FeCoCrCuNi高熵合金裂纹及孔洞结构的力学与微观构象演化的分子动力学模拟研究

本文采用分子动力学模拟研究了FeCoCrCuNi高熵合金裂纹和孔洞结构在不同轴向拉伸速率下的力学与微观结构演化机理.结果表明:应变速率越高FeCoCrCuNi裂纹结构对应更高的过冲应变和过冲应力,其主要原因是高拉伸速率会导致高强度的BCC结构及孪晶结构的生成,而BCC结构及孪晶结构的产生进而会抑制应力的下降,通过应力-应变曲线,可知FeCoCrCuNi裂纹模型在轴向应力作用下表现为塑性形变.对于不同尺寸的孔洞FeCoCrCuNi裂纹模型的应力结构分析,可以得出:孔洞尺寸越大, FeCoCrCuNi裂纹结构对应的过冲应变和过冲应力越小,其主要原因是大尺寸的孔洞造成孔洞之间产生裂纹的,进而会影响这个材料的屈服应变和屈服强度.     

孪晶对Be材料冲击加-卸载动力学影响的数值模拟研究

在高压、高应变率加载条件下,孪晶变形对材料的塑性变形具有重要的贡献,而目前孪晶对金属材料的动态屈服强度、冲击响应等的影响还没有被充分揭示.为此,本文考虑孪晶变形和晶粒碎化,针对铍(Be)材料在高应变率加载下的动态力学响应发展了含孪晶的热弹-黏塑性晶体塑性模型.经过和实验结果的对比,发现该模型可以更准确地预测Be材料在动态加载下,尤其是高压动态加载下的屈服强度.进一步,基于该塑性模型研究了Be材料在冲击加载下的准弹性卸载行为,结果表明剪切波速随着压力和剪应变的变化而发生变化是材料产生准弹性卸载现象的主要原因.此外,研究了冲击波卸载过程中Be材料孪晶的演化过程,发现Be材料卸载过程中也伴随着孪晶的产生.     

孔洞对FeCrNiCoCu高熵合金拉伸性能影响的分子动力学研究

孔洞是FeCrNiCoCu高熵合金在制备过程中常见的缺陷，为此本文利用分子动力学模拟方法构建含孔洞的FeCrNiCoCu模型进行单轴拉伸模拟，探究了孔洞位置、孔洞半径和变形温度对其力学性能的影响.研究发现，在Z轴晶向为[111]的晶体中和晶界处的孔洞会显著降低模型的屈服应变和屈服强度，但对模型的杨氏模量影响不大.随着晶界处孔洞半径的增大，在弹性阶段，孔洞半径增大使应力集中面积增大，有利于位错形核，模型的力学性能随之降低.在塑性变形阶段，随着孔洞半径的增大，初始位错更倾向于向Z轴晶向为[001]的晶体中扩展.在中、低温条件下(T<800K)，模型保持良好的力学性能;在高温条件下，力学性能显著降低.在高温塑性变形阶段，模型中的总位错线长度较低，平均流变应力也较低.     

单晶Cu(001)薄膜塑性变形的分子动力学模拟

使用分子动力学方法,模拟研究了单晶Cu(001)薄膜在双向等轴拉伸应变下的塑性变形行为.当应变超过一定值时,样品通过产生位错、层错及孪晶而发生塑性变形.当应变相对较低时,不全位错首先在薄膜表面形核并在密排面上滑移,留下堆积层错;当应变增加时,位错在表面与内部同时成核生长,层错数量也随之增加.分析了相邻滑移面上的位错之间相互作用形成孪晶的微观过程.材料内部形成大量堆积层错及孪晶后,较大孪晶的密排面上的原子也会发生滑移,形成孪晶内部的层错结构以释放残余应力.     

NiTi形状记忆合金形变机制的应变率相关性研究

利用高速拉伸实验机在宽的应变率范围内(0.001–1200 s^-1), 研究了NiTi形状记忆合金的宏观力学性能随应变率的变化规律, 并借助透射电子显微镜深入研究了微观结构在不同应变率下的演变机制. 研究发现: NiTi合金马氏体(B19’相)孪晶的解孪晶应力随应变率的升高而近乎线性增大, 表明NiTi 合金解孪晶应力具有正向应变率相关性. 在拉伸应变率为10 s^-1的样品微观结构中发现了大量的解孪晶区域, 而当应变率进一步增大到100 s^-1和 1200 s^-1时, 在样品中没有发现解孪晶区域的存在, 样品微观组织以孪晶形式存在. 该结果表明, NiTi合金的马氏体解孪晶速率应在 10–100 s^-1范围内. 在高应变率下(≥qslant10 s^-1)均发现了热引发奥氏体相(B2)的存在, 表明随应变率的增加, 拉伸过程由等温过程逐渐变为绝热过程. 此外, 在1200 s^-1 的样品差示扫描热量曲线中还发现了一个小肩峰, 表明相变过程由一步相变变为两步相变. 关键词： NiTi形状记忆合金 高速拉伸 应变率相关性 透射电子显微镜     

基于准连续介质方法模拟纳米多晶体Ni中裂纹的扩展

通过准连续介质方法模拟了纳米多晶体Ni中裂纹的扩展过程.模拟结果显示:裂纹尖端的应力场可以导致晶界分解、层错和变形孪晶的形成等塑性形变,在距离裂纹尖端越远的位置,变形孪晶越少,在裂纹尖端附近相同距离处,层错要远多于变形孪晶.这反映了局部应力的变化以及广义平面层错能对变形孪晶的影响.计算了裂纹尖端附近区域原子级局部静水应力的分布.计算结果表明:裂纹前端晶界处容易产生细微空洞,这些空洞附近为张应力集中区,并可能促使裂纹沿着晶界扩展.模拟结果定性地反映了纳米多晶体Ni中的裂纹扩展过程,并与相关实验结果符合得很好     

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

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

镁金属孪晶变形的实验和理论模型研究进展

孪晶变形作为密排六方(HCP)镁金属的重要变形机制,对镁金属的塑性硬化、破坏和织构演变等具有重要影响。影响孪晶变形的因素较多,有取向织构、晶粒尺寸、应变率、温度、晶界和应力状态等。首先重点介绍了前3种因素对镁金属孪晶变形的影响,孪晶的启动不再单一地考虑与取向相关的Schmid定律,需结合与临近晶粒间的应变兼容,晶粒尺寸对孪晶的影响同样可以采用Hall-Petch关系描述,只是关系式的斜率比滑移更大,提高应变率对孪晶成核和成长都有一定的促进作用;然后分析了现有常见的孪晶理论模型,最后展望了孪晶变形在实验和理论模型方面的发展方向。     

原位中子衍射研究两相NiTi合金的微力学相互作用和相变机理

NiTi合金的形状记忆效应与其微观结构特征密切相关,中子衍射技术可以在力学加载过程中原位观察块体NiTi合金的相变、晶间应变以及孪晶再取向等演化特征．结合两相NiTi合金宏观应力一应变曲线呈现的四种阶段性变形特征,利用原位中子衍射技术对其变形过程中的微观结构演化进行了分析．奥氏体初始体积份额约22％,在低应变硬化阶段,晶面（110）B2和（002）B19,的应变分别突然减小和增大表明出现了应力诱发马氏体相变,奥氏体体积份额迅速减小,产生了（011）II型孪晶;同时初始马氏体也开始发生再取向,随着应变量的增加,开始出现新的{20i}型马氏体孪晶,这种孪晶引起的应变卸载时不能回复．在高应变硬化阶段孪晶变形起主导作用,衍射峰半高宽变化较小;而在应变硬化饱和阶段则以滑移机制为主,大量位错的产生使衍射峰半高宽显著增加．     

晶粒尺寸对纳米多晶铁变形机制影响的模拟研究

利用分子动力学模拟方法研究了拉伸荷载作用下晶粒尺寸对纳米多晶铁变形机制的影响.研究结果表明杨氏模量随着晶粒尺寸的减小而减小.当晶粒尺寸小于15.50 nm时,纳米多晶铁的峰值应力和晶粒尺寸之间遵循反常的Hall-Petch关系,此时晶粒旋转和晶界迁移是其塑性变形的主要变形机制;随着晶粒尺寸的增大,变形孪晶和位错滑移在其塑性变形过程中逐渐占据主导地位.裂纹的形成是导致大晶粒尺寸模型力学性能降低的主要因素.纳米多晶铁在塑性变形中会出现孪晶界的迁移和退孪晶现象.此外还研究了温度对纳米多晶铁变形机制的影响.     

Investigation on dislocation-based mechanisms of void growth and coalescence in single crystal and nanotwinned nickels by molecular dynamics simulation

Abstract

Molecular dynamics simulations were conducted to elucidate dislocation mechanisms of the void growth and coalescence in single crystal and nanotwinned nickels subjected to uniaxial tension. The simulation results reveal that twin boundary is capable of decreasing the critical stress, suppressing the emission of dislocations and reducing the overall stiffness of the crystal. A size-scale dependence of critical stress is definitely illustrated through stress–strain response, where the larger void size leads to the lower critical stress and strain. It is the successive emissions of leading partials and the subsequent trailing partials that cause the atoms on the void surfaces to escape from the void surfaces continually, and consequently the voids grow to be larger and larger with increasing strain. The voids in the nanotwinned nickel coalesce earlier than those in the single crystal nickel even though the initiation of dislocations in the former is later than that in the latter. Void fraction remains a constant during elastic deformation, while it presents a linear increase with increasing strain during plastic deformation. Evolution of void fraction during void growth and coalescence is independent on void size.     

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.     

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.     

Multiscale simulation from atomistic to continuum – coupling molecular dynamics (MD) with the material point method (MPM)

A new multiscale simulation approach is introduced that couples atomistic-scale simulations using molecular dynamics (MD) with continuum-scale simulations using the recently developed material point method (MPM). In MPM, material continuum is represented by a finite collection of material points carrying all relevant physical characteristics, such as mass, acceleration, velocity, strain and stress. The use of material points at the continuum level provides a natural connection with the atoms in the lattice at the atomistic scale. A hierarchical mesh refinement technique in MPM is presented to scale down the continuum level to the atomistic level, so that material points at the fine level in MPM are allowed to directly couple with the atoms in MD. A one-to-one correspondence of MD atoms and MPM points is used in the transition region and non-local elastic theory is used to assure compatibility between MD and MPM regions, so that seamless coupling between MD and MPM can be accomplished. A silicon single crystal under uniaxial tension is used in demonstrating the viability of the technique. A Tersoff-type, three-body potential was used in the MD simulations. The coupled MD/MPM simulations show that silicon under nanometric tension experiences, with increasing elongation in elasticity, dislocation generation and plasticity by slip, void formation and propagation, formation of amorphous structure, necking, and final rupture. Results are presented in terms of stress–strain relationships at several strain rates, as well as the rate dependence of uniaxial material properties. This new multiscale computational method has potential for use in cases where a detailed atomistic-level analysis is necessary in localized spatially separated regions whereas continuum mechanics is adequate in the rest of the material.     

Effect of stress state on deformation and fracture of nanocrystalline copper:Molecular dynamics simulation

Deformation in a microcomponent is often constrained by surrounding joined material making the component under mixed loading and multiple stress states. In this study, molecular dynamics(MD) simulation are conducted to probe the effect of stress states on the deformation and fracture of nanocrystalline Cu. Tensile strain is applied on a Cu single crystal,bicrystal and polycrystal respectively, under two different tension boundary conditions. Simulations are first conducted on the bicrystal and polycrystal models without lattice imperfection. The results reveal that, compared with the performance of simulation models under free boundary condition, the transverse stress caused by the constrained boundary condition leads to a much higher tensile stress and can severely limit the plastic deformation, which in return promotes cleavage fracture in the model. Simulations are then performed on Cu single crystal and polycrystal with an initial crack. Under constrained boundary condition, the crack tip propagates rapidly in the single crystal in a cleavage manner while the crack becomes blunting and extends along the grain boundaries in the polycrystal. Under free boundary condition, massive dislocation activities dominate the deformation mechanisms and the crack plays a little role in both single crystals and polycrystals.     

Molecular dynamic simulations of elastomer structure and its influence on anisotropic stress under time-varying strain

We investigate the evolution of polymer structure and its influence on uniaxial anisotropic stress under time-varying uniaxial strain, and the role of external control variables such as temperature, strain rate, chain length, and density, using molecular dynamics simulation. At temperatures higher than glass transition, stress anisotropy in the system is reduced even though the bond stretch is greater at higher temperatures. There is a significant increase in the stress level with increasing density. At higher densities, the uncoiling of the chains is suppressed and the major contribution to the deformation is by internal deformation of the chains. At faster rates of loading stress anisotropy increases. The deformation mechanism is mostly due to bond stretch and bond bending rather than overall shape and size. Stress levels increase with longer chain length. There is a critical value of the functionality of the cross-linkers beyond which the uniaxial stress developed increases caused primarily by bond stretching due to increased constraint on the motion of the monomers. Stacking of the chains in the system also plays a dominant role in the behaviour in terms of excluded volume interactions. Low density, high temperature, low values of functionality of cross-linkers, and short chain length facilitate chain uncoiling and chain slipping in cross-linked polymers.     

α-Fe裂纹的分子动力学研究

通过分子动力学方法，模拟了α-Fe裂纹的单轴拉伸实验中的形变过程.研究了不同晶体取向裂纹的形变特点和断裂机理，观察到各种形变现象，如位错形核和发射，位错运动，堆垛层错或孪晶的形成，纳米空洞的形成与连接等.计算结果表明，裂纹扩展是塑性过程和弹性过程相结合的过程，其中塑性过程表现为由裂尖发射的位错导致的原子切变行为，而弹性过程的发生则是由无位错区中的原子断键所导致.同时还研究了α-Fe裂纹的形变特点和断裂机理与温度场和应力场的依赖关系.     

(001)面任意方向单轴应变硅材料能带结构

首先计算了(001)晶面单轴应变张量,在此基础上采用结合形变势理论的K ·P微扰法建立了在(001)晶面内受任意方向的单轴压/张应力作用时,应变硅材料的能带结构与应力(类型、大小)及晶向的关系模型,进而分析了不同单轴应力(类型、大小)及晶向对应变硅材料导带带边、价带带边、导带分裂能、价带分裂能、禁带宽度的影响.研究结果可为单轴应变硅器件应力及晶向的选择设计提供理论依据. 关键词： 单轴应变硅 K ·P法 能带结构     

Transition of deformation mechanisms in nanotwinned single crystalline SiC

ABSTRACT

The ability to experimentally synthesise ceramic materials to incorporate nanotwinned microstructures can drastically affect the underlying deformation mechanisms and mechanics through the complex interaction between stress state, crystallographic orientation, and twin orientation. In this study, molecular dynamics simulations are used to examine the transition in deformation mechanisms and mechanical responses of nanotwinned zinc-blende SiC ceramics subjected to different stress states (uniaxial compressive, uniaxial tensile, and shear deformation) by employing various twin spacings and loading/crystallographic orientations in nanotwinned structures, as compared to their single crystal counterparts. The simulation results show that different combinations of stress states and crystal/twin orientation, and twin spacing trigger different deformation mechanisms: (i) shear localised deformation and shear-induced fracture, preceded by point defect formation and dislocation slip, in the vicinity of the twin lamellae, shear band formation, and dislocation (emission) avalanche; (ii) cleavage and fracture without dislocation plasticity, weakening the nanotwinned ceramics compared to their twin-free counterpart; (iii) severe localised deformation, generating a unique zigzag microstructure between twins without any structural phase transformations or amorphisation, and (iv) atomic disordering localised in the vicinity of coherent twin boundaries, triggering dislocation nucleation and low shearability compared to twin-free systems.     

Effect of void size and Mg contents on plastic deformation behaviors of Al-Mg alloy with pre-existing void: Molecular dynamics study

The plastic deformation properties of cylindrical pre-void aluminum-magnesium (Al-Mg) alloy under uniaxial tension are explored using molecular dynamics simulations with embedded atom method (EAM) potential. The factors of Mg content, void size, and temperature are considered. The results show that the void fraction decreases with increasing Mg in the plastic deformation, and it is almost independent of Mg content when Mg is beyond 5%. Both Mg contents and stacking faults around the void affect the void growth. These phenomena are explained by the dislocation density of the sample and stacking faults distribution around the void. The variation trends of yield stress caused by void size are in good agreement with the Lubarda model. Moreover, temperature effects are explored, the yield stress and Young's modulus obviously decrease with temperature. Our results may enrich and facilitate the understanding of the plastic mechanism of Al-Mg with defects or other alloys.