开关磁阻电机设计MATLAB程序.docx
- 文档编号:6367184
- 上传时间:2023-01-05
- 格式:DOCX
- 页数:17
- 大小:20.84KB
开关磁阻电机设计MATLAB程序.docx
《开关磁阻电机设计MATLAB程序.docx》由会员分享,可在线阅读,更多相关《开关磁阻电机设计MATLAB程序.docx(17页珍藏版)》请在冰豆网上搜索。
开关磁阻电机设计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;
globalNqNrNsNphDaDsggilFela
globalmiu0sigma
globalwvuyhqpmn
globall_delta
globalalphthetabelta
globalwvuhbfmmiu0gammaghhmsdeltasigma
globalt1t2t3m1m2n1b1b2S5S6S6pn_theta1
globalp1p2_1piep2_2piep2_3piep2p3p4p5p6p7sum_p
globallFgFL_2DL1
globalbpsbprhcshcr
globalSgStsStrScsScrhmrltsltrlglcslcr
globalBtsBtrBgBcsBcrBps
globalHtsHtrHgHcsHcr
globalfai_2Faisum_Fi_s
globaltheta_1
globalL1LaoLuURp
globalW_pieTemTem_pie
globalfai_1mfai_sfai_mkmomiga
globalIIm_starImi_mi_s
globala0a1
globalSwkskiJ_pieSa_pieSaScuJI
globalUNPNnNitaA_pieB_delta_piePemTemTem_pie
globaltheta_ontheta_offtheta_ctheta_utheta_hrtheta_0theta_n1theta_2theta_3theta_4theta_abeta_sbeta_rbeta_sfbeta_rf
globallambdaDa_ratio_Ds
globalDids
globaltal_rB_deltaIi_mIm_starSwSaScuksJlavlRpGcuVsFeVrFeGFeAPcuK
globalalphtal_r_degreetal_s_degree
globalwvuyhqpmnmiu0
globalP1P2P3P4P5sum_PgFsigmalF
globalLuL2D
globalLuikPasi_uPasi_1Pasi_a
globalfai_msfai_mi_msi_maisoPasi_so
globalNphNsDsDabpslFeggi
globalhmrhmshms_piehcshcr
globalBts1Bts2Bts3Bts4Btr1Btr2Btr3Btr4BgBcs1Bcs2Bcr1Bcr2B_minB_maxB_aminB_amax
globalSts1Sts2Sts3Sts4Str1Str2Str3Str4Scs1Scs2Scr1Scr2ScsScrSg
globalHts1Hts2Hts3Hts4Htr1Htr2Htr3Htr4Hcs1Hcs2Hcr1Hcr2Hg
globallts1lts2lts3lts4ltr1ltr2ltr3ltr4lcs1lcs2lcr1lcr2lg
globalFai_amfai_0fai_2amfai_ami_amSum_Fsum_Fami_shrfai_shrmiu0Lao_hrgi_piedelta_betakdkmsigma_m
globalAcBcCckaPasi_1utheta_1utheta_xPasi_thetaxPasi_a_mtrixPasiapasiaPasihrpasihrPasi_hr_mtrixPasi_u_mtrixPasi_1_mtrixMatrix
globalPasi_ahrtheta_ahrAc_pieBc_pieCc_pieTorquecurrentLao_hrPasi1k1Pasihrk1ik1ik2ik3
globalk1k2k3k4k5k6k7k8k9k10LengthNsteph
Nq=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;%额定功率Kw
nN=1500;%额定转速1500r/min
ita=0.88;%额定效率
Pem=PN*(1+ita)/(2*ita);%5.电磁功率Kw
lambda=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=(Da/2)*sin((tal_r_degree-beta_rf)/2)-(bps/2);%已阅
%%%q=v/2;%已阅
%%%p=bps/2;
%%%m=u+q;%已阅
%%%n=y+q;%已阅
%%%P1=alph*(m^2)/(4*(w^2)*((2*v+u)^2));
%%%P2=(2/pi)*(log(n/h)+2*(n-h)*y/(w*v)-((n^2-h^2)/(4*((w*v)^2)))*(pi*w*v-2*(y^2))-(n^3-h^3)*y*pi/(6*((w*v)^2))+(n^4-h^4)*(pi^2)/(64*((w*v)
^2)));
%%%P3=(2/tal_r_degree)*log((2*tan(tal_r_degree-beta_rf)+pi-(tal_r_degree-beta_rf))/(2*tan(tal_r_degree-beta_rf)+pi-2*(tal_r_degree-beta_rf)-
beta_rf));
%%%P4=(2/(pi-tal_r_degree))*log(2*gi/(h*(pi-tal_r_degree)));
%%%P5=(p+h)/gi-(2/(pi-tal_r_degree));
%%%sum_P=P1+P2+P3+P4+P5;
%%%gF=(1/5)*(pi*n/2-pi*h/2+((pi/2-(tal_r_degree-beta_rf)/2))*h+g+gi);
%%%sigma=(2/pi)*(atan(2*n/gF)-(gF/(4*n))*log(1+((2*n)/gF)^2));
%%%lF=la+2*n*(1-sigma);
%%%miu0=0.4*pi*(10^(-6));%H/m
%%%L2D=(Nph^2)*miu0*lFe*(10^(-3))*sum_P;%此处乘以10^(-3)将铁心lFe:
mm->m
%%%Lu=L2D*(2*(lF/la)-1);%23.theta_u不对齐位置电感
Theta_u_Lu();
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%-theta_1位置磁化曲线
%%%h=(Da/2)*sin(tal_r_degree-beta_rf-(beta_sf/2))-(Da/2+g)*sin(beta_sf/2);
%%%b=(Da/2+g)*cos(beta_sf/2)-(Da/2)*cos(tal_r_degree-beta_rf-beta_sf/2);
%%%%b=(Da/2-g)*cos(beta_sf/2)-(Da/2)*cos(tal_r_degree-beta_rf-beta_sf/2);
%%%f=sqrt(h^2+b^2);
%%%theta=atan(h/(gi+g-b));
%%%%beta1=pi/2-beta_sf/2;
%%%beta1=(pi/2)-(beta_sf/2)-(beta_rf/2);%平行转子齿
%%%%beta2=pi/2-(tal_r_degree-beta_rf-beta_sf/2);
%%%beta2=(pi/2)-(tal_r_degree-beta_rf-(beta_sf/2))-(beta_rf/2);
%%%gamma=theta-(tal_r_degree-beta_rf-beta_sf/2)-beta_rf/2;
%%%iff>=gi
%%%gh=(Da/2)*(cos(beta_sf/6)-cos(beta_sf/2))+g;
%%%else
%%%gh=(Da/2+g)*(cos(beta_sf/2-beta_sf/7)-cos(beta_sf/2))+g;
%%%end
%%%%w=(Da/2+g)*tan(pi/Ns)-bps/2;
%%%%alph=pi/2-(asin((w+bps/2)/(Ds/2-hcs))+asin((bps/2)/(Ds/2-hcs)))/2;
%%%%v=sqrt((Ds/2-hcs)^2-((Da/2+g)^2)*(tan(pi/Ns)^2))-(Da/2+g);
%%%hms=(Ds/2)-hcs-(Da/2+g);
%%%t1=(gi-gh)/beta1;
%%%t2=(h/sin(theta))-gi+g;
%%%t3=(gi-f)/beta2;
%%%m1=u+(v/2);
%%%m2=u+(2/3)*v;
%%%n1=hms-m1;
%%%%delta=pi/2+beta_sf/2-acos(n1/Da);
%%%delta=(pi/2)-acos(n1/Da)+(beta_sf/2)*alph;
%%%p1=alph*(m^2)/(4*(w^2)*((2*v+u)^2));
%%%p2_1pie=(1/(pi/2+delta))*log((gh+(pi/2+delta)*n1)/gh);
%%%p2_2pie=(pi/(2*w*v))*((n1^2)/(2*(pi/2+delta))-gh*n1/((pi/2+delta)^2)+(gh^2)/((pi/2+delta)^3)*log((gh+(pi/2+delta)*n1)/gh));
%%%p2_3pie=((pi^2)/(16*(w^2)*(v^2)))*(n1^4/(4*(pi/2+delta))-(gh/(pi/2+delta))*(n1^3/(3*(pi/2+delta))-gh*(n1^2)/(2*((pi/2+delta)^2))+n1*
(gh^2)/(pi/2+delta)^3-(gh^2/(pi/2+delta)^4)*log(((pi/2+delta)*n1+gh)/gh)));
%%%p2=p2_1pie-p2_2pie+p2_3pie;
%%%p3=(1/beta1)*log((beta1*t1+gh)/gh);
%%%iff>=gi
%%%b1=sqrt((gi+t2)^2-gi^2);
%%%S5=(w*v-t2*b1/4)/(w*v);
%%%b2=sqrt((gi+hms-m2)^2-(gi^2));
%%%S6=(w*v-v*b2/12-t2*b1/4)/(w*v);
%%%p4=(bps-t1)/gi;
%%%p5=(S5^2)/(theta-gamma)*log(theta*(gi+t2)/(theta*gi+gamma*t2));
%%%%p6=(S6^2)/theta*log(theta*(gi+hms-m2)/(theta*(gi+t2)));
%%%p6=(S6^2)*(1/theta)*log((gi+hms-m2)/(gi+t2));
%%%gF=(1/8)*(((pi/2)*n1+gh+delta*n1)+gh+(gh+beta1*t1)+2*gi+(gi*theta+gamma*t2)+(gi+t2)*theta+(gi+hms-m2)*theta);
%%%else
%%%S6p=1-v/(36*w*tan(pi/Ns));
%%%p4=(bps-t1-t3)/gi;
%%%p5=(1/beta2)*log((f+beta2*t3)/f);
%%%p6=(2/pi)*(S6p^2)*log((f+(pi/2)*(hms-m2))/f);
%%%gF=(1/8)*(((pi/2)*n1+gh+delta*n1)+gh+(gh+beta1*t1)+2*gi+(f+beta2*t3)+f+(f+(pi/2)*(hms-m2)));
%%%end
%%%p7=alph*(m2^4)/(4*(w^2)*((2*v+u)^2));
%%%sum_p=p1+p2+p3+p4+p5+p6+p7;
%%%n_theta1=n1+gh;
%%%%xi=(pi/2)*(atan(2*n/gF)-(gF/(4*n))*log(1+((2*n/gF)^2)));
%%%lF=la+2*n_theta1*(1-sigma);
%%%L_2D=2*miu0*(lFe*(10^(-3)))*((Nph/2)^2)*sum_p;
%%%L1=L_2D*(2*lF/la-1);
Theta_1_L1()
%%-Theta_a对齐位置时的磁化曲线
%%%Sg=(bps+2*g)*lFe;%气隙截面积/mm^2
%%%Sts=bps*lFe;%定子极截面积/mm^2
%%%Str=bpr*lFe;%转子极截面积/mm^2
%%%Scs=hcs*lFe;%定子轭截面积/mm^2
%%%Scr=hcr*lFe;%转子轭截面积/mm^2
%%%hmr=gi;%第二气隙
%%%lts=2*(Ds/2-hcs-Da/2-g);
%%%ltr=2*hmr;
%%%lg=2*g;
%%%lcs=(Ds-hcs)*pi/2;
%%%lcr=(Da/2-hmr-hcr/2)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 开关 磁阻 电机 设计 MATLAB 程序