Chemithermomechanical Core Viscoelastic Compliance Modeling for Packaging Substrates

Strata

Multilayer paperboard achieves flexural rigidity through spatial separation of high-modulus outer skins separated by a low-density, high-bulk central layer. Bleached chemical pulp outer plies carry in-plane tension and compression during bending. Chemithermomechanical pulp placed in the middle ply acts as a structural spacer, resisting transverse shear while minimizing total mass per unit area.

Sourcing folding boxboard for pharmaceutical, cosmetic, or food packaging balances this stiffness efficiency against time-dependent deformation under sustained loads.

Lignin retention in chemithermomechanical processing preserves the native cellular structure of wood tracheids, yielding a bulk value reaching 1.8 to 2.2 cm³/g compared to 1.1 to 1.3 cm³/g for solid bleached sulphate boards under ISO 534 testing. The stiff, uncollapsed fibres form an open network with low relative bonded area. Under rapid mechanical loading, this network exhibits high elastic compliance.

When subjected to extended duration static stresses during warehouse stacking or high-speed automated carton erection, the amorphous polymers in the lignin-hemicellulose matrix begin to yield, initiating viscoelastic creep.

A twenty percent increase in middle-ply bulk reduces total fibre furnish costs while lowering long-term compressive resistance by nearly a third.

Board structures operating in multi-tier pallet arrangements experience continuous out-of-plane and in-plane stresses. Elastic beam equations fail to predict carton sidewall deflection over time because they assume constant stiffness. In a folding boxboard laminate, the effective bending stiffness decays as a function of duration, temperature, and ambient moisture.

Engineers sizing caliper for secondary packaging must determine both instantaneous elastic compliance and time-dependent compliance parameters to prevent container collapse under distribution conditions.

Constitutive Core Mechanics

The middle ply behaves as an orthotropic, porous continuum where mechanical properties vary along the machine direction, cross direction, and out-of-plane thickness direction. Chemithermomechanical pulp networks demonstrate significant directional anisotropy arising from fibre alignment during forming on the multiply fourdrinier or cylinder vat wire. Machine direction elastic modulus routinely exceeds cross direction modulus by a factor of 2.0 to 3.2, measured under ISO 1924-2 standard atmospheric conditioning at 23 °C and 50 percent relative humidity.

Transverse shear modulus governs the shear deformation of the thick core during out-of-plane bending. Thin board formulations follow classical Euler-Bernoulli beam kinematics where shear deformation remains negligible. In packaging boards exceeding 350 µm in caliper, thick-plate Timoshenko kinematic models become necessary to account for out-of-plane shear compliance.

The effective bending stiffness per unit width follows an expanded formulation incorporating core shear:

S_b = (E_skin (h^3 – h_core^3) / 12) + (E_core h_core^3 / 12)

When shear compliance is included, the total deflection under three-point loading increases by a shear correction factor directly tied to core transverse shear modulus G_xz and core thickness h_core. Low transverse shear stiffness in loose chemithermomechanical plies degrades the overall flexural resistance during score line folding and side-seam gluer compression.

Fibre stiffness dominates short-term resistance, while the amorphous lignin binder controls long-term structural integrity.

Relaxation

Polymers inside the middle ply respond to sustained stress through molecular rearrangement. Linear viscoelastic compliance models characterize this behavior across temporal scales spanning milliseconds during die-cutting to months in storage racking. The creep compliance function D(t), defined as the time-dependent strain divided by constant stress, captures both immediate elastic response and delayed viscous elongation.

Generalized Maxwell-Wiechert and Kelvin-Voigt models provide mathematical representations of this response. A discrete Prony series expansion calculates the transient compliance through a sum of decaying exponential terms paired with characteristic relaxation times:

D(t) = D_0 + Σ + t / η_0

D_0 represents the instantaneous elastic compliance. Coefficients D_i assign compliance weights to discrete retardation time constants τ_i, while η_0 defines long-term steady-state Newtonian flow viscosity. For chemithermomechanical pulp networks, five to seven Kelvin-Voigt elements spaced logarithmically across time decades accurately fit dynamic mechanical analysis data across frequency sweeps from 0.01 Hz to 100 Hz.

Time-Temperature Superposition Principles

Thermal acceleration alters the relaxation kinetics of paperboard constituents. Hemicellulose exhibits an apparent glass transition between 60 °C and 90 °C depending on water content, while native unaltered lignin softens between 120 °C and 150 °C under dry conditions. Moisture plasticization depresses these thermal transitions toward ambient room conditions.

The time-temperature superposition principle enables construction of a single master compliance curve over extended temporal horizons from short-term mechanical tests conducted at elevated temperatures.

The Williams-Landel-Ferry relationship calculates the horizontal logarithmic shift factor a_T for operational ranges above the glass transition temperature:

log10(a_T) = -C1 (T – T_ref) / (C2 + (T – T_ref))

Empirical fit parameters C1 and C2 for softwood-derived chemithermomechanical pulp cores differ markedly from petroleum polymers. Below the thermal transition, Arrhenius activation energy formulations describe the shift behavior with activation energies typically falling between 80 kJ/mol and 140 kJ/mol for unbleached high-yield fractions.

• Instantaneous glass compliance bounds the immediate elastic response under high-speed scoring anvils operating within ten-millisecond dwell windows.
• Retardation spectra distribution dictates the progressive load shedding from the compliant chemithermomechanical middle layer onto the stiff chemical pulp surface liners.
• Steady-state secondary creep governs irreversible carton bulging during long-term storage in temperature-controlled pharmaceutical distribution centres.
• Tertiary rupture threshold identifies the critical compliance inflection point where micro-buckling of individual core fibre walls initiates total panel delamination.

How do altered chemical pre-treatments in high-yield pulping shift the fundamental retardation spectrum of the finished board matrix?

Swell

Atmospheric moisture alters paperboard compliance. Cellulose microfibrils and hemicellulosic binders absorb ambient water molecules, disrupting inter-fibre hydrogen bonds. As relative humidity rises from 30 percent to 90 percent, equilibrium moisture content within a chemithermomechanical core climbs from 5.5 percent to over 13.0 percent, measured under ISO 187 conditioning protocols.

Water acts as an efficient plasticizer, lowering the storage modulus and accelerating viscoelastic relaxation across every time decade.

Under cyclic humidity conditions, the rate of deformation accelerates beyond values predicted by static high-humidity tests. This mechano-sorptive creep phenomenon induces rapid structural degradation under non-stationary environmental transit. Moisture sorption gradients create localized internal stresses that drive micro-scale sliding between uncollapsed wood fibres in the high-yield furnish.

Hygro-expansion strains further distort package geometry during ambient transitions. Transverse swelling coefficients in chemithermomechanical plies exceed longitudinal values by an order of magnitude due to radial fibre expansion. Differential swelling between the chemical pulp outer liners and the high-yield core generates out-of-plane panel warpage, altering the eccentricity of vertical stacking loads.

Under ISO 535 testing, a water absorption Cobb60 value above 45 g/m² accelerates transient shear compliance failure along carton creases.

Mathematical modeling of mechano-sorptive compliance requires the addition of a sorption-rate-dependent strain tensor. The total strain increment combines elastic, viscoelastic, hygro-expansion, and mechano-sorptive components:

dε_total = dε_elastic + dε_viscoelastic + β dM + μ |dM| σ

Parameter β represents the hygro-expansion coefficient, M denotes moisture ratio, μ represents the mechano-sorptive compliance coupling factor, and σ is the applied stress tensor. Omitting the mechano-sorptive term causes package life prediction models to overestimate warehouse stacking survival times by factors of two to four in maritime shipping lanes.

Ignoring hygro-expansion gradients during board specification leads to unrecoverable panel buckling and automatic sorting line stoppages in automated fulfillment hubs.

Stack

Predicting top-load box compression test performance over distribution cycles requires translating ply-level compliance tensors into structural container mechanics. The classical McKee formula calculates short-term ultimate peak load based on caliper, edgewise compression resistance (ISO 3037), and flexural stiffness (ISO 5628). This static formulation omits all time-dependent compliance parameters, creating dangerous blind spots for supply chain qualification.

Under continuous static pallet stacking loads, the vertical sidewalls of a corrugated or folding carton behave as thin plates under uniform edge compression. Stress redistribution occurs continuously across the panel width. The central region of the panel sheds stress toward the stiffer corner creases as out-of-plane deflection progresses.

The high-yield core plies at the carton corners experience intense transverse shear stresses alongside localized axial compression.

Worked Stacking Life Construction

Evaluate a standardized retail carton fabricated from 350 g/m² folding boxboard comprising a 220 µm chemithermomechanical core and two 50 µm solid bleached kraft outer liners. Total board caliper equals 320 µm. The package geometry measures 150 mm in length, 100 mm in width, and 200 mm in height.

Top-load static palletization imposes a continuous dead load of 180 N per carton at 23 °C and 75 percent relative humidity.

Instantaneous board properties determine initial baseline states. The edgewise compression resistance is 4.2 kN/m in the cross direction. Flexural stiffness in the machine direction equals 24.5 mN·m, while cross-direction stiffness equals 11.2 mN·m, yielding a geometric mean stiffness of 16.6 mN·m.

The short-term McKee static failure load evaluates as:

BCT_static = 5.876 CD_ECT^0.746 (MD_stiffness CD_stiffness)^0.127 Perimeter^0.492

Perimeter equals 500 mm. Substituting values produces a static failure capacity of 585 N. The initial applied stress ratio equals 180 N / 585 N, or 0.308. Under purely elastic assumptions, a safety factor of 3.25 suggests indefinite stacking survival.

Viscoelastic compliance modeling alters this calculation. The transient core compliance at 75 percent relative humidity causes geometric mean bending stiffness to degrade over time according to the Prony series fit:

S_b(t) = S_b0

Time t is expressed in hours. Concurrently, out-of-plane panel creep deflection increases following the viscoelastic compliance tensor. Critical elastic-viscoelastic bifurcation buckling occurs when out-of-plane creep displacement at panel center reaches 1.5 times the total board caliper, corresponding to 480 µm.

1. Hour 0 marks instantaneous elastic loading, creating an initial central panel deflection of 85 µm with zero structural damage.
2. Hour 24 exhibits progressive core compliance relaxation where bending stiffness drops to 14.1 mN·m, increasing central panel deflection to 195 µm.
3. Hour 120 shows secondary creep acceleration as middle-ply shear compliance drives out-of-plane displacement to 340 µm, transferring forty percent of face stresses onto the corners.
4. Hour 288 reaches the critical stability threshold where center deflection exceeds 480 µm, triggering sudden corner delamination and sidewall buckling at 12 days of warehouse dwell.

Paper mills often claim that transient laboratory humidity spikes do not compromise carton integrity once atmospheric equilibrium returns.

Yield

Switching from solid bleached sulphate board to a chemithermomechanical-cored folding boxboard alters sheet economics and converting efficiency. High-yield core fibres deliver equal bending stiffness at lower basis weights. A 300 g/m² folding boxboard matches the flexural rigidity of a 350 g/m² solid bleached sulphate sheet, creating a 14.3 percent reduction in total furnish tonnage requirements for identical carton specifications.

Yield calculations govern landed packaging unit economics. Buying by the metric ton and converting into individual cartons means area per ton dictates actual unit cost. A sheeted lot priced at 1,450 EUR per metric ton for 300 g/m² folding boxboard provides 3,333 m² of usable substrate per ton.

Solid bleached sulphate at 350 g/m² priced at 1,520 EUR per ton yields only 2,857 m² per ton. The raw substrate cost per thousand 0.15 m² carton blanks drops from 114.00 EUR to 65.25 EUR.

These material savings carry downstream converting constraints. Lower z-directional internal bond strength (ISO 16260) in chemithermomechanical layers limits high-speed die-cutting performance. Creasing rules must be configured with wider female channel profiles to prevent delamination during 180-degree pre-breaking on multi-point gluing machines.

A minimum Scott bond value of 130 J/m² maintains score line integrity during high-speed packaging line erection.

Tooling geometry requires adjustments when processing high-bulk substrates. Deeper creasing scores compress the bulky middle ply, permanently breaking hydrogen bonds within the high-yield furnish. Over-creasing induces localized viscoelastic fracture, creating spring-back resistance variations that jam cartoner feeder arms.

Procurement agreements specify DIN 55437-1 testing for score line bending resistance to ensure converting lines run without torque overloads.