skip to main content


Title: Elastic symmetry with beachball pictures
The elastic map, or generalized Hooke’s Law, associates stress with strain in an elastic material. A symmetry of the elastic map is a reorientation of the material that does not change the map. We treat the topic of elastic symmetry conceptually and pictorially. The elastic map is assumed to be linear, and we study it using standard notions from linear algebra—not tensor algebra. We depict strain and stress using the “beachballs” familiar to seismologists. The elastic map, whose inputs and outputs are strains and stresses, is in turn depicted using beachballs. We are able to infer the symmetries for most elastic maps, sometimes just by inspection of their beachball depictions. Many of our results will be familiar, but our versions are simpler and more transparent than their counterparts in the literature.  more » « less
Award ID(s):
1829447
NSF-PAR ID:
10247794
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Geophysical Journal International
ISSN:
0956-540X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract An elastic map $\mathbf {T}$ T describes the strain-stress relation at a particular point $\mathbf {p}$ p in some material. A symmetry of $\mathbf {T}$ T is a rotation of the material, about  $\mathbf {p}$ p , that does not change  $\mathbf {T}$ T . We describe two ways of inferring the group $\mathcal {S} _{ \mathbf {T} }$ S T of symmetries of any elastic map $\mathbf {T}$ T ; one way is qualitative and visual, the other is quantitative. In the first method, we associate to each $\mathbf {T}$ T its “monoclinic distance function” "Equation missing" on the unit sphere. The function "Equation missing" is invariant under all of the symmetries of  $\mathbf {T}$ T , so the group $\mathcal {S} _{ \mathbf {T} }$ S T is seen, approximately, in a contour plot of "Equation missing" . The second method is harder to summarize, but it complements the first by providing an algorithm to compute the symmetry group $\mathcal {S} _{ \mathbf {T} }$ S T . In addition to $\mathcal {S} _{ \mathbf {T} }$ S T , the algorithm gives a quantitative description of the overall approximate symmetry of  $\mathbf {T}$ T . Mathematica codes are provided for implementing both the visual and the quantitative approaches. 
    more » « less
  2. Abstract

    In linear elasticity, universal displacements for a given symmetry class are those displacements that can be maintained by only applying boundary tractions (no body forces) and for arbitrary elastic constants in the symmetry class. In a previous work, we showed that the larger the symmetry group, the larger the space of universal displacements. Here, we generalize these ideas to the case of anelasticity. In linear anelasticity, the total strain is additively decomposed into elastic strain and anelastic strain, often referred to as an eigenstrain. We show that theuniversality constraints(equilibrium equations and arbitrariness of the elastic constants) completely specify theuniversal elastic strainsfor each of the eight anisotropy symmetry classes. The corresponding universal eigenstrains are the set of solutions to a system of second-order linear PDEs that ensure compatibility of the total strains. We show that for three symmetry classes, namely triclinic, monoclinic, and trigonal, only compatible (impotent) eigenstrains are universal. For the remaining five classes universal eigenstrains (up to the impotent ones) are the set of solutions to a system of linear second-order PDEs with certain arbitrary forcing terms that depend on the symmetry class.

     
    more » « less
  3. 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
  4. There seems to be a basic misconception in several recent papers concerning the material symmetry of bodies in configurations that are pre-stressed. In this short paper we point to the source of the error and show that the material symmetry that is possible depends on the nature of the pre-stress. We also extend the results for material symmetry which have been well known within the context of simple elastic solids to the general class of simple materials. This generalization has relevance to the material symmetry of biological solids that are viscoelastic. 
    more » « less
  5. Abstract

    Numerical simulations using non-linear hyper-elastic material models to describe interactions between brain white matter (axons and extra cellular matrix (ECM)) have enabled high-fidelity characterization of stress-strain response. In this paper, a novel finite element model (FEM) has been developed to study mechanical response of axons embedded in ECM when subjected to tensile loads under purely non-affine kinematic boundary conditions. FEM leveraging Ogden hyper-elastic material model is deployed to understand impact of parametrically varying oligodendrocyte-axon tethering and analyze influence of aging material characteristics on stress propagation. In proposed FEM, oligodendrocyte connections to axons are represented via spring-dashpot model, such tethering technique facilitates contact definition at various locations, parameterize connection points and vary stiffness of connection hubs. Two FE submodels are discussed: 1) multiple oligodendrocytes arbitrarily tethered to the nearest axons, and 2) single oligodendrocyte tethered to all axons at various locations. Root mean square deviation (RMSD) were computed between stress-strain plots to depict trends in mechanical response. Axonal stiffness was found to rise with increasing tethering, indicating role of oligodendrocytes in stress redistribution. Finally, stress state results for aging axon material, with varying stiffnesses and number of connections in FEM ensemble have also been discussed to demonstrate gradual softening of tissues.

     
    more » « less