姚小贰 liujingzhou629
关注数: 43 粉丝数: 26 发帖数: 1,047 关注贴吧数: 17
用lammps拉伸一直出问题有大神能看下我的程序吗 units metal dimension 3 boundary p p p atom_style full read_data Fe.data lattice fcc 3.59 pair_style eam/alloy pair_coeff * * FeNiCr_Bonny_2013_ptDef.eam.alloy Cr Ni Fe neighbor 2.0 bin #截断半径 neigh_modify delay 1 check yes #邻近原子列表更新速度 # define groups----------------定义组 region 1 block INF INF INF INF INF 10 #定义一个叫1的区域 group lower region 1 #此区域里的原子叫lower region 2 block INF INF INF INF 103 INF #定义一个叫2的区域 group upper region 2 #此区域里的原子叫upper group boundary union lower upper #定义lower+upper = boundary group mobile subtract all boundary #定义mobile=all- boundary #initial velocities-------------设置原子初速度 velocity lower set 0.0 0.0 0.0 #设置left区域原子初速度为0 velocity upper set 0.0 0.0 0.3 #设置right区域z方向原子速度为0.3 velocity mobile ramp vz 0.0 0.3 z 3 17 sum yes min_style sd #能量最小模式 minimize 1.0e-6 1.0e-6 1300 10000 #驰豫 fix 1 lower setforce 0.0 NULL NULL #固定下边的原子 fix 3 upper setforce 0.0 NULL NULL fix 2 all npt temp 1300 1300 1 x 0 0 1 y 0 0 1 z 0 0 1 drag 1 thermo 200 #每200步输出一次结果 timestep 0.003 #时间步长0.003ps run 3000 #运行5000步 unfix 2 fix 3 all npt temp 1300 1300 1 x 0 0 1 y 0 0 1 drag 1 fix 4 all temp/rescale 1300 1.0 1.0 5 1.0 fix 5 all deform 1000 z erate 0.002 units box #设置沿z方向拉伸,应变率为0.002(1/ps) dump 1 all atom 100 dump1.txt #输出结果到dump.txt min_style sd #能量最小模式 minimize 1.0e-6 1.0e-6 1300 10000 #驰豫 run 200000 #运行10000步
1 下一页