求助福建单相电机程序操作说明

xuefengexpt · 发表于 2012-9-2 11:13

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？立即注册

×

本人看了福建单相电机程序，不清楚操作，求助福建单相电机程序操作说明。

xuefengexpt · 发表于 2012-9-2 11:20

[福建单相电机程序实例，

xuefengexpt · 发表于 2012-9-2 11:22

刚才没有附上实例：福建单相电机程序实例，

xuefengexpt · 发表于 2012-9-2 11:28

福建单相电机程序:实例。

wyx · 发表于 2012-9-11 07:34

应该是广科所的单相程序，本人一直在用

wenzee · 发表于 2012-9-13 20:28

这个程序。本人可以说很精通。而且还改过多处。
我在福安呆了6年。

lzemg007 · 发表于 2012-9-15 13:56

在福安呆了6年？谁啊？

wenzee · 发表于 2012-9-16 07:17

lzemg007 发表于 2012-9-15 13:56

登录/注册后可看大图

在福安呆了6年？谁啊？

杨心文。
设计单相是一直用这程序。

10 DIM A(100),B(100),C(100),D(100),E(102),B1(44),B2(44),B3(44)
20 DIM X(200),Y(200),H(100)
30 DIM A$(60),B$(96),C$(93),D$(93),E$(102)
'定义数组
40 RESTORE 40000 '读取第40000行后的58+95+91+93+100个输出结果字符串变量
42 FOR J=1 TO 58:READ A$(J):NEXT J
44 FOR J=1 TO 95:READ B$(J):NEXT J
46 FOR J=1 TO 91:READ C$(J):NEXT J
48 FOR J=1 TO 93:READ D$(J):NEXT J
50 FOR J=1 TO 100:READ E$(J):NEXT J
'定义char类型变量,输出变量的变量名
60 RESTORE 18000 '读取第18000行的输入数据：额定数据，冲片数据
65 READ NA$ 'NA$ 型号
70 FOR J=1 TO 31:READ A(J):NEXT J
'读取输入的额定数据及冲片数据
80 READ A(51),A(52),A(44),A(37)
'读取端环数据A1, B1, B2及A(37)-hs1
90 FOR J=1 TO 22:READ B(J):NEXT J
'读取主绕组数据
100 FOR J=1 TO 22:READ D(J):NEXT J
'读取副绕组数据
110 READ B(69),C(42) 'B(69)绝缘层厚度，C(42)叠压系数
115 PRINT "k2=? K3=? K7=? K10=?":INPUT K2,K3,K7,K10
'K2定子槽型， K3转子槽型， K7计算方法：1手动输入方式，2按18010、18020的数据自动计算，
'K10电机类型:1-YU, 2-YC, 3-YY
120 PI=3.14159 '定义π
125 IF K10<>3 THEN 190 '非电容运转跳至190行
130 PRINT "C=":INPUT E(17) '输入E(17)-Cr
140 E(18)=5*10^5/3.14159/A(3)/E(17) 'E(18)-Xc = 1 / (2*π*f*Cr), A(3)-frequcney, E(17)-Cr
'150 RESTORE 15000 '读取15000行的17个电容值
'160 U=E(17) 'U=E(17)-Cr
'170 GOSUB 7000
180 E(19)=V 'Rc =SUB7000 (Cr)插值求得,应该是 Rc=Ec*0.0022
190 A(29)=A(27)-2*A(28) 'D2=Di1-2*g A(29)-D2, A(27)-Di1, A(28) -g
200 A(32)=PI*A(27)/A(4) 'τp=π*Di1/P A(32)-τp 极距
210 A(33)=PI*A(27)/A(10) 'T1=π*Di1/Q1 A(33)-T1 定子齿距
220 A(34)=PI*A(29)/A(11) 'T2=π*D2/Q2 A(34)-T2转子齿距
230 A(35)=A(10)/A(4) 'Qp1=Q1/P A(35)-Qp1 毎极定子槽数
240 A(36)=A(11)/A(4) 'Qp2=Q2/P A(36)-Qp2毎极转子槽数
245 IF A(17)=0 THEN 260 'If (α=0) 定子槽斜肩角度=0 then 260
250 A(37)=(A(14)-A(13))*TAN(A(17)*PI/180)/2 'A(37)-hs1 = (bs1-bs0)*tg(α*π/180)/2
260 A(38)=PI*(A(27)+2*(A(16)+A(37)))/A(10)-A(14)
' bt11=π*(Di1+2*(hs0+hs1)/Q1 A(38)-bt11
270 A(39)=PI*(A(27)+2*(A(16)+A(37)+A(18)))/A(10)-A(15)
' bt12 =π*(Di1+2*(hs0+hs1+hs2)/Q1-bs2 A(39)-bt12
280 IF A(38)<A(39) THEN 288 'if(bt11 < bt12) then 288
285 A(40)=A(39):GOTO 290 'bt1 = bt12 then 290 A(40)-bt1
288 A(40)=A(38) 'bt1 = bt11
290 A(40)=A(40)+ABS(A(38)-A(39))/3 'bt1 = bt1+ABS(bt11-bt12)/3
300 A0=A(16)+A(37)+A(18) 'hs = hs0+hs1+hs2 A0-hs定子槽高
310 A(41)=(A(26)-A(27))/2-A0 'hc1 = (D1-Di1)/2-hs A(41)-hc1定子轭高, 平底槽
315 IF K2=2 THEN 320 'if(定子槽型为平底槽) then 320
317 A(41)=A(41)-2*A(19)/3 'hc1=hc1-2*R/3 圆底槽
320 A(49)=2*A(21) 'br1 =2*r1 A(49)-br1
330 A(42)=PI*(A(29)-2*(A(23)+A(24)))/A(11)-A(49)
' bt21=π*(D2-2*(hr0+hr1))/Q2-br1 A(42)-bt21
340 A(50)=2*A(22) 'br2=2*r2 A(50)-br2
350 A(43)=PI*(A(29)-2*(A(23)+A(24)+A(25)))/A(11)-A(50)
' bt22=π*(D2-2*(hr0+hr1+hr2))/Q2-br2 A(43)-bt22
360 IF A(42)<A(43) THEN 368 'if (bt21 < bt22) then 368
365 A(45)=A(43):GOTO 370 'bt2 = bt22 A(45)-bt2
368 A(45)=A(42) 'bt2 = bt21
370 A(45)=A(45)+ABS(A(42)-A(43))/3 ' bt2 = bt2 + ABS(bt21-bt22)/3
380 A(46)=A(23)+A(24)+A(25):IF K3=2 THEN 390
' hr=+hr0+hr1+hr2 A(46)-hr, if(转子槽型为平底槽) then 390
385 A(46)=A(46)+A(50)/3 'hr = hr+br2/3
390 A(47)=(A(29)-A(30))/2-A(46)+A(30)/4 'hc2= (D2-Di2)/2-hr+Di2/4 A(47)-hc2
395 A(48)=(A(20)+A(49))*A(24)/2+(A(49)+A(50))*A(25)/2:IF K3=2 THEN 410
' S2= (bo2+br1)*hr1/2+(br1+br2)*hr2/2 平底槽 if(转子为平底槽) then 410 A48)-s2
400 A(48)=PI*(A(49)^2+A(50)^2)/8+(A(49)+A(50))*A(25)/2:IF K3=1 THEN 410
' S2= π*(br1*br1+br2*br2)/8+(br1+br2)*hr2/2 圆底槽 if(转子为圆底槽) then 410
405 A(48)=(A(20)+A(49))*A(24)/2+(A(49)+A(50))*A(25)/2+PI*A(50)^2/8
' S2 = (bo2+br1)*hr1/2+(br1+br2)*hr2/2+π*br2*br2/8 平底槽
410 A(53)=A(51)*(A(52)+A(44))/2 'Sr = A1*(B1+B2)/2 A(53)-Sr
420 A(54)=A(29)-2*.05 'Dr0=D2-2*0.05 A(54)-Dr0
430 A(55)=A(54)-2*A(52) 'Dri=Dr0-2*B1 A(55)-Dri
440 A(56)=A(54)-A(52) 'Dr=Dr0-B1 A(56)-Dr
450 A(57)=A(55)/A(54) 'Karing = Dri/Dr0 A(57)-Karing
460 A(58)=A(4)*(1-A(57))/2*(1+A(57)^A(4))/(1-A(57)^A(4))
' Kring = (P*(1-karing)/2)*(Pow((1+karing), P)/(Pow((1-karing), P)
560 N4=0:N3=0 '
570 B(23)=0 'Ymav=0 B(23)-Ymav
580 FOR I=1 TO B(22)
590 B(23)=B(I)*B(10+I)+B(23) 'Ymav=Ymav+Ym(i)*Km(i)%
600 NEXT I
610 PRINT "Bg'= Ks'= Ke'= Jm'=":INPUT B(24),B(25),B(29),B(62)
'input B(24)-Bg', B(25)-ks', B(29)-ke', B(62)-Jm'
615 PRINT "a=? Ka=? s=?":INPUT D(23),D(53),E(32)
'Input D(23)-a, D(53)-Ka, E(32)-S'
617 PRINT "a1=? a2=?":INPUT A1,A2 'Input主副绕组并联支路数 a1, a2
620 U=B(25) 'B(25)-ks'饱和系数
630 RESTORE 10000
640 GOSUB 7000
650 B(26)=V 'αi 由Sub7000(ks')插值求得,数据为10000，100010 B(26)-αi极弧系数
660 GOSUB 7000
670 B(27)=V 'fb由Sub7000(ks')插值求得,数据为10020，100030 B(27)-fb波形系数
675 IF K7=1 THEN 735 'if (手动) then735
680 B(28)=B(26)*B(24)*A(32)*A(31) 'PHI'=αi*Bg'*Tp*L1 B(28)-PHI'毎极磁通
690 B(30)=A1*B(29)*A(2)*50*10^6/B(27)/B(28)/A(4)/A(3)/B(21)
' Zmp' = A1*ke'*U1*50*Pow(10,6)/(fb*PHI'*P*f*kym) B(30)-毎极导体数
700 B(30)=INT(B(30))
710 FOR I=1 TO B(22)
720 B(40+I)=INT(B(30)*B(10+I)/2+.5) 'Zm(i)=Int(Φ*Km(i)%/2+0.5)
730 NEXT I
732 IF K7=2 THEN 740 'If 手动 then740
735 FOR I=1 TO B(22):PRINT "Zm";I;"=?":INPUT B(40+I):NEXT I
'Input Zm(i)
740 B(32)=0 'Zmp=0 B(32)-Zmp
750 FOR I=1 TO B(22)
780 B(32)=B(32)+B(40+I)*2 'Zmp=Zmp+Zm(i)*2
790 NEXT I
800 B(33)=A(4)*B(32)/A1 'Zm=P*Zmp/A1 B(33)-Zm
810 B(61)=0 'kym=0 B(61)-kym
820 FOR I=1 TO B(22)
830 B(50+I)=SIN(B(I)*PI/2/A(35)) 'kym(i)=Sin(Ym(i)*π/(2*Qp1)), B(50+i)-kym(i)
840 B(61)=B(40+I)*B(50+I)+B(61) 'Kym=kym+Zmi*kymi B(61)-kym
850 NEXT I
851 B(61)=B(61)/B(32)*2 'Kym=2*kym/Zmp
854 B(63)=A(1)/A(2)/A(5)/A(6):IF K10<>3 THEN 860
' I1'=P2/(U1*η*cosφ) if(不是电容运转) then860 B(63)-I1'
856 B(63)=B(63)/1.414 'I1'=I1'/1.414
860 IF K7=1 THEN 960 'if (手动) then 960
861 RESTORE 13000
862 FOR J=1 TO 44:READ B1(J):NEXT J '44个线径
863 FOR J=1 TO 44:READ B2(J):NEXT J '44个含绝缘层的线径
864 FOR J=1 TO 44:READ B3(J):NEXT J '44个导线截面积
870 B(65)=SQR(4*A(1)/A(2)/A(5)/A(6)/PI/B(62)/A1):IF K10<>3 THEN 880
' Dm'=SQR(4*P2/(U*η*cosφ)*PI*A1) B(65)-Dm'
875 B(65)=B(65)/SQR(SQR(2)) 'Dm',
880 B(68)=0:N1=1 'Sm=0,
890 FOR I=1 TO 44
900 IF B1(I)>=B(65) THEN 920 'B1(j)>=Dm'
910 NEXT I
915 B(65)=B(65)/SQR(2):N1=2:GOTO 890
920 B(66)=B1(I-1)
930 B(67)=B2(I-1)
940 B(68)=B3(I-1)*N1
950 IF K7=2 THEN 980
960 PRINT "Dm=? Dms=? Sm=? N1=?":INPUT B(66),B(67),B(68),N1
'Input B(66)-Dm, B(67)-Dms, B(68)-Sm, N1
980 B(91)=(A(19)+A(14)/2)*(A(18)+A(37)-.15)+1.57*A(19)^2:IF K2=1 THEN 990
' S1=(R+bs1/2)*(hs2+hs1-0.15)+π*R*R/2 梨形槽 B(91)-S1 If(定子槽为圆底槽) then 990
985 B(91)=(A(14)+A(15))*(A(18)+A(37)-.15)/2
' S1=(bs1+bs2)*(hs2+js1-0.15)/2 平底槽S1
990 B(92)=B(69)*(2*(A(18)+A(37))+PI*A(19)):IF K2=1 THEN 1000
'Si=Ci*(2*(hs2+hs1)+π*R) 圆底槽B(92) -Si If If(定子槽为圆底槽) then 1000
995 B(92)=B(69)*(2*(A(18)+A(37))+A(15))
'Si=Ci*(2*(hs2+hs1)+bs2) 平底槽B(92) -Si
1000 B(93)=B(91)-B(92) 'As= S1-Si B(93)-As
1010 FOR I=1 TO B(22)
1020 B(70+I)=B(40+I)*B(67)^2/B(93)*N1 'Sfm(i)=N1*Zm(i)*Dms*Dms/As B(70+i)-Sfm(i)
1030 B(80+I)=A(31)+4.5*(A(27)+A(16)*2+A(37)*2+A(18)+A(19))*B(I)/A(10)
'Lm(i)=L1+4.5*(Di1+2*hs0+2*hs1+hs2+R)*Ym(i)/Q1 B(80+i)-Lm(i)
'la(i)=L1+4.3*(Di1+2*hs0+2*hs1+hs2+R)*Ya(i)/Q1
1040 NEXT I
1050 B(94)=0 'B(94)-Lm
1060 FOR I=1 TO B(22)
1070 B(94)=B(94)+B(40+I)*B(80+I)/B(32)*2 'Lm= Lm+2*Zm(i)*Lm(i)/Zmp B(94)-Lm
1080 NEXT I
1090 B(95)=2.13*B(94)*B(33)/B(68)/10^4/A1 'R1m=2.13*Lm*Zm/(10000*Sm*A1)
1100 GOTO 2460 '跳转至2460-2930计算副绕组
1150 A4=A(33)*(4.4*A(28)+.75*A(13)) 'T1*(4.4*g+0.75*bo1)
1160 C(1)=A4/(A4-A(13)^2) 'kc1=A4/(A4*A4-Bo*bo), C(1)-Kc1定子卡式系数
1170 C(2)=A(34)*(4.4*A(28)+.75*A(20)):C(2)=C(2)/(C(2)-A(20)^2):C(3)=C(1)*C(2)
' Kc2 = 转子卡式系数 kc=kc1*kc2
1180 C(4)=C(1)*C(2)*A(28) 'ge=Kc1*Kc2*g C(4)-gc有效气隙
1190 C(5)=A(16)/A(13)+2*A(37)/(A(13)+A(14)) 'λu1=hs0/bo1+2*hs1/(bo1+bs1)
1200 C(6)=A(18)/A(15) 'β1=hs2/bs2
1210 C(7)=A(14)/A(15) 'λ1=bs1/bs2
1220 C(8)=1/3-(1-C(7))*(1/4+1/3/(1-C(7))+1/2/(1-C(7))^2)/4-(1/(1-C(7))^3+LOG(C(7))/(1-C(7))^4)/4*(1-C(7))
' Kr11
1225 IF K2=2 THEN 1235 'if 转子为平底槽 then 1235
1230 C(9)=(2*PI^3-9*PI)/1536/C(6)^3+PI/16/C(6)-PI/8/(1-C(7))/C(6)-(PI^2/64/(1-C(7))/C(6)^2+PI/8/(1-C(7))^2/C(6))*LOG(C(7)):GOTO 1240
'Kr12
1235 C(10)=4*C(6)/(1+C(7))^2*C(8):GOTO 1250 'λL1=4*β1* Kr11/Pow((1+λ1),2)
1240 C(10)=C(6)*(C(8)+C(9))/(PI/8/C(6)+(1+C(7))/2)^2
1250 C(11)=C(5)+C(10) 'λs1=λu1+λL1
1260 C(12)=PI*(2*A(3)/100)*(B(33)*B(61)/1000)^2 'Kx=π*(2*f/100)*Pow(Zm*Kym/1000,2)
1270 C(13)=2.51*A(31)*C(11)*C(12)/A(10) 'Xs1=2.51*L1*λs1*Kx/Q1 定子槽漏抗
1280 C(14)=(A(31)+A(12))/2 'L=(L1+L2)/2
1290 C(15)=A(33)-A(13) 'Bt10=T1-bo1
1300 C(16)=A(34)-A(20) 'bt20=T2-bo2
1310 C(17)=.21*C(14)*(C(15)+C(16))^2*C(12)/A(10)/A(28)/(A(33)+A(34))
' Xz=0.21*L*Pow((bt10+bt20),2)*Kx/(Q1*g*()T1+T2) 总曲折漏抗
1320 C(18)=A(11)*C(17)/(A(10)+A(11)) 'Xz1=Q2*Xz/(Q1+Q2) 定子曲折漏抗
1330 C(19)=PI*(A(27)+A0)*B(23)*C(12)/2.54/A(10)/A(4)
' Xe1=π*(Di1+hs)*Ymav*Kx/(2.54*Q1*P) 定子端部漏抗
1340 C(20)=(C(13)+C(18)+C(19)) 'X1=Xs1+Xz1+Xe1 定子漏抗
1350 C(21)=.254*C(14)*A(32)*C(12)/C(4)/A(4)
' Xm0=0.254*L*τp*Kx/(ge*P) 不考虑饱和的激磁漏抗
1360 C(22)=(B(33)*B(61)/1000)^2 'K=Pow(Zm*Kym/1000,2)
1370 C(23)=(8.399999*(A(12)+A(51))/A(48)/A(11))*C(22)+5.3*A(56)*A(58)/A(4)^2/A(53)*C(22):C(23)=C(23)*1
'R2=(8.4*(L2+A1))*K/(S2*Q2)+5.3*Dr*kring*K/(P*P*Sr)
1380 IF A(20)=0 THEN 1392 'If (转子闭口槽) then 1392
1390 C(24)=A(23)/A(20):GOTO 1400 'λu2=hr0-bo2
1392 C(24)=B(33)*B(61)*B(63)*A(6)/A(11) 'λu2'=Zm*Kym*I1'*cosφ/Q2
1394 C(24)=.1+1.12*A(23)/C(24)*10^4 'λu2=0.1+1.12*hr0/(λu2'*10000)
1400 C(25)=A(49)/A(50) 'λ2=br1/br2
1410 C(26)=A(25)/A(50) 'β2=hr2/br2
1420 C(27)=.25+.3333333/(1-C(25))+.5/(1-C(25))^2+1/(1-C(25))^3+LOG(C(25))/(1-C(25))^4
'Kr21'=0.25+1/3/(1-λ2)+0.5/Pow(1-λ2,2)+1/Pow(1-λ2,3)+Log(λ2/Pow(1-λ2,4))
1422 C(27)=1/3-.25*(1-C(25))*C(27):IF K3=2 THEN 1452
' Kr21 = 1/3-0.25*(1-λ2)*Kr21' If (转子平底槽) then 1452
1430 C(28)=.02196/C(26)^3+.19635/C(26)-.3927/(1-C(25))/C(26)-(.1542/(1-C(25))/C(26)^2+.3927/(1-C(25))^2/C(26))*LOG(C(25))
'Kr22 = 0.02196/Pow(β2,3)+0.19635/β2-0.3927/((1-λ2)*β2)-0.1542/((1-λ2)*Pow(β2,2))+0.3927*Log(λ2)/(Pow(1-λ2,2)*β2)
1435 IF K3=3 THEN 1455 'If (转子圆形槽) then 1455
1440 C(29)=PI/4/C(26)*(PI/8/C(26)*(1-C(25)^2)+(1+C(25))/2)^2+(4+3*PI^2)*C(25)^2/32/C(26)^2*(PI/8/C(26)*(1-C(25)^2)+(1+C(25))/2)+(14*PI^2+39*PI)/1536*C(25)^4/C(26)^3
'Kr23=
1450 C(30)=C(26)*(C(27)+C(28)+C(29))/(PI/8/C(26)*(1+C(25)^2)+(1+C(25))/2)^2
' λL2=β2*(Kr21+Kr22+Kr23)/Pow((π*(1+Pow(λ2,2))/(8*β2)+Pow(1+λ2,2), 2) 圆底槽
1451 GOTO 1460
1452 C(30)=2*A(24)/(A(20)+A(49))+4*C(26)*C(27)/(1+C(25))^2:GOTO 1460
' λL2=2*hr1/(bo2+br1)+4*β2*Kr21/Pow(1+λ2,2) 平底槽
1455 C(30)=C(26)*(C(27)+C(28))/(PI/8/C(26)+(1+C(25))/2)^2+2*A(24)/(A(20)+A(49))
' λL2=β2*(kr21+kr22)/(π/8/β2+Pow((1+λ2)/2,2)+2*hr1/(bo2+br1)
1460 C(31)=C(24)+C(30) 'λs2=λu2+λL2
1470 C(32)=2.51*A(12)*C(31)*C(12)/A(11) 'Xs2=2.51*L2*λs2*Kx/Q2 转子槽漏抗
1480 C(33)=A(10)*C(17)/(A(10)+A(11)) 'Xz2=Q1*Xz/(Q1+Q2) 转子槽漏抗
1490 C(34)=C(21)/B(25) 'Xms=Xm0/Ks' 考虑饱和系数的激磁漏抗
1500 C(35)=C(34)+C(20) 'X01=Xms+X1 开路漏抗
1510 C(36)=C(34)/12/C(35)*(PI*A(4)/A(10))^2 'S1=Xms/(12*X01)*Pow(π*P/Q1,2)
1520 C(37)=(C(34)-C(32)-C(33))/C(35) 'S2=(Xms-Xs2-Xz2)/Xm0
1530 C(38)=C(36)^2+4*C(37) 'S12=S1*S1+4S2
1540 C(39)=(SQR(C(38))-C(36))/2 'Kp=(SQR(S1)-S1)/2 漏磁系数
1550 C(40)=C(39)^2 'Kr=Kp*Kp
1560 C(41)=C(39)/(2-C(40)) 'Kg=Kp/(2-Kr)
1620 C(50)=B(29) 'Ke=Ke' 压降系数
1630 C(51)=C(50)*A(2)*50*10^6/B(33)/B(27)/B(61)/A(3)
' PHI= Ke*U1*50*1000000/(Zm*fb*Kym*f) 毎极磁通
1640 C(52)=C(42)*A(31) 'Lfe1=Kfe*L1 定子有效长度
1650 C(53)=C(42)*A(12) 'Lfe2=kfe*L2 转子有效长度
1660 C(54)=A(40)*C(52)*A(35) 'St1=bt1*Lfe1*Qp1 定子齿部截面积
1670 C(55)=A(45)*C(53)*A(36) 'St2=bt2*Lfe2*QP2 转子齿部截面积
1680 C(56)=A(41)*C(52) 'Sc1=hc1*Lfe1 定子轭部截面积
1690 C(57)=A(47)*C(53) 'Sc2=hc2*lfe2 转子轭部截面积
1700 C(58)=A(32)*A(31) 'Sg=τp*L 气息面积
1710 C(59)=C(51)/B(26)/C(54) 'Bt1=PHI/(αi*St1) 定子齿磁密
1720 C(60)=C(51)*C(41)/B(26)/C(55) 'Bt2=PHI*Kg/(αi*St2) 转子齿磁密
1730 C(61)=C(51)*C(40)/B(26)/C(58) 'Bg=PHI*Kr/(αi*Sg) 气隙磁密
1740 C(62)=C(51)/2/C(56) 'Bc1=PHI/(2*Sc1) 定子额部磁密
1750 C(63)=C(51)*C(41)/2/C(57) 'Bc2=PHI*Kg/(2*SC2) 转子额部磁密
1760 RESTORE 11000 '读取11000行的磁密数据,DR510读11000行， DW470读12000行
1770 U=C(59) 'Bt1
1780 GOSUB 7000
1790 C(64)=V 'Sub7000(Bt1)查询att1
1800 RESTORE 11000 '读取11000行的磁密数据,DR510读11000行， DW470读12000行
1810 U=C(60) 'Bt2
1820 GOSUB 7000
1830 C(65)=V 'Sub7000(Bt2)查询att2
1840 C(66)=C(64)*A0 'Att1=att1*hs 定子齿部所需安匝数
1850 C(67)=C(65)*A(46) 'Att2=att2*hr 转子齿部所需安匝数
1860 C(68)=.8*C(4)*C(61) 'Atg=0.8*Ge*Bg 气隙所需安匝数
1870 C(69)=(C(66)+C(67)+C(68))/C(68) 'Ks=(Att1+Att2+Atg)/Atg 齿饱和系数
1900 IF ABS(C(69)-B(25))>.001 THEN 2130 'If Abs(Ks-Ks')>0.001 then 2130
1910 RESTORE 11000 '读取11000行的磁密数据,DR510读11000行， DW470读12000行
1920 U=C(62) 'Bc1
1930 GOSUB 7000
1940 C(70)=V 'atc1=Sub7000(Bc1)
1950 RESTORE 11000 '读取11000行的磁密数据,DR510读11000行， DW470读12000行
1960 U=C(63) 'Bc2
1970 GOSUB 7000
1980 C(71)=V 'atc2=Sub7000(Bc2)
1990 C(72)=(A(26)-A(41))*PI/2/A(4) 'Lc1=(D1-hc1)*π/(2*P)
2000 C(73)=(A(30)+A(47))*PI/2/A(4) 'Lc2=(di2+hc2)*π/(2*P)
2010 RESTORE 21000 '读取21000行的17个额部磁密
2020 U=C(62) 'Bc1
2030 GOSUB 7000
2040 C(43)=V 'Psi1=Sub7000(Bc1) 定子额部修正系数
2050 RESTORE 21000
2060 U=C(63) 'Bc2
2070 GOSUB 7000
2080 C(44)=V 'Psi2=Sub7000(Bc2) 转子额部修正系数
2090 C(74)=C(70)*C(72)*C(43) 'Atc1=atc1*Lc1*Psi1 定子额部所需安匝数
2100 C(75)=C(71)*C(73)*C(44) 'Atc2=atc2*Lc2*Psi2 转子额部所需安匝数
2110 C(76)=(C(66)+C(67)+C(68)+C(74)+C(75))/C(68)
'Ks0=(Att1+Att2+Atg+Atc1+Atc2)/Atg 总饱和系数
2120 GOTO 2210 'goto 2210
2130 B(25)=(B(25)+C(69))/2 'ks'=(Ks'+Ks)/2
2140 RESTORE 10000 '读取21个饱和系数值
2150 U=B(25) 'ks'
2160 GOSUB 7000
2170 B(26)=V 'αi=Sub7000(ks') 极弧系数
2180 GOSUB 7000
2190 B(27)=V 'fb=Sbu7000(αi) 波形系数
2200 N3=N3+1:PRINT "N3=";N3 'N3 += 1, 计算饱和系数用的次数
2205 GOTO 1490 'goto 1490
2210 C(77)=A(4)*C(54)*A0 'Vt1=P*St1*hs 定子齿部面积
2220 C(78)=A(4)*C(56)*C(72)*2 'Vc1=P*Sc1*Lc1*2 定子额部面积
2230 RESTORE 11200 '读取12000行的磁密数据,DR510读112000行， DW470读12200行
2240 U=C(59) 'Bt1
2250 GOSUB 7000
2260 C(79)=V 'pt1=Sub7000(Bt1) 定子齿部单位铁耗
2270 RESTORE 11200 '读取12000行的磁密数据,DR510读112000行， DW470读12200行
2280 U=C(62) 'Bc1
2290 GOSUB 7000
2300 C(80)=V 'pc1=Sub7000(Bc1) 定子额部单位铁耗
2310 C(81)=C(79)*C(77)/10^3 'Pt1=pt1*Vt1/1000 定子齿部铁耗
2320 C(82)=C(80)*C(78)/10^3 'Pc1=pc1*Vc1/1000 定子额部铁耗
2322 RESTORE 11200 '读取12000行的磁密数据,DR510读112000行， DW470读12200行
2324 U=C(60):GOSUB 7000 'Bt2
2326 C(45)=V:RESTORE 11200 '读取12000行的磁密数据,DR510读112000行， DW470读12200行
'pt2=Sub7000(Bt2) 转子齿部单位铁耗
2328 U=C(63):GOSUB 7000 'Bc2
2330 C(46)=V:C(47)=A(4)*C(55)*A(46)*C(45)/10^3
' pc2=Sub7000(Bc2) 转子额部单位铁耗 Pt2=P*St2*hr*pt2/1000 转子齿部铁耗
2332 C(48)=2*A(4)*C(57)*C(73)*C(46)/10^3 'Pc2=2*P*Sc2*lc2*pc2/1000 转子额部铁耗
2334 C(49)=3*C(47)+2.5*C(48):C(90)=2*C(81)+1.7*C(82)
' Pfe2=3*Pt2+2.5*Pc2 转子铁耗 Pfe1=2*Pt1+1.7*Pc1 定子铁耗
2340 C(83)=1.5*(2*C(81)+1.7*C(82)+3*C(47)+2.5*C(48))*(A(3)/50)^1.3
' Pfe=1.5*(2*Pt1+1.7Pc1+3*Pt2+2.5*Pc2)*(Pow(f/50,1.3)
' 铁耗应该去掉转子铁耗Pt2和Pc2, 这样才与理路和实际相符
2350 C(84)=C(21)/C(76) 'Xum=Xm0/ks0 考虑饱和系数的激磁电抗
2360 C(85)=C(84)+C(20) 'X0=Xum+X1 考虑总饱和系数的开路电抗
2370 C(86)=C(36)*C(39)*C(85) 'Xsk=S1*Kp*X0 斜槽漏抗
2380 C(87)=(C(32)+C(33)+C(86))*.8536 'X2=(Xs2+Xz2+Xsk)*0.8536 转子漏抗
2390 C(88)=C(20)+C(87) 'X=X1+X2 总漏抗
2400 C(89)=2*A(2)/(C(85)+C(88)) 'I0=2*U1/(X0+X) 空载电流
2405 E(59)=(6/A(4))^2*(A(26)/10)^4 'Pfw=Pow(6/P,2)*Pow(D1/10,4) 机械耗
2407 C(91)=C(89)^2*B(95)+C(89)^2*C(40)*C(23)/2+C(83)+E(59)
'P0=I0*I0*R1m+I0*I0*Kr*R2/2+Pfe+Pfw 空载损耗
2410 GOTO 3350 'goto 3350
2460 D(28)=0 'Yaav=0 D(28)-Yaav
2470 FOR I=1 TO D(22)
2490 D(28)=D(28)+D(I)*D(10+I) 'Yaav += Ya(i)*ka(i)%
2500 NEXT I
2510 IF K7=1 THEN 2564 'If (手动计算) then 2564
2530 D(24)=D(23)*B(32)*B(61)/D(21)/A1*A2
' Zap' = a*Ymav*Zmp*Kya'*A2/A1
2540 FOR I=1 TO D(22)
2550 D(30+I)=INT(D(24)*D(10+I)/2+.5)
' Za(I)= INT(Zap'*ka(i)%/2+0.5)
2560 NEXT I
2562 IF K7=2 THEN 2570
2564 FOR I=1 TO D(22):PRINT "Za";I;"=?":INPUT D(30+I):NEXT I
'Input Za(i) D(30+i)-Za(i)
2570 S0=0 'S0=0
2580 FOR I=1 TO D(22)
2590 S0=S0+D(30+I) 'S0 += Za(i)
2600 NEXT I
2604 D(91)=D(70+I) 'Sfmax = Sf(i) D(91)-Sfmax
2610 D(27)=S0*2 'Zap=2*S0 D(27)-Zap
2620 D(51)=A(4)*D(27)/A2 'Za=P*Zap/A2 D(51)-Za
2625 S0=0
2630 FOR I=1 TO D(22)
2640 D(40+I)=SIN(D(I)*PI/2/A(35)) 'Kya(i)%= Sin(π*Ya(i)/(2*Qp1))
2650 S0=S0+D(40+I)*D(30+I) 'S0 += Kya(i)*Za(i)
2660 NEXT I
2670 D(52)=S0/D(27)*2:IF K7=1 THEN 2768 'Kya=2*S0/Zap D(52)-Kya if (手动计算) then 2768
2675 IF K10<>3 THEN 2690
2680 D(53)=D(23) 'ka=a
2690 D(54)=B(68)/D(53) 'Sa'=Sm/ka
2700 D(57)=0:N2=1 'Sa=0
2710 FOR I=1 TO 44
2720 IF B3(I)>=D(54) THEN 2740
2730 NEXT I
2735 D(54)=D(54)/SQR(2):N2=2:GOTO 2710
2740 D(57)=B3(I)
2750 D(55)=B1(I)
2760 D(56)=B2(I)*N2
2765 IF K7=2 THEN 2770
2768 PRINT "Da=? Das=? Sa=? N2=?":INPUT D(55),D(56),D(57),N2
'Input D(55)-Da, D(56)-Das, D(57)-Sa, N2
2770 FOR I=1 TO D(22)
2775 D(91)=0 'D(91)-Sfmax=0
2780 D(60+I)=D(30+I)*D(56)^2/B(93)*N2:NEXT I 'Sfa(i)=2*Za(i)*Das*Das/As D(60+i)-Sfa(i)
2782 IF B(1)=A(35) THEN 2786 'If (Ym1=QP1) then 2786
2784 IF B(1)<>A(35) THEN 2792 'if (Ym1<> Qp1) then 2792
2786 D(71)=B(71)*2 'Sf1=2*Sfm1
2788 FOR I=2 TO B(22):D(70+I)=B(70+I)+D(60+B(I)-A(35)/2+I):NEXT I
' I = 2 to Nm, Sf(i)=Sfm(i)+Sfa(j) D(70+i)-Sf(i)
2790 GOTO 2821
2792 FOR I=1 TO B(22):D(70+I)=B(70+I)+D(60+B(I)-A(35)/2+I):NEXT I
'I = 1 to Nm, Sf(i)=Sfm(i)+Sfa(j) D(70+i)-Sf(i)
2821 D(91)=0 'Sfmax=0
2822 FOR I=1 TO D(22)+B(22) ' I = 1 to Nm+Na
2823 IF D(70+I)<D(91) THEN 2825 ' If sf(i) < Sfmax then 2825
2824 D(91)=D(70+I) 'Sfmax = Sf(i)
2825 NEXT I
2840 FOR I=1 TO D(22) ' i = 1 to Na
2850 D(80+I)=A(31)+4.3*(A(27)+A(16)*2+A(37)*2+A(18)+A(19))*D(I)/A(10)
'la(i)=L1+4.3*(Di1+2*hs0+2*hs1+hs2+R)*Ya(i)/Q1
2860 NEXT I
2870 S0=0
2880 FOR I=1 TO D(22) ' i = 1 to Na
2890 S0=S0+D(30+I)*D(80+I) 'So += Za(i)*La(i)
2900 NEXT I
2910 D(92)=S0/D(27)*2 'La = 2*S0 /Zap
2920 D(93)=2.13*D(92)*D(51)/(D(57)*10^4)/A2 'R1a=2.13*La*Za/(Sa*10000*A2) 副绕组电阻
2930 GOTO 1150 '跳回转至1150
2980 E(1)=B(95)+C(40)*C(23) 'Rmst=R1m+Kr*R2 主绕组启动总电阻
2990 E(2)=SQR(E(1)^2+C(88)^2) 'Zmst'=SQR(Rmst*Rmst+X*X) 不记饱和时主绕组启动总阻抗
3000 E(3)=A(2)/E(2) 'Imst'=U1/Zmst' 不记饱和时主绕组启动电流
3010 E(4)=1.615*C(14)*((C(11)-(A(16)+A(37))/A(13))/A(10))+1.615*C(14)*(C(31)-C(24))/A(11)+C(19)/C(12)/2
' λ1s=1.615*L*((λs1-hs1/bo1)/Q1+(λs1-λu2/Q2))+Xe1/Kx/2 不受饱和影响的磁漏导
3030 E(5)=C(88)/C(12)/2-E(4):N6=0 'λ2s=X/(2*Kx)-λ21 受饱和影响的磁漏导
3040 E(6)=4*SQR(2)*E(5)*B(33)*B(61)*E(3) 'PHI0=4*SQR(2)*λ2s*Zm*Kym*Imst' 受饱和影响部分如不作薄荷计算时的漏磁
3050 E(7)=25000*C(14)*(A(33)+A(34)-2*A(13)-2*A(20))^2/4/(A(33)+A(34))+25000*C(14)*(A(16)+A(37)+ A(23))
'Φm=25000*L*(Pow(T1+T2-2bo1-2*bo2,2)/(4*(T1+T2))+hs1+hr0 受饱和影响部分容许漏磁最大值
3060 E(8)=E(6)/E(7) 'f0m=PHI0/Φm
3070 RESTORE 14000 '读取14个f0m值
3080 U=E(8)
3090 GOSUB 7000
3100 E(9)=V 'fem=Sub7000(fom) 饱和漏抗曲线
3110 E(10)=E(5)*E(9)/E(8) 'λ2s'=λ2s*fem/f0m 受饱和影响部分在饱和时的漏磁导
3120 E(11)=2*C(12)*(E(4)+E(10)) 'Xmst'=2*Kx*(λ1s+λ2s') 受饱和影响时的启动总电抗
3130 E(12)=(1+C(23)^2/C(88)/C(85))/(1+(C(23)/C(85))^2) 'Cr=(1+R2*R2/(X*X0))/(1+Pow(R2/X0,2)) 启动漏抗系数
3140 E(13)=E(12)*E(11) 'Xmst=Cr*Xmst' 视在启动总电抗
3150 E(14)=SQR(E(1)^2+E(13)^2) 'Zmst=SQR(Rmst*Rmst+Xmst*Xmst) 主绕组启动总阻抗
3160 E(15)=A(2)/E(14) 'Imst=U1/Zmst 主绕组启动电流
3170 IF ABS(E(15)-E(3))<=.025 THEN 3200 'If Abs(Imst-Imst')<=0.025 then 3200
3180 E(3)=(E(15)+E(3))/2:N6=N6+1:PRINT "N6=";N6 'Imst'=(Imst+Imst')/2 N6+=1 goto 3040
3190 GOTO 3040
3200 E(16)=ATN(E(13)/E(1))*180/3.14 'Qmst=ATN(Xmst/Rmst)*180/3.14 主绕组起动电流相位角
3205 IF K10=1 THEN 3265 'if 电阻启动 then 3265
3210 PRINT "C=":INPUT E(17) 'input C
3220 E(18)=5*10^5/PI/A(3)/E(17) 'Xc=5*100000/(PI*f*C)
'3230 RESTORE 15000 '读取15000行的电容值数据
'3240 U=E(17) 'C
'3250 GOSUB 7000
3260 E(19)=E(18)*0.0022:IF K10<>1 THEN 3268 'Rc=Sub7000(C) if 非电阻启动 then 3268
3265 E(18)=0:E(19)=0 '电阻启动 Xc=0， Rc=0
3268 D(23)=D(51)*D(52)/B(33)/B(61) 'a=Za*Kya/Zm/Kym
3270 E(20)=D(93)+D(23)^2*C(40)*C(23)+E(19) 'Rast=R1a+a*a*Kr*R2+Rc 副绕组起动总电阻
3280 E(21)=D(23)^2*E(13)-E(18) 'Xast=a*a*Xmst-Xc 副绕组起动总电抗
3290 E(22)=SQR(E(20)^2+E(21)^2) 'Zast=SQR(Rast*Rast+Xast*Xast) 副绕组起动总电抗
3300 E(23)=A(2)/E(22) 'Iast=U1/Zast 副绕组起动电流
3310 E(24)=ATN(E(21)/E(20))*180/3.14 'Qast=ATN(Xast/Rast)*180/3.14 副绕组起动电流相位角
3320 E(25)=1.77*D(23)*C(40)*C(23)*E(15)*E(23)*SIN(E(16)*3.14159/180-E(24)*3.1415
9/180)/A(1) 'Tst=1.77*a*Kr*R2*Imst*Iast*Sin(Qmst-Qast)/P2n 起动转矩倍数
3330 E(26)=E(15)/E(22)*SQR((E(1)+E(20))^2+(E(13)+E(21))^2)
' I1st=Imst*SQR(Pow(Rmst+Rast,2)+Pow(Xmst+Xast,2))/Zast 起动电流
3340 E(27)=E(23)*SQR(E(19)^2+E(18)^2) 'Ecst=Iast*SQR(Rc*Rc+Xc*Xc) 电容器起动电压 goto 3790
3345 GOTO 3790
3350 E(28)=C(40)*C(23)/2 'M1=kr*R2/2
3360 E(29)=C(23)/C(85) 'M2=R2/X0
3370 E(30)=C(39)*C(85)/2 'M3=Kp*X0/2
3380 E(31)=C(39)*C(87)/2 'M4=Kp*X2/2
3395 N5=0
3400 E(33)=2-E(32) '2-S'=2-S'
3410 DEF FNA(X)=E(28)/X/(1+(E(29)/X)^2) 'FNA(X)=M1/X/(1+Pow(M2/X,2))
3420 DEF FNB(X)=(E(30)*(E(29)/X)^2+E(31))/(1+(E(29)/X)^2) 'FNB(X)=(M3*Pow(M2/X,2)+M4)/(1+Pow(M2/X,2))
3430 H(1)=E(32) 'H(1)=S'
3440 H(2)=E(33) 'H(2)=2-S'
3450 GOSUB 8000 'goto sub8000
3455 IF K10<>3 THEN 4000 'If 不是电容运转 then 4000
3460 FOR J=3 TO 24
3470 E(31+J)=H(J) 'E(34-43)= H(3-12) set Rf, Rb, Xf, Xb, Rt, Xt, Rta, Xta, Rf-Rb, Xf-Xb
'E(44-55)=H(13-24) set R3, X3, R4, X4, R5, X5, R6, X6, Im, Jm, Qm, P2
3480 NEXT J
3490 E(60)=SQR(E(50)^2+E(51)^2) 'Ia=SQR(R6*R6+X6*X6)
3500 E(56)=E(60)/D(57) 'Ja=Ia/Sa
3510 E(57)=ATN(E(51)/E(50))*180/PI 'θa=ANT(X6/R6)*180/PI
3520 E(58)=SQR((E(48)+E(50))^2+(E(49)+E(51))^2) 'I1=SQR((Pow(R5+R6,2))+Pow(X5+X6,2))
3530 E(61)=ABS((A(1)-E(55))/E(55)) 'P2%=Abs(P2-P2n)/P2
3540 IF E(61)<=.003 THEN 3600 'if P2%<=0.003 then 3600
3550 U9=E(32)-.1*E(61)*SGN(E(55)-A(1)) 'U9=S'-0.1*P2%*SGN(P2-P2n)
3560 E(32)=E(32)/2+U9/2 'S'= s'/2+U9/2
3580 N5=N5+1:PRINT "N5=";N5 ' N5+=1
3590 GOTO 3400
3600 E(62)=A(2)*(E(48)+E(50))+C(83)/2+.025*A(1) 'P1=U1*(R5+R6)+Pfe/2+0.025*P2n
3610 E(63)=E(55)/E(62) 'η=P2/P1
3620 E(64)=COS(ATN((E(49)+E(51))/(E(48)+E(50)))) 'cosφ=COS(ATN(X5+X6)/(R5+R6))
3630 E(65)=(E(52)*B(33)*E(53)+E(60)*E(56)*D(51))/(PI*A(27)) 'AJ=(Im*Zm*Jm+Iz*Za*Ja)/(PI*Di1)
3640 E(66)=(1-E(32))*120*A(3)/A(4) 'n=(1-S')*120*f/P
3642 EX=E(52)*B(95)*SIN(E(54)*3.14/180)+E(52)*C(20)*SIN((90+E(54))*3.14/180)
' EX=Im*R1m*SIN(Qm*3.14/180)+Im*X1*Sin((90+Qm)*3.14/180)
3644 EY=A(2)-E(52)*B(95)*COS((90+E(54))*3.14/180)-E(52)*C(20)*COS(E(54)*3.14/180)
' EY=U1-Im*R1m*Cos((90+Qm)*3.14/180)-Im*X1*Cos(Qm*.3.14/180)
3650 E(67)=SQR(EX^2+EY^2)/A(2):A7=ABS(E(67)-C(50)):IF A7<=.001 THEN 3660
' Ke=SQR(EX*EX+EY*EY)/U1 压降系数 A7=ABS(Ke-ke') if A7<=0.001 then 3660
3652 B(29)=(B(29)+E(67))/2:N4=N4+1:PRINT "N4=";N4:IF K7=2 THEN 680
' ke'=(ke'+ke)/2 N4+=1 if 自动计算 then 680
3654 GOTO 1620
3660 E(68)=E(60)*E(18) 'Ec=Ia*Xc
3670 E(69)=C(23)/C(88) 'R2/X=R2/X
3680 RESTORE 16000 '读取16000行的16个R2/X值
3690 U=E(69)
3700 GOSUB 7000
3710 E(70)=V:E(71)=2-E(70) 'Smax=Sub7000(R2/X), 有R2/X查询最大转矩是的转差率值
3720 H(1)=E(70) 'H(1)=Smax
3730 H(2)=E(71) 'H(2)=2-Smax
3740 GOSUB 8000
3750 FOR J=3 TO 24
3760 E(69+J)=H(J) 'E(72-93)=H(3-24), set Rfm, Rbm, Xfm, Xbm, Rtm, Xtm, Rtam, Xtam, Rfm-Rbm, Xfm-Xbm,
'R3m, X3m,R4m, X4m, R5m, X5m, R6m, X6m, Imm, Jmm, Qmm, P2m
3770 NEXT J
3780 E(94)=(1-E(32))*E(93)/(1-E(70))/E(55) 'Tm=(1-S')*P2m/(1-Smax)/P2
3785 GOTO 2980
3790 E(95)=7.85*A(26)^2*A(31)*C(42)*10^-3 'Gfe=7.85*D1*D1*L*Kfe/1000 硅钢片质量
3800 E(96)=8.899999*B(94)*B(33)*B(68)/10^5 'Gcum=8.9*Lm*Zm*Sm/100000
3810 E(97)=8.899999*D(92)*D(51)*D(57)/10^5 'Gcua=8.9*La*Za*Sa/100000
3820 E(98)=2.7*(A(12)*A(11)*A(48)+2*A(53)*PI*A(56))/10^3 'Gal=2.7*(L2*Q2*S2+2*Sr*PI*Dr)/1000
3821 E(100)=E(96)+E(97) 'Gcu=Gcum+Gcua
3825 PRINT "Bg=";C(61),"Ks=";C(69),"Ke=";E(67),"Sfmax=";D(91),"P1=";E(62),"P2=";E(55),"EAT=";E(63),"COS PHI=";E(64),"AJ=";E(65),"Ec=";E(68),"S=";E(32),"Smax=";E(70),"Qm=";E(54),"Qa=";E(57),"I1=";E(58)
3826 PRINT "Jm=";E(53),"Ja=";E(56),"a=";D(23),"Ka=";D(53),"Tst=";E(25),"Ist=";E(26),"Tmax=";E(94),"P2m=";E(93),"Zmp=";B(32),"Dm=";B(66),"Zap=";D(27),"Da=";D(55)
3827 PRINT "Bg=";C(61),"Bt1=";C(59),"Bt2=";C(60),"Bc1=";C(62),"Bc2=";C(63),"Kso=";C(76)
3828 IF K10<>3 THEN 3833 'IF(非电容运转) then 3833
3829 AA=(C(20)+2*E(36))/(B(95)+2*E(34)):XC=(AA^2+1)*(C(20)+2*E(36)):C=10^6/2/3.14159/A(3)/XC
' AA=a'=(X1+2*Xf)/(R1m+2*Rf) XC=(Aa*a+1)*(X1+2*Xf) C=1000000/(2*PI*f/XC)
3830 PRINT "a'=";AA,"C'=";C:PRINT "c''=":INPUT CJJ 'input C'
3831 XCC=10^6/2/3.14/CJJ/A(3):AJJ=SQR(XCC/(C(20)+2*E(36))-1):PRINT "C''=";CJJ,"a''=";AJJ:STOP
' XCC=C''=1000000/(2*PI*f*C) AJJ=a'=SQR(XCC/(X1+2*Xf)-1)
3832 GOTO 3837 'goto 3837
3833 IF K10=1 THEN 3837 'if 电阻起动 then 3837
3834 RFST=E(28)/(1+E(29)^2):XFST=(E(30)*E(29)^2+E(28))/(1+E(29)^2)
' RFST=M1/(1+M2*M2) XFST=(M3*M2*M2+M1)/(1+M2*M2)
3835 AA=(C(20)+2*XFST)/(B(95)+2*RFST):XC=(AA^2+1)*(C(20)+2*RFST)
' AA=a'=(X1+2XFST)/(R1m+2*RFST) XC=C'=(a'*a'+1)*(X1+2*RFST)
3836 C=10^6/(2*PI*A(3)*XC):PRINT "a'=";AA,"C'=";C:STOP ' C=1000000/(2*PI*f*C')
3837 STOP: PRINT " -----------------------";NA$;"--------------------"
3838 PRINT "Bg=";C(61),"Ks=";C(69),"Ke=";E(67),"Sfmax=";D(91),"P1=";E(62),"P2=";E(55),"EAT=";E(63),"COS PHI=";E(64),"AJ=";E(65),"Ec=";E(68),"S=";E(32),"Smax=";E(70),"Qm=";E(54),"Qa=";E(57),"I1=";E(58)
3839 PRINT "Jm=";E(53),"Ja=";E(56),"a=";D(23),"Ka=";D(53),"Tst=";E(25),"Ist=";E(26),"Tmax=";E(94),"P2m=";E(93),"Zmp=";B(32),"Dm=";B(66),"Zap=";D(27),"Da=";D(55)
3840 PRINT "Bg=";C(61),"Bt1=";C(59),"Bt2=";C(60),"Bc1=";C(62),"Bc2=";C(63),"Kso=";C(76)
3841 PRINT "a'=";AA,"C'=";C,,"a''=";AJJ,"C''=";CJJ:PRINT :PRINT :STOP
3845 FOR J=1 TO 10:PRINT A$(3*J-2);A(3*J-2),A$(3*J-1);A(3*J-1),A$(3*J);A(3*J)
3850 NEXT J :PRINT :PRINT :STOP
3851 FOR J=11 TO 20:PRINT A$(3*J-2);A(3*J-2),A$(3*J-1);A(3*J-1),A$(3*J);A(3*J)
3852 NEXT J :PRINT :PRINT :STOP
3860 FOR J=1 TO 11
3870 PRINT B$(3*J-2);B(3*J-2),B$(3*J-1);B(3*J-1),B$(3*J);B(3*J)
3880 NEXT J:PRINT :PRINT :STOP
3881 FOR J=12 TO 22
3882 PRINT B$(3*J-2);B(3*J-2),B$(3*J-1);B(3*J-1),B$(3*J);B(3*J)
3883 NEXT J:PRINT :PRINT :STOP
3884 FOR J=23 TO 32
3886 NEXT J:PRINT :PRINT :STOP
3890 FOR J=1 TO 10:PRINT C$(3*J-2);C(3*J-2),C$(3*J-1);C(3*J-1),C$(3*J);C(3*J)
3900 NEXT J:PRINT :PRINT :STOP
3901 FOR J=11 TO 20:PRINT C$(3*J-2);C(3*J-2),C$(3*J-1);C(3*J-1),C$(3*J);C(3*J)
3902 NEXT J:PRINT :PRINT :STOP
3903 FOR J=21 TO 31:PRINT C$(3*J-2);C(3*J-2),C$(3*J-1);C(3*J-1),C$(3*J);C(3*J)
3904 NEXT J:PRINT :PRINT :STOP
3910 FOR J=1 TO 11
3920 PRINT D$(3*J-2);D(3*J-2),D$(3*J-1);D(3*J-1),D$(3*J);D(3*J):NEXT J:PRINT :PRINT :STOP
3921 FOR J=12 TO 20
3922 PRINT D$(3*J-2);D(3*J-2),D$(3*J-1);D(3*J-1),D$(3*J);D(3*J):NEXT J:PRINT :PRINT :STOP
3923 FOR J=21 TO 31
3924 PRINT D$(3*J-2);D(3*J-2),D$(3*J-1);D(3*J-1),D$(3*J);D(3*J):NEXT J:PRINT :PRINT :STOP
3930 FOR J=1 TO 11
3940 PRINT E$(3*J-2);E(3*J-2),E$(3*J-1);E(3*J-1),E$(3*J);E(3*J)
3950 NEXT J:PRINT :PRINT :STOP
3951 FOR J=12 TO 22
3952 PRINT E$(3*J-2);E(3*J-2),E$(3*J-1);E(3*J-1),E$(3*J);E(3*J)
3953 NEXT J:PRINT :PRINT :STOP
3954 FOR J=23 TO 34
3955 PRINT E$(3*J-2);E(3*J-2),E$(3*J-1);E(3*J-1),E$(3*J);E(3*J)
3956 NEXT J:PRINT :PRINT
3960 PRINT "I1'=";I1J,"T=";T,"I1m=";I1M,"N1=";N1,"N2=";N2,"am=";A1,"aa=";A2
3970 END
3980 GOTO 2980
4000 FOR J=3 TO 12:E(31+J)=H(J):NEXT J 'E(34-43)=H(3-12) set Rf, Rb, Xf, Xb, Rt, Xt, Rta, Xta, Rf-Rb, Xf-Xb
'非电容运转的参数计算
4010 U8=COS(ATN(E(39)/E(38))):I1J=U8*A(2)/E(38) 'U8=Cos(ATN(Xt/Rt)) I1J=U8*U1/Rt
4020 E(58)=I1J+C(90)*E(38)/A(2)/SQR(E(38)^2+E(39)^2) 'I1=I1J+Pfe1*Rt/U1/SQR(Rt*Rt+Xt*Xt)
4030 E(53)=E(58)/B(68)/A1:E(65)=E(58)*E(53)*B(33)/PI/A(27) 'Jm=I1/Sm/A1, AJ=I1*Jm*Zm/Pi/Di1
4040 E(55)=I1J^2*E(42)*(1-E(32))-(E(59)+C(83)/2) 'P2=I1J^2*(Rf-Rb)*(1-S')-Pfw-Pfe/2
4050 E(61)=ABS(A(1)-E(55))/E(55) 'P2%=(Pn-P2)/P2
4060 IF E(61)<=.002 THEN 4100
4070 U9=E(32)-.1*E(61)*SGN(E(55)-A(1)) 'U9=S'-0.1*P2%*Sgn(P2-Pn)
4080 E(32)=(E(32)+U9)/2 'S'=(S'+U9)/2
4090 N5=N5+1:PRINT "N5=";N5:GOTO 3400
4100 E(62)=A(2)*E(58)*U8+C(90)
4110 E(63)=E(55)/E(62):E(64)=E(62)/A(2)/E(58):E(54)=-ATN(SQR(1-E(64)^2)/E(64))*180/3.14
4112 EX=E(58)*B(95)*SIN(E(54)*3.14/180)+E(58)*C(20)*SIN((90+E(54))*3.14/180)
4114 EY=A(2)-E(58)*B(95)*COS((90+E(54))*3.14/180)-E(58)*C(20)*COS(E(54)*3.14/180)
4120 E(67)=SQR(EX^2+EY^2)/A(2)
4130 A7=ABS(E(67)-B(29)):IF A7<=.001 THEN 4150
4135 B(29)=(E(67)+B(29))/2:N4=N4+1:PRINT "N4=";N4:IF K7=2 THEN 680
4140 GOTO 1620
4150 E(66)=(1-E(32))*120*A(3)/A(4)
4160 T=97.5*E(55)/E(66):E(69)=C(23)/C(88)
4170 RESTORE 16000:U=E(69):GOSUB 7000
4180 E(70)=V:E(71)=2-E(70)
4190 H(1)=E(70):H(2)=E(71):GOSUB 8000
4200 FOR J=3 TO 12:E(69+J)=H(J):NEXT J
4210 I1M=A(2)/SQR(E(76)^2+E(77)^2)
4220 E(93)=I1M^2*E(80)*(1-E(70))-(E(59)+C(49))*(1-E(70))
4230 E(94)=(1-E(32))/(1-E(70))*E(93)/E(55)
4240 GOTO 2980
7000 REM BASIC '插值
7010 READ N
7020 FOR J=1 TO N
7030 READ X(J)
7040 NEXT J
7050 FOR J=1 TO N
7060 READ Y(J)
7070 NEXT J
7080 FOR K=1 TO N-2
7090 IF U<=X(K+1) THEN 7130
7100 NEXT K
7110 K=N-2
7130 IF K=1 THEN 7160
7140 IF U-X(K)>=X(K+1)-U THEN 7160
7150 K=K-1
7160 X0=X(K):X1=X(K+1):X2=X(K+2)
7170 P=(U-X1)*(U-X2)/(X0-X1)/(X0-X2)
7180 Q=(U-X0)*(U-X2)/(X1-X0)/(X1-X2)
7190 R1=(U-X0)*(U-X1)/(X2-X0)/(X2-X1)
7200 V=P*Y(K)+Q*Y(K+1)+R1*Y(K+2)
7210 RETURN
8000 REM SUB
8010 H(3)=FNA(H(1)) 'H(3)=Rf=M1/S'/(1+Pow(M2/S',2)) 正序视在电阻
8020 H(4)=FNA(H(2)) 'H(4)=Rb=M1/(2-S'))/(1+Pow(M2/(1-S')),2) 副序视在电阻
8030 H(5)=FNB(H(1)) 'H(5)=Xf=(M3*Pow(M2/S',2)+M4)/(1+Pow(M2/S',2)),FNB(X)=(M3*Pow(M2/X,2)+M4)/(1+Pow(M2/X,2))
8040 H(6)=FNB(H(2)) 'H(6)=Xb=(M3*Pow(M2/(2-S')),2)/(1+Pow(M2/(2-S')),2)
8050 Z1=H(3)+H(4) 'Z1=Rf+Rb
8060 H(11)=H(3)-H(4) 'H(11)=Rf-Rb
8070 H(12)=H(5)-H(6) 'H(12)=Xf-Xb
8080 H(7)=B(95)+Z1 'H(7)=RT=R1m+Rf+Rb 总等值电阻
8090 H(8)=C(20)+H(5)+H(6):IF K10<>3 THEN 8310 'H(8)=XT=X1+Xf+Xb,总等值电抗,If不是电容运转 then8310
8100 H(9)=D(93)+D(23)*D(23)*Z1+E(19) 'H(9)=Rta=R1a+a*a*(Rf+Rb)+Rc
8110 H(10)=D(23)^2*H(8)-E(18) 'H(10)=Xta=a*a*XT-Xc
8120 H(13)=H(7)*H(9)-H(8)*H(10)-D(23)^2*(H(11)^2-H(12)^2) 'H(13)=R3=RT*Rta-XT*Xta-a*a*Pow(Rf-Rb,2)-Pow(Xf-Xb,2)
8130 H(14)=H(7)*H(10)+H(8)*H(9)-2*D(23)^2*H(11)*H(12) 'H(14)=X3=RT*Xta+Xt*Rta-2*a*a*(Rf-Rb)*(Xf-Xb)
8140 H(15)=A(2)*H(13)/(H(13)^2+H(14)^2) 'H(15)=R4=U1*R3/(R3*R3+X3*X3)
8150 H(16)=H(15)*H(14)/H(13) 'H(16)X4=R4*X3/R3
8160 Z4=H(9)-D(23)*H(12) 'Z4=Rta-a*(Xf-Xb)
8170 Z5=H(10)+D(23)*H(11) 'Z5=Xta+a*(Rf-Rb)
8180 H(17)=H(15)*Z4+H(16)*Z5 'H(17)=R5=R4*(Rta-a*(Xf-Xb))+X4*(Xta+a*(Rf-Rb))
8190 H(18)=H(15)*Z5-H(16)*Z4 'H(18)=X5=R4*(Xta+a*(Rf-Rb))-X4*(Rta-a*(Xf-Xb))
8200 Z6=H(7)+D(23)*H(12) 'Z6=Rt+a*(Xf-Xb)
8210 Z7=H(8)-D(23)*H(11) 'Z7=Xt-a(Rf-Rb)
8220 H(19)=H(15)*Z6+H(16)*Z7 'H(19)=R6=R4*Z6+X4*Z7
8230 H(20)=H(15)*Z7-H(16)*Z6 'H(20)=X6=R4*Z7-X4*Z6
8240 H(21)=SQR(H(17)^2+H(18)^2) '(H21)=Im=SQR(R5*R5+X5*X5)
8250 H(22)=H(21)/B(68) 'H(22)=Jm=Im/Sm
8260 Q5=H(18)/H(17) 'Q5=X5/R5
8270 H(23)=ATN(Q5)*180/PI 'H(23)=θm=ATN(X5/R5)*180/PI
8280 H(24)=H(11)*(H(17)^2+H(18)^2+D(23)^2*(H(19)^2+H(20)^2))
' H(24)=P2=(Rf-Rb)*R5*R5+X5*X5+a*a*(R6*R6+X6*X6)
8290 H(24)=H(24)+2*D(23)*(H(17)*H(20)-H(18)*H(19))*Z1
' P2=P2+2*a*(R5*X6-X5*R6)*(Rf+Rb)
8300 H(24)=H(24)*(1-H(1))-E(59)-.5*C(83) 'P2=P2*(1-S')-Pfw-Pfe/2
8310 RETURN
10000 DATA 21,1.0,1.05,1.1,1.15,1.2,1.25,1.3,1.35,1.4,1.45,1.5,1.55,1.6,1.65,1.7,1.75,1.8,1.85,1.9,1.95,2
' 21个饱和系数值
10020 DATA .635,.649,.66,.671,.682,.691,.7,.711,.72,.726,.733,.74,.745,.753,.76,.765,.77,.776,.782,.788,.792
' 21个极弧系数
10030 DATA 21,1.0,1.05,1.1,1.15,1.2,1.25,1.3,1.35,1.4,1.45,1.5,1.55,1.6,1.65,1.7,1.75,1.8,1.85,1.9,1.95,2
' 21个饱和系数值
10040 DATA 1.11,1.104,1.1,1.097,1.093,1.088,1.086,1.083,1.08,1.077,1.074,1.072,1.069,1.067,1.064,1.063,1.06,1.058,1.056,1.053,1.052
' 21个波形系数
10999 REM 131 Ge Chi Mi, 5000-18000
11000 DATA 131
11001 DATA 5000,5100,5200,5300,5400,5500,5600,5700,5800,5900
11002 DATA 6000,6100,6200,6300,6400,6500,6600,6700,6800,6900
11003 DATA 7000,7100,7200,7300,7400,7500,7600,7700,7800,7900
11004 DATA 8000,8100,8200,8300,8400,8500,8600,8700,8800,8900
11005 DATA 9000,9100,9200,9300,9400,9500,9600,9700,9800,9900
11006 DATA 10000,10100,10200,10300,10400,10500,10600,10700,10800,10900
11007 DATA 11000,11100,11200,11300,11400,11500,11600,11700,11800,11900
11008 DATA 12000,12100,12200,12300,12400,12500,12600,12700,12800,12900
11009 DATA 13000,13100,13200,13300,13400,13500,13600,13700,13800,13900
11010 DATA 14000,14100,14200,14300,14400,14500,14600,14700,14800,14900
11011 DATA 15000,15100,15200,15300,15400,15500,15600,15700,15800,15900
11012 DATA 16000,16100,16200,16300,16400,16500,16600,16700,16800,16900
11013 DATA 17000,17100,17200,17300,17400,17500,17600,17700,17800,17900,18000
11099 REM DR510-50 Chi Hua Qu Xian, 5000-18000
11100 DATA 1.58,1.60,1.62,1.64,1.66,1.69,1.71,1.74,1.76,1.78
11111 DATA 1.81,1.84,1.86,1.89,1.91,1.94,1.97,2.00,2.03,2.06
11112 DATA 2.10,2.13,2.16,2.20,2.24,2.28,2.32,2.36,2.40,2.45
11113 DATA 2.50,2.55,2.60,2.65,2.70,2.76,2.81,2.87,2.93,2.99
11114 DATA 3.06,3.13,3.19,3.26,3.33,3.41,3.49,3.57,3.65,3.74
11115 DATA 3.83,3.92,4.01,4.11,4.22,4.33,4.44,4.56,4.67,4.80
11116 DATA 4.93,5.07,5.21,5.36,5.52,5.68,5.84,6.00,6.16,6.33
11117 DATA 6.52,6.72,6.94,7.16,7.38,7.62,7.86,8.10,8.36,8.62
11118 DATA 8.90,9.20,9.50,9.80,10.1,10.5,10.9,11.3,11.7,12.1
11119 DATA 12.6,13.1,13.6,14.2,14.8,15.5,16.3,17.1,18.1,19.1
11120 DATA 20.1,21.2,22.4,23.7,25.0,26.7,28.5,30.4,32.6,35.1
11121 DATA 37.8,40.7,43.7,46.8,50.0,53.4,56.8,60.4,64.0,67.8
11122 DATA 72.0,76.4,80.8,85.4,90.2,95.0,100,105,110,116,122
11199 REM 131 Ge Chi Mi, 5000-18000
11200 DATA 131
11201 DATA 5000,5100,5200,5300,5400,5500,5600,5700,5800,5900
11202 DATA 6000,6100,6200,6300,6400,6500,6600,6700,6800,6900
11203 DATA 7000,7100,7200,7300,7400,7500,7600,7700,7800,7900
11204 DATA 8000,8100,8200,8300,8400,8500,8600,8700,8800,8900
11205 DATA 9000,9100,9200,9300,9400,9500,9600,9700,9800,9900
11206 DATA 10000,10100,10200,10300,10400,10500,10600,10700,10800,10900
11207 DATA 11000,11100,11200,11300,11400,11500,11600,11700,11800,11900
11208 DATA 12000,12100,12200,12300,12400,12500,12600,12700,12800,12900
11209 DATA 13000,13100,13200,13300,13400,13500,13600,13700,13800,13900
11210 DATA 14000,14100,14200,14300,14400,14500,14600,14700,14800,14900
11211 DATA 15000,15100,15200,15300,15400,15500,15600,15700,15800,15900
11212 DATA 16000,16100,16200,16300,16400,16500,16600,16700,16800,16900
11213 DATA 17000,17100,17200,17300,17400,17500,17600,17700,17800,17900,18000
11299 REM DR510-50 Sun Hao Qu Xian, 5000-18000
11300 DATA 5.15,5.35,5.55,5.76,5.99,6.18,6.38,6.50,6.78,7.00
11301 DATA 7.22,7.43,7.63,7.84,8.03,8.27,8.49,8.70,8.90,9.13
11302 DATA 9.36,9.56,9.76,9.98,10.2,10.4,10.6,10.8,11.0,11.3
11303 DATA 11.5,11.7,12.0,12.2,12.4,12.6,12.8,13.1,13.3,13.5
11304 DATA 13.8,14.0,14.3,14.5,14.8,15.1,15.3,15.6,15.9,16.2
11305 DATA 16.4,16.8,17.1,17.4,17.8,18.1,18.4,18.8,19.2,19.6
11306 DATA 20.0,20.4,20.8,21.2,21.7,22.2,22.6,23.0,23.5,24.0
11307 DATA 24.5,25.0,25.4,26.0,26.4,27.0,27.5,28.0,28.5,29.0
11308 DATA 29.5,30.0,30.5,31.0,31.6,32.1,32.6,33.1,33.6,34.2
11309 DATA 34.7,35.2,35.7,36.2,36.7,37.2,37.8,38.3,38.8,39.4
11310 DATA 39.8,40.4,40.9,41.4,41.9,42.4,42.9,43.5,44.0,44.5
11311 DATA 45.0,45.6,46.1,46.6,47.1,47.7,48.2,48.8,49.2,49.7
11312 DATA 50.2,50.7,51.3,51.8,52.3,52.9,53.5,54.1,54.8,55.5,56.1
11999 REM 131 Ge Chi Mi, 5000-18000
12000 DATA 131
12001 DATA 5000,5100,5200,5300,5400,5500,5600,5700,5800,5900
12002 DATA 6000,6100,6200,6300,6400,6500,6600,6700,6800,6900
12003 DATA 7000,7100,7200,7300,7400,7500,7600,7700,7800,7900
12004 DATA 8000,8100,8200,8300,8400,8500,8600,8700,8800,8900
12005 DATA 9000,9100,9200,9300,9400,9500,9600,9700,9800,9900
12006 DATA 10000,10100,10200,10300,10400,10500,10600,10700,10800,10900
12007 DATA 11000,11100,11200,11300,11400,11500,11600,11700,11800,11900
12008 DATA 12000,12100,12200,12300,12400,12500,12600,12700,12800,12900
12009 DATA 13000,13100,13200,13300,13400,13500,13600,13700,13800,13900
12010 DATA 14000,14100,14200,14300,14400,14500,14600,14700,14800,14900
12011 DATA 15000,15100,15200,15300,15400,15500,15600,15700,15800,15900
12012 DATA 16000,16100,16200,16300,16400,16500,16600,16700,16800,16900
12013 DATA 17000,17100,17200,17300,17400,17500,17600,17700,17800,17900,18000
12099 REM DW470 CHI HUA XI XING,5000-18000,131
12100 DATA 0.65,0.66,0.67,0.68,0.69,0.70,0.71,0.72,0.73,0.74
12101 DATA 0.75,0.755,0.76,0.77,0.78,0.79,0.80,0.81,0.82,0.83
12102 DATA 0.84,0.85,0.86,0.87,0.88,0.89,0.91,0.93,0.94,0.95
12103 DATA 0.96,0.97,0.98,1.00,1.02,1.05,1.06,1.07,1.08,1.10
12104 DATA 1.12,1.14,1.16,1.18,1.20,1.22,1.24,1.26,1.28,1.30
12105 DATA 1.34,1.35,1.40,1.42,1.44,1.49,1.55,1.58,1.60,1.65
12106 DATA 1.70,1.75,1.80,1.85,1.90,1.95,1.98,2.00,2.10,2.20
12107 DATA 2.30,2.40,2.50,2.60,2.70,2.80,2.85,2.90,3.00,3.20
12108 DATA 3.40,3.60,3.80,4.00,4.10,4.30,4.60,4.80,5.00,5.40
12109 DATA 5.90,6.20,6.60,7.00,7.80,8.50,9.00,9.80,10.0,11.0
12110 DATA 12.0,13.0,14.0,15.0,16.0,18.0,20.0,22.0,24.0,25.0
12111 DATA 28.0,30.0,32.0,34.0,38.0,40.0,44.0,47.0,50.0,54.0
12112 DATA 58.0,60.0,64.0,70.0,73.0,76.0,80.0,84.0,90.0,95.0,100.0
12199 REM 131 Ge Chi Mi, 5000-18000
12200 DATA 131
12201 DATA 5000,5100,5200,5300,5400,5500,5600,5700,5800,5900
12202 DATA 6000,6100,6200,6300,6400,6500,6600,6700,6800,6900
12203 DATA 7000,7100,7200,7300,7400,7500,7600,7700,7800,7900
12204 DATA 8000,8100,8200,8300,8400,8500,8600,8700,8800,8900
12205 DATA 9000,9100,9200,9300,9400,9500,9600,9700,9800,9900
12206 DATA 10000,10100,10200,10300,10400,10500,10600,10700,10800,10900
12207 DATA 11000,11100,11200,11300,11400,11500,11600,11700,11800,11900
12208 DATA 12000,12100,12200,12300,12400,12500,12600,12700,12800,12900
12209 DATA 13000,13100,13200,13300,13400,13500,13600,13700,13800,13900
12210 DATA 14000,14100,14200,14300,14400,14500,14600,14700,14800,14900
12211 DATA 15000,15100,15200,15300,15400,15500,15600,15700,15800,15900
12212 DATA 16000,16100,16200,16300,16400,16500,16600,16700,16800,16900
12213 DATA 17000,17100,17200,17300,17400,17500,17600,17700,17800,17900,18000
12299 REM DW470 SUN HAO XU XING,5000-18000,131
12300 DATA 4.84,4.99,5.07,5.15,5.30,5.46,5.62,5.77,5.93,6.08
12301 DATA 6.24,6.40,6.55,6.71,6.86,7.02,7.18,7.33,7.49,7.64
12302 DATA 7.80,7.96,8.11,8.27,8.42,8.58,8.97,9.36,9.44,9.52
12303 DATA 9.59,9.83,10.1,10.3,10.4,10.5,10.8,10.9,11.2,11.4
12304 DATA 11.7,11.9,12.1,12.3,12.5,12.6,12.9,13.3,13.5,13.7
12305 DATA 14.0,14.4,14.8,15.0,15.1,15.4,15.6,16.0,16.4,16.5
12306 DATA 16.8,17.0,17.2,17.4,17.6,17.9,18.3,18.6,18.7,19.1
12307 DATA 19.5,19.9,20.3,20.7,21.1,21.5,21.8,22.2,22.6,23.0
12308 DATA 23.4,23.8,24.2,24.4,24.6,25.0,25.4,25.7,26.0,26.6
12309 DATA 26.7,26.9,27.7,28.1,28.5,28.7,29.3,29.6,30.0,30.4
12310 DATA 31.2,32.0,32.8,33.2,33.5,34.3,34.7,35.1,35.9,36.7
12311 DATA 37.4,37.8,38.2,38.5,38.8,39.0,39.8,40.6,41.3,41.7
12312 DATA 42.1,42.9,43.7,44.5,45.2,46.0,46.8,47.6,48.4,48.8,49.1
13000 DATA .07,.08,.09,.1,.11,.12,.13,.14,.15,.16,.17,.18,.19,.20,.21,.23,.25,.28,.31,.33,.35,.38,.4,.42,.45,.47,.5,.53,.56,.6,.63,.67,.71,.75,.8,.85,.9,.95,1,1.06,1.12,1.18,1.25,1.3
'44个线径B1(44)
13010 DATA .09,.1,.11,.125,.135,.145,.155,.165,.180,.190,.200,.210,.220,.230,.24,.265,.29,.32,.35,.37,.39,.42,.44,.46,.49,.51,.54,.58,.61,.65,.68,.72,.76,.81,.86,.89,.96,1.01,1.08,1.14,1.2,1.26,1.33,1.38
'44个含绝缘层的线径B2(44)
13080 DATA .00385,.00503,.00635,.00785,.0095,.0113,.0133,.0154,.0179,.0201,.0227,.0254,.0283,.0314,.0346,.0415,.0491,.0616,.0755,.0855,.0962,.1134,.1265,.1385,.159,.1735,.1963,.2206,.2463,.2827,.3117,.3525,.3959,.4418,.5026
13090 DATA .5675,.6362,.7088,.7854,.8825,.9856,1.0936,1.2272,1.3273
'44个导线截面积
14000 DATA 14,0,.25,.5,.75,1,1.25,1.5,1.75,2,2.25,2.5,2.75,3,3.25 '14个fom
14020 DATA 0,.3,.5,.66,.77,.85,.905,.97,1,1.03,1.06,1.07,1.08,1.1 '14个饱和漏抗曲线值fem
15000 DATA 17,1,2,3,4,5,6,8,10,75,100,150,200,300,400,600,800,1200
'电容器17个容值
15020 DATA 20,16,15.5,15,14.5,14,13,12,4.3,3.2,2.2,1.6,1.1,.8,.55,.4,.275
'电容器17个电阻值
16000 DATA 16,.1,.2,.3,.4,.5,.6,.7,.8,.9,1,1.1,1.2,1.3,1.4,1.5,1.6 '16个R2/X值
16020 DATA .05,.092,.133,.159,.18,.2,.213,.228,.236,.244,.252,.259,.266,.272,.277,.28 '16个Smax曲线值
18000 DATA AO290-2,2200,230,50,2,.80,.85,27,0.4,1.8,24,18,9.5,.25,.568,.788,.075,0,0.84,.394,0,.25,.171,.02,.25,.45,14.5,7.7,.030,7.60,2.7,9.5,.7,1.70,1.80,.075
'NA$, A(1-30),A(51),A(52),A(44),A(37)
18010 DATA 11,9,7,5,0,0,0,0,0,0,.299,.278,.240,.183,0,0,0,0,0,0,.854,4
'B(1-22)主绕组跨距Ym、Km%,kym, Nm
18020 DATA 11,9,7,5,0,0,0,0,0,0,.518,.482,0,0,0,0,0,0,0,0,.854,4,.025,.95
'D(1-22)副绕组跨距Ya、ka%, kam, Na, B(69)-0.025绝缘层厚度， C(42)-0.95叠压系数
20188 *********************************************
21000 DATA 17,3000,4000,5000,6000,7000,8000,9000,10000,11000
21004 DATA 12000,13000,14000,15000,16000,17000,18000,18700 '17个额部磁密
21006 DATA .65,.645,.64,.63,.62,.6,.58,.55,.52
21008 DATA .475,.42,.365,.33,.3,.28,.265,.25 '额部长度修正系数 Psi1, Psi2
38845 PRINT B$(3*J-2);B(3*J-2),B$(3*J-1);B(3*J-1),B$(3*J);B(3*J)
40000 DATA "P2=","U1=","f=","P=","ETA=","COS FHI=","Ist=","Tst=","Tmax=","Q1=",
"Q2=","L2=","b01=","bs1=","bs2=","hs0=","ALPHA=","hs2=","R=","b02=",
"r1=","r2=","hr0=","hr1=","hr2=","D1=","Di1=","g=","D2=","Di2=",
"L1=","Tau=","T1=" '3, A$(1-33)
40010 DATA "T2=","Qp1=","Qp2=","hs1=","bt11=","bt12=","bt1=","hc1=","bt21=","bt22=",
"B2=","bt2=","hr=","hc2=","S2=","br1=","br2=","A1=","B1=","Sr=",
"Dr0=","DrI=","Dr=","Karing=","Kring=" ' A$(34-58)
40030 DATA "Ym1=","Ym2=","Ym3=","Ym4=","Ym5=","Ym6=","Ym7=","Ym8=","Ym9=","Ym10=",
"Km1%=","Km2%=","Km3%=","Km4%=","Km5%=","Km6%=","Km7%=","Km8%=","Km9%=","Km10%=",
"Kym'=","Nm=","Ymav=","Bg'=","KS'=","ALPHA I=","FB=","PHI'=","Ke'=","Zmp'=" 'B$(1-30)
40040 DATA "Wmp'=","Zmp=","Zm=","**","**","**","**","**","**","**",
"Zm1=","Zm2=","Zm3=","Zm4=","Zm5=","Zm6=","Zm7=","Zm8","Zm9=","Zm10=",
"Kym1=","Kym2=","Kym3=","Kym4=","Kym5=","Kym6=","Kym7=","Kym8=","Kym9=","Kym10=",
"Kym=" 'B$(31-61)
40050 DATA "Jm'=","I1'=","Sm'=","Dm'=","Dm=","Dms=","Sm=","Ci=","**","Sfm1=",
"Sfm2=","Sfm3=","Sfm4=","Sfm5=","Sfm6=","Sfm7=","Sfm8=","Sfm9=","Sfm10=","Lm1=",
"Lm2=","Lm3=","Lm4=","Lm5=","Lm6=","Lm7=","Lm8=","Lm9=","Lm10=","S1=",
"Si=","As=" ' B$(62-93)
40060 DATA "Lm=","R1m=","Kc1=","Kc2=","Kc=","Ge=","LAMDA U1=","BATA 1=","ALPHA 1=","Kr11=",
"Kr12=","LAMDA L1=","LAMDA S1=","Kx=","Xs1=","L=","BT10=","BT20=","Xz=","Xz1=",
"Xe1=","X1=","Xm0=","K=","R2=","LAMDA U2=","ALPHA 2=","BATA 2=" '8, B$(94-95), C$(1-26)
40070 DATA "Kr21=","Kr22=","Kr23=","LAMDA L2=","LAMDA S2=","Xs2=","Xz2=","Xms=","X01=","S1=",
"S2=","S12=","Kp=","Kr=","Kg=","Kfe=","Psi1=","Psi2=","pt2=","pc2=",
"Pt2=","Pc2=","Pfe2=","Ke=","PHI=","Lfe1=","Lfe2=","St1=","St2=","Sc1=",
"Sc2=","Sg=","Bt1=", 'C$(27-60)
40080 DATA "Bg=","Bc1=","Bc2=","att1=","att2=","ATT1=","ATT2=","ATg=","Ks=","atc1=",
"atc2=","Lc1=","Lc2=","ATc1=","ATc2=","Ks0=","Vt1=","Vc1=","pt1=","pc1=",
"Pt1=","Pc1=","Pfe=","Xum=","X0=","Xsk=","X2=","X=","I0=","Pfe1=",
"Po=","Ya1=","Ya2=","Ya3=" 'C$(61-91), D$(1-3)
40090 DATA "Ya4=","Ya5=","Ya6=","Ya7=","Ya8=","Ya9=","Ya10=","Ka1%=","Ka2%=","Ka3%=",
"Ka4%=","Ka5%=","Ka6%=","Ka7%=","Ka8%=","Ka9%=","Ka10%=","Kya'=","Na=","a=",
"Zap'=","Wap'=","Zai=","Zap=","Yaav=","**","**","Za1=","Za2=","Za3=",
"Za4=" 'D$(4-34)
40100 DATA "Za5=","Za6=","Za7=","Za8=","Za9=","Za10=","Kya1=","Kya2=","Kya3=","Kya4=",
"Kya5=","Kya6=","Kya7=","Kya8=","Kya9=","Kya10=","Za=","Kya=","Ka=","Sa'=",
"Da=","Das=","Sa=","**","**","**","Sfa1=","Sfa2=","Sfa3=","Sfa4=",
"Sfa5=" 'D$(35-65)
40110 DATA "Sfa6=","Sfa7=","Sfa8=","Sfa9=","Sfa10=","Sf1=","Sf2=","Sf3=","Sf4=","Sf5=",
"Sf6=","Sf7=","Sf8=","Sf9=","Sf10=","La1=","La2=","La3=","La4=","La5=",
"La6=","La7=","La8=","La9=","La10=","Sfmax=","La=","R1a=","Rmst=","Zmst'=",
"Imst'=" ' D$(66-93) E$(1-3)
40120 DATA "LAMDA 1S=","LAMDA 2S=","PHI 0=","PHI m=","fom=","fem=","LAMDA 2S'=","XMst'=","Cr=","Xmst=",
"Zmst=","Imst=","Qmst=","C=","Xc=","Rc=","Rast=","Xast=","Zast=","Iast=",
"Qast=","Tst=","I1st=","ECst=","M1=","M2=","M3=","M4=","S'=" 'E$(4-32)
40130 DATA "2-S'=","Rf=","Rb=","Xf=","Xb=","RT=","XT=","RTa=","XTa=","Rf-Rb=",
"Xf-Xb=","R3=","X3=","R4=","X4=","R5=","X5=","R6=","X6=","Im=",
"Jm=","Qm=","P2=","Ja=","Qa=","I1=","Pfw=","Ia=","P2%=","P1=",
"ETA=","COS PHI=","AJ=","n=","Ke=" 'E$(33-67)
40140 DATA "EC=","R2/X=","Smax=","2-Smax=","Rfm=","Rbm=","Xfm=","Xbm=","RTM=","XTM=",
"RTAM=","XTAM=","Rfm-Rbm=","Xfm-Xbm=","R3M=","X3M=","R4M=","X4M=","R5M=","X5M=",
"R6M=","X6M=","Imm=","Jmm=","Qmm=","P2m=","Tm=","Gfe=","Gcum=","Gcua=",
"GAL=","IMst%=" 'E$(68-99)
40150 DATA "Gcu=","S=","2-S=","RF=","RB=","XF=","XB=","RT=","XT=","I1'=",
"I1=","JM=","P2=","P2%=","P1=","ETA=","COS PHI=","Ke=","RPM=","T=",
"R2/X=","SM=","2-SM=","RFM=","RBM=","XFM=","SBM=","RTM=","XTM=","I1M=",
"P2M=","TM=" 'E(100)
40160 DATA "END"

复制代码

wenzee · 发表于 2012-9-16 07:19

这个程序，运行环境为微软的 GWBasic.
程序有3个子程序：
三相异步、单相异步、三相发电机

lzemg007 · 发表于 2012-10-9 09:10

wenzee 发表于 2012-9-16 07:17

登录/注册后可看大图

杨心文。
设计单相是一直用这程序。

根本无法运行啊！

wenzee · 发表于 2012-10-13 19:20

lzemg007 发表于 2012-10-9 09:10

登录/注册后可看大图

根本无法运行啊！

说明一下：老的Basic程序是这样运行的：
1. 运行程序如GWBasic, Qbasic,
2. 然后 load“文件名"
3. F5 Run

lzemg007 · 发表于 2012-10-14 09:14

哦，谢谢！我试一试

motors0709 · 发表于 2015-5-18 16:12

wenzee 发表于 2012-9-16 07:17
杨心文。
设计单相是一直用这程序。

楼主可否提供一下完整的basic程序，这些代码运行是错误的

zazhpe · 发表于 2019-7-4 13:44

研究的十分通透啊，程序都懂。

yinduhua · 发表于 2019-7-16 13:35

三相发电机的程序那位发一下

8075555 · 发表于 2019-7-24 09:56

yinduhua 发表于 2019-7-16 13:35
三相发电机的程序那位发一下

嗯，这是一个好提议！

lzemg · 发表于 2019-7-25 10:53

我有这个程序的使用说明

as224041358 · 发表于 2019-11-13 14:57

lzemg 发表于 2019-7-25 10:53
我有这个程序的使用说明

请问怎么使用啊，可以大概讲一下吗？

lzemg · 发表于 2019-11-15 09:08

as224041358 发表于 2019-11-13 14:57
请问怎么使用啊，可以大概讲一下吗？

无法使用，程序有问题

czvmkl · 发表于 2019-11-15 11:28

谁编译好程序吗？

		自动登录	找回密码
密码			立即注册

Flux电机有限元分析教程西莫团购入口 \| 奖励入口	当当网购物入口 \| 奖励入口	欢迎使用！西莫论坛App开放下载	Motor-CAD电机多物理域设计教材购买入口 \| 奖励入口
★新会员论坛须知★	《西莫电机技术》第39期发售火热进行中	Flux电机电磁阀有限元分析教程团购入口 \| 奖励入口	论坛微信公众平台欢迎入驻
西莫电机及相关产品供需交流群开放邀请	★ 论坛VIP会员申请 ★	Motor-CAD.MANATEE电磁热振动噪声教程 \| 奖励入口	西莫团队欢迎您的加盟！
宣传推广合作请联系QQ：25941174	西莫电机论坛微信群正式开放	Flux变压器与电抗器有限元分析团购入口 \| 奖励入口	西莫电机论坛技术版区QQ群汇总

求助福建单相电机程序操作说明

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

相关帖子

评分