本期給大家帶來lammps模擬帶缺陷鎳板剪切變形的in文件及其詳細注釋。初始模型如圖一所示：

? ? ? ? ? ? ? ?? 圖1 生成的初始模型?

in文件及注釋如下：

#利用eam勢函數模擬帶缺陷鎳板的剪切

#模型構成——上下鎳板夾可動鎳塊，鎳塊中有圓柱形缺陷，移動上鎳板使其發生剪切變形

# 3d metal shear simulation

units?????????????????? metal????

#單位制為metal

Boundary ??????? s s p?????

#剪切和厚度方向為非周期性浮動邊界，寬度方向為周期性邊界條件

atom_style?????? atomic

lattice??????????????? fcc 3.52??

#面心立方晶格，晶格常數3.52

region??????????????? box block 0 16.0 0 10.0 0 2.828427????????

#盒子總體大小，長16埃米，寬10，厚2.828427

create_box?????? 3 box

lattice??????????????? fcc 3.52 orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1

?#從原點開始，沿上述的晶體取向生長

create_atoms 1 box

#在盒子內填滿原子

pair_style ?????? eam

pair_coeff???????? * * Ni_u3.eam

#調用eam勢函數

neighbor? 0.3 bin

neigh_modify?? delay 5

region??????????????? lower block INF INF INF 0.9 INF INF

?#定義底板區域

region??????????????? upper block INF INF 6.1 INF INF INF

#定義頂板區域

group???????????????? lower region lower

group???????????????? upper region upper

group???????????????? boundary union lower upper

#頂板+頂板設為邊界群

group???????????????? mobile subtract all boundary

#除去邊界的地方為可動群

set???????????? group lower type 2

#設置底板原子類型為2

set???????????? group upper type 3

?#設置頂板原子類型為3

# void

region??????????????? void cylinder z 8 3.535534 2.5 INF INF

#以(8，3.535534)為軸線，2.5埃米為半徑定義出圓柱區域

delete_atoms? region void

#刪除圓柱區域的原子

# temp controllers

compute?????????? new3d mobile temp

#定義溫度的計算new3d(mobile區域統計平均)

compute?????????? new2d mobile temp/partial 0 1 1

#定義溫度的計算new2d,忽略x方向速度

# equilibrate

velocity??? mobile create 300.0 5812775 temp new3d

#設置原子初始速度

fix????????????? 1 all nve

#nve系綜

fix????????????? 2 boundary setforce 0.0 0.0 0.0

#固定邊界原子，鋼化原子，便于加載

fix????????????? 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0

#直接調溫法，溫度控制為300K+-10K之內

fix_modify???????? 3 temp new3d

?#fix 3 調溫時采用new3d的計算方案

thermo????????????? 25

?#每25步輸出一次熱力學統計量

thermo_modify??????? temp new3d

?#每次輸出的溫度按照new3d方案來

timestep? 0.001

#時間步長0.001ps

run??????????? 100

#弛豫100步

# shear 剪切

velocity??? upper set 1.0 0 0

#使頂板原子以1埃/ps的速度沿x正向運動

velocity??? mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes

#mobile原子的初始速度從0到1線性變化

unfix????????? 3

#結束弛豫時的溫度約束fix 3

fix????????????? 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0

#開始新的控溫，直接調溫法，溫度控制為300K+-10K之內

fix_modify???????? 3 temp new2d

#fix 3 調溫時采用new2d的計算方案

dump???????????????? 1 all atom 100 dump.shear.void

#每隔100步輸出一次原子信息

thermo????????????? 100

thermo_modify??????? temp new2d

reset_timestep??????? 0 #重新設當前為第0步

run??????????? 3000 #總共運行3000步

?最后，發生剪切后模型如下所示：

? ? ? ? ? ? ? ? ?圖2 發生剪切后的模型

最后，如果您有分子動力學相關需求，可以通過微信公眾號聯系我們呦~

微信公眾號：320科技工作室

創作挑戰賽新人創作獎勵來咯，堅持創作打卡瓜分現金大獎

總結

以上是生活随笔為你收集整理的剪切文件_lammps模拟带缺陷镍板剪切变形（in文件及注释）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。