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.     

Homogenization of inelastic solid materials at finite strains based on incremental minimization principles. Application to the texture analysis of polycrystals

The paper presents new continuous and discrete variational formulations for the homogenization analysis of inelastic solid materials undergoing finite strains. The point of departure is a general internal variable formulation that determines the inelastic response of the constituents of a typical micro-structure as a generalized standard medium in terms of an energy storage and a dissipation function. Consistent with this type of finite inelasticity we develop a new incremental variational formulation of the local constitutive response, where a quasi-hyperelastic micro-stress potential is obtained from a local minimization problem with respect to the internal variables. It is shown that this local minimization problem determines the internal state of the material for finite increments of time. We specify the local variational formulation for a distinct setting of multi-surface inelasticity and develop a numerical solution technique based on a time discretization of the internal variables. The existence of the quasi-hyperelastic stress potential allows the extension of homogenization approaches of finite elasticity to the incremental setting of finite inelasticity. Focussing on macro-deformation-driven micro-structures, we develop a new incremental variational formulation of the global homogenization problem for generalized standard materials at finite strains, where a quasi-hyperelastic macro-stress potential is obtained from a global minimization problem with respect to the fine-scale displacement fluctuation field. It is shown that this global minimization problem determines the state of the micro-structure for finite increments of time. We consider three different settings of the global variational problem for prescribed displacements, non-trivial periodic displacements and prescribed stresses on the boundary of the micro-structure and develop numerical solution methods based on a spatial discretization of the fine-scale displacement fluctuation field. Representative applications of the proposed minimization principles are demonstrated for a constitutive model of crystal plasticity and the homogenization problem of texture analysis in polycrystalline aggregates.     

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     

On the use of homogeneous polynomials to develop anisotropic yield functions with applications to sheet forming

This paper investigates the capabilities of several non-quadratic polynomial yield functions to model the plastic anisotropy of orthotropic sheet metal (plane stress). Fourth, sixth and eighth-order homogeneous polynomials are considered. For the computation of the coefficients of the fourth-order polynomial an improved set of analytic formulas is proposed. For sixth and eighth-order polynomials the identification uses optimization. Simple constraints on the optimization process are shown to lead to real-valued convex functions. A general method to extend the above plane stress criteria to full 3D stress states is also suggested. Besides their simplicity in formulation, it is found that polynomial yield functions are capable to model a wide range of anisotropic plastic properties (e.g., the Numisheet’93 mild steel, AA2008-T4, AA2090-T3). The yield functions have then been implemented into a commercial finite element code as constitutive subroutines. The deep drawing of square (Numisheet’93) and cylindrical (AA2090-T3) cups have been simulated. In both cases excellent agreement with experimental data is obtained. In particular, it is shown that non-quadratic polynomial yield functions can simulate cylindrical cups with six or eight ears. We close with a discussion on earing and further examples.     

On the non-oscillation criterion for multiplicative anisotropic plasticity at large simple shear deformation

The criterion for non-oscillatory stresses under monotonic large simple shear deformation in the context of multiplicative anisotropic plasticity is discussed. In particular, evolving anisotropy combined with a Hill type of yield criterion is considered. It is shown that a sufficient, but not necessary, criterion for a non-oscillatory stress is ellipticity of the first Piola–Kirchhoff stress. Loss of ellipticity corresponds to a critical value h_cr of the generalized plastic modulus. Similarly, the absence of limit points on the stress–strain relation motivates an alternative criterion in terms of a critical value h_sh  h_cr. Finally, this criterion is demonstrated analytically as well as numerically for an important class of models with evolving anisotropy of the saturation type.     

On the evolution of lattice deformation in austenitic stainless steels—The role of work hardening at finite strains

In this work, a three dimensional crystal plasticity-based finite element model is presented to examine the micromechanical behaviour of austenitic stainless steels. The model accounts for realistic polycrystal micromorphology, the kinematics of crystallographic slip, lattice rotation, slip interaction (latent hardening) and geometric distortion at finite deformation. We utilise the model to predict the microscopic lattice strain evolution of austenitic stainless steels during uniaxial tension at ambient temperature with validation through in situ neutron diffraction measurements. Overall, the predicted lattice strains are in very good agreement with those measured in both longitudinal and transverse directions (parallel and perpendicular to the tensile loading axis, respectively). The information provided by the model suggests that the observed nonlinear response in the transverse {200} grain family is associated with a competitive bimodal evolution of strain during inelastic deformation. The results associated with latent hardening effects at the microscale also indicate that in situ neutron diffraction measurements in conjunction with macroscopic uniaxial tensile data may be used to calibrate crystal plasticity models for the prediction of the inelastic material deformation response.     

A nonlinear finite element model for the stress analysis of soft solids with a growing mass

The paper presents a new finite element (FE) model for the stress analysis of soft solids with a growing mass based on the work of Lubarda and Hoger (2002). Contrary to the traditional numerical methods emphasizing on the influence of growth on constitutive equations, an equivalent body force is firstly detected, which is resulted from the linearization of the nonlinear equation and acts as the driver for material growth in the numerical aspect. In the algorithm, only minor correction on the traditional tangent modulus is needed to take the growth effects into consideration and its objectivity could be guaranteed comparing with the traditional method. To solve the resulted equation in time domain, both explicit and implicit integration algorithms are developed, where the growth tensor is updated as an internal variable of Gauss point. The explicit updating scheme shows higher efficiency, while the implicit one seems to be more robust and accurate. The algorithm validation and its good performance are demonstrated by several two-dimensional examples, including free growth, constrained growth and stress dependent growth.