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

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 

      PP1 = ((aa / (dd - aa) + mu1 * LL ^ 2 / 2 / (dd - aa) ^ 2) - 2 * aa / dd * dd / 2 / (dd - aa) - LL ^ 2 / 2 / (dd - 
aa) ^ 2) ^ 2 
      PP2 = ((LL ^ 2 / 2 / (dd - aa) ^ 2 + 1)) ^ 2 
      PP3 = ((aa / (dd - aa) + mu1 * LL ^ 2 / 2 / (dd - aa) ^ 2) - 2 * aa / dd * dd / 2 / (dd - aa) + 1) ^ 2 
      Sig1 = PP * SQR(.5 * (PP1 + PP2 + PP3)) 
      PRINT Sig1 
      Sig2 = PP * SQR(aa ^ 2 / cc1 ^ 2 * (((bet ^ 2 + 1) / (bet ^ 2 - 1)) ^ 2 + 1.8 * ((bet ^ 2 + 1) / (bet ^ 2 - 1)) + 
1.56)) 
 'Допускаемые напряжения для материала двухслойного цилиндра, Па 
      Sigd1 = 566000000: Sigd2 = 434000000 
      IF INQ = 2 THEN GOTO 100 
      IF INQ = 3 THEN GOTO 200 
      'ограничений в виде равенств нет 
      GOTO 300 
      'ограничения в виде неравенств 
100   
      R#(1) = 1# - Sig1 / Sigd1 
      R#(2) = 1# - Sig2 / Sigd2 
      R#(3) = X#(1) / .025 - 1# 
      R#(4) = 1# - X#(1) / .032 
      R#(5) = X#(2) / .03 - 1# 
      R#(6) = 1# - X#(2) / .037 
      R#(7) = X#(3) / .042 - 1# 
      R#(8) = 1# - X#(3) / .052 
      R#(9) = X#(4) / .008 - 1# 
      R#(10) = 1# - X#(4) / .015 
      GOTO 300 
      'целевая функция 
200   R#(11) = 3.14 * RO * (2 * (dd ^ 2 - aa ^ 2) + (cc1 ^ 2 - dd ^ 2) + 2 * (bb ^ 2 - cc1 ^ 2)) / 1000 
300   X#(1) = dd: X#(3) = bb: X#(2) = cc1: X#(4) = LL 
      END SUB