CC BY
33
АРХИТЕКТУРА И СТРОИТЕЛЬСТВО
УДК 624.04
НЕКОТОРЫЕ ВОПРОСЫ СВОБОДНЫХ НЕЛИНЕЙНЫХ КОЛЕБАНИЙ
БАЛОЧНЫХ КОНСТРУКЦИЙ
Докт. техн. наук, проф. БОСАКОВ С. В., инж. ЩЕТЬКО Н. С.
Белорусский национальный технический университет, РУП «Институт БелНИИС»
В прошлом веке были досконально исследованы многие вопросы линейной динамики элементов строительных конструкций, в том числе и задачи линейных колебаний стержней и балок [1-4]. Однако эти исследования, несмотря на богатство в содержании, зачастую весьма ограничены в возможностях качественного изучения нелинейных систем.
Строительные конструкции являются во многом нелинейными системами. Это обстоятельство указывает на необходимость выполнения нелинейного анализа, отправной точкой для которого могут быть решения, полученные в рамках линейной теории.
Нелинейным системам свойственны некоторые общие отличительные особенности, такие как: неизохронность свободных колебаний, неоднозначность зависимости амплитуды вынужденных колебаний от частоты гармонической вынуждающей силы, возможность существования автоколебательных режимов, реализация которых зависит от начальных условий [5-7]. Эти явления имеют место при практических расчетах, а также в работе реальных конструкций и не могут быть объяснены в рамках линейной теории.
Вовлечение нелинейностей в модели конструкций позволяет более точно проанализировать систему, запроектировать конструкцию либо оценить работу существующих конструкций и систем.
Традиционно принято ограничивать учет нелинейностей в динамике элементов строительных конструкций рассмотрением так называемой геометрической нелинейности протяженных конструкций. На это ориентированы
и программные комплексы, используемые для расчетов зданий и сооружений [8-11]. В отечественной и зарубежной литературе вопросы учета физической нелинейности строительных конструкций, подверженных динамическим воздействиям, не освещены в должной мере.
Задачи физически нелинейной динамики конструкций, рассмотренные в других областях науки [12, 13], не могут быть полностью приемлемы для расчетов строительных конструкций, но могут быть полезными при выборе направления и методов исследований и позволить уточнить известные решения с учетом нелинейных свойств, присущих строительным конструкциям. Отмеченное направление, по мнению авторов, может представлять определенный теоретический и практический интерес.
В качестве объектов исследования выбраны однопролетные балочные конструкции. Материал - один из основных конструкционных традиционно применяющихся в строительстве: металл, железобетон, древесина.
Балка является одним из основных конструктивных элементов, способным обеспечить адекватное моделирование многих объектов, таких как колонны, ригели, здания, мосты.
Важной задачей прикладной теории нелинейных колебаний является задача определения начальных условий, при которых система будет совершать устойчивые колебания. В случае превышения найденных величин эти же условия будут являться условиями неограниченного роста (срыва) амплитуд [14]. Цель настоящих исследований - определение особенностей нелинейных колебаний, обусловленных физической нелинейностью материала конструкций,
а именно: начальной энергии системы, необходимой для существования устойчивых свободных нелинейных колебаний балок, зависимости периода от начальной энергии однопролетных балок из металла, железобетона, древесины при характерных условиях опирания, а также определение степени влияния нелинейности на характеристики колебаний.
Для описания работы элементов изгибаемых конструкций актуально применение уравнений, содержащих мягкую упругую силу, что соответствует известному в практике понижению жесткости конструкции с ростом деформаций.
При решении уравнений с мягкой нелинейностью, описывающих колебания балок [15-18], авторами отмечена чувствительность решения к начальным условиям, что выражалось в резком неограниченном росте амплитуд в определенных интервалах изменения начальных условий даже при незначительных вариациях последних. Эти явления объясняются наличием нелинейных членов в разрешающем уравнении и отмечены рядом исследователей вопросов теории нелинейных колебаний [6, 19].
Устойчивость нелинейных колебаний однопролетных балок. При анализе поведения нелинейных колебательных систем удобно воспользоваться фазовыми диаграммами (плоскостями), объединяющими совокупность фазовых траекторий, отражающих состояние автономной системы, определяемой обобщенными координатами и обобщенными скоростями, т. е. /(у, у). Оперируя понятием фазовой плоскости, можно получить важные заключения относительно качественного характера движения [20].
Как уже отмечено в [15, 16], задачи нелинейной (с учетом физической нелинейности материала) динамики строительных конструкций с одной степенью свободы сводятся к решению уравнения Дюффинга с мягкой нелинейностью. Данное уравнение более точно описывает процессы деформирования, происходящие в системе, чем используемые в настоящее время линейные модели [1, 2].
Уравнение имеет следующий вид [18]:
у + ю2у - а2у3 = 0,
(2)
Му + гу -vy3 = 0. Либо после преобразования
где
Юо =
а =
Интеграл энергии уравнения (1) имеет вид
[21]
1 2
-у2 + и = С, 2
(3)
где и - потенциальная энергия; С - полная энергия системы.
В (1) гу -vy3 = Я(у) - восстанавливающая упругая сила. Ее работа определяет количество потенциальной энергии системы
и =
2 4
_гУ , уу
(4)
2 4
Тогда из выражения (3) определяем у как у = 2(С - и). (5)
Схемы графиков изменения и от у и у от у представлены на рис. 1. Отметим, что график зависимости у от у подобен фазовому портрету и отражает фазовые траектории, соответствующие предельно допустимым значениям начальных условий.
Сепаратриса
Область
неустойчивых у колебаний
Рис. 1. Схема зависимости энергии системы от начального отклонения и эллипс «критических» значений начальных условий, вызывающих нелинейные колебания одиночной массы с нарастающей амплитудой
Сепаратриса отделяет область начальных условий, которым соответствуют периодиче-
ские движения, от области начальных условий, которым соответствуют неограниченные движения.
Уравнение эллипса экстремальных значений начальных условий в нашем случае имеет вид
+=1
■2 2
Уо Уо
(6)
Критическое значение начальная потенциальная энергия приобретает при начальном отклонении
Уо =
(7)
Критическое значение начальной скорости
консервативной системы уо =
М V
находим
из условия
му2_
итах.
(8)
Точка с координатами I ./—, — I на графиН V 2v )
ке зависимости и от у получила название вершины энергетического барьера.
Для иллюстрации приведенных рассуждений рассмотрим шестиметровую однопролет-ную балку из основных конструктивных материалов - металла, железобетона, древесины.
Закон деформирования материала балки представим в виде
4 Е3
о(е ) = Ее---- е3,
27 о2
(9)
где с - предел прочности; Е - начальный модуль упругости материала балки.
Размеры сечений и характеристики материалов:
1. Металлическая двутавровая балка 30Б2: с Е = 2 • 1011 Па; с = 4 • 108 Па.
2. Железобетонная балка прямоугольного поперечного сечения.
Примем допущение, что рабочее армирование балки симметричное, т. е. жесткость в противоположных направлениях вертикальной плоскости одинакова. Это позволит избежать
расчетных сложностей, вызванных особенностями нелинейных колебаний железобетонных изгибаемых конструкций (часть из которых обозначена в (4)), проследить некоторые общие закономерности и сделать соответствующие выводы. Для железобетонной балки приняты следующие характеристики: Е = 3,5 • 1010 Па; с = 2 • 107 Па; Ь = 0,5 м; И = 0,22 м.
3. Балка из деревянного бруса квадратного поперечного сечения: Е = 1,1 • 1010 Па; с = = 9 • 106 Па; Ь = 0,36 м; И = 0,36 м.
В качестве начальных условий примем отклонение у0 = 0,02 м и скорость у0 = 0,1 м/с. Для всех балок распределенная масса приведена к сосредоточенной, обозначенной как М.
Размеры сечений обеспечивают равное значение начальной изгибной жесткости каждой из трех балок на уровне В = 1,54107 Н •м .
Проанализируем три характерных типа опи-рания балки (рис. 2).
у ^ 0
1/2
1/2
У !<} 0?
М1
/
_ I _
в
У ^ 0^М
I
х
г"
Рис. 2. Три типа опирания однопролетной балки с одной степенью свободы
Уравнения движения составляются с использованием принципа Даламбера. Задачи решаются энергетическими методами. После нахождения потенциальной энергии изгиба балки на основании теоремы Кастилиано [2] определяется нелинейная сила упругости. Ход решения подробно освещен в [15-18, 22].
2
Г
а
2
б
X
В качестве координатных функций используются фундаментальные балочные функции [18, 23], удовлетворяющие всем граничным условиям. Предварительно ограничимся рассмотрением изгибной формы колебаний.
Графическое отображение полученных решений представлено на рис. 3.
Полученные значения параметров г и V сведены в табл. 1.
Из представленных данных (рис. 3, табл. 1) следует, что при повышении общей вертикальной жесткости балок, обусловленной условиями закрепления (при переходе от свободного опирания к защемлению), прослеживается тенденция снижения значений начальных «критических» условий. Это соответствует меньшему
количеству потенциальной энергии и скорейшему преодолению энергетического барьера, что легко проследить по графикам зависимости и от у.
Для металлической балки область предельно допустимых значений начальных условий находится за пределами допустимых по жесткости и условиям нормальной эксплуатации [24, 25]. Высота энергетического барьера металлической балки значительно превышает величину, полученную для железобетонной балки и балки из деревянного бруса, что можно объяснить большим значением погонной жесткости двутавровой балки по сравнению с двумя другими.
а
и
б и
/ 1
/2
.1 -0. 05 0. 05 ' 0
з \
/1 -15 / 2 3 \
-0 0 1
--15-
Рис. 3. Графики зависимостей потенциальной энергии однопролетной балки от начального отклонения и эллипс критических значений начальных условий, вызывающих нелинейные колебания одиночной массы с нарастающей амплитудой: а - металлическая балка; б - железобетонная балка; в - балка из деревянного бруса; 1 - показано решение для свободно опертой балки; 2 - то же для опертой с одной и защемленной с другой стороны; 3 - для защемленной балки
в
Таблица 1
Параметры нелинейных колебаний балок с одной степенью свободы
Тип опирания (рис. 2) а б в
Металлическая балка. Двутавровое поперечное сечение 30Б2
г, 1/с2 3,32Е + 06 7,77Е + 06 1,35Е + 07
V, 1/(м2-с2) 1,43Е + 08 8,80Е + 08 2,93Е + 09
Железобетонная балка. Прямоугольное поперечное сечение 0,22x0,5 м
г, 1/с2 3,50Е + 06 8,19Е + 06 1,43Е + 07
V, 1/(м2-с2) 6,50Е + 08 4,00Е + 09 1,33Е + 10
Балка из деревянного бруса. Квадратное поперечное сечение 0,36x0,36 м
г, 1/с2 3,47Е + 06 8,12Е + 06 1,41Е + 07
V, 1/(м2-с2) 8,42Е + 08 5,18Е + 09 1,73Е + 10
Для железобетонной балки и балки из деревянного бруса значения погонной жесткости отличаются не более чем на 15 %, что подтверждается графиками бив рис. 3. Завышенное (на 45 %) значение «критической» скорости балки из деревянного бруса по сравнению с железобетонной балкой является следствием различий в массовых характеристиках материалов.
В численных экспериментах обнаружено также (рис. 4), что с ростом начальных условий, а, следовательно, и напряжений в конструкции возрастает вклад нелинейной составляющей (9) в разрешающем уравнении (1) в колебания балок. Это, в свою очередь, приво-
дит к увеличению погрешности линейного решения.
На основании полученных зависимостей проанализируем влияние количества начальной энергии системы, определяемое начальными условиями, на период свободных нелинейных колебаний однопролетных балок.
Зависимость периода свободных колебаний от начальной энергии системы. Выражение начальной энергии уравнения Дюффинга (1) в безразмерном виде имеет следующий вид:
е = ^л/уМ . (10)
г
А
0,06 0,04 0,02
-0,02 -0,04 -0,06
4
( У \ / Л
\\ 1 / \
\\ / / к \
0. 01 32 0. 33 0. 34 0. 35 3.
\\ i > ш V \
V * 1 / \ / \
\
36
0,10 0,05
-0,05 -0,10
у > \ \ 4 /
\ \ \ 1 ■ Л ■ / / А \
0. 31 0. \ 32 0. 3? 0. г 34 0. 35 0.0
1 \ У \ ч
б
а
I
Рис. 4. Свободные колебания (металлической свободно опертой) балки: сплошной линией показаны нелинейные колебания, пунктирной - линейные; а - у0 = 0,044 м, v0 = 10 м/с; б - 0,04 м, 15 м/с
По данным [5], период собственных колебаний уравнения (1) определяется выражением
Т =
4*Л
К
- 2в1
еу/2
1 + Л/Т-2ё2
(11)
где К(т) - полный эллиптический интеграл [26].
Формула (11) отражает такое свойство нелинейных колебательных систем, как неизохронность.
Согласно (11) переход к апериодическому движению произойдет, когда Т = да или
ёл/2 / \
= 1, так как К (1) = 1п да .
1 + 4\-2в£
После ряда преобразований находим корни аргумента полного эллиптического интеграла К(т):
>( 4е - 2^2 ) = 0; е1 = 0; е2 =
42 2
Тогда
г 2
(12)
(13)
На рис. 5 представлен график зависимости безразмерного периода колебаний металлической свободно опертой балки от безразмерной энергии. Нулевому значению энергии соответствует величина периода 2п. Максимальное
42
значение энергии--— .
Рис. 5. График зависимости Т - е для собственных колебаний металлической свободно опертой балки
Выражение (13) описывает движение по сепаратрисе. Отметим, что (11) не очень удобная формула, так как учитывает только влияние у0 -
начальной скорости отклонения, однако может быть полезна при решении определенного класса задач.
Известно также следующее решение для периода свободных колебаний [27]:
_у(С ) ,
Т (С) = ш Г -7Jy^=.
- у{с )4С-Щу)
(14)
Зависимость (14) определена через значения начальных отклонений, но не учитывает скорость.
В [18] авторами получены точные решения уравнения (1) для периода
Т =
4л/2К (т)
д/®2 +д/®<Т- 2а2 С
(15)
и амплитуды
А =
1
2СХ
+у]®4 - 2а2С
(16)
ю2 - ^/ю4 - 2а2С где т =-\ в которых учитыва-
ют + ^ю4 - 2а 2С1
ются значения обоих параметров, определяющих начальные условия и благодаря чему, по мнению авторов, их можно признать более корректными по сравнению с известными. Решения выражены через эллиптические функции. В выражениях (15), (16) введен коэффициент Сь который определяется с учетом коэффициента С2 при заданных начальных условиях решением системы уравнений
2С
®4 - 2а2С
-£п(и(0), т) = 20;
1
4С\Сп(и(0), т)0п(и(0), т) = у0.
(17)
Коэффициент С2 представлен в качестве множителя при величине и эллиптических функций Якоби [26] Сп(и, т), Бп(и, т) системы (17).
Начальные приближения для коэффициентов определяются как:
С1 = 1 ® -1 а22021;
С =-
a42K (m)
4 д/- 2а2Q
(18)
и уточняются решением системы (17).
Для графиков, представленных на рис. 4, по формулам (15)—(18) определены следующие значения:
а) С = 146,726351 м2/с2; С2 = 0,003511 с; T = 0,039876 с; А = 0,074178 м;
б) С = 271,726351 м2/с2; С2 = 0,002503 с; T = 0,046769 с; А = 0,110312 м.
Значение периода линейных колебаний для анализируемой системы составляет T = = 0,036129 с.
Погрешность вычисления относительно линейного решения для периода составляет от 9 % (для варианта а) до 23 % (для варианта б).
Все исследования проведены с использованием средств символьной математики и компьютерной алгебры Mathematica, Maple [20, 29], благодаря чему получены точные решения обозначенных выше задач.
Вклад нелинейности при рассмотрении свободных колебаний балок как консервативных систем повышается по мере увеличения значений начальных условий, что соответствует большим значениям напряжений в конструкции и, как следствие, росту погрешности линейного решения.
Факт снижения жесткости по мере роста напряжений в материале конструкции является известным в механике (статике) строительных конструкций, и полученные закономерности для колебаний балок имеют логическое объяснение.
Деградация жесткости конструкции может являться следствием не только нелинейности поведения материала под нагрузкой, но также деструктивных процессов, свойственных каждому из материалов: коррозии металла, гниения древесины, трещинообразования в растянутой зоне и дробления сжатой зоны железобетона.
На обозначенных свойствах нелинейных колебаний основаны получившие в последнее время широкое распространение методы диагностики конструкций [30], базирующиеся на анализе изменения колебаний конструкций.
Полученные авторами данные могут быть использованы и учтены при проектировании новых конструкций, а также для оценки динамических параметров и технического состояния существующих строительных конструкций зданий и сооружений.
Приведенные в настоящей работе решения справедливы для консервативных систем. Для более подробных практических расчетов элементов строительных конструкций необходимо учесть диссипативные свойства материалов, для чего полезно воспользоваться результатами, приведенными в [16].
В Ы В О Д Ы
1. Построены графики начальных условий, определяющих область устойчивых нелинейных колебаний балок с одной степенью свободы с различными опорными закреплениями для характерных конструкционных материалов.
2. Определены значение количества начальной энергии, необходимого для существования устойчивых колебаний, а также зависимость безразмерного периода от величины безразмерной энергии балки.
3. Получены точные значения амплитуды и периода нелинейных колебаний с учетом двух параметров начальных условий по предложенным авторами формулам.
4. Получены новые сведения о качественном характере колебаний элементов строительных конструкций, и представлена количественная оценка влияния нелинейности на колебательный процесс.
Л И Т Е Р А Т У Р А
1. Киселев, В. А. Строительная механика: спец. курс «Динамика и устойчивость сооружений» / В. А. Киселев. - М.: Мир, 1980. - 548 с.
2. Ржаницын, А. Р. Строительная механика / А. Р. Ржа-ницын. - М.: Высш. шк., 1991. - 438 с.
3. Справочник по динамике сооружений / под ред. Б. Г. Коренева, И. М. Рабиновича. - М.: Стройиздат, 1972. - 511 с.
4. Тимошенко, С. П. Колебания в инженерном деле / С. П. Тимошенко, Д. Х. Янг, У. Уивер. - М.: Машиностроение, 1985. - 472 с.
5. Моисеев, Н. Н. Асимптотические методы нелинейной механики / Н. Н. Моисеев. - 2-е изд., перераб. - М.: Наука, 1981. - 400 с.
6. Найфэ, А. Методы возмущений / А. Найфэ. - М.: Мир, 1976. - 455 с.
7. Пановко, Я. Г. Введение в теорию механических колебаний / Я. Г. Пановко. - М.: Наука, 1991. - 255 с.
8. Вычислительный комплекс SCAD / В. С. Карпи-ловский [и др.]. - М.: АСВ, 2004. - 529 с.
9. Программный комплекс для расчета строительных конструкций на прочность, устойчивость и колебания STARK ES. Версия 4.2 (2006). Руководство пользователя. - М.: Еврософт, 2006. - 383 с.
10. Симбиркин, В. Н. Статический и динамический расчет железобетонных монолитных каркасов зданий с помощью программного комплекса STARK ES: учеб. пособие / В. Н. Симбиркин, С. О. Курнавина; под ред. Ю. П. Назарова. - М.: ФГУП «НИЦ "Строительство"», ООО «Еврософт», 2007. - 158 с.
11. Moaveni, S. Finite Element Analysis: Theory and Application With ANSYS Prentice-Hall, Upper Saddle River, NJ, 1999. - Mode of access: http://www.proquest.com/ - Date of access: 2.08.2007.
12. Kuiyin, Mei. Nonlinear vibration of beam structures: analytical and explicit dynamics FEA evaluations // А dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the University of Kentucky, 2000. - Mode of access: http://www.proquest.com/ -Date of access: 15.09.2007.
13. Mahmood, M. Tabaddor. Nonlinear vibration of beam and multibeam systems // Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Engineering Mechanics, Sept. 24, 1996. - Mode of access: http://www.proquest.com/ - Date of access: 4.12.2007.
14. Вибрации в технике: справ. в 6 т./ ред. совет: В. Н. Челомей (пред.). - М.: Машиностроение, 1979 - Т. 2: Колебания нелинейных механических систем / под ред. И. И. Блехмана. - М., 1979. - 351 с.
15. Босаков, С. В. К теории свободных колебаний балок из физически нелинейного материала / С. В. Босаков, Н. С. Щетько // Вестник БНТУ. - 2006. - № 1. - С. 10-14.
16. Босаков, С. В. Исследование свободных затухающих колебаний балок из физически нелинейного материала // С. В. Босаков, Н. С. Щетько // Вестник БГТУ. -2006. - № 1. - С. 134-138.
17. Босаков, С. В. Об одном свойстве зависимости «момент - кривизна» для балок и его использовании в инженерных расчетах // С. В. Босаков, Н. С. Щетько // Строительная наука и техника. - 2006. - № 1. - С. 5 8-61.
18. Босаков, С. В. К теории нелинейных колебаний балок с одной степенью свободы / С. В. Босаков, Н. С. Щетько // Вестник БНТУ. - 2007. - № 3. - С. 5-9.
19. Бидерман, В. Л. Прикладная теория механических колебаний / В. Л. Бидерман. - М.: Высш. шк., 1972. -416 с.
20. David, Charles Hill. Time delayed dynamical systems and the Duffing Equation. A dissertation submitted to the faculty of the university of North California at Charlotte in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Applied Mathematics. Charlotte 2005. -Mode of access: http://www.proquest.com/ - Date of access: 6.02.2008.
21. Камке, Э. Справочник по обыкновенным дифференциальным уравнениям / Э. Камке. - 6-е изд. - М.: Наука, 1976. - 576 с.
22. Босаков, С. В. Метод Ритца в примерах и задачах по строительной механике и теории упругости / С. В. Босаков. - Минск: БГПА, 2000. - 144 с.
23. Фаддеева, В. Н О фундаментальных функциях оператора xIV / В. Н. Фаддеева // Труды математического института имени В. А. Стеклова. - 1949. - Т. XXVIII. -С. 157-159.
24. Нагрузки и воздействия: СНиП 2.01.07-85* / Минстрой России. - М.: ГП ЦПП, 1996. - 44 с.
25. ASCE STANDART/SEI 7 - 05. Minimum Design Loads for Buildings and Other Structures. USA, 2006. - Mode of access: http://www.proquest.com/ - Date of access: 26.11.2007.
26. Ленг, С. Эллиптические функции / С. Ленг; пер. с англ. - М.: Наука, 1984. - 312 c.
27. Ландау, Л. Д. Серия «Теоретическая физика». Т. 1 / Л. Д. Ландау, Е. М. Лившиц. - М., 1973. - 208 с.
28. Гандер, В. Решение задач в научных вычислениях с применением Maple и MATLAB / В. Гандер, И. Гржеби-чек. - М.: Вассамедиа, 2005. - 520 с.
29. Кулешов, А. А. Уравнения математической физики в системе Mathematica / А. А Кулешов. - Минск: БГУ, 2004. - 294 c.
30. Maeck, J. Damage Assessment of Civil Engineering Structures by Vibration Monitoring. Ph. D. Dissertation, Department of Civil Engineering, K. U. Leuven, Belgium, 2003. - Mode of access: http://www.proquest.com/ - Date of access: 19.06.2007.
Поступила 3.03.2008
