Response of reinforced concrete slabs to low-velocity projectile impact investigated using upper bound method
-
摘要: 基于不可压缩刚塑性材料模型和滑移线场理论,获得了单一容许速度场条件下刚性弹低速侵彻半无限介质的阻力函数。在此基础上,基于多速度容许场得到了刚性弹侵彻有限厚度靶的三阶段阻力曲线,并提出了震塌与贯穿的临界条件,通过与实验结果、UMIST公式及古比雪夫的对比,验证了本文方法在钢筋混凝土板低速撞击问题中的适用性,分析了弹头形状、冲击因子和钢筋阻力系数等参数对临界震塌(贯穿)厚度的影响。Abstract: Based on the incompressible-rigid-plastic material assumption and the slip line field theory, the resistance function of a rigid projectile penetrating a semi-infinite target at a low velocity was obtained with a single admissible velocity field. A three-stage resistance curve of a rigid projectile impacting on a thin target was analyzed under multiple velocity fields, where the critical conditions for scabbing or perforation were calculated. The methods and formulae for local effects on reinforced concrete slab under low-velocity impact were further verified using comparative analysis of the results from the experiments, the UMIST formulae, the Kuibyshev formulae, and the present paper's calculations. The relationships between the normalized critical scabbing/perforation thickness, and the nose-shape factor, the impact factor and the reinforcement factor were examined to present potential guide to experimental studies.
-
Key words:
- rigid-plasticity limit analysis /
- low velocity penetration /
- concrete /
- scabbing /
- perforation
-
液舱结构在弹丸侵彻下的动态响应问题近年来受到工程界广泛关注。战斗部接触爆炸会产生大量高速破片,对舰船舷侧防护液舱结构的侵彻会造成严重毁伤。针对高速弹体侵彻液舱结构的破损问题已开展了大量研究工作。J.H.McMillen等[1]早在20世纪40年代就针对弹体入水冲击波传播等问题通过实验进行了分析。D.Townsend等[2]则结合实验对高速破片侵彻薄壁水箱结构引起的水锤效应做了细致的研究,并提出了衰减冲击波压力的初步构想。随后,R.J.Disimile等[3]通过实验分析了钢质和铝质球形弹丸穿透背水薄铝板后在水中引起的水锤效应。D.Varas等[4]则从压力变化和管壁变形等方面对球形弹丸侵彻注水方管引起的水锤效应进行了研究,并分析了注水液面的影响。E.Deletombe等[5]比较了密闭和敞开水箱两种情形下水锤效应引起的冲击波压力差异。沈晓乐等[6]、徐双喜等[7]、李典等[8]针对高速破片侵彻液舱问题也开展了大量研究。通过分析可知,目前针对破片侵彻液舱问题的研究大多集中于穿透背水钢板后引起的水锤效应,对于破片穿透后的剩余特性研究较少,且研究对象大多为垂直背水薄板。本文中通过弹道实验,分析破片高速侵彻中厚背水钢板后的剩余特性,从弹体和靶板的破坏模式以及破片穿透背水钢板后的瞬时余速等方面,探讨破片高速侵彻倾斜和垂直背水钢板后的剩余特性的差异,力图揭示中厚背水钢板的抗高速侵彻机理,为后续的理论分析和数值模拟研究提供数据支撑。
1. 实验设计
实验采用15 mm口径滑膛弹道枪,通过火药发射。破片初速采用激光测速系统(1#和2#光幕靶)得到。实验弹丸采用质量3.3 g、边长7.5 mm的立方体破片。破片材料为经淬火处理的45钢。实验靶板为5 mm左右的Q235船用低碳钢。为模拟破片垂直和斜侵彻背水钢板的情形,制作了垂直90°和倾斜60°两种水箱结构,分别如图 1~2所示。
2. 弹道实验结果
采用高速摄影机对破片穿透背水钢板后在水中的运动过程进行了观测,得到背水钢板情形下破片余速即弹丸穿透钢板入水瞬间的速度。表 1中给出了弹道实验结果及主要参数。
表 1 实验结果及主要参数Table 1. Experimental results and main parameters工况 靶板厚度/mm 靶后介质 入射角/(°) 破片初速/(m·s-1) 破片余速/(m·s-1) 破坏情况 1 4.98 水(水箱Ⅰ) 90 1 105.0 229.3 穿透 2 4.96 水(水箱Ⅰ) 90 1 231.2 347.1 穿透 3 4.98 水(水箱Ⅱ) 60 1 058.1 167.8 穿透 4 4.94 水(水箱Ⅱ) 60 1 290.3 310.5 穿透 5 4.92 空气 90 1 029.0 384.0 穿透 6 4.96 空气 60 1 029.9 224.1 穿透 3. 实验结果及分析
3.1 弹体破坏模式分析
图 3给出了实验后收集的破片变形破坏形貌。由图可知,立方体破片不管是侵彻背水钢板还是背空钢板,头部均发生了严重的镦粗变形。对于垂直侵彻情形,破片头部近似呈圆形;而对于60°斜侵彻情形,破片头部虽然也近似呈圆形,但圆形中心向入射速度反方向偏移。进一步观察可知,立方体破片基本上是以面接触的形式侵彻靶板,且能够较好地保持侵彻姿态。
3.2 靶板破坏模式分析
由于钢板背后有液体存在,使得背液钢板的抗侵彻过程及耗能机制与背空钢板[9]差异较大。图 4中给出了中厚背水钢板的典型破坏形貌。由图 4可知,高速侵彻下,中厚背水钢板的穿甲破坏模式主要是剪切冲塞破坏。对于垂直和60°斜侵彻两种侵彻情形,主要区别在于穿孔形状。进一步观察图 5可知,背水钢板在厚度方向存在一定的绝热剪切现象。这说明高速侵彻背水钢板过程中,绝热剪切效应不容忽视。而未考虑绝热剪切效应的弹体侵彻背水钢板的能量耗散理论计算[7]是不符合实际情况的。
3.3 破片在水中运动轨迹分析
破片在水中运动的过程主要包括拖曳和空化两个阶段[3-4]。破片在水中运动过程中,垂直方向则受到水动升力和重力及低压效应的共同作用[10]。
结合图 6~9可知,破片初速较大时,即穿透背水钢板后的瞬时余速较大时,破片在水中的运动轨迹向上偏,如工况2(图 7)和工况4(图 9);相反,则向下偏,如工况3(图 8);若破片速度能使水动升力和重力及低压效应保持平衡,则破片的水中运动轨迹近似成水平,如工况1(图 6)。进一步比较可看出,破片在水中运动轨迹的偏转主要还是受入水初速的影响,钢板倾斜角度对轨迹偏转的影响较小。
结合图 6~7可看出,破片入水初期即产生了明显射流,射流方向与水平方向夹角均为约45°。这与弹体低速入水深闭合所产生的铅垂方向的射流[10]有较大差异。进一步分析图 8~9可看出,虽然破片在穿透倾斜背水钢板后也产生了上下方向的射流,但由于倾斜钢板的影响,向上射流的方向平行于钢板平面,与水平方向成60°夹角;向下射流的方向仍与水平方向成45°夹角。进一步结合图 6~9可知,破片入水初期形成的射流强度与破片入水初速有关,破片入水初速越大,射流越强。
进一步分析图 6~9中空泡的相关特性可知,破片穿透垂直背水钢板后的空泡形状比较对称,而穿透60°倾斜背水钢板后形成的空泡则往下偏。从空泡的大小来看,主要还是取决于破片穿透背水钢板后的余速即破片入水初速。
3.4 破片穿透背水钢板后的余速分析
通过分析高速摄影图片得到破片穿透背水钢板后的瞬时余速。通过拟合求导得到各工况破片穿透背水钢板后的瞬时余速,如表 1所示。由表 1通过计算可知,随着破片初速提高,背水钢板单位面密度吸能值相应提高。破片初速越高,提高值越大。相对于背空钢板而言,相同破片初速下,背水钢板的单位面密度吸能要稍大。
3.5 初始冲击波压力特性分析
大量研究结果表明[2, 5, 8],弹体高速侵彻水介质时形成的冲击波压力峰值较大。因此,本节主要分析初始冲击波压力特性。假设稀疏波以对称于自由液面虚拟镜像波源为起点向外传播,则初始冲击波和稀疏波的传播过程如图 10所示。以初始冲击波产生时刻为起点,得到初始冲击波传播至压力测点的时间t1为
t1=√L2h+B2+(D+h−H)2/cw (1) 式中:Lh为压力测点离撞击点的水平距离,B为压力测点离轨迹线的水平横向距离,D为撞击点离自由液面的高度,h为压力测点离水底高度,H为自由液面高度,cw为水中声速。
稀疏波到达压力测点所需的时间t2为
t2=√L2h+B2+(D+H−h)2/cw (2) 由于H>h,因此t2始终要大于t1,但随着水平距离Lh的增大,t2与t1的差值变小。
考虑到初始冲击波有一定的持续作用时间Δt,若稀疏波在初始冲击波持续作用时间内到达,则需考虑稀疏波对初始冲击波的影响;反之,稀疏波对初始冲击波无影响。由此得到以下判据:
{t2<t1+Δt稀疏波对初始冲击波有影响t2≥t1+Δt稀疏波对初始冲击波无影响 (3) 对于线性可压缩介质而言,波传播具有频率恒定的特性[11]。因此,初始冲击波在水中传播时,距离上的几何衰减仅对冲击波载荷幅值即峰值压力大小产生影响,对载荷的持续作用时间没有影响[4]。因而某一压力测点上,初始冲击波持续作用时间Δt可按下式计算[12]:
Δt=3Rv0w/(2cw) (4) 式中:R为被压缩水介质的区域半径,v0w为弹丸(破片)入水初速。
若压力测点位于水平弹道轨迹上,则稀疏波和初始冲击波传播至该点的时间差ΔT为
ΔT=t2−t1=(√L2h+4D2−Lh)/cw (5) 假设该压力测点的初始冲击波峰值压力与距离的平方成反比,且近似呈三角形脉冲形式,则考虑稀疏波的影响下,该压力测点的比冲量I近似等于:
I≈12P(Lh)ΔT=P(L0)(L0Lh)2√L2h+4D2−Lh2cw (6) 式中:P(L0)为距离撞击点L0(L0 < Lh)处的峰值压力值。
将上式对D求导得到:
∂I∂D=P(L0)(L0Lh)4Dcw√L2h+4D2 (7) 由上式可以看出,撞击点水深D对压力测点比冲量I的影响与该点离撞击点距离Lh的3次方成反比,即稀疏波对初始冲击波的影响程度随距离的3次方倒数下降,说明距离越远,稀疏波的影响越小。
图 11中给出了典型工况(工况2)中压力测点1~2的初始冲击波压力时程曲线。采用式(1)~(4)计算得到:工况2中初始冲击波到达测点2和测点1的时间差约为49 μs,而实验得到两测点峰值时间差约为51 μs,两者吻合较好。结合式(3)的判据可得,对于测点1,稀疏波无影响;而对于测点2,稀疏波有影响。而实际情况是,测点2受到了稀疏波影响,测点1未见明显影响。因此,考虑破片后续运动和空化效应影响,本文中对测点1~2中稀疏波是否影响的预测判别还是比较符合实验结果的。
表 2中给出了各工况下压力测点1~2各自的初始冲击波峰值压力大小。由表 2可知,各工况中测点1和测点2处的初始冲击波峰值压力大小顺序是相同的,即工况2>工况4>工况1>工况3。由于初始冲击波峰值压力大小与破片的入水初速近似呈线性增大[2, 4, 8, 12],因而由峰值压力大小顺序可反推得到破片入水初速即破片穿透背水钢板后的瞬时余速的大小顺序应该相同,这与上一节的分析结果是一致的。
表 2 破片侵彻背水钢板各工况测点1~2的峰值压力测量结果Table 2. Measured peak pressures of test points 1 and 2 for cases of water-bakced steel plates penetrated by fragments工况 破片初速v0/(m·s-1) 入射角/(°) 破片余速vr1/(m·s-1) 峰值压力P1max/MPa 峰值压力P2max/MPa 1 1 105.0 90 286.0 6.184 4.580 2 1 231.2 90 462.7 7.707 6.132 3 1 058.1 60 202.9 4.322 3.894 4 1 290.3 60 385.9 6.432 5.563 4. 结论
通过弹道实验,从弹体和靶板的破坏模式、破片穿透背水钢板后的瞬时余速等方面,分析了破片高速侵彻中厚背水钢板的剩余特性,得出主要结论如下:
(1) 破片高速侵彻中厚背水钢板过程中,头部会发生严重镦粗变形,而镦粗变形主要发生在侵彻背水钢板阶段,但背部水介质的“动支撑”作用不可忽视;
(2) 破片高速侵彻下,中厚背水钢板的主要破坏模式为剪切冲塞,迎弹面出现靶材的反向飞溅,而背弹面穿孔则存在绝热剪切现象;
(3) 破片穿透背水钢板初期,会产生空泡和射流,空泡大小和射流强度主要与破片入水初速有关,而空泡形状和射流方向则受背水钢板的倾斜角度影响;
(4) 破片穿透背水钢板后,在水介质中的运动轨迹会发生偏转,偏转方向主要与破片入水初速有关,破片初速较大时,运动轨迹向上偏,反之则向下偏;
(5) 由于水介质的“动支撑”作用和动能耗散效应,破片穿透背水钢板后的瞬时余速较穿透背空钢板后的余速要低,动能损失要大;
(6) 由于有自由液面的存在,水介质中初始冲击波压力特性需考虑稀疏波的影响,而稀疏波影响的可能性随离撞击点的距离的增大而提高,但影响程度会大大降低;
(7) 相同初速条件下,破片侵彻垂直背水钢板引起的初始冲击波的峰值压力大小较侵彻倾斜背水钢板要大,但峰值压力的提高幅度随破片初速的增大或离撞击点距离的增大均相应减弱。
-
实验 L/2a fc /MPa ft/MPa τs/MPa IBL vBL/(m·s-1) 实验 本文计算 T-1 2 25.0 2.6 4.0 10.5 27.0~35.7 47 T-2 2 25.2 3.1 4.4 10.5 41.7~56.8 49 T-3 2 161.9 7.3 17.2 10.5 34.7~58.5 97 T-4 2 175.3 13.8 24.6 10.5 76.0~104.0 116 表 2 本文计算结果与普通混凝土实验[11]的比较
Table 2. Comparison between experimental results [11] of normal strength concrete and present method
实验 L/2a fc /MPa ft/MPa τs/MPa IBL φ/mm @/mm fs/MPa μ vBL/(m·s-1) 配筋 实验 本文计算 D-1-1 2.0 35.0 3.0 5.12 15.9 2.5 34 382 0.43 165 103 D-1-2 2.0 35.0 3.0 5.12 15.4 3.0 25 183 0.40 222 101 D-1-3 2.4 35.0 3.0 5.12 20.8 3.0 25 183 0.40 232 117 D-1-4 2.4 35.0 3.0 5.12 127.1 5.0 27 473 2.68 236 291 D-1-5 2.0 34.0 3.4 5.38 53.6 3.25 21 600 1.76 210 193 D-1-6 2.0 34.0 3.4 5.38 34.8 2.5 20 650 1.19 162 156 表 3 本文计算结果与高性能混凝土实验[12]的比较
Table 3. Comparison between experimental results[12] of high performance concrete and present method
实验 L/2a fc /MPa ft/MPa τs/MPa IBL φ/mm @/mm fs/MPa μ vBL/(m·s-1) 配筋 实验 本文计算 D-2-1 4 40.0 4.0 6.3 66.5 8 100 400 0.64 204~245 187 D-2-2 4 108.0 10.8 17.1 38.3 8 100 400 0.24 273~276 233 D-2-3 4 102.0 10.2 16.1 39.1 8 100 400 0.25 281~289 229 D-2-4 4 104.0 10.4 16.4 38.9 8 100 400 0.24 287~291 231 D-2-5 4 113.0 11.3 17.9 37.7 8 100 400 0.22 262~289 237 D-2-6 4 106.0 10.6 16.8 38.6 8 100 400 0.24 291~307 232 D-2-7 4 101.0 10.1 16.0 39.3 8 100 400 0.25 286~292 229 D-2-8 4 93.0 9.3 14.7 40.7 8 100 400 0.27 292~313 223 D-2-9 4 94.0 9.4 14.9 40.5 8 100 400 0.27 313~314 224 D-2-10 4 102.0 10.2 16.1 39.2 8 100 400 0.25 287~289 229 D-2-11 4 103.0 10.3 16.3 39.0 8 100 400 0.25 287~292 230 -
[1] Hill R. The mathematical theory of plasticity[M]. Oxford: Clarendon Press, 1950. [2] Ravid M, Bodner S R. Dynamic perforation of visco-plastic plates by rigid projectiles[J]. International Journal of Engineering Science, 1983, 21(6):577-591. doi: 10.1016/0020-7225(83)90105-2 [3] Amini A, Anderson J. Modeling of projectile penetration into geologic targets based on energy tracking and momentum impulse principles[C]//Proceedings of the Sixth International Symposium on Interaction of Nonnuclear Munitions with Structures. Germany, 1993. [4] Tirosh J, Tylis A, Davidi G. Foreign object damage: Penetration of rigid projectiles into solids[J]. Mechanics of Materials, 2009, 41(5):535-544. doi: 10.1016/j.mechmat.2009.01.024 [5] 陈士林, 王明洋, 潘越峰.锥形端部弹体在岩石(混凝土)介质层中侵彻实用计算方法[J].爆炸与冲击, 2004, 24(1):7-15. http://www.bzycj.cn/article/id/9912Chen Shilin, Wang Mingyang, Pan Yuefeng. The method of calculation for penetration of conical nosed projectile in rock (concrete) layers[J]. Explosion and Shock Waves, 2004, 24(1):7-15. http://www.bzycj.cn/article/id/9912 [6] 王明洋, 施翠英, 陈士林.事故型撞击混凝土板的临界震塌与贯穿厚度计算方法[J].工程力学, 2009, 26(11):238-246. http://www.cnki.com.cn/Article/CJFDTotal-GCLX200911042.htmWang Mingyang, Shi Cuiying, Chen Shilin. Method of calculating critical spalling and penetration thickness of concrete slab of block under accident impact[J]. Engineering Mechanics, 2009, 26(11):238-246. http://www.cnki.com.cn/Article/CJFDTotal-GCLX200911042.htm [7] 冯淑芳, 王明洋, 任广昊, 等.深地下坑道衬砌结构抗局部冲击计算方法研究[J].岩石力学与工程学报, 2011, 30(6):1150-1156. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201106008Feng Shufang, Wang Mingyang, Ren Guanghao, et al. Research on calculation method for local impact of deep tunnel lining structure[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(6):1150-1156. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201106008 [8] 咸玉席, 文鹤鸣.平头弹侵彻半无限混凝土靶的工程模型[J].防护工程, 2012, 34(2):35-38. http://cdmd.cnki.com.cn/Article/CDMD-10358-1014432193.htmXian Yuxi, Wen Heming. An engineering model for the penetration of flat-nosed projectiles into semi-infinite concrete targets[J]. Protective Engineering, 2012, 34(2):35-38. http://cdmd.cnki.com.cn/Article/CDMD-10358-1014432193.htm [9] Li Q M, Ye Z Q, Ma G W, et al. Influence of overall structural response on perforation of concrete targets[J]. International Journal of Impact Engineering, 2007, 34(5):926-941. doi: 10.1016/j.ijimpeng.2006.03.005 [10] Tai Y S. Flat ended projectile penetrating ultra-high strength concrete plate target[J]. Theoretical and Applied Fracture Mechanics, 2009, 51(2):117-128. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0214989790/ [11] Dancygier A N. Effect of reinforcement ratio on the resistance of reinforced concrete to hard projectile impact[J]. Nuclear Engineering and Design, 1997, 172(1):233-245. http://cn.bing.com/academic/profile?id=723b6881b18fa7ba33bf200ca91ae81f&encoded=0&v=paper_preview&mkt=zh-cn [12] Dancygier A N, Yankelevsky D Z, Jaegermann C. Response of high performance concrete plates to impact of non-deforming projectiles[J]. International Journal of Impact Engineering, 2007, 34(11):1768-1779. doi: 10.1016/j.ijimpeng.2006.09.094 [13] 卡恰诺夫L M.塑性理论基础(第二版)[M].周承倜, 译.北京: 人民教育出版社, 1982: 332-385. [14] 冯淑芳.深埋地下结构抗爆动力计算方法研究[D].南京: 解放军理工大学, 2011. [15] 徐秉业, 刘信声.应用弹塑性力学[M].北京:清华大学出版社, 1995:510-512. [16] Wen H M, Xian Y U. A unified approach for concrete impact[J]. International Journal of Impact Engineering, 2015, 77:84-96. doi: 10.1016/j.ijimpeng.2014.11.015 -