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

книги / Механика сплошных сред (теоретические основы обработки давлением композитных материалов с задачами и решениями, примерами и упражнениями)

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
9.02 Mб
Скачать

 

3.1. ПРИМЕНЕНИЕ РАЗРЫВНЫХ КВ ПОЛЕЙ СКОРОСТЕЙ

где K = 1+ tg2 α

; α – угол захвата металла валками. Учитывая, что элемент

2

 

контактной поверхности d = KdE2 , а контактные напряжения трения

τот = const и τоп = const, находим значения первого слагаемого баланса мощности (1.3.22)

0

 

τ

 

 

h K 2

 

h

τопV

 

d

= τопV0

0

 

ln

1

 

2tg

α

hн

н

 

 

 

 

 

 

 

 

 

 

 

2

 

 

и четвертого слагаемого этого баланса

н

 

τ

 

 

h K 2

 

h

τотV

 

d

= −τот V0

0

 

ln

н

,

 

 

α

 

д

 

 

 

 

2tg

 

h0

 

 

 

 

 

 

2

 

 

 

где hн – высота проката в нейтральном сечении (Е2 = – н ). Отсюда сумма этих слагаемых при условии, что τот = τоп = τn, имеет вид

 

0

 

 

н

 

 

 

 

 

h K 2

 

h h

τn

τ V τd

τ

 

V τd

 

= −τn V

0

 

 

0 1

.

 

 

α

2

 

оп

 

 

от

 

 

0

 

 

н

 

 

д

 

 

 

 

 

2tg

2

 

hн

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переходим к расчету мощности внутренних сил Int с учетом (3.1.19):

 

 

 

 

α 0

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

4 E2

α

 

 

 

 

Int = 4τ

 

v h tg

 

 

 

 

 

 

1

+

1

tg2

 

 

dE dE

 

=

т

 

 

h2

 

 

 

2

 

0 0

 

2

 

 

 

h2

2

 

 

1

 

 

 

 

 

 

д

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 tg

α

 

 

 

 

 

 

 

 

 

 

 

= 2τ

v

h

 

0

 

2

1 + u2 dudE ,

 

 

 

 

 

т

0

0

h

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

0

 

 

 

 

 

 

 

 

 

371

3. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ОМД МЕТОДАМИ МСС

гдепроизведеназаменаинтегрированияпоE1 интегрированиемпоu: 2 E1 tg α2 h = u .

Тогда после интегрирования по u имеем

 

α

 

α

 

0

 

dE

2

 

Int = τт v0 h0 K tg

 

+ ln tg

 

+ K

 

 

.

2

2

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

 

При интегрировании по E2 выполним замену dE2 = −0,5dhtg α2 c помощью (3.1.3). Тогда

 

 

 

 

K tg

α

 

 

α

 

 

 

 

 

 

+ ln tg

 

+ K

Int = τ

т

v

h

 

2

 

 

2

 

ln

h0

.

 

 

 

 

 

 

 

 

0

0

 

 

2tg

α

 

 

 

h1

 

 

 

 

 

 

2

 

 

 

 

 

Окончательно из баланса мощности (3.1.23) получаем значение контактного касательного напряжения при плоской листовой прокатке:

 

 

α

 

α

 

h

+ tg2

α

 

α

(σзадн − σпер )

 

τт K tg

 

+ ln tg

 

+ K ln

0

 

 

+ tg

 

2

2

h1

2

2

τn =

 

 

 

 

 

 

 

 

. (3.1.23)

 

 

 

 

K 2 ln

h0 h1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hн2

В этой формуле высота нейтрального сечения hн определяется из формулы (3.1.3) при Е2 = – н . Положение нейтрального сечения н в очаге деформации определим из уравнения равновесия всех действующих на прокатываемый металл поверхностных сил с учетом принятого условия τот = τоп = τn:

 

 

τn ( д

н )

 

 

 

 

 

 

pn

 

 

 

 

 

 

 

 

 

2

 

τn н

cos α

− 2

д

sin

α h

σ

 

+ h σ

 

= 0.

 

 

 

 

 

 

 

задн

пер

 

 

α

 

 

 

α

 

2

 

 

α

2

0

 

1

 

 

 

cos

 

 

cos

 

 

cos

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

372

3.1. ПРИМЕНЕНИЕ РАЗРЫВНЫХ КВ ПОЛЕЙ СКОРОСТЕЙ

Отсюда

 

1

 

 

n

p

n

 

α

д h0

σзадн + h1

 

 

 

н =

 

 

τ

 

 

tg

2

 

σпер

,

(3.1.24)

 

 

 

 

2τn

 

 

 

 

 

 

 

 

 

 

 

где с учетом абсолютного обжатия h = h0 h1 геометрическая длина очага деформации (рис. 92)

 

 

 

 

д =

h

=

R

h

 

 

h2

.

 

 

(3.1.25)

 

 

 

 

2tg

α

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Подстановкой (3.1.24) в (3.1.3) при E2 = – н

и h = hн находим высоту нейт-

рального сечения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

= h

+

tg 2

 

τn pn tg

α

 

h

 

σ

 

 

+ h

σ

.

 

 

 

д

 

задн

н

1

 

τn

 

 

 

 

0

 

1

 

пер

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

Применяемые обычно в теории ОМД законы трения Э. Зибеля (2.2.29)

τn = 2 μт τт

и Т. Кармана (2.2.30)

τn = μp pn

позволяют записать высоту нейтрального сечения в виде

 

 

α

 

1

 

α

 

 

 

h σ

задн

h σ

пер

 

 

h

= h + tg

1

tg

 

 

0

1

.

(3.1.26)

2

μp

2

 

 

4μтτт

 

н

1

 

 

 

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

μт

2ε

 

μp =

μт mχ

 

 

;

 

,

(3.1.27)

(2 − ε)m

1+ mμт

373

3. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ОМД МЕТОДАМИ МСС

где относительное обжатие

ε =

 

h

;

(3.1.28)

h

 

параметр прокатки

0

 

 

 

 

 

 

 

 

 

m =

д

;

(3.1.29)

h

 

 

 

 

 

 

ср

 

 

средняя высота проката hср = 0,5(h0 + h1); коэффициент, учитывающий влия ние внешних зон, – χ = 0,4 при m < 1 и χ = 0 при m ≥ 1.

Формула (3.1.23) позволяет рассчитать момент прокатки М, необходимый для вращения рабочих валков радиуса R. При известной средней ширине про ката bср. момент прокатки

М = τnR2bср(α – γ),

(3.1.30)

где нейтральный угол определяется по формуле С. Экелунда с учетом заднего и переднего натяжений:

sinγ =

н

=

sinα

1 cosα

+

h0

σзадн h1

σпер

.

 

2

4μp

 

8 Rμтτт

 

 

R

 

 

 

 

 

Центральная сила P, действующая на валок со стороны прокатываемого ме талла, с учетом (2.2.30) и (3.1.27) также определяется контактным касательным напряжением (3.1.23)

P = pn Rb α =

 

τn

 

 

 

Rb α,

(3.1.31)

 

 

ср

 

 

 

ср

 

 

 

μp

 

где угол захвата определяется соотношением (3.1.4) при замене H на

д :

α = 2arctg

 

h

.

(3.1.32)

 

 

 

2 д

 

Самостоятельно. Используя кинематические параметры разрывного поля скоростей (3.1.2) при изменении текущей высоты проката h по окружности (3.1.6), по компонентам вектора скорости (3.1.7) найти компоненты тензора ско ростей деформаций (3.1.8), где

374

3.1. ПРИМЕНЕНИЕ РАЗРЫВНЫХ КВ ПОЛЕЙ СКОРОСТЕЙ

h′ =

2 E2

;

 

(3.1.33)

R2 E22

 

 

 

 

 

 

h′′ =

 

2 R2

 

;

(3.1.34)

 

(R2 E22 )3

интенсивность сдвиговых скоростей деформаций (3.1.9)

 

h

h2

 

E2

 

2 h2 2

H = 2V

 

0

+

1

h′′−

 

 

 

2

 

 

0

h

 

4

 

h

 

 

 

 

 

 

и определить контактное касательное напряжение τn.

Задача 3.1.2.2. Используя интенсивность сдвиговых скоростей деформаций (3.1.19) задачи 3.1.2.1, определить приращение температуры, вызванное деформационным разогревом при прокатке несжимаемого, идеального жесткопластичного материала.

Решение. В уравнении теплопроводности (1.4.62) слагаемое, характеризующее изменение температуры θ вследствие деформационного разогрева, имеет вид

d θ

=

Tσ Tξ

,

dt

cρJм

 

 

где c – удельная теплоемкость; ρ – плотность деформируемой сплошной среды; Jм – механический эквивалент тепла.

Самостоятельно. Показать, что для несжимаемых сред скалярное произведение тензора напряжений Тσ на тензор скоростей деформаций Тξ в (3.1.25) равно произведению интенсивности касательных напряжений Т на интенсивность сдвиговых скоростей деформаций Н: Тσ Тξ =ТН.

Далее, учитывая, что для идеальных жесткопластичных материалов Т = τт = const, при постоянных значениях плотности и удельной теплоемкости из (3.1.25) получаем

τt

θ= θ0 + cρтJм t0 Ηdt ,

375

3. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ОМД МЕТОДАМИ МСС

где θ0 – начальная температура деформируемого тела. В соответствии с формулой (1.2.162) здесь интегральное выражение точно представляет собой степень деформации сдвига Λ. Поэтому деформационный разогрев сплошной несжимаемой, идеальной жесткопластичной среды при постоянных теплофизических параметрах оценивается по формуле

θ = θ0 +

Λτт

.

(3.1.35)

 

 

cρJм

 

Для вычисления Λ в стационарных процессах воспользуемся формулой (1.2.163). Тогда, применяя формулы для V2 из (3.1.5) и (3.1.19), имеем

 

α 0

1

 

4 E

2

 

α

 

Λ = 4tg

2

 

 

1+

1

tg2

 

dE ,

 

 

 

 

 

 

h

 

h2

 

 

2

2

 

 

 

 

 

 

 

 

д

 

 

 

 

 

 

где принимаем во внимание соотношение (3.1.3) между текущей высотой про-

α

катываемой полосы h и эйлеровой координатой E2: dE2 = –dh/2tg 2 . Тогда

h1

1

 

4 E

2

 

α

h0

1

 

4 E

2

 

α

 

Λ = −2

 

1+

 

1

tg2

 

dh = 2

 

1+

 

1

tg2

 

dh.

h

h

2

2

h

h

2

2

h

 

 

 

 

h

 

 

 

 

 

0

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

В связи с тем, что при h = h1 уровень аргумента E1 (0 ≤ E1 ≤ 0,5h1) характеризует уровень линии тока на выходе из очага деформации, введем безразмерную величину u = 0,5 E1/h (0 ≤ u ≤ 1). Тогда после интегрирования по h находим значение высотного распределения накопленного к концу прокатки значения степени деформации сдвига:

Λ = 2 1+ u2 tg2 αln

h0

срез.

(3.1.36)

 

2 h1

 

 

Для разрывных КВ-полей скоростей часть внутренней энергии рассеивается на поверхностях разрыва вектора скорости. В данном случае это поверхности

E2 = 0 (h = h1) и E2 = – д (h = h0), для которых компоненты вектора скорости

имеют вид (3.1.21), а мощность среза на половине обеих поверхностей h = h1; h = h0 равна τтv0h0u2tg(α/2), где 0 ≤ u ≤ 1. В числителе второго слагаемого форму-

376

3.1. ПРИМЕНЕНИЕ РАЗРЫВНЫХ КВ ПОЛЕЙ СКОРОСТЕЙ

лы (3.1.36) стоит удельная работа Λτт. Исходя из теории размерностей, эквивалентная срезу удельная работа равна мощности среза, деленной на поток V0h0. Тогда эквивалентная скачку вектора скорости степень деформации сдвига в (3.1.36) Λсрез = u2tg(α/2), и окончательно эта формула принимает вид

Λ = 2 1+ u

2

 

2 α

h0

+u

2

 

2 α

 

 

tg

2 ln

 

 

tg

2 .

(3.1.37)

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Подстановкой (3.1.37) в (3.1.35) получим оценку распределения температуры θ, получаемой на выходе из очага деформации вследствие деформационного разогрева металла при плоской листовой прокатке:

 

τ

т

 

 

2

 

2

α

h

 

2

 

2

α

 

 

 

 

 

 

 

0

 

 

 

 

(3.1.38)

θ = θ0 + cρ J

 

1+ u

 

tg

 

2 ln h

+u

 

tg

 

2

.

 

 

 

 

 

 

 

 

 

м

 

 

 

 

1

 

 

 

 

 

 

 

Интегрированием (3.1.37) по u от 0 до 1 получим среднее значение степени деформации сдвига по сечению прокатанной полосы

 

 

 

 

 

 

α 1

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

2

α

2

 

h

1

 

α

 

Λ

ср

=

 

1+ tg

 

+

 

 

ln tg + 1+ tg

 

−1 ln

0

+

 

tg

 

(3.1.39)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

h1

 

 

 

 

 

 

 

 

 

 

2

 

tg

 

2

 

2

 

3

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с соответствующим средним значением деформационного разогрева прокатываемого металла

 

 

 

 

 

θ

ср

= θ

0

+

 

τт

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

ρ Jм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

1+ tg2

α

 

1

 

 

α

 

 

 

+ tg2

 

h

 

1

 

α

 

 

+

 

ln tg

 

+ 1

 

 

−1 ln

0

+

 

tg

. (3.1.40)

 

 

 

 

 

 

 

 

2

 

tg α

 

 

2

 

 

 

 

 

2

 

 

h1

 

3 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для оценки деформационного разогрева (3.1.35) вместо формулы (3.1.36) можно использовать формулу (1.2.165), полученную методом суперпозиции гар-

монических полей скоростей (задача 1.2.7.6) при D = h20 ; d = h21 :

377

3. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ОМД МЕТОДАМИ МСС

 

 

πh

 

πh

 

 

Λпов = 2ln

ctg

1

tg

0

.

(3.1.41)

4 H

 

 

 

4 H

 

 

При этом осевые значения степени деформации сдвига Λос (в первом случае u = 0, во втором d = D = 0) в обоих случаях совпадают и представляются в виде

Λос = 2ln

h0

.

 

 

(3.1.42)

 

 

 

 

 

 

 

 

h1

 

 

 

 

 

В этих случаях из (3.1.35) имеем осевое значение температуры

 

θ = θ

0

+

2τт

 

ln

h0

.

(3.1.43)

 

 

ос

 

cρJм

 

 

h1

 

 

 

 

 

 

В теории прокатки соотношение (3.1.43) получило название формулы С. Финка. Поверхностное значение Λ получается из (3.1.37) при u = 1:

Λпов = 2 1 + tg2 α ln

h0

+ tg

α .

(3.1.44)

 

2 h1

 

2

 

В связи со сложностью формулы (3.1.41) значение Λср из решения, получен$ ного методом суперпозиции гармонических полей скоростей, можно опреде$ лить только численным интегрированием.

Сравнение значений Λпов, рассчитанных по формуле (3.1.44) и по формуле (3.1.41), полученной методом суперпозиции гармонических полей скоростей, а также значений Λос (3.1.42) и Λср (3.1.39) в зависимости от относительного обжа$ тия ε (3.1.28) и параметра q = 0,5 h/R приведено в табл. 12. Для справки: в этой таблице дано значение параметра m (3.1.29), который связан с ε и q формулой

 

2ε

− ε2

 

m =

q

 

.

(3.1.45)

2 − ε

 

Из табл. 12 видно, что обе формулы (3.1.44) и (3.1.42) показывают одинако$ вую качественную тенденцию деформационного разогрева прокатываемого ме$ талла, связанного по формуле (3.1.35) с Λ: поверхностные слои разогреваются больше, чем осевые. Такая тенденция известна на практике, когда отдача тепла в окружающую среду пренебрежимо мала. Правильность количественной оцен$ ки поведения формул (3.1.41) и (3.1.44) необходимо проверять эксперименталь$

378

3.1. ПРИМЕНЕНИЕ РАЗРЫВНЫХ КВ ПОЛЕЙ СКОРОСТЕЙ

Таблица 12. Влияние параметров прокатки на значение Λ

q

ε

m

Значение Λпов по формуле

Λос

Λср

 

 

 

(3.1.41)

(3,1.44)

(3.1.42)

(3.1.39)

 

 

 

 

 

 

 

0,4000

0,1000

0,3684

1,1431

0,3557

0,2107

0,2591

 

0,2000

0,5443

1,4858

0,6596

0,4463

0,5174

 

0,3000

0,6985

1,7875

0,9884

0,7133

0,8051

 

0,4000

0,8478

2,0992

1,3600

1,0217

1,1346

 

0,5000

1,0000

2,4518

1,7946

1,3863

1,5227

0,3000

0,1000

0,4265

1,0247

0,3357

0,2107

0,2524

 

0,2000

0,6318

1,3274

0,6290

0,4463

0,5072

 

0,3000

0,8130

1,5964

0,9470

0,7133

0,7913

 

0,4000

0,9895

1,8787

1,3064

1,0217

1,1166

 

0,5000

1,1706

2,2042

1,7261

1,3863

1,4997

0,2000

0,1000

0,5237

0,8817

0,3123

0,2107

0,2446

 

0,2000

0,7778

1,1427

0,5937

0,4463

0,4954

 

0,3000

1,0035

1,3823

0,9002

0,7133

0,7756

 

0,4000

1,2247

1,6429

1,2468

1,0217

1,0967

 

0,5000

1,4530

1,9532

1,6517

1,3863

1,4748

но. При этом следует отметить два существенных момента: 1) с уменьшением параметра q и увеличением относительного обжатия ε (смещение области параметров толстолистовой прокатки, когда m < 1, в область тонколистовой прокатки, когда m ≥ 1) значения степени деформации сдвига Λ, рассчитанные по формулам (3.1.41) и (3.1.44), сближаются; 2) поверхностные значения Λпов, рассчитанные по формуле (3.1.44), незначительно отличаются от осевых значений Λос, представленных формулой (3.1.42), что иногда в упрощенных расчетах оправдывает оценку деформационного разогрева при прокатке особо тонких изделий по формуле С. Финка (3.1.43), а в более точных расчетах – по формуле (3.1.35)

сиспользованием значения Λср, определяемого по формуле (3.1.39).

3.1.3.Оценка технологических параметров процессов прессования и

волочения круглых прутков

Задача 3.1.3.1. Используя поле скоростей (3.1.14) задачи 3.1.1.2 с изменением текущего радиуса R по линейному закону (3.1.12), определить давление, необходимое для прессования круглого прутка радиуса R1 из заготовки радиуса R0

379

3. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ОМД МЕТОДАМИ МСС

 

 

(рис. 90), полагая, что

 

моделью деформируе-

 

мого металла является

 

идеальная жесткопла-

 

стичная среда.

 

Решение. Используя

 

граничные условия,

 

показанные на рис. 93,

 

составим баланс мощ-

 

ности (1.4.51) внутрен-

Рис. 93. Схема прессования прутка

них и поверхностных

сил без учета массовых

 

 

и инерционных сил.

Сначала рассмотрим мощность внутренних сил Int, которая в данном случае

состоит из двух составляющих: Int = Int1 + Int2. Первая составляющая мощнос-

ти внутренних сил Int1 связана с диссипацией энергии внутри конуса вследствие

отличия от нуля вектора скорости искажения dV в конусе.

Для непрерывных полей скоростей вторая составляющая Int2 равна нулю. Од-

нако в задаче используется разрывное КВ-поле скоростей. Поэтому необходимо

учитывать мощность среза Int2 на поверхностях разрыва вектора скорости.

Для определения мощности внутренних сил Int1 с помощью поля скоростей

(5.19) по формуле Дж. Стокса (2.23) сначала определим компоненты тензора

скоростей деформаций

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

3R2

E

2

 

 

[

rr

[

ΜΜ

0,5[

zz

V

0

tg ; [

rz

[

z r

V

0

r

tg

 

,

(3.1.46)

 

 

 

 

 

 

 

0 R3

 

 

0 2R4

 

 

 

а затем по преобразованной для осесимметричного течения несжимаемых сред формуле (1.2.161)

H 2 3[rr2 [r2z

(3.1.47)

найдем интенсивность сдвиговых скоростей деформаций

 

R2

 

3E2

 

2

 

 

H 2 3V0

0

tg 1

r

tg

 

.

(3.1.48)

R3

4 R2

 

 

 

 

 

 

 

Тогда первая составляющая мощности внутренних сил движения идеальной жесткопластичной среды (Т = Ωт = const) в конусе (рис. 93)

380