A Formulation of Anisotropic Elastic Damage Using Compact Tensor Formalism

An anisotropic elastic-damage model for initially-isotropic materials is presented. The model is based on a pseudo-logarithmic second-order damage tensor rate. To derive the complete expression of the tangent stiffness entering the rate constitutive law, various tensor operations and derivatives of tensor functions must be developed. Such derivations have been performed in compact form. Some useful tensor derivatives and a table of tensor algebra operations are given in Appendix. This note should interest engineering researchers involved in the development of constitutive models through tensor formalism. This revised version was published online in July 2006 with corrections to the Cover Date.     

各向加权异性张力有限元法及其应用

本文利用各向异性张力有限元的特点设计了灵活而有效的各向加权异性有限元算法，采用这种方法，计算了离心叶轮内部准三元流场。     

Equivalent Diffusion in a Thin Film Flow with a Non-One-Dimensional Velocity Field and Anisotropic Diffusion Tensor

For describing the mass transfer processes in channels, Taylor's dispersion theory is widely used. This theory makes it possible, with asymptotic rigor, to replace the complete diffusion (heat conduction) equation with a convective term that depends on the coordinate transverse to the flow by an effective diffusion (dispersion) equation with constant coefficients, averaged over the channel cross-section. In numerous subsequent studies, Taylor's theory was generalized to include more complex situations, and novel algorithms for constructing the dispersion equations were proposed. For thin film flows a theory similar to Taylor's leads to a matrix of dispersion coefficients.In this study, Taylor's theory is extended to film flows with a non-one-dimensional velocity field and anisotropic diffusion tensor. These characteristics also depend to a considerable extent on the spatial coordinates and time. The dispersion equations obtained can be simplified in regions in which the effective diffusion coefficient tensor changes sharply.     

Viscous Fingering Instability in Porous Media: Effect of Anisotropic Velocity-Dependent Dispersion Tensor

The viscous fingering of miscible flow displacements in a homogeneous porous media is examined to determine the effects of an anisotropic dispersion tensor on the development of the instability. In particular, the role of velocity-dependent transverse and longitudinal dispersions is investigated through linear stability analysis and nonlinear simulations. It is found that an isotropic velocity-dependent dispersion tensor does not affect substantially the development of the instability and effectively has the same effect as molecular diffusion. On the other hand, an anisotropic velocity-dependent dispersion tensor results in different instability characteristics and more intricate finger structures. It is shown that anisotropic dispersion has profound effects on the development of the fingers and on the mechanisms of interactions between neighboring fingers. The development of the new finger structures is explained by examining the velocity field and characterized qualitatively through a spectral analysis of the average concentration and an analysis of the variations of the sweep efficiency and relative contact area.     

Drained and Undrained Poroelastic Properties of Healthy and Pathological Bone: A Poro-Micromechanical Investigation

Poro-micromechanics allows for the quantification of poroelastic properties such as the Biot and Skempton coefficients, once a continuum micromechanics model for the material under consideration has been developed and validated. Employing such a model for the transversely isotropic elasticity of cortical and trabecular bone, we determine the tensors of Biot and Skempton coefficients as functions of the volume fractions of mineral, collagen, and the micropore space (Haversian and Volkmann canals, and the inter-trabecular space). Increase of microporosity, as experienced in osteoporosis, as well as decrease of mineral content, as experienced in osteomalacia, lead to an increase of Biot and Skempton coefficients, i. e. to magnification of the mechanical role of the marrow filling the micropore space. For quantification of the marrow pressure rise upon downfall, undrained conditions are appropriate, as can be shown by model predictions of non-destructive impact experiments.     

基于Oda裂隙结构张量的流纹岩各向异性试验研究

岩石内天然存在长度、倾角和形态不同的裂隙,造成岩石的各向异性特征。为揭示岩石内天然随机裂隙发育特征对岩石物理力学特性的影响规律,以泥巴山隧址区采集裂隙性流纹岩为研究对象,首先对试样裂隙进行素描统计分析;然后基于Oda裂隙结构张量,获得天然随机分布裂隙的几何统计参数;最后对裂隙性流纹岩试样分别进行单轴和常规三轴压缩试验,得到不同应力路径下流纹岩的应力-应变曲线及物理力学参数。分析Oda裂隙结构张量定义的各向异性参数与试验获得的力学参数之间的规律,研究结果表明:(1)Oda裂隙结构张量适用于天然随机分布裂隙的几何统计分析,各向异性参数A(F)越大,裂隙优势方向越明显;(2)单轴压缩下,随着各向异性参数I1和A(F)的增大,流纹岩各向异性程度增大,弹性模量减小,泊松比增大;(3)常规三轴压缩下,流纹岩弹性模量和泊松比随各向异性参数改变的规律较不明显,Oda裂隙结构张量不再适用。     

A direct representation of the rotation tensor in polar decomposition of deformation gradient and its applications

Through a variational approach, an explicit connection between the additive and the polar decompositions of deformation gradient has been established. An exact formula for determining the rotation tensor in polar decomposition is obtained. The formula is fundamental in continuum mechanics and can be used to separate the rotation and the pure strain in deformation, by which various approximate expressions can be easily obtained.     

Affine Transforms of Three-Dimensional Anisotropic Media and Explicit Formulas for Fundamental Matrices

The concept of classes of algebraically equivalent anisotropic three-dimensional media is introduced: elastic fields in such media are related by simple algebraic expressions. An explicit formula is obtained for a fundamental matrix with ten free constants in the elastic modulus tensor (rather than for five free constants, as in the famous case of transverse isotropy). A hypothesis is formulated and several questions are posed, which are related to the notion of algebraic equivalence under discussion. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 2, pp. 95–102, March–April, 2006.     

A Concise Proof of the Representation Theorem for Fourth-Order Isotropic Tensors

We present a new proof of the representation theorem for fourth-order isotropic tensors that does not assume the tensor to have major or minor symmetries at the outset.     

Anisotropic Elastic Materials that Uncouple All Three Displacement Components, and Existence of One-Displacement Green's Function

It was shown in an earlier paper that, under a two-dimensional deformation, there are anisotropic elastic materials for which the antiplane displacement u ₃ and the inplane displacements u ₁, u ₂ are uncoupled but the antiplane stresses σ₃₁, σ₃₂ and the inplane stresses σ₁₁, σ₁₂, σ₂₂ remain coupled. The conditions for this to be possible were derived, but they have a complicated expression. In this paper new and simpler conditions are obtained, and a general anisotropic elastic material that satisfies the conditions is presented. For this material, and for certain monoclinic materials with the symmetry plane at x ₃ = 0, we show that the unnormalized Stroh eigenvectors a _k for k = 1, 2, 3 are all real. The matrix A =[a ₁, a ₂, a ₃] is a unit matrix when the material has a symmetry plane at x ₂ = 0. Thus any one of the u ₁, u ₂, u ₃ can be the only nonzero displacement, and the solution is a one-displacement field. Application to the Green's function due to a line of concentrated force f and a line dislocation with Burgers vector v in the infinite space, the half-space with a rigid boundary, and the infinite space with an elliptic rigid inclusion shows that one can indeed have a one-displacement field u ₁, u ₂ or u ₃. One can also have a two-displacement field polarized on a plane other than the (x ₁, x ₂)-plane. The material that uncouples u ₁, u ₂, u ₃ is not as restrictive as one might have thought. It can be triclinic, monoclinic, orthotropic, tetragonal, transversely isotropic, or cubic. However, it cannot be isotropic. This revised version was published online in August 2006 with corrections to the Cover Date.     

Similarity Laws for the Velocity Distribution and the Reynolds Tensor Components in the Wall Region of a Turbulent Boundary Layer with Injection and Suction

Similarity laws of the distributions of the average velocity, tangential stress, and mean-square transverse velocity fluctuation are established in an intermediate zone of a turbulent boundary layer with injection and suction. This zone is located in the neighborhood of the wall outside the viscous sublayer. The similarity relationship for the velocity profile is a generalization of the well-known logarithmic law to include the case of the presence of a mass flow at the wall.     

Determination of the Rotation Tensor in the Polar Decomposition

An explicit representation for the rotation tensor which contains the lower powers of deformation gradient is proposed and used to evaluate the angle and axis of the rotation tensor. Some related equations about the rotation tensor are established. Through the approximate analysis, the relation between the S-R decomposition and the polar decomposition is examined. This revised version was published online in July 2006 with corrections to the Cover Date.     

Objective Tensor Rates and Applications in Formulation of Hyperelastic Relations

Following Ogden, a class of objective (Lagrangian and Eulerian) tensors is identified among the second-rank tensors characterizing continuum deformation, but a more general definition of objectivity than that used by Ogden is introduced. Time rates of tensors are determined using convective rates. Sufficient conditions of objectivity are obtained for convective rates of objective tensors. Objective convective rates of strain tensors are used to introduce pairs of symmetric stress and strain tensors conjugate in a generalized sense. The classical definitions of conjugate Lagrangian (after Hill) and Eulerian (after Xiao et al.) stress and strain tensors are particular cases of the definition of conjugacy of stress and strain tensors in the generalized sense used in the present paper. Pairs of objective stress and strain tensors conjugate in the generalized sense are used to formulate constitutive relations for a hyperelastic medium. A family of objective generalized strain tensors is introduced, which is broader than Hill’s family of strain tensors. The basic forms of the hyperelastic constitutive relations are obtained with the aid of pairs of Lagrangian stress and strain tensors conjugate after Hill (the strain tensors in these pairs belong to the family of generalized strain tensors). A method is presented for generating reduced forms of the constitutive relations with the aid of pairs of Lagrangian and Eulerian stress and strain tensors conjugate in the generalized sense which are obtained from pairs of Lagrangian tensors conjugate after Hill by mapping tensor fields on one configuration of a deformable body to tensor fields on another configuration.      

On the Wrinkling of Anisotropic Elastic Membranes

The existence and uniqueness of a relaxed energy function for anisotropic elastic membranes is established. It is shown that the extra strains during wrinkling abide by a normality rule akin to that of metal plasticity, and that throughout such wrinkling processes both the strain energy and the Piola-Kirchhoff stresses remain constant. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Remarkable Nature of Cylindrically Anisotropic Elastic Materials Exemplified by an Anti-Plane Deformation

A material is cylindrically anisotropic when its elastic moduli referred to a cylindrical coordinate system are constants. Examples of cylindrically anisotropic materials are tree trunks, carbon fibers [1], certain steel bars, and manufactured composites [2]. Lekhnitskii [3] was the first one to observe that the stress at the axis of a circular rod of cylindrically monoclinic material can be infinite when the rod is subject to a uniform radial pressure (see also [4]). Ting [5] has shown that the stress at the axis of the circular rod can also be infinite under a torsion or a uniform extension. In this paper we first modify the Lekhnitskii formalism for a cylindrical coordinate system. We then consider a wedge of cylindrically monoclinic elastic material under anti-plane deformations. The stress singularity at the wedge apex depends on one material parameter γ. For a given wedge angle α, one can choose a γ so that the stress at the wedge apex is infinite. The wedge angle 2α can be any angle. It need not be larger than π, as is the case when the material is homogeneously isotropic or anisotropic. In the special case of a crack (2α=2π) there can be more than one stress singularity, some of them are stronger than the square root singularity. On the other hand, if γ < there is no stress singularity at the wedge apex for any wedge angle, including the special case of a crack. The classical paradox of Levy [6] and Carothers [7] for an isotropic elastic wedge also appears for a cylindrically anisotropic elastic wedge. There can be more than one critical wedge angle and, again, the critical wedge angle can be any angle. This revised version was published online in August 2006 with corrections to the Cover Date.     

Stress Singularity in an Elastic Cylinder of Cylindrically Anisotropic Materials

The issue of stress singularity in an elastic cylinder of cylindrically anisotropic materials is examined in the context of generalized plane strain and generalized torsion. With a viewpoint that the singularity may be attributed to a conflicting definition of anisotropy at r=0, we study the problem through a compound cylinder in which the outer cylinder is cylindrically anisotropic and the core is transversely isotropic. By letting the radius of the core go to zero, the cylinder becomes one with the central axis showing no conflict in the radial and tangential directions. Closed-form solutions are derived for the cylinder under pressure, extension, torsion, rotation and a uniform temperature change. It is found that the stress is bounded everywhere, and singularity does not occur if the anisotropy at r=0 is defined appropriately. This revised version was published online in July 2006 with corrections to the Cover Date.     

Measurement of the Gruneisen Coefficient of Some Anisotropic Carbon Materials

Results of measurement of the Gruneisen coefficient of two anisotropic materials (4KMS carbon-carbon composite and UPV-1 pyrolytic carbon) are presented. The values of the Gruneisen coefficient are calculated on the basis of the measured amplitude of mechanical stresses arising in the material with absorption of electron-beam energy. The results for the 4KMS composite include not only a negative value of the Gruneisen coefficient in the range of electron energy density of several tens of joules per square centimeter but also its anomalously low value as compared to other materials. Significant decay of the stress wave propagating over this material is noted. The Gruneisen coefficient of the UPV-1 sample varies depending on the sample orientation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 171–179, November–December, 2005.     

The Onset of Convection in an Anisotropic Porous Layer Using a Thermal Non-Equilibrium Model

The stability of a horizontal fluid saturated anisotropic porous layer heated from below and cooled from above is examined analytically when the solid and fluid phases are not in local thermal equilibrium. Darcy model with anisotropic permeability is employed to describe the flow and a two-field model is used for energy equation each representing the solid and fluid phases separately. The linear stability theory is implemented to compute the critical Rayleigh number and the corresponding wavenumber for the onset of convective motion. The effect of thermal non-equilibrium and anisotropy in both mechanical and thermal properties of the porous medium on the onset of convection is discussed. Besides, asymptotic analysis for both very small and large values of the interphase heat transfer coefficient is also presented. An excellent agreement is found between the exact and asymptotic solutions. Some known results, which correspond to thermal equilibrium and isotropic porous medium, are recovered in limiting cases.     

Coriolis Effect on the Stability of Centrifugally Driven Convection in a Rotating Anisotropic Porous Layer Subjected to Gravity

We investigate natural convection in a fluid saturated rotating anisotropic porous layer subjected to centrifugal gravitational and Coriolis body forces. The Darcy model (including the centrifugal, gravitational and Coriolis terms; and permeability anisotropy effects) and a modified energy equation (including the effects of thermal anisotropy) is used in the current analysis. The linear stability theory is used to evaluate the critical Rayleigh number for the onset of convection in the presence of thermal and mechanical anisotropy. It is shown that the preferred solution comprises roll cells aligned parallel to the vertical z-axis. As a result, it is found that the Coriolis acceleration (or Taylor number) and the gravitational term play no role in the stability of convection.