Constitutive Mechano Sorptive Creep Modeling of Multi Ply Recycled Board Delamination under Dynamic Warehouse Stacking Loads

Swell

Fiber networks fabricated from recovered pulp display pronounced dimensional instability when exposed to fluctuating relative humidity. During humidity elevation, cellulose microfibrils absorb water molecules into their amorphous regions, pushing hydrogen-bonded chains apart. In a multi-ply recycled board structure, this moisture-driven expansion varies dramatically across the caliper.

The outer plies, typically manufactured from bleached virgin kraft or high-grade deinked pulp, feature tighter fiber orientation distributions and lower hygro-expansion coefficients. The core plies, composed of coarse mechanically pulped post-consumer fibers and filler stock, possess lower initial density, higher equilibrium moisture content, and significantly greater out-of-plane expansion potential.

Hygro-Expansion Anisotropy across Recycled Plies

Differential strain between adjacent board layers generates severe internal shear stresses long before external stacking forces are applied. Fiber orientation profiles established at the paper machine wet end dictate the principal axes of this expansion. The machine direction exhibits minimal length change during moisture uptake, whereas the cross-machine direction expands by factors three to five times larger.

Z-directional thickness increase is even more severe, often exceeding ten percent during a shift from fifty percent to ninety percent relative humidity. When these dimensional shifts occur unevenly through the sheet caliper, the boundaries separating internal plies absorb the resulting shear strain.

Starch adhesive interfaces applied during multi-ply web combination must transmit these hygro-expansion forces across the ply boundaries. Recycled fiber plies present reduced specific surface area for hydrogen bonding due to fiber hornification, a permanent loss of swelling capacity caused by repeated drying cycles. Converting operations compensate by applying native or cationic starches at the ply bond station.

While starch films bond well under dry ambient conditions, elevated moisture plastifies the starch matrix, lowering its glass transition temperature and reducing its elastic modulus. The combined action of fiber swelling and starch softening initiates micro-cracking within the interfacial bond zone.

Transient Moisture Gradients and Starch Adhesion Breakdown

Transient environmental cycles create steep moisture gradients through the board thickness. Relative humidity shifts do not penetrate the internal plies instantaneously. Water vapor diffuses inward from the exposed top and bottom surfaces, producing a moving front of elevated moisture content.

As the outer layers swell while the interior core remains dry, internal tensile stresses develop in the plane of the sheet, balanced by compression in the core. When the ambient air turns dry, the process reverses: outer plies contract rapidly, placing the core in Z-directional tension. Repeated humidity cycling drives mechanical fatigue directly into these adhesive planes.

Cross-directional hygro-expansion coefficients in recycled middle plies reach double the magnitude observed in virgin kraft top liners during humidity transients.

Substrate suppliers frequently defend stack failures under dynamic humidity by asserting that total board caliper remained within contract tolerances. This position ignores the transient strain history developed within internal ply boundaries. A board lot can pass standard room-temperature caliper audits while hiding severe micro-fractures along its internal starch lines.

The structural breakdown stems from cyclic moisture gradients rather than static dimensional overload.

Matrix

Constitutive modeling of mechano-sorptive creep requires a continuum mechanics formulation that couples viscoelastic stress relaxation with moisture-driven mechanical strain acceleration. Standard linear viscoelastic models predict that creep deformation under constant mechanical load follows a smooth asymptotic curve toward equilibrium. Paper and board substrates break this expectation completely when subjected to transient relative humidity.

Under changing moisture conditions, the rate of deformation accelerates beyond the sum of independent mechanical creep and hygro-expansion components. This non-linear phenomenon, known as the mechano-sorptive effect, accelerates structural collapse in stacked corrugated and solid board packaging.

Formulating Mechano-Sorptive Creep Tensors

Mathematical representation of the orthotropic strain tensor breaks total strain rate into four distinct constituent components. Total strain rate equals the elastic strain rate plus the linear viscoelastic strain rate plus the pure hygro-expansion strain rate plus the mechano-sorptive strain rate. The mechano-sorptive strain rate tensor depends directly on the absolute value of the moisture content change rate, multiplied by a material constitutive tensor and the current stress tensor.

This absolute value dependency creates asymmetric creep behavior: both wetting and drying cycles induce irreversible plastic deformation in the fiber network.

The constitutive equations for multi-ply recycled board incorporate distinct parameters for each ply layer due to varying fiber chemistry and recycled content ratios. Recycled fibers exhibit higher mechano-sorptive compliance than virgin fibers because hydrogen bonds in recycled networks undergo rapid breakage and reformation under stress during moisture transport. The constitutive relation takes the form:

dE_ms / dt = M |dM / dt| S

where dE_ms / dt represents the mechano-sorptive strain rate tensor, M is the dimensionless mechano-sorptive material matrix, dM / dt is the rate of moisture content change within the ply layer, and S is the Cauchy stress tensor. The material matrix M contains coupling coefficients that account for fiber alignment anisotropy across machine, cross, and thickness directions.

Stress Integration across Laminate Interfaces

Numerical implementation of this constitutive law within finite element codes requires implicit stress integration algorithms to maintain numerical stability during rapid humidity transitions. Plane stress assumptions valid for thin paper sheets fail when applied to thick multi-ply boards subjected to stacking loads. A full three-dimensional stress state must be integrated at each material point, tracking out-of-plane shear stresses and Z-directional tensile stresses at the ply boundaries.

Stacking loads impose multi-axial stress states on the lower containers in a warehouse stack. The compressive load applied to box corners creates vertical compression combined with lateral bulginess and out-of-plane shear along the vertical creases. When relative humidity cycles between fifty percent and ninety percent in twelve-hour storage facility cycles, the stress tensor within the internal plies rotates continuously.

This rotating stress field accelerates shear slip between starch-bonded fiber layers.

Transient humidity changes drive mechano-sorptive deformation regardless of whether the paperboard is absorbing or desorbing moisture.

Interlayer shear stresses concentrate near free edges and score lines. The constitutive matrix reveals that mechanical creep compliance increases by an order of magnitude when moisture flux exceeds two percent moisture content change per hour. Packaging lines using recycled board must account for this compliance shift during structural stack calculations.

• Interfacial Shear Discontinuities emerge at boundaries where adjacent plies possess mismatched mechano-sorptive compliance coefficients.
• Transient Compressive Buckling occurs when outer ply contraction forces internal core layers into localized out-of-plane bending.
• Plastic Strain Accumulation builds progressively with each humidity cycle, preventing full elastic recovery of box height.
• Stress Relaxation Loss reduces the load-carrying capacity of vertical box corners during prolonged high-humidity dwell periods.

Recycled board plies experiencing rapid moisture cycling lose structural stiffness faster than virgin grades of identical grammage.

Fracture

Delamination in multi-ply recycled packaging boards follows cohesive zone fracture mechanics principles. Under dynamic stacking loads and ambient humidity cycles, interface failure initiates at microscopic void spaces in the starch adhesive layer. Cohesive Zone Modeling (CZM) simulates this progressive damage by defining a traction-separation law along the internal ply interfaces.

The interface transfers stress linearly up to a critical penalty stiffness threshold, after which damage accumulation degrades the stiffness tensor until complete separation occurs.

Interlaminar Traction Separation Law Calibration

Calibration of the bi-linear traction-separation curve for recycled board plies requires measuring Mode I normal opening fracture energy and Mode II shear fracture energy. Mode I fracture energy governs Z-directional tensile splitting, while Mode II fracture energy governs interlaminar shear failure under compressive side wall bulging. Standard Scott Bond testing provides an empirical measure of total energy absorption during rapid z-direction rupture, but fails to isolate the critical displacement parameters needed for constitutive finite element modeling.

Double Cantilever Beam (DCB) and End-Notched Flexure (ENF) tests performed inside climate-controlled test chambers establish true constitutive fracture parameters. Test data shows that moisture content elevation drastically reduces both peak interfacial shear strength and critical strain energy release rates. At fifty percent relative humidity, a standard recycled middle-ply interface exhibits a Mode I fracture energy of 140 Joules per square meter.

Elevating relative humidity to ninety percent drops this value below 60 Joules per square meter, making the interface vulnerable to low-amplitude dynamic stress pulses.

Scott Bond Degradation Kinetics under Mechanical Strain

Dynamic warehouse stacking subjects bottom-tier packaging to cyclic mechanical excitation while static load rests on the structural corners. Vibrations transmitted through warehouse floors from heavy forklift traffic create periodic vertical acceleration peaks ranging from 0.25 g to 0.80 g at frequencies between 5 Hertz and 30 Hertz. These mechanical pulses impose high-frequency dynamic shear stresses across the board’s internal interfaces.

Consider a mathematical comparison between pure viscoelastic creep failure and mechano-sorptive dynamic delamination for a standard three-tier pallet stack loading condition:

Assume a static compressive stress of 1.2 Megapascals applied to the corner posts of a bottom-tier recycled board box. Under static conditioning at fifty percent relative humidity, linear viscoelastic creep produces a micro-strain rate of 1.5 micro-strain per hour. The interface retains ninety percent of its initial Mode I fracture energy after 500 hours, resulting in zero structural delamination.

Now introduce cyclic environmental conditions (relative humidity cycling between fifty percent and ninety percent every eight hours) combined with a low-amplitude mechanical vibration of 0.35 g at 12 Hertz. The mechano-sorptive strain rate accelerates to 48 micro-strain per hour during the wetting phases. The combined mechanical fatigue and moisture-driven softening degrade the interfacial damage parameter D from 0.0 to 1.0 within 72 hours.

Micro-cracks coalesce along the middle ply starch interface, causing total z-direction delamination and catastrophic corner post strain failure.

Mode I strain energy release rates drop by over sixty percent when board equilibrium moisture shifts from standard laboratory condition to high-bay storage ambient.

Ignoring mechano-sorptive fatigue parameters during box strength specification leads directly to premature stack collapse under real storage transit conditions. Standard Box Compression Test (BCT) formulas like McKee’s equation overestimate pallet stack life because they rely on static short-span compressive strength measured under uniform dry atmospheres.

Dynamics

Warehouse operations generate a complex dynamic stress environment that combines static payload mass, dynamic floor vibrations, and ambient environmental transport cycles. Automatic Guided Vehicles (AGVs), reach trucks, and high-bay conveyor systems induce mechanical oscillations that propagate through stacked pallets. When these dynamic forces interact with mechano-sorptive ply softening, failure modes transition from simple vertical crease crushing to extensive internal ply separation.

Can Dynamic Stacking Frequencies Accelerate Interlayer Rupture?

Resonance phenomena within multi-ply paperboard stacks compound structural degradation rates. A full pallet stack acts as a multi-degree-of-freedom spring-mass-damper system. The vertical stiffness of the system depends directly on the out-of-plane elastic modulus of the individual packaging boxes.

As moisture content cycles upward, the Z-directional modulus of recycled core plies drops, reducing the natural frequency of the entire pallet stack.

If the natural frequency of the stack degrades until it aligns with dominant warehouse floor vibration frequencies (typically 8 Hertz to 15 Hertz), resonance occurs. Dynamic amplification factors can double or triple the effective vertical compressive load acting on bottom-tier packaging. Under these resonant conditions, internal shear stresses at the starch interface exceed the moisture-degraded Mode II shear strength, initiating sudden delamination fronts that spread from box corners toward the panel centers.

Superposition of Vibration Harmonics and Sorptive Fatigue

The interaction between dynamic vibration frequencies and moisture absorption kinetics can be evaluated using a standardized qualification procedure for high-bay storage containers.

1. Mount fully loaded representative packaging units onto a multi-axis servo-hydraulic vibration table inside an environmental test chamber.
2. Condition the test chamber to fifty percent relative humidity at 23 degrees Celsius until package mass stabilizes.
3. Apply a static pre-load equivalent to maximum specified pallet stack height minus one tier payload mass.
4. Initiate sinusoidal vibration sweeps from 2 Hertz to 50 Hertz at a constant acceleration level of 0.5 g to identify initial stack natural frequencies.
5. Begin relative humidity cycling between fifty percent and ninety percent RH on a four-hour ramp schedule while maintaining continuous vibration at the identified natural frequency.
6. Monitor vertical displacement, structural stiffness loss, and acoustic emission signals to record exact time-to-delamination events.

Standard static box compression testing overstates real warehouse load capacity by up to forty percent when dynamic floor vibrations coincide with humidity shifts.

Unresolved questions remain regarding how multi-frequency dynamic excitation interacts with non-linear mechano-sorptive compliance matrices across varied recycled fiber pulp blends. Current finite element models struggle to predict exact crack initiation times when dynamic force amplitude varies stochastically alongside ambient relative humidity shifts.

Audit

Quality assurance protocols for recycled paperboard must move beyond basic grammage and dry burst strength testing to verify long-term resistance to mechano-sorptive delamination. Traditional acceptance testing relies on TAPPI T 541 for Z-directional tensile strength (Scott Bond) under static room conditions (23 degrees Celsius, 50 percent relative humidity). This snapshot measurement fails to expose structural vulnerabilities that emerge under real-world warehouse stacking conditions.

Standardized Creep Testing under Cyclic Relative Humidity

Comprehensive lot qualification requires dynamic laboratory testing that simulates combined mechanical load and environmental cycling. ISO 187 defines standard atmosphere conditioning, but advanced packaging audits demand transient dynamic conditioning per ASTM D4332 modified for cyclic relative humidity ramps. Test apparatus must apply continuous static compression to test specimens while varying relative humidity between forty percent and ninety percent every six hours.

Dynamic Mechanical Analysis (DMA) instruments equipped with relative humidity generators measure the complex storage modulus and loss modulus of individual board plies over time. A sharp drop in storage modulus during the first wetting ramp signals high sensitivity to mechano-sorptive softening. Boards showing a loss tangent (tan delta) spike above 0.15 during moisture absorption contain weakly bonded recycled fibers prone to rapid creep deformation.

Quality Verification Criteria for Recycled Boxboard Acceptance

Establishing clear technical metrics within procurement contracts safeguards buyers against sub-standard board deliveries that risk stack failure.

• Interfacial Adhesive Coverage must exceed ninety-five percent across all ply interfaces when evaluated by ultrasonic micro-tomography.
• Maximum Permitted Creep Index limits total Z-directional thickness reduction to less than eight percent after forty-eight hours under static load at high humidity.
• Delamination Resistance Factor requires retention of at least seventy percent of baseline Scott Bond strength after five consecutive moisture condensation cycles.
• Starch Gelatinization Uniformity verified through iodine staining analysis across cross-sectional microtome cuts of the unbonded sheet.

Including explicit cyclic humidity creep limits in procurement contracts shifts financial liability for stack collapse back to board mills delivering low-bond recycled grades.

Contractual specifications must mandate that board lots failing the cyclic relative humidity creep threshold trigger immediate rejection of the delivered run. Rejection clauses based strictly on dry burst strength fail to protect buyers against real warehouse delamination events.

Settlement

Mitigating delamination risk in automated high-bay warehouses requires balancing recycled fiber content against structural insurance costs. Incorporating high percentages of post-consumer recycled pulp reduces material input costs and aligns with corporate sustainability targets. Uncontrolled mechano-sorptive failure, however, incurs massive financial consequences including damaged inventory, automated warehouse downtime, safety hazard remediation, and customer penalty fees.

Commercial Cost Implications of Stack Safety Margins

Engineers historically compensated for recycled fiber degradation by over-specifying board grammage. Increasing sheet caliper by twenty percent adds direct material cost to every manufactured box while increasing shipping tare weight. Constitutive modeling demonstrates that simply increasing grammage does not prevent mechano-sorptive delamination if the internal starch bond remains sensitive to moisture cycling.

A heavy board with poor ply adhesion fails via delamination almost as fast as a lighter board with identical bond chemistry.

Targeted chemical modification of ply interfaces yields higher commercial returns than global grammage increases. Applying wet-strength resin additives (such as polyamide-epichlorohydrin) or cross-linked bio-starches at the ply bonding station increases Mode I fracture energy under humid conditions by up to eighty percent. This chemical intervention adds approximately three to five percent to raw substrate cost while allowing conversion engineers to maintain lower overall sheet grammage.

Recycled Content Penalties in Automated High-Bay Storage

Automated storage and retrieval systems (ASRS) enforce tight dimensional tolerance limits on stored pallet loads. Deflections exceeding fifteen millimeters at the top of a stacked pallet break optical safety curtains, triggering automated shut-downs of entire aisle cranes. Mechano-sorptive creep in bottom-tier boxes leads directly to structural lean, causing palleted loads to tip beyond allowable operational envelopes.

Chemical fortification of internal ply starch lines yields greater stack stability per dollar spent than indiscriminately increasing overall paperboard grammage.

Logistics operators routinely apply safety safety factors when establishing maximum stack heights for recycled packaging containers. Where virgin kraft boxes utilize a stack safety factor of 2.0 under standard storage conditions, recycled board configurations often require safety factors between 3.0 and 3.5 to account for unmodeled mechano-sorptive creep. This requirement reduces warehouse volumetric efficiency by limiting vertical pallet tiering, directly raising long-term storage overhead per stock keeping unit.

Substrate selection, adhesive formulation, and dynamic environmental controls must operate as an integrated system. Converting facilities that master constitutive modeling of multi-ply recycled board performance can optimize starch application rates, target wet-strength additive doses precisely where interfaces experience peak shear, and supply guaranteed storage ratings to warehouse operators.