Расчет обмоток трансформатора в MatLab
Добрый день.
Возникла проблема с контрольной работой. В задании необходимо найти индуктивности первичной и вторичной обмоток L1 и L2. Нашёл на просторах код для задания расчёт как и нужно методом объединённых матриц. Только, как я понял, расчитана общая индуктивность КЗ ZKZ катушек. Подскажите, как расчитать индуктивность L1 и L2?
%<<const>>
E=300; % напряжение подаваемое на первичную обмотку трансформатора
F=50; % частота питающей сети
R1=1; % активное сопротивление на 1-ой обмотке
R2=1; % активное сопротивление на 2-ой обмотке
W1=1200; % количество витков на 1-ой обмотке
W2=110; % количество витков на 2-ой обмотке
LM=0.1; % длина магнитопровода
SM=0.01; % площадь поперечного сечения магнитопровода
k=W2/W1; % коэффициент трансформации
r2=R2/(k); % приведённое активное сопротивление 2-ой обмотки
M=310; % относительная магнитная проницаемость
M0=pi*4e-7; % абсолютная магнитная проницаемость
Rn=440; % сопротивление нагрузки
RN=1e10; % активное сопротивление нагрузки в режиме холостого хода
LN=0.004; % индуктивность нагрузки
GEE=[1 0; 0 1];%
GMM=[1];
EV=[E; 0]; % матрица ЭДС
WEM=[1200; 110];
ZEVxx=diag([R1 R2+RN]); % матрица сопротивлений при холостом ходе
ZEVkz=diag([R1 R2]); % матрица сопротивлений при коротком замыкании
j=sqrt(-1); % мнимая единица
ZVM=[LN/(M*SM*M0)]; % ветвевивые магнитные сопротивлений
% Режим короткого замыкания
IVkz=1; % ток вторичной обмотки
ZEK=GEE*ZEVkz*GEE'; % контурное электрическое сопротивление
ZEMK=j*F*GEE*WEM*GMM'; % контурное электромагнитное сопротивление
ZMEK=-GMM*WEM'*GEE';
ZMK=GMM*ZVM*GMM'; % контурное магнитное сопротивление
EK=GEE*EV; % контурные ЭДС
DELTAZEK=-ZEMK*inv(ZMK)*ZMEK; % Внесенное сопротивление
IK=inv(ZEK+DELTAZEK)*EK; % матрица контурных токов
IV=GEE'*IK; % матрица ветвевых токов
UV=ZEVkz*IV+EV;
ZKZ=UV(1)/IVkz
Комментарии
Боюсь соврать, но помнится там половина бралась от общей.
Здравствуйте. Я где-то вот такую формулу нашёл L1=(1.26*M*W1^2*SM)/(10^8*0.1). Не знаю на сколько она применима к моей работе. В этом случае получается число. ZKZ MathLab выдаёт комплексным числом. Да и вообще, почему-то много за ним пересчитывать нужно. Коэффициент усиления при w1/w2 он считает как 0,09, хотя это 10,9.
Может, всё-таки нужно проверить, что вы сами написали?
Вы написали k=W2/W1. Матлаб вам и посчитал k=110/1200=0.9.
Вы правы. Уже совсем глазками не вижу.
@Centurio
Два вопроса:
Это имеете ввиду, что L1 = 0.91L, L2 = 0.09L ?
Это какие индуктивности, первичные/вторичные?
а может вообще рассеяния. а так L1=N^2*L2
А L2-? N - это сопротивление нагрузки, или количество витков?
В схеме же замещения Lm L1 L2 , N коэф трасформации.
Какая-то каша
L1 = Lkz(N^2/(N^2+1))
L2 = Lkz(1/(N^2+1))
так должно быть?
Вроде как да. Проверьте по учебнику. Вам проще он у вас под рукой (теоретически). Не у всех имеются под рукой книги по электрическим машинам. Тем более в нынешней пандемической ситуации, кто-то на удаленке.
Спасибо за помощь. Книги под рукой у меня нет.