Compensadores digitais

Ver Todos os Exercícios

Teoria

Introdução

Quando estávamos estudando o controle analógico, nós estudamos um importante elemento dos sistemas de controle que são os compensadores.

Os compensadores, são usados para melhorar algum aspecto da resposta do nosso sistema. Eles podem ser usados para melhorar a resposta transitória, deixando o sistema mais estável e menos oscilatório. Outro uso muito comum dos compensadores é para reduzir o erro em regime permanente.

Lembra como era a função de transferência padrão dos compensadores? Não? Pois ela tá aqui embaixo:

G c s = K c s + a s + b

Os compensadores podem ser do tipo avanço de fase ou atraso de fase.

Os compensadores de avanço de fase possuem o zero maior que o polo. Esses compensadores são usados para melhorar a resposta transitória do sistema.

a > b

Os compensadores de atraso de fase, por sua vez, possuem os zeros muito menores que o polo. Outra característica desses compensadores é que tanto os polos quanto os zeros são valores bem próximos de 0 . Os compensadores de atraso de fase são usados principalmente para reduzir o erro em regime permanente do sistema.

a < b

Os compensadores digitais, por sua vez, vão manter as principais características dos compensadores analógicos. A função de transferência geral dos compensadores digitais será análoga à dos compensadores analógicos, apenas mudando o s por z .

G c z = K c z - a z - b

Como os polos e zeros dos compensadores digitais podem ser positivos, nós vamos usar a expressão geral com “menos” e não com “mais”.

Além disso, os compensadores digitais têm a mesma função no sistema que os compensadores analógicos: melhorar a resposta transitória (avanço), reduzir o erro em regime permanente (atraso). E os métodos para projetar os compensadores digitais são análogos aos métodos para o projeto analógico. Vamos ver melhor como fazer esse projeto?

Compensadores de avanço

Para projetar compensadores de avanço analógicos, nós precisamos de especificações da resposta transitória. Essas especificações geralmente são dadas na forma do overshoot máximo permitido (ou fator de amortecimento, ζ ) e do tempo de decaimento, T s .

O S = e - ζ π 1 - ζ 2

T s ( 2 % ) = 4 ζ ω n

Com ζ e ω n nós devemos calcular o polo dominante, e todo o projeto se dá em torno desse polo dominante. O polo dominante tem a seguinte fórmula:

s = - ω n ζ ± i ω n 1 - ζ 2

Para projetar um compensador digital, nós também devemos seguir todos esses passos, e após achar o polo dominante no domínio s , devemos transformar esse polo para o domínio z .

A demonstração para essa transformação é um pouco chata, então vou pular direto para o que interessa. Para achar o polo dominante no domínio z , devemos encontrar o módulo e o ângulo de fase do polo no domínio z . Isso pode ser calculado da seguinte fórmula:

z = e - T ω n ζ  

∠ z = T ω n 1 - ζ 2

Onde T é o período de amostragem do sistema.

Sendo assim, após isso teremos o polo dominante no domínio z em coordenadas polares:

z = z ∠ z

Por exemplo, se tivermos os polos dominantes no plano s :

s = - 3 ± 1 , 5 i

Sendo assim, teremos:

- ω n ζ = - 3

ω n 1 - ζ 2 = 0,5

Podemos calcular então onde esse polo vai se encontrar no plano z utilizando as fórmulas ali de cima. Vamos considerar o período de amostragem t = 0,1 :

z = e - 3 T = 0,74

∠ z = 1,5 T = 015   r a d

Ou seja:

z = 0,74 ∠ 0,15

Mas isso está em coordenadas polares, significo que o polo tem módulo 0,74 e faz um ângulo de 0,15   r a d . Graficamente, podemos representar isso da seguinte forma:

Precisamos passar isso para coordenadas retangulares, ou seja, precisamos encontrar σ e ω para conseguir escrever o polo da seguinte maneira:

z = σ + ω i

Podemos fazer isso usando trigonometria:

ω = z × cos ⁡ ∠ z

σ = z × sen ⁡ ∠ z

Para o nosso exemplo teremos:

ω = 0 , 7 3

σ = 0 , 11

Portanto:

z = 0,73 ± 0,11 i

Dica: essa conversão também pode ser feita usando a função “Rec(“ da sua calculadora.

Em seguida, o projeto do compensador de avanço de fase deve seguir normalmente.

Para nossa felicidade, as regras que usamos no projeto analógico ainda funcionam aqui. Portanto:

∠ G z C ( z ) = - 180 º

G z C ( z ) = 1

Sendo assim, podemos arbitrar um zero e encontrar um zero que garanta que o lugar das raízes passe pelo polo dominante através da regra do ângulo de fase.

Após encontrar o polo e o zero, podemos encontrar o ganho do compensador através da regra do módulo.

Se você não lembra muito bem como utilizar essas regras, seria legal você dar uma olhada na aula de compensadores de avanço. Lá tem tudo explicadinho.

Compensadores de atraso

Para projetar os compensadores de atraso, nós devemos pensar principalmente no erro em regime permanente.

Quando a gente fazia o projeto para o caso analógico, a gente supunha o ganho do compensador como 1 e usava o teorema do valor final para encontrar uma relação entre o zero e o polo:

lim s → 0 ⁡ s + a s + b G s = k p

e ∞ = 1 1 + k p

Onde e ∞ é o erro em regime permanente desejado e k p é a constante de erro.

Ao aplicar esse teorema encontramos uma relação entre o polo e o zero, por exemplo:

b = 10 a

O próximo passo seria supor um valor bem pequeno para o zero, e calcular o valor do polo usando a relação encontrada acima.

Vamos pensar como podemos adaptar isso tudo para o mundo digital agora.

Primeiramente vamos lembrar da definição da variável z :

z = e T s

Sendo assim:

s → 0

z → 1

Portanto, para usar o teorema do valor final no domínio z , ao invés de fazermos a variável tendendo a 0 , ela deverá tender a 1 .

Portanto, para uma entrada rampa usaremos:

e ∞ = 1 1 + k p

k p = lim z → 1 ⁡ G O L z

E, para uma entrada rampa:

e ∞ = 1 k v

k v = 1 T lim z → 1 ⁡ z - 1 G O L z

Onde G O L ( z ) é a função de transferência do compensador e do processo:

G O L z = G C z G p ( z )

Sendo assim, após aplicar esse teorema teremos algo do tipo:

1 - a 1 - b = k

Onde k é uma constante qualquer.

Por exemplo, se temos um processo cuja função de transferência é:

G z = z z - 0,5 z - 0,2

E queremos um compensador que faça com que o erro em estado permanente seja 0,05 , deveremos primeiramenrte encontrar o k p que gera esse erro:

e ∞ = 1 1 + k p

k p = 19

Com k p poderemos encontrar a relação que queremos:

k p = lim z → 1 ⁡ G O L z

k p = lim z → 1 ⁡ z - a z - b z z - 0,5 z - 0,2

Isso vai nos dar:

2,5 1 - a 1 - b = 19

1 - a 1 - b = 7,6

E com isso conseguimos continuar o projeto do compensador de atraso normalmente.

Agora só precisamos arbitrar um valor qualquer para o zero e calcular o polo usando a relação que conseguimos através do teorema do valor final acima.

Quando estávamos fazendo o projeto analógico, nós usávamos um valor para o polo bem próximo de 0 . Aqui, como fizemos a mudança para o sinal digital, nós vamos querer que o zero seja o mais próximo de 1 possível. Geralmente 0,99 é um bom palpite para o zero.

Após isso, basta calcular o ganho do compensador através da regra do módulo:

G z C ( z ) = 1

Compensadores de avanço

Compensadores de atraso

Exercícios Resolvidos

Exercício Resolvido #1

Para entradas em degrau e rampa obtenha o erro em regime permanente para a função de transferência em malha aberta dada abaixo:

G z = 10 T z - 1 - 1 + z - 1 z - e - T

Passo 1

Vamos fazer primeira para a entrada degrau.

As fórmulas para o erro da entrada degrau são:

k p = lim z → 1 ⁡ G O L z

e ∞ = 1 1 + k p

Vamos calcular a constante do erro k p então:

k p = lim z → 1 ⁡ 10 T z - 1 - 1 + z - 1 z - e - T = 10 ∞ - 1 + 0 = ∞

Calculando o erro:

e ∞ = 1 1 + k p = 1 1 + ∞ = 0

Passo 2

Vamos agora calcular o erroa para a entrada rampa.

As fórmulas de erro para a entrada rampa são:

k v = 1 T lim z → 1 ⁡ z - 1 G O L z

e ∞ = 1 k v

Calculando primeiramente o k v :

k v = 1 T lim z → 1 ⁡ 10 z - 1 T z - 1 - 1 + z - 1 z - e - T = 10 T lim z → 1 ⁡ T - z - 1 + z - 1 2 z - e - T

k v = 10 T × T = 10

O erro será portanto

e ∞ = 1 k v = 1 10 = 0,1

Resposta

Entrada rampa:

e ∞ = 0

Entrada degrau:

e ∞ = 0,1

Exercício Resolvido #2

Para entradas em degrau e rampa obtenha o erro em regime permanente para o sistema mostrado abaixo:

G 1 s = 20 s + 3 s + 4 s + 5

Passo 1

O primeiro passo será encontrar a função de transferência do sistema no plano z .

No plano s temos que:

G O L s = 1 - e - T s s 20 s + 3 s + 4 s + 5

Precisamos aplicar a transformada Z na função de transferência, sendo assim:

G O L z = Z 1 - e - T s s 20 s + 3 s + 4 s + 5

Manipulando um pouco:

G O L z = Z 1 - e - T s Z 20 s + 3 s s + 4 s + 5

Como pela definição temos que:

z = e T s

Podemos fazer:

Z 1 - e - T s = 1 - z - 1 = z - 1 z

Portanto:

G O L z = z - 1 z Z 20 s + 3 s s + 4 s + 5

Passo 2

Vamos agora expandir 20 s + 3 s s + 4 s + 5 em frações parciais para encontrar a transformada Z desse termo:

20 s + 3 s s + 4 s + 5 = 20 0,15 s + 0,25 s + 4 - 0,4 s + 5

Consultando uma tabela de transformada Z , teremos:

Z 20 0,15 s + 0,25 s + 4 - 0,4 s + 5 = 20 0,15 z z - 1 + 0,25 z z - e - 4 T - 0,4 z z - e - 5 T

Logo:

G O L z = 20 z - 1 z 0,15 z z - 1 + 0,25 z z - e - 4 T - 0,4 z - e - 5 T

G O L z = 3 + 5 z - 1 z - e - 4 T - 8 z - 1 z - e - 5 T

Passo 3

Vamos fazer primeira para a entrada degrau.

As fórmulas para o erro da entrada degrau são:

k p = lim z → 1 ⁡ G O L z

e ∞ = 1 1 + k p

Vamos calcular a constante do erro k p então:

k p = lim z → 1 ⁡ 3 + 5 z - 1 z - e - 4 T - 8 z - 1 z - e - 5 T = 3

Calculando o erro:

e ∞ = 1 1 + k p = 1 1 + 3 = 0,25

Passo 4

Vamos agora calcular o erroa para a entrada rampa.

As fórmulas de erro para a entrada rampa são:

k v = 1 T lim z → 1 ⁡ z - 1 G O L z

e ∞ = 1 k v

Calculando primeiramente o k v :

k v = 1 T lim z → 1 ⁡ z - 1 3 + 5 z - 1 z - e - 4 T - 8 z - 1 z - e - 5 T

k v = 1 T lim z → 1 ⁡ 3 ( z - 1 ) + 5 z - 1 2 z - e - 4 T - 8 z - 1 2 z - e - 5 T = 0

O erro será portanto

e ∞ = 1 k v = 1 0 = ∞

Resposta

Entrada rampa:

e ∞ = 4

Entrada degrau:

e ∞ = ∞

Exercício Resolvido #3

Dado a função de transferência em malha aberta abaixo:

G O L z = 5 z - 0,65 z - 0,3 z + 0,1

É desejado implementar um compensador de avanço de fase que reduza o erro em regime permanente para uma entrada degrau desse sistema em 20 vezes. Para que isso aconteça, qual deve ser a relação entre o polo e o zero do compensador.

Passo 1

Esse exercício quer que a gente determine uma boa relação entre o zero e o polo do compensador de avanço que será implementado para que ocorra a redução de erro em regime permanente desejada para uma entrada degrau.

Para fazer isso, vamos com calma. Primeiramente, vamos calcular o erro em regime permanente para uma entrada degrau desse sistema sem o compensador.

As fórmulas para o erro da entrada degrau são:

k p = lim z → 1 ⁡ G O L z

e ∞ = 1 1 + k p

Ou seja:

k p = lim z → 1 ⁡ 5 z - 0,65 z - 0,3 z + 0,1 = 5 × 0,35 0,7 × 1,1 = 2,27

e ∞ = 1 1 + 2,27 = 0,31

Passo 2

O objetivo é que a compensação reduza o erro em 20 vezes, portanto, o novo erro deverá ser:

e * ∞ = 0,31 20 = 0,065

E o novo k p * :

1 1 + k p * = 0,065

k p * = 14,4

Passo 3

A função de transferência do sistema com compensação será:

G O L * z = z - a z - b 5 z - 0,65 z - 0,3 z + 0,1

Portanto, podemos achar a relação entre os polos, calculando o novo k p * :

k p * = lim z → 1 ⁡ z - a z - b 5 z - 0,65 z - 0,3 z + 0,1 = 14,4

Portanto:

1 - a 1 - b × 2,27 = 14,4

1 - a 1 - b = 6,37

Resposta

1 - a 1 - b = 6,37

Exercício Resolvido #4

Dados os polos dominantes no plano s :

s = - 2 ± 1,5 i

Encontre o polo dominante equivalente no plano z . Assuma T = 0,1

Passo 1

Bora lá.

Nesse exercício nós temo um par de polos no plano s e precisamos definir quais serão os polos equivalentes no plano z . Para isso, encontraremos primeiramente o polo z em coordenadas polares e em seguida converteremos esse polo para coordenadas polares. Sendo assim, primeiramente o polo dominante em z terá a cara:

z = z ∠ z

Sendo que:

z = e - T ω n ζ  

∠ z = T ω n 1 - ζ 2  

Passo 2

A fórmula geral de um polo dominante no plano s é:

s = - ω n ζ ± i ω n 1 - ζ 2

E o nosso polo é:

s = - 2 ± 1,5 i

Comparando os dois teremos:

ω n ζ = 2

ω n 1 - ζ 2 = 1,5

Portanto, podemos calcular o módulo e o ângulo do polo no plano z :

z = e - 2 × 0,1 = e - 0,2 = 0,82

∠ z = 1,5 × 0,1 = 0,15   r a d

Portanto, o polo no plano z em coordenadas angulares será:

z = 0,82 ∠ 0,15

Passo 3

Para encontrar esse polo em coordenadas retangulares, vamos usar a calculadora. Não precisa ser nenhuma calculadora superpoderosa para fazer essa conversão, aquela basiquinha que todo aluno de engenharia tem consegue fazer isso perfeitamente.

Para fazer a conversão você deve apertar o botão “shift” e em seguida o botão “Pol(” ativando assim, a função “Rec(“ da calculadora. Em seguida, digite o valor do módulo que calculamos, aperte o botão da vírgula, digite o valor do ângulo, e feche o parênteses. Certifique-se também de que sua calculadora está no modo “radianos” e não no modo “graus”.

Após fazer isso, aperte igual e aparecerá o valor da parte real no visor, para o nosso caos:

σ = 0 , 81

Para encontrar a parte imaginária, pressione a tecla “ALPHA” e em seguida a tecla “tan”. Aparecerá um “F” no visor da calculadora, aí basta apertar igual e aparecerá o valor da parte imaginária no visor. Para o nosso exemplo, esse valor é:

ω = 0 , 13

Portanto, os nossos polos dominantes no planos z serão:

z = 0 , 81 ± 0,13 i

Resposta

z = 0,81 ± 0,13 i

Exercício Resolvido #5

Projete um compensador de atraso que faça com que o erro em regime permanente do sistema abaixo seja de 0,02 para uma entrada degrau.

G z = 5 ( z - 0,5 ) z ( z - 0,3 )

Passo 1

Vamos primeiro descobrir qual o k p que faz com que o sistema produza o erro desejado. Para isso, vamos usar a fórmula que relaciona o erro com o k p :

e ∞ = 1 1 + k p

Para o erro ser 0,05 :

0,05 = 1 1 + k p

k p = 19

Passo 2

Agora, com o k p precisamos encontrar a relação entre o polo e zero do compensador usando a fórmula do próprio k p :

k p = lim z → 1 ⁡ G O L z

Onde G O L ( z ) representa a função de transferência em malha aberta que é o produto da função de transferência do processo com a função de transferência do compensador que estamos projetando. Ou seja:

k p = lim z → 1 ⁡ z - a z - b 5 ( z - 0,5 ) z ( z - 0,3 )

Obs: Nessa etapa, podemos considerar o ganho do compensador como sendo 1 .

Como o k p deve ser 49 , teremos:

lim z → 1 ⁡ z - a z - b 5 ( z - 0,5 ) z ( z - 0,3 ) = 19

Fazendo o limite:

1 - a 1 - b 5 × 0,5 0,7 = 19

1 - a 1 - b = 5,32

Passo 3

Agora que já temos a relação entre o polo e o zero do compensador, podemos arbitrar um dos dois e encontrar o outro.

Para isso, geralmente arbitramos o polo e encontramos o zero. O melhor a se fazer é arbitrar um valor para o polo bem próximo de 1 para que o zero não fique muito longe do polo.

Quanto mais próximos o zero e o polo estão, menor é o impacto do compensador na resposta trtansitória do sistema. Dessa forma diminui-se o erro em regime permanente sem alterar ζ e o tempo de acomodação.

Vamos então arbitrar b = 0,99 . Assim teremos:

1 - a 1 - 0,99 = 5,32

1 - a 0,01 = 5,32

a = 0,95

Dessa forma, a expressão para o compensador será:

G C ( z ) = 1 - 0,95 1 - 0,99

Como o polo e o zero estão bem próximos, não há necessidade de se calcular um ganho para o compensador uma vez que este será bem próximo de 1 .

Caso esses valores não fossem tão próximos, ou o rigor do projeto fosse maior, seria necessário encontrar os polos dominantes do sistema sem compensação e depois usar a lei dos módulos para encontrar o ganho do compensador.

Resposta

G C ( z ) = 1 - 0,95 1 - 0,99

Exercício Resolvido #6

Projete um compensador de atraso que faça com que o erro em regime permanente do sistema abaixo seja de 0,1 para uma entrada rampa.

G z = ( z - 0,8 ) ( z - 1 ) ( z - 0 , 2 )

Passo 1

Vamos primeiro descobrir qual o k v que faz com que o sistema produza o erro desejado. Para isso, vamos usar a fórmula que relaciona o erro com o k v :

e ∞ = 1 k v

Para o erro ser 0 , 1 :

0,1 = 1 k v

k v = 1 0

Passo 2

Agora, com o k v precisamos encontrar a relação entre o polo e zero do compensador usando a fórmula do próprio k v :

k v = lim z → 1 ⁡ ( z - 1 ) G O L z

Onde G O L ( z ) representa a função de transferência em malha aberta que é o produto da função de transferência do processo com a função de transferência do compensador que estamos projetando. Ou seja:

k p = lim z → 1 ⁡ [ z - 1 ] z - a z - b ( z - 0,8 ) ( z - 1 ) ( z - 0,2 )

Obs: Nessa etapa, podemos considerar o ganho do compensador como sendo 1 .

Como o k v deve ser 10 , teremos:

lim z → 1 ⁡ [ z - 1 ] z - a z - b ( z - 0,8 ) ( z - 1 ) ( z - 0,2 ) = 1 0

Fazendo o limite:

1 - a 1 - b 0 , 8 0,2 = 1 0

1 - a 1 - b = 2,5

Passo 3

Agora que já temos a relação entre o polo e o zero do compensador, podemos arbitrar um dos dois e encontrar o outro.

Para isso, geralmente arbitramos o polo e encontramos o zero. O melhor a se fazer é arbitrar um valor para o polo bem próximo de 1 para que o zero não fique muito longe do polo.

Quanto mais próximos o zero e o polo estão, menor é o impacto do compensador na resposta trtansitória do sistema. Dessa forma diminui-se o erro em regime permanente sem alterar ζ e o tempo de acomodação.

Vamos então arbitrar b = 0,99 . Assim teremos:

1 - a 1 - 0,99 = 2,5

1 - a 0,01 = 2,5

a = 0,9 7 5

Dessa forma, a expressão para o compensador será:

G C ( z ) = 1 - 0,9 7 5 1 - 0,99

Como o polo e o zero estão bem próximos, não há necessidade de se calcular um ganho para o compensador uma vez que este será bem próximo de 1 .

Caso esses valores não fossem tão próximos, ou o rigor do projeto fosse maior, seria necessário encontrar os polos dominantes do sistema sem compensação e depois usar a lei dos módulos para encontrar o ganho do compensador.

Resposta

G C ( z ) = 1 - 0,9 7 5 1 - 0,99

Exercício Resolvido #7

Dado o sistema realimentado com T = 0,2 s , projete um compensador digital que atenda às especificações:

G z = 0,16 z - 0,67

T s 5 % = 2 s

ζ = 0,5

Passo 1

O primeiro passo para fazer o projeto no plano z é encontrar os polos dominantes no plano s e depois convertê-lo ao plano z .

O polo dominante no plano s terá forma:

s = - ω n ζ ± i ω n 1 - ζ 2

O valor de ζ nós já temos que foi dado como especificação de projeto, ζ = 0,5 . Para achar ω n vamos usar a especificação do tempo de assentamento. Para uma faixa de 5 % a fórmula para o tempo de decaimento é a seguinte:

T s ( 5 % ) = 3 ω n ζ

Substituindo os valores:

2 = 3 0,5 × ω n

ω n = 3 0,5 × 2 = 3

Os polos dominantes no plano s serão portanto:

s = - 1,5 ± 2,6 i

Passo 2

Agora precisamos determinar onde estarão esses polos no plano z atravé das seguintes fórmulas:

z = e - T ω n ζ  

∠ z = T ω n 1 - ζ 2

Portanto, podemos calcular o módulo e o ângulo do polo no plano z :

z = e - 1,5 × 0,2 = e - 0,3 = 0,74

∠ z = 2,6 × 0,2 = 0,52   r a d

Portanto, o polo no plano z em coordenadas angulares será:

z = 0,74 ∠ 0,52

Fazendo a conversão de coordenadas polares para retangulares teremos:

z = 0,64 ± 0,37 i

Passo 3

Vamos representar os polos dominantes e o polo da função de transferência do processo no diagrama abaixo:

O ponto preto indica o polo dominante em 0,64 ± 0,37 i e o “xis” vermelho representa o polo da função de transferência do processo. Agora precisamos encontrar o polo e o zero do compensador. Para isso, podemos simplesmente arbitrar um polo qualquer e então encontrar onde vai estar o zero correspondente.

Como o nosso compensador será de avanço sabemos que o polo tem que estar mais à direita que o zero. Sendo assim, vamos tentar arbitrar o nosso polo um pouco mais à direita do nosso diagrama.

Eu vou escolher 0,85 simplesmente porque eu quero, mas você poderia ter escolhido qualquer valor para o polo. Sendo assim, vamos representar nosso diagrama agora:

O “xis” azul representa o polo do compensador e o círculo azul o zero. O zero do compensador foi colocado em um lugar arbitrário, sendo assim, precisamos determiner qual é o ponto exato do zero considerando que o polo foi colocado em 0,85 .

Passo 4

Para determinar o ponto exato do zero, usaremos a regra dos ângulos que diz que a soma dos ângulos dos polos menos a soma dos ângulos dos zeros deve ser igual a um múltiplo ímpar de 180 ° .

Sendo assim, teremos:

θ 1 + θ 2 - θ 3 = 180 ° ( 2 r + 1 )

Vamos calcular os ângulos θ 1 e θ 2 :

θ 1 = 180 - tan - 1 ⁡ 0,37 0,85 - 0,64   = 118,6 °

θ 2 = 180 - tan - 1 ⁡ 0,37 0,67 - 0,64 = 94,6 °

Sendo assim, podemos determinar θ 3 :

118,6 + 94,6 - θ 3 = 180

θ 3 = 33,2 °

Podemos agora determinar o ponto exato do zero do compensar, já que encontramos o ângulo que ele faz com o polo dominante:

tan ⁡ θ 3 = 0,37 0,64 - b

tan ⁡ 33,2 ° = 0,37 0,64 - b

Isso vais nos dar:

b = 0,07

Portanto, nosso compensador por enquanto está assim:

G c z = K c z - 0,07 z - 0,85  

Falta apenas calcular o ganho.

Passo 5

Para calcular o ganho, usaremos o critério do módulo. Para isso, precisaremos calcular as distâncias entre os polos dominantes e os polos e zeros do sistema:

Precisamos calcular h 1 , h 2 e h 3 :

h 1 = 0,37 2 + 0,85 - 0,64 2 = 0,43

h 2 = 0,37 2 + 0,67 - 0,64 2 = 0,37

h 3 = 0,37 2 + 0,64 - 0 , 07 2 = 0,70

A regra do módulo diz que:

K c × K × h 3 h 1 × h 2 = 1

Sendo K c o ganho do compensador e K o ganho do processo. Sendo assim:

K c = h 1 × h 2 K × h 3

K c = 0,43 × 0,37 0,16 × 0,70 = 1,42

Sendo assim, calculamos finalmente a função de transferência do nosso compensador de avanço:

G c z = 1,42 z - 0,07 z - 0,85

Resposta

G c z = 1,42 z - 0,07 z - 0,85

Exercício Resolvido #8

Considere o sistema de controle mostrado abaixo. No plano z , projete um compensador digital que apresente polos dominantes em malha fechada com uma constante de amortecimento ζ = 0,5 e tempo de acomodação para uma faixa de 2 % de 2 s . O período de amostragem utilizado será 0,2 s .

G z = 0,017 ( z + 0,88 ) ( z - 1 ) ( z - 0,67 )

Passo 1

O primeiro passo para fazer o projeto no plano z é encontrar os polos dominantes no plano s e depois convertê-lo ao plano z .

O polo dominante no plano s terá forma:

s = - ω n ζ ± i ω n 1 - ζ 2

O valor de ζ nós já temos que foi dado como especificação de projeto, ζ = 0,5 . Para achar ω n vamos usar a especificação do tempo de assentamento. Para uma faixa de 5 % a fórmula para o tempo de decaimento é a seguinte:

T s ( 2 % ) = 4 ω n ζ

Substituindo os valores:

2 = 4 0,5 × ω n

ω n = 4 0,5 × 2 = 4

Os polos dominantes no plano s serão portanto:

s = - 2 ± 3,5 i

Passo 2

Agora precisamos determinar onde estarão esses polos no plano z atravé das seguintes fórmulas:

z = e - T ω n ζ  

∠ z = T ω n 1 - ζ 2

Portanto, podemos calcular o módulo e o ângulo do polo no plano z :

z = e - 2 × 0,2 = e - 0,4 = 0,67

∠ z = 3,5 × 0,2 = 0,70   r a d

Portanto, o polo no plano z em coordenadas angulares será:

z = 0,67 ∠ 0,70

Fazendo a conversão de coordenadas polares para retangulares teremos:

z = 0,51 ± 0,43 i

Passo 3

Vamos representar os polos dominantes e o polo da função de transferência do processo no diagrama abaixo:

O ponto preto indica o polo dominante em 0,51 ± 0,43 i , os “xis” vermelho representam os polos e os círculos vermelhos representam os zeros da função de transferência do processo. Agora precisamos encontrar o polo e o zero do compensador. Para isso, podemos simplesmente arbitrar ou o polo ou o zero e depois encontrar o outro usando a regra dos ângulos.

Uma estratégia muito usada nessas situações é arbitrar um polo/zero de forma que ele cancele um polo/zero já existente, facilitando assim o projeto do compensador. Aqui por exemplo, podemos arbitrar o zero do compensador em 0,67 de forma que cancele com o polo do processo que já existe.

Sendo assim, nosso diagrama ficará:

Repare que o polo em 0,67 sumiu porque ele foi cancelado com o zero do compensador. O “xis” azul representa o polo do compensador que precisamos determiner em que ponto ficará.

Passo 4

Para determinar o ponto exato do polo, usaremos a regra dos ângulos que diz que a soma dos ângulos dos polos menos a soma dos ângulos dos zeros deve ser igual a um múltiplo ímpar de 180 ° .

Sendo assim, teremos:

θ 1 - θ 2 + θ 3 = 180 ° ( 2 r + 1 )

Vamos calcular os ângulos θ 1 e θ 3 :

θ 1 = tan - 1 ⁡ 0,43 0,88 + 0,51   = 17,9 °

θ 3 = 180 - tan - 1 ⁡ 0,43 1,0 - 0,51 = 138,7 °

Sendo assim, podemos determinar θ 2 :

- 17,9 + θ 2 + 138,7 = 180 °

θ 2 = 59,2 °

Podemos agora determinar o ponto exato do polo do compensar, já que encontramos o ângulo que ele faz com o polo dominante:

tan ⁡ θ 2 = 0,43 0,51 - b

tan ⁡ 59,2 ° = 0,43 0,51 - b

Isso vais nos dar:

b = 0,25

Portanto, nosso compensador por enquanto está assim:

G c z = K c z - 0,67 z - 0,25  

Falta apenas calcular o ganho.

Passo 5

Para calcular o ganho, usaremos o critério do módulo. Para isso, precisaremos calcular as distâncias entre os polos dominantes e os polos e zeros do sistema:

Precisamos calcular h 1 , h 2 e h 3 :

h 1 = 0,43 2 + 0,88 + 0,51 2 = 1,45

h 2 = 0,43 2 + 0,51 - 0,25 2 = 0,50

h 3 = 0,37 2 + 1,0 - 0,51 2 = 0,61

A regra do módulo diz que:

K c × K × h 1 h 2 × h 3 = 1

Sendo K c o ganho do compensador e K o ganho do processo. Sendo assim:

K c = h 2 × h 3 K × h 1

K c = 0,50 × 0,61 0,017 × 1,45 = 12,4

Sendo assim, calculamos finalmente a função de transferência do nosso compensador de avanço:

G c z = 12,4 z - 0,67 z - 0,25

Resposta

G c z = 1,42 z - 0,07 z - 0,85