基于小波包技术的爆破地震波特征提取及预报
doi: 10.11883/1001-1455(2004)03-0261-7
-
摘要: 在微风化花岗岩石场地进行了孔径为76mm、孔深为5m的单孔和单段多孔爆破实验研究,获得了该类场地爆破地震波加速度衰减规律,加速度衰减系数为k=3698,=2.046,并研究了该场地基频特征。运用了小波包技术对爆破地震波测试信号特征量提取,分析了试验所测爆破地震波不同频带下小波包系数的衰减规律。在此基础上建立了基于不同频带小波包系数的爆破地震波预报模型,该模型在爆破地震波波形、峰值及频率预报方面均具有较高的准确度。
-
钢筋混凝土因其优异的性能在军事和民用建筑领域得到了广泛应用。在现代战争中,为了更好地摧毁目标或发挥防护功能,对于钢筋混凝土在侵彻打击下动态响应的研究具有十分重要的意义。
目前,国内外学者对于弹体侵彻素混凝土的研究已相对成熟。钢筋混凝土方面,Young[1-2]在经历了数千次实验后,提出了弹体侵彻岩石、土、混凝土、钢筋混凝土的统一经验公式,可较好地预测侵彻深度,但受限于经验模型本身的特点,无法给出侵彻阻力等中间变量和运动参量随时间变化的关系。Luk等[3]、Forrestal等[4]认为侵彻过程中钢筋只起到阻碍混凝土裂纹开展的作用,忽略了其提供的直接阻力,并基于此建立了侵彻模型,但其计算结果与实验数据相差20%以上。黄民荣等[5]通过实验发现钢筋的破坏模式有两种:弯曲拉伸破坏和弯曲剪切破坏,并据此建立了两种失效模式下弹体侵彻钢筋混凝土的阻力模型,然而只适用于弹体低速撞击靶体,导致模型应用受限。欧阳春等[6]、周宁等[7]、穆朝民等[8]在空腔膨胀理论的基础上,考虑了钢筋对弹体侵彻提供的直接阻力,认为钢筋发生弯曲变形破坏,建立了阻力模型;刘志林等[9]认为钢筋发生弯曲拉伸耦合断裂失效,分析其对弹体的阻力作用。此外,他们注意到侵彻后弹体接触到的钢筋发生断裂,未接触到的钢筋则完好无损。
大量研究表明,根据不同速度条件下弹体响应的机制,可将侵彻划分为4个阶段[10],即:(1)刚性弹侵彻阶段(侵彻速度小于800 m/s);(2)变形弹侵彻阶段(侵彻速度介于800~1 500 m/s);(3)半流体侵彻阶段(侵彻速度介于1 500 ~3 000 m/s);(4)流体侵彻阶段(侵彻速度大于3 000 m/s)。本文中将研究撞击速度为800 m/s以下的刚性弹侵彻问题,并从钢筋弯曲剪切破坏方面分析,且基于以下假设:(1)钢筋与混凝土的粘结为理想粘结;(2)忽略钢筋与弹体之间、钢筋与混凝土之间的摩擦力;(3)弹体与钢筋接触时,只在接触点发生相互作用;(4)弹体侵彻钢筋混凝土靶体时,受力均匀且弹道不发生偏转。以钢筋网眼中心、钢筋交汇处、钢筋网格一边中点等3个特殊位置进行研究,建立弹体侵彻钢筋混凝土靶体的阻力模型,将计算结果与Young公式进行对比,结果吻合良好,表明本文中提出的理论模型能够较为合理地预测侵彻深度及其他侵彻过程量并给出弹体各运动参量随时间变化关系,进一步得到最不利的撞击位置及在此条件下的侵彻深度计算公式。
1. 刚性弹体侵彻钢筋混凝土的理论模型
弹体侵彻钢筋混凝土厚靶的弹道轨迹包含了两部分,即在入靶面形成的锥形前坑区域以及与弹体等直径的隧道区域。对于平头弹,Li等[11]利用滑移线理论得出前坑深度为Hrc=0.707d0,其中d0为弹体直径。
1.1 隧道阶段
1.1.1 基于刚塑性假设下单根钢筋的动态分析模型
单根钢筋的受力分析如图1所示(取一半为研究对象),A、H为塑性铰位置,
λ 为塑性铰A、H间的水平距离,q为混凝土提供给钢筋的压力,近似简化为均布荷载,Mp、Qp为钢筋塑性极限弯矩、剪力,Fz为弹体轴向侵彻阻力。设钢筋具有有限的剪切强度Qp,并忽略在极限状态时弯矩M和剪力Q之间的相互作用,则在塑性铰A处有MA=Mp且QA=Qp。
在小变形阶段,对AH段列竖直方向的动量定理:
ddt2(12ρλivi0)=Qp−qλi (1) 式中:
ρ 为钢筋线密度;λi为弹体撞击到第i层钢筋时钢筋上塑性铰A、H间的水平距离;vi0 为弹体撞击到第i层钢筋到钢筋失效间的钢筋截面A的任意时刻速度;Qp=πd21σs4√3 为钢筋极限剪力,q=d1Yc(2+π) [12]为混凝土对钢筋的阻力,d1为钢筋直径,σs为钢筋屈服应力,Yc为混凝土抗剪强度;t2为以弹体撞击到第i层钢筋为起始点的时间。弹体和与其直接接触部分钢筋的总质量
m′=m0+d0ρ (m0为弹体质量),由动量定理得:ddt2(m′vi1)=−2Qp−qd0−Fzi (2) 式中:
vi1 为弹体撞击到第i层钢筋到钢筋失效间的弹体任意时刻速度;Fzi 为弹体撞击到第i层钢筋时弹体受到的轴向侵彻阻力。将式(1)、(2)对时间积分得:
12ρλivi0=Qpt2−qλit2 (3) 2Qpt2+qd0t2+Fzit2=m′vi−m′vi1 (4) 式中:
vi 为弹体第i次碰到钢筋时的速度。AH段绕H点作转动,列动量矩定理:
ddt2(16ρλ2ivi0)=2Mp−qλ2i2 (5) 式中:为钢筋塑性极限弯矩。
式(5)对时间积分得:
16ρλ2ivi0=2Mpt2−qλ2i2t2 (6) 联立式(3)和(6)得:
qλ2i+2Qpλi−12Mp=0 (7) 则:
λi=−Qpq+√(Qpq)2+12Mpq (8) 由于塑性极限弯矩Mp和剪切强度Qp、混凝土对钢筋的阻力q都是常数,式(8)表明,在此阶段塑性铰H是一个驻定铰。
由式(3)得塑性铰A所在的钢筋端部的横向速度:
vi0=2Qpt2−2qλit2ρλi (9) 弹体的速度与塑性铰A所在的钢筋端部的横向速度并不相同,由式(4)可知:
vi1=vi−2Qpt2+qd0t2+Fzit2m′ (10) 加速度为:
˜ai=−2Qp+qd0+Fzim′ (11) 联立式(3)、(4)、(9),截面A两侧的相对速度为:
vi1−vi0=vi−(1+ρλim′)(2Qpt2−2qλit2ρλi)−2qλit2+qd0t2+Fzit2m′ (12) 将相对速度对时间积分可得A截面两侧的相对位移:
S1i=t2∫0(vi1−vi0)dt2=vit2−(1+ρλim′)(Qpt22−qλit22ρλi)−2qλit22+qd0t22+Fzit222m′ (13) 根据余同希等[13]的研究成果,S1i的极限值取为0.3d1且当钢筋剪切失效时的弹体侵彻位移:
Si=vit2−2Qpt22+qd0t22+Fzit222m′ (14) 1.1.2 多根钢筋的动态分析模型
由1.1.1节知,当弹体撞击到一根钢筋时,塑性铰H是驻定铰,则当弹体撞击到n根钢筋时,每根钢筋的塑性铰位置都为
λi 。弹体和与其直接接触部分的n根钢筋的总质量为
m=m0+n∑j=1mj (mj为与弹体接触的第j根钢筋的质量),列动量定理:ddt2(mvi1)=−2nQp−qn∑j=1bj−Fzi (15) 式中:bj为第j根钢筋与弹头部分重合的长度。
对其积分得:
2nQpt2+qt2n∑j=1bj+Fzit2=mvi−mvi1 (16) 可得,弹体的速度:
vi1=vi−2nQpt2+qt2n∑j=1bj+Fzit2m (17) 加速度:
˜ai=−2nQp+qn∑j=1bj+Fzim (18) 由式(3)、(9)、(19),每根钢筋截面A两侧的相对位移为:
S1i=vit2−(nρλim+1)(Qpt22−qλit22ρλi)−qt22n∑j=1bj+Fzit22+2nqλit222m (19) 弹体头部接触钢筋,且钢筋剪切失效时的弹体侵彻位移:
Si=vit2−2nQpt22+qt22n∑j=1bj+Fzit222m (20) 1.1.3 混凝土的动态分析模型
Forrestal等[4]利用空腔膨胀理论得出弹体在隧道区受到混凝土的阻力为:
Fzi=πd204(Sf′c+Nρ0vi(t)2) (21) 式中:
S=82.6f′−0.544c ,f′c 为混凝土无侧限抗压强度;N=1 为平头弹弹头形状系数,ρ0为混凝土密度,vi(t)为弹体任意时刻t的速度。任意时刻的加速度、速度、位移表达式为:
ai(t+Δt)=−Fzim0 (22) vi(t+Δt)=vi(t)−Fzim0Δt (23) si(t+Δt)=vi(t)Δt−Fzi2m0Δt2 (24) 1.2 模型算法流程
图2为弹体侵彻多层钢筋的示意图,其中,n为钢筋层数,h1为首层钢筋的埋置深度,hn为钢筋网间距。在隧道阶段,依据弹体是否与钢筋直接撞击,分为两种工况:(1)弹体撞击到钢筋,侵彻阻力包括钢筋与混凝土两部分的作用;(2)弹体未撞击到钢筋或钢筋失效,此时弹体仅受到混凝土提供的侵彻阻力。
图3为弹体侵彻半无限厚靶的算法流程图。其中,首先输入的参数包括(1)弹体参数:弹体质量、直径、长度、侵彻速度;(2)钢筋参数:钢筋密度、屈服强度、直径、首层钢筋埋置深度、钢筋网间距;(3)混凝土参数:混凝土的无侧限抗压强度、密度等。根据首层钢筋埋置深度与开坑深度的大小关系,判定弹体在进入隧道阶段时,是否撞击到钢筋,之后不断循环钢筋与混凝土阶段,直到弹体速度为零,其中在每一个阶段,取适当时间步长△t计算每个阶段的参量,不断循环,直到达到其边界条件,跳出此阶段进入下一个阶段。
2. 模型验证与分析
2.1 经典公式
Young[1-2]通过开展大量实验研究,得出了侵彻土、岩石、钢筋混凝土的统一公式:
H={0.0008KRN(M/A)0.7ln(1+2.15v2s×10−4),vs<61m/s0.000018KRN(M/A)0.7(vs−30.5),vs≥61m/s (25) 式中:M为弹体质量;A为弹体的截面面积;vs为弹体侵彻靶体的初始撞击速度;K为缩尺效应系数,当M<182 kg时,K=0.46M0.15;当M≥182 kg时,K=1。N为弹头性能系数,对于平头弹N=0.56;R为阻力系数,对混凝土R =0.085
Kc(11−P)(tchc)−0.06(35/σc)0.3 ,其中σc为实验时混凝土的无侧限抗压强度;Kc=(F/W1)0.3 ,W1为靶体宽度与弹体直径的比值,F=20,如果W1>F,Kc=1;P为混凝土中按体积计算的含钢百分率;tc为混凝土的凝固时间;hc为混凝土目标的厚度。在没有足够的数据无法计算R值时,采用R=0.9。由于缺乏大口径平头弹侵彻钢筋混凝土厚靶的公开实验数据,本文中将不同弹体直径和质量(80、90、100 kg作为工况1~3)、不同撞击速度(600、410、780 m/s作为工况3~5)等5种工况的计算结果与Young公式的计算结果对比来验证模型的有效性。选取的靶体的参数:靶体尺寸为2 m×2 m×2 m,无侧限抗压强度fc'=35 MPa,钢筋密度为7 850 kg/m3、钢筋屈服强度为400 MPa,钢筋直径d1=8 mm,网眼尺寸60 mm,首层钢筋埋置深度为200 mm,体积配筋率2%,钢筋层间距为160 mm;弹体参数见表1。
表 1 计算侵彻深度Table 1. Calculating penetration depth工况 弹体参数 弹体与钢筋的相对位置 侵彻深度/mm 误差/% 直径/mm 长度/mm 质量/kg 撞击速度/(m·s−1) 1 80 532 13 600 钢筋交汇处 649 0.6 钢筋网格一边中点 649 0.6 钢筋网眼中心 649 0.6 Young公式 653 − 2 90 599 18 600 钢筋交汇处 712 2.47 钢筋网格一边中点 711 2.6 钢筋网眼中心 711 2.6 Young公式 730 − 3 100 666 25 600 钢筋交汇处 801 3.84 钢筋网格一边中点 800 3.96 钢筋网眼中心 800 3.96 Young公式 833 − 4 100 666 25 410 钢筋交汇处 513 7.57 钢筋网格一边中点 513 7.57 钢筋网眼中心 512 7.75 Young公式 555 − 5 100 666 25 780 钢筋交汇处 1 046 4.56 钢筋网格一边中点 1 045 4.65 钢筋网眼中心 1 045 4.65 Young公式 1 096 − 此外,由于弹体撞击到钢筋的位置具有随机性,本文中,如图4所示,以3种典型位置的计算结果与Young公式进行比较,比较结果见表1。表1中计算结果与Young公式吻合良好,进一步验证了模型的可靠性。
2.2 分析不同着靶点
2.1节的分析表明:当钢筋直径较小时,不同着靶位置对于侵彻深度的影响较小。本节中,将弹靶参数改为:弹体长度为666 mm,质量m0=25 kg,弹体撞击速度为304 m/s;靶体尺寸为2 m×2 m×2 m,无侧限抗压强度
f′c =35 MPa,钢筋直径d1=30 mm,网眼尺寸60 mm,首层钢筋埋置深度100 mm,体积配筋率2%,钢筋层间距为50 mm。并分别对弹体直径为100 mm和140 mm进行分析,着重研究不同着靶点对于侵彻深度的影响,计算结果见图5~6。从图5中可以看出,弹体撞击到钢筋交汇处侵彻深度最大;从图6中可以看出,弹体撞击到钢筋网眼中心侵彻深度最大。2.3 平头弹侵彻钢筯混凝土靶的工程计算公式
根据防护角度的最不利位置,平头弹侵彻钢筋混凝土厚靶的工程计算公式推导如下:
由2.2知,弹体直径与网眼尺寸的比值不同,最大侵彻深度的着靶点不同。由于无法预估撞击到的钢筋层数,且钢筋对于弹体的阻滞作用较小,因此,通过计算弹体侵彻素混凝土靶体的侵彻深度,可根据钢筋层间距估算弹体撞击的钢筋层数和本文提出的解析模型,得到最终的侵彻深度。
以在侵彻阶段首次撞击到钢筋时为计算起始点,此时速度为v1。
(1)根据Forrestal等[4]计算侵彻素混凝土靶体的侵彻深度s:
s=2m0πd20ρ20ln(1+ρ0Sf′cv2)+0.707d0 式中:
v2=(2m0v20−πd30Sf′c)/(2m0+πd30ρ0) 。(2)估算弹体撞击的整钢筋层间距的层数x:
x=[shi] ,其中[]为取整函数,hi为钢筋层间距;(3)得到弹体撞击到x+1层钢筋时的速度vx+1。
①由混凝土阻力引起速度的减小,从弹体首次撞击到钢筋到第x层钢筋,依据动能定理得:
−xhi∫0Fzdx=vx+1∫v112m0dv2,v2x+1=e−2m0xhiB(A+Bv21)−AB 式中:
Fz=πd204(Sf′c+Nρ0v2)=A+Bv2 ,A=πd204Sf′c ,B=πd204Nρ0 。②由钢筋阻力引起的速度减小量:
Δv2=2(2nQp+qn∑j=1bj)m0×0.3d1x ,则:第x+1层钢筋时弹体速度为:
v2x+1=e−2m0xhiB(A+Bv21)−AB−2(2nQp+qn∑j=1bj)m0×0.3d1x (4)隧道阶段的侵彻深度
S′1 为:{S′1=xhi−v2x+12˜ax+1,若v2x+1>0且Sx+1=−v2x+12˜ax+1≤0.3d1S′1=xhi+0.3d1−m02BlnAA+B(2˜ax+1×0.3d1+v2x+1),若v2x+1>0且Sx+1=−v2x+12˜ax+1>0.3d1 当
v2x+1<0 时,此时弹体撞击到第x层,因此将式(2)中的x取x=x−1带入式(3)、(4)中计算。当弹体直径小于钢筋网眼尺寸时,弹体撞击到钢筋网眼中心时只受到混凝土提供的阻力,相比其他情况,可知此时受到的阻力越小,侵彻深度越大,越为不利。
综上所述,总侵彻深度L的公式为:
{L=−m02BlnAA+Bv20+0.707d0,若k<1L=xhi−v2x+12˜ax+1+h1=xhi+Sx+1+h1,若k>1且Sx+1=−v2x+12˜ax+1≤0.3d1L=xhi+0.3d1−m02BlnAA+BSx+1(2˜ax+1×0.3d1/Sx+1−2˜ax+1)+h1若k>1,且Sx+1=−v2x+12˜ax+1>0.3d1 式中:k为弹体直径与钢筋网眼尺寸的比值,h1为首层钢筋埋置深度。
v2x+1=e−2m0xhiB(A+Bv21)−AB−2(2nQp+qn∑j=1bj)m0×0.3d1×x; v21=e−2m0(h1−0.707d0)B(A+Bv20)−AB,x=[shi],若v2x+1<0,则x=[shi]−1; ˜ax+1=−2nQp+qn∑j=1bj+Fzx+1m0; n∑j=1bj=min{∑bn1,∑bn2,∑bn3} ,n取与n∑j=1bj 对应的ni值,对应如下:钢筋交汇处:
n1=4[k2]+2 ,∑bn1=(8[k2]∑i=1√14−i2k2+2)d0 ,钢筋网眼中心:
n2=4[k2−12]+4 ,∑bn2=8n24∑i=1√14−(i−12)21k2d0 ,一根钢筋中点:
n3=2[k2]+2[k2−12]+3 ,∑bn3=(4[k2]∑i=1√14−i2k2+4[k2−12]+1∑i=1√14−(i−12)21k2+1)d0 。3. 结 论
本文基于弹体侵彻素混凝土的经典理论,将钢筋的失效模式考虑为弯曲剪切破坏后,建立平头弹侵彻钢筋混凝土的阻力模型。将计算所得的结果与经典Young公式进行对比,结果吻合良好,表明了模型的可靠性,并能够较为合理地预测侵彻深度及其它侵彻过程量并给出弹体各运动参量随时间变化关系,通过分析不同着靶点位置对于侵彻深度的影响,发现当弹体直径小于钢筋网眼尺寸时,弹体撞击到钢筋网眼中心时,弹体受到的阻力较小,侵彻深度越大;当弹体直径大于钢筋网眼尺寸时,最不利着靶点位置视弹体直径与网眼尺寸的比值而定。最后,根据防护角度的最不利位置提出了平头弹侵彻钢筋混凝土厚靶的工程计算公式。
期刊类型引用(0)
其他类型引用(1)
-
计量
- 文章访问数: 2275
- HTML全文浏览量: 147
- PDF下载量: 62
- 被引次数: 1