книги / Устойчивость и колебания пластинок и оболочек с отверстиями
..pdfГлава 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
Изменение потенциальной энергии деформации срединной по верхности пластинки с размерами сторон 2а и 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 (gi2— 032)------------ |
|
(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 ) . |
где