流固耦合法模拟岩石爆破时耦合范围的确定

为研究岩石爆破问题模拟中的流固耦合算法以及提高计算效率,文章采用LS-DYNA软件开展大量计算,分别讨论了不同装药半径(r)、不耦合系数(K)和网格精度条件下流固耦合区域半径(R)对单元峰值应力衰减曲线的影响,并结合应力云图和应力时程曲线分析,确定出合理的流固耦合范围。结果表明:随着耦合装药半径r的增大,最小合理流固耦合半径R_(min)与r比值呈减小趋势;相同装药半径下随着K增大,R_(min)/r逐渐增大;网格精度对流固耦合半径取值影响较小;相同条件下,当R较小时,应力云图不稳定且应力时程曲线后期会迅速衰减,适当增大流固耦合半径会逐渐消除该现象。对于一般的岩石爆破问题,建议R/r取5～10。