⎜⎜
⎜⎝⎜⎜ ⎟
Заметим, что в разобранном примере математическая модель, представленная системой дифференциальных уравнений (8), адекватна физическому объекту только в пределах деформаций, удовлетворяющих закону Гука, т.е. ес-
ли фазовые координаты x1,x2,x3 достаточно малы по абсолютной величине.
К дифференциальным уравнениям вида (2) можно прийти и в результате линеаризации исходных нелинейных дифференциальных уравнений движения объекта. Опишем процедуру линеаризации.
Пусть математической моделью управляемого динамического объекта служит система нелинейных дифференциальных уравнений
y =Y t y v( , , ) , t ∈[t T0, ], y ∈ Rn, v ∈ Rr . (9) Относительно функции Y :[t0 ,T]× Rn+r → Rn предполагается существование непрерывных частных производных не ниже второго порядка включительно по каждому из аргументов.
Допустим, что некоторой функции v∗ :[t T0, ]→ Rr отвечает решение y∗()⋅ = y∗(⋅,t0 , y0 ,v∗()⋅ ) дифференциального уравнения (9), удовлетворяющее начальному условию y t( 0)= y0 . Предположим, что именно эта функция y∗ (⋅) является требуемым законом движения для управляемого объекта. Однако при физической реализации указанного управления v∗ (⋅) закон движения y()⋅ реального динамического объекта вследствие ряда факторов (неадекватность математической модели, наличие неконтролируемых возмущений, невозможность в точности удовлетворить начальным условиям и др.) будет отличаться от идеального движения y∗ ()⋅ . Для реализаций управляющих воздействий и отвечающих им движений примем следующее представление:
Обозначая
A( )t =
∂y ∂v
и учитывая, что
y∗( )t = Y (t, y∗(t), v∗(t)), t ∈ [t0 ,T ],
приходим к уравнениям (2), в которых C t( )= 0, t ∈[t T0, ].
Пример 3*. На горизонтальный плоскости находится двухзвенный механический манипулятор, каждое звено которого представляет собой абсолют-
li ,i =1,2. Первое звено соединено с неподвижным основанием манипулятора вращательной парой O1 , а со вторым звеном – враща-
Рис. 3 тельной парой O2 . Масса схвата
манипулятора – m , центр масс i-го звена находится в середине стержня – точке Ci , его масса – mi , момент инерции i-го звена относительно своего центра масс – Ii ,i =1,2. В соединительных парах могут развиваться управляющие вращательные моменты, соответственно, v1 и v2 ,
На горизонтальной плоскости, в которой расположен манипулятор, введем прямолинейную ось O1 x. Обозначим через ϕi угол, образованный i-м звеном манипулятора, i =1,2, с осью O1 x. Запишем дифференциальные уравнения движения манипулятора в форме уравнений Лагранжа второго рода, в которых в качестве обобщенных координат берутся углы ϕi ,i =1,2 . Кинетическая энергия манипулятора определяется по формуле
T = T1 + T2 + Tc , (13)
где Ti – кинетическая энергия i-го, i =1,2, звена, а Tc – кинетическая энергия схвата манипулятора. Последовательно вычисляем
T1 = I1ϕ1 + m1vC =4I1ϕ1 + m1l1ϕ1 = ϕ1 ,
28
T2 =
Tc =
Подставляя найденные величины энергий составных частей манипулятора в (13), находим
8 8
Введем обозначения
4 4
Тогда выражение для кинетической энергии манипулятора принимает вид
T =
Справедливы равенства
∂T d ∂T
∂ϕ1 dt ∂ϕ1
∂T d ∂T
∂ϕ2 dt ∂ϕ2