土体硬化模型参数试验研究及其在南昌地区基坑工程的应用

2021, Vol. 43

Issue (6): 38-47 doi: 10.11835/j.issn.2096-6717.2020.192

PDF

Article Options
	PDF

	Abstract

	Figures

	References


扩展功能

	把本文推荐给朋友



	Email Alert

	RSS

	RIS（文献管理工具）
本文作者相关文章
	刘伟煌

	朱怀龙

	贺斯进

	言建标

	徐长节

土体硬化模型参数试验研究及其在南昌地区基坑工程的应用

刘伟煌 ¹, 朱怀龙 ², 贺斯进 ³, 言建标 ⁴, 徐长节 ^1,5

1. 华东交通大学江西省岩土工程基础设施安全与控制重点实验室, 南昌 330013;
2. 江西交通职业技术学院, 南昌, 330013;
3. 南昌轨道交通集团有限公司, 南昌, 330038;
4. 浙江杭海城际铁路有限公司, 杭州, 314000;
5. 浙江大学滨海和城市岩土工程研究中心, 杭州 310058

收稿日期：2020-08-04

基金项目：国家杰出青年科学基金（51725802）；国家自然科学基金（51878276）

作者简介：刘伟煌(1995-), 男, 主要从事基坑工程研究, E-mail: 824098351@qq.com.

通信作者：贺斯进(通信作者), 男, 高级工程师, E-mail: 1697296226@qq.com.

摘要：基坑开挖数值分析的关键在于土体本构的选择和计算参数的合理选取，硬化模型（HS模型）是应用较为广泛的模型之一。通过采用GDS空心扭剪试验仪（SS-HCA）和常规固结仪完成大量室内试验，获取了南昌地区软土层的硬化模型参数，对各层土参数中模量之间的比例关系进行探讨，并将试验得到的结果与相关文献的结果进行比较和总结。运用PLAXIS有限元数值计算软件，对南昌某地下室深基坑工程进行数值模拟，土体本构模型采用程序内置的HS模型。结果表明，围护桩深层水平位移的实测数据与数值计算结果基本吻合，表明获取的HS模型参数及方法适用于南昌地区基坑开挖的数值分析。

关键词：土体硬化模型基坑工程空心扭剪试验仪室内试验

Experimental study on parameters of hardening soil model for soils and its application in foundation pit engineering in Nanchang

LIU Weihuang ¹, ZHU Huailong ², HE Sijin ³, YAN Jianbiao ⁴, XU Changjie ^1,5

1. Jiangxi Key Laboratory of Infrastructure Safety Control in Geotechnical Engineering, East China Jiaotong University, Nanchang 330013, P. R. China;
2. Jiangxi V & T College of Communications, Nanchang 330013;
3. Nanchang Rail Transit Co., Ltd., Nanchang 330038, P. R. China;
4. Zhejiang Hanghai Intercity Railway Co., Ltd., Hangzhou 314000, P. R. China;
5. Research Center of Coastal and Urban Geotechnical Engineering, Zhejiang University, Hangzhou 310058, P. R. China

Received: 2020-08-04

Foundation item: National Science Fund for Distinguished Young Scholars (No. 51725802); National Natural Science Foundation of China (No. 51878276)

Author brief: LIU Weihuang (1995-), main research interest: foundation pit engineering, E-mail: 824098351@qq.com.

Corresponding author: HE Sijin (corresponding author), senior engineer, E-mail: 1697296226@qq.com.

Abstract: The significance of numerical analysis of foundation pit excavation lies in the selection of soil constitution and reasonable input parameters.The hardening model(HS model) is one of the soil constitutive models that are widely used.In this paper, the GDS Hollow Torsion Shear Tester (SS-HCA) and consolidation apparatus were used to complete a series of laboratory tests.The parameters of the hardening model of the soft soil layer in Nanchang were obtained through the tests. The proportional relationship between the modulus of each layer of soil parameters has been discussed.The experimental results were compared with the result from literatures. Then, the PLAXIS finite element software was used to conduct numerical simulation of a certain basement deep foundation pit project in Nanchang. The soil constitutive model adopted the programmed HS model. The results show that the calculated value of the deep horizontal displacement of the retaining pile is basically consistent with the measured value. The HS model and the determination method for parameters have been proved to be applicable in numerical analysis of foundation pit excavation engineering in Nanchang area.

Keywords: hardening soil model foundation pit engineering hollow torsional shear tester laboratory test

随着经济的发展，南昌地区的基坑数量也随之增多，基坑规模也越来越大，加之城市用地日益紧张，基坑周边常遇到地下管线、地铁隧道、高架桥等建筑。在这种复杂情况下，基坑工程除满足自身的变形要求外，还需要保证其周边环境的安全性，这对基坑工程的设计提出了更高的要求。目前，分析基坑开挖对周边环境影响的方法有解析法和相关规范法^[1-2]。但考虑问题的复杂性，传统的解析法和规范法难以真实地反映基坑开挖期间对周边环境的影响。随着有限元和计算机的发展，数值分析法^[3-5]越来越成为深基坑变形分析中常用且有效的方法，而数值分析法的关键问题之一在于采用合适的土体本构模型和选取合理的计算参数^[6]。

土体硬化本构模型(HS模型)适用于基坑开挖的数值分析，它能得到较合理的围护桩(墙)位移、支护结构内力和坑外地表沉降^[6]。但由于HS模型包含的计算参数较多，要获得完整的HS模型参数难度很大。一些学者已对获取HS模型参数进行了一些研究，Calvello等^[7]通过反演分析法得到了芝加哥地区软土层的HS模型参数，并探讨了该方法影响其结果准确性的因素。李连祥等^[8]运用位移反分析方法获得济南典型土层的HS模型参数。刘畅^[9]通过室内三轴试验对天津地区某深基坑软土的硬化模型参数进行测定，获取了其部分模量参数。梁发云等^[10]通过大量的室内试验得到上海地区较为完整的土体HS模型参数。陈尚荣等^[11]通过室内试验获得了上海临港粉质砂土的HS模型参数。由于不同地区土体的硬化模型参数值存在一定差异，直接运用已有研究成果进行数值分析会对结果的准确性产生一定偏差。目前还没有对南昌地区土体硬化模型参数获取的相关报道。为此，针对南昌地区软土进行了室内土工试验，以获取其硬化模型参数c′、φ′、R_f、E₅₀^ref、E_oed^ref和E_ur^ref，并对获取的HS模型参数各种模量之间的比例关系进行探讨。随后，运用PLAXIS有限元数值计算软件对南昌某地下室深基坑工程进行数值模拟，将围护桩水平位移的计算结果与实测结果进行对比，从而验证获取的HS模型参数及方法在南昌地区基坑开挖工程中的适用性。

1 HS模型及参数介绍

HS模型，即土体硬化本构模型，是由Schanz等^[12]提出来的一种以塑性理论为基础的双屈服面模型，塑性部分采用了各向同性硬化准则。HS模型共包含11个参数：有效黏聚力c′、有效内摩擦角φ′、刚度应力水平相关幂指数m、参考应力p^ref、三轴排水剪切试验的参考加、卸载模量E_ur^ref、固结试验中的参考割线模量E_oed^ref、三轴排水剪切试验的参考割线模量E₅₀^ref、破坏比R_f、静止侧压力系数K₀、土的剪胀角ψ、加卸载泊松比v_ur。其中，参数K₀、ψ、p^ref、v_ur、m可参照已有的研究成果来确定，本文主要通过室内试验确定其他6个参数。

2 试验过程

2.1 现场取样

试验共采集了南昌地区某深基坑工地的4种不同的土样，每层土样取土深度也不相同，其中，各层土样的取土深度为3.0、6.2、7.5、10.1 m，对应的土层分别为黏土层②、淤泥质粉质黏土层③、粉质黏土层④和砂质粉土层⑤。其基本物理力学参数指标如表 1所示。

表 1 土层物理力学参数 Table 1 Mechanical parameters of soils

2.2 试验方法及试验设备

针对4种不同的土样，分别通过3种不同类型的室内试验来测定其参数c′、φ′、R_f、E₅₀^ref、E_oed^ref和E_ur^ref。1)通过常规三轴固结排水剪切试验获得参数c′、φ′、R_f、E₅₀^ref；2)通过三轴固结排水加、卸载试验获得加载模量参数E_ur^ref；3)通过标准固结试验获得切线模量参数E_oed^ref。

根据3种不同的试验分别采用不同的试验设备，其中，三轴固结排水剪切试验和三轴固结排水加、卸载试验采用的是GDS空心扭剪试验系统(SS-HCA)，见图 1，它包含常规GDS系统三轴仪的所有功能，由于GDS空心扭剪试验仪的试样为空心圆柱试样；而三轴试验试样为圆柱试样，因此，试验前需要将其底座换成常规GDS系统三轴试验的圆柱底座，见图 2。标准固结试验所采用的仪器为常规固结仪。

图 1 GDS空心扭剪试验仪器图 Fig. 1 Pictuers of hollow torsional shear tester GDS

图 2 装样过程示意图 Fig. 2 Diagrams of soil sample loading process

2.3 试验方案及步骤

2.3.1 常规三轴固结排水剪切试验

为了获得土样②~⑤的HS模型参数E₅₀^ref、R_f、c′和φ，进行了三轴固结排水剪切试验，在试验中，需对4种不同土样设置3种不同的围压，其中，考虑到土层②埋深较浅，因此，取70、100、200 kPa 3种不同围压，其余土样(③~⑤)所取的围压为100、200、300 kPa。

常规三轴固结排水剪切试验包括以下步骤：1)试样制备：将现场取回的原状土用切土盘制备成圆柱试样，试样尺寸为100 mm×50 mm(高度×直径)，然后用游标卡尺量取试样的实际直径和高度，并称其质量；2)预饱和：为了使试样充分饱和并节省时间，将制备好的试样放入饱和缸进行真空抽气饱和，饱和时间为24 h；3)反压饱和：压力室充满水后，打开围压、反压以及孔压阀门，围压设置为110 kPa, 反压设置为100 kPa，时间设置为3 h。4)B值检测：设置围压140 kPa，反压不变并保持关闭，测孔隙水压力系数B，若B=Δu/Δσ＞95%，认为试样达到饱和。5)固结：打开反压阀门，设置围压σ₃，对试样进行等向固结，时间设置为48 h。6)排水剪切：排水剪切采用等应变速率控制，排水剪切速率设置为0.005 2 mm/min。当试样应变值达到其最大应变值的20%时，停止试验。

2.3.2 三轴固结排水加载-卸载-再加载试验

三轴固结排水加、卸载试验包含7个步骤，其中，步骤1)~6)与三轴固结排水剪切试验相同，步骤7)进行轴向加载-卸载-再加载试验。选择GDS-ttAdvanced Loading试验模块，采用轴向应力进行控制，将围压设置为100 kPa，打开反压阀门设置反压为0。首先，轴向加载到试样预计破坏偏应力的40%，当达到目标值时，立即轴向卸载到0，然后再进行轴向加载到试样预计破坏偏应力的60%。

2.3.3 标准固结试验

为获取固结试验中的参考割线模量E_oed^ref，进行了标准固结试验，试验给土样②~⑤施加5种不同等级的荷载(50、100、200、400、800 kPa)。试验步骤依据土工试验标准进行。

3 试验结果分析

3.1 三轴固结排水剪切试验

图 3为土样②~⑤偏应力q(q=σ₁-σ₃)与轴向应变的关系曲线，其参考围压均为σ₃=100 kPa。

图 3 ②~⑤层土样三轴试验应力应变曲线 Fig. 3 Strain-stress curves of triaxial CD test on soil layers(②~⑤)

由图 3可知，轴向应变在0~10%时，各层土的偏应力随着轴向应变的增加而增大，当轴向应变达到一定范围时，偏应力也逐渐保持平缓或者稍微有所下降。当轴向应变达到15%时，曲线所对应的纵坐标(偏应力值)即为破坏值q_f^[13]，由此可得，②~⑤层土样的q_f分别为220.23、240.95、258.21、308.51 kPa。连接曲线的原点和0.5q_f所对应的点，其直线斜率即为E₅₀^ref，因此，可以得到围压为100 kPa时②~⑤层土样的参数E₅₀^ref分别为4.8、2.8、4.9、9.1 MPa。另外，取曲线稳定段对应的偏应力值作为渐近线值q_a，通过计算可得②~⑤层土样破坏比R_f(R_f=q_f/q_a)分别为0.95、0.90、0.92和0.96。将上述参数值E₅₀^ref和R_f汇总于表 2。

表 2 ②~⑤层土体E₅₀^ref和R_f值 Table 2 Values of E₅₀^ref and R_f of soil layers(②~⑤)

为获取各层试样的有效黏聚力c′和内摩擦角φ′，将各层土样分别进行其他两种围压的三轴固结排水剪切试验，图 4为各层土样的摩尔应力圆，由图 4可知，②~⑤层试样的摩尔圆均相切于同一条直线，由此可得各层土的有效黏聚力c′和内摩擦角φ′，并将其汇总于表 3。

图 4 ②~⑤层土体摩尔应力圆曲线 Fig. 4 Mohr circle of soil layers (②~⑤)

表 3 ②~⑤层土体c′和φ′值 Table 3 Values of c′ and φ′ of soil layers(②~⑤)

3.2 三轴固结排水加载-卸载-再加载试验结果

图 5为各层土样加、卸载试验的偏应力与轴向应变关系曲线。

图 5 ②~⑤层土三轴加、卸载试验应力-应变曲线 Fig. 5 Strain-stress curves of triaxial CD loading unloading-reloading rests on soil layers (②~⑤)

由图 5可知，当初次加载时，轴向应变随着荷载的增大而逐渐增加，卸载时，偏应力-应变曲线陡降，并且轴向应变稍有减小，再次加载，曲线开始变陡而随后逐渐变缓。各层土样曲线在加载-卸载-再加载试验过程中均表现为一个滞回圈。由参考文献[14]可得，将曲线滞回圈的两端点进行连接，该直线斜率表示围压100 kPa下试样的参考卸载再加载模量E_ur^ref。综上所述，②~⑤土层的E_ur^ref值分别为38.6、36.2、39.3、59.1 MPa。

3.3 标准固结试验结果

图 6为各土层在不同轴向荷载下与轴向应变关系曲线，图中的曲线通过拟合函数所获得，②~⑤层土样轴向荷载与轴向应变曲线拟合函数分别为p=3.017 7ε_a²+6.721 9ε_a、p=1.043 1ε_a²+5.631 3ε_a、p=4.937 6ε_a²+19.914 6ε_a、p=-13.448 8(1-1.656 5ε_a)，其中，每层土样曲线的拟合函数R²均为0.99。由图 6可知，各层土样的轴向荷载与轴向应变曲线变化规律基本相同，当轴向荷载较小时，曲线变化较为平缓，随着轴向荷载的逐级加大，曲线在逐渐变陡，曲线的斜率也随之增大。参考割线模量E_oed^ref为轴向荷载为100 kPa时曲线的斜率，由此可根据其曲线的拟合函数获得E_oed^ref，通过计算可得，各层土样的参考割线模量E_oed^ref分别为3.8、2.3、5.2、8.0 MPa。

图 6 各土层固结试验的轴向荷载与应变关系曲线图 Fig. 6 Relationships between load and strain of consolidation test

为了得到各层土的压缩模量E_s1-2，将固结试验的孔隙比与轴向荷载的关系曲线绘制于图 7，由图 7可知，各层土的孔隙比随着轴向荷载的增大而逐渐减小，当轴向荷载达到800 kPa时，各土层的孔隙比趋于平缓，通过计算可得，土层(②~⑤)的压缩模量E_s1-2分别为6.1、2.5、6.8、10.6 MPa。

图 7 各土层固结试验的孔隙比与轴向荷载关系曲线图 Fig. 7 Relationships between load and void ration of consolidation test

4 试验结果对比分析

通过3种不同的室内试验，获得南昌地区软土层(②~⑤)的HS模型参数c′、φ′、E_oed^ref、E₅₀^ref、E_ur^ref和R_f。将上述参数的取值汇总于表 4。从表 4可知，南昌地区土层③的c′、φ′、E_oed^ref、E₅₀^ref、E_ur^ref和R_f均小于土层②和④。砂质粉土层⑤的c′值均小于土层②~④，φ′值均大于②~④层土样。②~⑤层土样的E₅₀^ref和E_oed^ref之间的比例关系基本相同(E₅₀^ref=0.94~1.3E_oed^ref)；对于E_ur^ref和E₅₀^ref两者之间的比例关系，②~⑤层土体较为接近(E_ur^ref=6.5~12.9E_oed^ref)，南昌地区砂质粉土层⑤的E₅₀^ref值是E_oed^ref的1.2倍，E_ur^ref值是E₅₀^ref的6.5倍；而对于E_oed^ref与E_s1-2两者之间关系，土层②~⑤的比例关系大致相同(E_ur^ref=0.62~0.76E_s1-2)。

表 4 土体HS模型参数对比表 Table 4 Comparisions of parameters of the HS model for soils

土体名称		c′/kPa	φ′(°)	E_oed^ref/MPa	E₅₀^ref/MPa	E_ur^ref/MPa	R_f	E_s1-2/MPa	模量比例关系
南昌地区(本次试验)	②黏土	8.6	29.2	3.8	4.8	38.6	0.95	6.1	E₅₀^ref=1.3 E_oed^ref E_ur^ref=8.0 E₅₀^ref E_oed^ref=0.62 E_s1-2
	③淤泥质粉质黏土	6.5	25.6	2.3	2.8	36.2	0.90	2.5	E₅₀^ref=1.2 E_oed^ref E_ur^ref=12.9 E₅₀^ref E_oed^ref=0.74 E_s1-2
	④粉质黏土	25.8	28.2	5.2	4.9	39.3	0.92	6.8	E₅₀^ref=0.9 E_oed^ref E_ur^ref=8.0 E₅₀^ref E_oed^ref=0.76 E_s1-2
	⑤砂质粉土	2.0	36.8	8.0	9.1	59.1	0.96	10.6	E₅₀^ref=1.2 E_oed^ref E_ur^ref=6.5 E_oed^ref E_oed^ref=0.75 E_s1-2
上海地区^[10]	②黏土	10.37	29	3.7	4.7	38.8	0.91	5.9	E₅₀^ref=1.2 E_oed^ref E_ur^ref=8.4 E_oed^ref E_oed^ref=0.36 E_s1-2
	③淤泥质粉质黏土	8.13	34	2.5	3.0	33.9	0.68	2.3	E₅₀^ref=1.1 E_oed^ref E_ur^ref=11.6 E_oed^ref E_oed^ref=1.06 E_s1-2
	④淤泥质黏土	4.96	27	3.6	3.9	36.6	0.72	4.3	E₅₀^ref=1.08 E_oed^ref E_ur^ref=9.4 E_oed^ref E_oed^ref=0.85 E_s1-2
	⑤粉质黏土	16.90	34	5.9	6.4	42.4	0.89	6.7	E₅₀^ref=1.08 E_oed^ref E_ur^ref=6.7 E_oed^ref E_oed^ref=0.87 E_s1-2
天津滨海软土地区^[9]	②粉质黏土			3.8	2.8			1.8	E₅₀^ref=0.7 E_oed^ref E_ur^ref=2.1 E₅₀^ref
	③淤泥质黏土			3.5	2.3			1.8	E₅₀^ref=0.7 E_oed^ref E_ur^ref=2.1 E₅₀^ref
	④黏土			5.0	5.3			2.8	E₅₀^ref=1.1 E_oed^ref E_ur^ref=1.8 E₅₀^ref
	⑤粉质黏土			5.6	2.9			4.0	E₅₀^ref=0.5 E_oed^ref E_ur^ref=1.4 E₅₀^ref
无锡地区^[15]	③粉质黏土			4.8	8.0	25.2			E₅₀^ref=0.7 E_oed^ref E_ur^ref=3.2 E₅₀^ref
北京地区^[16]	粉质黏土			4.0	3.9	17.8	0.8	3.6	E₅₀^ref= E_oed^ref E_ur^ref=4.5 E₅₀^ref E_oed^ref=1.1 E_s1-2
LacustrineClay^[17]				37.6	37.6		150.4		E₅₀^ref= E_oed^ref E_ur^ref=4 E₅₀^ref
UpperBlodgett^[18]				1.6	2.4		10	0.7	E₅₀^ref= E_oed^ref E_ur^ref=4.3 E₅₀^ref

表 4 土体HS模型参数对比表 Table 4 Comparisions of parameters of the HS model for soils

为了对比分析，表 4中还列出了其他地区土体的参数。南昌地区各土层E₅₀^ref值为E_oed^ref的0.9~1.3倍，这与上海地区^[10]的1.02~1.2倍、天津滨海软土地区的0.5~1.1倍、无锡地区粉质黏土的1.7倍、北京地区粉质黏土的1倍、奥地利Lacustrine Clay^[17]的1倍以及美国Upper Blodgett^[18]的1.5倍比例比较接近。南昌地区各土层E_ur^ref值为E_oed^ref的6.2~12.9倍，这与上海地区^[10]的6.7~11.2倍、Lacustrine Clay^[17]软土地区的4倍以及美国Upper Blodgett^[18]的4.5倍比例比较接近，大于天津滨海地区软土的1.4~2.1倍。

5 工程验证

5.1 工程概况

南昌某地下室深基坑项目位于南昌市西湖区孺子路与船山路交汇处，整个基坑平面呈“L”形，场地北侧为孺子路，西侧为船山路，南面及东面为居民区，基坑平面位置如图 8所示。基坑包括B01和B04两地块，并将其两个地块整体设计，基坑周长约470 m²，面积约5 900 m²。B01和B04地块基坑开挖顺序不同，先施工B01地块，后施工B04地块，开挖深度分别约为11.45、10.85 m，工程采用钻孔灌注桩作为围护结构，B01和B04地块的围护桩桩长分别为19.45、21.85 m，桩径1 m，中心距1.2 m，在围护桩桩顶布置冠梁，冠梁截面尺寸为1 200 mm×800 mm(长×宽)；在B01地块竖直方向设置两道钢筋混凝土支撑，分别架设于桩顶以及距桩顶3.0 m位置处，在B04地块桩顶设置一道钢筋混凝土内支撑，内支撑横截面均为1 000 mm×800 mm(长×宽)，图 9为支护体系剖面图。

图 8 基坑平面位置示意图 Fig. 8 Pictures of plane position of foundation pit

图 9 B04地块支护体系剖面图(单位：mm) Fig. 9 Pictures of profile of area B04 support system(unit: mm)

5.2 计算模型与参数

5.2.1 计算模型

采用PLAXIS 3D有限元软件对该工程进行数值模拟，有限元模型尺寸为400 m×300 m×50 m(长×宽×深)，土体采用10节点四面体实体单元模拟，本构模型采用软件程序自带的HS本构，内支撑结构、冠梁和立柱桩均采用梁单元模拟。为了计算方便，围护结构采用6节点板单元模拟，并将排桩等效为地下连续墙，等效厚度采用文献[19]公式计算，此外，采用12节点界面单元模拟土和结构的相互作用。土体四周约束其法向位移，底部约束其竖向及水平位移。为了既能节约模型计算成本又能得到较为精确的计算结果，模型网格划分精度选择中等粗细程度，并且对围护墙、支撑以及立柱桩附近的网格进行局部加密，整个计算模型共包含15 462个计算单元, 计算模型网格划分见图 10。支撑、围护桩、立柱桩的物理力学参数见表 5。

图 10 有限元网格 Fig. 10 Finite element mesh

表 5 内支撑、围护桩、立柱桩计算参数 Table 5 Parameters for the support, retaining pile and column pile

5.2.2 土体参数的选取

模型中土层(②~⑥)采用HS模型，土层(②~⑤)参数c′、φ′、E_oed^ref、E₅₀^ref、E_ur^ref和R_f根据试验结果取值，其他参数根据表 6中已有文献成果取值。土层⑥为圆砾层，其基本参数可根据工程勘察报告取得：c′=0、φ′=37.3°，重度γ=20.5 kN/m³，泊松比υ=0.2, 压缩模量E_s1-2=35 MPa。对于南昌地区圆砾层，目前还没有对E_oed^ref、E₅₀^ref和E_ur^ref研究的相关报道，其参数取值可参考文献[20]获得，由此可得本工程场地②~⑥土层HS模型参数，具体见表 7。土层⑦为泥质粉砂岩层，土体本构采用摩尔库伦模型(MC模型)，基本参数：c′=35 kPa，φ′=200°，泊松比υ=0.2，重度γ=23.5 kN/m³, 弹性模量E=0.36 MPa，综上，可得模型中全部土层参数。

表 6 HS模型部分参数取值 Table 6 Evaluation of partial parameters of HS model

表 7 土体(②~⑥)HS模型参数指标 Table 7 Parameters of HS model of soil layers

土层	平均厚度/m	c′/ kPa	φ′/(°)	γ/(kN·m^-3)	K₀	ψ	m	p^ref/MPa	υ_ur	E₅₀^ref/MPa	E_ur^ref/MPa	E_oed^ref/MPa	R_f
②	2.7	8.6	29.2	18.00	0.51	0	0.8	100	0.2	4.8	38.6	3.8	0.95
③	2.6	6.5	25.6	17.10	0.57	0	0.8	100	0.2	3.8	34.2	2.3	0.90
④	1.5	25.8	28.2	20.30	0.53	0	0.8	100	0.2	4.9	39.3	5.2	0.92
⑤	2.1	2.0	36.8	20.50	0.40	0	0.8	100	0.2	9.1	59.1	8.0	0.96
⑥	4.3	0	37.3	20.50	0.39	7.3	0.8	100	0.2	26.3	131.5	26.3	0.90

表 7 土体(②~⑥)HS模型参数指标 Table 7 Parameters of HS model of soil layers

5.3 计算结果分析

为了监测基坑开挖期间基坑围护结构的变形，保证施工过程中基坑的安全性，在围护桩内埋设了多个测点，如图 8所示，现取AB和CD两段，其对应的两个测点分别为CX6和CX3，分析基坑开挖过程中围护桩沿深度方向的变形规律。图 11为基坑开挖期间围护桩水平位移实测值与计算值的对比曲线。由图 11可知，围护桩的深层水平位移曲线的分布状态随着基坑的逐次开挖在发生变化，并且在基坑开挖面以上的水平位移值逐渐增大，而基坑开挖面以下的水平位移值变化较小，当开挖至基坑底时，桩体的水平位移达到最大，且最大值出现在桩顶。桩体测点CX3和CX6最大水平位移的实测值分别为16.09、17.21 mm，数值模拟的最大水平位移值分别为17.08、17.91 mm，测点CX3和CX6的实测数据与数值模拟结果变形规律基本一致，其中，对于测点CX6，基坑第3次开挖至坑底时计算值稍大于实测值，总体而言，吻合较好。

图 11 围护桩水平位移对比曲线图 Fig. 11 Comparisions of horizontal displacement of piles

6 结论

针对南昌地区软土进行了室内土工试验，获得了土体的HS模型参数值，并将所得到的参数运用到南昌地区的基坑工程的数值分析中，得到如下结论：

1) 通过室内试验获得了南昌地区土体HS模型的主要参数c′、φ′、E_oed^ref、E₅₀^ref、E_ur^ref和R_f的取值及其模量之间的比例关系。

2) 南昌地区③层土样的c′、φ′、E_oed^ref、E₅₀^ref、E_ur^ref和R_f均小于②和④层土样。②~⑤层土样的E₅₀^ref和E_oed^ref之间的比例关系基本相同(E₅₀^ref=0.9E_oed^ref~1.3E_oed^ref)，E_ur^ref和E₅₀^ref两者之间的比例关系较为接近(E_ur^ref=6.5E_oed^ref ~12.9E_oed^ref)。

3) 南昌地区砂质粉土层⑤的c′值小于②~④层土样，φ′值大于②~④层土样，E₅₀^ref值是E_oed^ref值的1.2倍，E_ur^ref值E₅₀^ref值是的6.5倍，E_oed^ref是E_s-2的0.75倍。

4) 运用PLAXIS有限元数值软件对南昌某地下室深基坑工程进行三维数值分析，有限元数值计算结果与实测数据吻合较好，验证了获取的HS模型参数及方法适用于南昌地区基坑开挖的数值分析。

参考文献

[1]	陈晓平, 闫军. 深基坑支护结构的三维杆系有限元分析[J]. 岩土力学, 2001, 22(3): 258-261. CHEN X P, YAN J. 3D pole system FEM analysis for bracing structure of deep foundation pit[J]. Rock and Soil Mechanics, 2001, 22(3): 258-261. (in Chinese) DOI:10.3969/j.issn.1000-7598.2001.03.005
[2]	建筑基坑支护技术规程: JGJ 120—1999[S]. 北京: 中国标准出版社, 1999. Technical specification for retaining and protection of building foundation excavations: JGJ 120-1999[S]. Beijing: Standards Press of China, 1999. (in Chinese)
[3]	胡海英, 张玉成, 杨光华, 等. 基坑开挖对既有地铁隧道影响的实测及数值分析[J]. 岩土工程学报, 2014, 36(Sup2): 431-439. HU H Y, ZHANG Y C, YANG G H, et al. Measurement and numerical analysis of effect of excavation of foundation pits on metro tunnels[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(Sup2): 431-439. (in Chinese)
[4]	林杭, 陈靖宇, 郭春, 等. 基坑开挖对邻近既有隧道变形影响范围的数值分析[J]. 中南大学学报(自然科学版), 2015, 46(11): 4240-4247. LIN H, CHEN J Y, GUO C, et al. Numerical analysis on influence of foundation pit excavation on deformation of adjacent existing tunnel[J]. Journal of Central South University (Science and Technology), 2015, 46(11): 4240-4247. (in Chinese) DOI:10.11817/j.issn.1672-7207.2015.11.035
[5]	林敏博, 刘润, 徐余, 等. 基坑开挖对抗拔桩承载特性影响的数值分析[J]. 岩土工程学报, 2013, 35(Sup2): 735-739. LIN M B, LIU R, XU Y, et al. Numerical analysis of bearing behaviors of uplift piles during pit excavation[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(Sup2): 735-739. (in Chinese)
[6]	宋广, 宋二祥. 基坑开挖数值模拟中土体本构模型的选取[J]. 工程力学, 2014, 31(5): 86-94. SONG G, SONG E X. Selection of soil constitutive models for numerical simulation of foundation pit excavation[J]. Engineering Mechanics, 2014, 31(5): 86-94. (in Chinese)
[7]	CALVELLO M, FINNO R J. Selecting parameters to optimize in model calibration by inverse analysis[J]. Computers and Geotechnics, 2004, 31(5): 410-424. DOI:10.1016/j.compgeo.2004.03.004
[8]	李连祥, 刘嘉典, 李克金, 等. 济南典型地层HSS参数选取及适用性研究[J]. 岩土力学, 2019, 40(10): 4021-4029. LI L X, LIU J D, LI K J, et al. Study of parameters selection and applicability of HSS model in typical stratum of Jinan[J]. Rock and Soil Mechanics, 2019, 40(10): 4021-4029. (in Chinese)
[9]	刘畅. 考虑土体不同强度与变形参数及基坑支护空间影响的基坑支护变形与内力研究[D]. 天津: 天津大学, 2008. LIU C. Analysis of deformation and stress due to deep excavation considering different deformation and strength parameters of soil and space effect of excavation and retaining structure [D]. Tianjin: Tianjin University, 2008. (in Chinese)
[10]	梁发云, 贾亚杰, 丁钰津, 等. 上海地区软土HSS模型参数的试验研究[J]. 岩土工程学报, 2017, 39(2): 269-278. LIANG F Y, JIA Y J, DING Y J, et al. Experimental study on parameters of HSS model for soft soils in Shanghai[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(2): 269-278. (in Chinese)
[11]	陈尚荣, 李通达, 梁发云, 等. 上海临港砂质粉土硬化土小应变模型参数研究[J]. 同济大学学报(自然科学版), 2020, 48(6): 841-846. CHEN S R, LI T D, LIANG F Y, et al. Parameters of hardening soil model with small strain of shallow sandy silt in Shanghai Lin-Gang special area[J]. Journal of Tongji University (Natural Science), 2020, 48(6): 841-846. (in Chinese)
[12]	SCHANZ T, VERMEER P A, BONNIER P G. The hardening soil model: Formulation and verification[M]//Beyond 2000 in Computational Geotechnics-10 years of PLAXIS. Routledge, 2019: 281-296.
[13]	土工试验方法标准: GB/T 50123—1999[S]. 北京: 中国计划出版社, 1999. Standard for soil test method: GB/T 50123-1999[S]. Beijing: China Planning Press, 1999. (in Chinese)
[14]	李广信. 高等土力学[M]. 北京: 清华大学出版社, 1999. LI G X. Advanced Soil Mechanics[M]. Beijing: Tsinghua University Press, 1999. (in Chinese)
[15]	陈峰. 无锡地铁基坑典型地层本构模型适应性研究[D]. 上海: 同济大学, 2011. CHEN F. Adaptability of constitutive model for typical soil layers in Wuxi metro excavation [D]. Shanghai: Tongji University, 2011.
[16]	周恩平. 考虑小应变的硬化土本构模型在基坑变形分析中的应用[D]. 哈尔滨: 哈尔滨工业大学, 2010. ZHOU E P. Application of hardening soil model with small-strain in deformation analysis for foundation pit [D]. Harbin: Harbin Institute of Technology, 2010. (in Chinese)
[17]	LVFTENEGGER R, SCHWEIGER H, SCHARINGER F. 3D finite element analysis of a deep excavation and comparison with in situ measurements[M]//Geotechnical Aspects of Underground Construction in Soft Ground, CRC Press, 2008: 193-199. .
[18]	ROBOSKI J F. Soil parameters for constitutive models of compressible Chicago glacial clays [EB/OL].
[19]	LIAO S M, WEI S F, SHEN S L. Structural responses of existing metro stations to adjacent deep excavations in Suzhou, China[J]. Journal of Performance of Constructed Facilities, 2016, 30(4): 04015089. DOI:10.1061/(ASCE)CF.1943-5509.0000845
[20]	刘志祥, 张海清. PLAXIS高级应用教程[M]. 北京: 机械工业出版社, 2015. LIU Z X, ZHANG H Q. PLAXIS Advanced Application Course[M]. Beijing: China Machine Press, 2015. (in Chinese)
[21]	BRINKGREVE R B J, BROERE W. Plaxis material models manual[M]. Delft, 2006.
[22]	GAO D Z, WEI D D, HU Z X. Geotechnical properties of Shanghai soils and engineering applications[M]//Marine Geotechnology and Nearshore/Offshore Structures. 100 Barr Harbor Drive, West Conshohocken, PA 19428-2959: ASTM International, 1986: 161-178.
[23]	BOLTON M D. Discussion: The strength and dilatancy of sands[J]. Géotechnique, 1987, 37(2): 219-226. DOI:10.1680/geot.1987.37.2.219
[24]	JANBU J. Soil compressibility as determined by oedometer and trixial test [C]//Proceedings of the 3^rd European Conference on Soil Mechanics and Foundation Engineering. Wiesbaden, 1963, 1: 19-25.
[25]	徐中华, 王建华, 王卫东. 主体地下结构与支护结构相结合的复杂深基坑分析[J]. 岩土工程学报, 2006, 28(Sup1): 1355-1359. XU Z H, WANG J H, WANG W D. Analysis of a complicated deep excavations supported by substructures[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(Sup1): 1355-1359. (in Chinese)