多回路方法流程

1、输入电机电磁设计参数 sector=7; %每分支绕组划分的段数 %sector=1;%适用于计算Ld,Lq,L0,Ll等参数 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;%空载或负载时的气隙放大倍数K_saturation=1.8960;%根据电机运行工况适当调整其数值dmax=air_gap_length_max*Kcarter*K_saturation;%气隙

2、有效长度，乘以卡氏系数与饱和系数dmin=air_gap_length_min*Kcarter*K_saturation;%气隙有效长度，乘以卡氏系数与饱和系数 z_stator=42;m_stator=6;%定子支路数N_stator=7;%定子一个支路所包含的线圈数与抽头数之和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;%参见兰州电机

3、厂12kW电机计算程序原始数据形式 tLab=27.5e-6;rs1=1.5e-1;xc_stator=1.0e-12;wfd=94;afd=1;lfdl=0.01346;%励磁绕组漏磁自感系数 rffd=0.6;rotor_core_length=0.13;% 转子铁心长lfdl_end=0.00205;%励磁绕组端部漏感nc=6;%每磁极下转子阻尼条数sl=0.415e-6;%阻尼条漏感el=0.00724e-6;% 端环漏感sr=0.595e-3;%阻尼条电阻er=0.00417e-3;% 端环电阻dal=-0.5803 -0.3482 0.0 0.3482 0.5803 1.5708;

4、%第 1磁极下各阻尼回路中心线领先于转子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磁极下各阻尼回路的端环长度系数N_stator0=7;%定子一个支路所包含的线圈数%转子磁极尺寸,mpole_bp=0.121;% 极靴宽度pole_sate=asin(pole_bp/(dsi-2*air_gap_length_max);%极 靴包角a1=2*tp*pole_sate/pi;% 极弧系数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;% 极身长pole_bm=0.065;% 极身宽pole_hm=0.0384;% 极身高pole_Djm=0.1097;%转子磁轭外径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;%极靴漏磁导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;%极身漏磁导%计算励磁绕组的端部漏感、漏磁自感系数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);%励磁绕组的比漏磁导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;%励磁绕组端部漏感系数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;%励磁绕组漏磁自感系数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;%根据实际转子的几何尺寸确定以下各个参数，气

10、隙最小长度未考虑卡氏系数，此处参照交流电机及其系统的分析附图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;%上、下层线棒分界线位置参数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);输入绕组连接信息(明确所选定的回路由那个或哪些线圈组成)计算发电机的气隙导磁系数附录B式B2,乘以卡氏系数与饱和系数式B5,标么值，转化为有名值%计算发电机的气隙导磁系数%lp(k)： k=115,极靴部分的气隙导磁系数%ldi(k)： k=115，当磁通势零值位于极间中心线时的极间部分的气隙导磁系数%lqi(k)： k=115，当磁通势最大值位于极间中心线时的极间部

13、分的气隙导磁系数%ld(k): k=115，当磁通势零值位于极间中心线时的气隙导磁系数，ld(k)=lp(k)+ldi(k)%lq(k): k=115,当磁通势最大值位于极间中心线时的气隙导磁系数，lq(k)=lp(k)+lqi(k)%ldq(k)： k=115,通常计算时使用的气隙导磁系数，ldq(k)=(ld(k)+lq(k)/2计算定子各个回路之间的互感系数a%计算定子绕组各个线圈之间的互感系数(与气隙主磁链有关)交流电机及其系统的分析式(1-3-38) (1-3-20)式(1-3-39),(1-3-21)b 计算定子绕组各个回路之间的互感系数(与气隙主磁链有关)根据绕组连接信息叠加c

14、计算定子绕组各个线圈之间的端漏互感附录Cd 计算定子绕组各个回路之间的端漏互感根据绕组连接信息叠加e定子绕组各个回路之间的槽漏互感%mab0(m_stator1,m_stator1)：定子各个回路之间的互感的常数项(包括漏感)%mab2(m_stator1,m_stator1)：定子各个回路之间的互感的二次谐波项的幅值%mabal(m_stator1,m_stator1)：定子各个回路之间的互感的二次谐波项的相角%定子第 i 号与第 j 号回路之间的互感为：mab0(i,j)+mab2(i,j)*cos(2*(wt+mabal(i,j)计算定子各个回路的电阻计算定子各个回路与励磁绕组之间的互感

15、系数交流电机及其系统的分析式(1-6-3)、根据绕组连接信息叠加计算定子各个回路与各个阻尼回路之间的互感交流电机及其系统的分析式(1-6-5)、根据绕组连接信息叠加计算励磁绕组与各个阻尼回路之间的互感交流电机及其系统的分析式(1-5-24)计算各个阻尼回路之间的互感交流电机及其系统的分析式(1-5-17)计算各个阻尼回路的电阻计算各个阻尼回路的漏感交流电机及其系统的分析1-8节计算励磁绕组自感系数交流电机及其系统的分析式(1-5-3)计算励磁绕组漏磁自感系数交流电机及其系统的分析1-8节以下采用数值算法迭代求解同步发电机正常、稳定运行状态，即数字仿真。发电机保持同步速ws( rad/s )旋转。设定各状态变量初值假设在0时刻同步发电机已由原动机拖动至同步速旋转，但尚未施加励磁，则除转子转速外的各状态变量 数值均为0。假设在0+时刻施加励磁电流，其数值为Ifd0。设定迭代步长tStep=0.001、迭代次数tIterration=2000、转子初始位置角gama