On the modeling of evolving anisotropy and large strains in pearlitic steel

A phenomenological model for deformation induced evolution of anisotropy at large strains in pearlitic steel is proposed. The modeled anisotropy is based on a homogenization of an ideal pearlitic microstructure. An areal affine type of reorientation is assumed for the individual grains. Furthermore, a yield criterion of the Hill type is proposed and motivated from the grain reorientation. In each pearlitic grain the cementite lamellas have a privileged direction. The symmetry group of each individual grain is therefore considered transversally isotropic. In a virgin material, the privileged directions of the different grains are randomly oriented, which allows for the interpretation that the material on the macroscopic length scale is initially isotropic. However, the cementite lamellas in the grains tend to align after large stretching or shearing deformation. The modeled evolution of anisotropy on the macroscopic length scale shows a saturation characteristics under large deformations.     

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     

Thermoviscoplastic modelling of asymmetric effects for polymers at large strains

Glassy polymers such as polycarbonate exhibit different behaviours in different loading scenarios, such as tension and compression. To this end a flow rule is postulated within a thermodynamic consistent framework in a mixed variant formulation and decomposed into a sum of weighted stress mode related quantities. The different stress modes are chosen such that they are accessible to individual examination in the laboratory, where tension and compression are typical examples. The characterisation of the stress modes is obtained in the octahedral plane of the deviatoric stress space in terms of the Lode angle, such that stress mode dependent scalar weighting functions can be constructed. Furthermore the numerical implementation of the constitutive equations into a finite element program is briefly described. In a numerical example, the model is used to simulate the laser transmission welding process.     

Mechanical and optical characterization of an anelastic polymer at large strain rates and large strains

Two constitutive relations have been determined from test results that characterize, respectively, the uniaxial and photomechanical behavior of a polyester-styrene copolymer for strain rates from 10^?5 to 3×10³ in./in./s and strains up to 40 percent. The high-strain-rate data were obtained by means of a split-Hopkinson-bar apparatus. Intermediate-strain-rate tests, performed with the aid of a drop tower, were reported in an earlier paper. Quasi-static experiments were conducted on a standard testing machine. A nonlinear, four-parameter, elastic-viscoplastic model was constructed which describes the mechanical behavior. The parameters were determined by a least-mean-squares curve-fitting procedure. The viscoplastic parameters were found to obey a power law in strain rate. The photomechanical model was found to be linear with strain well into the plastic-deformation region, while the slope of the strain-birefringence curve for each strain rate also varied by strain rate to a power.     

Mathematical modeling of mechanical properties of metals and alloys at large strains

We discuss problems in mathematical modeling of the mechanical behavior of metals and alloys at large strains. Attention is mainly paid to the analysis of the stress-strain state of specimens and structural fragments made of highly plastic materials with the effect of stability loss under tensile stresses taken into account. We discuss the methods for determining the true property diagram at strains exceeding the ultimate uniform strain. We process experimental data and determine the true property diagrams for AMg6, AMg6M, and 1201 aluminum alloys and BrKh08 alloy. To calculate the load-carrying capacity of structural members, one often uses the conventional ultimate strength σ _b accepted in regulations as a material characteristic. But it follows from the method for experimentally determining this characteristic that it depends on the properties of the specimen viewed as a structure. As a result, a formal use of fracture criteria recommended in regulations leads to a discrepancy between design and experimental values of fracture loads. Nowadays, the finite element method is widely used in practical strength analysis. This method permits one to study the elastoplastic strained state of geometrically complicated structures in detail, take into account physical nonlinearity at large strains, determine damage boundaries, and improve experimental methodology. The wide capabilities of this method allow one to use test results more completely.     

An analytical micro-macro model for textured polycrystals at large plastic strains

An analytical micro-macro model of evolving plastic anisotropy is presented that is suitable for numerical simulation of forming processes. The model is based on the combination of a polycrystal model and different analytical procedures for writing anisotropic plastic potentials, expressing their coefficients in terms of texture coefficients, and updating the texture coefficients as function of the (tensorial) strain increment. The use of a fourth-order dual plastic potential (“C4”) in the analytical micro-macro model is studied, and this use is compared with that of Hill's [1948] yield criterion and also with the usual run of the Taylor model. The coefficients of the C4 potential depend linearly on the texture coefficients, which are updated using a variational polycrystal model. The analytical operation of this updating lies on the method first proposed by Eslinget al. [1984] and is described and checked in some detail. The predictions of the analytical micro-model compare well with measurements of the Lankford coefficient, provided the C4 potential is used. The predicted texture evolution is also in a good experimental agreement: a better one than with the Taylor model, which in some cases, gives a poor updating. The theoretical stress evolution during biaxial or plane-strain tension is experimentally consistent too, although in that case the C4 potential, closer to Taylor's model, makes no improvement as compared with Hill's quadratic criterion.     

Plane strain pure bending of sheets with damage evolution at large strains

An analysis of plane-strain bending at large strains for the rigid/plastic incompressible material model including arbitrary strain-hardening and damage evolution laws is performed. The fracture criterion is based on a critical value of the damage parameter. Numerical treatment is reduced to the system of two partial differential equations written in characteristic coordinates. The through-thickness distribution of the principal stresses and damage parameter as well as the variation of the bending moment with the radius of curvature of the concave surface are found for Swift’s hardening law and one specific damage evolution law. General tendencies in solution behaviour are in agreement with physical expectations.     

A thermodynamics framework and numerical aspects for transformation-induced plasticity at large strains

In this work, we present a macroscopic material model for simulation transformation-induced plasticity, which is an important phenomenon in metal forming processes. The model is formulated within a thermodynamic framework at large strains. In order to account for both, phase transformation and plasticity, yield functions are related to these effects. Then, applying the concept of maximum dissipation evolution equations are obtained for the inelastic strains, the transformation strains, a hardening variable and the volume fraction of martensite. Furthermore the numerical implementation of the constitutive equations into a finite element program is described. In a numerical example we investigate the austenite-to-martensite phase transformation in a shaft subjected to thermo-mechanical loading in a hybrid-forming process.     

Damage evolution in an expanding/contracting hollow sphere at large strains

A generalization of one of the classical problems of plasticity theory, expansion/contraction of a hollow sphere, is proposed assuming that the conventional constitutive equations for rigid plastic, hardening material are supplemented with an arbitrary ductile damage evolution law. No restriction is imposed on the hardening law in the analytic part of the solution. The initial/boundary value problem is reduced to two equations in characteristic coordinates. A numerical scheme to solve these equations is proposed. An illustrative example is given.     

The effect of anisotropy on the integral representation of a cylindrical pulse

Summary The infinite medium Green's function for a two dimensional anisotropic scalar wave equation is obtained in closed form using a technique developed by De Hoop¹). The effect of anisotropy on the complex contour integral representation of this Green's function is explicitly exhibited.Publication 367, Institute of Geophysics and Planetary Physics, University of California, Los Angeles.     

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.     

A theoretical and computational framework for anisotropic continuum damage mechanics at large strains

The main objective of this work is the formulation and algorithmic treatment of anisotropic continuum damage mechanics at large strains. Based on the concept of a fictitious, isotropic, undamaged configuration an additional linear tangent map is introduced which allows the interpretation as a damage deformation gradient. Then, the corresponding Finger tensor – denoted as damage metric – constructs a second order, internal variable. Due to the principle of strain energy equivalence with respect to the fictitious, effective space and the standard reference configuration, the free energy function can be computed via push-forward operations within the nominal setting. Referring to the framework of standard dissipative materials, associated evolution equations are constructed which substantially affect the anisotropic nature of the damage formulation. The numerical integration of these ordinary differential equations is highlighted whereby two different schemes and higher order methods are taken into account. Finally, some numerical examples demonstrate the applicability of the proposed framework.     

Thermodynamically consistent phase field approach to dislocation evolution at small and large strains

A thermodynamically consistent, large strain phase field approach to dislocation nucleation and evolution at the nanoscale is developed. Each dislocation is defined by an order parameter, which determines the magnitude of the Burgers vector for the given slip planes and directions. The kinematics is based on the multiplicative decomposition of the deformation gradient into elastic and plastic contributions. The relationship between the rates of the plastic deformation gradient and the order parameters is consistent with phenomenological crystal plasticity. Thermodynamic and stability conditions for homogeneous states are formulated and satisfied by the proper choice of the Helmholtz free energy and the order parameter dependence on the Burgers vector. They allow us to reproduce desired lattice instability conditions and a stress-order parameter curve, as well as to obtain a stress-independent equilibrium Burgers vector and to avoid artificial dissipation during elastic deformation. The Ginzburg–Landau equations are obtained as the linear kinetic relations between the rate of change of the order parameters and the conjugate thermodynamic driving forces. A crystalline energy coefficient for dislocations is defined as a periodic step-wise function of the coordinate along the normal to the slip plane, which provides an energy barrier normal to the slip plane and determines the desired, mesh-independent height of the dislocation bands for any slip system orientation. Gradient energy contains an additional term, which excludes the localization of a dislocation within a height smaller than the prescribed height, but it does not produce artificial interface energy. An additional energy term is introduced that penalizes the interaction of different dislocations at the same point. Non-periodic boundary conditions for dislocations are introduced which include the change of the surface energy due to the exit of dislocations from the crystal. Obtained kinematics, thermodynamics, and kinetics of dislocations at large strains are simplified for small strains and rotations, as well.     

A finite element model for shape memory alloys considering thermomechanical couplings at large strains

In the present work we propose a new thermomechanically coupled material model for shape memory alloys (SMA) which describes two important phenomena typical for the material behaviour of shape memory alloys: pseudoelasticity as well as the shape memory effect. The constitutive equations are derived in the framework of large strains since the martensitic phase transformation involves inelastic deformations up to 8%, or even up to 20% if the plastic deformation after the phase transformation is taken into account. Therefore, we apply a multiplicative split of the deformation gradient into elastic and inelastic parts, the latter concerning the martensitic phase transformation. An extended phase transformation function has been considered to include the tension–compression asymmetry particularly typical for textured SMA samples. In order to apply the concept in the simulation of complex structures, it is implemented into a finite element code. This implementation is based on an innovative integration scheme for the existing evolution equations and a monolithic solution algorithm for the coupled mechanical and thermal fields. The coupling effect is accurately investigated in several numerical examples including pseudoelasticity as well as the free and the suppressed shape memory effect. Finally, the model is used to simulate the shape memory effect in a medical foot staple which interacts with a bone segment.     

Computational modeling of inelastic large ratcheting strains

A framework for phenomenological hyperelasto-plasticity with initial anisotropy, kinematic hardening as well as anisotropic damage is presented in [Menzel et al., Int. J. Plasticity (2004), in press]. In this contribution, we exploit and extend this framework to include several back-stresses in order to capture the ratcheting response of polycrystalline metals subjected to cyclic stress with non-zero mid-value. The evolution equations for kinematic hardening resemble a linear combination of the multiple-Armstrong–Frederick and the Burlet–Cailletaud models, which are extended to the large strain setting. The capability of the model to capture various phenomenological characteristics, in particular multi-axial ratcheting, is illustrated by numerical examples. Comparisons with uni-axial and bi-axial experimental ratcheting results for carbon steel are given. Finally, the finite element analysis of a simplified railway turnout component subjected to cyclic loading is presented.     

A constitutive model for cold deformation of aluminium at large strains and high strain rates

A temperature-dependent constitutive model for viscoplastic deformation of aluminium based on a single, scalar internal variable is presented. The model is designed particularly for the strains, strain rates, and temperatures important for cold forging. Special attention is paid to the underlying physical processes that determine the flow stress in the metal. The kinetic constitutive equation is based on thermal activation of dislocations over an average potential barrier from various kinds of obstacles. Strain hardening is modelled through the internal variable which represents the increasing height of these barriers. The model is generalized to three dimensions, and it has been implemented in the finite-element code ABAQUS. Simulations of simple forging operations are presented.     

Evaluating large plastic strains in shock-loaded beams

Conclusions The simple relations obtained make it possible to satisfactorily evaluate large plastic strains (deflections) of uniform and nonuniform beams with fixed ends in the case where the beam is made of a strain-rate-sensitive material and is subjected to static or purely impulsive transverse loads. With allowance for the method in [2], these relations can be used to determine combinations of the parameters P and I of different forms of shock load corresponding to a given level of large plastic strains in beams.Moscow. Translated from Prikladnaya Mekhanika, Vol. 22, No. 3, pp. 66–71, March, 1986.