Прокудин А.Н., Фирсов С.В. Упругопластическое деформирование вращающегося полого цилиндра с жестким внешним покрытием // Вестник Пермского национального исследовательского политехнического университета. Механика. 2019. № 4. С. 120-135. DOI: 10.15593/peim.mech/2019.4.12
Prokudin A.N., Firsov S.V. Elastoplastic deformation of a rotating hollow cylinder with a rigid casing. PNRPUMechanics Bulletin, 2019, no. 4, pp. 120-135. DOI: 10.15593/perm.mech/2019.4.12
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 4, 2019 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
Б01: 10.15593/регш.шееЬ/2019.4.12 УДК 539.374
УПРУГОПЛАСТИЧЕСКОЕ ДЕФОРМИРОВАНИЕ ВРАЩАЮЩЕГОСЯ ПОЛОГО ЦИЛИНДРА С ЖЕСТКИМ ВНЕШНИМ ПОКРЫТИЕМ
А.Н. Прокудин, С.В. Фирсов
Институт машиноведения и металлургии ДВО РАН, Комсомольск-на-Амуре, Россия
АННОТАЦИЯ
Рассматривается вращающийся полый цилиндр с закрепленными торцами, внутренняя поверхность которого свободна от напряжений, а внешняя зафиксирована от радиальных перемещений. Предполагается, что цилиндр изготовлен из идеального изотропного упругопластического материала, а деформации в нем являются малыми и представляют собой сумму упругих и пластических деформаций. Напряжения связаны с упругими деформациями законом Гука. Пластические деформации определяются с помощью условия Треска-Сен-Венана и ассоциированного с ним закона пластического течения. Скорость вращения цилиндра сначала монотонно возрастает до максимального значения, а затем снижается до нуля.
С помощью упругого решения найдена зависимость для критической скорости вращения, при которой начинается пластическое течение. Установлено, что в зависимости от толщины цилиндра и коэффициента Пуассона пластическое течение может начаться как на внутренней, так и на внешней поверхности цилиндра. При этом на стадии нагрузки в цилиндре появляются три пластические области, а на стадии разгрузки - четыре. Эти области соответствуют двум граням и двум ребрам призмы Треска.
Для каждой пластической области найдено точное аналитическое решение определяющей системы уравнений. Приведены системы условий на границах между областями, обеспечивающие непрерывность полученных решений во всем цилиндре. Рассмотрены два случая с пластическим течением, впервые начинающимся на внутренней поверхности цилиндра и на внешней поверхности цилиндра. Получены аналитические выражения для скоростей вращения, при которых появляются новые области. Установлена зависимость между скоростями зарождения повторного и первичного пластического течения. Также найдено значение максимальной скорости вращения, достаточное для полного перехода цилиндра в состояние повторного пластического течения. Выявлено, что введение жесткого внешнего покрытия позволяет значительно увеличить ресурс эксплуатируемой детали.
©ПНИПУ
О СТАТЬЕ
Получена: 12 сентября 2019 г. Принята: 07 декабря 2019 г. Опубликована: 30 декабря 2019 г.
Ключевые слова:
вращающийся цилиндр, полый цилиндр, плоская деформация, упругопластичность, пластическое течение, закон Гука, условие Треска-Сен-Венана, повторное пластическое течение, остаточные напряжения.
© Прокудин Александр Николаевич - к.т.н., в.н.с., e-mail: [email protected], : 0000-0002-5156-424Х. Фирсов Сергей Викторович - м.н.с., e-mail: [email protected], : 0000-0001-7446-6231.
Prokudin A. Nikolaevich - CSc in Technical Sciences, Leading Researcher, e-mail: [email protected], : 0000-0002-5156-424X.
Firsov S. Viktorovich - Junior Researcher, e-mail: [email protected]), : 0000-0001-7446-6231.
В , .-Ч ' Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International (д/ Qy License (CC BY-NC 4.0) _ -W!H This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
ELASTOPLASTIC DEFORMATION OF A ROTATING HOLLOW CYLINDER WITH A RIGID CASING
A.N. Prokudin, S.V. Firsov
Institute of machinery and metallurgy FEB RAS, Komsomolsk-na-Amure, Russian Federation
ABSTRACT
A rotating hollow cylinder with fixed ends is considered, the inner surface of which is free of stresses, and the outer one is fixed from radial movements. It is assumed that the cylinder is made of an ideal isotropic elastoplastic material, and the deformations in it are small and represent the sum of elastic and plastic deformations. Stresses are associated with elastic deformations by Hooke's law. Plastic deformations are determined using the Tresca - Saint-Venant condition and the plastic flow rule associated with it. The cylinder rotation speed first monotonical-ly increases to a maximum value, and then decreases to zero.
By using the elastic solution, the dependence is found for the critical rotation speed at which the plastic flow begins. It is established that, depending on the thickness of the cylinder and the Poisson's ratio, plastic flow can begin, either on the inner or on the outer surface of the cylinder. In addition, 3 plastic regions appear in the cylinder at the loading stage, and 4 plastic regions appear at the unloading stage. These regions correspond to two faces and two edges of the Treska prism.
For each plastic region, an exact analytical solution of the determining system of equations is found. The system of conditions at the boundaries between the regions providing continuity of the obtained solutions throughout the cylinder is given. Two cases are considered, i.e. the case with a plastic flow which starts first on the inner, and then on the outer surface of the cylinder. Analytical expressions are obtained for rotational speeds at which new regions appear. The relationship between the nucleation rates of the secondary and primary plastic flow is established. The value of the maximum rotation speed sufficient for a complete transition of the cylinder to the state of the secondary plastic flow was also found. It has been revealed that the adding of a rigid casing can significantly increase the resource of an exploited part.
©PNRPU
ARTICLE INFO
Received: 12 September 2019 Accepted: 07 December 2019 Published: 30 December 2019
Keywords:
rotating cylinder, hollow cylinder, plain strain, elastoplasticity, plastic flow, Hooke's law, Tresca - Saint-Venant condition, secondary plastic flow, residual stresses.
Введение
Расчет напряжений во вращающихся дисках и валах относится к классическим задачам механики деформируемых твердых тел. Опубликовано множество научных трудов, посвященных решению данных задач в различных постановках. Обычно рассматривается три варианта кинематики задачи. Это цилиндры с закрепленными концами, свободными концами и диски различных профилей. Решение подобных задач входит в различные монографии и учебники [1-9]. Так, С.П. Тимошенко впервые [4] предложил рассчитывать распределения напряжений в форме решения соответствующей краевой задачи теории упругости. Для упругих тел данную задачу можно считать принципиально решенной после того, как А.И. Лурье получил общее решение этой задачи [5] с произвольной зависимостью упругого потенциала от инвариантов тензора конечных деформаций. В работах [10-13] рассмотрены задачи упругого деформирования вращающихся цилиндров из функционально-градиентных материалов. В работе [14] приведено решение задачи вязкоупругого деформирования сплошных и полых цилиндров из функционально -градиент-ных материалов. Вращение гиперупругих цилиндров приводится в работах [15-17].
Необратимое деформирование материала вращающегося цилиндра в рамках жесткопластического анализа рассмотрел впервые А. Надаи [6]. Изучению дефор-
мирования вращающихся дисков в условиях ползучести посвящена работа Ю.Н. Работнова [2]. Также можно отметить работы [18-21], посвященные изучению деформации ползучести вращающихся цилиндров. Упру-гопластическое деформирование впервые было рассмотрено в [22]. Но построенное решение, как позже было показано в [23], имело неисправимый изъян, заключающийся в том, что не обеспечивалась непрерывность перемещений на упругопластической границе. В отечественной учебно-методической литературе последнее обстоятельство называют «ошибкой К. Шевченко» [24]. Проявляется такая ошибка в попытках получить решение упругопластических задач в скоростях без должного отслеживания выполнения требования непрерывности перемещений на упругопластических границах. В работах [23, 25] приводится корректное решение упругопластической задачи для вращающегося цилиндра с закрепленными концами. Остаточные напряжения и деформации в сплошном цилиндре после его остановки получены в работе [26]. Полый цилиндр с закрепленными концами исследовался в [27].
Сплошной и полый вращающиеся цилиндры со свободными концами рассматривались в работах [28] и [29] соответственно. В работе [30] приводится сравнение результатов численных расчетов вращающихся сплошного и полого дисков для условий пластичности Треска и Мизеса. Полуаналитическое решение для сплошного диска при использовании условия Мизеса
приведено в [31]. В работе [32] рассматривался вращающийся сплошной цилиндр из линейно-упрочняющегося упругопластического материала как с закрепленными, так и со свободными концами. Вращающийся цилиндр с нелинейным законом упрочнения рассматривался в работе [33]. Упругопластический анализ сплошных и полых цилиндров из функционально-градиентных материалов рассмотрен в работах [34-37]. Вязко-пластическое течение вращающегося полого цилиндра со свободными и закрепленными торцами рассматривалось в работе [38]. В последнее время приобрела актуальность задача нахождения остаточных напряжений в трубе после вращения. Решение данной задачи используется для описания процесса автофретирования труб за счет вращения [39-43].
Из обзора материалов, посвященных вращающимся дискам и цилиндрам, можно видеть, что большинство работ посвящены рассмотрению сплошных и полых тел. Случаев, когда в качестве граничных условий выступают перемещения, гораздо меньше. Среди таких работ можно отметить [44-46], где рассматривался вращающийся диск с жестким включением, и [47] - с полностью зафиксированным диском.
Целью данной работы является получение точного решения задачи об упругопластическом деформировании вращающегося полого цилиндра с зафиксированной радиальной компонентой перемещений на внешней границе. Полученное в итоге решение позволит вычислить напряженно-деформированное состояние (НДС) в таком цилиндре в любой момент времени. Также оно позволит выяснить, какое влияние на НДС окажет введение жесткого ограничения на радиальную компоненту вектора перемещений на внешней граничной поверхности.
1. Постановка задачи
Рассмотрим бесконечно длинную цилиндрическую среду с упругими и пластическими свойствами. Цилиндр ограничен двумя поверхностями г = а и г = Ь (Ь > а). Граничная поверхность г = Ь полагается жестко зафиксированной, г = а - свободной. Положим, что такой цилиндр вращается с угловой скоростью ю (/). В начальный момент времени ю (0) = 0. В дальнейшем скорость вращения постепенно увеличивается до значения ю = ютах. После достижения этого значения скорость начинает плавно снижаться вплоть до полной остановки вращения. В цилиндрической системе координат г, 6, г для изотропного материала получим одномерный процесс деформирования относительно координаты г . Для общности приведем ее к безразмерному виду:
Р =5 =a; 0 <5<р<1. b b
Деформации в материале будем считать малыми и состоящими из упругих er, eee, ez и пластических
Prr, Pee, Pzz . Ненулевыми компонентами деформаций в нашем случае будут только компоненты drr и dee тензора полных деформаций:
du u
drr = err + prr = T7T; dee = eee + pee = 7T;
SP P
dzz = ezz + Pzz = °
(1)
где и = иг / Ь - безразмерные перемещения.
Напряжения в деформируемом слое обусловливаются упругими деформациями в соответствии с законом Гука:
E 1
ст„„ =
ст y ( 1 + v )(1 - ■2v )
E 1
СТ y ( 1 + v )(1 - -2v )
E 1
(1 + v )(1 - 2v )
((1 -v) er + veee + vezz ); (verr +(1 -v) eee+vezz ); (2)
(verr + veee + (1 - v) ezz )>
где Е - модуль Юнга; V - коэффициент Пуассона; сту - предел текучести.
Соотношения, обратные к (2), имеют вид
е = СТ- (ст - уст™ - уст ); е„„ = СТ- - уст - уст );
гг ц, \ гг 66 гг / > 66 ^ \ 66 гг гг ; '
E
E
e =-£- (ст -
zz т-т V zz
E
vct -vctc
О.
(3)
Будем считать угловое ускорение достаточно малым, чтобы можно было пренебречь силой инерции, определяемой вращательным ускорением. Это позволяет оставаться в рамках одномерности. В таком случае единственное уравнение равновесия примет вид
dCTrr I CTrr-CTee. = _Qp; 0 = ^ > 0,
dp
в
(4)
где р - плотность материала.
Для учета пластических свойств деформируемого слоя воспользуемся теорией пластического течения. Для задания поверхности нагружения возьмем условие пластичности Треска-Сен-Венана, в нашем случае принимающее вид
f = max(|стгг -
-CTzzl )-1 = 0. (5)
стгг - ст66 I , |стгг - стгг I , |ст66 - стгг
Ассоциированный закон пластического течения запишется в виде
Фгг ; ЛР6>6>=^-Г—; Фгг =Х~Т— , (6)
ёстгг dст66 d стгг
где X - положительный коэффициент.
CTrr =
ст
У
Использование закона (6) вместе с условием (5) приводит к пластической несжимаемости. Из чего следует, что объемная деформация является чисто упругой. Исходя из этого и соотношений (1), (3) получим
ди + и = 5: (1 -2у)(ст„. + стее + ст ) . (7)
2. Построение математической модели
Граничные условия для полого цилиндра с внешней зафиксированной границей запишутся в виде
ст„ (5) = 0, и (1) = 0.
Разобьем процесс деформирования на стадии. При росте скорости - упругое деформирование, упругопластиче-ское и полная пластичность. При снижении скорости - упругое деформирование, появление повторного пластического течения. Рассмотрим более подробно данный процесс.
2.1. Упругое деформирование
Из соотношений (1), (2), (4), с учетом равенства нулю пластических деформаций, получим дифференциальное уравнение, решение которого примет вид
« = 1АР+ D - i QP3.
2 Р 8 E 1 - v
(8)
где Д, D2 - константы интегрирования. Напряжения примут значения
_ E D,p2 -2D2 (1 -2v) 13-2v
оу 2Р2 (1 + v)(1 -2v) 8 1 -v
E gp2 + 2D2 (1 - 2v) - 11 + 2v oy 2P2 (1 + v)(1 -2v) 8 1 -v
p2Q,
P2Q;
О = v(0rr +°00)-
(9)
Для текущих граничных условий константы интегрирования Д, Б2 рассчитываются следующим образом:
D =
(1 + v)(1 - 2v) 1 - 2v + ( 3 - 2v)54
1 E 4 (1 -v) S2 +1 - 2v
oy (1 + v)(1 -2v) 1 -(3-2v)S2 E 8 (1 -v) S2 +1 - 2v
Q;
d2 =
S2Q.. (10)
Полученное упругое решение (8)-(10) справедливо для скоростей О < . Выход на поверхность нагруже-ния может происходить как на внутренней граничной поверхности Р = 5 , так и на внешней р = 1. Найдем значение скоростей, при которых это происходит. На внутренней границе напряжения удовлетворяют неравенству стее > стгг > ст^, на внешней - стее = стгг > стгг. Условие пластичности (5) на внутренней поверхности примет вид
а на внешней
оее-огг =1,
оаи -о„ = 1; о„ - о = 1.
Используя полученные условия пластичности и упругое решение (9), (10), получим следующие скорости начала пластического течения:
= 4
Q{ = 4
1 -v
1 -v S2 +1 -2v
1 - 2v
(1 -8- )'
S2 +1 -2v
1 -
2v (1 -S2)(1 -2v +S2 (3-2v))
(11)
Сравнив их, придем к выводу, что при 5 = / (2 -V) пластическое течение начнется одновременно на двух граничных поверхностях (рис. 1). При больших значениях 5 - на внешней граничной поверхности. При меньших - на внутренней. Отметим, что при 5 > ¡/>/3 пластическое течение всегда будет начинаться на внешней границе.
а б
Рис. 1. Условие начала пластического течения: а - зависимость между 5 и v , при которых пластическое течение начнется на внутренней границе (ниже графика), на внешней - (выше графика); б - зависимость скорости начала пластического течения на внутренней и внешней граничных поверхностях от значения 5 при v = 0,3
Fig. 1. The condition for the start of plastic flow: a - the relationship between 5 and v , under which the plastic flow begins on the inner boundary (below the graph), on the outer (above the graph); b - the dependence of the velocity of the start of plastic flow on the inner and outer boundary surfaces on the value 5 at v = 0.3
2.2. Упругопластическое деформирование
При скорости вращения Q = Q1 в цилиндре выполняется условие пластичности Треска (5) на внутренней граничной поверхности р = 5. Вследствие этого появляется область пластического течения, которую обозначим как I. Напряжения в области I (5<Р<Р1) удовлетворяют неравенству стее > стгг > arr. В оставшейся части цилиндра (Р1 < р <1) продолжается процесс упругого деформирования.
оее =
ст
y
С увеличением скорости вращения до О = О2 на внешней граничной поверхности Р = 1 выполняется условие пластичности (5) и зарождаются две области пластического течения, которые обозначим как II и III. Область II (Р2 < Р < Р3) соответствует грани призмы Треска стее > стгг > ст.., а область III (р3 < Р < 1) - ребру стее = стгг >стгг. В цилиндре по-прежнему сохраняется область упругого деформирования (< р < Р2).
При достижении скорости вращения О = О3 упру-гопластические границы Р1 и Р2 совмещаются между собой и упругая область исчезает. Так как области I и II принадлежат к одной грани призмы Треска, они сливаются вместе. Получившуюся область обозначим как II. При дальнейшем увеличении скорости вращения, вплоть до О = Ошах, граница Р3 будет постепенно продвигаться в направлении внутренней граничной поверхности р = 5 . Рассмотрим каждую область.
Упругая область. Полученное ранее упругое решение (8), (9) остается справедливым, однако значение констант интегрирования Д, Вг необходимо определять из условий на новых границах упругой области. При О1 < О < О2 это границы р = Р1 и р = 1. Затем при О2 < О < О3 это границы в = и Р = Р2.
Область I
Напряжения соответствуют грани призмы Треска стее > стгг > ст... Условие пластичности (5) принимает вид
стее-стгг =!.
где С2 - константа интегрирования.
Пластические деформации можно найти из (1), с учетом (3), (12), (13), (14):
^ =- ^ + Е (!-V2)-4 Е (! + у)(!-2у)ОР2;
Рее=-Р..; Рг = 0. (15)
Область II
Область II соответствует той же грани призмы Треска, что и область I. Соответственно, решения (12)-(15) справедливы и для области II с той поправкой, что константы интегрирования С1, С2 заменятся на С3, С4.
Область III
Напряжения в этой области соответствуют ребру призмы Треска стее = стгг > ст... Условие пластичности (5) примет вид
стш - ст =1; ст -ст =1.
ее гг > гг гг
С помощью данного условия и уравнения (4) получим значения напряжений
1 ,
ст.. = С5--Р2О + 1пР; стее=агг =1 + ст„, (16)
где С5 - константа интегрирования.
Подставив полученные напряжения (16) в (7), получим дифференциальное уравнение для перемещений
|+Р=Е(1 -2" )(2+3С5- 3 Р2О+31п Р),
С помощью данного условия и уравнения (4), получим решение которого имеет вид
1 2
drr =Ci - 2 ß2^ + ln ß; aee=1 + arr, (12)
где C - константа интегрирования.
Используя ассоциированный закон пластического течения (6), найдем
Prr = - Pee; Pzz = 0. Из (1) следует, что ezz = 0 . В результате осевое напряжение azz определяется так же, как и в случае упругости:
a zz =v(arr +a ee) . (13)
Подставив полученные напряжения (12), (13) в (7), получим дифференциальное уравнение для перемещений
dß + ß = E(1 + V)(1 -2v)(1 + 2Ci-ß2^ + 2lnß) ,
3 a C
u = (1 -2v) C5ß + Сб -2 EK ' ß
решение которого имеет вид ay
С
u =E (1 + v )(1 -2v ) C1ß^-ßL -
-T^y (1 + v )(1 - 2v) Qß3 + (1 + v )(1 - 2v )ß ln ß, (14)
4 E E
- - E (1 - 2v)Qß3+4 E (1 - 2v)(1+6ln ß)ß' (17)
где C6 - константа интегрирования.
Пластические деформации находятся из (1), (3), (16), (17):
prr=2 Oy (1 - 2v) с- С"+
+ 4 E (7 - 6V + 2 (1 - 2v) ln ß)-8 °Ey(1 - 2v)Qß2;
Pee = 2 E (1 - 2V) C5 + F -
- 1 E ( 3-2v-2 (1 - 2v) ln ß)+1 E (1 - 2v)Qß2;
Pzz =-Prr - Pee =- 'E(1 - 2v) C5 -
(1 -v + (1 - 2v) ln ß)+1' (1 - 2v)Qß2. (18)
E 2 E
2.3. Разгрузка
После достижения значения Q = ^max скорость вращения начинает уменьшаться. Это приводит к уменьшению напряжений и, соответственно, упругому деформированию цилиндрического слоя. На процесс упругости оказывают влияние накопленные пластически деформации. Обозначим их как
Prr = Prr i^-max ) ' рвв = Рвв max ) ' pzz = Pzz max ) •
Преобразуя уравнение равновесия (4) с учетом (1), получим уравнение равновесия в перемещениях для стадии разгрузки
_ (1 + v)(i-2v)
- + -——- u = —^-^-Qß + p;
du 1 du 1
ст„
dß2 ßdß ß2 E 1 - v
P =
1 1-2v
dprr , V I dp66 , dp.
1-V (prr pp66)+ dß + 1_v[ dß + dß J-
ß 1-,
dß dß
В общем случае решение данного неоднородного уравнения примет вид
В 1
и = 8^+^ + 2 в| Р d р-
- -1 fß2 p dß-^- (1 + V)(1 -2v) ß3 2ßJ E 1 -v
Q,
где В1, В2 - константы интегрирования.
Рассмотрим области с накопленными пластическими деформациями. Область I
и = В.Р+ 8-1Е )«
X р3 -1 (1 + у)(1 -2у)(1 -21пр)р.
2 Е
Область II
и = В3Р + В--1 Е + У)^ + 2 (1 -2у)^тах )х
х в3 - 4 (1 - 2у)(1 - 21п в)в.
2.4. Повторное пластическое течение
При дальнейшем снижении скорости вращения в цилиндре появляются области повторного пластического течения. В таких областях пластические деформации будут состоять из накопленных (ргг, р66, ргг)
и вторичных (згг, 566,) пластических деформаций:
p = p + s ;
г rr с rr rr '
1+ s66 ; pzz = pzz + Szz -
После снижения скорости вращения до значения О = О4 в цилиндре на внутренней граничной поверхности р = 8 зарождается область повторного пластиче-
ского течения. Обозначим эту область III 5 < ß < ß4. Ей соответствует грань призмы Треска arr > стгг > стее. В оставшейся части цилиндрического слоя (ß4 < ß < ß3 и ß3 < ß < 1) продолжается упругое деформирование.
При уменьшении скорости вращения до Q = Q5 зарождаются две пластические области на внешней граничной поверхности ß = 1. Обозначим их IV и V. Области IV (ß5 < ß < ß6) соответствует грань призмы Треска стгг >стгг >стее, а области V (ß6 < ß < 1) - ребро призмы Треска стгг > стее = стгг. Упругая область, соответственно, будет занимать оставшуюся часть деформируемого слоя (ß4 < ß < ß3 и ß3 < ß < ß5).
При дальнейшем уменьшении скорости, вплоть до Q = 0 , границы областей повторного пластического течения постепенно продвигаются внутрь деформируемого слоя в направлении поверхности ß = ß3. Рассмотрим подробней получившиеся области повторного пластического течения.
Область III
Напряжения соответствуют грани призмы Треска стгг > стгг > стее. Условие пластичности (5) принимает вид
CTrr -СТ66 =1-
Из данного условия и уравнения (4) получим 1 2
ц = A -2ß2Q-lnß; ^66=CTrr-1,
(19)
(20)
где А1 - константа интегрирования.
Используя ассоциированный закон пластического течения (6), найдем
5гг =-566; ^г = (21)
Учитывая накопленные в этой области пластические деформации (15), получим ргг = 0. Из (1) следует егг = 0 . В результате осевое напряжение стгг определяется так же, как и в случае первичного пластического течения (13). Подставив напряжения (13), (20) в (7), получим дифференциальное уравнение
дР + Р = Е (1 + у)(1 -2у )(2 А,-Р2О-21п р-1),
решение которого имеет вид
СТ у,
A
u =E (1 + v )(1 -2v ) A1ß^-ßL -
- 7ПТ(1 + v)(1 -2v)Qß3 (1 + v)(1 -2v)ßlnß, (22)
4 E E
где A2 - константа интегрирования.
Используя (1), (3), (13), (15), (20), (21), (22), получим
prr =- A (1 -v2)-1 (1 + v)(1 -2v)Qß2;
ß2 E v >4 E
p66=-prr; pzz = 0-
(23)
Область IV
Область IV соответствует той же грани призмы Треска, что и область III. Соответственно, для нее справедливо то же условие пластичности (19). Так же справедливо (20) с заменой A1 на A3. Однако накопленные пластические деформации принимают другие значения (18). Что в свою очередь, с учетом (21), приводит к следующей зависимости:
3 a A
u = (1 -2v) A5ß+A -2 ' ß
-3Е-2у)°Р3 -4Е-2у)(:+б1пp)p, (29)
где А5, А6 - константы интегрирования.
Из уравнений (1), (3), (28), (29) можно найти пластические деформации:
' zz =V(arr + ' ee)--Pzz .
CT,,
(24)
Подставив (20), (24) в (7), получим дифференциальное уравнение в перемещениях:
du u ' ,, 4 /, „ 4 — + - = —(1 + v)(1 - 2v) x 5ß ß E A '
x( 2 A3-Qß2-2ln ß-1)-(1 -2v ) P zz,
Prr =1 ay (1 - 2v) A5 - A -
2 E
ß2
-4E (7 - 6V + 2 (1 - 2v) ln ß)-85 EU - 2v)Qß2;
Pee = 2 E (1 - 2V) A5 + A + +4 E (3 - 2v- 2 (1 - 2V) ln ß)+8 E (1 - 2V)q^
решение которого при накопленных пластических деформациях (18) имеет вид
и = ^ (1 -2у )(1 + v ) ^3Р + А +1 ^ (!-2^ )х
E
ß 4 E
1 'h
х(1 - 2 (1 + 4у) 1п Р)Р~ -У (1 - 2у )х
8 Е
х( 2 (1 + V )О + (1 - 2у ) О шах )Р3 + 2 Е (! - 2у )2 С5Р, (25)
где А4 - константа интегрирования.
Из (1), (3), (20), (24), (25) можно найти пластические деформации:
Р.. =-А-4Е(! + 6у-8у2 -2(!-2^5 - 1пР))-
-1 у -2У)(2 + У)О + (3-2У)Ошах )Р2;
Рее = - Ртт ; Ргг = Ргг =- ^ ^Р +С5 )(1 - 2v)) +
+ 2 E (1 - 2v)ß2Qmax.
(26)
Область V
Напряжения соответствуют ребру призмы Треска ст . > ст ее = стгг. Условие пластичности (5) принимает вид
a „ - a „„ = 1; a „ - a „ = 1.
rr zz
(27)
По аналогии с областью III первичного пластического течения получим следующие значения напряжений:
1 2
ст.. = А--Р2О - 1пР; ст ее = ст гг = ст.. -1 (28)
2
и перемещений
Pzz = -Prr - Pee = -E(1 - 2v)A5 +
+ (1 - v + (1 - 2v) ln ß) +1 °U (1 - 2v)Qß2. (30) E 2 E
3. Случай начала пластического течения на внутренней границе
Рассмотрим следующие физические и геометрические параметры цилиндра:
5 = 0,4; v= 0,3; ay = 250 МПа;
E = 21 ГПа; Qmax = 13.
Для текущих параметров 5 <^v/ (2 -v) . Следовательно, пластическое течение впервые начнется на внутренней граничной поверхности (см. рис. 1). Из (11) получим Q1 = 50/9 .
На внутренней граничной поверхности зарождается пластическая область I. Для нахождения значений постоянных интегрирования для двух областей решим следующую систему уравнений:
Or (5) = 0; ue (1) = 0; u' (ß1 ) = ud (ß1); 'гг (ß1 ) = t (ß1),
где u',a Гг - перемещения и напряжения для области I; uel , rerl - перемещения и напряжения для упругой области. Далее верхний индекс будет использоваться для обозначения области.
Из данной системы уравнений можно найти значения констант интегрирования D1, D2, C1, С2. Для нахождения значений ß1 и Q2, при которых появятся пластические области II и III, решим добавочные уравнения
ст66(в1 )-стег (Р1 )=1; ст66(1 )-стег (1 )=1.
Решая второе уравнение при найденных константах интегрирования, получим
= 4
1 -v
1 -
(1 -2v)(l + 2Р2 ln| ] + P? 2^(1 -2v)(l + P4) + 2p? (1 -S2 (1 -v))
• (31)
Подставив полученную скорость в первое уравнение, получим
чения констант и скорости (35) и приравняв его нулю, получим биквадратное уравнение. Его решение можно записать в виде р2 = a±\[b . Так как нас интересует момент, когда пластические деформации равны нулю только в одной точке, данное уравнение распишется на два: b = 0 и р2 = a . Они примут вид
(1 - 2v)(p4 (W-(1 - 2v))-p2 + + (1 - 2v) S2p2 +6(p2 logS-S2 log Рз)) +
(1 -2v)(1 + pf)lnР1 -v(1 -Pf ) + 2(1 -v)S2 = 0. (32) +3W-S2 (5-4v) = W (1 -v)(1 -v2)(p-S2)2;
В связи с нелинейностью данного уравнения по Р1 для его решения будем использовать численный метод Ньютона. В итоге искомые величины примут значения
p1 = 0,42; = 5,95.
(33)
Далее, при О > О2 в цилиндре присутствуют три области пластического течения и одна - упругого деформирования. К неизвестным добавятся четыре константы интегрирования С3, С4, С5, С6 и две границы Р2, Р3. Для определения всех неизвестных воспользуемся системой уравнений
П.
. (S) = 0; иш (1) = 0;
p = p1 : и = ue p = p2 : uel = и' p = p3 : и11 =u
CTrr = CT„. ; CTee-CTrr = 1;
_el _II _el _el л
О = О ; Gaa - С = 1;
rr rr' ee rr '
(34)
II III II III II III
_" a=a ; a=a •
rr rr ' zz zz
Для нахождения момента О3, когда произойдет схлопывание областей I и II, данных уравнений недостаточно, даже при введении дополнительного условия Р1 = Р2. Чтобы найти этот момент, сначала рассмотрим случай, когда О > О3. В этом случае в цилиндре остаются только две области пластического течения II и III и, соответственно, пять неизвестных С3, С4, С5, С6, Р3. Итоговая система уравнений примет вид
стг (5) = 0; и111 (1) = 0; и11 (Р3 ) = и111 (Р3); ^(Рэ(Рэ); стГг (Рэ) = оГ(Р3).
Из первых четырех уравнений найдем значения констант. Подставив найденные константы в пятое уравнение, можно найти зависимость между скоростью вращения и положением границы:
" = 2
1 -v+(1 -2v)( log p3-log S) (1 -2v)(p32-S2)
(35)
В момент О = О3 на поверхности Р = Р1 =Р2 необратимые деформации равны нулю. Подставив в уравнение для пластических деформаций (15) найденные зна-
p2 =p2 =W- (1 -v)(p2 -S2);
W = 1 -v+(1 -2v)( log p3-log S) • (36)
Из (35), (36) получим следующие значения: "3 = 7,6894; p1 = p2 = 0,6746; p3 = 0,9106 .
Подставив в (35) значение " = "max, получим следующее положение границы:
"max = 13; p3 = 0,721.
На (рис. 2) представлены распределение напряжений в указанные моменты, итоговые пластические деформации и движение границ областей пластического течения.
После достижения максимального значения скорость вращения начинает постепенно уменьшаться. При этом в цилиндрическом слое присутствуют две области упругого деформирования I и II, разделенные граничной поверхностью p3. Она не меняет своего положения. Найдем значение скорости вращения, при которой появятся области повторного пластического течения. Для этого разрешим условия пластичности на внутренней (19) и внешней (27) граничных поверхностях. В итоге получим
Q' = " -
4 max
1 -v S +1 -2v
1 - 2v
= "max - 8
(1 -S2 )2
1 - v S2 +1 -2v
1 -
2v (1 -S2 )(1 - 2v + S2 (3 - 2v))
. (37)
Сравнив полученные формулы с (14), получим О4 = Отах - 2О; и О4 = Отах - 2О?. Из этого можно сделать вывод, что повторное пластическое течение впервые начнется на той же граничной поверхности, на которой началось первичное пластическое течение. Из соотношений (37) также можно найти минимальное значение Отах, при котором появится область повторного пластического течения. В нашем случае это Отах > 100/9 и 11,1111. При Отах = 13 повторное пластическое течение начнется в момент О, = 1,8889 .
д е
Рис. 2. НДС цилиндра при росте скорости вращения Q : а - напряжения при Q = Q1; б - напряжения при Q = Q2; в - напряжения при Q = Q3 ; г - напряжения при Q = Qmax ; д - пластические деформации при Q = Qmax ; е - развитие областей пластического течения
Fig. 2. Stress-strain state of the cylinder with increasing angular speed Q : a - stresses at Q = Q1; b - stresses at Q = Q2; c - stresses at Q = Q3 ; d - stresses at Q = Qmax ; e - plastic strains at Q = Qmax ; f - evolution of primary plastic regions
При О < О 4 в цилиндрическом слое находятся одна область пластического течения и две области упругого деформирования. Следовательно, имеется семь неизвестных В1,В2,В3,В4,А1,А2,Р4. Составим следующую систему уравнений для их нахождения:
ст? (5) = 0; и1 (1) = 0;
Р = Р4 : и111 = и'; Р = Р3: и1 = и'';
СТ = СТ ;
rr rr >
ст! = СТ1'.
ст'гг -СТ ее = 1; (38)
Для нахождения скорости О5, при которой на внешней граничной поверхности Р = 1 зародятся две
области повторного пластического течения IV и V, введем дополнительное уравнение
ст.. (1)-ст ее(1)=1.
Из данного уравнения найдем значение О5:
Q =Qmax - 8
1 -v
(1 -2v)fl + 2p4ln ^ W
1 -
2v (1 -2v)(1 + P4 ) + 2p4 (1 -52 (1 -v))'
Сравнив его с (31), заметим, что О5 =Ошах - 2О2, если заменить Р4 на Р1. Подставив полученную скорость в условие пластичности в системе (38), получим
г
(1 -2у)(1 + Р4) 1п в4-у(1 -р^ ) + 2 (1 -у)52 = 0. (39)
Как видим, полученное уравнение полностью совпадает с (32), только в качестве неизвестной стоит Р4. Соответственно, при решении получаем то же положение упругопластической границы:
Р4 = 0,42; Q5 =Qmax -2Q2 = 1,1.
(40)
При О < О5 добавятся неизвестные А3, А4, А5, А6,Р5,Р6. Для их нахождения внесем следующие изменения в систему (38):
0,6 0.7
в
uV (1) = 0;
ая = aV ; rr rr > a" rr "aee =
^¡V „V . a = a ; rr rr ГГ1V aee = aIV. zz
Р = Р6: и17 = и7;
Из получившейся системы найдем положения границ при условии О = 0 :
Р4 = 0,4545; Р5 = 0,9329; Р6 = 0,9691. (42)
На (рис. 3) представлены распределение напряжений в указанные моменты, итоговые пластические деформации, движение границ областей пластического течения, а также распределения перемещений в цилиндре при разных скоростях вращения.
д е
Рис. 3. НДС цилиндра при снижении скорости вращения Q : а - напряжения при Q = Q4; б - напряжения при Q = Q5 ; в - напряжения при Q = 0; г - пластические деформации при Q = 0; д - развитие областей пластического течения; е - перемещения при разных Q
Fig. 3. Stress-strain state of the cylinder with decreasing angular speed: a - stresses at Q = Q4; b - stresses at Q = Q5; c - stresses at Q = 0 ; d - plastic strains at Q = 0; e - evolution of secondary plastic regions; f - displacement at different Q
а
4. Случай начала пластического течения на внешней границе
Рассмотрим по аналогии случай, когда пластическое течение начинается впервые на внешней граничной поверхности. Для этого возьмем следующие параметры:
5 = 0,5; v= 0,3; сту = 250 МПа; Е = 21 ГПа; О^ = 18.
При таких параметрах, согласно (11), пластическое течение впервые начнется при скорости
01 = 91/15 и 6,0667 . Система уравнений для получения значения констант и нахождения положения границ областей будет выглядеть схоже с системой (34), за исключением уравнений для области I и границы Р1.
По аналогии с (31) аналитически вывели зависимость
02 от значений Р2 и Р3 (в связи с ее громоздкостью не
будем приводить ее здесь). В итоге получили следующие параметры, при которых пластическое течение начинается на внутренней граничной поверхности:
О2 = 7,7602; Р2 = 0,8075; Р3 = 0,922 . (43)
По аналогии с предыдущим случаем получим
О3 = 8,3845; Р! = Р2 = 0,6461; Р3 = 0,8985 (44)
Q max = 18; ß3 = 0,6934.
При снижении скорости вращения момент начала повторного пластического течения определяется по формулам (37). В нашем случае это Q4 = 88/15 » 5,8667 . Для дальнейшего решения используются системы (38), (41), только без уравнений для области III и границы ß4. Также удается вывести аналитическую зависимость скорости вращения, при которой начинается повторное пластическое течение на внутренней границе. Сравнив полученную формулу с формулой для Q2 , получим Q5 = Qmax - 2Q2, если заменить ß4, ß5 на ß2, ß3. В итоге получим
Q5 =Qmax - 2Q2 = 2,4796; ß5 = 0,8075; ß6 = 0,922. (45)
Можно заметить, что, как и в предыдущем случае, положения границ областей в момент начала повторного пластического течения на другой стороне цилиндра (40), (45) совпадают с таковыми при первичном пластическом течении (33), (43).
При дальнейшем уменьшении скорости для описания среды будем использовать систему (38), (41). Однако при скорости Q6 произойдет переход области IV через границу ß3 и она разобьется на две подобласти: IV и VI. Область VI будет идентична области III повторного пластического течения, только с заменой ко-
эффициентов интегрирования с А1, А2 на А7, А8. Это произойдет при следующих значениях:
О6 = 1,3294; Р4 = 0,5992; Р 5 = Р3 = 0,6934; Р6 = 0,9004.
Для продолжения расчета в системы (38), (41) внесем изменения:
VI.
J VI.
ст ... = ст „ .
ß = ß5 : u = u
VI IV VI IV ß = ß3 : u =u ; CTrr = CT rr •
CT 'rr -CT 6 6 =1;
(46)
При дальнейшем уменьшении скорости произойдет слияние областей III и VI. Для нахождения этого момента воспользуемся теми же соображениями, что и ранее для случая первичного пластического течения. Системы уравнений (38), (41), (46) преобразуются к виду
CT? (5) = 0; uV (l) = 0;
ß = ß3 : uVI =uIV; ß = ß6 : uIV = uV;
(47)
IV V
ст„ = ст„.
IV IV CTaa = CT „ •
Найдя из первых шести уравнений значения констант и подставив их в последнее уравнение, выразим скорость через положения границ:
Q = Q - 4
max
1 -v + (l-2v)( log ß5-log 5) (1 -2v)(ß2-52)
Сравнив полученное выражение с (35), получим Q7 = Qmax - 2Q3, если заменить Р5 на Р3. Выполняя те
же преобразования, что и ранее - только вместо равенства нулю пластические деформации должны равняться остаточным пластическим деформациям p" = p'r, - получим те же уравнения (36), только с заменой рр Р2 на Р4, Р5 и Р3 на Р6. Заметим, что сама граница Р3 не входит в полученные уравнения. Из этих уравнений получим
Q7 =Qmax - 2Q3 = 1,231; р4 =р5 = 0,6461; Р6 = 0,8985.
При остановке вращения граница между оставшимися двумя областями принимает положение Р6 = 0,8734. На (рис. 4) представлены графики распределения напряжений и пластических деформаций, характеризующие ключевые стадии процесса.
На основе полученных результатов можно составить следующие зависимости, которые с учетом ранее введенного обозначения для скорости вращения (4) примут вид
Q =Qmax - 2Q ; = <x - 2®?; Pr = Pl , (47)
где Qr, ®r - скорости вращения на стадии разгрузки (уменьшения скорости вращения); Ql, ®l - скорости вращения на стадии нагрузки (увеличения скорости вращения); Pr, Pl - положения границ имеющихся пластических областей на стадиях разгрузки и нагрузки соответственно.
и
0,5 0.6 0,7 0,8 0,9 1,0 0,5 0,6 0,7 0,8 0,9 1,0
0 0,5 1,0 1,5 0 12 3 4 5 6
д е
Рис. 4. НДС цилиндра при начале течения с внешней границы: а - напряжения при Q = Qmax ; б - пластические деформации при Q = Qmax ; в - напряжения при Q = 0; г - пластические деформации при Q = 0; д - движение границ при увеличении Q ;
е - движение границ при уменьшении Q Fig.4. Stress-strain state of the cylinder with plastic flow first started at outer border: а - stresses at Q = Qmax ; b - plastic strains at Q = Qmax ; c - stresses at Q = 0; d - plastic strains at Q = 0; e - evolution of primary plastic regions; f - evolution of secondary
plastic regions
Полученное соотношение (47) говорит о том, что мы можем на основе данных, полученных на стадии увеличения скорости вращения, проанализировать поведение цилиндрической среды на стадии уменьшения скорости вращения. А именно: где и какие области пластического течения образуются. Это позволяет, зная скорость, при которой впервые появляется пластическое течение в цилиндре, найти такое максимальное значение скорости, при которой при разгрузке в теле не
возникнет повторного пластического течения. Либо в конструкциях, работающих на критических скоростях, позволяет определить границу, до которой можно опускать скорость вращения, чтобы не происходило излишнее пластическое деформирование металла.
Однако данная зависимость (47) может не выполняться для других граничных условий. Что можно увидеть, например, в работах [21, 22] для сплошного цилиндра.
Также следует отметить, что добавление граничных условий на перемещения сильно увеличивает значение скоростей вращения, при которых начинается пластическое течение. Так, в работе [38] для полого цилиндра со свободными концами скорости, при которых пластическое течение появится и распространится на всю деформируемую среду, составляют Q1 = 1,1859 и Q3 = 2,9384 соответственно. Из той же работы [38] для полого цилиндра с закрепленными концами получим Q1 = 1,2609 и Q3 = 3,8717 . По формулам (11), (36) для полого цилиндра с жестким внешним покрытием с закрепленными концами для тех же параметров среды получили Q1 = 3,342 и Q3 = 8,5139 . С практической точки зрения это означает, что добавление ограничения на перемещения значительно повысит ресурс эксплуатируемой детали. В данном случае добавление жесткого ограничения на внешнюю граничную поверхность позволяет увеличить значение критической скорости вращения ю1 чуть более чем в полтора раза.
Выводы
Найдено точное аналитическое решение определяющей системы уравнений для задачи вращения полого цилиндра с жестким внешним покрытием.
Библиографический список
1. Толоконников Л.А. Механика деформируемого твердого тела. - М.: Высшая школа, 1979. - 318 с.
2. Работнов Ю.Н. Ползучесть элементов конструкций. -М.: Наука, 1966. - 752 с.
3. Работнов Ю.Н. Механика деформируемого твердого тела. - М.: Наука, 1979. - 744 с.
4. Тимошенко С.П., Гудьер Дж. Теория упругости. - М.: Наука, 1979. - 576 с.
5. Лурье А.И. Нелинейная теория упругости. - М.: Наука, 1980. - 512 с.
6. Nadai A. Theory of Flow and Fracture of Solids, Volume One. 2nd Edition. - New York: McGraw Hill, 1950. - 572 p.
7. Chakrabarty J. Theory of Plasticity: Third Edition. - Elsevier, 2005. - 896 p. DOI: 10.1016/C2009-0-05774-2
8. Sadd M.H. Elasticity: theory, applications, and numerics. -Elsevier, 2014. - 582 p. DOI: 10.1016/C2012-0-06981-5
9. Timoshenko S. Theory of elasticity. - New York: McGraw Hill, 2010. - 567 p.
10. Zenkour A.M. Analytical solutions for rotating exponentially-graded annular disks with various boundary conditions // International Journal of Structural Stability and Dynamics. - 2005. -Vol. 05. - No. 04. - Р. 557-77. DOI: 10.1142/S0219455405001726
11. Eraslan A.N., Akis T. On the plane strain and plane stress solutions of functionally graded rotating solid shaft and solid disk problems // Acta Mechanica. - 2006. - Vol. 181. - No. 1. -P. 43-63.
12. Akis T., Eraslan A.N. The stress response and onset of yield of rotating FGM hollow shafts // Acta Mechanica. - 2006. -Vol. 187. - No. 1. - P. 169-187.
Найдена зависимость для критической скорости вращения, при которой начинается пластическое течение. Установлено, что в зависимости от толщины цилиндра и коэффициента Пуассона пластическое течение может начаться как на внутренней, так и на внешней поверхности цилиндра.
Получены аналитические выражения для скоростей вращения, при которых появляются новые пластические области. Установлена зависимость между скоростями зарождения повторного и первичного пластического течения. Найдено значение максимальной скорости вращения, достаточное для полного перехода цилиндра в состояние повторного пластического течения.
Установлено, что добавление жесткого ограничения на внешнюю поверхность позволяет в полтора раза повысить значение критической скорости, при которой впервые появится пластическое течение.
Благодарность
Работа выполнена в рамках государственного задания № 075-00414-19-00.
Acknowledgment
The work was carried out within the framework of the state task No. 075-00414-19-00
13. Sburlati R. Analytical elastic solutions for pressurized hollow cylinders with internal functionally graded coatings // Composite Structures. - 2012. - Vol. 94. - No. 12. - P. 3592-3600.
14. Zenkour A.M., Elsibai K.A., Mashat D.S. Elastic and viscoelastic solutions to rotating functionally graded hollow and solid cylinders // Applied Mathematics and Mechanics. - 2008. -Vol. 29. - No. 12. - P. 1601-1616. DOI: 10.1007/s10483-008-1208-x
15. Horgan C.O., Saccomandi G. Large deformations of a rotating solid cylinder for non-gaussian isotropic, incompressible hyperelastic materials // J. Appl. Mech. - 2000. - Vol. 68. - No. 1. -P. 115-117. DOI: 10.1115/1.1349418
16. Bagheri A., Taghizadeh D., Darijani H. On the behavior of rotating thick-walled cylinders made of hyperelastic materials // Meccanica. - 2016. - Vol. 51. - No. 3. - P. 673-692. DOI: 10.1007/s11012-015-0233-x
17. Chadwick P., Creasy C.F.M., Hart V.G. The deformation of rubber cylinders and tubes by rotation // The ANZIAM Journal. - 1977. - Vol. 20. - No. 1. - P. 62-96. DOI: 10.1017/S0334270000001454
18. Немировский Ю.В., Янковский А.П. Оценка предельной скорости вращения искривленных металлокомпозит-ных дисков в условиях ползучести // Вестн. Чуваш. гос. пед. ун-та им. И.Я. Яковлева. - 2010. - № 2 (66). - С. 102-111.
19. Прокудин А.Н., Фирсов С.В. Расчет ползучести вращающегося цилиндра со свободными концами // Вестн. Чуваш. гос. пед. ун-та им. И.Я. Яковлева. Сер.: Механика предельного состояния. - 2018. - № 1 (35). - С. 63-73.
20. Банщикова И.А. Напряженно-деформированное состояние и длительность до разрушения вращающихся дисков при ползучести // Вестник Пермского национального исследо-
вательского политехнического университета. Механика. -2018. - № 4. - С. 20-32. DOI: 10.15593/peim.mech/2018.4.02
21. Бегун А.С., Ковтанюк Л.В. Расчет напряжений, деформаций и перемещений во вращающемся диске в условиях ползучести // Вестн. Чуваш. гос. пед. ун-та им. И.Я. Яковлева. Сер.: Механика предельного состояния. - 2019. - № 1 (39). -С. 84-93. DOI: 10.26293/chgpu.2019.39.1.011
22. Hodge P.G., Balaban M. Elastic-plastic analysis of a rotating cylinder // International Journal of Mechanical Sciences. - 1962. -Vol. 4. - No. 6. - Р. 465-476. DOI: 10.1016/S0020-7403(62)80008-3
23. Gamer U., Sayir M. Elastic-plastic stress distribution in a rotating solid shaft // Zeitschrift fur Angewandte Mathematik und Physik ZAMP. - 1984. - Bd. 35. - Nr. 5. - S. 601-617. DOI: 10.1007/BF00952107
24. Ивлев Д.Д. Три дискуссии по механике // Вестн. СамГУ. - 2007. - Т. 54, № 4. - С. 115-123.
25. Gamer U., Mack W., Varga I. Rotating elastic-plastic solid shaft with fixed ends~// International Journal of Engineering Science. - 1997. - Vol. 35. - No. 3. - Р. 253-267. DOI: 10.1016/S0020-7225(96)00085-7
26. Lindner T., Mack W. Residual stresses in an elastic-plastic solid shaft with fixed ends after previous rotation // ZAMM Zeitschrift fur Angewandte Mathematik und Mechanik. - 1998. -Bd. 78. - Nr. 2. - S. 75-86. DOI: 10.1002/(SICI) 1521-4001(199802)78:2<75::AID-ZAMM75>3.0.m;2-V
27. Gamer U., Lance R.H. Stress distribution in a rotating elastic-plastic tube // Acta Mechanica. - 1983. - Vol. 50. -No. 1-2. - Р. 1-8. DOI: 10.1007/BF01170437
28. Mack W. The rotating elastic-plastic solid shaft with free ends // Tech. Mech. - 1991. - No. 12. - Р. 119-124.
29. Mack W. Rotating Elastic-plastic tube with free ends // International Journal of Solids and Structures. - 1991. - Vol. 27. -No. 11. - Р. 1461-1476. DOI: 10.1016/0020-7683(91)90042-E
30. Rees D.W.A. Elastic-Plastic Stresses in Rotating Discs by von Mises and Tresca // ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik. - 1999. - Vol. 79. - No. 4. - Р. 281-88. DOI: 10.1002/(SICI)1521-4001(199904)79:4<281: :AID-ZAMM281>3.0.CO;2-V.
31. Aleksandrova N. Exact Deformation Analysis of a Solid Rotating Elastic-Perfectly Plastic Disk // International Journal of Mechanical Sciences. - 2014. - Vol. 88. - Р. 55-60. DOI: 10.1016/j.ijmecsci.2014.07.003
32. Eraslan A.N. On the linearly hardening rotating solid shaft // European Journal of Mechanics - A/Solids. - 2003. -Vol. 22. - No. 2. - Р. 295-307. DOI: 10.1016/S0997-7538(02)00002-5
33. Eraslan A.N. Von Mises' yield criterion and nonlinearly hardening rotating shafts // Acta Mechanica. - 2004. - Vol. 168. -No. 3-4. - Р. 129-144. DOI: 10.1007/s00707-004-0088-z
34. Akis T., Eraslan A.N. Exact solution of rotating FGM shaft problem in the elastoplastic state of stress // Arch Appl Mech. - 2007. - Vol. 77. - No. 10. - Р. 745-765. DOI: 10.1007/s00419-007-0123-3
References
1. Tolokonnikov L.A. Mekhanika deformiruemogo tverdogo tela [Mechanics of deformable solids]. Moscow, Vysshaia shkola, 1979, 318 p.
2. Rabotnov Iu.N. Polzuchest' elementov konstruktsii [Creep of structural elements]. Moscow, Nauka, 1966, 752 p.
3. Rabotnov Iu.N. Mekhanika deformiruemogo tverdogo tela [Mechanics of deformable solids]. Moscow, Nauka, 1979, 744 p.
35. Argeso H., Eraslan A.N. A Computational study on functionally graded rotating solid shafts // International Journal for Computational Methods in Engineering Science and Mechanics. -2007. - Vol. 8. - No. 6. - Р. 391-399. DOI: 10.1080/15502280701577842
36. Eraslan A.N., Arslan E. Plane strain analytical solutions to rotating partially plastic graded hollow shafts // Turkish Journal of Engineering and Environmental Sciences. - 2007. - Vol. 31. -No. 5. - Р. 273-287.
37. Nejad M.Z., Fatehi P. Exact elasto-plastic analysis of rotating thick-walled cylindrical pressure vessels made of functionally graded materials // International Journal of Engineering Science. -2015. - Vol. 86. - P. 26-43. DOI: 10.1016/j.ijengsci.2014.10.002
38. Прокудин А.Н., Фирсов С.В. Вязкопластическое течение вращающегося полого цилиндра // Дальневост. мат. журн. - 2018. - Т. 18, № 2. - С. 242-260.
39. Yuan R., Qian L., Gu K. Method of rotating autofrettage and analysis of stress for thick-walled cylinder Nanjing Li Gong Daxue Xuebao // Journal of Nanjing University of Science and Technology. - 1997. - Vol. 21. - No. 1. - P. 13-16.
40. Zare H.R., Darijani H. A novel autofrettage method for strengthening and design of thick-walled cylinders // Materials & Design. - 2016. - Vol. 105. - P. 366-374. DOI: 10.1016/j.matdes.2016.05.062
41. Zare H.R., Darijani H. Strengthening and design of the linear hardening thick-walled cylinders using the new method of rotational autofrettage // International Journal of Mechanical Sciences. - 2017. -Vol. 124-125. - P. 1-8. DOI: 10.1016/j.ijmecsci.2017.02.015
42. Kamal S.M. Analysis of Residual Stress in the Rotational Autofrettage of Thick-Walled Disks // Journal of Pressure Vessel Technology, Transactions of the ASME. - 2018. - Vol. 140. -No. 6. - P. 061402-061402-10. DOI: 10.1115/1.4041339
43. Kamal S.M., Perl M., Bharali D. Generalized plane strain study of rotational autofrettage of thick-walled cylinders-Part I: Theoretical analysis // Journal of Pressure Vessel Technology, Transactions of the ASME. - 2019. - Vol. 141. - No. 5. -P. 051201-051201-11. DOI: 10.1115/1.4043591
44. Güven U. The fully plastic rotating disk with rigid inclusion // Z. Angew. Math. Mech. - 1997. - Vol. 77. - No. 9. -P. 714-716. DOI: 10.1002/zamm.19970770912
45. Güven U. Elastic-plastic stress distribution in a rotating hyperbolic disk with rigid inclusion // International Journal of Mechanical Sciences. - 1998. - Vol. 40. - No. 1. - P. 97-109. DOI: 10.1016/S0020-7403(97)00036-2
46. Güven U. Elastic-plastic rotating disk with rigid inclusion // Mechanics of Structures and Machines. - 1999. - Vol. 27. -No. 1. - P. 117-128. DOI: 10.1080/08905459908915691
47. Güven U., Parmaksizoglu C., Altay O. Elastic-plastic rotating annular disk with rigid casing // ZAMM - Journal of Applied Mathematics and Mechanics. Zeitschrift für Angewandte Mathematik und Mechanik. - 1999. - Bd. 79. - Nr. 7. - S. 499503. DOI: 10.1002/(SICI)1521-4001(199907)79:7<499::AID-ZAMM499>3.0.CO;2-C
4. Timoshenko S.P., Gud'er Dzh. Teoriia uprugosti [Theory of elasticity]. Moscow, Nauka, 1979, 576 p
5. Lur'e A.I. Nelineinaia teoriia uprugosti [Nonlinear theory of elasticity]. Moscow, Nauka, 1980, 512 p.
6. Nadai A. Theory of Flow and Fracture of Solids, Volume One. 2nd Edition. New York, McGraw Hill, 1950, 572 p.
7. Chakrabarty, J. Theory of Plasticity: Third Edition. Elsevier, 2005, 896 p. DOI: 10.1016/C2009-0-05774-2
8. Sadd M.H. Elasticity: theory, applications, and numerics. Elsevier, 2014, 582 p. DOI: 10.1016/C2012-0-06981-5
9. Timoshenko S. Theory of elasticity. New York: McGraw Hill, 2010, 567 p.
10. Zenkour, A.M. Analytical solutions for rotating exponentially-graded annular disks with various boundary conditions. International Journal of Structural Stability and Dynamics, 2005, vol. 05, no. 04, pp. 557-577. DOI: 10.1142/S0219455405001726.
11. Eraslan A.N., Akis T. On the plane strain and plane stress solutions of functionally graded rotating solid shaft and solid disk problems. Acta Mechanica, 2006, vol. 181, no. 1. pp. 43-63.
12. Akis T., Eraslan A.N. The stress response and onset of yield of rotating FGM hollow shafts. Acta Mechanica, 2006, vol. 187, no. 1, pp. 169-187.
13. Sburlati R. Analytical elastic solutions for pressurized hollow cylinders with internal functionally graded coatings. Composite Structures, 2012, vol. 94, no. 12, pp. 3592-3600.
14. Zenkour, A.M., Elsibai K.A., Mashat D.S. Elastic and Viscoelastic Solutions to Rotating Functionally Graded Hollow and Solid Cylinders. Applied Mathematics and Mechanics, 2008, vol. 29, no. 12, pp. 1601-1616. DOI: 10.1007/s10483-008-1208-x
15. Horgan C.O., Saccomandi G. Large Deformations of a Rotating Solid Cylinder for Non-Gaussian Isotropic, Incompressible Hyperelastic Materials. J. Appl. Mech, 2000, vol. 68, no. 1, pp. 115-117. DOI: 10.1115/1.1349418
16. Bagheri A., Taghizadeh D., Darijani H. On the behavior of rotating thick-walled cylinders made of hyperelastic materials. Meccanica, 2016, vol. 51, no. 3, pp. 673-692. DOI: 10.1007/s11012-015-0233-x
17. Chadwick P., Creasy C.F.M., Hart V.G. The deformation of rubber cylinders and tubes by rotation. The ANZIAM Journal, 1977, vol. 20, no. 1, pp. 62-96. DOI: 10.1017/S0334270000001454
18. Nemirovskii Iu.V., Iankovskii A.P. Otsenka predel'noi skorosti vrashcheniia iskrivlennykh metallokompozitnykh diskov v usloviiakh polzuchesti [Estimation of the limiting rotation speed of curved metal composite discs under creep conditions]. I. Yakovlev Chuvash State Pedagogical University Bulletin, 2010, no. 2 (66), pp. 102-111.
19. Prokudin A.N., Firsov S.V. Creep Analysis of Rotating Cylinder with Free Ends. Bulletin of the Yakovlev Chuvash State Pedagogical University. Series: Mechanics of Limit State, 2018, no 1 (35), pp. 63-73.
20. Banshchikova I.A. The stress-strain state and duration until fracture of rotating disks in creep. PNRPU Mechanics Bulletin, 2018, no. 4, pp. 20-32. DOI: 10.15593/perm.mech/2018.4.02
21. Begun A.S., Kovtanuk L.V. Calculation of Stresses, Strains and Displacements in a Rotating Disk under Creep Conditions. Bulletin of the Yakovlev Chuvash State Pedagogical University. Series: Mechanics of Limit State, 2019, no. 1 (39). pp. 84-93. DOI: 10.26293/chgpu.2019.39.1.011
22. Hodge P.G., Balaban M. Elastic-plastic analysis of a rotating cylinder. International Journal of Mechanical Sciences, 1962, vol. 4, no. 6, pp. 465-476. DOI: 10.1016/S0020-7403(62)80008-3
23. Gamer U., Sayir M. Elastic-plastic stress distribution in a rotating solid shaft. Zeitschrift Fur Angewandte Mathematik Und Physik ZAMP, 1984, vol. 35, no 5, pp. 601-617. DOI: 10.1007/BF00952107
24. Ivlev D.D. Tri diskussii po mekhaniki [Three discussions on mechanics]. VestnikSamGU, 2007, vol. 54, no 4, pp. 115-123.
25. Gamer U., Mack W., Varga I. Rotating elastic-plastic solid shaft with fixed ends. International Journal of Engineering
Science, 1997, vol. 35, no. 3. pp. 253-267. DOI: 10.1016/S0020-7225(96)00085-7
26. Lindner T., Mack W. Residual stresses in an elastic-plastic solid shaft with fixed ends after previous rotation. ZAMM Zeitschrift fur Angewandte Mathematik und Mechanik. 1998, vol. 78, no. 2, pp. 75-86. DOI: 10.1002/(SICI)1521-4001(199802)78:2<75::AID-ZAMM75>3.0.C0;2-V
27. Gamer U., Lance R.H. Stress distribution in a rotating elastic-plastic tube. Acta Mechanica, 1983, vol. 50, no. 1-2. pp. 1-8. DOI: 10.1007/BF01170437
28. Mack W. The rotating elastic-plastic solid shaft with free ends. Tech. Mech., 1991, no. 12, pp. 119-124.
29. Mack W. Rotating Elastic-plastic tube with free ends. International Journal of Solids and Structures, 1991, vol. 27, no. 11, pp. 1461-1476. DOI: 10.1016/0020-7683(91)90042-E
30. Rees, D.W.A. Elastic-Plastic Stresses in Rotating Discs by von Mises and Tresca. ZAMM - Journal of Applied Mathematics and Mechanics, 1999, vol. 79, no. 4, pp. 281-288. DOI: 10.1002/(SICI)1521-4001(199904)79:4<281: :AID-ZAMM281>3.0.CO;2-V.
31. Aleksandrova, N. Exact Deformation Analysis of a Solid Rotating Elastic-Perfectly Plastic Disk. International Journal of Mechanical Sciences, 2014, vol. 88, pp. 55-60. DOI: 10.1016/j.ijmecsci.2014.07.003
32. Eraslan A.N. On the linearly hardening rotating solid shaft. European Journal of Mechanics - A/Solids, 2003, vol. 22, no. 2, pp. 295-307. DOI: 10.1016/S0997-7538(02)00002-5
33. Eraslan A.N. Von Mises' yield criterion and nonlinearly hardening rotating shafts. Acta Mechanica, 2004, vol. 168, no. 3-4, pp. 129-144. DOI: 10.1007/s00707-004-0088-z
34. Akis T., Eraslan A.N. Exact solution of rotating FGM shaft problem in the elastoplastic state of stress. Arch Appl Mech., 2007, vol. 77, no. 10, pp. 745-765. DOI: 10.1007/s00419-007-0123-3
35. Argeso H., Eraslan A.N. A Computational Study on Functionally Graded Rotating Solid Shafts. International Journal for Computational Methods in Engineering Science and Mechanics, 2007, vol. 8, no. 6, pp. 391-399. DOI: 10.1080/15502280701577842
36. Eraslan A.N., Arslan E. Plane strain analytical solutions to rotating partially plastic graded hollow shafts. Turkish Journal of Engineering and Environmental Sciences, 2007, vol. 31, no. 5, pp. 273-287.
37. Nejad, M.Z., Fatehi P. Exact elasto-plastic analysis of rotating thick-walled cylindrical pressure vessels made of functionally graded materials. International Journal of Engineering Science, 2015, vol. 86, pp. 26-43. DOI: 10.1016/j.ijengsci.2014.10.002.
38. Prokudin A.N., Firsov S.V. Viazkoplasticheskoe techenie vrashchaiushchegosia pologo tsilindra [Viscoplastic flow in a rotating hollow cylinder]. Far Eastern Mathematical Journal, 2018, vol. 18, no. 2, pp. 242-260.
39. Yuan R., Qian L., Gu K. Method of rotating autofrettage and analysis of stress for thick-walled cylinder. Journal of Nanjing University of Science and Technology, 1997, vol. 21, no. 1, pp. 13-16.
40. Zare, H. R., Darijani H. A novel autofrettage method for strengthening and design of thick-walled cylinders. Materials & Design, 2016, vol. 105, pp. 366-374. DOI: 10.1016/j.matdes.2016.05.062.
41. Zare H.R., Darijani H. Strengthening and design of the linear hardening thick-walled cylinders using the new method of rotational autofrettage. International Journal of Mechanical Sciences, 2017, vol. 124-125, pp. 1-8. DOI: 10.1016/j.ijmecsci.2017.02.015
42. Kamal S.M. Analysis of Residual Stress in the Rotational Autofrettage of Thick-Walled Disks. Journal of Pressure Vessel Technology, Transactions of the ASME, 2018, vol. 140, no. 6. pp. 061402-061402-10. DOI: 10.1115/1.4041339
43. Kamal S.M., Perl M., Bharali D. Generalized plane strain study of rotational autofrettage of thick-walled cylinders-Part I: Theoretical analysis. Journal of Pressure Vessel Technology, Transactions of the ASME, 2019, vol. 141, no. 5, pp. 051201051201-11. DOI: 10.1115/1.4043591
44. Guven U. The fully plastic rotating disk with rigid inclusion. Z. Angew. Math. Mech., 1997, vol. 77, no. 9, pp. 714-716. DOI: 10.1002/zamm. 19970770912
45. Guven U. Elastic-plastic stress distribution in a rotating hyperbolic disk with rigid inclusion. International Journal of Mechanical Sciences, 1998, vol. 40, no. 1, pp. 97-109. DOI: 10.1016/S0020-7403(97)00036-2
46. Guven U. Elastic-Plastic Rotating Disk with Rigid Inclusion. Mechanics of Structures and Machines, 1999, vol. 27, no. 1, pp. 117-128. DOI: 10.1080/08905459908915691
47. Guven U., Parmaksizoglu C., Altay O. Elastic-47. Plastic Rotating Annular Disk with Rigid Casing. ZAMM - Journal of Applied Mathematics and Mechanics, 1999, vol. 79, no. 7, pp. 499-503. DOI: 10.1002/(SICI)1521-4001(199907)79:7<499:: AID-ZAMM499>3.0.CO;2-C