Robust Control of Distributed Parameter Systems with Demonstration in Casting Technology and MATLAB/Simulink/DPS Blockset Software Support

© 2012 Belavý et al., licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Robust Control of Distributed Parameter Systems with Demonstration in Casting Technology and MATLAB/Simulink/DPS Blockset Software Support


Introduction
Most of the dynamical systems analysed in engineering practice have the dynamics, which depends on both position and time. Such systems are classified as distributed parameter systems (DPS). The time-space coupled nature of the DPS is usually mathematically described by partial differential equations (PDE) as infinite-dimensional systems. However, from point of view of implementation of DPS control in technological practice, where a finite number of sensors and actuators for practical sensing and control is at disposal, such infinite-dimensional systems need to be approximated by finite-dimensional systems. There are many dimension reduction methods, which can be used to solve this problem.
In the first mathematical foundations of DPS control, analytical solutions of the underlying PDE have been used (Butkovskij, 1965;Lions, 1971;Wang, 1964). That is the decomposition of dynamics into time and space components based on the eigenfunctions of the PDE. Continuous and approximation theories aimed to control of parabolic systems presents monograph (Lasiecka & Triggiani, 2000). Methodical approach from the view of time-space separation with model reduction is presented in (Li & Qi, 2010). Variety of transfer functions for systems described by PDE are illustrated by means of several examples in (Curtain & Morris, 2009). Well-known reduction methods based on finite difference method (FDM), or finite element method (FEM), spectral method require an accurate nominal PDE model and usually lead to a high-order model, which requires unpractical high-order controller.
An engineering approach for the control of DPS is being developed since the eighties of the last century (Hulkó et al., 1981(Hulkó et al., , 1987(Hulkó et al., , 1998(Hulkó et al., , 2009a(Hulkó et al., , 2009b. In the field of lumped parameters system (LPS) control, where the state/output quantities x(t)/y(t) -parameters are given as finite dimensional vectors, the actuator together with the controlled plant make up a controlled LPS. In this sense the actuators and the controlled plant as a DPS create a controlled lumped-input and distributed-parameter-output system (LDS).
In this chapter the decomposition of dynamics of controlled LDS into time and space components is introduced. Based on this decomposition a methodical framework of control synthesis decomposition into space and time tasks will be presented. In the space domain, approximation problems are solved. In the time domain, synthesis of control is performed by lumped parameter control loops, where robust controllers are used.
The casting technology is a typical case of the DPS. There in order to obtain the desired solidification structure, the casting process requires a specific temperature field of the mould, which is defined on complex-shape 3D definition domain. Modelling, simulation and evaluation of real-time experiments in this area is now widely accepted as an important tool in product design and process development to improve productivity and casting quality. For analysis of the casting process dynamics as DPS, especially temperature fields in the casting mould and control synthesis purposes, the benchmark casting plant with steel mould of complex-shape was designed at Faculty of Mechanical Engineering STU in Bratislava.
The main emphasis of this chapter is to present an engineering approach for the robust control of DPS with demonstration in the casting technology along with software support in the MATLAB & Simulink programming environment. This approach opens a wide space for novel applications of the toolboxes and blocksets of the MATLAB & Simulink software environment. For the software support of modelling, control and design of DPS, given on 1D-3D definition domains, the Distributed Parameter Systems Blockset for MATLAB & Simulink (DPS Blockset), a Third-Party Product of The MathWorks Companywww.mathworks.com/products/connections/ has been developed at the Institute of Automation, Measurement and Applied Informatics, Faculty of Mechanical Engineering STU, (Hulkó et al. 2003(Hulkó et al. -2010. Also a web portal named Distributed Parameter Systems Control -www.dpscontrol.sk has been created for those, who are interested in solving problems of DPS control (Hulkó et al., 2003(Hulkó et al., -2007. This web portal contains application examples from different areas of engineering practice, such as the control of technological and manufacturing processes, mechatronic structures, groundwater remediation, etc. In addition, this web portal offers the demo version of the DPS Blockset with the Tutorial, Show, Demos and DPS Wizard for download, along with the Interactive Control service for the interactive solution of model control problems via the Internet. In this chapter, for the control synthesis purpose, LDS models of temperature fields in the casting mould were created by means of evaluation of real-time experiments. Robust control synthesis based on internal model control (IMC) structure in the time domain has been done. Designed robust controllers were used for the robust control of preheating casting mould in the real-time experiment in accordance to casting technology requirements. Identification, uncertainty analysis of the models, robust control synthesis and experiments were performed with the software support of DPS Blockset and toolboxes of the MATLAB & Simulink, especially the System Identification Toolbox, Control Systems Toolbox, Robust Control Toolbox, Optimization Toolbox, System Identification Toolbox, Real-Time Windows Target and Simulink Design Optimization.

LDS/HLDS representation of DPS
In general, DPS are systems whose state or output quantities,     , , , / , , , X x y z t Y x y z t are distributed quantities or fields of quantities, where   , , x y z are spatial coordinates in 3D. These systems are often considered as systems whose dynamics is described by PDE. In the inputoutput relation, PDE define distributed-input/distributed-output systems (DDS) between distributed input,   , , , U x y z t and distributed output quantities,   , , , Y x y z t , at initial and boundary conditions given. Distributed parameter systems are very frequently found in various technical and non-technical branches with limited number of manipulated input quantities, or actuators. These lumped input quantities by means of interaction of fields and quantities generate distributed output of real DPS. Representation of such DPS is either in the form of LDS, Fig. 1 a), or in the form of LDS with zero-order hold unit H (HLDS), when discrete-time lumped input quantities are used, Fig. 1. b) , (Hulkó et al., 1981(Hulkó et al., , 1987(Hulkó et al., , 1998.

Dynamics of LDS
Distributed output of the linear LDS from zero initial conditions either in continuous, or in discrete-time (DT) is in the form: Dynamics of LDS/HLDS is decomposed to the time and space components:

Feedback control loop based on HLDS dynamics
Decomposition of dynamics enables also to decompose the control synthesis (CS) to time synthesis (TS) and space synthesis (SS) tasks in the feedback control loop of the distributed parameter system, Fig. 2.
Let us consider a step change of distributed parameter control quantity The goal of the control synthesis is to generate a sequence of control inputs   U k in such manner, that in the steady-state, for k   , the control error   , E x k will approach its minimal value   , E x   in the quadratic norm: First, in the SS blocks, the approximation both of sampled distributed controlled quantity

 
, W x k is assumed, the space control synthesis is performed in each time step k, which gives parameters. Graphical interpretation of the approximation problem (9) for HLDS defined on 1D space x is on Fig. 3.  in single-parameter control loops. During the control process, for k   the control task (7) is accomplished.
Finally, we may state as a summary, that in the feedback control of DPS with dynamics represented in the form of HLDS, the control synthesis is performed as: Space Tasks of Control Synthesis -as approximation tasks. Time Tasks of Control Synthesis -on the level of lumped parameter control loops.

Robust control system
In general, a mathematical model for the plant dynamics is the basis for analysis and design of control systems. Also for LDS representation of DPS lumped and distributed models are used. However, in practice, no mathematical model describes exactly a physical process. It is obvious, that although no model represents the process exactly, some of them will do so with greater accuracy than others.
The theory of the robust control represents one of the possible approaches to the control system design in the presence of uncertainty. The goal of the robust system design is to retain a good quality of system performance in spite of model inaccuracies and changes. For the design techniques, the following requirements are supposed to be fulfilled: formulation of nominal plant model, different plant uncertainty models and requirements for both, robust stability and performance.

Sources of uncertainties in the LDS structure
LDS representation of DPS means decomposition of dynamics to space and time components. Uncertainties may occur in both, time and space components.
In distributed parameter control system, according to Fig. 2, single-input, single-output control loops in the block TS are tuned as closed feedback control loops using usual methods. In these loops, as models of the controlled system, transfer functions in the z-domain are used. These transfer functions describe the dynamics between sequences In this case, the sources of uncertainties are given by:  procedure of dynamics modelling and possible change of parameters in models (4), (5)  solution of approximation problem (8), (9), where lumped quantities are obtained In order to treat uncertainties, it will be further assumed that the dynamic behaviour of a plant is described not by a single linear time invariant model, but by a family of linear time invariant models, i  . This family in the frequency domain, e.g. for models is the nominal plant model. Any member of the family i  fulfils the conditions: where   ai L j is an additive uncertainty and   ai L  is the bound of additive uncertainty. If we wish to work with multiplicative uncertainties, we define the relations: where   mi L j is a multiplicative uncertainty and   mi L  is the bound of multiplicative uncertainty.

Design of IMC robust controllers
A robust control system for HLDS will be designed using the Internal Model Control (IMC) strategy (Morari & Zafiriou, 1989) with the general structure depicted in Fig. 4. a). It is possible to transform this structure to the classical feedback control loop, Fig. 4. b) and to incorporate it into the TS block of the DPS feedback control system. The relationship between the classical feedback controller C and the IMC controller Q for the nominal model P  of the controlled process P is as follows, and vice-versa: It is well known that IMC strategy has the following properties: 1. Dual stability: Assume a perfect model,   P P   and if the controller and process are stable, then the IMC structure guarantees the closed-loop stability.
2. Perfect control: Assume a perfect model,   P P   and the closed-loop system is stable, while 1 Q P    , then there is no output steady-state error for set-point variance and disturbances.
The IMC structure thus provides the following benefits with respect to classical feedback: better dynamic response, system stability and robustness. One can search for Q instead of C without any loss of generality.
Structure of the distributed parameter feedback robust control system DPS with HLDS dynamics and IMC controllers is on Fig. 5.
are obtained from solution of the following minimization problem: Resulting IMC controller with filter   Fi Q z is in following form: Parameter of the filter i  is the only one tuning parameter to be selected by the user to achieve the appropriate compromise between performance and robustness and to keep the action of the manipulated variable within bounds. It must be chosen with respect to both, robust stability and robust performance condition: and low-pass filter (19)

Benchmark casting plant
The casting mould is one of the key components of a casting. It is well known, that the quality of the castings is affected strongly by the surface quality and the distribution of temperature in the mould, which has both time and space dependence. For study of the physical phenomena occurring during the casting solidification, from a DPS control point of view, control system development as well as mathematical model validation, a benchmark of the casting processes was designed (Belavý et al., 2009). At the study of casting processes and design of experimental plant, simulation studies in virtual software environments ProCAST and COMSOL Multiphysics were used.

Construction of the benchmark casting plant
Scheme of the benchmark casting plant is depicted in Fig. 6. The core item is the two-part steel mould of a complex-shape mounted in the frame of the ejector mechanism, Fig. 7. This mechanism is hinge-mounted to the main frame, to enable tilting of the mould for optimal filling and metal flow. Further, a hydraulic cooling circuit, which consists of an array of induction motor driven roller vane pumps, a bunch of hoses, a flow divider, a collector with built-in check valves, a plate heat exchanger and an expansion tank. The main cooling circuit is divided into five independent circuits thus enabling the control of heat extraction from the casting via the chills. The coolant flow is controlled by means of frequency converters, since volumetric pumps are used. The main heating circuit is also divided into five independent circuits.  Inside of the casting mould are built-in 26 electric heating elements, each with maximal heating power 400 W. Heating elements are grouped to 5 zones and their heating power is actuated by the input voltage range of (0 -10) V. In the body of the mould is also placed 7 water-cooled copper chills and 11 thermocouples, Fig. 8., Fig 9. Coordinates of measuring points of thermocouples in x-y plane are given in Tab. 1 and z -coordinate is -0.05 m. Location of built-in elements has been carefully designed based on simulation studies in software environments ProCAST and COMSOL Multiphysics, in order to have the possibility of preheating the mould in 5 zones achieving desired temperature profile as well as directional solidification of the casting by means of active heat removal. The temperature field in the mould-casting system is possible to estimate through interpolation of data, measured by thermocouples.

Measurement and control scheme in MATLAB & Simulink
The measurement and control task of temperature field in the permanent casting mould was performed in the MATLAB & Simulink environment, where a mould_exp_robust.mdl scheme was setup, Fig. 10. The scheme is composed of three main subsystems, namely: SENSING, HEATING and COOLING. Utilizing these subsystems and communication interface between process and computer, it is possible to measure dynamical characteristics of temperature field at zone heating. It is also possible to execute the experiment of controlled preheating of permanent mould before casting operation and controlled cooling during casting solidification.

Communication interface is performed by means of Simulink Real-Time Windows Target Toolbox and the communication is triggered by Connect to Target icon.
The SENSING subsystem enables the temperature measurement by thermocouples No. 1 to 11, which are permanently located in the bottom part of the steel mould, Fig. 9. These are marked on the scheme as P1 to P11. The subsystem also enables temperature measurements during casting solidification by temporal thermocouples W1 to W9, located in the casting domain. Thermocouple sensing junctions are located in the middle of the arm cross section, right above permanent thermocouples "P", except the node of the casting, where one temporal thermocouple points to the center of the node. Sensors S1 to S7 measure the temperature of cooling water in embedded chills. Display blocks of the sensors "P" and "W" in the scheme correspond to the real positions of sensors. Measured temperatures are saved to the data file and continuously displayed on the Scope block.
The HEATING subsystem has two basic operational regimes, which are activated by Switch2 block. In the manual control regime, it is possible to set heating performance in range (0÷10) Volt and measure the temperature transient characteristics with P1 to P11 thermocouples in given locations of the mould. The second regime activated by Switch2 block enables to perform controlled preheating of the casting mould to desired temperature profile W defined in 11 points, where thermocouples P1 to P11 are located. The control task is performed by DPS Robust IMC Control Synthesis block.  Fig. 11. , (Hulkó et al., 2003(Hulkó et al., -2010.

Distributed Parameter Systems Blockset for MATLAB & Simulink
The library of DPS Blockset shows Fig. 12. Blocks HLDS and RHLDS serve for modelling of distributed parameter systems as lumped-input/distributed-output systems with zero-order hold units. The block DPS Control Synthesis provides feedback to distributed parameter controlled systems in control loops with blocks for discrete-time PID, Algebraic, State-Space and Robust Control. The block DPS Input generates distributed quantities which can be used as distributed control quantities or distributed disturbances, etc. DPS Display presents distributed quantities with many options including export to AVI files. The block DPS Space Synthesis performs space synthesis as an approximation problem.
The block DPS Wizard in step-by-step operation, by means of several model examples with default parameters on 1D-3D definition domains, gives an automatized guide for arrangement and setting distributed parameter control loops. The block Demos contains examples oriented to methodology of modelling and control synthesis. The block Show contains motivation examples such as: Control of temperature field of 3D metal body (the controlled system was modelled in the virtual software environment COMSOL Multiphysics); Control of 3D beam of "smart" structure (the controlled system was modelled in the virtual software environment ANSYS); Adaptive control of glass furnace (the controlled system was modelled by Partial Differential Equations Toolbox of the MATLAB ), and Groundwater remediation control (the controlled system was modelled in the virtual software environment MODFLOW). The block Tutorial presents methodological framework both for formulation and solution of control tasks for distributed parameter systems.

Dynamics of the temperature field in the casting mould
For control synthesis purpose, LDS/HLDS models of temperature fields in the casting mould have been created by means of evaluation of real-time experiments. The measurement of temperatures fields in the casting mould was performed by MATLAB & Simulink scheme mould_exp_robust.mdl, Fig. 10.

Experimental identification of transfer functions
In the casting mould lumped inputs   1,5 i i U  are heating elements which act on sub-domains , (Zone 1 -5). Distributed output is the temperature field of the casting mould.
Temperatures in the casting mould were measured by 11 thermocouples as a time-response to the step change of heating power, which was activated by the input voltage step from 0 to 2,5 V, for heating elements separately in each zone. Results of measurements for Zone #1, #3 and #5 are depicted in Fig. 13, where are time and x-y space dependences of temperatures and Fig. 14 presents temperature profiles in steady-states.    (24) for identification has been chosen, see Fig. 15., where are also results of identification from actuating in zone #1. Comparison of measured and identified model output in zone #1 is presented in Fig. 16. Identified parameters of transfer functions Table 3. Continuous transfer functions with structure (24) are converted by means of function zpk to zero-pole-gain format (ZPK). Then, for control synthesis purposes, they are transformed to discrete transfer functions

Uncertainty analysis in the space and time domain
The proposed structure of the DPS control systems, Fig. 2, Fig. 5., are significant by decomposition of the control synthesis into the space and time subtasks. In the SS blocks, the approximation both of distributed controlled quantity  

Robust control process
Robust control of temperature fields of the casting mould as distributed parameter system was performed with MATLAB & Simulink and DPS Blockset software support. Robust controllers designed by the IMC control strategy were first optimised through their tuning parameters   i i  and then implemented to the control scheme for the real-time control of temperature fields of the casting mould in the benchmark casting plant.

Optimization of tuning parameters
In the MATLAB & Simulink environment, by means of the DPS Blockset, distributed parameter system of robust control, mould_robust_DZPK.mdl was arranged, see Fig. 21. It is DPS feedback control loop, where the DPS Robust Control Synthesis block includes both, time and spatial part of the control synthesis, see Fig. 22 a). In this case, the control system consists of five single parameter control loops, each for one zone of the mould, where discrete robust controllers based on IMC structure are used.
Parameters of filters were optimized in the presence of constraints of the criterion function in order to assure nearly aperiodic course of the quadratic norm of the distributed control error with respect to the robust stability and robust performance conditions, see Fig. 23 and optimization progress presents Fig. 24   Sensitivity with respect to uncertain element ... 'Tp1' is 100%. Increasing 'Tp1' by 25% leads to a 25% decrease in the margin. 'Tp2' is 54%. Increasing 'Tp2' by 25% leads to a 14% decrease in the margin.

Real-time robust control of preheating in the casting mould
Real-time robust control of the temperature fields of the casting mould in the benchmark casting plant was performed by means of Simulink designed block scheme mould_exp_robust.mdl. The structure of the block DPS Robust IMC Control Synthesis is depicted in Fig. 25, where are blocks Space Synthesis Y and Space Synthesis W, and block Robust controllers based on IMC with the same structure as in Fig. 22 b). Control process was performed for the reference variable -temperature profile given in 11 positions, where thermocouples are embedded. Results of the real-time robust control process are on Fig. 26. The quality of control both in the time and space domain is given by the quadratic norm of the distributed control error.

Conclusion
The aim of this chapter was to present the engineering approach for the robust control of DPS, which opens a wide space for novel applications of the toolboxes and blocksets of the MATLAB & Simulink software environment. This approach is based on the general decomposition of controlled DPS dynamics, represented by transient and impulse characteristics, into time and space components. Starting out from this dynamics decomposition a methodical framework was presented for the decomposition of control synthesis into the space and time subtasks. In the space domain an approximation problems were solved, while in the time domain the control synthesis was performed by lumped parameter SISO control loops, where various wellknown methods for design of controller is possible to utilize. The advantage of this approach is the relatively simple LDS model of DPS, which is directly suitable for control purposes and can be easily identified from input-output data by means of classical techniques.
Currently, it is interesting to formulate and solve tasks of control in various engineering branches, including the casting technology, by means of methods and tools of distributed parameter systems. Methodical approach presented in this chapter demonstrates simple possibilities, how to exploit the distributed dynamical characteristics on complex definition domains, obtained by evaluation of measured data for robust IMC control synthesis of DPS with respect to uncertainty of models and the real-time control according to technological requirements. Cyril Belavý, Gabriel Hulkó and