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

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

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

Коэффициентом Ck обозначим отношение скорости истечения через канал ©ь к средней скорости истечения:

ücp

(II-184)

 

Легко видеть, что

при прессовании в один канал

(п = 1 ) с\ = 1. В общем

случае эта величина отлична от

единицы и определяется соотношением между потоком металла, вытекающим из канала со&, и площадью этого канала:

Q k l _

Qk_

(11-185)

Фо

“ft

 

ѵо ~Z~

 

 

Для определения параметров

воспользуемся ва­

риационными принципами механики сплошных сред.

Зададим

кинематически возможное поле скоростей

и, удовлетворяющее уравнению несжимаемости divo = = 0, кинематическим граничным условиям и описываю­ щее течение металла при прессовании.

Распределение потоков будем искать из вариацион­

ного уравнения

 

 

Т(и,Н,Г) 8HdW -h jfx 6 v TdS=0

(11-186)

при дополнительном условии

 

 

+ <72 -------1~<7п1 •

 

 

Здесь Т —интенсивность

касательных

напряжений;

Н —интенсивность скоростей

деформации

сдвига; т —

касательные силы трения на боковой поверхности кон­

тейнера и

матрицы;

ѵѵ — касательная

составляющая

скорости;

W — объем

области течения;

5 — контактная

поверхность; и — температура.

В монографии [19] эта задача была решена с приме­ нением разрывных полей скоростей для идеально пла­ стического материала.

Использование разрывных полей скоростей, а также идеализация свойств материала (жестко-пластический), очевидно, не приводят к существенным погрешностям при анализе энергосиловых параметров процесса.

В то же время можно ожидать, что сложные свой­ ства реальных металлов и сплавов, в частности упроч­

100

нение, вязкость, а также неоднородность температурно­ го поля, могут существенно повлиять на кинематические условия процесса, исследование которых при прессова­

нии профилей сложной формы или при

многоочковом

прессовании является задачей особой важности.

Для более полного анализа откажемся

от гипотезы

о существовании разрывов и построим

непрерывное

X

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

кинематически возможное поле скоростей. Воспользу­ емся методом суперпозиции точечных источников и сто­ ков [20].

В работе [21] приводится решение задачи о неста­ ционарном температурном поле бесконечного цилиндра, внутри которого в точке с координатами р', Ѳ', z '= 0 расположен мгновенный тепловой источник производи­

тельности 2Q (рис. 16).

Это решение для случая тепловой изоляции боковой

поверхности цилиндра

(ди/дп= 0) после перехода к ста­

ционарному режиму принимает следующий вид:

и = 2Qf (р, Ѳ, z, р'.Ѳ'),

(IM 87)

где функция f записывается в виде двойного ряда:

Уcos т (Ѳ — Ѳ') X

101

 

fe|*° Jm[ h

£ jS „

 

 

(II-188)

 

(ßl + ^ 2) sm (ßk)

Здесь

Jm— функция

Бесселя первого рода m-го поряд­

ка; ßft — корни трансцендентного уравнения

 

=

(ІМ 89>

Из

симметрии решения относительно плоскости г —

= 0 следует, что условие ди/дп= 0 выполняется также на плоскости 2= 0.

Очевидно, фундаментальная функция

Ф = Q f ( P ,0 ,z ,p ',0 ')

(II-190)

для области 2> 0 удовлетворяет уравнению

Д<р = 0,

(ІИ 91)

на границе области течения

(поверхности цилиндра

и плоскости 2= 0) условию

 

^ = 0

(II-192)

дп

 

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

денного

стоком

производительности Q,

находящимся

в точке ЛІ(р', Ѳ', 2= 0) плоскости 2= 0.

прессования

В рассматриваемой задаче теории

в каждой точке Р области

 

(О= (Ö! +

С02 -1-------ь С0„

 

скорость истечения металла равна

 

ѵ(Р)= К>

если

/>£<»)*

(Н-193)

 

(0,

если

P £a)ft.

 

Это позволяет построить кинематически возможное безвихревое поле скоростей, порожденное системой сто­ ков, распределенных в области ю с удельной производи­ тельностью ц(р', Ѳ'). Потенциал такого течения, соглас­ но принципу суперпозиции, может быть получен интег-, рированием фундаментального решения по области со:

102

для каждой итерации, считаем распределение Т/2Н за­ данным предыдущим решением, т. е. решаем задачу о течении вязкой среды с переменным по объему коэф­ фициентом вязкости.
После выполнения условий сходимости уточняется температурное поле, делается новый расчет поля ско­ ростей и т. д.
Исследуем полученное решение. С этой целью из на­ чала координат на плоскости матрицы опишем окруж­ ность минимального радиуса г0, включающую в себя
103
Т (и, Н, Г)
(ІИ 95)
Т (и, Н, Г) gjj2
со ско­
J J wk
ветствующий истечению металла из канала ростью, равной единице.
Итак, потенциал поля скоростей представлен в виде линейной функции варьируемых параметров qh. Анало­ гичный вид будут иметь составляющие скорости, компо­ ненты скоростей деформаций.
Используя модифицированный метод Ритца [22] и решая задачу методом итераций, на каждом этапе решения будем иметь линейную систему уравнений от­ носительно параметров q^.
Практическая реализация описанного алгоритма сводится к следующему. В области течения И? и на гра­ нице 5 строится сетка, в узлах которой вычисляются и запоминаются значения потенциалов
Соответствующие значения компонентов вектора скорости и скорости деформаций находятся методом конечных разностей и также записываются в память ЭВМ.
Приняв в первом приближении скорости истечения vh равными средней скорости (сь=1), решаем уравне­ ние теплопроводности.
Далее осуществляется переход к блоку модифици­ рованного метода Ритца. Используя соотношения
где ф^=
f(p, Ѳ, z, p', Ѳ')р', dp'dQ' — потенциал, соот­
(IM 94)
Cd
«Po (p.ѳ**) -
JJ 0 (p'tѳ') / (p. Ѳ» 2. P', ö') p' dp' ctö'

область со. При выполнении условия r0I R o < .< L \ на до­ статочном удалении от стенок контейнера, мы приходим к задаче об истечении металла из полупространства 2 > 0 , ограниченного плоскостью z = 0, через каналы со*0

(см. рис. 17).

Наряду с цилиндрической системой координат будем использовать декартову систему координат (х, у, z), расположенную, как указано, на рис. 17.

Роль фундаментального решения в этом случае вы­

полняет

потенциал

стока, находящегося

в точке

с ко­

ординатами

(g, г|)

плоскости 2 = 0 . В точке М(х,

у, г)

этот потенциал равен:

 

 

ф = ------.

-

.. Q

---------- .

(II-196)

2яѴ (ж — £)* + (у — л)* + **

 

 

Прессование круглого прутка

В частности, если сток находится в начале коорди­ нат ( |= г | = 0), потенциал равен:

Ф = ..... .

(II-197)

2я у х2+

у2+ z2

104

Тем самым, переходя к сферической системе коор­ динат (г, %, Ѳ), получаем потенциал радиального тече-' ния металла:

Ф =

_0_

(II-198)

2лг

Составляющие скорости равны:

ѵ' = — £ г , '

ѵ%= ѵе - ° -

(И'«!»)

Компоненты

скорости

деформа­

ции находятся по формулам:

-

_ дѵг _ Q

 

 

г

дг

я г3

 

Н = — =

^

-

(11-200)

х

г

2яг3’

t

ІѴ_ __ _Q__

 

г

2яг3'

 

%с = \ ѳ = Ѵ = °- j

Интенсивность скоростей дефор­ мации сдвига равна:

н = 2 К е + Е Л + 3

Кз(з

я/-3

 

(11-201)

Рис. 18. Прессование прутка из круглого контейнера

Рассмотрим в качестве примера разрывное решение задачи о прессовании круглого профиля (рис. 18), ре­ зультаты которого используем в дальнейшем.

Описанное выше поле скоростей соответствует ради­ альному течению металла в области r0<.r<.R0— в по­ лой полусфере, наружный радиус которой равен радиу­ су контейнера, а внутренний — радиусу прессуемого прутка. Вне этой области металл перемещается как аб­ солютное жесткое тело.

Подсчитаем мощность внутренних сил, выделяемую в этой области при прессовании идеально пластическо­ го материала:

105

 

 

Я„Я/22я

/-о —

 

Л0 = JJJTs HdW = I

j* I*

— - y -■r2 sin %drd%dQ =

 

r

 

n

0 o

 

 

= Qaf In (&.)*.

 

 

(II-202)

Учитывая, что из условия несжимаемости

 

ѵ0R20 =

у2яг2,

 

 

 

получаем

 

 

 

 

Ä0 — Qosln — .

 

 

(11-203)

 

 

°0

 

 

 

Вычислим мощность, развиваемую силами трения

pGs на поверхности матрицы:

 

л»

 

 

 

(11-204)

Ах = j

j* pos |or| rdrdQ = — Qpos In ~ .

r0

0

 

 

 

 

Общая мощность составляет:

 

A = Ло +

Ä = Q pa(l

+ ^ -) l n ^ - .

(II-205)

т. е. при фиксированных свойствах материала и услови­ ях трения она пропорциональна потоку Q и логарифму

отношения скоростей In ѵ/ѵ0.

Следует отметить, что

в этом в общем известном [5]

решении не учитывается

дополнительная мощность, рассеиваемая на поверхно­ сти разрыва.

Рассмотрим теперь течение металла в несколько круглых каналов соь, расположенных на плоскости

г — 0.

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

Q_

Фо —

где

Я і

,

Я а

,

+ Яп_

(И-206)

ж

+

к

+

Rn

 

R k = V ( x - l k f

+

{ y - 4 \ k ? +

z \

в свою очередь §*, щ — координаты центра тяжести ка-

НЯ Л Я (і)£.

106

На рис. 19 методом Максвелла построены линии то­ ка на плоскости £ = 0 для некоторых случаев течения. Для сравнения на рис. 20 приведены траектории частиц

â

Рис. 19. Картина течения на плоскости 2=0:

a — q, = q2; 6 ~ q , - 2 q 2

металла, движущегося по поверхности прессовой матри­

цы {qx = q2).

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

107

 

_1_

0(1.11) dUr\

 

 

 

 

 

Фо

ЯV (x - ly + ü - W + z*

 

 

 

(П-207) -

 

 

где

 

 

 

 

 

 

ѵ(Р)

=

vk, если

Р £(ofc

 

 

О, если

Р £ (дк.

 

 

 

 

 

 

 

 

Аналогично

(II-194)

 

 

получаем:

 

 

 

 

Фо

=

»CP .ai Фі +

 

 

 

f

Фг +

••• + ^ - ф я

 

 

а2

 

 

 

а„

 

 

Рис. 20. Картина течения металла

 

 

на

поверхности

матрицы (qt= q 2)

где

Ф * = *

________ßdn_________

(Н-208)

 

V{х-І)* + (у-ц)* + г*

Вкачестве иллюстрации решим задачу о прессова­ нии тонкой полосы шириной 2 а и толщиной 2 6. В этом случае п — 1, <7і= а і= 1 -

сіУц

Фо

8яай

l)2+ (j/-ri)2 + z2

 

 

 

 

 

 

 

 

 

 

—а —b

 

 

Ввиду соотношения 6<Са можно приближенно запи­

сать:

 

 

 

Ф = 4іш

dl

 

1)2 + У2 + г*

 

 

 

 

 

 

Q

in * + а + У(х + а)2 + У2+ г2

(11-209)

 

 

х — а + Ѵ (х — а)2 + J/a + г2

 

 

 

108

Асимптотическое разложениё

 

Найдем

асимптотическое разложение потенциала

Ф о . С этой

целью в области z > 0 , заполненной

прессуе­

мым металлом, построим полусферу радиусом

г0с цент­

ром в начале координат. Как мы уже указывали, вели­

чина

радиуса выбрана так, чтобы область (Ü=

CÜI+ G)2+

+...+<ön целиком лежала внутри окружности L.

Подынтегральная

 

функция в выражении

(П-207)

может быть представлена в виде

 

________ ____________

т Ѵ а 2 — 2acosy-[-l

(II-210)

V ( X -

1)2 + - Г))* +

 

Z2

 

где

 

 

 

 

а =

cos у

хі + уп

 

 

 

 

 

Рг

г■—радиус-вектор точки Р (|, ц );

ррадиус-вектор точки М(х, у, z).

Вне полусферы

выполняется неравенство

а < 1

и справедливо следующее разложение [23]:

 

V а 2 — 2а cos у + 1

т£— О ат Рт(cos у).

(11-211)

Здесь Pm(cos у ) — шаровые функции (полиномы Ла­ гранжа) :

Р0 (cos у) =

1;

 

Рх(cosy) =

cosy;

(II-212)

Р2 (cos у)= ——(3 cos2 у — 1)

 

и т. д.

Подставив эти соотношения в (П-207), после инте­ грирования получаем:

 

/

ч

Q

f,

, хІс + УЧс

I

 

 

= - —

 

+

+

- 1-

[л- (2 7 ,— 7 (;)

— Д ) +

 

zp

 

 

 

 

 

+

6xyJxy -

z2(Jx +

7У)} + • • •} •

(Н-213)

109

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