Алгоритмическая и программная реализация.
Перед изложением алгоритма решения укажем какими входными параметрами мы оперируем. Входными параметрами здесь являются: линейный размер исследуемого участка акватории моря LA и сплоченность ледяного поля NA; интервал времени функционирования системы Ts; радиус rc и форма опоры сооружения; удельная энергия разрушения льда E0; плотность r, прочность Rice и толщина льда h; скорость движения льдины на свободной воде v и ее радиус rice. Прочность льда на смятие Rb оценивается согласно правилам [1]. Легко заметить, что начальное расстояние между льдинами l0 определяется выражением
(12)
Предполагается, что до опоры гидротехнического сооружения льдины следуют одна за другой на одинаковом расстоянии l0 между собой. Если перед опорой уже имеются льдины, то расстояние l(t) между крайней льдиной затора и вновь приближающейся льдиной определяется выражением
, (13)
где T - время окончания процесса прорезания и l(0)=l0. Случай T=0 соответствует отсутствию внедрения: а) либо кинетическая энергия льдины не достигла критического значения Ec, б) либо потенциальная энергия остановившегося ледяного поля недостаточна для реализации процесса прорезания. В этом случае оцениваем время t1 за которое приближающая льдина достигнет крайней в скопившемся перед опорой и остановившемся ледяном поле
, (14)
Далее системное время модели увеличивается на t1, число льдин в системе увеличивается на единицу, l(t) полагается равным l0 и цикл повторяется. Если l(t)³0, то в систему дополнительных льдин не поступает; в противном случае число льдин увеличивается на единицу.
Алгоритм функционирования системы «лед-опора ледостойкого основания».
П1. Открываем файл входных значений и считываем эти значения (OpenFilePar0)
p.s_Region –
p.R_ICE –
p.V_ICE -
p.D_ICE –
p.H_ICE –
p.Res_ICE –
p.N_ICE –
v.dt –
p.R_Con –
p.k_Con –
p.X_Con –
p.k1_Con –
p.T_s –
П2. Загрузка диалога (Load frmTask0)
П3. Вызов основной процедуры (Calculate)
а) процедура формирования промежуточных входных значений – первоначальное расстояние между льдинами (InitResult0)
б) печать входных значений (RezultRecFile)
Print p.s_Region – печать размер региона;
Print p.k_Con – печать коэффициент формы опоры;
Print p.D_ICE – печать плотность льда;
Print p.R_ICE – печать прочность льда;
Print p.H_ICE – печать толщина льда;
Print p.N_ICE – печать сплоченность льда;
Print p.T_s – печать время сезона;
Print p.L0 – расстояние между льдинами;
Print p.Ar0 – печать площадь льдины в) организация циклов расчета по времени и сам расчет:
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.