/structural-computation

title: Structural Computation
description: Finite element analysis fundamentals, form-finding methods, shell and gridshell structures, topology optimization, structural optimization, material-aware computation, and computational structural tools for AEC
version: 1.0.0
tags: [structural, FEA, form-finding, topology-optimization, shell, gridshell, Karamba, Kangaroo, funicular]
auto_activate: true
user_invocable: true
invocation: /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 Type | Dimensionality | DOF per Node | Captures | Typical Use | |---|---|---|---|---| | **Truss** | 1D (axial only) | 2 (2D) or 3 (3D) translations | Axial force only | Trusses, cables, bracing | | **Beam** | 1D (axial + bending) | 6 (3 translations, 3 rotations) | Axial, shear, bending, torsion | Frames, columns, beams | | **Shell** | 2D (membrane + bending) | 5-6 per node | In-plane + out-of-plane | Floors, walls, shells, slabs | | **Plate** | 2D (bending only) | 3 (1 translation, 2 rotations) | Out-of-plane bending | Floor slabs | | **Solid** | 3D | 3 translations | Full 3D stress state | Connections, nodes, foundations | | **Cable** | 1D (tension only) | 3 translations | Tension, large displacements | Cables, tendons | | **Spring** | 0D/1D | variable | Stiffness in specific DOFs | Supports, connections |

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:

• **As