Size Effects in Random Fiber Networks Controlled by the Use of Generalized Boundary Conditions

Evaluation of the parallel coupling constitutive model for biomaterials using a fully coupled network-matrix model

In this article we discuss the effective properties of composites containing a crosslinked athermal fiber network embedded in a continuum elastic matrix, which are representative for a broad range of biological materials. The goal is to evaluate the accuracy of the widely used biomechanics parallel coupling model in which the tissue response is defined as the additive superposition of the network and matrix contributions, and the interaction of the two components is neglected. To this end, explicit, fully coupled models are used to evaluate the linear and non-linear response of the composite. It is observed that in the small strain, linear regime the parallel model leads to errors when the ratio of the individual stiffnesses of the two components is in the range 0.1-10, and the error increases as the matrix approaches the incompressible limit. The data presented can be used to correct the parallel model to improve the accuracy of the overall stiffness prediction. In the non-linear large deformation regime linear superposition does not apply. The data shows that the matrix reduces the stiffening rate of the network, and the response is softer than that predicted by the parallel model. The correction proposed for the linear regime mitigates to a large extent the error in the non-linear regime as well, provided the matrix Poisson ratio is not close to 0.5. The special case in which the matrix is rendered auxetic is also evaluated and it is seen that the auxeticity of the matrix may compensate the stiffening introduced by the network, leading to a composite with linear elastic response over a broad range of strains.

Micromechanics of fibrous scaffolds and their stiffness sensing by cells

Mechanical properties of the tissue engineering scaffolds are known to play a crucial role in cell response. Therefore, an understanding of the cell-scaffold interactions is of high importance. Here, we have utilized discrete fiber network model to quantitatively study the micromechanics of fibrous scaffolds with different fiber arrangements and cross-linking densities. We observe that localized forces on the scaffold result in its anisotropic deformation even for isotropic fiber arrangements. We also see an exponential decay of the displacement field with distance from the location of applied force. This nature of the decay allows us to estimate the characteristic length for force transmission in fibrous scaffolds. Furthermore, we also looked at the stiffness sensing of fibrous scaffolds by individual cells and its dependence on the cellular sensing mechanism. For this, we considered two conditions- stress-controlled, and strain-controlled application of forces by a cell. With fixed strain, we find that the stiffness sensed by a cell is proportional to the scaffold's 'macroscopic' elastic modulus. However, under fixed stress application by the cell, the stiffness sensed by the cell also depends on the cell's own stiffness. In fact, the stiffness values for the same scaffold sensed by the stiff and soft cells can differ from each other by an order of magnitude. The insights from this work will help in designing tissue engineering scaffolds for applications where mechanical stimuli are a critical factor.

Emergence of an apparent yield phenomenon in the mechanics of stochastic networks with inter-fiber cohesion

In this work we investigate the contribution of inter-fiber cohesion to defining the mechanical behavior of stochastic crosslinked fiber networks. Fibers are athermal and store energy primarily in their bending and axial deformation modes. Cohesion between fibers is defined by an interaction potential. These structures are in equilibrium with the inter-fiber cohesive forces before external load is applied and their mechanical behavior is probed in uniaxial tension. Two types of configurations are considered: a state with high initial free volume in which contacts between fibers are scarce, and a state with low free volume and large number of fiber contacts. While in the absence of cohesion the response is hyperelastic, we observe that a yield point-like phenomenon develops as the strength of cohesion increases in both network types considered; we refer to this as an 'unlocking phenomenon'. The small strain stiffness increases as cohesion becomes more pronounced. The stiffness and unlocking stress are expressed in terms of network parameters and cohesion strength through a product of two functions, one dependent on network parameters only, and the other is a function of the cohesion strength. While the small strain response is controlled by cohesion, the large strain behavior is shown to be largely controlled by the network. Therefore, varying the strength of cohesion has no effect on strain stiffening. These observations provide a physical basis for the unlocking observed in both athermal and thermal network materials and are expected to facilitate the design of soft materials with novel properties.

Stiffening mechanisms in stochastic athermal fiber networks

Stochastic athermal networks composed of fibers that deform axially and in bending strain stiffen much faster than thermal networks of axial elements, such as elastomers. Here we investigate the physical origin of stiffening in athermal network materials. To this end, we use models of stochastic networks subjected to uniaxial deformation and identify the emergence of two subnetworks, the stress path subnetwork (SPSN) and the bending support subnetwork (BSSN), which carry most of the axial and bending energies, respectively. The BSSN controls lateral contraction and modulates the organization of the SPSN during deformation. The SPSN is preferentially oriented in the loading direction, while the BSSN's preferential orientation is orthogonal to the SPSN. In nonaffine networks stiffening is exponential, while in close-to-affine networks it is quadratic. The difference is due to a much more modest lateral contraction in the approximately affine case and to a stiffer BSSN. Exponential stiffening emerges from the interplay of the axial and bending deformation modes at the scale of individual or small groups of fibers undergoing large deformations and being subjected to the constraint of rigid cross-links, and it is not necessarily a result of complex interactions involving many connected fibers. An apparent third regime of quadratic stiffening may be evidenced in nonaffinely deforming networks provided the nominal stress is observed. This occurs at large stretches, when the BSSN contribution of stiffening vanishes. However, this regime is not present if the Cauchy stress is used, in which case stiffening is exponential throughout the entire deformation. These results shed light on the physical nature of stiffening in a broad class of materials including connective tissue, the extracellular matrix, nonwovens, felt, and other athermal network materials.

Probing soft fibrous materials by indentation

Indentation is often used to measure the stiffness of soft materials whose main structural component is a network of filaments, such as the cellular cytoskeleton, connective tissue, gels, and the extracellular matrix. For elastic materials, the typical procedure requires fitting the experimental force-displacement curve with the Hertz model, which predicts that f=kδ1.5 and k is proportional to the reduced modulus of the indented material, E/(1-ν2). Here we show using explicit models of fiber networks that the Hertz model applies to indentation in network materials provided the indenter radius is larger than approximately 12lc, where lc is the mean segment length of the network. Using smaller indenters leads to a relation between force and indentation displacement of the form f=kδq, where q is observed to increase with decreasing indenter radius. Using the Hertz model to interpret results of indentations in network materials using small indenters leads to an inferred modulus smaller than the real modulus of the material. The origin of this departure from the classical Hertz model is investigated. A compacted, stiff network region develops under the indenter, effectively increasing the indenter size and modifying its shape. This modification is marginal when large indenters are used. However, when the indenter radius is small, the effect of the compacted layer is pronounced as it changes the indenter profile from spherical towards conical. This entails an increase of exponent q above the value of 1.5 corresponding to spherical indenters. STATEMENT OF SIGNIFICANCE: The article presents a study of indentation in network biomaterials and demonstrates a size effect which precludes the use of the Hertz model to infer the elastic constants of the material. The size effect occurs once the indenter radius is smaller than approximately 12 times the mean segment length of the network. This result provides guidelines for the selection of indentation conditions that guarantee the applicability of the Hertz model. At the same time, the finding may be used to infer the mean segment length of the network based on indentations with indenters of various sizes. Hence, the method can be used to evaluate this structural parameter which is not easily accessible in experiments.

Random Fiber Network Loaded by a Point Force

This article presents the displacement field produced by a point force acting on an athermal random fiber network (the Green function for the network). The problem is defined within the limits of linear elasticity, and the field is obtained numerically for nonaffine networks characterized by various parameter sets. The classical Green function solution applies at distances from the point force larger than a threshold which is independent of the network parameters in the range studied. At smaller distances, the nonlocal nature of fiber interactions modifies the solution.