Skip to main content

Optimal design and performance evaluation of a flow-mode MR damper for front-loaded washing machines

Abstract

It is well known that the vibration of washing machines is a challenging issue to be considered. This research work focuses on the optimal design of a flow-mode magneto-rheological (MR) damper that can replace the conventional passive damper for a washing machine. Firstly, rigid mode vibration of the washing machine due to an unbalanced mass is analyzed and an optimal positioning of the suppression system for the washing machine is considered. An MR damper configuration for the washer is then proposed considering available space for the system. The damping force of the MR damper is derived based on the Bingham rheological behavior of the MR fluid. An optimal design problem for the proposed MR brake is then constructed. The optimization is to minimize the damping coefficient of the MR damper while the maximum value of the damping force is kept being greater than a required value. An optimization procedure based on finite element analysis integrated with an optimization tool is employed to obtain optimal geometric dimensions of the MR damper. The optimal solution of the MR damper is then presented with remarkable discussions on its performance. The results are then validated by experimental works. Finally, conclusions on the research work are given and future works for development of the research is figured out.

Background

It is well known that the vibration of washing machines is a challenging issue to be considered. The vibration of the washing machine is mainly due to the unbalanced mass of clothes distributed in the washing drum. This occurs most frequently in the spin-drying stage, because the drum spins at a relatively high speed causing the clothes to be pressed against the inner wall of the spin drum, and these can become a large unbalanced mass until the end of the stage. Particularly, in a front-loaded washing machine (drum-type washing machine), the unbalanced mass of clothes easily occurs and very severe due to the effect of gravity. The vibration of the washing machine is transferred to the floor causing noises, unpleasant feeling for humans, and failure of the machine.

There are many researches on vibration control of washing machines which can be classified into two main approaches. The first approach is based on the control of the tub balance to eliminate the source of vibration [1, 2]. In this approach, one type of dynamic balancer is used to self-balance the tub dynamics. A typical dynamic balancer is the hydraulic balancer containing salt water, which is attached to the upper rim of the basket. The liquid in the balancer moves to the opposite side of unbalance automatically due to the inherent nature of fluids when the rotational speed is higher than the critical speed of the spinning drum [1]. Another dynamic balancer that counteracts vibrations is to use two balancing masses. In this method, two balancing masses move along the rim of the basket. The rotation plane of the balancing masses can be easily chosen to be wherever judged suitable, always targeting at the reduction of the induced moments [2]. It is proved that the vibration of the washing machine can be significantly reduced by using a dynamic balancer. However, the complicated structure, high cost of manufacturing, and maintenance are a big obstacle for the wide application of this approach. In the second approach, the vibration of the washing machine is suppressed based on damping control of a suspension system [3]. It is noted that during the spinning process, the washing machine usually experiences the first resonance at quite low frequency, around 100 to 200 rpm. This results from the resonance of the washing drum due to the unbalanced mass. When the rotating speed exceeds 1,000 rpm, the side and rear panels of the frame may experience resonances which cause noises and vibration transferred to the floor. If a passive damper is used to reduce the vibration of the drum at low frequency, it will cause the vibration of the washing machine at high frequencies more severe. The reason is that more excitation force from the drum is transferred to the frame via the passive damper. Therefore, in order to effectively reduce the vibration of the washing machine at low frequency while the vibration of the machine at high frequencies is insignificantly affected, a semi-active suspension system such as a magneto-rheological (MR) damper should be employed.

Although there have been several researches on the design and application of MR dampers to control the vibration of washing machines [4–6], the optimal design of such MR dampers was not considered. The main objective of this study is to achieve the optimal design of the semi-active suspension system for washing machines employing MR dampers. Firstly, rigid mode vibration of the washing machine due to an unbalanced mass is analyzed and an optimal positioning of the suppression system for the washing machine is considered. An MR damper configuration for the washer is then proposed. The damping force of the MR damper is then derived based on the Bingham model of the MR fluid. An optimal design problem for the proposed MR brake is then constructed considering the damping coefficient and required damping force of the MR damper. An optimization procedure based on finite element analysis integrated with an optimization tool is employed to obtain optimal geometric dimensions of the MR damper. The optimal solution of the MR damper is then presented with remarkable discussions. The results are then validated by experimental works.

Methods

Vibration control of washing machine using MR damper

The washing machine object of this work is a prototype based on the LG F1402FDS washer manufactured by LG Electronics (Seoul, South Korea). A three-dimensional (3D) schematic diagram of the washer is shown in Figure 1. It is characterized by a suspended tub (basket) to store the water for washing linked to the cabinet by two springs and two dampers. The rotor is directly connected with the drum which rotates against the tub while the stator is fixed on the back of the tub. When the drum is rotating, the unbalanced mass due to the eccentricity of laundry causes the vibration of the tub assembly. The vibration of the tub assembly is then transmitted to the cabinet and the bottom through the springs and the dampers. In Figure 2, a two-dimensional (2D) simplified schematic of the machine is depicted. From the figure, the following governing equation of the washing machine can be derived:

m ü + c u ˙ sin 2 φ + β 2 + sin 2 φ − β 1 + ku sin 2 φ + α 1 + sin 2 φ − α 2 = F u t
(1)

where m is the mass of the suspended tub assembly including the drum, laundry, shaft, counter weight, rotor, and stator. For the prototype washing machine, m is roughly estimated about 40 kg. c is the damping coefficient of each damper, and k is the stiffness of each spring which is assumed to be 8 kN/m in this study. φ is the angle of an arbitrary direction (u direction) in which the vibration is considered. F u is the excitation force due to unbalanced mass in the u direction, F u  = F0cosωt = m u ω2R u cosωt, in which m u and R u are the mass and radius from the rotation axis of the unbalanced mass. From Equation 1, the damped frequency of the suspended tub assembly is calculated by

ω d = ω n 1 − ξ 2
(2)

where ω n = k sin 2 φ + α 1 + sin 2 φ − α 2 m and ξ = c sin 2 φ + β 2 + sin 2 φ − β 1 2 mk sin 2 φ + α 1 + sin 2 φ − α 2 .

Figure 1
figure 1

3D schematic of the prototype washing machine.

Figure 2
figure 2

2D simplified schematic of the prototype washing machine.

It is seen that the damped frequency and natural frequency of the tub assembly in the u direction are a function of φ. Therefore, in general, in a different direction of vibration, the tub assembly exhibits different resonant frequency. This causes the vibration to become more severe and hard to control. In the design of the suspension system for the tub assembly, the frequency range of the resonance in all directions should be as small as possible. From the above, it is easy to show that, by choosing α1 + α2 = 90° and β2 + β2 = 90°, Equation 1 can be simplified to yield

m ü + c u ˙ + ku = F u
(3)

In this case, the damped frequency and natural frequency of the tub assembly do not depend on the direction of the vibration. We have

ω d = ω n 1 − ξ 2
(4)

where ω n = k m and ξ = c 2 mk .

An inherent drawback of the conventional damper is its high transmissibility of vibration at high excitation frequency. In order to solve this issue, semi-active suspension systems such as ER and MR dampers are potential candidates. In this study, two MR dampers are employed to control the vibration of the tub assembly.

Figure 3 shows the schematic configuration of a flow-mode MR damper proposed for the prototype washing machine. From the figure, it is observed that an MR valve structure is incorporated in the MR damper. The outer and inner pistons are combined to form the MR valve structure which divides the MR damper into two chambers: the upper and lower chambers. These chambers are fully filled with the MR fluid. As the piston moves, the MR fluid flows from one chamber to the other through the annular duct (orifice). The floating piston incorporated with a gas chamber functions as an accumulator to accommodate the piston shaft volume as it enters and leaves the fluid chamber.

Figure 3
figure 3

Schematic configuration of the MR damper.

By neglecting the frictional force and assuming quasi-static behavior of the damper, the damping force can be calculated by [7]

F d = P a A s + C vis x Ë™ p + F MR sgn x Ë™ p
(5)

where P a , cvis, and FMR, respectively, are the pressure in the gas chamber, the viscous coefficient, and the yield stress force of the MR damper which are determined as follows:

P a = P 0 V 0 V 0 + A s x p γ
(6)
C vis = 12 ηL π R d t d 3 A p − A s 2
(7)
F MR = 2 A p − A s 2.85 L p t d τ y
(8)

In the above, τ y is the induced yield stress of the MR fluid which is an unknown and can be estimated from magnetic analysis of the damper and η is the post-yield viscosity of the MR fluid which is assumed to be field independent. R d is the average radius of the annular duct given by R d  = R − t h  − 0.5t d . L and R are the overall length and outside radius of the MR valve, respectively. t h is the valve housing thickness, t d is the annular duct gap, and L p is the magnetic pole length.

In this work, the commercial MR fluid (MRF132-DG) made by Lord Corporation (Cary, NC, USA) is used. The post-yield viscosity of the MR fluid is assumed to be independent on the applied magnetic field, η = 0.1 Pas. The induced yield stress of the MR fluid as a function of the applied magnetic field intensity (H) is shown in Figure 4. By applying the least squares curve fitting method, the yield stress of the MR fluid can be approximately expressed by

Ï„ y = C 0 + C 1 H + C 2 H 2 + C 3 H 3
(9)
Figure 4
figure 4

Yield stress of MR fluid as a function of magnetic field intensity.

In Equation 9, the unit of the yield stress is kilopascal while that of the magnetic field intensity is kA/m. The coefficients C0, C1, C2, and C3 are respectively identified as 0.044, 0.4769, −0.0016, and 1.8007E-6. In order to estimate the induced yield stress using Equation 9, first the magnetic field intensity across the active MRF duct must be calculated. In this study, the commercial FEM software, ANSYS, is used to analyze the magnetic problem of the proposed MR damper.

Optimal design of the MR damper

In this study, the optimal design of the proposed MR damper is considered based on the quasi-static model of the MR damper and dynamic equation of the tub assembly developed in the ‘Vibration control of washing machine using MR damper’ section. From Figure 2 and Equation 3, the force transmissibility of the tub assembly to the cabinet can be obtained as follows:

TR = 1 + 2 ξr 2 1 − r 2 2 + 2 ξr 2
(10)

where r is the frequency ratio, r = ω/ω n . The dependence of the force transmissibility on excitation frequency is presented in Figure 5. As shown from the figure, at low damping, the resonant transmissibility is relatively large, while the transmissibility at higher frequencies is quite low. As the damping is increased, the resonant peaks are attenuated, but vibration isolation is lost at high frequency. This illustrates the inherent trade-off between resonance control and high-frequency isolation associated with the design of passive suspension systems. It is also observed from the figure that when the damping ratio is 0.7 or greater, the resonant peaks are almost attenuated. Thus, the higher value of damping ratio is not necessary. It is noted that in Equation 5, the third term FMR is much greater than the other and the behavior of the MR damper can is similar to that of a dry friction damper. By introducing an equivalent damping coefficient Ceq such that the work per cycle due to this equivalent damping coefficient equals that due to the yield stress damping force of the MR damper, the following equation holds [8]:

C eq = 4 F MR Xωπ
(11)

or

F MR = Xωπ C eq 4 = Xωπξ km 2 = kXπξr 2
(12)
Figure 5
figure 5

Transmissibility of the tub assembly.

In the above, X is the magnitude of the tub vibration which is determined by

X = F 0 k 1 1 − r 2 2 + 2 ξr 2 = m u r 2 R u m 1 1 − r 2 2 + 2 ξr 2
(13)

From Equations 12 and 13, the required value of FMR can be determined from a required value of the damping ratio ξ as follows:

F MR = kπξ m u r 3 R u 2 m 1 1 − r 2 2 + 2 ξr 2
(14)

In this study, it is assumed that the spring stiffness is k = 10 kN/m, the mass of the suspended tub assembly is m = 40 kg, and the equivalent unbalanced mass is m u  = 10 kg located at the radius R u  = 0.15 m. With the required damping ratio ξ = 0.7, at the resonance r = 1 − ξ 2 , the required value of FMR can be calculated from Equation 14 which is around 150 N in this study.

Taking all above into consideration, the optimal design of the MR brake for the washing machine can be summarized as follows: Find the optimal value of significant dimensions of the MR damper such as the pole length L p , the housing thickness t h , the core radius R c , the width of the MR duct t d , the width of the coil W c , and the overall length of the valve structure L that

Minimize the viscous coefficient objective function , OBJ = c vis = 12 ηL π R d t d 3 A p − A s 2
Subjected to : F MR = 2 A p − A s 2.85 L p t d τ y ≥ 150 N

In order to obtain the optimal solution, a finite element analysis code integrated with an optimization tool is employed. In this study, the first-order method with the golden section algorithm of the ANSYS optimization tool is used. Figure 6 shows the flow chart to achieve optimal design parameters of the MR damper. Firstly, an analysis ANSYS file for solving the magnetic circuit of the damper and calculating the objective function is built using ANSYS parametric design language (APDL). In the analysis file, the design variables (DVs) such as the pole length L p , the housing thickness t h , the core radius R c , the width of the MR duct t d , the width of the coil W c , and the overall length of the valve structure L must be coded as variables and initial values are assigned to them. The geometric dimensions of the valve structure are varied during the optimization process; the meshing size therefore should be specified by the number of elements per line rather than the element size. Because the magnetic field intensity is not constant along the pole length, it is necessary to define paths along the MR active volume where magnetic flux passes. The average magnetic field intensity across the MR ducts (Hmr) is calculated by integrating the field intensity along the defined path then divided by the path length. Thus, the magnetic field intensity is determined as follows:

H mr = 1 L p ∫ 0 L p H s ds
(15)

where H(s) is the magnetic flux density and magnetic field intensity at each nodal point on the defined path.

Figure 6
figure 6

Flow chart to achieve optimal design parameters of the MR damper.

From the figure, it is observed that the optimization is started with the initial value of DVs. By executing the analysis file, first the magnetic field intensity is derived. Then the yield stress, yield stress damping force, and objective function are respectively calculated from Equations 9, 7, and 8. The ANSYS optimization tool then transforms the optimization problem with constrained design variables to an unconstrained one via penalty functions. The dimensionless, unconstrained objective function f is formulated as follows:

f x = OBJ OBJ 0 + ∑ i = 1 n P x i x i
(16)

where OBJ0 is the reference objective function value that is selected from the current group of design sets. P x i is the exterior penalty function for the design variable x i . For the initial iteration (j = 0), the search direction of DVs is assumed to be the negative of the gradient of the unconstrained objective function. Thus, the direction vector is calculated by

d 0 = − ∇ f x 0
(17)

The values of DVs in the next iteration (j + 1) is obtained from the following equation:

x j + 1 = x j + s j d j
(18)

where the line search parameter s j is calculated by using a combination of the golden section algorithm and a local quadratic fitting technique. The analysis file is then executed with the new values of DVs, and the convergence of the objective function is checked. If the convergence occurs, the values of DVs at this iteration are the optimum. If not, the subsequent iterations will be performed. In the subsequent iterations, the procedures are similar to those of the initial iteration except for that the direction vectors are calculated according to the Polak-Ribiere recursion formula as follows:

d j = − ∇ f x j + r j − 1 d j − 1
(19)
where r j − 1 = ∇ f x j − ∇ f x j − 1 T ∇ f x j ∇ f x j − 1 2 .
(20)

Results and discussion

In this study, an optimal design of the MR damper for the washing machine is performed based on the optimization problem developed in the ‘Optimal design of the MR damper’ section. It is assumed that the piston part and housing of the damper are made of commercial silicon steel, and the coil wires are sized as 21 gage (diameter = 0.511 mm) whose allowable working current is 2.5A. In the optimization, the applied current is assumed to be 2A. Figure 7 shows the optimal solution of the MR damper in case the pole length L p , the housing thickness t h , the core radius R c , the width of the MR duct t d , the width of the coil W d , and the overall length of the valve structure L are considered as design variables. From the figure, it can be found that with a convergence tolerance of 0.5%, the optimal process is converged after 14 iterations and the solution at the 14th iteration is considered as the optimal one. The optimal values of L p , t h , R c , t d , W d , and L, respectively, are 7.4, 2.5, 5.2, 2, 1.8, and 25 mm. It is noted that the optimal values of t d and t h are equal to their lower limits in this case. These limits are posed considering the stability and manufacturing cost of the damper. At the optimum, the viscous coefficient is significantly reduced up to 10 Ns/m from its initial value (290 Ns/m). When no current is applied to the coil, the damping ratio at the optimum is around 0.01 which is very small so that a very small transmissibility of vibration from the drum to the cabinet at high excitation frequency can be achieved. It is also observed from the figure that the yield stress force of the damper at the optimum is greater than the required value (150 N). The yield stress and power consumption as functions of the applied current of the optimized MR damper are shown in Figure 7c. From the figure, it is observed that the yield stress force almost increases proportionally to the applied current as the current varies from 0 to 2.5A.

Figure 7
figure 7

Optimal solution of the MR damper for the washing machine. (a) Design variables. (b) Viscous coefficient and yield stress force. c) Braking torque and power consumption vs. applied current.

In order to validate the above optimal results, experimental results of the optimized MR damper are obtained and presented. Figure 8 shows the experimental setup to test the performance of the optimized MR damper. In the figure, a crank-slider mechanism is employed to convert the rotary motion of the motor into the reciprocal motion of the damper shaft. The DC motor with a gearbox controlled by the computer is used to rotate the crank shaft at a constant angular speed of 0.5 rad/s. The damping force is measured by a load cell. The output signal from the load cell is then sent to the computer via the A/D converter for evaluation. Once the experiment process is stated, a step current signal from the computer is sent to the current amplifier. The output current from the amplifier, a step current of 2A, is applied to the coil of the damper.

Figure 8
figure 8

Experimental setup to test the performance of the optimized MR damper.

Figure 9 shows the step response of the prototype MR damper. It is noted that in this case, the angular velocity of the motor is kept constant at 0.5Ï€ rad/s and the step response of the MR damper is measured on the downward motion of the shaft. In the figure, the solid line presents the calculated damping force while the dash line presents the measured one. The results also show that, at the steady state, the average damping force is around 157 N which is a little greater than the calculated one (150 N). The difference may result from friction which is not taken into account in the modeling of the damper.

Figure 9
figure 9

Experimental step response of the optimized MR damper.

Conclusions

In this research, optimal design of a flow-mode MR damper to suppress the vibration of front-loaded washing machines was undertaken. After a brief introduction, the general governing equation of the washing tub assembly is derived. The position of the two springs and the two dampers to suppress the machine vibration was then optimally determined based on the governing equation. A configuration of the MR damper for the washing machine was then proposed, and the damping force was obtained based on the Bingham rheological model of the MR brake. The optimization problem for the damper was constructed such that the viscous coefficient of the damper is minimized and the yield stress force of the damper is greater than a required value that attenuates almost the resonant peak of the tub mechanism. The optimal design of the MR damper for a prototype washing machine was obtained based on the ANSYS finite element analysis of the MR damper magnetic circuit and first-order optimization method. The results show that performances of the MR damper are significantly improved with the proposed optimal design. A prototype of the optimized MR damper was manufactured, and experimental results on the performance of the prototype damper were obtained and presented. It was shown that the experimental results well agreed with the modeling one obtained from finite element analysis. As the second phase of this study, dynamics of the whole prototype washing machine equipped with the MR damper will be obtained and experimental results on the prototype washer will be conducted to evaluate the effectiveness of the optimized MR damper.

References

  1. Bae S, Lee JM, Kang YJ, Kane JS, Yun JR: Dynamic analysis of an automatic washing machine with a hydraulic balancer. Sound Vib 2002,257(1):3–18. 10.1006/jsvi.2001.4162

    Article  Google Scholar 

  2. Papadopoulos E, Papadimitriou I: Modeling, design and control of a portable washing machine during the spinning cycle. In Proceedings of the 2001 IEEE/ASME international conference on advanced intelligent mechatronics systems. Como, Italy; 2001. 8–11, July 2001

    Google Scholar 

  3. Lim HT, Jeong WB, Kim KJ: Dynamic modeling and analysis of drum-type washing machine. Int J Precis Eng Manuf 2010,11(3):407–417. 10.1007/s12541-010-0047-7

    Article  Google Scholar 

  4. Michael JC, David Carson J: MR fluid sponge devices and their use in vibration control of washing machines. In Proc. SPIE 4331. CA, USA: Newport Beach; 2001. March 4, 2001

    Google Scholar 

  5. Cristiano S, Fibio P, Sergio MS, Giuseppe F, Nicolas G: Control of magnetorheological dampers for vibration reduction in a washing machine. Mechatronics 2009, 19: 410–421. 10.1016/j.mechatronics.2008.09.006

    Article  Google Scholar 

  6. Aydar G, Evrensel CA, Gordaninejad F, Fuchs A: A low force magneto-rheological (MR) fluid damper: design, fabrication and characterization. J Intell Mater Syst Struct 2007,18(12):1155–1160. 10.1177/1045389X07083138

    Article  Google Scholar 

  7. Nguyen QH, Choi SB, Lee YS, Han MS: An analytical method for optimal design of MR valve structures. Smart Mater Struct 2009,18(9):385–402.

    Article  Google Scholar 

  8. Silva CW: Vibration: fundamentals and practices. New York: CRC; 2000.

    Google Scholar 

Download references

Acknowledgements

This research is funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 107.04.2011.07 and partially supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MEST) (No. 2010-0015090).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Quoc Hung Nguyen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

QHN carried out the research background of the study, modeling of the system and optimal design of MR dampers. NDN carried out the experiment. SBC participated in configuration proposal of the MR damper, providing experimental data of the MR fluid and review of the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nguyen, Q.H., Nguyen, N.D. & Choi, S.B. Optimal design and performance evaluation of a flow-mode MR damper for front-loaded washing machines. Asia Pac. J. Comput. Engin. 1, 3 (2014). https://doi.org/10.1186/2196-1166-1-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2196-1166-1-3

Keywords