trabalho de controle

Prévia do material em texto

IFMA
CAMPUS SÃO LUÍS – MONTE CASTELO
DEPARTAMENTO DE ELETROELETRÔNICA - DEE
ENGENHARIA ELÉTRICA INDUSTRIAL
CONTROLE ANALÓGICO
 					MOTOR CC 
Prof. Reginaldo Miranda
São Luís
2017
Representação dinâmica do Motor CC é representada pelo diagrama de blocos abaixo:
Implementando esse diagrama de blocos em um programa de simulação (SIMULINK/MATLAB), podemos verificar o comportamento dinâmico do motor. Seja um motor de corrente contínua de 50V; Ra= 2,903Ω; La=2,58 mH; Kt=0,055 Nm/A; 
Kv=0,0567 Vs/rad. ; O diagrama de blocos no SIMULINK é desenhado na forma:
Programação no Matlab para gerar a função de transferência do motor com entrada Va e saída W, é a seguinte:
clc
clear all
s1 = tf(1,[0.00258 2.903]);
s2 = tf(0.055,1);
s3 = tf(1,[2.62*(10^-5) (10^-5)]);
s4 = tf(0.0567,1);
G = series(s1,s2);
H = series(G,s3);
sys = feedback(H,s4,-1)
sys =
 
 0.055
 -------------------------------------
 6.76e-08 s^2 + 7.608e-05 s + 0.003148
 
>> pole(sys)
ans = 1.0e+03 *
 -1.0826
 -0.0430
Programação no Matlab para gerar a função de transferência do motor com entrada TL e saída W, é a seguinte:
s1 = tf(1,[2.62*(10^-5) (10^-5)]);
s2 = tf(0.055*0.0567,[0.00258 2.903]);
sys = feedback(s1,s2,-1)
sys =
 
 0.00258 s + 2.903
 -------------------------------------
 6.76e-08 s^2 + 7.608e-05 s + 0.003148
 
Continuous-time transfer function.
>> pole(sys)
ans = 1.0e+03 *
 -1.0826
 -0.0430
>> pzmap(sys)
Programação no Matlab para gerar a função de transferência do motor com entrada Va e saída Ia, é a seguinte:
s1 = tf(1,[0.00258 2.903]);
s2 = tf(0.055*0.0567,[2.62*(10^-5) (10^-5)]);
sys = feedback(s1,s2,-1)
sys =
 
 2.62e-05 s + 1e-05
 -------------------------------------
 6.76e-08 s^2 + 7.608e-05 s + 0.003148
 
Continuous-time transfer function.
>> pole(sys)
ans = 1.0e+03 *
 -1.0826
 -0.0430
>> pzmap(sys)
Programação no Matlab para gerar a função de transferência do motor com entrada TL e saída W, é a seguinte:
s1 = tf(1,[2.62*(10^-5) (10^-5)]);
s2 = tf(0.0567,[0.00258 2.903]);
sys = feedback(s1,s2,-1)
sys =
 
 0.00258 s + 2.903
 ------------------------------------
 6.76e-08 s^2 + 7.608e-05 s + 0.05673
>> pole(sys)
ans = 1.0e+02 *
 -5.6279 + 7.2285i
 -5.6279 - 7.2285i
>> pzmap(sys)