电机学Matlab仿真代码

最新推荐文章于 2024-03-06 21:58:28 发布

PANTHE

最新推荐文章于 2024-03-06 21:58:28 发布

阅读量1.7k

点赞数 3

分类专栏：仿真文章标签： matlab 开发语言算法

本文链接：https://blog.csdn.net/jiashuchampion/article/details/127980312

版权

仿真专栏收录该内容

3 篇文章 0 订阅

订阅专栏

仿真代码目录

资源链接在文章底部

一、变压器仿真

1.1电压供电时气隙大小对线电流和磁涌的影响

%气隙大小对线圈电流的影响
%clear all
%close all
clc
a=0.03;
b=0.03;
f=50;
u=500;
N=500;
l=0.15*4;
u0=4*pi*10^(-7);
uFe=5000*u0;
x=10^(-4):10^(-6):6*10^(-4);
phi1=u /(4.44*N*f);
y=((l-x)./(a.*b.*uFe)+x./(a.*b.*u0)).*phi1/N;
y1=((l-x)./(a.*b.*uFe)+x./((x+a).*(x+b).*u0)).*phi1/N;
figure
plot(x,y,'-k',x,y1,'--k');
xlabel('δ / m')
ylabel('I / A')
legend('忽略边缘效应','不忽略边缘效应')
grid on

仿真结果

%气隙大小对线圈电流的影响
clear;
a=0.03;
b=0.03;
f=50;
u=500;
N=500;
l=1;
u0=4*pi*10^(-7);
uFe=5000*u0;
x=10^(-4):10^(-6):6*10^(-4);
z=u./(4.44.*N.*f);
y=(l./(a.*b.*uFe)+x./((x+a).*(x+b).*u0)).*z/N;
y1=(l./(a.*b.*uFe)+x./((x+a).*(x+b).*u0)).*z/N;
plot(x,y,'-k',x,y1,':k');
xlabel('δ/m')
ylabel('I/A') 
legend('y','y1') 
grid on
pause
clf

%气隙大小对线圈磁通的影响
clear;
I=20;               %电流20A
U=500;              %电压500V
N=500;              %线圈匝数500匝
u0=pi*4e-7;         %真空磁导率μ0
uFe=5000*u0;        %铁芯磁导率μFe
S1=3*3e-4;          %铁芯面积3*3e-4 平方米
S2=3*3e-4;          %气隙面积3*3e-4 平方米
l=0.150*4;          %铁芯磁路长度150mm*4
Rm1=l/(uFe*S1);     %铁芯磁阻
x=1e-4:1e-6:6e-4;   %设置气隙长度变化范围，步长为1e-6 m
Rm2=x/(u0*S2);      %气隙磁阻
F=N*I;              %磁动势
Rm=Rm1+Rm2;         %铁芯和气隙磁阻之和
phi=F./Rm;          %线圈磁通φ
grid on
N=500;              %线圈匝数
I=20;               %线圈电流
u0=4*pi*10^-7;      %真空磁导率
u=5000*u0;          %铁心的磁导率
L=0.15;             %铁芯每段长度
Delta=0.0001:0.000001:0.0006;   %气隙长度
S1=0.03*0.03;             %铁芯截面积
S2=(0.03+Delta).*(0.03+Delta);         %考虑边缘效应的修正气隙处截面积
Rm1=(4*L)/(u*S1);   %除气隙外其余各部分的磁阻
Rm2=(Delta)./(u0.*S2);%气隙处磁阻
Rm=Rm1+Rm2;         %总磁阻
Flux=(N*I)./Rm;      %磁通
plot(x,phi,'k',Delta,Flux,':k');        %绘制图形
legend('phi','Flux')
axis([1e-4 6e-4 0.015 0.055]);%设置坐标范围
xlabel('δ/m');
ylabel('φ/Wb');
grid on

仿真结果