إستخدام octave لحساب العلاقة بين العزم والسرعة في الموتور الحثي

السلام عليكم ورحمة الله وبركاته

http://www.ece.ualberta.ca/~knight/electrical_machines/induction/basics/trq_speed.html

إذا يوجد أي إستفسار مرتبط بالموتور الحثي فأهلا وسهلا

وهذا هو الكود

https://bitbucket.org/xmonader/torque-speed-ccs/src/3c3923a3df63831a82d29689948da0ba7ee8c442/tsccs.m?at=default

 
 

%% Author  : xmonader [at] gmail.com
%% License : GPL v3
%% Site URL: 
%% moreinfo: http://www.ece.ualberta.ca/~knight/electrical_machines/induction/basics/trq_speed.html 

%Calculating Torque
function t = Torque(Vth, Rth, R2, s, Xth, X2, Ws)
	 t=(3*Vth*Vth*R2)/(( (Rth+R2/s)^2 + (Xth+X2)^2 )*s*Ws)
endfunction

%Calculating regular parameters <Vth,Rth,Xth>
function [Vth, Rth,Xth] = calcparams(Vph,r1,x1,Xm)

	Vth=Vph*(Xm/sqrt(r1^2 + (x1+Xm)^2));
	Zth=( (j*Xm)*(r1+j*x1))/(r1+j*(x1+Xm));
	Rth=real(Zth);
	Xth=image(Zth);
endfunction

%Generic function for plotting the curve
function plotcurve(Nm, t)
        figure(1)
        plot(Nm, t, 'k-', 'LineWidth', 2.0);
        xlabel('S');
        ylabel('T');
        title('IM T-S c/cs');
endfunction

%The main function
%Use it with the motor data
function plotfor(r1,r2,x1,x2,Xm,Vph, Ns,Ws)

        [Vth, Rth,Xth]= calcparams(Vph, r1, x1, Xm);
         
        %T-S curve for many slips
        s = (1:50)/50; 
        s(1) = 0.001;
        Nm = (1-s)*Ns;
        t = zeros(1,50);
        for  i = 1:50
                t(i) = Torque(Vth,Rth,r2,s(i),Xth,x2,Ws); 
        end

        plotcurve(Nm, t)

endfunction


%%a few examples
function test()
        r1=0.200;
        r2=.120;
        x1=x2=0.410;
        Xm=15;
        Vph=208/sqrt(3);
        Ns=3600;
        Ws=377;
        plotfor(r1,r2,x1,x2,Xm,Vph, Ns, Ws)
endfunction

function p1()

        r1 = 0.641;                 
        x1 = 1.106;                 
        r2 = 0.332;                 
        x2 = 0.464;                 
        Xm = 26.3;                  
        Vph = 460 / sqrt(3);    
        Ns = 1800;              
        Ws = 188.5;             
        plotfor(r1,r2,x1,x2,Xm,Vph, Ns, Ws)
endfunction

function p2()
        r1 = 0.641;               
        x1 = 0.750;                
        r2 = 0.300;                
        x2 = 0.500;              
        Xm = 26.3;                 
        Vph = 460 / sqrt(3);    
        Ns = 1800;             
        Ws = 188.5;
        plotfor(r1,r2,x1,x2,Xm,Vph, Ns, Ws)
endfunction