Моделирование колебательных процессов в машинах. Определение параметров виброизоляции
.pdf
МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ
Белорусский национальный технический университет
Кафедра «Теория механизмов и машин»
МОДЕЛИРОВАНИЕ КОЛЕБАТЕЛЬНЫХ ПРОЦЕССОВ В МАШИНАХ. ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ВИБРОИЗОЛЯЦИИ
Методическое пособие к лабораторным работам
Минск
БНТУ
2 0 1 5
МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ Белорусский национальный технический университет
Кафедра «Теория механизмов и машин»
МОДЕЛИРОВАНИЕ КОЛЕБАТЕЛЬНЫХ ПРОЦЕССОВ В МАШИНАХ. ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ВИБРОИЗОЛЯЦИИ
Методическое пособие к лабораторным работам
Рекомендовано учебно-методическим объединением по образованию в области машиностроительного оборудования
и технологий
Минск
БНТУ
2 0 1 5
УДК 60.001.11:531.8(076.5) ББК 34.41я73
М74
Авторы:
В.В. Кудин, Э.И. Астахов, М.В. Кудин, А.М. Авсиевич, А.А. Сухоцкий
Рецензенты:
В.М. Сурин, А.М. Тареев
Моделирование колебательных процессов в машинах. ОпредеМ74 ление параметров виброизоляции : методическое пособие к лабора-
торным работам / В.В. Кудин [и др.]. – Минск : БНТУ, 2015. – 38 с. ISBN 978-985-550-311-9.
В методическом пособии описаны схемы виброизоляции машин при силовом и кинематическом возмущениях, приведен расчет параметров виброизоляции и критерии ее эффективности.
Издание предназначено для студентов инженерно-технических специальностей вузов по дисциплине «Колебания в машинах», а также может быть использовано по отдельным работам в учебных курсах «Теория механизмов и машин», «Динамика машин» и другим родственным дисциплинам.
УДК 60.001.11:531.8(076.5) ББК 34.41я73
ISBN 978-985-550-311-9 |
© Белорусский национальный |
|
технический университет, 2015 |
1 ОСНОВНЫЕ ТЕОРЕТИЧЕСКИЕ ПОЛОЖЕНИЯ ВИБРОИЗОЛЯЦИИ МАШИНЫ
Виброзащита машин – это совокупность методов и средств уменьшения вибраций машины от внешних или внутренних источников вибровозмущения. Одним из распространенных методов виброзащиты является виброизоляция, заключающаяся в том, что между защищенным объектом 1 и источником вибровозмущения 2 устанавливается упругий элемент – виброизолятор 3 (Рис.1).
а) |
б) |
|
Рис. 1 |
Различают две схемы виброизоляции: а) с динамическим или силовым вибровозмущением (рис. 1, а), когда источник вибровозмущения 2 создает переменную силу F(t) , передаваемую через виброизолятор 3 на защищаемый объект 1; 2) с кинематическим вибровозмущением (рис. 1, б), когда вибровозмущение z 2 (t) источника 2
через виброизолятор 3 передается на защищаемый объект 1. Виброизолирующее устройство 3 представляет важнейшую
часть виброзащитной системы, назначение которой состоит в создании такого режима движения, инициируемое заданными возмущениями, при которых реализуется цель виброзащиты объекта.
3
Параметры виброизолятора, как и любого упругого звена, определяются по его динамической характеристике Rz (z, z) , представляющей зависимость силы реакции Rz виброизолятора от его деформации z и скорости z деформирования. В общем случае динамическая характеристика R(z, z) является нелинейной функцией.
Однако при анализе малых колебаний объекта вблизи положения равновесия можно производить линеаризацию и ограничится ли-
нейной зависимостью реакции Rz от перемещения z и скорости z .
Rz cz z bz z |
(1) |
Виброизолятор с линейной характеристикой (1) называется линейным, параметры которого ( сz – коэффициент жесткости, bz – коэффициент сопротивления) являются постоянными. В случае bz 0 зависимость (1) описывает характеристику линейного идеального упругого элемента (пружины), а при cz 0 – характери-
стику линейного вязкого демпфера. Таким образом, модель виброизолятора 3 с характеристикой (1) можно представить в виде па-
раллельного соединения пружины с c z и демпфера с bz . Рассмотрим эти две схемы виброизоляции объекта.
2 ВИБРОИЗОЛЯЦИЯ С ДИНАМИЧЕСКИМ (СИЛОВЫМ) ВИБРОВОЗМУЩЕНИЕМ
2.1 Малые колебания в системе с шестью степенями свободы
Объект, как твердое тело (корпус станка, прибор и т.п.) установленный на упругих амортизаторах , опирающихся на неподвижное основание (рис. 2), представляет собой колебательную систему, обладающую шестью степенями свободы и шестью частотами собственных колебаний, которые во избежание резонанса с частотами
4
возмущения сил F(t) должны быть заранее определены расчетным путем.
|
Рис. 2. |
|
Малые колебания рассматриваемой системы характеризуются: |
а) |
тремя линейными смещениями центра масс объекта xs , ys , zs ; |
б) |
1 , 2 , 3 – тремя углами поворота относительно осей x, y, z |
неподвижной системы координат с началом в центре масс S. При этом упругие свойства виброизолятора характеризуются соответственно коэффициентами жесткости cx , cy , cz в направлении осей x, y, z.
Если пренебречь массой виброизоляторов, их демпфирированием и гироскопическими эффектами, возникающими при колебаниях, то дифференциальные уравнения собственных колебаний такой системы можно записать при помощи уравнения Лагранжа 2-го ро-
да [1].
Потенциальную энергию системы как функцию координат смещения центров жесткости виброизоляторов представится в виде
5
|
1 |
2 |
2 |
|
2 |
|
|
U |
|
( cx 1 |
cy 1 |
cz |
1 ), |
(2) |
|
2 |
|||||||
|
|
|
|
|
|
||
где 1 , 2 , 3 – составляющие смещения центра Ok |
жесткости |
||||||
виброизолятора, имеющего координаты x, y, z. |
|
|
|||||
Воспользуемся известными кинематическими соотношениями между смещением произвольной точки и центра масс, тогда пере-
мещения |
, |
, |
3 |
центра |
жесткости |
i-го |
виброизолятора с |
|
1 |
2 |
|
|
|
|
|
начальными координатами x,y,z |
|
|
|||||
|
1 |
xs |
|
3 y |
2 z – по оси |
x , |
|
|
2 |
уs |
1z |
3 x – по оси |
y , |
(3) |
|
|
3 |
zs |
|
2 x |
1 y – по оси |
z. |
|
Подставив (3) в уравнение (2), получим выражение для потенциальной энергии системы
U |
|
1 |
( cx (xs |
3 y |
|
2 z)2 |
|
c y ( ys |
1 z |
3 x)2 |
|
||||
2 |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
2 |
1 |
|
|
2 |
|
2 |
2 |
2 |
|
|
|
cz (zs |
2 y |
1 y) |
|
|
|
( |
11 xs |
12 ys |
33 zs |
44 1 |
|||
|
|
|
2 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
55 32 ) |
15 xs 2 |
16 xs 3 |
|
24 ys |
3 |
34 zc 1 |
|
|||||||
|
35 zc 2 |
45 1 2 |
|
56 2 3 |
|
46 1 3 , |
|
|
|||||||
где
– суммирование производится по всем виброизоляторам:
11 |
cx , 22 |
cy , |
33 |
cz , |
44 |
|
(cz y2 |
cy z 2 ), |
|
55 |
(c z2 |
c x2 ), |
66 |
(c y2 |
c |
y |
x2 ), |
|
|
x |
z |
x |
|
|
|
|
|||
16 |
cx y, |
24 |
cy z, |
26 |
cy x, |
34 |
cz y, |
||
35 |
cz x, |
45 |
cz x y, 46 |
|
|
|
cy z x, |
||
56
cx
y
x.
6
Кинетическая |
|
|
энергия |
системы |
Т ТП ТВ , |
где |
||||||
|
m(x2 |
y2 |
z |
2 ) |
|
|
|
|
|
|
|
|
ТП |
s |
s |
|
s |
|
– кинетическая энергия при поступательном |
||||||
|
2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
движении центра масс; Т |
|
I |
2 |
– кинетическая энергия при вра- |
||||||||
|
|
|
||||||||||
В |
2 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
щательном движении системы относительно мгновенной оси, проходящей через центр масс объекта; I – момент инерции объекта относительно мгновенной оси;
– мгновенная угловая скорость вращения; m – масса объекта.
Обозначим углы 1, 2 и 3 – углы образованные мгновенной
осью вращения с осями x, y, z.
Проекции угловой скорости на оси координат
|
|
|
|
1 |
|
|
|
|
|
cos |
|
1 , 2 |
|
|
|
cos |
|
2 , |
|
3 |
|
|
|
|
cos |
3 . |
|
|
|
|
||||||||||
Момент инерции относительно мгновенной оси вращения |
|
|
|
|
||||||||||||||||||||||||||||||||||||
I |
|
|
I x cos2 |
|
|
1 |
I y cos2 |
|
|
2 |
|
I z cos2 |
|
3 |
2I xy cos |
1 cos |
|
2 |
||||||||||||||||||||||
|
2I yz cos |
|
|
2 cos 3 |
2I xz |
cos |
3 cos |
1 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
где Ix , I y , Iz , Ixy , I yz , Ixz |
– моменты инерции и центробежные мо- |
|||||||||||||||||||||||||||||||||||||||
менты инерции объекта относительно осей x, y,z. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Т |
|
|
|
1 |
(I |
|
|
2 |
I |
|
2 |
I |
|
|
2 ) |
I |
|
|
|
|
|
|
I |
|
|
|
|
|
I |
|
|
|
|
|
|
||||
|
В |
|
|
x |
y |
z |
xy |
1 |
2 |
xz |
1 |
3 |
yz |
2 |
3 |
|||||||||||||||||||||||||
|
|
2 |
|
|
1 |
|
2 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Кинетическая энергия объекта |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
T |
|
1 |
|
|
(x |
2 |
|
y |
2 |
|
z2 ) |
1 |
( |
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
2 |
11 |
|
|
s |
|
|
s |
|
s |
2 |
|
44 |
1 |
|
|
|
55 |
|
2 |
|
66 |
|
3 |
|
45 |
|
|
1 |
|
2 |
||||||||
|
56 2 3 |
|
|
|
46 3 1, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7
11 |
m, 44 |
Ix , 55 |
I y , 66 |
Iz , |
45 |
Ixy , 56 |
I yz , |
46 |
Ixz . |
С помощью уравнения Лагранжа 2-го рода получим шесть взаимосвязанных дифференциальных уравнений собственных колебаний объекта виброзащиты:
11 xs |
11 xs |
15 |
2 |
16 |
3 |
11 ys |
22 ys |
24 |
1 |
26 |
3 |
11 zs |
33 zs |
34 |
1 |
35 |
2 |
0
0
0
(4)
44 |
1 |
45 |
2 |
46 |
3 |
24 ys |
34 zs |
|
44 |
1 |
45 |
2 |
46 |
3 |
|
0 |
45 |
1 |
55 |
2 |
56 |
3 |
15 xs |
45 1 |
45 |
1 |
55 |
2 |
56 |
3 |
|
0 |
|
46 |
1 |
56 |
2 |
66 |
3 |
16 xs |
26 ys |
|
46 |
1 |
56 |
2 |
66 |
3 |
|
0 |
Решение системы (4) запишем в виде |
|
|
|
|
|
|
|
|
||||||||
xs |
A1 cos( |
t |
|
), ys |
A2 cos( |
t |
|
), zs |
|
A3 cos( |
t |
|
), |
|||
|
1 |
A4 cos( |
t |
|
), 2 |
A5 cos( |
t |
|
), |
3 |
A6 cos( |
t |
), |
|||
где A1, A2, A3, A4, A5, A6 – значение амплитуд изменения |
соответ- |
|||||||||||||||
ствующего параметра;
k – круговая частота собственных колебаний;
– начальная фаза колебаний.
Подставляя эти решения в систему (4) и предполагая, что система допускает решения отличное от нуля, то определитель системы будет равен нулю:
8
|
|
2 |
0 |
|
0 |
|
0 |
|
|
|
|
|
|
11 |
|
11 |
|
|
|
15 |
|
16 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
|
0 |
|
|
2 |
0 |
|
|
|
0 |
|
|
|
|
|
22 |
|
11 |
|
24 |
|
|
26 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
|
0 |
|
0 |
|
|
2 |
|
|
|
|
0 |
|
|
|
|
33 |
|
33 |
2 |
|
2 |
|
0 |
(5) |
|||
|
|
|
|
|
|
|
2 |
||||||
|
|
|
|
|
34 |
|
35 |
|
|
|
|
||
|
0 |
|
24 |
|
34 |
44 |
44 |
45 |
45 |
46 |
46 |
|
|
|
|
|
|
|
|
||||||||
|
|
|
0 |
|
|
|
2 |
|
2 |
|
2 |
|
|
|
15 |
|
|
35 |
45 |
45 |
55 |
55 |
56 |
56 |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
0 |
|
2 |
|
2 |
|
2 |
|
|
|
16 |
|
26 |
|
46 |
46 |
56 |
56 |
66 |
66 |
|
|
|
|
|
|
|
|
|
Развернув определитель (5), получим уравнение шестой степени
относительно 2 , решение которого возможно только численным методом с помощью ЭВМ. Рассмотрим некоторые частные случаи, важные в практическом отношении.
2.2 Объект устанавливается без перекоса виброизоляторов, имеющие одинаковые жесткости с z
Необходимо стремиться к тому, чтобы в состоянии покоя объект виброизоляции не имел перекосов, т.е. чтобы осадка всех виброизо-
ляторов была одинаковой. Это возможно, если |
34 |
35 |
0, |
т.е |
|
|
|
когда центр масс объекта (т. S) и центр жесткости виброизолятора (т. Oк) лежат на одной вертикальной линии.
Если коэффициенты жесткости сz всех виброизоляторов выбра-
ны одинаковыми, то третье уравнение в (4) станет независимым от других, т.е.
11zs
33zs 0 .
Следовательно, можно определить одну из собственных частот в направлении оси z
1
33 .
11
9
