Full Text

Введение Оболочечные конструкции широко применяются в различных областях техники, и их исследование важно для большого числа различных прикладных задач [1-3]. Например, оболочки, ослабленные по технологическим причинам большим числом вырезов (перфорированные оболочки), часто применяются в машиностроении и строительстве [4]. Такие конструкции подвергаются воздействию на них различных нагрузок [5, 6], в связи с этим возникает необходимость исследования их устойчивости. По сравнению с другими видами оболочечных конструкций перфорированным оболочкам посвящено сравнительно мало публикаций [4-16], и их исследование является достаточно актуальной задачей. Обзорных работ в данной области также немного [17-19]. Основная сложность при исследовании оболочек, ослабленных вырезами, состоит в необходимости решения задачи для многосвязной области. Для некоторых простых случаев возможно, применяя аппарат конформных отображений, свести исходную задачу к задаче для односвязной области. Если рассматривать оболочки, ослабленные вырезами, как оболочки ступенчато-переменной толщины [20], то в дискретно расположенных местах жесткость оболочки принимается нулевой. При этом на границе вырезов, как показано в работе [20], выполняются краевые условия свободного края, которые выполняются тем точнее, чем точнее решается сама краевая задача для таких конструкций. Главное преимущество такого подхода состоит в том, что область интегрирования уравнений равновесия становится односвязной. Существуют различные математические модели и алгоритмы исследования ослабленных вырезами оболочек. Так, в основе модели могут быть использованы уравнения в смешанной форме [4] или функционал полной потенциальной энергии [13]. В работах [6, 14, 15] исследование проводится на основе метода конечных элементов, а в работах [12, 13] - с применением метода Ритца. Перфорированные оболочки могут быть выполнены как из изотропных [4, 10, 11, 14], так и ортотропных материалов [6, 11]. В работе [4] рассматриваются статические и динамические задачи для пологих оболочек, перфорированных большим числом вырезов различной формы. Получены аналитические решения на основе сочетания асимптотического подхода и метода гомогенизации. В работе [15] проводится сравнительный анализ результатов расчета тонких упругих оболочек с вырезами с использованием моментной схемы конечных элементов и программных комплексов «ЛИРА» и SCAD. На примере пологих изотропных панелей показано влияние различных ослаблений на процесс нелинейного деформирования и потерю устойчивости. В работе [6] методом конечных элементов исследована устойчивость густо перфорированных цилиндрических оболочек под действием внешнего давления. Задача решена на основе теории оболочек с использованием конструктивно ортотропной модели. Получены коэффициенты снижения жесткости при различных значениях пористости. Проведено сравнение численно полученных результатов с аналитическими оценками Э.И. Григолюка и Л.А. Фильштинского [8]. Проведено исследование устойчивости упругой цилиндрической перфорированной оболочки под действием внешнего давления для двух вариантов граничных условий. Получены значения критических нагрузок и соответствующие им формы потери устойчивости в зависимости от длины оболочки и степени перфорации. 1. Постановка задачи Целью данной работы является анализ устойчивости перфорированных пологих оболочек двоякой кривизны и обоснование использования при проведении таких расчетов метода конструктивной анизотропии. 2. Вид рассматриваемых конструкций Будем рассматривать пологие оболочки прямоугольного плана, ослабленные вырезами, в геометрически нелинейной постановке в рамках теории оболочек Кирхгофа-Лява. Параметры Ляме и параметры кривизны пологой оболочки соответственно: Оболочка по контуру закреплена шарнирно-неподвижно, находится под действием равномерно-распределенной поперечной нагрузки и ослаблена большим числом сквозных вырезов со сторонами, параллельными сторонам оболочки (рис. 1). Рис. 1. Схематичное изображение пологой перфорированной оболочки 3. Математическая модель Математическая модель деформирования пологой оболочки, ослабленной вырезами, состоит из трех групп соотношений: 1) геометрических соотношений; 2) физических соотношений; 3) функционала полной потенциальной энергии. Геометрические соотношения в срединной поверхности оболочки, связывающие перемещения и деформации, имеют вид где - перемещения точки срединной поверхности оболочки вдоль осей соответственно. В слое, отстоящем на от срединной поверхности оболочки, деформации принимают вид где - функции изменения кривизн; - функция изменения кручения, Физические соотношения - связь напряжений и деформаций - принимают вид Здесь - модуль упругости и коэффициент Пуассона материала оболочки. Места расположения вырезов оболочки можно задать функцией где - толщина оболочки; индексы и указывают номер прямоугольного в плане выреза, стороны которого параллельны осям и ; - количество вырезов вдоль осей и соответственно; - единичные столбчатые функции, которые представляют собой разности двух единичных функций . Эти функции равны единице при а в остальных точках поверхности оболочки равны нулю. Здесь - ширина вырезов в направлении осей и соответственно. Таким образом, толщина оболочки в точке будет равна . Учитывать наличие вырезов будем двумя способами: дискретно и методом конструктивной анизотропии [13]. В зависимости от этого выражения для усилий и моментов будут иметь разный вид. Интегрируя напряжения по в пределах от до , получим вначале усилия и моменты, приходящиеся на единицу длины сечения и приведенные к срединной поверхности оболочки, при дискретном вводе вырезов: где - жесткостные параметры вырезов, При использовании метода конструктивной анизотропии усилия и моменты, приходящиеся на единицу длины сечения и приведенные к срединной поверхности оболочки, примут вид [13] где Функционал полной потенциальной энергии рассматриваемой оболочки будет иметь вид [21] При дискретном введении вырезов более подробно функционал можно записать в виде (1) где А при введении вырезов методом конструктивной анизотропии он запишется следующим образом: (2) 4. Алгоритм Для минимизации функционала полной потенциальной энергии применяется метод Ритца, при этом искомые функции аппроксимируются в виде где - неизвестные числовые параметры; - различные сочетания аппроксимирующих функций , удовлетворяющих краевым условиям при шарнирно-неподвижном закреплении контура, в качестве которых можно использовать тригонометрические функции В соответствии с методом Ритца разложения искомых функций подставляются в функционал (1) или (2), находятся производные по неизвестным числовым параметрам, которые затем приравниваются к нулю. В результате получается система нелинейных алгебраических уравнений, которая решается методом Ньютона. Также после перехода к безразмерным параметрам [22] для решения полученной системы можно применить метод продолжения решения по наилучшему параметру [23]. В качестве критерия оценки прочности принят критерий Мизеса: где - предел текучести материала; - интенсивность напряжений. Данный алгоритм был реализован в среде аналитических вычислений Maple 2015. 5. Результаты расчетов Рассматривались три тонкие пологие оболочки, квадратные в плане, входные параметры которых заданы в табл. 1. Рассматриваемые конструкции были ослаблены сквозными вырезами и находились под воздействием равномерно-распределенной поперечной нагрузки Оболочки выполнены из стали с модулем упругости МПа и коэффициентом Пуассона , МПа. Контур конструкции закреплялся шарнирно-неподвижно. Расчеты проводились при удержании членов разложения искомых функций в методе Ритца. Таблица 1 Входные параметры рассматриваемых вариантов оболочек Параметр Вариант 1 Вариант 2 Вариант 3 Линейный размер а, м 5,4 10,8 54 Главный радиус кривизны R, м 20,25 40,5 135,9 Толщина оболочки h, м 0,09 0,09 0,09 Распределение вырезов по оболочке проводилось двумя способами. В случае I (рис. 2, а) оболочка в плане разбивается на участки, число которых равно числу вырезов. Центр участка считается центром соответствующего выреза. Абсцисса i-го выреза вдоль оси имеет вид а ордината -го выреза вдоль оси - В случае II вырезы распределяются по оболочке равномерно (рис. 2, б). Абсцисса центра -го выреза вдоль оси имеет вид , а ордината -го выреза вдоль оси равна а б Рис. 2. Схема распределения 4 вырезов по оболочке Цель исследования - подтвердить сходимость метода конструктивной анизотропии. Предполагается, что сходимость происходит при большом числе вырезов. Покажем это следующим образом. Сначала фиксируется коэффициент относительной площади вырезов , после чего увеличивается число вырезов. Расчеты проводились для числа вырезов , притом рассматриваются только такие наборы вырезов, которые не содержат вырез в центре конструкции, так как анализ прогибов проводился в точке . На всех графиках будем обозначать кривую, полученную с использованием метода конструктивной анизотропии, как , а кривые, полученные путем дискретного ввода вырезов, как I и II для способов распределения вырезов I и II соответственно. Для оболочки варианта 1 был принят коэффициент . На рис. 2 показаны распределения вырезов по оболочке для способов I и II при 4 вырезах, а на рис. 3 - соответствующие кривые «нагрузка-прогиб». На рис. 4 представлены распределения 24 вырезов по оболочке (слева способ I, справа - II), а на рис. 5 - соответствующие кривые «нагрузка - прогиб». Рис. 3. Кривые «нагрузка-прогиб» для оболочки варианта 1 при 4 вырезах Результаты расчетов для оболочки варианта 1 при разном числе вырезов сведены в табл. 2. В первом столбце находится общее число вырезов , где квадратные скобки обозначают целочисленное деление; во втором - относительный размер выреза ; в третьем - критическая нагрузка потери устойчивости , полученная методом Рис. 4. Схема распределения вырезов при Рис. 5. Кривые «нагрузка-прогиб» для оболочки варианта 1 при 24 вырезах конструктивной анизотропии; в четвертом и пятом - критические нагрузки потери устойчивости и , полученные путем дискретного ввода вырезов по способам их распределения I и II соответственно; в шестом и седьмом - и Если и попадают в интервал [0,95; 1,05], то считаем, что метод конструктивной анизотропии сходится к дискретному вводу вырезов. На рис. 6 показан график зависимости относительной критической нагрузки потери устойчивости от общего числа вырезов для способов распределения вырезов I и II. Таблица 2 Результаты для оболочки варианта 1 при Критические нагрузки потери устойчивости , МПа , МПа , МПа 4 0,1667 2,7599 2,5699 2,2899 0,9312 0,8297 6 0,1361 2,5799 2,3899 0,9348 0,8659 8 0,1179 2,8799 2,6499 1,0435 0,9601 12 0,0960 2,6200 2,5400 0,9493 0,9203 16 0,0833 2,6399 2,6299 0,9565 0,9529 20 0,0745 2,7299 2,6700 0,9891 0,9674 24 0,0680 2,8299 2,7499 1,0254 0,9964 36 0,0555 2,7599 2,7100 1,000 0,9819 48 0,0481 2,7899 2,7399 1,0109 0,9927 64 0,0417 2,7859 2,7478 1,0094 0,9956 На основании рис. 6 и данных в табл. 2 можно сделать вывод, что метод конструктивной анизотропии для оболочки варианта 1 сходится при числе вырезов . В табл. 3 показаны аналогичные данные для оболочки варианта 2, а на рис. 7 - график зависимости относительной критической нагрузки от общего числа вырезов для способов распределения вырезов I и II. Рис. 6. График зависимости относительной критической нагрузки потери устойчивости от числа вырезов для оболочки варианта 1 при Рис. 7. График зависимости относительной критической нагрузки потери устойчивости от числа вырезов для оболочки варианта 2 при Таблица 3 Результаты для оболочки варианта 2 при Критические нагрузки потери устойчивости , МПа , МПа , МПа 4 0,2000 0,7138 0,6438 0,6288 0,9019 0,8809 6 0,1633 0,5688 0,5388 0,7968 0,7548 8 0,1414 0,5813 0,5463 0,8144 0,7653 12 0,1155 0,6075 0,5538 0,8511 0,7758 16 0,1000 0,6263 0,6925 0,8774 0,9702 20 0,0894 0,6550 0,6700 0,9177 0,9387 24 0,0816 0,6475 0,6325 0,9072 0,8862 30 0,0730 0,6763 0,6675 0,9475 0,9352 36 0,0667 0,7038 0,6850 0,9860 0,9597 48 0,0577 0,7088 0,6925 0,9930 0,9702 64 0,0500 0,7100 0,6957 0,9947 0,9748 Как видно из графика (см. рис. 7) и табл. 3, метод конструктивной анизотропии для оболочки варианта 2 сходится при Пользуясь полученными данными, предположим, что для оболочки варианта 3 метод конструктивной анизотропии сойдется уже при Рис. 8. Кривая «нагрузка - наибольший прогиб» для варианта оболочки 3 при способе введения вырезов II На рис. 8 показан график «нагрузка - наибольший прогиб» для оболочки варианта 3 при способе задания вырезов II. При этом то есть метод конструктивной анизотропии здесь также сходится. Также интерес могут представлять поля прогибов и напряжений рассматриваемых конструкций. На рис. 9, 10 показаны состояния конструкции для оболочки варианта 3 при до и после местной потери устойчивости. Цветом показана интенсивность напряжений Часто при отображении оболочки с наложенным полем прогибов, особенно до потери устойчивости, значения прогибов настолько малы, что их достаточно трудно заметить. Тем не менее даже малозаметные изменения в конструкции могут нести в себе важную информацию о ее состоянии. В связи с этим значения прогиба умножаются на масштабирующий коэффициент Рис. 9. Вариант оболочки 3 до потери устойчивости Рис. 10. Вариант оболочки 3 после потери устойчивости Из данных на рис. 10 видно, что интенсивность напряжений превышает значения предела текучести МПа. Таким образом, потеря прочности происходит во время местной потери устойчивости. Заключение В ходе исследования были рассмотрены две математические модели деформирования пологих оболочек, ослабленных вырезами. В первой модели вырезы вводятся дискретно с помощью единичных столбчатых функций, задающих места расположения вырезов. Во второй модели используется метод конструктивной анизотропии, позволяющий «размазать» нулевую жесткость вырезов по полю оболочки. Использование математической модели оболочки ступенчато-переменной толщины, когда область интегрирования является односвязной, не вызывает в дальнейшем существенных трудностей при проведении расчетов. Были исследованы три варианта оболочек при разном числе вырезов. В результате было продемонстрировано, что при увеличении числа вырезов дискретность их ввода теряется и становится возможным использование метода конструктивной анизотропии. При способе задания вырезов I потеря устойчивости наступает позднее, чем при способе задания вырезов II. Критерием использования метода конструктивной анизотропии для оболочек, ослабленных вырезами, являются два фактора: отношение общей площади вырезов к площади оболочки ( ) и отношение ширины выреза к линейному размеру оболочки ( ). При решение, полученное по методу конструктивной анизотропии, практически совпадает с решением, полученных при дискретном введении вырезов при . При должно выполняться условие . Таким образом, чем меньше тем больше может быть , но не превосходить 0,09.

About the authors

I V Kamenev

Saint Petersburg State University of Architecture and Civil Engineering

A A Semenov

Saint Petersburg State University of Architecture and Civil Engineering

References

Statistics

Views

Abstract - 200

PDF (Russian) - 60