•  栏目 研究性文章

三层靶板耗散子弹侵彻自适应FEM-SPH耦合算法仿真

作者

摘要

复合装甲因其显著的性能优势,在多个领域得到广泛应用。通过爆炸焊接技术,能够将铝合金与陶瓷材料牢固结合并制造反弹装甲。本研究利用ANSYS/LS-DYNA有限元程序对三层3mm厚材料(外侧两层为7039铝合金,中间层为氧化铝陶瓷)进行数值模拟,探讨该装甲在刚体侵彻后的动态响应。仿真采用自适应FEM-SPH耦合算法,特别适用于高速冲击工况下引起的网格破碎、飞溅等情况,具有良好的收敛性。研究的主要目的是评估爆炸焊接工艺对装甲防护性能的影响,具体目标包括:1) 探讨不同内应力条件下装甲的抗穿透能力;2) 分析子弹动能消耗和速度变化规律;3) 验证自适应FEM-SPH耦合算法在处理复杂材料模型中的有效性。仿真结果显示,随着内应力数值的增加,子弹的动能消耗显著增加,最大可达18.99%。击穿第三层靶板时子弹的速度差异为29.64%,表明内应力对子弹动能和速度有显著影响。此外,装甲的最大变形量从5.22mm减少到3.29mm,进一步证明了内应力对装甲变形的重要影响。自适应FEM-SPH耦合算法能够有效模拟氧化铝陶瓷的破坏过程,并且在考虑爆炸焊接产生的内应力后,仿真模型的精度得到了显著提高。

引言

复合装甲由于其出色的性能优势,在现代防护装备领域占据了不可或缺的地位。它由金属、陶瓷和纤维复合材料等构成,这些材料的巧妙结合赋予了复合装甲多项显著特性。

首先,高强度和高韧性是复合装甲最突出的特点之一。这种装甲能够有效抵御高速弹丸以及爆炸产生的破片冲击,确保在极端条件下也能提供可靠的保护。其次,轻量化设计是复合装甲的另一大亮点。通过采用铝合金、陶瓷和纤维复合材料等轻质但强度极高的材料,不仅能够在重量上大幅减轻,同时还能保持甚至提升防护效能,这对于提高车辆的机动性和燃油效率至关重要。此外,良好的能量吸收能力也是复合装甲的一大特色。其设计能够有效地分散并吸收冲击能量,比如当受到攻击时,陶瓷层会破碎以消耗部分能量,而金属层则负责阻止任何可能穿透的碎片或弹头。

更进一步,多层结构设计使得复合装甲的功能更加多样化。每一层都具有特定的功能,如有的层专门用于分散能量,有的层则是为了增强整体结构的坚固性。这种设计不仅能提高对多种威胁的适应性,还大大增强了整体的防护性能。因此,复合装甲被广泛应用于军事装备中,例如坦克和装甲车,它们可以有效抵挡子弹和炮弹碎片的攻击,极大提升了作战车辆的生存能力。而在警用和民用市场,复合装甲也被用于制造个人防护装备、防爆盾牌和防弹车辆等,为公共安全和个人防护提供了强有力的支持。总之,复合装甲凭借其卓越的性能,在不同领域的应用前景广阔,成为现代防护技术的重要组成部分。

对于ALE算法[1]虽然能够较好的处理由于计算过程中网格畸变引起的不收敛的问题,但是在针对陶瓷类脆性材料计算时,会出现网格极端变型和断裂的情况,就需要采用新的算法。鉴于本次研究对象为采用爆炸焊接工艺[2]制造的三层复合材料,中间层为氧化铝陶瓷材料[3,4],仿真将采用自适应FEM-SPH算法[5],自适应FEM-SPH耦合算法可以根据问题的复杂程度自动调整网格和粒子分布,提高计算效率和精度。尤其是当模型中出现极端变形和断裂问题。在装甲中采用陶瓷[6-11]等脆性材料的情况,子弹侵彻仿真采用该算法具有更好的效果。

1 计算模型

该模型设置为三层装甲,单层的尺寸为40mm×40mm×3mm,子弹的材料为4340钢材,以500 m/s的速度垂直侵彻三层金属靶板的中心位置。外侧两层板的材料采用7039铝合金,中间板的材料采用氧化铝陶瓷材料。鉴于该装甲为爆炸焊接加工制造,单层之间没有间隙,但在模型中建立接触对,在模拟中可以进行滑动。

Figure 1.

Figure 2.

图1模型几何参数示意图

爆炸焊接,通过爆炸产生的冲击能量转换为装甲版之间的内应力将材料间进行牢固结合,爆炸焊接工艺原理图见下图2所示。不同的爆炸焊接速度,则产生不同的爆炸焊接内应力。为探究不同爆炸焊接速度下,不同内应力(图3所示)对于模型结果的影响。由于内应力的测量较为困难,本次研究分别假设100MPa、200MPa、300MPa、400MPa和500MPa进行仿真研究,来探究装甲间不同内应力对于仿真结果的影响。

Figure 3.

图2 爆炸焊接工艺过程示意图

Figure 4.

图3 仿真模型示意图

Figure 5.

本次仿真采用Johnson Cook参数如下表1所示。

表1 模型中材料Johnson Cook参数表

材料 A/MPa B/MPa C n m
4340 钢板 792 510 0.26 0.014 1.03
7039 铝合金板 337 343 0.41 0.010 1.00
Table 1.

Johnson-Cook (JC) 本构模型是一种用于描述材料在高应变率、高温条件下塑性行为的经验模型。它特别适用于模拟金属材料在高速冲击或爆炸等极端条件下的响应。JC模型通过一组特定的参数来描述材料的流动应力随应变、应变率和温度的变化。以下是JC模型中主要参数的含义:

A-初始屈服应力。这是材料在零应变率和室温下的初始屈服强度

B-硬化常数。这个参数与材料的塑性硬化有关,反映了随着塑性应变增加,材料流动应力增大的程度。

n-硬化指数。这是一个无量纲的参数,用来描述材料硬化速率的快慢。n的值决定了流动应力随塑性应变增长的速度。

C-应变率敏感系数。这也是一个无量纲的参数,用来表示流动应力对应变率变化的敏感程度。较高的C值表明材料的流动应力随应变率的增加而显著增大。

m-热软化指数。同样是无量纲的参数,它描述了温度对材料流动应力的影响。随着温度的升高,材料的流动应力会下降,m值越大,这种下降的趋势越明显。

内层氧化铝陶瓷材料,当陶瓷出现微裂纹损伤破坏时,在压缩载荷的作用下呈现出渐进损伤状态,而这种大压强、高应变率条件下的脆性材料失效建模普遍采用损伤演化的Johnson−Holmquist(JH−2)材料模型,其中JH−2模型的强度模型为:

Figure 6.

损伤材料的归一化等效应力为:

Figure 7.

Figure 8.

完全断裂损伤材料的归一化等效应力为:

表2 模型中陶瓷Johnson−Holmquist参数

ρ/(k·gm -3 G/ GPa A B C M N β
3350.0 84.6 0.887 0.29 0.00377 0.53 0.748 1.0 1.0
/GPa HEL/GPa PHEL/ GPa D1 D2 K1/ GPa K2/ GPa K3/ GPa
0.01 3.84 1.84 0.01 -0.325 1.063 0 0
Table 2.

2 结果与讨论

在现代装甲设计中,相对于中间有间隙的双层装甲系统,通过爆炸焊接工艺制造的三层装甲展示了更高的防护性能和结构稳定性。这种三层装甲通常由外层金属、中间陶瓷层以及内层金属构成,其中间层材料如氧化铝陶瓷等脆性材料在抵御高速弹丸侵彻方面表现出了优异的特性。为了准确模拟这类复杂材料的行为,仿真模型的选择至关重要。在Ansys Lsdyna中,Johnson-Holmquist (JH-2) 材料模型因其能够有效描述陶瓷材料在高应变率下的破坏行为而被选用,其适用性在这个背景下显得尤为突出。

爆炸焊接是一种独特的工艺,它利用炸药爆炸产生的冲击动能将不同材料牢固地结合在一起,从而在装甲各层之间形成稳定的内应力。这些内应力不仅增强了装甲的整体强度,还提高了其抗冲击能力。然而,由于实际测试条件的限制,精确测量这些内应力的具体数值非常困难。因此,在本次研究中,我们采用Load_Segment_Set卡片技术对不同装甲间的应力进行施加,以探索这些内应力对子弹速度及靶板耗散动能的影响。

通过一系列仿真实验发现,随着装甲间内应力的增加,子弹的速度显著降低,同时靶板耗散的动能显著增加。这表明,装甲在承受冲击时能够更有效地吸收和分散能量,从而提高防弹性能。具体而言,当内应力从100MPa增加到500MPa时,子弹动能的耗散比例分别增加了29.64%和18.99%,显示出内应力对提升装甲防护能力的重要作用。仿真结果进一步证实了内应力在增强装甲性能方面的关键作用。较高的内应力不仅减少了子弹穿透的可能性,还大大降低了子弹对后方结构的潜在威胁。因此,在设计和评估此类装甲系统时,必须仔细考虑并优化内应力参数,以确保装甲具备最佳的防护效果。在涉及爆炸焊接工艺的三层装甲仿真中,合理设置和控制装甲间的内应力对于提高仿真精度和装甲防护性能的重要性。通过精确模拟内应力对子弹速度和靶板耗散动能的影响,研究人员能够更好地理解装甲材料在极端条件下的响应机制,并为开发更加高效、可靠的防护装备提供坚实的理论基础和技术支持。这一研究成果不仅推动了装甲材料科学的发展,也为未来防护技术的进步提供了宝贵的参考。

Figure 9.

Figure 10.

Figure 11.

Figure 12.

Figure 13.

Figure 14.

图4 击穿装甲后不同内应力下位移结果图

根据图4中的位移仿真结果图,可以观察到装甲在被击穿过程中,网格发生破损,由于陶瓷材料的本身延展性较差,在击穿过程中以网格以碎裂的方式破坏;通过击穿后的装甲的Y方向最大变形量的数值可以发现,随着内应力数值的增加,装甲变形量从最大的5.22mm逐步减少到3.29mm。说明内应力对于装甲的变形有很大影响,内应力值越大,装甲变形量越小。

表3 不同板间内应力下的子弹残余速度和耗散动能百分比

应力 / MPa 弹体初速度 (m/s) 击穿第一层板 击穿第二层板 击穿第三层板 靶板弹体耗散动能 %
历时 / μs 残余速度 /m · s -1 历时/ μs 残余速度 /m · s -1 历时/ μs 残余速度 /m·s -1
100 500 10.0 445.0 19.9 389.0 39.9 307.0 62.2
200 10.0 442.0 19.9 378.0 39.9 279.0 68.8
300 10.0 440.0 19.9 371.0 39.9 261.0 72.8
40 0 10.0 437.0 19.9 362.0 39.9 236.0 77.7
500 10.0 433.0 19.9 356.0 39.9 216.0 81.2
Table 3.

结合表3和图5的结果,可以发现,子弹击穿各层装甲的时间差异不大,但是随着内应力的提升,子弹在击穿各装甲后的速度有较为明显的降低,在击穿第三层装甲后,弹体速度由307.0m/s降低到216.0m/s且靶板弹体耗散动能也从62.2%降低到81.2%。

3 结论

在Ansys Lsdyna中,自适性FEM-SPH(有限元法-光滑粒子流体动力学)耦合算法被广泛应用以模拟复杂材料的动态行为,尤其是在多层装甲的破坏过程研究中显示出独特的优势。本研究聚焦于基于Johnson Cook模型的金属材料和Johnson-Holmquist (JH-2) 模型的陶瓷类材料,尤其是氧化铝陶瓷在三层装甲结构中的破坏特性。

鉴于爆炸焊接工艺的特点,该工艺过程中产生的内应力对仿真结果有着显著影响。因此,在构建仿真模型时,特别考虑了这一因素,通过引入内应力参数来提高模型的准确性。具体而言,通过调整不同的内应力数值进行了一系列子弹侵彻仿真实验。实验结果显示,尽管击穿各层装甲所需的时间差异不大,但装甲的最大变形量从5.22mm减少到了3.29mm,表明内应力对装甲的变形有显著抑制作用。进一步分析发现,在100MPa和500MPa的内应力条件下,子弹的速度以及靶板耗散的动能存在明显差异。随着内应力的增加,速度和动能损耗也随之增大,分别达到了29.64%和18.99%的变化率。这说明内应力不仅影响装甲的变形,还显著改变了子弹的侵彻效率及其能量损失模式。

由此得出结论,在进行涉及爆炸焊接工艺的侵彻仿真时,必须合理设定内应力参数。忽略或错误估计这些参数可能导致仿真结果与实际情况产生较大偏差,从而影响研究结论的有效性和可靠性。研究人员需要充分认识到装甲间内应力对于仿真结果的巨大影响,并在实验设计阶段就对其进行精确控制,以确保最终获得的数据具有较高的准确性和实用性。这不仅有助于更深入地理解装甲材料在极端条件下的响应机制,也为优化装甲设计提供了重要的理论依据和技术支持。

参考文献

[1] 常向阳,王自力. 爆炸成型弹丸侵彻钢靶的ALE 算法[J]. 解放军理工大学学报:自然科学版,2004,5(3):70-73.

[2]朱桂春,刘中枢,梁国峰,等.不锈钢06Cr18Ni11Ti与铸钢20Mn的爆炸焊接研究[J/OL].爆破,1-16[2024-11-03].

[3]赵宇.基于SPH方法爆炸焊接板材的数值模拟研究[D].安徽理工大学,2024.

[4]张文斌,杨海娟,刘翠荣,等.爆炸焊接N6/45#复合板的数值模拟与实验分析(英文)[J].稀有金属材料与工程,2024,53(06):1592-1600.

[5]杜忠华,赵国志,王晓鸣,等.双层陶瓷复合靶板抗弹性的研究[J].航空学报,2002,(02):147-150.

[6]秦溶蔓,朱波,乔琨,等.复合结构碳纤维防弹板的防弹性能仿真[J].工程科学学报,2021,43(10):1346-1354.

[7]申志强,蒋志刚,曾首义,等.穿甲子弹偏心入射陶瓷/钢复合靶板试验[J].国防科技大学学报,2007,(06):34-38.

[8]潘林,陈叶娣,黄敏高,等.基于CAE技术的便携式无人机干扰设备壳体注射模开发及应用[J].模具制造,2025,25(03):20-22.

[9]张见升,景建斌,孙浩,等.爆炸荷载下预应力T型梁毁伤效应数值仿真分析[J].兵工学报,2024,45(S2):193-198.

[10]陈斌,李亚,张绍伦,等.爆炸焊接制备C276/Q345R复合板及其热处理研究[J].铁合金,2025,56(01):29-35.DOI:10.16122/j.cnki.issn1001-1943.2025.01.007.

[11]周嘉楠,罗宁,梁汉良,等.TC4/Al 6063/Al 7075爆炸焊接复合板的微观结构分析:多尺度模拟与实验研究(英文)[J].稀有金属材料与工程,2025,54(01):27-38.

基金项目:浙江省教育厅一般科研项目(Y202351814)、浙江工商职业技术学院科研年度项目(KYND202203、KYND202420)、浙江工商职业技术学院科研创新团队项目(KYTD202101)

作者简介:杨少增(1991-),男,汉族,山东德州人,讲师,研究方向为模具设计、电火花加工、爆炸焊接等方面。杨少增 男 汉族 籍贯:山东省德州市 职称:讲师 研究方向:模具设计、电火花加工、爆炸焊接等方面

如何引用

杨少增, 马长越, 刘剑, & 蔡永根. (2025). 三层靶板耗散子弹侵彻自适应FEM-SPH耦合算法仿真. 工程与技术创新, 1(1), 79–84. 取读于 从 http://8.219.207.158/index.php/JETI/article/view/379
下载引用

参考

常向阳,王自力. 爆炸成型弹丸侵彻钢靶的ALE 算法[J]. 解放军理工大学学报:自然科学版,2004,5(3):70-73.

朱桂春,刘中枢,梁国峰,等.不锈钢06Cr18Ni11Ti与铸钢20Mn的爆炸焊接研究[J/OL].爆破,1-16[2024-11-03].

赵宇.基于SPH方法爆炸焊接板材的数值模拟研究[D].安徽理工大学,2024.

张文斌,杨海娟,刘翠荣,等.爆炸焊接N6/45#复合板的数值模拟与实验分析(英文)[J].稀有金属材料与工程,2024,53(06):1592-1600.

杜忠华,赵国志,王晓鸣,等.双层陶瓷复合靶板抗弹性的研究[J].航空学报,2002,(02):147-150.

秦溶蔓,朱波,乔琨,等.复合结构碳纤维防弹板的防弹性能仿真[J].工程科学学报,2021,43(10):1346-1354.

申志强,蒋志刚,曾首义,等.穿甲子弹偏心入射陶瓷/钢复合靶板试验[J].国防科技大学学报,2007,(06):34-38.

潘林,陈叶娣,黄敏高,等.基于CAE技术的便携式无人机干扰设备壳体注射模开发及应用[J].模具制造,2025,25(03):20-22.

张见升,景建斌,孙浩,等.爆炸荷载下预应力T型梁毁伤效应数值仿真分析[J].兵工学报,2024,45(S2):193-198.

陈斌,李亚,张绍伦,等.爆炸焊接制备C276/Q345R复合板及其热处理研究[J].铁合金,2025,56(01):29-35.DOI:10.16122/j.cnki.issn1001-1943.2025.01.007.

周嘉楠,罗宁,梁汉良,等.TC4/Al 6063/Al 7075爆炸焊接复合板的微观结构分析:多尺度模拟与实验研究(英文)[J].稀有金属材料与工程,2025,54(01):27-38.

作者简介