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.     

Anisotropy tensor of the potential model of steady creep

The Kelvin approach describing the structure of the generalized Hooke’s law is used to analyze the potential model of anisotropic creep of materials. The creep equations of incompressible transversely isotropic, orthotropic materials and those with cubic symmetry are considered. The eigen coefficients of anisotropy and eigen tensors for the anisotropy tensors of these materials are determined.     

General invariant representations of the constitutive equations for isotropic nonlinearly elastic materials

This paper develops general invariant representations of the constitutive equations for isotropic nonlinearly elastic materials. Different sets of mutually orthogonal unit tensor bases are constructed from the strain argument tensor by using the representation theorem and corresponding irreducible invariants are defined. Their relations and geometrical interpretations are established in three dimensional principal space. It is shown that the constitutive law linking the stress and strain tensors is revealed to be a simple relationship between two vectors in the principal space. Relative to two different sets of the basis tensors, the constitutive equations are transformed according to the transformation rule of vectors. When a potential function is assumed to exist, the vector associated with the stress tensor is expressed in terms of its gradient with respect to the vector associated with the strain tensor. The Hill’s stability condition is shown to be that the scalar product of the increment of those two vectors must be positive. When potential function exists, it becomes to be that the 3 × 3 constitutive matrix derived from its second order derivative with respect to the vector associated with the strain must be positive definite. By decomposing the second order symmetric tensor space into the direct sum of a coaxial tensor subspace and another one orthogonal to it, the closed form representations for the fourth order tangent operator and its inversion are derived in an extremely simple way.     

Electromechanical response of (3–0, 3–1) particulate,fibrous, and porous piezoelectric composites with anisotropic constituents: A model based on the homogenization method

An analytical framework based on the homogenization method has been developed to predict the effective electromechanical properties of periodic, particulate and porous, piezoelectric composites with anisotropic constituents. Expressions are provided for the effective moduli tensors of n-phase composites based on the respective strain and electric field concentration tensors. By taking into account the shape and distribution of the inclusion and by invoking a simple numerical procedure, solutions for the electromechanical properties of a general anisotropic inclusion in an anisotropic matrix are obtained. While analytical forms are provided for predicting the electroelastic moduli of composites with spherical and cylindrical inclusions, numerical evaluation of integrals over the composite microstructure is required in order to obtain the corresponding expressions for a general ellipsoidal particle in a piezoelectric matrix. The electroelastic moduli of piezoelectric composites predicted by the analytical model developed in the present study demonstrate excellent agreement with results obtained from three-dimensional finite-element models for several piezoelectric systems that exhibit varying degrees of elastic anisotropy.     

A yield function for anisotropic materials Application to aluminum alloys

A phenomenological yield function is proposed to represent the plastic anisotropy of aluminum sheets. It is an extension of the functions given by Barlat et al. [Int. J. Plasticity 7 (1991) 693] and Karafillis and Boyce [J. Mech. Phys. Solids 41 (1993) 1859]. The anisotropy is represented by 12 parameters in the form of two fourth order symmetric tensors. Four other parameters influence the shape of the yield surface uniformly. The role of each parameter is described in detail. The convexity of the yield surface is proved. The implementation of the proposed yield function is done in the 3D general case in an object-oriented finite element code. It is used to represent the anisotropy of a 2024 aluminum thin sheet and the adjustment is excellent. Other anisotropic materials from the literature are also well described by the proposed yield function.     

Solutions of inhomogeneity problems with graded shells and application to core-shell nanoparticles and composites

This paper first presents the Eshelby tensors and stress concentration tensors for a spherical inhomogeneity with a graded shell embedded in an alien infinite matrix. The solution is then specialized to inhomogeneous inclusions in finite spherical domains with fixed displacement or traction-free boundary conditions. The Eshelby tensors in the infinite and finite domains and the stress concentration tensors are especially useful for solving many problems in mechanics and materials science. This is demonstrated on two examples. In the first example, the strain distributions in core-shell nanoparticles with eigenstrains induced by lattice mismatches are calculated using the Eshelby tensors in the finite domains. In the second example, the Eshelby and stress concentration tensors in the three-phase configuration are used to formulate the generalized self-consistent prediction of the effective moduli of composites containing spherical particles within the framework of the equivalent inclusion method. The advantage of this micromechanical scheme is that, whilst its predictions are almost identical to the classical generalized self-consistent method and the third-order approximation, the expressions for the effective moduli have simple closed forms.     

The General Mathematical Theory of Plasticity and the Il’yushin Postulates of Macroscopic Definability and Isotropy

The physical laws characterizing the relation between stresses and strains are considered and analyzed in the general modern theory of elastoplastic deformations and in its postulates of macroscopic definability and isotropy for initially isotropic continuous media. The fundamentals of this theory in continuum mechanics were developed by A.A. Il’yushin in the mid-twentieth century. His theory of small elastoplastic deformations under simple loading became a generalization of Hencky’s deformation theory of flow, whereas his theory of elastoplastic processes which are close to simple loading became a generalization of the Saint-Venant–Mises flow theory to the case of hardening media. In these theories, the concepts of simple arid complex loading processes arid the concept of directing form change tensors are introduced; the Bridgman law of volume elastic change and the universal Roche–Eichinger laws of a single hardening curve under simple loading are adopted; and the Odquist hardening for plastic deformations is generalized to the case of elastoplastic hardening media for the processes of almost simple loading without consideration of a specific history of deformations for the trajectories with small arid mean curvatures. In this paper we discuss the possibility of using the isotropy postulate to estimate the effect of forming parameters in the stress-strain state appeared due to the strain-induced anisotropy during the change of the internal structures of materials. We also discuss the possibility of representing the second-rank symmetric stress and strain tensors in the form of vectors in the linear coordinate six-dimensional Euclidean space. An identity principle is proposed for tensors and vectors.     

On the representation of Prandtl-Reuss tensors within the framework of multiplicative elastoplasticity

This article presents some aspects of the formulation of finite strain elastoplasticity based on the multiplicative decomposition of the deformation gradient. A “canonical” structure of multiplicative elastoplasticity is discussed characterized by a geometrical setting relative to the intermediate configuration in terms of mixed-variant tensors and the exploitation of fundamental dissipation principles. The symmetric fourth-order elastoplastic moduli (so-called ‘Prandtl-Reuss-Tensors’ of the associative theory) appear as a consequence of the assumed metric-dependence of the flow criterion function in a characteristic structure which seems to be typical for large strain multiplicative elastoplasticity. Particular representations of “Prandtl-Reuss tensors” are outlined for isotropic response as well as for decoupled volumetric-isochoric stress response.     

Micro-mechanical behaviors of SMA composite materials under hygro-thermo-elastic strain fields

An analytical procedure to evaluate the behavior of shape memory alloy (SMA) composite under hygrothermal environment is presented. The SMA wires are considered as inclusions embedded in a homogeneous matrix medium of the composite. The inhomogeneity associated with the phase transformation and thermal strains in the SMA wire as well as the hygrothermal strain in the matrix is homogenized using Eshelby’s equivalent inclusion method. In the present work, a similar approach adopted for SMA composites by Marfia and Sacco [Marfia, S., Sacco, E., 2005. Micromechanics and homogenisation of SMA-wire-reinforced materials. J. Appl. Mech. 72 (2), 259–268.] is considered in order to validate the response of SMA composite subjected to thermo-elastic strain field. However, in the present approach, certain modifications and new derivations for the inelastic strain tensors is carried out. First, the constitutive laws for the SMA wire and matrix are expressed in terms of the average strain in the composite. The evolutionary equations used to characterize the pseudoelastic (PE) behavior of the SMA wire are redefined in terms of the eigen strains (phase transformation and thermal strains) occurring in the SMA wire, which are then expressed in terms of the average strain in the composite. Further, the SMA composite constitutive law under coupled hygro-thermo-elastic strain fields is proposed. The generic homogenized hygric and thermal inelastic composite tensors required for the proposed hygro-thermo-elastic constitutive law are derived. Finally, the SMA composite lamina is characterized using Eshelby’s equivalent inclusion method. Using the proposed modifications and derivations, the analytical results are validated for the case of thermo-elastic strain fields and the procedure is then extended to evaluate the SMA composite behavior under hygro-thermo-elastic strain fields. The results include the effect of thermo-elastic and hygro-thermo-elastic strains on the transformation stresses and the nature of hysteresis due to hygric and thermo-elastic strains.     

使用主值空间表示的各向同性塑性本构方程

针对各向同性材料,在内变量为标量的假定下,应用张量函数表示定理给出了其塑性应变增量的不变性表示.它的3个不可约基张量取决于应力张量、相互正交且共主轴.建立3个基张量构成的张量子空间与三维主值空间的对应关系,将共主轴的张量采用笛卡尔坐标系中的矢量描述,矢量在不同坐标系下的分量均为张量的一组不可约不变量.定义塑性应变增量对应的矢量为内变量增量,使用张量函数表示理论得到,内变量演化方程除取决于应力对应的矢量和内变量本身外,还取决于应力增量在张量子空间中的投影,该投影就是应力对应矢量的增量,因此,本构方程归结为确定主值空间中矢量之间的关系.最后表明,三维主值空间与张量子空间中的流动法则是等价的.     

A continuum damage mechanics framework for modeling micro-damage healing

A novel continuum damage mechanics-based framework is proposed to model the micro-damage healing phenomenon in the materials that tend to self-heal. This framework extends the well-known Kachanov’s (1958) effective configuration and the concept of the effective stress space to self-healing materials by introducing the healing natural configuration in order to incorporate the micro-damage healing effects. Analytical relations are derived to relate strain tensors and tangent stiffness moduli in the nominal and healing configurations for each postulated transformation hypothesis (i.e. strain, elastic strain energy, and power equivalence hypotheses). The ability of the proposed model to explain micro-damage healing is demonstrated by presenting several examples. Also, a general thermodynamic framework for constitutive modeling of damage and micro-damage healing mechanisms is presented.     

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

Algebra of Transversely Isotropic Sixth Order Tensors and Solution to Higher Order Inhomogeneity Problems

In this paper we provide a complete and irreducible representation for transversely isotropic sixth order tensors having minor symmetries. Such tensors appear in some practical problems of elasticity for which their inversion is required. For this kind of tensors, we provide an irreducible basis which possesses some remarkable properties, allowing us to provide a representation in a compact form which uses two scalars and three matrices of dimension 2, 3 and 4. It is shown that the calculation of sum, product and inverse of transversely isotropic sixth order tensors is greatly simplified by using this new formalism and appears to be appropriate for deriving new various solutions to some practical problems in mechanics which use such kinds of higher order tensors. For instance, we derive the fields within a cylindrical inhomogeneity submitted to remote gradient of strain. The method of resolution uses the Eshelby equivalent inclusion method extended to the case of a polynomial type eigenstrain. It is shown that the approach leads to a linear system involving a sixth order tensor whose closed form solution is derived by means of the tensorial formalism introduced in the first part of the paper.     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

磁致伸缩材料的非线性本构关系

给出了磁致伸缩材料的两个非线性本构关系,即标准平方型和双曲正切型。在确定一维问题的本构系数时,基于已有的实验结果,引进一个材料函数,用来描述磁致伸缩材料的压磁系数随预应力变化的关系。将非 线性本构关系的理论模型计算结果与实验曲线对比,结果表明标准平方型本构关系在中低磁场下能精确地模拟实验曲线,而双曲正切型本构关系在高磁场时能反映材料的磁致应变饱和现象。讨论了在标准平方型本构的一般三维情形,给出了确定本构系数的方法。     

Basic Issues Concerning Finite Strain Measures and Isotropic Stress-Deformation Relations

It is indicated that the commonly-used Rivlin–Ericksen representation formula for isotropic tensor functions exhibits some properties that might be undesirable for its reasonable and effective applications. Towards clarification and improvement, a set of three mutually orthogonal tensor generators is introduced to achieve an alternative representation formula for isotropic symmetric tensor-valued functions of a symmetric tensor. This representation formula enables us to express the unknown representative coefficients in terms of simple, explicit tensorial inner products of the argument tensor and the value tensor without involving their eigenvalues. In particular, the tensorial interpolation expressions thus obtained assume a unified form for the three different cases of coalescence of the eigenvalues of the argument tensor. Moreover, each summand in the alternative representation formula is shown to inherit the continuity and differentiability properties of the represented isotropic tensor function. These results are used to study some basic issues concerning finite strain measures and stress-deformation relations of isotropic materials, such as continuity and differentiability properties of the representation, determination of the representative coefficients in terms of experimental data for stress and deformation tensors, and computations of finite strain measures. This revised version was published online in June 2006 with corrections to the Cover Date.     

Analysis and simulation for tensile behavior of anisotropic open-cell elastic foams

Based on the elongated Kelvin obtained to investigate the tensile behavior Kelvin model＇s periodicity and symmetry in model, a simplified periodic structural cell is of anisotropic open-cell elastic foams due to the whole space. The half-strut element and elastic deflection theory are used to analyze the tensile response as done in the previous studies. This study produces theoretical expressions for the tensile stress-strain curve in the rise and transverse directions. In addition, the theoretical results are examined with finite element simulation using an existing formula. The results indicate that the theoretical analysis agrees with the finite element simulation when the strain is not too high, and the present model is better. At the same time, the anisotropy ratio has a significant effect on the mechanical properties of foams. As the anisotropy ratio increases, the tensile stress is improved in the rising direction but drops in the transverse direction under the same strain.