03b Modeling

Download as pdf or txt
Download as pdf or txt
You are on page 1of 28

Mathematical Modeling

Review
Chapter 2

MECH 436 Mathematical modeling Slide 1


Cascaded, parallel, & feedback
transfer functions (Matlab)
a) Cascaded
b) Parallel-connected
c) Feedback-connected

Find the transfer


function:
C (s)
R( s)

MECH 436 Mathematical modeling Slide 2


Automatic controllers

• An automatic controller:
– Compares the value of the plant output with the reference
input,
– Determines the difference,
– Produces a control signal that reduces this difference to 0.
• The manner in which the automatic controller
produces the control signal is called the control
action.

MECH 436 Mathematical modeling Slide 3


Automatic controllers

Signal at very low e.g., electric motor,


power level hydraulic motor,
pneumatic motor

Converts the output variable into a


suitable variable (e.g., displacement,
pressure, voltage)

MECH 436 Mathematical modeling Slide 4


Classifying controllers (power)

• Industrial controllers can be classified based on the


kind of power employed:
– Electric controllers.
– Pneumatic controllers (pressurized air).
– Hydraulic controllers (pressurized oil).

MECH 436 Mathematical modeling Slide 5


Classifying controllers (action)
• Industrial controller can also be classified based on their
control action:
– 2 position (on-off) controllers.
– Proportional controllers (P).
– Integral controllers (I).
– Proportional + Integral controllers (PI).
– Proportional + Derivative controllers (PD).
– Proportional + Integral + Derivative controllers (PID).

MECH 436 Mathematical modeling Slide 6


2 position (on-off) control action

• The actuating element has 2 fixed positions (on or off).


• With output u(t) and error e(t):

u (t ) = U1 , for e(t)  0
= U 2 , for e(t)  0

MECH 436 Mathematical modeling Slide 7


Classifying controllers (action)
• Two-position controllers are usually electrical devices,
where an electric solenoid operated valve is widely
used.
• Pneumatic two-position controllers also exist and are
called pneumatic two-position controllers.

MECH 436 Mathematical modeling Slide 8


Proportional control action
• The value of the controller output u(t) is proportional to
the actuating error signal e(t):

u (t ) = K p e(t )

U ( s)
Transfer function: = Kp Proportional gain
E ( s)

E(s) U(s)
+ Kp
-

MECH 436 Mathematical modeling Slide 9


Integral control action
• The value of the controller output change (du(t)/dt) is
proportional to the actuating error signal e(t):
t
= K i e(t )  u (t ) = K i 0 e(t )dt
du(t )
dt

U ( s) K i
Transfer function: =
E ( s) s

E(s) Ki U(s)
+
- s

MECH 436 Mathematical modeling Slide 10


Proportional + Integral control

• The control action is defined by:


Kp t
u (t ) = K p e(t ) +
Ti  e(t )dt
0

Integral time
U ( s)  1 
Transfer function: = K p 1 + 
E ( s)  Ti s 

E(s) K p (1 + Ti s ) U(s)
+
- Ti s

MECH 436 Mathematical modeling Slide 11


Proportional + Derivative control

• The control action is defined by:

de(t )
u (t ) = K p e(t ) + K pTd
dt
Derivative time
= K p (1 + Td s )
U ( s)
Transfer function:
E ( s)

K p (1 + Td s )
E(s) U(s)
+
-

MECH 436 Mathematical modeling Slide 12


Proportional + Integral + Derivative
• This combination has the advantages of each of the
three individual control actions.

Kp t de(t )
Controller equation: u (t ) = K p e(t ) +
Ti 0
e(t )dt + K pTd
dt
U ( s)  1 
Transfer function: = K p 1 + + Td s 
E ( s)  Ti s 

+
E(s)
(
K p 1 + Ti s + TiTd s 2 ) U(s)

- Ti s

MECH 436 Mathematical modeling Slide 13


PID controller using Op-amps
Z1
Z2
C1
R2 C2
R4
R1
R3
-
-
Ei(s) +
E(s) + Eo(s)

E ( s) Z2 R1 R C s +1
=− where Z1 = ; Z2 = 2 2
Ei ( s) Z1 R1C1s + 1 C2 s

MECH 436 Mathematical modeling Slide 14


PID controller using Op-amps

E ( s)  R2C2 s + 1  R1C1s + 1 
= −  
Ei ( s)  C2 s  R1 

Eo ( s) R4
We also have: =−
E ( s) R3

Eo ( s)  Ti 
Multiplying: = K p 1 + + Td s 
Ei ( s)  s 

Eo ( s) R4 ( R1C1 + R2C2 )  1 R1C1R2C2 


= 1 + + s
Ei ( s) R3 R1C2  (R1C1 + R2C2 )s (R1C1 + R2C2 ) 

MECH 436 Mathematical modeling Slide 15


Disturbance effect

• Since this is a linear system, superposition holds:


– Treat the reference input and the disturbance independently
– Add the outputs to each of these inputs to get the complete
output.

MECH 436 Mathematical modeling Slide 16


1- Disturbance alone

• Assume the reference input is zero and find CD


CD ( s) G2 ( s)
=
D( s) 1 + G1 ( s)G2 ( s) H ( s)

MECH 436 Mathematical modeling Slide 17


2- Reference input alone

• Assume the disturbance is zero and find CR

CR (s) G1 (s)G2 (s)


=
R(s) 1+ G1 (s)G2 (s)H (s)

MECH 436 Mathematical modeling Slide 18


Total response
(reference + disturbance)
• Add both response together:

C ( s) = CR ( s) + CD ( s) =
G2 ( s)
G1 (s) R(s) + D(s)
1 + G1 ( s)G2 ( s) H ( s)

MECH 436 Mathematical modeling Slide 19


Analysis
• What happens when |G1(s)H(s)| >>1 and
|G1(s)G2(s)H(s)| >>1?

• The response becomes independent of the disturbance:

CD ( s) G2 ( s)
= 0
D( s) 1 + G1 ( s)G2 ( s) H ( s)

• The closed loop transfer function is independent of the gains:

CR ( s) G1 ( s)G2 ( s) 1
=
D( s) 1 + G1 ( s)G2 ( s) H ( s) H (s)

MECH 436 Mathematical modeling Slide 20


Drawing a block diagram
• To draw a block diagram:
1) Write the equations that describe the dynamic
behavior of each component.
2) Take the Laplace transform of these equations,
assuming zero initial conditions.
3) Represent each Laplace-transformed equation
individually in block form.
4) Assemble the blocks.

MECH 436 Mathematical modeling Slide 21


Example – RC circuit
Start by writing the equations:

i=
ei − eo
eo =
 idt
R C
Next take the Laplace transforms:
Represent each equation in blocks:
Ei ( s) − Eo ( s)
I (s) =
R

I ( s)
Eo ( s) =
Cs

MECH 436 Mathematical modeling Slide 22


Example – RC circuit

Finally assemble the block into a complete block diagram:

MECH 436 Mathematical modeling Slide 23


Block diagram reduction

• Blocks can be connected in series only if the output of


one block is not effected by the next following block.
• If there are any loading effects between the
components, they must be combined into a single
block instead of cascading.
• Any number of cascaded blocks representing non-
loading components can be replaced by a single
block. In this case the resulting transfer function is the
product of the individual transfer functions.

MECH 436 Mathematical modeling Slide 24


Block diagram reduction

• We can reduce a complicated block diagram by step-


by-step rearrangement.
• The net effect is to reduce the subsequent mathematical
analysis.
• Unfortunately, the transfer functions in the new blocks
become more complex.

MECH 436 Mathematical modeling Slide 25


Example
H2

R - C
+- + G1 + G2 G3
+

H1

H2
G1
R - C
+- + + G1 G2 G3
+

H1

MECH 436 Mathematical modeling Slide 26


H2
G1
R - G1G2 C
+- + G3
1 − G1G2 H1

R G1G2G3 C
+-
1 − G1G2 H1 + G2G3 H 2

R G1G2G3 C
1 − G1G2 H1 + G2G3 H 2 + G1G2G3

MECH 436 Mathematical modeling Slide 27


Example

• The numerator for the closed loop transfer function is


the product of the transfer functions of the
feedforward path.
• The denominator for the closed loop transfer function
is:

1 −  product of the transfer functions around each loop =

1 − G1G2 H1 + G2G3 H 2 + G1G2G3

MECH 436 Mathematical modeling Slide 28

You might also like