破片群作用下复合材料层合板近场动力学损伤模拟*

2022-04-11陈志鹏马福临杨娜娜姚熊亮鞠金龙

爆炸与冲击 2022年3期

陈志鹏，马福临，杨娜娜，姚熊亮，鞠金龙

（哈尔滨工程大学船舶工程学院，黑龙江哈尔滨 150001）

反舰导弹的战斗部在爆炸后形成空爆冲击波以及破片群，舰船上层建筑防护结构不仅会受到大量破片群以及冲击波的耦合伤害，也会出现破片群单独侵彻作用一定时间之后冲击波再开始作用。破片群密集侵彻结构的破坏机理复杂多样，而纤维增强复合材料因自身力学性能复杂，其破坏机理的分析较金属靶板的破坏机理分析更困难。由于复合材料比强度与比模量较高，应力波在复合材料结构中的传播速度高于在金属结构中的传播速度，将更快地产生应力波叠加效应，从而产生更广泛的破坏范围。在应力波的叠加效应下，靶板的被撞击区域内能量密度与能量持续时间大幅增加，导致靶板在破片侵彻过程的抵抗能力大幅下降，破片的侵彻与穿透能力大幅提高，产生侵彻能力增强效应。

目前开展的实验主要集中于单个破片的侵彻作用，对多破片的实验研究较少。主要是因为开展复合材料破片侵彻的实验研究比较困难，尤其是多破片的侵彻实验，难点在于无法精确控制破片弹道轨迹和冲击速度。本文中，将从数值研究角度对复合材料破片群侵彻作用进行定量研究。针对脆性材料毁伤问题，包括裂纹扩张、分层损伤和高速侵彻等不连续问题时，采用有限元方法等传统方法仍具有局限性，为此，Silling于2000 年提出了近场动力学理论（peridynamic, PD），作为一种基于非局部作用的新兴粒子法，近20 年受到了广泛关注，尤其是在解决结构不连续性问题上有更好的应用前景。

本文中，基于键的复合材料近场动力学冲击模型，对破片群侵彻作用下纤维增强复合材料结构的损伤进行研究，分析复合材料结构破坏过程及损伤特性，探讨破片群的增强效应，并对破片速度、破片数量和破片群间距对侵彻能力增强效应的影响进行总结。

1 PD 基本理论

在近场动力学方法中，物质点x的作用范围由近场半径δ 定义，称为近场范围。近场范围内的所有物质点称为物质点x的近场域H，在近场域范围δ 之外的粒子认为不与物质点x产生相互作用力。物质点x在其近场域内物质点x的互相作用下产生变形，物质点x与x分别产生位移u和u，因此，物质点x与x之间的伸长率可以定义为：

2 复合材料层合板损伤模式与破坏过程

2.1 破片群侵彻模型

层合板模型几何尺寸为200 mm×100 mm×16 mm，采用的铺层为 [45/0/-45/90]，每一层厚度为1 mm，边界条件为四周固支，具体材料参数如表1 所示。本文中，采用自行开发的近场动力学程序对复合材料层合板进行建模，模型的离散物质点的宽度 Δ=1 mm，长度方向上设置200 个物质点，宽度方向上设置100 个物质点，厚度方向设置16 个物质点，可得物质点体积=1.0 mm；近场半径 δ =3.015Δ，时间步大小取为 Δ=2×10s，其中近场半径和时间步长已通过验算满足模型计算稳定性要求。

表1 层合板材料性能参数Table 1 Material properties of the laminate

基于上述模型，本文中研究破片数量、破片间距和破片速度对于板的损伤模式的影响效应。设定破片间距=1～2 mm，破片速度=200～1 500 m/s，破片的数量=2～9，以研究破片数量在侵彻作用中的影响效应，破片的分布形式如图1～2 所示。本文中采用球形破片，破片定义为刚性体，半径为8 mm，密度为7 850 kg/m。

图1 不同数量破片群排列示意图Fig. 1 Schematic diagram of the arrangement of fragment groups with different numbers

图2 破片群冲击工况示意图Fig. 2 Schematic diagram of the impact condition of the fragment group

2.2 损伤模式

根据近场动力学理论对于物质点损伤的定义，物质点的损伤情况可以由一个取值[0,1]的区间表示，=0 表示物质点未损伤，=1 表示物质点与其近场域内的物质点对点力全部消失。为在后处理中可以直观表达模型的损伤情况，用色条表示材料在该损伤模式下的损伤情况，如图3 所示。

图3 模型损伤程度Fig. 3 Damage degree of the model

破片群的侵彻破坏模式与破片间距、速度以及破片数量有关。由表2 可知，当破片群间距较小（=1 mm）和破片速度较低（=300 m/s）时，破片群侵彻较慢，板在初期形成鼓包变形，随后逐渐在迎弹面形成大变形和基体损伤，背弹面为纤维拉伸破坏。当速度较高时（=800 m/s），破坏模式为集团冲塞纤维层拉伸破坏，迎弹面沿面内纤维方向发生基体损伤与剪切损伤，在一定范围内产生变形，背弹面主要损伤表现为剪切损伤，损伤范围相对迎弹面更广。当破片速度达到高速时（=1 200 m/s），破坏模式依然表现为集团冲塞破坏，迎弹面及背弹面的主要损伤为剪切破坏，且背弹面的损伤范围较迎弹面大。

表2 破片间距 S=1 mm 时，在破片群侵彻下层合板的损伤模式Table 2 Damage modes of the laminate subjected to fragment group penetration when the fragment spacing S=1 mm

由表3 可知，在破片间距较大的条件下（=10 mm），当破片速度较低（=300 m/s）时，由于间距较大，层合板未出现集团冲塞破坏，但在破片冲击位置之间出现了裂纹的连接产生附加效应，在迎弹面出现大范围剪切破坏，损伤范围沿水平以及纤维方向，背弹面的剪切损伤范围相对较小。当破片速度较高（=800 m/s）时，损伤模式接近于独立穿孔的损伤模式，迎弹面冲击位置范围内产生损伤并在相互之间产生部分断裂情况，背弹面由于破片的高剩余速度产生大面积的拉伸断裂破坏。随着破片速度的提高，达到高速（=1 200 m/s）时，层合板的损伤模式主要是迎弹面的剪切穿孔，以及在冲击位置附近的部分基体损伤，由于穿透速度过快，损伤范围较小，背弹面产生纤维拉伸破坏。

表3 破片间距 S=10 mm 时，在破片群侵彻下层合板的损伤模式Table 3 Damage modes of the laminate subjected to fragment group penetration when the fragment spacing S=10 mm

由以上分析可知，单破片与破片群侵彻的复合材料层合板，其破坏模式有较大区别：单破片侵彻下，层合板的损伤模式以冲击部位的穿孔破坏以及小范围的剪切与基体损伤为主，而破片群侵彻下层合板的破坏模式主要是集团冲塞破坏、独立穿孔、大面积剪切拉伸断裂和基体损伤等。破坏程度与破片数量、间距和速度有关，其中破片间距主要影响集团冲塞破口损伤，随着间距的增大，逐渐转变为独立穿孔-附加效应。破片速度主要影响层合板中损伤的范围比例，随着速度的提高，层合板上下表面的损伤范围逐渐减小，层间损伤程度逐渐增大。

2.3 破片群侵彻过程

图4 为破片间距=1 mm、破片速度=800 m/s 的侵彻工况下，为了便于观察，以层合板的半剖模型表示层合板破坏的具体过程。由图4 可知，当=4 µs 时，破片群开始侵彻层合板，破片撞击后层合板将产生层间的压缩波以及面内应力波，压缩波在冲击位置的厚度方向产生局部压缩变形直至破坏失效，而应力波开始逐渐在各破片着弹点之间产生剪切破坏，其宏观表征体现为各破片着弹点附近形成浅坑状凹陷；当=8 µs 时，破片群之间开始产生剪切的延伸连接损伤，随着时间的延长逐渐扩大。由于间距较小，多个破片的动能会汇聚叠加，冲击位置能量密度增加形成大量集团剪切破口。当=18 µs 时，随着破片群侵彻过程中速度的不断衰减，每一层单层板的剪切破坏逐渐转变为拉伸破坏，集团剪切破口进一步形成；=38 µs 时，集团冲塞破口形成，破片群与冲塞块依靠惯性继续飞出。

图4 破片群侵彻下层合板的破坏过程（S=1 mm, v=800 m/s）Fig. 4 Failure process of the laminate subjected to fragment group penetration (S=1 mm, v=800 m/s)

图5 为破片群=10 mm、=800 m/s 的侵彻工况下，以板的半剖模型表示层合板破坏的具体过程。由图5 可知，当=4 µs 时，破片群开始侵彻层合板，冲击位置由于破片的高速冲击作用产生局部压缩变形，板发生剪切破坏，开始形成弹坑。当=8 µs 时，破片群继续侵彻作用，由于间距较大，冲击位置之间传播的应力波更加均匀分散，破坏较小，迎弹面主要在冲击位置产生独立剪切穿孔。同时，压缩波传播到背弹面产生拉伸破坏。当=16 µs 时，破片群冲击位置之间逐渐有一定损伤出现，体现为破片冲击的附加效应。当=42 µs 时，背弹面进一步拉伸，产生剪切破口，破片群穿透。

图5 破片群侵彻下层合板的破坏过程（S=10 mm, v=800 m/s）Fig. 5 Failure process of the laminate subjected to fragment group penetration (S=10 mm, v=800 m/s)

3 破片群侵彻下复合材料层合板损伤关键因素的影响

层合板在破片群侵彻过程中，多个破片分别在冲击位置造成损伤变形，并在板内产生应力波，当破片穿透层合板的时间＞/，其中为应力波传播速度，使冲击位置的能量密度与能量持续时间大大增加，提高破片的侵彻能力，即应力波的相互叠加会对破片的侵彻产生增强效应。当破片穿透层合板的时间＜/，破片穿透层合板的时间内，应力波未开始相互叠加，此时不会产生破片群的侵彻增强效应，侵彻能力接近于单破片侵彻能力。由此可知，破片间距和破片速度是破片群侵彻能力的主要影响因素，前者影响应力波的叠加发生时间，后者影响破片穿透的过程，两者共同决定了破片群的侵彻能力。此外，破片的数量也是影响破片群侵彻能力的重要因素，当破片数量过少，叠加效应不足以产生明显影响。

采用控制单一变量方法，开展层合板破片群侵彻关键因素的影响规律研究。为将破片群侵彻能力量化分析，采用破片群穿透后剩余速度最大的破片表征该工况下破片群的侵彻能力。而根据图3 中破片群的排列方式以及应力波叠加原理，中心破片将会受到主要的加强作用，其侵彻能力的增强效应最大，所以可以直接采用中心破片的侵彻状态表征相应工况下破片群的侵彻能力。

3.1 破片群数量的影响

保持破片群速度=800 m/s 和破片间距=1 mm 不变，仅改变破片数量获得不同工况下破片群侵彻的剩余速度，如图6 所示。单破片侵彻时，破片剩余速度为472 m/s。破片数量为2 枚时，两破片的应力波叠加效应开始产生一定影响，但由于破片数量较少，剩余速度为541 m/s。随着破片数量的增加，应力波的叠加效应对于中心破片的影响不断增强，破片剩余速度逐渐增高。

图6 破片数量对破片群侵彻能力的影响Fig. 6 Influence of the fragment number on the penetration ability of the fragment group

当破片数量增加至5 枚以上时，速度增高幅度相对较小，当破片数量到达9 枚时，破片数量对于破片群侵彻能力的影响已经微乎其微。

3.2 破片间距的影响

根据3.1 节的分析结果，固定破片数量=9以及破片速度=800 m/s，研究破片间距对破片群侵彻能力的影响。图7 为不同破片间距下破片群侵彻剩余速度与单破片侵彻剩余速度对比。当=0 mm时（设置的极端工况，每个破片都紧密地，无间隔地聚集在一起），破片群侵彻剩余速度比单破片侵彻剩余速度高50.4%，随着破片间距的增大，破片群剩余速度不断降低。当=15 mm 时，破片群与单破片侵彻剩余速度基本相等，可以认为破片群效应已退化为独立破片侵彻。这种破片群效应出现的主要原因在于，应力波的叠加所需时间随着破片冲击位置之间距离的增大而不断增长，直至出现破片已经开始穿透层合板而应力波叠加效应仍未作用，导致破片间距对破片群效应的不断减弱直至退化消失。

图7 破片间距对破片群侵彻能力的影响Fig. 7 Influence of the fragment spacing on the penetration ability of the fragment group

图8 为破片速度=800 m/s 和破片间距=5 mm 的破片群侵彻下，层合板的损伤过程。由图8 可知，当=4 µs 时，破片群开始侵彻层合板，冲击位置开始具有法向速度并产生弹坑；当=12 µs 时，随着应力波的传播，迎弹面冲击位置之间开始产生变形损伤，背弹面开始产生局部断裂；当=38 µs 时，由于层间损伤导致背弹面纤维层的拉伸损伤断裂，破口形成，损伤失效材料与破片一起飞出。

图8 破片群侵彻下层合板的破坏过程（S=5 mm, v=800 m/s）Fig. 8 Failure process of the laminate subjected to fragment group penetration (S=5 mm, v=800 m/s)

图9 为破片速度=800 m/s 和破片间距=20 mm 的破片群侵彻下，层合板的损伤过程。由图9 可知，当=4 µs 时，破片群着靶，层合板冲击位置获得法向速度产生弹坑，应力波开始向四周传播；当=12 µs 时，由于破片间距较大，应力波叠加效应已经可以忽略，破片在背弹面形成局部变形拉伸；当=42 µs 时，破片完全穿透层合板，在背弹面形成大范围剪切拉伸损伤，非冲击区域也产生一定的剪切拉伸损伤。

图9 破片群侵彻下层合板的破坏过程（S=20 mm, v=800 m/s）Fig. 9 Failure process of the laminate subjected to fragment group penetration (S=20 mm, v=800 m/s)

综上所述，破片间距较小时应力波叠加效应明显，应力波在破片冲击位置以外造成损伤，尤其是层间损伤情况，使得连接处产生拉伸断裂破坏进而失去抗冲击能力。当破片间距较大时，应力波叠加效应影响较小，非冲击区域的位移较小，不足以让背弹面提前产生局部拉伸破坏，背弹面仍具有一定的抗冲击能力，类似于单破片侵彻情况。

3.3 破片冲击速度的影响

固定影响因素破片间距=3 mm 和破片数量=9，研究破片速度=300, 800, 1200, 1500 m/s 等4 种工况下破片速度与破片群效应的关系。图10 为该破片群与相同速度单破片工况下剩余速度对比情况。由图10 可知，当=300 m/s 时，破片群剩余速度是单破片的1.78 倍。随着破片速度的升高，破片群与单破片的剩余速度差距不断减小，即群效应不断减弱。当=1 200 m/s 时，破片群的侵彻增强效应退化消失为单破片侵彻情况。这是由于在破片间距不变条件下，穿透时间随着破片速度升高而缩短，应力波叠加效应影响时间缩短，即＜/，因此破片群侵彻能力增强效应随破片速度的升高而不断退化直至消失。产生较大的变形，层间产生剪切破坏；当=12 µs 时，冲击位置之间产生大量剪切损伤，背弹面受到层间压缩波的作用产生剪切损伤，破片最下端已到达层合板底部位置；当=22 µs 时，破片完全穿透层合板。由于破片的高速侵彻，应力波叠加效应不足以影响破片群侵彻能力，损伤模式已退化为单破片模式。