ISSN 1992-6502 (Print)_
2018. Т. 22, № 1 (79). С. 58-66
Вестник УГАТУ
ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru
УДК 533.657
Математическая модель качки гидросамолета на режиме дрейфа
12 3
а.в. Власов , а.н. Варюхин , м.а. Овдиенко
1 vav-94@yandex.ru, 2 a.varyukhin@yandex.ru, 3 m.a.ovdienko@gmail.com
1 ФГУП «Центральный аэрогидродинамический институт имени профессора Н. Е. Жуковского» (ЦАГИ) 1-3 ФГАОУ ВО «Московский физико-технический институт (государственный университет)» (МФТИ)
Поступила в редакцию 28.03.2018
Аннотация. Разработана математическая модель динамики пространственной качки гидросамолета в условиях регулярного и нерегулярного волнения на режиме дрейфа. Модель основана на подходах, применяемых при решении задачи качки судов, но с учетом существенных различий в геометрических соотношениях между судном и гидросамолетом. Учитываются только вертикальные колебания и колебания по углу дифферента и крена. С целью верификации матмодели были выполнены расчеты динамики бортовой и килевой качки модели гидросамолета в условиях регулярного волнения и выполнено сопоставление результатов расчета и эксперимента в опыто-вом бассейне. Получены некоторые результаты поведения гидросамолета в условиях нерегулярного волнения, позволяющие спрогнозировать предельно допустимый уровень волнения.
Ключевые слова: гидросамолет самолет-амфибия; нерегулярное волнение; качка; высота волны; мореходность; линейная модель волнения; опытовый бассейн; матмо-дель; метод плоских поперечных сечений.
ВВЕДЕНИЕ
Работа посвящена построению математической модели качки гидросамолета на режиме дрейфа в условиях волнения. Для самолетов-амфибий и гидросамолетов этот режим очень важен, т.к. значительную часть времени они проводят на воде, а морская поверхность крайне редко бывает спокойной. Основными характеристиками для гидросамолета в этом режиме являются амплитуды и частоты вертикальных и угловых колебаний.
Задача о колебании гидросамолета на воде несколько отличается от качки корабля, хотя и имеет много общего. Это связано в первую очередь с существенными различиями в
геометрических соотношениях между кораблем и гидросамолетом. Несмотря на большое число работ, посвященных поведению судов на волнении, качка
гидросамолетов обделена вниманием. Рис. 1. схема наоегающеи волны
Особенно это касается опубликованных экспериментальных данных.
В работе рассматривается гидросамолет стандартной схемы, типа летающая лодка. Учитываются только вертикальные колебания и колебания по углам дифферента и крена. Продольная скорость самолета считается около нулевой, т.е. моделируется режим дрейфа и маневрирования. Волна набегает на самолет под углом к продольной оси (рис. 1). Далее будет отдельно рассмотрена модель морского волнения.
МОДЕЛЬ МОРСКОГО ВОЛНЕНИЯ
Нерегулярное морское волнение в линейной постановке можно представить в виде бесконечной суммы гармонических волн [1] и [3]. Мы же, приближенно, будем считать эту сумму конечной
n(t,х,z) = aik cosjQilkt-
=1 k=1
ш,.
--— [ х соб(у - ) + 2 -
£
)]+ Ф* }
где у-^ш - разница между генеральным
направлением распространения волн и
направлением распространения заданной
составляющей, щк - частота, ф— - фаза,
a
амплитуда составляющей волны.
aik =V2 S (®ik, Wuk )АШ-А^1
При этом
ш -ш Аш=- max min
N..
А^1 =
^1max -Vhm
N
Существуют различные модели спектра морского волнения. В данной работе мы будем использовать выражение, предложенное в работе [2].
S(ш,V1) = -9,43DfeO е" п ш I ш J
cos у1,
где В - дисперсия волнения, Ю - средняя частота волнения, а ют = 0,77 Ю - частота на которой двумерный спектр волнения £(ю,щ) имеет максимум (рис. 2).
В отечественной практике принято описывать нерегулярное волнение высотой волны 3% обеспеченности Л3%. Упрощенно это можно характеризовать тем, что при
наблюдении за волнением 3% волн будут иметь высоту не менее, чем Л3%.
Дисперсия и средняя частота зависят от высоты волны 3% обеспеченности следующим образом
D = 0,143
гк
"3%
2п
Ю =- ,-
3,3\ ^3%
Выполнив математические выкладки, можно перейти от двумерного энергетического
спектра к одномерному
^ Ю2
= ^ а ^ cos{югí---[хcos(y-
¿=1 у=1 ё
j) + z sin(V - j)] + 9i} cos2 j,
где
2 _
ai = — yj2 • S1 (ш i )Аш • Ау1
п
Аш =
2ш m
N...
А^1 =
п
N
V1
при этом
^(ш) = 9,43 Df^m| е" ш I ш J
0.05 0.04 0.03
S1 Ü.02 0.01 о
w/wm
Рис. 2. Одномерный энергетический спектр волне-
ния h3%=1 м
В данной модели фаза ¿-ой составляющей волны ф { берется случайным образом в пределах от 0 до 2п. Случайная частота Юг. ¿-ой составляющей имеет пределы от 0,5 ю до 2,5 ют, а угол распространения у у-ой составляю-
щей меняется от
п п
до 0 относительно угла
2 ди 2
генерального распространения у. Все случай-
ные величины берутся исходя из равномерного распределения.
УРАВНЕНИЯ ДИНАМИКИ
Основой для расчета динамики движения гидросамолета при качке на воде является гипотеза поперечных плоских сечений, смысл которой заключается в том, что гидродинамические силы, действующие на каждое отдельное сечение, не зависят от других сечений. Рассмотрим режим плавания гидросамолета, одно из поперечных сечений которого схематично представлено на рис. 3.
Одним из основных допущений является то, что мы пренебрегаем поперечным моментом, возникающем на фюзеляже, по сравнению с моментом от поплавков. Это обосновано тем, что расстояние от плоскости симметрии самолета до поплавка значительно больше ширины фюзеляжа, соответственно, и возникающий на поплавке поперечный момент намного превосходит поперечный момент фюзеляжа.
Вертикальные размеры поплавков и расстояние между уровнем спокойной воды и килем поплавка соизмеримы с высотой волн и размахом вертикальных колебаний. Поэтому колебания гидросамолета по крену
нельзя рассматривать без учета вертикальных колебаний. Уравнения движения по крену существенно нелинейны. При рассмотрении колебаний на значительном волнении, когда высота волны соизмерима с осадкой гидросамолета, необходимо учитывать нелинейность в уравнениях, определяющих вертикальную и килевую качку.
Вертикальная сила, действующая на заданное сечение фюзеляжа, определяется следующим выражением
ё (т •$)
/($, $,$) = /ге + /д +-
ж
где $ = п- У ~ погружение сечения, /гс -гидростатическая составляющая,
/д=2АрЪ(£)£ - демпфирующая составляющая, пропорциональная ширине погружен-
ТГ
ной части 2Ъ($), т = — рЪ2 ($) - присоединенная масса сечения. Суммарная вертикальная сила 7ф, действующая на фюзеляж
гидросамолета, вычисляется интегрированием по смоченной длине
'к
Гф = | / (х)ёх.
(1)
н
1 —1
1 нсг X г ^
Суммарный момент:
-к
мф = { f(x) ■xdx=
(2)
где lK и lH - длины кормовой и носовой части.
Величина силы f (x), действующей на поперечное сечение, и присоединенной массы m определяется различными случаями погружения фюзеляжа в воду:
1) киль фюзеляжа в воздухе
£(x) < y (x),
f (x) = 0, m = 0;
2) киль в воде, скула в воздухе
yk (x) <£(x) < yc (x),
f (x) = (£-Ук )[ £ (£-Ук ) +
tgP
+(A +—i )£ + — n ] -v tg^ 2tgp ,J
(3)
pn
tg 2P
(£-Ук )2(y +x Q),
m = p
n
2tg 2p
(£- Ук )2;
(4)
3) скула в воде, верхняя часть сечения фюзеляжа в воздухе
Л(х) <Е(х) < Уф (х),
/ (х)=Ръ [ ^ ус))+
+4 + П ЬС1Л ]-рПъ2(У+ хэ), (5)
П и 2
m = p- ¿с;
(6)
4) сечение фюзеляжа полностью погружено в воду
x) > Уф(x), f (x) = pbc[g(bctgp + 2(Уф -ус)) +
П , .... П
+A4 + - bcfj ] -p- bc2(yy + xQ),
П A 2
m = p- bc.
(7)
(8)
Вертикальные силы, действующие на подкрыльевые поплавки, вычисляются аналогично силам, действующим на фюзеляж. В то же время, т.к. обычно длины волн существенно больше длины поплавков, можно использовать приближенную методику и
уйти от необходимости суммирования составляющих сил в каждом сечении. Будем аппроксимировать контур сечения поплавка Ъ(Еп) синусоидой.
.П
(9)
Ъ(Е ) = Ьп в1и(-(4 - У,)) К
при укп <£п < ув, Ъ(£ ) = 0 в остальных случаях
где ук и ув - координаты нижней и верхней точек поплавка по оси У.
Получаются следующие выражения для определения сил, действующих на поплавки:
1) поплавок в воздухе
Е < ук , У = 0
^ п У кп > п
2) часть поплавка в воде
Укп <Еп < Ув ,
Y =pln[gVn + bi Ai„ +Kbbin + + П ¿ft]-pf b12ln (У + xn Q + ^пY)
(10)
где Vn = 2bnh П
П
1 - cos(-(£n - У к)) h
П
b = bn sm(—(^n -Ук)),
hn
b =-T bn c°s(f (£n - У к ))£ n ,
hn hn
h
bi =b пРи £n < Ук + у
b1 = bn при ^n ^ У к +
h
3) поплавок полностью погружен в во-
ду
£ ^ У ,
->п У в 5
h
Уп =pln [4gbn — + bn A£ n +
П
П
П
(11)
+^ b12ii ] -P ^ Kh (y + xnQ+)
Обозначим за Упл и Уппр силы, действующие на левый и правый поплавок соответственно. Б - вес самолета, К - высота
" цм
его центра масс. Далее из (3)-(11) переносим члены, содержащие вторые производные у, 0 и у, в правую часть, в результате чего получим систему уравнений, которую
H
окончательном виде можно записать следующим образом
а^у + а12ё + а1зУ = у \ + у '„л+ у ^ - £,
а2\У + а22ё +а2зУ = М \ф + (У 'пл + ■У 'ппр )Хп ,
а31у + аз2ё + аззУ = (У -У К + ^ -кцм sin у ,
где
Б ?
а,, =—+ I Ат(£)йх + Атш + Ат]
8 I
а22 = + | Ат^) х2 йх + х^ (АтпЛ + т,X
-г.
8
Б 2 2
азз =— К + Атпл +Атппр) 8 р
-г,
а12 а21
I Ат(^)хйх - хп (Атпл + ^ппр )
а,, = а,, = г (Ат -Ат ),
1з з1 п \ пл ппр / '
а„ = а„ = х г (Ат - Ат ).
2з з2 п п\ пл ппр /
Причем штрихованные силы и моменты для фюзеляжа будут определяться по (1)-(2) почти так же, как в (з)-(7), только без членов, содержащих вторые производные у, 9 иу
1) при х) < ук (х)
I (х) = 0;
2) при Ук(х) <£(х) < Ус(х)
I ( х) = (£-Ук )[ 8 (£-Ук ) + (А + tgв
п п ....
+—£)£ +—п];
3) при Ус(х) <£(х) < Уф(х)
I ( х) = рЪс [ 8 ФМ +
+2(£-ус)) + А^ + П ЪсП];
4) при £(х) > Уф (х)
I ( х) = рЪс[ 8 (ЪМ +
п
+2(Уф -Ус))+А£ +2ЪП
Аналогично для поплавков из (9)-(11) получим:
1) пРи £ п < Ук
У\ = 0 ;
2) пРи У, <£„ < У,
У'п = р1п [ 8Уп + Ь А£п + п 2
к. • 2' к
Ъ1 = ъп пРи £п ^ У, + у;
+пЪЪ£п Ъ12п];
Ъ1 = Ъ при ^п <У, +-
з) при £п ^ У У' =р/
п г }
4ёЪп — + ЬпА£п +п Ъ:2Л п 2
ДОПОЛНЕНИЕ К УРАВНЕНИЯМ ДВИЖЕНИЯ
Можно учесть присоединенную массу воздуха горизонтального оперенья. Пусть 1к - длина крыла, а кк - его средняя аэродинамическая хорда, тогда для оценки можно пренебречь формой крыла и считать его прямоугольным. Очевидно, что крыло может дать поправку только на момент крена, не делая при этом вклада в продольный момент.
Распишем силу и момент, создаваемые крылом
к п
рк = | р 2 к1 •(у + х у) йх =РпЧ1к у
С п г
Мк = | х-Р-к • (у+ ху)йх =Рпк^у. г 2 з
Забегая вперед, даже грубые оценки силы Ек показывают, что она значительно меньше остальных сил, действующих на самолет, и не может оказать существенного влияния на колебания самолета. Аналогичная сила, создаваемая задним горизонтальным опереньем, существенно меньше ввиду его малых размеров. Поэтому основная поправка, которую оно может внести, это поправка на продольный момент.
п
п
к
У, м
0.3 0.2 0.1 о
У, м 0.4 Г
0.2
-О. В -0.6 -0.4 -0.2 О 0.2 0.4 0.6 0.8
Рис. 4. Визуализация волнения и положения гидросамолета в продольной плоскости
О
-02
а — ■ - —4-
й ■- -с?. и. г;: ',:ъ х, м
Рис. 5. Визуализация волнения и положения гидросамолета в поперечной плоскости
Приближенно можно записать
Мго =9ПИ2го1го ■ ^ (у + Гг£\
где гоо - проекция расстояния от центра
масс самолета до заднего горизонтального оперенья на горизонтальную плоскость самолета.
Численные расчеты показали, что добавление этих моментов не оказывает сколько-нибудь заметного влияния на итоговое решение, поэтому в численных расчетах можно обойтись без них.
Система уравнений решается численно методом Рунге-Кутта 4-го порядка, для ее решения была написана вычислительная программа в среде МайаЬ. В программе предусмотрена возможность визуализации (рис. 4 и 5), а также вывод графиков изменения вертикального положения у и углов
дифферента 0 и крена у (рис. 5-11).
СРАВНЕНИЕ РАСЧЕТНЫХ ДАННЫХ С ЭКСПЕРИМЕНТОМ
Для верификации данной модели была проведена серия экспериментов в опытовом бассейне ЦАГИ. Модель самолета-амфибии в масштабе 1:16 подвергалась воздействию
поперечного и продольного волнения. В расчете данное волнение считалось регулярным и задавалось в виде синусоиды с заданной амплитудой. Сравнение результатов эксперимента и расчета показали хорошую точность данной модели (рис. 6-11). На рис. 6-9 представлены два опыта и расчеты для продольной качки, а на рис. 10-12 два опыта и расчеты для поперечной качки. Небольшие отклонения могут быть связаны с ошибками в измерениях параметров модели самолета, а также с тем, что данная математическая модель подразумевает свободно плавающий гидросамолет, что не совсем верно для проведенного эксперимента.
Рис. 6. Вертикальные колебания
о
-5 -10 -15
5
5-20 -25 -30 -35 -40
Рис. 7. Колебания по углу дифферента
Рис. 8. Вертикальные колебания
Рис. 9. Колебания по углу дифферента
2 4 к ; 1 0 12
а . А д
\ 1 Ч \
1 /
-Эксперимент _Расчёт_
С с
2,5
лДААЛГ
-Экспе
--Расчё
-Эксперимент -Расчёт
10 12
Рис. 11. Колебания по углу крена.
Рис. 12. Вертикальных колебаний.
КОЛЕБАНИЯ В УСЛОВИЯХ НЕРЕГУЛЯРНОГО ВОЛНЕНИЯ
Были проведены численные расчета модели самолета-амфибии в масштабе 1:16 на нерегулярных волнах с высотами 3% обеспеченности И3%о = 0,0625, = 0,1,
= 0,125 и углами генерального распро-
п п
странения волны у = 0, У = ^, ^=2"
(рис. 13-16). Нерегулярная волна считалась в соответствии с формулами из п. «Модель морского волнения», N=50, N=10.
Рис. 10. Вертикальные колебания
Рис. 13. Колебания волны к3%А = 0,1, у = 0
0,14
Рис. 15. Колебания по углу дифферента к3%А = 0,1, у = 0
0,25
0
О 2 4 6 8 10 12
I, С
Рис. 16. Колебания по углу крена Иш = 0,1, у = 0
ДОПУСТИМЫЕ ВОЛНЫ ДЛЯ МОДЕЛИ САМОЛЕТА-АМФИБИИ
Наиболее важный результат, который позволяет получить данная программа - это высоты волн, при которых происходит замыв волной частей самолета: хвоста или крыла. Этот параметр определяет максимально допустимое волнение, а как следствие, и мореходность гидросамолета.
В табл. 1 представлены максимально допустимые высоты волн 3% обеспеченности для модели современного самолета-амфибии при разных направлениях генерального распространения волн.
Таблица 1 Максимально допустимые высоты волн 3% обеспеченности для модели современного самолета-амфибии
Угол генераль- Макс. допус- Макс. допус-
ного распро- тимая волна тимая волна
странения вол- 3% обеспечен- 3% обеспечен-
ны ности для мо- ности для на-
дели туры
у = 0 0,18 м 2,9 м
П 0,18 м 2,9 м
у = -4
П 0,22 м 3,5 м
у = у
ЗАКЛЮЧЕНИЕ
Построена математическая модель волнения и качки самолета и выполнена ее верификация на экспериментальных данных. Эксперимент показал, что данная математическая модель достаточно точно описывает поведение самолета на волне, что подтверждает правильность допущений, сделанных при разработке модели. Также была выполнена серия расчетов качки гидросамолета на нерегулярном волнении с разными высотами волн 3% обеспеченности и разными направлениями генерального распространения.
В ходе анализа этих расчетов было выявлено, что с увеличением угла генерально-
« п
го распространения от 0 до увеличивается амплитуда колебаний по крену и уменьшается амплитуда колебаний по углу дифферента. На вертикальные колебания угол генерального распространения волн практически не влияет. Увеличение высоты полны приводит к увеличению всех видов колебаний. Эти закономерности вполне очевидны и могут выступать лишь в роли еще одного подтверждения корректности математической модели, построенной в данной работе.
СПИСОК ЛИТЕРАТУРЫ
1. Справочник по теории корабля / Войткунский Я. И. и др. в 3 т. Л.: Судостроение, 1985 г. [Y. I Vojtkunskiy «Handbook of ship theory", in Sudostroenie, 1985]
2. ОСТ 5.1003-80. Методика расчета качки водоизме-щающих кораблей и судов. [Электронный ресурс]. URL: https://www.twirpx.com/file/1642476/ (дата обращения 12.03.2018) [ OST 5.1003-80. (2018, Mar. 12) Method of calculation of pitching of displacement ships and vessels, (2018, Mar. 12) (in Russian) [Online]. Available: https://www.twirpx.com/file/1642476/ ]
3. Бородай И. К., Нецветаев Ю. А. Мореходность судов // Судостроение, 1982. [I. K. Boroday, U. A. Necvetaev "Seakeeping of vessels", Sudostroenie, 1982]
towing tank; mathematical model; method of cross sections
About authors:
VLASOV, Artur Vadimovich, student of the Department "Theoretical and applied Aero-hydrodynamics" (MIPT, 2018), hydrodynamics of marine aircrafts, planing, fast entry of bodies into the water.
VARYUKHIN, Anton Nikolaevich, senior researcher of the laboratory of Autonomous systems (MIPT 2005), hydrodynamics of marine aircrafts, high-speed hydrodynamics.
OVDIENKO, Maxim Aleksandrovich, post-graduate student of the Department "Theoretical and applied hydrodynamics" (MIPT 2012), hydrodynamics of marine air-crafts.
ОБ АВТОРАХ
ВЛАСОВ Артур Вадимович, студ. каф. «Теоретическая и прикладная аэрогидродинамика» (МФТИ(ГУ), 2018). Иссл. в обл. гидродинамики морских летательных аппаратов, скоростной гидродинамики, быстрого входа тел в воду.
ВАРЮХИН Антон Николаевич, канд. техн. наук, ст. науч. сотр. «Лаборатории автономных систем» (МФТИ(ГУ), 2005). Иссл. в обл. гидродинамики морских летательных аппаратов, скоростной гидродинамики.
ОВДИЕНКО Максим Александрович, асп. каф. «Теоретическая и прикладная аэрогидродинамика» (МФТИ(ГУ), 2012). Иссл. в обл. гидродинамики морских летательных аппаратов.
METADATA
Title: Mathematical model of the pitching seaplane on the mode drift
Authors: A. V. Vlasov 1, A. N. Varyukhin2, M. A. Ovdienko3
Affiliation:
1 Central Aero hydrodynamic Institute, Russia.
1 -3 Moscow Institute of Physics and technology, Russia.
Email: 1 vav-94@yandex.ru, 2 a.varyukhin@yandex.ru,
3 m.a.ovdienko@gmail.com
Language: Russian.
Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), ), vol. 22, no. 1 (79), pp. 58-66, 2018. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).
Abstract: The mathematical model of the dynamics of pitching seaplane under conditions of regular and irregular excitement on the mode drift is developed. The model is based on the ways that is used in solving the problem of pumping ships, but it considers the significant differences in the geometry between ships and seaplanes. Only vertical vibrations and oscillations of the angle of pitch and roll is considered. For verification of matmodel calculations of the dynamics of rolling and pitching motion of the model of a seaplane in conditions of the regular waves were made and the results of calculation and experiment in the testing pool was compared. Some results of the behavior of the seaplane under conditions of irregular excitement, allowing to predict the maximum permissible level of excitement, are obtained.
Key words: seaplane; amphibious aircraft; irregular waves; pitching; wave height; seakeeping; linear wave theory;