在之前的专栏中，曾介绍过动态统计某一个固定区域内原子数量的方法：
也介绍过动态统计某一个原子周围原子数量的方法：
下面介绍第三种类型：动态统计某一个原子周围不同类型原子数量的方法。
以小球的随机碰撞为例，原子类型为1的红色小球（id也为1）周围分布原子类型为2、3的蓝色和黄色小球。
分别统计红色小球周围蓝色原子和红色原子的数量，需要用到count(type,region)命令，改命令可统计区域region内type类型原子的数量。
对应的代码为：
fix 1 all nvegroup type2 type 2group type3 type 3variable x1 equal x[1]variable y1 equal y[1]variable z1 equal z[1]region ss sphere v_x1 v_y1 v_z1 15group mobile dynamic all region ss every 10variable number_all equal count(mobile)variable number_type2 equal count(type2,ss)variable number_type3 equal count(type3,ss)thermo_style custom step v_number_all v_number_type2 v_number_type3thermo 100run 10000
统计结果如下图所示：
报名一对一培训课程可领完整版《lammps中文教程》，配套学习更轻松！
扫码领取粉丝专属优惠券↓↓↓
