预应力锚杆作用机制的轴对称理论分析.docxVIP

预应力锚杆作用机制的轴对称理论分析 1 有限元计算的不足 预测生根设计的关键是确定固体衬底和稳定岩体之间的侧阻力分布模式。在现有的相关设计规范(《建筑边坡工程技术规范》(GB50330–2002)和《铁路路基支挡结构设计规范》(TB10025–2001))中，都假设锚固段侧阻力分布模式为均匀分布，而大量的研究结果表明：锚固段侧阻力并不是均匀分布的，而是在其前端形成峰值，逐步向末端减少并最终趋近于0，按照均匀分布模式计算是不合理的，有时甚至是不安全的。 有限元计算及模型试验结果表明：锚束体与浆体材料之间的结合应力在中等坚硬及坚硬地层中呈指数分布，其分布形态取决于荷载大小、锚束体弹性模量与地层、浆体材料弹性模量间的比例，比值越小，锚固段外端应力越集中，比值越大，侧阻力分布越均匀。 Hawkers,Evans(1951)以及后来的Phillips给出了一个钢材与浆体结合应力的计算公式： 式中：τx为距离锚固端外端x处的结合应力，τ0为锚固端外端处的结合应力，d为锚束体直径，A为锚束体与浆体间的结合应力与主应力相关的常数。 依此为基础，进而假定锚索锚固段浆体与岩体间也存在类似的关系式。但该计算过于简单，不能反映锚束体与围岩体相对刚度对侧阻力变化规律的影响。I.W.Farmer和A.Holmberg给出的锚杆锚固段侧阻力计算公式也为指数形式，即 其中， 式中：νr为岩体的泊松比，νg为灌浆材料的泊松比，dg为钻孔直径，d0为岩体影响直径。 式(2)应用较为普遍，特别在国外的相关文献中应用较多。本文在考虑锚束体与围岩体界面本构模型的基础上，将界面分为完全黏结和脱黏两种情况对界面的剪应力分布规律进行分析，给出了一个计算锚杆剪应力、轴力分布的理论计算模型。 2 锚束体-剪切位移间的线性动态本构模型 Y.Cai等采用2个阶段线性函数来描述预应力锚索锚固段岩体与灌浆材料之间接触面上的剪应力–剪切位移关系。第1阶段对应于弹性阶段，锚束体与围岩体之间完全耦合，接触面上剪应力与剪切位移呈比例变化，在此阶段，接触面处于无损状态；第2阶段对应于接触面上的残余强度，此时锚束体与围岩体之间产生相对运动，接触面之间只有残余强度存在。锚束体与围岩体界面本构模型如图1所示。 相应于2个阶段线性模型的剪应力–剪切位移关系可表示为 式中：u为剪切位移；K，ξ均为参数，对应于不同阶段的剪应力–剪切应变关系可分别表示为 式中：K1为界面刚度系数，τ1为锚固段极限黏结强度，τ2为锚固段残余黏结强度，umax为锚固段极限黏结强度对应的剪切位移。 3 理论方法的导出和实现 3.1 锚束体微段力学模型 当锚束体与围岩体接触面上的剪应力小于抗剪强度、界面处于完全黏结状态且没有相对位移发生时，两者之间应满足变形和应力协调关系，于是可采用弹性力学及其剪切滞模型的相关理论来研究界面特性。 受预应力轴向荷载作用的锚杆锚固段及影响范围内的围岩体是典型的弹性力学轴对称问题(见图2)。 考虑微元体轴向力的平衡条件，有如下方程成立： 式中：σzz为轴向应力，τrz为剪应力。 对式(6)沿a,b积分，有 式中：a为锚束体半径，b为围岩体内剪应力影响范围。 图3为锚束体及其影响岩体微段模型。考虑到τrz(b,z)=0，τrz(a,z)=τ(a,z)=τ(z)的条件，式(7)可进一步整理为 式中：为岩体内沿径向的平均轴向应力。 同样对于锚束体，存在如下关系： 式中：为锚束体平均轴向应力。 将式(8)代入式(6)可得 式(10)给出了剪应力沿径向的分布规律，其曲线见图3(b))。岩体本构关系为 式中：Gm为围岩体剪切模量，w为围岩体轴向位移。 将式(10)代入式(11)，并沿a,r积分可得 在此，引入平均位移差的概念，即 将式(12)代入式(13)积分并经过整理可得 其中， 将式(14)代入式(9)并进行微分可得 式中：Ef为锚束体复合体的弹性模量，可表示为 式中：Eg,Eh分别为钢筋和灌浆材料的弹性模量；ri为钢筋半径。 考虑到锚束体微段截面力的平衡条件(见图3): 于是有如下关系成立： 其中， 式(18)的通解为 相应的边界条件为 式中：P为作用在锚束体端部的外加荷载，l为锚杆锚固段长度。 求解式(19)可得 于是，式(19)可用下式表达： 相应地，锚束体与围岩体界面剪应力分布公式为 3.2 界面剪应力沿锚杆分布特性 一旦作用在锚束体与围岩体界面上的剪应力超过界面抗剪强度，界面就会发生破坏(脱黏)，根据界面本构关系，脱黏段的剪应力就为界面残余强度。作用在脱黏段上的剪应力就会重新分配和调整，于是可以计算出脱黏段的长度(见图4)。 如果不考虑界面脱黏情况，按照锚束体与围岩体完全胶结情况获得的剪应力沿锚杆分布特性见图4(a)。根据界面本构方程，界面剪应力并不可能按此规律分布。当界面剪应力超过界面抗剪强度时，