Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

双钢板混凝土组合板抗爆性能分析

赵春风 何凯城 卢欣 潘蓉 王静峰 李晓杰

张娅, 李晓彬, 彭帅, 施锐. 绝热指数γ对平面爆轰过程中不同复合波区的参数特性影响分析[J]. 爆炸与冲击, 2021, 41(10): 102102. doi: 10.11883/bzycj-2020-0458
引用本文: 赵春风, 何凯城, 卢欣, 潘蓉, 王静峰, 李晓杰. 双钢板混凝土组合板抗爆性能分析[J]. 爆炸与冲击, 2021, 41(9): 095102. doi: 10.11883/bzycj-2020-0291
ZHANG Ya, LI Xiaobin, PENG Shuai, SHI Rui. Influences of adiabatic index γ on the parameters of different complex wave zones in a planar detonation[J]. Explosion And Shock Waves, 2021, 41(10): 102102. doi: 10.11883/bzycj-2020-0458
Citation: ZHAO Chunfeng, HE Kaicheng, LU Xin, PAN Rong, WANG Jingfeng, LI Xiaojie. Analysis on the blast resistance of steel concrete composite slab[J]. Explosion And Shock Waves, 2021, 41(9): 095102. doi: 10.11883/bzycj-2020-0291

双钢板混凝土组合板抗爆性能分析

doi: 10.11883/bzycj-2020-0291
基金项目: 安徽省自然科学基金(2008085UD12);工业装备与分析国家重点实验室基金(GZ19106)
详细信息
    作者简介:

    赵春风(1983- ),男,博士,副教授,zhaowindy@hfut.edu.cn

  • 中图分类号: O383

Analysis on the blast resistance of steel concrete composite slab

  • 摘要: 钢-混凝土-钢组合板是一种新型的组合结构,与传统钢筋混凝土板相比,具有抗剪强度高、延性大、耗能能力强等特点,目前已经被广泛应用于核反应堆安全壳、海洋平台及储油罐等结构。本文中,设计并制作了缩尺的普通钢筋混凝土板和钢-混凝土-钢组合板,开展了在接触爆炸荷载作用下的实验研究,通过损伤分析、跨中最大挠度对比研究不同板的抗爆性能。基于ANSYS/LS-DYNA非线性有限元程序,研究了钢-混凝土-钢组合板的损伤模式、跨中最大挠度等,并与实验结果进行了对比分析,验证了有限元分析模型的准确性和适用性。参数化分析了炸药量、混凝土强度和钢板厚度等参数对钢-混凝土-钢组合板抗爆性能的影响规律。利用多参数回归分析方法,提出钢-混凝土-钢组合板跨中挠度的预测公式。结果表明:提高混凝土强度可以降低结构的塑性损伤, 增加钢板厚度可以有效降低钢-混凝土-钢组合板的跨中最大挠度。相对于普通钢筋混凝土板,钢-混凝土-钢组合板保持了良好的整体性,且具有继续承载的能力。拟合公式能够较好地预测钢-混凝土-钢组合板跨中挠度与药量和钢板厚度的关系。
  • 爆炸冲击波的传播与衰减一直是爆炸领域研究的重点问题[1]。当炸药爆轰后,爆炸气体产物内部形成的稀疏波会迅速衰减爆炸高压,其中不同方向的稀疏波会在中心汇聚,形成相互作用的稀疏波复合波区,稀疏效应的叠加加快了爆炸高压的衰减速率,且穿过复合波区的稀疏波会快速追赶上空气冲击波阵面,对该阵面的压缩特性进行衰减[2]。因此,在很多冲击波问题的研究中,为在实验段得到稳定的冲击超压,通常会在激波管实验中通过加长爆轰驱动段[3],来减缓反射稀疏波到达激波阵面的时间。可以看出,该稀疏复合波区的特性直接影响了爆炸过程的衰减特征。除稀疏复合波区外,随着流场运动,爆炸区域还会形成其他类型的相交复合波区,且每个波区的传播和衰减特征均不同。

    在平面爆轰流场中,该类复合波区均可以简化为两相向扰动波相交作用的特性耦合。在应用特征线法的爆炸理论中,已提出了两相向扰动波交汇后流场的一般解形式[4],但由于相交特性复杂,并没有通用的解析解,仅在假定爆炸气体绝热指数γ=3时有确定的特性,即相交后两簇波在x-t平面上的特征线仍均为直线。张守中[5]基于该γ=3时的特殊解形式,对多个典型爆轰产物流场中稀疏复合波区的数值解进行了计算,并将计算解代入到后续的流场分析。可以看出,对于复合波区而言,该特殊解形式简单,方便用于流场参数的快速估算。因此,当对计算精度要求不高时,对高密度凝聚炸药常近似取γ=3.0。

    但是在实际过程中,爆炸气体内部这种从极稠密的高压状态变化到较稀疏的低压状态,过程是极其复杂的,尤其是过程中爆炸气体的绝热指数γ是持续变化的。常用高密度凝聚炸药的初始γ0测定值一般在2.2~3.50之间[6-8],在根据γ律等熵状态方程进行膨胀衰减计算时,常采用在高压γ0到低压(1.1~1.4)分段取值计算的方法[9-10],所以γ≠3是这些复合波区的常态。

    为研究γ=3和γ≠3不同条件下爆炸不同复合波区的特性差异,本文中基于特征线法和理想气体的γ律等熵状态方程,对一平面爆轰过程中不同复合波区的波系相交特性进行理论规律分析。由于在一般解形式中,γ≠3的复合波区内相交的特征线为迹线未定的曲线,难以推导理论的位移-时间公式(即x-t解),本文中利用MATLAB对其平面爆轰过程进行流场模拟,验证并分析不同复合波区流场内的参数变化特性。

    在如图1所示的截面积为1 m2的圆形无限长管内,T0时刻在x(−1,1)区间内设有TNT固体装药。基于瞬时爆轰模型,在炸药原体积空间内,爆炸气体产物的密度取为固体装药的密度,初始压力取pCJ爆压乘以等效系数0.434 8[11],初始质点速度为零。

    图  1  平面爆轰过程中复合波区的形成及相对位置(+右行,−左行)
    Figure  1.  Formation and relative positions of complex wave zones in plane detonation (+ showing a going-right wave, − showing a going-left wave)

    利用特征波系的运动分析该爆炸流场,当t>0时,阶段1:空气受压缩形成冲击波阵面,在内部中心稀疏波的作用下爆炸气体压力迅速衰减;阶段2:原爆炸气体产物(③区)全部被稀疏扰动,两相向中心稀疏波相遇并相交作用,形成稀疏波-稀疏波复合波区(⑤区);阶段3:稀疏波波头完全扰动④区,继续向前扰动②区(扰动后形成稀疏⑥区),在到达与空气的交界面后,由于爆炸气体的波阻抗大于空气的波阻抗,会在空气中形成透射稀疏波,并向爆炸气体中反射一个压缩波,该压缩波又与原稀疏波相交作用,形成稀疏波-压缩波复合波区(⑦区)。

    为分析阶段2、阶段3中生成的两复合波区的参数特性,首先利用特征线法确定流场的初始波动特征(阶段1)。为便于分析流场中的能量流动,计算中假定高压气体和空气均为无黏理想气体,满足理想气体状态方程和γ律等熵方程。绝热指数γ不同,整个流场的参数分布不同。

    对于TNT炸药的pCJ爆压测定值一般在19~21 GPa,以pCJ=21 GPa时为例,爆热e为4 187 kJ/kg,则γ=2.34。根据类激波管问题的特征线理论,可以得到图1中阶段1时段无限长管内各区域的状态参量,见表1,且①区前端的冲击波阵面速度vs=5 266.14 m/s。当左右两稀疏波在原点相遇时,阶段1结束,记此时刻为T1,此时无限长管内各区域的状态参量分布,如图2所示,图2(a)(b)(d)为对称半流场,图2(d)Ed为线能量密度,Ed-intEd-kin分别代表内能密度、动能密度,沿管长度方向积分可得该时刻系统的内能及动能。

    表  1  阶段1各区域的状态参量
    Table  1.  Characteristic parameters of each region in stage 1
    区域压力p/
    MPa
    密度ρ/
    (kg·m−3
    质点速度u/
    (m·s−1
    气体声速c/
    (m·s−1
    绝热指数γ
    ③区9130.8 163003 620.5002.34
    ④区介于③区和②区之间,参数按中心稀疏波规律衰减2.34
    ②区 28.29 138.054 730.127 692.5102.34
    ①区 28.29 7.204 730.1272 345.5671.40
    ⓪区 0.101 1.2250 340.2941.40
    下载: 导出CSV 
    | 显示表格
    图  2  两稀疏波在原点相遇时流场内的特征参数分布
    Figure  2.  Distribution of characteristic parameters when two rarefaction waves meet at the origin

    由特征线理论可知,两个相向行波相交满足黎曼不变特性,即右行波的每一条C+特征线,沿该线的黎曼不变量J+为定值;左行波的每一条C特征线上,沿该线的黎曼不变量J为定值。对于理想气体,其J±=u±2c/(γ−1)。由此可得任一条C+C特征线相交后,交点位置Muc参数满足:

    u=12(J++J),c=γ14(J+J) (1)

    由于中心稀疏波特有的中心发散特性,其运动(④区)有通解方程。此时若一条C+C特征线穿过常态区与一相向中心稀疏波相遇,其迹线的解析方程为:

    x=γ+1γ1xc(cconsxct)3γγ+1±2γ1cconstx0 (2)

    式中:xc为中心稀疏波发散中心,x0为相向特征线出发点,ccons为常态区气体声速。

    图3图1所示流场的x-t平面示意图,图中⑤区就是两稀疏波特征线的相交区域,在交汇前各单波系的特征线运动均是直线。根据公式(2)可以得出⑤区边界线即图3中左右波头特征线C+0C−0x-t迹线。当γ=3时,式(2)变化为一直线形式。当γ≠3时,该x-t解为一明显的多项式曲线形式,可以看出当γ>3时C+0C−0的波速是不断减小的,曲线向内凹;γ<3时C+0C−0的波速是不断增大的,曲线向外凹。

    图  3  两中心稀疏波相交的x-t平面
    Figure  3.  x-t plane of intersection of two central rarefaction waves

    而对于复合波区内部任意两条相交特征线的x-t方程,其一般解形式[2]为:

    {C+:x=(u+c)t+F+C:x=(uc)t+F (3)

    式中:F+F分别为与特征线波速u+cu−c相关的系数。由式(1)可以看出,特征线相交会影响其运动波速,在复合波区内一条特征线会受N条相向特征线的影响,其x-t迹线变为斜率不断变化的曲线,所以式(3)很难得到直接的解析解。但由于在γ=3条件下,dx/dt=u±c=J+(或J)=常数,所有特征线均沿直线向前推进,所以可以通过左右行波的初始条件确定F+F

    图3中两中心稀疏波相遇时,由右行稀疏波t=0时,x=−1;左行稀疏波t=0时,x=1;可得F+=−1,F=1。可解得复合波⑤区任一交点的uc参数:

    u=xt,c=1t (4)

    对于理想气体而言,该点的压力、密度均可通过该气体声速得到(其中p3ρ3c3为③区参数):

    p=p3(cc3)2γγ1=p3(1c3t)2γγ1, ρ=ρ3(cc3)2γ1 = ρ3(1c3t)2γ1 (5)

    所以在γ=3时复合波⑤区内各处气体的声速c、压力p、密度ρ都是一个均布值,与气体分布位置无关,仅随运动时间t变化,并且随着t增长,cpρ逐渐减小,而气体的质点速度u同时受位置和时间的影响。可以看出该特殊解形式简单,且方程可解析。但当γ≠3时,复合波区内的各特征线均为变化曲线,任一交点的uc参数无通解。图4显示了近似利用γ=3的c值均布特性求解后原图2γ=2.34时流场的特性分布变化,其中E为能量,Eint为内能,Ekin为动能。

    图  4  应用c值均布特性后系统流场特性变化
    Figure  4.  Changes of flow field characteristics after the application of c-value uniform distribution

    图4(a)可以看出,复合波⑤区内除气体质点速度外,内部的其他参数都是均布状态。图4(b)计算了两稀疏波从初始T0、相遇T1到④区完全被扰动为复合波区的T2时刻过程中各区域的能量变化,可以发现,系统内的总能量在复合波⑤区形成后逐渐减少,到T2时刻总能量已减少了37.53%。这说明,当γ≠3时其复合波区的特性应与γ=3时有较大差异。

    从波系相交特性的角度分析,如图5所示,当γ≠3时选取3条右行C+特征线与左行中心稀疏波相交,C+0线按已有解析方程式(2)画出。由于图1中③区为常态区,根据特征线理论,该左行波C族线的右向黎曼不变量J+C+0线的J+0相同,所以在C+0线通过后,各C线的uc不变,仍按原斜率u-c运动,此时左行波仍保持中心发散特性。

    图  5  中心稀疏波的发散特性
    Figure  5.  Divergence characteristics of central rarefaction waves

    因此当C+1线(J+1J+0)与运动而来的C族线相交时,其特征迹线方程仍可按C+0的迹线求解方式得到,形式与式(2)类似。但由式(1)可知,与C+1线相交后,各C线的运动斜率(uc=3γ4J++γ + 14J)都将发生变化,将其延伸到x轴,可以看出此时的C族线并没有相交到一点,即中心稀疏波在受任一相向运动的特征线(其同向黎曼不变量与中心稀疏波的不同)影响后,将不再具有中心特性。

    所以当C族线再与C+2线相交时,C+2线已不能够利用式(2)方式求得迹线方程,只能通过与各C线逐一相交迭代计算。同理可推,当γ≠3时在图3中相向运动的两个中心稀疏波相交后,C+族线和C族线都不再具有中心特性,变为非中心稀疏波,即其波系发散特性发生了改变,因此难以确定其各特征线相交的具体位置。

    由上可知在x-t平面上无法直接对比其特性差异,改从能量角度分析,理想气体的内能与压力p相关,动能与密度ρ、质点速度u相关,结合式(5)可知,该系统的能量特性主要受uc变化的影响,所以为分析其特性差异,可以从复合波区的u-c特征进行切入。

    由特征线理论可知,对于每一个不变量J+/J,在u-c平面上必有一条曲线Γ+/Γ与之对应,且曲线Γ+Γ的交点处的(u,c)值与C+C特征线相应交点的(x,t)值具有一一对应关系。由于理想气体的等熵特性,u-c平面上的Γ+/Γ曲线均为直线。

    (1)复合波⑤区的u-c平面特性

    图6为两同等强度中心稀疏波相交的x-t平面和u-c平面。假定tk时刻,x-t平面上左行波波头C−0C+k线相交(A点),由对称可知,此时右行波波头C+0Ck线相交(B点),线AB为复合波区的长度,对应到u-c平面,可以找到相应的复合波区边界点AB,当γ=3时,x-t平面上特征线为直线,复合波区中点位置为P点,u=0。根据式(4)已知c与分布位置无关,为均值,因此在u-c平面,AB之间的路径为直线A-P-B

    图  6  两同等强度稀疏波相交的平面特性
    Figure  6.  Plane characteristics of the intersection of two rarefaction waves with equal intensity

    (2)复合波⑦区的u-c平面特性

    图1中复合波⑦区的类型是稀疏波和压缩波的相互作用。为分析⑦区的特性,首先了解两同等强度的压缩波相遇的复合波特性。在图7(a)中分别在长管内−x0x0处设置活塞,管内气体静止,初始声速c0。利用活塞相向匀加速运动,加速度为a,压缩管内气体,生成两相向等强度压缩波。两压缩波运动到中点(x=0)处相遇,形成压缩波-压缩波的复合波区。

    图  7  两同等强度压缩波相交的平面特性
    Figure  7.  Plane characteristics of the intersection of two compression waves with equal intensity

    同样利用式(3)进行图7复合波特性的求解。当γ=3时,由于dx/dt=u±c为常数,所有特征线仍均沿直线向前推进,活塞压缩生成的每条特征线的初始参数为(x+i, t+i, u+i, c+i)和(xi, ti, ui, ci),其中下标+为右行波,−为左行波。该区相交特征线的解为:

    {x=(u+c)t+F+x=(uc)t+Fc=FF+2t (6)

    由和常态区相邻的简单波系的黎曼特性可得:

    {u+i2c+i/2c+i((γ1)=2c0/2c0((γ1)ui+2ci/2ci(γ1)(γ1)=2c0/2c0(γ1)(γ1)u+i+c+i=2u+i+c0uici=2uic0 (7)

    此时假定tk时刻,C−0C+k线相交(A点),C+0Ck线相交(B点),线AB为复合波区的长度。为确定此时AB线上所有交点的声速c的特性规律,任取一交点MxM, tk)的方程展开:

    {xM=32at2+M+(2atkc0)t+M+c0tkx0xM=32at2M(2atkc0)tMc0tk+x0 (8)

    通过式(8)可得到线AB上相交两条特征线初始参数t+MtM的关系式:

    2(c0tkx0)=32a(t2M+t2+M)(2atkc0)(tM+t+M) (9)

    可以看出变量(tM+t+M)和因变量(t2M+t2+M)满足S=S1yS2x形式的直线方程,同时式(9)可以改写为:

    2c0tk+2atk(tM+t+M)=32a(t2M+t2+M) + c0(tM+t+M)+2x0 (10)

    式(10)等号右侧正好是式(6)解中FF+的展开表达式,将M点和B点的声速c进行对比可得:

    cMcB=FMF+MFkF+0=2c0tk+2atk(tM+t+M)2c0tk+2atk(tk+t+0) (11)

    从式(11)可以看出cM/cMcBcB1,因此在该压缩波复合波区的u-c平面上AB路径并不是一条直线。将ac0tk赋值后发现,(tM+t+M)在中心P点最大,在两端AB处最小,因此cPcB。由于在u-c平面上压缩波和稀疏波波头波尾走向相反,因此该波区的u-c平面如图7(b)所示,其AB路径为曲线A-P-B

    将该分析方法延伸至图1复合波⑦区的类型分析中,在长管内利用图8(a)中活塞向左匀加速运动生成不同强度(a1a5)的压缩波与迎面稀疏波相交,两相向波在x-t平面上的初始位置分别为(x0,0)、(x0,0)。当γ=3时结合式(4)和式(6)的分析,可以得到图8波区交点的F+F-解为:

    图  8  不同强度压缩波和稀疏波相交的平面特性
    Figure  8.  Plane characteristics of the intersection of a compression wave and a rarefaction wave with different intensities
    F+i=x0,Fi=xi(uici)ti (12)

    Fi展开,由此推出该波区任意时刻的声速解c

    c=FiF+i2t=32at2i+c0ti+2x02t (13)

    由式(13)可以看出,在指定a和时刻t后,该波区范围内c仅与各压缩线的初始形成时间ti有关,因此t时刻由压缩波波头到压缩波波尾c逐渐增大。但对于单条压缩线来讲,在复合波区内其c值随t的增长逐渐减小,这就是稀疏波对压缩波稀疏效应的体现。

    对于该路径上u的特性,根据稀疏波方程x=(u+c)t+F+可得:

    u=x+x0tc (14)

    由于在t时刻,(x+x0)/t项随x增大而增大,−c值会随x增大而减小,而c值变化速率主要受a的影响,因此该路径上u的特性也主要受活塞加速度a的影响,即受压缩波强度的影响。假定tk时刻,该复合波区边界交点为AB,改变活塞初始加速度(a1a5),以式(13)~(14)对应计算图8(b)u-c平面上的AB路径。计算得出,当a较小(a1a3)时,−c项的影响较小,AB路径上u基本呈单调递增的状态,但当a继续增大(a4a5)时,−c项的影响开始增大,AB路径上u开始有明显的先减小后增大的趋势,若继续增大a,可预见递减会成为主要趋势。

    综上分析可以看出,当γ=3时各类复合波区的c值均只跟时间有关,但在不同类型复合波区中分布形式不同,而u值与xt均相关,在同类型波相交的复合波区内是单调递增的,若相交波为等强度波,对称面u=0。在不同类型波相交中u值特性与参与波的强度直接相关。uc特性的组合就是复合波区的主要特性。

    虽然γ≠3时依据式(3)流场难以求解,但对于具有对称特性的复合波⑤区,仍可近似借鉴γ=3时的直线特性,假定其特征线也是沿直线运动,同样得到式(4)的c值均布特性(当然此假定过程并不等熵),然后放在原γ≠3的x-tu-c平面上进行对比和特征规律分析,如图9所示。

    图  9  γ≠3时两同等强度稀疏波相交的平面特性示意图
    Figure  9.  Plane characteristics of intersection of two rarefaction waves with equal intensity at γ≠3

    图9中假定tk时刻,A0-P0-B0为假定特征线沿直线运动时的c值均布直线。而其实际情况由式(2)及图3可知,当γ≠3时,x-t平面上复合波区内的特征线均为曲线,此时需分两种情况:

    (1)当γ>3时,波头C+0/C−0迹线的波速是不断减小的,所以波头C+0/C−0会比直线运动时更晚与Ck/C+k线相交(tk1tk),x-t平面上线A1B1为复合波区的长度,中点位置为P1点,u=0。但在u-c平面上复合波区边界点仍为A0B0两点,而P1点在x-t平面上高于P0点,也就是在u-c平面上其相交特征线的编号高于P0点,因此当γ>3时,在u-c平面上A1B1之间的路径为曲线A0-P1-B0,即此时复合波⑤区的特性参数变化规律为:c(或pρ)在复合波区中心处最小,向着复合波区边界逐渐增大。u为矢量,在对称中心处均为0,数值上向着复合波区边界逐渐增大。

    (2)当γ<3时,波头C+0/C−0迹线的波速是不断增大的,所以波头C+0/C−0会比直线运动时更早与Ck/C+k线相交(tk2tk),x-t平面上线A2B2为复合波区的长度,中点位置为P2点,u=0。但在u-c平面上复合波区边界点仍为A0B0两点,而P2点在x-t平面上低于P0点,也就是在u-c平面上其相交特征线的编号低于P0点,因此当γ<3时,在u-c平面上A2B2之间的路径为曲线A0-P2-B0,即此时复合波⑤区的特性参数变化规律为:c(或pρ)在复合波区中心处最大,向着复合波区边界逐渐减小。u为矢量,在对称中心处均为0,数值上向着复合波区边界逐渐增大。

    从上述分析可以看出,当γ>3和γ<3时复合波⑤区的u-c路径特性是相反的。

    为验证上述结论,也为得到非对称复合波⑦区在γ≠3时的u-c平面特性,本文利用MATLAB基于特征线法对图1所示的爆轰流场进行模拟,首先假定左右中心稀疏波④区各有N条特征线从中心发出,并在每个区域预先设置特征空矩阵,以便于下一区域直接调用上一区域的边界作为初始条件。为保证计算精度,计算中采用尽可能多的特征线充满流场,取N=10 000。

    图1中阶段3有2种不同的复合波形式,对于稀疏波相交的复合波⑤区,以公式(1)为迭代准则,假定在每个迭代胞格内特征线为直线,计算胞格各交点的(x,t)及(u,c),再推进到下一个胞格,如图10(a)所示;对于在介质分界面处的复合波⑦区,其入射波、反射波和透射波的关系[12]按下式计算:

    图  10  对复合波区的迭代思路
    Figure  10.  The iterative method for the complex wave zone
    u7u2=2ψexplψexpl+ψair(u6u2),p7p2=2ψexplψexpl+ψair(p6p2) (15)

    式中:ψexplψair分别为爆炸气体和空气的波阻抗。由接触面两侧的连续性可知,入射前p2=p1u2=u1;入射后p7=p8u7=u8。计算过程中入射特征线要先与前入射特征线生成的反射特征线相交后才能到达介质分界面,生成新的反射、透射特征线,如图10(b)所示,该过程中分界线DL也会变为曲线,要在迭代中要不断更新其位置,并对压缩波特征线增加判断语句,若压缩波特征线开始相交,则计算报错。由于压缩波波线汇聚处会形成间断面,属于波形畸变,此时特征线法不再适用。

    为充分讨论γ≠3的各种情况,本文中通过调节初始爆热,得到7组不同γ值:3.5、3.2、2.7、2.34、2.0、1.7、1.4。由于该流场程序具有通用性,因此增加γ=3工况作为对比。为验证该流场程序的合理性,以γ=2.34时为例,MATLAB计算生成的特征线流场如图11所示。每个特征线上节点位置处都包含了6个特征参数:压力、密度、质点速度、气体声速、内能密度、动能密度。选定时刻,即可确定该时刻内流场各处的参数分布。计算过程在Tn=5.93 ms时刻报错,显示⑦区N7=173与N7=172反射压缩特征线在26.569 5 m处相交。

    图  11  MATLAB计算生成的特征线流场
    Figure  11.  Characteristic line flow field generated by MATLAB calculation

    首先选定3个特征时刻进行流场对比:T2时刻(右行稀疏波波头刚穿出复合波⑤区)、T3时刻(右行稀疏波波头刚穿出②区,即到达介质分界面)、Tn时刻(计算结束时刻)。流场参数分布如图12所示,由于特征线流场左右对称,只显示右半流场。

    图  12  T2T3Tn时刻的流场参数分布
    Figure  12.  Distribution of flow field parameters at times T2, T3 and Tn

    图12(d)各区域随时间的能量变化曲线可以看出,Tn时刻之前流场的总能量保持守恒状态,说明该MATLAB迭代流场是合理可行的,若忽略Tn时刻报错继续计算,流场总能量开始增大。

    (1)复合波⑤区的u-c平面特性

    图11x-t平面可以看出,中心稀疏波相交作用的复合波⑤区在T1时刻开始形成。图13(a)为该区与图11x,t)交点对应的u-c平面,从图中不同时刻的u-c路径及图12(a)(c)中不同时刻⑤区内pρc的分布规律,可以看出在该波区内c(或pρ)均在波区中点处最大,向着两侧波区边界逐渐减小。这与上述γ<3时复合波⑤区的u-c特性分析是一致的。图13(b)γ不同时各流场中复合波⑤区在各自T2时刻的u-c路径,从对比中可以看出:不同流场的u值特性是一致的,均是沿着波区长度方向逐渐增大;而其c值特性则随着γ增大,其沿波区的路径逐渐由上凸曲线(γ<3,中心c值最大)-直线(γ=3,c值均布)-下凹曲线(γ>3,中心c值最小)变化,且离γ=3值越远时,曲线特征更明显。当γ值越接近3时,其曲线路径越接近直线特性。

    图  13  爆轰流场中复合波⑤区的u-c特性
    Figure  13.  u-c characteristics of complex wave zone ⑤ in the post-detonation flow field

    由于在爆轰流场中,复合波⑤区范围最广,且前期衰减速率及衰减跨度大,如对比γ=2.34时流场中图2图12,波区中点p值从T1时刻9 130.8 MPa快速衰减到了T2时刻的44.07 MPa,该时段内中心ρ值从1 630 kg/m3衰减到了166.8 kg/m3,因此该c值特性不同影响的pρ参数对系统的内能、动能均产生了较大影响。

    (2)复合波⑦区的u-c平面特性

    对于在T3时刻开始形成的复合波⑦区,通过图11Tn时刻压缩波特征线交点的位置可以看出,相交位置位于靠近介质分界面,也就是靠近反射压缩波波尾处。若将介质分界面DL看成一个可透射的活塞,若活塞向左做均匀加速、且波前静止的情况,间断出现的位置一般在波头上,而该⑦区流场中间断率先发生在了波尾处,说明后生成的反射压缩线追赶速率不断增大,也就是可透射的活塞DL向左加速度是不断增大的,间断将由波后向波前推进。在其他7组不同γ的流场中均出现了类似的间断位置,而且随着γ增大间断位置出现的越早。因此可以得出,该平面爆轰过程中的复合波⑦区是由右行稀疏波(其中γ≠3时为非中心稀疏波,γ=3时为中心稀疏波)和左行变加速反射压缩波相交形成的。

    图14(a)γ=2.34时该波区与图11x,t)交点对应的u−c平面,从c值特性看,在不同时刻该波区范围内c值均是沿压缩波波头向波尾(介质分界面)逐渐增大的,这与γ=3时的规律特性是一致的。从u值特性看,在T7时刻前u值沿着波区主要呈递增的趋势,但随着时间增长,u的增幅越来越小,且结合图12(c)中⑦区的u-x曲线,可以看出在Tn时刻该波区的u值出现了先减小后增大的变化趋势,由此也可以看出该反射压缩波的强度随生成时间是不断增大的。不过从整体上该波区内各位置的uc还是随时间不断衰减的。

    图  14  爆轰流场中复合波⑦区的u-c特性
    Figure  14.  u-c characteristics of complex wave zone ⑦ in the post-detonation flow field

    图14(b)γ不同时各流场中复合波⑦区在各自Tn时刻的u-c路径,从对比中可以看出:不同流场的c值特性是一致的,均是沿着波区长度方向逐渐增大,而其u值特性随γ增大逐渐由单调递增变化为先减小后增大,该先减小的特性也间接表现了右行稀疏波的强度,说明随γ增大,流场生成的稀疏波强度越大。

    可以看出,在该平面爆轰过程中不同区域都有其特定的衰减特性,根据其特性可以快速了解整个流场的参数分布情况,如复合波⑤区的边界点是流场压力最低的位置。同时若将图11模型继续计算,在Tn时刻之后约1 ms的时间,⑧区的右行透射稀疏波将与冲击波波阵面交汇,还会形成稀疏波和冲击波的复合阵面,此时空气超压和冲击波速开始被衰减。

    综上分析,uc平面特性的差异就是绝热指数γ不同时复合波区衰减特性不同的主要体现。

    (1)当γ≠3时,中心稀疏波在受任一相向运动的特征线(其黎曼不变量与该中心稀疏波不同)影响后,不再具有中心特性;

    (2)对于等强度中心稀疏波相交的复合波⑤区,主要表现在c值特性不同:当γ=3时,c(或pρ)在复合波区内为均布值,仅与时间相关;当γ>3时,c(或pρ)在复合波区中心处最小,向着复合波区边界逐渐增大;当γ<3时,c(或pρ)在复合波区中心处最大,向着复合波区边界逐渐减小;

    (3)对于稀疏波与压缩波相交的复合波⑦区,主要表现在u值特性不同:沿波区长度方向u值特性随γ增大逐渐由单调递增变化为先减小后增大。该u值特性受增长项和衰减项的共同影响,与相交波强度相关。

  • 图  1  RCS的几何尺寸及配筋方式

    Figure  1.  Dimensions of RCS and reinforcement layout

    图  2  SCS的几何尺寸和结构形式

    Figure  2.  Dimensions and structural style of SCS

    图  3  实验装置

    Figure  3.  Experimental setup

    图  4  测点布置

    Figure  4.  Arrangement of measure points

    图  5  数值模拟模型

    Figure  5.  Numerical model

    图  6  收敛性分析

    Figure  6.  Convergence analysis

    图  7  RCS试件破坏的实验结果

    Figure  7.  Experimental results of RCS damage

    图  8  RCS试件破坏的实验和数值模拟结果

    Figure  8.  Experimental and numerical results of RCS damages

    图  9  RCS试件钢筋变形的实验和数值模拟结果

    Figure  9.  Experimental and numerical results of RCS’s rebar deformation

    图  10  RCS试件钢筋的最大挠度

    Figure  10.  Maximum deflection of RCS’s rebar

    图  11  RCS试件测点的位移曲线

    Figure  11.  Displacement curves of RCS

    图  12  RCS试件测点的加速度曲线

    Figure  12.  Acceleration curves of RCS

    图  13  SCS试件破坏的实验结果

    Figure  13.  Experimental results of SCS damage

    图  14  SCS试件的栓钉拔出

    Figure  14.  Studs pull out of SCS

    图  15  SCS试件钢板的跨中挠度

    Figure  15.  Deflection in midspan of SCS’s steel plate

    图  16  SCS试件破坏的实验和数值模拟结果

    Figure  16.  Experimental and numerical results of SCS damages

    图  17  SCS试件核心混凝土的数值模拟结果

    Figure  17.  Numerical results of SCS’s concrete core

    图  18  SCS试件测点的位移曲线

    Figure  18.  Displacement curves of SCS

    图  19  SCS试件测点的加速度曲线

    Figure  19.  Acceleration curves of SCS

    图  20  不同炸药量时SCS中混凝土的有效塑性应变

    Figure  20.  Effective plastic strains of concrete in SCS with different explosive charges

    图  21  不同炸药量时SCS的跨中位移曲线

    Figure  21.  Mid-span displacement curves of SCS with different explosive charges

    图  22  不同炸药量时SCS的跨中最大位移

    Figure  22.  Maximum displacements of SCS with different explosive charges

    图  23  不同混凝土强度时SCS混凝土的有效塑性应变

    Figure  23.  Effective plastic strains of concrete in SCS with different concrete strengths

    图  24  不同混凝土时SCS的跨中位移曲线

    Figure  24.  Mid-span displacement curves of SCS with different concrete strengths

    图  25  不同混凝土时SCS的跨中最大位移

    Figure  25.  Maximum displacement of SCS with different concrete strengths

    图  26  不同钢板厚度时SCS混凝土的有效塑性应变

    Figure  26.  Effective plastic strains of concrete in SCS with different thickness of steel plate

    图  27  不同钢板厚度时SCS的跨中位移曲线

    Figure  27.  Mid-span displacement curves of SCS with different thicknesses of steel plates

    图  28  不同钢板厚度时SCS的跨中最大位移

    Figure  28.  Maximum displacements of SCS with different thicknesses of steel plates

    图  29  SCS跨中挠度与炸药量、钢板厚度的关系

    Figure  29.  Mid-span deflections of SCS versus explosive charges and thicknesses of steel plates

    表  1  材料力学性能参数

    Table  1.   Mechanical properties of materials

    材料型号弹性模量/GPa抗压强度/MPa屈服强度/MPa抗拉强度/MPa
    混凝土C30 3030
    钢筋HRB335200341472
    钢板Q235200235370
    焊钉A2-50200210500
    下载: 导出CSV

    表  2  试件的损伤和挠度

    Table  2.   Damages and deflections of specimens

    试件迎爆面爆坑尺寸/mm背爆面爆坑尺寸/mm实验跨中挠度/mm数值跨中挠度/mm是否发生贯穿破坏整体性是否能继续承载
    RCS360×300410×4005046.2一般
    SCS280×1803527.4
    下载: 导出CSV

    表  3  混凝土-钢-混凝土组合板跨中挠度经验公式拟合结果

    Table  3.   Fitting results of empirical formula for mid-span deflection of SCS

    炸药量w/g钢板厚度t/mm实际挠度γ0/mm预测挠度γ/mm误差/%
    1003.010.9010.623 142.54
    1503.012.0813.136 228.74
    2003.014.8813.372 4810.13
    2503.015.0516.003 816.34
    3003.027.1526.762 191.43
    3002.032.6232.594 100.08
    3002.529.1329.235 850.36
    3003.524.7624.870 360.45
    3004.023.3123.281 850.12
    下载: 导出CSV
  • [1] 王威, 张龙旭, 苏三庆, 等. 波形钢板剪力墙抗震性能试验研究 [J]. 建筑结构学报, 2018, 39(5): 36–44. DOI: 10.14006/j.jzjgxb.2018.05.005.

    WANG W, ZHANG L X, SU S Q, et al. Experimental research on seismic behavior of corrugated steel plate shear wall [J]. Journal of Building Structures, 2018, 39(5): 36–44. DOI: 10.14006/j.jzjgxb.2018.05.005.
    [2] 聂建国, 樊健生, 黄远, 等. 钢板剪力墙的试验研究 [J]. 建筑结构学报, 2010, 31(9): 1–8. DOI: 10.14006/j.jzjgxb.2010.09.015.

    NIE J G, FAN J S, HUANG Y, et al. Experimental research on steel plate shear wall [J]. Journal of Building Structures, 2010, 31(9): 1–8. DOI: 10.14006/j.jzjgxb.2010.09.015.
    [3] ZHAO C F, LU X, WANG Q, et al. Experimental and numerical investigation of steel-concrete (SC) slabs under contact blast loading [J]. Engineering Structures, 2019, 196: 109337. DOI: 10.1016/j.engstruct.2019.109337.
    [4] ZHAO C F, WANG Q, LU X, et al. Blast resistance of small-scale RCS in experimental test and numerical analysis [J]. Engineering Structures, 2019, 199: 109610. DOI: 10.1016/j.engstruct.2019.109610.
    [5] SOHEL K M A, LIEW J Y R. Behavior of steel-concrete-steel sandwich slabs subject to impact load [J]. Journal of Constructional Steel Research, 2014, 100: 163–175. DOI: 10.1016/j.jcsr.2014.04.018.
    [6] LIEW J Y R, WANG T Y. Novel steel-concrete-steel sandwich composite plates subject to impact and blast load [J]. Advances in Structural Engineering, 2011, 14(4): 673–687. DOI: 10.1260/1369-4332.14.4.673.
    [7] ZHAO C F, CHEN J Y. Damage mechanism and mode of square reinforced concrete slab subjected to blast loading [J]. Theoretical and Applied Fracture Mechanics, 2013, 63: 54–62. DOI: 10.1016/j.tafmec.2013.03.006.
    [8] ZHAO C F, WANG Q, LU X, et al. Numerical study on dynamic behaviors of NRC slabs in containment dome subjected to close-in blast loading [J]. Thin-Walled Structures, 2019, 135: 269–284. DOI: 10.1016/j.tws.2018.11.013.
    [9] 赵春风, 王强, 王静峰, 等. 近场爆炸作用下核电厂安全壳穹顶钢筋混凝土板的抗爆性能 [J]. 高压物理学报, 2019, 33(2): 025101. DOI: 10.11858/gywlxb.20180598.

    ZHAO C F, WANG Q, WANG J F, et al. Blast resistance of containment dome reinforced concrete slab in NPP under close-in explosion [J]. Chinese Journal of High Pressure Physics, 2019, 33(2): 025101. DOI: 10.11858/gywlxb.20180598.
    [10] 赵春风, 卢欣, 何凯城, 等. 单钢板混凝土剪力墙抗爆性能研究 [J]. 爆炸与冲击, 2020, 40(12): 121403. DOI: 10.11883/bzycj-2020-0058.

    ZHAO C F, LU X, HE K C, et al. Blast resistance property of concrete shear wall with single-side steel plate [J]. Explosion and Shock Waves, 2020, 40(12): 121403. DOI: 10.11883/bzycj-2020-0058.
    [11] 汪维, 杨建超, 汪剑辉, 等. POZD涂层方形钢筋混凝土板抗接触爆炸试验研究 [J]. 爆炸与冲击, 2020, 40(12): 121402. DOI: 10.11883/bzycj-2020-0180.

    WANG W, YANG J C, WANG J H, et al. Experimental research on anti-contact explosion of POZD coated square reinforced concrete slab [J]. Explosion and Shock Waves, 2020, 40(12): 121402. DOI: 10.11883/bzycj-2020-0180.
    [12] YAN C, WANG Y H, ZHAI X M, et al. Low velocity impact performance of curved steel-concrete-steel sandwich shells with bolt connectors [J]. Thin-Walled Structures, 2020, 150: 106672. DOI: 10.1016/j.tws.2020.106672.
    [13] 中华人民共和国住房和城乡建设部. 钢板剪力墙技术规程: JGJ/T 380–2015 [S]. 北京: 中国建筑工业出版社, 2016.
    [14] HALLQUIST J O. LS-DYNA keyword user’s manual [Z]. Livermore: Livermore Software Technology Corporation, 2007.
    [15] MALVAR L J, CRAWFORD J E, MORILL K B. K&C concrete material model release Ⅲ: automated generation of material model input: Technical Report TR-99-24.3 [R]. Glendale: Karagozian and Case Structural Engineers, 2000.
    [16] BISCHOFF P H, PERRY S H. Compressive behaviour of concrete at high strain rates [J]. Materials and Structures, 1991, 24(6): 425–450. DOI: 10.1007/BF02472016.
  • 期刊类型引用(0)

    其他类型引用(1)

  • 加载中
图(29) / 表(3)
计量
  • 文章访问数:  707
  • HTML全文浏览量:  366
  • PDF下载量:  122
  • 被引次数: 1
出版历程
  • 收稿日期:  2020-08-24
  • 修回日期:  2020-12-30
  • 网络出版日期:  2021-08-20
  • 刊出日期:  2021-09-14

目录

/

返回文章
返回