| 查看: 2898 | 回复: 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 |
» 猜你喜欢
有多少人是今天查系统知道结果的?
已经有19人回复
找导师
已经有7人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
有没有仍没收到信息的
已经有7人回复
国社科又开始会评了,不知道这次命运如何
已经有4人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
» 本主题相关价值贴推荐,对您同样有帮助:
求助小裂纹模拟软件或者如何模拟
已经有11人回复
溶剂热烧出的晶体总是很小,又有裂纹,不知道咋办?感谢!
已经有15人回复
大家知道哪个学校可以给三点弯曲样品做预制裂纹吗?价格呢?
已经有4人回复
求助关于包晶反应钢裂纹的问题
已经有5人回复
有关氧化铝陶瓷烧结后裂纹的现象,请大家分析原因
已经有23人回复
CVD涂层后表面裂纹的问题
已经有6人回复
请问碳化硅陶瓷烧结产生裂纹是什么原因?
已经有14人回复
【求助】ABS透明管打压会有小裂纹
已经有4人回复
【讨论】什么是小裂纹
已经有7人回复
【原创】ALOF软件-三维疲劳裂纹扩展模拟分析的实例
已经有20人回复
【交流】关于材料脆性断裂裂纹的问题
已经有11人回复
【求助】铸造裂纹及其他问题
已经有23人回复
【讨论】有没有做疲劳裂纹扩展速率研究的啊
已经有18人回复
» 抢金币啦!回帖就可以得到:
云南大学丁军桥团队2027年博士招生
+1/1985
面上中了,散金币。
+5/1825
八年奋战,面上中了,同时一个读了8年的博士(中途接手)答辩通过了,都是涉险过关
+2/470
散金币祈福
+5/358
货币银行学通俗教材版(一部超越哈耶克的经济学)
+1/252
【毕业前夕,畅谈科研】文集
+1/164
原子层沉积(ALD)磁控溅射PECVD等微纳代工服务:18817872921
+1/97
厦门大学崔明旭课题组招收2027届推免硕士/直博与申请考核博士生(金属有机、高分子)
+1/33
无穷无尽的科研生活焦虑
+4/12
招收CSC博士生(2027年9月入学)
+1/9
香港科技大学(广州) 强化学习、机器人、具身智能方向等方向硕博-27 spring fall
+1/8
南京航空航天大学激光焊接与再制造技术研究所招收2027年博士研究生
+1/7
2026年工程管理与智慧城市国际会议 (EMSC 2026)
+1/4
厦大化院海优团队招收2027级【博士】和【推免硕士】研究生各2名
+1/3
专利审查意见答复,是AI辅助生成中很难的一个功能...
+1/2
寻求合作【出版代理】兼职推广|无需坐班|一单5000元
+1/2
北航-国家级人才团队-招2027年9月博士
+1/1
中国科学院杭州医学研究所-科研助理,博士后及2027级博士招聘
+1/1
寻求合作【大学教师】兼职出版顾问|无需坐班|一单5000元
+1/1
上海大学优青团队招聘硕士、博士研究生、博士后、实验员以及相关教职人员
+1/1
2楼2012-09-05 10:11:52
3楼2012-09-27 15:19:20
4楼2012-09-27 19:51:35
5楼2012-09-28 08:49:08
★ ★ ★ ★
小木虫: 金币+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
7楼2012-10-05 10:37:35
8楼2014-11-12 14:07:34
9楼2014-11-12 15:17:47
10楼2015-04-01 23:46:19
11楼2018-10-19 16:21:54
12楼2018-10-19 16:22:39
★
小木虫: 金币+0.5, 给个红包,谢谢回帖
小木虫: 金币+0.5, 给个红包,谢谢回帖
|
本帖内容被屏蔽 |
13楼2018-10-22 12:51:03









回复此楼