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

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

' "расчет распределенной нагрузкии, сечения шнека, момента инерции и момента временного сопротив-
ления"
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
15 My = (-1) ^ (I + 1) * (dP * (TAN(BE) * (2 * R1 * X#(1) - X#(1) ^ 2) * 3.14 + (R1 - X#(1)) * tt) + 2 * dP * I
* FF * (2 * R1 * X#(1) - X#(1) ^ 2) * 3.14) * tt / 3.14
Mys0 = Mys0 +My: Mys1 = Mys1 + My * XI: Mys2=Mys2 + My * (XI) ^ 2
NEXT I
Ne = 3.14 ^ 2 * EE * J2 / ((2 * LL) ^ 2)
Np = PP * (3.14 * R1 ^ 2 + 2 * FF * (2 * R1 * X#(1) - X#(1) ^ 2) * 3.14)
IF (Ne) <= 0 THEN Ne = .0001: NpNe = (1 - Np / Ne)
IF (NpNe) <= 0 THEN NpNe = .5
ETA = 1 / (NpNe)
Wmax1 = q * LL ^ 4 * (3 / 4 + LL0 / LL - 1 / 4 * (LL0 / LL) ^ 4)
Wmax2a = (2 * LL0 * LL * Qysx1 + 3 * LL * Qysx2 - Qysx3) ^ 2
Wmax2b = (2 * LL0 * LL * Mys0 + 6 * LL * Mys1 - 3 * Mys2) ^ 2
Wmax2 = (Wmax2a + Wmaxb) ^ .5
Wmax = 1 / (6 * EE * J2) * (Wmax1 + Wmax2): Wpr = Wmax * ETA
Mmax = q * LL ^ 2 / 2 + ((Mys0) ^ 2 + Qysx1 ^ 2) ^ .5
Mkr = PP * 2 * 3.14 * (R1 ^ 3 - (R1 - X#(1)) ^ 3) / 3 * (TAN(BE) + FF * LL / tt)
SIG = ((Mmax / WW0 + Np / FFP + Np * Wpr / WW0) ^ 2 + (Mkr / WW0) ^ 2) ^ .5
PRINT "I="; I;"SIG="; SIG; "Wpr="; Wpr; "Wmax=";Wmax;"ETA="; ETA
PRINT "Mkr1="; Mkr1; "Mkr="; Mkr; "Mmax="; Mmax; "q="; q; "R#(9)="; R#(9)
PRINT #1, "I="; I; "SIG="; SIG; "Wpr="; Wpr; "Wmax="; Wmax; "ETA="; ETA
PRINT #1, "Mkr="; Mkr; "Mmax="; Mmax; "q="; q; "R#(9)="; R#(9)
IF INQ = 2 THEN GOTO 100
IF INQ = 3 THEN GOTO 200
'ограничений в виде равенств нет
GOTO 300
'ограничения в виде неравенств
100 R#(1) = 1 - SIG / SIGD
R#(2) = 1 - Wpr / WD
R#(3) = X#(1) / .001 - 1
R#(4) = 1 - X#(1) / .005
R#(5) = X#(2) / .001 - 1
R#(6) = 1 - X#(2) / .005
R#(7) = X#(3) / .001 - 1
R#(8) = 1 - X#(3) / .007
GOTO 300
'целевая функция

200 R#(9) = RO * gg * 3.14 * ((R1 - X#(1)) ^ 2 * LL0 + (R1 - X#(1)) ^ 2 * LL + 2 * (R1 - X#(1) / 2) /
COS(BE) * X#(1) * X#(2) * nk - X#(3) ^ 2 * (LL0 + LL))

300 END SUB

ПРОГРАММА 7

ПРОГРАММА РАСЧЕТА НА ПРОЧНОСТЬ ДВУХСЛОЙНЫХ ЦИЛИНДРОВ С ВИНТОВОЙ

КАНАВКОЙ В МЕСТЕ СОПРЯЖЕНИЯ. ЯЗЫК БЕЙСИК

‘ CalcT-LC.bas
' ПРОГРАММА РАСЧЕТА НА ПРОЧНОСТЬ ДВУХСЛОЙНЫХ ЦИЛИНДРОВ
' С ВИНТОВОЙ КАНАВКОЙ В МЕСТЕ СОПРЯЖЕНИЯ
' ввод исходных данных в системе измерения [СИ]
DATA .0225, .00005, .00045, .00055, .3, .0225, .0325, 566E6, 434E6, .6, 2E11, 150E6, .3
READ A, D, C, H0, HK, L, R0, RK, SIGMAT1, SIGMAT2, NT, E, PMAX, MU
' определение контактных давлений
Q1 = PMAX * 2 * A / D: Q2 = PMAX * 2 * A / C
' определение толстостенности наружного цилиндра
A1 = (SIGMAT2 * C / NT / A / PMAX) ^ 2 – .75
A2 = SQR(A1) + .1: A3 = SQR(A1) – 1.9: BETA = SQR(A2 / A3)
REM “определение наружного радиуса цилиндра”: B = C * BETA
' определение величины натяга
DELT = Q2 * C * RK * R0 / (R0 * HK + RK * 2 * H0) + Q2 * C * (1 – MU + BETA ^ 2 + BETA ^ 2 * MU) / 2 /
(BETA ^ 2 – 1) – PMAX * A * 2 * RK * R0 / (R0 * HK + RK * 2 * H0): DELT = DELT / E
' определение допускаемого контактного давления
QDOP = PMAX * (A / HO – (.5 – MU) * L ^ 2 / 2 / HO ^ 2)
A1 = (A / H0 – (.5 – MU) * L ^ 2 / 2 / H0 ^ 2 + .5) ^ 2 – (A / H0 + MU / 2 * L ^ 2 / H0 ^ 2) ^ 2 + (A / H0 + MU /
2 * L ^ 2 / H0 ^ 2) * (1 – L ^ 2 / 2 / HO ^ 2) + L ^ 2 / 2 / HO ^ 2 * (1 + L ^ 2 / 2 / HO ^ 2) + 1
A1 = A1 * PMAX ^ 2 + (SIGMAT1 / NT) ^ 2: A1 = SQR(A1)
QDOP = QDOP + A1: QDOP = QDOP * 2 * HO / D
' проверка выполнения условия по контактным давлениям
IF Q1 > QDOP THEN
PRINT "НЕ ВЫПОЛНЕНО УСЛОВИЕ ПО КОНТАКТНЫМ ДАВЛЕНИЯМ"
PRINT "НЕОБХОДИМО УВЕЛИЧИТЬ ПАРАМЕТРЫ:"
PRINT "d, c, H0, Hk, R0, Rk":

ELSE

PRINT "ИСХОДНЫЕ ДАННЫЕ: "
PRINT "РЕЗУЛЬТАТЫ РАСЧЕТА: "
'D,C,HO.HK,RO,RK')(X(5),A);

END IF: END

ПРОГРАММА 8

Программа минимизации массы конструкции двухслойного материального цилиндра с канавка-

ми в месте сопряжения. Язык Бейсик.

minMT-LC

'   Минимизация целевой функции методом скользящего допуска
'   Минимизация массы конструкции двухслойного материального цилиндра с канавками в местах со-
пряжения
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 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 "minVT-LC.res" FOR OUTPUT AS #1
PRINT #1, "   Расчет двухслойного цилиндра"
PRINT "-----------------------------------------------------------"
PRINT #1, " ------------------------------------------------------------"
FDIFER# = 1#
NX = 4 'общее число переменных
NC = 0 'общее число ограничений в виде равенств
NIC = 10 'общее число ограничений в виде неравенств
SIZE# = 5# 'опред. размер деформ. многогранника
CONVER# = .001# 'используется для окончания поиска
ALFA# = 1#
BETA# = .5#
GAMA# = 2#
10 STEP0# = SIZE#
X#(1) = .025# 'координаты
X#(2) = .033# 'стартовой
X#(3) = .042# 'точки
X#(4) = .008# '
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)
'Внутренний радиус внутреннего цилиндра, м
aa = .0225
'Плотность материала цилиндра, кг/м^3
RO = 7800
'Коэффициент Пуассона для стали
mu1 = .3
'Давление внутри цилиндра, Па
PP = 150000000
dd = X#(1): bb = X#(3): cc1 = X#(2): LL = X#(4)
PRINT cc1
bet = bb / cc1

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