发明名称 基于全局最优化的边坡稳定性分析条分法
摘要 本发明公开了一种基于全局最优化的边坡稳定性分析条分法,属于边坡稳定性分析领域,本发明将边坡滑体划分为竖直条块或斜条块,将条块底面的作用力(包括法向力、剪切力)以及条块间的作用力(包括法向力、剪切力)作为未知变量,以边坡的稳定安全系数作为目标函数,结合条块的平衡方程、条块底滑面的屈服条件、条间接触面的屈服条件等约束条件,建立边坡稳定性分析的非线性数学规划模型,并使用优化算法求解全局最优的最大稳定安全系数。使用本发明方法计算得到边坡安全系数具有下限解的性质。本发明方法具有概念明确、计算精度高等特点,可将其应用于土质边坡或岩质边坡的稳定性分析。
申请公布号 CN103485353B 申请公布日期 2015.08.12
申请号 CN201310436127.3 申请日期 2013.09.24
申请人 昆明理工大学 发明人 李泽;魏久坤;武亮
分类号 E02D17/20(2006.01)I 主分类号 E02D17/20(2006.01)I
代理机构 代理人
主权项 一种基于全局最优化的边坡稳定性分析条分法,其特征在于:将边坡滑体划分为垂直条块或斜条块,基于潘家铮最大值原理,将条块底面的作用力以及条块间的作用力作为优化变量,以边坡的稳定安全系数作为目标函数,结合条块的平衡方程、条块底滑面的屈服条件、条间接触面的屈服条件,建立边坡稳定性分析的非线性数学规划模型,并使用数学规划优化算法求解边坡稳定安全系数的最大值;上述方法按以下步骤进行:(1)确定边坡的计算参数根据边坡的实际情况,确定计算参数为:几何参数、材料参数、荷载参数,其中材料参数包括容重、凝聚力、摩擦角;(2)将边坡滑体划分为条块将边坡滑面以上的滑体划分为垂直条块或斜条块,条块数量为<i>n</i>;(3)建立求解边坡安全系数的条分法非线性数学规划模型目标函数为边坡稳定的安全系数,以条块底滑面的法向力、剪切力以及条间接触面的法向力、剪切力为优化变量,以边坡稳定安全系数为目标函数,约束条件为条块的平衡条件、条块底滑面的屈服条件、条间接触面的屈服条件,形成求解边坡最大安全系数的非线性数学规划模型:A、边坡条块的静力受力分析边坡滑体划分为条块,条块数量为<i>n</i>,取其中一个条块<i>i</i>进行受力分析,建立<img file="dest_path_image001.GIF" wi="29" he="18" />坐标系,竖直向下为<img file="dest_path_image002.GIF" wi="16" he="18" />轴正向,滑体滑动水平方向为<img file="dest_path_image003.GIF" wi="14" he="16" />正向;条块<i>i</i>的宽度为<img file="dest_path_image004.GIF" wi="29" he="25" />,底滑面的长度为<img file="dest_path_image005.GIF" wi="13" he="25" />,底滑面的倾角为<img file="dest_path_image006.GIF" wi="18" he="25" />,条间交界面的长度为<img file="dest_path_image007.GIF" wi="16" he="25" />,条块<i>i</i>的形心为<img file="dest_path_image008.GIF" wi="20" he="25" />点;条块<i>i</i>的形心<img file="374264dest_path_image008.GIF" wi="20" he="25" />上作用有水平力<img file="dest_path_image009.GIF" wi="20" he="25" />、竖向力<img file="dest_path_image010.GIF" wi="20" he="25" />;条块底滑面形心上作用为法向力<img file="dest_path_image011.GIF" wi="21" he="25" />、抗剪力<img file="dest_path_image012.GIF" wi="16" he="25" />;条块左侧面形心上作用为法向力<img file="dest_path_image013.GIF" wi="20" he="25" />、剪切力<img file="dest_path_image014.GIF" wi="21" he="25" />;条块右侧面形心上作用为法向力<img file="dest_path_image015.GIF" wi="29" he="25" />、剪切力<img file="dest_path_image016.GIF" wi="32" he="25" />;其中底面或条间法向力的方向规定为拉正压负;B、将边坡底滑面的强度储备系数作为目标函数,并求其最大值,强度储备系数为<img file="dest_path_image017.GIF" wi="141" he="21" />,其中<img file="dest_path_image018.GIF" wi="13" he="14" />,<img file="dest_path_image019.GIF" wi="14" he="17" />分别为边坡滑面原始的凝聚力和摩擦角,<img file="dest_path_image020.GIF" wi="16" he="18" />,<img file="dest_path_image021.GIF" wi="18" he="21" />分别为边坡滑面进行强度折减以后的凝聚力和摩擦角;C、约束条件约束条件包括:条块的平衡条件、条块底滑面的屈服条件、条间接触面的屈服条件;① 条块的平衡方程条块<i>i</i>,其受到外力、底滑面上以及左右侧面上的作用力而保持平衡,其平衡方程为:水平方向力的平衡方程:<img file="dest_path_image022.GIF" wi="232" he="25" />,<img file="dest_path_image023.GIF" wi="81" he="22" />,竖直方向力的平衡方程:<img file="dest_path_image024.GIF" wi="236" he="25" />,<img file="475425dest_path_image023.GIF" wi="81" he="22" />;② 条块底滑面的屈服条件当滑体处于滑动临界状态时,滑面应满足Mohr‑Coulomb屈服条件式,土条底滑面上的屈服条件可以写为:<img file="dest_path_image025.GIF" wi="134" he="25" />,<img file="317479dest_path_image023.GIF" wi="81" he="22" />,求解强度储备系数<i>K</i>时,将<img file="775006dest_path_image017.GIF" wi="141" he="21" />带入上式可得:<img file="dest_path_image026.GIF" wi="181" he="25" />,<img file="651695dest_path_image023.GIF" wi="81" he="22" />;③ 条间接触面的屈服条件当滑体处于滑动临界状态时,条间接触面应满足Mohr‑Coulomb屈服条件式,条间接触面上的屈服条件可以写为:<img file="dest_path_image027.GIF" wi="148" he="26" />,<img file="dest_path_image028.GIF" wi="100" he="22" />,求解强度储备系数<i>K</i>时,将<img file="100256dest_path_image017.GIF" wi="141" he="21" />带入上式可得:<img file="dest_path_image029.GIF" wi="196" he="26" />,<img file="429606dest_path_image028.GIF" wi="100" he="22" />;D、求解边坡最大安全系数的非线性数学规划模型求解边坡稳定性的非线性数学规划模型以安全系数为目标函数,以平衡条件、底滑面的屈服条件和条间接触面的屈服条件,数学模型具体表达式为:<img file="dest_path_image030.GIF" wi="320" he="150" />;(4)求解全局最优安全系数:采用序列二次规划法求解全局最优的最大安全系数,以及对应的底滑面和条间的法向力、剪切力。
地址 650093 云南省昆明市五华区学府路253号