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

книги / Устойчивость и колебания пластинок и оболочек с отверстиями

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

точности совпадают с результатами, полученным* с использованием «сплошной* модели с неоднородной жесткостью и метода Бубнова — Галеркина, то, очевид­ но, правомерно распространить метод сплошной модели в сочетаний с методом Бубнова — Галеркина на другие задачи.

К настоящему времени сведений о критических нагрузках для пластинок с отверстиями явно недостаточно для практики. В известной мере это объясняет­ ся тем, что применять обычные формы записи энергетического критерия устой­ чивости можно только, выполнив соответствующие предварительные условия [45, 46]. Та<к, для использования критерия в форме Брайена [46] необходимо ре­ шать плоскую задачу теории упругости, что само по себе не является простым для неодносвязных областей. Критерий в форме Н. А. Алфутова и Л. И. Балабука [1, 2] требует построения статически допустимого напряженного состоянчя. Использование критерия в форме С. П. Тимошенко [Э4] связано с необходимо­ стью предварительного определения бифуркационных перемещений точек контура пластинки. В связи с этим Л. М. Куршшым и К. А. Матвеевым [45—47] для решения задач устойчивости пластинки с отверстием предложен критерий в фор­ ме, не требующей выполнения каких-либо предварительных условий. Не при­ водя здесь вывода для критериального соотношения, отметим, что в основе его лежит условие: задачам вариационного типа соответствует не один, а целое мно­ жество функционалов. Переход от одного функционала к другому осуществля­ ется с помощью неопределенных множителей Лагранжа и некоторых преобразо­ ваний. При этом предварительные условия, которым должны удовлетворять иско­ мые функции, переходят в естественные условия для некоторого преобразованно­ го функционала. Вариационная задача для нового функционала дает, в частно­ сти, естественные условия для определения вводимых при преобразованиях функ­ ций— множителей Лагранжа. Сделанное ранее сравнение результатов вычисле­ ний по зависимостям Л. М. Куртина и К. А. Матвеева с решениями, приведен­ ным* в других работах и опытными данными, изложенными в предыдущем раз­ деле, показывает, что решение задач устойчивости для пластинок с круговым отверстием на оонове критерия работы [45] приводит к значительно занижен­ ным критическим параметрам, особенно с увеличением размера выреза.

Глава 4

УСТОЙЧИВОСТЬ ПРЯМОУГОЛЬНЫХ ПЛАСТИНОК С ОТВЕРСТИЯМИ ПРОИЗВОЛЬНОЙ ФОРМЫ

4.1.АНАЛИТИЧЕСКИЕ МЕТОДЫ ИССЛЕДОВАНИЯ

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

(ЛГ,р>+ЛГха ? )-5 -= О о г - £ = £ l —

_D°iiujbAr( e+

~

v н ’ (4-ц

2

т '‘Ы

 

/=1 '/=1

 

 

 

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

£ - = (а 2- г №

b*ij=y2ij -

уир

Q = m ujb*tjV + m2lJa*jа;

 

mXij =

sin 2 a.xuj — sin 2ax2ij\

(4 .2)

m.2,7 =

sin 2 $yw — sin 2 py2iJ;

 

S = a b ;

/- i ] - 1

Здесь через Si обозначена площадь, занимаемая вырезами.

Конфигурация вырезов сказывается на методе определения ко­ ординат вершин прямоугольников Хщ; х2ц\ унз\ Ущ, с помощью

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

Если пластинка перфорирована только вырезами в форме полу­ круга (рис. 4.1), то координаты вершины прямоугольников для этого случая находим по формулам

у щ = У 1— Ri 4 - АУ(Л

У г и = У ш + Д&;

 

Д

R 2i — ^ и~ У- — у ^

(4 .3 )

X2lJ= Xt\ У— 0,

1, 2 , . . . , г.

 

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

Далее рассмотрим пластинки с вы­ резами эллиптической формы. Пред­ полагаем, что они ориентированы та­ ким образом, когда оси эллиптическо­ го контура параллельны внешним сто­ ронам пластинки.

Параметрами эллипса являются координаты центра Х{ и yit а также малая — Ь* и большая а* полуоси. Кроме того, форму эллипса характе­ ризует ее эксцентриситет

у\ Ny

шшшшгшшхшш

NA !

м w

Vi’Ус Ny

Рве. 4.1

e = V { a * ) 2- {b * fla * .

(4 .4)

Координаты прямоугольников, представляющих эллиптический' вырез, будут определяться по формулам

У\)=Ш — £ *+ Д ^ £

*; Д ^ = 2

6 */г ;

 

 

k * = 0, 1, 2 , . . . , (

г - 1 ) ; / =

3 , 4

, . . . , г,

(4 .5)

V (^*)2 — 0>25 (y2j-\-yij — 2yjf.

Наконец, если контуры прямоугольного выреза не параллель­ ны наружному контуру пластинки, то в этом случае необходимые для вычисления соотношения значительно усложняются. Ниже бу­ дем рассматривать этот более общий случай. Для упрощения за­ писи индекс i опускаем. Вершинам прямоугольников присвоим по­ рядковые номера 1, 2, 3 и 4. Они имеют координаты соответственно (*ь У\)> (хь Уг), (*з. Уз), (*4, Уа) •Предположим, что последние под­ чиняются условиям * 1 < * 2 < * з < * 4 ; Уч<У\<У\<Уз-

Прямоугольный вырез для рассматриваемого случая представ­ ляется серией малых прямоугольников с контурами, параллельны­ ми сторонам пластинки. Высота этих прямоугольников равна

Ьу=1\1п Г=3, 4,...,

где

Y

[(*2 — Х& -

(У2~

[(ха — xzf + (у2у&]

1 г = ----------------------------------

 

 

■*з---------------------------------------

- ;

 

 

 

 

Х1

 

х 'з= х 1 +

{х з -

Xi) {УА—Ух)1{уг— t/x).

 

В этом случае ординаты прямоугольников будут:

 

У ц = У 2+

Ьу&;

k * = 0, 1,

2 , . . . , (г — 1);

Уч]=У ц+ьу;

</ср.у=0,5 [yl} -f у2}).

Если У2< У ц < У и то абсциссы

этих

прямоугольников находятся;

по следующим формулам:

 

 

х И = Х г- f 2— X j

(ycpJ y.)l(y2 — уf);

XV = X 2 + C*4- X2) (yc[)j -

y2)/(y2 - 1/4).

Когда yx< y 2 < .y A, TO l2=y-& y2 2lx, L y '= l 2jr.

В этом случае

y2j = y i j + b y r;

y'cPj= 0 ,5 (уу+у'у);

—*i) (У ы — уд/УЛ —Уй Х у = Х2+ {хА— лгх) (г/ср j —уЖ & . — 0.1).

Наконец, когда Уа< У 2) < У з, то

уЬ = t /4 + Ду#*; г/2у= г/1/ + Дг/;

*/ср/=0,5 (^1у+г/2;);

* y = * i + ( * а —* 1 ) (t/cp/— £/i)/(i/3— i/,):

^ 27= ^ 4 - f ( * 8 x 4) ( y cpJ - г/4)/(£/з - У а ) ‘

Для более ясного представления о степени влияния формы вы­ реза на устойчивость прямоугольных пластинок были выполнены вычисления по зависимостям (4.1) и (4.2) с учетом значений ко­ эффициентов для вырезов соответствующей формы.

На рис. 4.2 показан характер изменения коэффициента при ком бинированном нагружении для прямоугольных пластинок с цент­

ральным круговым вырезом

(сплошные кривые) и для

пластино.с

с вырезом в форме полукруга

(пунктирные).

 

 

На рис. 4.3 показан

характер изменения коэффициента

дня

квадратных пластинок с

эллиптическим центральным

вырезом в

функции от эксцентриситета.

 

 

 

В табл. 4.1 приведены значения коэффициента kx для пластин­

ки (см. рис. 4.1)

с различным числом вырезов в форме полукруга.

Координаты Xi/a

и y ja ,

где / = I, 2, 3......... 16, определяются

следу-

к

3ft

2,0

1,0

О 0,2 0,4 0,6 0,8 е

Рис. 4,3

ющими правилами: X i/ a = 0,2j ( / = 1

.для 1 = 1 , 5, 9, 13; } = 2 для 1 = 2 ,

6,

10,

14; j = 3

для i = 3, 7, 11, 15;

j = 4

для i= 4 , 8, 12,

16); г/,/а=

=

0,2г]

(т]=1

для 1=1, 2, 3, 4; т)—2 для

1 = 5 ,

6, 7, 8; rj =

3 для

1 =

= 9 , 10, 11, 12 и т] = 4 для 1 = 1 3 ,

14, 15, 16).

Все вычисления

про­

водились с точностью, равной 1C-7. Время расчета одного вариан­ та лежит в пределах 1 . . . 2 мин.

Т а б л и ц а 4.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значения \ = N y/ N x

 

 

 

 

 

 

 

 

Номера отверстий

0

0,2

0.4

0,6

0,8

1,0

 

 

 

 

 

 

 

 

 

 

 

 

 

1, 6,

 

11,

16

 

 

 

 

 

 

 

3,74

3,12

2,68

2,34

2,08

1,87

1,

16

 

 

 

 

 

 

 

 

 

 

 

3,8о

3,21

2,75

2,40

2,14

1,92

1

 

 

 

 

 

 

 

 

 

 

 

 

3,90

3,25

2,78

2,44

2,16

1,95

1, 6,

 

16

 

 

 

 

 

 

 

 

 

3,80

3,16

2,71

2,37

2,11

1,90

9,

10, 13, 14

 

 

 

 

 

 

3,74

3,12

2,68

2,34

2,08

1,87

9,

10,

13,

14,

15

 

 

 

3,70

3,08

2,64

2,31

2,05

1,85

5,

9,

 

10,

13,

14,

15

 

3,64

3,04

2,60

2,28

2,03

1,82

1, 5, 6, 9, 10, И ,

 

13— 16

3,45

2,87

2,46

2,15

1,91

1,72

1, 4, 5, 6, 9— 11,

 

13— 16

3,40

2,83

2,42

2,12

1,89

1,70

1— 16

 

 

 

 

 

 

 

 

 

3,32

2,76

2,37

2,07

1,84

1,66

1, 6,

7,

10,

11,

16

 

 

3,65

3,04

2,60

2,28

2,03

1,82

1, 4, 6, 7,110,

 

11,

13,

16

3,55

2,96

2,53

2,22

1,97

1,77

4.2.ЭКСПЕРИМЕНТАЛЬНЫЕ ДАННЫЕ

ИСОПОСТАВЛЕНИЕ ИХ С РАСЧЕТНЫМИ

Экспериментальному изучению поведения прямоугольных пла­ стинок с вырезами овальной формы посвящено, насколько извест­ но автору, лишь одно исследование В. Г. Налоева [58]. Целью этой работы было выявление влияния размеров выреза на величину ко­ эффициента устойчивости k0. На специально сконструированной установке [57] были испытаны в условиях одноосного сжатия рав­ номерно распределенной нагрузкой четыре прямоугольных пластин­ ки с отношением сторон 1 : 2 ; 2 : 1. Опытные пластинки изготовля­ лись из органического стекла. Выбор последнего основывался на следующих соображениях: во-первых, органическое стекло легко поддается обработке, что немаловажно, так как в каждом образце при исследовании производилось последовательно по три выреза; во-вторых, органическое стекло — линейно упругий материал с низ­ ким модулем Юнга, что позволило ограничиться при испытаниях невысокими нагрузками и в то же время испытать достаточно толстые образцы 5 мм) с малым начальным прогибом.

Для определения модуля Юнга материала экспериментальных образцов из оргстекла той же партии была изготовлена отдельная

Mktmt

1

w

0,5

6)

00,5

Pmc. 4 .4

 

 

1

i.

 

 

 

i

 

сl________

 

1

N

i

J

 

1

 

 

 

 

I

 

 

 

i

i

1________ J_________

1,0 1,5 c/(f

Ряс. 4.5

балка. Ее испытания показали, что £ = 3 , 0 6 ГЛа. Коэффициент Пу­ ассона считался равным 0,35.

Соотношения линейных размеров экспериментальных образцов, изготовленных согласно рис. 4.4, а, были следующие:

Пластинки I и II: а/Ь = 0,5; с/о=;0,5; dfc= 1;. 1,'5; 2 Пластинки III и IV: а/Ь=2; d/b= 0,5; c/d= 1; 1,5; 2

Через а здесь обозначены стороны, перпендикулярно которым при­ кладывалась внешняя нагрузка, через Ъ обозначено расстояние между опорными ножами. Длина сторон а опытных образцов у пластинок I— II была равна 448 мм, а у пластинок III и IV — 896 мм.

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

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

Каждый из четырех опытных образцов, отобранных для основ­ ных исследований, испытывался четыре раза: сначала как оплош­ ная пластинка, а затем три раза с вырезами, для высверливания которых образец вынимался всякий раз из установки.

Результаты испытаний образцов без выреза вполне удовлетво­ рительно согласовались с вычисленными значениями критической силы по известной классической (методике Саусвелла — Доннелла для оплошной пластинки. Форма потери устойчивости испытанных образцов также соответствовала расчетной. Наличие выреза не из­ меняло формы потери устойчивости широких (I, II) пластинок, но, как видно из рис. 4.5 (точки залитые), существенно снижало зна-

чение коэффициента устойчивости. Что касается длинных пластинок типа III, IV (точки не залитые), то вопрос о влиянии вьгреза на их устойчивость не удалось изучить достаточно детально, он тре­ бует дальнейшего изучения. Это вызвано тем, что появление выре­ за вносило существенные изменения в форму выпученной пластин­ ки при малых прогибах.

На рис. 4.5 нанесены не только экспериментальные точки, но и построены теоретические сплошные кривые, для чего использова­ лись необходимые соотношения предыдущего раздела. Кривая 1 — это вычисления для пластинок типа I, II, а кривая 2 — типа III, IV.

Помимо этого здесь же пунктиром нанесено решение (кривая <?), полученное методом конечных элементов для образцов типа I, II. Как видно из сопоставления, оба метода теоретического исследо­ вания и опытные данные дают критические параметры, близкие друг к другу. Пластинка и ее схематизация с помощью конечных элементов показаны на рис. 4.4 и б). Исследование устойчиво­ сти пластинок с центральным овальным вырезом на основе метода конечных элементов было проведено В. Г. Налоевым [59]. Благо­ даря двойной симметрии автор рассматривал только четверть пла­ стинки. В отличие от широко распространенной процедуры мето­ да конечных элементов В. Г Налоев рассматривал треугольный элемент как часть прямоугольного с теми же обобщенными коор­ динатами и тем же числом степеней свободы. Задача решалась в i два этапа. Сначала определялись компоненты напряжений, дейст­ вующих в пластинке в докритическом состоянии, в зависимости, ог нагрузки на наружном контуре, а затем найденные значения на­ пряжений использовались для решения задачи устойчивости, из ко­ торой находилась критическая сжимающая нагрузка. При решении плоской задачи теории упругости закон изменения компонентов перемещения произвольной точки конечного элемента применялся в форме, обеспечивающей линейный закон распределения нормаль­ ных и касательных напряжений в пределах элемента. Для упроще­ ния решения задачи устойчивости напряженное состояние в преде­ лах каждого конечного элемента предполагалось однородным.

При решении задачи устойчивости все элементы были разде­ лены на две группы: внутренние и внешние, примыкающие к на­ ружной кромке пластинки.

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

Аппроксимирующие функции прогиба для внутренних элемен­ тов задавались полиномом Мелоша:

w = a t-\-a^xJr a^y-f- алх 2 + а5ху - f а^у2 - f Д7* 3+

а^ 2У+

-f- а<рсу2 -]- а 10у3 - f а п х * у - \ - a l2x y * t

(4 .6 )

где х и у — безразмерные прямоугольные координаты. Для внеш­ них элементов, примыкающих к наружному контуру пластинки, ап-

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

(1 -

3 3 ? + 23?);

 

( 1 - 1 , 5 ? + 0,53?);

 

_

_

_

_

_

_

(4 .7 )

( х — 1,5л2+ 0 ,5 .? )

( х - 2 х 2 -\-х3).

 

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

Для сетки, показанной на рис. 4 .4 ,6 , решение задачи осущест­ влялось на ЭВМ «Минск-22 » с помощью стандартных программ,

входящих в математическое обеспечение машины.

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

можностью появления пластических деформаций

в

пластинке в

районе выреза. Этот вопрос может быть рассмотрен

по аналогии

с работой [59].

 

 

 

Условие, при котором определяется напряженное состояние пла­

стинки, может быть записано в виде

 

 

 

и)3э < <3пш

 

(4 * 8)

где

сгпп — предел пропорциональности материала

пластинки;

аэ

— критическое напряжение для пластинки по

Эйлеру; <о — ко­

эффициент концентрации напряжений в пластинке на границе вы­ реза.

Полагая модуль Юнга £ = 1 ,9 6 МПа и коэффициент Пуассона ц=С',3, можно получить из (4.8) предельное значение отношения толщины пластинки h к длине стороны а, при котором напряжения в пластинке могут достигнуть предела пропорциональности

fi/a < 0 ,7 4 3 - Ю-з 1 /о ппАсо/<э).

Значения коэффициента ю для некоторых типов пластинок, рас­ смотренных в этом разделе, получены в работе [58] для Ь{а = 0,5 ;

d / c = l ,0 и

d j b = 0,2;

С-,3; 0,4;

0,5; 0,6

соответственно равны

1,90;

2,36;

2,87;

3,47; 4,19,

а для

Ь /а = 0,5;

d j b = 0,5 и dfc = 0,5; 0,6;

0,7;

0,8;

0,9; 1,0 . Коэффициент

со

имеет значения 8,93; 6,49; 5,21;

4,41;

3,87;

3,47,

 

 

 

 

 

 

В заключение необходимо, по-видимому, обратить внимание на тот факт, что решение, приведенное в разд. 4.1, было получено без учета неоднородности докритического напряженного состояния. Однако, как мы убедились из Сопоставления с опытными данными и решением на основе метода конечных элементов, результаты вы­ числения но зависимостям, указанным в разд. 4 . 1, получаются

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

Глава 5

УСТОЙЧИВОСТЬ МНОГОСЛОЙНЫХ АНИЗОТРОПНЫХ ПРЯМОУГОЛЬНЫХ ПАНЕЛЕЙ И ПЛАСТИНОК С ВЫРЕЗАМИ

ИПЛАСТИНОК С ПОДКРЕПЛЕННЫМИ ВЫРЕЗАМИ

5.1.ОРТОТРОПНЫЕ ПРЯМОУГОЛЬНЫЕ ПЛАСТИНКИ

СВЫРЕЗАМИ ПРЯМОУГОЛЬНОЙ ФОРМЫ

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

Рассматриваем пластинки, выполненные из стекловолокнисто­ го армированного анизотропного материала (СВАМ ), считая его ортотропным материалом, так же, как это принято в работе [92]. Упругие свойства тако-го материала характеризуем четырьмя неза­ висимыми величинами: модулями упругости В\ и £ 2-по двум взаим­

но перпендикулярным направлениям х и у, модулем сдвига G и ко­ эффициентом Пуассона р, отвечающим поперечной деформации вдоль оси у.

С помощью импульсивной функции нулевого порядка можно записать переменные параметры жесткости сплошной пластинкимодели, эквивалентной перфорированной конструкции. Предпола­

гаем, что в последней имеется k прямоугольных вырезов,

контуры

которых

расположены

в

пределах х н < * < * 2г! У н < у < У 2и

t= 1,

2 ,

,, k. В этом случае

 

 

 

 

 

 

 

 

A = -EIOY (* .

Л1

Е 2= Е 2^ ( х ,

у)\ G = G y (x , у)у_

 

(5 .1 )

;где

Ею,

E 2Q и G — .модули упругости

реального

материала

пла­

стинки;

 

 

 

 

 

 

 

 

_

ft

 

 

*

[1 у

Уц)

 

 

У (-*1

У) — 1

 

У Уц)~\~^^*^0 (-Я

 

 

 

/=»1

 

 

1*-|

 

 

 

 

 

fr

 

ft

 

 

 

 

 

 

+ У г 0( * - * 1, ; ’у —

 

 

 

(5 -2 )

1-1