侵彻弹体结构纵向振动频率特性分析*

刘 波,杨黎明,李东杰,欧阳科

(中国工程物理研究院电子工程研究所，四川 绵阳 621999)

为了更好地使智能化武器达到最大的毁伤效果，就希望引信能达到“计层”与“计深”的功能，通常要求通过加速度传感器配合信号处理组件来测试处理弹体侵彻过程中的加速度信号[1-2]。此外，评价弹体内部的电子组件、装药等能否经受得住侵彻过程的高冲击环境，也需要对弹体穿靶过程中的加速度历程进行现场测试记录[3]。在现有的侵彻加速度现场测试技术中，较多采用弹载存储测试技术，由加速度传感器、记录电路模块和电池总共3 部分组成。其中，加速度传感器的敏感元件将自身受到的力信号转换为电信号，再由记录电路模块将电信号记录下来，电池为加速度传感器和记录电路模块提供电源[4-5]。

加速度传感器在测试弹体穿靶过程的加速度历程，不仅测得由于目标阻碍引起的弹体过载加速度，还会测到由于应力波传播引起的弹体结构振动的加速度[6]。在高冲击加速度测试中，弹体结构振动信号甚至强于弹体过载加速度信号。而在引信的“计层”与“计深”算法设计中，需要从测试的加速度信号中提取弹体的过载加速度信号。要准确测试弹体的过载加速度，就必须对弹体的振动特性有所了解。Hao等[7]、黄家蓉等[8]对弹体的振动特性和弹体侵彻加速度信号进行了分析，但都缺乏与简单解析解的对比，而且侵彻加速度信号也只是记录侵彻单层靶的数据。

本文中，将弹体简化为简单的均匀长直圆杆模型，对弹体的固有频率进行理论分析。此外，利用有限元软件ANSYS，对简化的弹体有限元模型进行结构模态分析、谐响应分析，并联系加速度传感器测试数据，进一步分析侵彻过程中侵彻弹体的振动特性[9]。

1 弹体振动理论分析

研究由于弹性应力波在侵彻弹体内传播所导致的结构振动时，可以忽略头部的局部弹性变形，将弹体简化为一个圆柱形均匀长直杆[10]。一维弹性应力波在杆内的传播方程为：

(1)

令U=X(x)T(t)，对一维应力波波动方程分离变量，得到：

Xn=Asin(ωnx)+Bcos(ωnx)

(2)

代入边界条件，得到简化模型的固有角频率为：

ωn=nπc0/L

(3)

弹体简化的长圆杆模型的纵向振动固有频率为：

fn=nc0/(2L)

(4)

式中：L为简析模型的长度。

2 弹体模态有限元模拟

利用ANSYS建立弹体有限元模型(见图1)，弹体全长为2.1 m，外径为0.38 m。弹体由弹壳、装药及压紧环3部分组成。弹体3部分材料都采用线弹性本构模型，材料参数见表1。

表1 材料模型参数Table 1 Material parameter in simulation

将弹体简化为均匀长直圆杆(杆材料采用弹壳的材料模型，杆长为弹长)，计算圆杆纵向振动的固有频率，并对比有限元模型模拟结果，见表2。然后提取有限元模型对应的弹体纵向的前四阶振型，如图2所示。由表2可以看出，有限元模型的前两阶固有频率模拟结果与简化的圆杆模型的近似值偏差都在5%以内，但是高阶固有频率的简析近似值与模拟结果还是有10%左右的差距。说明利用圆杆纵向振动波动方程求解纵向振动固有频率是较准确的。

表2 模态分析与谐响应分析结果与理论值的比较Table 2 Comparison of simulation results with theoretical values

3 弹体结构谐响应模拟分析

利用ANSYS建立有限元模型，结构与材料参数同上。对有限元模型进行瞬态分析，约束弹尖节点Y、Z方向的位移，在弹体头部施加冲击加速度载荷。设定谐响应频率范围为0～5 000 Hz，分为250个子步进行加载。

实际的侵彻测试实验中，一般由引信里的加速度测试记录装置现场测试侵彻过程中的加速度历程。而引信一般装于弹体尾部的压紧环内，提取模拟分析得到压紧环中心处对应节点的频率曲线，如图3所示。理论计算得到的、模态分析得到的和谐响应分析得到的弹体纵向振动固有频率，见表2。从表2可以看出，加速度测试装置在频率为1 200、2 360、3 460、4 380 Hz的简谐载荷下，位移响应最大。谐响应分析得到的前3阶固有频率与简析近似解偏差均不超过7%，再次证实了利用波动方程求解弹体振动固有频率的准确性。

4 实测加速度信号分析

图4为侵彻5层钢筋混凝土靶板的加速度传感器实测得到的穿靶过程中的加速度信号，利用MATLAB对加速度信号进行Fourier 变换，得到加速度信号的功率谱，如图5所示。从图5可以看出，加速度信号功率谱在1 114 Hz能量达到最大值，与前面模态分析的1 198 Hz比较接近，实际的测试数据也说明了利用波动方程求解弹体振动固有频率的有效性。在实际工况中，压紧环与弹体之间采用螺纹连接，而不是建模中直接接触的方式，实际的螺纹连接的刚度相对建模理偏小，导致实测的加速度信号的频率谐振点偏小。

但是测试值与模态分析、谐响应分析结果还是有8%左右的误差。这可能是：在实际的侵彻加速度历程测试中，测试装置一般螺纹连接装于压紧环内；而在有限元建模中，进行了简化处理，直接将测试装置与压紧环建模为一个整体，增加了测试装置与压紧环之间的连接刚度，导致有限元模型的结果偏高。

为了进一步了解实测加速度信号的频率特性，对图4中传感器曲线进行滤波分析。利用MATLAB中的小波分析工具对过载曲线进行简单的滤波分析，分析它在不同频段内的信号成分。对信号进行5级分解，图6为低通滤波结果，表示信号接近程度，而图7为带通滤波结果，表示信号细节部分。

从图6可以看出，信号在第4层(见图6(e))最理想，对实测加速度信号进行625 Hz低通滤波处理后，加速度信号明显显露出5个尖峰，对应弹体侵彻5层混凝土靶的加速度冲击；而在第5层(见图6(f))，对实测加速度信号进行312.5 Hz低通滤波处理后，加速度信号不能显示出明显的5层。从图7可以看出：在第1、2级高频噪声频段，幅值比较小，且信号比较集中，贯穿整个侵彻过程，在进行信号分析中可以忽略不计；在第3、4级频段内的信号幅值比第1、2级频段内的信号大，且其信号充满整个侵彻过程，符合弹丸振动信号特征。因此，可认为在第3、4级频段内的信号即为弹丸振动信号。图6中第4层与图7中第3、4级幅度的比较，验证了受应力波影响严重的信号中弹丸的刚体过载淹没在高频振荡信号中。

对信号进行滤波处理就是提取不同频段内的信号，则各层分解就是对信号进行低通或带通滤波的过程。根据上面的分析结果，可以确定图7中信号高频振荡信号的频率范围在625～2 500 Hz，也可以确定信号低通滤波处理的最佳截止频率，与理论计算得到1 235 Hz比较接近。

5 结 论

通过将弹体简化为圆柱形长直杆，理论分析了弹体振动的频率特性。建立有限元模型，对弹体进行了模态分析与谐响应分析，弹体振动在1 2 00Hz出现峰值，与简化模型的理论值吻合较好，前两阶固有频率的误差都在5%以内。并对实际测试的加速度传感器信号进行Fourier 变换和小波分析，信号功率谱在1 114 Hz能量最大，与模拟分析得到的固有频率有8%左右的差距。模拟分析中将压紧环建模为一个整体，增加了测试装置与压紧环之间的连接刚度，导致有限元模型的模拟结果偏高。加速度信号低通滤波结果也和理论计算值很好吻合。这些说明，利用波动方程求解侵彻弹体纵向振动固有频率具有有效性。

参考文献：

[1] 赵海峰,张亚,李世忠,等.侵彻弹体频率特性分析及过载信号处理[J].中国机械工程,2015,26(22):3034-3039.

ZHAO Haifeng, ZHANG Ya, LI Shizhong,et al. Frequency characteristics analysis and overload signal process of penetration projectile[J]. China Mechanical Engineering, 2015,26(22):3034-3039.

[2] 宋萍,李科杰.硬目标侵彻武器高冲击试验和高过载传感器技术国外发展概况[J].测控技术,2002,21(增刊1):30-32.

SONG Ping, LI Kejie. Development of the technique of high shock test on hard-target penetration weapons and high-gaccelerometer at broad[J]. Measurement & Control Technology, 2002,21(suppl 1):30-32.

[3] RASHKOVSKA A, TROBEC R, DEPOLLI M, et al. Time and frequency analysis of non-uniform sampling[C]∥MIPRO, Proceedings of the 33rd International Convention. 2010:236-240.

[4] 徐鹏,祖静,范锦彪.高g值侵彻加速度测试及其相关技术研究进展[J].兵工学报,2011,32(6):739-745.

XU Peng, ZU Jing, FAN Jinbiao. Research development of highgpenetration acceleration test and its correlative technology[J]. Acta Armamentarii, 2011,32(6):739-745.

[5] ZHANG Wendong, CHEN Lujiang, XIONG Jijun. Ultra-high g deceleration-time measurement for the penetration into steel target[J]. International Journal of Impact Engineering, 2007,34(3):436-447.

[6] 赵小龙,马铁华,徐鹏,等.弹丸侵彻混凝土加速度信号测试及分析[J].爆炸与冲击,2014,34(3):347-353.

ZHAO Xiaolong, MA Tiehua, XU Peng, et al. Acceleration signal test and analysis for projectile penetrating into concrete[J]. Explosion and Shock Waves, 2014,34(3):347-353.

[7] HAO H, LI X, SUN Y. Projectile structural response frequency characteristics analysis method in penetration process[J]. Journal of Vibration Measurement & Diagnosis, 2013,1997(16):7-11.

[8] 黄家蓉,刘瑞朝,何翔,等.侵彻过载测试信号的数据处理方法[J].爆炸与冲击,2009,29(5):555-560.

HUANG Jiarong, LIU Ruichao, HE Xiang, et al. A new data processing technique for measured penentration overloads[J]. Explosion and Shock Waves, 2009,29(5):555-560.

[9] ZHANG Q. Applications of FEM and lanczos algorithm in modal analysis of slender missile[J]. Journal of Projectile Rocket & Guidance, 2007,27(4):61-63.

[10] 王秋萍,张虎,吴迎春.金属杆中纵固有频率的研究[J].声学与电子工程,2005,24(8):144-147.

WANG Qiuping, ZHANG Hu, WU Yingchun. Study of longitudinal vibrating natural frequency of metallic rod[J]. Acoustics and Electronics Engineering, 2005,24(8):144-147.