分子动力学模拟的Lammps实现

《分子动力学模拟的Lammps实现》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、B、使用L-J势模拟裂纹的扩展L-J无量纲化计算中的单位B、使用L-J势模拟裂纹的扩展裂纹lowerupperleftupperleftlowerxyLammps计算输入文件#2dLJcracksimulation(问题的基本初始化)dimension2#2维分子动力学模拟boundaryssp#x,y方向不存在周期边界条件，z方向周期边界条件。atom_styleatomic#原子类型（金属）Neighbor0.3bin#建近邻表参数（rc+0.3）bin表示为近邻表类型。neigh_modifyde

2、lay5#间隔多少载荷步重新形成近邻表Lammps计算输入文件#creategeometry创建初始几何构形Latticehex0.93#指定晶格类型（二维hex）和晶格常数Regionboxblock0100040-0.250.25#定义一个区域create_box5box#在指定区域建立一个simulationbox,5表示原子类型的种类数create_atoms1box#在simulationbox中创建类型为1的原子（原子位置初始化）Mass11.0Mass21.0Mass31.0Mass41.

3、0Mass51.0#指定单个原子的质量Lammps计算输入文件#LJpotentials(指定原子作用势)pair_stylelj/cut2.5#指定lj势，截断半径为2.5pair_coeff**1.01.02.5#指定lj势参数#definegroups（便于加载）Region1blockINFINFINF1.25INFINFGrouplowerregion1#定义lower组（便于施加外加速度）Region2blockINFINF38.75INFINFINFGroupupperregion2#定义

4、upper组（便于施加外加速度）Groupboundaryunionlowerupper#定义总边界组Groupmobilesubtractallboundary#定义可动原子组（便于统计温度）Lammps计算输入文件regionleftupperblockINF2020INFINFINFregionleftlowerblockINF20INF20INFINFgroupleftupperregionleftuppergroupleftlowerregionleftlower#定义左上、左下原子组（便于指

5、定裂纹的存在）setgroupleftuppertype2setgroupleftlowertype3setgrouplowertype4setgroupuppertype5#指定原子类型（便于指定裂纹的存在）Lammps计算输入文件#initialvelocities初始化速度computenewmobiletemp#定义温度的计算（可动区域内统计平均）computenew2mobilestress/atom#定义原子应力的计算（整个区域）Velocitymobilecreate0.01887723t

6、empnew#按指定的温度（0.01）计算方法，初始化原子的速度Velocityupperset0.00.30.0#upper原子组y方向的速度为0.3Velocitymobilerampvy0.00.3y1.2538.75sumyes#mobile原子的速初始度从0到0.3线性变化#fixes施加约束fix1allnve#nve系综的积分算法fix2boundarysetforceNULL0.00.0#边界boundary上力条件，钢化原子，便于加载！！Lammps计算输入文件#run运行计算time

7、step0.003#时间间隔步Thermo200#每200步输出热动力学统计量thermo_modifytempnew#计算温度通过new指示的方法计算neigh_modifyexcludetype23#原子2，3之间作用取消（也就是通过不使他们在近邻表中出现实现）Dump1allatom500dump.crack#每隔500步将原子信息写入文件dump.crackDump2mobilecustom500dump2.cracktagxyzc_new2[2]run5000#进行5000步的模拟C.使用EA

8、M势模拟Ni的循环行为uuuzxy常温30K条件下金属材料模拟中Lammps的单位#3dmetal拉伸模拟#模拟条件的初始化Unitsmetal#指定模拟中的单位类型boundarysss#指定模拟的边界条件atom_styleatomic#指定原子类型，原子的属性Latticefcc3.52#指定材料的晶格类型和晶格常数regionboxblock-55-55-1515#xlo,xhi,ylo,yhi,zlo,zhi=区域box的上下限c