Обобщенной силой нестационарной голономной одностепенной

механической системы называется отношение виртуальной мощности системы сил (мощности системы сил на виртуальных скоростях) к обобщенной скорости:

Q = P

*

&

*

 

 

 

 

*

(1.50)

 

 

 

 

 

/ q, при P = Fi

(vi , ri ,t) vi

+ M j ω j

 

В общем случае Q есть функция трех параметров: Q =Q(q&, q,t). В (1.50) силы зависят от действительных (не виртуальных) скоростей. Уравнение Лагранжа для нестационарной одностепенной системы имеет такой же вид, как и для стационарной системы:

d

T

 

 

T

Q = 0

(1.51)

 

 

&

 

q

,

 

dt

q

 

 

 

т.е. полная производная по времени от обобщенного импульса Tq&, уменьшенная на величину Tq равна обобщенной силе системы

приложенных сил, определяемой по формуле (1.50) через мощность реальных сил на виртуальных скоростях.

Выполнив процедуры дифференцирования трехчленного выражения кинетической энергии (1.48), указанные в уравнении (1.51), получим окончательно динамическое уравнение в форме ОДУ второго порядка вида:

&&

 

1

&2 a2

&

a2

 

a1

 

a0

 

 

(1.52)

+b

 

+

Q

= 0

,

q

2

q

q

+ q

t

q

 

 

 

 

t

 

 

 

 

 

 

 

 

при ak

= ak (q,t),

k = 0,1,2, b = a21

 

 

 

 

 

К закону динамики нестационарной механической системы вида

(1.52) присоединяют начальные условия, начальную фазу движения

(1.53)

 

 

q(0) = q0 , q(0) = q0 или [q(0), q(0)] =[q0 , q0 ]

 

 

 

 

 

&

 

&

 

 

 

 

 

&

 

&

 

 

 

 

 

Математическую модель движения рассматриваемой механической

системы (1.52)-(1.53), можно записать в нормальной форме Коши:

(1.54)

&

 

a0

 

a1

 

a2

 

1

 

2 a2

 

&

*

 

1

u = b Q +

q

t

u

t

 

u

q

 

, q = u, при Q = P

 

/ u,

b = a2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(0) =u0 , q(0) = q0

 

 

 

(1.55)

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

21

Пример. Электродвигатель с вращающимся кривошипом.

Рис. 1 Вибрация подпружиненного электродвигателя с вращающимся кривошипом

Электродвигатель массой m1 , включающий массу корпуса, платформы

и четырех направляющих колес, установлен на пружинах общей жесткости с2 , с массой пружин m2 . Радиусы инерции каждого из четырех

направляющих колес i3 , радиусы колес r3 , масса колес m3 . На оси О двигателя установлен неуравновешенный маховик массой m4 с расстоянием до центра масс r4 , радиус инерции маховика относительно

оси О равен i4 , маховик равномерно вращается согласно кинематическому уравнению φ= ω4 t , при ω4=2. Задан общий момент вязкого трения на осях M3= b3 ω3 . Составить динамическое уравнение устройства.

Решение Уравнением ϕ= ω4t =2t задана нестационарная связь.

Пружины не создают связей, они лишь создают силы и имеют инерционность, которую согласно методу Рэлея учтем, прибавив к массе m1 двигателя третью часть массы пружин m2 . Имеем нестационарную

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

&

На рис.1 показано

Фазовый вектор системы: V=[v,y] при v = y .

положительное кинематическое состояние системы,

когда v>0 ,

y>0 .

Получим динамическое уравнение. Выразим силы

и скорости

через

фазовые координаты v , y . Имеем:

 

 

22

ω

 

=

v

,

v

= v +v

при v

=ω r

 

 

 

 

 

 

3

, v = vj

, v

OC,

 

 

 

r3

 

 

C

CO

 

 

CO

4 4

 

 

 

 

CO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.56)

ω

3 =ω3k , F2 = −c2 y j ,

M3 = −b3ω3 k = −(b3v / r3 )k ,

 

G1 = −m1g j , G4 = −m4 g j

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

 

+

1

m2

 

2

2

2

2

+ m4 (v

2

2

2

)

2T = m1

3

+ 4m3 v

 

+ 4m3i3 v

 

/ r3

 

+ 2v vCO +i4ω4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где v vCO = −vω4r4 sin(ω4t) .

Отсюда получаем кинетическую энергию в форме квадратного полинома относительно v :

 

 

 

T =

1 a v2 + a v + a

,

 

 

 

(1.57)

при обозначениях

 

 

2

 

2

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.58)

a2 = m1 + m2 / 3 + 4m3 + 4J3 / r32 + m4 ,

 

 

h4 = m4r4ω4 ,

a = −h sin(ω

t),

a

 

= m i2ω2

/ 2,

ω

4

= const

 

1

4

4

 

0

 

4 4

4

 

 

 

 

 

Частные производные от коэффициентов, входящие в (1.52)

(1.59)

a2

= 0,

a1

= −h ω

4

cos(ω

t),

a2

= 0, a0 = 0

y

 

t

4

4

 

 

t

 

 

y

 

 

 

 

 

 

 

 

 

 

 

Обобщенная сила Q по обобщенной скорости v = y& находится через вычисление виртуальной мощности системы всех приложенных сил на виртуальных скоростях точек приложения сил. При этом в выражениях сил зависящих от скоростей не следует заменять действительные скорости на виртуальные. Не принимаем во внимание крутящий момент электродвигателя вместе с моментом трения в подшипнике O , осуществляющих нестационарную связь – программное движение кривошипа, поскольку их мощность на виртуальном движении при фиксированном в связях параметре t (т.е. при ϕ = const ) равна нулю.

Мощность силы G4 подсчитываем на виртуальной скорости vC* = v . Центр

масс пружин амортизатора имеет скорость v / 2 . Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P* =G

V

+G V / 2 +G V + F V +4M

3

ω

3

= −(G +G / 2 +G + F + 4M

3

/ r )V

1

2

 

 

4

2

 

 

 

 

 

 

1

2

4

2

4

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

b3

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.60)

 

 

 

*

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q = P

/ v = −

y +c2 y +G1 +G2 / 2 +G4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в уравнение Лагранжа (1.52) выражения (1.59),(1.60) получаем динамическое уравнение устройства:

23

&&

 

b3

&

 

 

 

 

1

 

+G2

/ 2 +G4

h4ω4 cos(2t)

= 0,

b = a2

y

+b

r3

y +c2 y +G1

 

 

 

 

 

 

 

 

Окончательно имеем линейное неоднородное ОДУ второго порядка с постоянно действующими периодическими возмущениями следующей

структуры:

 

 

2

y +h0

h cos( pt) = 0, 0 n < k

(1.61)

&&

&

 

y +

2ny +k

 

 

при обозначениях

 

 

 

 

 

n =b /(2r3a2 ),

k = c2 / a2 , h0 = (G1 +G2 / 2 +G4 ) / a2 ,

h = h4ω4 / a2

Общее решение уравнения (1.61) состоит из частного периодического решения, имеющего гармоническую и постоянную составляющие, и общего решения однородного уравнения с двумя константами С и β :

y = H cos( pt

α) h0 / k

2

+Ce

nt

cos(k1t β), k1 = k

2

n

2

(1.62)

 

 

 

 

 

H =

 

h

 

,

cos(α) = h1H

(k2 p2 )

 

 

 

 

 

 

 

 

 

 

 

 

(1.63)

 

(k2 p2 )2 + 4n2 p2

 

 

 

 

 

 

 

 

 

 

Константы

С

и β

 

определяются

через начальные

условия

y(0) = y0 , y(0) = v0 y

посредством подстановки их в общее решение (1.62) и

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в производную от него.

 

 

 

 

 

 

 

 

 

 

 

Согласно

уравнению (1.62)

 

система

совершает

 

 

двухчастотные

колебания около

положения

y = −h / k2 с

собственной

 

циклической

частотой k1 и вынужденной частотой p . При этом амплитуда собственных

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

y = H cos( pt α) + yСТ при yСТ = −h0 / k2

(1.64)

Иными словами, система совершает гармонические колебания согласно формуле (1.64), асимптотически устойчивые относительно начальных возмущений в том смысле, что любой переходный процесс, начинающийся

с состояния

y(0) = y0

, y(0) = v0 , асимптотически переходит

в

 

 

&

 

периодические колебания относительно статического положения y = yСТ

с

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

Замечание В случае нестационарной одностепенной механической

системы уравнение Лагранжа и

уравнение энергии T& = P

приводят

к

различным дифференциальным

уравнениям. В данном

примере

из

уравнения T& = P получается дифференциальное уравнение, содержащее мощность крутящего момента. Такое уравнение можно использовать для

24

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

Замечание Имеется другой способ решения одностепенной нестационарной задачи. Можно обобщенную координату q переобозначить q1 и в уравнении нестационарных связей время t переобозначить q2 , вводя

условно дополнительную координату, дополнительную степень свободы. Получаем стационарную двухстепенную механическую систему для которой составляем два уравнения Лагранжа. Затем в этих уравнениях полагаем q2 =t , первое из уравнений есть искомое ОДУ, второе можно

отбросить или использовать для реакций нестационарных связей. В рассмотренном примере в качестве q2 можно принять угол поворота кривошипа и в конце принять q2 =ω4t

25

ЛАБОРАТОРНАЯ РАБОТА №1

ОПРЕДЕЛЕНИЕ КИНЕТИЧЕСКОЙ ЭНЕРГИИ МНОГОЗВЕННОГО МЕХАНИЗМА С ДВУМЯ СТЕПЕНЯМИ СВОБОДЫ

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

Цель работы

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

1. Описание установки

Рис.1.1. Окно виртуальной лабораторной установки

26

Изменяемые начальные параметры механической системы:

m11,m12 ,m21,m22 массы соответствующих звеньев, остальные звенья установки считаются невесомыми

k1,k2 коэффициент жесткости верхней и нижней пружины

k3 коэффициент вязкого трения в демпфере

k4 , k5 коэффициенты вязкого трения в подшипниках каждого колеса и в подшипниках маятника

i2 радиус инерции диска маятника, с увеличением радиуса инерции увеличиваются геометрические размеры диска

r3 радиус эксцентрика электродвигателя

r расстояние от оси подвеса маятника до центра диска массой

x(0), x(0),ϕ(0),ϕ(0) заданные начальные значения, фазовых

 

&

&

координат механической системы

Неизменяемые параметры механической системы:

r1 ,i1 радиус и радиус инерции каждого колеса платформы, их значения фиксированы и составляют r1 = 0.6 , i1 = 0.424

l длина стержня маятника, звено является однородным тонким стержнем, значение параметра фиксировано и равно l = 6.48

Параметры моделирования:

n шаг дискретизации, этот параметр влияет на скорость и

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

27

Рис. 1.2. Параметры системы с двумя степенями свободы

2. Математическая модель

Общая трехчленная формула для двойной кинетической энергии

плоского движения тела имеет вид

2T = mv2

+ J

 

ω2 + 2mv

 

v

 

,

v =ω OC,

v OC,

 

 

 

 

 

oz

o

co

v =ω

×OC ,

(1.1)

o

 

 

 

 

co

co

co

 

где vo - скорость полюса тела, OC - расстояние от полюса O до центра массы С,Oz – ось вращения, Joz - момент инерции тела относительно оси, проведенной через полюс перпендикулярно плоскости движения, vco -

относительная скорость вращения центра масс тела С вокруг полюса O .

vo vco =vo vcocos(vo vco )

В случае, когда за полюс принят центр масс тела, получаем из (1.1) двухчленную формулу:

2T = mv2

+ J

сz

ω2

(1.2)

с

 

 

 

28

Кинетическая энергия механической системы равна сумме кинетических энергий всех ее звеньев:

 

T =3T11 +T12 +T2

при T2 =T21 +T22

(1.3)

Здесь

 

 

T11

- кинетическая энергия каждого из трех колес массы m12

 

T12

- кинетическая энергия платформы массой m12

 

T21

- кинетическая энергия

стержня маятника массой

m21

совершающего плоское движение

T22 - кинетическая энергия диска массой m22 , совершающего плоское

движение

T2 - кинетическая энергия маятника, состоящего из стержня и диска.

T2 =T21 +T22

(1.4)

При обозначении (1.4), окончательная формула (1.3) принимает вид кинетической энергии механической системы:

T =3T11 +T12 +T2

(1.5)

2.1. Кинетическая энергия платформы и колес

Двойная кинетическая энергия платформы, движущейся поступательно определяется формулой

2T

= m v2

(1.6)

12

12

 

Кинетическая энергия трех колес с общей массой 3m11 , совершающих

плоское движение состоит из кинетической энергии поступательного движения вместе с полюсом О и кинетической энергии вращения вокруг О в поступательной системе отсчета

T11 =

3

(m11v2

+ Jcω2 ),

при Jc = m11i12 ,

ω = v / r1

или,

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

i2

 

2

 

 

 

 

T11 =

 

m11

1

+

1

v

 

(1.7)

 

 

 

2

2

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

1

 

 

Коэффициент инерции системы трех колес, приведенный к скорости

a11 =3m11(1+i12 / r12 )

2.2.Кинетическая энергия маятника

В данном механизме маятник совершает плоское движение. Найдем его кинетическую энергию T2 двумя способами.

29

Первый способ

Из выражения (1.2) следует, что кинетическая энергия маятника (Рис.6) состоит из энергии стержня и диска, совершающих плоское движение, из четырех слагаемых,

 

2T = m v2

+ J стω2 + m v2

+ J дискω2 ,

(1.8)

 

2

21 С1

С1

22 С2

С2

 

где vC1 -

скорость центра масс C1

стержня,

vС2 - скорость центра масс

C2 диска,

JСст1 и JСдиск2

- моменты инерции стержня и диска относительно

их центральных осей, перпендикулярных плоскости движения.

Рис. 1.3. Маятник

Скорости центров масс находятся как векторные суммы скорости

полюса и относительной скорости движения вокруг полюса:

vC1 = v + vC1_ O ,

vC 2 = v + vC 2 _ O

v = x

 

,

 

i

 

&

 

 

 

(1.9)

vC1_ O = vC1_ O e,

vC 2 _ O = vC 2 _ O e

e OC1

(1.10)

Здесь i -орт линейной координаты x , e - сопровождающий касательный орт угловой координаты ϕ направленный против часовой стрелки, O -

точка подвеса маятника, v - скорость точки O , равная скорости платформы, vC1_ O и vC 2 _ O - относительные скорости вращения центров

масс стержня и диска маятника вокруг точки подвеса O , vC1_ O и vC 2 _ O -

30