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

UptoLike

ВУЗ: 

ТГТУ | Тамбов

Составители: 

Рубрика: 

ПРИЛОЖЕНИЯ
ПРОГРАММА 1
' 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
'ограничений в виде равенств нет

Страницы