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

книги / Общая теория анизотропных оболочек

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
21.14 Mб
Скачать

 

ЛИТЕРАТУРА

341

26. О г и б а л о в

П. М., К о л т у н о в М. А ., Оболочки и

пластины,

Изд-во МГУ, 1969.

оболочек,

Здесь имеется

достаточно полная библиография по теории

§ 12. Результаты, приведенные в этом параграфе, интересны с точки зре­ ния расчета анизотропных слоистых оболочек со слабыми механическими характеристиками в поперечном направлении и, в частности, с точки зрения расчета оболочек, изготовленных из армированных пластиков. Необходимые сведения об этом материале с точки зрения теории оболочек можно найти, например, в работах:

27.А м б а р ц у м я н С. А ., Специфические особенности теории оболочек из современных материалов. Известия АН АрмССР (механика), т. 21, № 4, 1968.

28.

О г и б а л о в П. М.,

С у в о р о в а

Ю. В.,

Механика армированных

29.

пластиков. Изд-во МГУ 1965.

 

В. П.,

Т е т е р е Г. А .,

Сопро­

М а л м е й с т е р

А. К. , Т а м у ж

30.

тивление жестких

полимерных

материалов.

Изд-во «Знание»,

1967.

Т а р н о п о л ь с к и й

Ю. М.,

Р о з е

А.

В., Особенности расчета

31.

деталей из армированных пластиков. Изд-во «Знание», 1969.

В. В .г

А л е к с а н д р о в

А.

Л. , Б о р о д и н

М.

Я. , П а в л о в

 

Конструкции с заполнителями из пенопластов. Оборонгиз, 1962.

 

В этих работах читатель найдет достаточно полную библиографию по­ этому вопросу.

Рассматриваемый параграф написан согласно следующим работам: [3], а также

32.А м б а р ц у м я н С. А ., О двух методах расчета двухслойных ортотропных оболочек. Известия АН АрмССР (ФМ науки), т. 10, № 2, 1957.

33.А м б а р ц у м я н С. А ., К расчету двухслойных ортотропиых оболо­

34.

чек. Известия АН СССР, ОТН, № 7, 1957.

 

А м б а р ц у м я н

С. А., К общей теории анизотропных оболочек.

35.

ПММ, т. 22, в. 2, 1958.

Д. В., К теории орто-

А м б а р ц у м я н

С. А., П е ш т м а л д ж я н

 

тропных оболочек

и пластинок. Известия АН

АрмССР (ФМ науки),

т. 12, № 1, 1959.

36.А м б а р ц у м я н С. А ., К теории изгиба анизотропных пластинок и пологих оболочек. ПММ, т. 24, в. 2, 1960.

37.X а ч а т р я н А. А ., К расчету трехслойной ортотропной оболочки. Известия АН АрмССР (ФМ науки), т. 12, № 5, 1959.

38.А м б а р ц у м я н С. А., Еще одна уточненная теория анизотропных оболочек. Механика полимеров, № 5, 1970.

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

параграфа читатель может пополнить, просматривая [3, 32-38 ], а также 39. А м б а р ц у м я н С. А., П е ш т м а л д ж я н Д. В., О нелинейной

теории пологих ортотропиых оболочек. Известия АН АрмССР (ФМ науки),

т. И, № 1, 1958.

40.А м б а р ц у м я н С. А., Теория анизотропных пластин. Физматгнз,. 1967.

§ 13. Основные положения термоупругости анизотропных тел можно найти, например, в книге 41. Витольд Н о в а ц к и й , Вопросы термоупругости, Изд-во АН СССР

1962.

§ 13, п. 1. Посвящен классическим задачам термоупругости анизотропных слоистых оболочек. Этот пункт написан па основании статьи:

342

Н А П РЯЖ ЕН Н О Е СОСТОЯНИЕ И ДЕФОРМ АЦИИ О БО Л О Ч ЕК

[ГЛ . XI

42. А м б а р ц у м я н С. А ., Температурные напряжения в слоистых анизо­ тропных оболочках. Известия АН АрмССР (ФМЕ и Т науки), т. 5, № 6, 1952.

§13, п. 2. В этот пункт вошла лишь часть статьи:

43.Д у р г а р ь я н С. М. , К температурному расчету тонких ортотропных оболочек вращения. Инженерный журнал, т. 2, в. 3, 1962.

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

§13, п. 3. Написан на основании статьи:

44.A m b a r t s u m i a n s . A., D u r g a r i a n S . М., Some Thermoelastic Problems of Anisotropic Shells and Plates. Non-Classical Shell Problems. Proceedings of the IASS Symposium, Warsaw, 1964.

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

45.К о в а л е н к о А. Д., Пластинки и оболочки в роторах турбомашин. Изд-во АН УССР, 1955.

46.Б и р г е р И. А ., Круглые пластинки и оболочки вращения. Оборонгиз, 1961.

47. П о д с т р и г а я Я. С., Я р е м а С. Я., Температурные напряжения

в оболочках. Изд-во АН УССР, 1961.

48.Ф е о д о с ь е в В. И., Прочность и теплонапряжениость узлов жидкост­ ных ракетных двигателей. Оборонгиз, 1963.

49. О г и б а л о в П.

М., Г р и б а н о в В. Ф., Термоустойчивость пла­

стин и оболочек.

Изд-во МГУ, 1968.

Г л а в а III

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

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

§ 1. Свободные колебания

Уравнения движения элемента оболочки da dp dj в ортого­ нальной криволинейной системе координат а, р, f, согласно (16), запишутся следующим образом:

 

< в д . „ + < * л я * + < а д д . т-

v * ... +

 

 

 

 

 

д-а

 

 

 

+

=

р

 

 

(tfM i* + (В Д ^ ),т +

*Ъ>.. -

° Н 1:Р+

 

 

 

+ v ^ . r + а д . = р 5 г а д .

(1. 1)

 

 

 

 

 

- ° „ а д . т- « , а д . т= р ^ г а д .

 

где

р — плотность материала,

t — время.

 

 

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

на

уровне уточненной

теории (см. гл. I, § 7), с учетом явлений

поперечных сдвигов, имеем

 

 

 

 

« . = ( 1 + * п ) « - ^ + г ( 1 + т | ) | Ч -

 

 

-

Г 3 ( l + Т т ) Т Ф 1»

ф 1 = а 55? +

 

 

Вр = (1 + * 2 Т )У - 5 - д ^ + т ( 1 + Т т ) т ф. -

( 1 - 2 )

 

 

 

 

 

ф2= а44ф -f- а45ср,

 

 

u^ = w(а, р, t),

 

 

 

 

где а=и (а, р, t), v~v (а, р, /), <р=<р (а, р, t), <|^=<|>(а, р, t).

-344

К О Л ЕБА Н И Я И УСТОЙЧИВОСТЬ О БО Л О Ч ЕК

[ГЛ . III

Умножая уравнения (1.1) на dy, а кроме того, первые два уравнения на у dy и интегрируя их в пределах от —hi2 до hi2 с учетом (1.2) и (1.1.15), получим следующие уравнения движе­ ния элемента оболочки:

(ВТ,), -

Д „7’1 +

(,l.S,)i[ +

, l iii.?, + АВк,,\\ =

 

 

 

А , ^ Т 1 + ( ^ l ) , « +

 

^ . а 5 2 + ^ ^ 2 ^ 2

 

 

 

,

„ Г , d*v

А З , ,

,

,

Ч 1

дЗш

,

А5 / 3

 

,

, Л

^ ф 2 1

?AB\i h — — u (kl +

k,) в

^ j- , +

ш \2

+

 

J•

- ( К Т г +

W

[

(

B

W ^

- f

(ЛЛГ2)

 

 

 

(1.3)

(fiA^i).. +

) , э +

Л „Я, -

£ аМ2-

ЛВЛГ, =

 

 

 

 

 

. чпГ*»,о».

,

ч<?2И

 

АЗ 1

<)3ц, ,

А5 *Ф,П

 

— рАВ

|^12

(2kt

-f- k j df2

 

12 Ada^ 2+

12l)

dt, J,

Л/2),р - f (Sffi),.+ ^.A -

^.A ABN2 =

 

 

 

 

 

Л о Г А З /о ,

 

, ,

4 d?l>

 

АЗ 1

^Зц;

A5

* Ф Л

 

— рЛ5 ^

(2A2- f

dt2

 

J2B

 

dt*

120 dt* J *

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

В частном случае классической теории, пренебрегая явлениями

сдвигового

происхождения

(Ф .= 0 )

и

ограничиваясь точностью

fc(h, вместо правых частей уравнений (1.3)

(левые части остаются

неизменными) получим

соответственно

 

 

 

 

7л п д2и

 

1 л

д2и

г, (^ w

 

 

PhABW '

PhABW>>

Ph ^dt*'’

(1.4)

 

>

3

D

d3w

 

A3

. d‘3w

 

 

 

 

P 12^ dx t><2*

 

Р12Л ^^2-

 

В случае многослойной оболочки, собранной из произволь­

ного числа анизотропных слоев (см. гл. I,

§ 11), в правых частях

уравнений

движения

классической

теории будем иметь

 

 

г

Л В

д2ц

v т> дЗш

 

 

 

LeAts

dt2 '

К ?В dx dt2

 

 

 

С АВ d2v

 

е, .

d$w

d*w

 

 

L?Atf

dt2

 

Л РА

d$dt2Le

dt2

(1.5)

 

 

 

 

 

 

 

 

 

 

 

V

 

я ж

-

D*B S

*

 

 

 

 

 

 

 

 

 

 

d*w

 

 

 

V

 

5

dt'2

 

Ifd i*’

 

§ 1] СВО БО ДН Ы Е КО Л ЕБА Н И Я :м:>

где для приведенных характеристик плотности слоистой оболочки

имеем (рис.

34)

С9=

Я1

 

2

p . ( 8 . - V i ) >

г

«— 1

 

( 1. 6)

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

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

ным уравнениям движения (1.3) с соответствующими

правыми

частями (1.3), (1.4), (1.5)

соотношения

упругости соответственно

(1.7.19),

(1.1.24)—(1.1.26), (1.11.2),

геометрические

соотно­

шения

(1.1.5) — (1.1.8),

граничные

условия (1.7.27)—(1.7.29),

(1.1.27)— (1.1.30), а также начальные условия в момент времени t = 0, получим полные системы, с помощью которых могут быть опреде лены динамические характеристики оболочки.

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

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

Рассмотрим некоторые частные задачи по определению свобод­ ных колебаний анизотропных оболочек. Эти задачи имеют целью

1546 КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК £гл. in

«светить специфические особенности колебаний анизотропных оболочек.

1. Определение частот колебаний шарнирно опертой по всему контуру ортотропной цилиндрической панели. Исходная система дифференциальных уравнений колебаний круговой ортотропной цилиндрической (7?1=оо, R2—R) оболочки с учетом лишь инерци­ онных сил, происходящих от нормального перемещения, согласно приведенным выше соображениям, системе уравнений (1.7.56), л также следующим представлениям:

 

W

Jo (a55^U

* 44^22)

 

a44a55 1(Х)(^п^гг 12)^J

 

? = - [ £

+ £ ^ 2 - « 4 4 £

 

£ ( ^ 2 2 - В Д ] V .

(1.7)

 

 

 

 

Ф =

^22 + ~а ^12 — «55 IQ

 

 

“ ^ 12)]

 

запишется следующим образом:

 

.

д*

г

 

 

 

1

д*Р ,

лз г ^

0

д* ,

 

 

 

 

 

А

да2

\j^ ь 11~гА"хглг bi

-f-—

ь«»

 

 

 

12

 

дад$

12

 

Г

 

 

 

 

 

 

 

w ( “ «

 

 

 

 

 

 

 

 

 

 

 

+

Той

Т О - (в «^ п + в м ^ +

 

 

 

 

 

 

(1.8)

 

 

 

 

I

 

 

 

 

 

Г 2

^

 

 

 

 

h i

 

/ Т

Т

 

П

 

 

 

“ Г* «44«55 100 1 ^ U “ 22

" u J J

<>f2

— *■'»

^

( ^

22- £ У

-

w

( fl^

H

+ a« L 22) +

 

 

 

 

+

«44«55 Щ

 

 

 

 

 

(?2ф-

: 0,

 

 

 

 

 

J 22

 

да2

где F = F

(а,

р, /), ЧГ= Ч Г (я,

р,

t) — искомые функции, р= T0/g

( f 0 — удельный вес материала оболочки, g — ускорение силы тя­ жести) ,

г

D

I D

Г __ /)

22^2

| 73 д2

 

 

Ь 22— ^

" Г ^ 12^ 2 .

 

 

 

 

 

(1.9)

^12 = (®12 ~Ь ®6б)

^ » °>9 = ® 11®22

® 12

К системе уравнений (1.8) следует присоединить граничные условия, которые запишутся следующим образом (рис. 59):

при я = 0 , я

w—-v — Ml = T 1= ^ = 0,

при р = 0, р = Ь

w = и— М2 = Т2 = <р= 0.

1]

 

 

 

 

 

СВО БО ДН Ы Е КО Л ЕБА Н И Я

 

 

347

Решение системы (1.8)

представим в форме

 

 

 

 

 

F(a, р,

t) == Fmtt sin Х;а sin Х2р cos o>mt,

|

 

 

 

 

W(a, p, t) =

 

sin

sin X2p cos o>mnt,

J

 

гДе Fmn,

 

— некоторые постоянные, \=тт./а,

\~пп/Ъ, а

длина панели,

 

Ь — длина дуги

поперечного сечения панели,

штп — частота

колебаний,

/ л и г а — це­

 

 

 

лые числа, показывающие число полуволн

 

 

 

по продольным и поперечным направ­

 

 

 

лениям

соответственно.

 

 

 

 

 

 

 

Легко

заметить,

что

представления

 

 

 

(1.10) полностью

удовлетворяют

гра­

 

 

 

ничным

условиям

шарнирного

опира-

 

 

 

ния по всему контуру панели.

 

 

 

 

 

 

Подставляя

 

значения

F

и

 

из

 

 

 

(1.10) в

 

систему (1.8),

получим

сле­

 

 

 

дующую

 

однородную

систему

алгеб­

 

 

 

раических

уравнений

относительно

коэффициентов

и

Fmn4

-

^

» #

 

{А Л 4+

2 (Я » +

2Я6в) XfXf - f В22Ц +

 

 

+ ^ ( а д + а д ) [ а д 66м +

 

 

 

 

 

 

+

(ВпВ22-

в\2-

2ВпВ№) Цк\+

В22В№Ц]} +

 

 

 

+

 

<

«

{ i +

- г о - [ « * ( а

д +

а д

) +

 

 

 

+ « 4 4 ( а д + а д я + « л т (д п а д +

 

(1. И )

 

 

-j- (ВпВ22 Bj2

 

 

 

 

 

 

 

 

25125 66) Х|Х|-J- Л^ЛевХЦ} =

0,

Вщп дй6ш0 [ВпВ$в^{ “I- (ВцВ22В\2

2В12Вт) ХЩ- f

 

 

+

В22В66Ц] +

 

 

 

j l

+ - £

[ав5 (ВпЦ +

Я ввХ2)

+

 

Ч~ аи(В22Ц4“ Bggbf)]4“ «44«55Щ ) [ВцВмЦ4 " (ВПВ22

 

 

 

 

 

 

-

В\2-

2BnB J XfXf + B22B6e\*]} =

0.

Приравнивая нулю определитель однородной системы (1.11), получим для частот собственных колебаний ортотролной панели

348

К О Л ЕБА Н И Я И УСТОЙ ЧИВОСТЬ О БО ЛО ЧЕК

[ГЛ . Ш

«ледующее

выражение:

 

 

 

 

 

* - • =

C

c

21‘)ф‘=. L +>1Ь'«' f

+t i

+

 

+ ■ §

(« « Ч

+

%Ч>

 

 

 

 

( 1.12)

где

^2

 

 

 

 

 

^4

 

 

 

 

« 4А

2) +

 

 

b mn = 1 + -JQ- («55^11 +

а 44а 55

(^11^22 &1г)>

 

bu =

X2Вп ~j~Х“В66,

Ъгг =

Х2В22-J- Х2566,

(1.13)

^12 =

^1^2 (В12+ -®6б)-

 

 

 

 

Формула

(1.12)

получена

на основании уточненной

теории

с учетом поперечных сдвигов. В частном случае классической тео­

рии,

полагая ам=айй=0,

получим

 

Ч „ =

Щ-0[SuM + 2 (В12+

2Вт) ХЩ+ В22Ц] +

 

 

gBбб^о______________________ Х|_______________________

(1.14)

 

ТоR 2 В п В^).\ + (Вп В 22

В \г — 2B 12Bee) ХЦХ| + В 22В МЦ *

 

 

Сравнивая формулы

(1 .12)

и (1.14), замечаем, что поправка

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

Вчастном случае трансверсально изотропной оболочки, когда

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

Вп--- ^ 2 2 --- Я 1 2 + 2 В 66 ---- 1 ___V2 ’

^ 12---- 1 __V2 ’

#66 =

G — 2(1 + V)’

5 12 + В бв =

2 ( 1 — V)’

 

ВцВ^в =

Я А

~2

■*“ Bf2

2В12Вт) =

(1.15)

 

 

 

_

 

£2

 

 

 

_

2 ( l + v ) ( l — V2)

1

а44 — Я55 Q>1 V1 — V2 — V-

§ П

 

 

СВО БО Д Н Ы Е КО Л ЕБА Н И Я

349

Тогда

из (1.12)

получим следующее

выражение для

частот

колебаний

трансверсально изотропной цилиндрической

панели:

 

g D

М + Х | ) 8

, g E

Ц

(1.16)

 

тп

70*

1 + * * « 2 (Ч +

>1) ' Г 7оR2 (Ч + *1)2 ’

 

 

где

 

 

 

1 А2

 

 

 

D

 

Е к З

 

(1.17)

 

12 (1 1/2 ) ’

10 а-2 (1 — V2)G '-

 

 

 

Параметр h* представляет поправку к классической теории от учета явлений, связанных с поперечными сдвигами. Полагая h* равным нулю, из (1.16) получим формулу для определения частот колебаний изотропной оболочки:

g#

Ч

(1.18)

70Д2 (Xf + X|)2-

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

зотропных (трансверсально изотропных) оболочек

(E/G' > 1).

Для наглядного представления закономерностей

изменения

частот колебаний панели в зависимости от изменения параметра А* и подъемности оболочки (a/R) в нижеследующей таблице при­ ведены значения отношения шт„/й>тп, соответствующие некоторым

тонам

колебаний

(тге=1 ,

n = 1 ; т = 2 ,

п = 2)

трансверсально

изотропной панели (а=Ъ,

v = 0), вычисленные по

формуле

(1.16).

Рассматривая

расчетные

формулы

(1.12),

(1.14),

(1.16),

>(1.18)

и

таблицу 7, замечаем: с увеличением параметра h*,

т. е.

 

 

 

 

 

 

 

Т а б л

и ц

а 7

 

 

 

Я= со

 

Я — 2а

Я = а

 

л*

 

«и

 

0)27

«п

<22

ши

«И

 

 

®н

 

 

 

“и

«м

0,00

 

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

0,002

0,9807

0,9293

0,9837

0,9302

0,9891

0,9326

0,005

0,9539

0,8466

0,9615

0,8486

0,9742

0,8543

0,01

 

0,9137

0,7475

0,9282

0,7510

0,9522

0,7609

0,05

 

0,7092

0,4494

0,7637

0,4599

0,8481

0,4886

■отношений типа Bik/BiS, (E/G'), частоты колебаний уменьшаются -но сравнению с частотами, найденными по классической теории ,{h*= 0); различие между частотами колебаний, найденными по ^классической теории и по уточненной теории, менее заметное для

350 К О Л ЕБА Н И Я И УСТОЙЧИВОСТЬ О БО ЛО ЧЕК [ГЛ . Ш

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

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

[с2 (Д - f I)2 + 1 — h*Д] (Д + 2) w =

где

а оператор Лапласа на сфере, описанной в географической системе координат а, |3(A = R , B = R sin а) (рис. 60), имеет вид

 

 

Очевидно, отыскание

общего

 

 

решения

уравнения

(1.19)

сопря-

 

__

жено с некоторыми

трудностями.

 

 

Поэтому

будем

ограничиваться

 

 

лишь некоторым классом решений,

 

 

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

 

 

изотропной оболочки. Представим

Рис.

60.

решение уравнения (1.19)

в форме

w(а,

Р, *) = 1Р (а, P)sinarf.

(1.22)

 

 

Подставляя

значение w из

(1.22)

в

(1.19),

для

определения

W получим следующее уравнение:

[с2 (Д —f- I)2 —|—1 —й*Д] (Д + 2)Ж —

_ о>2 * g 2 - ( l — А» д ) ( Д 1 — v)W= 0.

(1.23)