为了攻击机库、高级指挥要地等各种不同形式的综合防护目标, “智能型”的钻地弹药相继出现^[1]。弹体高速侵彻硬目标的加速度-时间历程曲线是一个非常重要的测试量, 对于考核钻地弹的结构强度, 检验弹体设计合理性, 考核战斗部装药的安定性, 研究弹载电子仪器的存活性、引信的工作可靠性, 评价防护工事的抗侵彻效果至关重要。在现有的侵彻加速度测试技术中, 弹载存储测试技术较先进^[2-4], 它能记录膛内、飞行穿靶几个过程的实时加速度。国外对于硬目标侵彻过程动态参数测试起步较早, 瑞士的武器系统和弹药实验中心设计了高g值弹道飞行数据记录器, 质量较轻(3.58 kg), 安装了3个6×10⁴ g量程的Endevco_7270A压阻加速度计, 成功测得最大加速度(9×10⁴ g)。国内中北大学电子测试技术国防科技重点实验室通过选购B & K公司的8309或兵器204所988加速度计, 并采用泡沫铝缓冲保护, 电路模块质量控制在150g以内, 测试范围最高可达2×10 ⁵ g。经过多年的研究, 高g值侵彻加速度现场测试技术已经基本解决了低速和中速条件下(800 m/s以下)“测不着”的问题, 正朝着提高高速侵彻条件下数据捕获率和进一步解决“测不准”问题的方向发展。弹载测试装置的工作环境恶劣, 侵彻原型实验需要耗费大量资源, 进行必要的数值模拟分析是对侵彻实验的必要补充。对于加速度信号测试的影响分析有利于提高数据捕获率, 能够进一步优化弹载测试装置的设计。

1.   测试原理及实验

加速度测试系统原理图如图 1所示, 由4部分组成:加速度传感器、存储记录部分、数据接口、数据处理软件。测试装置如图 2所示。

图  1  加速度测试系统原理图

Figure  1.  Schematic diagram of acceleration testing system

下载: 全尺寸图片 幻灯片

图  2  加速度测试装置

Figure  2.  Acceleration test equipment

下载: 全尺寸图片 幻灯片

实验时, 将加速度过载测试仪与配重块放入子弹中, 用压筒的一端压紧测试仪下端盖的上部平面, 另一端加配合垫片, 并用尾翼组合件拧紧到位, 保证加速度过载测试仪与弹体内腔前部平面完全贴和, 尽可能保证刚性连接。

弹丸侵彻混凝土靶板后靶板正面破坏情况如图 3所示, 可以看到靶体表面成坑、靶体表面的径向和环向裂纹破坏现象, 证明弹体在撞击靶体时具有良好的姿态。图 4是弹体侵彻混凝土靶过程的加速度曲线图。图 5为以3 k Hz截止频率滤波后的加速度曲线图, 峰值加速度为4.133×10⁴ g。对侵彻加速度曲线二次积分得到侵彻深度曲线(见图 6), 积分位移为0.562 m, 与实测侵深0.600 m相差6.3%。

图  3  弹丸侵彻混凝土靶板的破坏情况

Figure  3.  Damage of concrete target penetrated by projectile

下载: 全尺寸图片 幻灯片

图  4  侵彻过程的加速度曲线

Figure  4.  Acceleration curve during penetration

下载: 全尺寸图片 幻灯片

图  5  3 k Hz滤波后的加速度曲线

Figure  5.  Acceleration curve after 3 k Hz filter

下载: 全尺寸图片 幻灯片

图  6  侵彻深度曲线

Figure  6.  Penetration depth curve

下载: 全尺寸图片 幻灯片

2.   弹体侵彻混凝土靶过程的数值模拟

在适量的实验基础上基于LS-DYNA进行数值模拟, 是对侵彻实验的重要补充^[5-6]。

弹丸壳体是按照实验弹尺寸建立的, 内部结构作了较大的简化, 将测试仪和套筒处理成一等重的质量块装在结构前部, 结构后部则同样用一等重的质量块来代替套筒和弹尾翼。由于侵彻计算中主要关心壳体和测试仪的刚体加速度以及侵彻最大行程, 不关心各结构的受力情况, 因此以上简化具有一定的合理性。为了模拟半无限靶, 将实验用的混凝土靶板半径取足够大, 除了对称面和自由面外, 模型的边界条件是将靶板的周身加以非反射边界条件(应力波通过此边界传播而不产生反射效应)以模拟实际目标的边界效应。测试仪视为各向同性弹塑性模型, 本构模型(强度模型)采用Von Mises准则。弹内尾部的套筒可等同为一个弹性体, 材料模型采用弹性模型即MAT_ELASTIC。混凝土结构采用JHC模型, 该模型能够很好地描述高应变率条件下混凝土的响应问题。考虑到弹丸垂直侵彻靶板为对称问题, 在进行有限元分析时选取1/2实际模型建立计算模型。有限元模型网格如图 7所示。

图  7  有限元网格划分图

Figure  7.  Finite element mesh for generation

下载: 全尺寸图片 幻灯片

图 8给出了弹丸侵彻靶板的数值模拟结果。从计算结果来看, 混凝土靶表面成坑漏斗形状、侵彻过程转弯等同实验现象较吻合。图 9为弹体轴向刚体加速度随时间的变化曲线, 与图 5相比, 幅值与持续时间基本相同, 与实验结果较接近, 表明了数值模型建立及参数选择的正确性。

图  8  t=2 500μs侵彻过程的Von-Mises应力云图

Figure  8.  Von-Mises stress nephogram of penetration process

下载: 全尺寸图片 幻灯片

图  9  侵彻加速度-时间曲线

Figure  9.  Penetration acceleration vs time

下载: 全尺寸图片 幻灯片

由于应力波的传播, 弹体的不同部位在侵彻过程中的加速度是不一样的, 应力波的波前是加速度的峰值之处^[7], 但波前过后峰值很快下降, 所以峰值处的脉宽一般很窄。图 10为弹壳体A、测试仪B及套筒C的侵彻加速度。通过对加速度信号的滤波, 将高频部分滤除, 最后只剩下各部分的刚体加速度, 而此刚体加速度各部分相同, 如图 11所示。

图  10  3个点的侵彻加速度

Figure  10.  Penetration acceleration at three points

下载: 全尺寸图片 幻灯片

图  11  3个点的侵彻加速度（滤波）

Figure  11.  Penetration acceleration at three points(filter)

下载: 全尺寸图片 幻灯片

在弹体强度设计、战斗部装药安定分析等时, 弹体的整体加速度, 即刚体加速度是大家非常关心的^[8], 加速度测试仪在弹体内的安装部位和安装方式不会对侵彻过程中的刚体加速度测量造成影响, 这可以通过采用适当的滤波截止频率进行滤波而得到所需要的加速度信号。

弹体内腔的前端面加工较粗糙, 为了使加速度计与弹体之间绝缘, 在加速度计与安装座间通常加绝缘纸垫, 当应力波从弹头传播到弹体与测试装置的接触面以及加速度计安装座时, 大部分应力波被反射, 所以加速度计基本不受应力波的影响。

3.   测试装置基础运动对实测加速度信号的影响分析

图 12的加速度计安装方式可以简化为图 13所示的模型, 图中c为阻尼器的阻尼系数, m为运动质量块的质量, k为弹性元件的刚度系数, x(t)为质量块所受激励, 它可被视为一个基础激励单自由度二阶系统。设加速度计的激励信号为基础运动:u=u₀ sinωt, 基础加速度:a_f(t)=-u₀ω²sinωt, 其相对基础的运动方程为

图  12  加速度测试仪安装图

Figure  12.  Installation drawing of acceleration recorder

下载: 全尺寸图片 幻灯片

图  13  加速度计力学模型

Figure  13.  Mechanical model for accelerometer

下载: 全尺寸图片 幻灯片

整理得

式中:ω_n为无阻尼固有角频率, ζ为阻尼比。若忽略阻尼, 则式(2)为

对杜哈梅积分逐次应用分部积分, 可把相对位移表示成关于 $\frac{1}{\omega}$的幂级数, 因加速度传感器的自振频率很高, 取其前3项为

因为 $-\frac{m a_{\mathrm{f}}(t)}{k}=-\frac{a_{\mathrm{f}}(t)}{k / m}=-\frac{a_{\mathrm{f}}(t)}{\omega_{n}^{2}}$, 上式右边第1项为加速度计对惯性力-ma_f(t)静力响应, 其余各项为加速度计对惯性力-ma_f(t)的动力响应, 它是以ω为频率的高频信号。式(4)中的后2项在加速度实测信号中反映出弹体和安装结构的高频振动的影响。

4.   测试装置安装结构刚度对实测加速度信号的影响

下面结合图 14中的测试装置安装结构, 分析铝压筒和加速度计安装座的刚度对实测加速度信号的影响。

图  14  加速度测试仪结构示意图

Figure  14.  Structural diagram of acceleration recorder

下载: 全尺寸图片 幻灯片

(1) 当弹体在炮膛内加速向前运动时, 测试装置因惯性向后运动压缩铝筒, 铝压筒的刚度将影响膛内加速度信号, 甚至导致测试失败。图 15中曲线出现剧烈振荡。计算机模拟可知铝压筒的固有频率为1.6 k Hz, 而图 16中在1.5 k Hz左右出现一个尖峰, 即膛内加速度信号中包含铝压筒刚度的影响。按1.5 k Hz滤波, 可得到弹体在膛内运动时的刚体加速度, 如图 17所示。

图  15  膛内加速度曲线

Figure  15.  Acceleration-time curve in bore

下载: 全尺寸图片 幻灯片

图  16  膛内加速度频谱

Figure  16.  Frequency spectrum of acceleration in bore

下载: 全尺寸图片 幻灯片

图  17  弹体在膛内的刚体加速度

Figure  17.  Rigid-body acceleration of projectile in bore

下载: 全尺寸图片 幻灯片

在图 16的频谱图中, 23 k Hz附近出现一个谐振峰, 表明在该频率附近结构发生剧烈振动。采用ANSYS对弹体进行模态分析, 取整个结构的四分之一进行分析。图 18显示弹体第23阶模态(除零频外)对应的振型, 主要为测试装置上、下端盖的弯曲振动, 其振动方向与加速度计的敏感方向相同, 并且该频率(22.0 k Hz)接近加速度计的谐振频率(25 k Hz)。可以判断, 正是第23阶模态(除零频外)对应的振动使加速度计输出过大, 可能导致测试在膛内阶段就失败。

图  18  弹体第23阶模态的振型

Figure  18.  Projectile mode shape of the twenty-third order mode

下载: 全尺寸图片 幻灯片

若将铝筒换为钢筒, 并不能明显提高测试装置上、下端盖的弯曲振动频率, 其值为23.37 k Hz。

由自由杆纵向振动频率公式(5)可知, l(长度)一定, E(材料弹性模量)、m (质量)均增大, 频率增大不明显。

为了增加铝压筒刚度, 使其上段为实心, 质量增加0.74 kg, 整体装配后弹体重心后移2.5 mm, 对弹的飞行稳定性影响很小。对应测试装置上、下端盖的弯曲振动的频率(第13阶)为28.46 k Hz。另外, 还可将铝筒和尾翼间的调整垫片改为钢片, 以增加预压力, 提高该振动频率。采取提高铝压筒刚度的措施后, 且膛内加速度曲线振荡比以前明显降低。这就从一个方面提高了数据的质量。

(2) 当弹体侵彻靶板时, 测试装置相对弹体有向前运动的趋势, 加速度计安装座压紧弹体, 因为安装座本身的刚度大, 而加速度计是通过钢螺栓安装在基座上, 所以此时整体安装刚度大, 可以认为加速度的信号是对弹体内腔前端安装面的响应。

5.   加速度计安装方式对实测信号影响的分析

由实验可知加速度测试仪在弹体内的安装部位和安装方式不会对侵彻过程中的刚体加速度测量造成影响, 但加速度计的输出信号是其对安装点(面)加速度激励的响应。所以不同的安装方式对加速度计的输出信号有很大影响。下面对目前已采用或可能的安装方式进行讨论。

(1) 将加速度计通过钢螺栓直接安装在弹体上。这种安装刚度最大, 频响特性最好, 但无论安装在弹体内空腔的前端还是尾翼盖上, 传感器数据线与记录电路的连接都难以实现, 所以目前国内外都采用类似于图 14的结构。

(2) 将电路模块与加速度计安装成一个整体, 或将加速度计直接灌封在电路模块中, 然后再对这一整体进行缓冲保护, 目的是减短连接到电路模块上的加速度计数据线的长度, 提高数据线的抗高冲击能力。这种安装方式可简化为图 19所示的力学模型。

图  19  整体安装的力学模型

Figure  19.  Mechanical model for overall installation

下载: 全尺寸图片 幻灯片

在侵彻过程中, 整体相对弹体以相对加速度a_r向前运动, 加速度计感知的绝对加速度为a₀=a+a_r, 实际上并不是弹体的加速度。

若从冲击响应的角度考虑, 这个问题可简化为一个基础激励的二阶系统, 如图 20所示, 整体的运动方程为

图  20  基础激励的2阶系统

Figure  20.  Two-order system in base excitation

下载: 全尺寸图片 幻灯片

令, $a=\ddot{x} / \ddot{x}_{\max }, \tau=t / T, \delta=-k y /\left(m \ddot{x}_{\max }\right),R=2 \pi \sqrt{m / k} / T$。其中 $\ddot{x}_{\max }$为弹体加速度最大值, 为弹体加速度瞬时值, T为弹体加速度脉冲持续时间。

对上式进行量纲一化, 得到

式(7)为相对于量纲一输入a的量纲一响应δ, 从式(7)可知, 当R很小时, 方程的前2项可以忽略, 于是δ=a; 当R变大时, 前2项发生影响, 第1项引起在a值上下的振荡, 第2项引起时间滞后。这将使弹体加速度波形发生变化, 由该加速度信号积分得到的着靶速度、侵彻深度与实测值之间的误差也随之加大。

所以, 采用类似于图 14的结构将加速度计刚性安装, 目前来看是比较理想的。

6.   结论

(1) 洛阳水利工程技术研究所进行了过载特性的理论计算^[9], 计算结果与实验测量数据符合较好。数值计算结果与实验结果较吻合, 表明此数值模型建立及参数选择的正确性。测试数据为实测信号的影响分析提供了事实依据。

(2) 由于应力波的传播, 弹体的不同部位在侵彻过程中的加速度是不一样的, 应力波的波前是加速度的峰值之处, 但波前过后峰值很快下降, 所以峰值处的脉宽一般很窄。在进行加速度测试时, 可根据不同的目的(如各部件的抗冲击能力等), 将加速度计安装在需要测试的位置。通过在加速度计安装座下加绝缘纸, 既可使加速度计与弹体绝缘, 又可大幅减轻应力波影响。

(3) 加速度测试仪在弹体内的安装部位和安装方式不会对侵彻过程中的刚体加速度测量造成影响, 可以通过采用适当的滤波截止频率进行滤波而得到所需要的加速度信号。加速度计的输出信号是其对安装点(面)加速度激励的响应, 不同的安装方式, 对加速度计的输出信号影响很大。在具体安装时, 加速度计应采用刚性连接, 同时应提高测试装置的安装刚度。