Математическое и компьютерное моделирование нелинейных распределенных механических систем тема диссертации и автореферата по ВАК РФ 05.13.18, доктор физико-математических наук Жигалов, Максим Викторович

• Жигалов, Максим Викторович
• доктор физико-математических наук
• 2013

Оглавление диссертации доктор физико-математических наук Жигалов, Максим Викторович

Введение. Исторический обзор.

Глава 1. Метод линеаризации и понижения порядка нелинейных дифференциальных уравнений математических моделей распределенных механических структур.

§1.1 Метод линеаризации для математических моделей, содержащих бигармонический оператор. Уравнения Кармана. Доказательство сходимости итерационной процедуры.

§1.2 Метод понижения порядка дифференциального оператора для математических моделей, содержащих бигармонический оператор.

§ 1.3 Метод линеаризации и понижения порядка для математических моделей на основе гипотез Кирхгофа, С.П. Тимошенко, Шереметьева - Пелеха.

§ 1.4. Метод линеаризации и понижения порядка для нелинейных уравнений, содержащих бигармонический или гармонический операторы.

§1.5 Численные примеры.

Выводы по главе.

Глава 2. Численные методы решения уравнений, содержащих бигармонический и гармонический операторы, для областей с криволинейной границей.

§ 2.1 Итерационная процедура сведения бигармонического уравнения к уравнению типа Пуассона для областей с криволинейной границей.

Доказательство сходимости итерационной процедуры.

§ 2.2 Численные процедуры методов конечных и граничных элементов для двумерных задач с криволинейной границей.

§ 2.3 Результаты решений методами конечных разностей, конечных и граничных элементов уравнения кручения стержня. Переход от уравнения Пуассона к уравнению Лапласа.

§ 2.4 Численное решение бигармонического уравнения. Сравнение результатов для процедуры понижения порядка, полученных различными численными методами.

Выводы по главе.

Глава 3. Нелинейная динамика однослойных балок, математической модели

Бернулли — Эйлера.

§3.1 Математическая модель нелинейной балки Бернулли - Эйлера.

§ 3.2 Сценарии перехода от гармонических колебаний к хаотическим.

Карты характеров колебаний.

§ 3.3 Исследование математической модели нелинейной динамики балок

Бернулли - Эйлера с помощью вейвлет-преобразования.

Выводы по главе.

Глава 4. Нелинейная динамика однослойных балок, математических моделей

С.П. Тимошенко и Шереметьева - Пелеха.

§4.1 Математические модели нелинейных балок С.П. Тимошенко и Шереметьева - Пелеха.

§ 4.2 Сценарии перехода от гармонических колебаний к хаотическим.

Карты характеров колебаний. Сравнение моделей.

§ 4.3 Исследование математических моделей нелинейной динамики балок С.П. Тимошенко и Шереметьева - Пелеха с помощью вейвлет преобразования.

Выводы по главе.

Глава 5. Нелинейная динамика многослойных балок.

§ 5.1 Математические модели многослойных балок с использованием гипотез

Бернулли - Эйлера, С.П. Тимошенко, Шереметьева - Пелеха.

§ 5.2 Анализ нелинейной динамики многослойных балок. Карты характеров колебаний.

§ 5.3 Математические модели контактных задач балок с использованием гипотез

Бернулли - Эйлера, С.П. Тимошенко, Шереметьева - Пелеха.

§ 5.4. Численное исследование балок с учетом трех типов нелинейности.

Выводы по главе.

Глава 6. Нелинейная динамика пластин и оболочек.

§6.1 Математические модели пластин и оболочек, основанные на гипотезе

Бернулли - Эйлера.

§ 6.2 Вейвлет-анализ при исследовании цилиндрических панелей и замкнутых цилиндрических оболочек.

§ 6.3 Анализ нелинейной динамики пластины произвольного очертания.

Выводы по главе.

Введение диссертации (часть автореферата) на тему «Математическое и компьютерное моделирование нелинейных распределенных механических систем»

При описании физического явления обычно вводят в рассмотрение гипотезы, строят функционал и из минимума его получают некоторую систему дифференциальных уравнений, обыкновенных или в частных производных, справедливую в определенной области, а также краевые и начальные условия. На этой стадии математическая модель замкнута, и для практического применения требуется получить решение, чаще всего численное, для конкретных значений числовых данных. Здесь, однако, возникают основные трудности, так как точному решению существующими аналитическими методами поддаются лишь уравнения самого простого вида внутри геометрически тривиальных границ.

В связи с этим возникла потребность в методах, позволяющих упростить исходные дифференциальные уравнения. Эти методы условно можно разбить на три большие группы: 1) линеаризация уравнения; 2) понижение размерности искомой функции; 3) понижение порядка дифференциального оператора.

Методы линеаризации

Существующие методы решения нелинейных задач в зависимости от уровня, на котором происходит линеаризация, можно разделить на две группы. Первая - линеаризация систем дифференциальных уравнений, вторая — линеаризация алгебраических уравнений, получающихся в результате применения к исходным дифференциальным методов дискретизации. Далее рассмотрены методы первой группы.

Одними из представителей этой группы являются методы Ньютона [1] и Ньютона-Канторовича [1]. Суть их в следующем: пусть задано нелинейное уравнение

Вд = 0, (1) с нелинейным, дифференцируемым по Фреше оператором, действующим из некоторого множества О банахова пространства Е\ в банахово пространство Е2. Если определен действующий из Е] в Е2 линейный оператор

Xм]}-1, (2) то для решения исходного дифференциального уравнения применима следующая итерационная процедура

Модификацией этой итерационной процедуры является следующая где х(0) - начальное приближение.

Метод Ньютона-Канторовича использован в статье Ю.Н. Санкина [2] для решения системы нелинейных дифференциальных уравнений, описывающих большие деформации оболочек вращения при симметричной нагрузке. При этом метод применен ко всей системе. Для нахождения матрицы функциональных производных предлагается использовать либо численную процедуру типа разностной, либо непосредственное дифференцирование операторов. Эффективность метода подтверждена решением задачи о нелинейном деформировании плоской мембраны. Тот же метод применен в статье В.Н. Мальгина [3] для создания алгоритмов, позволяющих решать задачи прочности, устойчивости и колебаний оболочек вращения, математических моделей, построенных на гипотезах типа С.П. Тимошенко. Метод Ньютона-Канторовича применяется последовательно к каждому из шести обыкновенных нелинейных уравнений. Линейные уравнения решаются методом ортогональной прогонки. Автором отмечается, что контрольные расчеты по всем программам показали хорошую сходимость алгоритмов и их устойчивость. В статье

Н.В. Валишвили [4] с помощью метода Ньютона исходная нелинейная краевая задача теории оболочек сводится к задаче Коши, для решения которой используется метод Рунге-Кутта. Что касается доказательств сходимости методов Ньютона и Ньютона-Канторовича для нелинейных уравнений, то они даны в статьях В.М. Вержбицкого [5] и А.В. Машкова [6]. Модификация метода Ньютона с использованием ряда Тейлора относится к другой группе и описана ниже.

Наиболее известный метод линеаризации — метод последовательных приближений. Существуют две модификации этого метода. В первой из нелинейного оператора выделяется линейный оператор со старшими производными. Все остальные составляющие исходного уравнения перебрасываются в правую часть и их формирование происходит за счет значений искомых функций, полученных на предыдущих итерациях.

Наиболее часто встречающимся является метод простой итерации. Суть его в следующем. Для решения нелинейной системы т=/, (5) предлагается представить оператор в виде

6) где Ьх - линейный оператор старшей степени, а Ь2 - линейный оператор младшей степени, Ь3 - нелинейный оператор. В результате получается следующая итерационная процедура:

Метод, предложенный в 1958 году В.З. Власовым и получивший название метода последовательных нагружений, заключается в представлении нагрузки в виде суммы отдельных ступеней. В пределах каждой ступени нагрузки считаются линейные теории рассматриваемых объектов (пластинки, оболочки и др.). Величины ступеней нагрузки выбираются так, чтобы, например, отношение прогиба к толщине пластинки укладывалось в рамки лимитной теории. Таким образом, решение системы нелинейных уравнений сводится к последовательному решению нескольких систем линейных уравнений. Этот метод был разработан В.В. Петровым в статье [7]. Этим методом решены как геометрически, так и физически нелинейные задачи, контактные задачи и задачи для массивных тел. В работах Э.И. Григолюка, В.И. Шалашилина, В.И. Кузнецова [8, 9] предложено уточнять решение, полученное методом последовательных нагружений, методом наискорейшего спуска, разработанным JI.B. Канторовичем, или методом Ньютона.

Одной из разновидностей методов последовательных приближений является метод квазилинеаризации, описанный в книге R.E. Bellman, R.E. Kalaba [10]. Этот метод представляет собой дальнейшую разработку метода Ньютона и его обобщенного варианта, предложенного JI.B. Канторовичем. Применение этих методов описано в статьях Я.М. Григоренко, А.Т. Василенко, H.H. Крюкова и многих других [11, 12].

Другой разновидностью методов последовательных приближений являются методы спуска, которые впервые были описаны G. Temple [13] и JI.B. Канторовичем [14]. Они состоят в том, что для решения функционального уравнения

А-ЯВ)Х-Р = 0, (8) отыскивают на каждом шаге итерационного процесса минимум одной из форм: Oí - потенциальная энергия системы, Ф2 — сумма квадратов левых частей уравнений или других на некотором подпространстве минимизации. Помимо геометрической интерпретации, метод спуска может быть истолкован «механически». Он имитирует движение системы, когда она освобождена в начальном положении, не представляющем положение равновесия. Для такой интерпретации достаточно номеру шага процесса присвоить название «дискретного времени». Таким образом, статическая система заменяется динамической с искусственно введенным временем, что обеспечивает ранее недостижимые вычислительные результаты. Классификация методов спуска и обзор работ по использованию их в задачах механики приведены в работах Д.В. Вайнберга и С.В. Симеонова [15-17].

Методы, принадлежащие третьей группе, используют представление искомых функций в виде суммы известного решения и малой добавки, играющей роль уточнения решения. Это представление дает возможность линеаризовать исходное уравнение. Одним из таких методов является метод, изложенный в статьях В.А. Пухлий, G.A. Thurston и других [18-21]. Суть его в следующем. Пусть дано нелинейное уравнение y\x) + f{y{x\x) = V. (9)

Будем искать решение в виде

10) где Jo - заданное решение, А - поправка. Подставляя в исходное уравнение и разлагая нелинейные члены в ряд Тейлора в окрестности точки (х0,у0), получаем + A" + /(70,x0) + /'(7o^o)0;i-J;o) + - = 0. (11)

Отбрасывая слагаемые, начиная с

0>,~Л)2=Д2, (12) приходим к линейному уравнению

A" + kA = F(y0,x0). (13)

Решая это уравнение, находим поправку, которую используем для нахождения первого приближения. Затем процесс повторяется до достижения заданной точности. Статья J. Griiters [22] посвящена применению аналогичного подхода к решению задачи для тела типа диска из анизотропного материала, нагруженного по краю. Нулевое приближение получается из решения задачи для изотропного тела.

Разновидностью методов, использующих малые добавки к решению, являются методы малого параметра. Д.Ф. Давиденко в [23, 24] предложен метод решения систем нелинейных уравнений, содержащих параметр, в случае, когда одно из значений параметра известно. Изложим кратко существо этого метода. Пусть дана система уравнений к(хх,х2,.,хп,Л) = Ъ £ = 1,2,.,и. (14)

Пусть также для л = л0 известны такие х;0, что

Ук(х10,х20,.,хп0,Л10) = к = 1,2,.,п. (15)

Функции /к определены и непрерывны в некоторой (и+1)-мерной области О изменения хх,х2,.,хп,Л и точка (х10,х20,.,хп0,Ад) е О. Пусть также в области (7 Якобиан

1(Х1,Х2,.,ХП,Л) = ^'/2'''''-(16) и\Х^,Х2,.,Хп)

Параметр Л принимается за независимую переменную, и хх,х2,.,хп считаются функциями Я. Исходная система дифференцируется по Л. В результате получается система обыкновенных уравнений, линейных относительно производных: к=1'2'-'п- <17>

1 дх1 аЛ дЛ

Так как якобиан этой системы отличен от нуля, она может быть решена

1х1 относительно производных —В итоге задача сводится к системе уравнении с1Л вида с!х

1- = Е({хх,.,хп,Л), / = 1,2,.72. (18) аЛ

К этой системе присоединяются граничные условия х,(Л0) = х109 / = 1,2,., л. (19)

Полученная задача является задачей Коши для переменной X. Интегрирование этой задачи любым из известных методов для значений на отрезке С даёт приближенное решение для исходной системы. Описанный метод нашел широкое применение для решения разнообразных задач механики, например в работах Э.И. Григолюка, М.И. Карасика и В.И. Шалашилина, А.И. Уздалева [25-29].

Описанные выше подходы линеаризуют исходную задачу, т.е. сводят ее к решению последовательности линейных задач. Однако все эти подходы не понижают порядок исходного дифференциального оператора, что приводит к дополнительным трудностям при численной реализации указанных подходов.

Остается еще одна важная проблема - размерность пространства исходной задачи. В связи с тем, что решение проводится в основном численными методами, размерность пространства существенно влияет на затраты машинного времени и точность полученного решения. Поэтому важным направлением являются методы понижения размерности пространства исходной задачи.

Методы понижения размерности

Численное решение трехмерных по пространственным координатам задач теории упругости даже на современном этапе развития вычислительной техники сопряжено со значительными трудностями.

Поэтому наиболее распространенным подходом трехмерных задач является понижение их размерности каким-либо приемом. На практике обычно используется один из следующих приемов.

Первый из них заключается в осреднении (интегрировании) по координате, по которой размер объекта мал по сравнению с двумя другими. Применимость такого приема связана с необходимостью знать (хотя бы приближенно) распределения параметров НДС по соответствующей координате. Если это распределение известно, например из геометрических соображений, то этот прием весьма эффективен и достаточно экономичен.

Второй из них применим для расчета тел вращения; суть его состоит в разложении параметров НДС в направлении угловой координаты (в цилиндрической системе координат) в ряды по тригонометрическим функциям.

После разложения в такие же ряды внешних наг система уравнений, описывающих объект, распадается на «К» независимых систем уравнений, где «К» - число гармоник, удерживаемых в разложении. Таким образом, вместо трехмерной задачи решаются «К» двумерных.

До середины 30-х годов XX столетия господствующим методом численного решения краевых задач был метод конечных разностей; достаточно упомянуть работы JI.B. Канторовича, В.З. Власова, A.A. Самарского и других [30-32]. И до сих пор к большинству задач математической физики, решаемых численно, применяются формулы метода конечных разностей.

Около 1935 г., однако, появляется несколько работ, посвященных решению уравнений в частных производных с двумя переменными, основной идеей которых является использование метода конечных разностей для получения решения по одной переменной, и, таким образом, решение уравнений в частных производных заменяется решением систем обыкновенных дифференциальных уравнений. Качественный тип этой системы зависит не только от вида уравнения, но и от типа граничных условий решаемой задачи. К таким работам относятся работы М. Г. Слободянского [33] для уравнений эллиптического типа и работа D. А. Hartree [34] для уравнения теплопроводности. Метод в дальнейшем стал называться методом прямых. На рубеже 60-70-х гг. XX века этот метод применялся для решения задач плоской теории упругости, теории пластин и оболочек, а также к задачам термоупругости.

Известно, что математические задачи теории упругости можно формулировать как вариационные задачи, то есть задачи о нахождении экстремума некоторых функционалов. Условия, при которых функционал имеет стационарное значение, определяются уравнениями Эйлера и естественными граничными условиями. Вариационные формулировки задач являются теоретической основой для построения прямых вариационных и вариационно-разностных методов, метода конечных и граничных элементов. Часто при решении задач применяются вариационные принципы Лагранжа и Кастильяно.

Рассмотренные методы понижения размерности позволяют свести дифференциальное уравнение в частных производных к обыкновенным дифференциальным уравнениям, которые, однако, остаются нелинейными и имеют высокий порядок дифференциального оператора.

Методы понижения порядка дифференциального оператора

Рассмотрим методы понижения порядка исходного дифференциального уравнения. Их можно разделить на три группы: введение новой функции, декомпозиции и аналогии.

Одна из статей XV. 1чГо-\уас1а [35], в которой рассмотрен изгиб полосы, свободно опертой по краям х=0 и х=а, нагруженной сосредоточенной силой в точке у=0. Дифференциальное уравнение для функции прогиба имеет вид ду4 дх2ду2 дх4 ' (20)

4 А Н где ¿г =—±,р =

Вводится новая функция ф по формуле + (21) ду дх

Подстановка функции (21) в уравнение (20) позволяет свести исходное уравнение четвертого порядка к уравнению второго порядка: д2со(х,у) „2 д2ф(х, у) '2 -0. (22) ту ex

При этом новые константы а и ß связаны со старыми следующими соотношениями: a2+ß2= 2ре2, a2ß2 = (23)

Аналогичный подход использован в работах Norman de Garrs, Е. Giangreco, L.S. Han [36-39].

Еще одним направлением среди методов понижения порядка исходного дифференциального уравнения являются методы расщепления или декомпозиции.

Суть подхода заключается в следующем: пусть требуется найти функцию и(р) точки р, удовлетворяющей в области Q уравнению

Аи{р~) = f(p), реП, (24) а на границе у краевым условиям:

Lju(s) = gj(s) sey, (25) здесь А - линейный дифференциальный оператор, /(р) - заданная в Q функция, Lj - линейные дифференциальные операторы, gj(s) - заданная на границе функция.

Предположим, что оператор А можно расчленить на несколько простейших линейных операторов, в сумме равных А, каждый из которых сравнительно легко обращается. Кроме того, соответствующее расчленение граничных условий также возможно.

Исходный оператор А представлен в виде

A = Z4, (26) к=i причем некоторые из слагаемых, входящие в операторы Ак, в исходный оператор могут и не входить. Обозначая h где f(p) = Yufk (р) » составим h вспомогательных задач относительно ик: к=1 причем граничное условие выбирают таким образом, чтобы в каждой фиксированной точке контура /оно выполнялось хотя бы для одного из ик.

Г. И. Пшеничным [40] отмечено, что в соответствии с предлагаемым методом формирование вспомогательных задач может быть выполнено бесчисленным количеством способов. Неоднозначность декомпозиции исходной задачи является следствием двух причин. Первая заключается в произволе выбора операторов ак, вторая - в произволе выбора краевых условий: вспомогательная задача может либо не содержать краевые условия, либо содержать их все, либо учитывать любую часть этих условий. Необходимо только, чтобы в каждой точке контура выполнялось краевое условие исходной задачи при объединении краевых условий вспомогательных задач. К работам в этом направлении относятся, например, статьи А.И. Уздалева, JI.A. Розина, H.D. Conway, A.W. Leissa [41-44].

Как известно, одним из свойств математического моделирования является его общность, т.е. одни и те же уравнения могут описывать различные процессы, происходящие как в природе, так и в обществе. Использование этого свойства заложено в упрощающих подходах в виде аналогий, например таких как аналогия Прандтля.

Особенный всплеск интереса к аналогии Прандтля виден с начала 50-х годов, когда для решения уравнения изгиба пластин начали использоваться методы, сводящие дифференциальное уравнение в частных производных четвертого порядка к системе уравнений второго порядка.

Au = fk(p)>

27)

Auh=fk(P)> PeQ

LjUk^gj{s) se/

28)

Ярким примером такого рода является статья Т. Nishihara и К. Tanaka [45], в которой задача об изгибе квадратной пластинки, находящейся под действием равномерно распределенной нагрузки, сводится к задаче об изгибе мембраны, т.е. уравнение четвертого порядка сводится к уравнению второго порядка. Авторами рассмотрены случаи шести типов граничных условий. Для каждого случая решены числовые примеры, показывающие достоверность предложенной методики. Использование различного рода аналогий описано в работах Ю.М. Почтмана, В.Д. Шайкевича, В. Roller, R.J. Duffin и других [46-51].

Описанные методы понижения порядка используются в основном для решения линейных задач, и важным вопросом является распространение этой методики на нелинейные задачи.

Приведенный исторический обзор методов понижения порядка и линеаризации дифференциальных уравнений в частных производных не ставит своей целью охватить весь спектр работ по данной тематике. Здесь отмечены только некоторые работы, которые отражают основную идею того или иного подхода. Более подробный обзор дан в статьях автора диссертации [52, 53].

Как видно из обзора, математические численные или аналитические методы, используемые для исследования нелинейных дифференциальных уравнений, являются более сложными и трудоемкими, чем методы линейного анализа, в связи с эти долгое время в различных областях механики изучались лишь простейшие нелинейные задачи. Развитию новых идей в понимании нелинейной механики способствовало появление компьютеров. Их использование для нелинейного анализа позволило исследовать сложное поведение механических систем. Одним из важнейших направлений является исследование хаотической динамики механических распределенных структур, причем интерес к нему, несмотря на длительную историю, не только не уменьшается, но и в последние годы возрастает, что свидетельствует о фундаментальности и актуальности этой проблемы. Исследованию нелинейной динамики в распределенных механических системах, таких как балки, пластины и оболочки, посвящено достаточно много работ, поэтому здесь приведены некоторые из них.

Стоит отметить, что работы по исследованию нелинейной динамики механических структур имеют не только научный, но и большой практический интерес, причем в различных областях. В работе Е. Kreuzer [54] изучаются нелинейные свойства морских добывающих платформ, кораблей и плавучих подъемных кранов в аэрогидродинамической среде акватории. Отмечается актуальность разработки методов нелинейного динамического расчета инженерных сооружений и технических конструкций в акватории. В работе S. Gabor [55] рассмотрена высокоскоростная обработка материалов резанием на токарном станке. В отличие от других работ, где использовался линейный анализ, исследована нелинейная вибрация в случае удвоения периода, проведено сравнение с известной бифуркацией Хопфа при вращении. Интенсификация глобального аттрактора при нестабильном резании приводит к противоречию между экспериментом и теорией, что требует более внимательного отношения к ограничениям применения метода малого параметра в условиях высокоскоростного резания. Для исследования системы мост-экипаж, движения двухосного автомобиля по шарнирно опертой балке, в работе S. S. Law [56] предложен новый метод идентификации подвижной нагрузки и статических напряжений, основанный на методе конечных элементов и вейвлет-анализе. Оба метода используются в работе для решения задачи идентификации режимов колебаний по измеренным реакциям моста (деформациям или ускорениям). Численное имитационное моделирование показало эффективность метода в условиях шума измерений и неровностей дороги. Аналогичной тематике посвящена работа F. Periard, J. Kalker [57], в которой рассмотрено решение во временной области задачи высокочастотного расчета трамвайных рельсов. Разнообразие тематики работ о нелинейной динамике подчеркивает важность исследований в этой области.

Одним из важнейших аспектов исследований в нелинейной динамике является изучение сценариев перехода к хаотическим колебаниям. Механизмы перехода к хаосу в колебаниях балок через бифуркации, под действием периодического нагружения рассмотрены в статье Т. Yagasaki [58]. Т. Yoshimura,

J. Hiño, T. Kamata, N. Ananthanarayana [59] рассмотрели хаотические сложные колебания нелинейной балки переменного поперечного сечения, шарнирно закрепленной, под действием нагрузки. Zhang Jian-wen, Cai Zhong-min, Yang Gui-tong [60] рассмотрели задачу по определению условий возникновения хаотических колебаний от вынуждающей нагрузки бесконечной по длине деформируемой балки, которая лежит на нелинейном упругом основании. Исследованию сценариев перехода в хаос для балок математических моделей Эйлера - Бернулли, С.П. Тимошенко, Шереметьева - Пелеха посвящены работы автора диссертации [61-66]. Исследованию сценариев для других механических структур посвящены работы В.А. Крысько и его учеников, например [67-74].

Для исследования нелинейной динамики используется математический аппарат, состоящий из преобразования Фурье, анализа показателей Ляпунова, фазовых и модальных портретов и автокорреляционной функции. В работе [75] P. J. Holmes, J. Marsden применили метод Мельникова в своих исследованиях хаотических колебаний однослойной балки при действии внешнего нагружения. Qiu Ping и др. [76] используют метод Мельникова для отыскания бифуркаций и хаотических колебаний в задаче о колебаниях пластины на нелинейном упругом основании. При этом предполагается, что сила реакции содержит линейную, изгибную и кубически нелинейную компоненты. Используя метод малых возмущений, были исследованы спектры Ляпуновских показателей, а также хаотические колебания упругой балки под действием периодической внешней силы в работе [77] A. Maewal. В работе Cardona А. [78] изложены некоторые аспекты реализации метода мультигармонического баланса (МНВ) с интенсивным использованием алгоритма быстрого преобразования Фурье на всех стадиях динамического расчета. Устойчивость полученных решений типовых задач нелинейных колебаний проанализирована апостериори по методу Флоке. Точность и эффективность изложенного подхода продемонстрированы на примерах численного расчета показателей колебаний осциллятора Дуффинга, защемленной балки с демпфером сухого трения и пологого висячего троса.

Появление вейвлет-анализа в начале 80-х годов прошлого века как синтеза научных идей из различных областей знаний, таких как математика, радиофизика, теория колебаний, привело к новому направлению исследований в нелинейной динамике. Кроме того, вейвлеты находят применение в обработке и синтезе сигналов, анализе изображений, полученных в различных областях, решении уравнений U. Lepik [79]; сжатии больших объемов информации.

Список публикаций по приложению вейвлет-анализа достаточно обширен; здесь отметим: во-первых, классические работы С.К. Chui, I. Daubechies, А. Grossman, S. Morlet [80-82]; во-вторых, работы, посвященные применению вейвлетов при анализе нелинейной динамики механических структур.

Анализ распространения волн в анизотропных многослойных пластинках был рассмотрен в [83], в [84] для исследования роторной системы, испытывающей вибрации, используются материнские вейвлеты Ньюлэнда. Работа Wongand Chen [85] посвящена применению вейвлета Морле для анализа механической системы, когда частота колебаний гармонических растет во времени. В последнее время появились работы по использованию вейвлетов совместно с нейронными сетями; так, в работе Gholizadeh S. [86] рассмотрена оптимизация конструкции с частотными ограничениями при помощи генетического алгоритма с применением нейронной сети вейвлетной радиальной базисной функции. Коэффициенты вейвлета Лапласа используются в качестве входного вектора для искусственной нейронной сети для анализа динамики в подшипниках качения [87].

Вейвлет-анализ для анализа колебаний механических систем использовался в работах научной группы В.А. Крысько [88-90]. Исследования автора по вейвлет-анализу опубликованы в [63, 64,91-93].

Вейвлеты позволили исследовать явление синхронизации для хаотических колебаний. Явление синхронизации колебаний известно достаточно давно и подробно описано в ряде книг, например [94-96].

Фазовая хаотическая синхронизация, т.е. захват фаз хаотических сигналов, при условии, что амплитуды сигналов остаются несогласованными, наблюдалась в радиотехнике, лазерной физике, медицине и других отраслях. Для получения фазы хаотического сигнала использовались различные подходы, например фаза рассматривалась как угол в полярной системе координат для траекторий хаотического аттрактора. Другим способом является использование преобразования Гильберта для сигнала. Третьим способом было использование сечения Пуанкаре. Однако все эти подходы дают результаты для достаточно простых систем. Для систем с плохо определенной фазой в работах A.A. Короновского и А.Е. Храмова [97, 98] предложено использовать непрерывное вейвлет-преобразование. Использование этого подхода позволяет определить фазу для каждого временного масштаба и проследить изменение фазовой синхронизации во времени. Для распределенных механических систем этот подход развит в работах В.А. Крысько, A.B. Крысько, И.В. Папковой, М.И. Коч, Т.В. Яковлевой и автора диссертационной работы [99-101].

Математические модели, описывающие многослойные конструкции, необходимо разделить на два класса. К первому относятся модели, в которых слои взаимодействуют между собой как одно целое и в процессе деформации конструкций не происходит разделения слоев между собой, такие конструкции принято называть спаянными. Второй класс - это пакеты балок, пластин и оболочек, соединенных через краевые условия, которые будем называть неспаянными. Чаще всего рассматриваются модели первого класса. Обзоры публикаций по исследованию многослойных конструкций, как спаянных, так и неспаянных в виде балок и пластинок, можно найти в статьях Б. Я. Кантора, M.L. Wilkins [102-104]. В статье П.С. Ланда [105] исследовано одно из практических приложений данной теории, имеющее междисциплинарный аспект -моделирование голосовых связок в горле, которые представлены в виде двух пластинчатых конструкций, соединенных пружинами со стенками трубы. Исследование установило, что под действием потока воздуха происходят хаотические колебания прикрепленных пластин. Достаточное число задач, в которых рассмотрено контактное взаимодействие стержней, приведено в книге

P. Wriggers [106], где имеется достаточно полный библиографический материал по контактному взаимодействию различных конструкций.

Исследования нелинейной динамики спаянных и неспаянных конструкций освещены в работах Sogabe Yuji [107], Ray К. [108], а также в работах В.А. Крысько, A.B. Крысько И.В. Папковой, М.И. Коч, B.C. Солдатова и автора диссертационной работы [109-112].

Целью диссертационной работы является построение нового метода решения нелинейных, дифференциальных уравнений в частных производных, сочетающего в себе понижение порядка дифференциального оператора и его линеаризацию, создание программных комплексов на основе этого метода, а также математическое и компьютерное моделирование нелинейной динамики распределенных механических конструкций с учетом геометрической и физической нелинейности и контактного взаимодействия.

Для достижения этой цели были решены следующие задачи:

• Анализ и оценка применимости известных методов линеаризации и понижения порядка нелинейных дифференциальных уравнений в частных производных.

• Разработка теоретических положений и области применимости нового метода понижения порядка и линеаризации дифференциальных уравнений.

• Развитие итерационных вычислительных схем численного решения нелинейных многомерных уравнений статики и динамики.

• Построение на основе созданных итерационных методов и алгоритмов комплексов программ для проведения численного исследования математических моделей нелинейной динамики механических структур в виде балок, пластин и оболочек.

• Решение на основе разработанных подходов и расчетных схем ряда задач для изучения новых явлений в нелинейной динамике механических структур в виде балок, пластин и оболочек.

Научная новизна положений, выносимых на защиту

1. Предложен и обоснован новый эффективный итерационный метод линеаризации и понижения порядка нелинейных дифференциальных уравнений в частных производных, основанный на сведении исходного уравнения к уравнению типа Пуассона на каждом шаге итерационной процедуры. Доказана сходимость предложенных итерационных процедур.

2. На основе предложенного метода разработаны алгоритмы и комплексы программ для различных типов граничных условий и вида областей. Использование предложенного метода позволило ускорить получение численного решения, уменьшить машинную погрешность.

3. Предложена итерационная процедура сведения бигармонического уравнения к уравнению типа Пуассона для областей с криволинейной границей, позволяющая разрешить парадокс Сапонджяна. Доказана сходимость предложенной итерационной процедуры.

4. Построены алгоритмы и разработаны комплексы программ для численного исследования нелинейной динамики балок, математические модели которых построены на основе гипотез Бернулли-Эйлера, С.П. Тимошенко, Шереметьева — Пелеха. Комплексы программ имеют адаптацию к различным видам нелинейности: геометрической, физической и конструктивной. Для программ получены охранные свидетельства.

5. Разработан универсальный программный комплекс для численного исследования и графического представления результатов на основе спектра Фурье, вейвлет-спектров (с различными материнскими вейвлетами), сечения Пуанкаре, показателей Ляпунова и автокорреляционной функции, фазового и модального портретов. Для программ получены охранные свидетельства.

7. Получены как классические, так и модифицированные сценарии перехода гармонических колебаний в хаотические для трех типов математических моделей балок и пластин. Проведены анализ и обобщение сценариев.

8. С помощью вейвлет-анализа впервые изучено явление потери устойчивости балок и пластин при действии поперечной знакопеременной нагрузки, а также замкнутых цилиндрических оболочек при действии радиальной знакопеременной и полосовой знакопеременной нагрузки.

9. На основании эвристического анализа спектра мощности, показателей Ляпунова, автокорреляционной функции построены 56 карт режимов колебаний балок для трех математических моделей, различных типов граничных условий, материалов балки и типов нагрузки, что позволило создать схему диагностики режимов колебаний. Исследовано влияние относительной толщины балки X = а/к на результаты решения задач статики и динамики.

10. Исследовано динамическое поведение многослойных балок на основе гипотез Бернулли - Эйлера, С.П. Тимошенко, Шереметьева - Пелеха, а также пакетов балок, соединенных между собой только через краевые условия, с учетом трех типов нелинейности. Получены новые эффекты: а) эффект фазовой хаотической синхронизации как для упругого материала, так и для физически нелинейного материала балок; б) явление «расслоения» пакета, когда колебания одной из балок после бифуркации происходят вокруг нового положения равновесия; в) явление «полной синхронизации», т.е. фазовой синхронизации и синхронизации сигналов.

Методы исследования

Используются общая методология математического моделирования, математический аппарат начально-краевых задач для нелинейных дифференциальных уравнений математической физики, методы нелинейной динамики, математический аппарат качественной теории дифференциальных уравнений и функционального анализа, численный метод Рунге-Кутты 4-го порядка, методы конечных разностей, граничных и конечных элементов, теория Фурье и вейвлет - анализа

Для создания программных комплексов на основе разработанных алгоритмов были использованы системы программирования С++ и FORTRAN.

Теоретическая и практическая значимость

Разработанный метод решения нелинейных дифференциальных уравнений в частных производных может быть применен при компьютерном анализе широкого класса математических моделей, описываемых дифференциальными уравнениями в частных производных.

Разработанные алгоритмы и программы могут быть использованы для анализа пространственно-временной динамики и частотных характеристик механических распределенных структур в проектной и расчетной практике конструкторских и научно-исследовательских организаций строительного, авиа-, судо- и машиностроительного профилей и приборостроения.

Полученные в диссертации результаты используются в учебном процессе при чтении курсов: «Методы линеаризации, понижения порядка и размерности», «Математические методы в нелинейной динамике», «Математическое моделирование динамических систем», «Проблемы хаоса и нелинейности. Синхронизация колебаний» для специальности «Прикладная математика и информатика», а также для других специальностей Саратовского государственного технического университета имени Гагарина Ю.А. С использованием полученных результатов издано учебное пособие «Математические модели и методы исследований сложных колебаний неклассических распределенных механических систем».

Методы и программы были внедрены и использованы для проектирования некоторых систем дугогасительных камер в ОАО «Контакт», г. Саратов.

Работа выполнялась в рамках госбюджетной темы кафедры «Математика и моделирование» - «Математическое моделирование нелинейных колебаний распределенных систем». Исследования проводились при финансовой поддержке Грантов Министерства образования РФ в рамках целевой программы «Развитие научного потенциала высшей школы на 2009-2011 годы», Российского фонда фундаментальных исследований (проекты №10-08-91153, №10-08-91332, №1208-00569), федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы, финансируемых за счет средств федерального бюджета, выделяемых по направлению расходов «НИОКР», мероприятию 1.2.1 «Проведение научных исследований научными группами под руководством докторов наук» (контракт № П321).

Достоверность результатов, представленных в диссертационной работе, обеспечивается корректностью применения математического аппарата, доказательством теорем сходимости решения по предложенным методам и алгоритмам, непротиворечивостью фундаментальным положениям методологии анализа нелинейной динамики распределенных механических систем, сравнением с признанными зарубежными и отечественными аналогами в области анализа, практическим использованием материалов диссертации и разработанных программных комплексов. Все основные результаты, представленные в диссертации, опубликованы в ведущих российских (из списка ВАК РФ) и зарубежных журналах и прошли рецензирование.

Апробация работы

Основные результаты работы докладывались и обсуждались на 11 международных и 22 всероссийских конференциях, симпозиумах и съездах:

9, 10, 11 Conference on Dynamical Systems: Analytical, Numerical Methods, Stability, Bifurcation and Chaos (Lodz, Poland, 2007, 2009, 2011);

IV Международной конференции по механике неоднородных структур (Тернополь, Украина, 1995);

1-й Международной научно-практической конференции

Дифференциальные уравнения и применения» (Санкт-Петербург, 1996);

XVIII Международной конференции по теории оболочек и пластин (Саратов, 1997);

8, 9, 11, 12 межвузовских и 1, 2, 3, 4, 7 всероссийских конференциях «Математическое моделирование и краевые задачи» (Самара, 1998-2011);

VIII, IX, X Всероссийских съездах по теоретической и прикладной механике (2001,2006, 2011);

Международном конгрессе «Нелинейный динамический анализ - 2007» (Санкт-Петербург, 2007);

Международном семинаре «Актуальные проблемы нелинейной механики оболочек» (Казань, 2008);

International Conference «Chaotic modeling and Simulation» - CHAOS 2009 (Technical University of Creet, Chania, Greece, 2009);

9th SSTA Conference «Shell Structures: Theory and Applications» (Gdansk-Jurata, Poland, 2009);

2 Международной конференции «Проблемы нелинейной динамики деформируемого твердого тела» (Казань, 2009);

VIII Международной научной конференции «Математические проблемы механики неоднородных структур» (Украина, Львов, 2010).

В законченном виде диссертация докладывалась на научных семинарах кафедры «Математика и моделирование» ФГБОУ ВПО «Саратовский государственный технический университет имени Гагарина Ю.А.» под руководством заслуженного деятеля науки и техники РСФСР, д.т.н., профессора В.А. Крысько (Саратов, 2012); на межкафедральном семинаре «Математическое моделирование, численные методы и комплексы программ» ФГБОУ ВПО «Саратовский государственный технический университет имени Гагарина Ю.А.» под руководством заслуженного деятеля науки РФ, д.ф.-м.н., профессора В.Б. Байбурина (Саратов, 2012).

Публикации

Результаты диссертационного исследования опубликованы в 45 печатных работах, из них 1 глава в монографии в издательстве Springer, 20 статей в ведущих иностранных журналах и входящих в Перечень ведущих рецензируемых журналов ВАК РФ, 10 авторских свидетельств о государственной регистрации программы для ЭВМ, а также 1 учебное пособие.

Личный вклад

Все основные результаты, на которых базируется диссертация, получены лично автором. Во всех совместных исследованиях автор принимал участие в выборе направления исследования и формулировке задач. Автору диссертации принадлежит ведущая роль в реализации численных методов и алгоритмов, проведении численных экспериментов и трактовке полученных результатов.

На защиту выносятся следующие основные результаты и положения

1. Обоснован теоретически и реализован в виде комплекса программ метод линеаризации и понижения порядка нелинейных дифференциальных уравнений в частных производных для различных типов граничных условий и вида областей, позволяющий ускорить получение решения, уменьшить вычислительную ошибку, решить парадокс Сапонджяна и упростить использование методов конечных и граничных элементов для сложных нелинейных задач.

2. Доказательство сходимости предложенных итерационных процедур.

3. Алгоритмы и комплекс программ для решения задач нелинейной динамики математических моделей балок на основе гипотез Бернулли, С.П. Тимошенко, Пелеха — Шереметьева с учетом трех типов нелинейности (геометрической, физической и конструктивной), а также алгоритмы и комплексы программ для анализа результатов на основе спектра Фурье, вейвлет-спектра, автокорреляционной функции, фазового и модального портретов.

4. Для различных математических моделей однослойных и многослойных балок впервые получен эффект динамической потери устойчивости при действии поперечной знакопеременной нагрузки. Для многослойных пакетов, связанных через краевые условия, это явление характеризуется «расслоением», т.е. балки перестают воздействовать друг на друга.

5. Построенные и проанализированные 56 карт характеров колебаний для трех математических моделей, ряда граничных условий, типов материала и видов нагрузки дают возможность построить схему диагностики режимов колебаний с возможностью предсказания последствий работы балочных структур. Для построения одной карты решено 9-104 задач и проведено исследование полученных данных с помощью методов нелинейной динамики.

6. Проведенный анализ сценариев перехода к хаотическим колебаниям для трех математических моделей балок позволил найти наиболее общий сценарий, заключающийся в формировании вокруг независимой частоты (или частоты бифуркации удвоения) парных частот, приводящих к зашумлению спектра и появлению хаотических колебаний.

7. Исследования фазовой хаотической синхронизации показали, что зоны синхронизации увеличиваются при смене граничных условий для одной из балок.

Структура и объем работы

Диссертация состоит из введения, шести глав, заключения и списка литературы из 155 наименований. Общий объем диссертации 335 страниц машинописного текста, включающего 126 рисунков, 17 таблиц.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Жигалов, Максим Викторович

Выводы по главе

1. Построены математические модели нелинейной динамики цилиндрических панелей, цилиндрических оболочек и пластин с криволинейной границей. Использовалась гипотеза Кирхгофа — Лява.

2. Для цилиндрических панелей с помощью Фурье и вейвлет анализа, а также анализа старшего показателя Ляпунова проведено исследование влияния на динамические характеристики числа членов ряда, используемых в методе Бубнова. Для одночленного приближения построены сценарии перехода в хаос. Получен классический сценарий Фейгенбаума для различных комбинаций поперечной знакопеременной нагрузки с продольных усилий. Получено новое явление - окно периодичности в вейвлет спектре во времени.

3. Для цилиндрических оболочек с помощью вейвлет-анализа изучена динамическая потеря устойчивости при действии знакопеременной поперечной нагрузки, что позволило качественно подойти к процессу исследования колебаний системы и выяснить характер колебаний в зонах до и после потери устойчивости динамической системы. *

4. С помощью метода установления изучены процессы сходимости динамических задач к статическим.

5. Построен сценарий перехода в хаос для пластинки со сложным очертанием границы под действием локальной поперечной знакопеременной нагрузки.

6. Создан комплекс программ для исследования нелинейной динамики цилиндрических панелей и оболочек, на который получен охранный документ.

317 Заключение

1. Предложен метод линеаризации и понижения порядка нелинейных дифференциальных уравнений в частных производных. Доказаны теоремы сходимости. На основе предложенного метода создан комплекс программ с использованием методов вариационных итераций, конечных разностей, конечных и граничных элементов, для разных типов граничных условий и видов областей.

2. Созданы алгоритмы и комплексы программ для исследования балок математических моделей на основе гипотез Бернулли, С.П. Тимошенко, Пелеха — Шереметьева, в которых учитываются три типа нелинейности (геометрическая, физическая и конструктивная).

3. Создан программный комплекс для численного исследования и графического представления результатов на основе спектра Фурье, вейвлет-спектров с различными материнскими вейвлетами, сечения Пуанкаре, показателей Ляпунова и автокорреляционной функции, фазового и модельного портретов.

4. Получены как классические, так и модифицированные сценарии перехода к хаотическим колебаниям. Проведен анализ и выявлен наиболее общий сценарий.

5. Исследовано влияние относительной толщины балки Л = а/И на результаты решения задач статики и динамики. Выявлено, что увеличение параметра приводит к сходимости результатов для моделей, построенных на основе гипотез Бернулли - Эйлера, С.П. Тимошенко, Шереметьева - Пелеха.

6. Изучено явление потери устойчивости балок и пластин при действии поперечной знакопеременной нагрузки и замкнутых цилиндрических оболочек при действии радиальной знакопеременной и полосовой знакопеременной нагрузок.

7. Исследовано явление синхронизации колебаний для математических моделей контактных задач для пакетов, состоящих из двух и трех балок. I 1

Список сокращений и условных обозначений а - длина структуры высота структуры

- прогиб структуры и(х, - перемещение срединной поверхности вдоль оси ОХ Р{х,у,Г) - функция усилия ух{х,() - поперечный сдвиг 8 - коэффициенты затухания для прогиба Я = я(х>У>0 - поперечная нагрузка Е — модуль Юнга - ускорение свободного падения Р - плотность материала у - объемный вес материала балки т а ~ коэффициент линейного температурного расширения в - разность температур у - коэффициент Пуассона 12ЕНЗ 12(1-~ ЦилиндРическая жесткость р - частота вынуждающей силы

Яо - амплитуда вынуждающей силы и к — зазор между слоями кх,ку- кривизны оболочки р, (0-фаза сигнала

- обозначение поперечной силы 6,3 - модуль сдвига

Для Главы 5:

И/ — толщина балки толщина балки в центре ^(х,/) - прогиб балок и1 (х, /) - перемещения в срединной линии Е1 - модуль Юнга материала балок G0^ - модуль сдвига р1 - удельная плотность материала У1 - коэффициент Пуассона еп - интенсивность деформаций <7п - интенсивность напряжений ез1 — интенсивность деформаций текучести <7з1 - интенсивность напряжений текучести К1 (х, г) - объемный модуль упругости £01 - объемная деформация у1 - коэффициент поперечной деформации

320

Список литературы диссертационного исследования доктор физико-математических наук Жигалов, Максим Викторович, 2013 год

1. Kantorovich L. V., Krylov V.l., V.l. 1958, Approximate Methods of Higher Analysis (The Netherlands, P. Noordhaff).

2. Санкин Ю.Н. Об одном численном методе в нелинейной теории тонкостенных упругих оболочек / Санкин Ю.Н. и др. // Тр. Ульяновск, политехи, ин-та. 1972. Т. 8. № 2. С. 191-202.

3. Мальгин В.Н. Алгоритмы решения задач прочности, устойчивости и колебаний оболочек вращения, основанные на уравнениях типа С.П. Тимошенко / В.Н. Мальгин // Методы решения задач упругости и пластичности: сб. ст. Горький, 1973. Вып. 7. С. 137-142.

4. Валишвили Н.В. Об одном алгоритме решения нелинейных краевых задач / Н.В. Валишвили // Прикл. матем. и механика. 1968. Т. 32. № 6. С. 1089-1092.

5. Вержбицкий В.М. О свободных от обращения вложенных итерациях Ньютона / В.М. Вержбицкий // Краевые задачи: сб. ст. Пермь, 1979. С. 8384.

6. Машков A.B. К вопросу использования метода Ньютона-Канторовича при решении задач нелинейной теории упругости эластомеров / A.B. Машков // Методы решения задач упругости и пластичности: сб. ст. Горький, 1990. Вып. 15. С. 141-145.

7. Петров В.В. К расчету пологих оболочек при конечных прогибах / В.В. Петров // Научные доклады высшей школы. Строительство. 1959. № 1. С. 27-35.

8. Григолюк Э.И. Проблемы нелинейного деформирования: Метод продолжения решения по параметру в нелинейных задачах механики твердого деформируемого тела / Э.И. Григолюк, В.И. Шалашилин. М.: Наука. Гл. ред. физ.-мат. лит., 1988. 232 с.

9. Шалашилин В.И. Метод продолжения решения по параметру и наилучшая параметризация в прикладной математике и механике / В.И. Шалашилин, Е.Б. Кузнецов. М.: Эдиториал УРСС, 1999. 222 с.

10. Григоренко Я.М. Численное решение задач о напряженном состоянии гибких некруговых цилиндрических оболочек / Я.М. Григоренко, А.Т. Василенко, H.H. Крюков // Прикладная механика. 1984. T. XX. № 3. С. 4044.

11. Temple G. The general theory of relaxation methods applied to linear systems / G. Temple // Proc. Roy. Soc., Ser. A, 1939.

12. Канторович JI.B. Об одном эффективном методе решения экстремальных задач для квадратичных функционалов (о градиентном методе наискорейшего спуска) / Л.В. Канторович // Докл. АН СССР. 1945. Т. 48. № 7. С. 345-357.

13. Вайнберг Д.В.Методы численного анализа в теории упругости / Д.В. Вайнберг, А.Л. Синявский // Тр 2-го Всесоюз. съезда по теор. и прикл. механике. 1964. М., 1964. С. 83-94.

14. Вайнберг Д.В. Метод спуска и программирование задач строительной механики пластин и оболочек / Д.В. Вайнберг, Е.С. Дихтярюк, А.Л. Синявский // ЭЦВМ в строит механике: сб. ст. М.: Стройиздат, 1966. С. 465-470.

15. Симеонов C.B. Некоторые методы решения нелинейных задач механики деформируемого твердого тела / C.B. Симеонов // Прикл. матем. и механ. 1964. Т. 28. № 2. С. 418-429.

16. Пухлий В.А. К расчету сопряженных оболочек переменной жесткости / В.А. Пухлий // Прикл. механика. 1989. Т. 25. № 11. С. 31-37.

17. Thurston G.A. A Numerical Solution of the Nonlinear Equations for Axisymmetric Bending of Shallow Spherical Shells / G.A. Thurston // Trans ASME. 1961. E. 28.№4. P. 557-562.

18. Mescall J. Numerical solutions of nonlinear equations for shells of revolution / J. Mescall // AIAA Journal. 1966. V. 4. № 11. P. 2041-2043.

19. Thurston G.A. Continuation of Newton's Method Through Bifurcation Points / G.A. Thurston // Trans ASME. 1969. E. 36. № 3. P. 425-430.

20. Grüters J. Iterative Lösung von Lastspanningsproblemen in anisotropen Körpern / J. Grüters // Z. angew. Math, und Mech. 1974. M. 54. № 4. T79-T80.

21. Давиденко Д.Ф. Об одном новом методе численного решения систем нелинейных уравнений / Д.Ф. Давиденко // Докл. АН СССР. 1953. Т. 88. № 4. С. 601-602.

22. Давиденко Д.Ф. О приближенно решении систем нелинейных уравнений / Д.Ф. Давиденко // Украинский мат. ж. 1953. Т. 5. № 2. С. 196-206.

23. Леньков В.Ф. Изгиб нелинейно упругого клина в квадратичной теории упругости / В.Ф. Леньков, В.М. Собин // Технол. машиностроения: сб. ст. Тула, 1969. Вып. 14. С. 113-116.

24. Карасик М.И. Об одном шаговом методе решения нелинейных уравнений теории упругости / М.И. Карасик, В.И. Шалашилин // Прочность конструкций: сб. ст. Уфа, 1976. № 1. С. 87-92.

25. Солянова О.Н. Изгиб квадратной пластинки с отверстием, внешний контур которой защемлен / О.Н. Солянова // Механика деформ. сред: сб. ст. Саратов: Изд-во Сарат. ун-та, 1973. Вып. 1. С. 5-10.

26. Уздалев А.И. Свободные колебания двухсвязанных пластин / А.И. Уздалев, Л.Н. Нагибин // Механика дефор. сред: сб. ст. Саратов: Изд-во Сарат. ун-та, 1973. Вып. 1. С. 136-142.

27. Самарский А. А. Введение в теорию разностных схем / А. А. Самарский. М.: Наука, 1971. 552 с.

28. Слободянский М. Г. Пространственные задачи теории упругости для призматических тел / М. Г. Слободянский // Ученые записки Московского государственного университета. Вып. 39. Механика. Г940

29. Hartree D. A. Method tor the Numerical or Mechanical Solution of Certain Types of Partial Differential Equations / D. A. Hartree // Proc. Roy. Soc. 1937. Vol. 161. P. 353.

30. Nowacki W. Az ortotrop lemezek elemeltenek nehany problemaja / W. Nowacki //Magyar tud. Akad.musz. tud. Oszt. Kozl. 1953. V. II. № 1-2. P. 123-141.

31. Giangreco E. Sur le calcul des dalles a rigidite variable // Non-Homogeity in Elasticity and Plasticity. London-New York-Paris-LosAngeles, Pergamon. Press, 1959. P. 63-76.

32. Han L.S. The buckling and deflection of isosceles-triangular plates / L.S. Han // Trans ASME. 1960. E. V. 27. №1. P. 207-208.

33. Пшеничнов Г.И. Метод декомпозиции решения уравнений и краевые задачи Г.И. Пшеничнов // Докл. АН СССР. 1985. Т. 282. №4. С.792-794.

34. Миндолин Ю.И., Уздалев А.И. О преобразовании уравнения в частных производных четного порядка к системе уравнений Лапласа и Пуассона / Ю.И. Миндолин, А.И. Уздалев // Известия РАН. МТТ. №2. 1993. С. 43-45.

35. Розин JI.A. Метод расчленения в теории оболочек / JI.A. Розин // ПММ. 1961. №25. С. 921-926.

36. Conway H.D. Application of the point-matching method to shallow-spherical-shell theory / H.D. Conway, A.W. Leissa // Trans ASME. 1962. E. V. 29. №4. P. 745-747.

37. Улитин B.B. Итерационные алгоритмы решения краевых задач механики на ЭВМ / В.В. Улитин. СПб: Изд-во СПб ун-та, 1991. 232 с.

38. Nishihara Т. Deflection of laterally loaded square plates under various edge conditions / Nishihara Toshio, Tanaka Kichinosuke // Mech. Fac. Engng, Kyoto Univ. 1953. V.15. № 4. P. 197-212.

39. Seth B.R. Bending of T-plate / B.R. Seth // Proc. 1st Congr. Theor. and APP1. Mech. 1955. Kharagpur, Indian Inst. Technol., s.a. p. 87-90.

40. Duffm RJ. On exterior boundary value problems in linear elasticity / R.J. Duffin, W. Noll // Arch. Ration. Mech. And Analysis, 1958. V. 2. № 2. P. 191196.

41. Roller B. Ortotrop lemezek alakvaltozasainak numerikas szamitasa / B. Roller //Epitoipari es Korlek. Miisz. Egyet, tud. Kozl. 1961. V. 7. № 3. p. 93-104.

42. Stussi F. Zur Prandtlschen Membranalogie der Torsion / F. Stiissi // Z. angew. Math und Phys. 1958. V. 96. № 5-6. P. 661-667.

43. Почтман Ю.М. Расчет шарнирно опертых пластин методом электрического моделирования / Ю.М. Почтман, В.Д. Шайкевич // Строит, механ. и расчет сооруж. 1965. № 2. С. 24-26.

44. Chang I-Dee, Finn R. On the solutions of a class of equations occurring in continuum mechanics, with application to the stokes paradox / I-Dee Chang, R. Finn // Arch. Mech. and Analysis. 1961. V. 7. №5. P. 388-40.

45. Жигалов M.B."Методы линеаризации дифференциальных уравнений механики деформированного твердого тела /°М.В. Жигалов, Т.В. Бабенкова //Вестник СГТУ. 2009. №2. Вып. 1. С. 24-38.

46. Жигалов М.В."Методы понижения порядка дифференциальных уравнений механики деформированного твердого тела (обзор) /"М.В.

47. Жигалов // Вестник Саратовского государственного технического университета. 2006. №1(11). Вып. 2. С. 13-32.

48. Kreuzer Е. Nonlinear dynamics in ocean engineering / E. Kreuzer // ICTAM 2004: 21st International Congress of Theoretical and Applied Mechanics, Warsaw, Aug.15-21, 2004: Abstracts and CD-ROM Proceedings. Warszawa, 2004. C. 5.

49. Nonlinear dynamics of high-speed milling-analyses, numerics, and experiments / Stepan Gabor, Szalai Robert, Mann Brian P., Bayly Philip V., Insperger Tamas, Gradisek Janez, Govekar Edvard // Trans. ASME. J. Vibr. and Acoust. 2005. 127. N 2. C. 197-203.

50. Law S. S. Moving load and prestress identification using wavelet-based method / S. S. Law, S. Q. Wu, Z. Y. Shi // Trans. ASME. J. Appl. Mech. № 2, 2008, T.75. C. 021014/1-021014/7.

51. Periard F. Time-domain solution of the high frequency dynamic behaviour ofttitram rails / F. Periard, J. Kalker //19 Int. Congr. Theor. and Appl. Mech., Kyoto, Aug. 25-31, 1996: Abstr. Kyoto, 1996. C. 689.

52. Takao Yagasaki. Bifurcations and chaos in quasi-periodically forced beam. Theory, simulation and experiment / Yagasaki Takao // J. Sound and Vibr. 1995. 183. №1. P. 1-31.

53. Random vibration of non-linear beam subjected to a moving load: a finite element method analysis / T. Yoshimura, J. Hino, T. Kamata, N. Ananthanarayana // J. Sound and Vibr. 1988. 12. № 2. P. 317-329.

54. Zhigalov M.V. Analysis of regular and chaotic dynamics of the Euler-Bernoulli using finite difference and finite element methods / M.V. Zhigalov, A. V. Krysko, O. A. Saltykova, J. Awrejcewicz // Modeling, Simulation andar

55. Жигалов M.B. Нелинейная динамика балок Эйлера-Бернулли. (Математическая модель, сценарии перехода колебаний из гармонических в хаотические) / М.В. Жигалов, В. А. Крысько, Е.Ю. Крылова // Известия вузов. Строительство. 2011. №2. С.15-21.

56. Жигалов М.В. Особенности нелинейных колебаний балок С.П. Тимошенко / М.В. Жигалов, А. В. Крысько, В.В. Солдатов // Изв. вузов. Строительство. 2009. № 5. С. 25-35.

57. Жигалов М.В. Нелинейная динамика балок Эйлера-Бернулли и типа Тимошенко / М.В. Жигалов, В.А. Крысько, O.A. Салтыкова // Известия вузов. Машиностроение. 2008. №6. С. 7-27.

58. Awrejcewicz, J. Feigenbaum Scenario Exhibited by Thin Plate Dynamics / J. Awrejcewicz, V. A. Krysko. // Nonlinear Dynamics. 2001. № 24. P. 373-398.

59. Крысько В.А. О сценариях перехода колебаний из гармонических в хаотические для гибких пластинок и их управление / В.А. Крысько, И.В. Папкова, Е.В. Салий // Вестник Саратовского государственного технического университета. 2009. № 4 (42).

60. Кузнецова Э.С. О влиянии температурного поля на сложные колебания замкнутых цилиндрических баллонов / A.B. Крысько, Я. Аврейцевич, Э.С. Кузнецова // Вестник Саратовского государственного технического университета. 2008. № 1(31). Вып. 2. С. 71-85.

61. Кузнецова Э.С. Исследование хаотических колебаний прямоугольных пластинок при действии поперечной знакопеременной нагрузки в температурном поле / В.А. Крысько, Э.С. Кузнецова, Н.Е. Савельева // Известия вузов. Машиностроение. 2006. №1. С. 3-9.

62. Крысько В.А. Управление хаотическими колебаниями гибких сферических оболочек / В.А. Крысько, И.В. Кравцова // Изв. РАН МТТ. 2005. № 1.С. 10-20.

63. Крысько В.А. Устойчивость замкнутых цилиндрических оболочек при неравномерном внешнем давлении / В.А. Крысько, Н.Е. Савельева // Известия вузов. Строительство. 2005. № 5. С. 3-14.

64. Крысько В.А. Стохастические колебания конических оболочек переменной толщины / В.А. Крысько, Т.В. Щекатурова // Известия вузов. Машиностроение. 2004. № 5. С. 3-13.

65. Крысько В. А. Колебания балки Эйлера-Бернулли на упругом основании по моделям Винклера и В.З. Власова / В. А. Крысько, А. С. Десятова. // Изв. вузов. Строительство. 2004. №7. С. 101-115.

66. Holmes, P. J. A partial differential equation with innately many periodic orbits: chaotic oscillations of a forced beam / P. J. Holmes, J. Marsden // Archives for Rational Mechanics and Analysis. 1981. 76. 135-166.

67. Bifiircation and chaos of the circular plates on the nonlinear elastic foundation / Qiu Ping, Wang Xin-zhi, Yeh Kai-yuan // Appl. Math, and Mech. Engl. Ed. -2003. 24, ~N 8. - C. 880-885.

68. Maewal A. Chaos in a harmonically excited elastic beam / A. Maewal // ASME Journal of Applied Mechanics. 1986. 53. 625-631.

69. Cardona A. Fast Fourier nonlinear vibration analysis / A. Cardona, A. Lerusse, M. Geradin // Comput. Mech. 1998. 22, 2. C. 128-142.

70. Lepik U. Application of the Haar wavelet transform to solving integral and differential equations / U. Lepik // Proc. Estonian Acad. Sci. Phys. Math. 2007. 56., 1.28-46.

71. Chui С. К. A cardinal spline approach to wavelets / С. K. Chui, and, J. Z. Wang // Proc. Amer. Math. Soc. 1991. 113. 785-793.

72. Grossman A. Decomposition of Hardy functions into square separable wavelets of constant shape / A. Grossman, S. Morlet // SI AM J. Math. Anal. 1984. Vol. 15. №4. P. 723.

73. Daubechies I. (). Painless nonorthogonal expansions / I. Daubechies, , A. Grossmann and Y. Meyer// J. Math. Phys. 1986. 27. 1271-1283.

74. Jeong H. Analysis of plane wave propagation in anisotropic laminates using a wavelet transform / H. Jeong // NDT&E Int. 2001. 34. 185-190.

75. Zheng Jibing. Application of wavelet transform to bifurcations and chaos study Jibing Zheng, Hangshan Gao, and Yinchao Guo // Appl. Math. Mech. (Shanghai, China). 1998. 19. 593-599.

76. Wong L. A. Nonlinear and chaotic behaviour of structural system investigated by wavelet transform technique / L. A. Wong, and, J.C. Chen // Int. J. NonLinear Mech. 2001. 36. 221-235.

77. Gholizadeh S. Structural optimization with frequency constraints by genetic algorithm using wavelet radial basis function neural network / S. Gholizadeh, E. Salajegheh, P. Torkzadeh // J. Sound and Vibr. 2008. N 1-2. T. 312. C. 316331.

78. Application of the Laplace-wavelet combined with ANN for rolling bearing fault diagnosis / F. Al-Raheem Khalid, Asok Roy, K. P. Ramachandran, D. K. Harrison, Steven Grainger // Trans. ASME. J. Vibr. and Acoust2008. N 5. T. 130. C. 051007/1-051007/9

79. Крысько B.A. Анализ нелинейных хаотических колебаний пологих оболочек вращения с помощью вейвлет преобразования / В.А. Крысько, И.В. Папкова, В.В. Солдатов // Известия РАН. Механика твердого тела. 2010. №1. С. 107-117.

80. Awrejcewicz J. On the wavelet transform application to a study of chaotic vibrations of the infinite length flexible panels driven longitudinally /

81. J. Awrejcewicz, A. Krysko, V. Soldatov // Int. J. of Bifurcation and Chaos in Applied Sciences and Engineering. 2009. Vol. 19. Issue 10. P. 3347-3371.

82. Аврейцевич Я. О памяти нелинейных дифференциальных систем в теории пластин / Я. Аврейцевич, Е. Ю. Крылова, И. В. Папкова, Т. В. Яковлева, В. А. Крысько // Вестник Нижегородского университета им. Лобачевского. 011. № 4. Ч. 2. С. 21-23.

83. Жигалов М.В. Анализ хаотических колебаний распределенных систем в виде балок Эйлера-Бернулли с помощью вейвлет преобразования / М.В. Жигалов, А. В. Крысько, В.В. Солдатов // Изв вузов. Авиационная техника. 2009. № 4. С. 21-24.

84. Жигалов М.В. Вейвлет-анализ колебаний замкнутых цилиндрических оболочек / М.В. Жигалов, В. А. Крысько, В.В. Солдатов, Э.С. Кузнецова // Вестник Саратовского государственного технического университета. 2009. №4. С. 24-30.

85. Жигалов М.В. О выборе типа вейвлета при изучении нелинейных колебаний балок с учетом поперечных сдвигов / М.В. Жигалов, А. В. Крысько, В.В. Солдатов, М.Н. Подтуркин // Вестник Саратовского государственного технического университета. 2009. №3. С. 14-22.

86. Пиковский А. Синхронизация: Фундаментальное нелинейное явление /

87. A. Пиковский, М. Розенблюм, Ю. Курте. М.: Техносфера, 2003. 496 с.

88. Блейхман И.И. Синхронизация в природе и технике / И.И. Блейхман. М.: Наука, 1981.352 с.

89. Синхронизация регулярных, хаотических и стохастических колебаний /

90. B. С. Анищенко, В. В. Астахов, Т. Е. Вадивасова, Г. И. Стрелкова. М.; Ижевск: Ин-т компьютерных исследований,. 2008. 144 с.

91. Короновский A.A. Анализ фазовой хаотической синхронизации с помощью непрерывного вейвлетного преобразования / A.A. Короновский, А.Е. Храмов // Письма в ЖТФ. 2004. Т. 30. Вып. 14. С. 2936.

92. Кантор Б. Я. Контактные задачи нелинейной теории оболочек вращения / Б. Я. Кантор. Киев: Наук, думка, 1990. 136 с.

93. Wilkins M.L. Calculation of elasto-plasic flow / M.L. Wilkins // Methods of Computational Physics / B. Alder, S. Fernbach and M. Rotenberg, editors. V.3. Academic Press, Boston, New York (1964).

94. Ланда П.С. Хаотические колебания в модели голосовых связок / П.С. Ланда // Известия вузов. Прикладная нелинейная динамика. 1998. Т. 6. №4. С. 57-67.

95. Wriggers P. Computational contact mechanics / P. Wriggers. 2nd ed. Springer-Verlag Berlin Heidelberg, 2006.

96. Newland D.E. Introduction to Random Vibrations, Spectral and Wavelet Analysis / D.E. Newland. Longman, N. Y., 1993.

97. Yuji S. Transient vibration analysis of elastically connected Timoshenko double-beam systems / Sogabe Yuji, Yokoyama Takashi, Tsuzuki Masayuki, Wu Zhiqiang // Nihon kikai gakkai ronbunshu. С N Trans. Jap. Soc. Mech. Eng. C. 1996. 62, 594. C. 429-437.

98. Ray K. The parametric instability of partially covered sandwich beams / K. Ray, R. C. Kar // J. Sound and Vibr. 1996. 197. 2. C. 137-152.

99. Жигалов M.B. Нелинейная динамика неоднородных структур / М.В. Жигалов, В. А. Крысько, А. В. Крысько // Математические проблемы механики неоднородных структур: тр. VIII Междунар. науч. конф. Украина, Львов, 14-17 сентября 2010 г. С. 14-21.

100. Жигалов М.В. Управление сложными колебаниями нелинейных многослойных балок / М.В. Жигалов, А. В. Крысько, О.А. Салтыкова // Известия вузов. Авиационная техника. 2008. № 3. С. 10-13.

101. Математические модели и методы исследования сложных колебаний неклассических распределенных механических систем: учеб. пособие / А.В. Крысько, М.В. Жигалов. Саратов: Сарат. гос. техн. ун-т, 2008.230 с.

102. Вольмир А.С. Нелинейная динамика пластин и оболочек / А.С. Вольмир. М.: Наука, 1972. 432 с.

103. Бочкарев В.В. Об одном подходе к решению геометрически нелинейных задач теории пластинок /В.В. Бочкарев, В.А. Крысько // Изв. вузов. Строительство и архитектура. 1981. №10. С. 30-34.

104. Ворович И.И. О существовании решений в нелинейной теории оболочек / И.И. Ворович / Изв. АН СССР. Сер. матем. 1955. Т. 19. № 4. С. 173-186.

105. Пшеничный Б.Н. Численные методы в экстремальных задачах / Б.Н. Пшеничный, Ю.М. Данилин. М.: Наука, 1975. 319 с.

106. Лионе Ж.-Л. Некоторые методы решения нелинейных краевых задач / Ж.-Л. Лионе. М.: Мир, 1972. 587 с.

107. Бублик Б.Н. Численное решение задач пластин и оболочек / Б.Н. Бублик. Киев: Изд-во Киев, ун-та,1969. 148 с.

108. Дородницын A.A. Асимптотическое решение уравнения Ван-дер-Поля / A.A. Дородницын // ПММ. Т. XI. 1947. С. 313-328.

109. Пальцев Б.В. О разложении решений задачи Дирихле и смешанной задачи для бигармонического уравнения в ряд по решениям распадающихся задач/ Б.В. Пальцев // ЖВМиМФ. 1966. Т. 6. № 1. С. 4351.

110. Шереметьев М.П. К построению уточненной теории пластин / М.П. Шереметьев, Б.Л. Пелех // Инженерный журнал. 1964. Т.4. Вып. 3. С. 3441.

111. Крысько А. В. Математические модели нелинейных распределенных систем в виде пластинчатых конструкций : дис. . д-ра физ.-мат. наук : 05.13.18, 01.02.04 / А. В. Крысько. Саратов, 2003. 343 с. РГБ ОД, 71:051/77

112. Ramberg W. Normal pressure tests of rectangular plates / W. Ramberg, A.E. Pherson, S. Levy // NACAT. № 849. 1942. 60 p.

113. Крысько В.А. Нелинейная статика и динамика неоднородных оболочек / В.А. Крысько. Саратов: Изд-во Сарат. ун-та, 1976. 214 с.

114. Баглай Р.Д. К обработке двумерных сигналов на ЭВМ / Р.Д. Баглай, К.К. Смирнов // ЖВМиМФ. 1974. № 1. С. 241-248.

115. Сьярле Ф. Метод конечных элементов для эллиптических задач / Ф. Сьярле. М.: Мир, 1980. 512 с.

116. Уздалев А.И. Температурные напряжения в пластинках ограниченных двухсвязанным контуром / А.И. Уздалев. Саратов: Изд-во СГУ, 1975. 176 с.

117. Лехницкий С.Г. Анизотропные пластинки / С.Г. Лехницкий. М.: Гостехиздат, 1957. 600 с.

118. Тимошенко С.П., Гудьер Дж. Теория упругости: пер. с англ. / С.П. Тимошенко, Дж. Гудьер. М.: Наука, 1975- 576 .

119. Стренг Г., Фикс Дж. Теория метода конечных элементов / Г. Стренг, Дж. Фикс. М.: Мир, 1977. 226 с.

120. Feigenbaum M.J. Quantitative universality for a class of nonlinear transformations / MJ. Feigenbaum // J. Stat. Phys. 1978. V. 1. № 1. P. 25-52.

121. Feigenbaum M.J. Quasiperiodicity in dissipative systems: A renormalization group analysis / M.J. Feigenbaum, L.P. Kadanoff, S.J. Shenker // Physica. 1982. V.D.5 P. 370.

122. Амбарцумян С.А. Теория анизотропных пластин / С.А. Амбарцумян. М.: Наука, 1967.

123. Reissner Е. On transverse vibration of thin shallow shells / E. Reissner // Quarterly of Appl. Vath. 13. №2 (1955). P. 169-170.

124. Timoshenko S.P. On the correction for shear of differential equation for transverse vibration of prismatic bar / S.P. Timoshenko // Philosophical Magazine. 41. 3 6 (1921). P. 744-746.

125. Ильюшин А.А. Пластичность. Основы общей математической теории пластичности / А.А. Ильюшин. М.: Изд-во АН СССР, 1963. 272 с.

126. Феодосьев В.И. Об одном способе решения нелинейных задач устойчивости деформируемых систем / В.И. Феодосьев // ПММ. 1963. Т. 27. Вып. 2. С. 265-274.

127. Биргер И.А. Некоторые общие методы решения задач теории пластичности / И.А. Биргер // Прикладная математика и механика. 1951. т. 15. Вып 6. С. 765-770.

128. Крысько В.А. Нелинейная статика и динамика неоднородных оболочек / В.А. Крысько. Саратов: СГУ, 1976. 216 с.

129. Бернштейн С.А. Сопротивление материалов: учеб. для машиностроит. специальностей вузов / С.А. Бернштейн. М.: Высш. шк., 1961. 464 с.

130. Osipov G.V., Pikovsky A.S., Rosenblum M.G., Kurths J. // Phys.rev. Lett 1997. V. 55. P. 2353.

131. Pikovsky A.S. Synchronization: a Universal Concept in Nonlinear Sciences / A.S. Pikovsky, M.G. Rosenblum, J. Kurths. Cambridge University Press, 2001.

132. Короновский А.А. Анализ хаотической фазовой синхронизации динамических систем с плохо определенной фазой / А.А. Короновский, А.Е. Храмов // Радиотехника и электроника. 2005. Т. 50. №7. С. 1-9.

133. Anishchenko V.S. Synchronization of Self-Oscillations and Noise-Induced Oscillations / V.S. Anishchenko, T.E.Vadivasova // Journal of Communications Technology and Electronics 2002. V. 47. № 2. P. 117-148.

134. Nonlinear Dynamics of Chaotic and Stochastic Systems. Tutorial and Modern Developments / V.S. Anishchenko, V. Astakhov, A. Neiman et al. Heidelberg: Springer Verlag, 2007. 461 p.

135. L'achaux J.P. et al. // Int. J. Bifurcation and Chaos. V. 10. P. 24-29.

136. Debnath L. Wavelet transforms and their applications / L. Debnath. Boston. Birkhâuser Verlag, 2002. 565 p.

137. Wave propagation and sampling theory. Part I: Complex signal land scattering in multilayer media / J. Morlet, G. Arens, E. Fourgeau, D. Giard // Geophysics. 1982. V. 47. P. 203-221.

138. Ободан Н.И. Устойчивость оболочек при неосемметричной деформации / JI.B. Андреев, Н.И. Ободан, А.Г. Лебедев. М.: Наука, 1988. 208 с

139. Баженов В.Г. Упругопластическое деформирование составных обол очечных конструкций при импульсных воздействиях. / В.Г. Баженов, А.П. Шинкаренко // Проблемы прочности. 1981. № 3. С. 25-29.5. />

140. Сагерлинд JI. Применение метода конечных элементов / Л. Сагерлинд. М.: Мир,: 1979. 393 с.

141. Krysko V.A. Chaotic vibrations of closed cylindrical shells in a temperature field / V.A. Krysko, J. Awrejcewicz, E.S. Kusnetsova, A.V. Krysko // International Journal Shock and Vibration. Part2. 2008. № 15. P.335-343.

142. Shian A.C. Dynamic Buckling of conical shells with Imperfection / A.C. Shian, T.T. Soong, D.S. Roth // AIAA Journal, 1974. Vol. 12. № 6. P. 24-30.

143. Кантор Б.Я. Нелинейные задачи теории пологих оболочек / Б.Я. Кантор. Киев: Наук, думка, 1971. 136 с.