UDK 621.74.047:66.011 ISSN 1580-2949 Original scientific article/Izvirni znanstveni članek MTAEC9, 46(4)325(2012) A FUZZY-BASED OPTIMAL CONTROL ALGORITHM FOR A CONTINUOUS CASTING PROCESS ALGORITEM, KI TEMELJI NA MEHKI LOGIKI, ZA OPTIMALNI NADZOR KONTINUIRNEGA PROCESA ULIVANJA Tomas Mauder, Cenek Sandera, Josef Stetina Brno University of Technology, Faculty Mechanical Engineering, Technicka 2, Brno, Czech Republic ymaude00@stud.fme.vutbr.cz Prejem rokopisa - received: 2011-10-19; sprejem za objavo - accepted for publication: 2012-02-24 Nowadays, continuous casting is used for processing almost one hundred percent of liquid steel into an intermediate shape. Steel with a poor structure and many defects is not acceptable to final customers, and therefore the producers must put strong emphases on its quality. Many serious steel defects can be eliminated by a preceding computer simulation, optimization and a proper control of the casting process. This paper describes an algorithm that optimizes the control parameters to ensure both a high production rate and a high quality of the products. These controlled parameters are the casting speed and all the cooling rates in the secondary cooling zone. The main principle of the algorithm is to get the surface and core temperatures to the prescribed values corresponding to the required ductility of steel. The whole optimization procedure consists of two separate parts, i.e., the numerical simulation of the process and the decision-making part based on fuzzy logic for modifying the control parameters. By incorporating the fuzzy logic into the optimization process, the algorithm has a very robust behaviour and an easy adaptation for different grades of steel. This algorithm runs in an off-line version and can be used as a preparation tool for a real casting process where the proper setting of the casting parameters is crucial for achieving high-quality products at an economical price. Keywords: fuzzy optimization, heuristic optimization, temperature field, continuous casting Postopek kontinuirnega ulivanja za pridobivanje vmesne oblike se danes uporablja pri obdelavi skoraj stoodstotnega deleža tekočega jekla. Jeklo s slabo strukturo in številnimi napakami ni sprejemljivo za končne kupce, zato morajo proizvajalci veliko pozornosti posvetiti kakovosti jekla. S predhodnimi računalniškimi simulacijami, z optimizacijo in ustreznim nadzorom procesa ulivanja se lahko odpravijo mnoge resne napake jekla. V tem prispevku je opisan algoritem za optimizacijo parametrov nadzora, ki zagotavlja visoko stopnjo proizvodnje ter visoko kakovost izdelkov. Med omenjene parametre nadzora spadata hitrost ulivanja in hitrost ohlajevanja v sekundarni coni ohlajanja. Glavno načelo algoritma je, da površina in temperatura sredice dosežeta določene vrednosti, ki ustrezajo zahtevani razteznosti jekla. Celoten postopek optimizacije obsega dva med seboj ločena dela - numerično simulacijo postopka in pa del, namenjen odločanju, ki temelji na mehki logiki za spreminjanje parametrov nadzora. Z integracijo mehke logike v postopek optimizacije algoritem doseže zelo robustno vedenje in enostavno prilagoditev različnim vrstam jekla. Ta algoritem deluje v različici »off-line« in se lahko uporablja kot pripravljalno orodje za dejanski proces ulivanja, kjer je ustrezna nastavitev parametrov ulivanja ključnega pomena za proizvodnjo visokokakovostnih izdelkov po ekonomični ceni. Ključne besede: mehka optimizacija, hevristična optimizacija, temperaturno polje, kontinuirno ulivanje 1 INTRODUCTION ferrostatic forces, a poor adjustment of the casting speed, roll gaps, bending and straightening of strand, etc. The Nowadays, continuous casting represents the main cracks caused by bending and straightening are method worldwide for forming molten steel into semi- influenced by the ductility of steel, and it is well finished products, such as slabs, blooms, and billets. The documented that the ductility of steel is reduced over quality of final products plays one of the most significant specific temperature ranges depending on the chemical roles for customers, and thus to satisfy their demands composition1. and maintain the highest possible productivity is a key Conducting industrial trials is very expensive, goal for each steelmaker. time-consuming and in some cases even impossible, Quality control is a prerequisite for continuous which makes a computer simulation the only sustainable casting and cannot be achieved without a proper option. There exist various ways of computer or matheknowledge of the main physical influences of the casting matical optimization and each of them has its advantages process, i.e., the solidification, the micro and macro and drawbacks. Therefore, to ensure the correctness of segregation, the crack formation, etc. During the the algorithm, it is necessary to validate the simulation processing of the steel in the continuous caster, there can output using real, measured data. arise many problems with its final quality. The solid Previous researches concerned with the optimal shell is permanently subjected to thermal and mechanical control of a continuous casting process were generally stresses and this can result in cracks or the breakout of based on simplified 1D or 2D temperature-field models1 liquid steel through the solid shell. The sources of the and were optimized by mathematical programming or mechanical stresses are the friction in the mould, heuristic methods, e.g., genetic algorithms2, the firefly algorithm3 or by neural networks. Many of these models are based on simplified assumptions, and therefore they describe the casting process very roughly and not satisfactory. We developed our original numerical model for the temperature field of the real caster and its results were validated by temperature measurements acquired by pyrometers. This more precise model simulates the process controlled by several numerical parameters and the goal is to find their values such that the resultant temperature field is optimal. This algorithm is inspired by our previous research3 and it is enhanced by a fuzzy logic inference mechanism, which makes its behaviour more robust and its setting easier to adjust. 2 PROBLEM DESCRIPTION AND AN OUTLINE OF THE ALGORITHM The natural inclination of all steelmakers is to cast as fast as possible but with a preservation of the required material quality. Thus, the problem of optimal casting control can be formulated as follows. The properties of the final material are highly dependent on its temperature courses emerging during the casting process. Therefore, we need to adjust the parameters of a caster in such a way that the final temperature field is optimal. We can control the process by the regulation of the casting speed and by all the cooling rates in the secondary cooling zone. Before the optimization starts, it is necessary to define the temperature field, which leads to the optimal material and mechanical properties. These optimal temperature courses are defined by experts (casting operators, material scientists, etc.) and they differ with the chemical composition of the steel. To deal with the high variability of the described problem, the heat and mass phenomenon must also be taken into account, and therefore the employed numerical model has to cover all the phase and structural changes. The search for the optimal cooling rates and the highest possible casting speed is left to the following fuzzy-logic heuristic algorithm. To describe the optimal temperature field, the experts define a set of points along the caster and prescribe their optimal temperature ranges. Then the algorithm randomly generates the values for all the input parameters and simulates the final temperature field. After comparing the computed temperatures with the prescribed ones, the algorithm decides which nozzles need to be adjusted and how to change the casting speed. Particular decisions for changing the values are taken by the fuzzy inference logic. These steps are iteratively repeated until the required temperatures are successfully reached. 3 MATHEMATICAL MODEL OF THE TEMPERATURE FIELD The developed numerical model of the 2D temperature field computes the formation of the solidification and the temperature distribution of the cast strand. This temperature field is modelled by the Fourier-Kirchhoff equation1-3, where the velocity component Vy(m/s) is considered only in the direction of casting: dH dr = k eff( T) A (T) + V, dH dy (1) Equation (1) describes a transient heat transfer, where keff/(W/m K) is the effective thermal conductivity, T/K is the temperature, H/(J/m3) is the volume enthalpy, r/s is the time and x and y are spatial coordinates. The phase and structural changes are included in the model by the latent-heat accumulation method, where the enthalpy is used as the primary variable and the temperature is calculated from a defined enthalpy-temperature relation- ship3: H-j p( ^)c(^) - p(^)AH dT (2) where AH/(J/kg) is the latent heat, p./(kg/m3) is the density, c/(J/kg K) is the specific heat capacity and fs is the solid fraction. In order to have a well-posed problem, the initial and boundary conditions must be stated. The boundary conditions include the heat fluxes in the mould and under the rollers, the forced convection under the nozzles and the free convection and radiation in the tertiary cooling zone. Equation (1) is discretized by the finite-difference method3,4 using an explicit formula for the time derivative. The mesh for the finite-difference scheme is non-equidistant and its nodes are adapted to the real rollers and nozzle positions. Our numerical model allows us to apply various enthalpy-temperature functions and thermal conductivity-temperature curves, and thus the temperature field can be calculated for various steels just by redefining the chemical composition. Thermo-physical parameters such as the thermal conductivity, the density, the specific heat capacity and the enthalpy and their temperature dependence are computed from a specific chemical composition of steel by using the solidification analysis package IDS. The numerical model is designed and verified for the radial slab caster operating in EVRAZ VITKOVICE STEEL, a. s. The caster contains 12 coolant circuits and the cross-section of the investigated slab is 1550 mm x 250 mm. 4 DESCRIPTION OF THE ALGORITHM The algorithm searches for the optimal values of the control parameters (cooling rates and the casting speed). The inputs to the algorithm are the required temperature ranges in the points along the caster and the maximum metallurgical length. For each cooling circuit, the maximum cooling rate is defined. The initial value of each control parameter is randomly chosen from the permitted range and the model is consequently evaluated by the numerical simulation. From the computed result, the algorithm extracts temperatures in the controlled points and by a comparison with the prescribed values determines their deviations (errors). Also, the metallurgical length (the length from the meniscus to the point where the steel is fully solidified) is computed. Using all this information the algorithm infers the modifications for all the controlled input parameters. The temperature courses at the controlled points are dominantly influenced by the two preceding coolant circuits. The closer circuit has a greater influence. Therefore, each coolant circuit defines a numerical value for each controlled point that describes how much it impacts on the temperature in the controlled point. These values are expertly estimated (in the range from 0 to 10) and they are strongly related to the distance from the controlled points. In the presented model there are two fuzzy inference mechanisms - one for the modification of the cooling rates and one for the control of the casting speed. The inference rules for the cooling rates have the following form: "IF error IS adj1 AND impact IS adj2 THEN modification IS adjS''. For each controlled point the adj1 describes the temperature error, the adj2 is the cooling impact of the given circuit and the ad;j3 is the modification inferred for the cooling rate. Table 1 shows all the rules in the matrix form. Table 1: The dependences of adj3 on adj1 and adj2. The abbreviations stand for Very Small, Small, Medium, Big and their fuzzy sets are equidistantly distributed along the corresponding universes. Tabela 1: Odvisnost označevalca 3 od označevalca 1 in označevalca 2. Okrajšave pomenijo VS - zelo majhen, S - majhen, M - srednji, B -velik, njihove mehke nastavitve pa so v enakih razdaljah razporejene po ustreznih splošnih parametrih. adi2 / adi1 VS S M B S VS VS S S M VS VS S M B VS S M B Sometimes, one circuit can get several different modifications and if it happens the algorithm takes the one with the highest absolute value. The defuzzification method is the standard centre-of-gravity function. If the maximum absolute value of all the temperature errors does not exceed a particular limit, the algorithm consequently computes the modification for the casting speed. The reason for introducing the limit is that if the maximum error is too big (i.e., the solution is far from the optimum), we have no information about whether the speed is high or not and at first it is better to stabilize the process and then to infer the speed modification. The rules for the modification of the casting speed (Table 2) are in the form: "IF maximum_error IS adj4 AND metalurgical_length IS iadjS THEN modification IS adj6'. The adj4 describes the maximum error over all the controlled points, the adj5 is the metallurgical length and the adj^> is the inferred modification of the casting speed. For a detailed explanation of the fuzzy logic and the fuzzy inference, see5. The values of the impacts for each cooling circuit to the controlled point are chosen to be 8 for the closest preceding circuit and 2 for the circuit placed before (except the last controlled point, where the distance to the second closest circuit is much longer, and therefore the values are 9 and 1). Table 2: The dependences of adj6 on adj4 and adj5. The abbreviations stand for Very Small, Small, OK, Big, Very Big and More, Little More, Little Less, Less, Nothing and their fuzzy sets are equidistantly distributed along the corresponding universes. Tabela 2: Odvisnost označevalca 6 od označevalca 4 in označevalca 5. Okrajšave pomenijo VS - zelo majhen, S - majhen, OK, B - velik, VB - zelo velik in MO - več, LM - malo več, LL - malo manj, L -manj, N - nič, njihove mehke nastavitve pa so v enakih razdaljah razporejene po ustreznih splošnih parametrih. adi4 / adi5 VS S OK B VB S MO LM LM LL L M LM N N N LL 5 RESULTS OF EXPERIMENT The examined grades of steel are the low-carbon steel S355J0H and the high-carbon steel X210CrW12. The expertly defined temperature ranges in the controlled points are depicted in Figure 1 and Figure 2 (the grey rectangles). The crucial reason for defining these temperature ranges in this way is that approximately the first half of caster is curved and therefore to decrease the mechanical stresses it is better to keep the temperature 1600 1400 1200 S! ■i 1000 casting speed 0.99 m/min ntc • |3«. m «01. 4M. 321.109, 270. 204. 293, 239. 1392, 1333] 800 600 400 \ \ , \ s N > h 1 v — top side temperature - - core temperature — bottom side temperature 10 15 20 distance below meniscus [m] 25 Figure 1: Final temperature distributions (S355J0H) Slika 1: Končne razporeditve temperature (S355J0H) 1600 1400 1200 casting speed 0.44 m/min ntc • (532. 539.4«9. 339, 300, 20«. 172. 175. 143. 149, 880, 897] U — 1000 800 600 400 200 K V V V ^ " " ^ ^ > H [\ — top side temperature - - core temperature 'irNiUf^ki i' — bottom side temperature ' "»Tf »'I 10 15 distance below meniscus [m] 20 Figure 2: Final Temperature distributions (X210CrW12) Slika 2: Kon~ne razporeditve temperature (X210CrW12) above a certain level (for S355J0H, 1000 °C, and for X210CrW12, 800 °C). The value of the maximum metallurgical length is 20 m (from the practice) and the casting temperature is 1550 °C. Reaching the temperature courses on surface ensures that the mechanical and material qualities of the final material will satisfy the demands. The optimal values computed by this algorithm (casting speed and heat-transfer coefficients) are shown in Figure 1 and Figure 2. Usually, the most important indicator characterizing the efficiency of iterative optimization algorithms is the number of evaluations of the model. The computation of the numerical model is very time-consuming, and therefore each repetition can significantly prolong the computations. Our algorithm is able to find the optimal input parameters in 50 evaluations, on average. The tests ran several times for different grades of steel and the number of evaluations never exceeded 65. The computation time takes a few hours. 6 CONCLUSION The problem of the optimization of a continuous casting process can be efficiently solved by the described algorithm. The algorithm based on heuristics and incorporating the fuzzy logic behaves in a very robust way and it is easily adaptable to any grade of steel and caster geometry. The number of evaluations of the included numerical model is low, and thereby the algorithm proves its high efficiency. Further research will be focused on extending the numerical model to 3D and a generalization of the used optimization constraints. Acknowledgement The authors gratefully acknowledge the financial support from the projects GA106/08/0606, GA106/09/ 0940 and GACR P107/11/1566, funded by the Czech Science Foundation, project ED0002/01/01 - NETME Centre and Specific research BD 13101001. 7 REFERENCES 1 J. P. Birat, et al.: The Making, Shaping and Treating of Steel: Casting Volume: 11th Edition. Alan W. Cramb, Pittsburgh, PA, USA: The AISE Steel Foundation, 2003, 1000 s 2 C. A. Santos, J. A. Spim, A. Garcia, Mathematical modeling and optimization strategies (genetic algorithm and knowledge base) applied to the continuous casting of steel, Engineering Applications of Artificial Intelligence, 16 (2003), 511-527 3 D. M. Stefanescu, Science and Engineering of Casting Solidification, Second Edition, New York, Springer Science, 2009, 402 4T. Mauder, C. Sandera, J. Stetina, M. Seda, Optimization of Quality of Continuously Cast Steel Slabs by Using Firefly Algorithm, Mater. Tehnol., 45 (2011) 4, 347-350 5 H. T. Nguyen, E. A.Walker, A First Course in Fuzzy Logic, CRC Press, 1999, 392