Проблемы машиностроения и надежности машин, 2019, № 5, стр. 10-17

ВЛИЯНИЕ ТИПА ПОДШИПНИКОВ НА ВИБРОАКТИВНОСТЬ ДВУХПРОЛЕТНОГО РОТОРА С КОСОСИММЕТРИЧНЫМ ОСТАТОЧНЫМ ДИСБАЛАНСОМ

О. А. Волоховская 1*, О. В. Бармина 1

1 Институт машиноведения им. А.А. Благонравова
г. Москва, Россия

* E-mail: olgaavol@yandex.ru

Поступила в редакцию 24.04.2019
Принята к публикации 25.06.2019

Полный текст (PDF)

Аннотация

Исследована зависимость виброактивности высокотемпературного отдела турбоагрегата, включающего в себя ротор высокого давления и ротор средне-низкого давления, при кососимметричном остаточном дисбалансе каждого из роторов от типа наиболее часто используемых подшипников. Рассмотрены шестиколодочные сегментные и эллиптические подшипники. При исследовании учтено наличие в системе неустранимого прогиба валопровода. Расчеты выполнены для турбины К-300-23.5 производства ЛМЗ. Приведен ключевой результат анализа.

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

Введение. Виброактивность турбоагрегата (ТА) определяется прежде всего величиной и распределением остаточного дисбаланса после балансировки роторов и дисбаланса, вызванного неустранимым прогибом, порожденным температурной нестабильностью роторов или несоблюдением условий пуска [13]. Кроме указанных факторов существенное влияние на уровень виброактивности оказывает тип подшипников, на которых установлен турбоагрегат [4]. Известно, что амплитуды колебаний в точках установки ТА на подшипники при рабочей скорости должны удовлетворять нормам ISO (International Standart Organisation): 1) $A \leqslant 90$ мкм; 2) в серединах пролетов не превышать величины зазоров в уплотнениях. Тип подшипников влияет на выполнение обоих условий ограничения амплитуд. Однако это влияние становится определяющим для выполнения первого условия для систем ротор высокого давления – ротор средне-низкого давления турбоагрегата (“РВД–РСНД” ТА) с кососимметричным остаточным дисбалансом каждого из роторов. При переходных колебаниях на выбеге турбоагрегата, возникающих на собственных частотах, которым соответствуют формы колебаний с близкими к середине пролетов узлами, наличие в системе кососимметричного дисбаланса способствует возрастанию амплитуд колебаний в точках установки системы на подшипники. При неблагоприятном взаимном расположении остаточных дисбалансов и дисбалансов от неустранимого прогиба для отдельных типов подшипников это может привести к выходу их из строя и необходимости аварийной остановки турбоагрегата для замены подшипников. Родственные проблемы рассмотрены в работах [510].

Для двух наиболее часто используемых типов подшипников – сегментных и эллиптических [4] сравним виброактивность системы “РВД–РСНД” турбоагрегата с кососимметричным остаточным дисбалансом на каждом из пролетов и неустранимым прогибом с наибольшим параметром формы [1].

Физическая модель и расчетная схема системы “РВД–РСНД”. Расчетная схема конструкции представлена на рис. 1. Она состоит из невесомого вала, имеющего изгибную жесткость, несущего на себе приведенные сосредоточенные массы mjj = 1, …, 9). Вал установлен на трех подшипниках. Оси Оx, Оy представляют собой главные оси жесткости наиболее нагруженного среднего подшипника: Оx – ось минимальной жесткости, Оy – ось максимальной жесткости.

Рис. 1.

Расчетная схема валопровода: $C_{{ii}}^{{(k)}}$, $B_{{ii}}^{{(k)}}$ – параметры жесткости и демпфирования подшипников, соответственно: k = 1, 2, 3 – номер подшипника слева направо, ii = 11 – для значений параметров в направлении оси Оx, ii = 22 – в направлении оси Оy.

Участок валопровода ротора высокого давления в процессе эксплуатации приобретает температурные изменения или изменения, обусловленные несоблюдением условий пуска турбоагрегата [13]. При этом оси РВД системы “РВД-РСНД” после установки в подшипники представляют собой плоские кривые, лежащие в плоскостях yОz или xОz. На рис. 2 представлена исходная форма оси, снятого с левой опоры валопровода системы “РВД-РСНД” (кривая О1ВС). Величина , где f – величина смещения левого конца РВД, ${\delta }$ – стрелка приобретенного прогиба РВД, может быть названа параметром формы кривой начального прогиба РВД.

Рис. 2.

Форма оси свободного от нагрузок валопровода системы “РВД-РСНД” с первоначально погнутым РВД, присоединенным муфтой к непогнутому РСНД.

Форма кривой прогиба не может быть предсказана заранее, так как зависит от большого числа различных и не известных факторов [1]. Из всех теоретически возможных форм неустранимого прогиба [1] для дальнейших расчетов выберем вариант, соответствующий наибольшему значению параметра формы α

(1)
$~w = f(\sin \pi z{\text{/}}2l - 1);\quad 0 < z < l,\quad \alpha = f{\text{/}}\delta = 4.75.$

Расчет вынужденных колебаний системы “РВД-РСНД” при выбеге. Значения физических параметров системы в расчетном примере соответствуют системе “РВД–РСНД” турбины К-300-23.5 производства ЛМЗ: длина каждого из пролетов l = 5.5 м; масса РВД М1 = 9600 кг, масса РСНД М2 = 17800 кг$;$ приведенные дискретные массы- m1 = 500 кг, m2 = 2800 кг, m3 = 3000 кг, m4 = 2800 кг, m5 = 1000 кг, m6 = 3000 кг, m7 = = 5000 кг, m8 = 8800 кг, m9 = 500 кг; приведенные изгибные жесткости – вала РВД $E{{I}_{1}} = 5.15 \times {{10}^{8}}$ Н/м, вала РСНД $E{{I}_{2}} = 18.025 \times {{10}^{8}}$ Н/м. Форма кривой неустранимого прогиба ротора высокого давления принята в соответствии с соотношением (1).

Для сравнения уровня виброактивности системы “РВД–РСНД” с кососимметричным дисбалансом выберем эллиптические и шестиколодочные сегментные подшипники скольжения, как наиболее распространенные в отечественной теплоэнергетике [4]. Схемы подшипников представлены на рис. 3а, б.

Рис. 3.

Основные схемы опорных подшипников: а – эллиптический, 1 – нижний вкладыш, 2 – верхний вкладыш, 3 – шип; б – сегментный шестиколодочный, 4 – колодки, 5 – корпус.

Значения параметров подшипников обоих типов для системы “РВД–РСНД” из рассматриваемого примера представлены в табл. 1.

Таблица 1.

Значения параметров подшипников

Вид Эллиптические Сегментные
№ подшипника, (№ массы) 1-(1) 2-(5) 3-(9) 1-(1) 2-(5) 3-(9)
${{С}_{{11}}}$, Н/м 0.022 × 109 0.157 × 109 0.114 × 109 0.127 × 109 0.172 × 109 0.164 × 109
${{С}_{{22}}}$, Н/м 0.685 × 109 1.208 × 109 1.041 × 109 0.343 × 109 1.368 × 109 1.038 × 109
${{B}_{{11}}}$, кг с–1 0.564 × 106 0.804 × 106 0.752 × 106 0.276 × 106 0.366 × 106 0.352 × 106
${{B}_{{22}}}$, кг с–1 3.256 × 106 5.605 × 106 4.944 × 106 0.440 × 106 2.453 × 106 1.917 × 106

На рис. 4 приведены линии неустранимого приобретенного прогиба системы “РВД–РСНД” с выбранным параметром формы прогиба РВД, установленной в подшипники двух рассматриваемых типов. При расчетах стрелка неустранимых прогибов в соответствии с нормами ISO на тепловую пробу была выбрана равной предельно допустимому значению ${\delta }$ = 20 мкм.

Рис. 4.

Линии неустранимого прогиба системы “РВД – РСНД” для рассматриваемого примера и варианты расположения неуравновешенностей: 1 – жесткие опоры; 2 – эллиптические подшипники, ось Ox; 3 – эллиптические подшипники, ось Oy; 4 – сегментные подшипники, ось Ox; 5 – сегментные подшипники, ось Oy. АБ – вариант расположения неуравновешенностей при наиболее виброактивном симметричном дисбалансе; Ф – вариант расположения неуравновешенностей при наиболее виброактивном кососимметричном дисбалансе.

Собственные частоты, проходимые турбоагрегатом при выбеге, и соответствующие им нормированные формы колебаний для системы с выбранными параметрами (табл. 1), были рассчитаны в [3]. У системы “РВД–РСНД” при обоих типах подшипников имеется по пять собственных частот, лежащих ниже рабочей скорости вращения ТА ω = 314.16 рад/с: три в направлении оси Оx и две в направлении оси Оy. Значения этих частот приведены в табл. 2 полужирным шрифтом.

Таблица 2.

Собственные частоты системы

Вид подшипников Эллиптические Сегментные
Порядковый номер собственной частоты по Х по Y по Х по Y
1-я, рад/с 100.8 182.7 110.9 178.7
2-я, рад/с 112.4 213.1 136.1 213.3
3-я, рад/с 242.8 538.7 295.0 510.1

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

(2)
$\begin{gathered} {\text{эллиптические}}\;{\text{подшипники}}\quad p_{{{\text{3x}}}}^{{{\text{эл}}}} = 242.8\;{\text{рад/с}}; \\ ~{\text{сегментные}}\;{\text{подшипники}}\quad p_{{{\text{3x}}}}^{{{\text{сег}}}} = 295.0\;{\text{рад/с}}. \\ \end{gathered} $

Соответствующие этим частотам нормированные формы колебаний представлены на рис. 5. Вынужденные колебания именно по этим формам могут возбуждаться при выбеге турбоагрегата при наличии кососимметричных дисбалансов обоих роторов и иметь при этом превышающие нормы ISO амплитуды в зоне подшипников. Это может привести к выходу подшипников из строя и аварийной остановке ТА.

Рис. 5.

Собственные формы колебаний, опасные для подшипников: 1${\varphi }_{{3x~}}^{{{\text{*эл}}}}$ – эллиптические подшипники, 2$~{\varphi }_{{3x~}}^{{{\text{*сег}}}}$ – сегментные подшипники.

Расчет амплитуд переходных колебаний системы “РВД–РСНД” в околорезонансных зонах частот (2) при кососимметричном дисбалансе проведен на основе методики, предложенной в [4] и ранее использованной в [1]

(3)
${{A}_{j}} = {{{\varphi }}_{{j(k)}}}({{p}_{{(k)}}}{\text{/}}{{H}_{{(kk)}}}){{{\chi }}_{{(k)}}},$
где Aj – амплитуда в точке z = zj расположения j-й сосредоточенной массы на оси ротора (рис. 1), ${{{\varphi }}_{{j(k)}}}$ значение k-й формы колебаний в точке z = zj (здесь и далее по индексам в скобках суммирование не проводится), p(k) – собственная частота колебаний, H(kk) – приведенный коэффициент демпфирования для колебаний по форме k, ${{{\chi }}_{{(k)}}}$$~$параметр возбуждения при колебаниях по форме k.

Приведенный параметр демпфирования H(kk) для модели “РВД–РСНД” выражается соотношением [4]

(4)
${{H}_{{(kk)}}} = {{b}^{1}}{\varphi }_{{1k}}^{2} + {{b}^{2}}{\varphi }_{{5k}}^{2} + {{b}^{3}}{\varphi }_{{9k}}^{2},$
где ${{{\varphi }}_{{1k}}}$, ${{{\varphi }}_{{5k}}}$, ${{{\varphi }}_{{9k}}}$ – значения главных форм колебаний ротора на трех опорах, расположенных в точках 1, 5, 9 на рис. 1 для каждой из плоскостей xОz и yОz, соответственно для каждого типа подшипников; ${{b}^{j}} = B_{{11}}^{j}$ – коэффициент демпфирования для подшипника с номером j каждого из рассматриваемых типов при колебаниях в плоскости xОz, ${{b}^{j}} = B_{{22}}^{j}$ – при колебаниях в плоскости yОz, (j = 1, 2, 3), (рис. 1).

Параметр меры возбуждения колебаний в системе “РВД–РСНД” при резонансе $~{{{\chi }}_{k}}$ включает в себя влияние неуравновешенностей, приобретенных за счет неустранимого прогиба (рис. 4), и остаточных неуравновешенностей масс с эксцентриситетами ${{e}_{j}}$j = 1, …, 9)

(5)
${{{\chi }}_{k}} = {\chi }_{k}^{f} + {\chi }_{k}^{\varepsilon };\quad {\chi }_{k}^{f} = f\sum\limits_{j = 1}^9 {{{m}_{j}}\tilde {u}({{z}_{j}})} {{{\varphi }}_{{jk}}};\quad {\chi }_{k}^{\varepsilon } = ~\sum\limits_{j = 1}^9 {{{m}_{j}}{{e}_{j}}{{{\varphi }}_{{jk}}}} {\text{.}}$
Здесь ${{{\varphi }}_{{jk~~}}}$ – собственные формы колебаний системы с соответствующим типом подшипников; $\tilde {u}({{z}_{j}})$ – значения функций неустранимого прогиба системы вдоль осей Ox и Oy (рис. 4) в точках крепления масс на валу ( j = 1, …, 9) для каждого типа подшипников. В соотношениях (4) и (5) для рассматриваемого примера следует положить k = 3.

В работах [1, 3] показано, что основной вклад в величину суммарного параметра возбуждения системы, заданного первым соотношением (5), дает параметр $~{\chi }_{k}^{{\varepsilon }}$, зависящий от распределения остаточных дисбалансов роторов по длине валопровода и вида соответствующей собственной формы колебаний. Из (5) следует, что наибольшее значение суммарного параметра возбуждения ${{{\chi }}_{k}}$ достигается, если ${\chi }_{k}^{{\varepsilon }}$ и ${\chi }_{k}^{f}$ имеют одинаковые знаки. Далее предполагаем, что для обоих типов подшипников вектора остаточных дисбалансов и кривые неустранимого прогиба лежат в одной плоскости. Тогда при неустранимом прогибе РВД, заданном уравнением (1), максимальное значение параметра ${{{\chi }}_{k}}$ достигается при расположении неуравновешенностей на пролетах, отвечающих варианту типа Ф (рис. 4). Здесь же приведено расположение неуравновешенностей для наиболее виброактивного на этих частотах симметричного точечного дисбаланса (вариант АБ) [3]. Из-за отсутствия норм на величину предельно допустимого момента инерционных сил, примем неуравновешенности EI (для первого пролета валопровода) и EII (для второго) при кососимметричном дисбалансе равными предельно допустимым их значениям при симметричном точечном дисбалансе [3]. Тогда для неуравновешенностей (индексы 2, 4, 6, 8 соответствуют номерам масс на валопроводе) получим

(6)
$\begin{gathered} \left| {{{E}_{2}}} \right| = \left| {{{E}_{4}}} \right| = {{E}_{{\text{I}}}} = {{M}_{1}}{\varepsilon ;}\quad \left| {{{E}_{6}}} \right| = \left| {{{E}_{8}}} \right| = {{E}_{{{\text{II}}}}} = {{M}_{2}}{\varepsilon }; \\ \Phi {\text{:}}\quad {{E}_{2}} = - {{E}_{1}},\quad {{E}_{4}} = {{E}_{{\text{I}}}},\quad {{E}_{6}} = {{E}_{{{\text{II}}}}},\quad {{E}_{8}} = - {{E}_{{{\text{II}}}}}; \\ \end{gathered} $
где $\varepsilon = \left| \varepsilon \right| = {\text{0}}{\text{.1}}g{\text{/}}{{{\omega }}^{2}}$ – модуль приведенных остаточных эксцентриситетов [10], M1, M2 – массы РВД и РСНД, соответственно.

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

(7)
$\chi _{{3x}}^{{\varepsilon (\Phi )}} = - {{E}_{{\text{I}}}}{\text{[}}\left| {{{\varphi }_{{3x}}}(z{}_{2})} \right| + \left| {{{\varphi }_{{3x}}}({{z}_{4}})} \right|{\text{]}} - {{E}_{{{\text{II}}}}}{\text{[}}\left| {{{\varphi }_{{3x}}}\left. {(z{}_{6})} \right| + \left| {{{\varphi }_{{3x}}}\left. {({{z}_{8}})} \right|} \right.} \right.{\text{],}}$
где ${{{\varphi }}_{{3x}}}$ – главные формы колебаний при частотах (2) для каждого типа подшипников.

В табл. 3 представлены значения амплитуд перемещений в контрольных точках валопровода системы с кососимметричным дисбалансом Ф (рис. 4) и неустранимым прогибом РВД (1), при использовании эллиптических и сегментных подшипников с параметрами, приведенными в табл. 1. Для сравнения показаны амплитуды колебаний в тех же точках при симметричном точечном дисбалансе варианта АБ, наиболее виброактивном симметричном дисбалансе на рассматриваемых частотах [3].

Таблица 3.

Амплитуды колебаний системы “РВД-РСНД” от прогиба и дисбалансов

Амплитуды в контрольных точках: подшипниках и серединах пролетов (мкм)
№ подшипника и/или массы z = 0 (1); m1 z = 0.5l m3 z = l (2); m5 z = 1.5 l m7 z = 2 l (3); m9
ЭЛЛИПТИЧЕСКИЕ ПОДШИПНИКИ: p3x = 242.8 рад/c
А f –8.3 0.7 4.0 0.7 –3.7
(Ф) Аε –51.8 –4.8 33.0 6.6 –32.4
(Ф) А f+ Аε –60.1 –4.1 37.0 7.3 –36.1
(АБ) А f+ Аε –14.0 1.3 7.0 1.1 –6.5
СЕГМЕНТНЫЕ ПОДШИПНИКИ: p3x = 295 рад/c
А f –9.6 –0.9 6.1 1.2 –6.0
(Ф) Аε –121.6 –11.3 77.4 15.4 –76.1
(Ф) А f+ Аε –131.2 –12.2 83.5 16.6 –82.1
(АБ) А f+ Аε –120.0 –11.1 76.0 7.0 76.0

На рис. 6 приведены графики амплитуд в контрольных точках валопровода для двух типов подшипников при кососимметричном дисбалансе типа Ф и симметричном дисбалансе АБ.

Рис. 6.

Амплитуды колебаний системы “РВД – РСНД”, установленной на пошипниках двух различных типов: 1, 2 – сегментные подшипники; 3, 4 – эллиптические подшипники; сплошные линии – кососимметричный дисбаланс Ф; пунктирные – симметричный дисбаланс АБ.

Основные результаты и выводы.

1. Для модели двухпролетного ротора высокотемпературной части валопровода турбины типа К-300-23.5 (рис. 1) рассчитаны амплитуды переходных колебаний в контрольных точках на частотах ${{p}_{{3x}}}$ при выбеге системы с наиболее опасным (с точки зрения виброактивности) типе кососимметричного дисбаланса варианта Ф. Рассмотрены наиболее распространенные в турбостроении типы подшипников – эллиптические и сегментные. Частоты $p_{{{\text{3x}}}}^{{{\text{эл}}}}$ и $p_{{{\text{3x}}}}^{{{\text{сег}}}}$ (2) выбраны потому, что им соответствуют собственные формы колебаний валопровода с узлами на каждом из пролетов (рис. 5), что в сочетании с кососимметричным остаточным дисбалансом в системе может быть опасно для сохранности подшипников.

2. Проведено сравнение амплитуд в контрольных точках системы при ее установке на эллиптические и сегментные подшипники на частотах $p_{{{\text{3x}}}}^{{{\text{эл}}}}$ и $p_{{{\text{3x}}}}^{{{\text{сег}}}}$ при дисбалансе типа Ф с соответствующими амплитудами при осесимметричном наиболее виброактивном дисбалансе типа АБ (рис. 4). Из табл. 3 и рис. 5 следует, что амплитуды системы “РВД–РСНД” в точках установки ее на подшипники при кососимметричном дисбалансе Ф превосходят соответствующие амплитуды при симметричном дисбалансе АБ: для эллиптических подшипников примерно в 4.3 раза, для сегментных – в 1.1 раза.

3. Амплитуды перемещений системы “РВД–РСРД” с сегментными подшипниками в точках ее установки на подшипники более, чем в 2 раза, превосходят амплитуды той же системы с эллиптическими подшипниками. При этом для системы с сегментными подшипниками амплитуды в точке А рис. 2 примерно в 1.5 раза превышают норму ISO (90 мкм) и являются опасными для сохранности подшипников. Для системы с эллиптическими подшипниками эти амплитуды примерно в 1.5 раза ниже верхней границы нормы ISO и для подшипников безопасны.

4. Для системы “РВД–РСРД” с выбранными физическими параметрами использование эллиптических подшипников приводит к снижению ее виброактивности более, чем в 2 раза (табл. 3, рис. 4) по сравнению с сегментными. Применение эллиптических подшипников является безопасным с точки зрения выполнения норм ISO и поэтому наиболее целесообразным.

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

  1. Volokhovskaya O.A., Barmina O.V. Evaluation of the Influence of the Initial Deflection Curve Shape for a Two-Span Rotor with Residual Imbalance on Its Vibroactivity during Rundown // Journal of Machinery Manufacture and Reliability. 2017. V. 46. № 2. P. 99.

  2. Volokhovskaya O.A., Barmina O.V. The Effect of Residual Imbalance Type in Vibration Activity of Double-Span Rotor with Different Curvature of the Axis under Rundown // Journal of Machinery Manufacture and Reliability. 2018. V. 47. № 5. P. 403.

  3. Kostyuk A.G., Volokhovskaya O.A. Vibration activity evaluation of double-span rotor at rundown caused by its initial curvature and residual unbalanced // Thermal Engineering. 2017. V. 64. № 1. P. 37.

  4. Костюк А.Г. Динамика и прочность турбомашин. М.: Изд. МЭИ. 2000. 479 с.

  5. Рыгин Р.Е. Влияние некомпенсированного моментного дисбаланса дисков на состояние гибкого ротора // Казанская наука. 2011. № 2. С. 42.

  6. Jozef H., Vojtech S., Ales Ja. Predictive control of radial rotor vibrations // Мир транспорта и технологических машин. 2014. № 1 (44). С. 83.

  7. Бранцевич П.Ю. Анализ причин изменения вибрационного состояния агрегата роторного типа // Актуальные вопросы машиноведения. 2013. Т. 2. С. 277.

  8. Урьев Е.В., Мурманский Б.Е., Бродов Ю.М. Концепция системы вибрационной диагностики паровой турбины // Теплоэнергетика.1995. № 4. С. 36.

  9. Миндрин В.И., Пачурин Г.В., Ребрушкин М.Н. Виды и причины вибрации энергетических машин // Современные наукоемкие технологии. 2015. № 5. С. 32.

  10. Клейманов Р.В., Коршунов А.В. Расчет амплитуды колебаний роторов турбомашин при разгоне // Компрессорная техника и пневматика. 2016. № 3. С. 24.

Дополнительные материалы отсутствуют.