IMC-optimization of a direct reduced iron phenomenological simulator
|
|
- Howard Richardson
- 5 years ago
- Views:
Transcription
1 IMCoptimization of a direct reduced iron phenomenological simulator André Desbiens Andrew A. Shook LOOP (Process Observation and Optimization Laboratory) Department of Electrical and Computer Engineering, Université Laval Pavillon AdrienPouliot, Quebec City (Québec), Canada G1K 7P4 desbiens@gel.ulaval.ca BHP Billiton Minerals Technology Off Vale Street, Shortland, NSW, Australia 237 andrew.a.shook@bhpbilliton.com Abstract This paper deals with the optimization of a direct reduced iron (DRI) phenomenological static simulator. The technique is based on an internal model controller (IMC) designed using a static approximate reduced (even linear) model. The cost function depends only on the reduced model rendering the minimization much easier. Appropriate corrections are calculated by the means of the IMC structure. Even in presence of a large model mismatch, constraints are respected. The convergence is usually much faster than a direct optimization based of the nonlinear process itself. The algorithm is successfully applied to a direct reduced iron phenomenological simulator, allowing to draw some metallurgical conclusions about the guidance of the process. 1 Introduction Several powerful optimization techniques exist. A comprehensive and uptodate description of these methods is detailed in [8]. When the model used to define the behaviour of the process to be optimized is linear, the minimization of the performance index is usually easy, even in presence of constraints. However if the model is nonlinear, more complex algorithms are needed. They are often hard to implement and setting their parameters may sometimes be difficult. The calculation time may also be long. The first aim of this paper is to present a static optimization algorithm which relies on an approximate reduced or even linear model of the nonlinear process. The proposed algorithm is simple to tune and to implement and converges rapidly. The solution may be suboptimal but it remains sufficiently close to the optimum for pratical applications as will be illustrated in this paper and in [1, 2]. Even in presence of a mismatch between the nonlinear process and the linear model, constraints are respected. The convergence is usually faster than a direct optimization based of the nonlinear process. The second objective of the paper is to apply the IMCoptimization to a direct reduced iron (DRI) phenomenological static simulator. The optimization consists in calculating the manipulated variables to obtain the desired metallurgical product properties while respecting production and safety constraints. 2 IMCoptimization The system to be optimized is static and nonlinear. It is represented by an operator H which acts on the inputs U R nu to produce the outputs Y = (HU) R ny. The cost function to be minimized is J = f (Y, U) (1) The optimization is subject to the following constraints V V max W = W eq (2) where V = (P U) R nx and W = (QU) R nw are some states of the static process, which are respectively obtained by applying the static and nonlinear operators P and Q to U. The operators H, P and Q being possibly highly nonlinear, the minimization of (1) subject to (2) may be difficult and time consuming. The IMC structure [3] depicted in Figure 1 facilitates and accelerates the optimization. The solution is obtained by simulating Figure 1 until convergence. The operators H, P and Q are only needed to calculate the corrections ɛ y, ɛ v and ɛ w through the IMC structure. The cost function and the constraints are based on simpler operators H M, P M
2 H M, P M, Q M W eq min J C U V max U e y e v e w (1ρ) 1ρ I U H H M P P M Q Q M Y M V M W M Y V W Under Assumption 1, the optimal solution is then U op = ( H T ΩH Ψ ) 1 ( H T ΩR Λ ) (7) Simulating until convergence Figure 1 using the error weight Ω IMC instead of Ω leads to U IMC = ( H T M Ω IMC H Ψ ) 1 ( H T M Ω IMC R Λ ) (8) The above equations are developed and analyzed only to gain a better understanding of the behaviour of the structure depicted in Figure 1. Comparing (8) to (7) shows that the IMCoptimization leads to an optimal result only if the IMC error weight is selected to respect the following equality Figure 1: IMCoptimization structure H T Ω = H T M Ω IMC (9) and Q M, which may even be linear. The cost function J C at step k of the simulation is J C = f (Y M ɛ y (k 1), U ) (3) and the constraints are V M V max ɛ v (k 1) W M = W eq ɛ w (k 1) (4) The time constant of the lowpass filter is defined by setting the parameter ρ ( ρ < 1). Theorem 1 (Respect of the constraints) If simulating Figure 1 converges to a steadystate when minimizing the objective function (3) subject to the constraints (4), then in the steadystate the constraints (2) are also respected. Proof. The following equations are obtained from Figure 1 ɛ v (k 1) = V (k 1) V M (k 1) ɛ w (k 1) = W (k 1) W M (k 1) Inserting (5) into (4) leads to (2) in steadystate. Assumption 1 The cost function is (5) J =.5 (R HU) T Ω (R HU).5U T ΨU Λ T U (6) where Ω R ny ny and Ψ R nu nu are positive definite weights. The vector Λ R nu is also a weight selected by the user. Also, to help the analysis of Figure 1, it is assumed that it converges to a steadystate, H and H M are matrices of linear gains and there is no constraint (they would be respected if present, see theorem 1). In practice, this is not critical since the weights are sometimes difficult to select and must often be fine tuned by trials and errors. The examples in Section 3.4 illustrate that, even when the process to be optimized is nonlinear, the proposed method leads to results certainly sufficiently close in practice to the optimal solution. Other examples confirming this property are detailed in [1, 2]. The stability (convergence) analysis of the system depicted in Figure 1 is a difficult topic. The tool to analyze the stability of discrete multivariable nonlinear systems is the Tsypkin criterion, which is the discrete analog to the Popov criterion used for analyzing continuous nonlinear systems. However the structural restrictions of even the more recent results [5, 6, 9] make impossible the stability analysis of Figure 1. To apply the Tsypkin criterion and its extensions, it is required to know the nonlinearities sector bounds which is difficult in practice when H consists of phenomenological equations (which is the case for the example presented in this paper). But even worse, H must be diagonal. Consequently, the stability analysis is restricted to the case where H is linear. In that case, the IMC structure is very robust and it can therefore be hoped that nonlinear systems can also be stabilized. For r r linear systems, it was demonstrated [4, 7] that, with an IMC controller, stability can be guaranteed by increasing the filter parameter ρ as long as Re ( [ ]) λ j HH 1 M > j = 1,..., r (1) where λ j [A] denotes the j th eigenvalue of A. Theorem 6.7 in [1] also confirms that a controller with integral action in all channel (which is the case for an IMCstructure) can be made stable even with large model errors by decreasing the performance (increasing ρ), provided that deth/deth M >, i.e. both H and H M have the same sign. Thus, ρ offers a
3 compromise between performance and robustness. To avoid unbounded constraints, it is also required to have detp /detp M > and detq/detq M >. Ore Gas to Recycle Loop When implementing the proposed structure, the simulation can be stopped when U U(k 1) 2 < η n u, where η >. The parameter η is therefore the desired precision for the elements of the solution U, equivalent to the termination tolerance for usual optimization algorithms. R4 R3 R2 R1 Recycle Gas HBI 3 Optimization of the direct reduced iron phenomenological simulator 3.1 Plant and simulator description BHP Billiton s Boodarie Iron direct reduced iron plant is located in Port Hedland, Western Australia. The plant uses FINMET technology to convert beneficiated iron ore fines to metallic iron in a series of pressurized fluidisedbed reactors. The iron ore feed is substantially hematite (Fe 2 O 3 ) particles smaller than 1 mm in diameter. The reactors are arranged in 4 parallel trains, each consisting of four sequential reactors (Figure 2). Within each reactor, the oxygen is removed from the ore by reaction with the gas, which contains large quantities of hydrogen. The gas passes upward through each reactor in the train, accumulating oxygen in the form of gaseous water. The ore passes down through each reactor in series under the influence of gravity. The reactors operate at temperatures between 4 and 8 C and at pressures of MPa. The extent of reduction of the iron ore is determined by a number of factors, including the intrinsic ore reducibility, the reducing gas composition, operating temperature and reactor residence time. The material from the final solids reactor (R1) is continuously formed into briquettes at elevated temperature. The plant is capable of producing 79 tonnes/hr of product from each of the four reactor trains. A fundamentallybased steadystate mathematical simulation of the plant was developed by personnel at BHP Billiton s Newcastle Technology Centre. This model applies laws of conservation of energy and mass together with descriptions of the thermodynamic and kinetic behaviour of the gas and ore to compute the temperatures and compositions occurring within a reactor train. A Gibbs free energy minimization technique is used to compute the gasphase thermodynamic equilibrium. Experimentallydetermined kinetic rates are applied to compute the degree of solids reaction. The mixing within each of the fluidised beds is computed using standard chemical engineering correlations. The composition of the reducing gas (containing H 2, Figure 2: DRI plant H 2 O, CH 4, CO, CO 2, N 2, H 2 S) is such that a large number of interacting gasgas and gassolid reactions are possible. As a result, the simulation is capable of predicting the significant changes to plant operations that occur due to alterations in gasphase equilibria. However, the complexity and nonlinear nature of these calculations mean that repeated simulations consume excessive amounts of computer time, thus rendering the online direct optimization very difficult. 3.2 The variables and the reduced models of the DRI simulator The manipulated variables (U) are the ore flow rate, the pressure at R1 inlet, the residence time of all four reactors, the reducing gas flow rate and temperature and the reformed gas flow rate. Since the recycle loop is closed, the recycle gas flow rate is obtained from a balance with other flow rates and therefore cannot be manipulated. The output of the simulator (Y ) is the percentage of metallization. The constrained variables (V ) are the level of all four reactors and the bed temperature of R1 and R2. In order to have signals with approximately all the same magnitude, they are normalized and the operation points are removed. Since the operators H and P (the DRI simulator itself) are complex and nonlinear, the optimization will rely instead on reduced models H M and P M and the structure depicted in Figure 1. The models were calibrated on data recorded from running the DRI simulator at various operating points. The operator H M predicting the normalized percentage of metallization is as follow y = 9 b i u i i=1 9 b ij u i u j i=1 j i 9 b ijk u i u j u k i=1 j i k j (11) Similar equations define the operator P M used to predict v 1 to v 6. Several bparameters in the reduced models are zero.
4 Percentage of metallization Set point NLop R1 bed temperature x Constraint NLop Percentage of metallization Set point NLop R1 bed temperature Constraint NLop R1 level NLop Ore flow rate NLop R1 level NLop Ore flow rate NLop Figure 3: Example The objective function and the constraints The objective function to be minimized by manipulating u 1 to u 9 is J = (y r) 2 λ 1 u 1 λ 7 u 7 λ 8 u 8 λ 9 u 9 (12) where r is the desired percentage of metallization. The scalars λ 1, λ 7, λ 8 and λ 9 are weights to be selected by the user. Minimizing J therefore corresponds to bringing the percentage of metallization y close to its set point (maintaining the quality of the product) while maximizing the ore flow rate u 1 (maximizing throughput) and minimizing the reducing gas flow rate u 7, the reducing gas temperature u 8 and the reformed gas flow rate u 9 (minimizing production cost). The minimization of J is subject to 9 inequality constraints on U and 1 inequality constraints on V. 3.4 Comparisons with direct optimization The IMCoptimization (IMCop) presented in Section 2 will be compared to two direct optimizations of (1) with respect to (2). The first direct optimization (NLop) uses the same initial guess as the IMCoptimization, which is zero for all variables. The initial guess for the second scheme () is the result of the optimization based on H M and P M, which is in fact the result of the first iteration of the IMCoptimization. To find the optimum, (NLop) and () rely on the Matlab function fmincon with the termination tolerance on U set to The parameter η in the IMCop scheme is set to the same value. No lowpass filter is used, i.e. ρ =. The minimum of (3) with respect to (4) is also found with fmincon. Example 1. The weights are λ 1 =.1 and λ 7 = λ 8 = λ 9 =. Figure 3 shows the stepbystep behaviour of some (normalized with the operating point removed) variables during the IMCoptimization. The lines labelled represent the parameter values of the DRI simulator while those labelled are the re Figure 4: Example 2 duced model states. In all graphs, the first point (iteration ) represent the simulator in openloop at its nominal operating point. The IMC structure starts to be effective only at the next iteration. Note that because all initial conditions are zero, the first iteration consists in doing an optimization based on H M and P M : the percentage of metallization predicted by H M is therefore equal to the set point at that first step. In the following iterations, corrections are calculated through the IMC structure to take into account the model mismatches. Since λ 1 was set to.1, the ore flow rate is increased to its maximum permissible value. To ease comparisons, the results obtained with NLop and are also plotted on the same graphs. The IMCop method needed to run the DRI simulator only 3 times (total calculation time: 28.3 s) while it was respectively run 23 times (191.9 s) and 11 times (625.4 s) for and NLop. All methods give almost equal terminal value of the cost function (12), IMCop being the best by.2% and.3%. All constraints are respected. Example 2. The weights are the same as in Example 1 but the desired percentage of metallization is larger. Again, IMCop needs to run the DRI simulator very few times (IMCop: 6 runs, 33.4 s; : 68 runs, s; : 178 runs, 92.9 s). IMCop leads to the smallest terminal value of the cost function by 1.5% and.4%. Figure 4 depicts some results. 4 Conclusion An interesting characteristic of the proposed IMCoptimization is to show the duality between control and optimization algorithms. The rapidity of convergence for optimization corresponds to the performance for control; divergence of the optimization is similar to unstability of the closedloop system. Hence, the
5 parameter ρ is related to the length of the steps for optimization and to the closedloop time constant for control, therefore allowing a tradeoff between performance and stability. The variable k in Figure 1 can be interpreted either as the optimization step from an optimization point of view or as the time when thinking in terms of control. The main advantages of the proposed IMCoptimization are: It can rely on a simple optimization algorithm and is easy to implement: If H M is linear, the minimum of the cost function can simply be found using quadratic programming. An analytical solution even exists in absence of constraints. If H M is nonlinear but simple, the optimization routine do not have to be robust and complex. For the DRI problem, when using a commercial optimization routine written in C instead of the Matlab function fmincon, the method NLop frequently diverges. Divergence has never occured using the same optimization routine combined with the proposed IMC approach. It is easy to tune: Because of the simplicity of the model appearing in the cost function, the default parameters of the minimization algorithm are usually adequate. The only parameters to tune are η and ρ, both with clear meanings. It converges rapidly: As illustrated in the examples detailed in Section 3.4 and in [1, 2], IMCop usually converges much faster than NLop when they both rely on the same optmization routine. If the behaviour of H differs significantly for various possible operating points, it is recommended to use a different H M for each of them. The IMCoptimization makes it possible to rapidly evaluate different situations according to the weights in the cost function, the desired percentage of metallization and the selected constraints. In most situtations, when the throughput is increased as much as possible, the quality (y) noticeably decreases and the production cost (u 7 u 8 u 9 ) is increased. Trying to minimize the production costs is translated into very large decreases of throughput and product quality. The tests allow to quantify the gains and losses and to determine the best tradeoff. Currently underway, the following step of the project is the online implementation of the IMCoptimization at the Port Hedland plant. Every time when needed, the simulator will be put in a state similar to the plant using online data. Then, the optimization will produce guidance advices to the operators by minimizing a cost function representative of the aimed metallugical product properties and constraints to be respected. Acknowledgements The authors are grateful to CIDAC (Centre for Integrated Dynamics and Control at the University of Newcastle), NSERC (Natural Science and Engineering Research Council of Canada) and BHP Billiton for their financial support and to BHP Billiton for its authorization to publish the results. References [1] D. Pomerleau, A. Desbiens, G.W. Barton, Realtime optimization of an extrusion cooking process using a first principles model, sumitted to CCA 23 Conf. on Control Appl., (Istanbul, Turkey), (23). [2] D. Pomerleau, A. Desbiens and D. Hodouin, Optimization of a simulated ironoxide pellets induration furnace, accepted for publication, 11th Mediterranean Conf. on Control, (Rhodes, Greece), 23. [3] C.E. Garcia and M. Morari, Internal model control 1. A unifying review and some new results, Ind. Engng Chem. Process Des. Dev., vol. 21, pp , [4] C.E. Garcia and M. Morari, Internal model control 2. Design procedure for multivariable systems, Ind. Engng Chem. Process Des. Dev., vol. 24, pp , [5] V. Kapila and W. Haddad, A multivariable extension of the Tsypkin criterion using a Lyapunovfunction approach, IEEE Trans. Autom. Control., vol. 41(1), pp , [6] M. Larsen and P. Kokotović, A brief look at the Tsypkin criterion: from analysis to design, Int. J. Adapt. Control Signal Process., vol. 15, pp , 21. [7] M. Morari, Response to comments on Internal model control 2. Design procedure for multivariable systems, Ind. Engng Chem. Res., vol. 26, pp. 633, [8] J. Nocedal and S.J. Wright, Numerical optimization, Springer, New York, [9] P. Park and S.W. Kim, A revisited Tsypkin criterion for discrete nonlinear Lur e systems with monotonic sectorrestrictions, Automatica, vol. 34(11), pp , [1] S. Skogestad and I. Postlethwaite, Multivariable feedback control: Analysis and design, John Wiley & Sons, England, 1996.
NONLINEAR PREDICTIVE CONTROL OF A ROUGHER FLOTATION UNIT USING LOCAL MODELS
NONLINEAR PREDICIVE CONROL OF A ROUGHER FLOAION UNI USING LOCAL MODELS A. Desbiens *, D. Hodouin ** and M. Mailloux * GRAIIM (Groupe de Recherche sur les Applications de l Informatique à l Industrie Minérale)
More informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
More informationIncorporating Feedforward Action into Self-optimizing Control Policies
Incorporating Feedforward Action into Self-optimizing Control Policies Lia Maisarah Umar, Yi Cao and Vinay Kariwala School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore
More informationMULTILOOP CONTROL APPLIED TO INTEGRATOR MIMO. PROCESSES. A Preliminary Study
MULTILOOP CONTROL APPLIED TO INTEGRATOR MIMO PROCESSES. A Preliminary Study Eduardo J. Adam 1,2*, Carlos J. Valsecchi 2 1 Instituto de Desarrollo Tecnológico para la Industria Química (INTEC) (Universidad
More informationDESIGN OF AN ON-LINE TITRATOR FOR NONLINEAR ph CONTROL
DESIGN OF AN ON-LINE TITRATOR FOR NONLINEAR CONTROL Alex D. Kalafatis Liuping Wang William R. Cluett AspenTech, Toronto, Canada School of Electrical & Computer Engineering, RMIT University, Melbourne,
More informationOn an internal multimodel control for nonlinear multivariable systems - A comparative study
On an internal multimodel control for nonlinear multivariable systems A comparative study Nahla Touati Karmani Dhaou Soudani Mongi Naceur Mohamed Benrejeb Abstract An internal multimodel control designed
More informationInternal Model Control of A Class of Continuous Linear Underactuated Systems
Internal Model Control of A Class of Continuous Linear Underactuated Systems Asma Mezzi Tunis El Manar University, Automatic Control Research Laboratory, LA.R.A, National Engineering School of Tunis (ENIT),
More informationA unified double-loop multi-scale control strategy for NMP integrating-unstable systems
Home Search Collections Journals About Contact us My IOPscience A unified double-loop multi-scale control strategy for NMP integrating-unstable systems This content has been downloaded from IOPscience.
More informationAsymptotic behaviour of Toeplitz matrix in multi-input multi-output model predictive control
23 European Control Conference ECC) July 7-9, 23, Zürich, Switzerland Asymptotic behaviour of Toeplitz matrix in multi-input multi-output model predictive control Quang Tran, Leyla Özkan, Jobert Ludlage
More informationONE can design optical filters using different filter architectures.
JOURNAL OF LIGHTWAVE TECHNOLOGY, VOL. 28, NO. 23, DECEMBER 1, 2010 3463 Comparison of Cascade, Lattice, and Parallel Filter Architectures Rohit Patnaik, Vivek Vandrasi, Christi K. Madsen, Ali A. Eftekhar,
More informationLMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System
LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz
More informationReduction of Random Variables in Structural Reliability Analysis
Reduction of Random Variables in Structural Reliability Analysis S. Adhikari and R. S. Langley Department of Engineering University of Cambridge Trumpington Street Cambridge CB2 1PZ (U.K.) February 21,
More informationIMPROVED CONSTRAINED CASCADE CONTROL FOR PARALLEL PROCESSES. Richard Lestage, André Pomerleau and André Desbiens
IMPROVED CONSTRAINED CASCADE CONTROL FOR PARALLEL PROCESSES Richard Lestage, André Pomerleau and André Desbiens Groupe de Recherche sur les Applications de l Informatique à l Industrie Minérale(GRAIIM),
More informationHANKEL-NORM BASED INTERACTION MEASURE FOR INPUT-OUTPUT PAIRING
Copyright 2002 IFAC 15th Triennial World Congress, Barcelona, Spain HANKEL-NORM BASED INTERACTION MEASURE FOR INPUT-OUTPUT PAIRING Björn Wittenmark Department of Automatic Control Lund Institute of Technology
More informationRobust control for a multi-stage evaporation plant in the presence of uncertainties
Preprint 11th IFAC Symposium on Dynamics and Control of Process Systems including Biosystems June 6-8 16. NTNU Trondheim Norway Robust control for a multi-stage evaporation plant in the presence of uncertainties
More informationRobust cross-directional control of paper making machines with saturating actuators
Robust cross-directional control of paper making machines with saturating actuators J. A. G. Akkermans, A. G. Wills and W. P. Heath March 2, 2004 1 Introduction Preface Technical Report EE03044 In this
More informationFeedback stabilisation with positive control of dissipative compartmental systems
Feedback stabilisation with positive control of dissipative compartmental systems G. Bastin and A. Provost Centre for Systems Engineering and Applied Mechanics (CESAME Université Catholique de Louvain
More informationRiccati difference equations to non linear extended Kalman filter constraints
International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R
More informationCourse on Model Predictive Control Part II Linear MPC design
Course on Model Predictive Control Part II Linear MPC design Gabriele Pannocchia Department of Chemical Engineering, University of Pisa, Italy Email: g.pannocchia@diccism.unipi.it Facoltà di Ingegneria,
More informationCircle Criterion in Linear Control Design
8 American Control Conference Westin Seattle Hotel Seattle Washington USA June -3 8 ThC8. Circle Criterion in Linear Control Design Marcel Heertjes and Maarten Steinbuch Abstract This paper presents a
More informationQuantitative Feedback Theory based Controller Design of an Unstable System
Quantitative Feedback Theory based Controller Design of an Unstable System Chandrima Roy Department of E.C.E, Assistant Professor Heritage Institute of Technology, Kolkata, WB Kalyankumar Datta Department
More informationReal-Time Implementation of Nonlinear Predictive Control
Real-Time Implementation of Nonlinear Predictive Control Michael A. Henson Department of Chemical Engineering University of Massachusetts Amherst, MA WebCAST October 2, 2008 1 Outline Limitations of linear
More informationComparative study of three practical IMC algorithms with inner controller of first and second order
Journal of Electrical Engineering, Electronics, Control and Computer Science JEEECCS, Volume 2, Issue 4, pages 2-28, 206 Comparative study of three practical IMC algorithms with inner controller of first
More informationThermodynamic and Stochiometric Principles in Materials Balance
Thermodynamic and Stochiometric Principles in Materials Balance Typical metallurgical engineering problems based on materials and energy balance NiO is reduced in an open atmosphere furnace by excess carbon
More informationIan G. Horn, Jeffery R. Arulandu, Christopher J. Gombas, Jeremy G. VanAntwerp, and Richard D. Braatz*
Ind. Eng. Chem. Res. 996, 35, 3437-344 3437 PROCESS DESIGN AND CONTROL Improved Filter Design in Internal Model Control Ian G. Horn, Jeffery R. Arulandu, Christopher J. Gombas, Jeremy G. VanAntwerp, and
More informationGAIN SCHEDULING CONTROL WITH MULTI-LOOP PID FOR 2- DOF ARM ROBOT TRAJECTORY CONTROL
GAIN SCHEDULING CONTROL WITH MULTI-LOOP PID FOR 2- DOF ARM ROBOT TRAJECTORY CONTROL 1 KHALED M. HELAL, 2 MOSTAFA R.A. ATIA, 3 MOHAMED I. ABU EL-SEBAH 1, 2 Mechanical Engineering Department ARAB ACADEMY
More informationGain Scheduling Control with Multi-loop PID for 2-DOF Arm Robot Trajectory Control
Gain Scheduling Control with Multi-loop PID for 2-DOF Arm Robot Trajectory Control Khaled M. Helal, 2 Mostafa R.A. Atia, 3 Mohamed I. Abu El-Sebah, 2 Mechanical Engineering Department ARAB ACADEMY FOR
More informationFormal verification of One Dimensional Time Triggered Velocity PID Controllers Kenneth Payson 12/09/14
Formal verification of One Dimensional Time Triggered Velocity PID Controllers 12/09/14 1: Abstract This paper provides a formal proof of the safety of a time triggered velocity PID controller that are
More informationDynamic Operability for the Calculation of Transient Output Constraints for Non-Square Linear Model Predictive Controllers
Dynamic Operability for the Calculation of Transient Output Constraints for Non-Square Linear Model Predictive Controllers Fernando V. Lima and Christos Georgakis* Department of Chemical and Biological
More informationROBUSTNESS COMPARISON OF CONTROL SYSTEMS FOR A NUCLEAR POWER PLANT
Control 004, University of Bath, UK, September 004 ROBUSTNESS COMPARISON OF CONTROL SYSTEMS FOR A NUCLEAR POWER PLANT L Ding*, A Bradshaw, C J Taylor Lancaster University, UK. * l.ding@email.com Fax: 0604
More informationADAPTIVE TEMPERATURE CONTROL IN CONTINUOUS STIRRED TANK REACTOR
INTERNATIONAL JOURNAL OF ELECTRICAL ENGINEERING & TECHNOLOGY (IJEET) International Journal of Electrical Engineering and Technology (IJEET), ISSN 0976 6545(Print), ISSN 0976 6545(Print) ISSN 0976 6553(Online)
More informationLinear-Quadratic Optimal Control: Full-State Feedback
Chapter 4 Linear-Quadratic Optimal Control: Full-State Feedback 1 Linear quadratic optimization is a basic method for designing controllers for linear (and often nonlinear) dynamical systems and is actually
More informationROBUST STABLE NONLINEAR CONTROL AND DESIGN OF A CSTR IN A LARGE OPERATING RANGE. Johannes Gerhard, Martin Mönnigmann, Wolfgang Marquardt
ROBUST STABLE NONLINEAR CONTROL AND DESIGN OF A CSTR IN A LARGE OPERATING RANGE Johannes Gerhard, Martin Mönnigmann, Wolfgang Marquardt Lehrstuhl für Prozesstechnik, RWTH Aachen Turmstr. 46, D-5264 Aachen,
More informationA METHOD OF ADAPTATION BETWEEN STEEPEST- DESCENT AND NEWTON S ALGORITHM FOR MULTI- CHANNEL ACTIVE CONTROL OF TONAL NOISE AND VIBRATION
A METHOD OF ADAPTATION BETWEEN STEEPEST- DESCENT AND NEWTON S ALGORITHM FOR MULTI- CHANNEL ACTIVE CONTROL OF TONAL NOISE AND VIBRATION Jordan Cheer and Stephen Daley Institute of Sound and Vibration Research,
More informationRobust Internal Model Control for Impulse Elimination of Singular Systems
International Journal of Control Science and Engineering ; (): -7 DOI:.59/j.control.. Robust Internal Model Control for Impulse Elimination of Singular Systems M. M. Share Pasandand *, H. D. Taghirad Department
More informationDesign and Stability Analysis of Single-Input Fuzzy Logic Controller
IEEE TRANSACTIONS ON SYSTEMS, MAN, AND CYBERNETICS PART B: CYBERNETICS, VOL. 30, NO. 2, APRIL 2000 303 Design and Stability Analysis of Single-Input Fuzzy Logic Controller Byung-Jae Choi, Seong-Woo Kwak,
More informationOptimization 2. CS5240 Theoretical Foundations in Multimedia. Leow Wee Kheng
Optimization 2 CS5240 Theoretical Foundations in Multimedia Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (NUS) Optimization 2 1 / 38
More informationDirect and Indirect Gradient Control for Static Optimisation
International Journal of Automation and Computing 1 (2005) 60-66 Direct and Indirect Gradient Control for Static Optimisation Yi Cao School of Engineering, Cranfield University, UK Abstract: Static self-optimising
More informationROBUST CONSTRAINED REGULATORS FOR UNCERTAIN LINEAR SYSTEMS
ROBUST CONSTRAINED REGULATORS FOR UNCERTAIN LINEAR SYSTEMS Jean-Claude HENNET Eugênio B. CASTELAN Abstract The purpose of this paper is to combine several control requirements in the same regulator design
More informationMachine Learning Applied to 3-D Reservoir Simulation
Machine Learning Applied to 3-D Reservoir Simulation Marco A. Cardoso 1 Introduction The optimization of subsurface flow processes is important for many applications including oil field operations and
More information80DB A 40DB 0DB DB
Stability Analysis On the Nichols Chart and Its Application in QFT Wenhua Chen and Donald J. Ballance Centre for Systems & Control Department of Mechanical Engineering University of Glasgow Glasgow G12
More informationCOMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE
COMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE Rafael ÅMAN*, Heikki HANDROOS*, Pasi KORKEALAAKSO** and Asko ROUVINEN** * Laboratory
More informationA NUMERICAL METHOD TO SOLVE A QUADRATIC CONSTRAINED MAXIMIZATION
A NUMERICAL METHOD TO SOLVE A QUADRATIC CONSTRAINED MAXIMIZATION ALIREZA ESNA ASHARI, RAMINE NIKOUKHAH, AND STEPHEN L. CAMPBELL Abstract. The problem of maximizing a quadratic function subject to an ellipsoidal
More informationRobust Adaptive MPC for Systems with Exogeneous Disturbances
Robust Adaptive MPC for Systems with Exogeneous Disturbances V. Adetola M. Guay Department of Chemical Engineering, Queen s University, Kingston, Ontario, Canada (e-mail: martin.guay@chee.queensu.ca) Abstract:
More informationx 2 x n r n J(x + t(x x ))(x x )dt. For warming-up we start with methods for solving a single equation of one variable.
Maria Cameron 1. Fixed point methods for solving nonlinear equations We address the problem of solving an equation of the form (1) r(x) = 0, where F (x) : R n R n is a vector-function. Eq. (1) can be written
More informationPROPORTIONAL-Integral-Derivative (PID) controllers
Multiple Model and Neural based Adaptive Multi-loop PID Controller for a CSTR Process R.Vinodha S. Abraham Lincoln and J. Prakash Abstract Multi-loop (De-centralized) Proportional-Integral- Derivative
More informationIterative Controller Tuning Using Bode s Integrals
Iterative Controller Tuning Using Bode s Integrals A. Karimi, D. Garcia and R. Longchamp Laboratoire d automatique, École Polytechnique Fédérale de Lausanne (EPFL), 05 Lausanne, Switzerland. email: alireza.karimi@epfl.ch
More informationOptimizing Control of Hot Blast Stoves in Staggered Parallel Operation
Proceedings of the 17th World Congress The International Federation of Automatic Control Optimizing Control of Hot Blast Stoves in Staggered Parallel Operation Akın Şahin and Manfred Morari Automatic Control
More informationNEW SUPERVISORY CONTROL USING CONTROL-RELEVANT SWITCHING
NEW SUPERVISORY CONTROL USING CONTROL-RELEVANT SWITCHING Tae-Woong Yoon, Jung-Su Kim Dept. of Electrical Engineering. Korea University, Anam-dong 5-ga Seongbuk-gu 36-73, Seoul, Korea, twy@korea.ac.kr,
More informationCHAPTER 6 CLOSED LOOP STUDIES
180 CHAPTER 6 CLOSED LOOP STUDIES Improvement of closed-loop performance needs proper tuning of controller parameters that requires process model structure and the estimation of respective parameters which
More informationAADECA 2006 XXº Congreso Argentino de Control Automático ITERATIVE LEARNING CONTROL APPLIED TO NONLINEAR BATCH REACTOR E.J.
ITERATIVE LEARNING CONTROL APPLIED TO NONLINEAR BATCH REACTOR E.J. ADAM (1) Institute of Technological Development for the Chemical Industry (INTEC), CONICET Universidad Nacional del Litoral (UNL). Güemes
More informationLecture XI. Approximating the Invariant Distribution
Lecture XI Approximating the Invariant Distribution Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Invariant Distribution p. 1 /24 SS Equilibrium in the Aiyagari model G.
More informationTwo-Link Flexible Manipulator Control Using Sliding Mode Control Based Linear Matrix Inequality
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Two-Link Flexible Manipulator Control Using Sliding Mode Control Based Linear Matrix Inequality To cite this article: Zulfatman
More informationLyapunov Stability of Linear Predictor Feedback for Distributed Input Delays
IEEE TRANSACTIONS ON AUTOMATIC CONTROL VOL. 56 NO. 3 MARCH 2011 655 Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays Nikolaos Bekiaris-Liberis Miroslav Krstic In this case system
More informationNumerical Methods. Root Finding
Numerical Methods Solving Non Linear 1-Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More informationAdvanced Aerospace Control. Marco Lovera Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano
Advanced Aerospace Control Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano ICT for control systems engineering School of Industrial and Information Engineering Aeronautical Engineering
More informationDesign of Decentralised PI Controller using Model Reference Adaptive Control for Quadruple Tank Process
Design of Decentralised PI Controller using Model Reference Adaptive Control for Quadruple Tank Process D.Angeline Vijula #, Dr.N.Devarajan * # Electronics and Instrumentation Engineering Sri Ramakrishna
More informationA Definition for Plantwide Controllability. Process Flexibility
A Definition for Plantwide Controllability Surya Kiran Chodavarapu and Alex Zheng Department of Chemical Engineering University of Massachusetts Amherst, MA 01003 Abstract Chemical process synthesis typically
More informationSparse, stable gene regulatory network recovery via convex optimization
Sparse, stable gene regulatory network recovery via convex optimization Arwen Meister June, 11 Gene regulatory networks Gene expression regulation allows cells to control protein levels in order to live
More informationFAULT-TOLERANT CONTROL OF CHEMICAL PROCESS SYSTEMS USING COMMUNICATION NETWORKS. Nael H. El-Farra, Adiwinata Gani & Panagiotis D.
FAULT-TOLERANT CONTROL OF CHEMICAL PROCESS SYSTEMS USING COMMUNICATION NETWORKS Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides Department of Chemical Engineering University of California,
More informationOnline monitoring of MPC disturbance models using closed-loop data
Online monitoring of MPC disturbance models using closed-loop data Brian J. Odelson and James B. Rawlings Department of Chemical Engineering University of Wisconsin-Madison Online Optimization Based Identification
More informationCONTROLLER PERFORMANCE ASSESSMENT IN SET POINT TRACKING AND REGULATORY CONTROL
ADCHEM 2, Pisa Italy June 14-16 th 2 CONTROLLER PERFORMANCE ASSESSMENT IN SET POINT TRACKING AND REGULATORY CONTROL N.F. Thornhill *, S.L. Shah + and B. Huang + * Department of Electronic and Electrical
More informationPassivity-based Adaptive Inventory Control
Joint 48th IEEE Conference on Decision and Control and 28th Chinese Control Conference Shanghai, P.R. China, December 6-8, 29 ThB.2 Passivity-based Adaptive Inventory Control Keyu Li, Kwong Ho Chan and
More informationCHEMICAL ENGINEERING LABORATORY CHEG 4137W/4139W. Reaction Kinetics Saponification of Isopropyl Acetate with Sodium Hydroxide
CHEMICAL ENGINEERING LABORATORY CHEG 4137W/4139W Reaction Kinetics Saponification of Isopropyl Acetate with Sodium Hydroxide Objective: The purpose of this experiment is to examine and determine the reaction
More informationTWO- PHASE APPROACH TO DESIGN ROBUST CONTROLLER FOR UNCERTAIN INTERVAL SYSTEM USING GENETIC ALGORITHM
International Journal of Electrical and Electronics Engineering Research (IJEEER) ISSN:2250-155X Vol.2, Issue 2 June 2012 27-38 TJPRC Pvt. Ltd., TWO- PHASE APPROACH TO DESIGN ROBUST CONTROLLER FOR UNCERTAIN
More informationModel-based PID tuning for high-order processes: when to approximate
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 25 Seville, Spain, December 2-5, 25 ThB5. Model-based PID tuning for high-order processes: when to approximate
More informationProblem Set III Quantization
Problem Set III Quantization Christopher Tsai Problem #2.1 Lloyd-Max Quantizer To train both the Lloyd-Max Quantizer and our Entropy-Constrained Quantizer, we employ the following training set of images,
More informationCHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER
114 CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 5.1 INTRODUCTION Robust control is a branch of control theory that explicitly deals with uncertainty in its approach to controller design. It also refers
More informationIntroduction to Model Predictive Control. Dipartimento di Elettronica e Informazione
Introduction to Model Predictive Control Riccardo Scattolini Riccardo Scattolini Dipartimento di Elettronica e Informazione Finite horizon optimal control 2 Consider the system At time k we want to compute
More informationA Method for PID Controller Tuning Using Nonlinear Control Techniques*
A Method for PID Controller Tuning Using Nonlinear Control Techniques* Prashant Mhaskar, Nael H. El-Farra and Panagiotis D. Christofides Department of Chemical Engineering University of California, Los
More informationIMC based automatic tuning method for PID controllers in a Smith predictor configuration
Computers and Chemical Engineering 28 (2004) 281 290 IMC based automatic tuning method for PID controllers in a Smith predictor configuration Ibrahim Kaya Department of Electrical and Electronics Engineering,
More informationEconomic Model Predictive Control Historical Perspective and Recent Developments and Industrial Examples
1 Economic Model Predictive Control Historical Perspective and Recent Developments and Industrial Examples Public Trial Lecture Candidate: Vinicius de Oliveira Department of Chemical Engineering, Faculty
More informationNavigation and Obstacle Avoidance via Backstepping for Mechanical Systems with Drift in the Closed Loop
Navigation and Obstacle Avoidance via Backstepping for Mechanical Systems with Drift in the Closed Loop Jan Maximilian Montenbruck, Mathias Bürger, Frank Allgöwer Abstract We study backstepping controllers
More informationImproved Damped Quasi-Newton Methods for Unconstrained Optimization
Improved Damped Quasi-Newton Methods for Unconstrained Optimization Mehiddin Al-Baali and Lucio Grandinetti August 2015 Abstract Recently, Al-Baali (2014) has extended the damped-technique in the modified
More informationMAGNETITE OXIDATION EXAMPLE
HSC Chemistry 7.0 1 MAGNETITE OXIDATION EXAMPLE Pelletized magnetite (Fe 3 O 4 ) ore may be oxidized to hematite (Fe 2 O 3 ) in shaft furnace. Typical magnetite content in ore is some 95%. Oxidation is
More informationGlobal stabilization of feedforward systems with exponentially unstable Jacobian linearization
Global stabilization of feedforward systems with exponentially unstable Jacobian linearization F Grognard, R Sepulchre, G Bastin Center for Systems Engineering and Applied Mechanics Université catholique
More informationAugust Progress Report
PATH PREDICTION FOR AN EARTH-BASED DEMONSTRATION BALLOON FLIGHT DANIEL BEYLKIN Mentor: Jerrold Marsden Co-Mentors: Claire Newman and Philip Du Toit August Progress Report. Progress.. Discrete Mechanics
More informationCHAPTER 2: QUADRATIC PROGRAMMING
CHAPTER 2: QUADRATIC PROGRAMMING Overview Quadratic programming (QP) problems are characterized by objective functions that are quadratic in the design variables, and linear constraints. In this sense,
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 12: Multivariable Control of Robotic Manipulators Part II
MCE/EEC 647/747: Robot Dynamics and Control Lecture 12: Multivariable Control of Robotic Manipulators Part II Reading: SHV Ch.8 Mechanical Engineering Hanz Richter, PhD MCE647 p.1/14 Robust vs. Adaptive
More informationLearning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System
Learning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System Ugo Rosolia Francesco Borrelli University of California at Berkeley, Berkeley, CA 94701, USA
More information3.1 Overview 3.2 Process and control-loop interactions
3. Multivariable 3.1 Overview 3.2 and control-loop interactions 3.2.1 Interaction analysis 3.2.2 Closed-loop stability 3.3 Decoupling control 3.3.1 Basic design principle 3.3.2 Complete decoupling 3.3.3
More information9. Two-Degrees-of-Freedom Design
9. Two-Degrees-of-Freedom Design In some feedback schemes we have additional degrees-offreedom outside the feedback path. For example, feed forwarding known disturbance signals or reference signals. In
More informationProcess Control, 3P4 Assignment 6
Process Control, 3P4 Assignment 6 Kevin Dunn, kevin.dunn@mcmaster.ca Due date: 28 March 204 This assignment gives you practice with cascade control and feedforward control. Question [0 = 6 + 4] The outlet
More informationStatic Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
More informationCHAPTER 10 Zeros of Functions
CHAPTER 10 Zeros of Functions An important part of the maths syllabus in secondary school is equation solving. This is important for the simple reason that equations are important a wide range of problems
More informationDESIGN AND CONTROL OF BUTYL ACRYLATE REACTIVE DISTILLATION COLUMN SYSTEM. I-Lung Chien and Kai-Luen Zeng
DESIGN AND CONTROL OF BUTYL ACRYLATE REACTIVE DISTILLATION COLUMN SYSTEM I-Lung Chien and Kai-Luen Zeng Department of Chemical Engineering, National Taiwan University of Science and Technology, Taipei
More informationON CHATTERING-FREE DISCRETE-TIME SLIDING MODE CONTROL DESIGN. Seung-Hi Lee
ON CHATTERING-FREE DISCRETE-TIME SLIDING MODE CONTROL DESIGN Seung-Hi Lee Samsung Advanced Institute of Technology, Suwon, KOREA shl@saitsamsungcokr Abstract: A sliding mode control method is presented
More informationControl Configuration Selection for Multivariable Descriptor Systems
Control Configuration Selection for Multivariable Descriptor Systems Hamid Reza Shaker and Jakob Stoustrup Abstract Control configuration selection is the procedure of choosing the appropriate input and
More informationProcess Model Formulation and Solution, 3E4
Process Model Formulation and Solution, 3E4 Section B: Linear Algebraic Equations Instructor: Kevin Dunn dunnkg@mcmasterca Department of Chemical Engineering Course notes: Dr Benoît Chachuat 06 October
More informationM. LeBreux, M. Désilets & M. Lacroix Faculté de génie, Université de Sherbrooke, Canada. Abstract. 1 Introduction
Advanced Computational Methods and Experiments in Heat Transfer XIII 517 Is the performance of a virtual sensor employed for the prediction of the ledge thickness inside a metallurgical reactor affected
More informationControl Configuration Selection for Multivariable Descriptor Systems Shaker, Hamid Reza; Stoustrup, Jakob
Aalborg Universitet Control Configuration Selection for Multivariable Descriptor Systems Shaker, Hamid Reza; Stoustrup, Jakob Published in: 2012 American Control Conference (ACC) Publication date: 2012
More informationNonlinear ph Control Using a Three Parameter Model
130 ICASE: The Institute of Control, Automation and Systems Engineers, KOREA Vol. 2, No. 2, June, 2000 Nonlinear ph Control Using a Three Parameter Model Jietae Lee and Ho-Cheol Park Abstract: A two parameter
More informationPrimal-Dual Interior-Point Methods for Linear Programming based on Newton s Method
Primal-Dual Interior-Point Methods for Linear Programming based on Newton s Method Robert M. Freund March, 2004 2004 Massachusetts Institute of Technology. The Problem The logarithmic barrier approach
More informationTime scale separation and the link between open-loop and closed-loop dynamics
Time scale separation and the link between open-loop and closed-loop dynamics Antonio Araújo a Michael Baldea b Sigurd Skogestad a,1 Prodromos Daoutidis b a Department of Chemical Engineering Norwegian
More information(Continued on next page)
(Continued on next page) 18.2 Roots of Stability Nyquist Criterion 87 e(s) 1 S(s) = =, r(s) 1 + P (s)c(s) where P (s) represents the plant transfer function, and C(s) the compensator. The closedloop characteristic
More informationChemical Reactions and Kinetics of the Carbon Monoxide Coupling in the Presence of Hydrogen
Journal of Natural Gas Chemistry 11(2002)145 150 Chemical Reactions and Kinetics of the Carbon Monoxide Coupling in the Presence of Hydrogen Fandong Meng 1,2, Genhui Xu 1, Zhenhua Li 1, Pa Du 1 1. State
More informationSTABILITY AND STABILIZATION OF A CLASS OF NONLINEAR SYSTEMS WITH SATURATING ACTUATORS. Eugênio B. Castelan,1 Sophie Tarbouriech Isabelle Queinnec
STABILITY AND STABILIZATION OF A CLASS OF NONLINEAR SYSTEMS WITH SATURATING ACTUATORS Eugênio B. Castelan,1 Sophie Tarbouriech Isabelle Queinnec DAS-CTC-UFSC P.O. Box 476, 88040-900 Florianópolis, SC,
More informationROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ]
ROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ] General control configuration with uncertainty [8.1] For our robustness analysis we use a system representation in which the uncertain perturbations are
More informationarxiv: v1 [cs.sy] 2 Oct 2018
Non-linear Model Predictive Control of Conically Shaped Liquid Storage Tanks arxiv:1810.01119v1 [cs.sy] 2 Oct 2018 5 10 Abstract Martin Klaučo, L uboš Čirka Slovak University of Technology in Bratislava,
More informationPID control of FOPDT plants with dominant dead time based on the modulus optimum criterion
Archives of Control Sciences Volume 6LXII, 016 No. 1, pages 5 17 PID control of FOPDT plants with dominant dead time based on the modulus optimum criterion JAN CVEJN The modulus optimum MO criterion can
More information