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

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

РГ8 ОД

НИРОДЖ ПАНДЕЙ

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

05.23.17 Строительная механика

АВТОРЕФЕРАТ

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

Москва - 1996 г.

Работа выполнена на кафедре Строительная механика Московского государственного строительного университета

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

профессор А.А.Амосов

Официальные оппоненты - доктор технических наук,

профессор Б.Ф.Власов кандидат технических наук, с.н.с. С.И.Трушин

Ведущее предприятие - Московский Государственный

Университет Путей Сообщения (МГУПС) Защита состоится 18 июня 1996 г. в I Ч час. Оп мин. на заседании диссертационного совета КОБЗ. 11.06 в Московском государственном строительном университете по адресу: Москва, Шлюзовая наб., д.8, ауд. N409.

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

Автореферат разослан " £ " О £ 1996 г.

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

Н.Н.Анохин

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

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

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

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

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

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

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

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

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

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

- на тестовых примерах уста* татав, получаемых по пред:.

- разработана методика расче и плит, основанная на си конструкции к однородной п

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

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

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

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

Апробация работы. Работа о£ нара кафедры Строительная мехаш Структура и объем диссерт?'

почения, содержащего сводку ос-литературы, включающего 78 наи-:ации составляет 125 страниц ма-ж 3 таблицы и 15 рисунков.

РЖАНИЕ РАБОТЫ

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

постановкой задачи расчета мно-контактной задачи сопряжения В.В.Болотина, Г.М.Амбарцумяна, Н.Москаленко, В.В.Парцевского, а также Я.М.Григоренко, В.Т.Ва-И.Андреева, А.А.Амосова и др. е некоторых основных работ дан-преимущества и недостатки дан-

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

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

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

В заключение главы на основании проведенного анализа формулируются основные цели и задачи диссертации.

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

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

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

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

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

с!Х(х)

= МХ + в (1)

с1х

-ь

где Х(х) - вектор искомых функций

Х(х) = X (б11, 612, б , U1. U2, U } (2)

^ Н Ü

И - матрица коэффициентов; В - вектор свободных членов.

Каждая из компонент вектора (2) в свою очередь представляет собой вектор, составленный из коэффициентов разложений функций в ряды по полиномам Лежандра

ff6ll.612.6 ,ui,u2,u) - if<°\f(1>.f(2>...f<N>>T (3) k 14 t,l

где fk =(k+l/2) | f (xi,X2)PkU)d£ ; (4)

4 = h/R ; -1 < С, < 1

h - полутолщина оболочки

PkU) - полиномы Лежандра k-го порядка.

В развернутом виде система уравнений (1) записывается следующим образом:

1

бц.1 = шцб12 + И126 " -V 1М - Fi ;

14 h 14

1

612.1 = m2l6u+m22U2+m23U +m24U2*+®254--V~% -F2 *,

4 4 h 24

1

б =шзабц +m32U2 + тззи + тз^г* + nfcsU--V-1M -F ;

14,1 4 4 h 14 4

111.1 = М41бц + Ш42И2 + И431^ + П144и

иг. 1 = 4151612 + ®52Ч1 ;

и = Иб1б + Ш62Ч1 + Ибзи1 ; 4,1 14

Коэффициенты этой системы уравнений определяются следующими формулами:

т 1 1 1

тц =-V-1; Ш12 --V-1 N1

I? Ь

Ло

т21= -т-V-1; т22---г-Г1^! - Н)У~1и2 +т2-~-Г1Г1;

И

«423 =--- ии + -ЮЗ.

№

фа фо

Ш24 ---- ЮГ1 ; Ш25 = -Ш - V-1 ;

№

1 *

таза ---- фгН) ;

I" 1 л л

тзг ---УЧфвУЧГг' + (ф3*1 + 4-4») V"1];

ИЪ

Изз=——-V-1 С(ФзА21 - «ЦАЮУ1 -(ф4А1 +Ф5Н)И]+ Ь

+8 1 = -^Г1 5 435 = Г1 (+4*1 - ;

- и -

4142 = - V"1 ; ГП4Э =

R

(Ф1ЛГ1 - ф2и);

ТП41 = (|>gl

h

Ш44 =

»2 h

m

sr 1

1451= +9 I

®52 =

R

1

®61 = Ф7 I

тез =

h

I

Здесь V, N, ü и Ü2 - квадратные матрицы N-го порядка, определенные в диссертации;

4>i - коэффициенты, зависящие от физических свойств материала;

М и Ui - векторы, определяемые видом заданных it

граничных условий на лицевых поверхностях оболочки, соответственно, статических или кинематических;

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

f (бц, б , Ui, и } = f(x) sin Шф 1 Н V

f Í6i2, U2> - f (X) COS П)ф

Введенные обозначения напряжений и перемещений представлены на рис.2. Индекс после запятой обозначает дифферен-

- IS -

дарование по соответствующей координате.

Далее рассматривается трехслойная оболочка со слоями заданной толщины V Л и V, где штрихи обозначают соответственно относительные толщины верхнего V и нижнего слоев У' ( X = h/R). Предполагается, что каждый из слоев определяется соответствующей матрицей жесткости материала. Показано, что формулировка контактной вадачи сопряжения слоев связана с построением матриц сопряжения Кь которые определяют вид векторов

М ', М , М " , UiVUi'-и Ui" и вависят от заданных усло-i 4 14 14

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

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

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

d X' Hit *t2 0 X' В'

- X SS Mai Mzz Ü23 X + в

dx X" 0 Ü32 Мзз X" В"

При этом в рамках разработанной методики для каждого из слоев может быть задан свой собственный порядок приближенной теории - И', N. Ж*. Таким образом, суммарный порядок системы (7) составляет 6 [ (Н'+1)+Ш+1)+Ш"+1) ].

Для каждого из слоев могут быть заданы соответствующие

граничные условия, определяемые условиями закрепления слоя.

СоХ

= 0; Сь X

х=0

- 0 (8)

х=1

где Со и Сь - некоторые прямоугольные матрицы, определяемые граничными условиями.

В совокупности для каждого из краев слоистой оболочки описанным образом может быть составлено 3 [ (Ы'+1)+(М+1)+ +(N"+1) ] граничных условий.

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

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

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

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

многослойной оболочки, могут быть разложены в сходящиеся ряды по полиномам Лежандра, причем в силу теоремы Нейманна для кусочно-непрерывной функции этот ряд имеет сумму 1/2[Г (4+0)+ + £ (4-0)3 в любой точке с, промежутка -1 < с,< 1.

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

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

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

Е б(к)(с*1,«2) РкСО « 0(0 ,П!к)(«Ь«2) Рк<0 (0) к=0 к=0

где бс1° и еС1° - векторы, составленные из коэффициентов разложении компонент тензоров напряжений и деформации:

t((К)^ (Ю^ сю}г НО { (к)\

^ и г с, а ' (10)

0(4)- симметричная матрица упругих постоянных, элементы которой рассматриваются как кусочно-непрерывные функции координаты 4 . Конкретный вид этих функций определяется упругими свойствами материала слоев.

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

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

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

,С0)

»Ш

010)+.

(И)

где 0ск5 представляют собой разложения матриц Б(4) по полиномам Лежандра.

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

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

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

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

В качестве тестового примера была выбрана многослойная цилиндрическая оболочка согласно рис.3.

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

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

С2 ---; (12)

к!(1-У1) (1+У1)(1-2У1)

где к2 и кз. - коэффициенты, определяющие параметры

------- --------------------модуля упругости слоев;

- vi и V2 - соответственно, коэффициенты Пуассона материала слоев; значения элементов матриц жесткости С Сг = 1 * £0) незначительно отличаются от матрицы жесткости, определяемой элементами г(0) - «(0) + 2в<0) И 8(0), где

а(0) = ь _ . а(0) = __м3ч

(l+vf0))(l-2v(0)) 2(1+ v(0)5 UJ

v(0) и g(0)_ соответственно значения обобщенных коэффициентов Пуассона и модулей упругости Юнга, определяемых по формуле

f(v(0\Et0>) = — ! fi fi+1dА (14)

2 1=1 Ai

где S - число слоев - рис.3.

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

Представленные на рис. 4-5 результаты расчета цилиндрической оболочки, подверженной внешнему давлению, позволяют заключить, что при х=0,1 расхождение результатов между данными, полученными по предлагаемой теории ( сплошная линия), и данными, получаемыми по теории С.А.Амбарцумяна, использующей гипотезы Кирхгоффа-Лява, не превосходит величины порядка 15%.

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

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

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

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

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

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

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

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

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

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

Рис. 2. Напряжения и перемещения в оболочке

Рис. 3. Схема пакета многослойной оболочки

iO

2,0

3,0

-to -2,0 -5,0

-5,0

\ //

V \\ 1 !

\\ \ \\ > ■■04 1 i i

/ и

----• ——.

Рис. 4. Нормальные перемещали в оболочке

Ч

Рис. 5. Осеше напряжения в оболочке