LIGGGHTS通过fix freeze命令固定颗粒保持不动

LIGGGHTS手册中关于fix freeze命令的用法和解释如下：fix freeze command — LIGGGHTS v3.X documentation (cfdem.com)

其主要作用时将颗粒上的力和力矩归零。

通过fix freeze命令结合group命令，可以实现让某些颗粒保持固定不动。比如：

#把颗粒group，以便freeze
group   fixParticle id 1:1000

# cfd coupling
fix     cfd all couple/cfd couple_every 10 mpi#50 mpi#
fix     cfd2 all couple/cfd/force

# apply nve integration to all particles that are inserted as single particles
fix         integr all nve/sphere

#对group以后的颗粒进行freeze
fix     fixPar fixParticle freeze

但是需要注意的是，在通过CFDEMcoupling进行耦合计算的时候，一定要注意，把fix freeze命令放在fix couple/cfd以及fix couple/cfd/force之后，否则freeze在一起的颗粒group在其自身受到的力和扭矩被清零之后，还会受到流体作用力，而一起成组的运动。