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

1200

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

Результаты вычислений по формулам (66) и (67) показаны на рис. 24. По оси абсцисс отложен параметр нагрузки

Р2 = К 12(1—v*j

по оси ординат — безразмерная дисперсия перемещения а^(г) в сфе­ рической оболочке. Параметр х0 = k jk ^ принимался равным 0,5, 1 , 2 и 4. Кроме того, на рис. 24 нанесена кривая, соответствующая пре­ дельному случаю х0->- оо (дельта-коррелированному полю начальных

отклонений). Правые ветви кривых соответствуют растягивающим уси­

лиям. При

р2 оо имеем

ф(0) ср0 (0).

Левые ветви соответствуют

сжимающим

усилиям. При

|32 - > —2, т.

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

к его критическому-значению, определяемому по линейной теории, дисперсия перемещений ^i(r) стремится к бесконечности.

§ 1 1.4. Краевые эффекты при докритических деформациях

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

цеских деформациях г

сечения будем обозначать через TWo(k2)• Вычислим эффективное вол­ новое число, соответствующее спектральной плотности TWo(k2)i

 

1/2

J

(^2) ^2

&2е

(71)

J

TwAkt) dk2

О

 

Волновое число k2e характеризует скорость изменения поля w0(г) в на­ правлении касательной к контуру оболочки. Пусть а2 — характерный размер срединной поверхности оболочки в этом направлении. Если выполняется условие

&2е^2 ^ 1 »

(72)

то при рассмотрении краевых эффектов в областях, достаточно удален­ ных от угловых точек, точек приложения сосредоточенных сил и т. п., можно отвлечься от учета граничных условий в направлении xt — const. Ищем решение уравнения (68) в виде канонического представления типа (1.180):

Щ (х1г х2)=

I

5

И70(kv k2) <р (лгх I klt k2) е1к‘ х>dki dk2;

 

ОО

00

(73)

 

 

Xi (xi>*2) =

J

J

W0 (kv k2) ф (xt I klt k2) elk‘ x‘ dkj dk2.

 

— CO

— 00

 

Здесь роль базисных функций выполняют выражения

ф1 = Ф(х1 |й1, k2) е1к‘Хг;

Фа = Ф ixi 1К k2)e‘k,x'.

Подставим разложения (73) в уравнения (68). В результате получим систему обыкновенных дифференциальных уравнений относительно функций <p(xi | ku k2) и Tjjfo I klt k2) (штрихами обозначено дифференци­ рование no *i):

D 1У—2k\ <p' -f ki ф) — Ф" — ~ ф^—()Va ф' — k\ Nt ф) =

 

= - { k \ N l + k lN 2)eiklXl\

(74)

—— (ф1У 2 й|ф -Ь^2 ф) + ф " —

=

Достаточные условия получим, используя метод из работы (151. Полагая у = —г2, представим уравнение (78) так:

yi + c1y3 + c2yi + c3y + ci = 0.

(79)

Для коэффициентов уравнения (79) получаем следующие выраже­ ния:

Ci = 4k\ +

С2 = +

(80)

cs = 4 k l+

D

k i+

- r f ^ - kl;

 

 

 

URi A2

_ (.8

,

t 6 I

Eh

u4

C4= «2

T---— «2 H------«2'

Решения типа краевого эффекта будут существовать, если среди корней уравнения (79) нет положительных действительных корней. Применяя к уравнению (79) известную теорему Декарта, приходим к достаточным условиям существования краевого эффекта

сх> 0; с2> 0; с3> 0; с4> 0.

(81)

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

ную равномерным внутренним или внешним давлением.

Тогда Rt =

= R 3 = R; Nx = N 2 = N. Достаточные условия принимают вид:

N > — 4kiD;

 

 

, / 4

2 k 4

(82)

N > — k lD U + ^

Здесь k * — параметр оболочки, определяемый по одной из фюрмул (58). С другой стороны, в данном случае нетрудно получить явные выраже­ ния для корней уравнения (78):

Из формул (83) вытекают необходимые и достаточные условия су­ ществования краевого эффекта:

N > — 2k\ D,

 

если

| /г2 1^ kt \

, /

Л4 \

если

(84)

k2 D \^ \ .- \ - - j\,

\k2\> k ^ .

Условия (82) и (84) проиллюстрированы на рис. 26. Граница, опре­ деляемая согласно (82), обозначена сплошной линией, а граница, оп­ ределяемая согласно (84), выделена штриховкой. При \k2\ = £* грани­ цы совпадают. Заметим, что соотно­

 

 

 

шение N = —2k^D соответствует до­

 

 

 

стижению критических

напряжений,

 

 

 

определяемых согласно линейной тео­

 

 

 

рии

упругой

устойчивости.

Таким

 

 

 

образом, краевой эффект существует

 

 

 

всюду в докритической

области. За­

 

 

 

метим также, что при

| У | <

2k\D

 

 

 

осциллирующий краевой эффект

[15].

 

 

 

При

N 2k2,D

— краевой

эффект

 

 

 

неосциллирующий.

 

 

 

 

 

 

 

Как видно из

формул (83), в

об­

 

 

 

ласти, достаточно

удаленной

от

гра­

 

 

 

ницы вырождения

краевого эффекта,

 

 

 

действительные

части

характеристи­

 

 

 

ческих корней имеют порядок

волно­

 

 

 

вого

числа k 2.

Таким

образом, в

вид /e2Oi

1.

Это значит,

этой области условие (77) принимает

что масштаб изменяемости

начальных

искривлений

в направлении

контура должен быть мал по сравнению

с характерным размером оболочки в направлении по нормали к кон­ туру. Этот вывод находится в соответствии с принципом Сен-Венана.

После того как функции q>(jci | к) и ф (л^ | к) найдены, дисперсии, корреляционные функции и т.п. характеристики вычисляются Извест­ ными методами (§ 1.12). В частности, для среднего квадрата перемеще­ ния аУх(г) имеем формулу

оо

< К (Г ) |2> = $ S ^ IO lT fo llO IM k .

(85)

— оо

 

Как и следовало ожидать, в рамках сделанных допущений

средний

квадрат не зависит от координаты х 2.

 

§ П.5. Растяжение пластины с начальными неправильностями

Пусть пластина с малыми начальными неправильностями испыты­ вает равномерное растяжение силами N (рис. 27). В этом частном слу­ чае в уравнениях (68) и (74) следует положить Ri~> оо; R 2-*• оо; ц =-

= N 2 = N. Второе уравнение системы (74) дает ф = 0, а первое vDaBнение принимает вид

D (<pIV—2k\ -J- ^2 ф) N —^2 <р) = — N -\-k\) е‘к,х‘.

Выпишем общее решение этого уравнения

 

Neikl х*

 

(86)

а = I

~i>№+ty+N

где га — корни характеристического уравнения

D { r * - k i r - N { r > - k l ) = 0.

Эти корни, очевидно, составляют

Г1,2 = ^ ^ 2 ' Л3.4 = ± V

k\ Jt-Xl,

где использовано обозначение

 

х2= "

(87)

D

Пока выполняется условие х2> —А:|,два характеристических кор­ ня — отрицательные действительные числа, т. е. концепция краево­ го эффекта применима в данной задаче. При положительных х2 (т е при растягивающих усилиях) эти корни по модулю имеют порядок k2 и более. Таким образом, можно ожи­

дать выполнения условия (77) и быст­ рого затухания краевого эффекта. При сжимающих усилиях краевой эффект затухает медленнее. Условие х2 = —k\,

как следует из (87), соответствует ра­ венству

N = N , ( 0, k2) = - k 2D.

Здесь N#(0, k2) — критическое усилие по линейной теории, соответствующее волновому числу k2 (при этом kx = 0). Итак, при достижении критического усилия краевой эффект вырождается.

Пусть х2 > — k\. Тогда решение (86), ограниченное на бесконеч­ ности, примет вид

<р = С1 е к,х'-\-С2е f 4 + * 2'i

х2 е‘к,х'

(88)

 

ft» + x 2

 

Как и ранее, используется обозначение А2 = k] Вычислим по­ стоянные Сх и С2для случая жесткого защемления пластины на линии

хх — 0. Требуя, чтобы выполнялись условия ф(0) — ф'(0) = 0- ле1*ко получим

х2 (К * !+ х 2+ /*,)

(*2 + х2) ( ] Л 2+ х 2- * 2) '

________ X2 (fe2+ tfei)_______

(*2 + х2) (]A jj + x2- * 2)

Выражение (88) после небольших преобразований записывается сле­ дующим образом:

Ф (*i | к) = ---------------------------------

[[ V k \ W + ikJ

*• -

{ k * w ) { V k i w - h )

 

(/:2+ ^ 1) Г ^ А2+х2

( V k \+ v ? —k ^ elklXl\ .

(89)

Применим формулы (85) и (89) для вычисления среднего Квадрата (дисперсии) перемещения при условии, что поле начальных неправиль­ ностей ш0(г) представляет собой узкополосное поле со спектральной плотностью, аппроксимируемой выражением:

CTQ

(90)

S * (к) = ^ - б ( | kx |- Л 0) б ( | k21- k 0).

Здесь ст„ и k0 — некоторые положительные постоянные. Подстановка выражения (90) в формулы (85) и (89) дает

< | w(г) |2> = о\

) 2 [(«о-1)2 + (2 + л) в-2***1+ 2е~2п°к°*>-

2 (1 + п0)е~ (|+я<,) *•

2ё~к°Хг (п01) (cos k0x1—sin£0Arx)-f

 

+ 2е~п° k°Xl (n0— 1) (cos k0 x1+

sin k0xx)— 2ne~k‘ Xt cos k0*x] (

(gi)

где использованы обозначения

 

 

 

n =

;

n0= ]/rl + n.

(92)

 

 

k 0 ™

 

 

Заметим, что правая часть формулы (91) остается конечной при я >

> —2 .

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

2

< 1^ ( 0°, х2) |2> =

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

кое значение.

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

ного преобразования. Рассмотрим вопрос о вычислении этих характе­ ристик в общей форме. Пусть т(г) — некоторое поле, связанное с по­ лями ^ (г) и %(г) соотношением

m = l1

Здесь

и / 2 — линейные операторы. Каноническое представление для

поля

т(г) имеет вид

оо

/я(г) =

$ №0(k)P(*ilk)e'*,Jr,A

 

(93)

 

— оо

 

 

 

где функция |х(хх | к) определяется следующим образом:

 

ц (х,|к)= [ф (*!|к)е1к

Хш] + /2 (ХХI к)е‘кгх‘] I е~1к

х\

Корреляционная функция поля ш(т) находится с учетом

соотно­

шений (70) и (93):

 

 

 

 

<m *(r)m(r')>= $

■Sa,0(к)

(дсх | k) ^ (jci |к)

_дг*)rfk.

— оо

 

 

 

 

В частности, полагая г = г \ получим формулу для среднего квадрата

оо

( Im (г) |2 > = $ •V ,(k)|p(*1 |k )|2dk.

(94)

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]