energie1.m

基于MATLAB电机磁场程序

%calcul de l'energie
[ux,uy]=pdegrad(p,t,u);
bx=uy;
by=-ux;
for i=1:length(bx)
    b(i)=(bx(i)^2+by(i)^2)^.5;
end
nu_e=1/(4e-7*pi);%inverse de la permeabilite de l'air
nu_c=nu_e;
nu_s=nu_e/500*ones(1,ntrg);%inverse de la permeabilite du fer
nu_r=nu_e/500*ones(1,ntrg);
h=sparse(1,ntrg);
ind_e=find(t(4,:)==2);
ind_s=find(t(4,:)==1);
ind_r=find(t(4,:)==3);
ind_c=find((t(4,:)>=4)&(t(4,:)<=27));
bstr=b(ind_s)'%matrice d'induction stator transpos閑(ligne)
brtr=b(ind_r)'%matrice d'induction rotor transpos閑(ligne)
betr=b(ind_e)'%matrice d'induction entrefer transpos閑(ligne)
bctr=b(ind_c)'%matrice d'induction conducteur transposee(ligne)
h(ind_e)=nu_e*b(ind_e);
h(ind_s)=nu_s(ind_s).*bstr;
h(ind_r)=nu_r(ind_r).*brtr;
h(ind_c)=nu_c*b(ind_c);
aire=pdetrg(p,t);
vol_trg=Long*aire;
energ_elem_e=0.5*h(ind_e).*betr.*vol_trg(ind_e)
energ_tot_e=sum(energ_elem_e);
energ_elem_s=0.5*h(ind_s).*bstr.*vol_trg(ind_s)
energ_tot_s=sum(energ_elem_s);
energ_elem_r=0.5*h(ind_r).*brtr.*vol_trg(ind_r)
energ_tot_r=sum(energ_elem_r);
energ_elem_c=0.5*h(ind_c).*bctr.*vol_trg(ind_c)
energ_tot_c=sum(energ_elem_c);
energ_total=energ_tot_e+energ_tot_s+energ_tot_r+energ_tot_c
enr=energ_total;