Многоцикловое и истирающее воздействия дрейфующего ледяного покрова на морские гидротехнические сооружения (Расчет конструкции ледостойких оснований. Особенности имитационного моделирования), страница 13

-  Время взыимодействия (T_InterAction)

Если v.v>0 да v.Tf=v.Fc/p.k1_Con/v.v; нет v.Tf=0 и рассматриваем либо процесс движения льда (Proc_MoveIce), либо процесс накопления льда (Proc_AddIce)

Если v.dt>0 да r.v=v.dx/v.dt; нет r.v=0 – расчетное (промежуточное) значение скорости

Если v.Tf>0 да v.Nf=v.t/v.Tf; нет v.Tf=0

Переход на

-  Печать промежуточных результатов и вычисление текущих переменных (OutResData)

Если v.dx=0 да Exit Sub нет

r.Tn=r.Tn+v.dt

r.n=r.n+1

Если r.v_Min>r.v да r.v_Min=r.v

Если r.v_Max<r.v да r.v_Max=r.v

Если r.dtMin>v.dt да r.dtMin=v.dt

Если r>.dtMax<v.dt да r.dtMax=v.dt

r.dtM=r.dtM+v.dt

Если r.dxMin>v.dx да r.dxMin=v.dx

Если r.dxMax<v.dx да r.dxMax=v.dx

r.dxM=r.dtM+v.dx

Если v.Fs=0 да Exit Sub

Если r.FsMin>v.Fs да r.FsMin=v.Fs

Если r.FsMax<v.Fs да r.FsMax=v.Fs

r.FsM=R.FsM+v.Fs

- Печать конечных результатов (EndShowResult)

Print T_end

Print

П5. К процедуре T_InterAction

- Процесс взаимодействия льда с опорой (Proc_MoveIce) вызывает функцию расета взаимодействия льда с сооружением П6.

- Процесс взаимодействия новых блоков льда (Proc_AddIce)

П6. Взаимодействие льда с сооружением (ContactIceConstr) вызывает:

П7. функцию расчета силы удара (FunForce)

и

П8. модель скорости движения льда (ModelVelosity)

Задаются входные параметры и начальные значения переменных; задаемся определенной моделью скорости и конкретным значением k; определяем массу отдельной льдины M₀=rhpr²_ice; вычисляем вектор A_i=A_c (x_i), где x_i=i×r_c/N (i=1¸N) и N заранее заданное целое число; полагается равными нулю число льдин пришедших в систему n_1i и ущедших из нее n_2i; индикатор текущего времени t полагается равным нулю.

П2. Если перед опорой имеет место свободная вода, то оцениваем глубину внедрения X и время прорезания T из соотношений (4) или (5). В противном случае, согласно предложенной выше мтодике, глубину внедрения оцениваем как корень уравнения (8) или (9).Если X > 0, то по (10) или (11) оцениваем T.

П3. Переходим к пункту П4 если X=0. Увеличиваем индикатор текущего времени t на T (t=t+T) и согласно (13) оцениваем l(t). Если линейный размер ближайшей к опоре льдины не превосходит определенное число диаметров сооружения, то число ушедших из системы льдин n_2i увеличивается на единицу. При l(t)>0 переходим к П5, в другом случае полагаем l(t) равным l₀ и число пришедших в систему льдин n_1i увеличиваем на единицу.

П4. Согласно (14) оцениваем t₁, число льдин в системе n_1i увеличиваем на единицу  и временной индикатор увеличиваем на t₁.

П5. Если t>T_s, то окончание работы, в другом случае - пересчитываем массу блока льдин перед опорой, выводим промежуточные результаты и уходим на П2.

П6. Окончание работы.

Программная реализация системы представляет собой открытый интерактивный комплекс и выполнен в рамках одного из самых мощных приложений Windows – программной оболочки Visual Basic 5. Посредством упорядоченного набора закладок(программных средств непосредственно самой оболочки) выполняется настройка и выбор математической модели, численных значений параметров и режимов функционирования системы. Поясним смысл вышесказанного.

При запуске системы посредством мыши предлагается выбрать ту или закладку. В каждой из них присутствует и/или модифицируется определенная информация об организации процесса вычислений, численных значений параметров и моделей. Так закладка «ДАННЫЕ» содержит численные значения входных параметров. Параметры, относящиеся к одному объекту упорядочены посредством их заключения в соответствующие контейнеры. Например, параметры для габаритов сооружения упакованы в контейнер «Параметры сооружения»,а маршруты файла параметров и файла результатов - в контейнер  «Ввод/Вывод».Модификация того или иного параметра.

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