run_winding.m

function run_winding()
% Generate the winding (litz wire) material data.
%
%    Map the different materials with a unique id.
%
%    (c) 2019-2020, ETH Zurich, Power Electronic Systems Laboratory, T. Guillod

% init
addpath(genpath('utils'))

% parse data
data = {};
data{end+1} = get_data('50um');
data{end+1} = get_data('71um');
data{end+1} = get_data('100um');
data{end+1} = get_data('200um');

% material type
type = 'winding';

% save material
save('data/winding_data.mat', '-v7.3', 'data', 'type')

end

function data = get_data(id)
% Generate the winding (litz wire) material data.
%
%    Parameters:
%        id (int): material id
%
%    Returns:
%        data (struct): material id and data

% get values
switch id
    case '50um'
        fill_litz = 0.47;
        d_strand = 50e-6;
        kappa_copper = 32.5;
    case '71um'
        fill_litz = 0.49;
        d_strand = 71e-6;
        kappa_copper = 23.5;
    case '100um'
        fill_litz = 0.51;
        d_strand = 100e-6;
        kappa_copper = 21.5;
    case '200um'
        fill_litz = 0.53;
        d_strand = 200e-6;
        kappa_copper = 19.0;
    otherwise
        error('invalid material id')
end

% maximum frequency (radius is equal to skin depth)
mu0_const = 4.*pi.*1e-7;
sigma = 5.8e7;
delta = 0.5.*d_strand;
f_max = 1./(pi.*mu0_const.*sigma.*delta.^2); % maximum frequency

% conductivity interpolation
material.interp.T_vec = [20 46 72 98 124 150]; % temperature vector
material.interp.sigma_vec = 1e7.*[5.800 5.262 4.816 4.439 4.117 3.839]; % conductivity vector

% assign param
material.param.fill_litz = fill_litz; % fill factor of the litz wire itself (not of the packing)
material.param.d_strand = d_strand; % strand diameter
material.param.f_max = f_max; % maximum frequency

% assign density
material.param.rho_copper = 8960; % volumetric density for copper
material.param.rho_iso = 1500; % volumetric density for insulation
material.param.kappa_iso = 5.0; % cost per mass for the insulation
material.param.kappa_copper = kappa_copper; % cost per mass for the copper

% assign constant
material.param.P_max = 1000e3; % maximum loss density
material.param.J_rms_max = 10e6; % maximum rms current density
material.param.P_scale_lf = 1.0; % scaling factor for LF losses
material.param.P_scale_hf = 1.0; % scaling factor for HF losses
material.param.T_max = 130.0; % maximum temperature
material.param.c_offset = 0.3; % cost offset

% assign
data = struct('id', id, 'material', material);

end