Dynamic networks containing multiple bond types within a continuous network grant engineers another design parameter – relative bond fraction – by which to tune storage and dissipation of mechanical energy. However, the mechanisms governing emergent properties are difficult to deduce experimentally. Therefore, we here employ a network model with prescribed fractions of dynamic and stable bonds to predict relaxation characteristics of hybrid networks. We find that during stress relaxation, predominantly dynamic networks can exhibit long-term moduli through conformationally inhibited relaxation of stable bonds due to exclusion interactions with neighboring chains. Meanwhile, predominantly stable networks exhibit minor relaxation through non-affine reconfiguration of dynamic bonds. Given this, we introduce a single fitting parameter, ξ , to Transient Network Theory via a coupled rule of mixture, that characterizes the extent of stable bond relaxation. Treating ξ as a fitting parameter, the coupled rule of mixture's predicted stress response not only agrees with the network model's, but also unveils likely micromechanical traits of gels hosting multiple bond dissociation timescales.
more »
« less
A network model of transient polymers: exploring the micromechanics of nonlinear viscoelasticity
Dynamic networks contain crosslinks that re-associate after disconnecting, imparting them with viscoelastic properties. While continuum approaches have been developed to analyze their mechanical response, these approaches can only describe their evolution in an average sense, omitting local, stochastic mechanisms that are critical to damage initiation or strain localization. To address these limitations, we introduce a discrete numerical model that mesoscopically coarse-grains the individual constituents of a dynamic network to predict its mechanical and topological evolution. Each constituent consists of a set of flexible chains that are permanently cross-linked at one end and contain reversible binding sites at their free ends. We incorporate nonlinear force–extension of individual chains via a Langevin model, slip-bond dissociation through Eyring's model, and spatiotemporally-dependent bond attachment based on scaling theory. Applying incompressible, uniaxial tension to representative volume elements at a range of constant strain rates and network connectivities, we then compare the mechanical response of these networks to that predicted by the transient network theory. Ultimately, we find that the idealized continuum approach remains suitable for networks with high chain concentrations when deformed at low strain rates, yet the mesoscale model proves necessary for the exploration of localized stochastic events, such as variability of the bond kinetics, or the nucleation of micro-cavities that likely conceive damage and fracture.
more »
« less
- Award ID(s):
- 1761918
- NSF-PAR ID:
- 10295769
- Date Published:
- Journal Name:
- Soft Matter
- ISSN:
- 1744-683X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Many biological materials contain fibrous protein networks as their main structural components. Understanding the mechanical properties of such networks is important for creating biomimicking materials for cell and tissue engineering, and for developing novel tools for detecting and diagnosing disease. In this work, we develop continuum models for isotropic, athermal fibrous networks by combining a single-fibre model that describes the axial response of individual fibres, with network models that assemble individual fibre properties into overall network behaviour. In particular, we consider four different network models, including the affine, three-chain, eight-chain, and micro-sphere models, which employ different assumptions about network structure and kinematics. We systematically investigate the ability of these models to describe the mechanical response of athermal collagen and fibrin networks by comparing model predictions with experimental data. We test how each model captures network behaviour under three different loading conditions: uniaxial tension, simple shear, and combined tension and shear. We find that the affine and three-chain models can accurately describe both the axial and shear behaviour, whereas the eight-chain and micro-sphere models fail to capture the shear response, leading to unphysical zero shear moduli at infinitesimal strains. Our study is the first to systematically investigate the applicability of popular network models for describing the macroscopic behaviour of athermal fibrous networks, offering insights for selecting efficient models that can be used for large-scale, finite-element simulations of athermal networks.more » « less
-
Toughness of soft materials such as elastomers and gels depends on their ability to dissipate energy and to reduce stress concentration at the crack tip. The primary energy dissipation mechanism is viscoelasticity. Most analyses and models of fracture are based on linear viscoelastic theory (LVT) where strains are assumed to be small and the relaxation mechanisms are independent of stress or strain history. A well-known paradox is that the size of the dissipative zone predicted by LVT is unrealistically small. Here we use a physically based nonlinear viscoelastic model to illustrate why the linear theory breaks down. Using this nonlinear model and analogs of crack problems, we give a plausible resolution to this paradox. In our model, viscoelasticity arises from the breaking and healing of physical cross-links in the polymer network. When the deformation is small, the kinetics of bond breaking and healing are independent of the strain/stress history and the model reduces to the standard linear theory. For large deformations, localized bond breaking damages the material near the crack tip, reducing stress concentration and dissipating energy at the same time. The damage zone size is a new length scale which depends on the strain required to accelerate bond breaking kinetics. These effects are illustrated by considering two cases with stress concentrations: the evolution of spherical damage in a viscoelastic body subjected to internal pressure, and a zero degree peel test.more » « less
-
more » « less
Abstract A continuum damage model was developed to describe the finite tensile deformation of tough double‐network (DN) hydrogels synthesized by polymerization of a water‐soluble monomer inside a highly crosslinked rigid polyelectrolyte network. Damage evolution in DN hydrogels was characterized by performing loading‐unloading tensile tests and oscillatory shear rheometry on DN hydrogels synthesized from 3‐sulfopropyl acrylate potassium salt (SAPS) and acrylamide (AAm). The model can explain all the mechanical features of finite tensile deformation of DN hydrogels, including idealized Mullins effect and permanent set observed after unloading, qualitatively and quantitatively. The constitutive equation can describe the finite elasto‐plastic tensile behavior of DN hydrogels without resorting to a yield function. It was showed that tensile mechanics of DN hydrogels in the model is controlled by two material parameters which are related to the elastic moduli of first and second networks. In effect, the ratio of these two parameters is a dimensionless number that controls the behavior of material. The model can capture the stable branch of material response during neck propagation where engineering stress becomes constant. Consistent with experimental data, by increasing the elastic modulus of the second network the finite tensile behavior of the DN hydrogel changes from necking to strain hardening.
-
Abstract 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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D1SM00753J
