Full Text

Введение Одной из задач теории упругости, решение которой имеет прикладное значение, является задача об изгибе шарнирно закрепленной многослойной панели с круговой осью, нагруженной равномерно распределённой нагрузкой. Использование в данном случае конечно-элементного метода при изменении геометрии слоёв и их упругих констант требует перестроения модели, что эквивалентно каждый раз решению новой задачи. Использование упрощённых моделей предусматривает введение гипотез относительно перемещений и упругих характеристик материалов слоёв, что сужает область использования таких моделей. Поэтому определённый интерес представляет получение аналитического решения, которое описывает изгиб шарнирно закрепленной многослойной панели с круговой осью без введения упрощающих гипотез. Решению статической задачи об изгибе многослойных пластин и оболочек при нагружении распределённой нагрузкой посвящено большое количество работ [1-27]. В исследовании [1] решена задача о цилиндрическом изгибе многослойной пластины при равномерно распределённой нагрузке. В статье [2] результаты работы [1] использованы для решения задачи об изгибе пластины в 3D-постановке. В работе [3] рассматривается цилиндрический изгиб пластины под действием локальной нагрузки. Автором работы [4] предложено точное решение 3D-задачи об изгибе пластины, где в качестве неизвестных функций использованы перемещения и поперечные напряжения. Использование смешанной формы закона Гука позволило краевые условия задачи сформировать в напряжениях и понизить порядок системы разрешающих уравнений до первого. Система разрешающих уравнений решена методом Навье. В работах [5; 6] рассмотрены пластины, упругие характеристики которых изменяются по толщине. В исследованиях [7; 8] авторы предложили решение для конструкций с заполнителем. Работы [9-13] посвящены решению плоской задачи теории упругости для прямолинейной многослойной балки. Аппроксимирующие функции в этих работах построены на основе полиномов Лежандра. В работе [14] предложен метод решения уравнений равновесия для пластин с поперечной анизотропией и модулем упругости, линейно изменяющимся по толщине пластины. В статьях [15-16] предложено замкнутое решение трёхмерной задачи о деформации шарнирно закреплённой прямоугольной пластины при действии механических и термических нагрузок для свободных и вынужденных колебаний шарнирно закреплённой пластины с поперечной анизотропией. В работе [17] предложено точное решение для прямоугольных в плане анизотропных композитных ламинатов с шарнирно закреплёнными кромками с упругими свойствами, переменными по толщине слоёв Значительно меньше работ посвящено приложениям теории упругости к расчёту оболочек, что связано со сложностью формирования и решения системы разрешающих уравнений. В работе [18] предложена процедура вывода уравнений деформирования многослойной оболочки с использованием в качестве аппроксимирующих функций полиномов Лежандра. В работе [19] приведено решение трёхмерной задачи теории упругости об изгибе тонкой многослойной оболочки под действием синусоидальной нагрузки. В работе [20] рассматривается трёхмерная модель для исследования свободных колебаний цилиндрической оболочки. Авторы статьи [21] для статического, динамического расчёта композитных цилиндрических панелей использовали систему дифференциальных уравнений, записанную в цилиндрической системе координат. В исследовании [22] разрешающие уравнения записаны в ортогональных криволинейных координатах для сферической оболочки с использованием теоремы Гамильтона - Кэлли (Cayley - Hamilton theorem). В качестве неизвестных использованы три перемещения и три компоненты поперечных напряжений. Наиболее близким к предлагаемому методу является подход, изложенный в работах [23-27]. Так, в работе [23] предложен метод приближенного решения задач изгиба многослойных пластин, цилиндрических и сферических оболочек при шарнирном закреплении кромок, согласно которому рассматриваемая область делится на слои. В пределах каждого слоя записывается система дифференциальных уравнений равновесия при постоянном радиусе кривизны. Слои связываются друг с другом условиями непрерывности перемещений и поперечных напряжений (условиями сопряжения) при переходе через границу раздела слоёв. Однако для получения необходимой точности в данном случае необходимо разбивать область на большое количество слоёв (порядка 100). Недостатком данного метода является его большая вычислительная нагрузка, поскольку для каждого слоя приходится записывать дифференциальные уравнения. Кроме того, все перечисленные методы связаны так или иначе с трудоёмкой, особенно для криволинейных оболочек, процедурой решения системы дифференциальных уравнений. Анализ вышеперечисленных работ показал, что решение задачи изгиба панели с круговой осью может быть выполнено либо приближенно или в численном виде, потому что: - во-первых, система уравнений равновесия Коши для этой задачи является нелинейной; - во-вторых, традиционный метод решения через обобщённую функцию напряжений Airy’s [28] при анализе многослойных панелей, слои которых имеют значительную поперечную анизотропию упругих характеристик при переходе через границу раздела слоёв, не пригоден в связи со сложностью формирования граничных условий смешанного типа и условий сопряжения слоёв. В связи с этим для решения задачи изгиба, как правило, используются приближённые методы. Анализ литературы показал, что в случае шарнирного закрепления кромок решение задачи значительно упрощается, а именно двухмерную задачу можно свести к одномерной, которая описывается системой нелинейных дифференциальных уравнений относительно одной переменной, измеряемой по толщине панели. В работе [29] предложен метод решения задачи об изгибе многослойной панели с прямолинейной осью для условий шарнирного закрепления её кромок. В данной работе этот метод применён к панели с круговой осью. При этом выражения, описывающие изменение окружных , радиальных перемещений и приложенной к ней нагрузки , записываются в виде функций: , ,, где функции и представляют собой полиномы степени ; k - номер слоя: , . (1) Коэффициенты полиномов (1) определяются из условий: - равновесия на наружном и внутреннем контуре модели; - равновесия на боковых гранях модели; - сопряжения слоёв; - удовлетворения системе дифференциальных уравнений равновесия Коши в нескольких точках по толщине каждого слоя. При этом решение задачи сводится к решению системы линейных уравнений относительно неизвестных коэффициентов и полиномов (1). От прямых вариационных методов, которые также сводятся к решению системы линейных уравнений, предлагаемый метод выгодно отличается тем, что: - во-первых, не требует выполнения сложных операций интегрирования при вычислении соответствующих функционалов; - во-вторых, система уравнений равновесия Коши в предлагаемом методе удовлетворяется точно, хотя бы в нескольких точках каждого слоя, в то время как в прямых вариационных методах эта система выполняется лишь приближённо. В работе [29] показано, что минимальная степень полиномов в предлагаемом методе должна быть не менее двух. При этом обеспечивается выполнение условий на контуре модели, условий сопряжения слоёв и системы дифференциальных уравнений в одной точке в середине толщины слоя. Дальнейшее увеличение степени полиномов возможно посредством увеличения количества точек, для которых выполняется система уравнений равновесия Коши. В случае использования полиномов первой степени задача имеет тривиальное решение: так, все неизвестные коэффициенты однозначно определяются граничными условиями и условиями сопряжения слоёв. При этом система дифференциальных уравнений равновесия Коши остаётся невыполненной. Поэтому в данной работе ставится задача решения плоской задачи изгиба многослойной панели с круговой осью при шарнирном закреплении его концевых сечений с использованием подхода [29]. 1. Однослойная панель Вначале рассмотрим в полярной системе координат однослойную панель с круговой осью. Радиус наружной поверхности панели , внутренней . Толщина панели в радиальном направлении равна Н, а его угловая длина . Концевые сечения имеют шарнирное закрепление. К наружной поверхности с радиусом приложена радиальная нагрузка . Необходимо отметить, что в случае нагрузки произвольного типа последняя раскладывается в ряд Фурье. Описанная методика применяется для каждой гармоники полученного ряда. Материал панели ортотропный с механическими характеристиками: - модули упругости в направлении координатных осей , ; - модуль поперечного сдвига в плоскости ; , - коэффициенты Пуассона. Упругие константы материала определялись по формулам: , , . Коэффициенты , полиномов (1) находились из условий: 1. Условие равновесия на наружном, внутреннем контуре при , , ; при 2. Условие равновесия на боковых гранях: . (2) 3. Система дифференциальных уравнений равновесия Коши в полярной системе координат: (3) 4. Условие равновесия слоёв при , где : . 5. Условие совместности деформации слоёв при , где : , где m - количество слоёв; k - номер слоя, - радиальная координата границы раздела слоёв; - радиальные нормальные, окружные нормальные и касательные напряжения. В формуле (2) выражение в правой части есть реакции на боковых гранях панели, которые найдены из условий его равновесия под действием нагрузки . Деформации связаны с окружными и радиальными перемещениями через соотношения Коши, записанные в полярной системе координат, а напряжения связаны с деформациями с помощью закона Гука [28]. Описанные выше условия 1, 2, 4 и 5 в предлагаемом методе удовлетворяются точно, в то время как система дифференциальных уравнений (3) выполняется точно только в нескольких точках по толщине каждого слоя. Степень нелинейности перемещений и напряжений по толщине панели определяется отношениями и . Чем меньше значения этих отношений, тем выше степень нелинейности указанных перемещений. Поэтому исследование влияния степени полинома на точность решения проводилось для различных значенийи . В качестве эталонной использовалась конечно-элементная модель, построенная в системе ANSYS, которая состояла из квадратичных элементов plane 183. Качество сетки оценивалось как с помощью стандартных процедур, которые используются в системе ANSYS, так и по условиям равновесия на контуре конечно-элементной модели и на границе раздела слоев [16]. В расчете использовались следующие параметры. Упругие характеристики: кгс/мм2, Геометрические размеры модели: мм, . Для анализа рассматривались модели с относительной толщиной Степень полиномов изменялась от 2 до 8. Интенсивность нагрузки, изменяющейся по закону синуса, равнялась 20 кгс/мм. Параметры напряженно-деформированного состояния однослойной панели, такие как радиальные, окружные и касательные напряжения и перемещения по толщине панели, сравнивались с аналогичными параметрами, найденными по эталонной конечно-элементной модели. На рис. 1 показана эпюра радиальных напряжений в кгс/мм2 по толщине однослойной панели. Точками показаны результаты, полученные с помощью эталонной конечно-элементной модели, а сплошной линией - с помощью предлагаемого метода. а b Рис. 1. Эпюры радиальных напряжений по толщине панели: а -; ; b - Fig. 1. Radial stress diagrams vs plate thickness: а - ; ; b - а b Рис. 2. Зависимость точности расчёта перемещений (а) и напряжений (b) от степени полиномов Fig. 2. Effects of polynomial degree on the calculation accuracy for displacement (a) and stress (b) Выбор степени полинома осуществлялся его повышением до достижения заданной точности. На рис. 2 представлены соответственно зависимости точности расчёта перемещений и напряжений от степени полинома (точки - результаты расчёта, линии - результаты по методу наименьших квадратов). Горизонтальная линия соответствует точности, принимаемой в инженерных расчётах, - 5 %. Например, при для расчёта перемещений с точностью 5 % достаточно степени полинома 3, а для расчёта напряжений с той же точностью необходима степень полинома 5. Таким образом, расчетные значения, полученные по конечно-элементной модели и по предлагаемому методу хорошо совпадают в широком диапазоне нелинейности от 2 до 20. Полученные результаты можно использовать в дальнейшем в качестве начального приближения при выборе степени полинома для слоёв многослойной панели. 2. Многослойные панели Предлагаемый метод расчета изгиба многослойной панели рассматривается на примере трехслойной панели, поскольку с увеличением количества слоев он остается неизменным. В этом случае методика остается неизменной, необходимо лишь добавить условия равновесия и условия совместности деформаций слоёв. В качестве примера рассмотрен расчет напряженно-деформированного состояния при изгибе трехслойной панели, конфигурация и основные геометрические размеры которой приведены на рис. 3. Радиусы наружной и внутренней поверхностей равны соответственно и Радиусы границ слоёв 1-2 и 2-3 соответственно равны и . Угловая длина панели равна . Упругие характеристики несущих слоёв и слоя заполнителя представлены в таблице. Рис. 3. Геометрия трёхслойной панели с круговой осью Fig. 3. Three layers circular panel geometry Степень нелинейности перемещений увеличивается с увеличением толщины слоя заполнителя. Поэтому рассматривались панели с несущими обшивками одинаковой толщины, равной 10 мм, но разной толщиной слоя заполнителя: 20, 40, 60 мм. При этом мм, α = π/3. На рис. 4 представлены эпюры окружных перемещений по толщине панели, которые были получены с помощью эталонной конечно-элементной модели (точки) и с помощью предлагаемого метода (сплошные линии) для толщин заполнителя 20 мм (а) и 60 мм (b). С ростом толщины заполнителя увеличивается степень нелинейности окружных перемещений заполнителя. На рис. 5-7 представлены эпюры касательных, нормальных радиальных и окружных напряжений, полученных с помощью эталонной конечно-элементной модели (точки) и с помощью предлагаемого метода (сплошные линии) для толщин заполнителя 20 мм (а) и 60 мм (b). Характеристики материала обшивок и заполнителя Skin and core material properties Характеристика Несущие слои Заполнитель Eq, кгс/мм2 2·104 0,01 Er, кгс/мм2 2·104 200 Gq, кгс/мм2 7690 2 mqr, mrq 0,3 0,3 а b Рис. 4. Эпюры окружных перемещений по толщине панели Fig. 4. Through thickness circular stress distribution а b Рис. 5. Эпюры касательных напряжений по толщине панели Fig. 5. Through thickness shear stress distribution а b Рис. 6. Эпюры нормальных радиальных напряжений по толщине панели Fig. 6. Through thickness radial stress distribution а b Рис. 7. Эпюры окружных напряжений по толщине панели Fig. 7. Through thickness circular stress distribution Рис. 8. Влияние на Δ для перемещений и напряжений Fig. 8. Effects on the Δ for displacement and stress Анализ эпюр показывает хорошее совпадение результатов, полученных с помощью эталонной конечно-элементной модели и с помощью предлагаемого метода. Погрешность не превышала 0,5 %. Дополнительно проведено сравнение результатов, полученных с помощью предлагаемого метода, и результатов, полученных с помощью модели, построенной на основе вариационного исчисления при линейном распределении продольных перемещений по толщине заполнителя [31]. На рис. 8 приведена зависимость разницы Δ между результатами вышеперечисленных расчётов перемещений и напряжений в зависимости от отношения Расчёты выполнены для трёхслойной панели с толщинами обшивок 2 мм. Толщина заполнителя равна: 15, 25, 35, 43, 55, 65 мм. Из рис. 8 следует, что при расчёте перемещений разница между результатами не превышает 5 %, если , а при расчёте напряжений - если . То есть для тонких панелей с можно использовать гипотезу о линейном распределение перемещений по толщине заполнителя. Заключение 1. Предлагаемый метод позволяет оценивать напряженно-деформируемое состояние при изгибе шарнирно закрепленных многослойных панелей с круговой осью с произвольным количеством слоев без введения упрощающих гипотез относительно перемещения и механических характеристик слоев. Точность метода сравнима с точностью метода конечных элементов при использовании квадратичных твёрдотельных элементов. 2. Новизна предлагаемого метода заключается в том, что полученное общее аналитическое решение справедливо для широкого диапазона изменения механических характеристик слоев, их геометрии и не связано с трудоёмкой процедурой решения системы разрешающих дифференциальных уравнений. 3. Целесообразность использования данного метода при проектировании обусловлена простотой его применения, высокой точностью и малыми временными затратами на анализ при переборе большого количества вариантов конструкции панели.

About the authors

S. A Bukaty

PJSC UEC-Saturn

N. V Osadchy

PJSC UEC-Saturn

References

Statistics

Views

Abstract - 287

PDF (Russian) - 158