1.1. Выбор гидроцилиндра.
По заданным параметрам выбираем гидроцилиндр ЦГ 140/100*1600
Основные параметры:
Диаметр поршня, мм 140;
Диаметр штока, мм 100;
Ход поршня, мм 1600;
Рабочее давление, МПА 25;
Расход рабочей жидкости, л/мин 120;
Нагрузка на гидроцилиндре, кг 1100;
Давление страгивания, МПа 0,25.
2. Выбор конструктивных элементов гидропривода.
2.1. Выбор насоса.
2.1.1. Расчет давления на выходе из насоса.
Рнрас = Ргц + ΔРзол + ΔРl, (2)
где: Ргц – рабочее давление в гидроцилиндре, Ргц=25 МПа (по заданию);
ΔРзол – потери давления в распределители;
ΔРl – потери давления в трубопроводе.
Определим диаметр трубопровода
dтр = (66,7·Qн / π·[Vж])1/2, (3)
где Vж – скорость течения жидкости в трубопроводе, Vж=2,5 м/с.
dтр = (66,7·120 / 3,14·2,5)1/2=31,9 мм
dтр = 32 мм
ΔРзол = 0,15 МПа [1]
ΔРl = 7,85·(L·Qн2 / d5) , (4)
где L – длина трубопровода, L = 10 м.
ΔРl = 7,85·(10·1202 / 325) = 0,03 МПа
Рн рас = 25 + 0,15 + 0,03 = 25,18 МПа
2.1.2. Расчет требуемой подачи насоса.
Qнрас = Qгц + ΔQгц + ΔQзол + ΔQпк, (5)
где: Qгц – рабочий расход в гидроцилиндре, Qгц = 20 л/мин;
ΔQгц – утечки в гидроцилиндре;
ΔQзол – утечки в распределителе;
ΔQпк – утечки в предохранительном клапане.
ΔQгц = ΔQ*гц ·Ргц / Р*гц, (6)
где ΔQ*гц – утечки в гидроцилиндре, рассчитанные для давления
Р*гц = 5 МПа, ΔQ*гц = 102 см3/мин.
ΔQгц = 102·25·106 / 5·106 = 510 см3/мин
ΔQзол = ΔQ*зол·Рзол / Р*зол, (7)
где ΔQ*зол – утечки в распределителе, рассчитанные для давления.
Р*зол = 5 МПа, ΔQ*зол = 300 см3/мин.
ΔQзол = 300·25·106 / 5·106 = 1500 см3/мин
ΔQпк = 0,1·Qн (8)
Qн рас = (120 + (510 + 1500)·10-3) / 0,9 = 134,5 л/мин
2.1.3. Выбор насоса по рассчитанным параметрам.
Исходя из условий Рн ≥ Рн рас
Qн ≥ Qн рас , выбираем насос.
Данным условиям соответствует аксиально-поршневой регулируемый насос 313.160.
Основные параметры:
Номинальный рабочий объем, см3 160;
Номинальная подача, л/мин 182;
Номинальное давление, МПа 20;
Номинальная частота вращения, об/мин 1200;
Максимальная частота вращения, об/мин 2650;
Минимальная частота вращения, об/мин 400;
Полный КПД, % 90.
2.2. Выбор электродвигателя.
2.2.1. Минимальная частота вращения вала электродвигателя:
nэд.min = nн.min = 400 об/мин.
2.2.2. Максимальная частота вращения вала электродвигателя:
nэд.max = nн.max = 2650 об/мин.
2.2.3. Расчет требуемой мощности электродвигателя:
Nэд рас = Nн / (ηм·ηн), (9)
где Nн – мощность насоса;
ηм – КПД муфты, ηм = 0,98;
ηн – КПД насоса, ηн = 0,9.
Nн = Рн рас ·Qн, (10)
Nн = 25,18·106·3,03·10-3 = 76 кВт
Nэд рас = 76·103 / (0,98·0,9) = 86 кВт
2.2.4. Выбор электродвигателя по рассчитанным параметрам.
Исходя из условий Nэд ≥ Nэд рас
nэд.min < nэд < nэд.max , выбираем электродвигатель.
Данным условиям соответствует асинхронный двигатель 4Ф280М6У3:
Мощность элегтродвигателя Nэд = 86 кВт;
Номинальная частота вращения вала электродвигателя nэд = 1000 об/мин.
2.3. Выбор предохранительного клапана.
Предохранительный клапан выбираем по номинальному расходу и давлению:
Ркл.ном ≥ Рн
Qкл.ном ≥ Qн - Qн рас
Данным условиям соответствует предохранительный клапан У462.8.5.0:
Номинальный расход Qкл.ном = 63 л/мин;
Номинальное давление на входе Ркл.ном = 20 МПа.
2.4. Выбор гидрораспределителя.
Распределитель выбираем по номинальному расходу и давлению:
Ррас.ном ≥ Рн
Qрас.ном ≥ Qн рас
Данным условиям соответствует распределитель типа МКРН.306.150.088:
Номинальный расход Qрас.ном = 135 л / мин;
Номинальное давление на входе Ррас.ном = 32 МПа;
Диаметр условного прохода Dу = 32 мм.
3. НЕЛИНЕЙНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ГИДРОПРИВОДА
Структурная схема нелинейной математической моделиданного гидропривода, построенная в программе MatLab-Simulink, приведена в приложении 1.
3.1. Нюансы построения модели.
В системе Simulink нежелательно использовать очень малые значения параметров блоков. Поэтому, для исключения ошибок счета, в некоторых частях схемы будем заменять единицы измерения (путем введения поправочных коэффициентов), т.е.:
60с = 1 мин;
1 м3/с = 6·107 см3/мин;
1кгс = 9,81 Н.
3.2. Описание звеньев нелинейной математической модели.
3.2.1. Входной ступенчатый сигнал.
Параметры входного ступенчатого сигнала в программе Simulink будут иметь следующие значения
Step time: 0;
Initial value: 0;
Final value: 1;
Sample time: 0.
3.2.2. Звено, описывающее распределитель.
Для данного распределителя с Dу = 32 мм:
Кр = 2;
Тр = 30 мс;
ζ = 0,8.
3.2.3. Коэффициент усиления К1.
На вход данного звена подается значение перемещения золотника распределителя, на входе имеем теоретический расход.
К1 = Кед.1·КQx, (11)
где: Кед.1 = 100 – коэффициент, учитывающий смену единиц измерения;
КQx = Qmax / xmax, (12)
где: Qmax – максимальный расход, проходящий через распределитель, Qmax = 135 л / мин;
xmax – максимальное перемещение золотника, xmax = 0,8 см.
КQx = 135·103 / 0,8 = 169·103 см2 / мин
К1 = 100·169·103 = 16,9*106 см3 / (мин·м)
3.2.4. Звено учитывающее нелинейность типа «насыщение по расходу».
На входе данного звена – теоретический расход, на выходе – реальный расход.
Параметры данного звена в программе Simulink будут иметь следующие значения:
- Upper limit: 13,45·104 (Qmax=Qн рас )
- Lower limit: -13,45·104 (Qmin=-Qн рас)
3.2.5. Звено, учитывающее сжимаемость рабочей жидкости.
На входе в данное звено – «расход сжимаемости», на выходе – давление в
поршневой полости гидроцилиндра.
Здесь Vц – объём гидроцилиндра
Vц = (dп2· /4) · h , (13)
где dп – диаметр поршня гидроцилиндра, см (см. задание)
h – ход поршня гидроцилиндра, см (см. задание)
Vц = (142·3,14)/4·160 = 24630 см3
- Е – модуль упругости, Е=14*10-5 см2 / кг
- Кqp – коэффициент утечек,
Кqp = Qзол / (рнрас – 0,5· рl ) (14)
Кqp = 1500 / (251,8 – 0,5·0,3) = 5,96 см5 / (кгс·мин)
3.2.6. Звено, учитывающее нелинейность типа «насыщение по давлению».
На вход в данное звено – теоретическое давление, на выходе – реальное давление.
Параметры данного звена в программе Simulink будут иметь следующие значения:
- limit: 250 (pmax = pгц)
- limit: -250 (pmin = -pгц)
3.2.7. Коэффициент усиления К2.
На вход данного звена подается значение реального давления в поршневой полости гидроцилиндра, на выходе имеем силу, создаваемую этим давлением.
К2 = Кед.2 · ·dп2/4, (15)
где Кед.2 = 9,81 – коэффициент, учитывающий смену единиц измерения
К2 = 9,81·3,14·142 / 4 = 1509,4 Н / кгс
3.2.8. Звено, описывающее механическую часть гидропривода.
Здесь m – приведенная масса на гидроцилиндре (см. задание);
fв.тр. – коэффициент вязкого трения.
f в.тр. = 0,036·lтруб · ·dп4/ dтруб3 , (16)
где lтруб – длина напорной линии, примем lтруб = 1 м;
dтруб – диаметр напорного трубопровода, примем dтруб = 0,032 м.
f в.тр = 0,036·1·3,14·0,144 / 0,0323 = 1,325 кг / с .
3.2.9. Звено, учитывающее нелинейность типа «сухое трение».
На входе в данное звено – теоретическая сила трения, на выходе – реальная сила трения. Максимальное значение силы сухого трения определяется следующим образом:
Рс.тр.max = рстраг · ·dп2/4, (17)
где рстраг – давление страгивания, рстраг = 0,25 МПа
Рс.тр.max = 0,25·106·3,14·0,142 / 4 = 3846,5 Н
Параметры данного звена в программе Simulink будут иметь следующие значения:
Upper limit: 3846,5 (Рс.тр.м = Р с.тр.maxрас)
Lower limit: -3846,5 (Рс.тр.м = -Р с.тр.maxрас)
3.2.10. Коэффициент усиления К4.
На вход данного звена подаётся значение скорости поршня гидроцилиндра, на выходе имеем силу сухого трения.
Для мгновенного срабатывания звена, описывающего сухое трение принимаем К4 = 1·106 кг/с.
3.2.11. Коэффициент усиления К3.
На вход данного звена подается значение скорости поршня гидроцилиндра, на выходе имеем геометрический расход в гидроцилиндре.
К3 = Кед.3 · · dп2 / 4, (18)
где: Кед.3 = 9,81 – коэффициент, учитывающий смену единиц измерения,
К3 = 9,81·3,14·142 / 4 = 1509,4 (см3·с) / (м·мин).
3.2.12. Интегрирующее звено.
На вход данного звена подается значение скорости поршня гидроцилиндра, на выходе имеем перемещение поршня гидроцилиндра.
3.2.13. Звено, ограничивающее перемещение поршня гидроцилиндра.
Параметры данного звена в программе Simulink будут иметь следующие значения:
Upper limit: 1,60 (xmax = h)
Lower limit: -1,60 (xmax = -h)
4. МОДЕЛИРОВАНИЕ ГИДРОПРИВОДА.
4.1. Переходный процесс по перемещению
Время полного переброса нагрузки: t переброса = 14,5с
4.2 Переходный процесс по скорости
Максимальная скорость поршня Vгц = 0,27 м/с
4.3. ЛАФЧХ.
Анализ частотных характеристик показал, что система обладает большими запасами устойчивости.