岩溶地层长锚索穿越溶洞锚固力计算方法研究_肖国庆.pdf

第 20 卷 第 2 期2023 年 2 月铁道科学与工程学报Journal of Railway Science and EngineeringVolume 20 Number 2February 2023岩溶地层长锚索穿越溶洞锚固力计算方法研究肖国庆，施成华，叶承敏，孙晓贺，宋银涛，郑可跃(中南大学 土木工程学院，湖南 长沙 410075)摘要：目前桩锚围护结构是基坑工程常用的支护方式，在岩溶地层进行施工时常会遇到长锚索穿越溶洞的情况。溶洞填充物的力学性质与一般地层相差较大，长锚索锚固力会有一定程度降低，进而对地铁基坑的稳定性产生不利影响。以徐州地铁2号线淮塔东路站基坑工程为背景，针对岩溶地层长锚索穿越溶洞时锚固力的计算问题，通过对比分析岩溶地层与一般地层长锚索的受力特征，考虑锚索穿越溶洞时剪应力的折减，建立长锚索穿越溶洞锚固力的计算模型。进而基于剪切位移法原理，采用荷载传递法，推导得到相应的长锚索锚固力计算方法。该计算方法得到的理论值与实测值吻合较好，从而验证了计算方法的可靠性。分析溶洞直径、锚固始端距溶洞距离及溶洞填充物剪切模量对长锚索锚固力的影响规律。研究结果表明：当溶洞直径不变时，溶洞距锚固始端距离越远，锚索锚固力越大，当溶洞分布在距锚固始端较近的位置(4 m内)时，锚固力的值受溶洞的影响很大，在此之外的溶洞影响较小；当溶洞与锚固始端距离不变时，溶洞直径越大，锚索锚固力越小；当溶洞直径及溶洞与锚固始端距离均不变时，溶洞填充物剪切模量越大，锚索锚固力也越大。当溶洞填充物的剪切模量较小时，采用注浆方式填充溶洞可有效提高长锚索的锚固力。关键词：岩溶地层；荷载传递法；剪切位移法；锚固力中图分类号：TD353 文献标志码：A 开放科学(资源服务)标识码(OSID)文章编号：1672-7029（2023）02-0641-10Study on calculation method of anchorage force of long anchor cable crossing karst cave in Karst stratumXIAO Guoqing,SHI Chenghua,YE Chengmin,SUN Xiaohe,SONG Yintao,ZHENG Keyue(School of Civil Engineering,Central South University,Changsha 410075,China)Abstract:Pile-anchor retaining structure is a common supporting method in foundation pit engineering,and long anchor cables often pass through karst caves during construction in karst strata.The mechanical properties of the cave filling are much worse than those of the general layer,so the anchoring force of the long anchor cable will be reduced,which will adversely affect the stability of the subway foundation pit.The foundation pit project of Xuzhou Metro Line 2 Huaita East Road Station was taken as the engineering background,and the calculation problem of the anchorage force of long anchor cable crossing karst cave in karst layer was considered.This paper established a calculation model of the anchorage force of long anchor cable crossing karst cave by comparing the 收稿日期：2022-02-28基金项目：国家自然科学基金资助项目(52178402)通信作者：施成华(1973)，男，安徽黄山人，教授，博士，从事隧道与地下工程方面的教学与科研工作；Email：DOI:10.19713/ki.43-1423/u.T20220352铁 道 科 学 与 工 程 学 报2023 年 2月stress characteristics of long anchor cable in karst layer and general layer and considering the reduction of shear stress in karst layer.Based on the principle of shear displacement method and load transfer method,the calculation method of anchorage force of long anchor cable was derived.The theoretical values obtained by this calculation method were in good agreement with the measured values,which verified the reliability of the calculation method.Finally,the influences of the diameter of the cave,the distance from the anchorage start,and the shear modulus of the cave filling on the anchorage force of the long anchor cable were analyzed.The results show that when the diameter of the cave is constant,the farther the cave is from the anchorage beginning of the anchorage,the greater the anchorage force is.When the cave is close to the beginning of the anchorage(within 4 m),the anchorage force is significantly influenced by the cave and the cave outside is less influenced.When the distance between the cave and the anchorage end is constant,the larger the cave is,the smaller the anchorage force is.When the diameter of the cave and the distance between the cave and the anchorage end remain unchanged,the long anchor cable gradually increases as the shear modules of the cave filling increases.When the shear modulus of the filler is not sufficient,filling the cavity with grout can effectively improve the anchoring force of the anchor cable.Key words:Karst stratum;load transfer method;shear displacement method;anchorage force 随着城市发展的需要，在岩溶分布较广地区建设地铁已不可避免。由于岩溶地层具有的复杂性，对地铁建设将产生很大的困难与挑战1。目前车站施工常采用明挖深基坑法，在基坑工程施工过程中存在一些难以预估的风险，如：地表沉降24、基坑塌陷58等。而这些风险归根结底还是与围护结构的稳定有关，当基坑围护结构采用桩锚结构时，锚索穿过溶洞时锚固力有一定程度的降低，影响基坑施工安全。现阶段关于锚索锚固力计算理论的研究，主要有2个方向。一是假设注浆体与围岩界面间剪应力沿锚固段全长均匀分布的平均剪应力理论。LITTLEJOHN等9假设黏结力沿锚固段全长均匀分布，提出了拉力型锚索锚固力设计公式。程良奎等10假定锚固体与周边岩土体界面的黏结强度为定值，提出了压力型锚索的设计方法。二是假设锚固力沿锚索轴线非均匀分布的非平均剪应力理论。PHILLIPS11认为拉力型锚索的锚固体剪应力沿锚固段呈指数分布，并由此得到了极限锚固力计算公式。LI等12提出了“锚固解耦段”的定义，认为部分解耦段与未解耦段面积组合的最大值即为注浆锚杆的极限锚固力。KIM13为得到锚固体与岩体分界面的荷载传递规律，采用t-z曲线法，取得了较好的研究效果。刘小斌等14通过对不同长度拉力型锚杆进行循环加载极限试验，发现锚固体轴向应力先升后降，黏结应力具单峰性，以及拉力型锚杆的锚固长度与极限抗拔承载力成正比。LIU等15为研究锚杆荷载传递规律及不同锚固条件下的有效性，考虑节理张开位移，将弹性地基梁模型和锚杆拉拔模型相结合，推导出单个节理面沿锚杆应力分布的解析解。HYETT 等16基于剪应力与剪切位移的线性关系，分析了裂隙围岩下锚杆的轴力分布规律，适用于实际的裂隙围岩条件。周炳生等17基于锚固体界面剪切位移分布曲线的统一形式，推导了无限长锚体界面的位移、轴力以及剪切应力的计算公式。CHEN等18分析了锚杆的系数、弹性模量和直径对锚杆荷载传递性能的影响，认为锚杆的轴向载荷从锚固端向自由端衰减，与拉拔荷载无关。刘超等19考虑压力型锚杆灌浆体受压产生径向膨胀的影响，推导了锚杆极限抗拔承载力与临界锚固长度的计算公式。综上所述，关于一般地层锚索锚固力的计算，国内外学者已进行了较为深入的研究，然而针对岩溶地层长锚索穿越溶洞时锚固力的计算，目前尚未有明确的计算方法，同时锚索锚固力受溶洞的影响因素也缺乏系统的认识。本文依托徐州地铁2号线淮塔东路站基坑支护工程，采用荷载传递法，基于剪切位移法非平均剪应力理论，推导锚索穿过溶洞的锚固力理论计算公式，同时642第 2 期肖国庆，等：岩溶地层长锚索穿越溶洞锚固力计算方法研究对锚固力的影响因素进行系统的分析，研究结果可为工程实际提供理论指导。1 长锚索穿越溶洞锚固作用机理1.1一般地层锚索锚固作用机理锚固单元体，主要由锚固剂、锚杆(索)、锚固剂-围岩界面、锚固剂-杆体界面、围岩组合而成，简称“三体两面”结构20。如图1所示，锚索由锚头、自由段和锚固段3部分组成。将锚固剂与锚索之间接触面称为第一界面，围岩与锚固剂接触面简称为第二界面。当长锚索发挥支护作用时，锚固单元体的应力状态如图2所示，其受到的径向和轴向应力分别以v和r表示。基于局部共同变形假设，把围岩体对锚索的剪切力简化为一系列独立作用的切向弹簧，同时锚固体微段上的应力与锚固体和围岩之间的相对位移是线性关系。对一般地层，锚索没有穿越溶洞，剪应力应处处连续，根据该假定的锚索锚固段力学计算简图如图3所示。1.2考虑岩溶锚索锚固作用机理以全填充型溶洞为例，当锚索锚固段部分穿过溶洞时，此时锚索锚固力有一定程度的降低，计算简图应按图4所示。此时把锚固段分成3段分图1锚固单元体结构组成示意图Fig.1Schematic diagram of the structural composition of the anchoring unit body图2锚固单元体应力分析Fig.2Stress analysis of anchorage unit图3不考虑岩溶锚固段力学计算简图Fig.3Mechanical calculation diagram without considering karst anchorage section图4考虑岩溶锚固段力学计算简图Fig.4Mechanical calculation of karst anchorage section643铁 道 科 学 与 工 程 学 报2023 年 2月开计算，值得注意的是，在岩层与溶洞交界处，应满足变形协调条件与静力平衡条件，即剪应力会发生突变，而正应力与位移连续。2 长锚索穿越溶洞锚固力计算方法2.1计算假定考虑到实际施工时锚索打入岩层时有一定角度，将锚索锚固体划分为若干子锚固段，如图5所示，并作如下假设21：1)同一子锚固段锚固体及周围岩土体均为均质弹性介质；2)同一子锚固段锚固体和岩土体之间力学参数为常值；3)锚固体横截面上正应力均匀分布，不考虑注浆体的破坏；4)锚固体与周围岩土体间无明显的相对滑移；5)溶洞填充物与周围岩层结合密切，与锚索锚固体也紧密接触。2.2锚索锚固力的推导根据荷载传递法，将锚固体视为连续弹性体，分成若干连续弹性的单元22，如图6所示，根据受力平衡，对z位置的锚固段dz有：(z)=-1CdP()zdz(1)而dz段的压缩量为：ds(z)=-P()zEAdz(2)联立式(1)和(2)可得：EAds2()zdz2-(z)C=0(3)式中：C为锚固体的周长；z为距锚固始端距离。设锚固体的半径为r0，弹性模量为Ep以及第i段锚固体与周围岩体剪切面上的剪切位移和剪应力分别为ui(r0)和i(z)。由式(3)得：d2ui()r0dz2-2r0Epi(z)=0(4)基于剪切位移法原理23，锚固体侧岩土体剪应力与剪切位移为线性关系。锚索锚固体侧岩土体剪应力和剪切位移可表示为：ui(r)=|i()z r0Gisln()rimrr0rrim0rrim(5)式中：ui(r)为距锚固体轴线水平距离r处的位移；Gis为周边第i段岩体的剪切模量；rim为第i段岩层中 变 形 可 以 忽 略 不 计 的 最 大 半 径，一 般rim=2.5li(1-is)；li为第i段连续岩层长度；is为第i段岩层的泊松比。联立式(4)和式(5)可得：d2ui()r0dz2-kiApEpui(r0)=0(6)式中：ki为第i段岩层的剪切刚度。ki=2Gisln()rimr0(7)求解式(6)可得：ui(r0)=C1eiz+C2e-iz(8)式中：C1和C2为待定常数；i为与第i段岩层剪切图5子锚固段分析模型Fig.5Analysis model of sub-anchorage section图6荷载传递法计算简图Fig.6Calculation diagram of load transfer method644第 2 期肖国庆，等：岩溶地层长锚索穿越溶洞锚固力计算方法研究刚度有关的常数。i=kiApEp(9)根据锚索锚固体轴力Pi(z)与应力应变之间关系可得：Pi(z)=-r20Epdui()r0dz(10)将式(8)代入式(10)和式(5)，可得连续锚固体的锚索轴力与剪应力的分布函数为：Pi(z)=-iApEp(Ci1eiz-Ci2e-iz)(11)i(z)=Gis()Ci1eiz+Ci2e-izr0ln()rimr0(12)式中：Ap为锚索锚固体的截面积；Ci1和Ci2为第i段微分待定常数。考虑计算方便，取计算图示为全填充型溶洞，如图7所示，图中1和3中2段为岩层，2段为全填充型溶洞，锚索在A和B两截面为岩层与溶洞的分界，岩层与溶洞填充物仅有剪切刚度的不同。设锚索1，2，3段剪切刚度分别为k1，k2，k3，则根据式(9)得：1=k1ApEp;2=k2ApEp;3=k3ApEp当z=0时，由式(11)，第1段锚固体有：P1(0)=-1ApEp(C11-C12)=P(13)当z=x时，对第1，2段锚固体，剪切位移与轴向力是相等的，由式(11)和式(12)：C11e1x+C12e-1x=C21e2x+C22e-2x(14)-1ApEpC11e1x-C12e-1x=-2ApEpC21e2x-C22e-2x(15)同理当z=x+s时，第2，3段锚固体有：C21e2()x+s+C22e-2()x+s=C31e3()x+s+C32e-3()x+s(16)-2ApEpC21e2()x+s-C22e-2()x+s=-3ApEpC31e3()x+s-C32e-3()x+s(17)当z=L时，由式(11)，第3段锚固体有：-3ApEp(C31e3L-C32e-3L)=0(18)当锚固体周边剪应力超过其弹性极限黏结强度时，界面将发生滑动破坏。故计算锚索穿越溶洞锚固力时，以锚固始端剪应力达到弹性极限黏结强度时锚固体的轴力值作为锚索锚固力值。根据式(12)，有：G1s()C11+C12r0ln()r1mr0=1(19)式中：1为第1段锚固体弹性极限黏结强度；G1s为第1段锚固体周边岩体的剪切模量。联立以上(13)(19)各式，可得最终锚索穿越全填充型溶洞锚固力P为：P=1-h1+h1ApEp1r0ln()rmr0G1s(20)其中：|h=1+f()1-f e21xf=2()ne22x-1ne22x+1n=m+2()2-m e22()x+sm=3e23()x+s-L-1e23()x+s-L+1(21)分析式(20)及(21)可知，长锚索穿越溶洞锚固力的计算相较一般地层增加溶洞与锚索的位置关系、溶洞直径以及溶洞填充物的力学参数。因此，确定这几项参数后便可以运用公式(20)进行长锚索穿越溶洞时锚固力的计算，定量分析在特定的岩溶条件下锚索锚固力是否能达到设计值。以上分析建立了锚索穿过全填充型溶洞时锚固力的计算方法，对于锚索穿过无填充型溶洞的情况，同样如图7所示，只需令溶洞段k2,2为0，图7全填充型溶洞锚固体剪应力分布图Fig.7Shear stress distribution diagram of anchor solid in fully filled karst cave645铁 道 科 学 与 工 程 学 报2023 年 2月仍按照全填充型溶洞的计算方法进行计算。当锚索穿过半填充性溶洞时，可将溶洞段进行细分，分为若干子锚固段，这样就可将溶洞段分为若干不同小段的全填充型和无填充性溶洞进行分别计算。限于篇幅原因，本文不过多赘述。上述公式是假设锚固体及其周围岩土体为弹性介质、且无相对滑移推导得到，因此式(20)及(21)适用于易出现溶洞的均质岩土体，如石灰岩、白云岩及泥灰岩等。2.3计算方法的验证2.3.1工程概况选取徐州地铁 2 号线淮塔东路站基坑里程K11+400.000断面进行验证。下伏基岩为中风化灰岩，岩层下含有直径为1 m的全填充型溶洞，充填物为硬可塑状褐黄黏土2425，夹灰岩、泥灰岩碎块。锚固体半径 r0为 0.075 m，弹性模量为 3.02104 MPa，锚固体弹性极限黏结强度 1为 150 kPa。锚固段岩层剪切模量G1s=G3s=45 MPa，泊松比1s=3s=0.3，溶洞填充物的泊松比2 s为0.36，剪切模量G2 s为 4 MPa。锚索选取第 2 道锚索，自由段 Lf=5 m，锚固段长度La=9 m，其中第1段锚固体长度为6 m。锚索选用3s15.2钢绞线，钢绞线标准强度1 860 MPa，锚索轴向拉力设计值180 kN。2.3.2实测与计算结果对比对第2道穿越溶洞的锚索进行抗拔试验，溶洞与锚索位置关系如图8所示。试验中采用千斤顶进行渐进逐级加载，如图 9 所示。本试验共进行 3次，当锚索锚固力趋于稳定时，读取得到锚索锚固力的值分别为 185.67，185.89 和 184.82 kN，平均值为185.46 kN。采用本研究所建立的计算方法，将现场相关地层与锚固体参数代入计算公式，得到锚固力值为213.88 kN。锚索的特征荷载包括锚索工作荷载与极限荷载，工作荷载指锚索正常工作时所承受的拉拔荷载，而极限荷载指锚索所能承受的最大荷载26。现场实测时，锚索未达极限值，尚属工作荷载。理论计算时，以锚固始端黏结强度达到其弹性极限黏结强度，锚固体与岩体即将发生相对滑移得到，理论计算结果为锚索弹性极限锚固力，属极限荷载。根据建筑基坑支护技术规程(JGJ1202012)27，锚杆(索)极限抗拔承载力与其轴向拉力标准值之比，即为抗拔安全系数，不应小于1.4。江冠等26根据现场6根锚索实测，得到极限荷载与工作荷载比值范围为1.11.3。当理论计算极限荷载213.88 kN时，根据规范计算理论值对应的工作荷载为152.77 kN，由江冠的研究，可得对应的工作荷载164.53194.44 kN。而现场实测值为 185.46 kN，与规范相比偏安全，而与江冠等现场实测数据吻合较好，因此本研究所提计算方法可为岩溶地层地铁车站长锚索穿越图8第2道锚索与溶洞位置关系Fig.8Position relationship between the second anchor cable and the karst cave(a)现场实况图；(b)工作原理图图9锚索拉拔试验示意图Fig.9Schematic diagram of anchor cable pulling experiment646第 2 期肖国庆，等：岩溶地层长锚索穿越溶洞锚固力计算方法研究溶洞锚固力的计算提供参考。3 锚索穿越溶洞锚固力影响因素分析为研究锚索锚固力受岩溶条件的影响程度，由以上锚索锚固力计算公式可知，共有 1，2，3，x及s 5个参数，由于1和3为岩层段与剪切刚度有关的量，一般没有变化，对最终结果影响不大，而2与溶洞的填充物剪切模量相关，故考虑剪切模量G2s，即实际影响因素为锚固始端距溶洞的距离 x，溶洞直径 s 及溶洞填充物的剪切模量G2s。其他参数同2.3.1节。3.1锚固始端距溶洞的距离的影响取溶洞直径 s分别为 1，2，3和 4 m，根据上述锚固力公式，绘图得到锚索锚固力随锚固始端距溶洞距离x变化的曲线如图10所示。单独分析每条曲线，即溶洞直径不变时，随着溶洞段距锚固始端的距离的增加，锚索锚固力的值逐渐增大，增长速率逐渐减小，最后趋于稳定。这与实际相符合，当x+时，即锚固体穿过岩层长度无限长，此时溶洞对锚固力影响不大。由图10易得，当溶洞直径s分别为1，2，3和4 m，锚固始端距溶洞距离x2.3，2.6，2.8和2.9 m时，锚索锚固力即可达到设计值180 kN。横向比较各条曲线，对于相同的锚固始端距溶洞距离，随着溶洞直径的增大，锚固力将有一定程度降低。当溶洞直径不超过 2 m，即s2m时，当溶洞距锚固始端的距离大于4 m，锚索的锚固力均不低于 200 kN；而s2 m时，锚索的锚固力不全大于200 kN。当溶洞段距锚固始端的距离为0.5 m时，锚索的锚固力均不超过110 kN。而当溶洞距锚固始端的距离大于4 m，锚索锚固力均不低于180 kN。故当溶洞分布在距锚固始端较近的位置(4 m 内)时，锚固力的值受溶洞的影响很大，在此之外的溶洞对锚索锚固力的影响较小。究其原因是锚索的剪应力值在锚固始端是最大的，剪应力主要分布为指数递减曲线，主要集中在始端附近，锚固始端较近的区域越长，剪应力分布范围越广，所得的锚固力值也就越大。3.2溶洞直径的影响取锚固始端距溶洞的距离x分别为1，2，3和4 m，得到锚索锚固力随锚固始端距溶洞距离x变化的曲线如图11所示。单独分析每条曲线，当溶洞段距锚固始端的距离不变时，随着溶洞直径增加，锚索锚固力的值逐渐减小。当溶洞段距锚固始端的距离分别为x=1 m，2 m，无论溶洞直径多大，锚索的锚固力均不能达到设计值180 kN的要求。而当溶洞段距锚固始端的距离x=3 m，溶洞直径只需s4.3 m，锚索锚固力即可以达到设计值180 kN的要求，这表明，对于特定的溶洞，可根据溶洞直径来确定最小锚固始端距溶洞的距离。而当溶洞距锚固始图10锚固力随锚固始端距溶洞的距离变化图Fig.10Anchoring force changes with the distance from anchorage beginning to the cave图11锚固力随溶洞直径变化图Fig.11Anchorage force changes with the diameter of karst cave647铁 道 科 学 与 工 程 学 报2023 年 2月端的距离x4 m，此时无论溶洞直径多大，锚固力变化不大。溶洞直径对锚固力影响较小，锚固力能满足设计180 kN要求。横向比较各条曲线，对于同一溶洞直径s，随着锚固始端距溶洞的距离x的增大，锚固力逐渐增大。且随着溶洞直径的增大，锚固力差值增加更加显著。因此，在实际工程中应保证锚固始端距溶洞距离尽量大，以确保基坑工程的稳定性。3.3溶洞填充物剪切模量的影响为分析溶洞填充物剪切模量G2s对锚索穿过溶洞时锚固力的影响，取溶洞段距锚固始端的距离x=2 m，溶洞直径 s=2 m，填充物剪切模量 G2s变化范围为145 MPa，其他参数与上述保持一致，得到锚索锚固力随溶洞填充物剪切模量G2s的变化曲线图，如图12所示。由图12可知：当溶洞直径和溶洞距锚固始端的距离不变时，随着溶洞填充物剪切模量的增大，锚索锚固力的值逐渐增大，大致趋势为线性增长。当溶洞填充物剪切模量大于25 MPa时，锚固力已经达到了180 kN的设计值。因此，当溶洞填充物剪切模量足够大时，锚索穿过溶洞时，不用对溶洞进行注浆填充，锚索锚固力也能达到设计要求。这与实际相符，当锚固力不满足设计要求时，常进行注浆加固，而注浆加固实质就是提高溶洞填充物的剪切模量。4 结论1)基于剪切位移法原理，采用荷载传递法，考虑锚固体穿过溶洞时，分界处的变形协调与力的平衡关系，建立岩溶地层车站基坑长锚索穿越溶洞时锚固力的计算模型，推导得到相应的长锚索锚固力计算方法，通过与实际现场测试的锚固力对比，理论计算方法与实际吻合较好，计算方法可为长锚索穿越溶洞锚固力的计算提供参考。2)当溶洞直径不变时，随着溶洞距锚索锚固始端距离的增大，长锚索的锚固力也逐渐增大，当溶洞距锚索锚固始端距离小于4 m内时，溶洞对锚索锚固力的影响很大，而在4 m外影响较小。3)当溶洞距锚索锚固始端距离不变时，随着溶洞直径增大，长锚索的锚固力逐渐减小，当溶洞直径较大且锚固始端距溶洞小于3 m时，必须对溶洞进行注浆加固，以使锚索的锚固力达到设计要求。4)当溶洞距锚索锚固始端距离和溶洞直径均不变时，随着溶洞填充物剪切模量的增大，长锚索的锚固力也逐渐增大。计算分析表明，在溶洞段距锚固始端的距离x=2 m，溶洞直径s=2 m该工况下，当溶洞填充物的剪切模量大于 25 MPa时，锚索锚固力达到设计值，可以直接穿过溶洞进行锚固施工。而当填充物的剪切模量较小时，采用注浆方式填充溶洞可有效提高锚索的锚固力。参考文献：1邹洋,彭立敏,张智勇,等.基于突变理论的岩溶隧道拱顶安全厚度分析与失稳预测J.铁道科学与工程学报,2021,18(10):26512659.ZOU Yang,PENG Limin,ZHANG Zhiyong,et al.Safety thickness analysis and stability prediction of tunnel roof in Karst region based on catastrophe theory.J.Journal of Railway Science and Engineering,2021,18(10):26512659.2郑映旭.岩溶地层中盾构隧道的影响特性研究D.广州:广州大学,2020.ZHENG Yingxu.Study on influence characteristics of shield tunnel in Karst stratumD.Guangzhou:图12锚固力随溶洞填充物剪切模量变化图Fig.12Anchoring force changes with shear modulus of cave filling648第 2 期肖国庆，等：岩溶地层长锚索穿越溶洞锚固力计算方法研究Guangzhou University,2020.3崔庆龙.岩溶地层盾构施工引起地表沉降的规律及预测研究J.现代隧道技术,2021,58(S2):110119.CUI Qinglong.Research on the law and prediction of ground settlement caused by shield construction in Karst strataJ.Modern Tunnelling Technology,2021,58(S2):110119.4WALTHAM A C,FOOKES P G.Engineering classification of Karst ground conditionsJ.Quarterly Journal of Engineering Geology and Hydrogeology,2003,36(2):101118.5堵紫薇.隐伏岩溶地层对半盖挖车站基坑稳定性影响及控制技术研究D.徐州:中国矿业大学,2019.DU Ziwei.Study on influence and control technology of concealed Karst stratum on stability of foundation pit of semi-covered excavation stationD.Xuzhou:China University of Mining and Technology,2019.6李结全.岩溶地区地铁隧道盾构掘进安全控制技术研究D.南宁:广西大学,2018.LI Jiequan.Research on safety control technology of metro shield tunnel construction in Karst areaD.Nanning:Guangxi University,2018.7张政.岩溶溶洞对盾构隧道掘进安全性的影响特征与处理方法研究D.济南:山东大学,2021.ZHANG Zheng.Study on the influence characteristics and treatment methods of Karst cave on the safety of shield tunnelingD.Jinan:Shandong University,2021.8李慎奎.武汉地区沙漏型岩溶塌陷数值分析与模型试验研究J.隧道建设(中英文),2020,40(7):981987.LI Shenkui.Numerical analysis and model test of hourglass-shaped Karst collapse in Wuhan areaJ.Tunnel Construction,2020,40(7):981987.9LITTLEJOHN G S,MITCHELL J M,JONES D L,et al.Half-day meeting.bs 8081 ground anchorages:a new British standard code of practiceJ.Proceedings of the Institution of Civil Engineers,1990,88(6):11171121.10 程良奎,刘启琛.岩土锚固工程技术的应用与发展:国际岩土锚固工程技术研讨会论文集(1996年:柳州)M.北京:万国学术出版社,1996:216221.CHENG Liangkui,LIU Qichen.Application and development of geotechnical anchorage engineering technology:proceedings of the international symposium on geotechnical anchorage engineering technology(Liuzhou,1996)M.Beijing:International Academic Press,1996:216221.11 PHILLIPS S H E,Factors affecting the design of anchorages in rockR.London:Cementation Research Ltd.,1970.12 LI C,STILLBORG B.Analytical models for rock boltsJ.International Journal of Rock Mechanics and Mining Sciences,1999,36(8):10131029.13 KIM N K.Performance of tension and compression anchors in weathered soilJ.Journal of Geotechnical and Geoenvironmental Engineering,2003,129(12):1138 1150.14 刘小斌,李兆锋,泰培,等.拉力型锚杆锚固体应力分布规律的现场试验研究J.地下空间与工程学报,2021,17(S1):6370.LIU Xiaobin,LI Zhaofeng,TAI Pei,et al.In-situ experimental investigation on stress distribution of grout body of tension-type ground anchorJ.Chinese Journal of Underground Space and Engineering,2021,17(S1):6370.15 LIU Junfeng,YANG Haiqing,WEN Haijia,et al.Analytical model for the load transmission law of rock bolt subjected to open and sliding joint displacementsJ.International Journal of Rock Mechanics and Mining Sciences,2017,100:19.16 HYETT A J,MOOSAVI M,BAWDEN W F.Load distribution along fully grouted bolts,with emphasis on cable bolt reinforcementJ.International Journal for Numerical and Analytical Methods in Geomechanics,1996,20(7):517544.17 周炳生,王保田,梁传扬,等.全长黏结式锚杆锚固段荷载传递特性研究J.岩石力学与工程学报,2017,36(S2):37743780.ZHOU Bingsheng,WANG Baotian,LIANG Chuanyang,et al.Study on load transfer characteristics of wholly grouted boltJ.Chinese Journal of Rock Mechanics and Engineering,2017,36(S2):37743780.18 CHEN Jianhang,HE Fulian,