Predicting Dynamic Interfacial Thermal Resistance and Viscoelastic Relaxation under Equibiaxial Stretch Gradients

Contact

Heat transfer coefficients across a pneumatic forming interface collapse by two orders of magnitude the moment a stretching thermoplastic membrane lifts away from the polished tool steel. In high-speed stretch blow moulding and heavy-gauge sheet forming, local contact pressures oscillate between 0.05 MPa and 1.20 MPa across complex corner radii. As the hot polymer sheet inflates against an aluminum or P20 tool cavity, equibiaxial extension thins the material while simultaneously flattening surface asperities.

This rapid planar extension alters the microscopic air gaps trapped between the tool surface and the polymer skin. When the polymer touches the cold cavity wall, instantaneous thermal contact resistance governs the local cooling rate, freezing molecular orientation into the part wall before viscoelastic stress relaxation completes.

Interfacial thermal conductance values span a wide experimental band under non-isothermal biaxial deformation. Static contact models routinely assume a uniform heat transfer coefficient between 1500 and 2500 W/(m²·K). Actual transient measurements during high-speed inflation reveal that conductance drops to 120 W/(m²·K) in regions where membrane tension bridges across microscopic tool toolmarks, then surges to 3800 W/(m²·K) once cavity pressure reaches peak packing levels.

The prediction of this boundary condition demands coupling the evolving surface micro-topography with the local contact pressure history.

Under an inflation pressure of 0.8 MPa against a 6061-T6 aluminum tool face, the measured interfacial heat transfer coefficient peaks at 3100 W/(m²·K) before dropping as thermal contraction pulls the frozen skin away.

Micro-asperity deformation at the boundary follows the local Deborah number of the surface layer. The characteristic relaxation time of the polymer melt competes directly with the millisecond contact duration. High strain rates increase the instantaneous yield resistance of the surface asperities, which inhibits intimate wet-out against the metal finish.

In contrast, slower inflation allows viscous flow into tool grind valleys, establishing continuous conduction paths. Tool venting directly modulates this interface by eliminating counter-pressure pockets that otherwise isolate the sheet from the chill surface.

Tooling manufacturers frequently explain corner thinning variations by asserting that vacuum draw or blow pressure simply failed to pack the outer geometry uniformly against the chill surface.

Strain

Planar extension fields generate severe thickness gradients across biaxially drawn components. In stretch-blow moulding of beverage containers or thermoforming of deep-draw industrial liners, the local draw ratio ranges from 1.5:1 on flat sidewalls to beyond 4.5:1 in spherical corner regions. Equibiaxial stretch gradients force polymer chains into high planar orientation, dramatically increasing mechanical anisotropy.

Because the true strain rate dictates the non-linear stress response of the macromolecular network, variable inflation velocities produce pronounced yield stress deviations across adjacent cavity sectors.

Constitutive modelling of this deformation relies on hyperelastic and visco-hyperelastic potential formulations. Standard Mooney-Rivlin or Ogden strain energy functions capture static equilibrium states. They fail, however, to predict the rate-dependent stress upturn observed at high extension ratios without incorporating dynamic relaxation functions.

During inflation, the rate of true strain varies continuously from 1 s⁻¹ at initiation to over 50 s⁻¹ immediately preceding cavity wall contact. High strain rates trigger strain-hardening regimes dominated by finite chain extensibility, which alters the mechanical energy stored within the relaxing network.

Local thermal history remains deeply coupled with this deformation field through strain-induced crystallization. In materials such as polyethylene terephthalate and polypropylene, rapid equibiaxial stretching above the glass transition temperature initiates nucleated crystal growth. The crystallization rate increases by several orders of magnitude under biaxial strain fields.

This crystalline development changes the specific heat capacity and thermal conductivity of the stretching membrane, dynamically modifying the heat dissipation rate into the tooling.

Non-uniform deformation introduces distinct physical failure mechanisms across the moulding:

• Corner webbing failure develops when excessive membrane tension prevents balanced lateral draw into female cavity pockets.
• Localized strain localization creates visible haze bands where excessive drawing forces necking in amorphous domains.
• Thermal boundary separation occurs when early contact freezing on convex tool features starves adjacent unchilled zones of material. The resulting wall section lacks structural rigidity.
• Molecular orientation splitting shatters structural ribs under modest impact loads. Internal residual stresses remain frozen along principal draw axes.

Whether secondary strain hardening can stabilize local thinning before micro-void coalescence triggers boundary rupture under extreme inflation speeds remains an active debate across commercial blow tool design.

Chill

Cooling circuit topology establishes the spatial dissipation rate of the moulding. Thermal management within high-cavitation tooling governs both cycle time and the final distribution of residual stresses. As molten polymer touches the chilled steel, a steep temperature gradient develops across the part thickness, often exceeding 120°C per millimeter.

This intense gradient sets up an internal skin-core morphology where the frozen surface layer retains non-equilibrium molecular extension while the molten core continues to undergo stress relaxation.

Thermal contact resistance acts as a dynamic valve throughout this freezing phase. Initially, high blow or vacuum pressure forces the soft polymer into microscopic tool irregularities, minimizing contact resistance. As the polymer cools below its crystallization or glass transition temperature, volumetric shrinkage takes place.

This shrinkage pulls the solidifying skin away from the tool steel, introducing a microscopic air gap. The thermal conductivity of this air gap is approximately 0.026 W/(m·K), which drastically reduces the local heat transfer rate and creates hot spots on thick-walled sections.

Cooling efficiency varies widely across common tool materials and interface conditions.

Transient heat transfer models must solve the coupled energy equations using time-dependent Robin boundary conditions. The heat flux across the moulding skin follows Fourier conduction modified by the interface conductance coefficient. When the tool temperature is kept too low, premature surface freezing locks in severe orientation stresses, promoting dimensional warpage after part ejection.

Conversely, elevated tool temperatures prolong cycle times, increasing landed part costs without guaranteeing improved dimensional capability.

Biaxial draw ratio directly affects thermal transport through the membrane. Stretching uncoils entangled polymer chains and aligns them parallel to the tool face, elevating in-plane thermal conductivity while reducing through-thickness conduction. In oriented polypropylene, through-plane conductivity drops to roughly half its isotropic value, slowing heat removal through the mold boundary.

Part ejection before the core cools below heat deflection temperature causes spontaneous warpage and dimensional collapse across unrestricted spans.

The sequence of thermal contact stabilization follows discrete mechanical stages during cavity inflation:

1. Initial membrane impact compresses macro-roughness peaks on the tool face within the first 15 milliseconds.
2. Pneumatic packing pressure forces molten asperities into micro-grooves, driving interface conductance to its operating maximum.
3. Solidification initiates at the skin layer, generating volumetric shrinkage forces that counter the applied cavity inflation pressure.
4. Micro-gap formation decouples the frozen skin from tool asperities, cutting heat flux by over sixty percent prior to mold opening.

Miscalculating the transient chill rate produces parts with massive frozen-in shear stresses that warp unpredictably upon secondary heat exposure, yielding high scrap rates during downstream filling and assembly operations.

Relaxation

Viscoelastic stress decay during forming controls the retained molecular orientation of the final moulding. When a thermoplastic sheet undergoes rapid biaxial drawing, entropy loss due to chain alignment generates substantial retractive forces. If contact cooling arrests molecular motion before these chains disentangle, the part retains high residual stresses.

These stresses manifest as post-moulding shrinkage, warp, and reduced environmental stress crack resistance. Predicting these phenomena requires multi-mode viscoelastic models capable of accounting for non-isothermal conditions.

Constitutive representations frequently employ the Kaye-Bernstein-Kearsley-Zapas or multi-mode upper-convected Maxwell models. The discrete relaxation spectrum represents chain dynamics through a distribution of relaxation times and relaxation moduli. The temperature dependence of these relaxation times obeys the Williams-Landel-Ferry relationship near the glass transition temperature, switching to an Arrhenius formulation at elevated melt regimes.

According to DIN 16742, dimensional tolerances on non-crystalline thermoplastic mouldings are evaluated 24 hours after production to account for ongoing viscoelastic relaxation and physical aging.

Stress relaxation proceeds at varying rates depending on local stretch ratios and instantaneous interface temperatures. In zones where high contact heat transfer chills the polymer instantly, molecular mobility freezes, locking in orientation. In thicker or poorly cooled zones, high core temperatures permit continued chain reptation, dissolving orientation and dissipating internal stresses.

This differential stress history across the part geometry generates unbalanced moments that warp the part upon ejection.

Thermal contact resistance directly dictates the boundary conditions of this relaxation equation. A high resistance maintains the polymer above its transition threshold longer, providing a wider window for stress dissipation. A sudden drop in interfacial resistance terminates relaxation abruptly.

Tool steel roughness values above 0.8 micrometers Ra trap insulating micro-voids that retard surface quenching.

A rigorous evaluation of tool feasibility hinges on comparing design options against process metrics.

Processing conditions that maximize relaxation also reduce the structural stiffness provided by molecular orientation. In container manufacturing, high orientation is necessary to achieve burst strength and gas barrier properties with minimal wall thickness. The process window must maintain orientation in critical load-bearing regions while dissipating residual stresses in sealing surfaces and base geometries.

Balancing these competing objectives requires accurate prediction of the local contact cooling history.

Rapid tool surface quenching preserves chain extension whereas delayed contact cooling permits complete entropic recovery.

Tolerance

Dimensional capability under DIN 16742 links back to the coupled effects of stretch gradients and cooling rates. Non-uniform shrinkage across equibiaxially drawn parts stems directly from localized differences in residual stress and crystallinity. In polyolefin thermoforming, unconstrained volumetric contraction reaches 1.8 to 2.5 percent, whereas heavily oriented sections shrink up to 3.5 percent in the draw direction while expanding slightly in thickness during relaxation.

Sourcing engineers face significant quote variance when toolmakers fail to account for anisotropic shrinkage in their cavity cutting coordinates.

Steel changes are expensive. Modifying an injection stretch blow tool or a multi-cavity thermoforming mold to correct post-mold warpage costs between 4,000 and 18,000 USD per cavity, depending on conformal water line placement. When differential shrinkage causes wall thickness variations outside of plus or minus 0.08 mm on high-speed filling lines, containers jam in guide rails and leak at capping stations.

Sourcing contracts must clearly specify dimensional inspection protocols, measurement holding fixtures, and environmental conditioning windows.

Process control at the press demands closed-loop monitoring of mold cavity pressure and mold wall temperature. Small shifts in incoming sheet temperature or pellet intrinsic viscosity alter the draw kinetics, shifting the peak contact point and the resulting thermal resistance profile. Setting tight tolerance limits on cooling water delivery temperature, maintaining turbulent flow in mold cooling channels with Reynolds numbers exceeding 10,000, and standardizing cycle timing guarantees stable boundary conditions across production shifts.

DIN 16742 Table 2 sets tolerance group TG4 for precision polyolefin mouldings, which allows buyers to reject production lots displaying uncontrolled warpage caused by unverified thermal contact variations.