Механика деформируемого твердого тела
DOI: http://www.dx.doi.org/10.24866/2227-6858/2019-4-2 УДК 539.374
А.Н. Прокудин, С.В. Фирсов
ПРОКУДИН АЛЕКСАНДР НИКОЛАЕВИЧ - к.т.н., ведущий научный сотрудник, AuthorID: 415431, SPIN: 6812-2451, ORCID: 0000-0002-5156-424X, ResearcherID: N-9344-2016, ScopusID: 35722777500, e-mail: [email protected] ФИРСОВ СЕРГЕЙ ВИКТОРОВИЧ - младший научный сотрудник, AuthorID: 783315, SPIN: 8267-2329, ORCID: 0000-0001-7446-6231, ResearcherID: D-1966-2018, ScopusID: 56976208300, e-mail: [email protected] Институт машиноведения и металлургии ДВО РАН Металлургов ул., 1, Комсомольск-на-Амуре, Россия, 681005
Упругопластическое деформирование вращающегося полого цилиндра с жестким покрытием на внутренней и внешней стенках
Аннотация: Рассматривается полый цилиндр из идеального изотропного упругопластического материала в условиях плоской деформации с ограничениями по радиальному перемещению на внутренней и внешней поверхностях. Предполагается, что полные деформации в цилиндре являются малыми и складываются из упругих и пластических деформаций. Напряжения и упругие деформации связаны законом Гука, пластические деформации определяются условием Треска и ассоциированным с ним законом пластического течения. Деформирование цилиндра происходит за счет центробежных сил, возникающих при вращении цилиндра с угловой скоростью, которая сначала увеличивается до выбранного максимального значения, а затем уменьшается вплоть до полной остановки цилиндра. Предполагается, что скорость вращения медленно изменяется со временем и угловое ускорение пренебрежимо мало. Установлено, что пластическое течение возникает одновременно на внешней и внутренней поверхностях цилиндра. Скорость вращения цилиндра, соответствующая началу пластического течения, найдена из решения упругой задачи. Для каждой пластической области найдено точное аналитическое решение определяющей системы уравнений. Приведены системы условий на границах между областями, обеспечивающие непрерывность полученных решений во всем цилиндре. Результаты расчетов проиллюстрированы графиками перемещений, напряжений и пластических деформаций.
Ключевые слова: вращающийся цилиндр, плоская деформация, упругопластичность, пластическое течение, закон Гука, условие Треска, повторное пластическое течение, остаточные напряжения.
Введение
Расчет напряженно-деформированного состояния вращающихся цилиндров и дисков относится к классическим задачам механики деформируемого твердого тела. В научной литературе опубликовано множество работ, посвященных решению данного класса задач для различных кинематических допущений и моделей материалов. Как правило, используется
© Прокудин А.Н., Фирсов С.В., 2019
О статье: поступила: 30.10.2019; финансирование: работа выполнена в рамках государственного задания ИМиМ ДВО РАН № 075-00414-19-00.
один из следующих вариантов кинематики: плоское напряженное состояние, плоская деформация и обобщенная плоская деформация. Исследуется как обратимое (упругое), так и необратимое (пластическое течение и ползучесть) деформирование вращающихся цилиндров и дисков. В последние годы значительное внимание уделяется неоднородным и анизотропным материалам.
Во многих случаях представляют интерес материалы, которые проявляют некоторую комбинацию перечисленных свойств, например упругопластические материалы. Вращающийся сплошной цилиндр из упругопластического материала впервые изучался в [8], где использовалось условие Треска и ассоциированный закон пластического течения, а также отдельно рассматривался случай конечных деформаций. Авторы предполагали, что в цилиндре возникает единственная пластическая область, в которой радиальное и тангенциальное напряжения совпадают и превышают осевое напряжения на величину предела текучести. Однако, как позднее было показано в [5], данное предположение приводит к разрыву перемещений на упругопластической границе. Корректное решение данной задачи опубликовано в [4, 5]. Авторами показано, что условие Треска впервые выполняется в центре цилиндра, что приводит к одновременному появлению не одной, а двух пластических областей, первая из которых соответствует ребру призмы Треска, а вторая - ее грани. Последующее увеличение скорости вращения приводит к развитию пластического течения и уменьшению упругой области. Решение [4, 5] является непрерывным во всем цилиндре на всех стадиях пластического течения. В [9] на основе результатов [4, 5] изучается распределение остаточных деформаций в сплошном цилиндре после его вращения. Авторами также установлено, что если скорость вращения цилиндра была достаточно высокой, то при уменьшении скорости в нем может возникнуть повторное пластическое течение.
Полый цилиндр с закрепленными концами исследовался в [3]. Авторы рассматривали как идеально-пластический, так и линейно-упрочняющийся материал. Установлено, что пластическое течение происходит только в одной области, которая зарождается на внутренней поверхности цилиндра и с увеличением скорости вращения распространяется на весь цилиндр.
Исследование упругопластического деформирования в цилиндрах со свободными концами осложняется дополнительным условием на суммарную осевую силу. В [11] получено решение для вращающегося сплошного цилиндра, а в [10] - для полого. Сравнение результатов [10, 11] и [4, 5, 9] показывает, что пластическое течение во вращающихся цилиндрах с закрепленными и свободными концами характеризуется различной конфигурацией пластических областей. Заметим, что в отмеченных работах [3-5, 8-11] использовалась модель идеального упругопластического материала. Влияние вязких свойств материала на необратимое деформирование вращающихся дисков и цилиндров изучалось в работах [1, 2].
Детальный обзор научной литературы, посвященной расчету упругопластических деформаций во вращающихся цилиндрах и дисках, показал, что наиболее часто рассматриваются два случая: это сплошной либо полый цилиндр (диск), на боковых поверхностях которого заданы напряжения. Другие виды граничных условий, при которых на одной или двух граничных поверхностях вращающегося полого цилиндра задано перемещение, рассматриваются значительно реже, несмотря на их значительный прикладной интерес. В [6] получено решение для вращающегося диска с жестким включением, а в [7] - для вращающихся дисков с жестким покрытием на внешней поверхности, при этом внутренняя поверхность диска могла быть как свободной, так и зафиксированной от радиальных перемещений.
Настоящая статья посвящена исследованию вращающегося полого цилиндра из идеального упругопластического материала с жестким покрытием на внутренней и внешней граничных поверхностях. Подобные граничные условия на практике могут возникнуть в тех случаях, когда исходная деталь состоит из композитного материала, в котором внешняя и внутренняя часть состоят из более жесткого материала, по сравнению с материалом центральной части. Предполагается, что цилиндр находится в плоском деформированном состо-
янии, а деформации в нем являются малыми. Для постановки задачи используется условие Треска и ассоциированный с ним закон пластического течения. Мы рассмотрим стадии нагружения и разгрузки, включая повторное пластическое течение. Ранее подобная задача не решалась.
Постановка задачи
Рассмотрим бесконечно длинную цилиндрическую среду с упругими и пластическими свойствами. Цилиндр ограничен двумя поверхностями г = а и г = Ь (Ь > а). Будем рассматривать случай, когда на двух граничных поверхностях наложено жесткое ограничение на радиальную компоненту перемещений. Положим, что такой цилиндр вращается с угловой скоростью а(). В начальный момент времени а(0) = 0. В дальнейшем скорость вращения постепенно увеличивается до значения а = атах. После достижения этого значения скорость начинает плавно снижаться, вплоть до полной остановки вращения. В цилиндрической системе координат (г, в, г) получим одномерный процесс деформирования относительно координаты г. Для большей общности приведем ее к безразмерному виду
р = Г, 8 = а, 0 <8<р< 1.
Ь Ь
Деформации в материале будем считать малыми и состоящими из упругих и пластических. Ненулевыми компонентами деформаций в нашем случае будут только компоненты ёгг и ёвв тензора полных деформаций:
7 ди , и ,
йгг = егг + ргг = —, йвв= евв + рвв=-, = егг + р22 = 0, (1)
ер р
где егг, евв, е22 - компоненты тензора упругих деформаций; ргг, рвв, р22 - компоненты тензора пластических деформаций; и = иг/Ь - безразмерные перемещения.
Напряжения в деформируемом слое определяются упругими деформациями в соответствии с законом Гука и в безразмерном виде примут значения
1 = Е--^-г ((1 - V) егг + Vе вв + ),
1 =
°0 ( 1 + V )(1 - -IV )
Е 1
°0 ( 1 + V )(1 - -IV )
Е 1
° (1 + V)(1 -IV)
'вв
+ ve,
:), (2)
(уегг + Veвв+(1 - V) егг),
где Е - модуль Юнга, V - коэффициент Пуассона, о0 - предел текучести.
Соотношения, обратные к (2), имеют вид
егг = 1 1гг - - ), ^ = у- - ° ), егг = у - 1 - ™вв ). (3)
Подставив (3) в (1), получим соотношения для определения пластических деформаций через перемещения и напряжения
ди °о ( \ и °о ( \
др & р Е ^
Ргг =-1 - - ™вв ) • Е
Будем считать угловое ускорение достаточно малым, чтобы можно было пренебречь силой инерции, определяемой вращательным ускорением. Это позволяет оставаться в рамках одномерности. В таком случае единственное уравнение равновесия примет вид
°вв =
( + (^( = _0Р, 0 = > о, (5)
дР Р (Го
где р - плотность материала, О = 0.(1) — параметр нагрузки.
Для учета пластических свойств деформируемого слоя воспользуемся теорией пластического течения. Для задания поверхности нагружения возьмем условие пластичности Треска
/ = таХ (|(г _ , \Ггт _(22\ , |) _1 = 0 . (6)
Ассоциированный закон пластического течения для граней призмы Треска примет вид
^ = , , ^ = , (7)
Л д( Л дгвв Л д(
где dX - положительный коэффициент.
Использование закона (7) вместе с условием (6) приводит к пластической несжимаемости. Из чего следует, что объемная деформация является чисто упругой. Исходя из этого и соотношений (4), получим
^ = (Ч1_>)(СТ„ (+(.). (8)
Граничные условия запишутся в виде
и (£) = 0, и (1) = 0. (9)
Весь процесс деформирования цилиндра состоит из нагрузки (рост скорости вращения и, соответственно, параметра нагрузки) и разгрузки (уменьшение скорости вращения и параметра нагрузки). При нагрузке сначала происходит упругое деформирование. Потом, при достижении напряжениями поверхности нагружения, в цилиндре образуются области пластического течения. Впоследствии пластическое течение распространяется на всю область деформирования. После достижения параметром нагрузки заданного максимального значения начинается стадии разгрузки. Во время этой стадии сначала происходит упругое деформирование. При дальнейшем уменьшении значения параметра нагрузки вследствие накопленных пластических деформаций образуются области повторного пластического течения. Рассмотрим более подробно указанные стадии.
Решение
Упругое деформирование
Из соотношений (1), (2), (5), с учетом равенства нулю пластических деформаций, получим дифференциальное уравнение, решение которого примет вид
1 А 1 Г (1 + ^)(1 _2у)
и = 1ДР + А_-^-1ПР3, (10)
2 Р 8 Е 1 _ V
где П2 - константы интегрирования. Напряжения примут значения
^ =Е ДР2 _ 2А (1 _ 2у) 13 _ 2у 2^
( Го 2Р2 (1 + У)(1 _2У) 8 1 -V '
Е ДР2 + 2Д (1 _ 2У) 11 + 2^
Г0в=— пп2(л^ у 0 Ч -Р ^ ( =У(Ггг + Гвв) •
Г0 2Р (1 + v)(1 _2У) 8 1 _у
Для текущих граничных условий (9) константы интегрирования А, П2 будут равны
д=_1а(1+м)(1-У)У
1 4Е 1-V 1 ' 8Е 1 -V
(12)
Полученное упругое решение (10)-(12) справедливо при О < О1. Пластическое течение начинается одновременно на двух граничных поверхностях. На внутренней граничной поверхности (в = Ь) напряжения удовлетворяют неравенству агг > авв = а22, а на внешней (в = 1) - авв = а22 > агг. Условие пластичности (6) примет вид
а -ааа = 1, а -а = 1
гг вв ' гг 22
на внутренней границе и
ааа-а = 1, а -а = 1
вв гг 22 гг
на внешней.
Используя одно из этих условий и упругое решение (11), (12), получим следующее значение параметра нагрузки, соответствующее началу пластического течения:
1 - V 4
^ = , _ - .
1 - 2М1 -У
На (рис. 1) формула (13) представлена графически при разных значениях V.
(13)
Рис. 1. Зависимость параметра нагрузки от 5 при разных V.
Упругопластическое деформирование
При О = О1 в цилиндре появляются 4 области пластического течения. Обозначим их как области I, II, III и IV. Напряжения в области I (Ь < в < в1) соответствуют ребру призмы Треска агг > авв = а22. В области II (в1 < в < в 2) - грани призмы Треска агг > авв > а22. В области III (в3 < в < в4) - грани авв > а22 > агг. В области IV (в4 < в < 1) - ребру авв = а22 > агг. В области в2 < в < вз продолжается упругое деформирование.
С увеличением значения параметра нагрузки до О = О2 на упругопластической границе в = в2 происходит выполнение равенства агг = авв. Это приводит к зарождению двух новых областей пластического течения V и VI. Область V (в2 < в < в5) соответствует ребру призмы Треска агг = авв > а22. Область VI (в5 < в < вб) соответствует грани авв > агг > а22. В области вб < в < в3 продолжается упругое деформирование.
При О = О3 происходит совмещение границ вб и в3. Одновременно с этим на месте их совмещения (в = вб = вз) выполняется равенство агг = а22. Это приводит к тому, что область упругого деформирования между областями VI и III исчезает и на ее месте появляется область пластического течения VII. Область VII (вб < в < вз) соответствует ребру призмы Треска авв > агг = а22.
При дальнейшем увеличении значения параметра нагрузки вплоть до О = Отах области пластического течения II, III, VI, соответствующие граням призмы Треска, постепенно
уменьшаются. В это время области пластического течения I, IV, V, VII, находящиеся между ними и соответствующие ребрам призмы Треска, постепенно расширяются.
Рассмотрим каждую область более подробно.
Упругая область
Полученное ранее упругое решение (10), (11) остается справедливым, однако значение констант интегрирования П2 необходимо определять из условий на новых границах упругой области. При 0у < О < 02 это границы в = в2 и в = в3. Затем при 02 < О < 03 это границы в = вб и в = в3.
Область I
Напряжения соответствуют ребру призмы Треска огг > овв = о22. Условие пластичности (6) принимает вид
— _— = 1, — _ — = 1
е гг вв ' гг 22
С помощью данного условия и уравнения (5) найдем значения напряжений
Ггг = С^Р^-МР, Гвв=(22 =Ггг"1 , (14)
где Су — константа интегрирования. В дальнейшем за С^ будем обозначать константы интегрирования для областей пластического течения на стадии нагрузки.
Подставив данные напряжения (14) в (8), получим дифференциальное уравнение для перемещений, решение которого имеет вид
и = 3—(1 _2^СР + С2_3—(1 _2^Р3 _ 1 —(1 _2*0(1 + 61пР)Р . (15)
Из уравнений (4), (14), (15) можно найти значения пластических деформаций.
Область II
Напряжения соответствуют грани призмы Треска огг > овв > о22. Условие пластичности (6) принимает вид
— _ — = 1
гг 22
Из ассоциированного закона пластического течения (7), с учетом (1) и данного условия пластичности можно выразить полные деформации через упругие. Из них, с помощью (3) и условия пластичности, можно найти напряжения
1 + 1 Е й„ + 2^вв __ Е Уйгг +(1 _у)йвв (16)
— = 2 + 2 ((1 )(1 _2у) ' Гвв=( (1 + ^)(1 _ 2^) , ^ =(гг _1 . (16)
Подставив полученные напряжения в уравнение равновесия (5), с учетом (1), получим дифференциальное уравнение, решение которого примет вид
= с Р+ с Р^ + — (1 + V) Р _ 2—(1 + V)(1 2у) ОР3.
3Р , С4Р , .1 , К Р 2 ОР . (17)
Е Е 7 + 2v
Из (1), (16), (17) можно получить итоговые выражения для напряжений. Из них и уравнений (4), (17) можно найти значения пластических деформаций.
Область III
Напряжения соответствуют грани призмы Треска овв > о22 > огг. Условие пластичности (6) принимает вид
°ee-°rr = 1
С помощью данного условия и уравнения (5) найдем напряжения 1
J = C5 - — ß2Q + ln ß, Oee= 1 + G„.
2 (18) Из ассоциированного закона пластического течения (7) найдем зависимости для компонент тензора пластических деформаций. С учетом полученных зависимостей из (1) следует, что ezz = 0. В результате осевое напряжение azz определяется так же, как и в случае упругости (11). Подставив полученные напряжения (18) в (8), получим дифференциальное уравнение для перемещений, решение которого имеет вид
u = J (1 + V)(1-2V) Cß + C6-—J (1 + v)(l-2v)Qß3 + E (1 + v)(l-2v)ßln ß
Пластические деформации можно найти из (4), с учетом (18)—(19). Область IV
Напряжения соответствуют ребру призмы Треска овв = azz > arr. По аналогии с областью I получим напряжения
1 9
E = C7--ß2n + ln ß, Oee=Oz =1 + E
2 (20)
и перемещения
u = 3 J (1-2v) C7ß + ^-3^l(1-2v )0ß +1^^(1-2v )(1 + 6ln ß)ß.
(21)
Область V
Напряжения соответствуют ребру призмы Треска arr = овв > azz. Условие пластичности (6) принимает вид
jrr=1 Jee-jzz =1. (22)
По аналогии с областью I получим напряжения
1 т
E = C9-- ß2n, aee=arr, E = E-1 (23)
и перемещения
u = 3 E (1-2v) Cß ß + ^l-lE (1-2v )ß-(1-2v )Qß3
2 E V ' ß 2 E V 8 E V ' . (24)
Область VI
Напряжения соответствуют грани призмы Треска овв > arr > azz. Условие пластичности (6) принимает вид
Eee-E =1.
По аналогии с областью II выразим напряжения через полные деформации, используя данное условие пластичности и соотношения (1), (3) и (7)
_Е_ (1 -v) drr +vdee _1 1 2vdrr + dee
E (1 + v)(1-2v) ' Eee" 2 + 2 j(1 + v)(1-2v)' E "Eee-1. (25)
Подставив напряжения (25) в уравнение равновесия (5), с учетом (1), получим дифференциальное уравнение относительно перемещений. Его решение примет вид
и = Сп р^ + С12 Р^ + — (1 + м)Р- 2—(1 + м)(1 - 2м) ор3
11 12 Ev ' Е 17 - 18М . (26)
Из (1), (25), (26) можно получить окончательные формулы для напряжений. Пластические деформации найдем, подставив полученные напряжения и перемещения в (4).
Область VII
Напряжения соответствуют ребру призмы Треска авв > агг = а22. Условие пластичности (6) примет вид
—вв-—гг = 1 —вв-—22 = 1 .
По аналогии с областью I получим напряжения 1 2
—гг =—22 = С13 -2Р°+ 1П Р —вв= 1 + —ГГ ,
(27)
и перемещения
и = 3—(1 - ) С13Р + С4 - 3—(1 - )Ор3 -1 — (1 - )(1 - 61п р)р
Разгрузка
После достижения значения О = Отах скорость вращения а начинает уменьшаться, что в свою очередь приводит к уменьшению значения параметра нагрузки О. В связи с этим уменьшаются напряжения и, соответственно, происходит упругое деформирование цилиндра. На процесс упругости оказывают влияние накопленные пластически деформации. Обозначим их как
Ргг = Ргг (0тах ), Рвв = Рвв (0тах ), Р22 = Р22 (0тах) (29)
Преобразуя уравнение равновесия (5) с учетом (1) и (29), получим дифференциальное уравнение в перемещениях для стадии разгрузки. В общем случае решение данного неоднородного уравнения примет вид
и = ДР + Р +1 р\ рр¿р-а— Р о,
р = -Рх)+%- + М (др«-др-^
р 1 V VГГ г*,, др 1 -М^др др) , (30)
где В], В2 - константы интегрирования.
Подставив значения накопленных пластических деформаций, для каждой из областей I—VII можно получить соответствующие формулы для перемещений. Напряжения можно найти из (2) с учетом (1), значений перемещений и накопленных пластических деформаций.
Повторное пластическое течение
При дальнейшем снижении скорости вращения в цилиндре образуются области повторного пластического течения. В таких областях пластические деформации будут состоять из накопленных (ргг, рвв, р22) и вторичных ^гг, 5вв, ^22) пластических деформаций
Ртт = Ргг + ^' Рев = Рвв + 5вв' Ргг = Ргг + 1^22 • (31)
При уменьшении параметра нагрузки до значения О = О4 в цилиндре зарождаются 4 области повторного пластического течения. Обозначим их VIII, IX, X, XI. Область VIII — в — в7) соответствует ребру призмы Треска овв = о22 > отт, область IX (в7 < в — в8) - грани о22 > овв > отт, область X (в9 — в < в 10) — грани отт > о22 > овв, область XI (в 10 — в — 1) - ребру
Отт > Овв = 022.
При дальнейшем уменьшении параметра нагрузки в момент, когда О = О5, упруго-пластическая граница в = в8 сравняется с предыдущей границей области в = в1. Соответственно, область пластического течения IX разобьем на две подобласти: IX (в7 < в — в1) и XII (в1 < в — в8). Различия между ними заключаются только в накопленных пластических деформациях.
Рассмотрим подробней получившиеся области повторного пластического течения.
Область VIII
Напряжения соответствуют ребру призмы Треска Овв = о22 > отт. Напряжения определяются из (20), перемещения из (21). Константы интегрирования С7, С8 заменяются на А1, А2. В дальнейшем для констант интегрирования областей повторного пластического течения будем использовать обозначения Аг-.
Область IX
Напряжения соответствуют грани призмы Треска о22 > Овв > отт. Условие пластичности (6) примет вид
а —а = 1
22 гг
По аналогии с областью II из ассоциированного закона пластического течения (7) найдем соотношения, связывающие компоненты тензора вторичных пластических деформаций. Из них и соотношений (1), (3) и (31) выразим значения полных деформаций через напряжения и накопленные пластические деформации. Из данной системы, с помощью соотношений (1) и условия пластичности, найдем значения напряжений
11 Е 1
<У„ = а —1 =---+
ГГ ~22 ' 2 2 <0 (1 + у)(1 -2у)ч
ди _ и _ „
— + 2у--Рг — 2,УРвв — Р22
дг г
Е 1 ( ди , ч и , ч „
у^Т + 0 — у)- —уР гг —(1 — у) Рвв—,Р 2.
<о (1 + у)(1 — 2у)
(32)
Подставив напряжения в уравнение равновесия (5), получим дифференциальное уравнение в перемещениях. В общем случае решение данного неоднородного уравнения примет вид
1 п42—2у 1 п—42—2у
и = А +А +1 I Р ар+1 1 ^^ Р ёр—
2у/2 — 2у ) 2^2 — 2у
((1+у)('—2у)
Е (1 — у)(7 + 2у) Е 42—2,
Р = 1(1 — 2у)( Р гг — 2Рвв+ Р 22 ) + % + 2^% + -дР 22
I /^гг ™ д| д| д| (33)
С учетом накопленных пластических деформаций получим
<вв =
u = + Atß^ + 31(1 -2v)Cß + C-21 (' + 'Г) пр -
- 5 ~F (1-2Г- п-f - 310 (3 + 2v+2 (1 - 2v) In ß)ß.
8 E 7 + 2v 4 E v ' (34)
Пластические деформации можно найти, подставив (32), (34) и накопленные пластические деформации в (4).
Область X
Напряжения соответствуют грани призмы Треска arr > azz > авв. Условие пластичности (6) принимает вид
1 -1ее=1.
С помощью данного условия и уравнения равновесия (5) найдем значение двух компонент тензора напряжений. Из ассоциированного закона пластического течения (7) найдем зависимости для компонент тензора вторичных пластических деформаций (srr = —see, szz = 0). Учитывая данные зависимости и значения накопленных пластических деформаций, найдем формулу для третьей компоненты тензора напряжений. В итоге получим
1 E
°rr = A5lnß, 1ее =°rr -' 1 = V(1rr + 1ее)--p
2 1 . (35)
Подставив напряжения (35) в (8), получим дифференциальное уравнение в перемещениях. Его решение, с учетом накопленных деформаций, имеет вид
u = A 1 (' + v)(1 -2v)ß + A + ^ 1 ('-2v)2 Cß-11(' + v)(1 -2v)Qß3 --i 1 (' - 2V)2 Qmaxß3 + \ 1 (' - 2v)(' - 2(' + 4v) in ß)ß.
(36)
Подставив (35), (36) в (4) и учтя накопленные пластические деформации, можно найти значения пластических деформаций.
Область XI
Напряжения соответствуют ребру призмы Треска arr > авв = azz. Напряжения принимают значения (14), перемещения - (15). Константы интегрирования C1, C2 заменяются на A7, A8.
Область XII
Данная область является продолжением области пластического течения IX. Соответственно перемещения определяются из формулы (33) с константами A9, A10 вместо Aj, A^. Из накопленных пластических деформаций в данной области следует, что p = 0. На основе этого формула (33) примет вид
u = Aß ß^ + Д0 ß^ -11 (' + v) ß - 2—(1 + v)(1 - 2v) Qß3
^ 10 EK f E 7 + 2V . (37)
Результаты расчетов
Расчеты производились при следующих физических и геометрических параметрах:
8 = 0,2; v = 0,3; 1 = 250 МПа; E = 21 ГПа; Q = 22
5 5 ' ' y ? ' max
Пластическое течение в цилиндре начинается на внешней и внутренней граничных поверхностях, согласно (13), при О1 = 175/24 ~ 7,2917. При О > О1, в цилиндре имеется 4 области пластического течения и одна область упруго деформирования. В результате получается 10 неизвестных констант интегрирования В1, П2, С1, С2, С3, С4, С5, Сб, С7, С8 и 4 границы областей в и в2, в3, в4. Для их определения составим систему уравнений вида
и1 = 0; и1У (1) = 0;
Д = Д: и1 = и11; Г = Г1; Гвв - ГГП ■ Г вв;
Д = Д2: и11 = ие1; < = Ггг ; Г гг — Г = гг
Д = Д3. ие1 = и111; < = Ггг ; Гвв — Г = гг
Д = Д4. и111 = и1У; г1:: = Ггг ; Г111 гг = ГУ, гг
(38)
где верхний индекс 1-ГУ обозначает принадлежность к области пластического течения под заданным номером, верхний индекс е1 — к области упругого деформирования.
Полученная система линейно зависит от коэффициентов интегрирования и нелинейно — от положения границ. Из первых двух уравнений в каждой строчке данной системы уравнений можно найти аналитические формулы для констант интегрирования. Границы областей можно найти численно из третьих уравнений в строчках со второй по пятую, подставив в них предварительно найденные значения коэффициентов.
Для нахождения значения параметра нагрузки О = О2, при котором появятся области пластического течения V и VI на упругопластической границе в2, разрешим одно из условий пластического течения (22) совместно с системой (38). Получим следующие значения:
П2 = 18,2737; Д = 0,259; Д = 0,4286; Д3 = 0,676; Д = 0,8338.
При О > 02 в цилиндре появятся еще две области пластического течения. То есть 4 константы интегрирования С9, С10, С11, С12 и две границы в5, вб. Для их определения в систему (39) внесем следующие изменения:
Д = Д. Д = Д. Д = Д6.
и11 = иУ; иУ = и*1; и*1 = ие1;
(Г =0 ;
гг гг ?
„У •
Г = Г ;
гг гг
гУ ^е1 ■
Г = Г ;
гг гг
Г вв Гвв;
* *1 Г = Г ;
гг гг ?
е1 е1 1
Гвв Г гг = 1
(39)
При О = 03 происходит совмещение упругопластических границ вб и в3. На месте их совмещения вместо области упругого деформирования появится область VII пластического течения. Следовательно, для нахождения данного значения можно воспользоваться либо условием равенства положения границ в5 и в3, либо равенства компонент напряжений огг и о22 на этих границах. В итоге получим
20,9841; Д = 0,2659; Д = 0,4088; Д = 0,4428; Д=Д= 0,6397; Д= 0,8192.
При дальнейшем росте скорости вращения в систему (38), (39) внесутся изменения
/? /? • „Л „.VII • ~-У1 ■ ~У1 ^У1 ■
Д = Д. и = и ; ГГг=Гуг ; Гу1=Гу1 ;
п п ■ „,У1 „.111- „.т. У11 111
Д = Д . и = и ; Г = Г ; Г = Г .
г ? гг гг ? гг гг
(40)
В итоге при заданном значении параметра нагрузки О = Отах получим следующие положения границ областей
Д= 0,2682; Д= 0,4031; Д= 0,4474; Д = 0,6351; Д= 0,6441; Д= 0,8146.
На (рис. 2) показаны графики распределения напряжений в цилиндре при найденных значениях параметра нагрузки, а также накопленные в итоге пластические деформации и диаграмма развития областей пластического течения.
0.6 в
- Ргг Рве " Р
А
Ргг
/?2 А
Рее
>04
0.03 0.02 0.01 0.00 -0.01 -0.02 0.03
V 1 1 \ 1 1 V 1 \ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^^
/1 ^^ 1 1 1 1 ЦпУ 1 1 1 1 VI V 1 1 \ 1 п ш ; 1у\
Р
0.2
0.4
0.6 д
1.0
0.6 г
д & — А
Ра — & — &
1.0 0.8 0.6 0.4 0.2
^ : ¡ш
Уугуп
/ VI
' V
! ! ¡п ——: !И
10
15
20'
Рис. 2. НДС цилиндра при нагрузке: а - напряжения при П = П?; б - напряжения при П = П2; в - напряжения при П = П3; г - напряжения при П = Птах; д - пластические деформации при П = Птах; е - развитие областей пластического течения.
При снижении скорости вращения полученные 7 областей пластического течения заменятся на 7 областей упругого деформирования. Неизвестными будут 14 констант интегрирования. Положения границ областей останутся неизменными. Составим следующую систему уравнений для нахождения констант интегрирования:
и
(¿) = 0; и" (1) = 0;
= |: и = и ; <1г = <г ; = Ре- и< = и< ; а« гг < = а<г ;
и = и«; < =<<«; и< = и ; а«1 гг Ш . = <гг ;
= Д: и« = и«; = <<г; = 14: и =и <; а111 гг =<«. гг
(41)
е
Для нахождения значения О = О4, при котором появятся области повторного пластического течения, совместно с данной системой разрешим одно из условий пластичности на граничных поверхностях цилиндра. В итоге получим О4 = 7,4167. Отметим, что численно получается О4 = Отах - 2 О1. К сожалению, затруднительно подтвердить данную зависимость аналитически, так как полученная формула для О4 является весьма громоздкой и нелинейно зависит от границ упругих областей.
При О < О4 появляются еще четыре области. В систему (41) внесутся следующие изменения:
ива (б) = 0; иш (1) = 0;
Д = Д: ивЛ = и«; VIII IX а = а ; гг гг ' авв = ^^ 22 ;
Д = Д8: иЖ = и; „.I . а = а ; гг гг аI -22 -а' = 1 гг
Д = Д9: и1в = ^; ^в —X. а = а ; гг гг аш гг -авв =
Д = Д10: ^ = -и*; „■X . а = а ; гг гг -X авв = а^.
(42)
Когда параметр нагрузки достигает значения О = О5, упругопластическая граница в8 достигает границы упругой области в1. Добавив в систему (41)-(42) равенство этих границ в качестве дополнительного условия, найдем данное значение параметра нагрузки
05 = 1,8167; Д = 0,2237; Д=Д = 0,2682; Д = 0,8638; Д0 = 0,9231
При дальнейшем уменьшении скорости образуется подобласть XII, что приводит к следующим изменениям в системе (41)-(42):
О О IX XII Ш XII
Д = Д: и = и ; агг = а гг ;
Г) Г) XII II XII II II II л
Д = Д: и =и ; агг =агг; а-агг =1; (43)
При остановке вращения (параметр нагрузки О = 0) из системы (41)-(43) получим следующие положения границ пластических областей
Д= 0,2295; Д= 0,2905; Д= 0,832; Д10 = 0,9066.
На (рис. 3) приведены графики напряжений, пластических деформаций и перемещений для найденных значений параметра нагрузки. Также показана диаграмма развития областей пластического течения.
Стоит отметить, что чем меньше значение 3, тем меньше скорость (и, соответственно, значение параметра нагрузки), при которой начнется пластическое течение, - и меньше значения накопленных к этому моменту деформаций. Так, получим следующие значения при 3 = 0,1:
0 = 7,0707; 02 = 15,2579; 03 = 18,9173; р„ (б, 03) = 0,0264 > и при 3 = 0,4:
Ц= 8,3333; 02 = 29,5646; 03 = 31,2925; р^ (б, 03) = 0,0455 .
При значении параметров 3 = 0,03 и Отах = 28 удается добиться появления новых областей повторного пластического течения. Одна область соответствует ребру призмы Треска с условиями а22 - агг = 1 и а22 - авв = 1, а вторая - грани с условием а22 - авв = 1. Получим следующие ключевые значения параметра нагрузки:
= 7,0063; 02 = 13,6929; 03 = 18,3197; Отах = 28; ргг (б,Отах) = 0,0479;
04 = 13,9874 = 0тах-20; 05 = 9,0434; 06 = 5,0417; 07 = 5,0225;
08 = 0,6142 = 0тах- 202,
ВЕСТНИК ИНЖЕНЕРНОЙ ШКОЛЫ ДВФУ. 2019. № 4(41)
где 06 - значение параметра нагрузки, при котором р9 = в2; 07 - значение, при котором в 10 = вз; 08 - значение, при котором появляются новые области повторного пластического течения.
Р\ р2 Рз
Офф
Рв Рз
- сТгг Овв - СГ:
2Х
Ра
Р1 Р%~Р\
—1-н-1- 1 и 1 1 и 1 1 и 1 и 1
1
1 1 1 | 1 1
1 1 1 1 I II 1 и 1 и 1 и 1 1 и 1 ^Ч^Ч. 1 и 1 V ут пт ! IV
Р
1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5
—1—1-1—1— -н- —|—I—I-
1 1 1 1 1
1 \ \ 1 1 1 ^^ 1 | 1 1 1 1 1
1 1 11 11 1 1 1 1
11 11 11 11 11 11 11 11 и \ I V 1 Г"4"», \ \ 1 \\ 1
11 11 11 11 1 \Д1 1 1 тч^
VIII1IX 1 1 1 1X 1 XI
Р
0.2
0.4
0.6 а
0.8
1.0
0.2
0.4
0.6 б
0.8
1.0
10
Д е
Рис. 3. НДС цилиндра при разгрузке: а - напряжения при П = П4; б - напряжения при П = П5; в - напряжения при П = 0; г - пластические деформации при П = 0; д - развитие областей повторного пластического течения; е - перемещения в цилиндре при различных П.
Если предположить, что для значения параметра нагрузки, при котором повторное пластическое течение распространится на весь цилиндр, выполняется то же соотношение, что и для значений, при которых появляются области повторного пластического течения О4 и О8, то получим, что максимальное значение параметра нагрузки должно принимать как минимум значение Отах = 203 = 36,6394. Однако для рассмотрения этого случая необходимо использовать модель конечных деформаций.
Можно отметить, что при упругости места, где компонента напряжений агг пересекает другие компоненты напряжений авв и а22, остаются неизменными при росте значения параметра нагрузки. Они соответственно принимают значения
В первом случае (при 3 = 0,2) эти точки принимают значения в = 0,4472 и в = 0,6177. Полученные расчеты свидетельствуют, что они довольно близки к положению мест, где зарождаются области пластического течения V ((в2(02) = 0,4286) и VII (в3(03) = 0,6397).
Для случая 3 = 0,03 получим значения в = 0,1732 и в = 0,5784. Места зарождения областей пластического течения соответственно в2(&2) = 0,1576 и в3(&3) = 0,6073. Также место зарождения области повторного пластического течения соответствует в8(&8) = 0,1576.
Заключение
В рамках модели малых деформаций было получено решение новой задачи упруго-пластического деформирования вращающегося цилиндра с жестким покрытием на внешней и внутренней стенках. Был полностью описан случай деформирования на стадии роста скорости вращения - от упругого деформирования до полного перехода цилиндра в пластическое состояние. Кроме того, частично рассмотрен случай повторного деформирования на стадии уменьшения скорости вращения. Рассмотрено упругое деформирование после того, как цилиндр полностью перешел в пластическое состояние. Рассмотрена стадия появления областей повторного пластического течения на внутренней и внешней граничных поверхностях и стадия зарождения дополнительных областей пластического течения (переход через ребро призмы Треска) внутри деформируемой среды. Полученное решение позволяет получить непрерывное распределение перемещений, напряжений и деформаций в цилиндре в зависимости от текущей скорости вращения и положения границ областей пластического течения.
При рассмотрении нескольких частных случаев были обнаружены определенные закономерности. Одна из них: связь скоростей, при которых зарождаются области повторного пластического течения, со скоростями зарождения областей первичного пластического течения. Другая обнаруженная закономерность связана с возможностью прогнозирования примерного местоположения области, где напряжения внутри среды будут выходить на ребро призмы Треска.
Сравнение полученных результатов с результатами работы [2], в которой рассматривался вращающийся полый цилиндр со свободными граничными поверхностями, позволяет отметить, что значение параметра нагружения О, при котором появляется пластическое течение и при котором весь цилиндр переходит в пластическое состояние, более чем в пять раз выше в цилиндре с жестким покрытием на граничных поверхностях. То есть нанесение жесткого покрытия на обе поверхности позволяет увеличить более чем в два раза значение скорости вращения а, при котором в цилиндре зарождается пластическое течение. Также в то время как в цилиндре со свободными граничными поверхностями пластическое течение распространится на весь деформируемый слой (в статье [2] это Ор = 3,8717), в цилиндре с жестким покрытием на двух граничных поверхностях пластическое течение еще даже не успевает начаться (в нашем случае пластическое течение появляется при ОI = 7,2917).
С практической точки зрения полученные результаты позволяют сказать, что изготовление труб и валов, подверженных высоким скоростям вращения, путем нанесения жесткого
(44)
покрытия или вставок позволит значительно увеличить ресурс данной детали. Также, как видно из представленных графиков пластических деформаций, мы можем заранее спрогнозировать, какие части данной детали будут подвержены значительным деформациям. В наших обозначениях это будут области I-IV. Границы максимального распространения этих областей мы можем приблизительно предсказать из результатов упругого решения, причем положение этих границ будет зависеть только от геометрических параметров детали, о чем нам говорит формула (44). Это позволит дополнительно увеличить ресурс детали, если подвергнуть данные области дополнительной обработке, в результате которой повысится сопротивление деформации данной области детали (повышение предела текучести).
СПИСОК ЛИТЕРАТУРЫ
1. Бегун А.С., Ковтанюк Л.В. Расчет напряжений, деформаций и перемещений во вращающемся диске в условиях ползучести // Вестник Чувашского гос. педагогич. ун-та им. И.Я. Яковлева. Серия: Механика предельного состояния. 2019. Т. 39, № 1. С. 84-93. DOI: https://dx.doi.Org/10.26293/chgpu.2019.39.1.011
2. Прокудин А.Н., Фирсов С.В. Вязкопластическое течение вращающегося полого цилиндра // Дальневосточный математический журнал. 2018. Т. 18, № 2. С. 242-260.
3. Gamer U., Lance R.H. Stress distribution in a rotating elastic-plastic tube. Acta Mechanica. 1983;50(1-2): 1-8. DOI: https://dx.doi.org/10.1007/BF01170437
4. Gamer U., Mack W., Varga I. Rotating elastic-plastic solid shaft with fixed ends. International Journal of Engineering Science. 1997;35(3):253-267.
DOI: https://dx.doi.org/10.1016/S0020-7225(96)00085-7.
5. Gamer U., Sayir M. Elastic-plastic stress distribution in a rotating solid shaft. Zeitschrift für Angewandte Mathematik und Physik ZAMP. 1984;35(5):601-617. DOI: https://dx.doi.org/-10.1007/BF00952107
6. Güven U. Elastic-Plastic Rotating Disk with Rigid Inclusion. Mechanics of Structures and Machines. 1999;27(1):117-128. DOI: https://dx.doi.org/10.1080/08905459908915691
7. Güven U., Parmaksizoglu C., Altay O. Elastic-Plastic Rotating Annular Disk with Rigid Casing. ZAMM. J. of Applied Mathematics and Mechanics. Zeitschrift für Angewandte Mathematik und Mechanik. 1999;79(7):499-503.
8. Hodge P.G., Balaban M. Elastic-plastic analysis of a rotating cylinder. Intern. J. of Mechanical Sciences. 1962;4(6):465-476. DOI: https://dx.doi.org/10.1016/S0020-7403(62)80008-3
9. Lindner T., Mack W. Residual stresses in an elastic-plastic solid shaft with fixed ends after previous rotation. ZAMM. Zeitschrift für Angewandte Mathematik und Mechanik. 1998;78(2):75-86. DOI: https://dx.doi.org/10.1002/(SICI)1521-4001(199802)78:2<75::AID-ZAMM75>3.0.CO;2-V
10. Mack W. Rotating Elastic-plastic tube with free ends. Intern. J. of Solids and Structures. 1991;27(11): 1461-1476. DOI: https://dx.doi.org/10.1016/0020-7683(91)90042-E
11. Mack W. The rotating elastic-plastic solid shaft with free ends. Tech. Mech. 1991;(12): 119-124.
FEFU: SCHOOL of ENGINEERING BULLETIN. 2019. N 4/41
Mechanics of Deformable Solids www.dvfu.ru/en/vestnikis
DOI: http://www.dx.doi.org/10.24866/2227-6858/2019-4-2 Prokudin A., Firsov S.
ALEXANDER PROKUDIN, Candidate of Engineering Sciences, Leading Researcher, e-mail: [email protected]
AuthorID: 415431, SPIN: 6812-2451, ORCID: 0000-0002-5156-424X,
ResearcherID: N-9344-2016, ScopusID: 35722777500
SERGEY FIRSOV, Junior Researcher, e-mail: [email protected]
AuthorID: 783315, SPIN: 8267-2329, ORCID: 0000-0001-7446-6231,
ResearcherID: D-1966-2018, ScopusID: 56976208300
Institute of machinery and metallurgy of Far Eastern Branch of RAS
1 Metallurgov St., Komsomolsk-on-Amur, Russia, 681005
Elastoplastic deformation of a rotating hollow cylinder with a hard coating on the inner and outer walls
Abstract: Consideration is being given to a hollow cylinder which is made of an ideal isotropic elastoplastic material in conditions of plain strain and with restriction on radial displacement on the inner and outer surfaces. It is assumed, that total strain in the cylinder is small and consists of elastic and plastic strains. Stresses are related with elastic strains by Hooke's law, plastic strains are determined by the Tresca conditions and the law of plastic flow associated with it. The deformation of the cylinder occurs due to centrifugal forces arising from the rotation of the cylinder with an angular speed, which first increases to the selected maximum value, and then decreases until the cylinder stops completely. It is assumed that the rotation speed slowly changes with time and angular acceleration is negligible.
It is established that plastic flow occurs simultaneously on both inner and outer surfaces of the cylinder. The rotation speed corresponding to the onset of plastic flow is found from the solution of the elasticity problem. Exact analytical solutions of the determining system of equations for each plastic field are found. The system of conditions at the boundaries between the regions providing continuity of the obtained solutions throughout the cylinder is given. The calculation results are illustrated by graphs of displacement, stresses and plastic deformations.
Keywords: rotating cylinder, plain strain, elastoplastic, plastic flow, Hooke's law, Tresca condition, secondary plastic flow, residual stress.
REFERENCES
1. 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;39(1):84-93. DOI: https://dx.doi.org/10.26293/chgpu.2019.39.L011
2. Prokudin A.N., Firsov S.V. Viscoplastic flow in a rotating hollow cylinder. Far Eastern Mathematical Journal. 2018;18(2):242-260.
For a complete list of references, see the previous page.