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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

MDeFF2012

新蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 21.1
散金: 85
紅花: 2
帖子: 238
在線: 221.9小時
蟲號: 2146432
注冊: 2012-11-25
性別: GG
專業(yè): 極端條件下使用的金屬材料

[求助] 碰撞模型弛豫出現(xiàn)變形

建立一個球形彈丸撞擊靶板的模型，弛豫時發(fā)生變形，請大家看看，有什么不妥？
log          log.equil

units          metal
dimension       3
boundary       s p p

atom_style       atomic
neighbor       2.0 bin
neigh_modify    every 5 delay 0 check yes

lattice       fcc 4.05
region       box block 0 20 -20 20 -0.5 0.5
create_box       2 box
create_atoms 1 box

mass          1 27
mass          2 27

pair_style    eam/fs
pair_coeff    * * Al_mm.eam.fs Al Al

timestep       0.001
thermo       100
thermo_style custom step etotal pe ke temp vol press

min_style    cg
minimize       0 1.0e-5 5000 20000

reset_timestep  0

velocity       all create 300.0 47289 dist gaussian units box
fix          1 all nvt temp 300.0 300.0 0.1

compute       peatom all pe/atom
compute       keatom all ke/atom
compute       c all cna/atom 3.4571
dump          1 all custom 100 *.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom
dump_modify    1 thresh c_c > 1

run          20000
unfix          1

region       1 sphere 2.0 0.0 0.0 2.0 side in
group          1 region 1
set             group 1 type 1

region       2 block 5 INF INF INF INF INF
group          2 region 2
set             group 2 type 2

group          model union 1 2
group          void subtract all model
delete_atoms group void

velocity       model create 300.0 47289 dist gaussian units box
fix          2 model nvt temp 300.0 300.0 0.1

dump          2 model custom 100 equil.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom
dump          3 1 custom 100 equil1.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom
dump          4 2 custom 100 equil2.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom
dump          5 model custom 100 *.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom
dump_modify    5 thresh c_c > 1

run          10000

unfix          2
undump       2
undump       3
undump       4
undump       5

write_restart restart.equil

回復此樓

» 猜你喜歡

085600材料與化工調(diào)劑已經(jīng)有10人回復
279分求調(diào)劑一志愿211 已經(jīng)有18人回復
0703化學求調(diào)劑已經(jīng)有5人回復
材料292調(diào)劑已經(jīng)有3人回復
一志愿陜師大生物學071000，298分，求調(diào)劑已經(jīng)有3人回復
306求0703調(diào)劑一志愿華中師范已經(jīng)有7人回復
333求調(diào)劑已經(jīng)有3人回復
335求調(diào)劑已經(jīng)有4人回復
336求調(diào)劑已經(jīng)有4人回復
303求調(diào)劑已經(jīng)有4人回復

1樓 2014-11-15 16:48:25

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

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

弛豫鐵電體的幾種模型已經(jīng)有10人回復

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 MDeFF2012 的主題更新