O efeito da adição de zeros e polos na função de transferência

Ver Todos os Exercícios

Teoria

Introdução

Agora que nós já vimos como se comportam e como são as respostas dos sistemas de primeira e segunda ordem, vamos continuar nossa jornada e tentar entender os sistemas de ordem superior.

Uma função de transferência normal pode ser descrita como:

Y s = K ∏ j = 1 m τ j s + 1 ∏ k = 1 n τ k s + 1 X ( s )

Onde n ≥ m .

Nós já vimos as funções de primeira ordem, onde m = 0 e n = 1 , as de segunda ordem, com m = 0 e n = 2 . Nessa seção vamos analisar mais alguns casos além desses e ver como poderemos simplificá-los.

Só uma coisinha antes de começar: as raízes do polinômio do numerador são chamadas de zeros, enquanto as raízes do polinômio do denominador são chamadas de polos.

Sendo assim, as funções de transferência de primeira ordem têm um polo e nenhum zero, as funções de transferência de segunda ordem têm dois polos e nenhum zero.

Nosso objetivo aqui então vai ser entender qual será o efeito da adição de polos e zeros na nossa função de transferência. Tentaremos entender como eles impactam a resposta do processo e vamos tentar procurar casos para poder simplificar essas funções desconsiderando polos ou zeros.

Bora entrar nessa viagem?

Adição de um polo

Vamos considerar uma função de transferência de terceira ordem do tipo:

Y s X s = K τ 2 s 2 + 2 τ ζ s + 1 s + α

Considerando que 0 < ζ < 1 os polos da função de transferência serão:

- ζ τ ± ζ τ i 1 - ζ 2

E

- α

Ou seja, teremos dois polos imaginários referentes à equação do segundo grau do tipo subamortecida e um polo real – α .

Quando tivermos mais de dois polos, como é o caso aqui, será interessante para nós procurar por polos dominantes.

“Mas o que raios são polos dominantes?

Calma, já te explico, mas antes vamos continuar analisando a nossa função de transferência.

Considerando uma entrada degrau unitário, teremos:

Y ( s ) = K τ 2 s 2 + 2 τ ζ s + 1 s + α 1 s  

A inversa dessa função, no domínio tempo, será algo do tipo:

y t = K A u t + e - ζ τ t B cos ⁡ ω d t + C sen ⁡ ω d t + D e - α t  

Os coeficientes A , B , C e D são resultantes da separação em frações parciais.

Não se assuste, você não precisa saber como fazer essa inversa, coloquei ela aqui só para a gente poder analisar o efeito dos polos na resposta final.

Vamos pensar juntos então. A adição do terceiro polo ( s + α ) resultou na nossa resposta na adição de um termo D e - α t . Esse por sua vez é um termo exponencial negativo, ou seja, quanto maior for o valor do expoente α t , menor será o valor do termo, e menos ele impactará na resposta final.

Sendo assim, quando o polo α é muito maior do que os outros dois polos, nós podemos desconsiderar o polo que está mais afastado, reduzindo a função de transferência em uma função de transferência normal de segunda ordem. Os outros dois polos, então, são os polos dominantes que eu havia falado agora há pouco.

Essa aproximação é muito boa, pois quando podemos fazê-la, nós conseguimos calcular todos os parâmetros como overshoot, tempo de subida, tempo de sedimentação da mesma forma que calculamos para uma função de segunda ordem comum.

”Ok, mas como eu vou saber que um polo é suficientemente grande para poder desconsiderar”

Boa pergunta!

Bom, isso vai depender do seu grau de exigência, se você precisa de uma aproximação mais exata, você vai precisar que esse polo seja muito grande para poder desconsiderá-lo, agora se seu grau de exigência não é muito grande, polos um pouco maiores que os outros já serão suficientes.

De modo geral, um valor seguro para fazer a aproximação é quando o maior polo for cinco vezes maior que os outros polos. Nesses casos a aproximação pode ser feita sem medo, visto que uma diferença de cinco vezes é bastante coisa. Mas como eu disse, se o seu grau de exigência for menor, a aproximação pode ser feita para valores menores.

Só para exemplificar melhor, a figura abaixo apresenta três respostas. y 1 ( t ) representa a resposta de uma função de transferência com apenas dois polos. y 2 ( t ) representa uma função com os mesmos dois polos de y 1 ( t ) e um terceiro polo muito maior que os outros dois. Por fim, a curva y 3 ( t ) representa uma função com os dois polos de y 1 ( t ) e um terceiro polo bem próximo dos outros dois.

Repare que como o terceiro polo de y 2 ( t ) é muito maior que os outros dois, a resposta é bem parecida com y 1 ( t ) . O mesmo não acontece com y 3 ( t ) cujo terceiro polo tem valor bem próximo dos outros dois.

Efeito do zero

Além da adição de polos à nossa função de transferência, também podemos ter casos em que são adicionados zeros. Por exemplo:

Y s X s = K s + α ( s + b ) ( s + c )

Os zeros, quando são adicionados em um sistema, afetam a amplitude, ou o overshoot, da resposta, mas não interferem em nada no tipo de resposta.

Portanto, se um sistema de segunda ordem é subamortecido sem nenhum zero, ao se adicionar um zero nesse sistema, ele continuará subamortecido.

Os zeros funcionam de maneira parecida com os polos. Quanto mais próximo for o zero dos polos, maior será o efeito sentido na resposta.

Quando os zeros se distanciam muito dos polos, seu efeito praticamente não é sentido e a função de transferência pode ser simplificada como uma função de transferência de segunda ordem comum.

Dessa forma, quando temos um zero próximo dos polos, a resposta transitória é levemente alterada, principalmente na amplitude e excesso (overshoot). Quando o zero se afasta muito dos polos, a resposta transitória é basicamente a mesma da função de transferência original, porém em uma versão em escala. Nesse caso, dizemos que o zero funciona como um fator de ganho.

Por exemplo, imagine uma função G s de segunda ordem.

G s = 1 s + b s + c

Se adicionarmos um zero a essa função:

s + α G s = s s + b s + c + α s + b s + c

Se o zero for muito maior que os polos, ou seja, se α for muito grande, a função de transferência pode ser simplificada como:

s + α G s ≅ α s + b s + c = α G ( s )

Portanto, quando o zero está muito longe dos polos, ele funciona como um termo multiplicativo da função de segunda ordem com dois polos, por isso que dizemos que ele funciona como um fator de ganho.

Existe ainda um outro caso em que α é negativo, ou seja, o zero é positivo.

Quando isso acontece, a resposta transitória tende a desviar valor final no início antes de se aproximar dele. O gráfico abaixo mostra esse comportamento:

Na engenharia elétrica esse tipo de sistema é chamado de sistema de fase não mínima. Os sistemas de fase não mínima sempre começam agindo contrário ao desejado antes de começarem a corrigir a variável controlada. Sendo assim, se um avião funcionasse como um sistema de fase não mínima, e você comandasse que ele virasse para a esquerda, primeiramente ele viraria para a direita antes de virar para a esquerda.

Cancelamento de polos com zeros

Além dos casos onde existe um polo ou um zero muito longe da origem, uma outra forma de fazer aproximações em uma função de transferência com múltiplos polos e zeros é cancelando um polo com um zero.

Esse cancelamento pode ser feito quando existe um zero bem próximo de um polo. Para avaliar se o cancelamento pode ser feito, devemos fazer o seguinte procedimento.

Vamos usar como exemplo a função de transferência hipotética abaixo:

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

Vamos supor que a gente considera que o zero em - a é bem próximo do polo em - d , e por isso queremos saber se é possível fazer o cancelamento. Para isso, devemos expandir a função de transferência em suas frações parciais:

G s = B s + b + C s + c + D s + d

Onde B , C e D são os coeficientes encontrados no processo de expansão em frações parciais. Para saber se o polo em - d pode ser cancelado com o zero em - a devemos avaliar o coeficiente referente a esse polo, ou seja, o coeficiente D .

Portanto, para cancelar o polo com o zero, o coeficiente D deve ser muito menor perto dos outros coeficientes ( B e C ). Se isso for verdade, podemos efetuar a simplificação e obteremos a função de transferência simplificada abaixo:

G s = K s + b s + c

Adição de um polo

Efeito do zero

Cancelamento de polos com zeros

Exercícios Resolvidos

Exercício Resolvido #1

Determine a validade de uma aproximação de segunda ordem para a função de transferência abaixo:

G s = 700 s + 15 s 2 + 4 s + 100

Passo 1

Para saber se podemos aproximar a função acima a uma função de transferência de segunda ordem, precisamos determinar o polos delas e ver se existe um que seja muito maior que os outros dois.

Um dos polos, é fácil de encontrar:

s = - 15

Passo 2

Vamos agora lembrar da nossa professora da oitava série e fazer Báskara para determinar os outrso dois polos:

s 2 + 4 s + 100 = 0

s = - 4 ± 4 2 - 4 * 1 * 100 2

s = - 4 ± - 384 2 = - 2 ± 9,79 i

Passo 3

Temos então um polo em - 15 e dois polos em - 2 ± 9,79 i . Para saber se podemos desconsiderar um desses polos, precisamos comparar a parte real desses polos.

Sendo assim, o polo em - 15 é mais de sete vezes maior que os outros dois polos (em módulo), o que faz com que esse polo possa ser desconsiderado.

Portanto, a aproximação é válida.

Resposta

Aproximação válida.

Exercício Resolvido #2

Determine a validade de uma aproximação de segunda ordem para a função de transferência abaixo:

G s = 360 s + 4 s 2 + 2 s + 90

Passo 1

Vamos determinar os polos da função acima para saber se podemos aproximar ela de uma função de segunda ordem.

O polo de ordem 1 é fácil determinar:

s + 4 = 0

s = - 4

Passo 2

Agora, vamos pegar o polo de segunda ordem e aplicar Báscara para determinar os outros dois polos da função de transferência.

s 2 + 2 s + 90 = 0

s = - 2 ± 2 2 - 4 * 1 * 90 2

s = - 2 ± - 356 2 = - 1 ± 9,43 i

Passo 3

Os nosso três polos são portanto:

s = - 4

s = - 1 + 9,43 i

s = - 1 - 9,43 i

Para saber se podemos aproximar precisamos ver quantas vezes o polo maior em módulo é maior que os outrs. Lembrando que para fazer essa comparação usamos sempre a parte real dos polos.

Sendo assim, podemos ver que o polo em - 4 é apenas quatro vezes maior que os outros dois polos, o que não produz uma aproximação muito boa para um modelo de segunda ordem.

Portanto, a aproximação não é válida.

Resposta

Aproximação não é válida.

Exercício Resolvido #3

Para a função de transferência dada abaixo, determine se há cancelamento entre o zero e o polo mais próximo.

G s = 26,25 ( s + 4 ) s s + 3,5 s + 5 ( s + 6 )

Passo 1

Vamo conferir se podemos cortar o zero em - 4 com o polo em - 3,5 (polo mais próximo do zero).

Para isso devemos fazer a expansão em frações parciais da função de transferência:

G s = A s + B s + 3,5 + C s + 5 + D s + 6

Passo 2

Vamos determinar os coeficientes A , B , C e D .

Se você não lembra da aula de cáluclo, vamos relembrar como fazemos isso. Para determinar um coeficiente a gente faz o limite da função antes da expansão sem o polo em questão tendendo ao valor do polo que foi retirado.

Assim ó:

A = lim s → 0 ⁡ 26,25 ( s + 4 ) s + 3,5 s + 5 ( s + 6 ) = 1

Ou seja, para determinar o A , a gente tirou o polo dele da função de transferência original, e fizemos que s fosse igual a esse polo.

Para os outros coeficientes:

B = lim s → - 3,5 ⁡ 26,25 ( s + 4 ) s s + 5 ( s + 6 ) = - 1

C = lim s → - 5 ⁡ 26,25 ( s + 4 ) s s + 3,5 ( s + 6 ) = - 3,5

D = lim s → - 6 ⁡ 26,25 ( s + 4 ) s s + 3,5 ( s + 5 ) = 3,5

Passo 3

Sendo assim, a expansão em frações parciais será:

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

Para avaliar se podemos cancelar o polo em - 3,5 precisamos ver o seu coeficiente nas frações parciais, ou seja, o coeficiente B .

Nos nossos cálculos o coeficiente B deu - 1 , o que é bem próximo dos outros coeficientes (precisamos avaliar o valor absoluto), portanto, o polo em - 3,5 não pode ser cacelado.

Resposta

Aproximação não é válida.

Exercício Resolvido #4

Para a função de transferência dada abaixo, determine se há cancelamento entre o zero e o polo mais próximo.

G s = 26,25 ( s + 4 ) s s + 4,01 s + 5 ( s + 6 )

Passo 1

O polo mais próximo do zero em - 4 é o polo em - 4,01 . Vamos fazer a expansão em frações parciais da função de transferência para avaliar se o cancelamento pode ser feito ou não.

G s = A s + B s + 3,5 + C s + 5 + D s + 6

Passo 2

Vamos determinar os coeficientes A , B , C e D .

Para determinar esses coeficientes, a gente pega a função de transferência e retira dela o polo que a gente deseja determinar o coeficiente. Então, basta fazer a variável s tender ao valor do polo que foi retirado.

A = lim s → 0 ⁡ 26,25 ( s + 4 ) s + 4,01 s + 5 ( s + 6 ) = 0,87

Ou seja, para determinar o A , a gente tirou o polo dele da função de transferência original, e fizemos que s fosse igual a esse polo.

Para os outros coeficientes:

B = lim s → - 4,01 ⁡ 26,25 ( s + 4 ) s s + 5 ( s + 6 ) = 0,033

C = lim s → - 5 ⁡ 26,25 ( s + 4 ) s s + 4,01 ( s + 6 ) = - 5,3

D = lim s → - 6 ⁡ 26,25 ( s + 4 ) s s + 4,01 ( s + 5 ) = 4,4

Passo 3

Sendo assim, a expansão em frações parciais será:

G s = 0,87 s - 0,033 s + 3,5 - 5,4 s + 5 + 4,4 s + 6

Repara que o coeficiente do polo - 4,01 foi muito pequeno, - 0,033 , que é menos de 5% do segundo menor coeficiente e por isso podemos fazer o cancelamento desse polo com o zero da função de transferência.

Sendo assim, a função de transferência aproximada seria:

G s = 26,25 s s + 5 s + 6

Resposta

Aproximação é válida.

G s = 26,25 s s + 5 s + 6

Exercício Resolvido #5

Determine a validade de uma aproximação de respost ao degrau de segunda ordem para a função de transferência abaixo:

G s = 185,71 ( s + 7 ) s + 6,5 s + 10 ( s + 20 )

Passo 1

Para avaliar se podemos aproximar a função de transferência acima a uma de segunda ordem comum, precisamos verificar se o polo mais próximo do zero em - 7 poderá ou não ser cortado.

Para isso, vamos fazer a expansão em frações parciais:

G s = A s + 6,5 + B s + 10 + C s + 20

Passo 2

Vamos determinar os coeficientes A , B e C .

Usando o método de cacular os coeficientes da expansão em frações parciais:

A = lim s → - 6,5 ⁡ 185,71 ( s + 7 ) s + 10 ( s + 20 ) = 1,96

B = lim s → - 10 ⁡ 185,71 ( s + 7 ) s + 6,5 ( s + 20 ) = 15,91

C = lim s → - 20 ⁡ 185,71 ( s + 7 ) s + 6,5 s + 10 = - 17,88

Passo 3

Sendo assim, a expansão em frações parciais será:

G s = 1,96 s + 6,5 + 15,91 s + 10 - 17,88 s + 20

Como podemos ver, o coeficiente do polo em - 6,5 , que é o polo mais próximo do zero, é cerca de 10% do valor dos outros polos. Apesar de ele ser bem pequeno, a aproximação para uma função de transferência ainda não é recoemendada.

Portanto, a aproximação não é válida.

Resposta

Aproximação não é válida.

Exercício Resolvido #6

Determine a validade de uma aproximação de respost ao degrau de segunda ordem para a função de transferência abaixo:

G s = 197,14 ( s + 7 ) s + 6,9 s + 10 ( s + 20 )

Passo 1

O polo mais próximo do zero aqui é o polo em - 6,9 , então precisamos decidir se ele é próximo o suficiente para fazer o cancelamento.

Para isso, vamos fazer a expansão em frações parciais:

G s = A s + 6 , 9 + B s + 10 + C s + 20

Passo 2

Vamos determinar os coeficientes A , B e C .

Usando o método de cacular os coeficientes da expansão em frações parciais:

A = lim s → - 6,9 ⁡ 197,14 ( s + 7 ) s + 10 ( s + 20 ) = 0,49

B = lim s → - 10 ⁡ 197,14 ( s + 7 ) s + 6,9 ( s + 20 ) = 19,08

C = lim s → - 20 ⁡ 1 97,14 ( s + 7 ) s + 6,9 s + 10 = - 19,56

Passo 3

Sendo assim, a expansão em frações parciais será:

G s = 0,49 s + 6 , 9 + 19,08 s + 10 - 19,56 s + 20

A expansão em frações parciais nos mostrou que o coeficiente do polo mais próximo do zero, o polo em - 6,9 , é muito pequeno perto dos outros dois coeficientes, portanto, é possível fazer o cancelamento do polo com o zero.

Portanto, a aproximação é válida.

G s = 197,14 s + 10 ( s + 20 )

Resposta

Aproximação é válida.

G s = 197,14 s + 10 ( s + 20 )

Exercício Resolvido #7

Avalie se é possível aproximar a função de transferência abaixo a uma função de transferência de segunda ordem. Se for possível, determine o overshoot e o tempo de assentamento (estabilização) para uma faixa de 5 % .

G s = 7,5 s + 20 ( s 2 + 0,6 s + 1 )

Passo 1

Primeiramente precisamos checar se é possível ou não fazer a aproximação. Como temos três polos, precisamos que um dos polos esteja muito mais afastado da origem que os outros dois. Os três polos são:

s = - 0,3 + 1,9 i

s = - 0 , 3 - 1,9 i

s = - 20

Como podemos ver o polo em - 20 está muito mais afastado que os outros dois. Ele é mais de seiscentas vezes maior que a parte real dos outros dois polos, o que faz com que ele seja uma ótima opção de cancelamento.

Sendo assim, a proximação é válida, e a função aproximada será:

G s = 7,5 s 2 + 0,6 s + 1

Passo 2

Precisamos então calcular o overshoot e o tempo de assentamento. Para isso precisamos determinar a constante de tempo, τ , e o fator de amortecimento, ζ da função de segunda ordem aproximada.

G s = 7,5 s 2 + 0,6 s + 6

Precisamos comparar essa função de transferência com a função de transferência padrão de segunda ordem:

G s = K τ 2 + 2 ζ τ s + 1

Sendo assim:

K = 7,5

τ = 1

ζ = 0,3

Passo 3

O cálculo do overshoot é dado pela fórmula abaixo:

O S =   e - π ζ 1 - ζ 2

Usando ζ = 0,3 teremos:

O S = 0,37

E tempo de assentamento para uma faixa de 5 % pode ser calculado pela fórmula abaixo:

t s = 3 τ ζ

Com τ = 1 e ζ = 0,3 :

t s = 10

Resposta

Aproximação é válida.

G s = 7,5 s 2 + 0,6 s + 6

O S = 0,37

t s = 10

Exercício Resolvido #8

Avalie se é possível aproximar a função de transferência abaixo a uma função de transferência de segunda ordem. Se for possível, determine tipo de amortecimento e o tempo de assentamento (estabilização) para uma faixa de 5 % para a função de segunda ordem aproximada.

G s = 25 ( s + 5 ) s + 4,9 ( s + 2 ) ( s + 3 )

Passo 1

Primeiramente precisamos checar se é possível ou não fazer a aproximação. Nesse casos temos um polo muito perto do nosso zero, então vamos expandir a função de transferência em suas frações parciais:

G s = A s + 4,9 + B s + 2 + C s + 3

Usando o método de cálculo dos coeficientes das frações parciais:

A = lim s → - 4,9 ⁡ 25 ( s + 5 ) ( s + 2 ) ( s + 3 ) = 0,45

B = lim s → - 2 ⁡ 25 s + 5 s + 4,9 s + 3 = 25,86

C = lim s → - 3 ⁡ 25 ( s + 5 ) ( s + 4,9 ) ( s + 2 ) = - 26,31

Podemos ver que o coeficiente A , que é relativo ao polo mais perto do zero, - 4,9 , é muito menor que os outros dois coeficientes. Portanto, a aproximação de cancelamento polo-zero pode ser feita.

Sendo assim, a proximação é válida, e a função aproximada será:

G s = 25 ( s + 5 ) ( s + 2 ) ( s + 3 )

Passo 2

Precisamos determinar o tipo de amortecimento e o tempo de assentamento da resposta aproximada. Para isso precisamos do tempo de amortecimento, ζ , e contstante de tempo, τ . Para isso, vamos primeiro expandir a nossa função de transferência:

G s = 25 s 2 + 5 s + 6

Precisamos comparar essa função de transferência com a função de transferência padrão de segunda ordem:

G s = K τ 2 s 2 + 2 ζ τ s + 1

Então precisamos fazer umas manipulações:

G s = 25 1 6 s 2 + 5 6 s + 1

G s = 150 s 2 + 5 6 s + 1

Sendo assim:

K = 150

τ = 1

ζ = 0,42

Passo 3

Para determinar o tipo de assentamento, podemos usar o valor de ζ e a figura abaixo:

Como o nosso fator de amortecimento está entre 0 e 1, a nossa resposta será subamorticida.

O tempo de assentamento para uma faixa de 5 % pode ser calculado pela fórmula abaixo:

t s = 3 τ ζ

Com τ = 1 e ζ = 2,5 :

t s = 7,2

Resposta

Aproximação é válida.

G s = 25 s 2 + 5 s + 6

Resposta subamortecida.

t s = 7 , 2