| 查看: 2892 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
有人做裂纹模拟的吗？看看我这个裂纹的in文件有什么问题
|
|||
|
# 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 |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有10人回复
-
基金系统什么内容也没有
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助小裂纹模拟软件或者如何模拟
已经有11人回复
- 溶剂热烧出的晶体总是很小，又有裂纹，不知道咋办？感谢！
已经有15人回复
- 大家知道哪个学校可以给三点弯曲样品做预制裂纹吗？价格呢？
已经有4人回复
- 求助关于包晶反应钢裂纹的问题
已经有5人回复
- 有关氧化铝陶瓷烧结后裂纹的现象，请大家分析原因
已经有23人回复
- CVD涂层后表面裂纹的问题
已经有6人回复
- 请问碳化硅陶瓷烧结产生裂纹是什么原因？
已经有14人回复
- 【求助】ABS透明管打压会有小裂纹
已经有4人回复
- 【讨论】什么是小裂纹
已经有7人回复
- 【原创】ALOF软件-三维疲劳裂纹扩展模拟分析的实例
已经有20人回复
- 【交流】关于材料脆性断裂裂纹的问题
已经有11人回复
- 【求助】铸造裂纹及其他问题
已经有23人回复
- 【讨论】有没有做疲劳裂纹扩展速率研究的啊
已经有18人回复
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zh1987hs: 金币+3, 谢谢 2012-10-01 11:02:11
小木虫: 金币+0.5, 给个红包，谢谢回帖
zh1987hs: 金币+3, 谢谢 2012-10-01 11:02:11
|
# 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 |
6楼2012-09-29 11:51:41
2楼2012-09-05 10:11:52
3楼2012-09-27 15:19:20
4楼2012-09-27 19:51:35
投票:
