Métodos de Ziegler-Nichols

Ver Todos os Exercícios

Teoria

Introdução

No capítulo de controladores nós aprendemos um pouco sobre esses equipamentos que são o grande cérebro das malhas de controle. Os controladores são os caras que tomam as decisões de como a malha deve reagir quando a variável controlada não está ajustada na especificação adequada.

Nós vimos que existem alguns tipos de controladores, cada um com sua própria função de transferência e com um propósito diferente para ser usado.

Um dos controles mais simples é o proporcional (P), ele corrige a variável controlada proporcionalmente ao erro. O controlador proporcional-integrativo (PI) corrige a variável proporcionalmente ao erro e ainda acrescenta uma parcela de integral do erro. O proporcional-integrativo-derivativo (PID) faz tudo isso e ainda tem um efeito derivativo que prevê a resposta do sistema. O PID por exemplo, tem função de transferência:

G C s = K c 1 + 1 τ i s + τ d s

Alguns livros preferem escrever da forma abaixo, mas dá na mesma:

G c s = K p + K i s + K d s

K p = K c

K i = k c τ i

K d = K c τ d

Na fórmula dos controladores existem parâmetros: o ganho do controlador ( K C ), a constante de tempo integral ( τ i ) e a constante de tempo derivativo ( τ d ). Esses parâmetros são fundamentais para o bom funcionamento do controlador e o sucesso do sistema de controle vai depender diretamente do quão bem ajustados eles estão.

Como o controlador é o centro de toda a malha de controle, nós precisamos dar uma atenção muito especial para os métodos de sintonização desses equipamentos. O ajuste fino dos parâmetros do controlador é o que vai definir não só as características da nossa resposta mas como também a estabilidade do sistema.

Os primeiros métodos de sintonia de controlador que veremos serão os métodos de Ziegler-Nichols. Os brabos propuseram dois métodos, um em malha aberta e um em malha fechada, para sintonizar controladores. Os métodos de Ziegler-Nichols também são chamados de métodos de decaimento de 1 / 4 porque já que eles foram desenvolvidos de modo que cada pico da resposta fosse 1 / 4 do pico anterior.

O método de Ziegler-Nichols de malha aberta

Esse método, que também é chamado de método da curva de reação, consiste em inserir uma entrada degrau unitário em um processo em malha aberta e observar a resposta desse sistema com o tempo.

Ziegler e Nichols perceberam que a maioria dos processos apresentam uma resposta parecida com a da imagem quando submetidos a uma entrada do tipo degrau unitário.

Sendo assim, eles elaboraram um método gráfico de forma que seja possível aproximar o sistema que tem uma resposta em “S” como a da figura acima a um sistema de primeira ordem mais tempo morto:

Y s X s = K τ s + 1 e - θ s  

O método consiste em definir o ganho, K , a constante de tempo, τ , e o tempo morto, θ , através do gráfico da resposta do sistema. Com esses três valores é possível calcular os parâmetros para o controlador.

Para fazer essa análise, é preciso traçar a linha tangente na curva no ponto onde a inclinação dessa linha é a maior possível (geralmente vai ser perto do meio).

Após traçar a tangente de maior inclinação o ponto em que ela cruza o valor inicial da resposta ( y ( 0 ) ) determina o tempo morto o tempo morto. A distância entre esse ponto e o ponto em que ela cruza a linha do valor da resposta em estado estacionário é τ . O ganho é a variação total da resposta dividido pela amplitude do degrau.

Os parâmetros para o controlador podem então ser calculados utilizando-se as correlações da tabela abaixo:

Esse método é útil quando não temos a função de transferência do processo, mas sabemos como sua resposta se comporta para uma entrada degrau unitário.

O método de Ziegler-Nichols de malha fechada

Esse método, também conhecido como método do ganho final, utiliza-se da equação característica de malha fechada de um processo para determinar os parâmetros do controlador.

Dado um processo qualquer G ( s ) , controlado por um controlador G C ( s ) , sua função de transferência em malha fechada será:

Y s X s = G C ( s ) G s 1 + G C ( s ) G s

O denominador dessa função de transferência é chamado de equação característica. Com essa equação conseguimos determinar, por exemplo, qual o maior ganho que o sistema pode ter antes de ficar instável, que é o ganho final.

O método de Ziegler-Nichols de malha fechada consiste em determinar o ganho e o período final de um processo e com esses valores sintonizar os parâmetros do controlador.

Para encontrar o ganho e período final, você deve considerar que está trabalhando com um controlador proporcional (mesmo que esteja sintonizando um PI ou um PID). Dessa forma, você deve determinar K c u e ω u que são o ganho e a frequência que tornam o seu sistema instável.

O período final pode ser encontrado através da frequência final:

T u = 2 π ω u

Existe uma série de métodos que você pode usar: substituição direta, critério de Routh, lugar das raízes. Fica a critério do cliente. Caso você não se lembre como aplicar esses métodos, pode ficar tranquilo que tem uma aula para cada um deles aqui no Responde Aí, só voltar lá e dar uma revisada.

Com os valores de K c u e T u em mãos, é só jogar na tabelinha abaixo e ser feliz.

Uma dica pessoal minha é usar o método da substituição direta para encontrar os valores de K c u e T u . Eu particularmente considero esse método mais fácil e direto. Entretanto, todos os métodos encontram respostas muito semelhantes e se você se sentir mais confortável usando o critério de Routh ou lugar das raízes tá tudo ótimo.

O método de Ziegler-Nichols de malha aberta

O método de Ziegler-Nichols de malha fechada

Exercícios Resolvidos

Exercício Resolvido #1

Um sistema pode ser aproximado por um modelo de primeira ordem mais tempo morto (FOPDT) conforme a função de transferência abaixo:

Y s X s = 2 3 s + 1 e - 0,5 s

Considerando-se essa aproximação, sintonize um controlador P, PI e PID usando as correlaçãoes de decaimento ¼ de Ziegler-Nichols.

Passo 1

Para fazer as sintonias, primeiro precisamos determinar os parâmetro do modelos FOPDT comparando a função de transferência do enunciado com o modelo padrão abaixo:

K τ s + 1 e - θ s

Como para o sistema que estamos trabalhando a função de transferência é:

2 3 s + 1 e - 0,5 s

Podemos tirar que:

K = 2

τ = 3

θ = 0,5

Para fazer as sintonizações, devemos utilizar as correlações da tabela abaixo:

Passo 2

Um controlador do tipo P tem apenas um parâmetro a ser sintonizado:

G C s = K C

Usando a correlação da tabela:

K C = τ K θ = 3 0,5 * 2 = 3

Portanto

G C = 3

Passo 3

Para um controlador PI, teremos que calcular dois parâmetros:

G C s = K C 1 + 1 τ i s

Usando as correlações da tabela:

K C = 0,9 τ K θ = 0,9   3 0,5 * 2 = 2,7

τ i = θ 3 = 0,5 0 , 3 = 1 , 6 7

A função de transferência do controlador será então:

G C s = 2,7 1 + 1 1,67 s

Passo 4

Por último, um controlador PID precisa que três parâmetros sejam sintonizados:

G C s = K C 1 + 1 τ 1 s + τ d s

Seguindo as correlações de Ziegler-Nichols para um controlador PID:

K C = 1,2 τ K θ = 1,2   3 0,5 * 2 = 3,6

τ i = 2 θ = 2 * 0,5 = 1

τ d = 0,5 * θ = 0,5 * 0,5 = 0,25

O controlador será então:

G C s = 3,6 1 + 1 s + 0,25 s

Resposta

Controlador P:

G C = 3

Controlador PI:

G C s = 2,7 1 + 1 0,17 s

Controlador PID:

G C s = 3,6 1 + 1 s + 0,25 s

Exercício Resolvido #2

Um certo sistema, ao ser perturbado com uma entrada degrau unitária, apresenta a seguinte curva de reação

Projete um controlador PID para esse sistema utilizando os parâmetros de Zielger-Nichols de malha aberta.

Passo 1

Para definir os parâmetros do controlador, precisamos aproximar o sistema da curva a uma função de transferência de primeira ordem com tempo morto (FOPDT). Para isso, precisamos traçar a maior tangente com maior inclinação possível no gráfico do enunciado.

Essa tangente, será mais ou menos assim (não presisa ser extao, lembre-se, é uma aproximação):

Passo 2

Agora precisamos encontrar os parâmetros da função de transferência de primeira ordem com tempo morto:

K τ s + 1 e - θ s

O ganho K , é a variação total da saída. Como a saída variou de 0 até 1 , o ganho é 1 .

K = 1 - 0 = 1

O tempo morto, θ , é a distância entre o início da curva e o ponto onde a reta tangente que a nós traçamos cruza o valor inicial da respota.

Pelo nosso esboço, podemos estimar o ponto em que isso acontece como mais ou menos 0,2 . Sendo assim:

θ = 0,2 - 0 = 0,2   s e g

Por fim, a constante de tempo, τ , é aproximada pela diferença entre o tempo que a reta tangente cruza o valor final e o tempo que a reta tangente cruza o valor inicial. Portanto:

τ = 2,2 - 0,2 = 2,0   s e g

Sendo assim, nossa função de primeira ordem com tempo morto (FOPDT) aproximada ficou com a seguinte cara:

1 2 s + 1 e - 0,2 s

Passo 3

Para projetar nosso PID precisamos encontrar três parâmetros: K c , τ i e τ d . Podemos encontrar as correlações para esses três parâmetros na tabela de Ziegler Nichols para malha aberta

Como nosso controlador é PID, vamos usar a última linha para definir nossos parâmetros. Então:

K c = 1,2 τ θ K = 1,2 2 0,2 = 12

τ i = 2 * 0,2 = 0,4

τ d = 0,5 * 0,2 = 0,1

A função de transferência do controlador será então:

G c s = 12   1 + 1 0,4 s + 0,1 s

Resposta

G c s = 12   1 + 1 0,4 s + 0,1 s

Exercício Resolvido #3

Um certo sistema, ao ser perturbado com uma entrada degrau unitária, apresenta a seguinte curva de reação

Projete um controlador PI para esse sistema utilizando os parâmetros de Zielger-Nichols de malha aberta.

Passo 1

Vamos começar encondtrando a função FOPDT que aproxima o sistema:

Pela figura acima, podemos ver que o ponto onde a tangente cruza o valor inicial é aproximadamente 2 , portanto:

θ = 2

Já o ponto que a tangente cruz a linha do falor final da resposta é aproximadamente 7 , portanto:

τ = 7 - 2 = 5

E o ganho do sistema é a variação da resposta dividida pela amplitude do degrau, a variação da resposta nós podemos estimar como 4 , e como o degrau foi um degrau unitário, o ganho é simplesmente:

K = 4

Sendo assim, a nossa aproximação ficará:

Y s X s = 4 5 s + 1 e - 2 s

Passo 2

Para projetar nosso controlador PI precisamos encontrar dois parâmetros: K c e τ i . Podemos encontrar as correlações para esses parâmetros na tabela de Ziegler Nichols para malha aberta

Como nosso controlador é PI, vamos usar a segunda linha para definir nossos parâmetros. Então:

K c = 0,9 τ θ K = 0,9 5 4 * 2 = 0,625

τ i = θ 0,3 = 2 0,3 = 6,67

A função de transferência do controlador será então:

G c s = 0,625   1 + 1 6,67 s

Resposta

G c s = 0,625   1 + 1 6,67 s

Exercício Resolvido #4

Um certo sistema, ao ser perturbado com uma entrada degrau unitária, apresenta a seguinte curva de reação

Projete um controlador PID para esse sistema utilizando os parâmetros de Zielger-Nichols de malha aberta.

Passo 1

Vamos começar encondtrando a função FOPDT que aproxima o sistema:

Pela figura acima, podemos ver que o ponto onde a tangente cruza o valor inicial é aproximadamente 1 , portanto:

θ = 1

Já o ponto que a tangente cruz a linha do falor final da resposta é aproximadamente 12, portanto:

τ = 12 - 1 = 11

E o ganho do sistema é a variação da resposta dividida pela amplitude do degrau, a variação da resposta nós podemos estimar como 5, e como o degrau foi um degrau unitário, o ganho é simplesmente:

K = 5

Sendo assim, a nossa aproximação ficará:

Y s X s = 5 11 s + 1 e - s

Passo 2

Para projetar nosso controlador PID precisamos encontrar três parâmetros: K c , τ i e τ d . Podemos encontrar as correlações para esses parâmetros na tabela de Ziegler Nichols para malha aberta

Como nosso controlador é PID, vamos usar a última linha para definir nossos parâmetros. Então:

K c = 1,2 τ θ K = 1,2 11 1 * 5 = 2,2

τ i = 2 * θ = 2 * 1 = 2

τ d = 0,5 * θ = 0,5 * 1 = 0,5

A função de transferência do controlador será então:

G c s = 2,2 1 + 1 2 s + 0,5 s

Resposta

G c s = 2,2 1 + 1 2 s + 0,5 s

Exercício Resolvido #5

Um sistema possui os seguintes valores para ganho ganho e período final:

K C u = 2,5

T u = 0,6

Sintonnize um controlador P, PI e PID para esse sistema usando as correlações de decaimento ¼ de Ziegler-NIchols em malha fechada.

Passo 1

Para fazer as sintonizações, usaremos a tabela abaixa que apresenta as correlações de Ziegler-Nichols em malha fechada:

Passo 2

Um controlador do tipo P tem apenas um parâmetro a ser sintonizado:

G C s = K C

Usando a correlação da tabela:

K C = 0,5 K C u = 0,5 * 2,5 = 1,25

Portanto

G C = 1,25

Passo 3

Para um controlador PI, teremos que calcular dois parâmetros:

G C s = K C 1 + 1 τ i s

Usando as correlações da tabela:

K C = 0 , 4 5 K C u = 0 , 4 5 * 2,5 = 1,125

τ i = T u 1,2 = 0,6 1,2 = 0,5

A função de transferência do controlador será então:

G C s = 1,125 1 + 1 0,5 s

Passo 4

Por último, um controlador PID precisa que três parâmetros sejam sintonizados:

G C s = K C 1 + 1 τ 1 s + τ d s

Seguindo as correlações de Ziegler-Nichols para um controlador PID:

K C = 0 , 60 K C u = 0 , 60 * 2,5 = 1 , 5

τ i = T u 2 = 0,6 2 = 0 , 3

τ d = 0,125 T u = 0,125 * 0,6 = 0 , 07 5

O controlador será então:

G C s = 1,5 1 + 1 0,3 s + 0,07 5 s

Resposta

Controlador P:

G C = 1 , 25

Controlador PI:

G C s = 1,125 1 + 1 0,5 s

Controlador PID:

G C s = 1,5 1 + 1 0,3 s + 0,075 s

Exercício Resolvido #6

Um sistema tem a seguinte função de transferência em malha aberta:

G O L = 20 s + 1 s + 2 s + 5

Sintonize um controlador PI usando o método de decaimento ¼ de Ziegler-Nichols em malha fechada.

Passo 1

Para aplicar o método de ZN em malha fechada, é preciso determinar o ganho e período final do sistema. Para isso, usaremos o método da substituição direta.

Dada a função de transferência em malha aberta, a equação característica da malha fechada sera´:

1 + G C G O L = 0

Sendo que para aplicar o método, devemos considerar a princípio um controlador apenas proporcional:

G C = K C

Então:

1 + 20 K C s + 1 s + 2 s + 5 = 0

s + 1 s + 2 s + 5 + 20 K C = 0

s 3 + 8 s 2 + 17 s + 10 + 20 K C = 0

Passo 2

Vamos então aplicar o método da substituição direta:

s = i ω u

K C = K C u

i ω u 3 + 8 i ω u 2 + 17 i ω u + 10 + 20 K C u = 0

- i ω u 3 - 8 ω u 2 + 17 i ω u + 10 + 20 K C u = 0

Devemos agora igualar tanto a parte real como a parte imaginária a 0 :

- i ω u 3 + 17 i ω u = 0

- 8 ω u 2 + 10 + 20 K C u = 0

Da parte imaginária temos:

- ω u 2 + 17 = 0

ω u = 17 = 4,12

Substituindo na parte real:

- 8 * 17 + 10 + 20 K C u = 0

K C u = 6,3

Além disso, sabendo que:

T u = 2 π ω u

T u = 1,52

Passo 3

Com os valores de K C u e T u calculados, é possível fazer a sintonia pelo método de ZN de malha fechada. Usaremos as correlações da tabela abaixo:

Como queremos um controlador PI, usaremos as correlações da segunda linha.

K C = 0,45 K C u = 0,45 * 6,3 = 2,835

τ i = T u 2 = 1,52 2 = 0,76  

Portanto, a função de transferência do controlador será:

G C s = 2,835 1 + 1 0,76 s

Resposta

G C s = 2,835 1 + 1 0,76 s

Exercício Resolvido #7

Considere o sistema mostrado na figura abaixo:

Calcule os ajustes de um PID para esse sistema utilizando o método de Ziegler-Nichols de malha fechada.

Passo 1

Para aplicar o método de ZN em malha fechada, é preciso determinar o ganho e período final do sistema. Para isso, usaremos o método da substituição direta.

A função de transferência dessa malha será:

Y s R s = G C G p 1 + G C G p

Onde G p é a função de transferência do processo (segundo bloco). O denominador dessa função de transferência é a nossa equação característica.

1 + G C G p = 0

Para fazer esse método, primeiramente consideramos o controlador como puramente proporcional e em seguida encontramos o ganho último pelo método da substituição direta:

G C = K C

G p = 0,5 1 + 2 s 1 + s 1 + 0,5 s

Portanto:

1 + 0,5 K C 1 + 2 s 1 + s 1 + 0 , 5 s = 0

1 + 2 s 1 + s 1 + 0,5 s + 0 , 5 K C = 0

s 3 + 3,5 s 2 + 3,5 s + 1 + 0,5 K C = 0

Passo 2

Vamos então aplicar o método da substituição direta:

s = i ω u

K C = K C u

i ω u 3 + 3,5 i ω u 2 + 3,5 i ω u + 1 + 0,5 K C u = 0

- i ω u 3 - 3,5 ω u 2 + 3,5 i ω u + 1 + 0 , 5 K C u = 0

Devemos agora igualar tanto a parte real como a parte imaginária a 0 :

- i ω u 3 + 3,5 i ω u = 0

- 3,5 ω u 2 + 1 + 0 , 5 K C u = 0

Da parte imaginária temos:

- ω u 2 + 3,5 = 0

ω u = 3,5 = 1,87

Substituindo na parte real:

- 3,5 * 3,5 + 1 + 0 , 5 K C u = 0

K C u = 22,5

Além disso, sabendo que:

T u = 2 π ω u

T u = 3,36

Passo 3

Com os valores de K C u e T u calculados, é possível fazer a sintonia pelo método de ZN de malha fechada. Usaremos as correlações da tabela abaixo:

Como queremos um controlador PID, usaremos as correlações da última linha.

K C = 0,60 K C u = 0 , 60 * 22,5 = 13,5

τ i = T u 2 = 3,36 2 = 1,68

τ d = 0,125 T u = 0,125 * 3,36 = 0 , 42

Portanto, a função de transferência do controlador será:

G C s = 13,5 1 + 1 1,68 s + 0,4 2 s

Resposta

G C s = 13,5 1 + 1 1,68 s + 0,42 s

Exercício Resolvido #8

Uma malha de controle de realimentação é representada pelo diagrama de blocos abaixo.

O ganho do processo é K = 2,5   % S T / % S C (porcentagem de saída do transmissor por porcentagem de saída do controlador) e as constantes de tempo são:

τ 1 = 5   m i n

τ 2 = 0,8   m i n

τ 3 = 0,2   m i n

Determine os parâmetros de sintonização do controlador para resposta de razão de decaimento de um quarto élo método de ganho final (ZN de malha fechada) para:

  1. Um controlador P
  2. Um controlador PI
  3. Um controlador PID

Passo 1

Para aplicar o método de ZN em malha fechada, é preciso determinar o ganho e período final do sistema. Para isso, usaremos o método da substituição direta.

A função de transferência dessa malha será:

Y s R s = G C G p 1 + G C G p

Onde G p é a função de transferência do processo (segundo bloco). O denominador dessa função de transferência é a nossa equação característica.

1 + G C G p = 0

Para fazer esse método, primeiramente consideramos o controlador como puramente proporcional e em seguida encontramos o ganho último pelo método da substituição direta:

G C = K C

G p = 2,5 5 s + 1 0,8 s + 1 0,2 s + 1

Portanto:

1 + 2,5 K C 5 s + 1 0,8 s + 1 0,2 s + 1 = 0

5 s + 1 0,8 s + 1 0,2 s + 1 + 2 , 5 K C = 0

0,8 s 3 + 5,16 s 2 + 6 s + 1 + 2,5 K C = 0

Passo 2

Vamos então aplicar o método da substituição direta:

s = i ω u

K C = K C u

0,8 i ω u 3 + 5,16 i ω u 2 + 6 i ω u + 1 + 2 , 5 K C u = 0

- 0,8 i ω u 3 - 5 , 16 ω u 2 + 6 i ω u + 1 + 2 , 5 K C u = 0

Devemos agora igualar tanto a parte real como a parte imaginária a 0 :

- 0,8 i ω u 3 + 6 i ω u = 0

- 5,16 ω u 2 + 1 + 2 , 5 K C u = 0

Da parte imaginária temos:

- 0,8 ω u 2 + 6 = 0

ω u = 7 , 5 = 2,74

Substituindo na parte real:

- 5,16 * 7,5 + 1 + 2 , 5 K C u = 0

K C u = 15,08

Além disso, sabendo que:

T u = 2 π ω u

T u = 2,30

Passo 3

Com os valores de K C u e T u calculados, é possível fazer a sintonia pelo método de ZN de malha fechada. Usaremos as correlações da tabela abaixo:

Primeiramente, vamos calcular os parâmetros para o controlador proporcional;

G C = K C

Pela tabela, teremos:

K C = 0,5 K C u = 0,5 * 15,08 = 7,54

G C = 7,54

Passo 4

Agora para o controlador PI:

G C = K C 1 + 1 τ i s

Pegando as correlações da tabela:

K C = 0,45 K C u = 0,45 * 15,08 = 6,79

τ i = T U 1,2 = 2,30 1,2 = 1,92

Portanto:

G C = 6,79 1 + 1 1,92 s

Passo 5

Por fim, vamos fazer o mesmo para o PID:

G C = K C 1 + 1 τ i s + τ d s

Olhando na última linha da tabela:

K C = 0,60 K C u = 0 60 * 15,08 = 9,05

τ i = T U 2 = 2,30 2 = 1,15

τ d = 0,125 T u = 0,125 * 2,30 = 0,29

Portanto:

G C = 9,05 1 + 1 1,15 s + 0,29 s

Resposta

  1. G C = 7,54
  2. G C = 6,79 1 + 1 1,92 s
  3. G C = 9,05 1 + 1 1,15 s + 0,29 s