Exercícios Resolvidos de Funções de Transferência no Simulink

Ver Teoria

Enunciado

Considere o motor DC abaixo:

Ele tem como entrada a tensão v a e como saída, o ângulo θ . Considerando as constantes:

  • Momento de inércia do motor e carga: J = 3,2284 ∙ 10 - 6   k g ∙ m ² ;
  • Coeficiente de atrito do motor e carga: B = 3,5077 ∙ 10 - 6   k g ∙ m 2 / s ;
  • Resistência de armadura: R a = 4   Ω ;
  • Indutância da armadura: L a = 2,75 ∙ 10 - 6   H ;
  • Constante elétrica: K e = 0,0274   V s ;
  • Constante de torque: K t = 0,0274   N m / A .

Considerando as condições iniciais nulas, o motor DC tem como função de transferência:

H s = θ ( s ) v a ( s ) = K t J L a s 3 + J R a + B L a s 2 + ( B R a + K e K t ) s

Usando o SIMULINK, simule a resposta ao degrau em um sistema em malha fechada.

Passo 1

O enunciado nos fornece a função de transferência de um motor DC e pede para construirmos o sistema em malha fechada no SIMULINK. Em seguida, simularemos a resposta ao degrau do sistema em malha fechada que construímos.

Nosso modelo no SIMULINK terá quatro blocos:

  • A entrada do modelo será um bloco “Step”, simulando um degrau unitário;
  • Usaremos um bloco “Sum” para construirmos o sistema em malha fechada;
  • A função de transferência H ( s ) será incluída no modelo por meio do bloco “Transfer Fcn”;
  • A saída do modelo será um bloco “Scope”, que mostrará a resposta do sistema ao degrau unitário.

Passo 2

Em novo modelo do SIMULINK, insira um bloco “Step”, que servirá como entrada no nosso modelo. Como o enunciado não deu maiores especificações sobre o degrau, não vamos alterar as propriedades de “Step”.

Em seguida, insira um bloco “Sum”, para começarmos a construir o sistema em malha fechada. Você encontrará esse bloco na biblioteca “Math Operations”:

O bloco “Sum” tem duas entradas e uma saída. Por padrão, cada entrada tem um sinal de + , indicando que o “Sum” vai somar dois sinais de entrada e a soma será fornecida pela saída do bloco. Contudo, em uma malha fechada, fazemos a subtração entre um sinal de referência (no caso, nosso bloco “Step”) e o resultado do produto entre esse sinal de referência e uma função de transferência (o sinal que sairá do bloco “Transfer Fcn” representando H ( s ) , como veremos mais pra frente).

Logo, temos que alterar os parâmetros do bloco. Dê um duplo-clique sobre “Sum” para acessar seus parâmetros. A propriedade que vamos alterar é a “List of signs”, destacada em vermelho abaixo:

Como o padrão do “Sum” é somar dois sinais, no campo “List of signs” temos | + + . Para que seja feita a subtração de sinais, troque por | + - . Clique em “OK” e você verá que o bloco vai ter esses sinais em seu interior:

Passo 3

O próximo passo é inserir H ( s ) no modelo. Vá à biblioteca “Continuous” e arraste um bloco “Transfer Fcn” para o modelo:

Para definir H ( s ) , dê um duplo-clique sobre o bloco “Transfer Fcn” que você acabou de incluir no modelo para ter acesso a seus parâmetros. Vamos mexer nos parâmetros “Numerator coefficients” (destacado em vermelho abaixo) e “Denominator coefficients” (destacado em azul).

Entramos vetores com coeficientes de polinômios, indo do elemento de maior expoente ao de menor expoente, em ambos os parâmetros. No “Numerator coefficients”, o vetor tem os coeficientes do polinômio do numerador de H ( s ) . Já o vetor “Denominator coefficients” contém os coeficientes do polinômio do denominador de H ( s ) .

O numerador de H ( s ) é K t = 0,0274   N m / A , logo, devemos trocar o [ 1 ] do campo destacado em vermelho por [ 0.0274 ] .

O polinômio do denominador é J L a s 3 + J R a + B L a s 2 + ( B R a + K e K t ) s . Considerando as constantes dadas no enunciado, temos:

J L a = 3,2284 ∙ 10 - 6 ∙ 2,75 ∙ 10 - 6 = 8,88 ∙ 10 - 12

J R a + B L a = 3,2284 ∙ 10 - 6 ∙ 4 + 3,5077 ∙ 10 - 6 ∙ 2,75 ∙ 10 - 6 = 1,29 ∙ 10 - 5

B R a + K e K t = 3,5077 ∙ 10 - 6 ∙ 4 + 0,0274 ∙ 0,0274 = 7,65 ∙ 10 - 4

Logo, o vetor a ser entrado no campo destacado em azul é:

8,88 ∙ 10 - 12   1,29 ∙ 10 - 5   7,65 ∙ 10 - 4   0    

Portanto, os parâmetros “Numerator coefficients” e “Denominator coefficients” de H ( s ) ficarão assim:

Clique em “OK” para registrar as mudanças e fechar automaticamente a janela com os parâmetros. Você notará que o bloco “Transfer Fcn” tem agora H ( s ) em seu interior. Aumente o tamanho do bloco, selecionando-o e arrastando um dos seus cantos. Dessa forma, você conseguirá ver H ( s ) por completo. Por enquanto, nosso modelo está assim:

Passo 4

Vamos inserir a saída do modelo: um bloco “Scope”, para vermos o gráfico da resposta ao degrau do motor DC. Após adicionar o bloco “Scope”, conecte os blocos, primeiramente como se fosse uma malha aberta:

Para fechar a malha, clique sobre a entrada de “Sum” que tem o sinal de menos e arraste. Aparecerá uma flecha vermelha tracejada. Conecte a extremidade “solta” dessa flecha na linha entre o bloco com H ( s ) e o “Scope”, para criar um ponto de derivação:

Pronto! Temos o sistema em malha fechada e já podemos executar o programa, para ver a resposta ao degrau.

Após encerrar a execução, clique duas vezes sobre o “Scope” para ver o gráfico. Use o botão “Autoscale”, destacado abaixo, para ajustar a escala do gráfico.

Vixe, ele ficou bem esquisito.

Passo 5

O motivo de o gráfico ter ficado esquisito é que o número de dados calculados na simulação ultrapassou o limite padrão do bloco “Scope”. Para corrigir esse problema, temos que acessar as propriedades do bloco “Scope”, clicando na engrenagem destacada abaixo (que está na barra superior da janela do “Scope”), e mudar a propriedade “Limit data points to last”.

Ao clicar na engrenagem, o “Scope parameters” se abrirá. Vá à aba “History”, onde tem a propriedade “Limit data points to last”, destacada abaixo:

Desmarque a caixinha à esquerda, dentro do destaque em vermelho. Dessa forma, não haverá limite de dados e o gráfico ficará bonitão.

Clique em “OK” para registar a alteração, feche a janela do “Scope” e rode a simulação novamente. Ao abrir o gráfico, você verá que ele terá um formato que estamos mais acostumados a ver em exercícios de controle com entrada degrau.

Resposta