Бесконечные системы линейных уравнений в задачах устойчивости прямоугольных пластин Текст научной статьи по специальности «Математика»

РУБРИКА «МАТЕМАТИКА»

БЕСКОНЕЧНЫЕ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ В ЗАДАЧАХ УСТОЙЧИВОСТИ

ПРЯМОУГОЛЬНЫХ ПЛАСТИН

Папков Станислав Олегович

канд. физ.-мат. наук, доц., зав. кафедрой «Высшая математика», ФГАОУ ВО Севастопольский государственный университет,

РФ, г. Севастополь E-mail: SOPapkov@sevsu. ru

Довгаленко Ольга Александровна

старший преподаватель кафедры «Высшая математика», ФГАОУ ВО Севастопольский государственный университет,

РФ, г. Севастополь

Чевын Майя Анатольевна

старший преподаватель кафедры «Высшая математика», ФГАОУ ВО Севастопольский государственный университет,

РФ, г. Севастополь

INFINITE SYSTEMS OF LINEAR EQUATIONS IN BUCKLING ANALYSIS

OF RECTANGULAR PLATES

Stanislav Papkov

candidate of Science, Head of department of Mathematics, Sevastopol State University

Russia, Sevastopol

Olga Dovgalenko

assistant professor of department of Mathematics, Sevastopol State University

Russia, Sevastopol

Maya Chevin

assistant professor of department of Mathematics, Sevastopol State University

Russia, Sevastopol

АННОТАЦИЯ

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

ABSTRACT

A new analytical solution for buckling of thin rectangular plate with clamped edges is presented. The problem has reduced to a quasi-regular infinite system of linear algebraic equations. Asymptotic behavior of unknowns is constructed for non-trivial solution of system. The accuracy and computational efficiency of the approach are demonstrated by examples.

Ключевые слова: устойчивость, аналитическое решение, бесконечная система линейных алгебраических уравнений.

Keywords: stability, analytical solution, infinite system of linear algebraic equations.

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

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

Рассмотрим тонкую пластину

(х,у) е{[-а; а]х [-6;6]} равномерно сжатую в двух направлениях под действием контурных нагрузок Ых и Ыу. Функция прогиба ^(х,у) должна удовлетворять линеаризованному уравнению устойчивости:

д2 w д2 w DAAw + N —w + N —w = 0 С1)

т дх2 ' ду2

На границе пластины Г заданы условия жесткого защемления:

I п dw wl = 0; —

|Г дп

= 0.

(2)

где D = -

Eh3

цилиндрическая жесткость, v ■

12(1 -v2)

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

Пусть, Q = N / D; P = N / D тогда общее решение уравнения (1) может быть представлено в виде четной функции по обеим координатам

да

w = X (An ch РшУ + Bn ch Рг.пУ) cos anx + n=1

да

+X (Cnch qh,nx+Dnch qi,nx) cos py

n=1

где an = л (n -1/ 2) / a; Pn = л (n -1/ 2) / b ; величины Р; n и „ (j = 1,2) являются корнями характеристических биквадратных уравнений

p4 - (2аП - P)p2 +а4п - Qa2n = 0, q4 - (2Д2 - QW +P4 - PPI = 0

для j = 1 выбираются со знаком «+», и для j = 2 со знаком «-»:

Pj,п= а -P/2 + V(Q-Р)а2„ + P2 / 4, q, = = .\pl -Q/2±7(P-Q)Pl + Q2/4

. (5)

Краевые условия на функцию прогиба пластины V можно выполнить точно:

да

Ч=+а = X (спсЬ +впсЬ ч2,„а) 008 РпУ =0;

И=1

ly=±b

если

= X( Anch Pl,nb + Bnch P2,nb)c0sanx = 0,

ch p b ch q, a

B =-A 'n ; D = -C. —. (6)

n n

ch Pl,nb'

nn

ch W2,na

Условия на нормальную производную, после подстановки (6) и разложения входящих в равенства функций в ряды Фурье, приводят к однородной бесконечной системе линейных алгебраических уравнений относительно неизвестных Хт, , связанных с коэффициентами Ст, Ат формулами:

(-1)n 9

X = Cn ^(ql - ql)ch q^a, Yn =

n

n+1

aPn

= An P2n -P2,n)chp^b

ba„

(7)

Бесконечная система алгебраических уравнений имеет вид

Xm A m =X

2аХ

n=1 Pm (аI + qtm )(a„2 + q22,m ) '

Y A у =X

m m / ,

2p Xn

; (8)

n=1 am (P + PL )(P + PL)

(m=1,2,...)

где

Am = q1,mth qt,ma - qi,mth qi,ma .

a pm (q2m - q{m) ;

Am _ P1,mth P1,mb - P2,m%2,mb

b am(Plm - Pl,m) .

Заметим, что при выводе бесконечной системы использовались соотношения теоремы Виета для уравнений (4), которые, в частности, дают тождество

(Д2 + рЬп )(Д2 + р2т) = («^+ ч2п )(«т + ч2,п)

n=1

Г

Далее, при помощи замены переменных

-1 = X,

= систему (8) можно записать в

каноническом виде

= Е(m =1, 2, — )

(9)

Бесконечная система линейных алгебраических уравнений называется регулярной, если для любого номера m:

^=Е1 мтп\< 1.

(10)

При выполнении условий регулярности (10) однородная бесконечная система (8) имеет нулевое главное решение, опираясь на теорему П.С.Бондаренко [2] можно показать единственность этого решения. Если условия (10) выполняются при т > > 1, то бесконечная система является квазирегулярной и ее исследование может быть сведено к конечной системе линейных уравнений порядка

Для системы (8) значения рядов (10) могут быть вычислены точно на основе известной формулы

Е

1

Я . = — Шяа.

1(п -1/2)2 + а2 2а

Откуда, учитывая, что при достаточно больших номерах m все коэффициенты бесконечной системы (8) положительны, можно найти, что = £2(П = 1.

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

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

*„ = X (1 + 8^), у = ¥я (1 + 81я)

После чего система (8) принимает вид

(11)

Хт дт=Е

(1 + 8т )2«п2 Уп

^ Рт (а1 + Ч21,т )(«п2 + ?22,т ) (1 + 81т К У

(12)

у д У = Е

т т / •

Рт («12 + Ч\,т )(«12 + ^т ) (1 + 8Хт )2Д2 Хп

=2 ат (Р2 + Р1т )(Р1 + Р{т )

(1 + 8\т )РРХ

«т (Р + Кт )(Р1 + Р2,т )

(т = 1,2,...)

Тогда оценка регулярности системы (12) дает:

Е

2а

"2 ДХРт («2 + Ч1т )« + 91т )

а

ДХ Рт («1 + 91,т Ж + 92,т ) = 1--. . . < 1;

Е

ДХРт («12 + Ч1т )(«12 + 42,т )

2Р2

=2 ДУ«т Р + р^ Ж + Р1т )

Р2 .

ДУ«т (Р12 + Р2т)(Р12 + Р2,т)

Р2

= 1 -

Д^т«т (Р12 + Р2т)(Р12 + Р2,т )

< 1,

при всех т > N. > 1, где Nr - номер, начиная с которого все коэффициенты системы положительны.

Из последней оценки следует, что бесконечная система (12) является квазирегулярной.

Для записанной в форме (9) системы (12) выразим через первые ^ = 2N неизвестных остальные неизвестные

= Е &

- т Е= (т > N ) .

1 = 1

(13)

После подстановки (13) в (9) при т > Nк получаем совокупность бесконечных систем относительно {& }

да

& =Е Мтп&п + мт (т >^;] = 1,2,...^) 14)

Системы (14) в силу квазирегулярности исходной системы будут регулярными, что означает для них существование единственного ограниченного решения для любых значений нагрузок P и Q. Первые неизвестные определяются как решение конечной системы

^ =ЕКтп?п (т = 1,2,..., N):

(15)

где Я = М + Е М &.

^ тп тп / 1 т)~ ]

!=NR +1

Система (15) на значениях критических нагрузок Pc и Qc имеет нетривиальное решение, которое позволяет построить собственную форму - форму потери устойчивости пластинки.

Для эффективного решения совокупности регулярных систем (14), заменой переменных

! = !Р2; ^ = &,«X (т > N) (16) регулярные системы (14) приводятся к виду,

+

7

г

т

п=1

п=1

п

+

п=1

ж

К= Еад+4 {т>МЛ,] = \,г...МЛ) (17)

приближенная формула [1], полученная как первый член приближения по методу Галеркина:

в котором удовлетворяют достаточному условию существования ненулевого предела у решения [2]

^ 22ш-1 = lim = < ■

(18)

Наличие положительного предела (18) у решений (17) дает возможность использовать для оценок первых неизвестных и предельного значения метод улучшенной редукции.

В качестве примера рассмотрим определение критических нагрузок и форму потери устойчивости для квадратной пластины. Для этой задачи известна

ЛГ а\Т ж2 Б

N +~т N. =—г

х Ь2 у 3Ь2

3Ь2 3а2 „ —^ + —^ + 2

(19)

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

wr

= (1 -cosж(х/а-1))(1 -cosж(y/Ь-1)).

При изменении Q е[0;2.5] (Р = 0) число нерегулярных уравнений в системе (12) остается одинаковым и равным Ык.= 2. В табл. 1 даны оценки методом улучшенной редукции первых неизвестных и предельного значения в системе (17) для 2гт ■

Таблица 1.

Оценки неизвестных в системе (16)

7} 23 2 4 2 5 <1

0.184113 0.899516 0.478141 0.603929

22 23 2 2 2 4 2 2 2 5 <2

5.28392 1.51362 1.89678 1.98505

Зная решения регулярных бесконечных систем (17), по формулам замены (16) вычислялись первые неизвестные в системах (14), предельные значения < дают возможность описать асимптотику решений этих систем при т > Ыь, откуда получаем асимптотику нетривиального решения бесконечной системы на критической нагрузке:

Ст Сг

Хт = ^ Ут = ~ (т ^да) ■ (20)

Рт «т

Таким образом, на критических нагрузках находятся все компоненты нетривиального решения. Критические нагрузки определяются как нуль определителя системы (15). На рис.1 можно увидеть, что Qc =2.5185(п/Ь)2 (при Рс = 0). В рассматриваемом примере порядок данного определителя Ыя.=2.

Для данной задачи значение критической нагрузки, полученное на основе (19), дается [1] как 0О = 2.667(ж/Ь)2, там же приводится уточненное значение 0 = 2.517(ж/Ь)2 ■

Сравнение с при

Рисунок 1. Определитель системы (15) как функция от 2

В следующей таблице дается сравнение найденных значений критической силы Qс согласно описанному в статье подходу и значений критической силы Qo согласно приближенной формуле (19)

Таблица 2.

:нным решением

Р / 0 0 0.25 0.5 0.75 1.0

Ос (Ь /ж)2 2.518 2.082 1.758 1.514 1.326

00 (Ь / ж)2 2.666 2.133 1.777 1.524 1.333

5.9 2.5 1.1 0.7 0.5

Заметим, что уже при отношении а/Ь > 1.1 формула (19) дает большую погрешность, и для дости-

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

Для определения формы потери устойчивости остается с учетом найденных значений критических нагрузок Pc и Qc выполнить численные оценки нетривиального ограниченного решения квазирегулярной бесконечной системы и воспользоваться аналитическим представлением прогиба в форме бесконечного ряда (3):

■=z

a(-1TxnPn

1(1+sln )(q2n - tin) n+1 . . (

^ch ъ,пх ch ch q1 na ch q2 na

cos РпУ +

+y b (-1)n+1 ynan

tí(1 + )(P2n - Pin )

chPi, „У ch Pi,ny

ch P1,nb ch Pi,nb

cosa x

На рис. 2 представлены формы потери устойчивости согласно (3) для случая одноосного сжатия пластины и для случая двухосного сжатия пластины равными силами.

a) b / a = 1/ 5,(b /n)2 Q = 1.779, Pc = 0

с) b / a = 1 / 2, (b / n)2 Q = 1.965, Pc = 0

b) b / a = 1 / 5, (b / n)2 Q = 0.986, Pc = Q

d) b / a = 1/2, (b / n)2 Q = 0.981, Pc = Q

Рисунок 2. Формы потери устойчивости

Как в случае одноосного растяжения, так и в случае двухосного сжатия форма потери устойчивости имеет тем большее число выпучин, чем сильнее

отличается форма пластины от квадрата. В случае одноосного сжатия это выражается более ярко с увеличением нагруженных сторон прямоугольника.

Список литературы:

1. Биргер И.А., Пановко Я.Г. (1968) Прочность, устойчивость, колебания в 3 т./М: изд-во Машиностроение, 1968. - 3т. - 567 с.

2. Канторович Л.В., Крылов В.И. (1962) Приближенные методы высшего анализа. М. -Л.: Физматгиз, 5-е изд. - 708 с.

3. Papkov S.O., Baneijee J.R. A new method for free vibration and buckling analysis of rectangular orthotopic plates. /Journal of Sound and Vibration 339, 2015. - P. 342 - 359