Semiactive Vibration Control for Horizontal Axis Washing Machine

A semiactive vibration control method is developed to cope with the dynamic stability problem of a horizontal axis washing machine.Thismethod is based on adjusting themaximum force values produced by the semiactive suspension elements considering a washingmachine’s vibration data (three axis angular position and three axis angular acceleration values in time). Before actuation signals are received by the stepmotors of the friction dampers, vibration data are evaluated, and then, the stepmotors start to narrow or expand the radius of bracelets located on the dampers. This changes the damping properties of the damper in the suspension system, and thus, the semiactive suspension system absorbs unwanted vibrations and contributes to the dynamic stability of the washing machine. To evaluate the vibration data, the angular position and angular acceleration values in three axes are defined in a function, and the maximum forces produced by semiactive suspension elements are calculated according to the gradient of this function. The relation between the dynamic stability and the walking stability is also investigated. A motion (gyroscope and accelerometer) sensor is installed on the top-front panel of the washingmachine because a mathematical model of a horizontal axis washingmachine suggests that the walking behavior starts around this location under some assumptions, and therefore, calculating the vibrations occurring there is crucial. Semiactive damping elements are located under the left and right sides of the tub. The proposedmethod is tested during the spinning cycle of washingmachine operation, increasing gradually from 200 rpm to 900 rpm, which produces the most challenging vibration patterns for dynamic stability. Moreover, the sound power levels produced by the washing machine are measured to evaluate the noise performance of the washing machine while the semiactive suspension system is controlled. The effectiveness of the proposed control method is shown through experimental results.


Introduction
During the operation of a washing machine, the main dynamic problem is an unbalanced load of unpredictable position and magnitude.As the drum rotates, clothes separate from each other, causing stability problems.A washing machine is dynamically stable if unbalanced load only produces the forces and torques in such a way as to create a small oscillation about the equilibrium point of washing machine.However, even for these small oscillation values, a washing machine may show walking instability behaviour because of washing machine's rigidness; thus, it can be said that just because a washing machine is dynamically stable does not mean that it has walking stability as well.Conrad and Soedel investigated the main reasons for walking instability in horizontal and vertical axis washing machines using rudimentary dynamic models [1].However, the washing machine model that they used did not contain elastic components, and therefore, the effects of passive suspension elements were neglected.Conrad later studied a more advanced dynamic mathematical model of horizontal and vertical axis washing machines for selected constraints and showed an approach to implementing dynamic elements such as a spring and a damper into the static model of such washing machines [2].These studies proved the importance of suspension elements.This led to the question of what the effect of semiactive suspension systems that have adjustable spring and damping coefficients would be if spring and damping elements having constant coefficients can contribute to the dynamic stability of washing machines.Researchers thus began focusing on solving the dynamic stability problem of washing machines using semiactive suspension systems.Semiactive suspension elements can successfully deal with the dynamic stability problem by manipulating the centrifugal forces acting on the cabinet by storing and dissipating them so as to safely transfer the centrifugal forces produced by clothes from the drum to the ground.Changing the magnitude and direction of transferred forces contributes to the stable dynamic behavior of washing machines [3].However, the question is how much of the force should be stored and how much should be dissipated.In this regard, the selection of the type of semiactive suspension is crucial.Magnetorheological (MR) dampers appear to be the best choice because their damping coefficients can be changed in wide variety; furthermore, they show fast response and are easy to control through an electronic interface [4].However, using them in the home appliance industry may be uneconomical owing to production constraints.Nevertheless, several studies have focused on using them for reducing vibrations in a semiactive suspension system, that is, for dynamic stability, and for noise reduction.Spelta et al. [5] improved an adaptive control method to deal with unwanted vibrations and noise using a semiactive suspension system (MR-controllable friction damper); this algorithm is based on changing the damping parameters in response to the vibrations.
In our study, the working structure of the semiactive suspension system is based on adjusting the maximum force produced under load.Step motors are installed on suspension elements located under the tub.For each actuation signal, the step motors narrow or expand bracelets on suspension elements made of polyurethane; therefore, the elastic properties of polyurethane start changing, and the stiffness and damping coefficients change accordingly.However, owing to the structure of the suspension system, the only known parameter is the maximum force that can be produced by suspension elements.In the manufacturing process of the suspension system, the maximum reaction forces under load for each elastic level of the suspension system are measured.No motion sensor is installed on the suspension system, and the actual stiffness and actual damping coefficients are unknown; therefore, the actual reaction force of the suspension system cannot be determined.Using this type of suspension system provides an economic advantage compared with MR semiactive suspensions; on the other hand, the actual stiffness and actual damping values of suspension elements are unknown during operation.
Some other studies have investigated vibration reduction using classic suspension systems for horizontal washing machines.Türkay et al. investigated parametric optimization methods for reducing the maximum orbit displacement.In their method, grid and sequential quadratic programming optimization methods were used [6].In another study, they used the Newton-Euler method to derive a nonlinear time-variant rigid body dynamic model of a classic suspension system of a horizontal axis washing machine for experimental assessment [7].Boyraz and Gündüz investigated the optimization of the vibration characteristics using a generic algorithm [8].Öztürk and Erol optimized the dynamic behavior of a horizontal axis washing machine using Adams (MSC), a commercial multibody simulation software [9].
Papadopoulas and Papadimitriou provided a mathematical model of a horizontal axis washing machine assuming that the machine is a rigid structure [10].In this model, the maximum angular velocity that can be reached (walking stability threshold) is determined according to other parameters such as the weight of clothes and friction between the ground and the feet.In our study, this mathematical model was first adapted to a washing machine with a semiactive suspension system; therefore, the assumption of rigidness is no longer completely valid, and the stiffness and damping coefficients can be considered in the equation for the maximum angular velocity that can be reached.This situation provides an opportunity for manipulation and proves that if the forces produced by semiactive suspension elements are changed appropriately, the maximum angular velocity can reach higher levels.However, the lack of modeling of some elastic components of the washing machine (cabinet, tub, etc.) and the assumption that the friction coefficient between the ground and the feet of the washing machine is constant may be weaknesses of the mathematical model.Sound performance of the washing machine must also be considered while the vibration performance is being improved [11].Therefore, sound power outputs of the washing machine for the proposed control method are measured during operation and compared with those while the washing machine is working with an uncontrolled suspension system.The similarities and differences of the outputs at different stages of the spin cycle are discussed.
Our study aims to solve the dynamic stability problem using low-cost equipment (suspension system that works based on adjusting the maximum force produced, motion sensor, and data acquisition card), and the success of the control method for walking instability is evaluated according to two criteria: Papadopoulas and Papadimitriou's walking stability threshold adapted to a horizontal axis washing machine working with a semiactive suspension system and reduction of the arithmetic mean of the position and acceleration values in three axes.

Theory
An unbalanced load owing to rotating clothes in the drum has unpredictable position and magnitude values.To simplify the calculations, some assumptions have to be made.The first assumption is that the radius of clothes in the drum,   , is constant.  and   are the components of the and -axes of the unbalanced force produced by rotating clothes in the drum: where   is the mass of the unbalanced load;  is the gravitational acceleration;  is the angle between the load and the -axis; and  is the angular velocity of the drum around the -axis.The second assumption pertains to the rigidness of the washing machine, with the forces produced by rotating clothes (  ,   ) and friction forces on feet (  ,   ,   , and   ) shown in Figure 1.
Points , , , and , respectively, indicate the front left, front right, back left, and back right feet of the washing machine.
If it is assumed that the washing machine has no suspension, the following moment of force equations are used to derive walking stability threshold.In the equations, moment of force parameters are calculated with respect to -axis shown in Figure 1: where  is the distance between the front and the back feet;  1 is the distance between the back feet and the center of gravity of the washing machine;  2 is the distance between the back feet and the center of gravity of the load; and   is the mass of the washing machine.When  2 is greater than half of , the sum of   and   on the -axis becomes greater than that of   and   on the -axis: When  1 is less than half of , the sum of   and   on the -axis becomes greater than that of   and   on the -axis: Considering these conditions, walking instability first occurs around the front feet because they are subjected to more force than the back feet.Thus, in our study, a motion sensor is installed on the middle of the front-top panel of the washing machine.
is the friction between the ground and the feet: Using ( 3) and ( 5), this equation can be rewritten as Substituting ( 1) and ( 2) into (10) gives To make the equation independent of , the derivative of (12) with respect to  has to be equal to 0: Then,  becomes All equations shown above hold under the assumption of the rigidness of the washing machine, and they change when the washing machine contains a suspension system.Because of the structure (spring and damping elements) of the semiactive suspension system and the passive suspension system,   and   given by ( 1) and ( 2) change, and their new forms are, respectively, given by where  1 is the stiffness coefficient of the left semiactive suspension element;  2 is the stiffness coefficient of the right semiactive suspension element;  3 is the stiffness coefficient of the left passive spring;  4 is the stiffness coefficient of the right passive spring;  1 is the damping coefficient of the left semiactive suspension element;  2 is the damping coefficient of the right semiactive suspension element;  1 is the angle between the ground and the left semiactive suspension element;  2 is the angle between the ground and the right semiactive suspension element; p1 is the displacement of the left semiactive suspension element; p2 is the displacement of the right semiactive suspension element;  3 is the angle between the top panel and the left passive spring;  4 is the angle between the top panel and the right passive spring; p3 is the displacement of the left passive spring; and p4 is the displacement of the right passive spring.Substituting these new forms into (3) and (5) gives Thus, the walking stability threshold changes as follows: As applied previously, to make the equation independent of , the derivative of (18) with respect to  has to be equal to 0: Changing  1 ,  2 ,  1 , and  2 changes the walking stability threshold of the washing machine; this change can be achieved using a semiactive suspension system.However, owing to the structural constraints of the semiactive suspension system used in this study, the only known knowledge about the suspension elements is the maximum force value that can be produced.Therefore, it can be said that  1 p1 +  1 ( p1 /) and  2 p2 +  2 ( p2 /) are the only known values in (19), and it is assumed that these values are always equal to the maximum forces produced by suspension elements.
To change the maximum force values produced by the semiactive suspension elements for reducing vibrations, the angular position and angular acceleration values produced by the washing machine have to be interpreted.This can be done using the function  that contains three-axis angular position (  ,   , and   ) and three-axis angular acceleration (  ,   , and   ) values, with the angular axis acceleration values being squared owing to the importance of acceleration in the vibration behavior at high frequencies: According to the gradient of function , the maximum force values produced by the suspension elements are chosen.
When the gradient of function  becomes higher, it indicates that the value of function  will increase dramatically.Therefore, the vibrations will reach higher levels, to prevent which the suspension elements must apply lower force values to the drum.Otherwise, the washing machine will show unstable vibration behavior: The intervals of function  for different maximum force values of suspension elements are shown in Table 1.
Each interval corresponds to a different vibration level of the washing machine, and the maximum forces applied by the semiactive suspension elements change to reduce the vibrations so as to deal with the dynamic stability problem that may occur during the spinning cycle.While the gradient of  increases, the maximum forces produced by the suspension elements decrease and vice versa.This process directly changes the effects of the suspension elements in the equation for the walking stability condition (19).Figure 2 shows the flowchart of the control algorithm used, for a better understanding.

Experimental Studies
Experiments were performed for the spinning cycle of a washing machine, increasing gradually from 200 rpm to 900 rpm, which produces the most challenging vibration patterns from the viewpoint of dynamic stability.The maximum weight capacity of the washing machine is 7 kg.The equation shown below is used as the first evaluation criterion (Papadopoulas and Papadimitriou's walking stability threshold adapted to a horizontal axis washing machine working with semiactive suspension system) for the success of the controlled semiactive suspension system. is the magnitude of  at each sampling time and  is the sum of samples: The higher the value of function   , the more stable the washing machine. values have to be calculated for each sampling time to obtain   values.As mentioned before, for each sampling time, vibration data are recorded, gradient of   function is calculated, and the maximum forces produced by suspension elements are changed to calculate .
The second criterion is the arithmetic mean of the position and acceleration values in three axes.This criterion can also be used for investigating the changing rates of each position and acceleration parameters.Function  V is defined as follows: where  is the magnitude of each signal (position or acceleration) at each sampling time and  is the sum of samples.The lower the value of function  V is, the less oscillation it shows from its initial value 0, indicating the success of the controller.The process of changing the maximum forces of suspension elements can only be achieved at the same time owing to the structure of the motor drive module.A positive edge PWM signal is sent from the DAQ card to the motor drive module through a circuit including a 3 k resistor and an NPN transistor, as shown in Figure 3.The connection between the resistor and the transistor is serial in nature.
A single PWM signal makes the suspension elements take only one step on the gradient line, and the number of PWM signals that have to be produced by the DAQ card is selected according to the gradient of function .For instance, if the gradient of function  is 3 at time  and 6 at time  +   , three PWM signals have to be produced to change the maximum forces of the suspension elements.  is selected as 300 ms for data acquisition and as 600 ms for control signal that drives step motors.
The MPU 6050 motion sensor is used in this study.It includes a 3-axis accelerometer and 3-axis gyroscope, and it has 16-bit analog-digital converters for digitizing the accelerometer and gyroscope outputs.The userprogrammable gyroscope scale is adjusted to ±2000 ∘ /s and the user-programmable accelerometer scale to ±16 g.Even though the DAQ card provides 5 V, the sensor has a voltage regulator for its working voltage of 3.3 V. Communication with the data acquisition card is performed using the I 2 C protocol at 400 Hz from the SCL and SDA channels.The communication software is programmed in MATLAB using the C programming language.As mentioned in (8), walking instability first occurs around the front feet because they are subject to more force than the back feet.Thus, the sensor is installed on the middle of the front-top panel of the washing machine, as shown in Figure 4.The Arduino Mega 2560 DAQ card is used in this study.It includes an ATmega microcontroller, and communication between the PC and the DAQ card is performed using the serial communication protocol.A positive edge PWM signal is sent from the digital pin of the DAQ card after receiving sensor data.The communication software is programmed in MATLAB using the C programming language.
The semiactive suspension system used in this study is shown in Figure 5. Inside the suspension elements, there is a cylindrical mass made of polyurethane.The data cable carries a signal from the motor drive module that drives the step motor.Depending on the type of signal, the bracelet starts to narrow or expand, causing a change in the elastic properties of the suspension element.
The experiments are performed in a semianechoic room, as shown in Figure 6.The sound power levels for uncontrolled and controlled suspension systems are measured according to ISO 3745.
The physical parameters of the washing machine are listed in Table 2.The trajectory of the passive spring elements is assumed to be 0.1 sin(5) and  1 ,  2 ,  3 , and  4 are assumed to be constant and equal to 55 ∘ .The friction coefficient between Step motor Data cable Bracelet Suspension element  the feet and the ground is assumed to be constant and equal to 5.

Results and Discussion
To show the success of the control method in damping vibrations, the measured position and acceleration values when the semiactive suspension system is uncontrolled (maximum force of left and right suspensions is 120 and 112, resp.) and controlled are shown in Figures 7-18.
As can be seen in Figure 7, for controlled system, variation between minimum and maximum values of   is considerably low.Minimum value of   is −6 ∘ and maximum value of   is +6 ∘ .On the other hand, for uncontrolled system shown in Figure 8, the value of   oscillates between −60 ∘ and +60 ∘ , which means that 90% reduction of minimum and maximum values of vibration for -axis is confirmed.Moreover, in Figure 7, it can be seen that the value of   settles around −1 ∘ ; the reason of this situation is the noise that corrupts sensor data.
In Figure 9, the value of   for controlled system settles around 0 ∘ with a low value of steady-state error, which is described with  V in Table 3; variation between minimum and maximum values of   is considerably low.Minimum value of   is −6 ∘ and maximum value of   is +6 ∘ .For uncontrolled system as shown in Figure 10, the value of   oscillates between −60 ∘ and +40 ∘ , which means that minimum value  of vibration in -axis is reduced by 90%, whereas maximum value of vibration in -axis is reduced by 85%.
In Figure 11, the value of   for controlled system settles around 0 ∘ ; peak values of   are high with respect to steadystate error of   , whereas the value of   settles around +2 ∘ for uncontrolled system as shown in Figure 12.
As can be seen from Figures 7-12, even though there are some peak values, vibrations in -, -, and -axes are considerably reduced for controlled system, especially for   and   parameters.Comparison for angular acceleration parameters is shown in Figures 13-18.In Figure 13, the value of   for controlled system settles around −0.4 ∘ ; peak values of   are high with respect to steady-state error of   and increase through the end of washing operation.For uncontrolled system shown in Figure 14, it can be seen that the value of   oscillates between −1 ∘ and +0.7 ∘ .
In Figure 15, the value of   for controlled system settles around 0 ∘ ; peak values of   are the highest ones among other angular acceleration parameters shown in Figures 13  and 17.However, this situation does not corrupt improvement percentage of  V as can be seen in Table 4.For uncontrolled system shown in Figure 16, it can be seen that the value of   oscillates between −1 ∘ and +0.7 ∘ and there are minimum and maximum peak values around −4 ∘ and +4 ∘ , respectively.
In Figure 17, the value of   for controlled system settles around 0 ∘ and varies between −2 ∘ and +2.5 ∘ .For uncontrolled system shown in Figure 18,   settles around 0 ∘ with high frequency change.As can be seen from Figures 13, 15, and 17, angular acceleration parameters for controlled system have peaks, just as angular position parameters shown in Figures 7, 9, and 11; however, in this case, the peak values reach beyond the boundaries of uncontrolled system's angular acceleration outputs.The main reason of the existence of these peak values is that, due to hardware constraints, the sampling time   is adjusted to 300 ms for data acquisition while it is adjusted to 600 ms for control signal that drives step motors.When the gradient of function  drastically changes (and the system becomes prone to be unstable), the control signal cannot quickly affect the system because of time delay between data acquisition and control signal.

Shock and Vibration
The  V values of each position and acceleration data for a controlled and an uncontrolled semiactive suspension system are listed in Table 3.
Table 4 suggests that the controller mostly succeeds in terms of the position   and acceleration   .

Shock and Vibration
The change in  over time is shown in Figure 19;   is calculated as 1080.2rpm.
Sound power levels of the washing machine for controlled and uncontrolled cases have also been measured according to ISO 3745 in a semianechoic room as shown in Figure 6 and the measurement results are presented in Figure 20.
Figure 20 shows that at the low spin speeds the sound power level is more stable when the semiactive suspension system is controlled; in comparison, the uncontrolled system shows more oscillatory behavior for the first spinning cycles.At the highest final spin speed, 900 rpm, Figure 20 shows that there is no major difference between the measured sound power levels.
In Figure 20, it can be seen that the value of sound power changes between 50 dBA and 63 dBA when the system is uncontrolled and the value of sound power changes between 54 dBA and 62 dBA when the system is controlled for the first spinning cycles.

Conclusions
This study proves that low-cost equipment including an MPU 6050 sensor, an Arduino Mega 2560 used as data acquisition card, and a semiactive suspension system working based on maximum force adjustment can be used to solve the walking instability problem of washing machines.When the system is controlled, angular position values in -, -, and -axes settle around 0 ∘ with low steady-state errors, which means that they are considerably reduced.Even though angular acceleration values have peaks during the operation, they have not any major effects on walking and dynamic stability.
Finally, the measured sound power levels prove that a washing machine with a controlled semiactive suspension system shows more stable sound performance than an uncontrolled one for the spinning cycle.Furthermore, even though there is an enormous difference between the controlled and the uncontrolled semiactive suspension systems from the viewpoint of vibration damping, the final sound power levels of the washing machine for both conditions are almost the same (72 dBA and 73 dBA), suggesting that there is no linear relation between vibrations causing walking instability and sound power levels produced by the washing machine in this case.
Creating J function Calculation of gradient of J function for each sample Changing max force of semiactive suspension element according to gradient of J function Measurement of acceleration and position values in three axes

Table 1 :
Different maximum force values of suspension elements.

Table 4 :
Improvement percentage of  V values for position and acceleration data compared with their uncontrolled values.