skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: 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
Award ID(s):
1824708 1561789
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Biomechanics and Modeling in Mechanobiology
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Intracranial aneurysms (ICAs) are focal dilatations that imply a weakening of the brain artery. Incidental rupture of an ICA is increasingly responsible for significant mortality and morbidity in the American’s aging population. Previous studies have quantified the pressure-volume characteristics, uniaxial mechanical properties, and morphological features of human aneurysms. In this pilot study,for the first time, we comprehensively quantified themechanical,collagen fiber microstructural, andmorphologicalproperties of one resected human posterior inferior cerebellar artery aneurysm. The tissue from the dome of a right posterior inferior cerebral aneurysm was first mechanically characterized using biaxial tension and stress relaxation tests. Then, the load-dependent collagen fiber architecture of the aneurysm tissue was quantified using an in-house polarized spatial frequency domain imaging system. Finally, optical coherence tomography and histological procedures were used to quantify the tissue’s microstructural morphology. Mechanically, the tissue was shown to exhibit hysteresis, a nonlinear stress-strain response, and material anisotropy. Moreover, the unloaded collagen fiber architecture of the tissue was predominantly aligned with the testingY-direction and rotated towards theX-direction under increasing equibiaxial loading. Furthermore, our histological analysis showed a considerable damage to the morphological integrity of the tissue, including lack of elastin, intimal thickening, and calcium deposition. This new unified characterization framework can be extended to better understand the mechanics-microstructure interrelationship of aneurysm tissues at different time points of the formation or growth. Such specimen-specific information is anticipated to provide valuable insight that may improve our current understanding of aneurysm growth and rupture potential.

    more » « less
  2. 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
  3. 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
  4. A comparative study is presented to solve the inverse problem in elasticity for the shear modulus (stiffness) distribution utilizing two constitutive equations: (1) linear elasticity assuming small strain theory, and (2) finite elasticity with a hyperelastic neo-Hookean material model. Assuming that a material undergoes large deformations and material nonlinearity is assumed negligible, the inverse solution using (2) is anticipated to yield better results than (1). Given the fact that solving a linear elastic model is significantly faster than a nonlinear model and more robust numerically, we posed the following question: How accurately could we map the shear modulus distribution with a linear elastic model using small strain theory for a specimen undergoing large deformations? To this end, experimental displacement data of a silicone composite sample containing two stiff inclusions of different sizes under uniaxial displacement controlled extension were acquired using a digital image correlation system. The silicone based composite was modeled both as a linear elastic solid under infinitesimal strains and as a neo-Hookean hyperelastic solid that takes into account geometrically nonlinear finite deformations. We observed that the mapped shear modulus contrast, determined by solving an inverse problem, between inclusion and background was higher for the linear elastic model as compared to that of the hyperelastic one. A similar trend was observed for simulated experiments, where synthetically computed displacement data were produced and the inverse problem solved using both, the linear elastic model and the neo-Hookean material model. In addition, it was observed that the inverse problem solution was inclusion size-sensitive. Consequently, an 1-D model was introduced to broaden our understanding of this issue. This 1-D analysis revealed that by using a linear elastic approach, the overestimation of the shear modulus contrast between inclusion and background increases with the increase of external loads and target shear modulus contrast. Finally, this investigation provides valuable information on the validity of the assumption for utilizing linear elasticity in solving inverse problems for the spatial distribution of shear modulus associated with soft solids undergoing large deformations. Thus, this work could be of importance to characterize mechanical property variations of polymer based materials such as rubbers or in elasticity imaging of tissues for pathology. 
    more » « less
  5. null (Ed.)
    Synopsis The adhesive toe pads of tree frogs have inspired the design of various so-called ‘smooth’ synthetic adhesives for wet environments. However, these adhesives do not reach the attachment performance of their biological models in terms of contact formation, maintenance of attachment, and detachment. In tree frogs, attachment is facilitated by an interconnected ensemble of superficial and internal morphological components, which together form a functional unit. To help bridging the gap between biological and bioinspired adhesives, in this review, we (1) provide an overview of the functional components of tree frog toe pads, (2) investigate which of these components (and attachment mechanisms implemented therein) have already been transferred into synthetic adhesives, and (3) highlight functional analogies between existing synthetic adhesives and tree frogs regarding the fundamental mechanisms of attachment. We found that most existing tree-frog-inspired adhesives mimic the micropatterned surface of the ventral epidermis of frog pads. Geometrical and material properties differ between these synthetic adhesives and their biological model, which indicates similarity in appearance rather than function. Important internal functional components such as fiber-reinforcement and muscle fibers for attachment control have not been considered in the design of tree-frog-inspired adhesives. Experimental work on tree-frog-inspired adhesives suggests that the micropatterning of adhesives with low-aspect-ratio pillars enables crack arresting and the drainage of interstitial liquids, which both facilitate the generation of van der Waals forces. Our analysis of experimental work on tree-frog-inspired adhesives indicates that interstitial liquids such as the mucus secreted by tree frogs play a role in detachment. Based on these findings, we provide suggestions for the future design of biomimetic adhesives. Specifically, we propose to implement internal fiber-reinforcements inspired by the fibrous structures in frog pads to create mechanically reinforced soft adhesives for high-load applications. Contractile components may stimulate the design of actuated synthetic adhesives with fine-tunable control of attachment strength. An integrative approach is needed for the design of tree-frog-inspired adhesives that are functionally analogous with their biological paradigm. 
    more » « less