Hooke's Law, Young's Modulus, and the FEM Stiffness Matrix

1. Hooke's Law: The Simplest Model of Elasticity

Hooke's Law — named after Robert Hooke who stated it in 1676 — is perhaps the most important equation in structural engineering: the force required to deform a spring is proportional to that deformation.

Where $F$ is the force (N), $k$ is the spring constant or stiffness (N/m), and $x$ is the displacement (m). The minus sign in the restoring force form ($F = -kx$) indicates that the spring force opposes the displacement.

Hooke's Law is the linear elastic model. Its validity is limited to small deformations in the elastic regime — before yield for metals, before significant geometric change, and before the material's microstructure changes. Remarkably, this simple law underpins virtually all linear structural FEM.

2. Spring Constant k: What It Means Physically

The spring constant $k$ has units of N/m (force per unit displacement). It represents how "stiff" the spring is — a high $k$ means you need a lot of force for a small displacement (stiff), while a low $k$ means the spring is compliant (soft).

For a cylindrical bar under axial load (the simplest structural element), the effective spring constant is:

Where $E$ is Young's modulus (Pa), $A$ is cross-sectional area (m²), and $L$ is length (m). This equation encodes everything: stiff material ($E$), large section ($A$), or short length all increase stiffness.

3. Springs in Series and Parallel

Understanding how springs combine is essential for modeling bolted joints, layered structures, and support systems.

Springs in Series (same force, displacements add)

The total stiffness is dominated by the weakest spring. A rubber gasket in series with a steel bolt: the total joint stiffness is close to the gasket stiffness, not the bolt.

Springs in Parallel (same displacement, forces add)

Parallel springs stiffen each other. Four bolts sharing a load are modeled as four springs in parallel: $k_{\text{joint}} = 4k_{\text{bolt}}$.

4. Young's Modulus: Hooke's Law for Solids

For a solid material under uniaxial stress, Hooke's Law in continuum form is:

Where $\sigma = F/A$ is stress (Pa), $\varepsilon = \Delta L / L$ is strain (dimensionless), and $E$ is Young's modulus (Pa). The bar spring constant follows directly:

5. Poisson's Ratio and Multi-Axial Hooke's Law

When you stretch a material axially, it contracts laterally. This is Poisson's effect:

For most metals, $\nu \approx 0.25$–$0.35$. For rubber (nearly incompressible), $\nu \approx 0.49$. For cork, $\nu \approx 0$ (designed to not change radially when compressed). The generalized Hooke's Law in 3D:

In matrix form, this defines the constitutive matrix $[\mathbf{D}]$ used in FEM element stiffness computation. For an isotropic solid, $[\mathbf{D}]$ depends only on $E$ and $\nu$.

Why Poisson's Ratio Matters in FEM

For near-incompressible materials ($\nu \rightarrow 0.5$, rubber, biological tissue), standard FEM elements "lock up" — they become artificially stiff because they can't correctly represent volume-preserving deformation. Specialized mixed formulations (hybrid elements in Abaqus: C3D8H, C3D8RH) or reduced integration are required. If you model rubber with standard C3D8 elements and $\nu = 0.499$, your stiffness will be off by an order of magnitude.

6. The FEM Stiffness Matrix [K]

For a single bar element (the simplest FEM element) with two nodes and axial DOF:

The $+1$ diagonal terms mean: applying unit displacement to one node while holding the other creates a restoring force $EA/L$ at each node (push-pull). The $-1$ off-diagonal terms represent the coupling between nodes — displacing node 1 creates a force at node 2.

For a full 3D solid element with 8 nodes × 3 DOF = 24 degrees of freedom, the element stiffness matrix is 24×24. Assembly of all elements into the global $[\mathbf{K}]$ (which might be 1,000,000×1,000,000 but very sparse) follows the same principle.

Stiffness Matrix Properties

• Symmetric: $[\mathbf{K}]^T = [\mathbf{K}]$ — from Maxwell's reciprocity theorem
• Positive semi-definite: $\mathbf{u}^T[\mathbf{K}]\mathbf{u} \geq 0$ — elastic energy is never negative
• Singular (before BCs): An unrestrained structure can translate and rotate freely — zero-energy rigid body modes
• Sparse: Each node only connects to neighboring nodes — most $K_{ij}$ entries are zero

7. Beyond Hooke's Law: Nonlinear Stiffness

Many real engineering scenarios involve nonlinear stiffness — the force-displacement relationship is not a straight line. Common sources:

For nonlinear problems, FEM replaces the single solve $[\mathbf{K}]\mathbf{u} = \mathbf{F}$ with an iterative Newton-Raphson loop: linearize about the current state, solve for a correction, update, repeat until convergence.

8. Spring-Back in Sheet Metal Forming

Spring-back is one of the most important practical applications of Hooke's Law in manufacturing. When a metal sheet is bent and the forming tool is removed, elastic recovery causes the bend angle to decrease (spring-back). The amount of spring-back depends directly on the elastic stiffness relative to yield strength:

Higher yield strength (like ultra-high-strength steel, UHSS, at 1500 MPa) relative to E (210 GPa) means more spring-back. UHSS parts can spring back 5–15°, requiring the die to be deliberately over-bent to compensate.

The Bauschinger Effect in Cyclic Loading

When a metal is plastically deformed in one direction and then loaded in the opposite direction, it yields at a lower stress than expected — the Bauschinger effect. This is important for spring-back predictions and fatigue analysis. Kinematic hardening models in FEM (Chaboche, Armstrong-Frederick) capture this behavior, while simpler isotropic hardening does not. For UHSS forming simulation, using the wrong hardening model can give spring-back predictions that are off by 2–3°.

9. Cross-Topics