Estabilidade e critério de Routh-Hurwitz

Ver Todos os Exercícios

Teoria

Introdução

Agora que nós já vimos um pouco melhor os diferentes tipos de função de transferência que nós podemos encontrar no nosso caminho, é hora de expandir um pouco nossos horizontes e olhar o nosso sistema de controle de uma forma mais ampla.

Às vezes podemos dar de cara com um cara desse tipo aqui:

Essa malha de controle vai ter função de transferência do tipo;

C s R s = 1 s + 1   . 3 2 s + 3 . 1 s + 1 s + 3 1 + 0,5 3 s + 1 . 1 s + 1 . 3 2 s + 3 . 1 s + 1 s + 3

E aí a gente vai esbarrar em um problema muito comum no controle de processos: essa malha é estável ou instável?

Em alguns casos, a nossa malha de controle pode ser instável, isso quer dizer que conforme o tempo passa, a resposta do processo ao invés de se aproximar de um valor determinado, vai se afastando cada vez mais desse valor. Quando isso acontece, nós falamos que o sistema é instável.

Bora aprender como a gente vê se a nossa malha é estável ou instável?

Controle feedback

Antes de a gente se aprofundar nos critérios de estabilidade, vamos dar uma relembrada em controle feedback.

O controle feedback é o mais clássico dos controles, é aquele que você tem que sacar tudo!

Nele, nós vamos ter duas variáveis principais: a variável controlada, e a variável manipulada. Basicamente a variável controlada é aquela que nós queremos manter em um determinado setpoint, enquanto a variável manipulada é a variável que nós vamos mexer para manter a variável controlada no nosso setpoint.

O diagrama de blocos simplificado do controle feedback é assim:

Vamos entender o que está acontecendo!

G c , G v , G p e G m são as funções de transferência do controlador, da válvula, do processo e do sensor-transmissor, respectivamente.

Presta atenção agora: após o nosso processo acontecer (em G p ) nós temos um valor para a nossa variável controlada C s , esse valor será medido por um sensor transmissor ( G m ) e comparado com um valor de setpoint ( R ( s ) ), após fazer essa comparação o controlador ( G c ) vai interpretar essa diferença entre variável medida e setpoint e tomar uma decisão. Quando ele toma essa decisão, ele envia um sinal para a válvula ( G v ) que altera a variável manipulada, fazendo que aconteça uma mudança no processo que altera o valor da variável controlada. E tudo isso acontece novamente. Incrível né?

Para o controle feedback tradicional, a função de transferência da malha fica:

C s R s = G c G v G p 1 + G c G v G p G m

Caso você não tenha entendido como chegamos a essa expressão, dá uma revisadinha na aula de álgebra de diagrama de blocos, você vai ver que não tem mistério.

Equação característica

Agora que já demos uma revisada em controle feedback, vamos começar a falar de fato de estabilidade.

Sempre que formos avaliar a estabilidade de um sistema, devemos olhar para a sua equação característica.

AH MAS O QUE DIABOS É EQUAÇÃO CARACTERÍSTICA? ”

Calma, não precisa gritar, já te explico.

Equação característica, nada mais é do que a expressão que fica no denominador da nossa função de transferência. Então, quando a gente estiver trabalhando com controle feedback, nossa equação característica vai ser:

1 + G c G v G p G m = 0

Para o exemplo que nós demos lá no início, ficaria assim:

1 + 0,5 3 s + 1 . 1 s + 1 . 3 2 s + 3 . 1 s + 1 s + 3 = 0

Rearranjando:

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

Na maioria dos casos a equação característica será um polinômio de ordem n . Nesse caso é um polinômio de ordem 5 .

Esse polinômio vai conter os segredos que vai nos dizer se o nosso sistema é estável ou instável.

O primeiro critério que a gente vai poder observar logo de cara é o critério dos sinais dos coeficientes.

O critério dos coeficientes diz que uma malha só é estável se todos os coeficientes da função de transferência tiverem o mesmo sinal.

Mas fique atento, esse critério é necessário, porém não é suficiente. Isso quer dizer que se o critério não for atendido, a malha necessariamente será instável, porém, se ele for atendido, não é possível determinar nada sobre a estabilidade do sistema.

Para garantir a estabilidade de um processo então, é necessário que todas as raízes da equação característica sejam números reais negativos ou números imaginários com a parte real negativa.

Em equações de até terceira ordem, talvez seja fácil determinar essas raízes, entretanto, em equações de ordem superior a coisa complica um pouco, e por isso precisamos de ferramentas para determinar de maneira mais fácil se o sistema é estável ou não.

Critério de Routh-Hurwitz

Uma ferramenta que temos para determinar a estabilidade de uma malha é o critério de Routh-Hurwitz.

Com esse critério, nós construiremos um diagrama, e esse diagrama será o que nos dirá se o sistema é estável ou instável.

Bora colocar a mão na massa então.

Imagina que a gente tenha uma equação característica de quinta ordem:

a 5 s 5 + a 4 s 4 + a 3 s 3 + a 2 s 2 + a 1 s 1 + a 0 = 0

Então, para construir nosso digrama nós colocaremos nossos coeficientes na matriz da seguinte forma:

Isso quer dizer que os coeficientes das duas primeiras linhas serão preenchidos pelos coeficientes da equação característica. Para o exemplo da equação de quinta ordem:

Após preencher as duas primeiras linhas com coeficientes da equação característica, devemos calcular os outros termos do diagrama até a linha n + 1 . Esses termos podem ser calculados da seguinte forma:

Para a linha 3 :

b 1 = a 4 a 3 - a 5 a 2 a 4

b 2 = a 4 a 1 - a 5 a 0 a 4

E assim sucessivamente.

Para a linha 4:

c 1 = b 1 a 2 - a 4 b 2 b 1

c 2 = b 1 a 0 - a 4 b 3 b 1

E assim sucessivamente em todas as linhas até todos os termos serem 0 .

Após montar todo o diagrama, o critério de Routh-Hurwitz diz que o sistema é estável se todos os coeficientes da primeira coluna forem positivos.

Cada mudança de sinal na primeira coluna, indica uma raiz da equação característica no semiplano direito (parte real positiva) e isso torna o sistema instável

Juro que na prática é mais fácil do que parece, bora praticar?

Controle feedback

Equação característica

Critério de Routh-Hurwitz

Exercícios Resolvidos

Exercício Resolvido #1

Defina a estabilidade do sistema representado pela malha feedback abaixo:

Passo 1

O primeiro passo, sempre que queremos determinar a estabilidade de uma malha, é encontrar a equação característica da malha. Para esse sistema, como estamos trabalhando com um feedback clássicozão, a função de transferência ficará assim:

C s R s =   3 2 s + 1 2 s + 1 s - 5 1 s + 1 s + 3 1 + 3 2 s + 1 2 s + 1 s - 5 1 s + 1 s + 3 3  

Caso você tenha dúvidas em como obter essa expressão, dá uma revisadinha na aula de álgebra de diagrama de blocos.

Nós precisamos da equação característica, e ela sempre será o denominador da função de transferência. Então:

1 + 3 2 s + 1 2 s + 1 s - 5 1 s + 1 s + 3 3 = 0

Passo 2

Vamos então deixar essa equação mais bonitinha pra gente poder pensar melhor sobre a estabilidade. Multiplicando toda a equação pelo denominador:

2 s + 1 s - 5 s + 1 s + 3 + 18 s + 1 = 0

Fazendo as multiplicações:

2 s 4 - s 3 - 35 s 2 - 29 s + 3 = 0

Passo 3

Agora que já temos a nossa equação característica, poderíamos determinar a estabilidade dela pelo critério de Routh-Hurwitz. Entretanto, só de olhar pra carinha dela, já conseguimos ter um spoiler da nossa resposta.

O critério dos coeficientes diz que para o sistema ser estável, a sua equação característica precisa ter todos os sinais iguais.

Nesse caso, o critério não foi atendido. Portanto o sistema é instável.

Obs: só lembrando que esse critério é necessário, porém não suficiente.

Resposta

O sistema é instável de acordo com o critério dos coeficientes.

Exercício Resolvido #2

Pelo critério de Routh-Hurwitz, determine se o sistema cuja função de transferência está representada abaixo é estável:

G s = 2 s 3 + 2 s 2 + 3 s + 4

Passo 1

A questão é bem direta: precisamos definir a estabilidade do sistema.

A nossa equação característica será:

s 3 + 2 s 2 + 3 s + 4 = 0

E precisamos aplicá-la no diagram de Roth-Hurwitz abaixo:

Passo 2

O primeiro passo seria então, distribuir os coeficientes nas duas primeiras linhas:

E então, podemos começr a calcular os outros termos.

Para a terceira linha:

b 1 = a 2 a 1 - a 3 a 0 a 2 = 2 3 - 1 4 2 = 1

b 2 = 2 0 - 1 0 2 = 0

Colocando no diagrama:

Passo 3

Calculando agora as linhas 4 e 5:

c 1 = 1 4 - 2 0 1 = 4

c 2 = 0

d 1 =   4 0 - 1 0 4 = 0

Colocando tudo no diagrama:

Resposta

O sistema é estável. Não ocorreu mudança de sinal na primeira coluna do diagarama.

Exercício Resolvido #3

Construa o diagrama de Routh para o sistema mostrado abaixo:

C s R s = 1000 s 3 + 10 s 2 + 31 s + 1030

Passo 1

Como já temos a função de transferência da malha, nossa função é simplesmente:

s 3 + 10 s 2 + 31 s + 1030 = 0

Devemos então calcular o diagrama de Routh do sistema, seguindo o modelo:

Passo 2

O primeiro passo, será preencher a primeira e a segunda linha com os coeficientes da equação característica:

Passo 3

Em seguida, iremos calcular os valores dos coeficientes da segunda linha:

b 1 = a 2 a 1 - a 3 a 0 a 2 = 10 31 - 1 1030 10 = - 72

b 2 = 10 0 - 1 0 10 = 0

Colocando os valores no diagrama:

Para a linha 4:

c 1 = - 72 1030 - 10 0 - 72 = 1030  

c 2 = 0

E para a linha 5:

d 1 = 0

Colocando tudo no diagrama:

Resposta

Exercício Resolvido #4

Determine se o sistema representado pela função de transferência abaixo é estável:

C s R s = 1 2 s 4 + 3 s 3 + s 2 + 2 s + 1

Passo 1

A nossa equação característica é:

2 s 4 + 3 s 3 + s 2 + 2 s + 1 = 0

E devemos montar nosso diagrama de Routh, seguindo o modelo para equação com 5 termos:

Passo 2

Primeiramente, vamos colocar os coeficientes da equação característica nas duas primeiras linhas:

Passo 3

Vamos então calcular os faores da linha 3:

b 1 = a 4 a 3 - a 5 a 2 a 4 = 3 1 - 2 2 3 = - 1 3

b 2 = a 4 a 1 - a 5 a 0 a 4 = 3 1 - 2 0 3 = 1

b 3 = 0

Colocando no diagrama:

Passo 4

Calculando os fatores da quarta linha:

c 1 = b 1 a 2 - a 4 b 2 b 1 = - 1 3 2 - 3 1 - 1 3 = 11

c 2 = b 1 a 0 - a 4 b 3 b 1 = - 1 3 0 - 3 0 - 1 3 = 0

E os termos da quinta e sexta linha:

d 1 = 11 1 - - 1 3 0 11 = 1

d 2 = 0

e 1 = 0

Colocando no diagrama:

O sistema não é estável pois ocorreu mudança de sinal na primeira coluna do diagrama de Routh-Hurwitz.

Resposta

Sistema instável.

Exercício Resolvido #5

Determine a estabilidade da função de transferência em malha fechada:

T s = 10 s 5 + 2 s 4 + 4 s 3 + 6 s 2 + 5 s + 3

Passo 1

O primeiro passo será determinar a nossa equação característica, como a função de transferência dada já é a da malha, a equação característica é:

s 5 + 2 s 4 + 4 s 3 + 6 s 2 + 5 s + 3 = 0

Em seguida devemos utilizar o diagrama e Routh-Hurwitz para determinar a estabilidade do sistema:

Passo 2

Vamos então dispor os coeficientes da equação característica nas duas primeiras linhas do diagrama de Routh-Hurwitz:

Passo 3

Em seguida, calculamos os coeficientes da linha 3:

b 1 = a 4 a 3 - a 5 a 2 a 4 = 2 4 - 1 6 2 = 1

b 2 = a 4 a 1 - a 5 a 0 a 4 = 2 5 - 1 3 2 = 3,5

b 3 = 2 0 - ( 1 ) ( 0 ) 2 = 0

Substtituindo os valores encontrados no diagrama:

Passo 4

Seguindo os cálculos agora, vamos para a linha 4:

c 1 = b 1 a 2 - a 4 b 2 b 1 = 1 6 - 2 3,5 1 = - 1

c 2 = b 1 a 0 - a 4 b 3 b 1 = 1 3 - 2 0 1 = 3

Não precisamos calcular os próximos termos dessa linha, pois eles serão 0 .

Colocando no diagrama:

Passo 5

Calculando as últimas linhas agora:

d 1 = - 1 3,5 - 1 3 - 1 = 6,5

d 2 = 0

e 1 = 6,5 3 - - 1 0 6,5 = 3

e 2 = 0

f 1 = 0

No diagrama:

Sendo assim, o sistema é instável, pois apresenta duas mudanças de sinal na primeira coluna do diagrama, indicando a presença de duas raízes no semi-plano direito (parte real positiva).

Resposta

Sistema instável.

Exercício Resolvido #6

Determine se a malha abaixo é estável:

Passo 1

Primeiramente devemos definir nossa equação característica. Para isso, precisamos da função de transferência da malha:

C s R s = G s 1 + G s

Com G ( s ) sendo a função de transferência do bloco.

A nossa equação característica será então:

1 + G s = 0

1 + 3 s 4 + 2 s 3 + 3 s 2 + s + 10 = 0  

Rearranjando:

s 4 + 2 s 3 + 3 s 2 + s + 13 = 0

Agora que já temos nossa equação característica, precisamos aplicá-la no modelo abaixo para conseguir determinar a estabilidade da malha:

Passo 2

Vamos então, dispor os coeficientes da equação característica nas duas primeiras linhas:

Vamos calcular agora os elementos da linha 3:

b 1 = a 3 a 2 - a 4 a 1 a 3 = 2 3 - 1 1 2 = 2,5

b 2 = a 3 a 0 - a 4 a - 1 a 3 = 2 13 - 1 0 2 = 13

b 3 = 0

Colocando no diagram:

Passo 3

Calculando agora os elementos da terceira linha:

c 1 = b 1 a 1 - a 3 b 2 b 1 = 2,5 1 - 2 13 2,5 = - 9,4

c 2 = b 1 a - 1 - a 3 b 3 b 1 = 2,5 0 - 2 0 0 = 0

c 3 = 0

Passando para o diagrama:

Passo 4

Agora para as duas últimas linhas:

d 1 = - 9,4 13 - 2,5 0 - 9,4 = 13  

d 2 = d 3 = 0

e 1 = 0

Completando o diagrama:

Resposta

O sistema é instável pois ocorre mudança de sinais na primeira coluna do diagrama de Routh-Hurwitz.

Exercício Resolvido #7

Determine se o sistema abaixo é estável:

Passo 1

Primeiro precisamos definir nossa equação característica. Já vimos que para o feedback unitário, a equação característica fica da seguinte forma:

1 + G 1 s G 2 s = 0

Onde G 1 s e G 2 ( s ) são os blocos do diagarama. Então:

1 + 2 s s + 3 10 s + 1 2 = 0

Rearranjando:

s s + 3 s + 1 2 + 20 = 0

s s + 3 ( s 2 + 2 s + 1 ) + 20 = 0

s s 3 + 5 s 2 + 7 s + 3 + 20 = 0

s 4 + 5 s 3 + 7 s 2 + 3 s + 20 = 0

Com a nossa equação característica em mãos, só precisamos aplicar o critério de Routh-Hurwitz para definir a estabilidade:

Passo 2

Vamos substituir os termos no diagrama então:

E vamos calcular os termos da terceira linha conforme as fórmulas:

b 1 = a 3 a 2 - a 4 a 1 a 3 = 5 7 - 1 3 5 = 6,4

b 2 = a 3 a 0 - a 4 a - 1 a 3 = 5 20 - 1 0 5 = 20

b 3 = 0

Vamos colocar esses termos no diagrama:

Passo 3

Vamos passar então para a quarta linha:

c 1 = b 1 a 1 - a 3 b 2 b 1 = 6,4 3 - 5 20 6,4 = - 12,625

c 2 = b 1 a - 1 - a 3 b 3 b 1 = 6,4 0 - 5 0 6,4 = 0

c 1 = 0

Então:

Passo 4

Para finalizar. Calculamos os coeficientes das duas últimas linhas:

d 1 = - 12,625 20 - 6,4 0 - 12,625 = 20

d 2 = 0

e 1 = 20 0 - - 12,625 0 20 = 0

Completando nosso diagrama:

Resposta

O sistema é instável porque existe mudança de sinal na primeira coluna do diagrama de Routh-Hurwitz.

Exercício Resolvido #8

Caracterize a malha abaixo quanto à sua estabilidade:

Passo 1

Primeiramente, vamos definir a função de transferência equivalente da malha, para que possamos encontrar a nossa equação característica.

Como esse é aquele diagrama feedback mais comum, a função de transferência será:

C s R s =   3 2 s + 1 2 s + 1 s + 2 1 s + 1 s + 3 1 + 3 2 s + 1 2 s + 1 s + 2 1 s + 1 s + 3 3  

Caso você tenha dúvidas em como obter essa expressão, dá uma revisadinha na aula de álgebra de diagrama de blocos.

Como a equação característica é sempre o denominador da função de transferência, nossa equação característica será:

1 + 3 2 s + 1 2 s + 1 s + 2 1 s + 1 s + 3 3 = 0

Vamos então rearranjar essa equação para conseguir usar o critério de Routh-Hurwitz. Multiplicando toda a equação pelo denominador:

2 s + 1 s + 2 s + 1 s + 3 + 18 s + 1 = 0  

Fazendo as multiplicações:

2 s 4 + 13 s 3 + 28 s 2 + 41 s + 24 = 0

Por fim, devemos substituir os termos no diagrama de Routh-Hurwitz, conforme o modelo:

Passo 2

Substituindo os termos nas duas primeiras linhas:

E vamos calcular os termos da terceira linha conforme as fórmulas:

b 1 = a 3 a 2 - a 4 a 1 a 3 = 13 28 - 2 41 13 = 21,7

b 2 = a 3 a 0 - a 4 a - 1 a 3 = 13 24 - 2 0 13 = 24

b 3 = 0

Colocando esses resultados no diagrama:

Passo 3

Vamos passar então para a quarta linha:

c 1 = b 1 a 1 - a 3 b 2 b 1 = 21,7 41 - 13 24 21,7 = 26,6

c 2 = b 1 a - 1 - a 3 b 3 b 1 = 21,7 0 - 13 0 21,7 = 0

c 1 = 0

Então:

Passo 4

Para finalizar. Calculamos os coeficientes das duas últimas linhas:

d 1 = 26,6 24 - 21,7 0 26,6 = 24

d 2 = 0

e 1 = 12 0 - - 6,375 0 12 = 0

Completando nosso diagrama:

Resposta

O sistema é estável porque não existe mudança de sinal na primeira coluna do diagrama de Routh-Hurwitz.