книги / Моделирование технологического оборудования
..pdfС учетом этого эквивалентная схема, отражающая динамику изгиба не весомого стержня постоянного сечения, показана на рис. 3.17.
|
|
■ A W |
|
|
|
|
|
|
С - 12 EJIl |
|
N |
|
|
ы D |
|
-Л Л Д г |
|
--- |
*TJ |
|
|
|
|
Ф; |
|
||
|
|
C = AEJI l |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
--- |
|
х |
PFyl Dr jМ |
Р |
kt - 6EJ/1 |
|
* № |
- |
|
|
•г |
|||||
|
к = 6EJI12 |
P F |
|
|
|
Рис. 3.17. Эквивалентная схема изгиба стержня
pF<i(j
Исходные данные для расчета изгиба стержня с помощью программы PAN имеют структуру, представленную в табл. 3.5.
Таблица 3.5
Исходные данные системы изгиба стержня
|
Ветвь |
|
Узел |
|
№ |
тип |
от |
ДО |
вли |
|
|
|
|
яния |
1 |
PF |
1 |
0 |
3 |
2 |
PF |
1 |
0 |
4 |
3 |
L |
1 |
2 |
|
4 |
PF |
2 |
0 |
3 |
5 |
PF |
2 |
0 |
4 |
6 |
PF |
3 |
0 |
1 |
7 |
PF |
3 |
0 |
2 |
8 |
PF |
3 |
0 |
4 |
9 |
L |
3 |
4 |
|
10 |
PF |
4 |
0 |
1 |
11 |
PF |
4 |
0 |
2 |
12 |
PF |
4 |
0 |
3 |
|
Значение |
C,R,L |
E,J,k |
к = 6£///j
к = 6EJ/I2
L = 1/(12EJIl1)
к = 6 EJli1 II a к =6EJtf к = 6 EJIl1 k\ = 6EJ/1
L = 1/(4EJIl)
к = 6 EJIl1 к = 6 EJli1 k\ = 6EJ/1
знак
-1 -1
-1
-1
4. МАТЕМАТИЧЕСКИЕ МОДЕЛИ И ПРИМЕРЫ МОДЕЛИРОВАНИЯ ОСНОВНЫХ УЗЛОВ И ПРОЦЕССОВ ТЕХНОЛОГИЧЕСКОГО ОБОРУДОВАНИЯ
4.1. Пример моделирования шпиндельного узла
Шпиндельный узел во многих случаях определяет динамическое качест во некоторых станков [19, 20, 26]. К таким станкам относятся токарные, коор динатно-расточные, расточные, внутришлифовальные и т.д. При моделирова нии шпиндельный узел рассматривают как ступенчатую балку на упругих опо рах с вязким демпфированием. Балку разбивают на стержни постоянного сече ния, разграниченные изменением диаметра, опорами, сосредоточенной массой, внешней сосредоточенной нагрузкой и т.д., которые имеют распределенные по длине свойства жесткости. Смонтированные на шпинделе шкивы, зубчатые ко леса, патрон, маховики представляются в виде сосредоточенных масс, распо ложенных на границах стержней. Опоры шпинделя отражают как линейную жесткость по соответствующим координатам с вязким демпфированием.
Для исследования реакций шпинделя на внешние воздействия в рамках линейной модели используют принцип суперпозиции.
Применение метода прямой аналогии для исследования динамических характеристик рассмотрим на примере шпинделя, изображенного на рис. 4.1.
Рис. 4.1. Принципиальная схема шпиндельного узла
Шпиндель разбит на три участка, каждый из которых отражается невесо мым стержнем постоянного сечения. На концах шпинделя расположены сосре доточенные массивы: шкив ременной передачи (узел 1) и патрон (узел 4); со стороны процесса резания действует сила F.
Правая опора шпинделя учитывается в модели пружиной, имеющей ра диальную и угловую жесткости (Су, Сф), и элементами демпфирования (Ry, /?ф).
Левая опора отражена пружиной, имеющей радиальную жесткость (Су), и элементом демпфирования (Ry).
Эквивалентная схема данного шпиндельного узла изображена на рис. 4.2.
Рис. 4.2. Эквивалентная схема шпиндельного узла
Для расчета динамики шпиндельного узла с помощью программы PAN структура исходных данных представлена в табл. 4.1.
Таблица 4.1
Исходные данные системы шпиндельного узла
|
Ветвь |
|
Узел |
№ |
тип |
от |
до |
1 |
С |
1 |
0 |
2 |
PF |
1 |
0 |
3 |
PF |
1 |
0 |
4 |
L |
1 |
2 |
5 |
R |
1 |
2 |
6 |
PF |
2 |
0 |
7 |
PF |
2 |
0 |
8 |
L |
2 |
0 |
9 |
R |
2 |
0 |
10 |
PF |
2 |
0 |
11 |
PF |
2 |
0 |
12 |
L |
2 |
3 |
13 |
R |
2 |
3 |
14 |
PF |
3 |
0 |
15 |
PF |
3 |
0 |
16 |
L |
3 |
0 |
17 |
R |
3 |
0 |
18 |
PF |
3 |
0 |
19 |
PF |
3 |
0 |
20 |
L |
3 |
4 |
21 |
R |
3 |
4 |
22 |
PF |
4 |
0 |
23 |
PF |
4 |
0 |
24 |
С |
4 |
0 |
25 |
С |
5 |
0 |
26 |
PF |
5 |
0 |
27 |
PF |
5 |
0 |
28 |
PF |
5 |
0 |
вли |
|
с д х |
|
Значение |
|
|
E J,k |
||||
|
|
|
|
|
|
||||||
яния |
|
|
|
|
|
|
|
|
|
|
|
О II £ |
|
|
|
|
* = 6EJH2 |
||||||
5 |
|
|
|
|
|
||||||
6 |
|
|
|
|
|
А |
|
= 6EJ/F |
|||
L |
= 1 / С |
, = |
1 /(ПЕЛ?) |
|
|
|
|
|
|||
«г II £ |
|
|
|
|
к |
= 6 £ |
///3 |
||||
5 |
|
|
|
|
|
||||||
6 |
|
|
|
|
|
А |
|
= 6 |
£ |
/ / / 2 |
|
г |
и |
3 |
|
|
|
|
|
|
|
|
|
04 |
и |
|
|
|
|
Т ^ б Ё Ш |
|||||
6 |
|
|
|
|
|
||||||
7 |
|
|
|
|
|
~кТШ Л ? |
|||||
L |
= МСУ= 1 |
/ ( 1 |
2 |
£ |
/ / / J ) |
|
|
|
|
|
|
е* |
II |
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
II |
|
|
1 |
|
7. |
|
|
|
|
|
А |
= 6 |
£ |
/ / / * |
||
L = 1/С/ |
|
|
|
|
|
|
|
|
|
||
R = Л / |
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
|
|
II |
|
|
1 |
8 |
|
|
|
|
|
к |
= 6EJ/12 |
||||
L |
= \/Су= 1 |
/ ( 1 |
2 |
E |
/ / / J ) |
|
|
|
|
|
|
II & |
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
к |
= 6EJ/12 |
||||
8 |
|
|
|
|
|
А |
= 6EJ/12 |
||||
II и |
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
А |
= |
6 |
£ |
/ / / 2 |
|
2 |
|
|
|
|
|
А |
= |
6 |
£ |
/ / / 2 |
|
6 |
|
|
|
|
|
.А: |
= 6EJH |
знак
-1 -1
-1 -1
-1 -1
-1
|
Ветвь |
|
Узел |
|
№ |
тип |
от |
ДО |
вли |
|
|
|
|
яния |
29 |
L |
5 |
6 |
|
30 |
R |
5 |
6 |
|
31 |
PF |
6 |
0 |
1 |
32 |
PF |
6 |
0 |
2 |
33 |
PF |
6 |
0 |
5 |
34 |
Pf |
6 |
0 |
2 |
35 |
PF |
6 |
0 |
3 |
36 |
PF |
6 |
0 |
7 |
37 |
L |
6 |
7 |
|
38 |
R |
6 |
7 |
|
39 |
PF |
7 |
0 |
2 |
40 |
PF |
7 |
0 |
3 |
41 |
PF |
7 |
0 |
6 |
42 |
L |
7 |
0 |
|
43 |
R |
7 |
0 |
|
44 |
PF |
7 |
0 |
3 |
45 |
PF |
7 |
0 |
4 |
46 |
PF |
7 |
0 |
8 |
47 |
L |
7 |
8 |
|
48 |
R |
7 |
8 |
|
49 |
PF |
8 |
0 |
3 |
50 |
PF |
8 |
0 |
4 |
51 |
PF |
8 |
0 |
7 |
52 |
C |
|
|
|
53 |
J |
4 |
0 |
|
Значение
C,R,L
L = 1/(4EJ/l)
R = ЯФ
L = l/(4£///) li 04
II и >5
L = 1/(4EJ/[)
IIP4
C = J 4
E,J,k
*= 6£///2
*= 6£///2
it = 6£///
к= 6EJ/l2
к= 6E///2 it = 6£///
к= 6£///-
*= 6E///1
к= 6£///
*= 6£///3
it = 6£///
it = 6£///~ II vO к = 6£///
F
_| знак ,
-1
-1
-1
-1
-1
-1
Инерционность материала шпинделя при необходимости может быть уч тена путем присоединения к концам стержней сосредоточенных масс по из вестной методике [27].
4.2. Пример моделирования механического привода станка
Привод вместе с электродвигателем рассматривается как цепная система [8,9,26,27]. Он представляет собой крутильную, упругую механическую систему, состоящую из последовательно соединенных посредством упругих и демпфирующих элементов (пружин и сопротивлений) сосредоточенных масс (массивов). Значение податливости упругих элементов определяется деформациями закручивания и изгиба валов, контактными деформация!ии в соединениях валов с колесами, зубчатых передачах, ременных передачах и т.д. [27]. Сопротивления отражают демпфирование в этих же элементах кинематических цепей привода и обычно включаются в расчетную цепь параллельно с упругими элементами.
Использование метода прямой аналогии для моделирования динамики механического привода рассмотрим на примере главного привода станка 6Н13ПБ (рис. 4.3) [27]. При этом применяются приемы моделирования рычажной системы, изложенные выше.
Рис. 4.3. Кинематическая схема привода станка 6Н13ПБ
Метод прямой аналогии не требует приведения системы к тому или иному валу. Механизм приведения скоростей органически заложен в схему составления расчетной схемы. Для этого достаточно разделить весь привод на подсистемы, которые включали бы в себя элементы с одинаковыми скоростями. А при установлении связей между подсистемами учитывается передаточное отношение кинематических пар. Исходя из этого, в приводе станка 6Н13ПБ целесообразно выделить 7 подсистем, взаимодействие между которыми отражено в механической цепи привода (рис. 4.4) элементами к,.
£|СЭ*И.З |
^ IВ+Ш+Ш |
L\B*U1+M |
^КЗ-ИЗ />4в+ш+ш |
Рис. 4.4. Механическая цепь привода станка 6Н13ПБ при л = 60 об/мин
Такой подход дает дополнительное удобство, так как при анализе математической модели в результате расчетов получаются действительные координаты, скорости и усилия в элементах исходной системы.
Элементами механической цепи отражены податливости и сосредоточен ные массы. При необходимости к ним могут быть добавлены элементы сопроти вления. На рис. 4.4 обозначены:
Ljв, 1 ш, LM, LHз, Z,K.з - податливости кручения соответствующих валов, приведенные к крутильной податливости контактные деформации шлицевых (шпоночных) соединений, муфт, изгибные и контактные деформации зубчатых передач;
Jp, Л,, */Шп - моменты инерции ротора электродвигателя, полумуфт, шпинделя;
J:ij - моменты инерции зубчатых колес, шкивов и т.п.
По механической цепи привода на основе второй системы аналогий получаем эквивалентную схему (рис.4.5). Фиктивные источники в эквивалентной схеме реализуют трансформаторный тип связи и отражают взаимодействие между подсистемами, учитывая передаточные отношения соответствующих передач.
Демпфирование в зубчатых зацеплениях мало [27], и им можно пренебречь. Демпфирование в других элементах привода, прежде всего в стыках соединений и материалов валов, учитываем в эквивалентной схеме элементами сопротивлений, включенных параллельно пружинам, отображающим жесткость соответствующих элементов исходной схемы.
Рис. 4.5. Эквивалентная схема привода станка 6Н13ПБ при п = 60 об/мин
Структура исходных данных для исследований динамики привода с помощью программы PAN приведена в табл. 4.2.
|
|
|
Исходные данные системы привода станка 6Н13ПБ |
Таблица 4.2 |
|||
|
|
|
|
||||
|
|
Ветвь |
|
Узел |
|
Значение |
знак |
[ № |
тип |
от |
Д О вли |
C,R,L |
E,J.k |
||
|
|
|
|
яния |
|
|
|
|
|
L |
1 |
0 |
и |
|
|
к |
3 |
1 R |
1 |
0 |
R = 1/бя |
|
|
| |
С |
1 |
0 |
с = / , |
|
|
|
Г |
4 |
L |
1 |
2 |
L = £пми+ш |
|
|
|----- |
R |
1 |
2 |
R = l/^uMii+ui |
|
|
|
|
5 |
|
|
||||
и |
|
С |
2 |
0 |
С = * /M + Z 2 7 |
|
|
Ветвь
№тип
7JLF
8ELF
9С
10L
11R
12С
13JLF
14ELF
15С
16JLF
17ELF
18С
19L
20R
21С
22JLF
23ELF
24С
25L
26R
27С
28JLF
29ELF
30С
31L
32R
133 С
34JLF
35ELF
36С
37L
38R
39С
40J
|
|
|
|
|
Окончание табл. 4.2 |
|
|
Узел |
|
|
Значение |
|
|
от |
до |
вли |
|
C,RJL |
EJ,k |
знак |
|
|
яния |
|
|
|
|
2 |
0 |
3 |
L = £к.з+и.з |
-it = -27/53 |
|
|
3 |
0 |
2 |
L = £кз+из |
-it = -27/53 |
|
|
3 |
0 |
|
С = Jz53 |
|
|
|
3 |
4 |
|
L = ^2В+Ш+Ш |
|
|
|
3 |
4 |
|
R = 1/62В+Ш+Ш |
|
|
|
4 |
0 |
|
О 1 |
N4 |
|
|
4 |
0 |
5 |
L = £к.з+и.з |
-А = -16/39 |
|
|
5 |
0 |
4 |
L = £к.з+и.з |
-it = -16/39 |
|
|
5 |
0 |
|
С = «/z39 |
-Л—-18/47 |
|
|
5 |
0 |
6 |
L = £кз+из |
|
||
6 |
0 |
5 |
L = ^к.з+и.з |
-it = -18/47 |
|
|
6 |
0 |
|
С = «/z47 |
|
|
|
6 |
7 |
|
L = Т#4в+ш+ш |
|
|
|
6 |
7 |
|
R —1/64В+Ш+Ш |
|
|
|
7 |
0 |
|
II и |
о |
-it =-20/75 |
|
7 |
0 |
8 |
L - |
£кз+из |
|
|
8 |
0 |
7 |
L = £кз+из |
-it = -20/75 |
|
|
8 |
0 |
- |
II CJ |
о |
|
|
8 |
9 |
|
L = ^5В+Ш+Ш |
|
|
|
8 |
9 |
|
R = 1/^5В+Ш+Ш |
|
|
|
9 |
0 |
|
С = Jz42 |
-it = -42/30 |
|
|
9 |
0 |
10 |
L = £кз+из |
|
||
10 |
0 |
9 |
L = Ькз+из |
-it = -42/30 |
|
|
10 |
0 |
|
II и |
о |
|
|
10 |
11 |
|
L - |
^бв+ш+ш |
|
|
10 |
11 |
|
R = Rбв+ш+ш |
|
|
|
11 |
0 |
|
II и |
о |
|
|
11 |
0 |
12 |
L - LKз+из |
-it = -80/60 |
|
|
12 |
0 |
11 |
L = ^кз+из |
-it = -80/60 |
|
|
12 |
0 |
|
п II |
Г |
|
|
12 |
13 |
|
L = Т/7В+Ш |
|
|
|
12 |
13 |
|
R = /?7В+Ш |
|
i |
|
13 |
0 |
|
С = */шп |
J = M |
|
|
13 |
0 |
|
|
|
-1 |
4.3. Пример моделирования гидропривода станка
Гидравлический привод широко используется в металлорежущих станках для создания движений формообразования, автоматизации и обеспечения работы некоторых узлов станков (гидравлических опор, устройств гидроразгрузки направляющих и т.п.). Метод прямой аналогии является эффективным способом получения математической модели таких гидрофицированных узлов станка [32]. Использование метода прямой аналогии рассмотрим на примере моделирования гидросуппорта ГСП-41, схема которого приведена на рис. 4.6.
Рис. 4.6. Гидравлическая схема гидросуппорта ГСП - 41