1200
.pdfРезультаты вычислений по формулам (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ё~к°Хг (п0— 1) (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) |