CC BY
17
УДК (539.3+624.07):534.1
Пространственные нелинейные колебания стержня с двумя жесткими заделками
Муницын А.И., канд. техн. наук, Крайнова Л.Н., соиск., Сабанеев Н.А., ст. преп.
Приводится аналитическое решение задачи о вынужденных изгибных колебаниях стержня с двумя заделками, имеющего близкие собственные частоты изгибных колебаний в двух взаимно перпендикулярных плоскостях вследствие близких значений главных осевых моментов инерции поперечного сечения. Учитывается геометрическая нелинейность, обусловленная изменением длины средней линии стержня при его пространственном движении. Исследуются колебания стержня в окрестности главного резонанса.
Ключевые слова: вынужденные колебания, собственные частоты, геометрическая нелинейность, главный резонанс.
Dimensioned Nonlinear Vibrations of Rod with Two Anchorages
A.I. Munitsyn, Candidate of Engineering, L.N. Krainova, Applicant, N.A. Sabaneev, Senior Teacher
The author presents an analytical solution of forced bending vibrations of a rod with two anchorages. The rod has close fundamental frequencies of forced vibrations in two mutually transverse plains due to close values of main axial torques of cross-section inertia. The authors consider geometrical nonlinearity, caused by length changes of the rod central line in its spatial motion. Rod vibrations are researched in the area of the main resonance.
Key words: forced vibrations, fundamental frequencies, geometrical nonlinearity, main resonance.
Пространственные колебания стержня с неподвижными шарнирными опорами и близкими значениями осевых моментов инерции исследованы в [1, 2]. При возбуждении колебаний в плоскости большей изгибной жесткости, в области резонанса обнаружены плоские и пространственные режимы колебаний стержня. При достаточно большом уровне возбуждения максимальные амплитуды реализуются именно на плоской форме колебаний. Этот участок амплитудно-частотной характеристики (АЧХ) изолирован и реализуется только при наличии внешних возмущений.
Постановка задачи. Пусть центральная ось стержня в недеформированном состоянии совпадает с осью х прямоугольной системы координат, главные оси инерции поперечного сечения параллельны осям у и г. Концам стержня соответствуют координаты х = 0 и х = Ь. Обозначим через V, w перемещения точек средней линии стержня. Пусть на стержень действуют внешние гармонические нагрузки во взаимно перпендикулярных направлениях q^(х, t) и q^(х, t), диссипацию учитываем по моделям внутреннего и внешнего трения. Уравнения пространственных нелинейных колебаний стержня в безразмерных переменных имеют вид [1, 2] a4V + vIV + р0й1У + а4ст0 V - 4у0е^" =
= qv (х )соэ(ц t + 01),
а4\л/ + cwIV + ср01/й1У + а4ст0 W - (1)
-4yoSoW" = qw (х )соэ(ц t + 02),
8о =11 (2 + W '2 )х.
0
Все перемещения и координата х отнесены к длине стержня £. Точками обозначено дифференцирование по безразмерному времени
t* = .
Здесь ю0 - первая собственная частота изгибных колебаний в плоскости Охг: ю0 = (а / L)2(Еиу)1/2(рР)-1/2; Е и р - модуль упругости и плотность материала стержня; Р -площадь поперечного сечения; иу, и; - моменты инерции сечения; е0 - деформация средней линии.
Введем безразмерные параметры:
р = р,Ю0, ст=-^-, qv = Еи3, qw = ,
®0 Еи у Еиу
где р. и ст. - коэффициенты в модели внутреннего и внешнего трения.
Последнее уравнение системы (1) получено из условия отсутствия относительного продольного смещения опор. Кручение не учитывается, поскольку в дальнейшем рассматриваются преимущественно изгибные колебания тонкого стержня в окрестности первой собственной частоты.
Система (1) должна быть дополнена граничными условиями, в частности, для стержня с двумя заделками: дv дw
и = V = w = — = — = 0 при х = 0,1.
дх дх (2>
Метод решения. Ограничимся случаем одномодового приближения и представим решение в виде
V(X,t) = 9oi(f)V(х), w(x,t) = 902(f)V(x),
где V(x) - первая форма колебаний стержня:
V (х) = Vo[(sin а - sha) (hax - cos ах) +
+ (cha - cos a) (hax - sin ax)],
a = 4,730, V0 =-0,0109.
Форма колебаний нормирована из условия, что прогиб в средней точке стержня равен единице. Подстановка в систему (1) и применение процедуры Бубнова-Галеркина приводит к системе с двумя степенями свободы:
ср 01 + п0ф01 +Ф01 + т(ф01 + ф02 )ф01 =
= f01cos(|i t + 01),
4502 + (п0 + в0(с - 1))ф02 + Сф02 + (3)
+Y (фС)1 + ф02) ф02 = f02 cos(^t + 02),
где у = 0,240у0 для стержня с двумя заделками,
П0 = Р0 + ст0 .
В уравнениях (3) функции p01(t), p02(t) описывают перемещения средней точки стержня в
двух плоскостях, отнесенные к длине стержня L. Выберем в качестве нового масштаба единицы длины амплитуду отклонений A, A<<L. Произведем замену переменных:
A A A A
ф01 = ф1 L ’ ф02 = ф2 L ’ f01 = f11 L ’ f02 = f12 L ! ^
В новых переменных уравнения (3) имеют следующий вид:
ср1 + snp1 + Ф1 + еу (2 + ф° |ф1 =
= ef1 cos(|it + 01),
ф2 + (ел + е2р8)ср 2 +(1 + е§)ф2 + (4)
/ 2 2 \ ()
+еу I ф2 + Ф2 IФ2 = ef2 cos(^t + 02 ),
( A f
где введен малый параметр е = I — I и следующие обозначения:
с =1 + ^ еп = ^ еР = Р0, ef1 = f11 = f01 "т = fr?,
A е
ef = f = f L = f02
2 = f12 = f02 A = е1/2 .
Амплитуды внешних нагрузок, параметры диссипации и разность собственных частот предполагаются малыми порядка е, что позволяет применить эффективные методы нелинейной механики. Отметим, что для стержней L >> r , где r -радиус инерции поперечного сечения r2 = Jz / F , и малость еу может быть достигнута выбором величины A. Решение системы уравнений (3) совпадает с решением уравнений (4) при значении е = 1 или A = L . Далее рассматриваются колебания в малой окрестности главного резонанса, т.е. при
ц = 1 + е^, е^ << 1.
Представим систему (4) в стандартной форме [4] при помощи перехода к медленным переменным:
(5)
Фк = ак соэ(цt + а к ), фк =-цакэ1п(цt + ак), к = 1,2, где а1, а2 - амплитуды парциальных колебаний; а1,а2 - фазовые добавки. При е = 0 эти величины постоянны и ц = 1. Для значений е> 0 медленные переменные изменяются со скоростью е. Произведя подстановку (5) в уравнения (4) и применяя метод усреднения, получим систему уравнений в медленных переменных:
1
ак = 2е^к Э1п(ак -0к) +
1 к 3—кж^/ \ 1
+ 8а1 а2 5|п2(ак-а3-к) -2епак,
ак = ^ е [8(к - 1) - 2^] -1 { I (\ \ 3 2
-----е ?к соэ(ак + 0к) +— еу ак +
(6)
2а,
8
1
+—еуа^к(2 + соэ2(а1 -а2)), к = 1,2.
8
Слагаемые 0(е2) в уравнениях (6) упущены. Далее будем рассматривать установившиеся вынужденные колебания стержня, что соответствует нулевым левым частям уравнений (6). Для численного решения системы нелинейных алгебраических уравнений воспользуемся методом продолжения решения по параметру. Если при некотором значении Хк известно приближенное решение системы уравнений г = (ак,а%,ак,а£)т , то для значения Хк+1 =Хк + А^к приближенное решение
представляется в виде гк+1 = гк +Агк. Подставляя в систему (6) и линеаризуя полученные уравнения, определяем приращения неизвестных из системы
(7)
Gk Ark = pk + Rk,
где Як - вектор невязки на предыдущем шаге решения; рк = (0,0Дк,Хк)т .
Знак действительной части всех собственных значений матрицы О позволяет сделать вывод об устойчивости решения. Поскольку матрица несимметрична, для решения проблемы собственных значений использовался ОР-алгоритм.
Результаты численного моделирования. В качестве примера рассмотрим задачу о вынужденных колебаниях латунной трубы длиной 1м, жестко закрепленной по краям. Внешний диаметр трубы с1 = 16 мм, внутренний - С0 = 14 мм. Радиус инерции се-
с Г. СТ
чения г = — = —, 1+ -0 = 5,3 мм, параметр
Р 4 V С2
нелинейности у = 0,24 J = 2125.. Параметр диссипации в системе уравнений (3) при-
нимался равным П0 = 5 • 10-3, параметр неравноже-сткости с = 1,012. Значение параметра диссипации определено по результатам экспериментальных исследований свободных колебаний. Возбуждение колебаний производилось электромагнитным вибровозбудителем, регистрация колебаний - датчиком виброускорений.
На рис. 1, 2 приведены амплитудно-
частотные характеристики а2(ц),а1(ц) для четырех значений амплитуды обобщенной силы /02, входящей в систему уравнений (3): /02 = 3,2*10-6 (кривая 1); /0)2 = 8-10"6 (кривая 2); /02 = 2,08*10-5 (кривые 3); /02 = 3,2*10-5 (кривые 4). Устойчивые решения показаны жирными линиями. Кривая 1 соответствует малым значениям амплитуды возбуждения, колебания трубы в двух плоскостях происходят независимо друг от друга. АЧХ практически совпадают с результатами, полученными по линейной теории. С увеличением амплитуды возбуждения (кривая 2) получаем две независимые нелинейные амплитудночастотные характеристики а1(ц) и а2(ц). С увеличением амплитуды возбуждения колебаний (кривые 3) происходит качественное изменение АЧХ. В области резонанса взаимодействие между формами колебаний приводит к существованию пространственных режимов движения стержня. На амплитудночастотных характеристиках четвертого случая нагружения (кривые 4) присутствуют все возможные режимы колебаний в области резонанса: два пространственных, с вращением в противоположные стороны и почти плоское движение в плоскости большей изгибной жесткости стержня. Для неравножесткого стержня при определенных значениях параметров задачи одновременно существуют несколько режимов колебаний.
Две пространственные формы движения отличаются направлением вращения сечений, причем одна из них реализуется на изолированном участке АЧХ. Максимальные амплитуды реализуются на колебаниях, имеющих форму эллипса, сильно вытянутого вдоль оси т, т. е. при дви-
жении практически в плоскости большей из-
гибной жесткости сечения стержня.______________
Рис. 1. АЧХ стержня в плоскости возбуждения колебаний
Рис. 2. АЧХ стержня в плоскости, ортогональной возбуждению колебаний
Список литературы
1. Муницын А.И. Пространственные нелинейные колебания стержня с неподвижными шарнирными опорами // Прикладная математика и механика. - 2006. - Т. 70. -Вып. 1. - С. 72-80.
2. Муницын А.И. Нелинейные колебания стержня с близкими значениями осевых моментов инерции поперечного сечения // Прикладная математика и механика. - 2009. - Т.73. - Вып. 3. - С. 427-438.
3. Каудерер Г. Нелинейная механика. - М.: Изд-во иностр. лит., 1961.
4. Найфэ А.Х. Введение в методы возмущений. -М.: Мир, 1984.
Муницын Александр Иванович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры теоретической и прикладной механики, телефон 8-910-667-11-48, e-mail: munitsyn@rambler.ru
Крайнова Лариса Николаевна,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», старший преподаватель кафедры теоретической и прикладной механики, телефон 8-910-981-68-57.
Сабанеев Николай Андреевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», старший преподаватель кафедры теоретической и прикладной механики, телефон 8-910-981-68-57.
