基于有限单元与极限平衡的土体稳定性评价方法技术

本发明专利技术属于岩土工程设计与施工中的土体稳定性评价技术领域，具体涉及一种基于有限单元与极限平衡的土体稳定性评价方法专利申请事宜。所述土体稳定性包括边坡稳定性、地基承载力和支挡结构土压力三类问题，具体包括：1）地形测绘、地质勘察与有限单元模型构建，2）总体约束方程构建，3）自变量，4）目标函数，5）评价与结论等步骤。本申请以有限单元与极限平衡理论为基础，能够适应土体性质的不连续、不均匀、在破坏时呈现的大变形等特性，且不需要引入最危险滑裂面假定、并能够同时适用于土体稳定性评价中的边坡稳定性、地基承载力和支挡结构土压力三类问题，其评价方法在理论上更为严格，结论与工程实践也更为符合。

【技术实现步骤摘要】
基于有限单元与极限平衡的土体稳定性评价方法
本专利技术属于岩土工程设计与施工中的土体稳定性评价
，具体涉及一种基于有限单元与极限平衡的土体稳定性评价方法专利申请事宜。
技术介绍
土体稳定性评价包括三类问题：边坡稳定性、地基承载力和支挡结构土压力。对于这三类问题，Columb和Rankine的方法被广泛应用到边坡稳定性和地基承载力评价中，属于一种极限分析的方法。在土体稳定性评价中，尽管还有有限元等其它方法，但极限平衡法因其理论简单、概念清晰、求解方便等优点而被广泛使用。经典极限平衡法的计算步骤可归纳为两步：首先是假定一个滑裂面，计算该面发生失稳的安全系数；然后是寻找最小的安全系数，也就是寻找临界滑裂面。为了对土体的稳定性进行评价，首先需要假定土体因某种扰动因素从目前的稳定状态进入极限平衡状态，此时，土体内出现一个假想的滑裂面，在该滑裂面上，每一点的法向应力和切向应力都满足Mohr-Columb强度准则。Fellenius于1927年提出了边坡稳定分析的圆弧滑动分析法，也就是瑞典圆弧法。该方法中，滑裂面假定为圆弧，土条底面的法向应力被看作是土条重力在土条底面法线方向的投影，因此对圆心取矩时不出现，使计算工作量大大简化。但该法忽略了条间力的作用，不能满足所有静力平衡条件。当滑裂面不是圆弧时，该法不再适用。Bishop法在瑞典圆弧法的基础上进行了发展，该法仍然假定滑裂面为圆弧，但考虑了条间力的影响。通过土条的竖向力平衡及力矩平衡推求出安全系数。Bishop法考虑了条间力，比瑞典圆弧法前进了一大步，对于均质土坡，目前普遍认为该法可以得到比较合理的解。但是，滑裂面为圆弧的假定限制了该法的应用范围。Janbu法假定滑裂面为任意形状，且假定条间力的作用点位置（即推力线位置）为己知。对每一土条，可列出两个静力平衡方程和一个力矩平衡方程，属于严格方法。但其推力线的假定必须符合条间力的合理性要求。尽管该法的滑裂面为任意形状，但必须事先假定，因此该法仅适用于滑裂面位置、形态均基本确定的边坡稳定分析。不平衡推力传递法是按照折线滑动面将滑动土体分块，假定块间力的合力与上一土块底面平行。该法对滑体仅进行分块，而没有分条，因而计算工作量大大减小。同样，该法也仅适用于滑裂面位置、形态均基本确定的边坡稳定性分析。上述所列仅是极限平衡方法范畴中用于边坡稳定性评价的几种典型方法。总体上，由上述方法可知：在经典极限平衡法中，无论是假定滑裂面为圆弧形、还是任意形状，也无论是严格方法、还是非严格方法，都必须首先假定滑裂面的形态（如Fellenius法和Bishop法），或者同时假定滑裂面的形态和位置（如Janbu法和不平衡推力传递法）。但这种假定在实际应用中并不完全准确，使得土体极限平衡分析法在实际应用中受到一定限制。随着计算机及工程实践的发展，更多的工程尝试采用有限元法对土体稳定性进行分析，与传统的极限平衡法相比，有限元法的优点有滑裂面的形状或位置不需要事先假定，并能够提供应力与变形的全部信息，理论体系相对严密，但有限元法需要引入土体中最不易弄清的本构关系，使得有限元法在实际应用时受到了很大限制。因此，设计一种不需要引入最危险滑裂面假定、不需要引入土体中最不易弄清的本构关系的土体稳定性评价方法，对岩土工程中边坡稳定性、地基承载力和支挡结构土压力工程的设计与施工都有着重要意义。
技术实现思路
本申请目的在于提供一种基于有限单元与极限平衡的土体稳定性评价方法，该方法不需要引入最危险滑裂面假定、不需要引入土体中最不易弄清的本构关系、并能够同时适用于土体稳定性评价中的边坡稳定性、地基承载力和支挡结构土压力三类问题，其评价方法在理论上更为严格，结论与工程实践也更为符合。本申请所采取的技术方案详述如下。一种基于有限单元与极限平衡的土体稳定性评价方法，所述土体稳定性包括边坡稳定性、地基承载力和支挡结构土压力三类问题，所述方法具体包括以下步骤：（一）地形测绘、地质勘察与有限单元模型构建针对待评价土体，进行地形测绘、地质勘察和土体物理力学参数测定，并建立土体稳定性评价的有限单元模型；需要简要解释强调的是：地形测绘应覆盖土体稳定性评价的范围；地质勘察应包括：查明评价范围内土体的地下水位线、具有不同物理力学性质的各类土体的空间分布等，并为评价范围内各类土体物理力学参数的测定提供必要的试样试件；需要测定的土体物理力学参数包括：土体内摩擦角、粘滞系数、天然容重、干容重、饱和容重、含水率等；根据各类土体的空间分布、地下水位线以及工程实践对计算精度的需求，将评价范围内的土体划分为有限个单元，形成土体稳定性评价的有限单元模型；实际划分中：对平面二维问题，可采用三角形和/或多边形单元进行划分；对空间三维问题，可采用四面体和/或多面体单元进行划分；（二）总体约束方程构建需要说明的是，为了便于直接采用弹性力学中的有关结论，本申请中节点应力的大小与方向均采用弹性力学中的有关规定，即正应力以拉应力为正，切应力以使微分体产生逆时针转动趋势为正；所述总体约束方程由屈服准则约束方程、节点平衡约束方程和应力已知边界条件约束方程共同构成；具体而言：（2.1）构建屈服准则约束方程根据测定出的土体物理力学参数及相应的土体强度理论，建立评价范围内各类土体以应力形式表达的屈服准则；以平面问题、Mohr-Columb屈服准则为例，设有限单元模型的节点编号依次为，第个节点所对应的节点应力分量分别为、和，该点在某一方向上的正应力和切应力分别为和，则有：，公式(1)；，公式(2)；进一步，设测定出的土体内摩擦角为，并记，根据Mohr-Columb屈服准则，其节点应力在任一方向上的正应力和切应力均应满足Mohr-Columb准则，可得到：，公式(3)；而公式(3)等价于，公式(4)；若记为：，公式(5)；，公式(6)；则可得到公式(3)等价于；进一步将公式(1)与公式(2)代入公式(5)的左边表达式，则可得到：，公式(7)；公式(7)中：；根据公式(7)，可得到公式(5)等价于：，公式(8)；根据Mohr-Columb准则，对任意的QUOTE，公式(8)均应成立，进一步可得到：，公式(9)；由于土体不能承受拉应力，因此可限定及，此时可得到公式(9)等价于，公式(10)；据上可得到，在限定及的条件下，公式(5)等价于公式(10)；类似地，将式(1)与式(2)代入公式(6)的左边表达式，并限定及，可得到公式(6)等价于，公式(11)；综合公式(10)与公式(11)，对平面问题，Mohr-Columb准则可综合地写为，公式(12)公式(12)中，；也即，公式(12)即为本申请的土体稳定性评价方法的屈服准则约束方程；需要说明的是，在公式(12)中，对于每个节点，Mohr-Columb准则有4个不等式约束方程；若节点总数为，则共有个屈服准则不等式约束方程；对三维空间问题、或其它以应力形式表达的屈服准则问题，均可采用上述方法，类似地建立屈服准则约束方程；（2.2）节点静力平衡约束方程构建对于平面问题，单个节点的静力平衡方程只包括和两个方向的平衡方程；对于三维空间问题，单个节点的静力平衡方程则包括、和三个方向的平衡方程；以平面问题为例，可采用如下方法建立节点静力平衡约束方程：如图1所示，设围绕节点的所有单元按逆时针本文档来自技高网...

【技术保护点】
1.一种基于有限单元与极限平衡的土体稳定性评价方法，其特征在于，所述土体稳定性包括边坡稳定性、地基承载力和支挡结构土压力三类问题，所述方法具体包括以下步骤：（一）地形测绘、地质勘察与有限单元模型构建针对待评价土体，进行地形测绘、地质勘察和土体物理力学参数测定，根据各类土体的空间分布、地下水位线以及工程实践对计算精度的要求，将评价范围内的土体划分为有限个单元，建立土体稳定性评价的有限单元模型；实际划分中：对平面二维问题，采用三角形和/或多边形单元进行划分；对空间三维问题，采用四面体和/或多面体单元进行划分；（二）总体约束方程构建所述总体约束方程由屈服准则约束方程、静力平衡约束方程和应力已知边界条件约束方程共同构成；（三）自变量以所有节点的应力分量为自变量；（四）目标函数对边坡稳定性问题，目标函数是判断满足步骤（二）的总体约束方程的自变量是否存在；若存在，判定该土体为稳定的；若不存在，判定该土体为不稳定的；对地基承载力问题，目标函数是判断满足步骤（二）的总体约束方程的自变量是否存在；若存在，判定该地基承载力满足要求；若不存在，判定该地基承载力不满足要求；对支挡结构的土压力问题，设待求边界上以节点应力为变量表示的土压力的大小为P；对极限最小土压力，目标函数是求解满足步骤（二）的总体约束方程条件下的P的最小值；对极限最大土压力，目标函数是求解满足步骤（二）的总体约束方程条件下的P的最大值；（五）评价与结论由步骤（二）的总体约束方程、步骤（三）的自变量、步骤（四）的目标函数共同构成稳定性评价方法的数学模型；采用最优化求解的方法，对上述数学模型进行计算求解；对边坡稳定问题，是判定数学模型是否有解；若数学模型的解存在，判定该土体为稳定的；若解不存在，判定该土体为不稳定的；对地基承载力问题，是判定数学模型是否有解；若数学模型的解存在，判定该地基承载力满足要求；若解不存在，判定该地基承载力不满足要求；对支挡结构土压力问题，是求出土压力P的最小或最大值；P的最小值即为极限最小土压力，P的最大值即为极限最大土压力。...

【技术特征摘要】
1.一种基于有限单元与极限平衡的土体稳定性评价方法，其特征在于，所述土体稳定性包括边坡稳定性、地基承载力和支挡结构土压力三类问题，所述方法具体包括以下步骤：（一）地形测绘、地质勘察与有限单元模型构建针对待评价土体，进行地形测绘、地质勘察和土体物理力学参数测定，根据各类土体的空间分布、地下水位线以及工程实践对计算精度的要求，将评价范围内的土体划分为有限个单元，建立土体稳定性评价的有限单元模型；实际划分中：对平面二维问题，采用三角形和/或多边形单元进行划分；对空间三维问题，采用四面体和/或多面体单元进行划分；（二）总体约束方程构建所述总体约束方程由屈服准则约束方程、静力平衡约束方程和应力已知边界条件约束方程共同构成；（三）自变量以所有节点的应力分量为自变量；（四）目标函数对边坡稳定性问题，目标函数是判断满足步骤（二）的总体约束方程的自变量是否存在；若存在，判定该土体为稳定的；若不存在，判定该土体为不稳定的；对地基承载力问题，目标函数是判断满足步骤（二）的总体约束方程的自变量是否存在；若存在，判定该地基承载力满足要求；若不存在，判定该地基承载力不满足要求；对支挡结构的土压力问题，设待求边界上以节点应力为变量表示的土压力的大小为P；对极限最小土压力，目标函数是求解满足步骤（二）的总体约束方程条件下的P的最小值；对极限最大土压力，目标函数是求解满足步骤（二）的总体约束方程条件下的P的最大值；（五）评价与结论由步骤（二）的总体约束方程、步骤（三）的自变量、步骤（四）的目标函数共同构成稳定性评价方法的数学模型；采用最优化求解的方法，对上述数学模型进行计算求解；对边坡稳定问题，是判定数学模型是否有解；若数学模型的解存在，判定该土体为稳定的；若解不存在，判定该土体为不稳定的；对地基承载力问题，是判定数学模型是否有解；若数学模型的解存在，判定该地基承载力满足要求；若解不存在，判定该地基承载力不满足要求；对支挡结构土压力问题，是求出土压力P的最小或最大值；P的最小值即为极限最小土压力，P的最大值即为极限最大土压力。2.如权利要求1所述基于有限单元与极限平衡的土体稳定性评价方法，其特征在于，步骤（一）中，地形测绘应覆盖土体稳定性评价的范围；地质勘察应包括：查明评价范围内土体的地下水位线、具有不同物理力学性质的各类土体的空间分布等，并为评价范围内各类土体物理力学参数的测定提供必要的试样试件；需要测定的土体物理力学参数包括：土体...

【专利技术属性】
技术研发人员：张金良，尹德文，吴昊，李希，
申请(专利权)人：黄河勘测规划设计研究院有限公司，
类型：发明
国别省市：河南,41