УДК: 62-762.445:677.017.442 DOI: 10.34759/trd-2021-116-04
К расчёту на герметичность фланцевого соединения трубопроводов с Z-образным металлическим уплотнением при воздействии внешней
осевой силы
Шишкин С.В., Бойков А.А.*
Московский авиационный институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, А-80, ГСП-3, 125993, Россия
*e-mail: a. boickov@yandex. ru
Статья поступила 17.12.2020
Аннотация
Рассматривается соединение трубопроводов с не контактирующими фланцами, которое герметизируется за счёт установки в узел Z-образного уплотнения. Ставится вопрос о необходимости расчёта соединения на герметичность в условиях нагружения продольной силой, возникающей в процессе эксплуатации объекта техники, в состав которого входит соединение. Выводятся зависимости, позволяющие рассчитать соединение, функционирующее условиях нагружения продольной силой, как на прочность, так и на герметичность.
Ключевые слова: Z- образное уплотнение, нагружение продольной силой, расчёт на герметичность.
Введение
Надёжность и эксплуатационная безопасность бортовых систем летательных аппаратов во многом обусловлена герметичностью соединений. В гидравлических и
пневматических системах трубопроводов авиационной и космической техники часто
используются фланцевые соединения трубопроводов с Z-образными
металлическими деформируемыми уплотнениями (рис. 1).
Рис. 1. Фланцевое соединение с 2-образным уплотнением Привлекательность данных соединений заключается в высокой плотности получаемого контакта между кромкой уплотнения и фланцем [1, 2, 3, 4]. При обеспечении герметичности нельзя допускать чрезмерного увеличения массы конструкции, поскольку оно приведёт к увеличению динамических и циклических нагрузок, действующих на летательный аппарат. Поэтому необходимо разработать методику, позволяющую рассчитать напряжённо-деформированное состояние деталей фланцевого соединения и получить величины, необходимые для расчёта соединения на прочность и герметичность в условиях действия внешней
эксплуатационной нагрузки
силы, действующей вдоль продольной оси
соединения.
Постановка задачи
Для расчёта на герметичность фланцевых соединений удобно использовать теорию Козени-Кармана [5], согласно которой, пространство между контактирующими поверхностями уплотнения и фланца представляется в виде эквивалентного пористого слоя, а утечка герметизируемой среды определяется по формулам [5]:
где: 0У(ж) и ру(г)—объёмная утечка соответственно жидкости и газа (отнесённая к атмосферному давлению); ^ - коэффициент динамической вязкости среды при температуре 1:; га, гь - внутренний и наружный радиус фактического контакта; р, ро -давление герметизируемой среды перед уплотнением и на выходе из него соответственно; К - коэффициент проницаемости пористого слоя; - толщина эквивалентного пористого слоя под нагрузкой.
Ь;л определяется по формуле:
где: у - угол наклона кромки уплотнения относительно оси острия; Я^, -шероховатость на контактирующих поверхностях уплотнения и фланца
^)= Пк (р _ р ). Отг)= !36 кисл (р -ро2)
Исл = (1 -е)( Щ + Я22)
в случае с Z-образным уплотнением:
К =(1 _*)( Щ + )х sin (у)
соответственно; в - относительная деформация шероховатости, которая
определяется по формуле:
в = ■
^ + Е22
где 5к - контактное сближение уплотнения и фланца за счёт деформации шероховатости, которое определяется по формуле [6]:
Я
Зк = 3,4 (Яаг + Яа2) (1)
где: q - давление на контакте уплотнения и фланца; НВ - твёрдость материала фланца.
Выражение (1) показывает, что герметичность соединения зависит в большей степени от высоты шероховатости контактирующих поверхностей уплотнения и фланцев, чем от давления на контакте. Поэтому, технологически обеспечивается высокая чистота контактирующих поверхностей фланцев и уплотнения. Однако, как видно из того же выражения (1), для осуществления расчёта на герметичность также необходимо знать и давление на контакте кромки уплотнения и фланца.
Расчёт соединения на этапе нагружения внутренним давлением рабочей среды был рассмотрен в работе [7]. Во время полёта летательного аппарата на соединение помимо внутреннего давления начинают действовать различные эксплуатационные нагрузки. Очень часто они выражены в виде силы, действующей вдоль оси соединения. Аналогично тому, что происходит с соединением на этапе подачи давления, при нагружении продольной силой, детали также испытывают дополнительные деформации и взаимные смещения [8]. Поэтому, давление на
контакте уплотнения и фланца изменяется по сравнению с тем, что было на этапе
подачи давления.
Необходимо разработать методику расчёта, позволяющую получить значение изменения давления на контакте уплотнения и фланца при нагружении продольной силой, чтобы рассчитать 5к по формуле (1), а затем утечку герметизируемой среды, после чего сравнить её с предельно допустимым значением и сделать вывод о герметичности соединения.
Расчётная модель и основные допущения Рассмотрим решение задачи при нагружении соединения растягивающей силой Тх. Будем полагать, что при растяжении узла внешнее усилие Тх передаётся к фланцам как равномерно распределённая кольцевая нагрузка на радиусе Так как разгрузка соединения является упругой, то целесообразно использовать метод суперпозиции и записывать уравнения, описывающие поведение соединения, в приращениях. При приложении внешней осевой силы болты дополнительно нагружаются усилием Ыб, их податливость служит причиной кинематического расхождения фланцев Л+, профили которых также получают приращения углов поворота 5фь Как и в выкладках, приведённых в работе [7], будем считать, что благодаря приспособляемости уплотнения его кромки не выходят из зацепления с продавленными канавками на поверхностях фланцев. Длину контакта кромок в уплотняемых стыках соединения также будем полагать неизменной 1]=сош! Решение задачи также будем искать в приращениях контактных усилий, используя принцип суперпозиции.
Методология исследования.
Вывод основных уравнений, описывающих поведение соединения.
Запишем уравнения равновесия деталей соединения вдоль оси вращения исходя из силовой схемы, приведённой на рис. 2.
П^1х = Мб
г2дР2х = г^-^Иб при этом:
Т
t = —
х 2 лг5
где: ЗР1х, 8Р2х - приращение осевой составляющей контактного усилия на первой и второй кромке уплотнения соответственно.
Имеем ситуацию, при которой неизвестных на 1 больше, чем количество
уравнений. Необходимо воспользоваться уравнением совместности осевых
перемещений точек деталей соединения [9, 10, 11]. Как в случае нагружения
соединения внутренним давлением, так и в случае нагружения его осевой силой,
происходит процесс упругой разгрузки соединения, приводящий к обратному
развороту уплотнения, кинематическому расхождению фланцев и их
дополнительному повороту. Общая принципиальная схема смещений и
соотношения между ними при нагружении соединения продольной силой Тх
аналогичны тем, что имеют место при воздействии на него внутреннего давления р
[12]. Следовательно, условие совместности осевых перемещений точек деталей
соединения при приложении к нему осевого усилия Тх определяется выражением,
6
аналогичным тому, что имело место на этапе нагружения соединения внутренним
давлением р [7]. При записи величин, как и в работе [7], будем обозначать:
индексом «1» - левый фланец, индексом «2» -- уплотнение, индексом «3» -- правый
фланец (см. рис 1). Таким образом, условие совместности перемещений точек
деталей соединения запишется следующим образом:
(Г1 - Г2) - ^ (Мз(Гх)(гб - г2) + М±(ТХ)(гб - Г!)) = Яб^б (3)
где: М1(Тх), М2(Тх), М3(Тх) - приращение осесимметричного изгибающего момента при приложении осевой силы Тх, действующего на левый фланец, уплотнение и правый фланец соответственно; 132(в) - характеристика жёсткости уплотнения, расчитанная на последней итерации вычислений, описанных в работе [13]; 13 - 3-я геометрическая характеристика сечений фланцев (одинаковая для обоих фланцев). Чтобы пользоваться данным уравнением, введём систему координат, аналогичную той, что была в работе [7]. Поэтому выберем положительные направления изгибающих моментов деталей так, чтобы они совпадали с положительными направлениями углов поворота деталей при выводе уравнения совместности перемещений в работе [7] (рис. 2). В итоге имеем 3 уравнения и 3 неизвестных. Теперь необходимо выразить изгибающие моменты, действующие на детали соединения, через силовые факторы, применённые для записи системы уравнений (2), подставить их в уравнение совместности перемещений (3) и найти неизвестные величины [14].
Вывод выражений для изгибающих моментов и решение системы уравнений
Запишем выражения для изгибающих моментов, действующих на детали
соединения, исходя из уравнений равновесия половинок колец согласно расчётной
схеме, приведённой на рис. 2.
Рис. 2. Расчётная силовая схема соединения, нагруженного осевой силой Т:
Труды МАИ. Выпуск № 116 Ьир://1гиёуша1. ги/
М2 (тх) = 8Р1ХГ? - 8Р2Хг*
Мг(Тх) = ^Гб + 8Р1хт2 - 1хг52 - 8QltrsH* - Зтиг5 - М1К(Тх) (4) Мз(Тх) = ^ + 8Р2хг22 - гхг} - SQзtrsH* - 6тз& - МзК(Тх)
где: 5шц, М1я(Тх) - приращения реактивной перерезывающей силы и
распределённого осесимметричного момента, действующих со стороны трубы, и момента, действующего вследствие изгиба болтов, соответственно. Следует отметить, что в записи выражений (4) М^Тх) не выглядят как функциональные зависимости изгибающих моментов от Тх. Однако, далее будет показано, что в конечном итоге изгибающие моменты являются функциями от Тх.
Выведем выражения для силовых факторов и 5ш^, действующих со
стороны трубы на фланец по методике, применённой в работах [7,15], исходя из
системы уравнений совместности перемещений сечения сопряжения трубы и
фланца (рис. 3). Данную часть задачи нельзя решать исключительно в приращениях,
иначе картина напряжённо-деформированного состояния будет некорректной.
Поэтому будем записывать систему уравнений в суммарных перемещениях с
последующим их раскрытием через суммарные силовые факторы:
'щ (0) = И^(Н*)
1 1 (5)
(0) = (РФ,
где: (0) - радиальное перемещение сечения трубы в месте сопряжения с фланцем; (Н*) - радиальное перемещение точки фланца в месте сопряжения с трубой; <рТ}_ (0) - угол поворота сечения трубы в месте сопряжения с фланцем, фФ}_— угол поворота фланца.
Рис. 3. Расчётная схема к выводу выражений для реакций, действующих со стороны
трубы на фланец 10
Введём следующие величины:
Мь = М + М(р) + М(ТХ), шг = т + 6т + 8тг, Qt = Q + SQ + SQt где: М, ш, Р - изгибающий момент фланца и распределённые реактивные момент и перерезывающая сила, действующие со стороны трубы на фланец соответственно на этапе затяжки соединения; М(р), 8т, SQ - приращения этих величин при подаче внутреннего давления; М(ТХ), 6т1, SQt - приращения тех же величин при нагружении соединения осевой силой.
Выразим перемещения сечения трубы через уравнения моментной теории оболочек [16]. Радиальное перемещение сечения сопряжения запишется в виде:
WTt(х) = е-Рх(А1 х бш(рх) + А2 х соб(рх)) + Ш где IV - частное решение дифференциального уравнения (далее - ДУ) упругой линии оболочки; в - параметр оболочки. Величина IV будет определяться выражением:
IV = (р-—)
где: ц - коэффициент Пуассона, Ет - модуль упругости материала трубы.
Так как данная величина не зависит от х, вторая и третья производные радиального перемещения по х равны тем же производным при отсутствии давления и осевой силы [17, 18, 19, 20]. Следовательно, вывод выражений для А1 и А2 , будет идентичен выкладкам работы [15]. Выражения для радиального перемещения и поворота сечения сопряжения трубы с фланцем примут следующий вид:
ш (о) = — + + — (6)
тс( ) 2(320 2$3В ЕтПт ЕтПт К '
^ (0)= (7)
Выразим перемещения фланца через зависимости теории осесимметричной деформации колец К.Б. Бицено [16]:
УТф, (Я*)= ^ + (8)
' Е!л Е13 к '
<Рф = ^ (9)
С учётом того, что хс-П)_ , исходя из рис. 3, определяется выражением
и при этом аг-П(: определяется выражением
N¿3
а = — ■
г-п* м^ получим, что
Пф, (Н' )= +
Так как осевая сила Тх не добавит принципиально новых силовых факторов, влияющих на образование N1, по сравнению с этапом нагружения соединения внутренним давлением, уравнение равновесия половинки кольца вдоль оси, перпендикулярной плоскости сечения, будет выглядеть аналогично тому, что имело
место при нагружении внутренним давлением [7]. Следовательно, выражения для N1
будут аналогичны выражениям для приведённым в работе [7]:
К = р(рщ + ^)-Qtlrs
(10)
где:
Нф — Щ + Н-кк
Таким образом, выражение (8) принимает вид для 1-го и 3-го фланцев соответственно:
И^ (н*)
мнн*
Е1з
+
Ей
(н*)
Е1з
* + 2р(рНф-^)
Е11
2Qt1rs Е11
2Qt3rs Е11
(11)
С учётом выражений (6), (7), (9), (11) система уравнений (5) приводится к следующему виду:
Для 1 -го фланца:
тН , , РГ? ^хг5 _ МНН* 2р(рнф++р^) 2(2нг3 --1—^—I--:---:— —--1-----
2Р2Б 2@3О ЕТНТ ЕтНт Е13
— Qt1 —
рэ 2р2Б Е13
Е11
Е11
(12)
Для 3-го фланца:
'Г1-Г2) ^зЪ
+ +
2Р2Б 203О ЕтНт ЕтНт
+
Е1з
РЭ 202О Е13
Е1Л
Е1-1
Решая системы уравнений (12) и (13), получаем выражения для Зш^
8Qt = ФМ(ТХ) + Л^
8т, = ХМ(Тх) + Г^
(14)
(15)
где:
_ 2в2хо . х в(1 + х)р . х_ 1 + 2ри*
2кт
ЕЕ
ЕЕ
1 +
Е1Х
Л = ^ • Г = Л!
Следует отметить, что, по сравнению с выкладками работы [7], выражения для 8Qt и 8т, оказались независимыми от номера фланца. Это произошло потому, что слагаемые, связанные с влиянием внутреннего давления, сократились при выводе выражений.
Реактивный момент от изгиба болтов определяется выражением:
Мт (Тх) = х (Тх)
(16)
Поскольку, как показывает практика, реактивный момент от изгиба болтов значительно меньше суммарного реактивного момента, действующего со стороны трубы (разница - 2 порядка), будем считать, что величина дб1 на этапе нагружения
соединения осевой силой равна аналогичной величине на этапе нагружения
внутренним давлением и на этапе затяжки соединения:
= % =
Если требуется повышенная точность расчёта, то его можно производить в несколько итераций как на этапе нагружения соединения внутренним давлением, так и осевой силой. Вывод выражения для дб приведён в работе [15].
Подставим выражения (14), (15), (16) в формулы изгибающих моментов. В результате получим:
Щ (Тх) = 0^Йб +
При этом:
= -2-Ь (Г1 - ъ); = ъ (г1 - ъ);
п
0 _ 2п(гб-г1) . ^ _ г5Г1-Г^2-Л{Г5Н*-Г{Г5; 1 1+Фг3Н*+Хг3+-вб{' 1* 1+Фг3Н*+Хг3+-дб}_ '
0 _ 2^.(Гб-Г2) . ^ _ Г3Г2-Г2-Л1Г3Н*-Т1Г3. 3* 1+Фг3Н*+Хг3+-дб3* 1+Фг3Н*+Хг3+дб/_ '
Подставляем выражения для изгибающих моментов в уравнение совместности перемещений (3). В результате получаем выражение для Ыб:
Е1
(гб - Г1) + ^ (гб - т-2) -у~Ъ\ ^ (Т1 - т-2)
К -:-—-X ^
Е13 {т0£) (г1 - Г2) - - 01'(г& - Г±) - 03<(Г& - Г2)
Подставляя полученное значение Ыб в уравнения равновесия фланцев вдоль
оси х (2), получаем значения 8Р1х, 8Р2х. Таким образом, становится возможным получить изгибающие моменты, действующие на детали соединения, и, как следствие, рассчитать их на прочность.
Суммарное контактное давление можно рассчитать таким же образом, как и на этапе нагружения соединения внутренним давлением [7]:
8Р]Х х 8ш(у) ^ = Ц
где: - суммарное контактное давление на соответствующей кромке уплотнения; цр - суммарное контактное давление, полученное на этапе нагружения соединения внутренним давлением.
Выводы
Исследование показало, что при решении герметологической задачи подход,
связанный с рассмотрением всех деталей соединения как колец большой кривизны,
подчиняющихся зависимостям теории осесимметричной деформации К.Б. Бицено,
удобен в использовании не только на этапе нагружения соединения внутренним
давлением [7], но и на этапе его нагружения осевой силой. Полученные зависимости
дают возможность провести расчёт соединения на герметичность и при этом не
имеют таких ограничений, как например зависимости, связанные с рассмотрением
фланцев как тонких осесимметричных круглых пластин [21], поскольку далеко не
каждый фланец в машиностроении можно моделировать таким образом вследствие
большого разнообразия их форм и размеров.
16
Библиографический список
1. Хилл Р. Математическая теория пластичности. - М.: Гостехиздат, 1956. -407 с.
2. Малинин Н.Н. Прикладная теория пластичности и ползучести. - М.: Машиностроение, 1975. - 400 с.
3. Давыдов Д.В., Мяснянкин Ю.М. О внедрении тел в жёсткопластическую среду // Вестник Воронежского государственного университета: Физика. Математика. 2009. № 1. С. 94 - 100.
4. Ивлев Д.Д., Максимова Л.А. О вдавливании индентора в идеальную жесткопластическую полосу // Известия РАН. Механика твердого тела. 2000. № 3. С. 131 - 136.
5. Бабкин В.Т., Зайченко А.А., Александров В.В. Герметичность неподвижных соединений и гидравлических систем. - М.: Машиностроение, 1977. -173 с.
6. Карагельский И.В., Алисин А.В. Трение, изнашивание и смазка. - М.: Машиностроение, 1978. - 400 с.
7. Бойков А.А., Шишкин С.В. К расчёту на герметичность фланцевого соединения трубопроводов с металлическим деформируемым /-образным уплотнением // Космонавтика и ракетостроение. 2020. № 5 (116). С. 86 - 99.
8. Биргер И.А., Иосилевич Г.Б. Резьбовые и фланцевые соединения. - М.: Машиностроение, 1990. - 368 с.
9. Пхон Х.К., Сысоев Е.О., Кузнецов Е.А., Сысоев О.Е. Прогнозирование
длительной прочности фасонных деталей гидрогазовых систем по параметрам акустической эмиссии // Труды МАИ. 2020. № 110. URL: http://trudymai.ru/published.php?ID=112828. DOI: 10.34759/trd-2020-110-4
10. Аунг Ч.Т., Бабайцев А.В. Исследование влияния геометрических параметров цилиндрической оболочки под давлением зажатой между абсолютно жесткими пластинами на ширину зоны контакта // Труды МАИ. 2020. № 113. URL: http://trudymai.ru/published.php?ID= 118135. DOI: 10.34759/trd-2020-113-18
11. Коровайцева Е.А. О некоторых особенностях решения задач статики мягких оболочек вращения при больших деформациях // Труды МАИ. 2020. № 114. URL: http://trudymai.ru/published.php?ID=118881. DOI: 10.34759/trd-2020-114-04
12. Биргер И.А., Мавлютов Р.Р. Сопротивление материалов. - М.: Наука, 1986. - 560 с.
13. Шишкин С.В., Бойков А.А., Колпаков А.М. К расчёту на прочность Z-образного металлического уплотнения с врезающимися элементами, деформируемого упруго-пластически // Труды МАИ. 2019. № 109. URL: http://trudymai.ru/published.php?ID=111374. DOI: 10.34759/trd-2019-109-9
14. Sadd Martin H. Elasticity. Theory, Applications and Numerics. Burlington, USA, Elsevier Inc, 2009, 536 p.
15. Бойков А.А., Шишкин С.В. К расчёту на прочность фланцевого соединения трубопроводов с металлическим деформируемым Z-образным уплотнением // Космонавтика и ракетостроение. 2020. № 4 (115). С. 45 - 56.
16. Бояршинов С.В. Основы строительной механики машин. - М.:
Машиностроение, 1973. - 456 с.
17. Zhiber A.V., Yur'eva A.M. On a Certain Class of Hyperbolic Equations with Second-Order Integrals // Journal of Mathematical Sciences, 2021, vol. 252, pp. 168 -174. DOI: 10.1007/s 10958-020-05151-y
18. Bagderina Yu. Yu. Eigenfunctions of Ordinary Differential Euler Operators // Journal of Mathematical Sciences, 2021, vol. 252, pp. 125 - 134. DOI: 10.1007/s10958-020-05147-8
19. Chuiko S. M., Nesmelova O.V. Nonlinear boundary-value problems for degenerate differential-algebraic systems // Journal of Mathematical Sciences, 2021, vol. 252, pp. 463
- 471. DOI: 10.1007/s10958-020-05174-5
20. Гражданцева Е.Ю. Интегральное исчисление функции одной переменной.
- Иркутск: Изд-во ИГУ, 2012. - 114 с.
21. Волошин А.А., Григорьев Г.Т. Расчёт и конструирование фланцевых соединений: справочник. - Л.: Машиностроение, 1979. - 128 с.
Tightness analysis method for flange connection of pipes with metal Z-shape seal during the influence of external axial force
Shishkin S.V., Boikov A.A.*
Moscow Aviation Institute (National Research University), MAI, 4, Volokolamskoe shosse, Moscow, A-80, GSP-3, 125993, Russia *e-mail: a. boickov@yandex. ru
Abstract
The most common problem in flange connections of pipes is to ensure their tightness during the functioning. At the same time it is unacceptable to design constructions of flange connections too strong and heavy, because it can lead to increasing of loads, acting on a flying vehicle. That's why it is necessary to have a method of tightness analysis. In this article connection with Z-shape seal is estimated during the phase of loading by external axial force. The problem is - to obtain relations, which allows to calculate increment of contact pressure between edge of the seal and surface of flange. It allows to calculate total contact pressure and use it in tightness analysis of connection by means of Kozeni-Karman's theory of equivalent porous layer. According to this theory, flanges and seal are estimated as rings with high curvature and modeled by means of theory of axisymmetric deformation of K.B. Bitseno. Problem is solved by using of system of equilibrium equations of flanges and equation of displacements' compatibility of contact points of flanges and seal. Influence of the pipe on flange was taken into consideration by means of using of system of equations of displacements' compatibility, applying to the profile of connection of flange and pipe. Displacements of pipe are estimated by means of moment theory of shells, displacements of flange are estimated by means of theory of axisymmetric deformation of K.B. Bitseno. As
a result, universal tightness analysis method was obtained, which can be more useful in this case, in comparison with relations, obtained by means of theory of axisymmetric thin plates, applied to flanges, because a lot of flanges in vehicles can't be modeled as axisymmetric thin plates because of their constructions.
Keywords: Z-shape seal, loading by axial force, tightness analysis.
References
1. Khill R. Matematicheskaya teoriya plastichnosti (The mathematical theory of plasticity), Moscow, Gostekhizdat, 1956, 407 p.
2. Malinin N.N. Prikladnaya teoriya plastichnosti i polzuchesti (Engineering theory of plasticity and creeping), Moscow, Mashinostroenie, 1975, 400 p.
3. Davydov D.V., Myasnyankin Yu.M. Vestnik Voronezhskogo gosudarstvennogo universiteta: Fizika. Matematika, 2009, no. 1, pp. 94 - 100.
4. Ivlev D.D., Maksimova L.A. Izvestiya RAN. Mekhanika tverdogo tela, 2000, no. 3, pp. 131 - 136.
5. Babkin V.T., Zaichenko A.A., Aleksandrov V.V. Germetichnost' nepodvizhnykh soedinenii i gidravlicheskikh system (Tightness of static connections and hydraulic systems), Moscow, Mashinostroenie, 1977, 173 p.
6. Karagel'skii I.V., Alisin A.V. Trenie, iznashivanie i smazka (Friction, deterioration and lubrication), Moscow, Mashinostroenie, 1978, 400 p.
7. Boikov A.A., Shishkin S.V. Kosmonavtika i raketostroenie, 2020, no. 5 (116), pp. 86 -99.
8. Birger I.A., Iosilevich G.B. Rez'bovye i flantsevye soedineniya (Thread and flange connections), Moscow, Mashinostroenie, 1990, 368 p.
9. Pkhon Kh.K., Sysoev E.O., Kuznetsov E.A., Sysoev O.E. Trudy MAI, 2020, no. 110. URL: http://trudymai.ru/eng/published.php?ID= 112828. DOI: 10.34759/trd-2020-110-4
10. Aung Ch.T., Babaitsev A.V. Trudy MAI, 2020, no. 113. URL: http://trudymai.ru/eng/published.php?ID=118135. DOI: 10.34759/trd-2020-113-18
11. Korovaitseva E.A. Trudy MAI, 2020, no. 114. URL: http://trudymai.ru/eng/published.php?ID= 118881. DOI: 10.34759/trd-2020-114-04
12. Birger I.A., Mavlyutov R.R. Soprotivlenie materialov (Strength of materials), Moscow, Nauka, 1986, 560 p.
13. Shishkin S.V., Boikov A.A., Kolpakov A.M. Trudy MAI, 2019, no. 109. URL: http://trudymai.ru/eng/published.php?ID= 111374. DOI: 10.34759/trd-2019-109-9
14. Sadd Martin H. Elasticity. Theory, Applications and Numerics. Burlington, USA, Elsevier Inc, 2009, 536 p.
15. Boikov A.A., Shishkin S.V. Kosmonavtika i raketostroenie, 2020, no. 4 (115), pp. 45 -56.
16. Boyarshinov S.V. Osnovy stroitel'noi mekhaniki mashin (The basics of structure mechanics of maschines), Moscow, Mashinostroenie, 1973, 456 p.
17. Zhiber A.V., Yur'eva A.M. On a Certain Class of Hyperbolic Equations with Second-Order Integrals, Journal of Mathematical Sciences, 2021, vol. 252, pp. 168 - 174. DOI: 10.1007/s10958-020-05151-y
18. Bagderina Yu. Yu. Eigenfunctions of Ordinary Differential Euler Operators, Journal of Mathematical Sciences, 2021, vol. 252, pp. 125 - 134. DOI: 10.1007/s10958-020-05147-8
19. Chuiko S. M., Nesmelova O.V. Nonlinear boundary-value problems for degenerate differential-algebraic systems, Journal of Mathematical Sciences, 2021, vol. 252, pp. 463 -471. DOI: 10.1007/s10958-020-05174-5
20. Grazhdantseva E.Yu. Integral'noe ischislenie funktsii odnoi peremennoi (Integral analysis of one-variable functions), Irkutsk, Izd-vo IGU, 2012, 114 p.
21. Voloshin A.A., Grigor'ev G.T. Raschet i konstruirovanie flantsevykh soedinenii: spravochnik (Analysis and design of flange connections. Reference book), Leningrad, Mashinostroenie, 1979, 128 p.