| 查看: 642 | 回复: 4 | ||
andw0wfm铁虫 (小有名气)
|
[求助]
模拟飞片和靶的的对撞问题
|
| 我想模拟一个对撞的拉伸破坏过程，但是对撞后压强和温度一直上不去，两个Velocity感觉只有一个起了作用（把两个Velocity命令换下行，效果居然差很多），In文件贴在2楼，麻烦大神帮忙看下。我知道Velocity是赋予初始速度，那这种过程需要SET FORCE来维持恒定速度嘛，如果要的话需要维持多久或者要SET FORE ALL吗？ |
回复此楼» 猜你喜欢
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
-
石河子大学（211、双一流）硕博研究生长期招生公告
已经有3人回复
-
一志愿西安交通大学材料工程专业 282分求调剂
已经有13人回复
-
北科281学硕材料求调剂
已经有6人回复
-
328求调剂，英语六级551，有科研经历
已经有10人回复
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
