A crystal plasticity model that incorporates stresses and strains due to slip gradients

This work is concerned with incorporating the kinematic and stress effects of excess dislocations in a constitutive model for the elastoplastic behavior of crystalline materials. The foundation of the model is a three term multiplicative decomposition of the deformation gradient in which the two classical terms of plastic and elastic deformation are included along with an additional term for long range strain due to the collective effects of excess dislocations. The long range strain is obtained from an assumed density of Volterra edge dislocations and is directly related to gradients in slip. A new material parameter emerges which is the size the region about a continuum point that contributes to long range strains.Using Hookean elasticity, the stress at a point is linearly related to the sum of the elastic plus the long range strain fields. However, the driving force for slip is postulated to be due only to the elastic stress so that the long range stress is a back stress in the constitutive relationship for plastic deformation. A consistent balance of the total deformation rate with the three proposed mechanisms of deformation leads to a set of differential equations that can be solved for the elastic stress, rotation and pressure which then implicitly defines the material state and equilibrium stress. Results from the simulation of a tapered tensile specimen demonstrate that the constitutive model exhibits isotropic and kinematic type hardening effects as well as changes in the pattern of plastic deformation and necking when compared to a material without slip gradient effects.     

高氮奥氏体不锈钢高温动态响应特性及本构关系

在293~873 K的环境下,采用分离式霍普金森杆装置对高氮钢试样进行了102~103 s-1应变率下的动态加载实验。结合准静态实验结果,分析了应变率和温度对材料塑性流动特性的影响。结果表明:高氮钢的动态力学行为具有很强的应变率敏感性和温度敏感性。当应变率达到400 s-1或更高时,流动应力随应变率的增加显著升高;在同一应变率下,流动应力随温度的降低明显升高。研究了温度和应变率耦合效应对材料塑性行为的影响,得出温度软化效应在高氮钢高温动态塑性变形中起主导作用。基于经典的Johnson-Cook（J-C）模型,通过对实验数据的分析,得出了高氮钢材料的修正J-C本构方程,经验证修正J-C方程预测结果与实验结果吻合。     

A dislocation-mechanics-based constitutive model for dynamic strain aging

Dynamic strain aging (DSA) is an important phenomenon in solute hardened metals and seriously affects the mechanical properties of metals. DSA is generally induced by the interaction between the moving dislocations and the mobile solute atoms. In this paper, only the interaction between moving dislocations and mobile solute atoms in a dislocation core area (core atmosphere) will be taken into account. To establish the constitutive model which can describe the DSA phenomenon, we improved the Zerilli-Armstrong dislocation-mechanics-based thermal viscoplastic constitutive relation, and added the effect of the interaction between the moving dislocations and core atmosphere. Because the constitutive relation established is based on the Zerilli-Armstrong relation, it can describe not only the DSA phenomenon, but also the mechanical behavior of metals over a broad range of temperatures (77K∼1000K) and strain rate (10⁻⁴∼10⁴ s⁻¹). The model prediction for tantalum fits well with the experimental data. Projected supported by the Chinese Academy of Sciences and the High Technical Project.     

第三型应变时效的提出与研究进展

第三型应变时效现象的发现使得传统的对于金属塑性流动行为的认识、位错的热激活理论以及常见的金属热粘塑性本构模型均需要进一步完善。为了系统地认识第三型应变时效,首先介绍了第三型应变时效现象区别于静态应变时效和Portevin-Le Chatelie动态应变时效的宏观特征,其次,对第三型应变时效的微观机理以及第三型应变时效与Portevin-Le Chatelier动态应变时效、蓝脆现象以及机械波谱的关联性进行了系统总结。最后,介绍了包含第三型应变时效的金属热黏塑性本构模型的发展。     

超高强度钢AF1410塑性流动特性及其本构关系

在本文中,为揭示超高强度钢AF1410的塑性流动性,并研究其塑性流动本构关系,利用CSS4410电子万能试验机和改进的Hopkinson拉压杆技术,对AF1410钢在温度从100K到600K,应变率从0.001/s到2000/s,塑性应变超过20%的塑性流动特性进行了试验研究。结果表明,拉伸加载下AF1410钢屈服强度低于压缩屈服强度,且随应变率增加,拉压屈服强度差值越来越大;该材料塑性流动应力对应变率敏感性低,而对温度较为敏感;随应变率的提高,该材料拉伸失效应变减小,但温度对失效应变无明显影响。最后基于位错的运动学关系,借助试验数据,获得了AF1410钢的塑性流动物理概念本构模型,并通过与经典J-C模型的结果对比对该物理概念本构模型进行了评估分析,表明该物理概念本构模型在较宽温度和应变率范围能较好的预测AF1410钢的塑性流动应力。     

金属热黏塑性本构关系的研究进展

金属材料的塑性流动行为依赖于温度和应变率,温度和应变率敏感性是金属材料塑性流动的最重要的本质特性之一,建立合适的热黏塑性本构关系来准确描述金属塑性流动行为的温度和应变率依赖性,是金属材料能被广泛应用的必要前提。为此,对金属热黏塑性本构关系的最新研究进展进行了综述,介绍了常见的几种金属热黏塑性本构关系并进行了详细讨论,给出了各本构关系的优势与不足,最后系统介绍了包含金属塑性流动行为中出现的第三型应变时效、或K-W锁位错结构引起的流动应力随温度变化出现的反常应力峰以及拉压不对称等行为的金属热黏塑性本构关系的研究进展。     

考虑微观结构特征长度演化的内变量黏塑性本构模型

在金属晶体材料高应变率大应变变形过程中,存在强烈的位错胞尺寸等微观结构特征长度细化现象,势必对材料加工硬化、宏观塑性流动应力产生重要影响。基于宏观塑性流动应力与位错胞尺寸成反比关系,提出了一种新型的BCJ本构模型。利用位错胞尺寸参数,修正了BCJ模型的流动法则、内变量演化方程,引入了考虑应变率和温度相关性的位错胞尺寸演化方程,建立了综合考虑微观结构特征长度演化、位错累积与湮灭的内变量黏塑性本构模型。应用本文模型,对OFHC铜应变率在10-4~103 s-1、温度在298~542 K、应变在0~1的实验应力-应变数据进行了预测。结果表明:在较宽应变率、温度和应变范围内,本文模型的预测数据与实验吻合很好;与BCJ模型相比,对不同加载条件下实验数据的预测精度均有较大程度的提高,最大平均相对误差从9.939%减小为5.525%。     

Size effects and idealized dislocation microstructure at small scales: Predictions of a Phenomenological model of Mesoscopic Field Dislocation Mechanics: Part I

A Phenomenological Mesoscopic Field Dislocation Mechanics (PMFDM) model is developed, extending continuum plasticity theory for studying initial-boundary value problems of small-scale plasticity. PMFDM results from an elementary space-time averaging of the equations of Field Dislocation Mechanics (FDM), followed by a closure assumption from any strain-gradient plasticity model that attempts to account for effects of geometrically necessary dislocations (GNDs) only in work hardening. The specific lower-order gradient plasticity model chosen to substantiate this work requires one additional material parameter compared to its conventional continuum plasticity counterpart. The further addition of dislocation mechanics requires no additional material parameters. The model (a) retains the constitutive dependence of the free-energy only on elastic strain as in conventional continuum plasticity with no explicit dependence on dislocation density, (b) does not require higher-order stresses, and (c) does not require a constitutive specification of a ‘back-stress’ in the expression for average dislocation velocity/plastic strain rate. However, long-range stress effects of average dislocation distributions are predicted by the model in a mechanistically rigorous sense. Plausible boundary conditions (with obvious implication for corresponding interface conditions) are discussed in some detail from a physical point of view. Energetic and dissipative aspects of the model are also discussed. The developed framework is a continuous-time model of averaged dislocation plasticity, without having to rely on the notion of incremental work functions, their convexity properties, or their minimization. The tangent modulus relating stress rate and total strain rate in the model is the positive-definite tensor of linear elasticity, and this is not an impediment to the development of idealized microstructure in the theory and computations, even when such a convexity property is preserved in a computational scheme. A model of finite deformation, mesoscopic single crystal plasticity is also presented, motivated by the above considerations.Lower-order gradient plasticity appears as a constitutive limit of PMFDM, and the development suggests a plausible boundary condition on the plastic strain rate for this limit that is appropriate for the modeling of constrained plastic flow in three-dimensional situations.     

Dynamic strain aging and related instabilities: experimental,theoretical and numerical aspects

Experimental data from uniaxial tensile tests on smooth and notched specimens of aluminium alloy 5083-H116 show that the material exhibits negative strain-rate sensitivity for strain rates within a certain range. The negative strain-rate dependence, which is attributed to dynamic strain aging, leads to serrated stress–strain curves, discontinuous plastic flow and propagating deformation bands during plastic straining (also denoted as the Portevin–Le Chatelier effect). Band analysis and linear perturbation analysis are performed using simple elastic-viscoplastic constitutive equations that include negative strain-rate sensitivity in a simplified manner. The negative strain-rate sensitivity allows for jumps in the plastic strain rate, which in turn permits the existence of localisation bands for the elastic-viscoplastic model. The simple elastic-viscoplastic constitutive model has been implemented in LS-DYNA, and non-linear finite element simulations of smooth and notched tensile test specimens are performed, allowing more detailed investigations into the effects of the negative strain-rate sensitivity on the material's behaviour.     

Three-dimensional constitutive relations of ideal plasticity and the flow on the coulomb-tresca prism edge

In the present paper, we consider basic relations of the mathematical theory of plasticity for the spatial state corresponding to the edge of the Coulomb-Tresca prism, which follow from the generalized associated flow law restricting the plastic flow freedom for the above states to the minimal possible extent. We found that the spatial relations of the theory of plasticity, formulated by A. Yu. Ishlinsky in 1946, can be derived from the above version of the theory of flow. We show that the A. Yu. Ishlinsky constitutive relations for states on the Coulomb-Tresca prism edge express the commutativity of the stress tensor and the tensor of plastic strain increments. We obtained one explicit form of the constitutive relation relating the stress tensor to the plastic strain increments for the stressed states corresponding to the Coulomb-Tresca prism edge.     

Constitutive modeling of pressure dependent plasticity and fracture in solder joints

A new constitutive model describing the pressure dependence of plasticity and fracture in solder joints is proposed. The pressure dependence of the solders flow stress is obtained from the Peierls stress necessary to move mobile dislocations. Conservation of volume during plastic deformation is achieved through a deviatoric associated flow rule. Fracture is considered as a result of damage accumulation depending on the amount of plastic strain and triaxiality of stress the sample has experienced. The constitutive model was implemented in the FEM code Abaqus using the user subroutine UMAT. The model predictions were compared to experimental results of tensile tests performed with Sn-3.5Ag solder joints of various thicknesses. Ultimate tensile strength and ductility showed a strong dependency on the solder joints geometry. It is demonstrated that the model can nicely be fitted to experimental results.     

FeCrNi合金静动态物理本构模型研究

以金属材料塑性变形的位错动力学为基础,将FeCrNi合金的流动应力分解为非热应力和热激活应力两部分.通过对该合金屈服应力随温度变化特性、屈服应力的应变速率特性、孪晶组织的温度特性及位错组态的应变速率特性进行分析,认为非热应力不只是应变的函数,还与温度和应变速率相关,因此对Johnson-Cook模型方程形式进行修正以描述非热应力. 同时认为影响热激活应力的微结构参数主要为位错阻碍间距\Deltal, 定义并推导出表征\Delta l演化的g函数的表达式,将其引入Kocks的热激活方程,从而建立FeCrNi合金的物理型本构模型.该模型初步实现了对FeCrNi合金从室温到高温、从准静态到动态塑性变形行为的描述.      

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.     

Al-Cu多晶合金中锯齿形屈服现象的时序统计研究

合金中溶质原子和可动位错之间的交互作用引起材料的固溶硬化,这一交互作用在合适的应变率和温度条件下可重复发生,称为动态应变时效(DSA),从而导致锯齿形的应力屈服,同时伴随着变形带的产生,被称为Portevin-Le Chartelier(PLC)效应。本文对Al-Cu多晶合金锯齿形屈服现象中应力跌幅、跌落时间等多个特征物理量时序演化规律进行了系统的研究。结果显示,跌落时间对应变不敏感;1mm和2mm厚度试件的应力跌幅和再加载时间都近似随应变线性增加,3mm厚试件的统计结果则显示出一个明显的转变点。最后本文结合动态应变时效原理和晶体位错学,分析了实验条件(应变率和试件尺寸)对各特征物理量的影响。     

Roles of plastic spin in shear banding

In this paper, the effects of plastic spin on shear banding and simple shear are examined systematically. Three types of plastic constitutive model with plastic spin are considered: (i) a non-coaxial model in which the direction of the plastic strain rate depends on that of the stress rate; (ii) a strain-softening model based on the J₂ flow theory; and (iii) the pressure-sensitive porous plasticity model. All the constitutive models are formulated in viscoplastic forms and in conjunction with non-local concepts that have been recently focused and discussed. First, behavior in simple shear is examined by numerical analysee with the aforementioned constitutive models. Moreover, some experimental evidences for stress response to simple shear are shown; that is, several large torsion tests of metal tubes and bars are carried out. Next, finite element simulations of shear banding in plane strain tension are performed. A critical effect of plastic spin on shear banding is observed for the noncoaxial model, while an almost negligible effect is observed for the porous model. The identical effects of plastic spin are observed, whether nonlocality exists or not. Finally, we discuss the relationship between the behavior in simple shear and the shear band formation. It is emphasized that this is a critical issue in predicting shear banding in macroscopic grounds.     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

6061-T6铝合金动态拉伸本构关系及失效行为

采用HMH-206高速材料试验机开展了6061-T6铝合金在0.001～100 s-1应变率范围内的静、动态拉伸力学性能实验,分析了其应力-应变响应特征和应变率敏感性,讨论了应变率对6061-T6铝合金流动应力和应变率敏感性指数的影响,并基于实验结果对Johnson-Cook本构模型进行了修正。结合缺口试件的实验结果和模拟数据,得到了材料的Johnson-Cook失效模型参数,并对模型的准确性和适用性进行了验证。结果表明,在拉伸载荷作用下,6061-T6铝合金表现出明显的应变硬化特征和应变率敏感性,其流动应力随应变率的升高而提高,修正的Johnson-Cook本构模型可以描述材料的动态塑性流动行为,建立的Johnson-Cook失效模型能够表征材料的断裂失效行为。     

Plastic deformation modeling of AL-6XN stainless steel at low and high strain rates and temperatures using a combination of bcc and fcc mechanisms of metals

Combination of physically based constitutive models for body centered cubic (bcc) and face centered cubic (fcc) metals developed recently by the authors [Voyiadjis, G.Z., Abed, F.H., 2005. Microstructural based models for bcc and fcc metals with temperature and strain rate dependency. Mech. Mater. 37, 355–378] are used in modeling the plastic deformation of AL-6XN stainless steel over a wide range of strain rates between 0.001 and 8300 s⁻¹ at temperatures from 77 to 1000 K. The concept of thermal activation analysis as well as the dislocation interaction mechanism is used in developing the plastic flow model for both the isothermal and adiabatic plastic deformation. In addition, the experimental observations of AL-6XN conducted by Nemat-Nasser et al. [Nemat-Nasser, S., Guo, W., Kihl, D., 2001. Thermomechanical response of AL-6XN stainless steel over a wide range of strain rates and temperatures, J. Mech. Phys. Solids 49, 1823–1846] are utilized in understanding the underlying deformation mechanisms. The plastic flow is considered in the range of temperatures and strain rates where diffusion and creep are not dominant, i.e., the plastic deformation is attributed to the motion of dislocations only. The modeling of the true stress–true strain curves for AL-6XN stainless steel is achieved using the classical secant modulus for the case of unidirectional deformation. The model parameters are obtained using the experimental results of three strain rates (0.001, 0.1, and 3500 s⁻¹). Good agreement is obtained between the experimental results and the model predictions. Moreover, the independency of the present model to the experiments used in the modeling is verified by comparing the theoretical results to an independent set of experimental data at the strain rate of 8300 s⁻¹ and various initial temperatures. Good correlation is observed between the model predictions and the experimental observations.