基于分子動力學模擬的金屬構件的彈-塑性分解方法
發(fā)布時間:2021-02-22 16:52
針對金屬多晶材料構件的分子動力學(MD)模擬,本文提出了一種新的彈-塑性分解方法.文章將MD運動軌跡分解為結構變形和熱振動,給出了計算結構變形的方法和近似公式;針對金屬多晶材料構件的當前構型,給出了基于FCC|BCC晶胞和四原子占位的四面體單元相組合的連續(xù)變形函數及計算變形梯度的算法;利用原子-連續(xù)關聯模型,發(fā)展了計算當前構型應力場和彈性張量的算法.分析了當構件承受過大載荷后在材料內部所產生的微觀缺陷,并將其分類標定為位錯、層錯、攣晶界、晶界和空位等;對于層錯和攣晶界討論了在彈性卸載過程中應滿足的剛體運動約束方程;利用極小勢能原理構造了基于當前構型的彈性卸載算法,進而給出了完整的基于MD模擬的計算彈-塑性應變的算法.最后,針對單晶銅納米線拉伸的MD模擬,計算了彈-塑性應變場,驗證了本文方法的合理性.本文提出的基于MD模擬的彈-塑性分解方法,為從微觀到宏觀的多尺度和多模型耦合計算提供了算法支撐.
【文章來源】:計算數學. 2020,42(03)北大核心
【文章頁數】:19 頁
【部分圖文】:
圖3彈-塑性分解(黑線和圓圈衷錄微缺陷)??
分,它們是在加載之前已經存在的原始缺陷,在加載過程中位于境界上的原子??一直隨著晶界滑移而變化其位置,故在彈性卸載中晶界上的原子仍應視為無約束的變形??狀態(tài).???空位是指至少具有一個晶格尺度且無原子占位的微觀缺陷,其勢能相對微弱.由于空位??周圍的原子各有其配位,在繼續(xù)加載和卸載過程中空位周邊原子均應視為無約束變形狀??態(tài).??為了使層錯和攣晶界等微觀缺陷子域在鞏4?的彈性卸載過程中保持形態(tài)不變,下面??給出它們在彈性卸載過程中應滿足的剛體運動約束方程:設缺陷4在卸載前后的位置如圖4??〔a_)所示,分別用虛線和實線標出,T為缺陷上的參考點,可取為W的形心,P為缺陷上任意??一點.(5x.lX2X3:?為全局坐標系,基向景為(ei肩。澹常;?;Txix2X3為平動坐標系,基向量與全局坐??標系一致;TxW3為固連坐標系,基向景為卸載前,固連坐標系IW2X’3與??平動坐標系rXlx2x3重合.卸載后,缺陷4移動到了新的位置,平動坐標系TxlX2x3僅隨W??發(fā)生平行移動,而固連坐標系IVlX’2x’3則不僅發(fā)生平移還發(fā)生了旋轉,見圖4(a)所示.??(4??(a)彈性卸載中缺陷的運動??0?)固連坐標系與平動坐標系??圖4缺陷的剛體運動:??在卸載前后,平動坐標系和固連坐標系的基向貴之間的變換關系為??e’i??ei??COS?^11?COS?汐?12?COS?汐?13??e’2??=M??e2??,M?=??COS?汐21?COS?〇22?COS?汐23??.e;3?_??_?e3?_??COS?沒3i?COS?0^2?COS?沒33??其中M為變換矩陣,%?(U'?=?表末向量?<?和e;?的夾角,如圖4?
3期?崔俊芝等:基于分子動力學模擬的金屬構件的彈-塑性分解方法?293??W初始構型?(b)中橫截面(010)??閣5單晶CJu納米線拉伸的數值構型??(a)加載過程與Y方向應變?(b)第40加載步后的橫截面?(e)第4(3加載步后的<r13??圖6第49加載步完成后的橫截面變形與應力ff13??(a)第49步加載后的縱截面?(b)第肋少加載后縱截面上關法點??圖7第49加載完成后的數值構型??表1初始構型的彈性常數??彈性??模量??Du??D22??D33??Di2(D2l)??Di3(D3l)??乃23(乃32)??-〇44??D55??Dqq??Wang??167.8??167.8??167.8??113.5??113.5??113.5??74.5??74.5??74.5??Kittel??168.4??168.4??168.4??121.4??121.4??121.4??75.4??75.4??75.4??/3〇-2??169.8??169.8??169.8??122.6??122.6??122.6??76.2??76.2??76.2??表2當前構型執(zhí)上4個點域的彈性常數(單位:GPa)??彈性??模量??Du??D22??D33??Di2(D2l)??Di3(D3l)??乃23(乃32)??-〇44??D55??Dqq??.K域2??152.5??169.6??165.2??112.7??143.1??106.4??62.5??56.2??92.6??K域3??152.4??164.3??164.2??114.4??141.5??100.4??63.6??49.4??90.4??K域6??152.
【參考文獻】:
期刊論文
[1]Macroscopic damping model for structural dynamics with random polycrystalline configurations[J]. Yantao Yang,Junzhi Cui,Yifan Yu,Meizhen Xiang. Acta Mechanica Sinica. 2018(03)
[2]基于原子模型的非局域連續(xù)模型[J]. 向美珍,崔俊芝,田霞. 中國科學:物理學 力學 天文學. 2011(03)
本文編號:3046286
【文章來源】:計算數學. 2020,42(03)北大核心
【文章頁數】:19 頁
【部分圖文】:
圖3彈-塑性分解(黑線和圓圈衷錄微缺陷)??
分,它們是在加載之前已經存在的原始缺陷,在加載過程中位于境界上的原子??一直隨著晶界滑移而變化其位置,故在彈性卸載中晶界上的原子仍應視為無約束的變形??狀態(tài).???空位是指至少具有一個晶格尺度且無原子占位的微觀缺陷,其勢能相對微弱.由于空位??周圍的原子各有其配位,在繼續(xù)加載和卸載過程中空位周邊原子均應視為無約束變形狀??態(tài).??為了使層錯和攣晶界等微觀缺陷子域在鞏4?的彈性卸載過程中保持形態(tài)不變,下面??給出它們在彈性卸載過程中應滿足的剛體運動約束方程:設缺陷4在卸載前后的位置如圖4??〔a_)所示,分別用虛線和實線標出,T為缺陷上的參考點,可取為W的形心,P為缺陷上任意??一點.(5x.lX2X3:?為全局坐標系,基向景為(ei肩。澹常;?;Txix2X3為平動坐標系,基向量與全局坐??標系一致;TxW3為固連坐標系,基向景為卸載前,固連坐標系IW2X’3與??平動坐標系rXlx2x3重合.卸載后,缺陷4移動到了新的位置,平動坐標系TxlX2x3僅隨W??發(fā)生平行移動,而固連坐標系IVlX’2x’3則不僅發(fā)生平移還發(fā)生了旋轉,見圖4(a)所示.??(4??(a)彈性卸載中缺陷的運動??0?)固連坐標系與平動坐標系??圖4缺陷的剛體運動:??在卸載前后,平動坐標系和固連坐標系的基向貴之間的變換關系為??e’i??ei??COS?^11?COS?汐?12?COS?汐?13??e’2??=M??e2??,M?=??COS?汐21?COS?〇22?COS?汐23??.e;3?_??_?e3?_??COS?沒3i?COS?0^2?COS?沒33??其中M為變換矩陣,%?(U'?=?表末向量?<?和e;?的夾角,如圖4?
3期?崔俊芝等:基于分子動力學模擬的金屬構件的彈-塑性分解方法?293??W初始構型?(b)中橫截面(010)??閣5單晶CJu納米線拉伸的數值構型??(a)加載過程與Y方向應變?(b)第40加載步后的橫截面?(e)第4(3加載步后的<r13??圖6第49加載步完成后的橫截面變形與應力ff13??(a)第49步加載后的縱截面?(b)第肋少加載后縱截面上關法點??圖7第49加載完成后的數值構型??表1初始構型的彈性常數??彈性??模量??Du??D22??D33??Di2(D2l)??Di3(D3l)??乃23(乃32)??-〇44??D55??Dqq??Wang??167.8??167.8??167.8??113.5??113.5??113.5??74.5??74.5??74.5??Kittel??168.4??168.4??168.4??121.4??121.4??121.4??75.4??75.4??75.4??/3〇-2??169.8??169.8??169.8??122.6??122.6??122.6??76.2??76.2??76.2??表2當前構型執(zhí)上4個點域的彈性常數(單位:GPa)??彈性??模量??Du??D22??D33??Di2(D2l)??Di3(D3l)??乃23(乃32)??-〇44??D55??Dqq??.K域2??152.5??169.6??165.2??112.7??143.1??106.4??62.5??56.2??92.6??K域3??152.4??164.3??164.2??114.4??141.5??100.4??63.6??49.4??90.4??K域6??152.
【參考文獻】:
期刊論文
[1]Macroscopic damping model for structural dynamics with random polycrystalline configurations[J]. Yantao Yang,Junzhi Cui,Yifan Yu,Meizhen Xiang. Acta Mechanica Sinica. 2018(03)
[2]基于原子模型的非局域連續(xù)模型[J]. 向美珍,崔俊芝,田霞. 中國科學:物理學 力學 天文學. 2011(03)
本文編號:3046286
本文鏈接:http://sikaile.net/kejilunwen/lxlw/3046286.html