УДК 539.32:519.873 ГРНТИ 90.27.32:30.19.15
ВЛИЯНИЕ ТЕПЛОФИЗИЧЕСКИХ ПРОЦЕССОВ НА РАСПРЕДЕЛЕНИЕ НАПРЯЖЕНИЙ И ДЕФОРМАЦИЙ ВО ВРАЩАЮЩЕМСЯ ПОЛОМ ЦИЛИНДРЕ
С.Ю. ЖАЧКИН, доктор технических наук, профессор
«Воронежский государственный технический университет» (г. Воронеж)
Н.А. ПЕНЬКОВ, кандидат технических наук
ВУНЦВВС «ВВА имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж)
А.С. ЧЕБОТАРЕВ, кандидат физико-математических наук
ВУНЦ ВВС «ВВА имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж)
Г.И. ТРИФОНОВ
ВУНЦ ВВС «ВВА имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж)
В статье рассмотрена задача определения напряжений, возникающих во вращающемся полом цилиндре, который находится в зоне действия стационарного теплового потока. Исследована модель однородного идеально упругого тела в рамках механики сплошных сред. Для определения связи напряжений и упругих деформаций используется закон Гука. Представлен анализ распределения главных компонент напряжений по толщине полого цилиндра. Сделаны выводы о полноте определения предела прочности металлического полого цилиндра.
Ключевые слова: напряжения, температура, деформации, закон Гука, гальваническое хромовое покрытие, цилиндр.
influence of thermophysical processes on the distribution of stresses and deformations in a rotating hollow cylinder
S.YU. ZHACHKIN, Doctor of Technical sciences, Professor
Voronezh State Technical University (Voronezh)
N.A. PENKOV, Candidate of Technical sciences
MESC AF «N.E. Zhukovsky and Y.A. Gagarin Air Force Academy» (Voronezh)
A.S. CHEBOTAREV, Candidate of Physical and Mathematical sciences
MESC AF «N.E. Zhukovsky and Y.A. Gagarin Air Force Academy» (Voronezh)
G.I. TRIFONOV
MESC AF «N.E. Zhukovsky and Y.A. Gagarin Air Force Academy» (Voronezh)
The article deals with the problem of determining the stresses arising in a rotating hollow cylinder that is located in the stationary heat flow action zone. A model of a homogeneous ideally elastic body is studied in the framework of continuum mechanics. Hooke's law is used to determine the relationship between stresses and elastic deformations. The analysis of the distribution of the principal stress components over the thickness of a hollow cylinder is presented. Conclusions are drawn about the completeness of determining the strength limit of a metal hollow cylinder.
Keywords: stresses, temperature, deformations, Hooke's law, electroplated chrome coating, cylinder.
Введение. Развитие конструкций машин происходит при постоянном стремлении к увеличению их производительности, что почти всегда сопровождается повышением механической и тепловой напряженности подвижных сопряжений деталей [1]. При этом
W U
ставятся задачи достижения высокой надежности и долговечности машины, снижения ее массы, сокращения расхода дефицитных материалов. Известно, что повышение долговечности машины даже в небольшой степени ведет к значительной экономии металла, уменьшению затрат на производство запасных деталей; сокращается число и объем ремонтов, а, следовательно, увеличивается количество фактически работающих машин [2].
Поэтому проблема повышения срока эксплуатации изделий и деталей, особенно в тех случаях, где их замена или ремонт (восстановление) по различным причинам невозможны или неоправданно дорогостоящи, является актуальной. Одним из эффективных направлений решения этой проблемы является нанесение защитных износостойких покрытий с целью упрочнения и восстановления деталей (рисунок 1).
Основным источником информации о параметрах и технологических особенностях существующих способах нанесения покрытия, в частности о физико-механических и теплофизических показателях формирования покрытия, остается эксперимент. Следовательно, при планировании использования технологии большое значение приобретает построение формальной модели при широком применении математических методов планирования.
При исследовании физико-механических и теплофизических характеристик детали, после нанесения функционального покрытия, применение математического моделирования весьма трудоемкий и сложный процесс.
В настоящее время в области теплофизики и механики сплошных сред при их аналитических и экспериментальных исследованиях появилась тенденция использования вычислительных комплексов для проведения инженерных расчетов. Системы уравнений, которые помогают оценить напряжения и упругие деформации детали (закон Гука), используют в качестве основного методологического подхода.
Рисунок 1 - Способы восстановления и упрочнения деталей
Следовательно, при проведении анализа напряжений, возникающих в детали, находящейся в зоне действия стационарного теплового потока, можно выделить специфические проблемы, характеризующие конкретную деталь и её функциональный слой. В частности, появляется необходимость моделирования систем уравнений по определению предела прочности детали.
м и
В процессе деформирования в наиболее нагруженных зонах деталей машин происходит накопление микроповреждений. Следствием такого процесса является повышение интенсивности макродеформирования нагруженных участков рабочих агрегатов [3]. Деформации такого типа оказывают существенное влияние на прочность и, как следствие, выносливость деталей.
Актуальность. Расчет на прочность деталей машин методами сопротивления материалов не позволяет в достаточной степени использовать их ресурс. В рамках допущений механики сплошных сред представлен расчет распределения напряжений по толщине нагружаемой втулки, позволяющий судить о положении опасного сечения и динамике развития главных компонент тензора напряжений в нём.
В представленной статье рассмотрен процесс деформирования металлической втулки, восстановленной хромированием в гальванической ванне [4, 5]. В расчетах учтено вращение детали и действие теплового потока во время нанесения покрытия. Наличие вращения, как показывают расчеты, оказывает непосредственное влияние на распределение напряжений и деформаций в формируемом слое хрома [6-8].
Постановка задачи. В приближении трехмерного напряженного состояния рассматривается задача о толстостенном цилиндрическом резервуаре, вращающемся вокруг своей оси с постоянной угловой скоростью и подверженном воздействию стационарного теплового потока [9, 10]. Выбирается цилиндрическая система координат рвг. Ось г проходит через центр цилиндра. Для проведения анализа напряженно-деформированного состояния используется модель изотропного идеального упругого тела. Внешняя объемная сила, действующая на цилиндр - сила инерции.
Определим напряжения и деформации в упругой области. Условие упругости запишем в
виде:
-аТ = Е ,-аТ = Е (ае
-аТ = Е (аг
- v (&, - V
)) )),
(1)
))
где £г - величина малой деформации в направлении главных осей тензора напряжений, & - главные значения тензора напряжений, а - коэффициент температурного расширения, V - коэффициент Пуассона [11].
Полагаем, что при соотношении длины цилиндра к его радиусу большем или равном
51 1 > 5 | перемещением вдоль оси цилиндра можно пренебречь: ^=0. Следовательно, е2=0 и ^а )
третье уравнение из (1) даёт:
= V(о-г +ав)-аЕТ.
(2)
Выражаем а2, подставляя (1) в уравнение (2). Система (1) преобразуется к виду:
(3)
ы
Система уравнений (3) позволяет определить искомые напряжения и деформации в цилиндре. Однако, имея в своем составе два уравнения, представленная система уравнений содержит четыре неизвестных.
Выбирая неинерциальную систему отсчета, связанную с вращающимся цилиндром, уравнение равновесия будет иметь вид [12]:
й&г а г ~ав 2 ^ —- + —-- + рч г = 0,
йг г
(4)
где р - плотность материала цилиндра, кг/м3, ч - угловая скорость вращения цилиндра вокруг собственной оси, 1/с.
Геометрические уравнения деформации сплошной среды (уравнения Коши) [13]:
йи и
£г = ~Г; ев = ~ ■
йг г
(5)
где и - радиальные перемещения при деформации.
Выражая Ог и ав из системы уравнений (3) и подставляя полученные значения в уравнение равновесия, получаем обыкновенное дифференциальное уравнение относительно модуля вектора перемещения:
йи ий2и 1 йи и Р4 (1" ) йТ
£г =—; ев=--- +----- +---- г = (1 + у)а —
йг г йг г йг г Е йг
(6)
В более удобной форме для интегрирования перепишем дифференциальное уравнение в
виде:
й Г1 йМ 1 = (1 + „ )айТ -Р (1 - "2)
йг I г йг I йг Е
г.
(7)
Для определения закона изменения перемещений необходимо задать поле температур в зоне рассматриваемого цилиндра. Положим, что в рассматриваемом случае имеет место стационарный тепловой поток [14]:
Т ( г ) = Т-Ь-Тг- 1п -,
1п
а
где температура внутри цилиндра - Та , температура окружающей среды - Ть. В принятых допущениях для перемещений будем иметь:
(8)
и = -
РЧ (1 - V2) 3 а(1 + у)(Тъ -Та)
8 Е
-г -
г 1п г
41п-
Л
а(1 + V)(ТЬ -Та) 41п ^ 2
С2
(9)
Определяем значения главных компонент напряжения путем подстановки (9) в (5) и (3):
ы
э
и
( „„.2
<У, =
Е
1 - V2
рм2 (1 - V2)(3 + V) 2 «(1 + V)(ТЬ -Та)
8Е
г --
1п г -
21п-
-(1 - V) С Ч- 1+1 с ■ -(1 + V )«Т
41п
(10)
<уа =
Е
1 - V2
РМ (1 - V2)(1 + 3v) 2 «(1 + V)(Т -Та)
8Е
г --
1п г -
21п-
Л 4^2 1 «(1 - V 2 )(Т - Та) 1 + /1.4 гр
-(1 -v)с :т —--—---г-с -(1+v)«T
41п-
(11)
Для осевой компоненты напряжения, с учетом (2) получим:
1 - V2
Л
р2('-V2)1п^201'пг{Т>-Т"' с.
2 Е
'п
'п
(12)
Постоянные интегрирования С1 и С2 определяются из граничных условий. При г=гь: Ог=0, поскольку на внешней границе нет внешних воздействий. При г=га: ов=0, так как нет трещин в месте контакта деталь-покрытие.
Таким образом, для постоянных интегрирования будем иметь:
С1 =
2 (гь2 + га2)
рм2 (1 -V)(1 + 3v-(3 + v)гь2) 2Е
«(2 (1 + V) 1п гь +1 - V )(Т - Та)
1п -
дм2 (1 - V)(1 + 3v)
4Е
(13)
С2 =
га гь (1+V)
4 (1 - V)(гЬ + га2)
рм2 (1 -V)(1 + 3v-(3 + V)г2)
2Е
«(2(1 + V) 1пгъ +1 - V)(Т - Та)
1п-
(14)
Результаты вычислений. Рассчитаем величину напряжений втулки. Зададим внутренний радиус цилиндра: а=49,7-10-3 м и внешний радиус Ь=50-10-3 м. Температура у внешней границы цилиндра равна температуре электролита хромирования: Та=36 °С. Температуру на внутренней поверхности примем равной Ть=0 °С (температура окружающей среды).
=
2
г
ь
Ы
Для рассматриваемого материала втулки значения модуля Юнга: Е = 200 ГПа, коэффициента температурного расширения а = 2,410-4, коэффициента Пуассона V = 0,3.
Подставляя значения постоянных интегрирования в выражения для главных компонент тензора напряжений окончательно получим:
1
аг = 4.62-10' • г2 - 0.88-1011 • 1п г - 77 • 1011 • — -10.08 • 1011
г
аа =-0.26•Ю8 • г2 + 0.88•Ю11 • 1пг + 3.08•Ю9 •1 -6.7Ы0
а =-21.97 •Ю6 • г2 + 0.52 •Ю11 • 1п г - 4.3 • 101
(15)
На рисунке 2 показаны графики зависимости главных компонент тензора напряжений от значения радиуса втулки.
-2x10
а- (г) аЭ (г)
02 (г) - 4x10
-6x10
8x10
Рисунок 2 - График зависимости главных компонент тензора напряжений как функций радиуса цилиндра
Выводы. Анализ полученного распределения главных компонент тензора напряжений в формируемом гальваническом композитном покрытии показал, что растягивающие напряжения создает угловая составляющая напряжения. Однако, её действие может быть скомпенсировано другими компонентами напряжений. Это может быть достигнуто соответствующим выбором режимов осаждения.
Для определения предела прочности металлической втулки предложенного подхода определения внутренних напряжений недостаточно. Необходимо рассмотрение зоны пластичности нагружаемой детали.
г
м и
СПИСОК ЛИТЕРАТУРЫ
1. Крагельский И.В., Алисин В.В. Трение, изнашивание, смазка. Справочник // М.: Машиностроение, 1978. Книга 1. 400 с.
2. Пузряков А.Ф. Теоретические основы технологии плазменного напыления: учеб. пособие / А.Ф. Пузряков. М.: Изд-во МГТУ им. Н Е. Баумана, 2008. 360 с.
3. Жачкин С.Ю., Пеньков Н.А. Восстановление гидрораспределителей сельскохозяйственной техники композитным гальваническим покрытием на основе хромовой матрицы // Труды ГОСНИТИ, 2012. Т. 109. Часть 2. С. 89-93.
4. Жачкин С.Ю. Холодное гальваноконтактное восстановление деталей. Воронеж: ВГТУ, 2002. 138 с.
5. Кудрявцев H.T. Электролитические покрытия металлами // Москва: Химия. 1979. 352 с.
6. Кутателадзе С.С. Основы теории теплообмена. Изд. 5-е перераб. и доп. М.: Атомиздат, 1979. 416 с.
7. Эрден-Груз Т. Основы строения материи. М.: Мир, 1976. 488 с.
8. Васильев Л.Л., Танаева С.А. Теплофизические свойства пористых материалов. Минск: Наука и техника, 1971. 265 с.
9. Тимошенко С.П., Дж. Гудьер. Теория упругости. Москва: Наука. 1975. 576 с.
10. Беляев Н.М. Сопротивление материалов. Москва: Наука. 1965. 856 с.
11. Биргер И.А. Сопротивление материалов. Москва: Наука. 1986. 560 с.
12. Седов Л.И. Механика сплошной среды. Т. 1. М.: Наука. 1976. 536 с.
13. Ландау Л.Д., Лифшиц Е.М. Теория упругости. М.: Наука. 1988. 248 с.
14. Ильющенко А.Ф., Шевцов А.И., Оковитый В.А., Громыко Г.Ф. Процессы формирования газотермических покрытий и их моделирование // Минск: Беларус. навука, 2011. 357 с.
REFERENCES
1. Kragel'skij I.V., Alisin V.V. Trenie, iznashivanie, smazka. Spravochnik // M.: Mashinostroenie, 1978. Kniga 1. 400 p.
2. Puzryakov A.F. Teoreticheskie osnovy tehnologii plazmennogo napyleniya: ucheb. posobie / A.F. Puzryakov. M.: Izd-vo MGTU im. N.E. Baumana, 2008. 360 p.
3. Zhachkin S.Yu., Pen'kov N.A. Vosstanovlenie gidroraspredelitelej sel'skohozyajstvennoj tehniki kompozitnym gal'vanicheskim pokrytiem na osnove hromovoj matricy // Trudy GOSNITI, 2012. T. 109. Chast' 2. pp. 89-93.
4. Zhachkin S.Yu. Holodnoe gal'vanokontaktnoe vosstanovlenie detalej. Voronezh: VGTU, 2002. 138 p.
5. Kudryavcev H.T. 'Elektroliticheskie pokrytiya metallami // Moskva: Himiya. 1979. 352 p.
6. Kutateladze S.S. Osnovy teorii teploobmena. Izd. 5-e pererab. i dop. M.: Atomizdat, 1979. 416 p.
7. Erden-Gruz T. Osnovy stroeniya materii. M.: Mir, 1976. 488 p.
8. Vasil'ev L.L., Tanaeva S.A. Teplofizicheskie svojstva poristyh materialov. Minsk: Nauka i tehnika, 1971. 265 p.
9. Timoshenko S.P., Dzh. Gud'er. Teoriya uprugosti. Moskva: Nauka. 1975. 576 p.
10. Belyaev N.M. Soprotivlenie materialov. Moskva: Nauka. 1965. 856 p.
11. Birger I.A. Soprotivlenie materialov. Moskva: Nauka. 1986. 560 p.
12. Sedov L.I. Mehanika sploshnoj sredy. T. 1. M.: Nauka. 1976. 536 p.
13. Landau L.D., Lifshic E.M. Teoriya uprugosti. M.: Nauka. 1988. 248 p.
14. Il'yuschenko A.F., Shevcov A.I., Okovityj V.A., Gromyko G.F. Processy formirovaniya gazotermicheskih pokrytij i ih modelirovanie // Minsk: Belarus. navuka, 2011. 357 p.
ы
© Жачкин С.Ю., Пеньков Н.А., Чеботарев А.С., Трифонов Г.И., 2020
Жачкин Сергей Юрьевич, доктор технических наук, профессор, профессор кафедры автоматизированного оборудования машиностроительного производства, «Воронежский государственный технический университет», Россия, 394026, г. Воронеж, Московский проспект, 14, zhach@list.ru.
Пеньков Никита Алексеевич, кандидат технических наук, начальник лаборатории научно-исследовательского центра (образовательных и информационных технологий), Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А, myth_np_nikit@mail.ru.
Чеботарев Андрей Сергеевич, кандидат физико-математических наук, доцент кафедры математики, Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А.
Трифонов Григорий Игоревич, младший научный сотрудник научно-исследовательского центра (проблем применения, обеспечения и управления авиацией Военно-воздушных сил), Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А, trifonov_gi@mail.ru.
W U