Подать статью
Стать рецензентом
Том 245
Страницы:
561-568
Скачать том:
RUS ENG

Численное моделирование двустенного шарового резервуара

Авторы:
М. Г. Каравайченко1
Л. И. Газалеев2
Об авторах
  • 1 — д-р техн. наук профессор Уфимский государственный нефтяной технический университет
  • 2 — инженер ЗАО «Нефтемонтаждиагностика»
Дата отправки:
2020-05-29
Дата принятия:
2020-09-16
Дата публикации:
2020-11-24

Аннотация

Обширным и важным классом многослойных конструкций оболочек являются трехслойные конструкции. В трехслойной конструкции важную роль играет жесткий заполнитель, за счет которого разнесены несущие слои, что придает пакету слоев высокую жесткость и прочность при относительно малом весе. Комбинируя толщины несущих слоев и заполнителя, можно добиться необходимых свойств трехслойной конструкции оболочки. По сравнению с традиционными одностенными трехслойная конструкция обладает повышенной жесткостью и прочностью, что позволяет уменьшить толщину и массу оболочек. С целью снижения металлоемкости конструкции шарового резервуара для хранения сжиженных газов в настоящей работе рассмотрена конструкция двустенного резервуара, в котором межстенное пространство заполнено армированным пенополиуретаном. Численное моделирование позволило определить параметры напряженно-деформированного состояния конструкции с погрешностью не более 5 %. Установлено на примере резервуара объемом 4000 м 3 , что пространственная конструкция стенки шарового резервуара позволяет снизить металлоемкость до 19 %. Областью применения результатов исследования может быть оценка напряженно-деформированного состояния шаровых резервуаров при их проектировании. Разработана методика построения конструкции двустенного шарового резервуара в программном комплексе SCAD, позволяющая произвести расчет напряженно-деформированного состояния (НДС) методом конечных элементов. Разработана численная модель двустенного шарового резервуара. Установлено, что для получения результатов расчета с погрешностью P ≤ 5 % размер конечного элемента не должен превышать 300×300×δ мм. В результате исследования конструкции двустенного шарового резервуара установлены конструктивные параметры, позволяющие обеспечить эксплуатационную надежность сооружения при снижении металлоемкости по сравнению с одностенным резервуаром на 19 %

Ключевые слова:
ортотропная оболочка шаровой двустенный резервуар численное моделирование двустенного шарового резервуара напряжения и перемещения
10.31897/PMI.2020.5.8
Перейти к тому 245

Введение. Основные подходы к изучению ортотропных оболочек были рассмотрены В.З.Власовым [12]. В работе [4] А.И.Голованов с соавторами рассмотрели проблемы построения моделей трехслойных оболочек методом конечных элементов (МКЭ). Разработали квадратичный конечный элемент, который распространяется на класс физически и геометрически нелинейных задач. Использован метод пошагового нагружения в форме модифицированной лагранжевой постановки.

В своей работе [9] Р.З.Киселева предложила алгоритм формирования матриц жесткости объемных конечных элементов (КЭ) с узловыми неизвестными в виде перемещений и их первых производных и алгоритм преобразования матриц жесткости и векторов узловых нагрузок КЭ, примыкающих к границе раздела слоев оболочек вращения с различными физико-механи¬ческими свойствами.

Ежегодно количество работ, исследующих трехслойные конструкции, увеличивается. Однако большая часть из них посвящена рассмотрению пластин для прямолинейной плоскости. Многие научные публикации посвящены многослойным оболочкам для плоских покрытий с различными характеристиками жесткости слоев [5, 15].

Наибольшее распространение в практических расчетах нашли дифференциальные теории: теория упругопластических процессов А.А.Ильюшина и теории пластического течения при комбинированном упрочнении, базирующиеся на концепции микронапряжений В.В.Новожилова.

В настоящее время существует множество вариантов моделей упругопластического деформирования металла, основанных на теории течения при его упрочнении. Предполагается, что металл однороден и начально изотропен. В процессе упругопластического деформирования в нем может возникнуть только деформационная анизотропия. В работе [1] представлен обзор работ, основанных на теории течения. Под научным руководством академика РАЕН В.С.Бондаря выполнена работа [11], в которой приведен вариант математической модели неупругого деформирования и накопления повреждений конструкционных материалов. Разработан программный комплекс оценки и прогнозирования ресурса конструкционных материалов, состоящий из модуля определения материальных функций по результатам базового эксперимента и модуля расчета кинетики напряженно-деформированного состояния и процесса накопления повреждений. В работе [14] предложен прикладной вариант теории упругопластических процессов и накопления повреждений материалов, основанный на теории В.С.Бондаря. Автор апробировал вариант теории на ряде конструкционных материалов, указал, что ошибка результатов расчета не превышает 10 % от результатов экспериментов, однако данную теорию можно применять только при малых деформациях. В теории оболочек под малыми деформациями понимают деформации, не превышающие толщины стенки.

Т.Р.Ищанов [8] анализировал напряженно-деформированное состояние тонких оболочек с учетом поперечного сдвига при различных вариантах аппроксимации угловых перемещений. На основе соотношений В.В.Новожилова разработан алгоритм исследования НДС при использовании скалярной и векторной интерполяционных процедур. В работе [3] приведен обзор работ, посвященных исследованию методов сведения трехмерной задачи теории упругости к двумерной теории пластин и оболочек. Н.А.Гуреева [6] выполнила сравнительный анализ эффективности использования конечных элементов треугольной и четырехугольной форм расчета оболочек вращения. Показано, что расчеты, выполненные с использованием треугольных КЭ, обладают существенно большей погрешностью по сравнению с конечно-элементными решениями, полученными на основе использования четырехугольных КЭ. Работа [10] посвящена качественному анализу метода конечных суперэлементов (МКСЭ) Федоренко. Получены априорные оценки погрешностей МКСЭ в шкале пространств Соболева на примере задачи Дирихле.

В химической, нефтехимической, нефтеперерабатывающей и других отраслях промышленности России для хранения жидкостей и газов широко применяют шаровые резервуары. За последние годы наблюдается увеличение производственных объемов промышленных предприятий, а с этим возрастает площадь территорий объектов. Но размеры промышленных территорий ограничены в своих параметрах, следовательно, необходимо увеличить единичную емкость шаровых резервуаров, что повлечет за собой уменьшение площади резервуарных парков, снижение металлоемкости на единицу объема хранимого продукта, снижение затрат на монтажные работы, на установку контрольно-измерительных приборов, сооружение технологических трубопроводов, а также повлияет на снижение потерь продукта от испарения.

Стенка резервуара является важной частью резервуарной конструкции. На оболочку резервуара действуют значительные нагрузки и воздействия. От выбора стали и параметров стенки зависит состояние резервуара в процессе строительства и дальнейшей эксплуатации. С увеличением вместимости резервуара для обеспечения прочности и устойчивости стенки необходимо либо применять сталь с улучшенными характеристиками, либо увеличивать толщину стенки, при этом учитывая экономическую выгоду.

В настоящее время при проектировании шаровых резервуаров в России руководствуются нормативными документами – ГОСТами* и сводами правил**, расчет внешних нагрузок производят по своду правил***. Актуальность данной работы заключается в том, что вопрос работоспособности шаровых резервуаров для хранения углеводородного сырья с учетом ортотропности физических свойств оболочек изучен недостаточно.

Конструкция двустенного шарового резервуара. В настоящей работе рассмотрен вариант двустенной конструкции (рис.1), состоящей из внешних слоев малой толщины и системы ферм, играющей роль слоя-заполнителя.

Рис.1. Общий вид трехслойного шарового резервуара объемом 4000 м3 1 – внутренняя оболочка резервуара; 2 – ферменный заполнитель; 3 – хранимый продукт; 4 – купол; 5 – холодильно-газовая машина; 6 – поплавковый указатель уровня; 7 – совмещенный узел для замера уровня, температуры нефти и отбора пробы; 8 – площадка обслуживания; 9 – внешняя оболочка резервуара; 10 – полка; 11 – опорная стойка; 12 – опорные связи; 13 – шахтная лестница; 14 – приемный и раздаточный патрубки; 15 – днище

Внутренний слой резервуара представляет собой оболочку, собранную из отдельных стальных листов, раскроенных на заводе из прокатных листов, а затем деформированных для придания сферической кривизны. Радиус кривизны листов равен радиусу оболочки. Внутренний слой выполнен из стали ОН9 с 9 %-м содержанием никеля, отличающейся высокой упругостью, необходимой для хранения криогенных жидкостей, с временным сопротивлением σв = 640 МПа и пределом текучести σT = 595 МПа.

Заполнитель представляет собой систему ферм, которая при выпучивании пластинок остается несгибаемой и создает для них дополнительные опорные закрепления. Деформации стенки внутри каждого такого поля рассматривают как самостоятельные. Ферменный заполнитель прикреплен к внешнему и внутреннему слою уголками, которые привариваются к оболочке в местах сварки листов внутренней оболочки.

Внешний слой шарового резервуара собирают из стальных листов длиной 6 и шириной 2 м так, чтобы место стыка стальных листов приходилось на балку жесткости, представляющую собой уголок. Таким образом, создано сварное соединение для односторонней сварки на подкладке, где подкладкой служит полка уголка. Внешняя оболочка резервуара опирается на опоры из труб, соединенные между собой связями. Внешний слой шарового резервуара изготовлен из стали класса прочности С345 марки 09Г2С. Особенностью этой марки стали является низкий порог хладноломкости, хорошая свариваемость и необходимая прочность.

Геометрические параметры внутреннего, наружного и срединного слоя стенки резервуара определяются расчетом в зависимости от заданных нагрузок [2].

Создание численной модели в программном комплексе SCAD [7]. Модели стенок резервуара разбили на конечные элементы – призмы размером 0,3×0,3×δi, где толщина листов стенки резервуара δi = 14 мм. Срединный слой принимали как стержневые элементы длиной конечного элемента l = 0,42 м. Опорные стойки приняли как стержни наружным диаметром 0,3 м с толщиной стенки 0,02 м.

Опорные столики, соединяющие стенку резервуара с опорными стойками, приняли как листовые элементы с дискретизацией на КЭ размером 0,1×0,1×0,02. В нижние части опорных стоек ввели связи по осям X, Y, Z ортогональной системы координат.

Сбор нагрузок и воздействий. Первыми шагами расчета трехслойного шарового резервуара является определение внешних нагрузок, построение расчетной схемы, а также установление характера (формы) эпюры приложения этих нагрузок. Можно выделить нагрузки, действующие на трехслойный шаровой резервуар: собственный вес металлоконструкции трехслойного шарового резервуара; снеговая нагрузка, рассматриваемая в трех вариантах; ветровая нагрузка; температурные климатические воздействия; нагрузка от избыточного давления; нагрузка от гидростатического давления. Расчет вариантов сочетания нагрузок* показал, что максимальная нагрузка, действующая на шаровой резервуар, составляет 232,5 кПа.

Рис.2. График зависимости относительных перемещений f/[f] и относительных напряжений σ/[σ] от толщин стенки оболочки δ
Рис.3. График зависимости относительных напряжений σ/[σ] и устойчивости стержней σкр от параметров стержня ферменного заполнителя (a) и шага фермы (б)

Для исследования применялись безразмерные величины: относительная величина напряжений σ/[σ], возникающая в сечении элемента под воздействием нагрузок и воздействий, где [σ] –допустимые напряжения, равные: для стали марки 09Г2С 190 МПа, для стали марки ОН9 – 315 МПа [11]

В модели резервуара принято, что внутренняя и наружная стенки имеют равные толщины. Исследование влияния толщины стенок на величины перемещений и напряжений показано на рис.2. По зависимостям σ/[σ] = F (δ) и f/[f] = F (δ) видно, что обе функции практически линейно убывающие. Исследуя поведение двух функций, приходим к выводу, что оптимальной толщиной оболочек с учетом запаса на коррозию δ = 14 мм.

В связи с тем, что в резервуаре будет храниться продукт с температурой Т = –160 oС, необходимо обеспечить теплоизоляцию. В качестве теплоизоляции принят слой пенополиуретана толщиной В = 400 мм с характеристиками: плотность ρ = 150 кг/м3, коэффициент теплопроводности λ = 0,3∙106 Вт/м∙К. Для исключения мостиков холода связи в фермах выполнены из стеклопластика марки ВПС-30 с характеристиками: плотность ρ = 1800 кг/м3, модуль упругости Е = 55∙103 МПа, коэффициент линейного расширения \(\ \alpha_t=8,3\cdot106{^oC^{-1}} \). На рис.3, а представлены зависимости прочности и устойчивости стержней ферменного заполнителя от их параметров. По зависимости можно сделать вывод об уменьшении относительных напряжений в ферменном заполнителе в связи с увеличением диаметра стрежня ферменного заполнителя. Исследуя поведение зависимости, приходим к выводу, что оптимальным параметром стрежня ферменного заполнителя будет являться 120×12 мм.

Из конструктивных соображений – рационального очертания фермы и удобства прикрепления – желателен угол, равный 45°. При малых углах ферма становится слишком вытянутой, при больших – высокой, что делает их громоздкими и неэкономичными. Исследуя зависимость относительных напряжений от угла между стенкой и стержнем ферменного заполнителя, можно утверждать, что угол, равный 45°, оптимален для ферменного заполнителя.

Размер шага, с каким следует устанавливать фермы, можно вывести согласно зависимости относительных напряжений и устойчивости стержней от шага расстановки ферм (рис.3, б). Можно утверждать, что оптимальным шагом будет являться расстояние не более 400 мм.

Рис.4. Распределение перемещений внешней стенки в двустенном шаровом резервуаре

Обсуждение результатов. Результатом расчета являются цветовые схемы распределения перемещений внешней стенки, усилий и напряжений в двустенном шаровом резервуаре. Результаты расчета представлены на рис.4, 5. На рис.4 видно, что максимальные перемещения во внешней стенке шарового резервуара не превышают 90 мм.

Предельно допустимые перемещения для стальных конструкций пролетом до 30 м, согласно своду правил*: D/200 = 20000/200 = 100 мм, где D – диаметр двустенного шарового резервуара, равный 20000 мм. Согласно расчету двустенного шарового резервуара, перемещения конструкции меньше 100 мм, что является допустимым.

Распределение главных и эквивалентных напряжений во внутреннем слое резервуара представлено на рис.5, а. Из рис.5, а видно, что в пластинах оболочки напряжения сжатия достигают 111,9 МПа, а в ребрах жесткости напряжения растяжения равны 173,5 МПа.

На рис.5, б представлено распределение главных и эквивалентных напряжений в срединном слое резервуара. Максимальные напряжения растяжения в срединном слое достигают 150,9 МПа, а максимальные напряжения сжатия 55,01 МПа. Причем преимущественно напряжения растяжения наблюдаются в нижнем поясе срединного слоя, где стержни примыкают к днищу оболочки.

В таблице представлены сравнительные результаты расчетов, полученные МКЭ с аналитическим методом расчета. Видно, что погрешность между результатами расчетов, полученных численным (МКЭ) и аналитическим методами, не превышают 4,35 %.

Сравнение полученных результатов

Аналитический расчет Численный расчет Погрешность, %
Максимальные перемещения общей конструкции, мм
93,1 89,58 4,07
Максимальные эквивалентные напряжения внутреннего слоя, МПа
180,78 173,56 4,16
Аналитический расчет Численный расчет Погрешность, %
Максимальные эквивалентные напряжения срединного слоя, кПа
150,9 157,46 4,35
Максимальные эквивалентные напряжения внешнего слоя, МПа
164,00 157,5 4,13

Для исследования металлоемкости трехслойного шарового резервуара необходим сравнительный объект – однослойный шаровой резервуар с идентичным объемом. Масса оболочки трехслойного шарового резервуара с рассчитанными параметрами равна 110 т, а масса оболочки однослойного шарового резервуара с толщиной стенки 40 мм, изготовленного в АО «Уралхиммаш», составляет 136 т. При сравнении масс двустенного и одностенного шаровых резервуаров масса разработанной конструкции меньше на 19 %.

Рис.5. Распределение главных и эквивалентных напряжений во внутреннем (а) и в срединном слое (б) двустенного шарового резервуара

Заключение. Представленные исследования проведены в рамках линейной постановки задачи. По предложенной методике расчета можно определить оптимальные сечения внутренней и наружной оболочки, сечения связей между оболочками, т.е. решить задачу в первом приближении. Результаты исследования показали перспективность многослойных конструкций в резервуаростроении.

Однако на практике невозможно изготовить оболочку с идеальной геометрией. Так в шаровом резервуаре для хранения углеводородного сырья имеет место конструктивная нелинейность, вызванная наличием врезок патрубков, люков, опорных элементов. Геометрическую нелинейность характеризуют остаточные сварочные деформации и ошибки монтажников при сборке резервуара из отдельных листов. В настоящее время существует методика лазерного сканирования крупных резервуаров [13], позволяющая оценить отклонения стенки от проектного положения с высокой разрешающей способностью. Физическая нелинейность в стенке резервуара вызвана формообразованием листовых оболочечных элементов и сваркой их между собой. Таким образом, научный и практический интерес вызывает вопрос оценки напряженно-деформированного состояния оболочек в условиях нелинейности при упруго-пластичном состоянии металла. Это станет следующим объектом наших исследований.

Выводы. Разработана методика построения конструкции двустенного шарового резервуара в программном комплексе SCAD, позволяющая произвести расчет НДС методом конечных элементов.

Разработана численная модель двустенного шарового резервуара. Установлено, что для получения результатов расчета с погрешностью P ≤ 5 % размер конечного элемента не должен превышать 300×300×δ мм.

Исследована конструкция двустенного шарового резервуара. Установлены конструктивные параметры, позволяющие обеспечить эксплуатационную надежность сооружения при снижении металлоемкости по сравнению с одностенным резервуаром на 19 %.

Литература

1. Абашев Д.Р. Развитие модели упругопластического деформирования, критериев усталости и методик идентификации материальных параметров конструкционных сплавов: Автореф. дис. … канд. физ.-мат. наук. Королев: Центральный научно-исследовательский институт машиностроения, 2016. 21 с.

2. Абдуллин И.Н. Моделирование ферменного заполнителя трехслойной конструкции // Поиск эффективных решений в процессе создания и реализации научных разработок в российской авиационной и ракетнокосмической промышленности: Международная научно-практическая конференция «Поиск эффективных решений в процессе создания и реализации научных разработок в российской авиационной и ракетнокосмической промышленности», Казань, 5-8 августа 2014. С. 307-312.

3. Аннин Б.Д. Неклассические модели теории пластин и оболочек / Б.Д.Аннин, Ю.М.Волчков // Прикладная механика и техническая физика. 2016. Т. 57. № 5. С. 5-14. DOI: 10.15372/PMTF20160501

4. Голованов А.И. Метод конечных элементов в статике и динамике тонкостенных конструкций / А.И.Голованов, О.Н.Тюленева, А.Ф.Шигабутдинов. М.: Физматлит, 2006. 392 с.

5. Горшков А.Г.Механика слоистых вязкоупругопластических элементов конструкций / А.Г.Горшков, Э.И.Старовойтов, А.В.Яровая. М.: Физматлит, 2003. 577 с.

6. Гуреева Н.А. Сравнительный анализ эффективности использования конечных элементов треугольной и четырехугольной формы в расчетах оболочек. Автореф. дис. … канд. техн. наук. Волгоград: Волгоградская государственная сельскохозяйственная академия, 2004. 28 c.

7. Зеливянский Е. Переход от архитектурного решения к расчетной схеме / Е.Зеливянский, В.Карпиловский, Э.Криксунов // SCADmaster, 2000. № 5. С. 40-43.

8. Ищанов Т.Р. Конечно-элементный анализ напряженно-деформированного состояния тонких оболочек с учетом поперечного сдвига при различных вариантах аппроксимации угловых перемещений: Автореф. дис. … канд. техн. наук. Волгоград: Волгоградский государственный аграрный университет, 2018. 21 с.

9. Киселева Р.З. Расчет слоистых пластин и оболочек вращения на основе трехмерных конечных элементов без предположения о деформировании нормали: Автореф. дис. ... канд. техн. наук. Волгоград: Волгоградская государственная сельскохозяйственная академия, 2010. 20 с.

10. Лазарева С.А. Аппроксимационные свойства метода конечных суперэлементов Федоренко // Вычислительные технологии. 2008. Т. 13. Специальный выпуск 4. С. 75-81.

11. Макаров Д.А. Математическое моделирование процессов неизотермического неупругого деформирования и накопления повреждений в конструкционных материалах. Автореф. дис. … канд. физ.-мат. наук. М., Московский автомеханический институт, 2005. 21 с.

12.Общая нелинейная теория упругих оболочек / С.А.Кабриц, Е.И.Михайловский, П.Е.Товстик, К.Ф.Черных, В.А.Шамина. СПб.: Изд-во Санкт-Петербургского университета, 2002. 388 с.

13. Сальников А.П. Оценка напряженно-деформированного состояния резервуаров по состоянию наземного лазерного сканирования. Автореф. дис. … канд. физ.-мат. наук. М.: Российский государственный университет нефти и газа им. И.М.Губкина, 2016. 22 с.

14. Семенов П.В. Прикладной вариант теории упругопластических процессов и накопления повреждений материалов. Автореф. дис. … канд. физ.-мат. наук. М.: Московский автомеханический институт, 2013. 21 с.

15. Abot J.L. Contact Law for Composite Sandwich Beams / J.L.Abot, I.M.Daniel, E.E.Gdoutos // Journal of Sandwich Structures and Materials. 2002. Vol. 4. № 2. P. 157-173. DOI: 10.1177/1099636202004002705

Похожие статьи

Методика расчета пневмокомпенсаторов для плунжерных насосов с погружным приводом
2020 Э. О. Тимашев
Пространственная нелинейность динамики метановыделения в подземных скважинах для устойчивого развития геотехнологий
2020 А. К. Джиоева, В. С. Бригида
Геологическое строение северной части Карского шельфа у архипелага Северная Земля по результатам последних исследований
2020 Е. А. Гусев, А. А. Крылов, Д. М. Урванцев, Ю. В. Горемыкин, П. И. Криницкий
Методика управления процессом бурения и экспериментальные исследования сил сопротивления при бурении долотами с резцами PDC
2020 В. В. Нескоромных, М. С. Попова, А. Е. Головченко, П. Г. Петенёв, Лиу Баочанг
Фальбанды Керетского архипелага Белого моря: характеристика состава пород и минералов, рудная минерализация
2020 Л. И. Салимгараева, С. Г. Скублов, А. В. Березин, О. Л. Галанкина
Изучение влияния алюминия на адгезию нержавеющей стали при газопламенном напылении
2020 С. Н. Степанов, Т. А. Ларионова, С. С. Степанов