Constitutive Multiphysics Modeling of Mechanosorptive Creep and Fiber Interface Hygro-Degradation in Secondary Containerboard

Mechanics

The mathematical description of creep in recycled paperboards under varying moisture states integrates transient mass transport with non-linear viscoelastic strain tensors. Secondary containerboard webs, produced predominantly from old corrugated containers, possess high structural heterogeneity, variable pore distributions, and damaged fiber walls resulting from successive repulping cycles. When exposed to dynamic environmental relative humidity, these paperboards exhibit strain rates significantly exceeding the sum of independent elastic, viscoelastic, and hygro-expansion strain components.

Moisture influx accelerates network strain. This phenomenon, known as mechanosorptive creep, undermines long-term compression stability in secondary containerboard packaging deployed across fluctuating climate zones.

Transient Moisture Transport and Stress Coupling

Water vapor diffusion through recycled paperboard webs follows non-Fickian kinetics governed by localized micro-structural swelling. Moisture transport inside the fibrous network occurs through bulk vapor phase movement within internal pore networks, surface diffusion along fiber walls, and bound water transport inside cellulosic cell structures. Localized moisture concentration shifts alter local mechanical stiffness, while applied stress fields simultaneously modify transport boundary dynamics by altering pore volume and pore tortuosity.

Mass transport modeling in secondary containerboard relies on an extended advection-diffusion equation incorporating strain-dependent diffusivity coefficients. The governing transport balance accounts for local transient moisture concentration within the sheet matrix:

dM / dt = grad( D(sigma, M) grad(M) ) + Q_sorp

Where M represents local dry-basis moisture content, D denotes the non-isotropic diffusion tensor modified by equivalent stress state sigma and moisture level M, and Q_sorp defines internal sorption kinetics between inter-fiber pores and cell wall micro-fibrils. Recycled fiber networks slip under load. Because secondary containerboard contains shorter fibers and higher micro-fines fractions than virgin kraft substrates, moisture sorption rates accelerate inside regions of high local porosity, generating steep moisture gradients across web calipers under rapid ambient transitions.

Constitutive Rheology under Dynamic Relative Humidity

Viscoelastic constitutive models extend standard Maxwell and Kelvin-Voigt element arrays by appending a mechanosorptive strain rate tensor. Total strain rate decomposition in secondary linerboard and corrugating medium under combined stress and moisture variation takes the following additive operational form:

d(epsilon) / dt = d(epsilon_e) / dt + d(epsilon_vp) / dt + d(epsilon_u) / dt + d(epsilon_ms) / dt

Elastic strain rate epsilon_e responds directly to stress updates divided by the moisture-dependent elastic modulus. Viscoelastic-viscoplastic strain rate epsilon_vp models intrinsic time-dependent rheology at constant moisture content. Unidirectional hygro-expansion rate epsilon_u reflects stress-free swelling or shrinkage driven by moisture changes.

Bound water acts as a plasticizer. The mechanosorptive strain rate tensor epsilon_ms couples applied stress directly to the absolute value of the moisture sorption rate:

d(epsilon_ms) / dt = m sigma | dM / dt |

Where m represents the constitutive mechanosorptive compliance coefficient tensor. Hornification permanently shrinks fiber pore volume. In recycled OCC furnishes, the magnitude of coefficient m increases substantially compared to virgin softwood pulp due to broken inter-fiber hydrogen bonding structures, pre-existing micro-cracks from converting operations, and variable hemicellulose distribution across damaged fiber surfaces.

Bond

Cellulosic fiber networks derive cohesive strength from hydrogen contacts at overlapping crossover points. In secondary containerboard, repulping cycles strip surface fibrils, shorten fiber lengths, and accumulate inactive inorganic filler loading. When relative humidity fluctuates, inter-fiber crossover contacts experience concentrated hygro-mechanical stresses.

Adsorbed water molecules penetrate cell walls, swell amorphous cellulose regions, and disassociate hydrogen bonds holding adjacent fibers together.

Micro-Mechanical Hygro-Degradation at Fiber Crossovers

Repeated recycling loops shorten individual wood fibers while introducing micro-fibrillar wall damage. During moisture adsorption, water molecules enter the amorphous matrix of cellulose and hemicellulose, forcing inter-chain spacing to expand. This molecular expansion induces micro-mechanical shear strains directly across the fiber-to-fiber bond area.

Fiber shortening lowers network cohesion. As moisture content rises above eight percent, hydrogen bond density inside the contact zone drops sharply, diminishing inter-fiber shear capacity.

Inter-fiber contact shear strength decreases by 42 percent when equilibrium moisture content rises from 6 percent to 14 percent at 23 degrees Celsius.

Lignin hinders swelling at interface contacts. High-yield secondary furnishes retaining residual lignin display distinct sorption kinetics compared to fully delignified kraft fibers. Non-uniform swelling between softwood fibers and recycled hardwood fines creates localized strain stress concentrations across crossover joints, prompting progressive micro-fracture under sustained stack loads.

Cohesive Zone Formulations for Secondary Furnish Interfaces

Numerical simulation of fiber separation employs traction-separation damage laws coupled to localized moisture density. A bilinear cohesive zone formulation captures both initial elastic stiffness and post-peak softening at fiber crossover joints. The maximum traction threshold drops explicitly as local moisture concentration increases:

T_max(M) = T_0 (1 – alpha (M – M_ref)) (1 – D_hygro)

Where T_0 defines pristine bond traction at reference moisture M_ref, alpha represents the moisture softening coefficient, and D_hygro tracks accumulated irreversible damage caused by cyclic sorption stress. Hydrogen bonds disassociate under vapor pressure. When local moisture concentration cycles repeatedly, D_hygro accumulates continuously, reducing the critical energy release rate required to propagate interface debonding across the secondary containerboard network.

Mill technical support frequently attributes localized wall collapse to inadequate wet-strength additive addition rather than baseline fiber degradation.

Climate

Distribution environments introduce dynamic temperature and vapor differentials during intermodal transport. Secondary containerboard boxes stacked inside ocean freight containers or unconditioned transit hubs undergo rapid relative humidity changes, driven by diurnal weather shifts and cold-chain temperature transitions. These environmental swings provoke continuous moisture adsorption and desorption within corrugated packaging structures, accelerating creep deformation far beyond levels predicted by static ambient testing.

Environmental Cycling Profiles in Cold-Chain Logistics

Standard conditioning protocols under ISO 187 specify constant 23 degrees Celsius and 50 percent relative humidity, whereas real shipping channels experience abrupt psychrometric oscillations. Relative humidity inside shipping containers regularly fluctuates between 35 percent and 92 percent over a twelve-hour period. Dew point transitions trigger liquid condensation directly onto corrugated flute tips, causing localized moisture contents to exceed 18 percent by weight.

Selecting secondary containerboard grades for variable logistics environments involves evaluating several structural parameters:

• Equilibrium moisture capacity evaluates the total water sorption potential across expected transit conditions.
• Hygral expansion coefficient dictates dimensional stability across machine and cross directions during humidity swings.
• Cyclic compliance ratio measures the factor by which transient moisture states accelerate creep strain beyond static levels.
• Internal sizing threshold determines liquid water resistance when dew point condensation strikes flute tips.

Structural Box Compression Loss in Transient Humidity

Top-to-bottom load resistance in corrugated containers decays rapidly when ambient humidity moves dynamically. Box compression performance under static load relies on the compressive strength of individual linerboard and corrugating medium plies, typically characterized by short-span compressive strength values under ISO 9863. Under transient humidity conditions, mechanosorptive creep induces wall bowing and flute flattening at loads well below nominal static limits.

Adherence to TAPPI T804 box compression parameters establishes static load limits but neglects the mechanosorptive safety factor required for maritime shipment.

Dew point shifts trigger rapid creep. Sorption cycles fatigue cellulosic bonds. As the outer linerboard ply adsorbs water, it expands hygrally while losing compressive stiffness, transferring vertical stack loads toward the inner liner ply and corrugating flute tips.

Subsequent desorption causes non-uniform structural shrinkage, trapping residual bending moments inside panel faces and triggering premature structural buckling. Section 4.2 of ISO 2233 dictates pre-conditioning regimes that alter the baseline moisture history, shifting acceptable compliance bands for secondary packaging by up to 18 percent.

Kinetics

Rate equations governing mechanosorptive strain accumulation account for time-dependent moisture sorption and applied stress vectors. In secondary containerboard, non-linear creep acceleration depends directly on the frequency and amplitude of moisture cycling. Rapid moisture fluctuations prevent the internal fiber network from reaching moisture equilibrium, generating steep moisture concentration gradients across the web caliper that amplify internal shear stresses and drive accelerated structural deformation.

Which Constitutive Parameters Govern Mechanosorptive Acceleration during Humidity Cycling?

The coupling constant linking moisture rate to strain velocity determines the magnitude of instantaneous creep acceleration. In multiphysics finite element formulations, primary constitutive parameters include the orthotropic elasticity matrix E(M), the hygro-expansion tensor beta, the viscoelastic relaxation spectrum tau_k, and the mechanosorptive compliance tensor m. In secondary OCC furnishes, the mechanosorptive parameter m exhibits pronounced stress nonlinearity when normalized compressive stress exceeds 35 percent of short-span compressive strength.

Failure mechanisms during dynamic moisture cycles progress through distinct structural stages:

• Flute wall delamination occurs when internal shear stress exceeds degraded inter-fiber bond capacity during desorption.
• Creep buckling failure strikes vertical panel walls once total non-linear strain surpasses critical stability thresholds.
• Scoreline cracking develops along fold lines where repeated swell cycles embrittle recycled fiber networks.
• Faceline blistering emerges when moisture gradients generate out-of-plane tensile forces between liner and medium.

Worked Calculation for Secondary Linerboard Deformations

Evaluating a 175 gram per square meter secondary testliner subjected to cyclic loading provides exact strain trajectory predictions. Assume a 175 gram per square meter testliner web containing 85 percent recycled OCC furnish, operating under a sustained compressive load of 2.5 MPa while relative humidity cycles between 50 percent and 85 percent every 4 hours at 23 degrees Celsius. Initial dry-basis equilibrium moisture content shifts between 6.5 percent (M_min) and 12.8 percent (M_max).

Initial elastic strain calculates directly from moisture-dependent elasticity at reference state:

epsilon_e = sigma / E(M_ref) = 2.5 MPa / 2800 MPa = 0.000893 (0.089 percent)

Viscoelastic strain after 100 hours at constant mean moisture content yields:

epsilon_vp = sigma J_0 (t / t_ref)^n = 2.5 0.000357 (100 / 1)^0.18 = 0.00204 (0.204 percent)

Mechanosorptive strain accumulation over 25 complete 4-hour cycles (100 hours total) sums the absolute moisture differentials per half-cycle, incorporating the secondary furnish compliance factor m = 0.0014 MPa^-1 per unit moisture change:

Delta_M_total = 25 2 (0.128 – 0.065) = 3.15

epsilon_ms = m sigma Delta_M_total = 0.0014 2.5 3.15 = 0.011025 (1.103 percent)

Total accumulated strain equals the sum of individual strain components:

epsilon_total = 0.000893 + 0.00204 + 0.011025 = 0.013958 (1.396 percent)

Rapid desorption strains internal networks. Cyclic sorption causes structural fatigue. Mechanosorptive strain accounts for 79 percent of total deformation under dynamic conditions, confirming that static viscoelastic models undercount cumulative creep in cycling humidity environments.

Rapid desorption cycles induce sharper internal strain gradients than slow equilibrium transitions, causing accelerated bond matrix fatigue.

Whether non-linear moisture transport models can accurately predict cumulative damage across twenty or more recycling cycles remains unresolved in current packaging science.

Specimen

Laboratory verification of containerboard performance relies on precise test geometry and controlled environmental chambers. Standard industrial quality protocols measure material strength at fixed moisture equilibrium states. Capturing mechanosorptive behavior demands dynamic multiaxial testing platforms capable of cycling environmental relative humidity while recording real-time dimensional changes and load relaxation profiles across paperboard samples.

Laboratory Calibration and Multiaxial Testing Protocols

Measuring mechanosorptive strain requires continuous axial load tracking inside environmentally regulated test enclosures. Standard short-span compressive strength tests under ISO 9863 provide instantaneous load values but fail to reveal time-dependent creep acceleration. Dynamic testing protocols combine constant tensile or compressive load applications with rapidly cycling humidity systems to isolate the mechanosorptive compliance tensor from static viscoelastic creep background effects.

Evaluating secondary containerboard reels for cold-chain or maritime logistics applications follows a sequential laboratory procedure:

1. Extract ten full-width sheet samples across the reel width immediately after unwinding the outer three protective wraps.
2. Condition test strips at 23 degrees Celsius and 50 percent relative humidity according to ISO 187 for 24 hours to standardize moisture baseline.
3. Mount test pieces in an automated tensile frame equipped with an environmental control cabinet operating at 2.5 MPa constant stress.
4. Cycle relative humidity between 30 percent and 85 percent at 2-hour intervals while logging continuous strain data for 48 hours.
5. Calculate the mechanosorptive creep coefficient by isolating moisture-driven strain rate components from static viscoelastic baseline creep.

Static compression values measured under equilibrium conditions fail to reveal the structural collapse triggered by humidity cycling.

Static tests miss dynamic field failures. Yield loss shifts mechanical performance. Laboratory procedures that evaluate containerboard solely under constant humidity conditions generate unrepresentative safety margins for packages destined for variable transit channels.

Selecting secondary linerboard based strictly on dry burst strength leads to premature box failure in humid logistics corridors.

Ledger

Financial metrics in corrugated packaging procurement depend on the yield balance between base paper weight and long-term performance under load. Secondary containerboard offers initial unit cost savings compared to virgin kraftliner substrates. Uncontrolled mechanosorptive creep introduces significant structural risks, forcing buyers to increase paper grammage or suffer product collapse penalties during ocean transit and long-term warehousing.

Commercial Sourcing and Downgrade Risk Metrics

Procurement teams evaluating secondary containerboard grades trade lower initial sheet costs against higher structural safety margins. Grammage compensation strategies frequently add ten to twenty percent higher basis weight when substituting recycled containerboard for virgin kraft, partially offsetting pulp cost savings. Grammage adjustments alter creep kinetics directly.

Optimizing substrate selection requires balancing landed raw material price against mechanosorptive safety factors. While standard testliner 3 yields a lower unit cost per tonne, its elevated mechanosorptive compliance requires increased sheet caliper to prevent box failure under transient maritime humidity. Miscalculating mechanosorptive creep allowances results in box collapse during warehouse storage, triggering uninsurable product damage claims and secondary repacking expenses across export supply networks.