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

книги / Остаточные напряжения теория и приложения

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

Рис. 5.19. Заннснмость коэффициента5 тепло­ отдачи ас от нременп

ыие оценок уровня остаточных напряжения по их текущим значе­ ниям (кривые 2 н 4) и по функционалу (3.42) (кривые 1 и 3) при естественном (кривые 1 и 2) и принудительном (кривые 3 и 4) охлаж­ дении. Видно, что оценка (3.42) выполняется с большой точно­ стью. Пики кривых в начальный период охлаждения обусловлены структурными превращениями исследуемой области. Эпюры ос­ таточных напряжений при естественном и оптимальном принуди­ тельном охлаждении содержатся на рис. 5.21. Сравнеппе обнару­ живает уменьшение максимальных сжимающих напряжений в центре стенки с —85 до —25 МПа и максимальных растягиваю­ щих напряжений в переходной от стенки к полкам области от +160 до 80 МПа. Заметим, что существовавшее опасепие образо­ вания при ускоренном охлаждении закалочных структур пе под­ твердилось. Из приведенной на рис. 5.22 термокинетической диа­ граммы распада переохлажденного аустепита [90] и кривых охлаждения для нескольких характерных точек исследуемой облас­ ти следует, что как при естественном, так и при принудительном охлаждении происходит перлитное превращение.

Исследование микрошлифов (при увеличении в 200 раз) так­ же показало, что при принудительном охлаждении по сравнению с естественным не изменяется тип структуры и лишь незначитель­ но изменяется размер зерна (структура делается несколько более мелкозернистой), что не ведет к сколько-нибудь существенному изменению физических или механических свойств.

Ввиду сложности технической реализации предлагаемого ре­ жима охлаждения была решена также несколько упрощенная за­ дача. Предположим, что тип системы охлаждения, расход и дав­ ление охлаждающей среды известны. Тем самым определен и ко­ эффициент’теплоотдачи как функция температуры охлаждаемой поверхности. Тогда определению подлежат только ширина в при­ нудительно охлаждаемой зоне и время охлаждения хг. Принуди­ тельное охлаждение осуществляется непрерывно в течение одного интервала времени до момента времени, с которого начинается увеличение значения функционала (3.42) [77].

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

81

мых двутавров балок примерно одинакова и составляет Ь =

=(0,5—0,6)2?, где В — ширина полок.

На рис. 5.23 представлены кривые изменения функционала

Ф (t) (3.42) при охлаждении широкополочной двутавровой балки 70HI1 от начальной температуры (температура конца прокатки) до температуры окружающей среды Г » = 20° С. Кривая 1 соот­ ветствует охлаждению профиля на воздухе. Продольные оста-

Рпс. 5.20. Уровень напряжений в исследуемой области, определяемый по напряжениям и из оценки (3.42)

J, 2 — при естественном охлаждении; 3, 4 —принудительном

Рис. 5.21. Эпюры продольных ос­ таточных напряжений

1 — естественное охлаждение; 2 — при­ нудительное

82

точные напряжения в характер­ ных точках поперечного сече­ ния (рис. 5.22) в этом случае составляют ах = — 170, ст2 = = - f 200, а 3 = — 160 МПа. Кри­ вые 2—4 соответствуют прину­ дительному охлаждению водой с удельным расходом q = 0,008 м3/м2 в течение 9,6, 18,5, 30 с соответственно.

Соответствующие кривой 2 (время принудительного охлаж­

дения

t j — 9,5

с) величины

остаточных

напряжений

рав­

ны:

ст* =

— 117,

а, =

145,

о3 = —169 МПа. Римскими циф­ рами на рисунке обозначены мо­ менты окончания принудитель­ ного охлаждения для соответст­ вующих кривых.

Кривая 3 описывает измене­

ние |Р |(т) для

оптимального

(с точки

зрения

наименьшего

уровня

остаточных напряже­

ний) времени охлаждения т* =

= 18,5 с. Остаточные напряже­ ния в этом случае составляют 21, 152, — 43 МПа соответст­ венно в точках 1—3.

Наконец, кривая 4 соответст­ вует времени принудительного охлаждения тх = 30 с, превы­ шающему оптимальное. Для этого режима после достиже­ ния некоторого минимального значения |Р |(т) вновь начи­ нает возрастать и, вообще гово­ ря, может превысить величину |Р |при естественном охлаж­ дении па воздухе. Напряжения, соответствующие кривой 4, рав­

ны:

= 160, о2 = 176, а3 =

= - 6

МПа.

На

рис. 5.24 приведены эпю­

ры остаточных напряжений для естественного и оптимального принудительного охлаждепия балки 70Ш1. Для других рас­ сматриваемых профилей получе-

Рнс. 5.22. Зависимость температуры от времени длясвободного (/) и принудительного {II) охлаждения (номера кривых 1,2, 3 соответствуют точкам 1, 2, 3, изображенным на рисунке внизу)

Рис. 5.23. Изменение величины функционала |Р |[ер (т)] в про­ цессе свободного (1) и принудитель­

ного {2—4) охлаждения для балки

70Ш1

Т а б л и ц а

5.4. Оптимальное время принудительного охлаждения

 

и остаточные напряжения и двутавровых профилях

 

 

 

Удельный

Оптимальное

Остаточные напряжения сх в точке, МПа

Тип балки

время прину­

расход воды а.

дительного

 

 

 

 

М*/(М*'С).

охлаждения

1

2

3

 

 

ti, с

30Б1

0,008

9,5

-1S0

183

-130

 

-21

153

-75

 

0,013

6,5

19

137

-73

ЗОШ

0,024

3,5

39

184

—80

_

_

-147

187

—153

 

0,008

9,5

29

130

—24

 

0,013

6,5

23

149

-19

60Б1

0,024

3,5

21

151

-11

_

_

-180

214

-139

 

0,008

15,5

-18

157

-31

 

0,013

12,5

-4

152

-34

70Б1

0,024

9,5

34

123

-23

_

_

-189

198

-140

 

0,008

18,5

40

.157

-45

 

0,013

12,5

3

156

-42

 

0,024

12,5

46

•150

-50

7 o m i

_

_

-168

200

-160

 

0,008

18,5

21

152

-43

 

0,013

12,5

-9

150

-74

№60

0,024

9,5

33

148

-98

_

-256

185

30

 

0,008

26,0

65

176

-42

 

0,013

20,0

81

168

-50

 

0,024

12,5

23

154

-20

ны аналогичные результаты. Значения остаточных напряжений в характерных точках поперечного сечения (см. рис. 5.4) для неко­ торых двутавровых балок сведены в табл. 5.4.

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

Полное решение задачи оптимизации, требующее решения тем­ пературной задачи при прокатке, задачи термоупругопластичности для исследования охлаждения профиля на холодильнике и задачи поиска управляющей функции ас (т), занимает 2,5 ч на БЭСМ-6. Результаты расчета показывают существенное снижение собст­ венных напряжений при использовании принудительного охлаж­ дения по сравнению с естественным. Однако эшоры напряжений в этом случае имеют более резкий характер, что связано с «жест­ костью» режима охлаждения.

Результаты использования данного упрощенного подхода для балки № 60 (материал — сталь 09Г2Д) были проверены экспери-

Рис. 5.24. Остаточные напряжения D широкополочпой балке 70Ш1 при ес­ тественном (7) и принудительном (2) охлаждении

ментально в промышленных условиях. Применялась охлаждающая установка с удельным расходом воды 8.3 -КГ3 м3/(м2-с) и щелевы­ ми форсунками. В условиях эксперимента испольаовались рас­ считанные параметры охлаждения: ширина охлаждаемой зоны b = 0,5 В, время охлаждения хг = 20 с. Полученное соответствие экспериментальных и теоретических данных можно признать удов­ летворительным: расчетные и экспериментальные значения на­ пряжений составляютдля центра стенки —114 и —100 МПа, для переходной области +81 и + 60 МПа, для кромки полки +20 и + 40 МПа соответственно. Напомним, что такой же уровень на­ пряжений достигается при отжиге в течение 1800 с при температу­ ре 450° С.

Предложенный алгоритм был испольэовап также для определе­ ния оптимальных режимов охлаждения несимметричных профилей для уменьшения остаточных напряжений и искривленности. В этом случае целевая функция имеет вид (3.20), а базис (3.31) ока­ залось достаточным ограничить членами, содержащими коорди­ наты в степени не выше первой (рис. 5.25). Результаты расчета подтверждают возможности рассмотренного в настоящей работе

Рис. 5.25.

Уровень на­

пряжений

в рессоре

45

X б,

определяемый

по

напряжениям (2), по

оценке с базисными функ­

циями

в виде

констант

(2), по

оценке

с

базис­

ными

 

функциями в

виде

линейных

функ­

ций

координат

 

(3)

85

подхода для разработки режимов охлаждения различных горя­ чекатаных профилей. Так, при использовании рассчитанных оптимальных режимов охлаждения рессорной полосы 45 х 6 мм (охлаждению подвергается пижняя кромка полосы, удельный рас­ ход воды 8,3-Ю "3 м3/(м2*с), время охлаждения = 4,12 с) оста­ точные напряжения снижаются в 2—3 раза, прогиб на единицу длины готового профиля уменьшается до 2,3 мм/м (при естествен­ ном охлаждении величина прогиба составляет 4,5 мм/м).

5.2. Решение задач расчета и снижения уровня остаточных напряжений при волочении проволоки

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

Рис. 5.26. Схеме четырех­ кратного волочильного стана

(а)

п одной из

се

клетей (б)

1

разматывающее

устройство;

2— 5 — волочильные

клети,

в — сматывающее

устройство;

7 — монолитная

волока;

8 — проволока,

 

 

9

— барабан

 

 

тая проволока находит применение в канатах, пружинах, в ка­ честве арматуры и т. д. Комплекс свойств, определяющих качество готовой проволоки, в частности уровень остаточных напряже­ ний, в основном формируется на последних этапах ее производ­ ства. Особенности деформации проволоки в монолитных воло­ ках приводят к появлению остаточных напряжений, достигающих 70—80% предела упругости металла [102].

Для исследования процессов волочения проволоки в основном применялись экспериментальные или аналитические методы [12, 13, 19, 37, 45, 46, 49, 85,118,128], причем последние с использова­ нием существенных допущений. Экспериментальные методы отли­ чаются большой трудоемкостью, привлечение же значительного числа допущений для аналитического решения дает возможность получить лишь интегральные характеристики, такие, как средняя температура, усилие волочения и т. д. Методика, излагаемая в дан­ ной работе, дает возможность рассчитать температурное поле и напряжения в течение всего времени пребывания проволоки на стане, в том числе и остаточные напряжения. Затем с помощью результатов, описанных в гл. III, проводится решение задачи снижения остаточных напряжений. Конкретные результаты ана­ лизируются для волочильного стана Белорецкого металлурги­ ческого комбината им. М. И. Калинина.

16

Кратко опишем технологию волочепия на этом стане. Стальная щоволока перед чистовым волочением находится на катушке. С разматывающего устройства 1 (рис. 5.26) проволока поступает на волочильный стаи, состоящий из 4 монолитных волок из победита в стальной обойме, имеющей внутреннее охлаждение водой. Температура волоки поддерживается в пределах 250— 300° С. Применяемая смазка — мыльный порошок. Через волоку проволока протягивается с помощью барабана, охлаждаемого из­ нутри водовоздушной смесыо, снаружи — естественным образом па воздухе. После четвертого барабана проволока подается на сматывающее устройство 6. Расстояние между волокой и осью барабана каждой клети равно 0,3 м; диаметр барабанов — 0,75 м;

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

Напряженно-деформированиое состояние и температурное по­ ле при волочении можно найти из решения связанной задачи термоупругопластичности (2.46) — (2.51). Отметим, что эффект свя­ занности здесь играет важную роль. При холодном волочении про­ волоки ив высокоуглеродистых сталей Ст. 70 и 80 со скоростью выхода из стана более 4 м/с температура поверхности из-за трения и деформационного разогрева поднимается уже в первой волоке до 200—270° С в зависимости от скорости волочения [62], а в послед­ них волоках непрерывных станов холодного волочения достигает 700° С [49], хотя средняя температура проволоки поднимается при этом лишь до 300—350° С [97]. Выравнивание температурного поля в поперечном сечении проволоки после волоки проходит за время порядка 1 мкс [97]. Существенная неоднородность температурного поля в волоке из-за неоднородности деформаций и влияния кон­ тактного трения, в свою очередь, влияет па механические харак­ теристики, прежде всего па предел упругости.

Задача исследования папряжешю-деформированного состоя­ ния рассматривается как стационарная для всего волочильного стана в целом; используется эйлерова система координат. Для зоны активного нагружения применяется приведенный в § 4.2

87

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

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

Конкретизируем далее граничные условия. При входе проволо­

ки в первую волоку действует противонатяжение Q(здесь и далее символ «Л» относится к заданным величинам) смотки проволоки с катушки (поверхность на рис. 5:27), здесь же задана темпе­

ратура входа на стан f . В промежутках между барабанами и во­ локами происходит естественное охлаждение проволоки на воз­ духе при отсутствии поверхностной пагрузкп (поверхность S2). На контактной поверхности в волоке (поверхность £ за) нормаль­ ная компонента скорости металла равна нулю из-за непроницае­ мости в волоку. В касательном направлении действуют силы тре­ ния. Анализ законов трения показал, что здесь паиболее целесо­ образно применить закон Амонтона — Кулона [85, 86, 1061, который обычно применяют при анализе развитых пластических течений, пока контактные касательные напряжения, определенные по нему, не превысят предела упругости на сдвиг т3, что соответ­ ствует нормальному давлению порядка ст„ [85]. Условие пластич­ ности для процесса волочения можно написать в виде ог — <хг = = <тв, где <зг — продольное, аг — радиальное напряжение. Так как кроме небольшого возможного участка трехосного сжатия на входе у поверхности волоки az > 0, оу < 0, следует, что оу < as. Ввиду малости угла наклона конусной части волоки <7Г ^ ап (оп — нормальное давлепие), и отсюда следует, что ап не достигает ве­ личины оу и оправдывает применение закона Амонтона — Куло­ на. При записи температурных граничных условий учтено, что теплообмен между волокой и проволокой ограничен [85], поэтому в расчетах им обычно пренебрегают. Использовано также ранее упомянутое утверждение о том, что в металл идет половина тепла от трения [71J. После выхода из волоки (поверхность S3Q) заданы скорость волочения в £-ой волоке vit равенство нулю продольного теплового потока и касательных напряжений в поперечном сече­ нии. Из решения определяется усилие волочения P t после i-й волоки.

Барабаны волочильного стана передают усилие волочепия на проволоку. Для исключения полного проскальзывания проволо­ ки в окружном направлении на поверхности барабана и обеспече­ ния тем самым устойчивой скорости волочения на барабанах воло­ чильных машин находится определенное число витков (около 10). Каждый новый виток при поступлении па барабан заставляет скользить остальные в продольном направлении барабана. После выхода из волоки при прохождении через барабан происходит

•в

разгрузка элементов проволоки от усилия Р $до величины проти-

вонатяжения QHl в следующей волоке. Охлаждение па барабане учитывается с помощью усредненного коэффициента теплоотдачи а (поверхность S4). Ha участке охлаждения перед смоткой (поверх­ ность iS5) отсутствуют поверхностные силы, задал теплообмен по закону Ньютона. Наконец, на торцевой поверхности перед смот­

кой (поверхность SQ) задана скорость волочения v, а тепловой по­ ток равен нулю.

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

вид (для решения

используются

цилиндрические

координаты)

02 =

О.

rrz — Hfc — Oi

T T,

x6E*?i;

 

<т =

0,

X (dT/dn) =

— a(T — T„),

x e

Si, St, Ss-

vn =

0,

Ti< = /o n,

X(dTldn) =

^ -W i,

x E ISM!

vz=

vit

Trz =

T0z= O ,

дГ/<?2 = 0,

x €E £ 36;

 

<*z =

Qi>

Trz =

T02 =

O,

T =

Ti,

хЕЕдУ3в;

 

Vz =

V,

T „ =

To2 =

0,

dT/dz= 0,

x e ^ c .

 

Здесь Wx — источник от

трения;

Т\ — начальная

температура

перед входом в i-ю волоку; / — коэффициент трения

(в расчетах

был принят для данных условий смазки равным 0,07 согласно [86]). Решение поставленной задачи термоупругопластичности про­ водилось методом конечных элементов (МКЭ) согласно алгорит­ мам, описанным в гл. IV. При волочении используется сталь Ст. 80. При анализе данного процесса можно пренебречь явлением ползучести, так как средняя температура проволоки не превышает к концу процесса 300—350° С, а поверхностная температура лишь на доли секунды в волоке поднимается до 600—700° С. Для опре­ деленного сечения проволоки длительность пребывания на воло­ чильном стане не превышает нескольких минут. По данным работ [87, 93 и др.], явления релаксации напряжений для сталей при температуре 500—700° С проявляются при выдержке в течение не­

скольких часов.

Скалярные свойства материала экспериментально изучены в отраслевой лаборатории прокатного производства Магнитогор­ ского горно-металлургического института. Зависимость сопротив­

ления

пластической

деформации

аа

от степени деформации

е и

температуры

Т

имеет

вид

о3 =

9,8 (80—0,057 Т +

-Ь 10,3 (100 е*)0»48 МПа.

Для выбора теории

пластичности, пра­

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

Сначала приведем траекторию для пятимерного пространства напряжений А. А. Ильюшина Sh [28]. В нашем случае вследствие осесимметричности задачи траектория лежит в трехмерном под-

ИПа ■fff~2

Рис. 5.28. Зависимость компонент вектора напряжений — 53 от продолжи­ тельности пластической деформации t при волочепии (1 — па поверхности, 2 — в центре)

Ряс. 5.29. Зависимость компонент вектора деформации 9X— Э3 от продолжи­ тельности пластической деформации t при волочении (1 — па поверхности,

.2 — 0 центре)

Рис. 5.30. Траектория деформирования для элемента в центре проволоки при

Рис. 5.31. Изменение радиуса кривизны траектории деформирования от вре­ мени при волочении (2 — в центре, 2 — па поверхности)

пространстве. Вводим трехмерный вектор напряжений S = £ iC i+

+ $гег +

*^зе 31 ГД0 е £(£ = 1, 2, 3) — орты. Здесь компоненты век­

тора $i определяются выражениями

St =

v \ [s r cos (Р - f

я/6) — so sin |3],

St =

/ 2 [sr sin (p +

я/6) - f s0 cos pi, S3 = j/J xrz,

где sr, so, Trz — компоненты девиатора напряжений. При анализе

принято р = 0.

ЯР