ШИНЫ. НЕКОТОРЫЕ ПРОБЛЕМЫ ЭКСПЛУАТАЦИИ И ПРОИЗВОДСТВА

Методика и алгоритм математического моделирования процесса выделения вулканизационных газов на ПЭВМ

Блок-схема алгоритма математического моделирования процесса обработки на ПЭВМ представлена на рис. 56.

Алгоритм имеет следующую структуру. Вначале вводятся исходные данные, константы, величина шага по времени при решении конечно-разностных уравнений. Затем следует под­программа инициализации графического режима. Перечислен­ные блоки составляют общую часть программы. Далее следует циклическая часть, которая многократно повторяется в процессе расчета. Она начинается с оператора-счетчика: К=К+1. Затем идет собственно расчет искомых величин. Расчет температур жидкой и газовой фаз в камере для двух рассматриваемых слу­чаев осуществляется по одним и тем же конечно-разностным формулам. После расчета температур и вывода результата на печатающее устройство алгоритм разветвляется на два направ­ления по признаку ветвления «X».

Значение переменной Х=1 соответствует процессам, про­текающим с бурным выделением газов. Значение Х==2 соот­ветствует процессам, протекающим с выделением небольшого количества газов. Значение X задается в блоке задания исход­ных данных. В первом случае вычисляются величины: парци­альных давлений, парциальных плотностей, парциальных по­токов и мольных долей жидких компонентов, входящих в со­став обрабатываемого изделия. Во втором случае вычисляются величины: парциальных давлений, парциальных концентраций, парциальных потоков и мольных долей жидких компонентов, входящих в состав обрабатываемого изделия. Каждый блок рас­чета заканчивается выводом результата на печать.

Далее следует общая для обоих случаев подпрограмма по­строения графиков. Затем проверяется условие достижения за­данного предела по времени Ктах. Если оно выполняется, то расчет заканчивается. Если нет, то уравнение передается к на­чалу циклической части алгоритма и осуществляется очеред­ной расчетный шаг.

Временной шаг выбирается в начале как 0,0001 доля от дли­тельности всего процесса в целом. Если такой шаг оказывается недостаточно мелким и приводит к численной неустойчивос­ти, то следует повторно запустить программу с шагом, умень­шенным в 2 раза. Эту операцию следует повторять до тех пор, пока решение не приобретет устойчивый характер.

Среднее время вычислений составляет 4-6 с на компьюте­ре РС АТ 80486.