Поэтому в первом приближении можно принять, что функция Теодорсена является амплитудно-фазовой характеристикой аэродинамических сил и моментов самолета, колеблющегося в несжимаемом потоке, а функция Сирса - амплитудно-фазовой характеристикой аэродинамических сил и моментов самолета, обтекаемого колеблющимся потоком.
Сложная зависимость функций Теодорсена и Сирса от приведенной частоты не позволяет использовать их непосредственно при исследовании возмущенного движения самолета в турбулентной атмосфере, где имеют место оба вида нестационарного обтекания.
Поэтому для целей динамики полета функции Теодорсена и Сирса необходимо аппроксимировать более простыми соотношениями, чем равенства (1.81) и (1.91).
Достаточно простые аппроксимирующие выражения для нестационарных аэродинамических сил и моментов можно получить, если исходить не из функций Теодорсена я Сирса, а из переходных функций нестационарных аэродинамических сил и моментов, а затем воспользоваться известными соотношениями, связывающими амплитудно-частотные характеристики и переходные функции [6].
Действительно [6], амплитудно-фазовая характеристика может быть получена из передаточной функции путем формальной замены параметра преобразования Лапласа р на j , где j 
1 и -угловая частота колебаний.
В свою очередь, между передаточной и переходной функциями существуют следующие соотношения:
L h t |
W p |
, |
(1.95) |
p |
|
|
|
|
|
|
где L - символ преобразования Лапласа; h(t ) - переходная функция;
W(р) - передаточная функция;
р - параметр преобразования Лапласа.
31
Переходные функции нестационарных аэродинамических сил и моментов являются решениями уравнений (1.50, 1.51), соответствующими скачкообразному изменению угла атаки и вертикальной скорости потока.
Как показано в [3], в случае несжимаемого потока переходные функции нестационарной подъемной силы, обусловленные скачкообразным изменением угла атаки и вертикальной скорости потока, приближенно можно представить в следующем, виде;
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
. |
|
У |
к |
t |
У |
ст |
0,64 |
0,36 1 |
Tk |
|||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
1 |
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
||
У |
т |
t |
У |
ст |
0,5 1 |
e T1T |
0,5 1 e T2T |
|||||
|
|
|
|
|
|
|
|
|
|
|
||
(1.96)
, (1.97)
где Ук - переходная функция нестационарной подъемной силы,
соответствующей скачкообразному изменению угла атаки; Уст - значение подъемной силы, определенное в
соответствии с гипотезой стационарности; |
|
|
|
|
|||||||
У т - |
переходная |
функция |
нестационарной |
подъемной |
|||||||
силы, |
соответствующая |
скачкообразному |
изменению |
||||||||
вертикальной скорости потока; |
|
|
|
|
|
||||||
Tк , T1T , T2T - постоянные времени. |
|
|
|
|
|||||||
|
Значения постоянных времени, входящих в (1.96) и |
||||||||||
.(1.9.7), определяются следующими соотношениями: |
|||||||||||
T |
1,32 |
вА |
T |
3,85 |
вА |
T |
0,5 |
вА |
(1.98) |
||
|
|
|
|
|
|
||||||
k |
|
V |
1T |
|
V |
2T |
|
V |
|
||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
32 |
|
|
|
|
|
|
где вА - средняя аэродинамическая хорда крыла;
V – скорость невозмущенного потока.
С учетом (1.17) и (1.18) равенства (3,96) и (1.97) –
можно представить в виде
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
t . |
|
(1.99) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
У |
к |
|
t |
|
У |
0,64 |
0,36 1 |
eTk |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
t |
|
|
|
|
|
У |
T |
t |
|
У |
|
0,5 1 |
e T1T |
0,5 1 |
e |
T2T |
T |
t |
, |
(1.100) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где У |
|
C |
y |
( |
)Skq ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(t) - угол атаки самолета; |
|
|
|
|
|
|
|
|
|
|
||||||||||
T (t) |
-приращение угла |
атаки |
|
самолета, |
вызванное |
|||||||||||||||
изменением вертикальной скорости потока. |
|
|
|
|
||||||||||||||||
|
Подставляя (1.99) и (1.100) в (1.95), после |
|||||||||||||||||||
преобразования |
получим |
следующие |
|
|
приближенные |
|||||||||||||||
выражения для передаточных функций нестационарной подъемной силы:
У к |
Р |
У |
|
0,64Т к Р |
1 . |
|
|
(1.101) |
||
|
Р |
|
Т к Р 1 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||
У т |
Р |
|
У |
|
0,5 Т1Т |
|
Т 2Т |
Р 1 |
|
(1.102) |
т Р |
Т1Т Т 2Т Р 2 |
Т1Т |
Т 2Т Р 1 |
|
||||||
|
|
|||||||||
33
Передаточные функции (1.101) и (1.102), строго говоря, справедливы только для крыльев больших удлинений, обтекаемых несжимаемым потоком.
Однако, их можно использовать и для самолетов с крыльями любой формы в плане, обтекаемыми как несжимаемым, так и сжимаемым потоком, если постоянные времени определять исходя из переходных функций, полученных для крыльев, имеющих определенную форму в плане.
Численные методы определения переходных функций крыльев различных форм в плане рассмотрены, например, в [3].
Исходя из сказанного, передаточные функции (1.101) и (1.102) целесообразно представить в следующем более общем виде:
|
|
|
|
Ук |
Р |
У |
|
в1к Р |
|
|
|
1 . |
|
(1.103) |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
Р |
|
|
1к Р |
1 |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
Ут |
Р |
|
У |
|
|
|
в1т Р |
1 |
|
|
(1.104) |
|||||||||||
|
|
т |
Р |
|
|
1т Р2 |
|
|
|
2т Р |
1 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
в1к |
квк |
|
|
|
|
|
к1к |
вА |
|
|
в1т |
|
квт |
2т . |
(1.105) |
|||||||||
1к |
|
1к |
V |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
к1к |
вА2 |
|
|
|
|
к2т |
вА |
|
|
|
|
|
|
|
|
||||||||
1т |
V2 |
|
|
2т |
|
V |
|
|
|
|
|
|
|
|
|
|||||||||
Коэффициенты К |
вк |
, |
|
К |
1к |
, К |
вт |
и К |
1т |
, |
входящие в |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
соотношения (1.105), подлежат определению по переходным функциям.
Из (2.99) следует, что для крыла бесконечного размаха, обтекаемого несжимаемым потоком,
34
Квк |
0,64 |
К1к |
1,32 Квт |
0,5 |
. |
К1т |
1,93 |
К2т |
4,35 |
|
|
|
|
Передаточным функциям (1.103) и (1.104) соответствуют следующие линейные дифференциальные уравнения, связывающие во временной области изменения нестационарной подъемной силы с изменением угла атаки, вызванного колебаниями самолета и колебаниями скорости потока:
1к |
У |
к |
t |
У |
к |
t |
У |
в |
|
t |
t . |
|
(1.I06) |
|
|
|
|
|
1к |
|
|
|
|
||||
|
|
t |
|
|
|
Ут |
|
У в1т t |
|
т t |
(1.107) |
||
1т Ук |
2т У t |
|
|
|
|||||||||
Так как запаздывание в изменении подъемной и боковой аэродинамических сил по сравнению с изменением углов атаки и скольжения определяется одними и теми же факторами, то очевидно, что изменение нестационарной боковой аэродинамической силы в зависимости от изменения угла скольжения может быть приближенно описано следующими уравнениями:
|
|
|
|
Zk t |
Z в1к |
|
t |
|
|
t |
. |
(1.108) |
|
|
|
1к Zk t |
|
|
|
|
|
|
|||||
1т Zk |
t |
2т Zk t |
Zт Z |
в1т |
|
t |
|
t , (1.109) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где Z |
K |
- |
нестационарная |
боковая |
|
аэродинамическая |
сила, |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
вызванная колебаниями самолета; |
|
|
|
|
|
|
|
|
|||||
ZR |
- |
нестационарная |
боковая |
аэродинамическая |
сила, |
||||||||
вызванная колебаниями боковой скорости потока,
35