Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2023. Т. 27, № 1. С. 102-118 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1955
EDN: ZHFNQY
УДК 539.374
О влиянии деформаций ползучести материала вращающегося цилиндра на последующее пластическое течение
С. В. Фирсов
Институт машиноведения и металлургии ДВО РАН Хабаровского Федерального исследовательского центра ДВО РАН, Россия, 681005, Комсомольск-на-Амуре, Металлургов, 1.
Аннотация
Изучается влияние деформаций ползучести на процесс пластического течения в материале на примере задачи вращения цилиндра с внутренней полостью (трубы), на внешнюю границу которого нанесено жесткое покрытие, предотвращающее его расширение в радиальном направлении. Задача решается в рамках теории малых деформаций. Для описания пластических свойств материала используется теория течения с ассоциированным с ним условием максимальных октаэдрических напряжений Мизеса, обобщенным на случай вязкопластического течения. Для описания вязких свойств используется широко применяемый степенной закон Нортона. В области пластического течения скорости необратимых деформаций складываются из скоростей пластической деформации и скоростей деформации ползучести. Из решения случая упругого деформирования получены зависимости для нахождения скорости вращения, при которой в материале цилиндра начнется пластическое течение. Составлена система интегро-дифференциальных уравнений для нахождения перемещений и напряжений в материале цилиндра при заданной скорости вращения и накопленных необратимых деформациях. По результатам численных расчетов получено, что наличие деформаций ползучести приводит к более позднему началу пластического течения, снижению скоростей пластических деформаций, а также к уменьшению области влияния пластического течения.
Ключевые слова: упругость, пластичность, ползучесть, вращение цилиндра, толстостенная труба, вязкопластичность, плоская деформация, малые деформации.
Получение: 20 сентября 2022 г. / Исправление: 15 февраля 2023 г. / Принятие: 17 марта 2023 г. / Публикация онлайн: 25 марта 2023 г.
Механика деформируемого твердого тела Научная статья
© Коллектив авторов, 2023 © СамГТУ, 2023 (составление, дизайн, макет) Q ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Фирсов С. В. О влиянии деформаций ползучести материала вращающегося цилиндра на последующее пластическое течение // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2023. Т. 27, № 1. С. 102-118. EDN: ZHFNQY. DOI: 10.14498/vsgtu1955. Сведения об авторе
Сергей Викторович Фирсов& https://orcid.org/0000-0001-7446-6231
младший научный сотрудник; лаб. проблем создания и обработки материалов и изделий;
e-mail: firsov.s.new@yandex.ru
Введение. На продвигающихся границах областей пластического течения в теории упругопластичности отмечается одномоментная смена в механизмах необратимых деформаций: в активном процессе — с вязкого (ползучесть) на пластический (течение); при разгрузке — наоборот. При этом вязким механизмом до начала пластического течения и при разгрузке, как правило, пренебрегают. Задача с одномоментным переходом при учете и деформаций ползучести, и деформаций пластического течения рассмотрена, по-видимому, впервые в [1]. В рамках теории больших деформаций [2] данная задача рассматривается позднее в [3-5]. В [6] отмечается, что подобная проблема возникает и в классической теории упруговязкопластичности при малых деформациях. Такое замечание делается на основе решения задачи о вращении цилиндра, изготовленного из упругопластического материала и вращающегося вокруг своей оси. Деформации ползучести и последующего вяз-копластического течения возникают в материале цилиндра за счет объемных центробежных сил инерции вращения.
Прочностные расчеты материалов вращающихся дисков и валов важны для ряда технических приложений [7] механики деформирования. Соответствующие краевые задачи для упругопластических и упруговязкопластиче-ских цилиндров и дисков рассматривались неоднократно [8-14]. Теоретически постановка таких задач обуславливает возможность получения точных решений [10,12,14-16]. Данная возможность определяется использованием кусочно-линейных пластических потенциалов, то есть классических условий пластического течения максимальных касательных напряжений (условий Треска-Сен-Венана) [10-14] или максимальных приведенных напряжений (условий Ишлинского-Ивлева) [15, 16]. Ситуация в этом вопросе такая же, как и в теории температурных напряжений в упругопластических телах [17-19].
1. Общие соотношения модели. Будем считать, что материал деформируемого тела обладает упругими, вязкими и пластическими свойствами, где под вязкими свойствами понимается накопление деформаций ползучести. При расчетах будем ограничиваться только малыми деформациями. Тогда полные деформации d можно представить суммой их обратимых (упругих) е и необратимых р составляющих:
d = е + р = 2(^и + Чти),
где и — вектор перемещений. Положим, что напряжения а зависят только от обратимых деформаций е. Для задания такой зависимости воспользуемся законом Гука
а = А Ьг(е)1 + (1)
Здесь Л и у — коэффициенты Ламе. Пластические деформации имеют место в материале в условиях принадлежности напряжений поверхности нагруже-ний (текучести) в пространстве напряжений /(а, сто) = 0 (сто — предел текучести). В условиях принятия принципа Мизеса [20] функция /(а, сто) оказывается пластическим потенциалом со следованием ассоциированного с поверхностью нагружения закона пластического течения:
" = Ж = & *> 0 (2)
В качестве условия пластического течения (поверхности нагружения) будем использовать следующее обобщение условия максимальных октаэдриче-ских напряжений Мизеса [21,22]:
Г3 1 1
¡(а,аа ) = у2(т - ЦО) •• (г - цО) - о*, т — ^ - - 1г ^ в — еР - - 1г (3)
где ц — коэффициент вязкого сопротивления пластическому течению.
Деформации ползучести в материале накапливаются схожим образом [7]:
— ^Ф(ЭД) — № (4)
<1а йТ <1а' ()
Для задания потенциала ползучести Ф воспользуемся степенным законом Нортона [23] с двумя параметрами В и п:
ёФ
ЛТ = (5)
Значение Т в соотношения (4) и (5) вводится при переходе от одномерного случая к трехмерному [7] и является некой мерой трехмерного напряженного состояния. Далее в качестве этой меры будет использоваться октаэдрическая мера напряжений Мизеса
Т2 — -т •• т.
2
2. Постановка задачи. Рассмотрим полый цилиндр, радиус внутренней полости которого К\. Боковая поверхность данного цилиндра радиусом покрыта жестким слоем, предотвращающим его радиальное расширение. Введем цилиндрическую систему координат (г,<р,х), ось г которой проходит через ось симметрии цилиндра. Деформации в цилиндре происходят за счет центробежных сил, возникающих при его вращении вокруг оси симметрии с угловой скоростью ш(Ь). Угловым ускорением пренебрегаем. Положим торцы цилиндра зафиксированными, то есть материал цилиндра находится в плоском деформированном состоянии. В этом случае компоненты вектора перемещений и тензора деформаций запишутся в виде
'иу — (т', , — 2, — 0, ^^^ — I — ,
-1 ' (6) ¿(¿хр — в^^ 1 Ррр — Т иг, — С хх 1 Рхх — 0-
Для удобства расчет будем производить в безразмерных переменных:
£ — г/В.2, и — иг/В,2, т — 1/1*, д^ — о^/ао
с безразмерными параметрами:
— Д1/Д2 — 0.2, а — Л/сто — 250, @ — — 195, X — Ва^* — 0.01, п — 4, ( — ^/М*) — 0.1.
Здесь ¿* — общая продолжительность процесса деформирования, т € [0; 1]. Для краткости безразмерные переменные в дальнейшем пишем без знака тильды. Уравнение равновесия запишется в виде
аггЛ + Г1 (^ = Ч^2, = ^ —. (7)
сто
Граничные условия задачи имеют вид
°гг(Со, т) = 0, и(1, г) = 0. (8)
3. Первоначальная упругость. Положим, что в начале деформирования напряжения малы и не возникает области пластического течения, а деформации ползучести незначительны и их можно отбросить. Иными словами, необратимые деформации р будут равны нулю. Тогда уравнение равновесия (7) с учетом закона Гука (1) и кинематических соотношений (6) приведется к виду
+ £ -£ 2и = -
е^2
а + 2,0'
Проинтегрировав данное дифференциальное уравнение, с учетом граничных условий (8) получим следующее решение задачи упругого деформирования [24]:
= 1 - е2 (2а + 30)Со4 - 2 + (а + 0)^(1 + С2) 2
В£ (а + 20)((а + 0) Со2 + 0)
= _ С2 - Со2 (2а + 30)(0(Со2 + С2) + (а + 0)^2) - 0(а + 0) ^ ^ 4£2 (а + 20)((а + 0) £2 + 0) ^
Со (а + 0) (0 - С4(2а + 0))
^ = 4£2(а + 2/?) (С2(а + 0)+0/ (9)
_ {4(2а + 30) (С2(а + 0) -0) , ^ 0{2 ((а + 0) - {2(2а + 0)) 4£2(а + 20) (£2(а + 0) + 0Г 4£2(а + 20) (£о2(а + 0) + 0) ^ ,
„ =1 а /(2а + 30) Со4 + 0 ^,,,2 ^ = 4 а + 201 (а + 0)£о2 + 0 24 Г '
С ростом угловой скорости будут увеличиваться и напряжения, и в определенный момент времени они выйдут на поверхность нагружения (выполнится условие £ = 1). Это приведет к появлению области вязкопластическо-го течения наравне с областью упругого деформирования. Из приведенного условия и полученных напряжений (9) найдем соответствующее значение угловой скорости фр. Полученная функция может иметь два локальных минимума при £ = £о и £ = 1. При некотором значении £о = Ср значения параметра нагружения в этих двух точках совпадают, то есть пластическое течение начинается одновременно на внутренней и боковой поверхностях. При меньших значениях £о < Ср пластическое течение впервые начинается на внутренней
поверхности (£ = {о ), при больших {о > Ср — на боковой (£ = 1 (а + 20 )((а + 0 ) Со2 + 0 )
ф2 =
С2 =
s»
4 ,_,
0(1 - {0)V3а2 + 6а0 + 402
(а + 20)(3а + 80 + 4^3а2 + 6а0 + 402)
40 д/3а2 + 6а0 + 402 :
2 (а + 20)((а + 0)Со2 + 0)
0(1 - Со2)((2а + 30)Со2 + 0)'
д/3а2 + 6а0 + 402 - 20 л/3а2 + 6а0 + 402 + 4а + 60 '
Со < СР
Со = Ср
Со > СР
(10)
Заметим, что при замене параметров Ламе на коэффициент Пуассона V и модуль Юнга Е:
А =
Ег/
(1 - 2г/)(1 + г/)'
^ =
Е
2(1 + )
скорости в (10) будут зависеть только от v и £о:
% (1 - ^)(£о2 + 1 - 2^)
А2 = s»
(1 - 2!/)(1 - Со2)2 VT-т+й'
(1 - V)(4 - + 4л/1 - V + г/2)
2(1 - 2г/)уТ-Т+/2 ' 4 (1 - Со2 + 1 - 2^)
(1 - 2г/)(1 - С2)((3 - 2г/)£2 + 1 - 2*)'
V1 - ^ + - 1 + 2г/ V1 - ^ + г^2 + 3 - 2 г/ '
Со < £р' Со = Ср' Со > Ср'
На рис. 1 графически представлены зависимости скорости вращения, при которой начинается пластическое течение фр, от положения внутренней полости £о при v = 28/89 (рис. 1, а), и места начала пластического течения как от £о, так и от v G (—1; 0.5) (рис. 1, b), где под «m» подразумевается случай начала пластического течения на внутренней полости £ = £о, под «oui» — на боковой поверхности £ = 1, а под «6oifa» —одновременное начало на двух поверхностях. На двух оставшихся графиках приведены распределения напряжений и перемещений в материале цилиндра при ф = фр и
Со = Д1/Д2 = 0.2' а = А/сто = 250' 0 = ^/^о = 195.
Значение угловой скорости при таких параметрах фр ~ 1.95263.
4. Пластическое (вязкопластическое) течение. При дальнейшем росте угловой скорости ф(т) > фр в материале цилиндра развивается область вязкопластического течения. Для описания механизма накопления пластических деформаций в этой области воспользуемся теорией пластического течения (2) с обобщенным потенциалом Мизеса на случай вязкопластического
фр 4
out
■ both х
7 1
1 / /
/ /
---
i !
0.2 0.4 0.6 0.8 a S Jrr ^ff ^ zz
1&
1 b-. 0.5 0
-0.5
-1 0.2
0.4
0.6
0.8
Co
1 0.8 0.6 0.4 0.2 0
both
, | , | , | , | , 1 1 1 1
■ on! -
, 1 , 1 , 1 , 1 , / in
i e
u, 10-3 0.4
0.3
0.2
0.1
0 0.2
— —1
0.8 -0.6 -0.4 -0.2 0 b
0.2 0.4 v
0.4
0.6 d
0.8
i e
Рис. 1. Деформирование упругого материала. Зависимость скорости вращения -фр, при которой начнется пластическое течение, от положения внутренней полости £ = £0 (а). Место начала пластического течения в зависимости от и v (b): in — на внутренней полости, out — на боковой поверхности, both — одновременно на внутренней полости и боковой поверхности. Распределение напряжений (с) и перемещений (d) в цилиндре в момент начала
пластического течения [Figure 1. Elastic material deformation. Variation of plastic flow onset angular velocity "фр with position of the internal cavity £ = £0 (a). Location of beginning of plastic flow depending on and v (b): in — on the inner cavity, out — on the lateral surface, both — simultaneously on the inner cavity and side surface. Distribution of stresses (c) and displacements (d) in the cylinder
at ф =
3
2
u
c
течения (3), которые примут вид [25]
Аутт — р — 1 Т - 1 )
фт — £гг — Т (2агг &<р<р ) ,
Лр^ф _ р _ 1 Т - 1 (2 )
Фт — е<р<р — т - °гг - ),
АУгх — р — 1 Т - 1 (2 )
фт — — з^ Т ( Стгг ,
1«/(„...._„. Л2 . СТ )2 , _ )2
л/2
В этом случае уравнение равновесия (7) сведется к виду
£ = ^V (orr — o4>lf) + (avv -azz) + (azz -arr) .
+C u,£ u = £ ^ + (prr - pipip) +
(11)
+ ^ + a+2ф (P+ P">*) - OW (12)
Представив левую часть данного уравнения (12) в форме [26]
-1 -2 А (1 А .
и,гг + Г и,г — Г и = — - — (ИГ) ,
аг V г аг /
проинтегрируем его и получим соотношение для нахождения компоненты вектора перемещений [27]:
и(е, г) = ^ме + ад г1 + ше, *)+с'ыс, *) —1 еф2
2 а + 20 у 4 у 8а + 20'
Д(£>т) = / ж 1(ргг(ж, г) — р^(ж, г)) Аж, (13)
Ло
/а + 0 а \
, т) = ж——— (ргг (ж, т)+ р^(ж, г)) + ———ргг (ж, гПАж. „/^0 \а + 20 а + 20 /
Ненулевые компоненты тензора напряжений примут следующий вид:
-2, о« а + 0 г 00,_1Т 12а + 30^,2
=(а+0)С1— 20С2е-2+20а+§ А — 20е-1/2— ^ е2^
^ = (а + 0)С1 + 20С2 Г2 + 20"ат~20 Л + 20С-1к —
а + 20 (14)
20 (аргг + 2(а + 0)р^) — 12а, +0 £2^2
а + 20 ^ 4 ! 4 а + 20-
_/ а а + 0 а \ 1а а2|2
=аС1+20(ат2071 — 2 — аг^М — 2 аг20е ф
Подставив полученные значения перемещения (13) и напряжений (14) в граничные условия (8), найдем
С1(т) = ((2а4+3+)20 +0 ^ — 20«а + 0X3 + 0)-1.
С2(т) = — ((2а + 30^.'+° —0'а + 0) «о¥2 + (а + 0)Й'зМ) «а + 0)Й + 0)
В итоге имеем следующие зависимости для перемещения:
1
и(е, г) = /1(е, ¿)+е-1д(е, *) — (0е + (а+0) ^-1) ад —
1 — £2 (2а + 30)Со4 — 0£2 + (а + 0)^(1 + £2) ,,2
8£ (а + 20)((а + 0) £°2 +0)
ф2 (15)
и напряжении:
— 11 - 2(£-212 - 2((а + ()^-^з -
е - е0 (2а + 3()(((Со2 + а + (а + ()е02е2) - ((а + (У1 о
4£2 (а + 2()((а + () + ()
-г
а + ( й + й ^ — 2(^+2/11 + -1/2 - 2((а + ()^т°1з-
2( ( + 2( +() ) С0(а + () ((- С4(2а + ())
(а^ + 2(а + Оррр) - тт^—, одл и2,—, , 'Л Ф2- (16)
а + 2( ^ ' 4 ' 4i2(а + 2() (£02(а + ()+()
С4(2а + 3() (С2(а + () -() ф2_ (С2 ((а + ()- {2(2а + ()) „ 4£2(а + 2() (С0(а + () + ()Ф 4£2(а + 2() (^(а + () + ()Ф '
а а + ( а
а +
а а + ( а
— Чаг^е71 - а/з - 2 - аг2(Ч +
+ 1 а / (2а + 3() СО + ( 2\Ф2
+ 4 а + 2(1 (а + ()£02 + ( ЧФ'
Используя полученные соотношения (15), (16), можно найти распределения перемещении и напряжений в материале цилиндра в заданный момент времени при известных значениях угловой скорости ф(т) и накопленных необратимых деформациях р, а с помощью соотношений (11) — прирост необратимых деформаций на следующем временном шаге. При расчетах принимаем, что коэффициент вязкости ( — 0.1, а угловая скорость
— к^Т - £ ВШ^Т), 0 < Г < V2,
|3.2, т> 1/2. ( )
Приведенные на рис. 2 результаты расчетов показывают следующее. Пластическое течение начинается у внутренней полости — о, затем у боковой поверхности — 1. При увеличении скорости вращения границы этих областей £ — т1(г) и { — т2(т) движутся навстречу друг другу и в определенный момент эти области объединяются в одну. При достижении угловой скоростью фиксированного значения Ф( ) — 3.2 уровень напряжений сокращается (уменьшается величина разницы Т - 1), что приводит к сокращению скоростей пластического течения вплоть до нуля. Хотя это и происходит практически одновременно, все же можно заметить, что скорости пластического течения сначала обращаются в нуль у боковой поверхности £ — 1, а затем граница обнуления скоростей пластического течения £ — т33 (г) быстро продвигается в направлении внутренней полости £ — £о. Когда т33 становится равным о, пластическое течение в среде полностью останавливается.
5. Учет ползучести до начала пластического течения. Рассмотрим деформирование материала, обладающего вязкими свойствами. Ползучесть, в отличие от пластического течения, начинает проявлять себя при наличии напряжений в материале, то есть еще на стадии упругого деформирования.
Е 1.01
1
0.99 1.01
0.99 1.01
0.99 1.01
0.99 1.001
0.999 1.001
pl ww Е
HSS5
0.999 0.2
u, 10-3 3
2
1
0 0.2
0.4
0.6
0.4
0.6 d
П = 0
3218
чЧЧЧЧЧЧ^
Т2 = 0.385'
т3| = 0.40754
......'....... I
: , : , : t4 = 0.4371
i Т5 = 0.49
1 t = r
0.8
0.8
е
mi
-Го2---тз
1 0.8 0.6 0.4 0.2
и-х-г
"I—-!-ГТ"
Ti | т2\ Т3 Т4 I Т5 ; i
0.3 0.36 0.42 0.48 0.54 Т b
• Grr & ^^ ........ & zz
s
_1_
0.2
р,10"3 10
0.4
0.6 c
0.8
i е
~р<р<р
о
-10
\ !
4
•
i е
0.2
0.4
0.(
e
Рис. 2. Деформирование упругопластического материала. Движение упругопластических границ в материале (b). Распределение интенсивности напряжений в определенные моменты времени (а), а также остаточные напряжения (с), перемещения (d) и необратимые
деформации (e)
[Figure 2. Deformation of an elastic-plastic material. Evolution of elastoplastic boundaries in a material (b). Distribution of stress intensity at certain points in time (a). Distribution of residual stresses (c), displacements (d) and irreversible strains (e)]
1
1
1
s
1
1
1
a
u
В этом случае необратимые деформации в материале цилиндра будут накапливаться с момента начала его вращения и их скорость накопления будет соответствовать скорости накопления деформаций ползучести. Для их нахождения воспользуемся теорией течения (4) со степенным законом Нортона (5) в качестве потенциала. После преобразований получим
dr
dp фф
d
dp Zz
: e",r _ X — n— 2 (2^rr - &ZZ),
ev _ ° V <Р X — n— 2 (f — - &ZZ),
v _ X — n— 2 1(2azz - — Vipip).
(18)
Данный случай математически схож со случаем упругопластического деформирования, то есть перемещения и напряжения в материале в текущий момент времени будут рассчитываться в соответствии с формулами (15) и (16) при известных угловой скорости ф и накопленных необратимых деформациях р, но прирост необратимых деформаций на следующем шаге будет определяться в соответствии с формулами (18).
На рис. 3 представлены результаты расчетов при росте угловой скорости согласно
1.95263 (10т - sin 20яг) , 0 ^ т ^ 0.1, 1.95263, т > 0.1.
ф(т)
Если сравнить полученные решения в момент времени т = 0.1 с аналогичными для случая упругого деформирования (рис. 1), то можно видеть, что ползучесть за время разгона успевает оказать достаточное влияние на распределение напряжений, в связи с чем их выхода на поверхность нагружения при ф(0.1) = фр не происходит. С течением времени, при зафиксированной скорости вращения ф(т) = фр напряжения продолжают релаксировать за счет значительного прироста необратимых деформаций. Это также приводит к росту перемещений в материале.
6. Учет ползучести при пластическом течении. Если же угловая скорость не остановится на значении фр и продолжит расти, то в определенный момент времени напряжения все-таки выйдут на поверхность нагруже-ния и в материале цилиндра начнется вязкопластическое течение. В области его влияния в образовании необратимых деформаций будут принимать участие уже два механизма. Положим, что скорости необратимых деформаций в области вязкопластического течения будут находиться в виде суммы скоростей деформаций ползучести и пластического течения:
dp rr
dP d
dp zz
d
_ e" + sp _
ev + ep _
_ e" + ep _
í 1 £ - 1 , Xyn-1 2
£ + 2—n ^(2arr — aVLp — oZz^
(3^—---+ X—n ^ - °rr - °zz),
í 1 — - 1 , X—n-1 2
(19)
+ X—n ^ (2^zz - @rr - ^^Lp)
V3C Е
Соответственно, получим, что перемещения и напряжения находятся в соответствии с формулами (15) и (16), а прирост необратимых деформаций
Е Е
1
0.5 0
-0.5
-1
0.2
0.4
0.6
0.8
1 е
u, 10-
0.4
0.2 -
0.2
p,io-3 0.2
0.4 0.6 0.8
c
Prr Vipip Pzz
0.1 о
-0.1
k...... V \..... \ -■ V Ч 1 ! 1 !
i , i ,
u, 10 0.6
0.4
0.2
0
1 f 0.2
' Р,Ю"3
2
----- !
■ .....
-----
\
!
2 0.4 0.6 0.8 b 3 u -
i e
0.4 0.6 0.8 1 f
d
Prr P<pip Pzz ........
1 0 -1
0.2
0.4
0.6
0.8
If
-2
_L_
0.2
0.4
0.6
f
0.8
If
Рис. 3. Деформирование вязкоупругого материала. Распределение в материале напряжений при т = 0.1 (а) и т = 1 (b), перемещений при т = 0.1 (с) и т = 1 (d) и необратимых деформаций при т = 0.1 (e) и т = 1 (f)
[Figure 3. Deformation of a viscoelastic material. Distribution of stresses at т = 0.1 (a) and t = 1 (b), displacements at т = 0.1 (с) and т = 1 (d), irreversible strains at т = 0.1 (e) and
г = 1 (f)]
определяется либо по формулам (19) в области пластического течения, либо по формулам (18) вне данных областей. Результаты подобного расчета при значениях угловой скорости, приведенных в соотношениях (17), графически показаны на рис. 4.
Как видно из представленных графиков, области пластического течения появляются позже, развиваются медленнее и сход напряжений с поверхности нагружения начинается еще до распространения пластического течения на весь материал цилиндра. При более тщательном изучении (уменьшении шага по времени) было замечено, что сход напряжений с поверхности на-гружения во второй области пластического течения у боковой поверхности начинается не со стороны упругопластической границы £ = т2(т) и не с бо-
Е
Е
а
а
а
а
а
а
w
a
u
0
e
1.01
0.99 1.01
0.99 1.0003
0.9997 1.0003
0.9997 1.0003
0.9997 1.001
pi ww S
0.999
К
= 0.3335
I I
■-Г1
i 1 1 T2 = = 0.385-
I3_
= 0.435'
1 I 1
7 i
; i ; I T4 = 0.43575
т = 1
e
i 0.8 0.6 0.4 0.2
0.3 0.36 0.42 0.48 0.54 т b
E -ar.
0.2
■ ^z
0.4
0.6
0.8
1 e
0.2
0.4
0.6
0.8
u, 10-3 6
4 -
2 -
0.2
0.4
0.6 d
0.8
i e
Рис. 4. Деформирование упругопластического материала с учетом вязких свойств. Движение упругопластических границ в материале (b). Распределение интенсивности напряжений в определенные моменты времени (а), а также остаточные напряжения (с), перемещения (d) и необратимые деформации (e)
[Figure 4. Deformation of an elastic-plastic material with viscous properties. Evolution of elastoplastic boundaries in a material (b). Distribution of stress intensity at certain points in time (a). Distribution of residual stresses (c), displacements (d) and irreversible strains (e)]
ковой поверхности £ = 1. На рис. 4 хорошо видно, что условие прекращения пластического течения £ < 1 впервые выполняется внутри области пластичности т2(г) < ( < 1. Она разбивается на две области т2(т) < ( < т4(т) и шз(т) < £ < 1. Первой исчезает область т2(г) < ( < т4(т), а затем граница т3 (г) достигает боковой поверхности £ = 1 и пропадает вторая область.
7. Сравнение двух решений с оценкой влияния первоначальной ползучести. Учет вязких свойств (ползучести) приводит к более позднему появлению области пластического течения, что наиболее отчетливо видно по движению границы т2(т). С ростом напряжений увеличивается скорость
S
1
1
1
1
E
1
1
1
c
a
u
0
e
деформаций ползучести, из-за чего увеличивается скорость релаксации напряжений. Это приводит к замедлению движения упругопластических границ, а с замедлением скорости роста угловой скорости напряжения начинают сходить с поверхности нагружения, и области пластического течения сокращаются. Однако пластическое течение заканчивается не во всей области деформирования, как было в случае упругопластического материала. Из-за перераспределения напряжений в окрестности внутренней полости поддерживается более высокий уровень напряжений. В связи с этим область пластического течения у внутренней полости, хотя и сокращается со временем, но все же сохраняется и при т = 1.
Конкурирующие интересы. Конкурирующих интересов не имею. Авторский вклад и ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.
Финансирование. Выполнено в рамках госзадания ХФИЦ ДВО РАН.
Библиографический список
1. Бегун А. С., Буренин А. А., Ковтанюк Л. В. Большие необратимые деформации в условиях изменяющихся механизмов их производства и проблема задания пластических потенциалов // ДАН, 2016. Т. 470, № 3. С. 275-278. EDN: WKDEDR. DOI: https://doi. org/10.7868/S0869565216270086.
2. Буренин А. А., Ковтанюк Л. В. Большие необратимые деформации и упругое последействие. Владивосток: Дальнаука, 2013. 312 с.
3. Begun A. S., Kovtanyuk L. V., Burenin A. A., Lemza A. O. On the mechanisms of production of large irreversible strains in materials with elastic, viscous and plastic properties // Arch. Appl. Mech., 2020. vol.90, no. 4. pp. 829-845. EDN: IPUUKS. DOI: https://doi.org/ 10.1007/s00419-019-01641-x.
4. Prokudin A. N., Firsov S. V. Antiplane strain of hardening elastoviscoplastic medium// J. Siberian Federal Univ. Math. Phys., 2018. vol.11, no. 4. pp. 399-410. EDN: XVATFR. DOI: https://doi.org/10.17516/1997-1397-2018-11-4-399-410.
5. Буренин А. А., Галимзянова К. Н., Ковтанюк Л. В., Панченко Г. Л. О согласовании механизмов роста необратимых деформаций полого шара при всестороннем сжатии // ДАН, 2018. Т. 482, №4. С. 403-406. EDN: YTJEZF. DOI: https://doi.org/10.31857/ S086956520003046-3.
6. Фирсов С. В., Прокудин А. Н., Буренин А. А. Ползучесть и пластическое течение во вращающемся цилиндре с жестким включением // Сиб. журн. индустр. матем., 2019. Т. 22, №4. С. 121-133. EDN: LZHISH. DOI: https://doi.org/10.33048/sibjim.2019. 22.412.
7. Работнов Ю. Н. Ползучесть элементов конструкций. М.: Наука, 2014. 752 с.
8. Gamer U., Lance R. H. Stress distribution in a rotating elastic-plastic tube// Acta Mech., 1983. vol.50, no. 1-2. pp. 1-8. DOI: https://doi.org/10.1007/BF01170437.
9. Gamer U., Mack W., Varga I. Rotating elastic-plastic solid shaft with fixed ends// Int. J. Eng. Sci., 1997. vol.35, no. 3. pp. 253-267. DOI: https://doi.org/10.1016/ S0020-7225(96)00085-7.
10. Antoni N. Contact separation and failure analysis of a rotating thermo-elastoplastic shrink-fit assembly// Appl. Math. Mod., 2013. vol.37, no. 4. pp. 2352-2363. DOI: https://doi. org/10.1016/j.apm.2012.05.018.
11. Mack W., Plochl M. Transient heating of a rotating elastic-plastic shrink fit // Int. J. Eng. Sci., 2000. vol.38, no. 8. pp. 921-938. DOI: https://doi. org/10.1016/S0020-7225(99) 00064-6.
12. Prokudin A. N. Exact elastoplastic analysis of a rotating cylinder with a rigid inclusion under mechanical loading and unloading // ZAMM, 2020. vol. 100, no. 3, e201900213. EDN: SYSSIM. DOI: https://doi.org/10.1002/zamm.201900213.
13. Прокудин А. Н., Буренин А. А. Упругопластическое деформирование вращающегося сплошного цилиндра из линейно-упрочняющегося материала // ПММ, 2021. Т. 85, №2. С. 172-192. EDN: EDKAMB. DOI: https://doi.org/10.31857/S0032823521020077.
14. Begun A. S., Burenin A. A., Kovtanyuk L. V., Prokudin A. N. Irreversible deformation of a rotating disk having angular acceleration// Acta Mech., 2021. vol. 232, no. 5. pp. 1917-1931. EDN: UPPSFS. DOI: https://doi.org/10.1007/S00707-021-02942-5.
15. Прокудин А. Н., Фирсов С. В. Упругопластическое деформирование вращающегося полого цилиндра с жестким внешним покрытием // Вестн. ПНИПУ. Механика, 2019. №4. С. 120-135. EDN: VUXFSF. DOI: https://doi.org/10.15593/perm.mech/2019.4.12.
16. Прокудин А. Н. Упругопластический анализ вращающегося сплошного цилиндра при условии максимальных приведенных напряжений // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2020. Т. 24, №1. С. 74-94. EDN: LJTYOU. DOI: https://doi.org/ 10.14498/vsgtu1737.
17. Буренин А. А., Ткачева А. В., Щербатюк Г. А. К расчету неустановившихся температурных напряжений в упругопластических телах // Выч. мех. спл. сред, 2017. Т. 10, №3. С. 245-259. EDN: ZHZVHJ. DOI: https://doi.org/10.7242/1999-6691/2017.10.3.20.
18. Буренин А. А., Ткачева А. В. Эволюция температурных напряжений в задаче Гадолина о сборке двухслойной упругопластической трубы // Вестн. ПНИПУ. Механика, 2020. №3. С. 20-31. EDN: HMUHTT. DOI: https://doi.org/10.15593/perm.mech/2020.3.03.
19. Буренин А. А., Ткачева А. В. Кусочно-линейные пластические потенциалы как средство расчетов плоских неустановившихся температурных напряжений // Изв. РАН. МТТ, 2020. Т. 6. С. 40-49. EDN: MFTXER. DOI: https://doi.org/10.31857/ S0572329920060057.
20. Быковцев Г. И., Ивлев Д. Д. Теория пластичности. Владивосток: Дальнаука, 1998. 528 с.
21. Мосолов П. П., Мясников В. П. Механика жесткопластических сред. М.: Наука, 1988. 208 с.
22. Ковтанюк Л. В., Шитиков А. В. О теории больших упругопластических деформаций материалов при учете температурных и реологических эффектов // Вестн. ДВО РАН, 2006. №4. С. 87-93. EDN: HZMRVR.
23. Norton F. H. The Creep of Steel at High Temperatures / Classic Reprint Series. London: Forgotten Books, 2017. 102 pp.
24. Фирсов С. В. Необратимые деформации вращающегося цилиндра // Изв. АлтГУ, 2018. Т. 102, №4. С. 114-117. EDN: YABMRN. DOI: https://doi.org/10.14258/izvasu(2018)4-21.
25. Прокудин А. Н., Фирсов С. В. Вязкопластическое течение вращающегося полого цилиндра// Дальневост. матем. журн., 2018. Т. 18, №2. С. 242-260. EDN: YUNRJB..
26. Буренин А. А., Ткачева А. В. О сборке двухслойной металлической трубы способом горячей посадки// Изв. РАН. МТТ, 2019. №3. С. 86-99. EDN: YPOULA. DOI: https://doi. org/10.1134/S0572329919030073.
27. Фирсов С. В. Пластическое течение и ползучесть в полом цилиндре с жестким внешним покрытием под действием внутреннего давления // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2021. Т. 25, №4. С. 696-715. EDN: TMRVNH. DOI: https://doi.org/10. 14498/vsgtu1877.
28. Прокудин А. Н., Фирсов С. В. Расчет ползучести вращающегося цилиндра со свободными концами// Вестн. ЧГПУ им. И.Я. Яковлева. Сер. Мех. пред. сост., 2018. №1 (35). С. 63-73. EDN: UVVBZG.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 1, pp. 102-118
d https://doi.org/10.14498/vsgtu1955
ISSN: 2310-7081 (online), 1991-8615 (print)
MSC: 74C10
The influence of creep deformations on the subsequent plastic flow in a material of rotating cylinder
S. V. Firsov
Institute of Metallurgy and Mechanical Engineering Far-Eastern Branch of RAS, Khabarovsk Federal Research Center of the Far Eastern Branch of RAS, 1, Metallurgov str., Komsomolsk-on-Amur, 681005, Russian Federation.
Abstract
The influence of creep deformations on the process of plastic flow in a material is studied by using the example of the rotational motion of a cylinder with an inner cavity (a pipe) that has a rigid coating on its outer boundary to prevent radial expansion. The problem is solved within the frameworks of the theory of infinitesimal deformations. The theory of plastic flow with the associated condition of maximum octahedral stresses of von Mises, generalized to the case of viscoplastic flow, is used to describe the plastic properties of the material. The Norton's power law is used to describe the viscous properties. In the plastic flow region, the irreversible deformation rates are composed of plastic deformation rates and creep deformation rates. The dependencies required to determine the rotational speed at which plastic deformation initiates in the cylinder material are derived from the elastic deformation solution. A system of integro-differential equations is compiled to find the displacements and stresses in the cylinder material for the specified rotational speeds and accumulated irreversible deformations. Numerical calculations show that the presence of creep deformations leads to a later initiation of plastic flow, a reduction in plastic deformation rates, and a decrease in the plastic flow influence area.
Keywords: elasticity, plasticity, creep deformation, rotating cylinder, thick-walled tube, viscoplasticity, plain strain, small strain.
Received: 20th September, 2022 / Revised: 15th February, 2023 / Accepted: 17th March, 2023 / First online: 25th March, 2023
Mechanics of Solids Research Article
© Authors, 2023
© Samara State Technical University, 2023 (Compilation, Design, and Layout) Q ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Firsov S. V. The influence of creep deformations on the subsequent plastic flow in a material of rotating cylinder, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 1, pp. 102-118. EDN: ZHFNQY. DOI: 10.14498/vsgtu1955 (In Russian). Author's Details:
Sergey V. Firsov © https://orcid.org/0000-0001-7446-6231
Junior Researcher; Lab. of the Problems of Creation and Processing of Materials and Products; e-mail: firsov.s.new@yandex.ru
Competing interests. I hereby declare that I have no conflict of interest in the authorship and publication of this article.
Author's Responsibilities. I take full responsibility for submitting the final manuscript
in print. I approved the final version of the manuscript.
Funding. The work was carried out within the framework of the state task of the KhFRC FEB RAS.
References
1. Begun A. S., Burenin A. A., Kovtanyuk L. V. Large irreversible deformations under conditions of changing mechanisms of their formation and the problem of definition of plastic potentials, Dokl. Phys., 2016, vol.61, no. 9, pp. 463-466. EDN: XFIORH. DOI: https://doi. org/10.1134/S102833581609007X.
2. Burenin A. A., Kovtanyuk L. V. Bol'shie neobratimye deformatsii i uprugoe posledeistvie [Large Irreversible Deformations and Elastic Aftereffect]. Vladivostok, Dal'nauka, 2013, 312 pp. (In Russian)
3. Begun A. S., Kovtanyuk L. V., Burenin A. A., Lemza A. O. On the mechanisms of production of large irreversible strains in materials with elastic, viscous and plastic properties, Arch. Appl. Mech., 2020, vol.90, no. 4, pp. 829-845. EDN: IPUUKS. DOI: https://doi.org/ 10.1007/s00419-019-01641-x.
4. Prokudin A. N., Firsov S. V. Antiplane strain of hardening elastoviscoplastic medium, J. Siberian Federal Univ. Math. Phys., 2018, vol.11, no. 4, pp. 399-410. EDN: XVATFR. DOI: https://doi.org/10.17516/1997-1397-2018-11-4-399-410.
5. Burenin A. A., Galimzyanova K. N., Kovtanyuk L. V., Panchenko G. L. Matching growth mechanisms of irreversible deformation of a hollow sphere under uniform compression, Dokl. Phys., 2018, vol.63, no. 10, pp. 407-410. EDN: GEIUPH. DOI: https://doi.org/ 10.1134/S1028335818100026.
6. Firsov S. V., Prokudin A. N., Burenin A. A. Creep and plastic flow in a rotating cylinder with a rigid inclusion, J. Appl. Industr. Math., 2019, vol. 13, no. 4, pp. 642-652. EDN: QGZHJN. DOI: https://doi.org/10.1134/S1990478919040070.
7. Rabotnov Yu. N. Creep Problems in Structural Members. Amsterdam, London, North-Holland Publ. Co., 1969, xiv+822 pp.
8. Gamer U., Lance R. H. Stress distribution in a rotating elastic-plastic tube, Acta Mech., 1983, vol.50, no. 1-2, pp. 1-8. DOI: https://doi.org/10.1007/BF01170437.
9. Gamer U., Mack W., Varga I. Rotating elastic-plastic solid shaft with fixed ends, Int. J. Eng. Sci., 1997, vol.35, no. 3, pp. 253-267. DOI: https://doi.org/10.1016/S0020-7225C96) 00085-7.
10. Antoni N. Contact separation and failure analysis of a rotating thermo-elastoplastic shrink-fit assembly, Appl. Math. Mod., 2013, vol.37, no. 4, pp. 2352-2363. DOI: https://doi.org/ 10.1016/j.apm.2012.05.018.
11. Mack W., Plochl M. Transient heating of a rotating elastic-plastic shrink fit, Int. J. Eng. Sci., 2000, vol.38, no. 8, pp. 921-938. DOI: https://doi.org/10.1016/S0020-7225C99) 00064-6.
12. Prokudin A. N. Exact elastoplastic analysis of a rotating cylinder with a rigid inclusion under mechanical loading and unloading, ZAMM, 2020, vol.100, no. 3, e201900213. EDN: SYSSIM. DOI: https://doi.org/10.1002/zamm.201900213.
13. Prokudin A. N., Burenin A. A. Elastoplastic analysis of a rotating solid shaft made of linear hardening material, Mech. Solids, 2021, vol.56, no. 7, pp. 1243-1258. EDN: UDGKGQ. DOI: https://doi.org/https://doi.org/10.3103/S0025654421070207.
14. Begun A. S., Burenin A. A., Kovtanyuk L. V., Prokudin A. N. Irreversible deformation of a rotating disk having angular acceleration, Acta Mech., 2021, vol. 232, no. 5, pp. 1917-1931. EDN: UPPSFS. DOI: https://doi.org/10.1007/S00707-021-02942-5.
15. Prokudin A. N., Firsov S. V. Elastoplastic deformation of a rotating hollow cylinder with a rigid casing, PNRPU Mechanics Bulletin, 2019, no. 4, pp. 120-135 (In Russian). EDN: VUXFSF. DOI: https://doi.org/10.15593/perm.mech/2019.4.12.
16. Prokudin A. N. Elastic-plastic analysis of rotating solid shaft by maximum reduced stress yield criterion, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2020, vol. 24, no. 1, pp. 74-94 (In Russian). EDN: LJTYOU. DOI: https://doi.org/10.14498/vsgtu1737.
17. Burenin A. A., Tkacheva A. V., Shcherbatyuk G. A. Calculation of the unsteady thermal stresses in elastoplastic solids, J. Appl. Mech. Tech. Phys., 2018, vol. 59, no. 7, pp. 1197-1210. EDN: LAVRJL. DOI: https://doi.org/10.1134/S0021894418070040.
18. Burenin A. A., Tkacheva A. V. Evolution of temperature stresses in the Gadolin problem of assembling a two-layer elastoplastic pipe, PNRPU Mechanics Bulletin, 2020, no. 3, pp. 20-31 (In Russian). EDN: HMUHTT. DOI: https://doi.org/10.15593/perm.mech/2020.3.03.
19. Burenin A. A., Tkacheva A. V. Piecewise linear plastic potentials as a tool for calculating plane transient temperature stresses, Mech. Solids, 2020, vol. 55, no. 6, pp. 791-799. EDN: EMMOEA. DOI: https://doi.org/10.3103/S0025654420060059.
20. Bykovtsev G. I., Ivlev D. D. Teoriia plastichnosti [Theory of Plasticity]. Vladivostok, Dal'nauka, 1998, 528 pp. (In Russian)
21. Mosolov P. P., Myasnikov V. P. Mekhanika zhestkoplasticheskikh sred [Mechanics of Rigidly Plastic Media]. Moscow, Nauka, 1988, 208 pp. (In Russian)
22. Kovtanyuk L. V., Shitikov A. V. On the theory of finite elastoplastic deformations of materials taking into account temperature and reological effects, Vestn. FEB RAS, 2006, no. 4, pp. 87-93 (In Russian). EDN: HZMRVR.
23. Norton F. H. The Creep of Steel at High Temperatures, Classic Reprint Series. London, Forgotten Books, 2017, 102 pp.
24. Firsov S. V. Irreversible Deformations of a Rotating Cylinder, Izv. Altai State Univ., 2018, vol.102, no.4, pp. 114-117 (In Russian). EDN: YABMRN. DOI: https://doi.org/10.14258/ izvasu(2018)4-21.
25. Prokudin A. N., Firsov S. V. Viscoplastic flow in a rotating hollow cylinder, Dal'nevost. Mat. Zh., 2018, vol. 18, no. 2, pp. 242-260 (In Russian). EDN: YUNRJB..
26. Burenin A. A., Tkacheva A. V. Assembly of a two-layered metal pipe by using shrink fit, Mech. Solids, 2019, vol.54, no. 4, pp. 559-569. EDN: BSSTLA. DOI: https://doi.org/ 10.3103/S0025654419040095.
27. Firsov S. V. Plastic and creep deformations of thick-walled cylinder with a rigid casing under internal pressure, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol.25, no. 4, pp. 696-715 (In Russian). EDN: TMRVNH. DOI: https://doi.org/10.14498/vsgtu1877.
28. Prokudin A. N., Firsov S. V. Creep analysis of rotating cylinder with free ends, Vestn. Chuvash. Gos. Ped. Univ. im I. Ya. Yakovleva. Ser. Mekh. Pred. Sost., 2018, no. 1 (35), pp. 63-73 (In Russian). EDN: UVVBZG.