Control Systems Lab: Submitted by

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 4

Control Systems Lab

P2

Submitted By:
Syed Muhammed Baqir Naqvi (170984)
Muhammed Inam Anwar (170966)
Qazi Ahasan Nawaz (171008)
Class/Sec: BEEP 5A
Submitted To: Engr. Hussain Asif
Instructor’s Remarks:

Submission Time: In Late in Working Days


Time 0 1 2 3 4 5

Marks: Total Marks Obtained Marks

Additional Remarks: _______________________________________


________________________________________________________

Instructor’s Signature: _____________


Complex Engineering Problem
Transfer Function Generator and PID
Controller
1. Objectives:
 To develop Transfer Function regarding the behavior of high tower buildings.
 Also to design PID controller to improve transient response.

2. Procedures:
 We need to develop response for following system.

 B1 represents the energy dissipation as the building moves.


 K1 is the lateral stiffness of the building structure.
 m2 is the mass of the sliding element (shaft and voice-coil formers).
 B2 is the viscous fiction coefficient describing the eddy-current losses in the voice coils.
 K2 is the stiffness of the wire spring.

We can summarize our system for the transfer function.

( m 1 s2 + B 1 s+ K 1 ) V 1 ( s )− ( B 2 s+ K 2 ) V 2 ( s )=Fw ( s ) … … .. ( i )
−( B 2 s+ K 2 ) V 1 ( s ) + ( m 2 s2 + B 2 s+ K 2 ) V 2 ( s ) =( s )!… … . ( ii )

m1 s 2+ B 1 s + K 1 −( B 2 s+ K 2 ) V 1( s) Fw (s )
[ −( B 2 s + K 2 ) 2
][ ] [ ]
( m 2 s + B 2 s+ K 2 ) V 2( s)
=
(s)!

 By using Crammer’s Rule we can evaluate the value of V1(s)

Fw ( s ) ( m1 s 2+ B 2 s + K 2 ) ( s )!(B 2 s + K 2)
V 1( s) = +
( a s +bs + cs + ds+e ) −( B 2 s+ k 2) ( a s +bs + cs2 +ds +e )−(B 2 s +k 2)2
4 3 2 2 4 3

 Whereas, a=m1m2, b=(m1B2+m2B1), c=(B1B2+m1K2+m2K1), d=(B1K2+B2K1), e=(K1K2)


Fw ( s ) ( m1 s 2+ B 2 s + K 2 )
V 1( s) =
a s 4 +bs 3+¿¿ ¿ ¿

 In order to get V1(s)/Fw(s) we assume there is no acting force on the system.

V 1 ( s ) ( m1 s2 + B2 s+ K 2 )
=
Fw ( s ) a s 4 +bs 3+¿ ¿¿ ¿
Let us have m1=1, B2=1, K2=1, a=1, b=1, c=1, d=1, e=1.

V 1 ( s ) ( s2 + s+1 )
=
Fw ( s ) s 4 + s3 +¿¿ ¿ ¿
Code:
a=[0 0 1 1 1]
b=[1 1 0 3 4]
y=tf(a,b)
sisotool(y)
Result:

 In order to get V1(s)/Fact(s) we assume there is no acting force on the system.

V 1( s) ( B2 s+ K 2)
= 4
( s )! a s +bs 3+¿ ¿¿ ¿
 Let us have B2=1, K2=1, a=1, b=1, c=1, d=1, e=1.
V 1( s) (s+ 1)
= 4 3 +¿¿¿
( s )! s +s ¿
Code:
a=[0 0 0 1 1]
b=[1 1 0 3 4]
y=tf(a,b)
sisotool(y)
Result:
Learning Outcomes:
 We were able familiarize ourselves with MATLAB’s SISO Design Tool.
 We saw PID high tower buildings require precise measurements regarding control systems accurate
approach to transfer functions is required.

You might also like