On constitutive modelling of porous neo-Hookean composites

In this paper a hyperelastic constitutive model is developed for neo-Hookean composites with aligned continuous cylindrical pores in the finite elasticity regime. Although the matrix is incompressible, the composite itself is compressible because of the existence of voids. For this compressible transversely isotropic material, the deformation gradient can be decomposed multiplicatively into three parts: an isochoric uniaxial deformation along the preferred direction of the material (which is identical to the direction of the cylindrical pores here); an equi-biaxial deformation on the transverse plane (the plane perpendicular to the preferred direction); and subsequent shear deformation (which includes “along-fibre” shear and transverse shear). Compared to the multiplicative decomposition used in our previous model for incompressible fibre reinforced composites [Guo, Z., Peng, X.Q., Moran, B., 2006, A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus. J. Mech. Phys. Solids 54(9), 1952–1971], the equi-biaxial deformation is introduced to achieve the desired volume change. To estimate the strain energy function for this composite, a cylindrical composite element model is developed. Analytically exact strain distributions in the composite element model are derived for the isochoric uniaxial deformation along the preferred direction, the equi-biaxial deformation on the transverse plane, as well as the “along-fibre” shear deformation. The effective shear modulus from conventional composites theory based on the infinitesimal strain linear elasticity is extended to the present finite deformation regime to estimate the strain energy related to the transverse shear deformation, which leads to an explicit formula for the strain energy function of the composite under a general finite deformation state.     

Conditions for Strong Ellipticity of Anisotropic Elastic Materials

We consider the problem of finding the transversely isotropic elasticity tensor closest to a given elasticity tensor with respect to the Frobenius norm. A similar problem was considered by other authors and solved analytically assuming a fixed orientation of the natural coordinate system of the transversely isotropic tensor. In this paper we formulate a method for finding the optimal orientation of the coordinate system—the one that produces the shortest distance. The optimization problem reduces to finding the absolute maximum of a homogeneous eighth-degree polynomial on a two-dimensional sphere. This formulation allows us a convenient visualization of local extrema, and enables us to find the closest transversely isotropic tensor numerically.      

Propagation of wave at the boundary surface of transversely isotropic thermoelastic material with voids and isotropic elastic half-space

The purpose of this research is to study the effect of voids on the surface wave propagation in a layer of a transversely isotropic thermoelastic material with voids lying over an isotropic elastic half-space. The frequency equation is derived after developing a mathematical model for welded and smooth contact boundary conditions. The dispersion curves giving the phase velocity and attenuation coefficient via wave number are plotted graphically to depict the effects of voids and anisotropy for welded contact boundary conditions. The specific loss and amplitudes of the volume fraction field, the normal stress, and the temperature change for welded contact are obtained and shown graphically for a particular model to depict the voids and anisotropy effects. Some special cases are also deduced from the present investigation.     

Theory of indentation on multiferroic composite materials

This article presents a general theory on indentation over a multiferroic composite half-space. The material is transversely isotropic and magneto-electro-elastic with its axis of symmetry normal to the surface of the half-space. Based on the corresponding half-space Green’s functions to point sources applied on the surface, explicit expressions for the generalized pressure vs. indentation depth are derived for the first time for the three common indenters (flat-ended, conical, and spherical punches). The important multiphase coupling issue is discussed in detail, with the weak and strong coupling being correctly revisited. The derived analytical solutions of indentation will not only serve as benchmarks for future numerical studies of multiphase composites, but also have important applications to experimental test and characterization of multiphase materials, in particular, of multiferroic properties.     

n the fourth order tensor valued function of the stress in return map algorithm

针对各向同性材料,基于一组相互正交的基张量,建立了一套有 效的相关运算方法. 基张量中的两个分别是归一化的二阶单位张量和偏应力张量,另一个则 使用应力的各向同性二阶张量值函数经过归一化构造所得,三者共主轴. 根据张量函数表示 定理,本构方程和返回映射算法中所涉及到的应力的二阶、四阶张量值函数及其逆都由这组 基所表示. 推演结果表明：这些张量之间的运算,表现为对应系数矩阵之间的简单 关系. 其中,四阶张量求逆归结为对应的3\times3系数矩阵求逆,它对二阶张量的变换 则表现为该矩阵对3times 1列阵的变换. 最后,对这些变换关系应用于返回映 射算法的迭代格式进行了相关讨论.     

An Area Modulus of Elasticity: Definition and Properties

A new modulus of elasticity is defined to be the ratio of an equibiaxial stress to the relative area change in the planes in which the stress acts. This area modulus of elasticity is intermediate in properties between Young's modulus and the bulk modulus. Expressions for the area modulus are computed in isotropic elasticity. A simple, convenient expression for the compliance tensor of transverse isotropy is found in terms of, amongst others, the longitudinal (axial) area modulus and this leads to a new, concise condition for positive definiteness of the compliance tensor. The limits of incompressibility, inextensibility and constant area are briefly considered. This revised version was published online in August 2006 with corrections to the Cover Date.     

On well-posedness of the dynamic problem for an anisotropic fluid-saturated solid

It is known that a high degree of anisotropy in the constitutive behaviour of a solid may result in the loss of hyperbolicity of the dynamic equations in the form of either complex-conjugate or purely imaginary characteristic wave speeds (flutter ill-posedness and shear band formation, respectively). In the present paper we investigate the characteristic wave speeds in the dynamic problem for a transversely isotropic fluid-saturated porous solid. Three cases are considered: a dry solid and a saturated solid under locally undrained and drained conditions. It is shown that, for given constitutive parameters of the solid skeleton, the dynamic problem for a drained solid may become ill-posed due to the flutter-type loss of hyperbolicity, while the dynamic equations for a dry and an undrained solids remain hyperbolic. For a given solid skeleton, the characteristic wave speeds are strongly influenced by the pore fluid compressibility which, in turn, is extremely sensitive to the presence of a small amount of free gas.     

On the Bending, Vibration and Stability of Laminated Rectangular Plates with Transversely Isotropic Layers

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｓｔａｔｉｃｓａｎｄｄｙｎａｍｉｃｓｏｆｒｅｃｔａｎｇｕｌａｒｐｌａｔｅｓｈａｖｅａｌｗａｙｓｂｅｅｎｏｆｐａｒｔｉｃｕｌａｒｒｅｓｅａｒｃｈｉｎｔｅｒｅｓｔｉｎｓｏｌｉｄｍｅｃｈａｎｉｃｓ.Ｉｎｐａｒｔｉｃｕｌａｒ,ｔｈｅｅｍｅｒｇｅｎｃｅｏｆｎｅｗｍａｔｅｒｉａｌｓａｎｄｎｅｗｓｔｒｕｃｔｕｒｅｓｈａｓｂｒｏｕｇｈｔｏｕｔｅｎｄｌｅｓｓｖｉｔａｌｆｏｒｃｅｉｎｔｏｔｈｅｓｔｕｄｙ .Ｉｔｉｓｗｅｌｌ＿ｋｎｏｗｎｔｈａｔｔｈｅｃｌａｓｓｉｃａｌｐｌａｔｅｔｈｅｏｒｙｂ…     

Determination of a Homogeneous Strictly Pseudoconvex Surface from the Coefficients of Its Normal Equation

Mathematical Notes -     

Axisymmetric wave propagation in a layered transversely isotropic medium

Summary In this paper, the method of numerical integration along bicharacteristics is generalized to the case of layered transversely isotropic medium for analysing the axisymmetric stress wave propagation. The stability of the present scheme is studied. The advantages and limitations of the method are discussed. Received 12 June 1996; accepted for publication 6 May 1997