版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

11B329

金蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 1714
散金: 119
紅花: 3
帖子: 554
在線: 109.8小時
蟲號: 2768861
注冊: 2013-10-31
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

[求助] 關于lammps里examples的問題（shear）已有1人參與

各位大神，希望能指點一下。
在equilibrate之前，lower與upper的溫度并沒有設定，只是設定了mobile的溫度，請問那其余兩組原子的初速度怎樣？還有就是原子的質量好像沒有給啊，不用給定嗎？
以下是shear的程序：
# 3d metal shear simulation

units metal

boundary s s p

atom_style atomic

lattice fcc 3.52

region box block 0 16.0 0 10.0 0 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
&
origin 0.5 0 0

create_atoms 1 box

pair_style eam

pair_coeff * * Ni_u3.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

set group upper type 3

# void

#region void cylinder z 8 5 2.5 INF INF

#delete_atoms region void

# temp controllers

compute new3d mobile temp

compute new2d mobile temp/partial 0 1 1

# equilibrate

velocity mobile create 300.0 5812775 temp new3d

fix 1 all 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

fix_modify 3 temp new3d

thermo 25

thermo_modify temp new3d

timestep 0.001

run 100

# shear

velocity upper set 1.0 0 0

velocity mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes

unfix 3

fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0

fix_modify 3 temp new2d

#dump 1 all atom 100 dump.shear

#dump 2 all image 100 image.*.jpg type type &

# axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0

#dump_modify 2 pad 4

#dump 3 all movie 100 movie.mpg type type &

# axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0

#dump_modify 3 pad 4

thermo 100

thermo_modify temp new2d

reset_timestep 0

run 3000

回復此樓

» 猜你喜歡

求調劑已經(jīng)有7人回復
275求調劑已經(jīng)有6人回復
招08考數(shù)學已經(jīng)有7人回復
303求調劑已經(jīng)有6人回復
306求調劑已經(jīng)有4人回復
求調劑已經(jīng)有7人回復
考研調劑已經(jīng)有4人回復
285求調劑已經(jīng)有6人回復
一志愿北京化工大學070300 學碩336求調劑已經(jīng)有7人回復
318求調劑已經(jīng)有6人回復

往前走吧！

1樓 2015-12-10 12:41:58

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

mgqqlwq

銅蟲 (著名寫手)

應助: 12 (小學生)
金幣: 12831.1
散金: 108
紅花: 14
帖子: 1282
在線: 207.6小時
蟲號: 2436148
注冊: 2013-04-25
專業(yè): 凝聚態(tài)物性I:結構、力學和

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
11B329: 金幣+10, ★★★很有幫助, 多謝！我一會試試！ 2015-12-11 07:42:15

我記得eam勢函數(shù)文件里有原子質量的信息，你可以看下pair_style eam的doc page。未設定初始溫度的group的溫度我不知道，不過是否可以試著把那兩個group的原子速度輸出出來看一下。希望對你有用。

贊一下(1人)

回復此樓

2樓2015-12-10 18:59:12

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 11B329 的主題更新