CC BY
9Математические модели деформирования оболочек переменной толщины с учётом различных свойств материалов
К.т.н. В.М. Жгутов*,
ООО «Архитектурно-строительная компания «Китеж»
Ключевые слова: оболочки переменной толщины; ребристые оболочки; ортотропия; линейная и нелинейная упругость; вязкоупругость; поперечные сдвиги; полные углы сдвигов; переменная жесткость; сдвиговая и крутильная жесткости ребер; функционалы полной энергии деформации оболочек
Введение
Оболочки как элементы разного рода конструкций широко применяются в различных областях техники и строительства. Тонкостенные элементы современных конструкций, представляющие собой оболочки, предназначены для работы под воздействием механических нагрузок, которые могут быть как статическими, так и динамическими.
Известно, что тонкие оболочки могут допускать прогибы, соизмеримые с их толщиной (даже под воздействием нагрузок, далёких от критических значений). Для придания в нужных местах большей жёсткости профиль тонких оболочек может иметь плавные утолщения. С целью повышения жёсткости тонкостенная часть оболочки может быть подкреплена дискретно расположенными рёбрами. В обоих случаях существенно повышается несущая способность конструкции при незначительном увеличении её массы.
Таким образом, всю конструкцию следует рассматривать как конструкцию переменной толщины. В зависимости от характера изменения толщины будем различать оболочки гладко-переменной и, соответственно, ступенчато-переменной толщины (ребристые оболочки).
Расчёты на прочность, устойчивость и колебания оболочечных конструкций играют важную роль при проектировании современных аппаратов, машин и сооружений. Тем не менее, поведение тонкостенных конструкций переменной толщины при котором проявляются геометрическая нелинейность, поперечные сдвиги, нелинейная упругость (пластичность, или, точнее, упругопластичность), ползучесть (вязкоупругость) и ортотропия материалов, а также переменность профиля исследовано недостаточно. Причины этого заключаются в сложности совместного учёта всех упомянутых факторов и необходимости решения громоздких нелинейных краевых задач.
Физические основы теории упругости изложены в энциклопедическом курсе Л.Д. Ландау и Е.М. Лифшица [1]. Прикладные аспекты теории упругости, пластичности и ползучести обстоятельно освещены в трудах Н.И. Безухова [2] и Н.Н. Малинина [3]. Анализ современного состояния теории оболочек, формулировка основополагающих принципов и построение модели упругих оболочек постоянной толщины приводится в весьма содержательной работе П.А. Жилина [4]. Современное состояние теории ребристых оболочек отражено в работах В.В. Карпова [5, 6], а также зарубежных учёных И. Бискова, Дж. Хансена, Б. Чробота, С. Фишера, С. Берта, В. Койтера, Ю. Джиунченга, Р. Лижо, Дж. Маковски, В. Петрашкевича, Х. Стампфа и др. [7-12].
Разработке геометрически нелинейных математических моделей упругих изотропных оболочек (преимущественно пологих) постоянной и переменной толщины посвящены работа В.В. Карпова и др. [5] (в которой использованы результаты исследований В.М. Жгутова [5, с. 8]), а также учебное пособие В.В. Карпова [6]. Однако в этих публикациях практически не рассматриваются оболочки гладко-переменной толщины и не учитывается эффект поперечных сдвигов. Иными словами, рассматриваются модели оболочек постоянной толщины и ребристых, основанные на гипотезе Кирхгофа-Лява.
Надо отметить, что в настоящее время ряд исследователей (Д.О. Астафьев, В.А. Гордон, В.М. Жгутов и др.) склонны считать неверной гипотезу Кирхгофа-Лява (гипотезу прямой нормали) для оболочки, подкреплённой рёбрами (массивными кольцами). Как показано и самим В.В. Карповым [5], поперечные сдвиги могут значительно влиять на напряжённо-деформированное состояние и устойчивость изотропных оболочек при (линейно) упругом их деформировании; с
увеличением же толщины конструкции и жёсткости подкреплений влияние эффекта поперечных сдвигов существенно возрастает. При больших перемещениях влияние поперечных сдвигов резко возрастает.
В ходе вычислительных экспериментов, проведённых В.М. Жгутовым, установлено, что в условиях нелинейной упругости или ползучести материала влияние поперечных сдвигов на напряжённо-деформированной состояние и устойчивость оболочек ещё более усиливается [13].
Геометрическая нелинейность в пособиях [5,6] описывается с помощью квадратичных членов от угловых перемещений в главных нормальных сечениях оболочки, характеризующих одновременно и сдвиги в этих плоскостях, что, очевидно, плохо согласуется с гипотезой Кирхгофа-Лява. По этой причине в задачах динамики при использовании гипотезы Кирхгофа-Лява возникают, как показано в работе В.М. Жгутова [14], некоторые математические некорректности, которые устраняются, исходя из физических соображений. Полученную в этом случае вариационным методом начально-краевую задачу В.М. Жгутов подвергнул корректировке [14]. В работе [14] показано также, что при учёте эффекта поперечных сдвигов математических некорректностей не возникает. Позднее такой же вывод сформулировал и В.В. Карпов с А.Ю. Сальниковым [15].
Кроме того, все результаты в работах [5, 6], относящиеся к ребристым оболочкам, получены с использованием функции Н, описывающей дискретное расположение рёбер по оболочке, их ширину и высоту, которая в общем случае, как показано В.М. Жгутовым [13, 16], оказывается некорректной и непригодной для вычислений.
В работах В.В. Карпова, Д.А. Барановой, Р.Т. Беркалиева и Т.В. Рябиковой [17 - 19] предпринята попытка создания теории ребристых оболочек «в единой системе координат» (на отсчётной поверхности оболочки), отличной от теории, изложенной в работах [5, 6]. Однако данная теория представляется ошибочной, потому что она основана:
1) на использовании «системы координат», которая в действительности таковой не является (введённые в работах [17 - 19] «координаты» не являются взаимно независимыми и оказываются неголономными);
2) на геометрических соотношениях, которые были некорректно выведены с использованием правила дифференцирования сложной функции (известного из курса математического анализа), не подлежащего применению в случае взаимозависимых переменных.
Проектирование и последующее создание легких, но вместе с тем прочных и надежных конструкций, требует дальнейшего совершенствования механических моделей деформируемых тел, а также разработки новых интегральных методов их расчета.
В связи с этим разработка более совершенных математических моделей деформирования оболочек представляется актуальной и важной задачей.
В настоящей статье для задач статики и динамики предложены единые геометрически нелинейные математические модели деформирования оболочек переменной толщины (как гладко-переменной толщины, так и ребристых), совместно учитывающие поперечные сдвиги, возможную ортотропию, нелинейную упругость и ползучесть материалов.
В случае ребристых оболочек на основе уточнённой функции Н [13, 15] учитывается также дискретность расположения рёбер, их ширина, высота, сдвиговая и крутильная жёсткости (в условиях возможного проявления указанных различных свойств материалов).
Постановка задачи
Рассматриваем оболочки общего вида, подразумевая достаточно широкий класс оболочек наиболее важных частных видов: пологие на прямоугольном плане, вращения (сферические, конические, цилиндрические, торообразные и т.д.) и многие другие оболочки.
Некоторую внутреннюю поверхность оболочки принимаем за отсчётную поверхность х3 = 0 . Координатные линии х1 и х2 криволинейной ортогональной системы координат (- а/2 < х1 < а/2 и - Ь/2 < х2 < Ь/2) направляем по линиям кривизны (параллелям и меридианам в случае оболочек вращения), а ось х3- по внутренней нормали отсчётной поверхности так, чтобы система координат х1, х2, х3 была правой. (Полагаем, что опредёленная
таким образом сеть координатных линий на отсчётной поверхности обеспечивает гладкость и регулярность её параметризации).
Элементы длин дуг координатных линий x1, x2 и оси x3 определяем по формулам:
dl1 = H1 • dx1, dl2 = H2 • dx2, dl3 = H 3 • dx3 = dx3,
где H1 = H1(x1,x2), H2 = H2(x1,x2) , H3 = 1 - метрические коэффициенты Лямэ. При этом H1 и H2 зависят от вида оболочки. Например: H1 = H2 = 1 для пологих оболочек и пластин; H1 = const и H2 = H2(x1) в случае оболочек вращения.
Переменную толщину оболочки h = h(xl7x2) задаём ограничивающими её гладкими (или ступенчато-гладкими) поверхностями zB = zB(x1,x2) e Сk и zH = zH(x13 x2) e Сk так, что h = zH -zB и zB < x3 < zH . (Принадлежность функции f (x1, x2) классу гладкости Ck означает, что функция имеет непрерывные частные производные до порядка k > 1 включительно; запись f (x1, x2) е C0 требует только непрерывности по совокупности аргументов). Полагаем, что
векторы (ковекторы) градиентов WzB и VzH отличны от нуля и коллинеарны
rang
Г г dzB / dx1 dzB / dx
v
Bi ^i B 2
dzH / dx1 dzH / dx2
= 1
У
H
в любой точке поверхности x3 = 0 .
Пусть К1 = К1(х1,х2) и К2 = К2(хх,х2) - главные кривизны отсчётной поверхности х3 = 0 оболочки в направлениях хх и х2 соответственно. Для любой точки отсчётной поверхности оболочки (как правило, не являющейся точкой уплощения) хотя бы одно из значений К1 и К2 отлично от нуля: (К1з К2) ^ (0,0). В случае точки уплощения (редком в теории оболочек) К К2) = (0,0) . Для пластин К1 = К2 = 0.
По определению главные радиусы кривизны отсчётной поверхности оболочки равны
Я1 = Я1(х1, х2) = и Я2 = Я1(х1, х2) = . Случаям К1 = 0 V К2 = 0 отвечают «бесконечно К1 К 2 большие» значения Я1 =x¡v Я2 =ю. Здесь V - оператор дизъюнкции предложений (логическое «или»).
Оболочки считаем тонкими, так что для любой точки отсчётной поверхности выполняется условие:
\ ~ ~ 1 1 А = шах^ —;—} < — ,
[ Я 11 20
где Я = ш1п(Я1, Я2) - наименьший из главных радиусов кривизны отсчётной поверхности данной оболочки в рассматриваемой точке; I = тт(^, 12) - наименьшая из длин 11 и 12 координатных
а/2 Ь/2
линий х1 и х2, проведённых в данной точке (11 = |Н1 • dx1 и 12 = |Н2 • dx2).
-а/2 -Ь/2
Как известно, область возможного применения теории тонких оболочек весьма велика. При
рассмотрении тонких оболочек всеми величинами, имеющими порядок малости И /Я (и выше), пренебрегают.
. И И п ,
В случае пластин имеем Д =— (в силу очевидного равенства — = 0), где I - наименьший
I К
из размеров 11 = а и 12 = Ь пластины в плане.
В случае ребристой оболочки за отсчётную поверхность х3 = 0 принимаем срединную поверхность обшивки толщиной И. Рёбра задаем с помощью ступенчато-гладкой функции Н = Н(х1зх2), характеризующей распределение рёбер по оболочке, их конечную ширину и высоту [13, 16]:
М _ N _ N М _ _
Н(х1, х2) = 2И2)5(х - х/) + 2ктб(х2 - х2) -225 - х/Жх2 - х2) ,
] =1 г =1 г=1 ] =1 ( )
где к\г>, 1 < г < N и И^, 1 < j <М - высоты рёбер, подкрепляющих оболочку в направлениях координатных линий х1 и х2 соответственно;
N и М - количества (числа) рёбер в направлениях х1 и х2 соответственно; И = шт(к(1), И2)) - высота фигуры, получающейся при пересечении г - го ребра в направлении х1 и j - го ребра в направлении х2;
5(х1 - х) и 5(х2 - х2) - единичные столбчатые функции, равные по определению единице в местах присоединения рёбер и нулю - вне таких мест.
При этом полагается, что
г (1)
• ширина г-го ребра в направлении х1 равна г(1) = ё -с{ (здесь с = х'2 —'— и
г Г(1)
ё = х'2 +-— , где х2 - ордината осевой линии прикрепления г - го ребра);
г 2 2 2
г (2)
• ширина j-го ребра в направления х2 равна г<(2) = Ь} - а}. (здесь а = х/ —^ и
г(2)
Ь = х/ , где х/ - абцисса осевой линии прикрепления j - го ребра),
поэтому для единичных столбчатых функций можно записать:
5 (х1 - х{ ) =
У , 1 , если aj < x1 < Ьу,
при любом другом x1;
- i fl , если c < x2 < d ,
^(X2 - X2)1 ^ 2 ^1
^ при любом другом X2 .
Таким образом, толщина ребристой оболочки равна h = h + H, причем zB = -h/2 и
zH = h/2 + H.
Отметим, что в формуле (1) высоты рёбер могут быть как постоянными (hf-1 = const и
h(2) = const), так и переменными величинами: hf1 = h(1)(x1,х'2) и h(2) = hj2)(xj,x2) ,
где x2 = const1 (1 < i < N), xj = const1 (1 < j < M); - a/2 < x1 < a/2, - b/2 < x2 < b/2.
Считаем, что рассматриваемые оболочки закреплены по контуру определённым способом и находятся под действием механической нагрузки (статической или динамической).
Учитываем совместно геометрическую нелинейность, поперечные сдвиги, переменную жёсткость, а также возможные проявления ортотропии, нелинейной упругости или вязкоупругости материалов. В случае ребристых оболочек учитываем, кроме того, дискретное расположение рёбер, их ширину, высоту, сдвиговую и крутильную жёсткости (в условиях возможного проявления указанных различных свойств материалов).
Математические модели деформирования рассматриваемых оболочек
Известно, что математическая модель деформирования оболочки состоит из геометрических соотношений (связь деформаций и перемещений), физических соотношений (связь напряжений и деформаций), функционала полной энергии деформации оболочки, из условия стационарности (минимума) которого следуют уравнения (устойчивого) равновесия или движения.
Геометрические соотношения Геометрические соотношения в отсчётной поверхности х3 = 0 оболочки получаются с
помощью операции абсолютного (ковариантного) дифференцирования векторного поля перемещений и с учетом геометрической нелинейности имеют вид
Ри, 1
е11 =- + -
д11 2
РЩ д,11
2
_ Ри2 1 ; е22 = ~д£ + 2
Ри
2
3 Ри2 Ри1 Ри3 Ри3
V д12 J
; Ум =у21 =—+—+■ • ,
12 21 д11 д\2 д11 д12 (2)
где е11 и е22 - деформации удлинения (сжатия) вдоль координатных линий х1 и х2 соответственно;
у12 =у21 - деформации сдвига в касательной плоскости (dx1,dx2) ;
и1 = и1(х1,х2), и2 = и2(х1,х2) и и3 = и3(х1, х2) - компоненты вектора перемещений (точек отсчётной поверхности) вдоль координатных линий х1 , х2 и оси х3 соответственно;
Р
-, 1 <а< 3 - операторы ковариантного дифференцирования (по направлениям 1а)
д1а
произвольных полей, в частности: скалярного поля а = а(х1,х2,х3), векторного поля а; = а{(х1, х2, х3), 1 < / < 3, поля тензора второго ранга ай = ай(х1, х2, х3), 1 < ¡,к < 3 и т.д.
Р
Операторы ковариантного дифференцирования-действуют по правилам [20, 21]:
д1а
Ра 1 да
а а
д1а Н а дха
Ра, 1 да; Д ^
а а-=----X, ак Тгка
д1а Н а дха к=1
и, соответственно,
Рагк 1 дагк ^
а а
гк д1а На дх
X (а1к Гг1а + аг1 Гк1а ) ,
1а а ила I=1
где Ггка = Т^^,х2,х3) - символы Кристоффеля (1-го рода), 1 <г,к,а<3 .
Как известно, символы Кристоффеля симметричны по крайним индексам при к Ф г, к Фа (Г,ка = Га) и антисимметричны по первым двум индексам (Ггка = -Гка), а потому величины Ггка с разными значениями индексов равны нулю (Ггка = 0 при ¡Ф к ,г Фа, к Фа). Это значит, что
в ортогональной криволинейной системе координат из 27 величин Ггка ненулевыми могут быть не более 12: г1кк = -ГЫк .
^ 1 дНк 1 д 1п Нк При этом Цкк =-Гкк Л = 77---^^ , а те 12 из 27 величин Ггка , к°т°рЫе в
НгН к дхг Нг дхг
ортогональной криволинейной системе координат могут быть отличными от нуля, имеют вид
Г =-Г = 1 дН 2 ; Г =-Г = 0 ; Г =-Г = 1 дН1 ; Г =-Г = 0 ;
-■-122 212 ^ ^ - ' Х133 х 313 и>А211 1121 ^ ^ - ' 1 233 1 323
Н,Н 2 дх1 Н,Н 2 дх2
Г =-Г = — = -Г ; Г =-Г = — дН2 =-К
311 131 гг Л1 ; 322 Х232 „ ~ Л2 .
Л1 С/Л3 Л 2 С/Л 2
Таким образом,
Бп, 1 дм, 1 дН, Бп, 1 дм, 1 дН2
—1 =---1 +---м2 - К,-п3; —1 =---1 +---2-п2;
д/, Н, дх, Н,Н2 дх2 д/2 Н2 дх2 Н,Н2 дх,
Бп2 1 дп2 1 дН, Бп2 1 дп2 1 дН2
-2 =---2 +---п,; —2 =---2 +---2-м -К2 -п3;
д/, Н, дх, Н,Н2 дх2 д/2 Н2 дх2 Н,Н2 дх, (3)
Бп3 1 дп3 Бп3 1 дм3
—3 =---3 + К, -п,; —3 =---3 + К2 -п2.
д/, Н, дх, д/2 Н 2 дх2
Введём обозначения:
©, = ©1(х1, х2) = ^ = ^ + К,-м1 ; ©2 =©2(х1. х2) = = ^ + К2 -п2 . (4) д/, Н, дх, д/2 Н2 дх2 (4)
Тогда с учётом выражений (3) и (4) геометрические соотношения (2) принимают вид:
1 дм, 1 дН, ^ 1 ©2
, =---1 +----1 • п2 - К, • п3 + — ©,;
11 н, ах, Н,Н2 ах2 2 1 3 2
1 дп2 1 дН2 1 2
е22 =---2 +---2-п, -К2 -п3 + — ©-,;
22 Н2 дх2 Н1Н2 дх1 1 2 3 2 2 (5)
1 дп2 1 дН, 1 дп, 1 дН2 ^ ^ Х12 = х21 =---2 +----1 • п, +----к +----2• п2 + ©, • ©2.
12 21 Н1 дх1 Н1 Н2 дх2 1 Н2 дх2 Н1 Н2 дх1 2 1 2
В ряде случаев можно полагать, что в процессе деформирования
1 дп3 1 дп3
---3 >> К, • п,;---3 >> К2 • п2,
Н1 дх1 1 1 Н2 дх2 2 2
а значит,
1 дп3 1 дп3
©, ~---; ©2 ~---.
Н1 дх1 Н2 дх2
Для пластин К, = К2 = 0 и Н, = Н2 = 1 (как было отмечено выше), а потому операторы
Б
ковариантного дифференцирования -, 1 <а< 3 совпадают с операторами обычного
д/а
д
дифференцирования-и геометрические соотношения (2), (5) максимально упрощаются:
дх„
=
ди1 1 дх1 2
/ \2 ди3
дх1
^22 =
ди2 дх2
1
+ — 2
2
ди3
V дХ2 У
ди2 ди1 ди3 ди
У12 = У21 = + ^ + ^
дх1 дх2 дх1 дх2
(6)
В соотношениях (2), (5) и (6) квадратичные члены характеризуют (приближённо) геометрическую нелинейность, которую следует учитывать в случаях, когда поперечные
перемещения и3 (прогибы) соизмеримы с толщиной оболочки И .
Подчеркнём, что тензор деформации при этом должен быть мал [1]; практически в нашем случае это означает требование и3 <<I (для любой точки отсчётной поверхности оболочки), т.е. прогибы должны быть малы по сравнению с продольными размерами оболочки.
Примечание. Следует помнить, что использование приближённых квадратичных выражений для деформаций может привести к ошибкам в практических расчётах. Но ошибки эти связаны, как правило, не с приближённостью исходных формул («приближенность» - всегда правило, а не исключение), а с тем, что при составлении уравнений равновесия зачастую «забывают» учитывать члены, обусловленные геометрической нелинейностью (которая квадратичными членами и описывается). В нашем случае, как указано ниже, уравнения равновесия (или движения) могут быть получены из функционала полной энергии деформации оболочки вариационным методом с помощью корректной математической процедуры вывода, что не приведёт к неучёту соответствующих важных членов. Для задач статики уравнения равновесия могут не составляться вообще, поскольку к функционалу полной энергии деформации оболочки могут быть применены «прямые» методы (метод Ритца). В этом случае в полученной системе алгебраических уравнений члены, описывающие геометрическую нелинейность, будут присутствовать.
Разумеется, было бы гораздо лучше осуществлять учёт геометрической нелинейности путём учёта изменения формы деформируемого тела на каждом шаге его нагружения. Но в случае оболочек это привело бы к чрезмерным математическим усложнениям.
Подход к решению геометрически нелинейных задач в теории оболочек, основанный на использовании квадратичных выражений для деформаций применён в работах В.В. Карпова и др. [5, 6]. Такой же подход приведён Л.Д. Ландау и Е.М. Лифшицем при изложении теории сильного изгиба тонких пластин в курсе [1].
Деформации поперечных сдвигов определяем по формулам:
У13 = с/(Х3 )ф1 ; У23 = с/(Х3 )Ф2 .
Здесь /(х3) - функция, характеризующая распределение напряжений т13 и т23 в главных
1 2н
нормальных сечениях оболочки, такая, что /(гВ ) = /(гн ) = 0, — | / (х3)ёх3 = 1,
И гВ
1 гн
— [/2(х3)ёх3 = Ис (с - константа); *
= 1/с (с
и.......
гВ
, , 3 . , Ои3 Ф1 = +--3 и Ф2 = *т2 +--3 - полные углы сдвигов, где ^ = и = tg^2, причём
д1х д12
= ¥\(Х\,х2) и у/2 = у2(хх,х2) - углы поворота отрезка нормали к отсчётной поверхности в соответствующих главных нормальных сечениях оболочки [19].
В качестве /(х3) используем квадратичную зависимость [19] (применимую как для оболочек гладко-переменной, так и ступенчато-переменной толщины):
/(Хь) = -14(х3 - гнХх3 - гв) = /о + /х3 + /2Х
И
6 гнгВ . 6(гн +гВ) ^ 6
где /о =--; / = ' н~2 В' и /2 =-^, и тогда с = 5/6.
И2 И И
Деформации в слоях x3 = const вычисляем по формулам [19]:
S11 ) = S11 + X3 'X11 ' ^22 = ^22 + X3 'X 22 ' Y12 = Y12 + X3 " 2X12 '
РФ, 1 дФ, 1 dH, ф где X11 =-L =--1 +--LФ2 ;
5/1 H1 dx1 H1H 2 dx2
r _ РФ 2 _ 1 дФ 2 1 dH2
L22 _ — +
dl2 H2 dx2 H1H2 dx1
РФ2 РФ, 1 5Ф2 1 5Ф, 1 2X,, _-2 +-1 —--2 + 1
12 dl1 dl2 H1 dx1 H2 öx2 HlH2
fdH, Ф dH2 Ф ^ —1 Ф1 +—2 Ф2 V öx2 dx1
поскольку для перемещений в слоях x3 = const, полагаем:
Mx3) _ Mj + x3 • Ф1, u(2x3) _ u2 + x3 • Ф2 и u3x3) _ u3.
1
Здесь xlk _ 2
fРФk РФ ^ -—+-:
V dli dlk J
, 1 < /, к < 2 - тензор изменения кривизны и кручения (записанный в
явно симметричном виде).
Физические соотношения
Физические соотношения для ортотропных линейно упругих материалов в соответствии с обобщенным законом Гука имеют вид [22 - 24]:
_ Gu (s^ + ^f) , ^22 _ G22 (s2f + ^if ), Т _ G^if ; Т _ G13Y13 , Т _ G23Y23 . (7)
Здесь G11 = Е1 /(1 - //) , G22 = Е2 /(1 - //); G12, G13, G23 - модули сдвига, где Е1, Е2 и /, /2- продольные модули Юнга и коэффициенты Пуассона (/Е2 = /2Ех) .
Для изотропных линейно упругих материалов физические соотношения следуют из формул (7), если положить Е1 = Е2 = Е, / = /2 = / и, следовательно, G11 = G22 = Е/(1 — /2) ,
= G13 = ^3 = G = Е/2(1 + /).
В случае нелинейно упругого изотропного материала физические соотношения в соответствии с деформационной теорией пластичности принимаем в виде [23, 24]:
ЕР ЕР ЕР ЕР ЕР
= — , ^22 = ^22 — ®22 , Т12 = Т12 — Т12 ; Т13 = Т13 — Т13 , Т23 = Т23 — Т23 , (8)
где линейно упругие составляющие тензора напряжений (отмечены индексом « Е ») определяются по формулам (7) (где Е1 = Е2 = Е и / = /2 = /), а нелинейно-упругие составляющие (отмечены индексом Р) вычисляются с помощью соотношений:
Р / \ Е Р / \ Е Р / \ Е Р / \ЕР /\Е
аи = )ап , ^22 = Ю(8г )а22 ; = ^ , Т13 = ®(ег )Т13 , Т23 = )Т23 . Здесь гэ(ег ) - функция А.А. Ильюшина; е - интенсивность деформации.
При учёте вязкоупругости (ползучести) материала физические соотношения в соответствии с теорией упруго-ползучего тела (применимой для полимерных материалов, стареющего и нестареющего бетона) имеют вид:
_ = _Е — _С _ = _Е — _С Т = _Е — _С ; = Е — С т = ТЕ — _С (9)
°11 _ °11 °11 , °22 _ 22 22 , П2 _ П2 П2 ; Т13 = Т13 — Т13 , Г23 _ Г23 Г23 , (9)
где линейно упругие (упруго-мгновенные) составляющие тензора напряжений (отмечены индексом « Е ») имеют вид соотношений (7), в которых все константы в общем случае являются функциями времени наблюдения t, а вязкоупругие составляющие (отмечены индексом « С ») вычисляются по формулам:
О = Я\^, ^ст* , а22 = т, ФЕ; тС = , Т = т, *К3, Тз = Я2(^ s)тE.
Здесь , 5) и 5) - функции влияния при растяжении (сжатии) и сдвиге соответственно, где 5 - время, предшествующее моменту наблюдения I.
Функционал полной энергии Функционал полной энергии деформации оболочки на данном отрезке времени [0,^] для задач динамики имеет вид [25, 26]:
W = - U + А)dt,
(10)
где К и и - кинетическая и потенциальная энергии оболочки; А - работа внешних сил.
Подчеркнём, что в задачах динамики подлежащие определению функции перемещений и1, и2, и3 и углов %, %2 являются не только функциями координат х, у (считающихся неподвижными), но и времени I: и1 = и1(х,у, t), и2 = и2(х,у, t), и3 = и3(х,у,I) и % х,у,t),
% =%(х,у,t) .
Для задач статики полная энергия деформации оболочки (функционал Лагранжа) может быть записана в виде [25, 26]:
W = U - А .
В соотношениях (10) и (11):
*=Р !!!
Q
HЛ
ди
(x3) Л2 (
dt
+
H 2
ди
(х3) Л2 (ди(x3) Л2
dt
+
dt
V У
dQ;
U = \ !!!(in3) +a22423) +T\2Yn) +т\зУиЪ) +т23у23) )dQ;
А = !!(Pxux + P2u2 + P3u3)dS.
(11)
(12)
(13)
(14)
Здесь р - плотность материала оболочки (р « const) ;
Q = [-a/2,a/2]х[-b/2,b/2]x[zB,zH] - компакт (замкнутое связное множество) в пространстве
(х, y,z);
S = [-a/2,a/ 2]х [- b/2,b/2] - компакт на плоскости (x,y);
dQ = dlldl2dx3 = HlH2dxldx2dx3 и dS = dl1dl2 = HlH2dxidx2 - элементы объёма и отсчётной поверхности оболочки соответственно;
Pl, P2 и P3 - компоненты внешней механической нагрузки в направлениях xl, x2 и x3 соответственно, действующей на единицу площади поверхности оболочки (в задачах статики
Px = P1(x,y) , P2 = P2(x,y) и P3 = P3(x,y); в задачах динамики Px = P1(x,y,t) , P2 = P2(x,y,t) и
P3 = P3( x, y, t) ).
Проинтегрируем по переменной x3 выражения (12) и (13) в функционалах (10) и (11).
0
S
В результате будем иметь соответственно:
* ii( 4
H ^ J2 +Г h 2 ^ 12 4ди3 12
dt
dt
+
+ 2 4
H 2 dUL + H 2 ди2 дФ 2
4 dt dt
dt dt
+ 4
H ^Y +fH дФ2^
4 dt
dt
ds
(16)
^ _ 1 Я[П + N22S22 + N12X12 + MnXn + 2M12X12 + M22X22 + + 023Ф2 ] ■
(17)
где
40 _ {dx3 _ h , 4 _ { x3dx3 _ h •
Zh ~ z + Z ^H ~
_ { x3dx3 _ h •—^—В и 42 _ {x32dx3 _ h •
2
2 + + 2 ^ + Z Я ZB + Z В
3
- погонные
^ J ^3
площадь, статический момент и момент инерции данного главного нормального сечения оболочки;
Zh
Zh
N11 = Jo"ndx3 и N22 = Ja"ndx3 - погонные усилия растяжения (сжатия) в направлениях x1 и x2;
Zb z в
Zh
N12 = jz12dx3 - погонное сдвиговое усилие в касательной плоскости (dx1, dx2) ;
Zb
Zh Zh
Mu = Jo^^d^ и M22 = Ja22x3dx3 - погонные изгибающие моменты в направлениях xx и x2;
Zb Zb
Zh
M12 = Jt12x3dx3 - погонный крутящий момент в касательной плоскости (dxl,dx2) ;
Zb
Zh Zh
Q13 = jr13dx3 и Q23 = jT23dx3 - погонные поперечные (перерезывающие) усилия в главных
Zb Zb
нормальных сечениях оболочки соответственно.
Заметим, что в случае симметричной по толщине оболочки (zB = -zH) величины
zH
Am-1 = Jzm-1dz = 0 , если m -1 нечётно (сравнимо с 1 по модулю 2: m -1 = 1(mod 2)), 1 < m < 3
zb
и соотношения для внутренних силовых факторов существенно упрощаются. Поэтому целесообразно (если при этом не нарушается условие гладкости отсчётной поверхности, как,
например, в случае ребристых оболочек) ввести новую отсчётную поверхность z0 = z0(x, y) так, чтобы Am-1 = 0, если m -1 = 1(mod2). Положение поверхности z0 можно определить из условия
Z" \т-\
z0 =(zH + zB)/2 [22] и тогда Am-1 = J (z*j dz*= 0 , где m -1 = 1(mod2), z* = z - z0,
* * ZЯ _ ZЯ — Z0 ' ZB _ ZВ — Z0 .
Компоненты тензора напряжений стп, сг22, т12 , т13 и т23 в выражениях для внутренних силовых факторов, действующих в оболочке (усилий , Ы22, Ыи, Q13, Q23 и моментов Мп,
М22,М12) вычисляем по формулам (7), (8) или (9) в зависимости от проявленных свойств материала (ортотропия, упругость, нелинейная упругость, вязкоупругость).
s
и
Z
Уравнения движения или равновесия оболочки могут быть получены, исходя из фундаментальных принципов наименьшего действия (в форме Гамильтона-Остроградского) или минимума потенциальной энергии (в форме Лагранжа) [25, 26], в соответствии с которыми:
т = 5|( - и+а) = о (18)
и
или, соответственно,
т = 5(и - А) = о,
где 5 - символ вариации.
Таким образом, построены единые математические модели деформирования оболочек переменной толщины с учётом различных свойств материалов (для задач статики и динамики).
Применение этих моделей позволяет совместно учитывать геометрическую нелинейность, поперечные сдвиги, возможную ортотропию, нелинейную упругость и ползучесть материалов. В случае ребристых оболочек учитывается также дискретность расположения рёбер, их ширина, высота, сдвиговая и крутильная жёсткости (в условиях возможного проявления указанных различных свойств материалов).
Начальные и граничные условия (отвечающие способу закрепления оболочки по контуру) предполагаются заданными.
В задачах статики для поиска минимума функционала (11), записанного с учётом выражений (15)-(17), может быть эффективно применён метод Ритца при разложении искомых функций перемещений и1, и2, и3 и углов % и Ч2 в ряды.
В задачах динамики для решения системы уравнений движения оболочки, эквивалентной вариационному уравнению (18), целесообразно последовательно применить методы Власова-Канторовича и Рунге-Кутта.
Если априори считать функции перемещений и1, и2, и3 компонентами потенциального
вектора (что, очевидно, почти всегда выполняется, поскольку возможные локальные повороты, нарушающие симметрию деформаций сдвига, не искажают картину деформированного и напряжённого состояния оболочки), то в этом случае мы имеем, что
Пиг Пи2 Ч Пиг Пи3 Ч Пи2 Пи3 Ч1 = "1 = , Ч 2 = '
д12 д11 дх3 д11 дх3 д12
а, значит, полные углы сдвигов:
Ф, = 2—3, Ф = 2- 3
и, кроме того,
д11 2 д12
х = 2 Ви3 х = 2 Ви3 = 2 Р\
Л,, — 2- , Л.— 2- , 2-Л-п — 2
11 д\х 22 д12 12 д11 д12 Таким образом, полученные математические модели существенно упрощаются.
Литература
1. Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. Т. VII. Теория упругости. М. : Физматлит, 2007. 264 с.
2. Безухов Н. И. Основы теории упругости, пластичности и ползучести. М.: Высшая школа, 1968. 512 с.
3. Малинин Н. Н. Прикладная теория пластичности и ползучести. М.: Машиностроение, 1975. 399 с.
4. Жилин П. А. Прикладная механика. Основы теории оболочек. СПб.: Изд-во Политехн. ун-та, 2006. 167 с.
5. Карпов В. В., Игнатьев О. В., Сальников А. Ю. Нелинейные математические модели деформирования оболочек переменной толщины и алгоритмы их исследования. М.: АСВ; СПб.: СПбГАСУ, 2002. 420 с.
6. Карпов В. В. Математическое моделирование, алгоритмы исследования модели, вычислительный эксперимент в теории оболочек. СПб.: СПбГАСУ, 2006. 330 с.
7. Byskov E., Hansen J. C. Postbuckling and imperfection sensitivity analysis of axially stiffened cylindrical shells with mode interaction // Journal of structural mechanics. 1980. № 2. Pp. 205-224.
8. Chrobot B. Mathematical models of ribbed shells // Studia Geotechnica et Mechanica. 1982. Vol. IV, №3-4. Pp. 55-68.
9. Fisher C. A., Berit C. W. Dynamic buckling of an axially compressed cylindrical shells with discrete rings and stringers // Trans. ASME. Ser. E. 1973. Vol. 40, № 3. Pp. 736-740.
10. Koiter W. T. General theory of mode interaction in stiffened plate and shell structures // WTHD Rep. 1976. № 590.
11. Wu Jiuncheng, Pan Lizhou. Nonlinear theory of multilayer sandwich shells and its application // Applied Math and Mechanics. 1997. Vol. 18, № 1. Pp. 19 -27.
12. Makowski J., Pietraszkiewicz W., Stumpf H. On the general form of jump conditions for thin irregular shells // Archives of Mechanics. 1998. Vol. 50, № 3. Pp. 483-495.
13. Жгутов В. М. Ответ профессору Карпову, Владимиру Васильевичу (о научном приоритете в методе конструктивной анизотропии для ребристых оболочек и на функционал, описывающий ползучесть их материала) // Инженерно-строительный журнал. 2011. № 3. С. 75 - 80.
14. Жгутов В. М. Нелинейные свободные колебания пологих оболочек ступенчато-переменной толщины: дисс. канд. техн. наук. СПб. , 2004. 177 с.
15. Карпов В. В., Сальников А. Ю. Вариационный метод вывода нелинейных уравнений движения пологих ребристых оболочек // Вестник гражданских инженеров. 2008. № 4. С. 121-124.
16. Жгутов В. М. Метод конструктивной анизотропии для ортотропных и изотропных ребристых оболочек // Инженерно-строительный журнал. 2009. № 8. С. 40 - 46.
17. Карпов В. В., Баранова Д. А., Беркалиев Р. Т. Программный комплекс исследования устойчивости оболочек. СПб. : СПбГАСУ, 2009. 102 с.
18. Карпов В. В. Уравнения равновесия для оболочки вращения для единой системы координат // Вестник гражданских инженеров. 2010. С. 173-179.
19. Карпов В. В., Рябикова Т. В. Оболочки вращения в единой системе координат // Вестник гражданских инженеров. 2010. № 2. С. 44-47.
20. Кочин Н. Е. Векторное исчисление и начала тензорного исчисления. М. : Наука, 1965. 428 с.
21. Акивис М. А., Гольдберг В. В. Тензорное исчисление. М. : Наука, Физматлит, 1969. 352 с.
22. Жгутов В. М. Геометрически нелинейные математические модели термоупругости оболочек переменной толщины // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. 2011. № 4. С. 46 - 56.
23. Жгутов В. М. Исследование прочности и устойчивости ребристых оболочек с помощью вычислительного эксперимента // Проблемы прочности материалов и сооружений на транспорте: Сборник докладов VII Международной конференции по проблемам прочности материалов и сооружений на транспорте 23-24 апреля 2008 года. СПб. : Петербургский государственный университет путей сообщения, 2008. С. 110 - 131.
24. Жгутов В. М. Математические модели и алгоритмы исследования устойчивости пологих ребристых оболочек при учёте различных свойств материала // Изв. Орловского гос. техн. ун-та. 2007. № 4. С. 20 - 23.
25. Компанеец А. С. Теоретическая физика. М. : Гос. изд-во техн.-теор. лит-ры, 1957. 564 с.
26. Жилин П. А. Теоретическая механика. Фундаментальные законы механики. СПб. : Изд-во СПбГТУ. 339 с.
*Владимир Михайлович Жгутов, Санкт-Петербург, Россия Тел. раб.: +7(812)378-20-83; эл. почта: abc_kitezh@mail.ru
Mathematical deformation models of variable thickness shells with calculation of different materials' behaviour
V.M. Zhgoutov
"Architectural and Engineering Company Kitezh" LTd.
+7(812)378-20-83; e-mail:abc_kitezh@mail.ru
Key words
variable thickness shells; linear and nonlinear elasticity; viscoelasticity, transversal movements; full shear angles; variable rigidity; shear and torsion rib rigidity; full deformation energy functional of shells
Abstract
Shells as elements of various structures are widely used in different fields of engineering and construction. Thin-walled elements of modern structures that represent shells designed to work under the influence of mechanical stress, which can be either static or dynamic. Calculations of strength, stability and vibrations of shell structures have an important role in the design of modern devices, machines and constructions.
Profile of shell can have smooth thickening to increase rigidity in some areas. There is also a variant of a thin-walled part of the shell strengthening by discretely spaced ribs.In both cases, carrying force of structures significantly increases with a slight increase of its mass.
In this article geometrically nonlinear mathematical deformation models of shells with variable thickness, in particular ribbed shells, (for static and dynamic problems) are proposed. Different material properties are taken into account (orthotropism, linear and nonlinear elasticity, viscoelasticity and creeping), as well as transversal movements, variable rigidity of ribbed shells besides its finite widths and height of ribs, shear and torsion rigidity.
References
1. Landau L. D., Lifshits Ye. M. Teoreticheskaya fizika. T. VII. Teoriya uprugosti [Theoretical Physics. Vol. VII. Theory of elasticity]. Moscow: FIZMATLIT, 2007. 264 p. (rus)
2. Bezukhov N. I. Osnovy teorii uprugosti, plastichnosti i polzuchesti [Fundamentals of the theory of elasticity, plasticity and creep]. Moscow: Vysshaya shkola, 1968. 512 p. (rus)
3. Malinin N. N. Prikladnaya teoriya plastichnosti i polzuchesti [Applied theory of plasticity and creep]. Moscow: Mashinostroyeniye, 1975. 399 p. (rus)
4. Zhilin P. A. Prikladnaya mekhanika. Osnovy teorii obolochek [Applied Mechanics. Fundamentals of the theory of shells]. Saint-Petersburg: Izd-vo Politekhn. un-ta, 2006. 167 p. (rus)
5. Karpov V. V., Ignatyev O. V., Salnikov A. Yu. Nelineynyye matematicheskiye modeli deformirovaniya obolochek peremennoy tolshchiny i algoritmy ikh issledovaniya [Non-linear mathematical models of variable thickness shells deformation and algorithms for their studies]. Moscow: ASV; Saint-Petersburg: SPbGASU, 2002. 420 p. (rus)
6. Karpov V. V. Matematicheskoye modelirovaniye, algoritmy issledovaniya modeli, vychislitelnyy eksperiment v teorii obolochek [Mathematical modeling, algorithms, research models, computer experiment in the theory of shells]. Saint-Petersburg: SPbGASU, 2006. 330 p. (rus)
7. Byskov E., Hansen J. C. Postbuckling and imperfection sensitivity analysis of axially stiffened cylindrical shells with mode interaction. Journal of structural mechanics. 1980. No. 2. Pp. 205-224. (rus)
8. Chrobot B. Mathematical models of ribbed shells. Studia Geotechnica et Mechanica. 1982. Vol. IV, No. 34. Pp. 55-68.
9. Fisher C. A., Berit C. W. Dynamic buckling of an axially compressed cylindrical shells with discrete rings and stringers. Trans. ASME. Ser. E. 1973. Vol. 40, No. 3. Pp. 736-740.
10. Koiter W. T. General theory of mode interaction in stiffened plate and shell structures. WTHD Rep. 1976. No. 590.
11. Wu Jiuncheng, Pan Lizhou. Nonlinear theory of multilayer sandwich shells and its application. Applied Math and Mechanics. 1997. Vol. 18, No. 1. Pp. 19 -27.
12. Makowski J., Pietraszkiewicz W., Stumpf H. On the general form of jump conditions for thin irregular shells. Archives of Mechanics. 1998. Vol. 50, No. 3. Pp. 483-495.
13. Zhgutov V. M. Magazine of civil engineering. 2011. No. 3. Pp. 75 - 80. (rus)
14. Zhgutov V. M. Nelineynyye svobodnyye kolebaniya pologikh obolochek stupenchato-peremennoy tolshchiny: diss. kand. tekhn. nauk [Nonlinear free vibrations of shallow stepwise-variable thickness shells: diss. Candidate. Technical. Science]. Saint-Petersburg: 2004. 177 p. (rus)
15. Zhgutov V. M. Magazine of civil engineering. 2009. No. 8. Pp. 40 - 46. (rus)
16. Karpov V. V., Baranova D. A., Berkaliyev R. T. Programmnyy kompleks issledovaniya ustoychivosti obolochek [The program complex for studying stability of shells]. Saint-Petersburg: SPbGASU, 2009. 102 p. (rus)
17. Kochin N. Ye. Vektornoye ischisleniye i nachala tenzornogo ischisleniya [Vector Calculus and the beginning of the tensor calculus]. Moscow: Nauka, 1965. 428 p. (rus)
18. Akivis M. A., Goldberg V. V. Tenzornoye ischisleniye [Tensor calculus]. Moscow: Nauka, Fizmatlit, 1969. 352 p. (rus)
19. Zhgutov V. M. Nauchno-tekhnicheskiye vedomosti Sankt-Peterburgskogo gosudarstvennogo politekhnicheskogo universiteta [Scientific and technical statements of the St. Petersburg State Polytechnical University]. 2011. No. 4. Pp. 46 - 56. (rus)
20. Zhgutov V. M. Problemy prochnosti materialov i sooruzheniy na transporte: Sbornik dokladov VII Mezhdunarodnoy konferentsii po problemam prochnosti materialov i sooruzheniy na transporte 23-24 aprelya 2008 goda [Problems of materials and structures strength for transport: Proceedings of the VII International Conference on the strength of materials and structures for transport on 23-24 April 2008]. Saint-Petersburg: Peterburgskiy gosudarstvennyy universitet putey soobshcheniya, 2008. Pp. 110 - 131. (rus)
21. Zhgutov V. M. Izvestiya Orlovskogo gosudarstvennogo tekhnicheskogo universiteta [Proceedings of the Oryol State Technical University]. 2007. No. 4. Pp. 20 - 23. (rus)
22. Kompaneyets A. S. Teoreticheskaya fizika [Theoretical physics]. Moscow: Gos. izd-vo tekhn. teor. lit-ry, 1957. 564 p. (rus)
23. Zhilin P. A. Teoreticheskaya mekhanika. Fundamentalnyye zakony mekhaniki [Theoretical Mechanics. The fundamental laws of mechanics]. Saint-Petersburg: Izd-vo SPbGTU. 339 p. (rus)
Full text of this article in Russian: pp. 79-90
