Influence of misfit stresses on dislocation glide in single crystal superalloys: A three-dimensional discrete dislocation dynamics study

In the characteristic $\gamma /\gamma \prime$ microstructure of single crystal superalloys, misfit stresses occur due to a significant lattice mismatch of those two phases. The magnitude of this lattice mismatch depends on the chemical composition of both phases as well as on temperature. Furthermore, the lattice mismatch of γ and $\gamma \prime$ phases can be either positive or negative in sign. The internal stresses caused by such lattice mismatch play a decisive role for the micromechanical processes that lead to the observed macroscopic athermal deformation behavior of these high-temperature alloys. Three-dimensional discrete dislocation dynamics (DDD) simulations are applied to investigate dislocation glide in γ matrix channels and shearing of $\gamma \prime$ precipitates by superdislocations under externally applied uniaxial stresses, by fully taking into account internal misfit stresses. Misfit stress fields are calculated by the fast Fourier transformation (FFT) method and hybridized with DDD simulations. For external loading along the crystallographic [001] direction of the single crystal, it was found that the different internal stress states for negative and positive lattice mismatch result in non-uniform dislocation movement and different dislocation patterns in horizontal and vertical γ matrix channels. Furthermore, positive lattice mismatch produces a lower deformation rate than negative lattice mismatch under the same tensile loading, but for an increasing magnitude of lattice mismatch, the deformation resistance always diminishes. Hence, the best deformation performance is expected to result from alloys with either small positive, or even better, vanishing lattice mismatch between γ and $\gamma \prime$ phase.     

The effect of inhomogeneous plastic deformation on the ductility and fracture behavior of age hardenable aluminum alloys

The role of alloy composition, grain structure, precipitate microstructure, and precipitate dislocation interactions on the plastic deformation characteristics and the resulting fracture behavior of two isotropic Al–Li–Cu–X alloys designated AF/C-458 (1.8 w/o Li) and AF/C-489 (2.1 w/o Li) was examined. Inhomogeneous deformation due to strain localization from the shearing of the δ′ (Al₃Li), θ′ (Al₂Cu), and T₁ (Al₂CuLi) precipitates lead to fine and coarse planar slip for the AF/C-458 and AF/C-489 alloys, respectively. The intensity of this planar slip was predicted through slip intensity calculations using precipitate density measurements, dislocation particle interactions, and grain boundary misorientation-slip continuity statistics. The slip intensity predictions were corroborated through atomic force microscopy (AFM) measured slip height offsets on the polished surface of single aged and 2% plastically strained tensile samples. Our results suggest that the low ductility of AF/C-489 in comparison to AF/C-458 is primarily due to the much larger slip lengths, i.e. grain size, which increased the strain localization and stress concentrations on grain boundaries, thus promoting low-energy intergranular fracture.     

A computational method for dislocation-precipitate interaction

A new computational method for the elastic interaction between dislocations and precipitates is developed and applied to the solution of problems involving dislocation cutting and looping around precipitates. Based on the superposition principle, the solution to the dislocation-precipitate interaction problem is obtained as the sum of two solutions: (1) a dislocation problem with image stresses from interfaces between the dislocation and the precipitate, and (2) a correction solution for the elastic problem of a precipitate with an initial strain distribution. The current development is based on a combination of the parametric dislocation dynamics (PDD) and the boundary element method (BEM) with volume integrals.The method allows us to calculate the stress field both inside and outside precipitates of elastic moduli different from the matrix, and that may have initial coherency strain fields. The numerical results of the present method show good convergence and high accuracy when compared to a known analytical solution, and they are also in good agreement with molecular dynamics (MD) simulations. Sheared copper precipitates (2.5 nm in diameter) are shown to lose some of their resistance to dislocation motion after they are cut by leading dislocations in a pileup. Successive cutting of precipitates by the passage of a dislocation pileup reduces the resistance to about half its original value, when the number of dislocations in the pileup exceeds about 10. The transition from the shearable precipitate regime to the Orowan looping regime occurs for precipitate-to-matrix elastic modulus ratios above approximately 3-4, with some dependence on the precipitate size. The effects of precipitate size, spacing, and elastic modulus mismatch with the host matrix on the critical shear stress (CSS) to dislocation motion are presented.     

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.     

An atomistically-informed dislocation dynamics model for the plastic anisotropy and tension-compression asymmetry of BCC metals

Atomistic simulations have shown that a screw dislocation in body-centered cubic (BCC) metals has a complex non-planar atomic core structure. The configuration of this core controls their motion and is affected not only by the usual resolved shear stress on the dislocation, but also by non-driving stress components. Consequences of the latter are referred to as non-Schmid effects. These atomic and micro-scale effects are the reason slip characteristics in deforming single and polycrystalline BCC metals are extremely sensitive to the direction and sense of the applied load. In this paper, we develop a three-dimensional discrete dislocation dynamics (DD) simulation model to understand the relationship between individual dislocation glide behavior and macro-scale plastic slip behavior in single crystal BCC Ta. For the first time, it is shown that non-Schmid effects on screw dislocations of both {110} and {112} slip systems must be implemented into the DD models in order to predict the strong plastic anisotropy and tension-compression asymmetry experimentally observed in the stress-strain curves of single crystal Ta. Incorporation of fundamental atomistic information is critical for developing a physics-based, predictive meso-scale DD simulation tool that can connect length/time scales and investigate the underlying mechanisms governing the deformation of BCC metals.     

Finite element analysis of grain-by-grain deformation by crystal plasticity with couple stress

Rigid–plastic crystal plasticity with the rate-sensitive constitutive behavior of a slip system has been formulated within the framework of a two-dimensional finite element method to predict the grain-by-grain deformation of single- and polycrystalline FCC metals. For that purpose, individual grains are represented by several numbers of finite elements to describe the sub-grain deformation behavior, and couple stress has been introduced into the equilibrium equation to be able to describe the size effect as well as to prevent mesh-dependent predictions. A modified virtual work-rate principle with an approximate interface constraint has been suggested to use a C ⁰-continuous element in the finite element implementation, and the couple stress work-rate has been formulated on the basis of an assumed constitutive behavior. Simulated plane-strain compressions of a single crystal cube show that the shearing and the deformation load are closely related to the imbedded lattice orientation of the crystal grain, and that the sub-grain deformation and the load magnitude can be controlled by the couple stress hardening. It is also confirmed that almost the same predictions are obtained for different mesh systems by considering the couple stress hardening. Simulated plane-strain compressions of a bi-crystal show considerably curved grain-by-grain surface profiles after large reduction for several combinations of the imbedded lattice orientation. The high couple stress hardening predicted around grain boundaries is supposed to be related to the grain size effect. It is also supposed that consideration of couple stress is necessary to predict the sub-grain or the grain-by-grain deformation, and the couple stress hardening may be used to describe the state of microstructures in grain.     

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.     

Transmission electron microscopy investigation of dislocation slip during superelastic cycling of Ni-Ti wires

Superelastic deformation of thin Ni-Ti wires containing various nanograined microstructures was investigated by tensile cyclic loading with in situ evaluation of electric resistivity. Defects created by the superelastic cycling in these wires were analyzed by transmission electron microscopy. The role of dislocation slip in superelastic deformation is discussed. Ni-Ti wires having finest microstructures (grain diameter <100 nm) are highly resistant against dislocation slip, while those with fully recrystallized microstructure and grain size exceeding 200 nm are prone to dislocation slip. The density of the observed dislocation defects increases significantly with increasing grain size. The upper plateau stress of the superelastic stress-strain curves is largely grain size independent from 10 up to 1000 nm. It is hence claimed that the Hall-Petch relationship fails for the stress-induced martensitic transformation in this grain size range. It is proposed that dislocation slip taking place during superelastic cycling is responsible for the accumulated irreversible strains, cyclic instability and degradation of functional properties. No residual martensite phase was found in the microstructures of superelastically cycled wires by TEM and results of the in situ electric resistance measurements during straining also indirectly suggest that none or very little martensite phase remains in the studied cycled superelastic wires after unloading. The accumulation of dislocation defects, however, does not prevent the superelasticity. It only affects the shape of the stress-strain response, makes it unstable upon cycling and changes the deformation mode from localized to homogeneous. The activity of dislocation slip during superelastic deformation of Ni-Ti increases with increasing test temperature and ultimately destroys the superelasticity as the plateau stress approaches the yield stress for slip. Deformation twins in the austenite phase ({1 1 4} compound twins) were frequently found in cycled wires having largest grain size. It is proposed that they formed in the highly deformed B19′ martensite phase during forward loading and are retained in austenite after unloading. Such twinning would represent an additional deformation mechanism of Ni-Ti yielding residual irrecoverable strains.     

Plane constrained shear of single crystal strip with two active slip systems

Within continuum dislocation theory the plane constrained shear of a single crystal strip with two active slip systems is considered. An analytical solution is found for symmetric double slip which exhibits the energetic and dissipative thresholds for dislocation nucleation, the Bauschinger translational work hardening, and the size effects. Comparison with discrete dislocation simulations shows good agreement between the discrete and continuum approaches. Numerical procedures in the general case of non-symmetric double slip are proposed.     

The role of superdislocations for modeling plastic deformation of lamellar TiAl

The deformation behaviour of lamellar TiAl-alloys is investigated by means of micromechanical finite element simulations. For the simulations, a lower and an upper bound periodic unit cell representing the lamellar substructure are used. Crystallographic slip is described by continuum crystal plasticity. For the simulations, 16 individual ordinary, super, and twinning systems divided into three morphological slip modes are taken into account. The strength and hardening parameters of the slip systems are determined by comparison with experimental data from room temperature compression tests of poly-synthetically twinned (PST) crystals of TiAl. The effect of the yield stress of the superdislocations is investigated by systematic parameter variations. The results show that the general anisotropy of PST-crystals can already be described with a purely morphological classification of the slip systems. Superdislocations cause a slight tension–compression anomaly if the yield stresses are higher than for the ordinary and twinning systems.     

A semi-implicit integration scheme for rate independent finite crystal plasticity

An efficient new numerical integration scheme is presented for rate independent crystal plasticity theory. A key feature of this approach is the ability to identify active slip systems prior to determining their shearing rate. Options are described for various cases of slip system hardening, including self hardening and latent hardening. Alternatives for the constitutive update are explored, including hyperelasticity based on the multiplicative decomposition of the deformation gradient as well as application of the consistency condition in a much more efficient hypoelastic formulation. Several conclusions are drawn concerning the influences of elastic and plastic properties on the activation of slip systems and their subsequent shearing rates. Key among these is the fact that, once activated, shearing rates are independent of the levels of shear flow resistance on the slip systems, provided that the plastic hardening moduli are much less in magnitude than the elastic moduli, as is usually the case. Determination of active slip systems and their shearing rates depend on the degree of elastic anisotropy of the crystal, but not on the magnitude of elastic stiffness.     

A continuum model for initiation and evolution of deformation twinning

Within continuum dislocation theory the plastic deformation of a single crystal with one active slip system under plane-strain constrained shear is investigated. By introducing a twinning shear into the energy of the crystal, we show that in a certain range of straining the formation of deformation twins becomes energetically preferable. An energetic threshold for the onset of twinning is determined. A rough analysis qualitatively describes not only the evolving volume fractions of twins but also their number during straining. Finally, we analyze the evolution of deformation twins and of the dislocation network at non-zero dissipation. We present the corresponding stress-strain hysteresis, the evolution of the plastic distortion, the twin volume fractions and the dislocation densities.     

A dislocation-based constitutive law for pure Zr including temperature effects

In this work, a single crystal constitutive law for multiple slip and twinning modes in single phase hcp materials is developed. For each slip mode, a dislocation population is evolved explicitly as a function of temperature and strain rate through thermally-activated recovery and debris formation and the associated hardening includes stage IV. A stress-based hardening law for twin activation accounts for temperature effects through its interaction with slip dislocations. For model validation against macroscopic measurement, this single crystal law is implemented into a visco-plastic-self-consistent (VPSC) polycrystal model which accounts for texture evolution and contains a subgrain micromechanical model for twin reorientation and morphology. Slip and twinning dislocations interact with the twin boundaries through a directional Hall–Petch mechanism. The model is adjusted to predict the plastic anisotropy of clock-rolled pure Zr for three different deformation paths and at four temperatures ranging from 76 K to 450 K (at a quasi-static rate of 10⁻³ 1/s). The model captures the transition from slip-dominated to twinning-dominated deformation as temperature decreases, and identifies microstructural mechanisms, such as twin nucleation and twin–slip interactions, where future characterization is needed.     

Strain hardening in 2D discrete dislocation dynamics simulations: A new ‘2.5D’ algorithm

The two-dimensional discrete dislocation dynamics (2D DD) method, consisting of parallel straight edge dislocations gliding on independent slip systems in a plane strain model of a crystal, is often used to study complicated boundary value problems in crystal plasticity. However, the absence of truly three dimensional mechanisms such as junction formation means that forest hardening cannot be modeled, unless additional so-called ‘2.5D’ constitutive rules are prescribed for short-range dislocation interactions. Here, results from three dimensional dislocation dynamics (3D DD) simulations in an FCC material are used to define new constitutive rules for short-range interactions and junction formation between dislocations on intersecting slip systems in 2D. The mutual strengthening effect of junctions on preexisting obstacles, such as precipitates or grain boundaries, is also accounted for in the model. The new ‘2.5D’ DD model, with no arbitrary adjustable parameters beyond those obtained from lower scale simulation methods, is shown to predict athermal hardening rates, differences in flow behavior for single and multiple slip, and latent hardening ratios. All these phenomena are well-established in the plasticity of crystals and quantitative results predicted by the model are in good agreement with experimental observations.     

单晶体和双晶体微观层次变形行为的有限元分析

从微观层次上研究金属材料的变形行为,将位错引入到本构关系中,用硬化函数描述材料的硬化规律,考虑了变形的率相关性,采用三维模型用大变形有限单元法对单晶体在单向拉伸载荷和循环载荷作用下的变形行为、双晶体在单向拉伸作用下滑移系的开动进行了模拟计算,得到了与实验一致的计算结果。     

Mechanical annealing under low-amplitude cyclic loading in micropillars

Mechanical annealing has been demonstrated to be an effective method for decreasing the overall dislocation density in submicron single crystal. However, simultaneously significant shape change always unexpectedly happens under extremely high monotonic loading to drive the pre-existing dislocations out of the free surfaces. In the present work, through in situ TEM experiments it is found that cyclic loading with low stress amplitude can drive most dislocations out of the submicron sample with virtually little change of the shape. The underlying dislocation mechanism is revealed by carrying out discrete dislocation dynamic (DDD) simulations. The simulation results indicate that the dislocation density decreases within cycles, while the accumulated plastic strain is small. By comparing the evolution of dislocation junction under monotonic, cyclic and relaxation deformation, the cumulative irreversible slip is found to be the key factor of promoting junction destruction and dislocation annihilation at free surface under low-amplitude cyclic loading condition. By introducing this mechanics into dislocation density evolution equations, the critical conditions for mechanical annealing under cyclic and monotonic loadings are discussed. Low-amplitude cyclic loading which strengthens the single crystal without seriously disturbing the structure has the potential applications in the manufacture of defect-free nano-devices.     

Self-consistent modeling of large plastic deformation, texture and morphology evolution in semi-crystalline polymers

A self-consistent model for semi-crystalline polymers is proposed to study their constitutive behavior, texture and morphology evolution during large plastic deformation. The material is considered as an aggregate of composite inclusions, each representing a stack of crystalline lamellae with their adjacent amorphous layers. The deformation within the inclusions is volume-averaged over the phases. The interlamellar shear is modeled as an additional slip system with a slip direction depending on the inclusion's stress. Hardening of the amorphous phase due to molecular orientation and, eventually, coarse slip, is introduced via Arruda-Boyce hardening law for the corresponding plastic resistance. The morphology evolution is accounted for through the change of shape of the inclusions under the applied deformation gradient. The overall behavior is obtained via a viscoplastic tangent self-consistent scheme. The model is applied to high density polyethylene (HDPE). The stress-strain response, texture and morphology changes are simulated under different modes of straining and compared to experimental data as well as to the predictions of other models.