CC BY
27ПРОБЛЕМЫ КОРАБЛЕСТРОЕНИЯ И ОКЕАНОТЕХНИКИ
УДК 629.572
А.А. Болотин
ПРИМЕНЕНИЕ МЕТОДА ДИСКРЕТНЫХ ВИХРЕЙ ДЛЯ ИССЛЕДОВАНИЯ ПОДВОДНЫХ КРЫЛЬЕВ
Нижегородский государственный технический университет им. Р.Е. Алексеева
Применение математического моделирования - необходимое условие для успешного проектирования судов на подводных крыльях в настоящее время. Представлено описание методики расчета гидродинамических характеристик подводных крыльев при помощи метода дискретных вихрей и приведены результаты расчетов.
Ключевые слова: подводные крылья, метод дискретных вихрей, гидродинамические характеристики, подъемная сила, лобовое сопротивление.
Идея поставить судно на крылья, безусловно, стала революционной в скоростном судостроении. Сочетание скорости, комфорта, экономичности в свое время сделали суда на подводных крыльях (СПК) очень популярным видом водного транспорта как у нас в стране, так и за рубежом. Долгие годы эксплуатации подтвердили эффективность этих судов, и возросший в последнее время интерес к СПК совершенно закономерен. Новый этап в развитии крылатых судов требует и новых методов в проектировании, в частности, более широкого внедрения методов математического моделирования.
Задача создания гидродинамического комплекса является основной при проектировании СПК. Выбор гидродинамической схемы и доведение ее до требуемых показателей - это результат многолетнего труда опытнейших специалистов. В основе этой работы в течение многих лет лежал модельный эксперимент, в частности, испытания буксируемых моделей на открытой воде. Использование математического моделирования, несомненно, способно во много раз сократить сроки и стоимость процесса создания подводных крыльев.
Первым этапом этой работы должно являться расчетное определение гидродинамических характеристик крыльев, как необходимых исходных данных для исследования динамики движения СПК, устойчивости и управляемости. Значения гидродинамических нагрузок и их распределение вдоль крыла необходимы для расчетов прочности.
В данной работе представлены результаты расчета гидродинамических характеристик подводных крыльев методом дискретных вихрей [1]. Вихревые методы расчета несущих поверхностей достаточно прочно вошли в процесс проектирования летательных аппаратов, но не очень активно используются для исследования подводных крыльев. Тем не менее, крыльевая система СПК является очень удобным объектом исследования. Тонкие несущие поверхности трапециевидной формы, отсутствие объемных элементов позволяют уверенно применять вихревые методы для расчета гидродинамических характеристик. Метод дискретных вихрей отличается относительной простотой и наглядностью, а также легко распространяется на сложные крыльевые системы, какими являются подводные крылья. Некоторая некорректность при определении распределения нагрузки вдоль профиля для подводного крыла является неактуальной, так как это сказывается на точности определения продольного мо-
© Болотин А.А., 2015.
мента крыла относительно его передней (или задней) кромки, а его величина ввиду большого расстояния крыла от центра тяжести судна незначительна.
В качестве объекта исследования выбрано носовое крыльевое устройство СПК «Комета», являющееся основой для крыльев большинства морских СПК. Схема крыла представлена на рис. 1. Профиль крыла - плоско-выпуклый сегмент с максимальной толщиной 5,5% .
Вид спереди
Вид сверху
Рис. 1. Схема крыльевого устройства
Метод дискретных вихрей предполагает разбиение крыла на площадки, размещение на площадке вихря с неизвестной циркуляцией, определение этих циркуляций и гидродинамических сил, действующих на каждую площадку, а затем определение сил, действующих на все крыло. Для моделирования рассмотренное крыльевое устройство разбито на 13 трапециевидных участков. Нумерация участков представлена на схеме. Профилированные участки крыла (с 1 по 8) заменяются криволинейными несущими поверхностями, расположенными вдоль средней линии профиля, стойки крыла - плоскими несущими поверхностями. Криволинейные несущие поверхности разбиваются на 6 частей по хорде, плоские на 5. Разбиение по размаху определяется удлинением рассматриваемого участка и меняется от 8 до 1. Граница раздела моделируется зеркальной вихревой системой с циркуляцией того же знака, что и у основной вихревой системы. Задача решается в линейной постановке, свободные вихри лежат в плоскости хорд крыла. Результатами расчета являются гидродинамические коэффициенты трех проекций гидродинамической силы: лобового сопротивления (С), подъемной
силы ( су ) и боковой силы ( С ), а также три проекции коэффициентов гидродинамического
момента (тх, т , тг). Коэффициенты связаны с гидродинамическими силами
соотношениями:
с г = г = = 2Ых = 2Му =_2М_
С р¥2Я ' Су р¥2Я ' С2 рУ2Я' тх рУ28Ь ' ту рУ28Ь ' т рУ28Ь '
где Ях, Я , Я2 - проекции гидродинамической силы;
Мх, Му, М2 - проекции гидродинамического момента;
V - скорость судна; р - плотность воды; £ - характерная площадь; Ь - характерный линейный размер.
Из составляющих лобового сопротивления в рамках представленного метода определяется только индуктивное сопротивление крыльевого устройства. Для определения полного
гидродинамического сопротивления необходимо определить составляющие вязкостного Сп
и волнового сопротивления Схвол.
Согласно [2], вязкостное сопротивление :
С_ = 2С..
1+(и -
где СШл=0А551(^ Яв)
2,58.
Яв=
УЬср/V;
и = 1 - е
-2(к ср)°
ср)°
Ки = 1 + (0,5 + с )е
г Ьр к ср = —,
Ьср
где с - относительная толщина профиля; V - кинематическая вязкость воздуха; Ьср - средняя хорда крыла; кср - среднее погружение крыла. Волновое сопротивление
С V 2 с =- у
Схвол
2Рг2
-2кс
, Гг1
1 -
2ж
Р7
-2кс
, Рг2
где Рг =
V
•М
- число Фруда.
ср
На рис. 2 приведены результаты расчета коэффициента подъемной силы в зависимости от угла атаки а и погружения крыльевого устройства. За величину погружения к принято отстояние задней кромки крыла в диаметральной плоскости от поверхности воды в долях хорды. В качестве характерной площади принята горизонтальная проекция всей несущей части крыльевого устройства.
На этом же рисунке приведены результаты модельных испытаний данного крыльевого устройства, полученные в опытовом бассейне ЦНИИ им акад. Крылова [3]. Испытания проводились на модели масштабом 1:10 при скорости буксировки 5,5 м/с.
Рис 2. Расчетные и экспериментальные значения коэффициента подъемной силы от угла атаки
С целью идентификации расчетных и экспериментальных данных расчет проводился при тех же числах Рейнольдса, что и эксперимент. Расчетные углы атаки и погружения выбраны такими же, как в эксперименте. Результаты расчета представлены на рисунке сплошными линиями. Экспериментальные данные показаны отдельными точками.
Представленные результаты свидетельствуют об удовлетворительной сходимости результатов расчета и испытаний, что подтверждает целесообразность использования метода дискретных вихрей для расчета подводных крыльев.
Несомненным достоинством предлагаемой методики является пространственная постановка задачи. Расчет гидродинамических характеристик при крене и дрейфе не требует введения дополнительных условий. В качестве примера на рис. 3, 4 представлены зависимости коэффициента кренящего момента тх рассмотренного крыльевого устройства от угла
крена у и от угла дрейфа в при различных погружениях. В качестве характерного линейного размера принят размах крыла. Зависимости, представленные на рис. 3, отслеживают вход в воду и выход из воды элементов крыльевого устройства с различными углами К-образности.
Рис. 3. Зависимость коэффициента кренящего момента от угла крена
Рис. 4. Зависимость коэффициента кренящего момента от угла дрейфа
Проведенные методические исследования свидетельствуют о высокой степени достоверности расчета боковых гидродинамических характеристик.
Для расчета гидродинамических характеристик с отклоненными органами управления достаточно изменить форму несущей поверхности в соответствии с углом отклонения управляющего элемента. При этом необходимо часть несущей поверхности, обслуживаемую органом управления, выделить отдельным участком, и разбиение по хорде провести таким образом, чтобы отклоненная часть несущей поверхности моделировалась отдельными площадками. Проведенные исследования показали, что расчетные результаты хорошо согласуются с экспериментальными в области углов отклонения органов управления, не превышающих 100.
Как указывалось ранее, получение гидродинамических характеристик производится путем интегрирования нагрузки, распределенной по площади крыла. Определение значений нагрузок, распределенных вдоль размаха крыла и вдоль стоек крыла, является отдельной задачей, достаточно эффективно решаемой при помощи данной методики.
Таким образом, предложенная методика позволяет с достаточной точностью получить полный спектр гидродинамических характеристик крыльевых устройств СПК, что значительно упрощает процесс проектирования судов. Полученные результаты являются исходными данными для математического моделирования динамики движения СПК, для решения задач получения кривой сопротивления [4], устойчивости, остойчивости и управляемости судна, а также получить исходные данные для расчета прочности.
Библиографический список
1. Белоцерковский, С.М. Тонкая несущая поверхность в дозвуковом потоке газа / С.М. Бело-церковский. - М.: Наука, 1965. - 244 с.
2. Егоров, И.Т. Гидродинамика быстроходных судов / И.Т. Егоров, В.Т. Соколов. - Л.: Судостроение, 1971. - 424 с.
3. Отчет по испытаниям в гидроканале моделей крыльевых устройств с поворотным закрылком. Технический отчет, ЦКБ по СПК, инв. № ОР 224-020-009, 1968.
4. Болотин, А.А. Методика расчета кривой сопротивления и посадки судна на подводных крыльях / А.А. Болотин // Проблемы прочности и пластичности: мевжуз. сб. / Нижегород. ун-т. - Н. Новгород, 2014. Вып. 76. С. 172-177.
Дата поступления в редакцию 02.06.2015
A.A. Bolotin
THE USE OF DISCRETE VORTICES METHOD FOR THE HYDROFOILS RESEARCH
Nizhny Novgorod state technical university n.a. R.E. Alexeev
The aim of this work is the development of methodology for hydrofoil hydrodynamic characteristics calculation. The method of discrete vortices, widely used in the aircraft industry, was chosen for calculation. The features of this method associated with its usage for hydrofoil calculation were distinguished. The hydrodynamic characteristics calculation for the fore wing device of the hydrofoil craft "Kometa" were done. The comparison of the calculation results and the results of towing tests of this wing device in Krylov Shipbuilding Research Institute isdemonstrated. The comparison concluded that theusageof this method for the hydrofoils calculations is feasible.The proposed method allows obtaining the raw data for mathematical modeling of the dynamics of hydrofoil craft motion, for the calculation of theresistance, landing, stability, controllability and strengthcurves.
Key words: hydrofoil, method of discrete vortices, hydrodynamic characteristics, dynamic lift, drag.
