Laboratorio de Sistemas de Controle 2

Disciplina Laboratório de Controle Discreto Trabalho em Equipe Considere o seguinte sistema analógico Este sistema deve ser controlado através de realimentação unitária por um controlador digital Dz com T01s e um segurador de ordem zero ambos em cascata com o sistema ver figura abaixo Os valores do ganho Kg e do polo pg da função de transferência Gps para cada equipe devem ser tirados da tabela no final deste texto Para todos os controladores calculados nos itens 2 a 5 devem ser traçados a resposta ao degrau unitário para a malha fechada com o controlador calculado e conferido as especificações desejadas em termos de tempo de subida sobrepasso e tempo de assentamento Também deve ser determinado os erros de regime 1 Discuta a estabilidade do sistema Gpz obtido a partir da avaliação de seus diagramas de Bode 1 ponto Mostre que o valor de ganho de estabilidade confere com o valor achado através do critério de Jury 2 Use a transformada bilateral w para projetar um controlador digital tal que a margem de fase seja de 60 com frequência de cruzamento de ganho maior que 1rads o erro estacionário seja reduzido em 10 vezes Justifique a escolha do controlador 2 pontos 3 Projete um controlador por alocação de polos que tenha desempenho de resposta transitória semelhante ao controlador do item 1 e erro ao degrau nulo Justifique seu projeto usando o diagrama de lugar das raízes 2 pontos 4 Projete um controlador do tipo PID que tenha desempenho de resposta transitória semelhante ao controlador do item 1 Quanto vale o erro de regime neste caso Justifique sua escolha de projeto 2 pontos 5 Projete um controlador em espaço de estado com desempenho semelhante ao controlador do item 1 Você pode usar como referência se necessário o par de polos encontrados no item 4 2 pontos 6 Sabemos que na prática o controlador discreto deve controlar o sistema contínuo Para comprovar o desempenho das malhas de controle projetados nos itens 13 4 e 5 monte os esquemas de controle usando o simulink Compare a resposta temporal para esses quatros esquemas de controle comente as semelhanças e eventuais diferenças Inclua os arquivos mxl desenvolvidos 1 ponto FORMATO DE ENTREGA Relatório um por equipe contendo todos os cálculos realizados para a solução das questões Matlab pode ser utilizado como apoio incluir comandos utilizados e figurasgráficos gerados Resposta apenas utilizando comandos do Matlab não serão aceitas Não é aceita resposta utilizando a toolbox SISOTOOLS nem PIDTUNE ALUNOS Ganho Kg Polo pg Cesar Batista 5 4 Na questao 2 o controlador que voces conseguiram não atinge o especificado o controlador acima não leva aos valores que voce indicam o que é comprovado pelo Bode que voce mostram logo a seguir Na questao 3 é pedido Projete um controlador por alocação de pólos que tenha desempenho de resposta transitória semelhante ao controlador do item 1 e erro ao degrau nulo Justifique seu projeto usando o diagrama de lugar das raízes 2 pontos A alocação de pólos deve ser feita no LUGAR DAS RAÍZES e não utilizando espaço de estado como vocês fizeram A questao de espaço de estado é a questao 5 Na questao 5 faltou a resposta ao degrau para atestar o desempenho Vocês não fizeram a questao 6 Análise em Tempo Discreto de Gps via ZOH Cesar Batista Matheus Diniz 8 de julho de 2025 Sumário 1 Exercício 1 2 11 Definições iniciais 2 12 Cálculo Analítico dos Coeficientes de Gpz 2 13 Estabilidade via Critério de Jury 3 14 Margem de Ganho pelo Diagrama de Bode 4 15 Lugar das Raízes 5 16 Resposta ao Degrau 6 17 Conclusões 7 2 Exercício 2 7 21 Especificações 7 22 Planta e ganho DC necessário 7 23 Escolha do compensador Lead 7 24 Resultados do Projeto 8 25 Análise do Diagrama de Bode Margem de Ganho e Margem de Fase 8 26 Análise da Resposta ao Degrau 9 27 Figuras 10 3 Exercício 3 10 31 Modelo Aumentado e Controlabilidade 11 32 Escolha dos Polos Desejados 11 33 Cálculo dos Ganhos 11 34 Desempenho em Malha Fechada 12 35 Justificativa pelo Lugar das Raízes 12 4 Exercício 4 12 401 Estrutura do PID Discreto 13 402 Especificações de Projeto 13 403 Cálculo Analítico do Ganho Integral 13 404 Resultados em Malha Aberta 14 405 Desempenho em Malha Fechada 14 406 Justificativa 14 1 Exercício 5 51 Modelo em Espaço de Estados 15 52 Escolha dos Polos Desejados 15 53 Cálculo de K e Nbar 16 54 Desempenho em Malha Fechada 16 55 Lugar das Raízes 16 1 Exercício 1 11 Definições iniciais Gps YsRs Kgs1spg Kg 5 pg 4 Ts 01 s 12 Cálculo Analítico dos Coeficientes de Gpz Fazendo frações parciais 5s1s4 As1 Bs4 A 53 B 53 gpt 53et 53e4t Em t kTs gk gpkTs 53e01k 53e04k Gpz sumk0 gkzk Os polos discretos são p1 e01 0904837 p2 e04 0670320 de modo que denz 1 p1 p2z1 p1 p2 z2 1 1575157 z1 0606530 z2 Em forma padrão Gpz b0 b1 z1 b2 z21 a1 z1 a2 z2 onde a1 1575157 a2 0606530 b0 0 b1 53p1 p2 0021157 b2 53p1 p2 53 p22 0017953 Em MATLAB num 0 00212 00180 den 1 15752 06065 13 Estabilidade via Critério de Jury Para realimentação unitária com ganho K o polinômio característico em z é 1 K Gpz 0 z2 α1 z α2 0 com α1 a1 K b1 α2 a2 K b2 onde numericamente a1 1575157 a2 0606530 b1 0021157 b2 0017953 As condições de estabilidade de 2ª ordem pelo critério de Jury são α2 1 1 1 α1 α2 0 2 1 α1 α2 0 3 Condição 1 α2 1 a2 K b2 1 1 0606530 0017953 K 1 4 0606530 0017953 K 1 0017953 K 0393470 K 03934700017953 21861 0606530 0017953 K 1 0017953 K 1606530 K 16065300017953 8950 Como K 0 a condição mais restritiva é K 21861 Condição 2 1 α1 α2 0 1 a1 K b1 a2 K b2 0 6 1 1575157 0021157 K 0606530 0017953 K 7 1 1575157 0606530 0021157 0017953 K 8 0031373 0039110 K 0 9 K 00313730039110 0802 10 Essa desigualdade é satisfeita para qualquer K 0 Condição 3 1 α1 α2 0 1 a1 Kb1 a2 Kb2 0 11 1 1575157 0021157 K 0606530 0017953 K 12 1 1575157 0606530 0021157 0017953 K 13 1968627 0003204 K 0 14 K 1968627 0003204 61470 15 Conclusão Das três condições de Jury a mais restritiva é K 21861 Portanto o ganho máximo que mantém o sistema estável em malha fechada é KmáxJury 21861 14 Margem de Ganho pelo Diagrama de Bode A margem de ganho Gm em malha aberta feedback unitário satisfaz GpejωTs Gm 1 em Gp 180 4 Figura 1 Mapa póloszeros do laço fechado Pz 0 Lendo o gráfico Gm 219017 20 log10 Gm 2681 dB 15 Lugar das Raízes O sistema contínuo Gps 5 s 1s 4 é discretizado por ZOH com Ts 01 s resultando em Gpz 0021157 z1 0017953 z2 1 1575157 z1 0606530 z2 Sejam p1 e01 09048 p2 e04 06703 os polos e z0 0848 o zero finito Assíntotas Número de assíntotas p z 1 ângulo ϕ 2k 1π 1 π centro em σa p1 p2 z0 1 24231 Ponto de breakaway Da condição Dzb 0 zb p1p2 2 07876 5 Conclusão Dois ramos partem de p1 e p2 Um ramo vai ao zero em z0 0848 O outro segue para o infinito ao longo do eixo real negativo A região de estabilidade em malha fechada é 0 K Kmax 2190 16 Resposta ao Degrau Figura 2 Resposta ao degrau de malha fechada 6 17 Conclusões Pelo Jury Kmax 21861 Pelo Bode Kmax 21902 268dB A coerência entre os métodos é melhor que 02 2 Exercício 2 Neste item projetamos um controlador digital Dz satisfazendo 1 Margem de fase desejada PMdes 60 2 Frequência de cruzamento de ganho ωcdes 1 rads 3 Redução do erro estacionário em degrau por fator 10 21 Especificações Margem de fase PMdes 60 Crossover de ganho ωcdes 1 rads Erro estacionário ess0 essd ess010 22 Planta e ganho DC necessário A planta discretizada por ZOH com Ts 01 s é Gpz b1 z1 b2 z21 a1 z1 a2 z2 Gp1 b1 b21 a1 a2 1252 Erro estacionário atual ess0 11 Gp1 0444 Para essd 00444 L0des 1 essdessd 2150 Kc L0des Gp1 1720 23 Escolha do compensador Lead Adotamos um compensador em avanço no domínio w Tustin Dw Kc 1 τ w1 α τ w com Kc 1720 Para PMdes 60 em ωc 12 φfalt 60 180 Gpj12 φmax φfalt 5 α 1 sin φmax1 sin φmax τ 1ωc α Convertendo via Tustin obtemos Dz n0 n1z1 1 d1z1 24 Resultados do Projeto Kc 1720 α 017 τ 180 PMobtido 5514 ωcobtido 130 rads ess 00444 25 Análise do Diagrama de Bode Margem de Ganho e Margem de Fase mostra as curvas de ganho acima e fase abaixo de Lz Dz Gpz Para verificar se as especificações foram atendidas identificamos Margem de Ganho GM valor de aumento de ganho antes de a fase atingir 180 GM 1200 dB Margem de Fase PM ângulo adicional de fase que o sistema ainda suporta quando Lejω 1 PMobtido 5514 próximo ao alvo de 60 8 Frequência de Cruzamento de Ganho ωc frequência em que Lejωc 1 ωcobtido 130 rads satisfazendo ωcdes 1 rads Interpretação 1 A curva de ganho cruza 0dB em ωc 2 Na curva de fase lêse L em ωc para calcular a PM 3 Uma maior folga de fase implica maior margem de robustez a variações de parâme tros 26 Análise da Resposta ao Degrau Tempo de Subida tr tr 130 s Sobressinal Mp Mp 028 Tempo de Assentamento ts ts 350 s Erro em Regime Permanente ess ess 00444 Observações O valor de tr e Mp reflete a margem de fase obtida maior PM geralmente reduz Mp e aumenta ligeiramente tr ts depende do amortecimento dominante com ζ 088 ts ficou em torno de 35s O erro estacionário confirma a redução de 10 em relação ao sistema sem compen sação 9 27 Figuras Figura 3 Diagrama de Bode de Lz DzGpz margens de ganho e fase Figura 4 Resposta ao degrau em malha fechada com Dz 3 Exercício 3 Objetivos Desempenho transitório similar ao controlador lead do Item 2 ts 35 s e Mp 03 10 Erro ao degrau nulo integrador interno 31 Modelo Aumentado e Controlabilidade A planta discreta Gpz em forma canônica controlável A a1 a2 1 0 B b0 b1 C 1 0 com valores numéricos a1 15752 a2 06065 b0 0 b1 00212 Estendemos com estado integrador xik 1 xik rk yk Ai 1 C 0 A Bc 0 B Br 1 0 0 Verificouse rank CAi Bc 3 sistema completamente controlável 32 Escolha dos Polos Desejados Para ts 35 s e Mp 03 adotamos amortecimento ζ 09 ωn 4 ζ ts 4 09 35 127 rads Os polos dominantes no contínuo são pc12 ζ ωn j ωn1 ζ2 114 j 058 Mapeando para z com Ts 01 s p12 epc12 Ts 089 j 006 Escolhemos o polo do integrador em p3 030 distante dos dominantes para não degradar o transitório 33 Cálculo dos Ganhos Alocando os polos p1 p2 p3 com realimentação de estados mais integrador Kaug Ki Kx1 Kx2 07823 326203 232551 Ou seja Ki 07823 Kx 326203 232551 34 Desempenho em Malha Fechada A resposta ao degrau unitário fornece tr 1090 228 s Mp 015 ts 2 390 s ess 0 degrau unitário Esses valores são muito próximos às metas ts 35 s e Mp 03 e o erro estacionário é nulo conforme desejado 35 Justificativa pelo Lugar das Raízes No lugar das raízes do sistema aumentado Figura 5 vemos que Dois polos complexos foram posicionados em 089j006 garantindo amortecimento elevado ζ 09 O polo real em 030 corresponde ao integrador interno deslocando o ganho DC para infinito e anulando o erro em regime Figura 5 Lugar das raízes do sistema aumentado Ai Bc com polos em 089 j006 030 4 Exercício 4 Neste item projetamos um controlador PID discreto que atenda simultaneamente a 1 Desempenho transitório semelhante ao controlador lead do Item 2 ts 35 s Mp 03 2 Erro de regime nulo para entrada degrau ess 0 12 401 Estrutura do PID Discreto Escolhemos a forma em z baseada em aproximação Tustin para a derivada e implementação exata para a integral Dz KP KI Ts 1 z1 KD 1 z1 Ts Aqui KP atua como ganho proporcional ajustando o ganho global do laço KI fornece ação integral pura garantindo ess 0 pois limz1 Dz KI Ts 1 z1 KD adiciona ação derivativa aumentando a margem de fase e melhorando a resposta transitória 402 Especificações de Projeto Para que o PID gere comportamento similar ao controlador lead projetado Frequência de cruzamento de ganho ωcdes 12 rads pois em controle digital ts 4 ζ ωc e com ζ 088 temos ts 35 s Margem de fase mínima PMdes 60 assegurando robustez similar Erro estacionário ess 0 para degrau unitário Período de amostragem Ts 01 s 403 Cálculo Analítico do Ganho Integral Para entrada degrau Rz 1 1 z1 a função de malha aberta é Lz Dz Gpz e o erro em regime para degrau é ess limz1 1 1 Lz limz1 1 z1 1 z1 1 z1 Dz Gpz Como 1 z1 Dzz1 KI Ts para ess 0 basta KI 0 O valor exato de KI será ajustado no passo de sintonia automática 404 Resultados em Malha Aberta O diagrama de Bode de Lz Dz Gpz Fig 6 mostrou PMobtido 7411 ωcobtido 120 rads Figura 6 Diagrama de Bode de Lz PIDz Gpz margens de ganho e fase 405 Desempenho em Malha Fechada Após aplicação de Dz a resposta ao degrau apresentou tr 1090 139 s Mp 475 ts 2 417 s ess 0 Embora ts e ωc estejam próximos das metas o sobressinal elevado 475 mostra que o PID é menos amortecido que o lead projetado no Item 2 406 Justificativa A ação proporcional e derivativa elevou a margem de fase para 74 garantindo robustez adicional A ação integral assegurou ess 0 O alto sobressinal pode ser reduzido ajustando KD para aumentar o amortecimento 14 5 Exercício 5 Objetivo projetar realimentação de estados discreta com précompensador Nbar de forma que o desempenho transitório seja semelhante ao controlador lead Item 2 ts 35 s e Mp 03 o erro em regime seja nulo ess 0 51 Modelo em Espaço de Estados Começamos pela representação em espaço de estados da planta discretizada via segurador de ordem zero ZOH com período de amostragem Ts 01 s A função de transferência discreta é Gpz b0 b1 z1 b2 z2 1 a1 z1 a2 z2 cujo modelo na forma canônica controlável é dado por xk 1 A x k B uk yk C xk com A a1 a2 1 0 B b0 b1 C 1 0 Substituindo os coeficientes obtidos numericamente a1 15752 a2 06065 b0 0 b1 00212 Dessa forma A 15752 06065 1 0 B 0 00212 C 1 0 Antes de prosseguir verificamos a controlabilidade do par A B calculando C B A B cujo posto é 2 igual à dimensão de A confirmando que todo estado pode ser controlado 52 Escolha dos Polos Desejados Para reproduzir no domínio discreto o mesmo desempenho transitório do controlador em avanço Item 2 que apresentava tempo de assentamento ts 35 s e sobressinal Mp 03 adotamos 1 Coeficiente de amortecimento aproximado ζ 09 2 Cálculo da frequência natural ωn 4 ζ ts 4 09 35 127 rads 3 Polos desejados no contínuo pc12 ζ ωn j ωn1 ζ2 114 j 058 4 Mapeamento para o domínio z via z esTs p12 epc12Ts 089 j 006 Desse modo requeremos que o polinômio característico do sistema em malha fechada tenha raízes exatamente em p1 e p2 53 Cálculo de K e Nbar Ganho de realimentação de estados K Aplicamos o comando K placeA B p1 p2 que retorna K K1 K2 103631 96453 Isso garante que a malha fechada Acl A B K possui exatamente os polos desejados Précompensador Nbar Para assegurar erro zero em regime permanente para entrada degrau unitário introduzimos um ganho de entrada Nbar A fórmula padrão é Nbar C I A B K1 B1 que numericamente resulta em Nbar 12188 Esse valor corrige o ganho do caminho de referência de forma que limk yk rk para rk 1 54 Desempenho em Malha Fechada Resposta ao degrau unitário tr 1090 205 s Mp 048 ts 2 329 s ess 0 Parâmetro Valor obtido tr 1090 205 s Mp 048 ts 2 329 s ess 0 55 Lugar das Raízes Na figura vemos o lugar das raízes de A BK com polos em 089 j006 que garantem amortecimento elevado ζ 09 e bos desempenho transitório Figura 7 Lugar das raízes de A BK com polos em 089 j006 17