A study of the design and calculation method of double-skin steel-concrete shield based on energy approach
-
摘要: 为了深入探讨双钢板混凝土结构抗强冲击时后附钢板的耗能方式变化,基于最小耗能原理和无量纲化分析,反推得到了不同材料和结构尺寸下后附钢板耗能状态。结合混凝土强度极限条件,建立了兼顾塑性变形和侵彻贯穿破坏形态的后附钢板综合耗能计算公式以及后附钢板最小临界厚度解析表达式。计算结果表明,薄钢板综合耗能可达仅考虑贯穿效应耗能的4~5倍。基于能量守恒原理,提出了双钢板混凝土遮弹层防贯穿设计六步法,给出了弹体临界贯穿速度和弹体余速计算公式。该计算公式与已有钢板混凝土结构侵彻试验结果吻合较好。Abstract: The calculation method of double-skin steel-concrete shield based on the energy method is discussed. Under the premise of a rigid projectile, the main structures that affect the energy dissipation of projectile penetration are front and rear steel plates, internal concrete, tied bars and distribution layer. In this paper, the energy consumption of steel plate in elastic state, plastic state and plastic membrane state is analyzed. Based on the principle of minimum energy consumption and dimensionless analysis, the energy consumption states of the rear steel plate of different materials and structural dimensions are inversely deduced. Combined with the strength limit condition of concrete, the formula for calculating the comprehensive energy consumption of the rear steel plate considering the plastic deformation and the penetration failure mode is proposed, and the analytical expression for the minimum critical thickness of the rear steel plate is given. The calculation results indicate that the comprehensive energy consumption of the sheet steel plate can reach 4−5 times that of merely considering the penetration effect. In view of the restraint effect of steel plates on both sides of concrete, the current relatively mature formula of concrete penetration energy consumption is modified. Based on the principle of energy conservation, a six-step method for designing double-skin steel-concrete shield is proposed, and the formula for calculating the critical penetration velocity of the projectile is given. The theoretical results are then compared with the existing test data. It indicates that the calculation formulas proposed in this paper are in good agreement with the test results and can provide a scientific guideline for the design of double-skin steel-concrete shield. The main principle of the six-step method is to take a full account of the plastic deformation energy consumption of the rear steel plate and the anti-penetration energy consumption of the front steel plate. This method is beneficial to reduce the thickness of the double-skin steel-concrete shield and improve the protection effect.
-
遮弹层是防护工程的重要组成部分,发挥着抵抗弹体侵彻贯穿的关键作用。为提高防护效能,遮弹技术研究方向主要包括提高材料性能增强侵彻阻力、使用组合结构以发挥材料性能、利用强度不对称迫使弹体偏转等[1-2]。双钢板混凝土结构通常由前/后钢板、混凝土、对拉钢筋或销钉组成,受力性能优良且施工便捷,在核电站安全壳和超高层建筑中应用广泛[3]。由于两侧钢板与混凝土相互约束,可以很好地发挥钢板抗侵彻性能,降低混凝土开坑和震塌效应,在遮弹层结构中也有广阔的应用前景。
在双钢板混凝土结构抗冲击性能研究方面,Remennikov等[4-5]、Zhao等[6]和安国青等[7]对双钢板混凝土板进行了落锤撞击试验和数值模拟研究,结果表明:栓钉和对拉钢筋等轴向约束可以有效发挥钢板膜力效应,提高钢板-混凝土的整体承载能力。在弹体侵彻研究方面,Hashimoto等[8]、Walter等[9]、Tsubota等[10]、Sugano等[11]、Barr等[12]、Mizuno等[13]、Abdel-Kader等[14]、Kim等[15]和Wu等[16]开展了一系列钢板混凝土、钢筋混凝土抗侵彻撞击性能对比试验,弹速涵盖100~730 m/s,获得了结构毁伤状态,以及弹体侵彻深度、剩余速度等宝贵的试验数据。
对于双钢板混凝土结构设计研究,主要方法为等效转换法和能量法。Abdel等[17]通过刚性弹侵彻素混凝土和低碳钢靶板的对比试验,提出了低碳钢板穿孔阻力与等效混凝土厚度计算公式。程毅等[18]通过将复合靶中钢板等效转换为混凝土,推导出了弹体剩余速度公式,并分析了不同转换公式预测方法的精确度。Bruhl等[19]提出的钢板混凝土结构设计方法主要包括3个步骤:(1)计算钢筋混凝土防贯穿临界厚度并取70%;(2)计算弹体-冲塞块剩余动能;(3)计算可以耗散剩余动能的后附钢板厚度。他们将计算结果与130个钢板混凝土靶抗侵彻试验结果进行了对比,发现该方法偏于保守,特别是试验中薄钢板的防护效能远超公式预测结果。王菲等[20]分析了钢板混凝土结构中各部分的耗能,给出了带对拉钢筋的双钢板混凝土墙防贯穿实用计算公式。上述研究中,钢板等效转换和能量耗散公式基于穿孔破坏形态[21]和空腔膨胀理论,与钢板所在位置无关。但通过观察试验现象可知,迎弹面钢板破坏形态与上述假设较吻合,但后附钢板部分为贯穿破坏,部分则发生了大挠度变形转变为拉伸撕裂破坏,两者破坏形态不同(见图1),因此需要进一步分析两者耗能区别。王明洋等[22-23]综合考虑了钢板-钢纤维混凝土遮弹层在接触爆炸时,钢纤维抗拔作用、钢板薄膜力作用和土介质反力等对能量耗散的影响,建立了钢板-钢纤维混凝土遮弹板接触爆炸下的极限设计分析方法。该方法适用于后附钢板较薄的工况,但当后附钢板达厘米级时,若仍按薄膜进行侵彻阻力分析则会产生较大误差。
本文中,基于最小耗能原理并结合无量纲分析,研究后附钢板材料参数和尺寸变化带来的耗能方式转变,建立兼顾塑性变形和侵彻贯穿破坏的后附钢板耗能计算公式。在此基础上,综合考虑迎弹面钢板、混凝土和分散层耗能,提出双钢板混凝土遮弹层防刚性弹侵彻贯穿的设计方法,给出计算公式,并与已有侵彻试验结果进行对比,验证其合理性。
1. 模型分析
弹体侵彻双钢板混凝土遮弹层的机理如图2所示,图中,
D 为弹体直径,Hsf 和Hsr 分别为迎弹面钢板和后附钢板厚度,Hc 为遮弹层中混凝土厚度,hcf 为混凝土前坑深度,hct 为弹体侵彻混凝土隧道深度,hcr 、Rcr 和φcr 分别为混凝土冲塞块高度、底面半径、锥体母线与侵彻方向的夹角,db 为拉筋或销钉间距。遮弹层下设分散层,当后附钢板较薄并产生大挠度变形时,分散层局部反力造成的耗能也应考虑。由于弹体设定为刚性,因此在侵彻双钢板混凝土遮弹层时,能量耗散主要包括迎弹面钢板耗能
Wsf 、混凝土耗能Wc 、拉筋或销钉耗能Wbar 、后附钢板耗能Wsr 和分散层耗能Wd 。遮弹层结构通常较厚,主要以销钉作为钢板与混凝土的加强连接方式,由于其长度较短[7],在抵抗冲塞块形成过程中耗能有限,保守计算时可以忽略。在弹体刚好贯穿遮弹层的临界情况下,弹体初始动能
E0 与遮弹结构耗散能量应相等,即:E0=Wsf+Wc+Wsr+Wd (1) 式中:
E0=mv20/2 ,m 为弹体质量,v0 为弹体着靶速度。本文公式中的各参数均采用国际单位制(质量:kg,长度:m,速度:m/s,材料强度和模量:Pa,等)。
遮弹层抗弹体侵彻过程中,钢板、混凝土等各个组成部分共同作用、相互影响。为简便分析,将各部分解耦进行单独研究。其中,迎弹面钢板耗能
Wsf ,混凝土耗能Wc ,以及分散层耗能Wd 等部分均采用较成熟理论,后附钢板耗能在1.1节中单独分析。1.1 迎弹面钢板耗能
参照C-L模型[24],通过弹体临界贯穿靶板速度公式,计算可得弹体贯穿迎弹面钢板时消耗的能量为:
Wsf=mN1σs2B0N2ρs[exp(πD2HsfB0N2ρs2m)−1] (2) 式中:
σs 为钢材的空腔膨胀准静态强度,ρs 为钢板密度,N1 和N2 为弹体头部几何结构参数;B0 为材料参数,根据文献[25-26],其取值范围变化较小,对于钢材可近似取为1.1。σs 可按下式计算[26]:σs=Ys√3[1+(Es√3Ys)n∫√3Es−2(1+υs)Ys√3Es0(−lnx)n1−xdx] (3) 式中:
Ys 、Es 和υs 分别为钢板屈服强度、弹性模量和泊松比,n 为钢材应变强化指数。N1 和N2 可通过下式计算[27]:N1={1+μmψ2(2π−4ϕo−2sin2ϕo)卵形弹1+μmcotα2锥形弹1+2μmψ2(2ϕb−sin2ϕb)钝头弹 (4) N2={8ψ−124ψ2+μmψ2(π2−ϕo−2sin2ϕo3−sin4ϕo12)卵形弹1+μmcot(α/2)1+cot2(α/2)锥形弹1−18ψ2+μmψ2(ϕb−sin4ϕb4)钝头弹 (5) 式中:
μm 为动摩擦因数,α 为锥形弹的顶部锥角,ψ 为卵形弹和钝头弹弹头部位曲径比(caliber radius head, CRH),ϕo 和ϕb 分别为卵形弹和钝头弹临界角度,表达式为ϕo=sin−1[1−1/(2ψ)] ,ϕb=sin−1[1/(2ψ)] 。弹体贯穿钢板后剩余速度
vsf 为:vsf=√2(E0−Wsf)m (6) 1.2 混凝土耗能
混凝土耗能可分为3个部分。
(1)在侵彻隧道形成阶段,根据文献[28]平均侵彻阻力模型,计算得到该阶段混凝土耗能为:
Wct=πD24(N1σc+μBN2ρcv2cf)hct (7) 式中:
σc 为混凝土的空腔膨胀准静态强度;μ 为弹体侵彻阻力系数,可取为0.4;B 为空腔膨胀动态阻力系数,混凝土材料可取为1.0;N1 和N2 分别通过式(4)和(5)计算;ρc 为混凝土密度;vcf 为弹体贯穿混凝土前坑以后的剩余速度;hct 为混凝土侵彻隧道深度,hct=Hc−hcf−hcr 。σc=82.6fc(fc106)−0.544 (8) 式中:
fc 为混凝土圆柱体试件抗压强度,其与立方体试件抗压强度对应关系参见GB/T 50081—2019《混凝土物理力学性能试验方法标准》[29]。vcf 的计算需考虑弹体在混凝土前坑形成和隧道侵彻阶段的阻力模型[25],可得作用在弹体头部的侵彻阻力F 为:F={c1z0≤z≤hcfπD24(N1σc+μBN2ρcv2cf)hcf<z≤hcf+hct (9) 式中:
c1 为待定系数,z 为弹体侵彻深度。在弹体侵彻混凝土前坑阶段(即
0≤z≤D ),根据加速度定律并考虑初始条件,计算可得:md2zdt2=−c1z⇒{z=vsfc2sin(c2t)c1=mc22 (10) 式中:
c2 为待定系数。假设在
t1 时弹体进入隧道侵彻阶段,根据式(9)~(10),可得此时同时满足:{mc2vsfsin(c2t1)=πD24(N1σc+μBN2ρcv2cf)vsfcos(c2t1)=vcfvsfc2sin(c2t1)=hcf (11) 根据式(10)~(11),可得:
c1=πD24hcf(N1σc+μBN2ρcv2cf)=mh2cf(v2sf−v2cf) (12) 进一步求解,可得:
v2cf=4mv2sf−πD2hcfN1σc4m+πD2hcfμBN2ρc (13) (2)根据弹体侵彻阻力方程(9),可得出在该阶段混凝土消耗的能量为:
Wcf=πD28(N1σc+μBN2ρcv2cf)hcf (14) (3)由于钢板包裹作用,混凝土冲塞块难以产生碎块飞溅情况,因此仍可对弹体侵彻产生阻力,参照前坑耗能方程,可得出在该阶段混凝土消耗的能量为:
Wcr=πD28(N1σc+μBN2ρcv2cf)hcr (15) 综上分析,混凝土总耗能
Wc 为:Wc=πD28(N1σc+μBN2ρcv2cf)(2Hc−hcf−hcr) (16) 遮弹层结构中混凝土厚度通常大于5倍弹径,因此混凝土冲塞块高度可取为
hcr=2.5D [28,30]。此外,考虑到钢板对混凝土的约束作用,参考文献[31],取hcf=D 。1.3 分散层耗能
考虑到分散层下方存在主体结构的实际情况,有限厚度的分散层动力响应可近似当做准静态问题处理[22]。此时,作用在遮弹层底部的介质反力可表示为:
Nd=KArωr (17) 式中:
K 为分散层动压缩模量,Ar 和ωr 分别为反力作用的面积和挠度位移。当后附钢板满足塑性耗能条件时,其响应区挠度沿半径呈抛物线分布,如图3所示。
积分可得分散层耗能计算公式:
Wd=πK∫Rp0r[ˉω0(1−r2R2p)]2dr=16πK(ˉω0Rp)2 (18) 式中:
Rp 和ˉω0 分别为后附钢板响应区半径和钢板中心点挠度位移最大值。2. 后附钢板耗能分析
2.1 模型建立
进行4个合理简化:
(1)由于钢板厚度有限且其下部的分散层弹性模量相对很小,垂直于钢板初始中面的应力在极短时间即可经历多次反射趋近于零,因此仅对平面内径向和环向应力进行分析;
(2)在模型试验、数值仿真和实弹打击时,为消除边界效应,遮弹层平面尺寸通常较大,同时由于混凝土对钢板的约束作用,钢板塑性响应区可视为固支圆板;
(3)弹体侵彻过程中,从混凝土弹性前驱波到达后附钢板时起,至弹体和结构停止运动时止,后附钢板一直承受冲击荷载作用,但由于侵彻前期混凝土结构下部变形极小,因此仅考虑侵彻末期混凝土冲塞块对钢板的作用,且根据动量守恒,假设加载应力为圆锥形;
(4)参考文献[32-33]中理想弹塑性固支圆板在均布荷载作用下的变形方程以及径向和环向应变相等的假定。
钢板承受冲击荷载示意图如图4所示,其中
σcr 为通过混凝土冲塞块作用在钢板上的荷载,Rp 为钢板塑性响应区半径。纪冲等[34]研究发现:混凝土冲塞块破坏面倾角
φcr 与混凝土单轴拉压强度和单轴抗拉强度的比ηc=fc/ft 有关。由此可得冲塞块底面(加载区)半径Rcr 为:Rcr=λcrD=D2(1+51.2ηc−1√3.6ηc−0.75) (19) 式中:
λcr 为加载区半径与弹体直径比值。在该区域,加载应力满足下列形式:
σcr=σ0(1−rRcr) (20) 式中:
σ0 为圆心处加载应力,且随时间变化。根据文献[22]关于钢板混凝土结构破坏时的后附钢板塑性响应区域研究,可得塑性响应区半径为:
Rp=λpD={Rcr+db2=D2(1+51.2ηc−1√3.6ηc−0.75+λb)带拉筋Rcr+hcr=D2(6+51.2ηc−1√3.6ηc−0.75)无拉筋 (21) 式中:
λp 为塑性响应区半径与弹体直径比值,λb 为拉筋/销钉间距与弹体直径比值。ω=ω0(1−r2R2p) (22) 式中:
ω0 为圆板中心挠度,随时间变化。2.2 屈服准则
将后附钢板视为理想弹塑性材料且率无关,满足Tresca屈服条件。加载时钢板应力与应变关系为:
σs={Esεε≤εpYsεp<ε<ˉε0ε≥ˉε (23) 式中:
Ys 为后附钢板屈服强度,ε 为钢板中的应变,εp 为钢板弹性极限应变,ˉε 为钢板承载力极限应变。当弯矩和膜力联合作用时,屈服条件表达式[33]为:
|MθM0|+(NθN0)2=1,|MrM0|+(NrN0)2=1 (24) 式中:
Mθ 、Mr 分别为单位长度环向和径向弯矩,Nθ 、Nr 分别为单位长度环向和径向轴力,M0 为钢板单位长度上的塑性极限弯矩M0=YsH2sr/4 ,N0 为钢板单位长度上的极限膜力N0=YsHsr 。可以看出,当膜力等于N0 时,钢板进入塑性膜阶段,此时弯矩为零。2.3 耗能计算
随着加载应力增加,圆板变形增大,并从弹性、弹塑性向塑性膜状态发展,钢板中也可能同时存在弹性区、塑性区、塑性铰线区和塑性膜区[36]。但由于各区域范围随挠度变化而变化,因此难以解析求解。参考文献[36]中耗能方程表达式,基于最小耗能原理,通过不同应力状态耗能对比,反推得出一定挠度变形下的钢板可能状态。
假设钢板中心产生单位位移
δω=1 。按弹性状态计算,钢板耗能为:
ΔWe=2π∫Rp02Esω0H3sr3R4p(1−υs)rdr+2π∫Rp04EsHsr(1−υs)ω30R4p(1−r2R2p)2rdr+2πEsω0H3sr3R2p(1−υs)=4πEsω0Hsr3R2p(1−υs)(H2sr+ω20) (25) 式中:
Es 和υs 分别为后附钢板弹性模量和泊松比。按弯矩轴力联合作用下的塑性状态计算,钢板耗能为:
ΔWp=2π∫Rp0YsH2srR2p[1−(11−υsEsYs)2ω40R4p(1−r2R2p)2]rdr+2π∫Rp04EsHsr(1−υs)ω30R4p(1−r2R2p)2rdr+πYsH2sr=2πYsH2sr−πYsH2sr(11−υsEsYs)2ω403R4p+4πEsω30Hsr3R2p(1−υs) (26) 按塑性膜状态计算,钢板耗能为:
ΔWf=2π∫Rp04YsHsrω0R2p(1−r2R2p)rdr+πYsH2sr=2πYsHsr(ω0+Hsr2) (27) 通过式(25)~(27)可知,钢板弹/塑性耗能均随着挠度增大而增加,即如果要保持钢板弹/塑性变形持续发展,加载应力需要不断增大。在考虑混凝土强度的情况下,可以计算得出当中心产生单位位移
δω=1 时荷载极限做功表达式:ΔWmax=πD2fc4(1−310R2crR2p) (28) 为便于分析,设
γ=Ys/fc ,λsr=Hsr/D ,λω=ω0/D ,并对上述耗能和极限做功方程进行无量纲化处理:ΔWeYsD2δω=4πλsrλω(λ2sr+λ2ω)3εpλ2p(1−υs) (29) ΔWpYsD2δω=2πλ2sr−πλ2srλ4ω3λ4pε2p(11−υs)2+4πλ3ωλsr3εpλ2p(1−υs) (30) ΔWfYsD2δω=2πλsr(λω+λsr2) (31) ΔWmaxYsD2δω=π4γ(1−310λ2crλ2p) (32) 同时,根据钢板极限应变条件,可知
ˉω0=√ˉεRp ,无量纲化后为:ˉλω=ˉω0D=√ˉελp (33) 分别取屈服强度
Ys=400, 700 MPa 的钢材作为后附钢板材料,混凝土选择fc=32 MPa 的普通混凝土和fc=80 MPa 的高强混凝土,销钉间距与弹体直径相当,具体参数详见表1,可计算得到不同厚度和挠度变形下钢板耗能与荷载做功计算结果,见图5。表 1 钢板混凝土结构常用材料和结构参数Table 1. Common material and structural parameters of steel plate concrete structureEs/GPa Ys/MPa υs ˉε fc/MPa ρc/(kg·m−3) γ λp ˉλω 带拉筋 无拉筋 带拉筋 无拉筋 210 400 0.3 0.2 32 2400 12.50 5.63 7.63 2.51 3.41 210 400 0.3 0.2 80 2400 5.00 5.63 7.63 2.51 3.41 210 700 0.3 0.15 32 2400 21.88 5.63 7.63 2.18 2.95 210 700 0.3 0.15 80 2400 8.75 5.63 7.63 2.18 2.95 从图5可以明显看出:弹性状态耗能基本上都小于等于塑性状态耗能,但随着钢板中心挠度增大,两者均呈指数级增长,超过塑性膜状态耗能。根据最小耗能原理,在变形量相同的情况下,结构总是按耗能少的方式运行,针对本文中所述情况,即在不考虑荷载极值条件的一般情况下,理想弹塑性钢板挠度较小时可以按弹性状态变形计算,荷载做功以弹性应变能储存在钢板中,卸载时恢复;钢板挠度较大时可以按塑性膜状态变形计算,荷载做功被塑性变形消耗。
当式(29)、(31)和(32)两两分别相等时,临界挠度值
λEF 、λEC 和λFC 分别满足:4λ3EF+[4λ2sr−6εpλ2p(1−υs)]λEF−3εpλ2pλsr(1−υs)=0 (34) λ3EC+λ2srλEC−3εpλ2p(1−υs)16γλsr(1−310λ2crλ2p)=0 (35) λFC=18γλsr(1−310λ2crλ2p)−λsr2 (36) 后附钢板累积塑性耗能方程为:
Wsr1={0λFC≤λECπYsD3λsr[λ2FC−λ2EC+λsr(λFC−λEC)]λEC<λFC<ˉλωπYsD3λsr[ˉλ2ω−λ2EC+λsr(ˉλω−λEC)]λFC≥ˉλω (37) 当
λFC=ˉλω 时,可计算得出钢板最小临界厚度:ˉHsr=[√ˉλ2ω+14γ(1−310λ2crλ2p)−ˉλω]D (38) 当钢板厚度小于最小临界厚度时,钢板破坏较早,无法完全发挥混凝土防护效能,因此在结构设计时后附钢板厚度应不小于此临界厚度。此时,后附钢板还可以承受弹体的侵彻冲击,按线性差值法计算弹体侵彻贯穿时的钢板耗能方程:
Wsr2=κmN1σs2B0N2ρs[exp(πD2HsrB0N2ρs2m)−1] (39) 式中:
κ 为折减系数。κ 按下式计算:κ={1λFC≤λEF1−λ2FC−λ2EFˉλ2ω−λ2EFλEF<λFC<ˉλω0λFC≥ˉλω (40) 综上所述,后附钢板可供耗散的能量为:
Wsr=Wsr1+Wsr2 (41) 在相同条件下,将式(41)除以式(2)的比值称为耗能增强因子。仍取表1的材料参数,弹体战斗部采用BLU-109/B(
m=897.9kg ,D=368.3mm [1]),2种耗能计算公式结果对比见图6。可以看出,在BLU-109/B弹体侵彻情况下,当钢板较薄时按式计算的后附钢板耗能最高可达仅考虑侵彻贯穿耗能的4~5倍,且钢板极限应变越大时比值越高;当钢板较厚时塑性耗能可忽略,可仅按侵彻贯穿破坏形态进行耗能计算。上述结论可以较好地解释文献[19]中薄钢板防护效能被低估的结果。图6中各曲线极大值点对应的钢板厚度,即为此种情况下后附钢板的最小临界厚度。
3. 遮弹层防贯穿设计计算方法
双钢板混凝土遮弹层设计计算总体原则是钢板总厚度一定时,充分发挥后附钢板塑性耗能作用和迎弹面钢板抗弹体侵彻贯穿作用,按照弹体动能计算-迎弹面钢板厚度估算-迎弹面钢板耗能计算-后附钢板厚度和耗能计算-分散层耗能计算-混凝土厚度计算共6个步骤进行。
(1)计算目标弹体初始动能
E0 ,选择双钢板混凝土遮弹层中钢板、混凝土、拉筋/销钉等的材料。(2)选择我国规范公式[37]计算弹体对混凝土侵彻深度,选取不贯穿系数
Kg [38],计算混凝土不贯穿厚度hg ,设定折减系数α (如0.3、0.4、……),计算混凝土折减厚度hc=αhg ,依据转换公式[18]估算与此相当的迎弹面钢板厚度Hsf 。其中:参考我国规范中混凝土侵彻深度计算公式[37],结合不贯穿系数Kg 计算hg :hg=Kgλ1λ2KqmD2v0 (42) 式中:
λ1 、λ2 和Kq 分别为弹形、弹径和混凝土侵彻系数。钢板与混凝土等效厚度的转换公式[18]为:hc=2896D116f0.375c(mρc)18Hsf916 (43) (3)利用式(2)和式(6)计算迎弹面钢板耗能
Wsf 和弹体穿透该钢板时的剩余速度vsf 。(4)利用式(38)计算后附钢板临界厚度,选取市面上厚度相近的钢板为后附钢板
Hsr ,利用式(41)计算后附钢板耗能Wsr 。(5)利用式(33)或式(36)计算后附钢板极限挠度变形,依据式(18)计算分散层耗能
Wd 。(6)利用式(1)计算需要混凝土耗散的能量
Wc ,依据式(16)计算混凝土厚度Hc 。如果双钢板混凝土遮弹层各参数已经给定,则可通过上述公式反算得到弹体临界贯穿速度
vbl :v2bl=πD2(N1σc−μBN2ρc4Wsf+πD3N1σc2m+πD3N2ρc)(2Hc−hcf−hcr)+8(Wsf+Wsr+Wd)4m−2πmD2μBN2ρc2m+πD3N2ρc(2Hc−hcf−hcr) (44) 如果弹体侵彻速度
v0 大于临界贯穿速度vbl ,则弹体余速vre 可用下式计算:vre=√v20−v2bl (45) 4. 验 证
文献[14]进行了钝头弹侵彻双钢板混凝土靶试验。试验弹体和靶标如图7所示,其中弹体质量m = 175 g;前后钢板密度ρs = 7 850 kg/m3,屈服强度Ys = 240 MPa,极限应变
ˉε=0.2 ;混凝土ρc = 2 400 kg/m3,fc = 26 MPa或fc = 50 MPa;销钉直径6 mm。由于结构参数已知,本文中通过计算临界贯穿速度理论值,并与试验结果[14]进行对比,验证理论分析的合理性,计算中钢材应变强化指数n = 0.1[19],μm= 0.02[24],弹性模量和泊松比取常用参数Es = 210 GPa,υs = 0.3。靶板背部临空,所以分散层耗能为零。根据式(44)计算得到的结果与试验结果对比如表2所示,可以看出本文计算公式可以准确预测双钢板混凝土遮弹层抗侵彻性能。
表 2 文献[14]侵彻试验结果与公式计算结果对比Table 2. Comparison between penetration test results in reference [14] and calculation results in this paper试验编号 迎弹面钢板厚度/
mm后附钢板厚度/
mm混凝土抗压强度/
MPa试验弹体着靶
速度/(m·s−1)试验结果 本文公式计算的临界
贯穿速度/(m·s−1)预测是否
准确St-1-1-A 1 1 26 316 未贯穿 344 准确 St-1-1-B 1 1 26 338 未贯穿 344 准确 St-1-1-C 1 1 26 349 贯穿 344 准确 St-1-2-A 1 2 26 320 未贯穿 354 准确 St-1-2-B 1 2 26 324 未贯穿 354 准确 St-1-2-C 1 2 26 331 未贯穿 354 准确 St-1-2 1 2 50 393 未贯穿 402 准确 St-2-1-A 2 1 26 335 未贯穿 354 准确 St-2-1-B 2 1 26 350 未贯穿 354 准确 St-2-1-C 2 1 26 431 贯穿 354 准确 文献[16]进行了卵形弹侵彻钢板-钢筋混凝土靶试验。试验弹体和靶标如图8所示,其中弹体质量m = 429 g;无迎弹面钢板,后附钢板密度ρs = 7 850 kg/m3,屈服强度Ys = 240 MPa,极限应变
ˉε=0.2 ;混凝土ρc = 2 400 kg/m3,无侧限圆柱体抗压强度fc = 41 MPa,内配4层直径6 mm的钢筋网(经计算,截面配筋率为0.56%)。钢板应变强化指数等参数的取值同上,迎弹面钢板和分散层耗能为零。由于配筋对混凝土的约束作用会增强其侵彻阻力和耗能,因此参考文献[39]中fc = 41 MPa混凝土中不同配筋率对弹体极限贯穿速度的影响,对式(16)等号右边乘以1.1进行修正。本文公式计算结果与试验结果的对比如表3所示,结可以看出:弹体余速公式计算结果与试验结果相近,相对误差均在11%以内;临界贯穿速度计算结果与试验结果吻合较好。
表 3 公式计算结果与侵彻试验结果[16]的对比Table 3. Comparison of calculation results with penetration test results[16]试验
编号后附钢板
厚度/mm试验弹体着靶
速度/(m·s−1)试验
结果试验弹体余速/
(m·s−1)本文公式计算结果 临界贯穿速度/(m·s−1) 预测是否准确 弹体余速/(m·s−1) 与试验结果的误差/% 1-1 1 641.5 贯穿 272 572 准确 290 6.6 1-2 1 540 未贯穿 0 572 准确 0 0 1-3 1 601 贯穿 180 572 准确 184 2.2 1-4 1 679 贯穿 329 572 准确 365 10.9 1-5 1 737 贯穿 422 572 准确 464 10.0 5. 结 论
对双钢板混凝土遮弹层主要耗能机制进行分析,提出了防侵彻贯穿设计方法和计算公式,并与试验结果进行了对比,得出以下主要结论。
(1)随着加载应力的提高,固支理想弹塑性钢板从弹性、弹塑性向塑性膜状态发展,但钢板弹性状态耗能基本上都小于等于塑性状态耗能。在不考虑荷载极值条件的一般情况下,理想弹塑性钢板挠度较小时可以按弹性状态进行变形计算,挠度较大时可以按塑性膜状态进行耗能分析。
(2)在考虑混凝土强度情况下,计算到了后附钢板最小临界厚度。在结构设计时,后附钢板厚度应不小于此临界厚度,以更好发挥混凝土防护效能。
(3)当后附钢板较薄时,需要考虑其塑性变形的耗能作用,经分析,其综合耗能最高可达仅考虑侵彻贯穿耗能的4~5倍,且钢板极限应变越大时比值越高。
(4)基于能量守恒原理,提出的双钢板混凝土遮弹层防贯穿设计六步法,以及临界贯穿速度和弹体余速计算公式,计算结果合理可信,可充分发挥后附钢板塑性耗能作用和迎弹面钢板抗弹体侵彻贯穿作用。
-
表 1 钢板混凝土结构常用材料和结构参数
Table 1. Common material and structural parameters of steel plate concrete structure
Es/GPa Ys/MPa υs ˉε fc/MPa ρc/(kg·m−3) γ λp ˉλω 带拉筋 无拉筋 带拉筋 无拉筋 210 400 0.3 0.2 32 2400 12.50 5.63 7.63 2.51 3.41 210 400 0.3 0.2 80 2400 5.00 5.63 7.63 2.51 3.41 210 700 0.3 0.15 32 2400 21.88 5.63 7.63 2.18 2.95 210 700 0.3 0.15 80 2400 8.75 5.63 7.63 2.18 2.95 表 2 文献[14]侵彻试验结果与公式计算结果对比
Table 2. Comparison between penetration test results in reference [14] and calculation results in this paper
试验编号 迎弹面钢板厚度/
mm后附钢板厚度/
mm混凝土抗压强度/
MPa试验弹体着靶
速度/(m·s−1)试验结果 本文公式计算的临界
贯穿速度/(m·s−1)预测是否
准确St-1-1-A 1 1 26 316 未贯穿 344 准确 St-1-1-B 1 1 26 338 未贯穿 344 准确 St-1-1-C 1 1 26 349 贯穿 344 准确 St-1-2-A 1 2 26 320 未贯穿 354 准确 St-1-2-B 1 2 26 324 未贯穿 354 准确 St-1-2-C 1 2 26 331 未贯穿 354 准确 St-1-2 1 2 50 393 未贯穿 402 准确 St-2-1-A 2 1 26 335 未贯穿 354 准确 St-2-1-B 2 1 26 350 未贯穿 354 准确 St-2-1-C 2 1 26 431 贯穿 354 准确 表 3 公式计算结果与侵彻试验结果[16]的对比
Table 3. Comparison of calculation results with penetration test results[16]
试验
编号后附钢板
厚度/mm试验弹体着靶
速度/(m·s−1)试验
结果试验弹体余速/
(m·s−1)本文公式计算结果 临界贯穿速度/(m·s−1) 预测是否准确 弹体余速/(m·s−1) 与试验结果的误差/% 1-1 1 641.5 贯穿 272 572 准确 290 6.6 1-2 1 540 未贯穿 0 572 准确 0 0 1-3 1 601 贯穿 180 572 准确 184 2.2 1-4 1 679 贯穿 329 572 准确 365 10.9 1-5 1 737 贯穿 422 572 准确 464 10.0 -
[1] 任辉启, 穆朝民, 刘瑞朝, 等. 精确制导武器侵彻效应与工程防护 [M]. 北京: 科学出版社, 2016: 1–19.REN H Q, MU Z M, LIU R Z, et al. Penetration effects of precision guided weapons and engineering protection [M]. Beijing, China: Science Press, 2016: 1–19. [2] 吴应祥, 秦有权, 徐翔宇. 矢量型抗侵彻结构研究进展 [J]. 防护工程, 2022, 44(6): 36–42. DOI: 10.3969/j.issn.1674-1854.2022.06.006.WU Y X, QIN Y Q, XU X Y. Progress on vector-based anti-penetration structure research [J]. Protective Engineering, 2022, 44(6): 36–42. DOI: 10.3969/j.issn.1674-1854.2022.06.006. [3] 樊健生, 丁然, 聂鑫, 等. 高性能双钢板混凝土结构研究与应用 [J]. 建筑结构学报, 2022, 43(9): 55–72. DOI: 10.14006/j.jzjgxb.2022.0076.FAN J S, DING R, NIE X, et al. Research and application of high-performance double steel-plate reinforced concrete structures [J]. Journal of Building Structures, 2022, 43(9): 55–72. DOI: 10.14006/j.jzjgxb.2022.0076. [4] REMENNIKOV A M, KONG S Y. Numerical simulation and validation of impact response of axially-restrained steel-concrete-steel sandwich panels [J]. Composite Structures, 2012, 94(12): 3546–3555. DOI: 10.1016/j.compstruct.2012.05.011. [5] REMENNIKOV A M, KONG S Y, UY B. The response of axially restrained non-composite steel-concrete-steel sandwich panels due to large impact loading [J]. Engineering Structures, 2013, 49: 806–818. DOI: 10.1016/j.engstruct.2012.11.014. [6] ZHAO W Y, GUO Q Q, DOU X Q, et al. Impact response of steel-concrete composite panels: experiments and FE analyses [J]. Steel and Composite Structures, 2018, 26(3): 255–263. DOI: 10.12989/scs.2018.26.3.255. [7] 安国青, 王蕊, 赵晖, 等. 双钢板混凝土组合板在撞击荷载下的动力响应 [J]. 东北大学学报(自然科学版), 2022, 43(8): 1192–1200. DOI: 10.12068/j.issn.1005-3026.2022.08.017.AN G Q, WANG R, ZHAO H, et al. Dynamic response of double-skin steel-concrete composite panel under impact loading [J]. Journal of Northeastern University (Natural Science), 2022, 43(8): 1192–1200. DOI: 10.12068/j.issn.1005-3026.2022.08.017. [8] HASHIMOTO J, TAKIGUCHI K, NISHIMURA K, et al. Experimental study on behavior of RC panels covered with steel plates subjected to missile impact [C]//Proceedings of the 18th International Association for Structural Mechanics in Reactor Technology. Beijing, China: IASMiRT, 2005: 2604–2615. [9] WALTER T A, WOLDE-TINSAE A M. Turbine missile perforation of reinforced concrete [J]. Journal of Structural Engineering, 1984, 110(10): 2439–2455. DOI: 10.1061/(ASCE)0733-9445(1984)110:10(2439). [10] TSUBOTA H, KASAI Y, KOSHIKA N, et al. Quantitative studies on impact resistance of reinforced concrete panels with steel liners under impact loading. part 1: scaled model impact tests[C]//Proceedings of the 12th International Association for Structural Mechanics in Reactor Technology. Stuttgart, Germany: IASMiRT, 1993: 169–174. [11] SUGANO T, TSUBOTA H, KASAI Y, et al. Local damage to reinforced concrete structures caused by impact of aircraft engine missiles. part 1: test program, method and results [J]. Nuclear Engineering and Design, 1993, 140(3): 387–405. DOI: 10.1016/0029-5493(93)90120-X. [12] BARR P, CARTER P G, HOWE W D, et al. Experimental studies of the impact resistance of steel faced concrete composites [C]//Proceedings of the 7th International Association for Structural Mechanics in Reactor Technology. Chicago, USA: IASMiRT, 1983: 395–402. [13] MIZUNO J, KOSHIKA N, SAWAMOTO Y, et al. Investigation on impact resistance of steel plate reinforced concrete barriers against aircraft impact. part 1: test program and results [C]//Proceedings of the 18th International Association for Structural Mechanics in Reactor Technology. Beijing, China: IASMiRT, 2005: 2566–2679. [14] ABDEL-KADER M, FOUDA A. Effect of reinforcement on the response of concrete panels to impact of hard projectiles [J]. International Journal of Impact Engineering, 2014, 63: 1–17. DOI: 10.1016/j.ijimpeng.2013.07.005. [15] KIM K S, MOON I H, CHOI H J, et al. A preliminary study on the local impact behavior of steel-plate concrete walls [J]. Annals of Nuclear Energy, 2017, 102: 210–219. DOI: 10.1016/j.anucene.2016.12.006. [16] WU H, FANG Q, PENG Y, et al. Hard projectile perforation on the monolithic and segmented RC panels with a rear steel liner [J]. International Journal of Impact Engineering, 2015, 76: 232–250. DOI: 10.1016/j.ijimpeng.2014.10.010. [17] ABDEL-KADER M, FOUDA A. Equivalent concrete thickness for perforation of mild steel plates [J]. Journal of Constructional Steel Research, 2017, 135: 213–229. DOI: 10.1016/j.jcsr.2017.04.016. [18] 程毅, 刘军, 郑雨晴, 等. 钢板混凝土复合靶抗贯穿性能的理论与数值分析 [J]. 应用力学学报, 2020, 37(4): 1441–1449. DOI: 10.11776/cjam.37.04.B056.CHENG Y, LIU J, ZHENG Y Q, et al. Theoretical and numerical investigation of perforation resistance of steel-concrete composite target [J]. Chinese Journal of Applied Mechanics, 2020, 37(4): 1441–1449. DOI: 10.11776/cjam.37.04.B056. [19] BRUHL J C, VARMA A H, JOHNSON W H. Design of composite SC walls to prevent perforation from missile impact [J]. International Journal of Impact Engineering, 2015, 75: 75–87. DOI: 10.1016/j.ijimpeng.2014.07.015. [20] 王菲, 刘晶波, 韩鹏飞, 等. 核工程钢板混凝土墙防撞击贯穿实用计算方法 [J]. 爆炸与冲击, 2020, 40(10): 105101. DOI: 10.11883/bzycj-2020-0020.WANG F, LIU J B, HAN P F, et al. A practical calculation method of steel plate concrete walls to resist perforation from missile impact in nuclear engineering [J]. Explosion and Shock Waves, 2020, 40(10): 105101. DOI: 10.11883/bzycj-2020-0020. [21] BØRVIK T, CLAUSEN A H, HOPPERSTAD O S, et al. Perforation of AA5083-H116 aluminium plates with conical-nose steel projectiles: experimental study [J]. International Journal of Impact Engineering, 2004, 30(4): 367–384. DOI: 10.1016/S0734-743X(03)00072-1. [22] 王明洋, 张胜民, 国胜兵. 接触爆炸作用下钢板-钢纤维混凝土遮弹层设计方法(I) [J]. 爆炸与冲击, 2002, 22(1): 40–45. DOI: 10.3321/j.issn:1001-1455.2002.01.008.WANG M Y, ZHANG S M, GUO S B. Design method of steel and steel-fiber concrete shelter plate under contact detonation [J]. Explosion and Shock Waves, 2002, 22(1): 40–45. DOI: 10.3321/j.issn:1001-1455.2002.01.008. [23] 王明洋, 钱七虎, 赵跃堂. 接触爆炸作用下钢板-钢纤维钢筋混凝土遮弹层设计方法(II) [J]. 爆炸与冲击, 2002, 22(2): 163–168. DOI: 10.3321/j.issn:1001-1455.2002.02.012.WANG M Y, QIAN Q H, ZHAO Y T. The design method for shelter plate of steel plate and steel fiber reinforced concrete under contact detonation [J]. Explosion and Shock Waves, 2002, 22(2): 163–168. DOI: 10.3321/j.issn:1001-1455.2002.02.012. [24] CHEN X W, HUANG X L, LIANG G J. Comparative analysis of perforation models of metallic plates by rigid sharp-nosed projectiles [J]. International Journal of Impact Engineering, 2011, 38(7): 613–621. DOI: 10.1016/j.ijimpeng.2010.12.005. [25] FORRESTAL M J, ALTMAN B S, CARGILE J D, et al. An empirical equation for penetration depth of ogive-nose projectiles into concrete targets [J]. International Journal of Impact Engineering, 1994, 15(4): 395–405. DOI: 10.1016/0734-743X(94)80024-4. [26] FORRESTAL M J, WARREN T L. Perforation equations for conical and ogival nose rigid projectiles into aluminum target plates [J]. International Journal of Impact Engineering, 2009, 36(2): 220–225. DOI: 10.1016/j.ijimpeng.2008.04.005. [27] CHEN X W, LI Q M. Deep penetration of a non-deformable projectile with different geometrical characteristics [J]. International Journal of Impact Engineering, 2002, 27(6): 619–637. DOI: 10.1016/s0734-743x(02)00005-2. [28] PENG Y, WU H, FANG Q, et al. A note on the deep penetration and perforation of hard projectiles into thick targets [J]. International Journal of Impact Engineering, 2015, 85: 37–44. DOI: 10.1016/j.ijimpeng.2015.06.013. [29] 中华人民共和国住房和城乡建设部, 国家市场监督管理总局. 混凝土物理力学性能试验方法标准: GB/T 50081—2019 [S]. 北京: 中国建筑工业出版社, 2019: 145–146.Ministry of Housing and Urban-Rural Development of the People’s Republic of China, State Administration for Market Regulation. Standard for test methods of concrete physical and mechanical properties: GB/T 50081—2019 [S]. Beijing, China: China Architecture and Building Press, 2019: 145–146. [30] PENG Y, WU H, FANG Q, et al. Residual velocities of projectiles after normally perforating the thin ultra-high performance steel fiber reinforced concrete slabs [J]. International Journal of Impact Engineering, 2016, 97: 1–9. DOI: 10.1016/j.ijimpeng.2016.06.006. [31] 张爽, 武海军, 黄风雷. 弹体侵彻钢筋混凝土靶开坑深度研究 [J]. 北京理工大学学报, 2018, 38(6): 565–571. DOI: 10.15918/j.tbit1001-0645.2018.06.003.ZHANG S, WU H J, HUANG F L. Investigation on crater depth of projectile penetrating into reinforced concrete target [J]. Transactions of Beijing Institute of Technology, 2018, 38(6): 565–571. DOI: 10.15918/j.tbit1001-0645.2018.06.003. [32] KARAGIOZOVA D, YU T X, SHI S Y, et al. On the influence of elasticity on the large deflections response of circular plates to uniform quasi-static pressure [J]. International Journal of Mechanical Sciences, 2017, 131/132: 894–907. DOI: 10.1016/j.ijmecsci.2017.07.032. [33] KYOHEI K, PIAN T H H. Large deformations of rigid-plastic circular plates [J]. International Journal of Solids and Structures, 1981, 17(11): 1043–1055. DOI: 10.1016/0020-7683(81)90012-3. [34] 纪冲, 龙源, 高振儒, 等. 弹丸冲击贯穿有限厚混凝土材料靶板的背面成坑效应 [J]. 振动与冲击, 2013, 32(4): 85–89,101. DOI: 10.13465/j.cnki.jvs.2013.04.003.JI C, LONG Y, GAO Z R, et al. Rear face crater-forming of a limited-thickness concrete target due to projectile penetrating [J]. Journal of Vibration and Shock, 2013, 32(4): 85–89, 101. DOI: 10.13465/j.cnki.jvs.2013.04.003. [35] CLOETE T J, NURICK G N. On the influence of radial displacements and bending strains on the large deflections of impulsively loaded circular plates [J]. International Journal of Mechanical Sciences, 2014, 82: 140–148. DOI: 10.1016/j.ijmecsci.2014.02.026. [36] 吴迪, 米国, 郭香华, 等. 空爆载荷作用下固支弹塑性圆板的动力学模型 [J]. 高压物理学报, 2022, 36(5): 054202. DOI: 10.11858/gywlxb.20220525.WU D, MI G, GUO X H, et al. Dynamic model of clamped elastoplastic circular plate under air blast loading [J]. Chinese Journal of High Pressure Physics, 2022, 36(5): 054202. DOI: 10.11858/gywlxb.20220525. [37] 王安宝, 邓国强, 张磊, 等. 混凝土侵彻公式的合理性分析 [J]. 防护工程, 2020, 42(6): 1–7. DOI: 10.3969/j.issn.1674-1854.2020.06.001.WANG A B, DENG G Q, ZHANG L, et al. Analysis on the rationality of concrete penetration formula [J]. Protective Engineering, 2020, 42(6): 1–7. DOI: 10.3969/j.issn.1674-1854.2020.06.001. [38] 龚自明, 方秦, 张亚栋, 等. 混凝土靶体侵彻不贯穿系数的试验研究 [J]. 兵工学报, 2009, 30(9): 1181–1186. DOI: 10.3321/j.issn:1000-1093.2009.09.006.GONG Z M, FANG Q, ZHANG Y D, et al. Experimental investigation into coefficients of concrete targets to prevent perforation [J]. Acta Armamentarii, 2009, 30(9): 1181–1186. DOI: 10.3321/j.issn:1000-1093.2009.09.006. [39] 陈小伟. 穿甲/侵彻力学的理论建模与分析(上册) [M]. 北京: 科学出版社, 2019: 249–264.CHEN X W. Modelling on the perforation and penetration (I) [M]. Beijing, China: Science Press, 2019: 249–264. 期刊类型引用(1)
1. 朱擎,李述涛,陈叶青,马上,魏万里,张生,陈嘉琳. 高强钢-钢筋混凝土复合防护结构厚度极限计算方法. 力学学报. 2024(07): 2077-2090 . 百度学术
其他类型引用(1)
-