CC BY
28
УДК 621.771.01
А.В. Сатонин, С.С. Настоящая, В.Г. Переходченко A.V. Satonin, S.S. Nastoyashaya, V.G. Perehodchenko
Донбасская государственная машиностроительная академия, г. Краматорск, Украина Donbass State Engineering Academy, Kramatorsk, Ukraine
А.Г. Присяжный A.G. Prisyagniy
Приазовский государственный технический университет, г. Мариуполь, Украина Azov State Technical University, Mariupol, Ukraine
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ТОНКИХ ПОЛОС ПРИ ПРОКАТКЕ
MATHEMATICAL SIMULATION OF STRESS-STRAIN STATE THIN STRIP IN THE ROLLING
На основе рекуррентного решения конечно-разностной формы условий баланса энергетических затрат, а также условия статико-динамического равновесия, рассматриваемых в рамках выделенных элементарных объемов, полученных путем разбиения зоны пластического формоизменения по ее длине, с использованием более корректных форм аналитических описаний условий внешнего контактного трения и условия пластичности, получили развитие численные математические модели локальных и интегральных характеристик напряженно-деформированного состояния металла при горячей прокатке относительно тонких лент и полос в чистовых рабочих клетях штрипсовых и широкополосных станов.
Ключевые слова: широкополосный стан, тонкая полоса, горячая прокатка, холодная прокатка, математическая модель, локальные и интегральные характеристики напряженно-деформированного состояния.
On the basis of recurrent decision finite difference form of terms balance of power expenses and the terms of static and dynamic balance, examined within the framework of the selected elementary volumes, got by breaking up of area of plastic forming on its length, with the use of more correct forms of analytical definition of terms of external contact friction and condition of plasticity, the numeral mathematical models of local and integral adjectives of the deflected mode of metal at the hot rolling in relation to thin strips and bars got development in the finishing working train skelp and wide-strip rolling mill.
Keywords: wide-strip rolling mill, thin bar, hot rolling, cold rolling, mathematical model, local and integral adjectives of the deflected mode.
Повышение требований к объемам и степени достоверности результатов математического моделирования различных технологических схем прокатного производства делает необходимым уточнение количественных оценок, а также использование достаточно строгих аналитических описаний граничных условий очага деформации, в качестве которых следует рассматривать механические свойства прокатываемых металлов, условия внешнего трения на контактных поверхностях рабочих валков и геометрические параметры зоны упругопластического формоизменения. Применительно к процессам прокатки относительно тонких полос в настоящее время достаточно широкое распространение получили основанные на конечно-разностных подходах численные математические модели, представленные в работах [1-3] и ряде других.
С целью расширения объемов и повышения степени достоверности предоставляемой информации разработана математическая модель напряженно-деформированного состояния металла при горячей прокатке относительно тонких полос. Непосредственно математическое моделирование заключалось в разбиении очага деформации на конечное n-е множество i-х элементарных объемов и в последующем рекуррентном решении конечно-разностной формы баланса энергетических затрат, рассматриваемого в рамках каждого из них. На рис. 1 представлена расчетная схема интегрального очага деформации. Зона пластического формоизменения металла протяженностью Ьпл была разбита на зону отставания протяженностью Ьот и зону опережения протяженностью Lon. В состав интегрального очага деформации входила и зона упругого восстановления прокатываемой полосы протяженностью ^п.
По аналогии с методикой работы [3] осуществлялось разбиение суммарного угла контакта ф02, определяемого как сумма углов контакта с зоной пластического формоизменения ф01 и зоной упругого восстановления, заведомо превышающих реальные значения данных углов, т.е.
ф02=ф01 + ф02, ф01 = arcsin (1,5л/RAH -Ah2 /4/R), ф02 = 0,3ф01.
Разбиение на конечное n-е множество i-х элементарных объемов, имеющих угловые характеристики для начального фа и для конечных ф, 2 граничных сечений осуществляли по следующей схеме решения (см. рис. 1, б):
Аф = фо2 / n, ф,1 = ф01 - Аф(, -1), ф, 2 = ф,1 - Аф = ф01 - Аф,,
а геометрические координаты, имеющие свое начало в вертикальной плоскости параллельной осям рабочих валков и соответствующие углу ф02, определяли по формулам
хп = R sin фй,
Xi 2 = R sin Фі 2.
Для первой итерационной процедуры рабочие валки принимались абсолютно жесткими, исходя из чего определялись текущие значения межвалко-вого зазора для начального hxi1 и конечного hxi2 граничных сечений каждого отдельного выделенного i-го элементарного объема:
hxi1 = h1 + 2R I1 - COs Фі1 )
Ki2 = h1 + 2R (1 - COs Фі2 ).
б в г
Рис. 1. Расчетная схема интегрального очага деформации применительно к численному конечно-разностному математическому моделированию напряженно-деформированного состояния металла при горячей прокатке относительно тонких полос (а), а также расчетные схемы элементарных объемов, выделенных в зоне упругого восстановления (б), в зоне опережения (в) и в зоне отставания (г) очага деформации
В рамках проведенного численного математического моделирования были приняты следующие допущения:
- пластическая деформация прокатываемой полосы является плоской и установившейся во времени;
- кинематика пластического течения металла в очаге деформации подчиняется гипотезе плоских сечений [4], а все граничные сечения очага деформации являются вертикальными;
- нормальные осевые напряжения ax и показатели удвоенного сопротивления сдвигу прокатываемого металла 2Kxi изменяются только по длине очага деформации;
- изменение текущих значений толщин hxi = hxi1...hxi 2, нормальных контактных pxi = pxi1...pxi2 и касательных контактных напряжений Txi =тх1...тх2 по длине каждого отдельного выделенного i-го элементарного объема линейно;
- аналитическое описание касательных контактных напряжений подчиняются закону трения Леванова и Колмогорова [5], в соответствии с чем ее численная интерпретация имеет вид:
Т»1 = 2 Kx# xi1 f1 - exp (-1,25 pX(i-1)1/a (.-1)1)
Txi 2 = 2Kx 2^ x 2 f1 - eXP (-1,25 p,i1/ °si1 ) (1)
где 2Kxl, 2Kxi 2 - удвоенные значения сопротивления сдвигу прокатываемого металла в начальном и конечном граничных сечениях выделенного i-го элементарного объема [5]; Vxl(2) - скорости перемещения прокатываемого металла в начальном и конечном граничных сечениях выделенного i-го элементарного объема [6], Vxl(2) = V1h1/hxl(2); pxi1, asi1 - нормальные контактные напряжения и напряжения текучести прокатываемого металла для начального граничного сечения i-го элементарного объема; ц xi 2 - текущее по длине очага деформации значение коэффициента пластического трения, определяемое в зависимости от геометрической координаты xi2 дифференцированно для зоны отставания и зоны опережения.
Принимая за основу полную форму записи условия пластичности [4] при плоской деформации и учитывая аналитическое описание величины касательных контактных напряжений Tx 2 (1) по отношению к нормальным осевым напряжениям ax 2, можно получить
a xi 2 = pxi 2 - ^ 4Kxi 2 -{{xi2 12Kxi2 f1 - exP ( 25 p^ / aSxi1 )]} =
= Рх2 - 2 К, 2ф-{{ [1 - ехр (-1,444рА/2КхЛ )
Рм 2 2 Кх12 ацхг 2 ,
ГДе ^2 = ^1 -{{2С1 - еХР(-1,444Рх11/2КхП )]} -
- вспомогательная пере-
менная, используемая для упрощения дальнейшей формы записи.
Приняв за основу уравнение баланса энергетических затрат [3], получим
где положительные значения нормальных осевых напряжений ам1 и ам 2 соответствуют напряжениям сжатия; Ах, - шаг разбиения очага деформации, Ах, = хг1 - хг2; Ув - окружная скорость вращения рабочих валков; ам - текущее значение угла контакта на рабочем валке, величина которого весьма незначительна, 008 ах = 1,0; КЛ - коэффициент немонотонности пластической деформации; рм, ах - плотность металла прокатываемой полосы и величина его ускорения, имеющего место в рамках выделенного г-го элементарного объема [3].
На основе зависимости (2) после соответствующих математических преобразований по отношению к искомой величине нормальных контактных напряжений рх12 в окончательном виде получим
Следует указать, что решение, аналогичное (3), может быть получено не только исходя из условия баланса энергетических затрат, но и исходя из условия статико-динамического равновесия всех сил, действующих в рамках выделенного г-го элементарного объема (см. рис. 1, в, г) [4]. Трансформированная аналитическая форма записи данного условия по отношению к искомой величине нормальных контактных напряжений рг 2 имеет вид
+
(2)
рхг2 ={ах;1 + 2Кхг 2а^г2 + 0,5 (т(1 / Ки + Тхг2 / 2 )Ах;
- 0,5КЛ(2 Кхп + 2 КХ12) (1/К 2 )-
- 0,25Рм (1 + НХ12(1/-1/^1) /)2}.
(3)
Рхі2 = {2[Схі1Кхі1 + 2Кхі22Кхі2 + рхі1 + Тхі2 ))і ] + Рхі1 (хі1 Кхі2 )
- рм (ж1 + Кхі2 ) ^і (Жі2 - УжіІ ) / ( 2Дхі )} / (ж1 + Кхі2 ) •
Последующий расчет напряженного состояния металла в зоне упругого восстановления прокатываемой полосы осуществляли на основе численного рекуррентного решения конечно-разностной формы условия статического равновесия выделенного г-го элементарного объема (см. рис. 1, б):
где величина нормальных контактных напряжений для конечного граничного сечения
где юп, Еп - коэффициент Пуассона и модуль упругости прокатывамой полосы.
Решив уравнение (4) относительно нормальных осевых напряжений ах! 2, получим основное уравнение рекуррентной схемы решения напряженного состояния металла в зоне упругого восстановления прокатываемой полосы:
По мере определения реального значения протяженности зоны опережения Ьоп (см. рис. 1, а) и соответствующих ей значений всех локальных характеристик напряженного состояния металла производили организацию еще одной итерационной процедуры, обеспечивающей учет упругого сплющивания рабочих валков. С этой целью на основе численной интерпретации методики И.Я. Штаермана [3] производили расчет упругих радиальных перемещений образующих поверхностей рабочих валков для начального и ко-
нечного 2 граничных сечений каждого отдельного і-го элементарного объема:
Ржі 2 =(хі 2 - Кхі1 )п / [К1 (1 -®2)
Схі2 = {СхіАі1 + 0,5 (Ржи + Ржі 2 )(2 - Кхі1 ) + ((і + Тхі2 )} < ^
2 2
8ІП (|СРг1 -Фіо|) АФ + К№Рхі1; (5)
х 1^8 1ф~2-ф«1 1 +(ЦМ1-3^1
2 2
^п(|ф,-2 -ф,-о|)ДФ+ К№Рхг2 , (6)
где ф,0 - усредненное значение угла контакта для данного г-го элементарного объема, по отношению к которому осуществляется численное интегрирование, ф,0 =(ф;1 + Ф,2 )/2; юв, Ев - коэффициент Пуассона и модуль упругости материала рабочих валков; Кк - коэффициент пропорциональности, являющийся упругой характеристикой контактных поверхностей рабочих валков
Учитывая упругие перемещения образующих поверхностей рабочих валков, количественные оценки 'Щ,2 которых получены на основе фор-
мул (5) и (6), принимая во внимание очевидные геометрические соотношения интегрального очага деформации (см. рис. 1), осуществляли перерасчет текущих значений толщин прокатываемых полос кхг1 и к, 2:
где 80102 - вводимая поправка первоначальной величины межцентрового расстояния 0102, определяемая исходя из условия кх ш1п = к1.
По мере определения новых значений кг1, кх,2 по формулам (7) и (8) производили расчет локальных характеристик напряженного состояния металла при данных текущих значениях толщин.
Рассматриваемая математическая модель включала в себя итерационный расчет силы, моментов и мощности прокатки. Помимо этого, был проведен расчет величины опережения, текущих и результирующих значений степени использования запаса пластичности, показателя напряженного состояния металла, а также величины напряжений переднего натяжения прокатываемой полосы и других параметров.
В качестве примеров результатов численной реализации полученной математической модели на рис. 2 представлены расчетные распределения локальных и интегральных характеристик напряженного состояния металла в зависимости от конечной толщины к1 и температуры прокатываемой полосы г.
[3].
К1 = к1 + 2К - ( ( - 1 ) 008 ф,1 - 8°1°2 ,
кхг 2 = к1 + 2К-( - 2 ) ф,- 2 -§°1°2 ,
(7)
(8)
1000
800
600
400
200
0
Рх V А
// и ч
К = 1,4 Я 1 1л' ч ч ч \ \ \ Ч\ V.
кх = 1,2 Гь г /г / г у/ ✓ к V/ >/'4 \ \ ^ \ \ \ч\ Л кх = 1,6
N 714 // V ч ч ч ч V'
Ы ъ—7 / н / V щ ч
Х12-Хц\1 =
-0,4 -0,2 0 0,2 0,4 0,6 О,
/ = Пг, I = щ
35
30
25
20
15
10
5
О
Р/В, кН/мм
М/В., кНм/мм
м/в
р/в^^ /г^ , = 900 °С
.к:: А... — 1 * *' *«.
У "А к-.. ^ »
*0 = 800 °С = 850 °С
0,7
0,6
0,5
0,4
0,3
0,2
ОД
О
1,2
1,3
1,4
1,5
1,6 мм
б
Рис. 2. Расчетные распределения локальных (а - ^ = 900 °С) и интегральных (б) характеристик напряженного состояния металла при горячей прокатке относительно тонких полос в зависимости от конечной толщины \ и температуры t0 (сталь 45,
к0 = 2,0 мм, Я = 350 мм, ^ = 0,35, У1 = 10 м/с). Р/В, М/В - приведенные к единице ширины проката сила и момент прокатки
а
Таким образом, распределение интегральных и локальных характеристик напряженного состояния носит довольно сложный характер, что подтверждает целесообразность использования численного подхода.
Список литературы
1. Мазур В. Л., Ноговицын A.B. Теория и технология тонколистовой прокатки (численный анализ и технические приложения). - Днепропетровск: PBA «Дніпро-VAL», 2010. - 500 с.
2. Caтонин A.B. Численная одномерная математическая модель процесса прокатки относительно тонких композиционных листов и полос, основанная на энергетическом подходе ІІ Удосконалення процесів та обладнання обробки тиском в металургії i машинобудуванш І ДДMA. - Краматорськ, 1998. -C. 36-41.
3. Федоринов B.A., Caтонин A.B., Грибков Э.П. Математическое моделирование напряжений, деформаций и основных показателей качества при прокатке относительно широких листов и полос І ДГMA. - Краматорск, 2010. - 243 с.
4. Целиков A.^, Никитин T.C., Рокотян C.E. Теория продольной прокатки. - М.: Металлургия, 1980. - 320 с.
5. Контактное трение в процессах обработки металлов давлением І A.ft Леванов [и др.]. - М.: Металлургия, 1976. - 416 с.
6. Caтонин A.B., Настоящая C.C., Переходченко B.A. Исследование напряженно-деформированного состояния металла при горячей прокатке относительно тонких полос ІІ Обработка материалов давлением: сб. науч. тр. І ДГМ^. - Краматорск, 2011. - № 3 (28). - C. 23-28.
Получено 1.11.2012
Сатонин Александр Владимирович - доктор технических наук, профессор, Донбасская государственная машиностроительная академия (84313, Украина, Донецкая обл., г. Краматорск, ул. Шкадинова, 72, e-mail: [email protected]).
Настоящая Светлана Сергеевна - аспирант, Донбасская государственная машиностроительная академия (84313, Украина, Донецкая обл., г. Краматорск, ул. Шкадинова,72, e-mail: [email protected]).
Нереходченко Виктор Александрович - аспирант, Донбасская государственная машиностроительная академия (84313, Украина, Донецкая обл., г. Краматорск, ул. Шкадинова, 72, e-mail: [email protected]).
Присяжный Андрей Григорьевич - старший преподаватель, Приазовский государственный технический университет (87500, Украина, Донецкая обл., г. Мариуполь, ул. Университетская, 7, e-mail: [email protected]).
Satonin Aleksandr Vladimirovich - Doctor of Technical Sciences, Professor, Donbass State Engineering Academy (84313, Ukraine, Donetsk region, Kramatorsk, Shkadinova st., 72, e-mail: [email protected]).
Nastoyashaya Svetlana Sergeevna - Graduate Student, Donbass State Engineering Academy (84313, Ukraine, Donetsk region, Kramatorsk, Shkadinova st., 72, e-mail: [email protected]).
Perehodchenko Viktor Aleksandrovich - Graduate Student, Donbass State Engineering Academy (84313, Ukraine, Donetsk region, Kramatorsk, Shkadinova st., 72, e-mail: [email protected]).
Prisyagniy Andrey Grigoryevich - Associate Professor, Azov State Technical University (87500, Ukraine, Donetsk region, Mariupol, Universitetskaya st., 7, e-mail: [email protected]).
