1、maxwell 電機(jī)氣隙磁密與用matlab進(jìn)行fft諧波分析1對電機(jī)進(jìn)行靜態(tài)場分析，分析完后，進(jìn)入后處理下載 (37.19 KB)2010-8-4 09:222需要在氣隙中間畫一條圓弧線。點開deometry菜單，點creat再選Arc如下圖所示。下載 (24.48 KB)2010-8-4 09:22然后輸入圓弧的中心（0，0）回車。在下一個界面輸入起始點坐標(biāo)。最后一個界面輸入這條弧線上的采樣點數(shù)（250），圓弧角90度，圓弧的分段數(shù)目（250），名字以及線的顏色，最后回車，就會得到下圖的圓弧了。下載 (13 KB)2010-8-4 09:22    下載 (18.6 K

2、B)2010-8-4 09:22    下載 (18.29 KB)2010-8-4 09:22下載 (6.56 KB)2010-8-4 09:223. 需要得到氣隙磁密。打開后處理計算器，依次選擇qtyB，即選擇磁密矢量。選擇geomlineairgap_line, 即選中剛才畫的那條弧線。選擇unit vect2d normal,求取圓弧線的徑向分量。選擇dot（點乘），求取圓弧線上的B的徑向分量。再選一遍那個圓弧線，然后點 2d plot,就會出現(xiàn)那個磁密分布圖了。下載 (30.96 KB)2010-8-4 09:22下載 (31.94 KB)2010-8-4 09:2

3、2下載 (25.33 KB)2010-8-4 09:22下載 (29.77 KB)2010-8-4 09:22下載 (30.82 KB)2010-8-4 09:22下載 (35.08 KB)2010-8-4 09:22下載 (17.12 KB)2010-8-4 09:224. 雖然maxwell本身也可以做fft分析，但小弟還是喜歡把數(shù)據(jù)導(dǎo)出來在matlab 中進(jìn)行分析，這樣更靈活一些。導(dǎo)出數(shù)據(jù)。點擊plot菜單save as2d plot。在彈出的對話框中輸入數(shù)據(jù)文件的名字。（小弟實在找不到更好的辦法導(dǎo)出數(shù)據(jù)了，如有哪位達(dá)人有更好的方法，請賜教。小弟在此謝謝了。）下載 (11.8 KB)20

4、10-8-4 09:22下載 (11.38 KB)2010-8-4 09:225.  對氣隙磁密進(jìn)行諧波分析。將第四步中生成的.dat文件拷出來放到一個文件夾中（保證matlab和數(shù)據(jù)文件的路徑相同）。然后將matlab文件也拷貝到這個文件夾中。打開這個m文件，輸入Ns=500(需要進(jìn)行分析的采樣點個數(shù)，由于我們在maxwell中只分析了一個磁極下的磁密，所以只有半個周期，我們需要通過鏡像生成后半個周期，這樣總采樣點個數(shù)為250*2=500)。Order是需要分析的諧波次數(shù)，輸入11就是分析到11次諧波。運行，就得到下面的兩個圖了，第一個是原始波形，基波分量以及各高次諧波；

5、第二個是個諧波分量的幅值大小柱狀圖。2010-8-4 09:22下載 (18.78 KB)2010-8-4 09:22這樣一個電機(jī)氣隙磁密諧波分析就完成了。clcclear all;format long; Ns=500;order=11; %*read the position and flux density*fid=fopen('B.dat','r'); %open the original filefidnew = fopen('b1.dat','w'); %write the new filewhile feof(fid

6、)=0 tline = fgetl(fid); %tline? if ischar(tline), break, end temp=abs(tline); Nlength=length(tline); isemptyline=0; %? if Nlength=0 isemptyline=1; end allspace=0; %? isspace=0; for i=1:Nlength T=temp(i); if T=32 isspace=isspace+1; end if isspace=Nlength allspace=1; break end end findalpha=0; %? for

7、j=1:Nlength T=temp(j); if (T>=65)&(T>=90)|(T>=97)&(T>=122) findalpha=1; break; end end if (findalpha)&(allspace)&(isemptyline=0) %? fprintf(fidnew,tline);fprintf(fidnew,'n'); end end fclose(fid);fclose(fidnew); fid1=fopen('b1.dat','r');flux_positio

8、n =fscanf(fid1,'%f',2,Ns);fclose(fid1);%*read file finish*flux_position=flux_position'pos=flux_position(:,1); flux=flux_position(:,2); figure;plot(pos,flux,'r');%plot origional waveformhold on;grid on;fft1=fft(flux,Ns);j=0;amp_har=zeros(1,(order+1)/2); for m=1:2:order j=j+1; fft1=fft(flux,Ns); fund_ele_front=fft1(m+1); fund_ele_back=fft1(Ns+1-m); amp_har(j)=(abs(fund_ele_front)/Ns*2; fft1=0*fft1; fft1(m+1)=fund_ele_front; fft1(Ns+1-m)=fund_ele_back; fft1=ifft(