The present invention relates generally to clothes washing machines, and more particularly, to a method and system for controlling the washing machine drum.
Residential and commercial clothes washing machines are well known. A generally cylindrical drum or basket for holding the clothing and other articles to be washed is rotatably mounted within a cabinet. Typically, an electric motor drives the drum. During is a wash cycle, water and detergent or soap are forced through the clothes to wash them. The detergent is rinsed from the clothes, then during one or more spin cycles the water is extracted from the clothes by spinning the drum.
One way of categorizing washing machines is by the orientation of the washing machine drum. Vertical-axis washing machines have the drum situated to spin about a vertical axis. Articles to be washed are loaded into the drum through a door, which is usually situated on the top of the washing machine. A vertical-axis washing machine drum includes an agitator situated therein, which cleans clothes by pushing and pulling them down into the water. A motor typically drives the agitator, in addition to spinning the vertically-oriented drum during spin cycles. The motor usually operates at a constant speed, and a series of gears or belts are configured to drive the proper component at the proper time during each washing machine cycle.
Horizontal-axis washing machines, having the drum oriented to spin about an essentially horizontal axis, do not include an agitator, and a variable-speed motor drives the drum. During wash cycles, the drum of the horizontal-axis washing machines rotates at a relatively low speed. The rotation speed of the drum is such that clothes are lifted up out of the water, using baffles distributed about the drum, then dropped back into the water as the drum revolves.
Both vertical and horizontal-axis washing machines extract water from clothes by spinning the drum, such that centrifugal force extracts water from the clothes. It is desirable to spin the drum at a high speed and extract the maximum amount of water from the clothes in the shortest possible time. Spin time is reduced, but more power is required to spin at a higher speed. The distribution of the clothes about the periphery of the drum affects the washing machine's ability to spin the drum at a high speed.
If the clothes in the drum are not distributed evenly, the drum will not be able to rotate at an acceptable spin speed without excessive drum vibration in the cabinet. Several factors contribute to drum load unbalance problems. For instance, with vertical-axis washing machines, when a wash or rinse cycle completes and the water is drained from the drum, the clothes are gathered at the bottom of the drum, and are not evenly distributed about the entire drum. Moreover, the drum typically is not perfectly cylindrical; but rather, includes a draft. When the drum spins, the clothes will “creep” up the sides of the drum. However, since a constant speed motor typically drives the vertically-oriented drum, the motor quickly ramps the drum up to the full spin speed. There is little chance for the clothes to distribute about the periphery of the drum, so they creep up the drum's sides in an unbalanced fashion.
The unbalanced, spinning drum vibrates within the cabinet. In conventional vertical-axis washing machines, if the vibration is too severe, the drum will trip a switch mounted inside the cabinet, stopping the drum's rotation and activating a drum-unbalance alarm. A user then manually redistributes the wet clothes within the drum, and restarts the spin cycle.
As discussed above, the drum in a horizontal-axis machine is driven by a variable speed motor. This allows the inclusion of a “distribution” cycle, wherein the drum is rotated faster than the rotation speed of a wash cycle, but slower than in a spin cycle. The drum rotation speed is gradually increased, until the clothes begin to “stick” to the sides of the drum due to centrifugal force. The slower rotation speed allows the clothes to more evenly distribute about the sides of the drum. Once the clothes have been distributed about the drum, the speed is increased to a full spin speed to extract the water from the clothes.
Horizontal-axis washing machines, however, are not immune to drum unbalance problems. If the clothes do not evenly distribute during the distribution cycle, the unbalanced load within the drum will cause unwanted vibrations as the drum rotates. Rather than applying all of the motor's power to spinning the drum at the highest possible speed, power is wasted in drum movement and cabinet vibrations. Detecting the amount of load unbalance allows the spin speed to be optimized to give an efficient wash, while at the same time minimizing the vibrations caused by centrifugal unbalance forces. In cases where the amount of unbalance detected is very high, the washing machine can be programmed to stop its spin cycle, shake the drum to redistribute the washing, and then restart the spin cycle.
Generally, washing machine spin cycles are pre-programmed to run for a fixed duration. In some machines, there may be several pre-programmed spin times stored in memory with the appropriate one being determined by the type of wash selected by the user. When the spin cycle is set to run for a pre-programmed duration, the spin time is likely to be set conservatively in order to ensure that the laundry is spun dry in most possible circumstances. As a result, the spin cycle may last longer than is actually required.
Therefore, it is desirable to automatically detect the end of a spin cycle allows the duration of the spin cycle to be optimized, thus saving time and energy.
In one aspect of the present disclosure, a clothes washing machine system includes a cabinet having a drum rotatably mounted therein. A motor is coupled to the drum for rotating the drum within the cabinet. A memory device is accessible by a processor, such as a digital signal processor (DSP). The memory device and processor may be separate devices or alternatively the memory device may be embedded with the processor itself, as is the case with a microcontroller. The memory contains program code for controlling operation of the washing machine, including automatically determining the end of the washing machines spin cycle. The motor spins the drum, and a predetermined parameter of the spinning drum is monitored. The drum is stopped in response to the monitored parameter.
In an exemplary implementation, the predetermined parameter is the inertia of the drum and laundry contained in the drum. The parameter, such as inertia, may be monitored discretely or continuously, and the end of the spin cycle may be detected based on the parameter reaching a predetermined value, or the rate of change of the parameter reaching a predetermined level, for example. The predetermined parameter, such as inertia, may be determined using a parameter estimator such as a recursive least squares estimator.
Other objects and advantages of the invention will become apparent upon reading the following detailed description and upon reference to the drawings in which:
While the invention is susceptible to various modifications and alternative forms, specific embodiments thereof have been shown by way of example in the drawings and are herein described in detail. It should be understood, however, that the description herein of specific embodiments is not intended to limit the invention to the particular forms disclosed, but on the contrary, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention as defined by the appended claims.
Illustrative embodiments of the invention are described below. In the interest of clarity, not all features of an actual implementation are described in this specification. It will of course be appreciated that in the development of any such actual embodiment, numerous implementation-specific decisions must be made to achieve the developers' specific goals, such as compliance with system-related and business-related constraints, which will vary from one implementation to another. Moreover, it will be appreciated that such a development effort might be complex and time-consuming, but would nevertheless be a routine undertaking for those of ordinary skill in the art having the benefit of this disclosure.
Referring back to
The memory 108 stores program instructions executed by the processor 109 for controlling operation of the washing machine system 100. Among other things, the processor 109 is programmed to control the speed of the motor 106 and thus, the drum 104. A typical speed control loop 200 is shown in
Speed/position detection and control of the rotor 106b may be implemented, for example, based on feedback from a tachometer on the motor shaft. In other washing machine systems embodying aspects of the present invention, position/speed of the rotor 106b is calculated or estimated from monitored parameters of the motor 106, such as the motor's voltage and current. These systems are often called “sensorless” systems since they do not use a physical transducer to determine the position and or speed. In such sensorless systems, the processor 109 receives indications of the motor current and voltage, and an algorithm stored in the memory 108 determines the position/speed of the rotor 106b based thereon.
As discussed in the Background section herein above, washing machines typically include a variety of operation cycles. Washing machines—particularly horizontal-axis machines—include one or more wash cycles, distribution cycles and spin cycles. Drum unbalance is rarely a significant problem during wash cycles, which, in a horizontal-axis machine, use a drum rotation speed of about 50 rpm to tumble the clothes in and out of the water. Distribution cycles typically operate at a drum rotation speed of about 55-110 rpm (clothes will begin to “plaster” or “stick” to the sides of the drum 104 at one G of centrifugal force).
In comparison, the minimum rotation speed that is normally considered a “spin cycle” speed is about 250 rpm. For example, a drum rotation speed of about 350-450 may be considered a “low” spin speed, a drum rotation speed of about 650-850 may be considered a “medium” spin speed, and a drum rotation speed of about 1,000 rpm may be considered a “high” spin speed. As discussed above, it is desirable to rotate the drum 104 at a high speed to extract the maximum amount of water from the clothes in the shortest possible time.
Detecting the amount of a load unbalance allows the spin speed to be optimized to give an efficient wash, while at the same time minimizing the vibrations caused by centrifugal unbalance forces. Accordingly, the memory 108 includes program code that when executed by the processor, implements methods to determine an unbalance mass in the rotating drum 104.
An exemplary method for determining the unbalance mass is summarized in
In block 250 of
In block 254 the inertial torque caused by any speed ripple is subtracted from the overall torque, which ensures that the remaining torque ripple is due entirely to gravity acting on the unbalanced mass. The inertial torque ripple is calculated from the inertia calculated in block 250 and the speed.
The drum position and compensated torque are then fed into a parameter estimator, which gives an estimate of the torque ripple in block 256. Blocks 252, 254 and 256 are repeated at a plurality of sampling points until the torque ripple estimate converges to a steady value, typically within five drum cycles. In block 258, the unbalance mass is determined based on the torque ripple. The method summarized in
During the washing machine distribution cycle the laundry gets plastered around the periphery of the drum 104 in a random, and essentially unpredictable, manner. The distribution of the laundry 110 will invariably be uneven as shown in
In accordance with aspects of the present invention, a mathematical model of the washing machine system is constructed. Thus, in order to detect drum unbalance it is necessary to identify the system variables that will be affected by drum unbalance. As noted above, these variables include drum/motor torque and drum/motor angular velocity. When unbalance is present, the drum/motor torque and angular velocity have a cyclically varying component at drum frequency (there may be some higher frequency components as well but these are essentially the result of second order effects or “dynamic unbalance”). The amplitude of this cyclically varying component, or ripple, gives a measure of the amount of unbalance present. The added advantage of using these variables to detect unbalance is that both can be estimated using sensorless control, thus eliminating the need for measurement instrumentation.
Accordingly, the dynamic model of the washing machine drum system obtained contains variables such as torque and speed as well as acceleration and drum position, which can be determined from speed and vice versa. The coefficients of the variables in the dynamic model are composite of various system parameters such as unbalance mass, drum inertia and friction coefficients.
Assuming that the variables of the dynamic model can be obtained by sensorless techniques, or even by using onboard instrumentation, a parameter estimator can be employed to determine the coefficients of the variables and hence the required system parameters such as unbalance mass. A suitable parameter estimator for this purpose is a recursive least squares (RLS) estimator.
There is a certain amount of unbalance inherent in an empty drum 104, so the equivalent point mass unbalance 112 will be that of both the empty drum 104 and the laundry 110 combined. There is no advantage in trying to separate drum unbalance 112 from the laundry 110 unbalance because it is both of these unbalances acting together that ultimately results in the observed vibrations of the washing machine.
An exemplary process of constructing the dynamic model of the drum 104 in accordance with one embodiment of the invention is further detailed herein. The first stage in constructing such a model of the drum is to assume that the axis of rotation of the drum is fixed—the drum suspension is ignored. The fixed axis model is then developed in three steps: the first step is to assume that the axis of the drum is horizontal and that drum friction and drag losses are zero; the second step is to include the effects of friction and drag; and finally the third step is to include drum inclination.
Drum position θd is measured in a counterclockwise direction between a reference point 120 on the drum 104 and a fixed datum 122, which is arbitrarily chosen to be at 3 o'clock. The location of the point mass unbalance 112 is given by φu with respect to the drum reference point 120 in a clockwise direction. The position of the point mass unbalance 112 with respect to the fixed datum 122 is therefore (θd−φu).
Assuming for the moment that there are no friction and drag losses associated with the drum 104, the drum torque τd in the direction of motion is given by:
τd=mu·g·ru·cos(θd−φu)+Jd·P(ωd) (1)
where Jd is the combined inertia of the drum and the unbalance mass and p is a differential operator, d/dt.
It can be seen from Equation 1 that if the velocity ωd were constant, i.e. uniform circular motion, the drum torque τd would vary cyclically with a mean value of zero. Therefore, for the fixed speed case Equation 1 reduces to:
τd=mu·g·ru·cos(θd−φu)
In practice there will be some variation in drum velocity ωd, depending upon the response of the system's speed controller 210 implemented by the processor 109. If the response of the speed controller 210 is very sluggish, the torque τd will be constant and zero. Therefore, for the fixed (and zero) torque case Equation 1 reduces to:
−mu·g·ru·cos(θd−φu)=Jd·p(ωd)
As the mean velocity of the drum ωd increases, the magnitude of the speed ripple decreases—as can be seen by considering the energy system:
where the left side of the equation is the potential energy gained due to gravity in moving the unbalanced mass 112 from the bottom to the top of the drum 104. The right side is the gain in kinetic energy in accelerating the unbalance mass 112 from ωd(min) to ωd(max), where ωd(min) is the angular velocity of the drum 104 when the unbalance mass 112 is at the top of the drum 104 and ωd(max) is the angular velocity of the drum 104 when the unbalance mass 112 is at the bottom of the drum 104.
In practice, both torque and velocity will vary, and so the observed behavior of the system will lie somewhere between the two extremes outlined above. The amount of torque and speed ripple will depend upon the response of the speed control loop 200, which in turn usually depends upon the proportional and integral gains of the speed controller 210. Therefore, the two cases considered are not appropriate for developing the system model, so the general equation given by Equation 1 is used from here onwards. However, it is important to consider the two extremes to understand that speed and torque ripple are affected by the response of the speed loop, and not just the unbalance mass 112.
In practice, a certain amount of torque is required to keep the drum running at a given speed in order to overcome forces due to drum friction and wind drag. Two types of friction will be considered: coulomb and viscous friction. Coulomb friction is a constant and is therefore independent of speed, although its polarity depends on the direction of rotation of the drum. Viscous friction, on the other hand, is proportional to the angular velocity of the drum. Both of these frictional forces take into account the resistance to motion due to the bearings, lubrication and wind drag, where the latter is assumed to be a type of viscous friction.
The coulomb friction is given by τcsignum(ωd), where signum(ωd) denotes that its polarity is the same as that of the drum's angular velocity, and the viscous friction is given by Bd·ωd, where Bd is the drum's viscous friction coefficient. Expanding Equation 1 to include these frictional forces gives:
τd=mu·g·ru·cos(θd−φu)+Jd·p(ωd)+Bd·ωd+τcsignum(ωd) (2)
The drum 104 in some versions of “horizontal” axis washing machines is not actually horizontal, but rather, is inclined at some angle to the horizontal.
Force=mu·g·cos(α)
The drum 104 illustrated in
Equation 2 can now be modified to include drum inclination as follows:
τd=mu·g·ru·cos(α)·cos(θd−φu)+Jd·p(ωd)+Bd·ωd+τcsignum(ωd) (3)
Note that if π=90°, the gravitational force on the unbalance mass 112 acts along the axis of the shaft 140 and therefore has no effect on shaft torque.
Referring back to
In developing the unbalance detection system, assumptions for the washing machine transmission 300 system include the following:
Based on
τm=Jm·p(ωm)+Bm·ωm+(Tmu−Tm1)·rmp (4)
where Jm is the motor pulley 312 inertia, Bm is the motor viscous friction coefficient, Tmu and Tm1 are the upper and lower belt tension, respectively, at the motor pulley 312, rmp is the radius of the motor pulley 312, and ωm is the motor shaft speed. The torque applied to the drum pulley 310, τd, is given by:
τd=(Tdu−Td1)·rdp (5)
where Tdu and Td1 are the upper and lower belt tension, respectively, at the drum pulley 310 and rdp is the radius of the drum pulley 310. As noted above in assumption 2, it is assumed that the top or bottom of the belt 314 is either fully tensioned or fully slack. This leads to the following:
Tm1=Td1=0 (6)
Tmu=Tdu (7)
Therefore, combining Equations 4 and 5 gives:
Assuming a perfect transmission system, the drum speed and motor speed are related as follows:
Equation 8 can therefore be written in terms of drum speed as follows:
From the fixed axis model, it can be seen that τd contains the terms Jd·p(ωd) and Bd·ωd, which are similar to the inertial and viscous terms as those in Equation 10. The drum and motor inertial and viscous terms combine to give:
Now,
Equation 10 simply reduces to:
The fixed axis model given by Equation 3 is referred to the drum, i.e. torque and speed are both drum quantities. However, it is more likely that the torque and speed variables that form the input to the model will be motor quantities, so Equations 9 and 11 will be required to refer motor speed and torque to the drum.
In order for the fixed axis model to be valid, the drum must not move significantly on its suspension. Therefore when using this model unbalance detection must be carried out when the drum is above minimum plaster speed, but less than the speed at which the suspension has a significant impact. Accordingly, in embodiments of the invention, the drum is run at a predetermined speed, such as 100 rpm, for some time period while the unbalance detection is implemented. Even if the drum is commanded to run at a predetermined constant speed, there will be some speed ripple that will result in a small amount of torque ripple due to the viscous friction term, Bd. However, if this torque ripple component is assumed negligible and the speed is otherwise substantially constant, the fixed axis model given by Equation 3 reduces to:
τd=mu·g·ru·cos(α)·cos(θd−φu)+Jd·p(ωd)+τmean (12)
In order for a unique set of parameters to be estimated for the fixed axis system, the input signal must be sufficiently exciting. This requirement will be met if the input signal excites all modes of the system; such a signal is said to be “persistently exciting.” For the system considered here, the input signal is effectively the drum speed, ωd, as the other variables, i.e. θd and p(ωd), are themselves determined by the drum speed and vice versa.
During unbalance detection the drum will be commanded to run at a constant speed. If the actual drum speed is perfectly constant during this time there will be no angular acceleration of the drum, i.e. p(ωd)=0. Under this mode of operation the inertia has no effect on the drum torque so it would not be possible for the parameter estimator to get an estimate of the drum inertia, Jd. However, in practice the speed controller 210 will not be capable of regulating the drum speed to a perfectly constant value and there will be some speed variation of the drum about the commanded speed 202.
Test results show that, at speeds pertaining to a fixed axis system, the speed ripple is approximately sinusoidal at a frequency corresponding to drum speed. At small levels of speed ripple, the unbalance torque will be approximately sinusoidal because cos(ωd·t−φu) ≈cos({overscore (ω)}d·t−φu), where {overscore (ω)}d is the mean speed. This causes a problem because when the frequency of the speed ripple is the same as the mean drum speed, the inertial torque will be sinusoidal at the same frequency as the unbalance torque. Under these circumstances numerous combinations of Jd, mu and φu will result in a small estimator error and as such, the estimator may have difficulty in uniquely identifying the inertial and unbalance torque components from the overall drum torque (At larger values of speed ripple this becomes less of a problem because cos(ωd·t−φu) ≠cos({overscore (ω)}d·t−φu) and so the unbalance torque contains harmonic components, whereas the inertial torque is still sinusoidal. The parameter estimator can then more easily distinguish between the inertial and unbalance components of torque).
If the drum inertia is known before unbalance detection is carried out, the problem described above can be eliminated by compensating the drum torque for inertia. Assume for the moment that the drum inertia is known before unbalance detection is carried out. This being the case, the constant speed fixed axis model can be written as:
τd
where τd
Therefore, if the drum inertia is known, the inertial torque ripple can be subtracted from the overall ripple so that the input to the estimator is just simply the torque ripple due to gravity acting on the unbalance mass, plus the mean. Expanding the τd
τd
The above equation can be written in the following form suitable for Parameter Estimators such as the RLS estimator:
and P1-3 are the parameters, or coefficients, estimated by the PE.
It can be seen from the parameter matrix, θ, that further manipulation of P2 and P3 is required to extract the peak torque ripple and the torque ripple phase lag. This is done as follows:
If the acceleration is calculated from drum position, rather than speed, it is essential that the drum position be filtered. The reason for this is that the acceleration is obtained by differentiating the drum position twice. A small amount of noise on the drum position signal could result in very noisy acceleration data, which could in turn affect the accuracy of the estimated parameters. (However, if the speed is known the problem is not as severe because only one differentiation is required to calculate the acceleration; drum position in this case is determined by integrating the speed.)
Being a sawtooth waveform, the drum position signal contains high frequency components not associated with noise. Therefore, applying a low pass filter directly to the drum position signal will distort it significantly. To avoid these problems the sine and cosine of the signal is filtered and then the filtered drum position is reconstructed as follows:
In the embodiment described above, the drum inertia is calculated prior to the unbalance determination. One method of calculating the inertia is to ramp the drum speed from one predetermined speed to another, for example, 100 rpm to 200 rpm, and feed the resulting data into a parameter estimator with Jd as a tractable parameter. The torque and speed profiles will then look as shown in
In accelerating from ωd1 to ωd2, it is likely that the suspension will start to affect the mean torque and torque ripple. However, these suspension effects are not taken into account by the fixed axis model given by Equation 2. In calculating the inertia only, the change in torque ripple due to suspension effects is not of much importance. However, the change in mean torque as the speed is increased will affect the inertia calculation, but can be modeled as an additional viscous friction term, Bs, if it is assumed that the change in mean torque between ωd1 to ωd2 is a reasonably linear function of speed. Moreover, this additional viscous term can be lumped together with the actual viscous friction term.
Assuming that the mean torque varies linearly with speed between ωd1 and ωd2, the profile shown in
τd=Jdp(ωd)+Bjωd+τJ (13)
Although torque ripple is not that important as far as estimating inertia is concerned—at least when using a speed profile such as that in
τd=Jdp(ωd)+BJωd+τJ+A·cosθd+B·sin θd (14)
where A and B are constants. Re-arranging Equation 14 into the form y=xTθ for parameter estimation yields:
P1-5 are the parameters estimated by the parameter estimator. If the drum moves on its suspension significantly the torque ripple will change and so, therefore, will the A and B coefficients. As a result the estimator parameters P4 and P5 will also change in an attempt to track the changes in A and B. This is acceptable because in this particular embodiment only the inertia parameter P1 is required and this is not affected significantly by a gradual change in torque ripple when using a speed profile such as that in
In another exemplary embodiment of the invention, an alternative inertia calculation is employed that requires the drum be linearly accelerated from a first predetermined speed to a second higher speed, but does not require the use of a parameter estimator. As noted above, at around minimum plaster speed there is very little suspension movement and one can assume a fixed axis model. However, as the drum is accelerated beyond the minimum plaster speed, the effects of the suspension will become noticeable.
If the ripple component of the torque waveform is removed, the resulting trace looks something like that shown in
The calculation of τd
The acceleration torque is then given by:
The processor implementation for an exemplary embodiment is detailed below. Discretising the integration process yields the following equations for a, b and c:
The acceleration torque is therefore given by:
The inertia (kg.m2) is then given by:
In some applications it may be possible to accelerate the drum from ωd1 to a suitably higher speed, ωd2, without exciting the suspension significantly—particularly if ωd1 is set to be the minimum plaster speed. If this is the case then another embodiment of the invention is possible whereby the inertia and unbalance mass can be calculated simultaneously, thus eliminating the need to calculate inertia and mass independently as described in the above embodiment. There is also no need to perform inertia compensation. In this embodiment the generalized fixed axis model given by equation 3 can be assumed. Expanding the cos term in equation 3 yields:
τd=mu·g·ru·cos(α)·cos(φu)·cos(θd)+mu·g·ru·cos(α)·sin(φu)·sin(θd)+Jd·p(ωd)+Bd·ωd+τcsignum(ωd) (15)
It will be noticed that Equation 15 is of the same form as the inertia model given by Equation 14, with A=mu·g·ru·cos(α)·cos(φu) and B=mu·g·ru·cos(α)·sin(φu). Essentially, this embodiment uses the same procedure as calculating the inertia using a Parameter Estimator in the previous embodiment, except that in addition to calculating inertia, the unbalance mass can be extracted from the A and B parameters because the axis of the drum essentially remains fixed. The method for calculating mu from A and B above is the same method as that described in the second half of the previous embodiment.
As noted above, peak inertial torque ripple will be determined by Jd for a given speed ripple, and the gravitational torque ripple will be determined by mu assuming that ru can be calculated before hand.
The inertia of the laundry 110 can be calculated by subtracting the empty drum inertia from Jd:
J1=Jd−Jempty
The inertia of the laundry will be approximately equal to:
where m1 is the laundry mass.
The volume of the laundry is approximately equal to:
V1=1π(rd2−ru2)
where ‘1’ is the length of the drum. Assuming a figure for the density of modern fabrics, D1, the mass of the laundry can be written as:
m1=D1V1
ru is therefore given by the following equation
(It is anticipated that D1 will depend upon the type of laundry to be washed, e.g. synthetics, woolens, etc. will each have its own density).
In another embodiment of the present invention, the problem of having an “insufficiently excited” system is solved by introducing a “dither” signal to sufficiently excite all modes of the system. In the context of the system under investigation, this involves adding a signal to the speed command signal 202 such that the actual drum speed fluctuates around a mean value in a deterministic manner.
One requirement of the dither signal is that it should be of reasonably small amplitude so that the drum speed does not deviate too much from its mean value. Another requirement of the signal is that it should not excite any system resonances that could result in excessive vibration.
As described above, one of the reasons why the inherent speed ripple may not sufficiently excite all modes of the system is because its frequency is equal to that of the mean drum speed. For low values of speed ripple, this means that the unbalance torque and inertial torque are both approximately sinusoidal at the same frequency and the parameter estimator could therefore have difficulty in distinguishing between the two torque components. However, if the drum ripple contained frequency components of a different frequency to that of the mean drum speed, the unbalance and inertial components of torque would also differ in frequency and would therefore be more easily identifiable by the parameter estimator.
One of the benefits of using a dither signal to excite the system over the embodiments described previously is that unbalance mass and inertia can be calculated simultaneously at a substantially constant drum speed. Therefore, during unbalance detection the mean drum speed can be fixed at a speed to allow the use of the fixed axis model given by Equation 3, or if torque ripple due to viscous friction is assumed negligible, Equation 12 may be used. No particular form of dither signal is required, though in certain implementations, a 1 Hz sinusoidal dither signal with an amplitude of about 0.2 rad/s was used with satisfactory results.
As noted herein above, detecting the amount of a load unbalance allows the spin speed to be optimized to give an efficient wash and also minimizes the vibrations caused by centrifugal unbalance forces. To further improve efficiency, it is desirable to optimize the length of the spin cycle. Rather than simply operate a spin cycle for some predetermined time period, the spin cycle lasts until the water is removed from the laundry contained in the drum.
Accordingly, in accordance with further aspects of the present disclosure, the memory 108 stores program instructions executed by the processor 109 for detecting the proper time to end the spin cycle. In certain embodiments, the spin cycle is initiated upon completion of a drum unbalance detection process as described herein above. The unbalance detection may be run near the minimum plaster speed (about 80-100 rpm, for example). Further, to minimize vibrations due to drum unbalances, the spin speed may be selected based on the amount of unbalance. Once the amount of drum unbalance has been quantified, the drum 104 is then accelerated up to the desired spin speed.
For example, in one system considered, the drum 104 is run at a maximum spin speed of about 1154 rpm if the unbalance mass detected is less than 0.7 kg. If the unbalance mass is more than 2.5 kg, the spin speed is clamped at around 650 rpm. If the amount of unbalance is too high to run at a suitable spin speed the processor 109 is programmed to shake the drum 104 in order to redistribute the laundry.
As water is removed from the laundry during the spin process, the mass of the laundry reduces, resulting in a system that is non-linear, i.e., it is characterized by a set of changing parameters. Towards the end of a spin cycle, the amount of water loss from the laundry will be minimal and the system will approach a linear state where all its parameters are virtually constant. Therefore, monitoring the changes of the system parameters during the spin cycle using a suitable parameter estimator provides a good indication of when the spin cycle can be stopped.
A suitable parameter to monitor in determining the end of a spin cycle (when to command the motor 106 to stop spinning the drum 104) is the total drum inertia, including the laundry contained therein. Tests on towels indicate that the mass of the “wrung” laundry is about 60% of its saturated mass. The drum inertia should therefore decrease significantly throughout the spin cycle, regardless of how the laundry is distributed around the periphery of the drum 104. Other parameters will also change during the spin cycle, such as the unbalance mass and radius, and parameters relating to the suspension, but the amount by which they change will depend upon the amount of unbalance, which could be negligible. Such parameters are therefore not suitable for tracking purposes with regard to detecting the end of a spin cycle.
One method of tracking parameter changes is to carry out parameter estimation during the spin cycle in a piecemeal fashion, whereby estimates are taken over short time frames at regular points during the spin cycle. During the short interval when an estimate takes place the system is assumed to be linear, i.e. constant parameters, which is a reasonable assumption if estimation takes place over just several drum cycles. The parameter estimates are then compared with the previous estimates to get a measure of the amount by which the parameters are changing.
Another method of tracking parameter changes is to run a parameter estimation scheme continually “on-line” during the spin cycle until it has deemed to have finished by the automatic detection process. Any recursive estimator employed to estimate parameters must therefore be adaptive to the changing system parameters during the spin cycle. In an exemplary implementation, parameter estimation during the spin cycle is be done at the constant, mean spin speed using the moving axis model. To get the recursive estimator to track parameter changes, a forgetting factor technique is employed. More specifically, a technique known as “directional forgetting” is used to avoid problems with “estimator wind-up” (a known problem associated with the use of forgetting factors) in the exemplary embodiment disclosed.
As shown in
The unbalance mass 112 will almost certainly decrease, but it will also be affected by a slight change in laundry distribution as centrifugal forces compress the laundry 110 more during the spin cycle. This may also result in a small change in the unbalance mass angle φu. A decrease in the unbalance mass 112 will result in a decrease in the unbalance inertial and centrifugal forces that cause the drum 104 to move on its suspension. However, since the total mass of laundry will also decrease, the total mass supported by the suspension will be less and will therefore, generally speaking, require less force to accelerate it.
During the spin cycle, the drum 104 will be commanded to run at a constant speed, which is determined by the amount of unbalance in certain embodiments. However, due to a combination of the unbalance torque ripple and the less than ideal response of the speed controller, there typically will be some speed variation of the drum 104 about its commanded or mean value.
Because inertia (unbalance mass) has to be estimated at a substantially constant drum speed, there may be insufficient inertial torque information for the estimator to get a good estimate of inertia Jd, particularly if there is very little speed ripple. Moreover, if the speed ripple is mainly composed of a component of the same frequency as the mean drum speed, the estimator may have difficulty in distinguishing between the inertial and unbalance torque components. In other words, when a constant drum speed is commanded, the system may not be sufficiently excited enough to allow an accurate estimate of inertia to be obtained.
As detailed herein above for detecting the drum unbalance, to sufficiently excite the system, a dither signal may be employed so that the estimator will find it easier to distinguish the inertial torque component from the unbalance torque. Various types of dither signals may be adopted; for example, the dither signal may be a sinusoidally varying signal that is added to the speed reference signal. Its purpose is to cause the drum speed to vary by a small amount about the commanded spin speed in a manner that persistently excites the system.
In the disclosure above regarding estimating the amount of unbalance 112, the drum torque, τd, of a “fixed axis” drum 104 (drum is fixed in space—does not move on its suspension) inclined at α to the horizontal is given by Equation 3, which is repeated below.
τd=mu·g·ru·cos(α)·cos(θd−φu)+Jd·p(ωd)+Bd·ωd+τcsignum(ωd) (3)
As noted above,
At spin speeds encountered in a typical horizontal axis washing machine, it is likely that any unbalance centrifugal force will be sufficient to cause the drum 104 to move on its suspension. For an approximately constant drum speed, the force, and hence the acceleration, of the center of the drum 104 will vary approximately sinusoidally at the drum frequency in the x and y plane as shown in
The moving axis model for an approximately constant drum speed is therefore of the form:
τd=mu·g·ru·cos(α)·cos(θd−φu)+Jd·p(ωd)+Bd·ωd+τcsignum(ωd)+B·cos(2·θd−φs)+τsusp (16)
Equation 16 is essentially Equation 3 with an additional second harmonic term and an additional constant term. As the laundry dries out in the spin cycle, the reduction in laundry mass will cause the overall suspended mass, i.e., the drum mass, to reduce and also the centrifugal force. Therefore, although the suspended mass will become lighter, making it easier to accelerate, the force causing the acceleration will also reduce.
The drum speed can therefore be split into two components as follows:
ωd={overscore (ω)}d+ωd
Using Equation 17, the viscous friction term in Equation 16 can be expanded to give:
τd=mu·g·ru·cos(α)·cos(θd−φu)+Jd·p(ωd)+Bd·ωd
where τd
In general the viscous friction coefficient is small and the speed controller is designed to keep the speed ripple down to a minimum—ideally no more than a few percent of mean drum speed. The torque ripple due to viscous friction can therefore be neglected in most cases to give
τd=mu·g·ru·cos(α)·cos(θd−φu)+Jd·p(ωd)+B·cos(2·θd−φs)+τd
For the purposes of parameter estimation Equation 19 is expanded to give:
τd=mu·g·ru·cos(α)·cos(φu)·cos(θd)+mu·g·ru·cos(α)·cos(φu)·sin(θd)+
Jd·p(ωd)+B·cos(φs)·cos(2·θd)+B·sin(φu)·sin(2·θd)+τd
Equation 20 can be written in the following form suitable for Parameter Estimators such as the RLS estimator:
and P1-6 are the parameters, or coefficients, estimated by the PE.
Assuming that the drum torque is given by Equation 20, the system parameters are given by Equation 23. During the spin cycle, these parameters will change as a result of water being removed from the laundry by centrifugal force. In order for a parameter estimator such as a recursive least squares estimator to track these parameter changes, the estimator must be adaptive to parameter variation. There are several techniques that enable a parameter estimator to become adaptive. A suitable technique is the use of a forgetting factor.
The forgetting factor gives the parameter estimator the ability to adapt to parameter changes by reducing the emphasis the estimator places on past information. The forgetting factor, λ, can be interpreted as a forgetting time constant with a data memory roughly equal to (1−λ)−1.
In general, the ability of the parameter estimator to accurately track parameter changes increases as the amount of forgetting is increased, i.e., as λ is reduced. However, this is not to say that one can arbitrarily set λ to a very low value to “be on the safe side” and assume that the estimator will have good tracking capability regardless of how the parameters change. It is therefore important to understand the factors influencing the choice of λ if satisfactory results are to be obtained.
If the estimator is only working on the last “data points remembered” (DPR) number of data points, then the data within that span of data points must not change too significantly if an accurate estimate is to be obtained. This means that less forgetting (high λ) is required for systems where the rates of change of system parameters are small. Conversely, if the rates of change of system parameters are high, then it is likely that accurate estimates are only obtainable with lower values of λ. The choice of forgetting factor is therefore determined by a combination of the sampling rate and the rate of change of parameters.
Given that the rate of change of parameters with time is the significant factor in the determination of λ, it is more convenient, and meaningful, to think of λ in terms of “time remembered” (TR), rather than the number of DPR. TR is simply DPR multiplied by the sampling time:
TR=τs·(1−λ)−1 (24)
If the TR is too small, then not enough information, with respect to time, is available for the estimator to get an accurate estimate.
In general, it can be stated that the end of the spin cycle occurs when the drum inertia has settled down to a reasonably constant value. This suggests that in using the inertia parameter to detect the end of the spin cycle, it is appropriate to monitor the rate of change of the drum inertia. In this way, the end of a spin cycle can be defined as being the point when the rate of change of inertia approaches a predetermined value close to zero.
The differentiation process is a good amplifier of noise, so when the inertia estimate 210 is differentiated, the result may contain so much noise that the rate of change of the general trend will not be identifiable. Filtering may be employed to overcome this problem. Because the rate of change of inertia is low when taken over the duration of a spin cycle, it is possible to heavily filter the signal without causing any significant phase distortion of the estimate's general trend. Any higher frequency components on the inertia estimate can be reduced significantly by employing filters with low cut-off frequencies and relatively high orders.
Alternatively a gain, G, can be incorporated in the differentiation process such that the rate of change of inertia, pJd, is approximated as follows:
At low frequencies the approximation given by Equation 25 is accurate, whereas at high frequencies the gain is limited to G. G is typically in the range 3 to 20. Equation 25 can be discretised for digital implementation by any of the standard methods such as the Tustin Transformation.
Another method of addressing the issue of a noisy estimate derivative is to effectively re-sample the inertia estimate at a lower rate before differentiation takes place. Because the inertia changes gradually over the course of a spin cycle, it is not necessary to calculate its derivative at the relatively high rates of sampling that are required by the estimator to get an accurate estimate of inertia Jd. For example, if the inertia is estimated at a rate of 1 ms, it may only be necessary to consider every 1000th inertia estimate when calculating p(Jd), which is effectively the same as resampling the data at 1 s. Moreover, this method could be used in conjunction with filtering.
In other embodiments, the spin cycle end is detected by using the raw inertia estimate only. Prior to the spin cycle being detected, an estimate of drum inertia and an approximate calculation of the unbalance radius is obtained as part of the unbalance detection process at around 100 rpm. Using a representative figure for the reduction in mass of the laundry, which can be obtained by experiment and stored in the memory 108, it is possible to calculate what the inertia should be at the end of the spin cycle. This value is then used as a “target” value for determining the completion of the spin cycle. When the inertia drops below the projected value and stays within specified limits for a specified period of time, the spin cycle is deemed to have completed.
The particular embodiments disclosed above are illustrative only, as the invention may be modified and practiced in different but equivalent manners apparent to those skilled in the art having the benefit of the teachings herein. Furthermore, no limitations are intended to the details of construction or design herein shown, other than as described in the claims below. It is therefore evident that the particular embodiments disclosed above may be altered or modified and all such variations are considered within the scope and spirit of the invention. Accordingly, the protection sought herein is as set forth in the claims below.
This application is a continuation-in-part of U.S. patent application Ser. No. 10/655,626, filed on Sep. 4, 2003, which claims the benefit of U.S. Provisional Patent Application No. 60/466,106, filed on Apr. 28, 2003, both of which are incorporated by reference herein.
Number | Date | Country | |
---|---|---|---|
60466106 | Apr 2003 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 10655626 | Sep 2003 | US |
Child | 11127634 | May 2005 | US |