挡土墙后EPS板减压性能模型试验研究

郑俊杰 邵安迪 谢明星 景丹

摘要：在剛性挡土墙后设置柔性垫层能有效减小墙背土压力。聚苯乙烯土工泡沫（EPS）是一种常见的柔性材料，具有受力易压缩的特性。通过开展室内挡土墙模型箱试验，对铺设不同弹性模量以及不同厚度的EPS板的各工况下墙背土压力分布规律进行了研究，详细分析了弹性模量和厚度对EPS板在不同填土表面荷载作用下的减压性能的影响;并通过FLAC3D建立有限差分数值模型，对不同压缩量情况下土压力分布规律进行了分析。研究结果表明：EPS板能有效地减小挡土墙墙后土压力;EPS板的弹性模量越小、厚度越大、墙后填土表面的外荷载越大，EPS板的减压效果越明显;EPS板厚增加到一定程度后，EPS板的减压效果不会再有明显提升，EPS板厚达到0.1倍墙高即可获得最佳减压效果。

关键词：挡土墙;聚苯乙烯土工泡沫;模型试验;数值模拟;土压力

中图分类号：TU432 文献标志码：A 文章编号：20966717（2020）04000109

收稿日期：20200305

基金项目：国家重点研发计划（2016YFC0800200）;国家自然科学基金（51678267、51878313、51808243）

作者简介：郑俊杰（1967 ），男，教授，博士生导师，主要从事岩土工程与隧道研究，Email：zhengjj@hust.edu.cn。

Received：20200305

Foundation items：National Key Research and Development Program （No. 2016YFC0800200）; National Natural Science Foundation of China （No. 51678267， 51878313， 51808243）

Author brief：Zheng Junjie （1967 ）， professor， doctorial supervisor， main research interests： geotechnical engineering and tunnel engineering， Email： zhengjj@hust.edu.cn.

Model test study on decompression performance of EPS boards behind retaining walls

Zheng Junjie， Shao Andi， Xie Mingxing， Jing Dan

（School of Civil Engineering and Mechanics， Huazhong University of Science and Technology， Wuhan 430074， P.R. China）

Abstract： Setting a flexible pad cusion behind the rigid retaining wall can effectively reduce the lateral earth pressure. Expanded polystyrene （EPS） geofoam is a common flexible material that has the characteristic of easy compression. In this paper， through the retaining wall model test in laboratory， the lateral earth pressure distribution of different EPS boards with different elastic modulus and thickness is studied; the influence of elastic modulus and thickness on the decompression performance of EPS boards under different loadings on the filling surface is analyzed in detail. The influence of compression value of EPS board on the lateral earth pressure on retaining wall is also analyzed by FLAC3D numerical modeling. The research results show that the EPS board can effectively reduce the lateral earth pressure behind the retaining wall; the smaller the elastic modulus of the EPS board， and the larger the thickness， as well as the external load on the surface of the backfill wall， the decompression effect of the EPS board becomes more significant. After the EPS board thickness increases to a certain extent， the decompression effect of the EPS board will not increase significantly. The best decompression effect can be obtained when the EPS board thickness reaches 0.1 times the wall height.

Keywords：retaining wall; EPS; model test; numerical simulation; earth pressure

挡土墙是土木工程、交通工程领域重要的结构物之一，其主要功能是抵抗侧向土压力，并防止墙后土体发生滑动或坍塌，目前已广泛应用于房屋建筑、水利工程、公路与铁路工程中。减小挡土墙墙背土压力对于保障工程安全、降低工程造价具有重要意义。迄今为止，诸多学者对挡土墙土压力进行了大量试验研究。Terzaghi[12]最早开展模型试验，证实了挡土墙土压力沿墙高呈非线性分布，该结论在岳祖润等[3]的离心试验中得到再次验证。对挡土墙土压力的研究很早便不局限于静止土压力，陈页开[4]和周健等[5]通过试验对不同挡土墙变位模式下的土压力开展了进一步研究。此外，也有学者采用有限差分软件FLAC3D对挡土墙土压力进行研究。谢明星等[6]分别对有限填土条件下挡土墙的土压力大小和分布进行了研究。加筋土挡墙作为极具应用价值的挡土墙，其土压力分布规律的研究更为必要，贾良[7]和陈建峰等[8]开展了数值模拟。Green等[9]对某混凝土悬臂式挡土墙在动态荷载作用下的侧向土压力进行了分析。

聚苯乙烯土工泡沫（EPS）是一种具有优良性能的轻型高分子材料，具有密度低、压缩变形大、缓冲隔振和隔热效果好的优点，可被用来解决许多岩土工程中的难题，如用于结构隔振、边坡稳定、软基处理、地埋结构减荷、冻土路基等诸多方面[10]。EPS的压缩变形特性可以有效减小刚性挡土墙的墙背土压力，学者们对此进行了大量试验研究。Tsukamoto等[11]的试验验证了在挡土墙和墙背填土之间沿竖向铺设压缩性EPS板能显著减小刚性挡土墙的静止土压力，Ikizler等[12]对其机理进行了深入探讨。EPS除了能减小静态荷载作用下的土压力，Gaskin[13]和Bathurst等[14]的振动台试验结果表明，在振动荷载作用下，土工泡沫材料对侧向土压力的减压作用仍然显著，其最大减小量可以达到60%。Bathurst等[15]也对EPS缓冲层减轻刚性基底和挡土墙地震负荷的效果进行了进一步论证。汪益敏等[16]基于考虑时间因素的有限差分原理，建立了水平条分法的数值模型，对挡土墙后的EPS土工泡沫缓冲层进行了分析，主要分析其对减小振动作用下挡土墙受力与变形的作用。

笔者基于直剪试验和无侧限抗压强度试验，测试砂土和EPS土工泡沫的力学性能及相关参数。在此基础上，开展挡土墙模型试验，采用有限差分软件FLAC3D进行数值模拟，研究不同上覆荷载作用下EPS板的弹性模量和厚度等参数对墙背土压力的影响程度，为实际工程中挡土墙后设置EPS板提供设计参考。值得注意的是，同时考虑挡土墙的移动与EPS板的压缩两个因素较为复杂，笔者的研究仅为挡土墙在静止状态下EPS板的减压性能。

1试验砂土和EPS材料特性测定

1.1试验砂土抗剪强度参数测定

对试验砂土进行直剪试验，以得到其黏聚力c和内摩擦角φ，由于砂土不需考虑排水固结的问题，试验采用快剪方法。试验采用ZJ型应变控制式直剪仪，装置构造如图1所示。

通过直剪试验可得到剪应力剪切位移关系曲线，如图2所示。剪应力剪切位移曲线上的剪应力峰值即为该砂土试样的抗剪强度。线性拟合得出抗剪强度垂直压力关系曲线，如图3所示。

土体的抗剪强度关系为

由式（1）可知：直线的截距为土样的黏聚力c值，直线的斜率即为土样内摩擦角φ的正切值。由图3可知，此组试验砂土的黏聚力为0 kPa，内摩擦角为37.29°。按同样的方式进行第2组直剪试验，得到砂土的黏聚力为0 kPa，内摩擦角为38.66°。将两组试验数据取平均值，即砂土的黏聚力为0 kPa，内摩擦角为37.98°。

1.2EPS材料弹性模量值测定

模型箱试验共采用5种规格的EPS材料，分别对应不同的EPS密度。可通过对各规格EPS材料特制成的相同尺寸的圆柱体试样进行质量测量并求得密度;对圆柱体试样开展无侧限抗压强度试验，根据绘制出的弹性应变阶段的轴向应力应变曲线得出EPS材料弹性模量值。试验采用南京宁曦土壤仪器有限公司生产的TSZ系列应变控制式三轴仪，仪器如图4所示。

图5为通过试验绘制出的5种规格EPS试样的轴向应力应变曲线。由图5可见，EPS材料的应力应变关系大致在0～2%的轴向应变范围内为线性，即在初始阶段表现为线弹性材料，这与Kim等[17]选取的EPS材料线弹性区段一致。取该区段的应力应变曲线进行线性拟合，由式（2）可求出这5种规格EPS材料的弹性模量E。

对每种规格的EPS另取一试样进行第2组无侧限抗压强度试验，取两组试验结果的平均值作为弹性模量E，见表1。由表1可以看出，随着EPS材料密度的增大，其弹性模量也相应地增大。

2挡土墙模型箱试验

2.1挡土墙模型箱

试验使用的自制模型箱长0.8 m、宽0.4 m、高0.48 m。箱体外部框架由30 mm方钢焊接而成，顶部无盖，模型箱侧板与底板的材质均为20 mm厚ABS塑钢，侧板与底板嵌固在外部不锈钢框架上，并使用密封胶对缝隙进行密封处理。挡墙面板为50 mm不锈钢板，底端固定、上端可通过铆钉改变其倾角。通过改变可移动挡板的位置来改变挡土墙后的填土宽度。EPS板均为定制，其面积与挡土墙面板尺寸一致。挡土墙模型箱如图6所示。

2.2试验步骤

1）模型箱就位，固定挡墙面板，将EPS板用玻璃胶固定在挡墙面板内壁。

2）在挡墙面板特定高度处（0、0.1、0.2、0.3、0.4 m），用玻璃膠将土压力传感器固定于EPS板表面，具体布置如图7所示。梳理接线并从挡墙面板侧边缘引出，将土压力传感器与应变采集仪连接，试验采用江苏东华测试技术有限公司的DH3816N静态应变测试系统，如图8所示。

3）对土压力传感器的白色通气管吹气，传感器的输出值由0开始有明显跃升，停止吹气输出值又逐渐减小为0，则表明土压力传感器安装正确。做好所有相关的准备工作后开始采样（采用5 s一次的定时采样）。

4）向模型箱内回填砂土，采用人工落雨法，使砂土从高于箱顶20 cm处自由下落，每回填10 cm高砂土静置30 min，分层铺设直至挡土墙顶部。待测试系统实时窗口中输出的土压力大小趋于稳定时测得的值即为无外加荷载时的墙背土压力。

5）稳定后开始逐级加载，如图9所示。在木板上均匀放置砝码来模拟施加在墙后填土表面的均布荷载，荷载均施加在墙后宽43 cm的填土表面。每增加一级荷载，需等待测量土压力值稳定后再增加下一级荷载。各级荷载大小见表2。

6）在砂土回填和砝码加载阶段始终保持静态应变采集仪的工作，测量和记录在各级荷载下，土压力沿挡墙高度的变化，并对数据结果进行分析。

7）卸掉荷载，清空砂土，取下土压力传感器，更换EPS板按同样的步骤进行下一工况的试验。完成每一工况试验后填写试验说明。试验包含8个工况，见表3。

2.3试验结果分析

2.3.1墙背土压力理论值计算

对于不含EPS板的工况1，理论上可计算静止土压力沿挡墙墙高的分布。对于本试验，试验砂样内摩擦角为φ=37.98°，重度γ=16 kN/m3，墙高H=0.48 m，静止土压力系数可根据式（3）求得。

由此计算出无荷载情况下墙背土压力的分布。

图10为无荷载情况下土压力实测值与理论值对比图。由图10可看出：由试验得到的实测土压力强度与理论值较为接近，在靠近墙踵处，实测值略小于理论值，但相差不大，证明了模型箱试验结果的可靠性。

2.3.2EPS板厚度的影響

图11为在无荷载、三级荷载、五级荷载和七级荷载条件下，工况1、2、3和4中土压力沿墙高的分布图。4个工况中EPS的弹性模量E相同，厚度t不同。对4种荷载条件下的土压力分布进行对比分析，结果表明：墙后填土表面荷载的增加对挡土墙靠近顶部部分影响较大，上部墙背土压力强度增加较为明显;放置EPS板后，墙背土压力相对于无EPS工况有明显减小，且在工况1下存在较大土压力强度对应的高度处减压效果更明显;从无荷载到七级荷载，墙背土压力越大，EPS的减压效果越明显;在同级荷载下，EPS板厚度越大时墙背土压力减小得越明显，这是由于随着EPS板厚度的增加，其压缩量也会相应增加，也就允许墙后土体发生更大的变形。

2.3.3EPS板弹性模量的影响

图 12为在无外荷载、三级荷载、五级荷载、七级荷载条件下，工况1、2、5、6、7和8中土压力沿墙高的分布图。除工况1未铺设EPS外，其余5个工况中EPS的厚度t相同，弹性模量E不同。

对4种荷载条件下的土压力分布进行对比分析，结果表明：墙后填土表面荷载的增加对挡土墙靠近顶部部分影响较大，墙背上部土压力强度增加较明显;放置EPS板后，墙背土压力相对于无EPS时有明显减小，且在工况1下存在较大土压力强度对应的高度处减压效果越明显的情况;从无外荷载到七级荷载，墙背土压力越大，EPS的减压效果就越明显;EPS材料弹性模量E越小，EPS板越柔，越容易发生压缩和变形，这使得采用E值更小的EPS板时墙背土压力减小的现象更为明显。

观察未施加荷载时无EPS（工况1）和有EPS工况下的土压力分布。明显看到工况1的土压力分布从上至下接近线性增大，而其余工况下的土压力强度均呈现从上至下先增大后减小的分布趋势。这是由于在有EPS的工况下，EPS产生的压缩量可视为挡土墙产生的位移，这与杨明辉等[18]通过试验得到的平动模式下的土压力强度分布规律相吻合。有无EPS工况下不同的土压力分布规律也表明EPS有均匀土压力分布的作用。

3数值模拟

3.1数值模型的建立

在试验过程中，由于EPS板的压缩量过小，用位移计无法测出，因此，采用数值模拟方法对EPS的压缩量与墙背土压力的分布之间的关系进行分析。

采用有限差分软件FLAC3D进行建模，所建数值模型尺寸均以模型箱的实际尺寸为标准。模型中，填土采用摩尔库仑模型，EPS材料、挡墙面板和模型箱箱壁均采用弹性模型。材料具体参数的选取见表4。根据挡土墙模型箱的实际情况，模型的边界条件为：模型底部用全约束限制所有方向的位移;可移动挡板以及挡墙面板对应的面则限制其水平方向的位移。建模过程中荷载的施加完全按照试验工况中的从无荷载逐级施加到七级荷载。

3.2数值模拟结果分析

3.2.1土压力强度模拟值与实测值对比

图13为工况1、2、7土压力强度模拟值与实测值的对比。由图13（a）可以看出，工况1无荷载情况下的土压力强度模拟值与实测值吻合良好，区别仅在墙踵处土压力强度模拟值减小至零。由图13（a）～（c）均可看到，在不同工况及不同荷载作用下，模拟值与实测值均能较好地吻合。将工况2、工况7的土压力强度分布情况与工况1进行对比，可以明显看到工况1中土压力强度从墙顶至墙踵大致呈线性增加，只是在接近墙踵处减小至零;而工况2及工况7均呈现挡墙中部土压力强度大、上下部较小的分布规律。如前文所述，由于EPS板压缩可视作挡土墙发生了移动从而导致了该现象。以上结果验证了数值模型的正确性。

3.2.2EPS压缩量分析

图14为工况2中EPS压缩量沿墙高变化分布图。由图14可以看出：EPS板沿墙高压缩量与土压力强度分布情况一致，尤其在无荷载作用时明显呈现出EPS板中部压缩量大、上下部压缩量较小的情况;随着荷载级数增加，EPS板的压缩量也增大，但由于荷载增加对墙体上部土压力影响更大，故EPS板上部的压缩量增加得更明显。这也导致施加外荷载后，EPS板的中上部压缩量均大于下部压缩量。

图15为七级荷载下A规格3种厚度的EPS板沿墙高的压缩量示意图。由图15可以看出：七级荷载作用下，EPS板中上部的压缩量大、下部压缩量小;厚度越大的EPS板，其压缩量也越大，允许土体产生更大的变形，因此，减小的墙背土压力也就越多。这也验证了板厚t越大EPS的减压效果越显著的结论。

图16为七级荷载下EPS厚度为4 cm时，不同弹性模量的EPS压缩量沿墙高变化曲线。由图16可以看出：在EPS板厚度t一定时，弹性模量E越小，其压缩量也就越大，减小的墙背土压力也就越多。弹性模量E差距较大的EPS板压缩量的差别在图中表现得尤为显著，如EPS（A）和EPS（E）。这也验证了弹性模量E越小EPS板减压效果越显著的结论。

3.3EPS静止土压力减压系数

挡土墙位移、墙高、墙后土体参数、EPS板厚度和密度等参数均对设置EPS板的挡土墙土压力有直接影响。该研究挡土墙位移为0，为进一步减少影响土压力的参数，定义EPS静止土压力减压系数

式中：En为铺设有EPS板的不同工况下静止土压力，E0为无EPS板工况下静止土压力。利用式（4）得出K值即可近似忽略墙高对减压效果的影响。针对试验中采用的砂土，基于数值模拟计算结果，对K值进行拟合。

EPS板密度ρ=6 kg/m3时，K随EPS板厚度t变化关系的拟合表达式为

EPS板厚度t=4 cm时，K随EPS板密度ρ的变化关系拟合表达式为

利用多项式拟合EPS静止土压力减压系数K关于密度ρ和厚度t的变化关系为

根据式（5）可得出K随EPS板厚度t的变化曲线，如图17所示。由图17可以看出：EPS板密度一定时，K均随板厚t增大而减小。如密度为6 kg/m3的EPS板板厚为1 cm时K值为0.712，板厚为6 cm时K值减小为0.547;其减小趋势呈非线性，随板厚t增加，在EPS板厚达到5 cm后，K值减小的趋势明显减缓。

这表明增大EPS板厚虽然能提高减压效果，但并非每单位厚度的EPS板都发挥了相同程度的减压作用，即部分EPS材料未能充分发挥作用。板厚增加到一定程度后，K值不会再有明显减小，但此时EPS板的压缩量仍在不断增大，故在实际工程中应根据经济效益和土体变形允许范围选择合适厚度的EPS板。由于EPS板压缩量可视为挡土墙位移，挡土墙主动土压力随挡土墙位移的增大而减小。根据上述分析，墙高H为0.48 m时，EPS板厚达到5 cm即为较佳状态，即EPS板厚达到0.1H即可获得最佳减压效果，此时，EPS板的平均压缩量约为0.000 1H。

不同EPS板厚度工况下，K随EPS板密度ρ的变化曲线如图18所示。由图18可以看出：EPS板厚度一定时，K均随密度ρ减小而减小，板厚为4 cm的EPS板密度为15 kg/m3时，K值为0.708，密度为6 kg/m3时，K值减小为0.565，且K值随密度ρ减小而减小的趋势大致呈线性，这与式（6）反映的变化规律一致。进一步证实了较低密度的EPS板能发挥较大的减压作用。K值越小，表明EPS的减压效果越好。K值最小值约为0.5，此时静止土压力仅为无EPS时的1/2，可见将EPS应用于挡土墙减压的研究极具价值。

4结论

为探究挡土墙后EPS板的减压性能，通过室内试验和FLAC3D数值模拟对铺设EPS板的挡土墙墙背土压力进行了研究，得到以下主要结论：

1）挡土墙模型箱试验和数值模拟结果均表明：有EPS工况下呈现出墙体中部土压力强度大，墙顶、墙踵处小的分布规律，这是由于EPS板产生的压缩量可视作挡墙移动;随着EPS弹性模量E减小，EPS材料的减压效果会增加;随着EPS板厚t增加，EPS材料的减压性能会提高;施加的外荷载越大，EPS的减压效果越好。

2）依托模型箱试验进行的数值模拟结果表明：EPS板沿墙高的压缩量与土压力强度分布情况一致，且EPS弹性模量E越小，板厚t越大，压缩量就越大。

3）根据模拟结果，拟合出试验土体的EPS静止土压力减压系数K关于EPS密度和厚度的計算公式，结果表明：K值随着密度ρ减小而减小，且变化趋势接近于线性;K值随着板厚t增大而减小，但随板厚t增加，K值减小的趋势明显减缓。EPS板厚达到0.1H可获得最佳减压效果。参考文献：

[1] TERZAGHI K. Large retaining wall tests IIpressure of saturated sand [J]. Engineering News Record， 1934， 112（22）： 259262， 316318， 403406， 503508.

[2] TERZAGHI K. Record earth pressure testing machine [J]. Engineering News Record， 1932， 109（29）： 365369.

[3] 岳祖润， 彭胤宗， 张师德. 压实粘性填土挡土墙土压力离心模型试验[J]. 岩土工程学报， 1992， 14（6）： 9096.

YUE Z R， PENG Y Z， ZHANG S D. Centrifuge model tests on lateral pressure on walls retaining compacted clayey backfill [J]. Chinese Journal of Geotechnical Engineering， 1992， 14（6）： 9096. （in Chinese）

[4] 陈页开. 挡土墙上土压力的试验研究与数值分析[D]. 杭州： 浙江大学， 2001.

CHEN Y K. Model test and numerical analysis of earth pressures on retaining wall [D]. Hangzhou： Zhejiang University， 2001.（in Chinese）

[5] 周健， 高冰， 彭述权. 不同位移模式下挡土墙的模型试验及数值模拟[J]. 岩石力学与工程学报， 2011， 30（Sup2）： 37213727.

ZHOU J， GAO B， PENG S Q. Model test and numerical simulation of retaining wall under different displacement modes [J]. Chinese Journal of Rock Mechanics and Engineering， 2011， 30（Sup2）： 37213727. （in Chinese）

[6] 谢明星， 郑俊杰， 曹文昭， 等. 有限填土路堤挡土墙主动土压力研究[J]. 华中科技大学学报（自然科学版）， 2019， 47（2）： 16.

XIE M X， ZHENG J J， CAO W Z， et al. Study of active earth pressure against embankment retaining wall of limited backfill [J]. Journal of Huazhong University of Science and Technology （Natural Science Edition）， 2019， 47（2）： 16. （in Chinese）

[7] 贾良. 加筋土挡墙的土压力试验与数值分析[J]. 深圳大学学报（理工版）， 2016， 33（1）： 5561.

JIA L. Soil pressure test and numerical calculation about reinforced earth retaining wall [J]. Journal of Shenzhen University Science and Engineering， 2016， 33（1）： 5561. （in Chinese）

[8] 陈建峰， 郑希， 李辉利， 等. 软土地基加筋土挡墙二维离散连续耦合模拟[J]. 土木工程学报， 2011， 44（Sup2）： 3336.

CHEN J F， ZHENG X， LI H L， et al. Coupling discontinuumcontinuum simulation of a geogridreinforced soil retaining wall（RSW） on soft soil [J]. China Civil Engineering Journal， 2011， 44（Sup2）： 3336. （in Chinese）

[9] GREEN R A， EBELING R M. Modeling the dynamic response of cantilever earthretaining walls using FLAC[C]//3rd International Symposium on FLAC： Numerical Modeling in Geomechanics， Sudbury， Canada， 2003.

[10] HORVATH J. The compressible inclusion function of EPS geofoam [J]. Geotextiles and Geomembranes， 1997， 15（1/2/3）： 77120.

[11] TSUKAMOTO Y， ISHIHARA K， KON H， et al. Use of compressible expanded polystyrene blocks and geogrids for retaining wall structures [J]. Soils and Foundations， 2002， 42（4）： 2941.

[12] IKIZLER S B， AYTEKIN M， NAS E. Laboratory study of expanded polystyrene （EPS） geofoam used with expansive soils [J]. Geotextiles and Geomembranes， 2008， 26（2）： 189195.

[13] GASKIN P A. Geofoam buffers for rigid walls： an investigation into the use of expanded polystyrene for seismic buffers [D]. Kingston， Ontario： Queens University， 2000.

[14] BATHURST R J， ZARNANI S， GASKIN A. Shaking table testing of geofoam seismic buffers [J]. Soil Dynamics and Earthquake Engineering， 2007， 27（4）： 324332.

[15] BATHURST R J， ZARNANI S. Earthquake load attenuation using EPS geofoam buffers in rigid wall applications [J]. Indian Geotechnical Journal， 2013， 43（4）： 283291.

[16] 汪益敏， BATHURST R J. 振動条件下EPS缓冲层对挡土墙受力与变形影响的水平条分法模型研究[J]. 土木工程学报， 2008， 41（10）： 7380.

WANG Y M， BATHURST R J. Horizontal slice method for force and displacement analysis of EPS geofoam seismic buffers for rigid retaining walls [J]. China Civil Engineering Journal， 2008， 41（10）： 7380. （in Chinese）

[17] KIM H， WITTHOEFT A F， KIM D. Numerical study of earth pressure reduction on rigid walls using EPS geofoam inclusions [J]. Geosynthetics International， 2018， 25（2）： 180199.

[18] 杨明辉， 戴夏斌， 赵明华， 等. 墙后有限宽度无黏性土主动土压力试验研究[J]. 岩土工程学报， 2016， 38（1）： 131137.

YANG M H， DAI X B， ZHAO M H， et al. Experimental study on active earth pressure of cohesionless soil with limited width behind retaining wall [J]. Chinese Journal of Geotechnical Engineering， 2016， 38（1）： 131137. （in Chinese）

（编辑章润红）