Viscoelastic Constitutive Modeling for High Speed Converting Operations on Board

Hysteresis

High Speed Creasing Dynamics under Dynamic Strain

Creasing paperboard at high production speeds subjects the fibrous network to strain rates exceeding 1000 per second. At these velocity regimes, the transient stress response deviates sharply from static tensile profiles. The polymer components within the fiber wall, primarily amorphous hemicellulose and lignin, exhibit pronounced time-dependent relaxation.

When a male creasing rule penetrates a multilayer board into a female counter matrix within a window of 1 to 5 milliseconds, the material behaves stiffer than its static tensile modulus indicates. This rate-dependent stiffening alters how energy dissipates across the plies. A high-speed die-cutter running solid bleached sulfate (SBS) at 8000 sheets per hour creates localized shear zones that fail to develop identical internal delamination geometries compared to a platen operating at make-ready speeds.

Mechanical response during this impact phase dictates fold formation. During the strike, fiber bonds rupture while other microstructural domains deform plastically or relax viscoelastically. Converting lines operating at industrial velocity encounter springback variations caused directly by the incomplete dissipation of stored elastic energy.

The fast strain rate shifts the effective glass transition of moist hemicellulose toward higher temperatures. Ambient converting conditions of 23 degrees Celsius and 50 percent relative humidity, conforming to ISO 187, preserve an equilibrium moisture content between 6 and 8 percent. Under rapid mechanical impact, moisture redistribution through the porous network lags behind the strain wave.

Water molecules act as plasticizers in steady-state loading, yet under microsecond impacts they remain locked in local hydrogen-bonded clusters. Internal stress peaks higher before ply separation initiates.

Creasing high-density bleached kraft at web speeds above 350 meters per minute elevates apparent compressive modulus by 35 to 55 percent compared to standard laboratory tensile tests performed under ISO 1924-2.

Failure to account for dynamic strain elevation shifts the shear plane toward the outer surface plies. In high-speed flatbed autoplaten die-cutters and rotary creasing systems, this shift produces surface checking or complete liner cracking on subsequent folding. When converters increase machine velocity from 4500 to 9000 impressions per hour without recalibrating matrix depth, the unrelaxed elastic stress stored in the outer tensile face ruptures the top coating layer.

The finished carton emerges with cracked score lines, compromising barrier properties and visual presentation. Tooling configurations derived from static hand-fed presses generate reject rates exceeding four percent when applied directly to high-speed automated runs.

Relaxation

Generalized Maxwell Formulations for Board Mechanics

Mathematical representation of dynamic paperboard deformation relies on parallel Maxwell branches coupled with an elastic equilibrium spring. High-speed folding and creasing require constitutive frameworks that capture multi-decade relaxation spectrums. A standard linear solid model fails to describe the continuous decay of stress during the hold period of a rotary folder-gluer belt.

The Generalized Maxwell model, formulated as a discrete Prony series, distributes the relaxation modulus across discrete characteristic relaxation times ranging from microseconds to several minutes.

The total Cauchy stress tensor in this formulation separates into volumetric and isochoric components to isolate matrix shear deformation from bulk fiber compression. Fiber networks show significant transverse anisotropy, requiring an orthotropic formulation of each elastic element in the series. The machine direction (MD), cross direction (CD), and out-of-plane thickness direction (ZD) display distinct stiffness coefficients, where MD tensile modulus typically exceeds CD modulus by a factor of two to three, and exceeds ZD modulus by more than two orders of magnitude.

Simulating the nip dwell of a scoring wheel operating on an 800 mm diameter cylinder at 400 meters per minute highlights the functional role of the shortest relaxation times. The contact duration between the steel rule and counter matrix lasts approximately 1.2 milliseconds. During this transient event, only Branch 1 and Branch 2 engage in viscous dissipation.

The remaining branches stay frozen in their elastic configurations, generating a high peak load on the creasing anvil. Once the scored board leaves the nip and transfers to the folder-gluer infeed over an interval of 200 to 500 milliseconds, Branch 3 relaxes. This relaxation relieves a portion of the internal bending moment before the folding sword initiates panel rotation.

Are Viscoplastic Thresholds Active during Nip Transit?

Irreversible ply delamination introduces plastic deformation that operates concurrently with viscoelasticity. A purely viscoelastic Maxwell series predicts complete strain recovery over infinite time. Converting operations rely on permanent deformation inside the crease channel to establish a mechanical hinge.

Viscoplasticity formulations employ a yield surface that expands or contracts as a function of plastic strain and plastic strain rate. In paperboard substrates, an orthotropic yield criterion, such as the Tsai-Wu or Hill criteria adapted for porous fiber matrices, defines the boundary between delayed reversible elasticity and unrecoverable fiber bond failure.

Compressing the board in the thickness direction causes irreversible cell wall collapse and shear banding between individual fiber plies. When the local stress state exceeds the rate-dependent yield surface, plastic flow occurs, governed by a hardening parameter tied to fiber compaction. Creasing models must combine non-linear viscoelasticity with viscoplastic softening: the viscous damper accounts for rate-dependent peak resistance during impact, while the plastic slider retains the indentation after tooling retraction.

Incorporating non-associated flow rules prevents an overestimation of volumetric expansion when shear stresses drive internal ply shear slip.

Finite element implementations of this dual behavior update the stress state incrementally using backward Euler integration algorithms. The numerical scheme calculates an elastic trial stress, evaluates the yield condition, and executes a radial return mapping step while relaxing the viscoelastic state variables. If the solver omits viscoplastic yielding, the model overpredicts crease opening moments by up to 80 percent, misleading tooling designers regarding rule penetration depths.

The standard ISO 5628 bending resistance test records static flexural behavior that systematically underestimates the peak resisting moments generated during commercial cartoning lines running at 100,000 packs per hour.

Delamination

Ply Separation and Cohesive Zone Formulations

Multilayer paperboards rely on controlled delamination of internal fiber plies to form functional hinges without breaking surface liners. Folding boxboard (FBB) pairs virgin bleached chemical pulp liners with a high-yield mechanical pulp middle layer. Solid bleached sulfate (SBS) comprises pure chemical pulp plies throughout.

In both structures, the crease creates high interlaminar shear stresses adjacent to the rule edges. These shear fields initiate localized delamination cracks that propagate parallel to the board surfaces. Successful creasing transforms a single monolithic plate of high flexural rigidity into several unbonded thin laminates that slide across one another, reducing bending stiffness by over 70 percent along the crease axis.

Cohesive zone modeling (CZM) serves as the primary continuum mechanics tool for tracking this progressive interface separation. The cohesive formulation dictates traction-separation laws governing interlaminar interfaces. Bilinear or exponential traction laws prescribe that interfacial traction increases linearly with opening or sliding displacement until reaching critical interface strengths, designated as interlaminar shear strength and out-of-plane tensile strength.

Beyond these peak values, a softening branch dictates damage evolution until the fracture energy is exhausted, creating a traction-free crack surface.

Standard ASTM D5528 double cantilever beam evaluations require dynamic correction factors when sizing adhesive open times on high-speed side-seam folder-gluers.

Converting machinery running at sustained speed generates localized frictional heating inside these sheared interfaces. The temperature of the middle plies can elevate by 4 to 12 degrees Celsius during successive scoring and pre-breaking operations. Elevated temperatures accelerate viscoelastic relaxation within secondary fiber walls, shifting local shear behavior.

When modeling cohesive zones under rapid converting, the constitutive formulation couples temperature and moisture fields directly to interface toughness parameters. Disregarding this coupling produces structural models that predict catastrophic liner fracture where production reality displays clean, compliant internal delamination.

Friction

Surface Interaction under Dynamic Belt and Guide Contact

Converting lines transport blanks using vacuum belts, rubber nip rollers, and steel folding swords. Surface friction dictates blank tracking, squareness, and registration accuracy into die-cutting and folding modules. Friction between paperboard surfaces and mechanical components cannot be treated as a constant Coulomb coefficient.

The dynamic friction coefficient changes as a function of sliding velocity, contact pressure, and surface temperature. The presence of surface finishes, including dispersion barriers, mineral pigment coatings, varnishes, and film laminations, introduces a distinct viscoelastic friction component.

Polymer barrier films and overprint varnishes exhibit adhesion-deformation friction regimes. When a polyurethane feed belt grips a carton blank, real contact area depends on the time-dependent deformation of microscopic coating asperities. Under higher line speeds, asperity contact times drop into the sub-millisecond range, reducing contact area and diminishing static grip.

Conversely, sliding friction over polished chrome folding rails generates localized shear in the coating layer. At typical folding velocities between 3 and 10 meters per second, the shear rate inside the top varnish film reaches extreme magnitudes, inducing thermal softening and shear-thinning behavior.

• Interfacial sliding velocity alters the shear strain rate within polymer topcoats, shifting friction coefficients from 0.28 under quasi-static testing to over 0.45 at high-speed processing.
• Contact normal pressure compresses surface asperities, transitioning the contact profile from elastic asperity contact to bulk plastic deformation of the mineral coating layer.
• Substrate surface roughness measured under ISO 8791-4 controls air entrapment pockets that lubricate the interface at linear board speeds above 300 meters per minute.
• Ambient relative humidity changes moisture content in the outer cellulose layers, increasing capillary adhesion forces and driving register drift across long converting runs.

Miscalculating friction coefficients at converting line speeds leads to misregistration inside the transfer sections of side-seam gluers. When blanks slip along the folding swords, pre-breaking fails to reach specified angles, commonly 180 degrees for crease one and 135 degrees for crease three. Asymmetrical folding generates out-of-square cartons that jam modern high-speed packing lines operating at 600 cartons per minute.

The production loss compounds when hot-melt or cold-set adhesives land off-target, fouling carrier belts and requiring line washdowns that halt output.

Rebound

Viscoelastic Recovery and Springback Dynamics

Following high-speed mechanical folding, carton panels exhibit time-dependent elastic recovery. Springback compromises the dimensional accuracy of folded packaging blanks during packing line erection. When a folding sword bends a scored panel to 90 or 180 degrees, the tensile strain on the outer surface and the compressive strain on the inner surface set up an internal resisting moment.

In a perfectly elastic medium, removing the tool returns the panel to its flat configuration. In an elastoplastic medium, the panel recovers by an amount proportional to its yield stress. In a viscoelastic paperboard, the recovery angle increases over seconds, hours, and days following the conversion event.

The time-dependent opening moment governed by stress relaxation in the folded hinge follows a decaying exponential trajectory. Directly after exiting the compression section of a folder-gluer, internal stresses remain concentrated within compressed fiber knuckles inside the crease. As the folded carton cures in the collection hopper under stack pressure, these peak stresses relax.

When blanks are extracted from packed shipping cases days later on the customer filling line, residual bending moments attempt to force the panels open against erecting fingers, increasing load on automated carton feeders.

To eliminate carton jamming, packaging specifications set strict boundaries on the ratio between board bending resistance and crease opening moment. The crease stiffness ratio, defined as crease opening resistance divided by uncut board stiffness, must land between 0.30 and 0.50 according to typical converting criteria. If the ratio climbs above 0.55, automated packaging line erecters bow the main carton panels instead of breaking the creases cleanly, causing machinery shutdowns.

When the ratio falls below 0.25, the crease lacks structural integrity, leading to panel bulging and package collapse during pallet stacking under compressive top loads.

The mill-side excuse for panel springback invariably points to environmental swings on the packaging floor rather than inconsistent cross-direction stiffness in the base board.

Tooling

Rule and Matrix Parameterization in High Speed Runs

Tooling geometry fixes the mechanical deformation field imparted to the board substrate. A creasing configuration comprises a hardened steel rule with an edge radius, paired with a counter matrix groove defined by its width and depth. For a given board thickness and furnish structure, classical empirical rules dictate channel width through linear equations adding rule thickness to a fixed multiple of board caliper.

These proportional rules fail when machine speeds exceed 6000 sheets per hour because dynamic strain hardening alters how the fiber matrix draws into the groove.

Rule penetration depth controls the level of out-of-plane shear stress. Under dynamic impact, insufficient penetration prevents shear separation of internal plies, leaving the board behaving as a monolithic sheet with an elevated bending moment. Excessive penetration crushes the bottom chemical pulp liner, creating tensile failure along the exterior matrix shoulders.

When high-speed converting lines experience micro-vibrations across the platen bed, actual penetration depths fluctuate within a band of plus or minus 25 microns across the sheet area.

1. Rule radius selection establishes the contact pressure distribution across the top coating, requiring a 2-point rule (0.71 mm width) with a precision-ground 0.35 mm radius to minimize local stress concentrations that initiate surface cracking on pigment coatings.
2. Channel depth calibration fixes the vertical shear displacement imparted to middle plies, set to exactly match compressed board caliper minus the density-corrected volume loss under high-speed impact.
3. Channel width optimization sets the unsupported shear span during indentation, where narrow matrices induce catastrophic tensile liner rupture and overly wide matrices produce diffuse, ill-defined folding lines with uncontrolled panel wander.
4. Counter matrix composition dictates viscoelastic damping of the tooling base, where steel counter plates maintain rigid dimensional boundaries while phenolic or prespan counters exhibit dynamic compliance that shifts effective penetration depths during extended production shifts.

Rotary die-cutting operations require compensating for tangential speed differentials between the upper rule cylinder and the lower anvil matrix cylinder. Because the board moves at linear web velocity while the rule tip traverses a circular arc, the entry and exit phases apply non-symmetric horizontal forces to the creasing channel. This dynamic scuffing distorts the internal shear field, creating an asymmetric delamination pattern that biases panel folding toward one direction.

Constitutive models must incorporate tool-substrate kinematic contact formulations that resolve these tangential shear components under microsecond dwell times.

Creasing matrix tolerances exceeding 15 microns across the working width violate standard quality criteria outlined in ECMA performance guidelines for automated folding cartons.

Adjusting make-ready tape beneath a steel counter plate alters local penetration dynamics by fractions of a millimeter. When operators add shims during machine setup, they modify both static depth and the acoustic impedance of the tooling stack. The dynamic shock wave generated as the upper platen strikes the matrix reflects differently through layered paper tape versus solid milled steel.

High-speed recording reveals that platen bounce increases local rule dwell variance by 12 percent, directly degrading crease uniformity from the gripper edge to the tail of the carton layout.

A supplier contract specifying commercial board convertibility remains unworkable unless it defines minimum dynamic delamination energy alongside static ISO 2493 bending resistance limits.