автореферат диссертации по , 05.00.00, диссертация на тему:Исследование устойчивости плоских пластинчатых систем

ШНЙСТЕРСТВО ВЫСШЕГО И СРЕДНЕГО СПЕЦИАЛЬНОГО ОБРАЗОВАНИЯ РСФСР

Тюменский индустриальный институт

ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПЛОСКИХ ПЛАСТИНЧАТЫХ СИСТЕМ

Диссертация

на соискание ученой степени кандидата технических наук

Научные руководители:

кандидат технических наук, профессор А.Н.КОСУХИН

кандидат технических наук, доцент В.И.КУЧЕИЖ

В.А.ЧУ РИЛОВ

Тюмень - 1975

г

ОГЛАВЛЕНИЕ

Стр.

ВВЕДЕНИЕ. ^............... 4

ГЛАВА I. РАСЧЕТ ИЗОЛИРОВАННЫХ ПРЯМОУГОЛЬНЫХ ПЛАОПШ НА УСТОЙЧИВОСТЬ ВАРИАЦИОННЫМ МЕТОДОМ В.З.ВЛАСОВА ....... 17

1.1. Основное дифференциальное уравнение 17

1.2. Фундаментальные балочные функции • • 21

1.3. Определение наименьшей критической 27 нагрузки ..............

1.4. Сходимость и оценки сходимости метода Канторовича-Власова в задачах устойчивости пластш*.........46

1.5. Пример расчета......................51

1.5,1. Алгоритм решения задачи на ЭВМ 54

ГЛАВА П. ПРИЛОЖЕНИЕ МЕТОДА ПЕРЕМЕЩЕНИЙ К

РАСЧЕТУ УСТОЙЧИВОСТИ ПЛАСТИН С ОСОБЕННОСТЯМИ В ОДНОМ НАПРАВЛЕНИИ. . • • 59

2.1. Основная система канонических уравнений метода перемещений и условие устойчивости. . ........... 59

2.2. Определение коэффициентов системы уравнений............. 62

2.3. Пример расчета......... . • 67

ГЛАВА Ш. ПРИЛОЖЕНИЕ ИОДА ПЕРЕМЕЩЕНИЙ К

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

3.1. Система канонических уравнений метода перемещений и условие устойчивости ................ 70

3.2. Определение коэффициентов системы уравнений

3.2.1. Определение реактивных усилий 74

3.2.2. Разрешающая система уравнений 77

3.3. Пример расчета........... 87

Стр.

ГЛАВА 1У. ЭКСПЕРИМЕНТАЛЬНОЕ ИССВДОВАНИЕ

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

4.1. Определение экспериментального значения критической нагрузки .... 94

4.2. Методы исследования деформаций и

'' напряжений, основанные на даровом

эффекте.......„............99

4.3. Определение аппроксимирующей функции

, формы потери устойчивости ..... НО

4.4. Примеры выполненных исследований» ♦ Н8

ЗАКЛЮЧЕНИЕ .................. 127

ЛИТЕРАТУРА .................. 129

т

л

ВВЕДЕНИЕ

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

Проблема устойчивости пластинчатых систем, находящихся под действием сжимающих внешних сил, насчитывает уже более 200 лет. В литературе по устойчивости пластин исключительно важное значение имеют работы русских и советских ученых. Можно сказать, что исследования С.П.Тимошенко, И.С.Бубнова, Б.Г.Галеркина, П.Ф.Папковича, А.А.Ильюшина и др. заложили основы устойчивости пластин, как раздела математической физики. В то же время ими решены многие практические задачи. К числу трудов, в которых сконцентрированы результаты различных исследований по устойчивости пластинок, относятся работы П.Ф.Папковича, А.Н.Динника, И.Я.Штаермана и А.А.Пиковского, Б.М.Броуде, А.С.Вольмира, А.Р.Рзканицына, В.В.Болотина и др.

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

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

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

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

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

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

сандрова и др. В указанных работах для изучения отдельной пластинки использовались решения Файлона (плоская задача) и М.Леви (изгиб пластинки), а условия сопряжения вдоль линий контакта пластин записывались на основе метода начальных параметров, сил и перемещений.

В задачах устойчивости указанное направление нашло отражение в работах А.С.Локшина (метод начальных параметров), В.А.Постнова, Ф.Блейха (метод сил). В.А.Александрова, В.Й.Климанова (метод перемещений) и др. В ряде работ П.Л.Амиро [б] , П.П.Абовского [х] » Ф.Бпейха [15], А.И.Косухина [45], Л.Анга и И.Ньюмака [98] и др. при стыковке пластин учитывались лишь изгибные деформации, причем условия сопряжения записывались как в форме метода сил, так и в форме метода перемещений.

В последнее время развитию метода перемещений способствовали работы А.М.Черняка, И.Э.Гольберга, В.И.Жесткого, Л.И.Хазовой и др., относящиеся к проверке прочности, устойчивости и изучению свободных колебаний пространственных систем, состоящих из прямоугольных изотропных и ортотропных пластин.

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

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

Проблемой устойчивости прямоугольных пластин, в том числе ослабленных прямоугольными вырезами, занимались А.А.Фельдман [эх] Г.П.Зиненко [зб], А.Баратов [ю^ и др. Их решения для пластинок с различными вариантами граничных условий были основаны на применении метода конечных разностей.

А.Н.1*узь [зз],подчеркивая общность задач устойчивости пластин и оболочек, ослабленных отверстиями, в связи с неоднородностью докритического напряженного состояния, подразделяет их на задачи трех типов, принимая в качестве критерия соотношение геометрических размеров пластины и выреза:

а) задачи первого типа: отверстия и вырезы имеют размеры одного порядка с минимальными внешними размерами пластин. Такие отверстия и вырезы существенным образом влияют на величину критической силы и форму выпученной поверхности. Эти задачи наиболее сложные и требуют для решения громоздких математических методов. Имеется немного исследований устойчивости прямоугольных пластин с большими прямоугольными отверстиями (решение проведено в основном методом конечных разностей);

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

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

Для исследования устойчивости изолированных элементов - пластинок различными исследователями применялись различные методы, которые можно свести к четырем основным:

а) методу тригонометрических рядов;

б) методу конечных разностей;

в) методу конечных элементов;

г) вариационно^ методу.

Метод тригонометрических рядов для определения критической нагрузки при расчете прямоугольных пластин на устойчивость применялся П.М.Огибаловым [бз] , С.П.Тиможенко [8б] , Р.А.Римским [72] , В.И.Климановым [44] , В.С.Ширмановым [э?] и др.

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

Метод конечных разностей применительно к задачам устойчивости пластин разработан в трудах П.М.Варвака [25] , А.Ш.Боженова и К.А.Турсунова [18] , А.А.Фельдмана [91] , А.Баратова [ю] , [12] , Ю.И.ВДузыченко [бб] и др.

Метод конечных разностей (метод сеток), в развитие и приложение которого значительный вклад внес П.М.Варвак, обладает несомненными достоинствами, но для задачи устойчивости метод не нашел пока еще широкого применения в связи с тем, что определение необходимых характеристик пластин в основной системе (перемещений и реакций на линиях сопряжений) связано с обращением матриц коэффициентов полной системы конечно-разностных уравнений для каздой пластины. Эти матрицы весьма высоких порядков, а методика расчета,включающая подбор критических параметров, требует их многократного обращения [эв] .

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

Методом сеток в уточненной постановке Ю.И.Лф'ЗЫченко [57] исследовал шарнирно опертую квадратную пластину, ослабленную квадратным отверстием, В результате точность удовлетворения граничным условиям не стала зависеть от густоты сетки.

Наряду с методом конечных разностей в последнее время интенсивно развивается метод конечных элементов.

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

Несмотря на целый ряд преимуществ, метод конечных элементов требует большой вычислительной работы.

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

Обзор способов решения задач устойчивости с использованием метода конечного элемента дан в работе Р.Галахера [юв] . Исследование устойчивости пластин МКЭ провел В.Г.Вдохов [1б] .

Метод конечных элементов в уточненной постановке использовал В.Г.Налоев [58] для решения ряда задач устойчивости прямоугольных пластин.

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

двух противолежащих краях пластины.

Примерно одновременно В.З.Власов и Л.В.Канторович предложили методику получения разрешающей системы уравнений с помощью приведения уравнений в частных производных к системе обыкновенных дифференциальных уравнений.

По сравнению с известными вариационными методами Еубнова-Гадеркина и Ритца-Тимошенко метод Канторовича-Власова обладает двумя существенными особенностями, делающими его весьма эффективным для решения поставленной задачи.

Во-первых, при решении уравнений в частных производных по методу Канторовича-Власова в качестве неизвестных принимаются не числовые параметры, а функции, зависящие от одной из двух независимых переменных. Благодаря этому, точность получаемых решений оказывается более высокой, а их сходимость более быстрой, чем это имеет место в вариационных методах Е*убнова-Галеркина или Ритца-Тимошенко.

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

Задача определения критических нагрузок аналитически в задачах устойчивости пластинчатых систем с достаточно высокой степенью точности сложна и громоздка, каким бы методом мы ни пытались ее решить»

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

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

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

Именно такому теоретико-экспериментальному методу, на наш взгляд, предстоит еще сказать свое слово,

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

Экспериментальные исследования устойчивости пластин с вырезом провел В.Г,Налоев [59] .

Б.М.Аллахвердов и др. [б] исследовали устойчивость анизотропных пластин ♦

Экспериментальные исследования напряженного состояния прямоугольной пластины после потери устойчивости провел В.Вальцак [l04] . Исследования были проведены при помощи тензометров омического сопротивления. Доказано достаточно точное совпадение результатов экспериментов с результатами теоретических вычислений, проведенных предварительно автором.

М.Скерли [ios] определял критическую нагрузку потери устойчивости квадратной пластины с круговым отверстием под действием равномерных сжимающих нагрузок. В работе приведены графики зависимости коэффициента устойчивости от размеров отверстия для двух видов закрепления сторон.

Вывод автора об увеличении жесткости пластины в связи с наличием отверстия по сравнению с отсутствием последнего при определенных значениях нагрузок и соотношениях размеров пластины и отверстия совпадает с результатами расчетов А.А.Фельдмана [91] Шлэк А.Л. в работе [Юб] изложил результаты экспериментального исследования устойчивости свободно опертых квадратных

пластинок с круговым центральным отверстием. Опыты проводились на пластинках с различным отношением диаметра отверстия к длине стороны пластинки.

И.Хейманн [107] провел экспериментальное исследование плосконапряженного состояния треугольной пластины с вырезом и без методом муаровых полос в сочетании с методом фотоупругости.

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

Метод муаровых полос получил развитие в работах И.П.Сухарева [82] , Б.И.Ушакова [90] , Б.П.Соколова [102] , Р.А.Дульнева [lOO] , В.В.Новицкого [62] , Е.Н.Андреевой [7] и др. Существенный вклад в развитие этого метода внесли П.Теокарис [*85j , Д.Дю-релли и В.Парке [33] , рассмотревшие решение значительного числа задач динамики, пластичности и изгиба пластин.

Экспериме