PSC_nomedoaluno_tipo_turma - pucrs
PSC_nomedoaluno_tipo_turma - pucrs
PSC_nomedoaluno_tipo_turma - pucrs
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Trabalho – Projeto de compensadores de avanço-atraso de fase.<br />
Método da Resposta em Freqüência<br />
Número Identificador<br />
Técnica a ser utilizada<br />
Data Limite (Ler as instruções deste item para enviar o trabalho)<br />
Roteiro<br />
Anexo<br />
I. Método da Resposta em Freqüência<br />
O método consiste na avaliação da resposta em freqüência do processo de forma a<br />
introduzir um compensador de avanço-atraso de fase que atenta adequadamente às<br />
especificações de projeto.<br />
II. Número Identificador<br />
Cada aluno de cada <strong>turma</strong> receberá um número identificador de acordo com a sua<br />
posição na lista de chamada. O aluno poderá consultar o seu número identificador na<br />
Tabela 8 em anexo nesse documento.<br />
III. Técnica a ser utilizada<br />
A técnica a ser utilizada na realização do trabalho será a mesma proposta na<br />
apostila Aula 6. O diagrama de bode do processo será gerado pela função projeto.p<br />
através da seguinte linha de comando no matlab<br />
>> projeto(<strong>tipo</strong>,<strong>turma</strong>,pontos)<br />
Maiores detalhes sobre esta função podem ser obtidos pelo help da função, isto é:<br />
>> projeto<br />
Observação: acrescentar 60, 65 ou 70, conforme a versão do Matlab<br />
utilizado.<br />
A função projeto.p acompanha todas as etapas do projeto do compensador,<br />
apresentando inicialmente as especificações de projeto (Figura 1), o diagrama de<br />
Bode (Figura 2) e as respostas temporais de malha fechada considerando<br />
realimentação unitária e sinais de entrada <strong>tipo</strong> degrau e rampa (Figura 3). Durante o<br />
desenvolvimento do projeto, a função pergunta realiza uma serie de perguntas que<br />
devem ser respondidas de forma adequada para que o aluno possa continuar o projeto<br />
(Figura 4).<br />
Professor: José Felipe Haffner 1
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Figura 1: Visualização da especificação de projeto.<br />
Figura 2: Diagrama de Bode de malha aberta do processo.<br />
Professor: José Felipe Haffner 2
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Resposta temporal em malha fechada e realimentação<br />
unitária para uma entrada <strong>tipo</strong> degrau.<br />
Resposta temporal em malha fechada e realimentação<br />
unitária para uma entrada <strong>tipo</strong> rampa.<br />
Figura 3: Respostas temporais do processo operando em malha fechada.<br />
Figura 4: Durante o desenvolvimento do projeto são realizadas perguntas que<br />
devem ser respondidas pelo aluno.<br />
Professor: José Felipe Haffner 3
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Os parâmetros de entrada da função projeto.p são mostrados abaixo:<br />
>> projeto(<strong>tipo</strong>,<strong>turma</strong>,pontos,A,B,C,D,E,F,G,H)<br />
Especificação dos parâmetros de entrada:<br />
<strong>tipo</strong> -> é o número identificador do aluno que determina o <strong>tipo</strong> de sistema (varia de 1 a 32)<br />
<strong>turma</strong> -> é o nº da <strong>turma</strong>: 460,470,480 e 490<br />
pontos -> numero de pontos do diagrama de Bode<br />
A -> ganho proporcional inserido no processo<br />
B -> numerador da função de transferência do compensador de atraso de fase kc1*[1 zc1]<br />
Observação: O compensador de atraso de fase inserido deve ter ganho DC unitário<br />
C -> denominador da função de transferência do compensador de atraso de fase [1 pc1]<br />
D -> numerador da função de transferência do compensador de avanço de fase kc2*[1 zc2]<br />
Observação: O compensador de avanço de fase inserido deve ter ganho DC unitario<br />
E -> denominador da função de transferência do compensador de avanço de fase [1 pc2]<br />
F -> numerador da função de transferência do segundo compensador de avanço de fase kc3*[1 zc3]<br />
Observação: O Segundo compensador de avanço de fase inserido deve ter ganho DC unitario<br />
G -> denominador da função de transferência do segundo compensador de avanço de fase [1 pc3]<br />
H -> Parâmetro quando colocado em 1, gera a nota do trabalho.<br />
Se H = 0, o modo nota não está ativado.<br />
A tabela 1 mostra o uso da função projeto.p durante a realização do trabalho. Inicialmente<br />
são empregados três parâmetros de entrada e depois na medida do desenvolvimento do<br />
trabalho, são introduzidos novos parâmetros relativos aos compensadores projetados.<br />
Comando a ser usado no MatLab Seqüência do projeto do compensador<br />
projeto(<strong>tipo</strong>,<strong>turma</strong>,pontos) Com três parâmetros de entrada, usado para verificar o<br />
diagrama de Bode do processo.<br />
Além disso, retorna a especificação desejada: constante de<br />
erro de velocidade, sobre-sinal e tempo de pico.<br />
projeto(<strong>tipo</strong>,<strong>turma</strong>,pontos,A) Com quatro parâmetros de entrada, usado para verificar o<br />
ganho proporcional:<br />
projeto (<strong>tipo</strong>,<strong>turma</strong>,pontos,A,B,C) Com seis parâmetros de entrada, usado para verificar o<br />
projeto do compensador de atraso de fase:<br />
projeto (<strong>tipo</strong>,<strong>turma</strong>,pontos,A,B,C,D,E) Com oito parâmetros de entrada, usado para verificar o<br />
projeto do compensador de avanço de fase:<br />
projeto (<strong>tipo</strong>,<strong>turma</strong>,pontos,A,B,C,D,E,H) Com nove parâmetros de entrada, usado para gerar a nota<br />
de um projeto usando um compensador de avanço de fase:<br />
projeto (<strong>tipo</strong>,<strong>turma</strong>,pontos,A,B,C,D,E,F,G) Com dez parâmetros de entrada, usado para verificar o<br />
projeto de um segundo compensador de avanço de fase:<br />
projeto (<strong>tipo</strong>,<strong>turma</strong>,pontos,A,B,C,D,E,F,G,H) Com onze parâmetros de entrada, usado para gerar a nota<br />
de um projeto usando dois compensadores de avanço de<br />
fase:<br />
Tabela 1: Seqüência do uso da função projeto durante o desenvolvimento do<br />
trabalho.<br />
Observação: Quando o modo nota é ativado a função necessita de um tempo de<br />
processamento para atribuir a nota ao aluno, portanto o modo nota só deve ser<br />
ativado no termino do projeto.<br />
Professor: José Felipe Haffner 4
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Para a geração da nota final do projeto do compensador são considerados o<br />
atendimento das especificações de projeto e as respostas das perguntas formuladas<br />
durante o projeto.<br />
IV. Data limite<br />
A data limite para entrega dos trabalhos é 15/05/2008. Não serão aceitos os<br />
trabalhos entregues depois da data prevista. Pois a partir do dia 16/05/2008 a nota<br />
será zerada. A entrega dos trabalhos será realizada da seguinte maneira<br />
• Enviar um e-mail ao professor da disciplina com os dados apresentados na<br />
tabela 1 em anexo:<br />
Email do professor: jfelipe@ee.<strong>pucrs</strong>.br<br />
• NO e-mail enviado deve ter o seguinte “subject”:<br />
<strong>PSC</strong>_<strong>nomedoaluno</strong>_<strong>tipo</strong>_<strong>turma</strong><br />
Observação: no campo “<strong>nomedoaluno</strong>” deve ser colocado o nome do aluno, no<br />
campo “<strong>tipo</strong>” deve ser colocado o numero atribuído ao trabalho e no campo<br />
“<strong>turma</strong>”deve ser colocado a <strong>turma</strong> do aluno. Exemplo: O aluno Paulo da <strong>turma</strong><br />
480 recebeu o trabalho de numero 20, portanto ao enviar o e-mail colocará o<br />
seguinte “subject”:<br />
<strong>PSC</strong>_Paulo_20_480<br />
• Não serão corrigidos arquivos mandados por e-mail que não contenham o<br />
“subject” normalizado.<br />
• A correção será feita utilizando as funções copy e paste para a transferência<br />
da função enviada por e-mail ao Matlab (para a verificação da nota final).<br />
Recomenda-se que o aluno teste este procedimento antes de enviar o e-mail.<br />
• Caso as informações fornecidas pelo aluno não funcionem adequadamente a<br />
partir do procedimento definido no item anterior, serão descontados 2.0<br />
pontos da nota final do trabalho do aluno.<br />
Professor: José Felipe Haffner 5
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
V. Roteiro<br />
a) Tarefa<br />
O aluno deve acessar o site www.ee.<strong>pucrs</strong>.br/~gacs e obter o material necessário<br />
para realização da tarefa (Ver Figura 5). O número identificador para a realização da<br />
tarefa deve ser obtido no mesmo endereço.<br />
Figura 5: Localização da Tarefa no site do GACS.<br />
b) Uso da função projeto.p<br />
A técnica a ser utilizada na realização do trabalho será a mesma proposta na<br />
apostila Aula 6. O diagrama de bode do processo será gerado pela função projeto.p<br />
através da seguinte linha de comando no matlab<br />
>> projeto(<strong>tipo</strong>,<strong>turma</strong>,pontos)<br />
Observação: <strong>tipo</strong> é o número identificador obtido anteriormente, <strong>turma</strong> = 490, e<br />
pontos=500 (número de pontos considerado para construir o diagrama de Bode).<br />
Maiores detalhes sobre esta função podem ser obtidos pelo help da função, isto é:<br />
>> projeto<br />
c) Avaliação do processo<br />
O aluno deve observar atentamente o diagrama de Bode gerado pela função<br />
projeto e preencher a tabela 2.<br />
Kp<br />
Kv<br />
ω0dB<br />
ω180º<br />
Margem de ganho<br />
Margem de fase<br />
Especificação Valor<br />
Tabela 2: Especificação do processo definido pela função projeto.<br />
Professor: José Felipe Haffner 6
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
d) Realizar a comparação entre as especificações da tabela 1 com as<br />
respostas temporais apresentadas pela função projeto<br />
Verificar a constante Kv encontrada no diagrama de Bode com a obtida diretamente<br />
da resposta temporal a rampa (ver figura 6).<br />
Verificar a margem de fase e a encontrada no diagrama de Bode com o sobressinal e<br />
o tempo de pico obtida diretamente da resposta temporal ao degrau (ver figura 2).<br />
Resposta temporal em malha fechada e realimentação Resposta temporal em malha fechada e realimentação<br />
unitária para uma entrada <strong>tipo</strong> degrau.<br />
unitária para uma entrada <strong>tipo</strong> rampa.<br />
Figura 6: Respostas temporais do processo operando em malha fechada.<br />
e) Avaliar as especificações de projeto<br />
Com as especificações obtidas da função projeto, preencher a tabela 3.<br />
Figura 7: Visualização da especificação de projeto.<br />
Professor: José Felipe Haffner 7<br />
Kv
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Especificações Valor<br />
Kv desejado<br />
Tempo de pico desejado<br />
Sobre-sinal desejado<br />
Margem de fase desejada<br />
ξ (obtido a partir do sobre-sinal)<br />
ωn (obtido usando o tp e ξ)<br />
ω0dB desejada<br />
Tabela 3: Especificações desejadas do projeto.<br />
f) Uso do ganho proporcional<br />
A primeira etapa do projeto será a escolha de um ganho proporcional para ajustar<br />
a ω0dB do processo para o valor desejado obtido na Tabela 3.<br />
>> projeto(<strong>tipo</strong>,<strong>turma</strong>,pontos,ganho)<br />
onde: ganho é um ganho proporcional inserido no sistema de controle realimentado<br />
com objetivo de ajustar a ω0dB no valor desejado.<br />
Kp<br />
Kv<br />
ω0dB<br />
Especificação Valor<br />
ω180º<br />
Margem de ganho<br />
Margem de fase<br />
Tabela 4: Especificação do Processo com um ganho proporcional definido pelo aluno.<br />
g) Realizar a avaliação do ganho proporcional escolhido<br />
Verificar as especificações da tabela 4 com as respostas temporais ao degrau e a<br />
rampa. Verificar quando se deveria aumentar o ganho proporcional para ajustar o Kv<br />
do processo ao valor desejado da tabela 3.<br />
h) Realizar o projeto de um compensador de atraso de fase<br />
A segunda etapa do projeto consiste no projeto do compensador de atraso de fase<br />
com objetivo de ajustar a Kv do processo para o valor desejado obtido na Tabela 3<br />
sem altera significativamente as demais especificações.<br />
>> projeto(<strong>tipo</strong>,<strong>turma</strong>,pontos,ganho,num1,den1)<br />
Professor: José Felipe Haffner 8
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
onde: num1 é o numerador da função de transferência do compensador de atraso de<br />
fase e den1 é o denominador da função de transferência do compensador de atraso de<br />
fase. O ganho DC do compensador deve ser unitário e qualquer ganho que deva ser<br />
inserido no sistema para ajustar o Kv deve ser colocado no parâmetro de entrada da<br />
função projeto ganho.<br />
Ao usar a função projeto serão realizadas perguntas sobre o desenvolvimento do<br />
projeto. Quando as perguntas forem respondidas de forma correta a função continua o<br />
processamento até terminar o projeto. Quando a função termina o projeto , é<br />
retornado o diagrama de Bode do sistema compensado bem com as respostas<br />
temporais do sistema atuando com o compensador. Preencher a tabela 5 e verificar se<br />
o projeto do compensador foi adequado. Nesse momento do projeto deve-se verificar<br />
se a especificação Kv foi plenamente atendida.<br />
Kp<br />
Kv<br />
ω0dB<br />
Especificação Valor<br />
ω180º<br />
Margem de ganho<br />
Margem de fase<br />
Tabela 5: Especificação do Processo com um ganho proporcional e um compensador<br />
de atraso de fase definido pelo aluno.<br />
i) Realizar o projeto de um compensador de avanço de fase<br />
A próxima etapa de projeto primeira etapa do projeto consiste no<br />
desenvolvimento do projeto de um compensador de avanço de fase para ajustar o<br />
sobre-sinal do processo para o valor desejado obtido na Tabela 3.<br />
Professor: José Felipe Haffner 9
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
VI) ANEXO<br />
Observação: Se o aluno utilizar somente um compensador de avanço recomenda-se usar o<br />
modelo apresentado pela tabela 6. Caso o aluno utilize dois compensadores de avanço<br />
recomenda-se usar o modelo apresentado pela tabela 7.<br />
Tabela 6: Dados que devem ser fornecidos ao Professor.<br />
%Nome do aluno:<br />
% Trabalho Individual – Projetos de Sistemas de Controle<br />
%Dados<br />
<strong>tipo</strong>= ; % número identificador do aluno na lista de chamada<br />
<strong>turma</strong>= ; %Turma do aluno<br />
Kp= ; %Ganho proporcional utilizado<br />
%Função de transferência do compensador de atraso de fase projetado<br />
n_atraso=[ ]; % vetor com os coeficientes do polinômio do numerador<br />
d_atraso=[ ]; % vetor com os coeficientes do polinômio do denominador<br />
%Função de transferência do compensador de avanço de fase projetado<br />
n_avanco1=[ ]; % vetor com os coeficientes do polinômio do numerador<br />
d_avanco1=[ ]; % vetor com os coeficientes do polinômio do denominador<br />
%<br />
%Respostas ao questionário do Projeto do compensador de atraso de fase<br />
W0dbs= ; %Entre com a nova freqüência de 0 dB esperada<br />
Kvs= ; % Qual a nova constante de erro esperada(Kv)<br />
%Respostas ao questionário do Projeto do compensador de avanço de fase<br />
fit1= ; %Entre com a introdução de fase esperada<br />
magt1= ; %Entre com a introdução magnitude do compensador [dB]<br />
W0dbt1= ; %Entre com a nova freqüência de 0 dB esperada<br />
pmt1= ; % Qual a nova Margem de fase esperada<br />
%Calculo da nota se foi utilizado um compensador de avanço<br />
respostas=[W0dbs Kvs fit1 magt1 W0dbt1 pmt1];<br />
projeto(<strong>tipo</strong>,<strong>turma</strong>,1000,Kp,n_atraso,d_atraso,n_avanco1,d_avanco1,1)<br />
Professor: José Felipe Haffner 10
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Tabela 7: Dados que devem ser fornecidos ao Professor.<br />
%Nome do aluno:<br />
% Trabalho Individual – Projetos de Sistemas de Controle<br />
%Dados<br />
<strong>tipo</strong>= ; % número identificador do aluno na lista de chamada<br />
<strong>turma</strong>= ; %Turma do aluno<br />
Kp= ; %Ganho proporcional utilizado<br />
%Função de transferência do compensador de atraso de fase projetado<br />
n_atraso=[ ]; % vetor com os coeficientes do polinômio do numerador<br />
d_atraso=[ ]; % vetor com os coeficientes do polinômio do denominador<br />
%Função de transferência do compensador de avanço de fase projetado<br />
n_avanco1=[ ]; % vetor com os coeficientes do polinômio do numerador<br />
d_avanco1=[ ]; % vetor com os coeficientes do polinômio do denominador<br />
%<br />
%Respostas ao questionário do Projeto do compensador de atraso de fase<br />
W0dbs= ; %Entre com a nova freqüência de 0 dB esperada<br />
Kvs= ; % Qual a nova constante de erro esperada(Kv)<br />
%Respostas ao questionário do Projeto do compensador de avanço de fase<br />
fit1= ; %Entre com a introdução de fase esperada<br />
magt1= ; %Entre com a introdução magnitude do compensador [dB]<br />
W0dbt1= ; %Entre com a nova freqüência de 0 dB esperada<br />
pmt1= ; % Qual a nova Margem de fase esperada<br />
%Função de transferência do segundo compensador de avanço de fase projetado<br />
n_avanco2=[ ]; % vetor com os coeficientes do polinômio do numerador<br />
d_avanco2=[ ]; % vetor com os coeficientes do polinômio do denominador<br />
%Respostas ao questionário do Projeto do Segundo compensador de avanço de fase<br />
fit2= ; %Entre com a introdução de fase esperada<br />
magt2= ; %Entre com a introdução magnitude do compensador [dB]<br />
W0dbt2= ; %Entre com a nova freqüência de 0 dB esperada:<br />
pmt2= ; % Qual a nova Margem de fase esperada:<br />
%<br />
%Calculo da nota se foi utilizado dois compensadores de avanço<br />
respostas=[W0dbs Kvs fit1 magt1 W0dbt1 pmt1 fit2 magt2 W0dbt2 pmt2];<br />
projeto(<strong>tipo</strong>,<strong>turma</strong>,1000,Kp,n_atraso,d_atraso,n_avanco1,d_avanco1,n_avanco2,d_avanco2,1)<br />
Professor: José Felipe Haffner 11
PROJETOS DE SISTEMAS DE CONTROLE<br />
PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO GRANDE DO SUL<br />
DEPARTAMENTO DE ENGENHARIA ELÉTRICA<br />
Tabela 8: Número Identificador para ser usado no trabalho<br />
Número Nome do Aluno<br />
1 ADRIANO TEIXEIRA DE CASTRO<br />
2 ALEX REZENDE DE QUADROS<br />
3 ALEX STEDILE CORREIA<br />
4 ALEXSANDER ECKERT<br />
5 BIANCO SANTANA<br />
6 CARLOS SPIEKERMANN<br />
7 CIBELE BORDINI DE CASTRO<br />
8 CLAYTON ROMERO DIAS<br />
9 DANIEL LONGONI<br />
10 DORA GLADYS RAMOS CARHUAPOMA<br />
11 EDISON LUIS RIBEIRO<br />
12 EZEQUIEL DA ROSA FERNANDES<br />
13 FABIO BANDA ROLAND<br />
14 FABRICIO COSTA PETRY<br />
15 FREDERICO PEREIRA LAYDNER<br />
16 GABRIEL ALVES DE SOUZA<br />
17 GABRIEL VIEIRA JURUA DA SILVA<br />
18 GEISON SCHEID WINCK<br />
19 GIOVANI TOMEDI<br />
20 GUSTAVO RODRIGUES STEIGLEDER<br />
21 IGUATEMI MACEDO GOULART<br />
22 IURI KIRIYAMA FORTE<br />
23 JONATAS MATTIAZZI<br />
24 LEANDRO R DE MORAES ARISI<br />
25 LORENZO BASSO BENEVENUTTI<br />
26 LUCIANO KRONHARDT DE LIMA<br />
27 LUIZA PIANCA<br />
28 MARCELO LOPES DIAS ALVES<br />
29 MARCIO FERNANDES FIRMINO<br />
30 MARCOS HECK<br />
31 MARCOS LUIZ MANGONI<br />
32 MARCUS VINICIUS PERIM NAUFEL<br />
33 MAYCON ROSA LARGURA<br />
34 RAFAEL GUEDES ABREU<br />
35 RAFAEL LIMA TOLEDO<br />
36 RENAN BIEGELMEYER DA SILVA<br />
37 ROBERT RAFAEL CORREA DA SILVA<br />
38 RODRIGO JOSE PIERRY<br />
39 RODRIGO SOUZA ROCHA<br />
40 SANDRO MARCELO SCHMAEDECKE<br />
41 WOLFGANG RICARDO ROWELL JUNIOR<br />
Professor: José Felipe Haffner 12