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

Влияние особенностей конструкционного демпфирования на нелинейные колебания виброзащитных систем Текст научной статьи по специальности «Физика»

CC BY
117
35
Поделиться
Ключевые слова
ВИБРОЗАЩИТНЫЕ СИСТЕМЫ / ГИСТЕРЕЗИС / АППРОКСИМАЦИЯ / УРАВНЕНИЕ ДВИЖЕНИЯ / КВАЗИ-ГАРМОНИЧЕСКИЕ КОЛЕБАНИЯ / VIBROPROTECTIVE SYSTEMS / HYSTERESIS / APPROXIMATION / EQUALIZATION OF MOTION / QUASIHARMONIC VIBRATIONS

Аннотация научной статьи по физике, автор научной работы — Лазуткин Г. В., Антипов В. А., Рябков А. Л.

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

Похожие темы научных работ по физике , автор научной работы — Лазуткин Г. В., Антипов В. А., Рябков А. Л.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

INFLUENCE OF FEATURES OF CONSTRUCTIONAL SUPPRESSION OF VIBRATION ON NONLINEAR FLUCTUATIONS OF SYSTEMS AT DISPERSION OF ENERGY

The task about the elementary form of approximation of change of the characteristics of devices on dispersion of energy with a constructional hysteresis is decided(solved) on the basis of minimization of a root-mean-square error of approach(approximation), and also equality of the areas of close loops of a hysteresis to required. With the help of the combined method of decomposition of the decision of the nonlinear differential equation of movement of systems at dispersion of energy and offered form of approximation the decision for approximate repeating fluctuations taking into account influence of features of isolation of a constructional hysteresis and overwhelming devices on the basic parameters of fluctuations is found.

Текст научной работы на тему «Влияние особенностей конструкционного демпфирования на нелинейные колебания виброзащитных систем»

УДК 620.318.6

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

© 2009 Г.В. Лазуткин, В.А. Антипов, А.Л. Рябков

Самарский государственный университет путей сообщения

Поступила в редакцию 05.12.2007

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

Системы конструкционного демпфирования (СКД) широко применяются при создании высокоэффективных виброзащитных систем (ВС) агрегатов и систем силовых установок, подвижного состава, радиоэлектронного оборудования, эксплуатируемых в условиях воздействия на них интенсивных динамических нагрузок. Однако, несмотря на многолетний опыт разработки и эксплуатации подобных ВС, до сих пор возникают трудности при создании методов расчета их динамических характеристик даже при простом гармоническом виде возбуждения. С одной стороны, это связано с неправомерностью традиционных подходов, основанных на предположении о гармонической, а не квазигармонической форме колебаний ВС [1]. С другой - отсутствием адаптированных к особенностям конструкционного демпфирования приближенных аналитических методов решения нелинейных дифференциальных уравнений движения ВС и корректного математического описания упругогистерезисных характеристик (УГХ) СКД.

Сложившийся подход к исследованию УГХ виброизолирующих и демпфирующих устройств, представляющих собой сложные СКД, основан на теоретическом изучении их простейших моделей [2, 3]. Они отображают основные особенности упругого и фрикционного взаимодействия с сухим трением множества контактирующих элементов, входящих в конструкцию устройств и имеющих определенную, в общем

Лазуткин Геннадий Васильевич, кандидат технических наук, докторант. E-mail: gamaba@list.ru. Антипов Владимир Александрович, доктор технических наук, профессор

Рябков Алексей Леонидович аспирант

случае, нелинейную геометрическую форму. Указанные обстоятельства позволяют выявить вид основных определяющих критериев подобия и с помощью методов теории подобия и анализа размерностей экспериментально определить обобщенные деформационные характеристики (УГХ) для изучаемого класса СКД.

Обычно УГХ устройств получают при их циклическом деформировании. Найденные при этом семейства петель гистерезиса преобразуют к обобщенному (безразмерному) виду с помощью афинных преобразований координат нагрузка-деформация. Причем для описания процессов деформирования произвольного нагружения СКД, свойственного переходным процессом неустановившихся колебаний ВС, определяют виды нелинейных преобразований процессов нагрузки и разгрузки петель гистерезиса, а также условия их сшивки в точках смены знака скорости деформации [2, 3]. При наличии полученных таким образом данных об УГХ устройств несложно получить известными численными методами (Рунге-Кутта и др.) решение дифференциального уравнения движения ВС, обладающее свойством точного. Последнее особенно важно при проведении поверочных расчетов динамических характеристик ВС.

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

стерезиса СКД. Оно должно достаточно полно и просто отображать важнейшие особенности УДХ устройств, связанные с нелинейностью их упругих и неупругих составляющих.

Представим описание контуров симметричных петель гистерезиса суммой двух функцией: ф(8, х, А) = Фу( х, А) + Фш(8, х, А), (1) где х, X, А - соответственно, скорость, перемещение (деформация) и ее амплитуда;

Фу=0,5(Ф1 + Ф2);Фт=0,5(Ф1- Ф2) - соответственно условные упругая и неупругая составляющие реакции СКД - И = Ф(8, х, А);

Ф1, Ф2 - функции, описывающие процессы нагрузки (X > 0) и разгрузки (X < 0) соответственно;

а = signx - ступенчатая функция, определяемая знаком скорости.

Воспользовавшись результатами анализа работ [2, 3, 4], отметим главные особенности поведения УГХ для различных СКД, выполненных на основе материала МР, тросов, многослойных пластин и т.п.:

- нелинейность функций Фу и Фт в наибольшей степени проявляется на концах отрезков X е [- А, А],где коэффициенты жесткости могут изменяться в десятки и сотни раз;

- форма неупругой составляющей Фт (гис-терезисной функции) в зависимости от амплитуды деформации А может быть близкой как к эллипсу, так и к параллелограмму с вертикальными стенками, характерному для сухого трения;

- упругая и неупругая составляющие имеющие сдвиг по фазе, между собой, равный п / 2.

Как показано в работе [2], аппроксимация упругой составляющей реакции СКД ортогональными полиномами Чебышева обеспечивает наилучшее приближение на концах отрезка X е [- А, А] ,в связи с чем запишем:

т

фу (X,А) (А)X9, ч е [1,з . . т]. (2)

9=1

Здесь Сц(А)- коэффициенты полиномов Чебышева степени ц, функционально зависящие от амплитуды деформации А, причем для каждого класса СКД они могут быть представлены по переменной А полиномами или другими видами аппроксимирующих функций.

Выберем вид аппроксимации гистерезисной функции с учетом указанных выше особенностей ее поведения. Предварительно заметим, что в самом простейшем виде она должна содержать по крайней мере два члена, один из которых представляет прямоугольник с половинами сторон Т(А), А, а другой - эллипс с полуосями рт(А), А. Это обеспечивает сдвиг по фазе, равный п / 2, между упругой и неупругой составляющей реакции СКД, а также скачкообразное

изменение жесткости в окрестности значений деформации х = А.

Таким образом можно записать:

I 1»

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Ф =а т

Т (А) + рт (А)л/А

«2 -X2

(3)

В дальнейшем коэффициенты аппроксимации Т(А), рт(А) и сц(А) будем записывать в виде Т, рт, сц, подразумевая при этом их функциональную зависимость от амплитуды А.

В дальнейшем коэффициенты аппроксимации Т(А), рт(А) и сц(А) будем записывать в виде Т, рт, сц, подразумевая при этом их функциональную зависимость от амплитуды А.

Определим коэффициенты аппроксимации Т, р, т и вид весовой функции /(X), где X = X / А, которые обеспечивали бы минимальную среднеквадратическую погрешность Q приближения

в = I | Фт -а

Т + Р„

41—X2

/(X)ОГ. (4)

дв дв

Найдем частные производные "т^ти-и,

дТ дРт

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

-

—11_

Фт —ат + р

/ (x)dX = 0;

Ф — а\

[т + Ртт11 — X2 ) л/1 — X2 /(1^ = 0.

= 0 (5)

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

1 г / Л

1 Ф —а т Т + Рт V! — X2

—1 [ )

^ = 0.

Разрешая с этим уравнением систему интегральных уравнений (5), получим два вида весовой функции /(X):

/( х)=1 /¿Т) = [!\—¥ I"1.

Тогда минимум функционала Q (4) в пространстве переменных Т и рт определяется с помощью двух различных систем интегральных уравнений:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

для / (X ) =1,

2

—1

2

Ф -\\T + p

yfl-X

dx = 0;

Фт-o\T + p

Vi - x2 j Vi - х2 dx = 0;

Для /2( x) = 1-1

Ф т -\\T + P

J

VT

x

dx

V

= 0;

Ф m -\\T + P

Vi-x2

dx = 0.

Разрешая эти системы уравнений относительно коэффициентов аппроксимации Т и рт,

получим: , _V

Для слУчая /2(х) = (] 1 - X2 ]' -

Рт =

T=

4

8 -п2

4

8 -ж2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

J Фт cosp, Л И - - sin р jdp

0 Фт cosp, Л) sin р ■- — jdp

(6)

где р = агссоБ х.

Аналогично можно получить выражения и для коэффициентов аппроксимации Т и рт при / (X) =1: т

Рт = П |ф„(\,С0Р 4sinpi--sinp|dP? Зп2 -320 V п ]

T=зТ-зФ^Р AsinpsinP- П^Р

(7)

Таким образом, с учетом выражений (2) и(3) реакция СКД, характеризуемая УГХ в виде семейства симметричных петель гистерезиса, может быть представлена в виде приближения:

(

R = Е с xq +\

q = i

T+p

i

Л

Л2 - x2

(8)

f

х+ Е с xq +о T+p

q=i

]/л2 -x2 j=ßco(t+а), (9)

где Рсо%{ю1 + а) - безразмерная возбуждающая нагрузка;

а - угол сдвига по фазе между возбуждающей нагрузкой и перемещением;

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

Отметим, что квазигармонические колебания ВС могут быть представлены разложением Фурье в виде:

t(t) = Е cosiat + Bi sin iat,

(10)

причем на характер закона движения при t=0 накладываются следующие ограничения:

п n

Е А = Л; ЕB = 0; i е [l,3,5. . .n]. (11)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

i=i i = i Найдем решение уравнения движения (9) с помощью комбинированного асимптотического метода (КАМ), основывающегося на сочетании методов малого параметра с вариационными принципами. Согласно КАМ построение асимптотических разложений уравнения движения (9) осуществляют по степеням малого парамет-

ра s

ее1: х ^s'xa; А ^s"Л

d , пу-

й=0 й=0 тем искусственного введения в уравнение (9) дополнительных функций: для реакции СКД - Ил = g2х, а для возбуждающей нагрузки -в л = всос08 cot. Коэффициенты g2 и Ьсо находят методом Бубнова - Галеркина, причем, как показано в работе [2], g2 является коэффициентом гармонической линеаризации по Чебышеву. Учитывая изложенное, преобразуем уравнение движения (9) к виду:

&

х +Y2 х = Еsd (ßcd cosat + ßsd sin ®t ^

где сч - коэффициенты полиномов Чебышева, для условной упругой составляющей реакции, а Т и рт - коэффициенты аппроксимации для неупругой составляющей, найденные с помощью выражений (6) или (7). Следует помнить отличие этих коэффициентов на величину А, А4 для случая переменной х по сравнению с х .

Воспользуемся полученным выражением для реакции СКД (8), и запишем дифференциальное уравнение движения ВС для установившихся квазигармонических колебаний:

+s

d=0

Y x-Е Cqxq -a\T + Pm

q=i

iAT-

22 2 x 2

(12)

f

причем

ß2 = Еßcd

V d=0 J

2

+

V d=0 J

(13)

tSa = Е ßd / Е ßcd. (14)

d=0 d=0

Ограничим q=1, m=3 и заметим, что по методу КАМ для уравнения (12) ß ^ = 0. Представим уравнение движения (12) в виде системы линейных дифференциальных уравнений по степеням порядка приближения sd:

т

т

т

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

т

т

i = l

т

q

т

т

q

т

eo &&0 +у2 x0 =pc0cosat;

(15)

xj + у1 x1 = РсХ cos mt + fisi sin mt + y1 x0 -

1 -Z cqxq -Jt + PnfA

q=J

(16)

при этом

X_Xo Xj

(17)

Порождающее решение уравнения (15) имеет гармоническую форму: x0 = A0 cosmt, (18)

где

A =-rCLY- (19)

Y -m x y

С учетом решения (18) представим неупругую составляющую реакции СКД, входящую в уравнение (16), в виде гармонического разложения:

a\T + p

2 2 2 x

4Т^ 1 . = — у -$тга + А0рт БтО1. п = i

В силу гармонической формы решения(18) и выбора для аппроксимации упругой составляющей реакции СКД (в том числе и дополнительной функции Ил) полиномов Чебышева, получим [2]:

з 1

у2x0 — у = —-с3А0 соБЗа.

9=1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4

Тогда уравнение (16) приобретает следующий вид:

x +y2x1 = вС1 cosmt - 4с3А0 cos3mt +

, 4T 4T " J . . + (fl* -PnA0--)sin mt--sinimt-

п

П 1=3 i

Чтобы обеспечить в порождающем решении (18) выбор полной амплитуды первой гармоники, потребуем отсутствия в поправочном слагаемом х1 членов, содержащих косинусы, следовательно, ЬС1 = 0.

Нахождение частного решения преобразованного вышеуказанным образом линейного дифференциального уравнения (16) не представляет затруднений и закон квазигармонических колебаний ВС (17) может быть представлен в следующем виде:

x{t ) =

вс

2 2 Y - m

-cos mt- -

4T

в j - P nA0--

4(Y - 9m 2) 4T

-cos 3mt +

22 Y - m

-sin mt -

п S i'tY2 -i2m2)

(20)

sin imt.

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

вс 0 = А0(у2 —а2);

О Л 4Т-А у2 —а2

Рл = Р тА0

n 0 2 2 2

п i=j Y -1 m

С учетом выражений (11), (13), (14) окончательно можно записать

4T^ Y2-m2

AV -m2)2 +|PnA0 + I =в2-

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1 п y -1 m

tga =

- + -

n 1

4T

А А л / i ') ') ')

Y -m п40 i=1 y -1 m A = A0 - Сз A

(21)

0 4(Y2 -9m2)

Сопоставляя его с полученным выражением (19) и, учитывая ограничения на синусоидаль-

Решения данной системы уравнений является основой для исследования влияния форм петель гистерезиса СКД и параметров возбуждения на характер квазигармонических колебаний ВС. При этом полученный закон колебаний ВС (20) вместе с системой уравнений (21) позволяет определить интенсивность гармонических составляющих и особенности поведения динамических характеристик: зависимости коэффициентов передачи и резонансных частот ВС от амплитуды возбуждающей нагрузки, затягивание резонансных режимов и т.п.

В качестве примера рассмотрим применения полученных результатов (см. выражения (20) и (21)) для исследования нелинейных колебаний ВС, в состав которых могут входить многослойные пластинчатые упругодемпфирующие опоры (УДО), тросовые виброизоляторы (рис. 1 а,б) [4]. УДО (см. рис. 1,а) состоит из двух пакетов многослойных пластин (1), прижатых друг к другу пружинной скобой (2). Тросовые виброизоляторы могут быть Г-образного типа или в виде гребенки с прямолинейными участками отрезков тросов (рис. 1,б).

На рис. 2 представлены обобщенные поля петель гистерезиса СКД [4] (вышеуказанных УДО и тросовых виброизоляторов). Выделяя из представленных петель гистерезиса СКД упругие и неупругие составляющие (гистерезисные функции) (1) и вычисляя значения для коэффициентов C (А) и Р^А), Т(А) при различных Ае [0,2;6,0] с помощью выражения (6) найдем значения вышеуказанных коэффициентов.

Анализируя полученные данные замечаем, что для больших амплитуд деформация Ар>2 свойственных резонансным режимам колебаний ВС, можно отметить, что значение величины р А<<^<1. Вследствие этого для амплитуд А>2

e

2 - x 2

00

2

С 3 A0

п

+

а) б)

Рис. 1. Системы конструкционного демпфирования: а) - УДО; б) - тросовые виброизоляторы

1 — а2

А 02(1 — а2)2 +-У — п 7=1 1 — г

4 Т

■=Р2

а =

п А,

п 1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

"У Г—7

г=1 1 г

А = А 0.

4Т ^

А, = — У

1

п 1=3 г (1 — г 2а2)

•2, 2>

тема уравнений (22) отображает приближенное решение для ВС с сухим трением, причем точное решение получено в работе [5] и может быть представлено в безразмерном виде:

л=.

О2

1 2 п

—Г '8 V .

а 2а

(24)

Рис. 2. Обобщенное поле упругогистерезисных петель:

-тросовых виброизоляторов с прямолинейными

упругогистерезиными элементами; ---УДО

можно считатьрт _ 0. Вместе с тем оказывается, что при А>2 величина С3А3<<С1А (д е [1,3] см. выражение (2)), что позволяет считать С3 _ 0, С=1

Тогда система уравнений (21) для значений амплитуд А>2 приобретает следующий вид

(22)

При этом высшие гармонические составляющие А. квазигармонических колебаний вычисляются с помощью выражения (20) по формуле

(23)

Заметим, что в этом случае полученная сис-

V (1—а2)2

Сравнение точного решения (24) с приближенным (22) показывает, что при исследовании колебаний ВС с высокой степенью конструкционного демпфирования крайне важно учитывать достаточно большое количество гармонических составляющих А. (23). При этом погрешность определения амплитуды ВС по сравнению с точным зависит не только от числа учитываемых гармоник в решении (22), но и значений Ь. При

О= 4

увеличении Ь до значений близких к р = — поп

грешности могут достигать до 30 % и более в зависимости от числа учитываемых гармоник. На рис. 3 представлена зависимость коэффициента

= А

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

передачи ^ = р от безразмерной частоты а для

различного числа учитываемых гармоник (рис 3,б) при Ь= 1,25. На рис. 3,а приведены значения погрешностей § приближенного решения (22) по сравнению с точным (24) в зависимости от числа учитываемых гармоник г и уровня возбуждающей нагрузки Ь.

На рис. 4 представлена зависимость (коэффициентов передачи на резонансе ^р = (А -амплитуда резонансных колебаний ВС) от уровня возбуждения при учете г =1001.

Как видно из приведенных данных коэффициент передачи на резонансных режимах работы виброизоляторов и демпферов с конструкционным демпфированием существенно зависит от уровня возбуждения и может изменяться от значений 1,2 до 5-10 и более. Причем, начиная с в > 1.265, интенсивность изменения коэффи-

I

a\

Рис. 3. Сравнение точного решения для сухого трения с приближенным: а - погрешность приближенного решения в зависимости от числа гармоник и уровня возбуждения, во е [1,25;1.265] б - характер влияния числа гармоник, учитываемых в приближенном решении, в0 = 1.25

циента передачи на резонансе резко увеличивается. Указанная особенность приводит к высокой чувствительности виброзащитных свойств СКД на резонансных режимах их работы даже в случае незначительного изменения сил сухого трения, особенно заметного при наработке в процессе эксплуатации ВС. Поэтому на практике при конструировании СКД под заданные технические требования необходимо стремиться к обеспечению таких начальных сил сухого трения (конструкционного демпфирования), чтобы относительная возбуждающая нагрузка в была меньше значения 1,26 в процессе всего ресурса ВС.

Рис. 4. Зависимость коэффициента передачи на резонансе от уровня возбуждения

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

СПИСОК ЛИТЕРАТУРЫ

1. Лазуткин Г.В. Вынужденные периодические колебания виброзащитных систем с конструкционным гистерезисом. Вибрационная прочность и надежность двигателей и систем летательных аппаратов. Куйбышев: КуАИ, 1986, с. 81- 86.

2. ЛазуткинГ.В. Виброизоляторы из материала МР (тип ДКУ). Деп. в ВИНИТИ - №1977 - 85 от 19.03.85. - б/о №771, библ. указ. №7 - 22 с.

3. Расчет и конструирование средств виброзащиты сухого трения / В.А. Антипов, Ю.К. Пономарев, А.И. Белоусов, и др. - Самара: СамГАПС, 2005. - 207 с.

4. Чегодаев Д.Е., Пономарев Ю.К. Демпфирование. -Самара: Издательство СГАУ. 1997 - 334 с.: ил.

5. Быховский ИИ. Основы теории вибрационной техники / И.И. Быховский. - М.: Машиностроение, 1968. - 362 с.

n p

3

2

INFLUENCE OF FEATURES OF CONSTRUCTIONAL SUPPRESSION OF VIBRATION ON NONLINEAR FLUCTUATIONS OF SYSTEMS AT DISPERSION OF ENERGY

© 2009 G.V. Lazytkin, V.A. Antipov, A.L. Ryabkov

Samara State University of Means of Communicationon

The task about the elementary form of approximation of change of the characteristics of devices on dispersion of energy with a constructional hysteresis is decided(solved) on the basis of minimization of a root-mean-square error of approach(approximation), and also equality of the areas of close loops of a hysteresis to required. With the help of the combined method of decomposition of the decision of the nonlinear differential equation of movement of systems at dispersion of energy and offered form of approximation the decision for approximate repeating fluctuations taking into account influence of features of isolation of a constructional hysteresis and overwhelming devices on the basic parameters of fluctuations is found. Keywords: vibroprotective systems, hysteresis, approximation, equalization of motion, quasiharmonic vibrations.

Gennady.Lazytkin, Candidate of Technics, Doctoral Candidate. E-mail: gamaba@list.ru. Vladimir Antipov, Doctor of Technics, Professor Alexey Ryabkov, Graduate Student