содержание .. 38 39 40 41 ..

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

170 A2 = (P1 * (DSPX / T)) ^ 2
180 DEL = SQR(A1 + A2 + 1)
190 F1 = T / P1 / D
200 B = (D2 – D1) / D1 / L: PRINT "B"; B
210 J0 = P1 * (D1 ^ 4) / 64: PRINT "J0"; J0
220 J1 = P1 * (D1 ^ 4) / 64
230 J2 = P1 * (D2 ^ 4) / 64
240 S = (PMAX * P1 * (D ^ 2)) / 4: PRINT "S1"; S
250 A1 = (((D – D1) / 4) * ((X ^ 2) / (L ^ 2) – 1)) – (((D2 – D1) / 6) * ((X ^ 3) / (L ^ 3) – 1))
260 NX = DEL * PMAX * COS(F1) * L: PRINT "A1"; NX
265 NX = NX * A1
270 NX = NX – ((PMAX * P1 * D * D) / 4)
275 REM **************Mx(X)************************************************
280 A1 = ((D ^ 2) / 2) * ((X ^ 2) / (L ^ 2) – 1)
290 A2 = .5 – (((2 / 3) * (X / L)) + ((X ^ 2) / (4 * (L ^ 2))))
300 A2 = A2 * ((X ^ 2) / (L ^ 2))
310 A2 = A2 – (1 / 12)
320 A2 = (D1 ^ 2) * A2
330 Mx = (DEL * PMAX * SIN(F1) * L) / 8: PRINT Mx
340 A3 = ((1 / 3) – (X / 4 / L)) * ((X ^ 3) / (L ^ 3)) – (1 / 12)
350 A3 = A3 * 2 * D1 * D2
360 A4 = D2 ^ 2 / 4 * (X ^ 4 / L ^ 4 – 1)
370 Mx = (-Mx) * (A1 – A2 – A3 – A4)
373 PRINT "NX"; NX, "Mx"; Mx
376 REM ***************MY(X)********************************************
380 A1 = (D – D1) * ((((-2) * COS(A * X)) / (A ^ 3)) – ((X * SIN(A * X)) / (A ^ 2)) + (2 / (A ^ 3)))
390 A2 = ((D2 – D1) / L) * (((6 / (A ^ 4)) * SIN(A * X)) – (((4 * X) / (A ^ 3)) * COS(A * X)) – (((X ^

2) / (A ^ 2)) * SIN(A * X)) + ((4 * X) / (A ^ 3)))

400 MY = ((DEL * PMAX * SIN(F1)) / (2 * L)) * (A1 – A2): PRINT "A1"; A1, A2
410 A1 = (COS(A * X) / (A ^ 2)) + ((X * SIN(A * X)) / A) – (1 / (A ^ 2))
420 A1 = ((D ^ 2) – (D1 ^ 2)) / L * A1
430 A2 = (2 * D1 * (D1 – D2)) / (L ^ 2)
440 A3 = (((2 * X * COS(A * X)) / (A ^ 2)) + (((X ^ 2) / A) – (2 / (A ^ 3))) * SIN(A * X) – ((2 * X) /

(A ^ 2)))

450 A2 = A2 * A3
460 A3 = ((3 * ((X ^ 2) / (A ^ 2))) – (6 / (A ^ 4))) * COS(A * X)
470 A4 = ((((X ^ 3) / A) – ((6 * X) / (A ^ 3))) * SIN(A * X)) + ((6 / (A ^ 4)) – ((3 * L * X) / (A ^ 2)))
480 A3 = A3 + A4
490 A3 = (((D2 – D1) ^ 2) / (L ^ 3)) * A3
500 MY = MY – ((DEL * PMAX * COS(F1)) / 8) * (A1 + A2 – A3)
502 PRINT "MY"; MY
505 REM ******************MZ(Z)*********************************************
510 MZ = -DEL * PMAX * SIN(F1) / 2 / L
520 A1 = (D – D1) * (((2 * SIN(A * X)) / (A ^ 3)) – ((X * COS(A * X)) / (A ^ 2)) + (X / (A ^ 2)))
525 C = -(((X ^ 2) / (A ^ 2)) * COS(A * X)) – (6 / (A ^ 4)) + ((L * X) / (A ^ 2))
530 A2 = ((D2 – D1) / L) * ((6 / (A ^ 4)) * COS(A * X) + (((4 * X) / (A ^ 3)) * SIN(A * X)) + C):

PRINT A1

540 A1 = A1 – A2
550 MZ = MZ * A1
560 A1 = (DEL * PMAX * COS(F1)) / 8
570 A2 = (((D ^ 2) – (D1 ^ 2)) / L) * ((SIN(A * X) / (A ^ 2)) – ((X * COS(A * X)) / A) + X / A)

580 A3 = 2 * D1 * (D1 – D2) / (L ^ 2)
590 A4 = ((2 * X * SIN(A * X)) / (A ^ 2)) – (((X ^ 2) / A) – (2 / (A ^ 3))) * COS(A * X) + ((L * X) /

A) – (2 / (A ^ 3))

600 A3 = A3 * A4
610 A4 = ((3 * ((X ^ 2) / (A ^ 2)) – (6 / (A ^ 4))) * SIN(A * X)) – (((X ^ 3) / A) – (6 * X / (A ^ 3))) *

COS(A * X) + (((L ^ 2) * X) / A) – ((6 * X) / (A ^ 3))

620 A4 = A4 * (((D1 – D2) ^ 2) / (L ^ 3))
630 MZ = MZ – A1 * (A2 + A3 – A4)
635 PRINT "MZ"; MZ
636 REM ******************VY(X)*****************************************
640 VY = (3 * DEL * (PMAX / 2) * SIN(F1) * ((D2 – D1) / (L ^ 2) / (A ^ 2))) – (2 * D1 * (D2 – D1) *

(DEL * (PMAX / 8) * COS(F1) / (L ^ 2) / A))

660 VY = VY * (((X ^ 4) / 12) – ((L / 6) * (X ^ 3)) + ((L ^ 3) * (X / 12))) – (4 * B) * (((X ^ 5) / 20) –

((L * (X ^ 4)) / 12) + ((L ^ 4) * (X / 30))) + (10 * (B ^ 2)) * (((X ^ 6) / 30) – ((L * (X ^ 5)) / 20) + (((L ^ 5) *
X) / 60)): PRINT "VY"; VY

670 A1 = ((X ^ 5) / 20 – (((L ^ 2) * (X ^ 3)) / 6) + ((14 * (L ^ 4) * X) / 120)) – (4 * B) * (((X ^ 6) / 30)

– (((L ^ 2) * (X ^ 4)) / 12) + ((18 * (L ^ 5) * X) / 360)) + (10 * (B ^ 2)) * (((X ^ 7) / 42) – (((L ^ 2) * (X ^
5)) / 20) + ((22 * X * (L ^ 6)) / 840))

680 A1 = A1 / E / J0: VY = VY / (J0 * E)
690 A1 = A1 * ((((DEL / 8) * PMAX * COS(F1)) / ((L ^ 3) * A)) * ((D2 – D1) ^ 2))
700 VY = VY – A1
703 PRINT "VY"; VY
704 REM ****************VZ(Z)***************************************
710 VZ = ((3 * (DEL / 8) * PMAX * COS(F1)) / ((L ^ 3) * (A ^ 2)) * ((D2 – D1) ^ 2)) – ((2 / 3) *

PMAX * SIN(F1) * (D – D1) * (1 / (2 * L * A)))

720 VZ = VZ * (((X ^ 4) / 12 – (L / 6) * (X ^ 3) + ((L ^ 3) / 12) * X) – (4 * B) * ((X ^ 5) / 20 – (L * (X

^ 4)) / 12 + ((L ^ 4) * X) / 30) + (10 * (B ^ 2)) * ((X ^ 6) / 30 – (L * (X ^ 5)) / 20 + ((L ^ 5) / 60) * X))

725 VZ = VZ / E / J0
730 A1 = (X ^ 5) / 20 – ((L ^ 2) * (X ^ 3)) / 6 + (14 * (L ^ 4) * X) / 120 – (4 * B) * ((X ^ 6) / 30 – ((L ^

2) * (X ^ 4)) / 12 + (18 * (L ^ 5) * X) / 360) + 10 * (B ^ 2) * ((X ^ 7) / 42 – ((L ^ 2) * (X ^ 5)) / 20 + (22 * X
* (L ^ 6)) / 840)

740 A1 = (1 / E / J0 / 3) * A1 * DEL * (PMAX * SIN(F1) / ((L ^ 2) * A)) * (D2 – D1)
750 VZ = VZ + A1
755 PRINT "VZ = "; VZ
760 MOX = SQR(MY ^ 2 + MZ ^ 2)
770 VOX = SQR(VY ^ 2 + VZ ^ 2)
772 AWD = J1 / J2
773 PRINT "J1/J2"; AWD
780 SE = (E / (L ^ 2)) * ETA * J2
790 VX = VOX / (1 – NX / SE)
800 FX = (P1 / 4) * DX
810 WOX = (P1 / 32) * (DX ^ 3)
820 GMAX = ABS(ABS(NX) / FX + MOX / WOX + (ABS(NX) / WOX) * VX)
830 A1 = (((L ^ 3) / 12) * (D – D1)) – ((L ^ 3) / 32) * (D2 – D1) – ((D – D1) * (L ^ 3)) / (16 * (P1 ^ 2))

+ ((D2 – D1) / (16 * (P1 ^ 2))) * L ^ 3

840 A1 = (A1 / L) * B * COS(F1)
850 A1 = A1 + ((P1 / 8) * L * (D ^ 2))
860 A2 = (L / 2) + B * (L ^ 2) + 6 * (B ^ 2) * (((L ^ 3) / 6) – ((L ^ 3) / (4 * (P1 ^ 2)))) + (4 * (B ^ 3)) *

(((L ^ 4) / 8) – ((3 / 8) * ((L ^ 4) / (P1 ^ 2))))

880 A3 = (B ^ 4) * ((L ^ 5) / 10 – (L ^ 5) / (2 * (P1 ^ 2)) + (3 / 4 / P1 ^ 4) * L ^ 5)
890 A2 = A2 + A3
900 A2 = ((A2 * E * J1) / (L ^ 2)) * (P1 ^ 2): PMAXP = A2 / A1

920 TMAX = ABS(Mx / 2 / WOX): GR = SQR((GMAX ^ 2) + 4 * (TMAX ^ 2))
940 NU = PMAXP / PMAX
950 PRINT "Внутpеннее усилие NX = "; NX, "MX = "; Mx, "MY = "; MY, "MZ = "; MZ
960 PRINT "Пpогибы VY = "; VY, "VZ = "; VZ, "Полный пpогиб VX = "; VX
970 PRINT "SIGMAMAX = "; GMAX, "Pmaxкp = "; PMAXP, "TAUmax = "; TMAX; "NU = "; NU;

"SIGMAR = "; GR

ПРОГРАММА 6

Программа минимизации массы конструкции шнека с разрывными витками. Язык Бейсик.

' minMSCRE
'Минимизация целевой функции методом скользящего допуска
' Расчет минимальной массы шнека с разрывными витками
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 "minMSCRE.res" FOR OUTPUT AS #1
PRINT #1, "Расчет минимальной массы шнека с разрывными витками "
PRINT " Расчет минимальной массы шнека с разрывными витками "
FDIFER# = 1#
NX = 3 'общее число переменных
NC = 0 'общее число ограничений в виде равенств
NIC = 8 'общее число ограничений в виде неравенств
SIZE# = 5# 'опред. размер деформ. многогранника
CONVER# = .001# 'используется для окончания поиска
ALFA# = 1#
BETA# = .5#
GAMA# = 1#
10 STEP0# = SIZE#
X#(1) = .003# 'координаты
X#(2) = .003# 'стартовой
X#(3) = .005# 'точки
' CALL prob(3)
'PRINT "ДЛЯ  ВВОДА  ИСХОДНЫХ  ДАННЫХ  НАЖМИТЕ F2 И  ВЫБЕРИТЕ  ПОДПРОГРАММУ
PROB"
IF FDIFER# < CONVER# THEN GOTO 9999
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

содержание .. 38 39 40 41 ..