罚单元与有限元耦合法计算巷道周边应力分布

将边界罚单元引入有限元方程来处理复杂边界条件问题，开巷初期在巷道周边设置一周边界弹簧元。通过２维元与边界罚单元（弹簧元）耦合分析，可计算出开巷初期巷道周边应力的分布。该方法简单易行，计算结果与实测位移接近，具有较高的精度，实现了巷道围岩变形过程的动态模拟。