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
- NSF-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
-
-
more » « less
Abstract Cell migration and traction are essential to many biological phenomena, and one of their key features is sensitivity to substrate stiffness, which biophysical models, such as the motor‐clutch model and the cell migration simulator can predict and explain. However, these models have not accounted for the finite size of adhesions, the spatial distribution of forces within adhesions. Here, we derive an expression that relates varying adhesion radius (
R s κ subE Y , for two‐dimensional cell cultures. Experimentally, we found that a cone‐shaped force distribution ( s = 1.05 -
more » « less
Abstract Tissues and engineered biomaterials exhibit exquisite local variation in stiffness that defines their function. Conventional elastography quantifies stiffness in soft (e.g. brain, liver) tissue, but robust quantification in stiff (e.g. musculoskeletal) tissues is challenging due to dissipation of high frequency shear waves. We describe new development of
finite deformation elastography that utilizes magnetic resonance imaging of low frequency, physiological-level (large magnitude) displacements, coupled to an iterative topology optimization routine to investigate stiffness heterogeneity, including spatial gradients and inclusions. We reconstruct 2D and 3D stiffness distributions in bilayer agarose hydrogels and silicon materials that exhibit heterogeneous displacement/strain responses. We map stiffness in porcine and sheep articular cartilage deep within the bony articular joint spacein situ for the first time. Elevated cartilage stiffness localized to the superficial zone is further related to collagen fiber compaction and loss of water content during cyclic loading, as assessed by independentT 2measurements. We additionally describe technical challenges needed to achievein vivo elastography measurements. Our results introduce new functional imaging biomarkers, which can be assessed nondestructively, with clinical potential to diagnose and track progression of disease in early stages, including osteoarthritis or tissue degeneration. -
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
-
Abstract Characterizing the elastic properties of soft materials through bulge testing relies on accurate measurement of deformation, which is experimentally challenging. To avoid measuring deformation, we propose a hydrodynamic bulge test for characterizing the material properties of thick, pre-stressed elastic sheets via their fluid–structure interaction with a steady viscous fluid flow. Specifically, the hydrodynamic bulge test relies on a pressure drop measurement across a rectangular microchannel with a deformable top wall. We develop a mathematical model using first-order shear deformation theory of plates with stretching and the lubrication approximation for the Newtonian fluid flow. Specifically, a relationship is derived between the imposed flowrate and the total pressure drop. Then, this relationship is inverted numerically to yield estimates of the Young’s modulus (given the Poisson ratio) if the pressure drop is measured (given the steady flowrate). Direct numerical simulations of two-way-coupled fluid–structure interaction are carried out in ansys to determine the cross-sectional membrane deformation and the hydrodynamic pressure distribution. Taking the simulations as “ground truth,” a hydrodynamic bulge test is performed using the simulation data to ascertain the accuracy and the validity of the proposed methodology for estimating material properties. An error propagation analysis is performed via Monte Carlo simulation to characterize the susceptibility of the hydrodynamic bulge test estimates to noise. We find that, while a hydrodynamic bulge test is less accurate in characterizing material properties, it is less susceptible to noise, in the input (measured) variable, than a hydrostatic bulge test.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1557/mrc.2018.98
