Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Моделирование технологического оборудования

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
4.93 Mб
Скачать

Гидросуппорт представляет собой следящую систему с однокромочным золотниковым распределителем. Основными узлами являются: гидростанция Г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)

где кР-

коэффициент силы резания,

 

е -

изменение толщины срезаемого слоя.

 

Соседние файлы в папке книги