matlab / отчет1

5

Завдання. Складіть структурну схему АД із короткозамкненим ротором у фазних координатах і виконайте моделювання за умови повної симетрії АД.

cur_calc.m function Io = cur_calc(u) global Lm Ess Err PSIo=u(1:6); tetar=u(7); c1=cos(tetar); c2=cos(tetar+2*pi/3); c3=cos(tetar+4*pi/3); L=[Ess -1/2 -1/2 c1 c2 c3; -1/2 Ess -1/2 c3 c1 c2; -1/2 -1/2 Ess c2 c3 c1; c1 c3 c2 Err -1/2 -1/2; c2 c1 c3 -1/2 Err -1/2; c3 c2 c1 -1/2 -1/2 Err]; L=2/3*Lm*L; PSIo=PSIo(:); Io=L\PSIo; tor_calc.m function M_em = tor_calc(u) global Zp Lm u=u(:)'; D=[0 1 -1; -1 0 1; 1 -1 0]; PSIr=u(1:3)'; Ir=u(4:6); M_em=Zp*Ir*D*PSIr/sqrt(3); init_param.m close all, clc global Lm Ess Err Zp Pn=2.2; Nn=875; Ix=4.6;J=0.028; Zp=3; fsn=50; Uslen=380; Rs=3.6; Xs=2.58; Is=7.2; Rr=5.7; Xr=2.63; Ir=4.5; R=diag([Rs Rs Rs Rr Rr Rr]); Usfen=Uslen/sqrt(3); wsn=2*pi*fsn; Xm=(Usfen-Ix*sqrt(Rs^2+Xs^2))/Ix; Lm=Xm/wsn; bs=Xs/Xm; br=Xr/Xm; Ess=1+3/2*bs; Err=1+3/2*br; Mn=Pn*1000/wsn;

plot_results.m figure(1) plot(t,[w,M]),grid on,hold on legend('w','M') xlabel('t,sec') ylabel('w,rad/sec, M,N*m') figure(2) plot(t,PSIs_ABC),grid on,hold on legend('faza A','faza B','faza C') xlabel('t,sec') ylabel('PSIs,Vb') figure(3) plot(t,Is_ABC),hold on,grid on legend('faza A','faza B','faza C') xlabel('t,sec') ylabel('Is,A') figure(4) plot(t,PSIr_abc),hold on,grid on legend('faza a','faza b','faza c') xlabel('t,sec') ylabel('PSIr,Vb') figure(5) plot(t,Ir_abc),hold on,grid on legend('faza a','faza b','faza c') xlabel('t,sec') ylabel('Ir,A')

Рисунок 1.2 - Структурна схема АД у фазних координатах

Статор

Ротор