上期我们讲述了双回路姿态控制器的设计方式,但是留下来了两个问题。
- 第一个问题是为什么采用协调转弯公式计算偏航角速率期望值 ψ ˙ 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致力于固定翼技术共享,汇聚固定翼领域技术极客,推动固定翼技术持续创新!