周期冲击载荷在采矿、土木等工程中较为常见.在周期冲击荷载作用下,巷道围岩会在极短时间内产生大量裂缝,块体会以猛烈方式脱离巷道围岩,造成岩爆、垮塌等灾害.自主开发的拉格朗日元和离散元耦合的连续-非连续方法主要包括4个计算模块:应力应变模块、节点分离模块、接触力求解模块和运动方程求解模块.鉴于单元畸变和局部自适应阻尼可能导致单元弹射出模型时速度失真,为了准确模拟弹射现象,在该连续-非连续方法的基础上,对弹射单元进行了刚化处理.该方法的实质是确保弹射单元在碰撞前后作刚体运动,这是通过对节点速度取平均、消除单元应力实现的.对周期冲击载荷作用下巷道围岩的变形—开裂—运动过程进行了模拟.模型尺寸为40 m×40 m,被剖分成160×160个正方形单元,巷道尺寸为6 m×6 m.模型的左、右侧面均为透射边界,即应力波经过时不会发生反射.通过在巷道顶板布置多个测点,获取了最大主应力在整个模拟过程中的演化规律.该过程包括巷道开挖前模型的平衡过程、巷道开挖后模型的平衡过程及应力波在模型中传播过程.研究发现,当应力波在模型中传播时,弹射单元任一节点的水平速度呈现上升—稳定—衰减的变化过程,这与不进行刚化处理的结果相比更符合实际;巷道顶板左、右对称线上大部分测点的最大主应力呈现近似正弦波动上升—衰减—稳定的变化过程.在某一节点分离前,最大主应力首先呈现有规律的波动;然后,波动幅度突然增大,并伴随着剧烈震荡.阐明了顶板的拉裂机理:当应力波传入模型之后,顶板某一位置的最大主应力的波峰和波谷随着时间的增加呈增加的趋势,反射的多个拉应力波的累积作用不断提升该位置最大主应力的波峰,致使拉裂.另外,还初步分析了周期冲击载荷频率对单元弹射初速度和剪、拉裂缝区段数目的影响.分析发现,随着频率的增加,脱离巷道围岩的单元的平均弹射初速度增加,剪、拉裂缝区段数目减小.