Home / Science / Materials science / Phase and thermal behaviour / Thermal mismatch stress calculatorTool #852

Thermal mismatch stress calculator: thin-film stress, Stoney wafer bow, bimetal strip and constrained bar

Work out the stress from different thermal expansion: the biaxial stress in a thin film on a substrate, film stress from wafer bow or curvature with Stoney's equation, the curvature and layer stresses of a bonded bimetal strip (Timoshenko), and the stress in a bar held between walls. The CTE can change linearly with temperature.

Keep it one tap away for everyday use.

Bookmark Thermal mismatch stress calculator?

Press Ctrl + D now. Your browser opens its own bookmark box, where you choose the name and folder.

  1. Materials
  2. ΔT
  3. Stress
  4. Bow
Nothing is uploadedWorks offline once open

1. What to work out

T₀ is where there is no stress, such as the deposition or bonding temperature. CTEs are in ppm/K (10⁻⁶ per K).

2. Film on a substrate

σ = Ef/(1 − νf) × (εs − εf), with ε = ∫α dT from T₀ to T₁: the film is held to the substrate's size in both directions. The substrate values give the bow from Stoney's equation. The example values are round made-up numbers.

Download

What it does

When two bonded materials expand by different amounts as the temperature changes, neither can take its free size, and stress builds up. This page covers four everyday cases:

  • Film on a substrate. A thin coating is held to the substrate’s size in both directions, so its stress is biaxial: σ = E_f/(1 − ν_f) × (ε_s − ε_f), where ε = ∫α dT is each material’s free thermal strain from the stress-free temperature T₀ (for example the deposition temperature) to the final temperature T₁. It also gives the wafer curvature, radius and bow from Stoney’s equation, κ = 6σh_f(1 − ν_s)/(E_s h_s²), and checks Stoney against a full two-layer beam solution.
  • Stress from wafer bow. Measure the bow (or radius) of a wafer before and after a film goes on, and Stoney’s equation gives the film stress, σ = E_s h_s² κ / (6(1 − ν_s) h_f), and the force per width σh_f.
  • Bonded strip (bimetal). Two layers bonded along their length bend when heated. Timoshenko’s formula gives the curvature, and a force-and-moment balance gives the stress at the faces and the interface of each layer, the tip deflection of a cantilever, and the rise of a free strip.
  • Bar held between walls. A bar with fixed ends gets σ = −E ∫α dT. With a gap it is stress-free until the gap closes; the page gives the closing temperature and compares the stress with the yield strength.

Tick CTE changes with temperature to use α(T) = α + slope × (T − T_ref) for every material; the thermal strain is then integrated exactly.

How to use it

  1. Choose what to work out, and type the stress-free and final temperatures.
  2. Type the moduli (GPa), Poisson’s ratios, CTEs (ppm/K) and sizes. Positive stress is tension.
  3. Read the stress, curvature and checks, and the plot of stress against temperature or through the thickness. Download everything as an image, PDF or Word file.

Worked example

A made-up 1 µm film (E_f = 70 GPa, ν_f = 0.33, α_f = 23 ppm/K) is deposited at 300 °C on a 500 µm wafer (α_s = 3 ppm/K, E_s = 130 GPa, ν_s = 0.28) and cooled to 20 °C. The mismatch strain is (3 − 23) × 10⁻⁶ × (−280) = 5600 µε, and the biaxial modulus is 104.5 GPa, so the film carries 585.1 MPa in tension. Stoney gives a curvature of 77.77 /km, a radius of 12.86 m, and a bow of 97.21 µm over 100 mm. With the CTE slopes on (0.02 and 0.004 ppm/K²) the stress rises to 650.6 MPa.

Going the other way, a bow of 50 µm over 100 mm on the same wafer and film means a radius of 25 m and a film stress of 300.9 MPa.

A made-up bimetal strip of two 0.5 mm layers (bottom E = 140 GPa, α = 1.5 ppm/K; top E = 110 GPa, α = 19 ppm/K) heated from 20 to 120 °C bends to a radius of 382.3 mm (curvature 2.615 /m); a 50 mm cantilever tip moves 3.265 mm. A steel-like bar (E = 200 GPa, α = 12 ppm/K) with fixed ends heated by 50 K carries 120 MPa in compression.

Tips

  • Only the stress-free state matters: if a film was deposited with an intrinsic stress, add it to the thermal stress.
  • Stoney assumes a film much thinner than the substrate, a bow smaller than the thickness, and an isotropic substrate; the page says when you are outside that.
  • Near the free ends of a bonded strip there are peel and shear stresses at the interface that beam theory leaves out; that is where layers come apart.

FAQ

Why E/(1 − ν) for a film? The film is held in both in-plane directions at once, which makes it stiffer than in a simple pull by 1/(1 − ν).

Which sign does the bow have? Positive bow here means the film side is concave (dished), which comes from a film in tension.

Is my input stored? No, all the sums happen on this page.

Where the numbers come from

  • Thin-film biaxial stress, Stoney's equation (Stoney 1909) and Timoshenko's bimetal formula (Timoshenko 1925, J. Opt. Soc. Am. 11, 233) in the textbook forms of Freund and Suresh, Thin Film Materials, ch. 2. fibub's code is checked in the tests against SciPy (integrate.quad, brentq) and a NumPy multilayer beam solution.
  • The example moduli, expansion coefficients, thicknesses and temperatures are made-up round values, not data for any real material.
  • Nothing you type is sent or stored.

Standard scientific values, such as constants, atomic weights and melting points, are facts found in every textbook; we name where ours were checked. If you are an author or publisher and want a value, a credit or a table changed or removed, tell us and we will act quickly.

Last updated 2026-11-09.

Related tools

Stress-strain curve analyser #846 Paste or open force-extension or stress-strain data from a tensile test. Get Young's modulus from a straight-line fit, the 0.2 % offset yield strength, tensile strength, uniform and fracture elongation, toughness, resilience, true stress-strain to the UTS and the strain-hardening exponent, with a plot and CSV. Laminate calculator #850 Get ply properties from fibre and matrix (rule of mixtures, Reuss or Halpin-Tsai), then use classical lamination theory for any stacking such as [0/±45/90]s: Q and Q-bar, the A, B and D matrices, effective laminate moduli, ply stresses and strains under in-plane forces and moments, and first-ply failure by maximum stress and Tsai-Wu. Diffusion calculator #851 Work out diffusion in solids with Fick's second law: the erfc profile under a constant surface concentration, a thin-film Gaussian source, or a diffusion couple. Get D from D₀ and Q at any temperature, the case depth for a target concentration, the time to reach a depth, the concentration profile plot, and an optional finite-difference check. AFM roughness analyser #861 Open AFM or profilometer height data as a text matrix, an x y z list or a 16-bit greyscale PNG or TIFF. Level it with a plane or polynomial (fitted to the whole image or a reference region) and line-by-line flattening, then get areal Sa, Sq, Ssk, Sku, Sz, profile Ra, Rq, Rsk, Rku, Rz, a dragged line profile, the step height between two regions, a height histogram and a colour map. 1D Schrödinger solver #923 Choose a textbook potential, type a formula or draw your own, and get the lowest energy levels and wavefunctions of the one-dimensional Schrödinger equation by the matrix (finite-difference) method, with the exact levels alongside for the box, harmonic oscillator, finite well, Morse, linear and Pöschl–Teller potentials. Adduct and neutral loss calculator #984 Type a formula or a neutral monoisotopic mass and get the m/z of 25 common electrospray adducts in positive and negative mode, the ions left after common neutral losses (water, ammonia, CO2, sugars, phosphate and your own), and match a list of observed m/z values within a ppm tolerance. With no molecule, it turns observed m/z values into candidate neutral masses.