CC BY
65
Алокова М.Х., Культербаев Х.П.
ИЗГИБНЫЕ КОЛЕБАНИЯ ВЕРТИКАЛЬНОГО СТЕРЖНЯ ПЕРЕМЕННОГО СЕЧЕНИЯ С СОСРЕДОТОЧЕННОЙ МАССОЙ
Alokova M.Kh., Kulterbaev Kh.P.
FLEXURAL VIBRATIONS OF THE VERTICAL ROD VARIABLE SECTION WITH A CONCENTRATED MASS
Аннотация. Рассматриваются изгибные колебания вертикального стержня переменного сечения с сосредоточенной массой. Изучаются колебания, возбуждаемые векторным гармоническим и случайным процессами с компонентами в виде кинематических и динамических возмущений. Определены амплитуды и среднеквадратические отклонения. Доказано, что влияние корре-лированности компонентов случайных возмущений на колебания аналогично влиянию сдвига фаз на гармонические колебания. Превосходство среднеквад-ратических отклонений над амплитуд объясняется наличием первой резонансной частоты в спектральной плотности сонаправленных возмущений.
Ключевые слова: колебания вертикального стержня, векторный процесс возмущений, спектральная матрица.
Abstract. Flexural oscillations of the vertical rod of variable section with the concentrated weight are considered. The oscillations excited vectorial harmonic and accidental by processes with components in the form of kinematic and dynamic perturbations are studied. Amplitudes and mean squared deviations are determined.
It is proved that influence of correlation of components of accidental perturbations on oscillations is similar to influence of phase shift on harmonic oscillations. Superiority of mean squared deviations over amplitudes is explained by existence of the first resonance frequency in spectral density the sonapravlennykh of perturbations.
Key words: oscillations of the vertical rod, vectorial process of perturbations, spectral matrix.
Введение. Высокие вертикальные гибкие стержни (рис.1) широко используются в строительстве в качестве опор линий передач, мачт, прожекторов, антенных устройств и т.д. Постановка краевой задачи о свободных поперечных колебаниях вертикальных стержней постоянного сечения с учётом собственного веса, но без учёта сил сопротивления была рассмотрена в [1]. При этом было отмечено, что найти собственные значения дифференциальных уравнений с переменными коэффициентами в замкнутой форме не удаётся.
В работах [2, 3] приближённым методом Релея -Ритца изучались родственные задачи для стержней постоянного сечения. Спектральная задача и вынуж-
денные поперечные колебания стержней постоянного и переменного сечений при отсутствии дискретных масс рассмотрена в [4, 5].
Постановка задачи. Вертикальный стержень (рис. 1а) длины I с погонной массой m(x), с осевым моментом инерции поперечного сечения J(x), с дискретной массой M0 и осевым моментом инерции 10, из материала с модулем упругости Е совершает плоские изгибные поперечные колебания, описываемые функцией и(х, t).
Рисунок 1 - Расчётная схема: а) общая схема, б) элемент стержня
Источником колебаний является векторный процесс возмущений
^)={^)Д2&), Н^)}.
Выделим элемент стержня (рис. 1б) с действующими силами: м - изгибающий момент, р - поперечная сила, N - продольная сила, Г4^)-поперечная равномерно распределённая нагрузка, р(х) = т(х - продольная неравномерно распределённая нагрузка, ё1 - дифференциал даламберовой силы инерции распределённой массы, ёЯ - дифференциал силы вязкого трения.
Следующие силовые факторы находятся в известных соотношениях с параметрами
Мх) = м0в +1 р(^, ё1 = тиёх, ёЯ = цти ёх .
(1)
Здесь и далее точкой над переменными обозначены производные по времени, ц - удельный коэффициент линейно -вязкого трения.
Пользуясь технической теорией изгибных колебаний балок и принципом Даламбера, получим дифференциальное уравнение частных производных гиперболического типа
Ъи"" + 2Ъ"и" + Би" +ри' + ти + цтИ = Г4, хе(0, I), 1 >-да, Ь = Б^ Б = Ъ" + N (2) Штрихи в верхних индексах означают дифференцирование по координате х. В математическую модель колебаний в дополнение к основному уравнению (2) должны включаться граничные и начальные условия. Для свободных и вынужденных установившихся колебаний, рассматриваемых ниже, начальные условия не требуются.
Свободные колебания. В этом случае внешние воздействия отсутствуют, в силу чего основное уравнение (2) становится однородным
0
Ъи"" + 2Ъ'и"' + Би" + ри' + ши + лти = 0, х е(0, /), г >-да (3)
Следует заметить, что условиям закрепления концов стойки, изображённым на рис. 1а, соответствуют следующие однородные граничные условия:
Ъ(0)и" (0,1>1ои'(0,1) = 0, [Ъ(0)и" (0,г)] ' ^(0,1) = о, (4)
и(/,1 ) = 0, и ' (/,1 ) = 0 (5)
Далее необходимо определить спектры собственных значений и функций задачи (3) - (5).
Пользуясь методом разделения переменных, представим и(х, г) как произведение
и(х,1 ) = Х(х)еХ1, (6)
где
Х = -р,+ (7)
- характеристический показатель, (ц, ш) - собственная пара, состоящая из коэффициента затухания и частоты свободных колебаний, ] - мнимая единица. Подстановка (6) в (4), (5) даёт
ЪХ"" + 2Ъ'Х" + БХ" + рХ ' + шХ(Х + л)Х = 0, х е(0, /) (8)
ъ ' (0)Х " (0)+Ъ(0)х "(0)-м0 х2х(0) = 0, Ъ(0)х " (0) - :0 х2х "(0) = 0, (9)
х' (/)=0, х(/ )=0. (10)
Определить X аналитическими методами не удаётся.
Выход состоит в применении метода конечных разностей [6] и переходе от аналитической модели (8) - (10) к численной.
Введём сетку с шагом Ь = // (п -1) и вместо функции непрерывного аргумента х(х) сеточную функцию дискретного аргумента у(х1 ) = у1, 1 - номер узла сетки. Таким образом, далее Х(х) будет представляться вектором У = {у1, у2,..., уп}. Замена производного с точностью о(ь2) в (8) по обычным правилам метода приведёт к алгебраическим уравнениям
г1У1-2 + с1У1-1 + ^У1 + е1У1+1 + ^+2 = 0, 1 = 3, 4,...,п - 2 , (11)
где обозначено
г1 = Ъ1 - ЬЪ1, с1 = -4Ъ1 + 2ЬЪ1 + Ь2Б1 - 0,5Ь3р1,
а1 = 6Ъ1 - 2Ь2Б1 + Ь4ш1Х(Х + л), е1 = -4Ъ1 - 2ЬЪ1 + Ь2Б1 + 0,5Ь3р1, = Ъ1 + ЬЪ1. Аналогичные замены произведём в граничных условиях (9), (10) и получим (12)
аУ1 +рУ2 + УУз +5У4 +еУ5 = 0, УУ1 +0у2 +куз + ^4 = 0, Уп-2 - 4уп-1 + 3уп = 0, Уп = 0
Здесь
а= 4Ъ1Ь - 5Ъ1 - 2М0 12Ь3, р = -10Ъ1Ь + 18Ъ1, у = 8Ъ1Ь - 24Ъ1, 5 = -2Ъ1Ь + 14Ъ1,
£=-3Ъ1, у = 4Ъ1 + 310Х2Ь, 0= -10Ъ1 - 4:0 12Ь, к = 8Ъ1 + 10Х2Ь, £ = -2Ъ1. Запишем (11), (12) как линейную однородную алгебраическую систему уравнений
о(х)у = 0. (13)
с(х) - квадратная разреженная матрица порядка п.
Элементы матрицы О являются функциями X, т. е. собственного значения задачи (7) - (9).
Условие существования нетривиального решения системы уравнений (13) даёт характеристическое уравнение
= 0, (14)
из которого определяются собственные значения, т. е. спектр
К, Х2,...}.
Уравнение (14) с учётом того, что его левая часть представляет комплексную функцию, можно переписать в виде ю)+jv2 (ц, ю) = 0.
Отсюда следует, что коэффициент затухания ц и частота свободных колебаний ю должны определяться из системы двух нелинейных алгебраических уравнений
^(ц, ю)= 0, V2(м,ю) = 0. (15)
Её решение будем находить с помощью метода координатного спуска [6].
Пример 1. Определим собственные значения и функции свободных колебаний стального стержня круглого поперечного сечения с параметрами: I = 18 м, п = 601, л = 0,02 с -1, E = 210 ГПа, M0 = 1000 кг, J0 = 200 кгм2, m(x) = 16 + 0,7х кг/м, I = ax + c м4, а = 0,0000002 м3, c = 0,0000007 м4.
Метод координатного спуска дал следующие значения собственных значений
ю = {8,90, 26,30 52,96} с , ц = 0,01 с .
Спектр собственных векторов, может быть найден с помощью системы уравнений (13). Первые три показаны на рис. 2.
О -Г^Г"-
4.------------:-----
е--------—
ю - - - - - -- -:----------/
Л2 --(■---------
14- - --Х-------_ _ __.
16----------
_ V _
М-1 о 1
Рисунок 2 - Первые собственные функции колебаний
Данные графики не похожи на приводимые в литературе для консольных балок и вертикальных стержней постоянного сечения, не несущих сосредоточенные массы. Причина в том, что инертная дискретная масса на верхнем конце создаёт эффект закрепления и собственные функции становятся подобными аналогичным кривым для стержня с неподвижными концами.
Вынужденные колебания при гармонических возмущениях. Теперь основное уравнение и граничные условия в совокупности образуют задачу
Ъи"" + 2Ъ'и'" + Би" + ри' + ти + лти = £4 , х е (0, /), г>-да. (16)
Ъ(0)и" (0,г)-д 0И "(0,г) = 0, [Ъ(0)и"(0,г)Г-М0и(0,г) = £1(1), и(/,г)=£2(1), и"(/,1)=£3(1). (17) Четыре возмущения являются компонентами векторного гармонического процесса
)£) £2(г) £3(1) £4(г)}, которые имеют частоты 0к и начальные фазы
£к (г) = аке-!(ок1к), к = 1,2,3, 4 . (18)
Здесь ак - вещественная амплитуда гармонических возмущений. Введём обозначение для комплексной амплитуды
Ак = аке-^к и упростим описание возмущений
£к (г) = Аке->°к1, к = 1,2,3,4. (19)
Выходной процесс и(х, г) в общем случае не будет гармоническим и, даже, периодическим. В то же время он будет суммой трёх гармоник с разными частотами 0к.
Периодическими такие колебания будут в том случае, если отношения 0к (] ф к) окажутся рациональными числами. Если 01 = 02 = 03 = 04 = 0, то выходной процесс будет гармоническим. Для этой задачи можно определять, кроме решения и(х, г), и амплитуду колебаний. Рассмотрим данный случай подробно.
Решение задачи (16), (17) ищется с помощью метода разделения переменных как произведение
и(х,г) = Х(х)еХг, (20)
где
Х = (21)
- характеристический показатель, х(х) - функция амплитуды вынужденных колебаний, подлежащая определению.
Подставим (20) в (16), (17) и получим
ЪХ + 2Ъ"X "" + Б1Х " + рХ " + тЛX = А4, Л = Х(Х + л), х е(0, /). (22)
ъ " (0)Х "(0)+Ъ(0)х ""(0)-м0х2х(0) = 0, Ъ(0)Х"(0)- 10х2х "(0) = А:, (23)
Х(/ )= А2, X " (/ )= А3. ( )
От дифференциального уравнения (22) перейдём к конечноразностной схеме. Можно показать, что будет получено конечноразностное уравнение
ЧУ1-2 + с1У1-1 + ФУ1 + е1 У1+1 + ^У1+2 = Я1, Я1 = А4ь4, 1 = 3, 4, ..., п - 2 , (24)
где коэффициенты имеют прежние значения, но с учётом (21). Граничные условия (23) приведём в соответствие с методом конечных разностей и запишем:
- на верхнем конце:
УУ1 +0У2 +КУ3 +^У4 = 0, ау1 +ру2 +УУ3 +5у4 +еу5 = Я2, Я2 = 2А^3. (25)
- на нижнем конце:
Уп-2 - 4Уп-1 + 3Уп = qn-1, Яп-1 = 2ЬА3 , Уп = Яп , Яп = А2 (26)
81
Уравнения (24) - (26) образуют неоднородную алгебраическую систему, которую можно переписать в матрично -векторной форме
сМУ = я, (27)
где - квадратная матрица порядка п,
я = {0,Я2 Яз.-Лп-^п-ъЯп}•
Система уравнений даёт вектор У и амплитуду
аи (х1 )= |У| •
Пример 2. Рассмотрим действие гармонических возмущений Г1 (1), ^(1), Г3 (1), Г4 (1) на вертикальный стержень с параметрами примера 1.
В качестве вещественных амплитуд возмущений возьмём
а1 = 50 Н, а2 = 0,05 м, а3 = 0,001 рад, а4 = 15 Н/м.
Изучим влияние возрастающих значений частоты возмущений
□ ={0 0,6 1 4,9 19 25} с-1
на амплитуду колебаний при отсутствии сдвига фаз, т. е. при
у =^2 = Уз =^4 = 0 •
Результаты счёта показаны кривыми на рис. 3. Кривая 1 соответствует статическому нагружению, кривые 2, 3 - возмущениям с весьма малыми частотами. Эти колебания происходят по форме статических прогибов.
Далее рост частот возмущений (кривая 4) приближает колебания к резонансным, амплитуды возрастают, мода колебаний становится подобной первой собственной форме на рис. 2. Далее колебания происходят уже по второй и третьей собственным функциям (кривые 5, 6).
Рисунок 3 - Амплитуды колебаний при возрастании частот возмущений
Пример 3. Рассмотрим колебания вертикального стержня при различных сдвигах фаз возмущений, но одинаковых частотах
ш = 4 с -1 •
Сходная задача рассматривалась в работе [7].
Изучим влияние начальных фаз у на амплитуду колебаний. Существенное значение будут иметь не величины начальных фаз, а разницы между ними (сдвиги фаз). Наиболее характерными будут три случая:
Т = [(0,0,0,0); (0, -л, -л, 0); (0, -л/2, -л/2, 0)
Результаты счёта по компьютерной программе представлены на рис. 4.
Рисунок 4 - Амплитуды колебаний при сдвигах фаз возмущений Первая из кривых соответствует синфазным возмущениям. Поэтому амплитуды имеют наибольшие значения. Во втором случае (кривая 2), кинематические возмущения нижнего конца направлены против динамических нагрузок, т. е. отстают от них на %. Поэтому амплитуды существенно уменьшились. При этом по очевидной причине нижний и средний участки стержня колеблются в противофазе. В третьем варианте (кривая 3), амплитуда принимает промежуточные значения.
Вынужденные колебания при случайных возмущениях. Зачастую источники колебаний являются случайными (сейсмика, атмосферная турбулентность и т.д.). Поэтому перейдём теперь от детерминистического варианта краевой задачи (16), (17) к стохастической. В таком случае векторный процесс возмущений ед будет стационарным, стационарно связанным центрированным случайным процессом с заданной спектральной матрицей, обладающей свойством эрмитовости
(
^ И =
811 ®12 813
Л
>21
5 22
5 23
(28)
V831 832 833)
Здесь и далее звёздочка означает переход к комплексно -сопряжённым величинам. Теперь и(х, 1:) будет пространственно-временным случайным полем.
Далее будем изучать стационарный случайный процесс колебаний, поэтому начальные условия не потребуются. Определим по заданной спектральной матрице (28) спектральную плотность 8и(х, ю) случайного поля отклонений стержня и дисперсию. Для определения спектральной плотности 8и(х, ю) воспользуемся вектором передаточных функций Н(х, ¿а), элементы которого суть реакции стержня на гармонические возмущения с единичной амплитудой. Например, Нл(х, ¿а) является комплексной амплитудой колебаний стойки при
единичном возмущении = е^1*. Спектральная плотность перемещений, а далее и дисперсия определяются по известным соотношениям
Би(х, ш) = Нт(х, ш) Б^ш )Н*, Эи(х) = {ви(х,ш) ёш.
(29)
Возьмём компоненты вектора в виде процессов со скрытой периодичностью (характерными частотами). В таком случае элементы спектральной матрицы можно представить в виде
) = 2 ^ Р?к1 ^
л[(ю2 -е2!)2 + 4а21 ю2]
ек/ - ак/ + Рк/:
к, / = 1, 2, 3, 4,
где ак, рк/ - параметры широкополосности и характерной частоты, рк -элементы неотрицательно определённой корреляционной матрицы, | рк | < 1, ркк=1; Ок, О/ - среднеквадратические отклонения возмущений.
Интегрирование в правой части (29) выполняется численными методами. Пример 4. Рассмотрим случайные колебания стержня из предыдущих примеров при следующих значениях параметров случайных возмущений в системе единиц СИ:
су - 0,01, в - {0 0,6 1 4,9 19 25}, % - {50 0,05 0,001, 15}, р - 1, к,1 - 1,2,3,4.
Заметим, что такие входные данные делают данную задачу стохастическим аналогом детерминистической задачи примера 2.
Результаты вычислений представлены на рис. 5. Сравнение кривых рис. 3 и рис. 5 обнаруживает аналогию между случайными и гармоническими колебаниями. Трансформации кривых при росте частот возмущений примерно одинаковые. Как и в гармонических колебаниях, рост характерной частоты возмущений приводит к понижению среднеквадратических отклонений (кривые 1, 2, 3). При малых значениях р к1 колебания происходят по форме статических прогибов. Далее при повышении частоты они сменяются колебаниями, по форме подобными второй и третьей собственным формам (кривые 4, 5, 6).
Рисунок 5 - Среднеквадратические отклонения при возрастании характерных частот случайных возмущений Сравнение двух рисунков обнаруживает характер влияния параметра широкополосности ак на случайные процессы возмущений. Из -за наличия шумовых составляющих, определяемых ими, в случайном процессе содержатся частоты, совпадающие с собственными частотами. Поэтому среднеквадратические отклонения получаются несколько большими, чем соответствующие амплитуды. По аналогии с изучением роли сдвига фаз возмущений в детерминистической задаче теперь рассмотрим влияние коррелированности компонентов вектора случайных возмущений на колебания. С этой целью зафиксируем параметры характерных частот и широкополосности на значениях
(1 1 1 ( 1 -1 -1 (1 0 0 1 ^
1 1 1 1 Р 2 N -1 1 1 -1 , Рз = 0 1 1 0
1 1 1 1 -1 1 1 -1 0 1 1 0
I1 1 1 ъ 1 1 -1 -1 ъ 11 0 0 1J
Р^ = 4 с-1, а^ = 0,01 с-1, к = 1, 2, 3, 4. Три случая сдвига фаз возмущений, имевшие место в детерминистической задаче примера 3 промоделируем с помощью трёх соответствующих нормированных корреляционных матриц
Р| =
Результаты вычислений представлены на рис. 6. Номера кривых соответствуют номерам корреляционных матриц. Сравнение кривых рис. 4 и рис. 6 показывает их сходство.
Влияние коррелированности компонентов случайных возмущений на колебания аналогично влиянию сдвига фаз на гармонические колебания.
Некоторое превосходство среднеквадратических отклонений над амплитудами (кривые 1, 3) объясняется наличием первой резонансной частоты в спектральной плотности сонаправленных возмущений.
Обратное имеет место при противодействии возмущений (кривые 2). Здесь таже причина создаёт усиление взаимопомех.
3 л
2
Рисунок 6 - Влияние коррелированности случайных возмущений на средне-
квадратические отклонения
Вывод:
1. Наличие дискретной массы на верхнем конце вертикального стержня создаёт эффект его «закреплённости» при колебаниях.
2. При низкочастотных колебаниях вертикальных стержней формы колебаний совпадают с формами статических прогибов.
3. При гармонических колебаниях формы и амплитуды колебаний стержня существенно зависят не только от частоты возмущений, но и от сдвига их фаз.
4. При случайных колебаниях формы и среднеквадратические отклонения прогибов колебаний стержня существенно зависят не только от доминирующей частоты возмущений, но и от степени коррелированности компонентов векторного процесса возмущений.
Библиографический список:
1. Коллатц Л. Задачи на собственные значения (с техническими приложениями). М.: Наука. 1968. 504с.
2. Kim C.S., Dickinson S.M. On the analysis of laterally vibrating slender beams subject to various complicating effects. J. Sound and Vibr.. 1988. 122, N 3, рр. 441455.
3. Auciello N. M. A study of the numerical convergence of Rayleigh-Ritz method for the free vibrations of cantilever beam of variable cross-section with tip mass. Eng. Trans.. 1996. 44, N 3-4, рр. 375-388.
4. Алокова М.Х. Свободные колебания вертикальной стойки переменного сечения. XI Всероссийский съезд по фундаментальным проблемам теоретической и прикладной механики: сборник докладов. (Казань, 20-24 августа 2015 г.). Казань: Изд-во Казан. ун-та, 2015. С. 135-137.
5. Культербаев Х.П. Колебания вертикальной стойки переменного сечения при гармонических и случайных векторных возмущениях. XI Всероссийский съезд по фундаментальным проблемам теоретической и прикладной механики: сборник докладов. (Казань, 20-24 августа 2015 г.). Казань: Изд -во Казан. ун-та, 2015. С. 2181-2184.
6. Вержбицкий В.М. Основы численных методов. М.: Высшая школа, 2002. 840 с.
7. Культербаев Х.П. Колебания дискретно -континуальной балки, кинематически возбуждаемые гармоническими и случайными векторными процессами // Изв. вузов. Сев.-Кавк. регион. Техн. науки. 2012. № 3. С. 54-59.
References:
1. L. Collatz is the Task on eigenvalues (with technical Annex-mi). Moscow: Nauka. 1968.504p.
2. Kim C. S., Dickinson S. M. On the analysis of laterally vibrating slender beams subject to various complicating effects. J. Sound and Vibr. 1988. 122, N. 3, pp. 441455.
3. Auciello N. M. A study of the numerical convergence of Rayleigh-Ritz method for the free vibrations of cantilever beam of variable cross-section with tip mass. Eng. Trans.1996.44, No.3-4, pp.375-388.
4. Alakova M. H. Free vibrations of a vertical rack of variable cross section. XI all-Russian Congress on fundamental problems of theoretical and applied mechanics: collection of papers. (Kazan, August 20-24, 2015). Kazan: Publishing house of Ka-zan.Universitypress, 2015. pp. 135-137.
5. Kulterbaev H. P. Upright Vibrations of variable cross section under harmonic and random vector perturbations. XI All-Russian Congress on fundamental problems of theoretical and applied mechanics: collection of papers. (Kazan, August 20-24, 2015). Kazan: Publishing house of Kazan. University press, 2015 - pp. 2181-2184.
6. Verzhbitsky V. M. Fundamentals of numerical methods. M.: Higher school, 2002. 840p.
7. Kulterbaev H. P. Oscillations of discrete-continuous beam is kinematically excited by harmonic and stochastic vector processes, Izv. universities. North.- Caucasus region. Tekhn. sciences.- 2012. No.3. pp. 54-59.
