Этот пример показывает вычисление и подтверждение нелинейной характеристики намагничивания ядра преобразователя. Начиная с основных значений параметров выведена базовая характеристика. Это затем используется в модели Simscape™ схемы тестирования в качестве примера, которая может использоваться, чтобы построить базовую характеристику намагничивания на осциллографе. Образцовые выходные параметры затем по сравнению с известными значениями.
Основные значения параметров, используемые в качестве основания для последующих вычислений:
Проницаемость свободного пространства,
Относительная проницаемость ядра,
Количество первичных поворотов,
Количество вторичных поворотов,
Эффективная длина магнитного сердечника,
Эффективная площадь поперечного сечения магнитного сердечника,
Базовое насыщение начинается,
Ядро, полностью насыщаемое,
mu_0 = pi*4e-7; mu_r = 3000; N1 = 23; N2 = 29; le = 0.032; Ae = 1.6e-5; B_sat_begin = 0.6; B_sat = 1.2;
Где:
Плотность магнитного потока,
Сила магнитного поля,
Линейное представление:
Нелинейное представление (включая коэффициент, 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; % Nonlinear representation H_nonlinear = 0:25:750; B_nonlinear = B_sat*tanh(a*H_nonlinear);
Вычисленные параметры могут теперь использоваться в модели Simscape. После того, как моделируемый, модель собирается вывести переменную журналирования Simscape, simlog, и некоторые сигналы с помощью выходных портов, вашего. Параметры схемы:
Исходное значение напряжения,
Исходная частота напряжения,
Исходное сопротивление напряжения,
Операционный усилитель ввел сопротивление,
Сопротивление обратной связи операционного усилителя,
Емкость обратной связи операционного усилителя,
% Circuit parameters Vs = 10; Freq_Hz = 60; R_Vs = 10; R_1 = 1e3; R_2 = 1e6; C_2 = 1e-6; % Open model and simulate modelName = 'ee_nonlinear_transformer'; open_system( modelName ); simOut = sim( modelName ); yout = get(simOut, 'yout'); simlog = get(simOut, 'simlog'); % Collect internal Simscape logging data for comparison I_simscape = simlog.Nonlinear_Transformer.Lm.i.series.values; phi_simscape = simlog.Nonlinear_Transformer.Lm.phi.series.values; % Collect model output data for comparison (as used for oscilloscope) I_primary = yout(:,1); int_V_secondary = yout(:,2);
Данные должны быть обработаны, чтобы обеспечить силу Магнитного поля и данные о плотности Магнитного потока для сравнения. Где:
Магнитодвижущая сила,
Магнитный поток,
Входное напряжение операционного усилителя,
Выходное напряжение операционного усилителя,
Уравнения, которые будут использоваться, следующие:
% Internal logging data H_simscape = I_simscape.*N1./le; B_simscape = phi_simscape./Ae; % Oscilloscope scaling and model output data H_measured = I_primary.*N1./le; phi_measured = (int_V_secondary.*R_1.*C_2)./N2; B_measured = phi_measured./Ae;
Эти три характеристики могут теперь быть наложены:
Заданная характеристика: вычисленный от основных параметров
Характеристика от журналирования: вычисленный от внутренних данных о регистрации Simscape
Характеристика от измерения: полученный измерением и вычислением с помощью электронной схемы тестирования
Из-за утечки и паразитных параметров, характеристика, полученная из электронной схемы тестирования, отличается к заданной характеристике. Однако схема тестирования и ее параметризация, как показывают, находят характеристику для данного преобразователя в подходящих допусках.
figure, plot( ... H_nonlinear,... B_nonlinear,... 'o',... H_simscape,... B_simscape,... H_measured,... B_measured,... '--'... ); 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( 'Defined characteristic', 'Characteristic from logging',... 'Characteristic from measurements', 'Location', 'NorthWest' );