Multiscale dislocation dynamics simulations of shock-induced plasticity in small volumes

Multiscale dislocation dynamics plasticity (MDDP) was used to investigate shock-induced deformation in monocrystalline copper. In order to enhance the numerical simulations, a periodic boundary condition was implemented in the continuum finite element (FE) scale so that the uniaxial compression of shocks could be attained. Additionally, lattice rotation was accounted for by modifying the dislocation dynamics (DD) code to update the dislocations’ slip systems. The dislocation microstructures were examined in detail and a mechanism of microband formation is proposed for single- and multiple-slip deformation. The simulation results show that lattice rotation enhances microband formation in single slip by locally reorienting the slip plane. It is also illustrated that both confined and periodic boundary conditions can be used to achieve uniaxial compression; however, a periodic boundary condition yields a disturbed wave profile due to edge effects. Moreover, the boundary conditions and the loading rise time show no significant effects on shock–dislocations interaction and the resulting microstructures. MDDP results of high strain rate calculations are also compared with the predictions of the Armstrong–Zerilli model of dislocation generation and movement. This work confirms that the effect of resident dislocations on the strain rate can be neglected when a homogeneous nucleation mechanism is included.     

Effect of Small Preliminary Deformation on the Evolution of Elastoplastic Waves of Shock Compression in Annealed VT1-0 Titanium

The evolution of an elastoplastic waves of shock compression in VT1-0 titanium in the as-annealed state and after preliminary compression is measured. A preliminary strain of 0.6% and the related increase in the dislocation density are found to change the deformation kinetics radically and to decrease the Hugoniot elastic limit. An increase in the preliminary strain from 0.6% to 5.2% only weakly changes the Hugoniot elastic limit and the compression rate in the plastic shock wave. The measurement results are used to plot the strain rate versus the stress at the initial stage of high-rate deformation, and the experimental results are interpreted in terms of dislocation dynamics.     

Effects of void shape and orientation on the elastoplastic properties of spheroidally voided single-crystal and nanotwinned copper

ABSTRACT

Void nucleation, growth and coalescence have been identified as the leading cause of ductile damage in metallic materials. To understand the underlying deformation and damage mechanisms, extensive theoretical, experimental and simulation efforts have been attempted on spherically voided metals. In this work, molecular dynamics simulations are performed to analyze the uniaxial straining deformation behaviours of both single-crystal and nanotwinned copper materials embedded with a preexisting spheroidal void. The coupling effects among twin boundary, spheroidal void aspect ratio and orientation on unidirectional elastoplastic behaviours are systematically examined. The dislocation-induced plastic deformation mechanism is also examined and compared with the one due to a perfectly spherical cavity. Simulation results show that elastic modulus increases with both spheroidal void aspect ratio and orientation. So do the yield stress, the first peak stress and the plasticity index. Another peak stress exists for most cases, except for a prolate void embedded in nanotwinned specimens. The slope between peak stresses decreases with both the spheroidal aspect ratio and orientation. The incorporation of a twin boundary results in lower elastic modulus, higher yield strength and smaller plasticity index. For an oblate void, the twin boundary gives rise to more severe strain softening behaviour. The dislocation extraction algorithm illustrates that the continuous nucleation, propagation and reaction of dislocations emanated from both the void front and twin boundary are responsible for the ductile damage of spheroidally voided crystals. The lower dislocation densities found in nanotwinned specimens indicate the desired suppression effects of twin boundary on dislocation activities.     

Modelling and simulation of dynamic recrystallisation based on multi-phase-field and dislocation-based crystal plasticity models

ABSTRACT

The mechanical properties of metals used as structural materials are significantly affected by hot (or warm) plastic working. Therefore, it is industrially important to predict the microscopic behaviour of materials in the deformation process during heat treatment. In this process, a number of nuclei are generated in the vicinity of grain boundaries owing to thermal fluctuation or the coalescence of subgrains, and dynamic recrystallisation (DRX) occurs along with the deformation. In this paper, we develop a DRX model by coupling a dislocation-based crystal plasticity model and a multi-phase-field (MPF) model through the dislocation density. Then, the temperature dependence of the hardening tendency in the recrystallisation process is introduced into the DRX model. A multiphysics simulation for pure Ni is conducted, and then the validity of the DRX model is investigated by comparing the numerical results of microstructure formation and the nominal stress–strain curve during DRX with experimental results. The obtained results indicate that in the process of DRX, nucleation and grain growth occur mainly around grain boundaries with high dislocation density. As deformation progresses, new dislocations pile up and subsequent nucleation occurs in the recrystallised grains. The influence of such microstructural evolution appears as oscillation in the stress–strain curve. From the stress–strain curves, the temperature dependence in DRX is observed mainly in terms of the yield stress, the hardening ratio, and the change in the hardening tendency after nucleation occurs.     

Nucleation of misfit dislocations and plastic deformation in core/shell nanowires

During fabrication of metal nanowires, an oxide layer (shell) that surrounds the metal (core) may form. Such an oxide-covered nanowire can be viewed as a cylindrical core/shell nanostructure, possessing a crystal lattice mismatch between the core and shell. Experimental evidence has shown that, in response to this mismatch, mechanical stresses induce plastic deformation in the shell and misfit dislocations nucleate at the core/shell interface. As a result, the mechanical, electrical and optoelectronic properties of the nanowire are affected. It is therefore essential to be able to predict the critical conditions at which misfit dislocation nucleation at the nanowire interface takes place and the critical applied load at which the interface begins deforming plastically. Two approaches are explored in order to analyze the stress relaxation processes in these oxide-covered nanowires: (i) energy considerations are carried out within a classical elasticity framework to predict the critical radii (of the core and shell) at which dislocation nucleation takes place at the nanowire interface; (ii) a strain gradient plasticity approach is applied to estimate the flow stress at which the interface will begin deforming plastically (this stress is termed “interfacial-yield” stress). The interfacial-yield stress, predicted by gradient plasticity, depends, among other material parameters, on the radii of the core and shell. Both approaches demonstrate how the geometric parameters of nanowires can be calibrated so as to avoid undesirable plastic deformation; in particular, method (i) can give the radii values that prevent misfit dislocation formation, whereas method (ii) can provide, for particular radii values, the critical stress at which interface deformation initiates.     

纳米多晶铜中冲击波阵面的分子动力学研究

冲击波阵面反映材料在冲击压缩下的弹塑性变形行为以及屈服强度、应变率条件等宏观量, 还与冲击压缩后的强度变化联系. 本文使用分子动力学方法, 模拟研究了冲击压缩下纳米多晶铜中的动态塑性变形过程, 考察了冲击波阵面和弹塑性机理对晶界存在的依赖, 并与纳米多晶铝的冲击压缩进行了比较. 研究发现: 相比晶界对纳米多晶铝的贡献而言, 纳米多晶铜中晶界对冲击波阵面宽度的影响较小; 并且其塑性变形机理主要以不全位错的发射和传播为主, 很少观察到全位错和形变孪晶的出现. 模拟还发现纳米多晶铜的冲击波阵面宽度随着冲击应力的增加而减小, 并得到了冲击波阵面宽度与冲击应力之间的定量反比关系, 该定量关系与他人纳米多晶铜模拟结果相近, 而与粗晶铜的冲击压缩实验结果相差较大.     

Selective excitation of two-wave structure depending on crystal orientation under shock compression

Single crystals subjected to shock compression exhibit responses with distinct two-wave structures for certain crystal orientations. However, little is known to date regarding how the shock response depends on crystal orientation, and especially why the two-wave structure depends on the crystal orientation. In this work, molecular dynamics simulations of shock compressions in copper single crystals are performed to investigate the orientation dependence of shock responses and the corresponding deformation mechanisms. Four copper single crystals with [001], [011], [012], and [123] crystal orientations along the depth direction are investigated. The [011], [012], and [123] crystal orientations of copper single crystals show distinct two-wave structures in their shock responses, while such a two-wave structure in the shock response is not seen for those orientations having a [001] crystal orientation. The potential causes are analyzed by considering the propagation velocities of both elastic and plastic waves. We develop a technique for identifying twin structures in face-centered cubic crystals and this technique can effectively identify the twin structure. The morphology of shock-induced defects(e.g., dislocations and twins) shows the significant dependence of crystal orientation and the mechanisms behind these are discussed in detail. Finally, the Johnson-Cook constitutive model describing dynamic deformations at high temperatures and high strain rates is used to analyze the relationships between the shock responses and microscopic defects. The predictions of the Johnson-Cook constitutive model are consistent with the results of the molecular dynamics simulations.     

Effect of hydrostatic pressure on the plasticity of copper

A new technique is designed to plot stress-strain curves during torsion at hydrostatic pressure up to 250 MPa. It is used to study the effect of pressure on the torsional strain to failure of copper. The experimental results demonstrate that the plasticity of the material increases in the pressure range 0–150 MPa and remains constant in the range 150–250 MPa. This reaction of the material is likely to be related to the fact that pressure can increase the dislocation density required for plastic flow.     

晶体塑性有限元方法研究辐照对多晶铜力学性能的影响

将辐照硬化理论与晶体塑性理论结合, 运用ABAQUS有限元分析软件模拟辐照后多晶铜的拉伸过程。分析辐照效应对材料屈服强度、硬化过程、晶体变形等力学性能的影响, 研究位错密度的演化及空间分布规律。数值模拟表明: 辐照效应提高多晶铜的屈服应力, 影响不同阶段的硬化和软化现象; 辐照剂量增大导致位错密度增殖总体变缓, 空间不均匀度增大; 晶体的塑性变形及晶体转动也受到辐照的影响, 在相同的应变条件下, 辐照剂量越大, 晶体塑性变形程度越小, 塑性变形分布不均匀度变大, 同时晶体转动程度及转动角离散度增大。     

非一维应变冲击加载下高纯铜初始层裂行为

提出了一种锥形靶层裂实验新方法,开展非一维应变冲击条件下高纯铜初始层裂行为实验研究,讨论了锥形靶内部损伤分布特征及其与自由面速度典型特征之间的内禀关系.结果显示:1)初始层裂的锥形靶内部出现了连续损伤区,损伤区扩展方向与锥面平行,从锥底到锥顶呈现了不同的损伤状态,从微孔洞独立长大到局部聚集,最后形成宏观裂纹,这种损伤状态分布特征归因于锥形靶内部拉伸应力幅值和持续时间的空间演化;2)通过锥形靶横截面损伤度定量统计分析,揭示损伤演化早期的微孔洞成核与早期长大过程是随机的,而损伤演化后期的微孔洞聚集过程具有显著的局域化特征;3)不同位置处实测的自由面法向粒子速度剖面呈现出典型的层裂Pull-back信号,但是通过与内部损伤分布特征对比,揭示基于Pull-back速度获得高纯铜层裂强度本质是微孔洞成核阈值应力,Pull-back回跳速度斜率反映了损伤演化速率,Pull-back回跳幅值与损伤度引起的应力松弛密切相关.     

Supersonic dislocation stability and nano-twin formation at high strain rate

Improved understanding of the plastic deformation of metals during high-strain-rate shock loading is key to predicting their resulting material properties. This paper presents the results of molecular-dynamics simulations which address two fundamental questions related to materials deformation: the stability of supersonic dislocations and the mechanism of nano-twin formation. The results show that aluminium plastically deforms by the subsonic motion of edge dislocations when subjected to applied shear stresses of up to 600?MPa. Although higher applied stresses initially drive transonic dislocations, this motion is transient, and the dislocations decelerate to a sustained subsonic saturation velocity. Slowing of the transonic dislocation is controlled by the interaction with excited Rayleigh waves. 800?MPa marks a critical shear stress at which dislocation glide gives way to nano-twin formation via the homogeneous nucleation of Shockley partial dislocation dipoles. At still higher applied stresses, additional dislocation dipole nucleation produces a mid-stacking fault transformation of the twinned material.     

单晶Cu冲击加载及卸载下塑性行为的微观模拟

使用分子动力学方法对室温下单晶铜沿［001］和［111］方向冲击加载及卸载下的塑性行为进行了模拟,得到了Hugoniot关系以及冲击熔化压力,与实验基本符合. 加载过程中,较高的初始温度有利于位错的形核与发展. 通过对冲击波在自由表面卸载过程的模拟和分析发现：卸载过程呈现“准弹性卸载行为”;沿［001］方向卸载后大量不全位错环与堆积层错消失,而沿［111］方向卸载后只有少量层错消失,部分层错甚至会发展扩大. 关键词： 分子动力学 冲击波 塑性     

A composite model of the plastic flow of amorphous covalent materials

A model based on the data available in the literature on the computer simulation of amorphous silicon has been proposed for describing the specific features of the plastic flow of amorphous covalent materials. The mechanism of plastic deformation involves homogeneous nucleation and growth of inclusions of a liquidlike phase under external shear stress. Such inclusions experience plastic shear, which is modeled by glide dislocation loops. The energy changes associated with the nucleation of these inclusions at room and increased temperatures have been calculated. The critical stress has been found, at which the barrierless nucleation of inclusions becomes possible. It has been shown that this stress decreases with an increase in temperature. According to the calculations, the heterogeneous (homogeneous) plastic flow of an amorphous material should be expected at relatively low (high) temperatures. Above the critical stress, the homogeneous flow is gradually replaced by the heterogeneous flow.     

The effect of the strain amplitude changes on the stress amplitude and resistivity of torsionally fatigued copper

Polycrystalline copper wires were cyclically strained in torsion in liquid nitrogen and the effect of the strain amplitude changes on the stress amplitude and resistivity was followed. It was found that both stress amplitude and resistivity are able to increase or decrease in dependence on the applied plastic strain amplitude. When the strain amplitude was decreased, the reversibility was not complete. By means of intermediate annealing it was found that both the dislocation density and the point-defect concentration follow the changes of the plastic strain amplitude.     

Localization of plastic strain under shock-wave loading of titanium alloy with a crack

The effect of shock-wave loading with an amplitude of 4 GPa and a duration of 10 μs on the localization of plastic flow and the dynamic fracture of VT-1 cast titanium alloy containing a disk-shaped crack is considered. It is found that the crack size decreases in the direction of shock wave propagation, the dislocation density grows and adiabatic shear bands appear near the crack, and secondary cracks originate in the areas of localized flow near the crack edges. The strain, strain rate, dislocation density, rate of dislocation generation, and dislocation velocity near and away from the healing crack are estimated.     

Dislocation Plasticity Effects on Interfacial Fracture

Analyses are reviewed where plastic flow in the vicinity of an interfacial crack is represented in terms of the nucleation and glide of discrete dislocations. Attention is confined to cracks along a metal-ceramic interface, with the ceramic idealized as being rigid. Both monotonic and fatigue loading are considered. The main focus is on the stress and deformation fields near the crack tip predicted by discrete dislocation plasticity, in comparison with those obtained from conventional continuum plasticity theory. The role that discrete dislocation plasticity can play in interpreting interface fracture properties in the presence of plastic flow is discussed.     

高应变率压缩下纳米孔洞对金属铝塑性变形的影响研究

本文利用分子动力学模拟方法研究了含纳米孔洞金属铝在[110]晶向高应变率单轴压缩下弹塑性变形的微观过程. 对比单孔洞和完整单晶的模型, 讨论了多孔金属的应力应变关系及其位错发展规律. 研究结果表明, 对于多孔模型的位错积累过程, 位错密度随应变的增加可大致分为两个线性阶段. 由同一个孔洞生成的位错在相互靠近过程中, 其滑移速度越来越小; 随着位错继续滑移, 源自不同孔洞的位错之间开始交叉相互作用导致应变硬化. 达到流变峰应力之后又由于位错密度增殖速率升高发生软化. 当应变增加到11.8%时, 所有孔洞几乎完全坍缩, 并观察到在此过程中有棱位错生成.     

Collective properties of defects,multiscale plasticity,and shock induced phenomena in solids

A statistically based approach is developed for the construction of constitutive equations that provides linkages between defect-induced mechanisms of structural relaxation, thermally activated plastic flow, and material response to extreme loading conditions. The collective properties of defects have been studied to establish the interaction of multiscale defect dynamics and plastic flow, and to explain the mechanisms leading to the universal self-similar structure of shock wave fronts. Pn explanation for structural universality of the steady-state plastic shock front (the four power law) and the self-similarity of shock wave profiles under reloading (unloading) is proposed. Structural characterization under transition from thermally activated dislocation glide to nonlinear dislocation drag effects is developed in terms of scaling invariants (effective temperatures) related to mesodefect induced morphology formed during the different stages of plastic deformation.     

Evaluating the effects of loading parameters on single-crystal slip in tantalum using molecular mechanics

This study is aimed at developing a physics-based crystal plasticity finite element model for body-centred cubic (BCC) metals, through the introduction of atomic-level deformation information from molecular dynamics (MD) investigations of dislocation motion at the onset of plastic flow. In this study, three critical variables governing crystal plasticity mediated by dislocation motion are considered. MD simulations are first performed across a range of finite temperatures up to 600K to quantify the temperature dependence of critical stress required for slip initiation. An important feature of slip in BCC metals is that it is not solely dependent on the Schmid law measure of resolved shear stress, commonly employed in crystal plasticity models. The configuration of a screw dislocation and its subsequent motion is studied under different load orientations to quantify these non-Schmid effects. Finally, the influence of strain rates on thermal activation is studied by inducing higher stresses during activation at higher applied strain rates. Functional dependence of the critical resolved shear stress on temperature, loading orientation and strain rate is determined from the MD simulation results. The functional forms are derived from the thermal activation mechanisms that govern the plastic behaviour and quantification of relevant deformation variables. The resulting physics-based rate-dependent crystal plasticity model is implemented in a crystal plasticity finite element code. Uniaxial simulations reveal orientation-dependent tension–compression asymmetry of yield that more accurately represents single-crystal experimental results than standard models.