多回路方法流程

1、輸入電機(jī)電磁設(shè)計(jì)參數(shù) sector=7; %每分支繞組劃分的段數(shù) %sector=1;%適用于計(jì)算Ld,Lq,L0,Ll等參數(shù) tp=2;1=0.12;dsi=0.23;dso=0.327;air_gap_1ength=6e-4;air_gap_1ength_max=6e-4;air_gap_length_min=6e-4;KCarter=1.172;a1=0.7096;K_saturation=1;%空載或負(fù)載時(shí)的氣隙放大倍數(shù)K_saturation=1.8960;%根據(jù)電機(jī)運(yùn)行工況適當(dāng)調(diào)整其數(shù)值dmax=air_gap_length_max*Kcarter*K_saturation;%氣隙

2、有效長(zhǎng)度，乘以卡氏系數(shù)與飽和系數(shù)dmin=air_gap_length_min*Kcarter*K_saturation;%氣隙有效長(zhǎng)度，乘以卡氏系數(shù)與飽和系數(shù) z_stator=42;m_stator=6;%定子支路數(shù)N_stator=7;%定子一個(gè)支路所包含的線圈數(shù)與抽頭數(shù)之和ytol_stator=9;wks=14;as_stator=2;M1=2;N1=4;N2=9;K1=10;Radius_bar=0.124;Radius_gap=0.1147;Re=2.525e-3;B_mendz=0.02;C_mendz=0.063;tLa=24.1e-6;tLb=45.4e-6;%參見(jiàn)蘭州電機(jī)

3、廠12kW電機(jī)計(jì)算程序原始數(shù)據(jù)形式 tLab=27.5e-6;rs1=1.5e-1;xc_stator=1.0e-12;wfd=94;afd=1;lfdl=0.01346;%勵(lì)磁繞組漏磁自感系數(shù) rffd=0.6;rotor_core_length=0.13;% 轉(zhuǎn)子鐵心長(zhǎng)lfdl_end=0.00205;%勵(lì)磁繞組端部漏感nc=6;%每磁極下轉(zhuǎn)子阻尼條數(shù)sl=0.415e-6;%阻尼條漏感el=0.00724e-6;% 端環(huán)漏感sr=0.595e-3;%阻尼條電阻er=0.00417e-3;% 端環(huán)電阻dal=-0.5803 -0.3482 0.0 0.3482 0.5803 1.5708;

4、%第 1磁極下各阻尼回路中心線領(lǐng)先于轉(zhuǎn)子d軸的電角度 bata_damp=0.07389 0.07389 0.1478 0.07389 0.07389 0.5567;%第 1 磁極下各阻尼回路的短距比ee=1 1 2 1 1 7.534;%第1磁極下各阻尼回路的端環(huán)長(zhǎng)度系數(shù)N_stator0=7;%定子一個(gè)支路所包含的線圈數(shù)%轉(zhuǎn)子磁極尺寸,mpole_bp=0.121;% 極靴寬度pole_sate=asin(pole_bp/(dsi-2*air_gap_length_max);%極 靴包角a1=2*tp*pole_sate/pi;% 極弧系數(shù)hp=0.5*(air_gap_length_ma

5、x-air_gap_length)*(dsi-air_gap_length)/(dsi/2-air_gap_length)-(dsi-air_gap_length_max)*cos(pole_sate);%極靴偏心距pole_Lm=0.13;% 極身長(zhǎng)pole_bm=0.065;% 極身寬pole_hm=0.0384;% 極身高pole_Djm=0.1097;%轉(zhuǎn)子磁軛外徑pole_a=(dsi-pole_Djm)/2-pole_hm-air_gap_length_min;%極 靴厚度pole_b=sqrt(dsi-2*air_gap_length_max)A2-pole_bpA2)/2-po

6、le_Djm/2-pole_hm;pole_hpm=(2*pole_a+pole_b)/3;pole_batam=asin(pole_bm/pole_Djm);pole_hm1=pole_hm+pole_Djm*(1-cos(pole_batam)/2;pole_gama=atan(pole_bp/(pole_Djm+2*pole_hm);pole_gama1=atan(pole_bm/(pole_Djm+pole_hm);pole_alfap=sqrt(pole_bpA2+(pole_Djm+2*pole_hm)A2)*sin(pi/2/tp-pole_gama);pole_alfakm=s

7、qrt(pole_bmA2+(pole_Djm+pole_hm)A2)*sin(pi/2/tp-pole_gama1);pole_Ap=5.1*pole_Lm*pole_hpm/pole_alfap+7*pole_hpm;%極靴漏磁導(dǎo)pole_d1=4e-3;%磁極壓板厚度pole_b1=6e-3;%磁極壓板寬度pole_Lm1=pole_Lm+2*pole_d1;pole_Am=2.6*pole_Lm1*pole_hm1/pole_alfakm+2.1*pole_hm1;%極身漏磁導(dǎo)%計(jì)算勵(lì)磁繞組的端部漏感、漏磁自感系數(shù)pole_dt=pole_a+air_gap_length_min-po

8、le_bpA2/(4*dsi);pole_cp=pi*dsi/(2*tp)-pole_bp-2*pi*pole_dt/(2*tp);pole_cm=pi*dsi/(2*tp)-pole_bm-pi/(2*tp)*(pole_hm+2*pole_a+2*air_gap_length_min);%勵(lì)磁繞組的比漏磁導(dǎo)lamda_f_1=4*(pole_dt/pole_cp-0.25)+1.75*(pole_bp/2-pole_bm/2)/pole_cp+0.2)-1.27*(pole_bp/2-pole_bm/2)/pole_cp-0.5)A2+1.15*pole_hm/pole_cm;lamda_

9、f_2=0.44*pole_bm/rotor_core_length;%勵(lì)磁繞組端部漏感系數(shù)lfdl_end=2*tp*4e-7*pi*wfdA2/afdA2*lamda_f_1*(rotor_core_length-l)+2*tp*4e-7*pi*wfdA2/afdA2*lamda_f_2*rotor_cor e_length;%勵(lì)磁繞組漏磁自感系數(shù)lfdl=2*tp*4e-7*pi*wfdA2/afdA2*(lamda_f_1+lamda_f_2)*rotor_core_length;rp=dsi/2-air_gap_length_min-hp;%根據(jù)實(shí)際轉(zhuǎn)子的幾何尺寸確定以下各個(gè)參數(shù)，氣

10、隙最小長(zhǎng)度未考慮卡氏系數(shù)，此處參照交流電機(jī)及其系統(tǒng)的分析附圖2 1sita0=a1*pi/2;%定子槽形尺寸(m)stator_b0=3e-3;stator_b1=8.8e-3;stator_h1=1e-3;stator_h2=1.5e-3;stator_R=5e-3;stator_hn=15.5e-3;stator_h3=8e-3;stator_t2=0.0133;stator_bs2=1.5e-3;stator_intersection=6.3e-3;%上、下層線棒分界線位置參數(shù)stator_tz=0.0172;stator_tz1=0.0176;stator_tz2=0.0188;sta

11、tor_bz1=0.0088;stator_bz2=0.0088;KCarter1=stator_tz*(4.44*air_gap_length+0.75*stator_b0);KCarter1=KCarter1/(stator_tz*(4.44*air_gap_length+0.75*stator_b0)-stator_b0A2);KCarter2=stator_t2*(4.44*air_gap_length+0.75*stator_bs2);KCarter2=KCarter2/(stator_t2*(4.44*air_gap_length+0.75*stator_bs2)-stator_b

12、s2A2);KCarter=KCarter1*KCarter2;tol=0.5*pi*dsi/tp;gama_stator=2*tp*pi/z_stator;bata_s=ytol_stator/(z_stator/(2*tp);輸入繞組連接信息(明確所選定的回路由那個(gè)或哪些線圈組成)計(jì)算發(fā)電機(jī)的氣隙導(dǎo)磁系數(shù)附錄B式B2,乘以卡氏系數(shù)與飽和系數(shù)式B5,標(biāo)么值，轉(zhuǎn)化為有名值%計(jì)算發(fā)電機(jī)的氣隙導(dǎo)磁系數(shù)%lp(k)： k=115,極靴部分的氣隙導(dǎo)磁系數(shù)%ldi(k)： k=115，當(dāng)磁通勢(shì)零值位于極間中心線時(shí)的極間部分的氣隙導(dǎo)磁系數(shù)%lqi(k)： k=115，當(dāng)磁通勢(shì)最大值位于極間中心線時(shí)的極間部

13、分的氣隙導(dǎo)磁系數(shù)%ld(k): k=115，當(dāng)磁通勢(shì)零值位于極間中心線時(shí)的氣隙導(dǎo)磁系數(shù)，ld(k)=lp(k)+ldi(k)%lq(k): k=115,當(dāng)磁通勢(shì)最大值位于極間中心線時(shí)的氣隙導(dǎo)磁系數(shù)，lq(k)=lp(k)+lqi(k)%ldq(k)： k=115,通常計(jì)算時(shí)使用的氣隙導(dǎo)磁系數(shù)，ldq(k)=(ld(k)+lq(k)/2計(jì)算定子各個(gè)回路之間的互感系數(shù)a%計(jì)算定子繞組各個(gè)線圈之間的互感系數(shù)(與氣隙主磁鏈有關(guān))交流電機(jī)及其系統(tǒng)的分析式(1-3-38) (1-3-20)式(1-3-39),(1-3-21)b 計(jì)算定子繞組各個(gè)回路之間的互感系數(shù)(與氣隙主磁鏈有關(guān))根據(jù)繞組連接信息疊加c

14、計(jì)算定子繞組各個(gè)線圈之間的端漏互感附錄Cd 計(jì)算定子繞組各個(gè)回路之間的端漏互感根據(jù)繞組連接信息疊加e定子繞組各個(gè)回路之間的槽漏互感%mab0(m_stator1,m_stator1)：定子各個(gè)回路之間的互感的常數(shù)項(xiàng)(包括漏感)%mab2(m_stator1,m_stator1)：定子各個(gè)回路之間的互感的二次諧波項(xiàng)的幅值%mabal(m_stator1,m_stator1)：定子各個(gè)回路之間的互感的二次諧波項(xiàng)的相角%定子第 i 號(hào)與第 j 號(hào)回路之間的互感為：mab0(i,j)+mab2(i,j)*cos(2*(wt+mabal(i,j)計(jì)算定子各個(gè)回路的電阻計(jì)算定子各個(gè)回路與勵(lì)磁繞組之間的互感

15、系數(shù)交流電機(jī)及其系統(tǒng)的分析式(1-6-3)、根據(jù)繞組連接信息疊加計(jì)算定子各個(gè)回路與各個(gè)阻尼回路之間的互感交流電機(jī)及其系統(tǒng)的分析式(1-6-5)、根據(jù)繞組連接信息疊加計(jì)算勵(lì)磁繞組與各個(gè)阻尼回路之間的互感交流電機(jī)及其系統(tǒng)的分析式(1-5-24)計(jì)算各個(gè)阻尼回路之間的互感交流電機(jī)及其系統(tǒng)的分析式(1-5-17)計(jì)算各個(gè)阻尼回路的電阻計(jì)算各個(gè)阻尼回路的漏感交流電機(jī)及其系統(tǒng)的分析1-8節(jié)計(jì)算勵(lì)磁繞組自感系數(shù)交流電機(jī)及其系統(tǒng)的分析式(1-5-3)計(jì)算勵(lì)磁繞組漏磁自感系數(shù)交流電機(jī)及其系統(tǒng)的分析1-8節(jié)以下采用數(shù)值算法迭代求解同步發(fā)電機(jī)正常、穩(wěn)定運(yùn)行狀態(tài)，即數(shù)字仿真。發(fā)電機(jī)保持同步速ws( rad/s )旋轉(zhuǎn)。設(shè)定各狀態(tài)變量初值假設(shè)在0時(shí)刻同步發(fā)電機(jī)已由原動(dòng)機(jī)拖動(dòng)至同步速旋轉(zhuǎn)，但尚未施加勵(lì)磁，則除轉(zhuǎn)子轉(zhuǎn)速外的各狀態(tài)變量 數(shù)值均為0。假設(shè)在0+時(shí)刻施加勵(lì)磁電流，其數(shù)值為Ifd0。設(shè)定迭代步長(zhǎng)tStep=0.001、迭代次數(shù)tIterration=2000、轉(zhuǎn)子初始位置角gama