Об одном решении Навье для частично нагруженной прямоугольной пластины Текст научной статьи по специальности «Математика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2009 Математика и механика № 1(5)

УДК 517.52:539.3

С.П. Сейранян

ОБ ОДНОМ РЕШЕНИИ НАВЬЕ ДЛЯ ЧАСТИЧНО НАГРУЖЕННОЙ ПРЯМОУГОЛЬНОЙ ПЛАСТИНЫ

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

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

Исследуемое решение получено Навье представлением неизвестного прогиба повторными синусными тригонометрическими рядами (ТР) и разложением нагрузки в повторный ряд Фурье. Путем подстановки таковых в уравнение равновесия Софи Жермен - Лагранжа [1] и формального повторного почленного дифференцирования по обоим символам суммирования ТР прогиба при применении бигармонического оператора в уравнении равновесия Навье приходит к уравнению, содержащему четыре повторных ТР. Далее Навье приводит сумму названных повторных ТР к единому повторному ТР с членом в виде алгебраической суммы членов предыдущих рядов и приравниванием коэффициентов полученного ряда нулю конечным числом корректных алгебраических действий определяет коэффициенты ТР прогиба.

В дальнейшем, однако, предложенное Навье решение требует строгого математического обоснования всех формально проведенных при его получении математических операций. Это операции четырехкратного почленного дифференцирования по обоим символам суммирования и операция обьединения четырех полученных таким образом повторных ТР в единый повторный ТР. Отмеченные действия как Навье, так и последующие поколения математиков и механиков средствами из классической теории функций оставляли необоснованными. Впервые попытка их обоснования в рамках классической теории функций принадлежит автору [2], где, однако, использовался искусственный подход.

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

ленной сходимостью рядов решений. Значительный вклад в построение решений для локально-нагруженных тонкостенных элементов конструкций в форме, доступной для эффективных вычислений, внесли А.Н. Крылов, С.П. Тимошенко, В.З. Власов, В.В. Новожилов, В.В. Васильев, В.М. Даревский, И.Ф. Образцов, Б.В. Нерубайло, В.П. Ольшанский, Э.М. Григолюк, В.М. Толкачев, Б.К. Михайлов, С. Лукасевич и др.

Обсуждая вычислительную способность решения Навье, предварительно отметим, что в [1] при сосредоточенной нагрузке относительно решения для данной пластины сообщается, что «двойные ряды этого решения непригодны для получения численных результатов, в особенности, если в них входят производные высших порядков от прогиба». В работе же [3], вообще говоря, отмечается, что «в случае сосредоточенных нагрузок приходится удерживать большое число членов, что повышает трудоемкость вычислений, не обеспечивая достаточной точности результатов. Далее, автор цитируемой книги, в которой обобщены известные результаты в области пластин и оболочек при локаьных нагрузках, приходит к выводу: «Поэтому в случае сосредоточенных нагрузок решений в виде рядов следует избегать. Здесь желательны решения в замкнутом виде».

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

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

(а) Известное неравенство [7]

ар р- + Ъсс[х > аЪ, а,Ъ > 0, р,q> 1, + с- = 1,

введением новых независимых переменных а, в и параметра у зависимостями

с применением предельных переходов при у ^ 0, у ^ 2 представляется в обобщенном виде

1. Некоторые вспомогательные математические результаты

а =

р ра2 , Ъ = ^с р2 , р = 2 у1 , а > 0, р > 0,

а2 +р2 > Суаур2-у , а> 0, р> 0,

= Г 2 у-0,5у(2-у)0,5у-1, при 0 < у< 2,

(1.1)

1, когда у=0, 2.

(Ь) Теорема. Пусть на отрезке А с (- да, да) \ {±2кп, к = 0, 1,...} задан косинусный ТР

ад

^ ат 008 тх, х еД , (1.2)

т=1

с коэффициентами, удовлетворяющими условию

Иш ат = 0. (1.3)

т^ад

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

Если же в принятых предположениях выполняется также условие

Иш тат = 0, (1.4)

т^ад

то ТР (1.2) почленно дифференцируется.

Доказательство. Заметим прежде, что исходный ТР удовлетворяет необходимым и достаточным условиям сходимости ТР Салема [8], откуда следует его сходимость на А. Применяя к ТР (1.2) преобразование Абеля и дифференцируя полученное выражение как дробь, при этом пользуясь, по условию, почленной дифференцируемостью ТР в числителе, получаем

Л ад Л 1 ^ , п сч 1

ат С°5 тх = — ^ Л 5Ш( У-0,5) х = ■

dxm=i m dx 2 sin0,5x “ v ’ 4 sin2 0,5x

ад ад

sin0,5xZАv (2v- 1)cos(v-0,5)x- cos0,5x Z Аv sin( v - 0,5) x

v=1 v=1

(1.5)

где А1 = -a1, Аv = av-1 -av , v = 2,3,... (1.6)

Отметим, что в (1.5) ТР в числителе сходятся, а заменатель - отличен от нуля. Значит, производная на А существует и конечна.

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

d ад 1 ад

— Z am cos mx =-------- ------Z А^, [ (v -1) sin vx - v sin( v -1) x]. (1.7)

dxm=j 4 sin 0.5 x v=j

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

d ^ N [ (v -1) sin vx-v sin( v -1) x]

— Z am cos mx = lim Z А.,-^--------------------------------------

dxm=1 m N^ад v=1 v (4 sin2 0,5 x)

1

lim

4 sin2 0,5 xN-

Z Аv (v -1) sin vx-Z Аvv sin( v -1);

.v=1 v=1

-----1----- lim

4 sin2 0.5xN^w

£ Av(v -1) sinvx -£ Л+1 (v +1) si

Sin vx

= lim

N ^w

N-1

АN (N-1) sinNx-2А2 sinx + £ [^v (v-1)-Ах,+1 (v +1)]

v=2

sin v x. (1.8)

Заметим теперь, что по условию почленно продифференцированный преобразованный ряд сходится на А. Поэтому в силу леммы Кантора - Гейне [9]

откуда

lim АN(N-0,5) = 0,

N ^w

lim АN(N-1) = lim [АN(N-0,5) (1 -0,5/(N-0,5))] = 0.

N

(1.9)

Следовательно, переходя в (1.8) к сумме пределов с учетом (1.9) и (1.6), приходим к доказательству первой части теоремы

d_

dx

£

m=1

1

am cos mx =

(2 a1 - a2 ) sin x - ^ [ (v -1) av-1 - 2 v av +

v=2

4sin2 0,5 x + (v+1) av+1 +(v+1) av+1 ] sin vx

(1.10)

Возвращаясь ко второй части теоремы, заметим, что при выполнении дополнительного условия (1.4) формально почленно продифференцированный исходный ТР в (1.10) оказывается связанным с ТР справа тождественным повторным преобразованием Абеля, что следует из необходимых и достаточных условий сходимости ТР Салема [8]. Теорема доказана.

2. О трехкратной почленной дифференцируемости решения Навье

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

Названное решение в прямоугольнике пластины G =[0, a] х [0, b] записывается в виде [1]

w w

w (^ У) = ££ wmn sinXmx sin V■nУ, m=1 n=1

где Wmn = B0 sinXm^ sin^пП sinXm Д^/2 sin ^n ДП/2/(Xm^n (X]« + ^П)2),

B0 = 16p (Dab Д^Дп), Xm = дал/ a,

|an = nn/b, x, y e G. (2.1)

С применением неравенства (1.1) имеем

I Wmn d S sin Xmxld xS d J sin Цny/d yJ I < B0Xm-1 VJn-' / ( X2 + Цn )2 <

< b^( x г1- цГ2у- j)=bJ( x m^s-j-r),

0 <y< 2, r = 2 y + 1 - s, - 2 < r < 5, 0 < s + j < 3. (2.2)

Здесь B1 > 0 - некоторая константа.

Оценка (2.2) при r = 1,5 приводит к неравенствам

\wmn д s sin х тх/ д Xs д1 sin ц„у/5 у1 \ <Vm\wmn д1 sin ц пу/д у1 |<

< в о хт-ч/ (х т+ц2)2 < Bi/(хт5цП'5_ s-1), 0 < s+j < 3.

С использованием (2.3) и (2.4) приходим к утверждениям

д1

-SSWmn SinХmX SinЦпУ =

(2.3)

да

0гкуда Wmn д S Sin Xmx/д xS д 1 Sin ЦпУ/д У1

п=1

да

<Ц Wmn д S Sin Кх/д xS д 1 Sin ЦпУ/д У]\ <

n=1

да да

<xm si wmnд 1 Sin цпу/5у1 | <вУ xm5 s ^ ц»5^ -1 <да; (2.4)

да да

S 1/x'm5 S 1 / цП'5 ~s -1 < да ,0 < s+1 < 3. (2.5)

n=1

I (2.4) приходим к утверждения!

да

„ , S wmn д SinXmXIд xS Sin ЦпУ =

д у n=i

да

= S wmn д Sin XmXlд xS д1 Sin ЦпУ/д У1 C C [G] (2.6)

n=1

дs+]w (x, у) = д s+1

д xs д у1 д xs д у1 m=1 П=

да да

= SSWmn 5S Sin Х mxlд xS д1 Sin ЦпУ/д У1 С С [G] , 0 < S + 1 < 3. (2.7)

m=1 п=1

Действительно, ряды

да

S Wmn д 1 Sin ЦпУ/д У1, 0 < 1 < 3, У е[0, b] (2.8)

П=1

имеют своими общими членами непрерывные функции переменной у е [0, b] и, что следует из оценок (2.3) и (2.4) при s = 0, мажорируются числовыми сходящимися рядами. Поэтому ряды (2.8) сходятся равномерно к непрерывным на [0, b] функциям [10] (признак Вейерштрасса, п0 430). Таким образом, при У е [0, b] имеем, что в (2.8) при 1 = 0 ряд сходится, а при 1 = 1 - сходится равномерно. Отсюда, с применением теоремы 7 ([10], п0 435) приходим к утверждению (2.6) при j = 1. Отправляясь далее от равенства (2.6) при 1 = 1 и опираясь на сходимость ряда (2.8) при 1 = 1, а также равномерную сходимость на [0, b] ряда в (2.8) при 1 = 2, приходим с использованием названной теоремы 7 к утверждению (2.6) и при 1 = 2. Если теперь при 1 = 2 воспользоваться равенством (2.6) и сходимостью ряда (2.8), а при 1 = 3 - равномерной сходимостью на [0, b] ряда (2.8), то с применением отмеченной теоремы 7 придем к завершению доказательства утверждения (2.6).

Рассмотрим теперь повторные ряды

да да

SSWmn 5S SinXmx/дxS д 1 Sin цп у/ду1, xу е G, 0 < s +1 < 3, (2.9)

как ординарные, в которых производится суммировании по индексу т. Общие члены данных рядов непрерывны в О как произведения непрерывных функций, а сами ряды мажорируются числовыми сходящимися рядами, что следует из (2.6), а также (2.4) и (2.5). Поэтому считая переменную х постоянной, полагая в (2.7) 5 = 0, у = 0 и применяя те же рассуждения, что и при доказательстве (2.6), придем с учетом равенства (2.6) и к утверждению (2.7) при 5 = 0. Пусть теперь фиксирована переменная у. Отправляясь в (2.6) от значений 5 = 0, у = 0, но при х е [0, Ь] повторяя рассуждения, примененные при доказательстве (2.6), приходим к утверждению (2.7) при у = 0. Аналогично, отправляясь от значения 5 = 0 в (2.7) при каждом у = 1, 2, приходим к доказательству (2.7) и при у = 1, 2 и прочих значениях I. Утверждение (2.7) доказано полностью.

3. О почленной дифференцируемости ряда решения Навье при вычислении частных производных до четвертого порядка в уравнении равновесия Софи Жермен - Лагранжа

С применением теоремы п.1 исследуется повторный ряд решения Навье для прогиба на частную почленную непрерывную дифференцируемость под обоими символами суммирования при вычислении производных в бигармоническом операторе в уравнении равновесия Софи Жермен - Лагранжа.

Введем обозначения

42 = В0 /(Кт ^п (+ ^)2К Хп(у) = вшц„у япвш Дп/2; (3.1)

да да

^ 5 (, у)=ЁЁ'^тп^ т цПг5 Хп(у) сов т1,5=1,3, (3.2)

т=1 п=1

и исследуем ряды (3.2) на сходимость. С применением (2.3) при 5 + у = 3 и 5 = 1, 3 получаем

\™тпк ^Г5Хп(у)сов т/| <| < вг/( х1,я5ц1п5); (3.3)

Ё ^Т5 Хп(у) сов т1

да да

•ё\^.пк,^ 1 <^дт5ё^п5,5=1,3. (3.4)

<

п=1 п=1

Таким образом, из (3.3) и (3.4) имеем, что каждый из ординарных рядов в повторном ряде (3.2) мажорируется числовым сходящимся рядом. Тогда оба ряда сходятся. Отсюда, подставляя (2.1) и (3.1) в (2.7) и полагая 5 = 1, 3 с применением формул преобразования произведения синусов, а также синуса и косинуса в их сумму и использованием сходимости в (3.2) ординарных рядов, последовательно получаем

д3 4 (х, у) =

5 х5 5 у 3~5

1 да да

= Т ЁЁ 'кт^Т5Хп (у) [ кт (% + Д%/2) - СОв кт (% - Д%/2)] С^ КтХ =

2 т=1 п=1

1 да да

= 44 ЁЁ'^т'п Кт Нп~5 Хп (у)[С05 Кт (Х + % + Д%/2) + С^ Кт (Х - % -Д%/2) -

т=1 п=1

- сов Кт (х + % - Д%/2) - сов кт (х - % + Д%/2)] =

1 да да

= 7 Z Z ^Г*Хи (У) cos К (х + % + А%/2) +

m=1 1_и=1

да

+Z 42 X m M^* Хи (>’) COS X m (X -%-А%/2) -

И=1

да

-Z 42 К^“*Хи (У) COS Xm (X + % - А%/2) -

И=1

да

-Z 42X *m^r* Хи (У ) COS Xm ( X -% + А%/ 2)

И=1

= 4 ^(-1)!+; ^* (Xi(x + (-1)г % + (-1)1 А%/2), У), * = 1,3. (3.5)

4 i, 1=0

Исследуем теперь ряды (3.5) на частную почленную непрерывную дифференцируемость по x под обоими символами суммирования. Из геометрических условий данной задачи имеем

0 < x < a, 0 <% + А%/2< a , А% > 0. (3.6)

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

X *% + А%/2, (3.7)

что равносильно исключению из областей определения функций Х1(х - £ + Д£/2)

нулей, с учетом неравенств (3.6), (3.7), получаем

X1(х + % + А%/2)е(0,2п), Х1(х-% + А%/2)е [-n,0)U(0,п]. (3.8)

Таким образом, множества значений Z(±), (±) функций X1 (х + £ +Д £ /2 ) переменной х е X(+) = [0, a] \ { х = £ +Д£/2} при каждом выборе знака не содержат точек +2кп, к = 0, 1,... Здесь в L(+),(+) левые знаки соответствуют таковым при £, а правые - при А£. Заметим, что L(+), (+) как множество значений линейной функции х связано с множеством ее определенния X(+) взаимнооднозначным соответствием, причем оба множества в зависимости от выбора знаков состоят одновременно из одного или двух непересекающихся промежутков.

Далее из (3.4) следует, что

да

lim Z^К^Т*Хи(y)cos ml = 0 . (3.9)

т^да

и=1

Таким образом, при ограничении (3.7) к ординарным при суммировании по m ТР в четырех повторных рядах в (3.5) допустимо применение преобразование Абеля [8].

Пусть теперь 1еА, где А - произвольный замкнутый отрезок, принадлежащий одному, но любому из множеств L(+), (+). Применяя на А преобразование Абеля, получаем

да да 1

Z Z42xm^*Хи(y)cos mi=-2. 0 5l m= и= 2sin0,51

X* sin ll2 Z Хи (y)

и=1

да / да да \

+ZIZ * Хи (y)-Z 4-1, их m,-^-* Хи (у) I sin(m - 0,5))

да да да

s = 1,3,

(3.10)

где

(3.11)

Здесь коэффициенты (3.11) преобразуются к виду

(3.12)

(3.13)

Далее с применением неравенств (1.1) и (2.2) для величин (3.12) и (3.13) приходим к оценкам

где В2 , В3 > 0 - некоторые константы.

Таким образом, при 5 = 1, 3 каждый ординарный ряд в повторном в числителе

(3.10), где производится суммирование по индексу п, имеет своим общим членом непрерывные функции переменной у и, с учетом (3.14) и (3.15), мажорируется числовым сходящимся рядом на отрезке [0, Ь] . Тогда каждый из них сходится равномерно [11] ( признак Вейерштрасса, п0 430 ) к непрерывной на [0, Ь] функции переменной у. Но при 5 = 1, 3 и каждый ординарный ряд в том же повторном ряде, образованный суммиованием по индексу т, а также ряд, полученный его формальным почленным частным дифференцированием по I, составлены суммированием по т непрерывных как произведение двух непрерывных функций функций переменной I е Д и переменной у е [0, Ь] функций, а также с учетом (3.14),

(3.15) мажорируется числовым сходящимся рядом. Поэтому и названные ряды сходятся равномерно к непрерывным на Д х [0, Ь] функциям I и у [10] (признак Вейерштрасса, п0 430 ) и, в силу теоремы 7 ( [ 11], п0 435 ), имеет место равенство

Очевидно и произведение в (3.10) sin / / 2 на ординарный ТР - непрерывную функцию переменной у непрерывно дифференцируется по /. Тогда числитель в

(3.10), рассматриваемый как ординарный ряд в повторном, образованный суммированием по индексу т = 1, 2,.. ., почленно непрерывно дифференцируется по /. Кроме того, с использованием (3.4) приходим к неравенствам

т=2 п=1 т=2 п=1

т=2 п=1

sin(m -0,5) /

д/

с C [Дх[0, b]], s = 1,3.

(3.16)

-Ё ^тП 1К ЦП^S Хп (у)sin ml

п=1

<( X MV+V3-s\ І/x, <

— I тп т г'п | / 1 _

п=1

Откуда

lim

да

< Bl/(Xl X0m5 ) 1 /ц)’5 , s = 1,3.

n =1

да

-Xwnl>)mXmц))^хп(у)sin ті |=o, s=1,3.

n=1

(3.17)

Таким образом, ординарный ряд в повторном ряду в (3.10), где т = 1, удовлетворяет всем условиям теоремы в п.1. Следовательно,

д^ (l,у) д і

= X X wilnX m ц))s х п (у)

д cos ml ді

: C[Ax[0,b]], s = 1,3. (3.1S)

Здесь непрерывная дифференцируемость функции ^s (l, у) вытекает из прежде доказанной непрерывной дифференцируемости числителя и непрерывности производной знаменателя алгебраической дроби (3.10).

Утверждение (3.18) доказано для значений ( l, у ) е Д х [0, b] при любом Д с Z(±), (±) и каждом из множеств L(±), (±). Покажем теперь, что оно остается верным и при ( l, у ) е L(±), (±) х [0, b] . Действительно, пусть l0 - произвольная точка из заранее выбранного одного, но любого из множеств L(±), (±). Каждое из L(±), (±) в зависимости от выбора знаков состоит из одного или двух непересекающихся промежутков. Поэтому l0 принадлежит промежутку. Но тогда найдется отрезок До, такой, что l0 е До, и одновременно Д0 содержится в том же промежутке, что и l0. Заметим, что на множестве Д0 х [0, b] утверждение (3.18) по доказанному верно. Значит оно верно и при ( l0 , у) е Д0 х [0, b]. Но ( l0 , у) - произвольная точка из одного, но любого заранее выбранного множества из множеств L(±), (±) х [0, b]. Следовательно, оно остается верным и в нем. Так как выбранное множество -любое, то (3.18) верно и в каждом из множест в L(±), (±) х [0, b]. Что и требовалось доказать.

Докажем теперь, что в повторном ряду в (3.18) на любом из множест в L(±), (±) х [0, b] допустимо изменение порядка суммирования. Так как множество L(±), (±) не содержит точек ± 2 кп, к = 0, 1, ..., и одновременно выполняется равенство (3.17), то к ординарному ТР с индексом суммирования m в повторном ряду (3.18) допустимо применение преобразования Абеля [9]. Применяя данное преобразование, получаем

m. s 3-s . . д cos ml X WmnKm^» Xn (у)-

д і

1

2 X, sin l/2

CJU CO ІЛ)

X wi(n)Xi+Vn“sXn(У)cos V2 + X Xв1п+1)ЦпХп(У)cos (m- V2)l

s = 1,3 .

(3.19)

Здесь

В^+1)ц = (Xs+1 w(m)-Xs+1, w(1), ) ц3^ =

mn r“n v m mn m-1 m-1,n' ^n

[(Xs+1w(')-X 1XX w(')) + (X 1X^w(1)-X^+'w(1), )] Ц3

LV^m mn m-1 m mn' v m-1 m mn m-1 m-1,n'J r^n

= X,Xs w(1) ц3~s +X ,B(s)

1 m mnr^n m-1 mn'

(3.20)

Далее с использованием (3.3), (3.14) и (3.15) оценим величину общего члена повторного ряда в (3.19):

Значит, замена общего члена в повторном ряду в (3.19) их абсолютными вели-

в повторном ряду в (3.19) справа допустима перестановка порядка суммирования. Кроме того, с учетом (3.3) имеем

что позволяет применить преобразование Абеля [8] к ординарному ТР с индексом суммирования т = 1, 2, ... в повторном ряду (3.19), если в последнем изменить порядок суммирования на обратный и обьединить оба сходящихся ряда с индексом суммирования п в единый сходящийся ряд. Выполняя данные преобразования, получаем

Далее равенство (3.5) дифференцируется по х е Х(±) , причем в правой части, в силу существования производных каждого из четырех рядов, переходим от дифференцирования суммы рядов к сумме производных от рядов. Диффенцируя теперь каждый ряд как сложную функцию переменной I , соответственно полагая I = (х ± ± ^ ± А ^ /2 ), используем равенство (3.18) на £(±), (±) х [0, Ь]. Затем последовательно вводим символ конечного суммирования по индексам / и] под символы бесконечных сумм сначала с индексом т, а затем и п. Данные действия допустимы ввиду сходимости каждого из ординарных бесконечных рядов как рядов, образующих сходящийся повторный ряд. В итоге имеем

(3.21)

чинами приводит к сходящемуся ряду. Поэтому в силу теоремы 3 ( [10] , п0 393 )

откуда

(3.22)

5 = 1,3.

(3.23)

т=1 п=1

да да

= XX н-™ С С [ G \{х = / ±А^/2}], * = = 3. (3.24)

m=1 n=i 3 х 5 у

Здесь мы использовали (3.18). Если же воспользоваться (3.23), то аналогично получаем

д 4 1 (Х у) XX ... 3s+1 sin X пх 53-* sin vmy _ С G \ у e.A/ }

= X mii ’""'^T*1-----------------------------------------------------------------dT*- c С [G \{x = =2}

* = 1,3. (3.25)

Заметим теперь, что и в исходном повторном ряду (2.1) допустима перестановка порядка суммирования, ибо ряд, составленный из модулей его общего члена, сходится, что следует из (2.3) при s = j = 0. Выполняя названную перестановку и повторяя рассуждения с точностью до обозначений, вместо (3.25) будем иметь

д 4|1'( х,у 1 =XX н д'" sinX"X d3'-* sn»ny С [ G \ {у = г, + Дч/2}] а х3-* д у**1 m=1 n=1 тп д у**1 а у3-* [ у v }

* = 1,3. (3.26)

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

Таким образом, согласно (3.24) и (3.25), все частные производные, входящие в уравнение равновесия Софи Жермен - Лагрангжа, рассматриваемые в совокупности, непрерывны на множестве

E = (G \ {х = / +А//2 })n(G \X у = п + Ап/2}) =

= G\({х = / + А//2}и{у = л + Ал/2}) (3.27)

и вычиляются почленным дифференцированием решения Навье (2.1) для прогиба под обоими символями суммирования.

Подставляя теперь повторные ряды (3.25) и (3.26) в уравнение Софи Жермен

- Лагранжа и используя в них, а также в повторном ряду для нагрузки сходимость ординарных рядов с индексом суммирования m и п, объединяем прежде всего, слева в уравнении, все четыре ряда в единый сходящийся ординарный ряд с индексом суммирования m. Данный ряд содержит в качестве общего члена сумму четырех прежних сходящихся ординарных ряда с индексом п. Поэтому, объединяя под индексом суммирования п и названные ряды в единый сходящийся ряд,

приходим на множестве E к уравнению, содержащему в левой части повторный ряд с нулевыми коэффициентами, а справа - нуль, что обращает уравнение в тождество.

Таким образом, решение Навье с применением разработанных нами средств методического характера из классической теории функций обосновано.

4. Об улучшении вычислительной способности решения Навье

Здесь автор настоящей работы, как и прежде в случае сосредоточенной нагрузки [3] - [5], ограничиваясь лишь повторными ТР, утверждает, что названные ряды - чаще уже полученные решения - могут стать весьма эффективным средством вычислений их сумм, если их предварительно преобразовать по предложенной им методике в быстросходящиеся ряды. В названной методике распространяется высокоэффективный метод А. Даду [12, 13] (п крат примененное преобразование Абеля), предназначенный для ускорения сходимости ординарных ТР,

на повторные ТР. В итоге полученное преобразование - точное аналитическое преобразование, которое сохраняет силу при произвольной нагрузке, в частности локально-распределенной. Его же вычислительная эффективность обосновывается как аналитическими оценками коэффициентов полученных рядов, так и вариантными вычислениями при различных значениях n с взаимным сравнением результатов.

С использованием известных формул для перерезывающих сил [1] и (2.7), (1.1), имеем

f ^ 1 D да да

Nx = -D —w= —---------XXdmn Xn (y)cos КX sinX„^ sinXm/2 =

d x d x d y

где

4 P = - ЬЛ.Л X (-1)г+; V(Xi(x + (-1)! | + (-iy A^2,y)), ab A^An t, j=o

dmn = 1I( Hn (Xm + нП )), 0 < dmn < B ) ] y=2-6 =

= в4/(x2my;e), 0<0<2,

да да

v (I, y) = dmn In

(y) cos ml.

m=1 n=1

(4.1)

(4.2)

(4.3)

Здесь В4 - некоторая константа.

Выражение (4.1) получено с использованием сходимости рядов (2.7), а также ряда (4.3), которая устанавливается с использованием оценки (4.2).

Далее преобразуем (4.3) с применением одно- и двукратного преобразования Абеля. Названные преобразования допустимы, так как

X dmnXn (y) c0s ml

n=1

<XI dmn Xn (y) c0s ml \ <XI dm

n=1

n=1

Откуда

Поэтому

< В4/x^ X^M,5 < да.

n =1

да

lim X dmnXn (y) c0s ml = °-

m^да

n =1

да да

V (l, y) = XX dmnXn (y) c0s m l =

(4.4)

m=1 n=1

да -j да да i да

x d1n Xn (y) +—xx ^ (y) sin (m - v2 )l =- - X dm Xn] ]+

n=1

sin l/ 2

m=2 n=1

' n=1

4sin21/ 2

x d 2n Xn (y) cos l+x x dmn)Xn (y) cos ml

(4.5)

где

dmn= dm-1,n dmn =X1(2 Xm / ( Hn (Xm-1 + M-и ) (Xm + Hn ) ), (46)

дШп - дт+\,п + 2 дтп Ат_1п - 2 Х1 (ци 3 Xт +Х])/

/(Цп (X т-1+ц 2) (X т+ц 2)(х т+1+ц2». (4.7)

Здесь применение неравенства (1.1) приводит к оценкам

0 < ^тп - о(хт/(Цп (Xт+цп))) < в^(ху)] -

- в^(Х3т5цп+5); (4.8)

1 4^ О ( Xт/ ( Цп (Xт + Ц^ )) < Вб /(XШУп-2У)] у-(4-5)/2 -

- в7(XШ-5Ц1n+5),0 <5< 4. (4.9)

Сравнение оценок (4.2), (4.8) и (4.9) показывает, что при первоначальном при-меннении преобразования Абеля улучшается сходимость каждого из ординарных рядов в повторном ряду в (4.5), при повторном же применении названного преобразования улучшается сходимость лишь ординарного ряда в повторном с индексом суммирования т.

Аналогичные результаты в силу симметрии выражения (2.1) получаются и для Ыу перестановкой в (4.1) - (4.9) символов х и у, а и Ь, £ и п, Д£ и Дп.

Далее приводятся численные результаты вычислений безразмерных перерезывающих сил Йх = Ых а / Р, Йу= Ыу а / Р. Расчеты производились при локальной нагрузке как без ускорения сходимости по формуле (4.1), так и с ускорением сходимости с применением формул (4.5) при неизменных значениях следующих параметров:

х/а - 0,49; у/Ь -|/а - п/Ь -0,5; Д|/а -Ап/а -0,005. (4.10)

В таблице при различных принятых значениях Ь/а содержаться результаты расчета безразмерной перерезывающей силы Их.

Ь/а Кх Йх / Йхт

I II III IV I II IV

1,0 17,6 15,6 15,8 15,9 1,11 0,99 1,01

1,1 17,6 15,6 15,8 15,9 1,11 0,99 1,01

1,2 17,5 15,6 15,8 15,8 1,11 0,99 1,00

1,4 17,4 15,6 15,8 15,8 1,10 0,99 1,00

1,6 17,3 15,6 15,8 15,8 1,09 0,99 1,00

1,8 17,1 15,6 15,8 15,8 1,08 0,99 1,00

2,0 17,0 15,6 15,7 15,7 1,08 0,99 1,00

3,0 16,1 15,1 15,2 15,2 1,06 0,99 1,00

5,0 14,0 13,4 13,5 13,5 1,04 0,99 1,00

В левой части таблицы в первом столбце приводятся значения, полученные без ускорения сходимости, во втором же и в третьем - с ускорением сходимости одно- и двукратным преобразованиями Абеля, а в четвертом - расчет при сосредоточенной нагрузке, полученный в [6] с применением трехкратного преобразования Абеля. Для сравнения результатов в правой части таблицы приводится отношение соответственных численных значений в столбцах I , II, IV к таковым, приведенным в столбце III. Сравнительный анализ показывает, что расчет без ускорения сходимости может привести к погрешностям, превосходящим точность классической теории пластин: относительная ошибка достигает 11 %. Примене-

ние же одно - и двукратного преобразования Абеля приводит к значениям, близким к таковым, полученным при сосредоточенной нагрузке: здесь разность относительных ошибок не превосходит 2 %.

Расчет же Ny при тех же значениях параметров во всех вышеупомянутых вариантах вычислений, связанных как с ускорением, так и без ускорения сходимости, показал, что порядок их величины не превышает 10-14, что обьясняется симметричным расположением нагрузки относительно оси OX и вычислением значений Ny в плоскости XOZ.

ЛИТЕРАТУРА

1. Тимошенко С.П., Войновский-Кригер С. Пластинки и оболочки: Пер. с англ. М.: Физ-матгиз, 1963. 635 с.

2. Сейранян С.П. Об обосновании одного решения Навье // Актуальные проблемы механики сплошной среды: Труды Международной конференции, посвященной 95-летию акад. НАН РА Н.Х. Арутюняна., г. Цахкадзор (Армения), 25 - 28 сентября, 2007. С. 391

- 395.

3. Лукасевич С. Локальные нагрузки в пластинах и оболочках: Пер. с англ. и польск. М.: Мир, 1982. 542 с.

4. Сейранян С.П. К решению В.З. Власова задачи изгиба прямоугольных в плане момент-ных пологих оболочек поперечной силой // Композиционные материалы и оптимальное проектирование: Тез. докл. Международной конференции. г. Агавнадзор, 25 - 28 сентября. Ереван: Гитутюн, 2006. С. 56 - 57.

5. Сейранян С.П. К задаче об изгибе прямоугольной пластины поперечной силой // Проблемы динамики взаимодействия деформируемых сред: V Международная конференция, г. Горис, 1 - 5 октября. Ереван: Гитутюн, 2005. С. 314 - 318.

6. Сейранян С.П. Об ускорении сходимости в задаче об изгибе прямоугольной пластины поперечной силой // Избранные вопросы теории упругости, пластичности и ползучести: Сб. статей, посвященный 75-летию акад. НАН РА М.А. Задояна. Ереван: Гитутюн, 2006. С. 266 - 273.

7. Беккенбах Э., Беллман Р. Введение в неравенства. М.: Мир, 1965. 165 с.

8. Salem R. Essair sur les sёries trigonomёtriques // Actual. Sci. Industr. № 862. Paris, 1940.

9. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. М.: Наука, 1970. Т. 3. 656 с.

10. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. М.: Физмат-гиз, 1959. Т. 2. 807 с.

11. Dadu A. Une mёthode d’accёrlёration de la convergence des sёries trigonomёtriques // Mathematica - Revue d’Analyse Numёrique et de Throne de rApproximation. 1980. T. 9. № 1. P. 27 - 33.

12. James E. Kiefer, George H. Weiss. A comparison of two methods for accelerating the convergence of Fourier series // Comput. Math. Appls. 1981. V. 7. No. 6. P. 527 - 535.

СВЕДЕНИЯ ОБ АВТОРЕ:

СЕЙРАНЯН Сурен Паруйрович - кандидат физико-математических наук, старший научный сотрудник Института механики Национальной академии наук Республики Армения. E-mail: [email protected]

Статья принята в печать 19.11.2008 г.