From Non-Linear Elasticity to Linearized Theory: Examples Defying Intuition

Material frame indifference implies that the solution in non-linear elasticity theory for a connected body rigidly rotated at its border is a rigid, stress-free, deformation. If the same problem is considered within linear elasticity theory, considered as an approximation to the true elastic situation, one should expect that if the angle of rotation is small, the body still undergoes a rigid deformation while the corresponding stress, though not zero, remains consistently small. Here, we show that this is true, in general, only for homogeneous bodies. Counterexamples of inhomogeneous bodies are presented for which, whatever small the angle of rotation is, the linear elastic solution is by no means a rigid rotation (in a particular case it is an “explosion”) while the stress may even become infinite. If the same examples are re-interpreted as problems in an elasticity theory based upon genuinely linear constitutive relations which retain their validity also for finite deformations, it is shown that they would deliver constraint reaction forces that are not in equilibrium in the actual, deformed, state. This furnishes another characterization of the impossibility of an exact linear constitutive theory for elastic solids with zero residual stress.      

Finite deformations of fibre-reinforced elastic solids with fibre bending stiffness

In the conventional theory of finite deformations of fibre-reinforced elastic solids it is assumed that the strain-energy is an isotropic invariant function of the deformation and a unit vector A that defines the fibre direction and is convected with the material. This leads to a constitutive equation that involves no natural length. To incorporate fibre bending stiffness into a continuum theory, we make the more general assumption that the strain-energy depends on deformation, fibre direction, and the gradients of the fibre direction in the deformed configuration. The resulting extended theory requires, in general, a non-symmetric stress and the couple-stress. The constitutive equations for stress and couple-stress are formulated in a general way, and specialized to the case in which dependence on the fibre direction gradients is restricted to dependence on their directional derivatives in the fibre direction. This is further specialized to the case of plane strain, and finite pure bending of a thick plate is solved as an example. We also formulate and develop the linearized theory in which the stress and couple-stress are linear functions of the first and second spacial derivatives of the displacement. In this case for the symmetric part of the stress we recover the standard equations of transversely isotropic linear elasticity, with five elastic moduli, and find that, in the most general case, a further seven moduli are required to characterize the couple-stress.     

On Isotropic Invariants of the Elasticity Tensor

Isotropic invariants of the elasticity tensor always yield the same values no matter what coordinate system is concerned and therefore they characterize the linear elasticity of a solid material intrinsically. There exists a finite set of invariants of the elasticity tensor such that each invariant of the elasticity tensor can be expressed as a single-valued function of this set. Such a set, called a basis of invariants of the elasticity tensor, can be used to realize a parametrization of the manifold of orbits of elastic moduli, i.e. to distinguish different kinds of linear elastic materials. Seeking such a basis is an old problem in theory of invariants and seems to have been unsuccessful until now. In this paper, by means of the unique spectral decomposition of the elasticity tensor every invariant of the elasticity tensor is shown to be a joint invariant of the eigenprojections of the elasticity tensor, and then by utilizing some properties of the eigenprojections a basis for each case concerning the multiplicity of the eigenvalues of the elasticity tensor is presented in terms of joint invariants of the eigenprojections. In addition to the foregoing properties, the presented invariants may also be used to form invariant criteria for identification of elastic symmetry axes.     

An atomistic and non-classical continuum field theoretic perspective of elastic interactions between defects (force dipoles) of various symmetries and application to graphene

Force multipoles are employed to represent various types of defects and physical phenomena in solids: point defects (interstitials, vacancies), surface steps and islands, proteins on biological membranes, inclusions, extended defects, and biological cell interactions among others. In the present work, we (i) as a prototype simple test case, conduct quantum mechanical calculations for mechanics of defects in graphene sheet and in parallel, (ii) formulate an enriched continuum elasticity theory of force dipoles of various anisotropies incorporating up to second gradients of strain fields (thus accounting for nonlocal dispersive effects) instead of the usual dispersion-less classical elasticity formulation that depends on just the strain (c.f. Peyla, P., Misbah, C., 2003. Elastic interaction between defects in thin and 2-D films. Eur. Phys. J. B. 33, 233-247). The fundamental Green's function is derived for the governing equations of second gradient elasticity and the elastic self and interaction energies between force dipoles are formulated for both the two-dimensional thin film and the three-dimensional case. While our continuum results asymptotically yield the same interaction energy law as Peyla and Misbah for large defect separations (∼1/rⁿ for defects with n-fold symmetry), the near-field interactions are qualitatively far more complex and free of singularities. Certain qualitative behavior of defect mechanics predicted by atomistic calculations are well captured by our enriched continuum models in contrast to classical elasticity calculations. For example, consistent with our atomistic calculations of defects in isotropic graphene, even two dilation centers show a finite interaction (as opposed to classical elasticity that predicts zero interaction). We explicitly find the physically consistent result that the self-energy of a defect is equivalent to half the interaction energy between two identical defects when they “merge” into each other. The atomistic, classical elastic and the enriched continuum predictions are thoroughly compared for two types of defects in graphene: Stone-Wales and divacancy.     

Nonlinear guided bulk waves in heterogeneous elastic structural elements

Guided nonlinear bulk strain waves are considered in heterogeneous composite materials, containing macro- or micro(nano)inclusions as the main sources of nonuniformity. Theory is based on the combination of the finite deformation and the 5-constant (Murnaghan) nonlinear elasticity, that leads under several assumptions to the only nonlinear doubly dispersive equation with variable coefficients (DDE) for a component of longitudinal strain in rod and in shell. Numerical simulations performed and experimental data are discussed, and lead to conclusions concerning dramatic influence of heterogeneity in dynamic behaviour of strain solitons in solids. We have shown how different smooth variations in elasticity may lead to either amplification or decay of strain solitons in both polymer rod and shell. These date can be useful in the NDT problems, and in problems of the solids integrity under intense elastic pulse loading.     

Some problems of linear elasticity for cylinders in micropolar orthotropic material

Some special problems for axisymmetric solids made of linearly elastic orthotropic micropolar material with central symmetry are dealt with. The first one is a hollow circular cylinder of unlimited length, subjected to internal and external uniform pressure. The second one is a hollow or solid circular cylinder of finite length, subjected to a relative rotation of the bases about its axis. In both cases, one of the axes of elastic symmetry is parallel to the cylinder axis; the other two are arbitrarily oriented in the plane of any cross-section of the solid. The elastic properties are invariant along the cylinder axis. It is shown that the two problems are governed by formally similar sets of ordinary differential equations in the kinematic fields (in-plane displacements and microrotations). In the general case, numerical solutions are derived. The solution for the cylinder subjected to radial pressure does not significantly differ from that obtained in classical elasticity, at least in terms of radial and hoop force stresses. In the case of a cylinder subjected to torsion the difference between the micropolar and the classical solutions is more pronounced. The torque induces twisting couple stresses about the cylinder axis of variable sign. Finally, size effects in terms of torsional inertia are pointed out.     

Theory of nonlocal asymmetric elastic solids

In this paper, a nonlinear theory of nonlocal asymmetric, elastic solids is developed on the basis of basic theories of nonlocal continuum fieM theory and nonlinear continuum mechanics. It perfects and expands the nonlocal elastic fiteld theory developed by Eringen and others. The linear theory of nonlocal asymmetric elasticity developed in [1] expands to the finite deformation, We show that there is the nonlocal body moment in the nonlocal elastic solids. The noniocal body moment causes the stress asymmetric and itself is caused by the covalent bond formed by the reaction between atoms. The theory developed in this paper is applied to explain reasonably that curves of dispersion relation of one-dimensional plane longitudinal waves are not similar with those of transverse waves.     

Waves in Constrained Linear Elastic Materials

This paper deals with the propagation of acceleration waves in constrained linear elastic materials, within the framework of the so-called linearized finite theory of elasticity, as defined by Hoger and Johnson in [12, 13]. In this theory, the constitutive equations are obtained by linearization of the corresponding finite constitutive equations with respect to the displacement gradient and significantly differ from those of the classical linear theory of elasticity. First, following the same procedure used for the constitutive equations, the amplitude condition for a general constraint is obtained. Explicit results for the amplitude condition for incompressible and inextensible materials are also given and compared with those of the classical linear theory of elasticity. In particular, it is shown that for the constraint of incompressibility the classical linear elasticity provides an amplitude condition that, coincidently, is correct, while for the constraint of inextensibility the disagreement is first order in the displacement gradient. Then, the propagation condition for the constraints of incompressibility and inextensibility is studied. For incompressible materials the propagation condition is solved and explicit values for the squares of the speeds of propagation are obtained. For inextensible materials the propagation condition is solved for plane acceleration waves propagating into a homogeneously strained material. For both constraints, it is shown that the squares of the speeds of propagation depend by terms that are first order in the displacement gradient, while in classical linear elasticity they are constant. This revised version was published online in July 2006 with corrections to the Cover Date.     

Bounds on the Effective Anisotropic Elastic Constants

Hill [12] showed that it was possible to construct bounds on the effective isotropic elastic coefficients of a material with triclinic or greater symmetry. Hill noted that the triclinic symmetry coefficients appearing in the bounds could be specialized to those of a greater symmetry, yielding the effective isotropic elastic coefficients for a material with any elastic symmetry. It is shown here that it is possible to construct bounds on the effective elastic constants of a material with any anisotropic elastic symmetry in terms of triclinic symmetry elastic coefficients. Similarly, it is then possible to specialize the triclinic symmetry coefficients appearing in the bounds to those of a greater symmetry. Specific bounds are given for the effective elastic coefficients of cubic, hexagonal, tetragonal and trigonal symmetries in terms of the elastic coefficients of triclinic symmetry. These results are obtained by combining the approach of Hill [12] with a representation of the stress-strain relations due, in principle, to Kelvin [25,26] but recast in the structure of contemporary linear algebra. This revised version was published online in August 2006 with corrections to the Cover Date.     

Non-linear theories of beams and plates accounting for moderate rotations and material length scales

The primary objective of this paper is to formulate the governing equations of shear deformable beams and plates that account for moderate rotations and microstructural material length scales. This is done using two different approaches: (1) a modified von Kármán non-linear theory with modified couple stress model and (2) a gradient elasticity theory of fully constrained finitely deforming hyperelastic cosserat continuum where the directors are constrained to rotate with the body rotation. Such theories would be useful in determining the response of elastic continua, for example, consisting of embedded stiff short fibers or inclusions and that accounts for certain longer range interactions. Unlike a conventional approach based on postulating additional balance laws or ad hoc addition of terms to the strain energy functional, the approaches presented here extend existing ideas to thermodynamically consistent models. Two major ideas introduced are: (1) inclusion of the same order terms in the strain–displacement relations as those in the conventional von Kármán non-linear strains and (2) the use of the polar decomposition theorem as a constraint and a representation for finite rotations in terms of displacement gradients for large deformation beam and plate theories. Classical couple stress theory is recovered for small strains from the ideas expressed in (1) and (2). As a part of this development, an overview of Eringen׳s non-local, Mindlin׳s modified couple stress theory, and the gradient elasticity theory of Srinivasa–Reddy is presented.     

Differential constitutive equations of incompressible media with finite deformations

A method is proposed for constructing a system of constitutive equations of an incompressible medium with nonlinear dissipative properties with finite deformations. A scheme of the mechanical behavior of a material is used, in which the points are connected by horizontally aligned elastic, viscous, plastic, and transmission elements. The properties of each element of the scheme are described with the use of known equations of the nonlinear elasticity theory, the theory of nonlinear viscous fluids, and the theory of plastic flow of the material under conditions of finite deformations of the medium. The system of constitutive equations is closed by equations that express the relation between the deformation rate tensor of the material and the deformation rate tensor of the plastic element. Transmission elements are used to take into account a significant difference between macroscopic deformations of the material and deformations of elements of the medium at the structural level. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 158–170, May–June, 2009.     

The dynamic response of an incompressible non-linearly elastic membrane tube subjected to a dynamic extension

The dynamic response of an isotropic hyperelastic membrane tube, subjected to a dynamic extension at its one end, is studied. In the first part of the paper, an asymptotic expansion technique is used to derive a non-linear membrane theory for finite axially symmetric dynamic deformations of incompressible non-linearly elastic circular cylindrical tubes by starting from the three-dimensional elasticity theory. The equations governing dynamic axially symmetric deformations of the membrane tube are obtained for an arbitrary form of the strain-energy function. In the second part of the paper, finite amplitude wave propagation in an incompressible hyperelastic membrane tube is considered when one end is fixed and the other is subjected to a suddenly applied dynamic extension. A Godunov-type finite volume method is used to solve numerically the corresponding problem. Numerical results are given for the Mooney-Rivlin incompressible material. The question how the present numerical results are related to those obtained in the literature is discussed.     

Initial stresses in elastic solids: Constitutive laws and acoustoelasticity

On the basis of the nonlinear theory of elasticity, the general constitutive equation for an isotropic hyperelastic solid in the presence of initial stress is derived. This derivation involves invariants that couple the deformation with the initial stress and in general, for a compressible material, it requires 10 invariants, reducing to 9 for an incompressible material. Expressions for the Cauchy and nominal stress tensors in a finitely deformed configuration are given along with the elasticity tensor and its specialization to the initially stressed undeformed configuration. The equations governing infinitesimal motions superimposed on a finite deformation are then used to study the combined effects of initial stress and finite deformation on the propagation of homogeneous plane waves in a homogeneously deformed and initially stressed solid of infinite extent. This general framework allows for various different specializations, which make contact with earlier works. In particular, connections with results derived within Biot's classical theory are highlighted. The general results are also specialized to the case of a small initial stress and a small pre-deformation, i.e. to the evaluation of the acoustoelastic effect. Here the formulas derived for the wave speeds cover the case of a second-order elastic solid without initial stress and subject to a uniaxial tension [Hughes and Kelly, Phys. Rev. 92 (1953) 1145] and are consistent with results for an undeformed solid subject to a residual stress [Man and Lu, J. Elasticity 17 (1987) 159]. These formulas provide a basis for acoustic evaluation of the second- and third-order elasticity constants and of the residual stresses. The results are further illustrated in respect of a prototype model of nonlinear elasticity with initial stress, allowing for both finite deformation and nonlinear dependence on the initial stress.     

Research on an orthogonal relationship for orthotropic elasticity

IntroductionAsymplecticsystematicmethodology[1- 3]forelasticitywasestablishedbyZhongWan_xie .Hepresentedcreativelythedualvectorsandthesymplecticorthogonalrelationshipandopenedaworkplatformparalleledtothetraditionalelasticity[4 - 9].AnewdualvectorandanewdualdifferentialmatrixLwerepresentedforasymplecticsystematicmethodologyfortwo_dimensionalelasticityandaneworthogonalrelationshipwasdiscoveredforisotropicplaneproblems[4 ]byLuoJian_hui.Theneworthogonalrelationshipisgeneralizedfororthotropicelas…     

THE INVERSE PROPOSITION ABOUT THE SYMMETRY OF LINEAR CONSTITUTIVE LAW

基于线弹性体本构方程系数矩阵的对称性,提出了其对应的逆命题问题,即若材料本构方程是线性且对称的,能否由此确定物体是完全弹性的?论文通过构造势函数的方法对该问题给出了肯定的回答,从而论证了对于符合线性本构关系的材料,其本构方程的对称性与物体的完全弹性相互蕴含,因而是相互等价的.     

Wave dispersion in gradient elastic solids and structures: A unified treatment

Analytical wave propagation studies in gradient elastic solids and structures are presented. These solids and structures involve an infinite space, a simple axial bar, a Bernoulli–Euler flexural beam and a Kirchhoff flexural plate. In all cases wave dispersion is observed as a result of introducing microstructural effects into the classical elastic material behavior through a simple gradient elasticity theory involving both micro-elastic and micro-inertia characteristics. It is observed that the micro-elastic characteristics are not enough for resulting in realistic dispersion curves and that the micro-inertia characteristics are needed in addition for that purpose for all the cases of solids and structures considered here. It is further observed that there exist similarities between the shear and rotary inertia corrections in the governing equations of motion for bars, beams and plates and the additions of micro-elastic (gradient elastic) and micro-inertia terms in the classical elastic material behavior in order to have wave dispersion in the above structures.     

全应变梯度挠曲电纳米梁有限单元法研究

挠曲电效应是一种存在于所有电介质材料中的特殊的力电耦合效应,本质上是应变梯度与电极化之间的线性耦合。然而,应变梯度会引入位移的高阶偏量,常给挠曲电问题的理论求解带来困难。且已有研究表明应变梯度弹性项会影响纳米结构中的力电耦合响应,但是现有的挠曲电研究大多忽略了应变梯度弹性的影响。因此,本文提出了一种既考虑应变梯度弹性,又考虑挠曲电效应的有效数值方法。基于全应变梯度弹性理论,建立了包含3个独立材料尺度参数的纳米欧拉梁的理论模型和有限元模型,提出了满足C₂弱连续的两节点六自由度单元。基于本文的有限单元法,以简支欧拉梁为例,通过分析讨论挠度、电势和能量效率,得到了挠曲电效应和应变梯度弹性项对梁的力电响应的影响。结果表明,挠曲电效应存在尺寸依赖性,且应变梯度弹性项在纳米电介质结构的挠曲电研究中的影响不可忽略。     

Periodic Homogenization and Material Symmetry in Linear Elasticity

Here homogenization theory is used to establish a connection between the symmetries of a periodic elastic structure associated with the microscopic properties of an elastic material and the material symmetries of the effective, macroscopic elasticity tensor. Previous results of this type exist but here more general symmetries on the microscale are considered. Using an explicit example, we show that it is possible for a material to be fully anisotropic on the microscale and yet the symmetry group on the macroscale can contain elements other than plus or minus the identity. Another example demonstrates that not all material symmetries of the macroscopic elastic tensor are generated by symmetries of the periodic elastic structure.