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

книги из ГПНТБ / Прессование алюминиевых сплавов. Математическое моделирование и оптимизация

.pdf
Скачиваний:
15
Добавлен:
25.10.2023
Размер:
13.98 Mб
Скачать

Кроме того, имеет место равенство тепловых потоков:

■м

■= — кс і—

(11-65)

 

' дп

s

'• дп

 

Уравнение состояния

 

Описывая

поведение реальных

металлов и сплавов

в процессах обработки металлов давлением, воспользу­ емся реологическими уравнениями [12].

Как обычно, предположим, что компоненты скорости деформации hk могут быть представлены в виде суммы

двух составляющих:

 

= +

66)

Упругие составляющие l eik линейно связаны с компо­

нентами скоростей напряжений от'-

\ik — Aiklm Olm .

Тензор упругих коэффициентов Ашт, очевидно, дол­ жен удовлетворять условиям симметрии по индексам г, k; I, т:

^ i k l m

^ k i l m

^ k i m l -

Вязко-пластические составляющие скоростей дефор­ маций ^vik удовлетворяют условию несжимаемости:

и связаны следующими соотношениями с компонентами девиатора напряжений:

№ 6?)

Интенсивность касательных напряжений Т=

/ ' sik sik есть некоторая функция температуры и, интенсивности скоростей вязко-пластических деформа­

ций Нг,= и интенсивности остаточных дефор­ маций сдвига:

t

Г11= fj№dt.

70

Таким образом, мы имеем уравнение состояния:

Т= Т(и, Н®, Г®).

Вдальнейшем, если противное не оговорено особо, будем предполагать упругие составляющие пренебрежи­ мо малыми по сравнению с вязко-пластическими, опу­ ская при этом значок ѵ в формулах, и получим уравне­ ние несжимаемости:

І и = div V = 0;

уравнение связи:

Е _ J L с .

ік’

уравнение состояния: Т = Т (и, И, Г).

Поскольку для несжимаемой среды

<*ik bk = Sik lik = ri 2l tk lik = TH,

уравнение энергии принимает следующий вид:

1a^vdS +

И Г pFt vt dW = (Т \ TUdW

 

D

1

dv2 ,w/

— p ----dW.

2

dt

(11-68)

(II-69)

(11-70)

Аналогично

 

®ik k 'Tк h

2|/*ö£< = T6H,

 

II

и вариационное уравнение (ІІ-51) можно записать в сле­ дующей форме:

 

a"8vdS +

f j j pF, öv, dW = f j j ' TöUdW

 

 

 

 

 

(ІІ-71)

 

Наконец, уравнение теплопроводности (II-57) сво­

дится к следующему:

 

 

du

( д2и

д2и

д2и \

(Н-72)

dt

X I — г

Н----- “

Н--------- I d- ѵТН.

дх7

дхі

дх.

 

71

Вероятность разрушения [13].

В теории разрушения металлов, развитой В. Л. Кол­ могоровым [14], предполагается, что можно ввести в

рассмотрение некоторую скалярную величину — «трещи­

новатость» и, которая характеризует пораженность эле­ ментарного объема деформируемого тела микродефек­ тами. _

Приращение to за время dt принимается пропорцио­ нальным приращению степени пластической деформации dT = Hdt и обратно пропорциональным пластичности Гр металла при реализуемом в данный момент времени на­ пряженном состоянии, показатель которого k = ocр/Т.

Моменту разрушения соответствует некоторая трещи­ новатость ар. Обозначив

(II-73)

приходим к условию разрушения м = 1 .

Воспользуемся методами математической теории на­

дежности [15]

для

расчета в е р о я т н о с т и

разруше­

ния металлов

при

пластической деформации.

Рассмот­

рим испытания на пластичность, проводимые при фикси­ рованной температуре и, скорости деформации Н, пока­ зателе напряженного состояния k.

Пусть образец разрушается при Г= у- В общем слу­

чае у — случайная величина с законом

распределения

Q ( T ) ^ P { y > T } .

(II-74)

Функция <2(Г) есть вероятность разрушения образ­ ца до деформации Г (рис. 8). Будем предполагать, что функция <2(Г) непрерывна и существует непрерывная плотность вероятности разрушения д(Г) = <2/(Г).

Наряду с функцией ф(Г) будем употреблять и дру­ гую функцию:

Р(Г) = 1 — <2(Г) = Р { у > Г } —

(И-75)

функцию надежности.

Рассмотрим некоторые числовые величины, характе­ ризующие функции Р и Q. Важнейшей из них является математическое ожидание деформации разрушения Му, совпадающее с пластичностью металла Гр при детерми­ нированном подходе к проблеме разрушения, развитом в работе [14].

72

Более логично считать, что

оо

 

Гр = Му = J уq (у) dy.

(11-76)

Другой характристикой

распределения вероятностей

является дисперсия деформации разрушения:

оо

Dy = М (Г — Гр) = Му2 — (Му)2= f у2 q (у) dy — Г*.

О

Рис. 8. Кривые распределения вероятности разрушения для различных зна­

чений вектора s (k, и, Н...)

Введем теперь широко используемую в теории на­ дежности «функцию опасности разрушения».

Поставим следующую задачу. Пусть образец не раз­ рушился до деформации Г. Какова вероятность того, что он разрушится на участке Г, Г+ДГ.

Обозначим эту вероятность через Q(T, Г+ДГ).

Имеем

 

 

Q(r, Г +

Д Г )= — у |^ -Д Г + 0(ДГ),

(Н-77)

где О (ДГ)— бесконечно мйлая величина

высшего по­

рядка. Введем обозначение

 

Я ( Г ) = —

Р(Г)

(И-78)

'

 

73

Тогда

 

 

г

Mr>dr

 

—j

 

Q (Г) = 1 — (1 — Qo) e *•

.

(II-79)

Эксперименты показывают [16, 17], что, как прави­ ло, имеет место нормальный закон распределения веро­ ятностей Q(F). Докажем это аналитически, следуя [15].

Предположим, что образец, вырезанный из слитка, прокатанной, прессованной или кованой заготовки, име­

ет «начальную трещиноватость» ©о, распределенную по

Рис. 9. Нормальный закон распределе-

Рис. 10. Функция опасности разру*

ния деформации разрушения V

шения

нормальному закону с небольшой дисперсией. Пусть при

нагружении образца параметр © изменяется детермини­ рованным образом (рис. 9):

© = / (Г, ©о),

причем функция © монотонна по Г.

По условию разрушение наступит тогда, когда ю

превысит некоторый критический уровень ©р. Значит, деформация разрушения определяется из уравнения

f (У, «о) =

НЛП

У = Ф («0, (Dp),

где ф — обратная к f по первому аргументу функция.

74 .

Поскольку coo имеет малую дисперсию, то, раскладывая

функцию ф по формуле Тейлора в ряд в точке a— M(nQ, где М —математическое ожидание, и пренебрегая чле­ нами второго порядка, получим:

у = ф (а, (Dp) + фа (а, ыр)((0о — а).

Но если о)о является нормальной величиной, то де­

формация разрушения у как линейная функция от wo также распределена по нормальному закону. График функции опасности разрушения при этом показан на рис. 10.

Перейдем к анализу вероятности разрушения элемен­ тарного объема в условиях сложного нагружения, когда вектор

s{k, и, II...),

компоненты которого характеризует весь комплекс па­ раметров (вид напряженного состояния, температуру, скорость деформации и др.) для произвольного момента времени t, не остается постоянным.

Воспользуемся формулой (11-78), записав:

d l n P = — Ks [ Tf ( o) Tfd ta .

Здесь Xs — функция опасности разрушения, Гр —

математическое ожидание деформации разрушения, вы­ численные для реализуемого в данный момент времени

вектора s;

 

 

 

 

< fe )= ^ =

H dt

 

 

(II-80)

Г?

r

sp

 

 

 

Математическое ожидание величины со, соответству­

ющей разрушению, равно единице.

 

 

Интегрируя (П-80), получаем

 

 

 

 

t

 

 

ln j r =

-

(r > , r f *

= - J Ы 1»

H«*.

(H-Sl)

 

co0

 

U

 

 

поэтому вероятность разрушения равна:

 

 

 

 

t

 

 

 

Q = 1 (1

- 1 М

гр“) ІШ

 

(11-82)

Qo)е

 

 

75

Легко видеть, что при монотонной деформации (s =

—const) из (П-82) будут получены распределения веро­ ятностей <3(Г), соответствующие испытаниям образцов на разрушение.

В общем случае формула (П-82) позволяет рассчи­ тать вероятность разрушения элементарного объема для немонотонных процессов.

6. МЕТОД КОНФОРМНЫХ ОТОБРАЖЕНИЙ В ПЛОСКИХ И ОСЕСИММЕТРИЧНЫХ ЗАДАЧАХ ТЕОРИИ ПРЕССОВАНИЯ

Плоское течение

Изучая решения, полученные в теоретической гидро­ механике [9], легко убедиться в том, что имеется глу­ бокая аналогия между движением идеальной жидкости и течением металла при обработке давлением. Можно ожидать в связи с этим, что методы гидродинамики при­ годны для построения приближенного «опорного» реше­ ния задач о пластическом течении.

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

Рассмотрим плоскопараллельное течение упрочняю­ щегося несжимаемого материала в области D — криво­ линейной полосе с границей 5 (рис. 11). Вместо коор­ динат Хі (( = 1, 2, 3) введем более удобные для плоско­ го случая переменные х = х ь у = х 2, причем комплексная

переменная

z — x-{-iy будет

обозначать

точки физиче­

ской плоскости,

а

ѵх = ѵх, ѵ2= ѵ у,

п3= 0 ,

£и = £хж, \22 =

~ і г / у I Іі2 =

5жг/і

0'11 =

ОХж, n 22 =

0' y y ,

G \ 2 z = O x y -

Основные уравнения в новых обозначениях принима­ ют следующий вид:

а) уравнения движения

76

д®хх

I доХу _

р dvx _

 

дх

ду

dt

(11-83)

д^ху_|_ дОуу __

.

 

дх

ду

dt

 

б)

условие несжимаемости

 

dvx I

дОу _

Q

(11-84)

d*

öl/

 

 

 

Рис. 11. Течение материала в м н о г о у г о л ь н о й области (прессование полосы)

77

в)

 

уравнения связи напряжения и скоростей дефо

маций:

 

 

 

 

^xx

"ср

2Т_ Щц '

 

II

дх

 

 

 

 

 

а ср

дѵу

 

(ІІ-85)

Ü !/iJ

Н ~ду

 

 

 

° х у =

_т_ /

дух .

 

 

 

II [

ду

 

 

 

где

 

 

 

 

 

т = | / + ^ ог«)

H= i/(¥+ f+W+¥j^

<іш>

2

 

®ср=

(11-88)

Будем предполагать шока, что упругие составляющие деформаций пренебрежимо малы по сравнению с пла­ стическими составляющими.

Для идеально пластического материала

Т =

T S =

const.

 

 

(ІІ-89)

 

В общем случае

 

 

Т =

Т(и,Н ,Г),

 

 

(11-90)

где

и — температура.

 

 

 

t

 

 

 

 

 

r =

jlM ;

 

 

 

(П-91)

 

о

 

 

 

 

 

 

г) уравнение

теплопроводности

 

 

 

 

 

 

 

да-92)

где х и V — константы;

 

du

ди

.

ди

,

ди

(ІІ-93)

dt

= --------

1--------

V .

-j---------

ду

dt

 

дх

'

 

Температурными

напряжениями для несжимаемой

среды также будем пренебрегать.

 

78

Обозначим компоненты поверхностных напряжений

оп, действующих на единицу площади, как опх и апу, имеем уравнение энергий в следующей форме:

f ( О п х ^ + OnyVy) d S = ^ T l l d x d y +

f f p

^ +

S

D

D

 

+

^ - v ^ d x d y .

 

(Н-94)

Обозначая р и —т, ѵп и ѵх соответственно нормаль­ ную и касательную составляющие поверхностной нагруз­

ки о п и скорости V, имеем

 

f (<Ѵ Ox + ояу Vy) dS = Â — Â T ,

 

s

 

 

где

 

 

 = I*pvn dS, Ax =

I w x dS,

(11-95)

s

 

 

так что

 

 

А — IJ THdx dy +

I TÜT dS -f

(11-96)

'D

S

D

Предположим, что граница S может быть разбита на две части: S b на которой заданы ѵп и т, и свободная по­ верхность S2, на которой р = т= 0 .

Тогда принцип виртуальной работы [уравнение (П-51)] можно записать в следующей форме:

JI* ТбНсД dy + J%ЬѵхdS + Jj p

8vx +

D

Si

D

 

+

^Ü Vy'jdxdy = 0.

(11-97)

 

Для идеально

пластичного

материала (T = Ts), ког­

да инерционными силами можно пренебречь, уравнения

(II-96) и (11-97) принимают следующий вид:

 

А = T S f f Нdx dy

%vxdS rnin.

(11-98)

D

 

 

Следует заметить, что непосредственное применение метода Ритца для решения описанной задачи связано

79

Соседние файлы в папке книги из ГПНТБ