An official website of the United States government
We propose and study a nonlinear elimination preconditioned inexact Newton method for the numerical simulation of diseased human arteries with a heterogeneous hyperelastic model. We assume the artery is made of layers of distinct tissues and also contains plaque. Traditional Newton methods often work well for smooth and homogeneous arteries but suffer from slow or no convergence due to the heterogeneousness of diseased soft tissues when the material is quasi-incompressible. The proposed nonlinear elimination method adaptively finds a small number of equations causing the nonlinear stagnation and then eliminates them from the global nonlinear system. By using the theory of affine invariance of Newton method, we provide insight into why the nonlinear elimination method can improve the convergence of Newton iterations. Our numerical results show that the combination of nonlinear elimination with an initial guess interpolated from a coarse level solution can lead to the uniform convergence of Newton method for this class of very difficult nonlinear problems.
more »
« less
Wrinkling instabilities for biologically relevant fiber-reinforced composite materials with a case study of Neo-Hookean/Ogden–Gasser–Holzapfel bilayer
Wrinkling is a ubiquitous surface phenomenon in many biological tissues and is believed to play an important role in arterial health. As arteries are highly nonlinear, anisotropic, multilayered composite systems, it is necessary to investigate wrinkling incorporating these material characteristics. Several studies have examined surface wrinkling mechanisms with nonlinear isotropic material relationships. Nevertheless, wrinkling associated with anisotropic constitutive models such as Ogden–Gasser–Holzapfel (OGH), which is suitable for soft biological tissues, and in particular arteries, still requires investigation. Here, the effects of OGH parameters such as fibers’ orientation, stiffness, and dispersion on the onset of wrinkling, wrinkle wavelength and amplitude are elucidated through analysis of a bilayer system composed of a thin, stiff neo-Hookean membrane and a soft OGH substrate subjected to compression. Critical contractile strain at which wrinkles occur is predicted using both finite element analysis and analytical linear perturbation approach. Results suggest that besides stiffness mismatch, anisotropic features associated with fiber stiffness and distribution might be used in natural layered systems to adjust wrinkling and subsequent folding behaviors. Further analysis of a bilayer system with fibers in the (x–y) plane subjected to compression in the x direction shows a complex dependence of wrinkling strain and wavelength on fiber angle, stiffness, and dispersion. This behavior is captured by an approximation utilizing the linearized anisotropic properties derived from OGH model. Such understanding of wrinkling in this artery wall-like system will help identify the role of wrinkling mechanisms in biological artery in addition to the design of its synthetic counterparts.
more »
« less
- PAR ID:
- 10191992
- Date Published:
- Journal Name:
- Biomechanics and Modeling in Mechanobiology
- ISSN:
- 1617-7959
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Soft biological tissues often function as highly deformable membranes in vivo and exhibit impressive mechanical behavior effectively characterized by planar biaxial testing. The Generalized Anisotropic Inverse Mechanics (GAIM) method links full-field deformations and boundary forces from mechanical testing to quantify material properties of soft, anisotropic, heterogeneous tissues. In this study, we introduced an orthotropic constraint to GAIM to improve the quality and physical significance of its mechanical characterizations. We evaluated the updated GAIM method using simulated and experimental biaxial testing datasets obtained from soft tissue analogs (PDMS and TissueMend) with well-defined mechanical properties. GAIM produced stiffnesses (first Kelvin moduli, K1) that agreed well with previously published Young's moduli of PDMS samples. It also matched the stiffness moduli determined via uniaxial testing for TissueMend, a collagen-rich patch intended for tendon repair. We then conducted the first biaxial testing of TissueMend and confirmed that the sample was mechanically anisotropic via a relative anisotropy metric produced by GAIM. Next, we demonstrated the benefits of full-field laser micrometry in distinguishing between spatial variations in thickness and stiffness. Finally, we conducted an analysis to verify that results were independent of partitioning scheme. The success of the newly implemented constraints on GAIM suggests notable potential for applying this tool to soft tissues, particularly following the onset of pathologies that induce mechanical and structural heterogeneities.more » « less
-
Over the past decades, the buckling instability of layered materials has been the subject of analytical, experimental, and numerical research. These systems have traditionally been considered with stress-free surfaces, and the influence of surface pressure is understudied. In this study, we developed a finite element model of a bilayer experiencing compression, and found that it behaves differently under surface pressure. We investigated the onset of buckling, the initial wavelength, and the post-buckling behavior of a bilayer system under two modes of compression (externally applied and internally generated by growth). Across a wide range of stiffness ratios, 1 < μf/μs < 100, we observed decreased stability in the presence of surface pressure, especially in the low-stiffness-contrast regime, μf/μs < 10. Our results suggest the importance of pressure boundary conditions for the stability analysis of bilayered systems, especially in soft and living matter physics, such as folding of the cerebral cortex under cerebrospinal fluid pressure, where pressure may affect morphogenesis and buckling patterns.more » « less
-
Abstract A simple equation modelling an inextensible elastic lining of an inner-lined tube subject to an imposed pressure difference is derived from a consideration of the idealised elastic properties of the lining and the pressure and soft-substrate forces. Two cases are considered in detail, one with prominent wrinkling and a second one in which wrinkling is absent and only buckling remains. Bifurcation diagrams are computed via numerical continuation for both cases. Wrinkling, buckling, folding, and mixed-mode solutions are found and organised according to system-response measures including tension, in-plane compression, maximum curvature and energy. Approximate wrinkle solutions are constructed using weakly nonlinear theory, in excellent agreement with numerics. Our approach explains how the wavelength of the wrinkles is selected as a function of the parameters in compressed wrinkling systems and shows how localised folds and mixed-mode states form in secondary bifurcations from wrinkled states. Our model aims to capture the wrinkling response of arterial endothelium to blood pressure changes but applies much more broadly.more » « less
-
Abstract Tension programmed shape memory polymer (SMP) fibers have been used as sutures for closing wide‐opened cracks per the close‐then‐heal strategy. However, the composite may be subjected to compression loading during service. These compression loads can reduce the amount of recoverable strain in these pre‐tensioned fibers, limiting their ability to close cracks. The purpose of this study is to investigate the effect of in‐service compression loading on the shape memory effect (SME) of composites consisting of SMP fiber and SMP matrix. To this end, pre‐stretched shape memory Polyethylene Terephthalate (PET) fibers were embedded into a shape memory vitrimer to obtain composite samples with different fiber volume fractions. The SME of both the PET fiber and the vitrimer was investigated. The effect of compression load on the SME of the composite was studied. It is found that, uniaxial compression on the composite along the fiber direction significantly reduced the shrinking ability of the embedded pre‐tensioned SMP fibers. Hence, this is a factor that needs to be considered when designing such types of self‐healing composites.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s10237-020-01345-0
