Unit cell debonding analyses for arbitrary orientation of plastic anisotropy

Debonding of rigid inclusions embedded in the elastic–plastic aluminum alloy Al 2090-T3 is analyzed numerically using a unit cell model taking full account of finite strains. The cell is subjected to overall biaxial plane strain tension and periodical boundary conditions are applied to represent arbitrary orientations of plastic anisotropy. Plastic anisotropy is considered using two phenomenological anisotropic yield criteria, namely Hill [Proceedings of the Royal Society of London A 193 (1948) 281] and Barlat et al. [International Journal of Plasticity 7 (1991) 693]. For this material plastic anisotropy delays debonding compared to plastic isotropy except for the case of Hill’s yield function when the tensile directions coincided with the principal axes of anisotropy. For some inclinations of the principal axes of anisotropy relative to the tensile directions, the stress strain responses are identical but the deformation modes are mirror images of each other.     

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.     

A hardening model for anisotropic materials

Results from proportionate-loading experiments under plane-stress states illustrate the existence of a field of uniform-hardening potentials for the yield and creep deformation behavior of isotropic, slightly anisotropic, aelotropic and orthotropic polycrystalline materials in the initially strain-free condition. For two different plane-stress states, it is shown that a linear functional relationship holds between the plastic-strain increment ratio and the stress ratio in these materials and that, consequently, the field is adequately modeled by a uniform-hardening anisotropic function that is quadratic in the components of deviatoric stress. Anisotropic plane-stress yield functions are formulated for any stage in the deformation process by combining the uniform-hardening function with the kinematic-hardening rule. The resulting functions, which correspond to rigid translations of initial yield loci according to Ziegler's rule, provide good agreement with experimental observations on a marked Bauschinger effect and an absence of cross hardening.     

A plane stress anisotropic plastic flow theory for orthotropic sheet metals

A phenomenological theory is presented for describing the anisotropic plastic flow of orthotropic polycrystalline aluminum sheet metals under plane stress. The theory uses a stress exponent, a rate-dependent effective flow strength function, and five anisotropic material functions to specify a flow potential, an associated flow rule of plastic strain rates, a flow rule of plastic spin, and an evolution law of isotropic hardening of a sheet metal. Each of the five anisotropic material functions may be represented by a truncated Fourier series based on the orthotropic symmetry of the sheet metal and their Fourier coefficients can be determined using experimental data obtained from uniaxial tension and equal biaxial tension tests. Depending on the number of uniaxial tension tests conducted, three models with various degrees of planar anisotropy are constructed based on the proposed plasticity theory for power-law strain hardening sheet metals. These models are applied successfully to describe the anisotropic plastic flow behavior of 10 commercial aluminum alloy sheet metals reported in the literature.     

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     

Finite strain analysis of simple shear using recent anisotropic yield criteria

The finite strain response of a rectangular block subjected to constrained simple shearing deformations is considered in order to evaluate the predictive capability of some recently proposed anisotropic yield functions. It is shown that, in the presence of plane anisotropy, the prediction of realistic second order normal stresses cannot be expected since every different initial orientation of the material axes relative to the loading axes results in a different response due to the rotation of the material axes with shear. Parametric studies are performed in order to determine possible limits on the material constants so that the predicted normal stresses remain second order with respect to the shear stress itself. Our numerical results indicate that, in particular, the commonly employed range of one parameter associated with grain related anisotropy renders results which no longer imply predicting a proper second order effect but rather introduces errors of the first order. The results suggest that the modelling of anisotropy with such phenomenological anisotropic yield functions should be limited to near-quadratic yield surfaces for applications involving stress states outside the biaxial tensile stress range.     

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.     

Failure prediction in anisotropic sheet metals under forming operations with consideration of rotating principal stretch directions

An approximate macroscopic yield criterion for anisotropic porous sheet metals is adopted to develop a failure prediction methodology that can be used to investigate the failure of sheet metals under forming operations. Hill's quadratic anisotropic yield criterion is used to describe the matrix normal anisotropy and planar isotropy. The approximate macroscopic anisotropic yield criterion is a function of the anisotropy parameter R, defined as the ratio of the transverse plastic strain rate to the through-thickness plastic strain rate under in-plane uniaxial loading conditions. The Marciniak–Kuczynski approach is employed here to predict failure/plastic localization by assuming a slightly higher void volume fraction inside randomly oriented imperfection bands in a material element of interest. The effects of the anisotropy parameter R, the material/geometric inhomogeneities, and the potential surface curvature on failure/plastic localization are first investigated. Then, a non-proportional deformation history including relative rotation of principal stretch directions is identified in a critical element of a mild steel sheet under a fender forming operation given as a benchmark problem in the 1993 NUMISHEET conference. Based on the failure prediction methodology, the failure of the critical sheet element is investigated under the non-proportional deformation history. The results show that the gradual rotation of principal stretch directions lowers the failure strains of the critical element under the given non-proportional deformation history.     

Crack growth resistance for anisotropic plasticity with non-normality effects

For a plastically anisotropic solid a plasticity model using a plastic flow rule with non-normality is applied to predict crack growth. The fracture process is modelled in terms of a traction–separation law specified on the crack plane. A phenomenological elastic–viscoplastic material model is applied, using one of two different anisotropic yield criteria to account for the plastic anisotropy, and in each case the effect of the normality flow rule is compared with the effect of non-normality. Conditions of small scale yielding are assumed, with mode I loading conditions far from the crack-tip, and various directions of the crack plane relative to the principal axes of the anisotropy are considered. It is found that the steady-state fracture toughness is significantly reduced when the non-normality flow rule is used. Furthermore, it is shown that the predictions are quite sensitive to the value of the maximum angle of deviation from normality in the non-normality flow rule.     

Anisotropic parameter identification using inhomogeneous tensile test

In this contribution, an inverse identification strategy of constitutive laws for elastoplastic behaviour is presented. The proposed inverse algorithm is composed on an appropriate finite element calculation combined with an optimisation procedure. It is applied to identify material anisotropic coefficients using a set up of easy performed laboratory tests. The used experimental data are the plane tensile test and the off axes tensile tests. The identified behaviour models are mainly based on Hill's quadratic yield criterion. Two cases of this yield criterion have been considered: the transverse isotropic and the orthotropic one under an associated and non-associated flow rule assumptions for each case. The yield surface has been assumed to expand isotropically (isotropic strain hardening law) as a function of the plastic work.In order to better describe anisotropic plastic properties of the studied materials, a recently planar anisotropic yield function is used. It is a non-quadratic yield criterion which takes account of anisotropic yield stresses as well as anisotropic strain ratios. It is subsequently shown that the agreement between inverse identification results and experimental measurements were improved.We prove also that the presented strategy is a good alternative to the simplified homogeneous tests assumption, especially for the plane tensile test.     

Theoretical considerations upon the MK model for limit strains prediction: The plane strain case,strain-rate effects,yield surface influence,and material heterogeneity

The paper presents a study of the Marciniak and Kuczynski (MK for short) model for the prediction of limit strains of orthotropic sheet metal under in-plane proportional biaxial stretching. In two particular cases analytical results can be obtained if the groove of the MK model is oriented along one of the in-plane symmetry axes. The first case is the plane strain loading mode. Necessary and sufficient conditions are derived for the MK-predicted plane strain limit strain to match exactly the experimentally measured limit strain. An example of material, the AA5182-O aluminum alloy, that does not satisfy these conditions is discussed. It is shown then that if a power-law strain rate sensitivity is included in the hardening law then the MK-model can match exactly any target plane strain limit strain. The second case is the non-hardening case for positive strain ratios. This case allows for an insight into the way the MK-predicted limit strains depend upon the yield function. Based on the theory developed for the plane strain case, material heterogeneity as a possible cause for unstable plastic flow is further discussed. It is shown that such heterogeneities can be modeled by perturbing the rate of deformation with an eigenstrain. This allows for an extension of the MK-model to sheets of uniform thickness.     

Incorporating the effects of fabric in the dilatant double shearing model for planar deformation of granular materials

The objective of this paper is to incorporate the effects of fabric and its evolution into the Dilatant Double Shearing Model [Mehrabadi, M.M., Cowin, S.C., 1978. Initial planar deformation of dilatant granular materials. J. Mech. Phys. Solids 26, 269–284] for granular materials in order to capture the anisotropic behavior and the complex response of granular materials in cyclic shear loading. An important consequence of considering the fabric is that one can have unequal shearing rates along the two slip directions. This property leads to the non-coaxiality of the principal axes of stress and strain rate, which is more appropriate for a material that exhibits initial and induced anisotropy. In addition, we employ a fabric-dependent elasticity tensor with orthotropic symmetry. The model developed in this paper also predicts one of the experimentally observed characteristics of granular materials: the gradual concentration of the contact normals towards the maximum principal stress direction.We implement the constitutive equations into ABAQUS/Explicit by writing a user material subroutine in order to predict the strength anisotropy of granular materials in a plane strain biaxial compression test and investigate the mechanical behavior of granular materials under the cyclic shear loading conditions. The predictions from this model show good quantitative agreement with the experiments of [Park, C.S., 1990. Anisotropy in deformation and strength properties of sands in plane strain compression, Masters Thesis, University of Tokyo; Park, C.S., Tatsuoka, F., 1994. Anisotropic strength and deformation of sands in plane strain compression. In: XIII ICSMFE, New Delhi, India; Okada, N., 1992. Energy dissipation in inelastic flow of cohesionless granular media. Ph.D. Thesis, University of California, San Diego].     

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.     

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.     

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.