Проблемы машиностроения и надежности машин, 2019, № 5, стр. 10-17
ВЛИЯНИЕ ТИПА ПОДШИПНИКОВ НА ВИБРОАКТИВНОСТЬ ДВУХПРОЛЕТНОГО РОТОРА С КОСОСИММЕТРИЧНЫМ ОСТАТОЧНЫМ ДИСБАЛАНСОМ
О. А. Волоховская 1, *, О. В. Бармина 1
1 Институт машиноведения им. А.А. Благонравова
г. Москва, Россия
* E-mail: olgaavol@yandex.ru
Поступила в редакцию 24.04.2019
Принята к публикации 25.06.2019
Аннотация
Исследована зависимость виброактивности высокотемпературного отдела турбоагрегата, включающего в себя ротор высокого давления и ротор средне-низкого давления, при кососимметричном остаточном дисбалансе каждого из роторов от типа наиболее часто используемых подшипников. Рассмотрены шестиколодочные сегментные и эллиптические подшипники. При исследовании учтено наличие в системе неустранимого прогиба валопровода. Расчеты выполнены для турбины К-300-23.5 производства ЛМЗ. Приведен ключевой результат анализа.
Введение. Виброактивность турбоагрегата (ТА) определяется прежде всего величиной и распределением остаточного дисбаланса после балансировки роторов и дисбаланса, вызванного неустранимым прогибом, порожденным температурной нестабильностью роторов или несоблюдением условий пуска [1–3]. Кроме указанных факторов существенное влияние на уровень виброактивности оказывает тип подшипников, на которых установлен турбоагрегат [4]. Известно, что амплитуды колебаний в точках установки ТА на подшипники при рабочей скорости должны удовлетворять нормам ISO (International Standart Organisation): 1) $A \leqslant 90$ мкм; 2) в серединах пролетов не превышать величины зазоров в уплотнениях. Тип подшипников влияет на выполнение обоих условий ограничения амплитуд. Однако это влияние становится определяющим для выполнения первого условия для систем ротор высокого давления – ротор средне-низкого давления турбоагрегата (“РВД–РСНД” ТА) с кососимметричным остаточным дисбалансом каждого из роторов. При переходных колебаниях на выбеге турбоагрегата, возникающих на собственных частотах, которым соответствуют формы колебаний с близкими к середине пролетов узлами, наличие в системе кососимметричного дисбаланса способствует возрастанию амплитуд колебаний в точках установки системы на подшипники. При неблагоприятном взаимном расположении остаточных дисбалансов и дисбалансов от неустранимого прогиба для отдельных типов подшипников это может привести к выходу их из строя и необходимости аварийной остановки турбоагрегата для замены подшипников. Родственные проблемы рассмотрены в работах [5–10].
Для двух наиболее часто используемых типов подшипников – сегментных и эллиптических [4] сравним виброактивность системы “РВД–РСНД” турбоагрегата с кососимметричным остаточным дисбалансом на каждом из пролетов и неустранимым прогибом с наибольшим параметром формы [1].
Физическая модель и расчетная схема системы “РВД–РСНД”. Расчетная схема конструкции представлена на рис. 1. Она состоит из невесомого вала, имеющего изгибную жесткость, несущего на себе приведенные сосредоточенные массы mj ( j = 1, …, 9). Вал установлен на трех подшипниках. Оси Оx, Оy представляют собой главные оси жесткости наиболее нагруженного среднего подшипника: Оx – ось минимальной жесткости, Оy – ось максимальной жесткости.
Участок валопровода ротора высокого давления в процессе эксплуатации приобретает температурные изменения или изменения, обусловленные несоблюдением условий пуска турбоагрегата [1–3]. При этом оси РВД системы “РВД-РСНД” после установки в подшипники представляют собой плоские кривые, лежащие в плоскостях yОz или xОz. На рис. 2 представлена исходная форма оси, снятого с левой опоры валопровода системы “РВД-РСНД” (кривая О1ВС). Величина , где f – величина смещения левого конца РВД, ${\delta }$ – стрелка приобретенного прогиба РВД, может быть названа параметром формы кривой начального прогиба РВД.
Форма кривой прогиба не может быть предсказана заранее, так как зависит от большого числа различных и не известных факторов [1]. Из всех теоретически возможных форм неустранимого прогиба [1] для дальнейших расчетов выберем вариант, соответствующий наибольшему значению параметра формы α
Расчет вынужденных колебаний системы “РВД-РСНД” при выбеге. Значения физических параметров системы в расчетном примере соответствуют системе “РВД–РСНД” турбины К-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а, б.
Значения параметров подшипников обоих типов для системы “РВД–РСНД” из рассматриваемого примера представлены в табл. 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 мкм.
Собственные частоты, проходимые турбоагрегатом при выбеге, и соответствующие им нормированные формы колебаний для системы с выбранными параметрами (табл. 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 амплитуды в зоне подшипников. Это может привести к выходу подшипников из строя и аварийной остановке ТА.
Расчет амплитуд переходных колебаний системы “РВД–РСНД” в околорезонансных зонах частот (2) при кососимметричном дисбалансе проведен на основе методики, предложенной в [4] и ранее использованной в [1]
где 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},$Параметр меры возбуждения колебаний в системе “РВД–РСНД” при резонансе $~{{{\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{.}}$В работах [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} $Параметр возбуждения от наиболее виброактивного кососимметричного дисбаланса типа Ф на частотах, заданных равенствами (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{],}}$В табл. 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 приведены графики амплитуд в контрольных точках валопровода для двух типов подшипников при кососимметричном дисбалансе типа Ф и симметричном дисбалансе АБ.
Основные результаты и выводы.
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 и поэтому наиболее целесообразным.
Список литературы
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.
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.
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.
Костюк А.Г. Динамика и прочность турбомашин. М.: Изд. МЭИ. 2000. 479 с.
Рыгин Р.Е. Влияние некомпенсированного моментного дисбаланса дисков на состояние гибкого ротора // Казанская наука. 2011. № 2. С. 42.
Jozef H., Vojtech S., Ales Ja. Predictive control of radial rotor vibrations // Мир транспорта и технологических машин. 2014. № 1 (44). С. 83.
Бранцевич П.Ю. Анализ причин изменения вибрационного состояния агрегата роторного типа // Актуальные вопросы машиноведения. 2013. Т. 2. С. 277.
Урьев Е.В., Мурманский Б.Е., Бродов Ю.М. Концепция системы вибрационной диагностики паровой турбины // Теплоэнергетика.1995. № 4. С. 36.
Миндрин В.И., Пачурин Г.В., Ребрушкин М.Н. Виды и причины вибрации энергетических машин // Современные наукоемкие технологии. 2015. № 5. С. 32.
Клейманов Р.В., Коршунов А.В. Расчет амплитуды колебаний роторов турбомашин при разгоне // Компрессорная техника и пневматика. 2016. № 3. С. 24.
Дополнительные материалы отсутствуют.
Инструменты
Проблемы машиностроения и надежности машин