高桩码头三维有限元分析.docxVIP

高桩码头三维有限元分析 岸壁的稳定一直是工程和学术界的研究主题。传统的圆弧滑动极限平衡方法使用简便,在工程中得到了广泛的应用,但是也有一定的局限性。首先,它将土视为刚性体,只能得到抗力分项系数,无法获得土体的应力和变形。其次,它不能计算存在桩基的岸坡稳定安全系数,在进行岸坡设计时,不考虑桩基的抗滑作用,这显然是偏保守的。此外,在高桩码头使用过程中,桩基的变位、桩和梁间的错位等都是经常出现的现象,极大地影响了岸坡稳定,而这些问题是传统的极限分析方法无法解决的。与极限平衡法相比,有限元强度折减法能够自动地求得任意形状的临界滑移面以及对应的最小抗力分项系数,还可以真实地反映坡体失稳及塑性区开展的过程,因此,在岸坡稳定分析中的应用备受重视。目前已有很多学者对有限元强度折减法进行了深入的研究,例如,Griffiths等,Dawson E.M.,Yuan-Liang Chang等,连镇营,张鲁渝、郑颖人等学者进行了研究工作,推导出了基于不同强度准则的计算方法。 笔者根据珠海某煤炭码头勘察报告提供的土体参数资料,建立了码头结构与岸坡三维有限元模型,并利用强度折减理论对其进行稳定性分析。 1 数值分析理论 1.1 抗力分值系数的计算 Duncan将安全系数(现行规范中称为抗力分项系数)定义为使岸坡达到临界破坏时对土的强度指标进行折减的程度,也就是将土体强度参数c,φ等同时除以折减系数F,得到一组新的c,φ值,然后作为新的材料参数输入有限元软件进行试算,当计算正好收敛时,即再稍大一些计算便不收敛,此时对应的F称为坡体的临界折减系数,此时土体达到临界状态,发生剪切破坏。可以用公式作如下表示: 有限元强度折减法无法直接用公式计算抗力分项系数,而需根据某种破坏判据来判定系统是否进入极限平衡状态,这样不可避免地会带来一定的人为误差。尽管如此,仍发展了一些切实可行的平衡判据,如:限定求解迭代次数(如500次),当超过限值仍未收敛则认为破坏发生;或利用可视化技术,当广义剪应变等值线自坡角与坡顶贯通则定义坡体破坏;或限定节点不平衡力与外荷载的比值大小。根据赵尚毅和郑颖人的研究,塑性区从坡角到坡顶贯通是岸坡破坏的必要条件而非充分条件。综合考虑,采用有限元计算不收敛作为岸坡失稳的判断依据,迭代次数10 000次,力收敛标准为l0-5,材料弱化系数步长为2.5%,依次降低弱化系数的数值直至计算不收敛时,则前一次收敛的数值即为临界弱化系数。 1.2 平面屈服准则 用折减系数法求解实际土坡稳定问题时,通常将土体假设成理想弹塑性体,其中本构模型常选用Mohr-Column准则、Drucker-Prager准则等。传统边坡稳定分析的极限平衡条分法采用的是Mohr-Coulomb准则,该准则较为可靠,在土体计算中应用十分广泛。但因为Mohr-Coulomb准则在三维应力空间中不是一个连续函数,是由6个分段函数所构成,在三维应力空间的屈服面为不规则的六角形截面的角锥体表面(图1),在π平面上的图形为不等角六边形,其存在尖顶和菱角,给数值计算带来困难。 Drucker-Prager屈服准则在主应力空间的屈服面为光滑圆锥面,在π平面上为圆形,不存在尖顶处的数值计算问题。因此,选用与Mohr-Coulomb匹配Drucker-Prager准则作为土的屈服准则,以利于数值计算。Drucker-Prager准则表示为: 式中:I1,J2分别为应力张量的第一不变量和应力偏张量的第二不变量。α,k是与岩土材料内摩擦角φ和粘聚力c有关的常数,不同的α,k在π平面上代表不同的圆。各屈服准则在π平面上的曲线如图2。 本文采用莫尔-库仑等面积圆D-P,此准则1990年由徐干成和郑颖人提出,并成功地运用于边坡稳定性分析中,换算关系如下: 本文采用关联流动法则,即膨胀角和内摩擦角取同值。 1.3 扩展d-p准则 ABAQUS中扩展的Drucker-Prager模型可以是理想塑性和强化塑性,其子午线形状可以是直线、双曲线和一般的指数函数,具有很强的适应性。本文采用子午线为线性性状的扩展D-P准则。有如下屈服函数: 式中:β为摩擦角;p=-1/3I1;凝聚力d和输入的硬化参数有关。使用材料单轴受拉屈服应力σt定义硬化时: 比较式(1)和式(2)可得: 则可以通过4步求出ABAQUS中等面积圆屈服准则的材料参数内摩擦角β和屈服应力σt: 1)首先将土体材料参数凝聚力c和摩擦角φ代入式(3),求得一组α、k; 2)然后将α,k代入到式(6)求得β和d; 3)将tanβ,d代入式(5)求出σt; 4)将求得的β和σt输入ABAQUS即可。 2 有限元模型及计算 选取珠海港某10万吨级煤炭码头建立模型。码头平台宽度为35.5 m,由桩基、及上部结构结构组成,码头通过引桥同后方陆域连接。排架间距为10 m,每个排