级联碰撞对层状珠光体力学行为的影响

铁素体合金钢是目前在核能工程界应用最为广泛的一种金属结构材料,以渗碳体和铁素体基体构成的层状珠光体是铁素体合金钢中常见的金相结构。深入理解辐照效应对层状珠光体力学性能的影响对高辐照条件下铁素体钢的材料设计与寿命评估有着重要的理论参考意义。基于以上考虑,本文采用分子动力学（MD）模拟,研究了连续低能铁原子级联碰撞对渗碳体/铁素体两相界面的破坏情况,探讨了经历不同程度级联碰撞的两相结构在单向拉伸以及压缩荷载下的初始屈服情况。通过对MD模拟结果的深入分析,得到了以下主要结论：a.辐照会破坏渗碳体/铁素体两相界面的失配位错结构,引起渗碳体的分解,并促进碳原子向铁素体的扩散;b.在单轴拉伸荷载作用下,级联碰撞会使初始屈服机制由{112}<111>位错滑移系的开动转变为间隙原子团簇附近位错环的形核与长大;c.在单轴压缩荷载作用下,级联碰撞会使初始塑性变形机制由{110}<111>滑移系的开动转变为{112}<111>滑移系的开动;d)无论在单轴拉伸还是压缩情况下,级联碰撞（及辐照效应）都会导致位错初始形核应力的提升。本文的研究结果为铁素体合金钢的辐照硬化和辐照脆化行为提供了新的微观解释,对于辐照条件下铁素体合金钢材料的优化设计有一定的参考意义。     

球化退火态重载车轮钢CL70磨损性能及组织演化

对CL70重载车轮钢进行球化退火获得球状珠光体组织. 将组织为片状珠光体和球状珠光体的CL70分别与U75V钢轨钢在滚动摩擦磨损试验机上以相同参数进行干摩擦纯滚动试验. 使用称重法测量磨损量、利用带电子背散射衍射附件的扫描电子显微镜及显微硬度计对两种组织形态的试样运行表面进行组织及硬度变化情况的观察与分析. 结果表明:球状珠光体组织磨损性能不及片状珠光体组织. 二者的磨损机制和强化机制不同，片状珠光体组织以疲劳磨损为主，通过塑性变形和晶粒不断细化至纤维状再到纳米晶，位错不断累积并达到良好的强化效果. 球状珠光体组织以黏着磨损为主，只有铁素体和少量渗碳体变形和碎化，硬度提升主要来自于渗碳体颗粒周围的位错集中，硬化效果较差.      

基于微结构损伤的珠光体材料本构行为研究

珠光体材料由许多随机取向的珠光体团组成，每一珠光体团又由铁素体片和渗碳体片交替 叠合而成. 实验表明具有较小片层间距的珠光体材料具有较高的强度和寿命. 通过分析构成 珠光体团的各相及界面的损伤特征，利用损伤耗散功的概念，得到了铁素体、渗碳体与界面 损伤的具有统一形式的演化规律，将其嵌入基于细片层微结构和非经典塑性理论的珠光体团 的弹塑性本构方程，利用Hill自洽方法，得到了珠光体材料的损伤本构描述. 所得到的损 伤本构描述显含表征微结构特征的片层间距，容易证明具有较小片层间距的珠光体材料具有 较好的综合力学性能. 对珠光体材料的拉伸和循环拉压特性进行了分析，得到了与实验一 致的结果.     

密排六方金属中的孪生及孪晶位错机制

密排六方晶体结构金属中可同时启动的滑移系少,孪生成为密排六方金属中重要的塑性变形形式．由于密排六方金属复杂的晶体结构,均匀切变不能保证所有晶格点都能与基体形成对称的晶体结构,因此密排六方金属的孪生通常为滑移和原子重组（shuffle）机制相结合．本文以密排六方金属中常见的{101 ̅2}、{101 ̅1}、{112 ̅2}及{102 ̅1}孪生为例,阐述不同类型孪生过程中的孪晶位错机制．分析表明,由于原子重组机制的参与,密排六方金属的孪生可以通过不同形式的孪晶位错实现．以上四种密排六方金属孪晶中,只有{112 ̅2}孪生中的一层孪晶位错是纯剪切机制,其余的孪生机制都需要原子重组的参与．孪生机制可以大致分为滑移主导、原子重组主导以及滑移-重组相结合的机制．当孪生类型确定时,即第一不畸变面（孪晶面）k_1（和孪晶剪切方向η_1）确定时,不同孪晶位错机制对应的孪晶剪切大小和方向均不同,第二不畸变面k_2和共轭剪切方向η_2也不相同,所导致孪晶的拉压性质也不同．不同剪切方向和大小的孪晶位错机制有可能在不同应力和温度条件下被激活,从而作为密排六方金属塑性的重要来源．     

金属材料辐照缺陷演化的分子动力学研究进展

辐照条件下，高能粒子在金属材料内部引入稠密的辐照缺陷，导致材料力学性能严重退化，缩短材料服役寿命，是辐照材料研究的关键问题。辐照缺陷多处在纳米尺度，故分子动力学方法是模拟辐照缺陷的有力工具，近年来被广泛用于研究辐照缺陷演化。本文总结了金属材料中辐照缺陷演化的分子动力学研究进展，介绍了级联碰撞、点缺陷、空洞、氦泡、Frank位错环、层错四面体等辐照缺陷，及其与位错、晶界等微结构的相互作用。分子动力学方法揭示的机制与模型，深化了学界对辐照效应的认识，有助于提高辐照材料力学性能和设计耐辐照材料。     

纳米压痕中位错运动的原子模拟

利用分子动力学方法模拟研究了金刚石压头压入Ni薄膜(111)晶面的纳米压痕过程中薄膜进入初始塑性后的纳观机制,采用中心对称参数(CSP)研究不同压入深度时薄膜内部的位错的萌生和生长情况.结果表明:压痕力-压痕深度曲线的每一次的剧烈的振荡,都是一次能量释放的过程,在薄膜内部的位错生长也最剧烈.加载过程中,压入深度为0.66nm时出现位错(层错),压入深度为0.93nm时出现明显的位错形核,随着压入深度的增加,多个位错形核相互作用形成梯杆位错.压入深度为1.4nm时,梯杆位错旁出现了棱柱形不全位错环,随着压入深度的增加,棱柱形不全位错环沿着{111}滑移面运动.在最大压入深度处,薄膜塑性形变达到最大.     

CrMnFeCoNi纳米晶温度相关的拉伸行为研究

摘要：高熵合金是一种由多种主元元素组成的新型合金。实验研究表明等原子比CrMnFeCoNi高熵合金在低温下具有比室温更高的拉伸强度和断裂韧性。本文针对这一现象，利用分子动力学模拟对平均晶粒尺寸为6 nm的CrMnFeCoNi纳米晶在300、200和77 K下分别进行拉伸模拟。模拟研究揭示了纳米尺度CrMnFeCoNi高熵合金力学行为的温度效应和强韧机理。微结构演化分析表明：低温下，塑性变形阶段，滑移系开动的较少，位错滑移所受的阻力越大，屈服强度和抗拉强度越大；模型破坏时，孔洞缺陷形核较慢，更多孔洞缺陷演化成断口，更多的断口分摊拉伸应变，使得高熵合金纳米晶的低温韧性更好。     

岩石拉压蠕变试验仪研制及应用

为了探索岩石在周期性恒定拉伸、压缩荷载作用下的蠕变行为,结合杠杆原理,设计研制了一种岩石杠杆式拉伸、压缩蠕变试验仪。该仪器具有挂重质量小、可方便切换拉、压荷载等特点。首先,通过标定好的数显式拉压荷重传感器对该试验仪拉伸、压缩荷载进行了校正,得出试验仪压缩、拉伸荷载杠杆扩力比分别为81.29倍与59.46倍,挂重质量与施加在岩样荷载呈线性关系(压缩作用下线性关系相关系数R2=0.99975,拉伸作用下R2=0.9991),荷载施加稳定。最后,采用该试验仪对红砂岩进行了单轴恒定拉、压循环荷载下的蠕变试验,探讨了受荷岩样拉压蠕变规律。上述成果丰富了岩石蠕变测试与研究内容,有助于岩石力学试验测试的发展。     

金属材料辐照缺陷演化的分子动力学研究进展

辐照条件下,高能粒子在金属材料内部引入稠密的辐照缺陷,导致材料力学性能严重退化,缩短材料服役寿命,是辐照材料研究的关键问题。辐照缺陷多处在纳米尺度,故分子动力学方法是模拟辐照缺陷的有力工具,近年来被广泛用于研究辐照缺陷演化。本文总结了金属材料中辐照缺陷演化的分子动力学研究进展,介绍了级联碰撞、点缺陷、空洞、氦泡、Frank位错环、层错四面体等辐照缺陷,及其与位错、晶界等微结构的相互作用。分子动力学方法揭示的机制与模型,深化了学界对辐照效应的认识,有助于提高辐照材料力学性能和设计耐辐照材料。     

CrMnFeCoNi纳米晶温度相关的拉伸行为研究

摘要：高熵合金是一种由多种主元元素组成的新型合金。实验研究表明等原子比CrMnFeCoNi高熵合金在低温下具有比室温更高的拉伸强度和断裂韧性。本文针对这一现象,利用分子动力学模拟对平均晶粒尺寸为6 nm的CrMnFeCoNi纳米晶在300、200和77 K下分别进行拉伸模拟。模拟研究揭示了纳米尺度CrMnFeCoNi高熵合金力学行为的温度效应和强韧机理。微结构演化分析表明：低温下,塑性变形阶段,滑移系开动的较少,位错滑移所受的阻力越大,屈服强度和抗拉强度越大;模型破坏时,孔洞缺陷形核较慢,更多孔洞缺陷演化成断口,更多的断口分摊拉伸应变,使得高熵合金纳米晶的低温韧性更好。     

Influence of single crystal orientation on homogeneous dislocation nucleation under uniaxial loading

Atomistic simulations are used to investigate how the stress required for homogeneous nucleation of partial dislocations in single crystal copper under uniaxial loading changes as a function of crystallographic orientation. Molecular dynamics is employed based on an embedded-atom method potential for Cu at 10 and 300 K. Results indicate that non-Schmid parameters are important for describing the calculated dislocation nucleation behavior for single crystal orientations under tension and compression. A continuum relationship is presented that incorporates Schmid and non-Schmid terms to correlate the nucleation stress over all tensile axis orientations within the stereographic triangle. Simulations investigating the temperature dependence of homogeneous dislocation nucleation yield activation volumes of ≈0.5- and activation energies of . For uniaxial compression, full dislocation loop nucleation is observed, in contrast to uniaxial tension. One of the main differences between uniaxial tension and compression is how the applied stress is resolved normal to the slip plane on which dislocations nucleate—in tension, this normal stress is tensile, and in compression, it is compressive. Last, the tension-compression asymmetry is examined as a function of loading axis orientation. Orientations with a high resolved stress normal to the slip plane on which dislocations nucleate have a larger tension-compression asymmetry with respect to dislocation nucleation than those orientations with a low resolved normal stress. The significance of this research is that the resolved stress normal to the slip plane on which dislocations nucleate plays an important role in partial (and full) dislocation loop nucleation in FCC Cu single crystals.     

On the evolution of lattice deformation in austenitic stainless steels—The role of work hardening at finite strains

In this work, a three dimensional crystal plasticity-based finite element model is presented to examine the micromechanical behaviour of austenitic stainless steels. The model accounts for realistic polycrystal micromorphology, the kinematics of crystallographic slip, lattice rotation, slip interaction (latent hardening) and geometric distortion at finite deformation. We utilise the model to predict the microscopic lattice strain evolution of austenitic stainless steels during uniaxial tension at ambient temperature with validation through in situ neutron diffraction measurements. Overall, the predicted lattice strains are in very good agreement with those measured in both longitudinal and transverse directions (parallel and perpendicular to the tensile loading axis, respectively). The information provided by the model suggests that the observed nonlinear response in the transverse {200} grain family is associated with a competitive bimodal evolution of strain during inelastic deformation. The results associated with latent hardening effects at the microscale also indicate that in situ neutron diffraction measurements in conjunction with macroscopic uniaxial tensile data may be used to calibrate crystal plasticity models for the prediction of the inelastic material deformation response.     

Measurement of Strain Localization Resulting from Monotonic and Cyclic Loading at 650 <Superscript>∘</Superscript>C in Nickel Base Superalloys

A methodology is presented for the use of the oxide scale that develops in polycrystalline Ni-base superalloys at service temperature, as a speckle pattern for μm-scale resolution strain measurements. Quantitative assessment of the heterogeneous strain field at the grain scale is performed by high-resolution SEM digital image correlation under monotonic and cyclic loading in polycrystalline Ni-base superalloys up to 650 ^°C. In the René 88DT superalloy, strain localization is observed near twin boundaries during low cycle fatigue (LCF) at intermediate temperatures, correlating with activation of {111} 〈110〉 and {111} 〈112〉 slip systems. A strong correlation between the microstructural configuration that promotes strain localization during monotonic loading and crack initiation at 650 ^°C in low cycle fatigue was observed.     

Texture development and plastic anisotropy of B.C.C. strain hardening sheet metals

A Taylor-like polycrystal model is adopted here to investigate the plastic behavior of body centered cubic (b.c.c.) sheet metals under plane-strain compression and the subsequent in-plane biaxial stretching conditions. The <111> pencil glide system is chosen for the slip mechanism for b.c.c. sheet metals. The {110} <111> and {112} <111> slip systems are also considered. Plane-strain compression is used to simulate the cold rolling processes of a low-carbon steel sheet. Based on the polycrystal model, pole figures for the sheet metal after plane-strain compression are obtained and compared with the corresponding experimental results. Also, the simulated plane-strain stress—strain relations are compared with the corresponding experimental results. For the sheet metal subjected to the subsequent in-plane biaxial stretching and shear, plastic potential surfaces are determined at a given small amount of plastic work. With the assumption of the equivalence of the plastic potential and the yield function with normality flow, the yield surfaces based on the simulations for the sheet metal are compared with those based on several phenomenological planar anisotropic yield criteria. The effects of the slip system and the magnitude of plastic work on the shape and size of the yield surfaces are shown. The plastic anisotropy of the sheet metal is investigated in terms of the uniaxial yield stresses in different planar orientations and the corresponding values of the anisotropy parameter R, defined as the ratio of the width plastic strain rate to the through-thickness plastic strain rate under in-plane uniaxial tensile loading. The uniaxial yield stresses and the values of R at different planar orientations from the polycrystal model can be fitted well by a yield function recently proposed by Barlat et al. (1997b).     

A dislocation density based crystal plasticity finite element model: Application to a two-phase polycrystalline HCP/BCC composites

We present a multiscale model for anisotropic, elasto-plastic, rate- and temperature-sensitive deformation of polycrystalline aggregates to large plastic strains. The model accounts for a dislocation-based hardening law for multiple slip modes and links a single-crystal to a polycrystalline response using a crystal plasticity finite element based homogenization. It is capable of predicting local stress and strain fields based on evolving microstructure including the explicit evolution of dislocation density and crystallographic grain reorientation. We apply the model to simulate monotonic mechanical response of a hexagonal close-packed metal, zirconium (Zr), and a body-centered cubic metal, niobium (Nb), and study the texture evolution and deformation mechanisms in a two-phase Zr/Nb layered composite under severe plastic deformation. The model predicts well the texture in both co-deforming phases to very large plastic strains. In addition, it offers insights into the active slip systems underlying texture evolution, indicating that the observed textures develop by a combination of prismatic, pyramidal, and anomalous basal slip in Zr and primarily {110}〈111〉 slip and secondly {112}〈111〉 slip in Nb.     

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.     

Micro-plasticity of surface steps under adhesive contact: Part II—Multiple-dislocation mediated contact hardening

The study of micro-plastic behavior of rough surface contacts is the critical link towards a fundamental understanding of contact, friction, adhesion, and surface failures at small length scales. In the companion paper (Yu, H.H., Shrotriya, P., Gao, Y.F., Kim, K.-S., 2007. Micro-plasticity of surface steps under adhesive contact. Part I. Surface yielding controlled by single-dislocation nucleation. J. Mech. Phys. Solids 55, 489–516), we have studied the onset of surface yielding due to single-dislocation nucleation from a stepped surface under adhesive contact. Here we analyze the contact hardening behavior due to multiple dislocations in a two-dimensional dislocation model. Continuum micro-mechanical analyses are used to derive the configurational force on the dislocation, while a modified Rice–Thomson criterion is used to model dislocation nucleation. Dislocations nucleated from the surface step are stabilized and pile up as a result of the balance between the resolved driving force and the non-zero lattice resistance in the solid. The dislocation pileup will exert a strong back stress to prevent further dislocation nucleation and thus lead to the contact hardening behavior, the degree of which depends on the slip-plane orientation. Particularly, we find that dislocation interactions between two slip planes can make the contact loading order-of-magnitude easy to nucleate multiple dislocations, which is thus named “latent softening”. A mechanistic explanation shows that the latent softening is closely related to the stress-concentration mode mixity at the surface step. Dislocation nucleation will modify the geometric characteristics of the surface step, so that the contact-induced stress state near the step, as described by the mode mixity, changes, which influences the subsequent dislocation nucleation. Our calculations show that the dislocation pileup on one slip plane can even cause the spontaneous dislocation nucleation on the other slip plane without further increase of the contact load. Furthermore, it is found that rough surface contacts at small length scale can lead to the dislocation segregation and the formation of a surface tensile sub-layer. The discrete-dislocation model presented here and in the companion paper provides novel insights in bridging the atomistic simulations and continuum plastic flow analysis of surface asperity contact.     

完整岩石拉压应力状态下的张裂破坏准则

一点的应力状态可由3个主应力$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示, 当规定主应力以压为正时, 沿最大主应力$\sigma_{1}$方向将产生收缩变形, 若中间主应力$\sigma_{2}$和最小主应力$\sigma_{3}$都远小于$\sigma_{1}$, 则沿$\sigma_{2}$和$\sigma_{3}$方向会产生横向扩张变形, 当横向扩张变形达到一定极限时, 将会在平行于$\sigma _{1}$的方向产生张裂破坏. 如何建立这种张裂破坏的强度准则目前尚缺乏研究, 最大拉应变理论(第二强度理论)有时被用来解释张裂破坏, 但最大拉应变理论难以应用于三向受力状态. 本文分别用$\varepsilon_{1}$, $\varepsilon_{2}$表示最大张应变和次大张应变, 则最大拉应变理论认为当$\varepsilon_{1}$达到单向拉伸屈服应变时, 材料将产生破坏. 而本文将根据$\varepsilon_{1}+\varepsilon_{2}$之和达到极限值$\varepsilon_u$来建立张裂破坏准则. 可以证明$\varepsilon_{1} +\varepsilon_{2}$所表示的是$\sigma_{1}$主平面的面积增长率. 当$\sigma_{3}<\sigma_{2} \ll \sigma_{1}$时, 大部分岩石都具有脆性破坏的特点, 所以可将破坏前的岩石视为满足广义胡克定律的线弹性材料, 这样用$\varepsilon_{1}$, $\varepsilon_{2}$表示的强度准则可通过$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示. 在这个过程中还可考虑岩石在拉伸和压缩时具有不同弹性参数和强度的特点, 并可通过单向拉伸和单向压缩的破坏状态来确定$\varepsilon_u$. 不管$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$是压应力, 还是拉应力, 或者$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$中有拉有压的情形, 基于$\varepsilon_{1} +\varepsilon_{2} =\varepsilon_u$都可建立相应的强度准则. 所建立的准则可以反映中间应力$\sigma_{2}$对强度的影响规律, 通过建立的强度准则还可以证明: 静水拉力能引起屈服, 而静水压力不能产生屈服; 压缩破坏能使塑性体积增大, 其结果比Mohr-Coulomb准则更能反映实际情形. 并通过拉压应力状态下的试验数据验证了所建立的强度准则, 所得理论计算结果和已有的试验数据吻合得很好. 通过提出的强度准则和圆盘劈裂的试验结果, 可获得更为可靠的岩石单轴抗拉强度.