Method for a pouring control and a storage medium for storing programs for causing a computer to carry out a process for controlling pouring

Abstract

To enable a ladle-tilting automatic pouring device to take less time for identification of the parameters and the device to pour highly precisely by sequentially updating pouring model parameters according to the pouring situation, the present pouring control method is based on a mathematical model of a process from input of control parameters to pouring of molten metal, the method including: identifying, using an optimization technique, a flow rate coefficient, a liquid density, and a pouring start angle that is a tilting angle of the pouring ladle when the flowing of the molten metal starts, which are the control parameters in the mathematical model, based on weight of liquid that flows out of the pouring ladle and tilting angle of the ladle that are measured during pouring, and a command signal that controls the tilting of the pouring ladle; and updating the control parameters to the identified control parameters.

Claims

1. A pouring control method for controlling pouring based on a mathematical model of a pouring process from input of control parameters to pouring of molten metal using a pouring ladle in an automatic pouring device with a tilting pouring ladle that pours the molten metal into a mold by tilting the pouring ladle that holds the molten metal, comprising: identifying, using an optimization technique, a flow rate coefficient, a liquid density, and a pouring start angle that is a tilting angle of the pouring ladle at which a flow out of the molten metal starts, wherein the flow rate coefficient, the liquid density, and the pouring start angle are the control parameters in the mathematical model, based on weight of liquid that flows out of the pouring ladle and tilting angle of the ladle that are measured during pouring, and a command signal that controls the tilting of the pouring ladle, and updating the control parameters to the identified control parameters, wherein the flow rate coefficient the liquid density, and the pouring start angle are identified by optimizing an evaluation function that is represented by a following equation,
{c.sub.id,.sub.sid,.sub.id}=arg min{.sub.0.sup.T(W.sub.Lex(t)W.sub.Lsim(t,c.sub.sim,.sub.ssim,.sub.sim)).sup.2dt+w.sub.1(c.sub.avgc.sub.sim).sup.2+w.sub.2(.sub.avg.sub.sim).sup.2}, where c.sub.id is an identified flow rate coefficient, .sub.sid is an identified pouring start angle, .sub.id is an identified liquid density, T is operating time required to pour molten metal into one mold, W.sub.Lex is data on outflow weight from the pouring ladle obtained from the automatic pouring device with a tilting-type ladle, W.sub.Lsim is outflow weight obtained by the simulation with the mathematical model using the lade tilting angle, c.sub.sim is a flow rate coefficient that was used in the simulation, .sub.ssim is a pouring start angle that was used in the simulation, .sub.sim is a liquid density that was used in the simulation, C.sub.avg is an average value of flow rate coefficients used until previous time, .sub.avg is an average value of liquid densities used until previous time, w.sub.1 is a weight coefficient used to control the variation of the flow rate coefficient for every pouring, and w.sub.2 is a weight coefficient used to control the variation of the liquid density for every pouring.

2. The pouring control method according to claim 1, wherein the flow rate coefficient and the liquid density are identified and updated every time one pouring is completed, and wherein an approximate function between the identified pouring start angle and a corresponding weight of liquid within the pouring ladle is calculated and updated after the consecutive pouring processes by the pouring ladle are completed.

3. The pouring control method according to claim 1, wherein the optimization technique is a Down-hill simplex method.

4. The pouring control method according to claim 2, wherein the optimization technique is a Down-hill simplex method.

Description

BRIEF EXPLANATION OF FIGURES

(1) FIG. 1 is a schematic perspective view that shows one example of the automatic pouring device with the tilting-type ladle.

(2) FIG. 2 is a block diagram of the pouring control method.

(3) FIG. 3 is a flowchart that shows the pouring control method for identifying and updating parameters.

(4) FIG. 4 is a schematic cross-sectional view of the pouring ladle.

(5) FIG. 5 is a schematic perspective view that shows the tip end of the lip of the pouring ladle.

(6) FIG. 6 is a schematic diagram that shows the result of a pouring experiment.

(7) FIG. 7 is a schematic diagram that shows the result of a pouring experiment.

(8) FIG. 8 is a schematic diagram that compares the result obtained from the shape of the pouring ladle with the approximate function with regard to the relationship between the pouring start angle and the weight of the liquid within the pouring ladle before pouring.

DESCRIPTION OF EMBODIMENTS

(9) The pouring control method of the present invention is explained by reference to the Figures.

(10) An example of the automatic pouring device with a tilting-type ladle that employs the pouring control method of the present invention is shown in FIG. 1. The automatic pouring device with a tilting-type ladle 1 (hereafter, automatic pouring device 1) comprises a pouring ladle 10 and servomotors 11, 12, and 13. The pouring ladle 10 carries molten metal. One of the servomotors is a servomotor 11 that tilts and also turns the ladle 10 around an axis . Another servomotor 12 moves the ladle 10 back and forth. The third servomotor 13 moves the ladle 10 up and down.

(11) Since the servomotors 11, 12, and 13 each have rotary encoders, the position and the tilting angle of the pouring ladle 10 can be determined. The servomotors 11, 12, and 13 are configured to be given a command signal from a computer. The computer in this specification denotes a motion controller such as a personal computer, a micro computer, a programmable logic controller (PLC), or a digital signal processor (DSP).

(12) A load cell is arranged at the lower end of a rigid structure that includes the pouring ladle 10 or at the lower end of the automatic pouring device 1, to measure the weight of the pouring ladle 10 that includes the liquid.

(13) By using the above-mentioned configuration, the automatic pouring device 1 can discharge the molten metal from the lip of the pouring ladle 10a, and pour the molten metal inside a mold 20 through a sprue of the mold 20a by controlling the servomotors 11, 12, and 13 to convey the pouring ladle 10 along a predetermined track.

(14) A mathematical model of the pouring process in the automatic pouring device 1 based on fluid mechanics will be derived here to build a model-based pouring control system that is a pouring control system based on the mathematical model. FIG. 2 shows a configuration example of the model-based pouring control system. Here, FIG. 2 shows a two-degree-of-freedom type pouring control system in which a feedforward control and a feedback control are incorporated.

(15) Once the computer 14 is given the desirable target outflow weight and the target pouring flow rate pattern, the computer 14 adjusts and outputs the command signal to the automatic pouring device 1 to achieve the target pouring flow rate and the target outflow weight, where the command signal may become the speed command and/or the position commands, depending on the control mode of the servomotors 11, 12, and 13. In addition, various aspects, such as voltage and pulses, can be adopted as the command signal.

(16) When pouring, the tilting angle of the ladle is measured by the rotary encoder, and the weight of the liquid within the pouring ladle is measured by a load cell provided on the automatic pouring device 1. The outflow weight of the liquid that outflows from the pouring ladle 10 can be measured by calculating the difference between the weight of the liquid within the pouring ladle before pouring and the weight of the liquid within the pouring ladle during the pouring.

(17) The measured tilting angle of the ladle and weight of liquid within the pouring ladle are output to the computer 14. The computer 14 controls the pouring operation based on them. Incidentally, the pouring control system in FIG. 2 may become a feedforward-type pouring control system by removing the feedback loop in FIG. 2.

(18) The computer 14 identifies and updates the model parameters based on the command signal, and the acquired tilting angle of the ladle and weight of the liquid within the pouring ladle. The pouring control system generates the command signals for the servomotors 11, 12, and 13, depending on the model parameters, by acquiring the command signal, the weight of liquid within the pouring ladle, and the tilting angle of the ladle that are detected through one pouring operation, by using these data and the mathematical model of the pouring process to identify the flow rate coefficient, the liquid density, and the pouring start angle, which are the model parameters of the pouring process, and by updating the model parameters within the pouring control.

(19) Next, based on the flowchart of FIG. 3, the identification and update process of the model parameters will be explained. At step 1, the initial model parameters and functional relationship (calibration curve) between the pouring start angle and the weight of the liquid within the pouring ladle are given the pouring control as parameters set for the pouring control. The pouring start angle is the tilting angle of the pouring ladle 10, at which a flow out of the molten metal begins. The initial model data as the initial model parameters include the shape of the pouring ladle, the liquid density, and the flow rate coefficient. The values that are employed for the pouring ladle design are used as data on the shape of the pouring ladle. The values that are considered to be appropriate through experiments and/or experience are used for the liquid density and the flow rate coefficient. The functional relationship between the pouring start angle and the weight of liquid within the pouring ladle can be obtained by calculating the volume of the liquid with which the pouring ladle is filled, which corresponds to the tilting angle of the ladle from data on the shape of the pouring ladle, multiplying the volume by the liquid density, and formulating the function. Incidentally, it is assumed that the pouring ladle 10 at this stage is already supplied with the molten metal and is ready to carry out the pouring operation.

(20) At step 2, the pouring machine is controlled based on the mathematical model discussed below, and the pouring from the pouring ladle 10 into the mold 20 is carried out.

(21) At step 3, the liquid density and the flow rate coefficient are identified as parameters to be updated by using an optimization technique explained later based on the outflow weight from the pouring ladle 10, the tilting angle of the ladle, and the command signal data that are acquired during a pouring operation from the pouring ladle 10 into a mold 20.

(22) At step 4, the identified pouring start angle and the weight of the liquid within the pouring ladle that were measured before pouring are stored as a set of data in the computer 14.

(23) At step 5, the liquid density and the flow rate coefficients that were input as the initial parameters to the pouring control and were used for the pouring control are updated online so that they are replaced by the liquid density and the flow rate coefficient, respectively, that were identified at step 3.

(24) At step 6, the computer 14 determines whether the pouring ladle 10 was supplied with the molten metal after or at step 2. If the pouring ladle 10 was not supplied with the molten metal (step 6: No), step 2 is carried out so that the pouring ladle 10 continues to pour the molten metal from the pouring ladle 10 to the mold 20. Thereby the liquid density and the flow rate coefficient are updated every time the pouring ladle 10 pours the molten metal.

(25) If the pouring ladle 10 is supplied with the molten metal (step 6: Yes), a cycle of pouring has been completed and then step 7 takes place.

(26) At step 7, the relationship between the pouring start angle and the weight of the liquid within the pouring ladle is represented by the approximate function based on a plurality of sets of data acquired from respective data sequences. The sequences are the identified pouring start angle and the weight of liquid within the pouring ladle that were measured before pouring that were acquired at step 4 every time the pouring ladle 10 pours the molten metal.

(27) At step 8, the approximate function of the previous pouring start angle and the weight of the liquid within the pouring ladle is updated to the approximate function obtained at step 7. At a new cycle of pouring, that approximate function is used for the pouring control.

(28) Repeating the above process allows for rapid handling of a change in the pouring environment, and for the pouring control with a high degree of accuracy, depending on the pouring condition.

(29) Below, the mathematical model of a pouring process based on fluid mechanics that is used when a parameter identification technique is built is shown. As the pouring control system based on such a mathematical model, the inventors propose the model-based pouring control system that is shown in Patent documents 4 and 5. First, a mathematical model from the command signal u [V] to the tilting angle [rad] of the ladle that is used at step 2 of the pouring control is shown in Equation (2).

(30) $\begin{array}{cc}\left[\mathrm{Math}.2\right]& \\ \frac{d}{\mathrm{dt}}\left[\begin{array}{c}\left(t\right)\\ \left(t\right)\end{array}\right]=\left[\begin{array}{cc}0& 1\\ 0& -\frac{1}{T_m}\end{array}\right]\left[\begin{array}{c}\left(t\right)\\ \left(t\right)\end{array}\right]+\left[\begin{array}{c}0\\ \frac{K_m}{T_m}\end{array}\right]u\left(t\right),& \left(2\right)\end{array}$

(31) where equation (2) shows the speed control mode, [rad/s] denotes the tilting angular speed of the ladle, T.sub.m [s] denotes a time constant of the motor system, and K.sub.m [m/s/V] denotes a gain constant. When the servomotors are in a position control mode, the equation is represented in the form of equation (2), to which the position feedback mechanism is added.

(32) The mathematical model from the tilting angular speed of the ladle to the pouring flow rate q.sub.c [m.sup.3/s] is represented by equation (3) and equation (4).

(33) $\begin{array}{cc}\left[\mathrm{Math}.3\right]& \\ \frac{\mathrm{dh}\left(t\right)}{\mathrm{dt}}=-\frac{q_c\left(h\left(t\right)\right)}{A\left(\left(t\right)\right)}-\left\{\frac{h\left(t\right)}{A\left(\left(t\right)\right)}\frac{A\left(\left(t\right)\right)}{\left(t\right)}+\frac{1}{A\left(\left(t\right)\right)}\frac{V_s\left(\left(t\right)\right)}{\left(t\right)}\right\}\left(t\right),\left(\left(t\right){}_s\right)& \left(3\right)\\ \left[\mathrm{Math}.4\right]& \\ q_c\left(t\right)=c{}_0^{h\left(t\right)}{L}_f\left(h_b\right)\sqrt{2{\mathrm{gh}}_b}{\mathrm{dh}}_b,\left(h\left(t\right)0\right)& \left(4\right)\end{array}$

(34) As is shown in FIG. 4, the symbol h[m] in equation (3) shows the level of the liquid above the lip of the pouring ladle. The symbol A [m.sup.2] denotes the surface area of the upper surface of the liquid within the pouring ladle, and V.sub.s [m.sup.3] denotes the volume of the part of the liquid that is lower than the lip of the pouring ladle. The symbol [rad] denotes the tilting angle of the pouring ladle. Equation (3) is useful when the upper surface of the liquid within the pouring ladle is located above the lower surface of the lip of the pouring ladle, and when the tilting angle [rad] is equal to or larger than the tilting angle .sub.s [rad] of the ladle when the liquid within the pouring ladle begins to flow out. The ladle tilting angle denotes the pouring start angle. Also, L.sub.f [m] of equation (4) represents the width of the lip of the pouring ladle at the depth h.sub.b [m] of the liquid in the pouring ladle from its surface as shown in FIG. 5. The symbol g [m/s.sup.2] denotes the acceleration of gravity. The symbol c denotes the flow rate coefficient. Equation (4) is useful when the height of the liquid within the pouring ladle is above the lower surface of the lip of the pouring ladle.

(35) Equation (5) shows the relationship between the outflow weight W [kg] and the flow rate q.sub.c [m.sup.3/s] of the molten metal.

(36) $\begin{array}{cc}\left[\mathrm{Math}.5\right]& \\ \frac{\mathrm{dW}}{\mathrm{dt}}=q_c,& \left(5\right)\end{array}$
where the symbol [kg/m.sup.3] shows the liquid density. The outflow weight W [kg] is measured by the load cell built in the automatic pouring device 1. The response delay in the load cell is represented using the first order lag of equation (6).

(37) $\begin{array}{cc}\left[\mathrm{Math}.6\right]& \\ \frac{{\mathrm{dW}}_L\left(t\right)}{\mathrm{dt}}=-\frac{1}{T_L}{W}_L\left(t\right)+\frac{1}{T_L}W\left(t\right)& \left(6\right)\end{array}$

(38) Where the symbol W.sub.L [kg] is the outflow weight measured by the load cell, and the symbol T.sub.L [s] denotes the time constant corresponding to the response in the load cell.

(39) Equations (2) to (6) are represented as a mathematical model of the automatic pouring device 1. The tilting angle [rad] of the ladle is detected by the rotary encoder, and the outflow weight W.sub.L [kg] is detected by the load cell. The pouring control system is built using the mathematical model of this automatic pouring device 1. When the feedforward-type pouring flow rate control is carried out using the inverse model, if the desirable pouring flow rate pattern q.sub.cref [m.sup.3/s] is given, the inverse function of equation (4) allows the height of the liquid h.sub.ref [m] to be obtained that can achieve the desirable pouring flow rate pattern shown in equation (7).

(40) [Math. 7]
h.sub.ref(t)=.sup.1(q.sub.cref(t))(7)

(41) Here, we can adopt a technique of obtaining the inverse function of equation (7) by applying polynomial approximation to the inverse function of equation (4) and/or by making equation (4) be adapted to finite dimensions and linearly-interpolating values between elements in order to derive equation (7).

(42) The tilting angular speed .sub.ref [rad/s] of the ladle that achieves a desirable pouring flow rate pattern q.sub.cref [m.sup.3/s] can be obtained by substituting the obtained height of the liquid h.sub.ref [m] into equation (8), derived from equation (3).

(43) $\begin{array}{cc}\left[\mathrm{Math}.8\right]& \\ {}_{\mathrm{ref}}\left(t\right)=-\left\{\frac{{\mathrm{dh}}_{\mathrm{ref}}\left(t\right)}{\mathrm{dt}}+\frac{q_{\mathrm{conf}}\left(t\right)}{A\left({}_{\mathrm{ref}}\left(t\right)\right)}\right\}{\left\{\frac{h_{\mathrm{ref}}\left(t\right)}{A\left({}_{\mathrm{ref}}\left(t\right)\right)}\frac{A\left({}_{\mathrm{ref}}\left(t\right)\right)}{\left(t\right)}+\frac{1}{A\left({}_{\mathrm{ref}}\left(t\right)\right)}\frac{V_s\left({}_{\mathrm{ref}}\left(t\right)\right)}{\left(t\right)}\right\}}^{-1}& \left(8\right)\end{array}$

(44) Reference tilting angle .sub.ref [rad] in equation (8) can be obtained from equation (9), where equation (2) is used. .sub.sref [rad] in equation (9) denotes the pouring start angle. It is the tilting angle of the ladle at which the liquid begins to flow out of the pouring ladle.

(45) [Math. 9]
.sub.ref(t)=.sub.0.sup.t.sub.ref(t)dt+.sub.sref(9)

(46) The tilting angular speed .sub.ref [rad/s] of the ladle obtained in equation (8) is realized by using the command signal u.sub.ref [V], which is derived using the inverse-model of the motor model shown in equation (2). The inverse-model of the motor model is shown in equation (10).

(47) $\begin{array}{cc}\left[\mathrm{Math}.10\right]& \\ u_{\mathrm{ref}}\left(t\right)=\frac{T_m}{K_m}\frac{d{}_{\mathrm{ref}}\left(t\right)}{\mathrm{dt}}+\frac{1}{K_m}{}_{\mathrm{ref}}\left(t\right)& \left(10\right)\end{array}$

(48) A feedforward-type pouring flow rate control can be built using equations (7) to (10). Here, in the feedforward-type pouring flow rate control, the height of the liquid h.sub.ref [m] is required to be twice-differentiable.

(49) When a two degree of freedom pouring flow rate control into which the feedforward control and the feedback control are incorporated is built, the two degree of freedom pouring flow rate control can be built as one technique, based on the flatness shown below. If the flat output F is the height of the liquid h, the feedback linearization mechanism of equation (11) is built based on equation (3).

(50) $\begin{array}{cc}\left[\mathrm{Math}.11\right]& \\ u\left(t\right)=-\frac{1}{K_m}\left\{A\left(\left(t\right)\right)\left(t\right)+q_c\left(F\left(t\right)\right)\right\}{\left\{\frac{A\left(\left(t\right)\right)}{\left(t\right)}F\left(t\right)+\frac{V_s\left(\left(t\right)\right)}{\left(t\right)}\right\}}^{-1}& \left(11\right)\end{array}$

(51) Here, assuming that the responsiveness of the motor is much better than that of the pouring process, u=K.sub.m can be represented without considering the dynamic characteristic of the motor. Thus, equation (11) can be obtained. Equation (11) allows the model from the new control input v to the height of the liquid h(=F) at the lip of the pouring ladle to be linearized as shown in equation (12).

(52) $\begin{array}{cc}\left[\mathrm{Math}.12\right]& \\ \frac{\mathrm{dF}\left(t\right)}{\mathrm{dt}}=\left(t\right)& \left(12\right)\end{array}$

(53) Thus, the feedback control mechanism in equation (13) is built for the new control input v.

(54) $\begin{array}{cc}\left[\mathrm{Math}.13\right]& \\ \left(t\right)=\frac{{\mathrm{dF}}^{*}\left(t\right)}{\mathrm{dt}}-K_p\left(F\left(t\right)-F^{*}\left(t\right)\right)-K_i\left(F\left(t\right)-F^{*}\left(t\right)\right)\mathrm{dt},& \left(13\right)\end{array}$
where the symbol F* denotes the desirable target height of the liquid (F*=h.sub.ref), and the symbols K.sub.p and K.sub.i are control parameters that adjust the performance of the following target value that makes the actual height of the liquid h follow the target height of the liquid h.sub.ref. The desirable pouring flow rate q.sub.cref is given. The height of the liquid h.sub.ref that achieves the desirable pouring flow rate can be obtained from equation (7). The two degree of freedom pouring flow rate control in equations (11) and (12) is carried out based on the height of the liquid h.sub.ref. Here, the height of the liquid h.sub.ref is required to be a once differentiable function when the two degree of freedom pouring flow rate control is carried out. Also, equation (11) is useful as well as the feedforward-type pouring flow rate control, when the tilting angle of the ladle is equal to or greater than the pouring start angle .sub.s.

(55) The two kinds of pouring flow rate controls shown in above are both model-based pouring flow rate controls, which are based on the mathematical model of the pouring process. Here, many of the model parameters are set depending on the shape of the pouring ladle. However, since the flow rate coefficient c depends on the characteristics of the liquid and the characteristics of the surface texture of the pouring ladle, the parameters need to be identified by experiments. Moreover, although the pouring start angle .sub.s can be obtained by deriving the volume of the liquid from the weight of the liquid within the pouring ladle before pouring and by using the volume of the liquid and the shape of the pouring ladle, a difference from the model due to the effects of the fluctuation of the shape of the ladle caused by accumulating slag could occur. Moreover, since the liquid density of the high-temperature molten metal is likely to fluctuate depending on the temperature, the molten metal is susceptible to the pouring environment. Then, as shown in FIG. 2, a technique of identifying the flow rate coefficient, the pouring start angle, and the liquid density can be built based on the outflow weight data on the liquid, data on the tilting angle of the ladle, and the command signal data, which are obtained by using the automatic pouring.

(56) The parameter identification at step 7 is carried out by minimizing the evaluation function in equation (14). Specifically, it is minimized by applying the Down-hill simplex method as an optimization technique to the evaluation function in equation (14). Here, when the Down-hill simplex method is used, the convergence of the parameter(s) is fast and the computational load can be small. Thus, the parameter update time can be preferably short. In addition, optimization techniques such as a genetic algorithm, or a sequential quadratic programming approach, can be adopted.

(57) [Math. 14]
{c.sub.id,.sub.sid,.sub.id}=arg min{.sub.0.sup.T(W.sub.Lex(t)W.sub.Lsim(t,c.sub.sim,.sub.ssim,.sub.sim)).sup.2dt+w.sub.1(c.sub.avgc.sub.sim).sup.2+w.sub.2(.sub.avg.sub.sim).sup.2} (14),

(58) where the symbol T [s] denotes the pouring motion time of the automatic pouring device 1 that pours the molten metal into one mold, W.sub.Lex [kg] denotes the weight data on the outflow from the pouring ladle that the automatic pouring device 1 obtains through the built-in load cell, W.sub.Lsim [kg] denotes the weight of the outflow that is obtained when the simulation is carried out through the mathematical model of equations (2) to (6) by using the command value sent to the motor and the ladle tilting angle that is measured by the rotary encoder. The symbols c.sub.sim, .sub.ssim, and .sub.sim denote the flow rate coefficient, the pouring start angle, and the liquid density, respectively, that were used in the simulation. The symbols C.sub.avg and .sub.avg denote averaged values of flow rate coefficients and liquid densities, respectively, that were used until the previous cycle, and are represented as equations (15) and (16), respectively.

(59) 0$\begin{array}{cc}\left[\mathrm{Math}.15\right]& \\ c_{\mathrm{avg}}=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}{c}_{\mathrm{id}}\left(k-i\right)& \left(15\right)\\ \left[\mathrm{Math}.16\right]& \\ {}_{\mathrm{avg}}=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}{}_{\mathrm{id}}\left(k-i\right),& \left(16\right)\end{array}$

(60) where the symbol k denotes the number of times a pouring is carried out, and N denotes the number of pourings to be averaged. When the flow rate coefficient and/or the liquid density of the liquid to be poured are constant, N can be set to the maximum number of the pourings. However, when high temperature molten metal is used, the flow rate coefficient and/or the liquid density, may vary, depending on the temperature characteristics. Thus, adjusting the number of N and deleing the identified data obtained by the past pouring allow the accuracy of the identified data to improve.

(61) The symbol w.sub.1 in equation (14) denotes a weight coefficient for controlling the variation of the flow rate coefficient for every pouring. The symbol w.sub.2 denotes a weight coefficient for controlling the variation of the liquid density for every pouring. Increasing these allows the variation of the flow rate coefficient and liquid density that are identified for every pouring to be low. Since an adjustment of the weight coefficient allows the effect on the flow rate coefficient and the liquid density to be adjusted, the parameter identification with a higher accuracy can be made possible and the accuracy of pouring can be improved. For example, when the effect of the temperature in the liquid density is significant, it is recommended that the value of w.sub.2 be set to be small.

(62) An identified pouring start angle .sub.sid [rad] is combined with the weight of the liquid within the pouring ladle W.sub.b [kg] before the pouring that is measured by the load cell to be a set, and is stored as a set of the identified pouring start angle and the weight of the liquid within the pouring ladle in the computer 14. The molten metal can generally be poured a plurality of times from the automatic pouring machine that is supplied with the molten metal once. The pouring start angle can be estimated from the weight of the liquid within the pouring ladle measured before the pouring by making the approximate function using the data sequence of the pouring start angles .sub.sid=(.sub.sid (1), .sub.sid (2), . . . .sub.sid (n)) that are identified for every pouring and the data sequence of the weight of liquid within the pouring ladle before pouring W.sub.b=(W.sub.b (1), W.sub.b (2), . . . W.sub.b (n)). The linear approximation and/or the polynomial approximation are often used as an approximate function.

(63) In addition, the present invention can be applied to the non-transitory medium. It is readable by a computer in which a pouring control program is stored. The program causes the computer to carry out the above-explained process. That is to say, the present invention can be applied to a non-transitory medium that is readable by a computer in which a program is stored. The program causes the computer to carry out a process for controlling pouring based on a mathematical model of a pouring process from the input of at least one control parameter to pouring of molten metal using a pouring ladle in an automatic pouring device with a tilting-type ladle that pours the molten metal into a mold by tilting the pouring ladle that holds the molten metal. The process comprises the following:

(64) identifying, using an optimization technique, a flow rate coefficient, a liquid density, and a pouring start angle that is a tilting angle of the pouring ladle at which a flow out of the molten metal starts, wherein the flow rate coefficient, the liquid density, and the pouring start angle are the control parameters in the mathematical model, based on the weight of the liquid that flows out of the pouring ladle and ladle tilting angle that are measured during pouring, and a command signal that controls the tilting of the pouring ladle, and updating the control parameters to the identified control parameters.
[Effects of the Embodiments]

(65) The pouring control method of the present invention includes a pouring control method for controlling pouring based on the mathematical model of the pouring process from the input of the control parameters to the pouring using the pouring ladle. As the method includes identifying and updating the flow rate coefficient, the liquid density, and the pouring start angle that are control parameters within the mathematical model using the optimization technique, the operation for identification of the parameters, which normally takes much time to complete, can take less time. And the control parameters can be updated to the value corresponding to the pouring condition. And the control can deal with changes in the pouring conditions. Thus, the accuracy of pouring can be improved.

(66) Further, since the mathematical model of the pouring process based on fluid mechanics has been derived, and a model-based pouring control system has been adopted that is a pouring control system based on the model, the automatic pouring devices with a tilting-type ladle, each of which devices has a pouring ladle with a different shape and/or a different kind molten metal, can share the common parameter(s). Thereby the system can be booted in a short time and can carry out the pouring process analysis.

(67) Further, the present invention can be applied to a non-transitory medium that is readable by a computer in which the pouring control program is stored, where the program causes the computer to carry out the above explained process.

Examples of Experiments

(68) We carried out the experiments of pouring to indicate the usefulness of the pouring control method of the present invention. The experiment conditions are the following:

(69) Shape of pouring ladle: Sector form pouring ladle

(70) Used liquid: water

(71) Target outflow weight: 1.55 kg

(72) Target pouring flow rate (stationary time): 510.sup.4 m.sup.3/s

(73) Pouring control: feedforward-type pouring flow rate control

(74) A weight coefficient w.sub.1: 3

(75) A weight coefficient w.sub.2: 0.01

(76) The experimental results are shown in FIGS. 6 and 7. FIG. 6 shows the result of the first time pouring experiment. The flow rate coefficient and the liquid density are given appropriately and the pouring start angle, corresponding to the weight of the liquid within the pouring ladle obtained from the drawing of the shape of the pouring ladle, is used. FIG. 7 shows the result of the fourth pouring experiment. The pouring control is carried out after the parameters are identified and updated. After being poured three times, the liquid is again supplied into the pouring ladle. FIG. 6 (A) and FIG. 7 (A) show the ladle tilting angle measured by the rotary encoder, and FIG. 6 (B) and FIG. 7 (B) show the outflow weight measured by the load cell. The solid line shows the experimental result, and the dashed line shows the simulation result obtained using the mathematical model of the pouring process.

(77) In the first experiment for pouring, shown in FIG. 6, with regard to the initial parameters used for the pouring control, the flow rate coefficient is 0.98, the liquid density is 110.sup.3 [kg/m.sup.3], and the pouring start angle is 21.70/180 [rad]. On the result of the parameters that are identified after the experiment of the first pouring, the flow rate coefficient is 0.98, the liquid density is 110.sup.3 [kg/m.sup.3], and the pouring start angle is 20.20/180 [rad]. A comparison before and after the parameter identification shows that the difference is small for the flow rate coefficient and the liquid density, but the difference is large for the pouring start angle.

(78) This difference of the pouring start angles affects the difference between the simulation result of the outflow weight and the experimental result shown in FIG. 6 (B). In the fourth pouring, where the pouring control was carried out after the parameter shown in FIG. 7 was identified and updated, the flow rate coefficient used for the pouring control was 0.99, the liquid density was 110.sup.3 [kg/m.sup.3], and the weight of the liquid within the pouring ladle was 5.58 kg. Thus, 30.86/180 [rad] was used as the estimated value of the pouring start angle. When the parameter was identified after the fourth pouring experiment, the flow rate coefficient that was used for the pouring control was 0.99, the liquid density was 110.sup.3 [kg/m.sup.3], and the pouring start angle was 30.90/180 [rad]. Since the flow rate coefficient, the liquid density, and the pouring start angle that were used for the pouring control were almost the same as those of the results of the parameter identification, and the parameters that are suitable for the pouring condition were used for the pouring control, it was confirmed that the result of the experiment matched that of the simulation, and that the liquid was poured with a high degree of accuracy.

(79) The relationship between the weight of the liquid within the pouring ladle before pouring and the pouring start angle is shown in FIG. 8. The dashed line shows the relationship between the weight of liquid within the pouring ladle and the pouring start angle that is obtained using the pouring ladle. The black circle mark shows the identified pouring start angle and weight of the liquid within the pouring ladle before pouring. The solid line shows the results of identification, which is approximately linear. The linearly approximated relationship between the weight of liquid within the pouring ladle and the pouring start angle is shown in equation (17).

(80) [Math. 17]
.sub.s=4.046W.sub.b+53.4332(17)

(81) In the fourth experiment of pouring, the pouring start angle is predicted using the linearly approximated relationship between the weight of the liquid within the pouring ladle before pouring and the pouring start angle. It is found from FIG. 8 that the pouring start angle obtained using the figure of the shape of the pouring ladle is much different than the pouring start angle obtained using the parameter identification. This difference is considered to be caused by an error in modeling due to the simplification of the shape when the pouring start angle is derived from the figure of the shape of the pouring ladle and the change over the years of the shape of the pouring ladle. The pouring control method of the present invention allows us to grasp the relationship between the accurate pouring start angle and the weight of the liquid within the pouring ladle before pouring, and to employ it for the pouring control.

(82) As shown above, it was confirmed that pouring with a high degree of accuracy can be achieved by using the pouring control method of the present invention.

DESCRIPTION OF THE REFERENCE NUMERALS

(83) 1 an automatic pouring device 10 a pouring ladle 10a a lip of the pouring ladle 11, 12, 13 servomotors 14 a computer 20 a mold 20a a sprue of the mold