Расчетная схема и математическая модель робота
Расчетная схема робота показана на рис. 1. Рассматриваемый в работе робот движется в неподвижной системе координат oxyz , причем предполагается, что робот находится все время на плоскости Oxy с координатой z=0.
Робот состоит из ведущего колеса 1, установленного на валу управляемого электродвигателя постоянного тока, который неподвижно закреплен на корпусе робота 2, и поддерживающего рояльного колеса 3. В дальнейшем принято, что элементы робота - недеформируемые твердые тела, а при взаимодействии ведущего колеса с шероховатой поерхностью возникает нормальная реакция, сила трения и момент трения качения. На рояльное колесо со стороны поверхности действует только нормальная реакция.
Рис.1 |
Пусть вес корпуса робота с электроприводом и рояльным колесом приложен в точке О2 .А вес ведущего колеса приложен в точке О1 , являющейся геометрическим центром колеса.
Для построения математической модели, описывающей движение робота, рассмотрим отдельно движение корпуса совместно с рояльным колесом и ведущего колеса. Пусть x2,y2, ф2 - обобщенные координаты, определяющие положение корпуса робота. Тогда дифференциальные уравнения, описывающие движение корпуса робота, можно представить в виде
M2x2 = RX2
M2 У2 =-m2g + N2 + R12
R1 |
Тродвигателя; m1 - масса колеса ф1 - |
J02 ф2 = МЭД + N2(l1 -12)-RX2(l3 -Г) + R12l2
Где m2,J02 - масса и момент инерции корпуса робота; m1- масса ведущего ко-
Леса; Rf2
12 - проекции силы, действующей со стороны колеса на корпус, N2 - нормальная реакция; l1, 12 - геометрические размеры корпуса; r - радиус колеса; Мэд - момент, действующий на корпус со стороны электродвигателя.
Пусть x1, y1 , ф1 - обобщенные координаты, определяющие положение ведущего колеса.
Дифференциальные уравнения движения колеса робота с учетом свойств управляемого электропривода примут вид
Ld + iR + сю (Р1 = U(t) dt
M1 Х1 = RX1 + F-П m1 У1 = Ry1 + N1
J01 ф1 =- N15 + Мэд + ф.
При этом предполагаем, что ф1 не
Равно нулю.
Здесь обозначено: L, R - индуктивность и активное сопротивление якоря электродвигателя; сю - постоянная элек - угловая скорость колеса; N1,Fj!^- нор-
Мальная реакция и сила трения, действующие на колесо со стороны опорной поверхности; Ry1, R^ - силы, действующие на колесо со стороны корпуса; x1, y1 , ф1 - обобщенные координаты, определяющие положение колеса; 5-коэф - фициент трения качения.
Будем считать, что момент электродвигателя зависит от тока якоря по следующему закону:
Myi = Mya (i) =
= M0ya + i |
RX _ nx 12 = - R21 • RX — и x 12 = - R21. |
.dM^ (8)
+ i3^^^^)3 +... di di
Ограничим этот ряд только вторым и третьим слагаемыми, то есть
Mya = 1 • C1e +13 • C3e,
Где^ C1e, C3e - экспериментально определяемые постоянные.
На первом этапе будем считать, что проскальзывание ведущего колеса по шероховатой поверхности отсутствует, кроме этого, колесо и опорная точка не отрываются от поверхности, то есть N1 > 0
N>0, тогда можно записать следующие уравнения связей:
TOC o "1-3" h z ф1г = x1. (9)
У1 =r. (10)
X2 = x1 +12 , У2 = l3 , ф2 = (11)
Запишем также очевидные равенства:
(12) (13)
Анализируя уравнения (1),(5),(7) в
Начальный момент времени, когда ф1 = 0 , с учетом условий (9)-(13) запишем следующее выражение, определяющее коэффициент трения качения в покое: 51 = Mya / N1.
Если коэффициент трения качения 51 меньше предельного 5 , то есть
51 = Mya / N1, (14)
Вращение колеса не происходит, и робот стоит на месте.
Только при выполнении равенства 51 = 5 (15)
Начинается движение.
Задача исследования динамических закономерностей пускового режима робота может быть сформулирована следующим образом. Пусть дано: m1, m2, l1, І2, І3, r, 5, J01Z, J02Z. L, R,croc1e2c3e.
Решая уравнения (1)-(7), с учетом условий (8)-(15) необходимо определить
Неизвестные: x1, y1, x2, y2, ф1, ф2, N1, N1, FsiE), Rx2, R|1, R12, Ry1, i, Mya для различных законов управляющего напряжения U(t).
Рассмотрим различные способы управляемого пускового режима. На рис.2 приведен график зависимости напряжения питания, поступающего на обмотки якоря электродвигателя, от времени [6-9].
Кусочно-постоянный 2-этапный алгоритм управления.
Данный алгоритм характеризуется 4- мя параметрами U1, U2, t, t2, изменение которых влияет на разгон робота. Преоб-
U, B
Разуем систему уравнений (1)-(7) с учетом условий (8)-(15) к виду
Ld + iR + сшср1 = U(t) dt
(m1 + m1)g = N1 + N2
J01 ф1 =-N15+ ic1e -
-(m1 + m1)r2 ф1 0 = ic1e + N2 (I1 -12) -
-m2r ф1 (13 - r) - (N1 - m1g)12.
Система уравнений (16)-(19) имеет 4 неизвестных ф1, N1, N2,i
Из (18) найдем нормальную реакцию N1:
N1 = (ic1e - (m1 + m2)r2 ф1- J01 ф1)/-5
Подставив это выражение в (17), получим
(m1 + m1)g =
= (ic1e - (m1 + m2)r2 ф1- J01 ф1)/ 5 + N2.
Отсюда определим величину реакции N2 в виде
N2 = (m1 + m1)g -
(iC1e - (m1 + m2)r2 ф1 - J01 ф1)
Из уравнения (19) найдем:
M2r2 ф1(13 - r) =
= ic1e + N2(11 -12) - (N1 - m1g)12.
U2
U1
T, c
Рис.2. Зависимость управляющего напряжения от времени
Ток для нахождения момента электродвигателя найдем из уравнения полного тока:
Ld + iR + cffl(Pl = U(t). dt
Таким образом, уравнения имеют три неизвестных: (j, N2,i.
Подставляя переменные, получим следующее:
[J02 + (m1 + m2)r2] ф = 1 (21) = [-(m1 + m2)g + - (m2gl2 - ic1e)]5 + ic1e. i1
TOC o "1-3" h z Ld + iR + cro((1 = U(t). (22)
Dt
Обозначим общую массу робота m: m1 + m2 = m, (23)
А J - приведенный момент инерции:
J = J01 + (m1 + m2)r2. (24
Тогда уравнение (22) принимает вид
J(1 = [m2gг - mg - —]§ + ic1e 11 i1
Или
I * S
J(1 = (m2gf - mg)5 + i(c1e - ic1e-). (25)
II 11
(26) |
Величина коэффициента трения качения - значительно меньше длины 12, тогда приближенно можно считать, что
- = в« 0. 1
Примем следующие обозначения:
12
A = M2Gf - mg. 11
Тогда с учетом (27) уравнение (25)
Запишем в виде
JcP1 = A-+i(c1e - c1e-);
11
(20) |
Ldi + iR + cro((1 = U(t). dt
Система дифференциальных уравнений (28)-(29) содержит две неизвестные величины ф1, i, определив которые можно найти зависимости всех параметров определяющих пусковой режим робота. Задавая различные законы управляющего напряжения, можно изучить закономерности пускового режима робота при отсутствии проскальзывания для любого закона управляющего напряжения. Полученные формулы позволяют изучать процесс запуска колесного робота из любых начальных условий, в том числе и нулевых, для различных стратегий управления. Наиболее удобной средой для выполнения вычисления является вычислительный пакет MATLAB SIMULINK.