Математическое моделирование колебательных движений крыла с гибким профилем

НАУЧНЫЕ СТАТЬИ И ДОКЛАДЫ • МЕХАНИКА

УДК 533.601.1

Ю. С. КАБАЛЬНОВ, Л. Ю. УРАЗАЕВА

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ КОЛЕБАТЕЛЬНЫХ ДВИЖЕНИЙ КРЫЛА С ГИБКИМ ПРОФИЛЕМ

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

ВВЕДЕНИЕ

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

Предпринятые исследования в данной области [8] свидетельствуют об относительно высокой эффективности подобных систем. В настоящее время на основе композитных материалов появляется возможность практической реализации движения гибких крыльев переменной геометрии с конечными амплитудами. Существующие методы решения задачи о нестационарном движении крыла в основном посвящены случаю бесконечно малых амплитуд, либо представляют собой чисто теоретические подходы, трудоемкие для практического применения, либо ограничиваются частными случаями движения, движением крыла вблизи поверхности. Весьма актуальна в связи с этим разработка метода, позволяющего находить удобное для практического применения и в то же время достаточно обоснованное решение для использования в практических расчетах.

Наибольшие трудности в нестационарных задачах доставляет учет влияния вихревого следа, тянущегося позади крыла. Это связано

с тем, что заранее неизвестны ни форма вихревого следа, ни вихревая плотность.

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

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

ИСПОЛЬЗУЕМЫЕ ОБОЗНАЧЕНИЯ

ХУИ — абсолютная неподвижная система координат;

хуг — относительная, неизменно связанная с движущейся пластинкой;

— скорость горизонтальная;

— скорость вертикальная;

— угловая скорость;

— угол атаки;

— полухорда;

— полуразмах крыла;

— относительное удлинение (достаточно велико, );

— частота колебаний;

— амплитуда колебаний;

— плотность (интенсивность) вихря;

— интенсивность присоединенных вихрей профиля;

— интенсивность вихрей следа;

— число Струхаля;

— циркуляция полная;

— циркуляция квазистационарная;

— циркуляция дополнительная;

, — проекции силы, вызываемой вли-

янием вихревого следа;

, — ядро интегро-дифференциально-

го уравнения для случая крыла бесконечного, конечного размаха.

1. ОБЩАЯ ПОСТАНОВКА ЗАДАЧИ

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

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

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

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

От части кромки непроницаемой поверхности, которую называют задней кромкой, сходит поверхность разрыва скоростей, причем скорость жидкости остается конечной на задней кромке. Постулат о конечности скорости на задней кромке был введен в теорию крыла Жуковским-Куттом. Данное условие отражает влияние вязкости реальных жидкостей в теории крыла в идеальной жидкости. На поверхности крыла выполняется условие отсутствия относительных нормальных скоростей. Это связано с тем, что при движении крыла частицы жидкости, находящиеся на нем, непрерывно соприкасаются с его поверхностью. В зависимости от требования от-рывности или безотрывности обтекания передней кромки имеют место разные классы задач. С точки зрения тех видов движения, которые встречаются в природе, наиболее эффективным является безотрывное обтекание передней кромки. Это условие будет выполняться: 1) за счет подбора параметров, обеспечивающих отсутствие отрыва с передней кромки; 2) за счет соответствующего изменения формы профиля крыла, гарантирующего плавный безударный вход. Проблемой является также полный учет влияния вихревого следа. При достаточно плавных движениях, которые как раз и являются объектом изучения в работе, согласно теоретическим и экспериментальным исследованиям влиянием деформации следа можно пренебречь в рамках точности рассматриваемой модели. Таким образом, рассматривается приближенное реше-

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

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

На основе анализа результатов эксперимента можно сделать вывод о возможности применения принципа суперпозиции при расчете нестационарных сил, при наложении определенных ограничений на параметры движения: ограничения на величину рассматриваемого угла атаки. Наиболее распространенным методом разложения решения на нестационарную и стационарную составляющие является метод гармонического анализа [1].

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

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

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

2. МАТЕМАТИЧЕСКАЯ ПОСТАНОВКА ЗАДАЧИ

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

здесь и — вектор допустимых управлений; ф — вектор пространства состояний; E(Q.) — пространство состояний системы; —

область, занятая жидкой средой; S — ограниченная поверхность, движущаяся в жидкости.

Операторное уравнение Т(ф,и) = О представляет собой уравнение гидродинамики идеальной жидкости. В этой задаче </> имеет смысл потенциала скоростей или потенциала ускорений.

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

В данном случае рассматривается колебательное движение крыла в несжимаемой идеальной жидкости.

Предполагается, что крыло обладает поступательной скоростью , угловой скоростью вращения ее относительно задней кромки , хорда крыла может изгибаться для обеспечения плавного схода следа с задней кромки. Считается, что амплитуда колебаний крыла не превышает 0,2 — таким образом, гарантируется, что обтекание крыла носит безотрывный характер. Относительное удлинение крыла , равное отношению , берется достаточно большим: 4^ I/а, здесь I полуразмах крыла, — полухорда. Отметим, что в природе удлинения в основном встречаются в диапазоне от 4 до 12 [6,7].

Пусть , , — неподвижная система ко-

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

уравнение Лапласа.

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

д2Ф д2Ф д2Ф _

о" о" о"" \2/

ОХ/ д)Г OZz

вне несущей поверхности и поверхности вихревого следа, образующегося при нестационарном движении крыла;

дФ

—— = (ис — шу) cos (пх) + (vc + сох) cos (пу);

на поверхности несущей поверхности (условие непротекания);

Р+ = 11

дФ

дп -

с1Ф

с1п

— условие непрерывности давления и нормальной скорости при переходе через вихревую поверхность;

Нт\7Ф = О,

где К = \/ж2~~

у2 + г2 при К —> ос — условие отсутствия возмущений в жидкости на бесконечном удалении от крыла;

с1Ф

с1п

< ОС

на задней кромке (постулат Кутта-Жуковского).

Задача является нелинейной и состоит в определении потенциала при за-

данных начальных условиях.

3. ПРЕДПОЛАГАЕМЫЕ ПОДХОДЫ К РЕШЕНИЮ

При решении прямой и обратной задачи встает проблема построения конструктивных моделей и соответствующих численных алгоритмов для решения.

Для исследований в теории крыла в основном применяется вихревой метод, позволяющий наиболее полно отразить влияние вихревого следа, а следовательно, и физическую картину обтекания крыла.

Фактически прямая задача теории крыла сводится к решению интегро-дифференци-ального уравнения, получаемого из условия непроницаемости.

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

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

метод дискретных элементов (вихревых или дипольных панелей);

метод распределенной сингулярности (функций ядра или разложения функций).

Метод дискретных элементов ввиду большей легкости в реализации применяется чаще, но он требует осторожности при применении, так как «по мере приближения контрольной точки к вихревой нити, индуцируемые в ней скорости стремятся к бесконечно большим значениям» [4, с. 74].

Второй способ является более надежным, так как при численных расчетах выделяется сингулярность, но и более трудоемким. Данный подход является распространением исследований Н. Н. Поляхова на случай колебаний с конечными амплитудами.

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

Вихревая модель крыла конечного размаха рассматривается для относительно больших удлинений (от 4 до 8).

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

В случае крыла с постоянной циркуляцией по размаху используется гипотеза плоских сечений.

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

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

• движение «ныряющей» пластинки; движение круговой дуги; движение параболической дуги; движение профиля переменной формы.

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

4. ОПРЕДЕЛЕНИЕ ЦИРКУЛЯЦИИ И СИЛ

Для описания движения вводятся две системы координат: — абсолютная неподвижная система и — подвижная систе-

ма координат, неразрывно связанная с хордой крыла.

Логично начать исследование с рассмотрения случая колебаний с конечными амплитудами бесконечно тонкой пластинки (крыло с бесконечным размахом, т. е. с удлинением больше 8). В этом случае имеем плоскую задачу. Предполагается, что пластинка совершает одновременно поступательное движение со скоростью Vc = (u0,vmcosut) и вращательное движение относительно задней кромки с угловой скоростью , которая определяется из кинематических условий. Нестационарное движение с переменной во времени циркуляцией сопровождается сходом вихревого следа с задней кромки пластинки. Принимается, что вихри следа располагаются по траектории задней кромки, уравнение траектории имеет вид

х,- !-т  (

У О = -----------81111/ ---------------

Щ \ »'0

Гладкий сход вихревого следа обеспечивается законом изменения угловой скорости , который обеспечивает в лю-

бой момент времени расположение пластинки по касательной к траектории задней кром-ки(«ныряющая пластинка»).

Уравнение непроницаемости для пластинки запишется в виде

Vі — V

V пр п — У ПЛ п

Vі

' sn;

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

На основе вихревой модели выписывается уравнение сохранения циркуляции:

Г = 2a J Flipn(l -тнв)(1в + о

7Г

о

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

Для данной модели движения, полагая в первом приближении, что получим

(г „(п , п ос ( Vm\2w\

Гі =

Го^С + тга-0,25

Т

/

Г2 =

/Го^

' 7Г

V Т

m 1 т= —

7Г

(aS+Tr)'2+a2(C + Tr ( 1^0,25 ( — V V «о

Я]

С= [ т*°{Х ~ Хо) I dX:

J О,

хо

оо

S — ИХ-,

А" О

Г0 = -iraaV^l - 0,5 ).

Щ

Vrr.

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

Ядро I = 1(ро, 9?) для прямолинейной пластинки в случае криволинейного следа является функцией — расстояния от задней кромки до вихря следа с координатами , ,

<р — угла наклона пластинки к оси ОХ, таким образом, в ядре учитывается влияние конечности амплитуды и положения пластинки в текущий момент времени.

Ядро для криволинейной пластинки можно представить в виде:

р — ^11 р

1 + 0 — V V «о

где — ядро в случае прямолинейного следа, который имеет место при колебаниях с бесконечно малыми амплитудами, дополнительное слагаемое свидетельствует о квадра-

'г,

тичной зависимости от величины

«о

Для пластинки с изменяющимся по времени профилем, аналогом изменяющейся формы тела рыбы при движении, ядро в уравнении сохранения циркуляции будет содержать дополнительно множитель иметь вид: (1 — ^),

где — координата перехода прямолинейной части профиля в изогнутую, в любой момент времени известную из закона изменения профиля.

Следующим этапом исследования являются исследования крыла конечного размаха при удлинении от 4 до 8, совершающего колебательные движения с конечными амплитудами. В данном случае вихревая система крыла должна отражать конечности размаха крыла, кроме конечности амплитуды колебаний крыла. При исследовании движения с учетом конечных колебаний крыла с постоянной и переменной по размаху циркуляцией при достаточно больших удлинениях строится вихревая схема крыла на основе предложенной Н. Н. Поляховым [1]. В вихревую модель кроме поперечных вихрей, вводятся также продольные вихри, являющиеся как бы продолжением присоединенных вихрей. Однако в отличие от случая нестационарного движения с малыми амплитудами, кроме кри-волинейности формы следа при движении с конечными амплитудами, необходимо учитывать и возросшую нестабильность вихревого следа. В связи с этим, в модели вихревого следа и в расчетах полагается, что на определенном конечном расстоянии от задней кромки вихревой след сворачивается в вихревую шахматную дорожку Кармана.

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

г = 2о / (V,:

V'

х

здесь — полная скорость крыла, , —

индуктивные скорости крыла, Ао — относительное удлинение крыла. Члены порядка А-2 возникают за счет учета конечности размаха присоединенных крыла и влияния продольных вихрей на участке от передней до задней кромки.

Часть циркуляции, вызываемая влиянием поперечных вихрей следа, имеет вид

Ао

Гпоп = [ + 1п) йх';

= ~

7Г

х-2

Ао

^ (1 — coяjв) (ро соя((р — д) + 1 + сой в) ^ _

о ' о + !)

= 100 + А11:),

)\/(^

здесь — расстояние от точки сече-

ния до вихря следа; — вид ядра для крыла бесконечного размаха с учетом конечности колебаний.

Часть циркуляции, отвечающая за учет продольных вихрей, имеет вид:

Ао

Гцрол = — (-^1 + £>1)Г с1Х\

где

в, = 2

7Г

X

X СС>8 р

------008 (7--------- 81119? X

«о а

гг

X <

О

0о2Ао2 + 1)':

где — расстояние от точки сечения до вихря следа.

Полученные формулы в случае крыла с изменяющимся сечением отличаются заменой величин

на

7Г

. . , ,

81111Р на 8111 р I 1---------------------------------

7Г

В случае крыла с сечением в форме круговой дужки формулы остаются такими же, как для ныряющего крыла, за исключением ква-зистационарной части, где учитывается изогнутость дужки.

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

Для «ныряющего крыла» и крыла с сечением в виде круговой дужки (если пренебречь кривизной профиля сечения) величины проекций , , силы, вызываемой влиянием

следа, будут равны

Хя =

-руйпт’х

00 і I,

У, = ращ 17,, о

СОЙ <р

где

Х' =

Х^Хй

1 =

ЭХ'

ОХ' 2у/Л

Р2

В = Ро СОЙ (<р — (/) + -у СОЙ 2(<р — (?) + Ро\/Р.

где

р = 1

Ро /

-р + Ро СОЙ (<р

?)

Формулы для Х5, в случае крыла с постоянной циркуляцией по размаху при достаточно большом удлинении содержат члены, отвечающие как за влияние конечности размаха, так и за влияние конечности амплитуды колебаний. На больших расстояниях от задней кромки вклад криволинейности следа и конечности амплитуды колебаний для крыла конечного размаха ведет себя как

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

Полученные выражения для сил Жуковского легко распространяются на случай крыла с изменяющимся сечением.

Анализ результатов показывает, что основной вклад в учет конечности амплитуды колебаний вносят величины, отвечающие за влияние поперечных вихрей следа. Учет влияния конечности размаха ослабляет влияние конечности амплитуды с уменьшением удлинения крыла. Также при учете конечности размаха сокращается длина участка вихревого следа, вносящего основной вклад в величины проекций силы Жуковского на оси подвижной системы координат по сравнению со случаем крыла бесконечного размаха.

5. ПРИМЕР

Численное экспериментирование открывает большие возможности по изучению нестационарного движения несущей поверхности, предоставляет возможности для проведения численного имитационного эксперимента, однако в любом случае окончательные выводы об адекватности той или иной предлагаемой модели можно делать только на основе сравнения полученных результатов с экспериментальными данными. Проведем расчеты для случая прямолинейной «ныряющей» пластинки, движущейся с конечными амплитудами. Исходные параметры ут/и0 = 0,15, а = 6, Хо = 0 (отметим, что расчеты выполнены для определенного положения пластинки).

Таблица

Влияние частоты колебаний на величину полной циркуляции

Частота колебаний Число Струхаля Полная циркуляция

1 0,600 1,044

2 1,200 9,273

3 1,800 18,314

4 2,400 17,282

5 3,000 20,034

6 3,600 22,194

Как и следовало ожидать, полная циркуляция нелинейно возрастает с увеличением числа Струхаля при всех остальных одинаковых значениях параметров, но скорость роста полной циркуляции замедляется, начиная с определенных значений числа Струхаля, дальнейшее увеличение числа Струхаля уже не так эффективно для увеличения циркуляции. Данный результат хорошо согласуется с выводами экспериментальных работ [6,7]. Результаты численного эксперимента указывают на существенное проявление нелинейности в рассматриваемом диапазоне частот. С

ростом величины '

«о

влияние конечности

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

ВЫВОДЫ

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Поляхов, Н. Н. Избранные труды. Аэродинамика / Н. Н. Поляхов. СПб. : Изд-во СПбГУ. 1997.380 с.

2. Уразаева, Л. Ю. Учет влияния конечности размаха и конечности амплитуды колебаний при расчете сил, действующих на крыло / Л. Ю. Уразаева // Прикладная механика : межвуз. сб. СПб., 1997. Вып. 10. С. 215-216.

3. Поляхов, Н. Н. Гармонические колебания крыла с постоянной циркуляцией по размаху / Н. Н. Поляхов, Л. Ю. Уразаева // Вестн. СПб. ун-та. Сер. 1. 1992. Вып. 3, № 15. С. 102-104.

4. Белоцерковский, С. М. Отрывное и безотрывное обтекание тонких крыльев идеальной жидкостью / С. М. Белоцерковский, М. И. Ништ. М.: Наука, 1978.352 с.

5. Панченков, А. Н. Теория оптимальной несущей поверхности / А. Н. Панченков. Новосибирск : Наука, 1983. 256 с.

6. Taylor, G. K. Flying and swimming animals cruise at a Strouhal number tuned for high power

efficiency / Graham K. Taylor, Robert L. Nudds, Adrian L. R. Thomas // Nature. Oct. 16, 2003. Vol. 425. P. 707-711.

7. Tobalske, B. W. Dial kinematics of flap-bounding flight in the zebra finch over a wide range of speeds / Bret W. Tobalske, Wendy L. Peacock, P. Kenneth // The J. of Experimental Biology. 1999. Vol. 202. P. 1725-1739.

8. Katz, J. Hydrodynamic propulsion by large amplitude oscillation of an airfoil with chordwise flexibility / J. Katz, D. Weihs // J. Fluid Mech. 1976. Vol. 74, pt. 1. P. 485-497.

ОБ АВТОРАХ

Кабальнов Юрий Степанович, проф., зав. каф. информатики. Дипл. инж. электронной техники (УАИ, 1971). Д-р техн. наук по управлению в технических системах (УГАТУ, 1993). Иссл. в обл. адаптивного и интеллектуального управления.

Уразаева Л. Ю., доц., каф. информатики. Дипл. математик (БГУ, 1981). Канд. физ.-мат. наук по мех-ке жидк., газа и плазмы (защ. в ЛГУ, 1985). Иссл. в обл. теории крыла, моделир. поведения сл. систем.