книги / Устойчивость и колебания пластинок и оболочек с отверстиями
..pdfОбщее решение уравнения (7.14) представляется в виде интеграль ных уравнений:
•йУ: |
г т { - [ / 2 (*) z" ~ Л (л) z lv] - |
|
/о (■*) |
— 2 |
I / 2 W ГУ4" /4 W гу] Л/з £Tix — |
— |
U 2 (*) г) + |
/ 4 |
(*) г)\ ( - |
l )y+1 |
t rJ* ^ (ле |
; (7.22) |
|
У-1 |
|
|
|
|
|
|
4 |
|
|
|
|
|
wl= z l + |
^ Ayoz-yerrv+ |
^ |
( - 1У+1 |
o e r^ ^ jxe г^ х ; |
|
|
|
ySt |
7 S1 |
|
о |
(7.23) |
|
|
|
|
|
|
|
|
|
|
|
( - 1 ) * 1-D ГУ ]Х\ **ГЧХС1Х |
|||
Здесь A jo— 'произвольные постоянные интегрирования |
( / = 1, 2, |
3, 4); функция z определяется линейным дифференциальным урав нением с постоянными коэффициентами:
aAzlv -\-(izZu -{~aQZ = 0. |
(7.24) |
||||
Параметр D есть определитель Вандермонда: |
|
||||
1 |
1 |
1 |
1 |
|
|
г х |
г2 |
|
г 3 |
г4 |
|
D = _2 |
2 |
|
2 |
2 = const ^-0, |
(7. 25) |
Гх |
Г 2 |
|
Г 3 |
Г 4 |
|
-3 |
_3 |
|
3 |
з |
|
Гх |
Г 2 |
|
Гз |
Г 4 |
|
где dj — определитель третьего порядка, получающийся из матри цы определителя D путем вычеркивания первой строки и /-го столбца.^
Функция [х определяется зависимостью вида
= |
б0 -j- |
—[—02^50^^ —f—бз'гед*11 —|— |
(7 .2 6 ) |
||
гд е |
в0= 0 ; |
в1= |
/ „ ( Л ) - * 1 |
_ 1 ; |
|
|
|
|
а0 |
|
|
9 _ _ / » (* )/о М ~ / о |
{X) /2(х) |
0L _ a L . |
|||
|
|
/о (•*) |
|
Л0 |
д0 |
|
|
|
|
|
(7 .2 7 ) |
a |
fi( x ) a 4— f 4{x)a2 . п |
fo(*) |
1 |
до /4(*) |
/о (*) во |
Исследуемую полоску берем по ширине, рав |
р |
||||
|
|||||
ной расстоянию между центрами вырезов вдоль |
|
||||
дуги. В этом случае уравнение |
изгиба полоски |
|
|||
запишется в форме (7.14). Его -решение ищем по |
|
||||
методу последовательных приближений Н. К. Ку |
-й -тз-й - |
||||
ликова. Оно записывается в форме (7.22). |
|
h It |
|||
Чтобы решить нашу задачу в нулевом |
приб |
|
|||
лижении, |
следует в (7.22) предположить, что |
ч Т г т Г П |
|||
р = 0 . Так как общее решение должно удовлетво |
|||||
|
|||||
рять граничным условиям задачи, то, подставляя |
Рнс. 7.1 |
||||
в них (7.22) |
и (7.29), .получим: |
|
|
||
|
|
|
|||
|
r i Л ю - { - Г 2Л 20 |
г з Л 3 о / ' 4 ^ 4 о = 0 ; | |
|
||
|
гie * Лю-Ьггв * Лго-Ь^зС^ Л3о-|-,‘4вг<*Л4о= 0;]3 |
||||
Л [/2 (0) ~ r \]iA\o-j- Л [/2 (0) — г1| Л2о+ |
гз [/2 (0) — rl]Лзо+З |
||||
|
+ г4 [/2 (0) — Л] Л40= 0 ; |
(7.28) |
Л[f2 [I)— r'\\ Л 10 4~ Л[/2 (/) — rl]ёг‘*Л2о + r3j[/2 (7)—r| j er,f Л 30-f-
—^4] er*^ 4 o = 0 .
Эта система имеет решение, отличное от нуля лишь в том слу чае, если определитель, составленный из коэффициентов при неиз вестных Aj[0| равен нулю. Раскрывая его и учитывая значения кор ней гь г2, г3 и г4, определяемые соотношениями (7.21), получим два -алгебраических уравнения:
г2 — г ? = 0 ; |
|
(7 .29) |
| Q—(ri4Ta)/ |
__Q |
(7 .30) |
|
|
Рассмотрим второе из них. Оно легко преобразуется, в результате имеем
|
cos {rt-f- г2) / — COS (rj — r 2) l = |
0 |
(7 .3 1 ) |
||||
или |
sin r x/ = |
0; |
sin r2l = 0 . |
|
(7 .3 2 ) |
||
Отсюда |
r it = n m |
m = 1, |
2, |
3 ,... |
|
(7 .3 3 ) |
|
|
г21 = л т |
m = 1, |
2, |
3 ,... |
|
(7 .3 4 ) |
|
Освобождаясь от корней, из зависимости |
(7.33) |
найдем |
|
||||
|
^ 2 = —r |
flo"fa2» |
|
(7 .3 5 ) |
|||
|
|
a* |
|
|
|
|
|
*где |
|
а==щ л//. |
|
|
|
|
Если число полуволн т достаточно велико, то можно записать условие минимума коэффициента a2t включающего внеш нюю нагрузку Р, приравнивая нулю производную от а2 по а:
— |
- |
а 0+ 2 а = 0 , |
|
da |
аз |
0 1 |
|
откуда |
|
4 — |
(7 .36) |
|
а = у а0. |
Подставляя (7.36) в (7.35), найдем выражение для критической силы
р _ 2 Уао
(7 .37)
«2
Такое же выражение получается из уравнения (7.34) и (7.2Э). Таким образом, мы получили приближенное выражение для критической нагрузки, действующей на полоску, выделенную из оболочки, ослабленной вырезами произвольной формы. Переходя от него к критическому напряжению для перфорированной оболоч
ки в целом, найдем
Эта формула получается из (7.37) путем деления на Л.
Решения (7.37) и (7.38) представляют собой лишь нулевое при ближение окончательного результата. В -случае необходимости за дача может -быть решена с более высокой точностью. Как показы вают решения других задач устойчивости [44], выполненные с по мощью данного метода, погрешность нулевого приближения, в сравнении с точным, получается небольшая, в пределах 10% . Для нашей задачи погрешность будет зависеть от характера изменения жесткости полоски от участка к участку, а также геометрических размеров оболочки и вырезов.
Рассмотрим, каким образом осуществляется более точное ис следование.
По окончании решения задачи в нулевом приближении ищется
выражение для функции р путем подстановки выражений до' и w"y найденных ранее, в отношение (7.26). При этом учитываются за висимости (7.21). Кроме того, определяются коэффициенты функ
ции р по (7.27) с использованием (7.19). Имея выражение для р,
можем записать общее решение уравнения (7.14) в первом при ближении. Удовлетворяя граничным условиям, вновь приходим к системе однородных уравнений относительно Аю, А 20, Лзо и Л40.
После приравнивания определителя системы, составленного из ко эффициентов при неизвестных А#, нулю, получим алгебраическое уравнение, из решения которого совместно с (7.21) находятся не известные величины г: и Р. По окончании расчета первого прибли жения вместо (7.38) получим новое выражение для критического напряжения:
ак р |
( 7 . 3 9 ) |
где | — коэффициент ( £ < 1), изменяющийся от приближения к при
ближению и стремящийся к своему постоянному значению. Расчеты второго и всех последующих приближений выполняются аналогич но первому.
Если предположить, что диаметр кругового выреза значительно меньше длины и диаметра оболочки, то справедливо будет счи тать действие реактивной силы отброшенной части оболочки оди наковым по всей длине выделенной полоски.
Далее рассмотрим некоторые частные случаи. Вначале иссле дуем устойчивость нагруженных осевой силой оболочек с опреде ленным числом вырезов квадратной формы. Прежде всего, подстав ляя в (7.38) зависимости (7.15) и (7.5), получим
( 7 . 4 0 )
где сткр определяется формулой (7.11) или (7.13),
да у |
( 7 . 4 1 ) |
i J и * ) |
|
f |
О |
Судя по зависимости (7.41), коэффициент к зависит от числа, фор мы и места расположения вырезов.
Рассмотрим устойчивость цилиндрической оболочки, ослаблен ной четырьмя вырезами, расположенными симметрично относитель но оси на средней параллели. Длина стороны выреза равна а. Тогда расстояние между центрами вырезов будет равно /i= 2jtR /4. Момент
инерции полоски шириною 1\, мысленно выделенной из оболочки, |
|
будет изменяться ступенчато вдоль образующей на трех участках. |
|
В пределах каждого участка он постоянен. На участках от нуля до |
|
(I — а) 12, а также от |
(1 + а )/2 до I жесткость полоски единичной |
ширины определяется |
по зависимости (7.16). На участке с коорди |
натами ( / — а )/2 . . . ( |
1 + а ) /2 жесткость полагается |
постоянной по |
всей ширине полоски, |
равной среднему значению, |
получающемуся |
к
О |
0,7 О,U |
0,в |
0,6 ф |
|
Рис. |
7.2 |
Рис. 7.3 |
делением истинной жесткости на ширину участка 1\. Тогда для по лоски единичной ширины она равна
(7.42).
Подставляя (7.16) и (7.42) в (7.41) и произведя интегрирование, получим
На рис. 7.2 приведены кривые, показывающие характер изме нения величины параметра k в функции от отношения длины квад ратного выреза к диаметру и к длине оболочки. О степени влия ния на устойчивость числа вырезов можно судить по графику, при веденному на рис. 7.3. Расчет был выполнен для оболочки с пара метрами a l(2 R )= 0,2 и а /1= 0,4 и различным числом квадратных вырезов, расположенных в среднем сечении.
Предлагаемый метод может быть применен для исследования устойчивости оболочек с подкрепленными вырезами. Подкрепле ние подбирается в соответствии с решением задачи прочности тако вым, чтобы оболочка получилась равножесткой во всех направле ниях. Если же подкрепление задано, то площадь его поперечного сечения можно равномерно распределить («размазать») до толщи ны стенки оболочки, уменьшая действительное отверстие. Форма фиктивного эквивалентного выреза аналогична действительной, а его контур эквидистантен.
По поведению вырезанных из оболочек полосок проводилась оценка устойчивости оболочек с вырезами в целом, результаты ко торой даны в [119]. Разница состоит лишь в том, в сравнении с дан ной работой, что авторы ее подбирали коэффициент, учитывающий наличие вырезов, из геометрических соображений, а не выводили с учетом изменяющейся жесткости полоски и переменности коэф фициента постели. В этой связи расчеты по зависимостям работы [119] дают несколько завышенные результаты в сравнении с рас сматриваемой работой.
7.3.ЦИЛИНДРИЧЕСКИЕ ОБОЛОЧКИ ПРИ ВНЕШНЕМ ДАВЛЕНИИ
Рассмотрим устойчивость длинной цилиндрической оболочки, 'ослабленной системой прямоугольных вырезов, расположенных на равном расстоянии друг от друга вдоль образующей. Предполага ется, что при потере устойчивости кривизна оболочки изменяется лишь в окружном направлении, а образующие остаются прямоли нейными (рис. 7 .4). Это предположение позволяет приближенно судить об устойчивости оболочки по устойчивости кольца, мыслен но выделенного двумя перпендикулярными оси сечениями. Анало гичные предпосылки были приняты в работе [68]. В отличие от нее здесь избран другой метол решения, позволивший получить окон чательные простые зависимости для определения величины крити ческого внешнего давления.
Кольцевой элемент, мысленно выделяемый из оболочки, при нимался по ширине равным расстоянию между центрами вырезов по образующей. Действующий на него изгибающий момент М, как известно [93], связан с изменением кривизны следующим соотно шением:
£ /у.= М , |
(7.44) |
тде Е — модуль упругости; I — момент |
инерции поперечного сече |
ния кольца, который для исследуемого случая является функцией угла 9.
Решение поставленной задачи ищем с учетом малых дефор маций.
Путем простых операций дифференциальное уравнение (7.44) преобразуется к виду
d2w |
W = MR2 |
(7.45) |
Пт |
El |
|
где Я — средний начальный радиус кривизны кольца.
Изгибающий момент появляется вследствие радиального упру гого смещения. Он равен M = N w , где N — нормальное усилие, возникающее в поперечном сечении кольца. Величину нормально
го усилия принимаем |
так |
же, |
как в работе |
«а |
]68], постоянной и равной 1среднвму значению: |
|
|||
N = p R l (* |
2^яГ )' |
<7 - 46> |
|
|
Здесь п — число прямоугольных вырезов в по |
|
|||
перечном сечении оболочки; |
5 — площадь се |
|
||
чения одного выреза. |
|
|
|
Рис. 7.4 |
Подставляя (7.46) в выражение для момента |
M — Nw, а по |
|||
следний в уравнение (7.45), получим |
|
|||
|
|
d2w |
k w = 0, |
(7 .47) |
|
|
d02 |
||
|
|
|
|
|
где |
k = k (B ) |
p& L Г i _ |
(7 .48) |
|
|
|
£ /(0 )l |
|
Функция прогиба, найденная из решения уравнения (7 .47), дол жна удовлетворять условиям симметрии изогнутой формы, выра жающимся в том, что
———= 0 при 0— 0 |
и 0 = — . |
(7.49) |
dd |
2 |
|
Линейное дифференциальное уравнение (7.48) с переменным коэффициентом k может быть решено, как и в предыдущем разде ле, методом последовательных приближений. Сравнивая уравнение (7.47) с (7.16), найдем, что
п = 2; / 2 (0 )= 1 ; / 2 (0)=О ; / 1(6 )= 0 ;
' • №- ' ^ |
[ 1- w |
] |
+ |
, - ' F ( e ,+ li |
(7- 50> |
||
fi( 0 )= p F '( e ) ; |
/ |
(0)— 0. |
|
||||
Характеристическое уравнение (7.18) |
будет иметь вид |
|
|||||
r 2-f-a1r - }- a 0= 0 . |
|
|
(7 .5 1 ) |
||||
Решая его относительно г, получим, что |
|
|
|||||
гх— — a J 2 -\-k] г2= |
— a J 2 — k, |
|
если |
k2= a i /4 — % > (); |
(7 .52) |
||
Г\= ~aJ2-\-qi\ г 2 = |
—aJ2—qi, если |
д2—а0—а ? / 4 > 0 . |
|||||
|
Коэффициенты уравнения (7.51) а0 и а\ пока неопределенные,
т\ФН Оот^О. Уравнение -(7.23) для нашей задачи имеет вид Кроме того, согласно (7.24) имеем
D = 1 1 |
— г2 — f i = co n st ф 0; ах= г 2; а2— гл. |
г2 |
|
Подставим далее (7.50) в зависимости (7.25) и (7 .27):
|
|L=A1-a// -j-X2'ZE//; |
|
xl = p m . - i ; x2 = - f i w _ i i |
(7 .5 3 ) |
|
«о |
Оо^(в) Ч |
|
Общее решение дифференциального уравнения (7.47) в соответст вии с (7.23) имеет вил
w |
0 Чъ) |
+ Г^ 2°е''’ + |
^ |
J ^ " >м ~ |
||
= ~ |
||||||
|
|
e r»8 |
^ (ie |
r*eof0l; |
|
|
|
|
ГЧ. |
|
|
J |
|
|
|
Оо |
|
|
||
|
|
а I |
|
*У2 |
£ er«e ^ jxe_ r‘flrf0 — |
|
|
■ш' = ЛюГ1ег,в+ ^ 2 о ^ е г»0 -| |
|
||||
|
|
|
|
П — гх\ |
|
|
|
|
— ег»в|‘ (1.е |
|
Га°^/вJ; |
(7.54) |
.^"=Л 10г1ег,°+ Л20г|ег , в |
Гrxe^e ^ [Te_ri8d6— |
1 |
о0 |
— г2ег»еf {Ге~Г1°^б1. |
|
Оо |
J |
Вычисляя ц с различной точностью и подставляя его в (7.51), по лучим общее решение в соответствующем приближении.
Рассмотрим вначале нулевое приближение. В этом случае jx=Of тогда решение преобразуется таким образом:
® = |
j (г1-4 1оег‘*+''2Ла)е'>1); |
|
-а!>'=Л10г1е'««+Л 20г2е'*«; |
(7.55) |
=Л10г?ег'° + Л2ог1ег»0.
Вполученное решение входят неопределенные параметры Сь До
инеопределенные корни Г\ и г2. Для нахождения а<> и а\ используем метол, применяемый при решении интегральных уравнений, суть которого состоит в осреднении соответствующих функций. На ос новании этого неопределенные параметры могут быть вычислены по следующей формуле:
/,(« )< » . 7 = 0 > ] > 2 . - м П. |
(7 .56) |
Осредняя функции /0(6) и f i(0) на рассматриваемом интервале по зависимостям (7.55), найдем
a' = 1 |
7 |
= |
^ /,(e )d e ; |
|
|
О» |
|
|
Оо |
0, |
|
|
|
|
|
||
а 0= — I— \ f ( Q ) d B = — £— ^ ( 6 ) ^ 0 . |
(7 .5 7 ) |
||||
ai— ао J |
|
|
®i— во |
«) |
|
во |
|
|
|
О0 |
|
Затем по зависимостям (7.57) |
и |
(7.52) |
определим |
неизвестны^ |
|
корни |
|
|
|
|
|
r x=qi\ r2 = — qi\ Ч= |
л |
/ |
^ ( е)^ е + |
(7 .5 8 ) |
|
|
|
Y |
|
Оо |
|
Решение (7.55) содержит |
произвольные постоянные А \0 и Л20. |
Их можно найти, удовлетворяя дополнительным условиям, которые следуют из условия симметрии выпученного кольца. Эти условия достаточно выписать для двух .поперечных сечений А к Б (см. рис.
7.4) с координатами 0 i= O и 02= |
jt/2: |
|
700,-0= 0; |
7О0,-«/2=О. |
(7 .5 9 ) |
Объединяя (7.5Э) и (7.55), получим систему алгебраических уравнений, с помощью которой можно определить постоянные Лю и Л 20*'
r 1 еГ1®' Л 10+ г 2ег»°» Л20= 0 ; гхe rifl* Л 10-f г2 ег*в» Л20= 0 . (7 .6 0 )
Эта система будет иметь решение, отличное от нудя, в том слу чае, если определитель, составленный из коэффициентов при неиз
вестных, равен нулю: |
|
|
гхег,0> |
г2 ег*в‘ = 0. |
(7 .6 1 ) |
гхег,в> |
г2ег*в» |
|
Учитывая соотношения (7.58) и раскрывая определитель, най |
||
дем, что sm«/i(02 — 0 i ) = 0 . Отсюда легко определим |
наименьший |
положительный корень, отличный от нуля: <7 (02 — 0 i)= jt или, под
ставляя вместо 02 и 0i их значения, найдем q = 2 . Найденное зна чение q подставим в (7.58), в результате получим формулу для вы числения критического давления в нулевом приближении:
3 (Bi — 0О)
(7 .6 2 )
0>
|>(б)<го
Рассмотрим решение конкретных задач. Прежде всего, если предположить, что мы имеем оболочку без отверстий, тогда из фор мулы (7.62) легко получается известное соотношение:
P i с р — 5 |
EI |
(7 .6 3 ) |
|
№ |
|||
|
|