Stability,bifurcation and post-critical behavior of a homogeneously deformed incompressible isotropic elastic parallelepiped subject to dead-load surface tractions

We study the equilibrium homogeneous deformations of a homogeneous parallelepiped made of an arbitrary incompressible, isotropic elastic material and subject to a distribution of dead-load surface tractions corresponding to an equibiaxial tensile stress state accompanied by an orthogonal uniaxial compression of the same amount. We show that only two classes of homogeneous equilibrium solutions are possible, namely symmetric deformations, characterized by two equal principal stretches, and asymmetric deformations, with all different principal stretches. Following the classical energy-stability criterion, we then find necessary and sufficient conditions for both symmetric and asymmetric equilibrium deformations to be weak relative minimizers of the total potential energy. Finally, we analyze the mechanical response of a parallelepiped made of an incompressible Mooney–Rivlin material in a monotonic dead loading process starting from the unloaded state. As a major result, we model the actual occurrence of a bifurcation from a primary branch of locally stable symmetric deformations to a secondary, post-critical branch of locally stable asymmetric solutions.     

Shear-band bifurcation of an isotropic compressible hyperelastic solid

This paper concerns shear-band bifurcations from the homogeneous finite plane deformation of an isotropic compressible hyperelastic solid. The governing equations for the incremental plane deformation superposed on the initial finite deformation are derived and then the equilibrium equations in terms of incremental displacements are classified into the elliptic type, parabolic type, etc. From this classification follows a restriction which should be placed on the strain-energy function in order that the equilibrium equations may be either elliptic or parabolic for all principal stretches. For the hyperelastic solid complying with this restriction, the condition for the shear-band bifurcation is obtained. Finally, the incremental displacement field of an infinite series of shear bands in a slab sandwiched between slippery rigid layers is established.     

Constitutive branching analysis of cylindrical bodies under in-plane equibiaxial dead-load tractions

Finite homogeneous deformations of hyperelastic cylindrical bodies subjected to in-plane equibiaxial dead-load tractions are analyzed. Four basic equilibrium problems are formulated considering incompressible and compressible isotropic bodies under plane stress and plane deformation condition. Depending on the form of the stored energy function, these plane problems, in addition to the obvious symmetric solutions, may admit asymmetric solutions. In other words, the body may assume an equilibrium configuration characterized by two unequal in-plane principal stretches corresponding to equal external forces. In this paper, a mathematical condition, in terms of the principal invariants, governing the global development of the asymmetric deformation branches is obtained and examined in detail with regard to different choices of the stored energy function. Moreover, explicit expressions for evaluating critical loads and bifurcation points are derived. With reference to neo-Hookean, Mooney-Rivlin and Ogden-Ball materials, a broad numerical analysis is performed and the qualitatively more interesting asymmetric equilibrium branches are shown. Finally, using the energy criterion, a number of considerations are put forward about the stability of the computed solutions.     

The Attainable Region of strain-invariant space for elastic materials

Within the theory of isothermal isotropic non-linear elasticity, the selection of the appropriate form for the strain energy function W in terms of the strain invariants is still an issue. The purpose of this paper is to introduce ideas and techniques which it is hoped will contribute to the task of finding an appropriate form for the strain energy. Three principal ideas are developed in this paper. Firstly, not all of invariant-space corresponds to real deformations. Constitutive equations only need to match real behaviour over a restricted part of invariant space, called the Attainable Region, bounded by states of deformation corresponding to uniaxial and equi-biaxial extension. Secondly, examples are given of how to exploit the fact that the Attainable Region is restricted. Mapping a deformation onto this region allows visualization of how close the deformation is to the well-understood uniaxial, equi-biaxial and simple shear deformations, and how this varies in space or time. Thirdly, acceptable strain invariants do not have to be obviously symmetric functions of the principal stretches. The ordered principal stretches are themselves invariants, and explicit unique algebraic expressions can be given through which the greatest, middle or least stretch can be calculated in terms of the usual invariants. Thus invariants can be chosen which are apparently non-symmetric functions of the ordered stretches.     

Volume changes associated with the deformation of rubber-like solids

Volumechanges accompanying the deformation of rubber-like solids are analysed on the basis of isotropic elasticity theory. In particular, a simple, but general, result relating the volume and the stretch in simple tension is obtained. This is achieved by the introduction of certain modified principal stretches which allow the dependence of the strain energy on the isochoric and the dilatational parts of the deformation to be considered separately. The fact that volume changes in rubber-like solids are typically of order 0.01 % is used to linearize the stress-deformation relations in the dilatation. This enables the dilatation to be given explicitly as a function of the stretch in simple tension. Specific results are obtained for certain classes of constitutive law and good agreement with the experimental data for simple tension is demonstrated.Results for equi-biaxial tension and pure shear are also given in anticipation of further experimental data becoming available. The need for volume-change data for a wide variety of types of strains is emphasized.     

Physical invariants for nonlinear orthotropic solids

Seven invariants, with immediate physical interpretation, are proposed for the strain energy function of nonlinear orthotropic elastic solids. Three of the seven invariants are the principal stretch ratios and the other four are squares of the dot product between the two preferred directions and two principal directions of the right stretch tensor. A strain energy function, expressed in terms of these invariants, has a symmetrical property almost similar to that of an isotropic elastic solid written in terms of principal stretches. Ground state and stress–strain relations are given. Using principal axes techniques, the formulation is applied, with mathematical simplicity, to several types of deformations. In simple shear, a necessary and sufficient condition is given for Poynting relation and two novel deformation-dependent universal relations are formulated. Using series expansions and the symmetrical property, the proposed general strain energy function is refined to a particular general form. A type of strain energy function, where the ground state constants are written explicitly, is proposed. Some advantages of this type of function are indicated. An experimental advantage is demonstrated by showing a simple triaxial test can vary a single invariant while keeping the remaining invariants fixed.     

Extremal paths of plastic work and deformation

Optimal paths of deformation between homogeneous states of finite strain are identified on the basis that the work expended should be an absolute minimum. The materials considered are classically rigid/plastic, isotropic and nonhardening, with any convex yield surface. A related minimum principle is derived for a class of functionals of the strain history alone, without reference to work or to material properties. Plane isochoric deformations are subsequently treated in particular detail, supported by simple parametric representations of the competing paths and associated kinematic data.     

Some New Implications of Certain Constitutive Inequalities in Plane Isotropic Elasticity

In plane isotropic elasticity a strengthened form of the Ordered–Forces inequality is shown to imply that the restriction of the strain-energy function to the class of deformation gradients which share the same average of the principal stretches is bounded from below by the strain energy corresponding to the conformal deformations in this class. For boundary conditions of place, this property (together with a certain version of the Pressure–Compression inequality) is then used (i) to show that the plane radial conformal deformations are stable with respect to all radial variations of class C ¹ and (ii) to obtain explicit lower bounds for the total energy associated with arbitrary plane radial deformations. For the same type of boundary conditions and together with a different version of the Pressure–Compression inequality, an analogous property in plane isotropic elasticity (established in [3] under the assumption that the material satisfies a strengthened form of the Baker–Ericksen inequality and according to which the restriction of the strain-energy function to the class of deformation gradients which share the same determinant is bounded from below by the strain energy corresponding to the conformal deformations in that class) is used (i) to show that the plane radial conformal deformations are stable with respect to all variations of class C ¹ and (ii) to obtain explicit lower bounds for the total energy associated with any plane deformation.     

On the failure of ellipticity of the equations for finite elastostatic plane strain

Summary In this paper we establish necessary and sufficient conditions, in terms of the local principal stretches, for ordinary and strong ellipticity of the equations governing finite plane equilibrium deformations of a compressible hyperelastic solid. The material under consideration is assumed to be homogeneous and isotropic, but its strain-energy density is otherwise unrestricted. We also determine the directions of the characteristic curves appropriate to plane elastostatic deformations that are accompanied by a failure of ellipticity.The results communicated in this paper were obtained in the course of an investigation supported by Contract N00014-75-C-0196 with the Office of Naval Research in Washington, D.C.     

Foam rheology: Relation between extensional and shear deformations in high gas fraction foams

A two-dimensional hexagonal foam cell model is used to derive analytic expressions for the bulk stress tensor and foam microstructure for any small homogeneous deformation. We show that calculations done for deformations where the principal axes of stress and strain coincide, such as in extension, are sufficient to provide all information about shear deformation. The stresses and foam structure for any given strain and initial cell orientation in shear bears a unique relation to a different strain and orientation in extension. Such a mapping is obtained using the assumption that the principal axes of strain and stress corotate with each other. This in turn implies that high gas fraction foams follow the Lodge-Meissner relation, i.e. the ratio of the normal-stress difference to the shear stress equals the shear strain. The spatially periodic structure of foam along with the fact that the cell centers move affinely with the bulk, makes the above assumption a justifiable one.     

Stress-Induced Phase Transformations of an Isotropic Incompressible Elastic Body Subject to a Triaxial Stress State

The present paper concerns the stable multiphase isochoric deformations for an isotropic elastic body subject to a surface traction of uniform Piola stress with two equal principal forces which are opposite to the third. To model the occurrence of such deformations, we consider a strain energy density function which depends on the first principal invariant of deformation through a non-convex function and which has an added linear dependence on the second invariant. We establish existence conditions for equilibrium multiphase deformations which give restrictions on the morphology of the connecting phases as well as on the orientation of the flat interfaces between the phases. Finally, by considering a special, but representative, form for the non-convex strain energy function, we show that there exists a “critical” value of the external load which allows for the emergence of stable coexistent deformation fields.     

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.     

Green's function for incremental nonlinear elasticity: shear bands and boundary integral formulation

An elastic, incompressible, infinite body is considered subject to plane and homogeneous deformation. At a certain value of the loading, when the material is still in the elliptic range, an incremental concentrated line load is considered acting at an arbitrary location in the body and extending orthogonally to the plane of deformation. This plane strain problem is solved, so that a Green's function for incremental, nonlinear elastic deformation is obtained. This is used in two different ways: to quantify the decay rate of self-equilibrated loads in a homogeneously stretched elastic solid; and to give a boundary element formulation for incremental deformations superimposed upon a given homogeneous strain. The former result provides a perturbative approach to shear bands, which are shown to develop in the elliptic range, induced by self-equilibrated perturbations. The latter result lays the foundations for a rigorous approach to boundary element techniques in finite strain elasticity.     

A semi-inverse shape optimization problem in linear anti-plane shear

The design of a semi-infinite fillet for efficient stress transmission is considered. The problem is treated within the context of anti-plane shear deformations of a homogeneous, isotropic, linearly elastic solid. Under a remote state of simple shear, it is desired to determine the shape of the traction-free lateral boundaries of a symmetric plane domain so that the shear stress distribution on the finite end is as uniform as possible. A semi-inverse approach for a particular class of semi-infinite profiles is used to examine this issue.     

New universal relations for nonlinear isotropic elastic materials

A nonlinear isotropic elastic block is subjected to a homogeneous deformation consisting of simple shear superposed on triaxial extension. Two new relations are established for this deformation which are valid for all nonlinear elastic isotropic materials, and hence are universal relations. The first is a relation between the stretch ratios in the plane of shear and the amount of shear when the deformation is supported only by shear tractions. The second relation is established for a thin-walled cylinder under combined extension, inflation and torsion. Each material element of the cylinder undergoes the same local homogeneous deformation of shear superposed on triaxial extension. The properties of this deformation are used to establish a relation between pressure, twisting moment, angle of twist and current dimensions when no axial force is applied to the cylinder. It is shown that these relations also apply for a mixture of a nonlinear isotropic solid and a fluid.     

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.     

Lateral buckling of beams with shear deformations – A hyperelastic formulation

The equilibrium and buckling equations are derived for the lateral buckling of a prismatic straight beam. A consistent finite strain constitutive law is used, which is based on a hyperelastic model for an isotropic material. The kinematics of the cross-sectional deformations are based on a Timoshenko type beam displacement of the cross-sectional plane using Euler angles and two shear finite rotations coupled with warping taken normal to the displaced plane. Also derived are the second order approximations to the displacements, curvatures, twist and internal actions. The constitutive relationships for the internal actions reveal new coupling terms between the bending moments, torsion and bimoment, which are functions of the cross-sectional warping and shear deformations. New Wagner type nonlinear torsion terms are derived which are functions of the warping of the cross-sectional plane, and are coupled to the twisting and shear deformations of the cross-section. Solutions are determined for the lateral buckling of a prismatic monosymmetric beam under pure bending and the flexural–torsional buckling under axial compression. For the flexural–torsional buckling problem it is found that the Euler type column buckling formula is consistent with Haringx’s column buckling formula while the torsional buckling formula is different to conventional equations. The second variation of the total potential is also derived. The effects of shear deformations are explored by examining the non-dimensional lateral buckling equation for a simply supported beam.     

Simple Shearing of Incompressible and Slightly Compressible Isotropic Nonlinearly Elastic Materials

The classical problem of simple shear in nonlinear elasticity has played an important role as a basic pilot problem involving a homogeneous deformation that is rich enough to illustrate several key features of the nonlinear theory, most notably the presence of normal stress effects. Here our focus is on certain ambiguities in the formulation of simple shear arising from the determination of the arbitrary hydrostatic pressure term in the normal stresses for the case of an incompressible isotropic hyperelastic material. A new formulation in terms of the principal stretches is given. An alternative approach to the determination of the hydrostatic pressure is proposed here: it will be required that the stress distribution for a perfectly incompressible material be the same as that for a slightly compressible counterpart. The form of slight compressibility adopted here is that usually assumed in the finite element simulation of rubbers. For the particular case of a neo-Hookean material, the different stress distributions are compared and contrasted.     

剪切变形对FGM圆板轴对称屈曲的影响

基于一阶剪切变形板理论，推导了功能梯度材料圆形板在边界面内均布压力作用下的轴对称屈曲方程。在推导过程中，忽略了前屈曲耦合变形。利用一阶板理论与经典板理论屈曲方程之间在数学形式上的相似性，得到了一阶板理论下功能梯度材料圆板与经典板理论下均匀圆板临界屈曲载荷之间的解析关系。利用这个解析关系，可以直接从已有的较为简单的经典理论的结果，获得一阶板理论下功能梯度材料板的临界屈曲载荷。