| 13 | 1/1 | 返回列表 |
| 查看: 2571 | 回復(fù): 12 | |||
[交流]
有人做裂紋模擬的嗎?看看我這個(gè)裂紋的in文件有什么問(wèn)題
|
|
# 3d Morse crack simulation units metal boundary s s p atom_style atomic # create geometry lattice fcc 4.05 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region box block 0 60 0 40 -0.25 0.25 create_box 5 box create_atoms 1 box mass 1 1.0 mass 2 1.0 mass 3 1.0 mass 4 1.0 mass 5 1.0 # Morse potentials pair_style morse 2.5 pair_coeff * * 1.0 1.0 2.5 neighbor 0.3 bin neigh_modify every 5 delay 0 check yes # define groups region 1 block INF INF INF 2 INF INF group lower region 1 region 2 block INF INF 38 INF INF INF group upper region 2 group boundary union lower upper group mobile subtract all boundary region leftupper block INF 20 20 INF INF INF region leftlower block INF 20 INF 20 INF INF region bonds block INF 20 19.75 20.25 -0.25 0.25 delete_atoms region bonds group leftupper region leftupper group leftlower region leftlower group halfupper1 region halfupper1 group halfupper2 region halfupper2 set group leftupper type 3 set group leftlower type 4 set group lower type 1 set group upper type 2 neigh_modify exclude type 3 4 neigh_modify exclude type 1 2 # 能量最小化 min_style sd minimize 1.0e-12 1.0e-12 10000 10000 # initial velocities 弛豫 compute new mobile temp velocity mobile create 0.001 887723 temp new timestep 0.001 thermo 200 run 10000 velocity upper set 0.0 0.01 0.0 velocity lower set 0.0 -0.01 0.0 velocity mobile ramp vy 0.01 0.01 y 20 38 sum yes # fixes fix 1 all nve fix 2 lower setforce 0.0 NULL 0.0 fix 3 upper setforce 0.0 NULL 0.0 # run thermo 200 thermo_modify temp new dump 1 all atom 1000 dump.crack.lammpstrj run 300000 |
» 搶金幣啦!回帖就可以得到:
+1/270
+1/194
+5/150
+1/84
+1/74
+1/64
+1/41
+1/34
+1/33
+1/31
+1/29
+1/15
+1/11
+1/10
+1/9
+1/7
+1/6
+1/5
+1/3
+1/3
榮譽(yù)版主 (職業(yè)作家)
|
# 3d Morse crack simulation units metal boundary s s p atom_style atomic # create geometry lattice fcc 4.05 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region box block 0 60 0 40 -0.25 0.25 create_box 5 box create_atoms 1 box mass 1 1.0 mass 2 1.0 mass 3 1.0 mass 4 1.0 mass 5 1.0 # Morse potentials pair_style morse 2.5 pair_coeff * * 1.0 1.0 2.5 neighbor 0.3 bin neigh_modify every 5 delay 0 check yes # define groups region 1 block INF INF INF 2 INF INF group lower region 1 region 2 block INF INF 38 INF INF INF group upper region 2 group boundary union lower upper group mobile subtract all boundary region leftupper block INF 20 20 INF INF INF region leftlower block INF 20 INF 20 INF INF region bonds block INF 20 19.75 20.25 -0.25 0.25 delete_atoms region bonds group leftupper region leftupper group leftlower region leftlower group halfupper1 region halfupper1 group halfupper2 region halfupper2 set group leftupper type 3 set group leftlower type 4 set group lower type 1 set group upper type 2 neigh_modify exclude type 3 4 neigh_modify exclude type 1 2 # 能量最小化 min_style sd minimize 1.0e-12 1.0e-12 10000 10000 # initial velocities 弛豫 compute new mobile temp velocity mobile create 0.001 887723 temp new timestep 0.001 thermo 200 run 10000 velocity upper set 0.0 0.01 0.0 velocity lower set 0.0 -0.01 0.0 velocity mobile ramp vy 0.01 0.01 y 20 38 sum yes # fixes fix 1 all nve fix 2 lower setforce 0.0 NULL 0.0 fix 3 upper setforce 0.0 NULL 0.0 # run thermo 200 thermo_modify temp new dump 1 all atom 1000 dump.crack.lammpstrj run 300000 |
禁蟲 (初入文壇)
|
本帖內(nèi)容被屏蔽 |
| 13 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +15 | minht 2026-03-26 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +10 | 幸運(yùn)的醬醬 2026-03-22 | 12/600 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 10/500 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過(guò) +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|