大型海洋平台基础结构有限元分析

1、16卷6期2007年12月自然灾害学报JOURNAL OF NAT URAL D IS ASTERS Vol .16,No .6Dec ,2007收稿日期:2007-06-10;修订日期:2007-09-20作者简介:徐秀丽(1963-,副教授，硕士，主要从事桥梁结构分析及振动控制研 究.E 2mail :ujxuxiuliyalioo .文章编号：100424574(20070620209207大型海洋平台基础结构有限元分析徐秀丽1，张元壮1，刘伟庆1，潘晓峰1，王仁贵2(1.南京工业大学土木工程学院，江苏南京210009; 2.交通部公路规划设计院，北 京 100010摘要:某大型海洋平台

2、基础采用平台板+桩基础组合形式，承台面为椭圆形，台面 尺寸为147.6m x9910m,平台板为混凝土板和钢梁组成的叠合梁板结构。该工程地处世界三大 强潮海湾之一,所受荷载与当地地质条件均较复杂。介绍了应用有限元对其进行结构分析的实施方法 及过程，采用弹性抗力法考虑了桩-土共同作用，重点探讨了钢梁加混凝土厚板组成的 叠合梁的有限元模拟方法。通过对比分析，提出了计算结果精度高、可操作性强的 叠合梁有限元模拟方案。有限元分析结果表明，该平台钢管混凝土桩与钢梁均满足 材料强度要求，但在施工阶段，桩内存在拉力。关键词:海洋平台基础;有限元法;钢-混 凝土叠合梁;结构分析中图分类号:U441;TU381

3、文献标识码:AF 1 n ite elem en t ana lyses of 1 aige 2si ze offshoie pl a tfbf m fbunda ti onXU Xiu 2h 1,ZHANG Yuan 2zhuang 1,L IU W ei 2qing 1,P AN Xiao 2feng 1,WANG Ren 2gm2(1.College of Civil Engineeimg.Naiy mg University of Teclmol ogy.Naiying 210009,Chma;2.H ighway Planning and Design I nstiniteM ims

4、tiy of Commuiucati ons of China,Beijing 100010,ChinaAbstract:A large 2size offshore p latf 01 m with asse mbled p latfbi m slab 2p lie f oundati on has an elli p tic shape p latf or m of 147.6m x99.0m di m ensi on .The p latf oi m slab is a super posed bea m 2slab stmctuie constructed by steel bea m

5、 and conciete slab .The gulf where the p latfbi m 1 ocates is one of the three biggest tide gulfs in the woildand the geol ogical conditi on there is comp lex .The stmctiue bears comp lex 1 oad .1 n tins paper,a p r ocess of static analysis fbf this stmctuie was inti oduced by finite ele ment method

6、 .The p lie -s oil inteiacti on p r oble m was s olved by elastic 2iesistance method,and the si m ulati on of composite steel bea m 2concrete slab was mamly discussed .The analysis method p i oposed in this paper is higli 2p lecise and high 2opeiable .The finite ele meat analysis result shows that t

7、hematenal stiength de mand of CFST （concrete filled steel tubes 2p lies and steel 2beam s can be satisfiedbut theie is pull in the p lies on the st ocks .Keywords:offshore p latfbf m f oundati on;fiinte ele nient method;steel 2concrete super posed bea m;stnicune analysis1工程概况某大型海洋平台是某跨海大桥工程的重要组成部分，通

8、过设置海中平台结构，一 方面可以作为大桥建设过程中施工、测量以及海上作业人员的工作基地，解决沿途 海域无可利用岛屿的困难，提高大桥施工的效率;另一方面在大桥竣工后则可成为一 个海中交通通讯服务、海事监控、海上紧急救援基地，同时也是一个绝佳的旅游观光平台（图1。海洋平台总面积约3.9万m2，共分为3个部分:桩基础及桩顶平台（停车、桩顶平VII 嘛叩 IAdf3OIO图1海洋平台建成后效果图Fig .1 Rendemig of the offshore p latfbi mt o be constnicted台以上屋面以下的上部结构（5层:餐饮、住宿、商务、屋面系统（网架结构。本工程所处建设场地自

9、然环境较为恶劣：（1地处世界三大强潮海湾之一,必须考 虑海流海浪对结构的影响;（2所在海域位置存在约16m的淤泥质可冲刷土层，地质条 件很差，海中平台采用桩基十桩顶平台的形式，桩-土相互作用显著;（3海洋平台所处海 湾海域辽阔,风力大，风期长，为台风影响和龙卷风出现多发区。在海中建造近4万m2的大型平台，加之所处场地的自然环境复杂，极大地增加了工程的建设难度，因此 必须对此平台进行全面的结构分析，以保证平台在各荷载工况下的安全性与使用性 能。海洋平台基础部分采用平台+桩基础组合形式，其中承台平面为椭圆形,台面标 高+11.8m （海平面518m，长轴半径147.6m,短轴半径99.0m,采用桩

10、柱一体钢管桩十 叠合梁结构并沿平台东西方向直桩部分设置交叉支撑，桩为直径1.6m的打入钢管桩 十填芯钢筋混凝土结构，共298根，其中直桩228根，斜桩70根;直桩、斜桩桩底标高均 为-80m,桩顶标高为+10.16m ,桩竖向投影长度为90.16m ,基础平台平面及剖面结构 见图2-4图2海中平台基础平面示意图（单位:cmFig .2Plan of the offshore p latf 01 m f oundati on (Umt:c mk3250J图3 A-A剖面示意图Fig .3 Pi ofile in A -Ap Ian图4 B-B剖面示意图Fig .4 Pi ofile in B -

11、B p Ian平台主梁与桩顶之间的连接采用可传递弯矩的刚节点，桩基础施工完毕后将节 点以及接长的钢管吊装焊接,并浇筑填芯混凝土。平台梁板结构由钢梁加铺预制混凝土板形成,其中平 台梁分为3类:A类梁,B类梁，次梁，均采用Q345钢材,具体形式见图5。为增加梁板 结构的整体性，除在钢梁上设置一定数量的连接件，还在各预制板间预留了 7cm的板 缝浇注现浇板带，最后在预制板板面浇筑8cm的整浇层（见图6,预制板混凝土及现浇混凝土均采用C60强度等级。考虑到该结构所处的特殊环境，需常年承担波流等 动荷载旦对钢筋混凝土的耐久性有较高的要求,在预制板的纵横向施加了双向预应 力。16012 自然灾害学报 卷图

12、 5 主梁与次梁截面（单位:nmiFig.5Ci oss seen ons of maui bea m and secondaiy.4类梁鳄图6叠合梁断面（单位:mmFig.6Ci oss seed on of super posed bea m（Umt:nun2有限元分析模型2.1钢管混凝土桩桩基础采用钢管混凝土桩形式，用梁单元模拟。建模时，根据混凝土与钢材的弹 性模量比值将钢筋混凝土填芯(C30与钢管桩换算为同种材料，钢管桩采用Q345级 钢，壁厚22mm,计算时取海水腐蚀后厚度17 mm。考虑桩与土的相互作用，采用桩基静力设计方法m值法，用水平弹簧模拟桩与土 的相互作用1。以-27.8m

13、为泥面标高(冲刷高度，计算出的水平弹簧刚度k如表1所 ZFo表1 m值法土弹簧刚度计算Table 1 Calculati on of s oil s p ring stiffness bymvahie method标高/m 土层厚度a/m计算宽度b p/ni抗力系数ni/(103kN - m-4 土层深度z/ni 土弹簧刚度 k/(103kN - m-1 -27.800-33.8006 2.3456421.2-39.8006 2.34512842.4-45.8006 2.345181263.6 -51.8006 2.345241684.8 -57.8006 2.3453021060-63.80

14、06 2.345362527.2 -69.8006 2.345422948.4 -75.8006 2.345483369.6 -80.0004.2 2.34552.22.2交叉支撑交叉支撑截面采用焊接工字形钢梁，钢材为Q345C级。根据交叉支撑的受力特 点，建模时用只受轴力的桁架单元模拟。2.3叠合梁根据结构设计方案，梁板结构为钢梁上加预应力混凝土板的叠合梁形式,本工程 中钢梁与混凝土翼板交接处设置一定数量的连接件传递接触面上的作用力,如果连 接件具有绝对刚性，则二者不会发生相对滑移，可称为完全共同作用的组合梁。由材 料力学可知，无共同作用和完全共同作用时,钢一混凝土组合梁截面的弹性抗弯刚度分

15、 别为2-5:E I=E c I c+E s I s(lEI=E I+EA 12(2 式中:EcIc,EsIs分别为混凝土翼板和钢梁截面的抗弯刚度;E I为无共同作用时，组合 梁”的弹性抗弯刚度工为混凝土翼板截面形心轴和钢梁截面形心轴之间距离;EA为截面组合抗拉 刚度，EA,A c,As分别为混凝土翼板和钢梁的截面面积;E I为完全共同作用组合梁截面的弹 性抗弯刚度。6期徐秀丽等:大型海洋平台基础结构有限元分析在建立有限元模型时，参照弹性理论为基础的容许应力设计体系6，采用如下假定:(1钢材和混凝土两种材料均认为是理想弹性体;(2钢梁与混凝土翼板之间有可靠的连接,相对滑移很小，可以忽略不计；(

16、3混凝土翼板按实体面积计算，不扣除其中受拉开裂的部分,为了简化计算.钢梁 托板面积忽略不计;(4不考虑混凝土板中的钢筋。如何有效模拟叠合梁是建立本有限元分模型的难题之一。本工程中钢梁均采用 两节点梁单元模拟，以Ti m osheiiko梁理论为基础，考虑梁剪切变形的影响,具有抗 拉、抗压、抗剪、抗弯、抗扭刚度;有限元分析中板有3种单元形式:薄板单元、厚 板单元、实体单元，由于厚板单元理论考虑了横向剪切应力的影响，因此较薄板单元 理论更为精确。由上述板、梁的模拟方法，在模拟叠合梁截面时，给出了如下5种方案，通过对比, 从中选择较优方案。方案1：使用实体单元建立模型；方案2:采用梁单元模拟工字钢梁

17、，混凝土板用厚板单元模拟，梁单元与板单元以 共用节点的方式连接;方案3:采用梁单元模拟工字钢梁;混凝土板用厚板单元模拟。 为了模拟板的真实位置，将板单元升至梁单元重心以上769.36mm (A类梁、 740.09mm (B类梁处，以梁单元的节点作为主节点，板单元对应节点作为从属节点进 行“刚体连接”；方案4:采用梁单元模拟工字钢梁;混凝土板用实体单元模拟。为了模拟板的真 实位置，将实体单元升至梁单元重心以上499.36nmi (A类梁、470.09nun (B类梁处， 以梁单元的节点作为主节点，对应的板底位置的实体单元节点作为从属节点，进行，刚 体连接；方案5:采用梁单元模拟工字钢梁;混凝土板

18、用实体单元模拟。为了模拟板的真 实位置，将实体单元升至梁单元重心以上499.36nmi (A类梁、470.09nun (B类梁处， 以梁单元的节点作为主节点，对应的板底与板顶位置的实体单元节点作为从属节点， 进行“刚体连接”。图7-9给出了方案3方案5中梁单元（A类梁与板单元或实体单元连接示意 图,B类梁类同650ViJKT*pA Hr)/ LnnV图7叠合梁模拟方案3Fig .7 Sche me 3of super posedbea mg 4 土 F扁I朱形心/650*301 (I| FVMMUK mnnA HBVft-Dvopue3 Imodelling图8叠合梁模拟方案4Fig .8 S

19、cheme 4of super posedbea m前66寸一刚性连接一主节点（梁形心）65()modelling图9叠合梁模拟方案5Fig .9 Sche me 5of super posedbeam modelliiig为了验证上述5种方案的可行性与有效性,本文给出如下算例进行比较2某两端固支叠合梁的跨度1 =6m ,混凝土板的计算宽度b e =1316nmi,混凝土板 的厚度h c =100mm ,混凝土强度等级为C 20;钢梁采用125a工字钢，钢号为Q 235。 梁在跨中承担竖向节点荷载100kN。按上述5种模拟方法分别对此叠合梁进行数 值分析（M 1-M 5分别代表方案1-5,分析

20、结果见图10-11。根据以上分析结果，并结合海洋平台基础的实际情况，对各模拟方案的优势和局限性作如下分析：(1方案1用实体单元建立模型，原理上应为精确的分析模型，但若采用此方案建 整个海洋平台基础模型，则建模过程复杂，计算量庞大，实现较为困难，且工程上梁的荷 载效应组合习惯采用内力,如果采用全部实体单元模型，则只能得到梁应力，无法获得 截面内力；(2方案2以梁单元与板单元共用节点的方式建立模型,建模过程简单，计算量小, 但此方案未能考虑板形心对梁形心相对距离所形成的空间刚度(以下简称“空间刚度作用”，致使计 算出的梁位移与应力结果偏大；212 自然灾害学报16卷(3方案3用梁单元与板单元建立模

21、型，通过模拟板与梁的实际位置考虑空间刚 度作用，建模过程简单计算量小，但是在台面不规则处,一些区域板单元已不符合板的尺寸比假定，故 此方案不能有效模拟形状不规则的叠合梁台面。对于形状规则旦厚度较小的板面， 建议使用此方案建立模型；(4方案4用梁单元与实体单元建立平台模型，通过模拟板与梁的实际位置考虑 了板的空间刚度作用，建模过程较方案2及方案3复杂，计算量也相应增大，但实体单 元能够较好地模拟平台面板不规则处区域，此方案中，钢梁与板之间的连接较实际情 况弱(实际情况为板面中心与梁中心相连计算结果偏安全，本文采用此方案建立平台 基础模型；(5方案5与方案4类似，但钢梁与板之间的连接较实际情况强，

22、计算结果略小偏 于不安全，平台基础也可采用此方案建模节点坐标/mm图io梁位移比较Fig .10 Compans on of dis p lace ments of bea.k(I -O-M4I.一.,./.01000 2000 3000 4000 5000 6000节点坐标/mmm图11梁底应力比较Fig .11 Compans on of stresses at bott om of bea m2.4钢筋混凝土板般在用实体单元模拟板时，实体单元是由面经过拖拉、旋转、偏移等方式生 成的,可先在原始面上生成板单元形式的网格，然后在生成体的同时自动生成三维实体网格。网格划分是 建模中较关键的一个

23、步骤，网格划分的好坏直接影响到计算的精度和速度一般应参照以下几点原则7-8(1四边形单元精度较三角形单元精度高，三角形单元建议只用作过渡；(2单元的纵横比不宜过大,如对于四边形，纵横比宜控制在14时,才能获得较好 的计算结果;(3四边形单元内角在90。左右，或至少在45。135。之间时，可获得较好的 计算结果;(4在集中荷载与形状突变或梁板相交处,网格应细分。3荷载般结构承担的荷载可分为5类：(1位移和约束位移和约束是用来定义模型固定的位置(零位移位置，也可以用非零位移来设定 个已知的偏移。本工程中将桩底的竖向自由度约束为零，以模拟桩支撑于大地。(2波流荷载及风荷载本工程中，波流荷载与风荷载需

24、要用线荷载形式加至相应的单元上。1波流荷载对于小尺度圆形构件(D /L 0.2,垂直于其轴线方向的单位长度的波浪力f,可按 Mons on 公式计算9-10:f=0.5C d pD |u |u +0.25C m p;rD 2u (3式中:C d ,C m分别为曳力系数和惯性系数;p为海水密度:D为圆形构件直 径;L为设计波长;u, u分别为水质点垂直于构件的速度和加速度。海流与波浪联合作用时,u为波浪水质点的速度矢量与海流速度矢量之和在垂直 于构件方向上的分矢量,根据海区的波浪参数选用不同的波浪理论，该平台波浪力用 司托克斯三阶波进行计算，并对不同入射方312 6期徐秀丽等:大型海洋平台基础结

25、构有限元分析214 自16卷然灾害学报向进行搜索，以得到在以某入射方向产生最大合波力及相应波峰位置。 根据计算的波浪荷载，在标高-27. 8 m （冲刷面以上对应梁单元以线荷载的形式施 加波流力，对应于实际情况下波浪对桩的作用力为东西向，在模型中为+ Y向（如图 2所示。2风荷载作用于平台上的风荷载按下式计算H:F = klk2pOA2（4式 中:kl,k2分别为风荷载体形系数及海上风压高度变化系数；pO为基本风压P0 = av ; a为风压系数；v, A分别为设计风速及受风面积。根据计算的风荷载（东西方 向，在标高+ 3. 7 m （海浪面以上至桩顶的相应单元上施加均布风荷载112kN/mo

26、 （3恒载及使用活荷载恒载为构件及建筑做法对应的重量，活荷载为由功能实施过 程中产生的荷载。（4温度荷载物体受温度荷载作用时，如果不受外界约束，则物 体只发生自由变形并无内部应力，只有当物体处于温度变化状态旦又有外界约束， 其变形受到限制时才会产生温度应力，故温度应力是由物体所受温度荷载作用与外 界约束共同作用的结果。本工程所施加荷载为：钢筋混凝土面板与钢梁体系升温 13. 3 C,体系降温21. 3 C。（5动力荷载波流力、地震力等，本文只考虑（1 （4 静荷载的作用，平台基础在其它荷载作用下的效应另行分析。4分析结果通过对 采用前述方案4建立的海洋平台基础有限元模型进行分析，得到平台基础在

27、恒载、 波流荷载及温度荷载作用下的内力（应力，见表2。根据分析结果，可对平台基础做 出如下评价：表2各荷载工况下海洋平台计算结果（方案4 Table 2 Numencal results of offshore p latfbnn under loaduig scenarios （ scheme 4 恒载 m in - 5 671. 13 - 33. 80 - 130. 36 - 917. 18 - 266. 79 - 44. 53 - 281. 98 - 456. 13 - 66. 05 -31. 17- 1. 09 - 0. 79 - 0. 14 - 2. 29 max 2 753. 13

28、 677. 92 648. 62 4 784. 54 4 983. 36 184. 29 1 220. 39 1 798. 10 269. 68 119. 88 4. 38 7. 46 9. 99 13.42 结构部位波流 min-2 610. 71 - 526. 54 -640. 97-4 727. 79 - 3 324. 22 - 184. 29 - 1 888. 13 - 1 731. 99 -254. 20- 119.81 -4. 79 - 7. 52-1. 14-13. 13 max 282. 33 267. 70 51. 57 372. 80 613. 12 37. 93 231.

29、45 413. 10 58. 79 18. 20 2. 08 3.31 3. 48 3. 48 温升 m m - 225. 81 - 155.52 - 27. 44 - 453. 84 - 874. 89 -37. 92 - 306. 88 - 324. 23 - 58. 79 - 19. 97-2. 14-2. 94- 1.32-5.31 max 361.64 249. 07 43. 94 726. 83 1 401. 13 60. 74 491. 46 519. 25 94. 15 31. 98 3.43 4. 71 8. 50 8.50 温 降 min-452. 15 -428. 72

30、- 82. 58 - 597. 04 - 981. 91 - 60. 74 - 370. 66 - 661. 59 - 94.15 -29. 15-3. 33 -5.30- 1.56-5. 57 内力 maxFx/kNFy/kNFz/kN237. 1429.65 桩 133. 56 721.98 266. 7944.53 My/(kN - mMz/(kN - inFy/kNFz/kN 305. 55 298. 06 66. 05 16. 58 1. 12 1. 05 1. 55 2. 20 梁 My/(kN mMz/(kN - m M zh / (N - nun S xx / (N - nun

31、 -2 -2 板 S yy / (N - nun 2 S pl / (N inm - 2 S pm ax / (N - nun -2:桩、注梁内力取值基于局部坐标系，沿杆件轴线方向为x轴，与其垂直的其余两轴分别为y轴、轴，M Zh为组合应力；混凝土板应力取值基Z 于全局坐标系O6期徐秀丽等：大型海洋平台基础结构有限元分析215(1波流荷载对平 台基础各构件(钢管桩、钢梁、混凝土板的影响是所有荷载中最大的，其次是温 降作用，恒载和风荷载影响相对较小。(2在波流荷载作用下，桩内会产生较大的 拉力；在自重荷载与施工荷载作用下，所有桩的轴力均为压力，但在桩基础的施工 过程中，由于上部结构尚未形成，恒载较小，需采取一定的压桩措施来防止桩被拔出 或波流冲倒。表3平台构件应力校核(3经对平台各构件在规范规定的各种组 合荷载 Table 3 Check of p latfoim stresses in components of p latfoi m 作用下的应 力进行验算，钢管混凝土桩与钢结构主梁均满足材料强度要求，见表3。(4由于本 文分析中未考虑混凝土板中预应力钢筋与普通钢筋的作用，表2中列出的混凝土板 的应力为名义应力,混凝土中的实际应力需另行分析。结构部位桩梁校核应力 /106 Pa 72. 6 164.