Exponenta Banner

Нелинейные характеристики индуктора

Открыть скрипт

Этот пример показывает сравнение нелинейного поведения индуктора для различной параметризации. Начиная с основных значений параметров выведены параметры для линейных и нелинейных представлений. Эти параметры затем используются в модели Simscape™ и симуляции сравненные выходные параметры.

Спецификация параметров

Основные значения параметров, используемые в качестве основания для последующих вычислений:

  • Проницаемость свободного пространства,

  • Относительная проницаемость ядра,

  • Количество обмотки поворотов,

  • Эффективная длина магнитного сердечника,

  • Эффективная площадь поперечного сечения магнитного сердечника,

  • Базовое насыщение начинается,

  • Ядро, полностью насыщаемое, 

mu_0 = pi*4e-7;
mu_r = 3000;
Nw = 10;
le = 0.032;
Ae = 1.6e-5;
B_sat_begin = 0.75;
B_sat = 1.5;

Вычислите данные о силе плотности и магнитного поля магнитного потока

Где:

  • Плотность магнитного потока,

  • Сила магнитного поля,

Линейное представление:

Нелинейное представление (включая коэффициент, a):

% Use linear representation to find value of H corresponding to B_sat_begin
H_sat_begin = B_sat_begin/(mu_0*mu_r);
% Rearrange nonlinear representation to calculate coefficient, a
a = atanh( B_sat_begin/B_sat )/H_sat_begin;

% Linear representation
H_linear = [-500 500];
B_linear = mu_0*mu_r*H_linear;

% Nonlinear representation
H_nonlinear = -5*H_sat_begin:H_sat_begin:5*H_sat_begin;
B_nonlinear = B_sat*tanh(a*H_nonlinear);

Отобразите плотность магнитного потока по сравнению с силой магнитного поля

Линейные и нелинейные представления могут быть наложены.

figure,plot( H_linear, B_linear, H_nonlinear, B_nonlinear );
grid( 'on' );
title( 'Magnetic flux density, B, versus Magnetic field strength, H' );
xlabel( 'Magnetic field strength, H (A/m)' );
ylabel( 'Magnetic flux density, B (T)' );
legend( 'B=\mu_0 \mu_r H', 'B-H curve', 'Location', 'NorthWest' )

Вычислите магнитный поток и текущие данные

Где:

  • Магнитный поток,

  • Текущий,

Линейное представление:

Нелинейное представление:

% Linear inductance
L_linear = mu_0*mu_r*Ae*Nw^2/le;

% Linear representation
I_linear = [-1.5 1.5];
phi_linear = I_linear.*L_linear/Nw;

% Nonlinear representation
I_nonlinear = le.*H_nonlinear./Nw;
phi_nonlinear = B_nonlinear.*Ae;

Отобразите магнитный поток по сравнению с текущим

Линейные и нелинейные представления могут быть наложены.

figure, plot( I_linear, phi_linear, I_nonlinear, phi_nonlinear );
grid( 'on' );
title( 'Magnetic flux, \phi, versus current, I' );
xlabel( 'Current, I (A)' );
ylabel( 'Magnetic flux, \phi (Wb)' );
legend( '\phi=I L/N_w', '\phi-I curve', 'Location', 'NorthWest' );

Используйте параметры в модели Simscape

Вычисленные параметры могут теперь использоваться в модели Simscape. После того, как моделируемый, модель собирается создать переменную журналирования Simscape, simlog.

modelName = 'ee_nonlinear_inductor';
open_system( modelName );
sim( modelName );

Заключение

Переменная состояния для всех представлений является магнитным потоком. Текущий, я и магнитный поток, можем быть получены из переменной журналирования Simscape, simlog, для каждого представления. Накладывание результатов симуляции от представлений разрешает прямое сравнение.

figure, plot( ...
    simlog.Linear_Inductor.inductor.i.series.values,...
    simlog.Linear_Inductor.inductor.phi.series.values,...
    simlog.B_vs_H.inductor.i.series.values,...
    simlog.B_vs_H.inductor.phi.series.values,...
    simlog.phi_vs_I.inductor.i.series.values,...
    simlog.phi_vs_I.inductor.phi.series.values,...
    'o'...
    );
grid( 'on' );
title( 'Magnetic flux, \phi, versus current, I' );
xlabel( 'Current, I (A)' );
ylabel( 'Magnetic flux, \phi (Wb)' );
legend( 'Linear (single inductance)', 'B-H characteristic', '\phi-I characteristic', 'Location', 'NorthWest' );

Документация Simscape Electrical
Поддержка