Stress relaxation in network materials: the contribution of the network
Stress relaxation in network materials with permanent crosslinks is due to the transport of fluid within the network (poroelasticity), the viscoelasticity of the matrix and the viscoelasticity of the network. While relaxation associated with the matrix was studied extensively, the contribution of the network remains unexplored. In this work we consider two and three-dimensional stochastic fiber networks with viscoelastic fibers and explore the dependence of stress relaxation on network structure. We observe that relaxation has two regimes – an initial exponential regime, followed by a stretched exponential regime – similar to the situation in other disordered materials. The stretch exponent is a function of density, fiber diameter and the network structure, and has a minimum at the transition between the affine and non-affine regimes of network behavior. The relaxation time constant of the first, exponential regime is similar to the relaxation time constant of individual fibers and is independent of network density and fiber diameter. The relaxation time constant of the second, stretched exponential regime is a weak function of network parameters. The stretched exponential emerges from the heterogeneity of relaxation dynamics on scales comparable with the mesh size, with higher heterogeneity leading to smaller stretch exponents. In composite networks of fibers whose relaxation time constant is selected from a distribution with set mean, the stretch exponent decreases with increasing the coefficient of variation of the fiber time constant distribution. As opposed to thermal glass formers and colloids, in these athermal systems the dynamic heterogeneity is introduced by the network structure and does not evolve during relaxation. While in thermal systems the control parameter is the temperature, in this athermal case the control parameter is a non-dimensional structural parameter which describes the degree of non-affinity of the network.
more »
« less
- PAR ID:
- 10326131
- Publisher / Repository:
- Royal Society
- Date Published:
- Journal Name:
- Soft Matter
- Volume:
- 18
- Issue:
- 2
- ISSN:
- 1744-683X
- Page Range / eLocation ID:
- 446 to 454
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Kristin S. Miller. (Ed.)Abstract Cardiomyocytes are viscoelastic and key determinants of right ventricle (RV) mechanics. Intracellularly, microtubules are found to impact the viscoelasticity of isolated cardiomyocytes or trabeculae; whether they contribute to the tissue-level viscoelasticity is unknown. Our goal was to reveal the role of the microtubule network in the passive anisotropic viscoelasticity of the healthy RV. Equibiaxial stress relaxation tests were conducted in healthy RV free wall (RVFW) under early (6%) and end (15%) diastolic strain levels, and at sub- and physiological stretch rates. The viscoelasticity was assessed at baseline and after the removal of microtubule network. Furthermore, a quasi-linear viscoelastic (QLV) model was applied to delineate the contribution of microtubules to the relaxation behavior of RVFW. After removing the microtubule network, RVFW elasticity and viscosity were reduced at the early diastolic strain level and in both directions. The reduction in elasticity was stronger in the longitudinal direction, whereas the degree of changes in viscosity were equivalent between directions. There was insignificant change in RVFW viscoelasticity at late diastolic strain level. Finally, the modeling showed that the tissue's relaxation strength was reduced by the removal of the microtubule network, but the change was present only at a later time scale. These new findings suggest a critical role of cytoskeleton filaments in RVFW passive mechanics in physiological conditions.more » « less
-
Abstract Entropy dynamics is a Bayesian inference methodology that can be used to quantify time-dependent posterior probability densities that guide the development of complex material models using information theory. Here, we expand its application to non-Gaussian processes to evaluate how fractal structure can influence fractional hyperelasticity and viscoelasticity in elastomers. We investigate how kinematic constraints on fractal polymer network deformation influences the form of hyperelastic constitutive behavior and viscoelasticity in soft materials such as dielectric elastomers, which have applications in the development of adaptive structures. The modeling framework is validated on two dielectric elastomers, VHB 4910 and 4949, over a broad range of stretch rates. It is shown that local fractal time derivatives are equally effective at predicting viscoelasticity in these materials in comparison to nonlocal fractional time derivatives under constant stretch rates. We describe the origin of this accuracy that has implications for simulating large-scale problems such as finite element analysis given the differences in computational efficiency of nonlocal fractional derivatives versus local fractal derivatives.more » « less
-
We present a study of the mechanical behavior of planar fibrous mats stabilized by inter-fiber adhesion. Fibers of various degrees of tortuosity and of infinite and finite length are considered in separate models. Fibers are randomly distributed, are not cross-linked, and interact through adhesion and friction. The variation of structural parameters such as the mat thickness and the mean segment length between contacts along given fibers with the strength of adhesion is determined. These systems are largely dissipative in that most of the work performed during deformation is dissipated frictionally and only a small fraction is stored as strain energy. The response of the mats to tensile loading has three regimes: a short elastic regime in which no sliding at contacts is observed, a well-defined sliding regime characterized by strain hardening, and a rapid stiffening regime at larger strains. The third regime is due to the formation of stress paths after the fiber tortuosity is pulled out and is absent in mats of finite length fibers. Networks of finite length fibers lose stability during the second regime of deformation. The scaling of the yield stress, which characterizes the transition between the first and the second regimes, and of the second regime's strain hardening modulus, with system parameters such as the strength of adhesion and friction and the degree of fiber tortuosity are determined. The strength of mats of finite length fibers is also determined as a function of network parameters. These results are expected to become useful in the design of electrospun mats and other planar fibrous non-cross-linked networks.more » « less
-
A balance between model complexity, accuracy, and computational cost is a central concern in numerical simulations. In particular, for stochastic fiber networks, the non-affine deformation of fibers, related non-linear geometric features due to large global deformation, and size effects can significantly affect the accuracy of the computer experiment outputs and increase the computational cost. In this work, we systematically investigate methodological aspects of fiber network simulations with a focus on the output accuracy and computational cost in models with cellular (Voronoi) and fibrous (Mikado) network architecture. We study both p and h-refinement of the discretizations in finite element solution procedure, with uniform and length-based adaptive h-refinement strategies. The analysis is conducted for linear elastic and viscoelastic constitutive behavior of the fibers, as well as for networks with initially straight and crimped fibers. With relative error as the determining criterion, we provide recommendations for mesh refinement, comment on the necessity of multiple realizations, and give an overview of associated computational cost that will serve as guidance toward minimizing the computational cost while maintaining a desired level of solution accuracy.more » « less
An official website of the United States government

