Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
48 views

Aircraft Controls Assignment - 2: by Pothula Shiva Harsh SR - No. 10197 Aerospace Engineering (M.E)

This document contains Matlab code and summaries of simulations analyzing the response of an aircraft's pitch rate controller to disturbances. It tests proportional (P), proportional with lag compensation, and proportional with fast lag compensation controllers. The P controller has high overshoot and slow response. Lag compensation reduces overshoot but not significantly. Fast lag compensation achieves a faster response time with much lower overshoot. Root locus and step response plots are generated to analyze stability and performance.

Uploaded by

Shiva Harsh
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
48 views

Aircraft Controls Assignment - 2: by Pothula Shiva Harsh SR - No. 10197 Aerospace Engineering (M.E)

This document contains Matlab code and summaries of simulations analyzing the response of an aircraft's pitch rate controller to disturbances. It tests proportional (P), proportional with lag compensation, and proportional with fast lag compensation controllers. The P controller has high overshoot and slow response. Lag compensation reduces overshoot but not significantly. Fast lag compensation achieves a faster response time with much lower overshoot. Root locus and step response plots are generated to analyze stability and performance.

Uploaded by

Shiva Harsh
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 10

Aircraft Controls Assignment -2

By
Pothula Shiva Harsh
Sr.No. 10197
Aerospace Engineering (M.E)

Simulink model:

Matlab code for P controller


clear all;
qc=1;
k=0.713;
sim('design_a.mdl');
figure;
plot(q.time,q.signals.values(:,1),q.time,q.signals.values(:,2),q.time,q.signa
ls.values(:,3),'linewidth',3);grid on;zoom on;
xlabel('\bf Time (sec)');
ylabel('\bf Response q_c, q_f')
title('\bf ');
xlabel('Time (sec)','FontName','default','FontSize',14,'FontWeight','Bold');
ylabel('Response q_c, q_f (deg/s)');
title('Pitch Rate response- P controller');
legend('Response','Gust','Command');
figure;
LTF=tf(17.4*[1 1.56],[1 1.568 11.5]);
rlocus(LTF);
sys = tf([17.4 87],[1 1.568 11.5]);
S = stepinfo(sys,'RiseTimeLimits',[0.05,0.95]);

Output:
RiseTime:
SettlingTime:
SettlingMin:
SettlingMax:
Overshoot:
Undershoot:
Peak:
PeakTime:

0.3084
4.8256
2.7171
6.0284
59.3724
0
6.0284
0.7636

Response:

Respectiverootlocus

GustResponseMatlabCode:
sys = tf([17.4 27.144], [ 1 1.568 11.5]);
sys1 = tf(.742);
f = tf (1);
sys2 = series(sys1,lag);
sys3 = series(sys2,f);
h = feedback(sys,sys3)
step(h);

Gust Response

Pcontrollerwithlagcompensator:

Matlabcode:
clear all;
qc=1;
k=0.713;
z1=0.6;
p1=0.02;
sim('design_b.mdl');
figure;
plot(q.time,q.signals.values(:,1),q.time,q.signals.values(:,2),q.time,q.signa
ls.values(:,3),'linewidth',3);grid on;zoom on;
xlabel('\bf Time (sec)');
ylabel('\bf Response q_c, q_f')
title('\bf ');
xlabel('Time (sec)');
ylabel('Response q_c, q_f (deg/s)');
title('Pitch Rate resoponse- Lag controller');
legend('Response','Gust','Command');
figure;
G=tf(17.4*[1 1.56],[1 1.568 11.5]);
C=zpk([-z1],-p1,k);
LTF=C*G;
rlocus(LTF);
sys = tf(LTF);
S = stepinfo(sys,'RiseTimeLimits',[0.05,0.95])

Response:
RiseTime: 0.3084
SettlingTime: 4.825
SettlingMin: 2.7171
SettlingMax: 6.0284
Overshoot: 59.3724
Undershoot: 0

Peak: 6.0284
PeakTime: 0.7636

Rootlocus

Response:

GustResponse:

Lag compensator with fast lag compensation


clear all;
qc=1;
k=0.713;
z1=4;
p1=0.14;
sim('design_b.mdl');
figure;
plot(q.time,q.signals.values(:,1),q.time,q.signals.values(:,2),q.time,q.signa
ls.values(:,3),'linewidth',3);grid on;zoom on;
xlabel('\bf Time (sec)');
ylabel('\bf Response q_c, q_f')
title('\bf ');
xlabel('Time (sec)');
ylabel('Response q_c, q_f (deg/s)');
title('Pitch Rate response- Lag controller');
legend('Response','Gust','Command');
figure;
G=tf([17.4 27.144],[1 1.568 11.5]);
C=zpk([-z1],-p1,k);
LTF=((C*G)/(1+(C*G)));
rlocus(LTF,[0:0.001:2]);
pause
sys = tf(LTF);
S = stepinfo(sys,'RiseTimeLimits',[0.05,0.95])
close all

Output:
RiseTime: 0.1323
SettlingTime: 1.8809
SettlingMin: 0.9016
SettlingMax: 1.0554
Overshoot: 7.7375
Undershoot: 0
Peak: 1.0554
PeakTime: 0.2428

Response

Rootlocus:

Matlabcodeforgustresponse
sys = tf([17.4 27.144], [ 1 1.568 11.5]);
lag= tf ([1 4],[1 .14]);
sys1 = tf(.713);
f = tf (1);
sys2 = series(sys1,lag);
sys3 = series(sys2,f);
h = feedback(sys,sys3)
step(h);

GustResponse:

You might also like