МЕТОДИКА ОЦЕНКИ ВЕСОВЫХ ХАРАКТЕРИСТИК КОНСТРУКЦИИ КРЫЛЬЕВ ТРАНСПОРТНЫХ САМОЛЁТОВ С ИСПОЛЬЗОВАНИЕМ МНОГОДИСЦИПЛИНАРНОГО ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ И МЕТОДОВ МНОГОКРИТЕРИАЛЬНОЙ ОПТИМИЗАЦИИ

УДК 629.7.01

МЕТОДИКА ОЦЕНКИ ВЕСОВЫХ ХАРАКТЕРИСТИК КОНСТРУКЦИИ КРЫЛЬЕВ ТРАНСПОРТНЫХ САМОЛЁТОВ С ИСПОЛЬЗОВАНИЕМ МНОГОДИСЦИПЛИНАРНОГО ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ И МЕТОДОВ МНОГОКРИТЕРИАЛЬНОЙ ОПТИМИЗАЦИИ

© 2020 О.У. Эспиноса Барсенас, О.Е. Лукьянов

Самарский национальный исследовательский университет им. академика С.П. Королёва,

г. Самара, Россия

Статья поступила в редакцию 27.10.2020

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

Ключевые слова: конструкция крыла, весовые характеристики, коэффициент силового фактора, математическое моделирование, метод дискретных вихрей, конечно-элементное моделирование, конструктивно-силовая схема, многокритериальная оптимизация, геометрические параметры. DOI: 10.37313/1990-5378-2020-22-5-120-127

ВВЕДЕНИЕ

Задача выбора основных геометрических параметров схемы самолёта на начальных стадиях проектирования формулируется в терминах нелинейного математического программирования [1, 2, 3, 4], что предполагает минимизацию (максимизацию) комплексной целевой функции-свёртки при удовлетворении геометрическим и функциональным ограничениям. Часто в качестве таких свёрток принимают взлётную массу самолёта, коэффициент топливной эффективности [1, 3] или более общие критерии, как себестоимость перевозок или прибыль от транспортной операции [5]. Важнейшим слагаемым таких целевых функций является масса конструкции планера самолёта, оказывающая непосредственное влияние на её итоговое значение. В связи с этим точность оценки массы конструкции планера самолёта, а в частности, крыла, несёт ключевое значение для достоверного определения основных лётно-технических и технико-экономических показателей будущего летательного аппарата (ЛА) ещё на начальных стадиях проектирования. Традиционные методики проектирования самолётов, сложившиеся

Эспиноса Барсенас Оскар Улисес, магистрант института авиационной техники. E-mail: oscar.espinosa.barcenas@gmail.com Лукьянов Олег Евгеньевич, кандидат технических наук,, доцент кафедры конструкции и проектирования летательных аппаратов. E-mail: lukyanovoe@mail.ru

в ведущих конструкторских бюро за длительное время, используют так называемые весовые формулы для оценки массы конструкции планера самолёта, которые базируются на результатах обработки и обобщения обширного количества статистического материала в совокупности с простыми аналитическими соотношениями сопротивления материалов, связывающими вес материала элементов конструкции и величину внутреннего силового фактора при типовых видах силовой работы. Качественная систематизация статистических данных и их скрупулёзный анализ позволили выстроить надёжный инструмент расчёта массы существующих конструкций с возможностью отслеживания влияния на её величину различных параметров ЛА, его особенностей, применяемых материалов, условий эксплуатации, предъявляемых требований и прочих факторов. Работа с такими математическими моделями требует минимальных временных и вычислительных затрат. По сути это чётко определённые аналитические зависимости, к которым могут быть применены аналитические непоисковые методы нелинейного программирования для использования их в задачах оптимизации параметров облика самолёта. Разработкой весовых формул занимались практически во всех ведущих инженерных коллективах по всему миру. Подробно с такими формулами можно ознакомиться в [6, 7, 8, 9, 10]. В то же время следует отметить, что математические модели, базирующиеся на информации о

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

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

В работе [11] представлены научные основы нового метода определения массы конструкций, базирующегося на их силовой работе. Авторы [11] и [12, 13] использует для расчёта теоретической массы силовой конструкции специфический интегральный критерий весового совершенства «силовой фактор» О, характеризующий величину и протяжённость действия усилий в конструкции:

О =

V

где (Гэкв - эквивалентные напряжения в конструкции, Па, V - объём конструкции, м-5

Этот критерий хорошо подходит для поиска и генерации новых силовых схем как с использованием аналитических, так и численных математических моделей механики твёрдого тела. Важнейшим преимуществом такого подхода является возможность оценки весового совершенства именно самой силовой схемы (вне зависимости от её сложности, а также материала конструкции). Критерий «силовой фактор» одновременно учитывает величину и протяжённость действия нагрузки в конструкции, а потому позволяет определять минимально потребный (теоретический) объём и, соответственно, массу равнопрочной полнонапряжённой силовой конструкции, необходимой для восприятия заданного спектра нагрузок с заданным уровнем напряжений. Для проведения весовых оценок нет необходимости подбора сечений силовых элементов: важна лишь силовая схема. В [11, 12] введено также понятие обезразмеренного варианта критерия О - коэффициента силового фактора СК, позволяющего производить сравнение двух или более конкурирующих силовых схем вне зависимости от их абсолютных размеров и величин действующих на них нагрузок. Величина критерия О пропорциональна полной массе исследуемой конструкции тк [11].

В работах [12, 13, 14] приведены примеры эффективного использования критериев О и СК для

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

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

1. УЧЁТ ВЛИЯНИЯ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ КРЫЛА И СИЛОВОЙ СХЕМЫ

НА ЕГО ВЕСОВЫЕ ХАРАКТЕРИСТИКИ

На рисунке 1 представлена общая схема функционирования методики, которая предлагается для определения величины критерия СК конструкции крыла с учётом влияния его силовой схемы и внешних геометрических параметров.

Предполагается, что в качестве исходных данных выступают следующие величины:

- геометрические: площадь крыла, удлинение, сужение, стреловидность, относительная толщина профилей, геометрическая крутка, прочие;

- полётные: высота полета и соответствующие ей параметры атмосферы, крейсерская и максимальная скорости полета, максимальная взлетная масса, максимальная эксплуатационная перегрузка.

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

После ввода исходных данных выполняется выбор силовой схемы крыла с определением количества и положения лонжеронов в крыле на основе

Оптимизация положения

Исходные данные лонжеронов

н ограничения

Г Я, ij, т0,Р0'Кф> | —* чКщвЖй:,;' . .......

1пэ  пэ -1 V y miw у тах • ••)

Рис. 1. Принципиальная схема функционирования методики определения весовых характеристик крыла с учётом его геометрических параметров с использованием математического моделирования

многокритериальном оптимизации, далее проводится математическое моделирование обтекания крыла на заданном режиме полёта для точного определения величины и характера распределения воздушной нагрузки по площади крыла. Затем выполняется расчёт напряжённого состояния конечно-элементной модели крыла под полученные воздушные нагрузки для расчёта критерия весового совершенства О и его безразмерного коэффициента СК, который выбран в данной работе в качестве эквивалента массы силовой конструкции крыла. Критерий СК, как и О, пропорционален полной массе конструкции Ск гс шк [11].

2. МАТЕМАТИЧЕСКИЕ МОДЕЛИ АЭРОДИНАМИКИ И МЕХАНИКИ ТВЁРДОГО ТЕЛА

В качестве математических моделей аэродинамики используется метод дискретных вихрей (МДВ) [15] в линейной стационарной постановке. Данный метод в выбранной постановке использует стационарное несжимаемое потенциальное течение, что, однако не позволяет производить учёт сил вязкости, но даёт возможность с высокой точностью вычислять распределение давления по площади крыла при минимальных вычислительных и временных затратах. Учёт силы вязкости и поправка на сжимаемость среды выполняются с помощью формул инженерной методики расчёта аэродинамических характеристик. Предполагается, что методика рассчитана на самолёты транспортной категории, которые оснащаются крыльями большого удлинения, обладающими в свою очередь линейной зависимостью коэффициента подъёмной силы от угла атаки в диапазоне углов атаки с безотрывным обтеканием.

Математическое моделирование по МДВ предполагает вычисление циркуляций в контрольных точках вихревой системы тонкой несущей поверхности (рис. 2 а), геометрически совпадающей со срединной поверхностью крыла, на основе численного решения системы дифференциальных уравнений, отражающих закон Био-Савара-Лапласа: [15, 16].

А = Г-У ,

где А - матрица аэродинамического влияния, 1/м, Г - матрица-столбец циркуляции вектора скорости, м2/с; У - матрица-столбец проекций скорости набегающего потока на нормаль к вихревым рамкам в контрольных точках, м/с.

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

ут = К,- п,

где Уп. - проекция скорости набегающего потока на г-ю нормаль; м/с, Ух - вектор набегающего потока, м/с; Ц - нормаль к г-й вихревой рамке; г - номер вихревой рамки.

На основе полученного распределения цир-куляций рассчитывается элемент нормальной силы в каждой контрольной точке с использованием теоремы Жуковского [15]. На рисунке 2 а представлен общий вид вихревой модели для реализации МДВ.

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

c = g

ya о

у тах

Р О

q

где суа - коэффициент подъёмной силы, потребный для обеспечения равенства подъёмной силы весу самолёта с учётом перегрузки, п у шах - максимальная эксплуатационная перегрузка, р0 - удельная нагрузка на крыло, qmax - максимальный скоростной напор.

При рассмотрении расчётного случая по минимальной эксплуатационной перегрузке в соотношении (4) множитель пэу шах меняется на п эу т!л.

Затем, путём расчёта по МДВ, выполняется поиск угла атаки а, на котором достигается суа, потребный для уравновешивания самолёта. Расчёт распределения воздушной нагрузки по крылу выполняется по МДВ для вычисленного угла атаки а.

Расчёт напряжённого состояния конструкции крыла и вычисление критерия G, выполняется методом конечно-элементного моделирования в линейной стационарной постановке, для чего автоматически выполняется построение расчётной сетки из оболочечных конечных элементов, воспринимающих растяжение/сжатие, сдвиг и изгиб. Расчёт силового фактора (1) здесь выполняется численно с конечно-элементной модели (КЭМ) (рис. 2 б) по соотношению:

О = I

а

'Г; ,

где i - номер элемента модели.

Полученный размерный критерий G обе-зразмеривается до коэффициента СК путём нормирования на характерную нагрузку (вес самолёта т^) и характерный линейный размер (корень квадратный из площади крыла 5):

Си =—.

то

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

3. ПЕРЕДАЧА НАГРУЗОК ИЗ ВИХРЕВОЙ МОДЕЛИ В КЭМ

Отдельной задачей является обеспечение корректного приложения к КЭМ воздушных нагрузок, которые получены аэродинамическим расчётом: размер и взаимное расположение вихревых рамок и конечных элементов в КЭМ могут не совпадать. Кроме того, КЭМ включает в себя нервюры и лонжероны, в совокупности представляющие собой систему жёстких в своей плоскости элементов - стенок, работающих на сдвиг. Однако, расположение некоторых контрольных точек вихревой модели может не совпадать с осью того или иного силового элемента на виде в плане. В таком случае нагрузка, полученная в такой контрольной точке, будет приходиться на нежёсткую из своей плоскости обшивку, что может привести возникновению локальных изгибающих моментов в ней, существенному повышению напряжений и значения силового фактора в этих местах и невозможности адекватно оценить массу конструкции крыла. Для преодоления этих проблем, в работе предложено соблюдать размеры ячеек двух моделей кратными друг другу, например, так, чтобы в одну вихревую рамку на виде в плане вписывались два или четыре элемента КЭМ при совпадении общих узлов. Сила, полученная в контрольной точке вихревой рамки, переносится вдоль линии действия напрямую в соответствующий узел КЭМ, расположенный на верхней поверхности крыла. Кроме того, при построении вихревой модели и КЭМ предварительно высчитывается количество разбиений несущей поверхности на вихревые рамки таким образом, чтобы контрольные точки совпадали с осями силовых элементов на виде в плане и приложение нагрузки приходилось непосредственно на них.

Рис. 2. Общий вид модели крыла: а) вихревая модель МДВ; б) конечно-элементная модель

4. ПОИСК ОПТИМАЛЬНОГО РАСПОЛОЖЕНИЯ ЛОНЖЕРОНОВ В КРЫЛЕ ПО ХОРДЕ

Задача выбора размещения лонжеронов в крыле по хорде формулируется в терминах нелинейного математического программирования с двумя критериями. В качестве критериев выступают строительная высота лонжеронов с одной стороны и межлонжеронная площадь поперечного сечения крыла - с другой. Строительная высота лонжерона зависит от его положения по хорде, определяется геометрией профиля крыла и подлежит максимизации с целью повышения его изгибной жёсткости. Площадь поперечного сечения крыла, заключённая между лонжеронами, представляет собой замкнутый контур, воспринимающий действие крутящего момента. Она подлежит максимизации для повышения крутильной жёсткости. Проектные переменные - координаты положения лонжеронов по хорде (рис. 3).

Рис. 3. Интервалы варьирования проектными переменными

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

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

У1( х,, х2) —

( Н3( х,) + Н 23( х2)у1

2

^ тт,

1

Л(^х2) = _ _ , ^ т1п

а>

X . < х, < х,

1шт 1 1шах

2шт 2 2шах

ч

где а>(х1, х2) = | (gв ^н )<Яс - площадь меж-

лонжеронной части поперечного сечения крыла; Н1

- строительная высота 2-го лонжерона; gg и gН

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

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

Решение поставленной задачи оптимизации предложено выполнять с использованием многоцелевого эволюционного алгоритма ЫБОЛ-П [17], относящегося к группе генетических алгоритмов.

5. РЕАЛИЗАЦИЯ АЛГОРИТМА

Изложенная методика алгоритмизирована и функционирует полностью в автоматическом режиме в виде программного обеспечения для ЭВМ на межплатформенной связке двух программных продуктов МЛТЬЛБ и Лшуз МесЬатса1 ЛРБЬ. При этом модуль расчёта аэродинамических характеристик по МДВ, модуль оптимизации положения лонжеронов в крыле, модуль автоматической генерации расчётных сеток для МДВ и КЭМ, а также общая управляющая программа реализованы на языке МЛТЬЛБ. Расчёт напряжённого состояния КЭМ и вычисление критерия «силовой фактор» - на языке ЛРБЬ. Блок-схема программы представлена на рисунке 4.

6. РЕШЕНИЕ ДЕМОНСТРАЦИОННОЙ ЗАДАЧИ

С целью апробации алгоритма и программного обеспечения выполнено решение прикладной задачи по оценке влияния геометрических параметров крыла гипотетического тяжёлого транспортного самолёта на коэффициент силового фактора его конструкции. Исходные данные: взлётная масса т0=350 т; высота крейсерского полёта Н=10000 м; крейсерская скорость V = 850 км/ч; максимальная эксплуатационная перегрузка пэушах = 2,5; минимальная эксплуатационная перегрузка п эу т1п= -1,5 удельная нагрузка на крыло р0 = 650 ДаН/м2. Геометрические параметры, принятые переменными: стреловидность х, сужение постоянные параметры: относительная толщина корневого профиля крыла с0=0,12, концевого - ск=0,08; удлинение X = 8,3; геометрическая крутка ф = 0°.

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

Из множества сочетаний для двухлонжерон-ного крыла выбраны относительные координаты положения лонжеронов х1 = 0,23 и х2 = 0,67. Данное сочетание проектных переменных обеспечивает оптимальность по Парето [18].

Рис. 4. Блок-схема алгоритма вычисления критерия СК

00

¿До & *

о

б)

0,2 0,22 0,24 0,26 Х1

Рис. 5. Кривые Парето: а) значения целевых функций; б) значения проектных переменных

На рисунках 6 а, б представлены распределения относительной циркуляции по размаху крыла при различных сужениях И и стреловид-

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

На рисунке 7 показана зависимость коэффициента силового фактора СК от этих двух переменных параметров - сужения И и стреловидности с, полученная на основе конечно-элементного моделирования под соответствующие распределения циркуляций.

Зависимости, представленные на рисунках 6 - 7, имеют логичный и ожидаемый качественный и количественный характер. Увеличение стреловидности повышает массу крыла из-за увеличения циркуляции в концевых сечениях крыла и уменьшения - в корневых (рис. 6 б). Увеличение сужения - снижает массу за счёт перераспределения воздушной нагрузки ближе к корневым сечениям и снижения изгибающего момента (рис. 6 а). Значения коэффициента силового фактора при соответствующих параметрах количественно и качественно совпадают с результатами, полученными в работах [19, 20], которые посвящены исследованию возможностей применения этого критерия для расчёта массы конструкции крыла транспортных самолётов.

Для получения одной точки на графике (рис.

7) по двум расчётным случаям (пуутах и пуутт

а)

0,4 0,6 0,8 Z 0,0

-3 ----4---5----6 б)

0,2 0,4 0,6 0,8 z -20° ----30----40°

Рис. 6. Распределение относительной циркуляции по размаху крыла: а) при разных сужениях крыла

Рис. 7. Зависимость критерия СК конструкции крыла в зависимости от сужения ц и стреловидности х

) потребовалось 25 секунд машинного времени на одном ядре в два потока на процессоре Core i5 с частотой 2,3 ГГц.

ЗАКЛЮЧЕНИЕ

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

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

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

1. Комаров, В.А. Многодисциплинарная оптимизация параметров крыла грузового самолета / В. А. Комаров, О.Е. Лукьянов // Общероссийский научно-технический журнал „Полет". - 2018. - № 3. -С.3-15.

2. Wunderlich, Tobias. Multidisziplinärer Entwurf und Optimierung von Flügeln für Verkehrsflugzeuge// Deutscher Luft- und Raumfahrtkongress 2009. 08. -

10. Sept. 2009. pp. 1-10

3. Кузнецов А.С. Алгоритм выбора рациональных параметров крыла с учетом аэродинамической и весовой эффективности. Известия Самарского научного центра Российской академии наук, т. 12, №1(2), 2010.

4. Ciampa P.D., Nagel B. AGILE Paradigm: the next generation collaborative MDO for the development of aeronautical systems // Progress in Aerospace Sciences. 2020 №119. С. 1-20

5. Арутюнов А.Г. Методика определения рационального облика коммерческого тяжелого рампо-вого грузового самолета на этапе концептуального проектирования из условия его прибыльности: Дисс. канд. техн. наук. Москва, 2017. 252 с.

6. Бадягин А.А., Егер С.М. и др. Проектирование самолетов. М.: Машиностроение, 1972. — 516 с.

7. Torenbeek E. Synthesis of Subsonic Airplane Design, Delft University Press 1983.

8. Raymer Daniel P. Aircraft design: A Conceptual Approach. Washington: Reston, VA : American Institute of Aeronautics and Astronautics, 5th ed., 2012. P. 227.

9. Roskam Jan. Airplane Design Part I : Preliminary Sizing of Airplanes. Design, Analisys and Research Corporation, Lawrence Kansas . 5-th printing, 2017. 206 p.

10. Шейнин В.М., Козловский В.И. Весовое проектирование и эффективность пассажирских самолетов. Т. 1. Весовой расчет самолета и весовое планирование. М., «Машиностроение», 1977, 344 с.;

11. Комаров В.А. Весовой анализ авиационных конструкций: теоретические основы// Общероссийский научно-технический журнал «Полёт». 2000. №1. С.31-39.

12. КомаровВ.А. Гуменюк А.В. Оценка весовой эффективности силовых схем несущих поверхностей // Вестник Самарского государственного аэрокос-

мического университета им. академика С.П. Королёва (национального исследовательского университета). 2002. № 1. С. 45-54.

13. Болдырев А.В., Комаров В.А. Использование высокоточного моделирования на ранних стадиях проектирования конструкций летательных аппаратов // Вестник Самарского государственного аэрокосмического университета им. академика С.П. Королёва (национального исследовательского университета). 2012. № 1(32). С 15-23.

14. Резниченко Г.А. Вычисление силового веса конструкций типа крыла с использованием метода конечных элементов. // Автоматизация проектирования авиационных конструкций. Межвузовский сборник. КуАИ: Куйбышев, 1979. С.77-83

15. Белоцерковский С.М. Тонкая несущая поверхность в дозвуковом потоке газа. Под ред. Н.И.Розальской, И.Ш. Аксельрод, Т.Д. Доверман. - М.: Наука, 1965. - 244 с.

16. Katz J., Plotkin A. Low-Speed Aerodynamics. М.: McGraw-Hill, Inc., 2001. 629с

17. Deb K., Agrawal S., Pratap A. A Fast Elitist Non-dominated Sorting Genetic Algorithm for Multi-objective Optimization: NSGA-II // International conference on parallel problem solving from nature. 2000. С. 849-858.

18. Ногин В.Д. Принятие решений в многокритериальной среде: количественный подход. - М.: Физматлит, 2002, 2005, 176 с.

19. Кузнецов А.С. Выбор геометрических параметров крыла с комплексным учётом аэродинамической и весовой эффективности: Дисс. канд. техн. наук. Самара 2011 г. 161 с.

20. Лаптева М.Ю. Разработка методики прогнозирования и учёта деформаций крыла на ранних стадиях проектирования с использованием модели тела переменной плотности: дисс. ... канд. техн. наук. Самара, 2012. 118 с.

METHODOLOGY FOR ASSESSMENT OF THE WEIGHT CHARACTERISTICS OF THE CARGO AIRCRAFT WINGS USING MULTIDISCIPLINARY NUMERICAL SIMULATION AND ANALYTICAL METHODS OF MULTI OBJECTIVE OPTIMIZATION

© 2020 O.U. Espinosa Barcenas, O.E. Lukyanov

Samara National Research University, Samara, Russia

The developed methodology determines the weight characteristics of a wing of high aspect ratio using the integral criterion - the carrying load factor coefficient and numerical mathematical modeling in the fields of aerodynamics and rigid body mechanics. The paper considers the application of a multi-objective optimization method based on genetic algorithms to find the rational position of the spars in the wing when choosing its structural scheme. The calculation of the air load distribution within the framework of the proposed method is based on mathematical modeling of flows by the vortex lattice method. The ideal weight criterion is calculated based on the stress state of the wing structure obtained by the finite element method. The efficiency of the method is shown by the example of solving a demonstration problem to assess the influence of several geometric parameters of the wing on the air load distribution and its weight characteristics.

Keywords: wing structure, weight characteristics, carrying load factor coefficient, mathematical modeling, vortex lattice method, finite element modeling, airframe, multi-objective optimization, geometric parameters.

DOI: 10.37313/1990-5378-2020-22-5-120-127 Oscar Ulises Espinosa Barcenas, Master's Student at Institute of Aeronautical Engineering. E-mail: oscar.espinosa.barcenas@gmail.com Oleg Lukyanov, Candidate Degree in Engineering, Associate Professor at Aircraft Construction and Design Department. E-mail: lukyanovoe@mail.ru