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

УДК 531.58

БОТ 10.18698/2308-6033-2016-11-1554

Моделирование баллистической стойкости двойной стенки космических аппаратов при высокоскоростном соударении

© Б.Т. Добрица1, Д Б. Добрица2

1МГТУ им. Н.Э. Баумана, Москва, 105005, Россия 2ФГУП «НПО им. С.А. Лавочкина», г. Химки, Московская область, 141400, Россия

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

Ключевые слова: баллистическое уравнение, высокоскоростное соударение, космический аппарат, экранная защита, экспериментальные исследования, численное моделирование.

Введение. Защита от угрозы воздействия высокоскоростных фрагментов космического мусора и метеорных частиц на элементы космического аппарата (КА) основана на процедуре оценки и снижения риска от таких воздействий. Если уязвимые компоненты КА идентифицированы, для них находят баллистические предельные уравнения (БПУ) — аналитические уравнения, определяющие порог непробоя высокоскоростными частицами для систем или компонентов КА. Для формулировки и проверки БПУ в рамках данной процедуры проведены тестирование высокоскоростных соударений на баллистических установках, численное моделирование, использованы полуэмпирические подходы и другие аналитические методы.

Основы механики деформируемого твердого тела систематизировано изложены в работе [1]. Алгоритм численного расчета двумерных упруго-пластических течений на основе метода конечных объемов рассмотрен в [2], основы численного моделирования задач из разных областей механики деформируемого твердого тела, жидкости и газа, в том числе трехмерных, — в работе [3]. Важность использования математического моделирования физических процессов и яв-

лений, а также применения методов численного моделирования для различных прикладных задач проиллюстрировано в [4]. Опубликованы расчетно-теоретические работы, посвященные численному моделированию динамических процессов деформирования композитных материалов [5, 6], в которых приведены примеры численной реализации разработанных математических моделей: они могут служить основой для выбора рациональных конструкций защитных преград. В работах [7, 8] представлены результаты расчетно-экспериментальных исследований по высокоскоростным взаимодействиям, проведенных для отработки противометеорной защиты.

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

Разработка инженерной методики моделирования баллистической стойкости конструкций при высокоскоростном соударении. Стойкость элементов конструкции КА к воздействию высокоскоростных твердых тел оценивается по БПУ, которые связывают критический диаметр ударника со скоростью воздействия у0, его углом, плотностью материала ударника и другими параметрами в виде множества точек (ёс, у0), характеризующих границу пробоя или отсутствия пробоя преграды для наиболее распространенных материалов и конфигураций.

Для построения баллистических уравнений используются обобщенные зависимости критического размера частиц ёс от параметров конструкции стенки, значения и направления скорости столкновения у0, а также удельного веса частиц [12]:

^ = /(Н V, 0,Рр ,Р*). (1)

Здесь ёс — критический диаметр ударника, см, который вызывает поражение конструкции на скорости воздействия; ^ — толщина внешнего слоя зашиты, см; S — расстояние между бампером и задней

стенкой, см; ^ — толщина преграды, см; у0 — скорость соударения, км/с; 9 — угол между вектором скорости частицы и нормалью к поверхности; р р — удельный вес частицы, г/см3; р ^ — удельный вес

материала стенки, г/см3. Поражение конструкции определяется как перфорация задней стенки или потеря материала (отделение осколка) от задней поверхности тыльной стенки конструкции.

Оценка стойкости элементов КА без учета экранной защиты (как одиночной стенки) к пробою при воздействии на них высокоскоростных метеорно-техногенных частиц проводится по БПУ для одиночной стенки [13]:

dc =

Г г—, ^2/3 V8/19

0,106022twНВ1'4 ^ Cw

VP p V v0 cos e

(2)

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

mc =р p pdC J 6,

где mc — критическая масса пробивающей частицы, г; HB — твердость по Бринеллю; cw — скорость звука в материале преграды, км/с.

Влияние экранно-вакуумной теплоизоляции (ЭВТИ) на проникающую способность высокоскоростных частиц при одиночной стенке можно оценить по формуле для увеличения эффективной толщины преграды при наличии ЭВТИ для скорости менее 7 км/с [14]:

t = t + K admli (3)

tweff = tw + KMLI , (3)

Pw

где tw — толщина монолитной преграды, см; KMLI = 4,5 — коэффициент, учитывающий влияние ЭВТИ; ADMLI — поверхностная плотность (масса, приходящаяся на единицу площади) ЭВТИ, г/см ; pw — плотность материала преграды, г/см3.

Для расчетов параметров пробоя высокоскоростными частицами (ударниками) разнесенной конструкции (экран + стенка) применяется уравнение Кристиансена и Керра [10]. Скорости соударения в этом уравнении разбиты на три диапазона в связи с различным характером баллистических кривых, которым соответствуют определенные функциональные зависимости.

Vh

Высокоэнтальпийный интервал скоростей V0 >-:

cos e

dc = khPp13 (V0 cos e)-23 p-1/9S12 (twpw)23 s 13 + Amli , (4)

где vH = 7 км/с; kh = 1,35 при большой толщине бампера,

- tb

kh - 7,451 + 0,411.

Условие малости толщины бампера:

< 0,126. (5)

tb

tfS113

Здесь s'h - 2,072 -10-3 aw — приведенное значение предела текучести (в оригинале задано в единицах ksi); cw — предел текучести материала

1

задней стенки (преграды), МПа; AML1 - kML1mML1 ^ ^^ j2 ; кми = 1,4; mML1 — поверхностная плотность ЭВТИ, г/см2.

Баллистический диапазон скоростей v0 < —Vl . , :

0 (cos е)1,5

dc - кг (tws f + CLtbPb ) p/'V" (cos е)-116, (6)

где ki = 1,9; aL - 3,626 -10-3 sw ; vL = 3 км/с, если ЭВТИ отсутствует; vL = 2 км/с, если ЭВТИ имеется; CL = 0,37 см3/г.

Диапазон дробления -— < V <- (область линейной ре-

(cos 0)1,5 cos 0

грессии):

d I v P-V3P-V9 oV2 (t P )23 _а/з + Л ^f Vq - vl(cos 0)-1'5 1 v 4 Vh (cos 0) - VL (cos 0) ' J

+

+ka ( t^f+CLtbPb ) pp0'5 (cos е)-5/6 ( Vh (cos е)-1 - V0

vh (cos е)-1 - Vl (cos е)-1,5

(7)

где кк1 = УнЪкИ; Л/г = у112ък1.

Отметим, что применение этих и других известных соотношений для двустенных конструкций к расчету пробоя двойных стенок, находящихся на сравнительно близком расстоянии (например, для тепловых труб под обшивкой), некорректно. Эти формулы разработаны для оптимальных параметров экранной защиты, в которых расстояние между стенками S > 15^С [9], а соотношение толщины переднего бампера и критического диаметра пробивающей частицы ^¿/^С = = 0,15. ..0,20.

Несложно убедиться, что в соответствии с уравнением (4) наращивание толщины бампера двустенной конструкции эффективно только до определенного предела, определяемого соотношением

tb max = 0,126tW/3S13, (8)

следующим из условия (5).

Примеры моделирования баллистической стойкости конструкций на основе предложенной методики. Рассмотрим конфигурацию двойной стенки. Такие конструкции регулярно встречаются при исследованиях стойкости КА к метеорно-техногенным воздействиям. Зададим следующие параметры двойной стенки (конфигурация 1): tb = 2 мм; tw = 2 мм; S = 4 см; материал бампера и задней стенки — сплав АМг6.

На рис. 1 приведена баллистическая предельная зависимость (БПЗ) для данной конфигурации, построенная по уравнению Кристи-ансена и Керра, нанесены экспериментальные точки, полученные в работе [7].

4,0 I 1 ------р

3,6 —\--------

3,2 -\--------

18 Л ___

4 2,0 - --"-^^----

1,2-------—^

0 2 4 6 8 10 12 14

Скорость, км/ч

Рис. 1. БПЗ конфигурации 1 с нанесением экспериментальных точек:

1 — пробой; 2 — на грани пробоя; 3 — непробой

Согласно соотношению (8), максимальная эффективная толщина бампера для конфигурации 1 ^ тах = 0,684 мм, следовательно, дальнейшее увеличение его толщины не приводит к улучшению характеристик БПЗ для высокоэнтальпийного интервала скоростей.

Предел роста толщины бампера для конфигурации 1 показан на рис. 2. Для расчета использован ударник из сплава Д16 (дюралюминий). Согласно графику на рис. 2, начиная с толщины бампера ^ = 0,7 мм кривые баллистической зависимости сливаются на высокоскоростном отрезке.

\

\

\

2 у v. , 3

0 2 4 6 8 10 12 14 Скорость, км/ч

2,2 2,0 1,8

г

| 1.6

1

£

1,2 1,0 0,8 0

1 3 5 7 9 11 13 15 17 19 Скорость, км/ч

Рис. 2. БПЗ конфигурации 1 при изменении толщины стенки бампера:

1-5 — Хъ = 0,5; 0,6; 0,7; 0,8 и 0,9 мм соответственно

Конфигурация 2, характерная для расчетов стойкости элементов КА к воздействию метеорно-техногенных частиц, имеет следующие параметры: Хъ = 1,5 мм, Хк = 1 мм, £ = 4 мм, материал бампера и задней стенки — алюминиевый сплав АД31.

Для расчета БПЗ в конфигурации 2 определим по соотношению (8) предельную толщину бампера. В данном случае рассматриваем упрощенную конструкцию, которая сводится к двум стенкам с расстоянием между ними Хъ тах = 0,2 мм, что значительно меньше имеющейся толщины 1,5 мм. Очевидно, что при увеличении толщины бампера расчет по формуле (2) для суммарной толщины бампер + задняя стенка даст значение стойкости конструкции выше, чем расчет по формуле (4) для двойной стенки на высокоэнтальпийном интервале скоростей. При произвольном увеличении толщины бампера Хъ >> Х™ БПЗ структуры будет стремиться к БПЗ одностенной конструкции, толщина которой равна Хъ + Хк. Очевидно также, что эффект бампера должен проявляться плавно при увеличении расстояния между стенками, начиная от нулевого. Если эффект бампера проявляется в расщеплении ударника на вторичное облако осколков, которые разлетаются в конический сноп [15], то площадь воздействия вторичного облака на заднюю стенку (т. е. эффект рассеяния вторичных осколков) будет увеличиваться пропорционально квадрату расстояния между стенками.

Для расчета БПЗ двухстенных конструкций с относительно близким расстоянием между стенками уравнение Кристиансена и Керра модифицировано с учетом этих соображений. Суть усовершенство-

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

Sid =

tb

V 0,126t/3 У

(9)

Затем для высокоэнтальпийного интервала скоростей Vo >

vh

cos 9

вво-

5 2

дим коэффициент пропорциональности Kprop = —^ (считаем, что со-

5id

отношение (4) корректно применять только при 5 =

Расчет БПЗ для V0 >

vh

cos 9

проводим на основе выбора при

vh

одного из трех вариантов, дающего наибольшее значение

Vo = О cos 9

диаметра пробивающей частицы: 1) модифицированное соотношение

dc = Kprop (khp-p13 (vo cos 9)"/3 p-19S12 (twpw)/3s13+Amlj ) +

Г , y/3 Y8/19

+ (1 K prop )

0,106022twHB14 /—

v0 cos 9

2) формула (2) для tb + tw (расчет как по одностенной конструкции);

3) формулы (4)-(5).

Расчет БПЗ для диапазона дробления

VL < V < VH

(cos 9)1'5 cos 9 • если используется 1-й вариант, то применение модифицированного соотношения дает формулу

dc =

K

prop

(khiPp Pb S (twPwУ s h +AMU) +

+ (1 -

f

s .__л 18/19 Л

xjt>1/4 Pwr„ \2/3

(1 - Kpop ) (VH )-2/3 0,106022^НВ1\1 ^ (cw )

P p

\

У У

x

V0 - VL (cos 9)-1,

1,5

VH (cos 9)-1 - VL (cos 9)-1,5

+ h

li (tws L + CLtbPb )x

x P-0'5 (cos 9)-5/6

Vh (cos 9)-1 - V0

VH (cos 9)-1 - VL (cos 9)"

1,5

ёс =

для 2-го варианта формула будет иметь вид

ч18/19 г

(V )

2/3 0,106022^НБ1/4

£(^

23

У0 -V (^ 0)

-1,5

Ун(^ 0)-1 0)"

1,5

+

+кй (^ + СьЬРь )р-0 5 (cos 0)-5/6

( Ун (cos 0)-1 - Уо Ун (cos 0)-1 - Уь (cos 0)-1,5,

• для 3-го варианта применяли формулу Кристиансена и Керра для двойной стенки.

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

Рис. 3. БПЗ для конфигурации 2, построенные по стандартному (1) и модифицированному (2) методам

Применение модифицированной формулы Кристиансена и Керра к конфигурации 1 дает результаты при определенном увеличении толщины стенки бампера (рис. 4). Без модификации все кривые совпадали бы с самой нижней кривой на участке У0 > 7 км/с.

Скорость, км/ч

Рис. 4. БПЗ конфигурации 1 при увеличении толщины стенки бампера по разработанному методу: 1-4 — tb = 2; 4; 6 и 8 мм соответственно

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

ЛИТЕРАТУРА

[1] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 4. Основы механики твердых сред. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2014, 624 с.

[2] Уилкинс М.Л. Расчет упругопластических течений. В кн.: Вычислительные методы в гидродинамике. Москва, Мир, 1967, с. 212-263.

[3] Wilkins M.L. Computer simulation of dynamic phenomena. Berlin-Heidelberg-New-York, Springer-Verlag, 1999, 264 p.

[4] Александров А.А., Димитриенко Ю.И. Математическое и компьютерное моделирование — основа современных инженерных наук. Математическое моделирование и численные методы, 2014, № 1 (1), с. 3-4.

[5] Димитриенко Ю.И., Димитриенко И.Д. Моделирование динамических процессов деформирования гибких тканевых композиционных материалов. Инженерный журнал: наука и инновации, 2014, вып. 5 (29). URL: http://engjournal.ru/search/author/40/page 1 .html

[6] Димитриенко Ю.И., Дзагания А.Ю., Беленовская Ю.В., Воронцова М.А. Численное моделирование проникания ударников в анизотропные упруго-пластические преграды. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2008, № 4, с. 100-117.

[7] Герасимов А.В., Пашков С.В., Христенко Ю.Ф. Защита космических аппаратов от техногенных и естественных осколков. Эксперимент и численное моделирование. Вестник Томского государственного университета. Математика и механика, 2011, № 4 (16), с. 70-78.

[8] Добрица Б.Т., Добрица Д.Б. Работоспособность защитных экранов с двойной сеткой при воздействии высокоскоростных частиц. Вестник Томского государственного университета. Математика и механика, 2015, № 4 (36), с. 64-70.

[9] Christiansen E.L. Design and performance equations for advanced meteoroid and debris shields. Int. J. Impact Engineering, 1993, vol. 14, pp. 145-156.

[10] Christiansen E.L., Kerr J.H. Ballistic Limit Equations for Spacecraft Shielding. Int. J. Impact Engineering, 2001, vol. 26, pp. 93-104.

[11] Добрица Д.Б. К вопросу расчета стойкости элементов конструкции космического аппарата при воздействии частиц космического мусора. Вестник ФГУП «НПО им. С.А. Лавочкина», 2012, № 5 (16), с. 53-58.

[12] Николаевский В.Н., ред. Высокоскоростные ударные явления. Москва, Мир, 1973, 528 с.

[13] Панасюк М.И., Новиков Л.С., ред. Модель космоса. В 2 т. Т. 2. Воздействие космической среды на материалы и оборудование космических аппаратов. Москва, кДу, 2007, 989 с.

[14] Ryan S., Christiansen E. Micrometeoroid and Orbital Debris Shield Ballistic Limit Analysis Program. NASA/TM -2009 214789.

[15] Солодов А.В., ред. Инженерный справочник по космической технике. Москва, Воениздат, 1977, 430 с.

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

Ссылку на эту статью просим оформлять следующим образом:

Добрица Б.Т., Добрица Д.Б. Моделирование баллистической стойкости двойной стенки космических аппаратов при высокоскоростном соударении. Инженерный журнал: наука и инновации, 2016, вып. 11.

http://dx.doi.org/10.18698/2308-6033-2016-11-1554

Добрица Борис Тимофеевич — канд. техн. наук, доцент кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор более 30 научных публикаций в области математики. e-mail: fs11@bmstu.ru

Добрица Дмитрий Борисович — канд. техн. наук, ведущий математик ФГУП «НПО им. С.А. Лавочкина». Автор более 30 научных публикаций в области проблематики защиты космических аппаратов от метеорно-техногенных воздействий. e-mail: dobrica@laspace.ru

Ballistic spacecraft double wall simulation at high speed collision

© B.T. Dobritsa1, D.B. Dobritsa2

:Bauman Moscow State Technical University, Moscow, 105005, Russia 2Federal State Unitary Enterprise NPO named after S.A. Lavochkin, Moscow, Region,

Khimki town, 141400, Russia

The article considers an engineering simulation method to calculate the ballistic limit depending on the double wall at the high-speed collision, which can be used in assessing space vehicles structural elements stability under the man-made meteor particles influence, looks at the method calculations results for two model variants equivalent to spacecraft design standard elements. We selected as models the fuel tank fragments with protection and a radiator cross-section with the built-in thermal tube, analyzed the proposed method application at increasing the bumper wall thickness for the fuel tank design that gave a positive result. The article shows various approaches validity while assessing spacecraft risks from the space debris or meteoroid impact damage - impact testing experimental results, numerical and engineering techniques.

Keywords: ballistic equation, high-speed collision, spacecraft, screen protection, experimental studies, numerical modeling.

REFERENCES

[1] Dimitrienko Yu.I. Mekhanika sploshnoy sredy. V 4 tomakh. Tom 4. Osnovy mekhaniki tverdykh sred [Continuum Mechanics. In 4 vols. Vol. 4. Fundamentals of solids]. Moscow, BMSTU Publ., 2014, 624 p.

[2] Wilkins M.L. Raschet uprugoplasticheskikh techeniy. Vychislitelnye metody v gidrodinamike [Calculation of elastic-plastic flow. Computational methods in fluid dynamics]. Moscow, Mir Publ., 1967, pp. 212-263 [in Russ.].

[3] Wilkins M.L. Computer simulation of dynamic phenomena. Springer-Verlag Publ., Berlin-Heidelberg-New-York, 1999, 264 p.

[4] Aleksandrov A.A., Dimitrienko Yu.I. Matematicheskoe modelirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2014, no. 1 (1), pp. 3-4.

[5] Dimitrienko Yu.I., Dimitrienko I.D. Inzhenernyy zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2014, no. 5 (29). Available at: http://engjournal.ru/search/author/40/page 1 .html

[6] Dimitrienko Yu.I., Dzaganiya A.Yu., Belenovskaya Yu.V., Vorontsova M.A. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki — Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2008, no. 4, pp. 100-117.

[7] Gerasimov A.V., Pashkov S.V., Khristenko Yu.F. Vestnik Tomskogo Gosudarstvennogo Universiteta. Matematika i Mekhanika — Bulletin of the Tomsk State University. Mathematics and Mechanics, 2011, no. 4 (16), pp. 70-78.

[8] Dobritsa B.T., Dobrica D.B. Vestnik Tomskogo Gosudarstvennogo Universiteta. Matematika i Mekhanika — Bulletin of the Tomsk State University. Mathematics and Mechanics, 2015, no. 4 (36), pp. 64-70.

[9] Christiansen E.L. Design and performance equations for advanced meteoroid and debris shields. International Journal of Impact Engineering, 1993, vol. 14, pp. 145-156.

[10] Christiansen E.L., Kerr J.H. Ballistic Limit Equations for Spacecraft Shielding.

International Journal of Impact Engineering, 2001, vol. 26, pp. 93-104.

[11] Dobritsa D.B. Vestnik FGUP NPO im. S.A. Lavochkina — Bulletin of FSUE NPO Lavochkin (Lavochkin Association), 2012, no. 5 (16), pp. 53-58.

[12] Nikolaevskiy V.N., ed. Vysokoskorostnye udarnye yavleniya [High Shock Phenomena]. Moscow, Mir Publ., 1973, 528 p.

[13] Panasyuk M.I., Novikov L.S., ed. Model kosmosa: Nauchno-informatsionnoe izdanie. V 2 tomakh. Tom 2. Vozdeystvie kosmicheskoy sredy na materialy i oborudovanie [Scientific-information publication. In 2 vols. Vol. 2. The impact of the space environment on materials and spacecraft equipment]. Moscow, KDU Publ., 2007, 989 p.

[14] Ryan S., Christiansen E. Micrometeoroid and Orbital Debris Shield Ballistic Limit Analysis Program. NASA / TM-2009-214789.

[15] Solodov A.V., ed. Inzhenernyy spravochnik po kosmicheskoy tekhnike [Engineering Guide on Space Technology]. Mosocw, Voenizdat Publ., 1977, 430 p.

Dobritsa B.T., Cand. Sci. (Eng.), Assoc. Professor, Department of Computational Mathematics and Mathematical Physics, Bauman Moscow State Technical University. Author of over 30 research publications in the field of mathematics. e-mail: fs11@bmstu.ru

Dobritsa D.B., Cand. Sci. (Eng.), Leading mathematician of Federal State Unitary Enterprise NPO. S.A. Lavochkin. Author of over 30 scientific publications in the field of spacecraft protection against orbital debris and micrometeoroids. e-mail: dobrica@laspace.ru