Predicting Dimensional Shrinkage Anisotropy in Sub-Core Laminates during Repeated High Temperature Lamination Passes
Predict anisotropic sub-core shrinkage by coupling glass yarn orientation with etched copper density across logarithmic decay models for each thermal pass.
Grain
Woven reinforcement in thin copper-clad laminates sets up a directional stiffness difference that drives all downstream thermal movement. Continuous-filament E-glass fabrics have different yarn counts, twist geometries, and mechanical tensions between the warp and weft directions. Warp yarns are held under continuous mechanical tension through scouring, silane finishing, and treater impregnation passes.
Weft yarns, shot across the loom during weaving, retain higher crimp and lower axial tension. Once the thermosetting resin climbs past its glass transition temperature during lamination, the lower modulus along the fill yarns allows more compressive relaxation than the taut warp can accommodate.
Sub-cores below 100 micrometres show especially sharp shrinkage asymmetry. Because the volumetric ratio of resin to glass runs much higher than in standard rigid cores, uncured or partially cross-linked polymer matrices dictate bulk dimensional behavior. As temperatures cross the glass transition threshold, thermal expansion in the resin shears against the yarn bundles.
Then, as the stack cools back down through vitrification, volumetric contraction pulls the crimped fill yarns into tighter compaction, leaving an immediate anisotropic shrinkage.
Warp shrinkage on standard 1080 E-glass cores remains below 0.04 percent at 210 degrees Celsius, while weft shrinkage exceeds 0.09 percent under identical vacuum press parameters.
Fabric style largely dictates how wide that directional split becomes. Square-weave fabrics like style 106 and 1080 pair tight yarn spacing with low filament counts, leaving large resin-rich windows that amplify contraction across the fill. Spread-glass fabrics such as style 1078 or 1067 flatten the filaments across both axes, lowering crimp angles and distributing resin more evenly.
That added mechanical restraint narrows the gap between warp and weft movement during high-temperature lamination.
| Glass Style | Nominal Dielectric Thickness (µm) | Resin Content (% by Weight) | Warp Shrinkage Mean (ppm) | Weft Shrinkage Mean (ppm) | Anisotropy Ratio (Weft to Warp) |
|---|---|---|---|---|---|
| 106 | 38 | 72.0 | -420 | -890 | 2.12 |
| 1080 | 65 | 64.0 | -310 | -680 | 2.19 |
| 1078 Flat | 68 | 62.0 | -260 | -410 | 1.58 |
| 2116 | 110 | 54.0 | -180 | -320 | 1.78 |
| 3313 Flat | 85 | 57.0 | -210 | -360 | 1.71 |
Resin chemistry also interacts with yarn mechanics. High-temperature hydrocarbon and polyphenylene ether systems feature glass transition temperatures above 200 degrees Celsius and moisture absorption values below 0.2 percent by weight under IPC-TM-650 Method 2.6.2.1. These non-polar systems undergo volumetric polymerization shrinkage as they cure, locking baseline stresses into the yarn intersections well before the bare core ever sees circuit etching.
Etching off the sub-core foil strips away the laminate’s primary mechanical restraint. Fully etched dielectric cores contract immediately once the copper is gone, releasing residual tensile stress locked in by electrodeposited foil bonding. Unconstrained, the sheet flattens out, and subsequent thermal cycles drive uninhibited relaxation along the more compliant weft axis.
Relaxation
Sequential build cycles in high-density interconnect production repeatedly push sub-cores back past their vitrification points. A four-plus-N-plus-four buildup sequence, for example, puts the central sub-core layers through five separate hydraulic press runs. Each cycle typically dwells between 195 and 225 degrees Celsius for 90 to 120 minutes under 1.8 to 3.2 megapascals of hydraulic pressure.
While the initial press achieves near-complete cross-linking density, subsequent thermal cycles trigger secondary physical aging and irreversible viscoelastic stress relaxation through the cross-linked network.
Polymer chains caught in non-equilibrium conformations during primary pressing undergo segmental rearrangement whenever reheated past the glass transition temperature. Frozen elastic strain energy locked into the twisted glass bundles dissipates through micro-creep over time, though that reorganization slows considerably with every additional cycle.
The initial lamination pass consumes most of the available residual strain. Secondary and tertiary cycles produce progressively smaller incremental movements, tracing a non-linear decay curve: thin sub-cores undergo a sharp initial contraction during the first sequential buildup pass, followed by roughly logarithmic decay across later cycles.

What Governs Asymmetry across Sequential Thermal Cycles?
Thermal degradation of silane coupling agents at the resin-glass interface accelerates dimensional drift across multiple heat runs. Near 220 degrees Celsius, thermo-oxidative stress weakens the bonds between the organofunctional silane sizing and the filament glass surface. As interfacial shear strength drops, the resin contracts more independently from the reinforcement, widening the dimensional split between the rigid warp filaments and compliant weft yarns.
| Lamination Pass | Incremental Warp Shift (ppm) | Incremental Weft Shift (ppm) | Cumulative Warp (ppm) | Cumulative Weft (ppm) | Net Anisotropy Delta (ppm) |
|---|---|---|---|---|---|
| Pass 1 (Core Prep) | -340 | -720 | -340 | -720 | 380 |
| Pass 2 (Inner HDI) | -160 | -310 | -500 | -1030 | 530 |
| Pass 3 (Outer HDI) | -70 | -130 | -570 | -1160 | 590 |
| Pass 4 (Cap Layer) | -30 | -50 | -600 | -1210 | 610 |
Press pads also degrade under continuous thermal cycling, altering how hydraulic pressure distributes across the working panel. Silicone rubber and cellulose buffers lose compliance over time, setting up pressure gradients toward the panel edges. Panels in the center of an opening also see higher effective platen temperatures than perimeter panels, creating radial thermal variations that warp the directional shrinkage profile across the working area.
Whether prolonged heat exposure ever brings the resin-glass matrix to true mechanical saturation remains unresolved in sequential laminate processing.

Vector
Predicting sequential shrinkage reliably requires an orthotropic constitutive model rather than flat percentage scaling. The global dimensional strain vector for a thin sub-core under sequential lamination combines thermal contraction, chemical cure shrinkage, resin content, and the directional restraint of patterned copper features. Net shrinkage along either axis works out to the bare composite baseline strain adjusted for residual copper volume and the directional layout tensor.
Copper distribution creates an uneven mechanical anchor across the sub-core surface. Electrodeposited copper has an isotropic thermal expansion coefficient near 17 parts per million per degree Celsius and an elastic modulus around 110 gigapascals, whereas the resin-glass dielectric has a lower modulus and expands anisotropically. Solid copper ground planes mechanically hold the underlying dielectric in place during thermal dwell, while heavily etched routing areas allow it to pull inward unconstrained.
Laminate sheets cut against the master roll grain direction reverse the predictable shrinkage axes and invalidate all phototool compensation files.
Calculating the directional dimensional change factor requires resolving warp and weft strain independently for each successive pass:
S_x(n) = S_0x exp(-k_x (n – 1)) (1 – C_x R_cu_x)
S_y(n) = S_0y exp(-k_y (n – 1)) (1 – C_y R_cu_y)
In these equations, S_x(n) and S_y(n) represent the incremental shrinkage strain in the warp and weft directions for pass number n. S_0x and S_0y denote the initial unrestrained shrinkage coefficients of the specific glass style and resin system. The terms k_x and k_y serve as the logarithmic decay constants for mechanical relaxation, while R_cu_x and R_cu_y represent the fractional copper density aligned parallel to the respective measurement axis.
C_x and C_y denote the copper coupling rigidity factors derived from foil thickness.
As a practical example, consider an 18-layer buildup using a 50-micrometre core on style 106 high-Tg FR-4, put through three consecutive lamination passes. Assume initial bare-core shrinkage factors of -450 parts per million in warp and -920 parts per million in weft. Copper retention runs at 75 percent along the x-axis (warp) because of bus routing, but drops to 20 percent along the y-axis (weft) across dense etched clearance arrays.
Decay constants are 0.85 for warp and 0.78 for weft, with a copper coupling coefficient of 0.65 for 18-micrometre foil.
Warp shrinkage on pass one calculates to -450 (1 – 0.65 0.75), yielding -230.6 parts per million. Weft shrinkage on pass one calculates to -920 (1 – 0.65 0.20), yielding -800.4 parts per million. For pass two, the warp increment drops to -450 exp(-0.85) (1 – 0.4875), giving -98.6 parts per million.
The weft increment drops to -920 exp(-0.78) (1 – 0.13), giving -366.1 parts per million. Summing the two-pass sequence yields a cumulative warp offset of -329.2 parts per million and a cumulative weft offset of -1166.5 parts per million. The resulting anisotropic spread reaches 837.3 parts per million between axes.
Stackup symmetry governs vertical stability across multi-pass assemblies. Several structural elements control how well sub-core layers hold their shape across repeated cycles:
- Glass fabric orientation alignment locks the warp yarns of every core layer strictly parallel to the nominal panel processing direction.
- Copper weight distribution balancing pairs equivalent foil thicknesses across opposing sub-core faces to suppress out-of-plane cylindrical warpage.
- Resin volume fraction matching maintains uniform prepreg ply selections across opposing sequential buildup stages.
- Etched pattern density compensation utilizes dummy copper thieving patches in clearance windows to equalize directional mechanical restraint.
Thicker copper foils suppress composite shrinkage more effectively than thin foils across all thermal cycles.

Register
Direct imaging phototools depend on accurate scaling offsets to align outer microvia drills with recessed sub-core target pads below. A 500-millimetre panel undergoing uncompensated 800 parts per million shrinkage along its weft axis loses 400 micrometres edge-to-edge. If 65-micrometre laser-drilled blind microvias target 150-micrometre landing pads, an uncompensated runout error past 35 micrometres destroys minimum annular ring integrity under IPC-6012 Class 3 requirements.
Dynamic laser direct imaging compensates for linear expansion and contraction by reading fiducial targets on every individual panel. Modern exposure systems apply four-point Cartesian scaling or multi-zone non-linear fitting to distort the nominal Gerber artwork dynamically. While these optical systems handle uniform stretching and basic linear trapezoidal distortions, complex local shifts caused by copper imbalance or uneven sub-core grain alignment quickly exceed what dynamic algorithms can correct.
Inner-layer misregistration exceeding forty micrometres causes immediate breakout on sub-core capture pads and generates microsection rejections under IPC-A-600 Class 3 evaluation.
Inner-layer prep operations apply pre-baking cycles to stabilize thin laminates before primary circuit etching. Baking bare cores unconstrained at 150 degrees Celsius for two to four hours relieves residual sheeting stresses, removing initial mechanical hysteresis prior to photoresist lamination and chemical etching. This pre-treatment step reduces the magnitude of unpredictable movement in subsequent vacuum press cycles.
Tooling hole punch units set the mechanical registration baseline for inner-layer panelization. Post-etch punch systems read optical targets etched directly onto the core copper rather than relying on raw mechanical edges. This optical alignment step centers the registration hole pattern on the actual etched circuit geometry, neutralizing bulk anisotropic shrinkage that occurred during the etching bake cycle.
Failing to predict directional shrinkage anisotropy across repeated press passes leads directly to radial layer misregistration, annular ring breakout, torn hole walls during drilling, and scrapped high-density panels.

