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

• Крысько, Антон Вадимович
• доктор физико-математических наук
• 2003

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

ВВЕДЕНИЕ (краткий исторический обзор по теме работы).

1.0 Методы сочленения элементов обол очечных конструкций.

2.0 Пространственные конструкции, элементы которых подчиняются гипотезам Кирхгофа-Лява, типа Тимошенко.

3. ° Комбинированные математические модели контактных задач теории пластин и оболочек.

4.° Контактные задачи для неспаянных, многослойных пластин и оболочек.

5.° Сложные колебания балок, пластин и оболочек.

6.° Вейвлет - преобразования.

ГЛАВА I. МАТЕхМАТИЧЕСКИЕ МОДЕЛИ ПРОСТРАНСТВЕННЫХ

КОНТАКТНЫХ ЗАДАЧ ТЕОРИИ ПЛАСТИН.

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

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

§ 1.3 Вариационный принцип и краевые задачи для пространственной пластинчатой конструкции на базе гипотез Кирхгофа.

§ 1.4 Вариационный принцип и краевые задачи для пространственной пластинчатой конструкции на базе гипотез Тимошенко.

§1.5 Математические модели пространственных пластинчатых конструкций в рамках трехмерной теории упругости и теории Кирхгофа.

§ 1.6 Математические модели пространственных пластинчатых конструкций в рамках трехмерной теории упругости и гипотез Тимошенко.

§ 1.7 Математические модели пространственных пластинчатых конструкций в рамках гипотез Кирхгофа и гипотез Тимошенко.

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

ГЛАВА И. НЕЛИНЕЙНЫЕ КОНТАКТНЫЕ СТАТИЧЕСКИЕ И

ДИНАМИЧЕСКИЕ ЗАДАЧИ ТЕОРИИ МНОГОСЛОЙНЫХ НЕСПАЯННЫХ ПЛАСТИН.

§2.1 Методы типа Канторовича-Власова для статических задач теории пластин.

§ 2.2 Контактные статические задачи для конструкций, состоящих из физически нелинейных неспаянных пластин.

§ 2.3 Построение коллокационных методов на основе хМетодов типа Канторовича-Власова.

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

§ 2.5 Численный анализ НДС статики двухслойных пластин МВИ и МКР.

§ 2.6 Динамические контактные задач теории многослойных неспаянных пластин.

§ 2.7 Математические модели и итерационные алгоритмы решения контактных задач для пространственных конструкций из неспаянных пластин в рамках гипотез типа Тимошенко и комбинированной модели.

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

ГЛАВА III. СЛОЖНЫЕ КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ 7 ОРТОТРОПНЫХ ГИБКИХ ПЛАСТИНОК

ПРИ ДЕЙСТВИИ ПРОДОЛЬНЫХ НАГРУЗОК

ПОСТОЯННЫХ ВО ВРЕМЕНИ.

§ 3.1 Устойчивость гибких ортотропных пластинок при действии продольных нагрузок.

§ 3.2 Последовательность бифуркаций Андронова-Хопфа в гибких ортотропных пластинках при действии продольных постоянных во времени нагрузок.

§ 3.3 Жесткая потеря устойчивости гибких ортотропных пластинок.

§ 3.4 Влияние начальных условий на последовательность бифуркации Андронова - Хопфа, хаос и жесткую потерю устойчивости гибких ортотропных пластинок.

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

ГЛАВА IV. ПАРАМЕТРИЧЕСКИЕ КОЛЕБАНИЯ ГИБКИХ

ИЗОТРОПНЫХ ПЛАСТИН.

§4.1 Методы сведения распределенных систем к сосредоточенным

MKP-0(h2)u0(h4).

§ 4.2 Расчет параметрических колебаний пластин

МКР-0(/г2) иМКР-0(/г4).

§ 4.3 Исследование влияния на сложные колебания пластинок некоторых параметров (начальная амплитуда возмущений, коэффициент демпфирования).

§ 4.4 Симметричные диссипативные параметрические колебания гибких пластин.

§ 4.5 Несимметричные параметрические колебания гибких пластин.

§ 4.6 Образование стоячих волн при диссипативных параметрических колебаниях гибких пластинок.

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

ГЛАВА V. МЕТОД БУБНОВА ПРИ ИССЛЕДОВАНИИ СЛОЖНЫХ

ПАРАМЕТРИЧЕСКИХ КОЛЕБАНИЙ ГИБКИХ ПЛАСТИНОК.

§5.1 Анализ систем со многими степенями свободы.

§ 5.2. Постановка задачи. Метод Бубнова-Галеркина.

§5.3 Численный эксперимент.

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

ГЛАВА VI. ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЯ В ТЕОРИИ

ПАРАМЕТРИЧЕСКИХ КОЛЕБАНИЙ ГИБКИХ ПЛАСТИНОК.

§ 6.1. Основы теории.

§ 6.2. Базисные функции вейвлет-преобразования.

§ 6.3. Свойства и возможности вейвлет - преобразования.

§ 6.4 Применение вейвлет - анализа к исследованию сложных параметрических колебаний гибких пластинок.

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

ВЫВОДЫ ПО ДИССЕРТАЦИИ.

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

Вводные замечания. Запросы практики уже в течении многих веков предопределили необходимость построения математических моделей основных элементов современных конструкций: стержней, балок, пластин, оболочек и различных их сочленений. История данной проблемы восходит своими корнями к работам JI. Эйлера, Я. Бернулли, С. Жермен, Ж.Л. Лагранжа, Л. М. А. К. Навье, А.Л. Коши, Г. Ламе, А.Ж К. Б. де Сен-Венана, Б.П.Э. Клай-перона, Г.Р. Кирхгофа, Р.Р.А. Клебша, И. Баушингера, Д.У. Релейя, Ф.Энгессера, А. Фёппля, И.Г. Бубнова, С.П. Тимошенко, Т. Кармана, В.В. Новожилова, Х.М. Муштари, К.З. Галимова, Н.А. Алфутова, С.А. Абарцумяна, А.Н. Андреева, В.Г. Баженова, А.Б. Богдановича, В.В. Болотина, А.С. Вольми-ра, А.Т. Василенко, В.В.Васильева, И.И. Воровича, Э.И. Григолюка, М.С. Танеевой, Я.М. Григоренко, Б.Я. Кантора, В.Г. Карнаухова, Ю.Г. Коноплева, М.С. Корнишина, В.И. Королева, В.А. Крысько, Н.Д. Кузнецова, Н.Ф. Морозова, Ю.В. Немировского, Ю.Н. Новичкова, И.Ф.Образцова, В.Н. Паймуши-на, Б.Л. Пелеха, В.В. Петрова, В.В. Пикуля, В.Г. Пискунова, А.П. Прусакова, Б.Е. Победри, И.Н. Преображенского, А.О. Рассказова, А.В. Саченкова, Э.Г. Терегулова, П.Е. Товстика, К. Ф. Черныха, Ю.В. Чеботаревского, Э. Рейс-нера и др.

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

Кирхгофа-Лява либо типа Тимошенко. Построенные математические модели позволят записать разрешающие уравнения и условия сопряжения (контакта). Это требует создания новых вариационных принципов.

Как отмечается в работе Э.И.Григолюка и В.М.Толкачева [53], "контактные задачи теории пластин и оболочек возникают при рассмотрении взаимодействия пластин и оболочек с жесткими и упругими телами (штампами), ребрами жесткости, при взаимном контакте пластин и оболочек". С позиции контактных задач могут рассматриваться слоистые пластины и оболочки, если вводить реакции взаимодействия между слоями и определять их из условий стыковки слоев. Механика контактного взаимодействия упругих и пластических тел представляет собой важную и интенсивно развивающуюся область современной механики деформируемого твердого тела.

Соответствующая контактная проблема для тонкостенных элементов (стержней, пластин, оболочек) ставилась и решалась, в основном, в последние три десятилетия. Полученные результаты в этом направлении систематизировали в монографиях В.М.Александрова [6] (1983 г.), Э.И.Григолюка и В.М.Толкачева [53] (1980г.), В.И.Моссаковского и В.С.Гундрамовича [131] (1978 г.), Б.Д. Пелеха и М.А.Сухорольского [158] (1980 г.), Б.Л.Пелеха, А.В.Максимуха и И.М.Коровачуха [157] (1988 г.), БЛ.Кантора [73] (1990г.), А.В. Крысько [86] (2000 г.).

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

7°. Методы сочленения элементов оболочечных конструкций. Тонкостенные оболочечные конструкции, состоящие из нескольких слоев и подкрепленные в различных направлениях ребрами, как правило, имеют сложную форму, представляя собой набор различных оболочек - цилиндрических, сферических, конических и более сложной конфигурации: с утолщениями, сварными швами и каналами охлаждения. Очень часто такие элементы выполнены из различных материалов (изотропные, анизотропные), что, естественно, накладывает свои ограничения на те или иные теории. В общем случае следует иметь законченную математическую теорию сочленения элементов конструкций, каждый элемент которой .описывается теми или иными кинематическими гипотезами (математическими моделями). В настоящее время общей единой теории нет.

Интенсивное развитие космической и авиационной техники привело к созданию ряда автоматизированных систем расчета осесимметричных оболочечных конструкций [78]. Оболочечные конструкции соединяются непосредственно между собой или с помощью шпангоутов. Использовалась теория гибких оболочек, кинематическая модель Кирхгофа-Лява, физическая нелинейность учитывалась по деформационной теории пластичности. Для шпангоутов использовалась теория тонкостенных стержней В.З.Власова [40]. Уравнения оболочек, колец получены в проекциях на оси, связанные с недеформированной системой координат кольца. Это позволило достаточно просто сформулировать задачу прочности' для оболочечных конструкций, состоящих из набора оболочек различного типа, которые состыкованы друг к другу по параллельным кругам непосредственно или с помощью упругих шпангоутов. Дискретизация -МКР, линеаризованная задача решалась методом ортогональной прогонки.

Получил распространение подход сопряжения оболочечных конструкций, основанный на схеме, когда подкрепляющие кольца рассматриваются как короткие оболочки вращения переменной толщины [52] (1990 г.).

Практический расчет сложных конструкций, составленных из разнородных элементов: стержней, пластин, оболочек, трехмерных тел - встречает большие трудности. Преодоление этих трудностей связано с идеей расчленения исходной задачи на ряд более простых. Первоначальный, традиционный подход состоял в' идеализации, огрублении связей между элементами. Расчет системы сводился к последовательному расчету элементов без учета их взаимного влияния. В настоящее время решение идет по двум направлениям: первый - метод суперэлементов (МСЭ), второй - итерационные методы решения контактных задач теории оболочек. Первый подход (МСЭ) развит в работах ряда авторов [44, 157, 217]. Второй подход - алгоритмы, реализующие итерационные методы (ИМ), как правило, значительно более просты. К их недостаткам относятся трудности исследования сходимости и, как следствие, - затраты машинного времени. Этот подход развивается в работах [2, 175].

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

Однако существует большое число нерешенных проблем. Еще в 1970 году в проблемной статье [258] была отмечена проблема МКЭ в расчетах оболочек. Это, прежде всего, несогласованность КЭ [45, 64, 134] (нарушение неразрывности первой производной от прогиба). В вариационных постановках для решения задач деформирования оболочек классическая теория оказалась малопригодной. Она была создана для сведения трехмерной задачи к двухмерной, но это достигалось ценой повышения порядка производной в искомых функционалах. Следует отметить, что функционалы теории упругости содержат лишь первые производные от искомых функций. Самой распространенной теорией оболочек, функционалы которой содержат только первые производные от перемещений, является теория оболочек и пластин типа Тимошенко (или Рейс-нера). В работе Р.Б. Рикардса [162] (1988 г.) на базе теории Тимошенко построены КЭ оболочек вращения и оболочек произвольной формы, лишенные недостатков классической теории.

Анализ имеющейся литературы дает возможность сделать вывод, что в настоящее время не построены КЭ для сочлененных оболочечных конструкций» когда отдельные элементы описываются различными математическими моделями (классическая теория, теория типа Тимошенко, трехмерная теория, теория стержней Кирхгофа-Клебша и др.). Построение такой теории, возможно только исходя из вариационных принципов. Вариационные формулировки создают предпосылки для развития прямых аналитических и численных методов. Созрела необходимость построения комбинированных, математических моделей контактных задач теории пластин и оболочек.

2°. Пространственные конструкции, элементы которых подчиняются гипотезам Кирхгофа - Ляеа, типа Тимошенко. Как было отмечено в предыдущем пункте, численный анализ позволил рассмотреть сложные задачи сопряжения оболочечных конструкций. В литературе имеются решения в рамках математических моделей теории оболочек, использующих гипотезы Кирхгофа-Лява.

Особенно интенсивно сочлененные конструкции стали рассчитываться за последние годы, в основном, рассмотрены соединения замкнутых цилиндрических оболочек с круглыми пластинками (работа Н.И. Фоминой [174]), соединение двух цилиндрических оболочек под прямым углом, (А.Л.Козак, Донг Хну Кун, Нгуен Ши Чан [80]), с конусом [67], с многослойным днищем [172]. В основном, вышеназванные конструкции рассчитывались методом конечных элементов.

В работе [239] Pierre Ladeveze дается обзор современного состояния вычислительной строительной механики и перспективы ее дальнейшего развития. Проанализировано более 180 публикаций в западных журналах. Выделяется, как перспективное, создание композитных сочлененных структур. Указываются преимущества, и недостатки. Группа исследователей из Казанского государственного университета создала комплекс программ по расчету тонкостенных конструкций как трехмерных тел [51] (1989 г.). Анализ НДС проводится совокупностью экспериментальных методов поляризационно -оптический, электротензометрия, голографическая интерферометрия) и МКЭ. Использовался искривленный конечный элемент оболочки, построенный на основе трехмерных уравнений теории упругости. Данный элемент, характеризующийся введением специальных кинематических гипотез типа Тимошенко, позволяет получить достоверные результаты, как для тонких оболочек, так и для конструкций средней толщины.

Сопряжению цилиндрических оболочек под прямым углом, обладающих малой сдвиговой жесткостью, посвящена работа И.Г.Стрельченко [170]. Отмечается, что при рассмотрении сопрягаемых оболочек в рамках кинематической модели Кирхгофа-Лява требуется восемь условий статического сопряжения, а для модели типа Тимошенко число условий статического сопряжения увеличивается.

В заключение данного параграфа следует выделить целый цикл работ В.Н.Паймушина и его учеников по созданию теории сопряжения элементов конструкций. Основополагающими работами здесь являются [141, 142]. Согласно этим работам тело разделяется на N элементов, через S^nk) обозначены поверхности сопряжения п- го элемента с k-ым элементом и, в соответствии с контактной формулировкой задачи, на Sfa) = S^ введены в рассмотрение неизвестные реактивные усилия взаимодействия . Используется третий закон Ньютона: q^nk) = - . Это приводит к формулировке удвоенного числа статических граничных условий, записанных для каждого элемента. Получен обобщенный вариационный принцип Лагранжа [3,156, 171].

Ряд работ В.Н.Паймушина, Ю.Я.Петрушенко [140, 143, 152] посвящен обобщению идей [141, 142] на более сложные задачи; применению обобщенного вариационного принципа Рейснера к построению линеаризованных задач механики; слоистых оболочек сложной геометрии; устойчивости и динамики составных оболочек; многослойных оболочек со слоями переменной толщины.

Анализ литературы показывает, что практически создана теория сочленения элементов конструкций, описываемых классической теорией пластин и оболочек. Интенсивно развивается теория сочленения оболочек, описываемых моделью типа Тимошенко (Рейснера). Но до настоящего времени она еще далека от совершенства.

3°. Комбинированные математические модели контактных задач теории пластин и оболочек. При исследовании НДС составных конструкций, включающих в себя тонкостенные (изотропные, трансверсально-изотропные, анизотропные) и массивные элементы, а также узлы их сопряжения, возникает необходимость в применении комбинированных математических моделей. Литературы данного направления существенно меньше, чем по направлениям, сформулированным в 2° данного обзора. Вопросам построения и применения таких моделей посвящен ряд работ [143, 269]. В этих работах для описания НДС в одних частях деформированного тела предлагается привлекать уравнение теории упругости, а в других - уравнения классической теории оболочек. В трех работах Я.Г.Савула [60, 164, 165] построены комбинированные математические модели трехмерной теории упругости и теории оболочек типа Тимошенко. Рассмотрены вариационные постановки краевых задач комбинированной модели. Приводится численный пример расчета рамы МКЭ. Построены простейшие КЭ сопряжений.

Развитию конечного элемента Варани (сопряжение тел, описываемых классической теорией оболочек и трехмерной теорией упругости) посвящена работа: Акимкин С.А. , Никишнов Г.П. Сопряжение трехмерных элементов и элементов тонкостенной оболочки для расчета напряженного состояния конструкции сложной формы //Деформация и разрушение материалов и элементов конструкций ядерных энергетических установок М, 1986. СЗ-9.

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

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

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

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

Наиболее проста линейная постановка для цилиндрических оболочек разной длины, установленных с натягом. Без учета обжатия, т. е. когда в решение входят сосредоточенные поперечные силы на границе зоны контакта, задача изучена авторами работ [25, 26, 54, 55], где решены дифференциальные либо интегральные уравнения. Обжатие по модели Винклера введено в работе [24], по модели упругого цилиндра и слоя — в [81, 82]. В двух последних работах контактное давление становится бесконечным на границах зон контакта. С помощью теории Тимошенко эта задача исследована в [144]. Решение такой же задачи получено [23] представлением контактного давления в виде суммы произведений неизвестных коэффициентов на заданные функции, ортонормиро-ванные на участке контакта. Коэффициенты вычисляются методом наименьших средних квадратов из кинематического условия контакта, граница зоны контакта уточняется итеративным путем. Этот подход позволяет существенно упростить расчеты, поскольку в нем не требуется решать дифференциальные или интегральные уравнения относительно контактного давления, результаты же полностью совпадают с данными [24, 25]. Такой же метод применен в работах [27, 28, 29] для анализа НДС двухслойного сильфона с податливым промежуточным кольцом.

Взаимодействие двух соосных цилиндрических оболочек разной длины с зазором между ними при нагружении внутренним давлением оболочки меньшего радиуса изучено в [178, 230]. Авторы работы [230] сопрягают аналитические решения уравнений равновесия оболочек в зоне контакта и вне ее, получают систему уравнений относительно произвольных постоянных, находят осевую координату границы зоны контакта, решая систему трансцендентных уравнений. Сочетание вариацинно-разностного метода с методом штрафной функции применено в [178]. Обжатие в обеих работах не учтено, использованы теории Кирхгофа — Лява, Тимошенко, Рейсснера.

Взаимодействие двух сферических оболочек рассмотрено по теории Кирхгофа — Лява в [170], причем обнаружено, что контактная реакция — распределенная по окружности сосредоточенная сила. Решение этой же задачи основано на теории Тимошенко в [152]. На границе зоны контакта получено контактное давление, равное нулю, хотя оно должно было бы принять конечное значение. Контакт двух круглых пластин, установленных с зазором при нагружении одной из них, изучен в [14] с использованием теории Жермен - Лагран-жа - Кирхгофа. На границе зоны контакта обнаружены сосредоточенные сила и момент. Теория Рейсснера позволила получить конечное значение контактного давления' qk на границе [209]. Задача о контакте между двумя прямоугольным пластинам решена вариационно-разностным методом в [59]. Перечисленные исследования исходят из линейной теории оболочек.

Геометрически нелинейная теория оболочек применена в работе [139] для изучения МКЭ контакта между слоями гофрированных мембран. Условия контакта здесь представлены специальными физически нелинейными элементами между узлами слоев, входящими в соприкосновение.

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

В случае если нормальные перемещения и напряжения на соприкасающихся поверхностях слоев совпадают, а касательные напряжения равны нулю, приходим к задаче об оболочке с проскальзывающими без трения слоями. Зоны контакта .при этом известны, что существенно упрощает задачу. В указанной постановке решены задачи статики слоистых цилиндрических [37] и сферических [142] оболочек. Метод последовательных приближений, основанный на «принципе поочередной непрерывности», в соответствии, с которым краевые задачи для слоев решаются на каждой итерации независимо, применен в [153, 174] для изучения слоистых цилиндров и цилиндрических оболочек. Задачи о контактном взаимодействии между тонкими пластинками и оболочками сложны, поскольку п£>и их решении приходится одновременно определять НДС и зоны контакта двух и более пластин и оболочек, в общем случае различной формы.

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

Контакт двух круглых пластин, установленных с зазором при нагружении одной из них, изучен Ю.П.Артюхиным [11] (1973 г.) с использованием теории Кирхгофа. На границе зоны контакта обнаружены сосредоточенные сила и момент. Использование теории типа Тимошенко позволило получить конечное значение контактного давления на границе [25] (1973 г.) решение осесиммет-ричной задачи.

Задача о контакте между двумя прямоугольными пластинками решена вариационно-разностным методом [59]. Данные решения выполнены в геометрически и физически линейной постановке.

Дискретный подход для пластин со слоями Тимошенко реализован в [149] с помощью предложенного авторами матричного метода, приводящего задачу к системе интегральных уравнений относительно контактного давления в неизвестных априори зонах. Здесь учтена возможность появления разрывов областей сопряжения слоев. Наиболее полно разработана дискретная теория в [149]. ■

5°. Сложные колебания балок, пластин и оболочек. Гибкие пластинки, подверженные интенсивному периодическому воздействию, представляют собой сложную динамическую систему, в которой в зависимости от шевеления (изменения) параметрами воздействия реализуются принципиально различные режимы колебаний. Процесс колебаний может сопровождаться большим разнообразием физических явлений, к которым относятся возникновение сложных резонансов; срыв" колебательного режима, приводящий к режиму изменения пространственно-временного состояния: стоячим либо бегущим волнам и др. явлениям^ например, появлению потери устойчивости по симметричной и несимметричной формам.

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

Широкомасштабные и планомерные исследования взаимосвязи хаоса и порядка ведутся относительно недавно. Они показали, что в поведении сложных нелинейных систем со многими степенями свободы при определенных условиях могут возникать регулярные пространственные и временные структуры [174], названные И. Пригожиным [50, 135, 157, 158] диссипативными. Наряду с этим возможна и обратная картина: из упорядоченного движения рождается хаос. Большой прогресс в концепции временного динамического хаоса достигнут в таких современных отраслях знаний, как физике плазмы [58], гидромеханике [124], электронике и радиофизике [261], в теории управления [158]. В интернете в разделе Bibliography on Stochastic Resonance американским исследователем Gammaitioni сделана подборка статей по указанному наименованию с 1980 по 1998 годы по следующим реферируемым журналам Physics Letters A; Physics Review Е; Optics Letters; Neutral Comput; Journal of Physical chemistig и др. Отмечается рост публикаций с 1996 г., в среднем в год публиковалось до 70 статей. На наш взгляд, это доказывает огромный интерес ученных к данной проблеме. Существенно скромнее успехи достигнуты в биологии [63] и химии [175]. В механике деформируемого твердого тела достижения в этом направлении еще скромнее. В нелинейных задачах теории пластин имеются лишь единичные публикации [244, 245]. Это связано с тем, что здесь приходится рассматривать пространственно-временную систему с бесконечным числом степеней свободы.

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

Явление турбулентности известно уже сотни лет, и многие исследователи пытались его изучить. В этом направлении опубликовано большое количество работ, но решить данную проблему почти не удаётся из-за чрезвычайной её сложности. Ниже мы приведём несколько моделей, описывающих переход механической системы к турбулентности. В настоящий момент известны следующие модели: Ландау [123], Рюэля-Тэкенса [262], Фейгенбаума [211], Помо-Манвиля [251]. Эти модели схематически приведены в табл. 0.1. Очень коротко опишем её.

Первоначально картина возникновения турбулентности, представленная Ландау [123], была основана на представлении об иерархии неустойчивости. При увеличении некоторого параметра, например, числа Рейнольдса, или числа Рэлея нелинейные колебания жидкости теряют устойчивость, и появляются всё новые и новые независимые частоты движения а> х,со 2,со При этом должно наблюдаться квазипериодическое движение с одной, двумя, тремя и т.д. основными частотами. Таким образом, мы приходим к последовательности бифуркаций, Андронова-Хопфа, т.е. к движению по поверхности некоторого тора возрастающей размерности. Движение выглядит всё более и более сложным, однако непрерывный спектр и хаотическое движение возникают лишь при бесконечном числе бифуркаций. Модель Ландау предложена для бесконечномерной модели. Модели же Рюэля-Тэкенса, Фейгенбаума и Помо-Манвиля связаны с конечномерными моделями.

Рюэль и Такенс [262] предложили другой механизм возникновения турбулентности, согласно которому сначала происходит две последовательные бифуркации Хопфа, как и в модели Ландау, затем нелинейность разрушает трёхчастотное движение и образуется «странный» аттрактор (табл. 0.1).

Некоторые экспериментальные данные, по видимому, подтверждают модель Рюэля-Тэкенса. Так, в спектрах мощности появляется сначала одна, затем вторая и, возможно, третья независимые частоты. На пороге появления третьей частоты внезапно возникает широкополосный шум, который свидетельствует о переходе к хаотическому движению. Экспериментально исследовались как вихри Тейлора в жидкости между вращающимися цилиндрами [212], так и конвекция Рэлея-Бенара [183]. Здесь обнаружена внутренняя синхронизация между частотами со х и со 2. В другом эксперименте по течению Куэтта [227] наблюдались, по крайней мере, четыре независимые частоты. Это указывает на то, что переход к турбулентности происходит не всегда после двух бифуркаций Анд-ронова-Хопфа.

Третья модель перехода к турбулентности, предложенная Фейгенбаумом, связана с последовательностью бифуркаций удвоения периода [211]. Переход начинается с бифуркации Андронова-Хопфа из устойчивого фокуса в предельный цикл с частотой со х. При дальнейшем увеличении параметра происходят последовательные бифуркации удвоения, приводящие к периодическому движению с частотами со Vt со v, со и и т.д. Эта последовательность сходится при

2 /4 '8 г некотором критическом значении параметра, при котором возникает странный аттрактор (см. табл. 0.1).

Таблица 0.1

Модель Механизмы

1. Ландау [25] Непод- виж. точка Бифур кация Хопфа Перио-дич. трае-ктор. Бифуркация Хопфа Квазипериод. траектория (две осн. частоты-) Бифуркация Хопфа Квазипериод. траектория (три осн. частоты) Бифур- кац Хопфа Турбул. движения

2. Рюэдя-Тэкенса [171 -«-«- -«-«- -«-«- Стран. Аттрактор

3.Фейгенбаума [26] -«-«(период Т) Период траект. (период 2T) Бифуркация удвоения Перио- дич. траект. (период 4T) Бифуркация удвоен. Странный атрак-тор

4. Помо-Манвиля [27] -«-«- ^ Обратная тан-ген, бифуркация Хаотич. движен. с перемежаемостью

Модель Фейгенбаума хорошо подтверждается численными экспериментами на простых моделях, таких, как аттрактор Рёслера, отображение Хенона, уравнение Дюффинга и др.

Четвёртый механизм возникновения турбулентности, лежащий в основе модели Помо-Манвиля [251] и связанный с переходом к хаотическому движению с перемежаемостью. В этой модели при увеличении некоторого параметра периодическая траектория непосредственно превращается в хаотическую с перемежаемостью в результате обратной тангенциальной бифуркации.

Численное моделирование квадратичного отображения подтверждает такое поведение. Оно существует и для модели Лоренца в некотором интервале параметров, в конвекции Рэлея-Бекара и так называемой химической турбулентности [259].

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

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

Общей трактовке указанных вопросов, посвящена, например, монография Капитаника [231], ориентированная на инженеров - практиков. Обзор результатов проведенных в 1980-е годы исследований свойств переходных процессов в колебаниях нелинейных систем опубликовал Капитаник в [234]. Построены отображения Пуанкаре и типы фазовых траекторий с различными вариантами неустойчивости. Хаотические эффекты сопоставляются с расчетными характеристиками нелинейности системы. В его же работе [232] проведен анализ условий перехода к хаотическому поведению в автономной самовозбуждающейся системе под действием периодического и случайного внешнего возмущения. Определены фазовые траектории с одной и двумя петлями и условия перехода к хаотическим фазовым траекториям. Хан, Жанг и Янг [222] рассматривали хаотические вынужденные колебания динамической системы второго порядка с квадратичной и кубической нелинейностями. Определялись условия возникновения хаоса по отображениям Пуанкаре, фазовым портретам и временным рядам.

В 1979 году Холмс [225] подробно исследовал хаотические движения слегка выпученного стержня, подвергающегося боковому синусоидальному возмущению, как системы с одной степенью свободы. Мун [253] установил, что гармонически вынужденное движение выгнутого стержня отчетливо демонстрирует хаотический характер: анализ сечений Пуанкаре показывает сложную, но устойчивую структуру. В работе предложен критерий установления порогового для возникновения хаоса значения амплитуды вынужденных колебаний как функции частоты. Тенг и Доуэлл [266] определили пороговое значение нагрузки, вызывающие хаотическое перемещение выгнутого стального консольного стержня, на основе анализа временного ряда. Подцар, Мун и Мухерджи [257] усовершенствовали расчетную модель шарнирно закрепленного упруго-пластического стержня с учетом геометрической и физической нелинейности. Развита численная процедура определения хаотических движений стержня при периодическом нагружении.

Хаотическим колебаниям деформируемых систем посвящен доклад Хана-гуда и Ашлани [223]. В частности, в нем для стержневых систем выбраны разные пути перехода к хаосу в колебательных процессах. А.С. Беломытцев и В.И. Карабан [19] для систем, моделирующих крутильные колебания силовых передач, обнаружили область странного аттрактора, возникающую в результате серии бифуркаций удвоения периода решения. Изучены квазипериодические колебания периодически возбуждаемой системы и многорежимные; периодические колебания. Хан, Ху и Янг [221] рассматривали стойку с жестко защемленными концами, к одному из которых динамически приложена сжимающая сила. Выявлены условия бифуркации форм движения: боковое выпучивание или осевые колебания прямолинейного стержня.

В 1985 г. Саймондс и Ю [264] теоретически обнаружили и экспериментально подтвердили, что при определенных условиях в колебаниях упруго пластической балки с закрепленными концами, нагруженной коротким поперечным импульсом, наблюдается аномально высокая чувствительность к изменению ведущих параметров. Детально изучая этот эффект, Ли, Саймондс и Бо-рино [241] нашли, что если в начальный период нагружения балки движение неизбежно носит хаотический характер, то и переходящий хаос регистрируется как во временном ряде, так и на фазовом портрете и спектрограмме. Установлена также экспоненциальная природа повышения чувствительности процесса к изменениям параметров. В серии статей, опубликованных Ю. Лепиком, эти наблюдения были продолжены и расширены. Например, в [245] анализируются нелинейные поперечные колебания выпученной при продольном сжатии балки под действием гармонического возбуждения. Дана оценка пороговой величины поперечной динамической нагрузки, при которой колебания упруго - пластической балки не переходят в хаотический режим. При этом оказалось, что для балок указанного типа установившиеся хаотические колебания в случае гармонического возбуждения гораздо более обычны, чем при импульсном нагружении.

Бифуркационные механизмы перехода к хаосу в сложных колебаниях балок под действием квазипериодического нагружения проанализированы в статье Ягасаки [277]. Абхьянкар, Холл и Ханагуд [183] предложили новый метод численного решения дифференциальных уравнений в частных производных, на основе которого провели компьютерные эксперименты для выявления закономерностей перехода к хаосу в колебаниях балочных структур. Иосимура, Хино, Камата и Анантанараяна [279] описали хаотические колебания шарнирно закрепленной балки переменного поперечного сечения под действием транспортной нагрузки. Жанг, Цай и Янг [280] решили задачу по выявлению условий возникновения хаотических вынужденных колебаний бесконечной деформируемой банки, лежащей на упругом нелинейном основании. Сложным колебаниям консервативных и диссипативных механических систем в виде многослойного пакета неспаянных балок посвящена работа В.А. Крысько, В.В. Боч-карева и Т.А. Бочкаревой [107]. Оказалось, что хаотические колебания в собственном смысле здесь не возникают, хотя и наблюдаются элементы переходного сценария Фейгенбаума.

Баженов В.А., Дехтерюк Е.С. и Петрина Ю. С. [18] численно исследовали бифуркации установившихся режимов вынужденных колебаний пластин и оболочек под действием периодических во времени нагрузок. Отмечен переход от регулярных (периодических и квазипериодических) колебаний к хаотическим. В работе [246] Ю. Лепик пытался выяснить возможность хаотических реакций в осесимметричных колебаниях упруго - пластических цилиндрических оболочек. В большинстве проведенных компьютерных экспериментов установившиеся колебания имели регулярный характер. Хан, Ху и Янг [220] провели анализ нелинейных колебаний упругой цилиндрической оболочки вращения, и нашли критические условия возникновения хаотического движения. Маэстрелло, Френди и Браун [248] изучали нелинейные колебания типовой панели фюзеляжа самолета. Для выбранного диапазона частот найдены линейные, квазилинейные при удвоении периода колебаний и хаотические динамические реакции панели при увеличении уровня акустического движения звуковой нагрузки.

Возникает необходимость изучения задач со многими степенями свободы. Хаотические движения для вынужденно изогнутой балки исследовались Тангом и Давелом [266] . В [174] изучались упругие хаотические колебания в изогнутых балках. Балка предполагалась начально-сжатой осевой нагрузкой и затем зафиксированной в сжатом состоянии. Колебания вызывались поперечной периодической нагрузкой или колебанием опор. Уравнения движения решались двумя методами: Бубнова - Галеркина и конечных разностей. Хаотический режим балки был установлен по фазовому изображению и по каскаду удвоений периода.

Исследованию нелинейных упруго-пластических колебаний изогнутых балок был посвящен ряд статей Симондса и его коллег [242, 264] . В 1985 году они рассмотрели следующий вопрос: балка с закрепленными концами подвергалась короткому интенсивному импульсу поперечной нагрузки, который производит пластический прогиб. Поскольку концы балки закреплены, то мембранные усилия должны также быть приняты в расчет. Решая уравнения движения, они нашли, что остаточные прогибы могут быть в направлении противоположном нагрузке. Это явление они назвали «аномальным» или «противоестественным» поведением балки.

В случае «противоестественного» поведения возникает вопрос: может ли быть реакция балки хаотической. Ли и Симондс [242] показали, что в случае модели с одной степенью свободы движение полностью определено, и там не может быть никаких хаотических колебаний. Что касается балки с двумя степенями свободы, то из исследования Симондса вытекает, что хаотические эффекты могут иметь место [241].

Динамическая реакция сплошных упруго-пластических балок при коротких импульсных „нагружениях была рассмотрена Лепиком [245, 246, 247]. В этих статьях уравнения движения интегрировались конечно-разностным методом или методом Бубнова, было установлено «противоестественное» поведение и хаотическая реакция балок.

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

Осесимметричные колебания упругих и упруго-пластических оболочек исследовались в нескольких статьях. Большинство из них посвящено динамическому выпучиванию под осевыми нагрузками или осевому удару массы по оболочке. В работе Флоранса и Гудьера [213] рассматривалось пластичное динамическое выпучивание цилиндрических оболочек под осевыми нагрузками.

Решение упрощается предположением, что оно разделяется на доминирующее движение и движение возмущенного типа. Применяется уравнение Прандтля -Рейсса; делается предположение, что в течение реакции не происходит упругое разгружение. В [182] представлены результаты для динамического поступательного выпучивания круговых трубок под осевым толчком. Для того чтобы исследовать динамический продольный изгиб цилиндрических оболочек, было разработано несколько численных схем (например, статья Морино и соавторов [254]). В статье [241], для интегрирования уравнений движения был предложен метод Бубнова, метод применим для произвольного числа степеней свободы; осевые силы инерции также принимаются во внимание. Составляющие уравнения основаны на теории Прандтля-Рейсса, рассматриваются упругое разгружение и обратное пластичное нагружение. Даются некоторые численные примеры и обсуждаются возможности хаотических колебаний.

Изучению явлений хаоса при колебаниях балок с различными краевыми условиями посвящен ряд работ. В [80, 255, 266] исследовались нелинейные поперечные- или продольные колебания балок, подверженных периодическому поперечному или продольному воздействию. Хаотические явления при колебаниях балки [96] были обусловлены нелинейностью краевых условий. Во всех перечисленных работах рассматривались однослойные балки.

В [18] приводятся результаты численного исследования нелинейных установившихся колебаний шарнирно-опертой квадратной в плане пластины при действии равномерно распределенного нормального давления, интенсивность, которого меняется во времени по гармоническому закону. Методом конечных элементов краевая задача сведена к системе обыкновенных дифференциальных уравнений. Установлено, что при переходе к хаотическим колебаниям реализуется серия бифуркаций удвоения периода.

Сценарий перехода диссипативных систем при воздействии на них гармонических нагрузок в различных отраслях современной науки, таких как радиофизика, радиоэлектроника, гидромеханика описан достаточно подробно в работах П. Берже, Н. Помо, К. Видаля [20]; А.С. Дмитриева, В.Я. Кислова [58]; Ю.И. Неймарка и П.С. Ланда [134]; А. Лихтенберга и К. Либермана [126], Ф. Муна [143] и др. Задачи же теории пластин существенно отличаются от задач, приведенных в указанных книгах, т.к. здесь мы по сути дела имеем многомерные системы: две пространственные координаты и время, и приходится рассматривать колебание системы с бесконечным числом степеней свободы во времени. 'В настоящей работе мы попытаемся описать сценарий перехода к пространственно-временному хаосу и особенности некоторых систем в развитом хаосе на примере гибких пластинок. Изучению сложных колебаний изотропных пластин при действии поперечных периодических нагрузок посвящена кандидатская диссертация Салий Е.В. [166], выполненная под руководством профессора Крысько В. А.

6°. Вейелет - преобразования. Термин «вейвлет» (дословный перевод -маленькая волна) появился сравнительно недавно - его ввели Гросс*ман и Мор-ле [219] в середине 80-ых годов в связи с анализом свойств сейсмических и акустических сигналов. В настоящее время семейство анализаторов, названных вейвлетами, начинает широко применяться в задачах распознавания образов; при обработке и синтезе сигналов (например, речевых); при анализе изображений различной природы (изображение радужной оболочки глаза, рентгенограмма почки, спутниковые изображения облаков или поверхности планеты, снимок минерала и т.п.); для изучения свойств турбулентных полей; для решения уравнений; для свертки (упаковки) больших объемов информации и во многих других случаях.

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

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

Далеко не случайно многие исследователи называют вейвлет-анализ «математическим микроскопом» - название прекрасно отражает замечательное свойство метода сохранять хорошее разрешение на разных масштабах. Способность этого «микроскопа» обнаружить внутреннюю структуру исследуемого объекта и изучить его локальные скейлинговые свойства уже продемонстрирована на многих примерах.

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

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

В первом разделе проводится аналогия между рядами Фурье и разложением в ряды по вейвлетам, определяется вейвлет-преобразование, описываются его свойства и некоторые следствия из них; материал основан, главным образом, на сборниках и монографиях [273, 274, 275] и прекрасных работах Ингрид Добечи [207, 208] и Мари Фарж [210]. Во втором разделе приводятся результаты применения вейвлет-преобразований к исследованию хаотических колебаний гибких пластин при продольных знакопеременных нагрузках.

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

Проведенный информационный анализ показал актуальность и необходимость исследования следующих проблем:

1. Построение математических моделей пространственных контактных задач теории пластин и оболочек в рамках однотипных и комбинированных моделей.

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

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

4. Не изучены статические и динамические задачи теории многослойных неспаянных пластин с учетом физической нелинейности и разномодульности материала конструкции.

5. Изучение хаоса в теории гибких изотропных и ортотропных пластинок, его структуры, сценариев развития, а также поиску новых периодических решений, их бифуркациям, в том числе жесткой потери устойчивости изучены недостаточно. Качественное и количественное сопоставление решений полученных различными методами (конечных разностей с аппроксимацией О (h2) и

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

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

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

Научная новизна работы заключается в следующем:

1) разработана единая методика, на базе метода множителей Лагранжа, получения условий сопряжения пластинчатых и оболочечных элементов в пространственных конструкциях, с учетом конечной площади зоны контакта, при этом для описания условий равновесия элементов конструкций могут быть применены произвольные гипотезы;

2) впервые дается новая классификация эволюционных задач теории неспаянных многослойных пластин - консервативно-диссипативные системы — когда один элемент рассматривают без учета диссипации, а другой - с ее учетом;

3) для стационарных задач разработаны новые итерационные алгоритмы с использованием метода вариационных итераций и доказана их сходимость для контактных задач теории многослойных неспаянных пластин с учетом обобщенных гипотез Винклера, созданы комплексы программ их расчета;

4) созданные комплексы программ для ПВМ позволили получить новые количественные данные о НДС, контактирующих пластинках в рамках обобщенных гипотез Винклера с учетом физической нелинейности и разномодуль-ности материала;

5) дается дальнейшее развитие метода коллокаций, предложены новые его модификации на основе восьми вариантов метода Канторовича-Власова (метод решетчатой коллокации, граничной коллокации и смешанной коллокаций);

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

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

8) впервые для анализа параметрических колебаний гибких пластин применяется вейвлет - анализ;

9) с помощью процедуры Бубнова на примере параметрических колебаний гибких пластин выявлены сценарии перехода колебаний из гармонических в хаотические в зависимости от числа мод.

Достоверность результатов. С одной стороны обеспечивается их полным соответствиям основным уравнениям термодинамики для конечных объемов сплошной среды, вариационными принципами, а с другой гарантируется рядом теорем и проведением численного эксперимента различными методами метод конечных разностей 0(h2), О (/г4), Бубнова в высших приближениях, Канторовича - Власова и его восьмью модификациями , быстрого преобразования Фурье, вейвлет - анализа, методами компьютерной динамики.

Практическая ценность полученных результатов заключается в использовании, при проектировании и оценке качества новых деформируемых конструкций, построение корректных моделей оболочечных конструкций, оболочек и алгоритмов. Некоторые из разработанных моделей многослойных неспаянных пластин, параметрических колебаний гибких изотропных и орто-тропных пластинок внедрены в расчетную практику НПЦ «Алмаз-Фазатрон», соответствующий документ прилагается к диссертации. Кроме того, результаты работы используются: в учебном процессе при чтении авторского спецкурса (Опубликовано учебное пособие по данной тематике [98]

Апробация работы Основные теоретические и практические результаты диссертационной работы докладывались и обсуждались на Международных, Всесоюзных и Всероссийских научно - технических конференциях: Ш Всесоюзной конференции «Механика неоднородных структур» (Львов, 1991, Украина); Ш Всесоюзном симпозиуме «Устойчивость и пластичность в механике деформируемого твердого тела» (Тверь, 1992, Россия); XVI symposium «Vibrations in physical systems» (Poznan, BlazeJevvko, 1994, Poland); The 3-d conference on dynamical systems «Theory and applications» (Lodz, 1995, Poland); III Международной конференции «Проблемы прочности материалов и сооружений на транспорте» (Санкт-Петербург, 1995, Россия); The 4-th Conference «Dynamical Systems-Theory and Applications» (Lodz, 1997, Poland); Седьмой Международной конференции «Математическое моделирование и краевые задачи» (Самара, 1997, Россия); XVII Международной конференции по теории оболочек и пластин (Саратов, 1997, Россия); VII Четаевской конференции (Казань, 1997, Россия); The 2-nd Belorussian Congress on Theoretical and Applied Mechanics, 'Me-chanics-99' (Minsk, 1999, Belorussia); The 5-th Conference «Dynamical Systems -Theory and Applications» (Lodz, 1999,Poland); XIX Международной конференции по теории оболочек и пластин (Н. Новгород, 1999, Россия); Международной конференции, посвященной 100-летию профессора Х.М. Муштари, 90-летию профессора К.З. Галимова и 80-летию профессора М.С. Корнишина «Актуальные проблемы механики оболочек» (Казань, 2000, Россия); I Международной конференции по нелинейной динамике механических и биологических систем (Саратов, 2000, Россия); Fourth International Colloquium on Computation of Shell & Spatial Structures, IASS-IACM 2000 (Chania-Crete, 2000, Greece); Chaos Control and Times Series (San-Paulo, 2000, Brasil); The Third International Conference on Thin-Walled Structures (Cracow, 2001, Poland); VIII Всероссийском съезде по теоретической и прикладной механике (Пермь, 2001,Россия); European Congress on Computational Methods in Applied Sciences and Engineering, ECCOMAS Computational Fluid Dynamics Conference 2001 (2001, Swansea, Wales, UK); The Seventh PAN American Congress of Applied Mechanics, РАСАМ VII (Temuco, 2002, Chile); Зимней школе по механике сплошных сред (тринадцатой) (Пермь, 2003, Россия); 13-й межвузовской конференции «Математическое моделирование и краевые задачи» (Самара, 2003, Россия). Подтверждением этих выступлений являются тезисы докладов, опубликованные в трудах этих конференций.

В целом работа докладывалась на научном семинаре кафедры «Теоретическая механика» СГТУ под руководством профессора Ю.В. Чеботаревского (Саратов, 2003).

На защиту выносятся:

1) методика построения, на базе множителей Лагранжа, замкнутых математических моделей контактирующих пластин и оболочек, описываемых в рамках одинаковых и различных статических и кинематических гипотез с учетом ширины зоны контакта;

2) постановка и методика численного решения, на базе итерационного метода и метода вариационных итераций контактных задач для пластин переменной толщины, различного рода нелинейностей и произвольного плана в рамках гипотез Винклера;

3) обоснование сходимости итерационных алгоритмов;

4) результаты численных экспериментов;

5) новая классификация эволюционных задач теории многослойных неспаянных пластин; сценарии развития стохастических параметрических колебаний гибких изотропных и ортотропных пластин; новый динамический критерий потери устойчивости пластин при продольных воздействиях;

6) результаты численного эксперимента по исследованию динамических задач теории гибких пластин, полученных методами конечных разностей с ап

О А. проксимацией 0(h ) и О (h ), методом Бубнова в высших приближениях и анализ периодических и хаотических колебаний пластин, полученных с помощью вейвлет - преобразования.

Публикации. По материалам диссертации опубликовано 35 работ [15, 62, 76, 84, 85, 86, 87, 88, 89, 91, 92, 93, 95, 98, 99, 100, 101, 106, 113, 114, 115, 117, 167, 168, 187, 191, 192, 193, 194, 195, 197, 198, 228, 229, 238,], в том числе 2 монографии [86, 106] и разделы в книгах [228, 229], вышедших в издательстве Naukowo - Techniczne Warszawa (2001г.), «Springer» (2003г.).

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

Структура и объем работы Диссертация состоит из введения, 6 глав, заключения и списка используемой литературы, содержащего 280 наименова

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

ЗАКЛЮЧЕНИЕ об использовании материалов докторской диссертации Крысько А В "Математические модели нелинейных распределенных систем в виде пластинчатых и оболочечных конструкций".

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

Разработанные в диссертации численные методы позволили повысить достоверность прогнозирования прочностных характеристик изделий, разрабатываемых в НПЦ

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

Заместитель Генерального директора по науке, -кандидат технических наук, Лауреаздосудар^венной премии СССР ка|^дидат физико-математических наук