Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Мехатроника и динамика мини-роботов

.pdf
Скачиваний:
0
Добавлен:
24.11.2025
Размер:
4.17 Mб
Скачать

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

Приближенное описание динамики сплошного деформируемого твердого тела на основе многомассовой системы может быть проведено с помощью лагранжевого формализма. Так как упругие деформации обычно малы, то линеаризованные уравнения движения представляются в виде (3.7), что для частотного анализа позволяет применять метод интегральных преобразований Лапласа. Это позволяет использовать методы анализа и синтеза систем, основанные на теории цепей и передаточных функций.

Применяя к уравнениям (3.7) преобразование Лапласа

p t e pt dt,

0

получим

 

Ap2

Bp c

 

p

 

 

p , p

 

d

 

 

 

 

M

 

.

 

 

 

 

 

 

 

 

 

dt

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

 

 

 

 

 

 

 

 

 

 

W p Ap2 Bp c 1

 

(3.8)

представляет собой матрицу операторов динамической податливости свободной системы, элементы которой обладают свойствами передаточных функций линейных систем и связывают преобразования Лапласа r-го входа и s-го выхода так, что решение (3.4) с учетом (3.5) в покомпонентном виде можно записать так:

n

s p Wrs p Mr p , s 1, 2, , n. (3.9)

r 1

90

Так как матрицы A, B, C симметричны, то и матрица W p также симметричная. Компоненты матрицы Wrs p вычисляются по формуле

W p

rs p

,

(3.10)

p

rs

 

 

 

 

 

где rs p – алгебраическое дополнение элемента r-й строки и s-го столбца этого определителя;

p det Ap2 Bp c – характеристический определитель

системы (3.10).

Обычно в машинах диссипация пренебрежимо мала, что позволяет в (3.8) положить B 0. Собственные формы и частоты систе-

мы удовлетворяют системе (3.8) при B 0, M 0.

Кинетическая и потенциальная энергии линеаризованной системы имеют вид

 

 

1

 

 

1

c , ,

 

 

T 2

A , ,

П

2

(3.11)

 

 

n

 

 

 

n

 

aik i k ;

c , cik i k .

 

A ,

 

 

 

 

 

 

 

 

 

 

 

i,k 1

 

 

 

i,k 1

 

Так как квадратичные формы определенно-положительные, то существует вещественная неособенная замена переменных

U

 

, detU 0,

 

1, , n ,

(3.12)

 

 

приводящая к суммеквадратов сразу обе квадратичные формы (3.11):

T

1

n 2

,

П

1

n

2

 

j

 

j j .

 

2 j 1

 

 

2 j 1

 

91

Обобщенные координаты j называются главными или нор-

мальными. Уравнения движения (3.6) в главных координатах запишутся в виде несвязанных n-уравнений второго порядка:

 

 

 

 

 

 

(3.13)

 

 

 

 

 

j j j 0, j 1, 2, , n.

Так как все j положительны, то решения (3.13) имеют вид

 

 

 

 

 

 

j c j sin jt j , j 1, 2, , n,

 

 

где j

j – частоты колебаний;

 

 

c j , j

– произвольные постоянные.

 

 

Тогда

 

 

представляется в виде

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

c ju j sin jt j .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

Здесь

u

u1, , un – вектор амплитуд.

 

 

Подставляя (3.10) в (3.7) и затем в (3.6), получим

 

 

 

 

 

 

u

 

 

c A

u

0.

 

 

Условие

0 приводит к уравнению

 

 

 

 

 

 

 

 

 

det c A 0.

(3.14)

Уравнение (3.14) определяет собственные частоты

i

i ,

i 1, , n.

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (3.6) в главных координатах имеет вид

 

 

 

 

 

 

 

 

 

Aг p2 Cг

 

Mг,

 

 

 

 

 

 

 

 

 

 

 

 

где Aг, Сг

 

– диагональные матрицы размера n 1 n 1 .

 

92

Запишем Aг в виде

Aг diag 0 , 1, , n ;

s as T Aus , s 0, , n;

n

0 Ir Ic ,

r 0

r diag 0 , 1, , n ,

s us T cu,

s 0, , n.

Представим компоненты вектора Мг в виде

MГs us T MГ us0M g usnMn .

Тогда (3.6) записывается в виде n 1 независимых уравнений:

Ic p2 0 M g Mn ;

ar p2 r r M g rnMn ,

(3.15)

r 1, 2, , n.

Решение системы (3.15) запишем в виде

 

 

0

 

1

 

M g Mn ,

 

 

 

 

Ic p

2

 

 

 

 

 

 

 

 

 

 

r

1

 

 

M g rnMn ,

r 1, , n.

a p2

 

r

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

Подставляя (3.12) в (3.6), получим выражения для -решения неоднородных уравнений (3.6), (3.7), описывающих динамику сис-

93

темы под действием внутренних сил. Покомпонентное выражение имеет вид

 

 

 

 

 

 

 

 

1

 

 

n

 

urs

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

M g

 

 

I

 

p2

 

 

 

 

 

 

 

 

 

 

 

 

c

r 1 a p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

r

n

 

 

1

 

 

 

n

 

 

u

u

rm

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

Mm

, s 0,1, , n.

I

 

p2

 

 

 

 

 

m 1

c

 

r 1 a p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

r

 

 

 

Выражения для компонентов передаточной матрицы запишутся в безразмерном виде:

 

 

 

 

1

 

 

n

 

m

 

 

 

 

Wrs p

 

 

 

 

rs

,

 

r, s 0,1, ,n;

Ic p2

 

 

 

 

 

 

 

 

m 1 m2 p2 1

 

 

 

 

 

 

m

m 1;

rs umrums m1.

Для установившегося процесса

p j ,

тогда динамические по-

датливости Wrs j

имеют вид

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

 

m

 

 

 

 

Wrs j

 

 

 

rs

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

r, s 0,1, , n.

 

 

 

2

 

2

2

 

 

 

 

 

Ic

m 1

1 m

 

 

Векторы

 

 

 

 

 

 

А, С, причем при любых

hs ортогональны в матрицах

s m выполняются условия

 

 

 

 

 

 

 

 

 

 

u

sT Aum 0,

u

sT cum 0.

Собственные частоты k находятся из уравнения

 

 

 

 

 

 

 

det С Ak2 0.

(3.16)

94

Выражение (3.16) является условием существования ненулевых амплитуд us , удовлетворяющих системе

С Ak2

u

s 0,

s 0,1, , n.

Собственные частоты выражаются через собственные формы (формула Рэлея):

2 uT Сus ks usT Сu ,

h s

причем считается, что uso 1, h0 T 1,1, ,1 . Применим к (3.9) обратное преобразование Лапласа

s t

s i

s p Mr p e ptdp,

 

s i

получим выражение для s t в виде

t

s t Wrs t Mr d .

r 10

В случае резонанса, что соответствует в (3.16)

1 k

m m

совпадению с одной из собственных частот системы в знаменателе Wrs j , он становится бесконечным. Это происходит, когда период

(частота) действия на систему, совпадает с собственной частотой.

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

asin t,

95

общее решение которого имеет вид

t csin( t ) * t ;

(3.17)

* t

a

 

 

t cos t.

(3.18)

2

Как следует из (3.18) функция * t растет, а потому линеаризо-

ванные уравнения вибрации должны заменяться более точными нелинейными.

Учет вязкости (диссипации) в системе также позволяет получать более реалистичные конечные режимы движения.

3.1.6. Антирезонансные частоты

Из анализа выражения (3.17) видно, что если закрепить массы в r-й и s-й материальных точках системы, величина Wrs j обра-

щается в 0 при совпадении с одной из собственных частот системы. Эти частоты называются антирезонансными. В отличие от резонансных частот, общих для всех Wrs , антирезонасные частоты

для каждой динамической податливости свои. Антирезонасные частоты диагональных элементов матрицы динамических податливо-

стей Wrr j располагаются между их резонансными частотами:

0 k1 k1 k2 kn kn ,

где k1 – антирезонасные, а ki – резонансные частоты.

При r s расположение антирезонансных частот зависит от перемен знака в ряду чисел

uoruos , u1ru1s , u2ru2s ,

(3.19)

Если знаки umrums и um 1,rum 1,s совпадают, то между

km и km 1

имеется антирезонансная частота km 1. Если знаки различные, то соответствующая антирезонансная частота отсутствует.

96

Число перемен знака в ряду (3.19) всегда равно s r.

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

3.2.Модели механических подсистем

иметод механических цепей

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

гими элементами, третьи — устройствами с трением. Отдельно представляют элементы, поставляющие энергию в механическую систему и возбуждающие ее движение: активные элементы или источники. Элементы, не имеющие независимых источников сил или кинематических величин, называют пассивными.

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

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

97

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

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

Переменные двухполюсника. Двухполюсник можно характеризовать силами F1 и F2 , приложенными к его полюсам 1 и 2, и пара-

метрами движения полюсов (рис. 3.6).

Рис. 3.6. Двухполюсник Д с полюсами 1 и 2

 

Для перемещений полюсов d1 и d2

 

имеем

 

 

 

 

 

 

d d x0

(x x )x0;

 

 

 

 

 

 

1

 

 

1

 

1

 

10

 

 

 

 

 

 

d

2

d

2

x0

(x x )x0

,

 

 

 

 

 

 

 

 

 

2

 

20

 

 

 

 

где xo – единичный вектор оси Ox;

 

 

 

 

 

 

x10 , x20 , x1, x2

– исходные и текущие координаты точек 1 и 2.

Скорости i и ускорения ai

полюсов определяют из выражений

i i x

0

 

 

 

0

;

ai ai x

0

i i x

0

,

i 1, 2.

 

di di x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Относительное движение полюсов удобно характеризовать с по-

мощью кинематических векторов двухполюсника d, и a (d

98

относительное перемещение, – относительная скорость; а – относительное ускорение полюсов; ниже эти векторы обобщенно обозначаются буквой k ):

d dx0

 

 

 

 

 

d

i

d

j

 

 

 

e ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

x0

 

 

 

 

 

 

e

 

;

(3.20)

 

 

j

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

a ax0

 

 

a

 

a

j

 

 

e

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

a

 

 

 

где d , и a кинематические переменные двухполюсника,

обобщенно обозначаемые буквой k;

ek k d, , a – единичные векторы, направленные от полюса

сбольшим алгебраическим значением кинематической переменной

кполюсу с ее меньшим алгебраическим значением.

Положительно направленным векторам d и отвечает сближе-

ние полюсов двухполюсника, поэтому при условии, что координата xi полюса i меньше координаты j полюса j(xi x j ):

k ki k j ;

k d, , a.

Так, для указанных на рис. 3.7 индексации полюсов и направления оси Ох

d d1 d2;

F1 F1x0 Fx0 F;

(3.21)

F2 F2 x0 Fx0 F,

где F воспринимаемая двухполюсником сила; F силовая переменная двухполюсника.

99