A constitutive model for orthotropic elasto-plasticity at large strains

Summary The paper presents a thermodynamically consistent constitutive model for elasto-plastic analysis of orthotropic materials at large strain. The elastic and plastic anisotropies are assumed to be persistent in the material but the anisotropy axes can undergo a rigid rotation due to large plastic deformations. The orthotropic yield function is formulated in terms of the generally nonsymmetric Mandel stress tensor such that its skew-symmetric part is additionally taken into account. Special attention is focused on the convexity of the yield surface resulting in the nine-dimensional stress space. Of particular interest are new convexity conditions which do not appear in the classical theory of anisotropic plasticity. They impose additional constraints on the material constants governing the plastic spin. The role of the plastic spin is further studied in simple shear accompanied by large elastic and large plastic deformations. If the plastic spin is neglected, the shear stress response is characterized by oscillations with an amplitude strictly dependent on the degree of the plastic anisotropy.accepted for publication 2 March 2004     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

A FE formulation for elasto-plastic materials with planar anisotropic yield functions and plastic spin

In this article a stress integration algorithm for shell problems with planar anisotropic yield functions is derived. The evolution of the anisotropy directions is determined on the basis of the plastic and material spin. It is assumed that the strains inducing the anisotropy of the pre-existing preferred orientation are much larger than subsequent strains due to further deformations. The change of the locally preferred orientations to each other during further deformations is considered to be neglectable. Sheet forming processes are typical applications for such material assumptions. Thus the shape of the yield function remains unchanged. The size of the yield locus and its orientation is described with isotropic hardening and plastic and material spin.The numerical treatment is derived from the multiplicative decomposition of the deformation gradient and thermodynamic considerations in the intermediate configuration. A common formulation of the plastic spin completes the governing equations in the intermediate configuration. These equations are then pushed forward into the current configuration and the elastic deformation is restricted to small strains to obtain a simple set of constitutive equations. Based on these equations the algorithmic treatment is derived for planar anisotropic shell formulations incorporating large rotations and finite strains. The numerical approach is completed by generalizing the Return Mapping algorithm to problems with plastic spin applying Hill’s anisotropic yield function. Results of numerical simulations are presented to assess the proposed approach and the significance of the plastic spin in the deformation process.     

Modelling inelastic behaviour of orthotropic metals in a unique alignment of deviatoric plane within the stress space

A finite strain constitutive model to predict the deformation behaviour of orthotropic metals is developed in this paper. The important features of this constitutive model are the multiplicative decomposition of the deformation gradient and a new Mandel stress tensor combined with the new stress tensor decomposition generalized into deviatoric and spherical parts. The elastic free energy function and the yield function are defined within an invariant theory by means of the structural tensors. The Hill’s yield criterion is adopted to characterize plastic orthotropy, and the thermally micromechanical-based model, Mechanical Threshold Model (MTS) is used as a referential curve to control the yield surface expansion using an isotropic plastic hardening assumption. The model complexity is further extended by coupling the formulation with the shock equation of state (EOS). The proposed formulation is integrated in the isoclinic configuration and allows for a unique treatment for elastic and plastic anisotropy. The effects of elastic anisotropy are taken into account through the stress tensor decomposition and plastic anisotropy through yield surface defined in the generalized deviatoric plane perpendicular to the generalized pressure. The proposed formulation of this work is implemented into the Lawrence Livermore National Laboratory-DYNA3D code by the modification of several subroutines in the code. The capability of the new constitutive model to capture strain rate and temperature sensitivity is then validated. The final part of this process is a comparison of the results generated by the proposed constitutive model against the available experimental data from both the Plate Impact test and Taylor Cylinder Impact test. A good agreement between experimental and simulation is obtained in each test.     

An anisotropic elastoplastic constitutive formulation generalised for orthotropic materials

This paper presents a finite strain constitutive model to predict a complex elastoplastic deformation behaviour that involves very high pressures and shockwaves in orthotropic materials using an anisotropic Hill’s yield criterion by means of the evolving structural tensors. The yield surface of this hyperelastic–plastic constitutive model is aligned uniquely within the principal stress space due to the combination of Mandel stress tensor and a new generalised orthotropic pressure. The formulation is developed in the isoclinic configuration and allows for a unique treatment for elastic and plastic orthotropy. An isotropic hardening is adopted to define the evolution of plastic orthotropy. The important feature of the proposed hyperelastic–plastic constitutive model is the introduction of anisotropic effect in the Mie–Gruneisen equation of state (EOS). The formulation is further combined with Grady spall failure model to predict spall failure in the materials. The proposed constitutive model is implemented as a new material model in the Lawrence Livermore National Laboratory (LLNL)-DYNA3D code of UTHM’s version, named Material Type 92 (Mat92). The combination of the proposed stress tensor decomposition and the Mie–Gruneisen EOS requires some modifications in the code to reflect the formulation of the generalised orthotropic pressure. The validation approach is also presented in this paper for guidance purpose. The \({\varvec{\psi }}\) tensor used to define the alignment of the adopted yield surface is first validated. This is continued with an internal validation related to elastic isotropic, elastic orthotropic and elastic–plastic orthotropic of the proposed formulation before a comparison against range of plate impact test data at 234, 450 and \({\mathrm {895\,ms}}^{\mathrm {-1}}\) impact velocities is performed. A good agreement is obtained in each test.     

On the implementation of anisotropic yield functions into finite strain problems of sheet metal forming

In this article the implementation of anisotropic yield functions into finite element investigations of orthotropic sheets with planar anisotropy is discussed within a plane-stress context. Special attention is focused on the proper treatment of the orientation of the anisotropic axes during deformation into the finite-strain range. As an example problem the hydrostatic bulging of a membrane is considered in conjunction with a recently proposed anisotropic yield function. It is shown that the aspects of the plane-stress assumption, which do not come into consideration in isotropic analyses, can play an important role on the accuracy of the solution when the rotation of the orthotropic axes enters the computation directly due to the presence of material anisotropy. When the material anisotropy is considered and when the deformation of the workpiece is not limited to the plane of the undeformed sheet (such as cup drawing, hydrostatic bulging, etc.), the numerical experiments indicate that the only correct formulation is the one based on numerically imposing the requirement that for the plane-stress application, the in-plane material axes have to remain in the plane of the sheet during the deformation.     

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.     

Modeling multi-axial deformation of planar anisotropic elasto-plastic materials,part I: Theory

In sheet metal forming processes local material points can experience multi-axial and multi-path loadings. Under such loading conditions, conventional phenomenological material formulations are not capable to predict the deformation behavior within satisfying accuracy. While micro-mechanical models have significantly improved the understanding of the deformation processes under such conditions, these models require large sets of material data to describe the micromechanical evolution and quite enormous computation expenses for industrial applications. To reduce the drawbacks of phenomenological material models under the multi-path loadings a new anisotropic elasto-plastic material formulation is suggested. The model enables the anisotropic yield surface to grow (isotropic hardening), translate (kinematic hardening) and rotate (rotation of the anisotropy axes) with respect to the deformation, while the shape of the yield surface remains essentially unchanged.Essentially, the model is formulated on the basis of an Armstrong–Frederick type kinematic hardening, the plastic spin theory for the reorientation of the symmetry axes of the anisotropic yield function, and additional terms coupling these expressions. The capability of the model is illustrated with multi-path loading simulations in ‘tension-shear’ and ‘reverse-shear’ to assess its performance with ‘cross’ hardening and ‘Bauschinger’ effects.     

正交各向异性薄板的弹塑性屈曲分析

本文应用混合硬化正交各向异性塑性理论和屈曲的能量法则，推导了正交各向异性薄板在面内压缩情况下的弹塑性稳定方程，计算了相应的临界荷载，并讨论几何形状、边界条件及诱导荷载比等对临界应力的影响。分析表明材料的各向异性性质对弹塑性临界应力的有较大的影响。     

Determination of Anisotropic Plastic Constitutive Parameters Using the Virtual Fields Method

The aim of the present study is to retrieve all the anisotropic plastic constitutive parameters from uniaxial loading. A complex geometry which can provide very heterogeneous stress states in a uniaxial tensile test was chosen for steel sheet specimens. A digital image correlation technique was used for the full-field heterogeneous strain measurement. The orthotropic Hill1948 yield criterion with Swift isotropic hardening was adopted as an elasto-plastic constitutive model. The virtual fields method (VFM) was employed as an inverse analytical tool to determine the constitutive parameters. All the parameters were successfully identified using the VFM by combining two tensile test results obtained in rolling and transverse directions.     

Damage in edge cracking of rolled metal slabs

Materials get damaged under shear deformations. Edge cracking is one of the most serious damage to the metal rolling industry, which is caused by the shear damage process and the evolution of anisotropy. To investigate the physics of the edge cracking process, simulations of a shear deformation for an orthotropic plastic material are performed. To perform the simulation, this paper proposes an elasto-aniso-plastic constitutive model that takes into account the evolution of the orthotropic axes by using a bases rotation formula, which is based upon the slip process in the plastic deformation. It is found through the shear simulation that the void can grow in shear deformations due to the evolution of anisotropy and that stress triaxiality in shear deformations of (induced) anisotropic metals can develop as high as in the uniaxial tension deformation of isotropic materials, which increases void volume. This echoes the same physics found through a crystal plasticity based damage model that porosity evolves due to the grain-to-grain interaction. The evolution of stress components, stress triaxiality and the direction of the orthotropic axes in shear deformations are discussed.     

The stress concentration near a rigid line inclusion in a prestressed, elastic material. Part I.: Full-field solution and asymptotics

A lamellar (zero-thickness) rigid inclusion, so-called ‘stiffener’, is considered embedded in a uniformly prestressed (or prestrained), incompressible and orthotropic elastic sheet, subject to a homogeneous far-field deformation increment. This problem is solved under the assumption of plane strain deformation, with prestress principal directions and orthotropy axes aligned with the stiffener. A full-field solution is obtained solving the Riemann-Hilbert problem for symmetric incremental loading at infinity (while for shear deformation the stiffener leaves the ambient field unperturbed). In addition to the full-field solution, the asymptotic Mode I near-tip representation involving the corresponding incremental stress intensity factor are derived and these results are complemented with the Mode II asymptotic solution. For null prestress, the full-field stress state is shown to match correctly with photoelastic experiments performed by us (on two-part epoxy resin samples containing an aluminum lamina). Our experiments also confirm the fracture patterns for a brittle material containing a stiffener, which do not obey a hoop-stress criterion and result completely different from those found for cracks. Issues related to shear band formation and evaluation of energy release rate for a stiffener growth (or reduction) are deferred to Part II of this article.     

塑性各向异性板料本构关系及冲压成形的数值模拟

推导了具有一般屈服函数形式的弹塑性速率型本构关系;给出了用于板料成形的Hill塑性各向异性屈服模型下本构关系的具体形式;用有限元动力显式计算程序MSC/DYTRAN模拟了金属板料的冲压成形;通过算例分析,考察了塑性各向异性对凸耳形成和大小以及对成形模拟结果准确性的影响;数值结果和实验结果表明：各向（厚向）异性本构模型比各向同性本构模型更真实地反映了板料的成形性。     

A framework for numerical integration of crystal elasto-plastic constitutive equations compatible with explicit finite element codes

In this paper we summarize the elements of a numerical integration scheme for elasto-plastic response of single crystals. This is intended to be compatible with large-scale explicit finite element codes and therefore can be used for problems involving multiple crystals and also overall behavior of polycrystalline materials. The steps described here are general for anisotropic elastic and plastic response of crystals. The crystallographic axes of the lattice are explicitly stored and updated at each time step. A plastic predictor–elastic corrector scheme is used to calculate the plastic strain rates on all active slip systems based on a rate-dependent physics-based constitutive model without the need of further auxiliary assumptions. Finally we present the results of numerous calculations using a physics-based rate- and temperature-dependent model of copper and the effect of elastic unloading, elastic crystal anisotropy, and deformation-induced lattice rotation are emphasized.     

On the load transfer from elastic inclusions to an infinite elastic orthotropic plane with cuts

The present paper deals with the problem of load transfer from elastic inclusions to an infinite elastic orthotropic plane with cuts located on one of the principal orthotropy directions. The constitutive system of equations of this problem is derived under the assumption that the inclusions are in a uniaxial stress state. The obtained system consists of a singular integro-differential equation and a singular integral equation for the jumps of the tangential stresses acting on the inclusion shores and for the derivative of the the cut opening function. The behavior of solutions of the system of constitutive equations at the endpoints of the inclusions and cuts is studied, and the solution of this system is constructed by the numerical-analytic discrete singularity method.     

Plastic spin: necessity or redundancy?

The plastic spin concept in large deformation anisotropic elastoplasticity theories with tensorial internal variables, is proved to be a necessary constitutive ingredient. Different inaccurate notions about the plastic spin are dispelled, and its presence in the theory is demystified as something very simple and straightforward. To this extent it is necessary to disassociate the plastic spin concept and the conjugate notion of constitutive spin from the foundation of kinematics, which caused confusion in the past, and define it only in relation to the constitutive equations of evolution of the tensorial internal variables. There, the plastic spin is related to the orientation aspect of such constitutive equations, and the multiplicity of the different internal variables suggests the necessity to have a different spin for each variable. In the process, a straightforward constitutive framework is developed which is based on classical hyperelasticity, yield criteria and invariance requirements of the constitutive functions under superposed rigid body rotation. Ad-hoc assumptions about stress corotational or convected rates and other fuzzy suggestions for different spins are not part of this development. Other topics such as the concept and simplifying effect of the spinless unstressed configuration and its comparison with the isoclinic configuration, some computational aspects, and the effect of small elastic strains are discussed, and all along the significance of plastic spin in the different equations is evaluated.     

Effect of Laminar Orthotropic Myofiber Architecture on Regional Stress and Strain in the Canine Left Ventricle

Recent morphological studies have demonstrated a laminar (sheet) organization of ventricular myofibers. Multiaxial measurements of orthotropic myocardial constitutive properties have not been reported, but regional distributions of three-dimensional diastolic and systolic strains relative to fiber and sheet axes have recently been measured in the dog heart by Takayama et al. [30]. A three-dimensional finite-deformation, finite element model was used to investigate the effects of material orthotropy on regional mechanics in the canine left ventricular wall at end-diastole and end-systole. The prolate spheroidal model incorporated measured transmural distributions of fiber and sheet angles at the base and apex. Compared with transverse isotropy, the orthotropic model of passive myocardial properties yielded improved agreement with measured end-diastolic strains when: (1) normal stiffness transverse to the muscle fibers was increased tangent to the sheets and decreased normal to them; (2) shear coefficients were increased within sheet planes and decreased transverse to them. For end-systole, orthotropic passive properties had little effect, but three-dimensional systolic shear strain distributions were more accurately predicted by a model in which significant active systolic stresses were developed in directions transverse to the mean fiber axis as well as axial to them. Thus the ventricular laminar architecture may give rise to anisotropic material properties transverse to the fibers with greater resting stiffness within than between myocardial laminae. There is also evidence that intact ventricular muscle develops significant transverse stress during systole, though it remains to be seen if active stress is also orthotropic with respect to the laminar architecture. This revised version was published online in July 2006 with corrections to the Cover Date.