Structural Computation

Structural Computation

1. Computational Structural Design Philosophy

The Unity of Form and Force

Structure is not an afterthought applied to a completed form. In the computational paradigm, structure is form. The geometry of a building, a bridge, a canopy -- each is a direct expression of the forces acting upon it. When we allow force to generate geometry rather than merely verifying geometry against force, we unlock an entirely different class of architectural expression: one that is simultaneously more efficient, more beautiful, and more materially honest.

This represents a fundamental epistemological shift: from analysis-after-design to analysis-as-design. The former treats structure as a constraint to be checked; the latter treats structure as a generative engine. In computational structural design, the distinction between "architect" and "engineer" dissolves. The designer works directly with force, curvature, material flow, and stress trajectories as primary design media.

Performance-Driven Geometry

Performance-driven geometry emerges when structural performance metrics (stress utilization, deflection, weight, embodied carbon) become the objective functions of a generative process. The designer defines boundary conditions, loads, material palettes, and performance targets. The computation produces geometry. This is not a loss of authorship -- it is an expansion of the design space far beyond what manual intuition can explore.

Key principles of performance-driven geometry:

• Axial action over bending -- Funicular forms carry load through compression or tension, avoiding costly bending. A catenary cable and an inverted catenary arch are the purest examples.
• Curvature as stiffness -- A flat plate is structurally weak; a curved shell resists load through membrane action. Double curvature provides bidirectional stiffness. The egg is the canonical example.
• Material where stress exists -- Topology optimization removes material where stress is absent and concentrates it along principal stress trajectories. The result resembles trabecular bone.
• Hierarchy and redundancy -- Efficient structures distribute load through hierarchical branching (trees, Gothic vaults, diagrid systems) with multiple load paths for robustness.

Historical Lineage

• Gaudi (1852-1926) -- Hanging chain models for the Colonia Guell chapel: physical form-finding producing pure compression vaults. The method is exact for catenary forms under self-weight.
• Isler (1926-2009) -- Hanging cloth models (fabric + plaster, inverted). Thin concrete shells at Deitingen, Norwich, Sicli: 40m+ spans at 70-90mm thickness. Physical analogue of dynamic relaxation.
• Frei Otto (1925-2015) -- Soap film experiments for minimal surface form-finding. Munich Olympic Stadium (1972), Mannheim Multihalle (1975). Founded the Institute for Lightweight Structures (IL) at Stuttgart.
• Philippe Block (b. 1980) -- Thrust Network Analysis (TNA): computational framework for compression-only funicular surfaces. Armadillo Vault (2016), NEST HiLo (2021).
• Contemporary -- Real-time structural feedback (Karamba3D, Kangaroo Physics), multi-objective optimization coupling structure with daylighting/thermal/fabrication, ML surrogates enabling population-based search over thousands of variants.

2. Finite Element Analysis Fundamentals

Overview

Finite Element Analysis (FEA) discretizes a continuous structure into a mesh of finite-sized elements, each governed by simple constitutive equations. By assembling element stiffness matrices into a global system, FEA solves for displacements, from which strains and stresses are derived. FEA is the universal method for structural verification in AEC.

Element Types

Element Formulations

• Linear (first-order) elements: Triangles (3-node), quadrilaterals (4-node), tetrahedra (4-node), hexahedra (8-node). Linear interpolation of displacements within each element. Require finer meshes but are computationally cheap per element.
• Quadratic (second-order) elements: Mid-side nodes added (6-node triangle, 8-node quad, 10-node tet, 20-node hex). Capture curved geometry and stress gradients more accurately. Preferred for stress analysis near holes, notches, and concentrated loads.
• Reduced integration: Uses fewer Gauss points than full integration (e.g., 1 point for a 4-node quad instead of 4). Faster but can exhibit hourglass modes (zero-energy deformation patterns). Remedied by hourglass control or selective reduced integration.

Mesh Requirements for FEA

Structural FEA meshes differ fundamentally from visualization/rendering meshes:

• Aspect ratio -- Elements should be as equilateral as possible. Aspect ratios above 5:1 degrade accuracy. Target < 3:1.
• Minimum elements per span -- At least 4-6 elements across any structural span; 8-12 for accurate stress recovery.
• Mesh grading -- Refine mesh near stress concentrations (holes, corners, load application points, support reactions) and use coarser mesh in low-gradient regions.
• Element quality metrics -- Jacobian ratio > 0.5, warpage < 15 degrees, skewness < 60 degrees for quads.
• Mesh transitions -- Avoid sudden jumps in element size (ratio of adjacent elements < 2:1).
• Compatibility -- Nodes must match at element boundaries. Hanging nodes require constraint equations or multi-point constraints (MPCs).

Boundary Conditions

Load Types

• Point load (concentrated) -- Force applied at a single node. Units: kN (SI), kip (Imperial).
• Line load (distributed) -- Force per unit length along a beam element. Units: kN/m, kip/ft.
• Area load (pressure) -- Force per unit area on shell/plate elements. Units: kN/m², psf.
• Self-weight -- Gravity applied to all elements based on material density. Acceleration: g = 9.81 m/s² (32.2 ft/s²).
• Wind load -- Pressure distribution from wind codes (EN 1991-1-4, ASCE 7-22 Ch. 26-31). Varies with height, exposure, building shape.
• Seismic load -- Equivalent lateral force (ELF) or response spectrum analysis per EN 1998 or ASCE 7-22 Ch. 11-23.
• Thermal load -- Temperature change causing expansion/contraction. Gradient through cross-section causes bending.
• Prestress -- Applied internal force (post-tensioning tendons, cable pretension).

Load Combinations

Structural design requires checking multiple load combinations per code. Key combinations:

• Eurocode ULS: 1.35G + 1.5Q (fundamental); with psi factors for multiple variable actions
• ASCE 7 LRFD: 1.2D + 1.6L (gravity dominant); 1.2D + 1.0E + L (seismic); 0.9D + 1.0W (uplift)

See references/fea-fundamentals.md for complete load combination tables with psi factors.

The Stiffness Matrix Concept

The fundamental FEA equation: [K]{u} = {F} where [K] is the global stiffness matrix (assembled from element stiffness matrices), {u} is the displacement vector (unknowns), and {F} is the force vector. The system is solved for {u}, from which strains and stresses are derived.

Key element stiffness matrices (truss: EA/L terms; beam: EI/L^3 terms; shell: Et and Et^3/(12(1-nu^2)) terms) are derived and tabulated in references/fea-fundamentals.md.

Result Interpretation

• Von Mises stress: sigma_vm = sqrt(sigma_x^2 - sigma_xsigma_y + sigma_y^2 + 3tau_xy^2). Combined equivalent stress for ductile material yield check.
• Principal stresses: sigma_1, sigma_2, sigma_3 (eigenvalues of stress tensor). Directions indicate natural force flow.
• Utilization ratio: U = demand/capacity. U < 1.0 = safe. Color: green (< 0.5), yellow (0.5-0.8), orange (0.8-1.0), red (> 1.0).
• Deflection: Check against limits (L/250 total, L/360 live load typical). See references/fea-fundamentals.md for complete deflection limits table.

Convergence Checking

• h-method -- Refine mesh (halve element sizes) and re-run. If results change by < 5%, mesh is adequate.
• p-method -- Increase element polynomial order (linear to quadratic) and re-run. More efficient for smooth solutions.
• Adaptive mesh refinement -- Automated refinement based on error estimators (Zienkiewicz-Zhu, residual-based).
• Target: Stress convergence within 2-5% between successive refinements.

Units (Quick Reference)

Full units table and cross-section property tables in references/fea-fundamentals.md.

3. Form-Finding Methods

Form-finding discovers structural geometry in static equilibrium under given loads, producing shapes that carry load through membrane action (tension or compression) rather than bending. These funicular forms are inherently material-efficient.

3.1 Hanging Chain / Catenary Method

A freely hanging chain adopts a catenary curve y(x) = a*cosh(x/a) where a = H/w. Inverting gives a pure compression arch. For uniform horizontal load, the shape is parabolic y = wx^2/(2H). Only valid for a single load case -- changing load changes the funicular shape.

3.2 Force Density Method (FDM)

Assigns a force-density ratio q = F/L to each element, linearizing the equilibrium equations (Schek, 1974): [C^T * diag(q) * C] * {x} = {p_x} - [C_f^T * diag(q_f)] * {x_f}. Linear system = direct solution, no iteration, guaranteed equilibrium. Force densities are prescribed (not forces), requiring experience to choose meaningful q-distributions. Best for cable nets, membranes, tensile canopies. See references/form-finding-methods.md for full derivation and worked example.

3.3 Dynamic Relaxation

A fictitious dynamic system oscillates under applied loads with artificial damping until static equilibrium is reached. Algorithm: initialize positions -> compute internal forces -> compute residuals -> update velocities with damping -> update positions -> check convergence (kinetic energy or residual force < tolerance) -> iterate.

Kinetic damping (Barnes, 1988) is the standard method: reset all velocities to zero when total kinetic energy peaks. Extremely robust, no parameter tuning. Time step must satisfy Courant condition: dt < sqrt(2*m_min/k_max). Best for cable nets, membranes, inflatables, gridshells. See references/form-finding-methods.md for pseudocode and damping strategies.

3.4 Thrust Network Analysis (TNA)

Block's method (2009) finds compression-only funicular surfaces using reciprocal force diagrams. A form diagram (horizontal projection of the network) has a dual force diagram where edge lengths represent horizontal force magnitudes. Vertical coordinates are solved from equilibrium. Scaling the force diagram changes the rise-to-span ratio without breaking the compression-only constraint.

Applications: unreinforced masonry vaults, stone shells, tile vaulting, 3D-printed compression structures. Software: RhinoVAULT, RV2, compas_tna. See references/form-finding-methods.md for reciprocal diagram theory.

3.5 Particle-Spring Systems (Kangaroo Physics)

Principle: Structural elements are represented as particles (nodes with mass) connected by springs (elements with stiffness). Loads, constraints, and geometric objectives are expressed as "goals" that apply forces to particles. The solver iteratively moves particles toward equilibrium.

Kangaroo 2 (Daniel Piker) uses a position-based dynamics approach:

• Each Goal calculates target positions for its participating nodes
• The solver averages target positions from all goals at each node
• Convergence occurs when goals agree (residual forces vanish)

Common goal types:

• Length -- Spring with rest length and stiffness
• Angle -- Angular spring between adjacent edges
• Planarize -- Forces a set of points toward coplanarity
• OnMesh/OnSurface -- Constrains points to lie on a reference geometry
• Load -- Applies constant force (gravity, wind)
• Anchor -- Fixes points with high stiffness
• Pressure -- Applies uniform pressure to a closed mesh (inflation)
• SoapFilm -- Minimizes area (minimal surface)
• Hinge -- Resists bending at mesh edges
• Floor/Plane collision -- Prevents penetration through a plane

Best for: Interactive form-finding, real-time feedback, multi-physics simulation (structure + fabrication constraints simultaneously), educational exploration.

3.6 Minimal Surfaces

Definition: A minimal surface has zero mean curvature (H = 0) at every point. Equivalently, it locally minimizes surface area for a given boundary. Soap films naturally form minimal surfaces.

Mathematical characterization: Mean curvature H = (kappa_1 + kappa_2) / 2 = 0, where kappa_1 and kappa_2 are principal curvatures. This means the surface curves equally in opposite directions at every point (anticlastic geometry).

Classical minimal surfaces:

• Catenoid -- Only minimal surface of revolution. Generated by rotating a catenary curve.
• Helicoid -- Ruled minimal surface. Generated by a line sweeping helically.
• Enneper surface -- Self-intersecting, algebraically simple.

Triply periodic minimal surfaces (TPMS):

• Schwarz-P (Primitive) -- Cubic symmetry, channel-like openings on faces.
• Schwarz-D (Diamond) -- Channels along body diagonals.
• Gyroid -- No planes of symmetry, no straight lines. Channel structure.
• I-WP (I-graph Wrapped Package) -- Complex interconnected channels.
• Applications in architecture: lightweight infill structures, acoustic panels, thermal exchangers, 3D-printed facade elements.

Structural implications: Minimal surfaces carry uniform tension under uniform pressure. They are structurally efficient for membrane structures (tents, canopies) but require edge cables or rigid boundaries to anchor the tension field.

Method Comparison Table

4. Shell Structures

Classification by Curvature

Membrane Theory vs. Bending Theory

Membrane theory assumes the shell is so thin that it carries load purely through in-plane forces (N_x, N_y, N_xy) with no bending moments. Valid when:

• Shell thickness t << radius of curvature R (t/R < 1/20)
• Loading is smooth and distributed
• Boundaries allow free in-plane movement (membrane-compatible supports)
• No concentrated loads or abrupt geometry changes

Bending theory (general shell theory) includes both membrane and bending actions. Required when:

• Near edges and supports (boundary layer effects, with decay length ~ sqrt(R*t))
• At geometric discontinuities (openings, ribs, thickness changes)
• Under concentrated loads
• When the shell is thick relative to its curvature

Buckling Analysis for Thin Shells

Thin shells are extremely sensitive to buckling. The classical critical buckling pressure for a sphere under uniform external pressure is:

p_cr = 2E / sqrt(3(1-nu^2)) * (t/R)^2

However, real shells buckle at 15-50% of the classical prediction due to geometric imperfections. Knockdown factors are essential:

• NASA SP-8007 knockdown factors for cylindrical and spherical shells
• Eurocode EN 1993-1-6 shell buckling assessment (LSA, MNA, GMNIA analyses)
• Nonlinear buckling analysis (GMNIA) -- Geometrically and Materially Nonlinear Analysis with Imperfections. The gold standard. Requires seeding geometric imperfections (first eigenmode shape scaled to fabrication tolerance, typically t/10 to t/2).

Shell Thickness Optimization

Shell thickness can vary spatially to match structural demand:

• Regions of high membrane force -- increase thickness
• Regions near supports -- increase thickness (bending zone)
• Central regions far from boundaries -- minimize thickness
• Parametric thickness variation: linear, stepped, or continuously graded
• Optimization methods: sensitivity-based gradient descent, evolutionary algorithms, or topology optimization of shell thickness field

Geometric Stiffness: Curvature = Strength

A fundamental principle of shell design: curvature provides stiffness. A flat plate of thickness t has bending stiffness D = Et^3/[12(1-nu^2)]. A curved shell of the same thickness resists load primarily through membrane action, with effective stiffness proportional to Et (not Et^3). This means a shell can be dramatically thinner than a flat plate spanning the same distance.

• Doubling curvature approximately halves the required thickness
• Gaussian curvature (K = kappa_1 * kappa_2) provides biaxial stiffness
• Zero Gaussian curvature (developable surfaces) allows deformation along the zero-curvature direction
• Corrugation adds effective curvature in one direction (folded plates)

Edge Conditions and Boundary Effects

• Free edge -- Membrane forces must vanish at the edge. Edge beam or thickening required to collect forces.
• Simply supported edge -- Vertical reaction provided but no moment restraint. Causes local bending concentration.
• Fixed edge -- Full moment connection. Highest bending demand at edge.
• Ring beam / edge beam -- Collects horizontal thrust from dome or vault. Sized for the horizontal component of membrane force.
• Boundary layer width -- The distance over which bending effects decay from an edge into the membrane zone is approximately L_b = sqrt(R * t) * pi / sqrt(12(1-nu^2)). For a concrete dome with R = 20m, t = 100mm: L_b is approximately 2.5m.

Ribbed Shells and Stiffening Strategies

When a pure shell is insufficient (due to buckling, openings, or concentrated loads), stiffening is added:

• Ribs -- Beam elements on the shell surface. Orthogonal grids or radial-circumferential patterns.
• Corrugation / folding -- Adds effective depth without adding material.
• Sandwich construction -- Two thin skins separated by a core (foam, honeycomb, lattice).
• Grid stiffening -- Continuous shell replaced by a lattice shell (gridshell) where elements act as ribs.

Historical Reference Shells

Material Considerations for Shells

• Reinforced concrete -- Classic shell material. Formwork cost is the main challenge; fabric-formed or CNC-milled foam formwork reduce cost.
• Timber -- CLT/LVL panels CNC-cut for segments. Good for compression shells. Joints are critical.
• Steel -- Welded plates for free-form shells. Expensive but precise.
• GFRP -- Lightweight, formable, corrosion-resistant. Limited by creep and fire performance.
• ETFE -- Pneumatic cushions. Transparent, lightweight, not structural in the traditional sense.
• Brick/stone -- Compression-only shells designed with TNA. Block Research Group projects.

5. Gridshell Structures

Definition and Classification

A gridshell is a structure with the shape and stiffness of a double-curved shell but made from a grid of linear elements (beams or laths) rather than a continuous surface. The grid carries load through a combination of membrane action (in-plane forces) and bending in individual elements.

Single-layer gridshells: One layer of beams. Structurally efficient when the form is funicular. Susceptible to buckling due to low bending stiffness in the normal direction. Requires in-plane bracing (diagonal cables or rods) to resist shear.

Double-layer (or multi-layer) gridshells: Two layers of beams separated by a core depth, connected by diagonal bracing or shear connectors. Much stiffer and more stable. Better for non-funicular or asymmetric loading.

Elastic Gridshells (Bent from Flat)

The defining innovation: a flat grid of flexible laths is assembled on the ground and then pushed/pulled into a doubly-curved shape. The laths bend elastically to accommodate the curvature. Once in shape, the grid is braced and the structure is stiff.

Process:

1. Fabricate flat grid (typically with sliding/rotating connections)
2. Lift or push grid into target shape using cranes, scaffolding, or actuators
3. Lock connections (tighten bolts, weld, or add bracing)
4. The stored bending energy in the laths provides initial stiffness

Key constraint: The lath material must accommodate the maximum bending strain without failure. For timber laths: epsilon_max = t / (2R_min) < epsilon_allow. This limits the minimum radius of curvature for a given lath thickness.

Material choices: Timber (oak, larch, bamboo), GFRP rods, steel tubes (for larger radii).

Rigid Gridshells (Assembled in Shape)

Beam elements are fabricated straight (or pre-bent to match local curvature) and assembled on scaffolding in the final curved shape. Connections are moment-rigid or pin-connected.

• No residual bending stress from erection
• Allow larger members and longer spans
• More complex fabrication and erection
• Examples: British Museum Great Court (Foster + Partners), Smithsonian Courtyard

Node Design and Connection Types

The node is the critical detail -- it must transfer forces between converging members while accommodating geometric complexity and remaining economically fabricable.

Connection types: Slotted plate nodes (simple, economical for orthogonal grids), cast steel nodes (3D-printed or investment-cast for complex geometry), bolted steel plates (adjustable, tolerance-forgiving), laminated timber finger joints (CNC-cut for all-timber gridshells), 3D-printed metal nodes (topology-optimized, SLM/DMLS, minimum material), and proprietary systems (MERO/KK ball joints for space frames).

Bracing Strategies

Single-layer gridshells require in-plane shear stiffness (bracing) to resist asymmetric loads:

• Diagonal cables -- Lightweight, nearly invisible. Pretensioned for load reversal. Most common.
• Diagonal rods -- Can resist both tension and compression. Heavier than cables.
• Rigid panels -- Glass or metal panels acting as shear diaphragms. Eliminates need for separate bracing.
• Third layer of members -- Triangulation. Converts the grid from quad-dominant (shear-flexible) to triangle-dominant (shear-stiff).

Form-Finding for Gridshells

• Compass method -- A kinematic method for elastic gridshells. Starting from a flat Chebyshev net (constant member spacing), the net is draped onto a target surface while maintaining constant member lengths. The resulting shape satisfies the fabrication constraint of using equal-length laths.
• Projection method -- A target surface is defined (from dynamic relaxation or other form-finding). A flat grid is projected onto this surface. Member lengths vary, requiring custom cutting.
• Inverse hanging model -- A cable net is hung under gravity, inverted, and the node positions define the gridshell geometry. Combined with dynamic relaxation for accuracy.
• Optimization-based -- Minimizing a combination of structural weight, bending stress in laths, and deviation from a target surface. Multi-objective optimization using genetic algorithms or gradient methods.

Panel Planarity

For glazed gridshells, quadrilateral panels must be planar or near-planar (tolerance: < L/500 to L/1000 of panel diagonal). Strategies:

• Conical meshes -- Normals at each node lie on a cone, guaranteeing planar quad faces
• Planarization -- Post-processing optimization of node positions (Kangaroo Planarize goal, Evolute Tools)
• Triangulation -- Always planar but increases member/node count and cost
• Curved panels -- Cold-bent or hot-bent glass. Expensive but increasingly feasible

Key Gridshell Projects

6. Topology Optimization

Overview

Topology optimization determines the optimal distribution of material within a design domain to maximize or minimize an objective function subject to constraints. Unlike sizing optimization (which changes member sizes) or shape optimization (which moves boundaries), topology optimization can create entirely new topologies -- holes, branches, and connections that were not present in the initial design.

SIMP Method (Solid Isotropic Material with Penalization)

The most widely used topology optimization method in AEC.

Design variable: Element pseudo-density rho_e in [0, 1] where 0 = void and 1 = solid.

Penalized stiffness: E_e = rho_e^p * E_0 where p is the penalization factor (typically p = 3). The penalization drives intermediate densities toward 0 or 1, producing a clear solid/void result.

Objective: Minimize compliance (= maximize global stiffness):

min C = {F}^T {u} = sum(rho_e^p * u_e^T * k_e * u_e)

Subject to:

[K(rho)]{u} = {F}     (equilibrium)
V = sum(rho_e * v_e) <= V_target   (volume constraint)
0 < rho_min <= rho_e <= 1          (bounds)

Solution: Sensitivity analysis (dC/d_rho_e) via adjoint method, followed by optimality criteria (OC) update or Method of Moving Asymptotes (MMA).

Level-Set Method

Represents the material boundary as the zero-level contour of a level-set function phi(x). phi > 0 = solid, phi < 0 = void. The boundary evolves by solving a Hamilton-Jacobi equation driven by shape sensitivities.

Advantages: Crisp boundaries (no intermediate densities). Natural for manufacturing. Disadvantages: Cannot easily nucleate new holes (requires seeding). More complex to implement.

ESO/BESO (Evolutionary Structural Optimization)

ESO (Xie and Steven, 1993): Iteratively removes elements with the lowest stress or strain energy density. Simple but can get stuck in local optima.

BESO (Bi-directional ESO): Allows both removal and addition of elements. More robust than ESO. Uses sensitivity filtering and stabilization.

Advantages: Simple to implement. Produces clear 0/1 solutions. Disadvantages: Heuristic (no formal convergence proof for ESO). Can produce mesh-dependent results without filtering.

Ground Structure Method

Starts with a highly connected truss network (the "ground structure") spanning the design domain. Optimization removes members with zero or near-zero force, leaving the optimal truss topology. Produces discrete, fabricable structures directly.

Best for: Long-span roof trusses, bridges, tower structures.

Design Domain, Loads, and Constraints

Setting up a topology optimization problem requires careful definition of:

• Design domain -- The 3D region within which material can be placed. Defined by architectural/spatial constraints.
• Non-design domains -- Regions that must remain solid (connection zones, load introduction points) or void (openings, passages).
• Loads -- Applied forces and their distributions. Multiple load cases handled by summing weighted compliances.
• Supports -- Boundary conditions (fixed, pinned, spring).
• Volume fraction -- Target material usage as a percentage of the design domain. Typical values: 20-50%. Lower fractions produce more organic, branching topologies.
• Symmetry -- Mirror symmetry planes, rotational symmetry, cyclic repetition.

Filtering

Without filtering, topology optimization produces checkerboard patterns (alternating solid/void elements) and mesh-dependent results (different meshes give different topologies).

• Sensitivity filter -- Smooths sensitivity values by averaging over a neighborhood of radius r_min. The most common approach.
• Density filter -- Smooths densities directly. Produces gray (intermediate density) boundaries.
• Heaviside projection -- After density filtering, applies a smooth step function to sharpen the 0/1 boundary. Controlled by a parameter beta that increases during optimization (continuation method).
• Minimum member size -- Enforced by choosing r_min >= minimum desired feature size.

Penalization Factor

The standard SIMP penalization p = 3 works well for stiffness-based optimization:

• p = 1: No penalization, result is full of gray.
• p = 2: Moderate penalization. Some gray remains.
• p = 3: Standard. Most intermediate densities penalized away.
• p = 4-5: Aggressive penalization. Can cause convergence issues.

Continuation approach: Start with p = 1 and gradually increase to p = 3 over iterations. Improves convergence and can find better optima.

Interpretation of Results

Raw topology optimization output requires post-processing:

• Thresholding -- Set a density cutoff (e.g., rho > 0.5 = solid). Produces a binary solid/void result.
• Smoothing -- Apply Laplacian or Gaussian smoothing to the boundary. Removes staircase artifacts.
• CAD reconstruction -- Fit NURBS surfaces or B-rep geometry to the smoothed result. Manual or semi-automated (e.g., nTopology, Altair Inspire).
• Structural verification -- Re-analyze the interpreted geometry with FEA to confirm performance. The post-processed shape may differ from the optimized result.

2D vs. 3D Topology Optimization

2D optimization is used for:

• Planar components (brackets, gusset plates, flat connection details)
• Cross-sections of beams or walls
• Pedagogical examples and initial exploration
• Computation time: seconds to minutes

3D optimization is used for:

• Spatial structural nodes
• Transfer structures (load paths through a 3D volume)
• Foundation layouts
• Full building volumes
• Computation time: minutes to hours (can be very large for fine meshes)

Tools for Topology Optimization

AEC Applications

• Structural nodes -- Topology-optimized steel nodes for space frames and gridshells. 3D-printed in metal. Up to 75% weight reduction vs. conventional welded nodes. Arup's 3D-printed steel nodes for a structural tree.
• Floor plates -- Ribbed slabs with topology-optimized rib patterns. ETH Zurich NEST HiLo project: 3D-printed concrete floor with 70% less material than flat slab.
• Facade brackets -- Connection elements between facade panels and primary structure. Optimize for multiple load cases (wind, self-weight, seismic).
• Foundations -- Material distribution in transfer beams and pile caps. Strut-and-tie models are a classical topology optimization concept.
• Furniture and pavilions -- Increasingly common for one-off or small-series production. Branch Technology, AI Build, MX3D.

7. Material-Aware Computation

Material Properties Table

Anisotropic Materials

Timber: Strongly anisotropic. Properties differ along grain (longitudinal), across grain (radial), and tangential directions. E_L : E_R : E_T is approximately 20 : 1.6 : 1. Computational models must account for grain direction. CLT (Cross-Laminated Timber) alternates grain direction for quasi-isotropic behavior in-plane.

Fiber composites (GFRP, CFRP): Properties depend on fiber orientation. Unidirectional laminates are strongly anisotropic. Quasi-isotropic layups ([0/+45/-45/90]s) provide balanced in-plane properties but are weaker than aligned laminates in any single direction. Classical Laminate Theory (CLT -- confusingly same acronym) governs composite analysis.

Masonry: Anisotropic due to mortar joints. Different stiffness and strength along bed joints vs. head joints vs. diagonal. Homogenized masonry models treat the assembly as an equivalent anisotropic continuum.

Material Behavior Models

• Linear elastic -- Stress proportional to strain (Hooke's law). Valid for most materials at service load levels. The basis for most FEA in AEC.
• Elasto-plastic -- Linear elastic up to yield, then plastic deformation at constant (or hardening) stress. Required for steel design at ULS. Bilinear or multilinear stress-strain models.
• Viscoelastic -- Time-dependent deformation under sustained load (creep). Critical for concrete and timber. Modeled as spring-dashpot combinations (Maxwell, Kelvin-Voigt, Burgers models).
• Nonlinear elastic -- Stress-strain curve is nonlinear but unloading follows the loading path. Rubber, some polymers.
• Brittle -- No plastic deformation before failure. Concrete in tension, glass, unreinforced masonry. Requires fracture mechanics or damage models.

Composite Action

When two materials act together (steel-concrete beams, timber-concrete floors), transformed section analysis converts to an equivalent single-material section using modular ratio n = E_1/E_2 (full composite). Partial composite action (interface slip) is modeled with interface springs. Shear connectors (headed studs, screws) transfer horizontal shear.

Material Efficiency Metrics

Timber and bamboo are exceptional: high strength-to-weight, low-to-negative embodied carbon. CFRP is phenomenal structurally but carries enormous environmental cost. Steel and concrete are the workhorses with moderate efficiency.

Digital Material Systems

Emerging paradigm of spatially varying material properties:

• Functionally Graded Materials (FGM) -- Continuous variation of composition (e.g., concrete with graded fiber density)
• Variable-density lattices -- 3D-printed lattice infill matched to stress fields from topology optimization (nTopology, Altair)
• Multi-material printing -- Simultaneous deposition of stiff and flexible materials for spatially tuned stiffness
• Programmable materials -- Shape-memory alloys, 4D printing, materials responding to stimuli (temperature, moisture)

8. Structural Computation Tools

8.1 Karamba3D (Grasshopper)

Capabilities: Real-time FEA within Grasshopper. Beam and shell elements. Linear and second-order analysis. Cross-section optimization. Utilization checking. Eigen-frequency and buckling analysis. Large deformation analysis (geometrically nonlinear).

Limitations: No material nonlinearity (no concrete cracking, no steel yielding). No dynamic time-history analysis. Not a code-checking tool (no automatic Eurocode/ASCE capacity checks -- requires manual setup). Accuracy depends on mesh quality.

Integration: Fully embedded in Grasshopper. Takes Rhino geometry (lines, meshes) as inpu

…(truncated)