热加工过程中动态再结晶现象的多相场研究

金属热加工过程中的动态再结晶引起的组织演化难以通过实验实时观察, 本文基于Ginzburg-Landau动力学方程, 构造多相场法与位错密度计算相耦合的物理模型, 模拟了热加工过程中的动态再结晶现象.研究了不同温度和不同应变速率下的动态再结晶过程, 阐述了应力-应变曲线由单峰形式转变为多峰形式的原因.此外, 本文利用多相场法对多阶段变形过程进行了系统模拟, 研究了静态回复对动态再结晶过程的影响, 分析了不同的热加工参数对动态再结晶动力学的影响, 发现在变形间断过程中, 晶粒尺寸不断增大, 较高的变形温度和较低的应变速率可以加速动态再结晶过程.     

On the homogeneous nucleation and propagation of dislocations under shock compression

The dynamic response of crystalline materials subjected to extreme shock compression is not well understood. The interaction between the propagating shock wave and the material’s defect occurs at the sub-nanosecond timescale which makes in situ experimental measurements very challenging. Therefore, computer simulation coupled with theoretical modelling and available experimental data is useful to determine the underlying physics behind shock-induced plasticity. In this work, multiscale dislocation dynamics plasticity (MDDP) calculations are carried out to simulate the mechanical response of copper reported at ultra-high strain rates shock loading. We compare the value of threshold stress for homogeneous nucleation obtained from elastodynamic solution and standard nucleation theory with MDDP predictions for copper single crystals oriented in the [0 0 1]. MDDP homogeneous nucleation simulations are then carried out to investigate several aspects of shock-induced deformation such as; stress profile characteristics, plastic relaxation, dislocation microstructure evolution and temperature rise behind the wave front. The computation results show that the stresses exhibit an elastic overshoot followed by rapid relaxation such that the 1D state of strain is transformed into a 3D state of strain due to plastic flow. We demonstrate that MDDP computations of the dislocation density, peak pressure, dynamics yielding and flow stress are in good agreement with recent experimental findings and compare well with the predictions of several dislocation-based continuum models. MDDP-based models for dislocation density evolution, saturation dislocation density, temperature rise due to plastic work and strain rate hardening are proposed. Additionally, we demonstrated using MDDP computations along with recent experimental reports the breakdown of the fourth power law of Swegle and Grady in the homogeneous nucleation regime.     

Analysis of flow stress and deformation mechanism under hot working of ZK60 magnesium alloy by a new strain-dependent constitutive equation

The present study investigates the variation of flow stress and microstructural evolution with strain for ZK60 magnesium alloy. A new constitutive equation was used to model the flow stress with excellent results. This constitutive analysis and the microstructural studies carried out on strained samples revealed the existence of two different regimes. At temperatures above 300 °C, moderate grain growth and intragranular dislocation activity. Yet, the calculated value of the activation energy and the marked increase in the equivalent strain to fracture indicated grain boundary sliding as a dominant mechanism in this regime of strain rate and temperature, with dislocation motion playing an ancillary role. At lower temperatures, deformation was exclusively governed by dislocation motion, with the extensive occurrence of dynamic recrystallization, which started at low strains, and absence of grain growth.     

双模晶体相场研究形变诱导的多级微结构演化

本文采用双模晶体相场模型,计算了双模二维相图;模拟了形变诱导六角相向正方相转变过程的多级微结构演化,详细分析了位相差、形变方向对位错、晶界、晶体结构、新相形貌的影响规律.模拟结果表明:形变方向影响正方相晶核的形核位置和生长方向,拉伸时正方相优先在变形带上形核,垂直于形变方向长大,而压缩时正方相直接在位错和晶界的能量较高处形核,平行于形变方向长大;位相差对形变诱发晶界甄没过程有显著影响,体现在能量峰上为,小位相差晶界位错的攀滑移和甄没形成一个能量峰,大位相差晶界位错攀滑移和甄没因分阶段完成而不出现明显的能量峰;形变诱导相变过程中各种因素相互作用复杂,是相变与动态再结晶的复合转变.     

Changes in the diffusion properties of nonequilibrium grain boundaries upon recrystallization and superplastic deformation of submicrocrystalline metals and alloys

The effect of an increase in the coefficient of the grain-boundary diffusion upon recrystallization and superplastic deformation of submicrocrystalline (SMC) materials prepared by severe plastic deformation has been studied. It is shown that the coefficient of the grain-boundary diffusion of the SMC materials is dependent on the intensity of the lattice dislocation flow whose value is proportional to the rate of the grain boundary migration upon annealing of SMC metals or the rate of the intragrain deformation under conditions of superplastic deformation of SMC alloys. It is found that, at a high rate of grain boundary migrations and high rates of superplastic deformation, the intensity of the lattice dislocation flow bombarding grain boundaries of SMC materials is higher than the intensity of their diffusion accommodation, which leads to an increase in the coefficient of the grain-boundary diffusion and a decrease in the activation energy. The results of the numerical calculations agree well with the experimental data.     

两相钛合金再结晶退火组织与织构演变的蒙特卡罗模拟

蒙特卡罗(MC)方法被广泛应用于模拟金属材料在退火过程中的静态再结晶行为. 在已有两相材料晶粒长大MC模型基础上, 引入形核阶段, 综合考虑再结晶晶粒吞并形变晶粒和再结晶晶粒竞争长大两种情况, 建立了退火时两相合金再结晶MC模型.结合电子背散射衍射所测 初始晶粒形貌、相成分、晶体学取向及应变储能相对值, 该模型被应用于TC11钛合金退火过程中的微观组织及织构演变模拟.结果表明, 所建模型能够较好体现退火过程中两相晶粒的形核及晶粒长大行为. 与β相相比较, α相具有较低的再结晶速率和较高的晶粒长大速率, 前者主要归结于α相较低的初始应变储能, 后者则体现了该条件下初始组织形貌、分布及两相比例对晶粒长大具有重要影响; 由于非均匀形核的影响, 模拟得到的再结晶速率变化与 假设均匀形核的Johnson-Mehl-Avrami-Kolmogorov 再结晶方程存在明显差异.同时, 两相的基本织构特征在退火过程中无明显变化, 但织构强度增加. 关键词： 两相钛合金 再结晶 蒙特卡罗方法 织构     

Mechanoluminescence induced by elastic and plastic deformation of coloured alkali halide crystals at fixed strain rates

Mechanoluminescence (ML) emission from coloured alkali halide crystals takes place during their elastic and plastic deformation. The ML emission during the elastic deformation occurs due to the mechanical interaction between dislocation segments and F-centres, and the ML emission during the plastic deformation takes place due to the mechanical interaction between the moving dislocations and F-centres. In the elastic region, the ML intensity increases linearly with the strain or deformation time, and in this case, the saturation region could not be observed because of the beginning of the plastic deformation before the start of the saturation in the ML intensity. In the plastic region, initially the ML intensity also increases linearly with the strain or deformation time, and later on, it attains a saturation value for large deformation. When the deformation is stopped, initially the ML intensity decreases at a fast rate; later on, it decreases at a slow rate. The decay time for the fast decrease of the ML intensity gives the relaxation time of dislocation segments or pinning time of the dislocations, and the decay time of the slow decrease of the ML intensity gives the diffusion time of holes in the crystals. The saturation value of the ML intensity increases linearly with the strain rate and also with the density of F-centres in the crystals. Initially, the saturation value of the ML intensity increases with increasing temperature, and for higher temperatures the ML intensity decreases with increasing temperature. Therefore, the ML intensity is optimum for a particular temperature of the crystals. From the ML measurements, the relaxation time of dislocation segments, pinning time of dislocations, diffusion time of holes and the energy gap between the bottom of the acceptor dislocation band and interacting F-centre level can be determined. Expressions derived for the ML induced by elastic and plastic deformation of coloured alkali halide crystals at fixed strain rates indicates that the ML intensity depends on the strain, strain rate, density of colour centres, size of crystals, temperature, luminescence efficiency, etc. A good agreement is found between the theoretical and experimental results.     

Atomistic simulation of the stacking fault energy and grain shape on strain hardening behaviours of FCC nanocrystalline metals

ABSTRACT

Ultra-fine grained copper with nanotwins is found to be both strong and ductile. It is expected that nanocrystalline metals with lamella grains will have strain hardening behaviour. The main unsolved issues on strain hardening behaviour of nanocrystalline metals include the effect of stacking fault energy, grain shape, temperature, strain rate, second phase particles, alloy elements, etc. Strain hardening makes strong nanocrystalline metals ductile. The stacking fault energy effects on the strain hardening behaviour are studied by molecular dynamics simulation to investigate the uniaxial tensile deformation of the layer-grained and equiaxed models for metallic materials at 300?K. The results show that the strain hardening is observed during the plastic deformation of the layer-grained models, while strain softening is found in the equiaxed models. The strain hardening index values of the layer-grained models decrease with the decrease of stacking fault energy, which is attributed to the distinct stacking fault width and dislocation density. Forest dislocations are observed in the layer-grained models due to the high dislocation density. The formation of sessile dislocations, such as Lomer–Cottrell dislocation locks and stair-rod dislocations, causes the strain hardening behaviour. The dislocation density in layer-grained models is higher than that in the equiaxed models. Grain morphology affects dislocation density by influencing the dislocation motion distance in grain interior.     

Shear-Banding Evolution Dynamics during High Temperature Compression of Martensitic Ti-6Al-4V Alloy

The isothermal compression dynamics of ternary Ti-6 Al-4 V alloy with initial martensitic structures were investigated in the high temperature range 1083-1173 K and moderate strain rate regime 0.01-10 s~(-1).Shear banding was found to still dominate the deformation mechanism of this process,despite its nonadiabatic feature.The constitutive equation was derived with the aid of Zener-Hollomon parameter,which predicted the apparent activation energy as 534.39 kJ/mol.A combination of higher deformation temperature and lower strain rate suppressed the peak flow stress and promoted the evolution of shear bands.Both experiments and calculations demonstrated that a conspicuous temperature rise up to 83 K could be induced by severe plastic deformation.This facilitated the dynamic recrystallization of deformed martensites,as evidenced by the measured microhardness profiles across shear bands.     

Development of dislocation-based unified material model for simulating microstructure evolution in multipass hot rolling

Dynamic recrystallization and recovery are two competing processes. Both may continue after hot deformation, such as during passes in multipass hot rolling processes, reducing dislocation density of materials and allowing larger plastic deformation to be achieved. The main objective of this research is to develop a set of mechanism-based unified viscoplastic constitutive equations which model the evolution of dislocation density, recrystallization and grain size during and after hot plastic deformation. This set of constitutive equations are determined for a C-Mn steel using an evolutionary programming (EP) optimization technique and implemented into the commercial finite element (FE) solver ABAQUS for process simulations. Numerical procedures to simulate multipass rolling are developed. FE analysis is carried out to simulate the evolution of grain size, dynamic/static recrystallization and recovery, and to rationalize their effects on the viscoplastic flow of the material in a two-pass hot rolling process.     

Residual stress field in steel samples during plastic deformation and recovery processes

An X-ray diffraction method was applied to measure residual stresses and stored elastic energy in deformed and annealed polycrystalline ferritic and austenitic steel samples. The orientation distribution of plastic incompatibility second-order stresses created during elastoplastic deformation was determined and presented in Euler space. Using deformation models, these stresses were correlated with different types of intergranular interactions occurring in the studied materials. An important decrease of the first- and the second-order residual stresses was observed during recovery and recrystallisation processes. Diffraction peak widths, related to dislocation density, were studied and correlated with stress variation during annealing process. Differences in stress relaxation between ferritic and austenitic samples were explained by different values of the stacking fault energy, which influences dislocation climb and cross-slip.     

Recrystallization of the superplastic Zn-Cd alloy

The initial fine-grained structure of superplastic alloys is attained by rather complicated thermomechanical treatment. Various processes can take place during annealing of such materials due to the release of deformation energy stored in specimens in the course of their preparation. This paper deals with the measurements of the electrical resistivity annealing curves of the superplastic Zn-0·25 wt. % Cd alloy. The pronounced electrical resistivity drop due to the recrystallization was observed between 320 and 360 K. The value of kinetic exponentn from Avrami's equation was found to be 1<n<2;n increases with increasing temperature of annealing. The influence of superplastic deformation on the electrical resistivity annealing curves was investigated. Due to the softening during superplastic deformation the recrystallization becomes less pronounced and shifts to higher temperatures. The kinetic exponentn decreases to unity and its temperature dependence disappears. The attempt was made to explain these results on the ground of literature findings on the kinetics of recrystallization after preceding hot deformation.     

An evolutional approach to the numerical construction of polycrystalline structures using the Voronoi tessellation

We developed a method to numerically construct a polycrystalline structure with a specified grain-size distribution by utilizing the genetic algorithm, which is suitable for global optimization of complex parametric space, and the Voronoi tessellation. The computational merits of our new method compared with the conventional optimization method are: (a) it is able to escape more easily from local minima of the penalty function; and (b) it is suitable for parallel processing. The method efficiently creates input data for numerical simulations of microstructural evolution such as recrystallization, grain growth, deformation, and fracture.     

Surface nanocrystallization mechanism of a rare earth magnesium alloy induced by HVOF supersonic microparticles bombarding

A nanostructured surface layer with a thickness up to 60 μm was produced on a rare earth Mg-Gd-Y magnesium alloy using a new process named HVOF-SMB (high velocity oxygen-fuel flame supersonic microparticles bombarding). The microstructural features of the treated surface at various depth of the deformed layer were characterized by optical microscopy (OM), transmission electron microscopy (TEM) and high-resolution transmission electron microscopy (HRTEM) with an aim to reveal the formation mechanism. Results showed that three steps during grain refinement process were found, i.e., twinning dominates the plastic deformation and divides the coarse grains into finer twin platelets at the initial stage, stacking faults are generated and a number of dislocation slip systems are activated leading to the cross slips with increasing strain and strain rate, eventually high-density dislocation networks, dislocation cells and dislocation arrays are formed, which further subdivides the twin platelets and residual microbands into sub-microstructures. As a result, homogeneous nanostructure with a grain size of about 10-20 nm is formed through dynamic recrystallization in the topmost surface layer. Based on the experimental observations, a grain refinement mechanism induced by plastic deformation with higher strain rate during the HVOF-SMB treatment in the rare earth Mg-Gd-Y alloy was proposed.     

Microstructural characterization and recrystallization kinetics of cold rolled copper

The recrystallization kinetics in pure copper deformed by cold rolling is mainly investigated by differential scanning calorimetry (DSC) under non-isothermal conditions. DSC curves show exothermic peaks corresponding to the stored energy released during recrystallization process. Variation of the heating rate and application of different methods allowed us to calculate two kinetic parameters of recrystallization: (i) the activation energy of the process was calculated using three different methods and (ii) the Avrami exponent was estimated using the Matusita method. On other hand, the microstructural evolution during recrystallization and grain growth of cold rolled copper were investigated at a temperature of 450 °C by scanning electron microscopy (SEM) and electron back scattered diffraction (EBSD).     

Microstructure evolution of zircaloy-4 during Ne ion irradiation and annealing: An in situ TEM investigation

The microstructural evolution of zircaloy-4 was studied, including the amorphization and recrystallization of Zr（Fe, Cr）2 precipitates, and the density of dislocations under in situ Ne ion irradiation and post annealing. The results show that irradiation at a relatively high temperature and dose induces the formation of nanocrystals in pre-amorphized Zr（Fe, Cr）2 precipitates. The recrystallized nanocrystals also have the structure of hcp-Zr（Fe, Cr）2. The formation of the nanocrystals is thought to be the consequence of competition between atomistic disordering and the recrystallization of precipitates under ion irradiation. The free energy of the nanocrystal is lower than that of the amorphous state, which is another reason for the recrystallization of the precipitates. With increased annealing temperature, the density of the nanocrystals is increased. The dislocation density sharply decreases with the increase in the annealing temperature, and its size increases.     

The relation between the polygonization and recrystallization of nickel

The recrystallization temperatures of 99·85% and 99·997% pure nickel were determined. It was found that up to a deformation of about 20% the temperature drop at the beginning and end of recrystallization was very steep and then much more gradual and linear. The greater the deformation and purity of the metal, the smaller was the degree of polygonization of the deformed material in which nucleation took place. It was also found that the growth rate of grains is slowed down at isothermic annealing. These phenomena are considered in relation to the non-uniformity of deformation and the decrease in the density of vacancies, produced in the process of deformation.     

相场再结晶储能释放模型与显微组织演变的模拟研究

在相场再结晶模型中提出了形式为 f(η_i^k,η_j^k) = E_s (η_i^de )² (1-(η_j^re)²)的冷变形储能项,并应用该模型模拟了 AZ31镁合金的再结晶过程,模拟结果和实验观测结果符合很好.研究表明, 引入储能释放模型可以实现再结晶形核物理过程的模拟; 模拟结果可以把合金在冷变形后退火的过程按照机理分为再结晶和热晶粒 长大两个阶段,模拟得出的理论再结晶时间是实验再结晶时间的2/3. 考察了冷变形应变大小对形变金属的亚晶粒尺寸和储能的影响机理和 试验结果,并将考察结果代入到改进后的再结晶模拟模型, 成功地再现了一个经典实验结果:随预先应变量的增加, 存在临界应变量对应的一个再结晶晶粒尺寸峰值.同时还给出了 这一经典实验结果的理论解释.     

Dislocation-mediated creep process in nanocrystalline Cu

Nanocrystalline Cu with average grain sizes ranging from ～ 24.4 to 131.3 nm were prepared by the electric brushplating technique.Nanoindentation tests were performed within a wide strain rate range,and the creep process of nanocrystalline Cu during the holding period and its relationship to dislocation and twin structures were examined.It was demonstrated that creep strain and creep strain rate are considerably significant for smaller grain sizes and higher loading strain rates,and are far higher than those predicted by the models of Cobble creep and grain boundary sliding.The analysis based on the calculations and experiments reveals that the significant creep deformation arises from the rapid absorption of high density dislocations stored in the loading regime.Our experiments imply that stored dislocations during loading are highly unstable and dislocation activity can proceed and lead to significant post-loading plasticity.     

Quantitative evaluation of dislocation density using minor-loop scaling relations

Relationships between minor hysteresis loops and dislocation density have been investigated at various temperatures from 10 to 600 K in polycrystalline nickel with tensile deformation. It was revealed that coefficients obtained from scaling relations between parameters of minor-loops are in linear proportion to stress at all measuring temperatures below its Curie temperature. Considering that dislocation density is generally in proportion to the square root of true stress, it is concluded that the coefficients are related with the square root of dislocation density. This method using minor hysteresis loops is useful for quantitative evaluation of dislocation density because of its very low measurement field.