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 
