Symmetry classes of the anisotropy tensors of quasielastic materials and a generalized Kelvin approach

The anisotropy matrices (tensors) of quasielastic (Cauchy-elastic) materials were obtained for all classes of crystallographic symmetries in explicit form. The fourth-rank anisotropy tensors of such materials do not have the main symmetry, in which case the anisotropy matrix is not symmetric. As a result of introducing various bases in the space of symmetric stress and strain tensors, the linear relationship between stresses and strains is represented in invariant form similar to the form in which generalized Hooke’s law is written for the case of anisotropic hyperelastic materials and contains six positive Kelvin eigen moduli. It is shown that the introduction of modified rotation-induced deformation in the strain space can cause a transition to the symmetric anisotropy matrix observed in the case of hyperelasticity. For the case of transverse isotropy, there are examples of determination of the Kelvin eigen moduli and eigen bases and the rotation matrix in the strain space. It is shown that there is a possibility of existence of quasielastic media with a skew-symmetric anisotropy matrix with no symmetric part. Some techniques for the experimental testing of the quasielasticity model are proposed.     

A constitutive theory for creep behavior of initially isotropic materials sustaining unilateral damage

The goal of the present work is to modify structure of the creep constitutive equations existing in the literature, and simultaneously to incorporate both damage induced anisotropy and unilateral damage into the constitutive model. The proposed nonlinear-tensor constitutive equation for creep together with the damage evolution equation take into account the secondary and tertiary creep of the initially isotropic materials. The material parameters of the model are determined using basic experiments. It is shown that the creep model is capable of describing available experimental data for the lateral creep responses under uniaxial compression.     

Applied creep theory for bodies with anisotropy due to plastic prestrain

Plastic strains in structures at the stages of manufacturing, testing, and approaching the operation regime cause anisotropic variations in the mechanical properties of materials, including creep strength. We consider the following special but practically important class of loading processes for originally isotropic materials: a simple active plastic strain is followed by a long-term steady-state loading within the elastic limits. To describe the second stage, we present the creep strain deviator in the form of an additive orthogonal decomposition in the directions of the repeated loading and the vector anisotropy. The coefficients in the decomposition are material functions of time, of the intensities of the preliminary and repeated loadings, and of the angle between the directions of these loadings. We obtain conditions on the material functions under which, at any given time instant, there is a one-to-one continuous correspondence between the stress and strain tensors for the model proposed and the boundary-value problem in the generalized statement has a unique solution; we also prove the convergence of the iteration method of elastic solutions used to find this unique solution. The model is identified according to the creep diagrams (under steady-state stresses of different values) determined for the material in the original state and after the plastic prestrain at an angle (zero, extended, and intermediate) to the direction of the repeated loading. We show that our results are in good agreement with the results available in the literature concerning experiments in this class of processes for stainless steel at high temperature. We propose an engineering version of the theory in which only the experimental data for uniaxial tension are used. We discuss the versions of the model for the cases in which the plastic preloading is cyclic (one-dimensional or circular) and the repeated loading is unsteady.     

Large strain creep analysis of thick-walled cylinders

The finite-strain theory has been used to study the creep behaviour of a thick-walled cylinder under large strains. The analysis is divided into two parts. In part 1 the creep deformation of a thick-walled cylinder of an anisotropic material subjected to internal pressure has been discussed. The effect of the anisotropy has been depicted graphically. It is found that the anisotropy of the material has a significant effect on the axial stress, strain and strain rate. Part 2 of the paper deals with the creep analysis of cylinders of either isotropic or anisotropic materials subjected to combined internal and external pressures. The effect of the anisotropy is found to be similar to that found in part 1. It is seen, however, that the introduction of external pressure results in decreasing the strain rate and thus increasing the life of the cylinder.     

State Equations of Unsteady Creep under Complex Loading

A mathematical model describing the unsteady creep of metals under complex loading is proposed. The results of numerical simulation of creep of St.304 steel in complex regimes of block multiaxial cyclic deformation are given. The numerical calculation results obtained are compared with the data of full-scale experiments. Creep is simulated in complex deformation processes accompanied by the rotation of main regions of stress, strain, and creep strain tensors.     

Phase field modeling of fracture in anisotropic brittle solids

A phase field model of fracture that accounts for anisotropic material behavior at small and large deformations is outlined within this work. Most existing fracture phase field models assume crack evolution within isotropic solids, which is not a meaningful assumption for many natural as well as engineered materials that exhibit orientation-dependent behavior. The incorporation of anisotropy into fracture phase field models is for example necessary to properly describe the typical sawtooth crack patterns in strongly anisotropic materials. In the present contribution, anisotropy is incorporated in fracture phase field models in several ways: (i) Within a pure geometrical approach, the crack surface density function is adopted by a rigorous application of the theory of tensor invariants leading to the definition of structural tensors of second and fourth order. In this work we employ structural tensors to describe transverse isotropy, orthotropy and cubic anisotropy. Latter makes the incorporation of second gradients of the crack phase field necessary, which is treated within the finite element context by a nonconforming Morley triangle. Practically, such a geometric approach manifests itself in the definition of anisotropic effective fracture length scales. (ii) By use of structural tensors, energetic and stress-like failure criteria are modified to account for inherent anisotropies. These failure criteria influence the crack driving force, which enters the crack phase field evolution equation and allows to set up a modular structure. We demonstrate the performance of the proposed anisotropic fracture phase field model by means of representative numerical examples at small and large deformations.     

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.     

Anisotropic polyconvex energies on the basis of crystallographic motivated structural tensors

In large strain elasticity the existence of minimizers is guaranteed if the variational functional to be minimized is sequentially weakly lower semicontinuous (s.w.l.s.) and coercive. Therefore, polyconvex functions which are always s.w.l.s. are usually considered. For isotropic as well as for transversely isotropic and orthotropic materials constitutive functions that are polyconvex already exist. The main goal of this contribution is to provide a new method for the construction of polyconvex hyperelastic models for more general anisotropy classes. The fundamental idea is the introduction of positive definite second-order structural tensors G=HH^T encoding the anisotropies of the underlying crystal. These tensors can be viewed as a push-forward of a cartesian metric of a fictitious reference configuration to the real reference configuration. Here the driving transformations H in the push-forward operation are mappings of the cartesian base vectors of the fictitious configuration onto crystallographic motivated base vectors. Restrictions of this approach are based on the polyconvexity condition as well as on the usage of second-order structural tensors and pointed out in detail.     

Laws of flow with a limiting gradient in anisotropic porous media

The equations of viscoplastic fluid flow through a porous medium are written for all types of anisotropy. It is shown that in anisotropic media the flows with a limiting gradient are characterized by two material tensors: the tensor of permeability (flow resistance) coefficients and the tensor of limiting gradients. A complex of laboratory measurements for determining the tensors of permeability coefficients and limiting gradients is considered for all types of anisotropic media. It is shown that the tensors of permeability coefficients and limiting gradients are coaxial. Conditions of flow onset and fluid flow laws are formulated for media with monoclinic and triclinic symmetries of flow characteristics.     

考虑颗粒转矩的接触网络诱发各向异性分析

颗粒材料的宏观力学行为与接触网络的组构各向异性密切相关, 根据接触点的滑动与否、转动与否和强弱力情况, 可以将颗粒间的接触系统分为不同的子接触网络. 一般而言, 不同的子接触网络在颗粒体系中的传力机制不同, 对宏观力学响应的贡献也有不同. 采用离散单元法(discrete element method, DEM)模拟了不同抗转动系数$\mu_r$下颗粒材料三轴剪切试验, 分析了剪切过程中不同子接触网络的组构张量的演变规律, 并探究了颗粒抗转动效应对子接触网络各向异性指标演变规律的影响. 研究发现: 剪切过程中转动、非转动接触的组构张量变化不是独立的, 受到颗粒间滑动与否的影响; 非滑动、强接触网络是颗粒间的主要传力结构, 非滑动接触网络的接触法向和法向接触力各向异性均随$\mu_r$的增大而增大, 其对宏观应力的贡献程度随$\mu_r$的增大而减小;强接触网络的接触法向各向异性随$\mu_r$的增大而增大, 但法向接触力各向异性随$\mu_r$的增大无明显变化, 强接触网络对宏观应力的贡献程度在不同$\mu_r$情况下均相同.      

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.     

On the modeling and simulation of induced anisotropy in polycrystalline metals with application to springback

The presence of initial, and the development of induced, anisotropic elastic and inelastic material behavior in polycrystalline metals, can be traced back to the influence of texture and dislocation substructural development on this behavior. As it turns out, via homogenization or other means, one can formulate effective models for such structure and its effect on the macroscopic material behavior with the help of the concept of evolving structure tensors. From the constitutive point of view, these quantities determine the material symmetry properties. Most importantly, all dependent constitutive fields (e.g., stress) are by definition isotropic functions of the independent constitutive variables, which include these evolving structure tensors. The evolution of these tensors during loading results in an evolution of the anisotropy of the material. From an algorithmic point of view, the current approach leads to constitutive models which are quite amenable to numerical implementation. To demonstrate the applicability of the resulting constitutive formulation, we apply it to the case of metal plasticity with combined hardening involving both deformation- and permanently induced anisotropy. Comparison of simulation results based on this model for the bending tension of aluminum-alloy sheet-metal strips with corresponding experimental ones show good agreement.     

Barnett-Lothe tensors and their associated tensors for monoclinic materials with the symmetry plane at x 3=0

The three Barnett-Lothe tensors S, H, L and the three associated tensors S(), H(), L() appear frequently in the real form solutions to two-dimensional anisotropic elasticity problems. Explicit expressions of the components of these tensors are derived and presented for monoclinic materials whose plane of material symmetry is at x ₃=0. We use the algebraic formalism for these tensors but the results are derived not by the straight-forward substitution of the complex matrices A and B into the formulae. Instead, we find the product –AB ^-1, whose real and imaginary parts are SL ^-1 and L ^-1, respectively. The tensors S, H, L are then determined from SL ^-1 and L ^-1. For S(), H(), L() we again avoid the direct substitution by employing an alternate approach. The new approaches require minimal algebra and, at the same time, provide simple and concise expressions for the components of these tensors. Although the new approaches can be extended, in principle, to monoclinic materials whose plane of symmetry is not at x ₃=0 and to materials of general anisotropy, the explicit expressions for these materials are too complicated. More studies are needed for these materials.     

A constitutive model for plastically anisotropic solids with non-spherical voids

Plastic constitutive relations are derived for a class of anisotropic porous materials consisting of coaxial spheroidal voids, arbitrarily oriented relative to the embedding orthotropic matrix. The derivations are based on nonlinear homogenization, limit analysis and micromechanics. A variational principle is formulated for the yield criterion of the effective medium and specialized to a spheroidal representative volume element containing a confocal spheroidal void and subjected to uniform boundary deformation. To obtain closed form equations for the effective yield locus, approximations are introduced in the limit-analysis based on a restricted set of admissible microscopic velocity fields. Evolution laws are also derived for the microstructure, defined in terms of void volume fraction, aspect ratio and orientation, using material incompressibility and Eshelby-like concentration tensors. The new yield criterion is an extension of the well known isotropic Gurson model. It also extends previous analyses of uncoupled effects of void shape and material anisotropy on the effective plastic behavior of solids containing voids. Preliminary comparisons with finite element calculations of voided cells show that the model captures non-trivial effects of anisotropy heretofore not picked up by void growth models.     

Anisotropic finite elastoplasticity with nonlinear kinematic and isotropic hardening and application to sheet metal forming

The paper discusses the derivation and the numerical implementation of a finite strain material model for plastic anisotropy and nonlinear kinematic and isotropic hardening. The model is derived from a thermodynamic framework and is based on the multiplicative split of the deformation gradient in the context of hyperelasticity. The kinematic hardening component represents a continuum extension of the classical rheological model of Armstrong–Frederick kinematic hardening. Introducing the so-called structure tensors as additional tensor-valued arguments, plastic anisotropy can be modelled by representing the yield surface and the plastic flow rule as functions of the structure tensors. The evolution equations are integrated by a new form of the exponential map that preserves plastic incompressibility and uses the spectral decomposition to evaluate the exponential tensor functions in closed form. Finally, the applicability of the model is demonstrated by means of simulations of several deep drawing processes and comparisons with experiments.     

On the forms of relationship between two noncoaxial second-order tensors (the case of plane strain or a plane stress state)

A relationship between noncoaxial tensors of stress and creep strain rate is established for the case of plane strain or a plane stress state. The basis is the experimentally substantiated hypothesis on the existence of a creep surface, which is a set of loading paths in the stress space that, at any time, ensure identical values of the creep intensity for a certain chosen measure and orthogonality of the creep strain rate vector to this surface. The relation obtained completely corresponds to available experimental data for complex loading. Lavrent'ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 2, pp. 135–140, March–April, 1998.     

Description of plastic anisotropy effects at large deformations. Part II: the case of transverse isotropy

In a previous paper (see Tsakmakis, 1999) a general thermodynamically consistent (visco-) plasticity theory has been developed, which accounts for anisotropy effects. For simplicity, isotropic hardening has not be regarded, while anisotropy arises from kinematic hardening and orientational evolution of the underlying substructure. In the present paper the capabilities of this theory are discussed for the study case of transverse isotropy. Anisotropy effects are elaborated in the free energy and the yield function by means of structural tensors. Characteristic features of the transversely isotropic model are illustrated for the case of homogeneous simple shear.     

The spherical vessel with anisotropic creep properties considering large strains

The effect of material anisotropy on creep of pressurized thick-walled spherical vessel has been discussed considering the large strain theory. It is found that the creep strain varies with varying anisotropy of the material. The results obtained for the anisotropic case have been compared with those obtained for the isotropic case and it is observed that the stress and strain distribution in the wall of the vessel is strikingly different for the two cases.