УДК 621.744.35 DOI: 10.14529/теМ70210
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ВИБРАЦИОННОГО ПРЕССОВАНИЯ ТРУБ
С.Р. Рахманов
Национальная металлургическая академия Украины, г. Днепропетровск, Украина
Разработана математическая модель динамики процесса прессования осесимметричной трубной заготовки на профильном прессе трубопрессовой линии. Согласно расчетной схеме и принятой динамической модели процесса прессования трубной заготовки получены аналитические выражения для соответствующих взаимодействующих границ пластической и упругой зон очага деформации. Получена аналитическая зависимость, характеризующая скорость распространения фронта пластической деформации в зависимости от скорости перемещения прессуемого металла. Определены динамические характеристики возникновения волновых явлений в прессуемом металле с учетом условий воздействия со стороны инструмента, что существенно дополняет ранее известные представления о напряженно-деформированном состоянии металла в очаге деформации. Установлены поля деформаций в зонах возмущенной области очага деформации с учетом воздействия технологических инструментов. Путем моделирования динамических процессов в прессуемом металле установлены особенности реализации высокоэффективной технологии управляемого вибрационного прессования труб на профильном прессе трубопрессовой линии. Выбором рациональных параметров (амплитудно-частотной характеристики) используемого активного процесса вибрационного прессования трубной заготовки достигнуто значительное повышение эффективности процесса и качества прошиваемых заготовок.
Ключевые слова: труба; пресс; матрица; металл; игла; калибр; очаг деформации; напряжение; деформация; модуль упругости; модуль сдвига; пластическая область; условия пластичности; упругая область; вибрация; частота; амплитуда; период.
Введение
Повышение эффективности процессов прессования бесшовных труб в условиях применения высокочастотных активных вибраций создает определенные условия для расширения технологических возможностей эксплуатируемого прессового оборудования. Разработка новых, высокоэффективных технологий производства бесшовных труб выдвигает, при этом, ряд актуальных задач. В мировой практике весьма перспективным направлением интенсификации режимов деформации при прессовании бесшовных труб считается использование активных высокочастотных ультразвуковых колебаний или механических вибраций. Применение передовых технологий при производстве бесшовных труб требует кардинально нового подхода к задаче вибрационного прессования, что, в свою очередь, связано с отработкой режимов деформации, прогнозированием показателей стойкости технологического инструмента (матриц, игл, контейнеров), оптимизацией качества бесшовных труб и повышением на-
дежности функционирования всего прессового оборудования [1-3].
Цель работы. Разработка математической модели и совершенствование современных вибрационных технологий производства горячекатаных труб. Постановка задачи. Современные технологические процессы прессования труб из малопластичных и трудноде-формируемых металлов на профильных прессах трубопрессовых установок с применением высокочастотных вибраций выдвигают ряд требований по совершенствованию как технологии производства труб, так и модернизации технологического оборудования [4, 5]. Последнее непосредственно связано с отработкой рациональных режимов деформации и прогнозированием показателей стойкости технологического инструмента, надежности функционирования оборудования для реализации самих высокочастотных вибраций [5] (рис. 1).
В ряде случаев с целью повышения эффективности технологические процессы прессования бесшовных труб предусматривают
Рис. 1. Схема расположения базового технологического инструмента на трубопрофильном прессе двойного действия: 1 - пресс-штемпель со встроенным вибратором; 2 - матрицадержатель со встроенным вибратором; 3 - пресс-шайба; 4 - игла (оправка) со встроенным вибратором; 5 - иглодержатель; 6 - контейнер; 7 - внутренная втулка; 8 - центрирующая втулка иглодержателя со встроенным вибратором; 9, 10 -опорные кольца пресса; 11 - опорное кольцо матрицы; 12 - промежуточная втулка; 13 - опорное кольцо матрицадержателя; 14 - матрица, 15 - трубная заготовка; 16 - экспандер или оправка
использование управляемых высокочастотных вибраций. Однако в силу недостаточной изученности, эти процессы не находят широкого применения в практике производства бесшовных труб. В частности, задачи по повышению эффективности процессов прессования труб на трубопрофильном прессе, в условиях применения высокочастотных вибраций, решаются без должного аналитического исследования и, как правило, экспериментально.
Выполнение операций прессования на профильном прессе преимущественно сопровождается значительными потерями рабочего усилия на преодоление сил сопротивления, характерного для процесса неравномерности деформаций по всему объему обрабатываемого металла.
Известно, что наиболее благоприятные условия деформации прессования труб достигаются в условиях использования высокочастотных вибраций, что получило некоторое свое отражение в ряде фундаментальных исследований, как в теоретическом плане, так и экспериментального характера.
Анализ литературных источников
Для обоснования и принятия решений по ведению прессования труб в активном управляемом режиме подробно проанализируем некоторые важные моменты известных процессов пластической деформации металла, моделируя возникающие волновые процессы в условиях применения высокочастотных активных вибраций.
Остановимся на некоторых моментах технологического процесса вибрационного
прессования труб. Процесс прессования труб на профильном прессе производится следующим образом: предварительно нагретая заготовка для прессования поступает в контейнер, имеющий температуру 380-420 °С. Вибри-руемая игла (оправка) принудительно, посредством стержня внедряется в предварительно сверленную или экспандированную трубную заготовку. При дальнейшем движении пресс-штемпеля и вибрируемой оправки металл заготовки уплотняется и занимает ограниченный объем пространства втулки контейнера вдоль оси прессования. Далее прессуемый металл трубной заготовки посредством вибрируемого пресс-штемпеля вдавливается в кольцевой зазор, образованный активно вибрируемой матрицей и иглой. Следует принимать во внимание то, что при этом в процессе вибрационного прессования в поперечном направлении очага деформации возникают волновые процессы, которые в дальнейшем сопровождаются динамическими эффектами, действующими на металл заготовки через технологические инструменты.
Установление некоторых скрытых возможностей существующих технологий и разработка новой технологии вибрационного прессования труб требует создания корректных математических моделей, адекватно отражающих условия реализации и особенности данного процесса на современных профильных прессах. В основу разработки ряда математических моделей и создания механизмов для вибрационного нагружения металла положены некоторые разработки В.М. Клименко, В Н. Потураева, Г.Б. Бухвалова, В.П. Сверден-
ко, В.И. Залесского, В.Н. Шаповала и других авторов; анализ ряда разработок обобщен и наиболее полно представлен в работах [6-11].
Большинством исследователей подтверждается тот факт, что применение высокочастотных вибраций при обработке металлов давлением приводит к уменьшению предела текучести материала и, как следствие, к значительному снижению энергосиловых параметров очага деформации. При этом важную роль играет фактор формирования рациональной объемной деформации металла и уменьшения внешних контактных и внутренних сил трения. Результатами исследований, приведенных в работах [12-16], установлена некоторая зависимость усилия прессования от частоты колебаний активных вибраций. При этом выявлено, что для достижения установившегося вибрационного режима деформирования заготовки необходимо определение зависимости между частотой, амплитудой и кинематическими параметрами технологического процесса прессования. Следовательно, амплитудно-частотная характеристика базовых технологических инструментов (пресс-штемпеля, матрицы, иглы) оказывает существенное влияние на эффективность процесса, поскольку при этом формируются некоторые условия для активизации вибрационного деформирования заготовки в процессе прессования.
Первые опытно-промышленные установки с использованием ультразвуковых колебаний были построены фирмами Uniform Tubes и Aero projects (США) в начале 60-х годов ХХ столетия. Затем на заводе фирмы Aero projects для производства бесшовных труб был создан пресс с подводом высокочастотных ультразвуковых колебаний непосредственно к матрице. Применение ультразвуковых колебаний позволило увеличить степень единичного обжатия и повысить скорость деформации заготовки, уменьшить число операций отжига и травления, а также количество отходов. Данные вопросы свое дальнейшее решение нашли на заводе фирмы Stainless Tube Division (США), где производство бесшовных труб из коррозионно-стойкой стали на профильных прессах осуществляют высокочастотными ультразвуковыми колебаниями технологического инструмента (матрицы) [17].
Промышленная технология прессования бесшовных труб из предварительно просверленных заготовок с использованием энергии
ультразвуковых колебаний освоена также в Англии и Японии. Фирмами Levy и Kolambia Jammeril (Англия) были созданы прессы усилием прессования до 40 МН, которые эксплуатировались в технологических линиях трубопрессовых установок для производства труб диаметром до 380 мм. Процессы прессования реализовались путем подвода к пресс-штемпелю ультразвуковых колебаний частотой в пределах 5-15 кГц [18, 19].
Аналогичный профильный пресс с подводом ультразвуковых колебаний непосредственно к игле, через стержень, работает в фирме «Кобэ» (Япония) [20, 21].
На трубопрокатных заводах России и Украины в настоящее время работают трубопрокатные линии, оснащенные современными прессами, выпуска фирмы SMS MEER (Германия) или отечественного производства. Однако по ряду необоснованных причин данное оборудование не снабжено устройством для реализации современной технологии вибрационного прессования труб.
Метод решения задачи
Для выработки научно обоснованных предложений рассмотрим несколько иную постановку и решение задачи о прессовании предварительно нагретой гильзы вибрируе-мой матрицей и иглой. Игла и матрица, используемые в качестве технологических инструментов профильного пресса, представляют осесимметричные абсолютно твердые тела. В ходе реализации высокочастотного вибрационного прессования труб одновременно вибрируемая игла и матрица воздействуют на прессуемый металл (однородную среду), занимающий некоторое полупространство нормально к поверхностям технологических инструментов. При этом возмущенное движение прессуемого металла в очаге деформации, вызванное вибрируемой матрицей и иглой, в процессе прессования, в радиальном направлении, предполагаем строго осесимметрич-ным. Выбираем следующую расчетную схему процесса прессования труб: используем цилиндрическую систему координат с началом в точке соприкосновения оправки с прошиваемой заготовкой; время прошивки отсчитываем от момента соприкосновения матрицы и иглы с металлом, причем ось x направляем вдоль оси прессования заготовки (рис. 2).
Решение задачи радиального волнового движения металла в процессе прессования,
л х = const
Рис. 2. Расчетная схема динамических процессов прессования и гильза с упруго-пластической зоной: 1 - пресс-штемпель; 2 - гильза; 3 - контейнер; 4 - игла (оправка); 5 - матрица; 6 - труба
основанное на использовании строгих нелинейных дифференциальных уравнений, в силу больших деформаций в настоящее время из-за недостаточного развития математического аппарата является затруднительным. Поэтому на данной стадии исследования наиболее приемлемы те методы, которые позволяют получить решение задачи в аналитической форме. Очевидно также, что данные методы хотя и обладают определенной погрешностью, но в основном охватывают физическую суть явлений, происходящих в процессе прессования предварительно просверленных или экспан-дированных заготовок в гильзы (трубы) заданных размеров.
Пусть H(t) - есть глубина «внедрения» матрицы в металл в момент времени t. Если уравнение образующей рабочей поверхности матрицы задано в виде r = f (x), (x > 0), то положение матрицы относительно выбранной системы координат в момент времени t > 0 согласно методу инверсии определяется r = f [ x + H (t)].
В некоторой фиксированной плоскости x = const сечение движущейся матрицы в каждый момент времени представляет собой сужающуюся окружность (кольцо) переменного радиуса r0(t) = f [t + H(t)], изменяющуюся в радиальном направлении со скоростью
Co(t) = ro(t) = df (x + H (t)) = dt
= f ' (x + H (t)K(t), (1)
где u0 (t) = H(t) - скорость прессования металла; f' (x + H (t)) - тангенс угла наклона
образующей рабочей поверхности матрицы к оси x на некоторой заданной глубине прессования.
Если образующая рабочей поверхности матрицы выполнена в виде эквивалентного конуса с углом полураствора , то f (x) = x tg(9^), то скорость сужения поперечного сечения матрицы в плоскости x = const будет равна c0(t) = u0(t)tg(9^).
Принимаем во внимание то, что технологический процесс прессования труб на профильном прессе ведется в квазиустановив-шемся режиме с постоянной скоростью. Поэтому можно в первом приближении считать, что U0 (t) = U0 = const. При этом каждое поперечное сечение матрицы при ее перемещении сужается с постоянной скоростью c0 = u0 tg(9^) = const, что дает основание рассматривать движение металла в радиальном направлении как автомодельное.
Предложенный подход к решению задачи вибрационного прессования состоит в следующем: считается, что каждый слой металла в некоторой фиксированной плоскости x = const движется зависимо от других слоев, причем движение начинается с того момента времени, когда рабочая поверхность матрицы достигает данной плоскости. Исходная задача сводится к задаче о сужении цилиндрической «каверны». Следовательно, данная задача представляет собой задачу о распространении одномерных радиальных волн в слое, перпендикулярном направлению движения металла в очаге деформации. Таким образом, можно сказать, что предлагаемый подход к решению
поставленной задачи основывается исключительно на методе плоских сечений, что позволяет развивать известные положения процессов деформирования заготовки технологическими инструментами (матрицей и иглой).
Результаты ряда экспериментальных исследований [2, 6] указывают на необходимость поэтапного исследования деформированного состояния прессуемого металла в стесненном очаге деформации между матрицей и иглой. Напряженно-деформированное состояние металла в очаге деформации показывает, что при существующих скоростях прессования в цилиндрической заготовке (металле) кроме зоны упругого состояния, преимущественно возникает зона пластического состояния, а в некоторых случаях и зона упругих деформаций. Поэтому первоначально рассмотрим отдельно указанные зоны и дадим решение динамической задачи для каждой из них. Пластическая область очага деформации расположена в зоне непосредственного контакта металла с матрицей. В рамках рассматриваемой задачи и выбранной динамической модели процесса в каждом сечении x = const данная зона представляет собой концентрическую область, ограниченную окружностью радиусом (поперечное сечение матрицы) r0 (t), которая пропорционально калибровке матрицы сужается со скоростью Со (t), и окружностью некоторого радиуса R(t), которая определяет границу пластической области со стороны упругой зоны и изменяется с неизвестной заранее скоростью с = r (t).
Область, примыкающая к пластической зоне по окружности r = R(t), как будет показано в дальнейшем, может рассматриваться как область, подверженная, в некоторых случаях, определенному виду разрушения.
В выбранной цилиндрической системе координат законы сохранения массы и количества движения для радиального перемещения металла в рассматриваемом слое очага деформации x = const имеют вид
др др (до о
— + о— = -р| — + - |; dt дг I дг r
да
Г | аг— а9
до до
= -р| — + о — дt дг
(2)
(3)
Зг г
где р(г, t) - плотность металла; о(г, t) - радиальная скорость частиц металла (положительная в сторону уменьшения размеров мат-
рицы); аг, а9 - соответственно радиальная и окружная компоненты напряжения (положительные при сжатии).
Уравнения (2) и (3) применимы как для пластической зоны, так и для упругой зоны. Пластичность металла, при рассматриваемом подходе, наиболее удобно описывать критерием текучести по Треску, который в нашем случае имеет вид
ar — а9 = ^0 , (4)
где т0 = const - предел текучести прессуемого металла.
Кроме того необходимо отметить, что относительное изменение объема, как известно из [23, 24], является упругой деформацией, которая обычно в теории пластичности считается пропорциональной давлению p = ~(cr +с9 +сх). В этом случае предполагаем, что между давлением p и изменением объема очага деформации существует линейная зависимость, которую представим в виде
p = K
f Р ^ 1 — ро
Р
(5)
где ро, р - соответственно начальная и текущая плотность прессуемого металла; К -объемный модуль упругости прессуемого материала.
Система четырех уравнений (2)-(5) не является замкнутой, так как содержит пять неизвестных функций: о, р, сг, се, сх.
Классическая теория пластичности при рассмотрении деформации металла с цилиндрической симметрией использует различные гипотезы и предположения относительно взаимосвязи напряжения сх с другими компонентами полного напряжения и деформаций. Одним из таких предположений в теории пластичности является предположение о линейности связи между девиатором тензора напряжений и девиатором тензора деформаций. Совместно с условием деформации указанное предположение приводит к следующей связи между напряжениями
а х = 2-К +а9) — G =
(6)
где г - объемная деформация материала.
Поскольку г величина достаточно малая, то в уравнении (6) часто пренебрегают по-
следним членом, в результате чего получаем
сх = ~(°г +се), что рекомендуется Хиллом
[23]. Однако результаты исследования [24] говорят о том, что вместо условия (14) в ряде случаев можно использовать условие сх = Се .
Необходимо указать, что оба последних условия принадлежат к классу более общего параметрического условия сх = 5Се + (1 - s)сr, соответственно при значениях параметра s = 1/2 и s = 1. Если использовать условие (6), то для
давления металла р = 1 [(2 - s)сr + (1 + я)се].
В силу изложенного, из (5), с учетом условия текучести, имеем
1-Р0. = 1
р
K
cr -3(1 + s)x0
(7)
Исключив из уравнений (2) и (9) Се при помощи условия (6) и р при помощи равенства (7), получим
^ + IKO + ^^.. I =
, 8о 8о .
:"р0 K|¥ + V I ;
(8)
8r r
(8c, 8c r 8r
1 + s
. —
8t
(9)
Для матрицы с произвольной криволинейной образующей г = / (х) и скоростью
движения пресс-штемпеля и0 (7), именно эти
два уравнения определяют пластическое движение металла в возмущенной инструментами области очага деформации г0 (7) < г < R(t).
Предположим, что пресс-штемпель трубную заготовку с постоянной скоростью «0 вдавливает в кольцевой зазор между матрицей и иглой. В этом случае граница контакта пластической зоны с матрицей в радиальном направлении движется с постоянной скоростью с0 = . В этом случае скорость распро-
странения пластической зоны с также постоянна и равна R (7) = ct. Следовательно, при
этом поставленная задача (8), (9) оказывается автомодельной, поскольку не имеет характерного размера. Перейдя в уравнениях (8) и (9) к
безразмерным величинам — = —; сг =—;
с К
т0 = , которые являются функциями только
г
автомодельной переменной в виде ^ = —,
М
получаем систему двух обыкновенных дифференциальных уравнений:
1
1 -Or + 3 (1 + s Ь
2 d о
= р2 dL ^ ;
1 -Cr + 1 (1 + s )Xq
fd Cr + О
dL L
( AT, гЛ
(10)
d о о
dL + L.
d c.
(L-о),
(11)
R С K
где P = —, cp =—. Cp P0
Система уравнений (10) и (11) описывает пластическое течение в области очага деформации c0)t < r < ct, в автомодельных переменных является интервалом L0 — L — 1. Здесь L0 - нижний предел автомодельной перемен" у С0 1
ной Lo = — < 1 соответствует границе контак-c
та металла с иглой.
Обратим внимание на тот факт, что граница r = R(t) заранее неизвестна, и в автомодельной задаче (из-за выбора в качестве характерной скорости величины c) неизвестным оказывается значение L0 на границе с матрицей. Таким образом, в подвижной области r0 (t) < r < R(t) с заранее неизвестной границей R(t) для матрицы с произвольной (криволинейной) образующей, решение задачи в пластической зоне сводится к решению системы нелинейных уравнений в частных производных (8) и (9).
Для случая прессования труб на матрице, выполненной в виде строго конического тела, с постоянной скоростью движения U0, задача в конечном итоге сводится к решению системы обыкновенных дифференциальных уравнений (10) и (11) в области с фиксированными границами L = L0 < 1 и L = 1 с заранее неизвестной величиной L0 . Из граничных условий при L = L0 известно значение скорости о = С0
0I
3
или u
= — = £0 , условия на границе £ = 1
' c
должны определяться из условия стыковки решения в пластической зоне с решением в зоне, которая примыкает к пластической. Отметим, что линеаризованная задача для пластической области базируется на основных предположениях линейной теории пластичности [23].
Результаты анализа задачи показывают, что хотя небольшое изменение плотности прессуемого металла происходит, но величина PolР ~ 1. Кроме того, величины и— и
дг
du
и— имеют незначимый второй порядок, по-
дг
этому, в линейном приближении задачи, вместо уравнений (8) и (9) воспользуемся уравнениями:
др
dt '
öOr
dr
■Р0 ^ du dr +u )=
r ■ = - du "р0 ¥ •
(12)
(13)
Исключив из уравнений (12) и (13) р, при помощи уравнения (7) придем к следующей системе уравнений [24]:
dcr т0 du
—-+— = -р0—•
dr r dt
dö + K idU + U| = o,
dt
dr r
(14)
(15)
которая является линеаризованным вариантом системы (8) и (9).
Отметим, что для матрицы с конической образующей, при постоянной скорости прессования м0, задача (30) и (31) является автомодельной и в безразмерных переменных приводится к системе обыкновенных дифференциальных уравнений:
dör + ^0 = ß2c • d £ £ d £'
du u „ dcr
— + - = £—L .
d£ £ d£
(16) (17)
на интервале переменной — £ — 1.
В дальнейшем ограничимся только решением линеаризованной системы уравнений (16) и (17), которая, как показывает сравнение с численным решением нелинейной системы (8) и (9), имеет удовлетворительную точность,
но обладает тем преимуществом, что линеаризованная система допускает замкнутое аналитическое решение. Из (16) и (17) имеем самостоятельное уравнение для определения й
(1 -Р2£2
d u и _
+Г-Т'
0'
(18)
Причем при £ = £о должно удовлетворяться граничное условие й|що = £0 . Решение уравнения (18), удовлетворяющее последнему условию, имеет вид
41 -Р2£2 =
>(£) = D
Р2£
Р2£
(19)
где D0 =
Р2£2 +^0 ф-Р2£2
Теперь из (16) получаем выражение для определения ör (£)
ö (£) = D>in
1 + 41 -Р2£2
Р£
+ D
(20)
где Dl - произвольная постоянная.
Постоянная Dl и безразмерная скорость
с
распространения пластической зоны р = —
должны быть определены из условия сопряжения пластической зоны с упругой зоной, примыкающей к пластической. Так как для внешней зоны могут применяться различные математические и физические модели, то здесь используется наиболее распространенная модель: модель упругой зоны в трубной заготовке.
Следует отметить, что при реализации некоторых режимов прессования в предварительно прошитой и нагретой трубной заготовке, наряду с пластической зоной, возникает некоторая упругая область (рис. 3). В связи с этим исследование динамических явлений в упругой области очага деформации является актуальным.
Упругая область очага деформации (в силу особенностей реализации технологического процесса прошивки) ограничена некоторой «жесткой» областью технологического инструмента (границами контейнера). В упругой области очага деформации имеем известные уравнения теории упругости [25], которые для радиального движения металла имеют вид:
а)
б)
в)
г)
Рис. 3. Осциллограммы процесса прессования экспандированой трубной заготовки на трубопрофильном прессе 50 МН фирмы SMS MEER (заготовка 0180 мм, материал 12Х18Н10Т, готовая труба 0105x10): а - ход траверсы пресса; б - сила прессования трубной заготовки; в - скорость прессования трубной заготовки; г - температура трубной заготовки перед прессованием
8а
8r
8 2u
r
8t
2
+ U V 2G 8U ;
8r r) 8r
i 8u u ^ u -а = а| —+ — | + 2G—; 8r r
8u u
-а = А | -+ —
8r r
(21) (22)
(23)
(24)
где u(r, t) - радиальное перемещение частиц
А vE металла; А=-
, ,-; Е,G - соответствен-
(1+у)(1-2у)
но модуль упругости Юнга и модуль сдвига прошиваемого металла; V - коэффициент Пуассона.
Напряжение по-прежнему считается положительным для сжатия, и в выражении (21) учтено, что в рамках рассматриваемого подхода осевое удлинение прошиваемой заготовки сравнительно мало.
Исключив из уравнения (21) напряжения сг и се согласно (22) и (23), после преобразований получим уравнение волновых процессов для возмущенной области очага деформации в виде
82u
1 8u + —■
8r2 r 8r r
a2 8t2 :
(25)
А + 2G
a =
P
Из уравнения (25) следует, что скорость распространения упругой цилиндрической волны в прошиваемой трубной заготовке равна скорости распространения продольных возмущений a . Момент времени, когда прессуемый металл касается матрицы и иглы, рассматриваемое сечение x = const следует считать нулевым. Очевидно, что к моменту времени t упругая волна в данном слое очага деформации захватит область радиусом r = at. Отметим, что если упругая зона очага деформации по исходному предположению задачи соседствует с пластической, то она расположена в кольце R(t) < r < at.
В случае, если калиброванный участок матрицы является коническим и поток металла движется с постоянной скоростью, то упругая зона расположена в кольце ct < r < at. При этом замкнутое решение задачи можно реализовать путем конструирования возмущенной области очага деформации.
Перейдем к детальному анализу математической модели возмущенной области очага деформации (25) и изучению волновых процессов в прессуемом металле.
С учетом определенных начальных и граничных условий задачи
«м> = f(r); drU = F(r); dr\r==R =0 (26)
dr dr
решение уравнения (25) представляем по методу Фурье в виде
u (r, t) = T (t)W (r). (27)
После подстановок (27) в дифференциальное уравнение волновых процессов (25) имеем
W "(r )T (t) +1W' (r )T (t) - -2 W (r )T (t) =
= ^ T"" (t )W (r).
(28)
Разделяя переменные задачи после некоторых преобразований (28) запишем
W (r) + 1w (r) - -2 W (r) T" (t)
W (r)
a2T (t)
= -X2. (29)
Из уравнения (29) имеем следующую систему независимых дифференциальных уравнений:
2 1
Ж " (г) + - Ж '(г) - (Я2 —-)Ж (г) = 0; (30)
1
— I
r
2„2п
Решения дифференциальных уравнений (30) и (31) согласно [26] соответственно запишем в виде
Ж(г) = СХ3Х (Хг) + С2Г1 (Яг); (32)
Т(0 = А cos(Хat) + В sin(Хat) (33)
где С1, С2, А, В - постоянные интегрирования задачи; ^ (Хг), У1 (Хг) - функции Бесселя соответственно первого и второго рода.
Дифференциальное уравнение (30) представляет известное уравнение Бесселя. Отметим, что в приведенном решении (32) при г ^ да второе слагаемое обращается в бесконечность 11(Хг) ^ да . Следовательно, при этом для удовлетворения начальных условий задачи постоянные интегрирования задачи соответственно определены в виде С = 1; С2 = 0 .
Тогда решение уравнения Бесселя в данном случае представим Ж(г) = ^(Хг).
Чтобы решение задачи удовлетворяло выбранному граничному условию справедливо —— |г=Я = 0 . Необходимо подчеркнуть, что dг
при этом из последнего условия находим требуемые корни (собственные числа) уравнения Бесселя
—( МХг) )| г=д = 0
или
Ji(XR) + J0 (XR) = 0.
XR
(34)
Из решения уравнения (34) с учетом некоторой замены ХЯ = ц находим все собственные числа (корни) уравнения Бесселя ц^, к = 1, 2, 3,...
Отметим, что собственным числам Хк = Цк / Я уравнения Бесселя соответствует собственная функция
W(r) = |.
(35)
Согласно методу Фурье решение уравнения (25) с учетом (32) и (35) представим в виде
<(r, t)=£
T"(t) + X2a T(t) = 0.
(31)
xJ
k=1
Mkat R
Ak cos
V-kat R
+ Bk sin
V-kat R
(36)
r
r
x
r
Из принятых начальных условий задачи запишем
f(r) = 1Л-Л [^j и F(r) = £BJi [^J .
Соответственно определяем постоянные интегрирования задачи Лк и Bk в виде:
Mtr
Лк =
Bk =
2
R2МЩ о 1 R
2
RaJ^T) о R
f rf (r) Ji( ^ )dr;
J rF (r) Ji( M- )dr.
(37)
Решение задачи при этом удобно представить в стандартной форме гармонических колебательных процессов
Ы(т, г) = ¿А яп( ^ + Фк )/1 ], (38)
k=1
R
R
где
Лк =у1Л,
k + Bk2 - амплитуды волновых
процессов в прошиваемом металле; ф^ =
Bt
начальные фазы волновых процессов.
Период волновых процессов в прошиваемом металле определяется
Т = ^
T = Kk '
где Kk =
R
(39)
частота распространения
упругих волн в прессуемом металле.
Для рассматриваемой скорректированной математической модели процесса прошивки и реологии металла сила осевого сопротивления оправки с учетом трения в очаге деформации, при известных предположениях [27, 28], определяется
F = 2лf ^ [f'(x) + ^о f(x)dx
(40)
где p - давление металла в очаге деформации при прошивке трубной заготовки; -коэффициент трения между прошиваемым металлом и оправкой в очаге деформации.
Предлагается следующая схема решения задачи. Сначала рассматривается задача в предположении, что возмущенная область состоит из пластической зоны и зоны упругости. При этом сначала решается уравнение для определения параметра Р (относительной
скорости распространения волны пластичности в металле), а затем по формулам (8) и (9) восстанавливаются истинные поля напряжений и скоростей. Далее проверяется выполнение условия се > -с0 во всей возмущенной области очага деформации. Если это условие выполняется, то указанным решением задачи ограничиваемся. При опровержении данного условия переходим к следующему этапу решения задачи. Далее задавшись определенными начальными и граничными условиями активных вынужденных гармонических колебаний технологического инструмента (оправки или иглы) в виде f (r) = r (t)sin(ra t) и
F(r) = c0 (t) cos(ra t) заданной частоты ю , математическим моделированием определяем эффективные параметры вибрационной прошивки трубной заготовки на этапе назначения технологического процесса.
Необходимо подчеркнуть, что поскольку процессы прошивки в рамках выбранной динамической модели очага деформации представляют волновые процессы, то путем приложения внешних активных управляющих воздействий, например, высокочастотных колебаний определенной частоты ю со стороны виб-рируемого технологического инструмента прошиваемому металлу (см. рис. 1), возможно добиться ускоренного внедрения оправки в металл и увеличение эффективности производства.
Отметим, что поскольку высокочастотный вибратор непосредственно встроен в иглу, то в процессе прессования трубы, кроме всего прочего, для защиты от тепловых воздействий он подвергается интенсивному охлаждению существующей системой охлаждения самой иглы. Под воздействием активных управляющих высокочастотных колебаний достигается значительное уменьшение предела упругости и текучести прессуемого металла, что обеспечивает снижение усилия прессования и приводит к повышению эффективности процесса прессования труб. При этом для стабилизации энергосиловых параметров очага деформации математическим моделированием могут быть определены требуемые рациональные амплитудно-частотные характеристики вынужденных колебаний технологических инструментов (матрицы и иглы). Типовые осциллограммы процесса прессования экспандированой трубной заготовки на трубопрофильном прессе 50 МН фирмы SMS MEER (заготовка 0180 мм,
k
материал 12Х18Н10Т, готовая труба 0105x10) представлен на рис. 3.
Стабилизацией параметров объемного напряженно-деформированного состояния металла, например, при прессовании трубной заготовки из стали 12Х18Н10Т диаметром 180 мм на трубопрофильном прессе 50 МН, применением активных высокочастотных вибраций частоты га = 7,5 кГц достигается снижение силы прессования с 17,9 до 15 МН (на 19 %) и за счет снижения значительных сил контактного трения повышение стойкости матриц и игл в два раза.
Необходимо особо подчеркнуть тот факт, что с применением высокочастотных активных вибраций, кроме всего прочего, под воздействием вынужденных колебаний системы значительно возрастает скорость расширения границы пластической зоны.
Полученные результаты показывают, что технологические процессы прессования труб из различных марок сталей и их сплавов с применением высокочастотной вибрации по рассматриваемой схеме достаточно корректно моделируются представленной математической моделью.
Разработанная математическая модель динамики процесса прессования труб определенно расширяет возможности анализа динамических явлений в очаге деформации. При этом решение задачи по прогнозированию напряженно-деформированного состояния очага деформации, как видно из выбранной динамической модели, базируется на реализации доступного математического аппарата. Полученные результаты позволяют реализовать управляемые режимы деформации в очаге и совершенствовать технологические процессы прессования труб на прошивном прессе. Уточненные динамические особенности очага деформации существенно дополняют существующие представления о напряженно-деформированном состоянии металла в условиях вибрационного прессования бесшовных труб.
Предложенная методика исследования технологического процесса прессования трубных заготовок в полые гильзы из различных марок сталей и их сплавов на профильном прессе по данной схеме заметно расширяет возможности анализа напряженно-деформированного состояния металла в очаге деформации. Полученные результаты указывают на высокую эффективность применения процесса вибрационного прессования и позволяют наметить
пути совершенствования технологического процесса производства бесшовных труб.
Выводы
Разработана динамическая модель процесса прессования труб на профильном прессе. Определены параметры возмущенной технологическими инструментами области очага деформации, соответствующей границам движения пластической и упругой зон, на границах взаимодействия матрицы и иглы в деформируемой трубной заготовке.
Составлена математическая модель распространения волновых процессов в прессуемом металле. Составлены дифференциальные уравнения для распространения фронтов пластической деформации и упругой зоны в возмущенной области очага деформации.
Получены аналитические зависимости для скорости распространения фронта пластической деформации и фронта распространения упругой зоны в возмущенной области очага деформации. Определены основные параметры границ сопряженных зон распространения волновых процессов в пластической и упругой областях соответственно. Анализом принятых моделей задачи о распространении волн в прессуемом металле определены динамические характеристики с учетом применения высокоэффективной технологии вибрационного процесса прессования труб (вибрация матрицы и иглы), что существенно дополняет ранее известные представления о напряженно-деформированном состоянии очага деформации. Выбором рациональных геометрических параметров образующей рабочей поверхности как самих матриц, так и амплитудно-частотных характеристик технологических инструментов (матрицы и иглы) при вибрационном прессовании на прессе достигнуто значительное повышение качества труб.
Литература
1. Перлин, И.Л. Теория прессования металлов / И.Л. Перлин. - М.: Металлургия, 1964. - 344 с.
2. Гуляев, Г.И. Прессование стальных труб и профилей / Г.И. Гуляев, А.Е. Притома-нов, О.П. Дробич. - М.: Металлургия, 1973. -192 с.
3. Медведев, М.И. Совершенствование процесса прессования труб / М.И. Медведев, Ю.Г. Гуляев, С.А. Чукмасов. - М. : Металлургия, 1986. -151 с.
4. Коликов, А.П. Новые процессы деформации металлов и сплавов / А.П. Коликов, П.И. Полухин, А.В. Крупин. - М.: Высшая школа, 1986. - 351 с.
5. Технология трубного производства: учеб. для вузов / В.Н. Данченко, А.П. Коликов, Б.А. Романцев, С.В. Самусев. - М.: Интермет Инжиниринг, 2002. - 640 с.
6. Боборыкин, Ю.А. Некоторые вопросы обработки металлов давлением с применением вибраций / Ю.А. Боборыкин // Материалы семинара «Применение вибраций для интенсификации процессов штамповки». - М.: МДННТП, 1961.
7. Иванова, В.С. Механизм пластической деформации при действии циклических нагрузок / В.С. Иванова //Металловедение и термическая обработка металлов. - 1960. - Vol. 2, no. 4. - P. 215-218.
8. Бухвалов, Г.Б. Обработка металлов в условиях вибрационного нагружения / Г.Б. Бухвалов, Н.П. Скрябин, В.П. Волегов // Доклады четвертой научно-технической конференции молодых специалистов. - М. : ВНИИМЕТМАШ, 1961.
9. Залесский, В.И. Исследование деформации при осадке в условиях вибрационного на-гружения / В.И. Залесский, И.П. Волков // Изв. ВУЗов. Черная металлургия. - 1965. - № 9.
10. Залесский, В.И. Вибрационное прессование с применением гидровинтового вибратора / В.И. Залесский, О.С. Мендыбаев // Изв. ВУЗов. Черная металлургия. - 1965. - № 11.
11. Потураев, В.Н. Некоторые результаты исследования вибропрессования / В.Н. Потураев, А.Ф. Миронюк. - М., 1966.
12. Стеблюк, В.И. Определение параметров вибрационных прессов для штамповой зачистки деталей / В.И. Стеблюк, И.П. Тар-таковский // Технология и автоматизация машиностроения. - Киев, 1968. - Вып. 3.
13. Зайцев, В.С. Геометрические условия вибрационной прокатки / В.С. Зайцев, И.М. Павлов, С.Ф. Буханов // Изв. ВУЗов. Черная металлургия. - 1969. - № 5.
14. Матвеев, И.Б. Выбор принципиальной схемы вибрационного пресса / И.Б. Матвеев // Кузнечно-штамповочное производство. -1973. - № 8.
15. Николаев, В.А. Давление и момент при прокатке полосы с вибрационным нагруже-нием / В.А. Николаев // Изв. ВУЗов. Черная металлургия. - 1971. - № 9.
16. Сверденко, В.П. Изучение влияния низкочастотных колебаний на силовые параметры при вибрационной осадке / В.П. Свер-денко, И.Г. Добровольский, В.Н. Булах // Изв. ВУЗов. Черная металлургия. - 1971. - № 7.
17. Thompson, F.C. The Fundamentals of Wire Drawing / F. C. Thompson // The Iron Age. -1938. - No. 15.
18. Patent Germany, no. 930201. The Wire Drawing. -1955.
19. Stancowic, H. Static-Dynamic Forging / H. Stancowic // Metal Treatment and Drop Forging. - 1958. - No. 159.
20. Solkowski, T. Étude du tréfilage avec filiers osillantes / T. Solkowski, A. Vilain // Metallurit (Belg.). - 1968. - Vol. 8, no. 4.
21. Winsper, C.T. The Superposition of Longitudinal Sonic Oscillations on the Wire Drawing Process / C.T. Winsper, D.H. Sansome // Proc. Inst. Mech. Engrs.. - 1968. - No. 25.
22. Клименко, В.М. Вибрационная обработка давлением / В.М. Клименко, В.Н. Ша-повал. - Киев: Техника, 1977. - 128 с.
23. Малинин, Н.Н. Прикладная теория пластичности и ползучести / Н.Н. Малинин. -М. : Машиностроение, 1975. - 400 с.
24. Безухов, Н.И. Основы теории упругости, пластичности и ползучести / Н.И. Безухов. - М. : Высшая школа, 1978. - 512 с.
25. Ландау, Л.Д. Теоретическая физика. В 10 т. Т. 7: Теория упругости / Л.Д. Ландау, Е.М. Лифшиц. - М. : Наука, 1987. - 248 с.
26. Кошляков, Н.С. Уравнения в частных производных математической физики / Н.С. Кошляков, Э.Б. Глинер, М.И. Смирнов // М. : Высшая школа, 1970. - 712 с.
27. Рахманов, С.Р. Динамические процессы при прошивке трубной заготовки на прессе / С.Р. Рахманов, О.Г. Гоман // Сталь. -2010. - № 6. - С. 73-78.
28. Рахманов, С.Р. Математическое моделирование процессов вибрационной прошивки трубной заготовки на прошивном прессе / С.Р. Рахманов // Вибрация в технике и технологиях. - 2010. - № 3. - С. 99-108.
Рахманов Сулейман Рахманович, канд. техн. наук, доцент кафедры теоретической механики и сопротивления материалов, Национальная металлургическая академия Украины, г. Днепропетровск; npfvostokplus@gmail.com.
Поступила в редакцию 20 марта 2017 г.
DOI: 10.14529/met170210
MATHEMATICAL MODELLING OF THE PROCESSES OF VIBRATIONAL PIPE PRESSING
S.R. Rakhmanov, npfvostokplus@gmail.com
National Metallurgical Academy of Ukraine, Dnepropetrovsk, Ukraine
The mathematical model of dynamics of process of pressing of an axial-symmetric pipe preparation on a parsing press of tube press installation is developed. For the accepted settlement scheme, dynamic model and process of pressing of pipe preparation analytical expressions of expansion of the borders corresponding to plastic and elastic zones in the deformation center are received. The analytical dependence characterizing the speed of distribution of the front of plastic deformation on speed of the pressed metal is received. Dynamic characteristics of emergence of the wave phenomena in the pressed metal in the conditions of influence from the tool are defined that significantly supplements earlier known ideas of the stress-strained condition of metal in the deformation center. Fields of deformations in zones of the disturbed area of the center of deformation taking into account influence of technological tools are established. By modeling of dynamic processes in the pressed metal features of realization of highly effective technology of the operated vibration pressing of pipes on a profile press of tube press installation are established. The choice of rational parameters (the amplitude-frequency characteristic) of the used active process of vibration pressing of pipe preparation has reached substantial increase of efficiency of process and quality of the pierced preparations.
Keywords: pipe; press; matrix; metal; needle; caliber; deformation center; stress; deformation; elasticity modulus; shear modulus; plastic area; plasticity conditions; elastic area; vibration; frequency; amplitude; period.
References
1. Perlin I.L. Teoriya pressovaniya metallov [Theory of Metal Pressing]. Moscow, Metallurgiya Publ., 1964. 344 p.
2. Gulyaev G.I., Pritomanov A.E., Drobich O.P. Pressovanie stal'nykh trub i profiley [Pressing of Steel Pipes and Profiles]. Moscow, Metallurgiya Publ., 1973. 192 p.
3. Medvedev M.I., Gulyaev Yu.G., Chukmasov S.A. Sovershenstvovanie protsessa pressovaniya trub [Improvement of Pipe Pressing Process]. Moscow, Metallurgiya Publ., 1986. 151 p.
4. Kolikov A.P., Polukhin P.I., Krupin A.V. Novye protsessy deformatsii metallov i splavov [New Processes of Deformation of Metals and Alloys]. Moscow, Vysshaya shkola Publ., 1986. 351 p.
5. Danchenko V.N., Kolikov A.P., Romantsev B.A., Samusev S.V. Tekhnologiya trubnogo proizvodstva [Technology of Pipe Production]. Moscow, Intermet Engineering Publ., 2002. 640 p.
6. Boborykin Yu.A. [Some Issues of Metal Forming Using Vibration]. Materialy seminara "Primenenie vibratsiy dlya intensifikatsii protsessov shtampovki" [Materials of the Seminar "Application of Vibration for Intensification of Stamping Processes"]. Moscow, MDNTP Publ., 1961. (in Russ.)
7. Ivanova V.S. Mechanism of Plastic Deformation in Cyclic Loading. Metal Science and Heat Treatment, 1960, vol. 2, no. 4, pp. 215-218. DOI: 10.1007/BF00681176
8. Bukhvalov G.B., Skryabin N.P., Volegov V.P. [Metal Forming Under Vibrational Loading]. Doklady chetvertoy nauchno-tekhnicheskoy konferentsii molodykh spetsialistov [Reports of the Fourth Scientific and Technical Conference of Young Professionals]. Moscow, VNIIMETMASh Publ., 1961.
9. Zalesskiy V.I., Volkov I.P. [Investigation of Deformation in Compression Under Vibrational Loading]. Izvestiya VUZ. Chernaya metallurgiya, 1965, no. 9. (in Russ.)
10. Zalesskiy V.I., Mendybaev O.S. [Vibrational Pressing Using a Hydro-Screw Vibrator]. Izvestiya VUZ. Chernaya metallurgiya, 1965, no. 11. (in Russ.)
11. Poturaev V.N., Mironyuk A.F. Nekotorye rezul'taty issledovaniya vibropressovaniya [Some Results of Investigation of Vibrational Pressing]. Moscow, 1966.
12. Steblyuk V.I., Tartakovskiy I.P. [Determination of Parameters of Vibrational Presses for Stamping Finish of Components]. Tekhnologiya i avtomatizatsiya mashinostroeniya [Technology and Automatization of Machine Building]. No. 3. Kiev, 1968. (in Russ.)
13. Zaytsev V.S., Pavlov I.M., Bukhanov S.F. [Geometrical Conditions of Vibrational Rolling]. Izvestiya VUZ. Chernaya metallurgiya, 1969, no. 5. (in Russ.)
14. Matveev I.B. [Choosing the Principal Scheme of a Vibrational Press]. Kuznechno shtampo-vochnoe proizvodstvo, 1973, no. 8. (in Russ.)
15. Nikolaev V.A. [Pressure and Moment in Band Rolling with Vibrational Loading]. Izvestiya VUZ. Chernaya metallurgiya, 1971, no. 9. (in Russ.)
16. Sverdenko V.P., Dobrovol'skiy I.G., Bulakh V.N. [Investigation of the Effect of Low-Frequency Oscillations on the Force Parameters in Vibrational Compression]. Izvestiya VUZ. Chernaya metallurgiya, 1971, no. 7. (in Russ.)
17. Thompson F.C. The Fundamentals of Wire Drawing. The Iron Age, 1938, no. 15.
18. The Wire Drawing. Patent Germany, no. 930201, 1955.
19. Stancowic H. Static-Dynamic Forging. Metal Treatment and Drop Forging, 1958, no. 159.
20. Solkowski T., Vilain A. Étude du tréfilage avec filiers osillantes. Metallurit (Belg.), 1968, vol. 8, no. 4.
21. Winsper C.T., Sansome D. H. The Superposition of Longitudinal Sonic Oscillations on the Wire Drawing Process. Proc. Inst. Mech. Engrs., 1968, no. 25.
22. Klimenko V.M., Shapoval V.N. Vibratsionnaya obrabotka davleniem [Vibrational Metal Forming]. Kiev, Tekhnika Publ., 1977. 128 p.
23. Malinin N.N. Prikladnaya teoriya plastichnosti i polzuchesti [Applied Theory of Plasticity and Creep]. Moscow, Mashinostroenie Publ., 1975. 400 p.
24. Bezukhov N.I. Osnovy teorii uprugosti, plastichnosti i polzuchesti [Fundamentals of the Theory of Elasticity, Plasticity and Creep]. Moscow, Vysshaya shkola Publ., 1978. 512 p.
25. Landau L.D., Lifshits E.M. Teoreticheskaya fizika. T. 7. Teoriya uprugosti [Theoretical Physics. Vol. 7. Elasticity Theory]. Moscow, Nauka Publ., 1987. 248 p.
26. Koshlyakov N. S., Gliner E. B., Smirnov M. I. Uravneniya v chastnykh proizvodnykh matematicheskoy fiziki [Partial Derivative Equations of Mathematical Physics]. Moscow, Vysshaya shkola Publ., 1970. 712 p.
27. Rakhmanov S.R., Goman O.G. [Dynamic Processes in Press Piercing of a Pipe Billet]. Stal', 2010, no. 6, pp. 73-78. (in Russ.)
28. Rakhmanov S.R. [Mathematical Modelling of Vibrational Piercing of a Tubular Billet with a Piercing Press]. Vibratsiya v tekhnike i tekhnologiyakh, 2010, no. 3, pp. 99-108. (in Russ.)
Received 20 March 2017
ОБРАЗЕЦ ЦИТИРОВАНИЯ
Рахманов, С.Р. Математическое моделирование процессов вибрационного прессования труб / С.Р. Рахманов // Вестник ЮУрГУ. Серия «Металлургия». -2017. - Т. 17, № 2. - С. 77-90. DOI: 10.14529/met170210
FOR CITATION
Rakhmanov S.R. Mathematical Modelling of the Processes of Vibrational Pipe Pressing. Bulletin of the South Ural State University. Ser. Metallurgy, 2017, vol. 17, no. 2, pp. 77-90. (in Russ.) DOI: 10.14529/met170210