Sistemas de Controle
Sistemas de Controle
Sistemas de Controle
Primeira Avaliao de
Sistemas de Controle
SERRA
2015/1
1
Sistemas de Controle
2
RESUMO
3
Sumrio
RESUMO................................................................................................................... 3
2. Questes............................................................................................................ 7
4
2.11 Projete um controlador pelo mtodo direto de Ragazzini
assumindo os mesmos plos de malha fechada do caso
contnuo. Use um perodo de amostragem aproximadamente
igual a 1/10 do perodo de oscilao da resposta de malha
fechada do caso contnuo ............................................................................ 44
REFERNCIAS ....................................................................................................... 54
5
1. O SISTEMA DE CONTROLE
( )
( )( )
( )
( ) ( )
( )
( )
( )
6
2. QUESTES
7
possvel determinar o zero do controlador assumido , cancelando
deste modo, o polo da planta.
( ) ( )
( )
( ) ( )
( )
( )
( )
( )( ) ( ) ( )
( )
( )( )( ) ( )
8
Com os parametros do controlador definidos:
( )
( )
( )
9
2.2 Verifique por simulao se o desempenho da malha fechada
do sistema contnuo est de acordo com o estabelecido
10
Para evidenciar ainda mais os resultados obtidos foi utilizado o Simulink. A
Figura 2.3 mostra a resposta do sistema em malha fechada ao degrau. Na
Figura 2.4 evidencia um sobressinal em torno de 8% como previamente
estabelecido.
11
Figura 2.5 Tempo de acomodao
12
2.3 Escolha um perodo de amostragem pela regra da margem de
fase. Apresente o diagrama de Bode para justificar a escolha.
13
Pela regra da regra da margem de fase:
( )
14
2.4 Converta o controlador contnuo em discreto. Use a tcnica
bilinear com prewarping e compare as respostas em frequncia
do controlador continuo e discreto
Por este motivo, pode ser utilizado o prewarping para compensar esta
distoro em uma frequncia determinada, usualmente a frequncia de corte. A
compensao da distoro em frequncia dada pela relao abaixo:
( )
15
Utilizando as funes c2dOptions e c2d do Matlab pode se obter o
seguinte controlador Bilinear com Prewarping:
16
2.5 Verifique por simulao o desempenho da malha fechada do
sistema amostrado em relao ao contnuo. Verifique se est de
acordo com o comportamento determinado pela regra da
margem de fase
( )
( )
( ) ( )
( )
( ) ( )
( )
17
Figura 2.10 Resposta ao degrau de H(s) e H(z)
18
A relao da margem de fase com o coeficiente de amortecimento
descrita pela seguinte equao:
eps = 0 : 0.0001: 1;
fi=atan(2./(((4+eps.^-4).^0.5)-2).^0.5);
plot(fi,eps);
19
Com o coeficiente de amortecimento obtido pelo grfico
possvel calcular o sobressinal do sistema com o controlador bilinear com
prewarrpping:
20
Figura 2.14 Sobressinal do sistema com o cotrolador bilinear com prewarpping
Como pode ser visto na Figura 2.14 o sobressinal foi em torno de 16,5%,
prximo do previamente calculado de 16,77% .
21
2.6 Encontre o equivalente ZOH do processo P(s) usando o
mtodo de fraes parciais. Mostre por simulao que a P(z)
encontrada a verso discreta de P(s)
( )
( ) ( ) { } ( ) { }
( )
22
Aplicando a tranformada Z a um perdo de amostragem de T = 0,15s
( ) ( )
{ } ( ) ( ) ( )
( )
{ }
( )
Logo:
( ) ( )( )
( )
( ) ( )
( )
( )
23
Figura 2.15 Resposta ao degrau de Pzoh encontrado via fraes parciais
Figura 2.16 Resposta ao degrau de Pzoh encontrado via fraes parciais (Zoom)
Pode ser observado nas Figuras 2.15 e 2.16 que o Pzoh encontrado pelo
mtodo das fraes parcias representa a verso discreta da planta a um
perodo de amostragem de 0,15 segundos .
24
2.7 Determine a funo de transferncia H(z). Para que valores de
k p a H(z) permanece estvel. Use o critrio de Jury. Teste as
restries encontradas para k p
( )
( )
( )
( )( )
( )
( ) ( ) [ ][ ]
( )( )
( ) ( )
( )( )
25
( )( )
( ) ( )
( )( )( )
( ) ( )
( ) ( )
( )
( ) ( )
( )
( )
( )
( ) ( )
26
Utilizando o critrio de Jury conforme Castrucci pginas 400 e 401, faz
necessrio atender quatro condies para uma equao caracterstica de 3
ordem:
( )
1) | |
2) ( )>0
3) ( )
4) | | | |
Onde
| |
1 Codio:
2 Codio:
3 Codio:
27
Obtidos os polinmios b2 e b0 em funo de kp foi plotado o grfico |b2|
e |b0| (Figura 2.17) para anlise do quarta condio onde |b2|>|b0|:
28
Analisando as Figuras 2.19 e 2.20 podemos concluir que pela 4
condio -12,7 < kp < 3,9.
29
Intercedendo as quatros condies temos que os valores de kp onde o
sistema estvel:
30
Figura 2.23 kp = 0,5
1
Figura 2.24 kp = 1,5
Figura 2.25 kp = 2
31
Figura 2.26 kp = 3
32
Figura 2.29 kp = 4,1
33
2.8 Encontre a y(k) a partir de Y(z)=H(z)*R(z) usando o mtodo
de fraes parciais. Mostre por simulao que a expresso
obtida para y(k) est correta.
Considerando:
( ) ( ) ( )
( )
( )
( )
( )
( )
[ ]
( )
34
Utilizando o comando residue do Matlab:
Portanto:
( )
( ) ( )
( ) ( )
( )
( ) ( )
( ) ( ) ( )
( )( )
( )( )
35
A Figura 2.30 mostra y(k) para 1 < k < 50 . O grfico se asselha
bastante a Figura 2.31 que a resposta de H(z) ao degrau unitrio.
36
2.9 Projete controladores deadbeat e deadbeat aumentado, ambos
para o equivalente ZOH do processo P(s). Use um perodo de
amostragem aproximadamente igual a 1/5 do tempo de
acomodao 2% da resposta de malha fechada do caso contnuo
Como a resposta em malha fechada do caso contnuo foi projetado para ter
um tempo de acomodao de 2 segundos foi calculado o novo perodo de
amostragem:
( )
( )
37
( )
( )
( )
( )
Onde ( ):
( )
E P( ):
( )
38
Logo a funo de controlador DeadBeat dada por:
( )
( )
( )
( )
( )
39
( )
( )
( )
( )
40
2.10 Verifique por simulao o desempenho da malha fechada do
sistema amostrado para os dois controladores deadbeat.
Compare os sinais de controle. Verifique se esto de acordo
com o esperado
41
Figura 2.33 Primeiro Passo DeadBeat
A Figura 2.33 mostra que o controlador Deadbeat deu o primeiro passo foi a
0,45. Enquanto o DeadBeatAumentado conforme a Figura 2.34 foi 25% do
DeadBeat (0,1125) .
42
Figura 2.35 Sinal de controle DeadBeat e do DeadBeat Aumentado
43
2.11 Projete um controlador pelo mtodo direto de Ragazzini
assumindo os mesmos plos de malha fechada do caso
contnuo. Use um perodo de amostragem aproximadamente
igual a 1/10 do perodo de oscilao da resposta de malha
fechada do caso contnuo
( ) ( )
44
Para se obter o controdor Ragazzini necessrio obedecer as seguintes
regras para escolha de H(z) .
I) Causalidade
H(1) = 1
III) Estabilidade
45
Ser utilizado os plos dominantes do sistema contnuo, com isso H(z) ser
de grau relativo 1 (2 plos e um zero), o mesmo grau relativo da planta (3 plos
e 2 zeros), o que atende o critrio de causalidade. Utilizando a definio de z
temos os seguintes plos de H(z):
( )
( )
( )
( )
( )
Logo
( )
( )
46
Tendo P(z) e H(z) possvel obter o controlador Ragazzini definido pela
seguinte expresso:
( )
( )
( ) ( )
47
2.12 Verifique por simulao o desempenho da malha fechada do
sistema amostrado usando o controlador obtido pelo mtodo
direto de Ragazzini. Compare com as simulaes obtidas no
item 5 e comente
Figura 2.34 Resposta ao degrau do processo com os controladores contnuo, bilinear com
prewarpping e ragazzini
Figura 2.35 Sobressinal do processo com os controladores contnuo, bilinear com prewarpping
e ragazzini
48
2.13 Projete um controlador pelo mtodo de Dahlin assumindo
como constante de tempo o tempo de acomodao 2% da
resposta de malha fechada do caso contnuo. Use um perodo
de amostragem aproximadamente igual a 1/10 dessa
constante de tempo
( )
( )
( )
Onde:
Portanto:
( )
( )
( )
49
Com o novo perodo de amostragem e utilizando o comando c2d do
Matlab foi possvel obter o equivalente zoh da planta:
( )
( )
( ) ( )
50
Usando o comando zpk temos:
Simplificando Dz:
( )( )
( )
( )( )
( ) ( )
( ) ( )
( )( )
( )
( )
51
2.14 Verifique por simulao o desempenho da malha fechada do
sistema amostrado usando o controlador obtido pelo mtodo
de Dahlin. Verifique se est de acordo com o esperado
52
Figura 2.37 Resposta ao degrau do processo com os controladores contnuo e Dahlin
53
REFERNCIAS
Controle Automtico - Plinio de Lauro Castrucci & Anselmo Bittar & Roberto
Moura Sales
54