详细内容

LAMMPS讲解33-水泥体系的钙矾石拉伸的LAMMPS模拟

使用ReaxFF研究水泥基体系是其重要的应用,可以研究其腐蚀老化等化学反应机理。本例发布钙矾石的拉伸性能模拟。in文件和data文件以及力场文件加微信lmp_zhushou获取。文件未经过测试,欢迎反馈问题。

首先是建模过程

打开MS,点击import

image.png 

在打开的对话框中,分别打开Structures->minerals->7-Sulfates->ettringite

image.png 

image.png 

image.png 

image.png 

执行操作,打开Redefine Lattice

image.png 

点击Redefine后关闭对话框

image.png 

 

打开Supercell对话框

image.png 

按需要修改扩胞数量,这里设置都是4,然后点击Create Spuercell

image.png 

点击补氢按钮

image.png 

打开Lattice Parametersγ修改为90.000a修改为原来值乘以cos30,修改前要确保取消Keep fractional。。。的勾选

image.png 

image.png 

建完的模型如下图,超出盒子的不是会由周期性边界条件显示的原因。

image.png 

将模型导出为cif格式

image.png 

用文本编辑器打开保存的cif文件拷贝盒子尺寸和坐标,在Excel中组合为data文件,记得要将坐标按照盒子三个方向上的尺寸进行缩放,然后将内容拷出到另一个文本文件建立data文件。下面就是建立的模型

image.png 

in文件内容

units    real

atom_style  charge

neighbor  2.0 bin

neigh_modify  every 1 delay 0 check no

 

read_data  system.data

 

 

pair_style  reaxff NULL

pair_coeff * * ff.reaxff C H O Ca Si Al S

 

timestep  0.1

velocity all create 300.0 4928459 rot yes dist gaussian

thermo 100

thermo_style custom step temp pe ke

thermo_modify flush yes

 

fix 2 all npt temp 300 300 10 iso 0 0 100

fix             3 all qeq/reaxff 1 0.0 10.0 1e-6 reaxff

 

fix 1 all reaxff/bonds 1000 bonds.reaxff

dump 1 all custom 1000 dump.lammpstrj id type q x y z

 

run    40000

 

unfix 2

fix 2 all npt temp 300 300 10 y 0 0 100 z 0 0 100

fix 3 all deform 1 x erate 0.001 remap x

 

run 1000000

 

 

感谢鲍路瑶老师的分享,内容来自于鲍老师分享出来的资料

如有需要添加微信:lmp_zhushou  进入微信群,帮助他人,共建社区

获取完整版lammps讲义可以加微信lmp_zhushou或加入QQ994359511

 

 


最新评论
请先登录才能进行回复登录
技术支持: CLOUD | 管理登录
seo seo