Нелинейная динамика трехслойных пластин при периодических и нестационарных воздействиях

Юрченко, Алевтина Анатольевна

Нелинейная динамика трехслойных пластин при периодических и нестационарных воздействиях тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

Юрченко, Алевтина Анатольевна АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Москва МЕСТО ЗАЩИТЫ

2012 ГОД ЗАЩИТЫ

01.02.04 КОД ВАК РФ

Диссертация по механике на тему «Нелинейная динамика трехслойных пластин при периодических и нестационарных воздействиях»

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

На правах рукописи

005009328

Юрченко Алевтина Анатольевна

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

Специальность 01.02.04 - механика деформируемого твердого тела

Автореферат

диссертации на соискание ученой степени кандидата физико-математических наук

2 ОЕЗ

Москва-2012

005009328

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Московский государственный технический университет «МАМИ»

Научный руководитель:

кандидат физико - математических наук, профессор

Коган Ефим Александрович

Официальные оппоненты:

доктор технических наук, профессор

Сухинин Сигизмунд Николаевич

доктор физико - математических наук, доцент

Гришанина Татьяна Витальевна

Ведущая организация

Институт проблем механики им. А.Ю.Ишлинского РАН

Защита состоится 29 февраля 2012 г. в 1600 на заседании диссертационного совета Д212.125.05 в ФГБОУ ВПО Московский авиационный институт (национальный исследовательский университет) по адресу: 125993, г. Москва, А-80, ГСП-3, Волоколамское шоссе, дом 4.

С диссертацией можно ознакомиться в библиотеке Московского авиационного института (национального исследовательского университета)

Автореферат разослан...-.........2012 г.

Ученый секретарь

диссертационного совета Д212.125.05 —' Федотенков Г.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

В ней отражены исследования, проводившиеся в 2004-2011 годах на-кафедре прикладной и вычислительной математики МГТУ «МАМИ» и в НИИ механики МГУ им. М.В.Ломоносова.

Целью диссертационной работы является:

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

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

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

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

Основные положения, выносимые на защиту:

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

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

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

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

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

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

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

- на семинаре по механике деформируемого твердого тела под руководством чл.-корр. РАН Э.И.Григолюка (2004, 2005 гг. и позднее в 2006-2010гг. в НИИ Механики МГУ им. М.В.Ломоносова); на «Девятой международной научно - технической конференции по динамике и прочности автомобиля», 15-17 марта 2005 года (М.: МГТУ «МАМИ»); на «Международном научном симпозиуме, посвященном 140-летию МГТУ «МАМИ», 23-24 марта 2005 года (Россия, М.: МГТУ «МАМИ»); на конференции - конкурсе молодых ученых НИИ Механики МГУ им. М.В.Ломоносова в октябре 2006 года; на научных конференциях «Ломоносовские чтения», секция механики. Апрель 2004

- 2008 годов. М.: МГУ им. М.В .Ломоносова; на Международном симпозиуме им. А.Г.Горшкова «Динамические и технологические проблемы механики конструкций и сплошных сред», Ярополец, Москва, 16-20 февраля 2009 года. Часть результатов работы получена и использовалась в рамках исследований по гранту РФФИ № 10-08-00258а.

Публикации. Список научных трудов по диссертационной работе составляет 8 публикаций.

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

Струюура и объем диссертации. Результаты изложены на 150 страницах машинописного текста, иллюстрированного 41 рисунком.

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

КРАТКАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

В.И.Мяченкова, Ю.В.Немировского, В.Н.Паймушина, В.И.Шалашилина, Л.А.Шаповалова, А.Р. ВогеБ1, С. СЫа, W.T. Койег'а, ЯЛ.МкоЫ, Н.Т. \Voer-п1е и др.

Важный резерв прочности и оптимизации конструкций был реализован, начиная с 40-х годов прошлого века, в связи с применением трехслойных и многослойных пластин и оболочек. Основной отличительной особенностью расчета трехслойных конструкций с маложестким промежуточным средним слоем является необходимость учета поперечного сдвига и поперечных нормальных напряжений и деформаций в заполнителе. Значительный вклад в развитие теории трехслойных и многослойных пластин и оболочек внесен А.Я.Александровым, С.А.Амбарцумяном, В.Н.Бакулиным, В.В.Болотиным, В.В.Васильевым, Н.К.Галимовым, Э.И.Григолюком,

B.Н.Кобелевым, Г.М.Куликовым, Ю.В.Немировским, В.Н.Паймушиным,

C.Н.Сухининым, ПЛ.Чулковым, Е.Н^зпег'ом и др. Различные аспекты исследований в этой области отражены в обстоятельных обзорах С.А.Амбарцумяна, В.В.Болотина с Э.И.Григолюком и Ю.Н.Новичковым, К.З.Галимова, Э.И.Григолюка с Е.А.Коганом, Г.М.Куликовым, И.Т.Селезовым и В.И.Мамаем, Я.М.Григоренко и А.Т.Василенко,

A.А.Дудченко, И.Ф.Образцова и А.СЛурье, Ю.В. Немировского,

B.В.Чепиги, А.Дуды, Л.М.Хэбипа, Ф. Плантема и др.

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

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

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

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

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

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

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

Поэтому на основании анализа состояния вопроса были сформулированы цели работы, указанные выше.

Во второй главе работы предложен единый подход к расчету нелинейных свободных колебаний и вынужденных колебаний трехслойных пластин конечного прогиба под действием периодически изменяющейся по времени поперечной нагрузки. Рассмотрены нелинейные колебания трехслойных пластин несимметричной структуры по толщине с жестким трансвер-сально изотропным заполнителем, податливым на поперечный сдвиг, и изотропными несущими слоями. Вынужденные колебания таких пластин описываются уравнениями Григолюка-Чулкова, в которых учтены также начальные неправильности формы координатной поверхности, поперечные инерционные силы и внешнее демпфирование. Эта система уравнений 10-го порядка в смешанной форме относительно разрешающих функций перемещений х и усилий /г, обобщающая уравнения однородных пластин Феппля-Кармана, имеет вид

У2У 2Р = ЕН

5Л

( Л2 О IV

2 Л

агуу + д\ д2\у д\ д2\у дг\ув

' дх2 ' ду2 + дхду' дхду дх2 ' ду2 8уг дх2

(1)

(д2Ъ 344

б*2 +'дх2

д2Г ду2 '

/-2— э2 Л О М> О IV.

ду2 ду2

а2^ а*2

+ 2

^ д2ц> д2™.

дхду дхду ,

= д-рИ

л52 5 —• + £■— 5/2 Э/

1-^У2

(2)

В уравнениях \¥(х,,у,0 = й'(*,.у,0 + (*,.)') - полный прогиб, и'е(х,у)- начальный прогиб, характеризующий отклонение пластины от идеальной формы, ъ/{х,у,Г) - дополнительный прогиб, выражающийся через разрешающую функцию перемещений % известным в теории трехслойных пластин соотношением:

м> =

1-^У2

. Р

1-у2 3 ЕЙ

Е =-У * 1 - осреднённый модуль упругости трехслойного пакета,

/г = - полная толщина пластины толщина к-го слоя, (к = 1,2,3),

*=1

к= 2с -толщина заполнителя), £> =

12(1-V2)

изгибная жесткость трех-

слойного пакета, р ■■

\201ф-уг) Ецх

параметр, характеризующий жесткость

заполнителя на поперечный сдвиг (<3модуль поперечного сдвига заполнителя (/ = 1,2)), 9 = ^^—— - параметр, характеризующий изгибную жест-

v ЕЛУк кость несущих слоев, V = 2^

ы 1-У

3 Г А

П-к2

- приведённый коэффициент

"к Ч"1 * 'I У

Пуассона, (Ек,ук - модуль упругости и коэффициент Пуассона к-то слоя ),

рИ = ^РкК (рк - удельная плотность материала к-то слоя), коэффициенты

»=1

т]х,т]2,1]} зависят от безразмерных жесткостных характеристики и без-

ЕА

1 -V

^ ЕЛ

-V2

размерных толщин слоев (А: = 1,2,3): х* :

Граничные условия для свободно опёртой по контуру пластины = У2/Г = 2' = У22' = У2У2^ = 0 удовлетворяются представлением разрешающей функции перемещений % и начального прогиба и виде

. . тях . пяу

X = f(t) sin-sin—S

a b

. . тях . пяу

w0=/„sin-sin—-

a b

В случае жесткого защемления по контуру трехслойной пластины граничные условия относительно функций F и х записываются в виде:

дх 8х}

п dF d'F При х = 0, а — = —- =

дх дх

. Р

При y = 0,b ^ - з

dF = d'F ду~дуг

д%_=д\х_ = 0

я ду ду1

(3)

Для построения координатных функций, удовлетворяющих условиям (3), задавался прогиб пластины в виде, удовлетворяющем граничным условиям относительно функции к и ее нормальной производной и1'

(4)

где

„ . . тях

Zm{ д;) = cos--cos

(

и соотношение vv =

w(x,y,t) = W(t)Za(x)Zn(y), (m + 2)ях

- плу (п + 2)яу

z„(y) = C0S^ ~C0S b '

(5)

1-*V

р .

X разрешалось относительно функции х- В

результате для этой функции получено выражение вида (6), содержащее амплитудное значение IV(/) и удовлетворяющее граничным условиям (4):

Xix,y,t) = W(t}

тях плу

cos-cos——

а_b

тлх (л + 2)лу

cos-cos --—

а Ъ_

1 + /С

+ п1

г V

1+*. [Ш +(п+2)1

(т + 2 )лх ппу

cos ---— cos-

а_Ь_

(т + 2)лх (п + 2 )яу

cos ---— cos --

а Ъ

1 + К„

т + 2

в котором

Кп =

я'И1

рь1'

+ п

Л = -

1+К,

т + 2

+ (и + 2)2

(6)

Для функции усилий F граничные условия жесткого защемления удовлетворяются, если принять F(x,y,t) = f(t)Zin{x)Z„(y),

Уравнения (1), (2) для шарнирно опёртых пластин интегрировались методом Бубнова - Папковича по пространственным координатам.

Для защемленных пластин решение краевой задачи (1-3) методом Бубнова - Папковича при выбранных координатных функциях (6) приводит к весьма трудоёмкой процедуре построения частного решения неоднородного

уравнения совместности деформаций (1). Поэтому метод ортогонализации Бубнова применялся непосредственно к системе уравнений (1), (2).

В результате для обоих вариантов граничных условий получено одно и то же по структуре нелинейное обыкновенное дифференциальное уравнение, описывающее в первом приближении вынужденные колебания трехслойной пластины под действием внешней поперечной нагрузки, изменяющейся по гармоническому закону q(x,y,t) = Q(x,y)cosQt и равномерно распределенной по поверхности пластины (Q{x,y) = Qa= const):

^f + +a2? +ал)=до(0, (?)

at at

В (7) <а02тп - квадрат частоты собственных малых колебаний трехслойной пластины. Для шарнирно опертых пластин параметры, фигурирующие в уравнении, отражающие геометрическую и механическую неоднородность трехслойного пакета по толщине, в предельном случае однослойной пластины при неучёте поперечного сдвига полностью совпадают с известными в литературе.

Для случая жесткого защемления при неучёте начальных прогибов аг= 0, а, =1,

4 f*? U,.

=«1* = -24\~у2^■ 1тп ■ [м2 JimJ!„ -m2(n + 2f J,mJAn -(rn + 2)V X Хщ

xJ*mJ3n +(f" + 2)2(n + 2)2JAmJ4n -J5mJ5n]-(NuJlmJu-Nl2J^J2„ -

- N2lJ2mJ]„ + N22J2mJ2„) ■

а интегралы метода Бубнова равны

JXm = j cos-Zk (x)dx, Jln = f cos -AZ, (y)dy, J2m = cos ^-'—Zk (x)dx,

о а о b о a

J2„=\cos^" + 'f>JI}-Zl(y)dy, J3m = Jcos —Zm{x)Zk(x)dx,

о ° о a

J3n = Jcos ^Z„(y)Z,(y)dy, JAm = jcos im±^Zm(x)Zk(x)dx, 0 " о a

J A„ = j cos -Zk(y)Z, {y)dy, Iqmn = }\q(x,y,t)Zk{x)Z,(y)dxdy,

a b oo

Jsm=\

2 . гтт „ . ... mm . (m + 2W 2

m sin--2m(m + 2) sin--sin^->— + (m + 2) x

a a a

. 2 (m + 2W

x sin ---—

Zk(x)dx.

^„=1 и2 ет2 -2п(п + 2)БШ•

Л Ь

,2 -2 ПТГУ

+(и+2)г8т2^±Щг,(у)4у. ь

Коэффициенты (/,7=1,2) зависят от относительных геометрических и

жесткостных характеристик слоев и параметров волнообразования.

Принимая далее прогиб пластины в виде х = ЛсобП/ и интегрируя уравнение (7) по полному периоду колебаний Г = 2я7П, получим уравнение, связывающее амплитуду колебаний А с относительной частотой нелинейных колебаний V = &/со0

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

Резонансные амплитуды А и частоты колебаний со^ трехслойной

пластины определяются как координаты точки пересечения амплитудно -частотной и скелетной кривых.

Расчеты выполнены для трехслойных пластин, геометрические и жест-костные параметры которых для основного варианта (несущие слои толщиной А, = И2 =0,15 см, заполнитель толщиной =5,7 см; модули упругости несущих слоев -Ет =7-104МПа, в заполнителе модуль упругости Ет = 40 МПа) соответствуют секции реального объекта. Рассматривались пластины с относительной толщиной И/Ь = 0,1 и пластины, имеющие достаточно большие размеры в плане {И/Ь и 0,033). Варьировались также относительное удлинение пластин Я = а/Ь, модуль упругости заполнителя от Е(1) = 40 МПа до Е01 =1500 МПа, собственная изгибная жесткость несущих слоев (изменением их толщины), параметр г, характеризующий демпфирование при колебаниях, начальные неправильности формы координатной поверхности и числа полуволн т и п.

Ниже приведены некоторые характерные графики, отражающие зависимость амплитудно - частотных характеристик от изгибной жесткости несущих слоев и от относительной толщины Л / Ь для шарнирно опертых и для защемленных по контуру трехслойных пластин. Как видно, изменение изгибной жесткости несущих слоев существенно сказывается на кривых А(у), а именно увеличение параметра 3 (от значения 3,28-10~\ что соответствует толщинам несущих слоев Л, = 0,2 см и Л2 = 0,15 см при общей толщине пакета слоев И = 6 см (кривые ¿>1 на рисунках 1,2) до значения 9 »1,78-10"3 (кривые >93) приводит к заметному снижению максимального значения амплитуды колебаний, но резонансная частота при этом практически не меняется. При этом графики 1а, 2а соответствуют значению И/Ь = 0,1, а графики 16,26 - значению И/Ь = 0,033.

06 А

0.! 04

0 3 02

- 1)1 ----1)2 ........ 1)3

Рис.2

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

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

В третьей главе работы приведена методика расчета и выполнено численное исследование напряженно - деформированного состояния трехслойных свободно опёртых прямоугольных пластин при кратковременном нестационарном нагружении. Вынужденные колебания трехслойных пластин конечного прогиба несимметричной структуры по толщине с жестким транс-версально изотропным заполнителем, податливым на поперечный сдвиг, и изотропными несущими слоями, описываются, как и в главе 1, уравнениями Григолюка-Чулкова (1), (2).

Рассматривались два случая нестационарного воздействия на пластину: кратковременная нагрузка, изменяющаяся по времени по кусочно-линейному закону (таким законом изменения нагрузки может моделироваться, например, воздействие наземного взрыва заряда взрывчатого вещества под днищем автомобиля) и ступенчатый (прямоугольный) импульс малой интенсивности, но относительно большой длительности по времени (таким нагружением .может моделироваться воздействие сильного взрыва на большом.удалении от объекта, в частности, на панели крыши кузовов - контейнеров, которые имеют достаточно большие размеры в плане).

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

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

„ п„ . тх . ¡лу . 1жх . )лу

р = Е Ъ(!)■■81П-81П , х = £ £ X,, С)!МП-МП V1.

/у а О ¡1 а Ь

^ X— . 1ЯХ . ]Яу = 2.1^00 МП—Iш —,

где /,_/- числа полуволн соответственно вдоль сторон аиЬ пластины.

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

и% ¿Хо

ат

т^1

1 I j к I

X - 2ук1 • /5,ий • 16,П]1 ]х (Д7 Х,/к> + ™оц ТХ-^-^ЪтпЧ = 0,

I } Ру

-,2

Я"

кт1г,щрч + Т^ШЕк/^'м;! х \Г11РиХчХи +

ЛЯ Л / ¡к I

' 1

+ + У2*2]- + Зй'о.и^)} = 0. (от,я = 1,2,...)

где /1т,.,...,/6п;/- интегралы метода Бубнова.

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

мировать правые части уравнений нормальной системы и находить значения IV в следующие моменты времени.

Выполнена оценка сходимости решения по сходимости частичных сумм рядов. При воздействии кусочно - линейного импульса на пластины с И/Ь = 0,1 решение получено до значений =16. Достоверность рас-

четов проверялась варьированием шага интегрирования дифференциальных уравнений. Некоторые характерные кривые изменения максимального прогиба в центре пластины в зависимости от безразмерного времени т = Сл/а (С. - скорость звука в пластине) при действии кусочно-линейного импульса для разных значений изгибной жесткости несущих слоев и жесткости заполнителя на сдвиг показаны на графиках 3, 4, для пластин с разными размерами в плане на рис. 5,6.

— р'ЧОШа — р'«150МПа......•1500МПа

Рис.4

0.1

0.0

-0,1

Л Л

Лл/ \л„ |Л„

■ма! ' \р

а=Ь=60 Рис.5

а=Ь=200 Рис.6

Максимальное значение г, которым ограничивались в расчетах гш =100, соответствует /»0,04 сек. Поэтому, как видно из графиков, максимальное значение прогиба достигается еще в период действия кусочно -линейного импульса, а затем пластина совершает свободные затухающие колебания.

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

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

При действии ступенчатого импульса малой интенсивности, но большой длительности по времени максимальные прогибы определены при значениях /ишх = ишх = 10. При этом относительная погрешность в определении Wm% составила около 4%. Максимальный прогиб при этом достигается еще до исчезновения силы, а затем пластина совершает колебания относительно деформированного состояния.

Сравнение с решением линейной задачи показывает, что при нагрузке, изменяющейся по времени по кусочно-линейному закону, влияние нелинейности на деформированное состояние пластин с И/а = 0,1 незначительно, но для относительно тонких пластин больших размеров в плане учёт нелинейности приводит к заметному снижению максимальных прогибов.

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

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

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

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

, w' + w2 z vv1 - W1 Z , ^ ^ Л

координаты w =-+--= w+—v, (-cszsc)

2 с 2 с

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

Применение принципа Гамильтона Остроградского приводит к вариационному уравнению вида

(, м-'=1 ,'ж1-'=1

2 ¡=1

/=17—1 ?

<2 Г"

+ д, - д2))&'}с1х1сЬс2)с11 + Щ[ТУ22Л2 + Л^й/, + Я22&2 + #12&, + (М22 2 + + 2 М12_, + N+^22у,2 +Л^12У„ +223 - М22йи',2+(/,22)2 + 2 ¿121 +

+ ^>,2+^,>.,+^22^.2+^12^1 ¿22^2 + Щ^П^! + ^12^2 +

I, и

+ Ни5ах + Я12&2 + (МП1 + 2Л/12д + Лгпм^+Л'12™,2+Л^.1+А^>2+бп)<5'1'--Л/и^„+(11и +21122 )сЛ —

В (8) Ь - действие по Гамильтону, Я - потенциальная энергия деформации пластины, К -кинетическая энергия системы, А - работа внешней поверхностной нагрузки,

Л/ , Нн, Цр 0„ - полные и обобщенные удельные усилия и моменты, выражающиеся через напряжения в слоях.

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

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

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

д2Щ | д21 - у д2Щ | я1 + у д2й2 | И й*о 2 ду1 2 8х0ду0 а,

- + Л

дхо дх2й 2 дх0 ду20

+ У дй> д2м

+ Я---

2 ду0 8х08у0

= 0,

¥о

,дги2 1-v д2и2 1 +V д2щ

'--2~+---

■ дх0

^ Л

-+я-

2 дх0ду0 а

, 8w 82w 1-v 8w d2\v

■---1----

SVo 9уо 2 ^o Sxl

1 + v dw 8 w 2 8x0 дх0ду0

- + Л2

= 0,

d2a¡. .21-v d2a¡ ,1 + v 82a¡ , д ÍL á( - . hdw^

— i J----^ д.--

8xq 2 dyl 2 дх0ду0 8x0

адх 0

= 0,

8% ■l-vd2á, ,1 + v д2а.

ду20 2 8x1 d2w , dw

2 8х0ду0 8y0 У, аду0

hdw^

= 0,

дт2

- + dn — + d2A0á0w-A: дт

h t _ (да, . да2 а V&0 ду0

д&о

\-vd%

2 дх.

2 дх08у0

2 ду02)дх,

, . , да, , да, дх0 дУ0

2д2й2 Л + v 52й, - + Я

Фо

2 &05v0

3w [ йн, . ди2

о У

—+

ЗУо V5xo дУо]ох0 *

ди, За, 13 w

+ Л Я^+v

,,, , Su, ди2) d2w

+ A(l-v) Я—L + v—- -

дуй дх0)дх0ду0

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

Граничные условия жесткого защемления удовлетворялись представлением компонент вектора перемещений в виде разложения в ряды по фундаментальным балочным функциям:

" a, =Zl4(r)Z:(x0)Zj(>;J)

> i ' у

¡ У ' /

' у

Здесь Z,(jt0) = sin Я,дг0 - íAA*0 - а, (cos Я д:0 - с/гЯх0)

2,00 = sin Я,}>, - лЛЯд, - а,(cos Я,у0 - сЛЯ^0).

Как и в главе 3, после интегрирования методом Бубнова полученная система обыкновенных нелинейных дифференциальных уравнений редуцировалась к нормальной системе и интегрировалась с учетом нулевых начальных условий методом Кутта-Мерсона. Интегралы процедуры Бубнова вычислялись численно методом Симпсона.

Проведено численное исследование полученного решения в зависимости от числа удерживаемых членов рядов. При воздействии кусочно - линейного импульса на пластины с И/Ъ = 0,1 решение получено до значений ттх=пш^= 10. Достоверность расчетов проверялась варьированием шага интегрирования дифференциальных уравнений. Расчеты выполнялись для трехслойных пластин с теми же геометрическими и жесткостными параметрами, что и в главе 3.

Некоторые характерные результаты численных расчетов приведены' ниже на графиках 7, 8, на которых показана зависимость й^ от изгибной жесткости несущих слоев (рис. 7) и жесткости заполнителя на сдвиг (рис. 8) при действии прямоугольного импульса.

0.1 V 0.0 '

(^ йу

ш 0

в/

¡Г

-¿•••••тейп»—1?"««смп».....£« «шип»

Рис.7

Рис. 8

Влияние геометрических параметров - относительного удлинения пластин Х = а!Ъ и тонкостенности пластин И1Ь качественно такое же, как и для шарнирно опёртых пластин: с ростом Я, как и для однослойных пластин при статическом изгибе, прогибы существенно возрастают. Отметим, что с увеличением линейных размеров при кратковременных нестационарных воздействиях как для шарнирно опертых, так и для жестко защемленных пластин, прогибы увеличиваются, но медленнее, чем при статическом нагру-жении, а именно примерно пропорционально квадрату линейного размера, что согласуется с известными данными для однослойных пластин.

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

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

В заключении работы даны общие выводы по полученным результатам.

Основные выводы и результаты работы

1. Предложен единый подход к решению задач о нелинейных колебаниях под действием периодически изменяющейся по времени поперечной на-

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

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

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

Основное содержание диссертации опубликовано в следующих работах:

Статьи в журналах из перечня ВАК:

1. Коган Е.А., Юрченко A.A. Нелинейные колебания защемленных по контуру трехслойных пластин // Проблемы машиностроения и надежности машин. М.: 2010, № 5. С. 25-34.

2. Юрченко A.A. Численное моделирование деформирования трехслойных панелей при кратковременных нестационарных воздействиях // Проблемы машиностроения и надежности машин. М.: 2011, № 4. С. 85-92.

Прочие публикации:

1. Григолюк Э.И., Коган Е.А., Юрченко A.A. О численном решении задач нелинейной динамики неоднородных пластин при нестационарном нагружении // Ломоносовские чтения. Секция механики. Научная конференция. Апрель 2004 года. Тезисы докладов. М.: Изд-во Моск. ун-та. С. 64- 65.

2. Коган Е.А., Юрченко A.A. Деформирование трехслойных пластин при импульсном нагружении // Материалы IX Международной научно - технической конференции по динамике и прочности автомобиля. 15-17 марта 2005 года. Под редакцией чл. - корр. РАН Э.И.Григолюка. М.: МГТУ «МАМИ», 2005. С. 172 -176.

3. Коган Е.А., Юрченко A.A. О колебаниях трехслойных пластин конечного прогиба при кратковременном импульсном нагружении // Международный научный симпозиум, посвященный 140 - летаю МГТУ «МАМИ» 2324 марта 2005 года. Материалы 49-ой Международной научно - технической конференции ААИ «Приоритеты развития отечественного автотракторостроения и подготовки инженерных и научных кадров». Секция 4. «Математические методы моделирования и оптимизации автотранспортных средств». Часть 2. М.: МАМИ, 2005. С. 23-26.

4. Юрченко A.A. Нелинейное деформирование трехслойных пластин при кратковременном динамическом нагружении // Труды конференции -конкурса молодых ученых 11-16 октября 2006 г. / Под ред. акад. РАН Г.Г.Черного, проф. В.А.Самсонова. М.: изд-во Моек ун-та, 2007, С. 304310.

5. Коган Е.А., Юрченко A.A. Нелинейное деформирование защемленных по контуру трехслойных пластин при импульсном нагружении // «Избранные проблемы прочности современного машиностроения». Сборник научных трудов, посвященный 85-летию чл.-корр. РАН Э.И.Григолюка. М.: Физматлит. 2008. С. 110-123.

6. Коган Е.А., Юрченко A.A. Построение амплитудно - частотных характеристик трехслойных пластин конечного прогиба // Динамические и технологические проблемы механики конструкций и сплошных сред. Международный симпозиум им. А.Г.Горшкова. Ярополец, Москва, 16-20 февраля 2009 года. Тезисы докладов М.: МАИ. 2009. С. 89-90.

Подписано в печать: 18.01.2012 Объем: 1,5 усл.п.л. Тираж: 100 экз. Заказ № 145 Отпечатано в типографии «Реглет» 119526, г. Москва, Страстной бульвар, д. 6,стр. 1 (495) 978-43-34; www.reglet.ru

Текст научной работы диссертации и автореферата по механике, кандидата физико-математических наук, Юрченко, Алевтина Анатольевна, Москва

61 12-1/465

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

ФГБОУ ВПО «МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ «МАМИ»

На правах рукописи

Юрченко Алевтина Анатольевна

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

Специальность 01.02.04 - механика деформируемого твердого тела

ДИССЕРТАЦИЯ

на соискание ученой степени кандидата физико - математических наук

Научный руководитель кандидат физико-математических наук, проф. Коган Е.А.

МОСКВА 2012

ОГЛАВЛЕНИЕ

Стр.

ВВЕДЕНИЕ..................................................................... 4

1. АНАЛИЗ СОВРЕМЕННОГО СОСТОЯНИЯ ИССЛЕДОВАНИЙ В ОБЛАСТИ НЕЛИНЕЙНОЙ

ДИНАМИКИ СЛОИСТЫХ ПЛАСТИН.............................. 11

Введение..................................................................... 11

1.1. Нелинейные колебания неоднородных пластин при периодических воздействиях...................................... 16

1.2. Динамическая реакция слоистых пластин при импульсном нагружении..................................................... 27

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

2.1. Построение амплитудно - частотных характеристик шарнирно опертых трехслойных пластин..................... 37

2.1.1. Решение методом гармонического баланса......... 46

2.1.2. Результаты расчетов..................................... 48

2.2. Построение амплитудно - частотных характеристик защемленных по контуру трехслойных пластин............ 58

2.2.1. Решение методом гармонического баланса......... 66

2.2.2. Результаты расчетов...................................... 67

2.3. Краткие выводы по главе......................................... 72

3. ИССЛЕДОВАНИЕ НАПРЯЖЕННО - ДЕФОРМИРОВАННОГО СОСТОЯНИЯ СВОБОДНО ОПЕРТЫХ ТРЕХСЛОЙНЫХ ПЛАСТИН КОНЕЧНОГО ПРОГИБА ПРИ КРАТКОВРЕМЕННОМ ДИНАМИЧЕСКОМ НАГРУЖЕНИИ................... 73

3.1. Динамическая реакция свободно опертых прямоугольных трехслойных пластин при воздействии взрыва 73 3.1.1. Постановка задачи........................................ 73

3.1.2. Решение начально - краевой задачи....................................77

3.1.3. Определение силовых факторов..............................................81

3.2 Воздействие кусочно - линейного импульса..................................87

3.3 Воздействие прямоугольного импульса..............................................98

3.4 Краткие выводы по главе..................................................................................100

4 НЕЛИНЕЙНОЕ ДЕФОРМИРОВАНИЕ ЗАЩЕМЛЕННЫХ

ПО КОНТУРУ ТРЕХСЛОЙНЫХ ПЛАСТИН ПРИ

НЕСТАЦИОНАРНОМ НАГРУЖЕНИИ............................................................101

4.1. О выборе аппроксимирующих функций..............................................101

4.2. Построение определяющих соотношений нелинейной динамики трехслойных ортотропных пластин............................104

4.3. Решение начально - краевой задачи........................................................119

4.4. Анализ результатов расчетов..........................................................................123

4.5. Краткие выводы по главе..................................................................................131

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ................................................................132

ЛИТЕРАТУРА......................................................................................................................................134

ВВЕДЕНИЕ

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

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

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

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

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

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

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

Настоящая работа и посвящена исследованию нелинейного поведения трехслойных пластин конечного прогиба при периодических и нестационарных воздействиях.

В ней отражены исследования, проводившиеся в 2004-2011 годах на кафедре прикладной и вычислительной математики ФГБОУ ВПО МГТУ "МАМИ" и в НИИ механики МГУ им. М.В.Ломоносова.

Целью диссертационной работы является:

Структура и объем диссертации. Результаты изложены на 150 страницах машинописного текста, иллюстрированного 41 рисунком.

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

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

На основании анализа состояния вопроса были сформулированы цели работы, указанные выше.

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

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

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

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

1). Кратковременная динамическая нагрузка, меняющаяся по времени по кусочно-линейному закону,

2) ступенчатый (прямоугольный) импульс малой интенсивности, но относительно большой длительности по времени.

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

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

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

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

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

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

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

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

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

щему колебания систем с вязким трением и с кубической упругой характеристикой.

Основные положения, выносимые на защиту:

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

Апробация работы. Основные результаты диссертационной работы доложены на научных конференциях и семинарах, в том числе: - на семинаре по механике деформируемого твердого тела под руководством чл.-корр. РАН Э.И.Григолюка (2004, 2005 гг., НИИ механики МГУ им. М.В.Ломоносова), на «Девятой международной научно - технической конференции по динамике и прочности автомобиля», 15-17 марта 2005 года (М.: МГТУ «МАМИ»), на «Международном научном симпозиуме, посвященном 140-летию МГТУ «МАМИ», 23-24 марта 2005 года (Россия, М.: МГТУ «МАМИ»), в трудах конференции - конкурсе молодых ученых. Октябрь 2006 года / Под ред. акад. РАН Г.Г.Черного, проф. В.А.Самсонова. М.: изд-во Моск. ун-та, 2007, на научных конференциях «Ломоносовские чтения», секция механики. 2004 - 2008 гг. М.: МГУ им. М.В.Ломоносова), на международном симпозиуме им. А.Г.Горшкова «Динамические и технологические проблемы механики конструкций и сплошных сред». Ярополец, Москва, 1620 февраля 2009 года.

Публикации. Список научных трудов по диссертационной работе составляет 8 публикаций.

1. АНАЛИЗ СОВРЕМЕННОГО СОСТОЯНИЯ ИССЛЕДОВАНИЙ В ОБЛАСТИ НЕЛИНЕЙНОЙ ДИНАМИКИ

СЛОИСТЫХ ПЛАСТИН

Введение

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

Геометрически нелинейная теория пластин стала интенсивно развиваться с начала XX столетия в связи с потребностями кораблестроения. Основополагающий вклад в эту теорию был внесен отечественным ученым - кораблестроителем И.Г.Бубновым [1,2]. Дальнейшее развитие теории было дано в трудах А.Фёппля [3], Т. Кармана [4] и К.Маргерра [5]. Наиболее существенные результаты, п