24小时热门版块排行榜

返回列表

andw0wfm

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 12.3
帖子: 59
在线: 33.5小时
虫号: 6282844
注册: 2017-04-10
性别: GG
专业: 金属材料的力学行为

[求助] 模拟飞片和靶的的对撞问题

我想模拟一个对撞的拉伸破坏过程，但是对撞后压强和温度一直上不去，两个Velocity感觉只有一个起了作用（把两个Velocity命令换下行，效果居然差很多），In文件贴在2楼，麻烦大神帮忙看下。我知道Velocity是赋予初始速度，那这种过程需要SET FORCE来维持恒定速度嘛，如果要的话需要维持多久或者要SET FORE ALL吗？

回复此楼

» 猜你喜欢

烦躁已经有5人回复
寻找合作团队已经有4人回复

1楼 2018-07-11 09:53:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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