cgdz.bas

Gwbasic下的程序源代码

5 PRINT "KJ岔管I(2#-6-2)锥展开计算程序座表值  A;半锥角 R1,R2,R3,R31,R41上下切口正锥径"
8 PRINT "L;母线实长;A' AA':展开图母线与Y轴夹角  CN:C点对应园心角 AN:母线对应园心角"
10 A = 3.20521593#: CN = (180 - 91.60135316#): R1 = 1712.68: R2 = 1707.33#: RM = 2801.524
15 K = 3.14159265# / 180: P = 2: P1 = 1
20 R3 = 2699.53#: R31 = 2733.14: B2 = 15.09453: B21 = 10.06222#
30 B1 = ATN((R2 - R1) / TAN(A * K) / (R1 + R2))
35 R4 = R3 * (1 + TAN(B2 * K) * TAN(A * K)) / (1 - TAN(B2 * K) * TAN(A * K))
36 R41 = R31 * (1 + TAN(B21 * K) * TAN(A * K)) / (1 - TAN(B21 * K) * TAN(A * K))
40 H1 = R2 / TAN(A * K): H2 = R4 / TAN(A * K): H21 = R41 / TAN(A * K)
45 RR1 = R2 / SIN(A * K): RR2 = R4 / SIN(A * K): RR21 = R41 / SIN(A * K)
50 PRINT "A="; A; "R1="; R1; "R2="; R2; "R3="; R3; "R31="; R31; "R4="; R4
51 PRINT "R41="; R41; "B1="; B1 / K; "B2="; B2; "B21="; B21
55 C0 = INT(CN / 5 + .5) * 5
58 FOR C = 0 * K TO CN * K STEP 5 * K
60 GOSUB 300
65 NEXT C
75 R3 = R31: R4 = R41: B2 = B21: RR2 = RR21: P = P1: H2 = H21
80 R4 = R3 * (1 + TAN(B2 * K) * TAN(A * K)) / (1 - TAN(B2 * K) * TAN(A * K))
85 FOR C = C0 * K TO 180 * K STEP 5 * K
90 GOSUB 300
95 NEXT C
100 C = CN * K
110 GOSUB 300
120 END
300 IF C = 90 * K THEN 330
320 F = ATN(H1 / R2 / COS(C)): GOTO 350
330 F = 90 * K
350 F1 = SIN(F) * SIN(B1) * R2 * (1 - COS(C)) / (COS(A * K) * SIN(F - B1))
360 X1 = (RR1 - F1) * SIN(C * SIN(A * K))
370 Y1 = (RR1 - F1) * COS(C * SIN(A * K)) - RM / SIN(A * K)
375 IF P = 2 THEN 385
380 C2 = C: GOTO 386
385 C2 = 3.14159265# - C
386 IF C2 = 90 * K THEN 390
388 F0 = ATN(H2 / R4 / COS(C2)): GOTO 400
390 F0 = 90 * K
400 F2 = SIN(F0) * SIN(B2 * K) * R4 * (1 - COS(C2)) / (COS(A * K) * SIN(F0 - B2 * K))
410 X2 = (RR2 - F2) * SIN(C * SIN(A * K))
420 Y2 = (RR2 - F2) * COS(C * SIN(A * K)) - RM / SIN(A * K)
425 PRINT TAB(0); "X1="; X1; TAB(15); "Y1="; Y1; TAB(29); "X2="; X2; TAB(45); "Y2="; Y2; TAB(60); "L="; RR2 - RR1 - F2 + F1; TAB(45); "AN="; C / K; TAB(60); "A'="; C * SIN(A * K) / K
430 CC = 3.14159265# - C
440 XX1 = (RR1 - F1) * SIN(CC * SIN(A * K))
450 YY1 = (RR1 - F1) * COS(CC * SIN(A * K)) / SIN(A * K)
470 XX2 = (RR2 - F2) * SIN(CC * SIN(A * K))
480 YY2 = (RR2 - F2) * COS(CC * SIN(A * K)) - RM / SIN(A * K)
490 PRINT TAB(0); "XX1="; XX1; TAB(15); "YY1="; YY1; TAB(29); "XX2="; XX2; TAB(45); "YY2="; YY2; TAB(60); "AA'="; CC * SIN(A * K) / K
500 RETURN