CC BY
49
DOI: 10.15593/2224-9982/2018.54.06 УДК 533.6.013.1
М.Ю. Егоров 1, Р.В. Мормуль2
2 Пермский национальный исследовательский политехнический университет, Пермь, Россия
1 ПАО НПО «Искра», Пермь, Россия
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЯЗКОУПРУГОГО ПОВЕДЕНИЯ ЗАРЯДА ТТ ПРИ АКУСТИЧЕСКОЙ НЕУСТОЙЧИВОСТИ ВНУТРИКАМЕРНОГО ПРОЦЕССА. ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ
Впервые прямым численным моделированием получен эффект «биения» (интермодуляции) первого рода (нелинейное наложение друг на друга колебаний газа с малой разностью частот) в камере сгорания ракетного двигателя твердого топлива (РДТТ) специального назначения.
Впервые определена оценка влияния акустической неустойчивости внутрикамерного процесса ракетного двигателя специального назначения с учетом пространственно-временной локализации динамической нагрузки на вибродинамическую прочность заряда ТТ при его температурно-временной вязкоупругой зависимости модуля упругости (ядра релаксации Максвелла) и реологии. Разработанный математический аппарат, созданный на его базе комплекс прикладных программ, разработанная методика расчета, проведенные методические исследования дают возможность существенно повысить надежность, улучшить энергомассовые, прочностные, эксплуатационные и другие характеристики РДТТ.
Ключевые слова: ракетный двигатель твердого топлива, акустическая неустойчивость, математическое моделирование, вычислительный эксперимент, вязкоупругое поведение, динамическая прочность, торцевое горение, камера сгорания, метод конечных элементов, определяющие соотношения.
M.Yu. Egorov1, R.V. Mormul2
1 Perm National Research Polytechnic University, Perm, Russian Federation 2 PJSC Research and Production Association "Iskra", Perm, Russian Federation
MATHEMATICAL MODELING OF THE VISCOELASTIC BEHAVIOR OF A CHARGE OF A SOLID PROPELLANT UNDER ACOUSTIC INSTABILITY OF AN INTER CHAMBER PROCESS. COMPUTATIONAL EXPERIMENT
For the first time, a intermodulation effect of the first kind (nonlinear superposition of gas oscillations with a small frequency difference) in the combustion chamber of a SPRM for a special purpose was obtained by direct numerical simulation.
The estimation of the influence of the acoustic instability of the intracameric process of a special-purpose rocket engine, taking into account the spatiotemporal localization of the dynamic load, on the vibrodynamic strength of the charge of the TT for its temperature-time viscoelastic dependence of the elastic modulus (Maxwell's relaxation kernel) and rheology is determined for the first time. The developed mathematical apparatus, the complex of applied programs created on its basis, the developed method of calculation, the conducted methodical investigations make it possible to significantly improve reliability, improve the energy mass, strength, operational and other characteristics of the SPRM.
Keywords: rocket engine of solid fuel, acoustic instability, mathematical modeling, computational experiment, viscoelastic behavior, dynamic strength, end combustion, combustion chamber, finite element method, determining relations.
Введение
Генерация акустических волн вихрями и обратное влияние акустики на образование вихрей приводят к возбуждению автоколебательного процесса в потоке. Энергия для этих колебаний берется из основного установившегося потока и ограничена нелинейными механизмами. Наличие обратной связи «вихри - звук» приводит к тому, что поток становится глобально неустойчивым. Сложная структура гомогенного потока в камере сгорания РДТТ специального
назначения обусловливает появление гидродинамических источников звука. Неустойчивость рабочих процессов в камере сгорания РДТТ объясняется неустойчивостью относительно первых продольных мод колебаний. Сложная геометрия камеры сгорания обусловливает вихревой поток продуктов сгорания. Выявленная крупномасштабность вихрей указывает на низкочастотный характер излучения звука в газодинамическом тракте двигателя. Резонансное взаимодействие акустического и вихревого полей предопределяет развитие акустической неустойчивости. Камера сгорания играет роль резонатора для колебаний звуковой природы и определяет обратную связь между вихревыми и акустическими колебаниями.
Неустойчивость процесса течения продуктов сгорания в РДТТ [1, 2, 4, 5, 7] определяется автоколебательным процессом изменения рабочих параметров, количественные показатели которых могут не соответствовать установленному диапазону изменения. Неустойчивость процесса течения в РДТТ способствует возникновению демаскирующих шумов. Камера сгорания РДТТ, как любой газовый объем, обладает характерными колебательными модами. В совокупности с прогрессивной по давлению плотностью тока продуктов сгорания газодинамическая неустойчивость приводит к установлению в камере сгорания автоколебательного процесса с конечной амплитудой.
Твердому топливу присуща явная физическая нелинейность механического поведения [3]. Исследование механических свойств ТТ показывает, что основной вклад в изменение вязкоупругих свойств материала наполнителя вносят максимальные деформации (напряжения) за время предыстории нагружения.
Для описания вязкоупругого поведения заряда ТТ при акустической неустойчивости внутрикамерного процесса использовано уравнение импульса с учетом конечно-элементной пространственной аппроксимации.
Левая часть системы уравнений (1) представляет собой согласованную матрицу масс, которая не является диагональной. Таким образом, для наилучшего численного разрешения системы (1) относительно импульсных нагрузок приведена диагонализация матрицы масс. Следует отметить еще одно допущение: тензор напряжений является известной величиной при аппроксимации по пространству, однако его можно записать через смещения и тем самым связать уравнение. Данное действие приводит к формированию глобальных матриц жесткости.
Для численного интегрирования системы (1) использовался восьмиузловой шестигранный элемент. Функция формы Ф . для данного элемента имеет вид
(1)
(2)
где ■, п ], С] принимают значения в своих узлах ±1. Силовые граничные условия записаны в виде
к =Т(г,^ ^{X
(3)
( -а- ) = 0.
Интегрирование уравнения импульса по времени
Запишем уравнение количества движения в следующем виде:
Мхп = <2" - ¥п,
(4)
где М — матрица масс; Qn — внешние нагрузки и объемные силы; ¥п — внутренние силы, определенные тензором напряжения. При этом для распределенных сил - это вклад (сумма) сил от ячеек, окружающих соответствующий узел, для сосредоточенных сил - узловая сила.
Аппроксимация системы уравнений (1) выполнена с помощью центральной конечно-разностной схемы «крест» [8]:
0п т^п ,1 1 ,1 ,1,1
. — Г. п +— п— п+— п +— п+—
Х; =-!-, х 2 = х 2 +Мп х", X,. 2 = Хп + X 2 М 2.
'
М:
(5)
Приведенная схема интегрирования имеет второй порядок точности по времени и обладает условной устойчивостью при выполнении условия Куранта т < к-^-.
Определяющие соотношения
Тензор напряжений Эйлера (Коши) для упругопластических материалов по времени интегрируется по времени следующим образом:
о а (г + йг) = о а (г) + о айг. Производная по времени может быть записана в следующем виде:
(6)
О =°ц +о,к Ц к + о аЦ.
(7)
Здесь Цк - тензор вращения, Цкк =
д— ду дх, дх
\
; Оц - тензор напряжений, вычислен-
а )
ный по Яуманну, с^ = Ст ек1; Сук1 - тензор 4-го ранга упругих постоянных, С т = о у ■■гш ; е,
ук1 ~ "у '-кП
тензор скоростей деформаций, е у =
дх; дх,
V 1 '
Тензоры скоростей деформации и вращения вычисляются в заданных точках элемента с использованием терминов матрицы деформация - смещение:
К = у к.
дх, к=1 дх, '
(8)
Таким образом, для вычисления тензора напряжений использованы следующие разностные формулы для схемы «крест»:
Щ п+1 ] п+1 п+—
о 1п+1 = о 1п + Г'1п + о'1 ^ 2^, оДг 2 = Сцк1 Дек1 2,
(9)
1 1 1 п+— п+— п+— Д£ 2 = £ 2 Дг 2 гп =
к1 Ьк1 ш ' 'а
1
п+-
опю ■2 + оIю.
'Р Р1 1Р Р'
V
1\ 1 п+—
Дг 2,
В любой момент времени напряжение может быть связано с деформацией интегралом свертки:
( й е ^
о(г) = | Е(г —;) •й£ йI, о V й ^ )
(10)
где о - напряжение; е - полная деформация (включает температурные деформации); - псевдовремя.
с
зв
Уравнение состояния применительно к режимам мгновенного и скоростного активного нагружений можно принять
для характеристики процессов релаксации напряжений (действие ДТ, Р) в виде
S¡j = | Я[ги (г), (г ,
о
для процессов ползучести (массовые силы) в виде
г
ву = | П[Ои (г), (г' — x"j)ds¡j,
(11)
(12)
где В ¡у, ву - компоненты девиаторов напряжений и деформаций; Ои, Еи - интенсивность напряжений и деформаций; Я, П - ядро релаксации и ползучести.
Влияние температуры на поведение материала описывается при помощи принципа соответствия времени - температуры. Математическое выражение данного принципа имеет вид
С« (Т, г) = Оа(Го, уг),
(13)
где у - функция смещения; уг - сокращенное время.
Уравнение WLF ("^Шат8-Ьапёе1-Реггу) используется для приблизительного определения функции
1п у =
( ст л V С2 + Т
1п10, Т =Т — Т0,
(14)
где с1 и с2 - материальные постоянные; Т0 - справочная температура, обычно называемая температурой перехода в прозрачное состояние.
Температурно-временная зависимость вязкоупругого модуля ТТ (баллистит) описывается соотношением
(
Е (г) = Е^ + (о — Е^) ехр
^ аг = А
г . у
л
V ®аг У
(15)
(1 1 >
V т То у
Константы ползучести ТТ определены экспериментально.
Корпус исследуемого РДТТ представляет собой оболочку вращения типа «полукокон». В качестве внутреннего теплозащитного покрытия корпуса использованы асбестовая прорезиненная ткань и смесь резиновая. В местах соединения с выходным блоком и дном передним в корпусе размещены соответственно фланец из титанового сплава и ряд штифтов из стали. Герметичность корпуса конструктивно обеспечивается введением в конструкцию корпуса гер-мослоя из резины, наносимого на ТЗП перед намоткой. С целью предотвращения растрескивания ТЗП при нагружении корпуса внутренним давлением в зоне вершины пера фланца имеется компенсатор из фторопластовой пленки, обеспечивающий раскрепление ТЗП от органопластика.
Цилиндрическая часть корпуса рассмотрена как многослойная оболочка радиусом Я, выполненная из композиционного материала симметричной структуры (т.е. структура слоя определяется парами одинаковых слоев, уложенных под углами ±ф) [10]. При этом обобщенные жесткости пакета материалов В у = 0, ¡Ф у, за исключением В12 = В21 Ф 0. Обобщенные жесткости
пакета материалов определяются формулами [10]
0
B12 - B21 -,
i=i
k
B11 - ^h (cos4 ф. + 2E1¡|12¡ sin2 ф. + E2¡ sin4 ф. + G12i sin2 2ф. ■,
i-i
k
Xh (( + E2¡ )sin2фг.cos2фг. + £1г|112г (sin4фг. + cos4фг.)-G12sin22фi
(16)
k
B22 - ^h (( sin4 ф. + 2E1 |12¡ sin2 ф. + E2¡ cos4 ф. + G12' sin2 2ф. ),
i-1
k
B33 - X ¿1 + E2 - 2É1K2)sin2 ф. cos2 ф. + G12 cos2 2ф.
i-1
Здесь h¡ - суммарная толщина i-го пакета материалов; Е/ - модуль упругости вдоль направления армирования; ¿2 - модуль упругости поперек направления армирования; G12 - модуль
сдвига i-го слоя; |2 - коэффициент Пуассона i-го слоя; ф. - угол армирования i-го слоя.
Средние модули упругости пакета материалов определены соотношениями [10].
Т}2 Т}2 о2
Bl - B22 " Bl -
< E >- B22 < E >- B11 < G >- B22 . и - B12 . и - BH (17)
< EX >--:-, < EY >--:-, < GXY >--:-; Иху - —; иух. (17)
h h h B11 B22
Жесткостные характеристики пакета цилиндрической части корпуса при растяжении, изгибе, сдвиге определяются соотношениями [10] < Ex > F - 2кгНЕ < Ex > - обобщенная
жесткость пакета при растяжении\сжатии; < Ex > I - nrh < Ex > - обобщенная жесткость пакета при изгибе; < Gxy > F - 2nrhz < Gxy > - обобщенная жесткость пакета при сдвиге;
< Gxy > •/ - 2nr3Й£ < G > - обобщенная жесткость пакета при кручении.
Описание нестационарного прогрева, воспламенения и последующего нестационарного и турбулентного горения заряда ТТ базируется на модели Мержанова-Дубовицкого с учетом (в рамках подхода Горохова-Липанова-Русяка) влияния газовой фазы на процесс горения в конденсированной фазе (k-фазе). Заряд ТТ моделируется как твердое тело, к которому применимы известные уравнения теплопроводности и химической кинетики.
Для математического описания процесса течения в камере сгорания РДТТ специального назначения использованы подходы механики сплошных сред [9].
Воздух, газообразные продукты сгорания заряда воспламенительного устройства и заряда твердого топлива считаем гомогенной смесью. Наличие незначительного по массе количества твердой фазы в такой смеси учитывается путем корректировки значений ее показателя адиабаты k и удельной теплоемкости при постоянном давлении по известным соотношениям.
Полная (трехмерная и нестационарная) система вихревых дифференциальных уравнений газовой динамики для гомогенного потока в камере сгорания (включая камеру воспламенителя), ракетного двигателя включает: уравнения неразрывности (сохранения массы), закон сохранения импульса по осям координат, закон сохранения полной удельной энергии газа. Система уравнений газовой динамики с учетом дополнительных соотношений интегрируется численно с помощью метода Давыдова (метода крупных частиц) - метода постановки вычислительного эксперимента, детально описанного в работах [1, 2, 5-7, 11]. Область интегрирования покрывается фиксированной в трехмерном пространстве (эйлеровой) равномерной расчетной сеткой, обеспечивающей изотропность и однородность вычислительного пространства.
Шаг интегрирования системы уравнений (1) по времени, определенный критерием Куранта, соотвествует шагу интегрирования по времени системы уравнений газовой динамики.
Сеточная топология конечно-элементной модели в задаче напряженно-деформированного состояния системы «корпус-элемент» отличается от кубической регулярной сеточной дискретизации вычислительного пространства при численном интегрировании системы уравнений газовой динамики. В связи с этим проведена пространственная интерполяция динамической нагрузки - волн давления в предсопловом объеме РДТТ.
Результаты численного моделирования по распределению параметров гомогенной смеси газа в пространстве полей течений, сопровождающемуся акустической неустойчивостью рабочего процесса в камере сгорания РДТТ, детально описаны в работах [1, 2]. Некоторые из них для безразмерного времени представлены на рис. 1-3.
р
0,0 -I----1-1----1-1-
/ 25,000 50,000 75,000 100.000 125.000
Рис. 1. Изменение во времени давления в камере сгорания (фрагмент, увеличено)
Р1(-х)-РЦ+х)-Р1(-у)- Р1(+у)
Рис. 2. Изменение во времени давления в 4 точках камеры сгорания (фрагмент, увеличено)
Рис. 3. Распределение удельной теплоемкости гомогенной смеси газа у поверхности горения ТТ для фиксированного момента времени
Результаты численного моделирования вязкоупругого поведения заряда ТТ при пространственно-временной локализации динамической нагрузки, сопоставленные с полученными МКЭ результатами вычислительного эксперимента при его статическом нагружении приведены на рис. 4-6. Максимальный уровень осевых напряжений сжатия реализуется в районе фаски торцевой поверхности контактного взаимодействия заряда ТТ с ТЗП передней крышки, и при динамическом нагружении его величина составляет 4,9 МПа, при статическом - 8,4 МПа. Динамический отклик ТТ относительно максимальных осевых перемещений сжатия при выходе РДТТ на основной режим работы практически полностью соответствует уровню осевых перемещений сжатия при статическом нагружении. Такая картина деформирования ТТ объясняется генерацией высокочастотной, но малоамплитудной газодинамической неустойчивости в камере сгорания РДТТ, а также низкой энергетикой самого баллиститного топлива.
У
1.1472 Мах
1.0452 0.94314 0.8411 0.73906 0.63702 0.53497 0.43293 0.33089 0.22885 0.12681 0.12681 0.0077272 Мт
Рис. 4. Карта радиальных перемещений системы «корпус-элемент» (мм) для момента выхода РДТТ на основной режим работы
б
Рис. 5. Карта распределения параметров НДС в осевом направлении системы «корпус-элемент»: а - перемещения (мм) - статическое нагружение камерным давлением 60 атм; б - напряжения (МПа) - статическое нагружение камерным давлением 50 атм
Верификация результатов численного моделирования. Результаты численного моделирования НДС армированного корпуса РДТТ с вкладным зарядом хорошо согласуются с показаниями тензодатчиков регистрации приращения периметра (ПК1 - ПК3), заднего днища (ЗД1 -ЗД6) и удлинения цилиндрической части корпуса (ПП1, III12) при проведении ОСИ (рис. 7-8).
б
Рис. 6. Карта осевых перемещений заряда ТТ (мм): а - статическое нагружение ТТ камерным давлением 38 атм; б - динамическое нагружение ТТ (выход ДТТ на режим)
а
Рис. 7. Схема установки тензодатчиков
ПК1 ПК2-ПКЗ
а
б
Рис. 8. Диаграммы перемещений корпуса РДТТ (мм), регистрируемые тензодатчиками
в процессе проведения ОСИ: а - приращения диаметра; б - удлинения цилиндрической части
Заключение
Проведена серия вычислительных экспериментов с использованием многопотоковой обработки информации (MPI - технология) для описания вязкоупругого поведения ТТ ракетного двигателя специального назначения. В результате проведенных численных расчетов установлено следующее:
1) в камере сгорания РДТТ реализуется высокочастотная, но малоамплитудная акустическая неустойчивость рабочего процесса;
2) численно и экспериментально подтверждена работоспособность системы «корпус-элемент» РДТТ специального назначения с позиции прочности при длительном воздействии комплекса термомеханических нагрузок;
3) результаты численного моделирования хорошо верифицируются с данными экспериментальной отработки двигателя на различных режимах его эксплуатации.
4) рассмотрены варианты повышения энергетики топлива за счет уменьшения времени работы двигателя.
Полученная расчетная информация (уровня постановки вычислительного эксперимента) может быть успешно использована при проектировании и отработке новых образцов ракетной техники и технических устройств на твердом топливе с высокими энергомассовыми, прочностными, шумовыми и другими эксплуатационными характеристиками.
Библиографический список
1. Численное моделирование нестационарных и нелинейных внутрикамерных процессов при срабатывании ракетного двигателя на твердом топливе специального назначения. Ч. 1. Постановка вычислительного эксперимента / М.Ю. Егоров, С.М. Егоров, Д.М. Егоров, Р.В. Мормуль // Вестник Пермского национального исследовательского политехнического университета. Аэрокосмическая техника. -2016. - № 47. - С. 53-72.
2. Численное моделирование нестационарных и нелинейных внутрикамерных процессов при срабатывании ракетного двигателя на твердом топливе специального назначения. Ч. 2. Результаты расчета / М.Ю. Егоров, С.М. Егоров, Д.М. Егоров, Р.В. Мормуль // Вестник Пермского национального исследовательского политехнического университета. Аэрокосмическая техника. - 2017. - № 48. - С. 26-34.
3. Алемасов В.Е., Дрегалин А.Ф., Тишин А.П. Теория ракетных двигателей. - М.: Машиностроение, 1980. - 533 с.
4. Численный эксперимент в теории РДТТ / А.М. Липанов, В.П. Бобрышев, А.В. Алиев [и др.] / под ред. А.М. Липанова. - Екатеринбург: Наука, 1994. - 301 с.
5. Егоров М.Ю., Егоров Я.В., Егоров С.М. Исследование неустойчивости рабочего процесса в двухкамерном РДТТ // Изв. вузов. Авиационная техника. - 2007. - № 4. - С. 39-43.
6. Егоров М. Ю. Метод Давыдова - современный метод постановки вычислительного эксперимента в ракетном твердотопливном двигателестроении // Вестник Пермского национального исследовательского политехнического университета. Аэрокосмическая техника. - 2014. - № 37. - С. 6-70.
7. Давыдов Ю.М., Егоров М.Ю. Численное моделирование нестационарных переходных процессов в активных и реактивных двигателях / под ред. Ю.М. Давыдова; НАПН РФ. - М., 1999. - 272 с.
8. Рихтмайер Р. Д., Мортон Х. Разностные методы решения краевых задач. - М.: Мир, 1972. - 420 с.
9. Лойцянский Л.Г. Механика жидкости и газа. - Изд. 6-е, перераб. и доп. - М.: Наука, 1987. - 840 с.
10. И.Ф. Образцов, В.В. Васильев, В. А. Бунаков. Оптимальное армирование оболочек вращения из композиционных материалов. - М.: Машиностроение, 1977. - 146 с.
11. Давыдов Ю.М. Крупных частиц метод // Математика. Большой энциклопедический словарь. -Изд-е 3-е. - М: Российская энциклопедия. - 1998/2000. - С. 303-304.
References
1. Egorov M.Yu, Egorov S.M., Egorov D.M., Mormul R.V. Chislennoye modelirovaniye nestatsionar-nykh i nelineynykh vnutrikamernykh protsessov pri srabatyvanii raketnogo dvigatelya na tvordom toplive spetsialnogo naznacheniya. Chast 1. Postanovka vychislitelnogo eksperimenta [Numerical simulation of unsteady and nonlinear intrachamber processes when triggering rocket engine on solid fuel for special purpose. Part 1. Statement of computational experiment]. PNRPUAerospace Engineering Bulletin, 2016, no. 47, pp. 53-72.
2. Egorov M.Yu., Egorov S.M., Egorov D.M., Mormul R.V. Chislennoye modelirovaniye nestatsion-arnykh i nelineynykh vnutrikamernykh protsessov pri srabatyvanii raketnogo dvigatelya na tvordom toplive spetsialnogo naznacheniya. Chast 2. Rezultaty raschota [Numerical simulation of unsteady and nonlinear in-trachamber processes when triggering rocket engine on solid fuel for special purpose. Part 2. Results of the calculation]. PNRPU Aerospace Engineering Bulletin, 2017, no. 48, pp. 26-34.
3. Alemasov V.E., Dregalin A.F., Tishin A.P. Teoriya raketnykh dvigateley [The theory of rocket engines]. Moscow, Mashinostroyeniye, 1980, 533 p.
4. Lipanov A.M., Bobryshev V.P., Aliyev A.V. Chislennyy eksperiment v teorii RDTT [Numerical experiment in the theory of solid propellant rocket]. Yekaterinburg, UIF «Nauka», 1994, 301 p.
5. Egorov M.Yu., Egorov Ya.V., Egorov S.M. Issledovanie neustoychivosti rabochego protsessa v dvukhkamernom RDTT [Study of working process instability in the two chamber solid propellant rocket engine]. Izvestiya vysshikh uchebnykh zavedeniy. Russian Aeronautics, 2007, no. 4, pp. 39-43.
6. Egorov M.Y. Metod Davydova - sovremennyy metod postanovki vychislitelnogo eksperimenta v raketnom tverdotoplivnom dvigatelestroyenii [Davydov method - modern method of statement of computing experiment in solid rocket motors industry]. PNRPU Aerospace Engineering Bulletin, 2014, no. 37, pp. 6-70.
7. Davydov Yu.M., Egorov M.Yu. CHislennoye modelirovaniye nestatsionarnykh perekhodnykh protses-sov v aktivnykh i reaktivnykh dvigatelyakh [Numerical simulation of unsteady transition processes in active and jet engines]. Moscow.: NAPS Russia, 1999, 272 p.8. Rikhtmayyer R.D., Morton KH. Raznostnyye metody resheniya krayevykh zadach. M.: Mir, 1972. 420 p.
9. Loytsyanskiy L.G. Mekhanika zhidkosti i gaza. Izd. 5-ye [Mechanics of liquid and gas. 5th ed.]. Moscow, Nauka, 1978, 736 p.
10. Obraztsov I.F., Vasilyev V.V., Bunakov V.A. Optimalnoye armirovaniye obolochek vrashcheniya iz kompozitsionnykh materialov [Optimum reinforcement of shells of rotation from composite materials]. Moscow, Mashinostroyeniye, 1977, 146 p.
11. Davydov Y.M. Krupnykh chastits metod [Large particle method]. In: Mathematics. Large Encyclopedic Dictionary. Publishing is 3rd ed. Moscow: Rossiyskaya entsiklopediya, 1998/2000, pp. 303-304.
Об авторах
Егоров Михаил Юрьевич (Пермь, Россия) - доктор физико-математических наук, профессор, профессор кафедры «Высшая математика» ФПММ, ФГБОУ ВО ПНИПУ (614990, г. Пермь, Комсомольский пр., д. 29, e-mail: egorov-m-j@yandex.ru).
Мормуль Роман Викторович (Пермь, Россия) - начальник патентного бюро ПАО «Научно-производственное объединение „Искра"» (614038, г. Пермь, ул. Академика Веденеева, д. 28), аспирант кафедры «Высшая математика» ФГБОУ ВО ПНИПУ (614990, г. Пермь, Комсомольский пр., д. 29, e-mail: rmormul@ yandex.ru).
About the authors
Mikhail Yu. Egorov (Perm, Russian Federation) - Doctor of Physical and Mathematical Sciences, Professor, Department of "Higher Mathematics", Perm National Research Polytechnic University (29, Komsomol-sky av., Perm, 614990, Russian Federation, e-mail: egorov-m-j@yandex.ru).
Roman V. Mormul (Perm, Russian Federation) - Design Engineer, Chief of Patent Office PJSC Research and Production Association "Iskra" (28, Academica Vedeneeva st., Perm, 614038, Russian Federation), PhD Student, Department of Higher Mathematics, Perm National Research Polytechnic University (29, Komso-molsky av., Perm, 614990, Russian Federation, e-mail: rmormul@yandex.ru).
Получено 08.08.2018
