2-4 气动力导数及导弹控制动力系数

本节定位:连接理论力学控制工程的核心枢纽。将风洞试验得到的非线性气动力数据,通过小扰动线性化转化为无量纲动力系数,使成熟的线性定常控制理论可直接应用于导弹设计。


一、知识图谱

气动力建模

线性化前提

气动力导数

动力系数

控制设计

小扰动假设

系数固定

升力导数

力矩导数

阻尼导数

转动系数

平动系数

滚转系数

频域设计

根轨迹法

最优控制

本节的核心任务是将非线性的气动力与气动力矩函数线性化。在导弹控制系统设计中,为直接利用最成熟的线性定常系统控制理论,必须引入小扰动线性化系数固定两大假设,从而将时变非线性微分方程组转化为小扰动线性定常方程组,进而方便地使用频率法、根轨迹、最优控制及鲁棒控制等方法进行设计。


二、总体结构图

非线性气动力

特征点选取

小扰动线性化

气动力导数

动力系数

控制设计

马赫数固定

攻角固定

泰勒展开

保留一阶

线性化过程遵循严格的数学与工程逻辑。首先,在导弹飞行弹道中选取若干特征点,每个特征点对应特定的马赫数、飞行高度与攻角。其次,假设导弹在特征点附近受到小扰动,气动力与力矩函数可在该点进行泰勒展开,仅保留一阶项,得到线性增量关系。最后,将气动力导数除以质量速度乘积或转动惯量,得到无量纲化的动力系数。


三、各模块设计图

平动义

俯仰义

滚转义

平动系

俯仰系

滚转系

运算

法向导

俯仰导

滚转导

来源

滚转力矩

俯仰力矩

升力侧力

滚舵偏导

滚转阻尼

攻角导数

舵偏导数

阻尼导数

法向攻导

除滚转惯

除俯仰惯

除质速积

c德尔塔

c欧米伽

a阿尔法

a德尔塔

a欧米伽

b阿尔法

b德尔塔

滚舵控制

滚转阻尼

静稳度量

舵效大小

阻尼强弱

机动效率

法向贡献

3.1 气动力导数定义

气动力导数是气动力与力矩对各运动变量的一阶偏导数,在特征点上取为常数。以滚转通道为例,滚转力矩 MxM_xMx 对滚转舵偏转角 δx\delta_xδx 的导数定义为
Mxδx=∂Mx∂δxM_x^{\delta_x} = \frac{\partial M_x}{\partial \delta_x}Mxδx=δxMx

小增量 Δδx\Delta \delta_xΔδx 引起的力矩增量为
ΔMx=Mxδx⋅Δδx\Delta M_x = M_x^{\delta_x} \cdot \Delta \delta_xΔMx=MxδxΔδx

该导数的值与特征点的选取密切相关。由于舵偏角增量通常仅为几度,可将 MxδxM_x^{\delta_x}Mxδx 在选定的特征点上看作常数。滚转阻尼导数 MxωxM_x^{\omega_x}Mxωx 的量纲为力矩除以单位滚转角速度,其符号恒为负,因为阻尼力矩总是阻止旋转运动。

对于俯仰与偏航通道,由攻角引起的升力 YYY 通常表示为
Y=12ρV2SCyY = \frac{1}{2} \rho V^2 S C_yY=21ρV2SCy
式中 ρ\rhoρ 为大气密度,SSS 为弹体特征面积(通常取横截面积),Cy(Ma,α,δz)C_y(Ma, \alpha, \delta_z)Cy(Ma,α,δz) 为升力系数,在给定马赫数下是攻角与舵偏角的函数。

定义以下升力导数:
Yα=∂Y∂α=Cyα⋅12ρV2SY^\alpha = \frac{\partial Y}{\partial \alpha} = C_y^\alpha \cdot \frac{1}{2} \rho V^2 SYα=αY=Cyα21ρV2S
Yδz=∂Y∂δz=Cyδz⋅12ρV2SY^{\delta_z} = \frac{\partial Y}{\partial \delta_z} = C_y^{\delta_z} \cdot \frac{1}{2} \rho V^2 SYδz=δzY=Cyδz21ρV2S

对于对称布局导弹,侧力系数 Cz(Ma,β,δy)C_z(Ma, \beta, \delta_y)Cz(Ma,β,δy) 与升力系数 CyC_yCy 相等,侧力导数 ZβZ^\betaZβZδyZ^{\delta_y}Zδy 的形式与升力导数完全对称。

俯仰力矩 MzM_zMz 的表达式为
Mz=12ρV2SLmzM_z = \frac{1}{2} \rho V^2 S L m_zMz=21ρV2SLmz
式中 mz(Ma,α,δz,ωz)m_z(Ma, \alpha, \delta_z, \omega_z)mz(Ma,α,δz,ωz) 为俯仰力矩系数,LLL 为弹体特征长度(一般取弹长)。

与力矩相关的导数定义为
Mzα=∂Mz∂α=12ρV2SLmzα,mzα=∂mz∂αM_z^\alpha = \frac{\partial M_z}{\partial \alpha} = \frac{1}{2} \rho V^2 S L m_z^\alpha, \quad m_z^\alpha = \frac{\partial m_z}{\partial \alpha}Mzα=αMz=21ρV2SLmzα,mzα=αmz
Mzδz=∂Mz∂δz=12ρV2SLmzδz,mzδz=∂mz∂δzM_z^{\delta_z} = \frac{\partial M_z}{\partial \delta_z} = \frac{1}{2} \rho V^2 S L m_z^{\delta_z}, \quad m_z^{\delta_z} = \frac{\partial m_z}{\partial \delta_z}Mzδz=δzMz=21ρV2SLmzδz,mzδz=δzmz
Mzωz=∂Mz∂ωz=12ρV2SLmzωz,mzωz=∂mz∂ωzM_z^{\omega_z} = \frac{\partial M_z}{\partial \omega_z} = \frac{1}{2} \rho V^2 S L m_z^{\omega_z}, \quad m_z^{\omega_z} = \frac{\partial m_z}{\partial \omega_z}Mzωz=ωzMz=21ρV2SLmzωz,mzωz=ωzmz

3.2 压心位置与静稳定性

MzαM_z^\alphaMzα 是气动力导数 YαY^\alphaYα 与质心到压心距离的乘积,该距离称为静稳定度。本质上它是气动力抗干扰项,用以表示导弹静稳定性的大小。

若压心在质心之后,导弹受扰产生攻角时,气动恢复力矩减小攻角,弹体称为静稳定弹体。反之,若压心在质心之前,扰动产生的攻角将导致更大攻角,弹体称为静不稳定弹体

一般来说,若质心距弹体头部约为弹长的百分之五十,则期望将压心置于距头部约为弹长百分之五十二至五十五的位置。在亚音速和低超音速时,压心位置比高马赫数情况下更靠前;低速飞行时,压心位置受攻角影响显著,这是因为弹体压心随攻角增大常往后移动,而控制面和弹翼的压心变化不大。

3.3 动力系数定义

为方便控制系统设计,通常将气动力导数除以质量和速度的乘积 mVmVmV,气动力矩导数除以各自的转动惯量。这就得到了与控制系统设计相关的各动力系数。

与转动有关的动力系数为 aα,aδ,aωa_\alpha, a_\delta, a_\omegaaα,aδ,aω;与平动有关的动力系数为 bα,bδb_\alpha, b_\deltabα,bδ;滚转动力系数为 cδ,cωc_\delta, c_\omegacδ,cω。各系数的代数符号与量纲见原书表2.4-1,主要表达式见原书表2.4-2。

动力系数的符号规定如下:

  • 对典型的气动静稳定、正常式控制导弹,所有气动动力系数 (aα,aδ,aω,bα,bδ,cδ,cω)(a_\alpha, a_\delta, a_\omega, b_\alpha, b_\delta, c_\delta, c_\omega)(aα,aδ,aω,bα,bδ,cδ,cω) 符号皆为正。
  • 对鸭式控制的导弹,只有 aδ<0a_\delta < 0aδ<0,其他皆为正。
  • 对静不稳定导弹,只有 aα<0a_\alpha < 0aα<0

3.4 各动力系数物理意义

aαa_\alphaaα 的物理意义为单位攻角产生的导弹俯仰旋转角加速度大小,单位为秒的负二次方,它反映了导弹静稳定性的大小。aαa_\alphaaα 越大,静稳定性越强,弹体自振频率越高。

aδa_\deltaaδ 为单位舵转角产生的导弹俯仰旋转角加速度大小,单位为秒的负二次方,它反映了舵对导弹转动控制的效率。

aωa_\omegaaω 为单位角速度产生的导弹俯仰旋转角加速度大小,单位为秒的负一次方,它反映了导弹气动阻尼的大小。

bαb_\alphabα 为单位攻角产生的导弹速度矢量旋转角速度大小,单位为秒的负一次方,它反映了导弹攻角产生的法向力转动速度矢量能力的大小,是表征导弹机动效率的重要气动导数。

bδb_\deltabδ 为单位舵转角产生的导弹速度矢量旋转角速度大小,单位为秒的负一次方。由于其值相对于 bαb_\alphabα 较小,因此对旋转速度矢量的贡献有限,但其大小对 aδa_\deltaaδ 有直接影响。

cδc_\deltacδ 为单位滚转舵转角产生的导弹滚转角加速度大小,单位为秒的负二次方,它反映了导弹滚转舵控制效率的大小。

cωc_\omegacω 为单位导弹滚转角速度产生的导弹滚转角加速度大小,单位为秒的负一次方,它反映了导弹滚转阻尼的大小。


四、协同设计图(数据流)

风洞试验

数据处理

非线性曲线

特征点选取

泰勒展开

气动力导数

除质量惯量

动力系数

控制律设计

驾驶仪设计

从风洞试验到控制设计的完整数据流如下:风洞试验给出不同马赫数、攻角、舵偏角下的气动力与力矩系数曲线;经数据处理后得到非线性函数关系;在特征工作点处进行泰勒展开,保留一阶项得到气动力导数;将导数除以 mVmVmV 或转动惯量,得到无量纲动力系数;最终代入状态方程,形成可直接用于控制律与驾驶仪设计的线性定常模型。


五、接口对接图

外部输入

马赫数

攻角

舵偏角

角速度

侧滑角

气动力计算

动力系数出

a阿尔法

a德尔塔

a欧米伽

b阿尔法

b德尔塔

c德尔塔

c欧米伽

动力系数计算模块的输入接口为导弹飞行状态参数与几何参数:马赫数决定气动环境,攻角与侧滑角决定弹体姿态,舵偏角与角速度决定控制输入。输出接口为七个核心动力系数,分别对应俯仰转动、平动与滚转三个通道。这些系数作为常数代入后续的状态空间方程或传递函数,实现从气动专业到控制专业的无缝衔接。


六、完整三层架构总览图

应用层

信号层

物理层

风洞试验

飞行试验

大气数据

弹体几何

非线性曲线

特征点

气动力导数

动力系数

状态空间

传递函数

频域分析

根轨迹

三层架构清晰划分了气动建模的工程边界。物理层提供试验数据与弹体几何参数,风洞试验给出不同马赫数与攻角下的气动力系数曲线,大气数据提供密度与声速,弹体几何提供参考面积与长度。信号层完成从非线性到线性的抽象,特征点选取将连续弹道离散化,泰勒展开得到气动力导数,除以质量与惯量后得到无量纲动力系数。应用层将动力系数组装为状态空间矩阵或传递函数,进而使用频域法、根轨迹法等成熟工具进行控制系统设计。


七、本节核心结论

核心结论

线性化三步

导数本质

系数桥梁

符号规定

特征点

泰勒展

保留一阶

偏导数

特征常数

除mV

除惯量

静稳全正

鸭式a德尔塔负

静不稳a阿尔法负

本节建立了从非线性气动力到线性控制模型的完整桥梁。核心成果可归纳为四点:

第一,小扰动线性化遵循三步流程:在弹道特征点处将非线性方程泰勒展开,保留一阶小量,假设过渡过程内系数不变,得到线性定常微分方程组。

第二,气动力导数本质上是力与力矩对运动变量的一阶偏导数,如 Yα,Mzα,Mzδz,MzωzY^\alpha, M_z^\alpha, M_z^{\delta_z}, M_z^{\omega_z}Yα,Mzα,Mzδz,Mzωz 等。这些导数在特征点上视为常数,其数值由风洞试验或飞行试验数据处理后得到。

第三,动力系数是导数除以 mVmVmV 或转动惯量的工程化结果,实现了量纲统一与控制理论的对接。七个核心动力系数 aα,aδ,aω,bα,bδ,cδ,cωa_\alpha, a_\delta, a_\omega, b_\alpha, b_\delta, c_\delta, c_\omegaaα,aδ,aω,bα,bδ,cδ,cω 分别表征静稳定性、舵效、阻尼、机动效率、舵面法向贡献、滚转舵效与滚转阻尼。

第四,符号规定直接反映气动布局特性:静稳定正常式导弹所有动力系数为正;鸭式布局仅 aδa_\deltaaδ 为负;静不稳定导弹仅 aαa_\alphaaα 为负。这一符号体系是后续第2.5节传递函数推导与第4章驾驶仪设计时判断系统稳定性的直接依据。

衔接下节:第2.5节将利用本节定义的七个动力系数,代入第2.2节简化的六自由度方程,通过拉普拉斯变换直接导出弹体俯仰通道与滚转通道的传递函数,完成从动力学方程到控制对象的最终转换。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

VectorShift

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值