Механика
Вестник Нижегородского университета им. Н.И. Лоб аче вского, 2013, № 1(3), с. 120-129
УДК 539.374
МОДЕЛЬ УПРУГОПЛАСТИЧЕСКОГО ТЕЧЕНИЯ ПРИ ПЕРЕМЕННОЙ СКОРОСТИ ДЕФОРМИРОВАНИЯ
© 2013 г. Л.А. Игумнов1, Т.Д. Каримбаев2, Ш. Мамаев3
1НИИ механики Нижегородского госуниверситета им. Н.И. Лобачевского ^Центральный институт авиационного моторостроения, Москва ^Московский физико-технический институт
Поступила в редакцию 30.11.2012
Многочисленные экспериментальные исследования свидетельствуют о заметном влиянии скорости деформации на пределы текучести, прочности, законы упрочнения и другие параметры процессов деформирования в различных материалах. Аналитические подходы, учитывающие указанные эффекты скорости деформаций, представлены в большом числе работ и продолжают разрабатываться. Однако, инженерные подходы, позволяющие в рамках экспериментальных и аналитических достижений описать сложные явления при переменных скоростях деформации, недостаточно развиты. Модель упругопластического течения при переменной скорости деформирования, разработанный в предположении, что скорость деформации является параметром, может в некоторой мере восполнить этот пробел. Он построен на совокупности экспериментальных зависимостей напряжений от деформаций, полученных при различных значениях скорости деформации.
Ключевые слова: упругость, пластичность, вязкость, текучесть, скорость деформации.
Наиболее значимые эффекты скорости деформаций
Влияние скорости деформации на механические характеристики материалов было установлено испытаниями на удар еще в конце XIX и в начале XX веков [1]. Изучению свойств материалов при динамических нагрузках посвящено большое число экспериментальных работ. Обзор экспериментальных исследований с анализом результатов испытаний и способов их проведения выполнен в работах [1—6].
На основании изучения опубликованных экспериментальных исследований, выполненных при различных скоростях деформирования, установлено, что наиболее чувствительными параметрами материалов к их изменению являются пределы текучести, прочности, а также характеристики упрочнения. Результаты испытаний труб из мягкой стали, приведенные в [7], характеризуют влияние скорости деформаций на пределы текучести и прочности (рис. 1). Пределы текучести аа и прочности а* возрастают с увеличением скорости деформации £. Динамический предел текучести аа растет быстрее, чем предел прочности а*, и с некоторых значений скоростей £ он (а* ) практически
а
достигает значений предела прочности аь . Эффект межзеренного скольжения сменяется
эффектом разрыва зерен. Для формирования скольжения необходим большой промежуток времени, чего нет при высоких скоростях деформирования и материал проявляет псевдо-хрупкие свойства. В [8] установлено, что динамический предел текучести а* мягкой стали в области скоростей деформации £ порядка 102 ^103с- достигает величины, в 2.5^3 раза превышающей статический предел текучести а* . Зоной повышенной чувствительности предела текучести а* к изменению скоростей деформации £ является диапазон скоростей от 1 до 102 с-1.
Рис.1. Зависимость пределов текучести и прочности (штриховая линия) для мягкой стали
В опытах [9] с чистым алюминием изучалось влияние на кривые деформирования а ~ е и упрочнение материала ёа / ёе изменений скорости деформации £ в диапазоне от 103 до 104 с-1. Эксперименты [10, 11] показали, что упрочнение материалов ёа / ёе при фиксированном значении деформации е существенно зависит от скорости деформации £ На рис. 2 приведены диаграммы растяжения стали 45, полученные при различных скоростях нагружения 9.8-1012 и 9.8 • 106Н/(м2 • с) [12]. С повышением скорости нагружения предел текучести а^ возрастает почти вдвое. Также существенно возрастает касательный модуль «упругости» ёа / ёе характеризующий степень упрочнения материала. В [13] построены кривые деформирования а ~ е для алюминиевых сплавов А7, АД0 и АД1 в состоянии поставки и после отжига (рис. 3). Упрочнения ёа / ёе при деформациях е, превышающих 3%, на всех кривых практически одинаковы и не зависят от скорости деформации с.
Рис. 2. Кривые деформирования стали 45 при скорости нагружения 9.8 • 1012Н/(м2 • с) (верхняя кривая) и
9.8 • 106Н/ (м2 • с) (нижняя кривая)
[9, 15] это влияние обнаружено при изменении скорости деформации £ всего на два порядка. На рис. 4 воспроизведены результаты испытаний [15] по растяжению алюминиевой проволоки с комбинациями двух скоростей деформации: ^ = 5.5 • 103 с-1 (нижняя кривая) и
<^2 = 4.5 • 10"1 с-1 (верхняя кривая).
0 0.1 0.2 екгОДЮ
Рис. 3. Кривые деформирования сплава АД1
В экспериментах [14], а также в [9, 15] было установлено влияние истории изменения скорости деформации £ на кривую деформирования материала ёа / ёе . Существенно, что в работах
ю г%
Рис. 4. Экспериментальные зависимости при скоростях деформирования: ^=5^10'4 с'1 (1), перехода с кривой на кривую ^=103 с'1 (2), ^=3*103 с'1 (3), ^=5^10 с'1 (4)
В точке А при е = 3.6% скорость растяжения была быстро снижена с £2 до £1. Пунктирная линия показывает, что в последующем процессе деформирования при £ = £1 кривая а = а(е) стремится к кривой а(е), соответствующей скорости % = %, и практически совпадает с ней
через небольшой участок запаздывания Де =
0.8%. Аналогичный переход был выполнен на другом образце в точке В при большей предварительной деформации (е = 19,5%). В этом случае стремление к кривой а = а(е) при % = % происходит несколько медленнее. В точках С (е = 3.6%) и Б (е= 19.5%) совершен переход от меньшего значения скорости деформации £= £1 к большему - £ =£2. В этих случаях также наблюдался переход от кривой а(е), соответствующей скорости деформации £1, к кривой а = а(е), соответствующей £2. Эти испытания показывают, что кривые а = а(е) построенные в опытах с постоянной, но изменяющейся от эксперимента к эксперименту скоростью деформации £ могут играть роль опорных кривых <т ~ е при построении математической модели динамического деформирования материалов. Начальный этап деформирования при изменении скорости деформации £, как следует из испытаний, проведенных в [16], осуществляется упруго. Только затем наблюдается нелинейное приближение к соответствующей кривой деформирования а = а(%;% + Д%). Накопленная на первом этапе нагружения необратимая пластическая деформация ер полностью сохраняется при любых изменениях скорости деформа-
ции £. Однако, момент появления нелинейных деформаций при изменении скорости деформации £ недостаточно экспериментально изучен. Неопределенность предела пропорциональности в зависимости а = а(е) приводит к определенным трудностям при моделировании процесса перехода с кривой, соответствующей одной скорости деформации £, на другую.
Проведенный краткий анализ результатов экспериментальных исследований о распространении механических возмущений подтверждает зависимость сопротивления материала от скорости деформации и проявление эффекта вязкости, что отражается в различных реологических моделях поведения материалов под нагрузкой [17-19]. Разработанные математические модели, описывающие явления при высокоскоростных нагружениях, обычно учитывают влияние скорости деформации на характер нелинейного деформирования. Работы Л. Мал-верна [20], В.В. Соколовского [21], П. Пэжина [3], В.Н. Кукуджанова [22, 23], С. Калиского
[24], Р.М. Нагди и С.А. Мерч [25] и многих других внесли существенный вклад в развитие теории упруго'вязко'пластического деформирования. Использующиеся обобщенные на упруго-вязко'пластические среды модели деформирования предоставляют широкие возможности для исследований [3, 19-24]. Однако обобщения, как правило, в меньшей мере охватывают отдельные экспериментально наблюдаемые явления и зачастую приводят к заметным осложнениям при решении конкретных инженерных задач. Это обстоятельство, быть может, является одной из причин относительно узкого применения указанных моделей деформирования в практике проектирования изделий. В связи со сказанным, продолжаются попытки создания математических моделей неупругого поведения материалов, чувствительных к изменениям скорости деформации [25-27]. Далее описана модель упругопластического течения при переменной скорости деформирования, которая, являясь развитием работы [25], позволяет на основе экспериментально полученных для различных скоростей деформации с зависимостей сг ~ е предложить алгоритм решения конкретных задач.
2. Вариант модели пластически деформируемых материалов, чувствительных к скоростям деформации £
Экспериментальные данные показывают, что пластическое состояние при динамическом нагружении достигается при более высоких уровнях нагрузок, чем при статическом нагру-
/(а.. ,ер,%.., х) = 0.
жении. Уровни пластических деформаций могут заметно зависеть от истории нагружения, в том числе, и от истории изменения скоростей деформации % . Последние эффекты проявляются особенно четко для материалов с хорошо выраженными пределами текучести а3. Отдельные аспекты влияния скорости деформации % на состояние тела могут быть учтены, если поверхность нагружения [25] записать в виде
р (1)
Здесь а., ер, %. являются компонентами тену у у
зоров напряжений, пластических деформаций и скоростей деформации, а X,- - скалярный параметр упрочнения, зависящий от истории изменения пластических деформаций ер . Параметр
упрочнения х,- является постоянным при фиксированных значениях ер , и остаточные деформации определяются его величиной. В качестве
величины X обычно принимают параметр Од-квиста [28]:
1/2
(2)
определяющий уровень накопленной пластической деформации е за всю историю нагружения (от 0 до ?). Влияние параметров вязкости на упругое и пластическое деформирование можно связать с зависящими от них изменениями характеристик упругости, а также поверхности нагружения (1). Указанные эффекты вслед за [25] можно аналитически записать, используя известные подходы теории пластического течения (см. [29]). В соответствии с экспериментальными выводами о незначительном влиянии скоростей деформаций £ на характеристики упругости принимается линейная зависимость шаровых составляющих тензоров напряжений а = (Гуду и деформаций е = е^ду ,
е = а/К (в дифференциальной форме
ёе = ёа/К); (3)
е = е8..+ эа =а8 + 5..,
У У У ’ У У У ’
где К - объемный модуль упругости, 8 -
символ Кронекера. Далее предполагается возможность представления компонент девиатора деформации в Зу виде суперпозиции упругой зУ
и пластической составляющих зр :
(4)
и линейная зависимость упругих з составляющих девиатора з деформации от соответству-
ющих составляющих девиатора напряжений ^ : зр = 5 / 2/и (в дифференциальной
форме ёзр = dsij / 2и). (5)
Появление приращений пластических составляющих девиатора деформаций ёзр связывают (см. [29, 30]) с поверхностью / нагружения (1), полагая их ортогональность к поверхности нагружения:
Ж
(6)
с1эр = / d X.
дs..
У
Здесь параметр X размерности напряжения пока неизвестен. Изменения пластических составляющих ёзр девиатора деформаций связаны с изменениями поверхности нагружения (1):
— = — сЬ„ +—ёзр + — + — йх. (7)
7 дstJ У дзр У д%в % дх Х
Условия дальнейшего активного нагружения или разгрузки следуют из определения поверхности нагружения (1) и характера ее изменения (7). В процессе разгрузки и нейтрального нагружения поверхность ^0 не изменяется, нет пластических деформаций (ёзр = 0) и их
накоплений (ёх = 0). Условие разгрузки и нейтрального нагружения в соответствии с (7) можно записать в виде
— < 0 или — = —^ +— < 0 . (8)
&у ' д%У '
Рост пластических деформаций ёзр * 0 и
появление вклада ёх = 0 для упрочняющегося тела связаны с расширением поверхности нагружения
(Г = 0и<Г = ^Д^>0. (9)
Ассоциированный с поверхностью нагружения /=0 (1) закон пластического течения (6) с учетом влияния скорости деформации £ принимает вид:
аэр
о, а/ < 0 и dsт +-^- «£ < о,
дЧ ро д£ ^ ’
дЧрЧ дЪрщ
д/ ^ ^ д/ д/
-^~аХ, а/ = 0 и -^—ds Л—-— «Е дч.. дч р д£
II р° ^ р°
(10)
>0,
вающего поля девиаторов напряжений ^ , деформаций з и их скоростей £, является более
глубоким и сложным. Принятая поверхность нагружения (1) зависит от скорости деформации £у, её форма и размеры изменяются в зависимости от скорости деформации £. Скорость деформации хотя присутствует в (1) как независимый параметр, тем не менее, она дифференциально связана с деформацией еу. В соответствии с определением (2) и зависимостями (6) имеют место соотношения
ёх = ёзр = (0.5йзр йзр )05,
ах= [0.5(д/ / )(д/ / )]05 ах,
(11)
где
эр = (0.5эрэр )05, 5. = (0.55.5.. )0
1 V 1111 ' 5 I V 1111 '
(12)
Скалярный множитель йХ, определяющий линейную зависимость между компонентами приращения девиаторов пластических деформаций ёзр и напряжений £ , находится из равенства й/ = 0. Подставляя в соотношение (7) выражение приращения девиатора пластической деформации ёзр (6) и приращения накопленной пластической деформации (11),
можно получить коэффициент йХ:
ах = -Ф
(13)
У ~У
Здесь введен новый скалярный параметр
ф1 = //л/
дх
дэр д$11
0.5
д*и д*и
(14)
Если использовать выражение (6), то приращение девиатора пластической деформации ёзр на активном участке нагружения можно
представить в виде
с1эр = -/ Ф
дч,.,.
д/
дч'
д/ л,
дїтп
(15)
где йХ - по-прежнему неизвестный параметр. По форме соотношения (10) не отличаются от традиционного ассоциированного закона пластического течения. Однако функция нагружения (1) является более сложной, чем в теории пластичности [28]. Фактическое содержание
выражения (—) ёХ, одновременно увязы-
Из выражений (14) и (15) следует, что изменения девиатора пластической деформации азр
могут происходить только на активном участке деформирования и они осуществляются как за счет изменения компонентов девиатора напряжений Ж у и изменения скоростей деформаций
, так и в силу изменения поверхности
нагружения/=0 от этих параметров. Вводя обозначения
= -Ф 111, , = -Ф1 (д//д%т), (16) выражения (15) для приращений пластических
деформаций dap можем записать в виде
dap = (dfl-s.. )(Fsds +F? dL ). (17)
ij \ J / ij mn mn mn ^mn / V /
В выражении (17) составляющая Kn (spq,L,z) приращения пластических деформаций dap обусловлена изменениями девиатора напряжения smn при постоянных значениях компонент тензора скоростей деформации
Lmn , а составляЮЩая FL (spq ,4pq , Z) вЫЗвана
изменениями тензора скоростей деформаций Lmn при постоянных значениях компонент девиатора напряжения smB. Все эти изменения рассматриваются при постоянной величине накопленной пластической деформации Z = sp = const. Выражение (17) для приращений компонентов девиатора пластических деформаций dap по форме подобно соотношениям неизотермической теории пластичности [31, 32]. Однако в отличие от зависимостей неизотермической теории пластичности, здесь переменные компоненты девиаторов напряжений
3, деформаций э и компоненты скоростей
деформации £j являются взаимно связанными величинами. В несвязанной неизотермической теории пластичности температура t является не зависящей от sff и э переменной. Влияние на
поверхность нагружения f таких параметров, как переменная скорость деформации £ij в соотношении (17) и температура t в несвязанной неизотермической теории пластичности определяется внутренним их качественным отличием. Несмотря на эти отличия, в [25] удалось показать нормальность вектора мгновенных приращений dap компонент девиатора пластических деформаций поверхности нагружения f=0 в форме (1).
Функция нагружения f=0 (1) записана в достаточно общем виде, что является причиной неопределенности модели. Кроме того, в совокупности принятых условий не выяснено каким образом меняются форма и размеры поверхности нагружения f=0 при изменении этих условий. Построение определяющих соотношений упруго-вязко-пластического деформирования исследуемой среды и разработки соотношений, пригодных для проведения непосредственных вычислений, завершается установлением конкретного вида функций f а также F^n и F^n, являющихся следствием принятой функциональной формы поверхности нагружения f.
Кроме того, необходимо установить способы экспериментального определения введенных функций. С целью конкретизации формы поверхности нагружения f=0 ниже используется один из частных случаев, соответствующий обобщенному на вязкое пластическое деформирование условию текучести Мизеса
f(siz ,L) = s2-c2(ZL) = 0. (18)
Здесь Sj, £ - интенсивности касательных напряжений и скоростей деформации, которые определяются соотношениями (11) и
[2 Ь -е -е..
Чз jj Чз е -е-. (19)
В (18) функция <Js(z Li) является мгновенным пределом текучести, равным значению интенсивности девиатора напряжения s\ при определенной интенсивности скорости деформации L, и накопленной пластической деформации Zi в текущий момент времени t. Начальная поверхность нагружения cs(L) = cs(0,L)
совпадает с обычным пределом текучести, зависящим от скорости деформирования £ (рис. 5). Последнее обстоятельство непосредственно охватывает экспериментальные результаты зависимости предела текучести материалов от скорости деформирования £. Геометрически зависимость (18) изображает в пространстве (si; эь £) поверхность деформирования образца (см. рис. 5) с отдельными кривыми деформирования, соответствующими различным скоростям деформации £. Процесс деформирования образца при заданной переменной скорости деформации £ = const можно изобразить в виде кривой s = g (э ) на этой поверхности. Совокупность
кривыгс s = g (э ) до предела текучести crs(0, £) образует линейчатую поверхность упругого деформирования. В сечениях э = const изображается зависимость интенсивностей напряжений Si, в том числе, пределов текучести c S от интенсивности скорости деформаций £i. Сечения s = const определяют зависимость интенсивности сдвиговых деформаций ai от скорости деформаций £■. Следует заметить, что модель построена в девиаторном пространстве напряжений и деформаций (si; эг), а используемый параметр £ в пространстве начальный деформаций. Обоснованием этому является то, что функция f=0 (18) при заданном £ не зависит [33] от гидростатического давления с
Всякая точка, изображающая пластически деформированное состояние материала, в любой
момент времени должна находиться на поверхности 5 =а (%,х), которая зависит от накопленной необратимой пластической деформации х и текущего значения интенсивности скорости деформации % (вязкость). С учетом принятой формы поверхности нагружения (18) основное соотношение (15) при продолжающемся нагружении можно записать так:
ёзр = (5у )[<*,■- (да/ д%)й% ]/[5у (да/ дх)] .(20)
В (20) были учтены равенства, вытекающие из (18):
д/ / д5УУ = 5У; д/ / д% = -2а!1да!1 / д% =
= -25 да5 / д%; = 25уА , (21)
д/ / дх = -— — / дх = -25 Уда5 / дх.
Соотношения (21) позволяют записать
1 ( 4 0,5
Ф =-(2ада / дх) {0.55.5.. } =
V 5 5 ( У У) (22)
= -252(да / дх).
Если ввести обозначения Г5 = 1/[(да /дх)], = Г5(даs /д%)], (23)
то ассоциированный с обобщенной поверхностью Мизеса (18) закон течения (20) с учетом влияния скорости деформации % можно записать в виде
ёзрУ = (У*, )[^й5У - ]. (24)
Соотношения (24) имеют место только на активном участке нагружения, когда в соответствии с (10) и (21) выполняется условие
^ - (да5 / д% )ё% > 0 , (25)
в противном случае ёзр = 0 . Соотношения (24)
дают возможность учесть известные экспериментальные факты [1]. Из выражения (24), (25) видно, что для предложенного соотношения условие активного нагружения и разгрузки, отличаются от существующих подходов. Способ получения соотношений (24) принципиально не отличается от того, что было использовано в
[25]. Точно так же, как и в аналогичных соотношениях работы [25], зависимости (24) продолжают содержать неопределенности, которые не позволяют их использовать в конкретных расчетах. Указанные неопределенности связаны с нерегламентированным (общим) характером изменения поверхности текучести а (х, % ) от
параметров х и %. Представляется, что конкретизация способов получения введенных соотношениями (23) функций & и выполнено впервые.
В соответствии с основными предположениями и соотношениями (2)-(5) и (24) на активном участке деформирования имеют место определяющие соотношения в форме
ёе = ёа8 / К + йъ.. / 2и +
У 8 У (26)
+5 / 5, )(Г5ё5, - ).
Далее предполагается, что поверхность пределов текучести Si = а5 (£ х) не зависит от вида напряженного состояния и её можно построить на основе результатов испытаний на одноосное растяжение, выполненных при различных, но постоянных скоростях деформации £=сош1. В
виртуальном испытании образца на растяжение начальное состояние определяется достигнутым уровнем растягивающих напряжений а11 и скоростей деформаций £:
а = аи/3; а = 0 (при i и у не равном 1 одновременно),
511 2стп/3, S22 S33 a11 /3,
512 = S13 = S23 = 0 Si ^пА/3,
S11 /Si = 2/V3, S22 / Si = S33 / Si = -1 / V3 S12/ Si = S13/ Si = S23/ Si = 0.
(27)
Если достигнутое равновесное состояние (27) приведено к новому равновесному состоянию с помощью дополнительного растягивающего напряжения dc ф 0, то
dc = dan /3, dsjj = 2dcn /3, ds22 = ds33 = -dc /3,
dst = dc / V3 dL = 0,
и в соответствии с (26), (27) приращения компонент тензора полной деформации выгчисля-ются из соотношений
dsjj = [K +1/ л + 2Fs ]ds, />/3,
dsl2 = dsl3 = ds23 = 0, (28)
ds22 = ds33 = [K - ds22 / 2л - Fs]ds( / \/3.
При этом приращение интенсивности сдвиговых деформаций находится из
dat = (1 / 2 л + Fs )ds . (29)
Из (29) можно найти одну из искомыгс функций
Fs = 1 /(ds. /da.)- 1/2л =
, (30)
= 1 / 2л -1 / 2л; 2л = d^( / da.
Здесь 2ju' — касательный модуль сдвига к экспериментальной кривой \ ~э . Таким образом, для экспериментального определения функции F следует рассмотреть простое растяжение образца при<^ = const. Из выфажения
(30) и условия /л' <л видно, что функция F всегда является положительной и обусловливает возрастание пластической деформации э? >- 0 при c/s) >- 0 и отсутствие её приращения (r/j/' = 0) при dst -< 0. Если достигнутое
начальное равновесное состояние (27) приведено к новому равновесному состоянию путем изменения интенсивности скорости деформации dL ф 0 , то приращения напряжений равны нулю dc = 0, dc = dsjj = ds22 = ds33 = ds = 0, и в соответствии с (26) приращения компонент тен-
зора деформации и интенсивности сдвиговых деформаций вычисляются из соотношений
deu =-(2/S)Fndn;
11 ' l . (31)
ds22 = ds33 = (1/v3)Fndn; dэiFndni.
В соответствии с (31) вторая неизвестная функция Fn определяется из выражения
Fn = dэ. / dn. (32)
По экспериментальным кривым деформирования s\ ~ э(, построенным при различных значениях скорости деформации П , должны быть построены сечения j ~ поверхности текучести при постоянном уровне интенсивности напряжений d = const. Касательные к кривой j ~ , позволяют определить эксперименталь-
ные выражения функции F (см.(31)) при текущем значении П . В соответствии с экспериментальными данными д^/д П >0 и соотношением (24) функция F принимает только отрицательные значения. Поэтому согласно равенству (24), рост пластических деформаций (cbf >- 0) можно наблюдать только при уменьшении скорости деформации d£ < 0. Этот эффект можно видеть и из соотношения (31) для основной компоненты деформации dsn. Таким образом, при возрастании скорости деформирования £ приращение пластических деформаций не наблюдается. Эти результаты качественно согласуются с экспериментальными наблюдениями. Подставляя выражения (30), (32) в соотношения (26), можно получить закон пластического течения с учетом влияния скорости деформации dє,.,. = daS„ / K + ds„ / 2u +
J J J +d / Si)H(FSdst - Fndn) = (33)
= dcrSy / K + dstj / 2u + (Sy / s )G.
Здесь
G = H(FSdsi - Fndn);
[0, если x < 0, (34)
[1, если x > 0.
Физические соотношения (33) могут быть записаны в виде
dєп = [dc - v(da22 + da33)]/ E + (sn / st )G, dєп = dsn / 2u + (s12 / d )G, dє22 = [dc - v(dc + da33)] / E + (s22 / st )G.
dє13 = dsl3 / 2u + (s13 / d )G, dє33 = [da - v(da22 + dan)] / E + (s33 / st )G. dє23 = ds23 / 2u + (s23 / d )G,
H (x) =
в котором V - коэффициент Пуассона, Е - модуль упругости.
При решении конкретных задач в каждый момент времени в каждой точке необходимо уточнить условия нагружения и разгрузки, которые связаны с переменными значениями скоростей деформации и уровнями интенсивностей напряжений 5 . Условия нагружения и разгрузки при простом растяжении и сжатии отличаются от исследованного в [31] варианта теории неизотермического течения. Ниже приводятся различные случаи нагружения и разгрузки. Для многих материалов с увеличением скорости деформации £ кривые деформирования повышаются [13, 34, 35]. Вследствие этого и соотношений (20), (24), (25) при >- 0 и
< 0 обязательно происходит нагружение (рис. 6а), так как согласно (25)
- (5с / 5^ > 0 . При < 0 и > 0
наблюдается процесс разгрузки (рис. 6б) и
< (5^ / 5£')й£,. На этих рисунках точка А обозначает начало рассматриваемого процесса нагружения, а точка В соответствует состоянию материала в его конце. Стрелка указывает направление протекания процесса для следующего этапа нагружения, соответствующего возрастанию напряжений.
Рис. 6. Зависимость кривых деформирования от скорости деформации: а - процесс нагружения (при йс > 0, < 0);
б - процесс разгрузки (при йс < 0, > 0)
Более сложные варианты нагружения показаны на рис. 7:
1. > 0, > 0, > (5с /
(нагрузка) (рис. 7а);
2. > 0, > 0, < (5с / 5^(раз-
грузка) (рис. 7б);
3. < 0, < 0, > (5с / 5^)й^
(нагрузка) (рис. 7в);
4. < 0, < 0, < (5с / 5£)й£ (раз-
грузка) (рис. 7г).
Рис. 7. Сложные варианты нагружения
На рисунке ds\, dЛ - соответственно изменения напряжений и скоростей деформации, (да I дП )d^i - изменение пределов текучести. Из условия (21) вытекает, что можно выбрать такие изменения скоростей напряжений ds и
деформаций dэf, при которых процесс при переменной скорости деформирования в определенном диапазоне будет все время упругим или, наоборот, непрерывно сопровождается накоплением пластических деформаций.
Заключение
Разработан вариант модифицированной теории течения, связывающей дифференциал пластической деформации с дифференциалами интенсивностей напряжений и скоростей деформации. Введены функции Fs (st, х ,Л),
F?(st ,X ,4i), определяющие вклад изменений напряжений и скоростей деформации в аккумулируемую пластическую деформацию. Функции Fs (s,., Хі Лі ), Fi(si, Хі ,4i) определяются из экспериментальных кривых деформирования, построенных при различных скоростях деформации. В рамках одноосного нагружения проанализированы условия активного нагружения и разгрузки при изменении основных параметров состояния - напряжения, скорости деформации. Такие экспериментальные факты, как повышение предела текучести, упрочнения, влияние истории изменения скорости деформации, плато деформаций, охватываются предложенной моделью упруго-вязко-пластического деформирования.
Список литературы
1. Бэлл Дж. Экспериментальные основы механики деформируемых твердых тел. Часть II. Конечные деформации. М.: Мир, 1984. 431 с.
2. Васин Р.А., Ленский В.С., Ленский Э.В. Динамические зависимости между напряжениями и деформациями II Механика. Новое в зарубежной науке, 5. Проблемы динамики упруго-пластических сред. M.: Мир, 1975. С. 1-38.
3. Пэжина П. Основные вопросы вязкопластич-ности. M.: Мир, 19б8. 175 c.
4. Кэмпбелл Дж. Эксперименты при высоких скоростях деформации II Механика: Сб. перев. и обзоров. Иностр. период. лит. 19бб. № 3.
З. Хольцер А.Дж. Обзор экспериментальных исследований в области динамической пластичности II Теоретические основы инженерных расчетов: Труды Амер. об-ва инж.-механиков. 1979. Т. 101, № 3. C. Зб-б1.
б. Каримбаев Т.Д., Мамаев Ш. Параметры
материалов, чувствительные к скоростям деформации // Новые технологические процессы и надежность ГТД: Научно-технический сборник
статей. М.: ЦИАМ, 2008. Вып. 8. С. 7-37.
7. Clark D.S., Duwez P.E. The influence of strain rate on some tensile properties of steeln // Proc. Amer. Soc. testing Materials. 1950. V. 50. P. 560-575.
8. Manjoine M. Influence of rate of strain and temperature on yield stresses of mild steel // J. Appl. Mech. 1944. V. 11. P. 211-218.
9. Васильев Л.И., Былина А.С., Загребеннико-ва М.П. О влиянии перемены скорости деформирования на пластическое растяжение // ДАН СССР. 1953. Т. 90. №5. C. 767-769.
10. Schweinghofer A., March W. Phys. stat Solids. 1968. V. 25. 2.
11. Lindholm U.S., Yeakley L.M. Dynamic deformation of single polyerystaline aluminium // J. Mech. and Phys. Solids. 1965. V. 13, 41. P. 41-53.
12. Орленко Л.П. Поведение материалов при интенсивных динамических нагрузках. М.: Машиностроение, 1964.
13. Брагов А.М., Ломунов А.К., Шахалов Г.И. Структура и механические свойства алюминия при высокоскоростной деформации // Прикладные проблемы прочности и пластичности. Методы решения / Горьков. университет. Горький, 1988. C. 129-134.
14. Tietz T.E., Dorn J.E. // Trans. Amer. Soc. Metals. 1949. 41A, №2.
15. Васильев Л.И., Еремина Л.И. О некоторых особенностях пластического растяжения с переменной скоростью // ДАН СССР. 1953. 93, №6. C. 10191023.
16. Малышев Б.М. Распространение догрузочных импульсов по натянутой проволoке // Изв. АН СССР. Механика и машиностроение. 1960. № 2. С. 120-124.
17. Клепачко Я. Анализ распространения фронта упругопластической волны в металлическом стержне // Механика: Сб. перев. и обзоров. Иностр. период. лит. 1971. №6. С. 112-130.
18. Шемякин Е.И. Динамические задачи теории
упругости и пластичности. Новосибирск: Изд-во
Новосибирского ун-тa, 1968. 336 с.
19. Cristescu N. Dynamic plasticity. Amsterdam: Horth-Holland Publ. Company, 1967. P. 608.
20. Малверн Л. Распространение продольных пластических волн с учетом влияния скорости деформации // Механика: Сб. перев. и обзоров. Иностр. период. лит. 1952. Т. 11, № 1. С. 153-161.
21. Соколовский В.В. Распределение упруго-вязко-пластических волн в стержнях // ПММ. 1948. Т. 12, № 3. С. 261-280.
22. Кукуджанов В.Н. Одномерные задачи распространения волн напряжении в стержнях. М.: ВЦ АН СССР, 1977. 55 с.
23. Кукуджанов В.Н. Распространение волн в упруговязкопластических материалах с диаграммой общего вида // Изв. РАН. МТТ. 2001. № 5. С. 96-111.
24. Kaliski S. On certain equations of dynamics of an elastic/viscoplastic body. The strain hardening properties
and the influence of strain rate // Bull. Acad. Polon. Sci., Serie Sci. Techn., 11, 7(1963).
25. Нагди Р.М., Мерч С.А. О механическом поведении вязко-упруго-пластических тел // Прикл. механика: Тр. Америк. об-ва инж.-механ. Сер.Е. 1963. T. 30, № 3. C. 3-12.
26. Эллин Ф., Ся З. Конститутивная модель неупругого поведения материала зависящего от скорости // Упругопластическое течение. Современное машиностроение. Сер. Б, 1991. № 8. C. 129-140.
27. Каримбаев Т.Д., Мамаев Ш. К исследованию поведения упруго-вязко-пластических тел при динамических нагрузках переменной скорости. В.кн: Актуальные проблемы механики деформируемого тела. Алма-Ата: Гылым, 1992. Ч. 2. C 53-68.
28. Ивлев Д.Д., Быковцев Г.И. Теория упрочняющегося пластического тела. М.: Наука, 1971. 231 с.
29. Качанов Л.М. Основы теория пластичности.
М.: Наука, 19б9. 420 с.
30. Новожилов В.В. Теория упругости. Л.: Судпромгиз, 1958. 370 с.
31. Биргер И.А., Демьянушко И.В. Теория пластичности при неизотермическом нагружении II Изв. АН СССР. МТТ. 19б8. № б. С. 10-11.
32. Шевченко Ю.Н. Термопластичность при переменных нагружениях. Киев: Наук. думка, 1918. 283 с.
33. Клюшников В.Д. Математическая теория пластичности. М.: Изд-во МГУ, 1979. 207 с.
34. Lindholm U.S. Some experiment with the split Hopkinson pressure bar // J. Mech. and Phys. Solids. 19б4. Vol. 12. №5. P. 311-333.
33. Викторов В.В., Шапиро Г.С. Об определении динамических диаграмм растяжения металлов при умеренно-высоких скоростях деформаций II Изв. АН СССР. МТТ. 19б8. №2. C. 184-181.
THE MODEL OF ELASTIC-PLASTIC FLOW AT A STRAIN VARIABLE RATE L.A. Igumnov, T.D. Karimbayev, Sh. Mamayev
Numerous experimental studies show a marked influence of strain rate on the yield stress, strength, hardening laws and other parameters of the deformation processes in various materials. Analytical approaches that take into account these effects of strain rate are presented in many researches, and continue to be developed. However, engineering approaches to within experimental and analytical advances describe complex phenomena at strain varying rates, are underdeveloped. Model of elastoplastic flow at strain variable speed, developed under the assumption that the deformation rate is a parameter, can fill this gap. It is built on a set of experimental stress-strain curves obtained at different strain rates.
Keywords: elasticity, plasticity, viscosity, fluidity, the strain rate.