Научная статья на тему 'Обзор применения численных методов для моделирования деформации и разрушения льда'

Обзор применения численных методов для моделирования деформации и разрушения льда Текст научной статьи по специальности «Физика»

CC BY
508
138
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Труды ВИАМ
ВАК
RSCI
Область наук
Ключевые слова
ЛЕД / РАЗРУШЕНИЕ / МОДЕЛИРОВАНИЕ / МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ / ГИДРОДИНАМИКА СГЛАЖЕННЫХ ЧАСТИЦ / МЕТОД ДИСКРЕТНЫХ ЭЛЕМЕНТОВ / ICE / FRACTURE / SIMULATION / FINITE ELEMENT METHOD / SMOOTHED PARTICLE HYDRODYNAMICS / DISCRETE ELEMENT METHOD

Аннотация научной статьи по физике, автор научной работы — Гриневич Д. В., Бузник В. М., Нужный Г. А.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Гриневич Д. В., Бузник В. М., Нужный Г. А.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

REVIEW OF NUMERICAL METHODS FOR SIMULATION OF THE ICE DEFORMATION AND FRACTURE

Рresents a review of various numerical methods use for modeling of the ice deformation and brittle fracture. The application of the finite element method, smoothed particle hydrodynamics, and the discrete element method is considered. The paper deals with modeling of sea and freshwater ice, impact, explosive loading, three-point and four-point bending, compression, contact interaction «ship-ice» and «support-ice», falling from a height, and the use of reinforcement etc.

Текст научной работы на тему «Обзор применения численных методов для моделирования деформации и разрушения льда»

УДК 52-17:539.42

Д.В. Гриневич1, В.М. Бузник1, Г.А. Нужный1

ОБЗОР ПРИМЕНЕНИЯ ЧИСЛЕННЫХ МЕТОДОВ

ДЛЯ МОДЕЛИРОВАНИЯ ДЕФОРМАЦИИ И РАЗРУШЕНИЯ ЛЬДА

DOI: 10.18577/2307-6046-2020-0-8-109-122

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

Ключевые слова: лед, разрушение, моделирование, метод конечных элементов, гидродинамика сглаженных частиц, метод дискретных элементов.

D. V. Grinevich1, V.M. Buznik1, G.A. Nuzhnyi1

REVIEW OF NUMERICAL METHODS FOR SIMULATION OF THE ICE DEFORMATION AND FRACTURE

Presents a review of various numerical methods use for modeling of the ice deformation and brittle fracture. The application of the finite element method, smoothed particle hydrodynamics, and the discrete element method is considered. The paper deals with modeling of sea and freshwater ice, impact, explosive loading, three-point and four-point bending, compression, contact interaction «ship-ice» and «support-ice», falling from a height, and the use of reinforcement etc.

Keywords: ice, fracture, simulation, finite element method, smoothed particle hydrodynamics, discrete element method.

"'Федеральное государственное унитарное предприятие «Всероссийский научно-исследовательский институт авиационных материалов» Государственный научный центр Российской Федерации [Federal State Unitary Enterprise «All-Russian Scientific Research Institute of Aviation Materials» State Research Center of the Russian Federation]; e-mail: admin@viam.ru

Введение

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

Моделирование необходимо для изучения сложных природных процессов, оно заключается в выделении основных особенностей и закономерностей, а также последующей идеализации изучаемого явления. Такое упрощенное представление об исследуемом объекте позволяет получить его описание с дальнейшим прогнозом поведения -например, деформации материала и разрушения. Адекватность моделей должна подкрепляться экспериментальными данными, а важным средством эффективного нахождения численного решения задач являются современные вычислительные методы и их реализация в программном обеспечении. Идеализация - сложный этап процесса изучения объекта, особенно материала сложной структуры - например, слоистой [1-3] или ячеистой [4]. Лед и композиционные материалы на его основе применяют в качестве конструкционного материала для построения инфраструктуры в арктических зонах

и в регионах с холодным климатом: дорог, взлетно-посадочных полос, площадок разгрузки и погрузки и др. [5-7]. При использовании льда желательно, чтобы его поведение было прогнозируемым. Хотя большое внимание уделяется его экспериментальному исследованию [8-11], для дальнейшего развития работ необходимо проводить моделирование льда, дающее возможность объяснять и предсказывать поведение конструкций в различных условиях эксплуатации.

Получение численных моделей льда - нетривиальная задача с учетом его варьируемых прочностных свойств, хрупкости, сильной чувствительности свойств льда к заморозке и условиям эксплуатации [11]. Лед - особый материал и на него нельзя напрямую распространить методы моделирования для остальных классов материалов. В композиционных материалах на основе льда изучение свойств усложняется вследствие взаимодействия разных по структуре фаз. Применение традиционных методов часто затруднено и требует предварительных исследований.

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

Метод конечных элементов

Метод конечных элементов (МКЭ) широко распространен для решения задач прикладной механики конструкций. Лед можно считать упругим только при очень малых деформациях, поэтому в этом случае допустимо использовать простейший вариант -модель линейно-упругого тела [12].

Моделирование более сложного процесса нелинейного деформирования и хрупкого разрушения требует использования расширенных моделей [13, 14], к которым можно отнести нелинейные модели МКЭ изотропных материалов Джонсона-Холмквиста, Леметра и др.

В работе [15] рассмотрены и проанализированы с помощью МКЭ различные модели материалов и проведено моделирование в программном комплексе LS-Dyna. Экспериментальные данные взяты из открытых научно-технических литературных источников и собственных данных автора по пресноводному поликристаллическому льду. Кроме применения различных критериев разрушения, автором для получения сходящегося решения применяется метод удаления элементов из модели при достижении предельного шага интегрирования.

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

В результате моделирования динамики контактного взаимодействия твердого тела с моделью льда получено, что наиболее согласующиеся с экспериментальными данными результаты демонстрируют модели упругопластического материала с кинематическим упрочнением и вязкопластический материал. В работе [15] для моделирования взята упругопластическая модель материала Прандтля-Рейса с изотропным упрочнением и критерием текучести Мизеса:

[а, - от=0,

где еР - компоненты приращения пластических деформаций; 'ksjj - компоненты девиатора напряжения; ог- - интенсивность напряжений; от - предел текучести.

В зоне ползучести эта модель может применяться без критерия разрушения. Для зоны разрушения более применима гидродинамическая модель поведения разрушенных элементов.

Рассмотренные модели являются наиболее универсальными, но требуют разделения их применения по исследуемым процессам деформации льда.

Применение изотропного упругопластического материала с изотропным упрочнением приводится в работах [16, 17] при моделировании ударного взаимодействия упругого цилиндра со льдом. Расчеты проводятся в программном комплексе LS-Dyna, а также приводятся результаты для авторской программы для двухмерной постановки задачи. Указывается и основной недостаток проведения подобного рода моделирования -точность определения свойств льда из экспериментальных данных. Так потребовалась корректировка методики обработки экспериментальных данных из-за завышенных значений кривой зависимости предела текучести при сжатии от скорости деформации. Получена удовлетворительная корреляция экспериментальных и расчетных данных в широком диапазоне скоростей удара (118-550 м/с). В настоящее время в работах [16, 17] указывается недостаток данных по деформированию разрушаемого льда, что необходимо для определения и расширения границ применимости физико-механических моделей динамического поведения материала.

В работах [ 18-20] подробно рассмотрено моделирование взрывного воздействия на лед. Коллективом авторов разработано собственное программное обеспечение для расчета такого воздействия на лед в двухмерной постановке задачи и используется модель упругопластического течения Прандтля-Рейса с критерием текучести Мизеса. Получена хорошая сходимость результатов расчетов с большим количеством собственных экспериментальных данных, накопленных авторами, а также взятых ими из научно-технических источников. Специфика работ авторов - взрывное воздействие на лед для расчистки русла рек.

Данные работы посвящены скоростному нагружению льда. Для изучения поведения ледяного покрова, ледяных переправ, льдин, припойного льда необходимо исследование механики квазистатического нагружения менее интенсивного характера. В отечественных работах [21-23] представлен опыт моделирования ледяной балки при четырехточечном изгибе. Лед описывается упругопластическим материалом с нелинейным критерием текучести Виллама-Варнке, разработанного для моделирования бетона. Выбор параметров для определения поверхности текучести по данному критерию представлен в работе [24]. Исследовано армирование, выполненное со стальным каркасом. Нагружение проведено только до начала разрушения, т. е. до активного растрескивания ледяной матрицы и наступления условной «пластичности» нагружаемой балки. Проведено моделирование и продемонстрирована хорошая корреляция с результатами эксперимента. На основе полученных моделей делается попытка вариативного исследования других материалов для армирования. Прочие типы нагружения не принимались во внимание. Данное исследование отличается тем, что рассматривается армированный лед, применение которого позволяет существенно повысить прочностные и деформационные свойства материала.

В зарубежной работе [25] также исследуется разрушение ледяной балки при четырехточечном изгибе. Проводится оптимизация модели материала с раздельным поведением при растяжении и сжатии, с критерием текучести по Мизесу. В качестве критерия разрушения выбран критерий для разрушения композиционных материалов Цая-Ву. При нагружении четырехточечным изгибом получена хорошая корреляция с экспериментальными данными, но применение модели потребовало проведения оптимизационного поиска параметров модели материала.

В работах [26, 27] представлено необычное применение метода когезии, который заключается в моделировании между конечными элементами промежуточных виртуальных интерфейсных элементов. При достижении предельных условий интерфейсные элементы разъединяются. Основная особенность метода заключается именно в заранее определенной зоне с интерфейсными элементами, т. е. возможное направление роста трещин определяется заранее. Традиционная область применения метода -механика композиционных материалов, граница сред, моделирование клеевых соединений, учет адгезии арматуры в железобетоне. Однако авторы адаптировали данный метод для хрупкого разрушения льда, для чего интерфейсные элементы были добавлены между всеми конечными элементами модели. Представлены результаты трехмерного моделирования для следующих типов нагружения: растяжение, сдвиг, растяжение при раскалывании (brazilian test), падение с высоты на упругую и жесткую поверхность. В целом получена удовлетворительная корреляция экспериментальных и расчетных данных, особенно в процессах с преобладанием одной трещины при растрескивании объемного тела. Но существует проблема контакта интерфейсных элементов и внешних тел - например, при моделировании падения образца на жесткую поверхность получены на порядок большие контактные силы. С увеличением количества элементов увеличилось количество степеней свободы, соответственно, значительно возросли требования к вычислительным ресурсам. Авторами работ [26, 27] получена зависимость результата от параметров разбиения объемного тела. Параметры контакта должны быть подобраны каждый раз в зависимости от сетки модели. Данные исследования заявляются первой попыткой в применении данного метода в механике хрупкого разрушения. Постановка задачи напоминает переход уже к бессеточным методам.

Применение МКЭ в моделировании хрупкого разрушения льда обобщено в работе [28]. В отличие от ранее указанных работ исследуется и пресноводный, и морской лед - менее прочный и более пластичный. Исследование проводилось на разномасштабных моделях, при этом демонстрируется неуниверсальность существующих моделей материалов.

Автор работы [28] использует разные типы моделей материалов МКЭ для хрупкого характера разрушения, в основном разработанных для таких материалов, как бетон и гранит (применен программный комплекс LS-Dyna). На основе анализа ряда работ приводится выбор подходящих для хрупкого разрушения моделей материалов. Проведено моделирование с применением следующих моделей материалов [29]:

- разрушаемый пеноматериал (MAT_063, MAT_CRUSHABLE_FOAM в LS-Dyna);

- модель материала с повреждаемостью Леметра (MAT_153, MAT_DAMAGE_3);

- модель материала с хрупким разрушением Говинджи-Кея-Симо (MAT_096, MAT_BRITTLE_DAMAGE);

- модель материала Джонсона-Холмквиста (MAT_111, AT_JOHNSON_ HOLMQUIST_CONCRETE).

Применение моделей пеноматериалов применительно ко льду рассматривалось также в работах [30-32]. Данная модель изотропная, одномерная. При растяжении происходит разрушение. Модуль упругости постоянен, напряжения при нагрузке для

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

Расчет напряжении после шага нагружения (п) с промежуточной оценкоИ (п+1/2) производится согласно формуле

п+1/2 п , 77* п+1/2 л ,,.+1/2 /ОЧ

°и =°и+ Ее п ^ ' (2)

где Оу - тензор напряжений; Е - модуль упругости; 8у - тензор скоростей деформации; А( - расчетный шаг времени.

О у <

Если значения главных напряжений выходят за поверхность текучести _п+1/2 то они нормируются до уровня поверхности:

п+1/2

п+1 О О, = О

'1 "у

п+1/2

О

Модель материала с хрупким разрушением Говинджи-Кея-Симо и модель бетона Джонсона-Холмквиста применялись ранее только для бетона и гранита, для льда они применены и адаптированы в работе [28].

Модель с хрупким разрушением Говинджи-Кея-Симо предполагает, что разрушение (накопление повреждений) в материале начинается, когда главные напряжения превышают некоторое пороговое значение (подробное описание дано в работе [29]). Повреждаемость учитывается путем уменьшения модуля упругости материала по экспоненциальному закону. Преимуществом данной модели является минимальный набор параметров, описывающих материал: массовые, жесткостные, вязкостные параметры, а также предельные напряжения.

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

Текущая прочность материала на расчетном шаге времени (п+1) с учетом повреждаемости равна (здесь и далее - знак * обозначает нормированный параметр)

* п+1 * п+1/2 гч/ * п+1/2 * п+1/2 \ (1\

о = о* - В {с,- -О^ ), (3)

где о* = Л(р* +Т) п [1+С-1п(в *р1)] - прочность неповрежденного материала; о*. = Б(р* )т [1+С-1п(в *р1)] - прочность для полностью разрушенного материала; г*р1 - скорость пластической деформации; р* - гидростатическое давление; Т, Б, С, т - константы материала.

Модель аккумулирует повреждения как от эквивалентной пластической деформации, так и от пластической объемной деформации:

п =-г-гг*ы> (4)

где 01, 02 - константы, описывающие характер разрушения конкретного материала.

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

Численные эксперименты по исследованию льда совместно с водой на основе модели повреждаемости Леметра можно найти в работе [34]. Совместные усилия авторов работ [25-27] и [34] представлены в международном проекте Евросоюза SAFEWIN для анализа судоходного движения судов во льдах Балтики. Как видно в публикации по проекту [35], авторами была выбрана и применялась именно модель Леметра с накоплением повреждаемости для анализа давления ледяного покрова на корпус судна. Проведенные испытания в специальном бассейне демонстрируют хорошую корреляцию моделирования с экспериментальными данными при данном типе разрушения.

Модель Леметра строится на концепции зависимости повреждаемости Б от величины пластической деформации. Влияние поврежденности среды на прочность выражается уравнением

* и+1/2

* и+1 _

Изменение меры поврежденности определяется соотношением

где ёр - эквивалентная скорость пластической деформации; и I - константы, характеризующие материал; У - плотность энергии деформации.

Полное описание модели приведено в работе [20].

В работе [28] рассмотрены две расчетные схемы - ударное нагружение бойком и контакт объемной пластины с твердым телом - аналог контакта «льдина-судно» или «льдина-опора моста» [36]. В отличие от работ [30-32] исследование в работе [28] показало, что модель пеноматериала с разрушением дает неудовлетворительный результат с завышенными на порядок напряжениями прежде всего для ударного нагружения, что можно объяснить более подробной моделью в работах [30-32], в которой использовалось разбивание ледового объекта на отдельные блоки с корректировкой характеристик для каждого блока.

Результаты для моделей, учитывающих повреждаемость, имеют желаемый физичный характер разрушения: кривая зависимости силы от перемещений для всех моделей материалов имеет пилообразный вид, адекватно описывает скачкообразно меняющийся характер разрушения льда. Порядок результирующих величин (реакций и перемещений) близок к экспериментальным данным. Однако следует отметить неустойчивость решения к изменению входных параметров моделей (варьирование приводит к несходящемуся решению) - определение параметров требует существенных сопутствующих материальных затрат на валидацию моделей, особенно для модели Джонсона-Холмквиста. Хрупкое разрушение с помощью приведенных моделей получить проще, чем пластичное деформирование льда. Моделирование потребовало значительных вычислительных ресурсов - проблемы со сходимостью решения в зоне разрушения приводили к вынужденному сгущению сетки. Отмечается недостаток экспериментальных данных с учетом вариативности свойств льда. В работе [28] указывается на возможное значительное влияние анизотропии морского льда на характер его деформирования, что пока подробно не исследовано. Отмечаются непростые условия применения МКЭ, приспособленного для континуальных сред при решении задач хрупкого разрушения. Отдельно предполагается, что предпочтительнее использовать дискретные бессеточные методы, такие как, например, гидродинамика сглаженных частиц.

Гидродинамика сглаженных частиц

Гидродинамика сглаженных частиц (Smoothed Particle Hydrodynamics - SPH) представляет собой бессеточный метод, в котором среда заменяется набором частиц, т. е. параметры среды распределяются между частицами [37]. Первоначально разработанный для астрофизики метод применяется во многих областях и прежде всего в вычислительной гидродинамике. Активно ведутся работы по адаптации метода для механики твердых тел [38-40].

На рис, 1, а приведена схема дискретизации области Q со смещением частицы i на вектор и. Для i-й частицы выделена область распространения функции W(r, h), которая представлена на рис. 1, б. Положение в пространстве меняется с r0 до rd (радиус-векторы i-й частицы), расстояние h называют сглаживающей длиной.

а) б)

б£1

а - схематичное представление дискретизации области О и ее деформации; б - схематичное представление нормированной сглаживающей функции ядра Ж(г, к)

Уравнение плотности тела в формулировке гидродинамики сглаженных частиц записывается через сглаживание параметров соседних частиц:

р(г)=^тгЖ(г - г, Ь). (7)

]

Уравнение состояния упругого тела, с учетом уравнения неразрывности среды, имеет вид [30]:

dva 1 д ab , a

dF=1 а?° + g ■

(8)

где уа - а-я компонента скорости; оаЬ - тензор напряжений; ga - объемные силы.

Уравнение движения для i-й частицы примет вид:

dW:

dva v

—т~ = > m 4

dt ^ 4

j

f ab ab

++П 5ab

2 2 4

dx,

ij a j+g ■

(9)

Выражение (9) дополнено искусственной вязкостью Пу для стабилизации реше-

ния (5ab - символ Кронекера).

Если в сеточных методах применительно к моделированию твердых тел существенен вопрос - как разделить элементы сетки при разрушении, то в бессеточных методах - как удержать частицы вместе.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Метод позволяет моделировать большие деформации и разрушение среды путем разделения скопления частиц. Это объясняет возросший интерес к данному методу в области моделирования хрупкого льда. В работе [41] представлено решение задачи разрушения образцов льда в двухмерной постановке задачи при четырехточечном изгибе и одноосном сжатии. Для пластической деформации используется критерий Друкера-Прагера, для валидации модели упругого тела применяется сравнительный расчет упругих колебаний консольно закрепленной балки.

Качественно характер разрушения точно совпадает с экспериментальными данными. Количественно же потребовалась корректировка модели, авторами работы [41] вводится интерполяция полей деформации по первой производной во времени для устранения неустойчивых колебательных процессов в частицах - данная особенность метода может вносить нефизичные эффекты в деформацию твердых тел.

В работе [42] метод получил развитие: проведено двухмерное моделирование разрушения ледяной балки при трехточечном изгибе, прогрессирующее разрушение балки об препятствие и разрушение равномерного ледяного покрова при прохождении судна с разными скоростями (но без моделирования воды). Эксперимент для сопоставления проводили в бассейне на масштабной модели. Получено хорошее соотношение экспериментальных и расчетных данных. В вычислительном эксперименте получена более консервативная оценка нагрузок, действующая на корпус модели корабля. Метод гидродинамики сглаженных частиц с точки зрения авторов работы [42] показывает удовлетворительную точность и устойчивость решения к изменению входных параметров. Предполагается дальнейшее развитие моделей, дополнение их водной средой и расчет системы «лед-судно-вода».

Метод гидродинамики сглаженных частиц больше подходит для моделирования водной среды. В ряде работ [43-46] лед рассматривается как упругий объект, движущийся на поверхности воды. Данный подход больше подходит для крупномасштабного исследования, анализа переноса больших масс льда, оценки давления жидкости на ледяной покров, без рассмотрения локальной деформации и разрушения.

Метод дискретных элементов

Другим бессеточным методом, применяемым для анализа хрупкого разрушения льда, является метод дискретных элементов (МДЭ). В нем физическое состояние системы складывается из физических состояний большого числа дискретных элементов. Набор элементов образует исследуемую среду, но каждый из них является независимым и взаимодействует с соседними элементами. Общее (макроскопическое) поведение системы является результатом взаимодействия отдельных частиц. В 1956 г. был разработан метод для моделирования частиц в молекулярной динамике [47], в 1971 г. метод был адаптирован для исследования динамики горных пород и назван методом дискретных элементов [48]. Теоретическое обоснование метода дано в работе [49] и показано, что его можно считать обобщением МКЭ. В настоящее время этот метод активно применяется для моделирования сыпучих сред [50].

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

частицами. Применение метода дискретных элементов для моделирования непрерывных сред потребовало выработки специальных подходов моделирования, нетипичных для исходного метода. Метод дискретных элементов был дополнен моделями связей в виде балок, соединяющих частицы. В зарубежной научно-технической литературе выделяют модели для объемных тел, подчеркивая их обособленность от исходного метода, - модель из связанных частиц (рис. 2). Связь между частицами представляет собой виртуальную балку, которая ограничивает тангенциальное I и нормальное п смещение частиц относительно друг друга (рис. 3).

(10)

= -упБпЛЫ ; = ^АЛ5?; ЬИп = -фnS¿5t;

5М = К/ 5t,

где ¥„, ^ - проекции результирующей силы в направлениях п и Мп - нормальный момент; ^ -

сдвиговый момент; А=лт62 - площадь сечения; гь - радиус жесткой связи; / = \ пг* - момент

инерции сечения; & - шаг по времени; «п и « - жесткости в нормальном и поперечном направлении; \п, - скорости в направлениях п и ю„, - угловые скорости вокруг соответствующих осей п и

Рис. 2. Модель связанных частиц тела: а - дискретные элементы; б - визуализация балочных связей между элементами

Разрыв связи между частицами происходит при условиях

°max< л + г rb'

A J (11)

F Мл.

Тmax < A + J rb ■

После разрыва связи частицы взаимодействуют как независимые твердые тела.

Эта модель была представлена в работе [51], разработана она для моделирования бетонных конструкций и скальных пород. На ее основе в настоящее время делаются попытки построить модели для льда. Одной из первых работ по оценке применения МДЭ для моделирования льда можно указать статью [52], в которой проведено моделирование одноосного сжатия и трехточечного изгиба простейших образцов морского льда прямоугольной формы в трехмерной постановке задачи. Полученные откалибро-ванные мезопараметры материала распространены на модели массивных льдин, воздействующих на опоры шельфовых установок, проведена оценка их воздействия [53]. Материал рассматривается как изотропный. Зависимость температуры, солености и анизотропии не исследуется, но предполагается в будущих исследованиях.

Метод дискретных элементов подходит для моделирования льдин с учетом того, что он адаптирован к независимому движению большого числа элементов. Причем геометрическая форма данных элементов не обязательно должна представлять собой сферу, что является простейшим случаем. В целом ряде работ [54-57] рассматривается подобное моделирование ледяного покрова, его раскалывание и воздействие на судно или опоры. Возможно также использование нескольких методов в единой модели: МДЭ - для хрупкого льда, МКЭ - для конструкций, на которые воздействует лед.

Как видно из данных приведенных работ, использование МДЭ перспективно в моделировании механики льда, но пока разработка этого направления только начинается.

Обсуждение и заключения

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

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

На основе проанализированных работ можно сделать некоторые обобщающие выводы.

Стандартизированная методика моделирования льда и композиционных материалов на его основе в настоящее время не разработана.

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

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

рассмотренная модель Джонсона-Холмквиста. Более простые модели с учетом повреждаемости демонстрируют более стабильный результат. Модель с повреждаемостью Леметра в нескольких исследованиях показала адекватный и удовлетворительный результат.

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

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

Метод дискретных элементов выглядит перспективным для применения в моделировании деформации и разрушения льда. Решаются задачи и в двухмерной, и в трехмерной постановке задачи. Имеется возможность построения мультифизичных моделей сложных систем «вода-лед-конструкция» (объединение вычислительной гидродинамики, МДЭ и МКЭ). При этом в методе имеется возможность применения неоднородных и несферических частиц. В дальнейшем также могут быть усовершенствованы модели связей частиц для более точного моделирования сложных процессов разрушения.

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

Исследование выполнено при финансовой поддержке Российского научного фонда (грант №18-13-00392).

Библиографический список

1. Орешко Е.И., Ерасов В.С., Лашов О.А., Подживотов Н.Ю., Качан Д.В. Численное исследование несущей способности слоистого материала // Все материалы. Энциклопедический справочник. 2019. №3. С. 16-21.

2. Орешко Е.И., Ерасов В.С., Луценко А.Н. Математическое моделирование деформирования конструкционного углепластика при изгибе // Авиационные материалы и технологии. 2016. №2 (41). С. 50-59. DOI: 10.18577/2071-9140-2016-0-2-50-59.

3. Орешко Е.И., Ерасов В.С., Подживотов Н.Ю., Луценко А.Н. Расчет на прочность гибридной панели крыла на базе листов и профилей из высокопрочного алюминий-литиевого сплава и слоистого алюмостеклопластика // Авиационные материалы и технологии. 2016. №1 (40). С. 53-61. DOI: 10.185.77/2071-9140-2016-0-1-53-61.

4. Яковлев Н.О., Гриневич Д.В., Мазалов П.Б. Математическое моделирование напряженно-деформированного состояния при сжатии сетчатой конструкции, синтезированной методом селективного лазерного сплавления // Вестник МГТУ им. Н.Э. Баумана. Сер.: Естественные науки. 2018. №6 (81). С. 113-127. DOI: 10.18698/1812-3368-2018-6-113-127.

5. Дынкин А.А., Верниковский В.А., Добрецов Н.Л., Бузник В.М., Каблов Е.Н. и др. Научно-технические проблемы освоения Арктики. М.: Наука, 2015. 490 с.

6. Бузник В.М., Каблов Е.Н. Состояние и перспективы арктического материаловедения // Вестник РАН. 2017. Т. 87. №9. С. 827-839.

7. Каблов Е.Н. Материалы нового поколения - основа инноваций, технологического лидерства и национальной безопасности России // Интеллект и технологии. 2016. №2 (14). С. 16-21.

8. Лавров В.В. Деформация и прочность льда. Л.: Гидрометеоиздат, 1969. 206 с.

9. Морской лед: справочное пособие. СПб.: Гидрометеоиздат, 1997. 402 с.

10. Каблов Е.Н. Материалы - основа любого дела // Деловая слава России. 2013. №2 (40). С. 4-9.

11. Паундер Э. Физика льда. M.: Мир, 1967. 189 с.

12. Гриневич Д.Б., Нужный ГА., Бузник Б.М., Яковлев Н.О., Гончарова Г.Ю., Разомасов Н.Д. Разрушение армированных ледовых композиционных материалов при изгибном механическом нагружении // Материаловедение. 2020. №1. С. 18-24.

13. Черепанов Г.П. Механика хрупкого разрушения. М.: Наука, 1974. 640 с.

14. Ерасов Ь.С., Орешко Е.И., Луценко AÄ Площадь свободной поверхности как критерий хрупкого разрушения // Aвиационные материалы и технологии. 2017. №2 (47). С. 69-79. DOI: 10.18577/2071-9140-2017-0-2-69-79.

15. Лобанов B.A. Моделирование льда в задачах с конечноэлементной постановкой // Дифференциальные уравнения и процессы управления: электрон. журн. 2008. №4. URL: https://diffjournal.spbu.ru/RU/numbers/2008.4/article.1.2.html (дата обращения: 20.04.2020).

16. Баландин B.B., Крылов С.Ь., Повереннов Е.Ю., Садовский B.B. Численное моделирование ударного взаимодействия упругого цилиндра со льдом // Проблемы прочности и пластичности. 2017. Т. 79. №1. С. 93-103. DOI: 10.32326/1814-9146-2017-79-1-93-103.

17. Баландин B.B., Брагов AM., Глазова Е.Г., Константинов A.^^, Кочетков A.B., Крылов С.Ь. Численно-экспериментальное исследование процессов ударного нагружения льда // Полярная механика: материалы Третьей Междунар. конф. Bладивосток, 2016. С. 293-305.

18. Глазырин Ь.П., Орлова Ю.Н., Орлов М.Ю., Повереннов Е.Ю. Расчет процесса внедрения крупногабаритного ударника в поликристаллический лед // Известия высших учебных заведений. Сер.: Физика. 2012. Т. 55. №9/3. С. 42-45.

19. Орлова Ю.Н., Орлов М.Ю. Исследование процесса взрывного нагружения льда // Bестник Томского государственного университета. Сер.: Математика и механика. 2015. №6 (38). С.81-89.

20. Глазырин Ь.П., Орлов М.Ю., Орлова Ю.Н. Исследование поведения льда при ударе и взрыве // Материалы, технологии и техника для освоения Aрктики и Сибири: сб. тез. III Междунар. конф. Томск: Изд-во Томск. гос. ун-та, 2019. С. 13.

21. Козин B.M., Bасильев A.C, Земляк Ь.Л., Ипатов К.И. Исследование предельного состояния ледяного покрова в условиях чистого изгиба при усилении армирующими элементами // Bестник Томского государственного университета. Сер.: Математика и механика. 2019. №61. С. 61-69.

22. Ипатов К.И., Bасильев A.C, Земляк Ь.Л., Лесков E.B. Несущая способность ледяного покрова при поверхностном армировании // Bестник ИШ ДBФУ. 2019. №3 (40). С. 3-12.

23. Лесков E.B., Ипатов К.И., Bасильев A.C, Земляк Ь.Л. Исследование применения базальтовых и углеродных композитов для поверхностного армирования ледовых переправ // Bест-ник ПГУ им. Шолом^лейхема. 2018. №1 (30). С. 59-63.

24. Тарануха НА., Bасильев A.G Aнализ критериев предельного состояния конструкций из композитных материалов // Ученые записки Комсомольского-на-Aмуре государственного технического университета. 2015. Т. 1. №3 (23). С. 81-87.

25. Ehlers S., Kujala P. Optimization-based material parameter identification for the numerical simulation of sea ice in four-point bending // Journal of Engineering for the Maritime Environment. 2013. No. 228. P. 70-80. DOI: 10.1177/1475090213486892.

26. Herrnring H., Kellner L., Kubiczek J., Ehlers S. Simulation of Ice-Structure Interaction with CZM-Elements. 2018. URL: https://www.researchgate.net/publication/328611406_Simulation_of_Ice-Structure_Interaction_with_CZM-Elements (дата обращения: 20.04.2020).

27. Herrnring H., Kellner L., Kubiczek J., Ehlers S. A cohesive model for ice and its verification with tensile splitting tests. 2019. URL: https://www.dynalook.com/conferences/12th-european-ls-dyna-conference-2019/adhesive-rivets/03_herrnring_tuhh.pdf (дата обращения: 20.04.2020).

28. Nisja H.A. Numerical Modeling of Brittle Failure in Ice Structures: Master thesis. Norwegian University of Science and Technology, 2014. URL: https://ntnuopen.ntnu.no/ntnu-xmlui/handle/11250/237581 (дата обращения: 20.04.2020).

29. Hallquist J.O. LS-DYNA Theory Manual. Livermore Software Technology Corporation, 2006. 680 p.

30. Gagnon R.E. A numerical model of ice crushing using a foam analogue // Cold Regions Science and Technology. 2011. No. 65 (3). P. 335-350. D01:10.1016/j.coldregions.2010.11.004.

31. Gagnon R.E., Derradji-Aouat A. First results of numerical simulations of bergy bit collisions with the CCGS Terry Fox Icebreaker // IAHR. 2006. No. 49 (3). P. 206-214. DOI: 10.1016/j.coldregions.2007.03.016.

32. Gagnon R.E., Wang J. Numerical simulations of a tanker collision with a bergy bit incorporating hydrodynamics, a validated ice model and damage to the vessel // Cold Regions Science and Technology. 2012. No. 81. P. 26-35.

33. Holmquist T.J., Johnson G.R., Cook W.H. A computational constitutive model for concrete subjected to large strains, high strain rates, and high pressures // The 14th International Symposium on Ballistic. 1993. P. 591-600.

34. Von Bock P.R., Ehlers S. Model scale ice - part B: Numerical model // Cold Regions Science and Technology. 2013. No. 94. P. 74-81.

35. Tomac T., Klanac A., Katalinic M. et al. Numerical simulations of ship resistance in model ice // Developments in Maritime Transportation and Exploitation of Sea Resources. 2013. P. 847-851. DOI: 10.1201/b15813-106.

36. Орешко Е.И., Ерасов В.С., Луценко А.Н. Особенности расчетов устойчивости стержней и пластин // Авиационные материалы и технологии. 2016. №4 (45). С. 74-79 DOI: 10.18577/2071-9140-2016-0-4-74-79.

37. Monaghan J.J. Smoothed particle hydrodynamics // Annual Review of Astronomy and Astrophysics. 1992. Vol. 30. Is. 1. P. 543-574.

38. Benz W., Asphaug E. Simulations of brittle solids using smooth particle hydrodynamics // Computer Physics Communications. 1995. Vol. 87. Is. 1-2. P. 253-265.

39. Deb D., Pramanik R. Failure process of brittle rock using smoothed particle hydrodynamics // Journal of Engineering Mechanics. 2013. Vol. 139. Is. 11. P. 1551-1565.

40. Gray J.P., Monaghan J.J., Swift R.P. SPH elastic dynamics // Computer Methods in Applied Mechanics and Engineering. 2001. Vol. 190 (49-50). P. 6641-6662. DOI: 10.1016/s0045-7825(01)00254-7.

41. Zhang N., Zheng X., Ma Q. Updated Smoothed Particle Hydrodynamics for Simulating Bending and Compression Failure Progress of Ice // Water. 2017. No. 9 (11). P. 3-24. DOI: 10.3390/w9110882.

42. Zhang N., Zheng X., Ma Q., Hu Z. A numerical study on ice failure process and ice-ship interactions by Smoothed Particle Hydrodynamics // International Journal of Naval Architecture and Ocean Engineering. 2019. Vol. 11. Is. 2. P. 796-808.

43. Цветков Д.О. Малые движения идеальной стратифицированной жидкости, частично покрытой упругим льдом // Вестник Удмуртского университета. Сер.: Математика. Механика. Компьютерные науки. 2018. Т. 28. Вып. 3. P. 328-347. DOI: 10.20537/vm180305.

44. Pan W., Tartakovsky A., Monaghan J. A smoothed-particle hydrodynamics model for ice-sheet and ice-shelf dynamics // Journal of Glaciology. 2012. No. 58 (208). P. 216-222. DOI: 10.3189/2012JoG11J084.

45. Staroszczyk R. SPH Modeling of Sea-ice Pack Dynamics // Archives of Hydro-Engineering and Environmental Mechanics. 2017. No. 64 (2). P. 115-137. DOI: 10.1515/heem-2017-0008.

46. Shen H.T., Su J.Sh., Lin L.W. SPH Simulation of River Ice Dynamics // Journal of Computational Physics. 2000. Vol. 165. Is. 2. P. 752-770.

47. Alder B.J., Wainwright T.E. Phase Transition for a Hard Sphere System // The Journal of Chemical Physics. 1957. No. 27 (5). P. 1208-1209. DOI: 10.1063/1.1743957.

48. Cundall P.A. The measurement and analysis of accelerations in rock slopes: PhD thesis. Imperial College London, 1971. 237 p.

49. Williams J.R., Hocking G., Mustoe G.G.W. The Theoretical Basis of the Discrete Element Method // Proceedings from Conference Numerical Methods of Engineering, Theory and Applications. 1985. P. 897-906.

50. Cundall P.A. A discrete numerical model for granular assemblies // Geotechnique. 1979. No. 29 (1). P. 47-65. DOI: 10.1680/geot.1979.29.1.47.

51. Potyondy D.O., Cundall P.A. A bonded-particle model for rock // Journal of Rock Mechanics and Minimg Sciences. 2004. No. 41. P. 1329-1364.

52. Ji S., Di S., Long X. DEM Simulation of Uniaxial Compressive and Flexural Strength of Sea Ice: Parametric Study // Journal of Engineering Mechanics. 2017. No. 143 (1). P. 4-11. DOI: 10.1061/(asce)em.1943-7889.0000996.

53. Ji S. Discrete element modeling of ice loads on ship and offshore structures // Proceedings of the 7th International Conference on Discrete Element Methods. 2017. No. 188. P. 45-54. DOI: 10.1007/978-981-10-1926-5_6.

54. Hisette Q., Alekseev A., Seidel J. Discrete Element Simulation of Ship Breaking Through Ice Ridges // Proceedings from the 27th International Ocean and Polar Engineering Conference. International Society of Offshore and Polar Engineers, 2017. P. 1-15.

55. Neto A.V. Improvement of Rubble Ice Generation in Numerical Simulation of Ice Ridge and Structure Interaction: Master thesis, Liège: Université de Liège, 2017. 77 p. URL: https://matheo.uliege.be/handle/2268.2/4432 (дата обращения: 20.04.2020).

56. Polojârvi A. Sea ice ridge keel punch through experiments: model experiments and numerical modeling with discrete and combined finite-discrete element methods: dissertation PhD. Aalto University, 2013. 60 p. URL: https://aaltodoc.aalto.fi/bitstream/handle/123456789/8957/ isbn9789526050829.pdf?sequence=1&isAllowed=y (дата обращения: 20.04.2020).

57. Tuhkuri J., Polojârvi A. A review of discrete element simulation of ice-structure interaction // Philosophical transactions. Series A: Mathematical, physical, and engineering sciences. 2018. No. 376 (2129). Article No. 20170335. DOI: 10.1098/rsta.2017.0335 URL: https://www.ncbi. nlm.nih.gov/pmc/articles/PMC6107613/ (дата обращения: 20.04.2020).

i Надоели баннеры? Вы всегда можете отключить рекламу.