15.2 Computational Methods: FEA & Numerical Analysis

PE Mechanical: Machine Design and Materials Exam

Track your progress as you read

Finished sections get a check mark.

1 Chapter 1: Basic Engineering Practice: Drawings, Economics & Design Methodology

1.1 Engineering Terms, Symbols, Drawing Interpretation & Conventions 1.2 Project Management & Engineering Economics 1.3 Design Methodology, Risk Assessment & Safety

2 Chapter 2: Physical & Mechanical Properties of Engineering Materials

2.1 Stress-Strain Relationships & Fundamental Mechanical Properties 2.2 Physical Properties, Temperature Effects, Hardness, Impact Toughness & Creep

3 Chapter 3: Mechanics of Materials: Fundamental Loadings & Deflections

3.1 Axial Loading, Normal Stress & Statically Indeterminate Members 3.2 Direct Shear, Bearing & Transverse Shear in Beams 3.3 Torsion of Shafts, Power Transmission & Thin-Walled Tubes 3.4 Beam Bending Stresses & Deflection Analysis

4 Chapter 4: Mechanics of Materials: Stability, Multiaxial Stress & Thermal Effects

4.1 Column Buckling & Elastic/Inelastic Instability 4.2 Combined Loading, Principal Stresses & Mohr's Circle 4.3 Thermal Stresses & Interference Stresses

5 Chapter 5: Failure Theories: Static Yielding & Fracture Mechanics

5.1 Static Failure Theories for Ductile Materials: von Mises & Tresca 5.2 Static Failure Theories for Brittle Materials & Fracture Mechanics

6 Chapter 6: Fatigue Analysis & Dynamic Failure

6.1 High-Cycle Fatigue, S-N Curves & Marin Modifying Factors 6.2 Fluctuating Stresses, Failure Criteria & Cumulative Damage

7 Chapter 7: Mechanical Attachments: Fasteners & Bolted Joints

7.1 Threaded Fasteners, Proof Strength & Preload Mechanics 7.2 Bolted Joint Stiffness, External Loading & Separation 7.3 Non-Threaded Fasteners: Keys, Pins, Clevises, Retaining Rings, Rivets, Lugs & Anchors

8 Chapter 8: Mechanical Attachments: Permanent Joints & Bonding

8.1 Welded Joints & Weld Group Stress Analysis 8.2 Brazing, Soldering & Structural Adhesive Bonding

9 Chapter 9: Power Transmission: Shafts & Gearing

9.1 Transmission Shaft Design, Critical Speeds & Deflection 9.2 Spur Gears: Geometry, Kinematics & AGMA Stress Equations 9.3 Helical, Bevel & Worm Gearing

10 Chapter 10: Power Transmission: Bearings, Flexible Drives, Friction Devices & Prime Movers

10.1 Rolling-Element Bearings & L10 Life Calculations 10.2 Hydrodynamic Journal Bearings & Lubrication Regimes 10.3 Belts, Roller Chains & Flexible Drives 10.4 Clutches, Brakes & Power Screws 10.5 Electric Motors, Engines & Flywheels

11 Chapter 11: Mechanical Components: Pressure Vessels & Structural Systems

11.1 Pressure Vessels & Piping: Thin and Thick Walls, ASME VIII and B31 11.2 Mechanical Springs: Helical, Belleville & Leaf Springs 11.3 Structural Trusses, Rigid Frames & Energy Methods

12 Chapter 12: Mechanical Dynamics, Vibrations & Fluid Power

12.1 Single-Degree-of-Freedom Vibrations & Resonance 12.2 Vibration Isolation, Multi-DOF Systems & Transmissibility 12.3 Fluid Power Systems: Hydraulics & Pneumatics

13 Chapter 13: Kinematics, Mechatronics & Control Systems

13.1 Mechanism Kinematics, Linkages, Cams & Simple Machines 13.2 Basic Mechatronics, Sensors & Closed-Loop Controls

14 Chapter 14: Supportive Knowledge: Fits, Tolerances & Manufacturing Methods

14.1 Limits, Fits & Geometric Dimensioning & Tolerancing (GD&T) 14.2 Manufacturing Methods, Machinability, Heat Treatment & Assembly

15 Chapter 15: Supportive Knowledge: Codes, Computation, Testing & Corrosion

15.1 Engineering Codes & Standards 15.2 Computational Methods: FEA & Numerical Analysis 15.3 Instrumentation, Testing, Inspection & Quality 15.4 Chemical Processes: Corrosion, Oxidation & Embrittlement

15.2 Computational Methods: FEA & Numerical Analysis

Key Takeaways Explain This Topic AI

In Finite Element Analysis (FEA), mesh convergence must be demonstrated for physical stress concentrations, whereas sharp re-entrant corners produce artificial mathematical stress singularities where calculated peak stress diverges toward infinity as mesh size approaches zero.

Boundary conditions cause most FEA errors: over-constraint adds false stiffness and stress peaks, while missing constraints leave a singular stiffness matrix; reactions must always balance the applied loads.

Newton-Raphson iteration x n + 1 = x n − f ( x n ) / f ′ ( x n ) x_{n+1} = x_n - f(x_n)/f'(x_n) x n + 1 ​ = x n ​ − f ( x n ​ ) / f ′ ( x n ​ ) converges quadratically near a root, while bisection is slower but guaranteed once the root is bracketed.

Simpson's 1/3 rule ∫ y d x ≈ ( h / 3 ) [ y 0 + 4 y 1 + 2 y 2 + ⋯ + y n ] \int y\,dx \approx (h/3)[y_0 + 4y_1 + 2y_2 + \dots + y_n] ∫ y d x ≈ ( h /3 ) [ y 0 ​ + 4 y 1 ​ + 2 y 2 ​ + ⋯ + y n ​ ] needs an even number of intervals and is more accurate than the trapezoidal rule for smooth data.

Last updated: October 2026

15.2 Computational Methods: FEA & Numerical Analysis

Modern machine design relies on computer-aided engineering (CAE), but the exam tests whether you understand what the software is doing: how a finite element model is built, when its answers can be trusted, and which numerical methods sit behind it. Hand-checkable numerical techniques also appear directly in problems.

Finite Element Analysis (FEA) Fundamentals

Finite Element Analysis discretizes continuous structures into finite subdomains (elements) to solve boundary-value stress, thermal, and dynamic problems.

Theoretical Formulation: The Stiffness Matrix Equation

In structural FEA, equilibrium is governed by the global stiffness equation: [ K ] { u } = { F } [K] \{u\} = \{F\} [ K ] { u } = { F }

where [ K ] [K] [ K ] is the global structural stiffness matrix ( N × N N \times N N × N , symmetric and banded), { u } \{u\} { u } is the vector of unknown nodal displacements, and { F } \{F\} { F } is the vector of applied nodal forces and boundary reactions.

The element stiffness matrix [ k ] e [k]^e [ k ] e is derived using the principle of virtual work: [ k ] e = ∫ V e [ B ] T [ D ] [ B ] d V [k]^e = \int_{V^e} [B]^T [D] [B] \, dV [ k ] e = ∫ V e ​ [ B ] T [ D ] [ B ] d V

where [ B ] [B] [ B ] is the strain-displacement matrix (spatial derivatives of shape functions N i N_i N i ​ ) and [ D ] [D] [ D ] is the material constitutive elasticity matrix (containing Young's modulus E E E and Poisson's ratio ν \nu ν ).

Element Types & Formulations

1D Elements :

Truss / Bar : 2 nodes, transmits pure axial force ( T x T_x T x ​ ).

Beam (Euler-Bernoulli / Timoshenko) : Transmits axial force, shear force, bending moments, and torsion. Timoshenko formulation incorporates shear deformation for deep beams.

2D Continuum Elements :

Plane Stress : Assumes zero out-of-plane normal and shear stresses ( σ z = τ x z = τ y z = 0 \sigma_z = \tau_{xz} = \tau_{yz} = 0 σ z ​ = τ x z ​ = τ yz ​ = 0 ). Applicable to thin plates loaded strictly in their plane.

Plane Strain : Assumes zero out-of-plane normal and shear strains ( ϵ z = γ x z = γ y z = 0 \epsilon_z = \gamma_{xz} = \gamma_{yz} = 0 ϵ z ​ = γ x z ​ = γ yz ​ = 0 ). Generates out-of-plane stress σ z = ν ( σ x + σ y ) \sigma_z = \nu(\sigma_x + \sigma_y) σ z ​ = ν ( σ x ​ + σ y ​ ) . Applicable to long, thick bodies of uniform cross-section (e.g., dams, tunnels, high-pressure pipes, retaining walls).

Axisymmetric : Models 3D bodies of revolution subjected to axisymmetric loading using a 2D planar section ( r , z r, z r , z ).

3D Solid Elements :

Tetrahedral (Tets) : Used for automated meshing of complex organic CAD geometry. First-order (4-node linear) tets must be avoided in stress analysis because they suffer from severe artificial shear locking (overly stiff). Second-order (10-node quadratic) tets possess mid-side nodes and provide excellent accuracy.

Hexahedral (Bricks) : 8-node linear or 20-node quadratic elements. Exhibit superior computational convergence, avoid shear locking, and are preferred for critical fatigue and contact stress zones.

• Mesh Convergence & The Stress Singularity Trap Peak Stress
• ^
• | /--- Stress Singularity (Sharp Corner r=0)
• | / (Stress -> Infinity as Mesh -> 0)
• | /
• | +---------------+------ Physical Stress Concentration (Fillet r>0)
• | / (Converges to True Peak Kt * sigma_nom)
• | /
• | /
• +---+-------------------------> Number of Elements (Refinement)

Mesh Convergence Study : The process of systematically refining element size ( h h h -refinement) or increasing element polynomial order ( p p p -refinement) until the output of interest (peak von Mises stress, displacement, or strain energy) changes by less than a specified threshold (typically < 2 % < 2\% < 2% to 5 % 5\% 5% ).

Physical Stress Concentration vs. Stress Singularity :

Physical Stress Concentration : Occurs at a geometric transition with a finite physical radius (e.g., a shaft fillet radius r = 3 mm r = 3\,\text{mm} r = 3 mm or an elliptical hole). As the mesh is refined, the calculated stress converges smoothly to the true physical stress: σ max ⁡ = K t σ nom \sigma_{\max} = K_t \sigma_{\text{nom}} σ m a x ​ = K t ​ σ nom ​ .

Stress Singularity : Occurs at idealized mathematical anomalies—such as a point load, a fixed boundary corner constraint, or a sharp re-entrant internal corner with zero radius ( r = 0 r = 0 r = 0 ) . Classical elasticity theory dictates that stress at a sharp crack or zero-radius corner is theoretically infinite ( σ ∝ r − λ → ∞ \sigma \propto r^{-\lambda} \rightarrow \infty σ ∝ r − λ → ∞ ).

Warning

In an FEA model with a stress singularity, mesh refinement will cause the calculated peak stress to increase without bound on every iteration . A design engineer must recognize a singularity, avoid using its diverging peak stress to calculate safety factors, and instead model the true manufacturing fillet radius or evaluate stresses several element rows away.

Verification vs. Validation (V&V in Computational Mechanics)

Verification : "Solving the equations right." The mathematical verification that the software correctly solves the underlying discrete model equations (e.g., checking mesh convergence, verifying reaction force equilibrium ∑ F = 0 \sum F = 0 ∑ F = 0 , confirming energy balance, comparing against closed-form beam solutions).

Validation : "Solving the right equations." The physical confirmation that the computational simulation accurately reflects real-world physics (e.g., comparing FEA strain and displacement predictions against physical strain gauge telemetry, photoelastic experiments, or full-scale destructive load tests).

Building a Trustworthy Model

Workflow : geometry (CAD) is simplified (defeatured), meshed, given material properties, loads and boundary conditions, solved, and post-processed.

Boundary conditions cause most errors. Over-constraining (fixing every degree of freedom on a face) adds false stiffness and creates artificial stress peaks; under-constraining leaves rigid-body motion and a singular stiffness matrix.

Symmetry : modeling one half or one quarter with symmetry constraints reduces size, but only if geometry, loads and supports are all symmetric.

Sanity checks : reaction forces must balance the applied loads, deformed shapes must make physical sense, units must be consistent (a model built in millimeters with stresses in MPa needs forces in newtons), and key results should be compared with a hand calculation.

• Numerical Methods Behind Engineering Software
• Root Finding

Newton-Raphson iteration solves f ( x ) = 0 f(x) = 0 f ( x ) = 0 by repeatedly following the tangent: x n + 1 = x n − f ( x n ) f ′ ( x n ) x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)} x n + 1 ​ = x n ​ − f ′ ( x n ​ ) f ( x n ​ ) ​

It converges quadratically near a simple root (the number of correct digits roughly doubles each step) but can diverge from a poor starting guess. Bisection is slower but always converges once the root is bracketed by a sign change.

Example: for f ( x ) = x 3 − 2 x − 5 f(x) = x^3 - 2x - 5 f ( x ) = x 3 − 2 x − 5 with x 0 = 2 x_0 = 2 x 0 ​ = 2 : x 1 = 2 − ( − 1 ) / 10 = 2.1 x_1 = 2 - (-1)/10 = 2.1 x 1 ​ = 2 − ( − 1 ) /10 = 2.1 ; x 2 = 2.1 − 0.061 / 11.23 = 2.09457 x_2 = 2.1 - 0.061/11.23 = 2.09457 x 2 ​ = 2.1 − 0.061/11.23 = 2.09457 ; x 3 = 2.094551 x_3 = 2.094551 x 3 ​ = 2.094551 , already correct to six digits. Nonlinear FEA (contact, plasticity, large deflection) uses the same idea, iterating each load step until the force residual is small.

Numerical Integration

For equally spaced data with step h h h : Trapezoidal rule: ∫ y d x ≈ h [ y 0 + y n 2 + y 1 + y 2 + ⋯ + y n − 1 ] \text{Trapezoidal rule: } \int y\,dx \approx h\left[\frac{y_0 + y_n}{2} + y_1 + y_2 + \dots + y_{n-1}\right] Trapezoidal rule: ∫ y d x ≈ h [ 2 y 0 ​ + y n ​ ​ + y 1 ​ + y 2 ​ + ⋯ + y n − 1 ​ ] Simpson’s 1/3 rule (even number of intervals): ∫ y d x ≈ h 3 [ y 0 + 4 y 1 + 2 y 2 + 4 y 3 + ⋯ + 4 y n − 1 + y n ] \text{Simpson's 1/3 rule (even number of intervals): } \int y\,dx \approx \frac{h}{3}\left[y_0 + 4y_1 + 2y_2 + 4y_3 + \dots + 4y_{n-1} + y_n\right] Simpson’s 1/3 rule (even number of intervals): ∫ y d x ≈ 3 h ​ [ y 0 ​ + 4 y 1 ​ + 2 y 2 ​ + 4 y 3 ​ + ⋯ + 4 y n − 1 ​ + y n ​ ]

Example: a spring test records force 0 0 0 , 2.4 2.4 2.4 , 4.2 4.2 4.2 , 5.4 5.4 5.4 and 6.0 kN 6.0\,\text{kN} 6.0 kN at 0 0 0 , 5 5 5 , 10 10 10 , 15 15 15 and 20 mm 20\,\text{mm} 20 mm . Simpson's rule gives the work W = ( 5 / 3 ) [ 0 + 4 ( 2.4 ) + 2 ( 4.2 ) + 4 ( 5.4 ) + 6.0 ] = ( 5 / 3 ) ( 45.6 ) = 76.0 J W = (5/3)[0 + 4(2.4) + 2(4.2) + 4(5.4) + 6.0] = (5/3)(45.6) = 76.0\,\text{J} W = ( 5/3 ) [ 0 + 4 ( 2.4 ) + 2 ( 4.2 ) + 4 ( 5.4 ) + 6.0 ] = ( 5/3 ) ( 45.6 ) = 76.0 J (since kN ⋅ mm = J \text{kN}\cdot\text{mm} = \text{J} kN ⋅ mm = J ); the trapezoidal rule gives 5 [ ( 0 + 6.0 ) / 2 + 2.4 + 4.2 + 5.4 ] = 75.0 J 5[(0 + 6.0)/2 + 2.4 + 4.2 + 5.4] = 75.0\,\text{J} 5 [( 0 + 6.0 ) /2 + 2.4 + 4.2 + 5.4 ] = 75.0 J .

Differential Equations and Linear Systems

Euler's method advances y ′ = f ( t , y ) y' = f(t, y) y ′ = f ( t , y ) by y n + 1 = y n + h f ( t n , y n ) y_{n+1} = y_n + h\,f(t_n, y_n) y n + 1 ​ = y n ​ + h f ( t n ​ , y n ​ ) ; it is simple but first-order accurate and can go unstable with large steps. The fourth-order Runge-Kutta method averages four slope estimates per step and is the usual general-purpose choice for dynamics and vibration simulations.

FEA ultimately solves the linear system [ K ] { u } = { F } [K]\{u\} = \{F\} [ K ] { u } = { F } by Gaussian elimination or a related factorization (or by iterative solvers for very large models). Explicit dynamic codes instead march in very small time steps limited by the time a stress wave takes to cross the smallest element.

Errors come from discretization (element size, time step), truncation of the numerical method and round-off ; refining the mesh or step and confirming that results stop changing is the universal check.