在做室内三轴试验时,通常需要保持较小的加载速率以保持试样的拟静态条件。那么在使用liggghts软件模拟三轴试验时,如何去控制试样保持所谓的“拟静态”呢?
主要思路就是通过颗粒的旋转动能来判断。由于所模拟的颗粒是球形颗粒,所以选用erotate/sphere命令去计算颗粒的旋转动能。(ligghts中还有针对不同颗粒的动能计算命令,如erotate/multisphere是用于计算颗粒簇动能,具体可以参考官网4. Commands — LIGGGHTS v3.X documentation)
在erotate/sphere命令中,计算球形颗粒旋转动能的公式为:
I表示转动惯量,ω表示粒子的角速度。
语法也非常简单:
compute ID group-ID erotate/sphere general_keywords general_values halfstep
其中,ID是用户给此计算命令所命名(可以随意),group-ID是用户先前所分配粒子组(通常选择all)。
示范:
compute 1 all erotate/sphere
这个命令的输出值为一个标量值KE,那么我们如何去输出它呢?只需要在thermo_style命令中去设置KE的输出即可。
示范:
thermo_style custom ke
这样,动能就会输出在thermo文件中了。有关thermo的命令可以参考:thermo_style command — LIGGGHTS v3.X documentation