沉入式大圆筒码头结构稳定性分析

沉入式大圆筒码头结构稳定性分析

大直径截面结构（以下简称大直径）是一种新型港口水工建筑物结构。成本低，施工周期长。适合于软土基和恶劣海况的优点，具有广阔的应用前景。然而,沉入式大圆筒结构受力条件及其作用机理复杂,目前尚缺乏广泛认可的理论体系与计算方法,稳定性分析方法无规范或标准可循。随着有限元技术的发展,利用有限元数值模型分析新型结构力学特性成为一种有效的途径。对于沉入式大圆筒等新型防波堤结构,由于承受的外载荷主要是波浪力,针对这一特点,近年来提出了新型防波堤结构稳定性分析的有限元加载系数法。其基本思想是,在有限元模型计算过程中,逐步增加波浪力,直至结构发生稳定性破坏,此时的波浪力即为结构的极限承载力。为表述清楚,定义一个表征荷载加载程度的加载系数:α=P/PDα=Ρ/ΡD式中:P为加载波浪力,PD为设计波浪力。当加载到结构极限承载力Pu时,加载系数即定义为结构稳定性安全系数K。文献利用加载系数法计算分析了沉入式大圆筒防波堤、箱筒型基础防波堤、半圆型防波堤、倒T型导管墙桩基防波堤等新型防波堤结构的稳定性,并在此基础上进一步分析土与结构相互作用关系,建立更为合理有效的简化计算方法等。然而,对于沉入式大圆筒等新型码头结构,其承受的外部荷载复杂,除了波浪力作用,还有墙后填料产生的土压力,筒前低水位和筒后地下水位差引起的剩余水压力,船舶荷载,堆货荷载、流动机械荷载等各种地面使用荷载以及码头面可变荷载作用产生的土压力等。此外,对码头结构进行稳定性分析时,必须考虑各种不利情况下的荷载组合,显然,利用加载系数法进行码头结构的稳定性分析变得十分困难。强度折减法是通过土体强度参数的降低来实现对失稳情况的模拟,近年来被广泛应用于边坡稳定性分析中。文献采用强度折减法进行边坡稳定性分析,以非线性有限元静力计算不收敛时的折减系数作为边坡稳定性安全系数,并得到边坡破坏时的滑动面。算例表明采用摩尔—库仑等面积圆屈服准则时,强度折减法求得的稳定安全系数与传统计算方法非常接近,证实了该方法的实际可行性。文献针对某一典型边坡算例,分别依据计算的收敛性、特征部位位移的突变性和塑性区的贯通性等3个失稳判据,采用强度折减法进行边坡稳定性分析。通过对比Spencer极限平衡法的计算结果,建议联合采用特征部位位移的突变性和塑性区的贯通性作为边坡的失稳判据。文献将强度折减法应用于边坡稳定性的三维分析,通过对比传统极限平衡法和二维分析的结果,验证了三维分析的可行性,认为采用有限元法对边坡稳定性进行分析时,二维分析难以较好地模拟边坡复杂的地质、地形、地势等实际条件,因此,在有条件的情况下,尤其是对大型边坡和重要工程,宜补充进行三维分析。文献认为土体强度参数c、φ在边坡失稳破坏时的发挥程度、发挥秩序以及衰减速度、衰减程度并非完全相等,提出了双系数折减法,即在强度折减过程中对c、φ采用不同的折减系数,而不同的折减方法将对应不同的安全系数和滑动带,使边坡的可能滑动位置变成一个可能滑动区域。沉入式大圆筒结构的抗滑、抗倾稳定性靠地基土的嵌固作用来维持,在极限状况下,认为外载荷所产生的土体实际剪应力与土体抵御外载荷所能发挥的最低抗剪强度相等,即,当不断折减土体的实际强度指标直至达到极限破坏状态时,外载荷所产生的剪应力与强度折减后所确定的土体抗剪强度相等,可通过弱化土体强度来实现对结构极限破坏状态的模拟。基于上述假设,并在沉入式大圆筒防波堤结构工程算例验证的基础上,建立了沉入式大圆筒码头结构稳定性分析的有限元强度折减法。结合工程实例,利用有限元软件ABAQUS,采用有限元强度折减法计算沉入式大圆筒码头结构的稳定性安全系数。1土体抗剪强度参数的计算有限元强度折减法是强度折减技术与弹塑性有限元方法的结合,其基本原理是:在有限元计算中,逐步折减土体的抗剪强度参数——粘聚力C和内摩擦角的正切值tanφ,直至土体达到极限破坏状态,定义强度折减法的安全系数Fv为:Fv=τ/τ′Fv=τ/τ′式中:τ为土体的初始抗剪强度,τ′为强度折减后极限状态下的土体抗剪强度。不同的屈服准则下,土体抗剪强度的表达式不同,以摩尔—库伦准则为例:τ=c+σtanφ。因此,强度折减过程可表示为:τ′=τFv=c+σtanφFv=c′+σtanφ′c′=cFvtanφ′=tanφFvτ′=τFv=c+σtanφFv=c′+σtanφ′c′=cFvtanφ′=tanφFv可见,强度折减法对于安全系数的定义与传统极限平衡法一致,亦是一种强度储备安全系数。在传统极限平衡法下,当土体达到极限状态(即安全系数为1)时,破坏面上的土体抗剪强度参数C和tanφ值降为c′和tanφ′。因此,强度折减法中的土体折减系数即为安全系数Fv。综上,采用有限元强度折减法进行稳定性分析,就是在有限元计算过程中,将土体的强度参数——粘聚力C和内摩擦角的正切值tanφ折减一个系数,生成一组新的土体参数,再代入有限元计算模型进行试算,以位移、广义剪应变等描述变形程度的某种物理量作为评判指标,判断土体是否达到极限破坏状态。反复进行上述步骤,当土体达到极限平衡时,其强度参数的折减系数,就是稳定性安全系数值。2项目总结2.1基桩式长桩方案某港码头结构方案之一是沉入式大圆筒结构。码头顶高程为9.50m,考虑构建5m高的胸墙,圆筒结构顶高程设为4.50m,筒底入土至-50.00m。采用钢质圆筒,直径26m,壁厚20mm,结构平面如图1所示。港区原泥面高程为-10.00m,码头前沿开挖至-19.50m。大圆筒筒内及筒后回填中粗砂,码头后方采用吹填砂。考虑到建筑物上部结构自重较大,门机等荷载也较大,该方案在码头前方胸墙及轨道梁下加设了部分长桩,将上部荷载直接传入下层地基,以减小结构变形。长桩底端入土至-60.00m,码头结构断面如图2所示。在结构稳定性分析中,由于基桩的存在有利于沉入式圆筒码头结构的稳定,为使计算偏于安全,计算中暂不考虑基桩的作用。简化后的大圆筒结构计算断面如图3所示。计算工况采用设计高水位+50年一遇波浪作用,设计波浪按S方向考虑;剩余水头采用1/3平均潮差;堆货荷载取用20kPa;筒顶承受胸墙等码头上部结构的自重作用。2.2土壤指数各层土性指标见表1。2.3墙前产生立波的影响设计水位和波浪条件如表2和表3所示。连续圆筒结构上的波浪力,按平面直墙波浪力计算方法确定,波峰作用时考虑0.9左右的折减系数,波谷作用时考虑0.95左右的折减系数。计算中为保守考虑取折减系数为1,且考虑最不利荷载——波吸力的作用。码头前沿泥面标高为-19.50m,根据波浪条件,码头前沿水深d=22.41m‚Tg/d−−−√=5.02<8d=22.41m‚Τg/d=5.02<8且d>H,故墙前产生立波。根据文献相应的计算公式,得静水面处波浪压强p0=0kPa,静水面以下深度H-hs处波浪压强ps=41.62kPa(hs为波浪超高),水底处波浪压强pd=18.90kPa,墙面波浪压强分布如图4所示,单位长度墙面上总波浪力P=639.68kN/m。3有限分析模型3.1墙后土体弹塑性整体模型的建立大圆筒码头结构由一系列圆筒薄壳连续排列构成,具有连拱形墙面,空间承载特性明显。采用有限元方法,利用有限元分析软件ABAQUS,建立大圆筒码头结构与土相互作用的三维弹塑性整体分析模型。为便于分析,假设波浪荷载的作用方向垂直于墙轴线,根据荷载及结构的空间对称性,在垂直于墙轴线方向取一个圆筒的一半和前后一定范围的土体作为建立有限元模型的区域,在边界上施加对称边界条件。其中,土体的计算域为:在垂直防波堤轴线方向,向两侧各取大圆筒直径尺寸的5倍;在深度方向,筒体以下取3倍的圆筒入土深度。整体计算域平面示意和立面示意分别如图5和图6所示。计算域边界条件为:地基表面为自由边界,底面为固定边界,前侧面和后侧面为侧限边界,左侧面和右侧面为对称边界。3.2地基土体本构模型由于大圆筒结构由钢材制成,结构强度和刚度远远大于土体强度和刚度,结构系统的位移和失稳破坏主要决定于地基土的变形和承载能力,故在有限元分析中大圆筒结构采用弹性模型,计算中大圆筒结构参数取用弹性模量E=210GPa,泊松比μ=0.3。土体本构模型采用C-M模型。在建立的三维有限元模型中,土体和结构均采用八节点三维减缩积分实体单元。3.3土体接触面的本构模型在结构物与地基相互作用问题中,由于和土体材料相比,结构材料具有很高的强度和刚性。因此,在应力作用下,接触面上土体和结构难以协调变形,容易产生接触面的分离和相对错动,需要在结构与土体接触的区域(包括大圆筒的内外壁及底部与土体接触的区域)引入接触面单元,建立主从接触面,以模拟在波浪力、剩余水压力、堆货荷载等外载荷作用下结构与周围土体的粘结、滑移和脱离现象。考虑到结构的弹性模量远大于土体的弹性模量,指定结构上的接触面为主接触面,土体上的接触面为从接触面。土体与结构的接触行为是一个高度非线性的问题,在接触面上需要确定不同时刻接触区域不同位置处的动接触状态以及接触行为的本构模型。在建立三维弹塑性耦合分析模型时,接触单元的本构模型切向采用库仑摩擦本构模型;法向采用硬接触方式,即假定相互作用的两个物体,相互挤压时不产生浸入变形,受法向拉力时在接触面上产生分离。土与结构相接触区域的不同位置处的接触状态是随着时间不断发生变化的,分为三种接触状态,即分离状态、粘结接触状态及滑动接触状态。当处于分离状态时,接触力为零;当处于粘结接触状态时,主从接触面相互挤压,但无相对滑动,摩擦力处于静摩擦状态。3.4稳定分析3.4.1有限元分析在土体强度参数确定中的应用基于强度折减法原理,采用有限元方法进行波浪力、剩余水压力、堆货荷载等外载荷作用下沉入式大圆筒码头结构稳定性分析的步骤如下:1)建立有限元分析模型。2)根据波浪力、剩余水压力、堆货荷载、上部结构自重等外载荷的分布,在结构上相应位置施加设计荷载,计算在设计荷载作用下土与结构体系的位移和应力。3)逐级降低有限元分析模型中全部土体的强度参数(包括粘聚力C和内摩擦角的正切值tanφ),计算相应土体强度下土与结构体系的位移和应力,直至有限元计算不收敛。4)绘制折减系数——结构位移关系曲线。5)确定依据折减系数——结构位移关系曲线判别结构失稳的准则。6)根据失稳判别准则,确定结构稳定性破坏对应的折减系数,将该系数定义为结构稳定性安全系数。3.4.2强度折减刑传统稳定性法虽然强度折减法的基本原理简单,近年来也在边坡稳定性分析中得到了广泛的应用,但目前尚缺乏统一的失稳评判标准。同时,由于边坡失稳与结构失稳的模式不同,现有的强度折减法边坡失稳判据也不能完全适用于判别大圆筒码头结构的稳定性。考虑到强度折减法的基本思想与传统的极限平衡法一致,均可称之为强度储备安全系数法,即主要考虑的是力和强度之间的关系问题。结合结构极限承载力判别标准,以折减系数——结构位移关系曲线(即FV-S曲线)斜率接近于零时对应的折减系数为结构的稳定性安全系数,此时即使增加微小的折减系数,结构也将发生非常大的变位,说明结构已经破坏。该判别准则精确、可靠,可避免由于其他原因造成有限元计算不收敛而导致的判别错误。4有限概率分析4.1层压高和结构稳定性为了验证本文所建立的结构稳定性分析有限元强度折减法的正确性,将其应用于某一沉入式大圆筒防波堤工程实例的稳定性分析,并将计算结果与防波堤稳定性分析的有限元加载系数法进行对比。沉入式大圆筒防波堤采用钢质圆筒,直径22m,壁厚16mm,筒顶标高8.50m,筒底入土至-49.00m。港区原泥面高程为-19.00m,筒内泥面以上回填块石。各层土性指标见表4。钢质圆筒材料参数取用弹性模量E=210GPa,泊松比μ=0.3。图7、图8分别给出了所选沉入式大圆筒防波堤结构在计算工况1:极端高水位+50年一遇波浪作用,结构顶高程为8.50m和计算工况2:设计高水位+10年一遇波浪作用,结构顶高程为4.00m条件下,强度折减法计算所得折减系数——结构转角关系曲线和加载系数法计算所得波浪力加载系数——结构转角关系曲线。根据所建立的结构稳定性分析强度折减法的失稳判别准则和有限元加载系数法常用的极限承载力判别标准,将强度折减法与加载系数法的结构稳定性分析计算结果列于表5。由对比结果可以看出,有限元强度折减法与加载系数法对于大圆筒防波堤结构的稳定性分析结果一致性良好,两种方法计算得出的结构稳定性安全系数基本相同。说明有限元强度折减法适用于结构的稳定性分析,且具有较高的精度,验证了所建立的有限元分析方法的正确性。4.3特征点折减系数—特征点的选取为便于分析,选取迎浪侧圆筒结构顶部端点及码头前沿泥面处筒体内壁上的结点作为分析特征点,如图9中A点、B点所示。通过特征点的折减系数—位移关系曲线可以方便地得到结构变位并以此判断结构的稳定性。由于结构整体刚度远大于地基刚度,可以看成刚体。定义A、B两点间水平位移差与两点间距的比值的反正弦值为结构转角θ,即:θ=arcsin[(UA−UB)/L]θ=arcsin[(UA-UB)/L]式中:UA、UB分别表示A点、B点的水平位移,L表示A点和B点之间的距离。4.3固结快剪指标强度图10和图11分别给出了本文所建立的沉入式大圆筒码头结构有限