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

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 
170 HN = H1 – C * X: D1 = D – 2 * HN 
190 D2 = D – 2 * H2: 200 HSR = (H2 + HN) / 2 
210 S1 = (PI ^ 2 * (D – HSR) * (D – 2 * HSR) * N) / (60 * HSR * SQR(PI ^ 2) * (D – 2 * HSR) ^ 2 + (T ^ 2)) 
220 J1 = ((PI – 2) * (D – 2) – (4 * (T * 2))) / (PI – 2) 
230 J2 = ((D + D2) – 3 – (D + DL) – 3) / (3 * (D2 – DD)) 
240 J3 = ((2.3 * (PI ^ 2) * (D – 5)) * (LOG(HN / H2)) / 2.3) / (((T ^ 2) + (PI ^ 2) * (D ^ 2)) * (HN – H2)): 250 J = 
J1 + J2 + J3 
260 PRINT "СКОРОСТЬ СДВИГА S1 = "; S1 
270 INPUT "ВВЕДИТЕ ЭФФЕКТИВНУЮ ВЯЗКОСТЬ ml = "; M1 
280 N11 = ((.0001) * (PI ^ 3) * (T – E) * LDOZ * J * M1 * (N – 2)) / (36 * T) 
290 N12 = ((.001) * A2 * DP * N) / 6: N1 = N11 + N12 

310 PRINT "МОЩНОСТЬ N1 = "; N1 
320 S2 = ((PI^2)*(D ^ 2) * N) / (60 * DL * SQR((PI ^ 2) * (D ^ 2) – (T ^ 2))) 
330 PRINT "СКОРОСТЬ СДВИГА S2 = "; S2 
340 INPUT "ВВЕДИТЕ ЭФФЕКТИВНУЮ ВЯЗКОСТЬ m2 = "; M2 
350 N2 = ((PI^3)*(D^3) * E * LDOZ * M2 * .0001 * (N ^ 2)) / (36 * DL * T) 
360 PRINT "МОЩНОСТЬ N2 = "; N2 
370 PRINT "ТЕХНОЛОГИЧЕСКАЯ МОЩНОСТЬ NT = "; N1 + N2 
380 END 

ПРОГРАММА 3 

 

Программа расчета технологической мощности и производительности двухшнекового смеси-

теля типа СН. Язык Бейсик. 

 

REM "Программа расчета технологической мощности и производительности" 
REM "      двухшнекового смесителя типа СН" 
 
OPEN "SNDPROB.res" FOR OUTPUT AS #1 
REM "Введите следующие параметры в системе измерения [СИ]" 
'FI = X#(1): HG = X#(2): e1 = X#(3): nw1 = X#(4): ek1 = X#(5): DG = X#(6) 
X#(1) = 4: X#(2) = .002: X#(3) = .001: X#(4) = 10 / 60: X#(5) = .008: X#(6) = .04 
' "Плотность набухшей массы клеевой композици, кг/м^3" 
RO = 850 
DG = X#(6) 
' "Вязкость перерабатываемого материала, Па*с" 
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) 
' "Индекс качества смеси, который задается и находится в интервале (0,85 - 0,95)" 
Jk1 = .95 
' "Расчет коэффициента сдвига методом последовательных приближений, об" 
 tauv1 = 10'- начальное приближение по времени, с 
15 GAM = tauv1 * X#(4) 
' "Коэффициент изменения качества смеси за один оборот шнеков" 
Nuu1 = -(-.06 * EXP(-.16 - .06 / 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) 
' "Сумма мощностей, затрачиваемых на отдельных участках шнеков"