Determination of Safety rotor Displacement zone For ACTIVE magnetic bearing reaction wheel

Full Text

Введение. В настоящее время активные магнитные подшипники (АМП) находят все большее применение в различных областях промышленности. АМП широко используются в газоперекачивающих агрегатах, в токарных, фрезеровальных, шлифовальных станках, в различных высокоскоростных роторных машинах [1-4]. Перспективно использование АМП и в космической технике, в частности, в качестве опор роторов электродвигателей-маховиков (ЭДМ) систем ориентации и стабилизации космических аппаратов [3; 5]. Основные преимущества АМП заключаются в отсутствии механического контакта между вращающимися и неподвижными элементами ЭДМ в штатном режиме работы, а также в возможности управления жесткостью и демпфированием системы АМП непосредственно в процессе работы ЭДМ. Бесконтактная левитация ротора в магнитном поле исключает механический износ элементов ЭДМ, снижает моменты трогания и сопротивления вращению, исключает необходимость в смазке, а также обеспечивает увеличение срока службы и возможность повышения рабочего диапазона частоты вращения ротора. Варьирование упруго-демпферными характеристиками магнитного подвеса позволяет значительно снизить вибрацию и связанные с ней помеховые силы и моменты сил, обусловленные остаточной динамической неуравновешенностью ротора ЭДМ. Однако, несмотря на достаточно высокую надежность систем АМП, при разработке ЭДМ с магнитными подшипниками необходимо принимать во внимание возможность их отказа. Исходя из этого, любые роторные машины с магнитными опорами содержат узел страховочных подшипников (СП). В соответствии с ГОСТ Р ИСО 14839-1-2011 [6] СП представляет собой вспомогательный подшипник в системе АМП, предназначенный для ограничения перемещений ротора и предотвращения его контакта с поверхностью статора АМП. Помимо предотвращения контакта ротора со статорными элементами ЭДМ в нештатных режимах, СП обеспечивают выбег вращающегося ротора либо до его полной остановки, либо до возобновления работы системы АМП. В качестве СП преимущественно используются подшипники качения либо скольжения [3; 7]. Одним из основных параметров СП является зазор - половина разности между внутренним диаметром радиального СП и внешним диаметром цапфы ротора в этом подшипнике (радиальный зазор) или расстояние между торцевой поверхностью упорного СП и заплечиком вала (осевой зазор) [6]. Обычно зазоры в СП составляют примерно половину зазоров в АМП [8; 9]: δр.СП = 0,5δр.АМП; δо.СП = 0,5δо.АМП, где δр.СП - зазор в радиальном СП (допустимое перемещение ротора в радиальном направлении, рис. 1); δр.АМП - зазор в радиальном АМП; δо.СП - зазор в осевом СП (допустимое перемещение ротора в осевом направлении); δо.АМП - зазор в осевом АМП. Однако при использовании системы АМП в ЭДМ введения узла СП по ряду причин недостаточно для обеспечения надежной работы и требуемого срока службы прибора. Во-первых, воздействие внешних механических нагрузок при транспортировке ЭДМ, выводе космического аппарата на орбиту вызывает интенсивные колебания ротора в СП, результатом которых являются ударные нагрузки, действующие на СП и приводящие к их разрушению. Во-вторых, в процессе орбитальной эксплуатации космического аппарата возникают нештатные ситуации, обусловленные внезапным отключением напряжения питания либо потерей устойчивости системы АМП. В таком случае ротор ЭДМ ляжет на СП, в которых будет осуществлять дальнейший выбег. Режим выбега ротора в СП зависит от силы трения между СП и вращающимся ротором [9], а также от отношения k динамических P и статических Q сил, действующих на ротор [10]: k = P/Q. Наиболее благоприятным является режим маятниковых колебаний, который соответствует k < 1, т. е. случаю, когда статические силы, действующие на ротор, превышают динамические. Вместе с тем при выбеге ротора в СП в процессе орбитального функционирования космического аппарата статические силы Q, действующие на ротор, практически отсутствуют (вследствие невесомости), следовательно, k > 1. При этом происходит полный обкат ротором внутренней поверхности СП в направлении, противоположном вращению. Данный режим крайне нежелателен, поскольку вызывает значительную нагрузку на элементы СП, повышенную вибрацию ротора, что в итоге приводит к интенсивному износу и повреждению ротора и СП [9; 10]. Кроме того, зазоры в СП снижают жесткость опорного узла, уменьшая значения собственных частот системы «ротор-СП», что усиливает вибрационные явления, особенно если в процессе выбега ротор проходит через критические частоты. Также вследствие отсутствия надежных методик расчета и оценки работоспособности СП при аварийных выбегах существуют значительные сложности при проектировании СП и моделировании их работы [3; 10-14]. Указанные факторы приводят к тому, что СП могут выдержать аварийный выбег ротора крайне ограниченное количество раз. Суммируя вышесказанное, необходимо отметить, что при использовании системы АМП в ЭДМ, помимо СП требуется устройство арретирования многократного действия, осуществляющее автоматическую механическую фиксацию ротора относительно корпуса ЭДМ как при его наземной транспортировке и выводе космического аппарата на орбиту, так и при орбитальном функционировании. Устройство арретирования. Для фиксации ротора относительно корпуса ЭДМ предложено устройство арретирования [15; 16], совмещенное с узлом СП, которые установлены на двух валах-фиксаторах, расположенных с обоих концов вала ротора. Валы-фиксаторы могут реверсивно перемещаться в осевом направлении на расстояние ±δф (рис. 2), обеспечивая либо жесткую механическую фиксацию ротора относительно корпуса ЭДМ, либо его бесконтактную левитацию в магнитном поле системы АМП. Преимуществом такого устройства является многократная автоматическая фиксация и центрирование ротора при возникновении нештатных режимов. В арретированном состоянии ротор вращается в СП аналогично традиционным ЭДМ с шарикоподшипниковыми опорами. Вибрации и нагрузки, действующие на СП в таком режиме, значительно меньше нагрузок, действующих в процессе выбега ротора в СП с зазором. Определение зоны допустимых перемещений ротора ЭДМ в АМП. С целью обеспечения безопасного режима работы ЭДМ с АМП необходимо определить требуемое осевое перемещение вала-фиксатора ±δф. С одной стороны, данное перемещение должно быть достаточным для того, чтобы ротор ЭДМ мог совершать колебания относительно своего номинального положения, вызванные остаточной динамической неуравновешенностью и различными внешними воздействиями. С другой стороны, должна быть полностью исключена возможность механического контакта и повреждения ротора и статорных элементов системы АМП и ЭДМ в целом. Перемещение вала-фиксатора ±δф определяется зоной возможных смещений ротора ЭДМ в активных устройствах системы АМП (радиальными δр.АМП и осевыми δо.АМП зазорами), расстояниями между данными устройствами и геометрической формой контактных поверхностей ротора и вала-фиксатора. Простейшим случаем смещения ротора от его номинального положения является линейное перемещение в осевом δZ и радиальном δX направлении, как показано на рис. 3. Рис. 1. Зазоры в СП и АМП: 1 - радиальный АМП; 2 - осевой АМП; 3 - ротор; 4 - радиальный СП; 5 - осевой СП C:\Поляков М.В\Статьи\2016\Вестник СибГАУ\Схема магнитного подвеса (2).png Рис. 2. Схема ЭДМ с системой АМП: 1 - ротор; 2, 3 - ротор и статор датчика положения ротора; 4, 5 - ротор и статор радиальной электромагнитной опоры; 6, 7 - статор и ротор осевой электромагнитной опоры; 8, 9 - датчики осевого и радиального зазора АМП; 10 - СП; 11 - вал-фиксатор Рис а б Рис. 3. Схема линейных перемещений ротора: разарретированное номинальное положение (а); радиальное δX и осевое δZ смещение ротора (б) В общем случае при одновременном перемещении по координатам Х и Z зазоры δр.СП, δо.СП в зависимости от перемещения δф и угла конуса контактной поверхности α определяются по выражению Из данного выражения можно получить частные случаи перемещения по одной из координат. Максимальное перемещение по координате Z (осевой зазор в СП δо.СП) при δХ = 0 δZ = δо.СП = δф. Максимальное перемещение по координате Х (радиальный зазор в СП δр.СП) при δZ = 0 . Осевые и радиальные зазоры в СП равны между собой только при α = 90°, в остальных случаях допустимое перемещение по одной из координат больше, чем по другой. Поскольку перемещение фиксатора δф является максимальным осевым зазором δо.СП, то целесообразно задавать его значение, выполняя условие δф < δо.АМП, и рассчитывать соответствующее допустимое радиальное перемещение δр.СП для заданного угла α. Например, зазоры δо.АМП и δр.АМП в системе АМП составляют 0,4 мм, диапазон перемещения фиксатора δф = 0,3 мм. Значения допустимого радиального перемещения δр.СП при различных углах α приведены на рис. 4, из которого можно сделать вывод, что углы α > 106º являются недопустимыми, поскольку радиальное перемещение в СП δр.СП превышает радиальный зазор в АМП δр.АМП. Следующим вариантом смещения ротора относительно своего номинального положения в системе АМП является его угловой поворот вокруг поперечной оси Y. В данном случае необходимо определить максимальный угол β, на который может повернуться ротор, а также соответствующее данному углу линейное смещение наружной поверхности ротора в радиальном направлении δXβ. Радиус окружности, по которой перемещается точка А, расположенная на контактной поверхности ротора на наибольшем удалении от его геометрического центра, определяется выражением , где L - расстояние от геометрического центра ротора до торца его вала (рис. 5); d - диаметр вала ротора. C:\Поляков М.В\Статьи\2016\Вестник СибГАУ\Рисунки\Рис. 4.png Рис. 4. Зависимость допустимого радиального перемещения δр.СП от осевого перемещения δо.СП при различных углах конуса контактной поверхности α C:\Поляков М.В\Статьи\2016\Вестник СибГАУ\Поворот.png Рис. 5. Угловой поворот ротора вокруг поперечной оси Уравнение прямой a, характеризующее контактную поверхность фиксатора, можно записать следующим образом: . (1) Координаты точки пересечения С контактной поверхности фиксатора и окружности, по которой перемещается точка А, определяются из решения уравнения . Данное уравнение приводится к квадратному, его решением является координата zС, соответствующая точке пересечения окружности, по которой перемещается точка А при повороте ротора, с прямой a, характеризующей контактную поверхность фиксатора. Координата xC точки пересечения рассчитывается по уравнению (1) подстановкой zC вместо z. Максимальный угол поворота можно определить по дуге окружности, ограниченной начальным положением точки А и точкой пересечения C. Вследствие малости перемещений дугу АС можно считать прямой, ее длина определяется из треугольника АВС: . Максимальный угол поворота (в градусах) . Радиальное смещение наружной поверхности вала ротора, обусловленное его поворотом на угол β, . Рассмотренные варианты линейных перемещений и угловых поворотов являются частными случаями отклонения ротора от номинального положения. Естественно, в процессе работы ротор может совершать одновременно как линейные, так и угловые смещения. Однако в таком случае максимальное отклонение ротора от номинального положения не превышает значений δр.СП, δо.СП, β, рассчитанных ранее для частных вариантов, что возможно получить, учитывая соответствующие радиальные и осевые смещения ротора в уравнениях для расчета углового поворота вокруг поперечной оси. Таким образом, результатом геометрического решения задачи определения зоны допустимых перемещений ротора являются максимальные линейные смещения в радиальном и осевом направлениях δр.СП, δо.СП, а также значение максимального угла β его поворота вокруг поперечной оси. Динамическое поведение вращающегося ротора в магнитном подвесе. В процессе вращения главная центральная ось инерции ротора неизбежно совершает линейные и угловые колебания, основной причиной которых является его остаточная динамическая неуравновешенность, являющаяся суммой остаточных статической и моментной неуравновешенностей. Статическая неуравновешенность приводит к линейным колебаниям центра масс ротора в радиальном направлении, моментная - вызывает угловые колебания главной центральной оси инерции ротора. Для обеспечения нормальной работы ЭДМ амплитуды данных колебаний не должны превышать максимально допустимых перемещений в страховочных опорах δр.СП, δо.СП, β, рассчитанных при геометрическом решении задачи определения зоны допустимых перемещений ротора. Линейные и угловые колебания симметричного ротора как жесткого тела по осям X и Y (см. рис. 2) описываются следующими уравнениями [3]: (2) где m - масса ротора; B - коэффициент демпфирования электромагнитной опоры; C - коэффициент жесткости электромагнитной опоры; e - эксцентриситет ротора; ω - угловая скорость ротора; Jэ, Jо - экваториальный и осевой моменты инерции ротора; R - коэффициент углового демпфирования электромагнитной опоры; D - коэффициент угловой жесткости электромагнитной опоры; γ - угол, характеризующий моментный дисбаланс ротора. Решением системы уравнений (2) являются амплитуды вынужденных линейных и угловых колебаний ротора по осям X и Y. Амплитуда вынужденных колебаний геометрического центра ротора по любому из радиальных направлений, возникающих вследствие его статической неуравновешенности [3], . Амплитуда вынужденных угловых колебаний геометрической оси ротора [3] Амплитуда результирующих колебаний точки ротора, максимально удаленной от его геометрического центра, . Расчетные максимальные амплитуды колебаний ASmax и Фmax (рис. 6) не должны превышать максимально допустимых перемещений δр.СП и β, рассчитанных при решении геометрической задачи. В противном случае необходимо либо увеличить зазоры в СП (увеличить перемещение фиксатора δф), либо изменить параметры системы магнитного подвеса. C:\Поляков М.В\Статьи\2016\Вестник СибГАУ\Рисунки\Рис. 6 (а, б).png а б C:\Поляков М.В\Статьи\2016\Вестник СибГАУ\Рисунки\Рис. 6 (в).png в Рис. 6. Графики зависимости амплитуды радиальных (а), угловых (б) и результирующих (в) колебаний ротора от его угловой скорости ω Недостатком такого расчета является рассмотрение ротора как абсолютно жесткого тела, при этом не учитываются его собственные частоты и изгибные формы колебаний, а также деформации его элементов под действием центробежных сил, возникающих в процессе его вращения. Анализ вынужденных колебаний ротора и определение перемещений его элементов с учетом изгибных мод колебаний целесообразно проводить с помощью систем конечно-элементного анализа. Перед этим необходим частотный (модальный) анализ исследуемого ротора с целью выявления значений его собственных частот и соответствующих им изгибных форм колебаний, а также определения областей ротора, в которых наблюдаются максимальные смещения элементов от номинального положения [17]. Например, по результатам частотного анализа ротора ЭДМ с АМП (рис. 7) установлено, что на первой собственной частоте максимальное отклонение от номинального положения наблюдается в его ободе (рис. 8). При проведении анализа вынужденных колебаний ротора целесообразно разместить на ободе в точках максимальных смещений виртуальные датчики, регистрирующие результирующее перемещение данной точки и перемещения по отдельным координатам. Кроме того, датчики можно установить в других точках, например в области контактной поверхности вала ротора, для получения актуальной информации о смещениях ротора в выбранных областях. Применение виртуальных датчиков позволяет по результатам анализа вынужденных колебаний автоматически формировать амплитудно-частотную характеристику, соответствующую выбранной точке ротора (рис. 9). Максимальные перемещения Amax точек ротора, определенные при анализе вынужденных колебаний, не должны приводить к его касанию со статорными элементами ЭДМ. В противном случае необходимо либо скорректировать конструкцию ЭДМ (увеличить зазоры между ротором и статорами ЭДМ, изменить форму ротора с целью повышения его жесткости, увеличения собственных частот и уменьшения амплитуды его изгибных колебаний), либо изменить параметры системы магнитного подвеса. Заключение. Таким образом, разработанная методика определения зоны допустимых перемещений ротора, устанавливающая зависимости между перемещением вала-фиксатора ±δф, зазорами в страховочных подшипниках δр.СП, δо.СП., зазорами в активных устройствах системы магнитного подвеса δр.АМП, δо.АМП, учитывает характер динамического поведения вращающегося ротора в системе АМП. Использование данной методики позволяет определить оптимальное значение перемещения фиксаторов δф, обеспечивающее бесконтактное вращение ротора в магнитном поле, создаваемом системой АМП в штатном режиме работы, и исключить возможность контакта ротора со статорными элементами системы АМП и другими элементами ЭДМ в аварийных ситуациях. G:\Для статей!!!\Статьи\2016\Вестник СибГАУ\Рисунки\Сборка маховика кэа 2.png а б Рис. 7. 3D-модель (а) и конечно-элементная модель (б) ротора ЭДМ C:\Users\бук\Desktop\частота.png C:\Users\бук\Desktop\частота 2.png Рис. 8. Форма колебаний ротора ЭДМ, соответствующая его первой собственной частоте Рис. 9. Зависимость перемещения в точках на ободе ротора от частоты вращения ротора

About the authors

M. V. Polyakov

JSC “Scientific and Industrial centre “Polyus”; National Research Tomsk Polytechnic University

Email: info@polus-tomsk.ru
56v, Kirov Av., Tomsk, 634050, Russian Federation; 30, Lenin Av., Tomsk, 634050, Russian Federation

References

Supplementary files