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

УДК 658.512:621.452

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

© 2011 Ф. В. Гречников, С. Ф. Тлустенко

Самарский государственный аэрокосмический университет имени академика С.П. Королёва

(национальный исследовательский университет)

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

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

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

Модель системы ТП в виде графа позволяет ввести индексные функции показателей производственных переходов в виде связей по дугам графа £(ху), которые могут иметь достаточно сложный вид. Функции задаются аналитически или алгоритмически.

Для процессов преобразований сборочного пространства при выполнении ТП операции представляются как п-шаговый процесс реализации решений с операторами

переходов А = / (я_1, ) = хг-1 + = хг. Ре-

зультат перехода определяется функцией

/ (^ и] ) = / {^ X )} .

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

стояния, к которому придет система в конце данного 1-го шага. При этом гарантируется, что управление, выбранное на любом шаге, не локально лучшее, а лучшее с точки зрения процесса в целом. Функция дискретного перехода соответствует содержанию задачи и адекватному переходу в системе автоматизированного проектирования технологических процессов (САПР ТП) сборки.

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

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

Рг-1 ( Х_1 ) =

,г = 1,2....,п, (1)

£ (х_1, х (и ))-(х (и}))

где Б - функция качества перехода и

р, (Хп ) = 0.

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

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

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

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

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

Рассмотрим параметрические формы

а1(А1,А2,t)и а2(А3,А4,t), описывающие состояния ТС:

) = А1 + А2 w;2(t);

w'(t) = А3 _ А4^),

(2)

где w

b

динамика выполнения операции;

ч - динамика подготовки производства для

выполнения операции; АI - параметр ресурсов; А2 - параметр эффективности; А3 - параметр, характеризующий индекс качества операции; А4 - параметр, характеризующий неустойчивость по некоторому критерию операции (фактор рисков). Уравнение, описывающее состояние ТС, представляет частный случай уравнения Риккати, которое в первом приближении решается в квадратурах, причем считаем, что ч^) и чф) - наблюдаемые фазовые координаты соответствующих векторов состояний в начале или в конце фиксированного промежутка времени

наблюдений за группой а1( А15 А2, t),

а2( А3, А4, t).

В качестве достаточных приближений значений для ч^) и чЬ^) в детерминированных условиях процесса выбираем:

w,

(t) = tg (yl a A2t.)

w* (t) = th(tJ A3 A41)

(3)

(4)

VA7

Решения уравнений (l)-(2), полученные для wb(t) и w(t), удовлетворяют аналогичным уравнениям в определённой метрике t.

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

При более точных условиях измерения фазовых координат wb(t) и w(t) в этом случае исходные уравнения имеют вид

wb (t) = A1 + A2 wb (t) + A5 sintwb (t)t +

+fb0 ]« *(t) - A5 C0SK (t)t + fb0] £*(t);

ч (t) = Л + А4 ) _ w;2(t) А6 +

+А7 8Ш[ ^ ^ + фЬ0а ) _ (6)

_А?С08[w*(t)t + фЫ)]Р ),

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

вающее характер колебаний по а ; р ^) -функция, характеризующая поставки для последующих процессов сборки при установленной их ритмичности, получаемая после подстановки ) и ч (t) в уравнение

динамики колебаний по Ь ; А5 - параметр, характеризующий устойчивость операции при межцеховых колебаниях объёмов и номенклатуры поставок материалов и комплектующих; А6 - параметр, характеризующий степень освоения в цехе указанных сборочных процессов; А7 - параметр, учитывающий возможные переналадки и переустройства рабочих мест исполнителей; фЬ0 -начальная фаза сезонных изменений.

Сгруппировав слагаемые ряда в уравнениях (5)-(6) по степеням t: t, t ,..., получаем несколько уровней классификации состояния системы ТП. Получение решения в виде ряда по степеням t определяет, какие параметры влияют на форму кривых

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

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

*

функции состояния ТС разложением в ряд базовой целевой функции по базисным векторам, причем коэффициенты разложения являются собственными значениями линейного оператора (Я1,1п), которым соответствуют вероятностные представления, рассчитываемые как плотности вероятностей такого состояния. Изображение оригинала заданной вероятности в комплексной области рассчитывается как

--------Г, (7)

(р +1)п+Г

где ^ - количество факторов и их выборочных значений в распределениях-композициях; р - комплексная величина; 1

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

Ортогональные тензоры Q: Е —® Е осуществляют повороты и зеркальные отображения (рис.1), при этом они сохраняют длины векторов, углы между ними и переводят один ортонормированный базис в другой. Определим понятие инвариантности тензора А: Е ——-® Е относительно ортогонального тензора 0: Е —2—» Е .

Л

А(зс) _

Г

Рис. 1 Преобразования технологического пространства сборки

Определение 1. Тензор А называется инвариантным относительно тензора Q, если А удовлетворяет равенству

А = Q о А о Q-Г. (8)

Уточним геометрический смысл этого понятия. Пусть Q - ортогональный тензор, осуществляющий поворот любого вектора в пространстве Е (ё1шЕ=3) на заданный угол

вокруг оси, ортогональной плоскости чертежа (рис.1).

Воспользуемся для наглядности искусственным построением: несмотря на то что О действует из пространства Е в пространство Е, будем считать, что О действует из пространства Е в пространство Е и «забудем» на время, что Е= Е. Е -«повернутое» пространство, в которое оператор

О: Е —Е переводит пространство Е. Обратный оператор О-1 : Е ——® Е возвращает пространство Е в исходное. Пусть на пространстве Е задан тензор А: Е —2® Е. Найдем в «повернутом» пространстве Е такой тензор А : Е —^® Е, который оказывал бы на любой вектор х є Е такое же действие, которое оказывает тензор А на прообраз вектора х (т.е. на вектор у = О-1(х)) в пространстве Е. Искомый тензор А должен удовлетворять равенству

"х є Е А(х) = О(А(у)) = О(А(д-1(х))), которое эквивалентно тензорному равенству

А = о о а о о-1.

Отождествим векторные пространства Е и Е (они в действительности неразличимы), в результате (3') переходит в (3). Рис.1 демонстрирует геометрический смысл понятия инвариантности тензора относительно ортогонального преобразования. Ортогональный тензор О, относительно которого тензор А инвариантен, принято называть симметрией тензора А. Понятие инвариантности тензора А относительно ортогонального тензора О обобщается и на тензорные функции.

Определение 2. Отображение Б: £(Е)^ £(Е) (называемое также тензорной функцией) инвариантно относительно ортогонального тензора О є £(Е), если оно удовлетворяет равенству

"А є £(Е),

О о Б(А) о О-1 = Б(О о А о О-1).

Определение 3. Отображение Б: £(Е)^ £(Е) инвариантно, если оно инвариантно относительно любого ортогонального тензора О є £(Е), т.е. для любого ортогонального

О є £(Е), для "А є £(Е)

О о Б(А) о О-1 = Б(О о А о О-1). (9)

Требование инвариантности отображения является настолько жестким, что в ряде случаев оно позволяет сразу выписать явный вид инвариантного отображения. Далее показано, что любое инвариантное отображение Б симметричных тензоров 8 е £(Е), &ше=3 может быть представлено в виде

Б(Б) = а1 + РБ + уБ2

5

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

Лемма 1. Если пространство Е раскладывается в прямую сумму

Е = ЕГ ФЕ2 Ф...ФЕк попарно ортого-

нальных подпространств, каждое из которых инвариантно относительно тензора А, и если отображение Б: £(Е)^ £(Е) инвариантно, то подпространства ЕГ....Ек инвариантны и относительно тензора Б(А).

Доказательство. Нетрудно видеть, что лемму достаточно доказать для двух взаимно ортогональных инвариантных относительно А подпространств Е1, Е2 таких, что Е = ЕГ Ф Е2.

Зададим два невырожденных отображения Q1:E ®Е и Q2:E ® Е:

1 на ЕГ

I I на Е 2

(12)

О1 Ч/"” О-1 = О1 на Е, (10)

| -1 на Е, | ,

022 Ч т Е 1 Н 021 = 02 на Е. (11) [ Т на Е 2 ]

Легко проверить, что они являются ортогональными. (Например, достаточно убедиться в выполнении необходимого и достаточного условия ортогональности 0'0*=Т). Из ортогональности отображения следует его линейность, таким образом, отображения 0і, 02 - ортогональные тензоры. Воспользуемся инвариантностью Б относительно любого ортогонального тензора, в частности, относительно ортогонального тензора Оі=10, и докажем инвариантность подпространства Е1 относительно тензора Б(А).

Инвариантность подпространства Е1 относительно тензора А и равенство Е = Ех Ф Е2 позволяют записать

(14)

"хєЦ Ах) єЦ Б(А)(х)=у+г

УєE, гєр.

Инвариантность отображения Б относительно О1 приводит к равенству

"А 01 о Б(А) о 01-1 = Б(01 о А о 01-1). (13)

Найдем значения левой и правой частей (13) для произвольного вектора

х є Ех: "х є Е1 :

01 о Б(А) о 01-1(х) = 01 о Б(А)(-Х1)=

= -01(у + г) = у - г.

Чтобы найти значение правой части

(13) от х є Ех, заметим, что тензоры АІр и

Е1

(0 о А о 01 1)р , индуцированные на инва-р1

риантном относительно тензора А и относительно тензора (01 о А о 01-1) подпространстве Ех, это позволяет записать "х є Ех

Б(01 о А о 01-1)(х) = Б(А)(х) = у + г. (15)

Из равенств (13)-(15) следует соотношение у - г = у + г (у є Ех, г є Е2), ко-

торое возможно только при г = 0, что равносильно утверждению ", х є Ех

Б(А)(х) = у є Е1, доказывающему инвариантность подпространства Е1 относительно тензора Б(А). Аналогично с помощью ортогонального тензора 02 (11) доказывается инвариантность подпространства Е2 относительно тензора Б(А).

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

тензора формулируется следующим образом.

Утверждение 1 (о симметричных тензорах). Тензор 8 : Е —^® Е симметричен тогда и только тогда, когда существует ор-тонормированный базис пространства Е, состоящий из собственных векторов Б.

Лемма 1 и это утверждение позволяют для любого симметричного тензора Б достаточно просто доказать симметричность тензора Б(Б) и его соосность с тензором Б для любого инвариантного отображения Б.

Определение 4. Симметричный тензор А є £(Е) соосен симметричному тензору В є £(Е), если любой собственный вектор

тензора В является и собственным вектором тензора А.

Замечание. Отношение соосности не является симметричным, т.е. из того, что тензор А соосен тензору В, еще не следует, что тензор В соосен тензору А. Например, единичный тензор в соответствии с опр.4 соосен любому симметричному тензору, обратное, очевидно, неверно.

Утверждение 2 (о соосности тензора Б(Б) тензору Б). Для любого симметричного тензора Б и для любого инвариантного отображения Б: £(Е)^ £(Е) тензор Б(Б) симметричен и соосен с тензором Б.

Доказательство. Ограничимся случаем трехмерного пространства (&шЕ=3). Пусть Б - симметричный тензор. По утв.1 существует базис (п1, п2, п3) пространства Е, состоящий из единичных собственных векторов тензора Б. Подпространства Е1, Е2, Е3 , натянутые на векторы (п1, п2, п3) соответственно, как известно, инвариантны относительно Б и, кроме того, попарно ортогональны и одномерны, это позволяет представить Е в виде их прямой суммы:

Е = Е1 Ф Е2 Ф Е3.

Из леммы 1 следует, что подпространства Е1, Е 2, Е 3 для инвариантного отображения Б также инвариантны относительно тензора Б(Б), а это означает, что векторы (п1, п2, п3) - собственные векторы тензора Б(Б), действительно

пк є Ек ^ Б(Б)(пк) = у є Ек ^ (16)

^ Б(Б)(пк) = РА. 1 >

Таким образом, для тензора Б(Б) найден ор-тонормированный базис (п1, п2, п3) пространства Е, состоящий из собственных векторов Б(Б), следовательно (утв.1), тензор Б(Б) симметричен. Докажем соосность симметричного тензора Б(Б) симметричному тензору Б. Для этого в соответствии с определением 3 достаточно доказать, что любой собственный вектор и тензора Б будет и собственным вектором тензора Б(Б). Пусть Ь - линейная оболочка вектора и, т.е. одномерное подпространство, натянутое на вектор и. Обозначим Ь* - ортогональное дополнение к Ь. Легко показать, что из того, что Б - симметричный тензор, следует, что подпространство Ь* инвариантно относи-

тельно Б, очевидно). Лемма 1 позволяет утверждать, что подпространства Ь и Ь* инвариантны также и относительно тензора Б(Б). Поскольку подпространство Ь по построению одномерно, то (см.(16))

и є Ь ® Б^Хи) = Ь и, т.е. и - собственный вектор и тензора Б(Б), что доказывает соосность тензора Б(Б) тензору Б. Доказанного утв.2 практически достаточно, чтобы выписать явное выражение тензора Б(Б) через тензор Б при сделанных предположениях. Особенно просто это в случае линейного отображения Б.

Разложение произвольного тензора А є £(Е) по диадному базису (§10 п.2) имеет вид

А = Ок ®■„ . (17)

к т

акт - коэффициенты матрицы тензора А в произвольном ортонормированном базисе (і1, і2, і3) пространства Е. Для симметричного тензора Б можно выбрать ортонормиро-ванный базис, состоящий из собственных векторов Б, в котором только три коэффициента в разложении (17) ненулевые. Действительно, для симметричного тензора Б существует ортонормированный базис (п1, п2, п3) из собственных векторов Б, матрица тензора Б в этом базисе диагональна, что приводит к следующему представлению этого тензора:

Б = X 1п1 ® п1 + X2п2 ® п2 + X3п3 ® п3 . (18)

Х1, X 2, X 3 - собственные числа Б, среди которых могут быть равные. Тензор Б(Б) со-осен с тензором Б согласно утверждению 2, следовательно, тензор Б(Б) имеет диагональную матрицу в базисе (п1, п2, п3) и для него возможно разложение по тем же трем диадам:

Б(Б) = СЦЇЇ1 ® п1 + а2п2 ® п2 + а3п3 ® п3 (19) В дальнейшем используется утверждение, представляющее и самостоятельный интерес.

Утверждение 3 (об ортогональном

преобразовании диады).

"а, Ь є Е, (іітЕ = 3, " ортогонального тензора 0 є £(Е)

О о а ® Ь о О-1 = 0(а) ® 0(Ь). (20)

Доказательство. "х е Е

^ о а ® Ь о Q"1 )(х) = ^ о а 0 Ь)^'1 (Х))=

= Q(a(b • Q-1(X))).

Учтем, что число, равное скалярному произведению Ь • Q-1(Х), можно представить в виде

Ь • Q-1 (Х) = Ь • Q*(Х) = Q(b) • Х.

Так как тензор Q ортогональный, сле-

1

довательно, Q- = Q . Используя это представление, получим равенство "Х е Е ^ о а 0 Ь о Q-1 )(Х) = рСа(р(Ь) • Х)) =

= ^(а) 0 Q(b))(Х), доказывающее утверждение 3.

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

Р(х, у) ^ Р(х, у, б). (21)

В общем случае п-местный предикат Р заменяется на п+1 - местный, всегда оставляя последнее место для последующей переменной состояния б, ассоциированной с этим предикатом Р. Содержательно это за-

ключается в том, что в некотором абстрактном сборочно-монтажном пространстве переменная состояния может принимать значения б1, б2, Бз... 810... Эти состояния

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

Библиографический список

1. Коротнев, Г.И. Топологические и тензорные методы описания производства летательных аппаратов [Текст] / Г.И. Коротеев

- М.: Полет.- 2003. - №4. - С. 17-23.

2. Ульянов, М.В. Ресурсно-эффективные компъютерные алгоритмы. Разработка и анализ. [Текст] / М.В. Ульянов - М.: Физ-матлит, 2008. - 304с.

3. Лазарсон, Э.В. Теория и методы решения многовариантных неформализованных задач выбора [Текст]: монография. - Пермь: Изд-во Перм. гос. техн. ун-та, 2008, 270 с.

THEORY REPRESENTATIONS MATHEMATICAL MODELS AND SEQUENCES OF PERFORMANCE OF TECHNOLOGICAL OPERATORS OF ASSEMBLAGE OF PLANES

©2011 F. V. Grechnikov, S. F. Tlustenko

Samara State Aerospace University named after academician S.P. Korolyov (National Research University)

Mechanical communications between collected elements of a design of the flying machine are modelled according to the scheme of spatial interrelation of assembly units of a design and initial, boundary conditions. The description of character of movement of an element of a product at assemblage is represented at quantitative level of a trajectory matrixes of conditions agrees the general count of assemblage of a product. Trajectories of movings are described by methods of analytical geometry and matrixes of movings, and also the operators, describing movement as vector size.

The operator, transformation, logic level, a contour, structure, interaction, bases.

Информация об авторах

Гречников Федор Васильевич, член - корреспондент Российской академии наук, доктор технических наук, профессор, заведующий кафедрой обработки металлов давлением, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Тел: (846) 267-43-02. Область научных интересов: обработка металлов давлением.

Тлустенко Станислав Федотович, кандидат технических наук, доцент кафедры обработки металлов давлением, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Тел: (846) 26746-01, 8-927-731-58-12. E-mail: [email protected]. Область научных интересов: обработка металлов давлением, сборка авиационных конструкций.

Grechnikov Fedor Vasilievich, term correlation of Sciences, Doctor of Technical Sciences, Professor, Head of the Department of Metal Forming, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). Phone: (846) 267-43-02. Area of research: metal forming.

Tlustenko Stanislav Fedotovich, Candidate of Technical Science, Associate professor of metal forming, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). Phone: (846) 267-46-01, 8-927-731-58-12. E-mail: [email protected]. Area of research: metal forming, assembly of aircraft structures.