pdfs

function [ N_turn , l_air_gap , B_max ] = Inductor_Winding_Calculation_Function( L , parameter , A_core , l_core , I_max , B_sat , method )
%% inputs
% L is the desired inductance
% A_core is the cross-section area of the core
% l_core is the length of the magnetic loop
% I_max is the maximum current
% B_sat is the saturation flux density
% parameter and method should correspond to each other:
% if the parameter is the mu_r of the core material,
% then method should be 'mu_r';
% if the parameter is the A_L of the core material,
% then method should be 'A_L'.
%% outputs
% N_turn is the result number of turns
% l_air_gap is the result air gap length
% B_max is the result maximum flux density
%%
mu_0 = 4 * pi * 1e-7;
switch method
    
    case 'mu_r'
        
        mu_r = parameter;
        
        % Magnetic Permeability
        mu = mu_0 * mu_r;
        % Reluctance
        R_m = l_core / mu / A_core;
                
    case 'A_L'
        
        A_L = parameter;
        R_m = 1 / A_L;
        mu = l_core / R_m / A_core;
        mu_r = mu / mu_0;
        
    otherwise
        
end
% Number of turns
N_turn = sqrt(L*R_m);
        
% Check Magnetic Saturation
B_max = N_turn * I_max / R_m / A_core;
if B_max > B_sat    % if saturated, then insert air gap
    N_turn = I_max * L / B_max / A_core;
    R_m = N_turn^2 / L;
    l_air_gap = mu_0 * A_core * R_m - l_core / mu_r;
else                % if not saturated, then no need to insert air gap
    l_air_gap = 0;
end
end