Научная статья на тему 'Разработка математической модели распределения сил натяжения по ширине полосы при горячей прокатке в непрерывной группе стана 2800'

Разработка математической модели распределения сил натяжения по ширине полосы при горячей прокатке в непрерывной группе стана 2800 Текст научной статьи по специальности «Физика»

CC BY
279
68
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГОРЯЧАЯ ПРОКАТКА / НЕПРЕРЫВНАЯ ГРУППА / МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ / МЕТОД ГАЛЕРКИНА / МЕЖКЛЕТЬЕВОЕ НАТЯЖЕНИЕ / HOT ROLLING / CONTINUOUS GROUP / FINITE ELEMENT METHOD / GALERKIN METHOD / BETWEEN MILL TENSION

Аннотация научной статьи по физике, автор научной работы — Оводенко А. М., Арышенский Е. В., Гречникова А. Ф., Беглов Э. Д.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Оводенко А. М., Арышенский Е. В., Гречникова А. Ф., Беглов Э. Д.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

DEVELOPMENTS MATHEMATICAL MODEL DISTRIBUTION TENSION FORCES ON WIDHT OF STRIP AT HOT ROLLING IN CONTINUOUS GROUP ROLLING MILL 28001Joint-Stock Company "Alcoa SMZ", Samara

In thе article is developed mathematical finite element model hot rolling in continuous group qvarto mill 2800. For the help created model it is possible to calculate distribution of tension on width of a strip on an exit from each mill of continuous group. The decision of this problem allows to select technological parameters hot rolling process so that they provided set flatness sheet.

Текст научной работы на тему «Разработка математической модели распределения сил натяжения по ширине полосы при горячей прокатке в непрерывной группе стана 2800»

УДК 621.771.22/24:669.2

РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ РАСПРЕДЕЛЕНИЯ СИЛ НАТЯЖЕНИЯ ПО ШИРИНЕ ПОЛОСЫ ПРИ ГОРЯЧЕЙ ПРОКАТКЕ В НЕПРЕРЫВНОЙ ГРУППЕ СТАНА 2800

© 2011 А.М. Оводенко1, Е.В, Арышенский2, А.Ф. Гречникова1, Э.Д Беглов1

1 ЗАО "Алкоа СМЗ", г. Самара 2 Самарский государственный аэрокосмический университет имени академика С.П. Королева

(национально исследовательский университет)

Поступила в редакцию 14.04.2011

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

Ключевые слова: горячая прокатка, непрерывная группа, метод конечных элементов, метод Галерки-на, межклетьевое натяжение.

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

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

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

Оводенко Алексей Максимович, директор прокатного производства. E-mail: Alexey.ovodenko@alcoa.com Арышенский Евгений Владимирович, кандидат технических наук, ассистент кафедры «Обработка металлов давлением». E-mail: omd@ssau.ru

Гречникова Анна Федоровна, инженер-технолог прокатного производства. E-mail: Anna_gr17@yahoo.com Беглов Эркин Джадатович, кандидат технических наук, ведущий инженер прокатного производства. E-mail: beg@mail.radiant.ru

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

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

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

б

Рис. 1. Форма полосы: а - полоса планшетная; б - полоса, на которой наблюдается непланшетность

Рис. 2. Механизм образования непланшетности: а - полоса в очаге деформации; б - схема развития деформации;

1 - валки; 2 - полоса

ления продольных сил натяжения по ширине полосы. Другими словами, данные о распределении величины удельного натяжения по ширине полосы равносильны информации о планшетности.

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

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

Профиль валка зададим следующей функцией

Х3 = И(х1, х2) =

(х1 ~ Ь Д )2

Я

+ И1г + ЗИ( хх, Х2), (1)

ЬД =^(Я01 - Ин )* Я ;

Рис. 3. Пятиклетьевая непрерывная группа прокатного стана "Кварто 2800"

Рис. 4. Схема очага деформации при прокатке в непрерывной группе

5И(х1,х2) - величина упругой деформации валков;

Н0г, Н1г - заданные входная и выходная номинальные толщины.

Введем допущение, что скорость полосы по х3 можно задать линейно изменяющейся и равной у3=0 на плоскости симметрии х1Ох2. Находим, её по следующему выражению

(

vз =

дИ

дИ

л

чдх1 дх2

И

где Я - радиус рабочего валка;

Ь - длина недеформированной дуги контакта

(2)

^ у ,,

где у1, у2, у3 - скорость полосы на осях 1,2 и 3 соответственно

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

х1Ох21 1(х,х2Д . Для того, чтобы найти их, разобьем плоскость х1Ох2 на прямоугольные девя-тиузловые конечные элементы [3] (рис. 5).При этом скорости ^ (х1з х2) и у2(х1зх2) будут определяться во всех девяти узлах конечного элемента а давление у0 (х1,х2) только в четырех узлах (1, 2, 3, 4) (рис. 6.)

Как следует из метода конечных элементов [4] непрерывная функция в узлах области решения аппроксимируется по следующей формуле.

¥ = [кг, , Кк ...ДО, (3)

где ^ -значенияискомой величиныв узлах элемента;

N - функция формы элемента с соответствующим номером узла @ = ^ ^ к и.т.д.

В нашем случае функции формы элемента в локальных координатах для девяти узлов записываются следующим образом

N1 = -Ыл2 -л),

N2 = 4Ы2 +£)(л2 -л),

48

36

24

12

Рис. 5. Разбиение области решения на конечные элементы

Рис. 6. Девятиузловой конечный элемент ( Т, £ - локальные оси координат в конечном элементе)

1

Nз = ^+£)]+]),

N4 = -£)(] + ]),

N5 = 4(1 -е)т2+т) ,

N 6 = 4(£2 +£)(1 -]2),

N 7=-4(1 -е)т2+т) ,

N8 = V2-£)(1 -Т2),

Зададим уравнения для определения неизвестных функций. Это будут два уравнения равновесия по осям х1 и х2 и условие неразрывности. При этом исходим из допущения что проекцией давления металла на валок на ось х2 можно пренебречь:

д(спИ) И) р дИ + Ау

■ -----р—--г

дх1

дх2

дх1

22 + у2

■т = 0 ;(5)

д(^21И) + д(^22И) _V

дх1

дх2

^Ау* + У22 1

1 +

^ дИ ^

дх1

д(с! !И) , д(Иу2 )

дх1 дх2

=0

т = 0 ;(6)

(7)

N9 =

(4)

где с11, с22, с 2 - компоненты тензора напряжений, с1 3 = с23 = 0 ;

р - величина нормального давления металла на валок.

Аи1 = УВ - и1, УВ - линейная скорость валка.

г - напряжение трения касательное к поверхности валка.

Принимаем следующий закон трения: г = та,,,

где т - показатель трения, а, - предел текучести материала.

Граничные условия для уравнений (3), (4), (5)

а11 = '

для ура

, ч . '1 (х2); (8)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

' (х2) - равномерное натяжение (известное из начальных условий) на поверхностях ABFE и DCGH,

Для того чтобы воспользоваться уравнениями (5, 6, 7) в методе конечных элементов, необходимо привести их к уравнениям следующего вида:

| N рЩ)ёЯ = 0, р = 1, к..

(9)

' Н^+^х+л

1И) + д(а 2И) дк, А дх дх2 дх Л/Ачг+

ё(а21к) + 0(а22к)

дх

дх2 д/ац2 +и22 1

w2dxdx

(10)

В выражении (7) применяем интегрирование по частям. Для второго и третьего интегралов применяем формулу Грина [4] чтобы исключить производные по компонентам тензора напряжений. Для удобства записи выражение (10) представим как сумму I = 11 +12 +13 +14 . С помощью ввода в выражение (10) интеграла 15 учтем граничное условие а11 = '1 (х2),. Через Г обозначена область интегрирования, представляющая собой пересечение прямоугольника DCGH (плоскости выхода металла из очага деформации, на котором мы задаем базовое значение натяжения) с плоскостью х1Ох2. С учетом описан-

ных выше преобразований выражение (10) запишется следующим образом:

1 = 11 + 12 + 13 + 14 + 15

где 11 =я

д(а11 И ) д(ки 2)

дх,

д х

м> рАх 1 Ах2 ;

12 =-

дм

■ЯИ*^

2 У

дм2 дм 2 +- 1

дх дх.

2 У

дм2 дх2

I з = -Я

дИ

дх

А V!

1 л/Ац1

22 + У

м1ёх1ёх; (11)

14 =Я"Т=

V А и1

А

2 , 2 + Ц 2

У

1 +

( дИ л 2

чдх1 у

гм 2 Ах 1 Ах 2;

15 = Ц к1 (х 2 V1А Г

где R - область интегрирования;

Ь( у ) - дифференциальное уравнение, определяющие у.

Затем мы аппроксимируем функцию у в узлах конечного элемента с помощью формулы (3), где ^ - узловые значения искомой величины.

Для того чтобы привести уравнения (5, 6, 7) к дискретному виду, воспользуемся методом взвешенных невязок с весовыми функциями мр (х1, х2) - для уравнения неразрывности и М (х1, х2), м2 (х1, х2) - для уравнений равновесия. Через А обозначена область интегрирования, представляющая собой пересечение плоскостью х1Ох2 .очага деформации (фигура ABCDEFGH на рис. 2). Тогда интеграл взвешенных невязок запишется следующим образом

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

*11 =а0 а22 =а0 +Л822; а12 = Л812;(12)

л = ~~— *ед - интенсивность девиатора напряжений,

3Ве

• eq

интенсивность тензора скоростей де-

формации. Зависимость *ВС1 от 8ес1 может быть получена по результам лабораторных испытаний на осадку или кручение.

Компоненты тензора скоростей деформации выражаются через скорости

щёхёх +

дц

ди2

1

8ц =——; 822 = ——; 812 =821 = —

дх1 дх2 2

(дц_ а^

дх2 дх1

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

(

*33 = *0 + л

дИ

дИ

\

чдх1

дх

; (13)

(14)

2 У

Используя (12, 13 и 14) и предполагая р = -а33 интегралы 12 и 13 из выражения (11) можно переписать так:

12 =Я/*(81( К8^ +|Ы УУЙ,(м ( |м¥х

А А V '=1 ,=1 У

(15)

131Ш*

л

А А гг 1 (Эи эи

2

2 У

хёх (16)

А

т

А

V

2

Г

2 У

V

А

и1 +

V

2

Тогда сумма интегралов т выражается через неизвестные функции у1 (х15 х2), у2 (х15 х2) и щ (х1, х2) и соответствующие им весовые функции , w1, .

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Теперь уравнения вида (11) приведенные к виду (9) аппроксимируем в соответствии с уравнением (3):

V ■ V1 -к и

N 0 Ы2 ... 0 Ы9 0

0 N 0 N. 0 N

(17)

И

щ■( ш-N, щ, щ, щ)

И

04 I

(18)

. \ д\[ 0 д\2 дЧ 0

дх дх дх

«22 -ВТ 0 N дх 0 .. а. х 0 ач х 1 У2 1 Iе

«2 1 д\ 1 N 1N 1N 1 N 1 N й

1 X дх 1/2 ¿X л/2 дх л/2 дх; л/2 ¿Х _ Г9?

[0] [К„ ]

Кш] К ]

У

(20)

[К. ] = Я

■24(1 [в])+( к>)(5г и

(23)

ит =

= Ш01>"Ш0Яр>; (24)

[У]

(25)

Граничные условия (8) учтены интегралом по контуру г, с использованием криволинейных аппроксимирующих, функций составляющих вектор .

Чтобы составленная с помощью матрицы жесткости система уравнений имела решение введем еще одно допущение. Предполагаем, что на поверхности DCGH действуют только напряжения у11 = 1(х2) и у33.

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

(18.а)

Для вычисления интеграла 15 по поверхности Г аппроксимируем натяжение на поверхности DCGH следующим выражением:

г (Х2 )■ N, N В, NB)T(t*2 ). (19)

Теперь, когда функции представлены дискретными величинами согласно методу конечных элементов [4], построим матрицу жесткости.

И-

г(х2) 0

IX

0 0

0

00 0 И

(26)

Тогда условие пластичности и формулу для определения нормального всестороннего сжатия у0 для указанных поверхностей можно записать как

г(х2 )2 - г(х2 Из + И ■ ат 2 /3; (27)

Ш0 ■ 3И +Шз) ,

(28)

здесь шт предел текучести. Разрешая (27) относительно Ш33 и подставляя полученное значение в выражение (28) получаем

(х2 )- 1

(

где

2

г (х2 )"

4Щ + г (х 2 / 3

Л

(29)

г 'д^ т

[Кщ ] *')

Й N ] [В]

(21)

[Кщ]-я( фги\))т-[§ 1

N

Модифицируем систему уравнений (21) следующим образом: удаляем все уравнения соответствующие узлам щ 0, лежащим на границе DCGH, Столбцы матрицы [Кущ ], соответствующие этим узлам, умножаем на значения полученные по (29) и переносим в правую часть, модифицируя век-

тор (У I. Таким образом, мы определяем силовые

условия в узлах, примыкающих к DCGH и система уравнений становится разрешимой.

2

0

2

0

т

е

А

2

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

Ф,

t 1V> tN 1 = v\ - j ,

(30)

где N - количество клетей непрерывного стана;

V1 - выходная скорость в узле , клети ' стана;

у°+1 , - входная скорость в узле , клети ' +1 ' - вектор узловых значений переднего натяжения для клети г

Искомым распределением натяжений будет такое при котором

Ф

i

t,

V V

,r р.. V N

= О . (31)

NeB

Z Í

1 г

Nb

t)dx2 = Ti.

(32)

5Ф / \ 5Ф..

- V,0)+£-

dVk

+ Ф..

10 10

1N

dt,

V 01...V0

(k -1° )+

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

шб т / ч

(-('°)}ь2 = 0,

q г

где к - номер приближения

ц - номер сечения в котором находится значение натяжения.

Частные производные вычисляем как конечные разности, для этого решаем систему (18) при малых изменениях только одного из узловых значений натяжений. Расчет заканчиваем по достижению сходимости по модулю вектора разности (')-и

Кроме того, интегралы удельного натяжения по поверхности DCGH должны быть равны полному межклетьевому натяжению Т', заданному как исходное данное в условиях задачи, в соответствии с (21)

Число неизвестных в уравнениях (31),(32) равно количеству узлов в которых заданы натяжения, во всех клетях плюс количество клетей. Распределение натяжений, которое обеспечивает выполнение все этих равенств является решением при форме межвалкового зазора заданного выражением (1). Система уравнений (31) (32) нелинейная. Для ее решения применим последовательные приближения. В качестве начального приближения можно принять равномерное распределение натяжения. Допустим, что на очередном шаге вычислений распределение натяжений по клетям задано последователь-

А0 -0 \

ностью векторов 1,....,tN >и скоростей

(у 10...У^М. Тогда следующее приближение найдем из линеаризованной формы (31), (32):

После каждого решения (18) производим пересчет упругих деформаций валков 5h( x1, x2), пользуясь известными методиками и аналитическими решениями для упругих деформаций полубесконечного тела под действием поверхностных нагрузок [5]. Таким образом, с помощью полученной математической модели можно определить распределение значений натяжений и давлений метала на валок по ширине прокатываемой полосы. Алгоритм реализации математической модели представлен на рис. 7.

ВЫВОДЫ

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

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

3. Предложен алгоритм компьютерной реализации математической модели.

СПИСОК ЛИТЕРАТУРЫ

1. Настройка, стабилизация и контроль процесса тонколистовой прокатки / Григорян Г.Г [и др.]. М.: Металлургия, 1975. 368 с.

2. Управление качеством тонколистового проката / В Л. Мазур [и др.]. Киев: Техника, 1997. 385 с.

3. P.M. Dixit, Dixit U.S. Modeling of metal forming and machining processes : by finite element method and soft computing . Springer Verlag London Limited , 2008. 590 p.

4. Сегерлинд Л. Применение метода конечных элементов [пер. с англ.]. М.: Мир, 1979. 392 с.

5. Формоизменение листового металла / Баркая Ф.Е [и др.]. М.: Металлургия, 1976. 260 с.

= 0;

(33)

t

T

k

kl

N

Рис. 7. Алгоритм реализации математической модели

DEVELOPMENTS MATHEMATICAL MODEL DISTRIBUTION TENSION FORCES ON WIDHT OF STRIP AT HOT ROLLING IN CONTINUOUS GROUP ROLLING MILL 2800

© 2011 A.M. Ovodenko1, E.V. Aryshenskij2, A.F. Grechnikova1, E.D Beglov1

'Joint-Stock Company "Alcoa SMZ", Samara 2 Samara State Aerospace University named after Academician S.P. Korolev (National Research University)

In the article is developed mathematical finite element model hot rolling in continuous group qvarto mill 2800. For the help created model it is possible to calculate distribution of tension on width of a strip on an exit from each mill of continuous group. The decision of this problem allows to select technological parameters hot rolling process so that they provided set flatness sheet.

Key words: hot rolling, continuous group, finite element method, Galerkin method, between mill tension.

Aleksey Ovodenko, Director of Rolling Division. E-mail: Alexey.ovodenko@alcoa.com

Evgenie Aryshenskij, Candidate of Technics, Assistant Lecturer at the Metal Working Department. E-mail: omd@ssau.ru Anna Grechnikova, Engineer-Technologist. E-mail: Anna_gr17@yahoo.com

Beglov Erkin Djavdatovich Candidate of Technics, Leading Engineer. E-mail: beg@mail.radiant.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.