第五讲_分子动力学模拟Lammps实现.ppt

B、使用L-J势模拟裂纹的扩展;B、使用L-J势模拟裂纹的扩展;Lammps计算输入文件;Lammps计算输入文件;Lammps计算输入文件;Lammps计算输入文件;Lammps计算输入文件;Lammps计算输入文件;C.使用EAM势模拟Ni的循环行为;金属材料模拟中Lammps的单位;# 3d metal 拉伸模拟 #模拟条件的初始化 Units metal #指定模拟中的单位类型 boundary s s s #指定模拟的边界条件 atom_style atomic #指定原子类型，原子的属性 Lattice fcc 3.52 #指定材料的晶格类型和晶格常数 region box block -5 5 -5 5 -15 15 #xlo,xhi,ylo,yhi,zlo,zhi =区域box的上下限 create_box 1 box #建立只有1中原子类型的simulation box create_atoms 1 box #在simulation box中创建类型为1的原子（原子位置初始化）;pair_style eam #定义原子之间相互作用势 pair_coeff * * Ni_u3.eam #定义作用势的参数，通过文件输入 Neighbor 2.0 bin #定义近邻表参数(r+rc),rc=2.0埃 neigh_modify delay 5 #定义更新近邻表的加载步间隔;region lower block INF INF INF INF INF -13 #生成lower区域，便于加载 region upper block INF INF INF INF 13 INF #生成upper区域，便于加载 group lower region lower #定义lower组 group upper region upper #定义upper组 group boundary union lower upper #定义boundary组 group mobile subtract all boundary #定义mobile组，便于计算统计温度（调温）;#定义计算和边界钢化 compute new2 all centro/atom #定义计算原子的适配度（针对所有的原子） compute new3 mobile temp 定义计算统计温度（针对可动原子） Fix 1 all nve #时间积分算法（微正则系综） fix 2 lower setforce NULL NULL 0.0 #使得lower group之间没有z方向的力，z方向钢化 fix 3 upper setforce NULL NULL 0.0 #使得upper group之间没有z方向的力，z方向钢化;#定义控温和原子速度初始化 fix 4 mobile temp/rescale 10 30.0 30.0 5.0 1.0 #使用直接调温法调整mobile原子的速度，每10步一调 #使得目标温度为30+-1K，温度在30+-5rescale之外调 fix_modify 4 temp new3 #调温时，计算温度采用news计算 Velocity mobile create 30 887723 temp new3 #初始化mobile group中原子的初速度 #通过温度30K初始化，温度通过new3计算，887723为随机数 ;dump 2 all custom 200 d_t.txt tag x y z c_new2 #定义输出，每两百步输出一次 #tag为原子标号 #x，y，z为原子坐标 #c_new2为new2所计算的???子适配度 #这一参数可以判断原子是否是理想晶格 #或者在位错和堆垛层错上 timestep 0.001 #时间步长间隔（ps=10e-12秒） #时间步长在飞秒量级10e-15秒 #小于原子自振周期的十分之一（0.1ps/10）;Lammps计算输入文件;#弛豫和循环加载 run 2000 #不加任何载荷计算2000步 velocity upper set 0 0 0.4 velocity lower set 0 0 -0.4 run 10000 velocity upper set 0 0 -0.4 velocity lower set 0 0 0.4 run 20000 velocity upper set 0 0 0.4 velocity lower set 0 0 -0.4 run 10000 #循环载荷模拟;D.使用EAM势模拟Ni的剪切行为;D.使用EAM势模拟Ni的剪切行为;# 3d metal shear simulation 3维金属剪切模拟 Units metal #采用金属材料单位