К проблеме разложения по собственным вектор-функциям в нестационарных начально-краевых задачах динамики оболочек вращения

Приложения дифференциальных уравнений

УДК 539.3: 517.956

Ю.Э. Сеницкий, А.Ю. Сеницкий

К ПРОБЛЕМЕ РАЗЛОЖЕНИЯ ПО СОБСТВЕННЫМ ВЕКТОР-ФУНКЦИЯМ В НЕСТАЦИОНАРНЫХ НАЧАЛЬНО-КРАЕВЫХ ЗАДАЧАХ ДИНАМИКИ ОБОЛОЧЕК ВРАЩЕНИЯ

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

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

1. Поскольку формальная процедура предлагаемого метода является новой [5], то сначала сформулируем в общем виде матричный структурный алгоритм биортогональных КИП применительно к линейной начально-краевой задаче, описывающей, в частности, нестационарное осесимметричное напряженно-деформированное состояние (НДС) оболочек вращения. Для этой цели в классе вектор-функций с интегрируемым квадратом в метрике сепарабельного вещественного пространства L2p (I) на сегменте I = [ü,01 ] области W : {I, x, t > 0} построено следующее биортогональное КИП [5]:

01

,t)

0

¥

u

i=1

Здесь u - искомая вектор-функция (матрица-столбец порядка mx1); K, G - две ядровые вектор-функции того же порядка, являющиеся ненулевыми решениями инвариантной и сопряженной однородных краевых задач, а , ßj - соответствующие им собственные значения; p(ß^)> 0 -

диагональная матрица

(k = 1,2...m).

Для упрощения дальнейшего рассмотрения включим p(ß) в ядровые вектор-функции, т.е.

К, =[Р(в)}',4K(m„в), С, =[р(#G(1,6), ul = [р(е)Р2и(Є„t),

и сохраним для К; , Оі , и1 прежние обозначения. Таким образом, имеем

6,

,і) =

0

а в (2) вместо скалярного произведения (К,, С,) следует считать (К,, с,).

Формула обращения (2) основана на таком соотношении биортогональности:

где 5- - символ Кронекера.

Было показано [5], [6], что при условии ограниченности трансформант (I) обеспечивается

~ т2 ~

единственность представлении и сходимость в метрике пространства Ьр разложении, определяемых формулой обращения (2).

Рассмотрим начально-краевую задачу, описываемую системой дифференциальных уравнении гиперболического типа:

А[и ]- Еи = /

и

і=0

ё,и + е,и = 0, 6 = 0;

ё2и + е2и = 0; 6 = 6, ,

где А [и ] - дифференциальный оператор вида

А[и] = аи + Ьи + си,

а, Ь, с, сіі2 , е 12 - квадратные матрицы, причем первые три из них функциональные, а - диагональная, Е - единичная и а, ЛІі2 - невырожденные матрицы. Все они имеют порядок (тхт).

/, и0, и0 - известные вектор-функции (тх1). Штрих и точка обозначают дифференцирование по пространственной переменной 9 и времени і. Соотношения (5), (6) представляют начальные и граничные условия.

К (4) приводится исходная система дифференциальных уравнений движения после деления на диагональную неособенную матрицу коэффициентов при инерционных членах. В случае,

когда оператор А сингулярный на одном из концов интервала, например 6 = 0, то соответст-

вующее равенство (6) заменяется на условие ограниченности (регулярности решения).

Если теперь подействовать интегральным оператором (1) на уравнения (4) и начальные условия (5), то в результате получаем

л2

После выполнения интегрирования по частям первого слагаемого (8) с учетом (7), удовле-

творяем обычным двум условиям структурного алгоритма метода КИП [2], т.е.

Ф*", о]0 =1"'Уатв - "т

Здесь А сопряженное А дифференциальное выражение

А* [О ] =

Равенство (10) представляет обращение в нуль билинейной формы на концах интервала, а (11) - операционное свойство.

Имея в виду обозначение трансформанты (11) и соотношения (10), (11), равенства (8), (9) образуют счетное множество задач Коши относительно (р(1,,,):

j(l >t) +12,t) = -F, t), i = 1,¥ ,

,0) = jo (1), j(l, t) = j 0 (li) > t=0

I t=0

Ее решение записывается в виде

t

j

где F,ф0,ф0 - соответствующие трансформанты (1х), т.е.

Ыя,,/)

^1 2 е1 2 ,

в =0, 0,

-'0 1 J

Из операторного равенства (11) следует такая система дифференциальных уравнений для компонентов вектор-функций О ядра КИП:

А* [О] = -Я2О .

Выражая из (6) производные и и подставляя их значения в (10), получаем следующие граничные условия для О:

1)атО + g1О = 0, в = 0; 2)атО + g 2 О = 0, в = в1.

Здесь

g1,2 = (а ^ + (й1,2 е1,2 ^ а — Ь . (17) Соотношения (15), (16) представляют сопряженную краевую задачу для определения компонентов Ок вектор-функции ядра О КИП. В случае, когда

ат\ = й,2 , (а Г — Ь

в=0, в1 1-2 !' '

причем матрицы йг<2 - симметричные, граничные условия (16) являются самосопряженными, т. е.

1)й1О + ехО = 0, в = 0; 2)й 2О + е2О = 0; в = в1 .

Составляем теперь скалярное произведение для членов равенства (15) с ядровой вектор-функцией А”. Имеем

Выполняя интегрирование по частям и используя аналогичные (10), (11) два условия

Ф ** [От, К ] 01 =[

из (19), (21) получаем соотношение биортогональности (3):

(Я — т )(о, К )= 0, , * ].

В операторном уравнении (21) А представляет сопряженное к А дифференциальное выражение [7], т.е.

А**[К] = А[К] = аК + ЬК + сК .

Если теперь воспользоваться операционным свойством (21), а затем равенством (20) с учетом граничных условий (16), то получаем инвариантную исходной однородную краевую задачу для компонентов Кк второй ядровой вектор-функции КИП К. Имеем:

А [К ] = — т? К,

1)й1 К + е1 К = 0, в = 0; 2)й2К + е2 К = 0 , в = в1.

После вычисления собственных вектор-функций О и К и соответствующих им собственных значений Яi, Л] из однородных систем обыкновенных дифференциальных уравнений (15),

(16) и (23), (24), общее решение и рассматриваемой начально-краевой задачи (4)-(6) находится по формуле обращения (2) с учетом (13), (14). Таким образом, возможность построения замкнутого решения задачи (4)-(6) определяется интегрируемостью систем обыкновенных уравнений (15) и (23).

i

i=1

где Gil2 =(G , Gt ) p - квадрат нормы в метрике L2p .

и (3) при этом

Замечание. Если К= G, то биортогональное КИП (1), (2) становится обычным КИП [2], причем трансформантой по-прежнему является скалярное произведение (1), а формула обращения принимает такой вид:

¥

,W„ „-2

i=1

112 =(G. G /р ____^......F________^

Соотношение биортогональности (3) при этом превращается в известное условие ортогональности

-X =

КИП (1), (25) предназначены для исследования самосопряженных начально-краевых задач (4)-(6), когда

A* [G ] = A[G ],

и справедливо тождество Лагранжа [7]

2. Будем рассматривать в сферических координатах (в, j) упругую оболочку вращения, имеющую толщину h и главные радиусы кривизны срединной поверхности R1 (в ), R2 (в ). Известно [8], что для поверхности вращения с произвольным меридианом коэффициенты ее первой квадратичной формы и соотношение Кодацци-Г аусса вычисляется по формулам

A1 = R1, B1 = R2 sin в ,

При учете поперечных сдвигов по Р. Миндлину [9], а также выражений (28), геометрические и физические уравнения для случая осесимметричной деформации определяются следующими равенствами [10]:

eв = ее + zCe , = ej + zcj, g0 = y + R1-1 ,

e° = R1-1 (и* + W*), ej = R2-1 (u*ctgв + W*), cl = R^y' , cj = R2-1yctgв ,

N* = cR-1

Nj = c[r2-1 (и*ctg в + W*) + vR— (и* + W*)], M j = d(r2-Vctg в + vR1-1y'),

Q*e = 0.5^1c(1 - v)[y + R1-1 (w*' - и*)].

Здесь и*, W* - тангенциальная и нормальная компоненты вектора перемещений; у - угол поворота в меридиальной плоскости; eв ,ej,g0z - относительные линейные деформации и угол сдвига, а св, Cj - соответствующие изменения кривизны, причем индекс «0» относится к срединной поверхности оболочки; z - координата в направлении толщины; k¡ - коэффициент поперечного сдвига; N*,Nj,Q* - нормальные и поперечное усилия; Mg,Mj - меридиальный и кольцевой изгибающие моменты;

C Eh D Eh3

C=(Tv) • D=T2FV^

E,v - модуль упругости и коэффициент Пуассона.

Для оболочки вращения загруженной произвольной осесимметричной динамической нагрузкой X* (в, t), Z * (в, t), M * (в, t) дифференциальные уравнения движения с учетом инерции поворота поперечных сечений записывается в виде [10]:

( U2 \ ..

(n* R2 sin в )'- Nj (R2 sin в ) + Q* R2 sin в - ph|

1 +- h

R1 R2 и* sin + X * R1 R2 sin 0 ;

Q* R2 sin в ) - Nj R1 sin в - N* R2 sin в - ph 1 +-------------R1R2 W* sin в + Z * R1R2 sin в =0;

12R R2

1 2 /

(m* R2 sin ^ ) _ Mj (R2 sin в ) - Q* R2 sin в

12

ph3 1 + (320)h 2

R1R2 y* sin в + M * R1R2 sin в = 0,

R1R2

где р - плотность материала оболочки.

Вводим безразмерные величины по формулам

и и* Rimax:

W = WR-1

" "*JV1max:

t = t* R11

^ imax

P(1 -v1

P1 = r1 R1-max ,

P2 = R2R^ax, X = X*R2

Ne = N*C 1 = p11 (u + W)+vp21 (uctge + W), Nj = NjC-1 = P2-1 (uctg в + W) +vP- (u + W) ,

*

'0 _ J,el

Г*(

' j - j

M в = M * R1max D ^ = РГУ + VP 2 Ctg в , Mj = Mj R1max D - = P2-Vctg в +VPlV ,

<Р

о„ = о„*с= к у + р,-1

притом Д|ш,х = тах Я,

Пусть замкнутая в вершине оболочка вращения опирается на упругий относительно углов поворота и линейных смещений контур, коэффициенты жесткости которого ¡л*, /3*,у*. Соответствующие этому случаю закрепления граничные условия формулируются следующим образом:

Ж* <¥ , и* <¥ , у < ¥ , 0 = 0,

N = Р*и*, о* = 7>*, М* = лУ , 0 = 0| .

Первые три неравенства являются условиями ограниченности в полюсе оболочки = 0 . После подстановки выражений (29) в уравнения (30) и условия (32) с учетом (31) и последнего равенства (28) приходим к начально-краевой задаче (4)-(6). Получающиеся при этом

sk II

c = \\c

sk

имеют следующий вид:

(в, t) = и (в, t)W(в, t)y (в, t)|T, f (в, t) = ||X(в, t)Z(в, t)M(в, t)|T ,

u,

(в )= ||U0 (в )W0 (в У0 (в )T , u 0 (в ):

u 0 (в )W 0 (в )y0 (в

P- 0 0 (VP2-1 ctg в1 - ß ) (P1-1 + vP2-1 J 0

d 2 0 kP1-1 0 , e2 - £ - k

00 P1-1 00 c t OQ 1 - )

a11 = a33 = k la22 = p12, ask = 0, (s ^ k), b11 = b33 = k lb22 = p12 (ctg в 1 -1 ~ 1

22

b12 = b21 = ГГ2 (1 + k + VP1 P—1

c21 = p

'22 ~ P1~ (ctg в + P2 P2‘ - P11P1-1 ) , l), b13 = b31 = 0, b23 = -a2 b32 = kp1-1 ,

cn + kp1-2 = C33 + ka~2 =-(p1P2 )-1 V + ctg2 в + P2 P-1 ctg в ),

12 = P1-2 [

-P1-2 [k(P1P2-1 + P^P1-1 + ^в)+ P1P2-1 v + P1P2-1 )ctgвJ, c13 = a2c31 = kPl-1,

kP1-1 (P2P1-1 + ctg в) , c32 = 0, c22 =-(P1P2 )-1 (P 2 P1-1 + P1P2-1 + 2v).

Элементы матриц (34) постоянные величины, т.к. P1 = P1 (в1 ), P2 = P2 (в1 ) .

В соответствии с (32) вместо первого равенства (6) следует иметь в виду условия регулярности решения, т.е.

u

c23 =

T

f (0 ) 0 0 Р1 p2 sin 0 00

Р(0 )= 0 f (0 ) 0 = 0 p1 р2 sin 0 0

0 0 a2 f (0 ) 0 0 a2p1 р2 sin 0

Соотношения (4)-(6), (36) при наличии (33)-(35) и известных матрицах f, u0, u0 представляют математическую формулировку рассматриваемой начально-краевой задачи в безразмерной форме. m - безразмерные коэффициенты жесткости опорного контура относительно

линейных смещений и углов поворота

b = ЬХ-С“1, у = у*RimaxC_1, m = rn*Rlmaxc-1.

3. Воспользуемся биортогональным КИП (1), (2), в котором весовая функция f (в) , зависящая от коэффициентов системы уравнений (4), (34), (35), определяется с точностью до постоянной такой квадратурой [2], [11]:

f q ) = exp(я- Jbkkde)= exp|p!2 q )JPl2 q )(ctgp + p\ q)p- q )- pj q )p- q ))deJ= Piq)p2 q)sin в , а ее матрица р(в ) соответственно

Если теперь подействовать интегральным оператором (1), (37) на уравнения (4), (34), (35) и начальные условия (5) в соответствии с алгоритмической процедурой пункта 1, получаем счетное множество задач Коши, разрешаемое относительно трансформанты <р(Яі, Ї) в виде (13), а также сопряженную краевую задачу (15) для вектор-функции ядра КИП 0(ЯІ, 0) и граничные условия (16), т.е. имеем:

а \\ (0 К + Ь„ (0 £; + оп (0 £, + Ъ„ (0 £; + 4 (0 )02 + с„ (0 )0, = -Я;0„

а 22 (0 К + »22 (0 £; + С 22 (0 )Сг + »21 0 К + 4 0 К + »23 (0 К + €„ (0 = -я? а,

а33

01 (о) < ¥, 02 (о) < ¥, G3 (о) < ¥, 0 = 0,

01 + р1 (ур201 — + (1 + уар—1 К = 0, 0 = 01,

К2 — К1 — /к 1 р1С2 + р1К3 = 0,

о3 + р1 ур—01 — ^ = о,

где 0(яг, 0 ) = 1^1 (Яі, 0 К (Яі, 0 К (Яі, 0 )|г

[(1 — г2 р—2 0 — р^ р;1 + (1 + к )р\р21 + (1 + 2£ )р11р1—1 ]

С12 = рГ21

с2*1 = -РГ2 tk -V + (2v + Орр-1 Jctg 0 + (k +v)pp-1 + кРІРі-1 + PiPp-2},

а в выражении (13) для <р(Яг-, t) соответственно имеем:

01

F (l, t) = J [Х (0, t )G1 (l, 0 ) + kZ (0, t )G2 (l, 0 ) + M (0, t )G3 (l, 0 )Jp1-1 p2 sin 6Ш ,

p1 p2 sin 0J0 .

j0 ) j0 (l) 01" = 0 0 0 0 0 Gi (Я, ,0 )+ W (0 ) W 0 (0 ) G 2 (Я0 )+ W00 ) v/00 І G3 (Я,0 І

Хотя при соответствиях и ~ С1, Ж ~ G2, у ~ G3 краевые условия (39) совпадают с (6), (34), (36), однако из равенств (35), (40), (41) следует, что

* * С12 * С12 , С21 * С21

Следовательно дифференциальный оператор А [и] для оболочек вращения с произвольным меридианом несамосопряженный, причем левая часть уравнений (38) совместно с (39) представляет А ^].

Применяя к дифференциальным уравнениям (38) аналогичное (1) КИП с весом (37) и ядровой вектор-функцией

к (т, в )=|к, (т, в )к (т, в )к, (ц,, в)г,

в соответствии с процедурой (19), (21) и учитывая (39), находим

L2 + C13 (ß )К 3 = Hi K1 ,

b23 (ß )К3 + C23 (ß )К3 = m K2 ,

ап (в К1' + ¿11 (в )К1 + Сц (в К + ¿12 (в )К2 + С12 (в )К2 °22 (в )К2 + ¿22 (в )К2 + С22 (в )К2 + ¿12 (в )К1 + С21 (в )К 1

а33 (в)к3 + Ь33 (в )к3 + с33 (в )к3 + с31 (в )к1 + Ь32 (в )к2 = -тг2К3.

Граничные условия для К(т, в ) при замене G1, 02 , 03 соответственно на К1, К2 , К3 тож-

дественно совпадают с (39), и вместе с (42) образуют оператор А**[К] = А[К]. В результате интегрирования однородных краевых задач (38), (39) и (42) с аналогичными граничными условиями получаем биортогональный базис G и К, а затем по формулам обращения (2) с учетом равенства (13) искомые функции перемещений и(в, ^), Ж (в, ^), у(в, ^). Таким образом, показана принципиальная возможность решения нестационарных задач динамики для оболочек вращения путем разложения в ряды по собственным вектор-функциям.

4. В качестве примера рассмотрим загруженную произвольной осесимметричной динамической нагрузкой упруго-защемленную относительно углов поворота непологую сферическую оболочку с конечной сдвиговой жесткостью. Для нее имеем

Я1

С учетом (43) матрицы коэффициентов (34) системы уравнений (4) п

эинимают такой вид:

1 0 0 ctg (1 +v + k) 0

a = 0 k 0 , b = -(1 + v + k) kctg k

0 0 1 0 - ka -2 ctga

с =

-

-

ka 0

k

kctg

-

Краевые условия (5), (33) сохраняются, а (6), (36), (32) для соответствующего закрепления на контуре формулируются следующим образом:

{u,W,y}<¥, 0= 0, и = W = 0, y +

Соотношения (45) получены из (6), (32) путем деления первых двух равенств на ß, у и предельных переходов ß  ¥ , g  ¥ .

Применяя к уравнениям (4), (44) и начальным условиям (5), (33) КИП (1) с весовой функцией f (0 ) = sin 0 по структурному алгоритму (8)-(24), находим A*[G\ = A**[K\ = A[G\.

Таким образом, начально-краевая задача (4), (5), (44), (45) является самосопряженной, т.е. G=K, и ядро преобразования есть нетривиальное решение следующей однородной краевой задачи (15):

L(G1) - (v + k + ctg2 0 )G + (1 + v + k )G2 + kG3 + A2 G1 = 0,

L(G2)-(2 +v)G2 -(1 + v + k)(g1 + G1ctg 0)+ k(g3 + G3ctg0)+ A2G2 = 0 L(G3)- (у + ka-2 + ctg2 0 )G3 + ka - (G1 - G2)+ A2G3 = 0,

{G1, G2, G3 }<¥, 0 = 0, G1 = G2, G3 + (vctg 01 - m)G3 = 0, 0 = 01.

Здесь Z(...) = —(-) + ctg в ("') . Трансформанта j(li,t), как и прежде, определяется равен-

de de

ствами (13) и (14), (33) с весом sinö.

Проинтегрируем систему дифференциальных уравнений (46). Для этой цели вводим потенциалы смещений U и усилий V по формулам

G1 = U', G3 = k-lV + U' - G-,

а также замечаем, что

_d_

d0

L(.-.)= L

d

d0

- sin -2

d

d0 '

После подстановки соотношений (48) в (46) с учетом (49), выражаем из второго равенства

G2, т.е.

G2 =-k„L(U + k 4iV).

В результате приходим к такой системе уравнений:

L(U)+кш U + k„,L(V)+ku,V = 0,

L2

исключая из которой U получаем разрешающее дифференциальное уравнение для V. Имеем

L

В процессе приведения (51) к (52) установлена связь U с V, т.е.

U = -к21,L2

В соотношениях (50)-(53) приняты такие обозначения:

к1, = 1 — V + 1 , к 2, = 1 + V , к3, = Ь к 4, =-(1 +V)-1. к5, = (1 + V)-1 [2(1 + v) — 12 J k6, = кЛ k7, = k-1 [(1 — v) — ka - + Я2 J, k8, = -k9l = ku , кш, = k1,k5, (k5, — k2, )-1 , k11, = k2,k4, (k5, — k2, ) , k12, = k3,k5, (k5, — k2, ) , k13, = k5, — k9, , k14, = k5, k8, ,

k15, = k5, k6, — k4,k9, , k16, = k5, k7, ,

k17, = к14, (к13, — k10, ) , k18, = (k4, — k11, Xk13, — k10, ) , k19, = (k15, — k12, Xk13, — k10, ) ,

k20, = k16, (k13, — k10, ) , k21, = k18, (k10, — k17, ) , k22, =(k11, — k19, )(k10, — k17, ) ,

k23, = (k12, — k20, Xk10, — k17, ) , k24, = k21, (k22, — k18, + k17,k21, ) ,

k25, = k21,

Вводим порождающее дифференциальное уравнение

L

В результате подстановки (55) в (52) формируется бикубическое уравнение относительно

xn:

Если учесть, что при x=cos9 (55) является уравнением Лежандра, и из (56) определить корни X , X, 2, X, 3, а затем сп, с, 2, с,3, то общее решение уравнения (52) можно представить в виде

V = Z (cos 0)+д„бсг„ (cos 0)J,

n=1

где Pc,n (c0S 0 ) , бсш ( cos 0) - функции Лежандра степени c,n первого и второго рода; C,n , D,n -произвольные постоянные.

Принимая во внимание условия регулярности (47), а также соотношения (53), (48), (50), (55), определяем компоненты вектор-функции ядра КИП:

G11, 0 )=—Z E,nC,nPL (cos 0 ), G2 (1, 0 ) = —Z F,nC,nPL (cos0 ),

n =1 n =1

G3

n =1

где P1n - присоединенная функция Лежандра первого порядка степени с ,n.:

En = k 21 ,i,n — k 22Л,n + k F,n = k[k 21 Л!, — k 22,X‘n +(k 23, — k 4, )E,n J.

H„ = k+ k 4,k -X —(1—k—л, E .

Дальнейшее решение очевидно, поскольку рассматривается задача на собственные значения. В результате подстановки равенств (57) в оставшиеся граничные условия (47) формируется однородная система алгебраических уравнений для C n:

Z CmEmPL (cos 01) = 0; Z CmFmPL (cos 01) = 0;

n=1 n=1

Z CmHm i Си(Си0+ ^ PL+1 (cos 01 ) + [(v — cm — 1)ctg01 — mPn (cos 01 )} = 0.

n=1 I sin 01 J

Из условия нетривиальности решения системы (58) составляется трансцендентное уравнение относительно 1 :

Det (l ) = 0 ; 1 = 1, ¥

и находятся постоянные Cin. Имеем:

сп = E, 2 Fil Pin (c°s qi )P„ 3 (cos qi)- Et 3F 2 p\3 (cos qi )PC12 (cos qi ),

C,2 = Ei3F1 PC'3 (c0S qi )Pc1 (c0S qi ) - Ei1F3PCl (c0S qi )Pc3 (c0S qi ) 5

C3 = Ei1 F,2PCn (c0S qi )Pc2 (c0S qi ) - Ei2FPC-2 (c0S qi )Pc1 (c0S qi ) .

Применяя теорему обращения к выражению (13) и учитывая равенства (57), определяются искомые компоненты вектора перемещений

¥

{«(в, t ),W (в, t )у(в, t )}=£j(l. t ){G,(1, в ),G2 (1, в ),GS (1, в )J|G,r3,

i=1

а затем по формулам (31) внутренние усилия в оболочке.

Полученные результаты справедливы и при исследовании неосесимметричных начальнокраевых задач динамики оболочек вращения с конечной сдвиговой жесткостью.

Широкий класс несамосопряженных линейных краевых задач динамики образуют также математические модели неоднородных или многослойных тонкостенных конструкций [5], [12], для решения которых эффективным представляется рассмотренный здесь метод биортогональ-ных КИП.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Улитко А.Ф. Метод собственных векторных функций в пространственных задачах теории упругости. Киев: Наукова думка, 1979. 264 с.

2. Сеницкий Ю.Э. Исследование упругого деформирования элементов конструкций при динамических воздействиях методом конечных интегральных преобразований. Саратов: Издательство СГУ, 1985. 176 с.

3. Мартыненко Н.А., Пустыльников Л.М. Конечные интегральные преобразования и их применение к исследованию систем с распределенными параметрами. М.: Наука, 1986. 393 с.

4. Сеницкий Ю.Э. Многокомпонентное обобщенное конечное интегральное преобразование и его приложение к нестационарным задачам механики // Известия высших учебных заведений. Математика. 1991, №4. С. 57-63.

5. Сеницкий Ю.Э. Обобщенные биортогональные конечные интегральные преобразования и их приложение к нестационарным задачам механики // Доклады РАН. 1995. Т. 341. №4. С. 474-477.

6. Сеницкий Ю.Э. Сходимость и единственность представлений, определяемых формулой обращения многокомпонентного обобщенного конечного интегрального преобразования // Известия высших учебных заведений. Математика. 1991 №9. С. 53-56.

7. НаймаркМ.А. Линейные дифференциальные операторы. М.: ГИТТЛ, 1954. 352 с.

8. ГольденвейзерА.Л. Теория упругих тонких оболочек. М.: Наука, 1976. 512 с.

9. Сеницкий Ю.Э., Еленицкий Э.Я. О физически непротиворечивой модели уточненной теории пластин и оболочек. // Доклады РАН. 1993. Т. 331. №5. С. 580-582.

10. Григолюк Э.И., Селезов И.Т. Неклассические теории колебаний стержней пластин и оболочек. // Механика деформируемых твердых тел. М.: ВИНИТИ, 1973. Т. 5. 272 с.

11. Сеницкий Ю.Э. О некоторых тождествах, используемых при решении краевых задач методом конечных интегральных преобразований // Дифференциальные уравнения. 1983. Т. 19, №9, с. 1636-1638.

12. Колесников С.В. Уточненная теория колебаний многослойной ортотропной пластины // ПММ. 1993. Т. 57. Вып. 5. С. 160-165.

Поступила 28.05.2004 г.