CC BY
86
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Оптимизация силовой части электрогидравлического
следящего привода с дроссельным регулированием
# 11, ноябрь 2013
Б01: 10.7463/1113.0623966
Замараев Д. С., Попов Д. Н.
УДК 62-52
Россия, МГТУ им. Н.Э. Баумана popov@bmstu.ru d.stalky@gmail.com
Введение
Согласно концепции оптимизации электрогидравлического следящего привода (ЭГСП) с дроссельным регулированием, изложенной в [1], в данной статье рассмотрен первый уровень оптимизации (Ур. 1) ЭГСП, исходя из предположения, что примененный в приводе электрогидравлический усилитель (ЭГУ) представлен пропорциональным звеном. Выполнению этого условия при втором уровне (Ур. 2) оптимизации будет посвящена следующая статья. Принятому здесь разделению ЭГСП на две оптимизируемых подсистемы соответствует широко применяемой блочной структуре ЭГСП (рис. 1).
Рис. 1. Блочная структура ЭГСП: БУ - блок управления, ЭМП - электромеханический
преобразователь, ГУ - предварительная ступень усиления, ЗР - золотниковый распределитель (основная ступень усиления), МОС - механическая обратная связь по положению золотника, ИМ - исполнительный механизм (гидроцилиндр), ДОС - датчик основной обратной связи, Иупр - сигнал управления, унагр - координата управляемого объекта, создающего нагрузку на выходное звено ЭГСП
В литературе описано большое количество методов многокритериальной оптимизации [2, 3, 4, 5]|[ каждый из которых имеет свои преимущества и недостатки. В связи с этим одним из этапов решения задачи многокритериальной оптимизации является выбор и обоснование метода.
Выбор параметров ЭГСП основан на том, что конструктор, имея техническое задание, предварительно проектирует его принципиальную схему. Пример такой схемы представлен на рис. 2. Основные конструктивные решения составляющих ЭГСП устройств широко распространены в приводах авиационных систем, станков с программным управлением, робототехнических систем и других современных систем управления [6, 7].
После второго уровня оптимизации, будет обращено внимание на обеспечение в структуре и характеристиках ЭГУ возможности применения методов корректирования динамики всего ЭГСП [7, 8].
Рис. 2. Принципиальная схема ЭГСП с ЭГУ типа «сопло-заслонка» с механической обратной связью по положению золотника (СЗ с МОС)
Процедура выбора параметров подсистемы Ур. 1 состоит из следующих этапов.
1) Составление математической модели привода и определение вектора варьируемых параметров.
2) Выбор критериев для минимизации отклонений характеристик проектируемого ЭГСП с различным сочетанием параметров от технического задания. Набор критериев
можно выразить в виде вектора целевых функций [№)ЛОО.....и ограничений в
виде gJ■ (х) > О
3) Как показывает практика [7], при условии, что частота пропускания ЭГУ СОэгу (частота, при которой фазочастотная характеристика ЭГУ пересекает уровень -90°) в 3.. .5 раз больше собственной частоты исполнительного механизма ^о , можно рассматривать подсистему Ур. 2 в составе всего ЭГСП, как пропорциональное звено с коэффициентом передачи
Qmax
к(Ц =
Ч.тах
Здесь Ртах - расход рабочей жидкости (РЖ) через окна золотникового распределителя при максимальном токе управления 1у.тах в обмотках ЭМП и равенстве нулю перепада давлений между полостями исполнительного механизма. Однако следует
заметить, что подсистема Ур. 2, при наличии в ней колебательного звена с большой постоянной времени, может создавать существенное, по сравнению с допускаемым техническим заданием, фазовое запаздывание. Таким образом, если одним из основных требований к ЭГСП является требования к значениям фазового запаздывания, замена подсистемы Ур. 2 пропорциональным звеном может оказаться некорректной.
В процессе оптимизации подсистемы Ур. 1 необходимо определить конструктивные параметры ЭГСП по критериям Б, (х), ^1.../, /<т , рассматривая подсистему Ур. 1 с подсистемой Ур. 2, представленной пропорциональным звеном.
1 Математическая модель ЭГСП и определение вектора варьируемых параметров
На этапе предварительного синтеза параметров ЭГСП для оценки его устойчивости и быстродействия в целом можно использовать линеаризованную математическую модель [9]. Однако для более подробного исследования характеристик ЭГСП с целью выявления условий, при которых может изменяться качественная природа решений уравнений, описывающих работу спроектированного привода (возникать катастрофы), необходимо затем применить нелинейную математическую модель [10]. Динамические характеристики реального ЭГСП обычно отличаются от характеристик его линейной модели вследствие:
- сухого трения в золотниковом распределителе и исполнительном механизме, или в устройствах, нагружающих выходное звено исполнительного механизма;
- вида функции, описывающей расходно-перепадную характеристику золотникового распределителя, с разрывом в нуле.
Используемые далее переменные состояния ЭГСП и их соответствие составляющим ЭГСП устройствам показаны на рис. 3.
Рис. 3. Схема подсистемы Ур. 1
Уравнение, описывающее электронную часть ЭГСП имеет вид
^у ^¿и ' С^Лп
пр "-ос
и )
"■пг I >
где ^ - ток управления в обмотках ЭМП, ^упр - сигнал управления, ^ос - сигнал от датчика обратной связи, ки - коэффициент передачи электронного усилителя.
Уравнение, описывающее подсистему Ур. 2 при принятом ранее допущении, можно записать в виде
где О - расход РЖ через окна золотникового распределителя при подаче тока управления
^ в случае равенства нулю перепада давлений между полостями исполнительного механизма.
Уравнения, описывающее подсистему Ур. 2, можно также записать в виде
^■ЗОУ) ^лсг ' ^у»
\хз I < X
з.тах>
Л,
где кх, = кд,- Д.тах - коэффициент передачи ЭГУ по ходу золотника, дтах - максимальный
ß
расход через ЭГУ при максимальном открытии окон золотникового распределителя хзтах.
Расход РЖ, протекающей через окно золотникового распределителя,
определяется известной из гидравлики формулой [7]
(} = Ъок, х31г?,р) ■ 7(рБХ-рБЫХ),
Где кз (р, Г, Ьок, хз, V, р) = кз(Г, хз, V, р)Ьок-^2/ р - полученная экспериментально
зависимость проводимости окна золотникового распределителя от плотности РЖ р, температуры РЖ 1, суммарной ширины окон золотникового распределителя Ьок, перемещения золотника Хз, направления потока жидкости V; р = рвых / рвх -
относительное противодавление, рвх и рвых - давления на входе и выходе окна
золотникового распределителя; к'з (1, Хз, V, р) - удельная проводимость окна
золотникового распределителя.
Уравнения, описывающие расход РЖ через окна золотникового распределителя, необходимо дополнить уравнениями, основанными на гипотезе неразрывности течения сжимаемой РЖ:
Здесь 5Л - площадь поршня гидроцилиндра (исполнительного механизма), уп -перемещение поршня гидроцилиндра, У\ и ^2 - объемы РЖ в каждой из полостей гидроцилиндра с присоединенными к ним гидролиниями, соединяющими золотниковый распределитель с гидроцилиндром, В - модуль объемной упругости РЖ, точкой над переменными в уравнении обозначена производная по времени.
При этом значения давлений в полостях исполнительного механизма ограничены неравенствами
Уравнение движения поршня гидроцилиндра с учетом нагрузки, действующей на его шток, имеет вид [9]
где сн - жесткость позиционной нагрузки, Ртр = Ятр + ктр-уп - трение в нагрузке на штоке
и в уплотнениях, состоящее из сухого трения величиной и вязкого трения с коэффициентом вязкого трения кр, тн - масса нагрузки, перемещаемой штоком
гидроцилиндра.
В соответствии с записанными выше уравнениями, варьируемыми параметрами для подсистемы Ур. 1 являются сила трения Р-ф, коэффициент передачи электронного
усилителя кщ, величина зоны нечувствительности ЭГУ. Также к варьируемым параметрам стоит отнести изменение передаточного отношения механизма, соединяющего выходное звено привода с управляемым объектом, в к раз.
2 Выбор критериев
Оптимизация силовой части ЭГСП осуществляется по следующим критериям [1]. 1) Целевые функции, характеризующие соотношение мощности и массы исполнительного механизма [11]. Для ЭГСП с поступательным движением выходного звена такой целевой функцией может быть безразмерное соотношение
Где т - масса ЭГСП; Ув - скорость поступательного движения выходного звена ЭГСП, Рв - максимальная сила, действующая на выходное звено ЭГСП, утах - максимальное перемещение выходного звена.
В техническом задании требования к предельной механической характеристике ЭГСП определяются из соображений обеспечения необходимых законов управления управляемым объектом и могут предъявляться в следующем виде.
а) Требования к максимальной развиваемой силе на выходном звене ЭГСП Рв и требования к максимальной скорости движения выходного звена ЭГСП Ув при отсутствии силы на выходном звене.
При изменении передаточного отношения механизма, соединяющего выходное звено ЭГСП с управляемым объектом, в к раз требования к ЭГСП изменятся следующим
образом: Рв' = Рв / к , Ув = ЬУв , Утах' = к'Утах .
Также для выбранной конструктивной схемы ЭГСП можно в первом приближении определить из соображений обеспечения прочности зависимость [12]
т = т(к, Ятр, рп).
На основе изложенного выше, можно записать целевую функцию для нового проектного варианта ЭГСП, как
VB2
где-в--постоянная величина.
P-y
в у max
б) Требования к скорости движения выходного звена ЭГСП Vp при заданной силе на выходном звене ЭГСП Pp и ограничение на значение максимальной развиваемой силы на выходном звене ЭГСП Рв.
Аналогично изложенному выше, можно получить целевую функцию в виде
однако в отношении —V— постоянной величиной остается ymax, а остальные величины
P-у
в У max
можно записать, как Vв(Pp, Sn, Рп), Рв =Рв(Рр, Sh, Рп).
2) Динамические показатели, выбираемые в соответствии с техническим заданием на ЭГСП. Например, для ЭГСП, используемых в авиационных системах управления, они представлены в виде [7]:
- требований на величину резонансного пика на логарифмической частотной характеристике Ьрез (как правило, не более 1...2 Дб), и частоту пропускания ^пп при логарифмическом значении амплитуд -3.-6 Дб для замкнутого контура привода;
- требований к фазочастотным характеристикам ЭГСП в виде ограничений на максимальный фазовый сдвиг перемещения выходного звена относительно входного сигнала в заданном диапазоне частот;
- кроме того, при необходимости демпфирования колебаний управляемого объекта в случае отсутствия специальных демпфирующих устройств могут предъявляться требования к динамической жесткости ЭГСП [8]. Динамическая жесткость определяется, как комплексное отношение возмущающей силы, действующей на выходное звено ЭГСП, к перемещению выходного звена, вызванному этим возмущением. Для обеспечения демпфирования управляемого объекта ЭГСП должен во всем диапазоне частот возмущающей силы иметь фазочастотную характеристику динамической жесткости максимально близкую к +90° и значение динамической жесткости не меньшее определенного в техническом задании.
В случае невозможности выполнения требований к динамическим характеристикам путем варьирования конструктивных параметров, целесообразно
вводить в структуру привода корректирующие устройства [7].
Поскольку работу ЭГСП описывает нелинейная математическая модель, при подаче управляющего гармонического сигнала, выходной сигнал имеет форму, отличную от гармонической, а частотные характеристики ЭГСП зависят от амплитуды входного сигнала. Частотные характеристики в таком случае определяют на основе первой гармоники разложения выходного сигнала в ряд Фурье для входных гармонических сигналов с различной амплитудой.
Расчет частотных характеристик ЭГСП целесообразно производить в следующем порядке.
а) Руководствуясь техническим заданием, выбрать:
- диапазон рабочих частот Wmax, Гц;
- шаг по частотной оси, представленной в логарифмическом масштабе A^iog (шаг по частоте
- набор амплитуд входного сигнала А = \А1,А1,..., Ап].
б) Для каждого сочетания амплитуды и частоты входного гармонического сигнала произвести расчет выходного сигнала.
в) Определить амплитуду и фазу первой гармоники разложения выходного сигнала в ряд Фурье (далее «первая гармоника»). Для сокращения времени вычислений процесс расчета выходного сигнала необходимо останавливать по условию одновременной малости первых и вторых производных амплитуды и фазы первой гармоники. В связи с чувствительностью результатов численного определения производных к дискретности процесса расчета, рационально производную заменить звеном с передаточной функцией
v 1
W (s) =-, dT < T <-,
\ / гт1 1 5 max 5
T'S + 1 ®max
где dTmax - максимальный шаг интегрирования, s = d.
dt
3) Энергетические характеристики.
Критерием, характеризующим энергетическую эффективность, может быть количество энергии, используемое для регулирования.
Также, с точки зрения обеспечения оптимальных энергетических характеристик, следует отметить случай, когда требования к предельной механической характеристике ЭГСП заданы скоростью выходного звена ЭГСП Vp при требуемой нагрузке Pp, а также минимальной Pe.min и максимальной Pe.max развиваемой ЭГСП силы. При таком сочетании
требований технического задания можно получить чрезмерно большой потребляемый расход РЖ в режиме холостого хода. Для его уменьшения следует повышать развиваемую силу. В случае Рв > Рвтах в структуру ЭГСП для ограничения развиваемой силы необходимо вводить нелинейную обратную связь по перепаду давления в гидродвигателе, либо применять на выходном звене гидропружину [12].
3 Решение многокритериальной задачи
Как показано в ряде работ, для решения задач многокритериальной оптимизации (МКО) ЭГСП целесообразно применять метод ЛП-поиска [3, 6]11или генетические алгоритмы (ГА) [13]. ГА универсальны и просты в применении. Многокритериальные модификации ГА [4] позволяют получать более полный фронт Парето по сравнению с методом ЛП-поиска при меньшем участии конструктора.
В качестве целевых функций при решении задачи МКО целесообразно выбирать два-три основных критерия из вышеперечисленных, остальные же следует отнести к ограничениям, накладываемых техническим заданием.
Для выбора оптимального проектного варианта в задаче МКО, как правило, необходимо сравнивать векторные критерии, элементы которых имеют различные размерности. В связи с этим, необходимо нормировать множество, составляющее полученный в результате применения выбранного алгоритма МКО фронт Парето. Если ни один из элементов фронта Парето не принимает нулевого значения, целесообразно нормировать пространство критериев следующим образом:
Здесь fmin - минимальное значение критерия f, fi(x) - значение критерия f при сочетании конструктивных параметров X, n - число критериев.
В качестве примера рассмотрим задачу оптимизации составной части ЭГСП -электрогидравлического агрегата, включающего в себя ЭГУ, ИМ и ДОС (рис. 1), со следующими параметрами:
развиваемая сила Ръ = 17000 Н,
максимальное перемещение выходного звена утах = 70 ■ 10~3 м, максимальная скорость выходного звена vB = 0,2 м/с ,
запаздывание по фазе в разомкнутом контуре ЭГСП не более относительно управляющего сигнала амплитудой, равной
половине от максимального, на частотах, соответственно, о) = [5 10 15 20] Гц.
Диапазон изменения варьируемых параметров и полученное оптимальное сочетание параметров ЭГУ приведены в табл. 1.
Таблица 1
Параметр, единица измерения Минимальное значение параметра Максимальное значение параметра Оптимальное сочетание параметров
к 0,7 1,3 0,72
к-ш' А/В 0,09 0,27 0,26
йгр, Н 100 600 165
/Срр, Н ■ с/м 50 300 300
рп,Па 15■106 21■106 19,5 ■ 106
В качестве критериев выбраны:
1) целевая функция, характеризующая отношение массы к мощности, в виде
2) целевая функция, характеризующая близость фазочастотной характеристики к оси абсцисс, в виде
При помощи многокритериальной модификации ГА получен аппроксимированный фронт Парето, затем все элементы фронта Парето нормированы. Результат оптимизации представлен на рис. 4.
1.4
1.2
о
0
1
а
о
д и
о
'Я
К
к:
ш
а
о
£
0 1) к
1
О)
(Г (Й к
п
(Г О Е
ас
ей рр
О
^0.2
0.8
0.6
0.4
Рч
о
И
Нормированный фронт Парето -1-1-1—
э .......•■•.......;.................
о
о
о
........../ч.....
О
Оптимальное решение
0 0.002 0.004 0.006 0.008 0.01 0.012
Нормированное значение критерия фаз о частотной характеристики
Рис. 4. Нормированный фронт Парето
Видно, что в результате варьирования конструктивных параметров не удалось существенно улучшить фазочастотную характеристику проектируемого ЭГСП. Таким образом, оптимальным проектным вариантом является сочетание параметров, при котором принимает минимальное значение целевая функция, характеризующая отношение массы к мощности.
Заключение.
В полученном решении изменение передаточного отношения механизма, соединяющего выходное звено ЭГСП с управляемым объектом в к раз, и коэффициент передачи электронного усилителя ки можно воплотить в реальной конструкции при отсутствии дополнительных неформализуемых ограничений (ограничения по элементной базе, ограничения связанные с компоновкой ЭГСП и системы, в которую он входит и т.п.).
Полученное значение давления питания рп можно реализовать, если разработчик ЭГСП проектирует и систему, в которую входит гидропривод. В таком случае оптимизацию необходимо проводить для всей системы. Иначе, полученное значение давления питания имеет рекомендательный характер для разработчика системы.
Сухое трение Rp и коэффициент вязкого трения kp зависят от конструкции уплотнительных узлов и типа нагрузки на выходном звене. Данные параметры чаще всего имеют некоторую степень неопределенности и сложно реализовать их при проектировании и производстве. Таким образом, при отсутствии отдельных требований к величине трения стоит отнести полученные значения параметров Rp и ктр к рекомендательным.
После определения действительной величины нагрузки трения, действующей на выходное звено ЭГСП, и выбора величины давления питания целесообразно повторно провести оптимизацию параметров ЭГСП.
Рассмотренный в статье первый уровень оптимизации ЭГСП с дроссельным регулированием соответствует концепции оптимизации, изложенной в [1]. Научная новизна заключается в постановке и решении задачи многокритериальной оптимизации ЭГСП.
Список литературы
1. Замараев Д.С., Попов Д.Н. Концепция оптимизации электрогидравлического следящего привода с дроссельным регулированием // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. № 6. DOI: 10.7463/0613.0569281
2. Карпенко А.П., Федорук В.Г. Аппроксимация функции предпочтений лица, принимающего решения, в задаче многокритериальной оптимизации. 3. Методы на основе нейронных сетей и нечеткой логики // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2008. № 4. Режим доступа: http://technomag.bmstu.ru/doc/86335.html (дата обращения 26.09.2013).
3. Статников Р.Б., Соболь И.М. Выбор оптимальных параметров в задачах со многими критериями. 2-е изд., перераб. и доп. М.: Дрофа, 2006. 175 с.
4. Deb Kalyanmoy. Multi-Objective Optimization Using Evolutionary Algorithms. John Wiley & Sons, 2001. 518 p.
5. Arora J.S., Marler R.T. Survey of multi-objective optimization methods for engineering // Struct Multidisc Optim. 2004. № 26. P. 369-395.
6. Боровин Г.К., Попов Д.Н. Многокритериальная оптимизация гидросистем: учеб. пособие. М.: Изд-во МГТУ им. Н.Э. Баумана, 2007. 94 c.
7. Попов Д.Н., Асташев В.К., Густомясов А.Н. и др. Машиностроение: энциклопедия. В 40 т. Т. IV-2 : Электропривод. Гидро- и виброприводы. Кн. 2: Гидро- и виброприводы / под ред. В.К. Асташева, Д.Н. Попова. М.: Машиностроение, 2012. 303 c.
8. Редько П.Г., Ермаков С.А., Константинов С.В. Электрогидравлические рулевые приводы систем управления полетом маневренных самолетов. М.: "Янус-К", 2006. 316 с.
9. Попов Д.Н. Механика гидро- и пневмоприводов: учебник для машиностроительных вузов. 2-е изд. М.: Изд-во МГТУ им. Н.Э. Баумана, 2002. 320 с.
10. Гилмор Р. Прикладная теория катастроф. 2-е изд. М.: Мир, 1984. 285 с.
11. Попов Д.Н. Оценка эффективности и оптимальное проектирование гидроприводов // Вестник машиностроения. 1986. № 9. С. 20-23.
12. Оболенский Ю.Г., Ермаков С.А., Сухоруков Р.В. Введение в проектирование систем авиационных рулевых приводов. М.: Изд-во ГУП г. Москвы "Окружная газета ЮЗАО", 2011. 344 с.
13. Andersson J., Krus P., Wallace D. Multi-objective optimization of hydraulic actuation systems // Proceedings of the 2000 ASME Design Automation Conference (DETC'00) (Baltimore, MD, 10-13 September 2000). DAC-14512. Режим доступа: http://cadlab.mit.edu/publications/00-andersson-asme-dt/dontindex/00-andersson.pdf
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Optimization of power subsystem of electrohydraulic servo drive with throttling governing # 11, November 2013 DOI: 10.7463/1113.0623966 Zamaraev D.S., Popov D.N.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
popov@bmstu.ru d.stalky@gmail.com
This paper describes optimization of the power subsystem of an electrohydraulic servo drive with throttling governing. The electrohydraulic servo drive was considered together with an electrohydraulic power amplifier replaced by a proportional link. Variable parameters were selected and objective functions for design alternatives were proposed in order to optimize the electrohydraulic servo drive. An example of two-criterion optimization of the power subsystem was considered in this paper. These criteria were nearness of phase-response characteristic to the axis of abscissae and the mass to power ratio. Multi-criterion modification of genetic algorithm was used to obtain the Pareto front. The described concept is based on generalization of optimization methods for the electrohydraulic servo drive known from literature references and includes new proposals on the structure of the optimization algorithm.
Publications with keywords: drosselling hydraulic actuator, multicriterion optimization, genetic algorithm, Electrohydraulic Servo actuator
Publications with words: drosselling hydraulic actuator, multicriterion optimization, genetic algorithm, Electrohydraulic Servo actuator
Referencrs
1. Zamaraev D.S., Popov D.N. Kontseptsiya optimizatsii elektrogidravlicheskogo sledyashchego privoda s drossel'nym regulirovaniem [Concept of optimization of an electrohydraulic servo drive with throttle regulation]. Nauka i obrazovanieMGTUim. N.E. Baumana [Science and Education of the Bauman MSTU], 2013, no. 6. DOI: 10.7463/0613.0569281
2. Karpenko A.P., Fedoruk V.G. Approksimatsiya funktsii predpochteniy litsa, prinimayushchego resheniya, v zadache mnogokriterial'noy optimizatsii. 3. Metody na osnove neyronnykh setey i nechetkoy logiki [Approximation of functions of preferences of decision
maker in the problem of multi-criteria optimization. 3. Methods on the basis of neural networks and fuzzy logic]. Nauka i obrazovanie MGTUim. N.E. Baumana [Science and Education of the Bauman MSTU], 2008, no. 4. Available at: http://technomag.bmstu.ru/doc/86335.html , accessed 26.09.2013.
3. Statnikov R.B., Sobol' I.M. Vybor optimal'nykhparametrov v zadachakh so mnogimi kriteriiami [Choice of the optimum parameters in problems with many criteria]. Moscow, Drofa, 2006. 175 p.
4. Deb Kalyanmoy. Multi-Objective Optimization Using Evolutionary Algorithms. John Wiley & Sons, 2001. 518 p.
5. Arora J.S., Marler R.T. Survey of multi-objective optimization methods for engineering. Struct. Multidisc. Optim, 2004, vol. 26 , no. 6, pp. 369-395. DOI: 10.1007/s00158-003-0368-6.
6. Borovin G.K., Popov D.N. Mnogokriterial'naia optimizatsiia gidrosistem [Multi-criteria optimization of hydraulic systems]. Moscow, Bauman MSTU Publ., 2007. 94 p.
7. Popov D.N., Astashev V.K., Gustomiasov A.N., et al. Mashinostroenie: entsiklopediia. V 40 t. T. 4-2: Elektroprivod. Gidro- i vibroprivody. Kn. 2: Gidro- i vibroprivody [Mechanical Engineering: Encyclopedia. In 40 vols. Vol. 4-2: Electric drive. Hydraulic and vibro drives. Book 2: Hydraulic and vibro drives]. Moscow, Mashinostroenie, 2012. 303 p.
8. Red'ko P.G., Ermakov S.A., Konstantinov S.V. Elektrogidravlicheskie rulevyeprivody sistem upravleniiapoletom manevrennykh samoletov [Electro-hydraulic steering drives of systems of control of flight of maneuverable aircraft]. Moscow, "Ianus-K", 2006. 316 p.
9. Popov D.N. Mekhanika gidro- ipnevmoprivodov [Mechanics of hydraulic and pneumatic drives]. Moscow, Bauman MSTU Publ., 2002. 320 p.
10. Gilmore R. Catastrophe Theory for Scientists and Engineers. New York, Dover Publications, 1981. (Russ. ed.: Gilmore R. Prikladnaya teoriya katastrof. V 2 t. T. 2. Moscow, Mir, 1984. 285 p.).
11. Popov D.N. Otsenka effektivnosti i optimal'noe proektirovanie gidroprivodov [Evaluation of effectiveness and optimum designing of hydraulic drives]. Vestnik mashinostroeniia, 1986, no. 9, pp. 20-23.
12. Obolenskii Iu.G., Ermakov S.A., Sukhorukov R.V. Vvedenie vproektirovanie sistem aviatsionnykh rulevykhprivodov [Introduction to the design of systems of aircraft steering gear]. Moscow, "Okruzhnaia gazeta IuZAO" Publ., 2011. 344 p.
13. Andersson J., Krus P., Wallace D. Multi-objective optimization of hydraulic actuation systems. Proceedings of the 2000 ASME Design Automation Conference (DETC'00), Baltimore, MD, 10-13 September, 2000, DAC-14512. Available at: http://cadlab.mit.edu/publications/00-andersson-asme-dt/dontindex/00-andersson.pdf , accessed 26.09.2013.
