Modelling & Design SimuLink
Modelling and Design of Feedback Controllers with SIMULINK
Problem 1: Cooling system
Consider the well-mixed cooling system shown below:
Cooling system design and operating conditions:
Jacket heattransfercoefficient: U = 1.225×106 J/hr.m2.C Area ofjacket/tankinterface: A = 31m2
Steady-stateflowrate: F = 2.83m3/hr
Fluiddensity: r= 993Kg/m3
Fluidheatcapacity: Cp = 4190 J/KgC
Mass of fluidintank: M = 18,200 Kg Steady-stateinlettemperature: T1 = 82C
Steady-stateoutlettemperature: T2 = 38 C (setpoint)
You may assume the fluid volume in the tank is constant and the fluid heat exchanger properties are constant but T1, T2, Tw, and F may change.
- Derive the transfer functions that relate T2 to T1, Tw, and F. Calculate the values for time constant andgains
- Sketch the instrumentation diagram for the process (transmitter, controller, controlvalve)
- With Simulink determine the PID tuning parameters using the closed- loop (constant oscillation) tuning method. Set point is constant. T1 changes –5 oC and F changes –0.4m3/hr
Additional information: Process delay = 0.5 hr
Transmitter output range = 4 – 20 mA Transmitter calibration range = 26 – 42 oC Transmitter time constant = 0
Valve gain = 1
Problem 2: Heat exchanger
Consider the process within a heat exchanger as given below:
Acetone enters at 120°C, heating up the incoming decane stream from 80°C to 95°C (at steady state). The outlet temperature of decane stream is maintained by manipulation of acetone inlet valve. At steady state operation, the flow rate of decane is 390 kmol/hr.
At time, t = 0, the decane flow rate is suddenly decreased to 250 kmol/hr, affecting the outlet decane temperature as shown in the response table shown below. For the second step testing, the acetone valve is suddenly opened from 20 to 30 %, giving a temperature response also tabulated in the table below.
Outlet Decane Temperature Responses to Step Tests
|Time (min)||Decane Temp|
|Time (min)||Decane Temp|
Decane 390 to250kmol/hr Acetone Valve 20% to30%
Valve gain is between 1 and 2 with zero time constant.
Transmitter has unity gain with negligible dead time and time constant
- Tune PI and PID controllers using the open loopequations
- Provide some conclusions
- Find the constant oscillation
This is my simulation block.
First I set PID coefficients Kp= 20, Ki = 0, Kd = 0
The result is
The oscillation is descending, so should increase Kp
When Kp = 40
The ouput is over oscillating. I need to decrease Kp
When Kp = 33
The output is constantly oscillating so from this value we can find the PID coefficients.
Kp = 16.5 / 2 * 1.7 = 14
Ti = 0.50 * Tk = 1.4
Td = 0.125 * Tk = 0.35
Kp = 14 Ki = 10 Kd = 4.9
Using this coefficient we get the following result.
Problem 2. Heat Exchanger
Using the data output of step input, we can plot the graph using MATLAB
This is the change of the acetone temperature versus the step input of the flow rate of decane.
With linear approximation we can get linear model.
The following MATLAB graph shows the output of the temperature when step input the acetone valve change.
From this graph we can get PID controller coefficients.