提供一个完整的lammps金属凝固in文件代码
时间: 2024-03-23 07:40:19 浏览: 145
以下是一个简单的金属凝固的LAMMPS输入文件代码:
```
#定义模拟体系
units metal
dimension 3
boundary p p p
atom_style atomic
#定义原子类型
lattice fcc 3.52
region box block 0 10 0 10 0 10
create_box 1 box
create_atoms 1 box
#定义相互作用势能
pair_style eam/alloy
pair_coeff * * Ni_u3.eam
#定义温度控制
velocity all create 300.0 12345
#定义时间步长
timestep 0.001
#定义运行时间
run 5000
#定义凝固过程
fix 1 all nvt temp 300.0 300.0 0.1
fix 2 all langevin 300.0 300.0 100.0 48279
compute 1 all temp
compute 2 all pe
compute 3 all ke
compute 4 all stress/atom NULL
variable c1 equal vol
variable c2 equal pxx
variable c3 equal pyy
variable c4 equal pzz
variable c5 equal c2+c3+c4
thermo_style custom step temp c_1 c_2 c_3 c_4 c_5 pe ke etotal press lx ly lz
thermo_modify format float %20.10f
#定义输出文件
dump 1 all atom 100 dump.lammpstrj
```
这个输入文件使用了EAM势能模型,模拟了一个简单的fcc结构的金属体系的凝固过程。通过定义原子类型,相互作用势能和温度控制等参数,并且通过定义fix、compute和thermo等命令来控制凝固过程的温度、压力、能量等参数,并输出到文件中进行分析。需要注意的是,具体的输入文件需要根据实际情况进行调整和修改。
阅读全文