在前面的推文中，介绍过lammps动态统计区域内原子数量、区域压力的方法。
本文介绍统计某一个区域内原子受力的方法。
group：原子组id dim：力的方向，可设置x、y、z region：区域id
对应的计算代码：
region mid block -20 20 -20 20 -20 20variable fxx equal fcm(all,x,mid)thermo_style custom step v_fxx
先设置中间区域mid，调用fcm()命令计算，最后用thermo命令输出结果。
全部in文件代码：
#气体粒子数量variable npart equal 1000#模拟基本参数设置units ljdimension 3atom_style atomicboundary p p pneighbor 0.5 binneigh_modify every 1 delay 0 check yes#三维box尺寸region box block -60 60 -20 20 -20 20#生成boxcreate_box 2 box#随机产生100个气体原子create_atoms 1 random ${npart} 324523 box#设置摩尔质量1mass * 1#力场设置，soft势pair_style soft 1.0pair_coeff * * 10.0#能量最小化minimize 1e-4 1e-4 1000 1000reset_timestep 0#温度初始化velocity all create 2.0 34234123 dist gaussian#保存轨迹文件dump img all atom 100 gas2.xyz#热力学输出thermo_style custom step temp ke pethermo 100#模拟步长timestep 0.001#设置nve系综fix 1 all nveregion mid block -20 20 -20 20 -20 20variable fxx equal fcm(all,x,mid)thermo_style custom step v_fxxrun 1000
---------- End -----------
报名一对一培训课程可领完整版《lammps中文教程》，配套学习更轻松！
扫码领取粉丝专属优惠券↓↓↓
