求助关于基于MATLAB的电机磁场建模与分析问题

thy562

请高手来解答一下，急需 。最好两天内，谢谢了。

: Y- p9 H2 j) U( n0 T
你问题太空洞了，别人怎么解答？？

thy562

请帮忙用有限元法做个同步电机电机磁场的模型相关程序（基于MATLAB）如下：% 1:air-gap, 2:culasse stator, 3: fer rotor,  
9 d+ F# X  w+ I! X( I% 4-7: stator conductor,  28 - 29: rotor conductor
5 q& d7 s- F. R% y. u/ n
3 F- J9 m/ }* _. [. `tic; clear;clc; close all;

5 I3 I: W) m. \  L  V. A4 _; P  q
$ W  N8 f4 u9 I3 x- r8 r
; o9 I4 s7 A( F, ~%     ------------- Draw the geometry of stator, rotor et air-gap --------------
9 x' p4 L* H& k5 T6 m$ Q7 L
, M  L1 ?* o2 r9 O) qentref=1e-3;    % air-gap depth
Rs=39.385e-3;        % radius of (rotor+airgap)
: d% K8 h& }6 ~( `$ n; W. B3 FRr=Rs-entref;        % radius of rotor
, m& v3 b, f/ z- K: g( ?) QRc=71.75e-3;         % radius of stator
; c8 E+ A6 t7 b4 `
Nb=140;         % number of turns per phase
9 \. y' f* X- Y" O& i8 ELong=0.125;     % longitude of machine in rotation axis direction
p=1;            % ??
f=50;           % ??

beta=35*pi/180; beta=beta/p;    % culasse rotor
beta1=50*pi/180; beta1=beta1/p; % epanouissement polaire

+ S+ K* D) V- v& U1 rdessinP1;                        % a sub-programme to draw the geometry of the machine section with the parameters defined above
%     --------------


%     ------------- Rotor positioning and initial mesh --------------

aa=-30*pi/180;         % initial angle of rotor
gr=groto(gr,(aa));     % rotating the rotor to the angle specified above using the sub-programme 'groto'
9 ~8 Y* A. @6 C0 Yg=[gsf gr];            % stator + rotor given the whole machine geometry
9 [/ s8 a% h4 z( b- M  w/ Mlimites;               % a sub-programme defining the bondary condition
[p,e,t]=initmesh(g);   % initial mesh for the partial differential equation resolution
4 u; |6 P: S! ^0 w- C( g4 O% g$ a
%   ------------- mesh refine if necessary ------------
; |( h, R9 e; u, v; r% [p,e,t]=refinemesh(g,p,e,t,[1 30]);
- O: I# U9 S- s& Q% pdegplot(g); axis equal, hold on ; pdemesh(p,e,t), hold off
% pause
%    -------------------------------------
4 t# `- W0 B5 m( m$ T

* @1 Y( o5 s) Z4 W, a( ]) A
, y4 j" E0 w6 j9 u%     -------------  specification of parameters required for calculation -------
    kexc=0;  
    sigmaexc=kexc*5e6;                       % conductivity of rotor coil conductor
: b" `0 t* c5 q; I) {    nuo=1/(4e-7*pi);                   % inverse of permeability of vacuum
. {7 k: l' f- _0 v  ]. A! ?    murs=500; murr=500;                % relative permeability of stator and rotor
  X3 Y* f9 n* `9 ?) B1 q- M" @                                    
    ks=1; kr=0;                        % ks=1 or 0: stator excited or not; kr=1 or 0: rotor excited or not.
   
1 Y  C+ S# r  B7 x0 }    Jex1=0;    Jex11=0;
    Jex2 = 10e6*ks;     Jex21 = -Jex2;
1 X9 [4 c! O6 c; d    Jex3 = 0;     Jex31 = 0;
%     Jex3=10e6*ks; Jex31=-Jex3;         % current density of phase a (stator)
9 @. g( `0 q* c4 c- G%     Jex1=-Jex3/2; Jex11=-Jex1;         % current density of phase b (stator)
% O- I. A- C, q# ]/ q7 \; `2 Y%         Jex2=-Jex3/2; Jex21=-Jex2;         % current density of phase c (stator)
    Jexr1=10e6*kr; Jexr2=-Jexr1;       % current density of rotor excitation

1 i( N7 n) }9 Z+ }5 `# G4 l    garnissageP1;                           % a sub-programme for the giving the current density and permeability at each mesh element
. w$ L) f! |' |) z9 G" K3 @: h3 R%     -------------
2 M# x8 e# K! U: V
9 f" z( z2 l: c! v, ]* y
- \3 z& b* r/ {) _5 N  A
( Y1 E( S2 Q6 c( N  |%     -------------   resolution of the partial differential equation and figure ploting ---------

2 Q/ \! k/ B# }) y, @u=assempde(bond,p,e,t,nu,0,J);         % resolution of the partial differential equation
% U=ASSEMPDE(B,P,E,T,C,A,F) assembles and solves the PDE problem -div(c*grad(u))+a*u=f

figure (1), pdemesh(p,e,t), axis equal, title ('Mesh of the machine section')
. d* s5 L* V: d7 M% ]$ Z7 Qfigure (2), pdegplot(g); axis equal; hold on; pdecont(p,t,full(real(u)),30), hold off, title('Magnetic field distribution'), % plots using 30 levels.
%    ------------------------------
. ]1 F+ f  }4 S" n


' {6 O' g4 I5 ^% -------------- calculation of the energy stocked in the machine (energy sum in the magnetic field) --------------------

. m) T+ d# E- Z/ e% z- d. n5 a% [ux,uy]=pdegrad(p,t,u);
: g0 \- H  `% M* j. E5 F% bx=uy;
% by=-ux;
) X) ?3 o% C: h# M$ }7 {; e% for i=1:length(bx);
1 [8 g( B; v  l4 ~3 _" {5 M%     b(i)=(bx(i)^2+by(i)^2)^.5;
% end
, D% a3 K& o1 S* ]% nu_e=1/(4e-7*pi);
; j6 C* |3 h) G0 {- Z% nu_c=nu_e;
7 Z+ x$ U2 X8 S) W! d6 r% nu_s=1/500*nu_e;
% nu_r=1/500*nu_e;
% h=sparse(1,ntrg);
% ind_e=find(t(4,:)==2);
1 Z0 \4 [. M" {2 u! q! m. X% ind_s=find(t(4,:)==1);
+ g% ~8 j' e' R% k, \1 I' c% ind_r=find(t(4,:)==3);
5 U4 @) S! C5 e, U6 t% ind_c=find((t(4,:)>=4) & (t(4,:)<=27));
) y& @& [3 \# h  ~- d0 v% h(ind_e)=nu_e*b(ind_e)';
8 V) I% w4 |7 Z0 \5 _2 M! ^2 c3 I% h(ind_s)=nu_s*b(ind_s);
% h(ind_r)=nu_r*b(ind_r);
7 P6 [  q! F7 j1 D  V% h(ind_c)=nu_c*b(ind_c)';
& p+ I& j* @1 u" ~* A) P4 Q% aire=pdetrg(p,t);
% vol_trg=Long*aire;
% energ_elem_e=1/2*(h(ind_e).*b(ind_e)').*vol_trg(ind_e);
% energ_tot_e=sum(energ_elem_e);
0 Q0 }* L- _1 m) C2 t% energ_elem_s=1/2*(h(ind_s).*b(ind_s)').*vol_trg(ind_s);
% energ_tot_s=sum(energ_elem_s);
  O4 r& p1 g: A2 E3 f  v5 g3 i% energ_elem_r=1/2*(h(ind_r).*b(ind_r)').*vol_trg(ind_r);
4 d+ n, l9 e& Q8 Y* ?% energ_tot_r=sum(energ_elem_r);
% ]! W; |# T9 \8 f% o% energ_elem_c=1/2*(h(ind_c).*b(ind_c)').*vol_trg(ind_c);
9 k! p! G: E. ?& p3 W5 S% energ_tot_c=sum(energ_elem_c);
% energ_total=energ_tot_e+energ_tot_s+energ_tot_r+energ_tot_c;
% enr=energ_total;
7 ?( O5 M: ^0 Y8 f& V7 d4 F0 @
toc;

zhoushaoxing

太深奥了，都没有看明白。。。