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

книги / Оболочки и пластины

..pdf
Скачиваний:
19
Добавлен:
19.11.2023
Размер:
71.66 Mб
Скачать

М а

— [(1

Vi©!/*) cos2 а -ф- (01/*

vx) sin2 а] In р ,

 

4 я

 

 

ра'

Р (е,А -уг)( 1 - е 1/.)

sin 2а In р.

4ле3/*

 

 

Компоненты напряженного состояния, как видно из этих формул, имеют ярко вы­ раженный эллиптический_характер относительно а — угла между осью £ и произволь­

ным радиусом-вектором г (£, <р).

Полагая в (3, 2, 58) S = E 2/E i= 1, v1= v 2=v, получаем известные асимптотические -формулы для изотропной цилиндрической оболочки:

 

 

 

 

PR2

 

 

 

 

 

 

 

 

ш==" 1

^

рЧпр’

 

 

 

 

 

Т

= — Т

=

---- — In

р ,

 

 

( 3 , 2 , 5 9 )

 

 

 

р

 

nR 4

 

 

 

 

 

 

М р = М а = — (1 -ф- v) In р.

 

 

 

5.

Примеры. Цилиндрическая и сферическая оболочки

 

Речь

идет

о математической

аналогии

в задаче

о действии

сосредоточенной

нагрузки

на

ортотропную

цилиндрическую

оболоч­

ку [55].

 

 

 

 

 

 

 

 

нагрузки на

1) Рассмотрим действие радиальной сосредоточенной

изотропную

цилиндрическую

оболочку. Уравнения равновесия

такой

оболочки получим из (3, 2, 23), полагая А=0 = 1, е=0:

 

 

 

у4? + m ~

+

dtp2

= i2b\RPb а -

£0,Ф),

 

 

 

 

dg2

 

 

 

 

 

 

 

 

 

2b2.= — V 12(1—V2).

 

 

(3,2,60)

 

 

 

 

h

 

 

 

 

 

 

Решение аналогичной задачи для ортотропной цилиндрической оболочки сводится к определению комплексной функции уравнения (при G= Go)

( +

0,/‘ Т гТ

? + i2b2 Ц

+ 0 Т Т i2m P b (S- So. Ф). (3.2,61)

\д Ъ 2

дф2 /

д12

дф2

Приведенный модуль G = Go позволяет установить следующую физиче­ скую аналогию между изотропией и ортотропией:

Е<->УЕгЕ2,

(3,2,62)

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

Если цилиндрическая оболочка подвергается осесимметричной де­ формации, то уравнения равновесия для изотропной и ортотропной оболочек совпадают по форме (независимо от G), если установить соответствие

2 Ь] 2Ь\

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

6 .~ в 'Ч ,

(3,2,63)

учитывая при этом (3, 2, 62) и следующее свойство 6-функции:

 

б ("л’ ’ <р) = л б ^ ’

 

Далее, устанавливая соответствие между величинами

 

Ы ^ b2/®1/*, Р, ~ P/0V<,

(3,2,64)

получим полное соответствие между уравнениями (3, 2, 60) и

(3, 2, 61).

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

ции F для изотропной и ортотропной оболочек тождественно совпадают, при этом любое решение задачи об изгибе ортотропной цилиндрической оболочки может быть получено из решения аналогичной задачи для изотропной оболочки. Пересчет результатов осуществляется по форму­ лам (3, 2, 62) —(3, 2, 64).

Иллюстрацией могут служить асимптотические соотношения (3, 2, 58) и (3, 2, 59), которые переходят одни в другие при помощи формул (3, 2, 62), (3, 2, 64).

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

мерного прогиба от безразмерной

длины

цилиндрической

оболочки,,

свободно опертой по торцам [40]:

 

 

 

w*0 = f( l 1*)-

 

 

Пользуясь формулами (3, 2, 54) — (3,

2, 62), получаем

такое соот­

ветствие:

= 0'Л^.

 

ш*0 = @>uWoy

 

Для проверки этих формул можно использовать числовые данные, при­ веденные в табл. 3.1 и 3.2.

 

 

Т а б л и ц а 3.1

 

 

Т а б л и ц а 3.2

1

11 880

41 170

63 370

97 270

1

1484

| 5147

7921

12150

а

1

10

20

40

а

0,5

5

10

20-

При действии двух самоуравновешенных сил на бесконечно длин­ ную цилиндрическую оболочку В. М. Даревским [32] выведена простая формула для максимального прогиба. Эта формула на основании ана­ логии для ортотропного материала будет иметь вид

(3,2,65)

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

Пусть на ортотропную оболочку с относительной длиной, удовлетворяю­ щей условию (3, 2, 66), действуют m равных по величине радиальных сосредоточенных сил, расположенных с одинаковыми интервалами на направляющей окружности. Тогда при

3 < m < - i

3 (1 — VjVg)

0

 

прогиб w0 будет равен

 

wп = пгя^0т г - ( 1<646 + 1,6:17 + 1,905 -L-j [3 (1 -

При значениях т , удовлетворяющих неравенству

пг >

з ( l - z w L f u

можно пользоваться формулой

W0 = — е- ^ 1(cos Ь\ -f sin b 11 \),

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

Для свободно опертой по торцам короткой цилиндрической оболоч­ ки при действии одной радиальной оси, приложенной в середине про­ лета, обобщается формула, полученная в статье [56]:

Юо = 1.6291/2-/2

(2я) '2 L ® J

Изложенный выше метод применим не только к задачам прочности,, но и к задачам устойчивости и колебаний ортотропных оболочек [55].

2) Рассмотрим действие нагрузки локального типа на сферическую оболочку.

Общий интеграл уравнений равновесия тонкостенной сферы был исследован в работе А. Л. Гольденвейзера [57], который доказал возможность разложения истин­ ного напряженного состояния сферической оболочки на три состояния — безмоментное, чисто моментное и смешанное, что значительно облегчает получение окончательного решения. В этой же статье точно решена задача для замкнутой сферы, находящейся под действием произвольно расположенной сосредоточенной силы, уравновешенной в полюсе. Основываясь на этих результатах, В. 3. Чернина [69, 70] построила прибли­ женные решения для изотропного сферического сегмента, находящегося под действием сосредоточенной силы или момента. Г. В. Чернышев [71] получил решение для сосре­ доточенной нагрузки, действующей на сферу переменной толщины и анизотропную сферу постоянной толщины. Лекки [64], используя теорию Хаверса, привел приближен­ ное решение для осесимметричных сосредоточенных нагрузок.

Точные уравнения равновесия для непологой сферы тбыли получены И. Н. Векуа* [58, 59] в комплексной форме и В. 3. Власовым [60] в действительной форме. Общее

решение уравнений В. 3.

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

в докторской диссертации

В. Г Рекача и приведено в статье [68]. Частные интегралы

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

даны в статье Н. Г. Гурьянова [79].

Задача об изгибе изотропной пологой сферической оболочки при различных видах осесимметричного нагружения подробно исследована М. Н. Ручимским [62]. В работах

[66 и 67] получено решение аналогичной задачи для

равномерно нагруженной оболочки,

а также для случая, когда на оболочку действует

сосредоточенная

сила, приложенная

в полюсе. Результаты сравниваются с экспериментом, проведенным

авторами. Частное

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

ного контура получено в работе [73].

изгибе

сферической

оболочки

Н. Г

Гурьянов [75, 76] дал решение задачи об

нагрузками, распределенными по площадке с центром

в произвольной точке,

а также

по участку

меридиана или параллели. Ю.^П. Артюхин

[39, 42]

исследовал

задачу о

прочности анизотропного однослойного и слоистого изотропного сферического сегмента, нагруженного сосредоточенной силой. Экспериментальному исследованию прочности и устойчивости сферического сегмента при действии на него сосредоточенной силы или нагрузки, распределенной по круглой площадке, посвящена работа [44].

Задача стыковки сферической оболочки с цилиндрическим патрубком бесконечной жесткости исследована в работе Бейларда [61], причем ось патрубка совпадает с осью вращения сферы. Джонсон [77] обобщил задачу Бейларда на случай наклонного дефор­ мируемого патрубка. Туус [63] опубликовал экспериментальные данные, подтверждаю­ щие теоретическое решение Бейларда.

Точное решение задачи о действии на сферическую оболочку произвольно распо­ ложенной и гармонически меняющейся во времени сосредоточенной силы с учетом по­ перечного сдвига и инерции вращения получено в [74], причем особо рассмотрен случай статического нагружения оболочки сосредоточенной силой.

В работе [47] решается задача о предельном равновесии сферического купол? пере­ менной толщины, на который действует внутреннее давление и сосредоточенная сила, приложенная в полюсе. Вместо обычного методаоценки напряженного состояния пред­

лагается находить критическую нагрузку по диаграмме «сила — перемещение»,

которая

строится с помощью шагового метода. Экспериментальная проверка

метода

осущест­

влена в [78].

 

 

 

'При изложении задачи будем опираться на статьи [73 и 76].

Уравнения В. 3. Власова для пологой сферической оболочки

 

- ^ Г V2V2(P + v2ay = °» V2(P — D R k \2v 2w +

-^ -Z =

0,

(3,2,67)

П.П

Rr

 

 

где

 

 

 

м _ 12(1-У 8)

 

 

 

RW

 

 

 

легко сводятся к одному уравнению относительно комплексной функции

F = Dk2w ---- — ср.

R

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

V2V2? - i V 2? = у Ы д г - б , 6 -Y ),

(3,2,68)

где

I, 0 — У) = Ь(х— I) £ 6(6 — у — 2пп).

П= — оо

Используя известные свойства б= функции (3, 2, 37) и разлагая функцию F в ряд

? = У. 6„Fnc o s/i(e — у),

ьп = {

1 , п

° ’

(3,2,69)

!ёо

^

2, п >

0,

 

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

 

 

 

 

V Svf. - Ф Л -

8<* - »•

 

{3,2,70)

Частный интеграл этого уравнения можно отыскать с помощью ин­ тегрального преобразования Ханкеля. Для этого (3, 2, 70) умножаем на xln(ux)dx и интегрируем в пределах от нуля до бесконечности.

Изображение Ханкеля для Fn записывается в таком виде:

 

 

оо

 

 

 

К

(и) = |

xln (их) Fn (X) dx.

 

Можно показать [81], что

 

о

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

j х 1 п (их) 4 2JFn (х) dx = ~ u2F'n (и),

 

О

 

 

 

 

 

 

00

 

 

 

 

 

 

1 xln (их) 4 nsj2nFn (х) dx = и4?* (и),

 

О

 

 

 

 

 

 

оо

 

 

 

 

 

 

р • J */„(ид:)6 ( x - t )

= - ? - I n(ul).

 

2я£

 

 

 

 

 

 

Отсюда получается формула для Fn

 

 

 

 

 

V

Р *

M “S)

 

 

 

 

а4

ш2

 

Переходим к интегралу обращения:

 

 

?„(*> = —

Г и1^ аШ ч ± йи,

 

пК

'

J

u2 (u2

О

 

Подставляя это выражение в (3, 2, 69) и используя формулу Ней­

мана для цилиндрических функций [79]

 

 

/о И ) = £

 

вя/„ (ug) /„ (их) cos п (6 — у),

 

п=0

 

 

 

 

 

где

 

 

| 2 _

2xl CO S ( 6 — у ) ,

(3,2,71)

1/2

=

* 2 +

получаем частное решение разрешающего уравнения, соответствующее сосредоточенной силе:

£ _Р_ Г /0 (иу) du __ iP f UIQ (uy) du

Xх' ~

2n .] и (в* 4- i)

2n \

и14- i

 

О

 

 

- ?

i

{К, (VTi/) + 1п!,).

J и )

Для более полной расшифровки этого решения можно использо­ вать теорему сложения для функции Макдональда [80]

/Со (1/Ту) = 2

V » (/Гб)/Сп (Vfx) cos п(6 -

б)

(*>

1),

(3, 2,73)

 

п

 

 

 

 

 

 

/Со (/Гу) =

£

6Д„ (/Гб) 1п( / Глг)cos п (9 -

Y)

(■* <

6)

& (3,2,74).

 

п

 

 

 

 

 

 

и формулу разложения для 1пу [73]

 

 

 

 

 

 

 

 

00

 

 

 

 

 

1п у =

In лт —

——

cos я (0 у)

(•«>!),

(3,2,75)

 

 

 

П=\

 

 

 

 

 

 

 

СО

,

 

 

 

 

 

lny = 1п| — ^ - i - ^ Y ^ ncosn(e — Y)

(*< б).

 

(3,2,76)

п = \

Чтобы получить решение для нагрузки, равномерно распределенной по окружности радиуса а с центром в полюсе, необходимо (3, 2, 72) про­ интегрировать по у от нуля до 2я, считая £= а. При этом следует поль­ зоваться формулами (3, 2, 73—3, 2, 76).

Для х> а получаем

Fi(x) = -тг- {/о(/Га)/Со(/Г*) + In*},

(3,2,77)

ZJt

 

а при * < а

 

?i(*) = ^ К о ( / Г а ) /0(/Г *) + 1п|}.

(3,2,78)

Чтобы получить частное решение уравнений (3, 2, 67), соответствующее нагрузке, равномерно распределенной по круговой площадке с центром в полюсе, необходимо взять двойной интеграл от (3, 2, 72) в пределах для у от 0 до 2л и для £ от 0 до а. Если х> а, интегрирование проводим, пользуясь формулами (3, 2, 73),(3, 2, 75). В итоге получаем

 

?2 = i[ a— (jfa

( / Й

+ -у - In*] ,

(3,2,79)

где q — интенсивность действующей на площадку нагрузки.

(3, 2, 79)

Если х < а , то х может быть и меньше и больше g, поэтому

приходится разбивать на два интеграла:

 

 

 

х

 

а

 

 

 

? = = ^ Л

d ' f "

[К° i V f y ) + 1 п у] d l + I “S " (* °

у Г у ) + l n ^

d y -

0

0

 

л:

 

 

 

Интегрирование

проводим с учетом

того,

что

в первом

интеграле

х>% и нужно пользоваться формулами (3, 2, 73),

(3, 2, 75), а во втором,

поскольку х<£, применяем формулы

(3, 2, 74) и

(3, 2, 76). Окончатель­

но получаем

 

 

 

 

 

 

Результат (3, 2, 8Й) получен с использованием

соотношений

(3, 2, 76)

ъ

 

_

;п

dzo(/Та)

 

 

 

 

 

_ ,•(,

Д.У<*)

 

f xZ, (VTx)ix=ia J h M J A .

 

J

 

 

 

da

 

db

 

где

 

 

 

 

 

 

 

>(УТх)

=

j

 

/ o ( V / i

x

) dz,{/tx)

=

'

1

 

X)

da

 

dx

 

In (V i* )

f{

iX)

-

Kn (V i x)

dIn (V ix )

 

 

dx

 

 

 

dx

 

 

Если выбрать новый полюс сферы в точке, отстоящей от старого на расстоянии g, и считать, что расстояние от нового полюса до произволь­ ной точки равно у, то, заменив в формулах (3, 2, 77) —(3, 2, 80) х на у, лолучим решение, соответствующее несимметричной относительно оси вращения нагрузке. В этом случае для сферы, нагруженной по окруж­ ности радиуса а с центром в точке g, будем иметь

 

? 1 ( У ) = - ^ { 1о(УТа)К0(УТУ) + \пу}

( у > а),

(3,2,81)

 

Ъ (У) = ^ - { K 0(VTa)I0(VTy) + \n l\ (у<а).

 

 

При

нагрузке,

распределенной по круглой площадке того

же

ра­

диуса с центром в g, приходим к такому результату:

 

 

 

 

 

F*(У) =

— °"(£ -‘а)- (Viy) +

In j/j

> а ),

 

 

F2 (у) =

<L[ 1 + a

/0 (/Т у) +

( in а —

L )

+ - f ]

(у <

а).

 

 

 

 

 

 

(3,2,82)

Поскольку х и у связаны зависимостью (3, 2, 71), можно пользоваться соотношениями (3, 2, 73) — (3, 2, 76). ^

Если теперь отделить вещественную и мнимую части в Fj (/=1, 2), получим выражения для прогиба и функции напряжений, через которые выражаются все компоненты напряженно-деформированного состояния

оболочки.

Общий интеграл однородных уравнений (3, 2, 67)

00

= У 6Л(Впхп+ Clnber„ X + с2пbei„х) cos п (0 — у), (3',2,83)

п

2nF*L - = у 6Л(Апхп + Clnbei„х — С2пЬегл х) cos п (0 — у),

PR

л=0

где W = /1 2 (1 — у») PR

2nEh*

Если к (3, 2, 83) добавим соответствующее-решаемой задаче част­ ное решение и найдем константы из граничных условий на внешнем контуре сегмента, то задачу можно считать решенной. Для нахождения констант частное решение нужно представить в таком же виде, что и (3, 2, .83). Для этого необходимо воспользоваться одним из соотноше­ ний (3, 2 ,73) — (3, 2, 76).

Если внешний контур х = х0 жестко защемлен, то должны выпол­ няться следующие условия:

w (х0, 0) =

(лг0, 6) =

и (х0, 0) =

V (х0, 0) =

0.

 

 

 

I

 

 

 

 

 

 

 

 

Находя перемещения из формул

 

 

 

 

 

 

 

 

и =-

"

+

£ 6 "

1 ^

т

г

д:”

+ ' С

0 Е

' 1 ( в

"

' |' ) ]

kR [

 

 

 

п=0

 

 

 

 

 

 

 

 

 

V = — у

---- —

----- Xn+1

s

i n

n ( 0

у )

+ k (1

- f

v )

kR

L l

(n+1)

 

 

 

 

 

Eh

 

(3,2,84)

дф

+ мдх ^

1 <?ф

7 ¥

и подставляя их значение в (3, 2, 84), приходим к системе четырех уравнений относительно констант 5 n, Cin, Сгп, А п. Эта система может быть сведена к системе трех уравненйи относительно £ n, Cin, Сгп и от­ дельному уравнению для Ап.

Для случая сосредоточенной силы, приложенной в точке, константы имеют вид [73]

А В п= IQ (/3/9---/4^10) ----^1 (^7^9

V ll)

^2 (^7^10---

Vll)»

ЛС1л =

X Q

(/7/9

^ 11)

^6 (^5^9

 

 

^2 ( ^ 1 1

 

Д ^ 2 Л =

(^3^11

^7^1о)

^1 (^6^11

^7^в)

~Ь ^6 (^6^10

^З^в)'

И

 

 

 

 

 

 

 

 

 

/IATQ -

=

~

^8^1/I

^4р1п +

^3^2Л “1" ^12>

 

где

 

 

 

 

 

 

 

 

 

Д —

(/3/9

Ij'io)

(W9

^4^в) *~f~^2 (^6^10

 

Здесь введены следующие обозначения:

 

 

 

=

Ьегл*0,

/2 —bei^^o»

^3 =

 

^4 =

 

/5 = пх*~\

/б = — (bei„£ кеглх0 +

Ьег„^ keirtл 0),

'» = -

(beiЛ кег>0 + Ьег, 6 kei>0),

I,

=

 

/9 = —Ьег' х0+

— Ьег„ х0,

l1Q= bei' х0---- - bei„ x0t

"

х0

xt

In = [b^r„ ^ к е г ^ 0 — 3L кег„ *0) —bei„| (kei'nx0 — j - kei„

l12= (ber„ l ker'nx0 - bei„ | kein*0) +

.

На рис. 3. 15 приводится график зависимости величины прогиба в точке приложения сосредоточенной силы от положения последней на меридиане.

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

§ 3. МЕТОД К. 3. ГАЛИМОВА ИНТЕГРИРОВАНИЯ НЕЛИНЕЙНЫХ УРАВНЕНИЙ РАВНОВЕСИЯ ПОЛОГИХ ОБОЛОЧЕК

Рассмотрим задачу определения больших прогибов пологой прямо­ угольной в плане цилиндрической панели, находящейся под действием равномерного поперечного давления q. Метод решения основан на раз­ ложении в ряды Фурье искомых функций и последующем точном удов­ летворении всех граничных условий и условия совместности деформа­ ций. Уравнение равновесия при этом удовлетворяется в среднем по, методу Бубнова — Галеркина. Этот метод интегрирования довольно общий и может быть успешно распространен на прямоугольные в плане пологие оболочки, подверженные действию различного рода локальных и переменных нагрузок.

1. Цилиндрическая панель

Функция напряжений Ф и функция прогиба w удовлетворяют систе­ ме нелинейных уравнений

 

у 2у2ф —Eh

d2w

\ 2

d2w

d2w

 

1

] ,

(3,3,1)

 

дхду

)

дх2

ду2

 

 

 

 

 

R *~дх7

 

D vW w =

д2Ф

d2w

, г д2Ф "/

(Pw ,

1 \) - 2о

д2Ф

 

d2w ' "Ь 9»

(3,3,2)

V V

ду2

дх2

 

дх2 \

ду2

R )

 

дхду

 

дх ду

 

где у 2 — оператор Лапласа, D — цилиндрическая жесткость, R — радиус

панели, Е — модуль

Юнга,

h — толщина

панели,

q — нормальное дав­

ление, действующее на поверхность панели, х — координата, направлен­

ная на образующей панели,

а // — координата,

взятая на

дуге панели.

Начало координат поместим в центре панели,

так

что

уравнениями

кромок будут х =±а , у=±Ь.

 

 

 

 

 

Граничными условиями задачи будут

 

 

 

7 \= 3*<t>—

P i , S

 

= 0 при х =

±

а,

(3,3,3)

а«/2

= Р1,

s =

-----дх ду

 

 

 

д2Ф

■- />„

S =

-

д2Ф

 

 

 

дх ду

 

 

 

дх2

 

 

 

 

 

 

и условия свободного опирания

w =

0,

= 0

при

х = ± а,

(3,3,5)

 

 

дх2

 

 

 

w =

0,

- f y = 0

при

у — ± Ь.

(3,3,6)

Здесь Р\ и Рг—заданные растягивающие — сжимающие усилия на кром­ ках х = ± а и у = ± Ь , а касательное усилие S принимается равным нулю.

Если взять w в виде ряда

0 0

оо

 

 

w = ^

2 Ann С0£? 2т~ 1

2П~ - ПУ-

(3>3>7)

ш=1 я=1

 

 

т п — постоянные), то граничные условия (3, 3, 5) и (3, 3, 6) удовле­ творяются. Чтобы удовлетворить статическим граничным условиям точно, необходимо использовать решение типа краевого эффекта. Поэто­ му функцию напряжения возьмем в виде ряда

 

 

 

 

 

Ф =

\

РгУ2+

- у Р**2+ Ф (*, У),

 

(3,3,8)

где

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф(*. у) = 2 [

C*Ek (*) cos -

f f

+

DkEl (у) cos - I f -

+

 

 

 

 

 

k = \

L

 

 

 

 

 

 

 

 

 

 

 

 

 

+

(X) cos

2k~

l ny +

Р*ф04_, (у) cos - 2k~

1 - Itxj +

 

OO

oo

 

mux

tiny

 

oo

oo

 

 

 

 

2n 1

,

V I

VT

r,

,

\^ l

1 П

r?

 

2 m — 1

 

+

2 J

2 J

 

C 0

S C

0

S

+

2

J

2 J

FmnCOS' '

2a

 

-----Yb-----n y ’

 

m = 0 n= 0

 

 

 

 

 

 

 

 

 

 

m — 1n — \

(3,3,9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

причем здесь введены обозначения:

 

 

 

 

 

 

 

 

 

 

 

 

и

 

knx

 

 

knx

 

 

 

 

 

 

 

Ек(х) =

sh—г ~

 

ch ■

f

+

 

 

(3,3,10)

 

 

 

knx

 

b

 

 

 

 

 

 

 

 

 

 

 

sh

kna

 

sh- kna

 

 

 

 

 

 

 

 

, 2k — 1

, 2 k 1

 

 

 

2 k — 1

2 k — 1

 

 

 

 

 

x sh — —— nx ch — ——

na — a sh — ——

na c h -----------nx

f 2k—I (* )

 

2b

 

 

2b

 

 

 

 

2b

 

2b

, (3,3,11)

 

 

 

 

 

 

,

2k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh -----------na

 

 

 

 

а функции

E\ (у)

и (flk—i (у)

получаются соответственно

из функций

Ек(х) и

f 2k-i(x) путем замены х на у ,,а не b на а.

 

 

 

Эти функции удовлетворяют условиям

 

 

 

 

 

 

 

 

Ека) = Ек(0) =

0,

Е°к (±Ь) = E°k (0) =

О,

 

 

 

 

 

 

/2t-i (i а ) — О Ф24-1 (i b) = 0.