Методика расчёта нагрузочных характеристик цилиндрических опор из материала МР с частичным охватом цапфы при прецессионном нагружении

Паровай Федор Васильевич; Пономарев Юрий Константинович; Котов Антон Сергеевич; Спивак Сергей Евгеньевич; Васюков Евгений Сергеевич

УДК 629.7

МЕТОДИКА РАСЧЁТА НАГРУЗОЧНЫХ ХАРАКТЕРИСТИК ЦИЛИНДРИЧЕСКИХ ОПОР ИЗ МАТЕРИАЛА МР С ЧАСТИЧНЫМ ОХВАТОМ ЦАПФЫ ПРИ ПРЕЦЕССИОННОМ НАГРУЖЕНИИ

1 Самарский государственный аэрокосмический университет 2Управляющая компания «Брянский машиностроительный завод»

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

Демпфер, материал МР, упругодемпфирующий элемент, характеристика, гистерезис, расчёт

Рассмотрим решение задачи о прецессионном нагружении цилиндрической опо-ры ротора, упругодемпфирующий элемент которой выполнен в виде незамкнутой ленты (рис. 1).

0

Рис. 1. Расчётная схема опоры с упругодемпфирую-щим- элементом- из материала МР в виде незамкнутой ленты

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

ный угол незаполненного упругодемпфи-рующим элементом пространства.

Кроме того, пусть центр прецессии смещён относительно начала координат, совпадающего с центром виброизолятора, на расстояния Х0 и Y0 вдоль осей OX и OY со -ответственно, что позволит решать задачи с предварительно нагруженной опорой посто -янно действующими силами.

Выделим двумя радиусными сечениями малый элемент с угловой протяженно -стью dj. Проекции суммарного смещения центра цапфы на оси координат найдутся в виде:

Г х = X0 + A ■ cos a;

[ y = Y0 + A ■ sin a.

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

q = q + A = A + ^JX0 + Y0 ■ cos(X — j) +

+A ■ cos(j— a)

(1)

(2)

где q| - переменная составляющая нормального смещения элемента dj.

Причём изменение угла X в зависимо -сти от положения центра прецессии описывается выражением:

*

§ =

X

arctan(—0), если Xo > 0 и Yo > 0;

Yo

p

— ,если X0 = 0 и Y0 > 0;

Xp

arctan(--) +— , если X0 < 0 и Y0 > 0;

Yo 2

p ,если X0 < 0 и Y0 = 0;

arctan( ) + —

X0 2

если X0 < 0 и Y0 < 0;

3р

— , если Х0 = 0 и 70 < 0;

X 3р

агСап(----- ) +---, если Х0 > 0 и 70 < 0.

^0 2

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

1, если — > 0;

df

dq

2, если —- < 0.

df

(4)

(5)

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

% = + А = А + у]Х02 + 702 • С08(Х -ф) +

+А • (-1)

о

где - амплитудное значение переменной составляющей нормального смещения элемента йц>.

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

(Я + Н)2 - Я2

d _____

p 6 ~ (R + H — q)2 — R2

— 1.

(6)

Аналогично амплитудное значение изменения относительной плотности 8-

p 0 •

d = (R + H )2 — R2

p 0 = (R + H — q0)2 — R2

— 1.

(7)

Сила сопротивления демпфера состоит из суммы нормальных dN и касательных сил dQ , действующих на бесконечно малом элементе dj:

dFS = dN + dQ . (8)

Для прецессионного нагружения суммарную силу сопротивления демпфера найдём в виде проекций на оси OX и OY:

[dFSX = dN ■ cosf - dQ ■ sin f;

[dFSY = dN ■ sin f + dQ ■ cos f.

Выражение (9) можно записать в виде: [dFSX = dN ■( cosf- f ■ sin f);

ldFSY = dN ■( sin f + f ■ cosf), где f - коэффициент трения, позволяющий выразить касательные силы dQ через нормальные dN.

Нормальная часть реакции dN может быть записана в виде:

(9)

(10)

т = о1 (8-,8-,у)• Цд •Я•Ь• ф. (11)

Здесь относительное обобщённое давление о*ф элемента dj определяется следующим выражением [1,2]:

1

-d-rd’V)=-d + 2• ИГ -dp)+

—k\d

d \

Р\

(12)

+(—1)v -(d)■ е ’

Po'

где 8- - текущее относительное изменение

р

плотности,

d

Po

- координата начала процесса про-

извольного процесса нагружения,

- (dp)=-(d)

P

- относительные сред-

ї д

ние напряжения цикла,

-(8 ) -*8)

- (8r ) = ■ r

P

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

относительные

ї д

обобщённые напряжения от сил трения,

а(8-) - величина «остаточной деформации» в точке начала процесса 8—,

У - параметр направления процесса нагружения, У =1 - для процесса нагрузки, У =2 - для процесса разгрузки,

Рпр - давление прессования ленты из материала МР [3].

При совместном рассмотрении (9-12) запишем выражение для проекций суммар-

ной силы сопротивления демпфера на оси ОХ и ОУ, действующей на бесконечно малом элементе dj :

сЛТ-х =аФф(8,8- ,у) • В,д ■ Я• Ь-(совф-/• Втф)ф;

-о •(13)

сЩу =афф(8,8 у) • Щд • Я Ь (тф+/ • акф)ф.

Интегрируя полученное выражение по углу ф, получим зависимость проекций силы сопротивления демпфера на оси ОХ и ОУ при прецессионном деформировании:

" 2р

рх = .К8- 8,у)' Щд ■ЯЬ'(соф-/ •

0 -(14)

2р__ 4 7

= Iстф°(8,8,У) • Щд • ЯЬ‘(втф+/ • совф)сф-

0

Из (9) получим выражения для проекций суммарной силы сопротивления демпфера на оси ОХ и ОУ, путём замены пределов интегрирования:

'

Цх = 1^(8 ,8у) • Щд ■ ЯЬ (совф-/ •sinф)dф;

в .(15)

г~

Цу = 1аф°(8,8-,у) • Щд ■ К-Ь (sinф+/ • сшф)ф.

в

На рис.2-4 представлены примеры расчётов петель гистерезиса для конструкции демпфера с частичным охватом вибратора упругодемпфирующим элементом. Основные параметры демпфера следующие: Я = 50 мм, Ь = 10 мм, Н = 30 мм, А = 5 мм, т = 250 г, в = 30°.

Очевидно, что наличие участков, не заполненных материалом МР, как показывают расчёты, проводит к искажению петель гистерезиса.

Рис. 2. Суммарный гистерезис и его составляющие в проекции на ось ОХ

с X о О

-20000 Смещение Вибратора, мм

Рис. 3. Суммарный гистерезис и его составляющие в проекции на ось ОУ

а0Х

' Ц £

Усилие сопротивления в проекции на оХ,

Рис. 4. Годографы силы сопротивления демпфера при прецессионном нагружении по круговой орбите с различными амплитудами

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

В

сп. ^

- относительный диаметр

1. d =

В

пр.

спирали; 2. - =

0

Ри

- относительная плотность уп -

-20000 Смещение Вибратора, мм

ругодемпфирующего элемента в свободном состоянии;

3. Ь = Н - относительная ширина УДЭ;

4. Я =

Я_

Н

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

тора;

5.

Т А -

Д =---относительный натяг;

Н

6. А = -

ап

- относительная амплитуда

деформирования( і = х,у);

7. / - коэффициент трения; Е • 3

к = -

- относительная жёст-

кость витка спирали, где 3 = -

ж • В,

64

g - ускорение свободного падения, равное 9,81 м/с2;

9. х =

- х

а

- относительная текущая без-

0 X

размерная деформация УДЭ вдоль оси ОХ, а0Х - «остаточная деформация» (рис.2);

10. у = - относительная текущая без-

а

07

размерная деформация УДЭ вдоль оси ОУ,

а,

07

«остаточная деформация» (рис.3);

— Г 11. Гх =

относительная сила сопро-

тивления вдоль оси ОХ,

ТХ - сила трения (см. рис.2.);

Г.,

12. = ■

£7

Т7

относительная сила сопро-

тивления вдоль оси О7,

Т7 -сила трения (см. рис.3);

Хп

13. Х0 = -

А

- относительное начальное

смещение центра траектории вдоль оси ОХ.

— 7

14 7 = Л

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

70 А

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

смещение центра траектории вдоль оси ОУ.

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

В качестве примера на рис.5 и рис.6 представлены безразмерные петли в проекциях на оси ОХ и ОУ, вычисленные для следующих критериев подобия: d =10, — 0 =

0.36, Ь = 0.33, Я = 1.66, А = 0.30, А = 1.8

/ = 0.15, к = 24001 30°.

X0 = 0, 70 = 0, в =

Рис.5. Зависимость относительной силы сопротивления демпфера в проекции на ось ОХ от амплитуды смещения вибратора в безразмерном виде

Рис. 6. Зависимость относительной силы сопротивления демпфера в проекции на ось О7 от амплитуды смещения вибратора в безразмерном виде

Как показывают расчёты, вырез в уп-ругодемпфирующем элементе опоры приводит к разным искажениям полей упругоги-стерезисных петель в проекциях на координатные оси ОХ и 07.

Выводы

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

4

препятствующей проворачиванию демпфера в опоре) имеет некоторое не заполненное материалом МР пространство, форма гистерезиса изменяется с одновременным сдвигом центра петли по оси силы в сторону не заполненного пространства. Опора при этом приобретает свойство саморазгрузки от веса ротора. Если статической нагрузки на опору нет, центр цапфы смещается в сторону незаполненного пространства; если такая нагрузка есть, то под её действием центр цапфы может возвратиться в геометрический центр опоры. Предложенная методика позволяет расчётным путём подобрать такое сочетание параметров (плотность материала МР, натяг в опоре, угол охвата коаксиального пространства, занятый материалом МР), при котором центр цапфы будет совпадать с центром опоры. Это является важным для обеспечения надёжности и долговечности роторных машин по уровню вибраций и износу уплотнений.

Библиографический список

1. Котов, А.С. Исследование количественных характеристик технологии изготовления изделий из материала МР/ А. С. Котов, Ю.К. Пономарев, Ф.В. Паровай // Вестник СГАУ «Проблемы и перспективы развития двигателестроения»: сб. науч. тр. - Самара: СГАУ, 2006.- ч.2.- С. 397-401.

2. Уланов, А.М. Основы проектирования систем виброзащиты с упругими элемента-

ми из материала МР / А.М. Уланов, Ю.К. Пономарев // Известия Самар. Науч. центра РАН, том 10, № 3, июль - сентябрь 2008 г. -с. 853 - 857.

3. Пономарев, Ю.К. Разработка математической модели виброизолятора с шаровым упругодемпфирующим элементом из материала МР / Ю.К. Пономарёв [и др.] // Вестник СГ АУ «Проблемы и перспективы развития двигателестроения»: сб. науч. тр. -

Самара: СГАУ, 2006.- ч. 1.- С. 392-397.

References

1. Kotov A.S. Investigation of quantitative characteristics of the technology of production of MR material/ A.S. Kotov, Y.K. Ponomaryov, F.V. Parovay // SSAU Bulletin. Problems and Perspective of Engine Designing Series. Samara, 2006, Pt.2. 397-401.

2. Ulanov, A.M. Fundamentals of vibration protection systems with elastic elements of the MR material / A.M. Ulanov, Y.K. Ponomaryov // Samara Science Center of RAS Bulletin, volume 10, № 3, July - September 2008. - p. 853

- 857.

3. Ponomaryov, Y.K. Development of mathematical model of vibration-isolator with spherical damping element of the MR material / Y.K. Ponomaryov [and others] // SSAU Bulletin. Problems and Perspective of Engine Designing Series: - Samara: SSAU, 2006.- Pt. 1.-p. 392-397.

METHODS OF CALCULATING THE LOAD CHARACTERISTICS OF RING DAMPERS OF THE MR MATERIAL WITH PARTIAL COVERAGE OF A BEARING WITH PRECESSIONAL LOADING

1Samara state aerospace university 2Bryansk Engineering Plant

The paper shows the method of calculating elastic and damping characteristics of ring damper made of the MR material with partial coverage of the coaxial space between a support and a cylindrical vibrating bearing. Shown that a precessional loading results in elliptical hysteresis loops at projections on coordinate axes of damper. The magnitude of reactions depends on the geometrical parameters of support, the density of the wired material MR, prior tightness of damper. Found that the volume of a gap in the ring damper deforms the elliptical shape of hysteresis loops and their orientation in the coordinate system of load characteristics.

Damper, hysteresis, material MR, calculating, damping characteristics

Информация об авторах

Паровай Федор Васильевич, кандидат технических наук, доцент, главный инженер отраслевой лаборатории вибрационной прочности и надёжности авиационных изделий Самарского государственного аэрокосмического университета. E-mail: [email protected]. Область научных интересов: теория и практика виброзащиты, контактные уплотнения из материала МР, автоматизация производства прессованных проволочных материалов.

Пономарев Юрий Константинович, доктор технических наук, профессор, заместитель главного конструктора отраслевой лаборатории вибрационной прочности и надёжности авиационных изделий Самарского государственного аэрокосмического университета. E-mail: [email protected]. Область научных интересов: теория и практика виброзащиты, математическое моделирование гистерезисных систем, конструкционное демпфирование.

Котов Антон Сергеевич, кандидат технических наук, ведущий специалист отраслевой лаборатории вибрационной прочности и надёжности авиационных изделий Самарского государственного аэрокосмического университета. E-mail: [email protected]. Область научных интересов: теория и практика виброзащиты, математическое моделирование гистере-зисных систем, конструкционное демпфирование.

Спивак Сергей Евгеньевич, аспирант Самарского государственного аэрокосмического университета. E-mail: [email protected]. Область научных интересов: теория гистере-зисных систем, автоматизация исследования гистерезиса, конструирование виброзащитных систем на базе материала МР, тросов, лент.

Васюков Евгений Сергеевич, технический директор управляющей компании «Брянский машиностроительный завод». E-mail: [email protected]. Область научных интересов -проблема шума и вибрации в транспортном машиностроении.

Parovay Fedor Vasilievich, candidat of technical Science, associate professor of Samara State Aerospace University. E-mail: [email protected]. Area of research: theory and practice of vi-broprotection, contact seals made of MR material, automation of production process of pressed wired materials

Ponomaryov Yuri Konstantinovich, Doctor of Technical Science, professor of Samara State Aerospace University. E-mail: [email protected]. Area of research: theory and practice of vibroprotection, mathematical modeling of systems with damping, structural damping.

Kotov Anton Sergeevich, Ph.D. in Technical Science, Leading Specialist of Research Laboratory of vibration strength and reliability of aviation products, Samara State Aerospace University. E-mail: [email protected]. Area of research: theory and practice of vibroprotection, mathematical modeling of systems with damping, structural damping.

Spivak Sergey Evgenievich, postgraduate of Samara State Aerospace University. E-mail: [email protected]. Area of research: theory of systems with damping, automation of hysteresis research, design of vibration protection systems made of MR material, wire rops and metallic ribbons.

Vasyukov Evgeniy Sergeevich, CTO of Bryansk Engineering Plant. E-mail: [email protected]. Area of research: problem of noise and vibration in transport engineering.