29.07.2014 Views

Chemical Process Control a First Course with Matlab

Chemical Process Control a First Course with Matlab

Chemical Process Control a First Course with Matlab

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

3 - 2<br />

are the two transfer functions for the two inputs X 1 (s) and X 2 (s).<br />

Take note (again!) that the characteristic polynomials in the denominators of both transfer<br />

functions are identical. The roots of the characteristic polynomial (the poles) are independent of<br />

the inputs. It is obvious since they come from the same differential equation (same process or<br />

system). The poles tell us what the time-domain solution, y(t), generally would "look" like. A<br />

final reminder: no matter how high the order of n may be in Eq. (3-4), we can always use partial<br />

fractions to break up the transfer functions into first and second order terms.<br />

3.1 <strong>First</strong> order differential equation models<br />

This section is a review of the properties of a first order differential equation model. Our Chapter 2<br />

examples of mixed vessels, stirred-tank heater, and homework problems of isothermal stirred-tank<br />

chemical reactors all fall into this category. Furthermore, the differential equation may represent<br />

either a process or a control system. What we cover here applies to any problem or situation as<br />

long as it can be described by a linear first order differential equation.<br />

We usually try to identify features which are characteristic of a model. Using the examples in<br />

Section 2.8 as a guide, a first order model using deviation variables <strong>with</strong> one input and <strong>with</strong><br />

constant coefficients a 1 , a o and b can be written in general notations as 1<br />

The model, as in Eq. (2-2), is rearranged as<br />

a 1<br />

dy<br />

dt +a o y = bx(t) <strong>with</strong> a 1 ≠ 0 and y(0) = 0 (3-5)<br />

τ dy<br />

dt<br />

+y = Kx(t) (3-6)<br />

where τ is the time constant, and K is the steady state gain.<br />

In the event that we are modeling a process, we would use a subscript p (τ = τ p , K = K p ).<br />

Similarly, the parameters would be the system time constant and system steady state gain when we<br />

analyze a control system. To avoid confusion, we may use a different subscript for a system.<br />

y<br />

1.2<br />

MK 1<br />

0.8<br />

0.6<br />

0.4 increasing τ<br />

0.2<br />

0<br />

0 2 4 6 8<br />

t<br />

10 12<br />

2<br />

K = 2<br />

1.5<br />

y<br />

M<br />

K = 1<br />

1<br />

0.5<br />

K = 0.5<br />

0<br />

0 2 4 6 t 8 10<br />

Figure 3.1. Properties of a first order transfer function in time domain. Left panel y/MK: effect of<br />

changing the time constant; plotted <strong>with</strong> τ = 0.25, 0.5, 1, and 2 [arbitrary time unit]. Right panel y/M:<br />

effect of changing the steady state gain; all curves have τ = 1.5.<br />

1 Whether the notation is y or y‘ is immaterial. The key is to find the initial condition of the<br />

problem statement. If the initial condition is zero, the notation must refer to a deviation variable.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!