Stroh formulation for a generally constrained and pre-stressed elastic material

The Stroh formalism is essentially a spatial Hamiltonian formulation and has been recognized to be a powerful tool for solving elasticity problems involving generally anisotropic elastic materials for which conventional methods developed for isotropic materials become intractable. In this paper we develop the Stroh/Hamiltonian formulation for a generally constrained and prestressed elastic material. We derive the corresponding integral representation for the surface-impedance tensor and explain how it can be used, together with a matrix Riccati equation, to calculate the surface-wave speed. The proposed algorithm can deal with any form of constraint, pre-stress, and direction of wave propagation. As an illustration, previously known results are reproduced for surface waves in a pre-stressed incompressible elastic material and an unstressed inextensible fibre-reinforced composite, and an additional example is included analyzing the effects of pre-stress upon surface waves in an inextensible material.     

An explicit expression for the surface-impedance matrix of a generally anisotropic incompressible elastic material in a state of plane strain

Surface and interfacial impedance matrices play an important role in the construction of Green's functions, the analysis of surface and interfacial waves and the stability assessment of pre-stressed half-spaces or joined half-spaces. This paper studies these matrices for generally anisotropic pre-stressed incompressible elastic materials. It is shown that the surface-impedance matrix satisfies a simple matrix equation which, for plane-strain deformations, can be solved exactly. As a result, explicit secular equations for surface and interfacial wave speeds and explicit wrinkling/buckling conditions for pre-stressed half-spaces and joined half-spaces are obtained. It is also shown that the plane-strain surface-wave problem is mathematically identical to the edge-wave problem for thin elastic plates. Thus, the uniqueness of surface-wave speed is settled by drawing upon a recent proof of the uniqueness of edge-wave speed. Examples are used to show that it is straightforward to solve the secular equations based on the given formulae either exactly (where possible) or numerically.     

A Representation Theorem for Material Tensors of?Weakly-Textured Polycrystals and?Its Applications in?Elasticity

Material tensors pertaining to polycrystalline aggregates should manifest also the influence of crystallographic texture on the material properties in question. In this paper we make use of tensors which form bases of irreducible representations of the rotation group and prove a representation theorem by which a given material tensor of a weakly-textured polycrystal is expressed as a linear combination of an orthonormal set of irreducible basis tensors, with the components given explicitly in terms of texture coefficients and a set of undetermined material parameters. Once the irreducible basis tensors that appear in the formula are determined, the representation formula, which is valid for all texture and crystal symmetries, will delineate quantitatively the effect of crystallographic texture on the material tensor in question. We present an integral formula and an orthonormalization process which serve as the basis for a procedure to determine explicitly the irreducible basis tensors required in the representation formula. For applications we determine a set of irreducible basis tensors for the elasticity tensor and a set for fourth-order tensors that define constitutive equations in incompressible elasticity and Hill’s quadratic yield functions in plasticity. We show that orientation averaging of a tensor can be done easily if we have in hand a set of irreducible basis tensors for the decomposition of the tensor in question. As illustration we derive a formula, which is valid for all texture and crystal symmetries, for the elasticity tensor under the Voigt model.     

An Invariant-Based Ogden-Type Model for Incompressible Isotropic Hyperelastic Materials

The Ogden model for an incompressible isotropic hyperelastic material is versatile enough to match complicated data for rubber-like materials at large deformations. However, the tensorial expression for the Cauchy stress in the Ogden model requires determination of the eigenvalues and eigenvectors of the left Cauchy-Green deformation tensor \(\mathbf{B}\). The objective of this paper is to propose an invariant-based Ogden-type model for isotropic incompressible hyperelastic materials. The strain energy function in this new model depends on classical invariants of \(\mathbf{B}\) and the Cauchy stress tensor can be expressed directly in terms of the tensor \(\mathbf{B}\) without need for its spectral form. Examples show that this new Ogden-type model retains the versatility of the original Ogden model in characterizing material response.     

Constitutive Models for Compressible Nonlinearly Elastic Materials with Limiting Chain Extensibility

Constitutive models are proposed for compressible isotropic hyperelastic materials that reflect limiting chain extensibility. These are generalizations of the model proposed by Gent for incompressible materials. The goal is to understand the effects of limiting chain extensibility when the compressibility of polymeric materials is taken into account. The basic homogeneous deformation of simple tension is considered and simple closed-form relations for the deformation characteristics are obtained for slightly compressible materials. An explicit first-order approximation is obtained for the lateral contraction and for the Poisson function in terms of the axial extension which is shown to be valid for each of two specific compressible versions of the Gent model. One of the main results obtained is that the effect of limiting chain extensibility is to stiffen the material relative to the neo-Hookean compressible case. Mathematics Subject Classifications (2000) 74B20, 74G55.     

On discontinuous strain fields in incompressible finite elastostatics

Piece-wise homogeneous three-dimensional deformations in incompressible materials in finite elasticity are considered. The emergence of discontinuous strain fields in incompressible materials is studied via singularity theory. Since the simplest singularities, including Maxwell’s sets, are the cusp singularities, cusp conditions for the total energy function of homogeneous deformations for incompressible materials in finite elasticity will be derived, compatible with strain jumping. The proposed method yields simple criteria for the study of discontinuous deformations in three-dimensional problems and for any homogeneous incompressible material. Furthermore the homogeneous stress tensor is also not restricted. Neither fictitious nor simplified constitutive relations are invoked. The theory is implemented in a simple shearing problem.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part I Non-isothermal flow

Tensor representation theory is used to derive an explicit algebraic model that consists of an explicit algebraic stress model (EASM) and an explicit algebraic heat flux model (EAHFM) for two-dimensional (2-D) incompressible non-isothermal turbulent flows. The representation methodology used for the heat flux vector is adapted from that used for the polynomial representation of the Reynolds stress anisotropy tensor. Since the methodology is based on the formation of invariants from either vector or tensor basis sets, it is possible to derive explicit polynomial vector expansions for the heat flux vector. The resulting EAHFM is necessarily coupled with the turbulent velocity field through an EASM for the Reynolds stress anisotropy. An EASM has previously been derived by Jongen and Gatski [10]. Therefore, it is used in conjunction with the derived EAHFM to form the explicit algebraic model for incompressible 2-D flows. This explicit algebraic model is analyzed and compared with previous formulations including its ability to approximate the commonly accepted value for the turbulent Prandtl number. The effect of pressure-scrambling vector model calibration on predictive performance is also assessed. Finally, the explicit algebraic model is validated against a 2-D homogeneous shear flow with a variety of thermal gradients. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University     

Stability of thick spherical shells

The pressure-radius relation of spherical rubber balloons has been derived and its stability behavior investigated before. In this work, we show that similar results remain valid for thick spherical shells of Mooney-Rivlin materials. In addition, we show that eversion of a spherical shell is possible for any incompressible isotropic materials if the shell is not too thick.     

对不可压黏弹材料非线性本构关系的探讨

运用张量分析,证明了适合描述不可压缩黏弹性材料非线性本构关系的3种理论：BKZ理论、Lianis理论和Chirstensen理论,运用到各向同性材料中它们的单轴拉伸具有一样的表现形式,同时也验证了3种理论在解释非线性本构关系时具有统一性．     

Superposition of Generalized Plane Strain on Anti-Plane Shear Deformations in Isotropic Incompressible Hyperelastic Materials

The purpose of this research is to investigate the basic issues that arise when generalized plane strain deformations are superimposed on anti-plane shear deformations in isotropic incompressible hyperelastic materials. Attention is confined to a subclass of such materials for which the strain-energy density depends only on the first invariant of the strain tensor. The governing equations of equilibrium are a coupled system of three nonlinear partial differential equations for three displacement fields. It is shown that, for general plane domains, this system decouples the plane and anti-plane displacements only for the case of a neo-Hookean material. Even in this case, the stress field involves coupling of both deformations. For generalized neo-Hookean materials, universal relations may be used in some situations to uncouple the governing equations. It is shown that some of the results are also valid for inhomogeneous materials and for elastodynamics.     

A class of universal relations in isotropic elasticity theory

A class of universal relations for isotropic elastic materials is described by the tensor equationTB = BT. This simple rule yields at most three component relations which are the generators of many known universal relations for isotropic elasticity theory, including the well-known universal rule for a simple shear. Universal relations for four families of nonhomogeneous deformations known to be controllable in every incompressible, homogeneous and isotropic elastic material are exhibited. These same universal relations may hold for special compressible materials. New universal relations for a homogeneous controllable shear, a nonhomogeneous shear, and a variable extension are derived. The general universal relation for an arbitrary isotropic tensor function of a symmetric tensor also is noted.     

Finite strain––isotropic hyperelasticity

This paper presents a strain energy density for isotropic hyperelastic materials. The strain energy density is decomposed into a compressible and incompressible component. The incompressible component is the same as the generalized Mooney expression while the compressible component is shown to be a function of the volume invariant J only. The strain energy density proposed is used to investigate problems involving incompressible isotropic materials such as rubber under homogeneous strain, compressible isotropic materials under high hydrostatic pressure and volume change under uniaxial tension. Comparison with experimental data is good. The formulation is also used to derive a strain energy density expression for compressible isotropic neo-Hookean materials. The constitutive relationship for the second Piola–Kirchhoff stress tensor and its physical counterpart, involves the contravariant Almansi strain tensor. The stress stretch relationship comprises of a component associated with volume constrained distortion and a hydrostatic pressure which results in volumetric dilation. An important property of this constitutive relationship is that the hydrostatic pressure component of the stress vector which is associated with volumetric dilation will have no shear component on any surface in any configuration. This same property is not true for a neo-Hookean Green’s strain–second Piola–Kirchhoff stress tensor formulation.     

Bounds and perturbation series for incompressible elastic composites with transverse isotropic symmetry

The effective elastic behavior of a transversely isotropic composite made from two incompressible elastic materials is examined. The set of all effective elasticity tensors for transversely isotropic finite rank laminar microstructures is described. The extremal property of this class of microstructures is used to derive a new more precise characterization of the set of effective shear moduli.The perturbation series for the effective elasticity tensor is considered. An explicit formula for the second order perturbation tensor is derived. We describe precisely the set of tensors that correspond to all second order perturbations consistent with transverse isotropy. We apply analytic methods [cf. 27] to show that all second order perturbation tensors are realized by finite rank laminar microstructures.Supported by NSF through Grant DMS-3907658.     

关于压电介质的守恒方程和路径无关积分

借鉴分析动力学中的Ｊａｃｏｂｉ积分和循环积分概念，以及电磁场理论中的能量矩概念，导出了压电介质在静态场中的守恒方程形式，由这些守恒方程即可得在位错，断裂力学和其他缺陷理论中应用广泛的路径无关积分。     

Plane strain slip line theory for anisotropic rigid/plastic materials

A slip line theory governing states of incipient plane flow at the yield point load is developed for anisotropic rigid/plastic materials which exhibit a reduced yield criterion, governing states of plane flow, that depends only on the deviatoric parts of the in-plane stress tensor. It is shown that every homogeneous, incompressible material which complies with the principle of maximum plastic work, but which is of otherwise arbitrary anisotropy, is of this class. The stress equilibrium requirements are seen to take a remarkably simple form expressing the constancy of the quantities mean in-plane normal stress plus or minus arc length around the governing yield contour in a Mohr stress plane along members of the two slip line families. Further, this generalization of the Hencky equations is valid for every material of the considered class. Some special features of yield contours containing corners and flat segments are discussed, and velocity equations are given for materials complying with the maximum work inequality. The theory is applied to obtain the solution for indentation of an arbitrarily anisotropic half-space with a flat-ended, rigid, frictionless punch. A simple, universal formula, involving only geometrical dimensions of the governing yield contour, is derived for the yield point indentation pressure.     

具有小密度差的两层流体中运动点源的二阶内波解

在具有自由面的两层流体中，运动点源生成的Kelvin船波存在两种模式，即表面波模式和内波模式。当上、下层流体密度比趋于1时，由内波模式计算的界面波幅趋于无穷大，这与实验事实相违背。为克服此困难，在自由面和界面作小波幅运动的假设，引入一个小密度差参数。研究了运动点源在无粘、不可压且具有小密度差的两层有限深流体中生成的高阶波动。首先利用摄动方法推导了各阶小参数满足的边值问题；其次，给出了小密度差情形下的可解性条件。证明了在密度比趋于1的极限情形，不存在导致界面波幅无穷大的内波模式；最后，利用Phillips的非线性共振相互作用理论，构造了具有自由面的两层有限深流体中Kelvin船波系的二阶一致有效波动解，并证明了该解在深水情形下退化为Newman关于均匀流体中自由面的二阶波动解。     

Finite amplitude acceleration waves in incompressible BKZ viscoelastic fluids— The propagation condition

A general propagation condition is derived which permits the calculation of the speed of propagation of a second-order acceleration wave passing through a particular non-linear incompressible viscoelastic fluid. The viscoelastic fluid is taken to obey the Bernstein, Kearsley and Zapas single integral constitutive model. The analysis is valid for arbitrary finite amplitude waves propagating through a medium undergoing an arbitrary large deformation. Three examples, rest history, steady simple shearing flow and steady simple extensional flow, are given to demonstrate the utility of the propagation condition.     

Two exact micromechanics-based nonlocal constitutive equations for random linear elastic composite materials

A Hashin-Shtrikman-Willis variational principle is employed to derive two exact micromechanics-based nonlocal constitutive equations relating ensemble averages of stress and strain for two-phase, and also many types of multi-phase, random linear elastic composite materials. By exact is meant that the constitutive equations employ the complete spatially-varying ensemble-average strain field, not gradient approximations to it as were employed in the previous, related work of Drugan and Willis (J. Mech. Phys. Solids 44 (1996) 497) and Drugan (J. Mech. Phys. Solids 48 (2000) 1359) (and in other, more phenomenological works). Thus, the nonlocal constitutive equations obtained here are valid for arbitrary ensemble-average strain fields, not restricted to slowly-varying ones as is the case for gradient-approximate nonlocal constitutive equations. One approach presented shows how to solve the integral equations arising from the variational principle directly and exactly, for a special, physically reasonable choice of the homogeneous comparison material. The resulting nonlocal constitutive equation is applicable to composites of arbitrary anisotropy, and arbitrary phase contrast and volume fraction. One exact nonlocal constitutive equation derived using this approach is valid for two-phase composites having any statistically uniform distribution of phases, accounting for up through two-point statistics and arbitrary phase shape. It is also shown that the same approach can be used to derive exact nonlocal constitutive equations for a large class of composites comprised of more than two phases, still permitting arbitrary elastic anisotropy. The second approach presented employs three-dimensional Fourier transforms, resulting in a nonlocal constitutive equation valid for arbitrary choices of the comparison modulus for isotropic composites. This approach is based on use of the general representation of an isotropic fourth-rank tensor function of a vector variable, and its inverse. The exact nonlocal constitutive equations derived from these two approaches are applied to some example cases, directly rationalizing some recently-obtained numerical simulation results and assessing the accuracy of previous results based on gradient-approximate nonlocal constitutive equations.     

Eshelby's problem of non-elliptical inclusions

The Eshelby problem consists in determining the strain field of an infinite linearly elastic homogeneous medium due to a uniform eigenstrain prescribed over a subdomain, called inclusion, of the medium. The salient feature of Eshelby's solution for an ellipsoidal inclusion is that the strain tensor field inside the latter is uniform. This uniformity has the important consequence that the solution to the fundamental problem of determination of the strain field in an infinite linearly elastic homogeneous medium containing an embedded ellipsoidal inhomogeneity and subjected to remote uniform loading can be readily deduced from Eshelby's solution for an ellipsoidal inclusion upon imposing appropriate uniform eigenstrains. Based on this result, most of the existing micromechanics schemes dedicated to estimating the effective properties of inhomogeneous materials have been nevertheless applied to a number of materials of practical interest where inhomogeneities are in reality non-ellipsoidal. Aiming to examine the validity of the ellipsoidal approximation of inhomogeneities underlying various micromechanics schemes, we first derive a new boundary integral expression for calculating Eshelby's tensor field (ETF) in the context of two-dimensional isotropic elasticity. The simple and compact structure of the new boundary integral expression leads us to obtain the explicit expressions of ETF and its average for a wide variety of non-elliptical inclusions including arbitrary polygonal ones and those characterized by the finite Laurent series. In light of these new analytical results, we show that: (i) the elliptical approximation to the average of ETF is valid for a convex non-elliptical inclusion but becomes inacceptable for a non-convex non-elliptical inclusion; (ii) in general, the Eshelby tensor field inside a non-elliptical inclusion is quite non-uniform and cannot be replaced by its average; (iii) the substitution of the generalized Eshelby tensor involved in various micromechanics schemes by the average Eshelby tensor for non-elliptical inhomogeneities is in general inadmissible.