Strain gradient plasticity analysis of transformation induced plasticity in multiphase steels

“To what extent do plastic strain gradients affect the strengthening resulting from the transformation of small metastable inclusions into hard inclusions within a plastically deforming matrix?” is the central question addressed here. Though general in the approach, the focus is on the behavior of TRIP-assisted multiphase steels. A two-dimensional embedded cell model of a simplified microstructure composed of a single metastable austenitic inclusion surrounded by a soft ferritic matrix is considered. The cell is inserted in a large homogenized medium. The transformation of a fraction of the austenite into a hard martensite plate is simulated, accounting for a transformation strain, and leading to complex elastic and plastic accommodation. The size of a transforming plate in real multiphase steels is typically between 0.1 and 2 μm, a range of size in which plastic strain gradient effects are expected to play a major role. The single parameter version of the Fleck–Hutchinson strain gradient plasticity theory is used to describe the plasticity in the austenite, ferrite and martensite phases. The higher order boundary conditions imposed on the plastic flow have a large impact on the predicted strengthening. Using realistic values of the intrinsic length parameter setting the scale at which the gradients effects have an influence leads to a noticeable increase of the strengthening on top of the increase due to the transformation of a volume fraction of the retained austenite. The geometrical parameters such as the volume fraction of retained austenite and of the transforming zone also bring significant strengthening. Strain gradient effects also significantly affect the stress state inside the martensite plate during and after transformation with a potential impact on the damage resistance of these steels.     

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.     

Modelling of laminated microstructures in stress-induced martensitic transformations

This paper is concerned with micromechanical modelling of stress-induced martensitic transformations in crystalline solids, with the focus on distinct elastic anisotropy of the phases and the associated redistribution of internal stresses. Micro-macro transition in stresses and strains is analysed for a laminated microstructure of austenite and martensite phases. Propagation of a phase transformation front is governed by a time-independent thermodynamic criterion. Plasticity-like macroscopic constitutive rate equations are derived in which the transformed volume fraction is incrementally related to the overall strain or stress. As an application, numerical simulations are performed for cubic β₁ (austenite) to orthorhombic γ₁′ (martensite) phase transformation in a single crystal of Cu-Al-Ni shape memory alloy. The pseudoelasticity effect in tension and compression is investigated along with the corresponding evolution of internal stresses and microstructure.     

含亚稳残余奥氏体HSLA TRIP钢的摩擦磨损性能研究

在SST-ST型销-盘摩擦磨损试验机上对特殊工艺热处理、成分为0.22C-1.50Si-1.65Mn的新型相变诱发塑性(HSLA Si-Mn系TRIP)钢的摩擦磨损性能进行研究,采用X射线衍射技术分析了摩擦过程中钢中残余奥氏体的变形诱发转变行为.结果表明,在本试验条件下,与相同成分并在两相区双相化处理的双相钢相比,TRIP钢具有较低摩擦系数和较好耐磨性.由于HSLA Si-Mn系TRIP钢中亚稳残余奥氏体在摩擦应力作用下的变形诱发相变作用,使得其磨损区的大量残余奥氏体向马氏体转变,一方面松弛了磨损表面的内应力,另一方面转变的马氏体使基体强化,从而使其具有较好摩擦磨损特性.     

A 3-D constitutive model for shape memory alloys incorporating pseudoelasticity and detwinning of self-accommodated martensite

A 3-D constitutive model for polycrystalline shape memory alloys (SMAs), based on a modified phase transformation diagram, is presented. The model takes into account both direct conversion of austenite into detwinned martensite as well as the detwinning of self-accommodated martensite. This model is suitable for performing numerical simulations on SMA materials undergoing complex thermomechanical loading paths in stress–temperature space. The model is based on thermodynamic potentials and utilizes three internal variables to predict the phase transformation and detwinning of martensite in polycrystalline SMAs. Complementing the theoretical developments, experimental data are presented showing that the phase transformation temperatures for the self-accommodated martensite to austenite and detwinned martensite to austenite transformations are different. Determination of some of the SMA material parameters from such experimental data is also discussed. The paper concludes with several numerical examples of boundary value problems with complex thermomechanical loading paths which demonstrate the capabilities of the model.     

Effect of phase-transformation on mechanical behavior of dual-phase steel plate

The mechanical behavior of dual phase steel plates is affected by internal stresses created during martensite transformation. Analytical modelling of this effect is made by considering a unit cell made of martensite inclusion in a ferrite matrix. A large strain finite element analysis is then performed to obtain the plane stress deformation state. Displayed numerically are the development of the plastic zone and distribution of local state of stress and strain. Studied also are the shape configuration of the martensite (hard-phase) that influences the interfacial condition as related to stress transmission and damage. Internal stresses are found to enhance the global flow stress after yield initiation in the ferrite matrix. Good agreement is obtained between the analytical results and experimental observations.     

Crystallographically based model for transformation-induced plasticity in multiphase carbon steels

The microstructure of multiphase steels assisted by transformation-induced plasticity consists of grains of retained austenite embedded in a ferrite-based matrix. Upon mechanical loading, retained austenite may transform into martensite, as a result of which plastic deformations are induced in the surrounding phases, i.e., the ferrite-based matrix and the untransformed austenite. In the present work, a crystallographically based model is developed to describe the elastoplastic transformation process in the austenitic region. The model is formulated within a large-deformation framework where the transformation kinematics is connected to the crystallographic theory of martensitic transformations. The effective elastic stiffness accounts for anisotropy arising from crystallographic orientations as well as for dilation effects due to the transformation. The transformation model is coupled to a single-crystal plasticity model for a face-centered cubic lattice to quantify the plastic deformations in the untransformed austenite. The driving forces for transformation and plasticity are derived from thermodynamical principles and include lower-length-scale contributions from surface and defect energies associated to, respectively, habit planes and dislocations. In order to demonstrate the essential features of the model, simulations are carried out for austenitic single crystals subjected to basic loading modes. To describe the elastoplastic response of the ferritic matrix in a multiphase steel, a crystal plasticity model for a body-centered cubic lattice is adopted. This model includes the effect of nonglide stresses in order to reproduce the asymmetry of slips in the twinning and antitwinning directions that characterizes the behavior of this type of lattices. The models for austenite and ferrite are combined to simulate the microstructural behavior of a multiphase steel. The results of the simulations show the relevance of including plastic deformations in the austenite in order to predict a more realistic evolution of the transformation process. This work is part of the research program of the Netherlands Institute for Metals Research (NIMR) and the Stichting voor Fundamenteel Onderzoek der Materie (FOM, financially supported by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO)). The research was carried out under project number 02EMM20 of the FOM/NIMR program “Evolution of the Microstructure of Materials” (P-33).     

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.     

单轴拉伸过程中双相不锈钢局部微观变形行为的细观力学有限元分析

针对双相不锈钢中奥氏体相和铁素体相分别展开了纳米压痕实验,并通过有限元反演得到两相各自的拉伸应力-应变关系,利用Voronoi Tessellation法生成代表性的微结构体积单元,对双相不锈钢的单轴拉伸行为进行了有限元仿真和模拟,研究了双相不锈钢在拉伸过程中的局部应力、应变分布和演化规律．结果表明,利用Voronoi Tessellation法建立单元模型,结合本文通过纳米压痕实验获取的两相力学性能参数,可以很好地模拟双相不锈钢的整体单拉行为,奥氏体比铁素体软,拉伸载荷下双相不锈钢的应变集中在奥氏体中,应力集中在铁素体中;局部应力应变的分布特征与两相分布特征和晶粒形状有关,最大应变值主要集中在奥氏体晶粒狭长且尖锐的区域,而最大应力则主要发生在铁素体晶粒狭长和尖锐的区域;对于奥氏体和铁素体晶粒占比相当的双相不锈钢,其虽然可以具有较为综合的宏观力学性能,但是其微观应力集中的区域和应力最大值相对较大．研究成果为进一步揭示双相不锈钢局部失效机理奠定了基础．     

A CONSTITUTIVE MODEL FOR TRANSFORMATION, REORIENTATION AND PLASTIC DEFORMATION OF SHAPE MEMORY ALLOYS

A constitutive model is developed for the transformation, reorientation and plastic deformation of shape memory alloys (SMAs). It is based on the concept that an SMA is a mixture composed of austenite and martensite, the volume fraction of each phase is transformable with the change of applied thermal-mechanical loading, and the constitutive behavior of the SMA is the combination of the individual behavior of its two phases. The deformation of the martensite is separated into elastic, thermal, reorientation and plastic parts, and that of the austenite is separated into elastic, thermal and plastic parts. Making use of the Tanaka’s transformation rule modified by taking into account the effect of plastic deformation, the constitutive model of the SMA is obtained. The ferroelasticity, pseudoelasticity and shape memory effect of SMA Au-47.5 at.%Cd, and the pseudoelasticity and shape memory effect as well as plastic deformation and its effect of an NiTi SMA, are analyzed and compared with experimental results.     

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.     

Phase proportioning in CuAlBe shape memory alloys during thermomechanical loadings using electric resistance variation

The microstructure of shape memory alloys changes with the thermomechanical history of the material. During thermomechanical loading, austenite, thermally-induced martensite or stress-induced martensite can be simultaneously present in the material. In applications integrating SMA parts, utilization conditions seriously affect the microstructure and can generate macroscopic strain or stress. Consequently, during thermomechanical loadings, it is important to be able to proportion the different phases and consequently to understand the kinetic transformation. This is very useful in the development of constitutive equations. This study shows, by a series of tests, that the proposed experimental method, based on the measurement of the variation of electric resistance of CuAlBe wires, permits to determine the volume fraction of the different phases present in the material (i.e., austenite, stress-induced martensite and thermally-induced martensite). The proposed method is applied to the most common thermomechanical behavior met in engineering applications of shape memory alloys: pseudoelasticity, pseudoplasticity, recovery-stress and stress-assisted two-way shape memory effect. The proportioning method based on a mixture law integrating the resistivity of pure phases present in the SMA is first performed on different two-phase mixture cases and then applied to a three phase mixture case.     

纳米单晶NiTi合金单程形状记忆效应的分子动力学模拟

采用基于第二近邻修正型嵌入原子势的分子动力学方法研究了纳米单晶NiTi合金的单程形状记忆效应,详细阐明了温度诱发马氏体相变和应力诱发马氏体重定向过程中纳米单晶的变形行为和微结构演化,进一步分析了加/卸载速率对NiTi合金单程形状记忆效应的影响。结果表明,NiTi纳米单晶在应力加载过程中发生马氏体重定向,卸载后存在残余应变;当加热到奥氏体转变结束温度以上时,马氏体逆相变为奥氏体相,残余应变逐渐减小,但未完全回复;随着应力加载速率的增加,重定向临界应力和模量逐渐增加;再次降温过程中不同加载速率下的原子结构演化各不相同。     

纳米单晶NiTi合金单程形状记忆效应的分子动力学模拟

采用基于第二近邻修正型嵌入原子势的分子动力学方法研究了纳米单晶NiTi合金的单程形状记忆效应,详细阐明了温度诱发马氏体相变和应力诱发马氏体重定向过程中纳米单晶的变形行为和微结构演化,进一步分析了加/卸载速率对NiTi合金单程形状记忆效应的影响。结果表明,NiTi纳米单晶在应力加载过程中发生马氏体重定向,卸载后存在残余应变;当加热到奥氏体转变结束温度以上时,马氏体逆相变为奥氏体相,残余应变逐渐减小,但未完全回复;随着应力加载速率的增加,重定向临界应力和模量逐渐增加;再次降温过程中不同加载速率下的原子结构演化各不相同。     

Estimation of deformation behavior of TRIP steels — smooth/ringed-notched specimens under monotonic and cyclic loading

Uniaxial tension tests were performed under a constant strain rate and various environmental temperatures from 77 to 373 K to identify the concrete form of the constitutive equation for TRIP steels. To elucidate the dependence of the martensitic transformation on the nonuniform deformation and to validate the proposed constitutive equation, the volume fraction of the martensite phase is predicted and measured using computational simulation and experimental procedures, respectively, for the uniaxial tension of bars with a ringed notch. The good correspondence between the local volume fraction of the martensite phase around the notch, obtained by both methods, verifies the validity of the proposed constitutive equation for the nonuniform deformation behavior. Subsequently, the computational simulations were performed to elucidate the deformation behavior of ringed-notched bars under compression and that of smooth/ringed-notched bars under cyclic loading.     

Experimental characterization and micromechanical modeling of superelastic response of a porous NiTi shape-memory alloy

Porous shape-memory alloys are usually brittle due to the presence of various nickel-titanium intermetallic compounds that are produced in the course of most commonly used synthesizing techniques. We consider here a porous NiTi shape-memory alloy (SMA), synthesized by spark-plasma sintering, that is ductile and displays full shape-memory effects over the entire appropriate range of strains. The porosity however is only 12% but the basic synthesizing technique has potential for producing shape-memory alloys with greater porosity that still are expected to display superelasticity and shape-memory effects. The current material has been characterized experimentally using quasi-static and dynamic tests at various initial temperatures, mostly within the superelastic strain range, but also into the plastic deformation regime of the stress-induced martensite phase. To obtain a relatively constant strain rate in the high strain-rate tests, a novel pulse-shaping technique is introduced. The results of the quasi-static experiments are compared with the predictions by a model that can be used to calculate the stress-strain response of porous NiTi shape-memory alloys during the austenite-to-martensite and reverse phase transformations in uniaxial quasi-static loading and unloading at constant temperatures. In the austenite-to-martensite transformation, the porous shape-memory alloy is modeled as a three-phase composite with the parent phase (austenite) as the matrix and the product phase (martensite) and the voids as the embedded inclusions, reversing the roles of austenite and martensite during the reverse transformation from fully martensite to fully austenite phase. The criterion of the stress-induced martensitic transformation and its reversal is based on equilibrium thermodynamics, balancing the thermodynamic driving force for the phase transformation, associated with the reduction of Gibbs’ free energy, with the resistive force corresponding to the required energy to create new interface surfaces and to overcome the energy barriers posed by various microstructural obstacles. The change in Gibbs’ free energy that produces the driving thermodynamic force for phase transformation is assumed to be due to the reduction of mechanical potential energy corresponding to the applied stress, and the reduction of the chemical energy corresponding to the imposed temperature. The energy required to overcome the resistance imposed by various nano- and subnano-scale defects and like barriers, is modeled empirically, involving three constitutive constants that are then fixed based on the experimental data. Reasonably good correlation is obtained between the experimental and model predictions.     

Numerical investigation of the interaction between the martensitic transformation front and the plastic strain in austenite

Phase-field simulations of the martensitic transformation (MT) in an austenitic matrix which has already undergone the plastic deformation are carried out. For this purpose the elasto-plastic phase-field approach of incoherent MT developed in a previous work [Kundin et al., 2011. A phase-field model for incoherent martensitic transformations including plastic accommodation processes in the austenite. J. Mech. Phys. Solids 59, 2082–2012] is used. The evolution equation for the dislocation density field is extended by taking into account the thermal and athermal annihilation of the dislocations in the austenitic matrix and the athermal annihilation at the transformation front. It is shown that the plastic deformation in the austenite caused by the MT interacts with the dislocation field and the MT front that leads to an inhomogeneous increasing of the total dislocation density. During the phase transformation one part of the dislocations in the austenite is inherited by the martensitic phase and this inheritance depends on the kinetics and the crystallography of MT. Another part of dislocations annihilates at the transformation front and decreases the dislocation density in the growing martensite. Based on the simulation results the specific type of phenomenological dependency between the inherited dislocations, the martensite phase fraction and the plastic deformation is proposed.     

A 2D finite element based on a nonlocal constitutive model describing localization and propagation of phase transformation in shape memory alloy thin structures

Here, the effects of localization and propagation of martensitic phase transformation on the response of SMA thin structures subjected to thermo-mechanical loadings are investigated using nonlocal constitutive model in conjunction with finite element method. The governing equations are derived based on variational principle considering thermo-mechanical equilibrium and the spatial distribution of the nonlocal volume fraction of martensite during transformation. The nonlocal volume fraction of martensite is defined as a weighted average of the local volume fraction of martensite over a domain characterized by an internal length parameter. The local version of the thermo-mechanical behavior model derived from micromechanics considers the local volume fraction of martensite and the mean transformation strain. A 4-noded quadrilateral plane stress element with three degrees of freedom per node accounting for in-plane displacements and the nonlocal volume fraction of martensite is developed. Numerical simulations are conducted to bring out the influence of material and geometrical heterogeneities (perturbations/defects) on the localization and propagation of phase transformation in SMA thin structures. Also, a sensitivity analysis of the material response due to the localization and the other related model parameters is carried out. The detailed investigation done here clearly shows that the localization of phase transformation has significant effect on the response of shape memory alloys.     

具有层状微观结构的NiTi单晶本构模型

建立了应力诱发的具有层状微观结构的NiTi单晶本构模型. 模型考虑母相和马氏体相弹性各向异性性质的差异, 以NiTi单晶相变过程中可能出现的24个马氏体变体为基础, 利用相变驱动力和理想界面的连续条件推导了马氏体相变的发生及发展过程, 以及单晶相变过程中宏微观应力应变的演化, 数值模拟了在不同加载方向材料的应力应变响应.结果表明, 对于不同的加载方向, NiTi单晶既存在强化也存在软化现象.