Q(x) x=0 = 0 cf + 23 mg = 0 f = 2c3 mg

Подставляя выражение статической деформации f в обобщенную силу получаем окончательно: Q = −cx . Динамическое уравнение колебаний имеет вид:

&&x + k

2

 

 

&

 

 

 

 

x = 0 , с начальным состоянием x(0) = l f , x(0) = v0 .

 

Решение полученной начальной задачи Коши имеет вид:

 

 

 

x = x(0)cos(kt) +

v0

sin(kt) = (l f )cos(kt) +

v0

 

sin(kt).

 

 

k

k

 

 

 

 

 

Пример 2. Решить предыдущую задачу при дополнительном условии, что наклонная поверхность вибрирует в вертикальном направлении с ускорением aO = Acos pt j1 , где орт j1 направлен вертикально.

Решение. В этом случае следует приложить в центре масс каждого звена переносную силу инерции Sek = −mk aO (наряду с силой Gk = mk g ), после этого систему отсчета Oxy можно условно считать неподвижной (инерциальной). Приложив в центре масс системы корпус-колеса силу GO = m(g aO )

(вместо силы G

= mg ), находим Q = P / x = −cx + 3 / 2 m(g aO ). В результа-

 

 

 

&

 

те получаем динамическое уравнение колебаний вида

 

&&

2

x h a1 cos pt = 0 при a1 = mA

3 2mпр ,

x + k

 

имеющее дополнительное гармоническое слагаемое.

5.10. Уравнение Лагранжа для голономной нестационарной одностепенной системы

Геометрическая связь в механической системе называется нестационарной, если она записывается в виде функционального уравнения, содержащего параметр t , т.е. если она изменяется с течением времени. Примеры нестационарных связей: математический маятник на нити переменной длины l=5t , прибор на платформе, совершающей заданное во времени программное движение, механизм с заданными кинематическими уравнениями движения одних звеньев относительно других. Голономная система называется одностепенной нестационарной, , если она подчинена хотя бы одной нестационарной связи наряду со стационарными связями, но при этом сохраняет одну степень свободы. Положение одностепенной механической системы определяется одной выбираемой обобщенной координатой q в каждое конкретное мгновение t .

В стационарной голономной системе радиус-вектор ri любой i-ой точки

системы можно найти на основании геометрического изображения объекта (желательно изображать положение q > 0) , как функцию обобщенной коор-

динаты ri = ri (q) . Скорость точки определяется функцией vi = ∂ri q q& .

88

Для нестационарных систем приходим к более сложным выражениям

ri = ri (q,t) , в связи с чем скорость v

 

содержит дополнительное слагаемое, не

 

 

 

i

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

содержащее q :

 

 

 

 

 

 

 

 

(5.47)

&

 

+∂ri

t

, где

 

= ∂ri

&

vi = ∂ri q q +∂ri

t = vi

vi

q q .

 

Составляющую vi скорости vi , вычисляемую при фиксированных значениях

параметра t в нестационарной связи, получаемую за счет изменения только обобщенной координаты q , назовем виртуальной скоростью.

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

 

1

&2

&

 

ak = ak (q,t),k = 0,1,2 .

(5.48)

T =

2 a2q

+ a1q + a0

при

 

Мощность системы приложенных к механической системе сил будем вычислять на виртуальных скоростях точек приложения сил, т.е. при условии, что в нестационарных связях параметр t закреплен в своем положительном значении. При этом не будем принимать во внимание не имеющие мощности идеальные реакции нестационарных и стационарных связей при закрепленном параметре t в связях. Имеем выражение виртуальной мощности системы сил:

P* = Fi

vi* +M

j ω

*j

(5.49)

Здесь силы Fi и моменты пар M j в общем случае могут зависеть от реальных скоростей vi и угловых скоростей ωj , а также - от координат и от времени.

Выразив в выражении (5.49) виртуальные скорости vi* , ω*j , а также действительные скорости и координаты через фазовые переменные q& и q , полу-

чим мощность P* (q&, q,t) ,в которой обязательно имеется общий множитель q&.

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

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

Q = P

*

&

*

 

 

 

 

*

(5.50)

 

 

 

 

 

/ q, при P = Fi

(vi , ri ,t) vi

+ M j ω j

 

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

Уравнение Лагранжа для нестационарной (и стационарной) одностепенной системы имеет вид:

d

T

 

 

T

(5.51)

 

 

&

 

q

= Q ,

 

dt

q

 

 

 

 

 

 

 

 

89

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

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

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

 

 

1

 

2

a2

 

a2

q +b

 

q

 

 

+ q

 

&&

 

2

&

 

q

&

t

 

 

 

 

 

при ak

= ak (q,t),

 

a1

 

a0

Q

(5.52)

+

= 0,

 

t

 

q

 

 

k = 0,1,2, b = a21

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

информацию о предистории движения

 

 

 

 

(5.53)

 

 

q(0) = q0 ,

q(0) = q0

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

 

 

 

 

&

 

 

&

 

 

 

 

 

 

&

&

 

 

 

 

 

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

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

(5.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

 

 

 

 

(5.55)

Вместо начальных условий (5.53), (5.55) могут быть заданы граничные

условия, или условия периодичности искомого решения. Пример

Рис. 18 Вибрация электродвигателя

90

Электродвигатель массой m1 , установлен на пружинах общей жесткости с2 , с массой пружины m2 . Радиусы инерции каждого из четырех направляющих колес i3 , радиусы колес r3 . На оси О двигателя установлен неуравновешанный маховик массой m4 с расстоянием до центра масс r4 , радиус инер-

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

Решение

Уравнением ϕ= ω4t =2t задана нестационарная связь. Пружины не соз-

дают связей, они лишь создают силы и имеют инерционность, которую со-

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

степенью свободы, в которой за обобщенную координату примем направлен-

ное вверх смещение от положения оси О, соответствующего недеформиро-

ванному состоянию пружины.

 

 

 

Фазовый вектор ситемы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

На рис.18 показано поло-

 

 

 

 

V=[v,y] при v = y .

 

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

когда v>0 ,

y>0 . Получим

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

через фазовые коорди-

наты v , y . Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

=

v

,

 

v = v

+v

 

 

при v

 

=ω r

 

 

 

 

 

 

 

 

 

3

 

 

 

, v = vj

, v

 

OC,

 

 

 

 

 

 

 

 

 

 

 

 

r3

 

 

 

C

 

CO

 

 

 

 

 

 

CO

 

 

4

4

 

 

 

 

 

 

 

 

CO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.56)

ω

3 =ω3k , F2 = −c2 y j , M3 = −b3ω3 k = −(b3v / r3 )k ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 = −m1g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

j

, G

4 = −m4 g

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Двойная кинетическая энергия устройства:

 

 

 

 

 

 

 

 

 

 

 

2T

 

 

1

m2

 

 

2

 

 

 

 

2

 

2

 

2

+ m4 (v

2

 

 

 

 

 

 

 

 

 

2

2

)

 

 

 

 

 

= m1 +

3

v

 

 

+ m3i3 v

 

/ r3

 

 

+ 2v vCO +i4ω4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

ма относительно v :

 

 

 

 

 

 

 

 

 

 

1 a v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.57)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T =

+ a v + a

,

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

1

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.58)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

h4 = m4r4ω4 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a = −h sin(ω

t),

 

a

= m i2ω2

/ 2,

 

ω

4

= const

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

4

 

 

 

 

4

 

 

 

 

 

0

4 4

4

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

a2 = 0,

 

 

 

a1

= −h ω

4

cos(ω

t),

a2

= 0,

a0 = 0

(5.59)

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

t

 

 

 

4

 

 

 

4

 

 

 

 

 

t

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обобщенная сила Q

по обобщенной скорости

 

&

 

 

 

 

v = y находится через

вычисление виртуальной мощности системы всех приложенных сил на вир91

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

сте с моментом трения в подшипнике O , осуществляющих нестационарную

связь – программное движение кривошипа, поскольку их мощность на виртуальном движении при фиксированном t (т.е. при ϕ = const ) равна нулю.

 

 

 

 

 

 

 

 

 

 

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

Мощность силы G

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

масс пружин амортизатора имеет скорость 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

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.60)

 

 

 

 

 

 

 

 

 

Q = P

*

 

 

 

&

 

 

 

+G1

+G2 / 2 +G4

 

 

 

 

 

 

 

 

 

 

 

 

 

/ v = −

y +c2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в уравненние Лагранжа (5.52) выражения (5.59),(5.60) по-

лучаем динамическое уравнение устройства:

 

 

 

 

 

 

 

 

 

&&

 

 

 

b3

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

+G2 / 2 +G4 h4ω4 cos(2t)

= 0,

b = a2

 

 

 

 

 

 

 

 

y

+b

y + c2 y +G1

 

 

 

 

 

 

 

 

 

 

 

r3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

постоянно действующими периодическими возмущениями следующей структуры:

&&

&

2

y +h0

h cos( pt) = 0,

0 n < k

(5.61)

y

+2ny + k

 

 

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

n =b /(2r3a2 ), k = c2 / a2 , h0 = (G1 +G2 / 2 +G4 ) / a2 , h = h4ω4 / a2

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

y = H cos( pt α) h0 / k

2

+Ce

nt

cos(k1t

β), k1 =

k

2

n

2

(5.62)

 

 

 

 

 

H =

 

h

,

α = arccos(h1H (k2 p2 ))

 

 

 

 

(k2

p2 )2 +4n2 p2

 

 

 

(5.63)

 

 

 

 

 

 

 

 

 

 

Константы

С и β

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

через

начальные

условия

y(0) = y0 , y&(0) = v0 y посредством подстановки их в общее решение (5.62) и в

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

Согласно уравнению (5.62) система совершает двухчастотные колеба-

ния около положения y = −h / k2 с собственной циклической частотой k1 и

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

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

92

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

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

(5.64)

Иными словами, система совершает гармонические колебания согласно формуле (5.64), асимптотически устойчивые относительно начальных возмущений в том смысле, что любой переходный процесс, начинающийся с состояния y(0) = y0 , y&(0) = v0 , асимптотически переходит в периодические колеба-

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

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

Примечание В случае нестационарной одностепенной механической системы уравнение Лагранжа и уравнение энергии T& = P приводят к различ-

ным дифференциальным уравнениям. В данном примере из уравнения T& = P

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

93