求解地基极限承载力上限有限元滑动面搜索法.docVIP

求解地基极限承载力上限有限元滑动面搜索法 　　摘 要：提出一种基于上限有限元的滑动面搜索法用于求解地基极限承载力。计算思路为：通过将破坏区域按照假定可能的破坏模式的方式进行网格划分并且网格参数化，然后利用上限有限元求解对应参数条件下的上限解，再从中选取最优上限解及其对应的滑动面，即达到利用上限有限元搜索最优滑动面的目的。该法综合了上限有限元计算灵活、刚体滑块上限法可获得直观滑动面的特点，可高效快速的获得一系列上限解并从中优选，具有计算成本低，且能保证计算精度等优势。通过条形基础地基极限承载力的计算分析，验证了该法的有效性，并可望应用于类似问题。 　　 关键词：极限分析上限有限元上限法刚体滑动面地基极限承载力 　　 　　 1 引言 　　在进行岩土工程稳定性分析时，采用极限分析上限法能直接快速获得上限解。传统的上限法（刚体滑块上限法）需预先假定破坏模式[1]，而上限有限元法通过将破坏区域离散成三角形单元，通过系统耗能最小原理直接搜索获得上限解[2,3]。 　　刚体滑块上限法对于特定问题求解简单，并可得到具体的滑动面；上限有限元法计算结果与网格单元存在依赖性，计算获得的滑动面形状不明显。 　　本文试图结合两者的优势，即采用假定破坏模式的方式划分上限有限元网格，将破坏区域的滑动面表示成带参数的多段线，调整滑动面参数并利用上限有限元进行计算获得一系列上限解，最后从中选取最优解。 　　本文以经典的条形基础地基承载力问题为例，说明采用上述思路的实现过程，计算结果与现有文献进行对比分析，论证基于上限有限元的滑动面搜索法对特定问题的可行性和有效性。 　　 2 极限分析上限有限元 　　 极限分析上、下限理论应用于岩土工程稳定性分析的原理和方法在Chen W F的著作中有详细的论述[1]。由于计算简便，极限分析上限法长期以来均采用与极限平衡法类似的假定刚体滑面模型的形式。 　　 然而，当破坏模式不易假定时，采用Sloan等提出的上限有限元模型就变得更加有效[2]。本文所采用的上限有限元原理和流程参考了Sloan等所作的工作[2,3]，区别在于网格划分按照假定???坏模式的方法进行。 　　 当采用线性规划模型时，上限有限元需对摩尔-库伦屈服准则进行线性化，之后模型即可转化为如下线性规划问题： 　　Minimize(1) 　　Subject to(2) 　　 式(1)、(2)中为目标函数(1)的系数向量；为决策变量，由单元节点速度（实域）、速度间断线辅助速度参数（非负值）和单元内部塑性乘子（非负值）组成；为等式约束系数矩阵，为等式约束右侧向量。 　　 由虚功率平衡方程获得目标函数(1)；等式 (2)由单元内部和速度间断线塑性流动约束条件、速度和应力边界约束条件组成。上限有限元基本原理及具体实现过程参见文献[3]。 　　 3 条形基础地基承载力 　　 3.1 上限有限元模型建立 　　 地基极限承载力可方便的转化为求解承载力系数,和，其中与自重相关的承载力系数的解答需通过数值计算获得。为减少篇幅，本文仅进行的计算。 　　 地基承载力上限有限元模型网格划分如图1所示。利用对称性，只考虑模型右侧的一半。可以看出，三角形单元的布置与现有的破坏模式近似。破坏模式由5个三角形组成，其中过渡区三角形个数；破坏模式可由图中的角度,,,唯一确定。为使刚体运动更加灵活，过渡区中每个三角形的角度,可取不同数值。以角度为决策变量，耗散能最小化为目标函数，对应的速度矢量闭合图为几何约束条件，即构成刚体滑块上限法的求解方法。 　　 刚体滑块上限法认为耗散能仅在间断线上发生[4]。采用上限有限元时，三角形单元也允许发生塑性变形，而不再是刚体。如图1，上限有限元求解承载力问题时，建立图示的坐标系，模型边界条件为：基础下方三角形单元左侧边节点向速度分量，上边向速度分量，向速度分量；下边界滑动面对应的边界条件（刚性边界）为，此约束并未将临近的三角形单元节点速度置零，而是在刚性边界添加虚拟节点，并将这些节点的速度分量均置零，然后在滑动面上施加间断线约束条件。 　　 　　图1 地基承载力上限有限元模型网格划分及边界条件 　　 由于基础下方的三角形向速度置零，意味着基础与地基之间完全粗糙。 　　 上限有限元计算前需确定三角形单元的数目()和摩尔-库伦屈服准则线性化对应的塑性乘子数目，之后将模型的几何和力学参数导入已编制的上限有限元程序，通过求解线性规划问题即可获得承载力系数的上限解。 　　 3.2 地基破坏时滑动面搜索策略 　　 按照图1模型设置的三角形单元较少，因此对于特定网格所获上限解精度不高。于是，将上述上限有限元的网格参数化，即由三个角度参数,,确定一系列的网格，计算对应的上限解，再从中选取最小值，与其对应的,,所确定的破坏模式和滑