Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Основы механики твердого деформируемого тела

.pdf
Скачиваний:
128
Добавлен:
14.05.2015
Размер:
26.55 Mб
Скачать

Глава 2

751

Эти равенства выражают так называемый ассоциированный закон течения, ибо они связаны (ассоциированы) с конкретным условием текучести (2.7).

Множитель , входящий в ассоциированный закон течения, остается неопределенным. Найти его значение, не конкретизируя физическую модель материала и не указывая историю нагружения, нельзя. И в самом деле, как об этом говорилось в п. 2.2, условие текучести (2.7) зависит от меры упрочнения q, а последняя может быть вычислена только при известной траектории (истории) нагружения и выбранной модели материала. Различные траектории нагружения с общими началом и концом приводят, вообще говоря, к различным условиям текучести и, как следствие, к различным решениям уравнений (2.12).

Если в формулу (1.36a), определяющую приращение интенсивности пластической деформации, ввести величины ip, которые даются формулой (2.12), то для множителя Лагранжа получится следующее выражение:

3

 

 

 

 

 

 

 

 

 

 

 

интp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.13)

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

∂f

 

∂f

2

∂f

 

 

∂f

2

∂f

 

∂f

 

2

 

 

 

 

 

 

+

 

 

 

 

+

 

 

 

 

 

 

 

 

 

∂σ

1

∂σ

2

∂σ

2

∂σ

3

∂σ

3

∂σ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ясно, что найти отсюда множитель нельзя, ибо величина интp неизвестна, но из равенства (2.13) можно установить неотрицательность этого множителя: dλ ≥0. Это наблюдение, в свою очередь, приводит к следующим утверждениям.

a) При активном нагружении тела

∂f

∂σi i > 0, f (σi) = 0, dλ > 0.

b) В случае пассивного нагружения (разгрузки)

∂f

∂σi i < 0, f (σi) = 0, dλ = 0.

c) При нейтральном нагружении

∂f i = 0, f (σi) = 0, dλ = 0. ∂σi

Если материал является идеально пластическим (нет упрочнения), то при его нагружении

∂f i = 0, f (σi) = 0, dλ > 0 ∂σi

752

 

Часть VII

и

∂f

 

 

i < 0, f (σi) = 0, dλ = 0

 

 

 

∂σi

вслучае разгрузки.

2.5.Теория течения. При помощи ассоциированного закона текучести (2.12) можно получить разные теории пластичности, отвечающие той или иной функции (2.7), которая определяет как поверхность текучести, так и закон упрочнения. Одну из таких теорий, называемую теорией течения, получают в предположении, что в процессе нагружения поверхность текучести равномерно расширяется (изотропное упрочнение). В этом случае двум равным по величине и обратным по знаку нагружениям отвечают напряжения одной и той же интенсивности, что упрощает построение теории. Однако изотропное упрочнение не описывает эффект Баушингера, т. е. вступает в противоречие с экспериментально установленным фактом. Это несущественно в случае, если история нагружения не содержит разгрузки из пластической зоны, но тогда, когда такая разгрузка была, результатами, полученными с помощью теории течения, надо пользоваться с осторожностью.

 

 

 

Φ инт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, поэтому 1p = p, 2p = 3p = 0, 5p. Равенства (1.26) и (1.36a) в этом случае дают

 

 

σинт = σ,

 

интp = p.

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

интp = p, Φ(ε p) = σ.

(2.14b)

 

 

Более конкретную информацию о функции Φ(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

 

 

 

 

(2.13a)

=

 

 

 

 

.

 

 

 

2

 

(σ1 −σ2)2 +(σ2 −σ3)2 +(σ3 −σ1)2

Эта формула с учетом равенства (1.26) записывается более компактно:

 

1

 

 

 

интp

.

=

 

 

 

 

 

 

2 σинт

При таком значении величины и производных (2.15) ассоциированный закон течения (2.12) принимает вид:

 

 

3

 

 

интp

 

 

p

=

 

(σi −σ0).

(2.16)

i

 

 

 

 

2

 

σинт

Отсюда следует, что компоненты тензора приращений пластических деформаций пропорциональны компонентам (σi −σ0) девиатора напряжений. Согласно закону Гука (I.6.6), упругая доля приращения деформации ε1 определяется равенством:

1(e) = E1 [1 − ν(2 + 3)].

Но σ2 +σ3 = 3σ0 −σ1 и тогда (см. формулу (1.37))

1

=

E (1+ν)1 3νdσ0

 

= 2G 1

1+ν 0

.

(e)

 

1

 

 

 

 

 

1

 

 

 

3ν

 

 

Стало быть,

 

 

 

1

i

 

3ν

.

 

 

 

 

 

 

i(e) =

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

2G

1+ν

 

 

 

754 Часть VII

Полное приращение i деформации εi складывается из величин i(e) и (2.16):

 

1

i

3ν

+

3

 

 

интp

 

i =

 

(σi −σ0).

 

 

0

 

 

 

 

2G

1+ν

2

 

σинт

Отсюда можно перейти к координатной форме представления приращений деформаций. Ниже такая форма приводится для базиса x, y, z, не совпадающего с главными осями:

 

1

 

 

 

3

 

 

интp

 

 

 

 

 

 

x =

 

 

[x − ν(y + z )] +

 

 

 

 

 

(σx −σ0),

E

2

σинт

 

 

 

· · ·

· · ·

· · ·

интp

 

· · ·

 

 

 

 

zx =

1

zx + 3

 

τzx.

 

 

 

 

 

 

 

 

 

 

 

 

G

 

σинт

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.17)

Полученные соотношения являются основными уравнениями теории течения. Для плоской деформации они были предложены Л. Прандтлем в 1924 г. Через 6 лет А. Рейс дал решение задачи в общем случае. При рассмотрении краевой задачи к шести физическим уравнениям (2.17) добавляются три условия равновесия Навье (I.2.6) и шесть геометрических соотношений. Последние берутся либо в форме уравнений Коши (I.5.7), либо в форме условий сплошности Сен-Венана (I.5.9).

Интегрирование полной системы уравнений теории течения возможно лишь при известной истории нагружения. Эта задача чрезвычайно сложна даже тогда, когда решение разыскивается численно. Поэтому желательно выделить те случаи, в которых возможны хотя бы небольшие упрощения. В этом отношении особый интерес представляет собой теория малых упругопластических деформаций, которая входит в число так называемых деформационных теорий пластичности.

2.6. Теория малых упругопластических деформаций. Выше уже говорилось, что величина достигнутой пластической деформации зависит не только от значений напряжений на данный момент времени, но и от всего пути нагружения, и что для каждого конкретного пути можно установить конечную связь между напряжениями и деформациями (т. е. проинтегрировать уравнения (2.17)). В этом случае (при условии, что нагружение активное) упрочняющееся тело аналогично нелинейно-упругому материалу в том отношении, что для обоих тел связь между напряжениями и деформациями взаимно однозначна. Для нелинейно-упругого тела физическим уравнениям может быть придана форма закона Гука с переменными упругими характеристиками, зависящими от деформаций. Замена такого рода конечными

σинт

Глава 2

755

соотношениями физических уравнений (2.17) и составляет основу деформационных теорий пластичности.

Получить искомые конечные соотношения между напряжениями и деформациями можно лишь при условии, что нагружение является простым: отказа от рассмотрения траекторий нагружения, не включающих разгрузку, еще не достаточно. При простом нагружении компоненты девиатора напряжений возрастают пропорционально некоторому параметру α:

σi − σ0 = α(σi − σ0) , 0 ≤ α ≤ 1.

Значком отмечены компоненты напряженного состояния в конце нагружения. Тогда (см. формулу (2.16))

ip = 32 (σi −σ0) pинт.

Следовательно,

ε p =

3

 

(σi −σ0)

 

 

 

p .

(2.18)

 

2

i

 

σинт

 

 

инт

 

После подстановки величины εip , даваемой этим равенством, в формулу (1.33a) для интенсивности пластических деформаций и выполнения ряда преобразований с использованием выражения (1.26) получится, что

εip = интp . (2.19)

Но тогда pинт = интp . Стало быть, в случае простого нагружения интенсивность интp приращения пластических деформаций равна приращению инт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

[σ

 

ν(σ

+σ

+σ

)]

12ν

σ

 

=

1

(σ

1

σ

).

 

 

E

 

E

 

2G

 

1

0

 

 

1

2

 

3

1

 

 

0

 

 

0

 

Таким образом,

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εie =

 

(σi

−σ0) + ε0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2G

 

 

 

 

 

 

 

+ εинтp ,

756

Часть VII

Эту упругую составляющую полной деформации εie с учетом зависимости σинт = 3инт, даваемой формулой (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 12ν σ ,

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

(N)2 +(M)2/3

(2.27)

Так как физические уравнения теорий течения и малых упругопластических деформаций содержат интенсивности напряжений, деформаций и приращений пластических деформаций, то без помощи выражений (2.26) и (2.27) задачу о распределении напряжений в упругопластическом теле не решить.

Согласно условию (2.4) начала текучести Губера – Мизеса, а именно оно принимается при решении рассматриваемой задачи, материал переходит в пластическое состояние, если

σинт = σT или τинт = τT.

(2.28)

Из формул (2.26) и (2.28) следует, что

 

σTT = σинтинт =

 

.

 

3

 

Значит,

 

σT =

 

τT,

(2.29)

3

σN2 + 3τM2 = σT2.

(2.26a)

Для несжимаемого материала ν = 1/2, а в этом случае, согласно формуле (1.37), E = 3G. Тогда в момент наступления текучести

σT = T, τT = T = 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,

 

 

= ()2 + ()2.

(2.33)

εT

σинт

 

 

 

εинт

 

 

 

 

интp

 

 

 

 

Теперь можно обратиться к непосредственному решению задачи. Сначала оно выполняется средствами теории течения.

2.7.1. Решение на основе теории течения. Из физических уравнений (2.17) этой теории нетривиальными остаются только два уравнения: одно для нормальных напряжений, другое – для касательных. При этом упругие составляющие деформаций отбрасываются (идеально жесткое пластическое тело):

 

 

 

3

 

интp

 

 

 

 

 

 

 

 

 

 

интp

 

 

 

 

N =

(σN −σ0), dγM = 3

τ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γ = τ

 

 

.

 

()2 + ()2

()2 + ()2

(2.34a)

Отсюда следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

τ =

 

dγ/dε

 

 

 

(2.35)

σ =

 

 

,

 

 

.

 

 

 

 

1+(dγ/dε)2

1+(dγ/dε)2