| 查看: 1120 | 回复: 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 |
回复此楼» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有11人回复
-
限项规定
已经有9人回复
-
免疫学博士有名额，速联系
已经有4人回复
-
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有5人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助，合金开坯表面总是出现裂纹........
已经有13人回复
- 【求助】晶面有裂纹 可以测试单晶不？
已经有3人回复
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
kite1755楼2013-01-04 10:02:32
2楼2013-01-03 17:06:56
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
3楼2013-01-03 19:46:13
4楼2013-01-04 09:34:08