Основы механики твердого деформируемого тела
.pdf
Глава 2 |
751 |
Эти равенства выражают так называемый ассоциированный закон течения, ибо они связаны (ассоциированы) с конкретным условием текучести (2.7).
Множитель dλ, входящий в ассоциированный закон течения, остается неопределенным. Найти его значение, не конкретизируя физическую модель материала и не указывая историю нагружения, нельзя. И в самом деле, как об этом говорилось в п. 2.2, условие текучести (2.7) зависит от меры упрочнения q, а последняя может быть вычислена только при известной траектории (истории) нагружения и выбранной модели материала. Различные траектории нагружения с общими началом и концом приводят, вообще говоря, к различным условиям текучести и, как следствие, к различным решениям уравнений (2.12).
Если в формулу (1.36a), определяющую приращение интенсивности пластической деформации, ввести величины dεip, которые даются формулой (2.12), то для множителя Лагранжа получится следующее выражение:
3 |
|
|
|
|
|
|
|
|
|
|
|
интp |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
dε |
|
|
|
|
|
|
|
|
|
|
(2.13) |
||||||
dλ = √ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
2 |
|
|
∂f |
|
∂f |
2 |
∂f |
|
|
∂f |
2 |
∂f |
|
∂f |
|
2 |
|||||||||||
|
|
|
|
− |
|
|
+ |
|
|
− |
|
|
+ |
|
|
− |
|
|
|
|
|
||||||
|
|
∂σ |
1 |
∂σ |
2 |
∂σ |
2 |
∂σ |
3 |
∂σ |
3 |
∂σ |
1 |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Ясно, что найти отсюда множитель dλ нельзя, ибо величина dεинтp неизвестна, но из равенства (2.13) можно установить неотрицательность этого множителя: dλ ≥0. Это наблюдение, в свою очередь, приводит к следующим утверждениям.
a) При активном нагружении тела
∂f
∂σi dσi > 0, f (σi) = 0, dλ > 0.
b) В случае пассивного нагружения (разгрузки)
∂f
∂σi dσi < 0, f (σi) = 0, dλ = 0.
c) При нейтральном нагружении
∂f dσi = 0, f (σi) = 0, dλ = 0. ∂σi
Если материал является идеально пластическим (нет упрочнения), то при его нагружении
∂f dσi = 0, f (σi) = 0, dλ > 0 ∂σi
752 |
|
Часть VII |
и |
∂f |
|
|
dσi < 0, f (σi) = 0, dλ = 0 |
|
|
|
|
|
∂σi |
|
вслучае разгрузки.
2.5.Теория течения. При помощи ассоциированного закона текучести (2.12) можно получить разные теории пластичности, отвечающие той или иной функции (2.7), которая определяет как поверхность текучести, так и закон упрочнения. Одну из таких теорий, называемую теорией течения, получают в предположении, что в процессе нагружения поверхность текучести равномерно расширяется (изотропное упрочнение). В этом случае двум равным по величине и обратным по знаку нагружениям отвечают напряжения одной и той же интенсивности, что упрощает построение теории. Однако изотропное упрочнение не описывает эффект Баушингера, т. е. вступает в противоречие с экспериментально установленным фактом. Это несущественно в случае, если история нагружения не содержит разгрузки из пластической зоны, но тогда, когда такая разгрузка была, результатами, полученными с помощью теории течения, надо пользоваться с осторожностью.
|
|
|
Φ dεинтp |
2 = σинт2 . |
(2.14a) |
Пусть условие текучести (2.7) определяется критерием Губера – Мизеса |
|||||
(2.4): |
|
|
|
|
|
f (σ1) ≡ |
1 |
[(σ1 |
−σ2)2 + (σ2 −σ3)2 |
+ (σ3 −σ1)2] − (Φ(q))2 = 0 |
(2.14) |
2 |
|||||
и мера q берется в форме (2.6) параметра Одквиста. Тогда (см. также равенство (1.26))
В конце п. 2.2 говорилось о независимости вида функции Φ(q) от вида напряженного состояния, что позволяет находить ее по диаграмме "σ–ε" простого растяжения, т. е. при σ1 = σ, σ2 = σ3 = 0. В состоянии текучести коэффициент Пуассона любого материала равен 1/2, поэтому dε1p = dεp, dε2p = dε3p = −0, 5dεp. Равенства (1.26) и (1.36a) в этом случае дают
|
|
σинт = σ, |
|
интp = dεp. |
|
||
|
|
dε |
|
||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
интp = dε p, Φ(ε p) = σ. |
(2.14b) |
|||
|
|
dε |
|||||
Более конкретную информацию о функции Φ(q) можно получить только при наличии диаграммы "σ–ε" (см. рис. 2.8).
Глава 2 |
753 |
Ассоциированный закон течения (2.12) содержит частные производные от функции f (σi) по напряжениям σi. Они же (производные) определяют множитель Лагранжа (2.13). Найти величины ∂f /∂σi можно, обратившись к соотношению (2.14). Например,
∂f
∂σ1 = (σ1 −σ2) − (σ3 −σ1) = 2σ1 −σ2 −σ3 −σ1 +σ1 = 3(σ1 −σ0),
где σ0 – среднее напряжение (1.21). Таким образом,
∂f |
= 3(σi − σ0) |
(2.15) |
∂σi |
и множитель Лагранжа имеет представление
1 |
|
|
|
|
интp |
|
||
|
|
|
dε |
(2.13a) |
||||
dλ = √ |
|
|
|
|
. |
|||
|
|
|
||||||
2 |
|
(σ1 −σ2)2 +(σ2 −σ3)2 +(σ3 −σ1)2 |
||||||
Эта формула с учетом равенства (1.26) записывается более компактно:
|
1 |
|
|
|
интp |
. |
dλ = |
|
|
dε |
|||
|
|
|
||||
|
2 σинт |
|||||
При таком значении величины dλ и производных (2.15) ассоциированный закон течения (2.12) принимает вид:
|
|
3 |
|
|
интp |
|
|
p |
= |
|
dε |
(σi −σ0). |
(2.16) |
||
dεi |
|
|
|
|
|||
2 |
|
σинт |
|||||
Отсюда следует, что компоненты тензора приращений пластических деформаций пропорциональны компонентам (σi −σ0) девиатора напряжений. Согласно закону Гука (I.6.6), упругая доля приращения деформации ε1 определяется равенством:
dε1(e) = E1 [dσ1 − ν(dσ2 + dσ3)].
Но σ2 +σ3 = 3σ0 −σ1 и тогда (см. формулу (1.37))
dε1 |
= |
E (1+ν)dσ1 − 3νdσ0 |
|
= 2G dσ1 − |
1+ν dσ0 |
. |
|||||||||
(e) |
|
1 |
|
|
|
|
|
1 |
|
|
|
3ν |
|
|
|
Стало быть, |
|
|
|
1 |
dσi − |
|
3ν |
. |
|
|
|
||||
|
|
|
dεi(e) = |
|
|
|
|
||||||||
|
|
|
|
|
|
dσ0 |
|
|
|
||||||
|
|
|
2G |
1+ν |
|
|
|
||||||||
Глава 2 |
755 |
соотношениями физических уравнений (2.17) и составляет основу деформационных теорий пластичности.
Получить искомые конечные соотношения между напряжениями и деформациями можно лишь при условии, что нагружение является простым: отказа от рассмотрения траекторий нагружения, не включающих разгрузку, еще не достаточно. При простом нагружении компоненты девиатора напряжений возрастают пропорционально некоторому параметру α:
σi − σ0 = α(σi − σ0) , 0 ≤ α ≤ 1.
Значком отмечены компоненты напряженного состояния в конце нагружения. Тогда (см. формулу (2.16))
dεip = 32 (σi −σ0) dεpинт.
Следовательно,
ε p = |
3 |
|
(σi −σ0) |
|
|
|
p . |
(2.18) |
|
|
dε |
||||||||
2 |
|||||||||
i |
|
σинт |
|
|
инт |
|
После подстановки величины εip , даваемой этим равенством, в формулу (1.33a) для интенсивности пластических деформаций и выполнения ряда преобразований с использованием выражения (1.26) получится, что
εip = dεинтp . (2.19)
Но тогда dεpинт = dεинтp . Стало быть, в случае простого нагружения интенсивность dεинтp приращения пластических деформаций равна приращению dεинтp интенсивности пластических деформаций.
Равенство (2.19), являющееся следствием простого нагружения, позволяет придать физическому уравнению (2.18) конечную форму (индекс теперь может быть опущен):
p |
3 |
|
p |
|
|
|
= |
|
|
εинт |
(σi − σ0). |
(2.20) |
|
εi |
|
|
|
|||
2 |
σинт |
|||||
Из формул (I.6.6)1, (1.38) и (1.37) следует (см. также вывод уравнений (2.17))
ε |
|
ε |
|
= |
1 |
[σ |
|
− |
ν(σ |
+σ |
+σ |
)] |
− |
1−2ν |
σ |
|
= |
1 |
(σ |
1 − |
σ |
). |
||
|
|
E |
|
E |
|
2G |
||||||||||||||||||
|
1 |
− |
0 |
|
|
1 |
2 |
|
3 |
1 |
|
|
0 |
|
|
0 |
|
|||||||
Таким образом, |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
εie = |
|
(σi |
−σ0) + ε0. |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
2G |
|
|
|
|
|
|
|
||||||||
756 |
Часть VII |
Эту упругую составляющую полной деформации εie с учетом зависимости σинт = 3Gεинт, даваемой формулой (1.41), можно представить в виде
e |
|
3εинтe |
|
εi |
= |
|
(σi −σ0) + ε0. |
2σинт |
|||
Тогда (см. соотношение (2.20))
e |
p |
|
3 |
εинтe +εинтp |
|
|
|
εi = εi |
+ εi |
= ε0 + |
|
|
|
(σi −σ0). |
(2.21) |
2 |
|
σинт |
|||||
При помощи формул (1.33), (1.33a) и εi = εei +εpi можно убедиться в том (формула для величины εинтe подобна зависимости (1.33a)), что
εинт = εинтe
т. е. интенсивности полных, упругих и пластических деформаций обладают теми же аддитивными свойствами, что и сами деформации. Следовательно,
εi − ε0 = |
3 |
|
εинт |
(σi −σ0), i = 1, 2, 3. |
(2.21a) |
2 |
|
σинт |
В осях 0x, 0y, 0z связь между напряжениями и деформациями описывается уже шестью уравнениями типа (2.21a):
εx − ε0 |
= 2 σинт (σx −σ0), |
|
||||||
|
3 |
|
εинт |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
· · · |
· · · |
|
· · · |
· · · |
|
|
||
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.22) |
γzx = 3 |
|
εинт |
τzx. |
|
|
|
||
|
|
|
|
|||||
|
|
|
σинт |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В теории малых упругопластических деформаций предполагается, что при пластическом течении объем тела не меняется, т. е. ε0 = 0. Для несжимаемого материала коэффициент ν поперечной деформации равен 1/2. Это значение коэффициента Пуассона, установленное при упругой работе материала, остается таким же и при малых пластических деформациях. Если же остаточные деформации велики, то величину ν вычисляют по формуле
ν = 1 − 1−2ν σ ,
2 2E ε
которая здесь приводится без вывода. В случае идеального пластического тела с ростом деформации ε отношение σ/ε стремится к нулю, а коэффициент ν – к 1/2.
Глава 2 |
757 |
Предположений о том, что нагружение является простым, история нагружения не содержит разгрузки и объем тела при пластическом деформировании не меняется, еще не достаточно для построения деформационной теории пластичности. Дело в том, что в уравнения (2.22) входит отношение εинт/σинт, вычисление которого требует знания связи между величинами σинт и εинт. Эту связь устанавливают по результатам испытаний материала на растяжение, что предполагает наличие единой кривой, о которой шла речь в п. 2.2. График зависимости σинт(εинт) называют диаграммой деформирования материала. Гипотеза о существовании единой кривой, т. е. о независимости функции σинт(εинт) от вида деформации, была выдвинута в начале XX века немецким исследователем П. Людвиком.
Из уравнений (2.22) вытекают две группы равенств:
ε1 |
−ε0 |
= |
ε2 |
−ε0 |
= |
ε3 |
−ε0 |
= |
k |
, |
γ1 |
= |
γ2 |
= |
γ3 |
= k; |
k = |
3εинт |
, |
|
|
|
|
|
−σ0 |
|
|
|
|
|
|||||||||
σ1 |
−σ0 |
σ2 |
−σ0 |
σ3 |
2 |
|
τ1 |
τ2 |
τ3 |
|
σинт |
||||||||
фиксирующих пропорциональность компонентов девиаторов напряжений и деформаций, а также главных касательных напряжений и главных сдвигов. При помощи первой группы этих формул можно следующим образом трансформировать коэффициент Лоде – Надаи (1.30):
µσ = |
2σ2 −σ1 −σ3 |
= |
2(σ2 −σ0)−(σ1 −σ0)−(σ3 −σ0) |
= |
|||
|
(σ1 −σ0)−(σ3 −σ0) |
||||||
|
σ1 −σ3 |
|
|||||
= |
2(ε2 −ε0)−(ε1 −ε0)−(ε3 −ε0) |
= |
2ε2 −ε1 −ε3 |
. |
|
||
(ε1 −ε0)−(ε3 −ε0) |
|
|
|||||
|
|
ε1 −ε3 |
|
||||
Стало быть (см. равенство (1.35)), |
|
|
|
|
|||
|
|
|
µσ = µε. |
|
|
|
(2.23) |
Таким образом, если теория малых упругопластических деформаций верна, то коэффициенты Лоде – Надаи по напряжениям и деформациям должны совпадать не только в случае простейших напряженных состояний (чистые растяжение, сжатие, сдвиг), но и при любых иных напряженных состояниях. Это обстоятельство используется для экспериментальной проверки деформационной теории пластичности (см. далее п. 4.4).
Уравнения (2.22) были получены в первой половине прошлого века академиком А. А. Ильюшиным. Им же был разработан метод их численного решения.
Теория малых упругопластических деформаций, основу которой составляют уравнения (2.22), является самым простым вариантом деформационной теории пластичности. В более общем случае связь между компонентами
758 Часть VII
девиаторов деформаций и напряжений дается равенством (ср. с формулой (2.21a))
εi − ε0 = (σi − σ0)ψ(σi, εi), |
(2.21b) |
где ψ – некоторая скалярная функция напряжений и деформаций. При ψ = 1 уравнение (2.21b) приобретает форму закона Гука для девиаторов, а при
ψ= 3 εинт 2 σинт
соотношение (2.21b) приводит к основному уравнению (2.21a) теории малых упругопластических деформаций. В случае, когда нагружение близко к простому, эта теория дает достаточно правильную картину распределения напряжений в теле.
2.7. Иллюстративный пример. Приводимый ниже пример предназначен не столько для демонстрации процедур решения задач теории пластичности методами, описанными в предыдущих пунктах, сколько для сопоставления напряженных состояний, к которым приводят теории течения и малых упругопластических деформаций. Дело в том, что на простом примере невозможно показать даже самые основные элементы вычислительных алгоритмов теории пластичности, а обращение к сложным примерам существенно загромоздило бы изложение.
На рис. 2.12 показана тонкостенная труба, выполненная из несжимаемого идеального упругопластического материала, по торцам которого приложены растягивающие силы и крутящие моменты. При δ = const R касательные напряжения кручения τM равномерно распределены и по толщине скорлупы, и вдоль ее средней линии (см. п. II.7.6). Равномерно распределены по поперечному сечению и нормальные напряжения σN растяжения. Кроме того, напряжения τM и σN не меняются вдоль оси трубы. Следовательно, напряженное состояние рассматриваемого тела однородно (к тому же оно является двумерным). На площадке dF , как видно по рис. 2.12, действуют только напряжения τM, σN, а на площадках, параллельных плоскости, проходящей через векторы напряжения τM и σN, напряжений нет вообще, поэтому в той же плоскости τMσN расположены и векторы главных напряжений σ1 и
Глава 2 |
759 |
σ2, значения которых находят по формулам (I.2.13c):
σ1 |
= 2 |
σN + σN2 +4τM2 |
, |
σ2 |
= 2 |
σN − σN2 +4τM2 |
. |
(2.24) |
||
|
1 |
|
|
|
|
1 |
|
|
|
|
Аналогичный вид имеют и формулы для главных деформаций:
ε1 |
= 2 |
εN + εN2 +4γM2 |
, |
ε2 |
= 2 |
εN − εN2 +4γM2 |
. |
(2.25) |
||
|
1 |
|
|
|
|
1 |
|
|
|
|
В двумерной задаче σ3 = ε3 = 0.
Напряжениям (2.24) отвечают следующие интенсивности (1.26) и (1.28) нормальных и касательных напряжений:
|
|
|
|
|
|
σинт = |
σN2 +3τM2, τинт = σN2/3+τM2. |
(2.26) |
|||
Аналогично, но уже при помощи формул (2.25), (1.33) и (1.36a), устанавливаются значения интенсивностей деформаций и приращений пластических деформаций:
εинт = |
|
, |
|
интp = |
|
. |
|
|
εN2 +γM2 /3 |
(dεN)2 +(dγM)2/3 |
(2.27) |
||||||
dε |
Так как физические уравнения теорий течения и малых упругопластических деформаций содержат интенсивности напряжений, деформаций и приращений пластических деформаций, то без помощи выражений (2.26) и (2.27) задачу о распределении напряжений в упругопластическом теле не решить.
Согласно условию (2.4) начала текучести Губера – Мизеса, а именно оно принимается при решении рассматриваемой задачи, материал переходит в пластическое состояние, если
σинт = σT или τинт = τT. |
(2.28) |
||||
Из формул (2.26) и (2.28) следует, что |
|
||||
σT/τT = σинт/τинт = √ |
|
. |
|
||
3 |
|
||||
Значит, |
|
||||
σT = √ |
|
τT, |
(2.29) |
||
3 |
|||||
σN2 + 3τM2 = σT2. |
(2.26a) |
||||
Для несжимаемого материала ν = 1/2, а в этом случае, согласно формуле (1.37), E = 3G. Тогда в момент наступления текучести
σT = EεT, τT = GγT = EγT/3
760 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Часть VII |
или (см. формулу (2.29)) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
= |
σT |
|
|
γT = √ |
|
|
εT. |
|
|
|
(2.30) |
||
|
|
εT |
, |
|
3 |
|
|
|
||||||||
|
|
|
|
|
|
|||||||||||
|
|
|
|
E |
|
|
|
|
|
|
|
|
|
|||
Дальнейшие выкладки удобнее выполнять при помощи обозначений |
||||||||||||||||
σ = |
σN |
, |
τ = |
τM |
, |
ε = |
εN |
, |
γ = |
γM |
, |
(2.31) |
||||
|
|
|
|
|||||||||||||
|
σT |
|
|
τT |
|
εT |
|
γT |
|
|||||||
которые позволяют записать условие начала текучести (2.26a) и формулы (2.26)1 и (2.27) для интенсивностей нормальных напряжений и деформаций в виде:
|
|
|
|
|
|
σ2 + τ 2 = 1, |
|
|
|
(2.32) |
||||
σT |
= σ2 + τ 2, |
εT |
= |
ε2 + γ2, |
|
|
= (dε)2 + (dγ)2. |
(2.33) |
||||||
εT |
||||||||||||||
σинт |
|
|
|
εинт |
|
|
|
|
dεинтp |
|
|
|
|
|
Теперь можно обратиться к непосредственному решению задачи. Сначала оно выполняется средствами теории течения.
2.7.1. Решение на основе теории течения. Из физических уравнений (2.17) этой теории нетривиальными остаются только два уравнения: одно для нормальных напряжений, другое – для касательных. При этом упругие составляющие деформаций отбрасываются (идеально жесткое пластическое тело):
|
|
|
3 |
|
интp |
|
|
|
|
|
|
|
|
|
|
интp |
|
|
|
|
||||
dεN = |
dε |
(σN −σ0), dγM = 3 |
dε |
τM. |
(2.34) |
|||||||||||||||||||
2 |
|
σинт |
σинт |
|||||||||||||||||||||
В предельном состоянии, что видно по формуле (2.28)1, σинт = σT. Кроме |
||||||||||||||||||||||||
того, при напряжениях (2.24) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
σ0 = |
1 |
|
(σ1 + σ2 |
+ σ3) = |
1 |
σN. |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
3 |
|
|
|
|
3 |
|
|
|
|
|
|
|
||||||||||||
Тогда (учитываются обозначения (2.31) и формулы (2.33)3, (2.30)2)): |
|
|||||||||||||||||||||||
dε = σ |
|
|
|
, dγ = τ |
|
|
. |
|
||||||||||||||||
(dε)2 + (dγ)2 |
(dε)2 + (dγ)2 |
(2.34a) |
||||||||||||||||||||||
Отсюда следует, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
1 |
|
|
|
|
τ = |
|
dγ/dε |
|
|
|
(2.35) |
|||||||||||||
σ = |
|
|
, |
|
|
. |
||||||||||||||||||
|
|
|
|
|||||||||||||||||||||
1+(dγ/dε)2 |
1+(dγ/dε)2 |
|||||||||||||||||||||||
