电动机双层叠绕电磁计算的VB程序.docx
- 文档编号:12097284
- 上传时间:2023-04-17
- 格式:DOCX
- 页数:33
- 大小:23.30KB
电动机双层叠绕电磁计算的VB程序.docx
《电动机双层叠绕电磁计算的VB程序.docx》由会员分享,可在线阅读,更多相关《电动机双层叠绕电磁计算的VB程序.docx(33页珍藏版)》请在冰豆网上搜索。
电动机双层叠绕电磁计算的VB程序
PrivateSubCommand1_Click()
EndSub
PrivateSubForm_Load()
ConstUN=380
Constf=50
ConstZ1=48
ConstZ2=44
Constπ=3.14159
ConstPN=18500
Constm1=3
Constη1=0.90
ConstCosφ=0.82
Constp=2
Constb01=3.8
Consth01=0.8
Consth11=1.2
Constb11=6.75
Consth21=17
Constr21=4.47
Constδ=0.55
Constb02=1
Consth02=0.5
Consth12=3
Consth22=28
Constr22=1.1
Constb12=2.2
ConstD1=290
ConstDi1=186
ConstDi2=70
ConstNs1=32
Consta1=2
ConstNi1=2
Constd11=15
Constd111=1.18
Constd222=1.3
ConstΔi=0.3
Constlt=200
ConstKFe=0.95
Constβ1=0.833
Consth=2
ConstKU1=0.875
ConstKL1=0.906
Constρcu=8.9*10^3
ConstρFe=7.8*10^3
Constq=4
Const∑s=0.0055
Constρcuw=0.0217*10^-6
α=p*2*π/Z1
UNφ=UN
Ikw=PN/(m1*UNφ)
Zp1=Z1/(2*p)
Zp2=Z2/(2*p)
D2=Di1-2*δ
τ=π*Di1/(2*p)
t1=π*Di1/Z1
t2=π*D2/Z2
bsk=t1
Nφ1=Ns1*Z1/(m1*a1)
N1=Nφ1/2
Ac11=0.985
hss=h21+h11
bt111=π*(Di1+2*h01+2*h11+2*h21)/Z1-2*r21
bt211=π*(Di1+2*h01+2*h11)/Z1-b11
bt11=(bt111+bt211)/2
Ass=(2*r21+b11)/2*(hss-h)+π*r21^2/2
Ai=Δi*(2*hss+π*r21)
Aef=Ass-Ai
sf=Ni1*Ns1*d222^2/Aef
lef=lt+2*δ
lfe=KFe*lt
MMM=Sin(q*α/2)
NNN=q*Sin(α/2)
Kd1=MMM/NNN
Kp1=Sin(π/2*β1)
kdp1=Kd1*Kp1
bt12=π*(D2-2*2/3*(h22+h12+h02))/Z2-b12
kE1=0.954
110:
E1=kE1*UNφ
Kss=1.195
50:
Dimkh(8,1)AsSingle
kh(0,0)=1:
kh(0,1)=1.1095
kh(1,0)=1.2:
kh(1,1)=1.09
kh(2,0)=1.5:
kh(2,1)=1.08
kh(3,0)=1.7:
kh(3,1)=1.072
kh(4,0)=2:
kh(4,1)=1.06
kh(5,0)=2.2:
kh(5,1)=1.057
kh(6,0)=2.5:
kh(6,1)=1.045
kh(7,0)=2.7:
kh(7,1)=1.04
kh(8,0)=3:
kh(8,1)=1.032
Dimk1AsInteger
k1=1
DoWhilek1<=8
IfKss
GoTo13
EndIf
k1=k1+1
Loop
13:
KNm=(Kss-kh(k1,0))*(kh(k1,1)-kh(k1-1,1))/(kh(k1,0)-kh(k1-1,0))+kh(k1,1)
φ=E1/(4*KNm*kdp1*f*N1)
At1=KFe*lt*bt11*Zp1
At2=KFe*lt*bt12*Zp2
hs1=r21+h21+h11+h01
hj11=(D1-Di1)/2-hs1+r21/3
hs2=h02+h12+h22+r22
hj22=(D2-Di2)/2-hs2+r22/3
Aj1=KFe*lt*hj11
Aj2=KFe*lt*hj22
Aδ=τ*lef
Dimxh(12,1)AsSingle
xh(0,0)=1:
xh(0,1)=1.57
xh(1,0)=1.05:
xh(1,1)=1.55
xh(2,0)=1.1:
xh(2,1)=1.53
xh(3,0)=1.15:
xh(3,1)=1.505
xh(4,0)=1.2:
xh(4,1)=1.485
xh(5,0)=1.25:
xh(5,1)=1.466
xh(6,0)=1.3:
xh(6,1)=1.45
xh(7,0)=1.35:
xh(7,1)=1.434
xh(8,0)=1.4:
xh(8,1)=1.42
xh(9,0)=1.45:
xh(9,1)=1.408
xh(10,0)=1.5:
xh(10,1)=1.398
xh(11,0)=1.55:
xh(11,1)=1.38
xh(12,0)=1.6:
xh(12,1)=1.382
DimiAsInteger
i=1
DoWhilei<=11
IfKss
GoTo1
EndIf
i=i+1
Loop
1:
Fs=(Kss-xh(i-1,0))*(xh(i,1)-xh(i-1,1))/(xh(i,0)-xh(i-1,0))+xh(i-1,1)
Bt1=Fs*φ/At1*10^6
Bt2=Fs*φ/At2*10^6
Bj1=φ/(2*Aj1)*10^6
Bj2=φ/(2*Aj2)*10^6
Bδ=Fs*φ/Aδ*10^6
Dimbh(57,1)AsSingle
bh(0,0)=1:
bh(0,1)=3.83
bh(1,0)=1.01:
bh(1,1)=3.92
bh(2,0)=1.02:
bh(2,1)=4.01
bh(3,0)=1.03:
bh(3,1)=4.11
bh(4,0)=1.04:
bh(4,1)=4.22
bh(5,0)=1.05:
bh(5,1)=4.33
bh(6,0)=1.06:
bh(6,1)=4.44
bh(7,0)=1.07:
bh(7,1)=4.56
bh(8,0)=1.08:
bh(8,1)=4.67
bh(9,0)=1.09:
bh(9,1)=4.8
bh(10,0)=1.1:
bh(10,1)=4.93
bh(11,0)=1.11:
bh(11,1)=5.07
bh(12,0)=1.12:
bh(12,1)=5.21
bh(13,0)=1.13:
bh(13,1)=5.36
bh(14,0)=1.14:
bh(14,1)=5.52
bh(15,0)=1.15:
bh(15,1)=5.68
bh(16,0)=1.16:
bh(16,1)=5.84
bh(17,0)=1.17:
bh(17,1)=6
bh(18,0)=1.18:
bh(18,1)=6.16
bh(19,0)=1.19:
bh(19,1)=6.33
bh(20,0)=1.2:
bh(20,1)=6.52
bh(21,0)=1.21:
bh(21,1)=6.72
bh(22,0)=1.22:
bh(22,1)=6.94
bh(23,0)=1.23:
bh(23,1)=7.16
bh(24,0)=1.24:
bh(24,1)=7.38
bh(25,0)=1.25:
bh(25,1)=7.62
bh(26,0)=1.26:
bh(26,1)=7.86
bh(27,0)=1.27:
bh(27,1)=8.1
bh(28,0)=1.28:
bh(28,1)=8.36
bh(29,0)=1.29:
bh(29,1)=8.62
bh(30,0)=1.3:
bh(30,1)=8.9
bh(31,0)=1.31:
bh(31,1)=9.2
bh(32,0)=1.32:
bh(32,1)=9.5
bh(33,0)=1.33:
bh(33,1)=9.8
bh(34,0)=1.34:
bh(34,1)=10.1
bh(35,0)=1.35:
bh(35,1)=10.5
bh(36,0)=1.36:
bh(36,1)=10.9
bh(37,0)=1.37:
bh(37,1)=11.3
bh(38,0)=1.38:
bh(38,1)=11.7
bh(39,0)=1.39:
bh(39,1)=12.1
bh(40,0)=1.4:
bh(40,1)=12.6
bh(41,0)=1.41:
bh(41,1)=13.1
bh(42,0)=1.42:
bh(42,1)=13.6
bh(43,0)=1.43:
bh(43,1)=14.2
bh(44,0)=1.44:
bh(44,1)=14.8
bh(45,0)=1.45:
bh(45,1)=15.5
bh(46,0)=1.46:
bh(46,1)=16.3
bh(47,0)=1.47:
bh(47,1)=17.1
bh(48,0)=1.48:
bh(48,1)=18.1
bh(49,0)=1.49:
bh(49,1)=19.1
bh(50,0)=1.5:
bh(50,1)=20.1
bh(51,0)=1.51:
bh(51,1)=20.2
bh(52,0)=1.52:
bh(52,1)=22.4
bh(53,0)=1.53:
bh(53,1)=23.7
bh(54,0)=1.54:
bh(54,1)=25
bh(55,0)=1.55:
bh(55,1)=26.7
bh(56,0)=1.56:
bh(56,1)=28.5
bh(57,0)=1.57:
bh(57,1)=30.4
DimjAsInteger
j=1
DoWhilej<=56
IfBt1
GoTo11
EndIf
j=j+1
Loop
11:
Ht1=(bh(j,1)-bh(j-1,1))/(bh(j,0)-bh(j-1,0))*(Bt1-bh(j-1,0))+bh(j-1,1)
DimkAsInteger
k=1
DoWhilek<=56
IfBt2
GoTo111
EndIf
k=k+1
Loop
111:
Ht2=(bh(k,1)-bh(k-1,1))/(bh(k,0)-bh(k-1,0))*(Bt2-bh(k-1,0))+bh(k-1,1)
DimlAsInteger
l=1
DoWhilel<=56
IfBj1
GoTo1111
EndIf
l=l+1
Loop
1111:
Hj1=(bh(l,1)-bh(l-1,1))/(bh(l,0)-bh(l-1,0))*(Bj1-bh(l-1,0))+bh(l-1,1)
DimmAsInteger
m=1
DoWhilem<=56
IfBj2
GoTo11111
EndIf
m=m+1
Loop
11111:
Hj2=(bh(m,1)-bh(m-1,1))/(bh(m,0)-bh(m-1,0))*(Bj2-bh(m-1,0))+bh(m-1,1)
Lt1=(h11+h21)+r21/3
Lt2=(h12+h22)+r22/3
Lj11=π*(D1-hj11)/(2*p*2)
Lj22=π*(Di2+hj22)/(2*p*2)
Kδ1=t1*(4.4*δ+0.75*b01)/(t1*(4.4*δ+0.75*b01)-b01^2)
Kδ2=t2*(4.4*δ+0.75*b02)/(t2*(4.4*δ+0.75*b02)-b02^2)
δef=Kδ1*Kδ2*δ
μ0=4*π/10^7
ft1=Ht1*Lt1/10
ft2=Ht2*Lt2/10
Dimch1(3,1)AsSingle
ch1(0,0)=1.24:
ch1(0,1)=0.61
ch1(1,0)=1.395:
ch1(1,1)=0.52
ch1(2,0)=1.55:
ch1(2,1)=0.4
ch1(3,0)=1.71:
ch1(3,1)=0.33
DimeAsInteger
e=1
DoWhilee<=3
IfBj1
GoTo7
EndIf
e=e+1
Loop
7:
Cj1=(ch1(e,1)-ch1(e-1,1))/(ch1(e,0)-ch1(e-1,0))*(Bj1-ch1(e-1,0))+ch1(e-1,1)
Dimch2(3,1)AsSingle
ch2(0,0)=1.24:
ch2(0,1)=0.44
ch2(1,0)=1.395:
ch2(1,1)=0.35
ch2(2,0)=1.55:
ch2(2,1)=0.25
ch2(3,0)=1.71:
ch2(3,1)=0.24
DimgAsInteger
g=1
DoWhileg<=3
IfBj2
GoTo8
EndIf
g=g+1
Loop
8:
Cj2=(ch2(g,1)-ch2(g-1,1))/(ch2(g,0)-ch2(g-1,0))*(Bj2-ch2(g-1,0))+ch2(g-1,1)
Fj1=Cj1*Hj1*Lj11/10
Fj2=Cj2*Hj2*Lj22/10
Fδ=Kδ1*Kδ2*δ*Bδ/(μ0*10^3)
kS=(Fδ+ft1+ft2)/Fδ
If(kS-Kss)/kS>0.01Or(kS-Kss)/kS<(-0.01)Then
Kss=kS-(kS-Kss)/3
GoTo50
EndIf
F0=Fδ+ft1+ft2+Fj1+Fj2
im=2*p*F0/(0.9*m1*N1*kdp1)
Imж=im/Ikw
Xmж=1/Imж
τy=π*(Di1+2*(h01+h11)+h21+r21)*β1/(2*p)*10^-3
lB=lt+2*d11
ss=(b11+2*r21)/(b11+2*r21+2*bt11)
cc=Sqr(1-(ss)^2)
Cs=τy/(2*cc)
lc=lB*10^(-3)+2*Cs
fd=Cs*ss
Cx=0.263*(N1*kdp1)^2*(lef/1000)*PN/(p*UNφ^2)/1000
λU1=h01/b01+2*h11/(b01+b11)
bb12=0.72
hb2=2.17
Dimhhh1(1,0)AsSingle
hhh1(0,0)=0.7:
hhh1(1,0)=0.8
If(hhh1(0,0) Dimhhh2(6,1)AsSingle hhh2(0,0)=2: hhh2(0,1)=1.1 hhh2(1,0)=2.1: hhh2(1,1)=1.15 hhh2(2,0)=2.2: hhh2(2,1)=1.2 hhh2(3,0)=2.3: hhh2(3,1)=1.25 hhh2(4,0)=2.4: hhh2(4,1)=1.3 hhh2(5,0)=2.5: hhh2(5,1)=1.35 hhh2(6,0)=2.6: hhh2(6,1)=1.4 DimlK1AsInteger lK1=1 DoWhilelK1<=5 If(hhh2(lK1,0) y5=(hhh2(lK1+1,1)-hhh2(lK1,1))/(hhh2(lK1+1,0)-hhh2(lK1,0))*(hb2-hhh2(lK1,0))+hhh2(lK1,1) EndIf lK1=lK1+1 Loop Dimhhh3(6,1)AsSingle hhh3(0,0)=2: hhh3(0,1)=0.98 hhh3(1,0)=2.1: hhh3(1,1)=1.02 hhh3(2,0)=2.2: hhh3(2,1)=1.06 hhh3(3,0)=2.3: hhh3(3,1)=1.11 hhh3(4,0)=2.4: hhh3(4,1)=1.16 hhh3(5,0)=2.5: hhh3(5,1)=1.2 hhh3(6,0)=2.6: hhh3(6,1)=1.24 DimlK2AsInteger lK2=1 DoWhilelK2<=5 If(hhh3(lK2,0) y6=(hhh3(lK2+1,1)-hhh3(lK2,1))/(hhh3(lK2+1,0)-hhh3(lK2,0))*(hb2-hhh3(lK2,0))+hhh3(lK2,1) EndIf lK2=lK2+1 Loop λL1=(y6-y5)/(hhh1(1,0)-hhh1(0,0))*(bb12-hhh1(0,0))+y5 EndIf AAA=KU1*λU1 BBB=KL1*λL1 λs1=AAA+BBB lt11=lt Xs1ж=2*m1*p*lt11*λs1*Cx/(Z1*kdp1^2*lef) Xδ1ж=m1*τ*∑s*Cx/(π^2*δef*kdp1^2*kS) XE1ж=1.2*(d11+0.5*fd)*Cx/lef Xσ1ж=Xs1ж+Xδ1ж+XE1ж λU2=h02/b02 λL2=2*h12/(b02+b12)+2.96 λs2=λU2+λL2 lt22=lt Xs2ж=2*m1*p*lt22*λs2*Cx/(Z2*lef) Dimhjj2(12,1)AsSingle hjj2(0,0)=3: hjj2(0,1)=0.097: hjj2(1,0)=4: hjj2(1,1)=0.053 hjj2(2,0)=5: hjj2(2,1)=0.034: hjj2(3,0)=6: hjj2(3,1)=0.025 hjj2(4,0)=7: hjj2(4,1)=0.017: hjj2(5,0)=8: hjj2(5,1)=0.013 hjj2(6,0)=9: hjj2(6,1)=0.01: h
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电动机 双层 电磁 计算 VB 程序