17 Equations That Changed the World-There’s More Than That!!!

Part 1

Author: Greg Caswell


 Download PDF



Mathematics has been a part of our lives forever and is used in numerous ways in our everyday lives. Recently, created by Ian Stewart, listed on Dr. Paul Coxon’s Twitter account and discussed on mathematics blogger Larry Philip’s site is a list of the “17 Equations that Changed the Course of History,” many of which have been mentioned on the Big Bang Theory TV series.

17 Equations Figure 1

 However, the list is incomplete. There are numerous equations that help to determine the reliability of electronics products that DfR believes should be included. These formulas are all integral to DfR’s Sherlock software and are the basis for the validation of the results obtained in an ADA analysis. Let’s take a look at them. 

Plated Through Hole Fatigue 

17 Equations Figure 1 a

Plated through holes (PTHs), also known as plated through vias (PTVs), are holes drilled through multilayer printed circuit boards (PCBs) that are electrochemically plated with a conductive metal (typically copper). These plated holes provide electrical connections between layers.

Because these plated holes are metallurgically bonded to annular rings on the top and bottom of the printed circuit board, they act like rivets and constrain the PCB. This constraint subjects the PTH to stresses when the PCB experiences changes in temperature.

PTH Fatigue is the circumferential cracking of plated through holes (PTHs) due to the differential expansion between the copper plating (~17 ppm) and the out-of-plane coefficient of thermal expansion (CTE) of the printed board (~45 to 70 ppm) during temperature variations. This failure mechanism was first reported by Bell Laboratories in 1976.

What are the drivers for PTH Fatigue?

 PTH Fatigue is influenced by maximum temperature, minimum temperature, PTH diameter, PTH copper plating thickness, copper plating material properties (ductility, yield strength), printed board thickness, printed board out-of-plane material properties (CTE, elastic modulus), and defects within the copper plating (voids, folds, etch pits, etc.).

How does the software assess PTH Fatigue? 

The software calculates a time to failure using the industry-accepted model published in IPC-TR-579, Round Robin Reliability Evaluation of Small Diameter Plated-Through Holes in Printed Wiring Boards. Life calculation for PTHs subjected to thermal cycling is a three step process. The first step is to calculate the stress being experienced by the copper barrel of the PTH. This is provided by the equations below,

 17 Equations Figure 2

 where a is coefficient of thermal expansion (CTE), T is temperature, E is elastic modules, h is the board thickness, d is the hole diameter, t is the plating thickness, and E and Cu correspond to board and copper properties, respectively.

 Once the stress is determined, the strain range is calculated by:

17 Equations Figure 3

 Where, Sy is the yield strength of copper. This strain is adjusted by two constants: a strain distribution (Kd) factor and a quality factor (KQ).

 17 Equations Figure 4

While the strain distribution factor tends to be set to a value of 1.6, KQ is dependent upon the quality of the PTH (i.e., the presence of defects such as voids, cracks, folds, etc.). The quality factor can range from 0 to 10 with the following delineations:

  • Extraordinary (KQ = 10)
  • Superior (KQ = 8.7)
  • Good (KQ = 6.7)
  • Marginal (KQ = 4.8)
  • Poor (KQ = 3.5)

Once the strain range is defined, the cycles to failure (Nf) can be calculated iteratively

 17 Equations Figure 5

 with Su being the ultimate tensile strength and Df being ductility of the plated copper.

 Reviewing the equations above, it can be seen that the designer and PCB manufacturer have the following controls over the reliability of PTHs

  • Out-of-plane CTE of the printed board
  • Plating thickness
  • Aspect ratio (hole diameter over board thickness)
  • Plating material properties (strength and ductility) and plating quality


Thermo-Mechanical Solder Joint Fatigue?

17 Equations Figure 6

 Solder joints, also known as interconnects, provide electrical, thermal, and mechanical connections between electronic components (passive, discrete, and integrated) and the substrate or board to which it is attached. Solder joints can be a first level (die to substrate) or second level (component package to printed board) connection. This module assesses the thermo-mechanical fatigue behavior of second-level solder joints.

 During changes in temperature, the component and printed board will expand or contract by dissimilar amounts due to differences in the coefficient of thermal expansion (CTE). This difference in expansion or contraction will place the second-level solder joint under a shear load. This load, or stress, is typically far below the strength of the solder joint. However, repeated exposure to temperature changes, such as power on/off or diurnal cycles, can introduce damage into the bulk solder. With each additional temperature cycle this damage accumulates, leading to crack propagation and eventual failure of the solder joint.

The failure of solder joints due to thermo-mechanical fatigue is one of the primary wearout mechanisms in electronic products, primarily because inappropriate design, material selection, and use environments can result in relatively short times to failure.

What are the drivers for Thermo-Mechanical Solder Joint Fatigue?

 Thermo-mechanical solder joint fatigue is influenced by maximum temperature, minimum temperature, dwell time at maximum temperature, component design (size, number of I/O, etc.), component material properties (CTE, elastic modulus, etc.), solder joint geometry (size and shape), solder joint material (SnPb, SAC305, etc.), printed board thickness, and printed board in-plane material properties (CTE, elastic modulus).

How does the software assess Thermo-Mechanical Solder Fatigue?

 The software calculates a time to failure using strain energy. The detailed methodology is provided by the equations below. The first equation, shown below, describes the force exerted on a solder joint during a thermal cycle.

 17 Equations Figure 7

where a is the CTE, T is temperature, L is one half component length, F is force, E is elastic modulus, A is effective solder joint area, G is the shear modulus, h is thickness and n is the Poisson ratio. The strain range induced in the solder joint during the thermal cycle is

 17 Equations Figure 8

The equations above are package-specific and account for the geometry, interconnect structure and material properties of the component and printed circuit board.

The stress on the solder joint is determined using the computed forces and this stress is combined with the strain to determine the energy dissipated by the solder during a thermal cycle through the following equation,

17 Equations Figure 9

The resulting strain energy is used to compute the number of cycles to failure for the component under temperature cycling using equations developed by Syed (see below).

 17 Equations Figure 10

What is the Shock and Vibration Module?

The Shock/Vibration module in Sherlock utilizes the finite element method to predict the circuit card assembly (CCA) response during mechanical shock and harmonic or random vibration events. The board response results are then used to make predictions on the robustness of the CCA to these events. In the case of vibration, high cycle fatigue predictions are made to determine the life of the interconnect (lead and solder joint). The shock analysis results are used to determine if a critical stress due to board bending is exceeded.

 Second level interconnects, provide electrical, thermal, and mechanical connections between electronic components (passive, discrete, and integrated) and the substrate or board to which it is attached. The vibration portion of the module assesses the high cycle fatigue behavior of secondlevel interconnects.

 How does the Software Assess Vibration?

 During exposure to vibration the circuit card assembly responds by cyclic deflecting in a manner that corresponds to its natural frequencies. These natural frequencies are dependent on the geometry, stiffness, mass and boundary conditions of the circuit card assembly. The software automatically generates the model and then utilizes the open source finite element analysis engine Calculix to conduct the shock and vibration response. Error! Reference source not found. shows vibration analysis of circuit board using Calculix. The Calculix FEA analysis engine was originally developed by Guido Dhondt1 and is open source covered under GNU General Public License2. The naming conventions and input style formats for CalculiX are based on those used by ABAQUS, a proprietary, general purpose finite element code developed and supported by Hibbitt, Karlsson & Sorensen, Inc (HKS) (http://www.abaqus.com). A discussion on CalculiX can be found at STRUCTURE (www.structuremag.org/article.aspx?articleID=1146), for validation of the CalculiX solver please visit bConverged (http://www.bconverged.com/benchmarks)

 17 Equations Figure 11

 The software generates a three layer shell model of the circuit card assembly with six-node triangular elements that define the top components, bottom components and printed wiring board. These six-node shell elements are expanded into 15 node three-dimensional wedge elements automatically by CalculiX during analysis.

The output of the finite element analysis is post-processed by the software and the maximum board level strains for every component are recorded for use in determining whether that component fails during exposure to the shock and vibration loads. There are four types of vibration analyses that the software conducts:

  1. Natural frequency extraction
  2. Single point harmonic vibration
  3. Swept harmonic vibration
  4. Random vibration

The first analysis, natural frequency extraction, is used to extract the fundamental frequencies of the circuit card assembly.

How does the Software Assess Vibration Fatigue using Board Strains?

The failure of solder joints due to vibration is based upon a technique similar to that developed by Steinberg. The main modification is converting from a displacement based criteria to a board level strain criteria. The Steinberg critical PCB deflection at which a component will survive 10 million cycles during harmonic vibration or 20 million cycles during random vibration is:

17 Equations 1 edit

The main issue with this approach is that it is limited to simple board geometries (since the maximum deflection is always assumed to be at the center of the PCB) and doesn’t account for board curvature. By utilizing finite element modeling the software eliminates some of the variables in the equation since they are accounted for in the FEA model. These include board length (B), position factor (r), board thickness (h), and the equation is simplified to:

17 Equations 1 edit3

The software uses this critical strain value and the FEA computed strain values to make fatigue predictions using the Basquin equation:

 17 Equations Figure 14

17 Equations 1 edit2

 Currently, the software does not make fatigue predictions based on a mechanical shock loading. Mechanical shock is instead viewed as an overstress event and is based on exceeding a predefined board level strain. This analysis is an extension to the type of information provided in IPC-9704 for BGA devices. The software utilizes these limits as the generic values for all components.

 17 Equations Figure 15

 If the strain level anywhere at the parts locations exceeds the acceptable amount the part is considered to fail the shock robustness portion of the analysis. Figure 3 shows plot of Maximum allowable principal strain vs. strain rate for circuit boards with variable thickness. The shock analysis utilizes the CalculiX based FEA to determine the circuit card assembly’s response and the part level strains which occur during the shock event.

What is the Coefficient of Thermal Expansion?

17 Equations Figure 16

Almost all materials exhibit a change in physical dimensions when subjected to temperature changes. The degree of expansion in response to a change in temperature is called the coefficient of thermal expansion.

The coefficient of thermal expansion (CTE) is critical because when two materials with different CTE are joined a stress is imparted because of the resultant displacement mismatch. This is the main driver of thermal-mechanical fatigue of electronic components.

During changes in temperature, the component and printed board will expand or contract by dissimilar amounts due to differences in the coefficient of thermal expansion (CTE). This difference in expansion or contraction will place the second-level solder joint under a shear load. This load, or stress, is typically far below the strength of the solder joint. However, repeated exposure to temperature changes, such as power on/off or diurnal cycles, can introduce damage into the bulk solder.

The calculation of the CTE of the printed circuit board is critical as input for determining thermomechanical fatigue of solder interconnects.

How is the PCB CTE Calculated

The Sherlock software uses lamina theory for calculating the CTE of printed circuit boards.

 Printed circuit boards (PCBs) are composed of alternating layers of glass-reinforced epoxy laminate / prepreg and copper foil

The mechanical properties of the laminate and copper foil are

  • Coefficient of thermal expansion, in-plane: CTExy
  • Coefficient of thermal expansion, out-of-plane: CTEz
  • Elastic modulus, in-plane: Exy
  • Elastic modulus, out-of-plane: Ez

The stackup tool allows the copper thicknesses to be in oz / mil / microns, where 1 oz = 35 microns = 1.4 mil

Term definitions

  • CTExypcb is the coefficient of thermal expansion, in-plane of the PCB
  • CTEzpcb is the coefficient of thermal expansion, out of plane of the PCB
  • Exypcb is the elastic modulus, in plane of the PCB
  • Ezpcb is the elastic modulus, out of plane of the PCB
  • CTExyln is the coefficient of thermal expansion, in-plane of layer n (the first and last layer will always be copper)
  • CTEzln is the coefficient of thermal expansion, out of plane of layer n
  • Exyln is the elastic modulus, in-plane of layer n
  • Ezln is the elastic modulus, out of plane of layer n
  • tln is the thickness of layer n
  • tpcb is the overall thickness of the PCB

Calculating mechanical properties of the printed circuit board takes the general forms.

 17 Equations Figure 17

The copper layers are computed assuming a combination of copper and unreinforced resin rich regions. This is a critical aspect of the CTE calculation especially when the copper weight increases above 2 oz.

This is just the beginning of the validation formulas that are integral to Sherlock. In Part II of this white paper DfR will address flex cracking of components, cyclic bending, Conductive Anodic Filament (CAF) and IC Wearout.

So, as previously stated, there are numerous equations that provide reliability validation. Look for the 2nd white paper to illustrate the rest that are currently in DfR’s Sherlock software and provide the basis for a Physics of Failure (PoF) approach to reliability assessments.

Don’t hesitate to contact DfR with any questions 


Plated Through Hole

  1. (http://www.slate.com/blogs/business_insider/2014/03/12/mathematical_equations_17_th at_changed_the_world.html
  2. IPC-TR-579, Round Robin Reliability Evaluation of Small Diameter Plated-Through Holes in Printed Wiring Boards, September 1988
  3. Goyal, H. Azimi, K.P. Chong and M.J. Lii, Reliability of high aspect ratio plated through holes (PTH) for advanced printed circuit board (PCB) packages, Proc. IRPS, 129-35 (1997).
  4. Kobayashi and S. Hayashida, A study on reliability modeling for through hole cracking failure in thermal enhanced PBGA laminate, Proc. ECTC, 1658-60 (2000)
  5. Smetana, Plated through hole reliability with high temperature lead-free soldering, The Board Authority, 4, 50-64 (2002)

Solder Joint Fatigue

  1. Wolf, J.P. 1988. Soil Structure Interaction Analysis in Time Domain, Prentice-Hall, Englewood Cliffs, NJ.
  2. Gazetas, G., “Formulas and Charts for Impedances of Surface and Embedded Foundations,” J. Geotech. Engng., 1991, ASCE, 117(9), 1363-1381.
  3. Engelmaier, W., "Chap. 17: Solder Attachment Reliability, Accelerated Testing, and Result Evaluation" in Solder Joint Reliability - Theory and Applications, edited by Lau, J. H., Van Nostrand Reinhold, New York, 1991, pp. 545-587.
  4. "Guidelines for Accelerated Reliability Testing of Surface Mount Solder Attachments," IPC-SM785, Institute for Interconnecting and Packaging Electronic Circuits, Lincolnwood, IL, July 1992.
  5. Jih, E. and Pao, Y. H., June 1995 “Evaluation of Design Parameters for Leadless Chip Resistors Transactions of the ASME Journal of Electronic Packaging, pp. 94-99, Vol. 117.
  6. Syed, A., “Accumulated Creep Strain and Energy Density Based Thermal Fatigue Life Prediction Models for SnAgCu Solder Joints,” ECTC 2004, pp. 737-746 - corrected.
  7. Qian Zhang, Abhijit Dasgupta, Dave Nelson, and Hector Pallavicini, “Systematic Study on Thermo-Mechanical Durability of Pb-Free Assemblies: Experiments and FE Analysis” Journal of Electronic Packaging, December 2005, Vol. 127, pp. 415 – 429
  8. Suhling, et. al., “Thermal cycling reliability of lead free chip resistor solder joints”,
  9. Soldering and SMT, vol. 16, no. 2, pp. 77–87, Jun. 2004.
  10. Woodrow, “Reliability and Leachate Testing of Lead-Free Solder Joints”, IPC
  11. Swan, et. al. ,”Development of Lead-Free Peripheral Leaded and PBGA Components to Meet MSL3 at 260° C Peak Reflow Profile”, IPC APEX 2001
  12. Unknown, RoHS Readiness, June 2004, update, Web-based
  13. Qi, et al., “Temperature profile effects in accelerated thermal cycling of SnPb and Pbfree solder joints”, Microelectronics Reliability (2005)
  14. Qi, et. al., ”Accelerated Thermal Fatigue of Lead-Free Solder Joints as a Function of Reflow Cooling Rate”, Journal of ELECTRONIC MATERIALS, Vol. 33, No. 12, 2004 9000 Virginia Manor Rd Ste 290, Beltsville MD 20705 | Phone: (301) 474-0607 | Fax: (866) 247-9457 | www.dfrsolutions.com
  15. Qi, et. al., ”Accelerated Thermal Cycling of Tin-Lead and Lead-Free Solder Joints”
  16. Dusek et. al., “Compatibility of Lead-Free Alloys with Current PCB Materials”, IMAPS 2002, pp. 110-115
  17. Dusek, et. al., “Effect of PCB Finish, Processing and Microstructure on Lead-Free Solder Joint Reliability”, NPL Report, September 2005
  18. Schubert et. al., “Lead-free Solder Interconnects: Characterization testing and Reliability”, Proceedings of the 3rd International Conference on Benefiting from Thermal and Mechanical Simulation in (micro)-Electronics (ESIME), Paris France, April 15-17, 2002, pp. 62-7

Coefficient of Expansion

  1. Woo Seok Chin, Dai Gil Lee, Laminating rule for predicting the dielectric properties of Eglass/epoxy laminate composite, Composite Structures, Vol. 77, Issue 3, Feb 07, pp 373-38 a. Density of E-glass (kg/m3), 2.54 × 103. Density of epoxy resin (kg/m3), 1.2 × 103
  2. wvu.edu/templates/showETD.cfm?recnum=3194 a. density of E-glass of 2.54 to. 2.59 g/cm3, which was obtained from ASTM D2734
  3. Smith, W. F., Principles of Materials Science and Engineering, McGraw-Hill, 1990 a. density of E-glass is 2.54