Influence of Different Material Constants Fitting Method on Predicting Warhead Impacting Metal Targets
-
摘要: 弹靶侵彻仿真中材料参数对计算结果有着至关重要的影响。为寻求一套适用于弹靶侵彻仿真计算的材料参数拟合方法,借助前期开展的靶板材料动态力学性能试验、靶板材料断裂试验,通过不同拟合方法依次得到不同的JC本构模型及失效模型参数,依据试验建立有限元计算模型,将数值计算结果与试验结果进行对比。结果表明:(1)对于同一材料的力学性能试验,采用不同的拟合方法可得到不同的JC本构、JC失效参数,二者会对弹靶仿真结果造成一定影响;(2)在不考虑温度软化项的前提下,采用高应变率作为参考应变率进行拟合能更加准确地表征材料在高应变率下的应力-应变关系,更加适用于弹靶侵彻强瞬态、高应变率作用过程仿真;(3)对于同一JC本构模型,采用平均应力三轴度拟合的JC失效模型较采用初始应力三轴度拟合的JC失效模型所得战斗部剩余速度计算结果偏小,仅采用拉伸试件结果拟合的JC失效模型较采用扭转、拉伸试件结果拟合的JC失效模型所得战斗部剩余速度计算结果偏小。Abstract: Material constants played a significant role in the numerical prediction of warhead impacting metal targets, in order to obtain a set of valuable material constants fitting methods, JC strength model and JC failure model constants were acquired by different fitting method based on material mechanical properties and failure experiments. Simulation results were compared with experiment outcome in which different nose warhead penetrated different thickness metal plates, and it indicated: (1) simulation results could be distinct using JC strength model and JC failure mode constants from different fitting methods for the same material mechanical property tests; (2) JC strength model constants fitting by stress-strain curves using high reference strain rate were suitable for the numerical predicting warhead penetrating metal plates; (3) JC failure model constants fitting by initial stress triaxiality or average stress triaxiality had few influence on calculation results. This research could provide assistance for the material constants fitting method in numerical simulation.
-
Key words:
- material constants /
- constants fitting /
- JC strength model /
- JC failure model /
- stress triaxiality /
- penetration
-
聚脲是由异氰酸酯组分和氨基化合物快速交联加聚反应而形成的一类高分子材料,因具有优异的力学性能[1-3],受到国内外专家和学者的广泛关注,并开展了爆炸防护研究。相关研究[4-6]表明,聚脲涂层能够减弱钢板在爆炸荷载作用下的损坏程度,增强整体结构的抗爆能力。而对于混凝土结构,无论是爆破拆除,还是反恐防爆,减少混凝土破碎导致的损害是目前研究的热点。在聚脲加固钢筋混凝土结构方面,美国空军实验室研究聚合物加固墙体时发现,相较于碳纤维、芳纶纤维等其他众多高分子聚合物材料,聚脲在墙体结构爆炸防护中表现出更好的抗爆效果[7]。Parniani等[8]在此基础上研究发现,聚脲能够提高钢筋混凝土梁的抗弯能力和延性。Raman等[9]和Ghaderi等[10]通过数值模拟发现,聚脲涂层可以提升钢筋混凝土板的位移控制和能量耗散,降低结构的应变能及动能,通过降低位移和能量吸收、减弱爆炸能量传递来提高抗爆能力。Iqbal等[11]通过聚脲涂覆混凝土爆炸实验,发现聚脲能够包覆混凝土爆炸后产生的碎片,减少碎片对周围人员及物体造成的伤害。王军国[12]、尚伟等[13]对聚脲加固墙体进行了爆炸试验,发现有聚脲涂层加固的墙体具有良好的抗爆能力,且背爆面涂覆墙体在爆炸作用下没有产生爆炸碎片,表明聚脲的包覆效果良好。聚脲涂层厚度对结构的抗爆能力有着巨大影响,为了深入研究聚脲涂层厚度对墙体抗爆能力的影响,Tao等[14]开展了不同工况下聚脲加固墙体的爆炸实验,发现喷涂3 mm厚聚脲的墙体在爆炸荷载下损伤较小,只有轻微的凹陷变形,而喷涂4 mm厚聚脲的墙体没有破坏,显示出良好的防爆性能。田力等[15]、Gu等[16]运用ANASYS/LS-DYNA软件模拟了在爆炸荷载作用下不同厚度聚脲涂覆墙体的破坏情况,结果表明,随着聚脲涂层厚度的增加,墙体中心点位移明显减小,聚脲加固后墙体的抗爆能力提高。
目前,人们已经开展了聚脲涂层的抗爆研究,并取得了良好的效果,但是对于聚脲涂层抗爆机理的研究尚不够完善,需要对爆炸过程中聚脲的力学性能及爆炸后涂层的微观及宏观破坏形式进行深入研究,以确定两者之间的关联。为此,本研究选择青岛理工大学功能材料研究所自主研发的抗爆型T26聚脲,分析不同应变率下T26聚脲的力学性能;通过广角X射线衍射(wide-angle X-ray diffraction,WXRD)实验,研究材料的微观结构,并利用差示扫描量热(differential scanning calorimetry,DSC)实验分析其热性能,研究T26聚脲在高温条件下的形态及其特点;在此基础上,对喷涂T26聚脲前后的钢筋混凝土板进行接触爆炸实验,通过1.4 kg当量TNT下混凝土板的宏观破坏形貌,分析其破坏模式,并借助电子显微镜观察涂层的微观形态,结合其热性能结果,探讨涂层的爆炸防护规律,分析其抗爆机理。
1. 聚脲涂层材料制备与性能
1.1 涂层试样制备
T26聚脲是由青岛理工大学功能材料研究所针对抗爆要求专门研发的高强抗爆型聚脲。喷涂设备包括PHX-40主机和AP-2喷枪(美国聚氨酯机械公司生产),喷雾温度为65 ℃,压力为17.24 MPa。由于材料力学性能是缓慢提升的,因此实验前需要将试样在常温下养护7 d,保证其力学强度稳定。
1.2 涂层性能实验
1.2.1 涂层力学性能
根据ASTM D412《硫化橡胶和热塑性弹性体拉伸试验方法》,将T26聚脲试样分别裁成无割口直角形与哑铃形。裁取哑铃形试样中间25 mm长区域作为测量区,在测量区内取两端及中间3点1、2、3(见图1(a))进行厚度测量,取平均值并输入计算机端。将待测哑铃形试样两端用万能试验机夹具夹住,使受力方向与试样测量区间保持重合,夹具间距离设置为70 mm,调零后运行。拉伸试样直至断裂,得到其拉伸强度和断裂伸长率。取直角形样片中心及两边3点进行厚度测试(见图1(b)),求平均值,进行撕裂强度测试。
1.2.2 涂层WXRD实验
选取尺寸为2 cm×2 cm的聚脲样片,将其固定在测角仪的载玻片上,采用RINT2000 vertical goniometer测角仪,对样片表面进行测试,得到衍射峰,并计算结晶度。扫描速度为10°/min,精度为 ± 0.02°,扫描范围为5°~80°。
1.2.3 涂层热性能
称取T26聚脲样品12.2 mg,利用NETZSCH DSC 204 F1差示扫描量热仪(见图2)对样品进行测量。为减小实验误差,在氮气条件下进行DSC测试,升温速率为10 ℃/min,温度范围为−100~800 ℃。
2. 涂层抗爆特性预估
2.1 涂层密度和力学性能
T26聚脲是传统聚脲配方优化调整后的高强抗爆型聚脲材料,T26聚脲的密度在1 g/cm3以下(详见表1)。对T26聚脲进行准静态拉伸实验,其应力-应变曲线及材料试件如图3(a)所示;对于该曲线的弹性阶段,采用一次函数y=0.493+79.729x进行拟合,得到材料的应力-应变曲线拟合情况,如图3(b)所示。所测得的T26聚脲的基本力学性能如表1所示。
表 1 T26聚脲与其他聚脲材料的物理性能对比Table 1. Physical properties comparison of T26 polyurea with other polyureasMaterial Density/
(g·cm−3)Elastic modulus/
MPaTensile strength/
MPaPeel strength/
(N·mm−1)Breaking elongation/
%T26 0.977 84.05 25.4 75.5 451.9 T1[13] 1.020 19.2 391.7 T2[17] 1.02−1.07 18.0 350.0 Note: (1) T1 is energy absorbing polyurea for strengthening masonry wall explosion test;
(2) T2 is common polyurea for explosion test of reinforced autoclaved aerated concrete plate.由准静态拉伸实验结果可知,T26聚脲的拉伸强度达到25.4 MPa,断裂伸长率为451.9%,相较于T1、T2两种聚脲,其抗拉强度和断裂伸长率均大幅增加。炸药产生的爆炸荷载可看作均布荷载作用到板表面,由于混凝土板的抗压强度远大于其抗拉强度,需要涂层承受较大的拉力作用,T26涂层的抗拉强度高,可以起到更好的保护作用;同时T26聚脲具有较高的断裂伸长率,说明T26聚脲涂层具有良好的柔性,可以更好地包裹混凝土板,降低碎片飞溅的速度。
2.2 涂层WXRD分析
图4为T26聚脲的WXRD图谱及分峰处理结果。从图4(a)中可以发现,在2
=20.52°附近出现较强的软段非晶漫散射峰,材料链段更趋向于规整有序排列或空间折叠结构,表明材料中的脲键作为分子链中的硬段部分,达到了远程有序,产生结晶。WXRD分峰处理结果如图4(b)所示,具体数据列于表2。T26聚脲微晶区位于38°附近,无定形区位于20°附近,结晶度(结晶峰面积与总的峰面积的比值)较高,达到24.11%,说明材料分子链中氢键较多。氢键具有易断裂、易形成的特点,并且氢键的作用力大于原子间作用力,小于化学键的作用力。随着氢键的断裂与形成,材料表现出优异的吸能效率,从宏观上看,有序规整的分子结构更有利于材料形成优异的力学性能。表 2 T26聚脲的WXRD分峰处理结果Table 2. WXRD peak separation results of T26 polyureaSpecimen Peak position/(°) Peak area Crystallinity/% Amorphous region Microcrystalline domain Amorphous region Microcrystalline domain T26 20.52 38.16 24042.85 7638.57 24.11 2.3 涂层热性能分析
图5为T26聚脲在−100~80 ℃范围内的DSC曲线。从图5中可以得到,材料的软段玻璃化转变温度(Tgs)为−44.9 ℃,硬段玻璃化转变温度(Tgh)为36.5 ℃,说明T26聚脲分子链为一定的微相分离结构,软段与硬段之间排列有序。在−100~80 ℃温度范围内,T26聚脲材料没有出现明显的吸热和放热峰,说明在较高的温度环境下未发生较大的变化,对高温环境的适应性较强。
由于爆炸产生的热量极高,因此需要对T26聚脲在更高的温度下进行测试,为此选择30~700 ℃的温度范围进行DSC实验,如图6所示。从图6中可以看出,随着温度的升高,375 ℃附近出现结晶熔融吸热峰,说明T26聚脲中存在结晶,与WXRD实验结果一致,因此材料可以在高温下通过硬段域内氢键的断裂吸收能量,产生结晶熔融吸热峰。材料分子链中的软段决定着材料的柔性,硬段决定着材料自身的强度。在爆炸过程中,炸药会瞬间释放能量,产生高温高热,T26聚脲在高温下会发生结晶熔融,通过涂层的结晶熔融吸收能量。同时,在高温环境下,T26聚脲材料的自由体积开始膨胀,分子链段运动变得活跃并逐渐舒展,材料更加柔软。在爆炸中,良好的柔韧性使T26聚脲能够包覆混凝土板表面,黏结爆炸碎片,减少爆炸碎片对周围环境造成的危害。
3. 抗爆实验验证
3.1 聚脲涂覆钢筋混凝土板的爆炸实验
3.1.1 实验方案及试样喷涂
按照GB50010—2010《混凝土结构设计规范》,混凝土强度等级为C40,板采用双层配筋,钢筋间距为150 mm,直径为10 mm,保护层厚度为15 mm,钢筋混凝土板的尺寸为1000 mm×1000 mm×300 mm,如图7所示。
T26聚脲涂层与混凝土基材良好的附着至关重要。在喷涂之前,需要将环氧改性聚氨酯底漆涂抹在基材表面上,这是由于:(1) 底漆的封闭能力较好,可以降低基材中的针眼缺陷;(2) 底漆可以将残留在基材表面的灰尘牢固地粘在基材表面,避免因灰尘导致的涂层缺陷;(3) 底漆的渗透能力较强,能同时与混凝土和聚脲涂层紧密黏结,增强涂层的附着力。底漆选用青岛理工大学功能材料研究所自主研发的环氧改性聚氨酯底漆,附着力强度可达到5 MPa。
利用喷涂设备将聚脲喷涂到钢筋混凝土板表面,采用全面喷涂方式,涂层厚度10 mm。将经T26聚脲喷涂处理后的钢筋混凝土板养护7 d后,安置在约束固定架上,背爆面侧挖深度为0.5 m、宽度为0.7 m的爆炸坑,炸药放置在钢筋混凝土板的上表面,进行有无涂层钢筋混凝土板的爆炸实验,TNT药量为1.4 kg,如图8所示。
3.1.2 爆炸后钢筋混凝土宏观形貌
根据实验方案(见表3),比较有无涂层钢筋混凝土板爆炸后迎爆面与背爆面的损伤形式。设置工况3为对照组,分析钢筋混凝土板迎爆面与背爆面T26涂层的防护效果。
表 3 宏观形貌分析方法Table 3. Analysis method of macroscopic morphologyExp. condition Coating Coating thickness/mm Failure plane Macroscopic morphology comparison 1 None 0 Contact blast face Failure mode and diameter Back burst face Failure mode and diameter 2 T26 10 Contact blast face Failure mode and diameter Back burst face Failure mode and diameter 3.1.3 爆炸后涂层的微观形貌
选取爆炸后不同爆炸区域的T26涂层,利用JSM-7500F扫描电子显微镜对爆炸后样品进行放大4000倍的微观形貌观察。该仪器的分辨率为1.0 nm(15 kV)和1.4 nm(1 kV),加速电压为0.1~30 kV,电子枪使用钨丝灯,放大倍数为2000。
将爆炸破坏后的T26聚脲涂层划分为3个区域,分别为断面区域Ⅰ、爆炸源中心区域Ⅱ及涂层边缘区域Ⅲ,如图9所示。分别裁取3个区域涂层,借助扫描电子显微镜进行微观形貌观察,研究爆炸冲击波在涂层爆炸核心区及边缘区域的传播规律以及对涂层微观结构的影响。
3.2 爆炸后形貌
3.2.1 宏观形貌
图10为无涂层钢筋混凝土板的破坏情况。通过对比可以发现,无涂层钢筋混凝土板迎爆面(图10(a))中心区域发生严重破坏,出现爆炸凹坑,混凝土破碎并裸露出钢筋,沿凹坑边缘向四周扩散大量裂缝;观察图10(b)发现,无涂层板背爆面表面混凝土发生大面积破碎,一些破片发生脱落,钢筋发生弯曲并凸出混凝土保护层,内部混凝土损坏成块,出现大量裂缝沿破坏区向四周扩散,无涂层板发生大面积震塌破坏,完全失去承载力。爆炸冲击波可近似看成由中心点向外扩散,在板迎爆面中心处产生极大能量,使板产生较大凹坑后,穿透整个混凝土板,导致板内部破碎,背爆面震塌。
图11为有涂层钢筋混凝土板背爆面破坏情况。从图11中可以发现,T26涂层钢筋混凝土板迎爆面(图11(a))中心区域发生破坏,凹坑边缘涂层向上凸起,破坏区域近似圆形,凹坑内混凝土破碎,并且从凹坑内飞溅出爆炸碎片。在破坏凹坑边缘,涂层发生剥离,涂层出现向上凸起的现象;而有涂层板背爆面完好(图11(b)),聚脲涂层没有出现破坏,涂层没有明显凸起,结构背爆面平整。相较于无涂层钢筋混凝土板,T26涂层钢筋混凝土板的凹坑破坏较规则,破坏程度减轻,并且背爆面的聚脲涂层能完全包覆住钢筋混凝土板的破坏碎片,极大地减少爆炸碎片飞溅造成的伤害。综合来看,对于迎爆面涂层,除了瞬间高温导致聚脲软化,爆炸反射波的稀疏拉伸作用也会使聚脲材料发生破坏,导致涂层材料发生撕裂;而对于背爆面涂层,由于聚脲涂层削弱了爆炸冲击拉伸波的作用,从而保护混凝土材料不被破碎,防止爆炸碎片飞溅。
为了深入探讨T26聚脲涂层的抗爆效果,对有无涂层钢筋混凝土板的破坏区域进行测量,结果见表4。从表4中可以得到:无涂层板迎爆面凹坑最大破坏区域直径达到41.0 cm,破坏深度为11.2 cm,而有涂层板迎爆面破坏直径为24.5 cm,破坏深度为9.4 cm。相较于无涂层板,有涂层板的迎爆面破坏程度大大减弱,破坏直径减小了40.24%,破坏深度减小了16.07%。无涂层板背爆面破坏直径达到60.0 cm,破坏凸起高度为5.0 cm,而有涂层板背爆面则表现良好,没有发生破坏和凸起。从碎片情况来看,有涂层板无论是迎爆面还是背爆面,碎片均减少,并且背爆面没有产生碎片。可见,T26聚脲涂层能够减弱钢筋混凝土板的破坏程度,背爆面涂层的抗爆能力比迎爆面更好,减少碎片的能力也更优。
表 4 有无涂层板的破坏结果对比Table 4. Comparison of failure results of coated plates with uncoated platesPlate Damage diameter Failure depth or bulge height Uncoating/cm Coating/cm Error/% Uncoating/cm Coating/cm Error/% Blast side 41.0 24.5 40.24 60.0 0 100 Rear side 11.2 9.4 16.07 5.0 0 100 Debris situation More and wide
dispersion rangeLess and more
concentrated
in pitsLarger, larger
splash rangeNo debris 3.2.2 微观形貌
图12显示了爆炸后T26聚脲涂层区域Ⅰ、区域Ⅱ和区域Ⅲ的局部微观形貌(放大2000倍)。从图12(a)可以看出,断面区域表面呈鳞片状分散式分布,表面排布错乱不平整,并且伴有较多撕裂层及少量的脆性破坏。可见,聚脲涂层部分区域因爆炸产生的高温而发生熔融,同时强大的冲击波使聚脲涂层产生撕裂破坏,爆炸结束后,温度瞬间降低,导致聚脲从玻璃态转变为黏弹态,断裂的聚脲固化成鳞片状附着在断面表层。观察图12(b)发现,聚脲涂层断面不均匀,组织致密,发生因较大塑性变形导致的明显的撕裂破坏,并且伴有较多裂缝,部分区域出现熔融现象。可见,聚脲在抵抗爆炸作用时受到高温和强大的冲击波作用,导致聚脲涂层表面发生不同程度的撕裂破坏及大量裂缝。从图12(c)可以看出,涂层边缘区域表面较光滑,没有发生熔融及撕裂破坏,说明边缘涂层受爆炸荷载作用较小,爆炸产生的大部分能量均在中心处释放,但是出现较多微小裂纹,并且裂纹由涂层凹陷处向四周延伸。图12(c)中凹陷处形状较规则,可能是由于原料中的气泡在喷涂过程中被破坏,导致涂层表面形成圆孔缺陷。在爆炸形成的冲击力下,圆孔处强度较低,导致裂纹从圆孔处向四周扩展,因此,在聚脲的喷涂过程中应尽量减少气泡的产生。
4. 抗爆机理验证分析
当爆炸冲击波从钢筋混凝土板迎爆面传播至背爆面时,板的破坏以爆炸凹坑、爆炸震塌、爆炸贯穿和爆炸冲切4种破坏模式为主[18-19]。本研究中钢筋混凝土板的破坏以迎爆面爆炸凹坑及背爆面爆炸震塌为主:无涂层板的迎爆面出现爆炸凹坑破坏,背爆面则完全震塌,混凝土成块脱落;有涂层板的迎爆面出现爆炸凹坑,而背爆面没有破坏,涂层完好无损。
爆炸能够瞬间产生大量的热量和强大的冲击能量,以波的形式向周围传播,可以穿透钢筋混凝土板造成其局部或整体破坏。观察爆炸后涂层的微观形貌可以发现,在强大的爆炸冲击波下涂层的破坏以撕裂破坏为主,并伴有少量脆性破坏,可见材料在高温高热下会转变为玻璃态,表现为刚性。通过对T26聚脲的WXRD和DSC分析可知,T26聚脲的分子链是一种微相分离结构,软段和硬段之间的分离程度较大并且排列有序,导致聚脲既有足够的力学强度,又能表现出类似于橡胶的柔软性。同时,随着温度的不断升高,T26聚脲会结晶和熔融,从而吸收热量,导致材料的自由体积逐渐膨胀,分子链段运动变得活跃并逐渐舒展;此外,T26聚脲的断裂伸长率较高,对于背爆面上的涂层来说,涂层承受的温度相对较低,仍保持其黏弹性行为,能够在爆炸荷载作用下通过自身的变形起到防护作用,包覆产生的爆炸碎片,达到“以柔克刚”的效果。
由爆炸应力波的破岩过程可知,在炸药爆炸初期,强烈的高温高压作用在迎爆面涂层,使涂层产生熔融和破坏,导致混凝土板形成粉碎区,此时主要表现为强冲击波的超高压对结构施加冲击压缩破坏效应[20]。由于粉碎区消耗冲击波的部分能量,使冲击波衰减成压应力波,压应力波继续传播使背爆面涂层主要承受拉应力,产生拉伸变形。因此,对于迎爆面涂层,需要承受更大的爆炸能量,由有涂层钢筋混凝土的接触爆炸实验可以知道,迎爆面涂层会出现穿孔破坏,通过自身的大变形吸收能量,首次减弱冲击波的强度,降低传播的能量。当爆炸冲击波经混凝土传播至背爆面涂层时,能量再次衰减,背爆面涂层能够承受剩余的能量,通过自身的变形使能量完全耗散,并能完全包裹混凝土破坏产生的碎片,起到良好的防震塌效果。迎爆面涂层和背爆面涂层的两次作用使钢筋混凝土板保持良好的完整性及稳定性,对钢筋混凝土板的抗爆能力具有重要意义。
综上研究可以发现,聚脲材料的抗爆机理与其自身热性能及力学性能密切相关。材料分子链段中硬段与软段的分离程度决定其力学性能,分离程度越高,力学性能越优异,抗爆效果越明显。此外,在抵抗爆炸荷载的过程中,迎爆面涂层和背爆面涂层分别依赖于其不同的力学特点:迎爆面涂层通过更高的强度和吸能效率吸收爆炸应力波,减弱反射波,从而抵抗爆炸能量;而背爆面涂层则需要更加柔软,通过橡胶态的性质包裹爆炸碎片,极大地减弱混凝土板的震塌损伤。
5. 结 论
(1) 对T26聚脲进行了准静态力学实验,得到T26聚脲的抗拉强度为25.4 MPa,撕裂强度为75.5 N/mm,断裂伸长率为451.9%,表现出更高的力学强度。通过WXRD和DSC得到,材料分子链是微相分离结构,软段与硬段之间排列规整有序,结晶度为24.11%;软段玻璃化转变温度为−44.9 ℃,硬段玻璃化转变温度为36.5 ℃。
(2) 对10 mm厚T26聚脲涂层进行1.4 kg TNT当量的接触爆炸,结果表明:无涂层钢筋混凝土板出现爆炸凹坑及爆炸震塌破坏,凹坑的最大直径为41.0 cm,最大深度为11.2 cm,背爆面破坏区域的最大直径为60.0 cm,板结构失去整体稳定性,并产生大量的爆炸碎片。T26聚脲涂覆钢筋混凝土板出现爆炸凹坑破坏,凹坑中心区的最大破坏直径为24.5 cm,最大深度为9.4 cm,背爆面涂层没有发生破坏和凸起,T26涂层能够良好地包裹爆炸碎片,降低碎片飞溅的伤害。综合来看,对于迎爆面涂层,除了瞬间高温导致聚脲软化,还有爆炸反射波的稀疏拉伸作用也使聚脲材料发生破坏,导致涂层材料发生撕裂;对于背爆面涂层,由于聚脲涂层削弱了爆炸冲击拉伸波的作用,从而保护混凝土材料不被破碎,防止爆炸碎片飞溅。
(3) 涂层的抗爆机理与其力学性能及热性能密切相关。T26聚脲具有优异的力学性能、良好的热稳定性以及较高的结晶度,能够在温度升高时产生熔融吸热峰,导致材料内部分子链断裂及重新组合,消耗部分能量,最终使其能够承受爆炸产生的荷载与能量的多方面耦合作用。迎爆面和背爆面涂层表现出不同的防护方式,迎爆面涂层通过变形吸收能量,背爆面涂层则展现出更好的包覆能力。
-
表 1 不同试件类型失效应变取值
Table 1. Failure strain of different specimens
Specimen Stress triaxiality Failure strain Initial stress triaxiality Average stress triaxiality Compression –0.333 Torsion 0.000 0.0006 1.340 Smooth 0.333 0.562 1.273 R=18 mm 0.413 0.663 1.140 R=8 mm 0.505 0.752 1.045 R=6 mm 0.556 0.805 0.990 R=2 mm 0.893 1.085 0.791 表 2 JC失效模型参数拟合结果
Table 2. JC failure model constants
表 3 尖头战斗部穿甲8 mm厚Q345B钢金属板仿真计算结果与试验结果对比
Table 3. Comparison between experimental results and simulation results for sharp nosed warhead penetrating 8 mm thick Q345B steel targets
Material model Impact velocity/(m·s–1) Residual velocity/(m·s–1) Strength model Failure model Experiment Simulation Low reference strain rate JC-F-1 208 185 180 Low reference strain rate JC-F-2 208 185 177 High reference strain rate JC-F-1 208 185 184 High reference strain rate JC-F-2 208 185 183 Low reference strain rate JC-F-3 208 185 180 Low reference strain rate JC-F-4 208 185 174 High reference strain rate JC-F-3 208 185 182 High reference strain rate JC-F-4 208 185 178 表 4 尖头战斗部穿甲15 mm厚Q345B钢金属板仿真计算结果与试验结果对比
Table 4. Comparison between experimental results and simulation results for sharp nosed warhead penetrating 15 mm thick Q345B steel targets
Material model Impact velocity/(m·s–1) Residual velocity/(m·s–1) Strength model Failure model Experiment Simulation Low reference strain rate JC-F-1 273 216 204 Low reference strain rate JC-F-2 273 216 202 High reference strain rate JC-F-1 273 216 212 High reference strain rate JC-F-2 273 216 212 Low reference strain rate JC-F-3 273 216 194 Low reference strain rate JC-F-4 273 216 181 High reference strain rate JC-F-3 273 216 201 High reference strain rate JC-F-4 273 216 194 表 5 钝头战斗部穿甲15 mm厚Q345B钢金属板仿真计算结果与试验结果对比
Table 5. Comparison between experimental results and simulation results for blunt warhead penetrating 15 mm thick Q345B steel targets
Material model Impact velocity/(m·s–1) Residual velocity/(m·s–1) Strength model Failure model Experiment Simulation Low reference strain rate JC-F-1 273 163 75 Low reference strain rate JC-F-2 273 163 50 High reference strain rate JC-F-1 273 163 174 High reference strain rate JC-F-2 273 163 172 Low reference strain rate JC-F-3 273 163 55 Low reference strain rate JC-F-4 273 163 0 High reference strain rate JC-F-3 273 163 171 High reference strain rate JC-F-4 273 163 156 -
[1] SCHEFFER D R, ZUKA J A. Practical aspects of numerical simulation of dynamic events material interfaces [J]. International Journal of Impact Engineering, 2000, 24(8): 821–842. doi: 10.1016/S0734-743X(00)00003-8 [2] ZUKA J A. Practical aspects of numerical simulations of dynamic events: effects of meshing [J]. International Journal of Impact Engineering, 2000, 24(9): 925–945. doi: 10.1016/S0734-743X(00)00012-9 [3] JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strains rates and high temperature [C]//Proceedings of the Seventh International Symposium on Ballistics. Netherland, 1983. [4] BØRVIK T, HOPPERSTAD O S, BERSTAD T, et al. A computational model of viscoplasticity and ductile damage for impact and penetration [J]. European Journal of Mechanics Solids, 2001, 20(5): 685–712. doi: 10.1016/S0997-7538(01)01157-3 [5] GUPTA N K, IQBAL M A. Experiment and numerical studies on the behavior of thin aluminum plates subjected to impact by blunt and hemispherical-nosed projectile [J]. International Journal of Impact Engineering, 2006, 32(12): 1921–1944. doi: 10.1016/j.ijimpeng.2005.06.007 [6] BØRVIK T, HOPPERSTAD O S. Numerical simulation of plugging failure in ballistic penetration [J]. International Journal of Solids and Structures, 2001, 38(25): 6241–6264. [7] BAO Y B, WIERZBICKI T. On fracture locus in the equivalent strain and stress triaxiality space [J]. International Journal of Mechanical Sciences, 2004, 46(12): 81–98. [8] GILIOLI A, WIERZBICKI T. Predicting ballistic impact failure of aluminium 6061-T6 with the rate-independent Bao-Wierzbicki fracture model [J]. International Journal of Impact Engineering, 2015, 76(15): 207–220. [9] TENG X, WIERZBICKI T. Evaluation of six fracture models in high velocity perforation [J]. Engineering Fracture Mechanics, 2016, 73(12): 1653–1678. [10] 陈刚. 半穿甲战斗部弹体穿甲效应数值模拟和实验研究[D]. 绵阳: 中国工程物理研究院, 2006.CHEN G. Numerical and experimental investigation on penetration effects of semi-armor-piercing warhead [D]. Mianyang: China Academy of Engineering Physics, 2006. [11] 姚熊亮, 徐小刚. 船用917钢抗冲击性能试验 [J]. 中国造船, 2004, 45(4): 35–41. doi: 10.3969/j.issn.1000-4882.2004.04.006YAO X L, XU X G. The experiment research of steel 917 under impact load [J]. Ship Building of China, 2004, 45(4): 35–41. doi: 10.3969/j.issn.1000-4882.2004.04.006 [12] 李营. 液舱防爆炸破片侵彻作用机理研究 [D]. 武汉: 武汉理工大学, 2014.LI Y. Study of anti-blast fragment penetration mechanism in liquid tank [D]. Wuhan: Wuhan University of Technology, 2014. [13] 伍星星, 孟利平, 刘建湖, 等. 大质量战斗部穿甲数值仿真对材料断裂极限参数确定分析[C]//第十五届战斗部与毁伤会议集. 北京, 2017.WU X X, MENG L P, LIU J H, et al. Parameter determination analysis of material fracture limit by numerical simulation armor piercing in mass warhead [C]//The 15th Warhead Damage Conference. Beijing, 2017. [14] BØRVIK T, HOPPERSTAD O S. On the influence of stress triaxiality and strain rate on the behavior of a structural steel. Part II. numerical study [J]. European Journal of Mechanics A-Solids, 2003, 22(1): 15–32. doi: 10.1016/S0997-7538(02)00005-0 [15] 陈继恩. 基于应力三轴度的材料失效研究 [D]. 武汉: 华中科技大学, 2012.CHEN J E. Research of material failure basic on stress triaxiality [D]. Wuhan: Huazhong University of Science and Technology, 2012. [16] 孟利平. 应变率和应力三轴度对船用钢变形和断裂的影响研究 [D]. 无锡: 中国船舶科学研究中心, 2016.MENG L P. Study on the influence of strain rate and stress triaxiality on deformation and fracture of marine steel [D]. Wuxi: China Ship Science Research Center, 2016. [17] 谢凡, 张涛, 陈继恩, 等. 应力三轴度的有限元计算修正 [J]. 爆炸与冲击, 2012, 32(1): 8–14. doi: 10.3969/j.issn.1001-1455.2012.01.002XIE F, ZHANG T, CHEN J E, et al. Updating of the stress triaxiality by finite element analysis [J]. Explosion and Shock Waves, 2012, 32(1): 8–14. doi: 10.3969/j.issn.1001-1455.2012.01.002 [18] 伍星星, 张伦平, 刘建湖, 等. 不同类型战斗部穿甲不同金属板试验与理论分析 [C]//第十五届战斗部与毁伤会议集. 北京, 2017.WU X X, ZHANG L P, LIU J H, et al. Experimental and theory analysis on different nosed warhead penetrating different thickness metal [C]//The 15th Warhead Damage Conference. Beijing, 2017. -