msp10.m

基于MATLAB电机磁场程序

%1:entrefer, 2:culasse stator, 3: fer rotor,  
%4?7: conducteurs stator,  28 et 29 conducteurs rotor 


entref=1e-3; %entrefer
Rs=39.385e-3;	%rayon int閞ieur stator
Rc=71.75e-3; 	%rayon culasse (ext閞ieur stator)
Rr=Rs-entref;		%rayon rotor
Nb=140; %nbre spires par phase
Long=0.125; %longueur utile machine
p=1;
f=50;
beta=35*pi/180; beta=beta/p; %culasse rotor
beta1=50*pi/180; beta1=beta1/p; % 閜anouissement polaire
dessinP1


%aa=0;
gr=groto(gr,(aa)); 

g=[gsf gr];

limites

[p,e,t]=initmesh(g);
%[p,e,t]=refinemesh(g,p,e,t,[1 30]);

%pdegplot(g); axis equal
%hold on ; pdemesh(p,e,t)
%pause

% D閒inition des constantes
    kexc=0;  
    sigmaexc=kexc*5e6;	
    nuo=1/(4e-7*pi); murs=500; murr=500;
    ks=1; kr=0;
    Jex1=1e6*ks; Jex11=-Jex1;
    Jex2=-Jex1/2; Jex21=-Jex2;
	Jex3=-Jex1/2; Jex31=-Jex3;
    Jexr1=1e6*kr; Jexr2=-Jexr1;

    garnissageP1	

    % Assemblage des matrices
u=assempde(bond,p,e,t,nu,0,J);