CN109753701B - 一种随机渗流场作用下的土质边坡可靠度分析下限法 - Google Patents
一种随机渗流场作用下的土质边坡可靠度分析下限法 Download PDFInfo
- Publication number
- CN109753701B CN109753701B CN201811568164.9A CN201811568164A CN109753701B CN 109753701 B CN109753701 B CN 109753701B CN 201811568164 A CN201811568164 A CN 201811568164A CN 109753701 B CN109753701 B CN 109753701B
- Authority
- CN
- China
- Prior art keywords
- soil slope
- soil
- slope
- random
- random number
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Abstract
本发明公开了一种随机渗流场作用下的土质边坡可靠度分析下限法,本发明方法为:步骤1、拟定土质边坡的计算参数;步骤2、建立土质边坡可靠度计算的极限状态函数;步骤3、生成土质边坡地下水水位的随机数;步骤4、生成土质边坡土体材料凝聚力、摩擦角的随机数;步骤5、使用有限单元离散土质边坡,得到有限单元网格;步骤6、进行土质边坡的随机渗流场计算;步骤7、建立土质边坡可靠度分析下限法线性规划模型;步骤8、求解土质边坡下限法线性规划模型;步骤9、统计土质边坡的可靠度指标。本发明可以获得土质边坡可靠度指标随地下水位随机变化的规律,并获得土质边坡的失效概率。
Description
技术领域
本发明涉及一种随机渗流场作用下的土质边坡可靠度分析下限法,属于边坡工程稳定性分析计算领域。
背景技术
土质边坡的稳定性问题是一个涉及人民生命、财产安全的重要问题。现阶段土质边坡的稳定性分析方法主要包括:确定性分析方法和可靠度分析方法,随着计算理论和计算机技术的巨大进步,土质边坡可靠度分析方法越来越受到重视。土质边坡的失稳破坏主要与两个因素息息相关,一是土体材料抗剪参数的不确定性,即土体材料凝聚力、内摩擦角均存在变异性;另一个是地下水位的随机变化,地下水位的随机变化导致土质边坡的渗流场是随机分布的,进而土质边坡内部各个点的孔隙水压力也是随机变化的。土质边坡的稳定性是一个与土体材料凝聚力、摩擦角和孔隙水压力变异性相关的随机问题。鉴于此,本发明基于国家自然科学基金项目(51564026)进行研究工作,提出一种能够同时考虑土质边坡地下水位随机变化和抗剪参数变异性的土质边坡可靠度分析方法。
发明内容
本发明提供了一种随机渗流场作用下的土质边坡可靠度分析下限法,以获得土质边坡的可靠度指标,为土质边坡可靠度分析提供一种新的计算方法。
本发明的技术方案是:一种随机渗流场作用下的土质边坡可靠度分析下限法,所述方法具体步骤如下:
步骤1、拟定土质边坡的计算参数;
步骤2、建立土质边坡可靠度计算的极限状态函数;
步骤3、生成土质边坡地下水水位的随机数;
步骤4、生成土质边坡土体材料凝聚力、摩擦角的随机数;
步骤5、使用有限单元离散土质边坡,得到有限单元网格;
步骤6、进行土质边坡的随机渗流场计算;
步骤7、建立土质边坡可靠度分析下限法线性规划模型;
步骤8、求解土质边坡下限法线性规划模型;
步骤9、统计土质边坡的可靠度指标。
所述土质边坡的计算参数包括土质边坡的几何参数,土质边坡的地层条件,土质边坡土体材料的凝聚力、摩擦角、容重和渗透系数,土质边坡地下水水位的信息。
所述土质边坡可靠度计算的极限状态函数为:
式中:Z是土质边坡的极限状态函数,Wn是土质边坡地下水水位随机数的数量,N是土质边坡土体材料凝聚力和摩擦角随机数的数量,kγ(tw,tm)是第tw个地水水水位作用下与第tm个土体材料凝聚力、摩擦角的随机数相关的容重超载系数的随机数;是与第tm个土体材料凝聚力、摩擦角随机数相关的极限容重随机数;γa是土质边坡土体材料实际的容重,cr(tm)是土质边坡土体材料凝聚力的第tm个随机数,是土质边坡土体材料摩擦角的第tm个随机数。
所述步骤3具体为:土质边坡的地下水水位在最高地下水水位与最低地下水水位之间随机变化,假设土质边坡地下水水位的变化规律符合截尾正态分布,土质边坡地下水水位的随机数按下式生成:
式中:tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的总数量,Hw(tw)是土质边坡地下水水位的第tw个随机数,μw是土质边坡地下水水位的均值,σw是土质边坡地下水水位的标准差,Random是正态分布随机数生成函数,Normal表示随机数符合正态分布;Hlb是土质边坡地下水水位的下界,Hub是土质边坡地下水水位的上界。
所述步骤4包括:
①假设土质边坡土体材料的凝聚力和摩擦角符合对数正态分布;
②确定土质边坡材料的凝聚力、摩擦角的均值、标准差;
③土质边坡土体材料凝聚力的随机量和摩擦角的随机量的随机数按下式生成:
式中:tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量,cr(tm)是土质边坡土体材料凝聚力的第tm个随机数,是土质边坡土体材料摩擦角的第tm个随机数,μc是土质边坡土体材料凝聚力的均值,是土质边坡土体材料摩擦角的均值,σc是土质边坡土体材料凝聚力的标准差,是土质边坡土体材料摩擦角的标准差,Random是正态分布随机数生成函数,lognormal表示随机数符合对数正态分布。
所述步骤5具体为:使用三角形三结点线性有限单元离散土质边坡,在总体坐标系(x,y)下,任意一个有限单元e中有三个结点,第e个单元的第i个结点具有x向有效正应力y向有效正应力和切应力同时,第e个单元的第i个结点作用有孔隙水压力其中,i=1,…,3。
所述步骤6具体为:包括:
①基于步骤5得到的有限单元网格,根据稳定渗流理论进行土质边坡的随机渗流场计算,二维稳定渗流公式如下:
式中:kx是土体材料x方向的渗透系数,ky是土体材料y方向的渗透系数,Hr是土质边坡内各点的随机水头函数;
②将土质边坡地下水水位随机数Hw(tw)作为已知边界条件带入上式,求解上式可获得所有有限单元结点的孔隙水压力值:其中:tw=(1,…,Wn),e=(1,…,Ne),Wn是土质边坡地下水水位随机数的数量,Ne是土质边坡中有限单元的数量;是第tw地下水水位随机数Hw(tw)作用下有限单元e中结点1的孔隙水压力值,是第tw个地下水水位随机数Hw(tw)作用下有限单元e中结点2的孔隙水压力值,是第tw个地下水水位随机数Hw(tw)作用下有限单元e中结点3的孔隙水压力值。
所述步骤7以有限单元结点的应力为决策变量建立求解土质边坡可靠的的下限法线性规划模型,具体为:①建立土质边坡有限单元的平衡方程;②建立土质边坡有限单元公共边的应力连续条件;③建立土质边坡有限单元的屈服条件;④建立土质边坡有限单元的静力边界条件;⑤建立目标函数;⑥根据土质边坡的极限状态函数、①、②、③、④和⑤建立的函数建立求解土质边坡可靠的下限法线性规划模型;
具体为:
①建立土质边坡有限单元的平衡方程为:
式中:e=(1,…,Ne),Ne是土质边坡中有限单元的数量;tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的数量;tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量;kγ(tw,tm)是第tw个地下水水位作用下与第tm个土体材料凝聚力、摩擦角的随机数相关的容重超载系数的随机数; 是分别是线性三角形有限单元e的6个形函数系数;分别是有限单元e的第1、2、3个结点沿x方向的有效正应力;分别是有限单元e的第1、2、3个结点沿y方向的有效正应力;分别是有限单元e的第1、2、3个结点的切应力;分别是第tw个地下水水位作用下有限单元e的第1、2、3个结点的孔隙水压力;Ae是有限单元e的面积;γe是有限单元e的容重;
②建立土质边坡有限单元公共边的应力连续条件:
式中:d=(1,…,Nd),Nd是土质边坡中有限单元公共边的数量; θd是公共边的倾角,逆时针为正;分别是公共边d的第1、2、3、4个结点沿x方向的有效正应力;分别是公共边d的第1、2、3、4个结点沿y方向的有效正应力;分别是公共边d的第1、2、3、4个结点的切应力;
③建立土质边坡有限单元的屈服条件:
式中:n=(1,…,Nn),Nn是土质边坡中有限单元结点的总数量;σn=[σ'xn,σ'yn,τxyn]T; Mp土体材料摩尔库伦屈服条件线性化正多边形的边数;σ'xn是第n个结点沿x方向的有效正应力;σ'yn是第n个结点沿y方向的有效正应力;τxyn是第n个结点的切应力;是土质边坡土体材料摩擦角的第tm个随机数;cr(tm)是土质边坡土体材料凝聚力的第tm个随机数;
④建立土质边坡有限单元的静力边界条件:
式中:b=(1,…,Nb),Nb是土质边坡中边界上的有限单元的数量; 分别是边界单元b中靠边界上的第1、2个结点的沿x方向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点的沿x方向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点的切应力;分别是边界单元b中靠边界上的第1、2个结点已知的沿法向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点已知的沿切向的切应力;分别是第tw个地下水水水位作用下边界单元b中靠边界上的第1、2个结点的孔隙水压力;θb是边界单元b靠边界上的第1、2个结点连线的倾角;
⑤目标函数:
在考虑土体的容重超载时,下限法目标函数为:
Maximize:kγ(tw,tm)
⑥建立求解土质边坡可靠度的下限法线性规划模型:
式中:Z是土质边坡的极限状态函数。
所述步骤8具体为:
S1、将土质边坡地下水水位随机数Hw(tw)从tw=1到tw=Wn循环带入渗流场计算公式,逐次计算土质边坡的Wn个渗流场,获得土质边坡中每个有限单元结点的孔隙水压力值:其中:tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的数量;e=(1,…,Ne),Ne是土质边坡中有限单元的数量;
S3、在tw=1到tw=Wn的每次循环中将cr(tm)、从tm=1到tm=N循环,将N个土体材料的凝聚力和摩擦角的随机数代入求解土质边坡可靠度的下限法线性规划模型,使用“单纯形”优化算法求解土质边坡可靠度的下限法线性规划模型,求解得到Wn×N个容重超载系数的随机数[kγ(tw,tm),(tw=1,…,Wn),tm=(1,…,N)]和对应的有限单元的正应力和切应力;同时,在tw=1到tw=Wn的每次循环中使用“二分法”迭代求解土质边坡的Wn×N个强度储备系数[km(tw,tm),(tw=1,…,Wn),tm=(1,…,N)]和与每个强度储备系数对应的有限单元的正应力和切应力。
所述步骤9具体为:
①、计算第tw个随机地下水水位随机数Hw(tw)作用下土质边坡N个强度储备系数的均值:
式中:tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的数量;tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量;是第tw个随机地下水水位随机数Hw(tw)作用下土质边坡的N个强度储备系数的平均值;
②、计算第tw个随机地下水水位随机数Hw(tw)作用下土质边坡N个强度储备系数的标准差:
③、计算土质边坡Wn×N个强度储备系数的均值:
式中:μkm是土质边坡Wn×N个强度储备系数的平均值;
④、计算土质边坡在所有随机地下水水位作用下的Wn×N个强度储备系数的标准差:
式中:σkm是Wn×N个强度储备系数的标准差;
⑤、计算土质边坡的失效概率
土质边坡的失效功能函数如下:
土质边坡的失效概率按下式计算:
式中:Pf是土质边坡的失效概率。
本发明的有益效果是:本发明以土质边坡为研究对象,考虑土质边坡地下水位的随机变化,将塑性极限分析下限定理、有限元离散技术、线性规划方法结合起来,以有限元结点的应力为决策变量建立土质边坡可靠度计算的下限法数学规划模型;通过该模型可以获得土质边坡可靠度指标随地下水位随机变化的规律,并获得土质边坡的失效概率。
附图说明
图1本发明的技术路线图;
图2有限单元示意图(标注三角形有限单元速度模式和孔隙水压力);
图3有限单元示意图(标注三角形有限单元不共节点模式);
图4下限法相邻有限单元之间的公共边示意图;
图5实施例1土质边坡的几何形状示意图(单位:m);
图6实施例1土质边坡随机地下水位直方图(单位:m);
图7实施例1土质边坡土体材料凝聚力随随机数分布图(单位:kPa);
图8实施例1土质边坡土体材料摩擦角随随机数分布图(单位:°);
图9实施例1土质边坡有限单元离散示意图;
图10实施例1土质边坡稳定渗流场孔隙水压力等值线(tw=30)(单位:kPa);
图11实施例1土质边坡tw=1时土质边坡强度储备系数的直方图;
图12实施例1土质边坡tw=30时土质边坡强度储备系数的直方图;
图13实施例1土质边坡tw=50时土质边坡强度储备系数的直方图;
图14实施例1土质边坡在50个地下水位作用下的50×1000个强度储备系数的直方图;
图15实施例1土质边坡tw=1~50时边坡强度储备系数概率密度分布图;
图16实施例1土质边坡tw=5、tm=293时总应力的第一主应力等值线(单位:kPa);
图17实施例1土质边坡tw=5、tm=293时总应力的第三主应力等值线(单位:kPa);
图18实施例1土质边坡强度储备系数的均值随地下水位的变化关系图;
图19实施例1土质边坡强度储备系数的标准差随地下水位的变化关系图;
图20实施例1土质边坡失效概率随地下水位的变化关系图。
具体实施方式
实施例1:如图1-20所示,一种随机渗流场作用下的土质边坡可靠度分析下限法,所述方法具体步骤如下:
步骤1、拟定土质边坡的计算参数;
步骤2、建立土质边坡可靠度计算的极限状态函数;
步骤3、生成土质边坡地下水水位的随机数;
步骤4、生成土质边坡土体材料凝聚力、摩擦角的随机数;
步骤5、使用有限单元离散土质边坡,得到有限单元网格;
步骤6、进行土质边坡的随机渗流场计算;
步骤7、建立土质边坡可靠度分析下限法线性规划模型;
步骤8、求解土质边坡下限法线性规划模型;
步骤9、统计土质边坡的可靠度指标。
流程图如图1所示。
进一步地,可以设置所述方法具体步骤如下:
步骤一、拟定土质边坡的计算参数:根据实施例土质边坡的时间情况确定计算参数,包括:确定土质边坡的几何参数,地层条件,土质边坡土体材料的凝聚力、摩擦角、容重和渗透系数,土质边坡地下水水位的信息。
土质边坡的剖面图如图5所示,具体参数如下:地层条件:土质边坡为均质边坡,土体材料为粘土;边坡几何参数:土质边坡顶宽20.0m、边坡高度10.0m、坡角为45°;土体材料的凝聚力取50.00kPa、摩擦角取10°、容重取18.0kN/m3;土体的渗透系数取7e-8m/s;地下水水位信息:土质边坡地下水位的最低水位是10.0m、地下水位的最高水位是20.0m,即地下水位Hw在10m至20m之间随机变化。
步骤二、建立土质边坡可靠度计算的极限状态函数。
土质边坡可靠度计算的极限状态函数为:
式中:Z是土质边坡的极限状态函数,tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的数量,tm=(1,…,N),N是土质边坡材料凝聚力/摩擦角随机数的数量,kγ(tw,tm)是第tw个地水位作用下与第tm个土体材料的凝聚力、摩擦角的随机数相关的容重超载系数的随机数;是与第tm个土体材料的凝聚力、摩擦角随机数相关的极限容重随机数;γa是边坡土体材料实际的容重,cr(tm)是边坡土体材料凝聚力的第tm个随机数,是边坡土体材料摩擦角的第tm个随机数。
本发明定义土质边坡的强度储备系数如下:
式中:km(tw,tm)是第tw个地下水水位作用下与第tm个抗剪参数的随机数相关的强度储备系数的随机数;c'r(tm)是边坡土体材料的凝聚力进行强度折减以后的第tm个随机数;是边坡土体材料的摩擦角进行强度折减以后的第tm个随机数。
具体的,将土质边坡土体材料的容重设为确定值;将土体材料的凝聚力和摩擦角设为相互独立的随机变量,并假设两者均符合对数正态分布;同时设边坡地下水位为随机变量,并符合正态分布;取边坡地下水位随机数的数量Wn=50,取土体材料凝聚力和摩擦角随机数的数量N=1000,并根据式(1)建立土质边坡可靠度计算的极限状态函数Z。
步骤三、生成土质边坡地下水位的随机数。
土质边坡在地下水位的最低水位至最高水位之间随机变化,本发明假设边坡地下水位符合截尾正态分布,边坡地下水位的随机数按下式生成:
式中:tw=(1,…,Wn),Wn是地下水水位随机数的总数量,Hw(tw)是土质边坡地下水位的第tw个随机数,μw是土质边坡地下水位的均值,σw是土质边坡地下水位的标准差,Random是正态分布随机数生成函数,Normal表示随机数符合正态分布;Hlb是土质边坡地下水位的下界,可以取边坡地下水位的最低水位;Hub是土质边坡地下水位的上界,可以取边坡地下水位的最高水位。
具体的,取边坡地下水位的下界Hlb=10.0m,取边坡地下水位的上界Hub=20.0m,取边坡地下水位随机数的数量Wn=50,取边坡地下水位的均值μw=16,取边坡地下水位的标准差σw=1.79,并根据式(3)生成边坡地下水位的50个随机数,这50个边坡地下水位随机数的直方图如图6所示。
步骤四、生成土质边坡材料的抗剪强度参数的随机数。
包括:①假设土质边坡土体材料的凝聚力和摩擦角符合对数正态分布,②确定土质边坡材料的凝聚力、摩擦角的均值、标准差;③土质边坡材料凝聚力的随机量和摩擦角的随机量的随机数按下式生成:
式中:tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量,cr(tm)是土质边坡土体材料凝聚力的第tm个随机数,是土质边坡土体材料摩擦角的第tm个随机数,μc是土质边坡土体材料凝聚力的均值,是土质边坡土体材料摩擦角的均值,σc是土质边坡土体材料凝聚力的标准差,是土质边坡土体材料摩擦角的标准差,Random是正态分布随机数生成函数,lognormal表示随机数符合对数正态分布。
具体的,设土质边坡土体材料的凝聚力、摩擦角符合对数正态分布;取土体的凝聚力的均值是50.00kPa、凝聚力的标准差是15.0,取土体的摩擦角的均值是10°、摩擦角的标准差是3.0;取土体材料的容重为18.0kN/m3,并设容重为确定值;取土体材料凝聚力和摩擦角随机数的数量N=1000;并根据式(4)分别生成土体的凝聚力、摩擦角的1000个随机数。实施例1土质边坡土体材料凝聚力随随机数分布如图7所示;实施例1土质边坡土体材料摩擦角随随机数分布如图8所示。
步骤五、使用有限单元离散土质边坡。
使用三角形三结点线性有限单元离散土质边坡,有限单元如图2、3所示,在总体坐标系(x,y)下,任意一个有限单元e中有三个结点,第e个单元的第i(i=1,…,3)个结点具有x向有效正应力y向有效正应力和切应力同时,第e个单元的第i(i=1,…,3)个结点作用有孔隙水压力两个相邻有限单元的公共边如图4所示,公共边上有编号为①、②、③、④的4个结点。
具体的,使用有限单元离散实施例土质边坡,土质边坡的有限元网格如图9所示,共离散得到的有限单元Ne=506个、有限单元之间公共边Nd=720条、有限单元之间结点Nn=1518条,每个有限单元结点具有x向有效正应力、y向有效正应力和切应力。
步骤六、进行土质边坡的随机渗流场计算。
包括:①基于步骤5得到的有限单元网格,根据稳定渗流理论进行土质边坡的随机渗流场计算,二维稳定渗流公式如下:
式中:kx是土体材料x方向的渗透系数,ky是土体材料y方向的渗透系数,Hr是土质边坡内各点的随机水头函数。
②将土质边坡地水水水位随机数Hw(tw),tw=(1,…,Wn)作为已知边界条件带入上式,求解上式可获得所有有限单元结点的孔隙水压力值:其中:tw=(1,…,Wn),e=(1,…,Ne),Wn是土质边坡地水位随机数的数量,Ne是土质边坡中有限单元的数量;是第tw地下水水位随机数Hw(tw)作用下有限单元e中结点1的孔隙水压力值,是第tw个地下水水位随机数Hw(tw)作用下有限单元e中结点2的孔隙水压力值,是第tw个地下水水位随机数Hw(tw)作用下有限单元e中结点3的孔隙水压力值。
具体的,将步骤三得到的50个地下水位随机数带入式(5),可获得土质边坡的50个稳定渗流场,并获得50个渗流场的506个有限单元结点的孔隙水压力值:tw=(1,…,50),e=(1,…,506)。图10是实施例1土质边坡第35个稳定渗流场的有限元结点孔隙水压力等值线(tw=35)。
步骤七、建立土质边坡可靠度度分析下限法线性规划模型。
该模型以有限单元结点的x向的有效正应力、y向的有效正应力以及切应力为决策变量,主要步骤包括:①建立土质边坡有限单元的平衡方程;②建立土质边坡有限单元公共边的应力连续条件;③建立土质边坡有限单元的屈服条件;④建立土质边坡有限单元的静力边界条件;⑤建立目标函数;⑥根据土质边坡的极限状态函数、①、②、③、④和⑤建立的函数建立求解土质边坡可靠的下限法线性规划模型。
①至⑥的步骤具体为:
①建立土质边坡有限单元的平衡方程为:
式中:e=(1,…,Ne),Ne是土质边坡中有限单元的数量;tw=(1,…,Wn),Wn是土质边坡地水位随机数的数量;tm=(1,…,N),N是土质边坡土体材料凝聚力和摩擦角随机数的数量;kγ(tw,tm)是第tw个地下水位作用下与第tm个抗剪参数的随机数相关的容重超载系数的随机数; 是分别是线性三角形有限单元e的6个形函数系数;分别是有限单元e的第1、2、3个结点沿x方向的有效正应力;分别是有限单元e的第1、2、3个结点沿y方向的有效正应力;分别是有限单元e的第1、2、3个结点的切应力;分别是第tw个地下水位作用下有限单元e的第1、2、3个结点的孔隙水压力;Ae是有限单元e的面积;γe是有限单元e的容重。
②建立土质边坡有限单元公共边的应力连续条件:
式中:d=(1,…,Nd),Nd是土质边坡中有限单元公共边的数量; θd是公共边的倾角,逆时针为正;分别是公共边d的第1、2、3、4个结点沿x方向的有效正应力;分别是公共边d的第1、2、3、4个结点沿y方向的有效正应力;分别是公共边d的第1、2、3、4个结点的切应力。
③建立土质边坡有限单元的屈服条件:
式中:n=(1,…,Nn),Nn是土质边坡中有限单元结点的总数量;σn=[σ'xn,σ'yn,τxyn]T; Ck=2sin(2kπ/Mp);Mp土体材料摩尔库伦屈服条件线性化正多边形的边数;σ'xn是第n个结点沿x方向的有效正应力;σ'yn是第n个结点沿y方向的有效正应力;τxyn是第n个结点的切应力;是土质边坡材料摩擦角的第tm个随机数;cr(tm)是土质边坡材料凝聚力的第tm个随机数;tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量。
④建立土质边坡有限单元的静力边界条件:
式中:b=(1,…,Nb),Nb是土质边坡中边界上的有限单元的数量; 分别是边界单元b中靠边界上的第1、2个结点的沿x方向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点的沿x方向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点的切应力;分别是边界单元b中靠边界上的第1、2个结点已知的沿法向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点已知的沿切向的切应力;分别是第tw个地下水位作用下边界单元b中靠边界上的第1、2个结点的孔隙水压力;θb是边界单元b靠边界上的第1、2个结点连线的倾角。
⑤目标函数:
在考虑土体的容重超载时,下限法目标函数为:
Maximize:kγ(tw,tm) (9)
式中:Maximize表示求“最大”;km(tw,tm)是第tw个地下水位作用下与第tm个抗剪参数的随机数相关的强度储备系数的随机数,tw=(1,…,Wn),Wn是土质边坡地水位随机数的数量;tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量。
⑥建立求解土质边坡可靠度的下限法线性规划模型:
步骤八、求解土质边坡下限法线性规划模型。
具体步骤为:
Step1、将Hw(tw)从tw=1到tw=Wn循环带入渗流场计算公式,逐次计算土质边坡的Wn个渗流场,获得土质边坡中每个有限单元结点的孔隙水压力值:其中:tw=(1,…,Wn),e=(1,…,Ne);
Step 3、在tw=1到tw=Wn的每次循环中将cr(tm)、从tm=1到tm=N循环,将N个土体材料的凝聚力和摩擦角的随机数带入求解土质边坡可靠的下限法线性规划模型,使用“单纯形”优化算法求解下限法线性规划问题,求解得到Wn×N个容重超载系数的随机数[kγ(tw,tm),(tw=1,…,Wn),tm=(1,…,N)]和对应的有限单元的正应力和切应力;同时,在tw=1到tw=Wn的每次循环中使用“二分法”迭代求解土质边坡的Wn×N个强度储备系数[km(tw,tm),(tw=1,…,Wn),tm=(1,…,N)]和与每个强度储备系数对应的有限单元的正应力和切应力。
计算结果是:获得了土质边坡50×1000=50000个强度储备系数和对应有限单元的正应力和切应力。
实施例1土质边坡tw=1时土质边坡强度储备系数的直方图如图11所示;实施例1土质边坡tw=30时土质边坡强度储备系数的直方图如图12所示;实施例1土质边坡tw=50时土质边坡强度储备系数的直方图如图13所示;实施例1土质边坡在50个水头作用下土质边坡50×1000个强度储备系数的直方图如图14所示;实施例1土质边坡tw=1~50时土质边坡强度储备系数概率密度分布图如图15所示。实施例1土质边坡tw=5、tm=293时土质边坡的第一主应力等值线如图16所示;实施例1土质边坡tw=5、tm=293时土质边坡的第三主应力等值线如图17所示。
步骤九、统计土质边坡的可靠度指标。
具体计算方法如下:
①、计算第tw个随机地下水位作用下土质边坡N个强度储备系数的均值:
②、计算第tw个随机地下水位作用下土质边坡N个强度储备系数的标准差:
③、计算土质边坡Wn×N个强度储备系数的均值:
式中:Wn是土质边坡地水位随机数的数量,N是土质边坡土体材料凝聚力和摩擦角随机数的数量,μkm是土质边坡Wn×N个强度储备系数的平均值。
④、计算土质边坡在所有随机地下水位作用下的Wn×N个强度储备系数的标准差:
式中:Wn是土质边坡地水位随机数的数量,N是土质边坡土体材料凝聚力和摩擦角随机数的数量,σkm是Wn×N个强度储备系数的标准差。
⑤、计算土质边坡的失效概率。当土质边坡的强度储备系数土质边坡的失效功能函数如下:
式中:tw=(1,…,Wn),tm=(1,…,N),Wn是土质边坡地水位随机数的数量,N是土质边坡土体材料凝聚力和摩擦角随机数的数量;I(tw,tm)是土质边坡的失效功能函数,km(tw,tm)是第tw个地下水位作用下与第tm个抗剪参数的随机数相关的强度储备系数的随机数,“if”是条件判断语句。
土质边坡的失效概率按下式计算:
式中:Wn是土质边坡地水位随机数的数量,N是土质边坡土体材料凝聚力和摩擦角随机数的数量,Pf是土质边坡的失效概率。
⑥、根据计算结果绘制土质边坡强度储备系数的直方图、概率密度曲线和累积概率密度曲线,以及绘制强度储备系数的均值和标准差、边坡失效概率随地下水位变化的关系图,绘制边坡的应力等值线。
计算结果是:实施例1土质边坡的50个地下水位作用下强度储备系数的均值标准差如表1所示;实施例1土质边坡的Wn×N个强度储备系数的均值μkm标准差σkm如表2所示;均值随地下水位变化如图18所示,标准差随水库上游库水位随机水头变化如图19所示;土质边坡的失效概率如表2所示,失效概率随地下水位的变化关系如图20所示。
土质边坡Wn×N个强度储备系数的均值是1.867、标准差是0.480,土质边坡的失效概率是0.928%。
表1第tw个地下水位作用下土质边坡N个强度储备系数的均值、标准差统计表
表2土质边坡Wn×N个强度储备系数的均值、标准差和失效概率统计表
项目 | 计算结果 |
均值μ<sub>km</sub> | 1.867 |
标准差σ<sub>km</sub> | 0.480 |
失效概率P<sub>f</sub>(%) | 0.928 |
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (9)
1.一种随机渗流场作用下的土质边坡可靠度分析下限法,其特征在于:所述方法具体步骤如下:
步骤1、拟定土质边坡的计算参数;
步骤2、建立土质边坡可靠度计算的极限状态函数;
步骤3、生成土质边坡地下水水位的随机数;
步骤4、生成土质边坡土体材料凝聚力、摩擦角的随机数;
步骤5、使用有限单元离散土质边坡,得到有限单元网格;
步骤6、进行土质边坡的随机渗流场计算;
步骤7、建立土质边坡可靠度分析下限法线性规划模型;
步骤8、求解土质边坡可靠度分析下限法线性规划模型;
步骤9、统计土质边坡的可靠度指标;
所述步骤7以有限单元结点的应力为决策变量建立求解土质边坡可靠的的下限法线性规划模型,具体为:①建立土质边坡有限单元的平衡方程;②建立土质边坡有限单元公共边的应力连续条件;③建立土质边坡有限单元的屈服条件;④建立土质边坡有限单元的静力边界条件;⑤建立目标函数;⑥根据土质边坡的极限状态函数、①、②、③、④和⑤建立的函数建立求解土质边坡可靠的下限法线性规划模型;
具体为:
①建立土质边坡有限单元的平衡方程为:
式中:e=(1,…,Ne),Ne是土质边坡中有限单元的数量;tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的数量;tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量;kγ(tw,tm)是第tw个地下水水位作用下与第tm个土体材料凝聚力、摩擦角的随机数相关的容重超载系数的随机数; 是分别是线性三角形有限单元e的6个形函数系数;分别是有限单元e的第1、2、3个结点沿x方向的有效正应力;分别是有限单元e的第1、2、3个结点沿y方向的有效正应力;分别是有限单元e的第1、2、3个结点的切应力;分别是第tw个地下水水位作用下有限单元e的第1、2、3个结点的孔隙水压力;Ae是有限单元e的面积;γe是有限单元e的容重;
②建立土质边坡有限单元公共边的应力连续条件:
式中:d=(1,…,Nd),Nd是土质边坡中有限单元公共边的数量; θd是公共边的倾角,逆时针为正;分别是公共边d的第1、2、3、4个结点沿x方向的有效正应力;分别是公共边d的第1、2、3、4个结点沿y方向的有效正应力;分别是公共边d的第1、2、3、4个结点的切应力;
③建立土质边坡有限单元的屈服条件:
式中:n=(1,…,Nn),Nn是土质边坡中有限单元结点的总数量;σn=[σ'xn,σ'yn,τxyn]T; Ck=2sin(2kπ/Mp);k=(1,…,Mp),Mp土体材料摩尔库伦屈服条件线性化正多边形的边数;σ'xn是第n个结点沿x方向的有效正应力;σ'yn是第n个结点沿y方向的有效正应力;τxyn是第n个结点的切应力;是土质边坡土体材料摩擦角的第tm个随机数;cr(tm)是土质边坡土体材料凝聚力的第tm个随机数;
④建立土质边坡有限单元的静力边界条件:
式中:b=(1,…,Nb),Nb是土质边坡中边界上的有限单元的数量; 分别是边界单元b中靠边界上的第1、2个结点的沿x方向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点的沿x方向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点的切应力;分别是边界单元b中靠边界上的第1、2个结点已知的沿法向的有效正应力;分别是边界单元b中靠边界上的第1、2个结点已知的沿切向的切应力;分别是第tw个地下水水水位作用下边界单元b中靠边界上的第1、2个结点的孔隙水压力;θb是边界单元b靠边界上的第1、2个结点连线的倾角;
⑤目标函数:
在考虑土体的容重超载时,下限法目标函数为:
Maximize:kγ(tw,tm)
⑥建立求解土质边坡可靠度的下限法线性规划模型:
式中:Z是土质边坡的极限状态函数。
2.根据权利要求1所述的随机渗流场作用下的土质边坡可靠度分析下限法,其特征在于:所述土质边坡的计算参数包括土质边坡的几何参数,土质边坡的地层条件,土质边坡土体材料的凝聚力、摩擦角、容重和渗透系数,土质边坡地下水水位的信息。
5.根据权利要求1所述的随机渗流场作用下的土质边坡可靠度分析下限法,其特征在于:所述步骤4包括:
①假设土质边坡土体材料的凝聚力和摩擦角符合对数正态分布;
②确定土质边坡材料的凝聚力、摩擦角的均值、标准差;
③土质边坡土体材料凝聚力的随机量和摩擦角的随机量的随机数按下式生成:
7.根据权利要求1所述的随机渗流场作用下的土质边坡可靠度分析下限法,其特征在于:所述步骤6具体为:包括:
①基于步骤5得到的有限单元网格,根据稳定渗流理论进行土质边坡的随机渗流场计算,二维稳定渗流公式如下:
式中:kx是土体材料x方向的渗透系数,ky是土体材料y方向的渗透系数,Hr是土质边坡内各点的随机水头函数;
8.根据权利要求1所述的随机渗流场作用下的土质边坡可靠度分析下限法,其特征在于:所述步骤8具体为:
S1、将土质边坡地下水水位随机数Hw(tw)从tw=1到tw=Wn循环带入渗流场计算公式,逐次计算土质边坡的Wn个渗流场,获得土质边坡中每个有限单元结点的孔隙水压力值:其中:tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的数量;e=(1,…,Ne),Ne是土质边坡中有限单元的数量;
9.根据权利要求1所述的随机渗流场作用下的土质边坡可靠度分析下限法,其特征在于:所述步骤9具体为:
①、计算第tw个随机地下水水位随机数Hw(tw)作用下土质边坡N个强度储备系数的均值:
式中:tw=(1,…,Wn),Wn是土质边坡地下水水位随机数的数量;tm=(1,…,N),N是土质边坡土体材料凝聚力/摩擦角随机数的数量;是第tw个随机地下水水位随机数Hw(tw)作用下土质边坡的N个强度储备系数的平均值;
②、计算第tw个随机地下水水位随机数Hw(tw)作用下土质边坡N个强度储备系数的标准差:
③、计算土质边坡Wn×N个强度储备系数的均值:
式中:μkm是土质边坡Wn×N个强度储备系数的平均值;
④、计算土质边坡在所有随机地下水水位作用下的Wn×N个强度储备系数的标准差:
式中:σkm是Wn×N个强度储备系数的标准差;
⑤、计算土质边坡的失效概率
土质边坡的失效功能函数如下:
土质边坡的失效概率按下式计算:
式中:Pf是土质边坡的失效概率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811568164.9A CN109753701B (zh) | 2018-12-21 | 2018-12-21 | 一种随机渗流场作用下的土质边坡可靠度分析下限法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811568164.9A CN109753701B (zh) | 2018-12-21 | 2018-12-21 | 一种随机渗流场作用下的土质边坡可靠度分析下限法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109753701A CN109753701A (zh) | 2019-05-14 |
CN109753701B true CN109753701B (zh) | 2022-05-24 |
Family
ID=66402942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811568164.9A Active CN109753701B (zh) | 2018-12-21 | 2018-12-21 | 一种随机渗流场作用下的土质边坡可靠度分析下限法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109753701B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113554259B (zh) * | 2020-04-26 | 2024-04-16 | 中国石油化工股份有限公司 | 一种渗流阻力与驱替压力双场匹配优化评价方法 |
CN112115530B (zh) * | 2020-08-05 | 2022-12-30 | 昆明理工大学 | 一种土质边坡的点失效概率的计算方法 |
CN112733392B (zh) * | 2020-10-30 | 2022-02-15 | 成都理工大学 | 基于二分类的边坡系统可靠度分析方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105224704A (zh) * | 2014-06-25 | 2016-01-06 | 中国科学院海洋研究所 | 一种海底不稳定性的评价方法 |
CN105279361A (zh) * | 2015-04-10 | 2016-01-27 | 中国电建集团成都勘测设计研究院有限公司 | 高土石不过水围堰边坡失稳风险率检测方法 |
CN106126796A (zh) * | 2016-06-20 | 2016-11-16 | 昆明理工大学 | 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法 |
CN107067333A (zh) * | 2017-01-16 | 2017-08-18 | 长沙矿山研究院有限责任公司 | 一种高寒高海拔高陡边坡稳定性监控方法 |
CN107330146A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑平动和转动效应的节理岩质边坡极限承载力分析上限法 |
CN108108556A (zh) * | 2017-12-20 | 2018-06-01 | 中国水利水电科学研究院 | 一种基于耗散-汇合结构的灌区水循环模型构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9506962B2 (en) * | 2013-05-08 | 2016-11-29 | Litecontrol Corporation | Resource control system with local reallocation |
-
2018
- 2018-12-21 CN CN201811568164.9A patent/CN109753701B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105224704A (zh) * | 2014-06-25 | 2016-01-06 | 中国科学院海洋研究所 | 一种海底不稳定性的评价方法 |
CN105279361A (zh) * | 2015-04-10 | 2016-01-27 | 中国电建集团成都勘测设计研究院有限公司 | 高土石不过水围堰边坡失稳风险率检测方法 |
CN106126796A (zh) * | 2016-06-20 | 2016-11-16 | 昆明理工大学 | 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法 |
CN107067333A (zh) * | 2017-01-16 | 2017-08-18 | 长沙矿山研究院有限责任公司 | 一种高寒高海拔高陡边坡稳定性监控方法 |
CN107330146A (zh) * | 2017-05-26 | 2017-11-07 | 昆明理工大学 | 一种同时考虑平动和转动效应的节理岩质边坡极限承载力分析上限法 |
CN108108556A (zh) * | 2017-12-20 | 2018-06-01 | 中国水利水电科学研究院 | 一种基于耗散-汇合结构的灌区水循环模型构建方法 |
Non-Patent Citations (2)
Title |
---|
"Improving the reliability of the frequency response function through semi-direct finite element model updating";Jae Hyuk Lim;《Aerospace Science and Technology》;20160731;第54卷;第59-71页 * |
"成都地区非饱和膨胀土边坡稳定性的有限元分析";吴宜峰;《中国优秀硕士学位论文全文数据库 基础科学辑》;20150930(第09期);第A011-11页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109753701A (zh) | 2019-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109753701B (zh) | 一种随机渗流场作用下的土质边坡可靠度分析下限法 | |
Karra et al. | Modeling flow and transport in fracture networks using graphs | |
CN108763697B (zh) | 一种岩质边坡中岩块失效概率的计算方法 | |
Millward-Hopkins et al. | Estimating aerodynamic parameters of urban-like surfaces with heterogeneous building heights | |
CN104933245A (zh) | 基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 | |
CN108170898B (zh) | 一种节理岩质边坡可靠度分析下限法 | |
CN112562078A (zh) | 一种三维地质分析预测模型构建方法 | |
CN112182731A (zh) | 非极限状态二维边坡稳定性评价方法 | |
CN107992687A (zh) | 基于蒙特卡洛方法的岩质边坡可靠度塑性极限分析上限法 | |
CN114297864B (zh) | 一种受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法 | |
Feldman et al. | Modeling the accumulation of wind-driven snow | |
Bahmanikashkooli et al. | Application of particle swarm optimization algorithm for computing critical depth of horseshoe cross section tunnel | |
CN110850477A (zh) | 一种基于松鼠搜索算法的电阻率溶洞识别方法 | |
CN109753700B (zh) | 一种土石坝可靠度分析上限法 | |
Chiozzi et al. | Fast kinematic limit analysis of FRP reinforced masonry vaults through a new genetic algorithm nurbs-based approach | |
Yazd et al. | Optimal designing of concrete gravity dam using particle swarm optimization algorithm (PSO) | |
CN109583135A (zh) | 针对s型沟槽填方场地边坡稳定性的极限平衡分析方法 | |
CN103268423A (zh) | 一种地理现象多点模拟空间尺度选择的方法 | |
CN115330981B (zh) | 一种边坡滑面搜索方法、系统、设备及可读存储介质 | |
CN113901552B (zh) | 一种堰体模型的建立方法及装置 | |
CN113435087B (zh) | 一种分析洞穴围岩局部稳定性的方法 | |
CN101221598B (zh) | 小流域淤地坝坝系规划用坝系拓扑视图的生成方法 | |
CN102968534A (zh) | 基于稳定分布的楼面活荷载标准值设计方法 | |
CN108446413A (zh) | 一种注浆成型扩底桩桩径的优化测定方法 | |
CN108319761A (zh) | 确定群桩位移的方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |