地表顶柱稳定性的砌体梁分析(Voussoir Analysis)
1 引言
地表顶柱的稳定性可以使用刚性分析【地表顶柱稳定性的刚性分析(Rigid Analysis)】和弹性分析【地表顶柱稳定性的弹性分析(Elastic Analysis)】，当侧向应力水平很低时并且跨度/高度比比较大时，顶柱在重力作用下可能会发生弯曲裂缝，一种非常经典的解析方法是砌体梁分析(Voussoir Analysis)。
2 破坏模式
砌体梁的破坏模式有三种：第一种是在柱墩出垂直滑动引起的剪切破坏，其计算原理与前两种方法相同；第二种是由于重力作用引起的弯曲破坏，这种破坏模式与弹性分析的弯曲类似，但梁会发生断裂；第三种是由于压缩产生的局部粉碎性破坏。
砌体梁分析作了如下假设：
(1) 顶柱不能承受拉应力，当发生弯曲时，主要通过形成抛物线压缩拱来实现自支撑；
(2) 砌体梁的变形发生在柱墩的滑动之前，抗滑稳定性在压力拱形成之后决定;
(3) 不考虑由原岩应力和/或开挖的几何形状引起的初始侧向应力，假定梁最初无应力；
(4) 横切结构的角度小于节理面假设的最小摩擦角。
3 输入参数
当使用Mohr-Coulomb强度准则进行刚性或弹性分析时，不需要原岩的单轴抗压强度UCS，但进行Voussoir分析，由于需要使用UCS计算局部挤压破坏的安全系数，因此无论采用哪一种强度准则，都必须输入UCS值。此外，如果使用Hoek-Brown准则，则需要输入岩石的弹性模量和泊松比。
4 计算指标
砌体梁分析有三个计算指标：
(1) 剪切安全系数
剪切安全系数的计算方法如下式所示:
其中
(2) 压缩安全系数
压缩安全系数FOS=UCS/产生的最大水平压应力
(3) 拱穿透概率
拱穿透概率(Arch Snap-Thru Probability)是砌体的主要破坏模式，计算使用的是Diederichs(1999)建议的算法[Diederichs, M. S. and P. K. Kaiser (1999). "Tensile strength and abutment relaxation as failure control mechanisms in underground excavations." 。
在Voussoir分析中计算的屈曲参数(buckling parameter)代表了由Voussoir分析迭代程序确定的在给定的几何形状和岩体模量下不稳定的拱形结构的百分比。Hutchinson等人(1996)确定出屈曲参数=35%作为超过屋顶应被视为不稳定的极限。这个35%的设计极限恰好对应于跨中挠度的变形梁，即为梁厚度的10%。因此，拱的稳定性也可以通过监测中跨的位移来评估，相对于无弯曲的状态。上述输入数据的计算结果如下：
可以看出，在剪切作用下破坏的失败的概率为0.6%，压缩作用下破坏的概率为0%，拱穿透下的破坏概率很低。如果逐渐减小矿柱的厚度，破坏概率逐渐增加，当矿柱厚度达到0.5m时，三种破坏模式同时达到，致使顶柱塌落。
5 结束语
随着数值模拟技术的发展，传统的砌体梁理论受到了严重挑战，特别是基于离散块体理论的技术能够更深入地了解顶柱的破坏机理，Hatzor(1998)[The stability of a laminated Voussoir beam: Back analysis of a historic roof collapse using DDA]使用DDA方法模拟了砌体梁；Oliveira(2014)[Revisiting the applicability of voussoir beam theory for tunnel design in Sydney]使用UDEC模拟了砌体梁；Abousleiman(2020)[Expanding Understanding of the Voussoir Beam Analog in Flat Roof Excavations Using the Discrete Element Method]也使用离散元模拟了砌体梁，更详细的讨论参看"voussoir beam theory"数据集。
