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

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

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

Дельта-функция двух переменных равна произведению двух б-функ- цнй одной переменной разных аргументов:

б (* — *о. У — Уо) = Ь(х — х 0) 6 (у — у 0).

(3,2,18)

Если на оболочку в точке ui = ai, (Z2 = a2 действует сосредоточенная сила

величины P j , то выражение для проекции нагрузки'на ось / можно запи­ сать так:

fi(at — aj, 02— 02)

(3,2,19)

/ = 1,2, 3.

AiAi

 

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

м г

Чя = А\А% 6a, (a i — a i> а 2 — а 2)-

При действии сосредоточенного момента Л42, направленного по оси аг,

М 2

 

 

а2 — а').

?3 =

2

ба, ( ° 1 — a i>

A ±A

 

3. Цилиндрическая ортотропная оболочка

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

Уравнения равновесия однослойной ортотропной цилиндрической оболочки в комплексной форме получим из общих соотношений (2 , 2 0 , 18), (2 , 20, 19), полагая в них

 

 

 

 

 

 

а г =

1 =

x/R , a2 =

ср = y/R:

 

 

 

 

 

 

 

 

 

д%

 

 

д<?

 

 

 

 

дё

+

^ г [ (

'1 + -

Л

')Т2 +

10

?!

 

 

д1

 

 

 

 

 

2Ьа

)

 

2Ь2

26*

 

 

1

 

 

Г

<э*г2

 

~

^ 1

 

я 1 д*7\

 

 

 

2

 

26*

 

L

ч* - +

0

5ф*

+

Ч V

+

 

 

 

 

 

 

 

 

'

a*fi

i

<з*т2

\~\

 

(3,2,20)

 

 

 

 

 

 

+ в(к

5|*

1

аФ*

JJ] = <7зД.

где

 

 

 

 

 

 

_ § L х

 

=

 

 

<

 

П

уА ' - ,* 8 S= h

 

e =

 

4G

 

Ei

E± 9

 

 

 

 

 

Ei

 

 

4G

 

Здесь 2Ь* =

V

1 2

0

( 1

v,v2) , R — радиус цилиндра.

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

Приведем систему (3, 2 , 20) к одному разрешающему уравнению, считая 9 i = <72= 0 . Первое из уравнений (3, 2, 20) удовлетворяется тож­

дественно, если положить

 

\ = — ( 1 + - £ - )

&F

 

 

 

асра ’

 

 

 

1

 

Да

\

26аба

)

 

(3,2,21)

 

 

 

 

 

 

Л \

*F

 

 

 

 

 

 

 

 

 

 

 

 

'= —- ( 1 + - О - }

 

 

 

 

 

 

 

 

Яа

V

 

2Ьа2ЬаJ 3£3<р

 

 

 

Второму уравнению

можно

удовлетворить

с

погрешностью

(А//?)2 по

сравнению с единицей, используя

(3, 2, 21) и следующее выражение для

Т2:

 

 

 

 

 

 

 

 

 

 

 

 

Т * = - - 4 -

d2F

 

/0

d*F

 

ie

d2F

(3,2,22)

\ д%

 

2Ь*

аср2

2b2

д%2

 

Я2

 

 

Выбранная таким образом комплексная функция F с такой же степенью

точности будет удовлетворять уравнению

 

 

 

 

 

+ 2 Х ^ -

+ в ™

- + 12Ь>

* F

 

0 ^ +

 

 

д|аафа

 

 

аф4

 

 

д£а

 

Эфа

 

+ е ( г - ^ - +

а2?

=

i2b2R3q3.

(3,2,23)

 

ч

а|а5фа

 

 

aga

 

 

 

 

 

 

Соотношения упругости могут быть записаны в комплексной форме:

М, =

ic (Т2 -

v2f ,),

М2=

£с0 (fx -

Vlf ,),

(3,2,24)

Я = —ic [bS +

(е + v,) 2]

с = tf/26a.

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

7\ =

7\ — £рк2,

Mi =

 

+ tpe2,

 

 

(i?2 )

 

 

 

 

 

 

 

 

 

 

 

2 = S + £рт,

Я =

Я — ip

,

 

(3,2,25)

 

 

р =

 

Ла / а

д

 

 

 

 

 

 

 

/12 (1 -

Vlv2) ’

 

 

 

 

 

 

 

 

 

 

_

1

а2ш

 

 

_

 

1

/

а2^

а» \

 

х *

^

' а|а

Хг

 

д а д

аФа

~ "аф ”)'

 

 

х = _1

/

а2ш

 

аи \

 

 

 

 

 

 

^~ V dldtp

~ ~ д Г )'

 

 

 

 

1

ди

 

 

 

1

 

dv

 

w

 

8l ~ ~R "

~df’

^ = ~R "д?

+ ~R '

 

 

 

СО

 

 

 

ао

 

ам

\

 

 

 

 

 

~

к

-

а |”

ь_аф J ‘

 

 

 

 

 

 

 

 

 

Пусть на свободно опертую замкнутую цилиндрическую оболочку в точ­ ке £= !о, <р = 0 действует нагрузка

Чз = Д2 б2я(| — io, ф), (3,2,26)

где

( I — £0> Ф) = £

6 (£ — | 0, ср — 2 я т ) .

(3,2,27)

Граничные условия будут удовлетворены, если комплексную функ­ цию F представить в виде двойного тригонометрического ряда:

 

 

2

£ fV > sln/mScos/iq>.

 

(3,2,28.

 

 

ш=1 л=О

 

 

 

 

 

 

 

Можно разложить в аналогичный ряд и периодическую 6-функцию,

используя свойства

(3, 2, 16)

и (3, 2, 18):

 

 

 

 

 

^2я (I — £0. Ф) =

£

£

(2 — ЭД Ятпsin /micos /ир.

(3,2,29)

Здесь

 

 

т—1л=0

 

 

 

 

 

 

п

_ siri /m§0

 

t

l

 

 

 

 

»

»

 

 

 

Чтп

 

n

«.

5l

 

 

 

 

 

 

2nil

 

 

Я

 

 

 

 

 

 

1,

n = 0, .

_

nut

 

 

 

 

 

n

 

~

 

Imn

 

~

 

 

 

 

 

0,

пфО,

 

 

ii

 

 

l — длина цилиндра. Подставляя

(3, 2, 26) — (3, 2, 29) в

уравнение-

(3, 2, 23), получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

0-тг\Рт п

bmnFmn — iCmn,

 

 

где

 

ЬщпРтп “Ь ^тп^тп ~

^ОтП,

 

 

 

 

amn =

а<1> — ia<?>

 

 

 

 

 

 

т

 

 

 

 

 

тп

 

 

т п

 

>

 

 

Ьтп =

е (2л2 -

1)Д ;

Cmn = 2b2RP (2 -

fig) qmn)

 

Omi = fm +

 

 

+

@rt2 (n2-

1); a{m =

2b2fm.

 

Из решения системы (3, 2, 29) следует

 

 

 

 

 

 

 

 

&mn

 

 

 

 

 

 

 

(3.2,30)

 

Д »

=

« ! > * +

> * - ч » .

 

Соотношения (3, 2, 24) и (3, 2, 26) дают возможность определить выражения для усилий и кривизн:

 

Tl

 

Тг =

[

Re

~ ±

 

1

 

* т

 

 

R2 Re*3*

---------1шА Ч ,

 

26*

 

<Эф*

J

0

►—1 3

*

1

е

 

 

---- г 1т д12

’ ’*” 26*’

аср2

 

дI2

 

 

 

 

(3,2,31)

а2? - A

i m

. **

1.

 

2Ь2

 

<?1<?ф J

 

 

[I m

a»F

e

 

d2F

 

1

-

 

2b2

Re -d<pa

 

 

1

 

d*F

X

Re

а » ?

 

 

Г I m

 

- +

 

д(р2

 

цЯ* L

 

acp*

 

 

 

 

d*F

X

 

d2F

(3,2,31)

X=

 

Re-

- 1 .

■ [ i m -

 

~2№

d m

рЯ2

 

dtd<(

 

j

Используя (3, 2, 30) и (3, 2, 31), находим прогиб

 

ОО ОШ

 

W = ___ 12Р**(1-У1У*).. у у (2 - 6оя )

sin/w6 cos пер, (3.2,32)

Exh3

i-J

 

 

m= 1 n=0

 

где

«тл (*) = almn Ч- 0 • п2, Ьтп(*) = 2еп2/ т .

Если сила Р приложена в середине пролета, то go = £i/2. Введем безраз­ мерную величину, связанную с максимальным прогибом:

wQ=

12 (1 — viv2) у3

£

2

(2 -e g )» ,

 

 

m=l,3,5 л=0,1...

 

 

^0

Е^тах^Я

 

(3,2,33)

 

P

 

 

 

 

 

Результаты вычислений по

формуле

(3,

2, 33) представлены на

рис. 3.11 и 3.121.

В качестве конкретных ортотропных материалов брались стеклопла­ стики КАСТ-В, ФН, СК-9Ф, СВАМ (1:1). Волокна основы располагались либо в продольном направлении, либо в поперечном. Для получения ре­ зультатов с более широким диапазоном изменения отношения Е^Е\ вво­ дился специальный модуль сдвига G0, равный

0 УЁ&

2(1+ / у л )

При этом в исходных уравнениях исчезали члены, содержащие комплекс­ но-сопряженную функцию. Вычисления показали, что при параметрах у>100, £i>3 введение G0 при определении w 0 дает погрешность не более 5%. Кривые на рис. 3.11 и 3.12 свидетельствуют о значительном влиянии ортотропии материала на прогибы оболочки. При фиксированных зна­ чениях у прогибы оболочки получаются тем меньше, чем больше отно­ шение Е2/Е1. Отсюда следует, что для замкнутых цилиндрических оболо­ чек любой длины выгоднее иметь большую жесткость в поперечном на­ правлении [29, 40].

1 Эти результаты опубликованы в работе [40].

Для изотропной цилиндрической оболочки система уравнений (3, 2, 20) сводится к известному уравнению В. В. Новожилова [51]

Д*У+ —

+ Й6* — = «2b]R (bq3+ ^ ~

(3,2,34)

<Э<р2

* д12

 

Функция Т (|, ф) связана определенными дифференциальными соотноше­

ниями с комплексными усилиями Ти Т2 и S. Рассмотрим конкретную задачу о действии на длинную замкнутую цилиндрическую оболочку двух равных по величине и противоположных по направлению радиальных

сосредоточенных сил, приложенных по концам одного диаметра (рис. 3.13). Для решения этой задачи применялись различные подходы [19, 32, 34, 48_50]. Мы будем придерживаться метода, изложенного в

работе [19].

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

такие компоненты:

 

 

qi = q2 = 0,

=

(3,2,35)

ГДе

 

 

5я (£ _ £ ', ф)=

2 6(Е — |',ф — тя).

(3,2,36)

т = — оо

Используя свойство (3, 2, 18), периодическую дельта-функцию (3, 2, 36) можно представить в виде

6я (Е-!',<Р) = 6( £ - П £ 6 ( ф - т я ) .

т = —оо

Заменяя бесконечную сумму, стоящую в правой части этого равен­ ства, эквивалентным тригонометрическим рядом [52]

 

оо

 

v

« ( » - » » > = т ( т + Е

с о в 2 т ф ) ’

10 п. М. Огибалов, М. А. Колтунов

будем иметь

6я(|-|',< р ) = ^ - £

( 1 ------

1-) ^ (S— £') cos т Ф, (3,2,37).

т = 0,2,...

где б? — символ Кронекера.

Уравнение (3, 2 , 34)

с учетом (3, 2 , 35) приобретает вид

 

 

Я(Г) =

- ^ ^ - Д

( б я),

(3,2,38)

где

 

 

R

 

 

 

 

 

 

 

 

 

W = A2 +

- ^ - + t2 ^ — ,

26* = ^

- / 1 2 (1

— v2)

 

aq>a

*

л

4

'

Рассмотрим сначала уравнение

 

 

 

 

 

М(Ф) =

6Я( £ - £ \Ф ) .

 

(3,2,39)

Нетрудно видеть, что если функция Ф(£, ср; |', 0) найдена, то решение

уравнения (3, 2, 38), т. е. основная комплексная

функция Т, сразу опре­

делится из соотношения

 

 

 

 

f = ^ j 2 ^ P _ д ( ф )

 

(3,2,40)

Уравнение (3, 2 , 39) будем интегрировать методом разделения пере­

менных, представляя искомую функцию Ф(£, ср,

0) в форме тригоно­

метрического ряда

оо

 

 

 

 

 

 

 

Ф(Е,Ч>,6 ' , 0 ) = —

V

Фт (£, £') cos mcp.

(3,2,41)

П

m =0 ,2 ,...

 

 

 

 

 

 

 

Такая запись продиктована разложением

(3, 2, 37). Подставляя

(3, 2,.37)

и (3, 2, 41) в уравнение (3, 2, 39)

и приравнивая коэффициенты при оди­

наковых косинусах, получаем уравнения для определения коэффициентов ряда (3, 2, 41)

_ 2 2 - / 6 а) + т * ( т 2 _ 1) ф т = (1 - ^ б (I - 1 '),

Это уравнение будем интегрировать с помощью преобразования Фурье

& <Фт } = ]

Фт (|, Г ) e r ^ - l ^ d H ,

(3,2,43)

00

 

учитывая при этом, что (52]

&{ Ь } = \.

Применяя преобразование (3, 2, 43) к уравнению (3, 2, 42) и пере­ ходя к интесралу обращения, получаем

 

2 — б?'

п

 

£ )Рdp

 

(3,2,44)

Ф т а д ' ) =

— —

С

 

 

 

4*

J

рП*4 -): 2 2 ib2) р2 + т 2 2 — ,1)

Обозначим подынтегральную функцию через f(p).

Эта

функция пред­

ставляет собой отношение двух конечных функций

 

 

 

 

 

flip)

 

 

(3,2,45)

 

 

f ( p ) =

 

 

 

 

Qm if)

 

 

где

 

 

 

 

 

 

 

fi ( p )

=

exp { i $ —

l ' ) p ) ,

 

 

Qm (p) = p4 + 2 (m2 ib2) p2 + m?(m2 — 1).

 

Полином Qm(p) имеет четыре комплексных корня

pjm — P/+2,m = (—

— l)'+1aJm-f ib.m (j =

1, 2).

Выражения

для ajm

и

bjm выписаны в

работе (19]. Согласно теории вычетов (53] интеграл (3, 2,

44) равен

^ f i P ) d p

2 n i (Вычр=р1т / (р) +

Вычр=р1т/ (р)],

£ > ! ' ,

2nt [Вычр=рзт f (р) +

Вычр=р4т / (р)],

| < I '.

 

В точках P = P JTO= —pj+2,m на плоскости комплексного переменного р подинтегральная функция (3, 2 , 45) будет иметь полюсы первого поряд­

ка. В этом случае ее вычеты будут определяться по формуле (53]

BbI4P=P*m f (Р) =

h

(Pbn)

’ k== 1’ 2 ’ 3 ’ 4 '

(3 )2 ,4 7 )

Подставляя (3, 2 , 47)

в

(3,

2 , 46)

и вновь возвращаясь к

интегралу

(3, 2, 44), окончательно можем записать

 

 

ф . ( 6 .

! ' ) =

- 4 =

^ - е х р ) - » . ( ! - O

I E - 1 ' l l .

 

 

 

ехР {ipjmII — S' 1}

- 2 , 4,

(3,2,48)

Ф .<6 , 1 ') = 4 - 2

 

Pim iPfm+ m2 — ib3)

 

 

/ = 1

 

 

 

 

 

 

 

Таким образом, решение уравнения (3, 2 , 39) построено в форме ряда (3, 2 , 41) с коэффициентами (3, 2 , 48) . Величины pjm затабулированы в

работе (19].-

Из соотношения (3, 2, 40) находим теперь основную комплексную функцию Т:

т а , ф)

iAblP

2 р фт

COS /Лф.

(3,2,49)

nR

 

m=0,2,...

 

 

 

 

 

 

Зная Г(£, ф), по известным формулам [51] можно определить усилия, моменты и перемещения. Остановимся на определении прогиба. Для этого используем зависимость

 

 

d2w

---- A(Ref),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

ю (£. Ф) =

 

A (Re Г) <^,<Ц +

иь (ф) + Щ

(ф).

(3,2,50)

Последние два члена

могут быть отброшены в силу условий «на беско­

нечности». Разделив

вещественную

и мнимую части в выражении (3, 2,

49) и-подставив результаты в (3, 2, 50), получим

 

 

 

 

w а , ф) =

4PR / 3

(1 — у2)

2

wm(£) cos ЩФ,

(3,2,51)

 

 

nh?E

 

m = 0,2,...

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w° &

=

O0* ^ ° s ^

+ sln

exP ( ~ bD>

 

 

2

 

 

 

 

 

 

 

 

 

 

 

(I) =

 

0 0 8

~

W)m Sin ajm D exp (— bjml),

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l >

o,

m = 2, 4,

 

 

 

 

Координата V положена равной п у л ю .

 

 

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

 

w'jm и

wim

приведены

 

 

 

в работе

[31].

 

вычислений по формуле

 

 

 

 

Результаты

 

 

 

(3, 2, 51)

представлены на рис. 3.14. Кри­

 

 

 

вые на этом графике характеризуют

 

 

 

изменение величины прогиба в различ­

 

 

 

ных точках срединной поверхности, об­

 

 

 

условленное

изменением

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

 

 

 

толщины

оболочки.

Сосредоточенные

 

 

 

силы приложены

по концам

вертикаль­

 

 

 

ного диаметра

направляющей

окружно­

 

 

 

сти (см. рис. 3.13). В этих точках имеет

 

 

 

место

максимальный

 

положительный

 

 

 

прогиб (верхняя кривая на рис. 3.14).

 

 

 

Наибольший

отрицательный прогиб до­

стигается в крайних точках горизонтального диаметра

 

(нижняя кривая

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

4.Асимптотические формулы

Впредыдущем пункте решения задач о действии сосредоточенных сил ьна цилиндрическую оболочку были получены в форме одинарных или двойных тригонометрических рядов. Для вычисления прогибов такие решения вполне приемлемы. Однако при вычислении усилий и моментов приходится сталкиваться с расходимостью рядов в точке приложения сосредоточенной нагрузки. Выход из этого положения лежит на пути вы­ деления особенности и установления так называемых асимптотических формул, которые характеризуют поведение усилий и моментов в окрест­ ности особой точки. В. М. Даревский [21], а затем А. С. Христенко [22—24] для вывода асимптотических формул использовали прием рас­ щепления тригонометрического ряда на два ряда с ограниченной и не­ ограниченной суммами, причем это расщепление производилось так, что сумму второго ряда удавалось представить в видё конечной комбина­ ции элементарных функций. Мы изложим здесь другой способ вывода асимптотических формул, развитый для ортотропных цилиндрических оболочек в работах Ю. П. Артюхина [54—55], а для изотропных оболочек

вработе Ю. П. Жигалко [38].

Рассмотрим уравнение

L(w) = 6(х, у),

(3,2,52)

где ь ( ---- , -----\ — линейный однородный дифференциальным Ъператор

\ дх

ду /

o n

эллиптического типа с постоянными

коэффициентами порядка 2 т ^ 2 .

Любое частное решение уравнения (3, 2, 52) называется фундаменталь­

ным.

Фундаментальное решение может быть определено лишь с точ­ ностью до слагаемого, являющегося решением однородного уравнения. Решение уравнения (3, 2, 52) можно найти методом разложения 6-функ­ ции на плоские волны. Это решение имеет вид [52]

(3,2,53)

где

т] = х cos Р + у sin Р — функция плоской волны.

С помощью подстановки

ц = р cos (Р — а),

р2 = х2 + у2\ а = arc tg

интеграл (3, 2, 53) преобразуем к виду

о

Второе слагаемое в правой части выражения (3, 2, 54) может быть отброшено, поскольку оно является регулярным решением уравнения (3, 2 , 52). Остается вычислить интеграл

I cos2* 1 2 (Р — a) dp

L (cosp, sin р)

Замена переменной / = ctgp переводит его в следующий интеграл:

 

2т - 2

dt

2

(/ cos а -ф- sin а)!

 

 

 

 

 

(',+ 1)”‘ (7 Т О ? '

/ г Ы

Этот интеграл можно вычислить с помощью теории вычетов и получить формулу для фундаментального решения [54]

 

.•2т —2 1

2

 

(t cos а ф

sin а)2 т —2

w (р, а) =

Iр

Ш р

Выч

/

t

1 \ (3,2,55)

л (2т — 2)!

 

\(t*

\)mL (

у

 

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

Вернемся к задаче о действии сосредоточенной нагрузки на ортотропную цилин­ дрическую оболочку. Допуская равенство G=G0 и исходя из того, что наивысшая осо­

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

д*Р

d*F

 

дI*

4 2 в '/г д£*а<р* - 4 ^ ~ - = ^ я р ба,ц>)-

(3,2,56)

По формуле (3, 2, 55) получаем

-

ib2RP

 

(3,2,57)

 

(в1^ cos2 а 4* sin2 а) pa In р э

 

4 я 0 3^4

 

 

 

р2 —

ф 2

 

Подставляя фундаментальное решение (3, 2, 57) в соотношения (3, 2, 31), полу­ чаем следующие асимптотические формулы для прогибов, усилий и моментов:

PR2 (б*/* cos2 а ф

sin2 а)

 

£i^3

8я£*цв,/4

Р п Р.

Dn

12 (1 _ ViVj) *

г р =

l ^ e V

T (0 ,/* sin* а +

cosa а>1п Р.

 

4nR6

 

 

 

Т„ =

РВ'и .

 

 

R

(в1/* sin* а -ф- cos* а) р* in р ,

(X

 

 

 

 

 

Р (1 _ 0 ’/*)

 

(3,2,58)

^

= —--------------sin 2а in р,

 

4 я Я 0 /‘

 

 

Р

[(в'/* 4- уг) cos*a 4r(l^rVa®1^*) sin*а] in р,

Мр=