matlab实现blob代码_多频铁耗系数的拟合含Matlab代码-CSDN博客

本文链接：https://blog.csdn.net/weixin_30950075/article/details/112404371

以有法为无法，以有限为无限，是为学术的最高境界。

目的和意义：

一般而言，为了在有限元计算过程中，更为精确地计算电机的铁耗（这只是理论上的假设，实际还是算不准），我们在有限元分析计算过程中可以采用如下公式来计及高频铁耗

其中，

、

和

为采用多频B-P曲线拟合后所得到的系数。

根据上述公式以及拟合思路，个人认为对上述三个系数可以有两种理解：

(1)从公式上来看，三个系数分别代表磁滞、涡流、杂散（非additional的杂散，此处为excessive）铁耗系数，是有物理意义的。

(2)从数值拟合角度来看，三个系数已经没有任何物理意义了，仅是拟合之后的数值解。

怎么理解见仁见智。

拟合方法：

多频B-P下的铁耗系数拟合一般有三种方法：1）ANSYS Maxwell自带的材料设置。2）Qucikfiled官方有一个Excel表格，可以方便计算（网址在这里，大家可以自行下载

Core loss coefficients calculatorquickfield.com

）。Core loss coefficients calculator）。3）Matlab自行编程。

这里我着重介绍一下Matlab自行编程，注意：这里我参考了SMEKLIB代码，大家可以自行查阅

https://github.com/AnttiLehikoinen/SMEKlib/blob/master/Examples/Core%20Loss%20Coefficient%20Fitting/fit_M330_35a.mgithub.com

下面放出我在其基础上修改过后的代码（非原创），个人认为应该比原代码更好理解一些，我在关键地方都给出了注释。

clear all
clc

B = [0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8]'; % flux density
Rho = 7600;  % density kg/m^3
P_init = [
0.02	0.05	0.12	0.33
0.08	0.2	0.48	1.27
0.17	0.41	1.02	2.69
0.28	0.67	1.68	4.49
0.4	0.97	2.47	6.66
0.53	1.3	3.37	9.19
0.68	1.68	4.39	12.11
0.84	2.1	5.54	15.44
1.02	2.56	6.82	19.22
1.22	3.07	8.25	23.54
1.44	3.64	9.86	28.48
1.69	4.29	11.6	34.12
2	5.07	13.7	40.62
2.4	6.06	16.3	48.24
2.94	7.4	19.6	57.86
3.67	8.86	23.2	70.24
4.32	NaN	NaN	NaN
4.73	NaN	NaN	NaN
]*Rho; % losses density, please remember the unit here is W/m^3, because W/kg * kg/m^3 = W/m^3.

freq = repmat( [50 100 200 400], [size(B),1]); % denfine the frequency, Hz

% P = (kh*f+kc*f^2)*B^2+ke*(f*B)^1.5, the following defintions do not
% contain kh, kc, ke. Usually, we can ignore the last term "excessive losses".
P_h = freq.*B.^2;
P_c = freq.^2 .* B.^2;
%P_e = freq.^1.5 .* B.^1.5; % This term can be added if needed.

% "" is called  mldivide, A*X = B, X = AB, it can be used for solving the
% least square, see "mldivide" in Matlab help document. Solve kh, kc, ke
% here.
k = [P_h(P_init>0) P_c(P_init>0) ]P_init(P_init>0 ); % solve the coefficients k

%recomputing losses
P_recomput = k(1)*P_h + k(2)*P_c; % in W/kg
P_recomput( P_recomput ==0 ) = 0;


figure('Position',[320, 230, 800, 450]);
b1 = plot (B(:,1),P_recomput);
set(b1,'LineWidth',1.5,'MarkerSize',8,'Marker','*');
hold on;
b2 = plot (B(:,1),P_init);
set(b2,'LineWidth',1.5,'MarkerSize',25,'Marker','.');
grid on;
set(gcf,'color','w');
set(gca,'FontSize',14,'FontName','Times New Roman');
xlabel('Flux density [T]');
ylabel('Core losses [W/m^3]');

运行结果如下：