圆形水工隧道应力和洞壁位移解析解.docxVIP

圆形水工隧道应力和洞壁位移解析解 0 压力隧道的荷载 地下工程中的坍塌钻是设计中常用的隧道结构，通常设计为圆形。隧洞在开挖前,岩体处于静止、平衡状态,开挖破坏了初始平衡状态,隧洞周围岩体应力发生重分布,达到二次应力状态。地下水是岩体赋存的重要环境因素之一,由于地下水的作用,常引起岩体应力重分布,从而影响岩体的变形与破坏,关系到隧洞工程的稳定和安全。对于水工压力隧洞,运行期均匀内水压力是主要荷载。外水压力是渗流过程中作用在围岩上的渗流体积力,且受隧洞施工和运行的影响,有自己的荷载历史。岩石、硬黏土、结构性土、密实性砂土和混凝土等材料在应力–应变关系上具有明显的峰值,峰值后应力随变形增大而降低,称为应变软化,之后趋于稳定的数值–残余强度。在体积应变与主应变曲线上,随着主应变的增大,体积应变先是减小,随后逐渐增大,称为剪胀。目前压力隧洞的研究多是集中在某一方面,如强度准则的选取、不同工况下主应力顺序的变化、岩石材料的应变软化和剪胀特性、渗流体积力的影响等,所得结果与实际情况存在差异。本文针对圆形水工压力隧洞,采用统一强度理论和弹脆塑性软化模型,考虑不同工况下主应力顺序、岩石应变软化、剪胀和渗流等综合影响,推导了隧洞应力和洞壁位移解析解,并通过工程算例分析本文方法与传统方法的不同,及各参数对隧洞切向应力和洞壁位移的影响。 1 岩相剪化特性理论参数 统一强度理论能合理地考虑中间主应力σ2效应,适用于各种拉压特性不同的材料。由于圆形隧洞纵向尺寸远大于其横向尺寸,且荷载和几何形状均是轴对称的,可认为隧洞处于平面应变状态,σz为中间主应力,且在塑性区σz=σ2=(σ1+σ3)/2。对于平面应变弹塑性问题,统一强度理论表达式为 式中c,?分别为材料的黏聚力和内摩擦角;tc,?t分别为材料的统一黏聚力和统一内摩擦角;σt,σc,τs和α分别为材料的抗拉强度、抗压强度、抗剪强度和拉压比;b为反映中间主剪应力及相应面上的正应力对材料屈服或破坏影响程度的参数,称为统一强度理论参数,0≤b≤1。 为突出反映岩土材料达到强度峰值后的应变软化特性,峰值至残余强度用应力跌落表示,称为弹脆塑性软化模型;用两折线分别模拟体积应变εv、第三主应变ε3与第一主应变ε1的关系,如图1所示。图中,ε1p,ε3p,εvp为峰值后局部坐标中的塑性应变,以拉伸为正,压缩为负;h,H为剪胀特性参数,且h=1+H。 采用平面应变状态下的统一强度理论,其塑性软化前、后的屈服准则为 式中cj,?j分别为材料的黏聚力和内摩擦角;ctj,?tj分别为材料的统一黏聚力和统一内摩擦角;j=0表示材料应变软化前,j=1表示材料应变软化后。 2 下挖深度持续型 岩体的渗流主要是地下水在裂隙(或孔隙)内的流动,一般情况下属层流,符合达西(Darcy)定律。设在隧洞施工期,距隧洞中心半径0R范围以外,外水压力与原始渗流场外水压力ip相同,如图2所示。 运行期由于隧洞内半径r1处作用渗流水压力pm,隧洞保持原始渗流场外水压力pi的半径将由施工期的R0调整为Rc,如图3所示。 施工期和运行期渗流体积力分别为 式中,pw,p′w分别为隧洞施工期和运行期的渗流体积力,R0,Rc可由钻孔实测或抽水试验确定。 3 工程期间，对围岩力学和洞壁位移的分析 3.1 围岩应力模型 设围岩为弹脆塑性软化材料,且在无穷远处受初始地应力q作用。隧洞在不同的工况下,第一主应力可能是径向应力,也可能是切向应力。在施工期,围岩由于开挖卸载进入塑性变形,此时径向应力σr为第一主应力。施工期含水围岩破坏时,剪应力τ与正应力σ+ηpw满足平面应变状态下的统一强度理论为 式中,η为渗流体积力的有效面积(单位面积上的有效面积)系数。为安全起见,在研究岩体的破坏和稳定时,一般取η=1。 施工期隧洞应力计算简图,如图4所示。图中,2r为施工期围岩塑性区外半径。 视含水围岩为两相介质体,由于渗流而形成渗流体积力场,于是平衡微分方程为 由式(5)(j=1)和边界条件解式(6)得 式中,ct1,?t1分别为围岩应变软化后的统一黏聚力和统一内摩擦角。 设围岩弹塑性交界处的径向压力为pr2,则围岩弹性区应力为 式中,pw为施工期围岩的渗流体积力,应取式(3)。 在r=r2处,围岩弹性区的应力满足初始屈服方程式(5)(j=0),整理得 式中,ct0,?t0分别为围岩初始的统一黏聚力和统一内摩擦角。 在r=r2处,根据弹塑性区径向应力连续,由式(7a)和式(9)得 3.2 围岩塑性区和弹性区应力 式中,E,μ分别为围岩的弹性模量和泊松比。 采用简化的应变关系曲线,见图1(b)和图1(c)。在塑性区, 由式(11)和式(13),并考虑到εre=-εθe,得 以式(12)为边界条件,求得式(14)通解,进而得隧洞洞壁位移 由式(10)计算出r2后,再由式(7)~