带主梁的简化模型与响应面联合的桥梁船撞易损性分析方法

樊伟　孙洋　申东杰　刘斌

摘   要：以往桥梁船撞动力分析以给定事件确定性分析为主，难以反映船撞作用的偶然性和概率性特征，以及不同能量撞击下的桥梁损伤演化特征. 为此，本文面向两类典型船舶，以桥梁墩柱受船舶撞击后的剩余承载能力作为损伤评估指标，较为系统地研究了桥梁船撞易损性. 首先，建立了受压RC墩柱受到侧向冲击后剩余承载能力的直接模拟方法，通过与试验结果进行对比，验证了该模拟方法的有效性. 然后，基于一座典型连续钢筋混凝土梁桥，建立了两种不同的有限元简化模型，并进行了比较和验证. 提出一种有限元简化模型与响应面代理模型联合的桥梁船撞易损性分析方法，獲得了两类典型船舶撞击下的桥梁易损性曲线. 结果表明：所建立的响应面具有良好的精度，可替代复杂的非线性有限元计算;两类船舶类型撞击下的桥墩剩余承载力的响应特征区别较大，在球艏船撞击下剩余承载能力随船速的增大而均匀减少，而在受驳船撞击时，剩余承载能力与临界船速密切相关，呈现出双折线的特征，在进行样本设计时需基于临界速度进行分段;在相同船速及质量的情况下，驳船撞击所造成的结构损伤以及失效的概率普遍要高于球艏船撞击，实际设计中应尤为关注.

关键词：船撞;易损性;剩余承载能力;响应面;简化模型

中图分类号：U447                                    文献标志码：A

Vessel-collision Vulnerability Analysis Method of Bridge Structures

Based on Simplified Model with Girders and Response Surface

FAN Wei1，2?，SUN Yang1，SHEN Dongjie1，LIU Bin1

（1. College of Civil Engineering，Hunan University，Changsha 410082，China;

2. Hunan Provincial Key Laboratory for Wind and Bridge

Engineering，Hunan University，Changsha 410082，China）

Abstract：The past studies primarily performed the deterministic dynamic analysis of vessel-bridge collision for a given collision event，which hardly reflected the contingency and probability characteristics of ship collision，as well as the damage evolution of bridge under impacts with various energy. For this reason，this study systematically investigated the vulnerability of a bridge under collisions of two typical vessels，where the residual capacity of bridge piers after vessel impacts was used as the damage assessment index. Firstly，a direct finite element （FE） simulation method was established to predict the residual capacity of axially-loaded RC column after lateral impact. The rationality of the simulation method was verified by the test results. Then，based on a typical continuous girder bridge，two different simplified FE model were established，compared and validated. By combining the validated simplified FE model and the response surface method，a vulnerability analysis method of the vessel-bridge collision was established，and the vulnerability curves of the bridge were obtained for the impacts of two kinds of typical ships. The results show that the response surface has good accuracy and is able to replace the complex nonlinear FE calculation;The response characteristics of the residual capacity of the pier are quite different for the impact of two types of vessels：the residual capacity of the pier under the impact of the bulbous-bow ship decreases uniformly with the increase of the ship speed，while the residual capacity is closely related to the critical ship speed under the impact of the barge and exhibits the bilinear characteristics so that the sample design needs to be segmented based on the critical speed;Under the same ship speed and quality，the probabilities of structural damage and failure caused by barge impacts are generally higher than those caused by bulbous-bow ship impacts，which should be paid special attention in the practical design.

Key words：vessel-bridge collision;vulnerability;residual load capacity;response surface;simplified model

船撞桥事故不但威胁着船舶通行的安全，也严重影响着桥梁的安全运营[1]，一旦发生船撞桥事故，桥梁结构可能承受巨大的侧向冲击荷载，因此在对通航水域的桥梁进行设计时，必须考虑船桥发生碰撞的可能性.

20世纪60年代末，人们开始研究船桥碰撞问题，著名的米诺尔斯基（Minorsky）理论[2]就在那时被提出，该理论为后人研究船桥碰撞奠定了基础. 目前，桥梁船撞安全问题已经受到了广泛关注. 相关研究主要集中于撞击力估算[3]、船撞桥数值模拟[4-6]、桥梁防撞设施研究[7-10]等. 但是，以往的这些研究大多都是基于确定性的有限元仿真计算进行分析，难以反映船撞作用的偶然性和概率性特征，以及不同能量撞击下的桥梁损伤演化特征. 对桥梁进行船撞易损性分析可以预测结构在不同类型的船舶以及不同船速撞击下发生各级破坏的概率，对结构的设计、加固和维修决策具有实际工程应用价值.

目前，桥梁的地震易损性研究已经受到国内外学者的广泛关注[11-15]. Singhal等[16]采用贝叶斯原理，分析1994年1月17日北岭加州地震混凝土框架建筑的地震损伤数据，采用Park-Ang地震损伤指数表示结构损伤，将地震易损性定义为在给定的地震强度下，损伤指标超过某一定值的条件概率. 目前，地震易损性的理论研究已经比较丰富.

与地震易损性研究相比，关于桥梁船撞易损性的研究成果却寥寥无几. 首先，船撞领域尚未有人提出一个广泛令人接受和信服的判断桥墩船撞损伤等级的指标. 张太磊[17]通过数值仿真计算，得出以墩底转角作为评价桥梁损伤等级的指标，并指出混凝土桥墩在受到船舶撞击后的破坏形式与地震作用下的破坏形式有所不同. 但是其仅仅对矩形墩柱在船舶正向撞击方面进行了研究，因此其损伤指标的适用性还有待进一步研究.

此外，前期研究[18]表明，桥梁船撞过程动力效应影响明显，对船撞下桥梁结构的影响进行合理的动力分析是必要的，因此，进行桥梁结构的易损性分析时，就需要大量样本的动力计算. 如果采用常规的精细化接触-碰撞有限元技术进行船桥碰撞的非线性显示动力分析，单个模型的计算将会消耗大量的计算时间，计算效率低. 若要进行大量样本计算，精细化全尺寸有限元模型必然会制约运算规模，因此，提高计算效率非常关键. 近年来也有学者[19-20]将响应面运用到了撞击下桥梁参数分析和可靠度分析中，极大地提高了计算效率. 其中就有张军等[19]提出运用响应面法进行船撞桥的可靠度计算，但是其并未进行非线性有限元动力碰撞计算，而是仅仅以规范的计算方法确定船撞力.

由此可见，高精度高效率的有限元简化模型与响应面代理模型的联合是实现大样本分析的前提条件. 因此，本文提出有限元简化模型与响应面代理模型联合的桥梁船撞易损性分析方法，为构建基于概率性的桥梁船撞设计与评估方法奠定基础.

1   墩柱剩余承载能力试验与直接模拟方法

1.1   墩柱剩余承载能力试验

Fan和Liu等[21-22]对一组不同参数的受压RC墩柱进行了落锤冲击试验，并对受到落锤冲击后的受压RC墩柱进行了轴向加载试验. 试验共对10根不同参数的RC墩柱进行了轴向加载试验，其中有2根未受损伤的试件，还有8根受落锤冲击后的受损试件. 这8根受损试件中，有1根被完全砸断，还有2根出现了“反拱”现象[21]. Fan等[21]指出，反拱的现象具有偶然性，存在一定的不可重复性. 因此，在本次模型验证中排除了这2根“反拱”试件以及1根完全砸断的试件，对剩余的7根试件开展了数值模拟.

1.2   墩柱剩余承载能力直接模拟方法

为了研究墩柱受冲击后剩余承载能力的直接模拟方法，采用LS-DYNA显式非线性动力有限元软件建立了如图1所示的精细化有限元模型.

其中，纵向和螺旋钢筋均采用 Hughes-Liu 积分梁单元模拟，采用*MAT_PIECEWISE_LINEAR_ PLASTICITY材料本构，该材料本构可以采用试验测得的实际应力-应变关系来定义. 混凝土采用六面体单点積分实体单元模拟，材料模型采用连续盖帽本构模型*MAT_CSCM，该材料本构已被证明可以合理地模拟冲击荷载下受压RC柱的动力响应[22-23]. 钢筋的梁单元与混凝土的实体单元之间使用了*MAT_GENERAL_NONLINEAR_1DOF_DISCRERE_ BEAM非线性弹簧单元相连接，来模拟钢筋与混凝土之间的粘结滑移关系. 非线性弹簧单元的加载卸载曲线的定义参照欧洲规范[24]中的规定.

为了模拟出落锤冲击后进行轴向加载的试验全过程，在建模过程中将三个试验阶段合并在一个计算模型中，如图1所示. 整个计算过程分为三个阶段，计算时间为0.14 s.

第一阶段为初始轴力的预加载，时间为0 s～0.022 s. 模型中采用关键字*CONSTRAINED_ NODAL_RIGID_BODY将预应力钢筋的两端与试件两端的钢板固定起来，预应力钢筋由索单元模拟，通过定义材料*MAT_CABLE_DISCRETE_BEAM将初始应力赋给了预应力钢筋;第二阶段为落锤冲击阶段，时间为0.022 s ～ 0.08 s. 在初始轴力完成加载后，落锤将以给定的初速度下落并与受压RC柱发生碰撞，在早于0.08 s的某个时间点完成全部碰撞过程，撞击力归零，锤头与试件完全分离;第三阶段为轴向加载阶段，时间为0.08 s ～ 0.14 s. 通过关键字*BOUNDARY_PRESCRIBED_MOTION_SET控制墩柱的两端钢板，在0.08 s时采用预先设定好的轴向位移开始压缩墩柱，即位移控制的轴向加载. 加载过程持续到墩柱失去轴向承载能力，结构失效.

0.2 ～ 0.5   （中度损伤）

0.5 ～ 0.8   （严重损伤）

0.8 ～ 1.0   （完全倒塌）       （4）

式中：Pr为桥墩受到损伤后的残余轴向承载能力;Pd为桥墩的初始轴向承载能力.

尽管采用简化分析模型能够提高计算效率，但需要大量计算的易损性分析而言，仍然是不够的. 因此，本文将用简化分析模型获得关键样本结果，采用响应面方法建立合理的响应面替代模型. 在此基础上，结合蒙特卡洛抽样的方法计算获得不同船撞情况下各等级损伤下的失效概率.

3.2   试验设计

Box-Behnken试验设计法（Box-Behnken Design，简称BBD），是响应面优化法常用的试验设计方法，以三因素为例，试验设计的试验点分布情况如图8所示.

合理的试验方法是获得良好的响应面代理模型的前提，因此，在进行响应面工况设计之前，对桥墩在分别受到驳船和球艏船不同船速撞击下的结构响应进行了数值模拟，结果如图9所示.

可以看出，桥墩在两种不同类型船撞击下的响应特征截然不同. 在球艏船的撞击下，桥墩的响应随着船速的增大近似线性上升，而在驳船的撞击下，桥墩的各响应呈现出双折线的形式. 当驳船的船速在0～0.8 m/s内上升时，桥墩的峰值响应迅速上升;当船速超过0.8 m/s后，峰值响应基本上不会再有很大变化. 产生这种现象的原因可能是，当撞击速度超过0.8 m/s后，驳船的船首在发生船桥碰撞的过程中前肋板屈曲，船首刚度突然下降至低于桥墩的侧向抗推刚度，驳船剩余的动能几乎全部由船首吸收，因此，桥墩响应未见上升而驳船船首变形严重. 这也是当船速超过0.8 m/s后，撞击力-时程曲线出现平台段的原因，如图10所示.

在上述工况分析中，可以认为0.8 m/s是使船首屈曲的临界船速v0. 根据试算可知，临界船速的大小并不是一个固定的数值，它与船舶总质量、混凝土强度、钢筋直径等因素息息相关.

考虑到桥墩在受到驳船和球艏船撞击下的不同响应特征，如果继续使用传统的Box-Behnken设计法，仅仅在每种因素（混凝土强度、钢筋直径、船速）的最大值、最小值和中心点进行试验，那么船速这一因素对驳船撞击下桥墩各响应的影响将会被错误地估计，如图11所示. 由这样的试验数据点拟合出来的响应面代理模型并不能反映桥墩受驳船撞击后响应的真实特征，代理模型外推能力较差.

对于船速这一因素来说，临界速度v0是一个重要的数据点. 在进行工况设计时，必须引入临界速度的影响. 因此，本文针对驳船撞击下的工况，提出了一种基于临界速度的分段BBD试验设计法，将原本[0，4]的船速区间分为[0，v0]和[v0，4] 两段. 对每一个区间分别进行BBD试验设计，如图11所示.

当桥梁受到船撞时，响应结果往往受到各种参数的随机性的影响，影响桥墩受冲击后剩余承载能力的主要参数有混凝土强度fc、纵向钢筋直径dl、箍筋直径ds、纵向钢筋屈服强度fy、箍筋屈服强度fys和船速v等. 根据试算可知，ds、 fy、 fys在其服从的随机分布范围内的变化对剩余承载能力的影响较小，故在本次研究中并未考虑这3个参数的随机性带来的影响，并在响应面代理模型的计算中将这3个参数设为固定值（ds = 16 mm，fy = 400 MPa，fys = 335 MPa）. 为了对桥梁船撞易损性进行分析，本节将建立如下响应面代理模型来预测桥梁受船撞后的残余承载能力：

Pr = f（ fc，dl，v）       （5）

参考现有文献的结论[28-30]，fc和dl均服从正态分布，统计参数如表1所示.

参考两种参数的概率密度分布情况，在进行试验设计时，将钢筋直径的变化范围设为28 ～ 36 mm，将混凝土强度的变化范围设为20 ～ 40 MPa. 设计工况如表2所示.

3.3   建立响应面代理模型

3.3.1   球艏船撞击下的桥墩响应代理模型

多项式响应面代理模型是响应面分析中常用的代理模型形式，本文首先尝试采用多项式代理模型对桥墩受驳船以及球艏船撞击下的剩余承载能力进行拟合.

代理模型采用多项式形式，其函数形式为：

Pr = p1 + p2 fc + p3 dl + p4 v + p5 fc dl + p6 fc v +

p7 dl v + p8 fc 2 + p9 dl 2 + p10 v2 + p11 fc 2dl +

p12 fc dl2 + p13 fc v2         （6）

式中：fc为混凝土强度;dl为钢筋直径;v为船速;p1～p13为多项式中各分项的待定系数.

根据试验结果拟合得到桥墩受球艏船撞击后的残余承载力响应面代理模型多项式的各项系数如表3所示. 二次多项式代理模型即可达到令人满意的拟合精度，可决系数（R2）为0.98. 各个参数对船撞后的残余承载能力的影响如圖12所示，图中坐标轴上未显示的参数均取中值：v = 2 m/s，fc = 30 MPa，dl = 32 mm.

然而，在拟合桥墩受驳船撞击下的残余承载能力响应面时，多项式代理模型始终无法达到一个令人满意的拟合精度，无论是二次多项式还是三次或者更高次多项式，最后的效果始终差强人意.

3.3.2   驳船撞击下的桥墩响应代理模型

为了提高响应面代理模型的拟合精度，针对桥墩在受到驳船撞击下的残余承载能力，本文提出了一种基于临界船速的多项式分段函数的代理模型，如式（7）所示：

Pr-barg e = f1（ fc，dl，v）       v∈（0，v0）

f2（ fc，dl，v）       v∈（v0，4）      （7）

根据试验结果拟合得到桥墩受驳船撞击后的残余承载力响应面代理模型多项式的各分项系数如表3所示. 分段拟合后，三次多项式代理模型可以达到令人满意的拟合精度，分段函数f1和f2的可决系数（R2）分别为0.999 9，0.999 3.

桥墩在受驳船撞击后的残余承载能力与各参数之间的关系如图13所示. 当船速小于0.8 m/s或大于0.8 m/s时，桥墩残余承载能力的响应特征有着明显差异.

由图13（c）中船速和混凝土强度对桥墩残余承载能力的影响可明显看出，当混凝土强度下降时，临界船速也会变小，具体表现为随着混凝土强度降低，响应面的转折点向船速低的方向移动.

值得一提的是驳船撞击下的桥墩剩余承载能力，如图13（b）（c）所示，当船速超过临界船速并继续提高至约2.5 m/s时，桥墩的剩余承载能力继续平缓下降;但是，当船速超过2.5 m/s时，桥墩的剩余承载能力反而呈现出上升的趋势，这一现象与荷载频谱特性、桥梁结构动力特性相关[25].

3.4   易损性分析

根据墩柱残余承载力的损伤指标，即式（4），进行桥墩船撞的易损性分析. 其中Pr可由相应的响应面代理模型计算得到. 将fc = 30 MPa、dl = 32 mm、v =0 m/s代入相应的响应面代理模型，则可得到Pd =

f（30，32，0）.

对于每一个确定的船速v∈（0，4）：对服从正态分布的随机变量fc和dl进行蒙特卡洛抽样，样本数量N = 107，统计出损伤指标D分别达到轻微损伤、中等损伤、严重损伤和倒塌的次数L1、L2、L3、L4 . 则在该船速下，桥墩达到相应损伤指标的失效概率为：

Pfn =        （8）

式中：n = 1、2、3、4分别表示损伤程度达到轻微损伤、中等损伤、严重损伤和倒塌. 对船速v从0 m/s到4 m/s，船速每间隔0.02 m/s计算一次结构的失效概率，将得到的数据点连线则可得到桥墩受船撞的易损性曲线，如图14所示.

对比驳船与球艏船的易损性曲线，可知在相同船速及质量的情况下，驳船撞击所引起结构失效的概率普遍比球艏船要大很多. 相比较于球艏船撞击工况，驳船撞击下的桥墩在较低的船速撞击下就能造成较大的损伤.

4   结论与展望

本文面向两类典型船舶，以桥梁墩柱受船舶撞击后的剩余承载能力作为损伤评估指标，提出了基于带主梁简化分析模型与响应面联合的桥梁船撞易损性分析方法，主要结论有：

1）所建立的高精度响应面模型可以替代需要进行复杂非线性计算的结构模型. 可以使用所建立的响应面代理模型在橋梁船撞易损性研究中进行大量的样本分析.

2）由于船首构造、外形等差异，不同船舶类型撞击下的桥墩剩余承载力的响应特征区别较大. 在球艏船撞击下，桥墩的剩余承载能力会一直随船速的增大而均匀减少;而在受驳船撞击时，桥墩的剩余承载能力与临界船速密切相关，呈现出双折线的特征. 因此，在进行样本设计时，为了反映桥墩受驳船撞击后响应的真实特征，需基于临界速度进行分段.

3）在相同船速及质量的情况下，驳船撞击所造成的结构损伤以及失效的概率普遍要高于球艏船撞击. 在进行实际风险评估时，应该尤为关注.

4）在本文算例分析中，无论是驳船还是球艏船都未能对桥墩造成严重及以上的损伤状态，但是实际工程中却屡见船桥碰撞的恶性事故，这是因为实际工程中的桥梁船撞事故往往还伴随着钢筋锈蚀以及下部结构冲刷等多种灾害同时作用. 此外桩土、水位变化、撞击角度、通航等级与船只吨位等各种不确定因素都会对事故的结果造成影响，这些都是值得进一步研究的问题. 在船撞易损性分析中考虑钢筋锈蚀以及冲刷等因素的影响将是本文的下一步研究方向.

参考文献

[1]    项海帆，范立础，王君杰. 船撞桥设计理论的现状与需进一步研究的问题[J]. 同济大学学报（自然科学版），2002，30（4）：386—392.

XIANG H F，FAN L C，WANG J J. State of art of ship collision design for bridges and future research[J].  Journal of Tongji University（Natural Science），2002，30（4）：386—392. （In Chinese）

[2]    MINORSKY V U. An analysis of ship collisions with reference to protection of nuclear power plants[R]. New York：Sharp （George G） Inc，1958.

[3]    ZHANG J F，CHEN X Z，LIU D J，et al. Analysis of bridge response to barge collision：Refined impact force models and some new insights[J]. Advances in Structural Engineering，2016，19（8）：1224—1244.

[4]   CONSOLAZIO G R，COWAN D R. Nonlinear analysis of barge crush behavior and its relationship to impact resistant bridge design[J]. Computers & Structures，2003，81（8/11）：547—557.

[5]    姜华，贺拴海，王君杰. 混凝土弹塑性损伤帽盖模型参数确定研究[J]. 振动与冲击，2012，31（15）：132—139.

JIANG H，HE S H，WANG J J. Parameters determination of elasto-plastic damage cap model for concrete materials[J].  Journal of Vibration and Shock，2012，31（15）：132—139. （In Chinese）

[6]    劉建成，顾永宁. 基于整船整桥模型的船桥碰撞数值仿真[J]. 工程力学，2003，20（5）：155—162.

LIU J C，GU Y N. Simulation of ship-bridge head-on collision based on finite element model of whole ship-bridge[J].  Engineering Mechanics，2003，20（5）：155—162. （In Chinese）

[7]    FAN W，GUO W，SUN Y，et al. Experimental and numerical investigations of a novel steel-UHPFRC composite fender for bridge protection in vessel collisions[J]. Ocean Engineering，2018，165：1—21.

[8]    FANG H，MAO Y F，LIU W Q，et al. Manufacturing and evaluation of large-scale composite bumper system for bridge pier protection against ship collision[J]. Composite Structures，2016，158：187—198.

[9]    许薛军，单成林. 夹层结构曲面环形浮式桥墩防撞套箱碰撞分析[J]. 湖南大学学报（自然科学版），2015，42（3）：106—111.

XU X J，SHAN C L. Impact analysis of the bridge pier anti-collision floating box sets made by sandwich structure with curved-shaped[J]. Journal of Hunan University（Natural Sciences），2015，42（3）：106—111. （In Chinese）

[10]  FAN W，YUAN W C，CHEN B S. Steel fender limitations and improvements for bridge protection in ship collisions[J]. Journal of Bridge Engineering，2015，20（12）：06015004.

[11]  赵国藩. 结构可靠度的实用分析方法[J]. 建筑结构学报，1984，5（3）：1—10.

ZHAO G F.A practical method for structural reliability analysis[J].  Journal of Building Structures，1984，5（3）：1—10. （In Chinese）

[12]  李立峰，吴文朋，黄佳梅，等. 板式橡胶支座地震易损性分析[J]. 湖南大学学报（自然科学版），2011，38（11）：1—6.

LI L F，WU W P，HUANG J M，et al. Research on the sesmic vulnerability analysis of laminated rubber bearing[J]. Journal of Hunan University（Natural Sciences），2011，38（11）：1—6. （In Chinese）

[13]  吕大刚，李晓鹏，王光远. 基于可靠度和性能的结构整体地震易损性分析[J]. 自然灾害学报，2006，15（2）：107—114.

L? D G，LI X P，WANG G Y. Global seismic fragility analysis of structures based on reliability and performance[J].  Journal of Natural Disasters，2006，15（2）：107—114. （In Chinese）

[14]  张菊辉.基于数值模拟的规则梁桥墩柱的地震易损性分析[D]. 上海：同济大学，2006：31—78.

ZHANG J H.Study on seismic vulnerability analysis of normal beam bridge piers based on numerical simulation[D].  Shanghai：Tongji University，2006：31—78. （In Chinese）

[15]  刘晶波，刘阳冰，闫秋实，等. 基于性能的方钢管混凝土框架结构地震易损性分析[J]. 土木工程学报，2010，43（2）：39—47.

LIU J B，LIU Y B，YAN Q S，et al.Performance based seismic fragility analysis of cfst frame structures[J].  China Civil Engineering Journal，2010，43（2）：39—47. （In Chinese）

[16]  SINGHAL A，KIREMIDJIAN A S. Bayesian updating of fragilities with application to RC frames[J]. Journal of Structural Engineering，1998，124（8）：922—929.

[17]  张太磊. 混凝土桥墩船撞损伤指标研究[D]. 重庆：重庆交通大学，2015：63—89.

ZHANG T L. Study on the concrete pier in ship knock injury index aspect[D]. Chongqing：Chongqing Jiaotong University，2015：63—89. （In Chinese）

[18]  FAN W，LIU Y Z，LIU B，et al. Dynamic ship-impact load on bridge structures emphasizing shock spectrum approximation[J]. Journal of Bridge Engineering，2016，21（10）：04016057.

[19]  张军，邓旭平，包龙生. 响应面法在船撞桥可靠度分析中的应用[J]. 石家庄铁道大学学报（自然科学版），2015，28（1）：28—33.

ZHANG J，DENG X P，BAO L S. Application of response surface method in reliability analysis of ship collision[J].  Journal of Shijiazhuang Tiedao University（Natural Science），2015，28（1）：28—33. （In Chinese）

[20]  FAN W，XU X，ZHANG Z Y，et al. Performance and sensitivity analysis of UHPFRC-strengthened bridge columns subjected to vehicle collisions[J]. Engineering Structures，2018，173：251—268.

[21]  FAN W，LIU B，CONSOLAZIO G R. Residual capacity of axially loaded circular RC columns after lateral low-velocity impact[J]. Journal of Structural Engineering，2019，145（6）：04019039.

[22]  LIU B，FAN W，GUO W，et al. Experimental investigation and improved FE modeling of axially-loaded circular RC columns under lateral impact loading[J]. Engineering Structures，2017，152：619—642.

[23]  FAN W，YUAN W C. Numerical simulation and analytical modeling of pile-supported structures subjected to ship collisions including soil-structure interaction[J]. Ocean Engineering，2014，91：11—27.

[24]  Model Code 2010-First Complete Draft （Volumn 1）： Fib Bulletin 55[S]. Lausanne，Switzerland：Fédération Internationale du Béton，2010：232—238.

[25]  FAN W，SHEN D J，HUANG X，et al. Reinforced concrete bridge structures under barge impacts：FE modeling，dynamic behaviors，and UHPFRC-based strengthening [J]. Ocean Engineering，2020，21（6）：108—116.

[26]  董銳.圆柱形桥墩在冲击荷载作用下的损伤指标试验研究[D]. 重庆：重庆交通大学，2016：11—13.

DONG R. Experimental study on damage mechanism of cylindrical bridge pier under impact load[D].  Chongqing：Chongqing Jiaotong University，2016：11—13. （In Chinese）

[27]  SHI Y C，HAO H，LI Z X.Numerical derivation of pressure-impulse diagrams for prediction of RC column damage to blast loads[J]. International Journal of Impact Engineering，2008，35（11）：1213—1227.

[28]  MIRZA S A，MACGREGOR J G. Variability of mechanical properties of reinforcing bars[J]. Journal of the Structural Division，1979，105（5）：921—937.

[29]  MIRZA S A，MACGREGOR J G，HATZINIKOLAS M. Statistical descriptions of strength of concrete[J]. Journal of the Structural Division，1979，105（6）：1021—1037.

[30]  BERTRAND D，KASSEM F，DELHOMME F，et al. Reliability analysis of an RC member impacted by a rockfall using a nonlinear SDOF model[J]. Engineering Structures，2015，89：93—102.