close

Вход

Забыли?

вход по аккаунту

?

Модель динамической системы шлифовального станка для обработки деталей подшипников транспортной техники и технологического оборудования..pdf

код для вставкиСкачать
Вестник СГТУ. 2013 №2 (71). Выпуск 2
УДК 621.923
В.А. Каракозова, А.И. Зорин, А.А. Игнатьев
МОДЕЛЬ ДИНАМИЧЕСКОЙ СИСТЕМЫ ШЛИФОВАЛЬНОГО СТАНКА
ДЛЯ ОБРАБОТКИ ДЕТАЛЕЙ ПОДШИПНИКОВ ТРАНСПОРТНОЙ ТЕХНИКИ
И ТЕХНОЛОГИЧЕСКОГО ОБОРУДОВАНИЯ
Предлагается уточненная модель съема припуска, используемая для построения
динамической системы станка.
Шлифование, подача круга, передаточная функция, динамические характеристики, модель динамической системы, съем припуска
V.A. Karakozova, A.I. Zorin, A.A. Ignatyev
MODEL OF DYNAMIC SYSTEM OF THE GRINDER FOR PROCESSING OF DETAILS
OF BEARINGS OF TRANSPORT EQUIPMENT AND PROCESSING EQUIPMENT
The specified model effort of cutting, used for creation of dynamic system of the machine is offered.
Grinding, supply range, the transfer function, dynamic characteristics, the model of a
dynamic system, effort of cutting
Для обеспечения качества колец при обработке на современных шлифовальных автоматах
необходимо учитывать, как достаточно известные методы технологии машиностроения, так и методы
теории управления. Последние позволяют не только эффективно выполнить анализ влияния
динамических характеристик станков на результаты обработки, но и учесть воздействие как
194
Вестник СГТУ. 2013 №2 (71). Выпуск 2
детерминированных, так и стохастических составляющих процессов в (технологической системе) ТС,
а также реализовать управление качеством на основе нескольких информационных параметров, в том
числе и по принципу обратной связи с применением приборов активного контроля [1, 2]. Поэтому очень
важно реализовать в производстве методики управления шлифованием на основе оценки динамического
состояния станка по стохастическим характеристикам и контролю качества обработки колец.
Разработка системы автоматического управления (САУ) процессами шлифования связана с
наличием достоверной математической модели процесса как объекта управления, который можно
описать дифференциальными или операторными уравнениями, передаточными, переходными или
частотными функциями [3]. На практике можно построить как теоретическую модель динамической
системы (ДС), так и экспериментально-аналитическую модель, применяя методы теории
идентификации [4].
В данной статье решена задача построения модели ДС в виде передаточной функции,
характерным отличием которой является разработка уточненной модели съема припуска,
Построению модели ДС предшествует составление функциональной схемы [5,6].
Для установившегося процесса шлифования сила резания представляется формулой
F = f1 ( k реж , q, v k , b, h ),
(1)
где kреж – коэффициент, характеризующий режущую способность круга; q – коэффициент,
характеризующий физико-механические свойства обрабатываемого материала; vk – скорость резания;
b – ширина шлифования, равная при врезном шлифовании ширине обрабатываемой детали; h –
глубина резания.
Кроме того, на динамику процесса резания при врезном шлифовании будут оказывать
влияние площадь S контактирующей (активной) режущей поверхности круга и детали (рис. 1),
которая определяется размерами ширины шлифования b и длиной дуги активного контакта l.
S = b⋅l ,
(2)
От колебаний величины S зависят
колебания усилия резания F, которые оказывают
влияние на температуру в зоне резания и скорость
съема припуска.
При постоянном значении b площадь S
контактирующей (активной) режущей поверхности
круга зависит от длины дуги контакта l.
Применяемые в шлифовании формулы для
определения длины дуги контакта, в основном,
эмпирического происхождения, не связывающие
размеры круга и отверстия детали, что снижает
точность вычислений по этим формулам, например
[7]
l = C1 ⋅ tФ0,5 ,
(3)
t
где ф – фактическая глубина шлифования.
Для выведения более точной формулы
длины дуги контакта l рассмотрим схему контакта
круга и детали при врезном шлифовании,
приведённую на рис. 1.
Найдём точки пересечения окружностей
абразивного круга и отверстия в детали. (A, C),
(рис. 1), для чего решим систему двух уравнений:
[X2 + (R2− r +2 h) )]2 + Y 2 = R 2 ,
X +Y =r
Рис. 1. Схема контакта круга и детали
при врезном шлифовании где: R – радиус
шлифуемого отверстия, r – радиус шлифовального
круга, X, Y – текущие координаты, l – длина дуги
активного контакта круга и детали, h – глубина
врезания круга в деталь за один оборот детали,
А, С – точки пересечения окружностей
абразивного круга и отверстия в детали
(4)
где X,Y – текущие координаты, R – радиус отверстия, r – радиус шлифовального круга.
Из совместного решения уравнений (4) следует:
XA =
C1 ⋅ (r − h )
,
C1 + h
YA =
2hrR
C1 + 2h
(5, 6)
195
Вестник СГТУ. 2013 №2 (71). Выпуск 2
где XA и YA координаты точки пересечения A( ХA, YA ) окружностей круга и детали, С1 = R − r .
При выведении формул (4) и (5) h2 было приравнено 0, как величина бесконечно малая
второго порядка по сравнению с другими входящими в формулу величинами.
В результате преобразований получена формула для расчета длины дуги контакта
абразивного круга с обрабатываемой деталью.
r
l = 2⋅
1 + (Y ') dx ,
∫
Y = r 2 − X 2 → Y'= −
2
C1⋅( r − h )
X
r2 − X 2
(7.,8)
C1 + h
где Y2 +X2 =r2.
Подставляя Y’ из формулы (8) в формулу (7), получим
X
l = r ⋅ arcsin
r
r
C1 ( r − h )
C
π
C ⋅ (r − h ) 
= r ⋅  − arcsin 1

r ⋅ (C1 + 2h ) 
2
(9)
Подставляя значения (8) в формулу (2) получим выражение для площади контактирующей
(активной) режущей поверхности шлифовального круга.
π
C ⋅ (r − h ) 
S = l ⋅ b = r ⋅ H ⋅  − arcsin 1

r ⋅ (C1 + 2h ) 
2
(10)
Формула (10) позволяет определить площадь контактирующей (активной) режущей
поверхности круга и детали (поверхности резания) при любых размерах круга и отверстия в детали,
полученных либо с датчиков, либо рассчитанных аналитически. Зная площадь поверхности резания
можно точно рассчитывать и контролировать усилие резания и адекватность процессу шлифования
величины поперечной подачи суппорта. В формуле величина h может приравниваться поперечной
подаче Vc .
Также одним из важных факторов, влияющих на аппаратный ход процесса управления
является мощность N, потребляемая двигателем привода вращения детали. Согласно [8], мощность N
N=V*F
(11)
где V – скорость (линейная), F – сила ( в нашем случае сила резания).
Уточним определение скорости V, которая в общем определяется траекторией перемещения
режущего элемента по шлифуемой поверхности детали и частотой (угловой скоростью) вращения
детали.
При выводе формул (8) и (9) за траекторию относительного перемещения активного
режущего элемента, как и в известной литературе по шлифованию, в т.ч. и по врезному, принималась
дуга окружности. В действительности, из-за непрерывности процесса врезного шлифования,
траекторией перемещения режущих элементов является кривая, отличная от окружности, больше
напоминающая спираль Архимеда.
С учетом формулы спирали Архимеда
ρ = K сп . ⋅ ϕ =
a
ϕ
2π ,
(12)
получена формула для определения величины текущего («мгновенного») радиуса шлифования r:
r = r0 + Vc ⋅ nt + К сп ⋅ ϕ = r0 + Vc ⋅ nt +
Vc
ϕ 

⋅ ϕ = r0 + Vc ⋅  nt +

2π
2π 

(13)
где r0 – радиус отверстия в заготовке (в момент начального соприкосновения круга и шлифуемой детали); a – шаг спирали Архимеда, принятый равным величине расчетной подачи суппорта Vc шлифовального станка за один оборот детали , nt – количество оборотов, сделанное деталью с момента
начального соприкосновения круга и шлифуемой детали до момента измерения радиуса отверстия.
Из формулы (13) видно, что текущий радиус r вращения режущего элемента непрерывно
увеличивается – за каждый оборот детали (φ = 2π) радиус вращения активного режущего элемента
будет увеличиваться на Vc. . Vc – поперечная подача суппорта на 1 оборот детали) благодаря врезанию
в деталь). С увеличением r будет непрерывно увеличиваться линейная скорость перемещения
абразивного режущего элемента по шлифуемой поверхности детали и связанная с ней потребляемая
196
Вестник СГТУ. 2013 №2 (71). Выпуск 2
мощность двигателя привода вращения детали. Согласно [8] линейная скорость при вращательном
движении материальной точки (в нашем случае абразивного режущего элемента
V =rω
(14)
Тогда для врезного шлифования линейная текущая (мгновенная) скорость перемещения
режущего активного элемента по шлифуемой поверхности детали будет равна с учетом формул (12) и
(13)

ϕ 

V = r ⋅ ω = r0 + Vc ⋅  nt +
 ⋅ω ,
2π 


где
(15)
ω – угловая скорость вращения детали.
Вводя в формулу мощности (11) коэффициент k шл – пропорциональности, зависящий от
специфики внутреннего врезного шлифования и размерности параметров, входящих в формулу и
учитывая формулу (14), получим формулу потребляемой мощности двигателя привода вращения
детали.

ϕ 

N = V ⋅ FY = K шл ⋅ r0 + Vc ⋅  nt +
 ⋅ FY ⋅ ω
2π 


(16)
Формула (2.20) позволяет по сравнению с известными формулами, например [8],
N заг = Pz ⋅ Vаг / (1000 ⋅ 60 ⋅ n2 )
(17)
более точно представлять и корректировать управление процессом врезного шлифования. В
практическом шлифовании принято считать, что для установившегося процесса, когда подача, усилие
резания и частота вращения детали не меняются, должна быть постоянной и мощность. В
действительности, согласно формуле (16), потребляемая мощность двигателем вращения детали
непрерывно возрастает от минимального до максимального значения в течение одного перехода
операции врезного шлифования.
В то же время, скорость резания при увеличении диаметра отверстия может быть принята для
данных размеров круга и детали при снятии припуска на чистовых переходах постоянной. Тогда при
врезном шлифовании при постоянной скорости резания, при неизменных режущей способности
круга и свойствах материала обрабатываемых деталей, при ширине шлифования b=const, колебания
площади и скорости контактирующей (активной) режущей поверхности круга и детали не
сказываются на глубине резания, сила резания F в установившемся режиме будет определяться лишь
фактической толщиной среза ад:
F = f 3 ( aд ).
(18)
Под толщиной среза ад здесь понимается некоторая приведенная (интегральная) толщина
срезаемой стружки металла, которая фактически обусловлена параметрами множества микростружек,
снимаемых элементарными зернами шлифовального круга в текущий момент времени.
В работе [6] автором рассматривается снимаемый с детали при врезном шлифовании припуск
(толщина среза) a0 ( p ) как сумма припусков a Д ( p ) – основного (расчетного) припуска, снимаемого
с детали, aИК (P ) – приращения припуска из-за износа круга, aУД ( p ) – упругой деформации
технологической системы и приводятся выведенные им формулы, т.е.
a0 ( p ) = a Д ( p ) + a ИК ( p ) + aУД ( p )
(19)
Несмотря на то, что съем припуска при врезном шлифовании рассматривается в схеме,
принятой Михелькевичем В.Н., происходящим по дуге окружности, а не по кривой Архимеда (в этом
случае расчетные формулы были бы более точными), основные закономерности, выведенные им для
определения толщины среза за один оборот детали, могут быть приняты для построения
динамической модели системы.
Так как при врезном шлифовании основная зависимость – это зависимость между силой
резания и скоростью поперечной подачи, то если за входную переменную шлифовального станка
xвх(t) принять функцию изменения скорости поперечной подачи vc(t), а за выходную переменную
xвых(t) – соответствующую ей функцию изменения радиальной силы Fy(t), то передаточная функция
ОУ будет равна:
197
Вестник СГТУ. 2013 №2 (71). Выпуск 2
Wc ( p ) =
Fy ( p )
Vc ( p )
.
(20)
С учетом обратной связи по упругой деформации системы и обратной связи по износу круга
получим обобщенную модель ДС при врезном шлифовании.
Для обобщенной структурной схемы процесса врезного внутреннего шлифования
доработаны базовые формулы [6], что позволило получить уточненную передаточную функцию ДС.
F ( p)
W ( p) = Y
=
Vc ( p )
(
K рез
⋅ 1− e
p
К рез
1+
р
(
⋅ 1− e
− pτ Д
− pτ Д
)⋅ (К
)
(21)
УС
⋅ р + К ИК )
После разложения в ряд Пада с точностью( 2/2) экспоненциальной функции ( 1 − e
гебраических преобразований формула (21) приняла вид:
W(p) =
pτ Д
12τ Д ⋅ К рез
) и ал-
(22)
τ Д2 р 2 + 6τ Д р + 12 + 12τ Д ⋅ K рез ⋅ (K УС ⋅ p + К ИК )
Учитывая, что ДС включает совокупность параллельно соединенных шпиндельный узел
(ШУ) детали и ШУ инструмента, следовательно WУС(р) = Wи(р) + Wд(p). Тогда обобщенная
структурная схема процесса врезного внутреннего шлифования будет иметь вид, приведенный на
(рис. 2), а соответствующая ей передаточная функция преобразуется к виду:
F ( p)
W ( p) = Y
=
Vc ( p )
K рез
p
1+
К рез
р
(
⋅ 1− e
− pτ Д
)⋅ [К
(
⋅ 1− e
УС
− pτ Д
)
(23)
⋅ р ⋅ (WИ ( р ) + W Д ( p ))+ K ИК ]
Рис. 2. Структурная схема динамической системы с учетом уточненной схемы съема припуска
Передаточные функции ШУ детали WД(р) и ШУ инструмента WИ(р) в общем случае имеют
сложную структуру, образованную совокупностью колебательных звеньев [6]. Выражения для них
имеют вид
h Дi
n
W Д ( р) = ∑
i =1
T p + 2γ Дi T Дi p + 1
2
Дi
2
hИj
m
,
WИ ( р) = ∑
j =1
T p + 2γ Иj TИj p + 1
2
Иj
2
,
(24),( 25)
где hДi , hИj – коэффициенты, обусловленные статической жесткостью отдельных элементов ШУ; ТДi ,
TИj – постоянные времени, обусловленные собственными частотами отдельных элементов ШУ; γДi ,
γИi – относительные коэффициенты демпфирования отдельных элементов ШУ.
198
Вестник СГТУ. 2013 №2 (71). Выпуск 2
Анализ формулы (23) с учетом выражений (24) и (25) достаточно сложен, поэтому с целью
упрощения передаточной функции WP(p) ШУ инструмента и ШУ детали рассматриваются как
колебательные звенья с одной основной частотой, а именно:
WИ ( р) =
hИ
,
2 2
TИ p + 2γ И TИ p + 1
WД ( p) =
hД
TД2 p 2 + 2 y Д T Д p + 1
(26)
Подставив значения (26) в (23), получим:
W ( p) =
12τ Д Kрез
1
(27)
τ Д2 p2 + 6τ Д p + 12 + 12τ Д Kрез 


h
h
Д
+ KИК ]
+ 2 2
КУС р 2 2 И

T
p
+
2
γ
T
р
+
1
T
p
+
2
y
T
p
+
1
И И
Д
Д Д
 И


Или в форме многочлена:
W ( p) =
a 3 p 3 + a2 p 2 + a1 p + a0
b4 p 4 + b3 p 3 + b2 p 2 + b1 p + b0
(28)
где аi – вj- – постоянные коэффициенты, вычисляемые по коэффициентам передаточных функций
звеньев ДС.
Отличительной особенностью врезного шлифования с точки зрения управления является
инерционность воздействия связанных с ним физических процессов на параметры обработки.
Например, после кратковременного действия дискретного источника тепла в любом звене системы,
вызываемые им температурные изменения в детали (увеличение размеров) могут идти сравнительно
длительное время уже после прекращения функционирования источника тепла. То же может быть
сказано и о реакции управляющей системы на колебания съема металла при неравномерном
припуске, которая может последовать также не сразу и т.д. Поэтому для инерционных процессов
обработки (а это относится не только к шлифованию, но и к любому управляемому
технологическому процессу металлообработки), режимы обработки рассчитываются, как правило, за
определенный промежуток времени, например, подача мм в минуту или мм за 1 оборот. Для учета
подобных допущений в моделях управляющих систем вводятся элементы запаздывания в структуру
ДС, что повышает точность и адекватность описания подобных систем управления и их
передаточных функций.
В нашем случае введение в структуру ДС передаточных функций WИ(р) и WД(р) повысило
порядок характеристического уравнения со второго до четвертого, что способствовало повышению
точности и адекватности описания ДС.
Построенная модель съема припуска используется для выбора режима шлифования при
обработке колец подшипников [9].
ЛИТЕРАТУРА
1. Игнатьев С.А. Мониторинг технологического процесса как элемент системы управления
качеством продукции / С.А. Игнатьев, В.В. Горбунов, А.А. Игнатьев. Саратов: СГТУ, 2009. 160 с.
2. Каракозова В.А. Экспериментально-аналитический метод определения передаточной функции динамической системы станка / В.А. Каракозова, А.А. Игнатьев // Автоматизация и управление в
машино- и приборостроении: сб. науч. тр. Саратов: Изд-во СГТУ, 2010. С. 90-92.
3. Егоров К.В. Основы теории автоматического регулирования. / К.В. Егоров. М.: Энергия,
1967. 648 с.
4. Игнатьев А.А. Основы теории идентификации объектов управления / А.А.Игнатьев,
С.А. Игнатьев Саратов: СГТУ, 2008. 44 с.
5. Бржозовский Б.М. Модель динамической системы шлифовального станка с учетом стохастичности процессов / Б.М. Бржозовский, С.А. Игнатьев // Автоматизация и управление в машино- и
приборостроении: сб. науч. тр. Саратов: СГТУ, 2001. С. 31-34.
6. Михелькевич В.Н. Автоматическое управление шлифованием / В.Н. Михелькевич. М.: Машиностроение, 1975. 34 с.
7. Зубарев Ю.М. Математическое описание процесса шлифования. Санкт-Петербургский институт машиностроения (ЛМЗ-ВТУЗ) Источник: Российский научно-технический журнал «Инструмент и технологии №17-18] 2004 г. С. 157-161.
199
Вестник СГТУ. 2013 №2 (71). Выпуск 2
8. Яворский Б.М. Справочник по физике / Б.М. Яворский, А.А. Детлаф. М., Наука, 1968. 940 с.
9. Каракозова В.А., Игнатьев А.А. Выбор подачи круга для обеспечения качества обработки
колец подшипников на основе идентификации динамической системы шлифовального станка / В.А.
Каракозова, А.А. Игнатьев // Вестник Саратовского государственного технического университета
2011. № 2 (56). С. 69-73.
Игнатьев Александр Анатольевич –
доктор технических наук, профессор,
заведующий кафедрой «Автоматизация
и управление технологическими процессами»
Саратовского государственного технического
университета имени Гагарина Ю.А.
Aleksandr A. Ignatyev –
Dr.Sc.Tech., the professor
the chair chief «Automation and management of
technological processes»
Gagarin Saratov state technical university
Каракозова Вера Алексеевна –
кандидат технических наук, доцент, заведующая
лабораторией кафедры «Автоматизация
и управление технологическими процессами»
Саратовского государственного технического
университета имени Гагарина Ю.А.
Vera A. Karakozova –
Ph.D., associate professor, the laboratory chief
the chair «Automation and management of
technological processes»
Gagarin Saratov state technical university
Зорин Анатолий Иванович –
инженер первой категории лаборатории
кафедры «Автоматизация и управление
технологическими процессами» Саратовского
государственного технического университета
имени Гагарина Ю.А.
Anatoliy I. Zorin –
Engineer of the first category of laboratory
the chair «Automation and management of
technological processes»
Gagarin Saratov state technical university
Статья поступила в редакцию 03.04.13, принята к опубликованию 30.04.13
200
1/--страниц
Пожаловаться на содержимое документа