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

LS-Sb88870

.pdf
Скачиваний:
10
Добавлен:
13.02.2021
Размер:
742.68 Кб
Скачать

 

 

 

 

 

P

=

 

1

 

0,0322 = 203 Н м 2 .

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходные данные для несинусоидальности напряжения:

 

I1= 5 А; I5= 3 А; I7 = 2 А; I11 = 1 А; I13 = 0,5 А.

 

 

Амплитуды основной волны магнитной индукции:

 

 

 

B = 1,35

 

 

10 0,9

 

 

4π 10 7 · 5 = 0,0318 · 5 = 0,159 Тл.

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2 1,2 0,0002

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично рассчитываются амплитуды магнитной индукции:

 

 

 

B5 = 0,0318 · 3

= 0,096 Тл; B7 = 0,0318 · 2 = 0,064 Тл;

 

 

B11 = 0,0318 · 1 = 0,032 Тл; B13 = 0,0318 · 0,5 = 0,016 Тл.

 

Амплитуды силовых волн:

 

 

 

 

 

 

 

 

 

 

 

 

P =

1

 

B2 =

 

107

 

0,1592 = 5032

Н м 2 ; P =

 

1

 

B2 = 1834 Н м 2

;

0

 

4 4π

 

0

1

1

 

 

 

 

 

 

 

 

 

 

 

 

5

 

5

 

 

P =

1

 

B2

= 815 Н м 2 ; P

=

1

B2

= 203 Н м 2 ;

 

 

 

0

 

0

 

 

7

 

 

 

 

7

 

 

 

 

 

 

 

 

 

11

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P =

 

 

1

 

B2

= 50 Н м 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

 

 

0

13

 

 

 

 

 

 

 

 

 

Значение

угловой

частоты:

ω1=

2 ω0 ;

ω5

=

5 ω0 ;

ω7 = 7 ω0 ; ω11=

= 11 ω0 ; ω13 = 13 ω0 . Исходные данные для расчета приведены в табл. 6.1.

Таблица 6.1

Величина

Размерность

 

 

 

 

Вариант

 

 

 

 

1

2

3

4

 

5

 

6

7

8

9

I01

 

А

10

10,5

11

9

 

9,5

 

10,5

11

10

9

I02

 

А

1,0

1,5

2,0

1,7

 

1,8

 

1,4

1,3

1,2

1,1

I1

 

А

5,0

5,1

5,2

5,3

 

5,4

 

5,5

4,9

4,8

4,7

I5

 

А

3,0

3,1

3,0

2,9

 

3,0

 

3,1

2,9

3,1

2,9

I7

 

А

2,0

2,1

2,2

2,1

 

2,0

 

1,9

2,0

2,1

2,2

I11

 

А

1,0

1,0

1,0

1,1

 

1,0

 

0,9

1,0

1,1

0,9

I13

 

А

0,5

0,5

0,6

0,5

 

0,6

 

0,6

0,5

0,4

0,4

Для всех вариантов

практического занятия № 6 принимаются одинако-

вые значения для

следующих параметров: w = 10; p = 2; δ = 0,2 мм;

Kоб =

= 0,9; Kc = 1,2.

 

 

 

 

 

 

 

 

 

 

 

 

31

Практическое занятие 7 РАСЧЕТ ОЖИДАЕМОГО УРОВНЯ ВИБРАЦИИ АД,

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

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

Исходные данные. Амплитуда пульсирующего момента M(t), Н·м; амплитуда радиальной электромагнитной силы, действующей на единицу площади, P, H; угловая частота пульсирующего момента ωm , рад/с; угловая ча-

стота радиальной силы ωp , рад/с.

Массогабаритные характеристики АД: масса статора Мst , кг; модуль

упругости материала статора Е, Н/м2 ; радиус внутреннего статора Rв, м; радиус наружный статора Rн, м; высота спинки статора h, м; длина статора l, м; средний радиус спинки статора Rc , м.

M(t)

J

Cτ

Ps (t)

MSt

Cr

Рис 7.1

Характеристики крепления АД к неподвижному фундаменту: коэффициенты жесткости виброизоляторов в направлении осей Y и Z – соответственно Cy и Cz , H/м; расстояние между амортизаторами Rа, м; ко-

32

эффициент механического сопротивления поступательному движению Rr ; количество амортизаторов n; угол между направлением действия силы с осью Z – ψ, … .

Требуется найти : Амплитуды колебания АД от действия приложенных момента и силы Zm , Z p и соответствующие уровни колебаний Lm , Lp .

Алгоритм расчета:

Расчет ожидаемых уровней колебаний от действия пульсирующего момента M(t). Момент инерции статора вместе с корпусом и присоединенными массами:

J = 0,5 Мst ( Rв2 + Rн2 ).

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

 

 

 

 

Cа =

CyCz

,

Czsin ψ 2

Сycos ψ 2

где Cy = nCy , Cz = nCz .

Коэффициенты жесткости и сопротивления виброизоляторов крутильным колебаниям соответственно:

C = Cа Rн2 ; K = α C J ,

где α = 0,1…0,2 – коэффициент демпфирования.

Амплитуда угловых виброперемещений статора АД при действии на не-

го, пульсирующего момента:

M t

γ =

 

.

C Jω2 2 K ω 2

Амплитуда линейных виброперемещений лапы АД в направлении оси Z :

Zm = 0,5 Rа γcos ψ.

Уровень ожидаемой вибрации лапы АД, возбуждаемой пульсирующем моментом на частоте ωm , дБ:

Lm = 20 lg

Zmωm2

.

2

3 10 4

 

 

Расчет ожидаемых уровней вибрации от действия радиальных электромагнитных сил. Колебания статора рассматриваются в плоскости сечения

33

кольца и характеризуются порядком колебаний r. Коэффициент жесткости кольца статора:

Cr =

r2

(r2 1)2 Elh2

r 2.

 

 

12R3

,

 

r2 1

 

 

 

c

 

 

Амплитуда радиальной электромагнитной силы, приложенной к полной поверхности статора:

Ps = 2π RвlP.

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

Z p =

Ps

 

.

Cr Мstω2p 2

 

 

Rrωp 2

Уровень колебаний лапы АД от действия радиальной силы на частоте ωp :

Lp = 20 lg

Z pω2p

.

2

3 10 4

 

 

Пример расчета. Данные по вынуждающим силам и моментам взяты из предыдущих расчетов. Численные значения исходных данных следующие:

M(t) = 2,03 H м; ωp = ωm = 2π 100 рад/с ; P = 20128 Н м 2 ; h = 0,018 м;

l = 0,08 м; Rc = 0,07 м; Мst = 40 кг; Е = 2

1011 Н м 2 ; Rв = 0,05 м; Rн =

= 0,1 м;

 

C

z

= 0,65 106 Н м 1;

C

y

= 0,35

106 Н м 1;

R = 0,2 м; n = 4;

 

 

 

 

 

 

 

 

 

 

а

r= 2p = 4;

 

Rr = 0; ψ = 0.

 

 

 

 

 

 

 

 

Расчет ожидаемых уровней вибрации от действия M(t).

 

J = 0,5

40( 0,052 0,12 ) = 0,25 Н м с2 ; C = 4

 

0,65 106 = 2,6 106 Н м 1;

 

 

 

 

 

 

 

 

z

 

 

 

 

Cа= Cz = 2,6 106 Н м 1; C

= 2,6

10 2

106 = 2,6

104 H м;

 

 

 

 

K = 0,1 2,60 104

0,25 = 8,1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ =

 

 

 

2,03

 

 

 

 

 

 

= 0,27 10 4 рад;

 

8,6 104 0,25 2π 100 2 2

 

 

 

 

 

 

8,12

100 2

 

34

Z

m

= 0,1 0,27

10 4 = 0,27

10 5 м; L

= 20 lg

0,27 10 5 2π 100 2

= 68 дБ.

 

 

 

 

 

 

 

m

 

 

 

2 3 10 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

Cr =

42(42 1)2

2π 2 1011 0,08 0,0183

= 2,87 1010 Н м 1;

 

 

 

42 1

 

 

12 0,073

 

 

 

 

Ps

= 2π 0,05 0,08

P = 2π 0,05 0,08 20128 = 505,6 H;

 

 

 

 

Z p =

 

 

505,6

 

= 1,76 10 8 м;

 

 

 

 

 

1010

 

 

 

 

 

 

2,87

40(2π 1002)

 

 

 

 

 

Lp = 20 lg

2π 100 21,76 10 8

= 24 дБ.

 

 

 

 

 

2 3 10 4

 

 

 

 

 

 

 

 

 

 

Варианты для расчета.

Данные по АД одинаковые для всех вариантов, а значения вынуждающих моментов и сил можно взять из предыдущих расчетов. Необходимо провести исследование влияния исходных параметров на уровень вибрации (например, масса и габаритные размеры АД, частоты и т.д.).

Практическое занятие 8 ВЫБОР И РАСЧЕТ ЭФФЕКТИВНОСТИ АМОРТИЗАТОРОВ

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

Исходные данные. Для полного расчета необходимо знать: массы отдельных частей агрегата Мi , кг; координаты центра тяжести этих частей относи-

тельно выбранной системы координат; неуравновешенные силы и моменты, действующие на агрегат, P, M; их частоты ωp , ωm , а также плечи, на кото-

35

рых приложены силы l, d, м. Кроме того, нужны справочные данные по амортизаторам: номинальная статическая нагрузка на сжатие Pн, статиче-

ские и динамические (для приближенного расчета) жесткости по осям Cxaj , Czaj ,Cyaj , которые располагаются согласно фронтальной проекции.

Требуется найти. Частоты собственных колебаний: поступательного движения вдоль осей fx , f y , fz и поворотные fвx , fвy , fвz , а также ча-

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

Алгоритм расчета. Находим центр тяжести всего агрегата, через который проводим новые оси системы координат XOY (на одной платформе, которая крепится на амортизаторах, может быть установлено несколько механизмов). Определим координаты центра тяжести отдельных частей агрегата и рассчитаем массу всего агрегата М = Mi . Заменяя отдельные части агрега-

та параллелепипедами, определим моменты инерции этих частей относительно новых осей и, суммируя, вычислим момент инерции массы всего аг-

регата: J = Ji . Крутящие моменты относительно осей X и Y

запишем в ви-

де: M x Pzd cosωpt ;

M y Pzl cosωpt , где Pz , Pzd , Pzl

амплитуды воз-

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

чество: Nmin М/Pн, где ( f0 / fz ) 2 номинальная статическая масса.

В зависимости от конструкции крепления агрегата к фундаменту выбираем фактическое число амортизаторов Nф Nmin , а также координаты установки

амортизаторов x j , y j , z j и угол их наклона θ. Динамические и статические

жесткости амортизаторов рассчитываются по формулам:

Czj = Cyaj cos2 θ+Cxaj sin2 θ;

Cyj = Cxaj cos2 θ+Czaj sin2 θ;

Cxj = Cyaj cos2 θ+Czaj sin2 θ.

36

Если угол наклона находится в плоскости ZOX, для всех амортизаторов

Cxj = Cxaj .

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

n

n

n

Cz = Czj

;Cy = Cyj

;Cx = Lви Cxj .

j 1

j 1

j 1

Поворотные жесткости рассчитываются по формулам:

n

n

n

Kx =

(Cyj Czj ) ; K y =

(Cxj Czj ) ; Kz = (Cxj Cyj ) .

j 1

j 1

j 1

При расчете зависимости частоты колебаний агрегата от действия неуравновешенных сил и моментов используются динамические жесткости амортизаторов; смещения от действия качки судна определяются по статическим жесткостям. Частота поступательных и поворотных колебаний вдоль и относительно оси Z, в Гц, находится следующим образом:

fz

1

Cz

;

fвz

1

Kz

.

М

Jz

 

 

 

 

Аналогично определяются частоты по осям X и Y. Связанные свободные колебания сдвига и поворота в плоскостях ZOY и ZOX могут быть опреде-

лены по формулам:

 

f

x1

 

 

 

1

 

1

C

x

 

K y

 

1

C

x

 

 

K y

2

C2Z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

;

 

fx2

 

 

 

2

 

М

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

J yМ

 

 

 

 

 

 

 

J y

 

 

М J y

 

 

 

f y1

 

 

 

1

 

1

Cy

 

K

x

 

 

 

1

Cy

 

 

K

x

 

 

2

Cy2Z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f y2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

J yМ

 

 

 

 

 

 

 

М J x

 

 

М J x

 

 

Амплитуда вынужденных колебаний агрегата вдоль оси Z = Zц. т для одномассовой системы определяется выражением:

 

P

 

f 2

 

 

P

Z

z

 

z

 

; =

z

,

 

 

 

 

fz2 f p2

 

Cz

 

 

Cz

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

37

приближенную формулу 250 . fz2

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

L

 

 

f p2

 

 

= 20 lg

1

 

.

 

ви

 

 

fz2

 

 

 

 

 

 

При f p fz растут амплитуды вынужденных колебаний, следовательно, действие амортизаторов имеет отрицательный эффект. Положительный эффект имеет место, если ( f p / fz ) > 2 . На практике, целесообразно, чтобы выполнялось условие ( f p / fz ) = 2,5…5. Эффективность поглощения энергии амортизации, %, определяется по формуле:

 

 

 

 

 

 

f02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П 1

 

 

 

 

 

 

 

 

 

100 ,

 

 

 

 

 

f 2

 

 

f 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

p

 

 

 

 

 

 

 

 

где f0 частота собственных колебаний (табл. 8.1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 8. 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эффективность

Относительная частота

f p

 

поглощения

 

 

f0

 

 

энергии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,4

 

2,0

 

 

3,0

 

4,0

 

5,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П, %

0

 

67

 

 

87,5

 

93

 

96

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

например, с частотой ходовой вибрации. При этом выполняется следующее

 

f pв

 

 

 

n

 

 

условие:

 

0,8, где f

pв

 

в

( n

число оборотов винта). В этом

fmin

 

 

 

60

в

 

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

38

с первой гармоникой частоты возмущения f p1. При выполнении условия

f p1

1,3 больших амплитуд колебаний на этой частоте не будет.

fmax

 

 

Пример расчета. Рассмотрим алгоритмы приближенного расчета, ко-

торый может быть рекомендован как ориентировочный для одномассовой

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

Zц. т= = 0,7 0,3 мм; Xст = 0,2 0,4 мм.

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

Z. Далее анализируем

отношение f p / fz

для частоты возмущения f p ,

соответствующей частоте питания агрегата.

Если это отношение меньше

двух, увеличиваем Zц. т

и добиваемся приемлемого значения f p / fz . Затем

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

оси

Cz по формуле

Cz / Zц. т. По справочнику, в соответствии со значение

Zц.т

выбираем Czj

и определяем минимальное количество соответствую-

щих амортизаторов Nmin Cz/Czj .

Пример расчета. Положение центра тяжести всего агрегата и направление осей считаем известными. Остальные данные следующие: М = 2280 кг;

J

x

= 260 кг м2 ;

J

y

= 1070 кг м2

; J

z

= 930 кг м2 ; P = 16000 H;

 

 

 

 

 

z

fp = 100 Гц.

 

 

 

 

 

 

 

 

Выберем для установки амортизатор типа АКСС-300И. Тогда минималь-

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

Nmin

2280 = 8. Справочные данные

 

 

 

 

 

 

 

 

300

по амортизаторам даны в табл. 8.2. По конструктивным соображениям принимаем Nф = 10; координаты расположения амортизаторов следующие

x1,2 = 900 мм; x3,4 = 700 мм; x5,6 = 500 мм; x7,8 = 200 мм; x9,10 = 0. По оси

Y амортизаторы расположены симметрично, поэтому y1,3,5,7,9 y2,4,6,8,10 = = 320 мм. Возвышение центра тяжести агрегата над плоскостью расположения амортизаторов принимаем z = 645 мм.

Динамическую жесткость рассчитываем при θ = 20º:

39

Czj = 395

104 0,942 + 270

104 0,342

= 37,7

105 H м 1 ;

Cyj = 270

104 0,942 + 395

104 0,342

= 28,4

 

105 H м 1 ;

 

Cxj = Cyaj =

11

105 H

м 1 ;

 

 

Cz = 10Czj = 377 105 H м 1 ; Cy

= 10 Cyj = 284

105 H м 1 ;

 

Cx = 10 Cxj = 110

105 H

м 1 .

 

 

Поворотные динамические жесткости рассчитываются по формулам:

Kx = (284 104 0,6452 + 377

104

0,322 ) = 156,5

105 H м 1 ;

K y = 377

104 ( 0,22 + 0,52 + 0,72 +0,92 ) 2 + 110 105

0,6452 =

 

= 165,6

105 H

м 1 ;

 

Kz = 110

104 0,322 10 + 284

10 4 (0,22 + 0,52 + 0,72 +0,92 ) 2 =

 

= 101

105 H

м 1 .

 

Частота поступательных и поворотных свободных колебаний относи-

тельно оси Z:

 

 

 

 

 

 

 

fz

1

377 105

= 20,5 Гц;

fвz

1

101 105

= 16,6 Гц.

2280

930

 

 

 

 

Связанно-свободные частоты колебания сдвига и поворота, вычисленные по приведенным выше формулам, равны:

 

f

x1

 

 

20,9

 

Гц;

 

f y1

 

 

 

42,1

Гц.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fx2

 

8,9

 

 

f y2

 

 

8,1

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, частоты собственных колебаний равны:

fz = 20,5 Гц;

fx1 = 20,7 Гц;

f y1 = 42,1 Гц;

fвz =16,6 Гц;

fx2 = 8,9 Гц;

f y2 = 8,1 Гц.

Амплитуда вынужденных колебаний агрегата:

 

P

 

f 2

 

= 1,1 10 4 м.

 

 

 

Zц. т = Z

z

 

z

 

Cz

fz2 f p2

 

 

 

 

 

40

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