Evolution of Force Chains Explains the Onset of Strain Stiffening in Fiber Networks

Fiber networks are the primary structural components of many biological structures, including the cell cytoskeleton and the extracellular matrix. These materials exhibit global nonlinearities, such as stiffening in extension and shear, during which the fibers bend and align with the direction of applied loading. Precise details of deformations at the scale of the fibers during strain stiffening are still lacking, however, as prior work has studied fiber alignment primarily from a qualitative perspective, which leaves incomplete the understanding of how the local microstructural evolution leads to the global mechanical behavior. To fill this gap, we studied how axial forces are transmitted inside the fiber network along paths called force chains, which continuously evolve during the course of deformation. We performed numerical simulations on two-dimensional networks of random fibers under uniaxial extension and shear, modeling the fibers using beam elements in finite element software. To quantify the force chains, we identified all chains of connected fibers for which the axial force was larger than a preset threshold and computed the total length of all such chains. To study the evolution of force chains during loading, we computed the derivative of the total length of all force chains with respect to the applied engineering strain. Results showed that the highest rate of evolution of force chains coincided with the global critical strain for strain stiffening of the fiber network. Therefore, force chains are an important factor connecting understanding of the local kinematics and force transmission to the macroscale stiffness of the fiber network.