Canonical moduli and general solution of equations of a two-dimensional static problem of anisotropic elasticity

Equations of a two-dimensional static problem of anisotropic elasticity are brought to a simple form with the use of orthogonal and affine transformations of coordinates and corresponding transformations of mechanical quantities. It is proved that an arbitrary matrix of elasticity moduli containing six independent components can be always converted by a congruent transformation to a matrix with two independent components, which are called the canonical moduli. Depending on the relations between the canonical moduli, the determinant of the matrix of operators of equations in displacements is presented as a product of various quadratic terms. A general presentation of the solution of equations in displacements in the form of a linear combination of the first derivatives of two quasi-harmonic functions satisfying two independent equations is given. A symmetry operator (i.e., a formula of production of new solutions) is found to correspond to each presentation. In a three-dimensional case, the matrix of elasticity moduli with 21 independent components is congruent to a matrix with 12 independent canonical moduli.     

Mathematical Models and Methods of the Mechanics of Stochastic Composites

The basic approaches used in mathematical models and general methods for solution of the equations of the mechanics of stochastic composites are generalized. They can be reduced to the stochastic equations of the theory of elasticity of a structurally inhomogeneous medium, to the equations of the theory of effective elastic moduli, to the equations of the theory of elastic mixtures, or to more general equations of the fourth order. The solution of the stochastic equations of the elastic theory for an arbitrary domain involves substantial mathematical difficulties and may be implemented only rather approximately. The construction of the equations of the theory of effective moduli is associated with the problem on the effective moduli of a stochastically inhomogeneous medium, which can be solved by the perturbation method, by the method of moments, or by the method of conditional moments. The latter method is most appropriate. It permits one to determine the effective moduli in a two-point approximation and nonlinear deformation properties. In the structure of equations, the theory of elastic mixtures is more general than the theory of effective moduli; however, since the state equations have not been strictly substantiated and the constants have not been correctly determined, theoretically or experimentally, this theory cannot be used for systematic designing composite structures. A new model of the nonuniform deformation of composites is more promising. It is constructed by performing strict mathematical transformations and averaging the output stochastic equations, all the constants being determined. In the zero approximation, the equations of the theory of effective moduli follow from this model, and, in the first approximation, fourth-order equations, which are more general than those of the theory of mixtures, follow from it     

Analytical procedure for determining the linear and nonlinear effective properties of the elastic composite cylinder

In this work we consider a cylindrical structure composed of a nonlinear core (inhomogeneity) surrounded by a different nonlinear shell (matrix). We elaborate a technique for determining its linear elastic moduli (second order elastic constants) and the nonlinear elastic moduli, which are called Landau coefficients (third order elastic constants). Firstly, we develop a nonlinear perturbation method which is able to turn the initial nonlinear elastic problem into a couple of linear problems. Then, we prove that only the solution of the first linear problem is necessary to calculate the linear and nonlinear effective properties of the heterogeneous structure. The following step consists in the exact solution of such a linear problem by means of the complex elastic potentials. As result we obtain the exact closed forms for the linear and nonlinear effective elastic moduli, which are valid for any volume fraction of the core embedded in the external shell.     

Exact Solutions for a Thick Elastic Plate with a Thin Elastic Surface Layer

A procedure has been developed in previous papers for constructing exact solutions of the equations of linear elasticity in a plate (not necessarily thin) of inhomogeneous isotropic linearly elastic material in which the elastic moduli depend in any specified manner on a coordinate normal to the plane of the plate. The essential idea is that any solution of the classical equations for a hypothetical thin plate or laminate (which are two-dimensional theories) generates, by straightforward substitutions, a solution of the three-dimensional elasticity equations for the inhomogeneous material. In this paper we consider a thick plate of isotropic elastic material with a thin surface layer of different isotropic elastic material. It is shown that the interface tractions and in-plane stress discontinuities are determined only by the initial two-dimensional solution, without recourse to the three-dimensional elasticity theory. Two illustrative examples are described.     

Integration of the Equilibrium Equations for Inhomogeneous, Transversely Isotropic Plates

A system of differential equilibrium equations for inhomogeneous transversely isotropic plates is derived based on the Fourier series in terms of Legendre polynomials. It is assumed that Poisson's ratios are constant and the elastic moduli are linear functions of the transverse coordinate. A method of finding the general solution to the system of equations derived is set forth     

拉压不同弹性模量薄板上圆孔的应力分析

采用弹性理论研究了拉压不同弹性模量薄板上圆孔的孔边应力集中问题．采用广义虎克定律推导出了拉压不同弹性模量薄板上圆孔边的应力平衡方程,并联合利用应力函数及边界条件得到了拉压不同弹性模量薄板上圆孔边的应力表达式．算例分析表明,当薄板材料的拉压弹性模量相差较大时,采用经典弹性理论研究薄板上圆孔的孔边应力是不合适的,当经典弹性理论与拉压不同弹性模量弹性理论的计算结果间的差别超过工程允许误差5%时,应该采用拉压不同弹性模量弹性理论进行计算．     

Evaluation of the second order elastic moduli

The expressions of the apparent linear elastic moduli and their first and second derivatives, with respect to hydrostatic pressure, are obtained according to the second order elasticity theory. As a particular case when the material is hyperelastic, formulae of the first derivatives of the linear elastic moduli reduce to those obtained by Seeger and Buck.     

Rigorous bounds on the effective moduli of composites and inhomogeneous bodies with negative-stiffness phases

We review the theoretical bounds on the effective properties of linear elastic inhomogeneous solids (including composite materials) in the presence of constituents having non-positive-definite elastic moduli (so-called negative-stiffness phases). Using arguments of Hill and Koiter, we show that for statically stable bodies the classical displacement-based variational principles for Dirichlet and Neumann boundary problems hold but that the dual variational principle for traction boundary problems does not apply. We illustrate our findings by the example of a coated spherical inclusion whose stability conditions are obtained from the variational principles. We further show that the classical Voigt upper bound on the linear elastic moduli in multi-phase inhomogeneous bodies and composites applies and that it imposes a stability condition: overall stability requires that the effective moduli do not surpass the Voigt upper bound. This particularly implies that, while the geometric constraints among constituents in a composite can stabilize negative-stiffness phases, the stabilization is insufficient to allow for extreme overall static elastic moduli (exceeding those of the constituents). Stronger bounds on the effective elastic moduli of isotropic composites can be obtained from the Hashin–Shtrikman variational inequalities, which are also shown to hold in the presence of negative stiffness.     

Influence of interfaces on effective properties of nanomaterials with stochastically distributed spherical inclusions

The method of conditional moments is generalized to include evaluation of the effective elastic properties of particulate nanomaterials and to investigate the size effect in those materials. Determining the effective constants necessitates finding a stochastically averaged solution to the fundamental equations of linear elasticity coupled with surface/interface conditions (Gurtin–Murdoch model). To obtain such a solution the system of governing stochastic differential equations is first transformed to an equivalent system of stochastic integral equations. Using statistical averaging, the boundary-value problem is then converted to an infinite system of linear algebraic equations. A two-point approximation is considered and the stress fluctuations within the inclusions are neglected in order to obtain a finite system of algebraic equations in terms of component-average strains. Closed-form expressions are derived for the effective moduli of a composite consisting of a matrix and randomly distributed spherical inhomogeneities. As a numerical example a nanoporous material is investigated assuming a model in which the interface effects influence only the bulk modulus of the material. In that model the resulting shear modulus is the same as for the material without surface effects. Dependence of the bulk moduli on the radius of nanopores and on the pore volume fraction is analyzed. The results are compared to, and discussed in the context of other theoretical predictions.     

FEA in elasticity of random structure composites reinforced by heterogeneities of non canonical shape

One considers linearly elastic composite media, which consist of a homogeneous matrix containing a statistically homogeneous random set of aligned homogeneous heterogeneities of non canonical shape. Effective elastic moduli as well as the first statistical moments of stresses in the phases are estimated. The explicit new representations of the effective moduli and stress concentration factors are expressed through some building block described by numerical solution for one heterogeneity inside the infinite medium subjected to homogeneous remote loading. The method uses as a background a new general integral equation proposed in Buryachenko, 2010a, Buryachenko, 2010b, which incorporates influence of stress inhomogeneity inside the inclusion on the effective field and makes it possible to reconsider basic concepts of micromechanics such as effective field hypothesis, quasi-crystalline approximation, and the hypothesis of “ellipsoidal symmetry”. The results of this reconsideration are quantitatively estimated for some modeled composite reinforced by aligned homogeneous heterogeneities of non canonical shape. Some new effects are detected that are impossible in the framework of a classical background of micromechanics.     

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.     

Bootstrap Elasticity: From Linear to Nonlinear Constitutive Representations

The normal modes of the 6×6 symmetric matrix of elastic moduli for linear anisotropic elasticity, also called Kelvin modes, provide orthogonal basis sets for the six dimensional space of symmetric, second order tensors in three dimensional Euclidean space. In turn the partitioning of the six space, induced by these bases and the multiplicity of each eigenvalue, provides the means for constructing six term minimal representations of nonlinear constitutive equations for materials of any symmetry from triclinic to cubic. The constructions also for the first time show clear connections to the linear elastic moduli, which through the eigenvalues set the scale for most, but not all, of the tensor generators. This approach also provides an alternate way to construct the well-known three term Rivlin-Ericksen representation for nonlinear isotropic materials.     

含柔性涂层的颗粒增强复合材料弹性模量估计

采用线弹簧型弱界面模型来模拟柔性涂层,研究柔性涂层对复合材料宏观弹性模量的影响。首先利用Mori-Tanaka方法和弱界面球形夹杂问题的弹性解,获得单夹杂内部的平均应力和平均应变,进而求得具有柔性涂层的复合材料的宏观弹性模量,并研究界面柔度对复合材料弹性模量的影响。     

Well-posedness of the fundamental boundary value problems for constrained anisotropic elastic materials

We consider the equations of linear homogeneous anisotropic elasticity admitting the possibility that the material is internally constrained, and formulate a simple necessary and sufficient condition for the fundamental boundary value problems to be well-posed. For materials fulfilling the condition, we establish continuous dependence of the displacement and stress on the elastic moduli and ellipticity of the elasticity system. As an application we determine the orthotropic materials for which the fundamental problems are well-posed in terms of their Young's moduli, shear moduli, and Poisson ratios. Finally, we derive a reformulation of the elasticity system that is valid for both constrained and unconstrained materials and involves only one scalar unknown in addition to the displacements. For a two-dimensional constrained material a further reduction to a single scalar equation is outlined.This paper is dedicated to Professor Joachim Nitsche on the occasion of his sixtieth birthday     

弹性地基上矩形薄板问题的Hamilton正则方程及解析解

利用辛算法求出弹性地基上矩形薄板问题的解析解,将弹性地基视为双参数弹性地基,直接从弹性矩形薄板的控制方程推导出了问题的Hamilton正则方程,为求出任意边界条件下问题的理论解奠定了基础,并且通过算例验证了文中所采用方法的正确性．     

The Relaxation of Two-well Energies with Possibly Unequal Moduli

The elastic energy of a multiphase solid is a function of its microstructure. Determining the infimum of the energy of such a solid and characterizing the associated optimal microstructures is an important problem that arises in the modeling of the shape memory effect, microstructure evolution, and optimal design. Mathematically, the problem is to determine the relaxation under fixed phase fraction of a multiwell energy. This paper addresses two such problems in the geometrically linear setting. First, in two dimensions, we compute the relaxation under fixed phase fraction for a two-well elastic energy with arbitrary elastic moduli and transformation strains, and provide a characterization of the optimal microstructures and the associated strain. Second, in three dimensions, we compute the relaxation under fixed phase fraction for a two-well elastic energy when either (1) both elastic moduli are isotropic, or (2) the elastic moduli are well ordered and the smaller elastic modulus is isotropic. In both cases we impose no restrictions on the transformation strains. We provide a characterization of the optimal microstructures and the associated strain. We also compute a lower bound that is optimal except possibly in one regime when either (1) both elastic moduli are cubic, or (2) the elastic moduli are well ordered and the smaller elastic modulus is cubic; for moduli with arbitrary symmetry we obtain a lower bound that is sometimes optimal. In all these cases we impose no restrictions on the transformation strains and whenever the bound is optimal we provide a characterization of the optimal microstructures and the associated strain. In both two and three dimensions the quasiconvex envelope of the energy can be obtained by minimizing over the phase fraction. We also characterize optimal microstructures under applied stress.     

The number of elastic moduli for anisotropic polynomial constitutive equations

This paper shows how to compute the number of elastic moduli for polynomial anisotropy constitutive equations. We investigate the algebraic structure of constitutive equations and we show that the knowledge of the number of elasticities for five types of anisotropy (described by cyclic groups) permits to determine the number of elasticities for all other symmetry groups relevant for elasticity theory. We compute the numbers of elastic moduli for polynomial constitutive relations of degree ≤6.     

弹性矩形板问题的Hamilton正则方程

为了采用辛算法求出弹性矩形板问题的解析解，中直接从弹性矩形板的控制方程出发推导了弹性矩形板，其中包括弹性矩形薄板和厚板问题以及弹性地基上矩形薄板和厚板问题的Hamilton正则方程，为利用辛几何方法求出任意边界条件下这类问题的理论解奠定了基础．     

聚氨酯泡沫塑料压缩杨氏模量的理论预测

通过微分法导出了泡沫塑料剪切模量和体积模量所满足的微分方程组，再利用联系泡沫塑料泊松比和孔隙比的Ｋｅｒｎｅｒ－Ｒｕｓｃｈ经验关系及泡沫塑料弹性常数间满足足的关系，在基体材料不可压缩的假设下，确定了泡沫塑料的杨氏模量。本文针对几种密度的泡沫塑料，分别对它们的杨氏模量进行了理论预测和实验测定，结果表明：理论预测的模量在较高密度下与实验符合的很好，在低密度下也给出相当好的近似值。此外，本文的结果同其他理     

State dependent linear moduli in ferroelectrics

A novel experimental method is used to measure the evolution of the linear elastic, dielectric and piezoelectric moduli of a soft ferroelectric ceramic during loading. The applied loading states are combinations of uniaxial compressive stress and electric field. Short pulses of electric field and stress are used to increment the remanent strain and polarization state of the material, while the rates of change of electric displacement and strain during unloading are used to assess the moduli. The remanent quantities are treated as state variables, with a view to expressing the moduli as functions of the material state. The piezoelectric moduli are found to vary approximately linearly with polarization, regardless of the remanent strain state, whilst the dielectric moduli and elastic compliances show more complex behaviour. A simple model of the state dependence of the moduli, based on varying the volume fractions of six crystal variants in the tetragonal system, is used to interpret the results.