Chinese Journal of Ship Research
钢/玻璃钢组合结构对高速弹丸的抗侵彻特性
张元豪，陈长海，朱锡430033海军工程大学 舰船工程系，湖北 武汉
摘 要：为探究钢与玻璃钢的组合结构形式对舰船舱壁复合装甲结构抗穿甲性能的影响，采用均质钢板前置和2后置玻璃钢来分别模拟舰船舱壁外设及内设复合装甲结构，结合高速弹道冲击实验，分析、比较 种结构形式组ANSYS/LS-DYNA合靶板的穿甲破坏模式和抗弹吸能能力。在此基础上，利用有限元分析软件 开展高速立方体弹丸侵彻组合靶板的数值模拟计算，分析组合靶板的侵彻过程，并与实验结果进行比较。结果表明，数值计算结果与实验结果较为吻合；2种组合靶板中复合装甲板的破坏模式均主要为钢板的剪切冲塞破坏和玻璃钢的纤维剪切断裂，后置组合靶板中玻璃钢背层伴随有纤维的拉伸破坏；前置组合靶板的抗弹吸能能力要稍大于后置组合靶板。关键词：舰船舱壁；复合装甲；抗弹性能；钢/玻璃钢组合结构
0引言
随着反舰导弹的不断发展，掠海飞行的半穿甲反舰导弹由于具有突防能力强以及侵入舷侧内爆等特点，已成为现代水面舰船面临的主要威胁。而战斗部内爆产生的高速破片会对舰船重要舱室造成二次毁伤，所以舰船舱壁的抗动能穿甲防护问题就显得尤为突出。目前，对于高速破片的动能穿甲作用，工程中主要是通过增加舱壁板厚或设置复合装甲的方式来尽可能减小对舱室的破坏。针对半穿甲战斗部内爆产生的高速破片的动能穿甲作用，现代舰船舱室广泛采用设置复合装Fiber Rein⁃甲的结构形式。纤维增强复合材料（ forced Plastics，FRP）以其高比强度和比刚度等优点，在舰船复合装甲防护结构中得到了广泛的应用［1-2］。但层合板结构及材料特性的复杂性导致其弹道冲击过程十分复杂，影响因素也较多，学者们致力于研究侵彻和穿透过程中复合靶板材料的动态力学变形机理。对于单一复合材料板或均质钢板的抗穿甲问题，国内外开展了大量的研究［3-5］。FRP的抗侵彻性能研究，Greaves［6-7］利用平头针对 S2 Glass /Phenolic弹撞击 厚层合板靶板，研究了其弹道穿透过程中的破坏机理；Zhu等［8-9］对锥头弹撞Kevlar -29/Polyester击 层合板的吸能机制进行了Wen 等［10-12 ］提出了针对不同形状弹体侵彻研究； FRP层合板预测侵彻深度和弹道极限的公式。覃FRP ［13］悦等 研究了卵形弹丸撞击下 层合板的穿透FRP性能；张颖军等［14］介绍了 层合板损伤性能的表征方法；谢恒等［15］采用有限元方法对不同破片Glass Fiber Reinforced Plas⁃模拟弹侵彻玻璃钢（ tics，GFRP）层合板的动态响应进行了分析。而对于纤维复合材料板与均质钢板组合形成的结构靶板的研究却并不多见。舰船舱壁结构相对于高速破片而言属于中厚板，所以高速破片对舰船舱壁结构的侵彻可以认为是中厚板结构抗高速穿甲的问题。为探讨舰船舱壁内设和外设复合装甲结构抗动能穿甲破坏机理的差异，比较内设与外设复合装甲结构的抗穿甲性能，本文将以均质钢板前置和后置玻璃钢分别模拟舰船舱壁外设与内设复合装甲结构，结合高速弹道冲击实验，分析、比较前置和后置组合靶板穿甲的破坏模式和抗弹性能。在此基础上，利ANSYS/LS-DYNA用有限元分析软件 开展高速立方体弹丸侵彻组合靶板的数值模拟计算，分析组合靶板的侵彻过程和破坏模式，并与相应的实验 结果进行比较。
1 实验及实验主要结果
14.5 mm实验采用 口径的滑膛弹道枪系统发射弹体，通过火药推进。为保证发射所必须的密封性和达到规定的速度，破片用特制的三瓣式铝合金弹托包覆，破片出膛后弹托通过弹托回收器分离，同时保持破片的运动轨迹。该系统包括激光测速装置。7.5 mm，设计质量 3.30 g采用边长 的立方体45#钢加工而成。弹，弹体由经淬火处理的Q235钢质靶板采用 钢，正方形靶板尺寸为400 mm × 400 mm × 5 mm 300 mm × ，暴露面积为300 mm 39 kg/m2 ，面密度为 ；玻璃钢靶板采用SW220 200 mm ×玻璃纤维布制作而成，尺寸为200 mm × 5 mm 12 kg/m2。二者之间采，面密度为用环氧树脂粘接形成组合靶板。实验采用靶架固定的形式，事先将靶架固定于滑道上，在实验过程中，再将靶板固定在靶架上。45# Q235 SW220钢、 钢及 玻璃纤维的材料参1 2数如表 和表 所示。
2 有限元仿真结果
LS-DYNA，建立高速弹体侵利用有限元软件彻组合结构靶板的数值仿真模型。弹体和靶板均
Lagrange采用八节点拉格朗日（ ）实体单元模拟， Lagrange采用 网格建模，弹体各边均八等分；在靶50 mm×50 mm，板中心选取正方形区域，其尺寸为 50对该撞击区域进行网格加密处理，各边平分 等分，远离该区域网格向靶板四角稀疏过渡；靶板厚10 1度方向 等分。网格模型如图 所示。 Plastic_弹体采用双线性弹塑性本构模型Kinematic Cowper-Symonds ，其应变率效应由 模型描述： EEh 1 N σd = σ0 + +è Dø ε 1 ae ε̇ ö E - E p h式中：σd 为动态屈服强度；σ0为静态屈服强度；E为弹性模量； E 为硬化模量； ε 为有效塑性应h p变；ε̇ 为等效塑性应变率；D，N为常数，对于低碳D=40.4 s-1，N=5。材料失效模型采用最钢，通常取大等效塑性应变失效准则。3弹体的材料参数［16］如表 所示。表3 弹体的材料参数Table 3 Material parameters of projectile /MPa /MPa D/s-1 σ0 E N失效应变 εf h 335 350 5 40.4 0.7 Johnson-Cook钢板材料采用 本构模型，该模型考虑了应变率强化、由绝热升温引起的软化效应，能反映材料在高应变率以及高温情况下材料性质的变化，其状态方程为m εp T - T0 pn σ = (A + Bε )1 + c ln 1aeçè y ε0 T -T m 0式中：σ 为钢板材料的动态屈服强度；A为静态屈y 235 MPa；B服强度，取为 为应变硬化模量，取为 300 MPa；n 0.26；c为应变硬化指数，取为 为应变率0.014；m 1.03；εp系数，取为 为热软化指数，取为为等效塑性应变； ε0 为参考塑性应变率，一般取1 s-1 1 793 K ；T 为材料的熔点，取为 ；T 为参考m 0 300 K。温度（取室温），取为由下述方程描述受冲击区域材料的失效： σh ε p ε = D1 + D exp D 1 + D4 · f 2 3 σeff ε̇ 0 T - T0 1 + D5 Tm -T 0 ~ Mises式中： D1 D 为材料常数； σeff 为 等效应5力；σ 为材料在三向应力状态下的静水压力；当h p破坏参数 D = åDε = 1时，材料发生失效。钢板εf 4的材料参数［16］如表 所示，其中G为剪切模量。表4 钢板的材料参数Table 4 Material parameters of steel plate G/GPa A/MPa B/MPa Tm/K T0/K n c m D1 D2 D3 D4 D5 80.8 235 300 0.26 0.014 1.03 1 793 300 0.4 0 0 0 0 Johnson-Cook CompositeDamage玻璃钢材料采用本构模型，该模型材料的应力应变关系为ε = 1 ( - 2) σ1 υ1σ 1 E1 ε = 1 ( σ - υ σ1) 2 2 2 E 2 1 τ + 2ε12 = ατ12 12 G12式中：ε1 ，ε ，ε12 分别为材料纵向、横向的拉伸应2变和面内剪应变；σ1 ，σ ，τ12 分别为材料纵向、横2向的拉伸应力和面内剪应力；υ1 ，υ 分别为纵向2和横向泊松比；E ，E ，G 分别为纵向、横向的1 2 12拉伸模量和面内剪切模量；α为非线性剪切应力参数。Chang-Chang 3种该模型采用的 失效准则的失效准则如下。1）基体开裂失效准则： σ2 - F = ( )2 + τ matrix S 2 >1当 F 时，产生基体开裂失效，材料常数matrix 0。E ，E ，υ ，υ都设为1 2 1 2 2）压缩失效准则： σ C σ - F =( )2 + [( )2 - 1] + τ 2 2 2 comp S12 2S12 C2 >1当 F 时，发生材料压缩失效，材料常数comp 0。，υ1 ，υ 都设为E 2 2
3）最终的失效模式是纤维断裂。σ1 - F = ( )2 + τ fiber S1 >1当 F 时，发生纤维断裂失效，材料常数fiber , , , , 0。玻璃钢的材料参数E1 E G12 υ1 υ 都设为2 2 5如表 所示。通过对实验结果的观察，表明弹体的侵彻为6垂直于靶板的正侵彻。表 所示为弹道实验与数
3 实验及计算结果分析 3.1 靶板侵彻过程及破坏模式分析 3.1.1 靶板类型Ⅰ的实验结果
1工况 中，前置玻璃钢和后置钢板的破坏形 3 3图 所示为对工况 的有限元模拟过程。由图可看出，弹体在侵彻玻璃钢的过程中使纤维层3（a））。弹体发生断裂，对其产生了剪切破坏（图继续侵彻到钢板，高速侵彻产生冲击波，冲击波反射造成弹体侵蚀，弹体因此进一步墩粗且侵彻速度下降明显，玻璃钢在侵彻区出现局部的纤维脱3（b））。钢板的破坏模式为剪切破坏，胶现象（图并在整个弹体侵彻过程中不断吸收弹体的侵彻能 表5 玻璃钢的材料参数Table 5 Material parameters of GFRP E11/GPa E22/GPa E33/GPa G12/GPa G13/GPa G23/GPa υ12 υ13 υ 23 18.22 18.22 6 0.12 0.3 0.3 6.75 6.75 3 Ⅰ为值仿真结果以及相关参数。表中，靶板类型5mm 玻璃钢+5 mm Ⅱ为钢板组合靶板，靶板类型5mm钢板+5 mm 玻璃钢组合靶板。 2貌如图 所示。由图可看出，前置玻璃钢在弹体侵彻下，呈现出纤维脱胶现象，弹体穿甲破坏模式主要为纤维的剪切断裂破坏，纤维的断裂面较为齐整；而后置钢板的穿甲破坏模式则主要是剪切冲塞破坏。由于弹体在穿透前置玻璃钢后速度较低，破片未完全穿透后面的钢板，剪切塞块未脱离钢板，导致钢板在破片穿甲部位形成局部的隆起。 量，弹体速度因此持续下降。在后期的侵彻过程中，钢板发生失效破坏，弹体也不断出现镦粗，直3 d至穿透整个后置钢板（图 （ ））。3工况 中弹体完全穿透了靶板，靶板的破坏4 3（e）与图4形貌如图 所示。由图 可看出，相比于1，前置玻璃钢的穿甲部位除发生纤维剪切工况断裂外，穿孔周围还出现了大量的纤维脱胶以及纤维断裂后的原纤维化现象。这主要是由于在速
度较高的情形下，弹体在剪切穿甲过程中与纤维的摩擦作用更剧烈，摩擦产生的热量更不容易传递，因而导致纤维被大量熔断。通过观察后置钢板的破坏形貌及穿孔周围靶材的颜色，并结合3（f）可知，其穿甲破坏模式主要是剪切冲塞，剪图切冲塞穿孔的周围也出现了一定程度的塑性变形。这主要是由于前置玻璃钢的存在降低了弹体 冲击后置钢板时的速度，同时断裂的纤维附着在弹头表面增大了弹体对后置钢板的作用面积。
3.1.2 靶板类型Ⅱ的实验结果
5 2图 所示为工况 中靶板的破坏形貌。由图可知，前置钢板在冲击区边缘存在明显的剪切冲塞痕迹，冲击区外围存在大量的变形。由此可得，其主要破坏模式为剪切冲塞破坏。钢板穿孔形状近似为圆形，穿孔周围有材料“反向溢出”现象，这是由弹体在穿甲过程中对近似流体状靶板材料的挤压作用造成的。而后置玻璃钢则未被穿透，主要穿甲破坏形式为迎弹面的纤维剪切断裂和背弹面的分层破坏，从后置玻璃钢的破坏形貌可看出，其背层在穿孔附近出现了纤维脱胶现象，而背层未穿透的纤维层形成了分层和局部变形。6 4图 所示为工况 的有限元模拟侵彻过程。由图可知，弹体接触钢板对其进行了剪切破坏，同6（a））。钢板失效后，将弹体时弹体出现墩粗（图的侵彻载荷传递到玻璃钢，在侵彻载荷作用下，玻璃钢背板纤维出现层裂且伴有纤维剪切与拉伸破坏，在整个侵彻过程中弹体侵彻能量不断被吸收， 6（b））。导致其速度的持续下降（图 而钢板冲塞块在弹体与玻璃钢之间进一步侵蚀弹体，背板纤维层完全破坏，失去防护能力，弹体携带钢板和纤6（d））。维板冲塞块穿透靶板（图
4 7工况 靶板的破坏形貌如图 所示，由图可看2出后置玻璃钢的破坏较工况 更为严重，结合6（f）分析表明，其主要体现在背层未被穿透的图纤维层的脱胶和纤维的原纤化现象上，并且伴随7（a 6（e）可有明显的纤维拉伸破坏。由图 ）与图知，前置钢板的破坏模式为剪切冲塞破坏。
3.2 剩余速度结果比较
6表 示出了实验与有限元模拟剩余速度结果，有限元模拟中弹体初速与实验工况一致。由3，4），实表可知，当初速较高时（工况 验剩余速度与有限元模拟剩余速度的相对误差较小，这说明3 4数值仿真结果是可靠的。工况 与工况 中的弹4体初速相差不大，而工况 弹体的剩余速度是工
3 2 1 000 m/s况 的近 倍。当弹体初速在 左右时， 2（b）可知，后置钢板凸包存在裂缝且仅有一由图小部分未脱离钢板，处于临界击穿状态；由5（b）观察到，弹体已经穿透至玻璃钢最后一层图的纤维层，所以通过分析实验后板的变形破坏及有限元模拟，近似认为靶板类型Ⅰ在 =1 050 m/s v0 Ⅱ =1 000 m/s时，靶板类型 在 v0 时均达到了弹道极限状态。
3.3 靶板抗侵彻性能对比分析
上节中提到本文认为靶板类型Ⅰ、靶板类型Ⅱ 1 050 1 000 m/s。因此，从的弹道极限分别约为 和弹道极限的角度来看，靶板类型Ⅰ要稍好于靶板类型Ⅱ。进一步比较工况3 4可知，2和工况 种组合靶板在弹体初速相差不大的穿透情形下，前者42.68（J·m ）/kg，后者的2的单位面密度吸能值为41.64（J·m ）/kg 2单位面密度吸能值为 。由此可Ⅰ Ⅱ见，靶板类型 的抗弹吸能较靶板类型 要大。这主要是因为弹体侵彻玻璃钢+钢板组合靶板时，一方面钢质背板对玻璃钢层有动支撑作用，增大了弹体在侵彻玻璃钢层过程中的局部惯性质量；另一方面，弹体在侵彻前置玻璃钢时，弹体侵彻所产生的压缩应力波传播至后置钢板，使得后置钢+玻板参与吸能并产生压缩功。而弹体侵彻钢板璃钢组合靶板时，弹体在侵彻前置钢板过程中所产生的压缩应力波使玻璃钢背层出现了轻微的层5（b）），这间脱胶现象（图 反而会降低后置玻璃钢的抗侵彻能力，从而在一定程度上降低钢板+玻璃钢组合靶板的整体抗侵彻能力。
4结论
本文通过高速弹道冲击实验并利用有限元分ANSYS/LS-DYNA析软件 的数值模拟计算，分析、2比较了 种结构形式组合靶板的穿甲破坏模式和抗弹吸能能力，得出如下主要结论： 1）在高速立方体弹侵彻下，钢板+玻璃钢组合靶板结构中钢板的破坏模式为剪切冲塞破坏，玻璃钢的破坏模式主要为纤维的剪切断裂破坏，玻璃钢背层伴随有纤维的层间脱胶及拉伸破坏。2）在高速立方体弹侵彻下，玻璃钢+钢板组合靶板结构中钢板的破坏模式为剪切冲塞破坏，玻璃钢的破坏模式为纤维的剪切断裂破坏。3）在高速立方体弹侵彻下，由于钢背板的支撑作用以及其压缩应力做功的影响，玻璃钢+钢板+玻璃钢的组组合靶板的抗弹性能要稍强于钢板合靶板。
参考文献：
1］朱锡，梅志远，刘润泉，等. ［ 舰用轻型复合装甲结构及其抗弹实验研究［J］. 爆炸与冲击，2003，23（1）： 61-66. ZHU X， MEIZY LIURQ ，et al. Warship's light com⁃ ， posite armor structure resistibility for ballistic impact ［J］. Explosion and Shock Waves， 2003， 23（1）： 61-66（in Chinese）. 2］朱锡，侯海量，谷美邦，等. ［ 抵御小口径火炮弹道侵J］. 2006，彻装甲防护模拟实验研究［ 爆炸与冲击， 26（3）：262-268. ZHU X， HOU H L GU M B ，et al. Experimental ， study on armor protection against ballistic impact of small caliber artillery ［J］. Explosion and Shock Waves，2006，26（3）：262-268（in Chinese）. 3］梅志远，朱锡，刘燕红，等. ［ 纤维增强复合材料层合板弹道冲击研究进展［J］. 力学进展，2003，33（3）： 375-388. MEI Z Y ，ZHU X， LIU Y H ，et al. The developments of fibre-reinforced composite laminates under ballistic impact［J］. Advances in Mechanics，2003，33（3）： 375-388（in Chinese）. 4］王晓强，朱锡，梅志远. ［ 纤维增强复合材料抗侵彻研究综述［J］. 玻璃钢/复合材料，2008（5）：47-56. WANG X Q，ZHU X，MEI Z Y. The development of fi⁃ ber-reinforced composites under ballistic impact［J］. Fiber Reinforced Plastics/Composites， 2008 （5 ）： 47-56（in Chinese）. 5王晓强，朱锡. ［］舰船用钢的抗弹道冲击性能研究进展［J］. 中国造船，2010，51（1）：227-236. WANG X Q，ZHU X. Review on ballistic impact resis⁃ tance of ship building stee［l J］. Shipbuilding of China， 2010，51（1）：227-236（in Chinese）. ［6］ GREAVES L J. Failure mechanisms in GFRP armour ［R］ .UK ：Unpublished UK DRA Report，1992. ［7］ GREAVES L J. Progress in modeling the perforation of GFRP by ballistic projectiles［R］. UK ：Unpublished UK DRA Report，1994. ［8］ ZHU G Q ，GOLDSMITH W，DHARAN C K H. Pene⁃ tration of laminated Kevlar by projectiles-I. experimen⁃ tal investigation［J］. International Journal of Solids and Structures，1992，29（4）：399-420. ［9］ ZHU G Q ，GOLDSMITH W，DHARAN C K H. Pene⁃ tration of laminated Kevlar by projectiles-II. analytical model［J］. International Journal of Solids and Struc⁃ tures，1992，29（4）：421-436. ［10］ WEN H M ，REDDY T Y，REID S R，et al. Indenta⁃ tion，penetration and perforation of composite lami⁃ nate and sandwich panels under quasi-static and pro⁃
jectile loading［J］. Key Engineering Materials，1998， 141/142/143：501-552. REDDY T Y， WEN H M ，REID S R，et al. Penetra⁃ tion and perforation of composite sandwich panels by hemispherical and conical projectiles［J］. Journal of Pressure Vessel Technology， 1998， 120 （2 ）： 186-194. REID S R， WEN H M ，SODEN P D，et al. Response of single skin laminates and sandwich panels to pro⁃ jectile impact［C］//WANG S S，WILLIAMS J J，LO K H. Composite materials for offshore operation-2. ［S.l. American Bureau of Shipping， 1999： ］： 593-617. 13］覃悦，文鹤鸣，何涛. FRP ［ 卵形弹丸撞击下 层合板J］. 2007，24（2）：的侵彻和穿透［ 复合材料学报， 131-136. QIN Y， WEN H M ，HE T. Penetration and perfora⁃ tion of FRP laminates under normal impact by ogi⁃ val-nosed projectiles［J］. Acta Materiae Compositae Sinica，2007，24（2）：131-136（in Chinese）. 14］张颖军，梅志远，朱锡. FRP ［ 层合板低速冲击损伤
92 （上接第 页）
海：上海科学技术文献出版社，1982. 5］ EH HC RH ，等.气动弹性［ 道尔 ，小柯蒂斯 ，斯坎伦力学现代教程［M］.陈文俊，尹传家，译.北京：中国宇航出版社，1991. ［6］ TSIEN H S.The poincaré-lighthill-kuomethod［J］. Ad⁃ vances in Applied Mechanics，1956，4：281-349. 7］ 张效慈，司马灿，吴有生. ［ 潜艇舵低速颤振现象及其预报［J］. 船舶力学，2001，5（1）：70-72. ZHANG X C，SIMA C，WU Y S. Low-speed flutter phenomenon of submarine rudder and its prediction ［J］. Journal of Ship Mechanics，2001，5（1）：70-72 （in Chinese）. ［8］ JEWELL D A，MCCORMICK M E. Hydroelastic 望［J］. 玻璃钢/复合材料，2011特性研究现状与展（1）：52-58. ZHANG Y J， MEI Z Y ，ZHU X. Current status and prospects of study on low velocity impact damaged characterization of fiber-reinforced composite lami⁃ nates［J］. Fiber Reinforced Plastics/Composites，2011 （1）：52-58（in Chinese）. 15］谢恒，吕振华. ［ 破片侵彻纤维复合材料板的有限元数值模拟［J］.清华大学学报（自然科学版），2012， 52（1）：96-101. XIE H，LV Z H. Finite element simulation of FRP plates impacted by fragments［J］. Journal of Tsinghua University（Science & Technology），2012，52（1）： 96-101（in Chinese）. 16］李茂，朱锡，侯海量，等. ［ 冲击波和高速破片对固J］.支方板的联合作用数值模拟［ 中国舰船研究， 2015，10（6）：60-67. LIM ，ZHU X， HOU H L ，et al. Numerical simula⁃ tion of steel plates subjected to the impact of both im⁃ pact waves and fragments［J］. Chinese Journal of Ship Research，2015，10（6）：60-67（in Chinese）. in-stability of a control surface：DTMB-TR-1442［R］. Carderok，MD：David Taylor Model Basin，1961. ［9］ WRIGHT J R，COOPER J E. Introduction to aircraft aeroelasticity and loads［M］. New York：John Wiley & Sons Ltd，2008. ［10］ LIM C W ，WU B S. A new analytical approach to the Duffing-harmonic oscillator［J］. Phys Lett：A，2003， 311（4/5）：365-373. 11］肖清，谢俊超，陈东阳. ［ 舵系统的颤振计算与分析［J］. 中国舰船研究，2016，11（5）：48-54. XIAO Q， XIE J C ，CHEN D Y. Flutter calculation and analysis of rudder system［J］. Chinese Journal of Ship Research，2016，11（5）：48-54（in Chinese）.