Проектирование и расчет валковых машин для полимерных материалов. Клинков А.С - 82 стр.

UptoLike

Программа 5
Расчет напряженно-деформационного состояния валка вальцев
99 CLS
100 WIDTH 80
101 DIM j(400), W(400), M(400), MK(400), J7(400), W7(400)
102 DIM SE(400), SI(400), SN(400), SS(400), SK(400)
103 DIM ST(400), ZX(400), X(400), BETA(400), M7(400)
104 DIM FP(400), XX(400), WW(400), JJ(400), MM(400)
105 DIM QS(400), MS(400), QQ(400), QM(400), MQ(400)
106 REM *Ввод исходных данных *
107 DATA 7,1,0.00,0.080,0.0001,0.04,0.00
108 DATA 1,0.00,0.095,0.04,0.11,.00
109 DATA 1,0.04,0.16,0.11,0.19,.04
110 DATA 4,0.04,0.16,0.19,0.35,.04
120 DATA 1,0.04,0.16,0.35,0.43,.04
130 DATA 1,0.04,0.95,0.43,0.394,.04
134 DATA 1,0.04,0.08,0.50,0.54,.04
140 READ N
150 FOR I = 1 TO N: FOR j = 1 TO 6
160 READ A(I, j)
165 NEXT j: NEXT I
170 DATA Сталь 40Л: READ MW$
180 DATA 0.54,0.12,0.42,36000,135e6 : READ L, L1, L2, Q, SD
190 DATA .024,.009,42,0.0525,0.016 : READ B1, T1, NP, ZP, DP
200 DATA 75,63,2.1E11,10.5E-6,.3,5.5,0.284
210 READ TV, TN, E, KT, KP, ND, V
220 REM * Расчет напряженного состояния валка *
230 RA = Q * (L2 - L1) / 2: PI = 3.14159: NV = V * 60 / PI / A(3, 3): KM = 9740 * ND / NV: Z = 0
235 hx1 = L / 100
240 FOR X = 0 TO L STEP hx1
250 Z = Z + 1: X(Z) = X * 1
260 I = 0
270 I = I + 1
280 IF A(I, 5) > X THEN 300
290 GOTO 270
300 RTF = TV - TN
310 P1 = E * KT * RTF / (2 * (1 - KP) * LOG(A(I, 3) / A(I, 5)))
320 P2 = LOG(A(I, 3) / A(I, 5))
330 ST(Z) = P1 * (1 - 2 * (A(I, 5) / 2) ^ 2 / ((A(I, 3) / 2) ^ 2 - (A(I, 5) / 2) ^ 2) * P2)
340 GOSUB 410: GOSUB 490: GOSUB 580: J7(Z) = j(Z): W7(Z) = W(Z): M7(Z) = M(Z)
350 ZX(M) = ZX(NP)
360 SI(Z) = M(Z) / W(Z): SS(Z) = 1.5 * SI(Z) + ST(Z): SK(Z) = MK(Z) / W(Z) / 2
370 SE(Z) = (1 - KP) / 2 * SS(Z) + (1 + KP) * SQR(SS(Z) ^ 2 + 4 * SK(Z) ^ 2) / 2
380 GOTO 870
390 REM * Подпрограмма *
400 REM * Определение крутящего момента в текущем сечении *
410 IF X < L1 THEN 440
420 IF X < L THEN 430
430 MK(Z) = KM: GOTO 460
440 MK(Z) = 0: GOTO 460
450 MK(Z) = KM * (X - L1) / L2