新型FRP桥墩防撞浮箱有限元仿真分析

桥墩防撞方面的优势，应用 FＲP 防撞浮箱能达到“既不伤墩、又少伤船、还少结
构自伤”的设计目的，实现长效防撞。
关键词： 纤维增强复合材料; 桥墩; 防撞浮箱; 船舶; 有限元仿真分析
中图分类号： U443． 22
文献标识码： A
重庆黄花园嘉陵江大桥主桥上部结构为 137． 16 m + 3 × 250． 00 m + 137． 16 m 的五跨预应力混凝土 连续刚构桥，其连续长度 1 024． 32 m。大桥于 2000 年建成交付使用后，大桥被船撞击事故时有发生。 为保证大桥安全，对船撞风险较大的 2 号主墩和 3 号主墩设计了围护在桥墩墩柱周围的纤维增强复合 材料（ FＲP） 防撞浮箱作为被动防撞保护结构，使桥墩受到不大于其设防船撞击力的失控船舶撞击时，经 FＲP 浮箱消能后传给桥墩的船撞击力小于桥墩的自身抗力，从而保证桥墩和全桥结构安全。
图 5 裸墩工况系统能量变化
图 6 带 FＲP 防撞箱工况系统能量变化
图 7 FＲP 防撞浮箱能量变化
图 8 裸墩工况桥墩受到的撞击力结果
第 23 卷第 1 期
杜柏松，等： 新型 FＲP 桥墩防撞浮箱有限元仿真分析
13
从能 量 角 度 分 析，对 于 裸 墩 而 言，船 舶 动 能 减 小而内能迅速增加，说明船舶变形极大； 而撞击 FＲP 浮箱时，船舶的内能增加较小，但 FＲP 浮箱的内能 增加极快，说明 FＲP 浮箱产生了较大的变形，吸收 了能量。桥 墩 两 端 约 束，无 法 大 幅 度 运 动，因 此 以 内能为主，并且在带有 FＲP 浮箱时，内能峰值出现 比裸墩时晚，船 直 接 撞 击 裸 墩 时，峰 值 出 现 时 间 没 有超过 0． 6 s，而采用浮箱保护后，撞击峰值出现的 时间都超过了 1． 6 s。从能量的角度可以看出，FＲP 浮箱作用是非常明显的，吸收了能量的同时，减小 图 9 带 FＲP 防撞浮箱工况桥墩受到的撞击力结果 了桥墩能量的峰值并且延长了桥墩能量峰值出现的时间。
表 4 FＲP 防撞浮箱的船撞动能转换及其消能比例计算结果
浮箱保护工 况船舶初始
动能 / MJ
7． 2
有浮箱保护工况船舶最大
撞击动能转换 / MJ
浮箱结
桥墩结
船舶保
构消能
构消能
留动能
3． 83
0． 01
0． 59
浮箱消能占 船舶最大撞 击动能比例 /%
53． 2
从桥墩受到船撞击力角度分析，桥墩受到的撞击力计算对比结果如表 5，可以看出带有 FＲP 浮箱撞 击工况下的桥墩应力相比裸墩受船撞击时的应力降低幅度极大，说明 FＲP 防撞浮箱保护桥墩效果 显著。
在裸墩工况下船撞整体有限元模型如图 3，船头直接与桥墩接触撞击，桥墩与船头刚度均较大，缓 冲距离小、撞击时间短、相互作用力大，对桥墩安全十分不利。
带 FＲP 防撞浮箱工况的船撞整体有限元模型如图 4，FＲP 防撞浮箱为桥墩提供了缓冲，吸收了船的 动能，延长了接触时间，使桥墩受力减少。
图 3 中心正撞裸墩头部
图 10 裸墩撞击工况船首变形
图 11 带 FＲP 防撞浮箱撞击工况 船首及浮箱变形图
14
河南城建学院学报
2014 年 1 月
5 结论 桥梁防撞结构经过多年发展，防护系统的种类已经很多，能够给予桥梁长期的稳定、安全、易于维护
的防撞保护逐渐成为可实现的目标，FＲP 防撞浮箱充分发挥新材料的特性，加以合理结构设计，在撞击 吸能方面有突出优势，结构的高度非线性使撞击时间大大延长，撞击过程中对船和桥墩应力时程曲线趋 于缓和，既保护了桥梁又保护了船舶。同时，浮箱结构的模块化组件，单构件质量轻，在发生撞击后部分 组件如发生破损可快速方便的更换组件，恢复整体结构可用性。本次计算中 FＲP 的材料优势明显，为 结构产生大变形非线性变形后产生尽可能少的结构自伤只有较好效果，FＲP 在其他结构类型防撞装置 中的应用也应进行探索尝试。
根据设计要求，计算采用内河航运常用的驳船模型，总质量 1 600 t，船舶材料参数如表 1、表 2 所 示，船舶采用壳单元划分，为了减小计算量，船头单元尺寸较小，平均尺寸 10 cm，而船尾在撞击过程中 变形极小，可以采用刚性单元，最 大 单 元 尺 寸 取 115 cm，全 船 共 划 分 单 元 20 262 个，有 限 元 模 型 如 图 1。
Vol． 23 No． 1 Jan． 2014
新型 FＲP 桥墩防撞浮箱有限元仿真分析
杜柏松，李亚迪，张锡祥，姜 旭
( 重庆交通大学土木建筑学，重庆 400074)
摘 要： 随着航道通航量增加，对桥墩的防护不断有新方法引入，基于在桥墩防撞结构
中首次引入 FＲP 新型材料的防撞浮箱有限元仿真分析，探索这种新型结构在
表 5 桥墩的船撞击力计算对比结果
桥墩承受的最大船撞击力 Pz / MN
无浮箱保护
有浮箱保护
模式 Pz1
模式 Pz2
39． 6
7． 5
桥墩承受的 船撞击力
降低比例 /% 81． 1
从船头受到的应力角度分析如图 10、图 11，采用 FＲP 浮箱保护后，船头的应力有大幅度降低； 船舶 撞击凹陷深度会有较大幅度的降低，撞击并未对船舶整体产生破坏性作用，仅局部产生变形，表明 FＲP 浮箱能够达到“少伤船”的目的，说明 FＲP 防撞浮箱保护船舶效果显著。FＲP 防撞浮箱的变形能力也使 得自身在发生变形后有较少损伤，损伤组件可方便更换，整体结构可继续使用。综上所述，就数值计算 结果而言，FＲP 防撞浮箱能达到“既不伤墩、又少伤船、还少结构自伤”的设计目的。
泊松比 0． 20
2． 4 单元退出计算的准则 由于 FＲP 材料的防撞浮箱材料为脆性材料，第一强度理论适用于这种脆性材料，该理论认为最大
拉应力是使材料发生断裂破坏的主要因素，即认为不论是什么应力状态，只要最大拉应力达到与材料性
12
河南城建学院学报
2014 年 1 月
质有关的某一极限值，材料就发生断裂。既然最大拉应力的极限值与应力状态无关，于是就可用单向应 力状态确定这一极限值。单向拉伸时只有 σ1 （ σ2 = σ3 = 0） ，当 σ1 达到强度极限 σb 时即发生断裂。故 据此理论得知，不论是什么应力状态，只要最大拉应力 σ1 达到 σb 就导致断裂。在模型试验中，FＲP 材 料的拉伸强度已经得到，横向为 300 MPa，纵向为 150 MPa。按弹性材料计算，碰撞完成后，应力超过 150MPa 的单元假设为破损单元，该方法结果为近似结果，计算速度较快。 3 碰撞及接触条件的定义
DOI:10.14140/ki.hncjxb.2014.01.013
第 23 卷第 1 期 2014 年 1 月
河南城建学院学报 Journal of Henan University of Urban Construction
文章编号： 1674 － 7046（ 2014） 01 － 0010 － 05
型参见图 2。
2． 3 桥墩模型的建立
本次计算选取受力不利的 2 号墩进行，墩和主梁为
固结体系，主梁质量相当大，在船 － 墩的碰撞过程中，碰
撞时间非常短，主梁提供了强大的惯性约束，因此，主梁
顶端可以简 化 为 固 结 模 式，墩 底 为 基 岩，亦 可 采 用 固 结
模式。
桥墩采用六面体单元划分，主要单元尺寸为竖向 50
一年当中嘉陵江水位是不断变化的，本文选取最低水位黄海高程 174 m，具有代表作用的工况进行 计算分析。船按照 1． 5 m 吃水深度取值，撞击吨位和撞击速度分别选择设计时考虑的最大撞击吨位 1 600 t 和速度 3 m / s。 1 有限元仿真分析目的
船桥碰撞是一个复杂的非线性动力学问题，是通过碰撞区域内的结构或者构件的变形和压溃，在短 时间内将巨大的冲撞能量吸收、转化和耗散的非线性动力响应过程。应用有限元仿真分析研究 FＲP 应 用在桥梁防撞结构中的耗能效果，通过有限元仿真计算为实际工程试验方案完善和试验参数选取提供 依据，更期望通过试验与有限元仿真计算结果对比，检验针对同一结构和相同荷载工况及边界条件进行 的结构计算的可靠性及计算结果可信度，从而为黄花园大桥的桥墩 FＲP 防撞浮箱工程施工图设计提供 既有科学依据、又有足够精度还偏于安全的计算方法。 2 有限元模型的建立 2． 1 代表船型选取和船舶有限元模型的建立
参考文献
［1］ 张锡祥，王智祥． 一种新型 FＲP 桥墩防撞浮箱结构［J］． 重庆交通大学学报： 自然科学版，2011，30（ 3） ： 388 － 393． ［2］ 肖波，鹿道凡． 船桥碰撞问题的有限元仿真分析［J］ ． 船海工程，2011，40（ 1） ： 143 － 146． ［3］ 陈国虞． 防御船撞桥的新装置及其机理研究［J］． 船舶工程，2007，29（ 4） ： 40 － 43． ［4］ 肖波，王爽． 桥墩防撞装置碰撞动力学分析［J］． 武汉理工大学学报： 交通科学与工版，2005，29（ 1） ： 49 － 51，67． ［5］ 汪宏，耿波． 菜园坝长江大桥船撞风险分析咨询报告［D］． 重庆： 重庆交通科研设计院，2008．
件。各构件的材料选择和截面尺寸，根据构件的抗撞受力要求和能量吸收、动量缓冲要求按复合材料力
学理论分析计算和设计确定。
FＲP 防撞浮箱采用壳单元划分，单元主要控制尺寸为 20 cm，为了提高计算速度，尽量采用四边形
单元进行划分，整个浮箱约划分 160 000 单元，整个系统共约划分 26． 4 万个单元，整个系统的有限元模
图 4 中心正撞 FＲP 防撞浮箱头部
碰撞过程中，碰撞系统的能量是守恒的，碰撞过程也即是能量的转换过程，从能量的状态变化可以 看出计算可靠性。在船舶撞击裸墩的过程中，动能向内能转换是主要出现的能量交换方式。裸墩及带 浮箱船撞计算的能量变化状态参见图 5 ～ 图 7，桥墩受撞击力变化参见图 8 ～ 图 9。
cm 长，截面尺寸为 20 cm，在船舶撞击部位，竖向尺寸加
密，碰撞区最大竖向尺寸为 20 cm，共约划分 84 000 个单
元。材料参数如表 3。
图 2 FＲP 防撞浮箱有限元模型
项目 桥墩
表 3 桥墩模型材料参数
材料特性 线弹性材料
容重 / ( kg / m3 )
2 600
弹性模量 / Pa
3． 45E + 10
容重 / ( kg / m3 )
7 800 7 800
弹性模量 / Pa
2． 10E + 11 2． 10E + 11
泊松比
0． 3 0． 3
项目 船头
屈服应力 / Pa
2． 35E + 08
表 2 船头补充参数
切线模量 / Pa
1． 18E + 09
应变率 参数 C 40． 4
应变率 参数 P
5
失效 应变 0． 34
2． 2 FＲP 防撞浮箱模型的建立
FＲP 桥墩防撞浮箱的总体结构形式，为围护在桥墩
四周的 FＲP 外围箱体结构与附着于外围箱体结构内表
面上的一排 FＲP 内衬八边形柱壳薄壁构件组成的组合
结构。FＲP 外围箱体结构，为浮箱的抗撞消能主体结构
和浮力平衡结构； FＲP 内衬柱壳构件，为浮箱与桥墩的弱
接触连接支承结构和撞击缓冲结构。
收稿日期： 2013 － 12 － 20 第一作者简介： 杜柏松( 1976 － ) ，男，湖北黄冈人，博士，重庆交通大学土木建筑学院副教授。
第 23 卷第 1 期
杜柏松，等： 新型 FＲP 桥墩防撞浮箱有限元仿真分析
11
项目
船头 船尾
表 1 桥墩及船舶主要构件的材料参数表
材料特性
弹塑性材料 刚性材料
图 1 撞击船舶有限元模型
黄花园嘉陵江大桥 FＲP 防撞浮箱材料及构件基本力学性能根据拉伸试验、压缩试验、弯曲试验、剪
切试验、冲击韧性试验等得出。FＲP 桥墩防撞浮箱新结构，为全 FＲP 结构，各构件分别为碳纤维增强复
合材料（ CFＲP） 或玻璃纤维增强复合材料（ GFＲP） 或芳纶纤维增强复合材料（ AFＲP） 杆系构件或薄壳构
接触是碰撞分析常用的一种约束方式，整个浮箱由 12 个子箱组成，子箱之间通过卡榫连接。从力 学模型上分析，是一种典型的接触受力模式。本计算模型中，船和浮箱或桥墩之间 LS-DYNA 中提供的 面—面接触； 浮箱与桥墩之间也采用面—面接触； 整个浮箱采用计算速度最快的自动单面接触模式。接 触参数： FＲP 材料之间的接触静摩擦系数取 0． 233，动摩擦系数取 0． 15； 船舶静摩擦系数取0． 30，动摩擦 系数取 0． 20； 船舶对 FＲP 箱或桥墩混凝土的静摩擦系数取 0． 30，动摩擦系数取 0． 20； FＲP 材料与混凝 土之间的静摩擦系数取 0． 30，动摩擦系数取 0． 20。荷载以速度的方式体现，在船舶行驶方向设定计算 速度 3 m / s，不计加速度及自重。计算时间取 3 s，试算发现，对于裸墩而言，计算时间超过 1s 即可； 对于 带 FＲP 防撞浮箱的结构，则采用 3 s。 4 有限元仿真结果及分析