CHAPTER 19: Single-Loop IMC Control

Similar documents
CHAPTER 10: STABILITY &TUNING

CHAPTER 13: FEEDBACK PERFORMANCE

CHAPTER 15: FEEDFORWARD CONTROL

Smith Predictor Based-Sliding Mode Controller for Integrating Process with Elevated Deadtime

2.9 Feedback and Feedforward Control

Solutions for Tutorial 10 Stability Analysis

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.010: Systems Modeling and Dynamics III. Final Examination Review Problems

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

Ch 12: Variations on Backpropagation

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

Chapter 8. Feedback Controllers. Figure 8.1 Schematic diagram for a stirred-tank blending system.

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon

Ph 20.3 Numerical Solution of Ordinary Differential Equations

The Wilson Model of Cortical Neurons Richard B. Wells

Architectural issues in the control of tall multiple-input multiple-output plants

Pseudo Dynamic Model Reference Nonlinear Control of a Continuous Bioreactor with Input Multiplicities

ma x = -bv x + F rod.

Torsion Experiment. Encoder #3 ( 3 ) Third encoder/disk for Model 205a only. Figure 1: ECP Torsion Experiment

Chapter 2: Introduction to Damping in Free and Forced Vibrations

Ufuk Demirci* and Feza Kerestecioglu**

COS 424: Interacting with Data. Written Exercises

Elastic Force: A Force Balance: Elastic & Gravitational Force: Force Example: Determining Spring Constant. Some Other Forces

Determining Controller Constants to Satisfy Performance Specifications

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. IX Uncertainty Models For Robustness Analysis - A. Garulli, A. Tesi and A. Vicino

Ocean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng

Autotuners based on the Smith predictor structures

Gary J. Balas Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, MN USA

Numerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term

Time-of-flight Identification of Ions in CESR and ERL

Efficient Filter Banks And Interpolators

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme

DRAFT. Memo. Contents. To whom it may concern SVN: Jan Mooiman +31 (0) nl

NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT

CHBE320 LECTURE XI CONTROLLER DESIGN AND PID CONTOLLER TUNING. Professor Dae Ryook Yang

A Novel Principle for Relay-Based Autotuning

Enhanced Single-Loop Control Strategies Chapter 16

Dynamics and PID control. Process dynamics

a a a a a a a m a b a b

Reducing Vibration and Providing Robustness with Multi-Input Shapers

Symbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm

Sensorless Control of Induction Motor Drive Using SVPWM - MRAS Speed Observer

Kinematics and dynamics, a computational approach

Non-Parametric Non-Line-of-Sight Identification 1

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

A Simple Regression Problem

CS Lecture 13. More Maximum Likelihood

Decentralized Adaptive Control of Nonlinear Systems Using Radial Basis Neural Networks

Feature Extraction Techniques

EMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS

16.30/31 September 24, 2010 Prof. J. P. How and Prof. E. Frazzoli Due: October 15, 2010 T.A. B. Luders /31 Lab #1

Analyzing Simulation Results

Pattern Recognition and Machine Learning. Artificial Neural networks

Principles of Optimal Control Spring 2008

Feedforward Networks

Figure 1: Equivalent electric (RC) circuit of a neurons membrane

Robustness Experiments for a Planar Hopping Control System

Dimensions and Units

Design of Sliding Mode Stabilizer for Wind Turbine Generator using Dynamic Compensation Observer Technique

Process Control & Design

The Weierstrass Approximation Theorem

Note-A-Rific: Mechanical

Data-Driven Imaging in Anisotropic Media

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics

All Excuses must be taken to 233 Loomis before 4:15, Monday, April 30.

lecture 36: Linear Multistep Mehods: Zero Stability

The Methods of Solution for Constrained Nonlinear Programming

A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)

Feedforward Networks

Pattern Recognition and Machine Learning. Artificial Neural networks

Lecture 8 Symmetries, conserved quantities, and the labeling of states Angular Momentum

Feedback Control of Linear SISO systems. Process Dynamics and Control

Kernel Methods and Support Vector Machines

Some Perspective. Forces and Newton s Laws

DESIGN, SIMULATION AND IMPLEMENTATION OF DECENTRALIZED PI CONTROLLERS FOR A MULTIVARIABLE FEED BLENDING SYSTEM AT THE FALCONBRIDGE SUDBURY SMELTER

Feedforward Networks. Gradient Descent Learning and Backpropagation. Christian Jacob. CPSC 533 Winter 2004

Donald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.

Hybrid System Identification: An SDP Approach

Curious Bounds for Floor Function Sums

This model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.

821. Study on analysis method for deepwater TTR coupled vibration of parameter vibration and vortex-induced vibration

LogLog-Beta and More: A New Algorithm for Cardinality Estimation Based on LogLog Counting

Handwriting Detection Model Based on Four-Dimensional Vector Space Model

Control Lab. Thermal Plant. Chriss Grimholt

Chapter 1: Basics of Vibrations for Simple Mechanical Systems

Sharp Time Data Tradeoffs for Linear Inverse Problems

Last Name Student Number. Last Name Student Number [5] Q1. Circuit analysis. [2] (a) For the following circuit, give the truth table.

Course support: Control Engineering (TBKRT05E)

Ştefan ŞTEFĂNESCU * is the minimum global value for the function h (x)

Finite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields

Using a De-Convolution Window for Operating Modal Analysis

Process Control, 3P4 Assignment 6

A Theoretical Analysis of a Warm Start Technique

System Design of Quadrotor

Easy Evaluation Method of Self-Compactability of Self-Compacting Concrete

Topic 5a Introduction to Curve Fitting & Linear Regression

A Simulation Study for Practical Control of a Quadrotor

Transcription:

When I coplete this chapter, I want to be able to do the following. Recognize that other feedback algoriths are possible Understand the IMC structure and how it provides the essential control features Tune an IMC controller Correctly select between PID and IMC

Outline of the lesson. Thought exercise for odel-based control IMC structure Desired control features IMC algorith and tuning Application guidelines

Let s quickly review the PID algorith E ( t ) SP ( t ) CV ( t ) 1 d CV MV ( t ) K c E ( t ) E ( t' ) dt ' Td T I 0 dt This is the de facto standard! I PID was developed in 1940 s PID is not the only feedback algorith PID gives good balance of perforance and robustness PID does not always give the best perforance Multiple PIDs are used for ultivariable systes

Let s look ahead to the IMC structure and algorith We will have another algorith to learn!!!! IMC was developed forally in 1980 s, but the ideas began in 1950 s IMC uses a process odel explicitly IMC involves a different structure and controller IMC could replace PID, but we chose to retain PID unless an advantage exists A single IMC can be used for ultivariable systes

Let s do a thought experient: 1. We want to control the concentration in the tank. 2. Initially, A 6 wgt% We want A 7 wgt% 3. Fro data, we know that A/ v 0.5 wgt%/% open Solvent Reactant A A wgt% 6% What do we do? Valve % open tie

Let s do a thought experient: 1. We want to control the concentration in the tank. 2. Initially, A 6 wgt% We want A 7 wgt% 3. Fro data, we know that A/ v 0.5 wgt%/% open Solvent Reactant A A wgt% 6% What do we expect to happen? Valve % open 2% tie

Let s do a thought experient: 1. We want to control the concentration in the tank. 2. Initially, A 6 wgt% We want A 7 wgt% 3. Fro data, we know that A/ v 0.5 wgt%/% open 7.3% A wgt% 6% Solvent Reactant A Should we be surprised? What do we do now? Devise 2 different responses Valve % open 2% tie Hint: We used a odel for the first calculation. What do we do if the odel is in error?

Let s do a thought experient: 1. We want to control the concentration in the tank. 2. Initially, A 6 wgt% We want A 7 wgt% 3. Fro data, we know that A/ v 0.5 wgt%/% open 7.3% A wgt% 6% 6.91% Solvent Reactant A We apply the feedback until we have converged. Valve % open 2% -0.6% tie

7.3% A wgt% 6% 5.91% Solvent Valve % open 2% -0.6% A tie Reactant D(s) d (s) SP(s) MV(s) C (s) v (s) P (s) - CV(s) Discuss this diagra and coplete the closed-loop block diagra for the IMC concept.

7.3% A wgt% 6% 5.91% Solvent Valve % open 2% -0.6% A tie Reactant D(s) d (s) SP(s) MV(s) C (s) v (s) P (s) - CV(s) This is a odel of the feedback process. (s) - Model error E (s)

D(s) d (s) SP(s) T p (s) MV(s) CP (s) v (s) P (s) - CV(s) What is CP?? E (s) (s) - Transfer functions CP (s) controller v (s) valve P (s) feedback process (s) odel d (s) disturbance process Variables CV(s) controlled variable CV (s) easured value of CV(s) D(s) disturbance E (s) odel error MV(s) anipulated variable SP(s) set point T p (s) set point corrected for odel error

What controller calculation gives good perforance? It is NOT a PID algorith Let s set soe key features and deterine what cp will achieve these features 1. Zero steady-state offset for step-like inputs 2. Perfect control (CVSP for all tie) 3. Moderate anipulated variable adjustents 4. Robustness to odel isatch 5. Anti-reset-windup

What controller calculation gives good perforance? SP(s) - T p (s) MV(s) CP (s) What is CP?? D(s) v (s) d (s) P (s) (s) - CV(s) E (s) 1. Zero steady-state offset: CV D( s) 1 CP (1 cp ) ( ) v p S d What condition for cp (s) ensures zero steady-state offset for a step disturbance? Hint: Steady-state is at t.

What controller calculation gives good perforance? SP(s) - T p (s) MV(s) CP (s) What is CP?? D(s) v (s) d (s) P (s) (s) - CV(s) E (s) 1. Zero steady-state offset: CV D( s) li CV ( t) t 1 CP li s 0 (1 cp ( ) v D (1 KcpK s s 1 K ( K p cp ) p ) K K Oh yeah, the final value theore! S d d 0 ) CONCLUSION K cp (K ) -1 Easily achieved because both are in coputer!

What controller calculation gives good perforance? SP(s) - T p (s) MV(s) CP (s) What is CP?? D(s) v (s) d (s) P (s) (s) - CV(s) E (s) 2. Perfect dynaic control: CV D( s) 1 CP (1 cp ) ( ) v p S d What is required for perfect dynaic control? Do you expect that we will achieve this goal?

What controller calculation gives good perforance? SP(s) - T p (s) MV(s) CP (s) What is CP?? D(s) v (s) d (s) P (s) (s) - CV(s) E (s) 2. Perfect dynaic control: CV D( s) 1 (1 cp ( s)) d CV D( s) 0 1 CP CP (1 ( ) v cp ( ) v p S p ) S d CONCLUSION cp (s) ( (s)) -1 Controller is inverse of odel! Can we achieve this?

2. Perfect dynaic control: 3. Moderate anipulated variable adjustents θs Ke ( τs 1) then cp ( τs 1) e K θs CONCLUSION cp (s) ( (s)) -1 Controller is inverse of odel! Can we achieve this? This is a pure derivative, which will lead to excessive anipulated variable oves This is a prediction into the future, which is not possible Conclusion: Perfect control is not possible! (See other exaples in the workshops.)

2. Perfect dynaic control: 3. Moderate anipulated variable adjustents CONCLUSION cp (s) ( (s)) -1 K( τ1s 1) ( τ s 1)( τ s 1) 2 3 then cp ( τ 2s 1)( τ 3s 1) K ( τ s 1) 1 Controller is inverse of odel! Can we achieve this? This has a second derivative, which will lead to excessive anipulated variable oves This could have an unstable controller, if τ 1 is negative. This is unacceptable. Conclusion: Perfect control is not possible! (See other exaples in the workshops.)

Let s begin our IMC design with the results so far. CONCLUSION K cp (K ) -1 Easily achieved because both are in coputer! CONCLUSION cp (s) ( (s)) -1 Controller is inverse of odel! We have loosened the restriction to a condition that can be achieved. Now, what does approxiate ean? How can we define the eaning of approxiate so that we have a useful design approach?

Separate the odel into two factors, one invertible and the other with all non-invertible ters. The invertible factor has an inverse that is causal and stable, which results in an acceptable controller. The gain is the odel gain, K. The non-invertible factor has an inverse that is noncausal or unstable. The factor contains odels eleents with dead ties and positive nuerator zeros. The gain is the 1.0. What do we use for the controller?

Separate the odel into two factors, one invertible and the other with all non-invertible ters. The IMC controller eliinates all non-invertible eleents in the feedback process odel by inverting - (s). cp [ ] -1 Looks easy, but I need soe practice.

F S solvent F A pure A AC Class exercise: We have two odels for the feedback dynaics for the 3-tank ixer. Deterine cp (s) for each. Epirical odel 5.5s 0.039e (1 10.5s) Fundaental odel 0.039 (1 5s) 3

Epirical odel 5.5s 0.039e (1 10.5s) Fundaental odel 0.039 (1 5s) 3 (1 10.5s) cp 0.039 ( 1 5s ) cp ( s ) 0. 039 3 Discuss these results. Do they ake sense? Are there any shortcoings? (Hint: Look at other desirable features.)

(1 10.5s) cp 0.039 ( 1 5s ) cp ( s ) 0. 039 3 First derivative Third derivative! 3. Moderate anipulated variable adjustents 4. Robustness to odel isatch To achieve these features, we ust be able to slow down the controller. We chose to include a filter in the feedback path.

D(s) d (s) SP(s) T p (s) MV(s) CP (s) v (s) P (s) - CV(s) (s) - f (s) E (s) f τ f 1 s 1 N f The value of N f is selected so that the product of the controller and filter has a polynoial in s with a denoinator order the nuerator order. The filter is designed to prevent pure derivatives and the filter constant is tuned to achieve robust perforance.

Class exercise: Deterine the structure of the filter for the two possible controllers we just designed for the three tank ixing process. Epirical odel 5.5s 0.039e (1 10.5s) Fundaental odel 0.039 (1 5s) 3 (1 10.5s) cp 0.039 ( 1 5s ) cp ( s ) 0. 039 3

3 ) 5 (1 0.039 ) ( s s ) 10.5 (1 0.039 ) ( 5.5 s e s s Epirical odel Fundaental odel 0.039 ) 10.5 (1 ) ( s s cp 0 039 5 1 3. ) ( ) ( s s cp Class exercise: Deterine the structure of the filter for the two possible controllers we just designed for the three tank ixing process. 1 1 1 ) ( s s f f τ 3 1 1 ) ( s s f f τ

Now, we ust deterine the proper value for the filter tie constant - this is controller tuning again! We know how to set the goals fro experience with PID. CV Dynaic Behavior: Stable, zero offset, iniu IAE MV Dynaic Behavior: daped oscillations and sall fluctuations due to noise. MV can be ore aggressive in early part of transient

We can tune using a siulation and optiization or Process reaction curve Solve the tuning proble. Requires a coputer progra. Apply, and fine tune. 1.5 1 0.5 0-0.5 0 5 101520253035404550 1 0.8 0.6 0.4 0.2 0 5 101520253035404550 v 1 TC v 2 COMBINED DEFINITION OF TUNIN First order with dead tie process odel Noisy easureent signal ± 25% paraeters errors between odel/plant IMC controller: deterine K f Miniize IAE with MV inside bound v 1 TC v 2 K 1 θ 5 τ 5 τ f??? We can develop tuning correlations. These are available for a lead-lag controller, cp (s) f (s).

0.55 IMC Tuning for a typical set of conditions and goals Scaled filter tie constant, τ f /(θτ) 0.50 0.45 0.40 0.35 0.30 0.25 0.20 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 Fraction dead tie, θ/(θτ)

Class exercise: Deterine the filter tie constant for the IMC controller that we just designed for the three tank ixing process. Epirical odel 5.5s 0.039e (1 10.5s) (1 10.5s) cp 0.039 f 1 τ f s 1 1

Class exercise: Deterine the filter tie constant for the IMC controller that we just designed for the three tank ixing process. Epirical odel 5.5s 0.039e (1 10.5s) (1 10.5s) cp 0.039 f τ f 1 s 1 1 1. The controller is lead-lag, so we can use the correlation. 2. θ/(θτ) 5.5/16.0 0.34 3. τ f /(θτ) 0.38 ; τ f (0.38)(15.5) 6.1 inutes 4. See the next slide for the dynaic response. 5. We can fine tune to achieve the desired CV and MV behavior.

F S solvent F A pure A AC 3.5 concentration 3 AC Discuss the control perforance. 2.5 0 10 20 30 40 50 60 70 tie 60 50 valve opening, % 40 30 20 10 v A 0 0 10 20 30 40 50 60 70 tie

Sith Predictor: A sart fellow naed Sith saw the benefit for an explicit odel in the late 1950 s. He invented the Sith Predictor, using the PI controller. The eleents in the box calculate an approxiate inverse D(s) d (s) SP(s) MV(s) C (s) v (s) P (s) - - CV(s) - (s) (s) - E (s) Notes: c (s) is a PI controller; - (s) is the invertible factor.

When do we select an IMC over a PID? 1. Very high fraction dead tie, θ/(θτ) > 0.7. 2. Very strong inverse responses. 3. Cascade priary controller with slow secondary (inner) dynaics 4. Feedforward with disturbance dead tie less than feedback dead tie. See the textbook for further discussion.

CHAPTER 19: IMC CONTROL WORKSHOP 1 Trouble shooting: You have deterined the odel below epirically and have tuned the IMC controller using the correlations. The closedloop perforance is not acceptable. What do you do? 1 Epirical odel K 1 θ 5 τ 5 IMC filter τ f 3.5 concentration anipulated flow 0.8 0.6 0.4 0.2 0-0.2-0.4 0 50 100 150 tie 0-0.5-1 -1.5-2 0 50 100 150 tie

CHAPTER 19: IMC CONTROL WORKSHOP 2 Design features: We want to avoid integral windup. 1. Describe integral windup and why it is undesirable. 2. Modify the structure to provide anti-reset windup. D(s) d (s) SP(s) T p (s) MV(s) CP (s) v (s) P (s) - CV(s) (s) - f (s) E (s)

CHAPTER 19: IMC CONTROL WORKSHOP 3 Design features: We want to obtain good perforance for set point changes and disturbances. The filter below affects disturbances. Introduce a odification that enables us to influence set point responses independently. D(s) d (s) SP(s) T p (s) MV(s) CP (s) v (s) P (s) - CV(s) (s) - f (s) E (s)

CHAPTER 19: IMC CONTROL WORKSHOP 4 Design features: We discussed two reasons why we cannot achieve perfect feedback control. Identify other reasons and explain how they affect the IMC structure design. D(s) d (s) SP(s) T p (s) MV(s) CP (s) v (s) P (s) - CV(s) (s) - f (s) E (s)

When I coplete this chapter, I want to be able to do the following. Recognize that other feedback algoriths are possible Understand the IMC structure and how it provides the essential control features Tune an IMC controller Correctly select between PID and IMC Lot s of iproveent, but we need soe ore study! Read the textbook Review the notes, especially learning goals and workshop Try out the self-study suggestions Naturally, we ll have an assignent!

CHAPTER 19: LEARNIN RESOURCES SITE PC-EDUCATION WEB - Tutorials (Chapter 19) The Textbook, naturally, for any ore exaples. Addition inforation on IMC control is given in the following reference. Brosilow, C. and B. Joseph, Techniques of Model- Based Control, Prentice-Hall, Upper Saddle River, 2002

CHAPTER 19: SUESTIONS FOR SELF-STUDY 1. Discuss the siilarities and differences between IMC and Sith Predictor algoriths 2. Develop the equations that would be solved for a digital ipleentation of the IMC controller for the three-tank ixer. 3. Select a feedback control exaple in the textbook and deterine the IMC tuning for this process. 4. Find a feedback control exaple in the textbook for which IMC is a better choice than PID. 5. Explain how you would ipleent a digital algorith to provide good control for the process in Exaple 19.8 experiencing the flow variation described.