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

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

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

Этот интеграл .выражается в общем случае 'через эллиптические функции I и II рода.

При <7=0 будем иметь гиперболоид вращения. БслиР = — щ(г\ -г\) (где гх — значение радиуса параллельного круга, при котором f (г) = 0), мы .получим тороидальную оболочку. Итак, формула (6, 8, 2 1) заключа­ ет в себе целый класс .поверхностей, .получающийся при различных соот­ ношениях .между давлением и осевой силон, приложенной к .полюсному отверстию.

3.Поверхность, нагруженная осесимметричным давлением

ицентробежными силами

ф[(г) = т ф 2( г ) , Фг(г)=0,

где m — масса единицы длины нити. Из

уравнения (6, '8, 3) d T = — rruahdr.

 

 

 

 

 

Интегрируя по г, получим

 

 

 

 

 

 

T =

- i - m

c o V

+

r 0 ( r 0) ,

 

( 6 , 8

где То — натяжение в точке .полюса нити. Подставляя (6, 8, 22) в (6, 8,

10), будем иметь

 

 

 

 

 

 

------/псо2р2J rein <р =

T Q х

г0.

(6,8,23)

Рассмотрим линии, удовлетворяющие уравнению '(6, 8, 23)

и имеющие

один и тот же полюс г0.при .различных значениях Т0. Из (6, 8, 23) имеем

sta 4. =

r0/ r

[ l _ ^

- p« ].

(6,8,24)

Та.к как sincp<^l, то 'Минимальное значение Т0«равно mcoVo

и линия

(6,

8, 23) совпадает с параллельным кругом радиуса г0. Если

Т 0 » ^ - Р

2,

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

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

Если г0 равно ,радиусу экваториального круга .поверхности, то -се­ мейство линий вырождается в одну линию — экваториальный круг, ко­ торый является одновременно и геодезической линией поверхности для любого значения Г0. Если г0= 0, то семейство линий (6, 8, 23) вырож­ дается в меридиан поверхности при любом значении Т0.

Уравнение (6, 8, 17) в .рассматриваемом случае будет иметь вид

I

 

 

 

W > ) - 1

±

т а зу г2] - 7 > 2а (r0) dr0 = %(г), (6,8,25)

 

а' (г)

где

 

 

 

 

 

 

 

1 (г) = яг j

<7{г) г dr -f rR (г).

 

 

Гх

 

Дифференцируя *(6, 8, 25) по г и введя обозначение

Ж= а ( г ) / / а ' ( г ) - 1 ,

получим

г

ri

 

Р ' ( ') X (г) * Р ( ' ) % ' ( ' )

(6,8,26)

 

Г

 

Пусть

 

Н

(г г0) , =

(6,8,27)

Уравнение (6, 8, 26) может быть 'сведено к интегральному уравне­ нию Вольтерра II рода .при некоторых 'Предположениях относительно функции Я(г, г0). Так, если «диагональ» Н(г, г) не обращается в нуль ни в одной точке интервала от г\ до г и Н(г, Го) .конечна в этом же ин­ тервале, то уравнение (6, в, 26) можно 1привести к уравнению II рода. Легко показать, что оба эти условия 'будут выполнены, если

Т0> mto2r2.

(6,8,28)

И т а к ,

 

Р' (г) %(') + Р (') %' (г) .

(6,8,29)

»

Г

Т0> то>2г2.

Умножая обе части (6, 8, 29) на {z^ rtyi, >интегрируя в пределах от

Г\ до z и применяя .правило Дирихле, получим

r . w d r . r

 

„ г

(6,8,29а)

J

J

PiP

J

Pi

Н* (г, О = J Н (г, г0) 2г dr

9i?

Го

где обозначено

Pi = (г2 —г2)’А.

Условия, наложенные на функцию Н(г, г0), 'Обеспечивают выполнение этих же условий для функции Я* (г, г0). Итак,

j Я* (z, г0) а (г0) drо = J 2(Р-^ -Х-)- dr.

(6,8,296)

Гх

Г0

Дифференцируя 'полученное уравнение по г, имеем

 

 

 

 

 

Я* (zz) a (z) +

Г — а (r.) dr0

dz

2 (Р'х + М dr

(6,8,29B)

 

 

 

 

 

J

дг

 

J

 

pi

 

 

 

 

 

 

 

r1

 

 

 

ri

 

 

 

 

Дифференцируя, имеем

 

 

 

 

 

 

 

 

 

 

 

 

дН*

_

2z

f

Hr (г, г) p

dr,

 

 

 

 

 

 

dz

 

2a — TQ J

p!

 

 

 

 

 

 

 

 

 

 

 

 

 

_a_

f

2(Р'Х + Рх')

d r _

P 2 Г

х ± Ж ± Р

 

 

 

ri

2 (PfX +

P x ') U ,t

 

Л r d r + ^ L .

 

 

dz J

Pi

 

J

 

zpi

 

 

 

 

Z

 

Pi

Тогда окончательно .получим

 

 

 

 

 

 

 

 

 

PIP

 

 

 

 

 

 

 

 

 

 

 

 

 

= Г 2(P"X 4

2PY +

PX") r jr

!

П

2 (P'X -fr

Px’) U r ,

 

 

 

J

 

г р !

 

 

 

 

 

 

 

 

 

Здесь

функция,

стоящая

в правой

части, терпит

разрыв в точке z = ri.

При решении уравнений в каждом конкретном случае необходимо иссле­ довать выражение

З Ф 'х Ф Р х ') U r ,

гPi

в точке Г\.

Уравнение (6, 8, 12а) равновесной'поверхности запишется:

^ P ( r ) r d r + j

V а ,г(г) — 1

 

 

 

(6 ,8 ,1 2 6 )

о'(г)

 

 

 

 

^ ' [ т . - ^ е у - т у ,

 

 

 

 

Г

 

 

 

 

 

nr^p(r)rdr +

у

г = %rz = f (г).

 

Г 2 (г)

 

 

 

хаМ

 

1 + Г

(г)

 

 

 

 

Г { г ) = ± г ( х )

/ ] / ' {

[ г

0-

тРаjЩ 2г*п* -

X2(г),

 

Г

 

 

X (г) d r

 

Z = / ( r ) = ± ^ . Y { [ Т“ -

 

 

 

^

Р

а] 2л2- Т^ о }« а - Х

2

-0 гран!ичи1вая'сь рассмотрением поверхностей, симметричных относи­

тельно ПЛОСКОСТИ Г = Гт а х ; /'(rm ax) = |СО

 

 

 

n = x(rmj

1 1/ [ T . - ^ f p l ^ J r L . - T l r l

 

обозначим

 

 

 

 

 

 

X (r max)

=

К ,;

 

 

 

-

 

 

 

гтах^о

 

 

 

 

тсо*

1

 

 

= сг.

 

2

(Гт а х - Го )^ Г

=

Й1-

 

Тогда

 

 

 

 

 

1 =

Кг

(1 — Пх)* — Cf = /с?.

(*)

 

Значит, параметры К\, С\, £2i, характеризующие нагрузки на поверх­ ность, и соотношения между щах и min отверстиями, не могут 'быть за­ даны произвольно, а связаны соотношением (*). В координатах К, С, У уравнение (*) представляет собой конус с вершиной в точке £2= 1.

При £2i='l Ki=0 и Ci =0. При £2i = 0,/Ci + С? = 1, т. e. все значе­ ния Ki и Ci лежат на окружности радиуса 1.

4.Формы поверхностей

Вработе [35] подробно исследуются возможные равновесные формы при действии постоянного давления и центробежных сил. Если в форму­

ле ('6, 8, 126) ©вести следующие параметры и переменные:

/ 2 •*

2праг\

m /QW3

£2 -

С = sin Ро, / с =

Та

иТо '

где Го— значение экваториального радиуса, Г0— 'натяжение в нити .на экваторе, Ро— угол нити *с меридианом *в .плоскости экватора, то она получит вид

dz

 

 

(1 _ с 2)1'* —

 

ТС(1 — Ла)

 

(6,8,12г)

 

 

 

 

 

 

 

dR

Г

1

12 Са

Г

1

/ 1

 

 

 

Кроме того, из выражения г ---- =

sin ср, формулы

(6, 8, 24)

и шредыду-

щего уравнения имеем

dsi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с_

 

 

 

 

 

 

 

 

R*

 

 

 

 

 

 

 

^

- [ О -

С*)’7’ -

\ К (1 -

V )] У''*

 

 

 

 

 

 

 

 

(6,8,12е)

Здесь ср — центральный угол траектории нити.

Легко видеть, что

(1 —

— «■)

 

 

(6,8,12г)

 

 

у

Дифференциальные уравнения

(6, 8, 12г), (6, 8, 12е) .полностью on-

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

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

у = [(7(2 _

Q2) *3 _ 4 {Q [! +

(Q/2)] +

К [(1 _ С2)1/, _ (ТС/2)]} *2 +

+

4 {(1 + Q/2)2- [(1 -

С2)1/, _

(К/2)]*} х - 4С2]1/»,

(6,8,30)

где x = R 2, в виде произведения. Один корень определяется выбором 2 = 0 в качестве плоскости 'симметрии и равен *1 = 1. Исключив этот корень, получим квадратный трехчлен

(Я 2 — Q2) х2 +

4 I"К (1 — С2)V. + Q — у - (К2— G2) х — 4С2

 

и два другие корня

 

 

 

* . = -

г & - Iя 2 + С2 ( * 2 -

;:

 

*, = -

- ^ - = Г {Н + [Я2 + с 2 (Я2 -

В Д /’};

(6,8,31)

я = К (1 — С2)’- + Q — — (Я2 — й2).

4

Следовательно, (6, 8, 30) можно записать в виде

у = [(Я2 — Q2) (1 —*)(* —х 2)

— *з]'л .

(6,8,32)

Возможны 24 варианта записи у, которые обеспечивают положи­

тельное значение .подкоренного выражения

(12 для /(>£2 и

12 для

/<<й). Однако большинство этих комбинаций не соответствует физиче­ ски реализуемым значениям £2, К и С для вращающихся оболочек при действии давления. Можно сформулировать следующие теоремы, кото­ рые определяют приемлемость данной комбинации корней.

Т е о р е м а

Б. Если К > £2 и *2Х ), то х3-<0 для

всех физически

реализуемых значений К , £2, С.

для всех физиче­

Т е о р е м а

;В. Если К<£2, х3> х 2, х2^>0, то *з>'1

ски реализуемых значений К, £2, С. Кроме того, из начального условия г = 0 при х= 1 вытекает требование, что * и 1 не могут быть .разделены

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

х2 и *3. Используя эти критерии, получим физически реализуемые ком­ бинации:

А /С > й , * * 0 20 < 1 ,

Б / O Q ,

* 3< l O < . t 2,

 

В 7C<Q,

*20 < 1 < * 3,

 

Г / С < й ,

1 < х < х ъ<.х3.

(6,8,33>

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

Из условия осевой симметрии форма поверхности .может быть .полу­ чена .вращением кривой в меридиональной плоскости <p = const вокруг оси z. Наклон меридиана задан уравнением (6, 8, 12в) и запишется а. виде

{[(1 - &)'!* -

(К/2)] + (К/2)х)

dz

 

dx

[(1 — С2)1/* — (К/2) (1 - х)]2}

{ * £ 1 + -J Q (! - х)]2 - С * -

Определяя у выражением (6, '8, 30), введем обозначение

(6,8,34)

X

Тогда уравнение меридиана можно записать в следующем виде:

Z ( x ) = ± { [(К/2) -

(1 -

С2)1/*] /ом -

(К/2) /, w .

Для случаев А а Г форма меридиана определяется выражением

'

{[/С (1 -

*8) - 2 (1 — С2)1/,] F, (фх, k,) -

- *

• (

1

(i MCx)}

(6,8,35)

где .Fi01>i, &i) — эллиптический интеграл

первого

рода, F2(ф^, kx) — эллип-

тичеекий интеграл второго

 

 

 

 

Для .случаев Б к В модуль k будет мнимым, и выражение для ме­ ридиана примет вид

Fi (ФА)1 + К (х2— х3) [F2(А,) — F2(ф2, А,)]},

(6,8,36)

где

 

Исключительным оказывается случай K = Q. Сделав предельный пере­

ход в (6, 8, 31), получим .следующие выражения для корней:

 

lim х2 = С2/[ 1 + (1 — С2)1/»]; Пт х3 = + оо.

(6,8,37)

Рассматривая решение в форме (6, 8, 35), найдем

 

lim [(/С2Й 2) ( 1 — *3)]'/‘ = 2 { Й [ 1 + (1 — С 2) ’/.] } * /.,

lim kx= 0.

К-Ю

K-*Q

Разлагая эллиптические интегралы в ряд и используя полученные выше оценки, запишем решение для .случая /<=Й в виде

Z =

 

X arc sln ( т 5 ^ ) ,Л - т 1(1 “ х) {х - х*)]',г}'

(6,8,38)

где х2 определяется выражением (6, 8, 37). Ори С= 0 (волокна .направ­ лены по меридиану) уравнение (6, 8, 38) дает

Z = {[1/(2й)‘/‘] ( l — -L й)агс sin (1 —*)*/. + ±

[2Q (1 — *)]■/.}.

(6,8,39)

Рассмотрим теперь условие существования точек возврата меридио-

налыной кривой. Обозначим через ха значение

х, при

котором

кривая

имеет нулевой найлон, и учтем, что х2— значение, при

котором произ­

водная обращается в бесконечность. Тогда точка возврата определяется равенством Хо=х2. Беря хо из уравнения (6, 8, 12в), а х2 из уравнения (6, 8, 31), условие х 0= х 2 после возведения в квадрат и группировки чле­ нов примет вид

Кй(1 — Q) + Й2(1 — С2) — (Л* + й2— 4Й) (m ) (1 — С2)1/. = 0.

Выражая из этого трехчлена й, имеем

2 = ________* _____/1

±

r i +

^ ' r ^ S / -

1 - ( / Э Д ( 1 - , < ? ) ' М

|,'-1

 

(1 — С8)1/2

I

 

L

К

\

1 — (2/Ю (1 — С*)'1'- 7J )

Чтобы получить й > 0,

берем нижний знак

и после преобразования ра­

дикала имеем

 

 

 

 

 

 

 

 

 

Й

*

 

Я

 

са________ Т ;* -

А

 

 

( 1 - С 2)1/»

Ц

1 — (2//С) (1 — С2)*^ J

Г

 

или.через

 

 

 

 

 

 

 

 

 

 

 

 

х0 = 1 - 4 о - с * ) * / .,

 

 

 

 

 

 

 

к

 

 

 

 

 

1— Х0 \ \ Х 0 /

 

)

1 — *0

 

 

Найдем выражение для центрального угла.

огибать ось Z при навивке.

Если

параметр

С Ф 0,

волокно будет

Центральный угол ср находится интегрированием уравнения

(6, 8, 12е).

Приняв за -переменную x = R2, преобразуем это уравнение >к виду

dtp =

+

 

 

 

 

С/2

 

 

2)' h

-dx

:{х 1

 

 

 

 

[(1 -< ? )* '.

 

 

 

 

 

 

 

*)

43*

Введя интегралы / ь определяемые уравнением (6, 8, 34), получим

1

 

< р = ± с Г - ^ - = ± CI1 (.V).

(6,8,40)

J dy

X

Интеграл /_i сводится к неполному эллиптическому интегралу третьего рода, который сводится .к табулированным тета-функциям [36]. 'Полный интеграл может быть представлен лямбда-функцией Неймана, которые также протабулированы. Через эллиптические интегралы первого и вто­ рого рода А0выражается следующим образом:

 

( - § - )

Л <>№ •

Ь) =

[Ft (ft) -

F, (A )] F, (ф, k')

+ F , 0k) F2 ( # ' ) .

(6.8,41)

где k'

= (1 А2)* /!.

Для случая А получим

 

 

 

 

 

 

 

Ф (**) =

±

 

яС

 

77

Ло (5l> ^i)>

(6,8,42)

 

 

■ * 2 * 3 ( К 2 -

Q 2 ) ]

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 — * 2

\ Ч г

 

Для случая Б

 

 

 

 

 

 

 

* 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф (*2) =

±

 

пС

 

гту^

^0 №2^2)*

(6,8,42а)

 

 

 

 

 

где

 

 

 

 

[— ***з(Я* — Й*)1 '*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin Ь, =

( - p

^

Y 7,

А, =

 

( i ^ - Y

/2

 

Для случая В

 

V 1 *3 J

"

 

\ *г — *з /

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф (*«) =

±

* - К 2)4 * \

 

2fl(ftt)

 

1

 

“ Т77 Л0 » ] ,

(6,8,426)

где

 

(Q

х3(х3 -х*)'/>

(хЛ ) ‘

 

 

 

 

 

 

 

 

 

 

1 — х2

У

 

 

 

sin £а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

случая Г

 

 

 

 

 

 

 

 

 

 

±

^

г

{

Р

^

.

+

1

^

где

 

 

 

 

 

 

 

 

 

• ' " И

^ Г

А, = ( - ^ 5 т У ' '

 

 

 

 

 

Для частного случая K = Q, неполные интегралы могут быть представле­ ны через элементарные функции. Для случая А, так как в пределе K->-Q

инаходим

^>i

П(Фх. а2, 0) = Г ----^ ---- = ----- L _ _ a rc tg [(1 — а2)1/, tg фД,

J 1 — а2 sin2 ф

( 1 - а 2)'/»

Подставляя «предельные оценки из (6, 8, 37) и др. в (6, 8, 40), най­ дем выражение для центрального угла

(6,8,43)

Эта формула «справедлива при /С = й для всех четырех случаев (6, '8, 33). Полагая здесь х-+х2, «получим ф-^:п;/2, т. е. при изменении радиуса от мак­ симального до «минимального значения волокно для случая K = Q и С Ф О делает четверть «оборота вокруг «о«си г.

Нам остается -построить формулу для длины волокна. Пространст­ венное положение волокна найдено — «оно полностью определяется функ­ циями Z(R) и cp(R). Однако для многих целей необходимо знать длину волокна. Например, для определения характера деформации витой обо­ лочки .при различных комбинациях давления и центробежных сил необ­ ходимо сохранить инвариантность длины волокна и приращения угла

(центральный угол на оборот). Из рис. 6Л(2 видно, что d L

cosp cosa = dZ.

Из уравнений (6, 8, 12в), (6, 8, 12) находим

 

dx

( 1 > 1/2Q) — l/2Qx

(6,8,44)

У

 

где у «определено уравнением

(6, '8, 30). Через интегралы

/ 0 и 1\ имеем

Для случаев А и 'Г, определяемых (6, 8, 33), найдем

 

где

 

 

Для случаев Б и В имеем

 

 

L = + --------------- ------------- гг {[2 £2 (1 xz) Fi 2&2) — ^

*^з) ^2 (Фг^)]}*

- [(К2 — Q2) (х2 — *3)] /’

 

(6,8,47)

 

 

где

Случай /С= Q может быть п-олучен из (6, 8, 46) путем перехода «к пре­ делу. При K-+Q «получим

L = ±

_

j[ , + - I п (1 - *s>] .rcsin(-l

-

(Q [1 -ф- (1 с 2)1^1

 

(6,8,48)

где

Анализ полученных результатов начнем с классификации форм ме­ ридиана. Классификация «осесимметричных оболочек, навитых из «воло­ кон, может быть представлена путем исследования совокупности «реше­ ний (6, 8, 35) —(6, 8, 37), «определяющих форму меридиана. Из этих вы­ ражений следует, «что форма меридиа-на зависит от .трех «параметров К,

Q, С, отределяющих 'соотношение между интенсивностью давления и центробежных сил -и углом соответственно. Будем .полагать К и Q неот­ рицательными, что соответствует внутреннему давлению и |растягивающим силам, действующим в волокнах. Тогда совокупность возможных решений можно разбить на два класса 1периодических -функций: «волно­ образный» (.рис. 6,13) и «петлеобразный» (рис. 6.14). Соответствующие поверхности вращения имеют вид гофрированных труб, или торов.

Эти два класса меридиональных кривых разделены переходной фор­ мой, имеющей точки возврата (рис. 6.13,д), и ограничены выродившейся

Z

 

Z

о (

 

 

 

а

К 1

1а :

1

Z

 

1

1

1

 

: кI

^6

1

1

 

 

Рис. 6.13

Рис. 6.14

периодической кривой, которая в случае волнообразного класса прини­ мает гиперболическую форму (рис. 6.10а), а в едучае петлеобразного класса представляет 'кривую с одной петлей « асимптотическими ветвя­ ми (рис. 6.146).

Каждый класс можно разбить на два подкласса. Кривые, принад­ лежащие волнообразному классу, могут быть «вогнутыми» в 'случае, если при Z = Z0имеет место минимум и R> 1 '(рис. 6Л0в), или «выпуклыми», если на экваторе реализуется максимум и R < \ (рис. 6.13d). -Переходная форма этих подклассов — прямая R — 1, соответствующая круговому ци­ линдру (рис. 6.10с). Кривые, принадлежащие петлеобразному классу, могут быть «прогрессивными» (рис. 6Л4а), когда последовательность значений R= 1 принимается ординатами кривой, начиная от R= 1, Z = 0, в точках, которым соответствует возрастание положительных Z, или «регрессивными», если точкам, где R —Л, соответствует возрастание от­ рицательных Z (рис. 6.14с). Переходной кривой между этими видами оказывается замкнутая петля, которая образует меридиан торовой обо­ лочки (рис. 6,14в). Изучая переходные формы, можно установить пре­ делы изменения К, £2 и С, соответствующие различным классам оболо­ чек. В результате для некоторых поверхностей получим 'условия для К, £2, С. Для гиперболоида необходимо, чтобы К = 0, £2=0, так как гипер­ болическая форма соответствует неограниченно возрастающим напряже­ ниям по сравнению е напряжениями от внутреннего давления или цент­ робежных сил.

Для цилиндра, определяемого постоянным значением радиуса, тре­

буется выполнение равенства х2= 1. Условие будет

 

£2 = С2— /((1 — С2)’/».

(6,8,49)