Моделирование процесса разрушения ледяного покрова под действием атмосферного давления и собственного веса тема диссертации и автореферата по ВАК РФ 01.02.04, кандидат физико-математических наук Сергеева, Анастасия Михайловна

• Сергеева, Анастасия Михайловна
• кандидат физико-математических наук
• 2006

Оглавление диссертации кандидат физико-математических наук Сергеева, Анастасия Михайловна

Введение

Глава 1 Способ разрушения ледяного покрова атмосферным давлением в плавающем контейнере

1.1 Описание способа

1.2 Расчет заполнения контейнера водой

Глава 2 Построение математической модели

2.1 Выбор численного метода решения задачи

2.1.1 Обоснование выбора численного метода

2.1.2 Численный метод решения дифференциальных уравнений механики твердого тела

2.2 Математическая модель

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

2.2.2 Построение численной схемы решения задачи

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

Глава 3 Результаты исследования

3.1 Разрушение льда в стационарном контейнере на воздушной подушке

3.2 Разрушение льда под действием гидростатического давления на стационарном контейнере

3.3 Разрушение льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере

3.4 Разрушение льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере с учетом его заполнения водой

Глава 4 Эволюция разрушения ледяного покрова под действием гидростатического атмосферного давления на подведенном плавающем контейнере с учетом его заполнения водой

Введение диссертации (часть автореферата) на тему «Моделирование процесса разрушения ледяного покрова под действием атмосферного давления и собственного веса»

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

В данной работе рассматривается ледяной покров рек и его разрушение под действием атмосферного давления и собственного веса. Ледяной покров рек неоднороден и характеризуется определенными особенностями строения, зависящими от условий льдообразования [1].

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

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

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

Чисто упругое поведение монокристалла обусловливается главным образом изменениями межмолекулярных расстояний под действием приложенного напряжения. Однако возбужденные напряжением движения дефектов (т. е. дислокации) также вносят свой вклад в деформацию. При движении дефектов к зонам равновесия твердое тело будет непрерывно деформироваться. Эта деформация будет не вполне упругой. Однако если напряжение прикладывается и снимается в течение достаточно короткого промежутка времени (например, при прохождении звуковой волны), дефекты не успевают участвовать в достаточной мере в деформации, которую в этом случае можно считать упругой. По этой причине константы упругости льда, получаемые при высокочастотных акустических измерениях более надежны, чем те же характеристики, получаемые из экспериментов, в которых измеряется деформация тела, испытывающего статически приложенную нагрузку. При температуре от -3 до -40 °С лед ведет себя как вполне упругое тело, которое подчиняется закону Гука, если приложенное напряжение не превышает определенного значения и продолжительность его воздействия достаточно коротка. Это происходит при напряжении сжатия до 0,1 МПа, скорости приложения нагрузки около 0,05 МПа/с и продолжительности воздействия напряжения менее 10 с.

Необходимо отметить, что коэффициент Пуассона уменьшается с понижением температуры примерно от 0,5 при -б°С до 0,38 при -40 °С. В то же время для монокристалла в диапазоне от 0 до -40 °С модуль Юнга и коэффициент Пуассона не зависят от температуры и равны соответственно 8,34*103 МПа и 0,35. Этот факт объясняется тем, что в деформировании поликристаллического льда существенную роль играет зависимость скольжения по границам зерен от температуры, а также возможное обратимое движение дислокации. Также, исследования показали, что структурно-текстурная анизотропия упругих свойств льда не превышает 10%.

Флюктуации плотности могут достигать 0,3% и более. Наименьшая измеренная плотность есть плотность «постаревшего» льда. Увеличение плотности на 0,18% соответствует увеличению расстояния между соседними атомами кислорода на 0,06%. [2].

Динамические исследования модулей упругости на 60 образцах глетчерного льда, извлеченных с разных глубин ледника с помощью термобура, были проведены В. А. Никитиным в районе ст. Новолазаревской (Антарктида). Толщина ледника составляла 460 м. Измерения скорости продольных и сдвиговых волн производились при температуре, близкой к температуре ледника на глубине, с которой был взят образец. Одновременно методом гидростатического взвешивания измерялась плотность льда образцов. Результаты показывают, что имеется тенденция к увеличению модулей упругости с увеличением плотности льда. При этом отклонение значений модулей от осредненных значений в большинстве случаев соответствует изменению плотности льда. Коэффициент Пуассона был равен 0,32-0,34. Модуль Юнга увеличивался по глубине ледника от 7,1*103 до 8,4*1 и МПа, модуль сдвига -от 2,7*10 до 3,1*10 МПа. Некоторые образцы, извлеченные с глубины 370 м и более, растрескивались, вероятно, вследствие снятия высоких напряжений, существующих в глубинных слоях ледника.

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

Природный лед содержит многочисленные повреждения - от субмикроскопических и микроскопических дефектов до крупных пор и магистральных трещин.

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

Вследствие многообразия факторов, определяющих физические свойства льда, данные по его прочностным характеристикам, полученные исследователями многих стран по результатам многих тысяч экспериментов за предшествующие 100 лет, имеют большой диапазон значений. Например, прочность на сжатие изменяется примерно от 4*10? до 130*105 Па (и даже больше), а прочность на изгиб - от 3*105 до 30*105 Па.

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

В настоящее время проблема прочности рассматривается с позиций механической и кинетической концепций.

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

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

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

Среди методов разрушения ледяных покровов особое место занимают методы, основанные на тепловом воздействии. Сравнивая затраты на разрушение единицы объема льда при различных видах воздействия, можно показать, что тепловое воздействие на лед, с целью его разрушения, не является наиболее выгодным с энергетической точки зрения - к примеру, в идеальных условиях, когда вся мощность идет на плавление льда, для полного растопления 1 м морского льда за 1 ч потребуется источник мощностью ~0,1 МВт. Тем не менее, тепловые методы, основанные на использовании дешевой энергии, например, отработанного тепла промышленных или судовых установок, а также атомных источников, могут явиться удобным и экономически выгодным, а в отдельных случаях и единственным средством преодоления ледовых затруднений. Естественно, рационально использовать на таяние льда и огромную даровую энергию Солнца.

Расчеты показывают, что за 1 мин при мощности 7,4 МВт 8,8 м3 льда могут быть нагреты на 18° (от -20 до -2°С). Следовательно, за 1 ч при расходе той же мощности можно нагреть уже около 530 м льда или же площадь 530 м2 при толщине льда 1 м. [2].

Нагревание льда приводит к соответствующему уменьшению его прочности. Дальнейшее нагревание льда Солнцем приводит к его разрушению.

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

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

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

Устройство, разработанное Э.Майером, для увеличения ледопроходимости ледокола, содержит ледоразрушающий агрегат, установленный на платформе, укрепленной в носу ледокола над ледяным полем. Имеет ширину, превышающую ширину ледокола. Ледоразрушающий агрегат выполнен в виде разнесенных по ширине платформы турбореактивных двигателей, от которых турбореактивная струя направлена на ледяной покров под углом в сторону кормы.

Устройство, разработанное Евдокимовым С.И. и Каштеляном В.Ш., для обеспечения ледопроходимости судов, имеет следующее строение: по бортам судна на уровне контакта корпуса со льдом устанавливаются сопла, подающие смесь горячего газа от судовой энергетической установки. Особенностью устройства являются струйные аппараты, связывающие сопла с энергетической установкой, подсасывающие воздух из атмосферы, причем соединительные трубопроводы сопел снабжен запорно-регулирующей аппаратурой.

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

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

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

При температуре -5°С 1 г бикарбоната калия может растопить 59 г льда, фторида натрия-33 г, сульфида натрия-27 г. При более низкой температуре рекомендуется применять другие соли: так, в диапазоне от -6 до -20 °С используют хлористые соли аммония, натрия и калия.

В месте контакта химического вещества со льдом образуется очаг рассола, который по мере повышения температуры углубляется. Наибольшее количество льда выплавляется при резких повышениях температуры. Хлористые соли натрия, калия и аммония плавят лед преимущественно по вертикали, тогда как у хлорида магния разрушительное действие проявляется и «в стороны».

Существует зависимость расхода химических веществ от температуры и толщины льда. При разрушении льда толщиной 10 см па. площади 1 га при температуре -10°С требуется израсходовать 3,5 т соли с кусками массой не менее 10 г. Для разрушения слабого речного предпаводкового льда той же толщины и на той же площади достаточно всего 0,5 т мелкозернистой соли.

Расход химикатов в теплую погоду с устойчивой положительной средней суточной температурой воздуха может быть уменьшен до 50%. При сплошном (без ледяных перемычек) протаивании льда, нормы расхода химических веществ обычно в 7 - 10 раз меньше массы льда, подлежащего растворению.

Ширина полосы химикатов при посыпке с автомашины составляет 0,30,4 м, с самолета 2,5-3,5 м. Применение авиации позволяет повысить степень механизации работ, но расход химических веществ увеличивается.

Опыт использования авиационно-химических способов при, очистке от льда подходов к порту показал, что при расходовании 50-100 т каменной соли на 1 га разрушение льда происходит на 15-20 суток раньше, а себестоимость составляет около 2 % стоитимости разрушения льда с помощью ледоколов. Химический метод как самостоятельный целесообразно применять на ограниченной площади для местного разрушения льда. Обычно этот метод следует рассматривать как вспомогательный и его рекомендуется применять совместно с работой ледорезных машин на участках рек с повышенной толщиной льда (толщина ) льда больше длины фрезы) и с включением большого количества бревен и других твердых предметов, а также в местах пересечения ледорезных трасс.

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

4-4,5 см через сутки после их нанесения на лед могут проникнуть в него на глубину от 20 до 70 см.

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

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

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

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

В США изучались возможности применения лазеров для разрушения льда [2]. В эксперименте по резанию образцов льда, полученных из дистиллированной и водопроводной воды путем смораживания крупной ледяной крошки диаметром до 1 см; использовался лазер непрерывного действия на двуокиси углерода мощностью 50 Вт. Пучок фокусировался зеркалом до плотности потока мощности 50 Вт/см . Применялся также импульсный лазер с энергией импульса 0,5 Дж и длительностью 1 мкс.

Лазер непрерывного действия резал лед с эффективностью 3-5 см /с. По оценкам американских ученых, лазер мощностью 50-100 Вт может резать ледяное поле на глубину 8 см со скоростью 0,514 м/с. Эффективность резания падает при горизонтальном положении образца, когда образующаяся при таянии вода не может вытекать из разреза.

Вследствие большой теплоемкости льда использование лазера для растопления больших объемов льда представляется нереальным.

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

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

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

Танкер-ледокол (пат. 3768427 США) предназначен для плавания как на чистой воде, так и в ледовых условиях. Он имеет широкий корпус и узкую надстройку. При плавании во льдах принимается водяной балласт с таким расчетом, чтобы ледорезы, установленные на рампе, обеспечили эффективное разрушение льда.

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

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

Для удаления ледяной чаши могут быть использованы потокоподогреватели [15]. Натурные опыты по удалению ее с помощью паровых потокоподогревателей проводились в 1971 г. на р. Северной Двине в г. Архангельске. С этой целью потокоподогреватели были установлены на плотиках перед носовой частью танкера, имеющего припай по ватерлинии и мощную ледяную чашу на подводной части корпуса. Однако из-за ограниченного времени и расхода водяного пара этот опыт не был завершен. За 42 ч работы ледяная чаша и припай были растоплены только на одной трети длины судна в его носовой части.

Наиболее эффективными способами разрушения ледяного покрова, как показала практика, являются способы, связанные с резонансными эффектами, или с использованием специальных судов, закачивающих под лед воздух, и последующее разрушение льда происходит под действием его массовых сил (под действием своей силы тяжести). Следует отметить, что хотя данные способы и считаются в настоящее время весьма эффективными, но являются весьма дорогими, да и эффективность их весьма относительна. Резонансный способ эффективен только при наличии толщины ледового покрова до 0,7 м [18]. Способ, связанный с ломкой льда при закачке под него воздушной подушки, также требует значительных затрат. Чтобы разрушить ледяной покров значительной толщины (1-3) м, необходимо построить плавающий компрессор стоимостью во много миллионов долларов, при этом необходимая для судоходства производительность не позволит реально его использовать для образования необходимого фарватера.

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

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

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

В качестве объекта исследования взят ледяной покров.

Предметом исследования явилось разрушения ледяного покрова.

В работе ставились и решались следующие задачи:

- построение пространственной математической модели разрушения ледяного покрова под действием собственного веса и атмосферного давления;

- разрушение льда в стационарном контейнере на воздушной подушке;

- разрушение льда под действием гидростатического давления на стационарном контейнере;

- разрушение льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере;

- разрушение льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере с учетом его заполнения водой;

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

- рассмотреть процесс развития трещины в ледяном покрове;

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

На защиту выносятся следующие основные положения:

-пространственная математическая модель разрушения ледяного покрова под действием атмосферного давления и собственного веса; -разрушение льда в стационарном контейнере на воздушной подушке; -разрушение льда под действием гидростатического давления на стационарном контейнере; -разрушение льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере; -разрушение льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере с учетом его заполнения водой;

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

Заключение диссертации по теме «Механика деформируемого твердого тела», Сергеева, Анастасия Михайловна

Основные результаты, полученные в диссертации

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

2. Построена пространственная математическая модель разрушения ледяного покрова под действием собственного веса и атмосферного давления;

3. Составлены блоки программы на языке Visual Fortran 6.1, которые дают возможность определять напряжения и перемещения при решении пространственной задачи;

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

5. Решена задача о разрушении льда под действием гидростатического атмосферного давления на стационарном контейнере;

6. Решена задача о разрушении льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере;

7. Решена задача о разрушении льда под действием гидростатического атмосферного давления на подведенном плавающем контейнере с учетом его заполнения водой;

8. Рассмотрен процесс развития трещины в ледяном покрове;

9. Выведена формула зависимости длины контейнера от толщины ледяного покрова и внешних характеристик его привода.

Список литературы диссертационного исследования кандидат физико-математических наук Сергеева, Анастасия Михайловна, 2006 год

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

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

3. Лурье А.И. Пространственные задачи теории упругости. М: Гостехиздат, 1955

4. Быковцев Г.И., Ивлев Д.Д Теория пластичности. Владивосток: Дальнаука. 1998

5. Прагер В. Введение в механику сплошных сред. ИЛ. 1963

6. Шевченко К.Н. Основы математических методов в теории обработки металлов давлением. М: Высшая школа. 1970

7. Лурье А.И. Теория упругости. М: Наука. 1970

8. Тимошенко С.П. Теория упругости. Л: Судпромиздат. 1985

9. Лейбензон Л.С. Курс теории упругости. М: Гостехиздат. 1947

10. Новожилов В.В. Теория упругости. Л: Судпромиздат. 1958

11. Ляв А. Математическая теория упругости. М: ОНТИ. 1935

12. Христианович С.А., Михлин С.Г., Девисон Б.Б. Некоорые вопросы механики сплошных сред. М: Изд.АН СССР. 1938

13. Мусхелишвили Н.И. Сигулярные и интегральные уравнения. М: Наука . 1968

14. Колмогоров В.Л., Макотра О.А., Моисеев Н.Я Математическая модель для численного решения нестационарных задач механики твердого тела модифицированным методом Годунова.// ПМТФ Т.45 №1. 2004 С. 66-72

15. Быков Л., Клюев В. Судоходные шлюза на морском пути Св. Лаврентия. Речной транспорт, 1976, № 2, с. 52-55

16. Тихомиров В.М., Суровин П.Г. Развитие усталостных трещин смешанного типа в образцах из стали.// ПМТФ Т.45 №1. 2004.С 135-142.

17. Богородский В.В., Гаврило В.П. Лед. Л.: Гидрометеоиздат, 1980.

18. Козин В.М., Онищук А.В.,Марьин Б.Н., Иванов Ю.Л., Повзык Н.Г., Шпорт В.И. Ледоразрушающая способность изгибно-гравитационных волн от движения объектов. Владивосток: Дальнаука. 2005.

19. Ильюшин А.А Ученые записи МГУ// Механика. Вып.39. 1949.

20. Гаврило В.Л., Никитин В.А. Попов И.К. К вопросу об определении прочности ледяного покрова на изгиб в натурных условиях// Труды Аркт. И Антарк. НИИ. 1984. Т. 386. С 39-43.

21. Гладков М.Г., Федеров Б.А. Исследование механических свойств морского льда методом одноосного сжатия // Труды ААНИИ. 1984, Т. 386.С. 44-52

22. Максименко ВН., Тягний А.В Расчетно-экспериментальный метод определения параметров разрушения конструкций с трещинами // ПМТФ Т.45 №4. 2004.С 168-175.

23. Лейбензон Л.С. Вариационные методы решения задач теории упругости. М: Гостехиздат. 1943

24. Михлин С.Г. Вариационные методы в математической физике. М: Наука. 1970

25. Михлин С.Г. Численная реализация вариационных методов. М: Наука. 1966

26. Ионов Б.П. Ледовое сопротивление и его составляющие. Л: Гидрометеоиздат. 1988.

27. Канн С.И. Морские льды. Л: Гидрометеоиздат. 1974

28. Тарновский И .Я., Поздеев А.А., Ганаго О.А. Деформация и усилия при обработке металлов давлением. Машгиз. 1959

29. Тарновский И.Я., Поздеев А.А., Вайсбург Р.А., Гунн Г.Я., Кошельников В.Л., Тарновский В.И., Скороходов А.Н., Коломогоров В.Л. Вариационные принципы механики в теории обработки металлов давлением. Металлургиздат. 1963

30. Моисеев Н.Н Численные методы в теории оптимальных систем. М: Наука. 1973

31. Черноусько Ф.Л., Баничук Н.В. Вариационные задачи механики и управления. М: Наука .1973

32. Шутов А.В. Анализ двух нелинейных моделей хрупкого разрушения твердых тел // ПМТФ Т.45 №6. 2004.С 95-102.

33. Самарский А.А. Тория разностных схем. 2-е изд. М: Наука. 1983

34. Марчук Г.И. Методы вычислительной математики. 2-е изд. М. 1980

35. Самарский А.А., Николаев Е.С. Методы решений сеточных уравнений. М: 1978

36. Стренг Г., Фикс Дж. Теория метода конечных элементов. М: Мир 1977

37. Леонтьев Л.В Вариационно-сеточный метод решения задач о собственных колебания упругих трехмерных тел, связанный с использованием ортогональных финитных функций// МТТ. № 3. 2002 С 117-126

38. Солодовников В.Н. К Алгоритму решения задач пластичности методом конечных элементов // ПМТФ Т.34 №6 1993

39. Легеньков А.П. Деформации дрейфующего льда в Северном Ледовитом океане. СПб.: Гидрометоеиздат. 1992

40. Маэно Н.Наука о льде: Пер. с яп. М: Мир. 1988

41. Жестка В.Д. Исследование напряженно-деформированного состояния ледяного покрова, находящегося под действием движущейся нагрузки, в условиях мелководья. // ПМТФ Т.41 №4 2000

42. Вильчинский А.В, Чугунов В.А. Моделирование течения льда в различных зонах ледников.//ПММ Т.65 вып.З 2001. С. 495-510

43. Бреббия К., Телес Ж., Вроубел Л. Методы граничных элементов. М: Мир. 1987

44. Зенкевич О.С. Методы конечных элементов в технике. М: Мир 1975

45. Колмогоров В.Jl. Численное исследование больших пластических деформаций и разрушения металлов //КШП. ОМД. №2 2003. С 4-16.

46. Одиноков В.И. О конечно-разностном представлении дифференциальных соотношений теории пластичности// Прикладная механика. 1985. Т.21 №1 С.97-102

47. Леонтьев В.Л. Метод конечных элементов теории упругости (смешанные вариационные формулировки). Ульяновск: Изд-во Средневолж. Науч. центра 1998

48. Песчанский И.С. Разрушение припайных льдов силами внешнего воздействия//Проблемы Арктики 1946 №1

49. Работнов Ю.Н. Механика деформируемого твердого тела. М: Наука 1988

50. Слепан Л.И. Механика трещин. Л: Судостроение. 1990

51. Смирнов В.Н. Некоторые вопросы натурного исследования деформаций и напряжений в ледяном покрове// Труды ААНИИ 1983 С.92-96

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

53. Качанов Л.М. Основы теории пластичности. М: Наука. 1969

54. Дж Мэйз Теории и задачи механики сплошных сред. М: Мир 1974

55. Тимохов Л.А. Хейсин Д.Е. Динамика морских льдов. Л: Гидромтеоиздат, 1987

56. Хейсин Д.Е. Динамика ледяного покрова. Л: Гидромтеоиздат. 1967

57. Хейсин Д.Е. Динамика ледяного покрова.// Механика и физика ледяного покрова. М.: Наука. 1983. С . 152-163

58. Храпатый Н.Г., Тахтеев В.А. О механизме образования трещин в ледяном покрове//Гидротехнические сооружения: Сб. науч. трудов. Владивосток: ДВГУ. 1981 С. 85-88

59. Полярус A.M., Романов Д.Ю. Об одном способе разрушения ледяного покрова // Проблемы механики сплошных сред и смежные вопросы технологии машиностроения: Сб. докладов второй конференции.

60. Владивосток, 31 августа 6 сентября 2003 г. Комсомольск-на-Амуре: ИМиМ ДВО РАН, 2003. С. 23-28.

61. Бердянников В.П. Изучение модуля упругости льда. Труды ГПИ, 1948. Вып. 7(61). С. 13-23.

62. Меркулов В.И., Одиноков В.И., Ловизин Н.С. Об одном подходе к численному решению задач упругопластического деформирования тел пространственной формы // Кузнецно-штамповочное производство. Обработка металлов давлением. 2001. № 6. С.12-19.

63. Макеранец Е.И., Одиноков В.И. Расчет пластического течения полых овальных цилиндров неограниченной длины // Известия АН СССР. Механика твердого тела. 1976. № 2. С. 103-110.

64. Одиноков В.И., Сергеева А.М Математическое моделирование одного нового процесса разрушения ледяного покрова // Прикладная механика итехническая физика 2006. - №2. - с. 139-146.

65. Сергеева A.M. Разрушение льда атмосферным давлением // Вестник ГОУВПО КнАГПУ.- 2004. с 93-96.

66. Odinokov V.I., Sergeeva. A.M. Mathematical Modeling for One New Method of Breaking Ice Cover // Journal of Applied Mechanics and Technical Physics. Springer New York. - Volume 47. №2 - 2006. - pp.266-273

67. Патент РФ № 2220878. Способ разрушения ледяного покрова / Одиноков В.И., Козин В.М. Бюл. № 1. Опубл. 10.01.04 г.

68. Сергеева A.M. Математическое моделирование процесса разрушения ледяного покрова под действием гидростатического атмосферного давления.// Деформирование и разрушение структурно-неоднородных сред и конструкций: Сборник докладов.- Новосибирск-2006. с-.

69. Свидетельство об официальной регистрации программы для ЭВМ № 2006610010. Математическое моделирование процесса разрушения ледяного покрова / Одиноков В.И., Сергеева A.M., Жигалкин К.А. Зарегистрировано в Реестре программ для ЭВМ 10 января 2006.