介原子分子動力學(xué)方法的發(fā)展與應(yīng)用
發(fā)布時間：2018-08-18 12:14
【摘要】：合金廣泛應(yīng)用于工業(yè)生產(chǎn),由此也得到學(xué)術(shù)界的廣泛關(guān)注與研究。但在經(jīng)典的嵌入式原子勢體系下,科學(xué)家們很難優(yōu)化得到可靠的合金勢函數(shù),這使得合金分子動力學(xué)研究嚴重滯后于實驗與理論的研究進展。為能夠利用分子動力學(xué)從原子層面研究合金的變形機制,我們基于嵌入式原子勢發(fā)展出用于合金分子動力學(xué)模擬的介原子分子動力學(xué)方法。所取得的研究進展如下：本文基于兩個基本假設(shè)“相似性假設(shè)”和“平均化假設(shè)”提出了介原子分子動力學(xué)方法。在這一方法中,所有的原子都不再代表合金中的不同元素的原子而是一個虛擬的原子—“介原子”。其勢函數(shù)的擬合也不再考慮不同元素之間的相互作用,轉(zhuǎn)而直接使用合金的實測材料屬性來優(yōu)化擬合。基于這一方法,本文自主開發(fā)了一套計算機程序,通過加權(quán)最小二乘法優(yōu)化模擬值與目標值的殘差,最后輸出滿足嵌入式原子勢格式的勢函數(shù)文件。為了降低單次優(yōu)化的時間并提升計算效率,本文提出了“幾何因子”和“白噪聲”等加速算法,這有效提升了程序優(yōu)化得到勢函數(shù)的求解能力。為了驗證本文所提出方法的可靠性,本文利用介原子分子動力學(xué)方法制作得到了一系列的銅鋁合金勢函數(shù)和虛擬可發(fā)生脆斷的面心立方勢函數(shù),并利用這些勢函數(shù)模擬研究層錯能對合金變形機制的影響和金屬的韌脆轉(zhuǎn)變效應(yīng)。本文的研究結(jié)果發(fā)現(xiàn)層錯能對于合金的變形機制有著重要影響：隨著層錯能的降低,金屬的變形機制中孿晶所占比例迅速提升。同時,自由表面能和不穩(wěn)定層錯能的比值直接決定了金屬的韌性或脆性斷裂屬性。這組勢函數(shù)的成功應(yīng)用證明了本方法的可靠性與普適性。為研究固溶體合金中的固溶強化效應(yīng),本文又對介原子分子動力學(xué)方法做了改進,提出了放縮因子介原子分子動力學(xué)方法。放縮因子的加入可以有效地在金屬中引入彌散分布的晶格畸變。并且,隨著晶格畸變的加大,當(dāng)超過晶格自身的穩(wěn)定性時,晶格就會垮塌成為非晶結(jié)構(gòu)。放縮因子的加入并不會對材料屬性有很大的影響,但又會使位錯滑動的臨界應(yīng)力Peierls應(yīng)力大幅度上升,明顯提升位錯滑動所面對的晶格阻力,進而影響材料變形機制。為了更好地推廣本文提出的介原子勢函數(shù),本文選擇了學(xué)術(shù)界近期的研究熱點“孿晶誘導(dǎo)塑形變形鋼”(TWIP鋼),并利用介原子分子動力學(xué)方法發(fā)展得到層錯能為19 mJ/m2的TWIP鋼介原子勢函數(shù),隨后利用這一勢函數(shù)系統(tǒng)研究低層錯能金屬形變孿晶的形核、生長以及位錯與孿晶的相互作用。研究發(fā)現(xiàn),形變孿晶與擴展位錯均在TWIP鋼的塑形變形過程中發(fā)揮著重要作用。而形變孿晶的孿晶界可以有效地細化晶粒。在阻礙位錯運動的同時,共格的孿晶界又提升了晶粒的容納位錯能力。在經(jīng)歷一定量塑性變形后,晶粒內(nèi)會出現(xiàn)大量點空穴或空穴管道。這些空穴主要來自于位錯割階的非保守運動或特定伯格斯矢量層錯的交叉。相對于傳統(tǒng)的嵌入式原子勢方法,本文所提出的介原子分子動力學(xué)方法可以有效降低合金勢函數(shù)的發(fā)展難度。我們相信,這一方法的提出和應(yīng)用必將對合金的分子動力學(xué)研究產(chǎn)生重要貢獻。
[Abstract]:However, under the classical embedded atomic potential system, it is difficult for scientists to optimize the reliable alloy potential function, which makes the research of alloy molecular dynamics seriously lag behind the experimental and theoretical progress. Based on the embedded atomic potential, we developed a mesoatomic molecular dynamics method for alloy molecular dynamics simulation. The research progress is as follows: Based on two basic hypotheses, similarity hypothesis and averaging hypothesis, a mesoatomic molecular dynamics method is proposed. In one method, all the atoms no longer represent the atoms of different elements in the alloy, but are a virtual atom called "mesoatom". The fitting of the potential function does not consider the interaction between different elements, but directly uses the measured material properties of the alloy to optimize the fitting. In order to reduce the time of single optimization and improve the computational efficiency, some acceleration algorithms such as "geometric factor" and "white noise" are proposed, which effectively improve the program optimization. In order to verify the reliability of the proposed method, a series of potential functions for Cu-Al alloys and virtual face-centered cubic potential functions for brittle fracture have been fabricated by using the mesoatomic molecular dynamics method. These potential functions are used to simulate the effect of stacking fault energy on the deformation mechanism of alloys and metals. The results show that stacking fault energy has an important effect on the deformation mechanism of alloys: with the decrease of stacking fault energy, the proportion of twins in the deformation mechanism of metals increases rapidly. Meanwhile, the ratio of free surface energy to unstable stacking fault energy directly determines the ductile or brittle fracture properties of metals. In order to study the solid solution strengthening effect in solid solution alloys, the mesoatomic molecular dynamics method has been improved and a scaling factor mesoatomic molecular dynamics method has been proposed. Moreover, with the increase of lattice distortion, the lattice collapses into amorphous structure when the lattice is more stable than the lattice itself. The addition of the scaling factor does not have a great influence on the material properties, but also makes the critical stress Peierls stress of the dislocation slip increase greatly, which obviously increases the lattice resistance faced by the dislocation slip, and then shadow. In order to generalize the mesoatomic potential function proposed in this paper, the twin-induced plastic deformation steel (TWIP steel) has been selected as a research hotspot in recent years, and the mesoatomic potential function of TWIP steel with stacking fault energy of 19 mJ/m2 has been developed by using the mesoatomic molecular dynamics method. The nucleation and growth of deformation twins in low stacking fault energy metals and the interaction between dislocations and twins are studied.It is found that deformation twins and extended dislocations play an important role in the deformation process of TWIP steel.The twin boundaries of deformation twins can effectively refine grains.The coherent twin boundaries are raised while the dislocation movement is hindered. After a certain amount of plastic deformation, there will be a large number of point holes or holes in the grains. These holes mainly come from the non-conservative motion of the dislocation cut order or the crossover of specific Burgers vector stacking faults. Compared with the traditional embedded atomic potential method, the mesoatomic molecular dynamics proposed in this paper is more effective. It is believed that the proposed method and its application will contribute greatly to the study of molecular dynamics of alloys.
【學(xué)位授予單位】：浙江大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2016
【分類號】：O56
本文編號：2189439
[Abstract]:However, under the classical embedded atomic potential system, it is difficult for scientists to optimize the reliable alloy potential function, which makes the research of alloy molecular dynamics seriously lag behind the experimental and theoretical progress. Based on the embedded atomic potential, we developed a mesoatomic molecular dynamics method for alloy molecular dynamics simulation. The research progress is as follows: Based on two basic hypotheses, similarity hypothesis and averaging hypothesis, a mesoatomic molecular dynamics method is proposed. In one method, all the atoms no longer represent the atoms of different elements in the alloy, but are a virtual atom called "mesoatom". The fitting of the potential function does not consider the interaction between different elements, but directly uses the measured material properties of the alloy to optimize the fitting. In order to reduce the time of single optimization and improve the computational efficiency, some acceleration algorithms such as "geometric factor" and "white noise" are proposed, which effectively improve the program optimization. In order to verify the reliability of the proposed method, a series of potential functions for Cu-Al alloys and virtual face-centered cubic potential functions for brittle fracture have been fabricated by using the mesoatomic molecular dynamics method. These potential functions are used to simulate the effect of stacking fault energy on the deformation mechanism of alloys and metals. The results show that stacking fault energy has an important effect on the deformation mechanism of alloys: with the decrease of stacking fault energy, the proportion of twins in the deformation mechanism of metals increases rapidly. Meanwhile, the ratio of free surface energy to unstable stacking fault energy directly determines the ductile or brittle fracture properties of metals. In order to study the solid solution strengthening effect in solid solution alloys, the mesoatomic molecular dynamics method has been improved and a scaling factor mesoatomic molecular dynamics method has been proposed. Moreover, with the increase of lattice distortion, the lattice collapses into amorphous structure when the lattice is more stable than the lattice itself. The addition of the scaling factor does not have a great influence on the material properties, but also makes the critical stress Peierls stress of the dislocation slip increase greatly, which obviously increases the lattice resistance faced by the dislocation slip, and then shadow. In order to generalize the mesoatomic potential function proposed in this paper, the twin-induced plastic deformation steel (TWIP steel) has been selected as a research hotspot in recent years, and the mesoatomic potential function of TWIP steel with stacking fault energy of 19 mJ/m2 has been developed by using the mesoatomic molecular dynamics method. The nucleation and growth of deformation twins in low stacking fault energy metals and the interaction between dislocations and twins are studied.It is found that deformation twins and extended dislocations play an important role in the deformation process of TWIP steel.The twin boundaries of deformation twins can effectively refine grains.The coherent twin boundaries are raised while the dislocation movement is hindered. After a certain amount of plastic deformation, there will be a large number of point holes or holes in the grains. These holes mainly come from the non-conservative motion of the dislocation cut order or the crossover of specific Burgers vector stacking faults. Compared with the traditional embedded atomic potential method, the mesoatomic molecular dynamics proposed in this paper is more effective. It is believed that the proposed method and its application will contribute greatly to the study of molecular dynamics of alloys.
【學(xué)位授予單位】：浙江大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2016
【分類號】：O56
【相似文獻】
相關(guān)期刊論文 前3條
1 ;光學(xué)概論[J];中國光學(xué)與應(yīng)用光學(xué)文摘;1999年02期
2 張剛;姜瑩;;另類視角話院士(二)[J];科學(xué)中國人;2007年08期
3 ;[J];;年期
相關(guān)博士學(xué)位論文 前1條
1 王鵬;介原子分子動力學(xué)方法的發(fā)展與應(yīng)用[D];浙江大學(xué);2016年
，本文編號：2189439
本文鏈接：http://sikaile.net/shoufeilunwen/jckxbs/2189439.html
最近更新
教材專著
