Full Text

Введение Задача о неосесимметричной потере устойчивости круглой пластины под действием радиальных сжимающих усилий, равномерно распределённых по её контуру, впервые поставлена и решена в работе Дж. Брайена (G. Bryan) 1891 г., в которой при выводе дифференциального уравнения для функции прогиба и граничных условий использован вариационный подход [1]. Там же получено первое критическое значение нагрузки для случая потери устойчивости по форме с одним узловым диаметром. Обзор последующих работ в этом направлении содержится в статье [2], монографии [3] и справочнике [4], где аналогичное исходное уравнение получено из условий равновесия для внутренних силовых факторов пластины. Ряд вопросов устойчивости прямоугольных и круглых пластин, в том числе при температурном воздействии, освещён в [5]. Среди современных работ данного направления отметим работу, посвященную анализу устойчивости круглых пластин при неравномерной сжимающей нагрузке, распределённой по окружности по косинусоидальному закону с изменением параметра от 0 до 1 (в последнем случае наблюдается переход к балочной форме потери устойчивости) [6]. В [7] исследованы формы потери устойчивости круглой пластины, частично опирающейся на упругое основание. Здесь обнаружена зависимость окружного числа волн в форме потери устойчивости пластины от жёсткости упругого основания. Эффект перехода от осесимметричной к неосесимметричной форме потери устойчивости круглой пластины с ростом действующего на неё равномерно распределенного давления исследован в [8-10]. Для описания этого эффекта в [8] выведена система уравнений равновесия несимметричных пологих оболочек в больших перемещениях. Решение этой системы отыскивалось методом Бубнова-Галёркина при задании формы потери устойчивости пластины радиусом а в полярных координатах в виде (1) где А и В - амплитуды осесимметричной и неосесимметричной складчатости формы потери устойчивости пластины; m - заранее неизвестное целое число, определяющее количество складок при обходе пластины по кругу при условии минимума интенсивности нагрузки для заданного уровня прогиба центра пластины. В результате было показано, что при закреплении края пластины, не препятствующего проскальзыванию в окружном и радиальном направлениях, неосесимметричная потеря устойчивости с минимальным значением критической нагрузки происходит по форме, имеющей 8 волн по окружности пластины. В [9] для определения критических нагрузок в неосесимметричных формах потери устойчивости использован метод малого параметра, с помощью которого минимальное значение сжимающей краевой критической нагрузки реализуется на форме, имеющей 6 волн по окружности пластины. Аналогичный метод использован в [10] в более общем случае, когда модуль упругости круглой изотропной пластины зависит от радиальной координаты. Исходное представление (1) для формы потери устойчивости использовано и в [11] при решении несколько иных уравнений, полученных как частный случай из уравнений для пологой сферической оболочки. Здесь также выведены условия неосесимметричной потери устойчивости пластины, не подкреплённые, правда, расчётным примером. Доказательство существования несимметричного решения при больших прогибах симметрично нагруженной пластины для конкретного вида нагрузки и нулевого коэффициента Пуассона дано в [12]. Наиболее подробный обзор, охватывающий в целом состояние работ по осесимметричной и неосесимметричной потере устойчивости пластин, содержится в [4], где, в частности, показано, что в тонких пластинах, подверженных изменениям температуры или влажности, возникают локальные сжимающие напряжения, которые способны вызвать продольный изгиб даже в отсутствие поперечной нагрузки. Среди экспериментальных результатов отметим испытания круглых пластин на устойчивость под действием равномерного давления, также обнаружившие неосесимметричную складчатость вблизи контура пластин [13]. Ниже представлены результаты экспериментов, показывающие краевую потерю устойчивости с большим числом волн по окружности пластины при изменении её температуры. Продемонстрированы возможности согласования классической теоретической модели [1] с наблюдаемыми формами потери устойчивости. 1. Экспериментальное проявление многоволновой краевой потери устойчивости круглой пластины Для проявления многоволновой краевой потери устойчивости круглой пластины проводился следующий эксперимент. В качестве объекта эксперимента использовались стандартные пластиковые крышки от упаковки пищевых продуктов, изготовленные из полипропилена (модуль упругости Е = 1,15 ГПа, коэффициент Пуассона n = 0,35, температура плавления 160 ºС) с диаметром плоской части 80 мм. Эти крышки надевались на обод цилиндрического сосуда (такого же диаметра) с кипятком. На части крышек в процессе нагревания образовывались неосесимметричные краевые формы потери устойчивости (рис. 1). а в г б Рис. 1. Фотографии многоволновых форм потери устойчивости круглых пластин: а - форма с 8 волнами по окружности ( - радиус окружности, проведенной через экстремумы формы, - радиус переходной линии); б, в - формы с 20 и 32 волнами по окружности; г - форма, воспроизведенная из [8] В зависимости от толщины пластин, которая изменялась для разных образцов от 0,18 до 0,4 мм, образовывались формы потери устойчивости с разным числом волн по окружности. Как правило, на пластинах минимальной толщины регистрировалась форма с большим числом волн, как на рис. 1, б, в, где видны формы с 20 и 32 волнами по окружности. Напротив, более толстые пластины теряли устойчивость по формам с 6-8 волнами по окружности (рис. 1, a). Для сравнения приведён также пример краевой потери устойчивости пластины из [8] (рис. 1, г) под действием равномерной нагрузки, не связанной с температурными деформациями. Отметим, что зарегистрированные формы многоволновой потери устойчивости похожи на формы краевых резонансов у выпуклых границ пластин и оболочек, отделённых от практически недеформированных внутренних областей переходными линиями [14]. Одна из таких линий показана на рис. 1, а штриховой окружностью радиусом . 2. Неосесимметричная краевая потеря устойчивости круглой пластины В качестве исходного уравнения для изучения потери устойчивости круглой пластины под действием нагрузок в её плоскости используется обычно уравнение [3] (2) в котором через sr, sj, t обозначены соответственно радиальная, окружная и тангенциальная компоненты плоского напряженного состояния пластины, схематично показанные для её элемента на рис. 2 (D, h - цилиндрическая жёсткость и толщина пластины). При сжатии пластины равномерно распределённой по контуру нагрузкой р возникающее в пластине напряжённое состояние, являясь частным случаем решения задачи Ламе для цилиндра (при отсутствии внутреннего отверстия) [5, 15] Рис. 2 Компоненты плоского напряженного состояния на элементе пластины толщиной h (3) позволяет свести уравнение (2) к линейному операторному уравнению четвёртого порядка [1, 3] (4) которое описывает в том числе температурное выпучивание пластины [4] (о связи нагрузки р с температурой см. следующий подраздел). Общее решение уравнения (4), регулярное в центре пластины, имеет вид [1, 3] (5) в котором функция Бесселя играет роль осциллирующей по радиусу компоненты, а описывает краевой эффект; с1, с2 - постоянные, определяемые из граничных условий на внешнем контуре пластины; m = 0, 1, 2, … - число волн в форме потери устойчивости по окружности пластины. Связь между постоянными с1, с2 и критическое значение нагрузки р находятся при подстановке выражения (5) в граничные условия. Так, в случае защемлённого края пластины получим, что значения критической нагрузки находятся через корни xm уравнения [1, 3]. Форма потери устойчивости при таких граничных условиях с точностью до постоянной имеет вид (6) где индекс m указывает на принадлежность данной формы выбранному числу волн по окружности пластины. Одна из таких форм ( ) изображена на рис. 3. Увеличение минимальной критической нагрузки Р с ростом числа волн по окружности в форме потери устойчивости пластины можно проследить по графику рис. 4, где по оси ординат отложено обезразмеренное значение Рис. 3 Расчетная форма потери устойчивости с 8 волнами по окружности пластины Рис. 4. Зависимость минимальной критической нагрузки Р от числа волн m в форме потери устойчивости пластины Внутренней границей осцилляции формы потери устойчивости является переходная линия , положение которой определяется из равенства . При функция прогиба пластины в трёхмерном представлении по радиальной и окружной координатам образует плато, что можно заметить как на графике теоретической формы потери устойчивости (рис. 3), так и на результатах визуализации экспериментальных форм потери устойчивости пластины (см. рис. 1). 3. Согласование экспериментов с расчётной моделью Наблюдаемые в экспериментах формы потери устойчивости качественно совпадают с расчётными формами (5), одна из которых отображена на рис. 3, и при отсутствии начального прогиба - с (1). Более полное согласование может быть достигнуто при совпадении таких характерных признаков, как расположение экстремумов формы потери устойчивости по отношению к центру пластины. Для нахождения этих экстремумов в выражениях (5) берётся производная по r, которая затем приравнивается к нулю. Решение получившегося неявного уравнения даёт координату экстремумов расчётной формы потери устойчивости при заданных граничных условиях по краю пластины. Так, при выбранных выше механических параметрах пластины с толщиной 0,37 мм, условиях жёсткого защемления и восьми волнах в форме потери устойчивости r1 » 30 мм; аналогичный расчёт при условиях шарнирного опирания даёт для r1 значение 33 мм. В то же время измерения на соответствующей экспериментально полученной форме показали, что этот радиус приближенно равен 32 мм. Для лучшего согласования результатов эксперимента с расчётной моделью наряду с условиями жёсткого защемления и шарнирного опирания было рассмотрено упругое опирание пластины с заранее неизвестной жёсткостью опоры. Формулировка условия упругого опирания для круглой пластины сводится к пропорциональности поперечного реактивного усилия [3]: и прогиба пластины на краю где - поперечная сила; , - радиальный и окружной изгибающие моменты, - крутящий момент, , - кривизны изогнутой поверхности в сечении вдоль диаметра и в перпендикулярном направлении; - кручение поверхности; с - коэффициент жёсткости опоры. Кроме того, при ставится условие Подстановка выражения (5) в эти условия приводит к характеристическому уравнению, из решения которого находится зависимость между параметром критической силы и коэффициентом жёсткости опоры . Полагая затем координату экстремумов прогиба совпадающей с экспериментальной, находим значения и с; в рассматриваемом примере при против для жёсткого защемления, а При этом радиус окружности, проведенной через экстремумы формы потери устойчивости, увеличился с 30 мм при жестком защемлении края пластины до 32,5 мм - при упругом. Приведенные данные иллюстрируют тенденцию сдвига экстремумов формы потери устойчивости к границе пластины при ослаблении краевых условий. Можно также оценить величину сжимающей нагрузки р при описанных выше условиях экспериментов в случае температурного расширения пластины. Для этого достаточно приравнять температурное изменение размеров круглой пластины при её нагревании на Т ºС к радиальному перемещению, являющемуся решением соответствующей задачи Ламе для сплошной пластины при действии сжимающей нагрузки р в плоскости пластины по её внешнему контуру. В результате , где at - температурный коэффициент линейного расширения материала пластины. Подобная формула для случая пластины, равномерно нагретой вместе с массивной обоймой, сделанных из материалов, имеющих разные температурные коэффициенты линейного расширения, приведена в [5]. Осесимметричная температурная задача устойчивости круглой пластины при местном нагреве рассмотрена в [16]. Задавшись значением коэффициента температурного расширения полипропилена at = 1,1·10-4 1/град [17] и температурой Т = 80 ºС, получим для величины температурной сжимающей нагрузки рс, вызывающей осесимметричный продольный изгиб пластины, значение 12,65 МПа, что намного выше максимального из рассчитанных критических значений этой нагрузки (3,17 МПа при m = 32) при многоволновой краевой потере устойчивости пластины (см. рис. 4). Таким образом, многоволновая краевая потеря устойчивости является естественной реакцией круглой пластины на действие сжимающей нагрузки при нагревании пластины в условиях граничного опирания, стесняющего её свободную деформацию. Заключение В работе представлены результаты экспериментов, показывающих возможность неосесимметричной потери устойчивости круглой пластины с большим числом волн, локализованных у края пластины при осесимметричной нагрузке, возникающей в процессе равномерного нагревания пластины. Показано, что с помощью простых аналитических моделей можно добиться достаточно полного описания экспериментально наблюдаемых особенностей эффекта краевого неосесимметричного выпучивания круглой пластины.

About the authors

R V Goldstein

Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences; Moscow Aviation Institute (National Research University)

A L Popov

Moscow State University of Civil Engineering

V M Kozintsev

Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences

D A Chelyubeev

Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences

References

Statistics

Views

Abstract - 176

PDF (Russian) - 145