Научная статья на тему 'Устойчивость движения неравножесткого ротора в упругих опрорах'

Устойчивость движения неравножесткого ротора в упругих опрорах Текст научной статьи по специальности «Физика»

CC BY
133
49
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УСТОЙЧИВОСТЬ ДВИЖЕНИЯ / РОТОР / УПРУГИЕ ОПОРЫ / АНИЗОТРОПИЯ УПРУГОГО ПОЛЯ / ROTOR / ANISOTROPY OF ELASTIC FIELD / A STABILITY OF A MOTION / ELASTIC SUPPORTS

Аннотация научной статьи по физике, автор научной работы — Клочков Б. Ф.

Исследована устойчивость движения неравножесткого ротора при трех вариантах конструкции его опор. Определена оптимальная схема ротора с невращающимися упругими опорами.

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

The article gives research of a stability of a motion of a various stiffness rotor in three versions of the construction of its supports. An optimal scheme of a rotor in nonrotating elastic supports is defined.

Текст научной работы на тему «Устойчивость движения неравножесткого ротора в упругих опрорах»

при годовой ставке 5 % — на 0.81 %, при годовой ставке 6 % — на 1.19 %, при годовой ставке 7 % — на 1.66 %, при годовой ставке 8 % — на 2.21 %, при годовой ставке 9 % — на 2.86 %.

В долларах сокращения размеров выплат составляют;

145505.92,

214487.95,

298569.12,

398447.29,

514774.80,

что свидетельствует о целесообразности реализации «мягкого» режима ежемесячных платежей при значительных размерах ставок и продолжительном времени выплат по инвестициям (рис. 1).

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

Список литературы

1. Крылов А. Н. О численном решении уравнения, которым в технических вопросах определяются частоты малых колебаний материальных систем / А. Н. Крылов // Изв. АН СССР. — 1931.

2. Уоткинс Д. С. Основы матричных вычислений / Д. С. Уоткинс. — М.; БИНОМ; Лаб. знаний, 2006. — 664 с.

3. Сахаров В. В. Применение матрицы Крылова для апериодического управления динамическими объектами / В. В. Сахаров, В. И. Королев // Журнал Университета водных коммуникаций. — 2011. — Вып. 1 (9).

4. Сахаров В. В. Совершенствование управления качеством сборки судовых механизмов в судоремонте / В. В. Сахаров, А. А. Кузьмин. — СПб.; Судостроение, 2012. — 202 с.

УДК 531 Б. Ф. Клочков,

канд. техн. наук, доцент, ГУМРФ имени адмирала С. О. Макарова

УСТОЙЧИВОСТЬ ДВИЖЕНИЯ НЕРАВНОЖЕСТКОГО РОТОРА

В УПРУГИХ ОПРОРАХ

STABILITY OF VARIOUS STIFFNESS ROTOR MOTION IN ELASTIC SUPPORTS

Исследована устойчивость движения неравножесткого ротора при трех вариантах конструкции его опор. Определена оптимальная схема ротора с невращающимися упругими опорами.

Выпуск 3

|Выпуск 3

The article gives research of a stability of a motion of a various stiffness rotor in three versions of the construction of its supports. An optimal scheme of a rotor in nonrotating elastic supports is defined.

Ключевые слова: устойчивость движения, ротор, упругие опоры, анизотропия упругого поля. Key words: a stability of a motion, a rotor, elastic supports, anisotropy of elastic field.

П

РИ вращении роторов в ряде случаев возникают колебания с частотами, кратными скорости их вращения. Известно, что одной из причин, вызывающих колебания с частотой, равной удвоенной скорости вращения, является неоднородность радиального коэффициента жесткости ротора (анизотропия упругого поля вала) [1; 2]. Примерами анизотропных роторов, то есть роторов двоякой жесткости, являются роторы двухполюсных электромашин; валы, сечения которых на некоторых участках в силу назначения машины прямоугольные.

Известные практические рекомендации — отстройка рабочих диапазонов частот от виброопасных, установка упругих опор для уменьшения анизотропии суммарного упругого поля между ротором и фундаментом [1; 3, с. 27-44] — оказываются часто недостаточными.

Возможен другой способ активной борьбы с указанным выше нежелательным явлением, заключающийся в применении дополнительного анизотропного упругого поля, вращающегося вместе с валом. Последнее может быть создано как на некоторых участках вала, так и с помощью специальных упругих опор [4]. Это обстоятельство открывает новый путь уменьшения ширины зон неустойчивости движения и амплитуд колебаний кратных частот.

Идею этого метода рассмотрим на примере простейшей роторной системы, состоящей из одного анизотропного вала с установленным в середине его и без перекоса диском (рис. 1).

Вводится система координат О^ц, вращающаяся вместе с ротором с постоянной угловой скоростью ш. В этой системе вал имеет жесткости c1 и с2 вдоль осей £, и ц и вращается в подшипниках, внутренние кольца которых также анизотропны с коэффициентами жесткости, равными соответственно c01 и с02 вдоль тех же осей. Диск с приведенной к его центру массой вала имеет массу m и удельный дисбаланс е с составляющими е1 и е2 вдоль осей £, и ц соответственно. Координаты геометрического центра диска и центров опор обозначим соответственно через £,, ц и 2Ю, ц0. При этом в силу симметрии считаем, что перемещения обеих опор одинаковы. Сила сопротивления движению принимается пропорциональной абсолютной скорости центра диска (коэффициент пропорциональности к).

Без учета сил неупругого сопротивления и массы опор уравнения движения имеют вид: = -к{4 - юл) - Ci - 4о) + 2/и сор + ми2(^ + ej + mg cos at,

mr\ = -A(f| + co4)-c2(p-Po)_ 2mat, + та2 (p + e2) - mg sin at,

0 = -2с0£0+С1£-Ы,

или после исключения £,0 и р0:

О = -2с02р0 + с2 (р - р0)

mt, + k(lI - юр) - 2/исор + (с® - та2 )t, = теуя2 + mg cos at,

(1)

/яр + А:(р + ш^) + 2та£, + (с(2) - та2 )р = те2а2 - mg sin at,

где С1 и С2) — эквивалентные коэффициенты жесткости, определяемые формулами:

11 1 11 1

с0) с + 2с ’ с(2) _ с + 2с '

°Э С1 ^С01 Э с2 ^С02

(2)

Если подобрать жесткости вала и опор так, что С1 = c(2), то полученные выше уравнения совпадают с уравнениями движения изотропного вала, поэтому зона неустойчивости между критическими скоростями стягивается в точку вследствие их равенства [1; 2]. Колебания с частотой 2ш, вызываемые весом, также не возникают.

Примем для сил неупругого сопротивления при деформации вала линейную зависимость от скорости деформации [5, с. 45-75]. Наряду с традиционно выполненными упругими опорами с невращающимся изотропным (однородным) упругим полем выше обоснована возможность применения опор с вращающимся анизотропным полем. Для сопоставления результатов параллельно проведем исследование устойчивости движения ротора двоякой жесткости в жестких опорах. С этой целью исследуем устойчивость движения ротора двоякой жесткости при трех вариантах установки (динамических схемах 1-3) (рис. 1).

Уравнения движения

Воспользуемся расчетной схемой, изображенной на рис. 1, и обозначениями, приведенными выше. При движении ротора возникает сила сопротивления движению (сила внешнего демпфирования), принятая пропорциональной абсолютной скорости центра О диска с коэффициентом пропорциональности к. При деформациях вала и упругих элементов опор в схемах 2 и 3 возникают силы неупругого сопротивления, принятые, согласно гипотезе Фохта, пропорциональными соответствующим скоростям деформаций с коэффициентом h.

Рассмотрим вначале схему 1. Для ротора, установленного по этой схеме, сила неупругого сопротивления F записывается непосредственно в проекциях на оси £, и ц, вращающиеся вместе с ротором: F' = -h£, F^ = -hip Так как сила внешнего демпфирования R пропорциональна абсолютной скорости центра диска, то она должна быть записана вначале в неподвижной системе координат Rx = -kx, Ry = -ку. По формулам преобразования координат:

x = £ cosat - ^ sin®t, у = £ cosat + ^ cosat (3)

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

7?! =-ki, + &юг|, Яц = -Ат) - kat,.

Таким образом в принятых обозначениях уравнения движения ротора по схеме 1 имеют вид [1]:

Выпуск 3

|Выпуск 3

т^ + (к + /г)4 -2тт\- tor) + (q - та2 )£ = mqco2 + mg cosat,

тц + (к + h)\\ +2та^ + ка£, + (c2 -т(й2)т\ = те2а2 - mg sin at (4)

Для обобщения результатов исследования полезно перейти к безразмерным (относительным) параметрам. Обозначим через ю0 = лАа + сД / 2т усредненную критическую скорость ротора в жестких опорах, р = а/а0 — безразмерную скорость вращения, 5 = к/та0 — безразмерный коэффициент внешнего демпфирования, к = Н/та0 — относительный коэффициент неупругого сопротивления, ц = (q - q)/(q + c2) — коэффициент анизотропии упругого поля вала и т = —

безразмерное время; тогда, обозначив производные d/dx штрихами, приведем уравнения движения (4) схемы 1 к виду:

if + (5 + - 2Р^' + (1 - р -р2)г| + Щ = е2р2 -^sinpx,

ю0

§' + (84- ky-' - 2pi!' + (И-|и-Р2)^-Р5л = qP2 + 4C0SPT- (5)

ш0

Обратимся теперь к схемам 2 и 3. Для этих схем силы неупругого сопротивления в вале пропорциональны относительным скоростям (£, - £0) и (т| - т|0). Силы неупругого сопротивления при деформациях упругих элементов опор для схемы 2 записываются сначала в неподвижной системе координат х, у, в которой они имеют вид -2hx0 и -2hy0, а затем, после преобразования по формулам (3), во вращающейся системе £,, ц: -2h(£0 - щц0) и -2h(p0 + ш2ю).

Для перехода к безразмерным параметрам обозначим через а = (q + q)/2q относительную податливость опор. Тогда уравнения движения ротора, установленного по схеме 2, примут вид:

Г + *(§' - 4о) + 8«' - Рп) + (1 + р)Й - §0) " 2РЛ' - РЧ = qP2 + 4C0S Р*>

ю0

л" + *(и' -1По) + 8(п' - Р5) + (! - Р)(П - По) + 2Р4' - Р2П = q>P2 - 4sinPT>

®о

К К S'- 2РПо) + 2^0 / а - (1 + - §о) = 0,

^(ЗПо - П' + 2Р^0) + 2п0 / а - (1 - Р)(П - Но) = °- (6)

Наконец, в схеме 3 упругое поле вращается вместе с ротором, поэтому силы неупругого сопротивления в опорах описываются непосредственно во вращающейся системе координат слагаемыми -2h£,0 и -2hp0. После приведения к безразмерному виду и некоторых преобразований уравнения движения ротора, установленного по схеме 3, принимают вид:

V + *£' - О + 8($' - Рп) + (1 + т, - $о) - 2Рн' - РЧ = qP2 + 4C0SPX>

ю0

п" + *(и' - По ) + 8(п' - Р£) + (1 - Р)(П - По ) + 2Р4' - Р2п = qp2 -4 sin Рт>

со0

*(3& - - 2РП0) + 2^о / а - (1 + р)($ - $о) = 0,

^(ЗПо - П' + 2р^0) + 2п0 / а - (1 - Р)(и - По) = °>

(7)

где а = (q + q) / (с + c ) — относительная усредненная податливость опор; v = (с - c ) /

(с + c ) — коэффициент анизотропии упругого поля опор. Отметим, что без учета сил неупругого сопротивления (к = 0) уравнения (7) принимают вид:

Г + 8(Г-РП)-2РП' +

2 (1 + Ц.)(1 + у) п2

а(1 + р) + 2(1 + v)

. (1-P)(1-V) п2

а(1 - р) + 2(1 - v)

4 = qP2 + 4C0SPX’

ап

n = e2P2+4sinPT-

а0

П' + 8(т1' + р§)-2рт,' +

(8)

Условие (2) в безразмерных параметрах записывается в виде

(1 + nXl + v) _ (1-nXl-v) а(1 + р) + 2(1 + v) а(1 - ц) + 2(1 - v)

Из квадратного относительно v уравнения (9) определяется значение

v° = [а(1 - р2) - ^а2(1-р2)2 +16р2] / 4р.

(9)

(10)

Анализ знаков показывает, что в формуле (10) перед корнем следует брать лишь «-». График зависимости v° = v°(a), построенный для примера при значении р = 0,2, изображен на рис. 2.

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

л' + 8(л' + р§) + 2К' + (-

----у - Р2)Л = е2Р2 + -%sin|k,

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

/ЖЙП Юп

V + «(§' - Рл) - 2Рл' + (^у - р2)^ = е,р2 + 4-cospT,

тап

оап

где с, = та:

2 2(l + p)(l + v) а(1 + р) + 2(1 + v)

■ та:

2 2(1-рХ1-у) а(1 - р) + 2(1 - v)

Эти уравнения совпадают с уравнениями движения изотропного вала, поэтому зона неустойчивости между критическими скоростями стягивается в точку вследствие их равенства. Колебания с частотой 2р, вызываемые весом, также не возникают. Динамика ротора, установленного по схеме 3, с учетом сил неупругого сопротивления исследуется ниже.

Устойчивость движения ротора в жестких опорах

Укажем предварительно, что в такой же постановке устойчивость этой схемы была исследована в [1], где получены следующие результаты.

1. Движение устойчиво во всем диапазоне скоростей, за исключением интервала между

критическими скоростями coj = ■Jcl / т и а2 = ^с2 / т или, в принятых здесь безразмерных параметрах, между Р1=-\Д-Й и Рг= \Д+Й.

2. С увеличением коэффициента внешнего демпфирования 5 зона неустойчивости сужается и сдвигается в сторону более низких скоростей. Однако вывод о неограниченной устойчивости анизотропного ротора за второй критической скоростью ш2 находится в противоречии с известным результатом, полученным в [1], о неустойчивости движения изотропного ротора в жестких опорах за критической скоростью. Поэтому проведем подробное исследование устойчивости дви-

Выпуск 3

|Выпуск 3

жения ротора, установленного по схеме 1. Применяя метод П-разбиений, решение однородной (с нулевыми правыми частями) системы уравнений (5) отыскиваем в виде

4 = A exp(ip0x), р = B exp(ip0 т), (11)

где А, В — произвольные постоянные; р0 — частота колебаний на границе устойчивости во вращающейся системе координат; i = V-I.

Подстановка решений (11) в однородные уравнения (5) приводит к характеристическому определителю

д= Р2-Р2-Ро+г‘(8 + £)А) — Р5 — 2/Rpo =0

р8 + 2фр0 Pi - Р2 - pi + г'(8 + к)р0

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

р40 -[2(1 + Р2) + (к + 8)2]/>о + (Pf -Р2)(Р| -Р2) + 82Р2 = 0, (12)

\-Ро -1 + Р2 - 28р2 I {к + 8)]/70 = 0.

(13)

Из уравнения (13) следует, что на границе зоны устойчивости частота колебаний, по форме которых происходит потеря устойчивости, равна 0. Подставив это значение р0 в (12), получим уравнения двух границ зоны неустойчивости

in,2

= 0,5

2-82 ± д/(2 — 82)2 -4(1-ц2)

(14)

Уравнения (14) совпадают с уравнениями, полученными в [1]. Однако уравнения (14) показывают, что на устойчивость движения схемы 1 никак не влияет внутреннее трение. Это свидетельствует о неполноте исследования. Возвратимся поэтому к уравнению (13), из которого найдем зависимость другой, ненулевой частоты собственных колебаний на границе устойчивости от скорости вращения

Pq =! — р2 — 2р28/ {к + 8).

Подставив это выражение в (12), получим уравнение третьей границы области устойчивости:

ргз = ^{гЪ + р -Ас)ТЩ, (15)

где a = (2k/5)2, b = k4/ 52 + 2k3/5 - 2k5 - 4(k/ 5 + 1)2, с = -(k/5 + 1)2 [ц + (k + 5)2].

Для получения действительных значении РГ3 в уравнении (15) перед корнем в подкоренном выражении для РГ3 следует брать знак «+».

Частота неустойчивых колебаний на третьей границе устойчивости определяется из уравнения (12) после подстановки в него Р = РГ3:

Роъ = ,

где d = -1 - Р23 - (k + S)2/ 2, f = (Р2 - Р23) (Р2 - Р23) + 52р2

ГЗ'

На рис. 3 сплошными жирными линиями показаны границы областей устойчивости движения анизотропного ротора в жестких опорах (штриховкой обозначены зоны неустойчивости). Как видно из рис. 3, вторая и третья границы пересекаются, откуда следует важный вывод: при малом внешнем демпфировании (в данном случае 5 = 0,02) движение неустойчиво во всем диапазоне за первой границей; причиной этой неустойчивости является внутреннее трение. При увеличении внешнего демпфирования движение становится вновь устойчиво за второй границей, но не во всем диапазоне скоростей, как это показано в [1], а лишь до третьей границы, после чего система окончательно теряет устойчивость.

Рис. 3. Границы областей устойчивости движения схемы 1 при р = 0,2 и к = 0,16

Итак, ротор, установленный в жесткие опоры, имеет предел устойчивости, отодвигающийся в сторону больших скоростей при увеличении демпфирования. Частота р2 развивающихся колебаний на первой и второй границах области устойчивости равна нулю во вращающейся системе координат, то есть равна частоте вращения в неподвижной системе, следовательно, на первой и второй границах система теряет устойчивость по форме прямой синхронной прецессии. График модуля частоты \р | (во вращающейся системе координат) изображен на рис. 3 штриховой линией. Этим частотам соответствуют в неподвижной системе координат частотыр31 = РГ3 -роъ и р32 = РГ3 + роу графики которых изображены сплошными линиями. Отметим, что при любом демпфировании рЪ1 = 1, то есть частота равна усредненной критической скорости рЪ1 = (Р1 + Р2) / 2.

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

Как следует из уравнений (14), в случае жестких опор ротора не существует практического способа подавления зоны неустойчивости [Рп и РГ2], так как требуемое для этой цели искусственное введение внешнего демпфирования влечет за собой применение упругодемпферных опор, то есть переход к динамическим схемам 2 или 3 анизотропного ротора в упругих опорах.

Устойчивость движения ротора в невращающихся упругих опорах

Устойчивость движения этой модели без учета внутреннего и внешнего трения была исследована в [1], где обнаружена единственная зона неустойчивости между критическими скоростями ш1 и ш2 или в безразмерных параметрах между

I 2(1-р) и = I 2(1+ р)

\ 2 + а(1 - р) ^ ^ 2 + а(1 + р)

(16)

Результаты численного исследования [6] представлены на рис. 4 (внутри изображенного объема движение устойчиво). Сечение граничной поверхности устойчивости-неустойчивости плоскостью 5 = 0,1 показано на рис. 5, на котором также нанесен график модуля частоты \ _р03| на третьей границе. Частотам ±р03 соответствуют в неподвижной системе две частоты: ру1 = Ргз - роз иру2 = Ргз + роз. График частоты р31 показан сплошной линией, график же частоты р3 2 ввиду больших ее значений не изображен.

Рис. 4. Граничная поверхность областей устойчивости движения схемы 2 при р = 0,16 и р = 0,2

Выпуск 3

Выпуск 3

Рис. 5. Границы областей устойчивости движения схемы 2 при 5 = 0,1

Как и в случае жестких опор, на первой и второй границах устойчивости во вращающейся системе координат р01 = р02 = 0 или в неподвижной системе р1 = Рп и р2 = РГ2, то есть потеря устойчивости на этих границах происходит по форме прямой синхронной прецессии. Частота р31 также является усредненной собственной частотой системы: р31 ~ 0,5(Р1 + Р2), где Р1 и Р2 — критические скорости без учета демпфирования, определяемые выражениями (16). Кривые Р1(а) и Р2(а) показаны на рис. 5 штрихпунктирными линиями. Как следует из рис. 4 и 5, при естественном внешнем демпфировании, под которым следует понимать демпфирование с коэффициентом 5 < 0,2, с ростом податливости опор а зона неустойчивости между критическими скоростями подавляется уже при жесткости каждой опоры, примерно равной жесткости вала. Ширина интервала устойчивой работы достигает максимума при а ~ 3. С увеличением внешнего демпфирования 5 зона устойчивости примерно линейно расширяется. Обнаружение третьей границы дополняет результаты исследования устойчивости, полученные в [1].

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

На рис. 6 изображены граничные поверхности устойчивости-неустойчивости в пространстве параметров а, Р и v, построенные по результатам численного исследования. Влияние коэффициента анизотропии v упругого поля опор на устойчивость движения показано на рис. 7, на котором изображено сечение граничных поверхностей плоскостью а = 1. Влияние коэффициента жесткости опор на устойчивость движения схемы 3 иллюстрирует рис. 8, на котором показаны сечения граничных поверхностей плоскостью v = 0,07. Наконец, влияние внешнего демпфирования 5 на ширину зон устойчивости показано на рис. 9.

Кроме того, на рис. 7-9 штрихпунктирными линиями показаны критические скорости Р1 и Р2, определяемые для схемы 3 без учета демпфирования формулами:

Рис. 6. Граничная поверхность области устойчивости движения схемы 3 при р = 0,2, к = 0,16 и 5 = 0,1

= 12 2(1 + ц)(1 + у) и = I 2(1-ц)(1-у)

1 V а(1 + (х) + 2(1 + v) Н2 V а(1 - ц) + 2(1 - v)

(17)

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

Рис. 8. Границы областей устойчивости движения схемы 3 при v = -0,07 и 5 = 0,1

Рис. 9. Границы областей устойчивости движения схемы 3 при а = 1 и v = - 0,9

Анализ приведенных результатов показывает, что выбор значения коэффициента анизотропии v вращающейся упругой опоры вблизи значения v°, определяемого выражением (10) (см. рис. 2), гарантирует подавление зоны неустойчивости между критическими скоростями без специального введения внешнего демпфирования. При этом даже значительное несовпадение этих значений, то есть «расстройка» по коэффициенту v не приводит к появлению неустойчивых колебаний в зоне [Р1, Р2]. Внешнее демпфирование усиливает этот эффект, и, кроме того, с ростом 5 зона устойчивой работы расширяется примерно по линейному закону. Частоты собственных колебаний на всех трех границах подчинены тем же закономерностям, что и в двух предыдущих схемах.

Выпуск 3

|Выпуск 3

Выводы

Прежде всего напомним, что количественные оценки верны для взятого в качестве примера случая р = (с1 - с2 ) / (с1 + с2) = 0,2 и h = 0,16 тю0, то есть к = 0,16.

Проведенный анализ устойчивости движения показал, что для всех трех схем установки вала двоякой жесткости применим общий способ подавления зоны неустойчивости между критическими скоростями путем увеличения внешнего демпфирования. В схеме 1, где этот способ единственный, указанная зона исчезает при 5 > 0,2. В схеме 2 увеличение податливости опор позволяет подавить неустойчивость при меньших значениях внешнего демпфирования, то есть для этой схемы зона неустойчивости между Р1 и Р2 может быть подавлена двумя способами — увеличением либо 5, либо а.

Что касается схемы 3, то она в указанном смысле «безразлична» к жесткости опор и внешнему демпфированию, так как при любом а и любом малом (но не нулевом, а реальном) 5 неустойчивость может быть подавлена путем выбора коэффициента анизотропии v упругого поля опоры, близким по значению к v°.

После подавления зоны неустойчивости между критическими скоростями возникает вопрос о ширине области устойчивой работы до окончательной потери устойчивости. В схеме 1 при 5 = 0,2 потеря устойчивости происходит при значении РГ3 = 2,25. С ростом внешнего демпфирования РГ3 также растет, но, по-видимому, значения 5 > 0,2 соответствуют специальным мерам по увеличению демпфирования (см. рис. 3). В схеме 2 с увеличением податливости опор РГ3 сначала увеличивается, достигая максимума при а ~ 3, затем уменьшается. Ориентировочно зона устойчивости в 3-4 раза превышает среднюю критическую скорость. Увеличение внешнего демпфирования примерно линейно расширяет зону устойчивости, причем эта тенденция имеет место при любых значениях податливости опор (см. рис. 5). Наконец, в схеме 3 увеличение податливости опор резко сужает область устойчивости, доводя ее ширину до полутора критических скоростей (5 = 0,1) (см. рис. 8). Увеличение внешнего демпфирования также почти линейно расширяет зону устойчивости при всех а. Максимальная ширина области устойчивого движения достигается при малых а в случае, если опоры близки к изотропным (v ~ 0); при большой податливости опор значение v в этом смысле безразлично.

Во всех трех схемах потеря устойчивости на первых двух границах происходит по форме прямой синхронной прецессии, окончательная же потеря устойчивости на третьей границе — по форме прямой несинхронной прецессии с частотами р31 = 0,5(Р1 + Р2) и рЪ2 = 2РГ3 - 0,5(Р1 + Р2 ), где РГ3 — скорость вращения на третьей границе.

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

Список литературы

1. Диментберг Ф. М. Изгибные колебания вращающихся валов / Ф. М. Диментберг. — М.: АН СССР, 1959. — 247 с.

2. Вибрации в технике: справ.: в 6 т. — М.: Машиностроение, 1980. — Т. 3: Колебания машин, конструкций и их элементов. — 544 с.

3. Цырлин А. Л. Динамика роторов двоякой жесткости / А. Л. Цырлин // Динамика гибких роторов. — М.: Наука, 1972.

4. Подшипник качения: А. с. 314009 СССР / Кельзон А. С., Лущик В. М., Циманский Ю. П. — Б. и., 1971. — № 27.

5. Келъзон А. С. Динамика ротора двоякой жесткости в упругих опорах / А. С. Кельзон, Б. Ф. Клочков // Динамика систем. — Горький, 1976. — Вып. 10.

6. Клочков Б. Ф. Вычислительные значения параметра на границе устойчивости решений системы обыкновенных линейных дифференциальных уравнений (алгоритм № П001219) / Б. Ф. Клочков // Алгоритмы и программы: информ. бюл. — М.: ВИНИТИ, 1975. — № 2.

УДК 620.197:629.12(4) В. Б. Леонов,

д-р техн. наук, профессор, Херсонская государственная морская академия;

И. И. Рублев,

аспирант,

Херсонская государственная морская академия

ИССЛЕДОВАНИЕ КОРРОЗИОННОЙ УСТОЙЧИВОСТИ МЕТАЛЛИЧЕСКИХ КОРПУСОВ СУДОВ,

ОБРАБОТАННЫХ ЭФФЕКТИВНЫМ НАПОЛНИТЕЛЕМ

STUDY OF CORROSION RESISTANCE OF METALLIC HULLS OF SHIPS TREATED BY EFFECTIVE FILLERS

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

In the article the analysis of existing methods for protecting metal hulls is carried out, a new technical solution for the effective filler having greater resistance to corrosive degradation in comparison with the conventional protective coatings is made. A new environmentally safe for marine, environmental filler is developed and tested. The practical uses of this coating, which as compared with the current standard coating is elastic are shown. This coverage does not have any existing analogues and has a relatively low cost because it is composed of industry and agriculture waste products. Covering having a high adhesion ability can keep on the metal surface with various parameters of the environment, as evidenced by the above experiments described in this article.

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

Key words: filler protective coating, filler component, crude rubber, aggressive environment, corrosion, corps, ship, steel, marine environment.

Введение

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

Выпуск 3

i Надоели баннеры? Вы всегда можете отключить рекламу.