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

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

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

Глава 3

УСТОЙЧИВОСТЬ ПРЯМОУГОЛЬНЫХ п л а с т и н о к С ВЫРЕЗАМИ КРУГОВОЙ ФОРМЫ

3.t. ПЛАСТИНКИ С ЦЕНТРАЛЬНЫМ КРУГОВЫМ ОТВЕРСТИЕМ

Задача о потере устойчивости свободно опертой квадратной пла­ стинки без отверстия при задании равномерных краевых перемеще­ ний в ее плоскости исследовалась С. П. Тимошенко и К .Маргером [94]. Они получили выражение для определения критического пере­ мещения сплошной квадратной пластинки со сторонами 2а. Оно имеет вид

 

 

я2А2

(3 .1 )

 

ик0= ----------------

 

кр

12(1+|л)<1

 

или

« кр= 0 ,6 3 2 6 7 /£ 2/<г

(3 .2 )

при

ц=.0,3. Этот результат

далее взят нами за основу

сравнения

с решениями, .полученными для пластинок с вырезами.

 

Одной Из первых работ, посвященных исследованию устойчиво­ сти прямоугольных пластинок с центральным круговым вырезом, была опубликованная в 1947 г. работа Леви, Воллея и Кроля [111]. В ней приведены результаты изучения поведения квадратных пла­ стинок с центральным круговым вырезом. Полученные в [111] за­ висимости были использованы практически лишь в 1951 г. Кумаи [110] для сопоставления со своими опытными и теоретическими данными. Работы [11С‘, 111] посвящены изучению поведения пря­ моугольных пластинок с круговым центральным отверстием в слу­ чае равномерного нагружения вдоль двух параллельных противо­ положных краев.

Рассмотрим такую же задачу, однако исследование будем про­ водить аналогично тому, которое было выполнено А. Л. Шлаком [119, 12С]. Он изучал влияние отдельных членов в функциях, ап­ проксимирующих перемещение, на окончательное значение крити­ ческого параметра нагружения. Теоретический анализ А. Л. Шлэк осуществлял на основе энергетического метода (метод Ритца). В результате им получено простое аналитическое выражение для определения критического параметра свободно опертой прямо­ угольной пластинки с центральным круговым отверстием при одно­ осном равномерном краевом перемещении. Рассмотрим этот же случай. Считаем, что прямоугольная пластинка имеет стороны дли­

ной 2а, 2Ь и центральное круговое отвер­ стие радиусом г. На рис. 3.1 показаны раз­ меры пластинки и область интегрирования S, занимающая 1/4 общей поверхности пла­ стинки.

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

Рис. 3.1

ческим центром прямоугольной пластинки.

Функции, аппроксимирующие нормаль­ ный прогиб w и перемещения в срединной поверхности пластинки и и v соответственно вдоль осей х и у, в

общем случае можно представить в следующем виде:

(1 - х 2/а2) (1 — у2/# 2) 2

А тп (х/а)т (у/Ь)п;

 

u = B x-{-h

(1 — * 2/а 2) (х/а) ^

А ш {х[а)т- {у}Ь)п\

(3 .3 )

v = C y + h

(1 — у2/* 2) (yjb) 2

Стп {х/а)т (у/Ь)п.

 

Однако аппроксимация перемещений функциями в такогл виде не учитывает особенностей в центре пластинки, появляющихся из-за концентрации напряжений. Перемещения и н о в форме (3.3) учитывают то обстоятельство, что наружные края пластинки пред­ полагаются свободно скользящими в опорах в ее плоскости.

Прежде чем рассмотреть устойчивость пластинки, изображен­ ной на рис. 3.1, определим необходимое число членов в функциях (3 .3 ), аппроксимирующих перемещения и, v и w для сплошной пластинки. Получающееся при этом выражение для критического краевого перемещения будем сопоставлять с решением (3 .2). При достижении определенной точности совпадения (в пределах 1% ) считаем, что функции, аппроксимирующие перемещения и, v и wr при которых достигнуто это совпадение, могут служить далее на­ чальной точкой для -приближенного исследования устойчивости пластинки с отверстием.

Рассмотрим сплошную упругую прямоугольную пластинку, на краях которой вдоль оси х задано равномерное перемещение. Оп­ ределим критическое значение краевого перемещения м1{р, при ко­ тором произойдет потеря устойчивости пластинки. Для анализа воспользуемся энергетическим методом, который особенно полезен, когда точное решение неизвестно и требуется найти только при­ ближенное значение для критических параметров нагружения. Кроме того, преимущество энергетического метода состоит еще и в том, а это особенно важно для исследования поведения тонкостен­ ных конструкций с отверстиями, что он не требует выполнения обя­ зательного условия точного метода вычисления критических пара­ метров нагружения при решении дифференциальных уравнений криволинейной формы равновесия. Суть этого требования — в не­ обходимости удовлетворения заданным краевым геометрическим я

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

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

Энергетический метод основан на рассмотрении потенциальной энергии выпучившейся пластинки.

Детальное обоснование нахождения критического значения на­ грузок с помошыо энергетического метода дано в работе П. Ф. Папковича [82].

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

равновесия в криволинейную равно:

 

 

 

 

д П = и + У —W,

 

 

(3.4)

где U — потенциальная энергия

изгиба и кручения

пластинки;

V — изменение потенциальной энергии деформации срединной по­

верхности пластинки при выпучивании; W — изменение потенциала

внешних сил, приложенных к пластинке, или же

потенциальная

энергия нагружения.

 

 

 

 

 

Потенциальной энергией деформации пластинки поперечными

силами Qxz и Qxv пренебрегаем из-за ее малости.

 

 

 

Потенциальная энергия изгиба и кручения

U

прямоугольной

пластинки с размерами сторон 2а и 2b равна:

 

 

 

и = _ ^ г г

2 ( 1 _

 

 

 

6(1 — f*2) J И дхР

fry* )

r L дх2 dyfi

V дхду

) ]|

 

 

 

 

 

(3 .5 )

Для нахождения V необходимо выписать выражения для линей­ ных и угловых деформаций срединной плоскости, обусловленных ее выпучиванием [82].

При одновременном наличии перемещений и и w искомая ли­

нейная деформация в направлении оси х :

 

еX

( 3 .6а)

Перемещения v не влияют на деформацию е*. Диалогично линей­ ная деформация срединной поверхности в направлении оси у:

(3.66)

При одновременном наличии перемещений а, и и w искомая уг­ ловая деформация срединной поверхности:

ди , dv , dw dw

Изменение потенциальной энергии деформации срединной по­ верхности пластинки с размерами сторон и 2b выражается че­ рез деформации следующим образом:

V = i ^ \ \ (е‘ + 4 + 2 н * А + - Ц ^ У%) dxdy.

(3 .7 )

Используем соотношение (-3.4), (3.5) и (3.7) для определения кри­ тических параметров нагружения. Применительно к нашему слу­ чаю требуется определить одноосное критическое перемещение внешних противоположных границ в направлении оси х. Критиче­ ское значение одноосного перемещения границ и* — это перемеще­ ние, при котором наряду с плоской формой равновесия пластинки возникает новая криволинейная форма равновесия.

Общим признаком равновесия материальной

системы является

экстремальность полной потенциальной энергии

П

системы.

В рассматриваемом нами случае равновесия

выпучившейся

пластинки полная потенциальная энергия:

 

 

П — П0-1~ ДП,

 

(3 .8 )

где П о— потенциальная энергия пластинки до выпучивания (не за­ висящая от амплитудных параметров прогиба Атп)\ ДП — изме­ нение потенциальной энергии при выпучивании, определяемое в виде

Д П = V — W .

(3 .9 )

Из условия экстремума энергии П следуют уравнения:

— = 0 = — 4 - — — —

(3. 10)

dXt ~ ~ d X i dXi dXt

где Xi — амплитудные параметры перемещений в аппроксимирую­

щих функциях (3 .3). Уравнение

(3.10) ведет

к целому ряду сов­

местных уравнений типа

ап

 

ап

 

0;

дБ, ■ = 0;

0;

(3 .1 1 )

дА/пп

тп

дСтп

 

Полученная таким образом система линейных однородных уравне­ ний (3.11) относительно неизвестных произвольных постоянных А тп, Втп, Стп . . имеет нулевые решения Атп = Втп= С тп= . . . =

=0 , которые соответствуют плоской форме равновесия пластинки,

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

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

Если члены ряда аппроксимирующих функций

(3.3)

совпадают

с истинным уравнением срединной

поверхности

пластинки

(при

весьма малых перемещениях), то

энергетический

метод

дает

точ-

ное значение критических параметров нагружения. В противном случае, что практически и имеет место [22], энергетический метод приводит к критическим параметрам, несколько превышающим их действительное значение.

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

® = Л 00/г (1-- х ^ а ? ) (1 У2/ Ь-)\

и = В х\

(3 .1 2 )

•и=0.

Подставляя (3.12) в выражение для потенциальной энергии изгиба (3.5), найдем

 

 

 

1

а )

■СцАоо,

(3. 13)

1 — (л*

1 5 а3

3 ab

~Fi

 

 

16

£А5

ь

 

 

 

(3 .14)

где

1 — [*2

{"5а3

 

 

 

9

 

 

 

 

Аналогичную подстановку проделаем для потенциальной энергии

деформации срединной плоскости (3.7). В итоге получим

 

V =

СцВ2 -f- Ci2^Aoo-}- ^i3-4oo»

(3 .15)

где

c u = 2£А

ab\

 

 

1 —

(х2

 

 

 

 

 

 

<ЗЛ 6>

 

1 0 2 4 _ ^ f _ /A

 

М

 

13

1575 1 — ^ 2 ^ 3

1

$3 1 7ab)

 

Рассмотрим потенциальную энергию нагружения.

Положим, что одна из кромок, например х = — а, неподвижна; при этом получим, что взаимное смещение кромок х — —а и х = + а должно иметь фиксированное значение. Обозначая ех как относи­ тельное сближение краев, можем записать:

е х = —й (^х -+д ^ х — а } '

( 3 . 1 7 )

Учитывая (3 .6), получим

ifh-T(^)2]dA:=const (ЗЛ8)

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

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

ь

P x = 2 h

Г

i_axdy для 0 < х < г ;

(3.19)

 

'

2

 

 

Ь

 

 

P x = 2 h

I* axdy для г < л ; < а ,

(3.20)

 

б

 

 

где

 

(£.v+ fte'/)-

 

Потенциальная энергия нагружения

 

 

 

<•-+-■> [---т е

dxdy.

(3.21)

 

 

о о

 

 

 

 

 

 

Подставляя в (3.21) необходимые соотношения из

(3.6) и функция

перемещения в форме (3 .12), получим:

 

 

 

 

 

 

W =

сзг5 2-(- С32ДЛ0О)

 

 

(3. 22)

где

 

 

ся =

1^а|х2Ь

 

 

(3. 23)

 

 

 

 

 

 

 

 

 

 

32

В »

( » +

f JL\.

 

 

 

32

45

1 — ц2( а

Ь )

 

 

 

Д ля заданного относительного сжатия ех

пластинки

неизвест­

ные постоянные В и А 00 в выражениях

(3.12)

можно найти из ус­

ловия, чтобы полная энергия деформации П

была

минимальной;

следовательно, учитывая соотношения (3 .13),

(3 .1 5 ), (3.22) и усло­

вие экстремума

(3 .10), имеем:

 

 

 

 

 

——= 0 — 2 В{Сц — ^31) —|—Лоо {С\2

£32)*

 

 

 

ов

 

 

 

 

 

 

(3 .24)

 

 

 

 

 

 

 

 

-----= 0 =

В (Cj2 — С32) + 2^00^13 "Ь ^21*

 

 

V-AQO

 

 

 

 

 

 

 

Реш ая систему

(3.24)

относительно неизвестных постоянных Л0о i

В, найдем, что

 

 

 

 

 

 

 

 

 

В = ______ £2i (gi2032)------------

 

(3. 25]

 

 

(<?ц — С31) — (<<12— С32)2

 

 

Д ля квадратной пластинки со сторонами, равными 2а, с помощью (3 .1 4 ), (3.16) и (3.23) можно установить, что

Выражение для критического перемещения «кр определим из соотношения

а

йкр= — ^ udx.

(3 .27)

о

Используя (3.12) и (3.25), из (3.27) получим

икр= li*h 2ja,

(3 .28)

где k * =0,70513.

Проведем аналогичные выводы для различных комбинаций Л тп,

используя для этого выражение для нормального прогиба до в фор­ ме (3.3), получим различные числовые коэффициенты в формуле для критического перемещения (3. 28).

 

 

Т а б л и о а 3.1

Коэффициенты членов,

 

 

учтенных в выражении

к*

Д, %

прогиба

 

 

Аоо

0,70513

11,45

AQO, A>2Q

0,67157

6,15

Л о о » А0 2

0,66685

5,40

AQOI ^ 2 0 » Л < )2

0,63555

0,46

Результаты такого исследования представлены в табл. 3.1. По­ лученные результаты иллюстрируют сходимость решения к точно­ му значению (3.2) в зависимости от вида функции, аппроксимиру­ ющей нормальный прогиб. Наибольшая точность (в пределах 1% ) получилась, когда в функции для до были использованы три члена с коэффициентами А00, А20 и Лог-

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

Однако до перехода к их рассмотрению целесообразно исследо­ вать вклад отдельных членов, используемых в функциях, аппрокси­ мирующих перемещения, в окончательный результат. Этот вопрос изучался различными авторами на первоначальном этапе поиска пригодных решений для прямоугольных пластинок с отверстиями. Так, Рао и Пикетт в работе [117] предложили ввести в аппроксимиру­ ющие функции член в виде логарифмической функции в дополне­ ние к тем, с помощью которых аппроксимировались перемещения для оплошной пластинки, чтобы получить более точное решение для пластинки с отверстием. Этот член использовался ими в зада­ че о колебаниях пластинки с отверстием. Принимая во внимание математическую схожесть в решениях задач по устойчивости с за ­ дачами по колебаниям, А. Л . Шлэк [119] по аналогии с работой

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

® = А ( 1 “

р ) [ А т + А 10' й

<а29а)

Так как учет перемещений и и v дает значительно

меньший вклад

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

u = B x + k (l -

-ff) -f- [Зх> + B *> ^ -\-В<,г^ + а Ь В а

 

 

 

 

 

 

(3 .2 9 6 )

v = C y -\ -h (l — — ) —

ГСоо+Сщ — + C 02— +<г*С0

^ -1 .

\

Ь-2 I b

\

'

(*2 +

02)2 j

 

 

 

 

 

(3. 29 B)

Сложность математических выкладок, встречающихся при ис­

следовании задачи на основе функций

смещения

в виде

(3 .29а),

(3.296) и (3 .29в ),

делает

использование

системы

(3.29) в качест­

ве решения весьма ограниченным [117]. В связи с этим А. Л . Шла­ ком использованы упрощенные приближенные решения в предпо­

ложении, что пластинка перфорирована лишь одним малым

(г /а <

< 0 ,3 ) центральным отверстием. Выводы осуществлены с

двумя

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

Вначале было найдено решение при использовании смещений и и о в форме (3.296) и (3.29в), функция же нормального прогиба w была упрощена до вида

( а 3 0 )

Предполагалось, что член с коэффициентом А\, взятый по анало­ гии с работой [117], должен учитывать наличие в пластинке от­ верстия. Д ля оценки его вклада вывод был проведен с учетом и без учета этого члена при одинаковых зависимостях для а и v. Резуль­ таты сведены в табл. 3.2, где дано сравнение числовых значений для критического параметра k*, позволяющего по формуле (3.28) определять величину критического одноосного омещения (С ==0) для пластинки с центральным круговым вырезом. Здесь же в про­ центах показано различие в результатах, определяемое по разни­ це в коэффициентах, в сравнении с решением при наличии членов с А т и А\.

 

Параметры, характери­

 

Значения у —г/а

 

 

 

зующие нормальный

 

 

 

 

прогиб по соотношению

0,01

0,1

0,3

 

 

(3.30)

 

 

 

 

^оо^О; А\^=0

0,70469

0,67а1э

0,62803

 

Доот^О; i4i = 0

0,70516

0,70888

0,75952

Расхождение

 

k\-k\

 

 

 

- --■ioo%

0,07

4,67

17,31

результатов

*

Вычисления выполнены для квадратных пластинок со стороной 2а и центральным круговым вырезом с радиусом г, характеризуе­ мым параметром y = rja .

Из табл. 3.2 видно, что логарифмический член с коэффициентом А\, дополняющим обычно используемый коэффициент Доо, влияет гораздо сильнее при увеличении размеров отверстия. Имея точное значение критического коэффициента k *= 0 ,6 3 2 6 7 для сплошной пластинки и сравнивая его с &2* из табл. 3.2, приходим к выводу, что упрощенные аппроксимирующие функции дают завышенные значения критического параметра для пластинок с очень малень­ ким отверстием, радиус которых близок к нулю, например, для 7 = 0 ,0 1 .

Второй тип функций смещения, на основе которых определялась величина критического одноосного перемещения противоположных торцов вдоль оси х ( С = 0 ) , имеет вид

и = В х + Л (1 - ^

f

[ад -Э ю у

+

-£];

(3.31)

v = C y + h(\ -

у

[Сю+ С 20у

+ с м

у }

 

На основании выведенных зависимостей получены числовые значения коэффициента k*, входящего в формулу (3.28) для крити­ ческого осевого смещения пластинки с центральным круговым от­ верстием. Результат таких вычислений:

у

о

0,1

0,3

k*

0,63555

0,62731 0,59335

Сравнение численных значений

для функций перемещения с

необычным (логарифмическим) членом (пунктирная кривая) и без него (сплошная кривая) дано на рис. 3.2. Кривые построены по результатам расчетов величины k* для пластинок с различным

 

 

значением у. Из их анализа видно, что функ­

 

 

ции смещения в форме (3.31) дают лучшие

 

 

результаты. Если диаметр отверстия стре­

 

 

мится к нулю, то решение приближается к

 

 

точному. Расхождение

в

этом случае со­

 

 

ставляет 0,5 %. Такой вывод для пластинок

 

 

с отверстием

основан

на

предпосылке О'

Рис.

3.2

том, что функции, использованные согласно-

выбранному методу исследования, дают за­

 

 

 

 

вышенные значения критического парамет­

ра. Метод

Ритца дает приближения

сверху

[117], сближающиеся,

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

w = - x 2Ja2) h — ~ ) Л о;

и В х —{—Л (1 — x^ja?) (х/а) [5 00-j- B 2QX21а2 -f- BQ2y2/b2); (3. 32)

v=C fif+A (1 - fib1) Ш [C00 + С20х2(а2-f- C02y2jb2].

Уточнение решения достигается прибавлением к решению (3.32) добавочных членов, например, типа логарифмического члена, ис­ пользованного в (З.ЗС).

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

На практике обычно за критерий потери устойчивости прини­ мается критическая нагрузка. В связи с этим следует сделать пере­ ход от критического перемещения к критической нагрузке. Чтобы сделать этот переход, необходимо напряжение вдоль края пластин­ ки выразить через функции перемещения'и затем проинтегрировать эти напряжения по площади пластинки. Используя зависимосги (3.19) и (3 .2 0 ), можно определить значение нагрузки как среднее значение Рх по длине пластинки. Интегрируя, получаем

где

Значение [P z]ср, соответствующее критическому перемещению мКр»- определяемому по зависимости (3 .28), будем считать критическим и обозначим его через Р*. С учетом двух последних формул

р * = 2 ё 1 .k xu * =

- 2Ef?

. k2,

1 — f i 2 1

a ( l —

, i 2 ) .

где