Инженерная оптимизация смесительного и валкового оборудования. Клинков А.С - 45 стр.

UptoLike

Рубрика: 

ПРОГРАММА 2
REM "имя файла - SNOPD01.bas"
REM "Минимизация целевой функции методом скользящего допуска"
REM "Минимизация технологической мощности двухшнекового смесителя при заданной его производительности "
DECLARE SUB FEAS ()
DECLARE SUB STAR ()
DECLARE SUB WRIT ()
DECLARE SUB SUMR ()
DECLARE SUB prob (INQ!)
DIM SHARED X#(50), X1#(50, 50), X2#(50, 50), R#(100), R1#(100), R2#(100), R3#(100), SUM#(50), F#(50), SR#(50),
ROLD#(100), H#(50), FLG#(10), A#(50, 50)
COMMON SHARED kik, NX, NC, NIC, STEP0#, ALFA#, BETA#, GAMA#, IN, INF, FDIFER#, SEQL#, k1, k2, k3, k4, k5,
k6, K7, K8, K9, FOLD#, L5, L6, L7, L8, L9, R1A#, R2A#, R3A#
OPEN "SNOpD01.res" FOR OUTPUT AS #1
PRINT #1, " Расчет технологической мощности двухшнекового смесителя"
PRINT " Расчет технологической мощности двухшнекового смесителя"
PRINT "-----------------------------------------------------------"
PRINT #1, " ------------------------------------------------------------"
FDIFER# = 1#
NX = 6 'общее число переменных
NC = 1 'общее число ограничений в виде равенств
NIC = 12 'общее число ограничений в виде неравенств
SIZE# = 5# 'опред. размер деформ. многогранника
CONVER# = .01# 'используется для окончания поиска
ALFA# = 1#
BETA# = .5#
GAMA# = 2#
10 STEP0# = SIZE#
'введите начальные значения варьируемых параметров в системе, [СИ]"
X#(1) = 4# ' FI координаты
X#(2) = .002# ' HG стартовой
X#(3) = .001# ' e1 точки
X#(4) = 10 / 60# ' nw1
X#(5) = .004# ' ek1
X#(6) = .05# ' DG1
'CALL prob(3)
' "ДЛЯ ВВОДА ИСХОДНЫХ ДАННЫХ НАЖМИТЕ КНОПКУ F2"
IF FDIFER# < CONVER# THEN GOTO 9999
PRINT " <<<<< *HAЧAЛO BЫЧИCЛEHИЙ* >>>>>"
PRINT #1, " <<<<< *HAЧAЛO BЫЧИCЛEHИЙ* >>>>>"
PRINT "-----------------------------------------------------------"
Далее с помощью кнопки “F2” войти в подпрограмму SUB PROB (INQ) и ввести с клавиатуры в тексте
программы значения параметров.
Затем запустить программу на счёт (F5)
SUB prob (INQ)
' "Заданная производительность, кг/ч"
QZ = 2
' "Плотность набухшей массы клеевой композиции, кг/м^3"
RO = 850
' "Диаметр рабочих органов (шнеков), м"
'DG = .1
DG = X#(6)
'FI = X#(1): HG = X#(2): e1 = X#(3): nw1 = X#(4): ek1 = X#(5): DG = X#(6)
FOR nno = 1 TO 6
IF X#(nno) < 0 THEN X#(nno) = -1 * X#(nno)
NEXT nno
' "Вязкость перерабатываемого материала, Па*с"
ETA1 = 27500: pm1 = .5
ETA2 = 27500: pm2 = .5
ETA3 = 27500: pm3 = .5
ETA4 = 27500: pm4 = .5
' "Заданный перепад давления, [Па] по длине рабочих органов и количество пар смесительных дисков"
PG1 = 5000000: nk1 = 45
' "Cредний коэффициент заполнения дисков на всей длине дисковой зоны"