Hydrogel mechanics below swelling equilibrium

Technical Summary: Hydrogel Mechanics Below Swelling Equilibrium

Problem Statement
Hydrogels are soft, water-rich polymer networks widely used in medicine, bioengineering, and soft robotics. While their mechanics are well-characterized in near-fully hydrated states, practical applications often involve dehydration, where environmental conditions cause water loss. Existing models typically focus on the fully swollen regime, where shape and volumetric deformations are tightly coupled: local stretching induces local changes in water content, and vice versa. However, it remains unclear whether drying merely causes hydrogels to shrink and stiffen or if it fundamentally alters their structure-property relationships. Specifically, there is a need to understand if the complex, fully coupled poroelasticity required for swollen gels can be simplified for drying gels without sacrificing accuracy.

Methodology
The authors combine experimental characterization with theoretical modeling to investigate the mechanics of drying hydrogels.

• Experimental Approach: The study utilizes poly(ethylene glycol) diacrylate (PEGDA) and polyacrylamide (PAM) hydrogels. The authors measure the drained compressive modulus ($M_{1D}$) and shear modulus ($G$) as the polymer volume fraction ($\phi$) increases due to drying. Compressive stiffness is measured via gel freezing osmometry (using cryosuction to drive water out), while shear stiffness is determined through indentation tests on partially dehydrated samples.
• Theoretical Framework: The authors employ standard poroelasticity theory, decomposing the total Cauchy stress ($\sigma$) into an elastic stress ($\sigma_{el}$) resisting network deformation and a mixing pressure ($\Pi_{mix}$) reflecting polymer-water affinity, balanced against interstitial water pressure ($p$).
  • Strongly Coupled Formulation: A reference model where volume change and shape deformation are solved simultaneously, treating $\sigma_{el}$ and $\Pi_{mix}$ as interdependent functions of displacement and polymer content.
  • Weakly Coupled Formulation: A proposed simplified model based on the observation that in drying gels, the mixing pressure dominates the volumetric response. This model decouples the problem into a hierarchy: first, the local polymer volume fraction is determined solely by the ambient pore pressure ($\phi \approx \Pi_{mix}^{-1}(p)$); second, the shape deformation is solved as an elastic problem subject to this prescribed volumetric constraint, using a Lagrange multiplier to enforce the volume change.

Key Contributions and Results

1. Separation of Stiffness Scales: Experiments reveal a stark disparity in mechanical response as hydrogels dry. The compressive modulus ($M_{1D}$) increases by two orders of magnitude (from $\sim1$ MPa to $\sim100$ MPa) as the polymer volume fraction rises, driven by the increasing difficulty of stripping water molecules (polymer-water affinity). In contrast, the shear modulus ($G$) remains relatively constant ($\sim0.1$ MPa), indicating that the polymer network elasticity changes little with dehydration.
2. Decoupling of Deformation Modes: This stiffness separation implies that in drying hydrogels, the resistance to volume change is dominated by polymer-water affinity, while the resistance to shape change is governed by the much weaker network elasticity. Consequently, the mixing pressure ($\Pi_{mix}$) effectively equals the pore pressure ($p$), regardless of mechanical loading or shape distortion.
3. Validation of the Simplified Model: The authors demonstrate that the "weakly coupled" formulation accurately predicts stresses and deformations under complex conditions, including anisotropic shrinkage in annular gels and extreme deformations in chiral spiral structures. The model's predictions match the "strongly coupled" reference solution across a wide range of uniaxial stretches and suction levels. Deviations are confined to a narrow regime of severe uniaxial compression combined with weak pore pressure.
4. Computational Efficiency: Beyond accuracy, the weakly coupled formulation offers significant computational advantages. By reducing the problem to a sequential hierarchy, it avoids the iterative root-finding required by the strongly coupled model. Numerical benchmarks show the simplified model reduces cumulative CPU time by approximately 25% in steady-state simulations and allows for much larger time steps in transient simulations due to improved numerical conditioning and robustness.

Significance and Claims
The paper claims to reveal a mechanical hierarchy in drying hydrogels: it is orders of magnitude harder to change their volume by squeezing out water than to deform their shape. This insight allows for the decoupling of the physics governing hydration from the physics governing shape deformation.

The authors assert that this simplified viewpoint provides a robust tool for modeling hydrogels in drying or swelling scenarios. Potential applications identified within the text include:

• Predicting tensile stresses and cracking in drying gels and epithelial tissues.
• Understanding swelling instabilities and volume changes in hydrogel sorbants.
• Designing hydrogel actuators with pre-programmed shape changes.
• Developing load-invariant sensing techniques, where hydrogel volume changes serve as a direct measure of water chemical potential (e.g., humidity, temperature, pore pressure) independent of mechanical loading.

The authors note that while this hierarchy is expected to be broadly applicable due to the inherent affinity of hydrogels for water, exceptions may exist for gels with strongly nonlinear elastic networks, those undergoing glass transitions or crystallization, or stimuli-responsive systems. They conclude that future research should verify the universality of this hierarchy across different hydrogel chemistries.