Method for filling a container with a fill product using a proportional valve
10112816 ยท 2018-10-30
Assignee
Inventors
Cpc classification
B67C3/286
PERFORMING OPERATIONS; TRANSPORTING
B67C3/007
PERFORMING OPERATIONS; TRANSPORTING
International classification
B65B1/04
PERFORMING OPERATIONS; TRANSPORTING
B67C3/00
PERFORMING OPERATIONS; TRANSPORTING
Abstract
Methods for filling a container with a fill product using a proportional valve are described. The method includes determining a function of the volume flow (q(t)) for a specified opening curve (H(t)) of the proportional valve or a specified required volume flow, and controlling the filling of the container based on the function of the volume flow (q(t)).
Claims
1. A method for filling a container with a fill product, comprising: determining a volume of fluid flowing through a proportional valve per unit time for each of a plurality of open positions of the proportional valve for a specified opening curve (H(t)) of the proportional valve to generate a curve of volume of the fluid flowing through the proportional valve as a function of time; and controlling the filling of the container based on the curve.
2. The method of claim 1, wherein determining the volume of fluid flowing through the proportional valve per unit time comprises determining a function (Kv(H)) of a flow factor (Kv) of the proportional valve for a plurality of open positions (H) of the proportional valve and a volume flow (q.sub.?(H)) in a steady state.
3. The method of claim 2, further comprising controlling subsequent filling processes using the determined volume of fluid flowing through the proportional valve per unit time or calculating a new volume flow as a function of time (q(t)) for a new filling process on the basis of the function (Kv(H)) of the flow factor.
4. The method of claim 2, wherein the volume flow (q.sub.?(H)) in the steady state is determined by:
5. The method of claim 2, wherein determining a function (Kv(H)) of the flow factor (Kv) of the proportional valve for the plurality of open positions (H) comprises measuring the flow factor (Kv) at discrete open positions H and determining a function of the flow factor (Kv(H)) from discrete flow factor values.
6. The method of claim 5, further comprising fitting the function (Kv(H)) of the flow factor to the discrete flow factor values using regression analysis, the least squares method, a fit algorithm, or determination of a polynomial equation.
7. The method of claim 6, wherein the polynomial equation is of the sixth order.
8. The method of claim 5, wherein measuring the flow factor (Kv) at an open position (H) of the volume flow (q) comprises determining a differential pressure (?p) and a density (?) of a medium flowing through the proportional valve.
9. The method of claim 8, further comprising calculating the flow factor (Kv) by:
10. The method of claim 8, wherein the differential pressure (?p) is determined using one or more of pressure sensors in a base of a supply vessel and addition of a static height to a valve outlet, level sensors in the supply vessel and addition of the static height to the valve outlet, and pressure sensors disposed at an outlet of a filling element.
11. The method of claim 1, wherein the curve is determined by a solution of a mechanical mesh equation.
12. The method of claim 11, wherein the solution of the mechanical mesh equation is solved for downward sloping open positions by
13. The method of claim 1, further comprising calculating a volume of fluid flowing through the proportional valve per unit time for a specified closing curve to end a filling process, calculating a lag volume by integration over the volume of fluid flowing through the proportional valve per unit time, and taking the lag volume into account when determining a time point at which the closing curve is traversed.
14. The method of claim 13, wherein the lag volume is calculated for a plurality of filling processes before an estimated end of filling.
15. The method of claim 1, further comprising determining a total fill volume of the fill product by integration of the curve over a course of a filling process, and determining an end of filling based on the determined total fill volume.
16. A method for filling a container with a fill product, comprising: determining a volume flow as a function of time (q(t)) for a specified opening curve (H(t)) of a proportional valve; calculating an actual volume flow (q.sub.actual-calculated) based on the determined volume flow as a function of time (q(t)) for an actual open position (H.sub.actual) of the proportional valve; and controlling a filling of the container with a required volume flow curve (q.sub.req(t)) as a reference variable and the calculated actual volume flow (q.sub.actual-calculated) as a controlled variable.
17. The method of claim 16, wherein determining the volume flow as a function of time (q(t)) comprises determining a function (Kv(H)) of a flow factor (Kv) of the proportional valve for open positions (H) of the proportional valve and a volume flow (q.sub.?(H)) in a steady state, wherein the volume flow (q.sub.?(H)) in the steady state is determined by:
18. The method of claim 17, further comprising calculating the flow factor (Kv) by:
19. The method of claim 16, wherein the volume flow as a function of time (q(t)) is determined by a solution of a mechanical mesh equation, and the solution of the mechanical mesh equation is solved for downward sloping open positions by
20. The method of claim 16, further comprising determining a total fill volume of the fill product by integration of the volume flow as a function of time (q(t)) over a course of a filling process, and determining an end of filling based on the determined total fill volume.
Description
BRIEF DESCRIPTION OF THE FIGURES
(1) Further embodiments and aspects of the present invention are more fully explained by the description below of the figures.
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
DETAILED DESCRIPTION
(10) Examples of embodiments are described below with the aid of the figures. In the figures, elements which are identical or similar, or have identical effects, are designated with identical reference signs, and repeated description of these elements is in part dispensed with in the description below, in order to avoid redundancy.
(11)
(12) Such a curve of the measured volume flow 1 results, for example, when a proportional valve is controlled with a specified opening curve, and accordingly differing open positions or ramp functions of the open position of the proportional valve are traversed during the filling process. The curve of the measured volume flow 1 that is shown in
(13) It can be seen in this example that an opening process 10 traverses a short plateau region 12 followed by an opening ramp 14. The main filling is then carried out in a plateau region 16 of high volume flow, after which a closing curve 18 is traversed. The example shown here is to be understood as an example only. It shows, however, a common behavior of the volume flow 1 when beverages are filled and the aim is to achieve as gentle a filling process as possible, with a low tendency to foaming and exactly the desired fill volume.
(14) It can also be recognized from the measured curve of the volume flow 1 that there is no sudden increase in the volume flow during the opening process 10, or sudden cessation of the volume flow when the closing curve 18 is traversed. Instead, the applicable curves have a finite gradient or a complex shape. This is because, first, a proportional valve needs a finite time for the opening and closing process, and secondly that the medium flowing through has a certain inertia, which prevents the medium from responding immediately to the open position of the proportional valve.
(15) The determination of the curve of the volume flow 1 by means of a specified opening curve of the proportional valve, which is described below, enables a volume flow curve 1 also to be determined without actual measurement of the volume flow. Instead, the volume flow curve 1 can be calculated or simulated by means of the Kv model that is described below. Accordingly, it is possible in normal filling operation to dispense with the measuring element in the form of a flow meter or a weighing cell that was described above, with the result that a filling device that is controlled in this manner can be designed in a simpler, more cost-efficient and more reliable manner.
(16) In order to enable the determination of the volume flow curve q(H(t)) from the opening position curve H(t) of the proportional valve, so that in actual filling operation it is possible to dispense with the use of a control system and in particular a measuring element, the following method is described.
(17) First a function of the flow factor Kv(H) of the proportional valve for every open position H of the proportional valve is determined. The flow factor Kv is also known as the flow coefficient. It is a measure of the attainable flow rate of a fluid or a gas through a proportional valve. It is here stated in the units ml/sec and can be interpreted as the effective cross-section. Each Kv value applies only to the associated open position H of the proportional valve.
(18) In order to determine the flow factor Kv, in an initial calibration process the proportional valve is controlled to adopt a specific open position H.sub.i, the fill product flow q(H) from the proportional valve is measured at this open position H.sub.i, and from this the flow factor Kv is determined in the steady state, for example by means of measurement using a measuring cell such as a weighing cell. This process is carried out for a multiplicity of discrete open positions H; of the proportional valve.
(19) The following relationship is obtained between the Kv value and the volume flow q.sub.? (volume flow in the steady state):
(20)
(21) where ?p is the differential pressure between the valve outlet and the static level of the fill product;
(22) where ? is the density of the fill product.
(23) Thus for the exact determination of the flow factor Kv, it is also necessary, besides the above-mentioned measurement of the volume flow at a specific open position, to determine the differential pressure ?p and the density ? of the medium flowing through the proportional valve.
(24) The differential pressure ?p between the valve outlet and the static level varies with the static level of the fill product, which can be accommodated for example in a supply vessel. In a beverage filling plant, for example, the following alternatives exist for determining the differential pressure.
(25) The total pressure difference can be determined by means of pressure sensors in the base of the supply vessel and addition of the static height to the valve outlet.
(26) The total fluid column bearing on the valve outlet can be determined by means of level sensors in the supply vessel and addition of the static height to the valve outlet. From this, the total pressure difference can be determined. This method of determining the total pressure difference may, however, have the disadvantage that it may not be possible to take account optimally of, for example, waves in the fill product in the vessel, and centrifugal effects brought about by differing filling rates.
(27) The differential pressure can be determined exactly by means of pressure sensors disposed directly at the outlet of the filling element. However, the pressure sensor may influence the flow behavior of the fill product at the filler valve outlet, and hence the filling process itself.
(28) The density ? of the fill product is usually known, or can be determined by known measurement methods. For water the density can be assumed to be approximately 1000 kg/m.sup.3, and then has no influence on the calculation.
(29) Accordingly, from the volume flow q measured for a specific open position H.sub.i, the specific differential pressure ?p and the specific ?, the Kv value for this open position can now be determined by:
(30)
(31) In order here to determine a function of the flow factor Kv(H) over the open positions H.sub.i, after the determination of all flow factors Kv(H.sub.i), a function of the flow factor over the open positions of the proportional valve is determined by means of the determination of a best-fit curve through the each of the flow factors Kv(H.sub.i). The best-fit curve can be determined for example by linear regression, the least squares method, a fit algorithm, or another known method for determining a best-fit curve through measured values. This determination and calculation is carried out for various discrete values of the open position H.sub.i.
(32) As the best-fit curve, a polynomial of the sixth order can for example be used, as shown for example in
(33) For a certain lift H of the proportional valve the resulting best-fit curve of the flow factor Kv is for example:
Kv(H)=c.sub.6*H.sup.6+c.sub.5*H.sup.5+c.sub.4*H.sup.4+c.sub.3*H.sup.3+c.sub.2*H.sup.2+c.sub.i*H+c.sub.7(3)
(34) where c.sub.1 to c.sub.7 are the respective coefficients for adjusting the function to the measured values.
(35) By means of the determination of the best-fit function, all intervening values of the open positions can then be taken into account in the filling process. The corresponding volume flow can thereby be calculated for steady states for every open position:
(36)
(37) In this, however, it must be taken into account that this function of the flow factor Kv(H) of the proportional valve for each open position relates to the corresponding volume flow in the steady state, i.e., after an extended period in which the open position of the proportional valve has remained constant. When the proportional valve opens, closes or travels from one open position to another, however, other dynamic influences become important.
(38) In order to address the dynamic influences, an analogy from the field of electrical engineering will first be used. In the table below certain electro-mechanical analogies are drawn:
(39) TABLE-US-00001 Electrical equivalent Mechanical equivalent ohmic resistance Kv value voltage differential pressure current volume flow inductance accelerated mass
(40)
(41) The power source, which provides a voltage U, thus represents the differential pressure which provides a corresponding pressure for the fluid. The current I corresponds to the volume flow, and the inductance L the accelerated mass.
(42) The dynamic behavior of this electrical circuit can be represented by means of a mesh equation, which involves an inhomogeneous differential equation of the first order:
(43)
(44) This inhomogeneous differential equation of the first order is satisfied by the following solution for increasing resistance values, i.e., continually added Rp1-n:
(45)
(46) where i.sub.0=current at the beginning of the settling process.
(47) It should however be noted that this settling process is constantly interrupted and replaced by a new settling process due to continually changing Rp.
(48) A similar solution of the mesh equation can be found for decreasing resistance values, i.e., continually deactivated Rp1-n:
(49)
(50) where i.sub.0=current at the beginning of the settling process.
(51) It should however be noted that this settling process too is constantly interrupted and replaced by a new settling process due to continually changing Rp.
(52) The electrical picture can be transferred to the mechanical equivalent, resulting in the following interrelationships for decreasing lifts H (corresponding to increasing resistance values):
(53)
(54) where q.sub.0=volume flow at the beginning of the settling process
(55) where m=constant.
(56) The following interrelationships result for increasing lifts H (corresponding to decreasing resistance values):
(57)
(58) where q.sub.0=volume flow at the beginning of the settling process
(59) where m=constant.
(60) These formulae are valid in each case for a defined incremental change in the open position. If the open position changes again during the settling process, the settling process is replaced by a new settling process with new parameters.
(61) For the curve of the volume flow 1 that is shown in
(62) Here it can immediately be seen that the measured volume flow 1 is substantially identical to the simulated volume flow 3. Accordingly, when the differential pressure, the density of the fill product and the flow factor are known, the simulated volume flow 3 can replace the continual measurement of the volume flow 1 during the filling operation.
(63) Thus if the described method is used, it is possible to dispense with a measuring element for controlling the filling process on the basis of the measured volume flow. Instead, after the calculation of the simulated volume flow 3, the proportional valve can be controlled directly in a simple manner to adopt a particular opening curve, and the desired fill volume can be reached. Thus a valve controller can be used in place of the control system.
(64) In order to determine the fill volume that is introduced into the container that is to be filled, the simulated volume flow curve 3 is then integrated.
(65) At this point a special situation at the end of the filling process is additionally considered. As already discussed above, the closing curve of a proportional valve is relatively prolonged, and the closing process can last as long as 500 milliseconds. For this reason, as well as the relative inertia of the volume flow, the volume flow is only slowly switched off, which results in a lag that must be taken into account in the calculation of the total fill volume.
(66) In this context,
(67)
(68) On the basis of the method described above for determining the simulated volume flow for a specified opening curve, the volume flow can also be calculated by similar means for the specified closing curve 40, i.e., a curve of the open positions H(t) of the proportional valve as it closes. By integration of the volume flow over the entire closing curve, the total lag volume thus can be determined.
(69) In
(70) In a variant shown in
(71) The calculation is performed as near to the expected switch-off time point as possible, in order to carry out the calculation of the lag volume in conditions that are as close as possible to those which actually obtain during the lag, for example as regards the fluid level in the supply vessel, the corresponding pressure in the vessel, or the product temperature, which in turn affects the viscosity of the fill product.
(72) In order to calculate the lag, the sequence of steps listed below is carried out. The calculation, in various embodiments, takes place at the point in time at which the proportional valve has reached the maximum volume flow, and enough time still remains before the end of the filling process:
(73) (1) Determination of the duration of the movement, from the acceleration of the motor and the maximum motor speed, in order to determine the closing curve H(t);
(74) (2) Determination of the prevailing static level in the region of the plateau 16 from:
(75) Kv(H) and the current volume flow q; and
(76) The static level is needed during the calculation, for the continual calculation of q.sub.?;
(77) (3) Simulated (calculated) incremental lowering of the open position H in accordance with the movement profile as a function of the time t;
(78) (4) Calculation of Kv(H) at the new open position;
(79) (5) Calculation of q.sub.? and ? with Kv(H);
(80) (6) Calculation of q(t)
(81) (7) .fwdarw.Go to step 3, or end the calculation if the duration of movement t has expired.
(82) By means of the calculation of the lag, it is possible to achieve greater accuracy in the case of fluctuations in the level of the fill product in the vessel, fluctuations in the pressure in the vessel, and fluctuations in temperature. It is thereby also possible to achieve optimum filling of the first container, and dispense with control of the lag using a measuring element.
(83)
(84) A required volume flow 50 is supplied as reference variable to the control system, which is embodied for example as a PID controller. The required volume flow 50 can also be a required volume flow curve q.sub.req(t), i.e., a volume flow which can vary over time, which is required to flow into a container that is to be filled.
(85) On this basis therefore, by means of the PID controller a required open position H.sub.req is determined, and passed as the required open position 52 to the actuator of the proportional valve. On the basis of this control command the proportional valve opens, and passes its actual open position H.sub.actual as the actual open position 54 to the Kv calculation model as the open position. Due to the inertia of the proportional valve, the actual open position 54 and the required open position 52 can diverge.
(86) Based on the actual open position H.sub.actual, the Kv model calculates a calculated actual volume flow q.sub.actual-calculated as the calculated actual volume flow 56, which is again supplied to the control system as the controlled variable, and compared with the reference variable of the required volume flow q.sub.req. Thus the controlled variable q.sub.actual-calculated is calculated based on the actual open position H.sub.actual, and is then used to control the required open position H.sub.req. In contrast to conventional control, it is possible in this manner to dispense with a measuring element for determining the controlled variable, i.e., the actual volume flow. The measuring element is instead replaced by the Kv model.
(87) As discussed and demonstrated above, by means of the Kv model the actual volume flow can be determined with good accuracy for a specific open position. The calculated actual volume flow q.sub.actual-calculated thus substantially corresponds to the actual volume flow q.sub.actual that is actually flowing through the proportional valve, which could be measured as the actual volume flow 58.
(88) To the extent applicable, all individual features described in the individual example embodiments can be combined with each other and/or exchanged, without departing from the field of the invention.