本期給大家?guī)?lái)lammps模擬帶缺陷鎳板剪切變形的in文件及其詳細(xì)注釋。初始模型如圖一所示：

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

in文件及注釋如下：

#利用eam勢(shì)函數(shù)模擬帶缺陷鎳板的剪切

#模型構(gòu)成——上下鎳板夾可動(dòng)鎳塊，鎳塊中有圓柱形缺陷，移動(dòng)上鎳板使其發(fā)生剪切變形

# 3d metal shear simulation

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

#單位制為metal

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

#剪切和厚度方向?yàn)榉侵芷谛愿?dòng)邊界，寬度方向?yàn)橹芷谛赃吔鐥l件

atom_style?????? atomic

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

#面心立方晶格，晶格常數(shù)3.52

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

#盒子總體大小，長(zhǎng)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

?#從原點(diǎn)開(kāi)始，沿上述的晶體取向生長(zhǎng)

create_atoms 1 box

#在盒子內(nèi)填滿(mǎn)原子

pair_style ?????? eam

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

#調(diào)用eam勢(shì)函數(shù)

neighbor? 0.3 bin

neigh_modify?? delay 5

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

?#定義底板區(qū)域

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

#定義頂板區(qū)域

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

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

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

#頂板+頂板設(shè)為邊界群

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

#除去邊界的地方為可動(dòng)群

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

#設(shè)置底板原子類(lèi)型為2

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

?#設(shè)置頂板原子類(lèi)型為3

# void

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

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

delete_atoms? region void

#刪除圓柱區(qū)域的原子

# temp controllers

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

#定義溫度的計(jì)算new3d(mobile區(qū)域統(tǒng)計(jì)平均)

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

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

# equilibrate

velocity??? mobile create 300.0 5812775 temp new3d

#設(shè)置原子初始速度

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

#直接調(diào)溫法，溫度控制為300K+-10K之內(nèi)

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

?#fix 3 調(diào)溫時(shí)采用new3d的計(jì)算方案

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

?#每25步輸出一次熱力學(xué)統(tǒng)計(jì)量

thermo_modify??????? temp new3d

?#每次輸出的溫度按照new3d方案來(lái)

timestep? 0.001

#時(shí)間步長(zhǎng)0.001ps

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

#弛豫100步

# shear 剪切

velocity??? upper set 1.0 0 0

#使頂板原子以1埃/ps的速度沿x正向運(yùn)動(dòng)

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

#mobile原子的初始速度從0到1線(xiàn)性變化

unfix????????? 3

#結(jié)束弛豫時(shí)的溫度約束fix 3

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

#開(kāi)始新的控溫，直接調(diào)溫法，溫度控制為300K+-10K之內(nèi)

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

#fix 3 調(diào)溫時(shí)采用new2d的計(jì)算方案

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

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

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

thermo_modify??????? temp new2d

reset_timestep??????? 0 #重新設(shè)當(dāng)前為第0步

run??????????? 3000 #總共運(yùn)行3000步

?最后，發(fā)生剪切后模型如下所示：

? ? ? ? ? ? ? ? ?圖2 發(fā)生剪切后的模型

最后，如果您有分子動(dòng)力學(xué)相關(guān)需求，可以通過(guò)微信公眾號(hào)聯(lián)系我們呦~

微信公眾號(hào)：320科技工作室

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

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

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