Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Problem Set 6

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

Harmattan 2021/2022 PROCESS DYNAMICS CHE 509

PROBLEM SET 6

Assigned out: 22/11/2022 Due date: 1 week

Goals: The main goals of this problem set are to equip students with the following skills sets:

1) ability to empirically model dynamical systems from process data,

2) ability to use MATLAB/SIMULINK system identification toolbox for process identification.

Time: This problem set will require you to put in an average of 3 hours. Don’t be intimidated by this.
It doesn’t have to be all at once. You can take break in between. Please have fun while solving the
problem. Best wishes!

[1] A single-tank process has been operating for a long period of


time with the inlet flow rate qi equal to 30.1 ft3/min. After the
operator increases the flow rate suddenly at t=0 by 10%, the liquid
level in the tank changes as shown in Table 6.1.
Assuming that the process dynamics can be described by a first-order
model, calculate the steady-state gain and the time constant using
three methods:
(a) From the time required for the output to reach 63.2% of the total
change

(b) From the initial slope of the response curve

(c) By fitting the data to the transient response equation of a first-order


system using nonlinear regression.

(d) Compare the data and the three models by simulating their step
responses.
Table 6.1
[2] Fig. 6.1 presents the response of a system to a unit step in the
input.
(a) Use these data to derive a FOPTD model of this system
(b) Plot the response of the model and compare with the data
(c) The FOPTD model does not capture all the features of the
original response. What is the feature that the FOPTD model fails
to capture? What causes this type of dynamic behavior?

Fig.6.1: Step response plot for problem [2]


[3]Consider the step response data y(k) in Table 6.2 below. At t=0, a unit
step change in u occurs, but the first output change is not observed until the
next sampling instant. Estimate the model parameters in the second order
difference equation from the input-output data using linear regression and
compare the results obtained by using the system identification toolbox
subroutine “arx” in MATLAB.

𝐲(𝐤) = 𝒂𝟏 𝒚(𝒌 − 𝟏) + 𝒂𝟐 𝒚(𝒌 − 𝟐) + 𝒃𝟏 𝒖(𝒌 − 𝟏) + 𝒃𝟐 𝒖(𝒌 − 𝟏)

Table 6.2: Step Response Data

k y(k)
1 0.060
2 0.222
3 0.300
4 0.359
5 0.450
6 0.480
7 0.520
8 0.600
9 0.692
10 0.700
11 0.771
12 0.833
13 0.888
14 0.900
15 0.910
16 0.925
17 0.934
18 0.950
19 0.965
20 0.980
N.B.: y(k) = 0 for k < 0; u(k) = 0 for k < 0

You might also like