laoyoung 发表于 2017-5-17 17:04

振动加速度频域分析

如题,我在做振动加速度的振动测试,得到的测试结果进行fft分析后加速度幅值最大值出现在800hz处,转速为3000转每分钟,8极9槽电机,SVPWM供电,空载。
一般在文献中看到的是2倍电频率处即400HZ,对电机模型电磁力进行fft变换,得到的是400hz出出现最大值,其次是800HZ。电磁力的 fft分析和振动fft分析如图

请问,为什么800HZ对应的振动幅值这么大

laoyoung 发表于 2017-5-17 22:15

现将自己做的展示部分内容
xx=rand(500,1001);yy=rand(250,500);
for i=1:500
   for j=1:1001
    xx(i,j)=rmf1000(1001*(i-1)+j);
   end
end
figure;
surf(xx);
y=abs(fft2(xx,500,1001))/500000*4;
z=flipud(y);
% y(1,1)=0;
% y=flipud(y);
figure;
surf(y);
% for i=1:500
%    for j=1:1001
%   if y(i,j)<1
%      y(i,j)=0;
%   end
%    end
% end
for i=1:250
   for j=1:500
      yy(i,j)=y(i+1,j+1);
   end
end
figure;
surf(yy);
% axis();
% set(gca,'YDir','reverse');
rmf1000是直接导入三维CSV文件
xx是电磁力空间和时间分布数组,yy是经过fft2变化后处理得到的最终结果
MAXWELL导出的电磁力csv文件,就是电磁力那一栏,是一个一维数组,我处理下,变成FFT2可以使用的数组
MATLAB里的FFT函数,对二维波形处理得到的图像,去掉直流分量是对称的,因此用fft2处理后,去掉第一行和第一列,取数组的左上四分之一,得到最终的结果,然而左上四分之一是与右下四分之一对称的,右上四分之一和坐下四分之一对称,左上和左下不相同,我不太懂三维FFT的原理,所以就如上处理了,有意见的请指正




njclark 发表于 2017-5-18 08:48

你可能搞混淆了,首先从你下面那张图,可得到信息:16th电磁引起的振动和1th转子引起的振动,两个频率上比较大。
其次,你提到“文献中看到的是2倍电频率处即400HZ”是不是指fft结果的对称性,通常只使用前半部分的结果,即使用采样频率一半的结果?

laoyoung 发表于 2017-5-18 09:02

前面那句话我没太明白,16th和1th分别是800HZ和50HZ处,后者是机械频率,在我的研究中,机械频率忽略掉,这里的问题是?
后面那句话,fft结果有对称性,确实采用的是一半的结果,不然的话,除了400HZ处还有另一处的幅值和400HZ处的幅值是相等的,二者呈对称状态,你的意思是说400HZ其实代表的是800HZ,400Hz处出现电磁力的最大值其实在maxwell自带fft分析结果里也是这样。
不好意思,我其实没太弄明白你的问题

laoyoung 发表于 2017-5-19 23:03

假如我要削弱51.6683那个径向电磁力,16阶力波有两种情况,定子4次和转子20次谐波相互作用,定子-20和转子4次相互作用,对于后者是不是只要减小转子永磁体四次谐波幅值就行?

lijian988998 发表于 2017-11-14 00:32

您好,请问 您的电磁力2维fft 做出来了么?

超级职业英雄 发表于 2017-11-14 01:11

lijian988998 发表于 2017-11-14 00:32
您好,请问 您的电磁力2维fft 做出来了么?

哈哈哈哈哈哈,看到你了

lijian988998 发表于 2017-11-14 01:49

超级职业英雄 发表于 2017-11-14 01:11
哈哈哈哈哈哈,看到你了

您好!请问您是?

超级职业英雄 发表于 2017-11-14 14:02

lijian988998 发表于 2017-11-14 01:49
您好!请问您是?

看我头像,这下认识了吗?哈哈哈
页: [1]
查看完整版本: 振动加速度频域分析