книги / Обобщенная теория анизотропных оболочек
..pdfОбобщенная теория анизотропных оболочек / Хо йя И. Ю.— Киев : Наук, думка, 1986,— 172 с.
Монография посвящена граничным задачам обобщенной теории нетонких анизотропных оболо чек и пластин. Редукция трехмерной задачи к двух мерной осуществляется при помощи метода разло жений по толщине с использованием полиномов Лежандра. Вариационным способом составляются линейные н нелинейные уравнения движения не тонких и тонких анизотропных оболочек, начальные и естественные граничные условия. Описываются методы построения общих регулярных решений уравнений равновесия трансверсально-изотропных оболочек и пластин, а также построения эквивалент ных уравнений. На основе данных решений рассмат риваются задачи концентрации напряжений около отверстий и строятся фундаментальные цртрццы ре шений. Рассматривается вопрос о- влиянии' попе речных модулей упругости и сдвига на напряженное состояние оболочек, частоту колебаний и величину критической силы.
Для научных и инженерно-технических работ ников, специализирующихся в области твердого де формируемого тела.
Ил. 25. Табл. 1. Библиогр.: с. 163— 168 (140 назв.).
Ответственный редактор Я. Ф. Каюк
Рецензенты К- И. Шнеренко, А. П. Мукоед
Редакция информационной литературы
X |
1703020000-275 |
194‘...........8в |
|
М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) |
ю