Научная статья на тему 'Методика идентификации струйно-золотниковой гидравлической рулевой машины'

Методика идентификации струйно-золотниковой гидравлической рулевой машины Текст научной статьи по специальности «Строительство и архитектура»

CC BY
261
63
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / ИДЕНТИФИКАЦИЯ / МЕТОДИКА / СТАТИСТИЧЕСКАЯ ОБРАБОТКА / ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Месропян Арсен Владимирович, Широкова Ксения Александровна, Целищев Владимир Александрович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Identification method of the jet-spool hydraulic steering machine

Mathematical modeling problems of the jet-spool hydraulic steering machine are observed. Principle circuit of the hydro-drive with the jet-spool hydraulic actuator, the mathematical model and computational modeling data, obtained in Maple application package, are presented. Identification methods of the jet-spool hydraulic steering machine are presented.

Текст научной работы на тему «Методика идентификации струйно-золотниковой гидравлической рулевой машины»

МАШИНОСТРОЕНИЕ • ГИДРАВЛИЧЕСКИЕ МАШИНЫ, ГИДРОПНЕВМОАГРЕГАТЫ

УДК 629.062-82

А. В. МЕСРОПЯН, К. А. ШИРОКОВА, В. А. ЦЕЛИЩЕВ

МЕТОДИКА ИДЕНТИФИКАЦИИ СТРУЙНО-ЗОЛОТНИКОВОЙ ГИДРАВЛИЧЕСКОЙ РУЛЕВОЙ МАШИНЫ

Рассмотрены вопросы математического моделирования струйно-золотниковой гидравлической рулевой машины. Представлена принципиальная схема гидропривода со струйно-золотниковым гидроусилителем, математическая модель и результаты численного моделирования в пакете прикладных программ Мар1е. Представлена методика идентификации. Математическая модель; идентификация; методика; статистическая обработка; численное моделирование

В настоящее время накоплена обширная база экспериментальных данных по снятию характеристик струйных гидравлических рулевых машин (СГРМ). Важной задачей при проектировании СГРМ является разработка оптимальной методики проектирования, учитывающей результаты обработки, обобщения и анализа экспериментальных данных по снятию статических и динамических характеристик СГРМ.

Критериями, определяющими оптимальность методики, могут выступать адекватные модели, временной фактор, универсальность методики, заключающаяся в возможности моделирования различных компоновочных схем СГРМ.

Применительно к СГРМ конечной целью разработки методики проектирования является потребность в инструментарии, позволяющем быстро и с высокой степенью точности и достоверности рассчитать характеристики и параметры одно- и двухкаскадных схем электрогидроусилителей.

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

Близость математической модели реальному объекту обуславливается ее сложностью и более неоднозначно взаимосвязанное влияние различных внешних и внутренних возмущающих факторов, влияющих на ее работо-

способность при исследовании динамических характеристик СГРМ.

Существует два способа формирования математических моделей. Первый способ состоит в том, чтобы разбить систему на такие подсистемы, свойства которых известны. Формальное математическое объединение этих подсистем позволяет получить модель всей системы. Такой подход и называется моделированием. Основной прием сводится к структуризации процесса в виде блок-схем, блоки которых состоят из более простых элементов.

При другом способе построения математической модели непосредственно используются экспериментальные данные. По существу, необходимо по некоторой совокупности параметров, измеряемых в процессе испытаний, воспроизвести недостающие характеристики струйно-золотниковой гидравлической рулевой машины (СЗГРМ) (или иного объекта), установить функциональные зависимости между входными и выходными параметрами, а также получить данные, необходимые для идентификации рулевой машины.

Задача идентификации может быть сформулирована следующим образом [1]. Задан объект, в процессе функционирования которого одновременно могут быть измерены его входная и выходная переменные. По результатам измерения и необходимо построить модель заданного объекта, т. е. найти оператор (точнее, его приближенное значение, оценку), ставящий в соответствие выходную у(£) и входную ж(£) функции. Если характеристикой рассматриваемого объекта является оператор Аг:

то задача идентификации заключается в определении его оценки А*г, т. е. (см. рис. 1.):

у* (г) = А1-х{і).

(2)

х(4

-Объект— -Модель

Рис. 1. Схема задачи идентификации

Пусть существует объект и его модель, на входе действует одна и та же переменная ж(£). Истинным выходом объекта является у(£), а на выходе модели .

Очевидно, что модель «беднее» реальной действительности, что и определяет существующие различия между у{Ь) и у*(і). Задача идентификации заключается в том, чтобы по реализациям ж(£) и у{Ь) найти оценку истинного оператора .

Таким образом, идентификация позволяет спроектировать СЗГРМ с параметрами, при которых достигается близость характеристик реальной установки и характеристик, построенных при помощи математического моделирования на ЭВМ с учетом экспериментальных данных. Этапы процедуры идентификации представлены на рис. 2.

Методику идентификации СЗГРМ можно представить как совокупность следующих этапов:

формулирование цели и задач исследования;

составление принципиальной схемы рулевой машины;

описание принципа действия СЗГРМ и основных физических процессов происходящих в рулевой машине;

выбор структуры модели; выбор параметров, которые задаются; выбор параметров, которые необходимо определить;

• расчет статических и динамических характеристик;

• статистическая обработка результатов экспериментов;

• построение стохастической математической модели;

• анализ полученных результатов и выдача рекомендаций по проектированию и доводки рулевых машин.

Основными элементами СЗГРМ являются: электромеханический преобразова-

тель (ЭМП), струйный гидрораспределитель (СГР), золотниковый гидрораспределитель (ЗГР) и гидродвигатель. На рис. 3 представлена принципиальная схема струйно-золотникового гидроусилителя. Электрический управляющий сигнал, преобразованный и усиленный в промежуточных устройствах автоматики, поступает в ЭМП. Якорь ЭМП смещается на угол , поворачивая в ту или иную сторону струйную трубку, сопло которой перекрывается с приемными окнами распределителя — приемной платы, что приводит к возникновению перепада давлений в золотнике. Когда сигнал задатчика будет скомпенсирован сигналом цепи обратной связи, струйная трубка вернется в исходное положение.

На рис. 4,5,6 представлены основные компоненты двухкаскадной схемы струйно-золотникового гидроусилителя (СЗГУ).

Объем полости струйной камеры между срезами питающего сопла и приемного канала обычно намного меньше суммарного объема жидкости в полостях гидропривода. Влияние локальных изменений скорости на участке между срезами сопла и канала также незначительно (рис. 5). Поэтому при анализе динамических свойств струйно-золотникового гидроусилителя можно пренебречь сжатием и инерционностью жидкости на участке между срезами питающего сопла и приемного канала.

На рис. 7,8 для большей наглядности процессов и явлений, происходящих в проточной части гидроусилителя, представлена область расчета СЗГУ.

Процедура идентификации

Выбор структуры

мате матич е ско и модели,

основанный на имеющейся

априорной информации об

исследуемых процессах,

происходящих в СЗГРМ.

Выбор критерия близости модели и объекта исследования, основанный на специфике задачи.

Определение параметров модели, оптимального с точки зрения выбранного критерия близости и построение характеристик.

Рис. 2. Этапы процедуры идентификации

плавающая

Рис. 3. Принципиальная схема СГРМ со струйно-золотниковым гидроусилителем

Рис. 4. Струйная трубка: а - в изометрии, б — в разрезе

Рис. 6. Золотниковый каскад управления: а — золотник, б — гильза золотника

Рис. 7. Проточная часть струйного каскада: 1 — струйная трубка, 2 — камера распространения струи, 3 — приемная плата

Рис. 5. Приемная плата струйного каскада: а — в изометрии, б — в разрезе

Построение в реальном масштабе элементов СЗГУ в пакете прикладных программ Solid Works (рис. 4-9) позволяет рассчитать массу и объемы полостей элементов, необходимых для математического моделирования.

Выбор структуры модели позволяет определить перечень параметров, задаваемых при идентификации и перечень параметров, определяемых при идентификации (см. табл. 1), после чего производится расчет статических и динамических характеристик СЗГРМ.

Рис. 9. Перемещение струйной трубки

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

Математическая модель СГРМ со струйно-золотниковым гидроусилителем включает следующие уравнения.

Уравнение моментов, действующих на якорь ЭМП, определяется формулой [4]:

J

сРа dt2

к„

а

. da

(3)

Уравнение электрической цепи ЭМП:

da ' ~dt

da

Kpe'її; (4)

и = к„ • (ие - Кос • у(*)).

Условием баланса является равенство расходов на выходе из СГР и расхода в ЗГР

А,

dt

Wp 2 • Е

dPd

dt

- т

= <

, г < z„;

(5)

Pd

Zr,-

Для золотникового гидрораспределителя уравнение сил и баланса расходов выглядят следующим образом [4]:

(I2

мх • (г//2-г‘- «(О) = Лии • г5 • <1-1-ип -

(6)

Ье ' ( ^х:юл (t) J С -Ї • Х:тл (t)

Р'і рх ( і; і. " і (j)

л chi , f

‘ y dt \2-Ed) dt

= ' х:ю:л (t) ' '

2 • (Pp - Pe) P

(7)

Уравнение движения нагрузки [4]:

R

си,

У

ъ

" dt

(8)

В пакете прикладных программ Мар1е были построены переходные процессы СЗГРМ: перемещение струйной трубки, золотника и выходного звена — гидродвигателя и их скорость (рис. 9-14). Построен фазовый портрет рулевой машины, описывающий поведение всей системы и наглядно определяющий положение элементов гидроусилителя (рис. 15,16).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Рис. l0. Скорость перемещения трубки

струйной

Таблица 1

Параметры, задаваемые при идентификации Параметры, определяемые при идентификации

электромеханический преобразователь

входное напряжение [/, индуктивность управляющей обмотки Ь, активное сопротивление управляющей обмотки коэффициент усиления усилителя по напряжению Кц, коэффициент противо-ЭДС Кре, коэффициент главной обратной связи Кос, угол поворота якоря ЭМП о, коэффициент вязкого трения Ъ, момент инерции якоря ЭМП J, сигнал рассогласования [/,,, максимальный момент якоря ЭМПЛ/ПШХ; управляющий ток ЭМП I, коэффициент момент-ной характеристики Кгг,л, коэффициент жесткости магнитной пружины Кта;

гидрораспределитель «струйная трубка-приемная плата»

теоретический расход через струйную трубку С)т, длина струйной трубки /, плотность рабочей жидкости р, перемещение струйной трубки, при котором наступает насыщение характеристик СГРМ г„; диаметр сопла струйной трубки й, максимальное перемещение струйной трубки 2тах, коэффициент расхода в СГУ рд, диаметр окна приемной платы йп, расстояние между струйной трубкой и приемной платой Я, коэффициент восстановления давления ;

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

величина позиционной нагрузки на золотнике Сх, масса золотника Мх, сила вязкого трения Ье_, сила сухого трения РіГХ, коэффициент расхода золотника (1х, давление питания Рр; давление слива Рс; диаметр золотника (!зол; ширина окна золотника Ьок'; площадь золотника Лзол; гидродинамическая сила С*гд; модуль объемной упругости рабочей жидкости в золотнике Е;

исполнительный гидродвигатель (гидроцилиндр)

коэффициент жесткости стенок гидроцилиндра Сук; сила сухого трения в гидроцилиндре Р"Г; коэффициент вязкого трения поршня гидроцилиндра Ь„; масса поршня, гидроцилиндра и рабочей жидкости тп; скорости, развиваемые СГРМ под нагрузкой и на холостом ходу; приведенная масса нагрузки М; статическая нагрузка Я; позиционная нагрузка С„; коэффициент жесткости проводки Лг„; коэффициент вязкого трения нагрузки Ь; усилие на штоке гидроцилиндра Л; эффективная площадь гидродвигателя А; сила сухого трения в гидроцилиндре Р^; модуль объемной упругости рабочей жидкости в гидродвигателе Еа;

Основные параметры гидроусилителя следующие: I = 0,25 А, М = 20 кг, давление питания в системе МПа, давление

питания первого каскада управления

, расход в системе ,

расход в струйной трубке ,

.

Результаты численного моделирования показывают, что представленные характеристики отвечают основным требованиям по устойчивости, точности, управляемости и быстродействию, предъявляемым к исполнительным механизмам систем управления летательными аппаратами, при этом статическая ошибка лежит в пределах 3 %, динамическая ошибка в пределах 10 %, что позволяет судить о высокой адекватности математической модели реальному объекту [3].

Анализ переходных процессов показывает, что характеристики имеют зоны неустойчивости, что связано с неоднозначным определением коэффициентов уравнений. В подоб-

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

Уточнение математического описания рулевых приводов целесообразно проводить с использованием результатов экспериментальных исследований, основной задачей которых является определение так называемых «недоопределенных» параметров СГРМ. По существу, необходимо по некоторой совокупности параметров, измеряемых в процессе испытаний, воспроизвести недостающие характеристики СГРМ, установить функциональные зависимости между входными и выходными параметрами, а также получить данные, необходимые для идентификации объекта.

Рис. 11. Перемещение золотника

Рис. 15. Фазовый портрет гидроусилителя

Рис. 12. Скорость перемещения золотника

Рис. 13. Перемещение гидроцилиндра

Рис. 14. Скорость перемещения гидроцилиндра

Рис. 16. Фазовый портрет гидроусилителя во времени

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

Статистическая обработка результатов экспериментальных исследований проводится следующими этапами (рис. 17). Основному этапу статистической обработки — проверке на соответствие нормальному закону распределения предшествует теоретический (качественный) анализ имеющихся данных с учетом процессов и явлений, влияющих на разброс параметров.

Рис. 17. Основные этапы статистической обработки результатов эксперимента

Основная схема проведения статистического анализа заключается в следующем: полученную в результате эксперимента совокупность значений группируют и на основании вида распределения выдвигают гипотезу о конкретном теоретическом законе распределения. Далее, пользуясь критерием Пирсона с определенной степенью надежности, проверяют эту гипотезу и определяют параметры гипотетического закона распределения рассеяния отклонений. Затем, пользуясь тестом Бартлета, проверяют гипотезу о постоянстве дисперсий и определяют характеристики смещения центра распределения. Однако следует отметить, что на практике критерий Пирсона применяется достаточно редко, так как при выборке небольшого объема он не дает полной картины и приходится использовать другие критерии для оценки распределения [2].

Использование классического подхода к обработке данных основано на принятии гипотезы о соответствии распределения этих данных какому-либо закону распределения. Как правило, таким законом является нормальный закон распределения, поскольку при стремлении количества данных к бесконечности все известные законы распределения стремятся к нему. Специфика производства и эксплуатации СГРМ не позволяет получать большие объемы экспериментальных данных (выборка не превышает 10-20 значений), поэтому о строгом соответствии распределения опытных данных какому-либо закону говорить затруднительно. В связи с этим возникает необходимость разработки методов обработки, нечувствительных к отклонениям от классических статистических моделей, т. е. таких, чтобы малая ошибка в математической

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

Робастные оценки должны удовлетворять двум основным требованиям:

1) мало уступать в эффективности оптимальным оценкам при выполнении основной модели, когда распределение действительно является гауссовским;

2) оставаться достаточно «хорошими» при отклонениях от основной модели, когда распределение отлично от гауссовского.

Робастные оценки строятся таким образом, чтобы их свойства оставались хорошими даже в том случае, когда истинное распределение экспериментальных данных отличается от предполагаемого. Так, для многих наборов данных более реалистичной является не строгая модель гауссовского распределения, а распределения с более «тяжелыми» хвостами (когда вероятность больших отклонений от среднего существенно больше, чем при гауссовском распределении).

Обработка, обобщение и анализ экспериментальных и теоретических исследований СГР, проводимых совместно ФГУП «Государственный ракетный центр — КБ им. акад. В. П. Макеева» позволили вывести эмпирические зависимости по расчету ряда стохастических коэффициентов, входящих в математическую модель рулевого привода, т. е. коэффициентов, имеющих случайный разброс в силу действия различных факторов.

Природа процессов и явлений, происходящих в проточной части СГРМ, описывается целым рядом сложных процессов, сопровождающих распространение высоконапорной скоростной струи в камере ограниченных размеров. В какой-то степени картину течения можно охарактеризовать коэффициентом расхода — отношением максимального расхода через гидродвигатель к расходу через струйную трубку. Он зависит от числа Рейнольдса Ие, отношения диаметра приемного сопла к диаметру нагнетающего сопла с!,, отношения диаметра приемного сопла к расстоянию между соплами Н, давлений питания и слива.

На рис. 18. представлена зависимость коэффициента расхода от геометрических параметров СГР, рассчитанная по выражению

0,769 - 0,023 х

__ _____ и

Н -4= (в \ л

х = + 0,329 -Н^ + 0,203- = .

(Р ' ' \н)

(9)

Рис. 18. Коэффициент расхода в СГУ

Погрешность между результатами расчета по эмпирической зависимости и экспериментами не превышает в среднем 3%. Эмпирическая зависимость коэффициента расхода от безразмерных конструктивных параметров позволяет оптимизировать выбор наиболее рациональных геометрических параметров СГР.

Экспериментальные данные ФГУП «Государственный ракетный центр — КБ им. акад. В. П. Макеева» по коэффициенту расхода представлены на рис. 19.

Рис. 19. Экспериментальные данные коэффициента расхода

Проведем статистическую обработку результатов экспериментальных исследований. Совокупность всех испытаний и результатов наблюдений, которые могут быть получены в данных условиях, называется генеральной совокупностью, а результаты, полученные при испытаниях, — выборкой из генеральной совокупности. С увеличением объема выборки плотность вероятностей оценки приближается к истинному значению математического ожидания , т. е. дисперсия оценки с ростом объема выборки стремится к нулю (рис. 20) [1]. В реальных условиях вы-

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

т

Рис. 20. Зависимость дисперсии оценки от объема выборки

Первичную обработку результатов измерений выполняют методами математической статистики, в результате чего получают статистические оценки (математическое ожидание, дисперсия, среднеквадратичное отклонение, плотность распределения вероятности и т. д.). Для оценки совпадения (согласия) эмпирических и теоретических статистик применяют различные критерии согласия. Наиболее распространенными и эффективными являются статистические критерии Пирсона, Колмогорова и Стьюдента [1,2].

Известные труды теории вероятности позволяют определить на основании экспериментальных данных математическое ожидание , дисперсию величины, среднее квадратичное отклонение , плотность распределения вероятности, с помощью функции Лапласа при заданной вероятности события определяются границы диапазона изменения, и характеристики СГРМ представляют уже не кривые, а области значений [1,2].

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

График плотности распределения вероятности коэффициента расхода в соответствии с нормальным законом распределения представлен на рис. 21.

Рис. 21. Плотность распределения вероятности коэффициента расхода и экспериментальные данные

В табл. 2 представлены основные методы проверки экспериментальных данных нормальному распределению.

Для оперативной оценки гипотезы нормальности распределения применялись следующие методы: проверка с помощью среднего абсолютного отклонения, проверка с помощью размаха варьирования, критерий, основанный на сравнении выборочных моментов с генеральными параметрами, критерий Пирсона и т. д. Использование различных критериев и экспресс-оценок не позволяет однозначно определить на примере одной конкретной выборки результатов экспериментальных исследований возможность и правомерность использования закона нормального распределения при статистической обработке.

Таблица 2

Наименование метода Выборка Основные уравнения Обработка коэффициента расхода при п = 15 Вывод

1 2 3 4

1 Проверка с помощью среднего абсолютного отклонения (С АО) С АО = где: п — объем выборки; х — среднее выборки. ^ ' Е (*-■ - *)2 £ х = ' г, ■ Гипотеза нормальности подтверждается, если - 0,79791 (^= С АО = 0,593 • 10_1) 5 = 0,267 0,7979 < 0,1633

Выборка: п >120

Вывод: условие нормальности распределения не подтверждается

2 Быстрая проверка гипотезы нормальности распределения осуществляется с помощью размаха варьирования Я Из экспериментальных данных выбираются максимальная и минимальная величины и вычисляется размах варьирования данной выборки И :=:: |я?тах Xтт | • среднее значение в выборке и среднеквадратичное отклонение Е ач 1 п х = ^; 5= Л-Е(х,--х)2. Гипотеза нормальности распределения подтверждается при [#],<[#]<[#]«> где нижнее и верхнее граничные значения берутся из таблицы [1] в зависимости от величины выборки п и уровня значимости р. Л = 0,1091, 5 = 0,2582, 0,43 < 0,4224 < 0,45

Вывод: условие выполняется не полностью или экспериментальные данные не подчиняются закону нормального распределения

Выборка 3 < п < 10000

Окончание табл. 2

1 2 3 4

3 Выборка 3 < п < 10000 Критерий, основанный на сравнении выборочных моментов с генеральными параметрами Вычисляются оценки математического ожидания, дисперсии и среднеквадратичные отклонения, выборочные коэффициенты эксцесса и ассиметрии п 71 = ^ ' Е (х; - р):1; 1=1 7-> = ^4 • Е (х; - Д)'4 - 3. 1=1 Распределения этих оценок сложны и малоизученны, обычно известны дисперсии этих величин -V. \ _ 6- (п — 1 ) ° П1] — („ + !).(„+3) ; -IV.. ч _ 24-(п-2)-(п-:Я (п + 1)2-(п+а)-(п+3)' Наблюдаемое распределение можно считать нормальным при |711 < 3 • \/сг-(')'!), 72 < 5 • \/а-(у2). р, = 0,87044, Б = 0,0011, а,, = = 0,033, 71 = —0,8916, 71 = —0,2255, а2( 71) = 0,54006, <г2(у2) = 0,20156, 0,8916 < 0,66317 0,2255 < 0,5269.

Выборка п > 20 Вывод: условие нормальности не подтверждается

4 Критерий согласия Пирсона Рассчитываются эмпирические частоты у = где р\ — вероятность выпадания данного значения, п — количество данных. Рассчитывается вероятность событий по правилу трех сигм, а также теоретические частоты и; = рг ■ Определяются к — число степеней свободы, г — число параметров предполагаемого распределения. По таблице для заданного уровня значимости и при известном числе степеней свободы определяется значение г’. Определяются значения П П о = Е ^ , х, = Е £ 1=1 1=1 Гипотеза принимается при х„ < \2 и прих; < к = 2, г = 12 р = 0.9973 Г’ = 21 х„ = 0,828 • 10_“х; = 13,864 0,8284 • 10-"' < 21, 13,864 < 21

Вывод: условие нормальности не подтверждается

Поэтому с помощью численного моделирования на ЭВМ удалось сравнить статические характеристики СГРМ с классической оценкой — средним арифметическим и характеристиками, построенными с помощью робастных методов, основанных на медианной оценке.

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

ных методов, основанных на медианной оценке (см. рис. 22).

Завершающим этапом идентификации СЗГРМ является расчет характеристик стохастической СЗГРМ

На первом уровне идентификации динамических характеристик СГРМ учитывается случайный разброс одного параметра — коэффициента расхода в СГУ. На следующих этапах разрабатываются модели, учитывающие случайный разброс двух стохастических параметров — коэффициента расхода и коэффициента восстановления давления в СГУ.

В результате построения стохастических математических моделей, характери-

стики рулевой машины представляют собой не кривые, а диапазон изменения значений с заданной степенью вероятности событий (см. рис. 23, 24).

стики рулевой машины представляют собой не кривые, а диапазон изменения значений с заданной степенью вероятности событий (см. рис. 23, 24).

Рис. 22. Сравнение обобщенных характеристик с экспериментальными данными

Рис. 23. Перемещение поршня гидроцилиндра

Рис. 24. Скорость перемещения струйной трубки

В результате построения стохастических математических моделей, характери-

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

• ряд параметров СЗГРМ не поддавался аналитическому описанию из-за отсутствия теоретических исследований в области проектирования высоконапорных рулевых машин, а также из-за сложных гидродинамических процессов, протекающих в проточной части рулевой машины. Поэтому в работе был использован метод идентификации, при котором математическая модель строится по результатам экспериментов или данным;

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

выбраны параметры, позволяющие рассчитать математическую модель и построить переходные процессы рулевой машины. Анализ переходных процессов показывает, что характеристики имеют зоны неустойчивости, что связано с неоднозначным определением коэффициентов уравнений. В подобной ситуации необходимо провести дополнительную статистическую обработку экспериментальных данных для доопределения неоднозначных коэффициентов;

выбран оптимальный метод обработки экспериментальных данных и проведена статистическая обработка результатов экспериментальных исследований по расчету стахо-стических коэффициентов;

разработана стохастическая математическая модель, учитывающая разброс параметров СЗГРМ с заданной степенью вероятности события.

СПИСОК ЛИТЕРАТУРЫ

ОБ АВТОРАХ

1. Валеева, О. В. Низкотемпературные твердотопливные газогенераторы: методы расчета рабочих процессов, экспериментальные исследования / О. В. Валеева, С. Д. Ваулин [и др.]. Миасс : ГРЦ «КБ им. академика В. П. Макеева», 1997. 268 с.

2. Гмурман, В. Е. Теория вероятностей и математическая статистика / В. Е. Гмурман. М. : Высшая школа, 1977. 479 с.

3. Кириллов, Ю. К. Струйные гидравлические рулевые машины (теоретические основы рабочих процессов) / Ю. К. Кириллов, А. М. Русак, Ю. Н. Скорынин [и др.]. Уфа : УГАТУ, 2002. 284 с.

4. Месропян, А. В. Численное моделирование характеристик струйно-золотниковых гидроусилителей / А. В. Месропян, К. А. Широкова // Гидропневмоавтоматика и гидропривод-2005 : сб. науч. тр. Т. 1. Ковров : КГТА, 2006. 326 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

5. Широкова, К. А. Использование статистических методов при проектировании струйных гидравлических рулевых машин / К. А. Широкова // Вопросы теории и расчета рабочих процессов тепловых двигателей : межвуз. сб. науч. статей. Уфа : УГАТУ, 2006. №21. С. 171-177.

Месропян Арсен Владимирович, доцент каф. прикл. гидромеханики. Дипл. инж.-мех. (УГАТУ, 1996). Канд. техн. наук по системам приводов (Пермск. гос. тех. ун-т, 2000). Иссл. в обл. гидроприводов систем управления ЛА.

Широкова Ксения Александровна, асп. той же каф. Дипл. бакалавр (УГАТУ, 2004) и магистр (там же, 2006) по энергомашиностроению. Готовит дис. в обл. систем автоматики двиг. установок.

Целищев Владимир Александрович, проф. той же каф. Дипл. инж.-мех. по гидравл. машинам (УГАТУ, 1982). Д-р техн. наук по тепловым двигателям (УГАТУ, 2000). Иссл. в обл. систем автоматики ЛА и двигательных установок.

i Надоели баннеры? Вы всегда можете отключить рекламу.