自适应网格加密求解洞室围岩塑性区

在弹塑性交界面上存在应力跳跃。 将应力跳跃值作为加密网格的准则, 不仅可使应力计算精度大幅度提高, 而且可求出塑性区的形状。 针对塑性区求解这一具体问题, 提出了一种修正的网格加密准则。 通过改变影响围岩塑性区的几个主要影响因素, 进行了大量的数值计算, 结果表明所提出的修正加密准则能正确地求出围岩塑性区的形状和大小。     

一种基于有限元应力场的边坡滑动面搜索方法

如何将有限单元法计算得到的结果是应力场、应变场和位移场等计算成果应用于边坡稳定安全系数及滑裂面位置的确定,是有限单元法边坡稳定分析中的一个重要研究课题.提出一种基于有限元应力场的最危险滑裂面的搜索方法,编制了搜索程序,采用澳大利亚边坡稳定分析程序的考核题进行对比验算,验证了该方法的正确性及程序的可行性.     

三维边坡最不利滑裂面的遗传算法搜索

将遗传算法引入到三维边坡稳定分析中。由于三维边坡滑裂面的复杂性,需对滑裂面作适当简化。认为滑裂面为椭球面,用7个控制参数来模拟生成滑裂面,再运用遗传算法来搜索滑坡的最不利的滑裂面。对于存在确定滑裂面的滑坡,7以将此滑裂面作为整体滑裂面,计算其整体稳定安全系数,然后在整体稳定的基础上,运用遗传算法搜索滑坡内部的最不利滑裂面,得出滑坡最小的稳定安全系数。     

自适应差分进化算法在边坡滑面搜索中的应用

边坡三维临界滑裂面的确定及其安全系数的求解是边坡稳定分析中的关键问题,而现有分析方法尚难以方便、可靠地获得三维滑裂面。在有限元数值应力场的基础上,计算边坡的安全系数。通过构造扩展椭球形滑裂面,利用自适应差分进化(SaDE)算法结合反向移动操作,实现了复杂边坡三维临界滑裂面的搜索。扩展椭球滑裂面构造方便,搜索灵活,比传统的球面或椭球面更能反映实际工程情况。整个求解过程为显式求解,不存在收敛性问题。算例考核结果表明,提出的方法具有良好的搜索效率与全局寻优能力,结果与刚体极限平衡法、有限元强度折减法吻合,计算结果合理可靠。通过实际工程边坡,进一步验证了该方法可方便、有效地应用于具有复杂岩性及环境条件的边坡三维稳定性分析。     

边坡稳定极限分析的单元集成法

提出的单元集成法为采用塑性力学上限定理分析边坡的稳定性提供了一种通用性较强的手段。它采用类似于有限元网格划分的方式离散边坡计算区域,并设定一个机动许可的滑动机构（包括滑裂面位置和速度场,它与网格划分方式无关）,在此滑动机构下可以很方便地计算每个单元贡献的外力功率和内能耗散率。把所有单元的能量相加就是滑坡体的总能量。然后,根据上限定理求得安全系数,并通过非线性数学规划方法找到最小值。对于直线滑裂面或对数螺旋滑裂面的单滑块机制,几个典型算例说明了该方法的有效性。     

基于有限元计算的边坡三维滑裂面搜索

提出一种基于有限元计算的边坡三维滑裂面搜索的方法,根据有限元计算的应力结果,计算可能滑裂面的安全系数,然后用遗传算法搜索出最危险滑裂面,并得出相应的安全系数作为边坡的安全度评价指标。该方法的优越性是只需进行一次有限元计算,不需要假定滑裂面,可以用于具有复杂地貌和地质构造的三维边坡。通过算例,验证该方法的可行性和程序的正确性,以大岗山水电站库区边坡工程实例验证该方法的实用性。     

平面剪切岩质边坡滑裂面的确定及稳定性分析

由于传统搜索方法对岩质边坡滑裂面的确定无法兼顾效率与精度,如何迅速准确确定潜在滑裂面仍然是个难题.极限平衡法在岩质边坡稳定性分析中备受认可,采用岩质边坡平面剪切滑动模型,以滑裂面的倾角来表征潜在滑裂面的位置;基于极值法,推导了极限平衡条件下平面剪切破坏型岩质边坡潜在滑裂面的解析解,并结合香港秀茂坪路边坡对其准确性进行了...     

基于统一强度理论的土石坝边坡稳定分析遗传算法

将非线性统一强度理论应用于土石坝边坡稳定分析中,考虑坝体各应力分量包括中间主应力对土质材料强度的影响,以区别于传统土石坝稳定分析中仅考虑剪切和拉伸极限强度的极限平衡法。提出了一个土石坝边坡稳定优化分析模型,并采用最优化遗传算法对土石坝边坡进行断裂面优化搜索。研究表明,统一强度理论中反映中间主应力效应系数值变化对坝体稳定安全系数有显著影响。利用遗传算法对土石坝边坡稳定进行分析时,不必事先假定滑裂面的形状,根据土石坝非线性统一强度理论的剪裂破坏准则进行滑裂面搜索。针对常规最优化方法容易陷入局部最优解的缺陷,基于MATLAB遗传算法工具箱对土石坝边坡稳定性进行分析,分析结果表明,统一强度理论结合遗传算法能克服常规最优化方法在寻优过程中一些缺陷,计算结果更自然合理。     

遗传算法在土坡整体稳定性分析中的应用

边坡稳定性分析的关键是如何确定最危险滑动面的位置并计算与之相对应的安全系数。由于传统的极限平衡分析方法很容易陷入局部极小值而不能找到真正的最危险滑裂面,因此采用瑞典条分确立土坡分析模型,用遗传算法搜索土坡最危险滑动面,进而求得土坡最小安全系数。该方法模拟了生物遗传进化的过程,克服了传统方法的局限性。通过和面积细分法所搜索的最危险滑动面和计算得到的土质边坡安全系数作对比,可得遗传算法在土质边坡稳定分析中具有较高的精度与可靠性。遗传算法可以很好地解决如何寻找土质边坡整体极值的问题。工程应用实例表明,遗传算法分析土质边坡的稳定性效果良好,具有很好的应用前景。     

用强度折减方法分析重力坝深层抗滑稳定性

目前,对重力坝深层抗滑稳定的分析方法尚无统一和明确的规定。强度折减方法在求解安全系数方面具有不用指定滑面、自动搜索滑裂面等优点。本文基于FLAC软件,采用强度折减方法对分属三种典型地质条件的多个算例进行分析。分析结果表明,无论采用相关联还是非相关联流动法则,采用强度折减方法得到的安全系数与Morgenstern-Price方法求解的安全系数都相差较小,误差小于5%,且采用相关联流动法则时得到的安全系数误差更是小于3.5%。根据在破坏时坝基塑性区的分布,提出了确定滑裂面的方法,该方法确定的滑裂面与M-P方法搜索到的滑裂面有很好的一致性。本文通过多个算例的分析证明基于FLAC软件的强度折减方法是一种可靠、有效的方法,可以应用于分析重力坝深层抗滑稳定性。通过对向家坝水电站泄12坝段的深层抗滑稳定分析表明该方法可以应用到重力坝深层抗滑稳定的工程分析中。     

有限元强度折减法计算边坡稳定的对比分析

采用目前边坡稳定性分析比较流行的强度折减法,对比研究了MIDAS/GTS、FLAC、ANSYS配合Drucker-Prager（简称D-P）屈服准则和Mohr-Coulomb（简称M-C）屈服准则时软黏土、硬黏土、弱膨胀土3种工况下计算结果的偏差。软黏土工况下D-P准则和M-C准则计算结果的偏差相对较小,当边坡土体为硬黏土时,采用D-P准则与采用M-C准则计算结果的偏差明显增加。3种软件2种屈服准则下的计算结果都反映出,硬黏土的滑动面比弱膨胀土和软黏土的滑动面浅,而且同等情况下MIDAS计算得到的滑动面比ANSYS计算得到的滑动面浅;坡度较小时FLAC（M-C）计算的安全系数比MIDAS（M-C）计算得到的大,坡度较大时则相反;坡度较小时计算过程中先出现塑性区贯通,后出现计算不收敛;坡度较大时计算过程中先出现计算不收敛,后出现塑性区贯通。坡度较小时计算不收敛时的折减系数与出现塑性区贯通时的折减系数差别较大;坡度较大时这一差别较小,甚至计算到不收敛时塑性区仍未贯通,在用MIDAS计算时这一现象反映得更加明显。     

边坡稳定强度折减有限元分析中的若干问题讨论

强度折减有限元法在边坡稳定性分析中得到了广泛应用,但仍存在一些问题未能很好地解决。首先探讨了不同单元类型对边坡稳定安全系数计算精度的影响,指出采用三角形二阶单元,能够较好地模拟土坡的渐近破坏过程,且相对网格大小控制在0.1～0.15之间时能获得较精确的解;其次讨论了边坡失稳的4个判据问题,指出在合理网格大小区间内,有限元强度折减至土坡破坏时,先后出现塑性区贯通、等效塑性应变带贯通、位移陡增与数值计算不收敛等4种表征,以此分别作为边坡失稳判据,所得到的安全系数具有较好的一致性;最后讨论了边坡破坏时滑动面深浅、滑出点位置与坡体材料的黏聚力、内摩擦角及坡角间的关系,指出黏聚力越大、内摩擦角越小或坡度越小,滑动面越深,且滑出点越远离坡脚。     

边坡稳定有限元分析中的三个问题

对大型有限元软件ABAQUS进行二次开发,提出了自动搜索安全系数的边坡稳定有限元分析模型。在此基础上进行大量变动参数研究,探讨了迭代不收敛、塑性区贯通及等效塑性应变贯通等3种边坡失稳判据的内在联系与适用性,其中迭代不收敛判据易于自动搜索的编程实现,且较少依赖研究者的经验。大多数情况下单元阶次不影响安全系数的确定,但一阶单元有时可能高估安全系数,建议采用二阶单元。基于更新拉格朗日格式实施边坡稳定大变形有限元分析,结果表明迭代不收敛准则不适用于大变形分析。     

物质点强度折减法及其在边坡中的应用

物质点法适用于模拟连续介质大变形,如边坡失稳全过程。在物质点法中应用强度折减法,用于边坡稳定性评价。与极限平衡法相比,二者安全系数计算值、滑动面位置结果基本一致;与有限元强度折减法相比,物质点法失稳评价标准的物理意义明确。利用物质点法大变形计算优势,评价边坡失稳后的破坏后果,通过算例说明其评价不同安全系数下的滑坡堆积形态及滑移距离的能力,尤其是评价滑坡对临近建筑物的影响程度的能力。物质点强度折减法可用于边坡稳定性评价及边坡破坏后果评价。     

浅埋隧道掌子面稳定性二维自适应上限有限元分析

为研究浅埋隧道掌子面稳定性及获取精细化的破坏模式,提出了一种上限有限元非结构化网格自适应加密策略。以单元耗散能权重指标作为网格自适应加密评判准则,该策略同时兼顾了单元尺度与塑性应变。应用高阶的6节点三角形单元并建立上限有限元线性规划模型,以多次反复计算和网格加密的方式实现了二维自适应上限有限元分析并编制了计算程序。利用条形基础地基极限承载力课题,从上限解精度和网格加密形态方面验证了该程序的有效性。针对浅埋隧道掌子面稳定性问题,展开多参数条件下的自适应上限有限元计算,分析了网格加密过程中单元总数与上限解精度的关系,列出不同隧道埋深和内摩擦角对应的隧道掌子面稳定性临界值的上限解,揭示出掌子面稳定性变化规律及精细化的破坏模式。     

基于M-P法的抗滑桩支护边坡稳定性分析

为评价采用抗滑桩支护后的边坡稳定性,从设计和验算两个角度着手,基于Morgenstern-Price（M-P）法建立了抗滑桩支护边坡的分析模型,进而得到了抗滑桩下滑力和边坡安全系数的表达式;通过引入自适应遗传优化算法,建立边坡稳定性分析优化模型,搜索采用抗滑桩支护边坡的非圆弧最危险滑动面;从确定安全系数求抗滑桩所受下滑力和确定抗滑桩抗力求边坡最小安全系数两种情况出发,探讨了抗滑桩位置对边坡稳定性的影响。结果表明：该方法可搜索出更加符合实际情况的边坡非圆弧最危险滑动面,并且能够得到抗滑桩的下滑力或边坡的最小安全系数,相同条件下抗滑桩应设置在边坡中部较为适合,才能最大发挥抗滑桩的加固作用。     

边坡稳定及加固分析的有限元强度折减法

将强度折减理论应用于土质边坡稳定及加固的弹塑性有限元分析中,提出了以某一幅值的总等效塑性应变区,从坡脚到坡顶贯通时为边坡破坏的标准,不但物理意义很明确,而且可以根据边坡临界破坏时的等效塑性应变区确定滑动面位置,并认为此前的折减系数即为边坡的安全系数。对某供水改造工程中边坡的稳定性分析及加固计算表明：该破坏判别标准是适宜的,能够在工程实际中加以运用。     

基于动力学和材料软化特性的边坡渐进破坏特征研究

边坡的失稳是一个从量变到质变的动态渐进破坏过程,此问题也是边坡领域研究的重点与难点之一。在考虑岩土材料软化特性和动力学求解的基础上,建立了边坡渐进破坏仿真的理论框架;利用ABAQUS软件的动力显式求解模块实现了边坡的渐进破坏仿真;根据塑性应变揭露了剪切带的扩展过程,由软化本构确定了滑面材料的分区演化规律,根据等效塑性应变确定了边坡的滑面,通过滑面位置将边坡分为滑体、滑带、滑床,并分别研究了边坡各分区内部特征点运动学变量的发展过程,从而揭示了边坡的渐进破坏过程;基于材料参数沿滑面的时空分布,利用矢量和法得到了边坡不同演化阶段的安全系数。对比该方法与Bishop法确定的滑面位置与安全系数,发现两种方法峰值和残余强度对应的安全系数比较接近,该方法搜索所得滑面位于Bishop法自动搜索的滑面之间,验证了此方法的合理性及可靠性。最后分析了材料软化特征对边坡稳定性的影响,在保持其他参数不变的条件下,增大残余黏聚力,边坡的滑面位置加深,安全系数的初始值减小,安全系数的快速减小阶段有所推迟,并且快速减小阶段经历的时间有所延长,稳定后的安全系数有所增大。保持其他参数不变,增大残余黏聚力对应的等效塑性应变阈值,边坡的滑面位置加深,安全系数的初始值减小,安全系数的快速减小阶段有所推迟,但快速减小阶段经历的时间基本不变,达到稳定的时间有所推迟,同时稳定后的安全系数略微有所增大。     

基于最大剪应变增量的边坡潜在滑动面搜索

滑动面位置和形状的判断是影响边坡地基稳定性分析合理性的主要因素,合理地确定滑动面的位置和形状是边坡地基稳定性分析首先要解决的问题。由岩土工程常用的Mohr-Coulomb度理论可知,岩土体破坏是由于某一面上的剪应力达到了岩土体的剪切强度,此时剪切面上必然发生较大的剪切变形,因此可以通过搜索最大剪应变增量的位置来确定边坡滑动面。基于以上思路,综合运用有限差分软件FLAC3D和曲线拟合技术,提出基于最大剪应变增量确定边坡地基潜在滑动面的方法。首先,利用有限差分软件FLAC3D计算极限状态时边坡内的应力、应变,在计算剖面中设置一系列的垂线,利用FISH语言编制程序搜索垂线上剪应变增量最大位置的坐标,以获得滑动面上的离散点坐标。然后,通过最小二乘的方法对离散点进行曲线拟合,得到滑动面的形状和位置。通过与极限平衡法的对比分析,验证文中方法的合理性。文中方法的判断标准物理意义明确,简便实用,且克服了对网格的依赖性。