开关磁阻电机设计MATLAB程序

%%%%-开关磁阻电机设计实例%%%%-Time:2013.08.20%%%%-BySunPeng(悟空)%%%------------------------------------------------------%%%-给定参数要求：%%功率：PN=11kw电源电压：UN=380V(线电压)额定转速：n=1500r/min额定效率：ita=0.88%%调速范围：200~2000r/min其中200~1500r/min为恒转矩特性1500~2000r/min为恒功率特性%%运行方式：连续运行绝缘等级：E电机转子承受的扭矩不得小于3500N.m%%%------------------------------------------------------%%%-相数、极数、绕组端电压:采用4相8/6极结构，功率变换器主电路选用每相一个主开关管的裂相式电路clc;clear;globalNqNrNsNphDaDsggilFelaglobalmiu0sigmaglobalwvuyhqpmngloball_deltaglobalalphthetabeltaglobalwvuhbfmmiu0gammaghhmsdeltasigmaglobalt1t2t3m1m2n1b1b2S5S6S6pn_theta1globalp1p2_1piep2_2piep2_3piep2p3p4p5p6p7sum_pgloballFgFL_2DL1globalbpsbprhcshcrglobalSgStsStrScsScrhmrltsltrlglcslcrglobalBtsBtrBgBcsBcrBpsglobalHtsHtrHgHcsHcrglobalfai_2Faisum_Fi_sglobaltheta_1globalL1LaoLuURpglobalW_pieTemTem_pieglobalfai_1mfai_sfai_mkmomigaglobalIIm_starImi_mi_sglobala0a1globalSwkskiJ_pieSa_pieSaScuJIglobalUNPNnNitaA_pieB_delta_piePemTemTem_pieglobaltheta_ontheta_offtheta_ctheta_utheta_hrtheta_0theta_n1theta_2theta_3theta_4theta_abeta_sbeta_rbeta_sfbeta_rfgloballambdaDa_ratio_DsglobalDidsglobaltal_rB_deltaIi_mIm_starSwSaScuksJlavlRpGcuVsFeVrFeGFeAPcuKglobalalphtal_r_degreetal_s_degreeglobalwvuyhqpmnmiu0globalP1P2P3P4P5sum_PgFsigmalFglobalLuL2DglobalLuikPasi_uPasi_1Pasi_aglobalfai_msfai_mi_msi_maisoPasi_soglobalNphNsDsDabpslFeggiglobalhmrhmshms_piehcshcrglobalBts1Bts2Bts3Bts4Btr1Btr2Btr3Btr4BgBcs1Bcs2Bcr1Bcr2B_minB_maxB_aminB_amaxglobalSts1Sts2Sts3Sts4Str1Str2Str3Str4Scs1Scs2Scr1Scr2ScsScrSgglobalHts1Hts2Hts3Hts4Htr1Htr2Htr3Htr4Hcs1Hcs2Hcr1Hcr2Hggloballts1lts2lts3lts4ltr1ltr2ltr3ltr4lcs1lcs2lcr1lcr2lgglobalFai_amfai_0fai_2amfai_ami_amSum_Fsum_Fami_shrfai_shrmiu0Lao_hrgi_piedelta_betakdkmsigma_mglobalAcBcCckaPasi_1utheta_1utheta_xPasi_thetaxPasi_a_mtrixPasiapasiaPasihrpasihrPasi_hr_mtrixPasi_u_mtrixPasi_1_mtrixMatrixglobalPasi_ahrtheta_ahrAc_pieBc_pieCc_pieTorquecurrentLao_hrPasi1k1Pasihrk1ik1ik2ik3globalk1k2k3k4k5k6k7k8k9k10LengthNstephNq=4;%1.相数Ns=8;%2.定子极数Nr=6;%3.转子极数UN=380;%额定电压Ud=3*sqrt(2)*UN/pi;%三相全波整流后直流电压U=1.1*(Ud/2);%4.绕组端电压（1.1系数为考虑电容滤波器存在导致的电压升高系数）%%%-主要尺寸选择PN=7.5;%额定功率KwnN=1500;%额定转速1500r/minita=0.88;%额定效率Pem=PN*(1+ita)/(2*ita);%5.电磁功率Kwlambda=1.2;%6.细长比（参照中小型交流电机经验数据，取0.5~3.0，对于SRM，典型取值1）A_pie=28000;%7.电磁负荷出选A_pie=28000A/m(15000~50000A/m);B_delta_pie=0.40T(0.29~0.55)B_delta_pie=0.4;ki=0.5;km=0.8;Da=(((6.1/(B_delta_pie*A_pie))*(ki/km)*(Pem*(10^3)/nN)*(1/(1.05*lambda)))^(1/3))*(10^3);%8.转子外径/mm/Da=Keep_decimal(Da,0);%la=lambda*Da;%9.贴心叠长la=Round_off_la(la);%%%-圆整铁心叠长l_delta=1.05*la;Da_ratio_Ds=Ratio_Da_Ds(Nq,Ns,Nr);%确定定转子外径比Ds=Da/Da_ratio_Ds;%10.定子外径Ds=Standardizing_Ds(Ds);%将定子外径规格化到已有机座New_Da_ratio_Ds=Da/Ds;%%%-其他结构尺寸及绕组匝数g=0.4;%11.气隙%beta_s=21;%12.定转子极弧，查表取经典值%beta_r=23;beta_s=Pore_arc_s(Nq,Ns,Nr);beta_r=Pore_arc_r(Nq,Ns,Nr);beta_sf=beta_s*pi/180;%定子极弧beta_rf=beta_r*pi/180;%转子极弧bps=(Da+2*g)*sin(beta_sf/2);%13.定转子极宽bps=Keep_decimal(bps,2);%保留两位小数bpr=Da*sin(beta_rf/2);bpr=Keep_decimal(bpr,2);%保留两位小数gi=bps/2;%14.第二气隙hcs=1.3*(bps/2);%15.定转子轭高hcs=Keep_decimal(hcs,2);%保留两位小数hcs_arter_roundoff=Round_off_h(hcs);%进行圆整hcr=1.4*(bpr/2);hcr=Keep_decimal(hcr,2);%保留两位小数hcr_after_roundoff=Round_off_h(hcr);%进行圆整Di=Da-2*(gi-g)-2*hcr;%16.轴径Di=Keep_decimal(Di,1);%保留一位小数ds=(Ds-Da-2*g-2*hcs)/2;%17.定子槽深ds=Keep_decimal(ds,2);%保留两位小数ds_after_roundoff=Round_off_h(ds);%进行圆整%[Di_aftergi_after]=Adapt_Di_gi(Di,Da,hcr,gi,g);%18.电机尺寸需根据工程设计特点进行圆整，取平行转子极结构[Digi]=Adapt_Di_gi(Di,Da,hcr,gi,g);kFe=0.93;lFe=kFe*la;%19.有效铁心长度lFe=Round_off_l(lFe);%进行圆整tal_r=(pi*Da)/Nr;%20.转子极距tal_r=Keep_decimal(tal_r,2);%保留两位小数tal_r_degree=2*pi/Nr;%%%%%%%%%%%%%%%%%%%%%%-以上程序段已经检查无误-%%%%%%%%%%%%%%%%%%%%%%%%%%%-控制参数theta_u=0;theta_on=0;%21.开通角theta_off=theta_u+((2*pi/Nr)-beta_rf)/2;%关断角theta_a=pi/Nr;theta_2=theta_a-abs(beta_rf-beta_sf)/2;theta_1=theta_a-(beta_sf+beta_rf)/2;%theta_0=-theta_1;%theta_0=theta_1-fai_s*theta_a/(12*fai_ms);%theta_0=theta_1-fai_s*theta_a/(12*fai_ms);theta_n1=(-1)*theta_1;theta_3=theta_a+(theta_a-theta_2);theta_4=tal_r_degree-theta_1;%theta_0=-10*pi/180;%theta_1=10*pi/180;theta_hr=theta_off;theta_c=theta_off-theta_on;%导通角%%%-每相绕组串联匝数Bps=1.6;B_delta=0.805*(bps/tal_r)*Bps;%重新校核Nph=(3.04*Nr*U*theta_c)/(nN*B_delta*(Da/(10^3))*(l_delta/(10^3)));%22.每相绕组串联匝数Nph1=(3.04*Nr*U*theta_c)/(nN*B_delta*(Da/(10^3))*(l_delta/(10^3)));Adapt_B_delta();%%%-对齐位置和不对齐位置磁化曲线计算%%-theta_u不对齐位置电感Lu计算%globalNrNsNphDaDsbpsggilFelahcsbeta_sfbeta_rfalphtal_r_degreetal_s_degree%globalwvuyhqpmnpi%%%tal_r_degree=2*pi/Nr;%转子极距角%%%tal_s_degree=2*pi/Ns;%定子极距角%%%%beta_sf=beta_s*pi/180;%定子极弧%%%%beta_rf=beta_r*pi/180;%转子极弧%%%%p=(Da/2+g)*sin(beta_sf/2);%%%%w=(Da/2+g)*tan(tal_s_degree/2)-p;%%%w=(Da/2+g)*tan(pi/Ns)-(bps/2);%%已阅%%%v=(Ds/2-hcs-(Da/2+g)/cos(pi/Ns))/cos(pi/Ns);%已阅%%%alph=(pi/2)-(beta_sf/2);%已阅%%%u=w/tan(alph);%已阅%%%%h1=sqrt((Ds/2-hcs)*(Ds/2-hcs)-p^2);%%%%h2=sqrt((Ds/2-hcs)*(Ds/2-hcs)-(p+w)^2);%%%%u=h2-h1;%%%%alph=atan(w/u);%%%%v=h1-(Da/2+g);%%%%y=(Da/2)+g-(Da/2)*cos(tal_r_degree-beta_rf)/2;%%%y=(Da/2)+g-(Da/2)*cos((tal_r_degree-beta_rf)/2);%已阅%%%%h=Da/2*sin((tal_r_degree-beta_rf)/2)-p;%%%h=