CC BY
32
УДК 519.6 DOI: 10.14529/mmpl80112
ОСОБЕННОСТИ ЗАДАНИЯ ГРАНИЧНЫХ УСЛОВИЙ ПРИ МОДЕЛИРОВАНИИ СЕТЧАТЫХ АНИЗОГРИДНЫХ КОНСТРУКЦИЙ
Т.В. Бурнышева, O.A. Штейнбрехер, А.Д. Ульянов
Новокузнецкий институт (филиал) Кемеровского государственного университета, г. Новокузнецк, Российская Федерация
В работе рассматривается влияние способа задания граничных условий при численном расчете напряженно-деформированного состояния сетчатых анизогридных конструкций. Численный расчет основан на применении метода конечных элементов в вариационной постановке. В работе приведены варианты задания граничных условий на нижней и верхней кромках конструкции. Приведены схемы закрепления в виде шарнирной заделки и жесткого закрепления для нижней кромки конструкции. Рассмотрены варианты задания осевого нагружения в виде распределенной нагрузки и «жесткого» узла дая верхней кромки конструкции. Приведена кинематика «жесткого» узла для конструкций с обшивкой и без обшивки. В качестве иллюстрации влияния вида закрепления на вычисления приведено сравнение численных расчетов и натурного эксперимента. Результаты работы показывают, что способ задания граничных условий оказывает значительное влияние на результаты расчета.
Ключевые слова: граничные условия; метод конечных элементов; дискретное моделирование; сетчатые оболочки.
Введение. Многие конструкции и сооружения, применяемые в строительстве и машиностроении, имеют вид многоэлементных сетчатых конструкций, образованных регулярной реберной структурой (рис. 1). В аэрокосмической отрасли, как правило, используют сетчатые оболочки, выполненные из композиционных материалов, которые являются силовыми и выступают составными элементами агрегатов космических летательных аппаратов [1, 2]. В этом случае набор структурных элементов конструкции должен обеспечивать прочность и жесткость оболочки под действующими нагрузками. В зависимости от наличия особенностей, таких как усиления и вырезы, различают анизогридные композиционные оболочки регулярной и нерегулярной структуры.
1. Дискретное моделирование. В работе [3] было показано, что для моделирования сетчатых конструкций, имеющих особенности, лучше использовать методы, основанные на дискретном моделировании, такие как метод конечных элементов (МКЭ) [4].
Одной из наиболее частых задач, связанных с данным классом конструкций, является расчет напряженно-деформированного состояния конструкции при осевом сжатии. Для этой задачи можно записать постановку в вариационном виде как минимум функционала
П = min П(и) (1)
Рис. 1. Сетчатые анизогридньте конструкции при граничных условиях в перемещениях
и\Г2 =0 (2)
и внешних нагрузках
P |А = Р, (3)
где и = \и, v, WT ,u £ U - варьируемое толе перемещений, U - пространство
Р
приложенная к верхней кромке оболочки), Г поверхность оболочки,Г ее верхняя кромка, Г2 - ее нижняя кромка, П(и) - потенциальная энергия как функционал, зависящий от перемещений.
Для большей согласованности результатов численного моделирования с данными натурного эксперимента требуется адекватное задание граничных условий задачи (1) (3). Известные CAD-системы, как правило, не предоставляют конечным пользователям (инженерам и конструкторам) возможности выбора способа задания граничных условий. Кроме того, алгоритмы и модели, заложенные в CAD-системе, закрыты для конечного пользователя. Возможность вариации задания граничных условий могут предоставить CAD-системы с открытым кодом. Такой системой является пакет программ «Композит-НК Анизогрид> [5], разработанный на базе Новокузнецкого института (филиала) Кемеровского государственного университета. Данный пакет программ имеет открытый код, что позволяет инженерам и конструкторам вносить изменения в алгоритм расчета. Данный пакет апробирован на расчете напряженно-деформированного состояния и устойчивости сетчатых анизогридньтх конструкциях. При этом сопоставления проводились как в сравнении с аналитическими расчетами, так и с известными CAD-системами. В настоящее время пакет программ используется при моделировании и проектировании сетчатых анизогридньтх конструкций АО «Центральный научно-исследовательский институт специального машиностроения> (г. Хотьково, Московская область). Дальнейшие рассуждения будут приведены в рамках алгоритмов, реализованных в пакете программ.
Как правило, CAD-системы моделирования основаны на использовании МКЭ. Тогда ребра анизогридной сетчатой оболочки можно представить в виде набора балок,
а обтпивка в виде набора пл&стин. В качестве конечного элемента для ребер можно использовать двухузловую балку Тимошенко (рис. 2 б), а в качестве элемента для обтттивки трехузловой треугольник Зенкевича (рис. 2 а).
Рис. 2. Конечные элементы: для пластины (а), для балки (б), где (в,1,и) местная система координат на элементе
Узлы конечно-элементой сетки верхней и нижней кромок могут принадлежать, как только одномерным конечным элементам балок, так и одновременно элементам балок и двумерным элементам обтттивки. в случае ее наличия.
2. Задание граничных условий. Для данного класса конструкций для задания граничного условия (2) на нижней кромке Г2 следует проводить закрепление в узлах перемещений и углов поворота либо шарнирным закреплением, либо жестким закреплением всех линейных перемещений и углов поворота узлов (заделкой). Шарнирное закрепление должно учитывать расположение точки крепления по толщине ребра (например, кольцевого) в физической модели и соотносить его с местом расположения в дискретной модели; в эту точку должен быть помещен узел конечно-элементной модели.
На рис. 3 представлена схема шарнирного закрепления перемещений щ в точке со сдвигом по толщине спирального ребра. Точка О расположена на нижней кромке оболочки и сдвинута по толщине кольцевого ребра на величину еп от узла О на поверхности приведения. Закрепление перемещений вдоль оси Ь в точке О будет соответствовать следующему условию щ + еп • 93 = 0, где еп - эксцентриситет для точки О. Если же выбрать узел в точке О, то граничное условие принимает канонический вид: щ = 0. Положение сечения относительно узла корректируется заданием эксцентриситета и введением дополнительного преобразования перемещений.
Рис. 3. Схема шарнирного закрепления в точке, отстоящей от поверхности приведения на нижней кромке оболочки: крупная штриховка кольцевое ребро, мелкая штриховка обтттивка
__
О'". е„ 'О
п
При решении задач на осевое сжатие, действующая сила (3), как правило, прикладывается к верхней кромке конструкции Г\. Используют два способа задания нагрузки: распределение сжимающей силы, в том числе и с учетом приложенного момента, по узлам верхней кромки и задание нагрузки через «жесткий> узел. Акту-эльным является исследование соответствия способа приложения нагрузки и способа ее моделирования для сетчатых анизогридньтх конструкций.
Рассмотрим кинематику «жесткого узла, расположенного в центре окружности, определяющей верхнюю кромку сетчатой оболочки (рис. 36).
«Жесткий> узел, со всеми присоединенными к нему элементами, может перемещаться только как жесткое тело, это позволяет привести нагрузку, распределенную по кромке произвольным образом, к действию сосредоточенных сил и моментов, приложенных в «жестком> узле.
«Жесткий> узел имеет шесть степеней свободы три поступательных и три вращательных, причем эти перемещения отсчитьтваются в главной системе координат этого узла. Удобно выбирать оси этой системы параллельными осям общей декартовой системы координат (У^,Х) (рис. 4).
Рис. 4. Схема «жесткого узла Степени свободы «жесткого узла определим как:
5
[пу, иг, их■ Оу, Ог, вх]Т ■ (4)
Перемещения всех точек элементов, присоединенных к жесткому узлу, выражаются через перемещения жесткого узла и эксцентриситеты (координаты точек в главной системе жесткого узла) следующим образом:
(5)
иV = их ■ сов(Х, в') + иу ■ сов(У, в') + иг ■ соа^, в') + Ох[—Zv ■ сов(У, в') +
+Уи ■ сов^, в')] + ОуZv ■ сов(Х, в') — ОгУ» ■ сов(Х, в'), = их ■ сов(Х, Ь') + иу ■ сов(У, Ь') + иг ■ сов^, Ь') + Ох [—Zv ■ сов(У, Ь') + +УV ■ со8^,г')] + Оу Zv ■ сов(Х,Ь') — О г Уи ■ сов(Х,г'), иП = их ■ сов(Х, п') + иу ■ сов(У, п') + иг ■ cos(Z, п') + Ох[—Zv ■ сов(У, п') +
+У, ■ cos(Z, п')] + ОуZv ■ сов(Х, п') — ОгУи ■ сов(Х, п'), эксцентриситеты узла элемента относительно «жесткого узла V, {в', Ь', п'} - главная система координат. Степени свободы узла V, расположенного на верхней кромке оболочки Г2 и входящего в одномерный балочный элемент спирального или кольцевого ребра, будут иметь вид:
где Yu, Zu, X,
Us
Ut'
Un
Os', Ot', 9n']T = Lv [uy, Uz, UX. 0Y, Oz, в
]T
X\
(G)
где - матрица перехода от главной системы «жесткого узла (У, Z, X) к главной системе [в',п'} узла. Степени свободы узла, расположенного на верхней кромке оболочки и входящего в треугольный конечный элемент обшивки, будут иметь вид:
П, щ>, пп!. О3>, 0г>] = [пу, uz, их. Оу, 9z, Ох] , (7)
где матрица - матрица перехода от главной системы «жесткого узла (У, Z, X) к главной системе [в',1', П} узла . Распределенная нагрузка задается или как равномерное распределение, или как доля осевой нагрузки с учетом приложенного момента. Таким образом, элементы верней кромки моделируются отдельно.
3. Результаты моделирования сетчатой анизогридной конструкции. Моделирование приложенной нагрузки через «жесткий> узел было использовано при численном расчете составной силовой конструкции (рис. 5), включающей стрингерную алюминиевую часть и сетчатую композитную часть с вырезами и усилениями [6, 7].
Рис. 5. Составная цилиндрическая оболочка: общий вид (а), схема расположения датчиков (б) и распределение осевых перемещений в композиционной части конструкции (в): 1 данные натурного эксперимента. 2 и 3 результаты численного расчета при задании нагрузки с использованием «жесткого узла и задании распределенной нагрузки
«
НК Анизогрид>. Натурный эксперимент на осевое сжатие проводился в АО «Центральный Научно-исследовательский институт специального машиностроения > пу-
тем нагружения сетчатой анизогридной оболочки через жесткую плиту. Результаты натурного эксперимента и численного расчета представлены в таблице.
Таблица
Результаты натурного эксперимента и численного расчета
Датчик Измеренные перемещения, мм Рассчитанные перемещения, мм
Нагружение через «жесткий> узел Распределенное нагружение
ВП1 0,86 0,963 0,6382
ВП2 1,363 1,3038 1,3091
ВПЗ 1,083 1,161 0,8045
ВП4 1,363 1,4071 1,4631
ВП5 1,155 0,8905 0,719
ВП6 1,053 1,1713 1,1756
ВП7 0,829 0,8884 0,6021
ВП8 1,195 1,0687 1,0225
Разница между расчетами с заданием нагрузки распределено и через «жесткий> узел в точках крепления датчиков составила 2 - 32%. Сравнение результата численных расчетов с натурным экспериментом (рис. 5 в) показало, что при задании нагрузки через «жесткий> узел погрешность расчета составила 3 - 12%, а при распределенной нагрузке - 4 - 37%.
Таким образом, моделирование через «жесткий> узел в случае с нерегулярной сеткой, нагруженной через плиту, имеет большую согласованность с натурным экспериментом, в отличие от распределенного задания нагрузки.
Вывод. Схемы задания граничных условий при численном моделировании напряженно-деформированного состояния силовых анизогридных конструкций оказывают значительное влияние на результаты расчетов. Результаты моделирования позволяют предположить, что задание нагрузки через «жесткий> узел хорошо согласуется с данными натурных экспериментов, в которых нагружение осуществляется посредством плиты. Тем не менее для каждой конструкции требуется проведение дополнительных исследований о корректности применения схем задания граничных условий.
Литература
1. Васильев, В.В. Композитные сетчатые балочные элементы для авиационных конструкций / В.В. Васильев, А.Ф. Разин // Вопросы оборонной техники. - 2006. - № 1/2. -С. 3-6.
2. Васильев, В.В. Апизогридные композитные сетчатые конструкции - разработка и приложение к космической технике / В.В. Васильев, В.А. Барынин, А.Ф. Разин, С.А. Петровский, В.И. Халиманович // Композиты и наноструктуры. - 2009. - № 3. - С. 38-50.
3. Бурнышева, Т.В. Сравнение дискретного и континульного подхода к расчету напряженного состояния сетчатых оболочечных конструкций при статическом нагружении / Т.В. Бурнышева, В.О. Каледин // Научно-технический вестник Поволжья. - 2011. -№ 4. - С. 113-116._
| .|2 Bulletin of the South Ural State University. Ser. Mathematical Modelling, Programming
& Computer Software (Bulletin SUSU MMCS), 2018, vol. 11, no. 1, pp. 137-144
4. Бате, К.-Ю. Методы конечных элементов / К.-Ю. Бате. - М.: Физматлит, 2010.
5. Каледин, В.О. Программная система для алгоритмизации численного решения задач механики сплошных сред / В.О. Каледин, Я.С. Крюкова, Н.В. Нагайцева, Е.В. Рав-ковская // Известия Алтайского государственного университета. - 2014. - № 1-1 (81). -С. 161-164.
6. Бурнышева, Т.В. Оценка взаимного влияния составных частей оболочечной конструкции при осевом сжатии / Т.В. Бурнышева, А.Ф. Разин, O.A. Штейнбрехер // Научно-технический вестник Поволжья. - 2012. - № 2. - С. 127-131.
7. Каледин, В.О. Интерпретация данных натурных испытаний оболочечной композитной конструкции при статическом осевом сжатии / В.О. Каледин, А.Ф. Разин, Т.В. Бурнышева, O.A. Штейнбрехер // Заводская лаборатория. Диагностика материалов. - 2015. -Т. 81, № 3. - С. 53-58.
Татьяна Витальевна Бурнышева, доктор технических наук, доцент, декан, факультет «Информационные технологии:», Новокузнецкий институт (филиал) Кемеровского государственного университета (г. Новокузнецк, Российская Федерация), [email protected].
Ольга Александровна Штейнбрехер, старший преподаватель, кафедра «Инфор-
:
(филиал) Кемеровского государственного университета (г. Новокузнецк, Российская Федерация), [email protected].
«
:
университета (г. Новокузнецк, Российская Федерация), [email protected].
Поступила в редакцию 18 сентября 2017 г.
MSC 74S05 DOI: 10.14529/mmpl80112
FEATURES OF THE BOUNDARY CONDITIONS IN THE MODELLING MESH DESIGNS ANIZOGRIDNYH
T.V. Burnysheva, O.A. Shteinbreher, A.D. Ulyanov
Novokuznetsk Branch Institute of Kemerovo State University, Novokuznetsk,
Russian Federation
E-mail: [email protected], [email protected], [email protected]
The influence of boundary conditions in the numerical calculation of the stress-strain construction of anisogrid structures is considered in the paper. The numerical calculation is based on the application of the finite element method in the variational formulation. The paper presents variants of the fulfillment of boundary conditions on the lower and upper edges of the construction. The schemes of fastening in the form of hinged closing and rigid fastening for the bottom edge of designs. The variants of specifying axial loading in the form of a distributed load and a "hard" node for the upper edge of the structures are considered. The article presents the kinematics of a "hard" node for lattice construction. As an example of the effect of fixing on calculations, a comparison of the number calculations and field experiment. The results show that the way to set the boundary conditions has a significant effect on the calculation results.
Keywords: boundary conditions; finite element method; discrete modelling; lattice shells.
T.B. Eypiibiiiieisa, O.A. HlTeiiHÖpexep, A./J. YjibHiiOB
References
1. Vasil'ev V.V. Mekhanika konstruktsiy iz kompozitnykh materialov [Mechanics of Constructions from Composite Materials]. Moscow, Mechanical engineering, 1988. (in Russian)
2. Vasil'ev V.V., Barynin B.A., Razin A.F., Petrakovskiy S.A., Khalimanovich V.l. [Anizogrid Composite Mesh Design - Development and Application of Space Technology]. Composites and Nanostructures, 2009, no. 3, pp. 38-50. (in Russian)
3. Burnysheva T.V., Kaledin V.O. [Comparison of Discrete and Continual Approach to Calculation of a Stressed State of Network Capsulate Designs at a Static Loading]. Nauchno-tekhnicheskiy vestnik Povolzh'ya, 2011, no. 4, pp. 113-115. (in Russian)
4. Bathe K.-Yu. Finite Elements Procedures. New Jersey, Prentice Hall, Pearson Education, 2010.
5. Kaledin V.O., Kryukova Ya.S., Nagaytseva N.V., Ravkovskaya E.V. [Software System for Algorithmization Numerical Solution of Problems of Continuum Mechanics]. Izvestiya Altayskogo gosudarstvennogo universiteta, 2014, no. 1-1 (81), pp. 161-164. (in Russian)
6. Burnysheva T.V., Razin A.F., Shteynbrekher O.A. [Evaluation of the Mutual Influence of the Parts Shell Structures under Axial Compression]. Nauchno-tekhnicheskiy vestnik Povolzh'ya, 2012, no. 2, pp. 127-131. (in Russian)
7. Kaledin V.O., Razin A.F., Burnysheva T.V., Shteynbrekher O.A. [Interpretation of These Field Tests of a Capsulate Composite Design at a Static Axial Compression]. Factory Laboratory. Diagnostics of Materials, 2015, vol. 81, no. 3, pp. 53-58.(in Russian)
Received September 18, 2017
