砂土地基上圆形浅基础三维破坏包络面的理论研究
开篇：润墨网以专业的文秘视角，为您筛选了一篇砂土地基上圆形浅基础三维破坏包络面的理论研究范文，如需获取更多写作素材，在线客服老师一对一协助。欢迎您的阅读与分享！
摘要：对砂土地基上圆形浅基础在竖向荷载V、水平荷载H及力矩M复合加载条件下的承载力进行了系统的三维有限元分析。在分析中，砂土假定为纯摩擦材料，遵循基于MohrCoulomb破坏准则的理想弹塑性本构关系。首先，对圆形浅基础的竖向承载力进行了有限元计算， 并与滑移线解法进行了对比，两种方法所得结果比较吻合。进而探讨了砂土内摩擦角对于基础在VH、VM荷载平面与VHM三维荷载空间内的破坏包络轨迹的影响。计算结果表明，与不排水情况下软黏土地基上基础破坏包络面相比，砂土地基上圆形浅基础的破坏包络面形状有较大差异， 但VH和VM平面内的破坏包络面形状仍具有较好的归一化特性。基于有限元计算结果，建立了圆形浅基础在VHM三维荷载空间内的破坏包络面方程， 该方程可用来合理评价复合加载条件下砂土地基上圆形浅基础的整体稳定性。
关键词：圆形浅基础；复合加载；砂土；破坏包络面；承载力
中图分类号：
TU311； P315.9文献标志码：A文章编号：16744764（2015）03006307
Abstract：
The bearing capacity of circular shallow footings subjected to combined loading composed of vertical component V， horizontal component H and moment M on sand is explored through three dimensional finite element numerical analyses. The sand is assumed to be purely frictional and obey elastoperfectly plastic relationship following the MohrCoulomb failure criterion. The capacity of circular footing under centrally vertical load is computed， and the calculated results agree well with those derived from slipline field theory. Then the influence of soil friction angle on the failure locus in the VH and VM load plane and VHM load space is investigated. The results show that the shape of failure envelopes of shallow foundations on sand is different from that on undrained clay， whilst normalization of the failure envelope on the VH and VM load plane by the pure vertical bearing capacity allow it to be generalized for varying conditions. Based on finite element analyses， the equation of failure envelope in the VHM three dimensional load space is proposed， which can be used to evaluate the general stability of circular footings on sand under combined loading condition.
Key words：circular footing； combined loading； sand； failure envelope； bearing capacity
对于圆形浅基础，传统的地基承载力分析理论主要有Terzaghi、Vesic、Hansen、Meyerhoff等提出的各种经验、半经验公式，对于倾斜与偏心荷载作用下的浅基础稳定性问题，一般是基于Terzaghi竖向承载力公式，分别通过引入荷载倾斜系数与Meyerhoff有效宽度假定来考虑水平荷载和偏心荷载对竖向承载力的影响，对此栾茂田等[1]已经进行了比较全面的评述。
随着石油、天然气和金属矿物等海洋资源的大力开发，各种海洋基础稳定性评价方面的问题得到了高度重视。与陆地上基础相比，海洋基础除了承受竖向荷载V以外，通常还要抵抗波浪和风暴等所引起的水平荷载H与力矩M。在这种复合加载情况下，海洋基础一般不会在单纯的竖向荷载作用下达到极限平衡状态，而是在竖向荷载、水平荷载与力矩的不同组合条件下发生失稳破坏。因此，将传统的地基承载力理论用于海洋浅基础稳定性评价时，可能出现较大偏差，从而不适合含有较大水平荷载和力矩的情况[2]。为了解决这个问题，部分学者[12]通过系统研究提出了破坏包络面理论，认为在复合加载条件下，地基达到整体破坏时各个荷载分量的组合在三维荷载空间（V，H，M）中将形成一个不依赖于加载路径的外凸曲面，其方程可由引起地基失稳时的各个荷载分量显式表达为f（V，H，M）=0。根据实际的受力状态与该破坏包络面之间的相对位置关系，可以直观评价设计荷载状态下海洋基础的整体稳定性。
目前所开展的研究工作大多针对软黏土地基，砂土地基上圆形浅基础破坏包络面方程的研究较少，并仅限于模型试验研究。Martin等[3]对于复合加载条件下黏土地基上纺锤形基础的力学响应进行了比较系统的室内小比尺模型试验。在此基础上，Gottardi等[4]、Cocjin等[5]分别对于密砂上圆形和条形浅基础开展了一系列模型试验，并对试验数据进行拟合，建议了砂土地基上浅基础的破坏包络面方程。Cassidy等[6]针对松散钙质砂地基开展了小比尺复合加载试验研究，主要探讨了宏单元模型中的硬化准则和流动法则。但这些试验工作都是针对某种相对密度或内摩擦角的砂土，因此，本文对于共面复合加载条件下砂土地基上圆形浅基础的承载力进行比较系统的有限元分析，探讨了砂土内摩擦角（对于圆形浅基础在VH、VM荷载平面与VHM荷载空间内的破坏包络轨迹的影响，并与已有的室内模型试验结果进行了对比。
1有限元模型
直径D=1 m的圆形浅基础位于砂土地基上。不考虑浅基础本身的变形，因此采用离散刚体单元模拟。地基模型的半径和深度都取为5D，经过试算，可以消除有限元模型中地基的边界效应影响。根据赵少飞等[7]的建议，土体单元类型选择8节点缩减积分实体单元（C3D8R），在浅基础附近的局部区域加密网格，单元数为26 240，如图1所示。
在分析中，砂土假定为纯摩擦材料，重度取为γ=20 kN/m3，采用基于MohrCoulomb破坏准则的理想弹塑性本构模型，变形模量E=100 MPa，泊松比v=0.3。很多实验结果已经表明，砂土剪胀角φ低于内摩擦角，但为了与常用的竖向承载力解进行对比，仍然采用了相关联流动法则，即取φ = 。土体为纯摩擦材料，基础与地基之间不能传递拉应力，因此，基础和地基间设为完全黏结约束，而未设置接触单元。
为了避免纯摩擦材料计算难以收敛的数值问题，同时又不影响材料真实的力学响应，需要在模型中人为赋予砂土一个很小的黏聚力值。根据Zhu等[8]利用ABAQUS软件计算矩形基础承载力时的建议，在本文考虑的砂土内摩擦角范围10°～30°内，当
2竖向承载力
对于内摩擦角=15°、20°、25°与30° 4种情况，分别进行了竖向承载力计算。在计算过程中，采用位移控制加载，当得到的基底中心处的节点反力V竖向位移v曲线斜率陡降或接近零时，与该状态相对应的基底反力即为地基的竖向承载力，如图2所示。
同时，也采用Martin[9]提出的滑移线方法和Terzaghi圆形浅基础公式对此问题进行了求解，几种方法所得竖向承载力（单位：kN）都列于表1，通过比较可以看到，在≥20°时有限元与滑移线解法所得结果都要比Terzaghi公式计算结果高，而在=15°时，有限元计算结果比其它两种方法偏低。有限元法与滑移线解法所得结果之间相差不大，最大误差为8%，从而说明本文有限元模型是基本合理的。
3复合加载数值试验方法
3.1swipe型加载方法
构建不同荷载平面上完整的破坏包络面，对于软黏土地基，只需要一条swipe加载路径[1011]，但对于砂土地基，则需要两条加载路径，如图3所示。加载路径I与软黏土地基相同，包括两个加载步骤：1）沿i方向从初始状态开始施加位移ui直到i方向反力达到极限值；2）固定i方向的位移，沿j方向施加位移uj直到j方向对应的反力Fj不随位移增大而改变，此时步骤2）所形成的加载轨迹可以近似作为ij平面上破坏包络面的一部分。加载路径II含有一个加载步骤，从初始状态开始约束i方向自由度，沿j方向施加位移uj直到j方向对应的反力Fj不随位移增大而改变，该步骤所形成的加载轨迹可作为包络面的另外一部分。
3.2probe型加载方法
probe型加载方法包括固定位移比加载、固定荷载比加载方法等，范庆来等[10]已经进行了比较详细的介绍，在此不再赘述。最近，赵少飞[7]建立了一种较为直观的probe型加载方法，该方法包含如下两个步骤：1）通过荷载控制，沿i方向（一般为竖向）在基础上直接施加一定的荷载分量Fi；2）保持所施加的荷载分量Fi不变，沿j方向进行位移控制加载，直到相应方向的反力Fj不再随位移增加而改变，由此可确定出破坏包络面上的一个点（Fi， Fj），如图3所示。通过多次加载，即可拟合一个完整的包络面。本文在构建VH、VM荷载平面上的包络面时，采用了该方法。
3.3荷载位移联合搜索方法
为了得到圆形浅基础在VHM荷载空间内的三维破坏包络面，需采用荷载位移联合搜索方法[1011]。这个方法包含如下3个步骤：1） 通过荷载控制，在基底中心施加一定大小的竖向荷载分量V；2） 保持该竖向荷载分量不变，进行HM荷载平面上的swipe型加载，得到破坏包络面的近似形式；3） 在此基础上，进行若干次probe型加载，确定最终的破坏包络面。
4VH平面上的破坏包络轨迹
联合采用swipe型与probe型两种数值加载方法，对于圆形浅基础在VH平面内的破坏包络轨迹进行研究。所得到的VH荷载平面内的破坏包络轨迹如图4（a）所示，采用竖向承载力Vult进行归一化后，可得VH荷载平面内的破坏包络面形状如图4（b）所示。
通过图4（a）可以看到，随着砂土内摩擦角的增大，VH荷载平面内的破坏包络面大小在不断增长。在内摩擦角=30°时，图4（a）还具体给出了swipe型加载路径I、II与probe型加载得到的3个数据点（分别是在V / Vult =0.3、0.5及0.7情况下得到的），可以看到swipe型加载路径I与路径II在V / Vult =0.5处汇合，构成了一个完整的包络面。probe型加载得到的数据点与swipe型加载路径基本吻合，考虑到数值计算误差，可以表明破坏包络面是不依赖于加载路径的。在其它内摩擦角情况下，也具有相同规律，在图4（a）中不再一一表达。
根据图4（b）可看出，破坏包络面形状类似于橄榄球形，基础所能承受的最大水平荷载Hmax大致出现在竖向荷载水平V/Vult=0.5处，而且对于不同内摩擦角情况下，Hmax≈0.13Vult。Gottardi等[4]针对内摩擦角=42.3°以及Bienen等[12]针对=34.3°的砂土所进行的模型试验也得到了基本一致的结论Hmax≈0.12Vult。Cassidy等[6]对于松散钙质砂也进行了试验，得到的结果表明Hmax≈0.15Vult。
因此，如图4（b）所示，不同内摩擦角情况下，采用竖向极限承载力Vult进行归一化后，包络面形状几乎完全重合，可采用式（1）进行描述。
HVult=4h0VVult1-VVult（1）
式中：h0=Hmax/Vult，其取值范围在0.12～0.15之间，对于本文有限元计算结果，h0=0.13。在图4（b）中还列出了Loukidis等[13]建立的条形浅基础破坏包络面方程，可见圆形与条形浅基础的破坏包络面形状存在一定差异。因此，在分析基础稳定性时，必须考虑其三维效应。
对于软黏土地基VH包络面，水平荷载最大值位于V=0，随着竖向荷载水平增大，基础承受水平荷载的能力不断下降[10]，这显然与砂土地基上基础破坏包络面特性有明显差异。
5VM平面上的破坏包络轨迹
采用类似数值加载方法，对于圆形浅基础在VM平面内的破坏包络轨迹进行了研究，其中力矩M是通过在基底中心处施加转角边界条件控制加载。所得到的VM荷载平面内的破坏包络轨迹如图5（a）所示。采用基础直径D与竖向承载力Vult之乘积DVult对于力矩M无量纲化，可得VM荷载平面内的归一化破坏包络面形状如图5（b）所示。
通过图5（a）可以看到，VM荷载平面内的破坏包络面大小也随着砂土内摩擦角的增大而在不断增长。归一化后的破坏包络面形状基本重合，也类似于橄榄球形，基础所能承受的最大力矩荷载Mmax大致出现在竖向荷载水平V/Vult=0.5处，Mmax=008DVult，对应着偏心距为e/D=1/625。Gottardi等[4]根据密砂的模型试验得到Mmax=01DVult，Bienen等[11]根据松砂上的试验结果得到的结论则是Mmax=0.075DVult，因此，可以认为，VM平面内砂土地基上圆形浅基础的归一化最大弯矩承载力m0=Mmax/（DVult）在0.075～0.1之间。不同内摩擦角情况下，包络面形状基本重合，如图5（b）所示，可采用如下抛物线方程式（2）进行描述。
MDVult=4m0VVult1-VVult
（2）
式中：m0取值范围在0.075～0.1之间，对于本文有限元计算结果，m0=0.08。与软土地基上基础VM包络面方程[10]进行比较，可以看到偏心荷载情况下砂土与黏土地基上破坏包络面特性也存在显著差异。
通过比较，可以看到式（4）的模拟结果较好，因此，采用如式（4）所示的三维破坏包络面方程来评价复合加载条件下砂土地基上圆形浅基础承载力是合理的。在实际应用中，只需根据竖向极限承载力公式或其它方法确定相应的竖向承载力Vult，进而根据方程（4），就可以得到破坏包络面的显式表达式f（V， H， M/D）=0，如图8所示。若浅基础设计荷载组合点（V，H，M）位于包络面上时，说明地基处于承载能力极限状态。当荷载组合点（V，H，M）处于该包络面内部，则浅基础是整体稳定的，反之，则将发生失稳破坏，在工程中应对基础进行重新设计[1416]。
7结论
1） 与不排水情况下软黏土地基上基础破坏包络面相比，砂土地基上圆形浅基础的破坏包络面形状有较大差异， 但VH和VM平面内的破坏包络面形状仍具有较好的归一化特性，可用抛物线方程进行表达。
2） 在一定竖向荷载水平下，不同内摩擦角情况下的HM破坏包络轨迹基本重合，形状为具有一定偏心度的椭圆。
3） 根据计算结果，提出了砂土地基上圆形浅基础三维破坏包络面方程。经过初步验证，该方程可以用来评价圆形浅基础在共面复合加载条件下是否稳定。
参考文献：
[1]
栾茂田，范庆来，王忠涛，等.海洋工程中若干土动力学问题的研究进展[C]//第七届全国土动力学学术会议论文集.北京：清华大学出版社，2006：2638.
Luan M T， Fan Q L， Wang Z T， et al. Astateoftheart review on some key issues of soil dynamics in ocean and offshore engineering [C]// The 8th National Conference on Soil Dynamics. Beijing： Tsinghua University Press， 2006： 2638. （in Chinese）
[2] Zhang Y， Bienen B， Cassidy M J， et al. Undrained bearing capacity of deeply buried flat circular footings under general loading [J]. Journal of Geotechnical and Geoenvironmental Engineering， ASCE， 2012， 138（3）： 385397.
[3] Martin C M， Houlsby G T. Combined loading of spudcan foundations on clay： laboratory tests [J]. Geotechnique， 2000， 50（4）： 325338.
[4] Gottardi G， Houlsby G T， Butterfield R. Plastic response of circular footings on sand under general planar loading [J]. Geotechnique， 1999， 49（4）： 453469.
[5] Cocjin M， Kusakabe O. Centrifuge observations on combined loading of a strip footing on dense sand [J]. Geotechnique， 2013， 63（5）： 427433.
[6] Cassidy M J， Byrne B W， Houlsby G T. Modelling the behaviour of circular footings under combined loading on loose carbonate sand [J]. Geotechnique， 2002， 52（10）： 705712.
[7] 赵少飞.复合加载条件下海洋地基承载力特性数值分析方法研究[D].辽宁 大连：大连理工大学，2006.
Zhao S F. A study on numerical methods for analyses of bearing capacity behavior of offshore foundations under combined loading [D]. Dalian： Dalian University of Technology， 2006. （in Chinese）
[8] Zhu M， Michalowski R L. Shape factors for limit loads on square and rectangular footings [J]. Journal of Geotechnical and Geoenvironmental Engineering， ASCE， 2005， 131（2）： 223231.
[9] Martin C M. Exact bearing capacity calculations using the method of characteristics [C]// Proceedings of the 11th International Conference on Computer Methods and Advances in Geomechanics. Turin： [s. n.]， 2005， 441450.
[10] Fan Q L， Luan M T， Gong X B. A unified equation of failure envelope for skirted foundations in normally consolidated clay [J]. Chinese Journal of Geotechnical Engineering， 2012， 34（10）： 19171924.
[11] 范庆来，栾茂田.各向异性软黏土地基上浅基础破坏包络面研究[J].岩石力学与工程学报，2010，29（11）：23622369.
Fan Q L， Luan M T. Study of failure envelope of shallow foundation on anisotropic soft clay [J]. Chinese Journal of Rock Mechanics and Engineering， 2010， 29（11）： 23622369. （in Chinese）
[12] Bienen B， Byrne B W， Houlsby G T， et al. Investigating sixdegreeoffreedom loading of shallow foundations on sand [J]. Geotechnique， 2006， 56（6）： 367379.
[13] Loukidis D， Chakraborty T， Salgado R. Bearing capacity of strip footings on purely frictional soil under eccentric and inclined loads [J]. Canadian Geotechnical Journal， 2008， 45（6）： 768787.
[14] 范庆来，郑静.复合加载模式下海上风机圆形浅基础亚塑性宏单元模型[J].土木建筑与环境工程，2014，36（3）：5963.
Fan Q L， Zheng J. A hypoplastic macroelement model for circular shallow foundations of offshore wind turbines under combined loading [J]. Journal of Civil， Architectural & Environmental Engineering， 2014， 36（3）：5963. （in Chinese）
[15] 范庆来，赵海涛，郑静，等.非共面复合加载条件下桶形基础稳定性研究[J].岩土力学，2013，34（12）：36433648.
Fan Q L， Zhao H T， Zheng J， et al. Stability of bucket foundations under noncoplanar combined loading [J].Rock and Soil Mechanics， 2013， 34（12）： 36433648. （in Chinese）
[16] 武科，马明月，栾茂田，等.扭剪荷载作用下桶形基础承载性能的弹塑性有限元数值分析[J].应用基础与工程科学学报，2012，20（5）：777786.
Wu K， Ma M Y， Luan M T， et al. Numerical analysis of bearing capacity behavior of bucket foundation subjected to torsion shear loading based on elastoplastic FEM [J].Journal of Basic Science and Engineering， 2012， 20（5）： 777786. （in Chinese）