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

UptoLike

Рубрика: 

ПРОГРАММА 4
'*************************** ВалЛОПТ ************************
'* Минимизация целевой функции методом скользящего допуска *
'* Расчет *
'*************************************************************
DECLARE SUB FEAS ()
DECLARE SUB STAR ()
DECLARE SUB WRIT ()
DECLARE SUB SUMR ()
DECLARE SUB prob (INQ!)
DIM SHARED X#(100), 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 "ВалЛОПТ.res" FOR OUTPUT AS #1
PRINT #1, " Минимизация объема конструкции валка вальцев лабораторных "
PRINT " Минимизация объема конструкции валка вальцев лабораторных"
PRINT "-----------------------------------------------------------"
PRINT #1, " ------------------------------------------------------------"
FDIFER# = 1#
NX = 3 'общее число переменных
NC = 0 'общее число ограничений в виде равенств
NIC = 10 'общее число ограничений в виде неравенств
SIZE# = 3 'опред. размер деформ. многогранника
CONVER# = .001# 'используется для окончания поиска
ALFA# = 1#
BETA# = .5#
GAMA# = .95#
10 'PRINT "введите начальные значения варьируемых параметров, [мм]"
STEP0# = SIZE#
X#(1) = 40# 'координаты
X#(2) = 40# 'стартовой
X#(3) = 40# 'точки
'PRINT "ДЛЯ ВВОДА ИСХОДНЫХ ДАННЫХ НАЖМИТЕ F2 и выберите PROB"
SUB prob (INQ)
'"ввод распорного усилия [Н]"
PP = 80000
'"ввод геометрических параметров, [мм]"
D1 = 150: D2 = 500: D3 = 130: DD = 160: dd2 = 95: dd3 = 80
ll1 = 320: ll2 = 110: ll3 = 115
'"ввод минимального диаметра осевого отверстия, [мм]"
dd0 = 40
'"ввод технологической мощности, [Вт]"
NTex = 5000
'"ввод частоты вращения, [об/мин]"
nobor = 30
'"ввод допускаемого прогиба валка, [мм]"
fdop = .25
'"ввод допускаемого напряжения, [МПа]"
SIG = 150
'"ввод модуля упругости материала валка, [МПа]"
EE = 200000
'"ввод плотности материала валка, [Н/мм^3]"
RRO = 78500 / 1000000000
' "ДЛЯ НАЧАЛА ВЫЧИСЛЕНИЙ НАЖМИТЕ F5"
gg1 = 3.14 / 4 * (DD ^ 2 - X#(1) ^ 2) * ll1 * RRO
gg2 = 3.14 / 4 * (dd2 ^ 2 - X#(2) ^ 2) * ll2 * RRO
qg = gg1 / ll1 + 2 * gg2 / ll2
qp = PP / ll1
Fr = D1 / D3: Mkr = 9550 * NTex / nobor
Mkr2 = Mkr / (1 + Fr): Mkr1 = Mkr2 * Fr
'PRINT "qp="; qp; "qg="; qg; "Fr="; Fr; "Mkr="; Mkr; "Mkr1="; Mkr1; "Mkr2="; Mkr2
'PRINT #1, "qp="; qp; "Fr="; Fr; "Mkr="; Mkr; "Mkr1="; Mkr1; "Mkr2="; Mkr2
PP1 = 2 * Mkr1 / D1: PP2 = 2 * Mkr2 / D2