?????????????? ?????? ? ????????? ?????????? ??????? ?????

• окончательно узловые напряжения определяются суммированием ранее
найденного вектора узловых напряжений и приращения dUo:
Uo ~ =Uo+dUo.
Найденные узловые напряжения в масштабе 1:1 соответствуют индукции
магнитного поля.
Исходные данные для расчета следующие:
parp – число пар полюсов; interval – число интервалов разбиения пазового
шага;
w – число витков фазной обмотки; chered – чередование фаз в пазах
индуктора;
Imax - амплитуда тока фазной обмотки; f – частота питания фазной
обмотки;
tau – полюсное деление фазной обмотки;aplast – толщина пластины;
ro – удельная проводимость пластины; delta – воздушный зазор;
h – половина расчетного пазового шага;
Vd – скорость движения пластины; Nt - число расчетных точек по времени
ta,he – входные и выходные узлы ветвей подсхемы;
C,A – топологические матрицы подсхем; Сobm – обмоточная матрица;
Сos – матрица цепи пересечений; Uo – вектор узловых напряжений
схемы;
B - вектор индукции в узлах решетки;
clear;
parp=4; interval=2; pazov=6*parp*interval; w=30/interval; chered=[2 1 3];
Imax=300; f=25; tau=0.2; aplast=5e-3; ro=1/20.0e6; mu0=4.0e-7*%pi; delta=0.01;
h=tau/(6*interval); Ki=2.*aplast*h/ro; v=15; o=11;
Vd=7; Nt=500; // задание скорости и числа расчетных точек по времени
ta= [3 4 5 6 7 8 9 10 11 2 3 5 7 9 11];
he=[1 3 4 5 6 7 8 9 10 11 12 12 12 12 12];
[C,A,g]=formc(ta,he);
//g('node_x')=[100 1100 200 300 400 500 600 700 800 900 1000 600];
//g('node_y')=[200 200 200 200 200 200 200 200 200 200 200 100];
//show_graph(g);
i=0;k=0;
for j=1:interval:pazov
k=k+1; i=i+1; if i==4 then i=1; end;
m=w*(-1)^k;
for n=0:interval-1
Cobm(j+n,chered(i))=-m; //формирование обмоточной матрицы
end;
end;
kp=int(2*pazov/5); oo=o*kp; ksv=kp+1;
Cos(2*kp+2,ksv)=0; //объявление размера матрицы цепи пересечений
dt=0.1*h/Vd; // вычисление расчетного шага по времени
for i=1:ksv
j=2*i; Cos(j-1,i)=-1; Cos(j,i)=1; // формирование матрицы цепи пересечений
if i