Mathematical modeling of the kinematics of metal forming in the screw calibers of a ball rolling mill

This article is available in Russian only.
Цитировать:
Математическое моделирование кинематики формообразования металла в винтовых калибрах шаропрокатного стана // Universum: технические науки : электрон. научн. журн. Каримов К.А. [и др.]. 2026. 9(150). URL: https://7universum.com/en/tech/archive/item/23440 (дата обращения: 17.10.2026).
Прочитать статью:
DOI - 10.32743/UniTech.2026.150.9.23440

УДК 621.771+519.6

Аннотация

В данной статье представлена математическая модель кинематики процесса формообразования металла при поперечно-винтовой прокатке. Модель основана на совместном описании геометрии винтовых калибров через функции переменного шага нарезки T(x) и радиуса валка R_{в}(x) . Разработаны системы дифференциальных уравнений первого порядка в матричной форме Коши для шести фазовых координат (перемещений, скоростей и ускорений заготовки по осевой, тангенциальной и радиальной составляющим) а также получены непрерывные функции распределения высоты, ширины, площади поперечного сечения и продольной скорости течения металла по всей длине очага деформации. Приведен алгоритм численного интегрирования системы методом Рунге-Кутты 4-го порядка в среде MathCAD с заданием начальных условий на входе заготовки в валки. Адекватность и математическая корректность модели верифицированы с помощью расчета закона постоянства секундного объема металла V_{сек} = const , где среднее отклонение не превышает 3–5 %.

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

Abstract

This article presents a mathematical model of the kinematics of the metal forming process during transverse helical rolling. The model is based on a combined description of the geometry of helical calibers using functions of the variable pitch T(x) and the roll radius R_{в}(x) . Systems of first-order differential equations in Cauchy matrix form have been developed for six phase coordinates (displacements, velocities, and accelerations of the workpiece along the axial, tangential, and radial components) and continuous distribution functions for the height, width, cross-sectional area, and longitudinal flow velocity of the rolled material have been obtained along the entire length of the deformation zone. An algorithm for the numerical integration of the system using the 4th-order Runge-Kutta method in MathCAD is presented, with initial conditions specified at the workpiece’s entry into the rolls. The adequacy and mathematical correctness of the model were verified by calculating the law of conservation of the instantaneous volume of metal, where the average deviation does not exceed 3–5 %.

The calculations revealed a parabolic nature of the vertical compression and a damped lateral expansion, as well as the determination of the critical section point, reflecting the transformation of the shape from circular to spherical. The calculated short-speed diagram illustrates a monotonic deceleration of the flow, determining a positive coefficient and confirming the stability of the kinematic rolling regime.

Ключевые слова: шаропрокатный стан; прокатные валки; винтовая линия; реборда; ось прокатки; математическая модель; кинематика; MathCAD 15.

Keywords: ball rolling mill; rolling rolls; screw line; flange; rolling axis; mathematical model; kinematics; MathCAD 15.

Введение

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

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

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

Материалы и методы

Поверхность винтового калибра в цилиндрической системе координат описывается переменным радиусом валка по вершине реборды и переменным шагом винтовой линии вдоль оси прокатки [5]. Длина очага деформации включает зоны захвата, интенсивного формообразования и калибровки (рис. 1).

Рисунок 1. Схема для определения проекций векторов скорстей заготовки
Рисунок 1. Схема для определения проекций векторов скорстей заготовки

При параллельной или скрещенной установке валков с межцентровым расстоянием A текущий радиус валка R_{в}(x) и текущий радиус шейки (перемычки) заготовки r_{п}(x) связаны соотношением:

R_{в}(x) = \frac{A}{2} - r_{п}(x) (1)

Для непрерывного дифференцирования геометрии в среде численного интегрирования профили валка R_{в}(x) , шейки заготовки r_{п}(x) и шага калибра T_{к}(x) задаются в виде кубических сплайнов или полиномиальных функций вида:

T_{в}(x) = T_{в0} - A_{t} \cdot x^{2} + B_{t} \cdot x^{3} (2)
r_{п}(x) = r_{п0} - \left( \frac{r_{п0} - r_{\min}}{l_{д}^{2}}()^{2} \right) (3)

где T_{к0} – начальный шаг винтовой линии в зоне захвата; r_{п0} – начальный радиус цилиндрической заготовки; r_{\min} – радиус перемычки в зоне отделения шара; A_{t} , B_{t} – коэффициенты апроксимации геометрии калибра.

Вектор полной скорости движения точек поверхности металлической заготовки в любой момент времени раскладывается на три основные компоненты [4]:

1. Осевая (поступательная) скорость \left( v_{x} \right) : задает скорость продвижения заготовки вдоль оси стана с учетом шага калибра T_{к}(x) и коэффициента осевого скольжения \eta_{x} :

v_{x} = \frac{\omega_{в}}{2\pi} \cdot T_{к}(x) \cdot \eta_{x} (4)

где \omega_{в} – угловая скорость вращения валков.

2. Тангенциальная (окружная) скорость v_{\tau} : определяет угловую скорость вращения заготовки \omega_{заг} с учетом тангенциального коэффициента скольжения \eta_{\tau} :

v_{\tau} = \omega_{заг} \cdot r_{заг}(x) = \omega_{в} \cdot R_{в}(x) \cdot \eta_{\tau} \Rightarrow \omega_{заг} = \omega_{в} \cdot \frac{R_{в}(x)}{r_{заг}(x)} \cdot \eta_{\tau} (5)

3. Радиальная скорость v_{r} : характеризует скорость обжатия металла внедряющимися ребордами валков:

v_{r} = \frac{dr_{заг}}{dt} = \frac{dr_{заг}}{dx} \cdot \frac{dx}{dt} = \frac{dr_{заг}}{dx} \cdot v_{x} (6)

Дифференцированием выражений скоростей по времени получены уравнения ускорений заготовки в очаге деформации:

Осевое ускорение:

a_{x} = \frac{dv_{x}}{dt} = \frac{\omega_{в} \cdot \eta_{x}}{2\pi} \cdot \frac{dT_{k}(x)}{dx} \cdot v_{x} (7)

Угловое ускорение (с учётом правила дифференцирования сложной функции):

\varepsilon_{заг} = \frac{d\omega_{заг}}{dt} = \omega_{в} \cdot \eta_{\tau} \cdot \left\lbrack \frac{1}{r_{заг}(x)} \cdot \frac{dR_{в}}{dx} \cdot v_{x} - \frac{R_{в}(x)}{r_{заг}(x)^{2}} \cdot v_{r} \right\rbrack (8)

Радиальное ускорение:

a_{r} = \frac{dv_{r}}{dt} = \frac{d^{2}r_{заг}(x)}{dx^{2}} \cdot v_{x}^{2} + \frac{dr_{заг}(x)}{dx} \cdot a_{x} (9)

Матричная форма системы дифференциальных уравнений (Форма Коши) Для численного решения кинематической задачи методом Рунге-Кутты 4-го порядка сформирован вектор фазовых координат, состоящий из осевого перемещения S_{x} , угла поворота \varphi_{заг} , текущего радиуса r_{заг} , а также их первых производных (скоростей) v_{x} , \omega_{заг} , v_{r} .

\left\{ \begin{matrix} \&\frac{d^{2}S_{x}}{dx^{2}} = \frac{\omega_{в} \cdot \eta_{x}}{2\pi} \cdot \frac{dT_{k}\left( S_{x} \right)}{dx} \cdot \frac{dS_{x}}{dt} \\ \&\frac{d^{2}\varphi_{заг}}{dt^{2}} = \omega_{в} \cdot \eta_{\tau} \cdot \left( \frac{1}{r_{заг}} \cdot \frac{dR_{в}\left( S_{x} \right)}{dx} \cdot \frac{dS_{x}}{dt} - \frac{R_{в}\left( S_{x} \right)}{{r_{заг}}^{2}} \cdot \frac{dr_{заг}}{dt} \right) \\ \&\frac{d^{2}r_{заг}}{dt^{2}} = \frac{d^{2}r_{заг}\left( S_{x} \right)}{dx^{2}} \cdot \left( \frac{dS_{x}}{dt} \right)^{2} + \frac{dr_{заг}\left( S_{x} \right)}{dx} \cdot \left( \frac{\omega_{в} \cdot \eta_{x}}{2\pi} \cdot \frac{dT_{k}\left( S_{x} \right)}{dx} \cdot \frac{dS_{x}}{dt} \right) \\ \end{matrix} \right.\ (10)

Результаты

Интегрирование системы дифференциальных уравнений (10) позволяет получить массивы кинематических параметров S_{x}(t) , \varphi_{заг}(t) , r_{заг}(t) , v_{x}(t) , \omega_{заг}(t) , v_{r}(t) , на основе которых построены основные расчетные графики формообразования.

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

Рисунок 2. Изменение радиуса перемычки заготовки по длине очага деформации S
Рисунок 2. Изменение радиуса перемычки заготовки по длине очага деформации S

Полученная кривая r_{п} = f\left( S_{x} \right) иллюстрирует непрерывный характер радиусной деформации заготовки при её движении сквозь очаг деформации длиной l_{д} = 200мм .

Кинематика вращательного движения металла и интенсивность скручивания заготовки в процессе винтовой перекатки исследованы на основе зависимости угловой скорости заготовки \omega_{заг} от её осевого перемещения S_{x} (рис. 3).

Рисунок 3. Зависимость угловой скорости вращения заготовки от осевого перемещения
Рисунок 3. Зависимость угловой скорости вращения заготовки от осевого перемещения

Анализ зависимости \omega_{заг} = f\left( S_{x} \right) показывает нелинейный гиперболический характер роста угловой скорости заготовки по мере её продвижения вглубь очага деформации.

Закономерности поступательного перемещения металла вдоль оси прокатки и динамика его вытяжки оценены по полученному графику изменения осевой скорости заготовки v_{x} от координаты очага деформации S_{x} (рис. 4).

Рисунок 4. Изменение осевой скорости движения заготовки вдоль оси прокатки S
Рисунок 4. Изменение осевой скорости движения заготовки вдоль оси прокатки S

Характер кривой v_{x} = f\left( S_{x} \right) отражает монотонное замедление осевого продвижения заготовки по длине бочки валка.

Снижение осевой скорости задается законом профилирования винтовой линии валка T_{k}(x) , шаг которой сужается к зоне калибровки до диаметра готового изделия 2R_{ш} . Такое замедление создает в очаге деформации эффект осевого подпора: идущие сзади слои металла с большей скоростью протяжки уплотняют металл в зоне интенсивного обжатия.

Выводы

1. Разработана аналитическая модель геометрии и кинематики, описывающая вращательно-поступательное движение заготовки в очаге деформации винтового калибра. Модель учитывает непрерывное изменение профиля шейки шара r_{заг}(x) , шага винтовой линии T_{k}(x) и коэффициентов скольжения металла \eta_{x} , \eta_{\tau} .

2. Сформирована система дифференциальных уравнений в матричной форме, пригодная для реализации в современных пакетах компьютерного моделирования (MathCAD). Представление уравнений в форме Коши позволяет рассчитывать полные кинематические поля (траектории, осевые v_{x} , угловые \omega_{заг} и радиальные v_{r} скорости, а также соответствующие ускорения заготовки) на любом этапе формообразования.

Список литературы

1. Аносов А.С., Голубчик Э.М. Прокатка помольных шаров в винтовых калибрах: теория и практика. С. 95⁠–⁠103. – Магнитогорск: Изд-во МГТУ им. Г.И. Носова, 2012. – 185 с.

2. Выдрин В.Н., Черных И.Н. Динамика и кинематика процессов поперечно-винтовой прокатки. С. 70⁠–⁠120. – Челябинск: ЮУрГУ, 2005. – 160 с.

3. Котенок В.И., Подобедов С.И. Создание эффективных калибровок валков шаропрокатных станов и расширение сортамента шаров на действующих и новых станках // Труды III конгресса прокатчиков, Липецк, 19⁠–⁠22 окт. – 2000. – С. 438⁠–⁠441.

4. Романенко В.П., Степанов А.И. Моделирование формоизменения металла в очаге деформации шаропрокатного стана // Обработка материалов давлением. – 2015. – № 2.

5. Филатов А.А. Кинематика течения металла при винтовой прокатке заготовок постоянного и переменного сечения // Известия вузов. Черная металлургия. – 2018. – Т. 61, № 4.

6. Целиков А.И., Барбар В.И., Василев Я.Д. Специальные станы для поперечно-винтовой прокатки. С. 185⁠–⁠210. – М.: Металлургия, 1971. – 248 с.

7. Целиков А.И., Барбарич М.В., Васильчиков М.В., Грановский С.П., Жукевич-Стоша Е.А. Специальные прокатные станы // М.: Металлургия. – 1971.

8. Granovski S.P. New Technologies and Mills for Rolling of Products in Screw Calibers. Moscow: Metallurgiya. – P. 62⁠–⁠87. – 1980.

9. Shakhobutdinov R.E., Karimova A.R. & Nosirov T.N. Tasks of mastering the production of import-substituting steel grinding balls products. Science and Innovative Development, – P. 36⁠–⁠44. – 2024.

10. Tofil, A. Overview of the research on roll forging process / A. Tofil, Z. Pater // Advances in Science and Technology Research Journal. – 2017. – Т. 11, № 2. – P. 72⁠–⁠86.

11. Tofil, A. Overview of the research on roll forging process / A. Tofil, Z. Pater // Advances in Science and Technology Research Journal. – 2017. – Т. 11, № 2. – P. 72⁠–⁠86.

References

1. Anosov A.S., Golubchik E.M. [Rolling of grinding balls in screw calibers: theory and practice]. Magnitogorsk, Publishing House of Nosov Magnitogorsk State Technical University, 2012. 185 P. (In Russ.)

2. Vydrin V.N., Chernykh I.N. [Dynamics and kinematics of cross-screw rolling processes]. Chelyabinsk, South Ural State University, 2005. 160 P. (In Russ.)

3. Kotenok V.I., Podobedov S.I. [Creation of effective roll calibrations for ball rolling mills and expansion of ball assortment on existing and new machines]. Trudy III kongressa prokatchikov, Lipetsk, October 19⁠–⁠22, 1999. Moscow, JSC «Chermetinformatsiya», 2000, P. 438⁠–⁠441. (In Russ.)

4. Romanenko V.P., Stepanov A.I. [Modeling of metal deformation in the deformation zone of a ball rolling mill]. Obrabotka materialov davleniem, 2015, no. 2, P. 112⁠–⁠119. (In Russ.)

5. Filatov A.A. [Kinematics of metal flow during screw rolling of billets of constant and variable cross-section]. Izvestiya vuzov. Chernaya metallurgiya, 2018, vol. 61, no. 4, P. 295⁠–⁠302. (In Russ.)

6. Tselikov A.I., Barbar V.I., Vasilev Ya.D. [Special mills for cross-screw rolling]. Moscow, Metallurgiya Publ., 1971. 248 P. (In Russ.)

7. Tselikov A.I., Barbarich M.V., Vasilchikov M.V., Granovsky S.P., Zhukevich-Stosha E.A. [Special rolling mills]. Moscow, Metallurgiya Publ., 1971. 336 P. (In Russ.)

8. Granovski S.P. New Technologies and Mills for Rolling of Products in Screw Calibers. Moscow, Metallurgiya. – 1980. – 87 P.

9. Import-substituting steel grinding balls products. Science and Innovative Development, 2024, P. 36⁠–⁠44. (In Russ.)

10. Tofil A., Pater Z. Overview of the research on roll forging process // Advances in Science and Technology Research Journal, June. – 2017. – P. 72⁠–⁠86.

11. Tofil A., Pater Z. Overview of the research on roll forging process // Advances in Science and Technology Research Journal, June. – 2017. – P. 72⁠–⁠86.

Информация об авторах
professor, Tashkent State Technical University named after Islam Karimov, Uzbekistan, Tashkent
Associate Professor, Tashkent State Technical University named after Islam Karimov, Uzbekistan, Tashkent
Associate Professor, Tashkent State Technical University named after Islam Karimov, Uzbekistan, Tashkent
senior Lecturer, Tashkent State Technical University named after Islam Karimov, Uzbekistan, Tashkent
senior Lecturer, Tashkent State Technical University named after Islam Karimov, Uzbekistan, Tashkent
ISSN 2311-5122. Article metadata is hosted on the eLIBRARY.RU platform.
Mass media registration cert.: EL No. FS77-91806 dated 17.06.2026
Journal founder: Universum LLC
Editor-in-Chief - Marina Yu. Zvezdina.
Top