Micromechanical modelling of the effect of plastic deformation on the mechanical behaviour in pseudoelastic shape memory alloys

Except for the recoverable strain induced by phase transformation, NiTi alloys are very ductile even in the martensite phase. The purpose of the present paper is to study the influence of permanent deformation, which results from plastic deformation of martensite, on the mechanical behaviour of pseudoelastic NiTi alloys. Based on phenomenological theory of martensitic transformation and crystal plasticity, a new three dimensional micromechanical model is proposed by coupling both the slip and twinning deformation mechanisms. The present model is implemented as User MATerial subroutine (UMAT) into ABAQUS/Standard to study the influences of plastic deformation on the stress and strain fields, and on the evolution of martensite transformation. Results show that with the increasing of plastic deformation the residual strain increases and the phase transformation stress–strain curves from the martensite to austenite become steeper and less obvious. Both characteristics, stabilisation of martensite and impedance of the reverse transformation, due to plastic deformation are captured.     

Coupling between martensitic phase transformations and plasticity: A microstructure-based finite element model

A microstructural finite element (MFE) model is developed to capture the interaction between martensitic transformations and plasticity in NiTi shape memory alloys (SMAs). The interaction is modeled through the grain-to-grain redistribution of stress caused by both plasticity and phase transformation, so that each mechanism affects the driving force of the other. A unique feature is that both processes are modeled at a crystallographic level and are allowed to operate simultaneously. The model is calibrated to pseudoelastic data for select single crystals of Ti–50.9at.%Ni. For polycrystals, plasticity is predicted to enhance the overall martensite volume fraction at a given applied stress. Upon unloading, residual stress can induce remnant (retained) martensite. For thermal cycling under load bias, plasticity is observed to limit the net transformation strain/cycle and increase the hysteretic width. Deformation processing, via plastic pre-straining at elevated temperature, is shown to dramatically alter subsequent pseudoelastic response, as well as induce two-way shape memory behavior during no-load thermal cycling. Overall, the model is suitable at smaller imposed strains, where martensite detwinning is not expected to dominate.     

Grain size effects in multiphase steels assisted by transformation-induced plasticity

The influence of the austenitic grain size on the overall stress–strain behavior in a multiphase carbon steel is analyzed through three-dimensional finite element simulations. A recently developed multiscale martensitic transformation model is combined with a plasticity model to simulate the transformation-induced plasticity effects of a grain of retained austenite embedded in a ferrite-based matrix. Grain size effects are included via a surface energy term in the Helmholtz energy. Tensile simulations for representative orientations of the grain of retained austenite show that the initial stability of the austenite increases as the grain size decreases. Consequently, the effective strength is initially higher for smaller grains. The influence of the grain size on the evolution of the transformation process strongly depends on the grain orientation. For “hard” orientations, the transformation rate is higher for larger grains. In addition, the phase transformation is partially suppressed as the grain size decreases. In contrast, for “soft” orientations, the transformation rate is lower for larger grains. The phase transformation is more homogeneous for smaller grains and, consequently, the effective transformation strain is larger. Nevertheless, in multiphase carbon steels with a relatively low percentage of retained austenite, the influence of the austenitic grain size on the overall constitutive response is smaller than the influence of the austenitic grain orientation.     

Computational simulation of the dependence of the austenitic grain size on the deformation behavior of TRIP steels

Due to the strain-induced martensitic transformation which occurs during plastic deformation, a transformation-induced plasticity (TRIP) phenomenon is generated. With the TRIP phenomenon, the TRIP steel possesses favorable mechanical properties such as high strength, ductility and toughness, and is frequently employed as a structural material. In the past, several researchers clarified experimentally that the strain-induced martensitic transformation and the deformation behavior of TRIP steel depend upon the austenitic grain size. In order to obtain the expected mechanical properties of TRIP steel through control of the austenitic grain size, prediction and control of the material characteristics in the deformation processes is essential. Here, the new strain-induced martensitic transformation kinetics model and constitutive equation of TRIP steels are proposed by considering the dependence of the austenitic grain size. Then, the deformation behavior of a type 304 austenitic stainless steel cylinder is simulated under different environmental temperatures with the various austenitic grain sizes by the finite-element method along with newly-proposed constitutive equations. Finally, the validity of proposed constitutive equations and the possibility of the improvement of the mechanical properties through control of the austenitic grain size are discussed.     

Couplings between plasticity and martensitic phase transformation: overall behavior of polycrystalline TRIP steels

Several couplings between plasticity and martensitic phase transformation are at the origin of remarkable properties of ductility and toughness in the case of TRIP steels. A micromechanical model is developed to predict the conditions of nucleation and growth of a martensitic microdomain inside an inhomogeneous plastic strain field. More explicit relations are developed in the case of a simple shear test where a heterogeneous plastic strain field leads to a significant decrease of the critical stress for martensitic transformation. The obtained results are combined with a kinetics and kinematics studies to derive the constitutive equation of an austenitic single crystal from which the overall behavior of a polycrystalline steel is deduced using the self-consistent scale transition method. Comparison with experimental data shows a good agreement.     

Micromechanics of plastic deformation and phase transformation in a three-phase TRIP-assisted advanced high strength steel: Experiments and modeling

The micromechanics of plastic deformation and phase transformation in a three-phase advanced high strength steel are analyzed both experimentally and by microstructure-based simulations. The steel examined is a three-phase (ferrite, martensite and retained austenite) quenched and partitioned sheet steel with a tensile strength of ~980 MPa. The macroscopic flow behavior and the volume fraction of martensite resulting from the austenite–martensite transformation during deformation were measured. In addition, micropillar compression specimens were extracted from the individual ferrite grains and the martensite particles, and using a flat-punch nanoindenter, stress–strain curves were obtained. Finite element simulations idealize the microstructure as a composite that contains ferrite, martensite and retained austenite. All three phases are discretely modeled using appropriate crystal plasticity based constitutive relations. Material parameters for ferrite and martensite are determined by fitting numerical predictions to the micropillar data. The constitutive relation for retained austenite takes into account contributions to the strain rate from the austenite–martensite transformation, as well as slip in both the untransformed austenite and product martensite. Parameters for the retained austenite are then determined by fitting the predicted flow stress and transformed austenite volume fraction in a 3D microstructure to experimental measurements. Simulations are used to probe the role of the retained austenite in controlling the strain hardening behavior as well as internal stress and strain distributions in the microstructure.     

Characterization of yielding behavior of polycrystalline metals with single crystal plasticity based on representative characteristic length

Single crystal plasticity based on a representative characteristic length is proposed and introduced into a homogenization approach based on finite element analyses, which are applied to characterization of distinctive yielding behaviors of polycrystalline metals, yield-point elongation, and grain size strengthening. The computational manner for an implicit stress update is derived with the framework of a standard multi-surface plasticity at finite strain, where the evolution of the characteristic lengths are numerically converted from the accumulated slips of all of slip systems by exploiting the mathematical feature of the characteristic length as the intermediate function of the plastic internal variables. Furthermore, a constitutive model for a single crystal reproduces the stress–strain curve divided into three parts. Using two-scale finite element analysis, the macroscopic stress–strain response with yield-point elongation under a situation of low dislocation density is reproduced. Finally, the grain size effect on the yield strength is analyzed with modeling of the grain boundary in the context of the proposed constitutive model and is discussed from both macroscopic and microscopic views.     

基于有限元模拟的SUS304奥氏体不锈钢摩擦耦合变形过程的应力分析

采用ANSYS/LS-DYNA软件对SUS 304奥氏体不锈钢薄板的摩擦耦合变形过程进行了数值模拟.采用隐式-显式序列求解法和分段线性塑性材料模型,分析了钢带摩擦耦合变形时的应力分布规律及载荷、下压量和滑动速度等因素对钢带剪应力、主应力及等效应力的影响.结果表明:摩擦耦合变形的试验参数显著影响钢带的应力分布,验证了钢带在低于其屈服强度的应力条件下发生塑性变形的摩擦诱发效应.对奥氏体不锈钢摩擦诱发马氏体转变行为的研究及其摩擦学性能的改善具有一定的指导意义.     

Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

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.     

On the strain hardening and texture evolution in high manganese steels: Experiments and numerical investigation

We present a systematic investigation on the strain hardening and texture evolution in high manganese steels where twinning induced plasticity (TWIP) plays a significant role for the materials' plastic deformation. Motivated by the stress–strain behavior of typical TWIP steels with compositions of Fe, Mn, and C, we develop a mechanistic model to explain the strain-hardening in crystals where deformation twinning dominates the plastic deformation. The classical single crystal plasticity model accounting for both dislocation slip and deformation twinning are then employed to simulate the plastic deformation in polycrystalline TWIP steels. While only deformation twinning is activated for plasticity, the simulations with samples composed of voronoi grains cannot fully capture the texture evolution of the TWIP steel. By including both twinning deformation and dislocation slip, the model is able to capture both the stress–strain behaviors and the texture evolution in Fe–Mn–C TWIP steel in different boundary-value problems. Further analysis on the strain contributions by both mechanisms suggests that deformation twinning plays the dominant role at the initial stage of plasticity in TWIP steels, and dislocation slip becomes increasingly important at large strains.     

Simulation of hydroforming of steel tube made of metastable stainless steel

The Olson–Cohen model for strain-induced deformation, further developed by Stringfellow and others, has been calibrated together with a flow stress model for the plastic deformation of metastable stainless steel. Special validation tests for checking one of the limitations of the model have also been carried out. The model has been implemented into a commercial finite element code using a staggered approach for integrating the stress–strain relations with the microstructure model. Results from a thermo-mechanical coupled simulation of hydroforming of a tube have been compared with corresponding experiments. The agreement between experimental results of radial expansion and martensite fraction and the corresponding computed results is good.     

The role of intergranular constraint on the stress-induced martensitic transformation in textured polycrystalline NiTi

The influences of grain boundaries and relative grain misorientations on stress-induced martensitic transformations in NiTi are studied using unique experiments and finite element modeling. Tensile and compressive mechanical tests reveal that polycrystalline NiTi with a dominant <111> fiber texture and single crystal NiTi oriented along the [111] direction exhibit nearly identical stress–strain curves during a stress-induced martensitic transformation. Micro-mechanical finite element simulations of fiber textured polycrystals and single crystals undergoing a multi-variant martensitic transformation confirm the relative indifference of the macroscopic transformation attributes to the presence of grain boundaries. On the microscale, the finite element simulations further reveal that the insensitivity of the transformation to intergranular constraint is linked to the local stress disturbance created by transforming grains. The transformation of grains that are favorably oriented with respect to the loading axis creates local stresses that invariably assist the transformation in neighboring grains, effectively lowering the influence of grain misorientations and boundaries on the macroscopic transformation behavior.     

Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation

A strain gradient dependent crystal plasticity approach is used to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. Material points are considered as aggregates of grains, subdivided into several fictitious grain fractions: a single crystal volume element stands for the grain interior whereas grain boundaries are represented by bi-crystal volume elements, each having the crystallographic lattice orientations of its adjacent crystals. A relaxed Taylor-like interaction law is used for the transition from the local to the global scale. It is relaxed with respect to the bi-crystals, providing compatibility and stress equilibrium at their internal interface. During loading, the bi-crystal boundaries deform dissimilar to the associated grain interior. Arising from this heterogeneity, a geometrically necessary dislocation (GND) density can be computed, which is required to restore compatibility of the crystallographic lattice. This effect provides a physically based method to account for the additional hardening as introduced by the GNDs, the magnitude of which is related to the grain size. Hence, a scale-dependent response is obtained, for which the numerical simulations predict a mechanical behaviour corresponding to the Hall-Petch effect. Compared to a full-scale finite element model reported in the literature, the present polycrystalline crystal plasticity model is of equal quality yet much more efficient from a computational point of view for simulating uniaxial tension experiments with various grain sizes.     

Anisotropic plasticity model coupled with Lode angle dependent strain-induced transformation kinetics law

A phenomenological macroscopic plasticity model is developed for steels that exhibit strain-induced austenite-to-martensite transformation. The model makes use of a stress-state dependent transformation kinetics law that accounts for both the effects of the stress triaxiality and the Lode angle on the rate of transformation. The macroscopic strain hardening is due to nonlinear kinematic hardening as well as isotropic hardening. The latter contribution is assumed to depend on the dislocation density as well as the current martensite volume fraction. The constitutive equations are embedded in the framework of finite strain isothermal rate-independent anisotropic plasticity. Experimental data for an anisotropic austenitic stainless steel 301LN is presented for uniaxial tension, uniaxial compression, transverse plane strain tension and pure shear. The model parameters are identified using a combined analytical–numerical approach. Numerical simulations are performed of all calibration experiments and excellent agreement is observed. Moreover, we make use of experimental data from ten combined tension and shear experiments to validate the proposed constitutive model. In addition, punch and notched tension tests are performed to evaluate the model performance in structural applications with heterogeneous stress and strain fields.     

Multi-scale finite element analyses of sheet metals by using SEM-EBSD measured crystallographic RVE models

In this study, two multi-scale analyses codes are newly developed by combining a homogenization algorithm and an elastic/crystalline viscoplastic finite element (FE) method (Nakamachi, E., 1988. A finite element simulation of the sheet metal forming process. Int. J. Numer. Meth. Eng. 25, 283–292; Nakamachi, E., Dong, X., 1996. Elastic/crystalline viscoplastic finite element analysis of dynamic deformation of sheet metal. Int. J. Computer-Aided Eng. Software 13, 308–326; Nakamachi, E., Dong, X., 1997. Study of texture effect on sheet failure in a limit dome height test by using elastic/crystalline viscoplastic finite element analysis. J. Appl. Mech. Trans. ASME(E) 64, 519–524; Nakamachi, E., 1998. Elastic/crystalline viscoplastic finite element modeling based on hardening–softening evaluation equation. In: Proc. of the 6th NUMIFORM, pp. 315–321; Nakamachi, E., Hiraiwa, K., Morimoto, H., Harimoto, M., 2000a. Elastic/crystalline viscoplastic finite element analyses of single- and poly-crystal sheet deformations and their experimental verification. Int. J. Plasticity 16, 1419–1441; Nakamachi, E., Xie, C.L., Harimoto, M., 2000b. Drawability assessment of BCC steel sheet by using elastic/crystalline viscoplastic finite element analyses. Int. J. Mech. Sci. 43, 631–652); (1) a “semi-implicit” finite element (FE) code and (2) a “dynamic explicit” FE code. These were applied to predict the plastic strain induced yield loci and the formability of sheet metal in the macro scale, and simultaneously the crystal texture and hardening evolutions in the micro scale. The isotropic and kinematical hardening laws are employed in the crystalline plasticity constitutive equation. For the multi-scale structure, two-scales are considered. One is a microscopic polycrystal structure and the other a macroscopic elastic plastic continuum. We measure crystal morphologies by using the SEM-EBSD apparatus with a unit of about 3.8 μm voxel, and define a three dimensional (3D) representative volume element (RVE) for the micro polycrystal structure, which satisfy the periodicity condition of crystal orientation distribution. A “micro” finite element modeling technique is newly established to minimize the total number of finite elements in the micro scale. Next, the “semi-implicit” crystallographic homogenization FE code, which employs the SEM-EBSD measured RVE, is applied to the 99.9% pure-iron uni-axial tensile problem to predict the texture evolution and the subsequent yield loci in the various strain paths. These “semi implicit” results reveal that the plastic strain induced anisotropy in the micro and macro levels can be predicted by our FE analyses. The kinematical hardening law leads a distinct plastic strain induced anisotropy. Our “dynamic-explicit” FE code is applied to simulate the limit dome height (LDH) test problem of the mild steel DQSK, the high strength steel HSLA and the aluminum alloy AL6022 sheet metals, which were adopted as the NUMISHEET2005 Benchmark sheet metals (Smith, L.M., Pourboghrat, F., Yoon, J.-W., Stoughton, T.B., 2005. NUMISHEET2005. In: Proc. of 6th Int. Conf. Numerical Simulation of 3D Sheet Metal Forming Processes, PART A and B(Benchmark), pp. 409–451) to estimate formability. The “dynamic explicit” results reveal that the initial crystal orientation distribution has a large affects to a plastic strain induced texture and anisotropic hardening evolutions and sheet formability.     

Microscale prediction of deformation in an austenitic stainless steel under uniaxial loading

In this study, the deformation behaviour of polycrystalline austenitic 316H stainless steel under uniaxial loading is investigated by means of in-situ neutron diffraction (ND) measurement and crystal plasticity-based finite element (FE) modelling. Data have been obtained for the macroscopic stress–strain response and the lattice strain evolution in the longitudinal and transverse direction relative to the uniaxial loading axis. Comparison between the model predictions and the ND measurements suggests that in most cases the FE model can predict the lattice strain evolution at the microscale and capture the general trends observed in the experiments. Both ND measurements and FE modelling simulations identify little effect of micromorphology effect on the longitudinal lattice strain evolution, while the transverse lattice strain response appears to be sensitive to the microstructure, in particular the initial crystallographic orientation of the material.     

Forming simulation of aluminum sheets using an anisotropic yield function coupled with crystal plasticity theory

This paper describes the application of a coupled crystal plasticity based microstructural model with an anisotropic yield criterion to compute a 3D yield surface of a textured aluminum sheet (continuous cast AA5754 aluminum sheet). Both the in-plane and out-of-plane deformation characteristics of the sheet material have been generated from the measured initial texture and the uniaxial tensile curve along the rolling direction of the sheet by employing a rate-dependent crystal plasticity model. It is shown that the stress–strain curves and R-value distribution in all orientations of the sheet surface can be modeled accurately by crystal plasticity if a “finite element per grain” unit cell model is used that accounts for non-uniform deformation as well as grain interactions. In particular, the polycrystal calculation using the Bassani and Wu (1991) single crystal hardening law and experimental electron backscatter data as input has been shown to be accurate enough to substitute experimental data by crystal plasticity data for calibration of macroscopic yield functions. The macroscopic anisotropic yield criterion CPB06ex2 (Plunkett et al., 2008) has been calibrated using the results of the polycrystal calculations and the experimental data from mechanical tests. The coupled model is validated by comparing its predictions with the anisotropy in the experimental yield stress ratio and strain ratios at 15% tensile deformation. The biaxial section of the 3D yield surface calculated directly by crystal plasticity model and that predicted by the phenomenological model calibrated with experimental and crystal plasticity data are also compared. The good agreement shows the strength of the approach. Although in this paper, the Plunkett et al. (2008) yield function is used, the proposed methodology is general and can be applied to any yield function. The results presented here represent a robust demonstration of implementing microscale crystal plasticity simulation with measured texture data and hardening laws in macroscale yield criterion simulations in an accurate manner.     

Modelling and experimental characterisation of the Lüders strain in complex loaded ferritic steel compact tension specimens

The use of 3D digital image correlation (DIC) has been used to capture the Lüders strains in a low carbon ferritic steel. Results were used to calibrate and compare with finite element (FE) results based on a constitutive plasticity model, capable of yield drop behaviour and therefore Lüders strains, by Zhang et al. (2001). Tensile tests were carried out at several strain rates to characterise the material behaviour. The results of these tests were used to fit parameters in the constitutive plasticity model. The FE model was then tested on a complex loading situation of in-plane compression of a compact tension (CT) specimen. The FE model predicts the shape and formation of the Lüders bands well. This FE model, using Zhang’s constitutive plasticity model, was used to predict the residual stress profile to compare with standard elastic–plastic isotropic hardening models with no yield point. The yield point reduced both the predicted peak tensile stress, at the notch root, and the amount of plastic strain. In regions where the plastic strain was of a similar size to the Lüders strain the stress profiles were perturbed from flat profiles predicted by the standard elastic–plastic hardening models.