Научный вестникНГТУ. - 2013. -№ 4(53) МЕХАНИКА
УДК 539.3
Применение глобальных координат в модели составной осесимметричной оболочки при анализе ее статического и динамического поведения
В.Е. ЛЕВИН, А.Н. ПЕЛЬ, Д.А. КРАСНОРУЦКИЙ, П.З. АЛКЖАЕВ
В статье получены разрешающие уравнения расчета параметров напряженно-деформированного состояния оболочек вращения при статическом и динамическом осесимметричном нагружении. При выводе применяется глобальная система координат срединной поверхности оболочки и векторов поворота точки на этой поверхности. В качестве примера приведены результаты расчета статического деформирования составной оболочки, состоящей из цилиндрической и конической части. В частном случае для круговой пластины из общих уравнений получаются известные уравнения.
Ключевые слова: оболочка, орт, система координат, кривизна.
ВВЕДЕНИЕ
Теория оболочек в ее традиционной форме [1] располагает хорошо развитым аппаратом, позволяющим решать многие задачи практики. Для вывода разрешающих уравнений, как правило, используются локальные системы координат, опирающиеся на геометрию самой оболочки - касательные и нормаль к срединной поверхности. В то же время такая параметризация не всегда удобна, например, при расчете составных оболочек. При решении задач деформирования оболочек, имеющих излом меридиана предполагается склейка решений, что накладывает определенные неудобства при получении численного результата.
Представленные в данной статье разрешающие уравнения построены с применением глобальных систем координат и вектора конечного поворота [2]. Этот подход удачно применялся для описания поворотов поперечных сечений при деформировании тонкого стержня 13—41. Вектор поворота используется для описания деформирования окрестности точки на срединной поверхности оболочки.
1. КИНЕМАТИЧЕСКИЕ СООТНОШЕНИЯ
Для вывода уравнений оболочки предполагаем, что выполняются гипотезы Кирхгофа -Лява. Их применение для описания деформирования оболочки позволяет перейти к задаче деформирования ее срединной поверхности.
Запишем векторное уравнение срединной поверхности оболочки вращения в следующем
виде:
г (5, ф) = г(5) 1Г + I" , (1)
где - орты глобальной цилиндрической системы координат, 5 - координата вдоль мери-
Статья получена 26 июня 2013 г. Работа выполнена при поддержке интеграционного проекта Министерства образования и науки РФ № 7.822.2011.
Возьмем на срединной поверхности некоторую точку О с координатами s = s{) и tp = ср0. Затем линиями jj = rtp0 , jj = г (cp0 + Лр), s2 = sq , s2= s20 + ds2 вырежем малый элемент поверхности ¿/.vi x ds2 (рис. 2). Поместим тройку ортов локального базиса % , е2, е3 в рассматриваемую точку О.
Локальный базис, связанный с поверхностью, образуется следующими единичными векторами:
Г = г'ф = - sin ф ij + cos ф i2, r2 = r,s Vr + ",s X = r,s cos Ф V1 + r,s sin Ф i2 + ",s h , (2)
q3 = " Jr - r,s I = z,s cos Ф V1 + ",s sin Ф i2 - r,s h •
^V2
С q
Рис. 1. Геометрия и деформирование оболочки вращения
Рис. 2. Деформирование элемента срединной поверхности оболочки Связь ортов локального и глобального базиса выражается соотношением:
í
rj = VjÁ • Р =
sin ф cos ф О r s cos ф r s sin ф z s z s cos ф z s sin ф -r s
\
(3)
Деформацию малого элемента оболочки представим следующим образом:
- перемещение элемента параллельно самому себе как жесткого целого и поворот относительно точки О ;
- изменение длин его сторон;
- изменение угла между сторонами.
При повороте малого элемента поверхности относительно некоторой оси, проходящей через точку О, тройка единичных векторов, ориентированная так же, как орты /|, /2, г3 , по-
вернется и займет положение ¡¡*, ь , ц . Положение ортов в пространстве после поворота будет определено, если известен закон их преобразования [5]:
i,
2 yjkvk , (#)
где Xjj;. - элементы матрицы поворота.
Компоненты матрицы поворота вычисляются по формуле
sin |й| 1 - cos |й|
yij 2 8 ij + eijk iq ЙХ + ekinejmn-¡—¡2-ЙтЙк ,
Pi Й
здесь c5 = c5(ji,í2) - вектор конечного поворота [2]; a>¿ -проекции этого вектора на оси глобальной системы координат; - символы Леви-Чевиты.
Вектор поворота нормали в случае осесимметричного деформирования (рис. 1) имеет вид
Й 2 -Apq 2 -ApPi¿4 2 Ap sin w ii - Ap cos w v2 . (5)
Таким образом,
к>1 2 Apsin w , Й2 2 -Ap cos w, Й3 2 0 , 2 Ap . (6)
В результате деформирования дойны отрезков ОА 1 2 dsi и \ОА 2 2 ds
изменятся и ста-
нут равными соответственно
О* А1
2 ds* и О* А! 2 ds" .Деформации S[ и s 2 в направлени-
ds-i - dsi ds* - ds2
ях íj и í2: si =—1-L, S2 --•
ds- ds2
Рассмотрим вектор OA 2 до поворота:
dr 2 e2ds2 2 p2kVkds2 . (7)
При повороте относительно точки О координаты этого вектора в поворачивающемся базисе не изменятся. Длина может измениться в результате деформации растяжения-сжатия. В итоге получим выражение для этого вектора после поворота и деформации
dr * 2 e2ds2 2 p2i4*ds2 2 (1 + s2)p2iyk/jds2 . (8)
Вместе с изменением длины при деформировании возможно изменение угла между ортогональными векторами ёу , . Обозначим конфигурацию векторов, связанных с деформированной поверхностью, ёу*, е", (рис.1.). Векторы ёу* и е" выражаются через векторы
q* q*
и е2 :
e** 2 e! cos у+e* sin у, e** 2 e* sin y+e2 cos у,
где у - угол, характеризующий нарушение ортогональности векторов еу и е2 (для ортогональных векторов у = 0).
^ q ** q * q* „ q ** ( cos У sin У ^ „ 1 ( cos У - sin У ^
Обозначим ex 2 a1yej , em 2 amnen ,a 21 • I,h 2-—I • I.
J J ^ sin у cos у) cos2y^-sin у cos у )
Соотношение (8) может быть записано в виде
7— ♦ —Ф Ф 7 Ф —Ф j ♦ л "Г" Ф j Ф / 4 \ л а "Г 7
dr = e ^ = aljejdsl = Рд^ ^ = (1 + s1)Pjkyknvndí1-
^ dU dr* dr
Поскольку-=---, где и - вектор перемещении срединнои поверхности, то
dSj dSj dSj
^ = [ aij (1 + si) Р jky kn -P in ] Vn> v = 1 >2- (9)
Соотношения (9) описывают кинематику деформирования срединной поверхности оболочки и связывают поворот элемента оболочки, растяжение его сторон и изменение угла между ними с перемещениями элемента. Метрика деформированной поверхности определяется шестью функциями координат срединной поверхности оболочки: е> |(.V| ), ю 2(^1,^2) > ю 3(^,^2), е iC?i,) > s2(íl'í2) и Y(ji> J2) • Теоретически эти шесть функций могут быть выражены через параметры перемещений - шесть частных производных. Ограничений на рассмотренные функции не накладывается. Эти соотношения являются нелинейными. В случае малых поворотов и растяжений возможна их линеаризация.
Найдем кривизны оболочки в направлении нормали к ней в двух сечениях для недефор-мированного и деформированного состояний:
, de3 q * 5e3* ¿** 5e3 ¿ * 5e3* ¿**
k1 = — ei> k1 =—*ej , k2 = — e2> k2 =—7e2 - (10)
3í2 ds2
В случае осесимметричного деформирования оболочки вращения координатная сетка на оболочке остается ортогональной, т. е. у = 0 .
Кривизны недеформированной срединной поверхности имеют вид
k1 =~~e = —(-Zj sin w ¿1 + z^ cos ф q — sin ф q + cos ф q) = —,
к2 2 Г2 = ( С08 ф ^ + "ф ^ " ^ ) ( ^ ф ^ + ф ^ + "^ ) =
2 — • (11)
Кривизны координатных линий на деформированной поверхности вычислим по формулам
* де3* ^ * 11 де3* ^ * 11/ . . . \
к, 2-е, 2-----г 2—---(„ соеАр + г, втАр)
1 1 г (1 + 8,) 5Ф 1 г (1 + 8,)) Л >я 1>'
,* де" ^* 1 ЭЙ* 1 / N
к2 2 —е2 2"--д е2 2"-(( " Г , яя" , -АХ,, ) (12)
дз* 1 + 82 дя 1 + 82 '
здесь е* 2 Р^'п 2 РкпУтп*т •
Запишем приращения кривизн
* 1 1 А*1 2 *1 - *1 2 —
Г
7~—-Г - — cos + r,s sin )--- >
-1 + s1 — ' ' ' r
Ak2 = - к2 = —— (z г s - г ssz s - Ap) - Z>ssr s - r z s . ( 1 3)
1 )
Произведем преобразование выражений для производных перемещений (9). Так как U = Urir + Uziz = Ur cosф iy + Ur sinф ¿2 + Uz /3 , то из (9) получаем
дб 1дб 1 ^ - 1
— = -— = -Uг гф = sjz , отсюда: е\=~Ur. (14)
oíi г ^Ф г г
^ dU dUr - dU -
Поскольку-= —-ir h--- iz , то можем записать:
ds ds ds
dUг = (1 + s2)( cosAp-z>s sinAp)-rs , ( 1 5)
ds
dU. ds
2 ( 1 + r2)(z,s cos AP + ^s sin AP) - z,s . ( 1 %)
Таким образом, в случае осесимметричного деформирования оболочки вращения связь угла поворота нормали к меридиану, деформации меридиана и окружной деформации с радиальными и осевыми перемещениями точек меридиана выражается с помощью соотношений (14), (15) и (16). Если углы поворота нормали и деформации меридиана малы, то соотношения (13), (15), (16) можно упростить:
¿иг 4
-~Т 2 8 2 - АР >—Т~ 2 £ 2 + АР '
ds ds
* Ар * 1
Аку 2 ¿1 - ¿1 2 _2-, АХ2 2 ¿2 - Х2 2 -Ар ^ ,81 2 — иг . (17)
г ' г
2. УРАВНЕНИЯ РАВНОВЕСИЯ ОБОЛОЧКИ
Рассмотрим равновесие элемента оболочки, находящийся под действием распределенной нагрузки (рис. 3). По его граням прикладываются внутренние погонные усилия и моменты, полученные в результате осреднений напряжений по толщине оболочки.
Направление погонной силы относительно нормали к грани, в которой она действует, может быть в общем случае произвольным (возможны нормальные, поперечные и касательные усилия). Вектор погонного момента лежит в плоскости элемента для обеих граней.
При выполнении гипотез Кирхгофа - Лява и введении погонных усилий и моментов будем рассматривать равновесие элемента срединной поверхности оболочки.
Представим дифференциалы параметров поверхности в следующем виде:
dsl = А^а^, ds2 = dа2,
где А\, А2 - параметры Ляме, а |, а 2 - безразмерные параметры на поверхности.
Для деформированной оболочки имеем ds 1= А \clai, = А^а2 ■ Уравнения равновесия сил и моментов, записанные для деформированного состояния, выглядят таким образом:
1 э( а*7) х э( 4*72)
--^-'----^-'-+4 = О, (18)
А{ а2 3«! A¡A! За 2
5(А2М1) 1 3(А1*М2) 1 3 г* q 1 3 г * q
-¿ + —--^-'- + х Т1 + —— х Т2 = 0. (19)
АГ а! 3«1 аГ а! За 2 аГ 3«1 а! За 2
Приложенные к граням вырезанного элемента распределенные усилия и моменты могут быть разложены на компоненты в локальном базисе. Локальный базис связан с деформированной срединной поверхностью оболочки.
Тем не менее с точки зрения получаемых соотношений более удобным является использование проекций моментов и сил на глобальные оси системы координат. Получаемые при таком проектировании уравнения оболочки являются свободными от кривизн исходной поверхности, что позволяет описывать оболочки с изломами меридиана.
Рассмотрим осесимметричное нагружение оболочки вращения. Введем параметры Ляме в виде следующих соотношений:
ds1 = rdw , ds2 = id' ,
здесь = s / i , i - характерный геометрический размер оболочки, 5 - длина меридиана.
В случае деформированной поверхности
dsl = r (1+ S1 )dw, dsr = i (1 + s 2 )d' .
В силу осевой симметрии деформирования векторы неизвестных погонных внутренних усилий и моментов имеют вид
Т1 = T1w vw , Т2 = T2rvr + Т2zvz ,
H = МьЪ = M1srsir + Mbzsiz , М2 = М2ф^ф . (20)
Полагаем, что деформации малы. В этом случае геометрия деформированной оболочки практически не отличается от геометрии недеформированной, поэтому можно принять допущение, что Al = А1 = г , А2 = А2 = I.
Вектор распределенной нагрузки представим так: q = qrir + qziz .
Запишем уравнение равновесия сил (18) в проекциях на оси цилиндрической системы координат ir, iz :
^df-t - еТц+ qrri = 0, ^^ + qzri = 0. (21)
d' T d'
Уравнение равновесия сил в проекции на ось /ф удовлетворяется тождественно.
Перепишем уравнение равновесия моментов (19) в проекции на ось : < (М 2Ф)
С'
■+емъгг, + ге-72гг,,) = о.
(22)
Приведем необходимые соотношения теории упругости.
Нормальные к граням элемента погонные силы Щ, А'2 и изгибающие моменты М\, М 2, приложенные к граням элемента срединной поверхности оболочки (соответствующие векторы лежат в срединной поверхности), связаны с изменениями геометрии срединной поверхности следующим образом:
ей ЕЙ
11 2 т1ф =-2 (г1 +уе2)> 12 2 т2гг^ + 2-2 (г 2
* 1 -V2 1 -V2
[-V >3
I — V ,3
ей3 ЕЙ3
М1 2 Мк 2 --- (( + VД¿2), М2 2 —М2ф 2 —-- (9X2 + УАХ1). (23)
12(1 -V2)
12(1 -V2)
Введем безразмерные геометрические параметры оболочки - 2 - / е , г 2 г / е ,наряду с £, = .<;/( , и неизвестные функции (соответствующие индексы не указываются):
и 2 и / е, т 2 ^т, М 212(1 -Г) ей ' ЕЙ
Примем в качестве разрешающих следующие функции:
иг ,и" , ДР , т2г , т2" , М2ф .
М .
(2#)
(25)
Остальные величины, которые описывают напряженно-деформированное состояние оболочки, выражаются через разрешающие функции.
Для разрешающих функций получаем систему дифференциальных уравнений:
лиг _ ■2 г,'
С '
сЮ.
т2гг,'+ т2г",' -V~UГ г
Юг
<т2г 1
С' г
т" 21 с' г
^2ф
<' 2 [ т2гг,' + т2- ^^ | + г,'ДР ,
/ — — _ „ иг агг (1 -V )
(V- 1)г,'т2г +VT2+ (1 ^^ -Ч г Е- ' ^ ^ г Е-
-г,'т2" -
(1 -V2)
ЕИ
с Ар
М
С '
2(у-1))£М2ф+ -(V2 -1) ^ Ар--('-т2"7л)
С' И
- л2
2ф г,' . - VДp,
12, т
(26)
Систему уравнений (26) дополняем необходимыми краевыми условиями. Для получения системы уравнений осесимметричных движений необходимо добавить в систему (26) инерционные составляющие: дг = рк(£)иг, д2 = , где - масса
оболочки на единицу площади; точки означают дифференцирование по времени.
3. НАХОЖДЕНИЕ ПАРАМЕТРОВ НДС СОСТАВНОЙ ОБОЛОЧКИ
Проведем расчет статического деформирования составной оболочки, представляющей собой усеченный конус и цилиндр (рис. 4) со следующими параметрами: h = h / R = 0.005 , радиус цилиндра R = 0.5 м, высота цилиндра - Н = 0.7 м, длина образующей конуса L = 0.3 м, модуль упругости материала оболочки Е = 20 МПа . Нагружение осуществим постоянным давлением Р = Р • Е.
Рассмотрим нагружение оболочки с двумя разными углами конусности ф.
Численное решение уравнений (26) будем получать методом пристрелки. На рис. 5 и 6 сплошной линией обозначена деформированная конфигурация оболочки, пунктиром - неде-формированная.
Рис. 4. Схема нагружения
0.8
0.75
0.7
0.65
\
\
1
II 11
1.35 0.4 0.45 0.5
0.55
0.705
0.702
0.7
0.633
\ \\
\\ V
.492 0.495 0.493 0.5 0.503 7
Рис. 5. Деформированная конфигурация оболочки при ф = 30°, .Р=3х10
0.75
0.45 0.475 0.5 0.525 — 8
Рис. 6. Деформированная конфигурация оболочки при ф = 90° , Р = 1 х 10
При традиционном решении задачи о нагружении оболочки, составленной из конуса и цилиндра, производится стыковка решений для цилиндрической и конической частей оболочки. В рамках предложенной методики результат получен решением краевой задачи для системы дифференциальных уравнений всей оболочки, записанных в глобальных переменных, которые не привязаны к исходной форме оболочки и заранее удовлетворяют уравнениям совместности. Стыковка решений в этом случае выполняется автоматически.
В случае ф = 90° коническая часть составной оболочки вырождается в кольцевую платину. Покажем, что используемые уравнения применимы и для расчета круглых платин как предельного случая формы конической оболочки.
4. ПЕРЕХОД К ТЕОРИИ ПЛАСТИН
Из уравнений (26) можно получить уравнения для круглой пластины. Геометрию оболочки-пластины зададим следующим образом:
Т(я) = ^ ; "(я) = 0.
Система уравнений (26) распадается на две подсистемы, описывающие изгиб круговой пластины (27) и ее растяжение в плоскости (28):
dUz -АР Мщ 1 9 ¿Г2г 1
-- = Ар , -- =--=-*--УзАр, —— = 3
d' d' Н г d' г
- д2г (1 -У2) -и 2"--=-
2" ЕН
-М2„ _ 1 _ _ 2 1 12 -= (у-1)!М2Ф+ Н(V2 -1)4Ар + Н —2" .
d' г ^ у2 п
(27)
dUг — у иг
—7Г = т2г -у—, d' г
2г =1
d' г
(У- 1)—2г + (1 -У2) Ц2 - дгЩ-У1
г Еп
(28)
Сведем систему уравнений (27) к уравнениям теории оболочек. Для этого из четвертого уравнения (27) выразим осевую силу:
т = Н —2" = 12
-М2ф . . 1 _ - 2 1 -(у- 1)М2ф -н(у2 - 1) —Ар
d' г ^ г2
(29)
Подставим выражение (29) в третье уравнение системы (27):
12
d 2 М-
2ф / 1 -М2ш т. 2 .. ±
--(V-2)г -2--п(V2 -1)-
d'2 г d' -
1_-ар
2 d'
д2 (1 -У2) ен '
(30)
Выразим из второго уравнения (27) изгибающий момент и подставим его в уравнение (30):
( ,2а.Л
12
-3ар d
+ 2!
^ Ар
d '2
1 "ар
г2 d' -3
д2 (1 -У2) ЕН
(31)
Подставляя вместо Ар первое уравнение системы (27), получим
1 аи„
d' #
- + 2
1 d3U"
1 d2U"
г -'3
-2 лк2
г d'
г" d'
Д
Д = -
ЕН 3
(32)
Уравнение (32) соответствует уравнению Соф* изгиб круговой симметрично нагруженной пластины.
12(1 -у2)
Жермен-Лагранжа [6], описывающему
ЗАКЛЮЧЕНИЕ
Применение глобальных координат при составлении уравнений деформирования оболочки позволяет уйти от проблем стыковки при расчете деформирования составных оболочек (с изломом меридиана). Получены уравнения осесимметричного деформирования срединной поверхности оболочки с использованием глобальных координат. Сделаны проверки на частных случаях. Проведены расчеты статического деформирования осесимметричной составной оболочки. Влияние излома в меридиане оболочки на характер деформирования подобно влиянию шпангоута.
СПИСОК ЛИТЕРАТУРЫ
[1] Биргер И.А. Прочность, устойчивость, колебания: справочник в трех томах. Т. 1. / Под ред. И.А. Биргера и Я.Г. Пановко. - М.: Изд-во «Машиностроение», 1968. - 832 с.
[2] Аргирис Дж. Современные достижения в методах расчета конструкций с применением матриц / Дж. Арги-рис. -1968,-М.: Изд-во лит-ры по стр-ву. - 242 с.
[3] Пустовой Н.В. Применение геометрически нелинейных уравнений стержня к расчету статики и динамики тросов. Часть 1 / Н.В. Пустовой, В.Е. Левин, Д.А. Красноруцкий // Научный вестник НГТУ. - 2012. - № 1 (46). -С. 83-92.
[4] Пустовой Н.В. Методика вычисления параметров больших поворотов поперечных сечений гибкого стержня при расчетах в рамках его дифференциальной модели. Часть 1 / Н.В. Пустовой, В.Е. Левин, Д.А. Красноруцкий // Научный вестник НГТУ. - 2013. - № 2 (51). - С. 155-165.
[5] Левин В.Е. Механика деформирования криволинейных стержней: монография / В.Е. Левин, Н.В. Пустовой. - Новосибирск: Изд-во НГТУ, 2008. - 208 с.
[6] Уманский A.A. Строительная механика самолета / A.A. Уманский. - М.: Оборонгиз, 1961. - 529 с.
REFERENCES
[1] Birger I.A. Prochnost', ustojchivost', kolebanija - Spravochnik v treh tomah. Tom 1. Pod red. I.A. Birgera i Ja.G. Panovko. - Moskva, izdatel'stvo "Mashinostroenie" , 1968. - 832 s.
[2] Argiris Dzh. Sovremennye dostizhenija v metodah rascheta konstrukcij s primeneniem matric. -1968.- M.: Izd-vo lit-ry po str-vu. - 242 s.
[3] Pustovoj N.V. Primenenie geometricheski nelinejnyh uravnenij sterzhnja k raschetu statiki i dinamiki trosov. Chast' 1 / N.V. Pustovoj, V.E. Levin, D.A. Krasnoruckij //Nauchnyj vestnikNGTU, 2012. -№ 1 (46). - S. 83-92.
[4] Pustovoj N.V. Metodika vychislenija parametrov bol'shih povorotov poperechnyh sechenij gibkogo sterzhnja pri ra-schetah v ramkah ego differencial'noj modeli. Chast' 1 / N.V. Pustovoj, V.E. Levin, D.A. Krasnoruckij // Nauchnyj vestnik NGTU, 2013.-№2(51).-S. 155-165.
[5] Levin V.E., Pustovoj N.V. Mehanika deformirovanija krivolinejnyh sterzhnej: monografija. - Novosibirsk: Izd-vo NGTU, 2008. - 208 s.
[6] Umanskij A.A. Stroitel'naja mehanika samoleta. - M.: Oborongiz, 1961. -529 s.
Левин Владимир Евгеньевич, доктор технических наук, доцент, зам. заведующего кафедрой прочности летательных аппаратов Новосибирского государственного технического университета. Основное направление научных исследований - динамика конструкций. Имеет более 40 публикаций, в том числе 1 монографию. E-mail: levin@craft.nstu.ru
Пель Александр Николаевич, кандидат технических наук, доцент, доцент кафедры прочности летательных аппаратов Новосибирского государственного технического университета. Основное направление научных исследований - деформирование оболочек. Имеет более 20 публикаций. E-mail: pell@ngs.ru
Красноруцкий Дмитрий Апександрович, кандидат технических наук, доцент кафедры прочности летательных аппаратов Новосибирского государственного технического университета. Основное направление научных исследований - механика деформирования криволинейных стержней. Имеет более 20 публикаций. E-mail: dakras@mail.ru
Алюкаев Павел Зарифович, аспирант кафедры техническая механика Сибирской государственной академии управления, инженер ОАО «ИСС» им академика Ф.М. Решетнева. Основное направление научных исследований - модальный анализ ЛА. Имеет 6 публикаций. E-mail: pavel-alyukaev@yandex.ru
V.E. Levin, A.N. Pel, D.A. Krasnorutskiy, P.Z. Alyukaev
Usage of global coordinates in model ofa compound axisymmetric shellfor the analys is of static and dynamic behavior
At this paper the equations for calculating the parameters of the stress-strain state of shells of revolution under static and dynamic axisymmetric loading are derived. The global coordinate system of the shell's middle surface and rotation vectors at points on this surface are used to derive these equations. A natural curvature of a shell is not used in these equations. It is possible to solve nonlinear boundary-value problem for any profile line of shell of revolution. As an example the results of calculation of the static deformation of compound shell of cylindrical and conic parts are presented. In contrast to the traditional approach when we have tojoin the solutions for parts , the presented method allows us to obtain solution for compound shell immediately. The shooting method is used to solve nonlinear boundary-value problem for derived system of equations. In a particular case for the round plate the known equations follow from the common equations.
Key words: axisymmetric shell, large displacements, shells of revolution, ort, global coordinate system, curvature.