| 查看: 596 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
andw0wfm铁虫 (小有名气)
|
[求助]
模拟飞片和靶的的对撞问题
|
|
| 我想模拟一个对撞的拉伸破坏过程，但是对撞后压强和温度一直上不去，两个Velocity感觉只有一个起了作用（把两个Velocity命令换下行，效果居然差很多），In文件贴在2楼，麻烦大神帮忙看下。我知道Velocity是赋予初始速度，那这种过程需要SET FORCE来维持恒定速度嘛，如果要的话需要维持多久或者要SET FORE ALL吗？ |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
心脉受损
已经有6人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
pxjdal
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 668.8
- 散金: 1
- 帖子: 91
- 在线: 100.9小时
- 虫号: 1017662
- 注册: 2010-05-13
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2018-07-12 16:10:52
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