本章定位:全书的被控对象数学模型。第4章驾驶仪设计、第8章制导律,全部建立在第2章的传递函数之上。掌握本章 = 掌握后续所有设计的根基。
一、第总架构图
原理点睛:本章只做一件事——把刚体弹体翻译成传递函数。
二、核心原理知识图谱:从物理到数学
原理点睛:
- 小扰动线性化:弹体在基准弹道附近微小偏离,非线性方程→线性定常方程。
- 气动力导数:力/力矩对 α, δ, ω 的偏导数,如 Yα=∂Y/∂αY^\alpha = \partial Y/\partial \alphaYα=∂Y/∂α。
- 动力系数:导数除以 mVmVmV 或 JJJ,统一量纲,得到 aα,bαa_\alpha, b_\alphaaα,bα 等工程化参数。
三、坐标系转换关系
原理点睛:
- 弹体系原点:质心;xbx_bxb 沿弹轴向前。
- 速度系:xvx_vxv 沿速度矢量。
- 攻角 α\alphaα:弹轴与速度矢量在俯仰面夹角;侧滑角 β\betaβ:偏航面夹角。
- 总攻角 αT\alpha_TαT 满足:tanα=tanαT⋅cosλ\tan\alpha = \tan\alpha_T \cdot \cos\lambdatanα=tanαT⋅cosλ,tanβ=tanαT⋅sinλ\tan\beta = \tan\alpha_T \cdot \sin\lambdatanβ=tanαT⋅sinλ(λ\lambdaλ 为攻角平面角)。
四、气动力导数 → 动力系数(核心公式链)
原理点睛:
- aα=−Mzα/Jza_\alpha = -M_z^\alpha / J_zaα=−Mzα/Jz:单位攻角产生的俯仰角加速度 → 静稳定性指标(aα>0a_\alpha>0aα>0 为静稳定)。
- bα=(P+Yα)/mVb_\alpha = (P+Y^\alpha)/mVbα=(P+Yα)/mV:单位攻角产生的速度矢量转角速度 → 机动效率指标。
- 静稳定度 x∗=−Mzα/Yα=aαJz/(bαmV)x^* = -M_z^\alpha/Y^\alpha = a_\alpha J_z / (b_\alpha mV)x∗=−Mzα/Yα=aαJz/(bαmV):压心到质心距离,一般取弹长 4%~8%。
五、俯仰通道传递函数全家福(重中之重)
原理点睛(必须死记):
所有以 δz\delta_zδz 为输入的俯仰通道传函,分母完全相同:
s2+(aω+bα)s+(aα+aωbα)s^2 + (a_\omega+b_\alpha)s + (a_\alpha+a_\omega b_\alpha)s2+(aω+bα)s+(aα+aωbα)
这是一个二阶振荡环节,无阻尼自振频率 ωm≈aα\omega_m \approx \sqrt{a_\alpha}ωm≈aα(因 aωbαa_\omega b_\alphaaωbα 极小)。
| 输出 | 分子特征 | 物理意义 |
|---|---|---|
| aya_yay | 二阶零点 A2s2+A1s+1A_2s^2+A_1s+1A2s2+A1s+1 | 法向过载(现代设计不可略去分子项) |
| θ˙\dot\thetaθ˙ | 同 aya_yay(因 ay=Vθ˙a_y=V\dot\thetaay=Vθ˙) | 弹道倾角速度 |
| ϑ˙\dot\varthetaϑ˙ | 一阶 (Tαs+1)(T_\alpha s+1)(Tαs+1) | 俯仰角速度,Tα=1/bαT_\alpha=1/b_\alphaTα=1/bα 叫攻角时间常数 |
| α\alphaα | 一阶 (B1s+1)(B_1s+1)(B1s+1) | 攻角,稳态增益 kα=aδ/aαk_\alpha=a_\delta/a_\alphakα=aδ/aα 叫操控比 |
关键结论:ϑ˙(s)=(Tαs+1)Vay(s)\dot\vartheta(s) = \frac{(T_\alpha s+1)}{V} a_y(s)ϑ˙(s)=V(Tαs+1)ay(s),说明角速度反馈相当于对过载的超前校正。
六、静稳定度 x∗x^*x∗ 对弹体特性的影响
原理点睛:
- x∗x^*x∗ 大 → 压心远离质心 → 恢复力矩大 → 自振频率高,但舵效相对变弱(增益小)。
- x∗x^*x∗ 过小 → 操控比剧增,但质心/压心微小偏移会引起增益剧烈变化 → 鲁棒性极差。
- 低静稳弹体必须引入伪攻角反馈(第4章)来人工增稳。
七、滚转通道(一阶惯性)
原理点睛:滚转通道远比俯仰/偏航简单,为一阶惯性:
γ˙(s)δx(s)=krTrs+1,kr=−cδcω, Tr=1cω\frac{\dot\gamma(s)}{\delta_x(s)} = \frac{k_r}{T_r s+1},\quad k_r=-\frac{c_\delta}{c_\omega},\ T_r=\frac{1}{c_\omega}δx(s)γ˙(s)=Trs+1kr,kr=−cωcδ, Tr=cω1
精简版本
第2章 弹道模型、气动导数、动力系数和弹体传函
本章定位:全书的被控对象数学模型。第4章驾驶仪设计、第8章制导律,全部建立在本章的传递函数之上。掌握本章 = 掌握后续所有设计的根基。
本章核心任务:把刚体弹体翻译成传递函数。
一、知识图谱
本章遵循严格的工程数学逻辑:从刚体假设出发,建立坐标系与六自由度方程,通过控制面分解与气动力线性化,最终导出俯仰与滚转两个核心通道的传递函数。所有后续控制设计均建立在最后这两个输出模型之上。
二、总体结构图
总体流程清晰呈现从物理对象到数学模型的转化链路。刚体弹体首先被置于统一的坐标系中,通过欧拉方程描述其空间运动;控制面规定将四舵输入分解为三通道控制量;气动力线性化得到无量纲动力系数;最终组装为俯仰通道的二阶振荡传函与滚转通道的一阶惯性传函。
三、各模块设计图
3.1 符号与坐标系模块
四个坐标系构成所有推导的基准框架。地面系为惯性参考,原点固定于地面。弹体系原点位于质心,与弹体固连,用于描述气动力与力矩。弹道系原点位于质心,轴沿速度矢量,轴位于铅垂面内,用于描述轨迹。速度系与弹道系区别仅在于轴位于弹体纵向对称面内而非铅垂面内。
攻角为弹轴与速度矢量在俯仰面夹角,侧滑角为偏航面夹角。总攻角满足
tanα=tanαT⋅cosλ,tanβ=tanαT⋅sinλ\tan\alpha = \tan\alpha_T \cdot \cos\lambda, \quad \tan\beta = \tan\alpha_T \cdot \sin\lambdatanα=tanαT⋅cosλ,tanβ=tanαT⋅sinλ
式中 λ\lambdaλ 为攻角平面角。
3.2 欧拉方程与简化模块
六自由度模型包含六个动力学方程与六个运动学方程。动力学方程描述状态量变化率受外力驱动,运动学方程描述位置与姿态的积分关系。
在弹道坐标系下,平动动力学方程为
mV˙=Fxt,mVθ˙=Fyt,−mVcosθψ˙V=Fztm\dot{V} = F_{x_t}, \quad mV\dot{\theta} = F_{y_t}, \quad -mV\cos\theta \dot{\psi}_V = F_{z_t}mV˙=Fxt,mVθ˙=Fyt,−mVcosθψ˙V=Fzt
在弹体坐标系下,经三条物理假设简化后,转动方程解耦为
Jxω˙x=Mx,Jyω˙y=My,Jzω˙z=MzJ_x\dot{\omega}_x = M_x, \quad J_y\dot{\omega}_y = M_y, \quad J_z\dot{\omega}_z = M_zJxω˙x=Mx,Jyω˙y=My,Jzω˙z=Mz
三条假设为:轴对称使惯积为零,小角速度使非线性乘积项可略,小侧向速度使哥氏耦合项可略。三轴解耦是后续三通道独立设计的数学根基。
3.3 控制面分解模块
四舵独立偏转角通过线性组合分解为三通道等效舵偏角。滚转取四舵平均,俯仰取上下舵面差值的一半,偏航取左右舵面差值的一半。
δx=14(δ1+δ2+δ3+δ4),δz=12(δ1−δ3),δy=12(δ4−δ2)\delta_x = \frac{1}{4}(\delta_1+\delta_2+\delta_3+\delta_4), \quad \delta_z = \frac{1}{2}(\delta_1-\delta_3), \quad \delta_y = \frac{1}{2}(\delta_4-\delta_2)δx=41(δ1+δ2+δ3+δ4),δz=21(δ1−δ3),δy=21(δ4−δ2)
沿坐标轴正向观察,舵面顺时针偏转为正。正常式布局下,正舵偏产生负方向俯仰与偏航力矩;鸭式布局因舵面位于质心之前,同向舵偏产生的力矩与正常式相反,但力方向保持一致。
3.4 气动力导数与动力系数模块
气动力导数是力与力矩对运动变量的一阶偏导数。升力与俯仰力矩的基本表达式为
Y=12ρV2SCy,Mz=12ρV2SLmzY = \frac{1}{2}\rho V^2 S C_y, \quad M_z = \frac{1}{2}\rho V^2 S L m_zY=21ρV2SCy,Mz=21ρV2SLmz
定义升力导数 Yα,YδzY^\alpha, Y^{\delta_z}Yα,Yδz 与力矩导数 Mzα,Mzδz,MzωzM_z^\alpha, M_z^{\delta_z}, M_z^{\omega_z}Mzα,Mzδz,Mzωz 后,将导数除以质量速度乘积或转动惯量,得到动力系数:
aα=−MzαJz,aδ=−MzδzJz,aω=−MzωzJza_\alpha = -\frac{M_z^\alpha}{J_z}, \quad a_\delta = -\frac{M_z^{\delta_z}}{J_z}, \quad a_\omega = -\frac{M_z^{\omega_z}}{J_z}aα=−JzMzα,aδ=−JzMzδz,aω=−JzMzωz
bα=P+YαmV,bδ=YδzmVb_\alpha = \frac{P+Y^\alpha}{mV}, \quad b_\delta = \frac{Y^{\delta_z}}{mV}bα=mVP+Yα,bδ=mVYδz
符号规定直接反映布局特性:静稳定正常式导弹所有系数为正;鸭式仅 aδa_\deltaaδ 为负;静不稳定导弹仅 aαa_\alphaaα 为负。
3.5 弹体传函模块
俯仰通道以舵偏角 δz\delta_zδz 为输入,法向过载 aya_yay 为输出,其传递函数为二阶振荡环节:
ay(s)δz(s)=ka(A2s2+A1s+1)Tm2s2+2μmTms+1\frac{a_y(s)}{\delta_z(s)} = \frac{k_a(A_2s^2+A_1s+1)}{T_m^2s^2+2\mu_mT_ms+1}δz(s)ay(s)=Tm2s2+2μmTms+1ka(A2s2+A1s+1)
分母决定系统固有特性,无阻尼自振频率 ωm≈aα\omega_m \approx \sqrt{a_\alpha}ωm≈aα。分子为二阶零点,现代导弹设计中绝对不可略去。攻角时间常数 Tα=1/bαT_\alpha = 1/b_\alphaTα=1/bα,表征导弹机动能力。
滚转通道为一阶惯性:
γ˙(s)δx(s)=krTrs+1\frac{\dot{\gamma}(s)}{\delta_x(s)} = \frac{k_r}{T_rs+1}δx(s)γ˙(s)=Trs+1kr
四、协同设计图
从风洞试验到控制设计的完整数据流如下:风洞试验给出不同马赫数、攻角、舵偏角下的气动力与力矩系数曲线;在特征工作点处进行泰勒展开,保留一阶项得到气动力导数;将导数除以 mVmVmV 或转动惯量,得到无量纲动力系数;最终代入状态方程,形成可直接用于控制律与驾驶仪设计的线性定常模型。
五、接口对接图
动力系数计算模块的输入接口为导弹飞行状态参数:马赫数决定气动环境,攻角决定弹体姿态,舵偏角决定控制输入。输出接口为核心动力系数,分别组装为俯仰通道的二阶振荡传函与滚转通道的一阶惯性传函,实现从气动专业到控制专业的无缝衔接。
六、完整三层架构总览图
三层架构清晰划分了建模的工程边界。物理层提供试验数据与弹体几何参数;信号层完成从非线性到线性的抽象,特征点选取将连续弹道离散化,泰勒展开得到气动力导数,除以质量与惯量后得到无量纲动力系数;应用层将动力系数组装为状态空间矩阵或传递函数,进而使用频域法、根轨迹法等成熟工具进行控制系统设计。
七、核心结论
本章建立了从刚体弹体到控制对象的完整数学桥梁,核心成果可归纳为四点:
第一,三轴转动完全解耦。基于轴对称、小角速度、小侧向速度三条假设,滚转、俯仰、偏航三通道的转动方程相互独立,为后续三通道独立驾驶仪设计奠定数学根基。
第二,气动力导数是连接风洞试验与控制理论的桥梁。在特征点处将非线性气动力泰勒展开,保留一阶项得到导数;导数除以 mVmVmV 或转动惯量得到动力系数,实现量纲统一与线性定常模型的构建。
第三,俯仰通道传递函数为二阶振荡环节,其无阻尼自振频率 ωm≈aα\omega_m \approx \sqrt{a_\alpha}ωm≈aα 直接由静稳定性决定;滚转通道为一阶惯性。两个传函构成了后续所有控制设计的被控对象。
第四,静稳定度 x∗x^*x∗ 是控制系统设计的关键权衡参数。一般取弹长百分之四到百分之八,过高导致增益小、阻尼弱,过低导致鲁棒性极差。低静稳弹体必须引入伪攻角反馈人工增稳。
八、掌握自检清单
- 能画出四个坐标系的转换关系,明确弹道系与速度系的轴区别
- 知道五个俯仰动力系数 aα,aδ,aω,bα,bδa_\alpha, a_\delta, a_\omega, b_\alpha, b_\deltaaα,aδ,aω,bα,bδ 的物理意义和量纲
- 能写出俯仰通道二阶振荡分母,并解释 ωm≈aα\omega_m \approx \sqrt{a_\alpha}ωm≈aα
- 知道法向过载传函的分子是二阶零点,且现代设计不可略去
- 理解攻角时间常数 Tα=1/bαT_\alpha=1/b_\alphaTα=1/bα 的物理意义
- 理解静稳定度 x∗x^*x∗ 对增益、频率、阻尼的三重影响
- 能写出四舵到三通道的分解公式
1064



被折叠的 条评论
为什么被折叠?



