содержание .. 42 43 44 45 ..

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

' "Для смесительных дисков, имеющих эксцентричную форму"
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

Mu4 = ETA4 * (gradV4) ^ (pm4 - 1)
tau4 = Mu4 * gradV4
NT4 = 2.75E-07 * 3.14 ^ 2 * DGs * AA * X#(4) * X#(5) * (ALF * 180 / 3.14) * tau4 * nk1
' "на выдавливание материала в зонах обратных витков и формующей головки"
NT5 = 2.75E-07 * PG1 * (DG ^ 2 - dGv ^ 2) * 3.14 * DGs * X#(4) * t1 * (ALF * 180 / 3.14)
' "Сумма мощностей, затрачиваемых на отдельных участках шнеков"
R#(14) = (NT1 + NT2 + NT3) * ZV + NT4 + NT5
NTS = R#(14)
'PRINT NT1; NT2; NT3; NT4; NT5
IF INQ = 2 THEN GOTO 100
IF INQ = 3 THEN GOTO 200
'ограничениe в виде равенств
R#(1) = (1 - Q / QZ)
CLS
LOCATE 2, 1: PRINT "Ждите! Идет процесс оптимизации"
LOCATE 4, 1: PRINT "Q="; Q; " [кг/ч]"; " P="; PG1; " [Па]"
LOCATE 5, 1: PRINT "tauv="; tauv; " [с]"; " LL="; LL; "[м]"
LOCATE 6, 1: PRINT "NTS="; NTS; " [кВт] "; "Nuu1="; Nuu1; "[дол. ед]"
LOCATE 7, 1: PRINT "FI - угол наклона винтовой нарезки шнека "; X#(1); "[град]"
LOCATE 8, 1: PRINT "h1- глубина винтовой нарезки шнека "; X#(2); "[м]"
LOCATE 9, 1: PRINT "DG - наружный диаметр шнека "; DG; "[м]"
LOCATE 10, 1: PRINT "nw1 - частота вращения шнека "; X#(4) * 6.28; "[рад/с]"
LOCATE 11, 1: PRINT "e1 - ширина гребня витка шнека "; X#(3); "[м]"
LOCATE 12, 1: PRINT "ek1 - ширина смесительного диска (кулачка)"; X#(5); "[м]"

LOCATE 2, 1: PRINT #1, "Ждите! Идет процесс оптимизации"
LOCATE 4, 1: PRINT #1, "Q="; Q; " [кг/ч]"; " P="; PG1; " [Па]"
LOCATE 5, 1: PRINT #1, "tauv="; tauv; " [с]"; " LL="; LL; "[м]"
LOCATE 6, 1: PRINT #1, "NTS="; NTS; " [кВт] "; "Nuu1="; Nuu1; "[дол. ед]"
LOCATE 7, 1: PRINT #1, "FI - угол наклона винтовой нарезки шнека "; X#(1); "[град]"
LOCATE 8, 1: PRINT #1, "h1- глубина винтовой нарезки шнека "; X#(2); "[м]"
LOCATE 9, 1: PRINT #1, "DG - наружный диаметр шнека "; DG; "[м]"
LOCATE 10, 1: PRINT #1, "nw1 - частота вращения шнека "; X#(4) * 6.28; "[рад/с]"
LOCATE 11, 1: PRINT #1, "e1 - ширина гребня витка шнека "; X#(3); "[м]"
LOCATE 12, 1: PRINT #1, "ek1 - ширина смесительного диска (кулачка)"; X#(5); "[м]"

GOTO 300
'ограничения в виде неравенств
100
R#(2) = X#(1) / 4 - 1
R#(3) = 1 - X#(1) / 8
R#(4) = X#(2) / (.003) - 1
R#(5) = 1 - X#(2) / (.015)
R#(6) = X#(3) / (.001) - 1
R#(7) = 1 - X#(3) / (.01)
R#(8) = X#(4) / (10 / 60) - 1
R#(9) = 1 - X#(4) / (100 / 60)
R#(10) = X#(5) / (.002) - 1
R#(11) = 1 - X#(5) / (.015)
R#(12) = DG / (.04) - 1
R#(13) = 1 - DG / (.1)
GOTO 300
'целевая функция (технологическая мощность)
200
R#(14) = NT1 + NT2 + NT3 + NT4 + NT5

300 X#(6) = DG':X#(1) = lam: X#(2) = HG: X#(3) = e1: X#(4) = nw1: X#(5) = ek1
'PRINT #1, "tauv="; tauv; "c"
'PRINT #1, "Q="; Q; "[м^3/ч]"
'PRINT #1, "LL="; LL; "[м]"
END SUB

ПРОГРАММА 1

ПРОГРАММА РАСЧЕТА МАКСИМАЛЬНОЙ ПРОИЗВОДИТЕЛЬНОСТИ ШНЕКОВОЙ

МАШИНЫ ПРИ ОПРЕДЕЛЕННОМ СОПРОТИВЛЕНИИ ФОРМУЮЩЕГО ИНСТРУМЕНТА И

ВАРЬИРОВАНИЕМ ГЛУБИНЫ НАРЕЗКИ ШНЕКА. ЯЗЫК БЕЙСИК

30 CLS
40 INPUT "ВЕЛИЧИНА РАДИАЛЬНОГО ЗАЗОРА DL, [м]"; DL
50 INPUT "КОЭФ-Т ГЕОМЕТРИЧЕСКОЙ ФОРМЫ ГОЛОВКИ, [м] К = "; K
60 INPUT "СКОРОСТЬ ВРАШЕНИЯ ШНЕКА, [об/c] N = "; N
70 INPUT "ВЫБРАННАЯ СТЕПЕНЬ СЖАТИЯ МАТЕРИАЛА I = "; W
80 INPUT "ДЛИНА ШНЕКА, [м] L = "; L
90 INPUT "ДИАМЕТР ШНЕКА, [м] D = "; D
100 INPUT "ГЛУБИНА ВИНТОВОГО КАНАЛА В ЗОНЕ ЗАГРУЗКИ ИЗМЕНЯЕТСЯ ОТ, [м]

HN = "; HN

110 INPUT ", [м] ДО HK = "; HK
120 INPUT "ШИРИНА ГРЕБНЯ ВИТКА ИЗМЕНЯЕТСЯ ОТ, [м] EN-"; EN
130 INPUT ", [м] ДО ЕК = "; EK
140 PI = 3.141592
150 FO = ATN(1 / PI)
160 R = D / 2
170 T = D
180 DL = DL * D
190 HN = HN * D
200 HK = HK * D
210 EN = EN * D
220 EK = EK * D
230 QMAX = 1E-10
240 FOR H1 = HN TO HK STEP .01 * D
250 H2 = .5 * (D – SQR(D – 2 – (4 * H1 * (D – H1)) / W))
260 FOR E = EN TO EK STEP .01 * D
270 J6 = (L / (H1 – H2)) * LOG(H1 / H2)
280 J7 = L / (H1 * H2)
290 J8 = L / 2 * ((HbH2) / (Hl – 2 * H2 – 2))
300 J9 = (L / (H1 – H2)) * LOG((2 * R – H2) / (2 * R – HD))
310 S = ((2 * R * J7 – J6) / (2 * R * TAN(FO)))

320 F = ((2 * R * J8 – COS(2 * FO) * J7) / (R * SIN(2 * FO))) + (TAN(FO)
* (J6 – J9)) / 4 * R – 2

330 B = (T – E) / (12 * F)
340 A = P1 * O * (T – E) * 8 / (2 * P)
350 F1# = ATN(T / (2 * PI * (R – H1)))
360 F2# = ATN(T / (2 * PI * (R – H2)))

370 J = 3 * 3.14 ^ 2 * D ^ 2 * DL – 3
380 I = (SIN(FI#)) ^ 3 * (SIN(F2#)) ^ 3 * (H1 – H2)
390 U = I * J
400 Z# = SIN(F1#)
410 Z2# = SIN(F2#)
420 'PRINT Z#, Z2#, F1#, F2#, H1, H2
430 V = 10 * L * T * E * (Z# – 3 – Z2# ^ 3)
440 G = U / V
450 Q = A * K * N / (K + B + G)
460 IF Q < QMAX THEN 480
470 QMAX = Q: H11 = H1: H21 = H2: A1 = A: B1 = B: G1 = G: E1 = E
480 NEXT E
490 NEXT H1
500 PRINT "МАКСИМАЛЬНАЯ ПРОИЗВОДИТЕЛЬНОСТЬ = "; QMAX; " [м^З/с]"
510 PRINT "ПРИ:"
520 PRINT "ГЛУБИНЕ ВИНТОВОГО КАНАЛА"
530 PRINT "В ЗОНЕ ЗАГРУЗКИ-"; H11; "[м]"
540 PRINT "В ЗОНЕ ДОЗИРОВАНИЯ-"; H21; " [м]"
550 PRINT "И ПРИ ПОСТОЯННЫХ"
560 PRINT "ПРЯМОГО ПОТОКА"; " ;A1;"; "(м^3);"
570 PRINT "ОБРАТНОГО ПОТОКА = "; B1; " [м^3]"
580 PRINT "ПОТОКА УТЕЧКИ-"; G1; " [м^3]"
590 PRINT "ШИРИНЕ ГРЕБНЯ ВИТКА ШНЕКА = "; E1; "[м]"
600 END

ПРОГРАММА 2

ПРОГРАММА РАСЧЕТА ТЕХНОЛОГИЧЕСКОЙ МОЩНОСТИ ШНЕКОВОЙ МАШИНЫ С

ПЕРЕМЕННОЙ ГЛУБИНОЙ НАРЕЗКИ ШНЕКА. ЯЗЫК БЕЙСИК

10 REM ПРОГРАММА РАСЧЕТА ТЕХНОЛОГИЧЕСКОЙ МОШНОСТИ
20 PI = 3.14
30 INPUT "ГЛУБИНА ВИНТОВОГО КАНАЛА ШНЕКА В 3ОНЕ ЗАГРУЗКИ, [м] H1 = "; H1
40 INPUT "В НАЧАЛЕ ЗОНЫ ДОЗИРОВАНИЯ, [м] HN'- = "; HN
50 INPUT "НА КОНЦЕ ШНЕКА, [м] Н2 = "; H2
60 INPUT "ШАГ ВИНТОВОЙ ЛИНИИ, [м] 1-"; T
70 INPUT "ДИАМЕТР ШНЕКА, [м] D = "; D
80 INPUT "СКОРОСТЬ ВРАШЕНИЯ ШНЕКА, [об/с] N = "; N
90 INPUT "ПЕРЕПАД ДАВЛЕНИЙ, [Па] DP = "; DP
100 INPUT "ШИРИНА ГРЕБНЯ НАРЕЗКИ, [м] Е = "; E
110 INPUT "ДЛИНА ЗОНЫ ДОЗИРОВАНИЯ, [м] Ldoz = "; LDOZ
120 INPUT "ДЛИНА ШНЕКА, [м] L"; L
130 INPUT "ПОСТОЯННАЯ ПРЯМОГО ПОТОКА, [м^3] А = "; A2
140 D1 = .003 * D: X = L – LDOZ: C = (H1 – H2) / L

содержание .. 42 43 44 45 ..