PX4固定翼控制器详解(二)——姿态控制器(续)

上期我们讲述了双回路姿态控制器的设计方式,但是留下来了两个问题。

  • 第一个问题是为什么采用协调转弯公式计算偏航角速率期望值 ψ ˙ s p \dot{\psi}_{sp} ψ˙sp
  • 第二个问题是空速缩放为什么 P I PI PI控制项的缩放系数 P I S c a l e r PI_{Scaler} PIScaler是平方比例计算,而 F F FF FF前馈项的缩放系数 F F S c a l e r FF_{Scaler} FFScaler计算是线性比例计算?

1. 为什么需要协调转弯?

侧滑角是机体系下的 x b x_b xb轴朝向与空速方向 V T V_T VT的夹角。
无风环境下,当飞行器沿着路径 C C C进行转弯时,机头朝向与飞行器飞行方向(飞行器在路径 C C C处的切线方向)之间存在的夹角为侧滑角。

当飞行存在侧滑角时,会有以下缺点:

  • 由于空速方向与机体纵轴不重合,会导致额外较大的阻力;
  • 机头朝向与飞行方向不一致,会引起乘坐不适;
  • 侧滑不利于飞机的操控和机动性;
  • 侧滑会导致机翼两侧受力不均匀,长期的侧滑会使飞机承受额外的应力,导致结构疲劳。

为了克服侧滑角的影响,需要实现协调转弯,控制机头朝向与速度方向一致。

我们先来看协调转弯下的受力分析图。升力为 L L L,升力的垂直分量为 L c o s ϕ Lcos\phi Lcosϕ,升力的水平分量为 L s i n ϕ Lsin\phi Lsinϕ。假设在等高度进行协调转弯,不难得出如下方程组:
{ m g = L c o s ϕ m a y = L s i n ϕ a y = V ψ ˙ \begin{cases}mg=Lcos\phi \\ma_y=Lsin\phi \\a_y = V \dot{\psi} \end{cases} mg=Lcosϕmay=Lsinϕay=Vψ˙
式中, a y a_y ay为转弯向心加速度, L L L为升力, V V V为空速。
求解改方程组可得协调转弯公式:
ψ ˙ = g V t a n ϕ \dot{\psi}=\frac{g}{V}tan\phi ψ˙=Vgtanϕ
若不在等高度下(匀速爬升)进行协调转弯,同样的受力分析如下图:

联立方程组:
{ m g = L c o s ϕ m a y = L s i n ϕ   c o s θ a y = V ψ ˙ \begin{cases} mg=Lcos\phi \\ma_y=L sin\phi\ cos\theta \\a_y = V \dot{\psi} \end{cases} mg=Lcosϕmay=Lsinϕ cosθay=Vψ˙
不难发现,不等高度下(匀速爬升)协调转弯,只需额外乘一个俯仰角 θ \theta θ的余弦值,得到协调转弯公式:
ψ ˙ = g V t a n ϕ   c o s θ \dot{\psi}=\frac{g}{V}tan\phi \ cos\theta ψ˙=Vgtanϕ cosθ
由此,我们看外环控制图中

期望的偏航角速度值 ψ ˙ s p \dot{\psi}_{sp} ψ˙sp,根据协调转弯公式可以得到:
ψ ˙ s p = g V T t a n ϕ s p c o s θ s p \dot{\psi}_{sp}=\frac{g}{V_T}tan\phi_{sp}cos\theta_{sp} ψ˙sp=VTgtanϕspcosθsp

2. 如何设计空速缩放?

内环的角速度控制中,空速缩放有两部分设计, P I PI PI空速缩放和 F F FF FF前馈空速缩放。

以滚转通道为例,滚转力矩计算公式在之前的气动力矩建模的文章中有详细解释,这里给出公式如下:
L A = C l ( 1 2 ρ V 2 ) S W b L_A=C_l(\frac{1}{2}\rho V^2)S_Wb LA=Cl(21ρV2)SWb
C l = C l β β + C l δ a δ a + C l δ r δ r + C l p ‾ p ‾ + C l r ‾ r ‾ C_l=C_l^{\beta}\beta+C_l^{\delta_a}\delta_a+C_l^{\delta_r}\delta_r+C_l^{\overline{p}}\overline{p}+C_l^{\overline{r}}\overline{r} Cl=Clββ+Clδaδa+Clδrδr+Clpp+Clrr
式中, p ‾ = p b 2 V \overline{p}=\frac{pb}{2V} p=2Vpb r ‾ = r b 2 V \overline{r}=\frac{rb}{2V} r=2Vrb

2.1. P I PI PI空速缩放——静态力矩补偿

p = 0 p=0 p=0时,阻尼项 C l p ‾ p ‾ C_l^{\overline{p}}\overline{p} Clpp为0,此时舵面偏转产生的瞬时力矩为:
l = q S b C l δ a δ a l=qSbC_l^{\delta_{a}}\delta_{a} l=qSbClδaδa
式中, q = 1 2 ρ V T 2 q=\frac{1}{2}\rho V_T^2 q=21ρVT2为动压。
提取出舵面偏转, k = 2 l ρ S b C l δ a k=\frac{2l}{\rho SbC_l^{\delta_{a}}} k=ρSbClδa2l为常数:
δ a = l q S b C l δ a = 2 l ρ V T 2 S b C l δ a = 2 l ρ S b C l δ a 1 V T 2 = k 1 V T 2 \delta_{a}=\frac{l}{qSbC_l^{\delta_{a}}} =\frac{2l}{\rho V_T^2SbC_l^{\delta_{a}}} =\frac{2l}{\rho SbC_l^{\delta_{a}}}\frac{1}{V_T^2} =k\frac{1}{V_T^2} δa=qSbClδal=ρVT2SbClδa2l=ρSbClδa2lVT21=kVT21
由此不难看出, P I PI PI控制下的舵偏 δ a \delta_a δa 与空速平方 V T 2 V_T^2 VT2 成线性比例关系。若要在不同空速下,产生同样大小的力矩,需要对 δ a \delta_{a} δa 进行对空速平方 V T 2 V_T^2 VT2进行比例缩放。

2.2 F F FF FF空速缩放——角速率补偿

速率控制器前馈的主要用途是补偿自然速率阻尼。假设当下匀速转动,滚转力矩为0,副翼产生的力矩应该精确补偿阻尼项产生的力矩:
− C l δ a δ a = C l p ‾ p ‾ = C l p ‾ p b 2 V T -C_l^{\delta_a}\delta_a=C_l^{\overline{p}}\overline{p}=C_l^{\overline{p}}\frac{pb}{2V_T} Clδaδa=Clpp=Clp2VTpb
提取出舵面偏转:
δ a = − C l p ‾ C l δ a p b 2 V T \delta_a=-\frac{C_l^{\overline{p}}}{C_l^{\delta_a}} \frac{pb}{2V_T} δa=ClδaClp2VTpb
由此可以看出, F F FF FF控制下的舵偏 δ a \delta_a δa V T V_T VT成线性比例关系。因此,该回路的舵偏 δ a \delta_{a} δa进行对空速 V T V_T VT的线性比例缩放。

3. 总结

本期文章介绍了协调转弯控制率,空速缩放设计原理。后续会陆续带来规划、控制相关的文章、仿真分享。

END

迅翼SwiftWing致力于固定翼技术共享,汇聚固定翼领域技术极客,推动固定翼技术持续创新!

  • 5
    点赞
  • 13
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值