Dynamic fragmentation of oxygen-free high-conducting copper under Mach stem loading
-
摘要: 为了深入了解金属材料在复杂加载下的动态破碎行为,在有限元模拟的基础上,设计了两类马赫杆加载实验,用于研究无氧铜在复杂加载下的动态破碎行为。实验中,采用火炮加载马赫透镜和激光粒子速度干涉仪测量自由面速度,实现了峰值压力分别为95.75和32.38 GPa的动态加载。结果表明,实验中成功实现了稳定的马赫杆加载,并且观察到马赫杆加载下无氧铜的2种不同近表面破碎模式,即高压下产生微层裂、低压下产生三角波层裂,且层裂区呈凸形分布。Abstract: To in-depth understand the dynamic fracture behaviors of metal materials under complex loading, based on the finite element simulation, two types of Mach stem loading experiments were designed and carried out to investigate the dynamic fragmentation of oxygen-free high-conducting copper (OFHC Cu) under complex loading. In the experiments, a powder gun was used to impact the Mach lens, and a laser particle-velocity interferometer was applied to measure the free surface velocity. And dynamic loadings with the peak pressures of 95.75 and 32.38 GPa, respectively, were achieved. Stable Mach stem loading was successfully generated, and the Mach stem-related features were consistent with the simulated ones. At the same time, two different near-surface fracture behaviors in the OFHC Cu were observed, namely the micro-spallation under high pressure and the triangular-wave spallation under low pressure, with the cracked area distributed in a convex shape. These findings have a certain value for further understanding the dynamic fracture behaviors of metal materials and can provide new experimental methods for understanding material failure under various complex loading conditions.
-
Key words:
- Mach stem /
- micro-spallation /
- OFHC Cu /
- dynamic fragmentation
-
为了攻击机库、高级指挥要地等各种不同形式的综合防护目标, “智能型”的钻地弹药相继出现[1]。弹体高速侵彻硬目标的加速度-时间历程曲线是一个非常重要的测试量, 对于考核钻地弹的结构强度, 检验弹体设计合理性, 考核战斗部装药的安定性, 研究弹载电子仪器的存活性、引信的工作可靠性, 评价防护工事的抗侵彻效果至关重要。在现有的侵彻加速度测试技术中, 弹载存储测试技术较先进[2-4], 它能记录膛内、飞行穿靶几个过程的实时加速度。国外对于硬目标侵彻过程动态参数测试起步较早, 瑞士的武器系统和弹药实验中心设计了高g值弹道飞行数据记录器, 质量较轻(3.58 kg), 安装了3个6×104 g量程的Endevco_7270A压阻加速度计, 成功测得最大加速度(9×104 g)。国内中北大学电子测试技术国防科技重点实验室通过选购B & K公司的8309或兵器204所988加速度计, 并采用泡沫铝缓冲保护, 电路模块质量控制在150g以内, 测试范围最高可达2×10 5 g。经过多年的研究, 高g值侵彻加速度现场测试技术已经基本解决了低速和中速条件下(800 m/s以下)“测不着”的问题, 正朝着提高高速侵彻条件下数据捕获率和进一步解决“测不准”问题的方向发展。弹载测试装置的工作环境恶劣, 侵彻原型实验需要耗费大量资源, 进行必要的数值模拟分析是对侵彻实验的必要补充。对于加速度信号测试的影响分析有利于提高数据捕获率, 能够进一步优化弹载测试装置的设计。
1. 测试原理及实验
加速度测试系统原理图如图 1所示, 由4部分组成:加速度传感器、存储记录部分、数据接口、数据处理软件。测试装置如图 2所示。
实验时, 将加速度过载测试仪与配重块放入子弹中, 用压筒的一端压紧测试仪下端盖的上部平面, 另一端加配合垫片, 并用尾翼组合件拧紧到位, 保证加速度过载测试仪与弹体内腔前部平面完全贴和, 尽可能保证刚性连接。
弹丸侵彻混凝土靶板后靶板正面破坏情况如图 3所示, 可以看到靶体表面成坑、靶体表面的径向和环向裂纹破坏现象, 证明弹体在撞击靶体时具有良好的姿态。图 4是弹体侵彻混凝土靶过程的加速度曲线图。图 5为以3 k Hz截止频率滤波后的加速度曲线图, 峰值加速度为4.133×104 g。对侵彻加速度曲线二次积分得到侵彻深度曲线(见图 6), 积分位移为0.562 m, 与实测侵深0.600 m相差6.3%。
2. 弹体侵彻混凝土靶过程的数值模拟
在适量的实验基础上基于LS-DYNA进行数值模拟, 是对侵彻实验的重要补充[5-6]。
弹丸壳体是按照实验弹尺寸建立的, 内部结构作了较大的简化, 将测试仪和套筒处理成一等重的质量块装在结构前部, 结构后部则同样用一等重的质量块来代替套筒和弹尾翼。由于侵彻计算中主要关心壳体和测试仪的刚体加速度以及侵彻最大行程, 不关心各结构的受力情况, 因此以上简化具有一定的合理性。为了模拟半无限靶, 将实验用的混凝土靶板半径取足够大, 除了对称面和自由面外, 模型的边界条件是将靶板的周身加以非反射边界条件(应力波通过此边界传播而不产生反射效应)以模拟实际目标的边界效应。测试仪视为各向同性弹塑性模型, 本构模型(强度模型)采用Von Mises准则。弹内尾部的套筒可等同为一个弹性体, 材料模型采用弹性模型即MAT_ELASTIC。混凝土结构采用JHC模型, 该模型能够很好地描述高应变率条件下混凝土的响应问题。考虑到弹丸垂直侵彻靶板为对称问题, 在进行有限元分析时选取1/2实际模型建立计算模型。有限元模型网格如图 7所示。
图 8给出了弹丸侵彻靶板的数值模拟结果。从计算结果来看, 混凝土靶表面成坑漏斗形状、侵彻过程转弯等同实验现象较吻合。图 9为弹体轴向刚体加速度随时间的变化曲线, 与图 5相比, 幅值与持续时间基本相同, 与实验结果较接近, 表明了数值模型建立及参数选择的正确性。
由于应力波的传播, 弹体的不同部位在侵彻过程中的加速度是不一样的, 应力波的波前是加速度的峰值之处[7], 但波前过后峰值很快下降, 所以峰值处的脉宽一般很窄。图 10为弹壳体A、测试仪B及套筒C的侵彻加速度。通过对加速度信号的滤波, 将高频部分滤除, 最后只剩下各部分的刚体加速度, 而此刚体加速度各部分相同, 如图 11所示。
在弹体强度设计、战斗部装药安定分析等时, 弹体的整体加速度, 即刚体加速度是大家非常关心的[8], 加速度测试仪在弹体内的安装部位和安装方式不会对侵彻过程中的刚体加速度测量造成影响, 这可以通过采用适当的滤波截止频率进行滤波而得到所需要的加速度信号。
弹体内腔的前端面加工较粗糙, 为了使加速度计与弹体之间绝缘, 在加速度计与安装座间通常加绝缘纸垫, 当应力波从弹头传播到弹体与测试装置的接触面以及加速度计安装座时, 大部分应力波被反射, 所以加速度计基本不受应力波的影响。
3. 测试装置基础运动对实测加速度信号的影响分析
图 12的加速度计安装方式可以简化为图 13所示的模型, 图中c为阻尼器的阻尼系数, m为运动质量块的质量, k为弹性元件的刚度系数, x(t)为质量块所受激励, 它可被视为一个基础激励单自由度二阶系统。设加速度计的激励信号为基础运动:u=u0 sinωt, 基础加速度:af(t)=-u0ω2sinωt, 其相对基础的运动方程为
m¨x+c˙x+kx=−maf(t) (1) 整理得
¨x+2ζωn˙x+ω2nx=u0ω2sinωt 式中:ωn为无阻尼固有角频率, ζ为阻尼比。若忽略阻尼, 则式(2)为
¨x+ω2nx=−af(t) (3) 对杜哈梅积分逐次应用分部积分, 可把相对位移表示成关于 1ω的幂级数, 因加速度传感器的自振频率很高, 取其前3项为
x(t)=−1ω2n[af(t)−af(0)cosωnt−1ωn˙af(0)sinωnt] (4) 因为 −maf(t)k=−af(t)k/m=−af(t)ω2n, 上式右边第1项为加速度计对惯性力-maf(t)静力响应, 其余各项为加速度计对惯性力-maf(t)的动力响应, 它是以ω为频率的高频信号。式(4)中的后2项在加速度实测信号中反映出弹体和安装结构的高频振动的影响。
4. 测试装置安装结构刚度对实测加速度信号的影响
下面结合图 14中的测试装置安装结构, 分析铝压筒和加速度计安装座的刚度对实测加速度信号的影响。
(1) 当弹体在炮膛内加速向前运动时, 测试装置因惯性向后运动压缩铝筒, 铝压筒的刚度将影响膛内加速度信号, 甚至导致测试失败。图 15中曲线出现剧烈振荡。计算机模拟可知铝压筒的固有频率为1.6 k Hz, 而图 16中在1.5 k Hz左右出现一个尖峰, 即膛内加速度信号中包含铝压筒刚度的影响。按1.5 k Hz滤波, 可得到弹体在膛内运动时的刚体加速度, 如图 17所示。
在图 16的频谱图中, 23 k Hz附近出现一个谐振峰, 表明在该频率附近结构发生剧烈振动。采用ANSYS对弹体进行模态分析, 取整个结构的四分之一进行分析。图 18显示弹体第23阶模态(除零频外)对应的振型, 主要为测试装置上、下端盖的弯曲振动, 其振动方向与加速度计的敏感方向相同, 并且该频率(22.0 k Hz)接近加速度计的谐振频率(25 k Hz)。可以判断, 正是第23阶模态(除零频外)对应的振动使加速度计输出过大, 可能导致测试在膛内阶段就失败。
若将铝筒换为钢筒, 并不能明显提高测试装置上、下端盖的弯曲振动频率, 其值为23.37 k Hz。
ωn=nπl√EAmn=0,1,2,⋯ (5) 由自由杆纵向振动频率公式(5)可知, l(长度)一定, E(材料弹性模量)、m (质量)均增大, 频率增大不明显。
为了增加铝压筒刚度, 使其上段为实心, 质量增加0.74 kg, 整体装配后弹体重心后移2.5 mm, 对弹的飞行稳定性影响很小。对应测试装置上、下端盖的弯曲振动的频率(第13阶)为28.46 k Hz。另外, 还可将铝筒和尾翼间的调整垫片改为钢片, 以增加预压力, 提高该振动频率。采取提高铝压筒刚度的措施后, 且膛内加速度曲线振荡比以前明显降低。这就从一个方面提高了数据的质量。
(2) 当弹体侵彻靶板时, 测试装置相对弹体有向前运动的趋势, 加速度计安装座压紧弹体, 因为安装座本身的刚度大, 而加速度计是通过钢螺栓安装在基座上, 所以此时整体安装刚度大, 可以认为加速度的信号是对弹体内腔前端安装面的响应。
5. 加速度计安装方式对实测信号影响的分析
由实验可知加速度测试仪在弹体内的安装部位和安装方式不会对侵彻过程中的刚体加速度测量造成影响, 但加速度计的输出信号是其对安装点(面)加速度激励的响应。所以不同的安装方式对加速度计的输出信号有很大影响。下面对目前已采用或可能的安装方式进行讨论。
(1) 将加速度计通过钢螺栓直接安装在弹体上。这种安装刚度最大, 频响特性最好, 但无论安装在弹体内空腔的前端还是尾翼盖上, 传感器数据线与记录电路的连接都难以实现, 所以目前国内外都采用类似于图 14的结构。
(2) 将电路模块与加速度计安装成一个整体, 或将加速度计直接灌封在电路模块中, 然后再对这一整体进行缓冲保护, 目的是减短连接到电路模块上的加速度计数据线的长度, 提高数据线的抗高冲击能力。这种安装方式可简化为图 19所示的力学模型。
在侵彻过程中, 整体相对弹体以相对加速度ar向前运动, 加速度计感知的绝对加速度为a0=a+ar, 实际上并不是弹体的加速度。
若从冲击响应的角度考虑, 这个问题可简化为一个基础激励的二阶系统, 如图 20所示, 整体的运动方程为
m¨y+c˙y+ky=−m¨x (6) 令, a=¨x/¨xmax,τ=t/T,δ=−ky/(m¨xmax),R=2π√m/k/T。其中 ¨xmax为弹体加速度最大值,
为弹体加速度瞬时值, T为弹体加速度脉冲持续时间。
对上式进行量纲一化, 得到
(R2π)2d2δdτ2+Rξπdδdτ+δ=a (7) 式(7)为相对于量纲一输入a的量纲一响应δ, 从式(7)可知, 当R很小时, 方程的前2项可以忽略, 于是δ=a; 当R变大时, 前2项发生影响, 第1项引起在a值上下的振荡, 第2项引起时间滞后。这将使弹体加速度波形发生变化, 由该加速度信号积分得到的着靶速度、侵彻深度与实测值之间的误差也随之加大。
所以, 采用类似于图 14的结构将加速度计刚性安装, 目前来看是比较理想的。
6. 结论
(1) 洛阳水利工程技术研究所进行了过载特性的理论计算[9], 计算结果与实验测量数据符合较好。数值计算结果与实验结果较吻合, 表明此数值模型建立及参数选择的正确性。测试数据为实测信号的影响分析提供了事实依据。
(2) 由于应力波的传播, 弹体的不同部位在侵彻过程中的加速度是不一样的, 应力波的波前是加速度的峰值之处, 但波前过后峰值很快下降, 所以峰值处的脉宽一般很窄。在进行加速度测试时, 可根据不同的目的(如各部件的抗冲击能力等), 将加速度计安装在需要测试的位置。通过在加速度计安装座下加绝缘纸, 既可使加速度计与弹体绝缘, 又可大幅减轻应力波影响。
(3) 加速度测试仪在弹体内的安装部位和安装方式不会对侵彻过程中的刚体加速度测量造成影响, 可以通过采用适当的滤波截止频率进行滤波而得到所需要的加速度信号。加速度计的输出信号是其对安装点(面)加速度激励的响应, 不同的安装方式, 对加速度计的输出信号影响很大。在具体安装时, 加速度计应采用刚性连接, 同时应提高测试装置的安装刚度。
-
表 3 设计参数
Table 3. Parameters of experimental design
实验编号 材料 飞片速度/(km·s−1) 厚度/mm 直径/mm 飞片 外圆 内圆 飞片 样靶 内圆 外圆 Mach-1 304不锈钢 LY12 铝 无氧铜 1.40 3.0 16.0 4.8 38.0 Mach-2 304不锈钢 TC4钛合金 无氧铜 0.50 12.0 26.0 14.0 45.0 -
[1] 陈华燕, 曾祥国, 朱文吉, 等. 爆炸荷载作用下桥梁动态响应及其损毁过程的数值模拟 [J]. 四川大学学报(工程科学版), 2011, 43(6): 15-19, 97. DOI: 10.15961/j.jsuese.2011.06.001.CHEN H Y, ZENG X G, ZHU W J, et al. Numerical simulation of dynamic response and damage process for bridge under blast loading [J]. Journal of Sichuan University (Engineering Science Edition), 2011, 43(6): 15–19, 97. DOI: 10.15961/j.jsuese.2011.06.001. [2] 朱建士, 胡晓棉, 王裴, 等. 爆炸与冲击动力学若干问题研究进展 [J]. 力学进展, 2010, 40(4): 400–423. DOI: 10.6052/1000-0992-2010-4-J2009-144.ZHU J S, HU X M, WANG P, et al. A review on research progress in explosion mechanics and impact dynamics [J]. Advances in Mechanics, 2010, 40(4): 400–423. DOI: 10.6052/1000-0992-2010-4-J2009-144. [3] 陈大伟, 王裴, 孙海权, 等. 爆轰波对碰驱动平面锡飞层的动力学及动载行为特征研究 [J]. 物理学报, 2016, 65(2): 024701. DOI: 10.7498/aps.65.024701.CHEN D W, WANG P, SUN H Q, et al. Loading characteristics and dynamic behaviors of the plane tin flying layer driven by detonation collision [J]. Acta Physica Sinica, 2016, 65(2): 024701. DOI: 10.7498/aps.65.024701. [4] 张凤国, 刘军, 何安民, 等. 强冲击加载下延性金属卸载熔化损伤/破碎问题的物理建模及其应用 [J]. 物理学报, 2022, 71(24): 244601. DOI: 10.7498/aps.71.20221340.ZHANG F G, LIU J, HE A M, et al. Modelling of spall damage evolution and fragment distribution for melted metals under shock release [J]. Acta Physica Sinica, 2022, 71(24): 244601. DOI: 10.7498/aps.71.20221340. [5] 程素秋, 陈高杰, 高鑫, 等. 不同装药战斗部壳体对水中兵器的爆炸威力 [J]. 爆炸与冲击, 2018, 38(6): 1372–1377. DOI: 10.11883/bzycj-2017-0127.CHENG S Q, CHEN G J, GAO X, et al. Estimation of underwater explosive energy for different charge warhead shells [J]. Explosion and Shock Waves, 2018, 38(6): 1372–1377. DOI: 10.11883/bzycj-2017-0127. [6] WALSH J M, SHREFFLER R G, WILLIG F J. Limiting conditions for jet formation in high velocity collisions [J]. Journal of Applied Physics, 1953, 24(3): 349–359. DOI: 10.1063/1.1721278. [7] ASAY J R, MIX L P, PERRY F C. Ejection of material from shocked surfaces [J]. Applied Physics Letters, 1976, 29(5): 284–287. DOI: 10.1063/1.89066. [8] ASAY J R. Material ejection from shock-loaded free surfaces of aluminum and lead: SAND-76-0542 [R]. Albuquerque, USA: Sandia National Laboratories, 1976. DOI: 10.2172/7136578. [9] ASAY J R. Thick-plate technique for measuring ejecta from shocked surface [J]. Journal of Applied Physics, 1978, 49(12): 6173–6175. DOI: 10.1063/1.324545. [10] SORENSON D S, MINICH R W, ROMERO J L, et al. Ejecta particle size distributions for shock loaded Sn and Al metals [J]. Journal of Applied Physics, 2002, 92(10): 5830–5836. DOI: 10.1063/1.1515125. [11] OGORODNIKOV V A, MIKHAĬOV A L, BURTSEV V V, et al. Detecting the ejection of particles from the free surface of a shock-loaded sample [J]. Journal of Experimental and Theoretical Physics, 2009, 109(3): 530–535. DOI: 10.1134/S1063776109090180. [12] ESCOBEDO J P, DENNIS-Koller D, CERRETA E K, et al. Effects of grain size and boundary structure on the dynamic tensile response of copper [J]. Journal of Applied Physics, 2011, 110(3): 033513. DOI: 10.1063/1.3607294. [13] FOWLES G R, ISBELL W M. Method for Hugoniot equation-of-state measurements at extreme pressures [J]. Journal of Applied Physics, 1965, 36(4): 1377–1379. DOI: 10.1063/1.1714313. [14] BROWN J L, RAVICHANDRAN G, REINHART W D, et al. High pressure Hugoniot measurements using converging shocks [J]. Journal of Applied Physics, 2011, 109(9): 093520. DOI: 10.1063/1.3590140. [15] THADHANI N N. Shock compression processing of powders [J]. Advanced Materials and Manufacturing Processes, 1988, 3(4): 493–549. DOI: 10.1080/10426918808953217. [16] 汤文辉. 冲击波物理[M]. 北京: 科学出版社, 2011: 174–181. [17] JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates, and high temperatures [C]// Proceedings of 7th International Symposium on Ballistics. The Hague, 1983: 541–547. [18] ANSYS AUTODYN Material [DB]. ANSYS, 2019. [19] 俞宇颖, 谭华, 胡建波, 等. 钽和LY12铝的高压声速测量 [J]. 爆炸与冲击, 2006, 26(6): 486–491. DOI: 10.11883/1001-1455(2006)06-0486-06.YU Y Y, TAN H, HU J B, et al. Measurements of sound velocities in shock-compressed tantalum and LY12 Al [J]. Explosion and Shock Waves, 2006, 26(6): 486–491. DOI: 10.11883/1001-1455(2006)06-0486-06. [20] WANG X M, SHI J. Validation of Johnson-Cook plasticity and damage model using impact experiment [J]. International Journal of Impact Engineering, 2013, 60: 67–75. DOI: 10.1016/j.ijimpeng.2013.04.010. [21] LEE S, BARTHELAT F, HUTCHHINSON J W, et al. Dynamic failure of metallic pyramidal truss core materials: experiments and modeling [J]. International Journal of Plasticity, 2006, 22(11): 2118–2145. DOI: 10.1016/j.ijplas.2006.02.006. [22] 侯日立, 周平, 彭建祥. 冲击波作用下LY12铝合金结构毁伤的数值模拟 [J]. 爆炸与冲击, 2012, 32(5): 470–474. DOI: 10.11883/1001-1455(2012)05-0470-05.HOU R L, ZHOU P, PENG J X. Numerical simulation of shock damage of LY12 aluminium alloy sructure [J]. Explosion and Shock Waves, 2012, 32(5): 470–474. DOI: 10.11883/1001-1455(2012)05-0470-05. [23] 辛春亮, 薛再清, 涂建, 等. 有限元分析常用材料参数书册 [M]. 北京: 机械工业出版社, 2020: 139. -