skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • A comparative study of two constitutive models within an inverse approach to determine the spatial stiffness distribution in soft materials
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. ; ; ; ( , Rendiconti di matematica e delle sue applicazioni)
    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
  2. ; ; ( , Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences)
    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
  3. ; ; ( , Rubber Chemistry and Technology)
    ABSTRACT

    Corners concentrate elastic fields and often initiate fracture. For small deformations, it is well established that the elastic field around a corner is power-law singular. For large deformations, we show here that the elastic field around a corner is concentrated but bounded. We conduct computation and an experiment on the lap shear of a highly stretchable material. A rectangular sample was sandwiched between two rigid substrates, and the edges of the stretchable material met the substrates at 90° corners. The substrates were pulled to shear the sample. We computed the large-deformation elastic field by assuming several models of elasticity. The theory of elasticity has no length scale, and lap shear is characterized by a single length, the thickness of the sample. Consequently, the field in the sample was independent of any length once the spatial coordinates were normalized by the thickness. We then lap sheared samples of a polyacrylamide hydrogel of various thicknesses. For all samples, fracture initiated from corners, at a load independent of thickness. These experimental findings agree with the computational prediction that large-deformation elastic fields at corners are concentrated but bounded.

     
    more » « less
  4. ; ; ( , Journal of Elasticity)
    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
  5. ; ; ( , Meccanica)
    Abstract

    The most widely-used representation of the compressible, isotropic, neo-Hookean hyperelastic model is considered in this paper. The version under investigation is that which is implemented in the commercial finite element software ABAQUS, ANSYS and COMSOL. Transverse stretch solutions are obtained for the following homogeneous deformations: uniaxial loading, equibiaxial loading in plane stress, and uniaxial loading in plane strain. The ground-state Poisson’s ratio is used to parameterize the constitutive model, and stress solutions are computed numerically for the physically permitted range of its values. Despite its broad application to a number of engineering problems, the physical limitations of the model, particularly in the small to moderate stretch regimes, are not explored. In this work, we describe and analyze results and make some critical observations, underlining the model’s advantages and limitations. For example, a snap-back feature of the transverse stretch is identified in uniaxial compression, a physically undesirable behavior unless validated by experimental data. The domain of this non-unique solution is determined in terms of the ground-state Poisson’s ratio and the state of stretch and stress. The analyses we perform are essential to enable the understanding of the characteristics of the standard, compressible, isotropic, neo-Hookean model used in ABAQUS, ANSYS and COMSOL. In addition, our results provide a framework for the parameter-fitting procedure needed to characterize this standard, compressible, isotropic neo-Hookean model in terms of experimental data.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email