Известия РАН. Теория и системы управления, 2023, № 6, стр. 166-176

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

К. С. Анисимов a, В. Н. Евдокименков b*, М. Н. Красильщиков b, К. И. Сыпало a, Н. Б. Топоров b

a ЦАГИ им. проф. Н.Е. Жуковского
Москва, Россия

b Московский авиационный институт (национальный исследовательский ун-т)
Москва, Россия

* E-mail: evdokimenkovvn@mai.ru

Поступила в редакцию 08.06.2023
После доработки 19.06.2023
Принята к публикации 31.07.2023

Полный текст (PDF)

Аннотация

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

Введение. Государственной программой Российской Федерации “Развитие авиационной промышленности на 2013–2025 годы” предусматривается повышение конкурентоспособности российской авиационной промышленности на мировом рынке авиационной техники и услуг. В свою очередь существующий уровень конкуренции на этом рынке требует повышения качества проектирования перспективных авиационных комплексов на основе современных информационных технологий.

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

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

Как полагают авторы статьи, этап функционального проектирования включает решение двух основных задач:

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

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

Метод решения первой из сформулированных задач подробно описан в [1, 2]. Показано, что эта задача может быть интерпретирована как классическая задача математического программирования, в которой множество допустимых значений летно-технических характеристик (ЛТХ), участвующих в расчете интегральных показателей, отражающих многодисциплинарные требования, предъявляемые к перспективному ЛА, формируется на основе накопленного опыта проектирования.

Как уже указывалось выше, задача оценки дефицита функциональных свойств и определения приоритетных научно-технических решений возникает в случае невозможности создания перспективного ЛА, удовлетворяющего всей совокупности многодисциплинарных требований, на современном уровне развития авиационных технологий. В формальной постановке эта задача определяется как обобщенная задача оптимизации [35], отличие которой от классической задачи математического программирования заключается в дополнительной операции оптимизации множества допустимых значений ЛТХ, поскольку влияние перспективных научно-технических решений проявляется, прежде всего, в изменениях формы и размеров этого множества, и как результат само множество становится дополнительным объектом оптимизации.

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

1. Постановка задачи оценки дефицита функциональных свойств перспективных ЛА. В рамках развиваемой концепции функционального проектирования [1, 2] приемлемым полагается облик перспективного ЛА, описываемый вектором проектных параметров r, при котором обеспечивается одновременное выполнение комплекса многодисциплинарных требований.

Формальный учет этих требований базируется на использовании интегральных показателей ${{I}_{j}}\left( r \right),j = \overline {1,N} $, количественно выражающих степень соответствия проектируемого ЛА заявленным требованиям:

(1.1)
$\left\{ \begin{gathered} {{I}_{{1~{\text{min}}}}} \leqslant {{I}_{1}}\left( r \right) \leqslant {{I}_{{1~{\text{max}}}}}, \hfill \\ {{I}_{{2~{\text{min}}}}} \leqslant {{I}_{2}}\left( r \right) \leqslant {{I}_{{2~{\text{max}}}}}, \hfill \\ {{I}_{{3~{\text{min}}}}} \leqslant {{I}_{3}}\left( r \right) \leqslant {{I}_{{3~{\text{max}}}}}, \hfill \\ \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .. \hfill \\ {{I}_{{N~{\text{min}}}}} \leqslant {{I}_{N}}\left( r \right) \leqslant {{I}_{{N~{\text{max}}}}}. \hfill \\ \end{gathered} \right.$

Предельные значения показателей ${{I}_{{j~{\text{min}}}}},{{I}_{{j~{\text{max}}}}}$ устанавливает специалист (эксперт), исходя из заданных требований к перспективному ЛА. Заметим, что система неравенств (1.1) позволяет в наиболее общем виде учесть требования к базовым функциям перспективного ЛА.

В дальнейшем комплекс показателей (1.1) будем представлять в виде вектор-функции $I\left( r \right) = {{\left( {{{I}_{1}}\left( r \right),{{I}_{2}}\left( r \right), \ldots ,{{I}_{N}}\left( r \right)} \right)}^{{\text{T}}}}$, для которого заданным полагается множество допустимых значений WI:

(1.2)
${{W}_{I}} = \left\{ {I\left( r \right):{{I}_{{j~{\text{min}}}}} \leqslant {{I}_{j}}\left( r \right) \leqslant {{I}_{{j~{\text{max}}}}},j = \overline {1,N} } \right\}.$

Компоненты вектора $r = {{\left( {{{r}_{1}},{{r}_{2}}, \ldots ,{{r}_{n}}} \right)}^{{\text{T}}}}$ в (1.1) в совокупности определяют ЛТХ, участвующие в вычислении интегральных показателей ${{I}_{j}}\left( r \right),j = \overline {1,N} $. В дальнейшем предполагается, что задано множество допустимых значений ЛТХ как n-мерный параллелепипед Wr определяемый совокупностью условий:

(1.3)
${{W}_{r}} = \left\{ {r:{{r}_{{i~{\text{min}}}}} \leqslant {{r}_{i}} \leqslant {{r}_{{i~{\text{max}}}}},i = \overline {1,n} } \right\}.$

Предельные значения ЛТХ ${{r}_{{i~{\text{min}}}}}$,$~{{r}_{{i~{\text{max}}}}}$, $i = \overline {1,n} $, в (1.3) задаются с учетом накопленного опыта проектирования ЛА, который базируется на доступных для использования научно-технических и технологических решениях.

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

В подобной ситуации использование перспективных научно-технологических решений, внедрение которых возможно в обозримой перспективе, является очевидным способом выполнения многодисциплинарных требований (1.1). Как уже говорилось, влияние этих перспективных решений выражается, прежде всего, в изменении размеров множества допустимых значений ЛТХ Wr .

Следовательно, задача устранения дефицита функциональных свойств проектируемого ЛА заключается в минимально возможном расширении размеров Wr с целью определения нового множества $W_{r}^{*}$, в пределах которого обеспечивается одновременное выполнение комплекса условий (1.1).

В результате построения нового множества $W_{r}^{*}$ будут определены приоритетные научно-технологические решения, потенциально обеспечивающие возможность расширения диапазонов предельных значений отдельных ЛТХ, для которых будет подтверждена такая необходимость. При этом формирование множества $W_{r}^{*}$ должно осуществляться с учетом временных и финансовых затрат, технологических сложностей и рисков, сопровождающих развитие любого перспективного научно-технологического решения и таким образом существенно влияющих на факт его выбора.

2. Оценка дефицита функциональных свойств перспективного ЛА как обобщенная задача оптимизации. Поскольку интегральные показатели ${{I}_{j}}\left( r \right),j = \overline {1,N} $, используемые для оценки многодисциплинарных требований (1.1), имеют разный физический смысл и выражены в разных единицах измерений, необходим переход к их безразмерным нормированным аналогам ${{\hat {I}}_{j}}(r),j = \overline {1,N} $, на основе линейного преобразования:

(2.1)
$~{{\hat {I}}_{j}}\left( r \right) = ~\frac{{{{I}_{j}}\left( r \right) - 0.5\left( {{{I}_{{j{\text{min}}}}} + {{I}_{{j{\text{max}}}}}} \right)}}{{0.5\left( {{{I}_{{j{\text{max}}}}} - {{I}_{{j{\text{min}}}}}} \right)}}.$

Это нормирующее преобразование отображает N-мерный параллелепипед WI в куб ${{\hat {W}}_{I}}$, задаваемый условиями:

(2.2)
${{\hat {W}}_{I}} = \left\{ {\hat {I}\left( r \right): - 1 \leqslant {{{\hat {I}}}_{j}}\left( r \right) \leqslant 1,j = \overline {1,N} } \right\}.$

Введем в рассмотрение скалярные функции ${{{{\varphi }}}_{j}}\left( r \right),j = \overline {1,N} $, количественно выражающие для любого фиксированного вектора ЛТХ r факт выполнения каждого из условий (1.1). Эти функции определим следующим образом:

(2.3)
${{{{\varphi }}}_{j}}\left( r \right) = \left\{ \begin{gathered} - ({{{\hat {I}}}_{j}}\left( r \right) + 1),~\quad {{{\hat {I}}}_{j}}\left( r \right) < - 1, \hfill \\ 0,\quad - 1 \leqslant {{{\hat {I}}}_{j}}\left( r \right) \leqslant 1, \hfill \\ {{{\hat {I}}}_{j}}\left( r \right) - 1,\quad {{{\hat {I}}}_{j}}\left( r \right) > 1. \hfill \\ \end{gathered} \right.$

Далее используем скалярную функцию Φ(r), количественно выражающую для любого фиксированного вектора ЛТХ r штраф за невыполнение совокупности условий (2.2). Эту функцию определим таким образом, чтобы признаком выполнения всего комплекса условий (2.2) было значение ${{\Phi }}\left( r \right) = 0{\kern 1pt} {\kern 1pt} {\text{:}}$

(2.4)
$\left\{ {\begin{array}{*{20}{c}} {r:\forall j~{{\varphi }_{j}}(r) = 0 \Leftrightarrow \Phi (r) = 0,} \\ {r:\exists j~{{\varphi }_{j}}(r) \ne 0 \Leftrightarrow \Phi (r) > 0.} \end{array}} \right.$

Примерами функций ${{\Phi }}\left( r \right)$, удовлетворяющих условиям (2.3), (2.4), могут быть следующие:

(2.5)
${{\Phi }}\left( r \right) = \mathop \sum \limits_{j = 1}^N {{{{\varphi }}}_{j}}\left( r \right),$
(2.6)
${{\Phi }}\left( r \right) = \mathop \sum \limits_{j = 1}^N {{\varphi }}_{j}^{2}\left( r \right),$
(2.7)
${{\Phi }}\left( r \right) = \mathop {\max }\limits_j {{{{\varphi }}}_{j}}\left( r \right).$

Выбор конкретного вида оптимизируемой функции ${{\Phi }}\left( r \right)$ из набора (2.5)–(2.7) осуществляется с учетом используемых численных методов оптимизации и оценки затрат машинного времени для проведения вычислений.

Из (2.3), (2.4) следует, что вектор ЛТХ r*, для которого обеспечивается вся совокупность условий (1.1), является решением оптимизационной задачи:

(2.8)
$r{\text{*}} = {\text{arg}}\mathop {\min }\limits_{r \in {{W}_{r}}} {{\Phi }}\left( r \right).$

В результате решения задачи (2.8) возможны следующие ситуации:

Вариант 1: ${{\Phi }}\left( {r{\text{*}}} \right) = 0.~$ В этом случае дефицит функциональных свойств проектируемого ЛА отсутствует, т.е. накопленный научно-технологический задел в области авиастроения, заключенный в множестве Wr (1.3), обеспечивает создание перспективного ЛА, отвечающего набору многодисциплинарных требований. Облик перспективного ЛА в первом приближении описывает вектор ЛТХ r*. Дальнейшее уточнение полученного таким образом предварительного облика ЛА проводится на стадиях структурного, обликового и компоновочного проектирования.

Вариант 2: ${{\Phi }}\left( {r{\text{*}}} \right) \ne 0~$. Создание перспективного ЛА, удовлетворяющего комплексу заявленных требований, на текущем уровне научно-технологического развития авиастроения, невозможно, т.е. существует дефицит функциональных свойств проектируемого ЛА. Этот дефицит может быть устранен за счет расширения границ ${{r}_{{i~{\text{min}}}}},{{r}_{{i~{\text{max}}}}},~,i = \overline {1,n} $, исходного множества допустимых значений ЛТХ Wr в результате внедрения новых научно-технологических решений. Иными словами, с учетом возможности внедрения передовых технологий формируется новое множество ЛТХ $W_{r}^{*}$, в пределах которого обеспечивается одновременное выполнение условий (1.1).

Рассмотрим эту задачу более подробно. Ранее мы уже упоминали о том, что, задача устранения дефицита функциональных свойств перспективного ЛА может быть интерпретирована как обобщенная задача оптимизации [1], в которой множество допустимых значений ЛТХ Wr становится дополнительным параметром оптимизации:

(2.9)
${{\Phi }}(W_{r}^{*},r{\text{*}}) = \mathop {\min }\limits_{{{W}_{r}} \in {{\Xi }}} \mathop {\min }\limits_{r \in {{W}_{r}}} {{\Phi }}\left( r \right) = 0,$
где Ξ – семейство всех (любых) множеств допустимых значений ЛТХ, которые могут быть сформированы с учетом внедрения перспективных научно-технологических решений.

Как подчеркивалось выше, внедрение любого перспективного научно-технологического решения требует финансовых и временных затрат. Учитывая это, определим затраты на устранение дефицита функционального свойства, связанного с i-й ЛТХ, на основе выражения:

(2.10)
${{D}_{i}}\left( r \right) = \left\{ \begin{gathered} d_{i}^{1}\left( {{{r}_{i}} - {{r}_{{i{\text{min}}}}}} \right),\quad {{r}_{i}} < {{r}_{{i{\text{min}}}}}, \hfill \\ 0,\quad {{r}_{{i{\text{min}}}}} \leqslant {{r}_{i}} \leqslant {{r}_{{i{\text{max}}}}}, \hfill \\ d_{i}^{2}\left( {{{r}_{i}} - {{r}_{{i{\text{max}}}}}} \right),\quad {{r}_{i}} > {{r}_{{i{\text{max}}}}}. \hfill \\ \end{gathered} \right.$

Функция $d_{i}^{1}\left( {{{r}_{i}} - {{r}_{{i{\text{min}}}}}} \right)$ в (2.10) является строго положительной и выражает финансовые или временные затраты, связанные с разработкой и внедрением научно-технологических решений, обеспечивающих расширение диапазона возможных значений ЛТХ ri путем уменьшения ее минимального ${{r}_{{i{\text{min}}}}}$ значения.

Соответственно функция $d_{i}^{2}\left( {{{r}_{i}} - {{r}_{{i{\text{max}}}}}} \right)$ также строго положительная и используется для оценки финансовых или временных затрат, связанных с разработкой и внедрением научно-технологических решений, направленных на увеличение максимального ${{r}_{{i{\text{max}}}}}$ значения ЛТХ ${{r}_{i}}$.

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

С учетом (2.10) обобщенные затраты на устранение дефицита функциональных свойств проектируемого ЛА по всему набору ЛТХ можно определить следующим образом:

(2.11)
$D\left( r \right) = \mathop \sum \limits_{i = 1}^n {{D}_{i}}\left( r \right).$

Заметим, что влияние перспективных научно-технических решений проявляется не только изменением размеров исходного множества ЛТХ Wr, но и его формы, поскольку внедрение этих решений корректирует в том числе характер связей между компонентами вектора ЛТХ.

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

Учитывая это, в [1] предложен алгоритм приближенного построения множества $W_{r}^{*}$ в классе параллелепипедов за счет расширения границ ${{r}_{{i~{\text{min}}}}},{{r}_{{i~{\text{max}}}}},~i = \overline {1,n} $, исходного параллелепипеда Wr. При этом, к сожалению, предложенный алгоритм не учитывает возможные затраты на реализацию перспективных научно-технических направлений, что может привести к заведомо нереализуемым проектным решениям. В настоящей статье описывается развитие алгоритма, предложенного в [1], позволяющее определить приоритетные направления научно-технологического развития с учетом возможных затрат на их развитие и внедрение.

Алгоритм объединяет ряд последовательных шагов, содержание которых приводится ниже.

Шаг 1. Формируется начальное приближение $W_{r}^{0}$ для множества допустимых значений ЛТХ, в качестве которого принимается исходный параллелепипед Wr:

(2.12)
$W_{r}^{0} = \left\{ {r:r_{{i~{\text{min}}}}^{0} \leqslant {{r}_{i}} \leqslant r_{{i~{\text{max}}}}^{0},i = \overline {1,n} } \right\},$
где $r_{{i~{\text{min}}}}^{0} = {{r}_{{i{\text{min}}}}}$, $r_{{i~{\text{max}}}}^{0} = {{r}_{{i{\text{max}}}}},i = \overline {1,n} $.

Дальнейшая процедура, обеспечивающая получение множества $W_{r}^{*}$, направленного на устранение дефицита функциональных свойств проектируемого ЛА, базируется на методе оптимизации нулевого порядка – методе покоординатного спуска [6]. Заметим, что описываемая ниже процедура во многом повторяет алгоритм параметрической оптимизации доверительных множеств в классе кубов, изложенный в [5].

Шаг 2. Последовательно варьируется значение ${{r}_{{1~{\text{min}}}}}$ в сторону уменьшения:

(2.13)
$r_{{1~{\text{min}}}}^{k} = r_{{1~{\text{min}}}}^{0} - k\delta {{r}_{1}},\quad k = 1,2,3 \ldots ,$
где ${{\delta }}{{r}_{1}}~$ – заданный шаг поиска по компоненте ${{r}_{1}}$. При этом предельные значения всех других компонент вектора r принимаются равными их исходным значениям: ${{r}_{{1~{\text{max}}}}} = r_{{1~{\text{max}}}}^{0},{{r}_{{i~{\text{min}}}}} = r_{{i~{\text{min}}}}^{0}$, ${{r}_{{i~{\text{max}}}}} = r_{{j~{\text{max}}}}^{0},i = \overline {2,n} $. Таким образом, для каждого нового значения $r_{{1~{\text{min}}}}^{k}$ имеем приближение $W_{r}^{k} \in {{\Xi }}$ для множества $W_{r}^{*}$, задаваемое совокупностью условий:

(2.14)
$W_{r}^{k} = \left\{ {\begin{array}{*{20}{c}} {r:r_{{1~{\text{min}}}}^{k} \leqslant {{r}_{1}} \leqslant r_{{1~{\text{max}}}}^{0},~r_{{j~{\text{min}}}}^{0} \leqslant {{r}_{j}} \leqslant r_{{j~{\text{max}}}}^{0},~j = \overline {2,n} ~} \end{array}~} \right\}$.

Шаг 3. Для множества $W_{r}^{k}$ решается задача оптимизации

(2.15)
$D\left( r \right) \to \mathop {{\text{min}}}\limits_{r \in W_{r}^{k}} ,$
где D(r) определяется в виде (2.11) и учитывается ограничение

(2.16)
${{\Phi }}\left( r \right) = 0.$

В результате находится вектор ${{r}^{k}}$, обеспечивающий выполнение многодисциплинарных требований, предъявляемых к перспективному ЛА с учетом внедрения перспективных научно-технических решений, и оценка затрат на их реализацию $D({{r}^{k}})$.

Шаг 4. Выполняется перебор значений $r_{{1~{\text{min}}}}^{k},k = 1,2,3 \ldots $, в результате чего находится, во-первых, граница $r_{{1~{\text{min}}}}^{*}$, во-вторых, соответствующее этой границе множество $W_{r}^{*} \in {{\Xi }}$, определяемое условиями (2.14), а также вектор ЛТХ $r{\text{*}}(r_{{1~{\text{min}}}}^{*})$, при которых достигается минимум затрат $D(r{\text{*}}(r_{{1~{\text{min}}}}^{*})) = d_{1}^{1}(r_{{1~{\text{min}}}}^{0} - r_{{1~{\text{min}}}}^{*})$ на разработку и внедрение новых научных и технологических решений, использование которых обеспечивает выполнение заявленных требований ${{\Phi }}(r{\text{*}}(r_{{1~{\text{min}}}}^{*})) = 0$.

Шаг 5. Оценивается эффективность использования перспективных технологий, направленных на увеличение значения $r_{{1~{\text{max}}}}^{0}$. Для этого последовательно варьируется значение ${{r}_{{1~{\text{m}}ax}}}$ в сторону увеличения:

(2.17)
$r_{{1~{\text{max}}}}^{k} = r_{{1~{\text{max}}}}^{0} + k\delta {{r}_{1}},\quad k = 1,2,3 \ldots $

При этом предельные значения всех других компонент вектора r принимаются равными их исходным значениям:

${{r}_{{1~{\text{m}}in}}} = r_{{1~{\text{min}}}}^{0},\quad {{r}_{{ij{\text{min}}}}} = r_{{j{\text{min}}}}^{0},\quad {{r}_{{j~{\text{max}}}}} = r_{{j{\text{max}}}}^{0},\quad j = \overline {2,n.} $

Таким образом, для каждого нового значения $r_{{1~{\text{max}}}}^{k}$ имеем приближение $W_{r}^{k} \in {{\Xi }}$ для множества $W_{r}^{*}$, задаваемое совокупностью условий:

(2.18)
$W_{r}^{k} = \left\{ {\begin{array}{*{20}{c}} {r:r_{{1~{\text{min}}}}^{0} \leqslant {{r}_{1}} \leqslant r_{{1~{\text{max}}}}^{k},~r_{{j~{\text{min}}}}^{0} \leqslant {{r}_{j}} \leqslant r_{{j~{\text{max}}}}^{0},~j = \overline {2,n} ~} \end{array}~} \right\}.$

Для каждого приближения $W_{r}^{k}$ решается задача условной оптимизации (2.15), (2.16). В итоге определяется вектор ЛТХ $r{\text{*}}(r_{{1~{\text{max}}}}^{*})$, при которых достигается минимум затрат $D(r{\text{*}}(r_{{1~{\text{max}}}}^{*}))$ = = $d_{1}^{2}(r_{{1~{\text{max}}}}^{*} - r_{{1~{\text{max}}}}^{0})$, связанных с внедрением новых научных и технологических решений, использование которых обеспечивает выполнение заявленных требований ${{\Phi }}(r{\text{*}}(r_{{1~{\text{max}}}}^{*})) = 0$.

Аналогичным образом оценивается возможность выполнения заданных требований, предъявляемых к перспективному ЛА, за счет изолированной коррекции предельных значений $r_{{j~{\text{min}}}}^{0},r_{{j~{\text{max}}}}^{0},~j = \overline {2,n,} $ всех других компонент вектора ЛТХ.

В результате работы описанного алгоритма определяется набор векторов $r{\text{*}}(r_{{i~{\text{min}}}}^{*})$, $~r{\text{*}}(r_{{i~{\text{max}}}}^{*})$, $~i = \overline {1,n,} $ описывающих в первом приближении варианты обликов перспективного ЛА, которые могут быть получены за счет внедрения новых научно-технических решений и удовлетворяют набору заданных требований:

${{\Phi }}(r{\text{*}}(r_{{~i{\text{min}}}}^{*})) = 0,{{\Phi }}(r{\text{*}}(r_{{~i{\text{max}}}}^{*})) = 0$
с оценками затрат $D(r{\text{*}}(r_{{i~{\text{min}}}}^{*})),D(r{\text{*}}(r_{{i~{\text{max}}}}^{*}))$ на их реализацию.

Решение принимается в пользу того варианта $\hat {r}$, реализация которого сопряжена с наименьшими затратами:

(2.19)
$D\left( {\hat {r}} \right) = \mathop {\min }\limits_{i = \overline {1,n} } \{ D(r{\text{*}}(r_{{i~{\text{min}}}}^{*})),D(r{\text{*}}(r_{{i~{\text{max}}}}^{*}))\} .$

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

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

Практическая реализация приведенного алгоритма требует многократного решения задач условной оптимизации, подобных (2.15), (2.16). Для этого может быть использован алгоритм на основе метода направленного случайного поиска [7]. Существенным преимуществом такого алгоритма является возможность распараллеливания вычислений. Работа алгоритма требует генерирования реализаций непрерывной равномерно-распределенной на интервале (0.1) случайной величины, для чего применен мультипликативный датчик равномерно-распределенных псевдослучайных чисел [8]:

(2.20)
${{{{\alpha }}}_{i}} = \left\{ {M{{\alpha }_{{i - 1}}}} \right\},$
где $i = 1,2, \ldots ,M = 1220703125,{{\alpha }_{0}} = {{2}^{{ - 32}}}$, фигурные скобки в (2.20) означают дробную часть числа.

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

$P\left( {X(\omega } \right) = {{x}_{k}}) = \frac{1}{K},\quad k = \overline {1,K} ,$
использован алгоритм:
(2.21)
${{\gamma }_{i}} = \mathop {\min }\limits_k \left( {\left( {{{\alpha }_{i}} - \frac{k}{K}} \right) < 0} \right),$
где ${{\alpha }_{i}}$ находится в соответствие с (2.20).

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

3. Пример оценки дефицита функциональных свойств перспективного самолета транспортной категории. Рассмотрим процесс формирования облика перспективного ЛА на примере классического транспортного рампового самолета, который должен обеспечивать беспосадочный перелет на дальность свыше 5000 км, перевозя при этом полезную нагрузку не менее 60 т.

Следуя [1], при формировании облика проектируемого самолета ограничимся только требованиями весовой и динамической эффективности. Это допущение никак не повлияет на реализацию описанного в разд. 2 алгоритма, но обеспечит более наглядное представление результатов.

Исходя из заданных требований к перспективному ЛА, определены предельные значения показателей ${{I}_{{j{\text{min}}}}},{{I}_{{j{\text{max}}}}},j = \overline {1,N} $ (табл. 1).

Таблица 1.

Состав показателей, определяющих требования к перспективному самолету

Наименование показателя, Ij Минимальное значение, ${{I}_{{j{\text{min}}}}}$ Максимальное значение, ${{I}_{{j{\text{max}}}}}$
Весовая отдача 0.28 0.8
Относительная масса конструкции 0.1 0.3
Крейсерская скорость, км/ч 300 900
Скорость отрыва, км/ч 0 100
Длина пробега, м 400 1200
Скорость посадочная, км/ч 0 100
Длина разбега, м 300 1200
Время разбега, с 0 50

Таким образом, приемлемым полагается облик проектируемого ЛА, для которого выполняется следующая совокупность условий:

(3.1)
$\left\{ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {0.28 \leqslant {{I}_{1}}\left( r \right) \leqslant 0.8,} \\ {0.1 \leqslant {{I}_{2}}\left( r \right) \leqslant 0.3,} \\ {300 \leqslant {{I}_{3}}\left( r \right) \leqslant 900,~} \end{array}} \\ {0 \leqslant {{I}_{4}}\left( r \right) \leqslant 100,} \end{array}} \\ {400 \leqslant {{I}_{5}}\left( r \right) \leqslant 1200,} \\ {0 \leqslant {{I}_{6}}\left( r \right) \leqslant 100,} \\ {300 \leqslant {{I}_{7}}\left( r \right) \leqslant 1200,} \\ {0 \leqslant {{I}_{8}}\left( r \right) \leqslant 50.} \end{array}} \right.$

В табл. 2 приведен перечень ЛТХ (компонентов вектора r), используемых для расчета интегральных показателей ${{I}_{j}}\left( r \right),j = \overline {1,8} $, представленных в табл. 1.

Таблица 2.

Вектор ЛТХ r, используемый для расчета интегральных показателей, отражающих соответствие проектируемого самолета предъявленным требованиям

Наименование ЛТХ, ri Минимальное значение, ${{r}_{{i{\text{min}}}}}$ Максимальное значение, ${{r}_{{i{\text{max}}}}}$
Коэффициент лобового сопротивления на посадке 0.33 0.34
Количество двигателей, шт. 4 4
Тяга двигателя максимальная, кгс 8290 8300
Максимальное значение коэффициента подъемной силы на взлете 1.9 1.95
Максимальное значение коэффициента подъемной силы на посадке 2.4 2.45
Вес перевозимой целевой нагрузки, кг 60 000 60 000
Дальность полета, км 5000 5000
Степень двухконтурности двигателя 4.4 4.5
Стреловидность крыла, град. 25 27
Площадь крыла, м2 200 210
Размах крыла, м 50 60
Коэффициент лобового сопротивления на крейсерском режиме 0.01 0.02
Время разбега, с 0 50

Конкретные значения этих ЛТХ получены на основе обобщения накопленного общемирового опыта проектирования транспортных рамповых самолетов и отражают текущий научно-технический уровень развития авиастроения. Значения, присутствующие в табл. 2, задают в 13-мерном пространстве исходное множество ЛТХ ${{W}_{r}}$.

Зависимости ${{I}_{1}}\left( r \right) - {{I}_{8}}\left( r \right)$, связывающие показатели, описанные в табл. 1, с компонентами вектора ЛТХ r, определены на основе фундаментальных положений теории полета [9, 10] с использованием соотношений, приведенных в табл. 3.

Таблица 3.

Соотношения, участвующие в расчете интегральных показателей ${{I}_{1}}\left( r \right) - {{I}_{8}}\left( r \right)$

Характеристики, участвующие в расчете интегральных показателей, Ij Выражение для расчета
Удлинение крыла $\lambda = \frac{{{{l}^{2}}}}{S}$
Эффективное удлинение крыла, ${{\lambda }_{{{\text{эф}}}}}$ ${{\lambda }_{{{\text{эф}}}}} = \frac{\lambda }{{1 + \pi ~\frac{\lambda }{{100}}~{\text{co}}{{{\text{s}}}^{2}}}}$
Коэффициент отвала поляры $A = \frac{1}{{\pi \lambda {{~}_{{{\text{эф}}}}}}}$
Коэффициент подъемной силы на крейсерском режиме, ${{C}_{{y~{\text{кр\;}}}}}$ ${{C}_{{y~{\text{кр\;}}}}} = \sqrt {\frac{{{{C}_{{x0}}}}}{A}} $
Коэффициент лобового сопротивления на крейсерском режиме, ${{C}_{{x~{\text{кр}}}}}$ ${{C}_{{x~\,{\text{кр}}}}} = {{с}_{{x0}}} + A~C_{{y~\,{\text{кр}}}}^{2}$
Статический удельный расход топлива номинальный ${{с}_{{р0}}} = \frac{{0.8}}{{1 + 0.5{\text{ }}\sqrt m }}$
Скорость горизонтального полета (крейсерская скорость), ${{V}_{{{\text{г}}{\text{.п}}}}}$ ${{V}_{{{\text{г}}{\text{.п}}}}} = \sqrt {\frac{{2~{{G}_{{{\text{взл}}}}}}}{{\rho ~S~{{C}_{y}}}}} $
Коэффициент подъемной силы в момент отрыва ${{C}_{{y{\text{отр}}}}} = 0.7~{{C}_{{y{\text{maxвзл}}}}}$
Коэффициент подъемной силы на посадке ${{C}_{{y{\text{пос}}}}} = 0.7~{{C}_{{y{\text{maxпос}}}}}$
Аэродинамическое качество на посадке $K = \frac{{{{C}_{{y{\text{пос}}}}}}}{{{{C}_{{x{\text{пос}}}}}}}$
Аэродинамическое качество на крейсерском режиме $K = \frac{{{{C}_{{y~{\text{кр}}}}}}}{{{{C}_{{x{\text{\;кр}}}}}}}$
Потребная тяга, ${{Р}_{{\text{п}}}}$, кгс ${{Р}_{п}} = \frac{{{{G}_{{{\text{взл}}}}}}}{K}$
Скорость отрыва, м/с ${{V}_{{{\text{отр}}}}} = \sqrt {\frac{{2~{{G}_{{{\text{взл}}}}}}}{{\rho ~S~{{C}_{{y{\text{отр}}}}}}}} $
Скорость посадочная, м/с ${{V}_{{{\text{пос}}}}} = \sqrt {\frac{{2~{{G}_{{{\text{пос}}}}}}}{{\rho ~S~{{C}_{{y{\text{пос}}}}}}}} $
Длина пробега, ${{L}_{{{\text{пр}}}}}$, м ${{L}_{{{\text{пр}}}}} = \frac{{{{{\left( {{{V}_{{{\text{пос}}}}} \mp W} \right)}}^{2}}}}{{g~\left( {\frac{1}{{{{K}_{{{\text{пос}}}}}}} + {{f}_{{{\text{тр}}}}}} \right)}}$
Длина разбега, ${{L}_{р}}$, м ${{L}_{р}} = \frac{{{{{\left( {{{V}_{{{\text{отр}}}}} \mp W} \right)}}^{2}}}}{{2~g~\left( {\frac{{{{Р}_{{\text{п}}}}}}{{{{G}_{{{\text{взл}}}}}}} - {{f}_{{{\text{тр}}}}}} \right)}}$
Удельный расход топлива на крейсерском режиме, ${{C}_{{{\text{уд}}}}}$, кг/кгсч ${{C}_{{{\text{уд}}}}} = {{C}_{{p0}}} + \frac{{0.4~M}}{{1 + 0.027~H}}$
Масса топлива, $~{{m}_{{\text{т}}}}$, кг ${{m}_{{\text{т}}}} = {{q}_{{{\text{км}}}}}l$
${{q}_{{{\text{км}}}}} = \frac{{{{q}_{{\text{ч}}}}}}{{{{V}_{{{\text{г}}{\text{.п}}}}}}}$
${{q}_{{\text{ч}}}} = {{C}_{{{\text{уд}}}}}~{{Р}_{{\text{п}}}}$

Обозначения, принятые в выражениях из табл. 3, имеют следующий смысл: $l,S,~{{\chi }}$ – соответственно размах, площадь и стреловидность крыла, ${{C}_{{x0}}}$ – коэффициент лобового сопротивления при нулевом угле атаки, $m$ – часовой расход топлива, ${{G}_{{{\text{взл}}}}}$ – взлетный вес, S – площадь крыльев, $~{{C}_{{y{\text{maxвзл}}}}}$ – максимальное значение коэффициента подъемной силы на взлете, ${{C}_{{y{\text{maxпос}}}}}$ – максимальное значение коэффициента подъемной силы на посадке, ${{{{\rho }}}_{0}}$ – плотность атмосферы на поверхности земли, ${{G}_{{{\text{пос}}}}}$ – посадочный вес, $~{{V}_{{{\text{пос}}}}}$ – посадочная скорость, W – скорость ветра, ${{K}_{{{\text{пос}}}}}$ – аэродинамическое качество на посадке, ${{f}_{{{\text{тр}}}}}$ – коэффициент трения при движении по взлетно-посадочной полосе (ВПП), ${{V}_{{{\text{отр}}}}}$ – скорость отрыва, ${{P}_{{{\text{ср}}}}}$ – среднее значение тяги двигателя в момент разбега самолета по ВПП, g – ускорение силы тяжести, H – крейсерская высота, M – число Маха на крейсерском режиме, ${{q}_{{{\text{км}}}}}$ – километровый расход топлива, ${{\Delta }}l$ – длина участка горизонтального полета, ${{q}_{{\text{ч}}}}$ – среднечасовой расход топлива на крейсерском режиме.

В результате решения задачи оптимизации (2.8) выяснилось, что в границах исходного множества ЛТХ ${{W}_{r}}$ невозможно определить облик самолета, удовлетворяющий совокупности условий (3.1), т.е. имеет место дефицит функциональных свойств. Следовательно, необходимо найти минимально возможное расширение границ множества ${{W}_{r}}$ допустимых значений ЛТХ и построить новое множество $W_{r}^{{\text{*}}}$ с использованием алгоритма, описанного в разд. 2.

Как указывалось выше, получение множества $W_{r}^{{\text{*}}}$ должно проводиться с учетом затрат на устранение дефицита функционального свойства, связанных с внедрением перспективных научных и технологических решений, которые определяются на основе (2.10), (2.11). В работе использованы линейные зависимости $d_{i}^{1}\left( {{{r}_{i}} - {{r}_{{i{\text{min}}}}}} \right),\,\,d_{i}^{2}\left( {{{r}_{i}} - {{r}_{{i{\text{max}}}}}} \right),\,\,i = \overline {1,13} $, выражающие время, которое, по мнению экспертов, необходимо для практического внедрения тех или иных технологий, позволяющих расширить диапазоны ${{r}_{{i~{\text{min}}}}},\,\,{{r}_{{i~{\text{max}}}}},\,\,i = \overline {1,13} ,$ соответствующих компонентов вектора ЛТХ r.

В результате реализации алгоритма (разд. 2) было определено новое множество допустимых значений ЛТХ $W_{r}^{{\text{*}}}$ и выделены возможные направления расширения границ исходного множества Wr  допустимых значений ЛТХ, представленные в табл. 4 темным цветом.

Таблица 4.

Множество допустимых значений ЛТХ $W_{r}^{*}$, построенное с учетом внедрения перспективных научно-технологических решений

Наименование ЛТХ, ri Минимальное значение, ${{r}_{{i{\text{min}}}}}$ Максимальное значение, ${{r}_{{i{\text{max}}}}}$
Коэффициент лобового сопротивления на посадке 0.33 0.34
Количество двигателей, шт. 4 5
Тяга двигателя максимальная, кгс 8290 8551
Максимальное значение коэффициента подъемной силы на взлете 1.9 1.98
Максимальное значение коэффициента подъемной силы на посадке 2.4 2.45
Вес перевозимой целевой нагрузки, кг 60 000 60 000
Дальность полета, км 5000 5000
Степень двухконтурности двигателя 4.4 4.7
Стреловидность крыла, град. 25 27
Площадь крыла, м2 200 210
Размах крыла, м 50 60
Коэффициент лобового сопротивления на крейсерском режиме 0.01 0.02
Масса взлетная, кг 80 000 82 113

Как следует из табл. 4, приоритетное значение для создания транспортного рампового самолета, отвечающего всему набору требований, сформулированных в табл. 1, имеют технологии, связанные с совершенствованием силовой установки. Удовлетворить заявленные требования удается за счет использования пяти двигателей, способных обеспечить максимальную суммарную взлетную тягу на уровне 8551 кгс и обладающих степенью двухконтурности, равной 4.7.

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

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

Заключение. Предложена формализация задачи оценки дефицита функциональных свойств проектируемого ЛА как обобщенной задачи оптимизации, в которой множество допустимых значений ЛТХ выступает в качестве дополнительного параметра оптимизации. Разработан алгоритм решения этой задачи с учетом накопленного научно-технического задела в области авиастроения, а также финансовых и временных ограничений.

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

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

Список литературы

  1. Анисимов К.С., Евдокименков В.Н., Красильщиков М.Н., Сыпало К.И., Топоров Н.Б. Оптимизация процесса совершенствования авиационных комплексов на базе концепции функционального проектирования // Изв. РАН. ТиСУ. 2022. № 1. С. 105–123.

  2. Evdokimenkov V.N., Toporov N.B., Vavilov D.S. Optimizing the Preliminary Design Profile of an Aircraft // Russian Engineering Research. 2022. V. 42. № 6. P. 599–602

  3. Демьянов В.Ф., Малоземов В.Н. Введение в минимакс. М.: Наука, Физматлит, 1976.

  4. Кибзун А.И., Лебедев А.А., Малышев В.В. О сведении задачи с вероятностными ограничениями к эквивалентной минимаксной // Изв. АН СССР. Техн. кибернетика. 1984. № 4. С. 73–80.

  5. Евдокименков В.Н., Динеев В.Г., Карп К.А. Инженерные методы вероятностного анализа авиационных и космических систем. М.: Физматлит, 2010.

  6. Васильев Ф.П. Лекции по методам решения экстремальных задач. М.: Изд-во МГУ, 1974.

  7. Аттеков А.В., Галкин С.В., Зарубин В.С. Методы оптимизации. М.: Изд-во МГТУ им. Н.Э. Баумана, 2003. 440 с.

  8. Вероятностные методы в вычислительной технике: Учеб. пособие / Под ред. А.Н. Лебедева и Е.А. Чернявского. М.: Высш. шк., 1986.

  9. Балакин В.Л., Лазарев Ю.Л. Динамика полета самолета. Расчет траекторий и летных характеристик: Конспект лекций. Самара: Самарск. гос. аэрокосм. ун-т, 2002. 56 с.

  10. Аэромеханика самолета / Под ред. А.Ф Бочкарева и В.В. Андреевского. М.: Машиностроение, 1985.

Дополнительные материалы отсутствуют.