Алгоритм приближенного расчета гидродинамических сил, действующих на гидросамолет при продольной качке Текст научной статьи по специальности «Физика»

Секция летательных аппаратов

УДК 629.735.33

С.Г. Муганлинский

АЛГОРИТМ ПРИБЛИЖЕННОГО РАСЧЕТА ГИДРОДИНАМИЧЕСКИХ СИЛ, ДЕЙСТВУЮЩИХ НА ГИДРОСАМОЛЕТ ПРИ ПРОДОЛЬНОЙ

КАЧКЕ

Расчет параметров продольной качки гидросамолета, находящегося на воде,

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

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

,

, -

лета, со скоростью V в идеальную несжимаемую жидкость с учетом сил тяжести. Пусть у = / (х) - уравнение боковой формы клина, а у1 = /1 (х) - уравнение формы поверхности жидкости.

Запишем граничные условия в неподвижной системе координат 01Х1У1, связанной с поверхностью невозмущенной жидкости.

Смоченную часть клина определяем величиной 2с, условиями: |х| < С, — к < у < / (с) — к, где к - величина погружения клина.

Моделируем смоченную часть контура 8А, и свободную границу 8 вихревой

пеленой [1], а касательную скорость жидкости представляем с помощью интеграла Коши [3]. Пусть закон движения плоского контура соответствует продольной качке гидросамолета на известной регулярной волне, а также - движение деформации смоченной части поверхности контура, вызванное начальным ударом о жидкость. Тогда величины и 8 являются функциями времени: ((), 8 (() . При этом на свободной границе за плоским контуром 8;, (() образуется вихревой гидродинамический след <г(( )с 8((). Граница 8;, и 8 имеет конечное число угловых точек. За положительные направления ортов нормали п и касательной 8 на границе принимаются соответственно направления внешней к жидкости нормали и положительного обхода полуплоскости, заполненной жидкостью. Орты п и 8 образуют правую плоскую систему координат. В начальный момент времени ось X совпадает с невозмущенной свободной границей. Скорость точек смоченной части

контура и 8 определятся еле дующим образом:

и5 ((, Г ) = У (() + У^ ((, г), (1)

где У0 () - скорость движения начала связанной с клином системы координат; г - радиус-вектор точек смоченной части контура; У^ - скорость деформационного движения точек контура.

Ставится задача определения скоростей частиц жидкости, вызванных движением контура и движением профиля свободной границы, обусловленное начальной волной и возмущенным движением жидкости от движущегося контура. Предполагается потенциальность поля скоростей во внутренних точках жидкости. Искомый потенциал скоростей Ф((, г) является гармонической функцией во всех внутренних точках жидкости и удовлетворяет на границе следующим условиям: непротекания смоченной поверхности; совместности движения; непрерывности давления

на свободной границе; конечности скоростей на смоченной поверхности (условие - );

; .

Введя комплексный потенциал скоростей жидкости, можно показать, что на смоченной поверхности клина и на свободной поверхности скорость жидкости определяется через логарифмический потенциал двойного слоя:

U±(Z) = ± 1Y(Z)-1 J r(Z)c0s^dZ Ze SK u(2) 2 2 J r

sr us

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

Uc (C) = Vos (Z) - srgn(Vos (C)ia- (( (C)- Von (Z)) " 2g( " У. ), (3)

где VOS (C),VOS (TO) - касательная составляющая скоростей в точке Ze S и в бесконечно удаленной точке невозмущенной свободной границы; VON (Z) - нормальная составляющая скорости V0 ; y(t), yx (t) - ординаты точек свободной поверхности возмущенного и невозмущенного движения.

С учетом выражения (2) уравнение (3) представляется в виде интегрального уравнения относительно плотности вихревого слоя границы SK U S ;

f

y(C) = 2

1 f YMc0s( >n )

2 J y((- Uc (Z) V S. Us '

Ze S. (4)

Вполне естественно различать плотность у^) вихревого слоя для различных элементов границы. Обозначим: Усм, усд, Усл - плотности вихревого слоя смоченной части клина, свободной поверхности и свободного вихревого следа соответственно. Допускаем, что плотности вихревого слоя Уш, усд, Усл не зависят , .

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

С учетом такого представления плотностей вихревого слоя использование математического выражения динамического условия на свободной границе условия непроницаемости на смоченной части контура, кинематического условия совместности движения на свободной границе и условия теоремы Томсона [3] при предположении о начале движения из состояния покоя, дает систему интегральных уравнений относительно плотностей вихревого слоя уси, усв, уа и скор ости возмущенного движения жидкости на свободной поверхности:

( -\ 1 Г ( ,^cos(r,п) „ 1 (V / Л ,п) „

Пе (, " ) =--Ьс ( I)-^ ^--Г(( ( I) + 7,г (I))-^ М +

п • Г П * "

ьк

+ (°°))' V,2 (()■

I ¡гЛ, л+1 I)+гЛ)))¡^ <П—V, )

П ' г к * г

\ 8к 8К У

+

+я С/1(г)—уЛ^1;

\С(/) — z\; е 8((). (5)

г / JчCos(r, 8) „ 1 г/ / л /JччCos(r, 8) „

¡ГсЛ*,/)—м/+- Г/)+ /Л/))—М/ =

i г ж 8 г (6)

= 2т(,")(( () + Уа(,"));г е 8„.

{г^, /)/ + )/ = 0; (7)

8, 8

§ = г)—2/)8—V() г = 8(). (8)

Участок смоченного контура определяется боковой формой погружающегося клина у = / (х) и передней точк ой этого участка, которую можно определить как точку пересечения контура клина у = / (х) с контуром брызговой струи у = /1 (х). Так как контур брызговой струи неизвестен, то для построения числового расчета погружения клина можно поступить следующим образом. В начальный момент погружения клина на основании имеющихся экспериментальных или опытных данных построить контур брызговой струи в виде кривой плавного сопряжения контура клина со свободной поверхностью жидкости. Таким образом, может быть найдено начальное положение передней угловой точки смоченной части контура клина. В дальнейшем, в каждый последующий рассматриваемый момент времени (( + Дt) положение передней кромки смоченного участка может

г

быть определено путем смещения передней точки вперед на величину, определяемую решением уравнения У0 At = / (х).

При продольной качке гидросамолета на воде скорости поступательного движения точек поперечных сечений лодки связаны выражением У0 = и1х, где и

- угловая скорость, а 1х - расстояние от рассматриваемой точки до центра масс .

воде можно задать ряд значений угловой скорости и (круговая частота колеба-), . -лее общим будет подход, при котором частота колебаний определяется из решения динамических уравнений продольного движения гидросамолета. Эти уравнения по существу должны дополнить систему интегральных уравнений (5) - (8), так как они содержат неизвестные гидродинамические силы, действующие на корпус гид.

Решение системы интегральных уравнений (5)-(9) каким-либо способом, позволяет осуществить определение плотности вихревого слоя и с помощью интегра-- :

Лр(( х у ) = р(, х у)- ро

-ч (х, у)

= Р (-((, х, у)-У0 ))2 - 2(у()-у „(())- 2 д ,г )г ' (9)

(хо, у0)

где (х0, у0) - критическая точка, в которой скорость равна нулю; (х, у) - текущая точка на контуре ; р(х, у) - давление в текущей точке смоченной части контура; ро - давление на свободной поверхности жидкости.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Логвинович Г.В. Гидродинамика течений со свободными границами. Киев: Наукова Думка. 1968.

2. Муганлинский СТ. Некоторые проблемы проектирования формы днища гидросамолетов // Сб. докладов IV научной конференции по гидроавиации «Гидросалон-2002», 2002.

3. Кочин Н.Е. и др. Теоретическая гидромеханика. Т.1 М.: Огиз. Гостехиздат. 1948.

УДК 623.746-519

С.Г. Муганлинский, ОЗ. Носко

ИССЛЕДОВАНИЕ ПРОФИЛЕЙ МИКРОБЛА СХЕМЫ "ЛЕТАЮЩЕЕ КРЫЛО"

В настоящее время в ряде стран мира ведутся разработки микробеспилотных летательных аппаратов (МикроБЛА) [1].

- 1 -диусом действия до 1 км [2].

В ходе исследований были определены летно-технические и массовые характеристики перспективного МикроБЛА. В настоящее время ведутся работы по соз-