Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

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

Optimal design and performance evaluation of a flow-mode MR damper for front-loaded washing machines nguyenhung2000vn@yahoo.com Department of Mechanical It is well known that the vibration of washing machines is a challenging issue to be Engineering, Industrial University of considered. This research work focuses on the optimal design of a flow-mode Ho Chi Minh City, Hochiminh-12 NVBao, Vietnam magneto-rheological (MR) damper that can replace the conventional passive damper Full list of author information is for a washing machine. Firstly, rigid mode vibration of the washing machine due to available at the end of the article 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. Keywords: Magneto-rheological; MR damper; Washing machine; Vibration control 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 © 2014 Nguyen et al.; licensee Springer. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 2 of 14 http://www.apjcen.com/content/1/1/3 tub balance to eliminate the source of vibration [1,2]. In this approach, one type of dy- namic 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 automatic- ally 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 vi- brations 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 ap- proach. 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 spin- ning 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 vi- bration 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 re- duce 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 op- timal 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 con- structed 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 sus- pended tub (basket) to store the water for washing linked to the cabinet by two springs Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 3 of 14 http://www.apjcen.com/content/1/1/3 Figure 1 3D schematic of the prototype washing machine. 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 as- sembly. 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) sim- plified schematic of the machine is depicted. From the figure, the following governing equation of the washing machine can be derived: 2 2 2 2 mu € þ cu _ sinðÞ φ þ β þ sinðÞ φ−β þ ku sinðÞ φ þ α þ sinðÞ φ−α ¼ F ðÞ t 1 2 u 2 1 ð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 esti- mated 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 is the excitation force due to unbalanced mass in the u direction, F =F cosωt=m ω R cosωt, in which m u 0 u u u and R 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 qffiffiffiffiffiffiffiffiffiffi ω ¼ ω 1−ξ ð2Þ d n qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 2 2 2 k½ sinðÞ φþα þ sinðÞ φ−α c½ sinðÞ φþβ þ sinðÞ φ−β 1 2 2 1 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi where ω ¼ and ξ ¼ : 2 2 2 mk sinðÞ φþα þ sinðÞ φ−α ½ 1 2 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 Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 4 of 14 http://www.apjcen.com/content/1/1/3 Figure 2 2D simplified schematic of the prototype washing machine. 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 α + α = 90° and β + β = 90°, Equation 1 can be simplified to yield 1 2 2 2 mu € þ cu _ þ ku ¼ F ð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 qffiffiffiffiffiffiffiffiffiffi ω ¼ ω 1−ξ ð4Þ d n qffiffiffi k c pffiffiffiffiffi where ω ¼ and ξ ¼ : 2 mk An inherent drawback of the conventional damper is its high transmissibility of vibra- tion 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. Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 5 of 14 http://www.apjcen.com/content/1/1/3 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 Figure 3 Schematic configuration of the MR damper. Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 6 of 14 http://www.apjcen.com/content/1/1/3 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. By neglecting the frictional force and assuming quasi-static behavior of the damper, the damping force can be calculated by [7] F ¼ P A þ C x_ þ F sgn x _ ð5Þ d a s vis p MR p where P , c , and F , respectively, are the pressure in the gas chamber, the viscous a vis MR coefficient, and the yield stress force of the MR damper which are determined as follows: P ¼ P ð6Þ a 0 V þ A x 0 s p 12ηL C ¼ A −A ð7Þ vis p s πR t 2:85L F ¼ 2 A −A τ ð8Þ MR p s y In the above, τ 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 is the average radius of the annular duct given by R = R − t − 0.5t . L and R are the overall length and outside d h d radius of the MR valve, respectively. t is the valve housing thickness, t is the annular h d duct gap, and L 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 in- dependent 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. Data Polynomial Curve Fit =0.044+0.4769H 10 2 3 -0.0016H +1.8007E-6H 0 50 100 150 200 250 300 Magnetic Field Intensity, H(kA/m) Figure 4 Yield stress of MR fluid as a function of magnetic field intensity. Yield Stress, (kPa) y Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 7 of 14 http://www.apjcen.com/content/1/1/3 By applying the least squares curve fitting method, the yield stress of the MR fluid can be approximately expressed by 2 3 τ ¼ C þ C H þ C H þ C H ð9Þ y 0 1 2 3 In Equation 9, the unit of the yield stress is kilopascal while that of the magnetic field intensity is kA/m. The coefficients C , C , C , and C are respectively identified as 0 1 2 3 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 devel- oped 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: sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 þðÞ 2ξr TR ¼ ð10Þ 2 2 ðÞ 1−r þðÞ 2ξr where r is the frequency ratio, r= ω/ω . 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 at- tenuated, 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 =0.05 =0.1 =0.2 =0.3 3 =0.4 =0.5 =0.6 =0.7 010 20 30 40 Excitation Frequency [Hz] Figure 5 Transmissibility of the tub assembly. Transmissibility Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 8 of 14 http://www.apjcen.com/content/1/1/3 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 F is much greater than the other and the behavior of the MR damper can is MR similar to that of a dry friction damper. By introducing an equivalent damping coeffi- cient C such that the work per cycle due to this equivalent damping coefficient equals eq that due to the yield stress damping force of the MR damper, the following equation holds [8]: 4jj F MR C ¼ ð11Þ eq Xωπ or pffiffiffiffiffiffiffi XωπC Xωπξ km kXπξr eq jj F ¼ ¼ ¼ ð12Þ MR 4 2 2 In the above, X is the magnitude of the tub vibration which is determined by sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi F 1 m r R 1 0 u u X ¼ ¼ ð13Þ 2 2 2 2 2 2 k m ðÞ 1−r þðÞ 2ξr ðÞ 1−r þðÞ 2ξr From Equations 12 and 13, the required value of F can be determined from a re- MR quired value of the damping ratio ξ as follows: sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi kπξm r R 1 u u jj F ¼ ð14Þ MR 2 2 2m ðÞ 1−r þðÞ 2ξr 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 = 10 kg located at the radius R = 0.15 m. With the required damping ratio ξ = 0.7, at the pffiffiffiffiffiffiffiffiffiffi resonance r ¼ 1−ξ , the required value of F can be calculated from Equation 14 MR which is around 150 N in this study. Taking all above into consideration, the optimal design of the MR brake for the wash- ing machine can be summarized as follows: Find the optimal value of significant di- mensions of the MR damper such as the pole length L , the housing thickness t , the p h core radius R , the width of the MR duct t , the width of the coil W , and the overall c d c length of the valve structure L that 12ηL Minimize the viscous coefficientðÞ objective function ; OBJ ¼ c ¼ A −A vis p s πR t 2:85L Subjected to : F ¼ 2 A −A τ ≥150 N MR p s y 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 , the housing thickness t , the core p h radius R , the width of the MR duct t , the width of the coil W , and the overall length c d c Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 9 of 14 http://www.apjcen.com/content/1/1/3 Figure 6 Flow chart to achieve optimal design parameters of the MR damper. 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 (H )is mr 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 ¼ HsðÞds ð15Þ mr where H(s) is the magnetic flux density and magnetic field intensity at each nodal point on the defined path. 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 cal- culated from Equations 9, 7, and 8. The ANSYS optimization tool then transforms the optimization problem with constrained design variables to an unconstrained one via Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 10 of 14 http://www.apjcen.com/content/1/1/3 penalty functions. The dimensionless, unconstrained objective function f is formulated as follows: OBJ fxðÞ ¼ þ P ðÞ x ð16Þ x i OBJ i¼1 where OBJ is the reference objective function value that is selected from the current group of design sets. P is the exterior penalty function for the design variable x.For the x i initial iteration (j = 0), the search direction of DVs is assumed to be the negative of the gra- dient of the unconstrained objective function. Thus, the direction vector is calculated by ðÞ 0 ðÞ 0 d ¼ −∇fx ð17Þ The values of DVs in the next iteration (j+1) is obtained from the following equation: ðÞ jþ1 ðÞ j ðÞ j x ¼ x þ s d ð18Þ where the line search parameter s is calculated by using a combination of the golden section algorithm and a local quadratic fitting technique. The analysis file is then exe- cuted 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 vec- tors are calculated according to the Polak-Ribiere recursion formula as follows: ðÞ j ðÞ j ðÞ j−1 d ¼ −∇fx þ r d ð19Þ j−1 ðÞ j ðÞ j−1 ðÞ j ∇fx −∇fx ∇fx where r ¼ : ð20Þ j−1 ðÞ j−1 jj ∇fxðÞ Results and discussion In this study, an optimal design of the MR damper for the washing machine is per- formed 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 , the housing thickness t , the core radius R , the width of the p h c MR duct t , the width of the coil W , and the overall length of the valve structure L are d d 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 , t , R , p h c t , W , and L, respectively, are 7.4, 2.5, 5.2, 2, 1.8, and 25 mm. It is noted that the opti- d d mal values of t and t are equal to their lower limits in this case. These limits are d h 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 Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 11 of 14 http://www.apjcen.com/content/1/1/3 L L W 10t p C d R t c h 2 4 6 8 10 12 14 Iteration design variables (a) Viscous Coefficient Yield Stress Force 2 4 6 8 10 12 14 Iteration (b) viscous coefficient and yield stress force 2.0 Yield Stress Force Applied Power 1.5 1.0 0.5 0 0.0 0.0 0.5 1.0 1.5 2.0 2.5 Applied Current (A) (c) braking torque and power consumption vs. applied current 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. 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 ob- served from the figure that the yield stress force of the damper at the optimum is Design Variables [mm] Yield Stress Force [N] Viscous Coefficient C [Ns/m] vis Yield Stress Force [N] Applied Power [W] Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 12 of 14 http://www.apjcen.com/content/1/1/3 Motor Driver DC Servo Translating Motor Computer Rod A/D Eccentric Cam Load Cell Current D/A Meter Amplifier Prototype MR damper Figure 8 Experimental setup to test the performance of the optimized MR damper. 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. 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 com- puter 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 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 0.0 0.5 1.0 1.5 2.0 Time [s] Figure 9 Experimental step response of the optimized MR damper. Damping Force [N] Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 13 of 14 http://www.apjcen.com/content/1/1/3 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. 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 gen- eral 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 ana- lysis of the MR damper magnetic circuit and first-order optimization method. The re- sults 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 model- ing one obtained from finite element analysis. As the second phase of this study, dy- namics 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 evalu- ate the effectiveness of the optimized MR damper. 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. 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 details Department of Mechanical Engineering, Industrial University of Ho Chi Minh City, Hochiminh-12 NVBao, Vietnam. Department of Mechanical Engineering, Inha University, Incheon 402-751, Korea. Received: 14 August 2013 Accepted: 9 December 2013 Published: 29 April 2014 References 1. Bae S, Lee JM, Kang YJ, Kane JS, Yun JR (2002) Dynamic analysis of an automatic washing machine with a hydraulic balancer. Sound Vib 257(1):3–18 2. Papadopoulos E, Papadimitriou I (2001) 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, 8-11, July 2001 3. Lim HT, Jeong WB, Kim KJ (2010) Dynamic modeling and analysis of drum-type washing machine. Int J Precis Eng Manuf 11(3):407–417 Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 14 of 14 http://www.apjcen.com/content/1/1/3 4. Michael JC, David Carson J (2001) MR fluid sponge devices and their use in vibration control of washing machines. In: Proc. SPIE 4331. Newport Beach, CA, USA, March 4, 2001 5. Cristiano S, Fibio P, Sergio MS, Giuseppe F, Nicolas G (2009) Control of magnetorheological dampers for vibration reduction in a washing machine. Mechatronics 19:410–421 6. Aydar G, Evrensel CA, Gordaninejad F, Fuchs A (2007) A low force magneto-rheological (MR) fluid damper: design, fabrication and characterization. J Intell Mater Syst Struct 18(12):1155–1160 7. Nguyen QH, Choi SB, Lee YS, Han MS (2009) An analytical method for optimal design of MR valve structures. Smart Mater Struct 18(9):385–402 8. Silva CW (2000) Vibration: fundamentals and practices. CRC, New York doi:10.1186/2196-1166-1-3 Cite this article as: Nguyen et al.: Optimal design and performance evaluation of a flow-mode MR damper for front-loaded washing machines. Asia Pacific Journal on Computational Engineering 2014 1:3. Submit your manuscript to a journal and benefi t from: 7 Convenient online submission 7 Rigorous peer review 7 Immediate publication on acceptance 7 Open access: articles freely available online 7 High visibility within the fi eld 7 Retaining the copyright to your article Submit your next manuscript at 7 springeropen.com http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Asia Pacific Journal on Computational Engineering Springer Journals

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

Download PDF
14 pages

Loading next page...
 
/lp/springer-journals/optimal-design-and-performance-evaluation-of-a-flow-mode-mr-damper-for-x5rx3FbtNO
Publisher
Springer Journals
Copyright
Copyright © 2014 by Nguyen et al.; licensee Springer.
Subject
Engineering; Theoretical and Applied Mechanics; Computational Science and Engineering; Classical Continuum Physics; Mathematical Applications in the Physical Sciences
eISSN
2196-1166
DOI
10.1186/2196-1166-1-3
Publisher site
See Article on Publisher Site

Abstract

nguyenhung2000vn@yahoo.com Department of Mechanical It is well known that the vibration of washing machines is a challenging issue to be Engineering, Industrial University of considered. This research work focuses on the optimal design of a flow-mode Ho Chi Minh City, Hochiminh-12 NVBao, Vietnam magneto-rheological (MR) damper that can replace the conventional passive damper Full list of author information is for a washing machine. Firstly, rigid mode vibration of the washing machine due to available at the end of the article 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. Keywords: Magneto-rheological; MR damper; Washing machine; Vibration control 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 © 2014 Nguyen et al.; licensee Springer. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 2 of 14 http://www.apjcen.com/content/1/1/3 tub balance to eliminate the source of vibration [1,2]. In this approach, one type of dy- namic 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 automatic- ally 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 vi- brations 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 ap- proach. 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 spin- ning 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 vi- bration 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 re- duce 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 op- timal 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 con- structed 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 sus- pended tub (basket) to store the water for washing linked to the cabinet by two springs Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 3 of 14 http://www.apjcen.com/content/1/1/3 Figure 1 3D schematic of the prototype washing machine. 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 as- sembly. 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) sim- plified schematic of the machine is depicted. From the figure, the following governing equation of the washing machine can be derived: 2 2 2 2 mu € þ cu _ sinðÞ φ þ β þ sinðÞ φ−β þ ku sinðÞ φ þ α þ sinðÞ φ−α ¼ F ðÞ t 1 2 u 2 1 ð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 esti- mated 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 is the excitation force due to unbalanced mass in the u direction, F =F cosωt=m ω R cosωt, in which m u 0 u u u and R 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 qffiffiffiffiffiffiffiffiffiffi ω ¼ ω 1−ξ ð2Þ d n qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 2 2 2 k½ sinðÞ φþα þ sinðÞ φ−α c½ sinðÞ φþβ þ sinðÞ φ−β 1 2 2 1 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi where ω ¼ and ξ ¼ : 2 2 2 mk sinðÞ φþα þ sinðÞ φ−α ½ 1 2 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 Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 4 of 14 http://www.apjcen.com/content/1/1/3 Figure 2 2D simplified schematic of the prototype washing machine. 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 α + α = 90° and β + β = 90°, Equation 1 can be simplified to yield 1 2 2 2 mu € þ cu _ þ ku ¼ F ð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 qffiffiffiffiffiffiffiffiffiffi ω ¼ ω 1−ξ ð4Þ d n qffiffiffi k c pffiffiffiffiffi where ω ¼ and ξ ¼ : 2 mk An inherent drawback of the conventional damper is its high transmissibility of vibra- tion 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. Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 5 of 14 http://www.apjcen.com/content/1/1/3 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 Figure 3 Schematic configuration of the MR damper. Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 6 of 14 http://www.apjcen.com/content/1/1/3 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. By neglecting the frictional force and assuming quasi-static behavior of the damper, the damping force can be calculated by [7] F ¼ P A þ C x_ þ F sgn x _ ð5Þ d a s vis p MR p where P , c , and F , respectively, are the pressure in the gas chamber, the viscous a vis MR coefficient, and the yield stress force of the MR damper which are determined as follows: P ¼ P ð6Þ a 0 V þ A x 0 s p 12ηL C ¼ A −A ð7Þ vis p s πR t 2:85L F ¼ 2 A −A τ ð8Þ MR p s y In the above, τ 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 is the average radius of the annular duct given by R = R − t − 0.5t . L and R are the overall length and outside d h d radius of the MR valve, respectively. t is the valve housing thickness, t is the annular h d duct gap, and L 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 in- dependent 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. Data Polynomial Curve Fit =0.044+0.4769H 10 2 3 -0.0016H +1.8007E-6H 0 50 100 150 200 250 300 Magnetic Field Intensity, H(kA/m) Figure 4 Yield stress of MR fluid as a function of magnetic field intensity. Yield Stress, (kPa) y Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 7 of 14 http://www.apjcen.com/content/1/1/3 By applying the least squares curve fitting method, the yield stress of the MR fluid can be approximately expressed by 2 3 τ ¼ C þ C H þ C H þ C H ð9Þ y 0 1 2 3 In Equation 9, the unit of the yield stress is kilopascal while that of the magnetic field intensity is kA/m. The coefficients C , C , C , and C are respectively identified as 0 1 2 3 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 devel- oped 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: sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 þðÞ 2ξr TR ¼ ð10Þ 2 2 ðÞ 1−r þðÞ 2ξr where r is the frequency ratio, r= ω/ω . 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 at- tenuated, 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 =0.05 =0.1 =0.2 =0.3 3 =0.4 =0.5 =0.6 =0.7 010 20 30 40 Excitation Frequency [Hz] Figure 5 Transmissibility of the tub assembly. Transmissibility Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 8 of 14 http://www.apjcen.com/content/1/1/3 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 F is much greater than the other and the behavior of the MR damper can is MR similar to that of a dry friction damper. By introducing an equivalent damping coeffi- cient C such that the work per cycle due to this equivalent damping coefficient equals eq that due to the yield stress damping force of the MR damper, the following equation holds [8]: 4jj F MR C ¼ ð11Þ eq Xωπ or pffiffiffiffiffiffiffi XωπC Xωπξ km kXπξr eq jj F ¼ ¼ ¼ ð12Þ MR 4 2 2 In the above, X is the magnitude of the tub vibration which is determined by sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi F 1 m r R 1 0 u u X ¼ ¼ ð13Þ 2 2 2 2 2 2 k m ðÞ 1−r þðÞ 2ξr ðÞ 1−r þðÞ 2ξr From Equations 12 and 13, the required value of F can be determined from a re- MR quired value of the damping ratio ξ as follows: sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi kπξm r R 1 u u jj F ¼ ð14Þ MR 2 2 2m ðÞ 1−r þðÞ 2ξr 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 = 10 kg located at the radius R = 0.15 m. With the required damping ratio ξ = 0.7, at the pffiffiffiffiffiffiffiffiffiffi resonance r ¼ 1−ξ , the required value of F can be calculated from Equation 14 MR which is around 150 N in this study. Taking all above into consideration, the optimal design of the MR brake for the wash- ing machine can be summarized as follows: Find the optimal value of significant di- mensions of the MR damper such as the pole length L , the housing thickness t , the p h core radius R , the width of the MR duct t , the width of the coil W , and the overall c d c length of the valve structure L that 12ηL Minimize the viscous coefficientðÞ objective function ; OBJ ¼ c ¼ A −A vis p s πR t 2:85L Subjected to : F ¼ 2 A −A τ ≥150 N MR p s y 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 , the housing thickness t , the core p h radius R , the width of the MR duct t , the width of the coil W , and the overall length c d c Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 9 of 14 http://www.apjcen.com/content/1/1/3 Figure 6 Flow chart to achieve optimal design parameters of the MR damper. 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 (H )is mr 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 ¼ HsðÞds ð15Þ mr where H(s) is the magnetic flux density and magnetic field intensity at each nodal point on the defined path. 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 cal- culated from Equations 9, 7, and 8. The ANSYS optimization tool then transforms the optimization problem with constrained design variables to an unconstrained one via Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 10 of 14 http://www.apjcen.com/content/1/1/3 penalty functions. The dimensionless, unconstrained objective function f is formulated as follows: OBJ fxðÞ ¼ þ P ðÞ x ð16Þ x i OBJ i¼1 where OBJ is the reference objective function value that is selected from the current group of design sets. P is the exterior penalty function for the design variable x.For the x i initial iteration (j = 0), the search direction of DVs is assumed to be the negative of the gra- dient of the unconstrained objective function. Thus, the direction vector is calculated by ðÞ 0 ðÞ 0 d ¼ −∇fx ð17Þ The values of DVs in the next iteration (j+1) is obtained from the following equation: ðÞ jþ1 ðÞ j ðÞ j x ¼ x þ s d ð18Þ where the line search parameter s is calculated by using a combination of the golden section algorithm and a local quadratic fitting technique. The analysis file is then exe- cuted 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 vec- tors are calculated according to the Polak-Ribiere recursion formula as follows: ðÞ j ðÞ j ðÞ j−1 d ¼ −∇fx þ r d ð19Þ j−1 ðÞ j ðÞ j−1 ðÞ j ∇fx −∇fx ∇fx where r ¼ : ð20Þ j−1 ðÞ j−1 jj ∇fxðÞ Results and discussion In this study, an optimal design of the MR damper for the washing machine is per- formed 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 , the housing thickness t , the core radius R , the width of the p h c MR duct t , the width of the coil W , and the overall length of the valve structure L are d d 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 , t , R , p h c t , W , and L, respectively, are 7.4, 2.5, 5.2, 2, 1.8, and 25 mm. It is noted that the opti- d d mal values of t and t are equal to their lower limits in this case. These limits are d h 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 Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 11 of 14 http://www.apjcen.com/content/1/1/3 L L W 10t p C d R t c h 2 4 6 8 10 12 14 Iteration design variables (a) Viscous Coefficient Yield Stress Force 2 4 6 8 10 12 14 Iteration (b) viscous coefficient and yield stress force 2.0 Yield Stress Force Applied Power 1.5 1.0 0.5 0 0.0 0.0 0.5 1.0 1.5 2.0 2.5 Applied Current (A) (c) braking torque and power consumption vs. applied current 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. 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 ob- served from the figure that the yield stress force of the damper at the optimum is Design Variables [mm] Yield Stress Force [N] Viscous Coefficient C [Ns/m] vis Yield Stress Force [N] Applied Power [W] Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 12 of 14 http://www.apjcen.com/content/1/1/3 Motor Driver DC Servo Translating Motor Computer Rod A/D Eccentric Cam Load Cell Current D/A Meter Amplifier Prototype MR damper Figure 8 Experimental setup to test the performance of the optimized MR damper. 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. 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 com- puter 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 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 0.0 0.5 1.0 1.5 2.0 Time [s] Figure 9 Experimental step response of the optimized MR damper. Damping Force [N] Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 13 of 14 http://www.apjcen.com/content/1/1/3 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. 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 gen- eral 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 ana- lysis of the MR damper magnetic circuit and first-order optimization method. The re- sults 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 model- ing one obtained from finite element analysis. As the second phase of this study, dy- namics 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 evalu- ate the effectiveness of the optimized MR damper. 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. 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 details Department of Mechanical Engineering, Industrial University of Ho Chi Minh City, Hochiminh-12 NVBao, Vietnam. Department of Mechanical Engineering, Inha University, Incheon 402-751, Korea. Received: 14 August 2013 Accepted: 9 December 2013 Published: 29 April 2014 References 1. Bae S, Lee JM, Kang YJ, Kane JS, Yun JR (2002) Dynamic analysis of an automatic washing machine with a hydraulic balancer. Sound Vib 257(1):3–18 2. Papadopoulos E, Papadimitriou I (2001) 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, 8-11, July 2001 3. Lim HT, Jeong WB, Kim KJ (2010) Dynamic modeling and analysis of drum-type washing machine. Int J Precis Eng Manuf 11(3):407–417 Nguyen et al. Asia Pacific Journal on Computational Engineering 2014, 1:3 Page 14 of 14 http://www.apjcen.com/content/1/1/3 4. Michael JC, David Carson J (2001) MR fluid sponge devices and their use in vibration control of washing machines. In: Proc. SPIE 4331. Newport Beach, CA, USA, March 4, 2001 5. Cristiano S, Fibio P, Sergio MS, Giuseppe F, Nicolas G (2009) Control of magnetorheological dampers for vibration reduction in a washing machine. Mechatronics 19:410–421 6. Aydar G, Evrensel CA, Gordaninejad F, Fuchs A (2007) A low force magneto-rheological (MR) fluid damper: design, fabrication and characterization. J Intell Mater Syst Struct 18(12):1155–1160 7. Nguyen QH, Choi SB, Lee YS, Han MS (2009) An analytical method for optimal design of MR valve structures. Smart Mater Struct 18(9):385–402 8. Silva CW (2000) Vibration: fundamentals and practices. CRC, New York doi:10.1186/2196-1166-1-3 Cite this article as: Nguyen et al.: Optimal design and performance evaluation of a flow-mode MR damper for front-loaded washing machines. Asia Pacific Journal on Computational Engineering 2014 1:3. Submit your manuscript to a journal and benefi t from: 7 Convenient online submission 7 Rigorous peer review 7 Immediate publication on acceptance 7 Open access: articles freely available online 7 High visibility within the fi eld 7 Retaining the copyright to your article Submit your next manuscript at 7 springeropen.com

Journal

Asia Pacific Journal on Computational Engineering – Springer Journals

Published: Apr 29, 2014

References

  • Dynamic analysis of an automatic washing machine with a hydraulic balancer
    Bae, S; Lee, JM; Kang, YJ; Kane, JS; Yun, JR
  • Proceedings of the 2001 IEEE/ASME international conference on advanced intelligent mechatronics systems
    Papadopoulos, E; Papadimitriou, I
  • Dynamic modeling and analysis of drum-type washing machine
    Lim, HT; Jeong, WB; Kim, KJ
  • Proc. SPIE 4331
    Michael, JC; David Carson, J
  • Control of magnetorheological dampers for vibration reduction in a washing machine
    Cristiano, S; Fibio, P; Sergio, MS; Giuseppe, F; Nicolas, G
  • A low force magneto-rheological (MR) fluid damper: design, fabrication and characterization
    Aydar, G; Evrensel, CA; Gordaninejad, F; Fuchs, A
  • An analytical method for optimal design of MR valve structures
    Nguyen, QH; Choi, SB; Lee, YS; Han, MS
  • Vibration: fundamentals and practices
    Silva, CW