Напряженно-деформированное состояние гибких слоистых пластин из анизотропных разносопротивляющихся материалов

Васильев, Николай Владимирович

Напряженно-деформированное состояние гибких слоистых пластин из анизотропных разносопротивляющихся материалов тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

Васильев, Николай Владимирович АВТОР

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

Тула МЕСТО ЗАЩИТЫ

2012 ГОД ЗАЩИТЫ

01.02.04 КОД ВАК РФ

Диссертация по механике на тему «Напряженно-деформированное состояние гибких слоистых пластин из анизотропных разносопротивляющихся материалов»

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

005015917

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

Васильев Николай Владимирович

НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ ГИБКИХ СЛОИСТЫХ ПЛАСТИН ИЗ АНИЗОТРОПНЫХ РАЗНОСОПРОТИВЛЯЮЩИХСЯ МАТЕРИАЛОВ

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

АВТОРЕФЕРАТ

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

3 [/(А.7 ¿0!Л

Тверь 2012

005015917

Работа выполнена в ФГБОУ ВПО «Тульский государственный университет».

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

доктор технических наук, профессор Трещев Александр Анатольевич

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

Охлопков Николай Леонидович доктор технических наук, профессор ФГБОУ ТГТУ, профессор каф. «Сопротивления материалов, теории упругости и пластичности»

Божанов Павел Валерьевич кандидат технических наук, доцент ООО «Инженерный Центр Промышленного Проектирования», генеральный директор

Ведущая организация ФГУП ГНПП «Сплав», г. Тула

Защита состоится «31» мая 2012 г. в_на заседании диссертационного совета Д 212.262.02 при ФГБОУ ВПО «Тверской государственный технический университет» по адресу: 170026, г. Тверь, наб. Афанасия Никитина 22, ауд. Ц-312.

С диссертацией можно ознакомиться в библиотеке Тверского государственного технического университета.

Автореферат разослан «_» _2012 года.

Ученый секретарь диссертационного совета

Гультяев Вадим Иванович

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

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

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

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

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

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

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

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

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

Задачи исследования. Для достижения поставленной цели необходимо решить следующие задачи:

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

2. Выполнить сравнение принятой модели деформирования анизотропных разносопро-тивляющихся материалов и наиболее известных моделей с экспериментальными диаграммами по деформированию этих материалов.

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

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

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

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

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

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

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

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

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

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

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

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

получением теоретических результатов строгими математическими методами, основанными на

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

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

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

Практическая значимость работы, выполненной в рамках госбюджетной НИР ТулГУ № 27.06 "Актуальные проблемы технологии строительных материалов и проектирования конструкций", заключается в следующих результатах:

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

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

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

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

Внедрение результатов работы осуществлено в организациях: ООО «Строительное проектирование» (г. Тула), ОАО «7УЛАОБЛГАЗ» {г. Тула). Программный продукт используется указанными предприятиями для экспертизы ресурса прочности конструкций при проведении проектных работ, НИР и ОКР.

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

• IV молодежная научно-техническая конференция студентов Тульского государственного университета «Молодежные инновации», Тула, ТулГУ, 2010 г.

• XIV Международная научно-практическая конференция "Современные технологии в машиностроении", Пенза, Приволжский Дом знаний, 2010 г.

• Международная научно-практическая конференция "Проблемы современного строительства", Пенза, Приволжский Дом знаний, 2010 г.

• XII Международная научно-техническая конференция «Актуальные проблемы строительства и строительной индустрии», Тула, ТулГУ, 2011 г.

• XII Международная научная конференция «Новые идеи нового века» ФАД ТОГУ, Хабаровск, Тихоокеанский государственный университет, 2012 г.

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

В полном объеме диссертация докладывалась на расширенном заседании кафедры «ССМиК» Тульского государственного университета 2 апреля 2012 года, 31 января 2012 года на научном семинаре по МДТТ при ГОУ ВПО «Тверской государственный технический университет», под руководством председателя совета доктора техн. наук, профессора В.Г. Зубчанинова, а также 14 марта 2012 года на научном семинаре по МДТТ им. Л.А. Толокон-никова при ГОУ ВПО «Тульский государственный университет», под руководством председателя совета доктора физ.-мат. наук, профессора A.A. Маркина.

Публикации. По теме диссертационной работы опубликовано 8 печатных работы. Основное содержание диссертации отражено в 8 публикациях, в том числе 3 работы в изданиях рекомендуемых ВАК РФ.

Структура и объем диссертации. Диссертация состоит из введения, трех разделов, заключения, списка литературы из 121 наименований и приложения. Диссертация содержит 117 страниц основного текста, в том числе 22 рисунка и приложение на 49 страницах, содержащего результаты и текст программы расчета круглой многослойной пластины, документы о внедрении. Общий объем работы -166 страниц.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

Первую группу составляют модели, в основу которых положена зависимость механических характеристик материала от знаков возникающих напряжений или развивающихся деформаций. В рамках данной модели выделены определяющие соотношения, предложенные в работах С.А. Амбарцумяна, A.A. Хачатряна, М.С. Саркисян, P.M. Джонса, К.В. Берта, Ф. Табаддо-ра, В.В. Петрова, И.Г. Овчинникова и других.

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

Экспериментальные исследования, связанные с деформированием анизотропных материалов проводились В.Г. Зубчаниновым, Н.Л. Охлопковым, A.B. Березиным, В.И. Строковым, Г.С. Писаренко, Е.В. Ломакиным, К.Н. Русинко, В.А. Паняевым, М.Я. Леоновым, P.M. Джонсом, Д.А.Р. Нельсоном.

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

• не учет влияния сложных видов напряженного состояния при определении жесткости материала;

• привлечение к расчету кусочных и непотенциальных зависимостей;

• большое количество констант, входящих в определяющие соотношения;

• наличие ограничений, накладываемых на некоррелируемые константы материалов;

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

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

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

Существенный вклад в развитие теории слоистых конструкций внесли А.Я. Александров, С.А Амбарцумян, В.В. Болотин, К.З. Галимов, Э.И. Григопюк, Я.М. Григоренко, Х.М. Муштари, В.В. Пикуль, A.A. Маркин и др.

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

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

теории Кирхгофа-Лява.

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

Напряженное состояние в точке анизотропного деформируемого тела в нормированном пространстве главных напряжений количественно предложено определить модулем вектора полного напряжения S = ^ak-ak, а качественно - косинусами направляющих углов ак =cosßt, которые связаны условием нормировки ак-ак=\,гр,ек='\ ,2,3.

Модуль вектора S можно рассматривать как норму векторного пространства главных напряжений, а направляющие косинусы ак - как главные нормированные напряжения.

При переходе от главных осей к произвольной ортогональной системе координат имеем следующие зависимости

5 = 'crv ■,ocIJ=(T,j/S\aIJ-aIJ=l. (1)

На базе введенных количественной и качественных характеристик напряженного состояния в работах Н.М. Матченко, З.В. Аркания и A.A. Трещева определяющие соотношения для анизотропных разносопротивляющихся материалов было предложено формулировать в виде потенциала деформаций W = A(av)S2, где А(ац) - стеленной полином от нормирован-

ных напряжений ai}, конкретизация которого производится в зависимости от вида структурной анизотропии.

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

Si = °Э1 = °12 S7 =СГ32 = = °21 • (2)

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

A(ak) = a0 + alai + akmatam/2-l W = BkSkS+BkmSkSm /2; (3)

а второй уровень - полиномом третьей степени

А(ак) = а0+ акак + актакат / 2+ahmakama„ 13;

W^B^SJl+B^SJJZS) (4)

где S = (SkSk)m \ak=S„IS\ акак = 1; (k,m,n = 1,2.....9).

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

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

^=QtaH)^;(iJ.k.m=1,2,3); (5)

№ Стк = АШк + Вша&; Сщ = AltJj + В,щ(аи + ам)\ i*j; (здесь и далее по индексам не суммировать);

Cll^Am + Buu<jau + ^2aIJ)\Cm = Avil + Bms[2av\yi С11А = + BiiJk(au + -j2aJky, j*k;

= Л*++; И;tem; Cijkm=сы/, Ajkm • Щкт ' компоненты тензоров четвертого ранга, определяющие физико-механические характеристики материалов из простейших экспериментов (в общем случае число независимых констант равно 42);

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

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

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

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

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

= (1 / е;+1 / е-у- о, 25[(i/е;+\/е; +

+1/Е;+1/ Е- - 20; / Е;+у- / Е;)\ ; (6)

В„ = V2(l / Е; -1 / Е-) - 0Д25л/2[(1 1 Е* +1 / £у+ -1 / EJ --ME-)- My], / е; - у- / £-)]; где Ef, vfj - модули упругости и коэффициенты поперечной деформации в направлениях соответствующих главных осей анизотропии; Щ - модули деформации в направлениях под углом 45° к соответствующим главным осям анизотропии.

Следует заметить, что результаты экспериментов могут представляться с использованием различных вариантов набора констант. Если вместо данных экспериментов по одноосному растяжению и сжатию под углом 45° к направлениям главных осей анизотропии имеются результаты испытаний на сдвиг в главных плоскостях х{х2, х2хъ, х]хг, то уравнения состояния для ортотропных материалов, связывающие касательные напряжения с деформациями сдвига можно упростить, представив следующим образом

еа = CwiTv где Ст = xnGv №1.2,3; И); (7)

{по индексам не суммировать).

Обоснование уравнений состояния также выполнялось путем сравнения экспериментальных диаграмм пропорционального деформирования ортотропного материала (стеклоткань - эпоксидная смола и графита AVCO Mod За) в направлении главных осей анизотропии и под углом 22,5°, 45° к ним с теоретически предсказанными по моделям A.A. Трещева, С.А. Амбар-цумяна и P.M. Джонса-Д.А.Р. Нельсона. При определении констант за основу были приняты данные экспериментов по одноосному растяжению, сжатию контрольных образцов вдоль главных осей анизотропии и результаты испытаний на чистый сдвиг в главных плоскостях анизотропии. Коэффициенты анизотропии при повороте осей координат вычисляются по обычным формулам преобразования компонентов тензора четвертого ранга, кроме этого пересчету подлежат нормированные напряжения.

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

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

Детально рассмотрена постановка задач, описана численная реализация и особенности расчета. Проведен анализ расчетов.

НДС прямоугольных пластин средней толщины

Рассмотрим упругое равновесие прямоугольной трехслойной пластины, отнесенной к декартовой системе координат (см. рис. 1). Внешние слои пластины одинаковы по своим свойствам и при этом в произвольной точке каждого слоя одна из плоскостей упругой симметрии параллельна срединной плоскости, а остальные две перпендикулярны к координатным линиям jc, = const, х, = const. Примем, что пластина нагружена нормально приложенной равномерно распределенной нагрузкой интенсивностью q. Вертикальную ось в сторону прогибов.

Е/,0 «№>.¡¡ = ¡2 —-i-!-Li-

/ Е/а> / \xî

Е,г<2) v/rZi

Рис. 1. Схема пластины

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

1) нормальное к срединной плоскости перемещение w не зависит от координаты х, (е, = 0) ;

2) нормаль к срединной плоскости после деформации поворачивается на угол <//, относительно оси и на угол i//2 относительно оси х2 ;

3) при определении параметров напряженного состояния влиянием нормальных напряжений ег3 пренебрегаем.

Деформированное состояние пластин определим компонентами перемещений точек срединной поверхности m, иг, ш = w.

Компоненты тензора деформаций выразятся через параметры деформации и кривизны % срединной поверхности (i, j = 1,2) согласно общим формулам

е.. =е..+х~х— (8)

У V У-и

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

и,,, + м-,, (9)

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

еп = )2 -*з</2„, еп = и2П +-^,2)2-х^,2 , (10)

у1г=и1,2+и1,1+и>,1м/,2-х3(м/1,1+1у1,2), Га=Гг + »п. Г23=^1 + и,.2-Для конкретизации структурной анизотропии материала пластин примем ортотропное тело. Тогда общие уравнения упругости для ортотропного разносопротивляющегося материала примем в соответствии с указанными выше гипотезами:

(11)

«11 ~ С Чш с М122 0 0 0

е22 с 1122 с 2222 0 0 0

е13 0 0 с И313 0 0

е23 0 0 0 ^-2323 0

,е12 0 0 0 0 ^1212

о-ц'

0-22

^23

где константы Ст имеют следующий вид:

~ А1Ш+ ВппСХп', Сп 22—^1122"

+ Втг(аи + а22); С22

2 — + В2222<*22 '

^1313 ^2323 — (^2323

в„,

^1212 (^1212 " где арст^/Б.

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

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

^114+^.2.2 = 0; ^2+^ = 0;

а„ +&>2 Н'=Н +^12.2 "'12 +Л'|2Ч +ип,2 = А/И!1 +Л/12,2 -а = 0; М22,2 +М12Ч -а = 0.

(12)

Решаемая задача изгиба пластинки имеет нелинейность двоякого характера, вследствие чего появляются определенные трудности при решении ее конечных разрешающих дифференциальных уравнений. Для решения поставленной задачи используется методика последовательных нагружений, предложенную В.З.Власовым и в последующем усовершенствованную В.В.Петровым в форме двухшагового метода последовательных возмущений параметров: Проведем лианеризацию с одновременным построением разрешающих уравнений. Пусть некоторое деформированное состояние пластинки будет начальным, а искомые функции будут определены на заданном уровне нагружения. Рассмотрим состояние, соответствующее некоторому малому возмущению нагрузки д + , которой будут соответствовать малые возмущения искомых функций + , и^Зи, и ч/1+8ц/1,у/г+ 8ц/г.

Опуская очевидные преобразования сразу выпишем уравнения статики в приращениях:

8Миа+8Мпп-8&= 0;

ЗМ11П +8Мпп -8<2г = 0;

¿21П +<*22>2 V,,, +<?ЛГ22,2 ЛГ11П

+ ^,22 ^22,2 +&У,а,2 +Щ211 М/,п +8щ,г Л\2,2 +Ы1т = ~8ц.

Для получения значений приращений внутренних усилий и собственно величин внутренних усилий на первоначальном этапе необходимо записать исходные геометрические соотношения в приращениях по вышеописанной методике, соответственно получим: 8еп = 8и1П +и<„ -хг(8у/1Ч), ^ = 8иг,г 8у>, 2 -х^{8у/1П),

8уп = 8щ, 2 +8и2Ч -ну,, 8м>,2 +м>,2 8мп +8™л <Ун>,2 -х, (<%,„ +8ц/гп), (14)

8уа=8ц/2+8ч>„\

ЗУп =81//1+8м>,1.

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

Зен §>|о> 8ап + 8еп да^ 8<У-а + деп дап 8ап + да, з 8аа + 2еи да2з 8агз

8еп = дап ■8ап + деп За-и ■8с22 + деп да12 ■8ап + ^22 даи ■8ап + &22 да,;

8е,2 деп 8ап + деп 8О22 + деп 8ап + де,2 8а,3 + де,, 8ап

д(Ти да22 дап дап 50*23

8е13 = Зеп дап 8ап + 8е,3 8ап 8ап + дап 8ап + де13 да,3 8а, з + 5е|3 дап 8аа

8вгз = де ц д<ти 8ан + &23 Эо"22 8а 12 + дап ■8а, 2 + дсгп ^13 + деп да23 ■8ап

или

8еп=Рт18ап + Рп228а22 + Рпп8ап + Р1П18аа + Рпа8ап = Рт18ап + Р22228а22 + РШ18ап + Рпп8ап + Рап8аи 5еП = Рццй*11 + Рп22^а2г + Рпп^п + + ^223 ■ (16)

<Цз = Рпп8ап + Р13228ап + Рип8ап + Рихз8ап + Рп138агз 5етз = ^зп^п + Р-аи5сгп + Ри\г8ап + Ртза^а + Рк-а8а-я

де

где = (и,к,т = 12,3). даи

Обращая указанные выражения получим:

[Cjfta]

8ern = Спп8еп + Cnl2Se12 + Cnn8en + Cln38el3 + Cila8ea 8an = Clm8en + C1222Sen + Cm28el2 + C22i38eu + C1222Sea (17)

Scrn = Ci2n8en+Cil228e22+CnnSen+Cm38et3 + Cn238ea 8an = Cml8en + CmlSe71 + Cml8en + Cim<5e13 + Cun8eu 8<t23 = C.im8en + C23228e22 + C23l28el2 +C23l3Sel3 +Cn23Se21 Приращения внутренних сил и моментов определяют традиционным образом, интегрируя выражения для приращений напряжений по толщине пластины, т.е. получим выражения для приращений усилий:

А/2 мг ыг

SNU = J Sandx3- SN12 = J Sa22dx3, SNn = J Scrndx3,

-h/2 -h/2 -hl г

h/2 h/2 SQl = J Srk3dx3, SQ2 = J Srk3dx3\ (18)

-h/1 -h/2 h/2 h/2 h/2 8MU = J 8апх^\ SM22 = J 8anx3dx3, 8Ma = J Sal2x3dx3.

-h/2 -h/2 -Ы1

Вычислив (18), учитывая (14) и подставив в (13) получим конечную линеаризованную систему разрешающих дифференциальных уравнений в приращениях перемещений. Окончательно система разрешающих уравнений имеет следующий вид:

8uim Dnn + ^8uinlDl2l2 + 8u2,n Dm2+^8u2,a Dlin -8щ,пGnn -8i//tnl Gm2 -

-8V/2,nGuu-8V/2,nGim=Fi(8w,8ul,8u 2,<Vi,M>w."i>m2>VpY/2)>

8ulnl Dml + ^8uimDml + Suln2Dnl2 + ^8ulm Dlin -8ц/„22 G2222 -S^,nGl2il -

-Sy/2n2 Glm - 8ц/г,п Gnn = F2(8W,8u1,8U2,8I//1,SV2,W,U1,U2,i//,,I//2),

-8wn Dm3 +8ulm Gim +^8щ,22 Gim + Su2,n Gn22 +-8u2n2Gl2i2 -

-SVmi Kmi -ЗУ^п Knil-Si//2,nKtltt-Si//2,22 КШ2 -|МД313 = <19) = F3 (Sw, 8u,, Su2, Stj/l, ¿>t//2, w, и,, «2, , У г), Sw,2 Dlm + SuV22 G2222 + ^8uml Gnn + Su2H2 Gim +^8и2Л2 Gim -

-8Ц/2Ч2 Km2 - 5Ц/2П2 K12l2 8Ц/у,22 ^цгг --<%,„ КШ2 ~8ч/,рпа = F^Sw^Su^SWiSw^w^u^y/^,

8wmD13l3 + 8w,22 D2323 + 2 D2за + ^!f2„ Din, -q =

= F,(8w,Зщ.Зщ,84/^8^,'»'.щ.щ,

Ау<™ ) = си<™ (х„х2,х])с1х3

^С^х^х^скз .

Правые части Р^.дщ.бщ^у/^г^.щ.щ.У^г) уравнений (19) имеют весьма

громоздкий вид, поэтому в тексте автореферата не приводятся.

Для полноты системы разрешающих уравнений (19) необходимо задать граничные условия. В частности для жесткого опирания они примут вид и, = 0; = 0; -и> = 0 и соответственно в приращениях 8щ =0; 84/,= 0; = 0, а для шарнирного опирания и, = 0; М„= 0; V = 0 и соответственно в приращениях 8и1 = 0; 8М„ =0; 8м = 0.

Начальные или накопленные к следующему этапу нагружения параметры системы обозначены буквами а 8^,8^,8^,8^ - текущие приращения, полученные

за очередной этап нагружения.

Вид переменных коэффициентов 0,]Ы, 0,]Ъп, Киш входящих в систему уравнений

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

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

Для апробации описанного выше подхода была решена модельная задача об изгибе жестко защемленной ортотропной пластины из слоистого композита (см. рис.1). Обкладки пластины выполнены из стеклопластика, а средняя часть пластины из графита марки Аи-Э [1] (см. табл. 1,2). Размеры в плане пластины 0,4x0,4 м; толщина пластины 0,06 м, толщина обкладок 0,01 м; равномерно распределенная нагрузка принималась от 0 до 650 кПа.

Материал Константы

Е,-МПа ЕгМПа ЕгМПа ЕгМПа Ез+ МПа Ез-МПа Ей* МПа Е,г-МПа

Стеклопластик 140000 70000 280000 140000 140000 70000 170000 875000

Графит Аи-Б 11850 10480 11850 10480 9450 7950 11040 9315

Таблица 2

Материал Константы

Е«* МПа Ei3-МПа МПа va МПа Кз МПа V* МПа V* МПа МПа

Стеклопластик 170000 875000 0,2 0,3 0,2 0,3 0,2 0,3

Графит ATJ-S 11040 9315 0,1 0,09 0,1 0,09 0,1 0,9

Для решения задачи разработан пакет прикладных программ на базе системы МАТ1.АВ рггооэ.

Было проведено сравнение полученных результатов с данными классической теории, причем разница по величинам перемещений составляла для различных величин нагрузки от 5 до 25 %. Для силовых факторов, в том числе моментов, разница колебалась в пределах от 6 до 40 %, и до 70 % процентов для максимальных напряжений.

На рис. 2-3 показаны графики изменения некоторых параметров напряженно-деформированного состояния пластины.

Величина прогиба W в центре плана пластины

- v/фЛЖ- - т - -

J&'f' \ \

- ----7 - -

0.02---J

—Ряекег по предяагммей мадепи «♦-С*-* Гйздсеичясюя пори» янинтдоми 5« ут» лодочных сдвигав С. А. А^бэрчУ*** P.M. Джонс-ДАР. Нельсон

Величина югруаки, Па

Рис. 2. Изменение прогиба W,mq центре пластины в зависимости от величины нагрузки

^ Напряжения ахх вдогъ диагонат пластины сверсу

—ф ~ Расчет по предлагаемой модели —О*- Классическая теория анизотропии

без учета поперечных сдвигов —GA. Аыбарцумян —В— P.M. Джонс -Д.А.Р. Нельсон

Расстояние вдоль диагонали пластины, м

Рис. 3. Распределение напряжений ах,<тх (сверху), Па вдоль диагонали пластины

НДС круглых прямоугольных пластин средней толщины

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

Решалась следующая прикладная задача об изгибе круглой пластины. Обкладки выполнены из стеклопластика, а средняя часть пластины из графита марки Аи-Э [1] (см. табл. 1, 2). Радиус пластаны 0,2 м; толщина пластины 0,06 м; толщина обкладок 0,0125 м; равномерно распределенная нагрузка принималась от 0 до 2 МПа.

Было проведено сравнение полученных результатов с данными классической теории, причем разница по величинам перемещений составляла для различных величин нагрузки от 4 до 20 %. Для силовых факторов, в том числе моментов, разница колебалась в пределах от 5 до 10 %, и до 60 % процентов для максимальных напряжений.

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

Рис. 4. Зависимость величины прогиба в центре пластины от величины равномерно распределенной нагрузки

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

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

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

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

Велтина прогиба УУ в цектрв плана пластины

Величина нагрузки. Па

Х10

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

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

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

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

4. Разработан алгоритм численного решения полученных систем дифференциальных уравнений на базе двухшагового метода последовательных возмущений параметров В.В. Петрова. Алгоритм реализован в виде программного комплекса в среде Matlab R2011а.

5. С использованием разработанного программного обеспечения решены следующие задачи по определению характеристик НДС:

- жестко защемленной по контуру и шарнирно опертой трехслойной квадратной пластины из ортотропных материалов: графита ATJ-S с внешними прокладками из стеклопластика под действием равномерно распределенной нагрузки, с учетом больших перемещений. Результаты расчета показали, что за счет учета эффекта разносопротивляемости и поперечных сдвигов удалось получить уточнение результатов, по сравнению с «классической теорией» и другими известными моделями в среднем до 25% для перемещений и углов поворота сечений и до 70% для максимальных напряжений;

- жестко защемленной по контуру и шарнирно опертой трехслойной круглой пластины из ортотропных материалов: графита ATJ-S с внешними прокладками из стеклопластика по действием равномерно распределенной нагрузки, с учетом больших перемещений. Уточнение результатов составило, по сравнению с «классической теорией» и другими известными моделями в среднем до 20% для перемещений и углов поворота сечений и до 60% для максимальных напряжений.

ОСНОВНЫЕ ПОЛОЖЕНИЯ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ В РАБОТАХ

(публикации в изданиях, включенных в перечень ВАК российских рецензируемых научных журналов, в которых должны быть опубликованы основные научные результаты диссертаций на

соискание ученой степени кандидата наук):

1. Васильев Н.В., ТрещёвА-А., Теличко В.Г. Моделирование напряженно-деформированного состояния гибких слоистых пластин из анизотропных разносо-противпяющихся материалов с учетом геометрической нелинейности II Известия Тул-ГУ. Технические науки. Вып. 2. Проблемы специального машиностроения. Тула: Изд-воТулГУ, 2011.-с. 541-547.

2. Васильев Н.В., Трещев A.A. Расчет напряженно-деформированного состояния трехслойной гибкой пластины из анизотропного разнсопротивляющегося материала с учетом больших прогибов II Строительная механика и расчет сооружений. Научно-технический журнал №6{239). М.: ЦНИИСК им. В.А. Кучеренко, 2011. - с. 42-48.

3. Васильев Н.В., ТрещёвА-А. Расчет напряженно-деформированного состояния трехслойной гибкой пластины из анизотропных разносопротивляющихся материалов с учетом больших прогибов II Известия ТулГУ. Технические науки. Вып. 2. Тула: Изд-во ТулГУ, 2012.-с. 179-188.

(публикации в остальных изданиях):

4. Васильев Н.В., Трещёв A.A. Напряженно-деформированное состояние гибких слоистых пластин из анизотропных, разносопротивляющихся материалов II IV-я молодежная научно-техническая конференция студентов Тульского государственного университета «Молодежные инновации»: сборник докладов. Часть 2 / под общей редакцией д-ра техн. наук, проф. Ядыкина Е.А. - Тула: Изд-во ТулГУ, 2010. - с. 33-35.

5. Васильев Н.В., Трещёв A.A. К расчету гибких слоистых пластин из анизотропных, разносопротивляющихся материалов II Современные технологии в машиностроении: сборник статей XIV Международной научно-практиктической конференции. - Пенза: Приволжский Дом знаний, 2010.-с. 271-272.

6. Васильев Н.В., Трещёв A.A. Поперечный изгиб гибких трехслойных пластин из анизотропных разносопротивляющихся материалов с учетом геометрической нелинейности II Проблемы современного строительства: сборник статей международной научно-практической конференции. - Пенза: Приволжский Дом знаний, 2011. - с. 33-37.

7. Васильев Н.В. К расчету гибких слоистых пластин из анизотропных, разносопротивляющихся материалов II XII Международная научно-техническая конференция с заочным участием «Актуальные проблемы строительства и строительной индустрии», Тула: Изд-во ТулГУ, 2011.-с. 5-6.

8. Васильев Н.В., Трещёв АЛ. Расчет напряженно-деформированного состояния трехслойной гибкой пластины из анизотропных разносопротивляющихся материалов II Новые идеи нового века - 2012 : материалы двенадцатой международной научной конференции ФАД ТОГУ: в 2 т. / Тихоокеанский государственный университет. - Хабаровск: Изд-во Тихооке-ан. гос. ун-та, 2012. - Т. 2. - с. 148-153.

Изд.лиц.ЛР № 020300 от 12.02.97. Подписано в печать 9.04.2012 г. Формат бумаги 60x84 'Лб. Бумага офсетная. Усл.печ. л. 1,1 Уч.изд. л. 1,0 Тираж 110 экз. Заказ 013 Тульский государственный университет. 300012, г. Тула, просп.Ленина, 92. Отпечатано в Издательстве ТулГУ. 300012, г. Тула, просп.Ленина, 95.

Текст научной работы диссертации и автореферата по механике, кандидата технических наук, Васильев, Николай Владимирович, Тула

61 12-5/2505

Тульский государственный университет

На правах рукописи Васильев Николай Владимирович

НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ ГИБКИХ СЛОИСТЫХ ПЛАСТИН ИЗ АНИЗОТРОПНЫХ РАЗН0С0ПР0ТИВЛЯЮЩИХСЯ МАТЕРИАЛОВ

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

твердого тела

Диссертация на соискание ученой степени кандидата технических наук

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

д. т. н., профессор Трещев A.A.

Тула - 2012

СОДЕРЖАНИЕ

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

1.ОБЗОР МОДЕЛЕЙ РАСЧЕТА ЭЛЕМЕНТОВ КОНСТРУКЦИЙ И МАТЕРИАЛОВ С УЧЕТОМ АНИЗОТРОПИИ И ЧУВСТВИТЕЛЬНОСТИ К ВИДУ НАПРЯЖЕННОГО СОСТОЯНИЯ ................................... 11

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

состояния ...................................................................................

1.2.Основные направления расчета многослойных

конструкций ..............................................................................37

2.ОПРЕДЕЛЯЮЩИЕ СООТНОШЕНИЯ ДЛЯ СТРУКТУРНО АНИЗОТРОПНЫХ РАЗН0С0ПР0ТИВЛЯЮЩИХСЯ МАТЕРИАЛОВ .................. 4 0

2.1.Вариант соотношений в девятимерном «пространстве», предложенный Н.М. Матченко,

А.А. Трещевым, З.В. Аркания ........................................4 0

2 .1.1. Ортотропное тело ...................................................45

2.1.2. Трансверсально изотропное тело .................48

2.2.Вариант определяющих соотношений для анизотропных разносопротивляющихся материалов ........ 54

2.2 .1. Общий случай анизотропии ............................... 54

2 . 2 . 2 . Плоскость упругой симметрии ........................55

2.2.3. Ортотропное тело ...................................................5 6

2.2.4. Трансверсально изотропное тело .................57

2.2.5. Определение констант уравнений состояния ..................................................................................57

2.2.6. Исследование принятых определяющих со-

отношений ................................................................... 61

3.ИССЛЕДОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ТРЕХСЛОЙНЫХ ПЛАСТИН СРЕДНЕЙ ТОЛЩИНЫ .....................7 8

3.1.НДС трехслойных круглых пластин средней толщины ...............................................................................................7 8

3 .1.1. Постановка задачи ................................................78

3 .1. 2 . Численная реализация .........................................8 6

3.1.3. Результаты расчета ..............................................8 6

3.2.НДС прямоугольных пластин средней толщины ......8 9

3 . 2 .1. Постановка задачи ................................................8 9

3.2.2. Результаты расчета ..............................................96

ЗАКЛЮЧЕНИЕ .................................................................................................98

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

ПРИЛОЖЕНИЯ .................................................................................................118

ВВЕДЕНИЕ

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

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

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

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

- б -

Для этой цели необходимо:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Внедрение результатов работы осуществлено в проектную практику проектно-сметного бюро на ОАО «ТУЛАОБЛГАЗ» (г. Тула) и ООО «Строительное проектирование» (г. Ту-

ла) . Использование результатов работы подтверждено актами о внедрении пакета прикладных программ.

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

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

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

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

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

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

1. ОБЗОР МОДЕЛЕЙ РАСЧЕТА ЭЛЕМЕНТОВ КОНСТРУКЦИЙ И МАТЕРИАЛОВ С УЧЕТОМ АНИЗОТРОПИИ И ЧУВСТВИТЕЛЬНОСТИ К

ВИДУ НАПРЯЖЕННОГО СОСТОЯНИЯ

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

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

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

- 12-

Модель С.А. Амбарцумяна.

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

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

„ _ орд^ , оРЯ0

(1.2.1)

где Б?? (податливости в системе координат pq) принимают

разные значения в зависимости от знаков главных напряжений :

при ор>0 и сгд > 0 : сРЯ _ оРЯ (1. 2 . 2)

при °р <о и сгд <0: оРЯ _ ЯРЯ 0íj (1. 2 . 3)

при °р >0 и а д <0: сРЯ сРЯ аРЯ оРЯ оРЯ оРЯ °11Ь' °12с' °22с' °61Ь' °62с' (1 .2 . 4

при °р <о и стд >0: оРЯ пРЯ аРЯ пРЯ оРЯ оРЯ °21с' °61с' °62Ьт (1. 2 . 5)

Податливости и связаны с податливостями

и в координатах главных осей анизотропии материала

обычными соотношениями теории упругости анизотропного тела: 511^ = 3111сС034(5 + + з±п2(5соз2(1 + 1п4/3,(1-2.6)

= Сс + С + - - Б™с) з1п2(1соз2(3, (1.2.7)

qpg 21tc

_ ^тп , , Qinn ~ ü21tc ^^

. QÜin _ 9 c'

lite + ü22tc ^ 21tc

,mn

- SgQtc) sin2¡icos2(5, (1.2.8)

qPq 22tc

S^tccos4p + (2Sj2tc + s66tc) sín2(3cos2(3 + S^tcsin4p,(1-2.9)

Sfftc =[S^tcsin^-S?>tccos^ + 0,5(2sr2tc + + s66tc) cos2/3cos22f3f3

(1.2.10

+ соз2(5соб22(5(5

Здесь нижний индекс Ьс означает, что каждый раз нужно брать соответствующие величины с индексом t или с; /3 — угол между главными направлениями напряженного состояния и главными направлениями материала (рис. 1.2.1) .

1.2.11

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

Податливости связаны с обычными «техническими»

постоянными формулами:

тп lite

Qmn _ qTtrn

b12tc ~ b21tc

тп /¡-.шш _ .nm /рлп

te ' te -~Vtc./ hte Г

mn 12tc

= 1/Ete,

qmn _ ? /r 66te ~ ±/ ^mntc'

(1.2.12)

(1.2.13)

(1.2.14)

(1.2.15)

где ~-£п/Ет, когда стш=сгс, а все остальные компоненты напряжений равны нулю. В соотношения (1.2.12 - 1.2.15)

входят восемь независимых постоянных материала: три — для растяжения, три — для сжатия и две — для сдвига.

Отметим, что в выражениях (1.2.1) урцФ0 , и, следовательно, главные оси напряжений не совпадают с главными осями деформаций. Кроме того, Амбарцумян связывал значения перекрестных податливостей и соответственно со знаками напряжений о и ад, причем условие =

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

координат) Б12ь совпадает с 312с только при чрезмерно жестких и физически нереальных условиях, наложенных на параметры ортотропного материала показанных в работе Исабе-кяна и Хачатряна. Следовательно, требование симметрии, не удовлетворяется. Это обстоятельство иллюстрируют диаграммы, показанные на рис. 1.2.2.

Рис. 1.2.2. Диаграммы податливостей