Attention! Your ePaper is waiting for publication!
By publishing your document, the content will be optimally indexed by Google via AI and sorted into the right category for over 500 million ePaper readers on YUMPU.
This will ensure high visibility and many readers!
• окончательно узловые напряжения определяются суммированием ранее найденного вектора узловых напряжений и приращения 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
end; k=0; m=0; j=1; Cdif(v,pazov,kp)=0; //объявление дифференцирующей матрицы Cdif(14,1,j)=0.5; Cdif(15,2,j)=0.5; k=0;j=j+1; for i=3:pazov+2 //формирование дифференцирующей матрицы k=k+1; if k==6 then k=1;j=j+1; end; if (i>pazov) then Cdif(10+k,i-2,j)=-0.5; else Cdif(10+k,i,j)=0.5; Cdif(10+k,i-2,j)=-0.5; end; end; for i=1:15 if i=2)&(npaz