• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

基于细观混凝土模型的刚性弹体正侵彻弹道偏转分析

邓勇军 陈小伟 姚勇 杨涛

王武, 杨军, 王安宝, 周布奎, 李晓军. 基于能量法的双钢板混凝土遮弹层设计计算方法研究[J]. 爆炸与冲击, 2023, 43(11): 114203. doi: 10.11883/bzycj-2023-0086
引用本文: 邓勇军, 陈小伟, 姚勇, 杨涛. 基于细观混凝土模型的刚性弹体正侵彻弹道偏转分析[J]. 爆炸与冲击, 2017, 37(3): 377-386. doi: 10.11883/1001-1455(2017)03-0377-10
WANG Wu, YANG Jun, WANG Anbao, ZHOU Bukui, LI Xiaojun. A study of the design and calculation method of double-skin steel-concrete shield based on energy approach[J]. Explosion And Shock Waves, 2023, 43(11): 114203. doi: 10.11883/bzycj-2023-0086
Citation: Deng Yongjun, Chen Xiaowei, Yao Yong, Yang Tao. On ballistic trajectory of rigid projectile normal penetration based on a meso-scopic concrete model[J]. Explosion And Shock Waves, 2017, 37(3): 377-386. doi: 10.11883/1001-1455(2017)03-0377-10

基于细观混凝土模型的刚性弹体正侵彻弹道偏转分析

doi: 10.11883/1001-1455(2017)03-0377-10
基金项目: 

国家自然科学基金项目 11225213

国家自然科学基金项目 11390361

国家自然科学基金项目 11390362

详细信息
    作者简介:

    邓勇军(1987—), 男, 博士研究生, 助理研究员

    通讯作者:

    陈小伟, chenxiaoweintu@yahoo.com

  • 中图分类号: O382

On ballistic trajectory of rigid projectile normal penetration based on a meso-scopic concrete model

  • 摘要: 为研究混凝土细观因素对刚性弹正侵彻弹道偏转的影响,基于骨料随机投放的思路建立混凝土三维细观几何模型,分析刚性弹正侵彻过程中发生弹道偏转的原因及可能影响因素,定量讨论混凝土细观因素对弹道偏转的影响。结果表明:混凝土细观数值模型可以较好地反映弹体正侵彻过程中弹道偏转等典型物理现象,且细观参数对于弹体弹道偏转有显著影响。刚性弹正侵彻细观混凝土时,存在一个弹体直径/骨料最大粒径比的特征比值。
  • 遮弹层是防护工程的重要组成部分,发挥着抵抗弹体侵彻贯穿的关键作用。为提高防护效能,遮弹技术研究方向主要包括提高材料性能增强侵彻阻力、使用组合结构以发挥材料性能、利用强度不对称迫使弹体偏转等[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  金属板不同破坏形式
    Figure  1.  Different failure modes of metal plates

    本文中,基于最小耗能原理并结合无量纲分析,研究后附钢板材料参数和尺寸变化带来的耗能方式转变,建立兼顾塑性变形和侵彻贯穿破坏的后附钢板耗能计算公式。在此基础上,综合考虑迎弹面钢板、混凝土和分散层耗能,提出双钢板混凝土遮弹层防刚性弹侵彻贯穿的设计方法,给出计算公式,并与已有侵彻试验结果进行对比,验证其合理性。

    弹体侵彻双钢板混凝土遮弹层的机理如图2所示,图中,D为弹体直径,HsfHsr分别为迎弹面钢板和后附钢板厚度,Hc为遮弹层中混凝土厚度,hcf为混凝土前坑深度,hct为弹体侵彻混凝土隧道深度,hcrRcrφcr分别为混凝土冲塞块高度、底面半径、锥体母线与侵彻方向的夹角,db为拉筋或销钉间距。遮弹层下设分散层,当后附钢板较薄并产生大挠度变形时,分散层局部反力造成的耗能也应考虑。

    图  2  刚性弹侵彻双钢板混凝土遮弹层的机理
    Figure  2.  Mechanism of a rigid projectile perforating on a double-skin steel-concrete shield

    由于弹体设定为刚性,因此在侵彻双钢板混凝土遮弹层时,能量耗散主要包括迎弹面钢板耗能Wsf、混凝土耗能Wc、拉筋或销钉耗能Wbar、后附钢板耗能Wsr和分散层耗能Wd。遮弹层结构通常较厚,主要以销钉作为钢板与混凝土的加强连接方式,由于其长度较短[7],在抵抗冲塞块形成过程中耗能有限,保守计算时可以忽略。

    在弹体刚好贯穿遮弹层的临界情况下,弹体初始动能E0与遮弹结构耗散能量应相等,即:

    E0=Wsf+Wc+Wsr+Wd
    (1)

    式中:E0=mv20/2m为弹体质量,v0为弹体着靶速度。

    本文公式中的各参数均采用国际单位制(质量:kg,长度:m,速度:m/s,材料强度和模量:Pa,等)。

    遮弹层抗弹体侵彻过程中,钢板、混凝土等各个组成部分共同作用、相互影响。为简便分析,将各部分解耦进行单独研究。其中,迎弹面钢板耗能Wsf,混凝土耗能Wc,以及分散层耗能Wd等部分均采用较成熟理论,后附钢板耗能在1.1节中单独分析。

    参照C-L模型[24],通过弹体临界贯穿靶板速度公式,计算可得弹体贯穿迎弹面钢板时消耗的能量为:

    Wsf=mN1σs2B0N2ρs[exp(πD2HsfB0N2ρs2m)1]
    (2)

    式中:σs为钢材的空腔膨胀准静态强度,ρs为钢板密度,N1N2为弹体头部几何结构参数;B0为材料参数,根据文献[25-26],其取值范围变化较小,对于钢材可近似取为1.1。σs可按下式计算[26]

    σs=Ys3[1+(Es3Ys)n3Es2(1+υs)Ys3Es0(lnx)n1xdx]
    (3)

    式中:YsEsυs分别为钢板屈服强度、弹性模量和泊松比,n为钢材应变强化指数。

    N1N2可通过下式计算[27]

    N1={1+μmψ2(4ϕo2sin2ϕo)1+μmcotα21+2μmψ2(2ϕbsin2ϕb)
    (4)
    N2={8ψ124ψ2+μmψ2(π2ϕo2sin2ϕo3sin4ϕo12)1+μmcot(α/2)1+cot2(α/2)118ψ2+μmψ2(ϕbsin4ϕb4)
    (5)

    式中:μm为动摩擦因数,α为锥形弹的顶部锥角,ψ为卵形弹和钝头弹弹头部位曲径比(caliber radius head, CRH),ϕoϕb分别为卵形弹和钝头弹临界角度,表达式为ϕo=sin1[11/(2ψ)]ϕb=sin1[1/(2ψ)]

    弹体贯穿钢板后剩余速度vsf为:

    vsf=2(E0Wsf)m
    (6)

    混凝土耗能可分为3个部分。

    (1)在侵彻隧道形成阶段,根据文献[28]平均侵彻阻力模型,计算得到该阶段混凝土耗能为:

    Wct=πD24(N1σc+μBN2ρcv2cf)hct
    (7)

    式中:σc为混凝土的空腔膨胀准静态强度;μ为弹体侵彻阻力系数,可取为0.4;B为空腔膨胀动态阻力系数,混凝土材料可取为1.0;N1N2分别通过式(4)和(5)计算;ρc为混凝土密度;vcf为弹体贯穿混凝土前坑以后的剩余速度;hct为混凝土侵彻隧道深度,hct=Hchcfhcr

    σc可用下式[25,28]表示:

    σc=82.6fc(fc106)0.544
    (8)

    式中:fc为混凝土圆柱体试件抗压强度,其与立方体试件抗压强度对应关系参见GB/T 50081—2019《混凝土物理力学性能试验方法标准》[29]

    vcf的计算需考虑弹体在混凝土前坑形成和隧道侵彻阶段的阻力模型[25],可得作用在弹体头部的侵彻阻力F为:

    F={c1z0zhcfπD24(N1σc+μBN2ρcv2cf)hcfzhcf+hct
    (9)

    式中:c1为待定系数,z为弹体侵彻深度。

    在弹体侵彻混凝土前坑阶段(即0zD),根据加速度定律并考虑初始条件,计算可得:

    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(v2sfv2cf)
    (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)(2Hchcfhcr)
    (16)

    遮弹层结构中混凝土厚度通常大于5倍弹径,因此混凝土冲塞块高度可取为hcr=2.5D[28,30]。此外,考虑到钢板对混凝土的约束作用,参考文献[31],取hcf=D

    考虑到分散层下方存在主体结构的实际情况,有限厚度的分散层动力响应可近似当做准静态问题处理[22]。此时,作用在遮弹层底部的介质反力可表示为:

    Nd=KArωr
    (17)

    式中:K为分散层动压缩模量,Arωr分别为反力作用的面积和挠度位移。

    当后附钢板满足塑性耗能条件时,其响应区挠度沿半径呈抛物线分布,如图3所示。

    图  3  后附钢板塑性响应变形示意图
    Figure  3.  Schematic diagram of plastic responsedeformation of rear steel plate

    积分可得分散层耗能计算公式:

    Wd=πKRp0r[ˉω0(1r2R2p)]2dr=16πK(ˉω0Rp)2
    (18)

    式中:Rpˉω0分别为后附钢板响应区半径和钢板中心点挠度位移最大值。

    进行4个合理简化:

    (1)由于钢板厚度有限且其下部的分散层弹性模量相对很小,垂直于钢板初始中面的应力在极短时间即可经历多次反射趋近于零,因此仅对平面内径向和环向应力进行分析;

    (2)在模型试验、数值仿真和实弹打击时,为消除边界效应,遮弹层平面尺寸通常较大,同时由于混凝土对钢板的约束作用,钢板塑性响应区可视为固支圆板;

    (3)弹体侵彻过程中,从混凝土弹性前驱波到达后附钢板时起,至弹体和结构停止运动时止,后附钢板一直承受冲击荷载作用,但由于侵彻前期混凝土结构下部变形极小,因此仅考虑侵彻末期混凝土冲塞块对钢板的作用,且根据动量守恒,假设加载应力为圆锥形;

    (4)参考文献[32-33]中理想弹塑性固支圆板在均布荷载作用下的变形方程以及径向和环向应变相等的假定。

    钢板承受冲击荷载示意图如图4所示,其中σcr为通过混凝土冲塞块作用在钢板上的荷载,Rp为钢板塑性响应区半径。

    图  4  后附钢板承受冲击荷载示意图
    Figure  4.  Schematic diagram of rear steel plate bearing impact load

    纪冲等[34]研究发现:混凝土冲塞块破坏面倾角φcr与混凝土单轴拉压强度和单轴抗拉强度的比ηc=fc/ft有关。由此可得冲塞块底面(加载区)半径Rcr为:

    Rcr=λcrD=D2(1+51.2ηc13.6ηc0.75)
    (19)

    式中:λcr为加载区半径与弹体直径比值。

    在该区域,加载应力满足下列形式:

    σcr=σ0(1rRcr)
    (20)

    式中:σ0为圆心处加载应力,且随时间变化。

    根据文献[22]关于钢板混凝土结构破坏时的后附钢板塑性响应区域研究,可得塑性响应区半径为:

    Rp=λpD={Rcr+db2=D2(1+51.2ηc13.6ηc0.75+λb)Rcr+hcr=D2(6+51.2ηc13.6ηc0.75)
    (21)

    式中:λp为塑性响应区半径与弹体直径比值,λb为拉筋/销钉间距与弹体直径比值。

    同时,假定圆板挠度沿半径呈抛物线分布[32,35]

    ω=ω0(1r2R2p)
    (22)

    式中:ω0为圆板中心挠度,随时间变化。

    将后附钢板视为理想弹塑性材料且率无关,满足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/4N0为钢板单位长度上的极限膜力N0=YsHsr。可以看出,当膜力等于N0时,钢板进入塑性膜阶段,此时弯矩为零。

    随着加载应力增加,圆板变形增大,并从弹性、弹塑性向塑性膜状态发展,钢板中也可能同时存在弹性区、塑性区、塑性铰线区和塑性膜区[36]。但由于各区域范围随挠度变化而变化,因此难以解析求解。参考文献[36]中耗能方程表达式,基于最小耗能原理,通过不同应力状态耗能对比,反推得出一定挠度变形下的钢板可能状态。

    假设钢板中心产生单位位移δω=1

    按弹性状态计算,钢板耗能为:

    ΔWe=2πRp02Esω0H3sr3R4p(1υs)rdr+2πRp04EsHsr(1υs)ω30R4p(1r2R2p)2rdr+2πEsω0H3sr3R2p(1υs)=4πEsω0Hsr3R2p(1υs)(H2sr+ω20)
    (25)

    式中:Esυs分别为后附钢板弹性模量和泊松比。

    按弯矩轴力联合作用下的塑性状态计算,钢板耗能为:

    ΔWp=2πRp0YsH2srR2p[1(11υsEsYs)2ω40R4p(1r2R2p)2]rdr+2πRp04EsHsr(1υs)ω30R4p(1r2R2p)2rdr+πYsH2sr=2πYsH2srπYsH2sr(11υsEsYs)2ω403R4p+4πEsω30Hsr3R2p(1υs)
    (26)

    按塑性膜状态计算,钢板耗能为:

    ΔWf=2πRp04YsHsrω0R2p(1r2R2p)rdr+πYsH2sr=2πYsHsr(ω0+Hsr2)
    (27)

    通过式(25)~(27)可知,钢板弹/塑性耗能均随着挠度增大而增加,即如果要保持钢板弹/塑性变形持续发展,加载应力需要不断增大。在考虑混凝土强度的情况下,可以计算得出当中心产生单位位移δω=1时荷载极限做功表达式:

    ΔWmax=πD2fc4(1310R2crR2p)
    (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γ(1310λ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 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
    下载: 导出CSV 
    | 显示表格
    图  5  钢板厚度和挠度变化对耗能的影响
    Figure  5.  Influence of deflection and thickness change on energy consumption of steel plate

    图5可以明显看出:弹性状态耗能基本上都小于等于塑性状态耗能,但随着钢板中心挠度增大,两者均呈指数级增长,超过塑性膜状态耗能。根据最小耗能原理,在变形量相同的情况下,结构总是按耗能少的方式运行,针对本文中所述情况,即在不考虑荷载极值条件的一般情况下,理想弹塑性钢板挠度较小时可以按弹性状态变形计算,荷载做功以弹性应变能储存在钢板中,卸载时恢复;钢板挠度较大时可以按塑性膜状态变形计算,荷载做功被塑性变形消耗。

    当式(29)、(31)和(32)两两分别相等时,临界挠度值λEFλECλFC分别满足:

    4λ3EF+[4λ2sr6εpλ2p(1υs)]λEF3εpλ2pλsr(1υs)=0
    (34)
    λ3EC+λ2srλEC3εpλ2p(1υs)16γλsr(1310λ2crλ2p)=0
    (35)
    λFC=18γλsr(1310λ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γ(1310λ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.9kgD=368.3mm[1]),2种耗能计算公式结果对比见图6

    图  6  钢板塑性变形对耗能的影响
    Figure  6.  Effect of plastic deformation of steel plates on energy consumption

    可以看出,在BLU-109/B弹体侵彻情况下,当钢板较薄时按式计算的后附钢板耗能最高可达仅考虑侵彻贯穿耗能的4~5倍,且钢板极限应变越大时比值越高;当钢板较厚时塑性耗能可忽略,可仅按侵彻贯穿破坏形态进行耗能计算。上述结论可以较好地解释文献[19]中薄钢板防护效能被低估的结果。图6中各曲线极大值点对应的钢板厚度,即为此种情况下后附钢板的最小临界厚度。

    双钢板混凝土遮弹层设计计算总体原则是钢板总厚度一定时,充分发挥后附钢板塑性耗能作用和迎弹面钢板抗弹体侵彻贯穿作用,按照弹体动能计算-迎弹面钢板厚度估算-迎弹面钢板耗能计算-后附钢板厚度和耗能计算-分散层耗能计算-混凝土厚度计算共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λ2Kq分别为弹形、弹径和混凝土侵彻系数。钢板与混凝土等效厚度的转换公式[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)(2Hchcfhcr)+8(Wsf+Wsr+Wd)4m2πmD2μBN2ρc2m+πD3N2ρc(2Hchcfhcr)
    (44)

    如果弹体侵彻速度v0大于临界贯穿速度vbl,则弹体余速vre可用下式计算:

    vre=v20v2bl
    (45)

    文献[14]进行了钝头弹侵彻双钢板混凝土靶试验。试验弹体和靶标如图7所示,其中弹体质量m = 175 g;前后钢板密度ρs = 7 850 kg/m3,屈服强度Ys = 240 MPa,极限应变ˉε=0.2;混凝土ρc = 2 400 kg/m3fc = 26 MPa或fc = 50 MPa;销钉直径6 mm。

    图  7  弹体及靶标示意图[14]
    Figure  7.  Overviews of the projectile and steel-concrete target[14]

    由于结构参数已知,本文中通过计算临界贯穿速度理论值,并与试验结果[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 准确
    下载: 导出CSV 
    | 显示表格

    文献[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%)。

    图  8  弹体及靶标示意图[16]
    Figure  8.  Overview of missile and steel-concrete target[16]

    钢板应变强化指数等参数的取值同上,迎弹面钢板和分散层耗能为零。由于配筋对混凝土的约束作用会增强其侵彻阻力和耗能,因此参考文献[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
    下载: 导出CSV 
    | 显示表格

    对双钢板混凝土遮弹层主要耗能机制进行分析,提出了防侵彻贯穿设计方法和计算公式,并与试验结果进行了对比,得出以下主要结论。

    (1)随着加载应力的提高,固支理想弹塑性钢板从弹性、弹塑性向塑性膜状态发展,但钢板弹性状态耗能基本上都小于等于塑性状态耗能。在不考虑荷载极值条件的一般情况下,理想弹塑性钢板挠度较小时可以按弹性状态进行变形计算,挠度较大时可以按塑性膜状态进行耗能分析。

    (2)在考虑混凝土强度情况下,计算到了后附钢板最小临界厚度。在结构设计时,后附钢板厚度应不小于此临界厚度,以更好发挥混凝土防护效能。

    (3)当后附钢板较薄时,需要考虑其塑性变形的耗能作用,经分析,其综合耗能最高可达仅考虑侵彻贯穿耗能的4~5倍,且钢板极限应变越大时比值越高。

    (4)基于能量守恒原理,提出的双钢板混凝土遮弹层防贯穿设计六步法,以及临界贯穿速度和弹体余速计算公式,计算结果合理可信,可充分发挥后附钢板塑性耗能作用和迎弹面钢板抗弹体侵彻贯穿作用。

  • 图  1  细观模型建立流程图

    Figure  1.  Flow diagram of meso-scopic modeling

    图  2  骨料投放示意图

    Figure  2.  Schematic diagram of aggregate distribution

    图  3  细观模型各组成部分有限元剖分

    Figure  3.  Each component's FEM of meso-scopic model

    图  4  靶板迎弹面破坏模式

    Figure  4.  Failure mode of the target on the impact side

    图  5  弹体尺寸

    Figure  5.  Projectile geometry

    图  6  混凝土靶板尺寸

    Figure  6.  Concrete target geometry

    图  7  不同尺寸的有限元模型弹道图

    Figure  7.  Ballistc trajectories of finite element models with different dimensions

    图  8  均匀模型混凝土靶的侵彻过程

    Figure  8.  Penetration process of uniform concrete models

    图  9  细观建模的混凝土靶的侵彻过程

    Figure  9.  Penetration process of mesoscopic concrete models

    图  10  侵彻过程中刚性弹运动示意

    Figure  10.  Motion of a rigid projectile during penetration

    图  11  弹体偏转角的定义

    Figure  11.  Definition of the deflexion angle of a rigid projectile

    图  12  弹体侵彻过程中各参数随时间的变化曲线

    Figure  12.  Motion parameter-time curves of the projectile during penetration

    图  13  靶体尺寸及侵彻位置示意图

    Figure  13.  Target dimension and penetration location

    图  14  骨料随机分布对弹体偏转的影响

    Figure  14.  Projectile deflection affected by aggregate random distribution

    图  15  骨料强度对弹体偏转的影响

    Figure  15.  Projectile deflection affected by aggregate strength

    图  16  砂浆强度对弹体偏转的影响

    Figure  16.  Projectile deflection affected by motar strength

    图  17  不同入射速度下,不同γ值时,弹体偏转角度的时程曲线

    Figure  17.  Deflection angle-time curves of the projectiles at different γ values and different impact velocities

    表  1  不同级配骨料粒径分布

    Table  1.   Distribution of aggregates with different gradations

    骨料粒径/mm 骨料类型 骨料级配 Vs:Vm:Vb:Vh
    5~20 Small Stone 1st gradation 1:0:0:0
    20~40 Medium stone 2nd gradation 5.5:4.5:0:0
    40~80 Big stone 3rd gradation 3:3:4:0
    80~150 Huge stone 4th gradation 2:2:3:3
    下载: 导出CSV

    表  2  弹体及混凝土材料参数

    Table  2.   Material parameters of projectile and concrete

    材料类型 ρ/(kg·m-3) E/GPa μc σt/MPa σc/MPa
    Aggregate 2 660 0.16 10 160
    Motar 2 280 0.22 4 15
    ITZ 2 000 0.16 2 10
    Concrete 2 440 0.2 5 48
    Projectile 8 020 210 0.3
    下载: 导出CSV

    表  3  剩余速度对比

    Table  3.   Contrast of residual velocities

    v0/(m·s-1) vr/(m·s-1)
    实验 模拟
    360 67 73
    381 136 157
    434 214 237
    606 449 470
    749 615 633
    1 058 947 963
    下载: 导出CSV
  • [1] Backmann M E, Goldsmith W.The mechanics of penetration of projectiles into targets[J].International Journal of Engineering Science, 1978, 16(1):1-99. doi: 10.1016/0020-7225(78)90002-2
    [2] 钱伟长.穿甲力学[M].北京:国防工业出版社, 1984:23-35.
    [3] Goldsmith W.Review:Non-ideal projectile impact on targets[J].International Journal of Impact Engineering, 1999, 22(2/3):95-395. http://d.old.wanfangdata.com.cn/Periodical/zgyxllx201711029
    [4] Corbett G G, Reid S R, Johnson W.Impact loading of plates and shells by free-flying projectiles:A review[J].International Journal of Impact Engineering, 1996, 18(2):141-230. doi: 10.1016/0734-743X(95)00023-4
    [5] Li Q M, Reid S R, Wen H M, et al.Local impact effects of hard missiles on concrete targets[J].International Journal of Impact Engineering, 2006, 32(1/2/3/4):224-284. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f836952d7e21bc3f5f0eb712e6968392
    [6] Warren T L, Fossum A F, Frew D J.Penetration into low-strength (23 MPa) concrete:Target characterization and simulations[J].International Journal of Impact Engineering, 2004, 30(5):477-503. doi: 10.1016/S0734-743X(03)00092-7
    [7] Forrestal M J, Frew D J, Hickerson J P, et al.Penetration of concrete targets with deceleration-time measurements[J].International Journal of Impact Engineering, 2003, 28(5):479-497. doi: 10.1016/S0734-743X(02)00108-2
    [8] 陈小伟.穿甲/侵彻问题的若干工程研究进展[J].力学进展, 2009, 39(3):316-351. doi: 10.3321/j.issn:1000-0992.2009.03.006

    Chen Xiaowei.Advances in the penetration/perforation of rigid projectiles[J].Advances in Mechanics, 2009, 39(3):316-351. doi: 10.3321/j.issn:1000-0992.2009.03.006
    [9] Neville A M.Properties of concrete[M].5ed.Prentice Hall, 2012:112-119.
    [10] 何翔, 徐翔云, 孙桂娟, 等.弹体高速侵彻混凝土的效应实验[J].爆炸与冲击, 2010, 30(1):1-6. http://www.bzycj.cn/CN/abstract/abstract8811.shtml

    He Xiang, Xu Xiangyun, Sun Guijuan, et al.Experimental investigation on projectiles' high-velocity penetration into concrete target[J].Explosion and Shock Waves, 2010, 30(1):1-6. http://www.bzycj.cn/CN/abstract/abstract8811.shtml
    [11] Wang Z M, Kwan A K H, Chan H C.Mesoscopic study of concrete I:Generation of random aggregate structure and finite element mesh[J].Computers and Structures, 1999, 70(5):533-544. doi: 10.1016/S0045-7949(98)00177-1
    [12] Fuller W B, Thompson S E.The laws of proportioning concrete[M]//The American Society of Civil Engineers.Transactions of the American Society of Civil Engineers: Vol.LIX, 1926: 67-172.
    [13] 刘光廷, 高政国.三维凸型混凝土骨料随机投放算法[J].清华大学学报(自然科学版), 2003, 43(8):1120-1123. doi: 10.3321/j.issn:1000-0054.2003.08.033

    Liu Guangting, Gao Zhengguo.Random 3-D aggregate structure for concrete[J].Journal of Tsinghua University (Science and Technology), 2003, 43(8):1120-1123. doi: 10.3321/j.issn:1000-0054.2003.08.033
    [14] Vervuurt A.Interface fracture in concrete[D].Delft University of Technology, 1997: 45-49.
    [15] Hanchak S J, Forrestal M J.Perforation of concrete slabs with 48 MPa (7 ksi) and 140 MPa (20 ksi) unconfined compressive strengths[J].International Journal of Impact Engineering, 1992, 12(1):1-7. doi: 10.1016/0734-743X(92)90282-X
    [16] 梁斌.弹丸对有界混凝土靶侵彻研究[D].绵阳: 中国工程物理研究院, 2004: 62-89. http://cdmd.cnki.com.cn/Article/CDMD-82818-2004124640.htm
    [17] 陈小伟, 李继承.刚性弹侵彻深度和阻力的比较分析[J].爆炸与冲击, 2009, 29(6):584-589. doi: 10.3321/j.issn:1001-1455.2009.06.005

    Chen Xiaowei, Li Jicheng.Analysis on the penetration depth and resistive force in the deep penetration of a rigid projectile[J].Explosion and Shock Waves, 2009, 29(6):584-589. doi: 10.3321/j.issn:1001-1455.2009.06.005
    [18] Chen Xiaowei.Dynamics of metallic and reinforced concrete targets subjected to projectile impact[D].Singapore: Nanyang Technological University, 2003: 50-67.
    [19] 李志康, 黄风雷.混凝土材料的动态空腔膨胀理论[J].爆炸与冲击, 2009, 29(1):95-100. http://www.bzycj.cn/CN/abstract/abstract8874.shtml

    Li Zhikang, Huang Fenglei.dynamic spherical cavity-expansion theory for concrete materials[J].Explosion and Shock Waves, 2009, 29(6):95-100. http://www.bzycj.cn/CN/abstract/abstract8874.shtml
  • 期刊类型引用(1)

    1. 朱擎,李述涛,陈叶青,马上,魏万里,张生,陈嘉琳. 高强钢-钢筋混凝土复合防护结构厚度极限计算方法. 力学学报. 2024(07): 2077-2090 . 百度学术

    其他类型引用(1)

  • 加载中
图(17) / 表(3)
计量
  • 文章访问数:  4817
  • HTML全文浏览量:  1562
  • PDF下载量:  603
  • 被引次数: 2
出版历程
  • 收稿日期:  2015-10-09
  • 修回日期:  2016-03-27
  • 刊出日期:  2017-05-25

目录

/

返回文章
返回