Euler–Bernoulli beam theory
Theory and Physics
The Most Basic Beam Theory
Professor, what is the "Euler-Bernoulli beam"? Euler also appeared in buckling theory, right?
It's the bending theory of beams established by Leonhard Euler and Daniel Bernoulli in the 18th century. It is the most fundamental theory in structural mechanics and also the starting point for beam elements in FEM.
Basic Assumptions
What assumptions does it make?
Three Basic Assumptions:
1. Plane Sections Remain Plane Assumption — Cross-sections that were plane before deformation remain plane after deformation.
2. Orthogonality Assumption — Cross-sections remain orthogonal to the beam axis (neutral axis) after deformation.
3. Small Deformation — Deformations are sufficiently small.
Assumption 2 seems important. If the cross-section doesn't tilt, that means...
Correct. It means shear deformation is ignored. The beam deforms only due to bending, and the tilt of the cross-section due to shear is zero. This is the most significant feature of the Euler-Bernoulli beam and also its greatest limitation.
Governing Equation
Please tell me the differential equation for bending.
A fourth-order ordinary differential equation concerning deflection $w(x)$:
Here, $EI$ is the bending stiffness, and $q(x)$ is the distributed load.
A fourth-order differential equation! Does solving it require four integrations?
Yes. The four integration constants are determined by four boundary conditions. Two conditions are needed at each end (displacement and rotation, or shear force and moment).
Physical meaning of each derivative:
| Derivative | Physical Quantity | Formula |
|---|---|---|
| $w$ | Deflection | |
| $w' = dw/dx$ | Rotation angle $\theta$ | |
| $w'' = d^2w/dx^2$ | Curvature $\kappa = M/(EI)$ | $M = EI w''$ |
| $w''' = d^3w/dx^3$ | Shear Force | $V = -EI w'''$ |
| $w'''' = d^4w/dx^4$ | Distributed Load | $q = EI w''''$ |
Differentiating deflection four times returns the load. That's elegant.
Memorizing this relationship is very useful for verifying FEM results. It forms a chain: deflection → rotation angle → curvature → moment → shear force → load.
Effect of Ignoring Shear Deformation
How much error occurs when shear deformation is ignored?
The error is smaller when the beam's span-to-depth ratio ($L/h$) is larger:
| $L/h$ | Contribution of Shear Deformation | Accuracy of Euler-Bernoulli |
|---|---|---|
| > 20 | < 1% | Sufficiently accurate |
| 10 to 20 | 1 to 5% | Practically acceptable |
| 5 to 10 | 5 to 20% | Caution required |
| < 5 | > 20% | Inaccurate. Should use Timoshenko beam |
So differences start to appear around $L/h < 10$.
Typical steel beams (H-shaped steel) have $L/h = 15 \sim 25$, so Euler-Bernoulli is sufficient. However, sandwich panels or short connecting beams ($L/h < 5$) require the Timoshenko beam theory.
Beam Elements in FEM
What does the Euler-Bernoulli beam element look like in FEM?
A 2-node beam element, with 3 degrees of freedom per node (in 2D): deflection $w$, rotation angle $\theta$, and axial displacement $u$.
An important feature is that the shape functions are Hermite polynomials (cubic polynomials). Unlike standard FEM elements (Lagrange polynomials), rotation angles are also nodal variables. This allows the accurate solution of the fourth-order bending differential equation with just 2 nodes.
Two nodes can solve a fourth-order equation! Is one element exact?
For a uniform cross-section and constant distributed load, a single element yields the exact solution. This is a major advantage of the Euler-Bernoulli beam element. Nodes need to be placed at points where concentrated loads act, but otherwise, a coarse mesh is sufficient.
Summary
Let me organize the Euler-Bernoulli beam theory.
Key points:
- Classical beam theory ignoring shear deformation — Cross-sections always remain orthogonal to the neutral axis.
- $EI w'''' = q$ — Fourth-order ordinary differential equation.
- Sufficiently accurate for $L/h > 10$ — Applicable to slender beams.
- Uses Hermite interpolation in FEM — Exact solution for uniform beams with 2 nodes.
- Use Timoshenko beam for $L/h < 5$ — Shear deformation cannot be ignored.
The deflection formula for a cantilever beam $\delta = PL^3/(3EI)$ solved in strength of materials classes is precisely the solution of this theory. FEM beam elements discretize this theory; the principle is the same.
The Birth of Euler Beam Theory
The Euler-Bernoulli beam theory originates from Leonhard Euler's 1744 work "De Curvis Elasticis." The assumption that "cross-sections remain plane and orthogonal to the beam axis after deformation" was initially controversial, but for slender beams (slenderness ratio L/h>10), it still provides accuracy within 1% error today.
Physical Meaning of Each Term
- Inertia Term (Mass Term): $\rho \ddot{u}$, i.e., "mass × acceleration". Have you ever experienced being thrown forward during sudden braking? That "feeling of being carried forward" is precisely the inertial force. Heavier objects are harder to set in motion and harder to stop once moving. Buildings shake during earthquakes because the ground moves suddenly while the building's mass "gets left behind." In static analysis, this term is set to zero, assuming "forces are applied slowly enough that acceleration is negligible." It absolutely cannot be omitted for impact loads or vibration problems.
- Stiffness Term (Elastic Restoring Force): $Ku$ or $\nabla \cdot \sigma$. When you stretch a spring, you feel a "force trying to return it," right? That's Hooke's Law $F=kx$, the essence of the stiffness term. Now a question—an iron rod and a rubber band, which stretches more under the same force? Obviously the rubber band. This "resistance to stretching" is the Young's modulus $E$, which determines stiffness. A common misconception: "High stiffness ≠ strong." Stiffness is "resistance to deformation," strength is "resistance to failure"—they are different concepts.
- External Force Term (Load Term): Body forces $f_b$ (e.g., gravity) and surface forces $f_s$ (e.g., pressure, contact forces). Think of it this way—the weight of a truck on a bridge is a "force acting on the entire volume" (body force), while the force of the tires pushing on the road surface is a "force acting only on the surface" (surface force). Wind pressure, water pressure, bolt tightening force... all are external forces. A common mistake here: getting the load direction wrong. Intending "tension" but applying "compression"—it sounds like a joke, but it actually happens when coordinate systems rotate in 3D space.
- Damping Term: Rayleigh damping $C\dot{u} = (\alpha M + \beta K)\dot{u}$. Try plucking a guitar string. Does the sound continue forever? No, it gradually fades. That's because the vibration energy is converted to heat by air resistance and internal friction in the string. Car shock absorbers work on the same principle—they intentionally absorb vibration energy to improve ride comfort. What if damping were zero? Buildings would continue swaying forever after an earthquake. Since that doesn't happen in reality, setting appropriate damping is crucial.
Assumptions and Applicability Limits
- Continuum Assumption: Treats material as a continuous medium, ignoring microscopic heterogeneity.
- Small Deformation Assumption (for linear analysis): Deformations are sufficiently small compared to initial dimensions, and the stress-strain relationship is linear.
- Isotropic Material (unless otherwise specified): Material properties are independent of direction (anisotropic materials require separate tensor definitions).
- Quasi-Static Assumption (for static analysis): Ignores inertial and damping forces, considering only the balance between external and internal forces.
- Non-Applicable Cases: Large deformation/large rotation problems require geometric nonlinearity. Nonlinear material behaviors like plasticity or creep require constitutive law extensions.
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Displacement $u$ | m (meter) | When inputting in mm, unify loads and elastic modulus to MPa/N system. |
| Stress $\sigma$ | Pa (Pascal) = N/m² | MPa = 10⁶ Pa. Be careful of unit inconsistency when comparing with yield stress. |
| Strain $\varepsilon$ | Dimensionless (m/m) | Note the distinction between engineering strain and logarithmic strain (for large deformations). |
| Elastic Modulus $E$ | Pa | Steel: ~210 GPa, Aluminum: ~70 GPa. Note temperature dependence. |
| Density $\rho$ | kg/m³ | In mm system: tonne/mm³ (= 10⁻⁹ tonne/mm³ for steel). |
| Force $F$ | N (Newton) | Unify to N in mm system, N in m system. |
Numerical Methods and Implementation
Beam Element Stiffness Matrix
Please tell me the stiffness matrix for the Euler-Bernoulli beam element.
For 2 nodes, with 2 bending degrees of freedom per node ($w_i, \theta_i$), a 4×4 stiffness matrix:
Where do these numbers $12, 6L, 4L^2$ come from?
They are obtained by differentiating the Hermite shape functions to create the $B$ matrix and integrating $\int_0^L B^T EI B \, dx$. Since Hermite shape functions are cubic polynomials, curvature (second derivative) becomes linear, allowing exact integration.
Adding axial force degrees of freedom ($u_i$) gives an independent 2×2 stiffness:
This is combined with the bending $[K_b]$ to form a 6×6 (2D) or 12×12 (3D) beam element stiffness matrix.
Equivalent Nodal Loads
How is a distributed load applied to a beam element?
Convert the distributed load $q$ into equivalent nodal loads. For a uniformly distributed load:
A force of $qL/2$ and moments of $\pm qL^2/12$ at each node... that's the fixed-end moment for a fixed-fixed beam, right!
Perfect understanding. Equivalent nodal loads are equal to the fixed-end reactions with the sign reversed. Knowing this correspondence allows intuitive verification of whether equivalent nodal loads are calculated correctly.
Element Names by Solver
| Element | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 2-node Beam (EB) | CBAR | B23 (2D), B33 (3D) | BEAM3 (2D), BEAM4 |
| 2-node Beam (Timoshenko) | CBEAM | B21 (2D), B31 (3D) | BEAM188/189 |
Are Nastran's CBAR and CBEAM different?
CBAR is the Euler-
Related Topics
なった
詳しく
報告