CC BY
9
№ 3.1 (72.1)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2020 г.
МОДЕЛЬ ДЕФОРМИРОВАНИЯ УПРУГОГО СЛОЯ ПЕРЕМЕННОЙ ТОЛЩИНЫ
Жолдасова Шарайна Амановна
канд. тех. наук, профессорАУНГ, Атыраускийуниверситет нефти и газа им. С. Утебаева
Республика Казахстан, г. Атырау E-mail: zh. sharaina@mail. ru
Коданова Шынар Кулмагамбетовна
канд. тех. наук, профессор АУНГ, Атырауский университет нефти и газа им. С. Утебаева
Республика Казахстан, г. Атырау E-mail: kodanova_s@mail. ru
Шабдиров Дарын Насипкалиевич
канд. физ.-мат. наук, профессор АУНГ, Атырауский университет нефти и газа им. С. Утебаева
Республика Казахстан, г. Атырау E-mail: [email protected]
MODEL OF DEFORMATION OF AN ELASTIC LAYER OF VARIABLE THICKNESS
Sharaina Zholdasova
Cand. technical Sciences, Professor AUNG, Atyrau University of Oil and Gas named after S. Utebayev,
Republic of Kazakhstan, Atyrau
Shynar Kodanova
Cand. technical Sciences, Professor AUNG, Atyrau University of Oil and Gas named after S. Utebayev,
Republic of Kazakhstan, Atyrau
Daryn Shabdirov
Cand. Ph. D., Professor AUNG, Atyrau University of Oil and Gas named after S. Utebayev,
Republic of Kazakhstan, Atyrau
АННОТАЦИЯ
Модель деформирования упругого слоя (пластин, оболочек) переменной толщины при малых удлинениях, сдвигах и углах поворота в общем случае представляют собой трехмерные уравнения линейной теории упругости. В работе [1] изложен метод понижения размерности задач упругого деформирования пластин и оболочек постоянной толщины с произвольными условиями для перемещений и напряжений на лицевых поверхностях. В работе [2] эта методика использовалась для построения уравнений плоского изгиба неоднородного упругого слоя переменной толщины. Ниже на основе этой же методики получены уравнения деформирования упругого слоя переменной толщины в произвольной криволинейной системе координат с произвольными условиями для напряжений и перемещений на лицевых поверхностях.
ABSTRACT
The model of deformation of an elastic layer (plates, shells) of variable thickness with small elongations, shifts, and rotation angles is generally a three-dimensional equation of the linear theory of elasticity. In [1], a method for reducing the dimension of problems of elastic deformation of plates and shells of constant thickness with arbitrary conditions for displacements and stresses on the front surfaces is described. In [2], this method was used to construct the equations of flat bending of a non-uniform elastic layer of variable thickness. The equations of deformation of an elastic layer of variable thickness in an arbitrary curvilinear coordinate system with arbitrary conditions for stresses and displacements on the front surfaces are obtained using the same technique below.
Ключевые слова: упругий слой, упругое деформирование, переменная толщина, теория упругости, ковари-антный базис, тензор напряжений, закона Гука, полиномы Лежандра, тангенциальное и поперечное перемещение.
Keywords: elastic layer, elastic deformation, variable thickness, elasticity theory, covariant basis, stress tensor, Hooke's law, Legendre polynomials, tangential and transverse displacement.
Библиографическое описание: Жолдасова Ш.А., Коданова Ш.К., Шабдиров Д.Н. Модель деформирования упругого слоя переменной толщины // Universum: Технические науки : электрон. научн. журн. 2020. N° 3.1(72.1). URL:
http://7universum. com/ru/tech/archive/item/9106
№ 3.1 (72.1)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2020 г.
1. Уравнения линейной теории упругости в векторной форме в произвольной системе координат ф состоят из уравнений равновесия.
d^ + f = 0
и обобщенного закона Гука [3]
(1.1)
Дифференцируя (2.1), получим выражения для локального базиса э системы координат %1 :
Эа= + 2 djö ^ ,Эз= 2 =ПП (Л3)
В каждой точке срединной поверхности So тройку векторов^} :
il = fgCv g , Сч = aiksirk х Э;, (1.2)
где
4а = эг
OS х rs), il = Jgothj.f = Jüf (1.3)
и приняты обозначения: э^ - ковариантный базис системы координат ^ а1'- компоненты тензора напряжений; и- вектор перемещений; f- вектор массовых сил; alksi - заданные функции координат, axb -прямое (тензорное) произведение. В (1.1)-(1.3) и в дальнейшем используется правило суммирования по повторяющимся индексам: латинские индексы пробегают значения 1,2,3, а греческие- 1,2.
Для простоты изложения ограничимся случаем следующих краевых условий на границе
Vi + q •и = G,
где VI -компоненты единичного вектора внешней нормали к 8^- кусочно-постоянная функция, принимающая значения либо 0, либо 1; С - заданная вектор-функция на 8. Таким образом, условия (1.4) обеспечивают задание на границе либо перемещений ^=0), либо напряжений ^=1).
2. Слой определим, как некоторую область V, для любой внутренней точки которой радиус-вектор / имеет вид
т = (2.1)
где
? е Ъ = ? е [-ЪИ (2.2)
_ х+ + х- —> _(
хп =-,Ах = х+ — х =
о 2
и приняты обозначения: И-половина толщины слоя в направлении^3 ( 2к = Ах),П -единичный вектор в направлении
Функции отображают плоскую об-
ласть с контуром на лицевые поверхности слоя Б±, а функция х0,(^отображает на срединную поверхность Б0 .
Пусть Е -боковая поверхность слоя, Ь - линия пересечения Е и 80. Согласно (2.1), Е-линейчатая поверхность, образованная семейством прямых линий, проходящих через точки контура Ь в направлении вектора П. Вектор Л может не быть нормалью к 80, соответственно, поверхности Еи 80 могут пересекаться не под прямым углом.
а^ = Л х (э^ х п) 'а3 = П'З^ = э3^3 = 0, (2.4)
будем рассматривать как новый локальный базис. Соответствующий биортогональный базис {а1 } имеет вид
а3 = э3 ,а3 = П'Э3 = эк|^3 = 0, (2.5)
где э1 -контравариантный базис системы координат . Из (2.3),(2.4) получаем связь э'1са'1:
(2.6)
Замечание. Поскольку в традиционном определении оболочки [4] вектор является нормалью к срединной поверхности 80 , а, следовательно, Е и 80 -ортогональны, представляется оправданным введение понятия слоя (соотношения (2.1)), которое следует рассматривать как обобщение на случай оболочки переменной толщины с наклонными боковыми поверхностями.
3. Разложения по полиномам Лежандра. В дальнейших рассуждениях, рассматривая функции (векторы, тензоры) Щ1 ,%3},будем предполагать, что они достаточно гладкие для равномерной сходимости рядов
где Рк -полиномы Лежандра; Р]к -моменты функции f порядкак , которые зависят от координат
ОТ1"?2"?3).
Сведение трехмерной задачи линейной теории упругости к двумерной будем проводить на основе аппроксимации неизвестных функций (векторов, тензоров) в виде отрезков рядов (3.1). Представим I1 и Й , используя (1.3), (2.5) в виде
k=0
US = ZT=0( <u)vK+ w^Pk (3.2)
CO
№ 3.1 (72.1)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2020 г.
где
да,-
i
i
мНу = | цц* а^Р^3, ® i = Pkdf3,
-i -i Q)iy = ö^[Ü]k .W(k) =
= п • [и]к,д0 = д = ^д/д0.
Величины м^уа будем называть моментами тан" 1 (к) генциальных напряжении, к -порядка, а а моментами поперечных (перерезывающих) сил к -порядка.
(к)
Можно показать, что м у а являются компонентами поверхностных тензоров [4] (или 8-тензоров [3])
1 + 2к
Аналогично и Ш(к) будем называть моментами тангенциальных и поперечных перемещении к -порядка.
4. Аппроксимация напряжении. Аппроксимация напряжении проводится в виде отрезков рядов (3.2). Для одних и тех же величин Ь1 используются два вида аппроксимации Т1и ТО , отличающихся только количеством удерживаемых членов в рядах (3.2):
ta (4.1)
2
Т0 = iM)Ya^y +
^а • n)rk,
г)Рк,
(о)
а • п
)Ро,
1
к=0
2
(к)3у •aYPk+ ^—-— q
к=0
I
1 + 2к
(к)
3 • пРк)
Го = •1
1 + 2к
к=0
2
Полагая для простоты / = 0 , условия равновесия усилиИ и моментов элемента слоя можно получить из уравнении (1.1)
п x
Г dt1
1 W
Pkd? =
= 0,(к = 0,1),п X f^Podt3 = 0 (4.2)
Для аппроксимации (4.1) уравнения (4.2) можно записать в форме
*xjfl = 0rfxd£=0. (4.3)
5. Аппроксимация деформации и перемещении. Аппроксимация деформации и перемещении осуществляется на основе энергетического тождества, которое является следствием равенства
xn)(uxn) + (ТО xn)(ux n)}dVf
ч
Г _ д ,_( _ _
= jiTL^('Tx<(uxn)) + T0
д
•^r(nx(uxn))}dV(,
dVj: = d^d^2d^3,
(5.1)
справедливого для любого вектора и и произвольных^, удовлетворяющих (4.3). Для аппроксимации вектора Й в виде отрезка ряда (3.2)
и~и = а?П=о щ Рк +пП=о (%Рк, (5.2)
используя формулу Гаусса-Остроградского и своиства полиномов Лежандра, равенство (5.1) можно привести к виду
[f J3üdS,](Ci + +fLhVoAn-=o>YfY+(y%]dL = f (<?•§+f3§)dV,, (5.3)
где V 0а = %•//; 1/0"- вектор нормали к£ в точках контура Ь и
ü0 = аУ И fc-О (ц^Рк + .
(5.4)
Обозначим правую часть равенства (5.3) через Е и подставим выражения Ь1 (1.3). Получим
Е=5уоЧ ЕцМ.М = J¿dí1dí2dí3,
где Бу - аппроксимации компонент тензора деформации £у
2eaß~2Eaß
= + (5.5)
2£аз~2Еаз = +
dÜo^
Р -I- düo
£33~^33 = (Э3 •
При этом отрезки рядов (5.2), (5.4) определяют две аппроксимации вектора перемещении, одна из которых отвечает производным от и по координатам (а =1,2), а другая- производнои по координате
Г3.
6. Граничные условия. Граничные условия (1.4) аппроксимируем в соответствии с левои частью равенства (5.3). На лицевых поверхностях имеем
(^)±(1-ч) + ^ = ^3 = ±1). (61)
Q
~ <
V.
1
-1
№ 3.1 (72.1)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2020 г.
На боковой поверхности £ получаем, соответственно
1 + 2к
2
(Vayv0aä;(1-q) + (JVäq
1(0^Уоа(1 — я) + ^Я = П[С]0. (6.2)
Равенства (6.1) должны выполняться в каждой точке лицевых поверхностей Б± , то есть при всех 1, £2) Е , следовательно, их необходимо рассматривать как дополнительную группу уравнений относительно неизвестных функций. С другой стороны условия (6.2) накладывают ограничения только в точках контура Ь, и таким образом являются граничными условиями для соответствующей дифференциальной задачи.
7. Закон Гука. Пусть напряжения а1' и деформации Бу , соответствующие (5.5), связаны законом Гука (1.2):
-пх ([G]fc хп)к- 0,1,
е-М*-f+^'f)
(7.1)
Соответственно, для коэффициентов разложений, исходящих в отрезки рядов (3.2), получим соотношения
где к=0,1; если 1=1,2 ;к= 0,1,2, если 1=3;
)Pkd^3, (7.2)
где к=0,1; 1=1,2,3, которые в общей сложности составляют 20 уравнений.
8. Двумерная краевая задача. Обобщая полученные выше результаты, сформулируем двумерную краевую задачу: найти 31 неизвестную функцию
(к)
1)ау (к - 0,1; а - 1,2; у
- 1,2); ™3Y(k - 0,1,2; Y - 1,2); qi(k- 0,1; i - 1,2,3); (k)ua(k - 0,1,2,3; а
(к).
- 1,2),™ (к - 0,1,2)
из 31 уравнения, в том числе 5-ти уравнений равновесия (4.1); 20-ти уравнений законов Гука (7.2) и 6-ти уравнений (6.1). Обозначим
Na - ^dö((0)aY ' К + (> ' п),Ма
j
Р±-Т3\
V -"Ja-аа+ w
S - Ыгу.ТГ wn, U - <а- aa,
e=±1,U±-Ü\ е=±1, (8.1)
и перейдем к новым неизвестным функциям N а,Ма,У,Р±,в,и± . При этом имеют место равенства
U;-v + 6-i3,
dU
2—- U+ - U- + 3P1(U+ + U- - 2V)
o^3
+ 5P2 (n х ((U+ - TT-) хп) - 2в)
(к)
.Г (dU; ~ ди\
-njicrß-w+cr3-w)p>
Сформулированная выше двумерная задача в новых неизвестных сводится к системе дифференциальных уравнений
1
U
+ пх(ИахЭ0^+Щ;х(Махп) + к-пх ((р+ — Р-)хп)-0,
dNc
+ Р+ -Р- -0,
№
ма - JFoäy t ß (äy-&-Щ + äyC~3 • g) Pidt3
'p± - JToäy tß(äy -c~ß-§ + äyc~3 • Ü)у±^3 +
+
-CW-jß + n- C3a3 • JJJ) <p±;df3,
dJß ■ " dJ3J~r o ■(1-q)+ Jßq - G±,
Uee33)
Где cp± - ZEtol-+r(±1)kPk,l - 0.1.
(8.2а)
(8.2б)
Назовем основными неизвестными величины ! ,М" ,7 , 8 , которые входят в уравнения (8.2) вместе со своими первыми производными по координате . Остальные неизвестные Р± , и± имеют смысл
напряжений и перемещений на лицевых поверхностях и входят в уравнения без производных, их будем называть дополнительными. Соответственно,
V
V
№ 3.1 (72.1)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
первую группу уравнении (8.2а) будем рассматривать как основную, а вторую (8.2б) - как дополнительную. Решая дополнительную линеиную систему относительно дополнительных функции, можно получить их выражения через основные неизвестные. Далее, подставляя наиденные Р±, и±, в первую группу уравнении, получим краевую задачу для основных неизвестных Ыа,Ма,7, в, которая будет
март, 2020 г.
иметь 10-ьш порядок. Соответствующие краевые условия определяются из соотношении (6.2).
Отметим особо справедливость энергетического равенства (5.3), которое является основои для обоснования корректности двумернои задачи, доказательства теоремы единственности, формулировки вариационных принципов.
Список литературы:
1. Иванов Г.В. Теория пластин и оболочек. Новосибирск. 1980.- 84 с.
2. Алексеев А. Е. Уравнения плоского изгиба упругого слоя переменной толщины / А. Е. Алексеев. - С .3-16 Динамика сплошной среды : сборник научных трудов / Акад. наук СССР, Сиб. отд-ние, Ин-т гидродинамики им. М. А. Лаврентьева. - Новосибирск : Наука, Сибирское отделение, Вып. 75 , 1986
3. Векуа И.Н. Некоторые общие методы построения различных вариантов теории оболочек.-М.: Наука, 1982.
4. Naghdy P.M., Foundation of Elastic Shell Theory.-Progress in Solid Mechanics 4, N2, p.1-90.
