УДК 621.73
КОНЕЧНО-ЭЛЕМЕНТНАЯ МОДЕЛЬ РАСЧЕТА ВЕЛИЧИНЫ НАКОПЛЕННОЙ ДЕФОРМАЦИИ В ПРОЦЕССЕ ГОРЯЧЕЙ РАСКАТКИ КОЛЕЦ
© 2009 Ф.В. Гречников1, Е.В. Арышенский1, Э.Д. Беглов2
1 Самарский государственный аэрокосмический университет 2 ОАО "Самарский металлургический завод"
Поступила в редакцию 13.02.2009
Разработана, конечно-элементная модель расчета степени накопленной деформации на различных этапах деформирования кольцевой заготовки. Сопоставление результатов моделирования и экспериментальных зависимостей подтверждает адекватность модели.
Ключевые слова: раскатка колец, макроструктура, рекристаллизация, накопленная деформация, метод конечных элементов, модель, матрица жесткости, равнопрочные вставки.
В практике производства ГТД широко используются кольцевые детали, имеющие многофункциональное назначение. К этим деталям предъявляются высокие требования по структуре и уровню механических свойств. Основным способом получения кольцевых деталей является горячая раскатка (рис.1). Особенностью этого процесса является наличие многократных актов локальной деформации заготовки в момент её нахождения в валках и сопутствующая многократная частичная рекристаллизация в междеформационных, паузах затрудняющая расчет общей (накопленной) деформации за процесс.
Это приводит к тому, что по сечению заготовки могут одновременно присутствовать различные по величине, в том числе, и критические степени [1] деформации. В свою очередь, критические степени деформации способствуют образованию крупного зерна при окончательном рекристаллизационном отжиге. В тоже время в местах, где деформация превысила критические значения, будет образовываться мелкозернистая структура. Таким образом, неоднородность деформации приводит к разнозернистости, т.е структурной неоднородности по сечению деталей и снижению уровня механических свойств. Чтобы избежать этого, необходимо знать на каждом этапе величину накопленной деформации, полученной металлом как на каждом локальном этапе деформирования, так и за весь период раскатки в целом. В связи с этим целью данной статьи является построение математической модели, позволяющей определять напряженно-де-
Гречников Федор Васильевич, доктор технических наук, профессор, член-корреспондент РАН, проректор по учебной работе. E-mail: teach@ssau.ru. Арышенский Евгений Владимирович, аспирант. E-mail: omd@ssau.ru.
Беглов Эркин Джавдатович, кандидат технических наук, ведущий инженер. E-mail: beg@mail.radiant.ru.
формированное состояние и величину степени накопленной деформации.
При разработке конечно-элементной модели учтено, что, благодаря симметрии, структура и свойства раскатанного кольца идентичны для всех сечений по окружности. Учитывая это обстоятельство, модель была построена не для всего кольца, а для сегмента, равного 6-ти длинам очага деформации. Сегмент разбивается на треугольные конечные элементы, как показано на рис. 2.
Угол р, определяющий положение элемента в области решения, находим по следующей формуле [2].
_ 12 • 1 ■ Кг
(2ЯН + 2ЯВ) , (1)
где ЯН, ЯВ - наружный и внутренний радиусы кольца;
К - средний радиус кольца в 1 обороте.
Ь - длинна дуги контакта с любым из валков. Для её определения применена формула [2]
Ь 1(2) АН , (2)
Рис. 1. Схема процесса горячей раскатки колец: 1 - заготовка, 2 -внутренний неприводной валок (дорн), 3 - внешний приводной валок, 4, 5 -направляющие ролики, 6 - конечный выключатель (контроль диаметра)
где Я2 - радиусы приводного и не приводного валков
А Ь - абсолютное обжатие Предварительно разбиваем область решения на четырехугольные секторы, каждый из которых соответствует двум соседним треугольным элементам. Имеется N рядов секторов по радиальному направлению и М - в тангенциальном направлении. Имеется 2 ■ N ■ М треугольных элементов и (М + 1) ■ (N + 1) узлов. Нумерация узлов показана на рис. 2. Координаты 1-го узла по осям 1 и 2 обозначим как хц, X'2 •
д
х„ = 1
д - д
ВН + _Н вн
2 ■ N
ЧМ )]• НМММ )| ;<3>
= 1 ЕВн.+Дн-Дн • то!± ^тоД
2
2 ■ N
{М ]] \М ум
. (4)
В процессе расчета координаты узлов в любой точке области расчета будут меняться на пе-
ремещения узлов йп , 2 . Для нахождения йп , 2 воспользуемся энергетическим методом [3]. Рассмотрим отдельный треугольный элемент 1 с узлами 1, 2, 3 на рис 3.
Предположим, что элемент вначале не напряжен,_ узловые силы равны 0. Затем силы А , У , /3 прикладываются к соответствующим узлам элемента. Новая конфи-
гур й =
ция узлов будет иметь смещение й 11, й'12, й, й22, й^, й32 . Верхний индекс относится к элементу, в дальнейшем его опускаем. Первый нижний индекс относится к узлу, а второй - к координате. Потенциальная энергия I новой конфигурации по отношению к исходной представляет собой разницу между энергией напряженного состояния, накопленной в элементе и и работой совершенной силами /2,/3 на векторе перемещений ё, [3].
I=и-Ж=2 |(п + ст22£22+^^ Уь-А 1й11 -
, (5)
а)
Рис 3. Задание граничных условий в задаче о деформировании сегмента
2
х
где ё12.......- перемещения в узлах элемента
по направлениям 1,2 соответственно;
/п...... /32 - силы, под действием которых
происходит смещение узлов в направлении 1,2 соответственно;
е11 е22 - нормальные, а е12 - касательный компоненты тензора деформации;
у11у22 - нормальные, у12 - касательный компоненты тензора напряжений.
Интегрирование проводится по объему ^ (в рассматриваемом случае плоской деформации -по площади элемента dF). Для удобства дальнейшего решения представим уравнение (5) в матричной форме.
-21'
I = — |а -е-ёГ-ё 2
= I ^ 2
- Г = 2\еТЩеГ - =
•и - — , (6)
Значения компонент вектора ё = |ё„ ■■■ ё32|| должны быть такими, чтобы потенциальная энергия I имела минимальное значение:
д1 п
■— = 0 ; Н1...3, . (7)
мц
После дифференцирования, в векторной форме получим:
-И -ИНг )-ё = f. (8)
Чтобы понять обозначения, ||в||, и ||и|| еще раз рассмотрим отдельный элемент, представленный на рис.3.
Если он треугольный как в нашем случае, и напряжения в нем меняются линейно, то в [3] рекомендуется связывать значения перемещения узлов элемента и его деформацию следующей формулой.
«12
•22 Х32 0
••31 - Х21
0 Х32 Х12
••31 — Х21 0
•ц -31 0
Х22 Х-32 Х11 Х31 Х32 Х12 Х21 Х11
•21 Хц •12 •22
2
2
2
2
2
(9)
В матричной форме выражение (9) запишем следующим образом:
е = \\Б\\ - ё. (9 а)
Как видно из (9) ||в|| выражает изменения координат узлов треугольного элемента при сохранении его площади и связывает перемещение в его узлах с накопленной деформацией.
В свою очередь ||и|| выражает связь между тензором деформации и тензором напряжений. Его значения различны для упругого и пластического состояния. Вывод ||И|| для обоих состо-
яний можно найти в [3]. Здесь приведены его значения, причем только для плоской деформации и энергетического подхода. Упругая деформация:
а = И - е,
(10)
И =
1 - V V
1- - 2v 1 - 2v
Е V 1 -V
1 + V 1- - 2v 1 - 2v
0 0
(11)
Пластическое состояние:
ёа = (
И
И1
)- ёе = |И| - ёе , (12)
где И1
Ие
для упругой части деформации, для пластической части деформации.
И1 =
20 Б
1 -V V 0
1 - 2у 1 - 2у
20 V 1 -V 0
1 - 2у 1 - 2у
0 0 1 2
а11 а11 а11 0 22 ^ а11 012
а22 а11" 0 22 0 22 0 22 а12
а12 а11 а12 0 22 а12 012
(13)
(14)
где модуль сдвига О =
Е
(15)
2(1 + у)' У '
8 - характеристический параметр упруго-пластического состояния
22 Б = —а
3
(
1 + 0 30
V
(16)
у
Данный параметр позволяет учесть зависимости напряжений от деформации и других параметров процесса, которые выражены через соотношение вида
0 = 0(е,е, Т , а в с), (17)
где е -накопленная деформация при одноосном сжатии (растяжении);
е - скорость деформации; Т - температура;
аоа а,в,с - эмпирически определяемые соотношения. Поиску таких соотношений посвяще-
0
0
е
ё
11
ё
12
0
ё
21
ь —
X
ё
22
ё
2
31
ё
32
но большое количество исследований. Нами использованы результаты [4] для сплавов, используемых при раскатки колец ГТД.
Вернемся к формуле (8), которая, как теперь понятно, выражает связь между усилием в элементе, с одной стороны, и напряжением, деформацией и перемещением - с другой. Исключив из формулы (8) перемещения, обозначим её левую часть следующим образом.
Щ = М-|И-B-dF- (18)
Щ - это матрица жесткости. В ней учтены все параметры деформации, приведенные выше. Если данная матрица приведена для одного треугольного элемента - она называется локальной. Глобальная матрица будет представлять собой матрицу правой части системы (M ++1) уравнений, формируемую как алгебраическая сумма локальных матриц каждого элемента.
Щ-d = f . (19)
Надо отметить, что нам уже известно напряжение <s на первом шаге из свойств материала, а на следующих шагах - из формулы (12). Система (19) с нулевым вектором f имеет бесконечное множество решений, соответствующих движению объекта моделирования как жесткого тела. Для получения решения, описывающего движение реально деформируемого тела, необходимо задать граничные условия. С вычислительной точки зрения, положение граничных условий состоит в построении вектора правых частей уравнения (19). С физической точки зрения - в учете контактного взаимодействия заготовки с приводным и неприводным валками. Для уравнений, описывающих узлы на свободных поверхностях заготовки, свободный член равен 0. Рассмотрим различные подходы к построению граничных условий на приводном и неприводном валке.
Для неприводного валка в первой половине дуги захвата силы направлены против направления движения металла, во второй - по направлению движения (рис. 3, б). Для каждого узла в контакте с валком направление действия сил известно. P - нормальное давление, т = juP -сила трения, j - коэффициент трения.
Рассмотрим уравнение (19), которое в развернутом виде для узла 9 запишется следующим образом (рис. 3,б).
k17,17 d91 + k17,18 d 92 + k17,19 d101 + k17,20 d102 +
+ k17,21 d111 + k17,22 d112 = f91 =
= jP cos (p3 - P sin (p3, (20)
k18,17 d91 + k18,18 d92 + k18,19 d101 + k18,20 d102 +
+ k18,21 d111 + k18,22d112 = f92 =
= P sin (p3 + /uP cos (p3. (21)
При решении уравнений (20) методом Гаусса, учтем условие непроникновения материала заготовки в неприводной валок:
d91 ■ sin (р3 = d92 ■ cos^3. (22)
Это условие позволит исключить из системы уравнений (19) d92 Данноепреобразование проделываем для всех уравнений, содержащих узлы, лежащие на поверхности неприводного валка.
На приводном валке известна скорость вращения, но неизвестно взаимное смещение поверхностей металла и валка. Применим следующий прием [3].
Введем фиктивный слой элементов. Покажем его на примере элемента, с узлами 7, 6 (рис 3а). Эти узлы движутся как жестко связанные с валком. Узлы контактного слоя металла 5 (рис. 3 а) движутся по поверхности валка. Матрица жесткости элемента K модифицируется с помощью показателя трения m. Элементы матрицы жесткости умножаются на m/m - ц. При
m, стремящемся к 0, элемент делается более жестким, моделируя низкое трение. При m ^ 1 моделируется "прилипание" материала к валкам. Элементы не моделируют слой смазки, но моделируют действие смазки. Каждый элемент фиктивного слоя создается на момент построения соответствующего реального элемента. Матрицы реального и фиктивного элементов могут быть сопоставлены и совместно, решены в уравнении (8). Перемещения фиктивных узлов известны, т.е они движутся как жестко связанные с валком.
Уравнения (19) для узла 5 (рис. 3 а) будут иметь следующий вид.
k9 3d 23 + k 9,4d 22 + k9,7 d41 + k9,8 d42 + k9,9 d51 + + k 9,10 d52 + k 9,15 d 81 + k9,16 d82 + k 9,13 d71 + + k 9,14d 72 + k 9,11 d61 + k 9,12 d62 = f51 , (23)
k10,3 d 21 + k10,4d 22 + k10,7 d41 + k10,8 d42 + k10,9 d51 + + k10,10 d 52 + k10,15 d 81 + k10,16 d 82 + k10,13 d71 + + k10,14d72 + k10,11 d61 + k10,12d62 = f52 . (24)
Так как усилие в узле 5 нормально к поверхности валка, то имеем:
f2Cos^2 = fs1sin (Р2, (25)
Условие непроникновения поверхности валка ds1 cos^2 = ds2 sin (p2, (26)
При составлении глобальной матрицы жесткости, преобразуя уравнения (23, 24) с учетом (25,
1 т^',.
о, -
Рис. 4. Схема расположения равнопрочных вставок в очаге деформации при раскатке. Н0 - толщина заготовки до попадания в валки; у, х - значения координат вставки;
а0,Ь0 и ах,Ьх
начальные и конечные размеры вставок соответственно
52, йЪ1, также можно восполь-
26), исключая /51, /5 зоваться при решении системы (19) методом Га-уссового исключения. В ходе решения находятся значения накопленной деформации, напряжений и перемещений, т.е напряженно деформированное состояние в очаге деформации.
Проверка адекватности модели осуществляется на основе экспериментальных исследований раскатки колец, приведенных в работе [5]. В данной работе был исследован очаг деформации кольца из алюминиевого сплава АМг6, в котором по-
ei й
слойно сверлились отверстия и заполнялись вставками из того же металла (рис 4). Раскатка колец с внешним диаметром 400 мм, внутренним 340 мм и толщиной 30 мм осуществлялась на кольцераскат-ном стане модели РМ1200 с диаметрами рабочих валков: верхнего приводного - 550 мм и нижнего неприводного - 200 мм; максимальная скорость подачи нажимного устройства составляла 16 мм/ сек.; скорость прокатки, предусмотренная конструкцией стана, соответствовала 1,5 м/сек. По результатам измерения вставок находились значения
> *
"ч Т| / ) / [>
^^-- у-----
13
0,5
ш
■о/
О— й-
___^ С.ГЧС1 ИГ I /1^1111.1С
¿■¡и нт I а
---V №|ен.нч I данные
5веп;рскс т;
-*
анспсро-.'а та
СгУ 1 ,и инмь ь?
С:ч:'ини 2 ^ I члк МЗДСЛ.-фЭБаМН!
■И л -I л и и е. 2 т.я 11 и. 7ВДШ1 V ■ДЙМ [-1
Рис. 5. Распределение интенсивности деформации по высоте очага деформации при раскатке кольцевого образца из сплава АМг6: е1 - степень накопленной деформации, у - координаты точки по оси у (причем Но /2 соответствует на оси ординат 1)
деформаций и напряжений, которые представлены на рис. 5. Представленные экспериментальные данные по раскатке кольца из сплава АМг6 были введены в разработанную конечно-элементную модель. На рис. 5 сопоставлены результаты моделирования и экспериментальные данные.
Как видно из графика, результаты эксперимента и моделирования практически идентичны (сходимость около 15 %).
ВЫВОДЫ
1. Для формирования в кольцевых деталях ГТД однородной макроструктуры и требуемого уровня механических свойств необходимо контролировать величину накопленной степени деформации на каждом этапе горячей раскатки заготовки.
2. Разработана, конечно-элементная, модель рас-
чета степени накопленной деформации на различных этапах деформирования кольцевых заготовок.
3. Сопоставление результатов моделирования и экспериментальных зависимостей подтверждает адекватность модели.
СПИСОК ЛИТЕРАТУРЫ
1. Лахтин Ю.М., Леонтьева В.П. Металловедение. М.: Машиностроение, 1980. 493 с.
3. Целиков А.И. Теория расчета усилия в прокатных станах. - М.: Металлугргиздат, 1962.
2. Finite-element plasticity and metalforming analysis / G.W. Rove., C.E.N. Sturgess, P. Hartly., Cambridge University Press, 2005. 296 c.
4 П.И. Полухин, Г.Я Гун, А.М. Галкин Сопротивление пластической деформации металлов и сплавов. , М. Металлургия, 1983, стр. 353
5 Костышев В.А., Шитарев И.Л. Раскатка колец. -Самара: СГАУ, 2000. С. 206.
THE FINAL-ELEMENT MODEL CALCULATION SIZE SAVED DEFORMATION IN THE PROCESS OF HOT ROLLING RINGS
© 2009 F.V. Grechnikov1, E.V. Aryshensky1, E.D. Beglov2
1 Samara State Aerospace University 2 Open Societies "Samara metal works"
It is developed, is final-element model of calculation degree the saved up deformation at various stages of deformation of ring preparation. Comparison of results of modelling and experimental dependences confirms adequacy of model.
Key words: rolling rings, macrostructure, recrystallization, the saved up deformation, method of final elements, model, a rigidity matrix, full-strength inserts.
Fedor Grechnikov, Doctor of Technics, Professor, Corresponding Member of Russian Academy of Sciences, Vice Rector for Academic Affairs. E-mail: teach@ssau.ru. Evgenie Aryshensky, Graduate Student. E-mail: omd@ssau.ru.
Erkin Beglov, Candidate of Technics, Leading Engineer. E-mail: beg@mail.radiant.ru