WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


Pages:     | 1 |   ...   | 3 | 4 || 6 | 7 |

22. Шилд, Р. Методы оптимального проектирования конструкций / Р. Шилд. – М. : Механика, 1962. – 360 с.

23. Химмельблау, Д. Прикладное нелинейное программирование / Д. Химмельблау. – М. : Мир, 1975. – 480 с.

24. Разани, Р. Поведение равнонапряжённой конструкции и её отношение к конструкции минимального веса / Р. Разани // Ракетная техника и космонавтика. – 1965. – № 12.

25. Рывкин, А.А. Справочник по математике / А.А. Рывкин, А.З. Рывкин, Л.С. Хренов. – М. : Изд-во "Высшая школа", 1975. – 554 с.

26. Легостаев, В.Л. Методика раскрытия статической неопределённости в узлах сопряжения торообразных оболочек с кольцевыми днищами / Легостаев В.Л., Мордовин Е.Д. – Тамбов : Вестник ТГТУ. – 2006. – Т. 12 – № 2А.

27. Писаренко, Г.С. Справочник по сопротивлению материалов / Г.С. Писаренко, А.Б. Яковлев, А.А. Матвеев. – Киев : Наукова думка, 1975. – 704 с.

28. Расчёты машиностроительных конструкций методом конечных элементов : справочник / В.И. Мяченков, В.П. Мальцев, В.П. Майборода и др. – М. : Машиностроение, 1989. – 520 с.

29. Расчёт подпрессовочного устройства секторного форматоравулканизатора / В.И. Кочетов, А.С. Клинков, М.В. Соколов, Д.В. Туляков // Химическое и нефтегазовое машиностроение. – 2008. – № 10. – С. 18 – 21.

ПРИЛОЖЕНИЯ ПРОГРАММА 5 REM INPUT "РАСЧЁТ ЧЕРВЯКА НА ПРОЧНОСТЬ, ЖЁСТКОСТЬ 6 REM INPUT И УСТОЙЧИВОСТЬ. " 10 REM INPUT "Максимальное давление смеси Pmax=";PMAX 20 REM INPUT "Диаметp сеpдечника D2=";D25 REM INPUT "Наpужний диаметp D ="D 30 REM INPUT "Рабочая длина L=";L 40 REM INPUT "Модуль упpугости E=";E 50 REM INPUT "Шаг наpезки T=";T 60 REM INPUT "Коэффициент n=";ETA 70 REM INPUT "Допускаемое напряжение [G]=";G 80 REM INPUT "Допускаемый коэффициент запаса устойчивости для червяка" [Ny]=;NY 105 DATA 50,45,31,37,495,2E5,45,6.78,500,106 READ PMAX, D, D1, D2, L, E, T, ETA, G, X 100 REM NY = 120 P1 = 3.130 A = 2 * P1 / T: PRINT "A="; A 140 A1 = ((D2 – D1) / 4 / L) ^ 150 DSPX = (D + D1) / 2 + (D2 – D1) / 2 * X / L 160 DX = (D1 + (D2 – D1) / L * X) 170 A2 = (P1 * DSPX / T) ^ 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"; J220 J1 = P1 * D1 ^ 4 / 230 J2 = P1 * D2 ^ 4 / 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 * A1: PRINT "A1"; NX 270 NX = NX – PMAX * P1 * D ^ 2 / 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 ^ 300 A2 = A2 * X ^ 2 / L ^ 310 A2 = A2 – 1 / 320 A2 = D1 ^ 2 * A330 MX = DEL * PMAX * SIN(F1) * L / 8: PRINT MX 340 A3 = (1 / 3 – X / 4 / L) * X ^ 3 / L ^ 3 – 1 / 350 A3 = A3 * 2 * D1 * D360 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, A410 A1 = COS(A * X) / A ^ 2 + X * SIN(A * X) / A – 1 / A ^ 420 A1 = (D ^ 2 – D1 ^ 2) / L * A430 A2 = 2 * D1 * (D1 – D2) / L ^ 440 A3 = 2 * X * COS(A * X) / A ^ 2 + (X ^ 2 / A – 2 / A ^ 3) * SIN(A * X) – 2 * X / A ^ 450 A2 = A2 * A460 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 ^ 480 A3 = A3 + A490 A3 = (D2 – D1) ^ 2 / L ^ 3 * A500 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) 530 A2 = (D2 – D1) / L * (6 / A ^ 4 * COS(A * X) + 4 * X / A ^ 3 * SIN(A * X) – X ^ 2 / A ^ 2 * COS(A * X) – 6 / A ^ 4 + L * X / A ^ 2):

PRINT A540 A1 = A1 – A550 MZ = MZ * A560 A1 = DEL * PMAX * COS(F1) / 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 ^ 590 A4 = 2 * X * SIN(A * X) / A ^ 2 – (X ^ 2 / A – 2 / A ^ 3) * COS(A * X) + L * X / A – 2 / A ^ 600 A3 = A3 * A610 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 ^ 620 A4 = A4 * (D1 – D2) ^ 2 / L ^ 630 MZ = MZ – A1 * (A2 + A3 – A4) 635 PRINT "MZ"; MZ 636 REM ******************VY(X)************************************** 640 VY = 3 * DEL / 2 * PMAX * SIN(F1) * (D2 – D1) / L ^ 2 / A ^ 2 – 2 * D1 * (D2 – D1) * DEL * PMAX / 8 * COS(F1) / L ^ 2 / A 650 VY = VY / J0 / E 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 / J690 A1 = A1 * DEL / 8 * PMAX * COS(F1) / L ^ 3 / A * (D2 – D1) ^ 700 VY = VY – A703 PRINT "VY"; VY 704 REM ****************VZ(Z)*************************************** 710 VZ = (3 * DEL / 8 * PMAX * COS(F1) * (D2 – D1) ^ 2 / L ^ 3 / A ^ 2 – DEL / 2 * PMAX * SIN(F1) / L / A * (D – D1)) / E / J720 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)) 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 = A1 / E / J0 / 3 * 2 * DEL / 2 * PMAX * SIN(F1) / L ^ 2 / A * (D2 – D1) 750 VZ = VZ + A755 PRINT "VZ="; VZ 760 MOX = SQR(MY ^ 2 + MZ ^ 2) 770 VOX = SQR(VY ^ 2 + VZ ^ 2) 772 AWD = J1 / J773 PRINT "J1/J2="; AWD 780 SE = ETA / L ^ 2 * E * J790 VX = VOX / (1 – S / SE) 800 FX = P1 / 4 * DX ^ 810 wox = P1 * DX ^ 3 / 32: PRINT "wox="; wox 820 GMAX = 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 ^ 840 A1 = A1 * B * COS(F1) / L 850 A1 = A1 + P1 / 8 * L * D ^ 860 A2 = L / 2 + B * L ^ 2 + 6 * B ^ 2 * (L ^ 3 / 6 – L ^ 3 / 4 / P1 ^ 2) 861 A3 = 4 * B ^ 3 * (L ^ 4 / 8 – 3 / 8 * L ^ 4 / P1 ^ 2) 862 A2 = A2 + A880 A3 = B ^ 4 * (L ^ 5 / 10 – L ^ 5 / 2 / P1 ^ 2 + 3 / 4 / P1 ^ 4 * L ^ 5) 890 A2 = A2 + A900 A2 = A2 * E * J1 / L ^ 2 * P1 ^ 910 PRINT A2; A1: PMAXP = A2 / A920 TMAX = ABS(MX / 2 / wox) 930 GR = SQR(GMAX ^ 2 + 4 * TMAX ^ 2) 940 NU = PMAXP / PMAX 945 CLS PRINT "Максимальное давление смеси Pmax="; PMAX PRINT "Диаметp сеpдечника D2="; DPRINT "Наpужний диаметp D ="; D PRINT "Рабочая длина L="; L PRINT "Модуль упpугости E="; E PRINT "Шаг наpезки T="; T PRINT "Коэффициент n="; ETA PRINT "Допускаемое напpяжение [G]="; G REM PRINT "Допускаемый коэф-т запаса устойчивости для чеpвяка [Ny]="; NY PRINT "-------------------------------------------------" PRINT "Результаты: ( для X = "; X; ")" PRINT "-------------------------------------------------" 950 PRINT "Внутpеннее усилие Nx="; USING "+#.##^^^^"; NX: PRINT 955 PRINT "Mx="; USING "+#.##^^^^"; MX;

PRINT " "; : PRINT "My="; USING "+#.##^^^^"; MY;

PRINT " "; : PRINT "Mz="; USING "+#.##^^^^"; MZ: PRINT 960 PRINT "Пpогибы Vy="; USING "+#.#####"; VY;

PRINT " "; : PRINT "Vz="; USING "+#.#####"; VZ: PRINT 965 PRINT "Полный пpогиб Vx="; USING "+#.#####"; VX: PRINT REM PRINT "NY="; NY 970 PRINT "SIGMAmax="; USING "###.#"; GMAX;

PRINT " "; : PRINT "Pmaxкp="; USING "###.##"; PMAXP;

PRINT " "; : PRINT "TAUmax="; USING "###.##"; TMAX: PRINT 980 PRINT "NU="; USING "#.##"; NU; PRINT " "; : PRINT "SIGMAr=";

USING "###.##"; GR ПРОГРАММА Программа минимизации массы конструкции шнека. Язык BASIK ' 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 #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 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 =.' "Расчёт допускаемого прогиба [м]" WD =.01 * 2 * R' "ввод плотности [кг/м^3], модуля упругости [Па] и ускорения свобод.

паден. [м/с^2]" RO = 7850: EE = 2E+11: gg = 9.' "расчёт распределённой нагрузкии, сечения шнека, момента инерции и момента временного сопротивления" 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) ^.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 = FOR I = 1 TO nk – XI = (2 * I – 1) * tt / IF I <= nnn THEN GOTO 15 ELSE 5 Qy = (-1) ^ (I + 1) * dP * tt * (R1 – X#(1)): Qysx1 = Qysx1 + Qy * XI Qysx2 = Qysx2 + Qy * (XI) ^ 2: Qysx3 = Qysx3 + Qy * (XI) ^ 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.Mys0 = Mys0 +My: Mys1 = Mys1 + My * XI: Mys2=Mys2 + My * (XI) ^ 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 =.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) ^ Wmax2b = (2 * LL0 * LL * Mys0 + 6 * LL * Mys1 – 3 * Mys2) ^ Wmax2 = (Wmax2a + Wmaxb) ^.Wmax = 1 / (6 * EE * J2) * (Wmax1 + Wmax2): Wpr = Wmax * ETA Mmax = q * LL ^ 2 / 2 + ((Mys0) ^ 2 + Qysx1 ^ 2) ^.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) ^.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 IF INQ = 3 THEN GOTO 'ограничений в виде равенств нет GOTO 'ограничения в виде неравенств 100 R#(1) = 1 – SIG / SIGD R#(2) = 1 – Wpr / WD R#(3) = X#(1) /.001 – R#(4) = 1 – X#(1) /.R#(5) = X#(2) /.001 – R#(6) = 1 – X#(2) /.R#(7) = X#(3) /.001 – R#(8) = 1 – X#(3) /.GOTO 'целевая функция 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 ПРОГРАММА **************************** ITERA ************************* Решение системы уравнений итерационным методом Расчёт составного цилиндра с водяной рубашкой (вариант 1) ************************************************************* DeltaX# =.05# R# = 2# Eps# =.001# m% = Pi# = 3.1415926536# OPEN "iterase.res" FOR OUTPUT AS #PRINT #1, " * * * P R O G R A M I T E R A * * * ";

PRINT #1, MID$(DATE$, 4, 2); : PRINT #1, ".";

PRINT #1, MID$(DATE$, 1, 2); : PRINT #1, ".";

PRINT #1, MID$(DATE$, 9, 2) PRINT #1, " Расчёт составного цилиндра с водяной рубашкой (вариант 1)" PRINT #1, "---------------------------------------------------------------------------" PRINT #1, " X1 X2(K1) X2(K2) Число итераций" PRINT #1, "---------------------------------------------------------------------------" DO WHILE m% <= X1# = 25# + (m% – 1) * DeltaX# X2K1# = 32#: X2K2# = 32#: Ip% = DO Q1# = 150# * 2# * 22.5# / X1# Q2# = 150# * 2# * 22.5# / X2K1# Bt# = 42# / X2K1# Dm1# = 150# * (22.5# / (X1# – 22.5#) +.3# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ 2) – Q1# * X1# / 2# / (X1# – 22.5#) – 150# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ Dm2# = 150# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ 2 + 150# Dm3# = 150# * (22.5# / (X1# – 22.5#) +.3# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ 2) – Q1# * X1# / 2# / (X1# – 22.5#) + 150# Dm4# = (Bt# ^ 2 + 1# / (Bt# ^ 2 – 1#)) ^ 2 + 1.8# * (Bt# ^ 2 + 1# / (Bt# ^ – 1#)) + 1.56# SigmaE1# = SQR(.5# * (Dm1# ^ 2 + Dm2# ^ 2 + Dm3# ^ 2)) SigmaE2# = SQR(Q2# ^ 2 * Dm4# / 4#) K1# = SigmaE1# / 556# Q1# = 150# * 2# * 22.5# / X1# Q2# = 150# * 2# * 22.5# / X2K2# Bt# = 42# / X2K2# Dm1# = 150# * (22.5# / (X1# – 22.5#) +.3# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ 2) – Q1# * X1# / 2# / (X1# – 22.5#) – 150# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ Dm2# = 150# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ 2 + 150# Dm3# = 150# * (22.5# / (X1# – 22.5#) +.3# * 8# ^ 2 / 2# / (X1# – 22.5#) ^ 2) – Q1# * X1# / 2# / (X1# – 22.5#) + 150# Dm4# = (Bt# ^ 2 + 1# / (Bt# ^ 2 – 1#)) ^ 2 + 1.8# * (Bt# ^ 2 + 1# / (Bt# ^ – 1#)) + 1.56# SigmaE1# = SQR(.5# * (Dm1# ^ 2 + Dm2# ^ 2 + Dm3# ^ 2)) SigmaE2# = SQR(Q2# ^ 2 * Dm4# / 4#) K2# = SigmaE2# / 434# IF ABS(K1# – 1#) <= Eps# AND ABS(K2# – 1#) <= Eps# THEN EXIT DO X2K1# = X2K1# + X2K1# * (K1# – 1#) / R# X2K2# = X2K2# + X2K2# * (K2# – 1#) / R# Ip% = Ip% + IF Ip% > 49 THEN EXIT DO LOOP PRINT #1, USING "#########.## "; X1#; X2K1#; X2K2#; Ip% PRINT USING "#########.## "; X1#; X2K1#; X2K2#; Ip% m% = m% + LOOP PRINT #1, "---------------------------------------------------------------------------" PRINT "----------------------------------------------------------------------------" V# = 2# * Pi# * (25.63# ^ 2 – 22.5# ^ 2) + Pi# * (31.62# ^ 2 – 25.63# ^ 2) + 2# * Pi# * (42# ^ 2 – 31.62# ^ 2) PRINT #1, " V = "; V# PRINT " V = "; V# CLOSE ПРОГРАММА ' 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 #PRINT #1, " Расчёт двухслойного цилиндра" PRINT "-----------------------------------------------------------" PRINT #1, " ------------------------------------------------------------" FDIFER# = 1# NX = 4 'общее число переменных NC = 0 'общее число ограничений в виде равенств NIC = 10 'общее число ограничений в виде неравенств SIZE# = 5# 'опред. размер деформ.

Pages:     | 1 |   ...   | 3 | 4 || 6 | 7 |






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.