Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Обобщенная теория анизотропных оболочек

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
12.55 Mб
Скачать

Обобщенная теория анизотропных оболочек / Хо­ йя И. Ю.— Киев : Наук, думка, 1986,— 172 с.

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

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

Ил. 25. Табл. 1. Библиогр.: с. 163— 168 (140 назв.).

Ответственный редактор Я. Ф. Каюк

Рецензенты К- И. Шнеренко, А. П. Мукоед

Редакция информационной литературы

X

1703020000-275

194‘...........

 

М221(04)-8в

© Издательство «Наукова думка», 1986

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

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

Существует ряд подходов построения теории оболочек и пластин при менее жестких допущениях, чем гипотеза нормального элемента. К ним относится теория типа С. П. Тимошенко, теория Э. Рейсснера, С. А. Амбарцумяна, Т. В. Койтера и др. Некоторые их приложения к задачам статики, динамики и устойчивости оболочек и пластин даны в монографиях С. А. Амбарцумяна [61, В. В. Болотина, Ю. Н. Новичко­ ва [12], Г. А. Ванина, Н. П. Семенюка, Р. Ф. Емельянова [15], А. Н. Гузя, И. С. Чер­ нышенко, Вал. Н. Чехова, Вик. Н. Чехова, К. И. Шнеренко [49], Б. Л. Пелеха [88],

атакже в работах А. Т. Василенко [14], А. П. Мукоеда [81], Ю. Н. Новичкова [85],

К.И. Шнеренко [128, 129]. Работы А. К. Галиньша [34], Э. И. Григолюка, Ф. А. Ко­ гана [44] носят обзорный характер.

Имеющиеся здесь подходы условно разделяются по классификации-И. И.-Воровн- ча [28] и А. Л. Гольденвейзера [43] на две группы. К первой относятся методы, содер­ жащие регулярный процесс замены решения трехмерной задачи последовательностью решении двухмерных, ко второй — методы, базирующиеся на принятии упрощаю­ щих гипотез без регулярного процесса.

Впервой группе выделяются два направления, связанные-с асимптотическими методами и методами разложения по толщине. Асимптотический способ'интегрирова­ ния уравнений трехмерной теории упругости описан в работах А. *Л. Гольденвейзе­ ра [42], К. О. Фридрихса [136], Э. Л. Рейсса [139], Л. А. Аголовяна [1], В. Л. Бер­ дичевского [8], В. А. Колос [65] и др. Иной подход к вскрытию исимптотики решений, основанный на однородных решениях, развивается в работах И. И. Воровича [29], О. К. Аксентяна, Т. В. Селезневой [5], А. Н. Ведерникова, Т. В. Виленской [19], Т. В. Виленской [26], О. К. Аксентяна [3], А. С. Космодамианского, В. А. Шалдырванз [66], Ю. А. Устинова [98], М. А. Шленева [127].

Метод разложений по толщине основан на представлении искомых функций в виде рядов по положительным степеням толщинной координаты. А. Коши и С. Пуас­ сон рассматравали задачу упругого равновесия пластины. Впоследствии Н. А. Киль­ чевский [64] этот метод распространил на задачи теории оболочек. Причем вместо

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

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

и оболочек описан в работах

А. И. Лурье 74J,

В. К. Прокопова [92], Ю. А. Груздева, В. К. Прокопова [47],

У. К. Нигула [84],

И. Т. Селезова [96].

и В. В. Понятовского [90, 91]

В работах И. Н. Векуа [20—23], П. Чикалы [132]

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

Лежандра. При этом для со­

ставления соответствующих граничных и начально краевых задач относительно коэф­ фициентов разложений (как функций координат срединной поверхности и времени) используются проекционный и вариационные методы. Проекционным способом в [20) И. Н. Векуа получены уравнения равновесия (движения) призматических оболочек переменной толщины. Дальнейшее развитие его и обобщение на произвольные тонкие пологие оболочки описано в работах [21—23]. При составлении уравнений равнове­ сия оболочек вращения П. Чикала [132] использует вариационный принцип возмож­ ных перемещений. В. В. Понятовский [90, 91] строит уравнения равновесия однослой­ ных и многослойных ортотропных пластин при помощи полуобратного метода СенВенана с применением вариационного принципа Кастильяно. В работах Ш. К. Гали­ мова [33], В. И. Гуляева, В. А. Баженова, П. П. Лнэунова [52], О. Ю Калекина [62], А. Солера [97, 99, 137], В. К. Чибирякова [126] и автора [116J эти методы распростра­ няются на нетонкие изотропные и анизотропные пластины п оболочки.

Изучению оператора теории оболочек И. Н. Векуа и решению граничных задач о напряженно-деформированном состоянии изотропных пластин и тонких пологих оболочек посвящены работы Т. С. Вашакмадзе [16— 18], Д. Г. Гордезианп [36—39|, Д. Г. Гордезиани, О. П. Комурджишвили [40, 41], А. Ф. Гюнтера [53,54], Г. В. Джа-

ианн [55], В. С. Жгенти [59, 60], Л. С.

Кнкнадзе |63|,

Т.

В.

Меупаргин

[75,

76|,

А. Р. Хволеса [100] н автора [101, 102,

104). В работах

[60,

61,

105-113,

117,

1221

излагается метод построения общих решений уравнений равновесия оболочек и плас­ тин. Обобщение метода И. Н. Векуа на анизотропные оболочки дано в |68, 77, 103,

112— 121]. Слоистые оболочки

и

пластины исследовали Н. К. Галимов [31, 32],

Б. Л. Пелех, В. А. Лазько [89],

В.

Е. Чепига [123, 124].

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

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

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

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

ГЛАВА 1

ОСНОВНЫЕ УРАВНЕНИЯ ОБОБЩЕННОЙ ТЕОРИИ АНИЗОТРОПНЫХ ОБОЛОЧЕК

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

§ 1. НЕКОТОРЫЕ СВЕДЕНИЯ ИЗ ДИФФЕРЕНЦИАЛЬНОЙ ГЕОМЕТРИИ И МЕХАНИКИ ДЕФОРМИРУЕМОГО ТВЕРДОГО ТЕЛА

Обозначим через R" евклидово пространство размерности п. Пусть V" — область из Rs, занимаемой упругим телом. Под телом (как обычно

в механике сплошной среды) подразумевается материальный конти­ нуум, непрерывно заполняющий данный объем *?/.

В дифференциальной геометрии свойство непрерывности играет важную роль в определении такого понятия, как дифференцируемое многообразие [9, 57]. По определению множество М называется п-

мерным многообразием, если для каждой

точки £ £ М существуют

связная открытая окрестность U, £ € U с

М и некоторый гомеомор­

физм ф : U -*■ ф (£/), переводящий U в открытое множество <р (U) — G

из п-мерного пространства Rn над полем

вещественных чисел R.

Если £ — точка из М, £ £ Uc, где i £ / (/ — некоторое фиксирован­ ное множество индексов), то координаты точки фг (£) £ Rn называются

локальными

координатами точки

£ £ М. Пара

(Uti qpf) образует

локальную

карту многообразия М,

а семейство

всех

локальных

карт —атлас.

 

 

базис в R3.

Пусть ik (k = 1, 2, 3) — некоторый фиксированный

Тогда любой точке £ £ *V ставится

в соответствие три вещественных

числа х*, называемых координатами точки £ относительно отобра­ жения ф,

ч> s e v + c p d )= **(?)< *s я 3.

Выбор

системы координат на V' имеет некоторый произвол. Если

(и {,

ф,) и

(£//, ф,) — две локальные

карты на *1/ и

Ut П £// Ф 0 . то

для

точки

£ £ ifi П Of

определены

две системы

координат: х* (£)

и хк>(£). Они связаны на

U{ П £7/ при помощи преобразования коор­

динат (рис. 1)

 

 

 

Ф/т: (*й) б % {Ut П Ui) -+■ x k’ = ф/ о ф Г 1(**) € Ф/ (U, П U /)•

Подробное

изложение дифференцируемых многообразий можно найти

в [57].

'

 

1.

 

Определение поверхности. Пусть Й — открытая ограниченная

в Я*I2 область с границей Г. Тогда поверхность S евклидового простран­ ства Я3 будет образом й при отображении г (рис. 2)

г£2) £ й - > г ( |\ i 2)£ .S c : Я3.

Данное отображение, сопоставляющее каждой точке £ = (£\ £2) £ Й радиус-вектор г (£ \ | 2), будем считать достаточно гладким. Эго озна-

 

 

 

Рис. 1.

 

 

 

чает,

что

г (I1, £2) — дифференцируемая

вектор-функция

параметров

I1,

| 2,

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

5.

 

 

Обозначим через аа базисные векторы на S:

 

 

 

 

 

&а — Го |г,<х = ~ЩаГ

а =

 

(1Л)

и по ним найдем единичный вектор внешней нормали

 

 

 

 

 

Яд —

fll X а2

 

 

 

 

 

3

|аа Х а 2|

 

 

где

знаком | • | обозначена норма евклидового пространства

Я3, а

X — векторное произведение. Наряду с ковариантным базисом аа (а =

=

1 ,2 )

рассмотрим контравариантный

(биортогональный)

базис

а^(Р =

1, 2), удовлетворяющий

условиям

 

 

 

 

 

аа •

= б£,

 

 

 

где 6а — символ Кронекера.

 

 

 

 

Векторные произведения этих векторов приводят к таким соотно­ шениям:

аа X ар =

с«ра3;

а® X а& =

с®Ра3;

а3 х ар =

£pvav;~

а3 X аР =

c&vav.

Здесь сар — дискриминантный тензор

V о.

€аЗ» 6а& — символы Леви — Чивита

[€«el = [£““| = [ _ ° ‘];

а = det

По повторяющимся индексам (как обычно в тензорном анализе' подразумевается суммирование, причем греческие буквы принимаю! значения I, 2, а латинские— 1, 2, 3.

Обозначим через аар компоненты первой квадратичной формы (или

метрического тензора) поверхности

ЯаР = Два = аа • ар = Г,а Гр.

(1.2)

При помощи тензора а«р вводится метрика на поверхности S, т. е

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

da = V a d ^ d l2.

(1.31

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

^ ]a , b[<=R-+ta = xa (t).

Тогда линейный элемент этой кривой имеет вид

d&=Va„»(x*Y (x»)'dt,

где (*v)' = dxVJ !) .

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

^аР = ^Ра =

^3,а * Яр = а3

.

/ а

OXL

Р

О * ’ )

ьр = а

&рх =

— аз.а -ар

 

Здесь а“х — компоненты контравариантиого метрического тензора. Дифференцируя базисные векторы аа, а* и а3 = а3 по координатам

£0, получаем вектор-функции

 

Яа,р — Гccf'Av -(- Ьа$а3 — Газ,х.а -f* £\хря3;

(1.5)

 

=

— Tpvav -f- b$а8;

(1.6)

 

аз,0 =

а3(з = — Ьрav,

(1.7)

где Гард и

= flvXrap,\ — символы Кристоффеля

первого и вто­

рого родов. Формулы (1.5) и (1.6) называются деривационными фор­ мулами Гаусса, а (1.7) — формулой Вейнгартена 1241.

Поверхность с заданными коэффициентами первой и второй ква­ дратичных форм будет неразрывной, если удовлетворяются уравнения Гаусса [241:

Rva&y == д^сср.у

dTgy.y

+ ГдуГXp>V

ГдзГ^у.у RCavCfrу (1*8)

 

а£0

 

 

и уравнения Петерсона — Кодацци

 

 

^y^ap

^Фау 0»

(1*9)

где Va — символ ковариантной

производной;

К = b\b\ b\b\>

Возьмем на поверхности S произвольную точку Р и вектор s, ле­ жащий в касательной к S плоскости и приложенный в точке Р. Плос­ кость, проходящая через s и вектор нормали m = а3 в Р, пересекает

поверхность 5 по кривой, определяемой отображениями

s £ ]й, Ь[ с : = xa (s).

Нормальная кривизна ks этой кривой выражается формулой [131]:

ааар(х*)' [х*У

Кривые, для которых ks принимает минимальное и максимальное зна­

чения, называются главными. Обозначим главные значения кривизн соответственно через kx и /г2.

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

1

- а Ь

( 1. 10)

аи -

 

Ьц k^A\\

где

м II

дт 1 36, | ~

дт

 

 

 

 

 

 

д12- I -

 

 

 

 

 

 

 

 

 

 

 

 

а символы

Кристоффеля

выражаются таким образом:

 

 

 

 

 

 

 

ЭЛ,

 

 

 

 

 

 

 

 

дАг

 

 

 

 

ЭЛ,

 

 

Гил =

А \

36,

 

 

1

Г 12,1 =

Г«ц,| = Л|

36а

1

Г22.1 =» - - 36,

1

(1.11)

 

 

 

 

dAt

 

Г12,2 = Г21,2 == Л2.

ЗЛ2

 

 

эл2

 

 

 

 

 

 

 

 

 

 

Г11,2 =

 

А - 3|2

I

36, • ;

Г22,2 :=

А

36а

 

 

и

 

1

дА 1

 

 

 

 

 

 

ал,

 

 

 

 

^2

зл2

,

 

Г',1 -

t

г |2 =

14 =

-д -

 

 

Ги =» —

 

А\

36,1

 

3|а

 

 

 

36, ’ ( 1.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г? ,_

 

А,

dAt

■*

Г?о

Го|

 

-

1

эл2

-

Гоа

1

ал2

 

 

1 и —

 

 

 

36а

^2

36,

1

J.22

А2

36а

 

 

 

 

 

 

 

»

* 12

1 21

 

 

 

Согласно (1.10) — (1.12)

уравнения

Гаусса

(1.8) и Петерсона — Ко­

дацци

(1.9) преобразуются

к виду

[87,

 

I25J

 

 

 

 

 

 

 

э

 

1

 

эл2

 

э

 

1

ал,

 

 

 

 

 

 

(1.13)

 

 

36, 1 Ах

agt

! ^

З6а

 

 

£-) = -

 

 

 

 

 

 

 

см

МГ-

см

 

 

 

 

 

 

 

 

 

 

 

^2^2 _ U

дА2 .

ЭЛ,/г,

 

ka

ЭЛ,

 

 

 

(114)

 

 

 

 

 

35,

 

 

36,

*

 

36а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2.

Определение геометрии оболочки.

Под оболочкой

обычно под­

разумевают трехмерный криволинейный слой, толщина которого на­

много меньше

двух других его

 

 

линейных размеров. Внутри слоя

 

 

между

его

граничными поверх­

 

 

ностями 5+

и

S- располагается

 

 

основная (координатная) поверх­

 

 

ность S. Расстояния между S и

 

 

S+, S~

обозначаются соответст-

 

 

 

 

с-Н

(—>

 

 

 

 

венно через

h

и h.

В

общем

Рис. 3.

 

случае

полутолщйны

(+)

(-)

 

 

h,

h яв-

(+)

(—)

ляются функциями координат

$3 £ 5. Если h (I1, £2) =

h (£\ £а),

то оболочка обладает симметрией относительно координатной поверх-

ности 5 и она

(+)

(—)

становится срединной поверхностью. При h =

h =

= h (h = const)

имеем оболочку постоянной толщины (рис. 3).

 

Отнесем оболочку к пространственной системе координат, нормаль­ но связанной с ее основной поверхностью S. Это означает, что семей-

ство координатных линий I1, | 2 расположено на поверхности 5, а семейство координатных линий | 3 направлено вдоль нормалей к S. Положение любой точки оболочки в данной системе координат опре­ деляется радиус-вектором

 

R =

г 4* £3а3.

(1.15)

Из

равенства (1.15) путем

дифференцирования

по координатам

£•' (i =

1 , 2 , 3) получаем базисные векторы

 

ga = j4a.Bv» g3 ~ a3,

где А'а. — тензор вида

А £ = а 1 - 1 * Ы

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

 

 

бар =

Aa.-A$-G\vt

баЗ =

0,

бзз “ 1*

 

( 1. 16)

Введенную таким

образом

параметризацию

слоя будем

называть

S g-параметризацией

[24).

 

 

 

 

 

 

 

 

 

Если система координат I1, £2 совпадает с линиями

кривизны, то

5 в-параметризация

 

оболочки

переходит

в

S

^ - параметризацию.

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

 

 

8 и =

^1*

622 =

^2»

бзз =

 

612 ^

бм = 823 =

( 1- 17)

где Hf — параметры Ламе;

 

 

 

 

 

 

 

 

 

H

^ A ^ l - k &

y ,

Н2 = А2(\ - й

2У ;

/ 7 , - 1 .

(1.18)

Элемент объема

представляется

формулой

 

 

 

 

 

 

 

 

 

 

dv = Bdadl3.

 

 

 

 

 

(1.19)

Здесь

da =

V ad^dl^,

0

=

 

= 1 — 2Н£3 -f

Kll,

причем Н и

К — средняя

и главная

кривизны

поверхности S.

 

 

3.

Деформация

оболочки.

В процессе

деформации

произвольная

точка М оболочки испытывает некоторое

перемещение. Пусть и (£\

I 2,

0 — вектор

перемещения.

Тогда

положение точки

М после

деформации определяется

радиус-вектором

 

 

 

 

 

 

 

 

 

 

 

R* = R 4- и.

 

 

 

 

 

(1.20)

Отсюда находим базисные векторы

 

 

 

 

 

 

 

 

 

ga = ga +

U,a

(a =

1, 2);

g3 = а3 +

u„

 

 

и компоненты метрического тензора

 

 

 

 

 

 

 

 

gcj =

ga + g/ • u,/+ g/ • ur/+

 

u,/.

 

(1.21)

Деформированное состояние среды измеряется разностью квадра­ тов линейных элементов после деформации и до деформации, т. е.

ds'2 — ds2 = 2eifd£dg,

(1.22)

ю