Автоматизированное проектирование и расчет шнековых машин - часть 50

Fmax03 = 1 / K1 ^ 2 * (RO1 * F1 * L / K1 - A1) * SIN(K1 * L)

Fmax = (Fmax01 - Fmax02 - Fmax03)/(E*J1)

IF Fmax1 < .0000001 THEN Fmax1 = 0
IF Fmax2 < .0000001 THEN Fmax2 = 0
PRINT "SIGekv="; SIGekv; SIG
99 IF INQ = 2 THEN GOTO 100
IF INQ = 3 THEN GOTO 200
‘ Ограничение в виде равенства

R#(1) = (1 - Q# / Q) * 1
R#(2) = (1 - (TD# - T0) / delT)
PRINT "Q#="; Q#; "dT="; (TD# - T0)

PRINT "N="; R#(15); "JB="; JB GOTO 300
100 ‘Интервалы изменения варьируемых параметров (ограничения в виде неравенств)
R#(3) = 1 - JB / JB1# ‘Левая граница
R#(4) = X#(1) - 16 ‘Левая граница

R#(5) = 21# - X#(1) ‘Правая граница

R#(6) = (X#(2) - .003#) * 10 ‘Левая граница
R#(7) = (7.000000000000001D-02 - X#(2)) * 100 ‘Правая граница
R#(8) = (X#(3) - 2.000000000000001D-02) * 10 ‘Левая граница
R#(9) = (9.000000000000001D-02 - X#(3)) * 10 ‘Правая граница
R#(10) = X#(4) - .1# ‘Левая граница
R#(11) = 1.5# - X#(4) ‘Правая граница
R#(12) = X#(5) - .2# ‘Левая граница
R#(13) = 1# - X#(5) ‘Правая граница
'PRINT "SIGekv="; SIGekv; SIG
R#(14) = 1 - SIGekv / SIG
GOTO 300
200 ‘Расчёт технологической мощности (целевой функции)
CS = COS(X#(1) * 3.1416 / 180)
SN = SIN(X#(1) * 3.1416 / 180)
NN1 = MU * W1 * (3.14 * X#(3) * X#(4)) ^ (1 + pm) / X#(2) ^ pm * (4 * SN ^ 2 + CS ^ 2) * X#(5)

* FZ / SN

NN2 = 3.14 * X#(3) * X#(4) * CS / 2 * W1 * X#(2) * dP

NN3 = (3.14 * X#(3) * X#(4)) ^ (pm + 1) * MU * E1 * X#(5) / TAN(X#(1) * 3.1416 / 180) / d1 ^ pm

* FZ

R#(15) = NN1 + NN2 + NN3
300 END SUB

ПРОГРАММА 10

Расчёт оптимальных технологических и конструктивных параметров двухшнекового смесите-

ля

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редний коэффициент заполнения дисков на всей длине дисковой зоны"
KKD = .8
' "Cредний коэффициент заполнения нарезной части сборного шнека"
KKN = .866
' "число витков в прямых насадках"
ZV = 10
'21 "Расчет среднего диаметра и межосевого расстояния"
DGs = DG - X#(2): AA = DGs
' "Расчет внутреннего диаметра шнеков (диаметра сердечника)"
dGv = DG - 2 * X#(2)
' "Расчет угла захвата"
ALF = 2 * (3.14 / 2 - ATN(AA / SQR(ABS(DG ^ 2 - AA ^ 2))))
' "Коэффициент уменьшения объема полости смешения за счет сопряжения шнеков"
cc1 = AA / DG
'PRINT "ДЛЯ НАЧАЛА ВЫЧИСЛЕНИЙ НАЖМИТЕ F5"
' "Расчет шага винтовой нарезки"
t1 = 3.14 * DG * TAN(X#(1) * 3.14 / 180): SS = t1
' "Расчет длины участка сборного шнека, занятая дисками"
LD = X#(5) * nk1
' "Расчет длины нарезной части сборного шнека"
LN = ZV * t1
' "Расчет общей сборного шнека"
LL = LD + LN
' "Расчет суммарного среднего коэффициента заполнения полости смешения"
KK1 = KKD * LD / LL + KKN * LN / LL
' "Расчет суммарного объема рабочей полости смешения"
' "Объем полости одного витка нарезки отдельно взятого шага"
VS = SS / 2 * (DG - X#(2)) * X#(2) * 3.14
' "объем полости одного витка нарезки отдельно взятого шага с учетом зазоров между шнеками в заце-
плении"
VSZ = (e1 * VS) / (e1 - .01 * DG)
' "Объем зазора на длине отдельно взятого шага между шнеком и корпусом"
VZ = .01 * 3.14 * DG ^ 2 * SS
' "полный действующий объем одного витка с учетом всех зазоров для отдельно взятого шага нарезки"
VZD = VZ / 2 + VSZ
' "объем рабочей полости смешения для насадок с одинаковым шагом нарезки"
VCN = VZD * ZV
' "Для смесительных дисков (кулачков), имеющих трехгранную форму"
'VCD = 3 * (LOG(AA * DG / 4 - .5238 * (AA ^ 2 - DG ^ 2 / 2))) / 2.3

' "Для смесительных дисков, имеющих эксцентричную форму"
VCD = 3.14 / 4 * (DG ^ 2 - DD ^ 2) * LD
' "Расчет суммарного объема рабочей полости смешения"
VC = 2 * cc1 * (VCN + VCD)
tauv2 = 3600 * KK1 * VC * RO / QZ
IF tauv2 < 0 THEN tauv2 = -1 * tauv2
' "Индекс качества смеси, который задается и находится в интервале (0,85 - 0,95)"
Jk1 = .95
' "Расчет коэффициента сдвига методом последовательных приближений, об"
tauv1 = 10'tauv2 / 2'- начальное приближение по времени, с
15 GAM = tauv1 * X#(4)
' "Коэффициент изменения качества смеси за один оборот шнеков"
Nuu1 = -(-.06 * EXP(-.16 - .06 / 1 * GAM))
' "Расчет среднего времени нахождения смеси в полости смешения"
IF Nuu1 <= 0 THEN Nuu1 = .0001
tauv = -LOG(1 - Jk1) / (Nuu1 * X#(4))
IF tauv < 0 THEN tauv = -1 * tauv
IF tauv > 1000 THEN tauv = 1000
'PRINT "tauv1="; tauv1 / 60; "[мин] "; "tauv="; tauv / 60; "[мин] "; "tauv2="; tauv2 / 60; "[мин] "
IF ABS((tauv - tauv1) / tauv) * 100 < 2 THEN 21 ELSE 19
19 tauv1 = tauv: GOTO 15
' "Расчет объемной производительности двухшнекового смесителя"
21 Q = 3600 * KK1 * VC * RO / tauv / 2

' "Расчет потребляемой мощности"
' "между корпусом и шнеками"
det1 = .01 * DG
gradV1 = 3.14 * DG * X#(4) / det1
IF gradV1 < 0 THEN gradV1 = -1 * gradV1
Mu1 = ETA1 * (gradV1) ^ (pm1 - 1)
tau1 = Mu1 * gradV1
NT1 = (.002 * (1 + ALF / (2 * 3.14)) * X#(3) * 3.14 ^ 2 * DG ^ 2 * X#(4) * tau1) / COS(X#(1) * 3.14 / 180)
' "в зазорах между боковыми поверхностями шнеков"
det2s = (det1 + t1 - 2 * X#(3)) / 2
gradV2 = 2 * 3.14 * DG * X#(4) / det2s
IF gradV2 < 0 THEN gradV2 = -1 * gradV2
Mu2 = ETA2 * (gradV2) ^ (pm2 - 1)
tau2 = Mu2 * gradV2
NT2 = .000033 * (3.14 * DG ^ 3 * SIN(ALF / 2) * (1 - COS(ALF / 2)) * X#(4) * tau2) / COS(X#(1) * 3.14 /
180)
' "в зазоре между вершиной и впадиной витков"
det2 = (DG - dGv) / (2 * COS(ALF / 2))
det3s = (det1 + det2) / 2
gradV3 = 3.14 * DG * X#(4) / det3s
IF gradV3 < 0 THEN gradV3 = -1 * gradV3
Mu3 = ETA3 * (gradV3) ^ (pm3 - 1)
tau3 = Mu3 * gradV3
NT3 = 5.4E-07 * (3.14 ^ 2 * X#(4) * X#(3) * DGs * (ALF * 180 / 3.14) * (DG - dGv) * tau3) / COS(X#(1) *
3.14 / 180)
' "между кулачками одной пары и корпусом"
det4 = .01 * DG
gradV4 = 3.14 * DG * X#(4) / det4
IF gradV4 < 0 THEN gradV4 = -1 * gradV4

содержание .. 48 49 50 51 ..