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

UptoLike

61
PRINT " <<<<< *HAЧAЛO BЫЧИCЛEHИЙ* >>>>>"
PRINT #1, " <<<<< *HAЧAЛO BЫЧИCЛEHИЙ* >>>>>"
Далее с помощью кнопки “F2” войти в подпрограмму SUB PROB
(INQ) и ввести с клавиатуры в тексте программы соответствующие пара-
метры.
Затем запустить программу на счет (F5)
SUB prob (INQ)
' "ввод перепада давления по длине шнека,[Па]"
nnn = 2' число витков (от начала винтовой нарезки до начала разрыва)
nk = 10 'число витков (суммарное)
PP = 7000000: dP = PP / nk
' "ввод геометрических параметров [м]"
LL0 = .016: R1 = .016: tt = .032: LL = tt * nk: BE = (17 * 3.14 / 180)
' "ввод допускаемого напряжения [Па] и коэффиц. трения"
SIGD = 325000000: FF = .2
' "Расчет допускаемого прогиба [м]"
WD = .01 * 2 * R1
' "ввод плотности [кг/м^3], модуля упругости [Па] и ускорения свобод.
паден. [м/с^2]"
RO = 7850: EE = 2E+11: gg = 9.81
' "расчет распределенной нагрузкии, сечения шнека, момента инерции и
момента временного сопротивления"
q = RO * gg * 3.14 / LL * ((R1 - X#(1)) ^ 2 * LL + 2 * (R1 - X#(1) / 2) /
COS(BE) * X#(1) * X#(2) * nk - X#(3) ^ 2 * LL)
FFP = 3.14 * ((R1 - X#(1)) ^ 2 - (X#(3)) ^ 2)
J2 = 3.14 * (R1 - X#(1)) ^ 4 / 4 * (1 - (X#(3) / (R1 - X#(1))) ^ 4)
WW0 = 3.14 * (R1 - X#(1)) ^ 3 / 2 * (1 - (X#(3) / (R1 - X#(1))) ^ 4)
'LL1 = (2 * 3.14 * (R1 - X#(1)) ^ 4 * WD * EE / q) ^ .25
PRINT "q="; q; "LL="; LL; "WD="; WD
PRINT "X#(1)="; X#(1); "X#(2)="; X#(2); "X#(3)="; X#(3)
PRINT #1, "q="; q; "LL="; LL; "WD="; WD
PRINT #1, "X#(1)="; X#(1); "X#(2)="; X#(2); "X#(3)="; X#(3)
' "ДЛЯ НАЧАЛА ВЫЧИСЛЕНИЙ НАЖМИТЕ F5"
Mys0 = 0: Mys1 = 0: Mys2 = 0: Qys1 = 0: Qys2 = 0: Qys3 = 0
FOR I = 1 TO nk - 1
XI = (2 * I - 1) * tt / 2
IF I <= nnn THEN GOTO 15 ELSE 5
5 Qy = (-1) ^ (I + 1) * dP * tt * (R1 - X#(1)): Qysx1 = Qysx1 + Qy * XI
Qysx2 = Qysx2 + Qy * (XI) ^ 2: Qysx3 = Qysx3 + Qy * (XI) ^ 3