Анализ динамики углового движения неуправляемого летательного аппарата методами качественной теории систем - page 5

β
=
=
(1+Δ
λ
1
λ
2
)(
β
0
+
ε
т
K
11
+
d
T
K
12
) +
μ
1
λ
(
ε
т
K
21
+
d
т
K
22
)
(1 + Δ
λ
1
λ
2
)
2
+ (
μ
1
λ
)
2
β
(
λ
);
(4)
α
=
=
(1+Δ
λ
1
λ
2
)(
ε
т
K
21
+
d
т
K
22
)
μ
1
λ
(
β
0
+
ε
т
K
11
+
d
т
K
12
)
(1 + Δ
λ
1
λ
2
)
2
+ (
μ
1
λ
)
2
Φ
1
(
λ
);
(5)
Φ
2
(
λ
)
α
=
ˉ
m
ω
x
ω
кр
ˉ
C
α
n
Δ
z
λ
ˉ
m
x
0
ˉ
C
α
n
Δ
z
ε
т
d
т
ˉ
C
α
n
Δ
z
P
т
J
x
sin(
ϕ
1
ϕ
2
)
,
(6)
где
C
11
=
ˉ
P
т
ω
кр
cos
ϕ
1
;
C
12
=
ˉ
m
ω
z
ˉ
m
α
z
ˉ
P
т
cos
ϕ
1
=
C
22
;
C
13
=
P
т
ˉ
m
α
z
x
т
J
sin
ϕ
1
;
C
21
=
ˉ
P
т
ω
кр
sin
ϕ
1
;
C
23
=
P
т
ˉ
m
α
z
x
т
J
cos
ϕ
1
;
K
12
=
P
т
ˉ
m
α
z
J
sin
ϕ
2
;
K
11
=
C
11
λ
+
C
12
+
C
13
;
K
21
=
C
21
λ
+
C
22
+
C
23
;
K
22
=
P
т
ˉ
m
α
z
J
cos
ϕ
2
.
Особые точки (
λ
,
α
,
β
)
системы уравнений (2) можно опреде-
лить, решая совместно уравнения (4), (5) и (6). Рассмотрим качествен-
ную картину изменения угла атаки в пространстве
α
λ
для случая
наличия заданных возмущений
β
0
. В соответствии с уравнениями (5)
и (6) при
λ
= 0
имеем следующие значения угла атаки:
Φ
1
(
λ
)
λ
=0
=
α
λ
=0
=
ε
т
(
C
22
+
C
23
) +
d
т
K
22
(1 + Δ
λ
1
)
Φ
1
(
λ
)
0
;
(7)
Φ
2
(
λ
)
λ
=0
=
α
λ
=0
=
ˉ
m
x
0
ˉ
C
α
n
Δ
z
ε
т
d
т
P
т
J
x
sin(
ϕ
1
ϕ
2
)
ˉ
C
α
n
Δ
z
Φ
2
(
λ
)
0
.
(8)
Из уравнений (7) и (8) очевидно, что
Φ
1
(
λ
)
0
зависит только от
возмущающих факторов эксцентриситета силы тяги
ε
т
и
d
т
, а функция
Φ
2
(
λ
)
0
зависит от произведения
ε
т
d
т
и возмущающего фактора
Δ
z
.
Угол наклона прямой
Φ
2
(
λ
)
в плоскости
α
λ
с учетом формулы
(6) определяется следующим соотношением:
d
2
(
λ
)] =
ˉ
m
ω
x
ω
кр
ˉ
C
α
n
Δ
z
.
(9)
Значение величины
λ
, при котором прямая
Φ
2
(
λ
)
пересекает ось
абсцисс в плоскости
α
λ
, определяется при
α
= 0
из уравнения (6)
как
λ
0
=
1
ˉ
m
ω
x
ω
кр
ˉ
m
x
0
+
ε
т
d
т
P
т
J
x
sin(
ϕ
1
ϕ
2
)
.
(10)
28 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2006. № 1
1,2,3,4 6,7
Powered by FlippingBook