8990 REM *** HP TIME-SHARED BASIC PROGRAM LIBRARY ********************** 8991 REM 8992 REM INDMTR: ANALYSIS OF A BALANCED POLYPHASE INDUCTION 8993 REM MOTOR 8994 REM 36154 (A820) REV A -- 7/71 8995 REM 8996 REM *** CONTRIBUTED PROGRAM *************************************** 8997 REM 9000 PRINT "INDMTR: BALANCED POLYPHASE INDUCTION MOTOR ANALYSIS" 9001 PRINT " PREPARED BY C.M.SIEGEL, DEC. 31, 1969." 9002 PRINT "IF YOU WANT BRIEFING REMARKS, TYPE 1, OTHERWISE 0." 9003 INPUT X 9004 IF X#1 THEN 9029 9005 PRINT "GIVES CURRENT, EFFICIENCY, TORQUE, LOSSES, ELECTRICAL" 9006 PRINT "POWER IN, AND MECHANICAL POWER OUT AS FUNCTIONS OF SPEED" 9007 PRINT "WHEN MACHINE PARAMETERS OR NO-LOAD TEST DATA ARE GIVEN." 9008 PRINT 9009 PRINT "R1 AND R2 ARE STATOR AND ROTOR WINDING RESISTANCES." 9010 PRINT "X1 AND X2 ARE STATOR AND ROTOR LEAKAGE REACTANCES." 9011 PRINT "R3 AND X3 ARE MAGNETIZING LOSS RESISTANCE AND MUTUAL REACTANCE." 9012 PRINT 9013 PRINT "WHEN YOU SEE THE WORD, CONTINUE?, TYPE 1 FOR AFFIRMATIVE." 9014 PRINT "TYPING 0 WILL LET YOU RESUBMIT DATA." 9015 PRINT "BASE SPEED (RPM) MEANS YOUR CHOICE OF 100% VALUE OF ROTOR SPEED." 9016 PRINT "START(%) AND STOP(%) ARE THE BEGINNING AND END OF THE" 9017 PRINT " SPEED RANGE OF INTEREST (% OF BASE SPEED)" 9018 PRINT "INCREMENT(%) IS THE SIZE OF INCREMENTS BETWEEN START AND STOP" 9019 PRINT "FOR WHICH ANALYSIS IS WANTED (% OF BASE SPEED)" 9020 PRINT "WHEN YOU ARE ASKED FOR: " 9021 PRINT "BASE SPEED (RPM), START(%), STOP(%), INCREMENT(%)?--" 9022 PRINT "ILLUSTRATION: PUTTING IN" 9023 PRINT " 1800, 80, 100, 2 GIVES MOTOR PERFORMANCE" 9024 PRINT " FOR THE RANGE 1440 TO 1800 RPM AT 36 RPM INCREMENTS." 9025 PRINT 9026 PRINT "ALL ELECTRICAL WATTS, VOLTS, AMPS, ETC ARE PER PHASE." 9027 PRINT "ALL MECHANICAL HP AND TORQUE ARE TOTAL." 9028 PRINT 9029 PRINT "TYPE 1 IF YOU WANT TO PUT IN MACHINE PARAMETERS, " 9030 PRINT " 2 IF YOU WANT TO PUT IN TEST DATA, OR" 9031 PRINT " 99 IF YOU ARE FINISHED." 9032 INPUT X 9033 IF X=2 THEN 9083 9034 IF X=99 THEN 9156 9035 IF X#1 THEN 9029 9036 PRINT "R1, R2, R3, X1, X2, X3?" 9037 INPUT R1,R2,R3,X1,X2,X3 9038 PRINT 9039 PRINT "FOR F(CPS) CORRESPONDING TO X1,X2,X3, TYPE VALUES OF" 9040 PRINT "F(CPS), NO. OF PHASES, NO. OF POLES." 9041 INPUT F,N1,N2 9042 PRINT 9043 LET W=2*3.14159*F 9044 LET M=X3/W 9045 LET L1=M+X1/W 9046 LET L2=M+X2/W 9047 LET G3=1/R3 9048 GOSUB 9152 9049 PRINT "OPERATING CONDITIONS:" 9050 PRINT "VOLTS(LINE-TO-NEUTRAL; RMS), F (CPS)?" 9051 INPUT V,F 9052 PRINT 9053 PRINT "BASE SPEED (RPM), START (%), STOP (%), INCREMENT (%)?" 9054 INPUT S,S1,S2,S3 9055 PRINT 9056 GOSUB 9152 9057 LET J=(S2-S1)/S3+1 9058 PRINT "RPM","TORQUE","CURRENT","P.F.","VARS IN" 9059 PRINT " ","WATTS IN","H.P. OUT","WATTS LOSS","EFFICIENCY" 9060 FOR Z=1 TO J 9061 LET S5=S1+(Z-1)*S3 9062 LET W=2*3.14159*F 9063 W2=W-S5*S*N2*3.14159/6000 9064 GOSUB 9136 9065 LET I=V*SQR(Y1*Y1+Y2*Y2) 9066 LET P=V*V*Y1 9067 LET P5=Y1/SQR(Y1*Y1+Y2*Y2) 9068 LET T5=W*M*F2*I*I/J3 9069 LET T=M*M*W2*R2*N1*N2*550*I*I/((R2*R2+(W2*L2)^2)*2*746*J3) 9070 LET P8=M*F2*(W-W2)*I*I/J3 9071 LET P9=P8*N1/746 9072 LET P6=P-P8 9073 IF Y1 <= 0 THEN 9076 9074 LET E=P8/P 9075 GOTO 9077 9076 LET E=P/P8 9077 LET P7=(-1)*V*V*Y2 9078 LET S9=S5*S/100 9079 PRINT S9,T,I,P5,P7 9080 PRINT " ",P,P9,P6,E 9081 NEXT Z 9082 GOTO 9048 9083 PRINT "BLOCKED-ROTOR QUANTITIES:" 9084 PRINT "TORQUE(FT-LBS), WATTS IN, VOLTS IN, CURRENT IN?" 9085 INPUT T1,P1,V1,I1 9086 PRINT 9087 PRINT "FREQ(CPS), NO. OF PHASES, NO. OF POLES?" 9088 INPUT F,N1,N2 9089 PRINT 9090 PRINT "NO-LOAD QUANTITIES:" 9091 PRINT "SLIP(RPM), WATTS IN, VOLTS IN, CURRENT IN?" 9092 INPUT S8,P2,V2,I2 9093 GOSUB 9152 9094 LET T5=T1*4*F*3.14159*746/(N2*N1*550) 9095 LET Q1=SQR(V1*V1*I1*I1-P1*P1) 9096 LET Q2=SQR(V2*V2*I2*I2-P2*P2) 9097 LET W=2*3.14159*F 9098 LET R2=T5/I1^2 9099 LET L5=Q1/(2*W*I1^2) 9100 LET R1=(P1-T5)/I1^2 9101 LET G3=P2/V2^2 9102 LET L1=Q2/(W*I2*I2) 9103 LET L2=L1 9104 PRINT "RESULT OF SUCCESSIVE ITERATIONS FOLLOWS:" 9105 PRINT 9106 PRINT "R1","R2","R3","X1 AND X2","X3" 9107 PRINT "L1-M","M","L1" 9108 LET K=0 9109 LET M=L1-L5 9110 LET W2=W 9111 GOSUB 9136 9112 LET R2=T5*J3/(F3*I1*I1) 9113 LET L5=(Q1*J3/(W*I1*I1)-M*(F3-2*F1+1))/(F3+1) 9114 LET R1=(P1-T5)/I1^2-H3*G3/J3 9115 LET M=L1-L5 9116 LET W2=S8*N2*3.14159/60 9117 GOSUB 9136 9118 LET G3=P2*K3/(H3*V2*V2)-R1*J3/H3-H1/H3 9119 LET L1=J3*Q2/(W*I2*I2)+M*F1 9120 LET L2=L1 9121 LET R3=1/G3 9122 LET X3=W*(L1-L5) 9123 LET X1=W*L5 9124 LET M=L1-L5 9125 PRINT R1,R2,R3,X1,X3 9126 PRINT L5,M,L1 9127 LET K=K+1 9128 IF K<5 THEN 9109 9129 PRINT 9130 PRINT "IF YOU WANT TO CONTINUE THE ITERATION, TYPE 2" 9131 PRINT "IF YOU ARE READY FOR CALCULATION OF BEHAVIOR, TYPE 1" 9132 INPUT X 9133 IF X=2 THEN 9104 9134 IF X=1 THEN 9049 9135 GOTO 9029 9136 LET F4=R2*R2+(W2*L2)^2 9137 LET F1=W2*W2*M*L2/F4 9138 LET F2=W2*M*R2/F4 9139 LET F3=F1*F1+F2*F2 9140 LET H1=W*M*F2 9141 LET H2=W*(L1-M*F1) 9142 LET H3=H1*H1+H2*H2 9143 LET J1=1+G3*H1 9144 LET J2=G3*H2 9145 LET J3=J1*J1+J2*J2 9146 LET K1=H1+R1*J1 9147 LET K2=H2+R1*J2 9148 LET K3=K1*K1+K2*K2 9149 LET Y1=(J1*K1+J2*K2)/K3 9150 LET Y2=(J2*K1-J1*K2)/K3 9151 RETURN 9152 PRINT "CONTINUE?" 9153 INPUT X 9154 IF X#1 THEN 9029 9155 RETURN 9156 END