24小时热门版块排行榜

返回列表

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 cuiziqi 的 39 个金币，回帖就立即获得 1 个金币，每人有 1 次机会

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cuiziqi

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 66
帖子: 70
在线: 7.3小时
虫号: 14359079

[交流] 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

三元模型结束.png