第47卷第4期
2021年4月
北京工业大学学报
JOURNAL OF BEIJING UNIVERSITY OF TECHNOLOGY
Vol.47No.4Apr.2021
圆墩爆炸易损性分析
闫秋实1，吕辰旭1，李述涛2
(1.北京工业大学城市与工程安全减灾教育部重点实验室，北京 100124;2.军事科学院国防工程研究院，北京 100036)
摘 要：对桥梁在不同爆炸荷载下的损伤情况进行评估研究可以为桥梁设计㊁桥梁灾后评估及保险投资等提供参考，以完成更好的防护设计以及达到更好的经济评估效果.采用非参数化方法，对某典型高速公路连续梁桥圆形截面墩柱开展了爆炸易损性分析.建立了爆炸荷载作用下圆墩高效数值分析模型，定义了墩柱损
伤评价指标及墩柱损伤等级，对爆炸荷载作用下墩柱的反应样本进行统计分析，形成墩柱在不同爆炸当量的爆炸荷载下的易损性曲线.结果表明，当爆炸当量超过200kg 时，墩柱轻度毁伤水平的超越概率基本已达到1；当爆炸当量超过300kg 时，墩柱发生超过中度水平的损伤概率已接近80%.界限值为0.8的曲线在0~500kg 爆炸当量范围内呈现持续上升趋势.易损性曲线可以预测墩柱在不同爆炸当量发生不同损伤等级的概率，以更加准确地对墩柱损伤进行评估.关键词：桥梁墩柱；爆炸荷载；非参数化方法；有限元分析；损伤等级；易损性分析中图分类号:U 443.22
文献标志码:A
文章编号:0254-0037(2021)04-0394-09
doi :10.11936/bjutxb2020100017
收稿日期:2020⁃10⁃28
基金项目：国家自然科学基金资助项目(51678018)；北京市自然科学基金资助项目(8182009)
作者简介：闫秋实(1983 )，男，副教授，主要从事结构抗爆㊁防护工程方面的研究,E⁃mail:yqs2011@bjut.edu
Vulnerability Analysis of Circular Section Pier Under Explosion Load
YAN Qiushi 1,LÜChenxu 1,LI Shutao 2
(1.Key Laboratory of Urban Security and Disaster Engineering of Ministry of Education,
Beijing University of Technology,Beijing 100124,China;
2.Institute of Defense Engineering,AMS,PLA,Beijing 100036,China)
Abstract :Explosion vulnerability analysis of bridge structure provides a reference for bridge design,disaster assessment and insurance evacuation.Explosion vulnerability of a typical expressway continuous bridge circular pier was analyzed to use non⁃parameter method in this study.An efficient numerical simulation model for bridge pier under explosive load was established,and damage feature and pier damage grade were defined.The vulnerability curve of pier under explosion load was obtained through statistics to analyze response samples of pier column under explosion load.Results show that when the explosion equivalent exceeds 200kg,the probability of exceeding the level of slight damage to the pier has reached 1;when the explosion equivalent exceeds 300kg,the probability of exceeding the moderate level of damage to the pier is close to 80%.The curve with damage level of 0.8shows a continuous
汽车爆炸upward trend in the range of 0-500kg explosive equivalent.The vulnerability curve can be used to predict the probability of various damage levels for piers at various explosion equivalents to assess the damage of the pier accurately.Key words :bridge pier;explosion load;non⁃parameter method;finite element analysis;damage level;vulnerability analysis
桥梁是交通运输系统中重要的枢纽结构，和平
时期是交通物流的重要保障，战争时期则是战略物
第4期闫秋实，等：圆墩爆炸易损性分析
资及武器运输的重要通道.桥梁在服役过程中可能遭受地震㊁爆炸等灾害作用，地震及爆炸荷载作用都存在很强的不确定性，同时桥梁结构本身由于其结构形式㊁施工工艺及自然环境等原因，也存在较大的差异性及不确定性.为了考虑工程结构本身及其所承受的地震荷载不确定性影响，使结构在不同水准地震作用下的破坏和损失都能为业主接受，国内外学者对桥梁结构进行了大量易损性分析，并形成了完整的研究理论与应用.随着恐怖活动㊁战争㊁偶然爆炸事故的增多，桥梁结构遭受恐怖爆炸袭击㊁偶然爆炸及军事精准打击等爆炸荷载作用的概率逐渐增大，而中国关于爆炸荷载作用下桥梁易损性研究相对较少.欧美等国军方研究部门已建立了关于爆炸荷载作用下桥梁结构易损性分析理论，但是因为保密原因一直未予公开.如何结合中国桥梁的实际情况，建立能够综合考虑多种不确定性的爆炸荷载作用下桥梁结构易损性
分析方法，是当前亟待解决的难题之一.
国内外学者对桥梁结构开展过一系列地震易损性研究.Yu等[1]㊁Jernigan等[2]基于弹性反应谱方法对美国不同地区的高速公路桥梁开展了易损性分析，研究发现桥梁在地震作用下的非线性行为对结构易损性分析有较大影响;Dutta[3]㊁Shinozuka等[4]在后续的桥梁研究中开始采用静力非线性时程分析来对结构进行模拟，以考虑结构在地震作用下的非线性行为.Hwang等[5]采用参数化方法对美国中东部的高速公路钢筋混凝土连续梁桥开展了地震易损性研究，综合考虑了桥梁结构自身材料的不确定性与地震波的随机性，给出了桥梁地震易损性基本研究方法;Zhang等[6]采用参数化易损性分析方法对公路桥梁支座进行了隔振效率以及优化设计的评估研究.Subrat等[7]采用无须假定易损性曲线形式的非参数化易损性分析方法，对三层钢结构进行地震易损性分析；刘阳冰等[8⁃9]采用非参数化方法对2种不同类型方钢管混凝土框架结构进行了地震易损性分析，综合考虑了结构材料与地震输入的不确定性，基于蒙特卡洛随机抽样方法共计算1280个反应样本，对结构的易损性进行评估和分析.郑鸿强[10]㊁尤志华[11]基于典型战斗部以及弹道学相关理论以及经验公式对斜拉桥㊁钢筋混凝土桥梁开展了一系列毁伤易损性研究，但是研究中未涉及概率和不确定性问题.参数化方法与非参数化方法是易损性分析中2种重要的易损性曲线形成方法，非参化方法基于大量的样本数据进行统计分析，得出易损性曲线的相关参数；参数化方法需事先假定易损性曲线的形式来进行后续的分析[12].
综上所述，目前关于桥梁结构或构件易损性研究仍多集中于地震作用下的易损性研究，而爆炸载荷作用
下的易损性研究相对较少.本文针对圆形截面桥梁墩开展了非参数化易损性分析研究，可用于桥梁结构在爆炸荷载下的性能设计以及桥梁结构面对爆炸荷载的损失预估评价.同时也为以后结构的爆炸易损性分析开展与深入研究提供了一种研究思路.
1 爆炸易损性不确定性分析与损伤指标确定1.1 易损性分析中的不确定性分析
不确定性问题是结构易损性分析中非常重要的一个问题，也是易损性引入概率分析的基础.不确定性存在于一切事物当中，工程结构也不例外.工程结构本身存在工程材料强度的随机性，例如钢筋和混凝土材料强度的随机性；工程结几何尺寸的不确定性，如桥梁结构跨径㊁梁高㊁墩高的误差不确定性，以及桥梁地震易损性研究中考虑到的桥梁伸缩缝尺寸的不确定性等[12].此外，还有学者在桥梁抗震加固易损性研究中考虑到加固方式的随机性对桥梁易损性的影响[13].结构本身存在一系列的不确定性，作用在结构上的荷载同样存在诸多的随机性，如地震易损性分析中的地震波随机性等.在易损性分析中可以只考虑结构自身存在的不确定性或外荷载的不确定性，也可以将两方面引起的不确定性同时进行考虑.越多的不确定性引入到分析当中，意味着研究分析的全面性与更高的精度，但同时会造成分析样本数量与参数分析工作量的大量增加，这会造成后续有限元分析所需的时间㊁资源大量的增加.
对于桥梁墩柱而言，存在面临汽车意外爆炸㊁以及军事打击带来的爆炸风险.这些爆炸荷载同样存在诸多的不确定性问题，例如汽车意外爆炸的位置，不同类型车辆的等效爆炸当量；位置以
及袭击武器的攻击烈度；军事武器打击精度㊁打击条件限制造成的与目标点的距离偏差等.在本研究中，通过引入打击精度造成的爆炸位置的不确定性，开展对桥梁墩柱的易损性分析.军事打击中常规武器对桥梁此类工程目标打击中多采用触地引信，墩柱属于竖向构件，如果爆炸发生在桥面，则主要是桥梁梁体受到损伤，对墩柱影响不大；接触地面而发生的爆炸，对墩柱造成的伤害会更加严重.因
593
北 京 工 业 大 学 学 报
2021年
此本研究基于蒙特卡洛随机抽样方法，以墩柱截面几何中心为基本点(原点)，考虑战斗部触地爆炸，进行样本点随机抽取.高伟亮等[14]在对地面建筑的毁伤评估研究中发现，常规导弹的落点散布服从二维正态分布，落点计算式为
X =X 0+0.8493218E CP -2ln r 1sin (2πr 2)
Y =Y 0+0.8493218E CP
-2ln r 1cos (2πr 2{
)
(1)
式中:(X 0,Y 0)为目标瞄准点坐标;E CP 为落点散布的圆概率偏差，其表示出现概率为50%的圆形误差范围的半径，与打击武器的精度相关[15]，在本研究中根据目标打击武器的性能指标取E CP =5m;r 1㊁r 2为[0,1]的2个随机数.基于上述随机落点计算公式利用Python 程序对炸点位置进行二维随机抽样，共考虑5种爆炸当量，分别为100㊁200㊁300㊁400㊁500kg .
每种当量抽取30个随机炸点样本，共计样本150个.炸点位置示意图及生成的一系列落点坐标
分别入如图1㊁2
所示，落点的坐标如表1所示，对于落点位置位于墩柱截面半径R 之内的情况，按照接触爆炸工况计算.
图1 落点示意图
Fig.1 Diagram of target point
图2 随机生成炸点样本示意图
Fig.2 Diagram of randomly generated explosion center
samples
表1 落点坐标
Table 1 Coordinates of target points
i x i /cm
y i /cm
i
x i /cm
y i /cm
123456789
101112131415-217.101.30
302.48-82.39
-612.03-331.52512.93
-
401.54-545.31-151.92-838.35-227.71198.11465.91461.50-14.45-255.68-2.4814.29352.71787.66-824.0932.11225.18-182.29-78.69-18.35-411.91-464.88526.8816171819
2021
222324252627
282930
438.92437.25126.47400.5268.23
680.67-393.60191.10-227.70
133.48183.78208.2042.41250.49-171.59
-715.97117.39
-383.34-326.98423.42-44.15-488.95
805.15188.68-830.17-269.41145.60
-49.13-129.58
-87.641.2 损伤指标确定
钢筋混凝土结构在爆炸荷载作用下产生的损坏主要是以冲击波作用下混凝土的剥落㊁向内侵蚀，结构或构件的局部变形等局部响应为主，尤其是在近爆荷载作用之下.而结构在地震作用下，主要以整体响应为主，因此以墩柱为例采用譬如位移延性比㊁截面曲率等损伤指标可以较好地反映出墩柱在地震作用下的整体损伤情况.但在爆炸分析中，采用此类指标难以准确反映构件真实的损伤情况，这是由爆炸荷载的特性所决定的.在建筑结构构件的爆炸荷载损伤评估研究中[16]，对于结构柱此类以轴向承载为主的构件，采用基于残余轴向承载力的损伤系数作为损伤指标是比较常见的.残余承载力相关指标的相关研究虽尚不完善，但具有坚实的物理背景，直接体现轴向承载力的状况，易于工程人员进行参考，能够反映出构件的整体损伤情况.参考建筑构件在爆炸荷载下的毁伤评估研究成果，基于钢筋混凝土墩柱剩余承载力判断准则，提出损伤系数D 作为墩柱损伤的量化指标，损伤系数D 的定义为
D =1-P r P N
(2)
式中:P r 为墩柱经受爆炸荷载后的残余承载力，通过
有限元模型求解得出;P N 为墩柱的初始承载力，可由公式㊁试验或数值模拟求出，本研究中采用有限元模型求解得出.
6
93
第4期闫秋实，等：圆墩爆炸易损性分析2 墩柱分析模型与损伤水平划分
2.1 墩柱有限元模型计算
墩柱模型根据国内某高速公路连续梁桥实际墩柱简化而来.墩柱高h =10m，截面半径0.9m，混凝土强度等级为C40.纵筋配置为
42
28，配筋率为
1.016%.箍筋采用普通箍筋，箍筋配置为B10@100/200，如图3所示.考虑到桥梁墩柱真实的端部约束情况，模型墩柱底端采用全固定约束，顶端采用铰接约束.在柱顶面设置加载面加载初始轴压力以及后续的位移加载，初始轴压比取0.
15.
图3 墩柱模型示意图
Fig.3 Diagram of pier model
运用LS⁃DYNA 显式动力分析程序对模型进行爆炸加载，爆炸荷载加载完成后进行残余承载力求解，计算过程如图4所示.图4位置考虑到桥梁圆墩模型尺寸较大，网格单元过小会造成计算效率降低，在有限元模型中，混凝土采用实体单元，钢筋采用梁单元.二者单元尺寸均为10cm，采取钢筋-混凝土单元共节点的方式建立分析模型.模型中混凝土材料采用LS⁃DYNA 材料库中Mat _Concrete _Damage_Rel3(MAT072R3)本构模型进行模拟.该本构模型是K&C 本构模型最新的升级改进版本，是主要针对混凝土在爆炸㊁冲击作用下的动力响应而开发的材料模型[17].其可通过给定混凝土单轴无侧限抗压强度a 0，自动生成其他参数，使用方便.072R3本构模型可以较好地模拟材料的损伤演化过
程和评估结构的损伤破坏程度.该材料还可以采用动力放大系数(dynamic increase factor,DIF)来考虑材料在冲击㊁爆炸等高应变率情况下的材料强度增强效应.本研究采用Malvar 等[18]
提出的动力放大
图4 计算过程示意图
Fig.4 Diagram of calculation process
系数计算方法来考虑应变率效应，抗压强度动力放大系数F CDI 和抗拉强度动力放大系数F TDI 计算公式分别为
F CDI =f dc f c =(ε㊃/ε㊃s )1.026α,ε㊃≤30s -1γ(ε㊃/ε㊃s )1/3,
ε㊃≥30s
{-1
(3)F TDI
=f dt
f t =(ε㊃/ε㊃
s )δ,ε㊃≤1.0s -1β(ε㊃/ε㊃s )1/3,
ε㊃
≥1.0s {
-1
(4)
式中:ε㊃
表示混凝土的应变率，混凝土在抗压时取值范围为10-6~160s -1，在抗拉时取值范围为3×
10-5~300s -1;ε㊃
s 表示混凝土的静态参考应变率，其
抗压时取值为3×10-5s -1，抗拉时取值为10-6s -1;f dc ㊁f dt 分别表示混凝土的动态抗压㊁动态抗拉强度;f c ㊁f t 分别为混凝土的静态抗压强度㊁静态抗拉强度;lg γ=6.156α-2,α=1/(5+0.9f c );lg β=6δ-2;δ=1/(1+0.8f c ).
钢筋采用MAT _PLASTIC _KINEMATIC (Mat _
003)材料模型，该模型通过给定钢材的屈服强度㊁切线弹性模量以及泊松比等参数来确定.在该模型中，可以通过Cowper and Symonds (C⁃S)应变率模型[19]来考虑应变率效应对钢材动态力学性能的影响，动力放大系数的计算式为
F DI (=1+
ε㊃
)
C
1/p
(5)
式中:ε㊃
是钢材的应变率;C 和p 是C⁃S 模型的应变率参数,C 和p 的取值与钢材种类有关.
本研究涉及到结构的爆炸荷载的加载,LS⁃
DYNA 是爆炸冲击领域具有高知名度的分析程序.
常用的混凝土结构爆炸荷载加载方式有流固耦合方法㊁简化荷载方法以及LS⁃DYNA 软件内置的LOAD_BLAST 加载方式等.LOAD_BLAST 加载方法设置简
单，仅需建立结构模型，无须建立实体与冲击波传播介质(默认空气)，内存占用小，计算速度快，但7
93
北 京 工 业 大 学 学 报2021年
无法考虑冲击波与结构相互作用发生的绕射.流固耦合加载方法可以较好地模拟出爆炸波的产生以及与结构的相互作用，更接近真实的爆炸反应情况，计算误差较小[20].但该方法需要消耗大量的计算内存资源与计算时间，对于本研究需要开展大量的工况计算来讲并不合理.耦合Load Blast 与ALE 加载方法(coupling Load Blast enhanced with ALE method,CLAM)[21]方法是一种在之前的爆炸模拟研究中使用相对较少的方法，其本质上也是一种流固耦合加载方法.CLAM 其方法通过在目标结构外层建立小范围的ALE 单元域，将ALE 区域冲击波入射面外层单元设置为接收单元(ambient element,AE)，接收由LS⁃DYNA 内置的Load Blast 程序生成的爆炸波压力数据从而实现冲击波与结构的相互作用而不必建立大规模的ALE 网格区域以减少计算
占用，对于本研究所面对的爆炸模拟工况来讲具有合理性.但是CLAM 方法与Load Blast 计算方法均受Load Blast 程序对于工况比例距离适用范围0.178~40.000m /kg 1/3的限制[22]，故在本研究中对于一般工况均采用CLAM 方法计算，对于接触爆炸等少量特殊工况采用流固耦合方法进行求解.本研究根据墩柱
构件尺寸建立如图5所示的CLAM 模图5 CLAM 模型示意图
Fig.5 Diagram of CLAM model
型，少量流固耦合工况模型基本同CLAM 模型，区别仅在于空气域的扩大以及建出实体.求解过程共分为P 1~P 4四个阶段，如图4所示，首先通过动力松弛算法(dynamic relax)加载初始轴压力P 0，之后加载爆炸荷载，爆炸荷载引起的构件震荡也是造成构件损伤的原因之一故待爆炸荷载加载完毕后允许墩柱自由震荡一段时间，之后通过重启动引入结构阻尼，使墩柱趋于平稳，开始在柱顶通过刚性加载板加载位移，求解残余轴向承载力，其中刚性加载板不参与流体耦合计算.墩柱的轴向承载力通过定义柱底位置的墩柱横截面并输出横截面轴向力曲线得到.圆形墩柱横截面为完全轴对称图形，可以将表1所示的不同落点的坐标通过计算转化为等效爆炸距离来进行加载，对于不同位置但距离相同的落点可以按照同一工况来进行计算以减少计算时间提高计算效率.表1中所示的共计30个落点可折算为如表2所示的27个工况.
表2 落点距离
Table 2 Distance of target point
j d i /m j
d i /m j
d i /m 123456789
0.60.81.92.22.32.42.52.62.8
101112131415161718
3.03.03.3
4.04.34.54.6
5.25.9
192021222324252627
6.36.66.8
7.07.1
8.38.48.5
9.7
由于计算工况较多，在此不再一一列举各工况计算所得的损伤情况，图6~14给出部分工况图6 100kg 接触爆炸损伤图
Fig.6 Damage of pier under a contact explosion of 100kg
的墩柱损伤图与轴力曲线.100kg 接触爆炸工况的墩柱损伤情况如图6所示，可见墩柱下部在爆轰产物的冲击下，已经完全损坏，完全丧失承载力，这与图7所示的轴向力曲线完全一致，墩柱在遭受爆炸荷载后无法继续承载，残余轴向承载力为零.如图8~10所示在等距离的爆炸工况下,
200kg 爆炸当量承载力测试曲线的峰值远高于100kg，前者的损伤情况也明显较后者严重，构件8
93
发布评论