МАШИНОСТРОЕНИЕ МАТЕРИАЛОВЕДЕНИЕ И ТЕРМИЧЕСКАЯ ОБРАБОТКА МЕТАЛЛОВ
УДК 539.3
В. М. ГРЕШНОВ, И. В. ПЯТАЕВА, В. Е. СИДОРОВ
ФИЗИКО-МАТЕМАТИЧЕСКАЯ ТЕОРИЯ ПЛАСТИЧНОСТИ И ПОЛЗУЧЕСТИ МЕТАЛЛОВ
Развивается теория необратимых деформаций металлов, основанная на ранее предложенной концепции синтеза основных положений и методов математической теории пластичности и физики прочности и пластичности. С единых позиций получены скалярные и тензорные определяющие соотношения, описывающие пластическую деформацию и ползучесть металлов в широком диапазоне температур, скоростей деформаций и напряжений в условиях сложного нагружения и с учетом истории нагружения. Пластичность; ползучесть; физико-феноменологическийподход; определяющие соотношения; физико-математическая теория
ВВЕДЕНИЕ
В математических теориях пластичности и ползучести металлов не полностью решены вопросы описания процессов деформирования в условиях сложного нагружения и с учетом истории нагружения [1,2]. И это, видимо, принципиальные трудности теории, связанные с чисто феноменологическим методом описания неравновесных процессов деформации (без учета структуры материала и ее эволюции). Поэтому расчетные методы, базирующиеся на данных теориях, являются, в отмеченных случаях, весьма приближенными и имеют известные ограничения. Кроме того, базовые экспериментальные методы определения входящих в определяющие уравнения материальных констант и функций трудоемки и сложны.
Со времени пионерской работы [3] предпринимаются попытки дальнейшего развития теории на основе учета современных положений физики пластической деформации и ползучести металлов [4-6]. Однако к настоящему времени ни одна из предложенных концепций не получила развития до создания расчетных методов, обеспечивающих решение прикладных задач.
В работах [7-9] предложен подход к развитию теории необратимых деформаций, основанный на синтезе методов и основных положений механики и физики пластической деформации. С использованием уравнений кинетического баланса применительно к плотностям дислокаций выводится скалярное (одноосное напряженное состояние) уравнение пластического течения, которое представляет
собой записанный в дифференциальной форме функционал напряжения течения от характеристик структуры, макропараметров деформации и их функций от времени. Полученное скалярное уравнение, в отличие от феноменологического функционала сопротивления металлов пластической деформации [10], раскрывает микромеханизмы основных процессов, протекающих в металлах при различных термомеханических условиях деформирования. Поэтому оно облегчает выбор функции нагружения, формулировку соответствующего принципа максимума и получение тензорного уравнения пластического течения с использованием основного метода математической теории пластичности.
В данной работе обобщаются и уточняются полученные ранее теоретические результаты и, в рамках разрабатываемого подхода в теории пластичности, выводятся уравнения дислокационной ползучести металлов.
1. ОСНОВНЫЕ ПОЛОЖЕНИЯ И УРАВНЕНИЯ ФИЗИКО-МАТЕМАТИЧЕСКОЙ ТЕОРИИ ПЛАСТИЧНОСТИ
В основу модели пластической деформации положены следующие положения физики прочности и пластичности:
1) доминирующим механизмом деформации является дислокационное скольжение в зернах (это положение, в том числе, определяет в основном и область применимости модели).
2) сопротивление деформации обусловлено торможением потока дислокаций барьерами различной природы, на которых они оста-
навливаются. Основными барьерами являются неподвижные дислокации (лес дислокаций). Причиной упрочнения является повышение плотности неподвижных дислокаций.
3) разупрочнение связано с уменьшением плотности неподвижных дислокаций за счет преодоления некоторыми из них потенциальных барьеров, а так же за счет аннигиляции, динамических процессов полигонизации и рекристаллизации, образования микронарушений сплошности (деформационной по-врежденности). При этом необходимую энергию они получают за счет термической активации и работы действующих напряжений. Микромеханизмы перечисленных процессов, в основном, имеют диффузионную природу, следовательно, высота барьера равна энергии активации самодиффузии.
4) при низких температурах (холодная деформация) дислокации преодолевают барьеры силовым способом, то есть только за счет работы действующих напряжений.
Система уравнений модели, при допущении достаточной производительности источников дислокаций, для жесткопластической среды включает [7]:
1) кинетическое уравнение баланса скалярных плотностей подвижных и неподвижных дислокаций
dfh
dt
— Рд^дя Ря^яд:
(1.1)
где и рд — усредненные по объему плотности неподвижных (сидящих на барьерах) и подвижных (вносящих вклад в деформацию) дислокаций соответственно; ид8 — частота превращения подвижных дислокаций в неподвижные; — суммарная частота исчезновения подвижных дислокаций за счет преодоления барьеров и превращения снова в подвижные, аннигиляции при переползании, протекания процессов динамической полиго-низации и рекристаллизации, образования зародышевых микротрещин (деформационной поврежденности).
При этом для термически активируемых необратимых процессов
активационный объем; т — фактор Тейлора для поликристаллов (при хаотической разо-риентировке зерен для разных типов кристаллических решеток т = 2,8 -г 3,1).
2) физические уравнения связи макрохарактеристик пластической деформации с характеристиками дислокационной структуры
Є = Єоexp
Є = Pgblh
U - crV/m кТ
(1.3)
(1.4)
где — усредненный по системам скольжения модуль вектора Бюргерса дислокаций (для металлов ; — средняя скорость
скольжения подвижных дислокаций.
Так как три из четырех перечисленных выше микромеханизмов исчезновения неподвижных дислокаций имеют диффузионную природу, а энергия активации дислокационного зарождения микротрещин близка к ней (энергия сублимации) [11], то и в (1.2) и (1.4) считается одинаковой и равной энергии активации самодиффузии и = /ЗОЬ3, где ¡3 = для разных металлов, — модуль сдвига. В качестве объема активации взят объем одноатомной цепочки длиной 1 / (сегмент дислокации между двумя барьерами), то есть V = I)2/т/р^. Частота щ равна частоте тепловых колебаний этого сегмента щ = уф^[р1, где VI) = 1012с-1 — частота Дебая. Установлено, что ¿о = е^Ъ^/р^, где е* = г'о [8].
Решением системы уравнений (1.1-1.4) является функционал напряжения течения
(т = /3(T)mG (Т) Ь [1-
kT{t) ^ ¿*by/fa'
/3 G(T)fe3 ¿(і)
¿(*) n /2„ bv
~r,-----Ps VDbX
bX
x exp
' p (T) G (T) b3 -kT(t)
1
/2
dt
(1.5)
Vsg = Щ CXP
U - crV/m kT
(1.2)
где — предэкспоненциальный множитель;
— усредненная по вышеперечисленным процессам энергия активации; — интенсивность напряжений; — постоянная Больцмана; — термодинамическая температура; —
где A = v/vgs — средняя длина свободного пробега подвижных дислокаций.
В случае е, Т = const из (1.5) следует
а = fi(T)mG(T)b х
кТ
/З (Т) G (Т) Ь3
In
£*Ьу/р1
VpI, (1.6)
где — накопленная за время деформирования или при степени деформации плотность неподвижных дислокаций.
Для учета истории нагружения, обусловленной зависимостями и , и численного интегрирования запишем (1.5) в виде следующего оператора. Из (1.6) получаем
<тГа)=Р(Т)тв(Т)Ьх
с1а{д} =
/3 (Т) тО (Т) Ь ткТд
^у/Р^д) ^2\/Ря(д)
¿*Ьу/Ря(д)
1 + 1п
Ч.д)
^Ра{д) ■
(1.7)
Из (1.1) имеем
д'Ря{д)
1
ЬХ
Р*(д-Ц
£(д)
/
х схр
¡3 (Т) С (Т) Ь3
а(д- Р
кЪ
(д)
ск
где ;
мер расчетного шага при численном расчете диаграммы , на котором интенсивность деформации получает приращение .
Уравнения (1.7) и (1.8) дополним соотношениями:
= 1.2.
,п
(1.8)
но-
Р я{д) Ря{д-1) + ^'Р.ч(д)1
Чд)
= е
(9-1)
ск
(д)'-
Чд)=Л£{д)/М{д): с^{д) = с^^{д)/'и{д)'-
Щд) = Н>
(д-1)
ни.
г(д)
(при испытании образцов осадкой);
а(.д)
~ °{д)
с!а.
(1.9)
1.10)
1.11)
1.12)
1.13)
1.14)
1.15)
\д)
кТ,
(д)
¡3 (Т) С (Т) Ь3
1п
£*Ьу/Р.ч(д-1)
£(д)
\/ Ря{д-1); (1.16)
где
1.д)
начальное напряжение течения на
расчетном шаге ; —
изменение размера образца при испытании или приращение перемещения рабочего органа деформирующей машины; — текущий размер образца при испытании; — скорость деформирования на шаге (скорость перемещения рабочего органа деформирующей машины).
Уравнения (1.11) и (1.15) записаны при допущении соосности тензоров , и , из-за малости и .
Исходное напряжение течения материала (начальный «предел текучести»)1 при д = 1 определяется из (1.6), как
(гт = /3 (Т) тС (Т) Ь х кТ{1}
1
1п-
Ряо
/3 (Т)С(Т)Ъ1 е(1)
\/~РйО',
(1.17)
где — исходная в материале (до нагрева и деформации) плотность неподвижных дислокаций.
При экспериментальной проверке модели установлено, что коэффициент ¡3 с ростом температуры увеличивается в соответствии с эмпирическим уравнением х Т + 0,3639. Это объясняется доминированием объемной диффузии (по телу зерен) при высоких и трубочной (по ядрам дислокаций) при низких. Последняя, как известно [11], имеет меньшую энергию активации.
Из изложенной модели как частные случаи вытекают модели холодной ,
теплой и горячей
деформации, где — температура плавления в . Установлено [8], что при холодной деформации дислокации преодолевают барьеры силовым способом, то есть энергия активации преодоления барьеров целиком обеспечивается работой действующих напряжений. Условия холодной деформации:
/3(7Ь3 = стЬ2 / т%ДГн\
(1.18)
1 Кавычки означают, что у вязкопластичного материала отсутствует предел текучести в смысле его общепринятого определения [9]. При уменьшении а материал не переходит в упругое состояние, деформация прирастает при меньшей £.
¿о = ¿*Ьу/р^ = є = і/0 = vvbyffhi = vsg.
(1.19)
С учетом (1.18) и (1.19) из (1.6) и (1.8) для холодной деформации имеем
а = fimGbyffTs,
dp.4 = ( - Ps ] de.
(1.20)
(1.21)
Интегрируя (1.21) при начальных условиях є = 0, ps = pso получаем
Ps =
(ЬХ) 1 [схр (е) - 1] + Pso схр (є)
а = ftmGb ■
(ЬХ) 1 [схр (є) - Г] + pso
схр (є)
(1.22)
(1.23)
Начальный предел текучести для жесткопластического тела из (1.23) при е = 0 или из (1.17) с учетом (1.19) определяется как
ат = pmGbyfp^o-
(1.24)
Теплая деформация протекает в условиях переменного активационного объема V = = b2/^p~s, который уменьшается с увеличением е и ps и эволюционное уравнение соответствует (1.8). Горячая деформация протекает в условиях V = Ь2/yfp~s = const и уравнение (1.8) в этом случае имеет вид
d'Ps(g)
х схр
1 Р.ч(д-І) , ,----------
ьх “
Є(9)
ft [Т)С{Т)ЪЛ
а(я- Рь~
ПК
кТ,
(9)
de
(д)-
(1.25)
Для расчета зависимостей а (є) по изложенной модели необходимо определение двух параметров и , которые являются характеристиками исходной дислокационной структуры металлов. Это можно сделать известными металлографическими методами. Возможно также их определение на основе экспериментально полученной при стандартных испытаниях цилиндрических образцов осадкой в условиях холодной деформации диаграммы деформирования по формулам [7,8]:
Р.чо = (сттР) / (ftmGb)
А =
ехр
Ь (ftmG)2 [схр (в) - 1] гг2 схр (є) — (ftmGb)2 pso
(1.26)
(1.27)
где — экспериментальный предел текучести материала при холодной деформации; и — интенсивность деформации из интервала и соответствующее ей на экспериментальной диаграмме деформирования значение интенсивности напряжений. Эти формулы получены из (1.24) и (1.23).
Из модели следует утверждение: физикофеноменологический функционал вида
при известной
диаграмме деформирования материа-
ла, полученной в условиях холодной деформации, позволяет определять его напряжение течения в температурно-скоростном диапазоне, в котором доминирующим механизмом деформации является дислокационное скольжение, а основным механизмом упрочнения — блокировка подвижных дислокаций барьерами дислокационного типа (дислокации леса, границы зерен и субзерен).
С целью экспериментальной проверки этого утверждения и одноосной модели пластической деформации были получены экспериментальные диаграммы деформирования качественных сталей 10, 20, 38хгНм, 20Г2Р, 30Г1Р в широком диапазоне температур и скоростей деформаций. Диаграммы строились по результатам стандартных испытаний образцов на осадку по методике, изложенной в [8].
На рис. 1 в качестве примера приведены экспериментальные (точки) и теоретические (сплошные кривые) диаграммы деформирования стали 10: 1 (ф) — скорость деформации 10_2 с-1; 2 (•) — 10_3 с-1; 3 (■) — 1(Т2 с"1; 4 (А) — 10"1 с"1. Значения 0{Т) брали из [12]. Характеристики стали: р30 = ;; ;; ;
. Отрезки на точках обозначают симметричное десятипроцентное отклонение от экспериментально определенной величины.
Проверка показала удовлетворительное совпадение экспериментальных и теоретических кривых .
а
690
590
490
190
90
а
Т=293 °К (20°С) '
( + ^ 5 . 1
< ' Т=723 4 і °К (450°С) ^ * і і 4
і А ІІ * ■ А 1—• 1 !3 2
і
1 Р 1
0 0,1 0,2 0,3 0,4 0,5 0,6
а
Рис.1
Для решения прикладных задач теории пластичности, например, при математическом моделировании технологических операций и процессов обработки металлов давлением, необходимы физико-феноменологические определяющие уравнения пластического течения для объемного напряженно-деформированного состояния. Необходимые предпосылки для их получения в рамках идеологии математической теории пластичности бы-
ли изложены в [9]. Однако некоторые формулировки нуждаются в уточнении.
В работе [19] было показано, что в случае вязкопластического тела (деформация металлов при повышенных температурах) уменьшение напряжения не переводит его из пластического состояния в упругое, а сопровождается уменьшением скорости деформации. Поэтому процесс разгрузки для него не имеет места. Дополнительно к процессам нагружения и нейтрального нагружения для вязкопластического тела были выведены понятия термических и динамических упрочнения и возврата [9]. Введены также понятия мгновенных на расчетном шаге функций нагружения, условия пластичности Губера-Мизеса и сформулирован новый для вязкопластического тела принцип максимума работы пластической деформации, учитывающий влияние истории нагружения.
При горячей деформации возникающая кристаллографическая текстура, обуславливающая деформационную анизотропию, «размывается» динамической рекристаллизацией и материал остается изотропным.
Возьмем мгновенную на расчетном шаге функцию нагружения в виде
аі-Яд) ^,(Т’Яд)
/9-І
фы
= 0 (1.28)
и, соответственно, мгновенное условие пластичности Губера-Мизеса
= ~ 1 8 ¡-Яд) с^8>Я.д)) \8і-Яд) с^8>Я.д)
/9-1
ф(9) (52^(9) + СІЄ{д)
= 0,
(1.29)
гд^'(9) и ^іЯд) — девиаторы тензоров а^(ду .
Из (1.29) следует
а(д) = а(д) + ^(9) = Ф(9) ^*(9) + (1Ч.9)^ •
Как видно, выбор условия пластичности в виде (1.29) означает, что интенсивность напряжения на шаге нагружения <7 определяется оператором (1.7)-і-(1Л6).
390
290
є
б
в
Ассоциированный закон течения вязкопластического тела будет иметь вид
д
аЬ{9) +
°Щд) + с1'°И
ЧІ.9)
()) ^,Х
(9)
х с!є.ІШ - <1,\{д)/{д)] = 0:
(ТЬ{д) + с^,(Т'Яд)
д) I ^'ч//(д)) ■ Множитель Лагранжа определяется как
^\д) \/ 2"
сіє
(д)
°ї(д)/0{(Т}лд)+(і'(Т>Яд)
СІЄ
(д)
°кд)
аЬ{д) + ^,(Т>Яд)
2 (а(д) + с1(Т^в)
Определяющие соотношения для вязкопластического тела выглядят следующим образом:
ск
(д)
. а
(д)
с1а
(д)
8Ь(д)
(1.30)
где и на каждом шаге нагруже-
ния рассчитываются по (1.16) и (1.7) соответственно.
Заметим, что из (1.30) как частный случай следуют определяющие соотношения классической теории течения [1]
3 ск
с''£и = 2~^3'Я
Действительно, приняв независимость от истории нагружения, получим
и
.
При решении прикладных задач вязкопла-стичности с применением численного метода конечных элементов, который реализует необходимый пошаговый алгоритм расчета, удобно использовать (1.30) в виде двух уравнений [9]
_ •* ^£($) і >Я9) 2 с1ст{д} ,Я9)'
З (Іє{д) 2 п!\
(д)
Подчеркнем также, что модель, включая расчет текущей плотности неподвижных дислокаций , позволяет прогнозировать текущее структурное состояние материала, в частности, линейный размер зерен и субзерен по известной в металлофизике формуле
= 1/у/р* [13].
При холодной направленной пластической деформации в металлах возникает деформационная анизотропия, одним из проявлений которой является эффект Баушинге-ра. Анизотропия определяет влияние на ст характера траектории нагружения. В работе [13] предложена физико-феноменологическая модель эффекта Баушингера, имеющая вид
а = /3тОЬ |
Ряс
Ае+ ехр(є) — 1
охр (є) 1,5АЬсхр(є)
(1.31)
Уравнение (1.31) описывает диаграмму деформирования материала с учетом предварительной деформации с обратным знаком . Для определения коэффициента предложен достаточно простой эксперимент, включающий волочение или прямое прессование (вид деформации — простое растяжение) прутка со степенью деформации из диапазона и последующую осадку
изготовленных из него образцов с определением предела текучести на сжатие . Значение находится по формуле
А =
О7!)2 - {/ЗтОЬ)2 ряо (/3тОЪ)2 е+
В общем случае сложного нагружения и деформирования, сопровождающегося сменой знака деформации, (1.31) можно записать в виде
Ряо
Ае
охр(є)
охр(є) — 1 I2 1,5АЬсхр(е) /
(1.32)
Для получения определяющих соотношений изотропного материала с анизотропным упрочнением примем условие текучести Губера-Мизеса в виде
2 Vз з гз гз
- ІФ. (е)]* - [Фг (е)]* - I*,,*,, - [Ф (е)]* - О,
(1.33)
где в соответствии с (1.32):
3 _ _ _ ,отГг: (Рио + Лє\ /2
2 ' ■ “СТ1“/3тСЧ^фМ
3 (2) (2)
2'И'И = °2 =
]_
а т. ( СХР (Е) ~ 1 ^ '2 ,3т6Чі,5ЬЛохр(£)і
I1) С2) I1) С2)
ЯгЗ > Ягз '.8ц— девиаторы тензоров 0\- , о\- , 0^
гг I1) , (2)
соответственно. При этом оц = а- + а- и,
К) •
Ассоциированный закон течения [1] д
да,
{аі;і<І: і;і - (Щ) = О,
*?
(1:1! = Зг/Л*;; = Зг/Л Как известно [1],
(1.34)
,, /3 сіє 1 сіє
ал = \ -- ,
2 / с>/ с>/
да-,да-,
2 а
(1.35)
Подстановкой (1.35) в (1.34) окончательно получаем
л,.. _ О1) . ,('2)
4 2 гг І 'і ч
(1.36)
где сг определяется для каждой по (1.32).
Видно, что по форме уравнения (1.36) совпадают с известными феноменологическими определяющими соотношениями теории пластичности изотропного материала с анизотропным упрочнением. Однако значение (1.36) состоит в том, что они получены в рамках единой концепции — физико-математической теории необратимых деформаций металлов.
2. ФИЗИКО-МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПОЛЗУЧЕСТИ МЕТАЛЛОВ
Современные расчетные методы феноменологической теории ползучести металлов (технические теории течения, упрочнения и старения) базируются, в основном, на фундаментальных результатах, развитых и обобщенных в [2]. Во введении были отмечены их недостатки. Проблема ползучести — это проблема вязкопластичности. Поэтому можно полагать, что на основе изложенной выше модели вязкопластичности могут быть получены основные уравнения дислокационной ползучести, контролируемой термодинамически активированными возвратом и скольжением дислокаций.
Для вывода уравнения, описывающего ползучесть при одноосном напряженном состоянии, из (1.6) выразим скорость деформации
р{Т)0{Т)Ъл
¿1д) = ¿*Ьу//3,П9_1)СХр
кТ,
(9)
а
(д)
р(Т)тО(Т)Ъ^Цу
(2.1)
где верхний индекс « », здесь и далее, обозначает макрохарактеристики ползучести.
Для введения в модель основной переменной — времени — в уравнение (1.8) подставим и изменение плотно-
сти неподвижных дислокаций при ползучести под действием приложенного напряже-
ния найдется как
^Ря(д)
- {р,Цд-1)У/г "1)Ьх
3/2 .
х охр
р(Т)С(Т)Ъл — а^Ь2 /Шу/рГо
кТ,
(9)
М(д}.
(2.2)
Для численного расчета кривых ползучести уравнения (2.1) и (2.2) дополним
очевидными соотношениями:
Р я(д) Ри(д-1) + (2.3)
СІЄ1в) = ё[9)М(9)- (2.4)
є(д) =є(9-1) +(1є(д): (2.5)
^{д) = ^{д-1) + ^'{д}- (2.6)
Оператор (2.1)-(2.6) описывает ползучесть металлов как в условиях Т,а° = const, так и T(t), crc(t). Следует заметить, что полученная модель объединяет современные технические теории ползучести: теорию старения , течения и
упрочнения .
Условием установившейся (стационарной) стадии ползучести при const будет и ползучесть протекает при
постоянной минимальной скорости ¿£„п. Для этой стадии из (2.2) имеем
— = pivDb х ЬХ
х exp
Р(Т)С(Т)ЪЛ -
kT
(2.7)
Подстановкой (2.1) в (2.7) получаем уравнение для определения стационарной плотности неподвижных дислокаций на установившейся стадии, которое после преобразований имеет вид
ln (Xpsi^vb/s*) =
o'
b2
mkT
1
1
■ (2.8)
Точное решение этого уравнения представляет собой сложное и громоздкое выражение. Поэтому ограничимся приближенным решением. Учитывая, что логарифмическая функция медленно изменяется с изменением аргумента при больших его значениях, примем
и приближенное решение уравнения (2.8) будет иметь вид
\/pS — y/Pso
ткТу/рГоЫ {ХряоииЬ1 е*
а
сЬ2
1 , (2.9)
где — стационарная плотность дислокаций на установившейся стадии.
Экспериментальными исследованиями влияния деформации ползучести на структуру металлов установлено, что общей закономерностью является формирование к началу установившейся стадии субзеренной структуры [14]. Для оценки линейных размеров субзерен предложено эмпирическое соотношение
d =
10,5Gb
(Tl
(2.10)
Уравнение (2.10) позволяет оценить d как
d =
1
— (\/ Pso) х
'М
mkT y/pZ In (Xpsoi/Db/£*
а
сЬ2
(2.11)
Видно, что (2.11) полнее отражает влияние различных факторов на структуру металла, чем (2.10). Подстановка (2.9) в (2.1) дает уравнение для определения минимальной скорости ползучести . При проверке адекватности модели использовали экспериментальные результаты работы [15]. На рис. 2 в качестве примера приведены экспериментальные точки и теоретические (сплошные линии) кривые ползучести стали 08 при
и (кривая
1) и ас = 75 МПа (кривая 2). Характеристики стали 08: рв0 = 2,12 • Ю10 см-2; А = 3,6 х см; МПа;
МПа. Общий результат проверки состоит в том, что модель весьма удовлетворительно описывает кривые ползучести различных материалов в широком диапазоне и .
Рис. 2
Одноосная модель ползучести (2.1)-(2.6), (2.9), (2.11) выведена из одноосной модели пластичности (1.6)—(1.16). Следовательно, для нее остаются в силе определяющие соотношения (1.30). Из этого следует, что в рамках предлагаемой модели постановка и решение краевых задач ползучести эквивалентны постановке и решению краевых задач пластичности при повышенных температурах.
£
t
ЗАКЛЮЧЕНИЕ
На основе объединения микро- и макропредставлений, методов микро- и макроописания необратимых деформаций в рамках единой концепции разработана модель дислокационной пластичности и ползучести металлов. В отличие от классических феноменологических теорий физико-математическая теория последовательно учитывает историю нагружения и позволяет прогнозировать структуру и, следовательно, свойства материалов после деформации. В связи со спецификой уравнений модели (дифференциальная форма), краевые задачи пластичности и ползучести целесообразно решать численным методом конечных элементов. Модель позволяет естественным образом осуществлять линеаризацию определяющих соотношений методом мгновенных касательных модулей при численном пошаговом расчете необратимых деформаций.
СПИСОК ЛИТЕРАТУРЫ
1. Ишлинский, А. Ю. Математическая теория пластичности / А. Ю. Ишлинский, Д. Д. Ивлев. М.: Физматлит, 2001. 704 с.
2. Работнов, Ю. Н. Ползучесть элементов конструкции / Ю. Н. Работнов. М. : Наука, 1966. 752 с.
3. Батдорф, С. Б. Математическая теория пластичности на основе концепции скольжения / С. Б. Батдорф, Б. А. Будянский // Механика. М.: Мир, 1962. № 1. С. 135-155.
4. Киселев, С. П. Модель упругопластического деформирования материалов на основе калибровочной теории дефектов с учетом диссипации энергии / С. П. Киселев // Прикладная механика и техническая физика. 2004. Т. 45, №2. С. 177-187.
5. Гриняев, Ю. В. Полевая теория дефектов на мезоуровне / Ю. В. Гриняев, В. Е. Панин // Доклады РАН. 1997. Т. 353, № 1. С. 37-39.
6. Лихачев, В. А. Структурно-аналитическая теория прочности / В. А. Лихачев, В. Г. Малинин. СПб.: Наука, 1993. 471 с.
7. Грешнов, В. М. Физико-феноменологическая модель сопротивления металлов
пластической деформации для расчета технологических процессов обработки металлов давлением. Сообщ. 1. Постановка задачи и вывод общего уравнения / В. М. Грешнов, Ф. Ф. Сафин, М. В. Грешнов // Проблемы прочности. 2002. № 6. С. 107-115.
8. Грешнов, В. М. Физико-феноменологическая модель сопротивления металлов
пластической деформации для расчета технологических процессов обработки металлов давлением. Сообщ. 2. Частные случаи модели и ее экспериментальная проверка /
В. М. Грешнов, Ф. Ф. Сафин, М. В. Грешнов // Проблемы прочности. 2003. № 1. С. 87-97.
9. Грешнов, В. М. Модель вязкопластического тела с учетом истории нагружения / В. М. Грешнов // Механика твердого тела. 2005. №2. С. 117-125.
10. Колмогоров, В. Л. Механика обработки металлов давлением / В. Л. Колмогоров. Екатеринбург : Изд-во Уральск. гос. техн. ун-та — УПИ, 2001. 836 с.
11. Штремель, М. А. Прочность сплавов. Ч. 1. Дефекты решетки. М. А. Штремель. М. : Металлургия, 1982. 280 с.
12. Сорокин, В. Г. Марочник сталей и сплавов / под общ. ред. В. Г. Сорокина. М. : Машиностроение, 1989. 640 с.
13. Грешнов, В. М. Прогнозирование разрушения металлов в процессах холодной пластической деформации. Сообщ. 2. Учет анизотропного упрочнения и экспериментальная проверка модели пластической деформации и разрушения / В. М. Грешнов [и др.] // Проблемы прочности. 1999. № 2. С. 74-84.
14. Пуарье, Ж. П. Высокотемпературная пластичность кристаллических тел / Ж. П. Пуарье. М.: Металлургия, 1982. 272 с.
15. Дегтярев, В. П. Пластичность и ползучесть машиностроительных конструкций / В. П. Дегтярев. М. : Машиностроение, 1967. 131 с.
ОБ АВТОРАХ
Грешнов Владимир Михайлович, проф., зав. каф. обр-ки металлов давлением. Дипл. инж.-мех. (УАИ, 1974). Д-р физ.-мат. наук по физ. тверд. тела (ИПМ АН Украины, 1992). Иссл. в обл. физ. и мех. пласт. деформ.
Патяева Ирина Владимировна, асп. той же каф. Дипл. магистр техн. и технол. по спец. «Машины и технология обработки металлов давлением» (УГАТУ, 2007). Готовит дис. в обл. физ.-мат. теории пластичности.
Сидоров Валерий Евгеньевич, ст. преп. той же каф. Дипл. инж.-мех. по спец. «Авиационные двигатели» (УАИ, 1979). Иссл. в обл. физ.-мат. теории ползучести.