A two-dimensional dislocation dynamics model of the plastic deformation of polycrystalline metals

Two-dimensional dislocation dynamics (2D-DD) simulations under fully periodic boundary conditions are employed to study the relation between microstructure and strength of a material. The material is modeled as an elastic continuum that contains a defect microstructure consisting of a preexisting dislocation population, dislocation sources, and grain boundaries. The mechanical response of such a material is tested by uniaxially loading it up to a certain stress and allowing it to relax until the strain rate falls below a threshold. The total plastic strain obtained for a certain stress level yields the quasi-static stress-strain curve of the material. Besides assuming Frank-Read-like dislocation sources, we also investigate the influence of a pre-existing dislocation density on the flow stress of the model material. Our results show that - despite its inherent simplifications - the 2D-DD model yields material behavior that is consistent with the classical theories of Taylor and Hall-Petch. Consequently, if set up in a proper way, these models are suited to study plastic deformation of polycrystalline materials.     

A modified Green–Lindsay thermoelasticity with strain rate to eliminate the discontinuity

Probing the mechanism of ultrafast thermoelastic processes is becoming increasingly important with the development of laser-assisted micro/nano machining. Although thermoelastic models containing temperature rate have been historically proposed, the strain rate has not been considered yet. In this work, a generalized thermoelastic model is theoretically established by introducing the strain rate in Green–Lindsay (GL) thermoelastic model with the aid of extended thermodynamics. Numerically, a semi-infinite one-dimensional problem is considered with traction free at one end and subjected to a temperature rise. The problem is solved using the Laplace transform method, and the transient responses, i.e. displacement, temperature and stresses are graphically depicted. Interestingly, it is found that the strain rate may eliminate the discontinuity of the displacement at the elastic and thermal wave front. Also, the present model is compared with Green–Naghdi (GN) models. It is found that the thermal wave speed of the present model is faster than GN model without energy dissipation, and slower than GN model with energy dissipation. In addition, the thermoelastic responses from the present model are the largest. The present model based upon GL model is free of the jump of GL model in the displacement distribution, and is safer in engineering practices than GN model. The present work will benefit the theoretical modeling and numerical prediction of thermoelastic process, especially for those under extreme fast heating.     

Dual-phase steel sheets under cyclic tension–compression to large strains: Experiments and crystal plasticity modeling

In this work, we develop a physically-based crystal plasticity model for the prediction of cyclic tension–compression deformation of multi-phase materials, specifically dual-phase (DP) steels. The model is elasto–plastic in nature and integrates a hardening law based on statistically stored dislocation density, localized hardening due to geometrically necessary dislocations (GNDs), slip-system-level kinematic backstresses, and annihilation of dislocations. The model further features a two level homogenization scheme where the first level is the overall response of a two-phase polycrystalline aggregate and the second level is the homogenized response of the martensite polycrystalline regions. The model is applied to simulate a cyclic tension–compression–tension deformation behavior of DP590 steel sheets. From experiments, we observe that the material exhibits a typical decreasing hardening rate during forward loading, followed by a linear and then a non-linear unloading upon the load reversal, the Bauschinger effect, and changes in hardening rate during strain reversals. To predict these effects, we identify the model parameters using a portion of the measured data and validate and verify them using the remaining data. The developed model is capable of predicting all the particular features of the cyclic deformation of DP590 steel, with great accuracy. From the predictions, we infer and discuss the effects of GNDs, the backstresses, dislocation annihilation, and the two-level homogenization scheme on capturing the cyclic deformation behavior of the material.     

A physically based model of temperature and strain rate dependent yield in BCC metals: Implementation into crystal plasticity

In this work, we develop a crystal plasticity finite element model (CP-FEM) that constitutively captures the temperature and strain rate dependent flow stresses in pure BCC refractory metals. This model is based on the kink-pair theory developed by Seeger (1981) and is calibrated to available data from single crystal experiments to produce accurate and convenient constitutive laws that are implemented into a BCC crystal plasticity model. The model is then used to predict temperature and strain rate dependent yield stresses of single and polycrystal BCC refractory metals (molybdenum, tantalum, tungsten and niobium) and compared with existing experimental data. To connect to larger length scales, classical continuum-scale constitutive models are fit to the CP-FEM predictions of polycrystal yield stresses. The results produced by this model, based on kink-pair theory and with origins in dislocation mechanics, show excellent agreement with the Mechanical Threshold Stress (MTS) model for temperature and strain-rate dependent flow. This framework provides a method to bridge multiple length scales in modeling the deformation of BCC metals.     

Modeling anisotropic strain hardening and deformation textures in low stacking fault energy fcc metals

The main issues and challenges involved in modeling anisotropic strain hardening and deformation textures in the low stacking fault energy (SFE) fcc metals (e.g. brass) are reviewed and summarized in this paper. The objective of these modeling efforts is to capture quantitatively the major differences between the low SFE fcc metals and the medium (and high) SFE fcc metals (e.g. copper) in the stress–strain response and the deformation textures. While none of the existing models have demonstrated success in capturing the anisotropy in the stress–strain response of the low SFE fcc metals, their apparent success in predicting the right trend in the evolution of deformation texture is also questionable. There is ample experimental evidence indicating that the physical mechanism of the transition from the copper texture to the brass texture is represented wrongly in these models. These experimental observations demonstrate clearly the need for a new approach in modeling the deformation behavior of low SFE fcc metals. This paper reports new approaches for developing crystal plasticity models for the low SFE fcc metals that are consistent with the reported experimental observations in this class of metals. The successes and failures of these models in capturing both the anisotropic strain hardening and the deformation textures in brass are discussed in detail.     

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.     

A mesoscale investigation of strain rate effect on dynamic deformation of single-crystal copper

A combined finite element (FE) simulation and discrete dislocation dynamics (DD) approach has been developed in this paper to investigate the dynamic deformation of single-crystal copper at mesoscale. The DD code yields the plastic strain based on the slip of dislocations and serves as a substitute for the 3D constitutive form used in the usual FE computation, which is implemented into ABAQUS/Standard with a user-defined material subroutine. On the other hand, the FE code computes the displacement and stress field during the dynamic deformation. The loading rate effects on the yield stress and the deformation patterning of single-crystal copper are investigated. With the increasing of strain rate, the yield stress of single-crystal copper increases rapidly. A critical strain rate exists in each single-crystal copper block for the given size and dislocation sources, below which the yield stress is relatively insensitive to the strain rate. The dislocation patterning changes from non-uniform to uniform under high-strain-rate. The shear stresses in the bands are higher than that in the neighboring regions, which are formed shear bands in the crystal. The band width increases with the strain rate, which often take places where the damage occurs.     

Modeling of thermally activated dislocation glide and plastic flow through local obstacles

A unified phenomenological model is developed to study the dislocation glide through weak obstacles during the first stage of plastic deformation in metals. This model takes into account both the dynamical responses of dislocations during the flight process and thermal activations while dislocations are bound by obstacle arrays. The average thermal activation rate is estimated using an analytical model based on the generalized Friedel relations. Then, the average flight velocity after an activation event is obtained numerically by discrete dislocation dynamics (DD). To simulate the dynamical dislocation behavior, the inertia term is implemented into the equation of dislocation motion within the DD code. The results from the DD simulations, coupled with the analytical model, determine the total dislocation velocity as a function of the stress and temperatures. By choosing parameters typical of the face centered cubic metals, the model reproduces both obstacle control and drag control motion in low and high velocity regimes, respectively. As expected by other string models, dislocation overshoots of obstacles caused by the dislocation inertia at the collisions are enhanced as temperature goes down.     

高速冲击表面处理对金属材料力学性能和组织结构的影响

高速冲击表面处理过程中的应变率对金属材料的宏观力学性能和微观组织结构都具有重要影响。根据当前应变率效应的研究成果,从宏观与微观相结合的角度出发,综述了高速冲击表面处理过程中应变率对金属材料强度和塑性的影响规律,并重点阐述了不同应变率下金属材料内部微观组织结构的演变规律,主要包括晶粒结构、绝热剪切带、相变、位错组态和析出相以及变形孪晶等。此外,还分析了组织结构随应变率的演化和微观变形机制的转变对材料力学性能的强化和弱化机理。最后,对高速冲击表面处理梯度组织的变形特点进行了总结。提出了不同组织结构对材料性能影响的综合效应模型,以期为应变率效应的深入研究奠定基础。     

Mechanical effects of the self-organization of dislocations and related critical characteristics

The effects of the dislocation pattern formed due to the self-organization of the dislocations in crystals on the macroscopic hardening and dynamic internal friction (DIF) during deformation are studied. The classic dislocation models for the hardening and DIF corresponding to the homogeneous dislocation configuration are extended to the case for the non-homogeneous one. In addition, using the result of dislocation patterning deduced from the non-linear dislocation dynamics model for single slip, the correlation between the dislocation pattern and hardening as well as DIF is obtained. It is shown that in the case of the tension with a constant strain rate, the bifurcation point of dislocation patterning corresponds to the turning point in the stress versus strain and DIF versus strain curves. This result along with the critical characteristics of the macroscopic behavior near the bifurcation point is microscopically and macroscopically in agreement with the experimental findings on mono-crystalline pure aluminum at temperatures around 0.5T _m . The present study suggests that measuring the DIF would be a sensitive and useful mechanical means in order to study the critical phenomenon of materials during deformation. The project supported by the National Natural Science Foundation of China under the Grand 19702019 & 19891180-4 and by the Chinese Academy of Sciences under the Grand KJ951-1-201.     

A variational multiscale constitutive model for nanocrystalline materials

This paper presents a variational multi-scale constitutive model in the finite deformation regime capable of capturing the mechanical behavior of nanocrystalline (nc) fcc metals. The nc-material is modeled as a two-phase material consisting of a grain interior phase and a grain boundary effected zone (GBAZ). A rate-independent isotropic porous plasticity model is employed to describe the GBAZ, whereas a crystal-plasticity model which accounts for the transition from partial dislocation to full dislocation mediated plasticity is employed for the grain interior. The constitutive models of both phases are formulated in a small strain framework and extended to finite deformation by use of logarithmic and exponential mappings. Assuming the rule of mixtures, the overall behavior of a given grain is obtained via volume averaging. The scale transition from a single grain to a polycrystal is achieved by Taylor-type homogenization where a log-normal grain size distribution is assumed. It is shown that the proposed model is able to capture the inverse Hall-Petch effect, i.e., loss of strength with grain size refinement. Finally, the predictive capability of the model is validated against experimental results on nanocrystalline copper and nickel.     

Multiresolution modeling of ductile reinforced brittle composites

Tungsten carbide-cobalt (WC-Co) is an important ductile reinforced brittle composite used in a range of important applications. The relationship between microstructure and mechanical properties of WC-Co is truly multiscale; micromechanical processes interact at different scales, resulting in permanent plastic deformation, damage accumulation and final failure of the composite. The goal of the current paper is to develop a continuum-based model, which captures the progressively finer scales of strain localization observed in WC-Co composites during plastic deformation and failure. This is achieved via a set of multiresolution governing equations; a microstress is introduced at each scale of strain localization, which represents the resistance to inhomogeneous strain localization at that scale. The extra constitutive models associated with these microstresses can be elucidated from the average response of separate computational cell models of a representative microstructure. The final multiresolution continuum model is capable of capturing the important length scales of deformation during the plastic stage of deformation without resorting to modeling microstructural scale features directly. The result is a more realistic continuum model; in particular the fracture toughness prediction is more physical when these length scales are incorporated compared to a conventional continuum approach.     

铝合金板材塑性性能的研究

铝合金在汽车工业中的广泛应用对于降低汽车重量、减少燃油消耗和汽车尾气的排放量具有十分重要的意义，但其室温塑性成形性能却受到了锯齿形屈服行为的影响，从而制约了铝合金进一步的推广应用。本文基于合金材料塑性变形过程中位错和溶质原子间相互作用的分析，建立了一个可用于描述锯齿形屈服现象的唯象本构模型。该模型将溶质原子对位错运动的钉扎效应和位错挣脱后的脱钉效应置于一个统一的框架内进行考虑，而这两个效应的相互竞争将决定材料宏观变形行为的发展演化。基于该模型的数值模拟结果和实验测试结果取得了良好的一致性，从而验证了理论和模型的有效性。     

基于晶体塑性理论研究铝材料高压高应变率下的强度特性

为了了解金属材料在极端加载下复杂动态响应过程中的多种机制和效应，重点针对Al材料在高压、高应变率加载下的塑性变形机制，在经典晶体塑性模型的基础上，对其中的非线性弹性、位错动力学和硬化形式进行改进，建立适用于高压、高应变率加载下的热弹-黏塑性晶体塑性模型。该模型可以较好地描述单晶铝和多晶铝材料屈服强度随压力的变化过程，相比宏观模型，用该模型还获得了多晶Al材料在冲击加载下的织构演化规律，揭示了织构择优取向行为和压力的关系。     

Constitutive modelling of the high strain rate behaviour of interstitial-free steel

A physically based modelling and experimental investigation of the work hardening behaviour of IF steel covering a wide range of strain rates including complex strain path and/or strain rate changes are presented. In order to obtain isothermal stress–strain curves at high strain rates, a procedure has been proposed with the aid of finite element analysis. The result reveals that the apparent excess of the flow stress after a jump in strain rate, which is frequently observed in bcc metals, is in fact due to the thermal softening at large strains, and that the flow stress after a jump in strain rate tends asymptotically to the values corresponding to the curve at the new strain rate. The strain rate affects not only the short-range stress but also the long-range stress via the strain-rate dependant evolution of dislocation structures. The proposed model is based on the dislocation model of intragranular hardening proposed by Teodosiu and Hu [Teodosiu, C., Hu, Z., 1995. Evolution of the intragranular microstructure at moderate and large strains: modelling and computational significance. In: Shen, S., Dawson, P. R., (Eds.), Proceedings of Numiform'95 on Simulation of Materials Processing: Theory, Methods and Applications. Balkema, Rotterdam, pp. 173–182] and extended to strain rate sensitive one with applying the results of the thermal activation analysis. A satisfactory agreement has been achieved between model predictions and experimental results.     

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

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

刃型位错理论模型的实验验证

刃型位错芯周围变形场的实验测量是多年来非常困难的研究任务,它导致目前有多种位错理论模型并存。为了检验刃型位错理论模型的适用性,使用透射电子显微镜直接观察并获得了多晶金中刃型位错的高分辨电子显微图像,并采用几何相位分析方法测量了刃型位错周围的位移场和应变场。将实验测量结果与线弹性理论位错模型、Peierls-Nabarro位错模型及Fore-man(a=4)位错模型进行了比较。结果表明,三种位错理论模型在远离位错芯的区域都能描述刃型位错变形场,但在距离位错芯较近的区域,Peierls-Nabarro模型是最适当的位错理论模型。     

A constitutive model accounting for strain ageing effects on work-hardening. Application to a C–Mn steel

One of the most successful models for describing the Portevin–Le Chatelier effect in engineering applications is the Kubin–Estrin–McCormick model (KEMC). In the present work, the influence of dynamic strain ageing on dynamic recovery due to dislocation annihilation is introduced in order to improve the KEMC model. This modification accounts for additional strain hardening rate due to limited dislocation annihilation by the diffusion of solute atoms and dislocation pinning at low strain rate and/or high temperature. The parameters associated with this novel formulation are identified based on tensile tests for a C–Mn steel at seven temperatures ranging from 20 °C to 350 °C. The validity of the model and the improvement compared to existing models are tested using 2D and 3D finite element simulations of the Portevin–Le Chatelier effect in tension.