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

3168

.pdf
Скачиваний:
0
Добавлен:
21.11.2023
Размер:
336.28 Кб
Скачать

N вар.

1

2

3

4

5

6

7

8

9

10

Схема №

1

1

1

1

1

2

2

2

2

2

М1

1

1

n+2

n+2

1

n+1

n+1

2

2

n+1

а1

4

2

4

2

4

2

2

8

8

8

f1

1

0

1

0

1

1

2

1

2

0

М2

n+2

n+2

1

1

n+2

2

2

n+1

n+1

2

а2

2

4

2

4

2

10

10

10

2

2

f2

0

1

0

1

0

0

1

0

1

1

К

1

1

1

1

1

1

1

1

1

1

а

3

2

1

2

3

1

1

1

1

1

P

50

50

50

50

50

10

10

10

10

10

Х1(0)

1

0

1

0

1

0

1

0

1

0

V1(0)

0

1

0

1

0

1

0

1

0

1

Х2(0)

1

1

0

0

1

1

0

0

1

1

V2(0)

0

0

1

1

0

0

1

1

0

0

N вар.

11

12

13

14

15

16

17

18

19

20

Схема №

3

3

3

3

3

4

4

4

4

4

М1

n+2

n+2

1

1

n+2

1

n+2

n+2

1

1

К1

1

1

1

2

2

2

2

1

3

3

а1

2

2

2

4

4

4

4

2

2

2

f1

1

0

1

0

1

0

1

0

1

0

М2

4

4

n+2

n+2

1

n+2

4

4

n+2

n+2

К2

2

2

1

1

1

3

3

2

1

1

а2

4

4

2

2

2

2

2

2

4

4

f2

0

1

0

1

0

1

0

1

0

1

K

 

 

 

 

 

1

1

1

1

1

а

 

 

 

 

 

1

1

1

1

1

 

 

 

 

 

20

 

 

 

 

 

P

 

 

 

 

 

2

2

1

3

1

Х1(0)

1

0

1

0

1

0

1

0

1

0

V1(0)

0

1

0

1

0

1

0

1

0

1

Х2(0)

1

1

0

0

1

1

0

0

1

1

V2(0)

0

0

1

1

0

0

1

1

0

0

 

 

 

 

 

 

 

 

 

 

 

N вар.

21

22

23

24

25

26

27

28

29

30

Схема №

5

5

5

5

5

6

6

6

6

6

М1

4

4

n+2

n+2

1

1

4

4

1

1

а1

2

2

2

4

4

4

4

2

2

2

f1

1

0

1

0

1

0

1

0

1

0

М2

n+2

n+2

1

1

n+2

n+2

n+2

n+2

n+2

n+2

а2

4

4

2

2

2

2

2

2

4

4

f2

0

1

0

1

0

1

0

1

0

1

а

4

4

4

4

4

1

1

1

1

1

P

20

20

20

20

20

16

16

16

16

16

Х1(0)

1

0

1

0

1

0

1

0

1

0

V1(0)

0

1

0

1

0

1

0

1

0

1

Х2(0)

1

1

0

0

1

1

0

0

1

1

V2(0)

0

0

1

1

0

0

1

1

0

0

 

 

 

 

 

 

 

 

 

 

 

N вар.

31

32

33

34

35

36

37

38

29

40

Схема №

7

7

7

7

7

8

8

8

8

8

М1

1

1

n+2

n+2

1

1

1

n+2

n+2

1

К1

 

 

 

 

 

0

0

1

1

1

а1

4

4

4

6

6

 

 

 

 

 

f1

1

0

1

0

1

1

0

1

0

1

М2

n+2

n+2

1

1

n+2

n+2

n+2

1

1

n+2

К2

 

 

 

 

 

1

1

0

0

0

а2

2

2

1

1

1

 

 

 

 

 

f2

0

1

0

1

0

0

1

0

1

0

К

1

1

1

1

1

1

1

1

1

1

а

2

1

2

1

2

 

 

 

 

 

P

10

10

10

10

10

 

 

 

 

 

Х1(0)

1

0

1

0

1

1

0

1

0

1

V1(0)

0

1

0

1

0

0

1

0

1

0

Х2(0)

1

1

0

0

1

1

1

0

0

1

V2(0)

0

0

1

1

0

0

0

1

1

0

21

 

Указания к разработке динамических уравнений упругих систем.

1.)

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

Т1

=

1

m × v2

,

Т2 =

1

J ×ω 2

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

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

для вращательного движения с угловой скоростью ω .

2.)

Потенциальная энергия упругого элемента системы при малых перемещениях х

вычисляется по формуле,

П

 

=

1

к × х2

, где к эффективная жесткость упругого

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

элемента. В случае если упругий элемент нагружен силой N0 , то его потенциаль-

ная энергия

 

П2 = N 0 × х .

Потенциальная энергия в поле силы тяжести на h

уровне П3

= m × g × h

 

 

 

 

 

 

3.)

Эффективная жесткость упругого элемента равна силе или моменту , приводя-

щих к единичному смещению элемента (эквивалент пружины Fу = -к × х ). Для

упругого защемленного стержня, нагруженного на свободном конце силой или моментом, эффективные жесткости будут:

при сжатии-растяжении

к =

ES

, при кручении

к =

JG

, при изгибе попереч-

l

 

l

 

 

 

 

 

 

 

 

 

 

ной силой к =

3ЕJ

при изгибе моментом к =

2ЕJ

 

 

 

,

 

 

.

 

l

3

l

2

 

 

 

 

 

 

 

 

 

 

 

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

 

 

» 1 +

1

α

 

 

 

1 - Cosα »

1

α 2

 

1 + α

,

Sinα ≈ α ,

,

 

 

2

 

 

 

 

2

 

 

 

 

- l = l(

 

-1) »

l

(α / l)2

 

 

l 2 + α 2

1 + (α / l)2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

5.) Функция влияния при изгибе шарнирно опертого стержня поперечной силой P

в точке ξ будет

u(x,ξ ) =

Px(l x)

(l 2 - x 2 - (l - ξ )2 ) .

 

 

 

6lEJ

22

3. ПРИМЕР ВЫПОЛНЕНИЯ РАСЧЕТНО-ГРАФИЧЕСКОЙ РАБОТЫ «АНАЛИЗ КОЛЕБАНИЙ УПРУГОЙ СИСТЕМЫ С ДВУМЯ СТЕПЕНЯМИ СВОБОДЙ»

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

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

m1 = 4, m2 = 1, k1 = 1, k2 = 0, k = 1,

 

~

~

= 0,

 

f1 = 1, f2

X 01

= 1

X 02

= 0

V01

= 0

V02

= 1

23

1. Составим уравнения динамической модели рассматриваемой механической системы.

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

х1 (t), x (t) .

Кинетическая энергия системы состоит из энергии двух движущихся точечных масс

Т(хɺ1 , хɺ2 ) = m1 xɺ12 + m2 xɺ22 ,

2 2

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

П(х , х

 

) =

к1

x 2

+

к2

x 2 +

к

(х

 

х )2

 

2

 

 

 

2

 

1

2

1

2

2

2

 

 

 

1 .

 

 

 

 

 

 

 

 

 

 

 

 

Используя уравнения Лагранжа 2-го рода

 

L

 

d

(

L

) + Fi

= 0 , где L = T П

 

 

 

 

 

 

хi

dt

ɺ

 

 

 

 

 

 

 

 

 

 

 

 

xi

 

 

функция Лагранжа, а Fi обобщенные силы, работа которых не включена в потен-

циальную энергию и в нашем примере состоит из заданных внешних сил, имеющих

гармонический характер Fi =

~

× CosWt .

fi

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

 

 

..

k1x1

+ k (x2 x1 )+ F1 = 0

 

m1 x1

 

 

..

k

 

 

k (x x ) + F = 0

m x

x

2

 

2

2

2

 

2

1

2

или разрешив уравнения относительно старших производных с учетом вида внешних сил , получим

..

 

 

k1 + k

 

 

 

 

 

 

 

 

~

 

 

+

 

k

 

 

 

=

f1

 

cos Ωt

x1

x1

 

x2

 

 

 

m

m

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

~ 1

 

..

 

 

k

 

 

 

k

 

 

 

 

 

f 2

 

.

x

x +

 

x

 

=

 

cos Ωt

 

 

 

 

 

2

 

 

1

 

m2

2

 

 

m2

 

 

 

 

 

 

m2

 

 

 

 

 

 

Запишем уравнения модели 4-го порядка в индексной стандартной форме:

24

..

x1 + с11 x1 + с12 x2 = f1 cos Ωt

..

x2 + с21 x1 + с22 x2 2 = f 2 cos Ωt

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

 

 

k1 + k

 

 

 

 

 

 

 

 

 

~

с11

=

; с12

= −

k

 

f1

=

 

f1

 

m1

m1

 

m1

 

 

 

 

 

 

 

 

 

2 + k

 

 

 

~

 

= −

k

 

=

k

 

=

 

f 2

с21

 

;

с22

 

 

 

f 2

 

 

 

 

m2

 

m2

 

 

 

m2

 

 

 

 

 

 

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

x1 (0) = X 01 = 1 x2 (0) = X 02 = 0

.

(0) = V01

.

(0) = V02

 

x1

= 0 x2

= 1

Запишем уравнения модели так же а векторной форме

..

X + СX = f CosΩt ;

 

 

 

 

(0)

=

 

 

 

.

 

(0) =

 

 

 

 

X

Х0

;

 

X

V0

Где введены следующие обозначения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x (t )

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

X

=

 

 

X (t) =

1

 

 

 

 

 

 

01

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

0

X 02

 

 

 

x2 (t )

 

 

 

 

 

 

 

 

 

(1)

, (2)

 

 

 

V

 

 

V

 

 

0

=

01

 

 

 

 

 

 

 

 

 

V02

 

с11

с12

 

 

 

f

 

 

 

 

С =

 

 

 

f =

1

 

-

 

с21

с22

- приведенная матрица жёсткости системы;

 

 

f2

 

 

 

 

 

 

 

приведенные амплитуды внешних сил..

Общее решение уравнений системы состоит из двух частей

X (t) = X собст. (t) + X вынужд. (t)

собственных колебаний, соответствующих движению системы без воздействия внешних сил и вынужденных колебаний , обусловленных именно действием внешних сил.

2. Будем искать собственные (свободные) колебания системы в форме:

25

 

 

 

 

 

λt

 

= l1e

λt

 

 

 

 

 

 

 

 

 

 

 

 

 

Х

собств = Le

x1

 

 

;

 

=

= l2 e

λt

 

 

 

 

 

 

x2

 

 

подставляя его в соответствующие однородные уравнения модели, получим:

 

 

2

 

 

 

 

 

 

 

 

с11 + λ

с12

 

 

 

 

 

 

 

 

 

 

× L = 0

(3)

 

с21

 

с22 +

λ

2

 

 

 

 

 

 

 

 

 

для наличия ненулевых решений необходимо выполнение условия

(λ) =

с11 + λ2

с12

= λ4 + (с + с

 

2 + (с с

 

с с

 

) = 0

,

 

с21

с22 + λ2

11

22

11

22

12

21

 

 

называемого характеристически уравнением системы. При отсутствии трения оно

будет биквадратным. Решив это биквадратное уравнение

получим чисто мнимые

корни: λ1,2 = ±iω1; λ3,4 = ±iω2 , где ω1 < ω 2

 

собственные частоты системы.

Собственные вектора, соответствующие собственным частотам ω1 2 .находим

из решения уравнений (3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

1

 

ω

 

 

 

2 =

 

1

 

 

~ L1

 

 

~ L

 

1

=

 

;

2

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

χ1

 

 

 

 

 

 

 

 

 

χ2

 

 

 

χ1 =

 

 

с21

 

 

;

χ 2 =

 

 

с21

 

 

 

 

 

 

ω 2

с

22

 

 

ω 2

с

22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

В рассматриваемом примере коэффициенты жесткости будут соответственно равны:

с

 

=

k1 +k

= 1+1

= 1

2

; с

12

= −

k

 

= −

1

 

 

f

1

=

 

f1

= 0,25

 

 

 

m1

 

 

 

11

 

m1

4

 

 

 

 

 

 

4

 

 

 

 

 

m1

с

21

= −

 

k

= −1;

 

 

 

 

 

с

22

=

 

k

= 1 f

2

= 0

 

 

 

 

 

 

 

 

 

 

 

m2

 

 

 

 

 

 

 

 

 

 

 

m2

 

 

 

 

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

λ4 + 3 λ2 + 1 = 0 ,

2 4

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

ω1 = 1 2 3 - 5 0,437 и ω2 = 1 2 3 + 5 1,144 .

Собственные вектора получаем:

26

ω

 

 

 

 

 

1

 

 

ω

 

 

 

 

 

1

 

 

~ L1

 

 

~ L2

1

=

 

 

;

2

=

 

 

 

 

 

 

 

χ1

 

 

 

 

 

 

χ2

, где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

χ1

=

− 1

 

 

≈ 1,24; χ 2

=

 

− 1

 

≈ −3,25

0,437

2

− 1

1,44

2

− 1

 

 

 

 

 

Тогда решение для собственных колебаний имеет следующий вид:

Х cобст = (C1 cosω1t + C2 sin ω2t )× L1 + (C3 cosω1t + C4 sin ω2t )× L2

или

 

 

собст (t ) = A1

 

cos(ω1t - ϕ1 ) + A2

 

cos(ω2t - ϕ2 )

 

 

 

 

Х

L1

L2

 

 

 

Где С1,C2 ,C3 ,C4 произвольные константы, а A1 1 , A2 2

выражаемые через

них амплитуды и фазы

 

 

 

 

 

 

 

 

 

 

 

 

 

C2

A2 =

 

tgϕ2

=

C4

 

 

 

A1 =

 

 

 

tgϕ1 =

C32 + C42

 

 

 

C12 + C22

 

 

 

C3

,

 

 

C1

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

1 (t ) = A1

 

cos(ω1t − ϕ1 ),

 

2 (t ) = A2

 

cos(ω 2 t − ϕ 2 )

Х

L1

Х

L2

называются собственными формами колебаний на соответствующих частотах Запишем общее решение (однородной) системы уравнений модели для собст-

венных колебаний по координатно:

х1cобст = (C1 cosω1t + C2 sin ω2t ) + (C3 cosω1t + C4 sin ω2t )×

х2cобст = (C1 cosω1t + C2 sin ω2t )χ1 + (C3 cosω1t + C4 sin ω2t )χ2 ×

или через амплитуды и начальные фазы:

x1собст = A1 cos(ω1t − ϕ1 ) + A2 cos(ω2t − ϕ 2 ) x2собст = A1χ1 cos(ω1t - ϕ1 ) + A2 χ 2 cos(ω2t - ϕ 2 )

Неизвестные константы С1, С2 , С3 , С4 или выражаемые через них неизвестные ам-

плитуды А1 А2 и начальные фазы ϕ1 2 собственных колебаний определяются из начальных условий (2)

3. Решение для вынужденных колебаний системы находим в виде правой час-

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

 

 

 

 

 

 

 

 

 

 

ти так как

F = f CosΩt , то и

Х(t) = BCosΩt

;

 

 

, где В =

1

 

 

 

 

 

 

 

 

 

 

 

 

 

В2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

27

Здесь В1 , В2 амплитуды вынужденных колебаний, а Ω их частота.

Подставляя это решение в неоднородные уравнения (1) получим:

 

с - W2

с

 

 

 

B

 

f

 

 

 

11

 

 

12

 

 

1

 

 

 

1

 

 

с21

с11

- W

2

 

×

 

 

=

f

 

 

 

 

 

B2

 

 

2

Ненулевое решение этих уравнений для неизвестных константВ1 , В2 таково:

 

 

 

f1

 

с12

 

 

 

 

 

с11 − Ω 2

f1

 

 

 

 

 

 

 

 

 

 

 

 

 

B1

=

 

f 2 с22 − Ω 2

 

 

B2

=

 

с

21

 

f 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Ω)

 

 

 

 

(Ω)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B1

=

 

f (с − Ω 2 )f с

B2

=

 

f (с − Ω 2 )f с

 

 

1

22

 

2

12

 

2

 

11

 

1

21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Ω)

 

 

 

 

 

 

 

 

(Ω)

 

 

где главный определитель совпадает по виду с характеристическим уравнением

(Ω) =

с

− Ω2

 

с

 

= Ω4 − (с

+ с 2 + (с с

 

с с )

с

11

с

 

 

12

 

 

 

 

− Ω2

11

22

11

22

12

21

 

 

 

 

 

 

 

 

 

21

 

22

 

 

 

 

 

 

 

 

Частоты внешних сил, при которых наблюдаются неограниченное нарастание амплитуд вынужденных колебаний называются резонансными частотами. Условием резонанса является (Ω) = 0 , из чего следует, что собственные частоты в системах без трения одновременно будут и резонансными. Частота внешней силы называется антирезонансной, если амплитуды вынужденных колебаний равны нулю т.е. В1=0 или В2=0

10 )2 =

f1 с22 f 2 с12

;

02 )2 =

f 2 с11 f1с21

 

f2

 

f1

 

В рассматриваемом числовом примере находим амплитуды:

28

(Ω) =

с11 − Ω2

 

 

 

с12

 

=

12 − Ω2 - 14

 

= (12 − Ω2 )(1 - Ω2 )14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с21 с11 − Ω2

 

 

 

− 1 1 - Ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

1

 

 

с

12

 

 

 

 

1

 

 

 

- 1

 

 

 

 

 

 

 

 

с − Ω2

f

1

 

 

 

 

 

1

− Ω2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

2

 

 

 

 

 

B =

 

 

f 2

с22 − Ω2

 

 

=

 

 

0 1 - Ω2

 

 

 

 

B =

 

 

 

c21

 

f 2

 

 

=

 

 

-1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Ω)

 

 

 

 

 

 

(Ω)

 

 

 

 

(Ω)

 

 

 

 

 

 

 

(Ω)

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B =

 

f

1

(с

22

− Ω2 )f с

=

1 − Ω2

B =

f

2

(с − Ω2 )f

с

21

=

− 1

 

 

 

 

 

2 12

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

(Ω)

 

 

 

 

 

 

 

 

 

 

(Ω)

 

2

 

 

 

 

 

(Ω)

 

 

 

 

 

 

 

 

 

(Ω)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и одну из антирезонансных частот Ω10

= 1, Ω02

= ∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

мы, для чего составим таблицу амплитуд вынужденных колебании

для различных

значений частоты внешней силы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

1

2

3

4

5

6

7

8

9

10

11

12

 

13

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,00

0,11

0,23

0,34

0,43

0,46

0,57

0,69

0,80

0,92

1,03

1,14

 

1,26

1,37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆( )

0,25

0,23

0,17

0,09

0,00

-0,02

-

-0,23

-

-

-

0,00

 

0,38

0,98

 

 

0,13

0,30

0,30

0,22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B1

4,00

4,28

5,44

10,12

-

-

-2,25

-

-

0,28

 

-

-

 

 

38,95

5,03

1,19

0,53

 

1,52

0,91

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B2

4,00

4,34

5,74

11,47

-

-

-4,26

-

-

-

 

2,61

1,02

 

 

49,27

7,47

3,33

3,28

4,63

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На основании таблицы построим по тачкам графики АЧХ

29

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]