maxwell电机气隙磁密与用matlab进行fft谐波分析文档格式.docx

maxwell电机气隙磁密与用matlab进行fft谐波分析文档格式.docx

《maxwell电机气隙磁密与用matlab进行fft谐波分析文档格式.docx》由会员分享，可在线阅读，更多相关《maxwell电机气隙磁密与用matlab进行fft谐波分析文档格式.docx（14页珍藏版）》请在冰豆网上搜索。

下载（18.6KB）

下载（18.29KB）

下载（6.56KB）

3.需要得到气隙磁密。

打开后处理计算器，依次选择qty—B，即选择磁密矢量。

选择geom—line—airgap_line,即选中刚才画的那条弧线。

选择unitvect—2dnormal,求取圆弧线的径向分量。

选择dot（点乘），求取圆弧线上的B的径向分量。

再选一遍那个圆弧线，然后点2dplot,就会出现那个磁密分布图了。

下载（30.96KB）

下载（31.94KB）

下载（25.33KB）

下载（29.77KB）

下载（30.82KB）

下载（35.08KB）

下载（17.12KB）

4.虽然maxwell本身也可以做fft分析，但小弟还是喜欢把数据导出来在matlab中进行分析，这样更灵活一些。

导出数据。

点击plot菜单—saveas—2dplot。

在弹出的对话框中输入数据文件的名字。

（小弟实在找不到更好的办法导出数据了，如有哪位达人有更好的方法，请赐教。

小弟在此谢谢了。

）

下载（11.8KB）

下载（11.38KB）

5. 

对气隙磁密进行谐波分析。

将第四步中生成的.dat文件拷出来放到一个文件夹中（保证matlab和数据文件的路径相同）。

然后将matlab文件也拷贝到这个文件夹中。

打开这个m文件，输入Ns=500（需要进行分析的采样点个数，由于我们在maxwell中只分析了一个磁极下的磁密，所以只有半个周期，我们需要通过镜像生成后半个周期，这样总采样点个数为250*2=500）。

Order是需要分析的谐波次数，输入11就是分析到11次谐波。

运行，就得到下面的两个图了，第一个是原始波形，基波分量以及各高次谐波；

第二个是个谐波分量的幅值大小柱状图。

下载（18.78KB）

这样一个电机气隙磁密谐波分析就完成了。

clc

clearall;

formatlong;

Ns=500;

order=11;

%**********************readthepositionandfluxdensity************************

fid=fopen（'

B.dat'

'

r'

）;

%opentheoriginalfile

fidnew=fopen（'

b1.dat'

w'

%writethenewfile

whilefeof（fid）==0

tline=fgetl（fid）;

%tline?

if~ischar（tline）,break,end

temp=abs（tline）;

Nlength=length（tline）;

isemptyline=0;

%

ifNlength==0

isemptyline=1;

end

allspace=0;

isspace=0;

fori=1:

Nlength

T=temp（i）;

ifT==32

isspace=isspace+1;

ifisspace==Nlength

allspace=1;

break

findalpha=0;

forj=1:

T=temp（j）;

if（（T>

=65）&

（T>

=90））|（（T>

=97）&

=122））

findalpha=1;

break;

if（~findalpha）&

（~allspace）&

（isemptyline==0）%

fprintf（fidnew,tline）;

fprintf（fidnew,'

\n'

fclose（fid）;

fclose（fidnew）;

fid1=fopen（'

flux_position=fscanf（fid1,'

%f'

[2,Ns]）;

fclose（fid1）;

%********************************readfilefinish*****************************************

flux_position=flux_position'

;

pos=flux_position（:

1）;

flux=flux_position（:

2）;

figure;

plot（pos,flux,'

%plotorigionalwaveform

holdon;

gridon;

fft1=fft（flux,Ns）;

j=0;

amp_har=zeros（1,（order+1）/2）;

form=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（fft1,Ns）;

fft1=real（fft1）;

plot（pos,fft1）;

holdon;

k=（1:

order）;

bar（k,amp_har）;

%peak_b=max（fft1）

%rms_b=0.707*peak_b