Анализ термоупругих напряжений в поковках с учетом эффекта релаксации Текст научной статьи по специальности «Физика»

УДК 669.14 : 621.78.082 : 51-74

АНАЛИЗ ТЕРМОУПРУГИХ НАПРЯЖЕНИЙ В ПОКОВКАХ С УЧЕТОМ ЭФФЕКТА РЕЛАКСАЦИИ

Д.А. Мирзаев, А.Д. Шабуров, А.О. Чернявский

Поскольку напряжения являются одним из основных факторов образования флокенов, важно понимать характер напряженного состояния непрерывно охлаждающейся цилиндрической поковки. Для этого был выполнен компьютерный расчет распределения внутренних напряжений, а также учтена возможность релаксации упругой деформации при охлаждении поковки. Расчет был произведен по нескольким вариантам при различных скоростях охлаждения для трех компонент напряжений - осевой, радиальной и тангенциальной. Расчет эффекта релаксации произведен на основании обработки и аппроксимации экспериментальных данных по ползучести стали марки 40ХГМ в температурном интервале охлаждения после противофлокенной термической обработки. Результаты расчета показали, что если при аналитическом расчете без учета релаксации абсолютная величина напряжений уменьшается и стремится к нулю при подходе к комнатной температуре, то расчет с учетом релаксации внутренних напряжений показывает, что после окончания охлаждения в поковке сохраняются остаточные напряжения, причем на поверхности они сжимающие, а в центре поковки растягивающие. Также было показано влияние данных напряжений на раскрытие зародыша флокена.

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

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

Ключевые слова: охлаждение на воздухе, внутренние напряжения, релаксация напряжений, ползучесть, расчет напряжений, моделирование охлаждения, поковки, флокенообразование, флокены, сталь 40ХГМ.

Введение

Появление флокенов, т. е. заполненных водородом трещин, при охлаждении поковок обусловлено выделением водорода в микротрещины вследствие уменьшения растворимости водорода в а-железе и их роста из-за непрерывного возрастания давления газа молизующегося водорода Н2 [1]. Заметное влияние на рост флокенов оказывают термические напряжения [2]. Температурный интервал образования флокенов составляет 200...20 °С. Поэтому анализ напряжений необходим, прежде всего, для низкотемпературной области на заключительной стадии охлаждения поковок.

В работе [3] был выполнен аналитический расчет напряжений в охлаждаемых поковках, основанный на известных уравнениях термоупругости [4]. Прежде всего, он подтвердил, что «амплитудные» величины напряжений сг, сг и се зависят от разности температур между поверхностью и осью цилиндра.

Максимальные значения растягивающих напряжений сг, и се наблюдаются у поверхности, тогда как сжимающих - в центре. Таким образом, при движении вдоль радиуса они изменяют свой знак. Напряжение сг только сжимающее. У поверхности оно равно нулю, а максимум модуля наблюдается у оси.

По мере возрастания длительности охлаждения разность температур между центром и поверхностью цилиндра уменьшается, соответст-

венно, снижаются величины всех напряжений, а при достижении осевой зоной комнатной температуры упругие напряжения должны исчезнуть. Именно в этом заключается основной недостаток предыдущих расчетов, поскольку имеются данные [5] о сохранении в поковках остаточных напряжений. Этот эффект связан с явлением релаксации на различном расстоянии от оси. В данной статье обсуждаются результаты расчёта термических напряжений с учетом их релаксации в охлаждаемых поковках их стали 40ХГМ, поставляемой ОАО «Уральская кузница» на экспорт, а потому технологические аспекты её тщательно изучаются.

Методика расчета и обсуждение результатов

В основе использованного метода анализа лежат следующие соображения [6, 7]. Допустим, что образец с исходной длиной 10 был упруго растянут так, что его удлинение составило Д10, а затем зане-волен. Упругие напряжения в таком образце будут

равны ст0 = Е —- = Ее0 , однако в ходе выдержки 1о

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

ст 1 дст(г, Т) так как теперь е0 = —ь еп или----- = -еп.

0 Е Е дг п

Рис. 1. Кривые ползучести стали марки 40ХГМ при различных параметрах испытания (температура-нагрузка): 1 - 560 °С, 150 МПа; 2 - 560 °С, 220 МПа; 3 - 600 °С, 150 МПа; 4 - 600 °С, 190 МПа; 5 - 600 °С, 230 МПа; 6 - 650 °С, 90 МПа;

7 - 650 °С, 150 МПа

0) х я

<D К SS

ч п

о я

Л

ч

(и н К о о X н

о

3)

10

Время,1

15

20

Рис. 2. Сравнение значений аппроксимирующей функции с экспериментальными данными для ползучести стали марки 40ХГМ: 1 - 560 °С, 150 МПа; 2 - 560 °С, 220 МПа; 3 - 600 °С, 150 МПа

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

1 dg(T, T) = -dsn, (1)

E

то есть развитие ползучести приводит к уменьшению (релаксации) напряжения. Если для данной стали известна деформация ползучести 8п как функция напряжения и температуры, то уравнение (1) позволяет установить связь дифференциалов напряжения и деформации ползучести, необходимую для численных расчетов. В справочных данных пакета конечно-элементного моделирования ANSYS для деформации ползучести приведены несколько аппроксимирующих функций еп(х, T), на основе которых построены программы численного интегрирования уравнения (1) для всех компонент тензора напряжений совместно с уравнениями теплопроводности и термоупругости.

С целью подбора функции еп(х, T) было про-

ведено экспериментальное исследование на ползучесть образцов стали 40ХГМ. От промышленной поковки стали 40ХГМ в кузнечно-прессовом цехе ОАО «Уральская кузница» были вырезаны заготовки 30x30x200 мм из середины радиуса поковки. Затем заготовки были подвергнуты отжигу по режиму: нагрев до 850 °С, выдержка 3 ч, охлаждение 2 ч с печью, далее на воздухе. Данный режим позволил получить в заготовках равновесную феррито-перлитную структуру, такую же, как после проти-вофлокенного отжига. Далее из заготовок были выточены стандартные образцы на ползучесть с диаметром рабочей части 7 мм и произведены испытания по ГОСТ 3248-811. Испытания на ползу-

1 Испытания были проведены на ОАО «Уральская кузница» (г. Чебаркуль) под руководством начальника лаборатории физико-механических испытаний ЦЗЛ Т.В. Ручьевой._

честь производились по нескольким вариантам при температурах 560-650 °С, соответствующих температурному интервалу охлаждения поковок после отжига, т. е. получены кривые ползучести (рис. 1, 2).

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

еп =— + схст 2т Е

С2тСз+1е-С„/Т .

1

Сз +1

+С стС е

С е-с7/ Т т

(2)

для которой неизвестные коэффициенты Сь С2, ..., С7 были определены путем последовательного сопоставления с экспериментальными зависимостями

(8п )т=ссп81 = /1 (ст) ; (8п )а=ссп81 = /2 (Т) и т. д. Для областей нестационарной и стационарной ползучести: С = 78; С2 = 2,17; С3 = -0,55; С4 = 14 000; С5 = 750,4; С6 = 4,55; С7 = 28 240. На рис. 2 (1-3) сравниваются экспериментальные кривые зависимости и их аналитическая аппроксимация; степень их согласия можно признать достаточно хорошей.

Расчет напряжений сг, сг, с0 ист при охлаждении цилиндрических поковок с учетом релаксации был проведен в пакете ANSYS путем совместного численного решения методом конечных элементов уравнения теплопроводности, возникновения упругих напряжений и уравнения вида (1) для релаксации.

Параметры расчета были приняты такими же, как в работе [3]. Рассматривается цилиндрическая поковка из стали 40ХГМ диаметром 450 мм.

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

1) на воздухе (рис. 3 и 4), при этом коэффициент теплоотдачи а, как и в работе [3], был принят равным 34 Вт/м2К ;

С учётом релаксации

Без учета релаксации

I

К

о

-од

---ОТ /

---ОТ ' / Дя

/

/У / ✓

- ----

—ое

--— ---<зт

---от

б)

г)

20 ю Н о

-10 -

-об

---02

---ОТ У/' я

1

— -

- - - '

д)

---Я

1 гц,ч. ч 1

-аб V

---стг

---от

-2 ^

-4

20 10 0 -10

-20 -30

ж)

Рис. 3. Результаты расчета внутренних напряжений при охлаждении на воздухе поковки диаметром 450 мм с учетом (а, в, д, ж) и без учета релаксации (б, г, е, з): а, б - через 0,5 ч; в, г - через 5 ч; д, е - через 7 ч; ж, з - через 15 ч после начала охлаждения

20 -

О -20 -40 -60 -80

I

К

С учётом релаксации

- -*■ ~ -----

/ - -ОН

! > -ах

1; -- -от

а)

г — г- ' ' Г"

X Л*

Л У/ ~ ! -об - (32

-от

Без учета релаксации

20 0

-120

1 / -| / [ / -ОН

/ / ---ах ---ог

б)

/ Л г ■ / /' . / ч -об

/ ---ах

---от

25 Время, ч

Рис. 4. Результаты расчета внутренних напряжений при охлаждении на воздухе поковки диаметром 450 мм в зависимости от длительности охлаждения с учетом (а, в, д) и без учета релаксации (б, г, е): а, б - на оси; в, г - в середине радиуса; д, е - на поверхности

поковки

2) при замедленном охлаждении (рис. 5) со скоростью ~13 °С/ч (а = 5 Вт/м2К ).

В случае малой длительности охлаждения напряжения с0 и сг резко уменьшаются по мере удаления от поверхности и меняют свой знак, оставаясь затем приблизительно одинаковыми, но малыми по величине. Различия в величине напряжений ещё не велико, но уже после 30-минутного охлаждения (рис. 3, а, б) с0 и сг у поверхности отличаются почти в 1,5 раза. После 30 мин охлаждения существенно возрастают и сжимающие напряжения с0 и с2 вблизи оси поковки. А это означает, что теперь релаксация напряжений будет заметно проходить не только у поверхности, но и в осевой зоне. Отметим также, что увеличение выдержки приводит к уменьшению г0, при котором с0 обращается в 0. После 30 мин охлаждения температура в середине поковки составляет 633 °С, а у поверхности - 573 °С. Фактически не охваченным релаксацией остается участок между осевой и приповерхностной зоной, в котором напряжения с0 и сг оказываются небольшими. Релаксация напряжений у поверхности (где напряжения растягиваю-

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

Расчет для поковки, охлаждающейся в течение 5 ч (рис. 3, в, г), показывает, что график радиального изменения напряжений без учета релаксации принципиально сохраняет свой вид, но величины напряжений у оси и на поверхности уменьшаются в силу снижения разности температур между поверхностью и центром (Тп = 193 °С, Тц = 213 °С). В этом случае эффект релаксации приводит к уменьшению с0 и с2 и в осевой, и в приповерхностной зонах. Вблизи г = R немного меняется характер изменения напряжений. Вместо восходящего участка появляется практически горизонтальный участок (рис. 3, в).

Для 7-часового охлаждения (рис. 3, д) характерно отчетливое формирование трех зон вдоль радиуса поковки. Дальнейшее уменьшение разности температур поверхности и центра (Тп = 129 °С, Тц = 142 °С) снижает термические напряжения.

Без учета релаксации

а)

&

ев

К

-4 -

-6 -

_- -- ~ " . — - —

—V 1 1 1

С /,'

//

■ / -ае

1 ---ох

---от

15 10

100 Время, ч

Рис. 5. Результаты расчета внутренних напряжений п диаметром 450 мм в зависимости от длительности учета релаксации (б, г, е): а, б - на оси; в, г - в сере

поковки

А поскольку поверхность образца из-за релаксации была увеличена, то теперь охлаждающаяся серединная зона начинает создавать в поверхностных слоях сжимающие напряжения ае. Кривые напряжений, при расчете которых релаксация не учитывалась (рис. 3, е), по-прежнему отражают растягивающие напряжения у поверхности и сжимающие напряжения на оси. В данном случае видно, что ае и аг определяются разностью температур поверхности и центра.

После 15-часового охлаждения отчетливо наблюдается (рис. 3, ж) возникновение около оси поковок растягивающих напряжений ае, аг и аг, а у поверхности - сжимающих напряжений ае и аг, величина которых достигает 20 МПа. Расчет в варианте отсутствия релаксаций (рис. 3, з) дает ~ 1,5 МПа на поверхности и ~ 2 МПа у оси.

Более подробно временной характер изменения напряжений ае, аг и аг на оси, при г = 1/2 R и у поверхности представлены на рис. 4 с учетом релаксации и без неё. Общий характер изменения напряжений оказывается следующим: в отсутствие релаксации напряжения у поверхности растяги-

вающие, а напряжения при г = 1/2 R и на оси -сжимающие. Причем в первых двух случаях (рис. 4, б, г) их графики проходят через минимум и затем начинают возрастать по модулю асимптотически, приближаясь к 0 при продолжении охлаждения. Аналогично и на поверхности, пройдя через максимум, ае и аг уменьшаются до 0 при дальнейшем охлаждении. При учете релаксации (рис. 4, а, в, д) наблюдается, во-первых, уменьшение уровня всех напряжений, может быть, за исключением очень малых времен охлаждения, причем при времени охлаждения более 7 ч графики напряжений пересекают 0 и, по существу, напряжения превращаются в остаточные напряжения растяжения около оси и сжатия у поверхности. Аналогично и для г = 12 R . При этом сжимающие напряжения будут препятствовать образованию флокенов, в которых трещина перпендикулярна действующему напряжению. Во-вторых, следует учитывать, что сжимающие напряжения вдоль е и г будут препятствовать образованию флокенов, плоскость которых перпендикулярна этим напряжениям, поэтому поверхность возникающего флокена

должна быть ориентирована перпендикулярно к радиусу, что в действительности довольно часто наблюдается в производстве и отмечается в публикациях [8].

Рассмотренные выше графики внутренних напряжений относились к охлаждению поковки диаметром 450 мм на воздухе. Теперь рассмотрим случай, когда поковки охлаждаются замедленно со скоростью ~ 13 °С/ч, что примерно соответствует скорости охлаждения поковки в термосе [9]. Зависимость внутренних напряжений от времени охлаждения в двух вариантах (с учетом и без учета релаксации напряжений) представлены на рис. 5. Характер изменения напряжений оказывается таким, как и в решении, рассмотренном в случае охлаждения на воздухе (см. рис. 4). В отсутствие релаксации сохраняется общий характер изменения напряжений по радиусу и по времени, но все напряжения оказываются существенно меньшими (рис. 5, б, г, е). В случае учета релаксации характер изменения напряжений также не меняется, но величина остаточных напряжений после охлаждения снижается приблизительно в 3 раза (рис. 5, а, в, д).

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

(Т -4- (Т _ 4- (Т

-. (3)

Dh = D0eRTe

3

Поэтому в выражении для коэффициента диффузии водорода появляется дополнительный сомножитель

-б ст¥м_

кт (4)

то есть энергия активации заменяется на Q ^ Q -ст¥М, где ¥М - парциальный молярный

объем растворенного водорода, который оценивается в 2 см3/моль. Если максимальное значение ст при охлаждении на воздухе у поверхности оцени-

ч -20-20 + 0 „ ^ вается (см. рис. 3, ж) ---= -13,3 МПа,

то дополнительный множитель, отражающий эффект уменьшения диффузии, будет равен

СТУМ -13,3106 210-6

е кт = е 8,314-з00 = 0,99, то есть влияние напряжений на диффузию водорода пренебрежимо мало.

Однако на самой поверхности водород, по-видимому, удаляется ещё в процессе термической обработки слитков, нагрева под ковку, самой ковки и охлаждения после неё, так как при ковке про-

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

Выводы

1. При охлаждении цилиндрических поковок в однофазной а-области возникают осевые (стг), радиальные (стг), и тангенциальные (ст0) напряжения. Если при аналитическом расчете без учета релаксации абсолютная величина напряжений уменьшается и стремится к нулю при подходе к комнатной температуре, то расчет с учетом релаксации внутренних напряжений показывает, что после окончания охлаждения в поковке сохраняются остаточные напряжения, причем на поверхности они сжимающие, а в центре поковки растягивающие.

2. Наличие на поверхности сжимающих напряжений препятствует раскрытию флокенов.

3. Показано, что влияние внутренних напряжений на диффузию водорода пренебрежимо мало.

Работа выполнена в рамках государственного задания Минобрнауки РФ и частично поддержана грантом РФФИ 14-03-00618.

Литература

1. Дубовой, В.Я. Флокены в сталях: моногр. / В.Я. Дубовой. -М.: Металлургиздат, 1950. - 332 с.

2. Склюев, П.В. Зависимость остаточных напряжений от скорости охлаждения после отпуска /П.В. Склюев, Б.Д. Петров // Проблемы металловедения и термообработки. - 1960. - Вып. 2. -М.; Свердловск: Машгиз. - С. 136-142.

3. Мирзаев, Д.А. Внутренние напряжения при охлаждении крупных поковок / Д.А. Мирзаев,

A.Д. Шабуров // Вестник ЮУрГУ. Сер. «Металлургия». - 2012. - Вып. 18, № 15 (274). - С. 112-115.

4. Тайц, Н.Ю. Технология нагрева стали / Н.Ю. Тайц. - М.: Металлургиздат, 1962. - 567 с.

5. Склюев, П.В. Термическая обработка крупных поковок / П.В. Склюев. - М.: Машиностроение, 1976. - 48 с.

6. Работнов, Ю.Н. Сопротивление материалов / Ю.Н. Работнов. - М.: Физматгиз, 1962. - 456 с.

7. Борздыка, A.M. Релаксация напряжений в металлах и сплавах / А.М. Борздыка, Л.Б. Гецов. -М.: Металлургия, 1978. - 256 с.

8. Штремель, М.А. Кинетика раскрытия внутренней зернограничной трещины водородом / М.А. Штремель, А.А. Князев // Физика металлов и металловедение. - 1986. - Т. 62, № 4. - С. 645-651.

9. Шабуров, А.Д. Теплообмен при охлаждении поковок в термосе / А.Д. Шабуров, Д.А. Мирзаев,

B.А. Смолко // Вестник МГТУ им. Г.И. Носова. -2011. - № 2. - С. 65-67.

Мирзаев Джалал Аминулович, д-р физ.-мат. наук, профессор кафедры физического металловедения и физики твёрдого тела, Южно-Уральский государственный университет (г. Челябинск); mirzayev@physmet. susu.ac.ru.

Шабуров Андрей Дмитриевич, и. о. начальника кузнечно-прессовой лаборатории, Челябинский филиал ОАО «Уральская кузница» (г. Челябинск); [email protected].

Чернявский Александр Олегович, д-р техн. наук, профессор, заведующий кафедрой прикладной механики, динамики и прочности машин, Южно-Уральский государственный университет (г. Челябинск); cher@sopro. susu. ac.ru.

Поступила в редакцию 5 сентября 2014 г.

Bulletin of the South Ural State University

Series "Metallurgy" _2014, vol. 14, no. 3, pp. 40-47

ANALYSIS OF THERMOELASTIC STRESSES IN FORGED PRODUCTS ON ACCOUNT OF THE RELAXATION EFFECT

D.A. Mirzaev, South Ural State University, Chelyabinsk, Russian Federation, mirzayev@physmet. susu. ac. ru.

A.D. Shaburov, Chelyabinsk Branch of the Urals Stampings Plant, Chelyabinsk, Russian Federation, [email protected].

A.O. Chernyavsky, South Ural State University, Chelyabinsk, Russian Federation, cher@sopro. susu. ac. ru.

Since stress is one of main factors in flake formation, it is important to understand the nature of the stress state of continuously cooling cylindrical forgings. For this purpose computer calculations of internal stress distribution were made as well as the possibility of elastic deformation relaxation during cooling of forging was taken into account.

Calculations were carried out for several options at different cooling rates for three stress components -axial, radial and tangential. Calculations of the relaxation effect were made by processing and approximating the experimental data of 40KhGM steel creep in a temperature range of cooling after anti-flake heat treatment.

Calculations showed that at analytical calculations without account of relaxation the absolute stress value decreases and tends to zero when approaching room temperature, whereas calculations with account of relaxation of internal stresses showed residual stresses at the end of forging cooling. Moreover, they are compressive on the surface and tensile in the center of the forging. The impact of these stresses on the flake disclosure was shown.

The second step was to analyze the influence of internal stresses on hydrogen diffusion and flake formation. The impact of external stresses on diffusion is usually limited to the decrease of the activation energy of diffusion to an additional value. Calculations have shown that this value has no effect on the activation energy.

The results of this article can be useful to specialists dealing with problems of flake formation in steel as well as to experts in computer modeling of technological processes in metallurgy.

Keywords: air cooling, internal pressure, stress relaxation, creep, stress analysis, simulation of cooling, forgings, flake formation, flakes, 40KhGM steel.

References

1. Dubovoy V.Ya. Flokeny v stalyakh [Flakes in Steels]. Moscow, Metallurgizdat Publ., 1950. 332 p.

2. Sklyuev P.V., Petrov B.D. [Dependence of Residual Stress on the Cooling Rate After Tempering]. Proble-my metallovedeniya i termoobrabotki [Problems of Metallurgy and Heat Treatment], 1960, vol. 2. pp. 136-142. (in Russ.)

3. Mirzaev D.A., Shaburov A.D. [Internal Stresses During Cooling of Large Forgings]. Bull. of the South Ural St. Univ. Ser. Metallurgy, 2012, no. 15 (274), issue 18, pp. 112-115. (in Russ.)

4. Tayts N.Yu. Tekhnologiya nagreva stali [Heating Technology of Steel]. Moscow, Metallurgizdat Publ., 1962. 567 p.

5. Sklyuev P.V. Termicheskaya obrabotka krupnykh pokovok [Heat Treatment of Large Forgings]. Moscow, Mashinostroenie Publ., 1976. 48 p.

6. Rabotnov Yu.N. Soprotivlenie materialov [Strength of Materials]. Moscow, Fizmatgiz Publ., 1962. 456 p.

7. Borzdyka A.M., Getsov L.B. Relaksatsiya napryazheniy v metallakh i splavakh [Stress Relaxation in Metals and Alloys]. Moscow, Metallurgiya Publ., 1978. 256 p.

8. Shtremel' M.A., Knyazev A.A. [Kinetics of Disclosure of Internal Grain Boundary Cracks by Hydrogen]. Fizika metallov i metallovedenie, 1986, vol. 62, no. 4, pp. 645-651. (in Russ.)

9. Shaburov A.D., Mirzaev D.A., Smolko V.A. [Heat Transfer During Cooling of Forgings in a Thermos]. Vestnik of Nosov Magnitogorsk State Technical University, 2011, no. 2, pp. 65-67. (in Russ.)

Received 5 September 2014