|
|
[交流]
LAMMPS建模后无法熔化是为什么?
请问各位大佬一个问题,我在用LAMMPS建了一个三元合金球形颗粒模型之后,给它加热到3000K他依然没有怎么熔化,这是为什么? 同样的合金势函数,如果是一个立方体系,它在2100K的时候就熔化了。但换成这种球形体系就熔化不了,这是为什么呢? in文件如下:
lattice fcc 4.05
region spher1 sphere 0 0 0 10
region spher2 sphere 20.5 0 0 10
region spher3 sphere 40.7 0 0 10
region spher union 3 spher1 spher2 spher3
create_box 3 spher
create_atoms 1 region spher1
create_atoms 2 region spher2
create_atoms 3 region spher3
thermo_modify lost ignore
pair_style eam/alloy
pair_coeff * * Farkas_Nb-Ti-Al_1996.eam.alloy Nb Ti A
timestep 0.0001
thermo 100
# ---------- Define Interatomic Potential -----------
neighbor 0.5 bin
neigh_modify every 1 delay 0 check yes
min_style cg
minimize 1.0e-6 1.0e-8 1000 1000
run 100
velocity all create 300 49284 dist gaussian loop local
min_style sd
minimize 1.0e-12 1.0e-12 10000 100000
fix 1 all nve
fix 2 all langevin 300 300 $(100.0*dt) 699483
dump 2 all custom 2000 simu1dump_* id type x y z c_coord
run 100000
unfix 1
unfix 2
undump 2
thermo_modify lost ignore
fix 3 all npt temp 300 3000 $(100.0*dt) y 50 50 1 z 50 50 1 drag 0.2
dump 3 all custom 2000 simu1dump_* id type x y z c_coord
run 20000
unfix 3
undump 3
thermo_modify lost ignore
fix 31 all npt temp 3000 3000 $(100.0*dt) y 50 50 1 z 50 50 1 drag 0.2
dump 31 all custom 1000 simu1dump_* id type x y z c_coord
run 20000
unfix 31
undump 31
![LAMMPS建模后无法熔化是为什么?]()
三元模型结束.png
![LAMMPS建模后无法熔化是为什么?-1]()
三元球形初始.png |
» 猜你喜欢
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有3人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
|