1 循环次数以及时间步设置
在进行第三步计算前,首先使用d.mo.isCrack=1与d.setStandarddT对裂隙记录和时间步长短进行设定。可适当扩大时间步以提高计算效率(小于4倍)。
对于位移或应力边界条件连续变化的数值模拟,需将该过程进行微分。设置变量totalCircle和stepNum两层循环嵌套进行计算。其中totalCircle为中间文件保存次数,而totalCircle*stepNum则为总的平衡迭代次数。在确定的位移或应力总变化量的情况下,totalCircle*stepNum的值越大,单步位移或应力变化量越小,数值模拟越精确。
而对于微波破岩来说,将半径增量(elementExpandRate)乘以单元的平均半径(d.mo.aR),并除以总循环次数(totalCircle*stepNum),得到单次迭代计算时单元半径膨胀量。单次膨胀量与劲度系数的乘积等于单步单元膨胀后产生的膨胀力增量,其值需小于单元间的抗拉力。