Structural Engineering: Mechanics, Finite Element Modeling, and Dynamics

Structural engineering is the branch of civil and mechanical engineering focused on designing structural frameworks capable of supporting applied static and dynamic loads without collapse, unacceptable deflection, or fatigue degradation.

This article provides the mathematical and computational foundations of structural analysis, from classical continuum mechanics and beam theories to matrix stiffness formulations, Finite Element Modeling (FEM), and seismic dynamic response spectra.

1. Static Equilibrium, Determinacy, and Stress States

Static Equilibrium in 3D

A structural system or substructure in static equilibrium satisfies the vector equations of force and moment equilibrium:

\sum \mathbf{F} = \mathbf{0}, \quad \sum \mathbf{M}_O = \mathbf{0}

In Cartesian coordinates, this expands into six scalar equilibrium equations:

\sum F_x = 0, \quad \sum F_y = 0, \quad \sum F_z = 0

\sum M_x = 0, \quad \sum M_y = 0, \quad \sum M_z = 0

Structural Determinacy and Indeterminacy

For a 2D planar frame or truss with $b$ members, $r$ support reaction components, and $j$ joints:

Static Determinacy ( $r + b = 2j$ for trusses, $r + 3b = 3j$ for frames): Unknown internal forces and reactions can be solved uniquely using equilibrium equations alone.
Static Indeterminacy ( $r + 3b > 3j$ ): Redundant internal load paths require compatibility equations and constitutive force-deformation relationships (e.g., Virtual Work or Matrix Stiffness methods).

   Determinate Simple Beam:                 Indeterminate Fixed-Fixed Beam:
     (Pinned)        (Roller)                 (Fixed)            (Fixed)
      /\              (O)                     |===|              |===|
     =====           =====                     |   |==============|   |
      ^               ^                        |===|              |===|
     [Rx, Ry]        [Ry]                     [Rx, Ry, Mz]       [Rx, Ry, Mz]
    (3 Unknowns = 3 Eqns)                    (6 Unknowns > 3 Eqns; Indeterminate)

2D Plane Stress and Mohr's Circle

For a 2D stress element subjected to normal stresses $\sigma_x, \sigma_y$ and shear stress $\tau_{xy}$ , the state of stress on a plane inclined at angle $\theta$ is:

\sigma(\theta) = \frac{\sigma_x + \sigma_y}{2} + \frac{\sigma_x - \sigma_y}{2} \cos 2\theta + \tau_{xy} \sin 2\theta

\tau(\theta) = -\frac{\sigma_x - \sigma_y}{2} \sin 2\theta + \tau_{xy} \cos 2\theta

The Principal Stresses $\sigma_1, \sigma_2$ (where shear stress $\tau = 0$ ) are the eigenvalues of the stress tensor:

\sigma_{1, 2} = \frac{\sigma_x + \sigma_y}{2} \pm \sqrt{\left(\frac{\sigma_x - \sigma_y}{2}\right)^2 + \tau_{xy}^2}

Mohr's Circle Representation:
          Shear Stress (τ)
                ^
                |         .(σ_avg, τ_max)
                |       / | \
                |     /   |   \  Radius R
     -----------+----+----+----+----------> Normal Stress (σ)
               O|   σ₂  σ_avg  σ₁
                |     \   |   /
                |       \ | /
                |         .

Yield Criteria for Multiaxial Stresses

Von Mises Yield Criterion (Distortional Energy Theory): Yielding occurs when the equivalent von Mises stress $\sigma_v$ equals yield strength $\sigma_y$ : $\sigma_v = \sqrt{\frac{1}{2}\left[(\sigma_1 - \sigma_2)^2 + (\sigma_2 - \sigma_3)^2 + (\sigma_3 - \sigma_1)^2\right]} \ge \sigma_y$
Tresca Yield Criterion (Maximum Shear Stress Theory): $\tau_{\max} = \frac{\sigma_1 - \sigma_3}{2} \ge \frac{\sigma_y}{2}$

2. Flexure, Shear, and Elastic Column Stability

Euler-Bernoulli Beam Theory

Euler-Bernoulli beam theory assumes that cross-sections remain plane and perpendicular to the deformed neutral axis during bending (valid for slender beams where length $L \gg$ depth $h$ ).

The governing 4th-order differential equation for beam deflection $w(x)$ under distributed load $q(x)$ is:

E I \frac{d^4 w}{dx^4} = q(x)

where $E$ is Young's modulus and $I = \int y^2 dA$ is the second moment of area.

Beam Kinematics and Stresses:
- Bending Moment:      M(x) = -EI · (d²w/dx²)
- Shear Force:         V(x) = dM/dx = -EI · (d³w/dx³)
- Flexural Stress:     σ(y) = - (M · y) / I
- Transverse Shear:    τ(y) = (V · Q) / (I · b)
                       where Q = ∫_y^top y' dA' (First moment of area)

Column Buckling and Euler Critical Load

For an idealized slender elastic column subject to axial compressive load $P$ , bifurcation buckling occurs at the Euler Critical Load $P_{\text{cr}}$ :

P_{\text{cr}} = \frac{\pi^2 E I}{(K L)^2} = \frac{\pi^2 E A}{\lambda^2}

where:

$K$ is the effective length factor ( $K = 1.0$ for pinned-pinned, $K = 0.5$ for fixed-fixed, $K = 0.7$ for fixed-pinned, $K = 2.0$ for cantilever).
$\lambda = \frac{K L}{r}$ is the Slenderness Ratio, with radius of gyration $r = \sqrt{I/A}$ .

Buckling Modes for Boundary Conditions:
  Pinned-Pinned (K=1.0)     Fixed-Pinned (K=0.7)     Fixed-Fixed (K=0.5)
        /\                          |===|                   |===|
        ( | )                      (  |  )                 |  |  |
        ( | )                       ( | )                 (   |   )
        ( | )                        (|)                   |  |  |
        /\                          /\                      |===|

3. Finite Element Method (FEM) and Matrix Stiffness Formulation

The Finite Element Method (FEM) transforms continuous differential equations into discrete matrix algebraic systems $K \mathbf{u} = \mathbf{F}$ .

2D Beam Element Stiffness Matrix

For a planar Euler-Bernoulli frame element with 2 nodes and 3 degrees of freedom per node (axial displacement $u_i$ , transverse displacement $v_i$ , and rotation $\theta_i$ ):

Local Degree-of-Freedom Vector:
u_e = [ u₁, v₁, θ₁, u₂, v₂, θ₂ ]ᵀ

The element stiffness matrix $k^e$ in local coordinates is:

k^e = \begin{bmatrix} \frac{EA}{L} & 0 & 0 & -\frac{EA}{L} & 0 & 0 \\ 0 & \frac{12EI}{L^3} & \frac{6EI}{L^2} & 0 & -\frac{12EI}{L^3} & \frac{6EI}{L^2} \\ 0 & \frac{6EI}{L^2} & \frac{4EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{2EI}{L} \\ -\frac{EA}{L} & 0 & 0 & \frac{EA}{L} & 0 & 0 \\ 0 & -\frac{12EI}{L^3} & -\frac{6EI}{L^2} & 0 & \frac{12EI}{L^3} & -\frac{6EI}{L^2} \\ 0 & \frac{6EI}{L^2} & \frac{2EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{4EI}{L} \end{bmatrix}

Coordinate Transformation and Global Assembly

Transforming element stiffness from local to global coordinates using direction cosines $c = \cos \theta, s = \sin \theta$ :

K^e = T^T k^e T

Global assembly sums element contributions at shared node DOFs:

\mathbf{K}_{\text{global}} = \sum_{e=1}^{N_e} A_e^T K^e A_e

Enforcing kinematic boundary conditions (supports where $u = 0$ ) removes matrix singularity, allowing direct solver solution $\mathbf{u} = \mathbf{K}^{-1} \mathbf{F}$ .

4. Structural Dynamics and Earthquake Engineering

Structures subjected to transient dynamic excitations (wind gusts, rotating machinery, earthquakes) require time-dependent equations of motion.

Single-Degree-of-Freedom (SDOF) Dynamic System

For a mass $m$ constrained by stiffness $k$ and viscous damper $c$ under base ground acceleration $\ddot{u}_g(t)$ :

m \ddot{u}(t) + c \dot{u}(t) + k u(t) = -m \ddot{u}_g(t)

Dividing by $m$ :

\ddot{u}(t) + 2 \zeta \omega_n \dot{u}(t) + \omega_n^2 u(t) = -\ddot{u}_g(t)

where:

Natural circular frequency: $\omega_n = \sqrt{k/m}$ (natural period $T_n = \frac{2\pi}{\omega_n}$ ).
Damping ratio: $\zeta = \frac{c}{2 m \omega_n} = \frac{c}{2 \sqrt{k m}}$ (typically $\zeta \approx 0.05$ for reinforced concrete and steel).

SDOF Mechanical Oscillator Model:
           Base Motion ü_g(t)
          ====================
             |            |
         [Damper c]  [Spring k]
             |            |
          +------------------+
          |      Mass m      | ---> Relative Displacement u(t)
          +------------------+

Response Spectrum Analysis (RSA)

An earthquake Response Spectrum plots the maximum peak response (spectral displacement $S_d$ , pseudo-velocity $S_v$ , or pseudo-acceleration $S_a$ ) of an SDOF oscillator as a function of period $T_n$ for a specified damping ratio $\zeta$ :

S_a(T_n, \zeta) = \omega_n^2 S_d(T_n, \zeta)

Design Elastic Response Spectrum (ASCE 7 / Eurocode 8):
  Spectral Accel Sa (g)
     ^
 S_DS|      +----------------+ (Constant Acceleration Plate)
     |     /                  \
     |    /                    \
 S_D1|   /                      \___ (Constant Velocity 1/T)
     |  /                            \___ (Constant Displ 1/T²)
     +--+----+------------------+----+--------------------> Period T (s)
        0   T_0                T_s   T_L

In equivalent lateral force procedures, total seismic base shear $V_b$ is:

V_b = C_s W = \frac{S_{DS}}{R / I_e} W

where $W$ is total seismic building weight, $R$ is the Response Modification Coefficient (accounting for non-linear inelastic ductility and energy dissipation), and $I_e$ is the structural importance factor.

5. Structural Systems Comparison

+---------------------------+-------------------+--------------------+------------------------+
| Structural Lateral System | Height Efficiency | Ductility (R factor)| Typical Failure Mode   |
+---------------------------+-------------------+--------------------+------------------------+
| Moment Resisting Frames   | Low/Mid (≤ 20 fl) | High (R = 8)       | P-Delta drift instability|
| Concentric Braced Frames  | Mid-Rise (≤ 40 fl)| Moderate (R = 6)   | Brace compressive buckl|
| Eccentric Braced Frames   | Mid/High (≤ 60 fl)| High (R = 8)       | Shear link plastic yield|
| Reinforced Concrete Walls | High-Rise (≤ 80 fl)| High (R = 5 - 7)   | Flexural crushing/shear|
| Outrigger & Mega-Columns  | Supertall (> 80 fl)| High               | Differential shortening|
+---------------------------+-------------------+--------------------+------------------------+

References

Hibbeler, R. C. (2018). Structural Analysis (10th ed.). Pearson.
Chopra, A. K. (2020). Dynamics of Structures: Theory and Applications to Earthquake Engineering (5th ed.). Prentice Hall.
Bathe, K. J. (2014). Finite Element Procedures (2nd ed.). Prentice Hall.
American Society of Civil Engineers. (2022). Minimum Design Loads and Associated Criteria for Buildings and Other Structures (ASCE/SEI 7-22). ASCE.
Gere, J. M., & Goodno, B. J. (2018). Mechanics of Materials (9th ed.). Cengage Learning.