Численное исследование процесса деформирования ледяного покрова ледокольным судном

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2013 Математика и механика № 3

УДК 539.3

А.М. Сергеева, А.В. Ткачева, В.И. Одиноков ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ПРОЦЕССА ДЕФОРМИРОВАНИЯ ЛЕДЯНОГО ПОКРОВА ЛЕДОКОЛЬНЫМ СУДНОМ

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

Ключевые слова: ледяной покров, ледокольное устройство, напряженно-деформированное состояние, численные методы, математическая модель.

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

Работа направлена на исследование напряженно-деформированного состояния (НДС) ледяного покрова, находящегося под воздействием ледокольного устройства [1]. Целью является нахождение оптимального расстояния между надводной и подводной частями устройства, при котором в ледяном покрове возникнут наибольшие разрушения.

Описание устройства

Рассмотрим ледокольное судно, изображенное на рис. 1, которое состоит из надводной 1, подводной 2 и двух боковых частей 3. В надводной части расположены система управления, пассажирский салон, снизу установлены резцы 4 и система фрез 5. Угол наклона форштевня р . Подводная часть имеет спереди резец 6, отсеки для воздуха, топлива, двигателя судна. Боковые части содержат систему коммуникаций, гидропривод перемещения в вертикальной плоскости надводной части относительно подводной, имеют спереди клинья для прорезания льда и вертикальные фрезы. Расстояние между надводной и подводной частями (И) должно быть И > 1,5^, где £ - толщина ледяного покрова. Эта величина устанавливается системой гидроцилиндров, расположенных в боковых частях, перед началом движения судна.

Работа устройства осуществляется следующим образом. Ледокольное судно набирает требуемую скорость на чистой воде и находит на ледяной покров. При

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

Постановка задачи

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

Согласно результатам экспериментальных исследований, описанных в работе [2], лед ведет себя упруго при температурах ниже -5 °С и нагрузках менее 10 с. Примем лед как упругое тело, по данным работы [2, 3] имеем: модуль Юнга -Е = (87,6-0,219-0,001792)-102 МПа ; 9 - температура льда; коэффициент Пуас-

сона - ц = 0,36 ; модуль сдвига - G = E/2(1 + ц); коэффициент объемного сжатия -K = (1 -2ц)/E. В качестве критериев нарушения сплошности ледяного покрова примем напряжения стррс и 4 МПа, ст™ и-8 МПа, приведенные в работе [2].

Используя основные уравнения теории упругости, строится математическая модель устройства. Материал, из которого изготовлено устройство, не рассматривается, поэтому трением пренебрегаем.

Уравнение движения

В уравнениях (1), (2) используется суммирование по повторяющимся индексам; р - плотность льда; G = G(9) - модуль сдвига льда; [стгу ] - тензор напряжений; [6^ ] - тензор деформаций; vi, Fi - проекции скорости перемещений и удельной объемной силы по координатным осям xi, i = 1,2,3 ; F2 = F3 = 0; X - ко-

мации, K - коэффициент объемного сжатия.

Для расчета температурного поля воспользуемся решением уравнения теплопроводности с учетом Х = Х0 (1 + a9), запишем его в конечном виде

При выводе (6) использовались граничные условия: при x1 = 0 9 = 0°С; при x1 = h 91 =-30°С, h - толщина льда, a = -0,0159 град-1.

Рассмотрим 1/2 области льда и 1/2 носовой части ледокольного судна, которые являются симметричными относительно плоскости Ox1 x2, проведенной вдоль

судна. На рис. 2. представлена схема нагружения ледяного покрова.

В плоскости, находящейся на расстоянии h ^ - толщина ледяного покрова) от плоскости Ox2x3, выделим три области: ^ и S4, S16, где S1 - поверхность взаимодействия льда с атмосферой, S4 - площадка контакта форштевня с льдом, S16 - область взаимодействия клина с атмосферой. В плоскости Ox2x3 определим

Закон Гука для малых деформаций

і = 1,2,3; ст=3kє;

Уравнение теплопроводности

эффициент теплопроводности; 9 - температура льда; т>- =—, т - время дефор-

дт

четыре области - £2, $3, £14 , £15 (рис. 2). Площадка £2 характеризует взаимодействие льда с несжимаемой жидкостью - водой. Участки льда £3 и $14 моделируют соответственно взаимодействие подводной части судна с нижней поверх-

нижней поверхности льда, ограниченный областью £15, моделирует взаимодействие клина с водой. Поверхности £6 и £7 расположены параллельно друг другу, причем £7 находится в плоскости симметрии 0x1 Х2, а $6 - на расстоянии Ь от нее. Поверхность £5 расположена в плоскости 0х1 х3 параллельно свободной кромке льда £8. Свободная кромка льда £8 находится в одной плоскости с поверхностью £9, на которой задаются граничные условия, моделирующие жесткое взаимодействие рассматриваемой области деформирования льда с неподвижным ледяным полем. На поверхностях £10, £11, £12, £13, £15, £16 задаются граничные условия, отражающие взаимодействие внедряющегося в лед клина ледокола. Площадки £10 и £12 характеризуют боковые части клина, причем £10 лежит параллельно плоскости 0x1 Х2 на расстоянии (Ь + Ь2) , а $12 находится под углом а к $10. Области £п и £13, расположены параллельно плоскости 0х1 х3, моделируют площадки острия и основания клина соответственно.

ностью льда и неподвижность льда в направлении оси 0x1, т.е. $ = 0 . Участок

а

к% Лед Клин

Рис. 2. Схема нагрузки ледяного покрова: а - в вертикальной плоскости, б - в горизонтальной

Граничные условия задачи (рис. 2):

Здесь р0 - атмосферное давление; у - удельный вес льда; И - толщина льда; V* - скорость передвижения судна; в - угол между носовым отсеком надводной части судна и льдом, принят равным 15°; а - угол острия клина (рис. 2), принят равным 25°.

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

Для решения системы уравнений (1) - (5) с учетом граничных условий воспользуемся алгоритмом, описанным в работе [5, 6] и программным комплексом [7]. Достоинствами данного метода является то, что на границе исследования области можно задать как скорости, так и напряжения, возможно провести расчета как в упругом, так и в пластичном теле, практически любую область можно представить в криволинейных ортогональных элементах, и нет привязки к системе координат. Этот метод довольно прост для представления в различных языках программирования.

Рассматриваемая система уравнений нелинейна. Перепишем в уравнениях (1) инерционный член I, в виде

где vk будет считаться известным, тогда система (1) - (5) будет линейной.

1. Исследуемая область деформирования разбивается на элементы ортогональной формы. Рассчитывается матрица длин дуг элементов.

2. Время исследования (деформирования) т разбивается на т временных промежутков.

3. Насчитывается поле температур по каждому элементу.

4. Насчитываются значения модуля сдвига и коэффициент объемного сжатия по каждому элементу.

5. Задаются начальные условия задачи.

6. Задается шаг по времени Дх.

7. Задаются граничные условия задачи.

8. Задается (V- )0 в (7).

9. Насчитывается матрица коэффициентов и свободных членов новой эквивалентной системы в соответствии с последовательностью, описанной в работе [7].

10. Решается система линейных уравнений методом Гаусса.

11. По каждому элементу (его граням) (//') насчитываются ст^-, V,, ^.

12. Уточняются значения V,. Производится сравнение (V- )п и (V-)п-1, п - номер итерации. Если заданная точность соответствия (V-- )п и (V- )п-1 удовлетворяется, то следует операция 12, в противном случае осуществляется операция 8.

Алгоритм решения задачи

13. Исследуются элементы, непосредственно контактирующие с деформирующим устройством, на наличие растягивающих напряжений; если такие элементы существуют, то на границе льда и устройства ($ - площадка контакта) следует переназначение граничных условий стгг-1£ = 0 ; -,] = 1,2,3 , следует операция 9, иначе 14.

14. Исследуются площадки элементов на наличия напряжений, превышающих предел прочности льда. Если такие элементы присутствуют, то граничные условия перестраиваются ст^ |£ = 0, -,] = 1,2,3 , следует операция 9, иначе 15.

15. Осуществляется шаг по времени х = х + Дхг-, если х<2Дх , то следует опе-

т

рация 5, если - нет, то - операция 16.

16. Окончание расчета.

Результаты исследований

Расчеты проводились с использованием программного комплекса [7]. Язык программирования Фортран, алгоритм линейный. Область деформирования ледяного покрова разбивалась на 2500 ортогональных элементов. Количество уравнений составило 7680. Временной интервал одной итерации, затраченный компьютером с процессором на базе ЮТЕЦК)Соге(ТМ)15СРи 760@2.80вШ, составил примерно 1 час. Для сходимости решения с установленной точностью 0,001 потребовалось до 25 итерационных процессов.

Геометрические параметры устройства: а - расстояние от подводной части судна до окончания исследуемой области деформирования льда в направлении, противоположном положительному направлению оси х2, принято равным 29,0 м; а1 - глубина внедрения клина в лед - принята равной 3,0 м; а2 - расстояние от подводной части судна до надводной - 4,0, 5,0, 6,0 м; а3 - ширина площадки контакта надводной части судна с ледяным покровом - 1,0 м; Ь1 - половина длины площадки контакта надводной части судна с ледяным покровом - 3,0 м; Ь2 - расстояние от надводной части судна до клина - 1,5 м; Ь3 - ширина клина - 0,8 м; Ь4 -длина свободной кромки льда от клина до закрепленной части льда 0,03 м; с -расстояние от свободной кромки льда до края рассматриваемой области деформирования ледяного покрова - принято равным 19,0 м, по всей этой длине лед жестко закреплен с припайным льдом; ё - ширина подводной части судна - 1,0 м; ё - половина длины площадки контакта подводной части судна с ледяным покровом - 0,5 м, ширина острия клина 0,03 м. Скорость устройства 10 км/ч, т = 0,4с , Дт = 0,1с .

Авторы ледокольного судна предполагают, что максимальная толщина льда, по которому способно пройти судна 1,0 м, а рабочая толщина составляет около

0,5 м. Поэтому толщину ледяного покрова принимаем равной 0,5 м.

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

Устойчивость получаемого решения проверялась численно при отклонении от начальных данных на малую величину 0,01 %. В результате данные с веденной погрешностью и без погрешности совпали.

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

тями равно 5,0 м. На рис. 3 представлена эпюра напряжений по толщине ледяного покрова.

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

Л-А„

100 МПа 1--------1

1н>:що,ч ч. НАДВОД. Ч,

ГГ

Б-Б

і____________________________________

и

подвод Ч.

НАДВОД. Ч.

ПОЛНОЛ. Ч.

а^2

/-/

100 МПа

клин [надвод. ч.

ГТТТТТТТТТТТТТТІ ГГП

подвод, ч. клин

|_ _ __І I________

п-п ________г

///-///

І І 11 І І 111 1

КЛИН НАДВОД. ч.

©

а

б

Рис. 3 (начало). Эпюра напряжений ледяного покрова в момент образования магистральной трещины и трещины с наклонного бока клина: а - схема, б - эпюра напряжений до образования магистральной трещины

сг,,

А-А_

20 МПа I-----1

НАДВОД. Ч. - — — 1 , КЛИН

1 1 I

Б-Б

ти

т п

_ подвод, ч.

гггггггп

г клин [нлдвод. ч

11 11 <

Рис. 3 (окончание). Эпюра напряжений ледяного покрова в момент образования магистральной трещины и трещины с наклонного бока клина: в - эпюра напряжений после образования магистральной трещины

Наибольшие из них находятся фронтально от клина и над подводной частью. Нормальные напряжения, перпендикулярные действию устройства, носят растягивающий характер, наибольшие из них расположены у наклонной части клина. Все напряжения с удалением от устройства сглаживаются. Как видно из рис. 3, а, напряжения на порядок превосходят критические, поэтому лед будет трещать и ломаться. Фронтально от клина пройдет основная трещина, ее называют магистральной, она разделит лед на две зоны - внутреннюю и внешнюю. На рис. 4 представлены распределения напряжений в ледяном покрове в момент возникновения магистральной трещины, она выделена белой линией, видно, что концентратор напряжений расположен в вершине трещины. Это свидетельствует о том, что трещина будет расти. От наклонной боковой поверхности клина пойдет трещина, на рис. 4 показано расположение напряжений в этой области, лед здесь разрушится на мелкие фрагменты. После образования магистральной трещины пойдет падение напряжений рис. 3. В ледяном покрове возникнет изгибающий момент, подтверждение этому - смена знака напряжения по толщине льда, расположенного в области над подводной частью судна. Во внутренней зоне ледяного покрова в верхнем слое льда напряжения носят растягивающий характер, а в нижнем слое -сжимающий (эпюра вид 1-1 и 11-11, рис. 3, б) произойдет изгиб льда.

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

в

ва а2 равном 4 м. С увеличением параметра а2 снижаются напряжения в данной области, при этом длина магистральной трещины увеличивается с 3 до 4 м. На рис. 4 трещины изображены прямыми белыми линиями.

Рис. 4. Распределение напряжений в верхнем слои ледяного покрова, при контакте с ледокольным судном, время контакта 0,3 с, параметр устройства а2, м: а - 4, б - 5, в - 6

Далее рассмотрим более подробно нагрузку на составные части судна. Наибольшая нагрузка придется на клин устройства: острие 0,992, 0,813, 0,689 тс, наклонная боковая площадка (S12) - 2,041, 2,034, 2,024 тс соответственно при а2 равном 4,0, 5,0, и 6,0 м. Нагрузка на подводную часть составит 0,0059, 0,0052, 0,052 тс соответственно при а2 равном 4,0, 5,0, и 6,0 м. Нагрузка на надводною часть увеличится с увеличением расстояния а2 составит 0,0086, 0,0074, 0,0067 тс. Видно, что при увеличении расстояния а2 нагрузка на составные части устройства уменьшается.

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

ЛИТЕРАТУРА

1. Пат. 2276037 РФ МПК В 63 В 35/08. Ледокольное судно / В.И. Одиноков, В. Л. Колмогоров. № 2004114052/11; Заяв. 06.05.2004; Опубл. 10.05.2006, Бюл. № 13.

2. Богородский В.В., Гаврило В.П. Физические свойства. Современные методы гляциологии. Л.: Гидрометеоиздат, 1980.

3. Войтковский К.Ф. Основы гляциологии. М.: Наука, 1999.

4. Одиноков В.И. Численное исследование процесса деформации материалов бескоорди-натным методом. Владивосток: Дальнаука, 1995.

5. Сергеева А.М., Марченко О.В. Деформация ледяного покрова под действием атмосферного давления и собственного веса льда // Труды Нижегородского технического университета им. Алексеева. 2010. № 4 (83). С. 157-165.

6. Одиноков В.И., Сергеева А.М., Захарова Е.А. Построение математической модели для численного анализа процесса разрушения ледяного покрова // Математическое моделирование. 2008. Т. 20. № 12. С. 15-26.

7. Одиноков В.И., Сергеева А.М., Ткачева А.В. Математическое моделирование процесса разрушения ледяного покрова ледокольной судном: свидетельство о государственной регистрации программы для ЭВМ № 2012614206. Зарегистрировано в Реестре программ для ЭВМ 12 мая 2012.

Статья поступила 12.12.2012 г.

Sergeeva А.М., Tkacheva A.V., Odinokov V.I. NUMERICAL STUDY OF THE DEFORMATION PROCESS OF THE ICE COVER UNDER THE ACTION OF AN ICEBREAKING SHIP. Using the theory of small elastic deformations, a mathematical model of an icebreaking device is built to study the effect of the device on the ice cover destruction.

Keywords: ice cover, ice-breaking device, the stress-strain state, numerical analysis, mathematical model.

Sergeeva Anastasiya Мikhailovna

E-mail: serg-nasty@mail.ru

Tkacheva Anastasiya Valerievna

E-mail: 4nansi4@mail.ru

Odinokov Valery Ivanovich (Institute for Machine Engineering

and Metallurgy of the Far-Eastern Branch of Russian Academy of Sciences)

E-mail: mail@imim.ru