A phase-field theory of dislocation dynamics, strain hardening and hysteresis in ductile single crystals

A phase-field theory of dislocation dynamics, strain hardening and hysteresis in ductile single crystals is developed. The theory accounts for: an arbitrary number and arrangement of dislocation lines over a slip plane; the long-range elastic interactions between dislocation lines; the core structure of the dislocations resulting from a piecewise quadratic Peierls potential; the interaction between the dislocations and an applied resolved shear stress field; and the irreversible interactions with short-range obstacles and lattice friction, resulting in hardening, path dependency and hysteresis. A chief advantage of the present theory is that it is analytically tractable, in the sense that the complexity of the calculations may be reduced, with the aid of closed form analytical solutions, to the determination of the value of the phase field at point-obstacle sites. In particular, no numerical grid is required in calculations. The phase-field representation enables complex geometrical and topological transitions in the dislocation ensemble, including dislocation loop nucleation, bow-out, pinching, and the formation of Orowan loops. The theory also permits the consideration of obstacles of varying strengths and dislocation line-energy anisotropy. The theory predicts a range of behaviors which are in qualitative agreement with observation, including: hardening and dislocation multiplication in single slip under monotonic loading; the Bauschinger effect under reverse loading; the fading memory effect, whereby reverse yielding gradually eliminates the influence of previous loading; the evolution of the dislocation density under cycling loading, leading to characteristic ‘butterfly’ curves; and others.     

Overall softening and anisotropy related with the formation and evolution of dislocation cell structures

In this work, a model, based on a representation of the dislocation cell microstructures by a non-local two-phase material with evolving microstructures, is proposed for the elastic–plastic behavior of metals under monotonic and sequential loading. The first phase represents the cell interior and the second one, the cell walls. The evolution of the microstructure is taken into account considering the cell-wall interfaces as free boundaries. Finally, the accumulation within walls of dislocations crossing the cells defines a non-local hardening process. Assuming a piecewise uniform plastic strain field and assuming ellipsoidal cells, the free energy of the system is calculated. The driving and critical forces associated with the plastic flow of the two-phases and the morphology of the cells are established. In a third part, numerical results are presented for monotonic and sequential loading. The results show an overall softening related to the destabilization of the dislocation microstructures which occurs in sequential as well as monotonic paths.     

A crystal plasticity theory for latent hardening by glide twinning through dislocation transmutation and twin accommodation effects

Imagine a residual glide twin interface advancing in a grain under the action of a monotonic stress. Close to the grain boundary, the shape change caused by the twin is partly accommodated by kinks and partly by slip emissions in the parent; the process is known as accommodation effects. When reached by the twin interface, slip dislocations in the parent undergo twinning shear. The twinning shear extracts from the parent dislocation a twinning disconnection, and thereby releases a transmuted dislocation in the twin. Transmutation populates the twin with dislocations of diverse modes. If the twin deforms by double twinning, double-transmutation occurs even if the twin retwins by the same mode or detwins by a stress reversal. If the twin deforms only by slip, transmutation is single. Whether single or double, dislocation transmutation is irreversible. The multiplicity of dislocation modes increases upon strain, since the twin finds more dislocations to transmute upon further slip of the parent and further growth of the twin. Thus, the process induces an increasing latent hardening rate in the twin. Under profuse twinning conditions, typical of double-lattice structures, this rate-increasing latent hardening combined with crystal rotation to hard orientations by twinning is consistent with a regime of increasing hardening rate, known as Regime II or Regime B. In this paper, we formulate governing equation of the above transmutation and accommodation effects in a crystal plasticity framework. We use the dislocation density based model originally proposed by Beyerlein and Tomé (2008) to derive the effect of latent hardening in a transmuting twin. The theory is expected to contribute to surmounting the difficulty that current models have to simultaneously predict under profuse twinning, the stress-strain curves, intermediate deformation textures, and intermediate twin volume fractions.     

Discrete dislocation dynamics modelling of mechanical deformation of nickel-based single crystal superalloys

Discrete dislocation dynamics (DDD) has been used to model the deformation of nickel-based single crystal superalloys with a high volume fraction of precipitates at high temperature. A representative volume cell (RVC), comprising of both the precipitate and the matrix phase, was employed in the simulation where a periodic boundary condition was applied. The dislocation Frank-Read sources were randomly assigned with an initial density on the 12 octahedral slip systems in the matrix channel. Precipitate shearing by superdislocations was modelled using a back force model, and the coherency stress was considered by applying an initial internal stress field. Strain-controlled loading was applied to the RVC in the [0 0 1] direction. In addition to dislocation structure and density evolution, global stress-strain responses were also modelled considering the influence of precipitate shearing, precipitate morphology, internal microstructure scale (channel width and precipitate size) and coherency stress. A three-stage stress-strain response observed in the experiments was modelled when precipitate shearing by superdislocations was considered. The polarised dislocation structure deposited on the precipitate/matrix interface was found to be the dominant strain hardening mechanism. Internal microstructure size, precipitate shape and arrangement can significantly affect the deformation of the single crystal superalloy by changing the constraint effect and dislocation mobility. The coherency stress field has a negligible influence on the stress-strain response, at least for cuboidal precipitates considered in the simulation. Preliminary work was also carried out to simulate the cyclic deformation in a single crystal Ni-based superalloy using the DDD model, although no cyclic hardening or softening was captured due to the lack of precipitate shearing and dislocation cross slip for the applied strain considered.     

Inspection of free energy functions in gradient crystal plasticity

The dislocation density tensor computed as the cud of plastic distortion is regarded as a new constitutive variable in crystal plasticity. The dependence of the free energy function on the dislocation density tensor is explored starting from a quadratic ansatz. Rank one and logarithmic dependencies are then envisaged based on considerations from the statistical theory of dislocations. The rele- vance of the presented free energy potentials is evaluated from the corresponding analytical solutions of the periodic two-phase laminate problem under shear where one layer is a single crystal material undergoing single slip and the second one remains purely elastic.     

基于微结构动态演化机制的单晶镍基高温合金晶体塑性本构及其有限元模拟

因其优异的高温力学性能,镍基单晶高温合金在航空航天和能源等领域得到了广泛的应用.镍基单晶高温合金优异的高温性能来源于其特有的两相微结构.基于代表体胞模型及分块均匀化方法,以位错密度为主要内变量,发展了一个包含两相微结构和位错演化信息的单晶镍基高温合金塑性行为的本构模型.该本构模型充分考虑了镍基单晶合金中位错在基体相和沉淀增强相中的多种演化机制,例如,基体位错八面体滑移、立方滑移、位错攀移、交滑移、位错弓出、位错切过沉淀增强相以及位错Kear-Wilsdolf(K-W)锁形成与解锁等.在商用有限元软件ABAQUS的框架下,编制了UMAT用户材料子程序.利用该用户子程序,对单晶和多晶镍基高温合金在不同温度、不同加载方向下的单调塑性、循环塑性、蠕变等典型行为进行了计算模拟.结果表明:该晶体塑性本构模型能"统一地"刻画镍基高温合金在不同温度、不同方向下的多种变形行为,并与实验结果具有良好的一致性.     

Quantifying dislocation microstructure evolution and cyclic hardening in fatigued face-centered cubic single crystals

Discrete dislocation dynamics simulations were performed to investigate the dislocation microstructure evolution and cyclic hardening during the early stages of fatigue loading in nickel single crystals. The effects of the crystal size and initial dislocation densities on both the mechanical response and the evolution of dislocation microstructure were quantified. Crystals having an initial dislocation density of 10¹²  m⁻² and diameter less than $2.0\mathrm{\mu }\mathrm{m}$ do not show any dislocation density multiplication or cyclic hardening. In contrast, crystals having the same initial dislocation density and diameters larger than $2.0\mathrm{\mu }\mathrm{m}$ show a significant dislocation density accumulation in the form of dislocation cell-like structures, even after only a few number of loading cycles. This dislocation density accumulation was also accompanied by considerable cyclic hardening. The dislocation cell size and its wall thickness increase with increasing crystal size. With increasing dislocation density the critical crystal size, at which dislocation cell-structures form, decreases. The information theoretic entropy is utilized as a metric to quantify the extent of dislocation patterning and the formation and evolution of dislocation cell structures over time. Cross-slip was found to play a dominant role in the dislocation cell-structure formation. Further insights on the mechanisms contributing to the observed behavior are presented and discussed.     

A multi-scale computational model of crystal plasticity at submicron-to-nanometer scales

Plastic flow in crystal at submicron-to-nanometer scales involves many new interesting problems. In this paper, a unified computational model which directly combines 3D discrete dislocation dynamics (DDD) and continuum mechanics is developed to investigate the plastic behaviors at these scales. In this model, the discrete dislocation plasticity in a finite crystal is solved under a completed continuum mechanics framework: (1) an initial internal stress field is introduced to represent the preexisting stationary dislocations in the crystal; (2) the external boundary condition is handled by finite element method spontaneously; and (3) the constitutive relationship is based on the finite deformation theory of crystal plasticity, but the discrete plastic strains induced by the slip of the newly nucleated or propagating dislocations are calculated by dislocation dynamics methodology instead of phenomenological evolution equations used in conventional crystal plasticity. These discrete plastic strains are then localized to the continuum material points by a Burgers vector density function proposed by us. Various processes, such as loop dislocation evolution, dislocation junction formation etc., are simulated to verify the reliability of this computational model. Specifically, a uniaxial compression test for micro-pillars of Cu is simulated by this model to investigate the ‘dislocation starvation hardening’ observed in the recent experiment.     

Polychromatic microdiffraction analysis of defect self-organization in shock deformed single crystals

A spatially resolved X-ray diffraction method – with a submicron 3D resolution together with SEM and OIM analysis are applied to understand the arrangements of voids, geometrically necessary dislocations and strain gradient distributions in samples of Al (1 2 3) and Cu (0 0 1) single crystals shocked to incipient spallation fracture. We describe how geometrically necessary dislocations and the effective strain gradient alter white beam Laue patterns of the shocked materials. Several distinct structural zones are observed at different depths under the impact surface. The density of geometrically necessary dislocations (GNDs) is extremely high near the impact and back surface of the shock recovered crystals. The spall region is characterized by a large density of mesoscale voids and GNDs. The spall region is separated from the impact and back surfaces by compressed regions with high total dislocation density but lower GNDs density. Self-organization of shear bands is observed in the shock recovered Cu single crystal.     

纳米晶体铝在单轴拉伸下的位错堆积对微裂纹扩展的影响

摘要：针对纳米晶体材料,研究了单轴拉伸载荷作用下纳米晶体铝中的裂纹与裂纹尖端发射的位错所形成的滑移面之间的相互作用。通过分布位错法,将裂纹和滑移面等效为均匀分布的连续位错,获得了裂纹面上应力场。并引入裂纹尖端的无位错区,研究了裂纹尖端无位错区对微裂纹的萌生和主裂扩展之间的影响。结果表明,不考虑裂纹尖端无位错区时,裂纹长度较短,会先在晶界处形成微裂纹,主裂纹较长时,主裂纹会直接穿晶扩展。滑移面与裂纹尖端夹角较大时,会增加裂纹尖端发射的位错个数,从而抑制主裂纹的扩展。考虑裂纹尖端无位错区时,无位错区先于晶界处出现微裂纹,通过主裂纹与微裂纹之间位错的相互发射,导致裂纹与尖端处微裂纹汇合,有效加速了主裂纹的扩展。     

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

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

A cyclic plasticity model for single crystals

The Armstrong–Frederick type kinematic hardening rule was invoked to capture the Bauschinger effect of the cyclic plastic deformation of a single crystal. The yield criterion and flow rule were built on individual slip systems. Material memory was introduced to describe strain range dependent cyclic hardening. The experimental results of copper single crystals were used to evaluate the cyclic plasticity model. It was found that the model was able to accurately describe the cyclic plastic deformation and properly reflect the dislocation substructure evolution. The well-known three distinctive regimes in the cyclic stress–strain curve of the copper single crystals oriented for single slip can be reproduced by using the model. The model can predict the enhanced hardening for crystals oriented for multislip, showing the model's ability to describe anisotropic cyclic plasticity. For a given loading history, the model was able to capture not only the saturated stress–strain response but also the detailed transient stress–strain evolution. The model was used to predict the cyclic plasticity under a high–low loading sequence. Both the stress–strain responses and the microstructural evolution can be appropriately described through the slip system activation.     

Dislocation pile-ups in bicrystals within continuum dislocation theory

Within continuum dislocation theory the plastic deformation of bicrystals under a mixed deformation of plane constrained uniaxial extension and shear is investigated with regard to the nucleation of dislocations and the dislocation pile-up near the phase boundaries of a model bicrystal with one active slip system within each single crystal. For plane uniaxial extension, we present a closed-form analytical solution for the evolution of the plastic distortion and of the dislocation network in the case of symmetric slip planes (i.e. for twins), which exhibits an energetic as well as a dissipative threshold for the dislocation nucleation. The general solution for non-symmetric slip systems is obtained numerically. For a combined deformation of extension and shear, we analyze the possibility of linearly superposing results obtained for both loading cases independently. All solutions presented in this paper also display the Bauschinger effect of translational work hardening and a size effect typical to problems of crystal plasticity.     

Fatigue-induced reversible/irreversible structural-transformations in a Ni-based superalloy

Cyclic loading and the subsequent fatigue-induced structural transformations have been investigated with in-situ neutron diffraction and thermal characterization for a single-phase, polycrystal nickel-based alloy. The lattice-strain evolution is compared with bulk parameters, such as the applied stress and thermal response as a function of the fatigue cycles. In-situ neutron-diffraction and temperature-evolution measurements identify the development of different stages in the fatigue-induced structural transformations, such as bulk hardening, softening, and eventual saturation. An increase in the dislocation density and the formation of planar-patterned dislocation structures are responsible for hardening during the early cycles. With further cyclic loading, the rearrangement of the dislocations results in cyclic softening. A transition is observed during the saturation cycles, which is characterized by the emergence of lattice-strain asymmetry in the loading and transverse directions. The dislocation density and dislocations-wall spacing are determined with diffraction-profile analyses and complemented by transmission-electron microscopy. The thermal behavior of the sample during deformation correlates with corresponding in-situ observation of the time-dependent dislocation structure. An anomaly during saturation cycles is believed to arise from dislocation self-organization – possibly during the formation of microcracks.     

Nonlocal continuum crystal plasticity with internal residual stresses

We derive a three-dimensional constitutive theory accounting for length-scale dependent internal residual stresses in crystalline materials that develop due to a non-homogeneous spatial distribution of the excess dislocation (edge and screw) density. The second-order internal stress tensor is derived using the Beltrami stress function tensor φ that is related to the Nye dislocation density tensor. The formulation is derived explicitly in a three-dimensional continuum setting for elastically isotropic materials. The internal stresses appear as additional resolved shear stresses in the crystallographic visco-plastic constitutive law for individual slip systems. Using this formulation, we investigate two boundary value problems involving single crystals under symmetric double slip. In the first problem, the response of a geometrically imperfect specimen subjected to monotonic and cyclic loading is investigated. The internal stresses affect the overall strengthening and hardening under monotonic loading, which is mediated by the severity of initial imperfections. Such imperfections are common in miniaturized specimens in the form of tapered surfaces, fillets, fabrication induced damage, etc., which may produce strong gradients in an otherwise nominally homogeneous loading condition. Under cyclic loading the asymmetry in the tensile and compressive strengths due to this internal stress is also strongly influenced by the degree of imperfection. In the second example, we consider simple shear of a single crystalline lamella from a layered specimen. The lamella exhibits strengthening with decreasing thickness and increasing lattice incompatibility with shearing direction. However, as the thickness to internal length-scale ratio becomes small the strengthening saturates due to the saturation of the internal stress.Finally, we present the extension of this approach for crystalline materials exhibiting elastic anisotropy, which essentially depends on the appropriate Green function within φ.     

DP钢板应变路径多步演变下的弹塑性力学行为

针对DP高强双相钢板在复杂载荷作用下的弹塑性力学特征,提出利用三步拉伸力学实验,对比分析单轴循环加载和非等轴加载下材料的各向异性硬化、永久软化和弹性模量衰减特性等力学行为,揭示应变路径多步演变下的弹塑性力学特性.研究结果表明:材料再加载初期的瞬态行为与应变路径有关,在初期瞬态阶段显示出明显的各向异性,且再加载角度、预应变量以及马氏体含量的增加都会加剧两种DP钢过渡期后的应力软化行为.同时对不同加载下DP钢弹性模量的演化进行了对比,发现弹性模量的衰减程度也随着再加载角度,预应变量和马氏体含量的增加而增加,弹性模量衰减程度在多步非等轴加载比单轴循环载荷作用下表现更明显.     

A statistical analysis of the elastic distortion and dislocation density fields in deformed crystals

The statistical properties of the elastic distortion fields of dislocations in deforming crystals are investigated using the method of discrete dislocation dynamics to simulate dislocation structures and dislocation density evolution under tensile loading. Probability distribution functions (PDF) and pair correlation functions (PCF) of the simulated internal elastic strains and lattice rotations are generated for tensile strain levels up to 0.85%. The PDFs of simulated lattice rotation are compared with sub-micrometer resolution three-dimensional X-ray microscopy measurements of rotation magnitudes and deformation length scales in 1.0% and 2.3% compression strained Cu single crystals to explore the linkage between experiment and the theoretical analysis. The statistical properties of the deformation simulations are analyzed through determinations of the Nye and Kröner dislocation density tensors. The significance of the magnitudes and the length scales of the elastic strain and the rotation parts of dislocation density tensors are demonstrated, and their relevance to understanding the fundamental aspects of deformation is discussed.     

Size effects under homogeneous deformation of single crystals: A discrete dislocation analysis

Mechanism-based discrete dislocation plasticity is used to investigate the effect of size on micron scale crystal plasticity under conditions of macroscopically homogeneous deformation. Long-range interactions among dislocations are naturally incorporated through elasticity. Constitutive rules are used which account for key short-range dislocation interactions. These include junction formation and dynamic source and obstacle creation. Two-dimensional calculations are carried out which can handle high dislocation densities and large strains up to 0.1. The focus is laid on the effect of dimensional constraints on plastic flow and hardening processes. Specimen dimensions ranging from hundreds of nanometers to tens of microns are considered. Our findings show a strong size-dependence of flow strength and work-hardening rate at the micron scale. Taylor-like hardening is shown to be insufficient as a rationale for the flow stress scaling with specimen dimensions. The predicted size effect is associated with the emergence, at sufficient resolution, of a signed dislocation density. Heuristic correlations between macroscopic flow stress and macroscopic measures of dislocation density are sought. Most accurate among those is a correlation based on two state variables: the total dislocation density and an effective, scale-dependent measure of signed density.     

Study of size effects in thin films by means of a crystal plasticity theory based on DiFT

In a recent publication, we derived the mesoscale continuum theory of plasticity for multiple-slip systems of parallel edge dislocations, motivated by the statistical-based nonlocal continuum crystal plasticity theory for single-glide given by Yefimov et al. [2004b. A comparison of a statistical-mechanics based plasticity model with discrete dislocation plasticity simulations. J. Mech. Phys. Solids 52, 279-300]. In this dislocation field theory (DiFT) the transport equations for both the total dislocation density and geometrically necessary dislocation (GND) density on each slip system were obtained from the Peach-Koehler interactions through both single and pair dislocation correlations. The effect of pair correlation interactions manifested itself in the form of a back stress in addition to the external shear and the self-consistent internal stress. We here present the study of size effects in single crystalline thin films with symmetric double slip using the novel continuum theory. Two boundary value problems are analyzed: (1) stress relaxation in thin films on substrates subject to thermal loading, and (2) simple shear in constrained films. In these problems, earlier discrete dislocation simulations had shown that size effects are born out of layers of dislocations developing near constrained interfaces. These boundary layers depend on slip orientations and applied loading but are insensitive to the film thickness. We investigate the stress response to changes in controlled parameters in both problems. Comparisons with previous discrete dislocation simulations are discussed.