Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

FORMATION OF VIBRATION IMPACTS TAKING INTO ACCOUNT RANDOMNESS IN THE PARAMETERS OF ELECTRO-HYDRAULIC ACTUATOR OF VIBRATION TEST STAND

Makarenkov A.M. 1 Tun Tun Kyaw 1 Tin Aye Kyaw 1
1 Bauman Moscow State Technical University (Kaluga Branch)
This article is devoted the problem of forming a random vibration impacts using an electro-hydraulic vibration test stand. The electro-hydraulic servo actuator is used as a drive of vibration test bench. The electro-hydraulic actuator of this type is characterized by a random instability of physical parameters, therefore it can introduce considerable distortion in the statistical characteristics of the generated vibration impacts. If the input signal of vibration test stand is a mixture of deterministic and random signals, the randomness of parameters of an electro-hydraulic actuator leads to significant fluctuations in the variance of generated vibration impacts, since the correlation function of the output signal of a stochastic system is determined not only by the correlation function of the input signal, as is the case for deterministic systems, but its mathematical expectation. The method is proposed for compensation of the randomness in the parameters of the electro-hydraulic servo drive, the input of which receives the amount of random and deterministic signals. This method is based on the proportional change in the variance of the random component of the input signal, which compensates for the influence of its deterministic component. The desired law of change of the mentioned variance is determined as a result of minimization of a quadratic functional, calculated using the averaged projective model of a stochastic system. A numerical example is given demonstrating the effectiveness of the proposed approach.
vibration test stand
random parameters
stochastic system
electro-hydraulic servo drive
mathematical model
projective model
1. Baranov V.N., Zakharov Yu.E. Elektrogidravlicheskie i gidravlicheskie vibratsionnye mekhanizmy [Electro-hydraulic and hydraulic vibration mechanisms]. Moscow, Mashinostroenie, 1977. 326 p.
2. Lapin S.V., Egupov N.D. Teoriya matrichnykh operatorov i ee prilozhenie k zadacham avtomaticheskogo upravleniya [The theory of matrix operators and its application to problems of automatic control]. Moscow, Izd-vo MGTU im. N.E. Baumana, 1997. 496 p.
3. Metody inzhenernogo sinteza slozhnykh sistem upravleniya: analiticheskii apparat, algoritmy prilozheniya v tekhnike [Methods of engineering synthesis of complex control systems: analytical apparatus, algorithms, applications in engineering]. Ed. by K.A. Pupkov and N.D. Egupov. Moscow, Izd–vo MGTU im. N.E.Baumana, 2012.
4. Popov D.N. Dinamika i regulirovanie gidro- i pnevmosistem [Dynamics and regulation of hydraulic and pneumatic systems]. Moscow, Mashinostroenie, 1977. 424 p.
5. Pupkov K.A., Egupov N.D., Makarenkov A.M., Trofimov A.I. Teoriya i kompyuternye metody issledovaniya stokhasticheskikh system [Theory and computer methods of investigation of stochastic systems]. Moscow, Fizmatlit, 2003. 400 p.

Повышение качества и надежности выпускаемых изделий приборо- и машиностроения обеспечивается в числе прочего механическими испытаниями на вибрационную нагрузку, которые выполняются как для разрабатываемых новых изделий, так и для серийных образцов. При вибрационных испытаниях объектов, имеющих большую массу, требуется достаточно мощный силовой привод, в качестве которого в виброиспытательных стендах часто используется электрогидравлический следящий привод (ЭГСП), отличающийся высокой удельной мощностью и быстродействием [1], [4]. При испытаниях на случайную вибрационную нагрузку необходимо получить на выходе виброиспытательного стенда испытательное воздействие (перемещение рабочего стола или платформы стенда с закрепленным изделием) в виде случайного процесса с требуемой корреляционной функцией и математическим ожиданием. Однако ввиду нестабильности параметров ЭГСП, обусловленных как внешними (температура окружающей среды), так и внутренними (технологический разброс значений параметров, содержание газовой фазы в рабочей жидкости) случайными факторами, имеет место искажение упомянутых статистических характеристик формируемого вибрационного воздействия. Например, если на вход виброиспытательного стенда подается аддитивная смесь детерминированного и случайного сигналов, то случайность параметров виброиспытательного стенда приводит к существенным колебаниям дисперсии формируемого случайного вибрационного воздействия, вызванным колебаниями детерминированной составляющей входного сигнала. Это является одной из особенностей стохастических систем, состоящей в том, что корреляционная функция, а значит, и дисперсия выходного сигнала системы со случайными параметрами определяется не только корреляционной функцией входного сигнала, как это имеет место для детерминированной системы, но и математическим ожиданием данного сигнала. Таким образом, указанные колебания дисперсии выходного сигнала ЭГСП, искажающие характеристики вибрационного воздействия, могут привести к снижению качества вибрационных испытаний и получению неадекватных результатов в случаях, когда требуется обеспечить точное соответствие формируемых вибрационных воздействий некоторым реальным. Следовательно, задача компенсации фактора случайности параметров ЭГСП является актуальной.

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

Компенсация фактора случайности параметров

Для компенсации фактора случайности параметров ЭГСП предлагается использовать подход, состоящий в формировании такого случайного входного сигнала, который обеспечил бы устранение указанного влияния детерминированной составляющей (математического ожидания) входного сигнала за счет соответствующей модуляции амплитуды его случайной составляющей. В терминах теории управления такой подход соответствует задаче синтеза программных управлений в статистической постановке и относится к классу обратных задач динамики [3]. Далее рассмотрим алгоритм синтеза программного управления механическими колебаниями штока гидроцилиндра ЭГСП при формировании случайного вибрационного воздействия, основанный на использовании усредненной проекционной модели и процедуры оптимизации квадратичного функционала.

Математическая модель ЭГСП в линейном приближении описывается следующим дифференциальным уравнением:

mak03.wmf, (1)

где y(t) – входной сигнал (электрическое напряжение), x(t) – выходной сигнал (перемещение штока поршня).

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

mak04.wmf, (2)

mak05.wmf, (3)

где A – квадратная матрица проекционной характеристики (матричный оператор) системы; mak06.wmf и mak07.wmf – квадратные матрицы проекционных характеристик корреляционных функций входного и выходного сигналов; mak08.wmf и mak09.wmf – вектор-столбцы проекционных характеристик функций математических ожиданий выходного и входного сигналов; mak10.wmf – оператор математического ожидания; T – знак транспонирования.

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

mak12.wmf,

mak13.wmf

mak14.wmf,

mak15.wmf.

Проекционная характеристика системы A строится с использованием матричных операторов интегрирования [2]. В выражение для нее входят коэффициенты уравнения (1), поэтому матрица A является случайной. Операция усреднения mak16.wmf реализуется с использованием приема разложения стохастического оператора A в матричный ряд с последующим аналитическим усреднением членов этого ряда [5]. Аналитическое усреднение состоит в выражении стохастических моментов (порядка выше второго) случайных коэффициентов исходного уравнения модели (1) через дисперсии и математические ожидания этих коэффициентов для каждого члена упомянутого ряда.

Операторная форма записи решений, характерная для проекционных методов, позволяет предложить следующий подход к компенсации фактора случайности параметров ЭГСП на основе анализа выражения (2). Из данного выражения видно, что влияние ненулевого математического ожидания входного сигнала обусловлено вторым и третьим слагаемым в его правой части. Идея компенсации этого влияния состоит во введении матрицы умножителя на некоторую корректирующую функцию f(t) (матричный оператор умножения Af) в первое слагаемое, то есть в преобразовании (2) к виду

mak17a.wmf

mak17b.wmf, (4)

что соответствует введению умножителя в цепь центрированной случайной составляющей входного сигнала системы. Далее скорректированная случайная и детерминированная составляющие поступают на вход сумматора, формирующего смешанный входной сигнал y(t).

Уменьшение среднего уровня дисперсии выходного сигнала системы, обусловленное действием умножителя Af, компенсируется добавлением к входному сигналу второй случайной составляющей mak18.wmf в виде центрированного случайного процесса, некоррелированного с первой случайной составляющей mak19.wmf и имеющего корреляционную функцию того же вида, но с дополнительным коэффициентом k. Например, если корреляционная функция случайной составляющей входного сигнала имеет вид mak20.wmf, то корреляционная функция второй случайной составляющей mak21.wmf будет определяться как

mak22.wmf. (5)

Тогда выражение (4) будет иметь вид

mak23a.wmf

mak23b.wmf, (6)

а сигнал на входе ЭГСП будет формироваться следующим образом:

mak24.wmf. (7)

Если детерминированная составляющая является гармоническим сигналом mak25.wmf, то функция f(t) может выглядеть так:

mak26.wmf, (8)

где коэффициент r определяет уровень компенсации, а ? – фазовый сдвиг, позволяющий получить «противофазный» сигнал f(t), гасящий колебания дисперсии выходного сигнала системы.

Для нахождения параметров r, ? и k воспользуемся приемом минимизации следующего функционала, выражающего близость корреляционных функций выходного сигнала стохастической системы (1) и детерминированной системы, также описываемой уравнением (1), в котором все случайные коэффициенты заменены своими математическими ожиданиями:

mak27.wmf, (9)

где mak28.wmf и mak29.wmf – корреляционные функции выходного сигнала детерминированной и стохастической системы соответственно, при этом последняя зависит от искомых параметров компенсации.

Проекционная аппроксимация модели (1) позволяет представить функционал (9) в следующей удобной для вычисления форме:

mak30.wmf, (10)

где mak31.wmf – элементы квадратной матрицы, вычисляемой как

mak32.wmf, (11)

где mak33.wmf и mak34.wmf – проекционные характеристики корреляционных функций mak35.wmf и mak36.wmf соответственно.

Проекционная характеристика mak37.wmf (r, ?, k) в (11) вычисляется по усредненной проекционной модели системы (1).

Пример компенсации случайности параметров электрогидравлического привода

В качестве примера решения задачи компенсации влияния случайности параметров ЭГСП, описываемого моделью (1), где n = 5 и m = 0, построенной на основе типовой модели ЭГСП [4], на дисперсию его выходного сигнала, рассмотрим отработку детерминированного сигнала mak38.wmf с наложенной на него центрированной случайной составляющей mak39.wmf с корреляционной функцией mak40.wmf. В качестве случайного параметра возьмем коэффициент вязкого трения на золотнике электрогидравлического усилителя, обозначив его как hЗ. Данный параметр будем считать гауссовой случайной величиной с математическим ожиданием mak41.wmfи дисперсией mak42.wmf. При этом в уравнении (1) случайными будут коэффициенты а1, а2, а3, а4, математические ожидания и дисперсии которых вычисляются через заданные mak43.wmf и mak44.wmf. В качестве ортогонального базиса Ф(t) будем использовать базис функций Уолша с l = 64.

makar1.wmf

Рис. 1. Дисперсия выходного сигнала ЭГСП: 1 – для детерминированной модели; 2 – для стохастической модели без компенсации влияния случайности параметра; 3 – для стохастической модели c компенсацией влияния случайности параметра

makar2.tif

Рис. 2. Корреляционная функция входного сигнала ЭГСП

Выполняя минимизацию функционала (10) методом Нелдера-Мида при начальных значениях r = 0,5, ? = 0,3, k = 0,5, находим следующие оптимальные значения параметров компенсации: r = 0,5267, ? = 0,2213, k = 0,4664.

На рис. 1 представлены графики дисперсии выходного сигнала ЭГСП x(t) до и после компенсации влияния случайности параметра hЗ.

На рис. 2 приводится график корреляционной функции скорректированного входного сигнала y(t), формируемого согласно (7).

Заключение

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и Правительства Калужской области (грант № 14-41-03071).