Перейти к:
Использование метода конечных элементов для моделирования резервуара высокого давления из полимера, армированного углеродным волокном
https://doi.org/10.23947/2687-1653-2022-22-2-107-115
Аннотация
Введение. За последнее десятилетие спрос в мире на резервуары высокого давления существенно возрос, особенно в таких областях промышленности как авиационная, космическая, химическая и нефтегазовая. Находясь под постоянным воздействием высокого внутреннего давления, стенки резервуаров испытывают повышенное напряжение, что может стать причиной их внезапного разрушения. Для устранения такой возможности и улучшения прочностных характеристик резервуары изготавливают в форме металлических цилиндров с внутренним покрытием из композитного материала, состоящего из смолы, армированной углеродными волокнами. Цель настоящей статьи заключалась в изучении влияния угла наклона углеродного волокна на цилиндрические резервуары и определении величины максимального разрушающего давления с использованием метода конечных элементов программы ANSYS.
Материалы и методы. С помощью программы ANSYS создана конечно-элементная модель резервуара, состоящего из центральной части в виде металлического цилиндра с внутренним покрытием из полимера, армированного углеродными волокнами. По торцам резервуара размещены полусферы со спиральной намоткой. Для моделирования слоистого композитного материала использована командная оболочка SHELL 99, для определения критерия разрушения резервуара высокого давления — теория Tsai-Wu.
Результаты исследования. Модель цилиндрического резервуара рассчитывалась для двух видов намотки волокон: кольцевой и спиральной при различных углах их наклона. Результаты анализа величины давления для различных углов наклона волокон показывают, что, начиная со значения угла 0° и до 45º оно увеличивается, а затем до значения 65º — уменьшается. Наибольшее давление, которое может выдержать резервуар, армированный углеволокном, составляет 207 МПа при угле наклона волокон ± 45º .
Обсуждение и заключения. Анализ исследований показал, что при угле наклона волокон ± 45º максимальное напряжение оказалась наименьшим, а максимально возможное разрушающее давление при том же угле составило 207 МПа. Из этого следует, что оптимальный угол ориентации волокон для обеспечения безопасной работы резервуара высокого давления составляет ± 45º , а резервуар из углепластика при том же угле намотки волокон имеет максимальную прочность.
Ключевые слова
Для цитирования:
Антибас И.Р., Дьяченко А.Г. Использование метода конечных элементов для моделирования резервуара высокого давления из полимера, армированного углеродным волокном. Advanced Engineering Research (Rostov-on-Don). 2022;22(2):107-115. https://doi.org/10.23947/2687-1653-2022-22-2-107-115
For citation:
Antipas I.R., Dyachenko A.G. Using the Finite Element Method to Simulate a Carbon Fiber Reinforced Polymer Pressure Vessel. Advanced Engineering Research (Rostov-on-Don). 2022;22(2):107-115. https://doi.org/10.23947/2687-1653-2022-22-2-107-115
Введение. Использование резервуаров высокого давления из углеродных волокон нашло широкое распространение в различных областях промышленности из-за своих отличительных свойств, таких как малый вес и высокая прочность. Поэтому в последнее время спрос на такие резервуары существенно возрос в тех случаях, когда вес играет большую роль [1][2].
Важнейшими областями применения резервуаров такого типа являются аэрокосмическая, авиационная и химическая техника. Кроме того, резервуары, армированные волокном, широко используются для увеличения напора при транспортировке нефти и газа. Во многих случаях резервуары такого назначения подвергаются воздействию высокого внутреннего давления, следствием чего может стать скачкообразное повышение напряжения на стенках сосуда и внезапное их разрушение, что наносит большой ущерб материальным и людским ресурсам [3][4].
В [5] проведен ряд исследований, направленных только на изучение напорного резервуара из многослойного композита, в ходе которых определена ожидаемая прочность на смятие. В [6] изучалось поведение вращающегося составного напорного резервуара под воздействием внутреннего давления и осевой нагрузки; в [7] проведены исследования по изучению влияния тепловых нагрузок, действующих на многослойный композитный напорный бак. В работе [8] исследовалось поведение многоугольного композитного напорного бака пяти различных форм под действием внутреннего давления различного режима, а в [9] — изучалась работоспособность составного напорного бака при воздействии поперечных нагрузок.
Проектирование резервуара из композитного материала является сложной проблемой, поэтому для проведения полного и точного анализа необходимо выбрать приоритетные факторы. Авторами изучался резервуар высокого давления, армированный несколькими слоями углеродных волокон и подвергаемый нагрузке от внутреннего давления. Цель работы — изучение влияния угла наклона углеродного волокна на цилиндрические резервуары и определение величины максимального разрушающего давления с использованием метода конечных элементов программы ANSYS.
Материалы и методы. Теоретическое исследование. Резервуар высокого давления из композитных материалов (CFRP) состоит из центральной части, представляющей собой металлический цилиндр с внутренним покрытием, состоящим из полимера, армированного углеродными волокнами. Торцевые поверхности резервуара имеют форму полусфер со спиральной намоткой (рис. 1).
Рис. 1. Продольное сечение резервуара
Резервуары, работающие при высоких давлениях и армированные углеродным волокном, изготавливаются методом филаментной намотки (рис. 2). Для получения необходимой стабильности армирования волокна отправляются на движущуюся тележку при тщательном выборе их координации, а затем наматываются на цилиндрическую поверхность. На стабильность координации волокон влияет несколько факторов: температура, форма и обработка поверхности, а также степень сцепления волокон с матрицей. Угол наклона намотки контролируется скоростью перемещения тележки и скоростью вращения цилиндрического барабана. С целью получения высоких эксплуатационных и прочностных свойств внутренняя цилиндрическая поверхность резервуара покрывается несколькими слоями волокон [10].
Рис. 2. Схема намотки волокон
Угол наклона волокон оказывает большое влияние на свойства резервуаров, поэтому важным вопросом является поиск соответствующего угла для каждой части резервуара. Угол ориентации волокон определяется необходимой величиной трения между волокнами и слоем композитного материала, как показано в соотношении:
(1)
где R — расстояние между центром и точкой слоя; R0— радиус центральной оси; Rt1 — радиус по касательной к поверхности цилиндра при [8].
Критерии разрушения Tsai-Wu. При исследовании и моделировании резервуара высокого давления, выполненного из композитного материала, использовались критерии разрушения по теории Tsai-Wu [11]. Реализация уравнения (1) необходима для исследования ожидаемого разрушения ортотропных материалов по теории Tsai-Wu:
Упругие свойства определяются четырьмя независимыми константами: E11, E22, G12, V12, представленными в таблице 1.
Таблица 1
Характеристики композитного материала LY5052/T300
Силы рассчитываются по следующим уравнениям:
где Xt — растягивающая сила в продольном направлении; Yt — растягивающая сила в поперечном направлении; Xc — сила давления в продольном направлении; Yc — сила давления в поперечном направлении; S — сила сдвига.
Максимальное разрушающее напряжение достигается при выполнении одного из следующих соотношений [12]:
Свойства материалов и моделирование методом конечных элементов. Композитный материал, использованный для исследуемого резервуара, представляет собой полимер, армированный углеродным волокном типа T300, а в качестве полимерного материала применяется эпоксид LY5052. Композитные материалы имеют ортотропную природу, поэтому процесс их моделирования конечными элементами сложнее, чем изотропных материалов, таких как алюминий и сталь.
На рис. 3 представлена конечно-элементная модель резервуара высокого давления из углепластика, внутренний слой которого состоит из алюминиевого сплава, армированного восемью слоями композитного материала Carbon/Epoxy T300/LY5052.
Модель имеет следующие размеры:
- длина резервуара —1200 мм;
- диаметр резервуара по центру — 300 мм;
- общая толщина — 64 мм;
- толщина одного слоя — 6,5 мм;
- толщина футеровки — 0,12 мм.
Рис. 3. Конечно-элементная модель резервуара CFRP
Для исследования резервуаров высокого давления из композитных материалов очень важно выбрать подходящий тип конечных элементов. Программа ANSYS содержит конечные элементы SHELL и SOLID, необходимые для моделирования слоистых композитных материалов. Для исследования использовали программу SHELL 99, которая ускоряет расчёт конструкции, имеющей до 250 слоев. Этот элемент представляет собой многоуровневую линейную структуру с восемью узлами и шестью степенями свободы и позволяет пользователю определять гибкость, наклон слоев и плотность каждого слоя.
При изучении композитных материалов формирование слоев является одним из основных вопросов, поскольку каждый слой имеет свой угол наклона, а волокна в каждом слое имеют различные углы наклона, поэтому свойства каждого слоя необходимо определять отдельно.
Таблица 2
Механические свойства CFRP и алюминия 6061
Создана модель с использованием девяти слоёв композитного материала с различными углами наклона волокон:
На рис. 4 показана последовательность слоев с наклоном волокон под углом ±45°, реализованная в программе ANSYS.
Рис. 4. Последовательность слоев с наклоном волокон под углом ±45°
Цифровое моделирование. Расчет напряжений и перемещений. Работа резервуара высокого давления проанализирована с использованием критериев разрушения по теории Tsai-Wu. Для расчета максимальных напряжений и перемещений применялось внутреннее рабочее давление 35 МПа. Наибольшее напряжение резервуара наблюдалось при угле наклона кольцевых волокон — 0° и различных углах наклона спиральных волокон:
На рис. 5 показано распределение эквивалентных напряжений и перемещений в резервуаре при угле наклона волокон ±45°. В таблицах 3, 4 приведены максимальные и минимальные значения перемещений и эквивалентных напряжений в направлениях координатных осей X, Y, Z.
Рис. 5. Эквивалентное распределение в резервуаре: а — напряжений; б — перемещений
Таблица 3
Максимальные и минимальные значения перемещений в направлениях осей X, Y, Z
Таблица 4
Максимальные и минимальные значения напряжений в направлении осей X, Y, Z
Определение ожидаемого разрушающего давления и оптимального угла наклона волокон. Нагружение резервуара осуществлялось постепенно увеличивающимся внутренним давлением, начиная с его рабочего значения — 35 МПа. Затем полученное расчетное максимальное напряжение внутреннего давления сравнивали с конечным допускаемым напряжением по условию: , где:
— расчетное максимальное напряжение;
— конечное допускаемое напряжение. Допускаемое напряжение для резервуаров высокого давления из углепластика составляет 1210 МПа.
Модель цилиндрического резервуара рассчитывалась для двух видов путей намотки волокон: кольцевого и спирального, а также при различных углах наклона волокон:
При кольцевом пути намотки угол наклона волокон к оси цилиндра составляет 0°. На рис. 6 показана конечно-элементная модель цилиндрического резервуара высокого давления.
Рис. 6. Конечно-элементная модель цилиндрического резервуара высокого давления
Результаты исследования. На рис. 7 показано изменение максимального напряжения, начиная со значения угла от 0º и до 45º оно уменьшается, а затем до значения 65º возрастает. Установлено, что максимальное напряжение в цилиндрическом резервуаре высокого давления, армированного углеродными волокнами, было наименьшим при угле наклона волокон ±45º.
Рис. 7. Зависимость максимального напряжения от угла наклона волокон
На рис. 8 представлены результаты анализа величины давления для различных углов наклона волокон, которые показывают, что, начиная со значения угла 0º и до 45º давление увеличивается, а затем до значения угла 65º уменьшается. Наибольшее давление, которое может выдержать резервуар, армированный углеволокном без разрушения, составляет 207 Мпа. Это давление имеет место при угле наклона волокон ±45º, т. о. оно и является ожидаемым значением разрушающего давления для резервуара.
Рис. 8. Зависимость давления от угла наклона волокон
На рис. 9 представлено распределение напряжений в резервуаре с углом наклона волокон ±45º, причем максимальное напряжение при приложении давления составляет 1213 МПа и оно выше, чем допускаемое напряжение для резервуаров высокого давления из углепластика, которое составляет 1210 МПа.
Рис. 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.5923.j.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
Об авторах
И. Р. АнтибасРоссия
Антибас Имад Ризакалла, доцент кафедры «Основы конструирования машин», кандидат технических наук, доцент
г. Ростов-на-Дону, пл. Гагарина, 1
А. Г. Дьяченко
Россия
Дьяченко Алексей Геннадьевич, доцент кафедры «Основы конструирования машин», кандидат технических наук
г. Ростов-на-Дону, пл. Гагарина, 1
Рецензия
Для цитирования:
Антибас И.Р., Дьяченко А.Г. Использование метода конечных элементов для моделирования резервуара высокого давления из полимера, армированного углеродным волокном. Advanced Engineering Research (Rostov-on-Don). 2022;22(2):107-115. https://doi.org/10.23947/2687-1653-2022-22-2-107-115
For citation:
Antipas I.R., Dyachenko A.G. Using the Finite Element Method to Simulate a Carbon Fiber Reinforced Polymer Pressure Vessel. Advanced Engineering Research (Rostov-on-Don). 2022;22(2):107-115. https://doi.org/10.23947/2687-1653-2022-22-2-107-115