| 查看: 680 | 回复: 4 | ||
andw0wfm铁虫 (小有名气)
|
[求助]
模拟飞片和靶的的对撞问题
|
| 我想模拟一个对撞的拉伸破坏过程，但是对撞后压强和温度一直上不去，两个Velocity感觉只有一个起了作用（把两个Velocity命令换下行，效果居然差很多），In文件贴在2楼，麻烦大神帮忙看下。我知道Velocity是赋予初始速度，那这种过程需要SET FORCE来维持恒定速度嘛，如果要的话需要维持多久或者要SET FORE ALL吗？ |
回复此楼» 猜你喜欢
我在等一个没有答案的答案
已经有3人回复
-
半夜喝咖啡
已经有5人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有6人回复
-
今年审到国自然15份，谈谈感受
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
andw0wfm
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 12.3
- 帖子: 59
- 在线: 33.5小时
- 虫号: 6282844
- 注册: 2017-04-10
- 性别: GG
- 专业: 金属材料的力学行为
|
# ---------------- Initialize Simulation --------------------- clear units metal dimension 3 boundary p p p atom_style atomic timestep 0.001 # ---------------- Create Atoms --------------------- lattice fcc 4.05 origin 0.0 0.0 0.0 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region sim_box block 0 45 0 45 0 120 units lattice create_box 1 sim_box create_atoms 1 region sim_box group sim_box region sim_box # ---------------- Define Interatomic Potential --------------------- pair_style eam/fs pair_coeff * * Al_mm.eam.fs Al dump 1 all cfg 1 atomchushi*.cfg id type xs ys zs mass dump_modify 1 element Pu run 0 undump 1 # ---------------- Equilibration --------------------- reset_timestep 0 velocity all create 5 12345 mom yes rot yes dist gaussian fix equilibration all npt temp 5 5 0.1 iso 0 0 1 drag 1 variable eq equal "press/10000" variable eq1 equal "pxx/10000" variable eq2 equal "pyy/10000" variable eq3 equal "pzz/10000" thermo 1000 thermo_style custom step temp etotal ke pe lx ly lz press v_eq v_eq1 v_eq2 v_eq3 dump 2 all cfg 10000 atomequil*.cfg id type xs ys zs mass c_myPE c_myKE dump_modify 2 element Al run 10000 undump 2 unfix equilibration # ---------------- Shock --------------------- change_box all boundary p p s reset_timestep 0 fix 1 all nve velocity right set 0 0 -15 units box velocity left set 0 0 15 units box variable eq equal "press/10000" variable eq1 equal "pxx/10000" variable eq2 equal "pyy/10000" variable eq3 equal "pzz/10000" thermo 1000 thermo_style custom step temp etotal ke pe lx ly lz press v_eq v_eq1 v_eq2 v_eq3 #Use cfg for OVITO dump 3 all cfg 1000 atom*.cfg id type xs ys zs mass c_myPE c_myKE dump_modify 3 element Al dump 4 yugongniu cfg 1000 atomquyu*.cfg id type xs ys zs mass c_myPE c_myKE dump_modify 4 element Al run 40000 |
2楼2018-07-11 09:55:40
andw0wfm
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 12.3
- 帖子: 59
- 在线: 33.5小时
- 虫号: 6282844
- 注册: 2017-04-10
- 性别: GG
- 专业: 金属材料的力学行为
3楼2018-07-11 10:22:40
andw0wfm
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 12.3
- 帖子: 59
- 在线: 33.5小时
- 虫号: 6282844
- 注册: 2017-04-10
- 性别: GG
- 专业: 金属材料的力学行为
4楼2018-07-11 14:35:41
pxjdal
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 668.8
- 散金: 1
- 帖子: 91
- 在线: 100.9小时
- 虫号: 1017662
- 注册: 2010-05-13
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2018-07-12 16:10:52
精华III:
20