An official website of the United States government
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
Characterization of the stiffness distribution in two and three dimensions using boundary deformations: a preliminary study
We present for the first time the feasibility to recover the stiffness (here shear modulus) distribution of a three-dimensional heterogeneous sample using measured surface displacements and inverse algorithms without making any assumptions about local homogeneities and the stiffness distribution. We simulate experiments to create measured displacements and augment them with noise, significantly higher than anticipated measurement noise. We also test two-dimensional problems in plane strain with multiple stiff inclusions. Our inverse strategy recovers the shear modulus values in the inclusions and background well, and reveals the shape of the inclusion clearly.
more »
« less
- Award ID(s):
- 1663435
- PAR ID:
- 10096287
- Date Published:
- Journal Name:
- MRS Communications
- Volume:
- 8
- Issue:
- 03
- ISSN:
- 2159-6859
- Page Range / eLocation ID:
- 893 to 902
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Intervertebral disc (IVD) degeneration is a significant health issue that can lead to severe complications. Recent research has highlighted the close relationship between disc degeneration and the biomechanical properties of the IVD. This study introduces an innovative approach—magnetic resonance imaging (MRI) elastography of the human IVD—using an explicit inverse solver to identify the non-homogeneous shear modulus map of the IVD. The key advantage of this explicit solver is its streamlined optimization process, focusing only on the shear moduli of the nucleus pulposus (NP), annulus fibrosus (AF), and their interface. This approach reduces the optimization variables, making it more efficient than traditional pixel-based approaches. To validate this method, we conducted a plane strain numerical example, observing a consistent underestimation of the AF/NP shear modulus ratio by a scaling factor of approximately 1.5. Further investigations included comprehensive sensitivity analyses to various noise levels, revealing that the proposed method accurately characterizes shear modulus distribution in the AF and NP regions, with a maximum relative error of the AF/NP shear modulus ratio remaining below 8%. In addition, applying this approach to real human IVDs underin vitrocompression or bending, demonstrated its effectiveness, yielding an AF/NP shear modulus ratio within a reasonable range of 6–15. In summary, the proposed method offers a promising direction for MRI elastography of the human IVD.more » « less
-
Abstract Auxetic materials have a negative Poisson’s ratio and are of significant interest in applications that include impact mitigation, membrane separations and biomedical engineering. While there are numerous examples of structured materials that exhibit auxetic behavior, the examples of engineered auxetic structures is largely limited to periodic lattice structures that are limited to directional or anisotropic auxetic response. Structures that exhibit a three-dimensionally isotropic auxetic response have been, unfortunately, slow to evolve. Here we introduce an inverse design algorithm based on global node optimization to design three-dimensional auxetic metamaterial structures from disordered networks. After specifying the target Poisson’s ratio for a structure, an inverse design algorithm is used to adjust the positions of all nodes in a disordered network structure until the desired mechanical response is achieved. The proposed algorithm allows independent control of shear and bulk moduli, while preserving the density and connectivity of the networks. When the angle bending stiffness in the network is kept low, it is possible to realize optimized structures with a Poisson’s ratios as low as −0.6. During the optimization, the bulk modulus of these networks decreases by almost two orders of magnitude, but the shear modulus remains largely unaltered. The materials designed in this manner are fabricated by dual-material 3D-printing, and are found to exhibit the mechanical responses that were originally encoded in the computational design engine. The approach proposed here provides a materials-by-design platform that could be extended for engineering of optical, acoustic, and electrical properties, beyond the design of auxetic metamaterials.more » « less
-
Bacteria experience substantial physical forces in their natural environment including forces caused by osmotic pressure, growth in constrained spaces, and fluid shear. The cell envelope is the primary load-carrying structure of bacteria, but the mechanical properties of the cell envelope are poorly understood; reports of Young’s modulus of the cell envelope of E. coli are widely range from 2 MPa to 18 MPa. We have developed a microfluidic system to apply mechanical loads to hundreds of bacteria at once and demonstrated the utility of the approach for evaluating whole-cell stiffness. Here we extend this technique to determine Young’s modulus of the cell envelope of E. coli and of the pathogens V. cholerae and S. aureus. An optimization-based inverse finite element analysis was used to determine the cell envelope Young’s modulus from observed deformations. The Young’s modulus of the cell envelope was 2.06±0.04 MPa for E. coli, 0.84±0.02 MPa for E. coli treated with a chemical known to reduce cell stiffness, 0.12±0.03 MPa for V. cholerae, and 1.52±0.06 MPa for S. aureus (mean ± SD). The microfluidic approach allows examining hundreds of cells at once and is readily applied to Gram-negative and Gram-positive organisms as well as rod-shaped and cocci cells, allowing further examination of the structural causes of differences in cell envelope Young's modulus among bacteria species and strains.more » « less
-
The extracellular matrix provides macroscale structural support to tissues as well as microscale mechanical cues, like stiffness, to the resident cells. As those cues modulate gene expression, proliferation, differentiation, and motility, quantifying the stiffness that cells sense is crucial to understanding cell behavior. Whereas the macroscopic modulus of a collagen network can be measured in uniform extension or shear, quantifying the local stiffness sensed by a cell remains a challenge due to the inhomogeneous and nonlinear nature of the fiber network at the scale of the cell. To address this challenge, we designed an experimental method to measure the modulus of a network of collagen fibers at this scale. We used spherical particles of an active hydrogel (poly N-isopropylacrylamide) that contract when heated, thereby applying local forces to the collagen matrix and mimicking the contractile forces of a cell. After measuring the particles’ bulk modulus and contraction in networks of collagen fibers, we applied a nonlinear model for fibrous materials to compute the modulus of the local region surrounding each particle. We found the modulus at this length scale to be highly heterogeneous, with modulus varying by a factor of 3. In addition, at different values of applied strain, we observed both strain stiffening and strain softening, indicating nonlinearity of the collagen network. Thus, this experimental method quantifies local mechanical properties in a fibrous network at the scale of a cell, while also accounting for inherent nonlinearity.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1557/mrc.2018.98
