Substituição direta

Ver Todos os Exercícios

Teoria

Introdução

No tópico anterior, nós vimos a técnica do lugar geométrico das raízes. Essa técnica nos mostrava em qual lugar do plano se encontrariam as raízes da equação característica, conforme nós variamos um parâmetro, o ganho por exemplo.

Agora nós veremos uma técnica parecida. O método da substituição direta vai ser uma ferramenta muito importante no nosso estudo do controle de processos. Esse método vai nos mostrar o limite superior do ganho da nossa malha.

Então por exemplo. Dado o seguinte sistema:

Nós vamos obter:

K u = 23,8   % S C % S T

T u = 28,7   s

Que são os valores de ganho e período final, ou seja, os valores de ganho e período que fazem com que a resposta do sistema esteja no limite da estabilidade, ou em outras palavras, que o sistema seja marginalmente estável. Isso quer dizer que se aumentarmos um pouquinho que seja o ganho, o sistema já se torna instável.

O índice u em K u e T u vem da palavra em inglês ultimate. E serve justamente para indicar o ganho ou período final (último).

Antes de mergulharmos de cabeça no método vamos relembrar umas coisas que nós vimos no nosso tópico de estabilidade. Vou destacar especialmente esses dois pontos:

  • Para ser estável, um sistema precisa que todos os coeficientes da sua equação característica tenha o mesmo sinal. Esse é um critério necessário, porém não suficiente, o que quer dizer que caso o sistema não atenda a esse critério ele necessariamente é instável, porém, se ele atender não podemos determinar nada sobre sua estabilidade a partir desse critério;
  • Para ser estável, todas as raízes da equação característica precisam necessariamente estar no semiplano complexo da esquerda, ou seja, terem a sua parte real negativa. Caso qualquer raiz da equação característica apresente parte real positiva, o sistema se torna automaticamente instável.

O método da substituição direta

Então vimos que para que um sistema seja estável, as raízes de sua equação característica precisam se encontrar todas no semiplano complexo da esquerda. Graficamente, isso quer dizer que elas precisam se encontrar na seguinte região:

Ou seja, a raiz precisa ser um número com a parte real negativa, ou um número puramente real negativo se ele se encontrar sobre o eixo dos reais.

Pelo gráfico podemos ver que quando a raiz for um número imaginário puro (estiver sobre o eixo dos imaginários) o sistema estará no seu limite de estabilidade, ou seja, qualquer variação para a direita o tornará instável.

Então o método da substituição direta vai forçar que a raiz da equação esteja sobre o eixo dos imaginários, levando o sistema ao seu limite e encontrando os valores dos seus parâmetros que causam essa situação.

Por exemplo, uma raiz qualquer tem a seguinte fórmula geral:

s = σ + i ω

Onde σ é a parte real, e ω é a parte imaginária.

Se quisermos encontrar o limite de estabilidade do nosso sistema, assumiremos que:

s = i ω u

Ou seja, a raiz se encontra sobre o eixo dos imaginários.

Após fazer as substituições devemos separar a parte real da parte imaginária e encontrar os valores de ganho e frequência final que satisfazem ao sistema.

Com o valor de ω u é possível encontrar T u através da fórmula:

T u = 2 π ω u

Para organizar tudo na sua cabeça eu fiz um passo a passo:

  • Obter a equação característica do seu sistema;
  • Fazer as substituições:
  • s = i ω u

    K = K u

  • Separar a parte real da parte imaginária e igualar ambas as partes a 0 ;
  • Resolver o sistema resultante e encontrar os valores de K u e ω u ;
  • Encontrar o valor de período final pela fórmula:

T u = 2 π ω u

Entendeu tudinho? Não entendeu? Acha que sim? Não faz pergunta difícil? Calma, vamos fazer um exemplo para tentar ver como tudo isso funciona na prática.

Botando a mão na massa

Vamos retomar então nosso exemplo lá do início:

O primeiro passo é determinar nossa equação característica. Como esse é um feedback comunzão, nossa equação característica vai ficar assim:

1 + K 1 10 s + 1 50 30 s + 1 0,016 ( 3 s + 1 ) = 0

Rearranjando:

10 s + 1 30 s + 1 3 s + 1 + 0,8 K = 0

Multiplicando os termos em parênteses:

900 s 3 + 420 s 2 + 43 s + 1 + 0,8 K = 0

Agora que já encontramos a nossa equação característica podemos fazer nossas substituições:

s = i ω u

K = K u

Substituindo na equação característica:

900 i ω u 3 + 420 i ω u 2 + 43 i ω u + 1 + 0,8 K u = 0

900 i 3 ω u 3 + 420 i 2 ω u 2 + 43 i ω u + 1 + 0,8 K u = 0

Sabendo que:

i 2 = - 1

i 3 = i 2 i = - 1 i = - i

Temos:

- 900 i ω u 3 - 420 ω u 2 + 43 i ω u + 1 + 0,8 K u = 0

Rearranjando:

- 900 ω u 3 + 43 ω u i - 420 ω u 2 + 1 + 0,8 K u = 0

Para que a igualdade acima seja verdadeira, tanto a parte real quanto a parte imaginário precisam ser iguais a 0 , sendo assim:

- 900 ω u 3 + 43 ω u = 0

- 420 ω u 2 + 1 + 0,8 K u = 0

Agora precisamos resolver esse sistema. Como a primeira equação só depende as frequência, podemos começar por ela

- 900 ω u 3 + 43 ω u = 0

- ω u ( 900 ω u 2 - 43 ) = 0

O resultado em que ω u = 0 sempre será irrelevante para nós, então:

900 ω u 2 - 43 = 0

ω u = 43 900 = 0,2186   r a d / s

Agora, com o valor de ω u , encontramos K u na segunda equação:

- 420 0,2186 2 + 1 + 0,8 K u = 0

K u = 23,8 % S C % S T

Aproximação de Padé

Em alguns casos, podemos encontrar processos com tempo morto.

O tempo morto é representado no domínio frequência como uma função exponencial:

e - t 0 s  

Onde t 0 é o tempo morto.

Quando temos esse termo, não conseguimos resolver a nossa substituição direta normalmente porque a exponencial bagunça tudo.

Então a gente, que não é bobo nem nada, dá uma “trapaceada” e usa uma aproximação conhecida como aproximação de Padé.

Essa aproximação, some com o termo exponencial e transfora ele em um termo mais amigável que é bem mais fácil de trabalhar:

e - t 0 s ≅ 1 - t 0 2 s 1 + t 0 2 s

Então se nós o bloco do transmissor do nosso exemplo tivesse um tempo morto de 10   s e g , ele seria escrito da seguinte forma:

G m = 0,016 3 s + 1 e - 10 s

Quando fôssemos aplicar o método da substituição direta, não poderemos fazer com a função do jeito que está, precisaremos então aplicar a aproximação de Padé. A função de transferência ficaria então:

G m = 0,016 3 s + 1 1 - 10 2 s 1 + 10 2 s = 0,016 3 s + 1 1 - 5 s 1 + 5 s

O método da substituição direta

Botando a mão na massa

Aproximação de Padé

Exercícios Resolvidos

Exercício Resolvido #1

Um processo pode ser definido pela seguinte função de transferência:

G 1 s = K 5 s - 1 τ V s + 1 τ T s + 1

Onde τ V e τ T são as constantes de tempo da válvula de controle e do transmissor, respectivamente. Considerando um controlador proporcional, encontre a faixa do ganho da malha, K K C para a qual a malha é estável se τ V = 0,1   m i n e τ T = 1,0   m i n

Passo 1

Primeiramente devemos identificar noss malha para podermos encontrar a equação característica dela.

Como a função de transferência dada no exercício já envolve as dinâmicas da válvula e do transmissor, a malha pode ser simplificada da seguite maneira:

Onde G 1 ( s ) é a função de transferência do processo, válvula e transmissor que foi dada no enunciada e G C ( s ) é a função de transferência do controlador.

Como o enunciado diz que o controlador [e do tipo proporcional:

G C = K C

E como a malha é uma malha do tipo feedback, a função de transferência equivalente fica:

C s R s = G C G 1 1 + G C G 1

Mas nós estamos interessados na equação característica, que nada mais é do que o denominador da função de transferência. Sendo assim, nossa equação de transferência para esse caso será:

1 + G C G 1 = 0

Passo 2

Com os valores das constantes dados, nossa função de transferência fica assim:

G 1 s = K 5 s - 1 s + 1 ( 0,1 s + 1 )

Substituindo na equação característica:

1 + G C G 1 = 0

1 + K C K 5 s - 1 s + 1 ( 0,1 s + 1 ) = 0

Multiplicando tudo pelo denominador:

5 s - 1 s + 1 0,1 s + 1 + K K c = 0

0,5 s 3 + 5,4 s 2 + 3,9 s - 1 + K K C u = 0

Passo 3

Uma vez que temos a equação característica bonitinha, o próximo passo é fazer a substituição:

s = i ω u

K K C = K K C u

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

Como:

i 2 = - 1

i 3 = - i

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

Agora devemos separar a parte real da imaginária e igualar as duas a 0 :

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

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

Da parte imaginária temos:

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

- 0,5 ω u 2 + 3,9 = 0

ω u = 3,9 0,5 = 2,79   r a d / s

Substituindo na parte real:

- 5,4 2,79 2 - 1 + K K C u = 0

K K C u = 43,12

Portanto:

K K C < 43,12

Passo 4

O enunciado pede que a gente encontre a faixa de estabilidade pra K K C , entretanto achamos apenas o limite superior.

Para achar o limite inferior, usamos o critério dos coeficientes, que diz que para a malha ser estável todos os coeficientes devem possuir o mesmo sinal.

Sendo assim:

- 1 + K K C > 0

K K C > 1

Portanto, a faixa de estabilidade para o ganho da malha é:

1 < K K C < 43,12

Resposta

1 < K K C < 43,12

Exercício Resolvido #2

Para uma malha de realimentação representada pelo diagrama de blocos abaixo:

Determine o ganho e período finais para um controlador proporcional para a função de transferência de processo abaixo:

G 1 s = 1 4 s + 1 2 s + 1 s + 1

Passo 1

Como a malha que estamos trabalhando é do tipo feedback, a função de transferência equivalente da malha é:

C s R s = G C G 1 1 + G C G 1

Mas nós estamos interessados na equação característica, que nada mais é do que o denominador da função de transferência. Sendo assim, nossa equação de transferência para esse caso será:

1 + G C G 1 = 0

Passo 2

Agora devemos substituir as funções de transferência do controlador e do processo na nossa equação característica:

G C = K C

G 1 s = 1 4 s + 1 2 s + 1 s + 1

Substituindo na equação característica:

1 + G C G 1 = 0

1 + K C 1 4 s + 1 2 s + 1 s + 1 = 0

Multiplicando tudo pelo denominador:

4 s + 1 2 s + 1 s + 1 + K c = 0

8 s 3 + 14 s 2 + 7 s + 1 + K C = 0

Passo 3

Finalmente então, podemos aplicar a nossa substituição para encontrar o ganho e período finais:

s = i ω u

K C = K C u

8 i ω u 3 + 14 i ω u 2 + 7 i ω u + 1 + K C u = 0

Como:

i 2 = - 1

i 3 = - i

- 8 i ω u 3 - 14 ω u 2 + 7 i ω u + 1 + K C u = 0

Separando a parte real da imaginária e igualando as duas a 0 :

- 8 i ω u 3 + 7 i ω u = 0

- 14 ω u 2 + 1 + K C u = 0

Da parte imaginária temos:

- 8 i ω u 3 + 7 i ω u = 0

- 8 ω u 2 + 7 = 0

ω u = 7 8 = 0,93   r a d / s

Substituindo na parte real:

- 14 0,93 2 + 1 + K C u = 0

K C u = 11,25   % S C % S T

Para determinar o períofo final, devemos aplicar a fórmula:

T u = 2 π ω u

Logo:

T u = 2 π 0,93 = 6,72   m i n

Resposta

K C u = 11,25 % S C % S T

T u = 6,72   m i n

Exercício Resolvido #3

Para uma malha de realimentação representada pelo diagrama de blocos abaixo:

Determine o ganho e período finais para um controlador proporcional para a função de transferência de processo abaixo:

G 1 s = ( 0,5 s + 1 ) 4 s + 1 2 s + 1 s + 1

Passo 1

Como a malha que estamos trabalhando é do tipo feedback, a função de transferência equivalente da malha é:

C s R s = G C G 1 1 + G C G 1

Mas nós estamos interessados na equação característica, que nada mais é do que o denominador da função de transferência. Sendo assim, nossa equação de transferência para esse caso será:

1 + G C G 1 = 0

Passo 2

Como estamos trabalhando com um controlador proporcional, a função de transferência para o controlador será:

G C = K C

E a do processo:

G 1 s = ( 0,5 s + 1 ) 4 s + 1 2 s + 1 s + 1

Colocando tudo na equação característica:

1 + G C G 1 = 0

1 + K C ( 0,5 s + 1 ) 4 s + 1 2 s + 1 s + 1 = 0

Multiplicando tudo pelo denominador:

4 s + 1 2 s + 1 s + 1 + K c ( 0,5 s + 1 ) = 0

8 s 3 + 14 s 2 + 7 s + 1 + K C 0,5 s + K C = 0

8 s 3 + 14 s 2 + 7 + 0,5 K C s + 1 + K C = 0

Passo 3

Fazendo aquela substituição marota:

s = i ω u

K C = K C u

8 i ω u 3 + 14 i ω u 2 + ( 7 + 0,5 K C u ) i ω u + 1 + K C u = 0

Como:

i 2 = - 1

i 3 = - i

- 8 i ω u 3 - 14 ω u 2 + ( 7 + 0,5 K C u ) i ω u + 1 + K C u = 0

Devemos então separar a parte real da imaginária e igualar as duas a 0 :

- 8 i ω u 3 + ( 7 + 0,5 K C u ) i ω u = 0

- 14 ω u 2 + 1 + K C u = 0

Da parte imaginária temos:

- 8 i ω u 3 + 7 + 0,5 K C u i ω u = 0

- 8 ω u 2 + 7 + 0,5 K C u = 0

ω u = 7 + 0,5 K C u 8

Substituindo na parte real:

- 14 7 + 0,5 K C u 8 2 + 1 + K C u = 0

- 14 7 + 0,5 K C u 8 + 1 + K C u = 0

- 14 7 + 0,5 K C u + 8 = 8 K C u = 0

- 98 - 7 K C u + 8 + 8 K C u = 0

K C u = 90   % S C % S T

Substituindo o valor de K C u :

ω u = 7 + 0,5 K C u 8

ω u = 7 + 45 8 = 2,55 r a d / m i n

Para determinar o períofo final, devemos aplicar a fórmula:

T u = 2 π ω u

Logo:

T u = 2 π 2,55 = 2,46   m i n

Resposta

K C u = 90 % S C % S T

T u = 2,46   m i n

Exercício Resolvido #4

Para uma malha de realimentação representada pelo diagrama de blocos abaixo:

Determine o ganho e período finais para um controlador proporcional para a função de transferência de processo abaixo:

G 1 s = e - 0,6 s 6 s + 1

Passo 1

Como a malha que estamos trabalhando é do tipo feedback, a função de transferência equivalente da malha é:

C s R s = G C G 1 1 + G C G 1

Mas nós estamos interessados na equação característica, que nada mais é do que o denominador da função de transferência. Sendo assim, nossa equação de transferência para esse caso será:

1 + G C G 1 = 0

Passo 2

Como estamos trabalhando com um controlador proporcional, a função de transferência para o controlador será:

G C = K C

Porém, aqui esbarramos com um problema. A função do processo tem um termo exponencial: o famoso tempo morto.

Sempre que precisarmos fazer substituição direta em uma função de transferência com tempo morto nós acionaremos o nosso megazord chamado aproximação de Padé.

Relembrando a fórmula da aproximação de Padé:

e - t 0 s ≅ 1 - t 0 2 s 1 + t 0 2 s

Para o nosso exercício então:

e - 0,6 s ≅ 1 - 0,6 2 s 1 + 0,6 2 s = 1 - 0,3 s 1 + 0,3 s

A nossa função de transferêncua para o processo será então:

G 1 s = 1 - 0,3 s 6 s + 1 1 + 0,3 s  

Colocando tudo na equação característica:

1 + G C G 1 = 0

1 + K C 1 - 0,3 s 6 s + 1 1 + 0,3 s = 0

Multiplicando tudo pelo denominador:

6 s + 1 1 + 0,3 s + K c ( 1 - 0,3 s ) = 0

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

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

Passo 3

Hora de aplicar nossa substituição:

s = i ω u

K C = K C u

1,8 i ω u 2 + ( 6,3 - 0,3 K C u ) i ω u + 1 + K C u = 0

Como:

i 2 = - 1

- 1,8 ω u 2 + ( 6,3 - 0,3 K C u ) i ω u + 1 + K C u = 0

Agora, já sabe né, devemos separar a parte real da imaginária e igualar as duas a 0 :

( 6,3 - 0,3 K C u ) i ω u = 0

- 1,8 ω u 2 + 1 + K C u = 0

Da parte imaginária temos:

6,3 - 0,3 K C u i ω u = 0

6,3 - 0,3 K C u = 0

K C u = 6,3 0,3 = 21 % S C % S T

Da parte real:

- 1,8 ω u 2 + 1 + K C u = 0

- 1,8 ω u 2 + 1 + 21 = 0

ω u = 22 1,8 = 3,50   r a d / m i n

Para determinar o períofo final, devemos aplicar a fórmula:

T u = 2 π ω u

Logo:

T u = 2 π 3,50 = 1,80   m i n

Resposta

K C u = 21 % S C % S T

T u = 1,80   m i n

Exercício Resolvido #5

Para uma malha de realimentação representada pelo diagrama de blocos abaixo:

Determine o ganho e período finais para um controlador proporcional para a função de transferência de processo abaixo:

G 1 s = ( 0,3 s + 1 ) s + 1 2 5 s + 1

Passo 1

Como a malha que estamos trabalhando é do tipo feedback, a função de transferência equivalente da malha é:

C s R s = G C G 1 1 + G C G 1

Para o método da substituição direta, o que tilizamos é a equação característica da malha, que nada mais é que o denominador da função de gransferência queivalente, ou seja, a nossa equação característica é:

1 + G C G 1 = 0

Passo 2

O nosso controlador é do tipo proporcional, portanto sua função de transferência será:

G C = K C

E a do processo:

G 1 s = ( 0,3 s + 1 ) s + 1 2 5 s + 1

Colocando tudo na equação característica:

1 + G C G 1 = 0

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

Multiplicando tudo pelo denominador:

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

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

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

Passo 3

Aplicando a nossa famosa substituição:

s = i ω u

K C = K C u

5 i ω u 3 + 11 i ω u 2 + ( 7 + 0,3 K C u ) i ω u + 1 + K C u = 0

Como:

i 2 = - 1

i 3 = - i

- 5 i ω u 3 - 11 ω u 2 + ( 7 + 0,3 K C u ) i ω u + 1 + K C u = 0

Separando a parte real da imaginária e igualando as duas a 0 :

- 5 i ω u 3 + ( 7 + 0,3 K C u ) i ω u = 0

- 11 ω u 2 + 1 + K C u = 0

Da parte imaginária temos:

- 5 i ω u 3 + ( 7 + 0,3 K C u ) i ω u = 0

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

ω u = 7 + 0,3 K C u 5

Substituindo na expressão da parte real:

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

- 11 7 + 0,3 K C u 5 + 1 + K C u = 0

- 11 1,4 + 0,06 K C u + 1 + K C u = 0

- 15,4 - 0,66 K C u + 1 + K C u = 0

0,34 K C u = 14,4

K C u 42,35   % S C % S T

Substituindo o valor de K C u :

ω u = 7 + 0,3 K C u 5

ω u = 7 + 12,7 5 = 1,98   r a d / m i n

Para determinar o períofo final, devemos aplicar a fórmula:

T u = 2 π ω u

Logo:

T u = 2 π 1,83 = 3,16   m i n

Resposta

K C u = 42,35 % S C % S T

T u = 3,16   m i n

Exercício Resolvido #6

Para uma malha de realimentação representada pelo diagrama de blocos abaixo:

Determine o ganho e período finais sabendo que o controlador é do tipo proporcional e que:

G v = 1 2 s + 1

G p = 1 0,5 s + 1

H = 1 s + 1

Passo 1

Essa é aquela malha feedback mais clássica de todas, cuja função de trsnferência equivalente é a seguinte

C s R s = G C G v G p 1 + G C G v G p H

Se você não se lembra como chegar nessa expressão, dá um revisadinha na aula de álgebra de diagrama de blocos, essa deduação tem que estar na ponta da língua.

Porém, para fazer a substituição direta precisamos mesmo é da equação característica, que é o denominador da função de transferência

1 + G C G v G p H = 0

Passo 2

Agora devemos substituir as funções de transferência do controlador e do processo na nossa equação característica:

G C = K C

G v = 1 2 s + 1

G p = 1 0,5 s + 1

H = 1 s + 1

Substituindo na equação característica:

1 + G C G v G p H = 0

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

Multiplicando tudo pelo denominador:

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

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

Passo 3

Finalmente então, podemos aplicar a nossa substituição para encontrar o ganho e período finais:

s = i ω u

K C = K C u

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

Como:

i 2 = - 1

i 3 = - i

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

Separando a parte real da imaginária e igualando as duas a 0 :

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

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

Da parte imaginária temos:

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

- ω u 2 + 3,5 = 0

ω u = 3,5 = 1,87   r a d / m i n

Substituindo na parte real:

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

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

K C u = 11,25   % S C % S T

Para determinar o períofo final, devemos aplicar a fórmula:

T u = 2 π ω u

Logo:

T u = 2 π 1,87 = 3,36   m i n

Resposta

K C u = 11,25 % S C % S T

T u = 3,36   m i n

Exercício Resolvido #7

Para uma malha de realimentação representada pelo diagrama de blocos abaixo:

Determine o ganho e período finais sabendo que o controlador é do tipo proporcional e que:

G v = 1 2 s + 1

G p = ( 0,1 s + 1 ) ( 0,5 s + 1 ) ( 1,5 s + 1 )

H = 2

Passo 1

Nesse exercício temos uma malha feedback com controlador, válvula, processo e transmissor. Essa malha tem a seguinte função de transferência equivalente:

C s R s = G C G v G p 1 + G C G v G p H

Se você não se lembra como chegar nessa expressão, dá um revisadinha na aula de álgebra de diagrama de blocos, essa deduação tem que estar na ponta da língua.

Sendo assim, nossa equação característica será:

1 + G C G v G p H = 0

Passo 2

Agora devemos substituir as funções de transferência do controlador e do processo na nossa equação característica:

G C = K C

G v = 1 2 s + 1

G p = ( 0,1 s + 1 ) ( 0,5 s + 1 ) ( 1,5 s + 1 )

H = 2

Substituindo na equação característica:

1 + G C G v G p H = 0

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

Multiplicando tudo pelo denominador:

2 s + 1 0,5 s + 1 1,5 s + 1 + 2 K c ( 0,1 s + 1 ) = 0

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

1,5 s 3 + 4,75 s 2 + ( 4 + 0 , 2 K C ) s + 1 + 2 K C = 0

Passo 3

Aplicando aquela substituição marota:

s = i ω u

K C = K C u

1,5 i ω u 3 + 4,7 5 i ω u 2 + ( 4 + 0 , 2 K C u ) i ω u + 1 + 2 K C u = 0

Como:

i 2 = - 1

i 3 = - i

- 1,5 i ω u 3 - 4,7 5 ω u 2 + ( 4 + 0 , 2 K C u ) i ω u + 1 + 2 K C u = 0

Separando a parte real da imaginária e igualando as duas a 0 :

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

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

Da parte imaginária temos:

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

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

ω u = 4 + 0 , 2 K C u 1,5

Substituindo na parte real:

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

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

Multiplicando toda equação por 1,5 :

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

- 19 - 0,95 K C u + 1,5 + 3 K C u = 0

- 17,5 + 2,05 K C u = 0

K C u = 17,5 2,05

K C u = 8,54   % S C % S T

Substituindo na expressão da frequência:

ω u = 4 + 0,2 K C u 1,5 = 4 + 1,71 1,5 = 1,95

Para determinar o períofo final, devemos aplicar a fórmula:

T u = 2 π ω u

Logo:

T u = 2 π 1,95 = 3,22   m i n

Resposta

K C u = 8,54 % S C % S T

T u = 3,22   m i n

Exercício Resolvido #8

Para uma malha de realimentação representada pelo diagrama de blocos abaixo:

Determine o ganho e período finais para um controlador proporcional para a função de transferência de processo abaixo:

G 1 s = e - s s + 1 2

Passo 1

Como a malha que estamos trabalhando é do tipo feedback, a função de transferência equivalente da malha é:

C s R s = G C G 1 1 + G C G 1

Mas nós estamos interessados na equação característica, que nada mais é do que o denominador da função de transferência. Sendo assim, nossa equação de transferência para esse caso será:

1 + G C G 1 = 0

Passo 2

Como estamos trabalhando com um controlador proporcional, a função de transferência para o controlador será:

G C = K C

Só que dessa vez apareceu um termo exponencial que indica um tempo morto, ou seja temos que pedir ajuda para o nosso amigo Padé e sua aproximação:

e - t 0 s ≅ 1 - t 0 2 s 1 + t 0 2 s

Para o nosso exercício então:

e - 0,6 s ≅ 1 - 1 2 s 1 + 1 2 s = 1 - 0 , 5 s 1 + 0,5 s

A nossa função de transferêncua para o processo será então:

G 1 s = 1 - 0,5 s s + 1 2 ( 1 + 0,5 s )  

Colocando tudo na equação característica:

1 + G C G 1 = 0

1 + K C 1 - 0,5 s s + 1 2 ( 1 + 0,5 s ) = 0

Multiplicando tudo pelo denominador:

s + 1 2 ( 1 + 0,5 s ) + K c ( 1 - 0 , 5 s ) = 0

s 2 + 2 s + 1 ( 1 + 0,5 s ) + K C - 0 , 5 K C s = 0

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

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

Passo 3

Fazendo aquela substituição esperta:

s = i ω u

K C = K C u

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

Mas:

i 2 = - 1

i 3 = - i

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

Agora, já sabe né, devemos separar a parte real da imaginária e igualar as duas a 0 :

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

- 2 ω u 2 + 1 + K C u = 0

Da parte imaginária temos:

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

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

ω u = 2,5 - 0,5 K C u 0,5 = 5 - K C u

Substituindo na parte real:

- 2 ω u 2 + 1 + K C u = 0

- 2 5 - K C u 2 + 1 + K C u = 0

- 2 5 - K C u + 1 + K C u = 0

- 10 + 2 K C u + 1 + K C u = 0

3 K C u - 9 = 0

K C u = 9 3 = 3 % S C % S T

Substituindo na expressão da frequência:

ω u = 5 - K C u

ω u = 5 - 3 = 2 = 1,41   r a d / m i n

Para determinar o períofo final, devemos aplicar a fórmula:

T u = 2 π ω u

Logo:

T u = 2 π 1,41 = 4,44   m i n

Resposta

K C u = 3 % S C % S T

T u = 4,44   m i n