考虑稳定渗流和波浪力作用下内河航道岸坡稳定性分析.docx

考虑稳定渗流和波浪力作用下内河航道岸坡稳定性分析 河流底板的不稳定和破坏是世界上所有河流边坡的现象。河流上建设枢纽后, 航道两岸塌岸更为严重。航道岸坡失稳引起沿岸地带的建筑物破坏, 造成财产损失;同时, 岸坡失稳将侵蚀两岸土地资源, 危害航道通航条件, 也极大增加了航道的维护费用。因此, 开展航道护岸稳定的研究具有十分重要的意义。 影响航道岸坡稳定性的因素很多, 包括渗流、降雨、波浪力等, 本文主要研究稳定渗流和波浪力对岸坡的影响。在渗流方面, 陈丽刚基于ABAQUS对边坡稳定渗流进行了流固耦合分析。张晓咏应用ABAQUS对坝体的渗流进行了研究, 并验证了强度折减法的可靠性。周群华、章广成等研究了水位变化对边坡稳定性的影响。在船行波方面, 项菁等对船行波波要素的特点以及现有的一些计算公式进行了总结和分析。卢无疆对高速双体客船和普通船的船行波进行了现场观测, 得到了各种波浪要素。简文彬对边坡对循环荷载的响应进行了研究。 虽然不少学者已经对边坡失稳的影响因素进行了研究, 但都只是针对单一影响因素研究, 缺乏对多种因素共同作用下边坡失稳的研究, 实际上边坡失稳是多种因素共同作用下的结果, 因此多因素共同作用下边坡的稳定性研究也是必不可少的。本文除了分别分析稳定渗流和波浪力这2个影响因素之外, 还将其组合起来, 研究这2种因素共同作用对边坡稳定性的影响。 本文采用大型有限元分析软件ABAQUS, 通过数值模拟对岸坡进行流固耦合分析, 并采用强度折减法分析岸坡的稳定性。在单独考虑渗流时, 先计算不同静水位下边坡的安全系数, 然后在数值模型上施加不同的水头差, 计算稳定渗流情况下的安全系数, 最后与极限平衡方法的计算结果相对比。在单独考虑波浪力时, 参考实测的波浪要素, 根据规范计算出波压力的分布, 并在模型上施加周期荷载, 研究不同波高、周期的船行波对边坡稳定的影响。在考虑共同作用时, 将之前计算的工况组合起来, 研究2种因素共同作用下对岸坡稳定性的影响, 并与单因素的计算结果进行对比分析。 1 基本原则 1.1 土体渗流方程的建立 ABAQUS可对多孔介质的渗流和变形进行耦合分析, 可进行非饱和土的渗流计算。ABAQUS中处理孔隙介质中的流体流动的方式是将孔隙体视为多相材料, 孔隙中的流体可包含两部分:一是液体, 通常认为压缩性相对很低;另一个则是气体, 认为是可压缩的。土体的体积包括两部分:土颗粒的体积和孔隙的体积, 孔隙的体积等于孔隙中液体的体积与气体体积之和。计算中有限元的网格固定在土骨架上, 气体或液体可流过网格, 但需要满足流体的连续性方程。土体的力学特性通过采用有效应力定义的本构模型来模拟, 液体的渗透采用Forchheimer渗透定律模拟, 常用的Darcy定律是它的简化。 式中:δv为虚速度场;δε=sym (坠δv/δx) 为虚应变;σ为真实 (Cauchy) 应力;t为单位面积上的表面力;为单位体积上的体积力。 渗流连续方程为 式中:vw为流体的平均流速;n为边界S的外法线;该方程对ρw0进行了归一化。 Forchheimer渗透定律表达式为 式中:g为重力加速度;为土体的渗透系数;β (e) 为“速度系数”。当β=0时, 式 (3) 即为Darcy定律。 ABAQUS能将渗流场和应力场直接藕合, 无需进行渗流场和应力场的反复迭代, 只要按时间过程连续求解就可得到全部结果, 即通过将节点位移和孔隙水压力作为节点自由度进行空间离散, 将应力平衡方程和渗流连续方程写成矩阵形式, 并对渗流连续方程引入时间积分, 得到耦合控制方程, 然后在每个时间步内求解方程, 并同时满足位移边界条件和渗流边界条件。 1.2 有限元分析方法 边坡稳定的分析方法主要可分为两大类:第一类可称作传统方法, 如极限平衡法、极限分析法、滑移线法等;第二类为基于计算机技术的数值分析方法, 如有限元法、离散元法、边界元法等。对于稳定渗流的土坡稳定分析, 传统的极限平衡条分法分成2种方法:一种是根据边界孔隙水压力和饱和的重度考虑渗流作用的常用方法;另一种是根据渗透力和浮重力考虑渗透作用的分析方法。前者需要绘制流网求得孔隙水压力, 而后者由于计算渗透力很不方便, 工作量较大, 在工程实践中较少采用。本文采用的有限元强度折减法, 能够很好地克服极限平衡法的缺点, 因而已经被广泛的应用到了有数值分析中。 抗剪强度折减系数 (SSRF:Shear Strength Reduction Factor) 的概念是由Zienkiewicz等1975年在《土工弹塑性有限元数值分析》一文中首次提出的, 其定义为:在外载荷保持不变的情况下, 边坡内土体所能提供的最大抗剪强度与外载荷在边坡内所产生的实际剪应力之比。在极限状况下, 外载荷所产生的实际剪应力与抵御外载荷所发挥的