книги / Моделирование технологического оборудования
..pdfГидросуппорт представляет собой следящую систему с однокромочным золотниковым распределителем. Основными узлами являются: гидростанция Г48-1; 1 -напорный трубопровод, подающий рабочую жидкость от станции к гидроцилиндру; 2,5 - задняя и передняя полости гидроцилиндра; 4 - отверстие в поршне, соединяющее полости; 5 - неподвижный поршень гидроцилиндра; б - подвижная часть суппорта с корпусом гидроцилиндра; 7 - канал, соединяю щий переднюю полость с однокромочным золотником; 8 - корпус следящего золотника; 9 -плунжер следящего золотника; 10 -рычаг обратной связи; 11 -ко пир; 12 - сливной трубопровод.
Для разработки математической модели необходимо разбить гидросуп порт на три подсистемы: гидравлическую, механическую и гидравлическую, которые взаимодействуют друг с другом через полости гидроцилиндра. Меха ническая цепь гидросуппорта изображена на рис. 4.7.
Элементы механической цепи отражают физические явления, происхо дящие в гидромеханической системе. Работа гидростанции представлена источ ником постоянного давления Рн, то есть считается, что гидроклапаны идеальны, а пульсацией давления от работы насоса можно пренебречь. Элемент RTотра жает сопротивление напорного трубопровода течению жидкости, элементы Си Сз учитывают податливость жидкости в трубопроводах и самих трубопроводов. Элемент /?л учитывает сопротивление отверстия в поршне перетеканию жидко сти из одной полости в другую, оказывающее важнейшее влияние на работу гидросуппорта, а элемент N - нелинейный закон истечения жидкости через ко льцевую щель золотника. Элемент SUM| представляет собой сумматор, где складываются движения корпуса и плунжера золотника. Результат этого сложе ния определяет величину открытия кольцевой щели золотника. Элемент М от ражает инерционность механической части суппорта, элемент F - внешнюю си лу, в качестве которой выступает сила резания, элементы S\, S2 - площади поло стей гидроцилиндра, через которые осуществляется взаимодействие между под системами.
На основе механической цепи строится эквивалентная схема (рис. 4.8).
Для получения |
эквивалентной |
схемы |
подсистем |
|
гидравлических |
втор3* |
|||
использована первая система аналогий, а механической подсистеМь1 |
TJjna |
|||
система аналогий. Фиктивные источники JF отражают связь гиратОрн |
|
|||
между подсистемами. |
основе метода |
узловых |
„магическая |
|
Полученная на |
потенциалов матем*** |
|
модель отражает динамические процессы в гидросуппорте ГСП-41 и м быть использована для анализа его работы. Структура исходных данных, подготовленная для анализа динамики гидросуппорта с помощью программы
PAN, |
представлена в табл. 4.3. |
|
|
|
|||
|
|
|
|
|
|
|
Таблица 4.3 |
|
|
Исходные данные системы гидросуппорта ГСП - 41 |
|
||||
Ветвь |
|
Узел |
|
|
Значение |
|
|
№ |
тип |
от |
до |
вли |
C,R,L |
E,J,k |
знак |
|
|
|
|
яния |
|
|
|
1 |
ER |
1 |
0 |
|
Лт |
р» |
|
2 |
С |
1 |
0 |
|
с, |
S, |
|
3 |
JF |
1 |
0 |
2 |
|
|
|
4 |
JF |
2 |
0 |
1 |
м |
S, |
-1 |
5 |
С |
2 |
0 |
|
F |
|
|
6 |
J |
2 |
0 |
|
|
-1 |
|
"7 |
JF |
2 |
0 |
3 |
|
s 2 |
-1 |
8 |
JF |
3 |
0 |
2 |
с2 |
& |
|
9 |
С |
3 |
0 |
|
|
|
|
10 |
NOTR |
3 |
0 |
4 |
N |
|
|
11 |
R |
1 |
3 |
|
R:, |
|
|
12 |
SUM |
4 |
0 |
2 |
к = 1 |
0 |
|
Рис. 4.10. Расчетная {а) и эквивалентная (5) схемы движения рабочего органа, отображающие демпфирование процесса трения
В расчетной схеме элемент Сх отражает жесткость связи, удерживающей стол в направлении оси х; Их- коэффициент сил сопротивления* по оси х, опре деляемый демпфированием в связи; Ит- коэффициент сопротивления по оси х, определяемый силами вязкого трения смазки; М - инерционность стола.
На основе второй системы аналогий получаем эквивалентную схему дви жения стола (рис. 4.10,6).
Используя метод узловых потенциалов, по эквивалентной схеме строим математическую модель, которая представляет собой дифференциальное урав
нение: |
|
М х+ (Их+ /?т )х+ Ctx = F. |
(4.1) |
Коэффициенты в уравнении описаны выше.
Наиболее часто, когда возникает задача исследования плавности движе ния рабочего органа станка, процесс трения отображается в модели как аперио дическое звено, имеющее передаточную функцию
fVAS) = *т |
(4.2) |
Г,5 + Г |
|
где кл - статическая характеристика контактного трения; Гт - постоянная времени процесса трения.
При моделировании методом прямой аналогии в этом случае процесс трения рассматривается как отдельная подсистема, взаимодействующая с меха нической подсистемой. Механическая цепь привода стола с учетом процесса трения, отраженного в виде апериодического звена, изображена на рис. 4.11.
Рис. 4.11. Механическая цепь привода рабочего органа с учетом процесса трения, отраженного в виде апериодического звена
Здесь элементом к отображено взаимодействие между механической под системой и процессом трения. Используя вторую систему аналогий, получаем эквивалентную схему привода исполнительного органа с учетом процесса тре ния как апериодического звена (рис. 4.12).
|
х |
|
т |
~ |
|
|
м - - F /k |
= 1 |
V |
1 |
PF |
Рис. 4.12. Эквивалентная схема привода рабочего органа с учетом процесса грсния. отраженного в виде апериодического звена
Взаимодействие между подсистемами в эквивалентной схеме учтено включением двух фиктивных источников: JF, PF, которые обеспечивают скоро стную связь:
PF =k F r,
(4.3)
JF = kT x.
Используя метод узловых потенциалов теории электротехники, получаем математическую модель, описывающую динамику привода рабочего органа с учетом процесса трения как апериодического звена:
|
М х+ hx х+ Схх - F + FT=0, |
(4.4) |
|
|
|
|
hTFT +Fr +kTx = 0, |
|
где М - |
масса стола; |
|
С, - |
жесткость связи рабочего органа по оси х\ |
|
Их- |
коэффициент демпфирования в связи по оси х; |
|
F - |
внешнее воздействие по оси х\ |
|
F, - |
сила трения; |
|
kt - |
статическая характеристика трения. |
|
В приводах рабочих органов прецизионных станков на точность обработ ки оказывает влияние как плавность движения, так и величина всплывания. По этому при перемещении рабочего органа необходимо исследовать влияние обо их этих параметров. В этом случае привод подач должен рассматриваться в пространстве двух координат. Расчетная схема привода в пространстве двух координат представлена на рис. 4.13.
Рис. 4.13. Расчетная схема привода рабочего органа в пространстве двух координат
При использовании метода прямой аналогии каждая из координат рас сматривается как отдельная подсистема. Вторая система аналогий отображения дает эквивалентную схему, показанную на рис. 4.14.
х |
у |
Рис. 4.14. Эквивалентная схема привода рабочего органа в пространстве двух координат
Фиктивные источники в эквивалентной схеме отражают взаимодействие между подсистемами, то есть перемещение рабочего органа по координатам л и у. Эта связь обусловлена процессами трения в направляющих. Выражения для фиктивных источников имеют следующий вид:
PF = кт- у,
(4.5)
JF = кд х.
Эти уравнения связи отражают условия смешанного трения.
При учете влияния жидкостного трения в эквивалентную схему необхо димо ввести дополнительный элемент. Сила жидкостного трения пропорцио нальна координате^ и скорости скольжения:
|
F = -* -К |
(4.6) |
где |
i t , - коэффициент пропорциональности, |
|
|
У0 - скорость скольжения. |
|
Если рассматривать У0 как постоянную величину, то структура эквива лентной схемы остается без изменений. Изменится только коэффициент в урав нении связи:
JF = (kJ+kxVa) y . |
(4.7) |
Для учета изменений скорости скольжения при перемещении рабочего органа по направляющим в эквивалентную схему необходимо ввести нелинейный элемент жесткости, который учитывает жидкостное сопротивление движе нию. Эквивалентная схема при таком подходе показана на рис. 4.15.
Рис. 4.15. Эквивалентная схема привода рабочего органа в пространстве двух координат с учетом жидкостного и смешанного трения
Полная математическая модель динамической системы привода подач в пространстве двух координат, полученная по приведенным выше эквивалент ным схемам, представляет собой систему двух обыкновенных дифференциаль ных уравнений:
М х+ Нхдг+ Схх + куу =О,
(4.8)
М у+ hy у+ Суу + kQх = О,
где М - масса рабочего органа; К к у - коэффициенты демпфирования в связи по оси х и стыке по оси у\
Сх>Су - жесткости связи по оси х и стыка по оси у;
kj. - коэффициент пропорциональности между перемещением по оси у и силой трения;
kQ - коэффициент пропорциональности между скоростями скольжения и усилием всплывания.
Исходные данные для программы PAN, соответствующие эквивалентной схеме (см. рис. 4.14), приведены в табл. 4.4.
В табл. 4.5 представлены исходные данные по эквивалентной схеме, изо браженной на рис 4.15, которая описывает динамику привода рабочего органа в пространстве двух координат.
Таблица 4.4 Исходные данные привода рабочего органа с учетом силы трения
|
Ветвь |
|
Узел |
|
Значение |
|
№ |
тип |
от |
до вли |
C,R,L |
E,J,k |
знак |
|
|
|
яния |
|
|
|
1 |
L |
1 |
0 |
L = 1/ С, |
2 |
R |
1 |
0 |
II 0* |
3 |
С |
1 |
0 |
с = м |
4 |
PF |
1 |
0 |
2 |
5 |
JF |
2 |
0 |
1 |
6 |
R |
2 |
0 |
R = l//jT=7V |
7 |
L |
2 |
0 |
L = 1 |
8 |
J |
1 |
0 |
|
k= \
и
Таблица 4.5 Исходные данные привода рабочего органа в пространстве двух координат
|
Ветвь |
|
Узел |
|
|
Значение |
|
№ |
тип |
от |
до |
вли |
C,R,L |
E,J,k |
знак |
|
|
|
|
яния |
|
|
|
i |
L |
1 |
0 |
|
L = 1 ICX |
|
|
2 |
R |
1 |
0 |
|
II * |
|
|
|
|
|
|
||||
3 |
С |
1 |
0 |
|
с = м |
|
|
4 |
PF |
1 |
0 |
2 |
|
к= \ |
|
5 |
JF |
2 |
0 |
1 |
|
и |
|
|
|
|
|||||
6 |
С |
2 |
0 |
|
с = м |
|
|
7 |
L |
2 |
0 |
|
ги |
|
|
8 |
R |
2 |
0 |
|
R = \!hy |
|
|
9 |
J |
1 |
0 |
|
|
1 |
|
Процесс резания оказывает определяющее влияние на динамические ха рактеристики и виброустойчивость динамической системы станка.
Так как всем рабочим процессам присущи общие черты [17], то для моде лирования процесса резания методом прямой аналогии используются приемы, аналогичные моделированию процесса трения.
Математические модели для процесса резания достаточно разработаны [13,17,19,21]. Вид и сложность модели определяется решаемой задачей. Для описания процесса резания используются: статическая модель, апериодическое звено, дифференцирующее звено, полная линейная математическая модель и т.д.
Динамические характеристики процесса резания и связь с другими звень ями динамической системы станка определяются геометрическими параметра ми срезаемого слоя (рис. 4.16).
Рис. 4.16. Геометрические параметры срезаемого слоя: a.h - толщина и ширина срезае мого слоя; I- глубина резания; <р - главный угол в плане; е - изменение толщины срезае
мого слоя
Если пренебречь инерционностью и демпфированием процесса резания, то модель будет представлять собой статическую характеристику вида
|
Р = кг гч |
(4.9) |
где кР- |
коэффициент силы резания, |
|
е - |
изменение толщины срезаемого слоя. |
|