📄 form1.frm
字号:
Begin VB.Label Label18
Caption = "Nt1"
Height = 375
Left = 4320
TabIndex = 34
Top = 2160
Width = 495
End
Begin VB.Label Label17
Caption = "a1"
Height = 375
Left = 3000
TabIndex = 32
Top = 2160
Width = 495
End
Begin VB.Label Label16
Caption = "Ar"
Height = 375
Left = 1680
TabIndex = 30
Top = 2160
Width = 375
End
Begin VB.Label Label15
Caption = "Dr"
Height = 375
Left = 360
TabIndex = 28
Top = 2160
Width = 375
End
Begin VB.Label Label14
Caption = "y1"
Height = 375
Left = 8520
TabIndex = 26
Top = 1440
Width = 375
End
Begin VB.Label Label13
Caption = "Tst"
Height = 495
Left = 7080
TabIndex = 24
Top = 1440
Width = 495
End
Begin VB.Label Label12
Caption = "Tm"
Height = 495
Left = 5760
TabIndex = 22
Top = 1440
Width = 495
End
Begin VB.Label Label11
Caption = "Ist"
Height = 495
Left = 4320
TabIndex = 20
Top = 1440
Width = 375
End
Begin VB.Label Label10
Caption = "cosΦ"
Height = 495
Left = 2760
TabIndex = 18
Top = 1440
Width = 495
End
Begin VB.Label Label9
Caption = "η"
Height = 375
Left = 1680
TabIndex = 16
Top = 1440
Width = 495
End
Begin VB.Label Label8
Caption = "lt"
Height = 255
Left = 360
TabIndex = 14
Top = 1440
Width = 495
End
Begin VB.Label Label7
Caption = "δ"
Height = 375
Left = 8520
TabIndex = 12
Top = 720
Width = 495
End
Begin VB.Label Label6
Caption = "Di1"
Height = 495
Left = 7080
TabIndex = 10
Top = 720
Width = 375
End
Begin VB.Label Label5
Caption = "D2"
Height = 375
Left = 5760
TabIndex = 8
Top = 720
Width = 255
End
Begin VB.Label Label4
Caption = "D1"
Height = 375
Left = 4320
TabIndex = 6
Top = 720
Width = 615
End
Begin VB.Label Label3
Caption = "P"
Height = 375
Left = 3000
TabIndex = 4
Top = 720
Width = 495
End
Begin VB.Label Label2
Caption = "UN"
Height = 375
Left = 1680
TabIndex = 2
Top = 720
Width = 495
End
Begin VB.Label Label1
Caption = "PN"
Height = 255
Left = 360
TabIndex = 0
Top = 720
Width = 375
End
End
Attribute VB_Name = "三相异步电动机验算程序"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Option Explicit
Public Sub cmdinput_Click()
PN = Val(Text1): UN = Val(Text2): p = Val(Text3): D1 = Val(Text4): Di2 = Val(Text5): Di1 = Val(Text6): δ = Val(Text7)
lt = Val(Text8): η1 = Val(Text9): cosΦ1 = Val(Text10): Ist2 = Val(Text11): Tm1 = Val(Text12): Tst1 = Val(Text13): Y1 = Val(Text14)
DR = Val(Text15): AR = Val(Text16): a = Val(Text17): Nt1 = Val(Text18): Nc = Val(Text19): H = Val(Text20): Δi = Val(Text21)
d11 = Val(Text22): Δ = Val(Text23): d = Val(Text24): d12 = Val(Text25): Ps1 = Val(Text43): Z1 = Val(Text44): Z2 = Val(Text45)
b01 = Val(Text26): h01 = Val(Text27): b11 = Val(Text28): h11 = Val(Text29): h21 = Val(Text30): α1 = Val(Text31): r21 = Val(Text32)
b02 = Val(Text33): h02 = Val(Text34): b12 = Val(Text35): b22 = Val(Text36): b32 = Val(Text37): b42 = Val(Text38): α2 = Val(Text39)
h12 = Val(Text40): h22 = Val(Text41): h32 = Val(Text42) '输入数值
J11 = 5: f = 50: KI = 0.95: KFe = 0.95: KNm = 1.1: ρ1 = 0.0000000217: ρ2 = 0.0000000434: ρcu = 8900
ρFe = 7800: c = 1.05: δ1 = 5: KB = 1.04: Pi = 3.1415926: μ0 = 0.4 * Pi / 1000000: m1 = 3: k1 = 2.5: k2 = 2
Picoutput.Print "η'="; η1, "cosΦ'="; cosΦ1, "Ist'="; Ist2, "Tm'="; Tm1, "Tst'="; Tst1, vbCrLf
Rem (一)额定数据和主要尺寸
IKW = PN * 1000 / m1 / UN
q1 = Z1 / 2 / m1 / p: Zp1 = Z1 / 2 / p: Zp2 = Z2 / 2 / p
lef = lt + 2 * δ: D2 = Di1 - 2 * δ: τ = Pi * Di1 / 2 / p: KE1 = 0.0108 * Log(PN) - 0.013 * p + 0.931
bt11 = Pi * (Di1 + 2 * h01 + 2 * h11 + 2 * h21) / Z1 - 2 * r21
bt12 = Pi * (Di1 + 2 * h01 + 2 * h11) / Z1 - b11 'bt11表示定子齿上部宽度，bt12表示定子齿下部宽度
bc1 = (bt11 + bt12) / 2
bt211 = Pi * (D2 - 2 * (h02 + h12)) / Z2 - b12
bt212 = Pi * (D2 - 2 * (h02 + h12 + h22)) / Z2 - b22 'bt211表示转子槽口上部齿宽，bt212表示槽口下部齿宽
If bt211 >= bt212 Then
bc21 = Pi * (D2 - 4 / 3 * (h02 + h12 + h22)) / Z2 - (2 * b22 + b12) / 3
Else: bc21 = Pi * (D2 - 2 / 3 * (h02 + h12 + h22)) / Z2 - (b22 + 2 * b12) / 3
End If
bt221 = Pi * (D2 - 2 * (h02 + h12 + h22)) / Z2 - b32
bt222 = Pi * (D2 - 2 * (h02 + h12 + h22 + h32)) / Z2 - b42 'bt221表示转子槽底上部齿宽，bt222表示槽底下部齿宽
If bt221 >= bt222 Then
bc22 = Pi * (D2 - 2 * (h02 + h12 + h22 + 2 / 3 * h32)) / Z2 - (2 * b42 + b32) / 3
Else: bc22 = Pi * (D2 - 2 * (h02 + h12 + h22 + 1 / 3 * h32)) / Z2 - (b42 + 2 * b32) / 3
End If '计算齿宽
t1 = Pi * Di1 / Z1: t2 = Pi * D2 / Z2: β = Y1 / Zp1: bsk = t1: Ns1 = 2 * Nc: NΦ1 = Ns1 * Z1 / m1 / a
N1 = NΦ1 / 2: d = d11 + Δ: Ac1 = Ac(d11) '调用自定义函数1求导线截面积
As1 = (2 * r21 + b11) / 2 * (h11 + h21 - H) + Pi * r21 ^ 2 / 2
Ai = Δi * (2 * (h11 + h21) + Pi * r21 + 2 * r21 + b11): Aef = As1 - Ai
sf = Nt1 * Ns1 * d ^ 2 / Aef '计算槽满率
Dim m As Integer
If sf > 0.8 Or sf < 0.7 Then
m = MsgBox("槽满率不合格", vbOKOnly & vbExclamation, "槽满率")
End If
α = p * 2 * Pi / Z1: Kd1 = Sin(q1 * α / 2) / q1 / Sin(α / 2): Kp1 = Sin(β * Pi / 2)
Kdp1 = Kd1 * Kp1
At1 = KFe * lt * bc1 * Zp1 / 1000000: At21 = KFe * lt * bc21 * Zp2 / 1000000: At22 = KFe * lt * bc22 * Zp2 / 1000000
hj11 = (D1 - Di1) / 2 - (h11 + h21 + h01 + r21) + r21 / 3
hj21 = (D2 - Di2) / 2 - (h02 + h12 + h22 + h32)
Aj1 = KFe * lt * hj11 / 1000000
Aj2 = KFe * lt * hj21 / 1000000: Aδ = τ * lef / 1000000 '计算齿轭部截面积和空气隙面积
Lt1 = (h11 + h21) + r21 / 3: Lt21 = h12 + h22: Lt22 = h32: Lj1 = Pi * (D1 - hj11) / 4 / p
Lj2 = Pi * (Di2 + hj21) / 4 / p: 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)
kδ = kδ1 * kδ2: δef = kδ * δ '计算齿轭部和空气隙磁路长度
sinα = (b11 + 2 * r21) / (b11 + 2 * r21 + 2 * bc1): cosα = Sqr(1 - sinα ^ 2)
τy = Pi * (Di1 + 2 * (h01 + h11) + h21 + r21) * β / 2 / p: Cs = τy / 2 / cosα
lb = lt + 2 * d12: lc = lb + 2 * Cs: fd = Cs * sinα
Cx = 4 * Pi * f * μ0 * (N1 * Kdp1) ^ 2 * lef * PN / (m1 * p * UN ^ 2)
'计算漏抗系数
λu1 = h01 / b01 + 2 * h11 / (b01 + b11): Call Lou1 '调用过程Lou1()求λL1(自定义函数5)
Ku1 = (3 * β + 1) / 4: KL1 = (9 * β + 7) / 16: λs1 = Ku1 * λu1 + KL1 * λL1
Xs1 = 2 * m1 * p * lt * λs1 * Cx / (Z1 * Kdp1 ^ 2 * lef) '计算定子槽漏抗
XE1 = 1.2 * (d12 + 0.5 * fd) * Cx / lef '计算定子端部漏抗
λu2 = h02 / b02: Call Lou2 '调用过程Lou2()求λL2
λs2 = λu2 + λL2
Xs2 = 2 * m1 * p * lt * λs2 * Cx / Z2 / lef '计算转子槽漏抗
XE2 = 0.7569 * DR * Cx / (lef * 2 * p) '计算转子端部漏抗
R1 = ρ1 * 2 * N1 * lc * 1000 / (Nt1 * Ac1 * a): R11 = R1 * IKW / UN
Gcu = c * lc / 1000 * Ns1 * Z1 * Ac1 / 1000000 * Nt1 * ρcu: GFe = KFe * lt / 1000 * (D1 + δ1) ^ 2 / 1000000 * ρFe
RB = ρ2 * KB * lt * 4 * m1 * 1000 * (N1 * Kdp1) ^ 2 / AB / Z2: RB1 = RB * IKW / UN
RR = ρ2 * DR * 4 * m1 * 1000 * (N1 * Kdp1) ^ 2 / (AR * 2 * Pi * p ^ 2): RR1 = RR * IKW / UN
R22 = RB1 + RR1 '计算定转子电阻和材料重量
Rem (二)磁路计算
Do While (1) '计算满载电势标幺值(循环1)
E1 = KE1 * UN
Ks = 0: Ks1 = 3.6
Do '计算饱和系数(循环2)
Ks1 = Ks - (Ks - Ks1) / 3
Call Fs1 '调用自定义函数2求波幅系数
αp = 1 / Fs: Φ = E1 / (4 * KNm * Kdp1 * f * N1): Bt1 = Fs * Φ / At1
Bt21 = Fs * Φ / At21: Bt22 = Fs * Φ / At22: Bj1 = Φ / Aj1 / 2: Bj2 = Φ / Aj2 / 2: Bδ = Fs * Φ / Aδ
Ht1 = cihua(Bt1): Ht21 = cihua(Bt21): Ht22 = cihua(Bt22)
Hj1 = cihua(Bj1): Hj2 = cihua(Bj2) '调用自定义函数3查磁化曲线
Ft1 = Ht1 * Lt1 / 10: Ft21 = Ht21 * Lt21 / 10: Ft22 = Ht22 * Lt22 / 10: Fδ = δef * Bδ / 1000 / μ0
Ks = (Fδ + Ft1 + Ft21 + Ft22) / Fδ: error2 = Abs(Ks - Ks1) / Ks
Loop While (error2 > 0.0005) '*********循环2结束
If Bt1 > 1.8 Or Bt21 > 1.8 Or Bt22 > 1.8 Or Bj1 > 1.8 Or Bj2 > 1.8 Or Bδ > 1 Then
m = MsgBox("磁密过高", 0 & vbExclamation, "磁通密度")
End If
If Ks > 1.55 Then
m = MsgBox("饱和系数过高", 0 & vbExclamation, "饱和系数")
End If
Cj1 = Cjs(p, Bj1, hj11 / τ) '调用自定义函数41
Cj2 = Cjr(p, Bj2, hj21 / τ) '调用自定义函数42查轭部校正系数Cj
Fj1 = Cj1 * Hj1 * Lj1 / 10: Fj2 = Cj2 * Hj2 * Lj2 / 10: F0 = Fδ + Ft1 + Ft21 + Ft22 + Fj1 + Fj2
Im = 2 * p * F0 / (0.9 * m1 * N1 * Kdp1): Im1 = Im / IKW: Xm = 1 / Im1
Rem (三)参数计算
Σs = s(q1, β) '调用自定义函数6求Σs
Xδ1 = m1 * τ * Σs * Cx / (Pi ^ 2 * δef * Kdp1 ^ 2 * Ks): Xσ1 = Xs1 + XE1 + Xδ1
ΣR = R(p, Z2) '调用自定义函数8求ΣR
Xδ2 = m1 * τ * ΣR * Cx / (Pi ^ 2 * δef * Ks): Xsk = 0.5 * (bsk / t2) ^ 2 * Xδ2
Xσ2 = Xs2 + Xδ2 + XE2 + Xsk: Xσ = Xσ1 + Xσ2 '计算漏抗
Rem (四)工作性能计算
I1p = 1 / η1
L3: '计算效率(循环3)
σ1 = 1 + Xσ1 / Xm: Ix = σ1 * Xσ * I1p ^ 2 * (1 + (σ1 * Xσ * I1p) ^ 2): I1Q = Im1 + Ix
KE = 1 - (I1p * R11 + I1Q * Xσ1): error1 = Abs(KE - KE1) / KE: KE1 = KE - (KE - KE1) / 8
If error1 <= 0.00025 Then Exit Do
Loop '*********循环1结束
KE0 = 1 - Im1 * Xσ1: K = KE0 / KE
Bt10 = K * Bt1: Bj10 = K * Bj1
I11 = Sqr(I1p ^ 2 + I1Q ^ 2): I1 = I11 * IKW: J1 = I1 / (a * Nt1 * Ac1): A1 = m1 * NΦ1 * I1 * 1000 / Pi / Di1
I21 = Sqr(I1p ^ 2 + Ix ^ 2): I2 = I21 * IKW * m1 * NΦ1 * Kdp1 / Z2: IR = I2 * Z2 / (2 * Pi * p)
JB = I2 / AB: JR = IR / AR '计算定转子电流及电密
Pcu11 = I11 ^ 2 * R11: Pcu1 = Pcu11 * PN * 1000: PAl21 = I21 ^ 2 * R22: PAl2 = PAl21 * PN * 1000
Ps = Ps1 * PN * 1000: Pfw = (3 / p) ^ 2 * (D1 / 1000) ^ 4 * 10000: Pfw1 = Pfw / PN / 1000
Gt = 2 * p * At1 * Lt1 * ρFe / 1000: Gj = 4 * p * Aj1 * Lj1 * ρFe / 1000
Phet = p1(Bt10): Phej = p1(Bj10) '调用自定义函数9查轭部损耗曲线
PFet = Phet * Gt: PFej = Phej * Gj: PFe = k1 * PFet + k2 * PFej: PFe1 = PFe / PN / 1000
Σp = Pcu11 + PAl21 + Ps1 + Pfw1 + PFe1: PN1 = 1 + Σp '计算各种损耗
η = 1 - Σp / PN1: error3 = Abs(η - η1) / η
If error3 > 0.00025 Then
η1 = η + (η - η1) / 5: GoTo L3
End If '*********循环3结束
cosΦ = I1p / I11: PFer = ((k1 - 1) * PFet + (k2 - 1) * PFej) / PN / 1000
Sn = PAl21 / (1 + PAl21 + PFer + Ps1 + Pfw1): nN = 60 * f * (1 - Sn) / p
Tm = (1 - Sn) / 2 / (R11 + Sqr(R11 ^ 2 + Xσ ^ 2))
Rem (五)起动过程
hB = h12 + h22 + h32: ξ = 0.001987 * hB / 1000 * Sqr(f / ρ2)
Ka = Ka1(b32 / b22) '调用自定义函数12求截面宽度突变系数Ka
KF = KFx(ξ, 1): Kx = KFx(ξ, 2) '调用自定义函数11查转子肌肤效应系数KF,Kx
hpx = (h12 + h22 + h32) * Kx * Ka
If hpx > (h12 + h22) Then
bpx = b42 + (b32 - b42) * (h12 + h22 + h32 - hpx) / h32: b42 = bpx
hx = hpx - (h12 + h22): h32 = hx
Else
hx1 = hpx - h12: h22 = hx1
bpx1 = b12 + (b22 - b12) * hx1 / h12: b22 = bpx1
End If
Call Lou2: λL2st = λL2 '调用过程Lou2()求λL2st
Ist1 = 24 * Tm * IKW: Ist = 0
Do '计算起动电流(循环4)
Ist1 = Abs(Ist + (Ist - Ist1) / 8)
Fst = Ist1 * Ns1 / a * 0.707 * (Ku1 + Kd1 ^ 2 * Kp1 * Z1 / Z2) * Sqr(KE0)
βc = 0.64 + 2.5 * Sqr(δ / (t1 + t2)): BL = μ0 * Fst / (2 * δ / 1000 * βc)
Kz = Kz1(BL) '调用自定义函数10查起动时漏磁饱和系数Kz
Cs1 = (t1 - b01) * (1 - Kz): Cs2 = (t2 - b02) * (1 - Kz)
Δλu1 = (h01 + 0.58 * h11) / b01 * Cs1 / (Cs1 + 1.5 * b01)
λs1st = Ku1 * (λu1 - Δλu1) + KL1 * λL1
Xs1st = λs1st * Xs1 / λs1: Xδ1st = Kz * Xδ1
Xσ1st = Xs1st + Xδ1st + XE1 '起动时定子漏抗
Δλu2 = h02 * Cs2 / b02 / (Cs2 + b02)
λs2st = (λu2 - Δλu2) + λL2st
Xs2st = λs2st * Xs2 / λs2 '起动时转子槽漏抗
Xδ2st = Kz * Xδ2: Xskst = Kz * Xsk
Xσ2st = Xs2st + Xδ2st + XE2 + Xskst '起动时转子漏抗
Xσst = Xσ1st + Xσ2st: R2st = KF * RB1 + RR1
Rst = R11 + R2st: Zst = Sqr(Rst ^ 2 + Xσst ^ 2) '起动时总阻抗
Ist = IKW / Zst: error4 = Abs(Ist - Ist1) / Ist
Loop While (error4 > 0.0015) '*********循环4结束
Ist3 = Ist / I1
Tst = R2st * (1 - Sn) / Zst ^ 2
η = Format(η, "##.####"): cosΦ = Format(cosΦ, "##.####"): Ist3 = Format(Ist3, "##.####")
Tm = Format(Tm, "##.####"): Tst = Format(Tst, "##.####") '格式化输出
Picoutput.Print "η="; η, "cosΦ="; cosΦ, "Ist*="; Ist3, "Tm*="; Tm, "Tst*="; Tst, vbCrLf
'输出性能指标
sf = Format(sf, "#.####"): Ks = Format(Ks, "#.####"): Bt1 = Format(Bt1, "#.####"): Bt21 = Format(Bt21, "#.####")
Bt22 = Format(Bt22, "#.####"): Bj1 = Format(Bj1, "#.####"): Bj2 = Format(Bj2, "#.####"): Bδ = Format(Bδ, "#.####")
Im1 = Format(Im1, "#.####"): Ix = Format(Ix, "#.####"): I1Q = Format(I1Q, "#.####"): KE = Format(KE, "#.####")
R11 = Format(R11, "#.####"): R22 = Format(R22, "#.####"): Cx = Format(Cx, "#.####"): Xσ1 = Format(Xσ1, "#.####")
Xσ2 = Format(Xσ2, "#.####"): Σp = Format(Σp, "#.####"): nN = Format(nN, "#.####"): J1 = Format(J1, "#.####")
Picoutput.Print "sf="; sf, "ks="; Ks, "R1*="; R11, "R2*="; R22, "Bt1="; Bt1, vbCrLf
Picoutput.Print "Bt21="; Bt21, "Bt22="; Bt22, "Bj1="; Bj1, "Bj2="; Bj2, "Bδ="; Bδ, vbCrLf
Picoutput.Print "Cx="; Cx, "Xσ1="; Xσ1, "Xσ2="; Xσ2, "Im*="; Im1, "Ix*="; Ix, vbCrLf
Picoutput.Print "I1Q*="; I1Q, "KE="; KE, "Σp*="; Σp, "nN="; nN, "J1="; J1, vbCrLf
End Sub
Private Sub Lou2() 'λL2子过程
AB1 = (b02 + b12) / 2 * h12: AB2 = (b12 + b22) / 2 * h22: AB3 = (b32 + b42) / 2 * h32
AB23 = AB2 + AB3: AB = AB1 + AB23
Kr11 = Kr(b12 / b02, 1): Kr12 = Kr(b12 / b02, 2): Kr13 = Kr(b12 / b02, 3)
Kr21 = Kr(b22 / b12, 1): Kr22 = Kr(b22 / b12, 2)
Kr23 = Kr(b22 / b12, 3): Kr31 = Kr(b42 / b32, 1) '调用自定义函数7求槽下部比漏磁导系数Kr
λh12 = (b02 * h12 ^ 3 * Kr11 + AB23 * h12 ^ 2 * Kr12 + AB23 ^ 2 * h12 * Kr13 / b12) / AB ^ 2
λh22 = (b12 * h22 ^ 3 * Kr21 + AB3 * h22 ^ 2 * Kr22 + AB3 ^ 2 * h22 * Kr23 / b22) / AB ^ 2
λh32 = (b32 * h32 ^ 3 * Kr31) / AB ^ 2: λL2 = λh12 + λh22 + λh32
End Sub
Private Sub cmdend_Click()
End
End Sub
Private Sub h11Lbl_Click()
b01 = Text26: b11 = Text28: α1 = Text31: Pi = 3.1415926
h11 = (b11 - b01) / 2 * Tan(α1 * Pi / 180): h11 = Format(h11, "#.##"): Text29.Text = Str(h11)
End Sub
Private Sub h12Lbl_Click()
b02 = Text33: b12 = Text35: α2 = Text39: Pi = 3.1415926
h12 = (b12 - b02) * Tan(α2 * Pi / 180): h12 = Format(h12, "#.##"): Text40.Text = Str(h12)
End Sub
Private Sub dLbl_Click()
d11 = Text22: Δ = Text23
d = d11 * 1 + Δ: Text24.Text = Str(d)
End Sub
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -