Способ определения статических и нестационарных аэродинамических производных моделей летательных аппаратов и устройство для его осуществления

Способ определения статических и нестационарных аэродинамических производных моделей летательных аппаратов и устройство для его осуществления

Иллюстрации

Показать все

Реферат

Изобретение относится к области экспериментальной аэродинамики летательных аппаратов.

Важность определения статических и нестационарных аэродинамических производных летательных аппаратов в широком диапазоне углов атаки и скольжения обусловлена, прежде всего, необходимостью обеспечения безопасности эксплуатации самолета, в частности предотвращением возникновения сваливания самолета и попадания его в штопор. По мировой статистике, более половины аварий и катастроф происходит именно по этой причине (Авиация общего назначения. Рекомендации для конструкторов. Под ред. д.т.н. проф. В.Г. Микеладзе, М., ЦАГИ, 2001, с.213).

Из существующего уровня техники известен способ определения статических и нестационарных аэродинамических производных моделей летательных аппаратов в аэродинамической трубе методом вынужденных колебаний, при котором модель посредством тензовесов и державки устанавливают в подвижное поддерживающее устройство, модель колеблют с фиксированной частотой и малой амплитудой поочередно относительно связанных осей Ox, Oy и Oz, измеряют во времени действующие на модель аэродинамические нагрузки и ее угловое положение (Г.С. Бюшгенс, Р.В. Студнев, Аэродинамика самолета. Динамика продольного и бокового движения. М., Машиностроение, 1979, с.32-34). Недостатком данного способа является невозможность определения продольных статических и нестационарных аэродинамических производных моделей летательных аппаратов и характеристик рыскания при углах скольжения, не равных нулю, в случае, если угол атаки модели не равен нулю, а также характеристик крена при углах скольжения, не равных нулю, во всем диапазоне углов атаки модели.

Наиболее близким к заявленному техническому решению является способ определения статических и нестационарных аэродинамических производных моделей летательных аппаратов, при котором испытания проводят как в потоке, так и без потока в аэродинамической трубе, модель посредством тензовесов и державки устанавливают в подвижное поддерживающее устройство, модель колеблют с фиксированной частотой и малой амплитудой поочередно относительно осей Ox, Oy и Oz, измеряют во времени действующие на модель нагрузки и ее угловое положение, производят вычитание инерционных и весовых нагрузок, а затем производят вычисления статических и комплексов нестационарных и вращательных аэродинамических производных (Экспериментальные исследования и математическое моделирование нестационарных аэродинамических характеристик моделей самолетов. Труды ЦАГИ - сборник статей. М., 2010 г., выпуск 2689, с.5-6). Недостатком указанного способа является то, что он также не позволяет определить статические и нестационарные аэродинамические производные моделей летательных аппаратов при углах скольжения, отличных от нуля.

Из существующего уровня техники также известно устройство для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов, содержащее поворотный круг в рабочей части аэродинамической трубы, стойку на поворотном круге, державку модели с тензовесами, смонтированную в подшипниковом узле стойки, привод, связанный тягой с державкой, датчик угла поворота державки (Г.С. Бюшгенс, Р.В. Студнев, Аэродинамика самолета. Динамика продольного и бокового движения. М., Машиностроение, 1979, с.32-34). Недостатком описанного устройства является невозможность проведения испытаний при колебаниях модели в продольном канале и по рысканию при углах скольжения, не равных нулю, в случае, если угол атаки модели не равен нулю, а также отсутствие возможности проведения испытаний при углах скольжения, не равных нулю, во всем диапазоне углов атаки модели при ее колебаниях по крену.

Наиболее близким к заявленному техническому решению является устройство для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов, содержащее поворотный круг в рабочей части аэродинамической трубы, стойку на поворотном круге, подвижную Г-образную раму, державку модели с тензовесами, смонтированную в подшипниковом узле рамы, привод, связанный тягами с рамой и державкой, датчики угла поворота рамы и державки (Экспериментальные исследования и математическое моделирование нестационарных аэродинамических характеристик моделей самолетов. Труды ЦАГИ - сборник статей. М., 2010 г., выпуск 2689, с.5-6). Недостатком этого устройства является то, что оно также не позволяет проводить испытания по определению статических и нестационарных аэродинамических производных моделей летательных аппаратов при углах скольжения, отличных от нуля.

Задачей и техническим результатом изобретения является создание способа и устройства для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов при наличии угла скольжения во всем диапазоне углов атаки, что увеличивает безопасность эксплуатации самолета, в частности, позволяет предотвратить возникновение сваливания самолета и попадание его в штопор.

Решение поставленной задачи и технический результат достигаются тем, что в способе определения статических и нестационарных аэродинамических производных моделей летательных аппаратов, заключающемся в том, что испытания проводят как в потоке, так и без потока в аэродинамической трубе, модель посредством тензовесов и державки устанавливают в подвижное поддерживающее устройство, модель колеблют с фиксированной частотой и малой амплитудой поочередно относительно осей Ox, Oy и Oz, измеряют во времени действующие на модель нагрузки и ее угловое положение, производят вычитание инерционных и весовых нагрузок, в процессе испытаний модель вместе с державкой и тензовесами поворачивают на фиксированные углы крена, результаты эксперимента при колебаниях относительно осей Oy и Oz обрабатывают совместно, при этом значения статических и нестационарных аэродинамических производных определяют по формулам

C α = C 1 ω z cos γ − C 1 ω y sin γ ,

C β = C 1 ω z sin s γ − C 1 ω y cos γ ,

C ω ¯ y   +   C α ˙ ¯   =   C 2 ω z   cos γ − l 2 b a C 2 ω y   sin γ ,

C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β = 2 b a l C 2 ω z sin γ + C 2 ω y cos γ ,

где C=C_y, C_x, m_x, m_y или m_z, C^α - статическая производная аэродинамического коэффициента от угла атаки, C^β - статическая производная аэродинамического коэффициента от угла скольжения, C ω ¯ y + C α ˙ ¯ - комплекс продольных вращательной и нестационарной аэродинамических производных; C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β - комплекс боковых вращательной и нестационарных аэродинамических производных, γ - угол поворота державки по крену, C 1 ω z - составляющая сигнала в фазе с углом колебаний модели при колебаниях относительно оси Oz, C 1 ω y - составляющая сигнала в фазе с углом колебаний модели при колебаниях относительно оси Oy, C 2 ω z - составляющая сигнала в фазе с угловой скоростью колебаний модели при колебаниях относительно оси Oz, C 2 ω y - составляющая сигнала в фазе с угловой скоростью колебаний модели при колебаниях относительно оси Oy, l - размах крыла модели, b_a - САХ модели.

Решение поставленной задачи и технический результат достигаются также тем, что в устройство для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов, содержащее поворотный круг в рабочей части аэродинамической трубы, стойку на поворотном круге, подвижную Г-образную раму, державку модели с тензовесами, смонтированную в подшипниковом узле рамы, привод, связанный тягами с рамой и державкой, датчики угла поворота рамы и державки, введен механизм дистанционной установки державки, состоящий из корпуса, установленного на подшипниковом узле рамы посредством подшипников, электрического мотор-редуктора и выходного вала с упругим диском-шестерней, установленных соосно державке и соединенных с ее свободным концом, пневмотормоза диска, установленного в корпусе и охватывающего упругий диск-шестерню, дополнительного датчика угла установки державки по крену, который вынесен от оси державки и связан с диском-шестерней, рычага-зажима, установленного на наружной поверхности корпуса механизма с возможностью фиксации корпуса относительно рамы и соединения тягой с приводом.

Сущность предложенных способа и устройства для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов поясняется на фиг.1-6, на которых представлены:

на фиг.1 - схема проведения испытаний для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов при колебаниях относительно различных связанных осей;

на фиг.2 - схема заявляемого устройства для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов;

на фиг.3 - схема механизма для дистанционной установки державки под заданным углом крена;

на фиг.4 - график зависимостей углов атаки α и скольжения β от угла поворота круга в рабочей части аэродинамической трубы θ и угла крена державки γ;

на фиг.5, 6, 7 - графики зависимостей статических и нестационарных аэродинамических производных модели треугольного крыла от угла атаки для различных значений угла скольжения.

Перечень используемых позиций в описании изобретения.

1. Поворотный круг в рабочей части аэродинамической трубы.

2. Подвижная Г-образная рама.

3. Стойка.

4. Тяга.

5. Привод.

6. Державка.

7. Тензовесы.

8. Механизм дистанционной установки державки.

9. Модель летательного аппарата.

10. Датчик угла поворота рамы.

11. Подшипниковый узел рамы.

12. Выходной вал.

13. Электрический мотор-редуктор.

14. Упругий диск-шестерня.

15. Датчик угла установки державки по крену.

16. Корпус.

17. Пневмотормоз диска.

18. Подшипник.

19. Рычаг-зажим.

20. Датчик угла поворота державки.

Способ определения статических и нестационарных аэродинамических производных моделей летательных аппаратов заключается в том, что испытания проводят как в потоке, так и без потока в аэродинамической трубе, модель посредством тензовесов и державки устанавливают в подвижное поддерживающее устройство, модель колеблют с фиксированной частотой и малой амплитудой поочередно относительно осей Ox, Oy и Oz, измеряют во времени действующие на модель нагрузки и ее угловое положение, производят вычитание инерционных и весовых нагрузок. Для этого устройство последовательно собирают в конфигурации для проведения исследований производных модели при ее колебаниях относительно одной из связанных осей (фиг.1). Модель летательного аппарата 9 устанавливают в подвижное поддерживающее устройство - державку 8 с тензовесами 7 (фиг.2). При помощи механизма дистанционной установки державки 8 модель поворачивают на фиксированный угол крена γ_i, при помощи поворотного круга в рабочей части аэродинамической трубы 1 модель 9 устанавливают под углом θ_i к оси потока аэродинамической трубы (фиг.2). Привод 5 настраивают на нужные величины амплитуды А и круговой частоты ω и колеблют подвижную Г-образную раму 2 по гармоническому закону

Δ θ ( t ) = A sin ω t .             ( 1 )

Величина амплитуды колебаний является малой и равна 2÷3°. Без потока в аэродинамической трубе в течение шестнадцати периодов колебаний производят измерение тензовесами 7 и запись пяти компонент нагрузок, действующих на модель Y, Z, M_x, M_y, M_z, а также показаний датчиков угла поворота державки 20 и рамы 10. Исследования проводят для нескольких значений углов θ_i. Определяют значение инерционных и весовых нагрузок. Затем включают поток в аэродинамической трубе и указанные операции повторяют. Из нагрузок, полученных в потоке, вычитают инерционные и весовые нагрузки. Затем получают зависимости статических и нестационарных аэродинамических производных модели с использованием алгоритма, описанного ниже.

Для проведения испытаний при заданных средних значениях угла атаки α_i и скольжения β_i модель необходимо устанавливать под соответствующими углами γ_i и θ_i. Связь между величинами указанных углов описывается системой уравнений

t g α = t g θ cos γ sin β = sin θ sin γ                   ( 2 )

При решении этой системы уравнений получают требуемые значения углов поворота круга θ_i и державки γ_i для заданных значений углов атаки и скольжения по следующим формулам

cos 2 θ = cos 2 α + cos 2 β               ( 3 ) sin 2 γ = 1 − cos 2 β 1 − cos 2 α cos 2 β

Графически приведенные зависимости представлены на фиг.4. При проведении исследований необходимые значения углов θ и γ по заданным углам α и β рассчитывают в реальном времени на компьютере.

При повороте модели с тензовесами по крену при колебаниях по тангажу и рысканию угловая скорость вращения модели имеет проекцию одновременно на две оси связанной системы координат модели - ω_z и ω_y (фиг.2). В связи с этим, результаты эксперимента при колебаниях относительно осей Oy и Oz обрабатывают совместно - только одновременная обработка результатов колебаний по тангажу и рысканию с последующим решением соответствующей линейной системы уравнений позволяет получить требуемые зависимости статических и нестационарных аэродинамических производных модели.

Рассмотрим малые колебания модели по тангажу. Пусть угол поворота рамы установки изменяется по закону

Δ θ = ∈ sin ω t ,                     ( 4 )

Где ∈₀<<1 - амплитуда колебаний, ω=2πf - частота. Тогда малые изменения углов атаки и скольжения модели относительно средних значений можно представить формулами

Δ α ( t ) = Δ θ ( t ) cos γ , Δ β ( t ) = Δ θ ( t ) sin γ .                 ( 5 )

Компоненты угловой скорости движения модели можно выразить следующим образом:

ω y = Ω ( t ) sin γ ω z = Ω ( t ) cos γ ,                     ( 6 )

где Ω = Δ θ = ∈ 0 ω cos ω t .

Следовательно, при наличии поворота модели по углу крена при колебаниях по тангажу в связанной системе координат возникают угловые скорости вращения модели ω_y и ω_z, в отличие от случая γ=0, когда ω_y=0.

Выражения для производных по времени углов атаки и скольжения можно найти с помощью кинематических соотношений:

α ˙ = ω z − ( ω x cos α − ω y sin α ) t g β , β ˙ = ω x sin α + ω y cos α .                         ( 7 )

С учетом (6), можно получить

α ˙ = Ω ( t ) ( cos γ + sin γ sin α t g β ) , β ˙ = Ω ( t ) sin γ cos α .                             ( 8 )

Следовательно, в линейном приближении, с учетом (5), (6) и (8), при колебаниях по тангажу аэродинамические коэффициенты можно представить в следующем виде

C = C 0 + C α Δ θ cos γ + C β Δ θ sin γ + C α ˙ ¯ b a V Ω ( sin α t g β sin γ + cos γ ) + + C β ˙ ¯ l 2 V Ω sin γ cos α + C ω ¯ y l 2 V Ω sin γ + C ω ¯ z b a V Ω cos γ = = C 0 ( C α cos γ + C β sin γ ) Δ θ + [ ( C ω ¯ z + C α ˙ ¯ ) ] b a V cos γ + ( C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β ) l 2 V sin γ ] Ω ,

где C=C_y, C_z, m_x, m_y или m_z.

Таким образом, каждый канал тензовесов при малых колебаниях по тангажу воспринимает некоторое среднее значение C₀ и сигналы в фазе с опорным сигналом Δθ(t) и в фазе с сигналом Ω(t)=Δθ(t). Составляющая сигнала в фазе с опорным сигналом пропорциональна

C 1 ω z = C α cos γ + C β sin γ .             ( 9 )

Составляющая сигнала в фазе с угловой скоростью (опережение по фазе на π/2 выражается следующим образом

C 2 ω z = ( C ω ¯ z + C α ˙ ¯ ) b a V cos γ + ( C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β ) l 2 V sin γ .     ( 10 )

Таким образом, при γ ≠ 0 в результате эксперимента при вынужденных колебаниях по тангажу могут быть выделены сложные комплексы производных, описываемые выражениями (9) и (10).

Проанализируем теперь кинематику модели при колебаниях по рысканию. В этом случае рама установки повернута так, что бы ее колебания изменяли угол рыскания модели. При этом ее малые колебания (4) приводят к следующим изменениям мгновенных значений углов атаки и скольжения относительно средних значений, определяемых выражением (2)

Δ α ( t ) = − Δ θ ( t ) sin γ , Δ β ( t ) = Δ θ ( t ) cos γ .                 ( 11 )

Компоненты угловой скорости имеют следующий вид

ω y = Ω ( t ) cos γ ω z = − Ω ( t ) sin γ                     ( 12 )

Выражения для производных от углов атаки и скольжения в случае колебаний по рысканию определяются из вышеприведенных формул (12) с учетом кинематических соотношений (7):

α ˙ = Ω ( t ) ( − sin γ + cos γ sin α t g β ) , β ˙ = Ω ( t ) cos γ cos α .                               ( 13 )

В линейном приближении коэффициенты сил и моментов, действующих на модель при колебаниях по рысканию, с учетом (11), (12) и (13) можно представить в виде:

C = C 0 ( − C α sin γ + C β cos γ ) Δ θ + [ − ( C ω ¯ z + C α ˙ ¯ ) ] b a V sin γ + ( C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β ) l 2 V cos γ ] Ω .

Таким образом, по аналогии со случаем колебаний по тангажу, при малых колебаниях по рысканию каждый канал тензовесов также воспринимает некоторое среднее значение C₀ и сигналы в фазе с опорным сигналом Δθ(t) и в фазе с сигналом Ω(t)=Δθ(t):

C 1 ω y = − C α sin γ + C β cos γ .             ( 14 )

С учетом выражений (9) и (14,) решая систему уравнений для нахождения производных C^α и C^β, получим

C 2 ω y = − ( C ω ¯ z + C α ˙ ¯ ) b a V sin γ + ( C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β ) l 2 V sin γ ,     ( 15 )

где C^α - статическая производная аэродинамического коэффициента от угла атаки, C^β - статическая производная аэродинамического коэффициента от угла скольжения, C 1 ω z - составляющая сигнала в фазе с углом колебаний модели при колебаниях относительно оси Oz, C 1 ω y - составляющая сигнала в фазе с углом колебаний модели при колебаниях относительно оси Oy. Коэффициент C 1 ω z получен методом линейной регрессии при колебаниях модели по тангажу для заданных средних значений углов атаки и скольжения, а коэффициент C 1 ω y - при колебаниях по рысканию.

Аналогично, с учетом выражений (10) и (15), можно найти комплексы производных

C ω ¯ y   +   C α ˙ ¯   =   C 2 ω z   cos γ − l 2 b a C 2 ω y   sin γ ,

C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β = 2 b a l C 2 ω z sin γ + C 2 ω y cos γ .

В этом выражении C ω ¯ z + C α ˙ ¯ - комплекс продольных вращательной и нестационарной аэродинамических производных; C ω ¯ y + C β ˙ ¯ cos α + C α ˙ ¯ 2 b a l sin α t g β - комплекс боковых вращательной и нестационарных аэродинамических производных, C 2 ω z - составляющая сигнала в фазе с угловой скоростью колебаний модели при колебаниях относительно оси Oz, C 2 ω y - составляющая сигнала в фазе с угловой скоростью колебаний модели при колебаниях относительно оси Oy, l - размах крыла модели, b_a - САХ модели.

Обрабатывая совместно результаты эксперимента при колебаниях по тангажу и рысканию, можно найти все необходимые комплексы стационарных и нестационарных производных при заданных средних значениях углов атаки и скольжения.

При колебаниях по крену для установочного угла γ≠0 отлична от нуля только одна компонента угловой скорости движения модели так же, как и в случае γ=0

ω x = Ω = Δ θ = ∈ 0 ω cos ω t .

Т.е. для таких колебаний начальный поворот на некоторый постоянный угол крена не приводит к дополнительным взаимосвязям между комплексами аэродинамических производных.

Выражения для производных углов атаки и скольжения имеют вид

α ˙ = − Ω ( t ) cos α t g β ,

β ˙ = Ω ( t ) sin α .

Следовательно, в случае малых колебаний модели по крену коэффициент пропорциональности C 2 ω x составляющей сигнала тензовесов, опережающий по фазе на π/2 опорный сигнал, может быть выражен следующим образом

C 2 ω x = C ω ¯ x + C β ˙ ¯ sin α − C α ˙ ¯ 2 b a l cos α t g β .

Поэтому исследовать зависимость демпфирования крена от угла скольжения можно, не решая дополнительных систем линейных уравнений.

Предлагаемый способ реализуется при помощи устройства для определения статических и нестационарных аэродинамических производных моделей летательных аппаратов, собираемого в одном из трех вариантов для колебаний модели относительно одной из связанных осей (фиг.1). Устройство содержит поворотный круг 1 в рабочей части аэродинамической трубы (фиг.2), стойку 3 на поворотном круге, подвижную Г-образную раму 2, державку 6 модели с тензовесами 7, смонтированную в подшипниковом узле рамы 11, привод 5, связанный тягами 4 с рамой или державкой (в зависимости от конфигурации устройства), датчики угла поворота рамы 10 и державки 6. При колебаниях модели по тангажу Г-образная рама 2 устанавливается в вертикальной плоскости (фиг.1), при колебаниях по рысканию подвижная Г-образная рама 2 переводится в горизонтальное положение, при этом используется стойка 3 большей высоты. При колебаниях модели по крену тяга от привода крепится не к Г-образной раме 2, а к рычагу-зажиму 19 механизма 8. Для дистанционно управляемого поворота державки на заданный угол крена в устройстве использован механизм дистанционной установк