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

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

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

На основании выражений (4.66) и (4.67) выводим формулы для от личных от нуля компонентов докритического напряженного состояния

 

а Г = - р " > ;

а.,

« М - М Й * .

 

 

а22

а11а22

а?2

 

 

 

 

 

 

 

р

= _________ -

Ра2ДМ М

~ Mft?)_______

f

(4.68)

 

уагг [а^Щ / (а(,^)2] + (1

— у) а $ (аиаг2 — afo

 

(1)_____________ - f t 4 . 1 a l W

~ ( ^ ) , l__________

 

 

 

Ya„ [аЦ>в$ -

(oii»)«J + (1 - Y) 4a («..««~ aU

 

Решения

уравнений (2.8)

(второй вариант теории

малых

началь*

пых деформаций] и (2.12) [третий вариант теории малых начальных деформаций] для статической плоской задачи сводится к решению уравнений (4.13), (4.14). Для внутреннего слоя в выражениях (4.13)...

...(4.15) все величины следует отмечать индексом (1) вверху.

Граничные условия при хг — const запишем

в таком

виде:

 

СТ22 1*5= ±V>! =

a 22J k =

±Vft;

®12 U ,=±V h =

°42 U =±V ftI

 

uxl*j=±vh =

u\ *

 

и2 k —±yb = H2IJ k=±v/il

(4.69)

CT22 |r,= ±.1

= Oi

 

=

0.

 

 

Исследуем общую форму

потери

устойчивости трехслойной

плас­

тинки.

 

 

 

 

 

 

 

Решение уравнения (4.13) при условии (4.68) выберем в виде

X = {A Xch — £iXt +

А2ch — С2*г +

А3sh — SiX2 -f-

 

 

+ Л4 sh - у S2x2) sin - у

хц

 

 

 

/

ch j - &"xz +

 

 

\

 

(4.70)

ХП) =

0 2 ch -J-

 

sin f

*г.

 

Заметим, что при таком выборе решений на торцах пластины

=* 0, /

обеспечиваем выполнение в интегральном смысле граничных условий шарнирного опирания (4.18).

Подставляя решения (4.70) в граничные условия (4.69), выводим трансцендентное характеристическое уравнение для определения кри­ тической нагрузки в виде (4.7), где элементы с,у определяем по форму­ лам:

Си — 1^22^12^-- аП (д11--- ®Р) +

а12 (а12

Gia)] SI shaSl*,

C\i ~ 1^22^12^2---Q22 (QU

®P)

ai2 (flu +

G12)I S2 sh ttSa!

C19 5=3 [а22^12^1 -- a22 (aii

®P)

ai2 (ai2 +

&lt)l Si Ch Д^1»

Си = MhfivA — an (a,, — bp) +

aH (aa +

G13)] t2 ch a£2;

cls = c10 = 0; c21=

Gn (au — 6p + a12Si) ch a ^ ;

(4.71)

toi

' Сгг = Gn (а» -

бр + о,2$ ch о£2; с23 = Gn (а,, -

бр +

sh

см =

G,2 (Qii — бр + а 12Сг) sh а£2; с25 =

ся» =

 

cal =

Ia22G,2d — Q23 (аГ — бр) +

а12 (а12 +

^ г)! S ish

 

с32 =

(a2,Gn ^ — а22 (а„ — бр) +

а1г (о,2 +

G,a)l Сгsh a vS?.

 

Cga= |a22Gl2d — a22(Qu — бр) + «и(aw+

 

& ch av^ ’

 

Рм = |O22G,2^ — Оаг (au — бр) +

а 12 (а12 +

^ 12)I С?ch а 7£з*

 

см -

-

f u M (S',0)2 -

<& (оГ/ -

6р">) +

 

аЙ («Й +

ОЙ» Ci*’ sh а ^ ' ’

cse =

-

1сС Ж (Й")1 -

atf (а!1,1- бр(,)) +

 

оЙ (аЙ +

0Й>1 Й" ^

«Т й'1;

 

 

 

 

c4J =

С,2 (Оц — бр -J- Q12S1) ch a Vd»

 

 

 

 

 

 

 

 

с42 = Gu (а„ — бр +

о12$

ch avS2;

 

 

 

 

 

 

 

 

с« =

G12 (а„ — бр +

а12$

sh ау£,;

 

 

 

 

 

 

 

 

с« =

G,2(flj, — бр +

аИЦ) sh ау£2',

 

 

 

 

 

 

<Ъ -

-

ОЙ* |а й -

бр"» +

оЙ (d")2l ch aT d,>;

 

 

 

 

c40 =

— G $ | O (II) — 6p(l) +

eff (£2 Vl chavCsK.

 

 

£51 =

(G,2SI -

au +

бр) ch avd;

c62 =

 

(G12^ -

a,i +

бр) ch avCrf

4>з =

(G,2SI -

a„ +

бр) sh ау&;

c„ =

 

(G,2C> -

a u +

бр) sh ау£*

 

 

 

C65 =

-

|Gi2 ( $ V -

all1+

 

6p("| Ch oTti0;

 

 

 

 

 

cse =

— |G}y (d’’)2 — oiV +

6p“)| ch avCs

 

 

 

 

 

=

(°i2+

Cl2) £, sh ау£,;

ce3 =

 

(a12 -f- G,2) £2 sh ауС2»

 

 

^вз =

(«и +

G,2) £, ch ayS,;

<V =

(«12 +

G]2) £2 ch aySaJ

 

 

 

 

 

с05 = _ ( о Й

+

ОЙ) d ,,shaT;l,);

 

 

 

 

 

 

 

cee = — la i2 + G12) £2*sh ау$*'.

 

a — ^ y - .

 

 

Для

тонкостенных

трехслойных пластин

(a =

nhll <

1),

исполь­

зуя приближенные формулы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh

» a d

+ - i . (ad)3;

ch ad- ®

1

+

4 “ В Д 2-

<4 7 2 )

с точностью до a* из уравнений

(4.7),

(4.71) выводим формулу

 

 

 

 

 

 

 

ай>аю ~

(а1г)2

 

 

(1 — <у)Э

ОцД22~

 

 

 

 

4

 

 

вИ»

 

 

 

 

 

3

 

 

о2а Ц. (4. 73)

В технических постоянных (4.25) формула (4.73) принимает вид

 

Р „ = 4 - а ‘[у>---- - J L ------ + (1 — V)3------1

^

1 .

(4-70

102

Для изотропных слоев

Е<"

- + п

(4.75)

- — * " * [* I - ( v {l)p

Выражения (4.73)...(4.75) совпадают с формулами для критической нагрузки, вычисленной с привлечением гипотезы Кирхгофа — Лява для всего пакета, что свидетельствует об асимптотической точности гипотезы в теории устойчивости длинных трехслойных пластин.

§ 32. Области применимости прикладных теорий в задачах устойчивости стержней и пластин из материалов с низкой сдвиговой жесткостью. Сравнение с экспериментальными данными

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

Существующие сравнения результатов по устойчивости стержней и пластинок на основании теории Эйлера — Бернулли — Кирхгофа и уточненных теорий типа Тимошенко не являются достаточными, так как в этих случаях приближенные решения сравниваются с другими, также приближенными решениями.

Рассмотрим сначала устойчивость бесконечно длинной пластинки толщиной 2/г (ось Ох-л) и шириной I (ось Охх) при цилиндрическом из­ гибе. Будем считать, что пластинка изготовлена из ориентированного композиционного материала и сжата вдоль оси 0хуусилиями интенсив­ ности р. Исследование проведем на основании кинематической гипо­ тезы Тимошенко (§ 13).

Линеаризированные уравнения равновесия бесконечно длинной пластинки получаем из системы уравнений (1.161):

Q x„ x, — pw,Xlxt = 0;

(4.76)

MXliXl — QXt = 0.

Из выражений (1.160) с учетом соотношений (1.158) выводим формулы для изгибающего момента и перерезывающей силы:

М х, = а ц ф ^ . х ,; Qx, = 2 hGls (фх, + к; *,). (4 .7 7 )

Используя соотношения (4.76) и (4.77), приходим к следующей си­ стеме уравнений:

2hGls (фх,.х, + о».*,*,) — pwlXtXl — 0;

(4.78)

4 “ вцЛЧ*,.*,*, — Gia (Ф*. + и».*,) = 0.

Граничные условия шарнирного опирания при х, = 0, / будут выпол­ нены, если принять

=Атsin -

фх, = Втcos -

юз

Эти выражения подставим в систему уравнений (4.78). Тогда

 

+

 

 

A j > ^ ) ‘ s i n ^ = 0;

 

 

 

 

 

 

(4.79)

1 _ . 2 ( тл У О „ „ /ЛЯХ, . г

(л М*1

mn*l

I О

тпХ1\ Л

[~Г1

Втcos —J-L +

Gl3

^Д, - J- COS x x

+ 5 mcos —J—J = 0.

Тривиальное решение системы уравнений (4.79) приводит к значе­

ниям Ат =

Вт = 0 или

до = ip*, =

0. Отбрасывая это решение, не­

обходимо положить:

 

 

 

 

 

Р— 2G13ft — 2G13/i

 

 

 

 

- о »

- с - - 4 - « и ( ^ )

= 0

 

 

или

 

 

 

 

 

 

X anh*Gi3(—

) — Р [^1э+ х

йц/i2 ( х / 1] — 0.

Критическое значение сжимающих усилий представим в виде

 

4 -

M

- f f

 

 

 

 

 

 

°13+ X а,,аг

 

Еха?

 

Eva?

н - 3Gis(l — v13v3l)

1+3G13(1— v,3v3l)

 

 

 

 

 

(4.80)

n

2£j/i3

 

 

жесткость; рм — критиче­

где Ox, =

X (i — v13v3I)-----цилиндрическая

 

 

ское значение сжимающих усилий по теории

 

 

Кирхгофа — Лява;

 

Если в

формулу для

перерезывающей

силы

ввести множителем

коэффициент сдвига к, то формула для критического сжимающего уси­

лия

(4.80)

примет вид

 

 

 

 

 

 

 

 

 

 

р« ' = --------------Pjb

(4-81>

При

к = 5/6 формула (4.81) совпадает с полученной в работе [791 на

основе уточненной теории пластин [139]. На рис.

11 представлены гра­

фики

изменения

величины р* = рКр/рзл,

вычисленные по

формуле

(4.81) [штриховые линии], в зависимости от параметра а =

nh/l при

k =

5/6. Кривым 1, 2, ..., 8 соответствуют значения параметра EJG,3=

=

6,

8,

10,

20,

30,

50,

70,

100.

 

 

 

 

Значения критических нагрузок, подсчитанные по формуле (4.81)

при

к — 5/6,

сравнивались с соответствующими значениями крити­

ческих нагрузок, полученными при решении характеристического трансцендентного уравнения (4.19) для различных механических ха­

104

рактеристик, которые приведены в табл. 5. В результате получено, что относительная погрешность не превышает 1...2% [79].

Формулу для критического сжимающего усилия сжатого вдоль оси 0х3 усилиями интенсивности р стержня кругового поперечного сечения, изготовленного из ориентированного композиционного мате­ риала, легко получить из формул (4.80) или (4.81), заменив DXf на E3J

и 2kGX3h на kG13nR 2

(J — момент инерции поперечного

сечения

стер­

жня):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рчр

—_

 

 

 

 

(4.82)

 

 

 

 

 

 

 

 

Gian/?2

 

 

 

где рэл

имеет вид

(4.26)

[эйлерова

критическая

нагрузка].

р* =

На

рис.

3...4

представлены

графики изменения величины

= Ркр/Рэл, полученные

по формуле

(4.82) [штриховые линии], в за­

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

= nRtl при k = 5/6 и различных отноше-

 

и/к

6 15/0 = 2; Е/Ег=

 

 

3;; Е/Е,=ш3

Таблица 7

 

 

2

С„/0 =

CIS/G= 4; Е/Ег= 4

 

 

 

(1 О'’) юо%

 

(1 -

р•) 100%

U - P '’>100%

 

 

1/12

20,11

17,45

26,53

23,99

32,00

29,64

 

 

1/10

26.04

23,37

33,58

31,34

39,73

37,80

 

 

1/8

35,67

32,23

44,36

41,59

51,01

48,71

 

ний E3lGl3.

Кривым

1,

2, ..., 7 соответствуют значения

E3/G13 = 6,

10, 20, 30,40,60, 100. Сплошные кривые соответствуют решениям транс­ цендентных уравнений (4.7), (4.10), полученных на основании трех­ мерной линеаризированной теории.

В табл. 7 приведены относительные погрешности, возникающие при определении критических нагрузок по классической теории (первый столбец) и по уточненной статической теории Тимошенко[3,116] для круговой трансверсально-изотропной пластинки.

Из табл. 7 следует, что погрешность классической теории достига­ ет для сравнительно толстых пластин около 50%. Погрешность же уточненной теории не превышает 4% для рассмотренных геометриче­ ских параметров и механических характеристик материала.

На рис. 17 приведены результаты по устойчивости стеклопласти­ кового стержня кругового поперечного сечения, армированного пре­ имущественно в продольном направлении. Сплошная кривая соответ­ ствует решениям уравнений (4.7), (4.10), а экспериментальные дан­ ные даны из работы [241. Как следует из графика, теоретическая кри­ вая хорошо согласуется с экспериментальными данными.

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

1. Значения критических нагрузок, подсчитанных с привлечением классических гипотез Кирхгофа — Лява и плоских сечений для тон­ костенных пластин и длинных стержней, совпадают с первыми членами асимптотических разложений выражений для критических нагрузок,

105

которые вычислены на основании строгой трехмерной липеаризиро- ■

ванной теории.

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

стигать 50% и более.

3. При расчете на устойчивость тонкостенных стеклопластиковых

(£ 3/G < 20) шарнирно-опертых стержней и пластин (яУ?//, ш ^

0,13)

наибольшие погрешности

классических тео­

рий не превышают

10...12%.

Значит,

для

указанных диапазонов изменения парамет­

ров £ 3/G, я R/1 и

со

можно

применять в

инженерных расчетах указанные выше тео­

рии. Однако для более толстых стекло­

пластиковых

(£ 3/G <

20) стержней и пла­§*

стин (0,13 <

я R/1,

со ^

0,42)

погрешности

классических

теорий

 

могут

превышать

40%.

г4. Погрешности уточненных теорий типа Тимошенко для исследо­ ванных задач (в случае однородных докритических деформаций) не превышают 4% для рассмотренных областей изменения механических и геометрических параметров (0,01 < £ / £ 3 < 0,8; я R/1, со < 0,42; £ 3/G ^ 100). Выражения для параметра со для различных типов плас­ тин приведены в § 30. Погрешности прикладных теорий устойчивости, использующих гипотезы плоских сечений и Кирхгофа— Л ява, для жестко защемленных стержней и пластин больше, чем в случае шар­ нирного опирания. Такой же вывод можно сделать и для двумерной теории типа С. П. Тимошенко (рис. 12... 14).

ГЛАВА 5. УСТОЙЧИВОСТЬ ОБОЛОЧЕК

§ 33. Цилиндрические оболочки. Постановка задач и методы решений

Рассмотрим задачу об устойчивости замкнутой кру­ говой цилиндрической оболочки толщиной 2ft, длиной I и радиусом R срединной поверхности. Оболочка сжата по оси усилиями интенсив­ ности р и нагружена внешним боковым равномерным давлением q. Материал ее будем считать цилиндрически-ортотропным.

Исследование проведем на основании трехмерной линеаризиро­ ванной теории при малых докритических деформациях (предполага­ ются малыми по сравнению с единицей удлинения и сдвиги, а начальное состояние определяется по геометрически линейной теории — второй вариант уравнений теории малых докритических деформаций). В этом случае трехмерные линеаризированные уравнения и граничные усло­ вия в лагранжевых координатах имеют вид (2.8) и (2.10), а выражения для компонентов тензора деформаций выражаются обычными форму­ лами (2.11).

106

Соотношения упругости в случае малых деформаций для цилиндри- чески-ор готронного тела представляются в виде (1.142)...(1.144), (4.25).

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

Решения линейных уравнений теории упругости ищем в виде

u? = u(r); и£ = 0; «5 = С % .

(5.1)

Тогда система уравнений Ляме сводится к одному уравнению

- + 4 - ^ н - « --5 -+в- ^ - с о - о.

(5.2)

Граничные условия на боковых поверхностях оболочки записываем следующим образом:

осевое сжатие

Оп |г=я±/1 = 0;

(5.3)

внешнее боковое давление

сV \r=R+n = — q\

= 0,

(5.4)

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

винтегральном смысле

 

 

 

л+л

 

 

 

 

 

 

 

 

 

 

 

a%rdr =

— 2phR.

 

 

(5.5)

 

 

 

R~lI

 

 

 

 

 

 

 

Общее решение уравнения

(5.2) имеет вид:

 

 

 

 

 

 

и, =

С,г* +

С2г~ь +

АС°г,

 

 

 

Составляющие напряжений определяем по формулам:

 

 

 

 

 

-£ -)(-£ -)

- ( * • - - ^ ) ( т г )

] :

* - c [ i

-

- £

) ( i - ) 1"

+ X ( * . - - ? ) ( *

) " ~

]; (5.7)

 

азз = С

[ 1

(fe3

$ . )

 

 

+

(* . -

 

 

] i

Здесь

 

 

a% =

0;

о?з =

0;

003 = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ki

(1 +

e)**1-

(1 -

e)x+l . .

=

(1 -6 )-^ + ' - (1 + 8 ) - * + ' .

U -4- e)2X-— (1 — e)2X

*

 

(1 — e)- ^ —(1+ e)2X

 

 

 

и

~ 1

+

Дгз A (gn + ai«) + °is .

p ____A ..

 

 

3

euA +

 

 

auA (a13

+

a33) + аэз’* ’

107

,

f. ам -Д и>.

A (g„ +

a12) +

a„

 

- д (1 + в)Ч-'

 

 

4

* a u X.— a,a

А (о1Я+

a2S) -f aS3 ’

1

(1 +

e)2* — (1 — e)2*-

a

a « и £ ± £ а _ ;

 

 

 

 

 

 

 

 

 

 

 

3

1

a1IX + 6u

 

 

 

 

 

 

 

 

 

 

 

 

C =

C° (Л (axi +

a12) +

a13J;

С =

С°|Л (a18 +

a„) + a33|.

(5.8)

 

 

- { - * * • - x f l

 

 

 

+ * f * - (I -

 

+

 

+ - A c

 

1(1 +

e),~l -

<> - «У-Ч} (2 M

+

<■»> + aMl E

-

И («и +

Ou) +

a j

[(1 +

e)*+* -

a -

e)*+«I - 2 4 + i ^

*L_

+

 

 

 

 

 

 

 

 

 

 

 

аПЛ“Г в1г

f'-r l

 

+

|A (eu +

an) + a13| [(1 +

e)1-* — (1 -

e )'- x| - Ж -

g“ .~

a|3M ~ ‘ .

 

 

 

 

 

 

 

 

 

 

i — к

апл — а1г j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 5 .9 )

Для трансверсально-изотропной оболочки (ось изотропии совпа­ дает с осью 0ха) имеют место равенства (2.38), (4.30). Поэтому получаем только при осевом сжатии распределение напряжений такое же, как

ив изотропном полом цилиндре (4.1).

Вслучае внешнего бокового равномерного давления даже для изо­

тропной оболочки докритическое напряженное состояние является не­ однородным.

§ 34. Трансверсально-изотропная оболочка

Рассмотрим устойчивость трансверсально-изотропной оболочки при осевом сжатии. Для этого построим решения уравнений (2.10) в ци­ линдрической системе координат. В случае трехмерной задачи для трансверсально-изотропного тела, ось симметрии которого совпадает с осью Ох3, при однородных докритических деформациях решение урав­ нений (2.10) сводится к решению уравнений (2.30). В зависимости от

соотношений между величинами возможны различные случаи пред­ ставления решений уравнений (2.30).

Рассмотрим случай вещественных и положительных I? при Еа ^13Решение можно выбрать в виде (3.26).

Величины для трансверсально-изотропной оболочки [соотноше­ ния упругости имеют вид (2.37)1 при условии (2.38) определяем по формулам:

£г.з =■ С ±

■|^/"^2

(Q331+ Ра°з)(033 (< + о°3)

^

G +_ O33

 

 

^

’» У =

Сп

 

а 11 (азз+ °зз) + О (Д + Р°э) — (а13 -I- Gf

(5.10)

2С =

*

 

 

an(J

 

108

Представим решения уравнений (2.30), отличные от (3.26). При ве­

щественных

и отрицательных £? и

в случае, когда

Сз ф £з, в реше­

нии (3.26) для функции X вместо функций 1„ (г) и Кп (г) следует брать

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

функции Бесселя Jn (г) и функции

Неймана N„ (z).

При действительных $ и комплексных Ц и £з для функции V полу­

чаем (£? > 0) такое же выражение,

как в (3.26), а для функции X со­

гласно выражению

(2.32) находим

 

 

X = [A;,L Re Jn UD + A ll Im Jn (y&) + Al„ Re HP (yt2r) +

 

+

А)пп Im НУ (Y£201 COS л0 COS ух3,

(5.11)

где H[n — функции

Ханкеля первого рода.

 

Для других случаев можно получить аналогичные представления решений.

Из соотношений (2.10) выводим при условии (4.1) формулы для нахождения составляющих поверхностных сил на круговой цилиндри­ ческой поверхности при г = const:

Проекции внешних усилий, приложенных к торцам оболочки, при х3 = 0, / имеют вид (4.3).

Выражения для перемещений в круговой цилиндрической системе координат в случае трансверсально-изотропного тела имеют вид (4.2).

Предполагая, что боковые поверхности ненагружены, из соотно­ шений (2.10), (4.1) получаем граничные условия на цилиндрических поверхностях:

Orr Ir=R±h = 0; а,о |г=я±й = 0; а гз |г=я±/« = 0.

(5.13)

Из выражений (4.2), (4.3) и (3.26) следует, что прих3 = 0, /выпол­ няются следующие граничные условия:

иг|*,=о./ = 0 ; ын|*,=о./ = 0; Р3|*,=о./ = 0-

(5.14)

Таким образом, при х3 = 0, / автоматически выполняются в интеграль­ ном смысле условия шарнирного опирания.

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

у щ / з Т Г = ^ > 10 и з о -

Учитывая представления для перемещений и соотношения (1.145), (2.38), выводим формулы для компонентов тензора напряжений,

выраженных через функции Y

и X:

„ w _ L _ £ _____L J L W - , Г а1Я(0 - р ) а»

Orr = (fln вю)*v т агав

г* д й +

-----5^+G ^ ***» +

аов = »a,i - a12) ~

— flia (°mН~ С)

I /

1

+ _ ? ! _

|v .

g.3+"G

_ Ш

____* ! _

г ^

г

ао2ахя

^

дгдхя

JA>

(7 - £

Т й г ) ¥ +

 

" £ r +

 

+ ( ■ £ & - °“) W + t o - «a) H + ^ k +

+ i s d ] x :

- { (°” ~ t o

 

л - &

 

+

 

 

( - g r +

 

.

1

 

w

.

I

 

 

 

(5.15)

 

+

" Г “5Г +

75-“

"®TJ

 

 

 

 

 

 

 

 

 

O K . - 0 . ,[ ( - f - |r

+ 4

" W

- ------+

 

 

- ■ £ - w - + * - w b - 2 ^ j y x } i

i

flu

G

(

аз_

-|

1__Ё!_ .

i а»

\ l Y i

Г

а1я +

\ дгЧ0 ^

 

г дгдй +

7 Г W

j | Л ) •

Подставляя решения (3.26) в граничные условия (5.13), из условия нетривиальности решения алгебраической системы уравнений (ис­ пользуя при этом рекуррентные соотношения для модифицированных цилиндрических функций) получаем трансцендентное уравнение для определения критической нагрузки:

Здесь

det| a //! =

0 (f, / = 1,2, . . . . 6).

(5.16)

 

 

 

 

 

«и (/.+ ,. О =

/«+. К , (х + е)1

К 1C, (х + «№

а п = a u ( -

К я+,, К„);

« и (/„+ „ /„, а =

/„+ . [£, (х +

«)] +

+ [£s +

Н--- (v j- eii ^1

1£г (х + е)1;

сси = аи (— К п + н K nt £2)•

(5.17)

но