| 5 | 1/1 | 返回列表 |
| 查看: 1153 | 回復(fù): 5 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
in 文件,裂紋拉不開.求助深入解讀。 已有1人參與
|
|||
|
# 3d Morse crack simulation units metal boundary p p p atom_style atomic neighbor 0.3 bin neigh_modify every 5 delay 0 check yes # 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 26.9 mass 2 26.9 mass 3 26.9 mass 4 26.9 mass 5 26.9 # Morse potentials pair_style morse 2.5 pair_coeff * * 1.0 1.0 2.5 # 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 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 # minimize min_style sd minimize 1.0e-12 1.0e-12 10000 10000 # initial velocities relax compute new mobile temp velocity mobile create 750.0 887723 temp new timestep 0.001 thermo 500 run 10000 velocity upper set 0.0 0.04 0.0 sum yes velocity lower set 0.0 0.0 0.0 sum yes velocity mobile ramp vy 0.0 0.04 y 2 38 sum yes # fixes fix 1 lower setforce 0.0 0.0 0.0 fix 2 upper setforce 0.0 0.0 0.0 fix 3 all nvt 750.0 750.0 0.01 fix 4 all temp/rescale 100 750.0 750.0 3 0.01 fix 5 mobile deform 1000 y erate 0.003 units box fix 6 mobile setforce 0.0 NULL 0.0 compute 1 all stress/atom compute 2 all temp compute 3 all reduce sum c_1[2] compute 4 all pe compute 5 all ke variable 2 equal c_3 variable 3 equal c_2 variable 4 equal c_4 variable 5 equal c_5 fix 7 all ave/time 1 100 1000 v_2 v_3 v_4 v_5 file a_stress.txt thermo_style custom step temp pe etotal press vol thermo 200 thermo_modify temp new dump 1 all atom 1000 dump.crack.lammpstrj min_style sd minimize 1.0e-12 1.0e-12 10000 10000 run 200000 |
木蟲 (正式寫手)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305求調(diào)劑 +7 | RuiFairyrui 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料學(xué)碩333求調(diào)劑 +10 | 北道巷 2026-03-24 | 10/500 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +5 | 舉個栗子oi 2026-03-24 | 6/300 |
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分數(shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開頭均可,過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|