1 模型平衡状态检查
1) 使用d.status.dispEnergy()查看当前系统中的能量,包括Kinetic Energy(动能)、Elastic Energy(弹性势能)、Gravity Energy(重力势能)、Total Energy(总能量)和Heat(热量)。当动能远小于弹性势能时候,模型平衡良好。
2) 通过d.show()命令显示整个数值模拟的状态,总共包含4张图,分别为StressZZ(Z方向应力分布)、Energy Conversion(能量转化曲线)、Normal force of boundary(边界受力曲线)、Heat generation(热量生成曲线)。若数值模型平衡良好,则在各类能量曲线与受力曲线的末端会有较长的水平段。由于应力波会在试样中反射传递,边界受力的平衡需要较长时间,当边界受力曲线都已呈完全水平时,模型达到平衡。
3) 在后处理模块中,观察模型中颗粒运动速度,当颗粒速度达到极小值(如1e-10m/s)时,可将单元视为静止,模型平衡良好。实际上单元很难达成完全静止。