skip to main content


Title: A comparative study of two constitutive models within an inverse approach to determine the spatial stiffness distribution in soft materials
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
Award ID(s):
1663435
NSF-PAR ID:
10057051
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
International journal of mechanical sciences
Volume:
140
ISSN:
0020-7403
Page Range / eLocation ID:
446-454
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The deformations of several slender structures at nano-scale are conceivably sensitive to their non-homogenous elasticity. Owing to their small scale, it is not feasible to discern their elasticity parameter fields accurately using observations from physical experiments. Molecular dynamics simulations can provide an alternative or additional source of data. However, the challenges still lie in developing computationally efficient and robust methods to solve inverse problems to infer the elasticity parameter field from the deformations. In this paper, we formulate an inverse problem governed by a linear elastic model in a Bayesian inference framework. To make the problem tractable, we use a Gaussian approximation of the posterior probability distribution that results from the Bayesian solution of the inverse problem of inferring Young’s modulus parameter fields from available data. The performance of the computational framework is demonstrated using two representative loading scenarios, one involving cantilever bending and the other involving stretching of a helical rod (an intrinsically curved structure). The results show that smoothly varying parameter fields can be reconstructed satisfactorily from noisy data. We also quantify the uncertainty in the inferred parameters and discuss the effect of the quality of the data on the reconstructions. 
    more » « less
  2. Abstract

    A hyperelasticity modelling approach is employed for capturing various and complex mechanical behaviours exhibited by macroscopically isotropic polydomain liquid crystal elastomers (LCEs). These include the highly non-linear behaviour of nematic-genesis polydomain LCEs, and the soft elasticity plateau in isotropic-genesis polydomain LCEs, under finite multimodal deformations (uniaxial and pure shear) using in-house synthesised acrylate-based LCE samples. Examples of application to capturing continuous softening (i.e., in the primary loading path), discontinuous softening (i.e., in the unloading path) and auxetic behaviours are also demonstrated on using extant datasets. It is shown that our comparatively simple model, which breaks away from the neo-classical theory of liquid crystal elastomers, captures the foregoing behaviours favourably, simply as states of hyperelasticity. Improved modelling results obtained by our approach compared with the existing models are also discussed. Given the success of the considered model in application to these datasets and deformations, the simplicity of its functional form (and thereby its implementation), and comparatively low(er) number of parameters, the presented isotropic hyperelastic strain energy function here is suggested for: (i) modelling the general mechanical behaviour of LCEs, (ii) the backbone in the neo-classical theory, and/or (iii) the basic hyperelastic model in other frameworks where the incorporation of the director, anisotropy, viscoelasticity, temperature, softening etc parameters may be required.

     
    more » « less
  3. Mascia, Corrado ; Terracina, Andrea ; Tesei, Alberto (Ed.)
    We study a model of dislocations in two-dimensional elastic media. In this model, the displacement satisfies the system of linear elasticity with mixed displacement-traction homogeneous boundary conditions in the complement of an open curve in a bounded planar domain, and has a specified jump, the slip, across the curve, while the traction is continuous there. The stiffness tensor is allowed to be anisotropic and inhomogeneous. We prove well-posedness of the direct problem in a variational setting, assuming the coefficients are Lipschitz continuous. Using unique continuation arguments, we then establish uniqueness in the inverse problem of determining the dislocation curve and the slip from a single measurement of the displacement on an open patch of the traction-free part of the boundary. Uniqueness holds when the elasticity operators admits a suitable decomposition and the curve satisfies additional geometric assumptions. This work complements the results in Arch. Ration. Mech. Anal., 236(1):71-111, (2020), and in Preprint arXiv:2004.00321, which concern three-dimensional isotropic elastic media. 
    more » « less
  4. 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
  5. Accurate characterization of the mechanical properties of the human brain at both microscopic and macroscopic length scales is a critical requirement for modeling of traumatic brain injury and brain folding. To date, most experimental studies that employ classical tension/compression/shear tests report the mechanical properties of the brain averaged over both the gray and white matter within the macroscopic regions of interest. As a result, there is a missing correlation between the independent mechanical properties of the microscopic constituent elements and the composite bulk macroscopic mechanical properties of the tissue. This microstructural computational study aims to inversely predict the hyperelastic mechanical properties of the axonal fibers and their surrounding extracellular matrix (ECM) from the bulk tissue's mechanical properties. We develop a representative volume element (RVE) model of the bulk tissue consisting of axonal fibers and ECM with the embedded element technique. A multiobjective optimization technique is implemented to calibrate the model and establish the independent mechanical properties of axonal fibers and ECM based on seven previously reported experimental mechanical tests for bulk white matter tissue from the corpus callosum. The result of the study shows that the discrepancy between the reported values for the elastic behavior of white matter in literature stems from the anisotropy of the tissue at the microscale. The shear modulus of the axonal fiber is seven times larger than the ECM, with axonal fibers that also show greater nonlinearity, contrary to the common assumption that both components exhibit identical nonlinear characteristics. Statement of significance The reported mechanical properties of white matter microstructure used in traumatic brain injury or brain mechanics studies vary widely, in some cases by up to two orders of magnitude. Currently, the material parameters of the white matter microstructure are identified by a single loading mode or ultimately two modes of the bulk tissue. The presented material models only define the response of the bulk and homogenized white matter at a macroscopic scale and cannot explicitly capture the connection between the material properties of microstructure and bulk structure. To fill this knowledge gap, our study characterizes the hyperelastic material properties of axonal fibers and ECM using microscale computational modeling and multiobjective optimization. The hyperelastic material properties for axonal fibers and ECM presented in this study are more accurate than previously proposed because they have been optimized using seven or six loading modes of the bulk tissue, which were previously limited to only two of the seven possible loading modes. As such, the predicted values with high accuracy could be used in various computational modeling studies. The systematic characterization of the material properties of the human brain tissue at both macro- and microscales will lead to more accurate computational predictions, which will enable a better understanding of injury criteria, and has a positive impact on the improved development of smart protection systems, and more accurate prediction of brain development and disease progression. 
    more » « less