Чудаков М.В. Построение поляр и расчет динамики полета дозвуковых транспортных и пассажирских самолетов
Подождите немного. Документ загружается.
141
Re
f
2.991 10
8
×=
Re
f
V
ras
L
f
ν
h
⋅:=
∆C
xa_f
0.008:=
Определение лобового сопротивления фюзеляжа и гондол
двигателя
C
xa_oper
0.009:=
S
oper
63.87:=
Определение коэффициента лобового сопротивления
оперения
C
xamin
0.0143=
C
xamin
C
xa_min_izol_kr
1k
int
S
kr_fuz
S
⋅−
⋅:=
C
xa_min_izol_kr
C
xpminizolkr
Σ∆C
xa
+:=
C
xpminizolkr
dvaC
f
0.93 2.8 c
otn
⋅+
()
⋅ 15c
otn
⋅+ M
ras
4
+
⋅:=
Σ∆C
xa
0.0025:=
k
int
0.25:=
dvaC
f
0.005:=
Определяем по графику 1.4 коэффициент 2Сf
Re 3.781 10
7
×=
Re V
ras
b
eqv
ν
h
⋅:=
Определение минимального коэффициента лобового
сопротивления крыла
м
b
eqv
5.379=
b
eqv
S
l
:=
Эквивалентная хорда крыла, м
142
C
x0
0.03=
C
x0
1.03
SC
xamin
⋅ S
mid
C
xa_f
⋅+ S
mid_GD
n
GD
⋅ C
xa_GD
⋅+ S
oper
C
xa_oper
⋅+
S
⋅:=
Площадь миделя гондол двигателя
S
mid_GD
1.58:=
Число гондол двигателя
n
GD
3:=
Минимальный коэффициент суммы сопротивлений самолета
C
xa_GD
0.04:=
C
xa_f
0.058=
C
xa_f
C
f
η
c
⋅η
m
⋅
S
pov
S
mid
⋅∆C
xa_f
+:=
S
pov
766.963=
S
pov
2.85 L
f
⋅ S
mid
⋅:=
η
m
1.2:=
Определяем по графику 1.7 коэффициент
ηµ
η
c
1.08:=
Определяем по графику 1.6 коэффициент
η
c
λ
fn
1.191=
λ
fn
L
fn
S
mid
π
:=
λ
f
5.962=
λ
f
L
f
4
π
S
mid
⋅
:=
C
f
0.002:=
Определяем по графику 1.5 коэффициент Сf
143
Cya
i
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
C
ya_max_kr
:=
i012..:=
Построение поляры самолета
C
ya_max_kr
1.172=
C
ya_max_kr
C
ya_max_seh
K
η
⋅ 1 cos χ
()
+
()
⋅
2
:=
K
η
0.92:=
Определение по табл. 1.3 коэффициента К
η
C
ya_max_seh
1.4:=
Определение прироста коэффициента вредных
сопротивлений при углах атаки, отличных от нулевого
λ
ef
λ
1
S
kr_fuz
S
+
:=
δ 0.03:=
Определение по рис. 1.8 поправки
δ
η 4.994=
η
b
0
b
k
:=
λ 6.98=
λ
l
2
S
:=
Определение индуктивного сорротивления самолета
144
Cxa
i
1 δ+
()
πλ
ef
⋅
Cya
i
()
2
⋅:=
C
ya_sred
i
Cya
i
C
ya_max_kr
:=
C
ya_sred
i
0
0
1
2
3
4
5
6
7
8
9
10
11
12
0
0.085
0.171
0.256
0.341
0.427
0.512
0.598
0.683
0.768
0.854
0.939
1
=
Определение по рис 1.10 коэффициент
∆
С
x
a_vr
∆C
xa_vr
i
0
0
0
0
0
0
0
0.001
0.003
0.004
0.012
0.02
0.033
:=
Определение Cxa_сам
C
xa_cam
i
C
x0
∆C
xa_vr
i
+ Cxa
i
+:=
C
xa_cam
i
0
0
1
2
3
4
5
6
7
8
9
10
11
12
0.03
0.031
0.033
0.035
0.039
0.044
0.05
0.058
0.068
0.078
0.096
0.115
0.137
=
145
0.02 0.032 0.044 0.056 0.068 0.08 0.092 0.1 0.12 0.13 0.14
0.15
0.3
0.45
0.6
0.75
0.9
1.05
1.2
1.35
1.5
Cya
i
C
xa_cam
i
Построение кривой коэффициента подъемной силы Cya=f(
α
)
a
0
0.1:=
a
χ
a
0
cos χ
()
⋅
1
a
0
cos χ
()
⋅
πλ⋅
+
:=
a
χ
0.082=
α
kr
C
ya_max_kr
1M
ras
2
−⋅
a
χ
α
0
+:=
α
kr
9.875=
j01..:=
α
j
0
7
:=
Линейный участок кривой Cya=f(
α
)
C
ya
j
a
χ
α
j
α
0
−
()
⋅
1M
ras
2
−
:=
146
C
ya
j
0.023
0.837
=
Массив шкалы значений
p06..:=
C
ya
p
0.023
0.837
0.9
1
1.1
1.172
1.168
1.6
:=
α
p
0
7
8
9
9.5
α
kr
10
11
:=
012345678910
0.12
0.24
0.36
0.48
0.6
0.72
0.84
0.96
1.08
1.2
C
ya
p
α
p
147
Определение максимальных коэффициентов лобового
сопротивления и подъемной силы крыла при посадке с
применением механизации крыла
Определение максимальных коэффициентов лобового
сопротивления
∆C
xa_mex
0.1:=
l
mex
1.87:=
S
kr_mexS.kr
2 η⋅ l
mex
η 1−
()
⋅−
η 1+
l
mex
l
⋅:=
h 5.64 10
6−
×:=
λ
fik
λ
ef
1 0.065
1
h
⋅+
⋅:=
Cxaimex
i
1 δ+
()
πλ
fik
⋅
Cya
i
()
2
⋅:=
C
xamex
i
1.3 C
x0
⋅∆C
xa_mex
S
kr_mexS.kr
⋅+∆C
xa_vr
i
+ Cxaimex
i
+:=
C
xamex
i
0.042
0.042
0.042
0.042
0.042
0.042
0.042
0.043
0.045
0.046
0.054
0.062
0.075
=
148
Определение максимальных коэффициентов подъемной силы
крыла при посадке с применением механизации крыла
∆C
ya_vid
0.32
c
otn
0.16
+
S
kr_mexS.kr
⋅:=
Cyamaxmex
i
Cya
i
∆C
ya_vid
cos χ
()
⋅+:=
Cyamaxmex
i
0.019
0.119
0.219
0.319
0.419
0.519
0.619
0.719
0.819
0.919
1.019
1.119
1.191
=
Поляра при посадке с применением механизации крыла
0.04 0.0435 0.047 0.0505 0.054 0.0575 0.061 0.0645 0.068 0.0715 0.075
0.12
0.24
0.36
0.48
0.6
0.72
0.84
0.96
1.08
1.2
Cyamaxmex
i
C
xamex
i
149
Построение поляр при М>Мкр
М=0,75
z03..:=
C
xavol_kr_χ
z
0.002
0.004
0.006
0.01
:=
C
xavol_kr
z
C
xavol_kr_χ
z
cos χ
()
3
⋅:=
C
xa_vol_f
0:=
C
xa_vol_op
0.002:=
C
xa_vol_cam
z
C
xavol_kr
z
C
xa_vol_op
S
oper
⋅ C
xa_vol_f
S
mid
⋅+
()
S
+:=
C
xa_cam
z
0.03
0.031
0.033
0.035
:=
C
xacam
z
C
xa_cam
z
C
xa_vol_cam
z
+:=
C
xacam
z
0.032
0.034
0.037
0.041
=
∆C
ya_vol
z
0
0.08
1.12
1.16
:=
150
C
yacam
z
0.023
0.837
0.954
1.07
:=
C
ya_cam
z
C
yacam
z
∆C
ya_vol
z
+:=
C
ya_cam
z
0.023
0.917
2.074
2.23
=
М=0,8
C
xavol_kr08_χ
z
0.008
0.01
0.012
0.019
:=
C
xavol_kr08
z
C
xavol_kr08_χ
z
cos χ
()
3
⋅:=
C
xa_vol_f08
0.1:=
C
xa_vol_cam08
z
C
xavol_kr08
z
C
xa_vol_op
S
oper
⋅ C
xa_vol_f08
S
mid
⋅+
()
S
+:=
C
xa_vol_cam08
z
0.025
0.026
0.027
0.031
=
C
xacam08
z
C
xa_cam
z
C
xa_vol_cam08
z
+:=
C
xacam08
z
0.055
0.057
0.06
0.066
=