Inductor Design

clc;clear all;
%% units & constant
mT = 10^-3;mm2 = 10^-6;mm = 10^-3;u0 = 4*pi*10^-7;mm3 = 10^-9;

% Core size/basic parameters
Bsat = 400*mT;
Ve = 3310*10^-3;%(cm3)
Ae = 58*mm2;
Le = 57*mm;
ue = 2000*u0;
As = 8; % (cm2)

%% Design parameters
N = 6;
Imax = 15;
Irms = 12;
Lg = 7*0.0889*mm;      % real gap distance = 2*Lg

%% calculation
Rc = Le/ue/Ae;
Rg = Lg*2/u0/Ae;
Rtotal = Rc+Rg;

% constraint #1
Phi = N*Imax/Rtotal;
Bcal = Phi/Ae;

% constraint #2
L = N^2/Rtotal*10^6;  % unit uH

%% Loss Estimation
% core loss
Cm = 2e-5;x = 1.8;y = 2.5;Ct2 = 0.77e-4;Ct1 = 1.05e-2;Ct0 = 1.28;
T = 80;
Ct = Ct0-Ct1*T+Ct2*T*T;
f = 500*10^3;
% Pcore_pV = Cm*Ct*f^x*Bcal^y;
 Pcore_pV = 400;
Pcore = Pcore_pV*Ve*0.001;

% winding loss
Rdc = N*(2*10^-8*4*(10/1000)/(pi*(0.04/2000)^2))/1162;
Rac = Rdc*1.5;
Pcu = Rac*Irms*Irms;

Ptot = Pcu + Pcore;
dTtot = (Ptot*1000/As)^0.833;