?? magckt2.m
字號:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% magckt2.m - Analysis of Fig 3.40 magnetic circuit.
% Calculates flux linkage & current arrays
% for specified value of air gap length x.
% Leakage flux is considered in calculations.
% Used by soldes.m.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [LamI] = magckt1(x);
soldata,
if x==0; x=1e-12; else; end; del=0.01*x;
% Air gap & leakage permeances
Pg=4e-07*pi*(pi*d^2/4/x+0.52*pi*d);
Pl=4e-07*pi^2*(d+db)*lp/(db-d)/2;
% Generate left & right flux-mmf curves
npts=25; phimax=Bmax*pi*d^2/4; phi=linspace(0,phimax,npts);
for i=1:npts
FL(i)=l1*Hcs(phi(i)/A1)+l2*Hcs(phi(i)/A2)+l3*Hcs(phi(i)/A3);
phiL(i)=phi(i)+FL(i)*Pl;
FR(i)=l6*Hcs(phi(i)/A6)+l7*Hcs(phi(i)/A7)+l8*Hcs(phi(i)/A8)+ ...
l9*Hcs(phi(i)/A9);
phiR(i)=phi(i)+FR(i)*Pl;
end
% Solve center loop looking up FL & FR at same values of flux
for i=1:npts
mmfR=interp1(phiR, FR, phi(i));
mmfL=interp1(phiL, FL, phi(i));
F(i)=l4*Hcs(phi(i)/A4)+l5*Hcs(phi(i)/A5)+(lp/2-x)* ...
Hcs(phi(i)/A10)+phi(i)/Pg+mmfR+mmfR;
LamI(1,i)=N*phi(i); LamI(2,i)=F(i)/N;
philR(i)=mmfR*Pl; philL(i)=mmfL*Pl;
end
% Activate to plot leakage & plunger fluxes
% Can handle from command window: magckt1(x)
% plot(F/N,phi,F/N,philL,F/N,philR); grid;
% ylabel('Flux, Wb'); xlabel('Coil current, A');
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -