Механизм насоса с качающейся кулисой
СОДЕРЖАНИЕ: Принципы работы, механизм и назначение насоса с качающейся кулисой. Структурный анализ и силовой расчет рычажного механизма. Особенности выполнения геометрического расчета зубчатой передачи. Синтез кулачкового механизма, порядок построения его профиля.Содержание
Введение
1 Синтез и анализ рычажного механизма
1.1 Структурный анализ механизма
1.2 Определение недостающих размеров
1.3 Определение скоростей точек механизма
1.4 Определение ускорений точек механизма
1.5 Диаграммы движения выходного звена
1.6 Определение угловых ускорений и скоростей
1.7 Определение ускорений центров масс звеньев механизма
1.8 Аналитический метод расчёта механизма
2 Силовой расчет рычажного механизма
2.1 Определение сил инерции
2.2 Расчет диады 4-5
2.3 Расчет диады 2-3
2.4 Расчет кривошипа
2.5 Определение уравновешивающей силы
2.6 Определение мощностей
2.7 Определение кинетической энергии и приведенного момента инерции механизма
3 Геометрический расчет зубчатой передачи. Проектирование планетарного механизма
3.1 Геометрический расчет зубчатой передачи
3.2 Определение передаточного отношения планетарной ступени и подбор чисел зубьев колес
3.3 Определение частот вращения зубчатых колес
4 Синтез и анализ кулачкового механизма
4.1 Диаграммы движения толкателя
4.2 Определение минимального радиуса кулачка
4.3 Построение профиля кулачка
Список использованных источников
Введение
Механизм насоса с качающейся кулисой применяется в нефтеперерабатывающей промышленности и предназначен для откачки жидкости нефтяных скважин.
Подача жидкости регулируется автоматически за счёт кулачкового механизма.
Поршень получает возвратно-поступательное движение в цилиндре от электродвигателя через планетарный редуктор и рычажный механизм О1АО2С.
При движении поршня вверх осуществляется рабочий ход, при движении поршня вниз - холостой.
При рабочем ходе на поршень 5 действует сила полезного сопротивления.
Механизм насоса с качающейся кулисой - одностороннего действия.
Кулачок 6 получает вращение посредством зубчатой передачи z5-z6.
1. Синтез и анализ рычажного механизма
Схема механизма:
Исходные данные:
Q=3450 H
H=240 мм
m3=42 кг
K=1,6
m5=35 кг
nкр=150 об/мин
O1O2=625 мм
nдв=1500 об/мин
1.1 Структурный анализ механизма
Степень подвижности механизма определим по формуле:
W=3n-2p1-p2;
Где n- число подвижных звеньев,
P1- число одноподвижных кинематических пар,
P2- число двуподвижных кинематических пар,
W=35-27-0=1
Разложение механизма на группы Ассура:
Формула строения механизма: I(0,1)II(2,3)II(4,5)
Механизм II класса, 2 порядка
1.2 Определение недостающих размеров
Угол размаха кулисы
b=180
= 180(1,6–1)/(1,6+1)=41,53
Длину кривошипа определяем по формуле:
lO1A=0102 sin
=0,625sin20,76=0,22м,
Длину lO2В определим по следующей формуле:
lO2В=
=0,24/20,354=0.338 м,
Выберем масштабный коэффициент
Строим 12 планов положений механизма, приняв за начало отсчета крайнее положение, соответствующее началу рабочего хода механизма.
1.3 Определение скоростей точек механизма
Определяем u точки А кривошипа:
uA=w1lO1A
где w1-угловая скорость вращения механизма, определяется по формуле:
w1=
=
рад/с,
Скорость точки А определим по формуле:
uA=1 O1A=15,710,22=3,46 м/с,
План скоростей строим в масштабе:
ku=
=3,46/69,2=0.05 мс-1/мм
Скорость точки A’ находим графически, решая совместно систему:
uA’= ku РuA’
По свойству подобия определяем скорость точки C’:
РVc’=13633/280=16мм
Абсолютное значение скорости точки
C’= ku puC’=0.05*17=0,85м/с
Составим систему уравнений скоростей для нахождения u в точке C:
uс’=uc’+ucc’
uc’=uc’c+uc’c
На плане puC=19мм. Абсолютное значение u в точке C:
uC= ku puC=0.05*19=0,95м/с.
Для остальных 11 положений скорости определяются аналогично, их значения приведены в таблице 1.1
Таблица 1.1 - Значения скоростей
| Скорости, м/с | Положение механизма | |||||||||||
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
| uА | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 | 3,46 |
| uА’ | 1,65 | 2,95 | 3,4 | 3,4 | 3 | 2,15 | 0,7 | 1 | 2,9 | 3,15 | 1,7 | 0 |
| uC | 0,95 | 1,2 | 1,23 | 1,24 | 1,18 | 0,96 | 0,38 | 0,7 | 2,23 | 2,48 | 1,18 | 0 |
| uC’ | 0,85 | 1,24 | 1,29 | 1,26 | 1,18 | 0,96 | 0,37 | 0,65 | 2,14 | 2,42 | 1,16 | 0 |
1.4 Определение ускорений точек механизма
Ускорение точки А направлено по кривошипу к центру вращения О1.
aA=w12lO1A=15,7120.22=54,3м/с2
Выбираем масштабный коэффициент ускорений:
ka=
=54,3/108,6=0.25 м*с-2/мм
На плане ускорений изображаем ускорение точки А отрезком 
а=108 мм. Ускорение точки А’ определяем, решая совместно систему
Для этого используем
aA’Ak=2w3 VA’A=2VA’O2VA’A/A’O2=2Pva’aa’c/A’O2
c- коэффициент перечета определяется по формуле:
с=
=0.052/(0.50.0025)=2
Ускорение точки С’ определим из соотношения:
Ускорение точки С найдём, решая совместно систему
где
аCC’k=23VCC’=2 РuA’СС’c/А’О2;
Ускорения всех точек механизма найдены. Ускорения для остальных положений механизма определяются аналогично. Значения ускорений сводим в таблицу 1.2.
Таблица 1.2- Значения ускорений.
| Ускорения, м/с2 | Положение механизма | ||||||
| 1 | 3 | 5 | 7 | 9 | 11 | 0 | |
| aA | 54,4 | 54,4 | 54,4 | 54,4 | 54,4 | 54,4 | 54,4 |
| aA’Ak | 14,37 | 7,74 | 12,02 | 7,55 | 24,7 | 23,56 | 0 |
| aA’O2n | 3,88 | 13,84 | 10,74 | 0,775 | 20,24 | 7,9 | 0 |
| ac’ | 15,3 | 5 | 7 | 24 | 40,6 | 48,1 | 32 |
| aCC’k | 0,75 | 0,5 | 1,25 | 0,00 | 3,45 | 1,9 | 0 |
| aC | 16,45 | 6,25 | 8,75 | 26,5 | 41,8 | 50,35 | 34 |
1.5 Диаграмма движения входного звена
Диаграмму перемещения S-t строим используя полученную из плана механизма траекторию движения точки С. Диаграммы скоростей V-t и ускорений а-t определяются из полученных 12-ти планов скоростей и планов ускорений.
Масштабные коэффициенты диаграмм
Кs=0.0025 (м/мм)
КV=0.05 (мс-1/мм)
Ка=0,5 (мс-2/мм)
1.6 Определение угловых скоростей и ускорений
1==3,14150/30=15,71 (рад/с)
3=
=3,46/2800,0025=4,94 (рад/с)
1=0 (рад/с-2)
3=
=14,88/2800,0025=21,26 (рад/с-2)
1.7 Определение ускорений центров масс звеньев механизма
Ускорения центров масс звеньев механизма определяем из планов ускорений.
aS=KaРаS3=0.5*45=22,5 м/с2
Определяем относительные угловые скорости.
10= 1=15,71 рад/с;
30= 3=4,94 рад/с;
12= 1- 3=10,77 рад/с;
45= 3=4,94 рад/с;
1.8 Аналитический метод расчёта механизма
Исходные данные:
= 625 мм;
= 15,71 рад/с;
= 220 мм;
= 41,530;
= 308 мм;
= 20,760;
Расчет ведется для первого положения кулисы:
;
В проекциях на координатные оси:
;
Разделим второе уравнение ED Equation.3
;
Берем производную от левой и правой части:
Угловая скорость кулисы:
Угловое ускорение кулисы:
Составим векторное уравнение:
Проектируем на оси координат:
;
;
;
;
;
Расчет скоростей и ускорений для первого положения механизма.
Угловая скорость кулисы:
Угловое ускорение кулисы:
;
рад/с2 ;
Скорость точки С :
;
;
м/с;
Ускорение точки C :
;
рад/с2;
2 Силовой анализ рычажного механизма
2.1 Определение сил инерции
Исходные данные:
w1=15,71 рад/с;
Q=3450 Н;
m5=35 кг;
m3’=12 кг;
m3’’=30 кг;
Определим силы инерции:
U5=-m5aC;
U5= m5PaCKa;
U5=3530,60,5=535,5 (Н);
U3’=-m3’aS3’;
U3’=m3’PaS’Ka;
U3’=12450,5=270 (Н);
U3’’=-m3’’aS3’’;
U3’’=m3’’PaS’’Ka;
U3’’=30170,5=255 (Н);
Определим веса звеньев:
G5=m5g;
G3’=m3’g;
G3’’=m3’’g;
G5=359,8=343,35 (Н);
G3’=129,8=117,72 (Н);
G3’’=309,8=294,3 (Н);
Сила полезного сопротивления Q=3450 Н.
Разбиваем механизм на группы Ассура в соответствии с формулой строения I(0,1)II(2,3)II (4,5). Начинаем силовой рассчёт самой удалённой от кривошипа диады.
2.2 Расчёт диады II (4,5)
Выделим из механизма диаду 4-5 и нагружаем её силами. Составляем уравнение равновесия диады 4-5:
Р(4,5)=, R50+Q+G5+U5+R43=0 (1)
Уравнение содержит два неизвестных- модули реакций R50 и R43, поэтому оно решается графически. Строим план сил по уравнению равновесия (1).
Для построения плана сил выбираем масштаб сил Кр
Кр=
=3450/172,5=20 н/мм
Из плана сил определяем реакции:
R50= R50 Кр=6620=1320 Н;
R43= R43 Кр=22120=4420 Н;
2.3 Расчёт диады II (2,3)
Выделим диаду 2-3 и нагрузим её силами. Действие отброшенных звеньев 1,0 на третье заменяем действием реакций связей R21 и R30, которые требуется определить. Реакцию R21 направляем перпендикулярно линии движения ползуна, модуль неизвестен. Реакция R30 в шарнире О2 неизвестна ни по модулю ни по направлению; на схеме направляем её произвольно. Действие отброшенного звена 4 на третье известно: Реакция R34 равна по величине и противоположно направлена реакции R43, которая уже определена из плана сил диады II (4,5). Силы тяжести G3’ и G3’’ наносим на диаду в центрах масс стержней S3’ и S3’’. Силы инерции U3’ и U3’’ прикладываем в точках К’ и К’’, расположенных на расстоянии 2/3 длин стержней. Силы инерции направляем противоположно ускорениям центров масс согласно плана ускорений.
Составляем условия равновесия диады II(2,3):
Р(2,3)=0, R21+G3’+U3’+G3’’+U3’’+R34+R30=0 (2)
Данное уравнение содержит три неизвестных: модуль реакции R21, модуль и направление реакции R30. Значит уравнение (2) графически не решается. Реакция R21 может быть определена аналитически из уравнения моментов сил относительно точки О2.
М О2 (зв.2,3)=0, R21AO2-U3’hu3’+G3’hg3’-U3’’hu3’’-G3hg3’’-R34O2C=0;
Откуда
R21= (U3’hu3’- G3’hg3’+ U3’’hu3’’+ G3hg3’’+ R34O2C)/ AO2
R21=(270233-117,7253+255102-294,374+4500132)/280=2539 Н
Теперь уравнение (2) содержит два неизвестных, а следовательно решается графически.
Строим план сил диады II(2,3) по уравнению (2). Считаем отрезки плана сил:
= U3’/Кр=270/20=13,5 мм.
= U3’’/ Кр=255/20=12,75 мм.
= R21/ Кр=2539/20=126,95 мм.
= G3’/ Кр=117,72/20=5,8 мм.
= G3/ Кр=294,3/20=14,7 мм.
Согласно уравнению (2) строим сумму векторов сил, откуда находим:
R30= 
Кр=27420=5480 Н.
2.4 Расчёт кривошипа
Силовой расчёт кривошипа состоит в определении реакции стойки на кривошип R10 и уравновешивающей силы Ру, имитирующей действие силы со стороны двигателя.
Реакция R21 известна, так как R12= R21. Величина Рур определиться из уравнения моментов сил относительно точки О1 кривошипа.
М О1 (зв.1)=0, РурАО1-R12hR12=0
Рур’= R12hR12/ АО1=2539 40/88=1154 Н
Реакция стойки на кривошип R10 определиться из условия равновесия кривошипа:
P(кр)=R21+Py+R10=0 (3)
По уравнению (3) строим план сил кривошипа, откуда определяем искомую реакцию R10
R10= R10Кр=11020=2200 Н.
2.5 Определение уравновешивающей силы методом Жуковского
Уравновешивающую силу можно определить с помощью план скоростей по методу рычага Жуковского.
Строим повёрнутый на 90план скоростей и приложим к нему все внешние силы, действующие на механизм. План скоростей рассматриваем как жёсткий рычаг с опорой в полюсе. Рычаг находится в равновесии под действием приложенных сил.
Составляем уравнение равновесия рычага в форме суммы моментов сил в форме суммы моментов сил относительно полюса плана скоростей.
МPv1=0
Pyp’Pva-(Q+U5+G5)PvC-U3’’hU3’’-G3’’hG3’’-U3’hU3’+G3’hG3’=0
Pyp’=((Q+U5+G5)PvC+ U3’’hU3’’+ G3’’hG3’’+ U3’hU3’- G3’hG3’)/ Pva
Pyp’=((3450+535,5+343,35)47+25533+294,324+27069-117,7213)/179
Pyp’=1173 Н
Сравниваем значения Pyp и Pyp’, найденные двумя способами
=( Pyp’- Pyp)/ Pyp’
=(1173-1154)100%/1173=1,62%
2.6 Определение мощности
Мгновенная потребная мощность привода насоса без учёта потерь мощности на трение определяется соотношением:
Npy=PypVA=11733,46=4058,58 Вт
Мощность привода, затраченная на преодоления только полезной нагрузки:
NQ=QVc=34500,95=3277,5 Вт
Потери мощности во вращательных кинематических парах:
N10=R10f’(1-0)rц=22000,13215,710,025=114,5 Вт
N12=R12f’(1-3)rц=25390,13210,770,025=90,2 Вт
N30=R30f’(3-0)rц=54800,1324,940,025=89,3 Вт
N45=R45f’(3-5)rц=44200,1324,940,025=72,05 Вт
Где rц-радиус цапфы вала, rц=0,025 м,
f’- приведенный коэффициент трения, f’=(1,2…1,5)f=0,132
Потери мощности в поступательных кинематических парах:
N23=R23f’VA’A=25390,1321,65=553 Вт
N34=R34f’VC’C=44200,1320,85=495 Вт
N50=R50f’VC=13200,1320,95=165,5 Вт
Суммарная мощность трения:
Nтр=Ni=N10+N12+N30+N45+N23+N34+N50
Nтр=114,5+90,2+89,3+72,05+553+495+165,5=1579,2 Вт
Мгновенная потребляемая мощность двигателя:
N=NРу+Nтр
N=4058,58+1579,2=5637,78 Вт
2.7 Определение кинетической энергии и приведенного момента инерции механизма
Кинетическая энергия механизма равна сумме кинетической энергии звеньев:
Тмех=Тi
Для механизма насоса с заданными параметрами кинетическая энергия звена равна:
Тi=Т3+Т5=
Где
JO2’=
=120,352/3=0,49 кгм2
JO2’’=
=300,1552/3=0,24 кгм2
Т3=(0,49+0,24)4,942/2=8,9 Дж
Т5=350,95/2=16,62 Дж
Тмех=8,9+16,62=25,52 Дж
За звено приведения обычно выбирают ведущее звено. Так как у исследуемого механизма ведущим звеном является кривошип, то кинетическая энергия определится по формуле:
Tпр=
Откуда находим приведенный момент инерции:
Jпр=
Jпр=225,52/15,712=0,2 кгм2
3 Геометрический расчет зубчатой передачи. Проектирование планетарного механизма
3.1 Геометрический расчет зубчатой передачи
Исходные данные:
Число зубьев шестерни Z5=11;
Число зубьев колеса Z6=25;
Модуль m=6 мм;
Нарезание проводится методом обкатки инструментом реечного типа, который профилируется на основе исходного контура по ГОСТ 13755-81 и имеет следующие значения: угол профиля 
; коэффициент высоты головки 
; коэффициент радиального зазора 
;
Определяем геометрические параметры эвольвентной передаче.
Определяем минимальный коэффициент смещения:
Z517и Z5+Z634, следовательно, передача равносмещенная,
x5=(17-Z5)/17=(17-11)/17=0,35 мм;
x6=-x5=-0,35 мм;
Определяем делительное межосевое расстояние:
а= 0,5m(Z5+Z6)= 0.56(11+25)=108 мм;
Определяем высоту зуба:
h=m(2ha*+c*)=6(21+0,25)=13,5 мм;
4) Делительная высота головки зуба:
ha=m(ha*+x);
ha5= m(ha*+x5)= 6(1+0,35)= 8,1 мм;
ha6=m(ha*+x6)=6(1-0,35)= 3,9 мм;
5) Делительная высота ножки зуба:
hf= m(ha*+C-x);
hf5= m(ha*+C-x5)= 6(1+0,25-0,35)= 5,4 мм;
hf6= m(ha*+C+x6)= 6(1+0,25+0,35)= 9,6 мм;
Диаметр делительной окружности:
d5= mZ5= 611= 66 мм;
d6= mZ6= 625= 150 мм;
Диаметр основной окружности:
db5= mZ5cos()= 611 cos(20)= 62,05 мм;
db6= mZ6cos()= 625 cos(20)= 147 мм;
Диаметр окружности вершин зубьев:
da=mZ+2m(ha*+x);
da5=mZ5+2m(ha*+x5)=611+26(1+0,35)= 82,2 мм;
da6=mZ6+2m(ha*+x6)=625+26(1-0,35)= 157,8 мм;
Диаметр окружностей впадин зубьев:
df=mZ-2m(ha*+C*-x);
df5=mZ5-2m(ha*+C*-x5)=115-26(1+0,25-0,35)=55,2 мм;
df6=mZ6-2m(ha*+C*-x6)=255-26(1+0,25+0,35)=130,8 мм;
10) Делительная окружная толщина зуба:
S=0,5m+2mxtg();
S5=0,53,146+260,35tg(20)= 10,9 мм;
S6=0,53,146-260,35tg(20)= 7,9 мм;
Делительный шаг:
P= m = 3,146=18,84 мм;
12) Основной шаг:
Pb= m cos()= 3,1460,94=17,7 мм;
13) Радиус кривизны галтели
=0,38m=2.28 мм;
14) Строим зубчатую передачу с масштабным коэффициентом Kl=0,00025 м/мм;
15) Проверяем коэффициент торцевого перекрытия
а) аналитический метод:
1,57
б) графический метод:
где 
– длина активной линии зацепления.
3.2 Определение передаточного отношения планетарной ступени и подбор числа зубьев колес
Исходные данные:
nкр=150 мин-1;
nдв=1500 мин-1;
Z5=11;
Z6=25;
знак передаточного отношения привода (-)
Составляем общее передаточное отношение механизма:
Рассчитаем передаточное отношение 
и 
через исходные данные:
Из исходного уравнения определяем передаточное отношение планетарной ступени:
;
Составляем формулу Виллиса для планетарной передачи:
;
;
Запишем через числа зубьев передаточное отношение обращенного механизма:
;
Подбираем числа зубьев:
; 
;
Z1+Z2=Z4-Z3;
Z1+Z2=30+30=60
Z3+Z4=85-25=60
Z1=30, Z2=30, Z3=25, Z4=85
По выбранным числам зубьев определяем размеры колес:
d=mZ;
d1=640=240 мм;
d2=640=240 мм;
d3=625=150 мм;
d4=685=510 мм;
d5=611=66 мм;
d6=625=150 мм
Масштабный коэффициент построения Кl=0,001 м/мм;
Для построения плана скоростей редуктора определяем скорость точки А:
м/с;
Строим план скоростей. Масштабный коэффициент плана скоростей
мс-1/мм;
3.3 Определение частот вращения зубчатых колес аналитическим методом
n1= nдв=1500 мин-1;
n6= nкр=150 мин-1;
;
мин-1;
;
мин-1;
мин-1;
Значения частот вращения получим графическим методом:
мин-1;
мин-1;
мин-1;
мин-1;
4 Синтез и анализ кулачкового механизма
4.1 Диаграммы движения толкателя
Исходные данные:
Максимальный подъём толкателя h=29 мм;
Фазовый рабочий угол =290;
Дезаксиал е=0 мм;
nкр=150 об/мин;
Z5=11;
Z6=25
Угол давления =25;
По заданному графику V-t графическим диференцированием получим график а-t, графическим интегрированием - S-t. Базы Н1=20 мм, Н2=25 мм. Методом исключения общего параметра t получим график V-S, a-S, a-V. Масштабные коэффициенты графиков:
Ks=
м/мм;
Kv=
мс-1/мм
Kt=
c/мм;
Ka=
мс-2/мм
4.2 Определение минимального радиуса кулачка
Минимальный радиус кулачка выбирается из условия выполнения угла давления. Для этого строим совмещённый график S’-V, где S’- текущее перемещение в стандартном масштабе КS’=0,0005 м/мм, V- аналог скорости.
На совмещённом графике на горизонтальных линиях откладываем аналоги скорости в масштабе КS’
x1=
мм
x2=
К совмещённому графику проводим две касательные под углом давления . Ниже точки пересечения касательных выбирается центр вращения кулачка и соединяется с началом совмещённого графика. Это и будет минимальный радиус кулачка.
R0’=R0’KS’=400,0005=0,02 м;
4.3 Построение профиля кулачка
Профилирование кулачка выполняется методом обращённого движения. Для этого строим кулачок в масштабе Кl=0,00025 м/мм. Проводим окружность радиусом R0’ и окружность радиуса е. Откладываем угол р=290. Делим его на 12 частей и через точки деления проводим оси толкателя в обращённом движении. Вдоль осей толкателя откладываем текущее перемещение толкателя от окружности R0’. Соединяя полученные точки получим центровой профиль кулачка. Радиус ролика выбираем из условия:
rp=(0,2…0,4)R0’=0,2540=10 мм
Минимальный радиус действительного профиля:
R0=R0’-rp=40-10=30 мм
Обкатывая ролик по центровому профилю получаем действительный профиль.
Public Sub kul()
Dim I As Integer
Dim dis1, dis2, R, a1, a2, arksin1, arksin2, BETTA, BET As Single
Dim R0, FIR, FI0, FII, SHAG, E As Single
Dim S(1 To 10) As Single
R0 = InputBox(ВВЕДИТЕ МИНИМАЛЬНЫЙ РАДИУС КУЛАЧКА RO)
FIR = InputBox(ВВЕДИТЕ РАБОЧИЙ УГОЛ КУЛАЧКА FIR)
FI0 = InputBox(ВВЕДИТЕ НАЧАЛЬНОЕ ЗНАЧЕНИЕ УГЛА ПОВОРОТА КУЛАЧКА FI0)
E = InputBox(ВВЕДИТЕ ДЕЗАКСИАЛ E)
For I = 1 To 10
S(I) = InputBox(ВВЕДИТЕ СТРОКУ ПЕРЕМЕЩЕНИЙ S( I ))
Next I
FIR = FIR * 0.0174532
SHAG = FIR / 10
FI0 = FI0 * 0.0174532
FII = FI0
For I = 1 To 10
dis1 = (R0 ^ 2 - E ^ 2) ^ (1 / 2)
dis2 = S(I) ^ 2 + R0 ^ 2 + 2 * S(I) * dis1
R = dis2 ^ (1 / 2)
a1 = E / R
a2 = E / R0
arksin1 = Atn(a1 / (1 - a1 ^ 2) ^ (1 / 2))
arksin2 = Atn(a1 / (1 - a2 ^ 2) ^ (1 / 2))
BETTA = FII + arksin1 - arksin2
BETTA = BETTA * 180 / 3.1415
Worksheets(1).Cells(I, 1) = R
Worksheets(1).Cells(I, 2) = BETTA
FII = FII + SHAG
Next I
End Sub
Список использованных источников
1. Машков А.А. Теория механизмов и машин. Мн., 1971.
2. Артоболевский И.И. Теория механизмов и машин. М., 1975.
3. Фролов К.В., Попов С.А., Мусатов А.К. и др. Теория механизмов и машин под ред. К.В. Фролова М., 1986.
4. Попов С.А., Тимофеев Г.А. Курсовое проектирование по теории механизмов и механике машин. М., 1998.