基于非侵入式随机有限元法的地下洞室可靠度分析

提出了地下洞室变形可靠度分析的非侵入式随机有限元法。介绍了随机多项式展开的基本原理,采用 GEOSLOPE 的 SIGMA/W 模块进行确定性有限元分析。提出了随机多项式展开与 SIGMA/W 模块接口方法及其流程图,从而实现了确定性有限元分析和随机分析一体化。最后研究了非侵入式随机有限元法在地下洞室变形可靠度分析中的应用。结果表明：非侵入式随机有限元法使得随机分析与确定性有限元分析互不耦合,其计算效率是传统的蒙特卡罗模拟方法无可比拟的,它是地下洞室变形可靠度问题分析一种有效的方法。采用二次衬砌支护是提高地下洞室可靠度有效的方法。此外,岩体变形模量的变异性对地下洞室变形可靠度有非常明显的影响,而岩体重度的变异性对可靠度基本上没有影响。因此,在地质勘查中要尽可能准确地确定岩体的变形模量,从而有效地提高地下洞室变形可靠度。     

考虑空间变异特性的边坡可靠度分析简化方法

介绍了一种衡量边坡稳定性的简化处理方法,采用蒙特卡洛法与极限平衡法,对某均质边坡进行了可靠度分析,并将简化方法的结果与已有结果进行了对比,结果表明:失效概率误差控制在20%以内,可以尝试采用简化方法对边坡进行可靠度分析。     

考虑地层变异性和土体参数空间变异性的边坡可靠度全概率设计方法

目前边坡稳定设计研究中大多数考虑了土体参数的空间变异性,但忽略了地层变异性的影响。为此,提出了一种同时考虑这两类变异性的边坡可靠度全概率设计方法。在全概率设计框架内,将广义耦合马尔可夫链模型与随机场模型相耦合用于同时表征地层变异性和土体参数空间变异性,给出了所提方法的计算流程。利用澳大利亚珀斯市钻孔资料,以某边坡为例进行可靠度设计,为说明在边坡可靠度设计中同时考虑地层变异性和土体参数空间变异性的重要性,分析了仅考虑土体参数空间变异的情况,进一步分析了同时考虑两类变异性的情况,并对二者进行了比较。结果表明:所提出的边坡可靠度设计方法能够有效地考虑边坡中存在的地层变异性和土体参数空间变异性。当仅考虑土体参数空间变异性时,边坡可靠度设计结果很大程度上取决于所采用地层的分布情况,特别是地层分布中抗剪强度较强土体类型占比高于真实情况时,将导致得到的最优设计方案偏于危险。反之,若地层分布中抗剪强度较弱的土体类型占比高于真实情况,得到的最优设计方案偏于保守。因此,为准确地得到最优设计方案,在边坡可靠度设计中应同时考虑地层变异性和土体参数空间变异性的影响。     

考虑时效特性的锚固岩质边坡变形可靠度分析

 提出基于非侵入式随机有限元法的岩质边坡变形可靠度分析方法。建立考虑腐蚀效应的锚杆锚固力随服役时间变化模型。给出边坡变形可靠度分析非侵入式随机有限元法的计算步骤。研究锚固边坡变形失效概率与最大允许变形值之间的关系,并基于参数敏感性分析提出边坡最大允许变形取值方法。以锚固岩质边坡为例,证明所提方法的有效性。结果表明：非侵入式随机有限元方法为考虑时效特性的边坡变形可靠度分析提供一条有效的途径。随着锚杆服役时间增加,锚杆腐蚀作用对锚固边坡变形可靠度的影响越明显。考虑时效特性的边坡变形失效概率常用对数 与边坡最大允许变形值间存在近似线性关系,而且随着边坡可靠度水平增加,该线性关系表现得越明显。     

考虑土体空间变异性的边坡最危险滑动面随机分析方法

现有边坡可靠度研究均未考虑土体空间变异性对边坡最危险滑面的影响。为此，提出了考虑土体空间变异性的边坡最危险滑动面随机分析方法。首先，采用谱表现法建立了表征土体空间变异性的随机场模型。在此基础上，提出了基于SIGMA/W和SLOPE/W的自动定位搜索最危险滑动面方法。其次，采用非侵入式随机分析方法研究了抗剪强度参数空间变异性对边坡最危险滑动面空间分布的影响。最后，采用算例验证了所提方法的有效性。结果表明：提出的边坡最危险滑动面随机分析方法能够有效地确定边坡最危险滑动面空间分布特征。土体抗剪强度参数的空间变异性对边坡最危险滑动面的空间分布特征有重要的影响，它直接决定了边坡最危险滑动面的位置和滑体规模。土体抗剪强度参数波动范围越大，最危险滑动面的空间分布范围越大。随着土体抗剪强度参数水平向和竖直向波动范围比值的增大，边坡上部发生局部滑动的可能性增大。抗剪强度参数的变异系数越大，最危险滑动面的空间分布范围越大，边坡发生小规模局部滑动的可能性越大。     

考虑参数空间变异性的无限长边坡可靠度分析

目前有关土体参数空间变异性对边坡稳定性影响的研究没有考虑抗剪强度参数随深度变化的影响。为此,提出了考虑土体抗剪强度参数均值随深度变化的无限长边坡稳定性概率分析方法。采用Karhunen-Loeve展开建立了表征土体空间变异性的随机场模型。探讨了考虑土体抗剪强度参数空间变异性时边坡失效概率和最危险滑动面的变化规律。最后,以无限长边坡稳定性概率分析问题为例验证了所提方法的有效性。结果表明：土体抗剪强度参数的空间变异性对无限长边坡失效概率有明显的影响,边坡失效概率随土体抗剪强度参数相关距离的增加而减小。对于不排水黏性边坡来说,边坡不排水抗剪强度随深度变化越明显,边坡失效概率越小。边坡最危险滑动面大部分都位于边坡底部。对于摩擦/黏性边坡来说,随着边坡抗剪强度随深度变化的增强,边坡失效概率有先减后增的趋势。抗剪强度参数随深度变化不同趋势对最危险滑动面分布规律有明显的影响。不考虑抗剪强度参数均值随深度变化将会导致最危险滑动面最可能出现在边坡底部。     

考虑地层变异性和土体参数变异性的边坡可靠度分析

现有边坡可靠度分析大多数只考虑土体参数的固有变异性,而忽略了地层变异性的影响。为此,提出了同时考虑地层变异性和土体参数固有变异性的边坡可靠度分析方法,利用耦合马尔可夫链模拟地层分布,采用基于乔列斯基分解的中点法离散土体参数随机场,采用有限元强度折减法计算边坡安全系数,通过蒙特卡洛法模拟进行边坡可靠度分析。利用澳大利亚珀斯市钻孔资料,以边坡可靠度问题为例阐明了同时考虑地层变异性和土体参数固有变异性的重要性,研究了钻孔布置方案对边坡可靠度的影响规律,结果表明:提出的方法能够有效地反映地层变异性和土体参数固有变异性对边坡可靠度的影响。当钻孔数目较少时,模拟的边坡土体类型分布与真实边坡土体类型分布相差较大,此时忽略地层变异性将导致边坡可靠度不精确的估计结果。钻孔布置方案对边坡失效概率和安全系数有明显的影响,钻孔应尽可能多的布置在边坡关键影响区域。边坡失效概率和安全系数统计量与钻孔数目并不呈单调关系,但是随着钻孔数目的增加,边坡失效概率和安全系数统计量逐渐收敛至"精确解"。     

考虑岩体空间变异性的隧洞围岩变形随机分析

准确和全面地预测隧洞开挖所引起的围岩变形是隧洞设计和施工中的一个关键问题。目前常用的确定性分析方法不能考虑岩体的空间变异性,而这种空间变异性对围岩变形计算结果的准确性和可靠性有着重要的影响。为此,提出采用非侵入式随机分析方法研究岩体空间变异性对围岩变形的影响。采用随机场模拟岩体的空间变异性,采用有限元计算围岩的变形,研究岩体弹性模量和侧压力系数的空间变异性对围岩变形特别是地表位移的影响。结果表明,非侵入式随机分析方法的优势在于有限元计算和随机分析不耦合。侧压力系数的空间变异性对围岩变形的影响远小于弹性模量的影响。由于低强度占优效应的影响,岩体的特征模量小于数学平均值。弹性模量随机场中横向和竖向波动范围参数对围岩变形有明显不同的影响。     

基于极限平衡法和有限元法的边坡协同式可靠度分析

极限平衡法和有限元法是两种常用的边坡稳定分析方法。基于极限平衡法的边坡可靠度分析计算效率较高,但需要假定失效模式,从而导致计算结果不准;与之相反,基于有限元法的边坡可靠度分析更为严格,但计算效率较低。为此,提出了一种新的基于随机模拟的边坡可靠度分析方法——边坡协同式可靠度分析方法。该方法包括初步可靠度分析和精细可靠度分析两步,可以同时利用极限平衡法和有限元法的优势,实现既高效又准确的边坡可靠度分析。通过一个考虑空间变异性的两层土坡算例验证了该方法的有效性,结果表明:协同式可靠度分析方法与基于有限元法的蒙特卡洛模拟或子集模拟相比,不仅具有一致的可靠度分析结果,而且显著提高了小概率水平下的计算效率,促进了基于有限元法的边坡可靠度分析在实际工程中的应用。该方法可以将基于极限平衡法的边坡可靠度分析成果合理纳入到基于有限元法的边坡可靠度分析中,从而获得大量的失效样本,以制定合理的边坡防治措施。该方法非常适用于高维可靠度问题,如考虑空间变异性的边坡可靠度问题。     

锦屏一级水电站左岸坝肩边坡施工期高效三维可靠度分析

 基于可靠度理论的边坡稳定性分析是一个重要的发展趋势,然而目前国内外关于复杂三维高边坡可靠度方面的研究几乎还是空白。以锦屏一级水电站左岸坝肩边坡为例,研究了施工阶段三维边坡可靠度问题。由于三维边坡稳定性分析有限差分强度折减法计算量较大,并且安全系数没有显式表达式,通过参数敏感性分析减少随机变量数目,提出三维边坡可靠度分析的非侵入式随机有限差分法。探讨了预应力锚索和深层混凝土抗剪洞2种主要加固措施失效对边坡变形、稳定及其可靠度的影响规律。研究了施工阶段边坡横河向位移、安全系数以及失效概率的变化规律。结果表明：参数敏感性分析可以有效地识别敏感性较大的随机变量,从而提高边坡可靠度计算效率。非侵入式随机有限差分法实现了概率分析与通用商业软件如FLAC3D的有机结合,极大地简化了可靠度分析过程,而且能充分利用参数敏感性分析中边坡稳定性计算结果,为解决三维高边坡可靠度问题提供一种有效的工具。此外,锦屏一级水电站左岸坝肩边坡所采取的预应力锚索和深层混凝土抗剪洞联合加固措施能够有效地控制边坡变形和确保边坡稳定性。施工阶段边坡开挖扰动作用对边坡稳定可靠度有着重要的影响。     

考虑参数空间变异性的边坡稳定可靠性有限元极限分析

当土性参数空间变异性较大时,极限平衡法得到的滑移面不尽合理。阐述了基于广义变分原理的有限元极限分析方法,采用混合有限元方法,构筑了线性应力三角形单元与线性速度三角形单元,结合强度折减法与线性规划算法,建立了边坡稳定安全系数上下限分析方法,分析了土的抗剪强度参数空间变异性对边坡稳定性的影响,并与3种典型极限平衡法进行了对比。结果表明,FELA方法可有效搜索边坡临界滑移面,并给出安全系数的严格上下限。对于简单均质边坡,有限元极限分析与极限平衡法结果接近,极限平衡法结果大多位于极限分析的上下限内;对于空间变异性较大的边坡,有限元极限分析法可以有效搜索可能的多种临界滑移面,而极限平衡法则存在显著偏差,且往往高估滑坡风险。强度参数的空间变异性还导致边坡安全系数分布形式变化显著,仅采用安全系数无法反应这一变化。根据安全系数的分布形式,给出了土性参数设计值建议。     

考虑堆石料空间变异性的土石坝地震反应随机有限元分析

探讨了一种将随机场模拟技术与有限元方法相结合的随机有限元法在土石坝地震反应分析中的应用途径。采用局部平均细分法模拟堆石料物理力学性质随机场,使得能够考虑堆石料的空间变异性。提出尽量利用建立材料物理性质与力学性质关系的计算模型,以降低获取材料力学性质统计特性信息的成本。以一实际工程为例,采用该随机有限元法计算了堆石料的干密度、孔隙比、不均匀系数、平均粒径、Duncan–Chang E–B模型模量系数以及初始摩擦角具有不确定性时土石坝的随机地震响应及永久变形。结果表明:筑坝堆石料的不确定性会使大坝响应及永久变形发生一定程度的离散;忽略材料的不确定性可能导致低估大坝的地震反应;表征大坝地震反应的各性能指标并不一定都符合正态或对数正态分布。     

考虑概率分布影响的低概率水平边坡可靠度分析

目前考虑土体参数空间变异性的边坡可靠度分析所研究的边坡几何尺寸相对较小。为有效地分析考虑参数空间变异性的几何尺寸相对较大的边坡可靠度问题,提出了基于多重响应面与子集模拟的边坡可靠度分析方法。以一个坡高为24 m的两层非均质黏土边坡为例验证了提出方法的有效性,并探讨了正态、对数正态、极值I型、Gamma和Beta这5种概率分布类型对边坡可靠度的影响。结果表明,提出方法具有以下优势:1可以有效地计算考虑多参数空间变异性的低概率水平边坡可靠度;2可以较好地解决几何尺寸相对较大的边坡可靠度问题;3有较高的参数敏感性分析计算效率,可为调查概率分布类型和波动范围等参数统计特征对边坡可靠度的影响提供技术支持。此外,参数概率分布类型对边坡可靠度具有重要的影响,常用的正态和对数正态分布分别用于表征参数概率分布特征时,可能会分别高估和低估边坡失效概率。     

边坡的弹塑性有限元可靠度分析

以弹塑性有限元理论和可靠度理论为基础,基于偏微分技术及增量初应力法,对采用Mohr-Coulomb屈服准则的边坡土体进行了弹塑性随机有限元可靠度分析,比较了基于强度折减法和基于滑面应力分析的弹塑性随机有限元分析法的异同:前者编程简单,可利用现有大型有限元计算软件计算边坡的整体可靠指标,但不能确定边坡中相应的潜在滑动面的位置,计算量较大,运算速度较慢;后者可以同时求出边坡的整体可靠指标及潜在滑动面的位置,运算速度快,但编程复杂。文中还研究了基于滑面应力分析的弹塑性随机有限元可靠度分析方法中滑面可靠指标的4种计算方法,比较了不同功能函数形式对边坡可靠指标及滑面位置的影响,指出采用考虑滑面方向的功能函数形式来求解可靠指标的方法3与基于强度折减法的可靠度分析本质相同,较为合理。