Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Амиро И.Я. Ребристые цилиндрические оболочки

.pdf
Скачиваний:
14
Добавлен:
24.10.2023
Размер:
11.82 Mб
Скачать

напряженное состояние оболочки практически не зависит от Bsnt (0).

Формулы (11.118) можно получить и непосредственно из системы уравнений (II. 1), если последнюю упростить в соответствии с та­

кими

предположениями:

 

1.

Перемещения можно представить в виде

и(£,Ѳ) = «0ß , Ѳ) +

+ ц, ß, Ѳ); Vß, Ѳ) = ü0 ß , Ѳ) + ü, ß , 0); w ß, 0) =

ш0 ß, Ѳ )+ш Д, Ѳ).

2.Изменяемость uQ, v0, w0 в окружном направлении такая же, как изменяемость внешней нагрузки.

3.Деформации срединной поверхности оболочки, определяемые перемещениями и0, v0, w0, значительно больше соответствующих де­

формаций, определяемых перемещениями ult ѵъ wv

4. Касательные усилия взаимодействия отшивки и ребер распреде­ лены по ширине панели между ребрами равномерно, влияние дис­ кретного размещения ребер полностью определяется радиальными усилиями взаимодействия.

d2wx

б2^

5 . at «

6Ѳ2

В соответствии с этими предположениями можно упростить и

формулы для

вычисления

усилий

в

обшивке

и

ребрах

(1.32) и

(1.33). Так, для оболочек, нагруженных

циклически симметрич­

ными нагрузками, после элементарных выкладок получаем

Nx *

 

Eh

du^

 

 

N

 

 

Eh

 

 

 

 

( 1 - v V

 

dl • — ѴШ„

 

 

( 1 - v V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Nху ÄJ N ух ‘0,

 

 

 

 

 

 

 

М .

 

 

Eh 2

dlw.

 

M

=

vM ,

M

 

 

 

 

 

 

12(1

— V2)/"2

i ,

xy

yx'

 

 

 

 

<502

''■ *

 

 

 

' "

 

 

N ,

E ciF ci

duo

 

 

 

E

•/ /

/

d?w

 

d?Wi

;

 

 

 

 

Z E U E L I

— ± - _ j -------------- !

 

 

 

 

 

d-l

M ‘ *

s = s

|

L

 

CL

 

 

 

 

 

' 2

\ d t

 

d t

 

 

K i

= ^

o.

 

Покажем, что предложенный выше приближенный метод позво­ ляет с достаточно высокой точностью определить напряженное со­ стояние ребристой оболочки. Сравнение с точным решением удобно выполнить на примере шарнирно опертой оболочки, нагруженной внешним давлением q. Поскольку точное решение этой задачи по­ лучено выше (§ 1 данной главы) в двойных тригонометрических рядах, приближенное решение также будем разыскивать в форме двойных тригонометрических рядов. Примем, что ус= бс = 0, тіс

а2. С учетом принятых допущений решение системы уравнений (IIЛ), упрощенной в соответствии с предположениями V—5, полу-

6—39

81

ч а е м в т а к о м в и д е *:

Л

і - Л

2v

 

 

 

» - 0 ,

и

Eh

q г 1 + 2 (1 -

 

 

 

4 m=

 

 

 

r* (1 — V*12)

 

 

Л і - ѵ 2)

 

w =

- q — 1+ A ilziZLs®

m

^ -

где

 

 

 

 

 

 

 

 

 

^

^ n - o

n .

 

 

 

Формулы для определения усилий в обшивке записываются так:

Л /,= 0,

N = qr —

22°

= qrT'

 

 

1+ 22° (1 - V2)

Му = - a2kV q -

 

 

blk cos IkB

22° (1 — V2) I

S

- akV qG '

0

1 +■--------------

/=1

 

Mx = vMy = — a k 2r2qG\.

Вычисления выполнены для оболочки2, у

которой с — 0,1258,

т)*= 2,40,

= 9,6,

blk— 1, V = 0,3. Результаты сопоставления

прогибов,

усилий и

моментов, полученных

с помощью точного3

и приближенного решений, приведены в табл.

11. Из анализа дан­

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

Рассмотрим граничные условия. Статические граничные условия после подстановки в них (11.118) и пренебрежения малыми слагае-

 

С О

1Поскольку при 0 < С К і в РЯД

qmsin dmg = 1, прогибы, а также при-

т—1

веденные ниже усилия не зависят от ^ в этом интервале. 2 Для этой оболочки неравенства (11.76) выполняются.

3В качестве «точного» решения использовались данные, полученные для

£= -|^. однако необходимо отметить, что слабая зависимость приведенных па­

раметров от I, полученная из приближенной теории, обнаружена и при табу­ ляции точного решения. Как в точном, так и в приближенном решении сумми­ ровался 41 член ряда но I.

82

Т а б л и ц а 11

 

 

 

 

к Ѳ

- 0

 

 

 

 

 

 

 

 

А’В=л

 

 

Р еш е­

 

*

 

 

 

 

 

*

 

*

 

 

 

»

 

*

 

ние

 

 

Т '2

 

 

 

 

ГІ

 

 

°2

 

W *

 

 

 

 

с і

С2

W *

 

Т 2

 

С І

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Точ­

0,01

0,04

—0,97

0,42

1,43

1,95

—0,05

—0,93

—0,23

—0,80

ное

Приб-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ли-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

жен-

0

0

—0,94

0,47

1,55

1,94

0

—0,94

—0,24

- 0 ,7 8

ное

мыми приобретают такой вид:

 

 

 

 

 

 

 

 

 

 

 

8—46'

0X„lSl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о:,

sin n ft) +

 

 

( 1 - л

 

 

 

 

 

{б»,п, к , «>* " ,в +

 

 

 

s = l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

Д

Д

' (bn COS пѲ) +

(bn sin пѲ)]} = Г 1Яі>

 

 

 

Eh

 

8 - 4 6 '

 

 

 

 

 

 

 

 

 

 

 

 

-

2

/ * ', t [ s , s ; , ( - c : |Sin»,e

+

 

 

 

 

( 1 - л

 

 

 

 

 

 

 

 

S=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

Д

 

cos п,Ѳ) + Snts (0)] =

5,Пі,

 

 

 

 

 

 

 

8—46'

 

 

 

 

 

 

 

 

 

 

 

 

- д

г Ѵ

 

 

I

с

 

<6А ,

(с ».cos “іѳ +

К

sin л.ѳ>+

 

 

 

 

 

S=1

 

 

 

 

 

 

 

 

 

 

 

+ Q ‘„ K ,

Ic ; ® r

№. « »

 

я0) +

о ; ® ?

№. Sin ВѲ )] +

 

і (Ѳ )> =

 

 

 

 

 

8—46'

 

 

 

 

 

 

 

 

 

 

 

 

 

------—~—2~

V

еХп'Л {ЬП* (Cs

cos n .Q + D ns sinn.0) +

 

 

 

] __ѵ 2

 

s=1

 

 

1

л,

п, \ Л,

1

1

Л,

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

Д

[С;ф ^‘ (Ьаcos пѲ) +

D \Фп2' (bn sin лѲ)] +

G^ (Ѳ)} =

Glv

Здесь

приняты

обозначения1:

 

 

 

 

 

 

 

 

Т п , =

Д , -

 

Ѵ П 1М п , -

v F n t -

K

S n , s L n , -

V , h K , -

 

 

n , = Т

с * . , ,

д , - д ? ч > - e j £ , д , - д г § ; ) +

 

 

 

 

 

 

 

 

 

Y Y

 

А

— 6 Y2 C ,

 

 

 

 

 

 

 

 

 

 

 

 

1

I C ^ / i j S

 

 

c ^ « | S / І 1 *

 

 

 

 

 

1 Как нетрудно проверить, 5Ді(Ѳ) мало по сравнению с 7 я>и S„*.

6*

83

S*H, =

 

 

 

+ ^sK) + a\s'hFSn, + a\ K

~

- ь ? [ ^ ( ^ + * п Л ) + а\* (пА

+ Mn)i.

К

(ö) = ■ Ц ^ х - .л (Ѳ ) + ß2^ , s[ d- T

- +

B l

(0)] >

< 4 = Л ц . K

. -

( 2 — V) ^

^

-

( 2 ■ - V)

 

-

 

- ■ K S( B i . -

(2 -

V) « ’ ] Г Пі -

 

(2 -

V) « , / Q } ,

 

 

 

 

 

 

 

 

 

 

dK

w

,

, 0

 

< (Ѳ )

<?; (Ѳ)= *

X ^ ^ , (Ѳ) + (2 — v) x„iS —^ t

+

^2

V^Xn,5

до

 

 

 

 

 

 

 

 

wo8

Qs

=

-n Y3

(^s

— A?‘T

) — 6 Y2

— A?'SZ,S) +

^я,

 

"c^njS *

rtj

 

1

^ r t i '

cЛ/ijS ' n,

 

1

л,/

~

 

 

 

 

+

Т ІсО І.Іц -

ö 'X l X -

 

 

 

 

Gl = «“ « X l -

4 ) П , -

 

v n ^

-

 

[(зі. -

v / i X

-

 

К

=

л .< . к

 

-

к

х

) -

6ex„lS cs; - а г а ;> +

 

 

 

 

+

\

X; A

, - ÖCX„1A 1.

 

 

 

 

с ; (Ѳ) =

а2

' 2

„ s

 

) +

d2C*

(Ѳ)

dB' (Ѳ)

 

x ; sc ; ( 0

v - ^

i i +

v - ^ i -

 

 

Л„

=

4

Af

 

У"і

' б '

1

I O

X ” iS

y « l

 

 

 

nl

 

1 — V

£2

2

а2*4V J +

<Л6

6 ’

 

_ 22"’

■~ Л 4“ '

Кинематические граничные условия можно получить, приравняв перемещения (11.116) и угол поворота к заданным. Очевидно, что трудности в решении задачи не исчерпываются выводом формул (11.116). Как нетрудно заметить, восемь или четыре (при ^ = 0 ; 1) произвольных постоянных в общем случае недостаточно для удо­ влетворения бесконечной системы граничных условий. В связи с этим необходимо ввести дополнительные физически оправданные предположения, позволяющие свести граничные условия к системе линейных алгебраических уравнений восьмого (% > 1) или че­ твертого (при пх= 0; 1) порядка. Рассмотрим два варианта пред­ положений,

84

1. В граничных условиях можно отбросить слагаемые, изменяе­

мость

которых

в окружном

направлении равна Ik ± 'h (і > 1)

(если

оболочка

подвержена

действию циклически симметричной

нагрузки, то этому варианту соответствует интегральное удовле­ творение граничных условий).

2. Граничные условия можно удовлетворять только на ребрах. Качественный анализ погрешностей, вносимых этими предпо­ ложениями, проведен на примере полубесконечной оболочки, на­ груженной циклически симметрично. Приведем выводы, получен­

ные в итоге указанного исследования.

Если на торце оболочки заданы усилия, то более близкие к точным значениям произвольных постоянных дает второй вариант предположений, причем, чем больше жесткость ребер и их число, тем ближе к точным значения произвольных постоянных, найден­ ных из приближенных граничных условий. При весьма большом числе ребер {аг& 1) первый вариант предположений также дает величины произвольных постоянных, достаточно близкие к точным.

Если на торце заданы перемещения, то второй вариант предпо­ ложений также приводит к значениям произвольных постоянных, более близким к точным, чем первый. Первый вариант предполо­ жений приводит к хорошим результатам при достаточно большом числе ребер.

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

Приближенное решение уравнений

равновесия

определялось

по

формуле (II. 118). В качестве

частного решения

неоднородной

си-

стемы уравнений равновесия

(II. 1)

принималось

W04

Для

рассмотренных задач был важен способ удовлетворения граничного условия относительно прогиба (неоднородное граничное условие). Вычисления выполнены для двух оболочек [ус=бс = 0, с = 0,1258; kt,2 = 9,6;ѵ = 0,3; Дс » 0; т)* = 240 (табл. 12) ит]* = 60 (табл. 13)]; усилия и прогибы определялись в среднем сечении оболочки. Ока­ залось, что при принятой величине с более близкие к точным1 значениям прогибов и усилий дает второй вариант предположений, В табл. 12 и 13 приведены данные, полученные на основе этого ва­ рианта предположений. Очевидно, что с увеличением с различие между числовыми данными, полученными для первого и второго вариантов предположений, будет уменьшаться.

Из сравнения данных табл. 11 и 12 следует, что для рассмотрен­ ных оболочек приближенное удовлетворение граничных условий не

1 В качестве точных значений усилий и прогибов для шарнирно-опертой оболочки принимались значения, полученные в результате табулирования двой­ ных тригонометрических рядов, для жесткозащемленной оболочки — одинар­ ных (см. § 3 этой главы).

85

 

 

 

 

 

 

 

 

 

Т а б л и ц а

12

Реше­

 

 

кВ=0

 

 

 

 

К ІІ ©

 

 

 

 

*

*

*

 

 

*

*

*

 

 

ние

ДО*

°2

ДО*

°2

 

г і

Г2

 

г і

Г2

°1

Приб-

 

 

 

 

 

 

 

 

 

 

 

ЛН-

 

 

 

 

 

 

 

 

 

 

 

жен-

0,02

0

—1,03

0,45

1,50

1,91

0

—1,03

—0,23

—0,77

ное

Точ­

0,01

0,04

—0,97

0,42

1,43

1,95

—0,05

—0,93

—0,23

—0,80

ное

 

 

 

 

 

 

 

 

 

Т а б л и ц а

13

Закреп­

 

 

*.ѳ=о

 

 

 

 

*Ѳ=я

 

 

 

 

 

 

 

 

 

 

 

 

ление

до*

Г2

 

ДО*

ті

Г2

С1

С2

торцов

г .

с і

С 2

Шар-

0,10

0

—1,03

0,43

1,42

1,56

0

—1,03

—0,22

—0,71

нир-

0,12

0,04

—0,98

0,38

1,29

1,88

—0,04

—0,93

—0,21

—0,74

ное

 

 

 

 

 

 

 

 

 

 

 

отти­

 

 

 

 

 

 

 

 

 

 

 

ра­

 

 

 

 

 

 

 

 

 

 

 

ние

 

 

 

 

 

 

 

—1,03

 

 

 

Же-

0,02

—1,03

0,46

1,54

1,91

—0,23

—0,77

сткое

0

 

—0,97

0,40

1,30

1,82

 

—0,94

—0,23

—0,76

за­

 

 

 

 

 

 

 

 

 

 

 

щем­

 

 

 

 

 

 

 

 

 

 

 

ление

 

 

 

 

 

 

 

 

 

 

 

П р и м е ч а н и е :

В числителе указаны числа, полученные

в результате табуляции

при­

ближенного решения, в знаменателе—полученные в результате табуляции точного решения.

вносит дополнительных существенных погрешностей в расчет. Сравнение данных, приведенных в табл. 12 и 13, показывает, что уменьшение жесткости ребер в четыре раза не привело к существен­ ному уменьшению точности приближенного расчета. Сравнение результатов точного и приближенного расчетов жесткозащемленной оболочки показывает, что и в этом случае получены достаточно близкие результаты. Однако здесь следует отметить, что точное ре­ шение в данном случае также будет приближенным, а его погреш­ ность такова, что учет отброшенных членов должен привести, на­ пример, к увеличению G* и Gf в точке Ѳ = 0.

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

86

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

Уравнения равновесия ребристой оболочки вращения в перемещениях можно вывести так же, как это сделано в главе 1 для цилиндрической оболочки, в связи с чем вывод уравнений равновесия не приводится.

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

Уравнения равновесия в перемещениях имеют вид

л

Т

Z i ( ц , о , w) +

— Л 4

- < ( Р ) I 4>с ( Р )

L ic (“ >« О d P*= °<

 

1

Л

 

 

 

~k

 

 

 

L2 (и, V, Е0 ) = О,

(И.119)

 

 

л

 

Ц (и, V, W)

АгА%~ ^ ( ß ) j ^ ( ß )

L3c w) dß == °-

л

~k

Здесь приняты такие обозначения: L1( Л2, L3— дифференциальные операторы уравнений равновесия неподкрепленных оболочек вращения в перемещениях [48];

Li >

 

 

 

/ 1 ди

 

w V

 

* д 1

dw ' J.

 

 

 

 

 

 

 

 

 

~

da

A\

 

J

 

 

д 1

 

д 1

д 1

dw

 

d_ 1 d l

du

L 3Ct'U’ W) =

R

\ g- ^

g~

g -

A ^

g~

'

da

da

\

da.

, б

 

 

 

ö

d

1

 

dw

 

Ус

( ди

 

 

c da

\

da

 

 

да

\

 

да

 

R ^ *

\ da

 

 

 

~

E/ r k

g*

 

hc ~

~

 

^yc

W

^

 

 

 

Yc_ 2лR*B

c ~

R* Yc’

 

_

 

2n (R*)3B

 

 

a, ß — ортогональные криволинейные координаты точки на срединной поверхности обшивки; Аи А2, Ri, R2— коэффициенты первой квадратичной формы этой поверхности и ее главные радиусы кривизны: R* (R* > /г) — характерный раз­ мер обшивки;

 

/ * (ß )

при

ß £

-

ße

-

,

ßc '

 

 

 

 

ß i

~2~ ’ ß l

'

2

J’

<

(ß ) =

 

 

ß

ß°

г

ß°

 

 

0

при

ß £

 

 

 

~2~ ’ ß j

 

 

 

 

 

 

 

 

 

 

 

2nA2/(a,ß)

где ß; ~k і, ßcA3

— ширина

ребра (t =

1, 2,

• ,

й). /*№)

2

87

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

действия по ширине ребра, / (a,ß) —

,

где Ф (ß)

— известная функция ß.

 

 

_

_

 

Аг

 

 

 

 

На торцах оболочки а : и а.,

статические граничные условия задаются в виде

 

 

.

 

л

 

 

 

 

 

 

 

 

Т

»

^

 

 

'"w

_

 

 

kBR

_ Г

(Ц'

 

N-a+ ' * -

<P)j

 

<ß )

^ =

Nä>

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

^ a ß

 

S a ß ’

 

 

 

 

,

kB (tf*)2

,,

 

 

дМ-с(и, иі)

 

( 11. 120)

 

 

2яЛ,Лг

®JЛ

■Фсф)

да

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

Л*й +

kV ü 7

< (ß)J

<

(ß)

(«. ®) dß=M_.,

 

 

 

 

Л

 

 

 

 

 

а кинематические

в

виде

 

~~k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и =

и, V = о, ю = к». фх =

Фі.

 

(11.121)

В уравнениях (11.120), (11.121) приняты такие обозначения: Ти Si, Qi, üi — соответственно нормальное и обобщенное касательные усилия, обобщенная пере­ резывающая сила и изгибающий момент на торце неподкрепленной оболочки [48]:

 

ди

w \

бсЯ*

д

1

dw

}

д/

т*

 

 

да

J

 

Л!

да

Ä!

ö«

 

1 1 ’ 2л.Л,

 

Mäc

V ?*

^

'

de

 

j

 

ди

до

 

Лх да Аг да

 

':і

да

 

 

 

 

 

 

 

 

 

_

Q, +

2яЛ2

i|)c(ß) Q0,

 

 

^ = G; + 2 ^ < ( P ) ( GC+ W

 

 

где Т*, S*. Q*, G* — внешние усилия,

заданные на торце обшивки;

и, v, w, ф]—

перемещения, заданные

на

торце

обшивки; Тс, Qc, Ga — усилия,

заданные на

торце ребра.

Предполагаем, что k достаточно велико и изменяемость напряженно-деформи­ рованного состояния в окружном направлении такова, что для расчета обшивки можно использовать уравнения «обобщенного краевого эффекта» [48], т. е. уравне­ ния, в которых отброшены слагаемые, пропорциональные /г2 в первом и втором и /г2«, №ѵ в третьем уравнениях равновесия (11.119). Соответствующие слагаемые отбрасываются и в статических граничных условиях (11.120).

88

Решение системы уравнений равновесия (11.119) будем искать в виде рядов

со

со

 

со

 

и = Л

Н; COS /Äß, V =

sin/&ß, Id =

wl COS lk$.

(11.122)

/=о

;=i

 

i= о

 

В результате подстановки (11.122) в (11.119) получаем однородную бесконечную систему обыкновенных дифференциальных уравнений.

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

0.25

 

0,25 _

,—0,25

 

гі

Ф(а)

Т)

Ф(а)

W, = w,eи с

< р ( а ) '

(11.123)

Здесь ф (а) — функция

изменяемости;

ul, v l,w l — функции интенсивности.

Асимптотические ряды для u[t vf w* строим так:

*

*S *

*c • V~>

= 1

“/ - vi = L "Л wi = \ wi .

s=0

s=0

s=0

Функции интенсивности в s-м приближении представляем в виде

,

0,25—0,5s

*c

—0,5s .

_

—0,5s

Щ =

ul I1! с)

 

'— *e

-*

 

.w /= ^ (T ls)

,Wt

------= ®; (11C

 

(11.124)

(11.125)

где т|с = T]cé4. Подставив (11.123)— (11.125) в уравнения равновесия, собрав и приравняв нулю1 коэффициенты при т)* в одинаковых степенях, получаем после­

довательность бесконечных систем обыкновенных дифференциальных

уравне­

ний, каждая из которых содержит только uf, of, wf(s'

$).Системы,

входящие

в указанную последовательность, можно записать в таком виде:

 

 

4 (“0' °’ w0>+

ЛЙ2

 

 

 

 

-<)= 0,

 

 

 

 

 

 

 

 

/,=o

 

 

 

 

 

 

 

L3 (UQ, 0, aDg)

h2

1

d

A2

d

1

 

d

A2 du)Q

2

 

12 AXA2 da Aj da AXA2 da

da

 

 

 

 

 

 

 

 

 

iR*f

Y

W

s c K

’ t y

- 0’

 

 

 

 

 

AIA2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f,=о

 

 

 

 

 

 

 

 

L\ {u]

\ of

wst

*)

1 — V Pfe2

 

2bthR*

bL .Jus~ \ ws~ ') =0,

 

 

___U/s + —Lk—

X ' b

 

 

 

 

 

 

Ai

1

 

AtA2

2 J

hh lc(

l'

 

 

 

 

 

 

 

 

 

 

 

ii=о

 

 

 

7 '/ „ s

s ...s, I

1 V

1

d

1

 

d

ІЛ ,,s—1\

 

 

 

 

 

 

 

 

 

 

 

 

 

М /Г 1) = c

При

этом

использовались

предположения т|е > ■

Л2

, бе=

:6Ч£'5, |б'К1.

89

_

h21

 

 

 

h2

d

A. d

1

d А dwf 2

 

’) + ü j 4 ( ^ ) 4 +

 

 

 

 

 

 

 

26,. (R*)2

 

 

 

 

 

 

 

(11.126)

где

+ 4 &

г - ь

V

L3cK> <> = °-

 

 

(,=o

 

 

 

 

 

 

 

Lj(<4, vsu m^) = Li(u5,

и $ -

 

1 — v/2£2

 

 

2

Л2“' ’

 

 

 

 

 

 

 

% (4 0*

^f) = L2 («?,4

 

-

Ц р

i- 4 -Л -d-

^ .

 

 

 

 

 

■*

л і da

л іЛ2

da

^ з і и\,ѵ% w])

~

 

 

 

h?

 

л

 

— L 3(«f,

u®, a if)+ -j^ 4 (ife )

wf 4"

 

/I2

d ^2

d

1

d

\

dw\

 

 

 

12AxAi

da 4

da

 

da 4

da

 

 

 

 

 

 

 

 

Я

 

 

 

«Г =°Г =wf —0 приs"< °. è№=-2іг J

<(ѳ)cos lkBde-

 

 

 

 

 

 

 

 

я

 

 

 

 

 

 

 

 

 

~~k

 

 

Для получения решений всех бесконечных систем, входящих в (11.126), до­ статочно знать решение системы нулевого приближения (s = 0), через которое можно выразить решения других систем.

Рассмотрим систему уравнений нулевого приближения подробнее. Опреде­ лив из второго уравнения (11.126) (s = 0) сумму по^ и подставив ее в пятое урав­ нение (11.126) (s = 0), получим

W-(/&)»

+ 2bIk [ У г + ж ) АгА2 da <Ѵ8) -

 

12А\

1 — V

d

»8 = о (t = 1,2,...).

4iy43

da

(11.127)

 

 

Далее, пользуясь (11.127), находим суммы по Іъ входящие в первое, второе, третье и пятое уравнения 1 (11.126) (s = 0). Подставив полученные суммы в пер­ вое и второе уравнение (11.126) (s = 0), получим систему 2 уравнений относи­

тельно ц° и Шц :

Lu “§ + LIA = °'

(IU28)

^ г и о + ^ < = °-

Решив эту систему, можно определить как ш°г из (11.127), так и о® из четвертого

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

1 Из четвертого уравнения (11.126) (s=0) следует, что ы?=0 при 1X^1. 2 Дифференциальные операторы Z-slt s2 (sA, s2= l ,2) здесь не приводятся.

90

Соседние файлы в папке книги из ГПНТБ