Транспортное машиностроение. 2023. № 7(19). С. 66-76. ISSN 2782-5957 (print) Transport Engineering. 2023. no. 7(19). P. 66-76. ISSN 2782-5957 (print)
Научная статья Статья в открытом доступе УДК 629.463.32
doi: 10.30987/2782-5957-2023-7-66-76
МОДЕЛИРОВАНИЕ ОСТАТОЧНЫХ ДЕФОРМАЦИЙ КОТЛА ВАГОНА-ЦИСТЕРНЫ НА ОСНОВЕ ТЕОРИИ ПОЛЗУЧЕСТИ
Артем Юрьевич Сурнинш, Сергей Валерьевич Беспалько2
1,2 Российский университет транспорт, Москва, Россия
1 mister.surnin@inbox.ru, https://orcid.org/0009-0006-5532-5445
2 besp-alco@yandex.ru, https://orcid.org/0000-0002-6027-6039
Аннотация
Цель исследования заключается в том, чтобы определить влияние остаточных деформаций на поведение котла при эксплуатации.
Задачи: разработать модель ползучести по теории течения и методику расчета деформаций ползучести котла железнодорожной цистерны, возникающих при длительной эксплуатации. Использованы методы теории упругости и ползучести, метод конечных элементов.
Новизна работы заключается в разработке методики определения остаточных деформаций в котле вагона-цистерны. Результатами исследования являются зависимости величины остаточных де-
Ссылка для цитирования:
Сурнин А.Ю. Моделирование остаточных деформаций котла вагона-цистерны на основе теории ползучести / А.Ю. Сурнин, С.В. Беспалько // Транспортное машиностроение. - 2023. - № 07. - С. 66-76. doi: 10.30987/27825957-2023-7-66-76.
Original article Open Access Article
MODELING OF RESIDUAL DEFORMATIONS OF A TANK CAR BOILER BASED ON CREEP THEORY
Artem Yuryevich Surnin1^, Sergey Valeryevich Bespalko2
1,2 Russian University of Transport, Moscow, Russia
1 mister.surnin@inbox.ru, https://orcid.org/0009-0006-5532-5445
2 besp-alco@yandex.ru, https://orcid.org/0000-0002-6027-6039
Abstract
The study objective is to determine the effect of residual deformations on the behavior of the boiler during operation.
Tasks: to develop a creep model based on flow theory and a methodology for calculating creep deformations of a railway tank boiler that occur during long-term operation. The methods of elasticity and creep theory, finite element method are used.
The novelty of the work is in the development of a methodology for determining residual deformations in the boiler of a tank car. The results of the
формаций от различных факторов, деформации котла в зонах опор, выборка параметров модели ползучести.
Выводы: в исследовании было смоделировано уравнение состояния ползучести, на основе которого проведены расчеты различных вариантов технических решений по уменьшению остаточных деформаций в котле. Проведено сравнение вариантов.
Ключевые слова: котел, цистерна, метод, элементы, ползучесть, состояние, теория, течение, деформация.
study are the dependences of the value of residual deformations on various factors, the deformation of the boiler in the support zones, the selection of parameters of the creep model.
Conclusions: in the study, the creep state equation is modeled, on which basis calculations of various technical solutions for reducing residual deformations in the boiler are carried out. The options are compared.
Keywords: boiler, tank, method, elements, creep, state, theory, flow, deformation.
66
© Сурнин А. Ю., Беспалько С. В., 2023
Reference for citing:
Surnin AYu, Bespalko SV. Modeling of residual deformations of a tank car boiler based on creep theory. Transport Engineering. 2023; 7: 66-76. doi: 10.30987/2782-5957-2023-7-66-76.
Введение
Перевозка грузов по железным дорогам занимает важное место в общем объеме грузооборота. Это связано в первую очередь с высокой экономической эффективностью железнодорожных цистерн. Для грузовых вагонов важным показателем является суточная доходность, зависящая от эффективности слива груза [1, 2]. Из-за длительного воздействия постоянной нагрузки от перевозимого груза в металле возникает и развивается явление ползучести. В процессе эксплуатации в нижней зоне котла накапливаются остаточные деформации, что приводит к образованию карманов (рис.1), которые препятствуют
полной разгрузке цистерны. Кроме того, растущая интенсивность использования вагонов, сокращение времени погрузки и разгрузки увеличивают интенсивность накопления деформаций, уменьшая срок службы котла [3, 4].
Чтобы уменьшить влияния остаточных деформаций на работу котла цистерны, необходимо рассмотреть поведение котла при длительном воздействии постоянных нагрузок с учетом явления ползучести. Также следует подобрать технические решения по улучшению прочностных свойств котла.
Рис. 1. Карманы, образующиеся из-за необратимых деформаций
Fig.1. Bags formed due to irreversible deformation
Методы
Формирование математической модели котла осуществляется на основе теории ползучести [5-7] и метода конечных элементов (МКЭ) [8-13].
Напряженно-деформированное состояние (НДС) котла опишем следующей системой уравнений [5, 14, 15]:
даг дт
ху
дх дт
ду
+
дту
ух
дх дт2Х дх
+
+
дОу
ду дт
+
дх дт
+ Х = 0,
уг
дг
+ У = 0,
гу
ду
+
до1 дг дих
дх ' диу
ду '
ди7
+ 2 = 0,
дх '
ди^. ди
,
Уху Ууг Угх =
ду
+
У
ди
у
дг ди7
+
+
дх ди2 ду
диу
дх дг ах = 2Сех + Хе, оу = 2Сеу + Хе, а2 = 2Сег + Хе,
= СУ: туг = & У; тгх = &Уг
1ху
ху> 'yz'
(1)
где ох,оу
®г, Тху
Т,
Ту
ух, '-хг, ^гх, ^уг, '-гу
компоненты тензора напряжений;
£х, Еу, ^г, Уху, Уух, Ухг, Угх, Ууг, Угу - компоненты тензора деформаций; Х,У,2 -проекции объемных сил на оси координат х, у, 2; Е - модуль упругости первого рода; V - коэффициент Пуассона; С - модуль
3
Т7Л7
упругости второго рода; е - объемная деформация; А - коэффициент Ляме.
Для описания явления ползучести материала примем теорию течения. Уравнение поверхности потенциала ползучести будет иметь вид:
Г=-511511-[РШ,1)Г = 0,
где: Бц - компоненты девиатора напряжений; - интенсивность скорости деформации ползучести; t - время.
Зависимость компонентов тензора скоростей деформаций ползучести от компонентов девиатора напряжений имеет вид:
3Е-
гС 3 >1 Г
(3)
где: - интенсивность напряжения; а^ -компоненты тензора напряжений; - параметр Кронекера.
Учтем в (2) следующее выражение:
2
(2)
(4)
В итоге получим:
а[ = Е((1с,С). (5)
Формулы (3) и (5) характеризуют явление ползучести по теории течения.
Рассмотрим метод конечных элементов, суть которого заключается в том, что непрерывная величина (скалярная функция или вектор) разбивается на совокупность кусочно-непрерывных функций.
Запишем полную потенциальную энергию системы:
П(?) = и(р) - Ш(р), (6)
£
х
Е
У
=
J
где: и - потенциальная энергия деформаций; Ш - работа внешних нагрузок; V - возможные перемещения.
Для решения уравнения (6) его необходимо минимизировать на множестве функций V. Функция, которая дает минимум, является искомой функцией перемещения, обозначим ее №. Для вычисления функции ш введем пробные функции (метод Релея-Ритца) :
V = В1,р = В2,... ,У = Вп, (7) где п - размерность дискретной системы.
Далее определяется комбинация, минимизирующая уравнение (6) по пространству функций В:
п
№ = (8) 1=1
и = 11 гайУ = ЦУ Щ1В1 У и]БВ] ■ (Ю)
I 1
Обозначим:
Ъ = 1рВг. (15)
Тогда (14) примет вид:
Ш = [и]т[Р], (16)
где: [Р]т = {Рг,Р2.....Рп].
Учитывая (13) и (16) в (6) и применяя условие (9), получим:
[К][и] = [Р], (17)
где: - матрица жесткости конструкции; [и] - вектор перемещений; [Р] - вектор нагрузки.
Система уравнений (17) является разрешающей. Расчет проведен в программном комплексе
где щ - узловые значения функции ш (параметры Ритца).
Условие минимума функционала П запишем в виде системы уравнений: д
— П(ш) = 0, (9)
дщ
где / = 1, ..., п.
Для вычисления энергии деформации и в (6) необходимо отыскать интеграл произведения функций напряжения и деформации. Для этого введем оператор дифференцирования Ь, переводящий перемещения в деформации, и оператор 5, в котором к оператору Ь добавлена матрица Гука:
) 2)
Индекс ] изменяется в тех же пределах, что индекс 1. Обозначим:
Кц = | (ЬВ^БВ;).
[и]т = {и1,и2, ...,ип}. (12)
Учитывая (11) и (12), запишем (10) в следующем виде: 1
(11)
и =~[и]т[К][и].
(13)
Работа внешних сил имеет вид:
Ш = I р У и^, (14)
= I рУи&>
где р - вектор, включающий в себя объемные, поверхностные и сосредоточенные силы.
2
Расчетная схема и допущения
Расчетная схема, изображенная на рис. 2, представляет собой конечно-элементную модель котла цистерны с лежнями и лапами. При формировании конечно-элементной сетки использовались тетраэдральные конечные элементы.
Опора котла на лапы и лежни смоделирована контактными конечными элементами типа «поверхность-поверхность» [16]. В качестве исходных данных принимались: изотропный материал, свойства которого соответствуют физическим ха-
рактеристикам стали марки 17ГС, уровень заполнения котла принимался равным 97 %, температура конструкции не превышала 100 °С.
Силовая схема смоделирована распределёнными на опоры силами, соответствующими реакциям опор на гидростатическую нагрузку, возникающую при наливе груза [17]. Вычисление реакций проведен отдельным расчетом (рис. 3). На рисунке конечные элементы имеют цвет, соответствующий величине реакции.
O.SOO 1.500
Рис. 2. Расчетная схема Fig. 2 The design model
П
Рис. 3. Реакция опор на вес груза Fig.3. The reaction of the supports to the mass of the goods
В табл. 1 приведены величины сил, приходящихся на опоры.
Таблица 1
Суммарные реакции опор
Table 1
Total reactions of supports
Тип опоры Суммарные реакции опор, Н
Центральные опоры 186046,02
Лежневые опоры 412509,48
Сформируем уравнение состояния ползучести согласно (3) и (5):
= C1ac2tc3,
(18)
где С1, С2, С3 - постоянные ползучести
1
(для сталей С1 > 0, С2 > 1,- < С3 < 1 ).
С1, С2, С3 найдем, используя справочные данные [18]:
$с = 2.4713 • 10-30а3-2467^-5. (19) В исследовании рассмотрены следующие технические решения:
- монтаж броневого листа в зонах опор котла;
Результаты
Предложена методика определения остаточных деформаций в котле цистерны, учитывающая явление ползучести материала. На основе методики проведен анализ зависимостей остаточных деформаций от различных факторов. Сравнение результатов проведено относительно базового расчета с толщиной листа котла 1 см без усиления шпангоутами, тягами и броневым листом.
На рис. 4, 5 изображены поля деформаций различных расчетных вариантов. Теплые оттенки цвета показывают деформации со знаком «+» относительно вертикальной оси, вектор деформаций направлен внутрь котла. Холодные оттенки показывают деформации со знаком «-» относительно вертикальной оси, вектор направлен наружу. Учитывая особенности силовой схемы, деформация в центральной зоне котла является величиной, на которую опускается котел относительно плоскости сливного прибора. В табл. 2 приведены данные о величинах деформаций и глубине кармана для рассматриваемых вариантов.
- усиление котла внутренними или внешними разрезными шпангоутами;
- усиление котла внутренними тягами;
-увеличение толщины листов котла.
На рис. 6 и 7 приведены зависимости глубины кармана от толщины листа и толщины броневого листа, соответственно. По графикам видно ожидаемое уменьшение деформаций и глубины кармана за счет возрастания жёсткости конструкции при увеличении толщины оболочки или же локального утолщения. При этом эффективность броневого листа на порядок ниже, чем увеличение толщины оболочки.
На рис. 8 приведена обобщающая гистограмма, наглядно показывающая эффективность рассмотренных решения по снижению остаточных деформаций в котле. Лучший результат показали варианты с усилением котла шпангоутами.
Увеличение толщины котла на 50 % уменьшает глубину кармана в 3 раза. Применение броневого листа в зоне лежней позволит уменьшить глубину кармана в 1,5 раза, что является менее эффективным чем увеличение толщины листов котла, но при этом более технологичным, благодаря простой технологии монтажа на уже работающие вагоны-цистерны.
Вариант Деформация в области лежневых опор, см Глубина кармана, см
Толщина листа котла 1 см (базовый расчет) 4,03 9.7
Толщина листа котла 1,1 см 2,9 7,6
Толщина листа котла 1,2 см 2,3 5,9
Толщина листа котла 1,3 см 1,9 4,5
Толщина листа котла 1,5 см 1,1 2,9
Котел с броневым листом 0,77 6,7
Котел с внутренними шпангоутами 0,2 0,8
Котел с внешними шпангоутами 0,27 0,6
Таблица 2
Величина деформации и глубина кармана
Table 2
Deformation amount and pocket depth
в)
г)
д)
Рис. 4. Деформация котла с броневым листом: а -толщина листа 1,1 см; б - толщина листа 1,2 см; в - толщина листа 1,3 см; г - толщина листа 1,5 см; д - толщина листа 1,8 см Fig.4. The deformation of a tank car shell equipped with an armor plate thickness of 1.1 cm: а - plate thickness of 1.1 cm; b - plate thickness of 1.2 cm; в - plate thickness of 1.3 cm; г - plate thickness of 1.5 cm; д - plate thickness of 1.8 cm
а) б)
Рис. 5. Деформация котла, укрепленным шпангоутами: а - внутренние шпангоуты; II - наружные разрезные шпангоуты
Fig.5. The deformation of a tank car shell with timber reinforcement. a - inside timbers; b - outside split timbers
0
1 1,1 1,2 1,3 1,4 1,5
Толщина котла, см The thickness of the tank car shell, cm
Рис. 6. Зависимость глубины кармана от толщины котла
Fig. 6. The dependence of the bag's depth on the thickness of the tank car shell
12
0
0 0,5 1 1,5 2 2,5 3 3,5
Толщина броневого листа, см The thickness of the armor plate, cm
Рис. 7. Зависимость глубины кармана от толщины броневого листа
Fig. 7. The dependence of the bag's depth on the thickness of the armor plate
Рис. 8. Эффективность рассмотренных технических решений
Fig.8. Efficiency of the considered technical construction
Заключение
Применение тяг уменьшает глубину кармана в 3 раза, но сложность установки тяг внутри котла снижает перспективу применения.
Усиление котла шпангоутами является самым эффективным решением из рассмотренных, позволяющая уменьшить
СПИСОК ИСТОЧНИКОВ
1. Герасимов К. В. Нагруженность кузова полувагона глухого типа при падении глыбы груза : специальность 05.22.07 «Подвижной состав железных дорог, тяга поездов и электрификация» : автореф. дисс. ... канд. техн. наук / Герасимов Кирилл Вячеславович; Брянский государственный технический университет. М., 2017. 20 с.
2. Маслов И. Г. Состояние котла цистерны при воздействии очага пламени в аварийной ситуации : специальность 05.22.07 «Подвижной состав железных дорог, тяга поездов и электрификация» : автореф. дисс. ... канд. техн. наук /
глубину кармана в 15 раз. При этом процесс установки внешних разрезных шпангоутов, наряду с броневым листом, технологичный.
В дальнейшем планируется разработать методику подсчета объема остающегося груза в карманах.
Маслов Илья Геннадьевич; Московский государственный университет путей сообщения Императора Николая II. М., 2016. 24 с.
3. Шишкин Д. М. Методы расчёта остаточных напряжений в поверхностно упрочнённых призматических деталях с концентраторами напряжений в условиях реологического деформирования : специальность 1.1.8 «Механика деформируемого твердого тела» : автореф. дисс. ... канд. техн. наук / Шишкин Дмитрий Михайлович; Самарский государственный технический университет. Самара, 2021. 20 с.
4. Архипов А. В., Смольянинов А.В. Концентрация напряжений в зоне геометрических отклонений формы котла цистерны. Транспорт Урала. 2007. №1 (12). С. 53-64.
5. Мейз Дж. Теория и задачи механики сплошных сред. М.: Изд-во ЛКИ, 2007. 318 с.
6. Зотеев В.Е., Макаров Р.Ю. Численный метод определения параметров модели ползучести в пределах первых двух стадий. Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2017. Т. 16, № 2.
C. 145-156.
7. Радченко В.П., Саушкин М.Н. Феноменологический метод расчета остаточных напряжений и пластических деформаций в полом поверхностно упрочненном цилиндрическом образце. ПММ. 2013. № 1(45). С. 168-177.
8. Рычков С.П. Моделирование конструкций в среде Femap with NX Nastran. М.: ДМК Пресс, 2013. 784 с.
9. Numerical methods in finite element analysis / Bathe, K.-J., Wilson E. L. NJ : Prentice Hall, 1976. 544 p.
10. Finite element procedures : 2nd edition / Bathe, K.J. NJ : Prentice Hall, 2014. 1043 p.
11. An introduction to finite element analysis / Norrie,
D. H. de Vries G. New York : Academic Press, 1978. 301 p.
12. The Finite Element Method in Structural and Coun-tinium Mechanics / Zienkiwicz O. C., Cheung Y. New-York : McGraw-Hill Book Company, 1967. 272 p.
13. Finite elements and approximation / Zienkiewicz, O.C., Morgan K. New York : Wiley, 1983. 328 p.
14. Козлов М. П. Определение напряжений в оболочках цилиндрических частей котлов цистерн. Транспорт Урала. 2009. № 3 (22). С. 71 - 76.
15. Hoger A. On the determination of residual stress in an elastic body. Journal of Elasticity. 1998;16:303-324.
16. Емельянов И.Г., Миронов В.И., Якушев А.В., Кузнецов А.В. Контактное взаимодействие котла вагона-цистерны с лежневыми опорами. Известия Самарского научного центра Российской академии наук. 2011. Т. 13, № 1(2). С. 436-439.
17. Kim S.H., Kim C.H., Oh H. [et al.] Residual stresses and viscoelastic deformation of an injection molded automotive part. Korea-Australia rheology journal. 2007;19(4). 183-190.
18. ГОСТ 5520-79. Прокат листовой из углеродистой, низколегированной и легированной стали для котлов и сосудов, работающих под давлением : технические условия. М. : Изд-во стандартов, 2003. 15 с.
REFERENCES
1. Gerasimov KV. Loading of the gondola car body when a block of cargo falls [abstract of dissertation]. [Moscow (RF)]; Bryansk State Technical University; 2017.
2. Maslov IG. The condition of the tank boiler when exposed to a flame in an emergency situation [abstract of dissertation]. [Moscow (RF)]; Moscow State University of Railway Engineering; 2016.
3. Shishkin DM. Methods of calculating residual stresses in surface-hardened prismatic parts with stress concentrators under rheological deformation [abstract of dissertation]. [Samara (RF)]; Samara State Technical University; 2021.
4. Arkhipov AV, Smolyaninov AV. Stress concentration in the zone of geometric deviations of the tank boiler. Transport of the Urals. 2007;1(12):53-64.
5. Maze J. Theory and problems of mechanics of solid media. Moscow: LKI Publishing House; 2007.
6. Zoteev VE, Makarov RYu. Numerical method for determining the parameters of the creep model within the first two stages. Vestnik of samara university. Aerospace and mechanical engineering. 2017;16(2):145-156.
7. Radchenko VP, Saushkin MN. Phenomenological method for calculating residual stresses and plastic deformations in a hollow surface-hardened cylindrical sample. The Journal of Applied Mathematics and Mechanics. 2013;1(45):168-177.
8. Rychkov SP. Modeling of structures in Femap environment with NX Nastran. Moscow: DMK Press; 2013.
9. Bathe K-J, Wilson EL. Numerical methods in finite element analysis. NJ: Prentice Hall; 1976.
10. Bathe K-J. Finite element procedures. 2nd ed. NJ: Prentice Hall; 2014.
11. Norrie DH, de Vries G. An introduction to finite element analysis. New York: Academic Press; 1978.
12. Zienkiwicz OC, Cheung Y. The finite element method in structural and countinium mechanics. New-York: McGraw-Hill Book Company; 1967.
13. Zienkiewicz OC, Morgan K. Finite elements and approximation. New York: Wiley; 1983.
14. Kozlov MP. Determination of stresses in the shells of cylindrical parts of tank boilers. Transport of the Urals. 2009;3(22):71 - 76.
15. Hoger A. On the determination of residual stress in an elastic body. Journal of Elasticity. 1998;16:303-324.
16. Emelyanov IG, Mironov VI, Yakushev AV, Kuz-netsov AV. Contact interaction of the boiler of a tank car with leaning supports. Proceedings of the Samara Scientific Center of the Russian Academy of Sciences. 2011. Vol. 13, No. 1(2). pp. 436-439.
17. Kim SH, Kim CH, Oh H. [et al.] Residual stresses and viscoe-lastic deformation of an injection molded automotive part. Korea-Australia rheology journal. 2007;19(4).183-190.
18. GOST 5520-79. Rolled sheet of carbon, low-alloy and alloy steel for boilers and vessels operating under pressure : technical conditions. M. : Publishing House of Standards, 2003. 15 p.
Информация об авторах:
Сурнин Артем Юрьевич - аспирант Российского университета транспорта (РУТ(МИИТ); тел. 89013473427.
Беспалько Сергей Валерьевич - доктор технических наук, профессор кафедры «Вагоны и вагонное хозяйство» Российского университета транспорта (РУТ(МИИТ); тел. 8-905-585-85-36.
Surnin Artem Yuryevich - Postgraduate Student of Russian University of Transport (RUT (MIIT); phone: 89013473427.
Bespalko Sergey Valeryevich - Doctor of Technical Sciences, Professor of the Department of Cars and Car Facilities at Russian University of Transport (RUT (MIIT); phone: 8-905-585-85-36.
Вклад авторов: все авторы сделали эквивалентный вклад в подготовку публикации. Contribution of the authors: the authors contributed equally to this article.
Авторы заявляют об отсутствии конфликта интересов. The authors declare no conflicts of interests.
Статья опубликована в режиме Open Access. Article published in Open Access mode.
Статья поступила в редакцию 15.05.2023; одобрена после рецензирования 15.06.2023; принята к публикации 27.06.2023. Рецензент - Антипин Д.Я., кандидат технических наук, доцент кафедры «Подвижной состав железных дорог», директор учебно-научного института транспорта Брянского государственного технического университета, член редколлегии журнала «Транспортное машиностроение».
The article was submitted to the editorial office on 15.04.2023; approved after review on 15.06.2023; accepted for publication on 27.06.2023. The reviewer is Antipin D.Ya., Candidate of Technical Sciences, Associate Professor of the Department of Railway Rolling Stock, Director of the Educational and Scientific Institute of Transport at Bryansk State Technical University, member of the Editorial Board of the journal Transport Engineering.