Математическое моделирование разрушения ледяного покрова от сжимающих усилий

Размещено на http://www.allbest.ru/

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАЗРУШЕНИЯ ЛЕДЯНОГО ПОКРОВА ОТ СЖИМАЮЩИХ УСИЛИЙ

Сергеева А.М., Одиноков В.И., Захарова Е.А., Марченко О.В.

Комсомольск-на-Амуре, Россия

Достаточно серьезной проблемой является обеспечение передвижения судов в условиях ледяной корки и значительного ледяного покрова. Существует не мало экспериментальных исследований, связанных с разрушением льда. Фундаментальные теоретические исследования по ломке ледяных полей имеют место, но, в большинстве случаев, лед моделируется как ледяная пластина. В данном случае проводится теоретическое исследование работы устройства в виде ледокольной приставки, защищенное патентом РФ [1]. Работа направлена на построение математической модели, с помощью которой проводится численный анализ процесса разрушения ледяного покрова. А также исследуется зависимость разрушения льда от геометрических параметров исследуемого устройства. Разрушение льда происходит как от растягивающих напряжений, появляющихся во льду при внедрении клина ледокольной приставки, так и за счет напряжений сжатия между двумя клиньями, реализующих потерю устойчивости ледяной пластины в поперечном направлении. Сжимающие усилия должны достигать таких величин, при которых ледяная пластина теряет устойчивость и разрушается. Для обеспечения более быстрой потери устойчивости используется специальная поперечная нагрузка.

Устройство для разрушения льда представлено на рис.1

В носовой части судна 1 при помощи креплений 2 устанавливают ледокольную приставку 3, перед которой в ледяном покрове 4 создают его свободную кромку 5. Приставку 3 выполняют в виде плавучей емкости, в диаметральной плоскости 6 которой устанавливают форштевень 7, наклоненный под углом к горизонтальной поверхности, а по бортам - боковые ножи 8 с углом в плане ? и наклоненные к горизонту под таким же, но отрицательным по отношению к форштевню углом. Благодаря такой геометрической форме форштевня и боковым ножам при их контакте со свободной кромкой 5 в ледяном покрове между ножами 8 будут возникать сжимающие усилия 9. При определенном расстоянии l и усилиях 9 это приведет к потере устойчивости формы участка льда 10 и его разрушению. Отломанный участок льда 10 от сплошного льда 4 притопится наклонным днищем приставки 3, и его обломки бортами 11 судна 1 раздвинутся под кромки образовавшегося канала 12.

деформация трещина лед напряжение

В работе решается пространственная задача о деформации ледяного покрова при одновременном внедрении в лед двух клиньев и вертикальной нагрузки, приложенной по оси симметрии. Для проведения теоретических исследований использовался численный метод [2]. Деформируемую среду (клин и лед) будем полагать упругой и изотропной, таким образом, решается двухкомпонентная система. Учитывая симметрию, рассматривается половина области деформирования. Массовые силы не учитываются. Используя теорию упругости для малых деформаций и Эйлерову систему координат, запишем систему дифференциальных уравнений.

(1)

(2)

(3)

Уравнение теплопроводности (только для льда)

(4)

Здесь G_t - модуль сдвига; t=1, 2; G₁ - модуль сдвига клина; G₂ - модуль сдвига льда; G₂ = G(q); ? - температура; k_t - коэффициент объемного сжатия; k₁ - клина; k₂ - льда; k₂ = k(q); [s_ij] - тензор напряжений; [e_ij] - тензор деформаций; u_i - проекции перемещений по координатным осям х_i, i = 1, 2, 3; ? - коэффициент теплопроводности.

Уравнение теплопроводности (4) записано для стационарного случая: Ледяной покров можно рассматривать как пластину конечной толщины, в которой При при Тогда уравнение (4) примет вид

(5)

Коэффициент ? изменяется от температуры по линейному закону [3]

(6)

Интегрируя уравнение (5) и учитывая, что , получим

(7)

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

(8)

где Р₁ - атмосферное давление; Р₂ - давление, оказываемое на нижнюю часть льда, контактирующую с водой; Р₃ - нагрузка, находящаяся на поверхности льда; h₀ - толщина льда; - угол наклона клиньев; b₁ - половина ширины нагрузки; l₁ - ширина рассматриваемого льда; l₂ - лед, контактирующий с неподвижной основой, u^* - заданное перемещение клина.

Для решения задачи разработан алгоритм:

1. Задаются начальные условия и шаг по времени Dt.

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

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

4. Насчитывается поле температур по каждому элементу, входящему в ледяное поле; поле температур по элементам инструмента задается как ?*.

5. Насчитываются значения и по каждому элементу (n - номер элемента); в инструменте и - задаются для стали.

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

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

8. По каждому элементу (его граням) (ij) насчитываются .

9. Уточняются граничные условия по поверхности S₆ в соответствие со следующими условиями. Если переназначение , образуется полость (рис.3) с соответствующими граничными условиями: Если операция 10.

10. По напряжениям в плоскостях I, II (рис. 2) исследуется наличие и величина магистральной трещины из условия - предел прочности льда при растяжении. Если это неравенство выполняется, то следует перестройка сетки, в которой строится конфигурация трещины, на соответствующих поверхностях которой нормальные напряжения и касательные напряжения

11. Если перестройка элементов и граничных условий имела место, то выполняется операция 6. В случае отсутствия изменений следует операция 12.

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

При теоретическом моделировании использовались свойства пресноводного льда. Следует отметить, что пресноводный лед прочнее морского.

При нагрузке, перпендикулярной направлению длинных осей кристаллов модуль Юнга (Е) равняется [3]

Коэффициент Пуассона ?

Коэффициент объемного сжатия (k)

Модуль сдвига (G)

Для инструмента

На рис. 4 представлены результаты численного решения при h₀ = 1м, b₂ = 20 м, b₀ = 14 м, l = 130 мм, b₁ = 2 м, a = 30⁰, р₃ = 2.5 атм; р₁ = 1 атм, р₂ = 1.092 атм. Радиус конца клина 100мм. Значение р₃ включает значение р₁. Эпюры на рис.4 приведены в сечении Показано на рис. 4 начало внедрения. Получаем огромные сжимающие напряжения по поверхностям S₆, S₇. На поверхности S₈ нормальное напряжение = - 300 МПа. Такие напряжения превосходят предел прочности при сжатии , и лед при внедрении будет ломаться. Кроме того, на магистральных поверхностях I, II (рис. 2) имеем на некотором расстоянии от кромки клина растягивающие напряжения , превосходящие предел прочности льда при растяжении Следует автоматическая перестройка сетки с изменением граничных условий.

На рис. 5,а (сплошными линиями) приведены эпюры перестроенной системы при Затемненные участки - пустоты. У основания клина пустота образовалась в результате появления растягивающих напряжений, получившихся при изменении граничных условий из-за наличия магистральной трещины. Как только появилась магистральная трещина - сжимающие напряжения по поверхностям клина сильно снизились. По толщине льда в разных сечениях приведены эпюры на рис. 5,б. Видим, что напряжения в области клина по высоте h₀ сжимающие с максимальным значением на поверхности (x₁ = h₀). На рис. 5,в приведены эпюры по толщине h₀ в различных сечениях. С удалением от конца клина по координате x₃ сжимающие напряжения уменьшаются, а в сечении 2 на поверхности x₁ = h₀ переходят в растягивающие. В сечении III-III напряжения по сечениям 1, 2, 3 - растягивающие (рис.5,г) с наибольшим значением в сечении 2, отстоящим от конца клина по координате x₃ на 1660 мм. Пунктирными линиями на рис. 5 показаны эпюры при внедрении клина на 200 мм. Видим увеличение по абсолютным значениям как сжимающих, так и растягивающих напряжений. Линией со звездочкой приведены на рис. 5,г эпюры при внедрении клина на 300 мм. При этом происходит появление поперечных трещин, так как растягивающие напряжения намного превосходят порог прочности льда.

На рис. 6 показана эволюция появления трещин при внедрении клина. На рис. 7 приведены эпюры напряжения , вблизи конца магистральной трещины при внедрении клина на 300 м (линии со звездочкой) и 400 мм (линии штрихпунктирные). На рис. 7,а - эпюры в сечении ; на рис. 5,б - эпюры по высоте х₁ в плоскостях, сопряженных с рис. 5,а и в сечениях I и II.

Видим, что магистральная трещина по концу клина развивается дальше на величину 3320 мм; имеют место и трещины, параллельные магистральной. Растягивающие напряжения превосходят предел прочности на растяжение по всей толщине льда (рис. 5,б). Образуются и поперечные трещины (растягивающие напряжения (рис. 7,в, г)) в основном с правой стороны клина на значительном расстоянии от него по координате x₃.

Литература

1. Патент РФ №2229415. Устройство для разрушения ледяного покрова/ Горкунов Э.С., Колмогоров В.Л., Козин В.М., Одиноков В.И. Бюл. №15. Опубл. 27.05.2004

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

3. Богородский В.В., Гаврило В.П., Недошивин О.А. Разрушение льда. Металлы, технические средства. Л.: Гидрометеоиздат, 1983. 232 с.

Размещено на Allbest.ru