lammps计算聚合物例子_lammps案例分享：聚乙烯拉伸过程分子动力学模拟

本文分享一个聚乙烯拉伸的案例代码​。复合物的拉伸过程和金属的拉伸类似，​差别在于势函数的设置。

金属原子直接没有键连接，因此，在进行金属材料拉伸时，只要指定势函数文件名和指定原子类型即可。

聚合物则需要设置键、键角、二面角等势函数参数​。

具体模拟过程如下：

（1）在MS中生成50个单体、10个链的聚乙烯​，设置cvff力场，导出car文件。

生成的聚乙烯如下图所示：

（2）运行msi2lmp pp -class I -frv cvff >data.pp命令，生成pp.data聚乙烯结构数据文件。

（3）在in.pp中设置势函数参数，读取pp.data，驰豫后进行拉伸​。

in.pp代码​如下：

#聚乙烯拉伸
# 模型初始化
units		real
boundary	p p p
atom_style	full

#设置键势类型
bond_style      harmonic
#设置键角势类型
angle_style     harmonic
#设置二面角势类型
dihedral_style	harmonic
#设置对势类型
pair_style	lj/cut 10.5

#读取聚乙烯结构数据
read_data	pp.data

#设置邻居参数
neighbor