Newton's Laws of Motion — The Foundation of Every FEM Simulation

1. Why Newton Still Rules CAE

Three hundred and fifty years after Newton published his Principia Mathematica, his three laws of motion remain the bedrock of every structural FEM code, every crash simulation, and every dynamic analysis run in industry today. Whether you're running Abaqus, LS-DYNA, NASTRAN, or Ansys, the solver's heart is still Newton — just expressed in matrix form and applied simultaneously to thousands or millions of nodes.

This article starts with the simplest possible statement of each law, builds your physical intuition with real engineering examples, and then shows you exactly how Newton's laws transform into the FEM equation of motion. By the end, reading a solver's documentation about "mass matrix", "damping matrix", and "stiffness matrix" will make complete sense.

2. First Law: Inertia and Static Equilibrium

Newton's First Law: An object at rest stays at rest, and an object in motion stays in motion at constant velocity, unless acted on by a net external force.

In engineering terms, this is the law of inertia. The critical insight for CAE is the special case of rest: if something isn't accelerating, the net force on it must be zero.

This is the entire basis of linear static FEM analysis. When you run a "static structural" simulation — checking whether a bracket will yield under a given load, verifying a weld's safety factor — you're applying Newton's First Law at every node: the forces must balance.

Inertia in Crash Simulations

The flip side of the First Law is that inertia resists changes in motion. In a crash event, the car body decelerates rapidly but the occupant (not yet in contact with anything) wants to keep moving at the original speed. This is why seatbelts, airbags, and crumple zones all exist — they provide controlled forces to decelerate the occupant gradually rather than abruptly.

In FEM, inertia appears in the mass matrix [M]. A heavy component resists acceleration; a light one accelerates easily. Getting the mass distribution right is critical for any dynamic simulation.

3. Second Law: F = ma and Vector Forces

Newton's Second Law: The net force on an object equals its mass times its acceleration.

The bold notation is important — force and acceleration are vectors. In 3D space, this single equation actually represents three scalar equations:

For a finite element model with $N$ nodes, each with 3 degrees of freedom, you have $3N$ equations of this form — all solved simultaneously. That's why FEM for real structures involves millions of equations.

Force Decomposition in Practice

A key skill in CAE pre-processing is resolving forces into components. Consider a bolt loaded at an angle — you need to find the tensile and shear components separately because they have different allowable limits and failure modes.

FEM handles this automatically through the stiffness matrix and coordinate transformations, but understanding the underlying force decomposition helps you set up boundary conditions correctly and interpret results intelligently.

4. Third Law: Action-Reaction and Contact Mechanics

Newton's Third Law: For every action there is an equal and opposite reaction.

This law is deceptively simple but has profound implications for FEM, especially in contact and interface modeling. When two surfaces come into contact, the contact force on body A from body B must exactly equal and oppose the contact force on body B from body A.

Contact in FEM: Where Third Law Gets Tricky

Implementing the Third Law in FEM contact algorithms is one of the hardest problems in computational mechanics. The challenge: you don't know in advance which nodes will be in contact, and the contact force must be found as part of the solution. Common approaches include:

• Penalty method: Apply a stiff spring between penetrating surfaces (violates the Third Law slightly but converges)
• Lagrange multiplier method: Enforce exact non-penetration; satisfies Third Law exactly but adds unknowns
• Augmented Lagrangian: Combines both; used in Abaqus Standard by default

5. Worked Example: Crash Test Simulation

Let's apply all three laws to the canonical CAE problem: a frontal crash test. A 1500 kg car strikes a rigid barrier at 50 km/h (13.9 m/s). The front crumple zone crushes over 0.4 m in 80 ms.

Step 1: Initial Momentum and Required Impulse

The car must go from 13.9 m/s to 0, so the impulse (change in momentum) is 20,850 N·s.

Step 2: Average Deceleration Force

Using impulse = F·Δt, with Δt = 0.080 s:

That's 261 kilonewtons — about 18 times the car's weight. This is why crash events are violent and why explicit FEM (LS-DYNA, Abaqus Explicit) is needed rather than static analysis.

Step 3: Peak Deceleration (g-level)

Structural members in the load path must survive ~18g of acceleration loading. The FEM model captures how this load distributes through the A-pillar, firewall, and floor pan.

Step 4: Energy Check (First Law of Thermodynamics)

All of this kinetic energy must be absorbed by plastic deformation of the crumple zone. The FEM analysis tracks internal energy absorption to verify that safety cells (cabin, B-pillar) absorb minimal energy while crumple zones absorb maximum.

6. From F=ma to the FEM Equation of Motion

Now let's make the connection explicit. Consider a single mass-spring-damper system — a model of a car suspension, or any vibrating structure:

Where:

• $m\ddot{u}$ — inertia force (Newton's 2nd Law: mass × acceleration)
• $c\dot{u}$ — damping force (viscous resistance, proportional to velocity)
• $ku$ — spring (stiffness) force (Hooke's Law: stiffness × displacement)
• $F(t)$ — external applied force (time-varying)

This is Newton's Second Law with all forces included — inertia equals all other forces. Now extend this to a structure with $N$ nodes and multiple degrees of freedom. Every node gets its own equation. Assembled into matrix form:

This is the FEM equation of motion — the heart of every dynamic structural analysis. Every term traces directly back to Newton's Second Law.

7. Understanding [M]ü + [C]u̇ + [K]u = F

Let's dissect each matrix in the equation of motion and understand its physical meaning:

Mass Matrix [M]

The mass matrix represents the inertia of the structure. For a lumped mass model, [M] is diagonal — each entry is just the mass associated with that node. For consistent mass formulation, off-diagonal terms appear, representing the coupled inertia between nodes within an element.

In explicit FEM (LS-DYNA), a diagonal [M] is critical because it allows the equation of motion to be solved node-by-node without matrix inversion — enabling the extremely fast explicit time integration used in crash simulation.

Stiffness Matrix [K]

Assembled from element stiffness matrices, [K] represents the structure's resistance to deformation. Each element contributes a local $[k_e]$ that is transformed to global coordinates and assembled. The result is typically sparse and symmetric:

For linear static analysis (no inertia, no damping), the equation collapses to $[\mathbf{K}]\mathbf{u} = \mathbf{F}$ — a large sparse system of linear equations.

Damping Matrix [C]

Physical damping is complex to model accurately. In practice, most FEM codes use Rayleigh damping:

Where $\alpha$ and $\beta$ are calibrated from measured structural damping ratios at two frequencies. This formulation preserves the mathematical convenience of modal analysis.

Time Integration Methods

8. Practical CAE Tips from Newton's Laws

Tip 1: Always Check Reaction Forces (3rd Law)

After any FEM analysis, sum the reaction forces at all fixed boundaries. They must equal the sum of all applied loads — Newton's Third Law. If they don't, something is wrong with your model setup. This 30-second check catches most boundary condition errors.

Tip 2: Know When to Use Static vs. Dynamic Analysis (1st Law)

If the loading rate is slow enough that inertia effects are negligible, use linear static. The key check: is the loading frequency much less than the structure's lowest natural frequency? A rule of thumb: if $f_{\text{load}} < 0.1 \times f_{n1}$, static analysis is adequate. For anything faster — impact, blast, seismic — you need dynamic.

Tip 3: g-Loading Simplification (2nd Law)

For structures subjected to constant acceleration (e.g., aircraft components under maneuver loads, components on vibrating machinery), you can apply the equivalent body force $\mathbf{b} = \rho\mathbf{g}_{\text{effective}}$ as a distributed load rather than running a full dynamic simulation. This is the "inertial relief" approach widely used in aerospace FEM.

Tip 4: Mass Scaling in Explicit FEM

In explicit FEM, the stable time step is $\Delta t \leq \frac{l_{\min}}{c}$, where $c$ is the wave speed in the material. Small elements force tiny time steps. A common workaround is "mass scaling" — artificially increasing the density of small elements to reduce wave speed and allow larger time steps. This is safe only when kinetic energy stays below ~1% of internal energy — monitor it carefully.

9. Cross-Topics