Dynamic response solutions of unsaturated soil foundation using soil-water characteristic curve considering deformation effects
-
摘要: 为研究移动荷载作用下的非饱和地基动力响应,在传统V-G土-水特征曲线模型的基础上,建立了考虑变形效应的土-水特征曲线模型,并基于该模型推导了新的非饱和土动力控制方程,从而完整描述非饱和土在受动力作用下的水-力耦合作用。进一步的,采用2.5维有限元法对控制方程进行求解,求解结果分别与单相弹性介质,双相饱和介质和三相非饱和介质的解析解进行对比,均验证了该求解方法的准确性;不同介质模型的计算耗时分析表明2.5维有限元法是目前求解多孔介质动力问题的一种优势算法。通过数值分析发现,采用传统的、未考虑变形效应的土-水特征曲线会低估非饱和地基的振动强度。Abstract: To investigate the dynamic response of unsaturated soils under moving loads, a new soil-water characteristic curve (SWCC) model is established considering the deformation caused by applied load based on the traditional V-G SWCC model. Using this modified SWCC model, a dynamic governing equation for unsaturated soils, which fully describes the water-force coupling effects of unsaturated soils under dynamic loading, is derived. The governing equation is solved using the 2.5-dimensional finite element method (2.5D FEM). The obtained solutions are compared with the analytical ones for single-phase medium, double-phase saturated medium and three-phase unsaturated medium, respectively, which all confirm the accuracy of the proposed solution method. The computational time analysis for different medium models demonstrate that the 2.5D FEM is an advantageous algorithm for solving the dynamic problems of porous media. The numerical analysis reveals that using the traditional SWCC without considering the deformation will underestimate the vibration intensity of unsaturated foundations.
-
0. 引言
近年来,随着城市化的进一步加深,基坑工程的建设朝着超大、超深的方向发展[1],深基坑降水和开挖将会引起地下水渗流场、土层中应力场的变化,只有正确认识基坑降水开挖过程中地连墙的受力变形机制,才能更好地指导深基坑的设计与施工。
目前,国内外已有许多学者采用多种方法对基坑开挖条件下地表沉降、地连墙变形进行研究。马险峰等[2]以上海某超深基坑为对象开展离心模型试验,研究了超深基坑开挖过程中地下连续墙侧移、坑外地层变形等规律;陈宝国等[3]通过模型试验方法,研究了内支撑体系、地连墙和墙后土体之间的协调变形特性,得出了墙后侧向土压力、坑外地表沉降等随支撑体系调节方式的变化规律;Fan等[4]假定围护结构变形曲线为分段二次多项式,推导出任意变形模式下的地表沉降解析解。
上述研究主要集中于基坑开挖卸荷引起的基坑变形及坑外地表沉降,忽略了基坑降水疏干对基坑施工造成的影响,当前基坑降水问题的研究重点在于坑内外水位变化及其对周围环境的影响。Goh等[5]通过有限元分析评估了影响地下水位下降引起地面沉降性能的关键参数,并开发了人工神经网络模型来估计最大地面沉降;Zeng等[6]基于现场抽水试验结果,采用流固耦合数值模型模拟基坑降水试验,揭示了地下水位下降和墙体移动共同作用下的地面沉降特征。
现阶段对于基坑降水对地连墙受力及变形的研究成果较少[7-8],且现有研究重点集中于基坑降水对地连墙的受力变形研究,缺乏降水开挖共同作用对地连墙受力变形的研究,而实际施工过程中基坑降水与开挖同时进行,地连墙会受到降水疏干及开挖卸荷的共同作用影响。现有研究大多采用数值模拟方法,其在考虑土与结构相互作用时所使用的假设条件过多,准确性有待进一步验证。而模型试验方法可以更好地还原工程条件,适合进行规律性研究,故常被用于基坑降水方面的研究。Xu等[9]进行了一系列室内试验研究,研究了地连墙的埋深和抽水井的过滤段长度对地下水位变化的影响;曾超峰等[10-12]开展了室内模型试验,揭示了基坑开挖前降水引发基坑变形的机制。
综上,当前大多数研究仅考虑基坑降水或基坑开挖单独作用对坑外水位变化、地表沉降及基坑变形的影响,而针对基坑降水与开挖共同作用下地连墙受力及变形规律的研究较少。本文为揭示基坑降水与开挖共同作用下地连墙受力及变形机制,开展了模型试验研究,模拟了基坑干砂开挖与降水开挖过程,通过对比不同条件下地连墙受力及变形的联系与区别,分析了基坑开挖卸荷与降水疏干对地连墙受力及变形的影响机制。
1. 模型试验系统
1.1 工程概况
如图 1所示为南昌地铁4号线某车站基坑平面布置及模拟区示意图,车站总长238 m,标准段宽22.7 m,标准段开挖深度约16 m,端头井段开挖深度约17.2 m。开挖范围属于强透水含水层,潜水位高程为4.10~7.04 m,故布置疏干井进行潜水疏干。地连墙采用C35钢筋混凝土,厚度800 mm,深度约为28 m。
本地铁车站基坑为狭长形基坑,仅针对如图 1所示的模拟区进行研究,即疏干井J7和观测井G8所在区域进行研究。所模拟区域为狭长形基坑中间段,不受到边角效应的影响,呈现出近似于平面应变的状态,因此,本模型试验为拟三维试验。
由于本模型试验是在1g缩尺条件下进行,难以满足所有物理力学相似条件,主要控制地连墙尺寸及地连墙抗弯刚度满足相似比。根据实验室场地条件及现有模型试验[10-12]取值经验,取相似比例为1∶40。值得注意的是,模型试验目的在于定性揭露地连墙的受力变形规律研究,而非对实际工程进行定量反推预测,研究重点也并非坑后土体的沉降特征,为简化试验,不对土体进行缩尺设计。
1.2 模型箱
如图 2所示为模型试验装置实物图,试验模型箱尺寸为260 cm×60 cm×120 cm(长×宽×高),为钢制无盖长方体箱体,长度方向为透明钢化玻璃观察窗,设置一道高强度槽钢加固箱体,限制箱体侧向变形。
如图 3所示为模型试验装置示意图,左侧为储水箱,储水箱与储土箱由透水隔板隔开,储水箱侧壁贴有刻度线,通过流量可调节抽水泵动态补给储水箱,严格控制储水箱内水位动态稳定,形成常水头边界。透水隔板全断面均匀打孔,开孔率为48%,与试验土层孔隙率相似,两侧粘贴过滤网布,防止砂土流失。模型箱右侧为基坑开挖区域,开挖区平面尺寸为30 cm×60 cm,最大开挖深度为40 cm。降水方式为坑内管井井点降水,选用外径为20 mm、内径为16 mm的PVC管模拟降水,井管长度为65 cm(实际工程降水井长度为26 m,缩尺40倍),过滤段长度为20 cm(实际工程过滤段长度为8 m,缩尺40倍),为保持过滤段强度以及足够的进水面积,开孔采用“梅花桩”型错位排列,为避免土粒在进水孔内堵塞,井管外包裹一层尼龙过滤网布进行反滤处理。采用SR400箱式蠕动抽水泵进行基坑降水,通过标定实验确定抽水泵最大流速为1.63 L/min,通过预降水试验检验降水井的降水效果,证明本模型降水井与水泵共同作用可以将坑内水位降至设计深度。根据Xu等[9]在中砂中开展的类似降水试验经验,本研究取降水影响半径与开挖影响半径均为最大开挖深度的4倍以上,取180 cm。
试验监测装置包括测压管、百分表、应变片、孔隙水压力传感器(型号:金诚测试BWK-S型,量程0~10 kPa,分辨率0.1 kPa,满量程测量误差0.01 kPa)和土压力传感器(型号:金诚测试BW型,量程为0~30 kPa,分辨率为0.1 kPa,满量程测量误差为0.15 kPa),测压管与箱底测压孔连接,通过对齐管内凹液面与刻度条进行读数,可达到浮漂监测方法[11-12]相同的数据读取效果,百分表安装于基坑外土体表面及地连墙顶部,应变片、孔隙水压力传感器和土压力传感器均按照图 3所示的空间排布粘贴于有机玻璃地连墙中线处。
1.3 试验土体
试验采用福建平潭ISO标准砂,填土厚度为110 cm,颗粒级配曲线如图 4所示,土体参数如下:d60=0.90 mm,d30=0.35 mm,d10=0.13 mm,计算可得:Cu= 6.92,Cc=1.05,根据《土的工程分类标准:GB/T50145— 2007》[13],可判定该砂级配良好,根据《建筑地基基础设计规范:GB50007—2011》[14],可判定该砂属于中砂。
模型箱中砂土填筑完成后,取样开展常规土工试验[15]得到砂土土性参数(均进行3次平行试验并取平均值)。土性参数如表 1所示,其中,湿密度根据环刀法得到,含水率根据烘干法得到,相对密度根据比重瓶法得到,孔隙比根据测得的湿密度、含水率、相对密度由公式计算得到,渗透系数根据常水头法得到,砂土内摩擦角及有效内摩擦角通过三轴压缩试验得到,压缩模量根据固结试验得到(如图 5所示)。
表 1 土性参数表Table 1. Physical parameters of model soil土层 ρ0/(g·cm-3) w0/% Gs e K/(m·d-1) φ/(°) φ′/(°) Es/MPa 砂土 2.01 29.6 2.67 0.72 18.68 22.8 31.6 11.67 1.4 支护结构
本地铁车站基坑施工区土质较好,坑顶钢筋混凝土支撑间距为12 m,基坑开挖模拟区长边长度为60 cm,根据相似比计算得支撑间距应为30 cm,故于基坑长边中心位置处布设一道支撑。根据实际工程支撑截面尺寸及现有模型试验经验[16],支撑材料选用PVC管材,弹性模量为3.5 GPa,泊松比为0.38,设计PVC管支撑的截面尺寸为20 mm×20 mm,支撑布设情况如图 3(a)所示。
根据现有模型试验经验[10-12, 17],地下连续墙材料选用有机玻璃板,地连墙按照抗弯刚度等效相似原则设计,地连墙有机玻璃板厚度计算公式为
dm=dpn(1−ν2m1−ν2p⋅EpnEm)13。 (1) 式中:E为弹性模量;d为地连墙厚度;ν为泊松比;下标p表示原型值;n为相似比(1∶40);下标m表示整体相似分析模型值。
实际工程中地下连续墙厚度为800 mm,长度为28 m,钢筋混凝土弹性模量为35 GPa,泊松比为0.17。本试验中,地连墙采用有机玻璃模拟,其弹性模量为3.37 GPa,泊松比为0.25,将上述参数代入式(1)中,可得有机玻璃板厚度为12 mm,深度为70 cm。
1.5 试验工况及步骤
如表 2所示,开展2个工况的模型试验,工况Ⅰ为降水开挖试验,降水开挖交替进行,先降水后开挖;工况Ⅱ为干砂开挖试验,分4层开挖。根据曾超峰等[10]在粗砂中进行的模型试验经验,当模型稳定降水10 min后,模型箱内部水位达到稳定、地表沉降和地连墙侧移达到稳定状态,故本模型试验每一降水、开挖步骤均稳定30 min。
表 2 模型试验工况Table 2. Programmes of model tests工况 开挖深度/cm 降水后水位/cm Ⅰ降水开挖 10,20,30,40 -15,-21,-31,-41 Ⅱ干砂开挖 10,20,30,40 — 注:初始地下水位为-15 cm,每层开挖前,将坑内地下水位降至开挖深度以下1 cm处。 以工况Ⅰ降水开挖试验为例,各工况试验均进行3次平行试验以减少试验误差,试验步骤如下。
(1)试验准备工作
首先对模型箱进行密封性试验,确保模型箱未出现渗水情况,将透水隔板和地连墙安装至设计位置,并在地连墙两侧用硬塑料薄膜做密封处理,砂土采用人工砂雨法填筑,控制砂土落距500 mm,落砂速度约120 g/s,每层铺约50 mm直至达到需要高度,填土过程中严格控制填土压实效果。当地连墙两侧同步填砂至地表以下3 cm时,架设支撑,继续填砂至地表高度。将各监测装置安装至设计位置,归零备用。
(2)注水加压,自重固结
连接水循环系统,调节供水蠕动泵转速使储水箱水头稳定在初始水位处,地下水流通过透水隔板进入储土箱,至水头达到初始水位处,停止注水,待各测压管水头稳定,可认为土层浸润饱和,静置24 h,待土层完成固结沉降,补充填土至设定标高。
(3)降水开挖试验
采用分层开挖的方式进行基坑开挖,各层开挖深度均为10 cm。本模型试验采用人工开挖,开挖前,记录各测压管读数为初始地下水位,维持储水箱内水位始终保持在-15 cm处,开启数据采集系统,通过抽水泵进行基坑降水,遵循“先降水,后开挖”的原则,降水至各开挖层低标高以下1 cm处,静止30 min,记录各监测装置读数,开挖至设定标高,静置30 min,记录各监测装置读数,重复上述过程,直至试验结束。
工况Ⅱ干砂开挖试验与上述步骤相似,区别在于干砂开挖试验不进行注水加压及基坑降水步骤。
2. 试验结果验证
Clough等[18]提出了基坑开挖“凹陷”型地表沉降经验曲线,可预测墙后不同距离处的地表沉降,Hsieh等[19]对地表沉降预测经验曲线进行了修正,并对沉降影响区域进行了划分。图 6为降水开挖与干砂开挖地表沉降归一化曲线与经验曲线对比图,其中δ为墙后任意点处地表沉降,δmax为墙后最大地表沉降,d为坑外测点与地连墙距离,Hmax为基坑开挖深度,取开挖至30 cm处结果进行验证。可以看出,模型试验结果与经验曲线吻合较好,干砂开挖工况下地表沉降影响范围与Clough经验曲线拟合效果更好;降水开挖工况下地表沉降影响范围与Hsieh经验曲线拟合效果更好,但最大沉降点出现位置更远,这是由于悬挂式地连墙条件下,降水开挖导致坑外地下水位下降,使得坑外土体产生附加固结沉降,造成了更大的地表沉降值和地表沉降影响范围。
3. 地连墙受力分析
监测点布置受模型试验尺寸所限,不能对研究剖面全线监测,故仅对数据监测点采集数据进行分析。土压力盒数据是在填土完成并初始平衡后进行采集的,因此获得的土压力监测数据是以静止土压力为初始点的变化数据,值得注意的是,墙侧土压力盒所监测到的数据均为侧向总土压力变化值,故本文后续侧向土压力均指侧向总土压力。静止土压力系数可通过Jaky[20]提出的计算公式获得:
K0=1−sinφ′。 (2) 叠加试验测得的侧向土压力数据,从而得到试验过程中的侧向土压力变化规律。无地下水时,基坑开挖后地连墙受到的侧向土压力是深度的函数:
p0=K0γz(干砂开挖)。 (3) 存在地下水时,基坑降水开挖后地连墙受到的侧向土压力需要分为水位线上、下两部分进行计算。对于砂土地基而言,水位线以下侧向土压力计算一般采用水土分算:
p0=K0γh+K0γ′(z−h)+γw(z−h)(降水开挖)。 (4) 式中:φ′为土的有效内摩擦角;γ为土的天然重度;γw为水重度;γ′为土的浮重度;h为水位埋深;z为监测点埋深,0<z<H;H为地连墙埋深;K0为侧向土压力系数。计算简图见图 7。
3.1 墙前/后侧向土压力(干砂开挖)
图 8为干砂开挖墙前/后侧向土压力沿深度分布图,其中z为监测点埋深,p为监测点处侧向土压力值,正值代表测点处侧向土压力指向坑外,负值代表侧向土压力指向坑内(若不做特殊说明,下文侧土压力方向也如此规定)。由于墙前侧向土压力变化显著,而开挖面以下监测点较少,故墙前侧向土压力图像采用散点图。由图 8(a)可以看出,在开挖深度范围内,墙前侧向土压力以减小为主,并随深度的增加,墙前侧向土压力减小幅度逐渐增大,几乎减小至0 kPa,其原因是基坑开挖导致开挖深度范围内作用于地连墙的初始侧向土压力消散;在开挖深度以下,墙前侧向土压力出现了显著的增长现象,这是由于开挖深度范围内墙前侧向土压力消散,地连墙受力平衡被打破,发生挠曲变形,使得开挖深度以下墙前侧向土压力形式由静止土压力变为被动土压力。
墙后侧向土压力亦将重分布,由图 8(b)可以看出,同一埋深处随着开挖深度的增加,墙后侧向土压力减小幅度逐渐增大;不同埋深处干砂开挖墙后侧向土压力以减小为主,并且侧向土压力减小幅度呈现出先增后减的趋势,减小幅度峰值均出现在埋深45 cm处。由于墙顶支撑的作用,地连墙顶部挠曲变形较小,墙后浅埋位置侧向土压力无明显变化,墙后埋深45 cm测点处侧向土压力减小幅度最大,与其位置处相对较大的挠曲变形有关。
图 9为干砂开挖墙前与墙后侧向土压力之差(△ pd=|pB|−|pA|)沿深度分布图,其中Δpd是干砂开挖地连墙侧向土压力之差,正值代表测点处Δpd指向坑内。可以看出,随着埋深的增加,Δpd的分布形式呈现出“S”型,以开挖深度40 cm为例,地连墙上部Δpd为正,地连墙受到指向坑内的土压力差,Δpd正峰值出现在地下30 cm处,为2.1 kPa;下部Δpd为负,地连墙受到指向坑外的土压力差,Δpd负峰值出现在地下60 cm处,为-3.8 kPa,Δpd负峰值约为正峰值的1.8倍。随着开挖深度的增加,Δpd正负峰值均逐步增大且出现下移趋势,这将导致地连墙挠曲变形峰值下移。
3.2 墙前/后侧向土压力(降水开挖)
图 10为降水开挖墙前/后侧向土压力沿深度分布图。由图 10(a)可以看出,同一埋深处随着降水开挖深度的增加,墙前侧向土压力减小幅度逐渐增大;不同埋深处降水开挖墙前侧向土压力均减小,开挖深度范围内墙前侧向土压力几乎减小至0 kPa,开挖深度以下墙前侧向土压力形式同样变为被动土压力。
墙后侧向土压力亦将重分布,由图 10(b)可以看出,降水开挖墙后侧向土压力变化形式与干砂开挖墙后侧向土压力变化形式相近,不同之处在于,降水开挖墙后侧向土压力最大减小幅度的埋深更深,墙后深埋位置处土压力减小幅度相较于干砂开挖更大。
降水开挖墙前/后侧向土压力分布与坑内降水导致的坑内外稳定的渗流场密切相关,将在3.4节中进行分析。
图 11为降水开挖墙前与墙后侧向土压力之差(△ pw=|pB|−|pA|)沿深度分布图,其中Δpw是降水开挖地连墙侧向土压力之差,正值代表测点处Δpw指向坑内。可以看出,随着埋深的增加,Δpw的分布形式呈现出“S”型,上部为正,地连墙受到指向坑内的土压力差,Δpw正峰值出现在地下45 cm处,为3.2 kPa;下部为负,地连墙受到指向坑外的土压力差,Δpw负峰值出现在地下60 cm处,负峰值为-4.9 kPa,Δpw负峰值约为正峰值的1.5倍。与干砂开挖情况不同的是,降水开挖条件下,随着开挖深度的增加,Δpw正峰值峰值位置变化较小。
3.3 墙前/后孔隙水压力
图 12为降水开挖作用于墙前/后的墙侧孔隙水压力沿深度分布图。根据试验工况可知,每层开挖前先将坑内水位降至开挖深度下1 cm处,图 12中,正值代表地连墙墙前测点处孔压,负值代表地连墙墙后测点处孔压(孔压正负值仅代表与地连墙的相对位置)。由图 12(a)可以看出,同一埋深处,随着降水开挖深度的增加,墙前孔隙水压力逐渐减小,在降水完成后地下水位以上各监测点处孔压值均减小至0 kPa左右;随埋深的增加,墙前孔压减小幅度呈现出先增后减的趋势,在地下45 cm处出现减小幅度峰值,这与坑内降水引起的坑内外稳定的渗流场有关,过滤段位于坑内埋深45~65 cm范围内,可视为渗流终点,渗流终点处孔压减小幅度最大[21],随着距离渗流终点距离的增加,孔压减小幅度越小。
坑内外稳定的渗流场亦将导致墙后孔压重分布,由图 12(b)可以看出,降水开挖墙后孔压变化形式与墙前孔压变化形式一致,墙后孔压减小幅度随开挖深度的增加而逐渐增大;但随埋深的增加,墙后孔压减小幅度逐渐增大,这是由于埋深越深,距离渗流终点越近。埋深70 cm测点处,墙后孔压略大于墙前孔压,说明存在从坑外流向坑内的渗流场。这与曾超峰等[10]的研究成果一致。
图 13为降水开挖作用于墙前/后的墙侧孔隙水压力之差(△ u=|uB|−|uA|)沿深度分布图,Δu为地连墙两侧孔压差,正值代表测点处孔压差Δu指向坑内。可以看出,随着埋深的增加,Δu分布形式为“单峰”型,地连墙主要受到指向坑内的孔压差,表明墙前孔压减小幅度普遍大于同埋深处墙后孔压减小幅度,Δu峰值出现在埋深30~45 cm处,这是由于此范围处墙前/后渗流路径最长,所以孔压差Δu最大。
3.4 降水开挖与干砂开挖土压力对比分析
为对比分析降水开挖与干砂开挖情况下地连墙受力情况的区别,取开挖深度为30,40 cm的侧向土压力变化值进行研究,图 14为降水开挖与干砂开挖墙前/后侧向土压力变化值沿深度分布对比图,其中δp为侧向土压力变化值,负值代表监测点位置处侧向土压力减小。
由图 14(a)可以看出,在开挖深度范围内,墙前侧向土压力均以减小为主,并随埋深的增加,减小幅度逐渐增大;在开挖深度以下,降水开挖与干砂开挖墙前侧向土压力变化形式不同,随着埋深的增加,干砂开挖情况下δp出现正值,表明墙前侧向土压力超过了初始静止土压力值;而降水开挖情况下δp未出现正值,表明墙前土压力始终未超过初始静止土压力。这与降水开挖情况下墙前孔压的变化情况有关,如图 12(a)所示,由于降水井降水作用(过滤段位于埋深45~65 cm范围内,可视为渗流终点),墙前开挖深度以下孔压大幅度降低,进而降低了墙前开挖深度以下的侧向总土压力,从而导致δp始终为负值(被动土压力值相对于干砂开挖工况更小)。
由图 14(b)可以看出,降水开挖与干砂开挖墙后侧向土压力变化形式相似,随着埋深的增加,墙后侧向土压力变化分布形式为“单峰”型,不同之处在于,侧向土压力变化幅度峰值出现位置及大小有所不同。以开挖深度为40 cm为例,干砂开挖墙后侧向土压力峰值出现在埋深45 cm处,为-3.9 kPa;降水开挖峰值出现在埋深60 cm处,为-4.9 kPa。这与基坑坑内降水导致的坑内外形成的稳定渗流场有关,从图 12(b)可以看出,墙后地下水存在自上而下的渗流,随着埋深的增加,孔压减小幅度逐渐增大,从而导致墙后侧向土压力减小幅度峰值位置较干砂开挖有所下移,且减小幅度峰值更大。
图 15为降水开挖与干砂开挖侧向土压力之差沿深度分布对比图,其中,Δp为地连墙受到的侧向土压力之差,正值代表测点处Δp指向坑内。可以看出,两种情况下Δp的变化形式相近,不同之处在于,Δp峰值位置和大小有所不同,以开挖深度为40 cm为例,Δpd正峰值出现在埋深30 cm处,Δpw正峰值出现在埋深45 cm处;两者负峰值出现在埋深60 cm处,Δpw正负峰值均更大,这将导致更大的地连墙弯矩及挠曲变形。
4. 地连墙变形分析
4.1 地连墙弯矩
根据马险峰等[2]的研究,地连墙弯矩可通过粘贴于地连墙中线上的应变片监测数据计算得到。计算公式如下:
M(z)=EID△ ε(z)。 (5) 式中:△ ε(z)为地连墙在埋深z两侧应变片的差值;D为地连墙厚度;EI为地连墙抗弯刚度;M(z)为z深度处的地连墙弯矩。
图 16为降水开挖与干砂开挖地连墙弯矩沿深度分布图,其中M为地连墙弯矩,指向基坑内部为正,相反为负。可以看出,降水开挖与干砂开挖地连墙弯矩变化趋势相同:当开挖深度为10 cm时,地连墙弯矩为负,且弯矩值较小;当开挖深度超过20 cm时,地连墙弯矩随深度变化形式改变,弯矩分布形式呈“S”型,上部地连墙弯矩指向坑内,M正峰值出现在地下15 cm处,由图 16(a),(b)可以看出,干砂开挖情况下M正峰值达到18.1 kN·m/m,降水开挖情况下M正峰值达到29.9 kN·m/m;下部地连墙弯矩指向坑外,M负峰值位置随开挖深度的增加而降低,以开挖至埋深40 cm处为例,两种情况下M负峰值均出现在地下60 cm处,干砂开挖M负峰值达到-19.4 kN·m/m,降水开挖M负峰值达到-27.8 kN·m/m。从图 15侧向土压力之差的对比分析可得出,降水开挖的侧向土压力之差大于干砂开挖的侧向土压力之差,地连墙的弯矩变化情况与相应的地连墙受力变化情况一致。
4.2 地连墙挠曲变形
根据现有模型试验经验[2, 22],地连墙挠曲变形可通过粘贴于地连墙中线的应变片监测数据计算得到,以地连墙底部为基点,将应变片监测数据换算为墙体的曲率分布,进而计算得到墙体挠曲变形随深度的变化,计算简图如图 17所示。挠度公式推导过程如下。
根据平面曲线的曲率公式的近似形式,可得每一微分段dL的曲率半径ρ及转角θ:
θ=tanθ=dLρ=εdL。 (6) 对于第1个测点,极坐标系与直角坐标系一致,则第1个测段的挠曲变形值为
dx1=−ρ1(1−cosθ1)。 (7) 对于第i个测点,极坐标系相当于初始直角坐标系转动角度为Θi,则各测段的挠曲变形值为
dxi=(dL−ρisinθi)sinΘi−ρi(1−cosθi)cosΘi 。 (8) 累加可得第i个测点处地连墙侧向挠曲变形值:
xi=i∑1 dxi=dx1+dx2+⋯+dxi。 (9) 式中:ρ为每一微分段的曲率半径;θ为相邻微分段之间的转角;dL为微分段长度;dxi为第i个微分段的挠曲变形;xi为第i个微分段的累计挠曲变形。
上述过程需满足以下条件:①地连墙侧移属于小变形范畴;②地连墙仅考虑挠曲变形;③地连墙为多段圆弧组成的光滑曲线,段内各点曲率相等。
图 18为降水开挖与干砂开挖地连墙侧移沿深度分布图,其中δv为地连墙挠曲变形,指向基坑内部的挠曲变形为正,相反为负。可以看出,降水开挖与干砂开挖地连墙挠度变化趋势相同:地连墙挠曲变形值随着开挖深度增大而增大,挠曲变形呈“内凸”型,这是由于坑顶支撑的作用限制了地连墙顶部的挠曲变形。随着开挖深度的增大,地连墙挠曲变形峰值位置逐步下移,且均出现在开挖深度附近,这与徐中华等[23]对上海软土地区的统计分析结果一致。对比图 18(a),(b)可知,降水开挖情况下的地连墙挠曲变形值均大于干砂开挖情况,如当开挖深度为40 cm时,降水开挖情况下地连墙侧移峰值达到0.024 mm,而干砂开挖情况下地连墙侧移峰值仅为0.014 mm。地连墙侧移趋势与地连墙弯矩变化及相应的地连墙受力变化趋势一致。
5. 结论
本文为揭示基坑降水与开挖共同作用下地连墙受力变形机制,开展了室内模型试验,模拟了基坑干砂开挖与降水开挖过程,通过对比分析不同试验工况得到以下4点结论。
(1)基坑开挖卸荷改变了基坑内外土体中的应力场,墙两侧侧向土压力均以减小为主,墙前部分监测点出现显著的增长现象(被动土压力区),侧向土压力之差分布形式呈“S”型,而降水开挖工况试验过程被动土压力则更小,这由坑内降水引起的墙前孔压大幅度减小导致。
(2)相较于干砂开挖工况,降水开挖工况下墙后侧向土压力减小幅度更大且峰值位置下移,侧向土压力之差峰值更大且峰值位置下移,这由坑外向坑内稳定的渗流导致。
(3)地连墙弯矩的变化较好地吻合了地连墙受力情况的变化,相对于干砂开挖试验结果,降水开挖工况在地连墙顶部支撑及侧向水、土压力的共同作用下,降水开挖条件下的地连墙正、负弯矩值均更大。
(4)基坑干砂开挖与降水开挖条件下,地连墙侧移发展规律相同,挠曲变形形式均呈“内凸”型,降水开挖情况下地连墙侧移值更大,与地连墙侧向土压力分布情况及弯矩分布情况吻合。
值得注意的是,本文研究中土层条件为单一砂层,地连墙的受力变形均较小,未达到极限状态,故无法对开挖后的主动/被动侧水土分算和水土合算进行对比分析。未来可开展基于砂土-黏土复杂地层基坑降水开挖诱发地连墙受力变形研究。
-
Ks/GPa Kw/GPa Ka/kPa G/MPa n ρs/(kg·m-3) ρw/(kg·m-3) ρa/(kg·m-3) Sr 35 2.25 145 25 0.125 0.001 2500 0.001 0.001 0.001 Kb/MPa Ks/GPa Kw/GPa Ka/kPa G/MPa n ρs/(kg·m-3) ρw/(kg·m-3) ρa/(kg·m-3) Sr kD/(m·s-1) 8.77 11 2.25 145 2997.6 0.125 0.3 2500 1000 0.001 1 1×10-6 Kb/MPa Ks/GPa Kw/GPa Ka/kPa Gs/MPa n ρs/(kg·m-3) ρw/(kg·m-3) ρa/(kg·m-3) 8.33 35 2.25 145 3.85 0.35 0.45 2650 1000 1.28 Sr Sa ηw/(Ns·m-2) ηa/(Ns·m-2) Sw0 κ/m2 β/(Pa-1) m k φ' 0.5,0.8 0.5,0.2 0.001 18×10-6 0.05 5.3×10-13 1.0×10-4 0.5 2 10° -
[1] 黄长生, 周耘, 张胜男, 等. 长江流域地下水资源特征与开发利用现状[J]. 中国地质, 2021, 48(4): 979-1000. HUANG Changsheng, ZHOU Yun, ZHANG Shengnan, et al. Groundwater resources in the Yangtze River Basin and its current development and utilization[J]. Geology in China, 2021, 48(4): 979-1000. (in Chinese)
[2] 吴庆华, 汪啸, 范越. 长江中下游地下水资源战略储备选址适宜性评价指标体系[J]. 长江科学院院报, 2022, 39(8): 145-151, 158. WU Qinghua, WANG Xiao, FAN Yue. Evaluation index system for the suitability of groundwater strategic reserve site in middle and Lower Yangtze River[J]. Journal of Yangtze River Scientific Research Institute, 2022, 39(8): 145-151, 158. (in Chinese)
[3] 池田俊雄. 地盤地質と鉄道土木の 50年[J]. 応用地質, 1998, 39(1): 4-10. IKEDA T. 50 years of geology and railway engineering [J]. Applied Geology, 1998, 39(1): 4-10. (in Japanese)
[4] Conrad, British Columbia. Canadian National Train No. Q-102-51-26 Derailment[R]. Canada: Transportation Safety Board of Canada, 1997: 1-24.
[5] EASON G. The stresses produced in a semi-infinite solid by a moving surface force[J]. International Journal of Engineering Science, 1965, 2(6): 581-609. doi: 10.1016/0020-7225(65)90038-8
[6] SHENG X, JONES C J C, PETYT M. Ground vibration generated by a load moving along a railway track[J]. Journal of Sound and Vibration, 1999, 228(1): 129-156. doi: 10.1006/jsvi.1999.2406
[7] 雷晓燕, 徐斌, 徐满清. 半无限弹性空间中移动荷载动力响应的频域-波数域比例边界有限元法分析[J]. 振动工程学报, 2017, 30(5): 798-805. LEI Xiaoyan, XU Bin, XU Manqing. Using the frequency-wave domain scaled boundary finite element method for the dynamic response of the elastic half space due to moving loads[J]. Journal of Vibration Engineering, 2017, 30(5): 798-805. (in Chinese)
[8] BIOT M A. Mechanics of deformation and acoustic propagation in porous media[J]. Journal of Applied Physics, 1962, 33(4): 1482-1498. doi: 10.1063/1.1728759
[9] LU J F, JENG D S. A half-space saturated poro-elastic medium subjected to a moving point load[J]. International Journal of Solids and Structures, 2007, 44(2): 573-586. doi: 10.1016/j.ijsolstr.2006.05.020
[10] CAI Y Q, SUN H L, XU C J. Response of railway track system on poroelastic half-space soil medium subjected to a moving train load[J]. International Journal of Solids and Structures, 2008, 45(18/19): 5015-5034. http://core.ac.uk/download/pdf/82345375.pdf
[11] GAO G Y, CHEN Q S, HE J F, et al. Investigation of ground vibration due to trains moving on saturated multi-layered ground by 2.5D finite element method[J]. Soil Dynamics and Earthquake Engineering, 2012, 40: 87-98. doi: 10.1016/j.soildyn.2011.12.003
[12] BIAN X C, HU J, THOMPSON D, et al. Pore pressure generation in a poro-elastic soil under moving train loads[J]. Soil Dynamics and Earthquake Engineering, 2019, 125: 105711. doi: 10.1016/j.soildyn.2019.105711
[13] 胡静, 唐跃, 张家康, 等. 高速列车荷载作用下饱和软土地基动力响应研究[J]. 岩土力学, 2021, 42(11): 3169-3181. HU Jing, TANG Yue, ZHANG Jiakang, et al. Dynamic responses of saturated soft soil foundation under high speed train[J]. Rock and Soil Mechanics, 2021, 42(11): 3169-3181. (in Chinese)
[14] 徐明江. 非饱和土地基与基础的动力响应研究[D]. 广州: 华南理工大学, 2010. XU Mingjiang. Study on Dynamic Response of Unsaturated Soil Foundation and Foundation[D]. Guangzhou: South China University of Technology, 2010. (in Chinese)
[15] MUALEM Y. A new model for predicting the hydraulic conductivity of unsaturated porous media[J]. Water Resources Research, 1976, 12(3): 513-522. doi: 10.1029/WR012i003p00513
[16] LU Z, FANG R, YAO H L, et al. Dynamic responses of unsaturated half-space soil to a moving harmonic rectangular load[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2018, 42(9): 1057-1077. doi: 10.1002/nag.2780
[17] FANG R, LU Z, YAO H L, et al. Study on dynamic responses of unsaturated railway subgrade subjected to moving train load[J]. Soil Dynamics and Earthquake Engineering, 2018, 115: 319-323. doi: 10.1016/j.soildyn.2018.08.037
[18] GAO G Y, YAO S F, YANG J, et al. Investigating ground vibration induced by moving train loads on unsaturated ground using 2.5D FEM[J]. Soil Dynamics and Earthquake Engineering, 2019, 124: 72-85. doi: 10.1016/j.soildyn.2019.05.034
[19] 李绍毅. 土体饱和度对移动荷载引起多层非饱和铁路地基振动的影响[J]. 岩土力学, 2021, 42(1): 151-159, 167. LI Shaoyi. Influences of soil saturation on the vibration of multi-layer unsaturated railway ground induced by moving load[J]. Rock and Soil Mechanics, 2021, 42(1): 151-159, 167. (in Chinese)
[20] Corey R H Brooks and A T. Hydraulic properties of porous media and their relation to drainage design[J]. Transactions of the ASAE, 1964, 7(1): 26-28. doi: 10.13031/2013.40684
[21] VAN GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated Soils1[J]. Soil Science Society of America Journal, 1980, 44(5): 892. doi: 10.2136/sssaj1980.03615995004400050002x
[22] FREDLUND D G, XING A Q. Equations for the soil-water characteristic curve[J]. Canadian Geotechnical Journal, 1994, 31(4): 521-532. doi: 10.1139/t94-061
[23] 孙德安. 非饱和土力学特性及本构模型[J]. 岩土工程学报, 2023, 45(1): 1-23. doi: 10.11779/CJGE20221450 SUN De'an. Mechanical behaviors and constitutive model for unsaturated soils[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(1): 1-23. (in Chinese) doi: 10.11779/CJGE20221450
[24] 蔡国庆, 田京京, 李舰, 等. 考虑变形及滞回效应影响的三维土-水特征曲面模型[J]. 土木工程学报, 2019, 52(11): 97-107. CAI Guoqing, TIAN Jingjing, LI Jian, et al. A three-dimensional soil water characteristic surface model considering deformation and hysteresis effect[J]. China Civil Engineering Journal, 2019, 52(11): 97-107. (in Chinese)
[25] 胡冉, 陈益峰, 周创兵. 考虑变形效应的非饱和土相对渗透系数模型[J]. 岩石力学与工程学报, 2013, 32(6): 1279-1287. doi: 10.3969/j.issn.1000-6915.2013.06.023 HU Ran, CHEN Yifeng, ZHOU Chuangbing. A relative hydraulic conductivity model for unsaturated deformable soils[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(6): 1279-1287. (in Chinese) doi: 10.3969/j.issn.1000-6915.2013.06.023
[26] 张雪东, 赵成刚, 刘艳, 等. 变形对土水特征曲线影响规律模拟研究[J]. 土木工程学报, 2011, 44(7): 119-126. ZHANG Xuedong, ZHAO Chenggang, LIU Yan, et al. Modeling study of the relationship between deformation and water retention curve[J]. China Civil Engineering Journal, 2011, 44(7): 119-126. (in Chinese)
[27] YANG Y B, HUNG H H. A 2.5D finite/infinite element approach for modelling visco-elastic bodies subjected to moving loads[J]. International Journal for Numerical Methods in Engineering, 2001, 51(11): 1317-1336. doi: 10.1002/nme.208
[28] 王勖成. 有限单元法[M]. 北京: 清华大学出版社, 2003. WANG Xucheng. Finite Element Method[M]. Beijing: Tsinghua University Press, 2003. (in Chinese)
[29] HALL L. Simulations and analyses of train-induced ground vibrations in finite element models[J]. Soil Dynamics and Earthquake Engineering, 2003, 23(5): 403-413. doi: 10.1016/S0267-7261(02)00209-9
-
其他相关附件