doi: 10.5862/MCE.48.8
Вибрация шельфовых сооружений при действии льда
Инженер-гидротехник Е.А. Салганик,
ЗАО «ГТ Морстрой»; д.ф.-м.н., профессор К.Н. Шхинек,
Санкт-Петербургский государственный политехнический университет
Аннотация. Целью статьи является изучение процесса динамического взаимодействия льда и шельфовых сооружений и определение условий возникновения вибраций. Для этого решены задачи по созданию полноразмерной конечно-элементной модели взаимодействия ровного льда и шельфовых сооружений, сравнению полученных результатов условий возникновения вибраций с результатами полевых исследований.
Моделирование проводилось с помощью программы конечно-элементного анализа ANSYS. Для создания модели льда с возможностью разрушения были решены задачи моделирования процесса трещинообразования, проверки полученной многокомпонентной модели с помощью базовых тестов по определению ее прочностных свойств, определению уровня нагрузок на сооружения различного типа.
Объектом исследования является роль упругих и инерционных параметров шельфовых сооружений, геометрических и динамических параметров ледового поля в процессе возникновения вибраций при взаимодействии сооружения и льда. В качестве предметов для моделирования сооружений были выбраны платформа «Моликпак» в Канаде и маяк «Норстрёмсгрюнд» в Швеции.
Исследование содержит апробацию результатов, их сравнение с данными полевых исследований, с величинами нагрузок, полученных с помощью рекомендаций стандарта ISO 19906 «Нефтяная и газовая промышленность. Сооружения арктического шельфа» и с результатами дискретно-элементного моделирования. Полученные результаты имеют приемлемый уровень погрешности по сравнению c результатами других исследований процесса возникновения динамических нагрузок при действии льда.
Ключевые слова: вибрация; ледовые нагрузки; шельфовые сооружения; механика разрушения; компьютерное моделирование; метод конечных элементов; ANSYS
Введение
Растущий объем деятельности человека в арктической зоне требует точной методики определения нагрузок на сооружения. Ледовые нагрузки в этой зоне превышают значения суммы всех остальных нагрузок. Обычно ледовые нагрузки рассматриваются со стороны ровного льда, торосов и айсбергов [1]. В данной статье рассматриваются динамические нагрузки со стороны ровного льда.
За последние несколько десятилетий значительное количество гидротехнических сооружений подвергалось влиянию динамических нагрузок со стороны ровного льда, превышающих расчетные [2]. Одним из наиболее значительных событий было возникновение вибрации нефтедобывающей платформы «Моликпак» в море Бофорта в апреле 1986 года [3].
Современные стандарты определения нагрузок на шельфовые сооружения со стороны ровного льда предполагают определение динамических нагрузок. Международный стандарт ISO 19906 [4] предлагает методы их нахождения. Однако величины нагрузок, полученные по этому стандарту, нередко бывают превышены на различных гидротехнических сооружениях [1]. Кроме того, методика нахождения этих нагрузок по стандарту ISO не учитывает многие значимые для динамического воздействия параметры, такие как скорость ледового поля, форма гидротехнического сооружения, его жесткостные и инерционные параметры.
Целью этого исследования является изучение процессов динамического взаимодействия льда и шельфовых сооружений на основе конечно-элементного моделирования. Анализ литературы показывает, что с помощью этого метода пока что не было описано возникновение вибраций шельфовых сооружений [5].
Методология исследования основана на использовании метода конечных элементов для моделирования динамического взаимодействия льда и сооружения с учетом образования внутренних трещин во льду. Для процесса трещинообразования использовалась модель связывающих элементов.
Статья состоит из описания исследования свойств модельного льда, их соотношения с реальными физическими параметрами и расчета вибрации сооружений. Соотношение модели с физическими параметрами льда состоит в проведении базовых тестов для используемой модели льда для определения достоверности основных характеристик материала: прочность на одноосное сжатие, прочность на одноосное растяжение, прочность на изгиб. В третьей части исследования определены глобальные статические и динамические нагрузки, которые были сравнены с нагрузками, определяемыми с помощью стандарта ISO 19906.
Обзор иностранной и отечественной литературы
Ранние исследования вибрации сооружений показывают, что при рассмотрении взаимодействия льда и сооружения следует учитывать несколько диапазонов скоростей, которые и определяют характер взаимодействия [1].
Бленкарн и Пейтон провели натурные эксперименты, которые показали, что колебания сооружения происходят при малых скоростях льда [6]. Также они получили зависимость прочности льда на одноосное сжатие от скорости деформации льда.
Карна [7] с помощью лабораторных экспериментов определил диапазон скоростей, при которых возникают вибрации. Он пришел к выводу, что при возникновении вибраций максимальная скорость сооружения превышает скорость льда на 40 %.
Уравнения, полученные Пейтоном, и модель вибраций Маттанена использовались для анализа динамических взаимодействий на сооружения [8].
С помощью полевых измерений на различных гидротехнических сооружениях были получены условия возникновения вибраций платформ в Бохайском заливе [9], заливе Кука [10], море Бофорта [3], Ботническом заливе [11].
Показано, что возникновение вибраций зависит от соотношения между скоростью льда, его прочностными и геометрическими параметрами, размерами сооружения, его жесткостью и массой [2].
Существенный шаг в решении вопроса возникновения вибраций сделан в работах [5, 12], описывающих процесс разрушения льда через образование трещин. Для описания связи усилий и деформаций в работе [13] использована модель Твергаарда и Хатчинсона, описывающая связь усилий и деформаций в предварительно созданных трещинах. Также использовались положения механики разрушения, разработанные Барренблаттом [14].
В области статического и динамического конечно-элементного моделирования работали Дайан [15] и Санд [16].
Для задачи описания процесса динамического взаимодействия сооружения и льда Жиленковым [17] и Шхинеком [18, 19] была создана дискретно-элементная модель.
Динамическая нагрузка в ISO 19906 определяется из условий воздействия на сооружение силы, величина которой определяется с помощью идеализированной кривой и зависит от времени [4]. Нагрузка задается в определенном диапазоне значений в виде повторяющихся треугольников с периодом, равным собственной частоте колебаний сооружения Fstr.
Постановка задачи
Рассматривается следующая задача: на податливое сооружение, расположенное на дне моря, действует ледяное образование. Задача решается путем изучения взаимодействия сооружения со льдом, определения нагрузок, действующих на сооружение, и возникающих при этом усилий и реакции сооружения.
В качестве программного комплекса для моделирования выбран комплекс ANSYS. Для корректного описания свойств льда должен быть проведен комплекс исследований. Эти исследования касаются вопросов развития трещин между элементами и связи параметров отдельных элементов модели льда и параметров конечно-элементной сетки с глобальными прочностными свойствами льда, такими как модуль упругости, коэффициент Пуассона, пределы прочности на растяжение, сжатие и изгиб в различных направлениях.
Для подтверждения правильности моделирования процесса необходимо сравнить полученные результаты с данными полевых наблюдений и рекомендациями нормативных документов.
Механические свойства льда
Морской лед в общем случае является анизотропным и нелинейным материалом. При динамических нагрузках он описывается основными параметрами: пределом прочности на сжатие ау и растяжение Т, модулем Юнга E, коэффициентом Пуассона V и коэффициентом кулоновского трения ^ [1].
При высоких значениях скорости деформации а, которым соответствует процесс взаимодействия льда и гидротехнических сооружений, лед разрушается путем образования трещин. При увеличении нагрузки лед сначала деформируется упруго, затем пластично при параллельно идущем процессе образования трещин.
В основе линейно-упругой механики разрушения лежит предположение о существовании доминирующих трещин, роль которых в моделировании льда выполняют связывающие элементы. Для характеристики трещины существует коэффициент интенсивности напряжений К.
Согласно экспериментам [20] зарождение трещин при высоких скоростях деформации происходит при растягивающих напряжениях до 1,2 МПа, развитие трещин - после этого предела.
1/2
Критический коэффициент интенсивности напряжений К|О составляет 100-250 кПа/м [4]. Большие значения соответствуют морскому льду, меньшие - лабораторным образцам.
Интенсивность выделения энергии О1 в нормальном направлении определяется по формуле
Ирвина [20]:
0 = = к (1 -V2) (1)
1 Е' Е ' где Е - модуль Юнга льда, V - коэффициент Пуассона льда.
Параметры конечно-элементной модели
Для моделирования принимаются следующие параметры льда: модуль Юнга E = 5 ГПа; коэффициент Пуассона V = 0,3; предел прочности на сжатие О = 1,5-3,0 МПа; предел прочности на растяжение Т = 0,3-1,1 МПа [22]. Для льда используются полноинтегрируемые 8-узловые объемные элементы и упруго-пластичный материал [23, 24] (табл. 1).
Таблица 1. Входные параметры для упруго-пластичного материала объемных элементов льда [25]
Параметр Единица измерения Значение параметра
Массовая плотность р кг/м3 910
Модуль Юнга Е ГПа 5,0
Коэффициент Пуассона V - 0,3
Предел текучести а у МПа 1,0-3,0
Модуль упругого упрочнения Е1ЭП МПа -6,7; 0; 1,0; 5,0; 10,0
Кривая зависимости напряжения от пластической деформации а = а(£р ) - Рисунок 1
При малых напряжениях материал ведет себя как упругий с линейной зависимостью напряжений и деформаций, характеризуемой модулем Юнга Е. Для перехода к пластичности напряжения в материале должны превысить значения предела текучести. Далее, если задана кривая зависимости напряжения а от пластической деформации ар (рис. 1), напряжения в
материале линейно уменьшаются (разупрочнение) при росте деформации до определенного значения напряжения (0,01 МПа), после чего материал становится идеально пластичным.
(Т. МПа
10
Пластиче разупрочн :кое гние
Вя-жид кая f юсть 1
0.2
0.4
0.6
0.8
Рисунок 1. Кривая зависимости напряжения а от пластической деформации ар
Если задан модуль упругого упрочнения Е13п, после преодоления предела текучести наступает фаза линейного упрочнения. В этом случае зависимость напряжения и деформации материала билинейная. Для корректного моделирования процесса образования трещин модель льда составляется из основных и связывающих элементов (рис. 3).
Рисунок 2. Схема образования трещины при использовании связывающих элементов: при росте деформации б растет тяга Т(х) до момента достижения деформацией
критического значения бс
При внедрении связывающих элементов необходимо задать закон связи тяги и деформации. В качестве этого закона удобно использовать закон Твергаарда и Хатчинсона [13] (рис. 4).
Рисунок 3. Разделение конечно-элементной модели льда на вертикальные и горизонтальные связывающие элементы и объемные элементы
Салганик Е.А., Шхинек К.Н. Вибрация шельфовых сооружений при действии льда 75
Т, МПа
1 í т шах
/
/
/ !
! \ - 1 А 2 - -
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Рисунок 4. Трилинейный закон тяги I и деформации ис связывающих элементов
В работе для моделирования использовался материал с произвольной связью тяги и деформации, которая была принята трилинейной (табл. 2).
Таблица 2. Параметры связывающих элементов
Параметр Вертикальные связывающие элементы Горизонтальные связывающие элементы
Массовая плотность р кг/м3 910 910
Деформация до предельной тяги А1 0,08 0,10
Деформация до ослабления материала Л2 0,45 0,55
Интенсивность нормального выделения энергии Н/м 2-5200 2-5200
Интенсивность тангенциального выделения энергии ¿и, Н/м 2-5200 2-5200
Предельная тяга Ттах, МПа 0,3-1,0 1,1
Статическое давление на сооружение от ледовой нагрузки р находят через значения напряжений на поверхностях конечных элементов на границе между сооружением и льдом через интегрирование локальных значений давления по поверхности контакта:
р = |ргйЛ . (2)
Глобальная ледовая нагрузка Р определяется по формуле:
^ = рЛ = рБИ, (3)
где й - ширина сооружения, И - толщина льда.
Реакция сооружения определяется по формуле:
^.^Г = кХ.1Г , (4)
где - коэффициент упругости пружины, определяющей податливость сооружения, -перемещение сооружения относительно центра своего равновесия.
Прочность на сжатие
При взаимодействии с сооружением лед подвергается значительным сжимающим нагрузкам, поэтому необходимо проверить реалистичность поведения модели на примере хорошо исследованных испытаний на одноосное сжатие [26]. Значения прочности льда на одноосное неограниченное сжатие находятся в пределах от 0,5 МПа до 12 МПа [1].
} -0.06- /
Рисунок 5. Схема определения прочности модели льда на одноосное сжатие
Расчетные образцы представляли собой призмы с квадратным сечением. Задавалась нагрузка Р©, линейно возрастающая со временем. Регистрировались нагрузка при разрушении
Ртах и площадь поперечного сечения А . Прочность льда на одноосное сжатие С равна:
P
q _ -1 max
" A
(5)
Была определена зависимость прочности модели на сжатие от прочностных параметров отдельных материалов: предела текучести <гу объемных элементов, предела прочности на
растяжение Т и сдвиг Б связывающих элементов, интенсивности выделения энергии при трещинообразовании С| и Си.
Прочность на сжатие С модели льда значительно зависит от формы и параметров кривой зависимости между напряжением и деформацией.
С, МП а
2.5
1.5
1.0
Ж
^"модель
• • аналитическое решение г—►
1.5
2.0
2.5
3.0
Рисунок 6. Зависимость прочности на сжатие С от предела текучести <х
для модели льда и для аналитического решения С(< ) = <
y
В интервале значений предела текучести 1,5-3,0 МПа разность между пределом текучести объемных элементов ау и пределом прочности на сжатие образца С не превышает 12 % (рис. 6).
С, МП а
3.5
3.0 2.5 2.0 1.5 1.0 0.5
А
\\ \\
—-Т = 1.С т П 1 МПа МПа гическое ре] -
• • аналгп -- пение - -►
► С, Н/м
1
10
100
1 ООО
10 000
Рисунок 7. Зависимость прочности на сжатие C от интенсивности выделения энергии при трещинообразовании G при различных пределах прочности на растяжение I для модели льда и для аналитического решения
Прочность льда на сжатие С зависит от параметров связывающих элементов. При изменении интенсивности выделения энергии при трещинообразовании в нормальном и касательном направлениях С и Сц в пределах 10-5200 Н/м прочность льда на сжатие изменяется в диапазоне 95-150 % от предела текучести объемных элементов ау (рис. 7).
Наилучшее соответствие с аналитическим решением модель льда показала при пределе текучести объемных элементов в диапазоне 1,5-2,5 МПа и интенсивности выделения энергии при трещинообразовании 5200 Н/м.
Прочность модели льда на сжатие:
прямо пропорционально увеличению предела текучести
объемных
• повышается элементов;
• повышается при понижении интенсивности выделения энергии трещинообразования;
• незначительно понижается при уменьшении предела прочности на растяжение и сдвиг вертикальных связывающих элементов.
Изгибная прочность
Экспериментальные значения изгибной прочности а у находятся в пределах 0,1-1,0 МПа. Изгибная прочность ау определяется прочностью льда на растяжение [27].
Для проведения тестов была создана модель, представляющая собой балку прямоугольного сечения на двух роликовых опорах (рис. 8). Разрушению отвечает момент удаления первого связывающего элемента модели. Нагрузка приложена к центру верхней поверхности балки вдоль оси 2.
Рисунок 8. Схема определения прочности модели льда на изгиб
Неизменными при моделировании изгиба балки остаются параметры горизонтальных связывающих элементов: предельная тяга в нормальном и касательном направлении
тког = тког = х,шш.
(Тр МП а
6.0 5.0 4.0 3.0 2.0 1.0 0.0
А
— Г 1.0 -■■Т = 0.3 МПа
МПа • 1 - ■
.у ♦ ж ♦ ж ♦ Ж
♦ / ♦ Ж ♦ Ж ♦ Ж
........ I—►
G, Н/м
1
10
100
1 000
10 000
Рисунок 9. Зависимость прочности на изгиб < от интенсивности выделения энергии
при трещинообразовании О при различных пределах прочности на растяжение Т
Предел прочности на изгиб увеличивается при возрастании интенсивности выделения энергии при трещинообразовании С (рис. 9) и находится в пределах 1,5-6,0 от величины прочности на растяжение Т вертикальных связывающих элементов.
Прочность на растяжение
Прочность льда на одноосное растяжение находится в пределах 0,1-2,0 МПа [1]. Для рассматриваемой модели льда на прочность образца при растяжении влияют параметры
вертикальных связывающих элементов: предельная тяга в нормальном направлении Тх и интенсивность выделения энергии в нормальном направлении О1.
При незначительных величинах интенсивности выделения энергии предел прочности на растяжение связывающих элементов тх совпадает с полученным пределом прочности на
растяжение модели льда т . При значительной интенсивности выделения энергии часть растягивающих усилий поглощается в связывающих элементах, совершая работу по относительному перемещению их противоположных плоскостей.
Горизонтальная нагрузка приложена ко всем узлам крайних граней балки, перпендикулярных направлению нагрузки, вдоль оси X (рис. 10).
Рисунок 10. Схема определения прочности модели льда на растяжение
В зависимости от значений О прочность на растяжение модели Т находится в пределах 1,0-1,7 от предельной тяги связывающих элементов Тх = \,0МПа и в пределах 0,9-5,2 от
предельной тяги связывающих элементов Тх = 0,3МПа.
По результатам моделирования свойств льда можно сделать следующие заключения.
1. Определено влияние основных параметров модели на прочностные характеристики модельного льда.
2. Прочностные параметры модельного льда, в общем, близки к параметрам реального льда. Однако следует учитывать некоторое влияние масштаба, особенно при растяжении. Следует отметить, что коэффициент интенсивности К1с является размерным, что приводит к масштабному эффекту в механике разрушения.
3. В зависимости от свойств льда по данным испытаний на сжатие, растяжение и изгиб можно откалибровать входные параметры для получения необходимых прочностных свойств модели.
Глобальная нагрузка на жесткие сооружения
Для уверенности в результатах расчетов вибрации предварительно рассмотрена более простая задача определения квазистатических ледовых нагрузок на сооружения цилиндрического и прямоугольного поперечного сечений и проведено сопоставление их с результатами, определенными по нормам.
Рисунок 11. Конечно-элементная схема взаимодействия льда и жесткого или податливого сооружения прямоугольного или цилиндрического сечения
При моделировании рассматриваются прямоугольное и цилиндрическое сооружения, жестко закрепленные в узлах поверхности, не контактирующей со льдом.
Модель льда соответствует модели, описанной в предыдущих разделах. Часть входных параметров модели льда была постоянной во время всех проведенных тестов: плотность льда
р = 910кг / м ; модуль Юнга Е = 5ГПа; коэффициент Пуассона у = 0,3 ; предел на растяжение
и сдвиг вертикальных и горизонтальных связывающих элементов Т^°г = Т^°г = 1,1МПа.
Для нахождения оптимальных параметров в процессе проведения экспериментов меняются следующие величины: скорость ледового поля V = 0,1-1,0 м/с; интенсивность выделения энергии при трещинообразован
ии СЗIII = 100-5200 Н/м; толщина льда h = 0,3-0,5 м. Поперечный размер сооружения й = 0,5-0,7 м. При различных входных параметрах модели находятся значения давления от ледовой нагрузки ртах и глобальной ледовой нагрузки Ртах.
Значения максимального статического давления ртах прямо пропорциональны значениям
предела текучести объемных элементов льда <Уу При увеличении значений интенсивности
выделения энергии при трещинообразован
ии СС||| = 100-5200 Н/м давление увеличивается
(рис. 12).
МПа
2.5
2.0
1.5
1.0
0.5
0.0
i
прямоугси 1ьное у/^f
- ™цилиндри чс с кое ^j« / / ^ / /
/ ✓ ✓
/
- -
G, H/M
10
100
1 000
10 000
Рисунок 12. Зависимость максимального статического давления р шах на жесткое
сооружение цилиндрического и прямоугольного сечения от интенсивности выделения энергии при трещинообразовании О при скорости льда V = 0,5 м/с
Максимальное давление ршах прямо пропорционально значениям предела прочности на
сдвиг вертикальных связывающих элементов Т^,еП. Давление на жесткие сооружения достигает своего максимума в момент появления первых трещин, дальнейшее дробление льда не вызывает значительных нагрузок.
Глобальная нагрузка на податливые сооружения
Податливое сооружение при взаимодействии с ледовым полем можно описать двумя основными характеристиками: коэффициентом упругости к^ и массой М . На процесс взаимодействия также влияет скорость ледового поля V.
Частота собственных колебаний сооружения ¥«г может быть определена по формуле:
F -_
Г str -
1
2п\
к
str
M
(6)
str
Для более точного понимания влияния податливости на статическую нагрузку в данном разделе рассматривается взаимодействие льда и сооружения с параметрами, аналогичными использованным при моделировании взаимодействия с абсолютно жестким сооружением: толщина льда И = 0,3 м; диаметр цилиндра й = 0,5 м; ширина сооружения прямоугольного сечения й = 0,5 м.
При моделировании постоянными оставались: плотность льда р = 910кг/м ; модуль Юнга Е = 5ГПа; коэффициент Пуассона у = 0,3; предел прочности вертикальных и
горизонтальных связывающих элементов
Т™п = Т™ = 1,0МПа и Т™г = Т™г = 1,1МПа; предел текучести объемных элементов льда а у = 2,0МПа. Также не изменялись параметры кривой связи пластических напряжений и деформаций объемных элементов льда.
Варьировались скорость ледового поля V = 0,05; 0,1; 0,2; 0,5; 1,0 м/с и частота собственных колебаний сооружения = 1,5; 7,5; 10 Гц за счет изменения массы сооружения М
коэффициента упругости сооружения к
str
str
Для цилиндрического сооружения частоты собственных колебаний изменялись за счет изменения массы М3(г в диапазоне 10-1000 т при постоянном коэффициенте упругости
к3(г = 40МН / м и за счет изменения коэффициента упругости к^г в диапазоне 10-40 МН/м при
постоянной массе сооружения М = 10 т.
По итогам моделирования были определены условия возникновения вибраций при различных жесткостных характеристиках сооружений (рис. 13 и 14). При изменении жесткости сооружения кя(г параметры условий возникновения вибраций изменяются менее значительно,
чем при изменении массы сооружения М^ .
Рисунок 13. График зависимости отношения реакции сооружения и ледовой нагрузки Rstr / р от скорости льда V при собственной частоте сооружения Fstr = 1; 5; 7,5; 10 Гц и
постоянной жесткости 10 МН/м
Увеличению массы сооружения соответствуют большие скорости диапазона возникновения вибраций. Для сооружения с собственной частотой колебаний Fstr = 5 Гц при увеличении массы в 4 раза с 10 т до 40 т отношение реакции сооружения Rstr и ледовой нагрузки Fice увеличилось в 1,4 раза.
Диапазон максимальных вибраций соответствует диапазону скоростей ледового поля 0,10,5 м/с. При возникновении вибраций максимальные значения реакции сооружения Rstr находятся в пределах 2,0-3,4 от соответствующих значений ледовой нагрузки Fice.
Рисунок 14. График зависимости отношения реакции сооружения и ледовой нагрузки Я3{г / р от скорости льда V при собственной частоте сооружения Р5(г = 2,5; 5; 7,5 Гц и
постоянной массе 10 т
Для сооружения прямоугольного сечения с шириной й = 0,5 м были проведены испытания при постоянной частоте собственных колебаний Р^,- при различных значениях коэффициента упругости кя(г и массы сооружения М^ .
Уровень нагрузок и реакций на сооружения прямоугольного сечения с аналогичной диаметру шириной й = 0,5 м превышает уровень нагрузок и реакций на цилиндрическое сооружение. Значительное воздействие на уровень вибраций сооружений прямоугольного сечения оказывает масса сооружения. При увеличении массы в 10 раз величина динамического коэффициента возрастает в 2 раза (рис. 15). Однако при дальнейшем увеличении массы при сохранении частоты собственных колебаний значения динамического коэффициента находятся в пределах 2,5-3,5.
Рисунок 15. График зависимости отношения реакции сооружения и ледовой нагрузки Rstr / р от скорости льда V при собственной частоте сооружения прямого сечения Fstr = 5 Гц в зависимости от массы сооружения
Аналитическое решение подобной задачи с помощью рекомендаций стандарта ISO 19906 не позволяет учитывать влияние массы на уровень вибраций, т. к. масса вместе с жесткостью сооружения определяет постоянную частоту его собственных колебаний.
Для аналитического решения задачи мы задаем сооружение материальной точкой с одной степенью свободы с соответствующей массой Mstr , на линейной пружине с жесткостью kstr , на
которую действует нагрузка, зависящая от времени, с максимальным значением Fmax и
амплитудой Ар = pmax - pmin = 9pmax ■
Коэффициент q принимается равным 0,3 и уточняется при нахождении максимальной
скорости сооружения при нагружении V^ = 0,084 • Pstr ■ Полученные для сооружений с
собственной частотой колебаний Fstr = 1; 5; 10 Гц нагрузки находятся в диапазоне 1,2-1,7 МН, что в 3,5-5,0 раз превышает нагрузки, полученные с помощью моделирования. Динамический коэффициент при аналитическом решении задачи находится в пределах значений 2,6-3,6.
Динамическая нагрузка на маяк «Норстрёмсгрюнд»
Маяк «Норстрёмсгрюнд» находится в северной части Ботнического залива на территории Швеции. Его высота составляет 42 метра, диаметр на уровне ватерлинии - 7,5 м. Тип сооружения - цилиндрический кессон.
Рисунок 16. Конечно-элементная схема взаимодействия маяка «Норстрёмсгрюнд» и
ледового поля
Во время исследовательского проекта по измерению нагрузок от тонкого льда (LOLEIF) [28] маяк «Норстрёмсгрюнд» был оборудован панелями для измерения ледовых нагрузок и устройствами по измерению толщины льда. Во время проекта STRICE были проведены измерения нагрузок на маяк в период между 2001 и 2003 годами. Частота снятия данных варьировалась в пределах 1-100 Гц для разных панелей.
Скорость льда в рассматриваемых 2001 и 2002 годах находилась в пределах 0-0,6 м/с, толщина ровного льда в пределах 0-0,7 м. В этот период наиболее вероятным было значение скорости льда 0,1 м/с. По данным за 2001 год максимальное эффективное давление на маяк находилось в пределах 0,4-1,6 МПа. Наибольшим значениям максимального давления (1,0-1,6 МПа) соответствовали толщина льда 0,2-0,3 м и скорость льда 0,15-0,17 м/с.
Полевые и модельные нагрузки необходимо сравнить с рекомендованными нагрузками на гидротехнические сооружения от действия ровного льда. Максимальная глобальная ледовая нагрузка Fmax и максимальное статическое давление pmax в соответствии со стандартом ISO равны:
= CRhn ( D J = 1,4 • 0,6
-0,38
V 0,6 у
-0,16
= 1,13МПа,
Ртах = РтахР^ = 1,13 • 7,5 • 0,6 = 5,09МН
Для определения зависимости уровня вибраций сооружения от параметров ледового поля было смоделировано несколько вариантов взаимодействия при скоростях и толщинах ледового поля, наблюдаемых в отчетах различных исследовательских проектов. Параметры модели маяка (табл. 3) включают в себя массу сооружения М3(г, жесткость сооружения , коэффициент демпфирования Л и диаметр й.
Таблица 3. Параметры модели маяка «Норстрёмсгрюнд»
Масса M str Жесткость kstr Коэффициент демпфирования Л Диаметр D
кг Н/м - м
7,13 ■ 106 8,33 ■ 108 0,03 7,5
Были определены величины нагрузок и реакций маяка при действии льда толщиной 0,3 м и 0,6 м при скоростях V = 0,05-1,00 м/с. Уровень вибраций определен с помощью отношения реакции и нагрузки при различных скоростях ледового поля (рис. 17). Максимальные динамические нагрузки наблюдаются при скорости 0,20 м/с, что соответствует полевым данным проекта STRICE.
Rsti 1 Fice
1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0
А
/ / /
■ df— ¥ * >
N
Ч
шш ■h = 0.6 м
—h = 0.3 м
-1-1- -►
V, м/с
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40
Рисунок 17. График зависимости отношения реакции маяка «Норстрёмсгрюнд» и ледовой
нагрузки R
str / Е„
от скорости льда V
Нагрузки были определены и аналитически с помощью рекомендаций стандарта ISO 19906. Уровень нагрузок (2,5-5,1 МН) и значения динамического коэффициента (2,6-5,0) при таком способе моделирования значительно превышают значения конечно-элементного моделирования и значения, полученные в ходе полевых экспериментов проектов LOLEIF и STRICE на маяке «Норстрёмсгрюнд».
При максимальном давлении на измерительные панели, установленные на маяке, в 1,6 МПа и среднем значении прочности льда на сжатие 1,1 МПа значения динамического коэффициента не превышают 1,5. Аналогичные значения отношения реакции и нагрузки получены с помощью моделирования в программном комплексе ANSYS.
Динамическая нагрузка на платформу «Моликпак»
Платформа «Моликпак» представляет собой сооружение кессонного типа, предназначенное для установки на подводную песчаную берму. 12 мая 1986 года на платформе возникли значительные вибрации, что привело к увеличению нагрузок на сооружение до 250 МН.
Таблица 4. Параметры максимальных нагрузок на платформу «Моликпак» [2]
Дата Скорость льда Толщина льда Нагрузка Прочность льда
- м/с м МН МПа
17 февраля 1986 0,03-0,06 0,5-0,7 60 1,7
7 марта 1986 0,05 2,8 300 1,3
12 апреля 1986 0,06 3,5 400 1,5
12 мая 1986 0,1 1,9-3,5 250 1,4
25 июня 1986 0,1 1,7-2,5 115 0,8
Частота собственных колебаний платформы Fstr составляет 1,8 Гц. По рекомендациям стандарта ISO 19906 вибрации на платформе должны возникать при скоростях, близких к значению 0,15 м/с. Коэффициент упругости и масса платформы Моликпак были определены через ее геометрические характеристики (табл. 5).
Таблица 5. Параметры модели платформы «Моликпак»
Масса Mstr Жесткость kstr Коэффициент демпфирования Л Диаметр D
кг Н/м - м
7,8107 1,0-1010 0,03 90,0
Рисунок 18. Конечно-элементная схема взаимодействия платформы «Моликпак» и ледового поля
Проведено моделирование взаимодействия платформы и ледового поля с толщиной h = 1 м и диапазоном скоростей V = 0,05-0,5 м/с (табл. 6). Максимальная реакция платформы соответствует скорости ледового поля V = 0,1 м/с при величине динамического коэффициента 2,2.
Таблица 6. Результаты моделирования взаимодействия льда и платформы «Моликпак»
h V т/max Vstr max xstr E tan Rstr / Flce Rstr Rstr / Fjce
м м/с м/с м МПа МПа МН МН -
1,0 0,05 0,03 28,2 1,4 - 173,4 264,8 1,53
1,0 0,1 0,05 42,0 1,4 - 173,6 382,0 2,20
1,0 0,2 0,09 22,6 1,4 - 173,0 225,0 1,30
1,0 0,3 0,28 35,2 1,4 - 176,0 281,6 1,61
1,0 0,5 0,27 35,5 1,4 - 180,6 282,4 1,56
Значения полученных нагрузок превышают значения полевых нагрузок из-за того, что в полевых условиях ледовое поле действовало не на всю платформу «Моликпак», а вектор скорости поля не был перпендикулярен поверхности сооружения. Диапазон максимальных динамических нагрузок при моделировании соответствует условиям вибрации при эксплуатации платформы.
Заключение
Целью работы являлось изучение процесса возникновения вибраций сооружения при действии ровного льда. Для этого решена задача по созданию конечно-элементной модели льда в программном комплексе ANSYS.
Для созданной модели льда проведены базовые тесты по определению прочностных характеристик материала: одноосное сжатие, растяжение и изгиб. Были уточнены и связаны с глобальными характеристиками льда входные параметры конечно-элементной модели. Также определены статические нагрузки на сооружения круглого и прямоугольного сечения, их значения сопоставлены с рекомендациями стандарта ISO 19906. Для моделей податливых сооружений при различных величинах массы и жесткости определены нагрузки и реакции при различной скорости и толщине ледового поля, найдены условия возникновения вибраций. Значения нагрузок сопоставлены со значениями, определенными аналитически согласно рекомендациям стандарта ISO 19906.
Апробация модели проведена на примере маяка «Норстрёмсгрюнд» и платформы «Моликпак», параметры взаимодействия сооружения и ровного льда соотнесены со значениями полевых измерений. Величины нагрузок и условия возникновения динамических нагрузок, полученных с помощью конечно-элементного моделирования, имеют приемлемые расхождения с полевыми данными.
Данный метод определения ледовых нагрузок от ровного льда на податливые и жесткие сооружения позволяет учитывать множественные параметры: прочность льда на сжатие, растяжение и изгиб, массу и жесткость сооружения, форму поверхности сооружения, скорость и толщину ледового поля, коэффициенты трения между льдом и материалом сооружения. Это обеспечивает его преимущества перед аналитическими и полуаналитическими моделями.
В данной модели использовался упрощенный способ учета плавучести льда, что обеспечило приемлемую скорость расчета. Для более детального расчета можно учесть влияние вязких свойств воды, для которой была создана конечно-элементная модель. Однако учет этого параметра требует использования значительных расчетных ресурсов. Также для большей точности расчета необходимо учесть зависимость прочностных свойств льда от скорости деформации, что может увеличить уровень возникающих вибраций.
Литература
1. Лосет, С., Шхинек, К.Н., Гудместад, О., Хойланд, К. Воздействие льда на морские и береговые сооружения. СПб.: Лань, 2010. 272 с.
2. Bjerkâs M. Ice action on offshore structures - with applications of continuous wavelet transforms on ice load signals: PhD Thesis. Norwegian University of Science and Technology, 2006. 176 p.
3. Jeffries M.G., Wright W.H. Dynamic response of "Molikpaq" to ice - structure interaction // Proceedings of the OMAE Conference. 1988. Pp. 201-220.
4. ISO/FDIS 19906: Petroleum natural gas industries - Arctic offshore structures, ISO TC 67/SC 7 N. International Organization for Standardization, 2010.
5. Gïirtner A. Experimental and Numerical Investigations of Ice-Structure Interaction: PhD Thesis. Norwegian University of Science and Technology, 2009. 182 p.
6. Peyton, H. Sea Ice Forces. Ice Pressure against Structures // National Research Council of Canada, Ottawa, Technical Memorandum No. 92. Pp. 117-123.
7. Karna T., Turunen R. A straightforward technique for analyzing structural response to dynamic ice action // OMAE'90. Proceedings of the 9th International Conference Offshore Mechanics and Arctic Engineering. 1990. Pp. 135-142.
8. Nandan H., Younan A. H., Deng L. Ice induced vibration implementation of Maattanen model and development of design supplements // Proceedings of the 21st International Conference on Port and Ocean Engineering under Arctic Conditions, Montreal, Canada, July 2011. Pp. 510-519.
9. Yue Q., Xianguing B., Sun B., Zhang T., Xiangyu C. Full scale ice measurements on JZ20-2 platform // Proceedings of the IAHR Ice Symposium. 1996. Vol. 1. Pp. 282-289.
10. Blenkarn K.A. Measurement and Analysis of the ICE Forces on Cook Inlet Structures // Proceedings of the 2nd Offshore Technology Conference, Houston, Texas. 1970. Vol. 2. Pp. 265-378.
11. Olsen A.I. Time-Variant Reliability of Dynamic Systems by Importance Sampling and Probabilistic Analysis of Ice Loads / Norwegian University of Science and Technology, 2006. 153 p.
12. GQrtner A., Bjerkas M., KQhnlein W., Jochmann P., Konuk I. Numerical simulation of ice action to a lighthouse: Proceedings of the ASME 28th International Conference on Ocean, Offshore and Arctic Engineering. 2009. Pp. 175-183. 2009.
13. Tvergaard V., Hutchinson J.W. The Relation between Crack Growth Resistance and Fracture Process Parameters in Elastic-Plastic Solids // Journal of Mechanics and Physics of Solids. 1992. Vol. 40. Pp. 1377-1397.
14. Barrenblatt G.I. The Mathematical Theory of Equilibrium of Cracks in Brittle Fracture // Advances of Applied Mechanics. 1962. Vol. 7. Pp. 55-129.
15. Daiyan H., Sand B. Numerical Simulation of the Ice-Structure Interaction in LS-DYNA // 8th European LS-DYNA Users Conference, Strasbourg, 2011.
16. Sand B. Nonlinear finite element simulations of ice forces on offshore structures: PhD Thesis, Lulea University of Technology, 2008. 241 p.
17. Zhilenkov A., Shkhinek K., Thomas G. The Vibration of Fixed Offshore Structures under Ice Action. // SPE Arctic and Extreme Environments Conference and Exhibition, Russia, 2011. Pp. 917-933.
18. Карулин Е.Б., Карулина М.М., Шхинек К.Н., Томас Г.А.Н. Математическое и физическое моделирование взаимодействия чертырёхопорной конструкции с ровным льдом // Труды ЦНИИ им. акад. А.Н. Крылова. 2011. №63. С. 81-92.
19. Шхинек К.Н. Ледовые нагрузки на арктические платформы // Гидротехническое строительство.1994. №3. С. 33-36.
20. Lemaitre J., Chaboche J.-L. Mechanics of solid materials. Cambridge university press, 1990. 584 p.
21. Морозов Е.М. ANSYS в руках инженера: Механика разрушения.М.: ЛЕНАНД, 2010. 456 c.
22. Kuehn G.A., Lee R.W., Nixon W.A., Schulson E.M. The structure and tensile behavior of first-2014 sea ice and laboratory-grown saline ice // Journal of Offshore Mechanics and Arctic Engineering. 1990. Vol. 112. Pp. 357-363.
23. Weimar K. LS-DYNA User's Guide / CAD-FEM GmbH, PWS, 2001. 219 p.
24. Hallquist J.O. LS-DYNA user's manual. Livermore software technology corporation, 2009. 3020 p.
25. Муйземнек А.Ю., Богач А.А. Математическое моделирование процессов удара и взрыва в программе Ls-Dyna. Пенза: ИИЦ ПГУ, 2005. 106 с.
26. Лавров В.В. О строительных нормах прочности льда на сжатие // Гидротехническое строительство. 1967. №2. С. 39-42.
27. Никитин В.А., Сухоруков К.К. Прочность морских ледяных полей и ледовые нагрузки // Метеорология и гидрология. 1998. №12. С. 88-95.
28. Jochmann P., Schwarz J. Ice Force Measurement at Lighthouse Norstromsgrund - Winter 1999 // Validation of Low Level Ice Forces on Coastal Structures: LOLEIF, Report No.:5, Hamburg, 2000. Pp. 749-753.
Евгений Александрович Салганик, Санкт-Петербург, Россия Тел. моб.: +7(952)2266134; эл. почта: salganikea@mail.ru
Карл Натанович Шхинек, Санкт-Петербург, Россия Тел. моб.: +7(921)9499683; эл. почта: karl-1303@mail.ru
© Салганик Е.А., Шхинек К.Н., 2014
doi: 10.5862/MCE.48.8
Ice induced vibrations of offshore structures
E.A. Salganik
CJSC GT Morstroy, Saint-Petersburg, Russia +79522266134; e-mail: salganikea@mail.ru
K.N. Shkhinek
Saint-Petersburg State Polytechnical University, Saint-Petersburg, Russia
+79219499683; e-mail: karl-1303@mail.ru
Key words
vibration; ice load; offshore structures; fracture mechanics; computer simulation; finite element method; ANSYS
Abstract
The aim of the article was to analyze the process of the dynamic interaction between ice and offshore structures and find conditions of vibration beginnings. The problems of creating a full-scale finite element model of ice and offshore structure interaction and comparison of the received results of vibration conditions with the field data were solved. Modeling was realized with the finite element analysis software ANSYS.
The problems of modeling cracks formation, checking the multilateral model using basic ice strength tests and finding ice loads level on different types of structures were solved in order to create a destructible ice model. The object of the investigation is the role of elastic and inertial offshore structure parameters and geometrical and dynamical ice field parameters in the process of ice-induced structure vibrations. Molikpaq drilling and oil production platform in Canada and Norstromsgrund lighthouse in Sweden were chosen as objects for modeling.
The research contains verification of the results, their comparison to field data, load magnitudes from ISO 19906 «Petroleum and natural gas industries - Arctic offshore structures» algorithm and results of discrete element modeling. The results are acceptable in comparison with the results of other ice-induced vibration research.
References
1. L0set S., Shkhinek K.N., Gudmestad O., H0yland K. Vozdeystviye lda na morskiye i beregovyye sooruzheniya [Actions from Ice on Arctic Offshore and Coastal Structures]. Saint-Petersburg: Lan, 2010. 272 p. (rus)
2. Bjerkas M. Ice action on offshore structures - with applications of continuous wavelet transforms on ice load signals: PhD Thesis, Norwegian University of Science and Technology, 2006. 176 p.
3. Jeffries M.G., Wright W.H. Dynamic response of "Molikpaq" to ice - structure interaction. Proceedings of the OMAE Conference. 1988. Pp. 201-220.
4. ISO/FDIS 19906: Petroleum natural gas industries - Arctic offshore structures, ISO TC 67/SC 7 N. International Organization for Standardisation, 2010.
5. Gurtner A. Experimental and Numerical Investigations of Ice-Structure Interaction: PhD Thesis. Norwegian University of Science and Technology, 2009. 182 p.
6. Peyton H.R. Sea Ice Forces. Ice Pressure against Structures. National Research Council of Canada, Ottawa, Technical Memorandum No. 92. Pp. 117-123.
7. Karna T., Turunen R. A straightforward technique for analyzing structural response to dynamic ice action. OMAE'90. Proceedings of the 9th International Conference Offshore Mechanics and Arctic Engineering. 1990. Pp. 135-142.
8. Nandan H., Younan A. H., Deng L. Ice induced vibration implementation of Maattanen model and development of design supplements. Proceedings of the 21st International Conference on Port and Ocean Engineering under Arctic Conditions, Montreal, Canada, July 2011. Pp. 510-519.
9. Yue Q., Xianguing B., Sun B., Zhang T., Xiangyu C. Full scale ice measurements on JZ20-2 platform. Proceedings of the IAHR Ice Symposium. 1996. Vol. 1. Pp. 282-289.
10. Blenkarn K.A. Measurement and Analysis of the ICE Forces on Cook Inlet Structures. Proceedings of the 2nd Offshore Technology Conference, Houston, Texas. 1970. Vol. 2. Pp. 265-378.
11. Olsen A.I. Time-Variant Reliability of Dynamic Systems by Importance Sampling and Probabilistic Analysis of Ice Loads / Norwegian University of Science and Technology, 2006. 153 p.
12. Gurtner A., Bjerkâs M., Kuhnlein W., Jochmann P., Konuk I. Numerical simulation of ice action to a lighthouse. Proceedings of the ASME 2009 28th International Conference on Ocean, Offshore and Arctic Engineering, 2009. Pp. 175-183.
13. Tvergaard V., Hutchinson J.W. The Relation between Crack Growth Resistance and Fracture Process Parameters in Elastic-Plastic Solids. Journal of Mechanics and Physics of Solids. 1992. Vol. 40. Pp. 1377-1397.
14. Barrenblatt G.I. The Mathematical Theory of Equilibrium of Cracks in Brittle Fracture. Advances of Applied Mechanics. 1962. Vol. 7. Pp. 55-129.
15. Daiyan H., Sand B. Numerical Simulation of the Ice-Structure Interaction in LS-DYNA. 8th European LS-DYNA Users Conference, Strasbourg, 2011.
16. Sand B. Nonlinear finite element simulations of ice forces on offshore structures: PhD Thesis, Luleâ University of Thechnology, 2008. 241 p.
17. Zhilenkov A., Shkhinek K., Thomas G. The Vibration of Fixed Offshore Structures under Ice Action. SPE Arctic and Extreme Environments Conference and Exhibition, Russia, 2011. Pp. 917-933.
18. Karulin E.B., Karulina M.M., Shkhinek K.N., Tomas G.A.N. Matematicheskoye i fizicheskoye modelirovaniye vzaimodeystviya chertyrekhopornoy konstruktsii s rovnym ldom [Mathematical and physical modeling of four foot construction and level ice interaction]. Trudy TsNII im. akad. A.N. Krylova [Proceedings of the Krylov State Research Centre]. 2011. No. 63. Pp. 81-92. (rus)
19. Shkhinek K.N. Ledovyye nagruzki na arkticheskiye platformy [Ice loads on arctic platforms]. Power Technology and Engineering. 1994. No. 3. Pp. 33-36. (rus)
20. Lemaitre J., Chaboche J.-L. Mechanics of solid materials. Cambridge university press, 1990. 584 p.
21. Morozov E.M. ANSYS v rukakh inzhenera: Mekhanika razrusheniya [ANSYS in the hands of the engineer: Fracture Mechanics]. Moscow: LENAND. 2010. 456 p. (rus)
22. Kuehn G.A., Lee R.W., Nixon W.A., Schulson E.M. The structure and tensile behavior of first-year sea ice and laboratory-grown saline ice. Journal of Offshore Mechanics and Arctic Engineering. 1990. Vol. 112. Pp. 357-363.
23. Weimar K. LS-DYNA User's Guide / CAD-FEM GmbH, PWS, 2001. 219 p.
24. Hallquist J.O. LS-DYNA user's manual. Livermore software technology corporation, 2009. 3020 p.
25. Muyzemnek A.U., Bogach A.A. Matematicheskoye modelirovaniye protsessov udara i vzryva v programme Ls-Dyna [Mathematical modeling of shock and explosion in the program LS-DYNA]. Penza: IITs PGU, 2005. 106 p. (rus)
26. Lavrov V.V. O stroitelnykh normakh prochnosti lda na szhatiye [Building codes for ice compressive strength]. Power Technology and Engineering. 1967. No. 2. Pp. 39-42. (rus)
27. Nikitin V.A., Sukhorukov K.K. Prochnost morskikh ledyanykh poley i ledovyye nagruzki [Sea ice floe strength and ice loads]. Russian Meteorology and Hydrology. 1998. No. 12. Pp. 88-95. (rus)
28. Jochmann P., Schwarz J. Ice Force Measurement at Lighthouse Norstromsgrund - Winter 1999. Validation of Low Level Ice Forces on Coastal Structures, LOLEIF, Report No. 5, Hamburg, 2000. Pp. 749-753.
Full text of this article in Russian: pp. 72-88