Similar Analysis of PELE Penetrating Metal Target Fragmentation Effect
-
摘要: 为了研究横向效应增强型侵彻体(penetrator with enhanced lateral effects, PELE)侵彻金属靶板破碎效应的相似规律,选取PELE的壳体破碎长度和靶后破片散布半径作为衡量PELE破碎效应的两个物理参量,基于量纲理论对PELE破碎效应问题进行相似分析,应用AUTODYN软件开展了4组相似模型数值模拟,并进行了两组相似模型验证试验。研究结果表明:通过相似理论分析,确定了PELE破碎效应满足严格的几何相似律。在800~2000 m/s撞击速度范围内,归一化处理的壳体破碎长度和靶后破片散布半径数值模拟结果及试验结果与几何尺寸无关,仅随撞击速度的提升呈线性增长,从而证明了PELE侵彻金属靶的破碎效应满足几何相似律。
-
关键词:
- 横向效应增强型侵彻体 /
- 侵彻 /
- 破碎效应 /
- 量纲分析 /
- 相似律
Abstract: To study similar law of the fragmentation effect of penetrator with enhanced lateral effects (PELE) penetrating a metal target plate, breaking length of the PELE jacket and scattering radius of the fragments behind target are selected as two physical parameters to measure the fragmentation effect of PELE. A similar analysis on the fragmentation effect of PELE was conducted based on dimensional theory. Four groups of scaling model simulations were carried out using AUTODYN software, and two groups of scaling model verification tests were done. It is determined that the fragmentation effect of PELE satisfies geometric similar law through the similar theory analysis, and in the range of 800–2000 m/s impact velocity, the normalized breaking length of jacket and dispersion radius of fragments both increase linearly with the impact velocity in the simulation and test results, which proves that the fragmentation effect of PELE penetrating the metal target satisfies the geometric similarity law. -
夹芯结构由于具有较好的吸能缓冲性能以及较高的比强度和比刚度,被广泛应用于航空航天、交通运输、航海等领域。近年来随着对夹芯板研究的不断深入,提出了各种桁架夹芯和褶皱型夹芯等形式。其中,层级结构思想引起了众多学者的关注。Lakes[1]通过研究发现,层级结构与非层级结构相比具有更加显著的刚度和强度,同时也具有足够的稳定性及可靠性。Zhang等[2-3]通过梯形波纹板夹芯结构的静态压缩和低速冲击实验,研究了不同芯层材料、不同载荷下结构的能量吸收和变形过程,结果表明不同芯层材料和冲击能量下结构的能量吸收性能不同。Hou等[4]对不同层数的梯形夹芯板结构进行了实验研究,发现不同层数结构的能量吸收性能不同。Meza等[5]对多层级桁架结构的演变模式进行了分析,提出了几种不同结构的破坏形式和变形模式。Kooistra等[6]研究了层级波纹板结构的6种失效模式,并给出了结构在不同相对密度及角度下的破坏机理图。Velea等[7]提出了一种制备二级夹芯结构的方法,将一级波纹板结构芯层替换为多孔填充芯层,通过实验和理论分析发现,二级结构有着良好的刚度和强度。Wu等[8]将金字塔型夹芯结构的芯层替换为层级夹芯结构,理论推导出了结构的等效刚度,通过实验研究了芯层几何参数对结构变形模式的影响,并总结出结构的6种不同变形模式。Wu等[9]设计了二级金字塔形夹芯结构,进行了理论分析和数值模拟,结果表明二级芯层结构具有较好的强度和稳定性。Liu等[10]将传统格栅夹芯结构芯层板替换为蜂窝夹芯结构,通过实验研究和数值模拟研究发现,在压缩载荷下二级芯层结构有着较好的能量吸收性能和稳定的变形性能,在轻量化和能量吸收方面具有更好的工程应用前景。
多孔结构及蜂窝结构有着良好的力学性能和能量吸收性能[11]。Jing等[12]通过实验研究了3种不同多孔芯层夹芯板结构在脉冲载荷下夹芯板的变形模式,结果表明不同芯层夹芯结构的变形模式不同,芯层和面板尺寸不同的夹芯结构的力学性能不同。Chen等[13-14]根据层级蜂窝结构的设计理论设计了不同芯层形式的层级蜂窝结构,并由理论分析得到结构变形模式,利用有限元对结构变形模式进行了验证,与传统蜂窝结构相比,层级蜂窝结构具有更好的能量吸收性能。Yin[15]、Qiao[16]等利用有限元软件对层级蜂窝结构进行了数值模拟,研究结果表明层级蜂窝结构的力学性能远优于传统蜂窝。
已有的研究主要集中于传统波纹板夹芯结构和层级蜂窝结构,关于多层级夹芯结构的力学性能与能量吸收的研究较少。本研究基于波纹板夹芯结构,设计了3种二级波纹板夹芯结构,通过理论分析和数值模拟相结合的方法研究一级和二级夹芯结构的变形模式与能量吸收性能,理论预测结构的临界失效载荷;建立不同芯层层数的二级波纹板夹芯结构的有限元模型,研究不同芯层厚度对多层级波纹板夹芯结构变形模式和能量吸收性能的影响,并与传统的一级夹芯波纹板结构进行对比,分析其比吸能与结构效率。
1. 层级夹芯结构的芯层设计和理论分析
1.1 芯层结构设计
根据层级蜂窝结构的设计理论,将一级梯形波纹板夹芯结构的芯层替换为三角形多孔夹芯结构,共设计了3种芯层结构,如图1所示。图1(a)为一级波纹板夹芯结构的芯层,斜板长度为la,顶板长度为l,斜板角度为
ω ,芯层板厚度为ta;二级结构的芯层如图1(b)、图1(c)和图1(d)所示,其中图1(b)为二级单层结构,图1(c)和图1(d)分别为二级双层和三层结构,二级结构芯层由大支撑面板和小支撑构成,水平方向小支撑为横向小支撑。二级结构中不同层数结构的尺寸参数相同,斜向大支撑面板长度为la,厚度为ta,小支撑长度为lb,厚度为tb,小支撑与外板角度均为θ ,外板与水平方向夹角为ω 。结构垂直于纸面方向的厚度为30 mm。1.2 理论分析
小支撑厚度与大支撑面板厚度相同的结构模型如图2所示,在外部载荷Fx和Fy作用下考虑横向小支撑的变形,基于小变形假设,由受力分析可以得到结构变形和力之间的关系,通过几何分析可以得到合力与变形之间的关系
{Fx2=Fncosω+FtsinωFy2=Fnsinω−Ftcosω (1) {δx=δncosω+δtsinωδy=δnsinω−δtcosω (2) 式中:Fx和Fy分别为沿x轴和y轴的外部载荷,
δx 和δy 为沿x轴和y轴的变形,Fn和Ft分别为夹芯夹层在外部载荷下的轴向压缩分量和剪切分量,δn 和δt 分别为夹芯夹层的轴向变形和剪切变形。当只有Fx作用时,
δy =0,δx ≠0,此时有Fn=EAδnl=Fx2cosω (3) δn=FnlEA (4) 式中:EA为夹芯夹层截面的抗拉刚度。由结构几何关系可知
δx=δncosω (5) 当只有Fy作用时,
δy ≠0,δx =0,此时有Fn=EAδnl=Fy2sinω (6) 同样有
δy=δnsinω (7) 当结构芯层为二级结构时,如图2(d)所示,(1)式~(7)式仍然成立。根据文献[17]给出的不同尺寸比下结构的失效机制,当
lb/la>√εy/πsinθ (其中εy 为材料屈服应变)时层级结构出现大支撑面板塑性屈服的失效模式。本研究中结构模型的尺寸满足上述条件,因此失效模式为大支撑面板发生塑性屈服,即大支撑面板达到材料屈服应力时结构失效。由(1)式和(2)式可得单胞结构大支撑面板在y方向上发生塑性屈服时结构的合力Fy′=(2n+2)σybtasinω (8) 式中:n为二级结构层数,b为结构垂直于纸面方向的厚度,σy为屈服应力。
根据结构变形模式,考虑横向小支撑的变形,如图3所示,Fb为横向小支撑在y方向的力,
δb 为δn 在y方向的分量,由此可知δb=Fbl3b3D (9) 弯曲刚度
D=Et3b12(1−ν2) (10) 式中:E为弹性模量,
ν 为泊松比。由此可以得到Fb=2nσyt3bsinω4(1−ν2)l2b (11) 式中:tb和lb分别为二级结构芯层小支撑厚度和长度。二级结构单胞在y方向的合力为
Fy=Fy′+Fb (12) 2. 有限元模型
有限元模型主要包括3部分:刚性压板、刚性支撑底板和波纹板夹芯结构。波纹板夹芯结构由芯层和上、下面板构成,上、下面板设置为刚体,厚度为2 mm。图4为二级双层结构有限元模型。波纹板夹芯结构和下部底板采用壳单元,刚性压板采用实体单元。刚性压板速度V=1 m/s,下部刚性底板固定,中部波纹板夹芯结构芯层与上、下两面板之间采用绑定约束。结构材料为铝合金6061-T6,材料参数如表1所示,其中
ρ 为密度。图5为材料应力-应变曲线[6]。结构单元类型采用四节点线性缩减积分壳单元(S4R),单元特征长度L=1 mm。芯层和上、下面板为自接触,刚性压板与上面板、底板与下面板为面面接触,静摩擦系数和动摩擦系数分别为0.20和0.15[18]。表 1 6061-T6铝合金材料性能Table 1. Material properties of aluminum 6061-T6Material ρ/(g·cm–3) σy/MPa E/GPa ν Al 6061-T6 2.700 251 69 0.33 2.1 准静态分析
利用LS-DYNA进行准静态分析时[19]:首先在模拟过程中结构动能要相对于内能较小,其次结构的力-位移曲线相对于速度的变化是稳定的。如图6(a)所示,结构在压板速度为0.5 m/s和1.0 m/s时动能相对于内能较小,结构整体内能受速度变化的影响较小;图6(b)为结构在压板速度为0.5 m/s和1.0 m/s时的力-位移曲线,可以看出两条曲线的重合度较高,速度变化对结构力-位移曲线的影响较小。通过对结构进行网格敏感性验证,并考虑到计算时间等因素,计算时取压板的速度为1.0 m/s,网格单元尺寸为1 mm。
2.2 能量吸收指标
为了评估结构的能量吸收性能[20],选取总能量吸收(EA)、平均载荷(Pm)、比吸能(
ξEA )、载荷效率(AE)4个评价指标。总能量吸收EA为结构在变形过程中吸收的能量
EA=∫δ0F(x)dx (13) 式中:
δ 为结构的变形位移(取δ =70 mm),F(x)为压缩过程中的瞬时载荷。平均载荷Pm定义为单位变形长度的能量吸收
Pm=EA/δ (14) 比吸能
ξEA 定义为结构单位质量的能量吸收ξEA=EA/m (15) 式中:m为结构总质量。比吸能越大,结构的能量吸收性能越好。
载荷效率定义为平均载荷与峰值载荷的比值
AE=Pm/Pk (16) 式中:Pk为峰值载荷。较好的吸能结构应该具有较大的载荷效率,并且力-位移曲线波动较小。
3. 模拟结果与分析
3.1 芯层板厚度的影响
建立一级和二级结构模型,其中一级结构芯层壁厚为2 mm,二级结构模型中ta=tb,即二级结构内部小支撑厚度与外部大支撑面板厚度相同。层级结构芯层板厚度分别为0.565、0.75、1.00、1.25 mm,结构几何参数如表2所示。图7给出了一级和二级层级结构在不同芯层板厚度下的力-位移曲线。由图7可知,二级结构平台力明显高于一级结构。一级结构在靠近芯层上部产生塑性铰,随着变形的增大未产生新的塑性铰,力-位移曲线变化平稳。二级结构随着变形增大发生折叠塑性变形,力-位移曲线会随之发生波动;在芯层厚度为1.25 mm时二级单层结构由于变形模式发生变化,力-位移曲线发生波动而随之上升;随着二级结构层数的增加,结构平台力增大;随着二级结构厚度的增大,其平台力增大。
表 2 模型尺寸Table 2. Size parameters of the modelType l/mm la/mm lb/mm ta/mm tb/mm θ/(°) ω/(°) First order 50 100 2 60 Second order-1 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 Second order-2 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 Second order-3 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 图8所示为不同芯层厚度下结构失效载荷理论值与模拟值的对比。对于数值模拟结果,当结构的承载力达到峰值时认为结构失效,即取力-位移曲线上的峰值载荷为结构失效载荷,可以看出理论值与模拟值吻合较好。当结构厚度为0.565 mm时,由于其厚度较小,在压缩载荷下的变形模式受屈曲影响较大,故理论值较模拟值偏大;当结构厚度为1.25 mm时,小支撑端部变形受边界条件的影响较大,理论值相对于模拟值偏小。
3.2 变形模式
图9所示为芯层板厚度t=2 mm时一级结构在不同应变下的变形模式,应变为结构在受载荷方向上的变形位移与结构上、下两面板之间原始距离之比。由图9可知:一级结构首先发生屈曲,在靠近芯层斜支撑上部产生塑性铰,随着变形的增大,结构发生失稳,从而产生非对称变形。图10所示为二级结构在应变为0.3时不同芯层厚度下的变形模式,可以看出:当厚度较小时,结构的变形发生在整体斜支撑与整体水平支撑的交界处,整体斜支撑内部的小支撑发生逐层折叠的塑性变形;当芯层厚度为1.25 mm时,结构斜支撑整体产生塑性铰而发生非对称变形;厚度变化对二级双层和三层结构的变形模式影响较小。在芯层厚度较小时,二级结构主要为小支撑整体逐层发生塑性变形,从而有更多能量转换为非弹性能,提高了结构整体的能量吸收能力。
3.3 能量吸收性能分析
二级结构变形模式主要为小支撑逐层折叠的塑性变形,相对于一级结构,其在变形过程中小支撑不断发生塑性变形产生塑性铰,从而有更多的能量转换为非弹性能,提高了结构的能量吸收性能。图11为不同芯层厚度结构的比吸能柱状图。由图11可知:层级波纹板夹芯结构的比吸能高于传统波纹板夹芯结构;随着结构芯层厚度增大,结构的比吸能不断提高;由于二级单层结构变形模式的变化,芯层厚度较小时二级单层结构的比吸能高于二级双层及三层结构,厚度较大时二级单层结构的比吸能低于二级双层及三层结构;二级双层结构的比吸能略高于二级三层结构。
图12为结构载荷效率随芯层厚度的变化关系。随着芯层厚度的增大,结构的载荷效率提高,二级单层结构由于变形模式的变化,载荷效率发生波动。由于二级单层结构在芯层厚度为1.25 mm时结构变形模式为斜支撑整体发生失稳,相对于小支撑发生折叠变形时其能量吸收性能降低,从力-位移曲线可以看出,结构平台力减小导致结构的载荷效率下降。二级双层和三层结构的载荷效率不断增加,二级三层结构的载荷效率高于二级双层结构。
4. 结 论
设计了不同层数的层级波纹板夹芯结构,利用数值模拟方法研究了在压缩载荷下层级夹芯结构的变形规律与能量吸收性能,理论推导了结构临界失效载荷公式;分析了结构参数对其变形模式和能量吸收性能的影响,并与一级结构进行了对比分析,得到如下研究结果。
(1)理论分析得到的结构失效载荷与数值模拟结果吻合较好。
(2)在准静态压缩载荷作用下,一级结构首先发生屈曲,随着变形增大,结构发生失稳而产生非对称变形;二级单层结构在厚度为1.25 mm时由于芯层整体产生塑性铰而失稳;二级双层和二级三层结构在整体斜支撑与整体水平支撑的交界处,整体斜支撑内部的小支撑发生逐层折叠的塑性变形。
(3)二级波纹板夹芯结构的比吸能显著大于一级结构;二级结构芯层小支撑发生逐层折叠的塑性变形时,结构能量吸收性能较好;随着芯层厚度的增大,二级结构的比吸能和载荷效率增加,芯层厚度较小时二级单层结构的比吸能高于二级双层和三层结构,二级双层结构的比吸能略大于二级三层结构。
-
表 1 PELE和靶板材料模型参数
Table 1. Material model parameters of PELE and metal target
Component ρ/(g·cm−3) Grüneisen coefficient C1/(km·s−1) S1 Shear modulus/GPa Yield
stress/GPaPlastic strain Jacket (D-180 K) 18.000 0 4.03 1.237 139.0 1.50 Filling (A-G3) 2.650 1.97 5.24 1.400 27.5 0.30 Steel target
(XC 48)7.823 0 4.57 1.490 77.0 0.80 1.2 Aluminum target (A-U4G) 2.800 2.0 5.20 1.360 26.7 0.40 1.2 Component Principal tensile
failure stress/GPaPrincipal tensile
failure strainCrack softening Fracture energy/(J·m−2) Stochastic failure Stochastic variance Erosion strain Jacket (D-180 K) 2.8 0.035 Yes 45 Yes 36.65 0.6 Filling (A-G3) 0.5 No No 0.8 Steel target
(XC 48)No No Aluminum target (A-U4G) No No 表 2 不同比例模型的几何参数
Table 2. Geometric parameters of different scale model
λ d1/mm d2/mm L/mm l/mm D/mm H/mm s/mm v/(m·s−1) 1 5 3 25 22.5 10 1.5 0.125 800–2000 2 10 6 50 45.0 20 3.0 0.250 3 15 9 75 67.5 30 4.5 0.500 4 20 12 100 90.0 40 6.0 0.750 表 3 PELE弹丸和靶板的几何尺寸
Table 3. Geometric dimension of PELE projectile and target plate
λ d1/mm d2/mm L/mm l/mm H/mm v/(m·s−1) 1 10 7 25 22.5 3 800–1600 2 20 14 50 90.0 6 -
[1] 程春, 杜忠华, 陈曦, 等. 壳体刻槽PELE冲击破碎数值模拟 [J]. 弹道学报, 2019, 31(3): 92–96.CHENG C, DU Z H, CHEN X, et al. Numerical simulation for impact fragmentation of jacket grooved PELE [J]. Journal of Ballistics, 2019, 31(3): 92–96. [2] 杜忠华, 宋丽丽. 横向效应增强型侵彻体撞击金属薄板理论模型 [J]. 南京理工大学, 2011, 35(6): 822–826.DU Z H, SONG L L. Theoretical model of penetrator with enhanced lateral effect impacting thin metal target [J]. Journal of Nanjing University of Science and Technology, 2011, 35(6): 822–826. [3] 朱建生, 赵国志, 杜忠华, 等. 小口径PELE作用薄靶板影响因素的实验研究 [J]. 实验力学, 2007, 22(5): 505–510. doi: 10.3969/j.issn.1001-4888.2007.05.010ZHU J S, ZHAO G Z, DU Z H, et al. Experimental study of the influence factors on small caliber PELE impacting thin target [J]. Journal of Experimental Mechanics, 2007, 22(5): 505–510. doi: 10.3969/j.issn.1001-4888.2007.05.010 [4] 朱建生, 赵国志, 杜忠华, 等. PELE垂直侵彻薄靶的机理分析 [J]. 爆炸与冲击, 2009, 29(3): 281–288. doi: 10.3321/j.issn:1001-1455.2009.03.010ZHU J S, ZHAO G Z, DU Z H, et al. Mechanism of PELE projectiles perpendicularly impacting on thin target plates [J]. Explosion and Shock Waves, 2009, 29(3): 281–288. doi: 10.3321/j.issn:1001-1455.2009.03.010 [5] 朱建生, 赵国志, 杜忠华, 等. 着靶速度对PELE横向效应的影响 [J]. 力学与实践, 2007, 29(5): 12–16. doi: 10.3969/j.issn.1000-0879.2007.05.003ZHU J S, ZHAO G Z, DU Z H, et al. The influence of impact velocity on the penetrator with enhanced lateral effects (PELE) [J]. Mechanics in Engineering, 2007, 29(5): 12–16. doi: 10.3969/j.issn.1000-0879.2007.05.003 [6] CHENG C, XU L Z, PANG Z J, et al. Fragmentation characteristics and damage capacity of radial layered PELE [J]. Latin American Journal of Solids and Structures, 2020, 17(8): e327. [7] 朱建生, 赵国志, 杜忠华, 等. 靶板厚度对横向效应增强型侵彻体作用效果的影响 [J]. 南京理工大学学报(自然科学版), 2009, 33(4): 474–479.ZHU J S, ZHAO G Z, DU Z H, et al. Influence of target thickness on lateral effect of PELE [J]. Journal of Nanjing University of Science & Technology (Natural Science Edition), 2009, 33(4): 474–479. [8] 蒋建伟, 张谋, 门建兵, 等. 不同内核材料PELE弹丸对多层靶穿甲实验研究 [J]. 北京理工大学学报, 2010, 30(9): 1009–1012. doi: 10.15918/j.tbit1001-0645.2010.09.009JIANG J W, ZHANG M, MEN J B, et al. Experimental study on multi-layered target penetration of PELE with different cores [J]. Transactions of Beijing Institute of Technology, 2010, 30(9): 1009–1012. doi: 10.15918/j.tbit1001-0645.2010.09.009 [9] PAULUS G, SCHIRM V. Impact behaviour of PELE projectiles perforating thin target plates [J]. International Journal of Impact Engineering, 2006, 33(1): 566–579. [10] VERREAULT J. Analytical and numerical description of the PELE fragmentation upon impact with thin target plates [J]. International Journal of Impact Engineering, 2015, 76: 196–206. doi: 10.1016/j.ijimpeng.2014.09.012 [11] VERREAULT J, HINSBERG N V, ABADJIEVA E. PELE fragmentation dynamics [C]//27th International Symposium on Ballistics. Freiburg, Germany: International Ballistics Society, 2013. [12] FAN Z J, RAN X W, TANG W H, et al. The model to calculate the radial velocities of fragments after PELE penetrator perforating a thin plate [J]. International Journal of Impact Engineering, 2016, 95: 12–16. doi: 10.1016/j.ijimpeng.2016.04.011 [13] 佚名. 聚能射流侵彻过程模型律 [J]. 力学学报, 1974, 1: 3–12.YI M. Model law of penetration process of shaped jet [J]. Chinese Journal of Theoretical and Applied Mechanics, 1974, 1: 3–12. [14] 陶为俊, 黄风雷, 浣石. 聚能射流侵彻过程的相似律及其在数值仿真中的应用 [J]. 广州大学学报 (自然科学版), 2008, 7(1): 29–33.TAO W J, HUANG F L, HUAN S. Analogue law of the penetration process of shaped charge jet and application on numerical simulation [J]. Journal of Guangzhou University (Natural Science Edition), 2008, 7(1): 29–33. [15] 程怡豪, 李干, 岳松林, 等. 混凝土超高速侵彻效应的相似规律 [J]. 力学与实践, 2019, 41(5): 543–549.CHENG Y H, LI G, YUE S L, et al. Similarity laws of hypervelocity penetration effects of concrete [J]. Mechanics in Engineering, 2019, 41(5): 543–549. [16] 彭永, 卢芳云, 方秦, 等. 弹体侵彻混凝土靶体的尺寸效应分析 [J]. 爆炸与冲击, 2019, 39(11): 1–10. doi: 10.11883/bzycj-2018-0402PENG Y, LU F Y, FANG Q, et al. Analyses of the size effect for projectile penetrations into concrete targets [J]. Explosion and Shock Waves, 2019, 39(11): 1–10. doi: 10.11883/bzycj-2018-0402 [17] CHENG C, DU Z H, CHEN X, et al. Damage of multi-layer spaced metallic target plates impacted by radial layered PELE [J]. Defence Technology, 2020, 16(1): 201–207. doi: 10.1016/j.dt.2019.05.002 [18] 高光发. 量纲分析理论与应用 [M]. 北京: 科学出版社, 2021.GAO G F. The theory and application of dimensional analysis [M]. Beijing: Science Press, 2021. [19] DING L L, ZHOU J Y, TANG W H, et al. Damage characteristics of PELE projectile with gradient density inner core material [J]. Materials, 2018, 11(12): 1–19. [20] DING L L, ZHOU J Y, TANG W H, et al. Research on the crushing process of PELE casing material based on the crack-softening algorithm and stochastic failure algorithm [J]. Materials, 2018, 11(9): 1–19. -