Особенности моделирования роторов различного назначения.Объектом исследований в работе выбран ротор любой возможной конфигурациис валом, состоящим из участков кусочно-постоянного круглого(кольцевого) сечения. К такому представлению может быть приведен валпрактически любой роторной машины (турбины, компрессора, детандера,генератора и т.д.) [3]. Навесные элементы могут быть любой конфигурации,но для них должно быть известно местоположение центра масс навалу, масса, экваториальный и полярный моменты инерции. Применяемыев настоящее время подшипники (качения, скольжения, магнитные) могутмоделироваться как жестко защемленные, шарнирные или упругодемпферныеопоры в зависимости от степени близости к этим вариантамхарактеристик применяемого типа подшипников.Кроме того, для проведения расчетных исследований динамики должныбыть известны такие параметры ротора, как свойства материала вала, остаточныедисбалансы навесных элементов, разгонная характеристика и диапазонрабочих скоростей вращения ротора.Теоретические положения численного анализа роторной динамики.При анализе процессов, происходящих в роторной машине, следует разделятьпонятия критической скорости и резонанса. Первое связано с потерейустойчивости вращающегося ротора под действием возбуждения из-за собственнойнеуравновешенности [4-7]. Но в многовальных турбомашинах колебанияодного ротора могут быть вызваны неуравновешенностью другого [5,с. 167-170]. Такие колебания, в отличие от критических, называют резонансными[5]. Кроме того, резонансные вибрации определяются динамическимсостоянием машины в целом и могут возбуждаться, например, кинематическиили периодическими силами постоянного направления различной природы[4]. При этом ротор может совершать движение типа прямой или обратной,синхронной или несинхронной прецессии.В работе используются численный метод расчета динамических характеристик– метод конечных элементов (МКЭ) [8]. Определение критическихчастот ротора, как и все другие анализы, выполняется с использованиемуравнений роторной динамики, которые в матричной форме имеют вид [8]:[ M][U] + ([C] + [G])[U] + ([K] + [B])[U] = [F] , (1)где M, C и K – соответственно матрицы масс, демпфирования и жесткости, U– вектор-столбец узловых перемещений, F – вектор-столбец динамическойнагрузки, приведенной к узлам конечно-элементной сетки, а G и B – гироскопическаяматрица и матрица демпфирования, связанного с вращательнымдвижением, которые позволяют учесть зависимость динамических характеристик(например, собственных частот) от угловой скорости вращения ротора.Последняя матрица изменяет общую матрицу конструкционной жесткости,что может привести к неустойчивости движения ротора.Далее при [F] = 0 этим уравнениям придается форма, соответствующаяISSN 2078-9130. Вісник НТУ «ХПІ». 2012. № 55 (961) 119
задаче о собственных значениях, и для дискретных значений угловой скоростивращения в заданном диапазоне находятся собственные частоты, соответствующиеизгибным (поперечным) формам колебаний. Это позволяет построитьчастотную диаграмму (диаграмму Кэмпбелла) и по ней определитькритические скорости вращения с учетом гироскопического момента.На ротор при вращении действуют центробежные силы, проекции которыхна оси системы координат (y и z), перпендикулярные оси вала, изменяютсяпо гармоническому закону в противофазе [8]:FFyz2= ω (Fcosαcosωt+ Fsin αsinωt);2= ω (Fcosαsinωt+ Fsin α cosωt);F = me,где m – неуравновешенная масса, e – эксцентриситет, α – фаза неуравновешенности,ω – угловая скорость вращения.В комплексном виде равенства (2) могут быть записаны как:2Fy= ω (Fa− iFb)e ;(3)2iωtFz= ω ( −Fb− iFa)e .При выполнении гармонического анализа с учетом таких сил решаетсяуравнение вида (1) с правой частью (3) и строится амлитудно-частотная характеристика(АЧХ), например, для центров масс навесных элементов илиопорных участков. Это позволяет найти резонансные режимы с определениемамплитуд колебаний, построить траектории движения ротора, соответствующиеэтим угловым скоростям и в конечном итоге оценить их опасность.Интегрированное программное средство и методика вычислений.Для расчета динамических характеристик роторной системы (многоопорногоротора), то есть проведения расчетов собственных частот и форм невращающегосяротора, получения зависимостей собственных частот от угловой скоростивращения в заданном диапазоне и построения частотной диаграммыКэмбелла, выполнения гармонического анализа и построения амплитудночастотныххарактеристик с визуализацией траекторий движения ротора быласоздана программа «Sol». Она интегрирована с многоцелевым пакетом проектированияи конечно-элементного анализа. Геометрические параметрыротора в виде длин и диаметров отдельных участков вала задаются пользователемв главном окне программы «Данные», вид которого показан на рис. 1.Кроме того, в этом же окне имеется возможность указания мест расположениясосредоточенных масс с определением их параметров (масса, экваториальныйи полярный моменты инерции, дисбаланс), радиальных и осевыхопор с выбором вида опирания и заданием их характеристик (коэффициентовжесткости и сопротивления), свойств материала вала, а также типа расчета(критические скорости, гармонический анализ) и его параметров (количествособственных частот и форм, диапазон угловых скоростей для построениядиаграммы Кэмбелла, диапазон частот для построения АЧХ и количество120 ISSN 2078-9130. Вісник НТУ «ХПІ». 2012. № 55 (961)iωt(2)
- Page 1 and 2:
ISSN 2078-9130ВІСНИКНАЦІ
- Page 3:
Вісник Національно
- Page 6 and 7:
принят кандидатом
- Page 9 and 10:
А. М. Журавлевой и О.
- Page 11 and 12:
ции (1976 г.), орденом
- Page 13 and 14:
ук.- Х.: 1955. - 12 с. 4. Бо
- Page 15 and 16:
следующие формулы:
- Page 18 and 19:
m+ 1 tε ω+( ) ( ) .1 0,5∫ gmt
- Page 20 and 21:
ты количества движ
- Page 22 and 23:
туды импульсов U с ,
- Page 24 and 25:
напряжения на конд
- Page 26 and 27:
напряжений построе
- Page 28 and 29:
∂ u 1ε ij= u i j+ u j i+ u k iu,
- Page 30 and 31:
сти (6) выбираем зна
- Page 32 and 33:
Введение. Одним из
- Page 34 and 35:
∂ψ ∂ψx y ∂ψ ∂ψψ xy = +
- Page 36 and 37:
На рис. 2, а показано
- Page 38 and 39:
c = 10 5 Н/м (кривая 5) е
- Page 40 and 41:
которой величина н
- Page 42 and 43:
Рисунок 11 - Одиннад
- Page 44 and 45:
симметричной конст
- Page 46 and 47:
Рисунок 24 - Третий в
- Page 48 and 49:
Рисунок 36 - Шестой в
- Page 50 and 51:
Рисунок 48 - Девятый
- Page 52 and 53:
Рисунок 60 - Двенадц
- Page 54 and 55:
Рисунок 72 - Пятнадц
- Page 56 and 57:
variables functions and ANSYS. By t
- Page 58 and 59:
числения были повт
- Page 60 and 61:
Ψ - угол между осью x
- Page 62 and 63:
Вычисление микро н
- Page 64 and 65:
Табл. 4 и 5 показываю
- Page 66 and 67:
вен разности двух д
- Page 68 and 69:
безопасную работу
- Page 70 and 71: Из предыдущего опы
- Page 72 and 73: Однако, изменения т
- Page 74 and 75: УДК 539.1С. Н. ИСАКОВ,
- Page 76 and 77: абРисунок 1 - График
- Page 78 and 79: k1 = -0,1; k2 = 8,0167; k3 = -13,75
- Page 80 and 81: Рисунок 5 - Распреде
- Page 82 and 83: УДК 519:539:534С.В. КРАСН
- Page 84 and 85: жидкости. Схема при
- Page 86 and 87: Максимальные велич
- Page 88 and 89: 2( x,t)∂ uu( x,t)= 0, = 0 при
- Page 90 and 91: Тогда для изображе
- Page 92 and 93: УДК 539.3О. О. ЛАРІН, к
- Page 94 and 95: Рисунок 1 - Схема ко
- Page 96 and 97: Під час досліджень
- Page 98 and 99: Окрім СКЗ вібропри
- Page 100 and 101: Список литературы:
- Page 102 and 103: стях деформации (пр
- Page 104 and 105: гдеcrεij - тензор ско
- Page 106 and 107: абРисунок 5 - Перера
- Page 108 and 109: Выводы. Для изучени
- Page 110 and 111: нутого образца. В с
- Page 112 and 113: pгде εi- интенсивнос
- Page 114 and 115: = H11εr + H12ε θ ;θ H ε 21 r+
- Page 116 and 117: Кинетический закон
- Page 118 and 119: личение давления а
- Page 122 and 123: расчетных точек в э
- Page 124 and 125: сеточной дискретиз
- Page 126 and 127: Рисунок 4 - Поперечн
- Page 128 and 129: абРисунок 8 - Амплит
- Page 130 and 131: ции такого вида пол
- Page 132 and 133: моделирующих навес
- Page 134 and 135: J Ψ + c ϕazy∗2⎡ 1my + byy+ c
- Page 136 and 137: Структурная схема
- Page 138 and 139: ключить влияние ве
- Page 140 and 141: Общие соотношения
- Page 142 and 143: Колебания прямоуго
- Page 144 and 145: s K - положительные к
- Page 146 and 147: 2 ⎧ 4G( 0,0, − sK) = f2( δ )
- Page 148 and 149: ленном крае пласти
- Page 150 and 151: чия жесткого тела н
- Page 152 and 153: Чтобы определить н
- Page 154 and 155: ( n+1)1εK= arctg, (21)( n)χ ⎡*
- Page 156 and 157: Рисунок 1 - Графики
- Page 158 and 159: С целью дальнейшей
- Page 160 and 161: упругой анизотропи
- Page 162 and 163: тигло 56 %, а σ Т упал
- Page 164 and 165: УДК 539.3С. Ю. СОТРИХИ
- Page 166 and 167: Рисунок 2 - Блок-схе
- Page 168 and 169: УДК 539.1А. В. СТЕПУК,
- Page 170 and 171:
сти деформирования
- Page 172 and 173:
УДК 534.1:539.3А. Н. ШУПИ
- Page 174 and 175:
Рассмотрим цилиндр
- Page 176 and 177:
Рисунок 2 - Кристалл
- Page 178 and 179:
стержни отжигались
- Page 180 and 181:
с. 12. Шупиков А. Н. Не
- Page 182 and 183:
- неравномерность п
- Page 184 and 185:
Меридиональное сеч
- Page 186 and 187:
Рисунок 4 - Расчетны
- Page 188 and 189:
ЗМІСТКедровская О.
- Page 190:
НАУКОВЕ ВИДАННЯВІС