| 查看: 2650 | 回复: 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 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助小裂纹模拟软件或者如何模拟
已经有11人回复
- 溶剂热烧出的晶体总是很小，又有裂纹，不知道咋办？感谢！
已经有15人回复
- 大家知道哪个学校可以给三点弯曲样品做预制裂纹吗？价格呢？
已经有4人回复
- 求助关于包晶反应钢裂纹的问题
已经有5人回复
- 有关氧化铝陶瓷烧结后裂纹的现象，请大家分析原因
已经有23人回复
- CVD涂层后表面裂纹的问题
已经有6人回复
- 请问碳化硅陶瓷烧结产生裂纹是什么原因？
已经有14人回复
- 【求助】ABS透明管打压会有小裂纹
已经有4人回复
- 【讨论】什么是小裂纹
已经有7人回复
- 【原创】ALOF软件-三维疲劳裂纹扩展模拟分析的实例
已经有20人回复
- 【交流】关于材料脆性断裂裂纹的问题
已经有11人回复
- 【求助】铸造裂纹及其他问题
已经有23人回复
- 【讨论】有没有做疲劳裂纹扩展速率研究的啊
已经有18人回复
» 抢金币啦！回帖就可以得到:
-
坐标广州，征女友
+2/96
-
欧洲“CO2捕集与加氢催化转化”博士招聘（捷克-德国-斯洛文尼亚三国项目）
+1/89
-
盐城师范学院生物与医药（生物技术与工程领域）还有调剂名额
+1/82
-
工学门类（前两位代码08）满足国家考研A区工科门类分数线要求的同学看过来
+1/49
-
江西科技师范大学086生物与医药
+1/37
-
盐城师范学院生物与医药（生物技术与工程领域）还有调剂名额
+1/36
-
江西科技师范大学086生物与医药
+1/36
-
中新绿色科技研究院（山东济南）招收机械设计方向实习生
+1/34
-
爱情短剧
+1/23
-
上海交通大学-电气工程学院-高温超导磁体-2026年秋季入学全日制博士生招聘
+1/21
-
求助请教各位:如果您的算法被顶尖院校学术不端怎么办？
+1/11
-
AIRS微纳机器人中心诚招研究型硕士生（Mphil）、全职科研助理（RA）、实习生（Intern�
+1/10
-
河南省医学科学院王宁利院士/苗彦彦研究员科研团队、伏羲慧眼云平台团队2026年博士、�
+1/9
-
北化大学材料科学与工程学院研究生调剂
+1/7
-
青岛科技大学高分子学院硕士招生
+1/7
-
香港科技大学用于先进LCD器件的偏振发射钙钛矿纳米棒增强薄膜
+2/6
-
材料与化工专业还有大量指标，欢迎大家调剂，调剂系统已经开放
+1/4
-
【每日皆可乐】探秘细胞因子：趋化因子家族解析
+1/1
-
申请临床医学博士、医学博士
+1/1
-
招07理学调剂生，剩余2个名额
+1/1
★ ★ ★ ★
小木虫: 金币+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
