压电双晶片在制导子弹中的应用

程罡; 陈超; 陈恒

doi:10.15892/j.cnki.djzdxb.2019.06.031

弹箭与制导学报 >

2019 , Vol. 39 >Issue 6: 142 - 145

DOI: https://doi.org/10.15892/j.cnki.djzdxb.2019.06.031

压电双晶片在制导子弹中的应用

程罡 ,
陈超 ,
陈恒

展开

南京航空航天大学机械结构力学及控制国家重点实验室, 南京 210016

程罡(1993-),男,河南信阳人,硕士研究生,研究方向:压电梁在飞行器中的应用。

收稿日期: 2018-11-20

网络出版日期: 2025-05-12

基金资助

国家自然科学基金(51575259)

中央高校基本科研基金(NS2015004)

收起

Application of Piezoelectric Bimorph in Guided Bullets

CHENG Gang ,
CHEN Chao ,
CHEN Heng

Expand

State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

Received date: 2018-11-20

Online published: 2025-05-12

Fold

摘要

为了实现对微小口径高旋弹药的精确控制,延伸压电双晶片的应用领域,基于变质心控制机制,提出了一种应用于微小口径高旋稳定弹的调姿态装置。该方案通过压电双晶片摆动质量块来产生质量矩的变化,基于离心场作用达到输出力矩的放大。建立了双晶片-质量块系统的结构动力学模型,并进行了有限元仿真分析。建立了弹头整机(弹头部分)的刚体动力学模型。理论分析结果表明:该调姿态系统能够满足姿态调整的需求,并能够实现输出力矩的放大。

关键词： 压电双晶片; 变质心; 高旋子弹; 姿态控制

本文引用格式

程罡 , 陈超 , 陈恒 . 压电双晶片在制导子弹中的应用[J]. 弹箭与制导学报, 2019 , 39(6) : 142 -145 . DOI: 10.15892/j.cnki.djzdxb.2019.06.031

Abstract

In order to realize the precise control of micro-caliber high-elastic ammunition,the application field of extended piezoelectric bimorph is proposed.Based on the variable centroid control mechanism,a tuning attitude device for small-caliber high-rotation stable missile is proposed.The solution generates a change in the mass moment by the piezoelectric bimorph oscillating mass,and the amplification of the output torque is achieved based on the centrifugal field action.The structural dynamics model of the dual-wafer-mass system was established and the finite element simulation analysis was carried out.A rigid body dynamics model of the warhead (the warhead part) was established.The theoretical analysis results show that the attitude system can meet the requirements of attitude adjustment and can amplify the output torque.

Key words： piezoelectric bimorph; moving mass; high-speed spinning bullets; adjustment control

0 引言

随着军事科技的发展,“精确打击”的需求已愈来愈大,传统非制导弹药也面临着信息化和智能化的升级,因此常规子弹药的智能化已成为了一种趋势,尤其是基于身管发射平台的具备高旋特性的微小口径子弹药,由于其飞行时间短,可利用空间小,其制导和控制的设计面临着不同于常规制导炮弹的技术难题和挑战,变质心控制技术作为再入飞行器的一种新兴控制技术,非自旋弹头的变质心控制理论研究已经比较成熟,但对于自旋的飞行器变质心控制研究成果相对较少^[1-5]。文献[6]提出了一种基于直线驱动装置的双质量块滑模系统的变质心调控方式,该方式虽然将变质心成功应用于高速自旋弹头上,但其滑动响应频率受到限制,并不能在超高速飞行弹药上满足相应要求。

因此研究一种控制装置使其拥有更高的响应及控制带宽并且能够保持高旋弹体的良好气动外形是具有重大意义的。压电双晶片具有响应快、带宽高、分辨率高的特点,文中采用了双晶片带动质量块系统的控制装置的配置方式,利用离心场的作用对输出力矩进行了放大,弥补了传统方案的不足。利用哈密尔顿原理和拉格朗日方程分别建立了双晶片-质量块动力学模型以及整机动力学方程。对双晶片质量块变质心控制系统进行了仿真并分析了控制装置的性能。

1 压电双晶片控制装置及工作原理

文中研究的基于压电双晶片-变质心系统的控制装置包括一套双晶片质量块驱动系统,可放入12 mm口径弹药的压电双晶片驱动器,图1为子弹整体示意图,可应用于一般枪管,无需另行设计配套枪管。弹头内部配置方式如图2所示。

显示原图|下载原图ZIP|生成PPT

图1 制导子弹整体示意图

显示原图|下载原图ZIP|生成PPT

图2 压电双晶片变质心调姿系统配置示意图

在图2中双晶片质量块装置通过轴承与外部结构连接,当外部高速旋转时,内部双晶片结构并不随着进行高速旋转,并可以独立控制双晶片结构旋转以达到全方位调控和快速的频响。高旋弹体本身(除去内部压电双晶片质量块装置)就是一个绕其自转轴高速旋转的刚体,由于旋转的质量块做变质心运动,产生的离心力会增大变质心力矩的输出,依靠陀螺效应能够保证其飞行稳定性,图3为该变质心调节系统工作原理。

显示原图|下载原图ZIP|生成PPT

图3 双晶片-质量块系统工作原理示意图

双晶片控制质量块m沿bd方向做往复运动,子弹以角速度k转动,子弹转动周期一致,因此一个周期内,子弹会沿某个方向做持续摆动,气动力会产生一个持续的升力来让子弹转向(图中t₁和t₃时刻产生升力)。

2 双晶片-质量块系统动力学建模及仿真^[7]

考虑到变质心运动会产生离心力,故建模时考虑了离心力做功。只建立双晶片结构的动力学模型,暂不考虑陀螺进动效应。

2.1 静力学微分方程的建立

将执行机构简化为图4所示的自由端连接一质量块的悬臂梁,图中F₁为等效离心力分布于梁的各个单元。要建立静力学微分方程,首先建立动力学方程。

显示原图|下载原图ZIP|生成PPT

图4 固支双晶片-质量块模型简图

质量块动能:

(1) T_m=

$\frac{1}{2}$

$\frac{\partial w}{\partial t}$

)²+

$\frac{1}{2}$

$\left(\frac{{\partial }^{2}w}{\partial x\partial t}|{ }_{x=l}\right)$

(2)

$\begin{array}{l}{W}_{r}={\rho }_{s}{A}_{s}{\int }_{0}^{w}{\int }_{0}^{l}w{k}^{2}dxdy+\\ 2{\rho }_{p}{A}_{p}{\int }_{0}^{w}{\int }_{0}^{l}w{k}^{2}dxdy+{\int }_{0}^{w { }_{l}}m{w}_{l}{k}^{2}dx\end{array}$

(3) δW=-

$\frac{1}{2}$

Mδ

$\frac{\partial w}{\partial x}\left|{ }_{0}^{l}+\rho A{k}^{2}{\int }_{0}^{l}w\delta wdx+m{w}_{l}{k}^{2}\right.$

由哈密顿原理:

(4) ${\int }_{{t}_{1}}^{t { }_{2}} \left(\delta T-\delta U+\delta W\right)$dt=0

得到运动方程:

(5)ρA $\frac{{\partial }^{2}w}{\partial {t}^{2}}$+EI $\frac{{\partial }^{4}w}{\partial {x}^{4}}$-ρAk²w=0

该模型边界条件:

(6)

$\left\{\begin{array}{l}J\frac{{\partial }^{2}w}{\partial {t}^{2}}+EI\frac{{\partial }^{2}w}{\partial {x}^{2}}+M=0, x=l　\\ m\frac{{\partial }^{2}w}{\partial {x}^{2}}+EI\frac{{\partial }^{3}w}{\partial {x}^{3}}+mw{k}^{2}=0, x=l\\ \frac{\partial w}{\partial x}=0, x=0　　　　　　　\\ w=0, x=0　　　　　　　\end{array}\right.$

式(1)~式(6)中:ρA为双晶片线密度;J为质量块绕质心转动惯量;EI为双晶片弯曲刚度;M为电压产生的力矩;w为双晶片挠度;k为双晶片-质量块系统的转动角速度;m为质量块质量。

令式(5)中

$\frac{{\partial }^{2}w}{\partial {t}^{2}}$

=0,得到静力学微分方程:

(7)EI $\frac{{\partial }^{4}w}{\partial {x}^{4}}$-ρAk²w=0

其通解为:

(8) w=Ae^λx+Be^-^λx+Ccosλx+Dsinλx

边界调节为:

(9)

$\left\{\begin{array}{l}EI\frac{{\partial }^{2}w}{\partial {t}^{2}}+{M}_{e}=0, x=l　\\ EI\frac{{\partial }^{3}w}{\partial {x}^{3}}+mw{k}^{2}=0, x=l\\ \frac{\partial w}{\partial x}=0, x=0　　　　\\ w=0, x=0　　　　\end{array}\right.$

2.2 有限元仿真分析

利用ANSYS对压电双晶片质量块系统进行了仿真分析,图5为加载电压为100 V,转速为150 rad/s时的位移云图。进而分析了压电双晶片在20 V、60 V、100 V电压环境下,随着转速的变化压电双晶片输出位移的变化,并与MATLAB数值解进行对比,得到结果如图6。

显示原图|下载原图ZIP|生成PPT

图5 压电双晶片有限元位移云图

显示原图|下载原图ZIP|生成PPT

图6 不同电压下转速-输出位移关系图

由图6可知,随着输入电压的增加,双晶片输出变形有所增加,但双晶片本身的变形小,高转速能够增加离心力进而增大双晶片的变形,从而增大变质心产生的力矩,改变外弹体所受到的升力。

3 整机动力学方程的建立^[8]

由于弹壳(即子弹弹头除去内部可移动质量块的剩余部分)与内部质量块之间的连接属于刚性连接,故该系统的约束是理想约束,可以利用拉格朗日第二类方程建立系统的运动方程。

建立多刚体系统的定坐标系O,弹体坐标系O₁,动坐标系O₂。其中O取地面坐标系;弹体坐标系O₁的坐标原点在弹壳质心处;取与地面坐标系完全平行,原点在弹壳质心处的坐标为动坐标系O₂。

子弹弹头的动能由两部分组成:

①弹壳部分动能:

(10) T=ρ∫_v(r^T)'rdv=T₁+T₂

式中:T₁为平均动能,T₂为转动动能。

(11)T₁= $\frac{1}{2}$M·((x²)'+(y²)'+(z²)')

(12)

$\begin{array}{c}{T}_{2}=\frac{1}{2}\left[{\omega }_{{x}_{1}}　{\omega }_{{y}_{1}}　{\omega }_{{z}_{1}}\right]·\left[\begin{array}{lll}{I}_{x}& & \\ & {I}_{y}& \\ & & {I}_{z}\end{array}\right]·\left[\begin{array}{l}{\omega }_{{x}_{1}}\\ {\omega }_{{y}_{1}}\\ {\omega }_{{z}_{1}}\end{array}\right]\end{array}$

②质量块动能:

(13)

$\begin{array}{l}{T}_{m}={\rho }_{1}{\int }_{{v}_{1}}({r}^{T})\text{'}{r}_{1}d{v}_{1}=\\ 　\frac{1}{2}m\frac{d}{dt}\left[{x}_{p}+{y}_{p}+{z}_{p}\right]·\frac{d}{dt}\left\{\left[\begin{array}{l}{x}_{p}\\ {y}_{p}\\ {z}_{p}\end{array}\right]+\left[\begin{array}{l}X\\ Y\\ Z\end{array}\right]\right\}+\\ \frac{1}{2}m\frac{d}{dt}\left[X,Y,Z\right]·\frac{d}{dt}\left\{\left[\begin{array}{l}{x}_{p}\\ {y}_{p}\\ {z}_{p}\end{array}\right]+\left[\begin{array}{l}X\\ Y\\ Z\end{array}\right]\right\}　\end{array}$

式(10)~式(13)中:M为弹壳质量;m为质量块质量;r为弹壳质心矢径;r₁为质量块矢径;(x,y,z)为弹壳质心在坐标系O的值;x₁、y₁、z₁分别为滚动轴、偏航轴、俯仰轴、(x_p,y_p,z_p)为m的质心在O₂中的值。

利用拉格朗日方程建立系统的运动方程:

(14)

$\begin{array}{l}M·\left[\begin{array}{l}X″\\ Y″\\ Z″\end{array}\right]+\frac{1}{2}m·\left[\begin{array}{l}D\text{'}{ }_{x}\\ D\text{'}{ }_{y}\\ D\text{'}{ }_{z}\end{array}\right]+m·\left[\begin{array}{l}X″\\ Y″\\ Z″\end{array}\right]=\\ {F}_{Ox}+{F}_{Oy}+\left(M+m\right)g+{F}_{Oz}\end{array}$

(15)M·X″+ $\frac{1}{2}$m·D'_x+m·X″=F_Ox

(16)M·Y″+ $\frac{1}{2}$m·D'_y+m·Y″=(M+m)g+F_Oy

(17)M·Z″+ $\frac{1}{2}$m·D″_z+m·Z″=F_Oz

(18)

$\begin{array}{c}\left(\frac{dA}{dt}·\left[\begin{array}{l}\stackrel{}{\phi ″}\\ \stackrel{}{\psi ″}\\ \stackrel{}{\gamma ″}\end{array}\right]+A·\left[\begin{array}{l}\stackrel{}{\phi \text{'}}\\ \stackrel{}{\psi \text{'}}\\ \stackrel{}{\gamma \text{'}}\end{array}\right]\right)-m·\frac{d}{dt}\{\left[x{\text{'}}_{p},y{\text{'}}_{p},z{\text{'}}_{p}\right]·U·R\}={M}_{\phi }\end{array}$

(19)

$\begin{array}{c}\left[\begin{array}{lll}0& 1& 0\end{array}\right]·\left(\frac{dA}{dt}·\left[\begin{array}{l}\phi ″\\ \psi ″\\ \gamma ″\end{array}\right]+A·\left[\begin{array}{l}\phi \text{'}\\ \psi \text{'}\\ \gamma \text{'}\end{array}\right]\right)-m·\frac{d}{dt}\{\left[\begin{array}{lll}x{\text{'}}_{p}& y{\text{'}}_{p}& z{\text{'}}_{p}\end{array}\right]·O·R\}={M}_{\psi }\end{array}$

(20)

$\begin{array}{c}\left(\frac{dA}{dt}·\left[\begin{array}{l}\phi \text{'}\\ \psi \text{'}\\ \gamma \text{'}\end{array}\right]+A·\left[\begin{array}{l}\phi ″\\ \psi ″\\ \gamma ″\end{array}\right]\right)-m·\frac{d}{dt}\left\{\right[x{\text{'}}_{p},y{\text{'}}_{p},z{\text{'}}_{p}]·Y·R\}={M}_{\gamma }\end{array}$

式(14)~式(20)中:(x'_p,y'_p,z'_p)为m的质心在坐标系O₁的位置;φ为俯仰角;ψ为偏航角;γ为滚动角。求出上述6个公式中广义力,因为相对速度坐标系的阻力Q,升力L和侧向力H易求出,因此通过速度坐标系向惯性坐标系转换可得对应于X、Y、Z的广义力F_Ox、F_Oy、F_Oz。除此之外,为求对应于φ、ψ、γ的广义力矩M_φ、M_ψ、M_γ,应先求相对广义弹体坐标的气动力矩,结果如下:

(21)

$\begin{array}{l}{M}_{Z}={M}_{Z0}+{M}_{Z}^{\alpha }·\alpha +{M}_{Z}^{{\omega }_{Z}}·{\omega }_{Z}+{M}_{Z}^{\stackrel{·}{\alpha }}·\stackrel{}{\stackrel{·}{\alpha }}+{F}_{X}·{\delta }_{Y}\\ {M}_{Y}={M}_{Y}^{\beta }·\beta +{M}_{Y}^{{\omega }_{Y}}·{\omega }_{Y}+{M}_{Y}^{\stackrel{·}{\beta }}·\stackrel{·}{\beta }+{F}_{X}·{\delta }_{Z}\\ {M}_{X}={F}_{Z}·{\delta }_{Y}+{F}_{Y}·{\delta }_{Z}\end{array}$

通过坐标旋转得到:

(22)

$\begin{array}{l}{M}_{\phi }={M}_{OZ}\\ {M}_{\psi }=-{M}_{OX}·sin\psi +{M}_{OY}·cos\phi \\ {M}_{\gamma }={M}_{OX}·cos\phi cos\psi +\\ {M}_{OY}·sin\phi cos\psi -{M}_{OZ}·sin\psi \end{array}$

至此,变质心控制子弹弹头整机的空间六自由度动力学方程建立完成。

4 总结

1)基于压电双晶片变质心系统提出了一种新型姿态调节方案,充分利用了压电双晶片响应快的特点,同时增大了输出位移和输出力矩,可放入12.7 mm口径弹体内。

2)结构动力学分析和仿真分析结果表明在相同输入电压下,随着变质心系统转速的增加,压电双晶片输出变形不断增加。变质心系统输出力矩增加,充分验证了其用于高速自旋飞行器姿态调整的可行性。

文中提出的新型变质心姿态调控装置可以为内部空间狭小,高速自旋,控制带宽高位特征的微小飞行器的姿态调控装置提供参考。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

[1]	邓以高, 田军挺, 王亚锋, 等. 飞行器姿态控制方法综述[J]. 战术导弹控制技术, 2006(2) :7-13.

[2]	魏泽宇, 冯秀莲, 刘山, 等. 制导子弹的发展现状与关键技术[J]. 飞航导弹, 2015(10):38-42.

[3]	BARRETT R, LEE G. Design and testing of piezoelectric fli ght control actuators for hard-launch munitions[C]// FLATAU A B.Proceedings of SPIE Vol.5390 (Smart Structures and Materials 2004:Smart Structures and Integrated Systems).[S.l.]:SPIE,2004:459-470.]

[4]	HU K M, WEN L H. Scheme design and mechanical analysis about new-type bullet based piezoelectric rudder actuator[J]. Journal of Mechanical Engineering, 2015, 51(21):104-112.

[5]	ZHAO H C, WANG T, FAN S L. Modeling and adaptive sliding mode control of varyin gcentroid spinning warhead[J]. Science & Technology Reviem, 2009(5):43-46.

[6]	ZHAO H C, SHI X J, FAN S L. The variable centroid control of a spinning warhead based on a new collocating mode[J]. Aerospace Control, 2011, 29(4):3-7.

[7]	Zhao Cs. Ultrasonic motors technologies and application[M]. Beijing: science press, 2007.

[8]	YI Y, ZHOU F Q, ZHOU J. Motion analysis on moving center of mass controlled missile[J]. Aerospace Control, 2000, 18(3):1-5.

Options

文章导航

模态框（Modal）标题

摘要

本文引用格式

Abstract

0 引言

1 压电双晶片控制装置及工作原理

图1 制导子弹整体示意图

图2 压电双晶片变质心调姿系统配置示意图

图3 双晶片-质量块系统工作原理示意图

2 双晶片-质量块系统动力学建模及仿真[7]

2.1 静力学微分方程的建立

图4 固支双晶片-质量块模型简图

2.2 有限元仿真分析

图5 压电双晶片有限元位移云图

图6 不同电压下转速-输出位移关系图

3 整机动力学方程的建立[8]

4 总结

参考文献

2 双晶片-质量块系统动力学建模及仿真^[7]

3 整机动力学方程的建立^[8]