Математические модели сложных колебаний балок в условиях знакопеременных и ударных нагрузок
На правах рукописи
Салтыкова Ольга Александровна
МАТЕМАТИЧЕСКИЕ МОДЕЛИ СЛОЖНЫХ КОЛЕБАНИЙ БАЛОК
В УСЛОВИЯХ ЗНАКОПЕРЕМЕННЫХ И УДАРНЫХ НАГРУЗОК
Специальности: 05.13.18 – Математическое моделирование, численные
методы и комплексы программ
01.02.04 – Механика деформируемого твердого тела
А в т о р е ф е р а т
диссертации на соискание ученой степени
кандидата физико-математических наук
Саратов 2008
Работа выполнена в Государственном общеобразовательном учреждении высшего профессионального образования «Саратовский государственный технический университет» (Россия) и Техническом Университете г. Лодзь (Польша)
Научные руководители – доктор технических наук, профессор
Крысько Вадим Анатольевич (Россия)
– доктор, профессор
Аврейцевич Ян (Польша)
Официальные оппоненты –доктор физико-математических наук, профессор Серазутдинов Мурат Нуриевич
(Казанский государственный технологический университет)
– доктор физико-математических наук, профессор
Землянухин Александр Исаевич
(Саратовский государственный университет
им. Н.Г. Чернышевского)
Ведущая организация - Институт проблем точной механики и управления РАН, г. Саратов
Защита состоится «26» ноября 2008г. в 13.00. на заседании диссертационного совета Д 212.242.08 при ГОУ ВПО «Саратовский государственный технический университет» по адресу: 410054, Саратов, ул. Политехническая, 77, Саратовский государственный технический университет.
С диссертацией можно ознакомиться в научно-технической библиотеке ГОУ ВПО «Саратовский государственный технический университет».
Автореферат разослан « 23 » октября 2008 г.
Ученый секретарь
диссертационного совета А.А. Терентьев
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. В связи с появлением сложных инженерных и технических сооружений в последние десятилетия к вопросам динамики конструкций проявляется особый интерес. Активно расширяется спектр применения балочных конструкций. Так, развитие авиационной, строительной и морской техники выдвинуло в число наиболее актуальных задач изучение нелинейного, в том числе хаотического поведения балок, их динамики и устойчивости при воздействии внешних нагрузок. Резко возрастают требования к оценкам прочности и экономичности таких конструкций. Отдельным вопросом изучения динамики механических систем является вопрос о переходе систем в состояние хаоса под действием различного вида нагрузок.
Значительный вклад в изучение этого вопроса внесли такие ученые, как П.С.Ландау, Е.А.Хопф, М.Фейгенбаум, Н.Помо, Ю.И.Неймарк, В.А.Крысько, Я.Аврейцевич, П.А.Ланда, П.Манневиль, Д.И.Трубецков, У.Лепик и др.
Исследованиям хаотических колебаний пологих, замкнутых цилиндрических оболочек, круглых и прямоугольных пластинок, бесконечно длинных панелей, а также балок, посвящены работы В.А.Крысько, Я.Аврейцевича, Ю.Г.Коноплева, А.В.Крысько, Т.В.Вахлаевой, Т.В.Щекатуровой, И.В.Папковой, Н.Е.Савельевой, Э.С.Кузнецовой, М.В.Жигалова, О.Н.Киреевой, Г.Г.Наркайтиса. Исследования поведения пластин и оболочек под действием ударных нагрузок можно найти в работах В.Г.Баженова, В.А.Крысько, А.М.Варыгина и др. Однако в работах этих авторов не достаточно изучены сложные колебания балок в условиях знакопеременных и ударных нагрузок. В связи с этим важной и актуальной является задача построения детерминированных математических моделей, позволяющих исследовать хаотические колебания балок при воздействии поперечной знакопеременной нагрузки и продольного удара груза.
Целью работы является построение математических моделей нелинейных колебаний сложных механических систем в виде балок с учетом гипотез Эйлера-Бернулли, С.П. Тимошенко, Пелеха-Шереметьева. Для достижения этой цели необходимо решить следующие задачи:
- Разработка математических моделей для сложных колебаний гибких балок по гипотезам Эйлера-Бернулли, С.П.Тимошенко, Пелеха-Шереметьева для некоторых типов граничных условий под действием знакопеременной и ударной нагрузок.
- Изучение сценариев перехода от гармонических колебаний к хаотическим для различных гипотез с учетом некоторых управляющих параметров.
- Разработка алгоритма и комплекса программ на ПЭВМ для качественного исследования сложных колебаний диссипативных систем в виде упругих балок с учетом различных гипотез при произвольных граничных условиях.
- Качественное исследование динамики гибких балок на основе нелинейной динамики в зависимости от изменения следующих параметров: краевых условий, амплитуды и частоты равномерно распределенной поперечной и продольной знакопеременных нагрузок, ударных нагрузок, величины диссипативных членов, угла поворота и искривления нормали.
Направление исследований данной работы является изучение динамики сложных механических систем в виде балок, изучение сценариев перехода к хаосу в этих системах, выявление точности используемой математической модели.
Основными методами исследований являются методы качественной теории дифференциальных уравнений; методы конечных разностей с аппроксимацией 
, конечных элементов с аппроксимацией в форме Бубнова-Галеркина по пространственным координатам, позволяющие сводить уравнения в частных производных к системам обыкновенных дифференциальных уравнений. Исследуется сходимость этих методов в зависимости от шага разбиения по пространственной координате; методы анализа показателей Ляпунова; методы типа Рунге-Кутта четвертого порядка точности, позволяющие решать обширные системы обыкновенных дифференциальных уравнений по времени. Для указанных методов, исследована устойчивость решения в зависимости от соотношения шагов по временной и пространственной координатам.
Достоверность и обоснованность обеспечивается корректной физической и математической постановкой задачи, применением методов математического моделирования, качественной теории дифференциальных уравнений и нелинейной динамики, сравнением результатов, полученных методом конечных разностей и методом конечных элементов по пространственной координате.
Результаты, полученные автором диссертации, согласуются с имеющимися физическими представлениями, основанными на экспериментах.
Научная новизна работы заключается в следующем:
- Построена математическая модель сложных колебаний гибких упругих балок с некоторыми краевыми условиями под действием поперечной знакопеременной и ударной нагрузок, с использованием гипотез Эйлера-Бернулли, позволяющая проводить качественный анализ динамического поведения балки при действии различных управляющих параметров (тип трения и значения коэффициентов диссипации энергии, учет различных граничных условий, относительная толщина балки, частота и амплитуда вынуждающих колебаний, скорость груза в момент удара, отношение массы балки к массе груза и др.).
- Построенная математическая модель сложных колебаний гибких упругих балок с произвольными краевыми условиями под действием поперечной знакопеременной нагрузки, с использованием гипотез С.П. Тимошенко учитывает инерционные составляющие, связанные с поворотом сечений. Реализация этой математической модели позволила провести качественный анализ динамического поведения балки при действии различных управляющих параметров (граничные условия, амплитуда и частот вынуждающих колебаний, величина геометрического параметра).
- Построенная математическая модель сложных колебаний гибких упругих балок с произвольными краевыми условиями под действием поперечной знакопеременной нагрузки с использованием гипотез Пелеха-Шереметьева, учитывающая инерционные составляющие, связанные с поворотом сечений, позволяет изучать влияние граничных условий, амплитуды и частоты вынуждающих колебаний, а также геометрического параметра на динамическое состояние балки.
- Сценарии перехода колебаний гибких балок от гармонических в хаотические, для различных граничных условий, с учетом различных управляющих параметров. Изучение сценариев позволяет провести их классификацию в соответствии с ранее изученными, а также прогнозировать поведение исследуемой динамической системы во времени при изменении управляющих параметров.
- Учет различных типов трения (кулоновское, нелинейное, линейное) и упругих оснований Винклера и В.З.Власова для гибкой и жесткой балки Эйлера-Бернулли позволило провести сопоставление влияния учета упругих оснований Винклера и В.З.Власова и типов трения, для математической модели балки Эйлера-Бернулли.
- Разработаны и реализованы алгоритмы, методика и комплекс программ анализа хаотических колебаний гибких диссипативных систем в виде гибких упругих балок с произвольными краевыми условиями, находящихся под действием поперечной знакопеременной и ударной нагрузок, с помощью которых проведен качественный анализ состояний балки Эйлера-Бернулли, С.П. Тимошенко, Пелеха-Шереметьева.
На защиту выносятся следующие результаты и положения:
- Построены математические модели нелинейной динамики гибких балок, подчиняющихся гипотезам Эйлера-Бернулли, С.П. Тимошенко, Пелеха-Шереметьева. Это позволило впервые провести сравнительный анализ состояния системы при использовании различных математических моделей, а также определять зоны хаотических колебаний и управляющие параметры (амплитуду вынуждающих колебаний, скорость груза и отношение массы груза к массе балки в случае ударной нагрузки, коэффициенты внешнего и внутреннего трения, геометрические параметры балки) для предотвращения негативных последствий.
- Разработан и реализован в виде пакета программ для ПЭВМ универсальный алгоритм расчета гибких балок при действии произвольной нагрузки, в том числе ударной, с учетом различных видов трения и оснований Винклера и В.З.Власова, с помощью которого проведен качественный анализ (сигналов
, фазовых портретов, сечения Пуанкаре, спектров мощности на основе Фурье анализа, автокорреляционных функций, знаков Ляпуновских показателей) хаотических колебаний гибких диссипативных систем в виде балок при различных типах граничных условий. - Впервые построенные карты зависимости характера колебаний системы от управляющих параметров
, под действием знакопеременной поперечной нагрузки вида 
для каждой модели с учетом некоторых типов граничных условий, при нескольких значениях геометрического параметра, дают представление о динамическом состоянии системы для каждого набора управляющих параметров. - Для статической и динамической задач установлено существенное влияние относительной толщины балки на напряженно - деформированное состояние системы.
- Показано, что переход колебаний из гармонических в хаотические для гибких балок при действии поперечной знакопеременной и ударной нагрузок на разных частотах вынуждающих колебаний может происходить по различным сценариям, таким как: сценарий Фейгенбаума, сценарий Рюэля, Такенса, Ньюхауза, модифицированный сценарий Рюэля, Такенса, Ньюхауза, модифицированный сценарий Помо - Манневиля.
Практическая ценность и реализация результатов. Предложенные математические модели позволяют решать широкий класс практических задач динамики нелинейных гибких упругих балок при произвольных краевых условиях. Разработанный алгоритм позволяет исследовать колебания механических систем в виде балок в зависимости от управляющих параметров (амплитуды и частоты возбуждающей нагрузки, краевых условий, геометрического параметра, диссипативных членов, упругого основания). Предложенные алгоритмы расчета динамики нелинейных упругих балок могут быть использованы при создании балочных систем для инженерных конструкций, в приборостроении при инженерных расчетах. Результаты, полученные в диссертации, использовались при чтении курса лекций «Математические модели и методы исследования сложных колебаний неклассических распределенных механических систем», а так же при написании учебных пособий по данной теме. Работа выполнена при финансовой поддержке гранта 2006-2008 гг. РФФИ № 06-08-01357 и гранта СГТУ 1.3.08.2008 г.
Апробация работы. Основные положения и результаты диссертации представлялись на III Международной конференции по теории нелинейной динамики механических и биологических систем (Саратов, 2004); XXI Международной конференции по теории оболочек и пластин (Саратов, 2005); 8th conference of Dynamical systems «Theory and Application, DSTA 2005» (Lodz, Poland, 2005); научном семинаре кафедр «Автоматика и биомеханика» и «Прочность материалов и конструкций» Технического университета г. Лодзь (Польша, 2006); третьей Международной научно-технической конференции молодых ученых и студентов (Украина, г. Донецк, 2007); 9th conference of Dynamical systems «Theory and Application, DSTA 2007» (Lodz, Poland, 2007); IV Всероссийской научной конференции с международным участием «Математическое моделирование и краевые задачи» (Самара, 2007); VII Межрегиональной научно-практической конференции студентов и аспирантов (Новокузнецк, 2007); Международном семинаре, посвященном памяти заслуженного деятеля науки ТАССР проф. А.В.Саченкова (Казань, 15-17 сентября 2008).
В законченном виде диссертация докладывалась на научных семинарах кафедры «Высшая математика» Саратовского государственного технического университета под руководством заслуженного деятеля науки и техники РФ, д.т.н., профессора В.А.Крысько, (Россия, 2008); кафедры «Автоматика и биомеханика» Технического университета г. Лодзь (Польша, 2008) под руководством профессора Я.Аврейцевича; на научном семинаре Lodzkiego Polskiego Towarzystwa Mechaniki Teoretycznej i Stosowanej г. Лодзь (Польша, 2008) под руководством профессора Катаржины Коваль-Михальской; на межкафедральном семинаре «Математическое моделирование, численные методы и комплексы программ» Саратовского государственного технического университета под руководством заслуженного деятеля науки РФ, д.ф.-м.н., профессора В.Б.Байбурина (Россия, 2008).
Публикации. Основное содержание и результаты диссертации опубликованы в 12 печатных работах, в том числе, 4 в журналах из перечня ВАК РФ.
Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения, списка использованной литературы. Работа содержит 164 страницы, в том числе 19 рисунков, 49 таблиц. Список использованной литературы включает 123 наименования.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность диссертационного исследования, дается исторический обзор результатов по математическому моделированию балочных систем, сформулирована цель работы, приводится краткое содержание диссертации.
В первой главе формулируются основные гипотезы и допущения для построения математической модели сложных колебаний балок Эйлера-Бернулли. Рассматриваются однослойные упругие балки, подчиняющиеся закону Гука. Проводится сравнительный анализ методов типа конечных разностей для сведения уравнений в частных производных к задаче Коши. Разработан алгоритм расчета динамики балок Эйлера-Бернулли методами конечных разностей с аппроксимацией 
по пространственной координате и конечных элементов с аппроксимацией по Бубнову-Галеркину. Обеспечивается достоверность получаемых результатов.
Как двумерная область балка определяется следующим образом: 
, 
. Рассматривается балка с прямоугольным поперечным сечением длиной 
, высотой 
, нагруженная распределенной на единицу длины поперечной знакопеременной нагрузкой 
, действующей в направлении оси oz. Колебания балки рассматриваются в рамках гипотезы Эйлера-Бернулли, т.е. в предположении, что поперечные сечения, перпендикулярные оси балки до изгиба, остаются плоскими и перпендикулярными изогнутой оси и не деформируются в своей плоскости. Система дифференциальных уравнений в перемещениях в безразмерном виде, описывает движение балки с учетом диссипации энергии (1), черточки над безразмерными параметрами, для простоты опущены:
 | (1) |
где 
, 
, 
- нелинейные операторы. Для сведения уравнений (1) к безразмерному виду использовались следующие безразмерные параметры: 
, 
, 
, 
, 
, 
, 
, 
.
К уравнениям (1) следует присоединить одно из краевых условий и начальные условия:
1. Заделка – заделка:
 | (2) |
2. Шарнир – шарнир:
 | (3) |
3. Заделка – шарнир:
 | (4) |
4. Заделка – свободный край:
 | (5) |
Начальные условия:
 | (6) |
На рассматриваемую балку действует знакопеременная поперечная нагрузка вида
 | (7) |
Бесконечномерная задача (1)-(6) с помощью метода конечных разностей с аппроксимацией 
, а также метода конечных элементов в представлении Бубнова-Галеркина, сводится к конечномерной системе обыкновенных дифференциальных уравнений. Системы обыкновенных дифференциальных уравнений решаем методом Рунге – Кутта четвертого порядка точности.
Важным вопросом при решении задач хаотической динамики является вопрос достоверности получаемых результатов. В связи с тем, что дифференциальные уравнения, описывающие поведение используемых математических моделей, существенно не линейны, получение аналитических решений для проверки численных результатов становится невозможным. Данный тип динамических задач решается впервые, таким образом, сопоставление с результатами других авторов невозможно. Единственным выходом является сопоставление с результатами, полученными с помощью разных численных методов. Поэтому используем сравнение решений по методу конечных разностей с аппроксимацией 
(МКР) и методу конечных элементов (МКЭ) в форме Бубнова-Галеркина.
Отметим, что для каждого метода, используя принцип Рунге, проводились исследования на предмет сходимости результатов по разбиениям: а) пространственной координаты; б) времени. Было установлено, что оптимальным разбиением по пространственной координате является 
, шаг по времени 
. На основании этих исследований выбирались параметры задачи, используемые в дальнейшем. Показана сходимость численных результатов с использованием МКР и МКЭ на примере не только сходимости сигнала 
, но и фазовых портретов, спектров мощности в представлении Фурье 
.
Для подтверждения достоверности результатов и следуя идее А.Пуанкаре о том, что лучше изучать все многообразие орбит, чем следить за какой-то конкретной, были построены карты динамических режимов для управляющих параметров 
, которые представляют собой графическое отображение результатов решения динамической задачи. Важным вопросом является наиболее полное отображение информации при минимальных затратах машинного времени. Для построения карт на область пространства 
была наложена сетка, в узлах которой производилась идентификация характера колебаний. Предварительно исследовался вопрос о сходимости решения при увеличении количества разбиений области 
. Расчеты показали, что 
является оптимальным. Для построения одной карты, с разрешением 
необходимо просчитать 
вариантов. Время счета одного варианта составляет примерно 2.5 часа.
| МКР | МКЭ | Условные обозначения |
 а) |  б) |  |
Рис. 1 Карты динамических режимов
Для консольной балки приведены карты динамических режимов, полученные двумя разными методами: конечных разностей (рис.1а) и конечных элементов (рис. 1б).
Условные обозначения применимы ко всем нижеприведенным картам. Как видно, карты для каждого набора управляющих параметров, построенные разными методами, практически полностью совпадают. Для всех граничных условий было проведено сравнение численных результатов по методам исследования (сопоставление сигналов, спектров мощности, карт), в результате чего сделан вывод о достоверности получаемых численных решений.
Во второй главе исследуются сценарии перехода гармонических колебаний в хаотические для гибких балок Эйлера-Бернулли в условиях поперечной знакопеременной нагрузки для различных граничных условий.
Так, исследован сценарий перехода системы к хаосу для граничных и начальных условий (3), (6), в результате чего выявлено, что система переходит к хаотическим колебаниям через последовательное появления трех частот, значения которых линейно не зависят друг от друга, то есть наблюдаем модифицированный сценарий Рюэля, Такенса, Ньюхауза, который был впервые предложен В. А.Крысько и И.В.Кравцовой. Для граничных и начальных условий (5), (6) переход системы к хаосу происходит по следующему сценарию: колебания системы на частоте возбуждения, появление одной линейно независимой частоты, обогащение спектра рядом частот, линейно зависящих от первых двух. Если рассматривать сценарий в случае заделки обоих концов балки (2) и начальных условий (6), то для частот возбуждающих колебаний 
и 
сценарии различны.
Проведено сравнение результатов расчета гибких балок Эйлера-Бернулли для прогиба и перемещений при некоторых граничных условий. На рис. 2 представлены карты динамических режимов при граничных условиях (4).
Исследуется влияние коэффициентов диссипации энергии 
на частотные характеристики по длине балки для 
. При 
и 
количество и значения частот по длине балки совпадает, то есть эти значения являются наиболее оптимальными при расчетах. Параллельно исследовалось влияние относительной толщины балки на характер колебаний гибких балок.
Для учета влияния Кулоновского, нелинейного и линейного трения, а также упругих оснований Винклера и В.З.Власова, уравнения движения, в безразмерном виде (модель Эйлера-Бернулли), запишем как
 | (8) |
Нелинейные операторы 
аналогичны приведенным выше. Дополнительные безразмерные параметры: 
, 
. Система уравнений (8) является дважды нелинейной. С одной стороны, она учитывает нелинейную зависимость между деформациями и перемещениями, а с другой - нелинейность сил трения от скорости, которая описывается моделью Кулона (при параметре m = 0). Для исследования колебаний при линейном трении следует положить m = 1, при нелинейном m = 2.
К уравнению (8) следует присоединить граничные (2) – (4) и начальные условия.
Дано сравнение динамического состояния балки Эйлера-Бернулли на основании Виклера и без его учета (рис. 3) для балки с шарнирным опиранием краев. Как видно, учет основания Винклера в математи-ческой модели приводит к существенному измене-нию характера колебаний балки.
Третья глава посвящена изучению колебаний гибких балок Эйлера-Бернулли при действии продольного удара груза массой 
. В момент удара на свободный край балки, согласно второму закону Ньютона, действует сила: 
. Приводя данное уравнение к безразмерному виду, получим 
.
Тогда, к уравнениям (1) присоединим уравнения на границе в безразмерном приведении:
 | (9) |
где 
, и начальные условия:
 | (10) |
где 
- скорость груза в момент удара.
Данная бесконечномерная задача с помощью метода конечных разностей с аппроксимацией 
сводится к конечномерной системе обыкновенных дифференциальных уравнений, которую также решаем методом Рунге – Кутта четвертого порядка точности.
При увеличении скорости груза V, система переходит к хаосу по сценарию Фейгенбаума. Получена первая константа Фейгенбаума, относительная погрешность которой составляет 0.005.
Анализ характера распространения волны 
по длине балки позволяет говорить о том, что характер волны 
не зависит от скорости груза, но ее амплитуда увеличивается. При ударе по торцу балки грузом массой 
со скоростью 
происходит процесс распространения продольных волн и их отражения. Время начала резкого увеличения амплитуды поперечных колебаний примем за 
.
Исследование вопроса о зависимости 
от геометрических параметров балки и от скорости груза в момент удара позволяет сделать вывод о том, что значение 
значительно больше в случае относительно толстых балок 
и напрямую зависит от скорости груза.
В четвертой главе формулируются основные гипотезы и допущения для построения математической модели гибкой балки С.П.Тимошенко.
Размерная система дифференциальных уравнений в перемещениях с учетом вышеприведенных безразмерных параметров и дополнительных 
, 
сводится к системе обыкновенных дифференциальных уравнений в перемещениях описывающих движения балки с учетом диссипации энергии:
 | (11) |
где 
- нелинейные операторы аналогичные используемым в уравнениях (1), а 
.
К системе дифференциальных уравнений (11) следует присоединить одно из нижеследующих граничных условий, и начальные условия.
1. Заделка – заделка:
 | (12) |
2. Шарнир – шарнир:
 | (13) |
3. Заделка – шарнир:
 | (14) |
4. Заделка – свободный край:
 | (15) |
Начальные условия:
 | (16) |
Бесконечномерную задачу с помощью метода конечных разностей сводим к конечномерной системе обыкновенных дифференциальных уравнений. Применяя процедуру метода Бубнова–Галеркина, получаем разрешающие уравнения МКЭ. Отметим, что аппроксимирующие функции для прогиба и перемещений были использованы исходя из свойств математической модели. Сборка матриц масс, демпфирования и формы производится так же как и для модели Эйлера-Бернулли.
На рассматриваемую балку действует знако-переменная поперечная нагрузка вида (7). Чтобы убедиться в достоверности получаемых численных результатов, решаем описан-ную систему уравнений методами конечных элементов и конечных разностей, при всех граничных условиях. В целях сокращения затрат времени, можно проанализировать характер колебаний системы, только для одного значения частоты, по всему интервалу нагрузки. Для этого строим шкалы зависимости характера колебаний от амплитуды внешней нагрузки. На рис.4 показаны шкалы и графики зависимости 
, для граничных условий (12).
На основании этих результатов, можно сказать, что при значениях амплитуды нагрузки, соответствующей нехаотическим колебаниям, значения 
по МКР и МКЭ совпадают полностью. После перехода к хаосу значения максимального прогиба отличаются, но несущественно.
В работе выявлено явление динамической потери устойчивости балки под действием поперечной знакопеременной нагрузки, которое характеризуется резким изменением значения максимального прогиба при незначительном изменении амплитуды вынуждающих колебаний. Явление динамической потери устойчивости хорошо видно при переходе системы от т. A к т. В и от т. С к т. D, что на шкале иллюстрируется сменой режима колебаний. При переходе от т. Е к т. F происходит обратный переход системы от хаоса к гармоническим колебаниям, что также можно наблюдать на шкале. В этом случае значения прогиба уменьшаются в 1,5 раза.
Исследование влияния граничных условий на характер поведения системы, в частности анализ карт динамических режимов, дает основание говорить об их существенном влиянии (рис. 5).
а) заделка-заделка б) шарнир-шарнир в) заделка-шарнир
Рис. 5 Карты динамических режимов для разных видов граничных условий
Для несимметричных граничных условий (14) исследованы сценарии перехода к хаотическим колебаниям для частот 
. Выявлено, что в области управляющих параметров 
нет единого сценария перехода колебаний из гармонических в хаотические. Существуют подобласти 
, в которых переход совершается по различным сценариям.
Исследование колебаний нелинейных конструкций на различных режимах (от гармонических до хаотических) предъявляет повышенные требования как к численным методам, так и к математической модели. Учет поперечного сдвига при построении математической модели нелинейных колебаний балки приводит к качественно иной картине характера колебаний системы, что хорошо видно при сравнении карт.
Анализ шкал колебаний, графиков 
и карт динамических режимов колебаний для моделей Эйлера-Бернулли и С.П.Тимошенко показывает, что результаты, полученные по этим моделям, существенно различаются.
Пятая глава посвящена исследованию сложных колебаний гибких балок Пелеха-Шереметьева. Приведены основные гипотезы и допущения, построена математическая модель гибкой балки Пелеха-Шереметьева.
Система обыкновенных дифференциальных уравнений в перемещениях, описывающих движения балки с учетом диссипации энергии в безразмерном приведении, имеет вид:
 | (17) |
Нелинейные операторы 
и безразмерные параметры, те же, что и выше. К системе дифференциальных уравнений (17), следует присоединить одно из нижеследующих граничных условий и начальные условия.
1. Заделка – заделка:
 | (18) |
2. Шарнир – шарнир:
 | (19) |
3. Заделка – шарнир:
 | (20) |
Начальные условия:
 | (21) |
Достоверность численного решения обеспечена сопоставлением результатов по методам конечных разностей и конечных элементов.
а)модель Эйлера-Бернулли б)модель С.П.Тимошенко в)модель Пелеха-Шереметьева
Рис. 6 Карты динамических режимов для различных моделей
Очень важным моментом является правильный выбор математической модели исследуемой распределенной системы. Рассмотрим этот вопрос на примере выбора той или иной расчетной схемы для однослойных изотропных балок 
(рис. 6).
Анализ этих карт позволяет сделать вывод, что учет поворота нормали приводит к существенному изменению режимов колебаний (сопоставление моделей Эйлера-Бернулли и С.П.Тимошенко, Эйлера-Бернулли и Пелеха-Шереметьева). В то время когда учет искривления нормали (сопоставление моделей С.П.Тимошенко и Пелеха-Шереметьева) приводит к изменениям в основном на высоких частотах.
В заключении приводятся основные результаты и выводы по работе.
Основные выводы по диссертации
- Построены общие теории и математические модели сложных колебаний гибких балок Эйлера-Бернулли, С.П.Тимошенко, Пелеха-Шереметьева, проведен качественный анализ динамического поведения рассматриваемых систем.
- Предложен эффективный алгоритм решения поставленных задач. Разработан и реализован комплекс программ анализа хаотических колебаний балок с некоторыми краевыми условиями, находящихся под действием поперечной знакопеременной и ударной нагрузок.
- Разработан комплекс программ для качественного исследования сложных колебаний балок с помощью метода конечных разностей с аппроксимацией
и метода конечных элементов с аппроксимацией по Бубнову-Галеркину. - Обоснован выбор типа метода конечных разностей для сведения уравнений в частных производных к задаче Коши.
- Проведено исследование сходимости метода конечных разностей в зависимости от числа разбиений по пространственной координате для балок Эйлера-Бернулли при действии поперечной знакопеременной нагрузки.
- Выявлены области сценария Фейгенбаума на картах
для консольной балки Эйлера-Бернулли при действии ударной нагрузки грузом массой 
, где происходило до 4 бифуркаций Хопфа, что позволило вычислить константу Фейгенбаума с относительной погрешностью 0.005. - Дается сопоставление решений, полученных методом конечных разностей и методом конечных элементов для каждой из моделей при действии поперечной знакопеременной нагрузки.
- Рассмотрено влияние некоторых типов трения (кулоновское, нелинейное, линейное) на балку Эйлера-Бернулли и исследовано поведение балок Эйлера-Бернулли на упругих основаниях Винклера и В.З.Власова.
- В соответствии с известными сценариями перехода колебаний балочных конструкций в хаос проведена классификация колебаний балок, находящихся под действием поперечной знакопеременной и продольной ударной нагрузки. Выявлены и исследованы сценарии Фейгенбаума, Рюэля, Такенса, Ньюхауза, модифицированные Рюэля, Такенса, Ньюхауза, Помо-Манневиля, характерные для колебаний исследуемых систем, и выявлены их области на картах динамических режимов.
- Для каждой рассматриваемой модели были отмечены явления динамической потери устойчивости при действии знакопеременной поперечной нагрузки, что подтверждается резким увеличением максимального прогиба при малом изменении амплитуды вынуждающих колебаний.
- Исследовано влияние геометрического параметра
на характер поведения балки для каждой математической модели. Выявлено, что по мере увеличения параметра 
результаты, получаемые по моделям Эйлера-Бернулли, С.П.Тимошенко, Пелеха-Шереметьева, сходятся. - Проведено качественной сравнение результатов, полученных для каждой математической модели, что позволило сделать вывод о пределах применимости каждой модели в зависимости от геометрических параметров балки, это позволит при расчете конструкций избегать ситуаций потери устойчивости системы.
Публикации по теме диссертации
Работы, опубликованные в ведущих научных журналах и изданиях, определенных ВАК РФ
- Салтыкова О.А. Нелинейная динамика балок Эйлера-Бернулли и типа Тимошенко / В.А. Крысько, М.В. Жигалов, О.А. Салтыкова // Известия вузов. Машиностроение. 2008. № 6. С. 7-27.
- Салтыкова О.А. Особенности сложных хаотических колебаний балок Эйлера-Бернулли и типа Тимошенко в зависимости от краевых условий / В.А. Крысько, М.В. Жигалов, О.А. Салтыкова // Известия вузов. Строительство. 2008. № 9. С. 4-10.
- Салтыкова О.А. Управление сложными колебаниями нелинейных многослойных балок / А.В. Крысько, М.В. Жигалов, О.А. Салтыкова // Известия вузов. Авиационная техника. 2008. №3. С. 10-13.
- Салтыкова О.А. Диссипативная динамика геометрически нелинейных балок Бернулли–Эйлера / В.А. Крысько, М.В. Жигалов, О.А. Салтыкова, А.С. Десятова // Известия РАН. Механика твердого тела. 2008. №6. С. 128 – 136.
Работы, опубликованные в других изданиях
- Салтыкова О.А. Нелинейные колебания балки Эйлера-Бернулли под действием продольного удара груза массой
. / О.А. Салтыкова, В.А. Крысько // Труды третьей Международной научно-технической конференции молодых ученых и студентов 11-13 декабря 2007 года, г. Донецк (Украина). Донецк, 2007. С. 489-492. - Салтыкова О.А. Нелинейные колебания гибкой балки модели С.П. Тимошенко / В.А. Крысько, А.М. Варыгин, О.А. Салтыкова // Нелинейная динамика механических и биологических систем: межвузовский научный сборник. Саратов, 2004. С. 205-212.
- Салтыкова О.А. Сложные колебания гибких балок при продольном ударе / В. А. Крысько, А. М. Варыгин, О. А. Салтыкова // Труды XXI Международной конференции по теории оболочек и пластин. Саратов. 2005. С.288-294.
- Салтыкова О.А. Сложные колебания гибких балок для некоторых типов краевых условий // VII Межрегиональная научно-практическая конференция студентов и аспирантов. Новокузнецк, 2007. С. 14-19.
- Салтыкова О.А. Математическая модель нелинейной динамики балок с учетом поперечных сдвигов / М.В. Жигалов, О.А. Салтыкова // Математическое моделирование и краевые задачи: труды четвертой Всероссийской научной конференции с международным участием. Самара, 2007. С. 228-231.
- Салтыкова О.А. Нелинейная динамика упругих балок моделей Бернулли-Эйлера, С.П.Тимошенко и Шереметьева-Пелеха / В.А.Крысько, А.В.Крысько, М.В.Жигалов, О.А.Салтыкова // Актуальные проблемы нелинейной механики оболочек: материалы Международного семинара, посвященного памяти заслуженного деятеля науки ТАССР проф. А.В.Саченкова. Казань, 15-17 сентября 2008. Казань, 2008. С. 81-84.
- Saltykova O.A. Vibration of Flexible Beam Subjected to a Longitudinal Impact / V.A.Krysko, J.Awrejcewicz, Yu.V.Chebotarevskiy, O.A.Saltykova // 8th Conference on Dynamical Systems Theory and Applications. December 12-15. 2005. Lodz, Poland, 2005. P. 719-727.
- Saltykova O.A. Analysis of Regular and Chaotic Dynamics of the Euler-Bernoulli using Finite Difference and Finite Element Methods / A.V. Krysko, M.V. Zhigalov, J. Awrejcewicz, O.A. Saltykova // 9th Сonference of Dynamical Systems Theory and Application. December 17-20. 2007. Lodz, Poland, 2007. P.657-669.
