Thermal Expansion — From Bridge Joints to Solder Fatigue

1. Thermal Expansion: Small Effects, Big Consequences

Thermal expansion coefficients are tiny numbers — typically $10^{-6}$ per kelvin. Yet these small numbers, multiplied by large structures, large temperature swings, and high stiffness, produce forces and stresses large enough to buckle railway tracks, shatter ceramic coatings, and fatigue solder joints on circuit boards through millions of thermal cycles.

In FEM, thermal loading is often overlooked in initial design phases but becomes critical in detailed design — particularly for structures exposed to temperature cycles (engines, electronics, bridges) where differential expansion between materials causes cumulative fatigue damage.

2. Linear Thermal Expansion: ΔL = αLΔT

The linear coefficient of thermal expansion (CTE) $\alpha$ characterizes how much a material lengthens per unit length per degree of temperature change:

For isotropic materials, the same expansion occurs in all directions. For anisotropic materials (carbon fiber composites, single crystals), CTE is a tensor with different values in different directions. For a 3D solid, the thermal strain in FEM is added to the mechanical strain:

The mechanical strain $\boldsymbol{\varepsilon}_{\text{mechanical}} = \boldsymbol{\varepsilon}_{\text{total}} - \boldsymbol{\alpha}\Delta T$ is what drives stress — thermal strain alone (in an unconstrained material) produces no stress.

3. CTE of Engineering Materials

Note the wide range — from Invar at 1.2 to aluminum at 23 ppm/K. Joining materials with very different CTEs creates stress when they're bonded and subjected to temperature change.

4. Constrained Thermal Stress: σ = EαΔT

When thermal expansion is fully constrained (zero displacement allowed), the entire thermal strain is converted to mechanical stress:

For steel heated by 100°C while fully constrained: $\sigma = 210 \times 10^3 \times 12 \times 10^{-6} \times 100 = 252$ MPa — already near yield of mild steel (235–355 MPa). A 200°C heating cycle would exceed yield, causing permanent plastic deformation.

Thermal Shock and Ceramics

Ceramics have low CTE (2–6 ppm/K) and high $E$ but very low fracture toughness. When a ceramic component experiences a sudden temperature gradient (thermal shock), the surface expands or contracts faster than the interior. The resulting thermal stress can exceed the ceramic's tensile fracture strength, causing catastrophic cracking. This is why refractory bricks, ceramic cutting tools, and glass cookware are designed to minimize thermal shock — and why thermal FEM with transient analysis is essential in high-temperature ceramic component design.

5. Bridge Expansion Joints and Railway Buckling

Bridge expansion joints accommodate the seasonal and diurnal thermal movement of bridge decks. Without them, thermal stress would cause either buckling (compression in summer) or tensile fracture (winter shrinkage).

Railway Track Stability

Continuous welded rail (CWR) has no expansion joints. Instead, it's laid at a specific "stress-free temperature" (typically 27°C in temperate climates, 35°C in hot climates) with a pre-tension applied during welding. The rail is stressed in tension at low temperatures and compression at high temperatures. At the maximum rail temperature (~60°C summer), the compressive stress:

If this compressive stress, combined with any horizontal mis-alignment or soft ballast, exceeds the lateral buckling resistance, the track buckles sideways — "sun kink." Railway infrastructure maintenance agencies monitor track temperature and issue speed restrictions on hot days based on exactly this FEM-validated calculation.

6. Bimetal Actuators and Thermostat Strips

A bimetallic strip bonds two metals with different CTEs. When heated, the metal with higher CTE expands more, causing the strip to bend toward the lower-CTE side. The relationship between temperature change and tip deflection for a bimetal strip of length $L$ and total thickness $h$:

For a 100mm strip, $\alpha_2 - \alpha_1 = 10 \times 10^{-6}$/K, $h = 1$ mm, $\Delta T = 50$ K: tip deflection = $3 \times 10^{-5} \times 50 \times 0.01 / 0.002 = 7.5$ mm — easily enough to open or close an electrical contact. Bimetal thermostats, circuit breakers, and HVAC actuators all work on this principle.

7. CTE Mismatch in Electronics Packaging

Electronic packaging suffers the most severe CTE mismatch problems in engineering. A silicon chip (CTE = 2.6 ppm/K) soldered to a copper PCB pad (CTE = 17 ppm/K) experiences repeated differential thermal strain each power cycle:

This 0.115% strain amplitude per thermal cycle, accumulated over 10,000–100,000 power cycles, fatigue-cracks the solder joint connecting chip to board. This is the dominant failure mode of flip-chip BGA (Ball Grid Array) packages in automotive electronics, where temperature swings from -40°C to +125°C are required for 15+ years.

Coffin-Manson Solder Fatigue Life

The Coffin-Manson relation predicts solder joint fatigue life $N_f$ from plastic shear strain amplitude $\Delta\gamma_p$:

FEM thermo-mechanical analysis computes the plastic strain distribution in each solder joint, feeding directly into this life prediction. The FEM-Coffin-Manson methodology is now the industry standard for BGA package qualification, with Ansys, Abaqus, and COMSOL all having validated implementations.

8. FEM Thermo-Structural Coupling

Thermo-structural FEM requires two analyses in sequence (or simultaneous coupling):

1. Thermal analysis: Solve the heat equation to get temperature field $T(x,y,z)$ or $T(x,y,z,t)$
2. Structural analysis: Apply the temperature field as a load; nodes expand according to their local temperature and CTE

In Abaqus, the thermal strain is automatically included in the constitutive model when *EXPANSION is defined:

*MATERIAL, NAME=STEEL
*ELASTIC
210000., 0.3
*EXPANSION
12.E-6,   ! CTE in ppm/K (coefficient per degree C)

*STEP
*STATIC
*TEMPERATURE, FILE=THERMAL_RESULTS.odb
*END STEP

Fully Coupled vs. Sequential Approach

For problems where structural deformation changes heat transfer (e.g., contact opens/closes under thermal loading), fully coupled thermal-structural analysis is needed. Use Abaqus coupled temp-displacement elements (C3D8T, S4T). For most engineering problems where the temperature field doesn't depend on the deformation, the sequential approach is adequate and significantly faster.

9. Cross-Topics