Pratica C1 - Sistemas Com Polos Reais
Pratica C1 - Sistemas Com Polos Reais
Pratica C1 - Sistemas Com Polos Reais
Escola de Engenharia
Laboratório de Automação e Controle - EMA 105
Prof. Willian Moreira Duarte
Modelagem do sistema
t t
V (t ) qi (t )dt qo (t )dt
0 0
Variáveis:
V(t) = volume de fluido no tanque
qi(t) = vazão de entrada
qo(t) = vazão de saída
Para obtermos a função de transferência, vamos considerar como entrada da função a vazão de
entrada qi(t) e como saída a vazão de saída qo(t).
Precisamos então obter uma equação diferencial que apresente como variáveis apenas a entrada e a
saída. Utilizaremos algumas relações conhecidas como:
V (t ) Ah(t )
em que h(t) é a altura do nível do fluido e A é área do tanque (constante). Assim:
t t
Ah(t ) qi (t )dt qo (t )dt
0 0
Porém, ainda temos na equação a variável h(t) que não é a entrada nem a saída do sistema.
Utilizamos então a relação
h(t )
qo (t )
R
representando que a vazão de saída é proporcional à altura do nível do tanque com uma constante
de proporcionalidade R. Essa constante depende de fatores como densidade e viscosidade do fluido,
perda de carga na tubulação e válvula etc. Após a substituição, teremos:
t t
ARqo (t ) qi (t )dt qo (t )dt
0 0
Agora temos na equação somente constantes e variáveis de entrada e saída, então podemos obter a
função de transferência. Aplicando a transformada de Laplace:
Qi ( s) Qo ( s )
ARQo ( s)
s s
Qo ( s ) 1
Qi ( s ) ARs 1
(a) Definir no Matlab a função de transferência que relaciona Qo(s) e Qi(s) para R=0,9 m/(m3/s) e
A=4m2.
Comandos:
R=0.9
A=4
num=[1]
den=[A*R 1]
G1=tf(num,den)
Resposta do Matlab:
R =
0.9000
A =
num =
den =
3.6000 1.0000
G1 =
1
---------
3.6 s + 1
(b) Traçar qo(t) quando qi(t) sofre uma variação degrau unitário, para as seguintes situações:
Comandos:
clear all
close all
R=0.9
A1=4
A2=2
A3=1
num=[1]
den1=[A1*R 1]
den2=[A2*R 1]
den3=[A3*R 1]
G1=tf(num,den1)
G2=tf(num,den2)
G3=tf(num,den3)
step(G1)
hold on
step(G2)
hold on
step(G3)
hold on
Gráfico:
Comentários:
O tempo de subida é tanto mais rápido quanto menor é a área da seção do tanque, estabilizando o
sinal mais rapidamente.
G A (s ) GB (s ) GC (s )
Utilizando o Matlab:
(a) Determinar a função de transferência GA(s) que relaciona Q2(s) e Q1(s);
Comandos:
R=0.9
A=4
num=[1]
den=[A*R 1]
GA=tf(num,den)
Resposta do Matlab:
R =
0.9000
AA =
AB =
AC =
num =
denA =
3.6000 1.0000
denB =
2.7000 1.0000
denC =
0.9000 1.0000
GA =
1
---------
3.6 s + 1
GB =
1
---------
2.7 s + 1
GC =
1
---------
0.9 s + 1
Q1 ( s ) Q2 ( s ) Q3 ( s ) 1 Q4 ( s )
1 1
AA RA s 1 AB RB s 1 AC RC s 1
G A (s ) GB (s ) GC (s )
G AGB (s )
Comandos:
clear all
close all
R=0.9
AA=4
AB=3
AC=1
num=[1]
denA=[AA*R 1]
denB=[AB*R 1]
denC=[AC*R 1]
GA=tf(num,denA)
GB=tf(num,denB)
GC=tf(num,denC)
GAGB=GA*GB
Resposta do Matlab:
R =
0.9000
AA =
AB =
AC =
num =
denA =
3.6000 1.0000
denB =
2.7000 1.0000
denC =
0.9000 1.0000
GA =
1
---------
3.6 s + 1
GB =
1
---------
2.7 s + 1
GC =
1
---------
0.9 s + 1
GAGB =
1
--------------------
9.72 s^2 + 6.3 s + 1
Comandos:
clear all
close all
R=0.9
AA=4
AB=3
AC=1
num=[1]
denA=[AA*R 1]
denB=[AB*R 1]
denC=[AC*R 1]
GA=tf(num,denA)
GB=tf(num,denB)
GC=tf(num,denC)
GAGBGC=GA*GB*GC
Resposta do Matlab:
R =
0.9000
AA =
AB =
AC =
num =
denA =
3.6000 1.0000
denB =
2.7000 1.0000
denC =
0.9000 1.0000
GA =
1
---------
3.6 s + 1
Continuous-time transfer function.
GB =
1
---------
2.7 s + 1
GC =
1
---------
0.9 s + 1
GAGBGC =
1
---------------------------------
8.748 s^3 + 15.39 s^2 + 7.2 s + 1
Q1 ( s ) Q2 ( s ) Q3 ( s ) Q4 ( s ) H C (s )
1 1 1
R
AA RA s 1 AB RB s 1 AC RC s 1
Comandos:
HC=GAGBGC*R
Resposta do Matlab:
HC =
0.9
---------------------------------
8.748 s^3 + 15.39 s^2 + 7.2 s + 1
(e) Utilizando o Matlab, considerar que a entrada q1(t) é um degrau unitário, traçar em um
único gráfico as respostas q2(t), q3(t), q4(t) e h3(t).
Comandos:
GAGBGC=GA*GB*GC
HC=GAGBGC*R
step(GA)
hold on
GAB=GA*GB
step(GAB)
hold on
GABC=GA*GB*GC
step(GABC)
hold on
step(HC)
hold on
ylim([0,1.05])
legend()
Gráfico: