边坡稳定可靠性的随机有限元分析

随机有限元法可以处理土性参数的变异性和空间相关性。对二阶摄动随机有限元法的摄动理论和程序进行了研究，提出了偏导矩阵的组集方法，采用正态分布随机变量的正交变换法来提高计算效率。考虑土性参数随机场作用和土性指标之间的互相关性，建立边坡局部抗剪和总体稳定性可靠度的随机有限元分析模型，对某实际土坡进行了可靠度计算，计算结果较为符合实际。     

顺层边坡稳定性及可靠度的随机有限元分析法

阐述了基于摄动法的随机有限元基本原理,采用以均值为中心的2阶摄动随机有限元,运用局部平均的空间离散方法,对顺层边坡的应力场进行了2阶统计均值、协方差和方差的分析,并结合顺层边坡沿软弱夹层滑动的破坏特征,推导了基于Mohr-Coulomb准则的顺层边坡整体稳定的安全系数和可靠度的计算表达式。结合典型顺层岩质路堑边坡工程实例,采用摄动随机有限元法进行了边坡稳定性和可靠度的分析计算。结果表明,文中介绍的摄动随机有限元法可以考虑多种随机变量对边坡系统的影响,可同时采用安全系数和可靠度作为边坡稳定性的评价指标,判定结果更为准确和全面。     

基于随机响应面法的边坡系统可靠度分析

提出了一套基于随机响应面法的边坡系统可靠度分析方法。该方法首先从大量潜在滑动面中筛选出代表性滑动面。针对每条代表性滑动面,采用Hermite多项式展开建立其安全系数与土体参数间的非线性显式函数关系（即随机响应面）。然后,采用直接蒙特卡洛模拟计算边坡系统失效概率。在蒙特卡罗模拟中,采用所有代表性滑动面的随机响应面计算每一组样本所对应的边坡最小安全系数。最后,以两个典型多层边坡系统可靠度问题为例验证了该方法的有效性。结果表明：文中提出的边坡系统可靠度分析方法能够有效地识别边坡代表性滑动面,具有较高的计算精度和效率,并且确定代表性滑动面时无需计算滑动面间的相关系数。同时该方法可以有效地计算低失效概率水平的边坡系统可靠度,为含相关非正态参数的边坡系统可靠度问题提供了一条有效的分析途径。此外,多层边坡可能同时存在多条潜在滑动面,基于单一滑动面（如临界确定性滑动面）或者部分代表性滑动面进行边坡系统可靠度分析均会低估边坡失效概率。     

边坡稳定可靠度分析的并行随机有限元方法

将Monte-Carlo随机方法与并行有限元方法结合,设计了基于PC机群的并行随机有限元算法,每一次随机模拟产生一个随机场,作为一个样本,代入有限元进行稳定性计算,发挥并行有限元计算量大、计算速度快的特点,克服了Monte-Carlo方法计算时间长的缺点。用C++语言结合MPI消息传递接口开发了相应的计算软件,计算了1 000个随机样本的边坡可靠度问题,验证了程序的可靠性及高效性。     

约束随机场下的边坡可靠度随机有限元分析方法

目前边坡可靠度中常用的简化分析方法,不考虑边坡土体的空间变异性,每次计算整个边坡都取用相同的强度参数,由离散点试样试验得到的土体参数统计特性只能反映点特性,而边坡的稳定性受滑面上平均抗剪强度特性控制,因此,需要考虑空间范围内的平均特性。描述空间变异性的随机场理论对变异性较高的土体,实际上高估了其空间变异性。把随机场理论和地质统计中的区域化变量理论结合起来,建立约束随机场,并在此基础上进行Monte-Carlo随机有限元分析。计算实例表明,在高变异性条件下约束随机场能有效降低完全随机场的模拟方差,得到更低的破坏概率。对比了随机有限元和简化法的计算结果表明,简化法在土体强度变异性很高时其结果并非偏于保守。另外也指出了可靠度分析中存在的边坡尺度效应和简化法的适用条件。     

隧道围岩弹塑性随机有限元分析及可靠度计算

在弹塑性有限元分析的初应力法基础上,推导了弹塑性随机有限元增量初应力法的计算迭代格式,结合可靠度理论,给出了围岩稳定的可靠度分析方法,在计算出可靠度的同时还给出了围岩和锚喷支护结构的应力特征值。并对一工程算例进行了详细的分析和讨论,得出了一些有益的结论。     

随机地震作用下滑坡稳定性可靠度分析

应用边坡可靠性分析中常用的蒙特卡罗法,并结合稳定性计算的剩余推力法和Sarma法,研究随机地震作用下边坡的可靠性。以三峡库区黄土坡滑坡为例,分析了地震效应对边坡稳定性的影响。计算结果表明：地震加速度对滑坡稳定性影响显著;竖向地震加速度小于水平地震加速度2/3时,对滑坡稳定性影响较小,大于水平加速度时,对滑坡稳定性影响明显。因此在一般工程设计中,可以只考虑水平向的地震作用,而对于设计烈度较大的大型工程,则应同时考虑水平向和垂直向的地震作用。     

薛家湾开源路高路堑边坡工程稳定性的三维有限元分析

针对准格尔旗薛家湾镇开源路高边坡欠稳工程,采用MIDAS/GTS有限元软件对其施工过程及边坡稳定性进行了三维数值模拟,得到了不同工况下的变形及支护完成后的安全系数,并与实际工程现场所测变形结果作了比较分析。结果表明,三维FEM计算结果与工程实际变形较为相似,其结果对类似高路堑边坡的设计及加固有一定的参考价值。     

边坡可靠度分析的非侵入式随机有限元法

提出基于非侵入式随机有限元法的边坡可靠度分析方法,并编写计算程序NISFEM。采用有限元滑面应力法计算边坡安全系数,将Hermite随机多项式展开与SIGMA/W和SLOPE/W模块有机结合实现边坡可靠度非侵入式随机分析。根据随机多项式展开系数,给出边坡安全系数前4阶统计矩（均值、标准差、偏度和峰度）和Sobol指标解析表达式,并采用Sobol指标进行边坡可靠度参数敏感性分析。最后,以均质土坡可靠度问题为例,证明该方法在边坡可靠度分析中的有效性。结果表明,边坡可靠度分析的非侵入式随机有限元法能够有效地考虑边坡变形对边坡可靠度的影响,计算效率远远高于蒙特卡罗模拟方法（MCS）,是解决复杂边坡可靠度问题一种有效地分析手段;黏聚力和内摩擦角变异性对边坡安全系数前四阶统计矩具有明显的影响,重度变异性对安全系数前4阶统计矩几乎没有影响;抗剪强度参数间负相关性对边坡安全系数均值几乎没有影响,但对安全系数标准差、偏度和峰度均有明显的影响。此外,随着抗剪强度参数间负相关性的增加,边坡安全系数由近似正态分布逐渐变为明显的非正态分布。     

随机有限元洁分析水工隧洞结构的可靠度

首次提出用平面等参随机有限元分析水工隧洞结构可靠度的方法，推导了有关公式，编写了计算程序。结合实例分析了衬砌和不衬砌隧洞的可靠度，并探讨了渗流场对隧洞结构可靠度的影响。     

心墙堆石坝渗透稳定可靠性分析的随机响应面法

将大规模渗流有限元计算与随机响应面法相结合,对双江口心墙堆石坝进行渗透稳定可靠性分析。在基于随机响应面法的可靠度分析框架内,堆石坝稳定渗流有限元计算过程和可靠度分析过程分开独立进行,通过对心墙渗透坡降较大区域的节点建立统一的渗透稳定功能函数,采用渗流有限元分析方法和随机响应面法,计算出该区域每个节点处的渗透破坏失效概率,并将最大失效概率作为心墙的失效概率。最后,分析了心墙渗透系数、覆盖层渗透系数、上游水位与心墙具有最大失效概率节点处渗透坡降的相关关系,以及心墙渗透系数和上游水位的变异性对心墙渗透破坏失效概率的影响。计算结果表明,随机响应面法3阶Hermite展开就能够保证良好的计算精度,且计算耗时较小;双江口堆石坝心墙具有最大失效概率节点处的渗透坡降与上游水位密切相关,而与心墙本身的渗透系数呈弱负相关关系,与覆盖层渗透系数的相关性不显著;随着上游水位变异性的增大,心墙失效概率急剧增大,而这种效应对于心墙渗透系数并不明显。研究成果为随机响应面法在实际工程中的应用奠定了一定的基础。     

基于可靠度的黄土斜坡稳定性分析

以延安市宝塔区虎头峁黄土斜坡为研究对象,在野外详细调查和勘查的基础上,分析确定斜坡可靠度的随机变量为滑带土土体容重r、内聚力c和内摩擦角φ。分别用统计矩法、一次二阶矩法、几何法和蒙特卡洛模拟法计算斜坡的可靠度和破坏概率,得出安全系数为1.00~1.15,可靠度指标为0.35~0.45,破坏概率为32%~37%,说明斜坡具有较大的失稳可能性。计算过程表明,用几何法计算所得的可靠度指标精度最高,一次二阶矩法次之,统计矩法计算精度最低,蒙特卡洛模拟法可作为一种校验方法,对其他方法进行校验和对比。在不需要获取传统的定值安全系数的情况下,可用几何法计算斜坡的可靠度指标,方法较为简便。当需要同时获取安全系数时,应使用统计矩法或一次二阶矩法计算。     

地基固结沉降随机有限元计算和可靠度分析

基于比奥（Biot）固结理论，采用邓肯-张本构模型，考虑随机场的影响，首次给出了地基固结沉降随机有限元分析方法，并基于随机有限元对点固结沉降和差异固结沉降的可靠度进行了分析。     

基于有限元强度折减法的边坡动力稳定性可靠性分析

传统有限元强度折减法在边坡稳定性数值分析中取得了一定的成功,但由于未考虑岩土体材料参数的变异性等不确定性因素,尚不能直接应用于边坡稳定性特别是动力稳定性可靠性分析问题。为此提出了基于有限元强度折减法的地震边坡动力稳定性可靠性分析方法。将有限元极限分析法、动力分析法和可靠性分析法三者耦合,分析求解边坡在地震作用下的动力稳定性可靠性问题,并将这一过程在数值计算程序中得以实现。在计算分析过程中,克服了原方法需不断人工试算才能得到边坡安全系数而无法量化处理问题,并对边坡动力失效准则进行了适用于程序化的改进,使其计算过程完全实现自主运行。结合典型算例分析结果表明,该方法显著的特点是能较全面地反映岩土体的动力特性和边坡岩土体材料强度参数的变异性及相关性,所得结果相对更加合理且更符合工程实际。该方法既是对有限元强度折减法的应用范围的有益推广,也为边坡动力稳定性可靠性问题研究提供了一条新的有效途径。     

地震加速度动态分布及对高土石坝坝坡抗震稳定的影响

《水工建筑物抗震设计规范》（DL5073-1997）建议的土石坝地震加速度动态分布系数适合于150 m以下的土石坝,而前许多土石坝的设计高度已远大于150 m。与低坝相比,高坝受地基刚性的约束减弱,坝体的自振周期延长。在坝体地震反应中,高阶自振周期与地震卓越周期耦合的机率增大,高阶振型的振动易被激发放大,从而导致坝体地震加速度沿坝高的分布与低坝有所不同。据此,采用有限元法研究了高土石坝的加速度分布,提出了高于150 m土石坝的地震加速度动态分布系数图示。在此基础上,利用堆石材料的非线性强度准则,对高度超过150 m的土石坝坝坡抗震稳定性作了进一步的分析,得出随着坝体地震加速度动态分布系数的降低,坝坡临界破坏的安全系数有所提高。