Imchar.m
上传用户:eighthdate
上传日期:2014-05-24
资源大小:270k
文件大小:1k
源码类别:

其他行业

开发平台:

Matlab

  1. % This program called by im produces the torque/speed curve
  2. % for an induction motor.
  3. %
  4. % Copyright (C) 1998 by H. Saadat
  5. slip1 = Smax:-.005:.0001;
  6. if Smax > 1, Smax =1;else,end
  7. slip2 = 1:-.005:Smax;
  8. slip=[slip2, slip1];
  9. np1 = length(slip1);
  10. np2 = length(slip2);
  11. R2d1 = R2./slip1;
  12. R2d2 = R2./slip2;
  13. Z2d1=R2d1+ j*X2*ones(1,np1);
  14. Z2d2=R2d2+ j*X2*ones(1,np2);
  15. I2d1 = Vth./(Zth*ones(1,np1)+Z2d1);
  16. I2d2 = Vth./(Zth*ones(1,np2)+Z2d2);
  17. Pgd1 = 3*R2d1.*abs(I2d1).^2;
  18. Pgd2 = 3*R2d2.*abs(I2d2).^2;
  19. Tdev1 = Pgd1/Ws;
  20. Tdev2 = Pgd2/Ws;
  21. nr1 = Ns*(ones(1,np1)-slip1);
  22. nr2 = Ns*(ones(1,np2)-slip2);
  23. if length(nr2) > 4
  24.    x1= 1.7*nr2(4);
  25.    else x1 = 1.7*nr1(2); end
  26. y1 = 0.86*Tst;
  27. x2 = 0.96*(1- Smax)*Ns; y2 = 0.89*Tmax;
  28. x3 = 0.95*Nr;  y3 =0.93*Td;
  29. fy = 0:Td/20:Td; nf=length(fy);
  30. fx = ones(1,nf)*Nr;
  31. figure
  32. plot(nr1, Tdev1, '-g', nr2, Tdev2, 'r', fx, fy,'b'),grid
  33. xlabel('Rotor speed - rpm')
  34. ylabel('Developed Torque  Nm')
  35. title('Motor torque-speed curve')
  36. text(x1, y1, 'Ts')
  37. text(x2, y2, 'Tmax')
  38. if Nr ~= 0
  39.   text(x3, y3, 'Td')
  40. else, end