浅谈几种常用本构模型在强夯沉降计算中的应用.docx

浅谈几种常用本构模型在强夯沉降计算中的应用

摘要：本论文以厦门航空港区港联动区与货运配套区围填工程，以Ⅱ-1b块强夯地基处理为例，结合分层总和法的地基沉降量计算理论对地基沉降量计算进行应用，并采用线弹性模型、非线性弹性模型（邓肯-张模型）、弹塑性模型（莫尔库仑、修正剑桥模型）等三种类型的四个模型进行强夯处理后场地工后沉降二维有限元分析，对各不同模型进行了分析对比，可为类似地基处理工程研究、设计提供了参考和借鉴。

关键词：强夯法，沉降，数值模拟

1本构模型研究

土作为天然地质材料在组成及构造上呈现出高度的各向异性、非均质性、非连续性和随机性，在力学性能上表现出强烈的非线性、非弹性和粘滞性，土的本构模型就是反映这些力学性态的数学表达式。一般认为，一个合理的土的本构模型应该具备理论上的严格性、参数上的易确定性和计算机实现的可能性。自Roscoe等(1958年～1963年)[70,71]创建剑桥模型至今，各国学者已发展数百个土的本构模型，其中包括不考虑时间因素的线弹性模型、非线弹性模型、弹塑性模型和近来发展起来的内时模型、损伤模型及结构性模型等，常用的模型只有极少数几个。

本文将结合土的本构模型发展历史，介绍了土的本构模型的一般建立方法、土的本构模型的发展进程和常用本构模型，最后结合厦门航空港区港联动区与货运配套区围填工程（二标段）工程对几种常用土的本构模型进行分析对比研究。

2几种常用本构模型在强夯沉降计算中的应用

2.1工程概况

现以厦门航空港区港联动区与货运配套区围填工程（二标段）为例说明一下这种方法的工程应用过程，

该工程为笔者负责设计的一填海造地项目，2007年10月本工程正式开工，施工单位进场，首先进行临时围堰的建设，随后进行吹填造地及地基处理，最后进行排水等配套设施的建设，截止2009年7月全部完工，工程实体满足设计要求。

厦门航空港区港联动区与货运配套区围填工程，位于厦门岛的北部，厦门大桥与集美大桥之间高崎国际机场西北侧的滩涂地带。本工程围填面积76.10×104m2；分为三个标段，分别为一标段（31.20×104m2）、二标段（32.40×104m2）和三标段（12.50×104m2）。工程主要内容为围填造地和软基处理。

地基处理方案为根据原泥面高程情况，Ⅱ－1a块与Ⅱ－1b块区回填中粗砂分别至高程2.0m与1.0m，再回填砂性土至设计标高，然后进行强夯处理。工程围填造地经地基处理后，地基承载力达到80kPa以上，工后残余沉降量小于25cm。强夯设计参数如下：

（1）强夯拟定点夯2遍，满夯1遍，点夯单击夯能要求2500kN.m；夯锤直径应不小于2.5m，夯点布置为正方形，第一遍强夯夯点与第二遍强夯夯点间距均为6m×6m，点夯每点击数不小于8击，每遍点夯的最后两击平均夯沉量小于10cm。第一遍强夯夯点与第二遍强夯夯点呈梅花型布置。

（2）普夯能量为1000KJ，夯锤直径应不小于2.5m，要求夯印搭接为1/3锤底直径，每点夯点总贯入量大于100cm时，需在原夯点上填料进行补夯，普夯最后两击平均夯沉量不大于5cm。

（3）锤重要求在200kN左右。

2.2几种常用土的本构模型

（1）摩尔一库仑(Mobr-Coulomb)模型

库仑(C．A．Coulomb)在土的摩擦试验、压剪试验或三轴试验的基础之上，于1773年提出了库仑破坏准则，即剪应力屈服准则，其准则方程为：

=C+atanφ??????????????????????????????(5.1)

其中，c为土的粘聚强度；φ为内摩擦角。如果已知三轴试件内破坏面与小主应力方向之间的倾角为βf，则由普通三轴试验的莫尔圆，将破坏面上的剪应力与法向应力代入库仑破坏准则，得到发生在某破坏面时主应力表达的破坏准则，即莫尔一库仑准则：

σ1-σ3=2ccosφ+(σ1+σ3)sinφ?????????(5.2)

试验研究表明，Mobr-Coulomb准则是符合岩土材料的屈服和破坏特征的。但是由于其屈服面在空间中的角度性质的影响，只要应力落在“脊梁”(棱角)附近，屈服函数沿曲面的外法线方向导数不易确定，则粘塑性的应变率不易确定。另外，在角锥顶点也存在不连续的问题，这就增加了使用上的困难，特别是限制了其在土工数值计算上的应用。

（2）Drucker-Prager(D--P)模型

经典的Tresca准则和Mises准则都没有考虑平均正有效应力对材料屈服性状的影响，为了考虑该影响条件，Drucker与Prager[5]于1952年提出了在应力空间中为一圆锥形屈服面的广义Mises屈服与破坏准则，或称为D-P屈服破坏准则。屈服函数为：

????(5.3)

其中，a，k为D-P准则材料常数，按照平面应变条件下的应力和塑性变形条件，Dmcker与Prager推导得出了，a、k与C