DOI： 10.11817/j.issn.1672-7207.2020.03.023

基于库仑摩擦原理的充填体强度力学模型

常宝孟^{1, 2},杜翠凤^{1, 2},魏丁一^{1, 2},张爱卿³

（1. 北京科技大学 金属矿山高效开采与安全教育部重点实验室，北京，100083；

2. 北京科技大学 土木与资源工程学院，北京，100083；

3. 北华航天工业学院 建筑工程学院，河北 廊坊，065000）

摘 要：

单元为研究对象，基于库仑摩擦原理建立摩擦充填体强度力学模型，采用理论分析与工程实例相结合的方法研究其变化规律，并应用模型预测值进行现场试验验证模型有效性。研究结果表明：摩擦充填体强度模型预测值随充填体高度的增加而逐渐增加，且充填体高度超过30 m后，由于静摩擦向动摩擦的转化，其增加幅度明显降低；强度预测值随充填体长度的增加变化不大，与宽度无关。摩擦充填体强度模型预测值随充填体尺寸的变化规律与高水自立强度模型的相似，但摩擦充填体强度模型更简易且应用范围更广泛。充填体揭露过程中仍保持稳定，未出现垮塌现象，证明该力学模型是可靠的，该模型满足了安全性和降低充填成本的双重目的。

关键词:库仑摩擦原理；充填体尺寸；充填体强度；力学模型

中图分类号:TD853.34          文献标志码:A

文章编号:1672-7207（2020）03-0777-06

Strength mechanics model of filling body based on Coulomb tribology principle

CHANG　Baomeng^{1, 2}, DU　Cuifeng^{1, 2}, WEI　Dingyi^{1, 2}, ZHANG　Aiqing³

(1. State Key Laboratory of High-Efficient Mining and Safety of Metal Mines， Ministry of Education，

University of Science& Technology Beijing, Beijing 100083， China；

2. School of Civil & Resources Engineering, University of Science& Technology Beijing, Beijing 100083, China；

3. College of Architecture Engineering, North China Institute of Aerospace Engineering, Langfang 065000, China)

Abstract: Based on Coulomb tribology principle and taking an infinitesimal unit as the research object, a strength mechanics model of friction filling body was established and the variation law was studied by theoretical analysis and engineering examples. The field test was carried out by applying the predicted values to verify the validity of the model. The results show that the predicted strength of strength mechanics model of friction filling body increases gradually with the height of filling body, and the increasing extent decreases significantly when the height of filling body is over 30 m because of the static friction transforming to dynamic friction. The predicted strength changes little with the length and has no correlation with the width of filling body. The predicted strength of friction filling body model varies with the increase of size of filling body similarly to that of the high water self-support strength model, but the strength mechanics model of friction filling body is simpler and has wider application range. The filling body remains stable during the process of uncovering in the field test, and no collapse occurs. It is proved that the mechanical model is reliable and satisfies the dual purpose of safety and reduces the filling cost.

Key words: Coulomb tribology principle; size of filling body; strength of filling body; mechanics model

充填体强度是影响矿山安全生产的重要因素。若充填体强度过低，则无法起到充填效果，引起采场顶板垮冒及地表塌陷等问题；若充填体强度过高，则充填成本会相应增加^[1-2]。充填体强度在满足安全的前提下，工程应用中必须要考虑其成本^[3] 。因此，对于充填体强度的研究非常必要。许多国内外专家对充填体强度开展了大量研究。基于对充填体与围岩相互作用关系的分析，目前已有一些常用的充填体强度计算模型，如Terzaghi强度模型、Thomas强度模型、蔡嗣经经验公式法及高水自立强度模型^[4-5]；王俊^[6]在考虑充填体与围岩的接触条件和结构尺寸的前提下，建立了空场嗣后胶结充填体强度模型；HUANG等^[7]基于薄板理论确定了充填体强度的计算方法；杨磊等^[8]利用极限平衡法得出充填体矿柱与岩壁接触时的强度模型；刘光生等^[9]基于Marston拱应力二维模型建立了充填体的三维拱应力解析模型；曹帅等^[10]基于弹性力学平面应变假设建立了空场嗣后充填体的力学模型并进行了理论求解；叶洪涛^[11]运用弹性力学半逆解法建立了充填体-矿柱的力学模型，揭示了不同尺寸充填体对矿柱水平和剪切应力的变化规律。但是，由于缺少对充填体与围岩间摩擦的考虑，目前常用的充填体强度力学模型在计算时导致强度相对偏大，且对充填体尺寸有相应的适用范围^[12]。因此，本文作者在分析前人对充填体强度理论研究的基础上，基于库仑摩擦原理，建立一个易于计算且适用范围广泛的充填体强度力学模型，并通过工程实例对其进行可行性验证。

1 常用充填体强度的力学模型

基于对充填体与围岩相互作用关系的分析，许多学者已建立一些常用的充填体强度计算模型。

1.1　Terzaghi强度模型

由于充填材料的力学结构特性与固结土相似，在分析充填体的应力分布时常用太沙基计算法^[13]：

(1)

式中：σ_ν为充填体设计强度，MPa；H为充填体高度，m；W为充填体的宽度，m；γ为充填体容重，N/m³；C为材料黏聚力，MPa；φ为材料内摩擦角，(°)；K为充填体侧压系数；；L为充填体长度，m。

1.2　Thomas强度模型

充填材料与矿壁围岩间的摩擦力导致充填体内成拱^[14]，其作用在充填体底部的垂直应力为

(2)

式中：ρ_b为充填料的散体堆积密度。

1.3　蔡嗣经经验公式

蔡嗣经半立方抛物线强度经验公式^[15]为

(3)

式中：B为经验系数，当H＜50 m时，B=600；当H＞100 m时，B=100。

1.4　高水自立强度模型

高水自立强度模型^[4-5]如式(4)所示：

(4)

式中：σ_α为充填体自立强度，MPa；ρ为充填体密度，g/cm³；α为充填材料坍塌角，(°)，α=45°+φ/2；φ₁为充填体与围岩的摩擦角，(°)；C₁为充填体与围岩的黏结强度，MPa。

上述模型中，Terzaghi强度模型、Thomas强度模型和蔡嗣经经验公式均未充分考虑充填体与围岩间的摩擦，高水自立强度模型考虑了充填体与矿壁间的关系，但是它们之间的黏结强度和内摩擦角并未明确取值，给公式应用带来不便，且受到充填体高度H与充填体宽度W的限制，只有当满足时，计算结果才有意义。

2 基于库仑摩擦原理的充填体强度力学模型

2个固体在法向力作用下被挤压在一起时的摩擦力表现为静摩擦和动摩擦2种性质^[16]，且动摩擦因数μ_k近似等于静摩擦因数μ_S。当静摩擦力达到最大值，即临界状态时，

(5)

式中：β为滑动面角度。

摩擦充填体强度力学模型的假设条件如下：1) 充填体内部属于各向同性；2) 只考虑滑动体与两侧围岩间的摩擦，不考虑滑动体与前后围岩间的摩擦；3) 充填体与围岩间的内摩擦角根据充填体的内摩擦角取值^[17]；4) 充填体与围岩间的摩擦因数与充填体内部的摩擦因数相同。

在高水自立强度模型的基础上，摩擦充填体强度力学模型基于库仑摩擦原理计算其充填体两边的摩擦力，取充填体中无限小水平截面单元进行分析，如图1所示。

根据力的平衡可得：

(6)

根据库仑定律，摩擦力为

(7)

将式(7)代入式(6)可得：

(8)

分离变量，积分可得：

(9)

式中：σ为作用于充填体内部应力；Κ_a为主动土压力系数，。

根据假设条件可知，达到临界状态时，摩擦充填体自立强度公式为

(10)

由式(10)可知，摩擦充填体自立强度主要与充填体的容重、高度、长度及围岩间摩擦因数有关。

图1　基于库仑摩擦原理的充填体强度力学分析

Fig. 1　Strength mechanics analysis of filling body based on Coulomb tribology principle

3 充填体强度力学模型计算结果对比

为验证摩擦充填体强度可以满足安全性的要求，对内蒙古某铅锌矿进行现场工业试验，该矿山矿房长×宽×高为35 m×15 m×50 m。

3.1　充填体高度对强度预测值的影响

当充填体长度为35 m，矿房宽度为15 m，充填体黏聚力为105 kPa，内摩擦角为33.2°时，5种充填体强度模型预测值随充填体高度变化的规律如图2所示。

图2　充填体高度对充填体强度的影响

Fig. 2　Influence of room height on strength of filling body

由图2可知：5个充填体强度模型的预测值均随充填体高度增加而增加。其中，蔡嗣经经验公式计算的强度最高；由于公式适用范围的限制，只有当充填体高度大于35 m后高水自立强度模型才有预测值。而摩擦充填体强度模型预测值介于Thomas强度模型和高水自立强度模型的预测值之间，且充填体高度超过30 m后，充填体强度预测值增加幅度明显降低，这可能是充填体与围岩接触面的静摩擦转化为动摩擦所致。

3.2　充填体宽度对强度预测值的影响

当充填体长度为35 m，矿房高度为50 m时，5种充填体强度模型预测值随充填体宽度的变化规律如图3所示。

图3　充填体宽度对充填体强度的影响

Fig. 3　Influence of width on strength of filling body

由图3可知：Terzaghi强度模型和Thomas强度模型的强度预测值随充填体宽度增加而增大；高水自立强度模型预测值随充填体宽度增加而略有下降，但变化幅度不大；蔡嗣经经验公式和摩擦充填体强度模型的预测值变化与充填体宽度变化无关；蔡嗣经经验公式预测值最大。

3.3　充填体长度对强度预测值的影响

当充填体宽度为15 m，矿房高度为50 m时，5种充填体强度模型预测值随充填体长度的变化规律如图4所示。

图4　充填体长度对充填体强度的影响

Fig. 4　Influence of length on strength of filling body

由图4可知：Terzaghi强度模型、Thomas强度模型和蔡嗣经经验公式计算的充填体强度变化与矿房长度变化无关，高水自立强度模型预测值随充填体长度增加而增加。摩擦充填体强度模型预测值介于Thomas强度模型和高水自立强度模型的预测值之间，充填体长度的变化对摩擦充填体强度影响较小。

综上所述，摩擦充填体强度模型预测值随充填体尺寸的变化规律与高水自立强度模型的变化规律基本一致，因为二者均考虑了充填体与围岩间的摩擦，但摩擦充填体强度计算结果略大于高水自立强度模型的计算结果，这是由于摩擦充填体强度并未考虑充填体内部摩擦力。与高水自立强度模型相比，摩擦充填体强度力学模型较简单，且不受充填体高度和宽度的限制。因此，在工程应用中摩擦充填体强度力学模型对充填体强度预测有重要意义。

4 工程实例验证

为了验证摩擦充填体强度力学模型的正确性，在内蒙古某铅锌矿进行现场工业试验。该矿产能为30万t/a，采矿方法为浅孔留矿法，充填方式为全尾砂结构流体胶结充填。目前主要开采920 m和860 m阶段，矿房和矿柱长×宽×高分别为35 m×15 m×50 m和35 m×9 m×50 m。

测定全尾砂的基本物化参数和充填料浆的流变参数等，同时进行配比后的强度和物理力学参数实验。现场选取充填配比为1:6，质量分数为76%的全尾砂胶结充填体，28 d强度及物理力学参数如表1所示。

表1　充填体的物理力学参数

Table 1　Physical and mechanical parameters of filling body

将全尾砂胶结充填体物理力学参数分别代入5个充填体强度力学模型中，其计算结果如表2所示。由表2可知：5种模型计算的充填体强度预测值从大到小依次为：蔡嗣经经验公式、Terzaghi强度模型、Thomas强度模型、摩擦充填体强度模型、高水自立强度模型，且摩擦充填体强度模型预测值介于Thomas强度模型和高水自立强度模型的预测值之间。现场实测充填体强度约为1.68 MPa，远大于摩擦充填体强度模型计算的理论强度。

为了验证摩擦充填体强度可以满足安全性的要求，在试验矿房采用充填配比为1:12、质量分数为74%和强度为255 kPa的充填体进行充填。充填后回采矿柱，发现充填体揭露过程中仍保持稳定，未出现垮塌现象，证明该强度可以保证充填体自立，摩擦充填体强度力学模型是可靠的^[18-19]。

不同充填配比的全尾砂胶结充填体材料用量对比如表3所示。由表3可知：充填配比为1:12的充填体与充填配比1:6的充填体相比，全尾砂和水的用量较多，水泥用量减少了97.17 kg/m³。当地水泥价格为18元/袋，充填体可节省的水泥费用约为35元/m³。因此，摩擦充填体强度力学模型既满足了矿房充填体安全性的要求，又达到了降低充填成本的目的。

表2　不同模型计算的充填体强度

Table 2　Strengths of filling body calculated by different models   MPa

表3　不同配比时充填体材料用量对比

Table 3　Comparison of filling material dosage with different ratios

5 结论

1) 在分析已有充填体强度力学模型的基础上，基于库仑摩擦原理计算充填体与围岩间的摩擦力，建立了摩擦充填体强度力学模型。

2) 摩擦充填体强度模型预测值随充填体高度的增加而逐渐增加，且充填体高度超过30 m后，由于静摩擦向动摩擦的转化，充填体强度预测值的增加幅度明显降低；充填体强度预测值随充填体长度增加变化不大，与充填体宽度无关。

3) 摩擦充填体强度模型随充填体尺寸的变化规律与高水自立强度模型的基本一致。由于摩擦充填体强度并未考虑内部摩擦力，其预测值略大于高水自立强度模型的预测值。但摩擦充填体强度力学模型较高水自立强度模型简易，且适用范围广。因此，摩擦充填体强度力学模型对工程中充填体强度预测有重要意义。

4) 摩擦充填体强度力学模型满足了安全性要求和降低充填成本的双重目的。

参考文献：

[1] 李夕兵, 彭定潇, 冯帆, 等. 基于中厚板理论的深部崩落转充填隔离矿柱稳定性分析[J]. 中国矿业大学学报, 2019, 48(3): 484-494.

LI Xibing, PENG Dingxiao, FENG Fan, et al. Stability analysis of horizontal insulating pillar in deep mining from caving to filling method on the basis of refined plate theory[J]. Journal of China University of Mining & Technology, 2019, 48(3): 484-494.

[2] 孙希奎, 赵庆民, 施现院. 条带残留煤柱膏体充填综采技术研究与应用[J]. 采矿与安全工程学报, 2017, 34(4): 650-654.

SUN Xikui, ZHAO Qingmin, SHI Xianyuan. Research and application on the technology of paste backfilling fully mechanized in residual strip pillar[J]. Journal of Mining & Safety Engineering, 2017, 34(4): 650-654.

[3] 刘浪, 辛杰, 张波, 等. 矿山功能性充填基础理论与应用探索[J]. 煤炭学报, 2018, 43(7): 1811-1820.

LIU Lang, XIN Jie, ZHANG Bo, et al. Basic theories and applied exploration of functional backfill in mines[J]. Journal of China Coal Society, 2018, 43(7): 1811-1820.

[4] 蔡嗣经. 矿山充填力学基础[M]. 北京: 冶金工业出版社, 1994.

CAI Sijing. Mine filling mechanics foundation[M]. Beijing: Metallurgical Industry Press, 2009.

[5] 刘建庄, 赵春景, 浑宝炬. 胶结充填体自立的受力分析与研究[J]. 有色矿冶, 2008, 24(4): 12-13, 16.

LIU Jianzhuang, ZHAO Chunjing, HUN Baoju. Analysis and research of cement filling body to self-support[J]. Non-Ferrous Mining and Metallurgy, 2008, 24(4): 12-13, 16.

[6] 王俊, 乔登攀, 李广涛, 等. 大空区尾砂充填体自重压密模型与应用[J]. 岩土力学, 2016, 37(S2): 403-409.

WANG Jun, QIAO Dengpan, LI Guangtao, et al. Dead weight compression model and application of tailings filling body in large goaf[J]. Rock and Soil Mechanics, 2016, 37(S2): 403-409.

[7] HUANG Yucheng. Study on artificial roof stability of downward drift mining with high-water solidifying backfill[C]//Proceedings of the 8th International Symposium on Mining with Backfill. Beijing: The Nonferrous Metals Society of China, 2004: 192-195.

[8] 杨磊, 邱景平, 孙晓刚, 等. 阶段嗣后胶结充填体矿柱强度模型研究与应用[J]. 中南大学学报(自然科学版), 2018, 49(9): 2316-2322.

YANG Lei, QIU Jingping, SUN Xiaogang, et al. Research and application on strength model of cemented backfill pillar for stage subsequent filling mining method[J]. Journal of Central South University(Science and Technology), 2018, 49(9): 2316-2322.

[9] 刘光生, 杨小聪, 郭利杰. 阶段空场嗣后充填体三维拱应力及强度需求模型[J]. 煤炭学报, 2019, 44(5): 1391-1403.

LIU Guangsheng, YANG Xiaocong, GUO Lijie. Models of three-dimensional arching stress and strength requirement for the backfill in open stoping with subsequent backfill mining[J]. Journal of China Coal Society, 2019, 44(5): 1391-1403.

[10] 曹帅, 杜翠凤, 谭玉叶, 等. 金属矿山阶段嗣后充填胶结充填体矿柱力学模型分析[J]. 岩土力学, 2015, 36(8): 2370-2376.

CAO Shuai, DU Cuifeng, TAN Yuye, et al. Mechanical model analysis of consolidated filling pillar using stage-delayed backfill in metal mines[J]. Rock and Soil Mechanics, 2015, 36(8): 2370-2376.

[11] 叶洪涛. 充填体-矿柱系统力学模型构建及破坏机制研究[D]. 唐山: 华北理工大学矿业工程学院, 2018: 20-50.

YE Hongtao. Study on mechanics model construction and fracture mechanism of filling body-pillar system[D]. Tangshan, China: North China University of Science and Technology. School of Mining Engineering, 2018: 20-50.

[12] 黄庆享, 李亮. 充填材料及其强度研究[J]. 煤矿开采, 2011, 16(3): 38-42.

HUANG Qingxiang, LI Liang. Research on stowing material and its strength[J]. Coal Mining Technology, 2011, 16(3): 38-42.

[13] 蔺港, 孔令刚, 詹良通, 等. 基于太沙基土拱效应考虑基质吸力影响的松动土压力计算模型[J]. 岩土力学, 2015, 36(7): 2095-2104.

LIN Gang, KONG Linggang, ZHAN Liangtong, et al. An analytical model for loosening earth pressure considering matric suction based on Terzaghi soil arch effect[J]. Rock and Soil Mechanics, 2015, 36(7): 2095-2104.

[14] 杜晓丽, 宋宏伟, 陈杰. 煤矿采矿围岩压力拱的演化特征数值模拟研究[J]. 中国矿业大学学报, 2011, 40(6): 863-867.

DU Xiaoli, SONG Hongwei, CHEN Jie. Numerical simulation of the evolution of the pressure arch during coal mining[J]. Journal of China University of Mining & Technology, 2011, 40(6): 863-867.

[15] 吴爱祥, 沈慧明, 姜立春, 等. 窄长型充填体的拱架效应及其对目标强度的影响[J]. 中国有色金属学报, 2016, 26(3): 648-654.

WU Aixiang, SHEN Huiming, JIANG Lichun, et al. Arching effect of long-narrow cemented paste backfill body and its effect on target strength[J]. The Chinese Journal of Nonferrous Metals, 2016, 26(3): 648-654.

[16] 瓦伦丁 L 波波夫. 接触力学与摩擦学的原理及其应用[M]. 李强, 雒建斌，译. 北京: 清华大学出版社, 2019: 25-45.

POPOV V L, Contact mechanics and tribology principle and its application[M]. LI qiang, LUO jianbin, trans. Beijing: Tsinghua University Press, 2019: 25-45.

[17] 朱志彬, 刘成平. 充填体强度计算及稳定性分析[J]. 采矿技术, 2008, 8(3): 15-17, 25.

ZHU Zhibin, LIU Chengping. Filling body strength calculation and stability analysis[J]. Mining Technology, 2008, 8(3): 15-17, 25.

[18] ZHOU Baojing, XU Jinhai, ZHAO Maosen, et al. Stability study on naturally filling body in gob-side entry retaining[J]. International Journal of Mining Science and Technology, 2012, 22(3): 423-427.

[19] 韩斌, 吴爱祥, 王贻明, 等. 低强度粗骨料超细全尾砂自流胶结充填配合比优化及应用[J]. 中南大学学报(自然科学版), 2012, 43(6): 2357-2362.

HAN Bin, WU Aixiang, WANG Yiming, et al. Optimization and application of cemented hydraulic fill(CHF) with low strength aggregate and extra fine grain full tailings[J]. Journal of Central South University(Science and Technology), 2012, 43(6): 2357-2362.

（编辑 赵俊）

收稿日期： 2019 -07 -12; 修回日期： 2019 -08 -23

基金项目(Foundation item)：国家重点研发计划项目(2018YFC0604605)；国家自然科学基金资助项目(51274023) (Project(2018YFC0604605) supported by the Key Research and Development Program of China; Project(51274023) supported by the National Natural Science Foundation of China)

通信作者：杜翠凤，博士，教授，从事矿井通风、粉尘抑制及采矿方法研究；E-mail：ducuifeng@126.com

摘要:以充填体无限小单元为研究对象，基于库仑摩擦原理建立摩擦充填体强度力学模型，采用理论分析与工程实例相结合的方法研究其变化规律，并应用模型预测值进行现场试验验证模型有效性。研究结果表明：摩擦充填体强度模型预测值随充填体高度的增加而逐渐增加，且充填体高度超过30 m后，由于静摩擦向动摩擦的转化，其增加幅度明显降低；强度预测值随充填体长度的增加变化不大，与宽度无关。摩擦充填体强度模型预测值随充填体尺寸的变化规律与高水自立强度模型的相似，但摩擦充填体强度模型更简易且应用范围更广泛。充填体揭露过程中仍保持稳定，未出现垮塌现象，证明该力学模型是可靠的，该模型满足了安全性和降低充填成本的双重目的。