MatLab

Engenharia Elétrica
6º Período
DISCIPLINA: Análise de Sistemas Lineares

PROFESSOR : Roselito de Albuquerque Teixeira

ATIVIDADE: Modelagem e Simulação de Sistemas
ALUNOS:

DATA: 12/06/2013

1 Sistema Suspensão de Automóvel
Considere o sistema suspensão de um automóvel descrito no Livro do Ogata 4ed. pag. 108, exercício A.3.16 e faça o que se pede:

1) Obtenha o modelo matemático do mesmo (função de transferência);

2) Considerando m = 355kg, b = 1860N/m/s e k = 14400N/m, obtenha a resposta deste sistema a um quebra-mola. Diferentes velocidades deverão ser testadas, quais sejam:
Substituindo os valores na função de transferência obtemos:

Engenharia Elétrica
6º Período
Para obter a resposta deste sistema a um quebra-mola utilizamos o software MATLAB cuja programação segue abaixo:

clear close all s=tf('s'); V = 5;
%Entrada de velocidade em KM/H v= (v*1000)/3600;
%KM/h para m/s t = ((1.4)/v);
[NUM,DEN] = PADE(t/2,25);
Ga =tf(NUM,DEN);
Gp = 0.2*(6.283*v)/1.4;
Gat=(0.2*(6.283*v)/1.4)*Ga;
%Entrada do sistema gm = s^2+(((6.283*v)/1.4)^2);
%Gs = ((1860*s+14400)/(355*s^2+1860*s+14400));
Gsu= ((684.3*s+14400)/(355*s^2+684.3*s+14400)) step(s*((Gp+Gat)/gm)*Gsu) %grid
%funçao para gerar o quebra-mola
% step(s*(Gp/gm))%gerar o quebra-mola v=3.6 hold on
%step(s*(Gat/gm))
hold on
%step(s*((Gp+Gat)/gm))
hold on
%step(s*((Gp+Gat)/gm)*Gs)
%grid hold on
%step(s*((Gp+Gat)/gm)*Gsu,'--r')
axis([0 10 -.3 .3])
Obs: Para cada valor de velocidade descrita teremos que alterar o valor da variável v. Segue abaixo a resposta demonstrada em gráfico de cada velocidade.
(a) 5 km/h
(b) 20 km/h
(c) 40 km/h.

Engenharia Elétrica
6º Período
(a) 5 km/h o sistema responde da seguinte forma:

(b) 20 km/h o sistema responde da seguinte forma:

Engenharia Elétrica
6º Período
(c) 40 km/h os sistema responde da seguinte forma:

Ao fazer as análises dos gráficos com as diferentes velocidades observamos que quanto mais a velocidade