Further study of rubber-like elasticity： elastic potentials matching biaxial data＊

By virtue of the rational interpolation procedure and logarithmic strain, a direct approach is proposed to obtain elastic potentials that exactly match uniaxial data and shear data for elastomers. This approach reduces the determination of multiaxial elastic potentials to that of two one-dimensional potentials, thus bypassing usual cumbersome procedures of identifying a number of unknown parameters. Predictions of the suggested potential are derived for a general biaxial stretch test and compared with the classical data given by Rivlin and Saunders(Rivlin, R. S. and Saunders, D. W. Large elastic deformation of isotropic materials. VII: experiments on the deformation of rubber. Phill. Trans. Royal Soc. London A, 243, 251–288(1951)). Good agreement is achieved with these extensive data.     

Convex Fung-type potentials for biological tissues

The anisotropic, non-linear elastic behavior of biological soft tissue is typically accounted for by the hypothesis of hyperelasticity, i.e., the existence of an elastic potential. Fung-type potentials, based on the exponential of a quadratic form in the components of the Green-Lagrange strain, have been widely used in soft tissue modeling, and have inspired potentials in which the exponential was replaced by other monotonically increasing functions. It has been shown that simple fitting of the parameters of a Fung-type potential to experimental stress-strain curves may lead to non-convexity, with undesirable effects on the reliability of the algorithms used in Finite Element simulations. In this paper, we prove that the necessary and sufficient condition for the strict convexity of a Fung-type potential is that the quadratic form in the exponential is positive definite. This result provides a clear physical meaning for the parameters featuring in the quadratic form, and their relationship with the small-strain elastic moduli. This consistency relationship must be respected in order to guarantee that the Fung-type potential correctly reduces to the quadratic potential of classic linear elasticity in the small-strain approximation. Furthermore, we show that, when the conditions of convexity and consistency with the linear theory are respected, Fung-type potentials become a one-parameter family, and we discuss the consequences of this result for when fitting experimental data. An erratum to this article can be found at     

On multiple circular inclusions in plane magnetoelasticity

A general series solution to the magnetoelastic problem of interacting circular inclusions in plane magnetoelasticity is provided in this paper. By the use of complex variable theory and Laurent series expansion method, the general expression of the magnetic and the magnetoelastic complex potentials for the circular inclusion problem is derived. Expanding the definition of the Airy’s stress function of pure elastic field into the magnetoelastic field and applying the superposition method, the general expression then can be reduced to a set of linear algebraic equations and solved in a series form. An approximate closed form solution for the case of two arbitrarily located inclusions is also provided. For illustrating the effect of the pertinent parameters, the numerical results of the interfacial magnetoelastic stresses are displayed in graphic form.     

Interacting elliptic inclusions by the method of complex potentials

An accurate series solution has been obtained for a piece-homogeneous elastic plane containing a finite array of non-overlapping elliptic inclusions of arbitrary size, aspect ratio, location and elastic properties. The method combines standard Muskhelishvili’s representation of general solution in terms of complex potentials with the superposition principle and newly derived re-expansion formulae to obtain a complete solution of the many-inclusion problem. By exact satisfaction of all the interface conditions, a primary boundary-value problem stated on a complicated heterogeneous domain has been reduced to an ordinary well-posed set of linear algebraic equations. A properly chosen form of potentials provides a remarkably simple form of solution and thus an efficient computational algorithm. The theory developed is rather general and can be applied to solve a variety of composite mechanics problems. The advanced models of composite involving up to several hundred inclusions and providing an accurate account for the microstructure statistics and fiber–fiber interactions can be considered in this way. The numerical examples are given showing high accuracy and numerical efficiency of the method developed and disclosing the way and extent to which the selected structural parameters influence the stress concentration at the matrix–inclusion interface.     

Fourth-order quasi-harmonic equations of state for solids of cubic and tetragonal symmetry

General expressions are given for the dependence of the pressure and the effective elastic moduli on deformation and temperature in the form of a Taylor series expansion with respect to elastic and thermal strains. The temperature dependence of these expressions is derived within the quasi-harmonic approximation of lattice dynamics. The expressions are developed in terms of the Lagrangian strain and an alternative strain measure identical with the Eulerian strain for a pure deformation. They are then used to obtain the third- and fourth-order equations of state for crystals of cubic and tetragonal symmetry and to relate the parameters entering these equations to quantities which are commonly (or may be potentially) measured experimentally. It is shown that available ultrasonic data are not completely sufficient to evaluate the parameters of fourth-order equations of state. For tetragonal symmetry, this problem is still in abeyance; while in the cubic case, it is possible to estimate the fourth-order parameters from shock-wave data and so to give illustrative numerical applications of our equations. Finally, the third- and fourth-order Hugoniots and isotherms of Cu and Ag are calculated in terms of both the Lagrangian and Eulerian strain measures.     

一种土的非线性弹性本构模型参数的反演方法

旨在提出一种土的非线性弹性本构模型参数反演的方法。以现今普遍实行的地基载荷试验为基础,依据遗传算法的组合优化理论,采用正演计算和遗传算法优化相结合的方式,建立了土层非线性弹性本构模型参数反演的方法;并依据某黄土场地地基载荷试验数据,实施了黄土土层非线性弹性本构模型参数反演的全过程。计算结果表明,所建立的方法可以实现土层非线性弹性本构模型中相互关联的多个参数的组合优化,并在对初始值要求较低的情况下,可以获得良好的参数反演精度。从而为土的变形特性分析和土与其中及相邻结构的共同作用分析,提供了较好的土体本构模型参数的确定方法。     

Elastic interaction of point defects in crystals with cubic symmetry

The energy of elastic mechanical interaction between point defects in cubic crystals is analyzed numerically. The finite-element complex ANSYS is used to investigate the character of interaction between point defects depending on their location along the crystallographic directions 〈100〉, 〈110〉, 〈111〉 and on the distance from the free boundary of the crystal. The numerical results are compared with the results of analytic computations of the energy of interaction between two point defects in an infinite anisotropic medium with cubic symmetry. The interaction between compressible and incompressible defects of general type is studied. Conditions for onset of elastic attraction between the defects, which leads to general relaxation of the crystal elastic energy, are obtained.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

The bending moment and springback in pure bending of anisotropic sheets

Finite deformation rigid plastic and elastic–plastic analyses of plane strain pure bending of a plastically anisotropic sheet is presented. An efficient method for finding the exact solution is proposed by extending the previously developed method to the stage of unloading. Using this method the solutions are obtained in closed form or reduced to a numerical treatment of ordinary integrals, or an ordinary differential equation, or transcendental equations. An effect of plastic anisotropy and elastic properties on the bending moment is analyzed. The distribution of residual stresses is illustrated and an effect of material and process parameters on springback is investigated.     

Wave processes in nanostructures formed by nanotube arrays or nanosize crystals

Basic parameters of wave processes in structures formed by arrays of parallel nanosize crystals or nanotubes grown in the direction normal to the substrate are obtained. This problem is considered in modeling the behavior of nanoelectromechanical systems, for instance, sensors that utilize such structures. The parameters obtained can also be used to determine the effective elastic characteristics of nanoobjects forming this structure.     

Complex unloading behavior: Nature of the deformation and its consistent constitutive representation

Complex (nonlinear) unloading behavior following plastic straining has been reported as a significant challenge to accurate springback prediction. More fundamentally, the nature of the unloading deformation has not been resolved, being variously attributed to nonlinear/reduced modulus elasticity or to inelastic/“microplastic” effects. Unloading-and-reloading experiments following tensile deformation showed that a special component of strain, deemed here “Quasi-Plastic-Elastic” (“QPE”) strain, has four characteristics. (1) It is recoverable, like elastic deformation. (2) It dissipates work, like plastic deformation. (3) It is rate-independent, in the strain rate range 10⁻⁴-10⁻²/s, contrary to some models of anelasticity to which the unloading modulus effect has been attributed. (4) To first order, the evolution of plastic properties occurs during QPE deformation. These characteristics are as expected for a mechanism of dislocation pile-up and relaxation. A consistent, general, continuum constitutive model was derived incorporating elastic, plastic, and QPE deformation. Using some aspects of two-yield-function approaches with unique modifications to incorporate QPE, the model was implemented in a finite element program with parameters determined for dual-phase steel and applied to draw-bend springback. Significant differences were found compared with standard simulations or ones incorporating modulus reduction. The proposed constitutive approach can be used with a variety of elastic and plastic models to treat the nonlinear unloading and reloading of metals consistently for general three-dimensional problems.     

Auxetic mechanics of crystalline materials

In the present paper, we analyze uniaxial deformation of crystals of different systems with negative Poisson’s ratios, known as auxetics. The behavior of auxetic crystals is studied on the basis of extensive knowledge on the experimental values of elastic constants of different crystals, gathered in the well-known Landolt-Börnstein tables. The competition between the anisotropy of crystal structures and the orientation of deformable samples results in the dependence of the elastic characteristics of deformation, such as Young’s modulus and Poisson’s ratio, on the orientation angles. In the special case of a single angle, a large number of auxetics were found among crystals of cubic, hexagonal, rhombohedral, tetragonal, and orthorhombic systems and the character of variations in their response due to changes in orientation was determined.     

Contact problems for elastic bodies with initial stresses (rigid punches)

Conclusions In this work, the results are presented of investigations into contact problems for prestressed solids with special reference to rigid punches. It should be mentioned that in the last 5–6 years further investigations into the contact problems for solids with inifial stresses with special reference to elastic punches using the formulation of the second approach, i. e., in the general united form for compressible and incompressible solids with elastic potentials of the arbitrary structure, were also carried out by the authors of this article.The authors believe that further progress in the development of the mechanics of contact interaction of prestressed solids (both for rigid and elastic punches) is determined by examination of more complicated problem groups (for elastic, viscoelastic, and plastic solids), and by carrying out experiments to determine the strength of the effect of the initial stresses on the main characteristics of contact interaction of the structural material. It is also evident that, from the practical viewpoint, it is interesting to carry out these investigations for nonhomogeneous initial states. However, this is associated with considerable mathematical problems, even in theoretical developments and, these problems are even greater in application in practice.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Khmel'nitsk Technological Institute of Household Maintenance. Translated from Prikladnaya Mekhanika, Vol. 25, No. 8, pp. 3–18, August, 1989.     

有限变形下多晶晶体塑性模型算法及应用

用Sanna和Zacharia^[1]所提出的延性单晶本构模型的积分算法和Taylor多晶模型假设研究了时间步长和硬化模型的选取对多晶集合体的应力应变响应和织构演化的影响。该算法是利用变形梯度乘法分解获得弹性变形梯度演化方程，用增量迭代法积分该方程，显式更新各滑移上的临界分切剪应力。算例的结果表明该算法具有时间步大，计算效率高的特点，另外，不同硬化模型的选取对多晶集合体应力应变响应的预测有明显的影响但对织构演化的预测影响不大。     

Crack edge displacement and elastic constant determination for an orthotropic material

A method for determining the elastic constants of an isotropic material, based on crack edge displacement data, is extended to an orthotropic material. Complex potentials are used to obtain the stresses and displacements for plane strain. Mode I crack problems in three mutually orthogonal planes are considered and solved. In particular, the expressions of crack edge displacements are obtained in an explicit form. An iterative statistical identification method, based on a Bayesan approach, is used to identify the elastic constants of an orthotropic medium from the Mode I crack displacements measured from the mid-point of the crack. Some graphics are displayed to illustrate the convergence of the pertinent parameters and the approach of the analytical displacements to their experimental values.     

任意轴对称弹性体吸附接触的广义Maugis模型

通过线性叠加Sneddon方法和Lowengrub-Sneddon方法分别给出的解, 得到了一个弹性半空间 轴对称混合边值问题的一般解，进而研究了两个一般轴对称弹性体的正向无摩擦吸附接触问 题. 考虑任意有效的表面形状(要求中心部分首先进入接触)和任意的表面吸附作用，推广 得到了广义Maugis模型. 该模型是一个半解析的模型，它可以分解成表面形状和表面吸附 作用的分别独立影响的两部分，以及一个关联变形和吸附作用的式子. 利用Dugdale模型近 似表面吸附作用，得到了具有任意有效的表面形状的广义M-D模型. 它在强吸附或软材料条 件下的极限形式是广义JKR模型，而在弱吸附或硬材料下的另一个极限形式是广义DMT模型.     

Pair vs many-body potentials: Influence on elastic and plastic behavior in nanoindentation of fcc metals

Molecular-dynamics simulation can give atomistic information on the processes occurring in nanoindentation experiments. In particular, the nucleation of dislocation loops, their growth, interaction and motion can be studied. We investigate how realistic the interatomic potentials underlying the simulations have to be in order to describe these complex processes. Specifically we investigate nanoindentation into a Cu single crystal. We compare simulations based on a realistic many-body interaction potential of the embedded-atom-method type with two simple pair potentials, a Lennard-Jones and a Morse potential. We find that qualitatively many aspects of nanoindentation are fairly well reproduced by the simple pair potentials: elastic regime, critical stress and indentation depth for yielding, dependence on the crystal orientation, and even the level of the hardness. The quantitative deficits of the pair potential predictions can be traced back: (i) to the fact that the pair potentials are unable in principle to model the elastic anisotropy of cubic crystals and (ii) as the major drawback of pair potentials we identify the gross underestimation of the stable stacking fault energy. As a consequence these potentials predict the formation of too large dislocation loops, the too rapid expansion of partials, too little cross slip and in consequence a severe overestimation of work hardening.     

ON LARGE DEFORMATION UNILATERAL CONTACT PROBLEM WITH FRICTION (I)—INCREMENTAL VARIATIONAL EQUATIONS

Based on the theory and technique of nonlinear geometric field theory of continuum, a more general incremental variational equation for elastic and plastic large deformation in co-moving coordinate is established in this paper. An expression for two and three-dimensional continua is derived, and the incremental variational equation for large deformation of changing boundary contact and the variational inequality in rate form are obtained, which provides the theoretical basis for the computation of elastic-plastic large deformation contact problem with friction.