原創(chuàng)?YJ學長?LAMMPS交流站javascript:void(0);?大家好，本文是小編根據(jù)自己所學的知識，因能力有限，將介紹一種高熵合金的lammps簡單拉伸過程，本教程只適合于新手學習，希望能給讀者提供具體工程意義還讀者結合實際同時希望專業(yè)前輩提出不足。

? ? 如下圖為本文建立的合金試件模型。

01

基本思路

本文采用deform命令對高熵合金進行一個拉伸操作

02

完整代碼

完整in文件如下：

#初始模擬系統(tǒng)設置units metal #單位boundary p p p #邊界條件atom_style atomic #原子類型timestep 0.001 #時間步長#構建模型lattice fcc 3.56 #設置晶格參數(shù)region box block 0 100 0 20 0 20 #生產(chǎn)長方形區(qū)域create_box 5 boxcreate_atoms 1 box#將type為1的原子轉換為其他原子set type 1 type/ratio 2 0.2 87393 #將1原子轉換為2原子set type 1 type/ratio 3 0.5 87393 #將1原子轉換為3原子set type 1 type/ratio 4 0.5 87393 #將1原子轉換為4原子set type 3 type/ratio 5 0.5 87393 #將1原子轉換為5原子#設置原子類型mass 1 27 #Al原子mass 2 56 #Fe原子mass 3 28 #Si原子mass 4 64 #Cu原子mass 5 24 #Mg原子#拉伸設置fix 1 all deform 2000 x erate 0.05 remap x #設置拉伸方式dump 1 all atom 500 lasheng.xyz #輸出文件run 100000

03

運行結果

? ??接著，在cmd窗口運行in文件，將生成的xyz文件拖到OVITO軟件中，并在OVITO軟件中進行渲染，就產(chǎn)生了如下結果：

接著

在加上使用輸出的應力應變數(shù)據(jù)（本文中in文件未給出）在origin中繪制應力-應變曲線，如下圖所示：

更多案例盡在公眾號：LAMMPS交流站

?

總結

以上是生活随笔為你收集整理的[lammps教程] 合金拉升模拟的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。