2.5 Анализ динамической модели гидросистемы
Динамическая модель описывает переходный процесс гидросистемы. В общем случае система дифференциальных уравнений, описывающих гидравлическую систему, имеет вид:
(48)
где А – матрица Якоби,
- вектор фазовых координат,
- вектор функции внешних воздействий.
С учетом произведенных ранее расчетов, запишем систему дифференциальных уравнений, представляющую динамическую гидросистему:
(49)
Для динамической модели матрицу Якоби можно записать аналогично статической модели:
(50)
Переходный процесс определяется в результате численного интегрирования системы (49), для чего необходимо произвести выбор ряда параметров.
Пусть переходный процесс оценивается как реакция системы, находящейся в состоянии покоя, на ступенчатое воздействие вида:
(51)
где u0 и uk – начальное и конечное значение функции воздействия u(t), причем u0 и uk – константы, (u0 ≠ uk):
(52)
Начальные (46) и конечные (47) значения всех фазовых координат определены при анализе статического режима (таблица 5).
=> (53)
Если система устойчивая, то через некоторый промежуток времени, система перейдет из состояния V0 в состояние Vk. Для численного интегрирования будем использовать неявный метод Эйлера. Вектор входных воздействий Vk при имеет вид:
(54)
Решим систему дифференциальных уравнений методом Эйлера. Формула численного интегрирования неявного метода Эйлера имеет вид:
(61)
Совместное преобразование двух последних выражений приводит к записи:
(62)
где - модифицированная матрица Якоби на k+1 шаге, которая формируется по следующему правилу: диагональные элементы матрицы Якоби на k-ом шаге пересчитываются по формуле:
(63)
Остальные элементы не изменяются. Для матрицы размерности 5х5 получаем:
(64)
- модифицированный вектор входных воздействий на k+1 шаге, определяемый по формуле:
(65)
Решение системы уравнений (61) дает значение фазовых координат на k+1 шаге, то есть в момент времени tk+1.
Алгоритм неявного метода Эйлера с постоянным шагом интегрирования:
-
задание шага интегрирования h (принимаем h = 0,5);
-
задание начальных значений фазовых координат при t0 = 0;
-
вычисление времени tk+1=tk+h, где k=0, 1, 2…;
-
вычисление модифицированных матриц и на k+1 шаге;
-
решение системы уравнений (61) с целью определения в момент времени tk+1;
-
переход к этапу (3) до тех пор, пока в случае устойчивой системы фазовые координаты не достигнут состояния конечного значения .
Начальные значения вектора определяются на основании входных воздействий системы. В качестве начальных значений фазовых переменных берем вектор начальных значений .
Рисунок 7 – Графики фазовых координат f(n)0, f(n)1, f(n)2, f(n)3
Рисунок 7 – Переходный процесс гидросистемы