CC BY
13
МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ MACHINE BUILDING AND MACHINE SCIENCE
УДК 678.549
https://doi.org/10.23947/2687-1653-2022-22-2-107-115
H) Check for updates
Научная статья
Использование метода конечных элементов для моделирования резервуара высокого давления из полимера, армированного углеродным волокном
И. Р. Антибас : , А. Г. Дьяченко
Донской государственный технический университет, Российская Федерация, г. Ростов-на-Дону, пл. Гагарина, 1 И Imad.antvpas@mail.ru
Аннотация
Введение. За последнее десятилетие спрос в мире на резервуары высокого давления существенно возрос, особенно в таких областях промышленности как авиационная, космическая, химическая и нефтегазовая. Находясь под постоянным воздействием высокого внутреннего давления, стенки резервуаров испытывают повышенное напряжение, что может стать причиной их внезапного разрушения. Для устранения такой возможности и улучшения прочностных характеристик резервуары изготавливают в форме металлических цилиндров с внутренним покрытием из композитного материала, состоящего из смолы, армированной углеродными волокнами. Цель настоящей статьи заключалась в изучении влияния угла наклона углеродного волокна на цилиндрические резервуары и определении величины максимального разрушающего давления с использованием метода конечных элементов программы ANSYS.
Материалы и методы. С помощью программы ANSYS создана конечно-элементная модель резервуара, состоящего из центральной части в виде металлического цилиндра с внутренним покрытием из полимера, армированного углеродными волокнами. По торцам резервуара размещены полусферы со спиральной намоткой. Для моделирования слоистого композитного материала использована командная оболочка SHELL 99, для определения критерия разрушения резервуара высокого давления — теория Tsai-Wu. Результаты исследования. Модель цилиндрического резервуара рассчитывалась для двух видов намотки волокон: кольцевой и спиральной при различных углах их наклона. Результаты анализа величины давления для различных углов наклона волокон показывают, что, начиная со значения угла 0° и до 45о оно увеличивается, а затем до значения 65о — уменьшается. Наибольшее давление, которое может выдержать резервуар, армированный углеволокном, составляет 207 МПа при угле наклона волокон ± 45о.
Обсуждение и заключения. Анализ исследований показал, что при угле наклона волокон ± 45о максимальное напряжение оказалась наименьшим, а максимально возможное разрушающее давление при том же угле составило 207 МПа. Из этого следует, что оптимальный угол ориентации волокон для обеспечения безопасной работы резервуара высокого давления составляет ± 45о, а резервуар из углепластика при том же угле намотки волокон имеет максимальную прочность.
Ключевые слова: резервуар высокого давления, компьютерная модель, угол намотки, композиционное покрытие, углеродное волокно, полимерное связующее.
Благодарности: авторы выражают благодарность редакции и рецензентам за внимательное отношение к статье и указанные замечания, которые позволили повысить ее качество.
Для цитирования: Антибас, И. Р. Использование метода конечных элементов для моделирования резервуара высокого давления из полимера, армированного углеродным волокном / И. Р. Антибас, А. Г. Дьяченко // Advanced Engineering Research. — 2022. — Т. 22, № 2. — С. 107-115. https://doi.org/10.23947/2687-1653-2022-22-2-107-115
<и К X <и
<и ю о X
к 3
eö
<и К X <и О CP
н о о X X
3
CÖ
© Антибас И. Р., Дьяченко А. Г., 2022
Original article
Using the Finite Element Method to Simulate a Carbon Fiber Reinforced Polymer Pressure Vessel
Imad Rizakalla Antipas ED, Alexey G. Dyachenko
Don State Technical University, 1, Gagarin sq., Rostov-on-Don, Russian Federation И Imad.antypas@mail.ru
Abstract
Introduction. Over the past decade, global demand for pressure vessels has increased significantly, specifically in such industries as aviation, space, chemical, and oil and gas. Being under the constant impact of high internal pressure, the walls of the tanks are under increased stress, which can cause their sudden destruction. To eliminate this probability and improve the strength characteristics, the tanks are made in the form of metal cylinders with an internal coating of composite material consisting of resin reinforced with carbon fibers. This article aimed at studying the effect of the angle of inclination of carbon fiber on cylindrical tanks and determining the maximum destructive pressure using the finite element method of ANSYS program.
Materials and Methods. Using the ANSYS program, a finite element model of a tank was created. It has a central part, which is a metal cylinder with an internal coating of composite material consisting of polymer reinforced with carbon fibers. At the ends of the tank, spiral wound hemispheres were placed. In these studies, SHELL 99 was used to model the layered composite material. The Tsai-Wu theory was used to determine the pressure tank failure criterion. Results. The cylindrical tank model was calculated for two types of fiber winding paths: annular and spiral, at different angles of their inclination. The results of the pressure value analysis for different fiber inclination angles showed that, starting from the angle value of 0° and up to 45°, it increased, and then, up to the angle value of 65°, it began to decrease. The critical pressure value for a carbon fiber reinforced tank was 207 MPa, which was obtained at a fiber angle of 45o.
Discussion and Conclusion. Analysis of the studies showed that at a fiber inclination angle of 45o, the value of the maximum stress turned out to be the smallest, and the maximum possible destructive pressure at the same angle was 207 MPa. It follows, that the optimal fiber orientation angle to provide safe operation of the high-pressure tank is ± 45°, and the carbon fiber tank, calculated at the same fiber winding angle, has the maximum strength value. Keywords: high-pressure tank, computer model, winding angle, composite coating, carbon fiber, polymer binder. Funding information: the authors would like to thank the editorial board and the reviewers for their advertency to the article and the comments, which contributed to its quality.
For citation: I. R. Antypas, A. G. Dyachenko. Using the Finite Element Method to Simulate a Carbon Fiber Reinforced Polymer Pressure Vessel. Advanced Engineering Research, 2022, vol. 22, no. 2, pp. 107-115. https://doi.org/10.23947/2687-1653-2022-22-2-107-115
Введение. Использование резервуаров высокого давления из углеродных волокон нашло широкое распространение в различных областях промышленности из-за своих отличительных свойств, таких как малый вес и высокая прочность. Поэтому в последнее время спрос на такие резервуары существенно возрос в тех случаях, когда вес играет большую роль [1, 2].
Важнейшими областями применения резервуаров такого типа являются аэрокосмическая, авиационная и химическая техника. Кроме того, резервуары, армированные волокном, широко используются для увеличения напора при транспортировке нефти и газа. Во многих случаях резервуары такого назначения подвергаются ^ воздействию высокого внутреннего давления, следствием чего может стать скачкообразное повышение
и напряжения на стенках сосуда и внезапное их разрушение, что наносит большой ущерб материальным и д
о людским ресурсам [3, 4]. тз
В [5] проведен ряд исследований, направленных только на изучение напорного резервуара из многослой-£ ного композита, в ходе которых определена ожидаемая прочность на смятие. В [6] изучалось поведение враща-
СЛ
^ ющегося составного напорного резервуара под воздействием внутреннего давления и осевой нагрузки; в [7] проведены исследования по изучению влияния тепловых нагрузок, действующих на многослойный композитный напорный бак. В работе [8] исследовалось поведение многоугольного композитного напорного бака пяти различных форм под действием внутреннего давления различного режима, а в [9] — изучалась работоспособность составного напорного бака при воздействии поперечных нагрузок.
Проектирование резервуара из композитного материала является сложной проблемой, поэтому для проведе-108 ния полного и точного анализа необходимо выбрать приоритетные факторы. Авторами изучался резервуар высокого давления, армированный несколькими слоями углеродных волокон и подвергаемый нагрузке от внут-
о, £ л
реннего давления. Цель работы — изучение влияния угла наклона углеродного волокна на цилиндрические резервуары и определение величины максимального разрушающего давления с использованием метода конечных элементов программы ANSYS.
Материалы и методы. Теоретическое исследование. Резервуар высокого давления из композитных материалов (CFRP) состоит из центральной части, представляющей собой металлический цилиндр с внутренним покрытием, состоящим из полимера, армированного углеродными волокнами. Торцевые поверхности резервуара имеют форму полусфер со спиральной намоткой (рис. 1).
Рис. 1. Продольное сечение резервуара
Резервуары, работающие при высоких давлениях и армированные углеродным волокном, изготавливаются методом филаментной намотки (рис. 2). Для получения необходимой стабильности армирования волокна отправляются на движущуюся тележку при тщательном выборе их координации, а затем наматываются на цилиндрическую поверхность. На стабильность координации волокон влияет несколько факторов: температура, форма и обработка поверхности, а также степень сцепления волокон с матрицей. Угол наклона намотки контролируется скоростью перемещения тележки и скоростью вращения цилиндрического барабана. С целью получения высоких эксплуатационных и прочностных свойств внутренняя цилиндрическая поверхность резервуара покрывается несколькими слоями волокон [10].
Угол наклона волокна
Цилиндрический барабан
Движущаяся тележка
Полимер
— Волокна
#
!с Он
Рис. 2. Схема намотки волокон
Угол наклона волокон оказывает большое влияние на свойства резервуаров, поэтому важным вопросом является поиск соответствующего угла для каждой части резервуара. Угол ориентации волокон определяется необходимой величиной трения между волокнами и слоем композитного материала, как показано в соотношении:
а
(Я ) = 51
51П
(
±8
я - я
Л"
V я - я0 У
(1)
где я — расстояние между центром и точкой слоя; я — радиус центральной оси; я — радиус по касательной к поверхности цилиндра при 8 = 0 [8].
Критерии разрушения Tsai-Wu. При исследовании и моделировании резервуара высокого давления, выполненного из композитного материала, использовались критерии разрушения по теории Tsai-Wu [11]. Реализация уравнения (1) необходима для исследования ожидаемого разрушения ортотропных материалов по теории Tsai-Wu:
Р1°11 + Р2°22 + Р6°12 + ^1^2 + Р66ТП + 2^2СТ11СТ22 = 1 .
«
К X
<и «
и ю о X
к 3
и
к
X
<и О
Н О
о X X
3
Упругие свойства определяются четырьмя независимыми константами: Е1Ь Е22, С12, У12, представленными в таблице 1.
Таблица 1
Характеристики композитного материала LY5052/T300
Свойства Углеродное волокно Т300 Эпоксид LY5052
Модуль упругости 230 ГПа 3,0 ГПа
Прочность на растяжение 3,5 ГПа 71,0 ГПа
Плотность 1760 г/см3 1,14 г/см3
Силы рассчитываются по следующим уравнениям:
F = ---; F =1 -1; F - о; fu = _J-
1 XX' Y Y ' XX,
; F66 = --±; Fi2 =-
где X — растягивающая сила в продольном направлении; Y — растягивающая сила в поперечном направлении; X — сила давления в продольном направлении; Y — сила давления в поперечном направлении; S — сила сдвига.
Максимальное разрушающее напряжение достигается при выполнении одного из следующих соотношений [12]:
^ > 1; ^ > 1; in > 1.
X
Y
S
Свойства материалов и моделирование методом конечных элементов. Композитный материал, использованный для исследуемого резервуара, представляет собой полимер, армированный углеродным волокном типа Т300, а в качестве полимерного материала применяется эпоксид ЬУ5052. Композитные материалы имеют ортотропную природу, поэтому процесс их моделирования конечными элементами сложнее, чем изотропных материалов, таких как алюминий и сталь.
На рис. 3 представлена конечно-элементная модель резервуара высокого давления из углепластика, внутренний слой которого состоит из алюминиевого сплава, армированного восемью слоями композитного материала СагЬоп/Ероху Т300/ЬУ5052.
Модель имеет следующие размеры:
- длина резервуара —1200 мм;
- диаметр резервуара по центру — 300 мм;
- общая толщина — 64 мм;
- толщина одного слоя — 6,5 мм;
- толщина футеровки — 0,12 мм.
й о тз
л
£
и
<и
fp £ -Й
Рис. 3. Конечно-элементная модель резервуара CFRP
Для исследования резервуаров высокого давления из композитных материалов очень важно выбрать подходящий тип конечных элементов. Программа ANSYS содержит конечные элементы SHELL и SOLID, необходимые для моделирования слоистых композитных материалов. Для исследования использовали
программу SHELL 99, которая ускоряет расчёт конструкции, имеющей до 250 слоев. Этот элемент представляет собой многоуровневую линейную структуру с восемью узлами и шестью степенями свободы и позволяет пользователю определять гибкость, наклон слоев и плотность каждого слоя.
При изучении композитных материалов формирование слоев является одним из основных вопросов, поскольку каждый слой имеет свой угол наклона, а волокна в каждом слое имеют различные углы наклона, поэтому свойства каждого слоя необходимо определять отдельно.
Для формирования слоев требуются характеристики материала, количество слоев, угол наклона волокон, толщина слоя и количество точек интегрирования в каждом слое. Приведем сравнительные сведения о свойствах применяемых материалов: CFRP и алюминия 6061. Плотность: CFRP — 1570 кг/м3, алюминия — 2750 кг/м3. В табл. 2 показаны механические характеристики применяемых материалов: модуль упругости E; предел прочности при сдвиге G; предел прочности при растяжении V.
Таблица 2
Механические свойства CFRP и алюминия 6061
Параметры в направлениях CFRP Алюминий 6061
координатных осей, ГПа
Ex 128 7070
Ey 10,5 70
Gxy 5 70
Gyz 5
Gzx 5
Vxy 0,27 0,3
Vyz 0,4 0,3
Vzx 0,02 0,3
Создана модель с использованием девяти слоёв композитного материала с различными углами наклона волокон:
[+25°/ -25° ] , [+30°/ -30° ], [+35°/ -35° ], [+40°/ -40° ], [+45°/ -45° ], [+50° / -50° ], [+55° / -55° ] .
На рис. 4 показана последовательность слоев с наклоном волокон под углом ± 45°, реализованная в программе ANSYS.
<и К X <и
<и ю о X
к 3
<и К X <и О
а н о о X
к 3
ей
Рис. 4. Последовательность слоев с наклоном волокон под углом ± 45°
111
Цифровое моделирование. Расчет напряжений и перемещений. Работа резервуара высокого давления проанализирована с использованием критериев разрушения по теории ТБа1^и. Для расчета максимальных напряжений и перемещений применялось внутреннее рабочее давление 35 МПа. Наибольшее напряжение резервуара наблюдалось при угле наклона кольцевых волокон — 0° и различных углах наклона спиральных волокон:
[+25°/ -25° ], [+30°/ -30° ], [+35°/ -35° ] , [+40°/ -40° ], [+45°/ -45° ], [+50° / -50° ], [+55°/ -55° ] .
На рис. 5 показано распределение эквивалентных напряжений и перемещений в резервуаре при угле наклона волокон ± 45°. В таблицах 3, 4 приведены максимальные и минимальные значения перемещений и эквивалентных напряжений в направлениях координатных осей X, У, 2.
а) б)
Рис. 5. Эквивалентное распределение в резервуаре: а — напряжений; б — перемещений
Таблица 3
Максимальные и минимальные значения перемещений в направлениях осей X, У, 2
Направление деформации Минимум, мм Максимум, мм
суммарная векторная 0 6,473
по оси X -3,914 3,914
по оси У - 6,473 0
по оси 2 -3,906 3,909
Таблица 4
Максимальные и минимальные значения напряжений в направлении осей X, У, 2
Направление деформации Минимум, МПа Максимум, МПа
суммарная векторная 11,886 463,756
по оси X - 81,852 528,485
по оси У -155,587 341,607
по оси 2 - 80,307 530,07
й о тз
л
£
и
и
£ Л
112
Определение ожидаемого разрушающего давления и оптимального угла наклона волокон. Нагружение резервуара осуществлялось постепенно увеличивающимся внутренним давлением, начиная с его рабочего значения — 35 МПа. Затем полученное расчетное максимальное напряжение внутреннего давления сравнивали
с конечным допускаемым напряжением по условию:
а <а
тах и
где: сттах — расчетное максимальное
напряжение; аи — конечное допускаемое напряжение. Допускаемое напряжение для резервуаров высокого давления из углепластика составляет 1210 МПа.
Модель цилиндрического резервуара рассчитывалась для двух видов путей намотки волокон: кольцевого и спирального, а также при различных углах наклона волокон:
[+25°/ -25° ], [+30°/ -30° ], [+35°/ -35° ] , [+40°/ -40° ], [+45°/ -45° ], [+50° / -50° ], [+55°/ -55° ] .
При кольцевом пути намотки угол наклона волокон к оси цилиндра составляет 0°. На рис. 6 показана конечно-элементная модель цилиндрического резервуара высокого давления.
Рис. 6. Конечно-элементная модель цилиндрического резервуара высокого давления
Результаты исследования. На рис. 7 показано изменение максимального напряжения, начиная со значения угла от 0о и до 45о оно уменьшается, а затем до значения 65о возрастает. Установлено, что максимальное напряжение в цилиндрическом резервуаре высокого давления, армированного углеродными волокнами, было наименьшим при угле наклона волокон ± 45о.
и «
I
к 1) о К л
180 160 140 120 100 80 60 40 20
0 30 35 40 45 50 55 60 65 Угол наклона волокна Рис. 7. Зависимость максимального напряжения от угла наклона волокон
На рис. 8 представлены результаты анализа величины давления для различных углов наклона волокон, которые показывают, что, начиная со значения угла 0о и до 45о давление увеличивается, а затем до значения угла 65о уменьшается. Наибольшее давление, которое может выдержать резервуар, армированный углеволокном без разрушения, составляет 207 Мпа. Это давление имеет место при угле наклона волокон ± 45о, т. о. оно и является ожидаемым значением разрушающего давления для резервуара.
и
к
X
<и «
и ю о X
к 3
и
к
X
<и О
Н О
о X
к 3
0
Ö О T3
л
£
м (U
fp Л
л в
а Д
Угол наклона волокна Рис. 8. Зависимость давления от угла наклона волокон
На рис. 9 представлено распределение напряжений в резервуаре с углом наклона волокон ± 45о, причем максимальное напряжение при приложении давления составляет 1213 МПа и оно выше, чем допускаемое напряжение для резервуаров высокого давления из углепластика, которое составляет 1210 МПа.
MODAL SOLUTION
STEP=1
SUB=1
TIME=1
AN
STWSI TOP RSYS=0 DMX=25.255 SMN=20.52 SMX=1213
(AVG)
30.62 293.186 666.86 818.616 1081
161.853 424.518 687183 949.847 1213
Рис. 9. Распределение напряжений в резервуаре при нагнетании давления и угле наклона волокон ± 45о
Обсуждение и заключения. Конечно-элементная модель резервуара, армированного углеродным волокном, спроектирована и исследована с использованием программы ANSYS с применением в процессе моделирования элемента SHELL 99. Создано несколько моделей с различными углами наклона волокон. Максимальное напряжение и разрушающее давление для резервуара высокого давления рассчитывались с использованием критериев Tsai-Wu. Установлено, что максимальное напряжение оказалась наименьшим при угле наклона волокон ± 45о, а максимально возможное разрушающее давление в этих же условиях составило 207 МПа. Это означает, что оптимальный угол наклона волокон для обеспечения безопасной работы резервуара высокого давления составляет ± 45о, а резервуар из углепластика, рассчитанный с тем же углом намотки волокон, имеет максимальную прочность.
Библиографический список
1. Sankar Reddy, S. Design, Analysis, Fabrication and Testing of CFRP with CNF Composite Cylinder for Space Applications / S. Sankar Reddy, C. Yuvraj, K. Prahlada Rao // International Journal of Composite Materials Structures. — 2015. — Vol. 5. — P. 102-128. https://doi.org/10.5923j.cmaterials.20150505.03.html
2. Subhash N. Khetre. Design and Analysis of Composite High Pressure Vessel with Different Layers using FEA / N.Subhash. P. T. Khetre, Nitnaware, Arun Meshram // International Journal of Engineering Research & Technology (IJERT). — 2014. — Vol. 3. — P. 1460-1464.
3. H. S. da Costa Mattos. Analysis of burst tests and long-term hydrostatic tests in produced water pipelines / H. S. da Costa Mattos, L .M. Paim, J. M. L. Reis // Engineering Failure Analysis. — 2012. — Vol. 22. — P. 128-140. https://doi.org/10.1016/j.engfailanal.2012.01.011
4. Goldin Priscilla, C.P. Effect of statistical scatter in the elastic properties on the predictability of first ply failure of a polymer composite pressure vessel / C. P. Goldin Priscilla, J. Selwin Rajadurai, A. Krishnaveni // Journal of Engineering Research. — 2021. — P. 128-140. https://doi.org/10.36909/jer.11529
5. Chang, R. R. Experimental and theoretical analyses of first-ply failure of laminated composite vessels / R. R. Chang // Composite Structures. — 2000. — Vol. 49 — P. 237-243. https://doi.org/10.1016/S0263-8223(99)00133-6
6. Levend Parnas. Design of fiber-reinforced composite vessels under loading condition / Levend Parnas, Nuran Katuc // Composite Structures. — 2002. — Vol. 58. — P. 83-95. https://doi.org/10.1016/S0263-8223(02)00037-5
7. Burst failure load of composite pressure vessels / Aziz Onder, Onur Sayman, Tolga Dogan, Necmettin Tarakcioglu // Composite Structures. — 2009. — Vol. 89. — P. 159-166. https://doi.org/10.1016/j.compstruct.2008.06.021
8. Stress analysis of non-conventional composite pipe / M. A. Wahab, M. S. Alam, S. S. Pang [et al.] // Composite Structures. — 2007. — Vol. 79. — P. 125-132. https://doi.org/10.1016/j.compstruct.2005.11.054
9. Guedes, R. M. Stress-strain analysis of cylindrical pipe subjected to a transverse load and large deflection / R.M. Guedes // Composite Structures. — 2009. — Vol. 88. — P. 188-194. https://doi.org/10.1016/j.compstruct.2008.03.031
10. Madhavi, M. Design and Analysis of Filament Wound Composite Pressure Vessel with Integrated-end Domes / M. Madhavi, K. V. J. Rao, K. Narayana Rao // Defence Science Journal. — 2009. — Vol. 59. — P. 73-81. http://dx.doi.org/10.14429/dsj.59.1488
11. Tsai, S. W. Composites Design. 3rd ed. / Stephen W. Tsai // Dayton, Ohio : Think Composites, 1987.
12. Pihong Xu. Finite element analysis of burst pressure of composite hydrogen storage vessels / Pihong Xu, Jinyang Zheng, P.F. Liu // Materials & Design. — 2009. — Vol. 30. — P. 2295-2301. https://doi.org/10.1016/j.matdes.2009.03.006
Поступила в редакцию 04.04.2022 Поступила после рецензирования 19.04.2022 Принята к публикации 20.04.2022
Об авторах: и
К
Антибас Имад Ризакалла, доцент кафедры «Основы конструирования машин» Донского государственного К
технического университета (344003, РФ, г. Ростов-на-Дону, пл. Гагарина, 1), кандидат технических наук, ^
доцент, ResearcherID, БсоршГО, ORCID, imad.antvpas@mail.ru §
к
к
Дьяченко Алексей Геннадьевич, доцент кафедры «Основы конструирования машин» Донского д
государственного технического университета (344003, РФ, г. Ростов-на-Дону, пл. Гагарина, 1), кандидат ^
технических наук, БсоршГО, ORCID, alexev-a2@mail.ru К
<о к
Заявленный вклад соавторов: щ
И. Р. Антибас — научное руководство, постановка задачи, определение методологии исследования, сбор и о
сР
анализ аналитических и практических материалов по теме исследования, критический анализ и доработка н
решения, компьютерная реализация решения задачи; А. Г. Дьяченко — анализ научных источников по теме §
исследования, критический анализ и доработка текста. д
й
Конфликт интересов.
Авторы заявляют об отсутствии конфликта интересов.
Все авторы прочитали и одобрили окончательный вариант рукописи.
