基于上限有限元法的下伏溶洞地基极限承载力分析_赵明华.pdfVIP

水 利 学 报 2015年6月 SHUILI XUEBAO 第46卷 增刊1 文章编号：0559-9350（2015）S1-0006-06 基于上限有限元法的下伏溶洞地基极限承载力分析 赵明华，卢晓明，张 锐 （湖南大学 岩土工程研究所，湖南 长沙 410082） 摘要：简要地介绍上限有限元法的基本原理与计算过程，结合岩溶区下伏溶洞地基受力特点，提出合理的计算假 定，建立了一个可考虑多种主要因素影响的数值计算模型，利用上限有限元法，计算出了多种工况下的地基极限 承载力上限解，并分析溶洞顶板厚度、溶洞宽度和内摩擦角对地基极限承载力的影响规律。最后，在对大量能量 耗散图和速度场图分析的基础上，提出三类典型的下伏溶洞地基破坏模式。研究结果表明，地基极限承载力随着 内摩擦角和溶洞顶板厚度的增大而增大，随着溶洞宽度的增大而减小；当溶洞顶板厚度较小时，内摩擦角对地基 承载力影响较小，当内摩擦角较大时，溶洞宽度和溶洞顶板厚度对地基承载力的影响才显著；岩溶区下伏溶洞地 基存在冲切破坏、冲切剪压破坏和Prandtl破坏三类破坏模式。 关键词：地基基础；溶洞；极限承载力；上限法；有限元分析 中图分类号：TU457 文献标识码：A doi：10.13243/ki.slxb.2015.S1.002 1 研究背景 在工程建设中不可避免地遇到下伏溶洞的地基，设计和施工时需特别注意溶洞对地基稳定性的 ［1］ 影响，若其处治不当可能给工程带来极大的安全隐患和经济损失 。目前对于岩溶区下伏溶洞地基承 载力问题的研究主要从试验研究、理论研究和数值分析三个方面进行。在试验研究方面，Kiyosumi ［2］ 等 对坚硬地层中空洞上方的条形基础地基承载力进行了模型试验，并对试验过程中的地基破坏特征 ［3］ 进行了描述和分析；刘庭金等 对黏土下伏空洞地基的破坏过程进行了室内模型试验，并与有限元模 ［4］ 拟的结果进行了对比。在理论研究方面，刘之葵等 从弹性理论出发，利用莫尔-库仑屈服准则，对 ［5］ 土洞周围土体稳定性进行计算判别；刘辉等 利用极限分析上限法推导出圆形空洞上方条形基础的地 ［6］ 基极限承载力解析公式。在数值分析研究方面，胡庆国等 采用ANSYS对溶洞上方条形基础地基极 ［7］ 限承载力进行计算，得到了不同溶洞跨度、顶板厚度时的地基极限承载力；阳军生等 采用ABAQUS ［8］ 对岩溶上方圆形基础地基极限承载力进行了计算，并对其影响因素进行了分析；彭芳乐等 利用 PLAXIS研究了空洞的位置、形状、大小以及地基土的类型等诸因素对浅基础承载力与沉降的影响。 然而，上述数值研究方法对地基破坏时的极限状态无法直接模拟，对地基极限承载力的确定主要是 依据特征点的荷载-位移曲线（P-S曲线）进行人为判别，而不能直接得到。 ［9］ 为此，赵明华等 采用一种新型岩土工程数值方法——上限分析有限元法，对考虑地下空洞影响 的岩溶区地基极限承载力进行了计算和分析，由于上限分析有限元法计算误差的主要来源于网格离 散误差，因此，本文在文献［9］的基础上，引入一种网格自适应的上限分析有限元法，该方法无需人 为划分网格，从而避免由网格划分带来的计算误差，同时还能在网格单元数相同的条件下，获得更 高的计算精度和计算效率，并针对文献［9］对下伏溶洞地基