Другие журналы

научное издание МГТУ им. Н.Э. Баумана

НАУКА и ОБРАЗОВАНИЕ

Издатель ФГБОУ ВПО "МГТУ им. Н.Э. Баумана". Эл № ФС 77 - 48211.  ISSN 1994-0408

Идентификация параметров математической модели электрогидравлического привода

# 06, июнь 2013
DOI: 10.7463/0613.0574706
Файл статьи: Scherbachev_Р.pdf (746.24Кб)
автор: Щербачев П. В.

УДК 62-522.2

Россия, МГТУ им. Н.Э. Баумана

scherbachev.p.v@gmail.com

 

Введение

Математическое моделирование на сегодняшний день является неотъемлемой частью проектирования и исследования гидравлических приводов. Наличие математической модели позволяет провести численный эксперимент задолго до появления опытного образца изделия, получить его ожидаемые характеристики, либо оптимизировать уже имеющееся устройство. В случае с электрогидравлическими приводами часто встают задачи синтеза регулятора, отработки специального закона управления. Эти задачи на первом этапе решаются на математических моделях приводов. В связи этим актуальным становится вопрос адекватности математической модели.

Физические процессы, протекающие при работе электрогидравлических приводов, в математической модели описываются системой алгебраических и дифференциальных уравнений. В зависимости от решаемой задачи необходимо выбирать степень приближения модели к действительности. Чем точнее требуется описать процесс или явление, тем больше численных значений параметров разного рода нужно знать, будь то характеристики материалов, точность изготовления отдельных элементов и т.д. В тех случаях, когда стоит задача оптимизации имеющегося устройства, большая часть параметров системы априори известна. Однако, даже имея готовый работоспособный агрегат, не всегда удается путем непосредственных измерений получить численные значения всех интересующих параметров. В электрогидравлических приводах одним из таких проблемных мест является описание электрогидравлических усилителей мощности (ЭГУ).

В следящих электрогидравлических приводах часто применяют двухкаскадные ЭГУ золотникового типа. В данной статье речь пойдет об электрогидравлическом приводе вращательного движения, схема которого изображена на рисунке 1. Коленчатый вал 1 приводится во вращение двумя гидроцилиндрами 3 через шатуны 4. Согласованная работа цилиндров обеспечивается подачей специальных сигналов на два электрогидравлических усилителя (на схеме не показаны), выходы которых соединены с полостями цилиндров. На качество характеристик такого привода большое влияние оказывает качество исполнения ЭГУ, в том числе точность изготовления золотника. При составлении математической модели необходимо качественно и количественно описать золотниковые пары ЭГУ.

 

 

Рисунок 1 - Схема электрогидравлического привода вращательного движения

 

В [1] предложен метод описания точности изготовления цилиндрического золотникового распределителя. Схема приведена на рисунке 2. Основными параметрами точности являются:

·       S – перекрытия кромок;

·       R – радиусы скругления кромок;

·       δ – радиальный зазор в паре золотник-гильза.

Параметры R и δ являются приблизительно одинаковыми для всех рабочих кромок. Параметр S может достаточно сильно варьироваться в силу технологии изготовления. Поэтому при составлении математической модели следует учесть индивидуальный закон открытия каждой из дросселирующих щелей золотникового распределителя.

 

 

Рисунок 2 – Схема цилиндрического золотникового распределителя.

 

Математическая модель

Расчетная схема, принятая для составления математической модели, приведена на рисунке 3.

 

Рисунок 3 – Расчетная схема.

 

Уравнения математической модели

1.     Уравнения, описывающие электрогидравлические усилители.

ЭГУ представлены дифференциальным уравнением первого порядка согласно [2].

 ,

где  – постоянная времени гидроусилителя;

и  – смещения золотников первого и второго ЭГУ;

 – коэффициент усиления;

 и  – входные сигналы ЭГУ:

где  – текущий угол поворота вала;

 – заданный угол поворота вала.

Методы формирования управляющих сигналов  и  описаны в [4].

Смещение золотника ограничивается неравенствами:

,

где  – максимальное смещение золотника.

Из-за наличия неточности изготовления золотников величина открытия дросселирующих щелей не равна смещению золотника [1].

,

где  – номер ЭГУ

 – открытие -ой щели -го распределителя (см. схему на рисунке 2);

 – радиус скругления кромок;

 – радиальный зазор в паре золотник-гильза;

 – перекрытия кромок.

 – вспомогательные величины:

2.     Расходы через кромки золотниковых распределителей.

где  – давление питания в системе;

 – давление в сливе системы;

 – давление в полости “A–го цилиндра;

 – давление в полости “B–го цилиндра;

 – расход из линии нагнетания в полость “A”  –го цилиндра;

 – расход из полости “A–го цилиндра в линию слива;

 – расход из линии нагнетания в полость “B”  –го цилиндра;

 – расход из полости “B–го цилиндра в линию слива;

 – проводимость щели золотникового распределителя:

где  – коэффициент расхода щели;

 – ширина окон в гильзе золотникового распределителя;

 – количество окон в гильзе золотникового распределителя у каждой кромки;

 – плотность рабочей жидкости.

3.     Расходы через обратные клапаны.

где  – расход из полости “A–го цилиндра в линию нагнетания через обратный клапан;

 – расход из линии слива в полость “A–го цилиндра через обратный клапан;

 – расход из полости “B–го цилиндра в линию нагнетания через обратный клапан;

 – расход из линии слива в полость “B–го цилиндра через обратный клапан;

 – проводимость обратного клапана в линию нагнетания;

 – проводимость обратного клапана из линии слива.

4.     Уравнения баланса расходов, приведенные к форме Коши.

где  – приведенный модуль упругости жидкости;

 – объем сжатия в полости “A–го цилиндра;

 – объем сжатия в полости “B–го цилиндра;

 – смещение поршня –го цилиндра;

 – площадь поршня.

5.     Сила, развиваемая поршнем гидроцилиндра.

где  – суммарная сила сухого и вязкого трения.

где  – сила Кулоновского трения;

 – сила страгивания;

 – коэффициент вязкого трения;

 – пороговое значение скорости.

Подробное описание принятой модели силы трения представлено в  [3] и [5].

6.     Движущий момент привода.

Момент от поршня одного гидроцилиндра рассчитывается исходя из кинематики механизма с учетом шатуннного влияния.

где  – длина кривошипа;

 – длина шатуна.

7.     Уравнение движения нагрузки (приведенное к форме Коши).

где  – приведенный к валу момент инерции;

 – внешний нагружающий момент.

8.     Уравнения кинематической связи.

Решение уравнений математической модели проводилось численным методом при помощи алгоритма, реализованного на языке программирования C++. Численные значения параметров в системе СИ приведены в таблице 1.

Таблица 1

Параметр

Значение

Размерность

20

МПа

0.1

МПа

0.06

с

1.59*10-4

м/В

3.0*10-4

м

9*10-6

м

2*10-6

м

3.41*10-4

м2/(с*Па0.5)

2.67*10-7

м3/(с*Па0.5)

2.67*10-7

м3/(с*Па0.5)

1.5*109

Па

3*10-4

м3

3*10-4

м3

1.15*10-3

м2

1000

Н

100

Н*с/м

10-6

м/с

600

Н

0.02

м

0.4

м

0.328

кг*м2

20

Н*м

 

Идентификация параметров

Величины перекрытий в паре «золотник-гильза», входящие в уравнения математической модели, изначально неизвестны. Непосредственно измерить их также не представляется возможным. В таком случае для их определения целесообразно воспользоваться предлагаемым методом идентификации. Для этого необходимо определяемые величины представить набором варьируемых параметров . Также нужно выбрать критерий несоответствия  математической модели и моделируемого объекта. Идентичность модели и объекта будет достигнута при обеспечении минимума критерия  на множестве параметров . Таким образом, метод идентификации параметров представляет собой задачу многомерной оптимизации. Блок-схема описанного алгоритма решения задачи представлена на рисунке 4. В качестве численного метода поиска минимума был использован метод «деформируемого многогранника», известный также как метод Нелдера-Мида.

Наибольшее влияние неточности изготовления золотниковых распределителей оказывают на характер изменения давлений в полостях цилиндров. Поэтому величина несовпадения  давлений в полостях цилиндров, полученных экспериментально и теоретически, выбрано в качестве вышеуказанного критерия .

 

 

Рисунок 4 – Блок-схема алгоритма идентификации.

 

В ходе эксперимента привод работал в режиме слежения по угловой скорости вращения выходного звена. Зависимость давлений в полостях от времени, полученная экспериментально, показана на рисунке 5.

 

Рисунок 5 – Экспериментальные данные.

 

Для иллюстрации влияния перекрытий на изменение давления на рисунке 6 показаны зависимости давлений в полостях цилиндров, полученные при расчете математической модели с нулевыми значениями перекрытий.

 

 

Рисунок 6 – Результаты моделирования.

 

Во время проведения эксперимента данные с датчиков давления фиксировались с частотой 1000Гц. Критерий совпадения характеристик был получен суммированием разниц между давлением, полученным экспериментально, и давлением, полученным при моделировании, в течение одного оборота с интервалом по времени 0.001с (см. рисунок 7). Математическое выражение сформированного критерия следующее:

 

 

Рисунок 7 – Формирование критерия.

 

Результаты

Полученные в результате процедуры идентификации значения перекрытий в золотниковых парах представлены в таблице 2.

Таблица 2

ЭГУ 1

ЭГУ 2

№ кромки

S, мкм

№ кромки

S, мкм

1

8.0

1

-0.6

2

-5.4

2

19.5

3

13.8

3

-0.1

4

-5.1

4

15.1

 

Следует отметить, что полученные цифры показывают не только неточность изготовления золотниковых пар, но также включают в себя смещение золотника от нейтрального положения при нулевом входном сигнале – дрейф нуля ЭГУ.

На рисунках 8 и 9 представлены графики изменения давлений в полостях цилиндров с учетом вычисленных значений перекрытий. Пунктирные линии отображают экспериментальные данные, сплошные – результаты моделирования.

 

Рисунок 8 – Сопоставление результатов (цилиндр 1).

 

 

Рисунок 9 – Сопоставление результатов (цилиндр 2).

 

Заключение

По полученным результатам, основываясь на хорошем совпадении экспериментальных данных и результатов моделирования, можно сделать вывод о том, что предложенный алгоритм идентификации параметров эффективен.

Для проведения поиска минимума критерия в предложенном алгоритме удобно использовать метод «деформируемого многогранника», т.к. в этом случае не требуется вычислять частные производные минимизируемого критерия по всем варьируемым параметрам.

Для получения адекватных математических моделей электрогидравлических приводов с дроссельным регулированием следует учитывать микрогеометрию кромок золотниковых распределителей.

Список литературы

1.               Фомичев В.М. Проектирование цилиндрических золотниковых распределителей следящих электрогидравлических приводов : учеб. пособие. M.: Изд-во МГТУ им. Н.Э. Баумана, 2009. 19 с.

2.               Фомичев В.М. Проектирование электрогидравлических усилителей следящих приводов : учеб. пособие. M.: Изд-во МГТУ им. Н.Э. Баумана, 2009. 42 с.

3.               Ломакин В.О., Шумилов И.С., Щербачев П.В. Математическое моделирование электрогидравлического следящего привода высокоточного регулирования вращательного движения // Известия ВУЗов. Машиностроение. 2011. № 10. C. 38-45.

4.               Щербачев П.В., Семенов С.Е. Электрогидравлический привод с дроссельным регулированием с повышенной энергоэффективностью  // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн.  2012. № 10. DOI: 10.7463/1012.0465528

5.               Руководство пользователя. Matlab 7.12.0. Режим доступа: http://www.mathworks.com/  (дата обращения 30.05.2013). 

Поделиться:
 
ПОИСК
 
elibrary crossref ulrichsweb neicon rusycon
 
ЮБИЛЕИ
ФОТОРЕПОРТАЖИ
 
СОБЫТИЯ
 
НОВОСТНАЯ ЛЕНТА



Авторы
Пресс-релизы
Библиотека
Конференции
Выставки
О проекте
Rambler's Top100
Телефон: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)
  RSS
© 2003-2022 «Наука и образование»
Перепечатка материалов журнала без согласования с редакцией запрещена
 Тел.: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)