CC BY
83
УДК 539.37: 621.771.07
Е. Б. Бочектуева
ЧИСЛЕННОЕ ОПРЕДЕЛЕНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ В ВАЛКАХ И УСИЛИЙ ПРОТИВОИЗГИБА В ЧЕТЫРЕХВАЛКОВОЙ КЛЕТИ ПРОКАТНОГО СТАНА
Разработаны расчетные методы и созданы программные средства для численного определения НДС и усилий противоизгиба холоднолистового четырехвалкового прокатного стана с учетом остаточных термонапряжений. Расчет термонапряжений осуществлен посредством решения задачи термоупруговязкопластично-сти для среды с нестабильной структурой. В основу решения объемной контактной задачи для определения рабочих напряжений от прокатки положен конечно-элементный комплекс ANSYS. Представлены результаты численной реализации применительно к реальному прокатному стану КВАРТО500/1300x1700. Установлено, что приближенное аналитическое решение, основанное на балочной схеме, дает заниженные значения усилий противоизгиба, а также что остаточные напряжения от термообработки валков значительно больше эксплуатационных напряжений от прокатки.
E-mail: [email protected]
Ключевые слова: листопрокатный стан, клеть кварто, усилия противоизгиба, объемная контактная задача, метод конечных элементов, остаточные термические напряжения, теплопроводность, термоупруговяз-копластичность, нестационарная структура.
Важной проблемой при изготовлении листового проката является обеспечение постоянности его толщины. Разнотолщинность листа возникает вследствие изгиба рабочих валков в процессе прокатки. В четы-рехвалковых станах (кварто) для уменьшения изгиба рабочих валков используются опорные валки, на которые передается основная часть усилия прокатки [1], кроме того, к рабочим валкам прикладываются усилия противоизгиба, выгибающие их в противоположном направлении и выравнивающие стрелу прогиба. Как правило, усилия противоизгиба выбирают посредством прикидочных расчетов и ориентировочных значений, полученных в процессе эксплуатации аналогичных прокатных станов.
В настоящей работе разработана основанная на МКЭ математическая модель прокатного стана кварто, позволяющая численно определять НДС в валках и усилия противоизгиба. Следует отметить, что при одних и тех же размерах валков и усилиях прокатки НДС в них будут различными в зависимости от того, какой термической обработке подвергались валки в процессе изготовления и какие остаточные напряжения сохранились в них к началу эксплуатации. Остаточные термонапряжения, как правило, превышают эксплуатационные напряжения,
Рис. 1. Схема нижней половины стана КВАРТО 600/1500x1700:
1 — прокатываемый лист; 2, 3 — рабочий и опорный валки; Q — усилие противоизгиба
возникающие при прокатке [2], поэтому для адекватного определения усилий противоизгиба необходимо учитывать в расчете значения этих остаточных напряжений.
Для обоснования целесообразности применения разработанной математической модели были сопоставлены значения усилий противо-изгиба, полученные путем численного эксперимента, и значения, найденные по традиционной методике [1], основанной на балочной схеме. В качестве объекта исследования выбран холоднолистовой четы-рехвалковый прокатный стан КВАРТО 600/1500х1700, состоящий из двух рабочих (сталь 90ХФ, диаметр 0,6 м, длина рабочей части 1,7 м) и двух опорных валков (сталь 90Х, диаметр 1,5 м, длина рабочей части 1,7 м). Усилие прокатки принималось равным 20 МН, а ширина прокатываемого листа — 1,4 м. Остальные геометрические размеры, используемые в расчете, приведены на рис. 1. В силу симметрии задачи относительно прокатываемого листа на рисунке изображена только нижняя часть клети.
Согласно традиционной методике [1] рабочий валок рассматривается как балка с постоянным поперечным сечением. Дифференциальное уравнение упругой линии у(г) балки в общем случае имеет вид
Е1хУ1У (г) = f (г), (1)
где f (г) — функция изменения поперечной погонной нагрузки по длине балки; Е — модуль Юнга; 1х — момент инерции поперечного сечения.
Общее решение этого дифференциального уравнения складывается из общего решения однородного уравнения при f (г) = 0 и частного
Рис. 2. Силовая схема рабочего валка решения неоднородного уравнения [3]:
EIx v(z) = Со + Cl z + c^— + C3 — +
2 6
(z - t)3 6
f (z) dt, (2)
где co , ci , c2 , c3 — постоянные интегрирования, определяемые граничными условиями.
Из силовой схемы, приведенной на рис.2, следует, что функция f (z) для рабочего валка имеет следующий вид:
Г 0, l < | z | < d;
f (z) = { -p(y), b< | z | < l; (3)
[ q(z) - p(z), 0 < | z | < b.
Рассмотрим граничные условия задачи. Симметрия рабочего валка относительно срединной плоскости приводит к тому, что в начале координат угол поворота равен нулю (v'(0) = 0); изгибающий момент в точках приложения усилий противоизгиба также равен нулю (v''(±d) = 0), а поперечные силы — равны этим усилиям (EJxv//'(±d) = Q). С учетом граничных условий выражение для прогиба валка (2) принимает вид
2 3 0 2 0
EJxv(z) = Co + Qz2d + Qj + j z^p(z ) dt + y q(z) dt+
-l -b
0,
+ <
(z -1) 6
l < | z | < d; p(z) dt, b < | z | < l;
-i
(4)
(z - t)3 . , , ? (z - t)3 . . , v ; -q(z) dt - v ^ p(z) dt, 0 < | z | < b.
6
6
-b
-i
z
2
3
Прогиб в уравнении (4) определен с точностью до постоянной интегрирования с0 — прогиба в начале координат. Последнее связано с тем, что при определении общего решения дифференциального уравнения четвертого порядка (1) имеется возможность записать только три граничных условия. Принимая для простоты давление ( прокатки и контактное давление р между рабочим и опорным валками постоянными, получаем формулу для определения прогиба:
, , й г2I2 г2 Ь2 г4 г4 Е1хь(г) = со + ф—---—р + — ( - —Р + ^^(-
- ^(Я^2^(Н(!), (5)
где Н(Ь), Н(I) — функции Хевисайда.
Разность прогибов рабочего валка посередине и у кромки листа с учетом того, что р1 = (Ь + ф, будет равна
т^тьт^ I4 (ЬУ( Ь2 7 Ь Л (Ь\2 ( Ь2 й 1 ВДО, Ь)=((- ^ -6У +1) + ф13{у) (йт - 2! + 4
(6)
Очевидно, что для равномерной толщины прокатываемого листа величина Ау в последнем выражении должна равняться нулю. Приравняв левую часть уравнения (6) к нулю, можно получить формулу для определения усилия противоизгиба:
^ Ь (7Ь Ь2 Л /( Ь2 й 1 \ ф = (4 \ 61 - 6Р - ^/{ш - 2 + 4) ■ (7)
Давление прокатки ( определяется как отношение усилия прокатки к длине рабочей части валка. Подставляя геометрические размеры Ь, I и й, равные согласно рис. 1 соответственно 0,7; 0,85 и 1,350 м, находим значение усилия противоизгиба 891 кН, что составляет 4,46% усилия прокатки. Из практики известно, что суммарное усилие противоизгиба (2ф) обычно составляет 10... 20% усилия прокатки, т.е. применение формулы (7) приводит к заниженным результатам. Это объясняется тем, что отношение длины рассматриваемого рабочего валка к диаметру бочки приблизительно равно 3 и определять усилия противоизгиба на основе балочной схемы можно только в первом приближении. Кроме того, остаточные напряжения в валках оказывают заметное влияние на НДС валков и, следовательно, на значения усилий противоизгиба.
Остановимся на численном методе расчета остаточных напряжений от термообработки. Для вычисления термонапряжений фактически нужно решить нелинейную нестационарную задачу теплопроводности, выполнить моделирование структурообразования при термообработке и собственно определить напряжения. Поскольку теплофизи-ческие коэффициенты зависят от температуры и структуры, кинетика
структурных превращений зависит от изменения температуры, а параметры уравнений состояния при пластичности и ползучести зависят от температуры и структуры, то все три задачи оказываются взаимосвязанными. Термообработку крупногабаритных прокатных валков обычно выполняют таким образом, что на центральные области валков в течение длительного времени действуют высокие температуры, при которых проявляются реономные свойства стали. В поверхностных слоях вследствие интенсивного охлаждения формируются пластические деформации. Поэтому для адекватного описания деформирования стали при термообработке необходимо использовать модель упруговязкопластической среды с нестабильной структурой [4].
Для решения нелинейной нестационарной задачи теплопроводности с граничными условиями третьего рода применяли МКЭ. Для осе-симметричного изотропного тела в случае переменных теплофизиче-ских коэффициентов эта задача описывается следующим дифференциальным уравнением [5]:
срд1 = 1 • д (Хг-^ + — (а—^ + (8)
дт г дг \ дг) дг \ дг) '
где ¿(г, г, т) — температура; г, г — координаты по радиусу и по оси соответственно; т — время; с — теплоемкость; А — коэффициент теплопроводности; р — плотность; (у — мощность удельных источников энерговыделения. Для аппроксимации в уравнении (8) производной по времени использовали безусловно устойчивую центральную разностную схему Кранка-Никольсона. В расчете принимали, что коэффициент теплопроводности зависит от температуры и структуры, а теплоемкость — только от температуры; тепловыделения при структурных превращениях учитывали посредством слагаемого — мощности удельных источников энерговыделения, включенного в уравнение теплопроводности.
При моделировании структурообразования при термообработке стали для описания изотермического распада аустенита в перлит использовали уравнение Авраами [5], а для бейнитного превращения — более сложное уравнение, учитывающее ускорение хода превращения и снижение предельной степени распада аустенита под нагрузкой. Переход от изотермической кинетики превращения аустенита к неизотермическим условиям проводили согласно теории изокинетических реакций на основе правила аддитивности. Уравнения для определения удельной доли новой структуры на п-м шаге имеют вид [4]
1п(1 - уТУУа)
K (tn)
1/n(tn)
КЫ = {1 - exp [-K(tn) (тп + Ar„)n(in)] } ^аб, (9)
где тп — приведенное время, которому соответствует удельная доля новой структуры, накопленная на предыдущем шаге в предположении, что превращение протекает при температуре данного шага £п; РА — удельная доля аустенита, сохранившегося к началу бейнитно-го превращения. Для перлитного превращения этот параметр равен единице.
Напряжения определяли, решая задачу термоупруговязкопластич-ности для материала с нестабильной структурой [4]. В основу решения был положен пошаговый метод определения дополнительных (начальных) деформаций [7]. В этом случае приращение тензора полной деформации можно представить в виде суммы
= Д4 + д4 + Д4 + Дет, (10)
где , Дер,-, Де^ — приращение упругой, пластической и деформации ползучести соответственно; ДеТ — приращение свободной деформации, учитывающей температурные и структурные изменения объема; 8^ — символ Кронекера.
Согласно методу дополнительных деформаций решение задачи термоупруговязкопластичности сводится к последовательному решению задачи термоупругости, при этом три последних слагаемых в уравнении (10) объединяются в одно, т.е.
= Д4 + Д4, (11)
где Де0 = ДеР + Де^ + 8^ Дет — дополнительная деформация.
Принцип расчета термонапряжений на каждом временном шаге заключался в итерационном уточнении приращений деформаций пластичности и ползучести, включенных в состав дополнительной деформации, при которых интенсивность напряжений, вычисленная в задаче термоупругости, не превышала бы значения, полученного по кривой растяжения для накопленной к данному шагу пластической деформации. Задачу термоупругости так же, как и температурную задачу, решали МКЭ. Изолинии наиболее опасных остаточных осевых напряжений в рабочем валке после индукционной закалки токами промышленной частоты (ТПЧ) с отпуском приведены на рис. 3.
При создании математической модели изгиба валков в процессе прокатки в целях снижения числа используемых конечных элементов учитывали симметрию задачи относительно продольного сечения и срединной плоскости (рис. 4). Для определения НДС в валках и выбора на основе условия постоянности прогиба рабочего валка усилий противоизгиба решали задачу о контактном взаимодействии двух цилиндрических поверхностей. Для этого в зоне контакта рабочего и опорного валков были введены специальные контактные конечные элементы.
Рис. 3. Изолинии остаточных осевых напряжений (МПа) в рабочем валке после индукционной закалки ТПЧ с отпуском
Математическая модель была создана в среде А^УБ 11.0. Конечно-элементную сетку формировали посредством "протягивания" плоской фигуры, разбитой на плоские четырехугольные 8-узловые вспомогательные конечные элементы МЕБИ200, вдоль продольной оси валков. Эта фигура имела внутренние области. Внешняя граница фигуры представляла собой поперечное сечение рабочего и опорного валков, внутренние области — поперечные сечения шеек; имелась также область более мелкого разбиения в зоне контакта валков. Для постепенного перехода от зоны мелкого разбиения к зоне обычного разбиения дополнительно задавали промежуточную область. В результате процедуры "протягивания" получали конечно-элементную модель, состоящую из 20-узловых объемных (3Б) элементов 80ЬГО95. Данный конечный элемент позволяет без потери точности использовать нерегулярную форму сетки, а также описывать модели с искривленными границами. Поверхности рабочего и опорного валков в зоне полосы контакта моделировали 8-узловыми поверхностными трехмерными контактными С0МТА174 и ответными ТАЯ0Е170 элементами, предназначенными для моделирования контактного взаимодействия между трехмерными деформируемыми поверхностями. Причем эти элементы располагали на поверхностях трехмерных объемных элементов 80ЬГО95. Для учета остаточных радиальных (аг), осевых (ог), окружных (а г) и касательных (тгг) напря-
Рис. 4. Конечно-элементная модель клети стана кварто
жений при решении контактной задачи эти напряжения задавали в А^УБ как начальные посредством использования функции inis.
В процессе численного счета задавали различные значения усилий противоизгиба и определяли погрешность толщины листа, которую рассчитывали по формуле
6 = 2А-Д, (12)
где А- — максимальная разность радиальных перемещений узлов, расположенных на рабочей поверхности валка, в зоне контакта с прокатываемым металлом; Ь — толщина листа (3 мм).
Вначале задавали усилие противоизгиба, рассчитанное по формуле (7). Результаты численного расчета показали, что разнотолщинность листа в этом случае составляет ~3 %. Далее увеличивали усилия про-тивоизгиба и анализировали значение разнотолщинности. Численный расчет показал, что наименьшая разнотолщинность листа возникает при усилиях противоизгиба, составляющих 7 % усилия прокатки. На рис. 5 приведены изолинии наиболее опасных осевых напряжений для этого случая. Из сопоставления рис. 3 и 5 следует, что остаточные термонапряжения в валках значительно больше эксплуатационных напряжений от прокатки.
Выводы. 1. Приближенное аналитическое решение, основанное на балочной схеме, дает существенно заниженные значения усилий противоизгиба.
Рис. 5. Изолинии суммарных осевых напряжений (МПа) от термообработки и прокатки в рабочем валке с рациональными значениями усилий противо-изгиба
2. Остаточные напряжения от термообработки валков значительно больше эксплуатационных напряжений от прокатки.
3. Корректное определение усилий противоизгиба в четырехвалко-вых клетях прокатных станов возможно только численными методами, при этом необходимо учитывать в расчете остаточные напряжения от термообработки.
Работа выполнена в рамках гос. контракта № 02.513.11.3487 от 8 октября 2009 г.
СПИСОК ЛИТЕРАТУРЫ
1. Целиков А. И., Никитин Г. С., Рокотян C. E. Теория продольной прокатки. - М.: Металлургия, 1980. - 319 с.
2. Покровский А. М. Оценка ресурса прокатных валков с учетом остаточных напряжений от термической обработки // Производство проката. - 2005. - № 9. -С. 26-31.
3. Смирнов В. И. Курс высшей математики. - М.: Физматгиз. - Т. II, 1958. -628 с.
4. Вафин Р. К., Покровский А. М., Лешковцев В. Г. Прочность термообрабатываемых прокатных валков. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. - 264 с.
5. К о з д о б а Л. А. Методы решения нелинейных задач теплопроводности. - М.: Наука, 1975.-228 с.
6. КристианДж. Теория превращений в металлах и сплавах: Пер. с англ. М.: Мир. -Ч. 1: Термодинамика и общая кинетическая теория, 1978. - 808 с.
7. Термопрочность деталей машин / Пол ред. И.А. Биргера, Б.Ф. Шорра. - М.: Машиностроение, 1975. - 456 с.
Статья поступила в редакцию 16.10.2009
Елена Баторовна Бочектуева окончила Томский государственный университет в 2003 г. Аспирантка кафедры "Сопротивление материалов" Восточно-Сибирского государственного технологического университета (Улан-Удэ). Автор четырех научных работ в области численных методов оценки прочности машин.
Ye.B. Bochektueva graduated from the Tomsk State University in 2003. Post-graduate of "Strength of Materials" department of the East-Siberian State Technological University (Ulan-Ude). Author of 4 publications in the field of numerical methods of machine strength estimation.
