Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

ASSESSMENT OF THE INTENSE DEFORMED CONDITION OF THE ICE COVER AT THE MOVEMENT OVER ICEBREAKING PLATFORMS ON THE HOVERCRAFT

Larin A.G. 1 Moskvicheva Yu.A. 1
1 Nizhny Novgorod State Technical University n.a. RE Alekseev
The movement of an icebreaking platform with a small speed over an ice cover is considered. Model of the ice field – the thin isotropic uniform plate of constant thickness lying on the elastic basis of hydraulic type. The intense deformed condition of the infinite, semi-infinite plate and a semi-infinite plate having cut in the form of the ice channel is defined.The load of the platform is considered evenly distributed on the area of a rectangle with the varied length relation to width.Influence of size of pressure, form and the area of load of the intense deformed state is investigated. Recommendations are given for designing platforms.
ice cover
stress-strain state
bending
icebreaking platform hovercraft

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

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

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

При разрушении ледяного покрова ЛПВП характеризуются значительными давлениями в ВП (до 10 кПа) и относительно низкими скоростями движения (до 12 км/час).

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

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

Рис. 1. Схемы взаимодействия ЛПВП со льдом

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

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

Расчет НДС ледяной пластины проводим методом конечных элементов в пакете COSMOSM. В качестве расчетных характеристик материала приняты: модуль упругости Е = 5000 Мпа, коэффициент Пуассона μ = 0,34. Предел прочности льда на изгиб принят по рекомендациям В.В. Лаврова [4], где приведен график σн = f(k) при отрицательных температурах воздуха. Этот график с достаточной точностью можно аппроксимировать зависимостью

Рис. 2. Зависимость σ = f(h) по В.В. Лаврову [4]_________ и аппроксимация • • • • •

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

В зависимости (2):

Е, μ – модуль упругости и коэффициент Пуассона для льда;

h – толщина льда;

ρв – плотность воды;

g – ускорение свободного падения

Для выбора размеров конечного элемента был выполнен расчет для разных его размеров и проведено сравнение с аналитическим решением. Элементы разных размеров Δ «бесконечной» пластины нагружались равным по величине суммарным усилием Q и полученный прогиб под ним сравнивался с прогибом от сосредоточенной силы P той же величины, вычисленным по формуле Герца [6] для бесконечной пластины в начале координат.

Где D – цилиндрическая жесткость пластины при изгибе,

Схема нагружения показана на рис. 3

Рис. 3. Расчетная схема при определении размера элемента

При относительном размере элемента решение Герца и численное решение по предложенной методике, относительная погрешность не превышает 0,5%.

Бесконечная пластина, нагруженная равномерной нагрузкой в форме квадрата

Нагрузка равномерно распределена в форме квадрата со стороной а. Относительный размер нагрузки a/r варьировался от 0,1 до 3,5 с шагом 0,1 и от 3,5 до 5,0 с шагом 0,5. Изгибные напряжения принимались по рис. 2. В качестве расчетного принято давление равное Р = ρgh, Па (при h = 0,1 м). Безразмерные величины определялись по формулам:

где W – прогиб пластины.

График безразмерных напряжений   в ледяной пластине при постоянном давлении приведен на рис.4

Как видно из рис. 4, максимальные напряжения в рассмотренном случае возникают при безразмерной стороне квадрата a/r = 3 ÷ 3,5.


Рис. 4. График безразмерных напряжений в ледяной пластине при постоянном давлении, распределенном по площади квадрата со стороной а

Напряжения для пластины могут быть вычислены с помощью формул [6]

Изгибающий момент от сосредоточенной силы в цилиндрических координатах можно определить [6]

где:– постоянная Эйлера ();

R – текущий радиус;

P – сосредоточенное усилие;

– радиальный и тангенциальный моменты.

Подставив (6) в (5) можно видеть, что связь между напряжением и нагрузкой носит линейный характер для заданной толщины льда, формы и площади нагрузки.

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

Бесконечная пластина, нагруженная равномерной нагрузкой в форме прямоугольника

Нагрузка распределена в форме прямоугольника со сторонами а и в. Отношение сторон нагрузки а/в варьируется от 0,5 до 2,0 с шагом 0,1 и от 2,5 до 4,0 с шагом 0,5. Относительная площадь варьируется от 0,5 до 1,5 с шагом 0,25 и от 2,0 до 4,0 с шагом 0,5.

В качестве расчетной принята нагрузка . Результаты расчета приведены на рис.5.


Рис. 5. Безразмерные напряжения в бесконечной пластине при постоянной нагрузке, распределенной по площади прямоугольника

Полубесконечная пластина, нагруженная равномерной нагрузкой в форме прямоугольника

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

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

Анализ полученных данных показывает, что при относительно небольших площадях (до ) напряжения, вызываемые одинаковой по величине суммарной нагрузкой в 1,5 раза больше для полубесконечной пластины, чем для бесконечной. По мере роста площади нагрузки это отношение уменьшается.

Полубесконечная пластина с вырезом в форме канала

Нагрузка распределена в форме прямоугольника со сторонами а и в. Отношение сторон нагрузки варьировалось от 0,5 до 2,0 с шагом 0,1 и от 2,5 до 4,0 с шагом 0,5. Относительная ширина канала варьировалась от 1,0 до 4,0 с шагом 1,0. Расчетной принята нагрузка h=0,1 м.

Отношение напряжений в «полубесконечной» пластине с вырезом при постоянной нагрузке приведены в табл. 1

Таблица 1

Отношение напряжений в «полубесконечной» пластине с вырезом при постоянной суммарной нагрузке

Отношение a/b

 

Относительная ширина канала

 

1

2

3

4

В углу канала

На оси симметрии

В углу канала

На оси симметрии

В углу канала

На оси симметрии

В углу канала

На оси симметрии

0,5

0,808

0,735

0,544

0,397

0,358

0,214

0,225

0,109

0,6

0,738

0,694

0,470

0,363

0,296

0,189

0,179

0,093

0,7

0,676

0,656

0,407

0,331

0,247

0,167

0,146

0,080

0,8

0,620

0,620

0,354

0,303

0,208

0,149

0,122

0,070

0,9

0,570

0,587

0,309

0,278

0,178

0,133

0,103

0,061

1,0

0,525

0,556

0,271

0,256

0,153

0,119

0,090

0,054

1,1

0,485

0,528

0,239

0,236

0,134

0,108

0,080

0,049

1,2

0,448

0,501

0,211

0,218

0,119

0,098

0,072

0,044

1,3

0,415

0,476

0,188

0,202

0,107

0,090

0,066

0,041

1,4

0,385

0,453

0,168

0,188

0,097

0,083

0,060

0,038

1,5

0,357

0,432

0,152

0,175

0,089

0,077

0,056

0,035

1,6

0,332

0,412

0,137

0,164

0,083

0,072

0,053

0,033

1,7

0,309

0,393

0,125

0,154

0,077

0,067

0,050

0,031

1,8

0,288

0,376

0,115

0,145

0,072

0,064

0,047

0,029

1,9

0,269

0,359

0,107

0,137

0,068

0,060

0,045

0,028

2,0

0,251

0,344

0,099

0,129

0,065

0,057

0,042

0,026

2,5

0,183

0,280

0,074

0,102

0,052

0,045

0,034

0,021

3,0

0,138

0,234

0,061

0,085

0,043

0,038

0,028

0,018

3,5

0,109

0,199

0,052

0,072

0,037

0,033

0,024

0,015

4,0

0,090

0,173

0,046

0,063

0,033

0,028

0,021

0,013

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

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

Следовательно, для расчетного случая при определении НДС следует выбирать нагружение полубесконечной пластины, т.к. он требует больших усилий.

Оптимальное соотношение сторон нагрузки при проектировании ЛПВП целесообразно принимать равным 1.

Рецензенты:

Зуев В.А., д.т.н., профессор, заведующий кафедрой кораблестроения и авиационной техники Нижегородского государственного технического университета им. Р.Е. Алексеева, г. Нижний Новгород;

Грамузов Е.М., д.т.н., профессор кафедры кораблестроения и авиационной техники Нижегородского государственного технического университета им. Р.Е. Алексеева, г. Нижний Новгород.