基于PDV模型对一维“锅盖效应”气态水迁移的研究

唐国航; 王乃东; 姚仰平

doi:10.11779/CJGE2022S2036

基于PDV模型对一维“锅盖效应”气态水迁移的研究 English Version

唐国航^{1, 2,},
王乃东^1, ,,
姚仰平¹

1.
北京航空航天大学交通科学与工程学院, 北京 100191
2.
同济大学地下建筑与工程系, 上海 200092

基金项目:

国家自然科学基金面上项目 51978028

详细信息

作者简介:
唐国航(1995—)，男，新疆乌鲁木齐人，博士研究生，主要从事岩土工程土工试验方面研究。E-mail: TangGuohang@tongji.edu.cn

通讯作者:
王乃东, E-mail：wangnd@buaa.edu.cn

中图分类号: TU43
计量
- 文章访问数: 129
- HTML全文浏览量: 19
- PDF下载量: 14
出版历程
- 收稿日期: 2022-12-05
- 网络出版日期: 2023-03-26
- 刊出日期: 2022-11-30

One-dimensional "pot cover effect" for migration of gaseous water based on PDV model

1.
Department of Civil Engineering, Beihang University, Beijing 100191, China
2.
Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China

摘要

摘要: “锅盖效应”是一种重要的土中水汽迁移现象，常发生于季节性冻土区域，带来许多工程病害的发生。在目前“锅盖效应”试验与理论的研究、非饱和土体水汽迁移模型的基础上，考虑多种影响因素，基于PDV模型提出并推导了一维“锅盖效应”气态水迁移方程，用于描述或预测“锅盖效应”现象，并分析了温度与饱和度对“锅盖效应”迁移通量的影响。
- 锅盖效应 /
- 水汽迁移 /
- 迁移通量 /
- 浓度梯度
Abstract: The "pot cover effect" is an important phenomenon of water vapor migration in soils, which often occurs in the seasonal permafrost regions and causes many engineering diseases. Based on the current "pot cover effect" experiments and theoretical studies and unsaturated models for migration of soil moisture, considering the multiple influence factors, a one-dimensional "pot cover effect" gaseous water transport equation is proposed by using the PDV model. It is used to describe or predict the "pot cover effect" phenomenon, and the influences of the temperature and degree of saturation of soils on the "pot cover effect" of migration flux are investigated.
- pot cover effect /
- water vapor migration /
- migration flux /
- concentration gradient

HTML全文

0. 引言

随着城市建设的快速发展和建筑技术的不断进步,城市地下空间得到大规模开发和利用,超深超大基坑工程不断涌现。近年来,基坑开挖面积在10×10⁴ m²以上,开挖深度达40 m以上的工程项目越来越多。例如,上海地铁4号线修复工程深基坑开挖深度接近41 m,上海世博地下变电站基坑开挖深度接近34 m,天津的高层建筑基坑最大平面尺寸已达548 m×187 m。同时,基坑工程中也出现了一些亟待解决的问题,坑底土体回弹问题就是其中之一^[1-4]。基坑开挖产生坑底回弹,同时周边围护结构变形,也会要造成基底的隆起,回弹再压缩变形往往是建（构）筑物沉降变形的主要组成部分。

已有学者对基坑回弹变形的特点及其对坑内桩基承载力的影响进行了研究。Iwasaki等^[5]最早认为坑底土体竖向回弹会影响坑底桩基,认为土体回弹会对桩产生向上的侧摩阻力从而导致桩身上拔。对于基坑开挖与桩基的相互作用,Finno等^[6]利用平面有限元对一基坑开挖引起临近桩基破坏的实例进行了分析。刘畅等^[7]采用有限元数值模拟分析结合现场实测数据,研究了采用逆作法施工的基坑回弹变形问题,分析了工程桩、支护结构、楼板对坑底回弹变形的影响。查甫生等^[8]通过有限元软件ABAQUS,以坑底无桩和坑底群桩两种基坑为研究对象,对比分析了有桩、无桩情况下,深基坑开挖卸载的变形特性,得出了工程桩可以使基坑周围沉降、基坑中心隆起、围护结构向坑内位移这几种变形明显减小。曹力桥^[9]利用ABAQUS软件分析了存在工程桩和不存在工程桩基坑开挖和降水下的三维模型,通过对比分析了基坑开挖降水过程中基坑隆起的基本规律,得出工程降水对深基坑土体的压密作用及工程桩对坑底变形有明显的抑制作用。冯虎等^[10]利用FLAC数值模拟软件研究了坑内工程桩对软土超深基坑抗隆起稳定的影响规律以及作用机理,结果表明,墙趾土层特性、地连墙插入深度、基坑宽度和潜在滑裂面之内的工程桩对基坑抗隆起稳定有着非常显著的影响。

本文结合江苏某隧道明挖基坑工程,利用PLAXIS 3D软件,采用小应变土体硬化（HSS）模型作为土层的本构模型,建立了太湖隧道第二仓基坑的1/4模型,应用该模型,分析研究了坑底工程桩的桩长、桩径、桩刚度对基坑回弹变形的影响规律。

1. 工程概况

江苏某湖底隧道工程,公路等级为双向六车道（全线紧急停车带）高速公路,设计速度100 km/h,隧道总长10709 m,净宽16.75 m,湖中最大开挖深度达15 m,建成断面示意图见图1。

图 1 隧道断面示意图

Figure 1. Schematic diagram of tunnel section

下载: 全尺寸图片幻灯片

本文研究的基坑地层分布较均匀,项目隧道工程场地主要为粉土及粉质黏土,局部夹软土层,设计采用明挖施工,主要采用放坡加垂直支护形式进行开挖,基坑工程规模大、施工时间长、施工工序和工艺流程复杂等特点,研究基坑开挖卸荷引起的坑底回弹变形是该项目长期变形控制的关键技术问题之一。

现场监测由施工方进行,由于在施工过程中隆起测点遭到施工破坏,所以图2只给出了第二仓基坑孔隙水压力、立柱隆沉、桩土深层水平位移、围护结构水平位移及支撑轴力的测点图。

图 2 第二仓基坑测点布置图

Figure 2. Survey point layout of the second warehouse foundation pit

下载: 全尺寸图片幻灯片

2. 数值模型的建立

数值模拟标段选用第二仓K25+135—K25+515标段,基坑尺寸为400 m×80 m,开挖深度为15 m,为了便于数值计算,选取1/4的基坑进行建模,基坑三维模型及网格划分如图3所示。

图 3 基坑三维模型及网格划分示意图

Figure 3. Three-dimensional model and grid drawing of foundation pit

下载: 全尺寸图片幻灯片

土层分层情况已进行适当简化,已在断面图4中标明。自上而下分别为：2-1粉质黏土,2-3粉土,2-4淤泥质粉质黏土,3-1粉质黏土,3-2粉质黏土,4-1b粉质黏土,4-1黏土,层底标高-90 m。

图 4 基坑断面示意图

Figure 4. Schematic diagram of excavation section

下载: 全尺寸图片幻灯片

2.1 土体本构模型

基坑开挖一般属于临时性工程,工期较短,所以按不排水条件进行分析,且不考虑开挖过程对土体扰动的影响；土体本构模型采用小应变硬化模型。通过标准固结试验获得土体参考切线模量 $E_{oed}^{ref}$ $E_{oed}^{ref}$ ,通过三轴固结排水剪切试验获得土体参考割线模量 $E_{50}^{ref}$ $E_{50}^{ref}$ 、破坏比R_f和土体强度参数c′,φ′值,通过三轴固结排水加卸载剪切试验获得参考加卸载模量 $E_{ur}^{ref}$ $E_{ur}^{ref}$ ,具体试验过程与试验数据处理不再赘述。HSS模型参数取值见表1。

表 1 土层小应变本构模型参数取值表

Table 1. Parameter value table of soil layer HSS model

土层	$E_{oed}^{ref}$ $E_{oed}^{ref}$ /MPa	$E_{50}^{ref}$ $E_{50}^{ref}$ /MPa	$E_{ur}^{ref}$ $E_{ur}^{ref}$ /MPa	$G_{0}^{ref}$ $G_{0}^{ref}$ /MPa	γ_0.7/10^-4	R_f	$c^{'}$ $c^{'}$ /kPa	$φ^{'}$ $φ^{'}$ /(°)
2-1	5.067	3.24	41.32	112.0	2.0	0.70	5.0	37.5
2-3	8.592	4.03	26.93	100.0	2.0	0.62	4.7	32.2
2-4	2.190	7.32	29.11	58.1	2.0	0.90	7.2	29.0
4-1	2.566	7.06	32.87	115.8	2.0	0.86	16.5	38.0

下载: 导出CSV

| 显示表格

2.2 围护结构及工程桩模型

基坑采用上部放坡与钻孔灌注桩结合的围护结构形式,如图4所示,围护结构及支撑均按弹性材料考虑,三维模型见图5。

图 5 支护体系模拟示意图

Figure 5. Schematic diagram of support system simulation

下载: 全尺寸图片幻灯片

（1）采用两道内支撑,第1道为钢筋混凝土支撑,间距为8 m,截面为800 mm×1000 mm,刚度为3×10⁷ kPa,重度为25 kN/m³；第2道支撑为钢管支撑,间距为4 m,截面直径为609 mm,厚度为16 mm,刚度为2.1×10⁸ kPa,重度为78.5 kN/m³。

（2）止水帷幕及连续墙采用板单元进行模拟,本文在计算时将连续墙的刚度按C30混凝土模量的80%取值,即为2.4×10⁷ kPa,泊松比取为0.2。

（3）被动加固区截面尺寸为5 m×5m,沿着基坑纵向满长布置,材料为水泥土,莫尔-库仑模型,刚度为1.5×10⁵ kPa,泊松比为0.25,抗剪强度为750 kPa。

（4）工程桩直径为600 mm,桩长为15 m,桩间距为8 m,桩的刚度为3×10⁷ kPa,桩的重度为7 kN/m³。

2.3 边界条件

对称侧面约束法向自由度及垂直于另外两个方向的转动自由度,非对称侧面仅约束法向自由度,底面约束所有自由度。

2.4 施工工况模拟

数值模拟开挖方案采取在模型竖直方向上开挖一层土,施工一道撑的做法,工况描述见表2。

表 2 计算工况

Table 2. Calculation conditions

施工阶段	工况描述
初始阶段	平衡初始地应力
1	施工围护结构与工程桩
2	放坡开挖至-4 m,降水至-4 m
3	开挖至-5 m,降水至-5 m
4	施工第一道内支撑
5	开挖至-10 m,降水至-10 m
6	施工第二道内支撑
7	开挖至-15 m,降水至-16 m

下载: 导出CSV

| 显示表格

2.5 模型结果验证

为保证有限元数值模拟计算结果的准确可靠,提取了基坑的砼支撑轴力监测数据与基坑围护墙水平位移数据,与PLAXIS模拟结果进行对比,结果见图6。可以看出钢支撑轴力模拟值的变化趋势与实测值的变化趋势一致,开挖至坑底时的误差为7.5%,可以看出钢支撑轴力曲线吻合度较高。

图 6 钢支撑轴力对比图

Figure 6. Axial force contrast diagram of steel bracing

下载: 全尺寸图片幻灯片

3. 工程桩对基坑回弹变形的影响研究

3.1 桩长对基坑回弹变形的影响分析

研究工程桩桩长对基坑回弹变形的影响,保持工程桩桩径（d=0.6 m）、刚度（E=3×10⁷ kPa）不变,分别取工程桩桩长为5,10,15,20,25 m以及无工程桩共6种条件。

由图7看出,无桩的回弹变形曲线与有桩的回弹变形曲线存在明显区别。无桩的回弹变形曲线呈现较为平滑的“凸”形,最大回弹变形量为33 mm,当桩即不施工工程桩时,基坑回弹增长较快,当距坑边约10 m,基坑回弹量不再增长,保持在约33 mm左右。

图 7 桩长-回弹曲线图

Figure 7. Pile length - rebound curves

下载: 全尺寸图片幻灯片

当工程桩存在时,有桩的回弹变形曲线呈现波浪型,波谷处即为打桩的位置,此时,基坑回弹量增长较慢,在距坑边7 m左右就不再增长,最大值为27 mm左右；并且施工工程桩的位置,基坑回弹出现显著的减小。桩长每增加5 m,基坑中心回弹量减小约8%。

3.2 桩径对基坑回弹变形的影响分析

研究工程桩桩径对基坑回弹变形的影响,保持工程桩桩长（L=25 m）、刚度（E=3×10⁷ kPa）不变,分别取工程桩桩径为d=0.6 m,L=0.8 m,d=1.0 m,d=1.2 m,d=1.4 m以及无工程桩共6种条件。

由图8可以看出,各个桩径工况下,随着桩径的增加,基坑回弹变形量逐渐减小,工程桩桩径每增加0.2 m,基坑中心的回弹量减小4%。

图 8 桩径-回弹曲线图

Figure 8. Pile diameter - rebound curve

下载: 全尺寸图片幻灯片

3.3 桩刚度对基坑回弹变形的影响分析

研究工程桩刚度对基坑回弹变形的影响,保持工程桩桩长（L=25 m）、桩径（d=0.6 m）不变,分别取工程桩刚度为E=2.5×10⁷ kPa,E=3.0×10⁷ kPa,E=3.5×10⁷ kPa,E=3.8×10⁷ kPa以及无工程桩5种工况。

由图9可以看出,改变工程桩桩刚度后,基坑回弹变形并未产生明显变化,结果表明,坑底工程桩的桩刚度对基坑回弹变形的影响最小。

图 9 桩刚度-回弹曲线

Figure 9. Pile stiffness-rebound curves

下载: 全尺寸图片幻灯片

4. 结论

（1）基坑开挖初期,开挖深度较浅,基坑回弹变形为弹性变形,呈现四周小中间大的特点；随着开挖深度的增加,侧向卸荷逐渐增大,基坑内部土体产生塑性变形,基坑回弹变形呈现四周小中间大的特点。

（2）由基坑回弹变形曲线可以看出,工程桩附近土体回弹量明显小于相邻土体,最终开挖完成时,有桩较无桩可以减小约20%的回弹变形。

（3）结合基坑回弹变形曲线可以看出,桩长每增加5 m,基坑回弹变形减小约5%左右；桩径每增加0.2 m,基坑回弹变形减小约2%左右；工程桩桩身刚度的变化对基坑回弹变形的影响不是特别明显。

图 1 “锅盖效应”的发生机理

Figure 1. Mechanism of "pot cover effect"

下载: 全尺寸图片幻灯片

图 2 “锅盖效应”的试验原理

Figure 2. Test principle of "pot cover effect"

下载: 全尺寸图片幻灯片

图 3 迁移通量J与温度T的关系曲线

Figure 3. Relationship between diffusion flux J and temperature T

下载: 全尺寸图片幻灯片

图 4 迁移通量J与饱和度S_r的关系曲线

Figure 4. Relationship between diffusion flux J and soil saturation S_r

下载: 全尺寸图片幻灯片

参考文献(26)

[1]	李强, 姚仰平, 韩黎明, 等. 土体的"锅盖效应"[J]. 工业建筑, 2014, 44(2): 69–71. https://www.cnki.com.cn/Article/CJFDTOTAL-GYJZ201402018.htm
[2]	LUO T, QU X, WANG N D, et al. Pot cover effect and its prevention: an experimental study in the field. Cold Regions Science and Technology, 2019, 167: 102845. doi: 10.1016/j.coldregions.2019.102845
[3]	华遵孟, 沈秋武. 敦煌机场跑道破坏原因与防治[J]. 西部探矿工程, 1994, 6(4): 16–19. https://www.cnki.com.cn/Article/CJFDTOTAL-XBTK404.009.htm
[4]	EIGENBROD K D, KENNEPOHL G J A. Moisture accumulation and pore water pressures at base of pavements[J]. Transportation Research Record: Journal of the Transportation Research Board, 1996, 1546(1): 151–161. doi: 10.1177/0361198196154600117
[5]	罗汀, 陈含, 姚仰平, 等. 锅盖效应水分迁移规律分析[J]. 工业建筑, 2016, 46(9): 6–9. doi: 10.13204/j.gyjz201609002
[6]	姚仰平, 王琳. 影响锅盖效应因素的研究[J]. 岩土工程学报, 2018, 40(8): 1373–1382. doi: 10.11779/CJGE201808002
[7]	宋二祥, 罗爽, 孔郁斐, 等. 路基土体锅盖效应的数值模拟分析[J]. 岩土力学, 2017, 38(6): 1781–1788. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201706028.htm
[8]	张升, 贺佐跃, 滕继东, 等. 非饱和土水汽迁移与相变: 两类"锅盖效应"的试验研究[J]. 岩土工程学报, 2017, 39(5): 961–968. doi: 10.11779/CJGE201705022
[9]	滕继东, 贺佐跃, 张升, 等. 非饱和土水气迁移与相变: 两类"锅盖效应"的发生机理及数值再现[J]. 岩土工程学报, 2016, 38(10): 1813–1821. doi: 10.11779/CJGE201610010
[10]	罗汀, 曲啸, 姚仰平, 等. 北京新机场"锅盖效应"一维现场试验[J]. 土木工程学报, 2019, 52(增刊1): 233–239. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2019S1030.htm
[11]	MILLY P C D. A simulation analysis of thermal effects on evaporation from soil[J]. Water Resources Research. 1984, 20(8): 1087–1098. doi: 10.1029/WR020i008p01087
[12]	MIYAZAKI T. Condensation and movement of water vapor in sand under temperature gradient[J]. Trans Jpn Soc Irrig Drain Rural Eng, 1976, 61: 1–8.
[13]	CAHILL A T, PARLANGE M B. On water vapor transport in field soils[J]. Water Resources Research, 1998, 34(4): 731–739. doi: 10.1029/97WR03756
[14]	ZHANG Y Z, MA W, WANG T L. Characteristics of the liquid and vapor migration of coarse-grained soil in an open-system under constant-temperature freezing[J]. Cold Regions Science and Technology, 2019, 165: 102793. doi: 10.1016/j.coldregions.2019.102793
[15]	SHERWOOD T K. The drying of soil-II[J]. Industrial & Engineering Chemistry, 1929, 21(10): 976–980.
[16]	CEAGLSKE N H, HOUGEN O A. Drying granular solids[J]. Industrial & Engineering Chemistry, 1937, 29(7): 805–813.
[17]	HENRY P S A. Diffusion in absorbing media[J]. Proc R Soc London, 1939, A171: 215–241.
[18]	PHILIP J R, DE VRIES D A. Moisture movement in porous materials under temperature gradients[J]. Transacations, American Geophysical Union, 1957, 38(2): 222–231.
[19]	DE VRIES D A. Simultaneous transfer of heat and moisture in porous media[J]. Transactions, American Geophysical Union, 1958, 39(5): 909.
[20]	FREDLUND D G, RAHARDJO H. Soil Mechanics for Unsaturated Soils[M]. New York: Wiley, 1993.
[21]	EWEN J, THOMAS H R. Heating unsaturated medium sand[J]. Geotechnique. 1989, 39(3): 455–470.
[22]	陈含. 冻融循环条件下锅盖效应机理及试验研究[D]. 北京: 北京航空航天大学, 2020.
[23]	PARTINGTON J R. An Advanced Treatise on Physical Chemistry[M]. London: Longmans, Green, 1949.
[24]	KOPONEN A, KATAJA M, TIMONEN J. Tortuous flow in porous media[J]. Physical Review E, 1996, 54(1): 406–410.
[25]	KRISCHER D. Wärme- und Stofftransport in Feuchten Gütern[M]. Trocknungstechnik. Berlin, Heidelberg: Springer Berlin Heidelberg, 1978: 269–286.
[26]	司娟娟. 不同级配下砾石土土水特征曲线研究[J]. 低碳世界, 2016(15): 248–249. https://www.cnki.com.cn/Article/CJFDTOTAL-DTSJ201615151.htm