Partial Uniqueness and Obstruction to Uniqueness in Inverse Problems for Anisotropic Elastic Media

We consider the inverse problem of identifying the density and elastic moduli for three-dimensional anisotropic elastic bodies, given displacement and traction measurements made at their surface. These surface measurements are modelled by the dynamic Dirichlet-to-Neumann map on a finite time interval. For linear or nonlinear anisotropic hyperelastic bodies we show that the displacement-to-traction surface measurements do not change when the density and elasticity tensor in the interior are transformed tensorially by a change of coordinates fixing the surface of the body to first order. Our main tool, a new approach in inverse problems for elastic media, is the representation of the equations of motion in a covariant form (following Marsden and Hughes, 1983) that preserves the underlying physics.In the case of classical linear elastodynamics we then investigate how the type of anisotropy changes under coordinate transformations. That is, we analyze the orbits of general linear, anisotropic elasticity tensors under the action by pull-back of diffeomorphisms that fix the surface of the elastic body to first order, and derive a pointwise characterization of parts of the orbits under this action. For example, we show that the orbit of isotropic elastic media, at any point in the body, consists of some transversely isotropic and some orthotropic elastic media. We then derive the first uniqueness result in the inverse problem for anisotropic media using surface displacement-traction data: uniqueness of three elastic moduli for tensors in the orbit of isotropic elasticity tensors. ^†Partially supported by an MSRI Postdoctoral Fellowship. Research at MSRI is supported in part by NSF grant DMS-9850361. This work was conducted while the first author was a Gibbs Instructor at Yale University. ^‡Partially supported by an MSRI Postdoctoral Fellowship, and by NSF grant DMS-9801664 (9996350).     

A novel internal dissipation inequality by isotropy and its implication for inelastic constitutive characterization

In the present work a novel inelastic deformation caused internal dissipation inequality by isotropy is revealed. This inequality has the most concise form among a variety of internal dissipation inequalities, including the one widely used in constitutive characterization of isotropic finite strain elastoplasticity and viscoelasticiy. Further, the evolution term describing the difference between the rate of deformation tensor and the “principal rate” of the elastic logarithmic strain tensor is set, according to the standard practice by isotropy, to equal a rank-two isotropic tensor function of the corresponding branch stress, with the tensor function having an eigenspace identical to the eigenspace of the branch stress tensor. Through that a general form of evolution equation for the elastic logarithmic strain is formulated and some interesting and important results are derived. Namely, by isotropy the evolution of the elastic logarithmic strain tensor is embodied separately by the evolutions of its eigenvalues and eigenprojections, with the evolution of the eigenprojections driven by the rate of deformation tensor and the evolution of the eigenvalues connected to specific material behavior. It can be proved that by isotropy the evolution term in the present dissipation inequality stands for the essential form of the evolution term in the extensively applied dissipation inequality.     

Crack propagation in an orthotropic medium with coupled elastodynamic properties

In this paper the Mode-I elastodynamic problem of a crack propagating in an orthotropic medium is studied under the condition that the matrix of elastodynamic coefficients has repeated eigenvalues. It is shown that the crack is constrained in an elastodynamic state which is defined through a compulsory condition coupling its velocity with the elastic parameters of the orthotropic medium. The dynamic stress and displacement components ahead of the crack tip as well as the energy release rate are explicitly obtained.     

Mindlin second-gradient elastic properties from dilute two-phase Cauchy-elastic composites Part II: Higher-order constitutive properties and application cases

Starting from a Cauchy elastic composite with a dilute suspension of randomly distributed inclusions and characterized at first-order by a certain discrepancy tensor (see part I of the present article), it is shown that the equivalent second-gradient Mindlin elastic solid: (i) is positive definite only when the discrepancy tensor is negative defined; (ii) the non-local material symmetries are the same of the discrepancy tensor, and (iii) the non-local effective behaviour is affected by the shape of the RVE, which does not influence the first-order homogenized response. Furthermore, explicit derivations of non-local parameters from heterogeneous Cauchy elastic composites are obtained in the particular cases of: (a) circular cylindrical and spherical isotropic inclusions embedded in an isotropic matrix, (b) n-polygonal cylindrical voids in an isotropic matrix, and (c) circular cylindrical voids in an orthotropic matrix.     

Static equilibrium of an elastic orthotropic medium with an arbitrarily oriented triaxial ellipsoidal inclusion

The stress state of an elastic orthotropic medium with an arbitrarily oriented triaxial ellipsoidal inclusion is analyzed. A solution is obtained using the triple Fourier transform and the Fourier-transformed Green’s function for an infinite anisotropic medium. The high efficiency of the approach is demonstrated by solving the problem for a transversely isotropic material with a spheroidal cavity for which the exact solution is known. A numerical analysis is conducted to study the stress distribution over the surface of the inclusion with different orientations in the orthotropic space. It is revealed that in some cases the orientation of the inclusion has a strong effect on the stress concentration __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 4, pp. 55–61, April 2007.     

An ‘extended’ volumetric/deviatoric formulation of anisotropic damage based on a pseudo-log rate

Following a framework of elastic degradation and damage previously proposed by the authors, an ‘extended’ formulation of orthotropic damage in initially isotropic materials, based on volumetric/deviatoric decomposition, is presented. The formulation is founded on the concept of energy equivalence and makes use of second-order symmetric tensor damage variables. It is characterized by fourth-order damage-effect tensors (relating nominal to effective stresses and strains) built from the underlying second-order damage tensors and decomposed in product-form in isotropic and anisotropic parts. The formulation is developed in two steps. First, secant relations are established. In the isotropic case, the model embeds a path parameter allowing to range between pure volumetric to pure deviatoric damage. With the two undamaged material constants this makes a total of three constant parameters plus an evolving scalar damage variable, giving rise to a four-parameter model with two varying isotropic material coefficients. In the anisotropic case, the model is still characterized by the same three material constants plus three evolving variables which are the principal values of a second-order damage tensor. This leads to a six-parameter restricted form of orthotropic damage. In the second step, damage evolution rules are formulated in terms of a pseudo-logarithmic rate of damage. This allows to define meaningful conjugate forces that constitute a feasible space in which loading functions and damage evolution rules can be defined. The present ‘extended’ formulation is closed by the derivation of the tangent stiffness.     

Laws of viscoplastic fluid flow in anisotropic porous and fractured media

In invariant tensor form, the laws of viscoplastic fluid flow are formulated for capillary and fractured media with a periodic microstructure that has orthotropic and transversely isotropic symmetry in the flow properties. An analysis of the laws of viscoplastic fluid flow in transversely isotropic and orthotropic porous and fractured media shows that in formulating the equations it is necessary to distinguish between the permeability tensor and the limiting gradient tensor, which may differ in the symmetry of the flow characteristics, and that the flow law is multivariant and admits one-, two-, and three-dimensional flows.     

含正交排列夹杂和缺陷材料的等效弹性模量和损伤

研究含正交排列夹杂和缺陷材料的等效弹性模量和损伤,推导了以Eshelby－Mori－Tanaka方法求解多相各向异性复合材料等效弹性模量的简便计算公式,针对含三相正交椭球状夹杂的正交各向异性材料,得到了由细观参量（夹杂的形状、方位和体积分数）表示的等效弹性模量的解析表达式．在此基础上,提出了一个宏细观结合的正交各向异性损伤模型,从而建立了以细观量为参量的含损伤材料的应力应变关系．最后,对影响材料损伤的细观结构参数进行了分析．     

On a full-strength juncture between an isotropic matrix and an orthotropic inclusion under a triaxial load

The possibility exists for the full-strength state of the juncture between an isotropic medium and an orthotropic inclusion when the phases are in ideal contact and the stresses are tensile (compressive). It is found that an ellipsoid, whose semiaxes depend heavily on the elastic characteristics of the matrix and inclusion and the parameters of external loading, is such a form. A system of transcendental equations is obtained to search for the unknown semi-axes. Previous data confirm the results obtained. Numerical investigations are performed, and interrelations are established between geometrical parameters and loading conditions and the properties of the phases for transversally isotropic and orthotropic inclusions. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev. Translated from Prikladnaya Mekhanika Vol. 36, No. 1, pp. 95–102, January, 2000     

Effects of thermal relaxations on thermoelastic interactions in an infinite orthotropic elastic medium with a cylindrical cavity

Thermoelastic interactions in an infinite orthotropic elastic medium with a cylindrical cavity are studied. The cavity surface is subjected to ramp-type heating of its internal boundary, which is assumed to be traction free. Lord–Shulman and Green–Lindsay models for the generalized thermoelasticity theories are selected since they allow for second-sound effects and reduce to the classical model for an appropriate choice of the parameters. The temperature, radial displacement, radial stress, and hoop stress distributions are computed numerically using the finite-element method (FEM). The results are presented graphically for different values of the thermal relaxation times using the three different theories of generalized thermoelasticity. Excellent agreement is found between the finite-element analysis and analytical and classical solutions.     

Nonlinear vibration analysis of laminated composite Mindlin micro/nano-plates resting on orthotropic Pasternak medium using DQM

The nonlocal nonlinear vibration analysis of embedded laminated microplates resting on an elastic matrix as an orthotropic Pasternak medium is investigated. The small size effects of micro/nano-plate are considered based on the Eringen nonlocal theory. Based on the orthotropic Mindlin plate theory along with the von Kármán geometric nonlinearity and Hamilton's principle, the governing equations are derived. The differential quadrature method (DQM) is applied for obtaining the nonlinear frequency of system. The effects of different parameters such as nonlocal parameters, elastic media, aspect ratios, and boundary conditions are considered on the nonlinear vibration of the micro-plate. Results show that considering elastic medium increases the nonlinear frequency of system. Furthermore, the effect of boundary conditions becomes lower at higher nonlocal parameters.     

A NOTE ON MORI-TANAKA＇S METHOD

Explicit expressions of Mori-Tanaka＇s tensor for a transversely isotropic fiber rein- forced UD composite are presented. Closed-form formulae for the effective elastic properties of the composite are obtained. In a 3D sense, the resulting compliance tensor of the composite is symmetric. Nevertheless, the 2D compliance tensor based on a deteriorated Mori-Tanaka＇s tensor is not symmetric. Nor is the compliance tensor defined upon a deteriorated 2D Eshelby＇s tensor. The in-plane effective elastic properties given by those three approaches are different. A detailed comparison between the predicted results obtained from those approaches with experimental data available for a number of UD composites is made.     

Basic Issues Concerning Finite Strain Measures and Isotropic Stress-Deformation Relations

It is indicated that the commonly-used Rivlin–Ericksen representation formula for isotropic tensor functions exhibits some properties that might be undesirable for its reasonable and effective applications. Towards clarification and improvement, a set of three mutually orthogonal tensor generators is introduced to achieve an alternative representation formula for isotropic symmetric tensor-valued functions of a symmetric tensor. This representation formula enables us to express the unknown representative coefficients in terms of simple, explicit tensorial inner products of the argument tensor and the value tensor without involving their eigenvalues. In particular, the tensorial interpolation expressions thus obtained assume a unified form for the three different cases of coalescence of the eigenvalues of the argument tensor. Moreover, each summand in the alternative representation formula is shown to inherit the continuity and differentiability properties of the represented isotropic tensor function. These results are used to study some basic issues concerning finite strain measures and stress-deformation relations of isotropic materials, such as continuity and differentiability properties of the representation, determination of the representative coefficients in terms of experimental data for stress and deformation tensors, and computations of finite strain measures. This revised version was published online in June 2006 with corrections to the Cover Date.     

Further improvement on fundamental solutions of plane problems for orthotropic materials

On the basis of the ecisting fundamental solutions of displacements, further improvement is made, and then the general fundamental solutions of both plane elastic and plane plastic problems for orthotropic materials are obtained. Two parameters based on material constants α₁ = α₁ are used to derive the relevant expressions in a real variable form. Additionally an analytical method of solving the singular integral for the internal stresses is introduced, and the corresponding results are given. If α₁ = α₁ = 1, all the expressions obtained for orthotropy can be reduced to the corresponding ones for isotropy. Because all these expressions and results can be directly used for both isotropic problems and orthotropic problems, it is convenient to use them in engineering with the boundary element method (BEM).     

Two-dimensional transient dynamic response of orthotropic layered media

In this study, two-dimensional transient dynamic response of orthotropic plane layered media is investigated. The plane multilayered media consist of N different generally orthotropic, homogeneous and linearly elastic layers with different ply angles. In the generally orthotropic layer, representing a ply reinforced by unidirectional fibers with an arbitrary orientation angle, the principal material directions do not coincide with body coordinate axes. The solution is obtained by employing a numerical technique which combines the use of Fourier transform with the method of characteristics. The numerical results are displayed in curves denoting the variations of stress and displacement components with time at different locations. These curves clearly reveal, in wave profiles, the scattering effects caused by the reflections and refractions of waves at the boundaries and at the interfaces of the layers, and also the effects of anisotropy caused by fiber orientation angle. The curves properly predict the sharp variations in the response at the neighborhood of the wave fronts, which shows the power of the numerical technique employed in the study. By suitably adjusting the elastic constants, the results for multilayered media with transversely isotropic layers, or layers with cubic symmetry, or isotropic layers can easily be obtained from the general formulation. Furthermore, solutions for some special cases, including Lamb’s problem for an elastic half-space, are obtained and compared with the available solutions in the literature and very good agreement is found. Preliminary version presented at the Second International Congress on Mechatronics (MECH2K3), Graz, Austria, July 14-17, 2003.     

Determining equations of two-phase flows through anisotropic porous media

A new interpretation of the concept of relative phase permeability is given. Relative phase permeabilities are represented in the form of fourth-rank tensors. It is shown that in the case of anisotropic porous media functions depending not only on the saturation but also on the anisotropy parameters represented in the form of ratios of the principal values of the absolute permeability coefficient tensor correspond to the classical representation of the relative phase permeabilities. For a two-phase flow in anisotropic porous media with orthotropic and transversely-isotropic symmetry a generalized two-term Darcy’s law is analyzed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 87–94, March–April, 1998. The work was carried out with support from the Russian Foundation for Fundamental Research (project No. 96-01-00623).     

Orthotropic elastic–plastic material model for paper materials

Paper and paperboard generally exhibit anisotropic and non-linear mechanical material behaviour. In this work, the development of an orthotropic elastic–plastic constitutive model, suitable for modelling of the material behaviour of paper is presented. The anisotropic material behaviour is introduced into the model by orthotropic elasticity and an isotropic plasticity equivalent transformation tensor. A parabolic stress–strain relation is adopted to describe the hardening of the material. The experimental and numerical procedures for evaluation of the required material parameters for the model are described. Uniaxial tensile testing in three different inplane material directions provides the calibration of the material parameters under plane stress conditions. The numerical implementation of the material model is presented and the model is shown to perform well in agreement with experimentally observed mechanical behaviour of paper.     

Vibrations of plates with initial geometrical periodical imperfections interacting with a periodic elastic foundation

Summary  The aim of this contribution is to formulate an engineering theory describing the effect of periodic structure of the elastic foundation on the dynamic behaviour of a plate with geometrical periodical imperfections, Fig. 1. It is shown that the obtained solutions depend on the structure length parameter l of the elastic foundation and the wave length of imperfections. The general results are illustrated by an example and compared with the known theories of orthotropic plates. Received 7 February 1999; accepted for publication 6 October 1999