Optimal design of a double coil magnetorheological fluid damper with various piston profiles

Similar documents
Research on Damping Characteristics of Magneto-rheological Damper Used in Vehicle Seat Suspension

A STUDY OF HYDRAULIC RESISTANCE OF VISCOUS BYPASS GAP IN MAGNETORHEOLOGICAL DAMPER

1576. Development of a variable-damping magnetorheological damper with multiple poles

Parameter optimisation design for a six-dof heavy duty vehicle seat suspension

Magneto Rheological Semi-Active Damper with External By-pass Circuit in Modular Structure

SMART FLUID SELF ADAPTIVE DAMPER SYSTEM (SFSADS)

A REVIEW OF TWO WHEELER VEHICLES REAR SHOCK ABSORBER

MAGNETORHEOLOGICAL FLUID DAMPER WITH EXTERNAL MAGNETIC CIRCUIT. Alexandru Boltosi, Adrian Chiriac, Liviu Bereteu, Ramona Nagy, Voichiţa Hule

Design of closing electromagnet of high power spring operating mechanism

Magnetorheological suspension damper for space application

Forced vibration frequency response for a permanent magnetic planetary gear

Performance analysis of a magnetorheological damper with energy harvesting ability

Fig Electromagnetic Actuator

Investigation & Analysis of Three Phase Induction Motor Using Finite Element Method for Power Quality Improvement

Influence of Cylinder Bore Volume on Pressure Pulsations in a Hermetic Reciprocating Compressor

Theoretical and Experimental Investigation of Compression Loads in Twin Screw Compressor

Analysis and measurement of damping characteristics of linear generator

DESIGN AND PERFORMANCE ANALYSIS OF MAGNETORHEOLOGICAL VALVE MODULE WITH ANNULAR-RADIAL CONCEPT

Experimental Investigation of Effects of Shock Absorber Mounting Angle on Damping Characterstics

Transient analysis of a new outer-rotor permanent-magnet brushless DC drive using circuit-field-torque coupled timestepping finite-element method

Study on the Controllable Characteristics of Electrorheological Valve Using Serial Multielectrode

Semi-Active Suspension for an Automobile

A STUDY OF THE CENTRIFUGAL COMPRESSOR DISCHARGE PIPELINE CONSTRAINED OSCILLATION. KIRILL SOLODYANKIN*, JIŘÍ BĚHAL ČKD KOMPRESORY, a.s.

Performance prediction of magnetorheological valves under various type of fluid and flux path

Performance of a MR Hydraulic Power Actuation System *

Structural Analysis of Pick-Up Truck Chassis using Fem

Research on the Structure of Linear Oscillation Motor and the Corresponding Applications on Piston Type Refrigeration Compressor

Available online at ScienceDirect. Procedia CIRP 33 (2015 )

Parametric Modeling and Finite Element Analysis of the Brake Drum Based on ANSYS APDL

KINEMATICAL SUSPENSION OPTIMIZATION USING DESIGN OF EXPERIMENT METHOD

Integrated active and semi-active control for seat suspension of a heavy duty vehicle

China. Keywords: Electronically controled Braking System, Proportional Relay Valve, Simulation, HIL Test

EXPERIMENTAL ANALYSIS OF MAGNETO-RHEOLOGICAL FLUID BASED FRONT FORK DAMPER USING ELECTRO-MAGNETIC CONCEPT

CHAPTER 4 MR DAMPER DESIGN. In this chapter, details of MR damper geometry and magnetic circuit design are provided.

Chapter 2 Dynamic Analysis of a Heavy Vehicle Using Lumped Parameter Model

Comparing FEM Transfer Matrix Simulated Compressor Plenum Pressure Pulsations to Measured Pressure Pulsations and to CFD Results

Optimization of Three-stage Electromagnetic Coil Launcher

INTERNATIONAL JOURNAL OF DESIGN AND MANUFACTURING TECHNOLOGY (IJDMT) CONSTANT SPEED ENGINE CONROD SOFT VALIDATION & OPTIMIZATION

Research on Optimization for the Piston Pin and the Piston Pin Boss

The Effects of Magnetic Circuit Geometry on Torque Generation of 8/14 Switched Reluctance Machine

EVS28 KINTEX, Korea, May 3-6, 2015

Multi Body Dynamic Analysis of Slider Crank Mechanism to Study the effect of Cylinder Offset

INWHEEL SRM DESIGN WITH HIGH AVERAGE TORQUE AND LOW TORQUE RIPPLE

Optimization of Scissor-jack-Damper s Parameters and Performance under the Constrain of Human Comfort

MARINE FOUR-STROKE DIESEL ENGINE CRANKSHAFT MAIN BEARING OIL FILM LUBRICATION CHARACTERISTIC ANALYSIS

Application of Simulation-X R based Simulation Technique to Notch Shape Optimization for a Variable Swash Plate Type Piston Pump

EFFECT OF SURFACE ROUGHNESS ON PERFORMANCE OF WIND TURBINE

MODELING SUSPENSION DAMPER MODULES USING LS-DYNA

ANALYSIS OF THE INFLUENCE OF HYDRAULIC CYLINDER DIAMETER TO THE TOTAL DAMPING FORCE AND THE GENERATED ELECTRICITY OF REGENERATIVE SHOCK ABSORBER

Development of an innovative diaphragm accumulator design and assembly process

The Effect of Friction between a Cylindrical Guide and Magnetic Tape on Lateral Tape Motion

1874. Effect predictions of star pinion geometry phase adjustments on dynamic load sharing behaviors of differential face gear trains

Modal Analysis of Automobile Brake Drum Based on ANSYS Workbench Dan Yang1, 2,Zhen Yu1, 2, Leilei Zhang1, a * and Wentao Cheng2

Chapter 5. Design of Control Mechanism of Variable Suspension System. 5.1: Introduction: Objective of the Mechanism:

The Optimal Design of a Drum Friction Plate Using AnsysWorkbench

Full Scale Experimental Evaluation for Cable Dampers

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

Static And Modal Analysis of Tractor Power Take Off (PTO) Gearbox Housing

Chapter 2. Background

Effect of concave plug shape of a control valve on the fluid flow characteristics using computational fluid dynamics

Aerodynamic Characteristics of Sedan with the Rolling Road Ground Effect Simulation System

1036. Thermal-hydraulic modelling and analysis of hydraulic damper for impact cylinder with large flow

CHAPTER 4 MODELING OF PERMANENT MAGNET SYNCHRONOUS GENERATOR BASED WIND ENERGY CONVERSION SYSTEM

Optimal Design of a New Magnetorheological Mount Using Flexible Valve-Type Structure

STABILITY ENHANCEMENT OF RAILWAY VEHICLE DYNAMICS PERFORMANCE IN LATERAL DIRECTION USING FUZZY BOGIE-BASED SKYHOOK CONTROL

Improving muffler performance using simulation-based design

Effect of Permanent Magnet Rotor Design on PMSM Properties

Investigating the Benefit of Using Magnetorheological Fluids in a Shock Absorber

ONLINE NON-CONTACT TORSION SENSING METHOD USING FIBER BRAGG GRATING SENSORS AND OPTICAL COUPLING METHOD. Yoha Hwang and Jong Min Lee

Vehicle Dynamic Simulation Using A Non-Linear Finite Element Simulation Program (LS-DYNA)

THE DYNAMIC CHARACTERISTICS OF A DIRECT-ACTING WATER HYDRAULIC RELIEF VALVE WITH DOUBLE DAMPING: NUMERICAL AND EXPERIMENTAL INVESTIGATION

Transverse Distribution Calculation and Analysis of Strengthened Yingjing Bridge

MULTI-PARAMETER OPTIMIZATION OF BRAKE OF PISTON

FINITE ELEMENT SIMULATION OF SHOT PEENING AND STRESS PEEN FORMING

Research and Development of Mechanically Adjustable Fluid Viscous Damper Dan-Feng SONG*, Yong-Jin LU

The operating principle and experimental verification of the hydraulic electromagnetic energy-regenerative shock absorber

Electromagnetic Field Analysis for Permanent Magnet Retarder by Finite Element Method

A Linear Magnetic-geared Free-piston Generator for Range-extended Electric Vehicles

Driver roll speed influence in Ring Rolling process

This chapter gives details of the design, development, and characterization of the

Load Analysis and Multi Body Dynamics Analysis of Connecting Rod in Single Cylinder 4 Stroke Engine

AN OPTIMAL PROFILE AND LEAD MODIFICATION IN CYLINDRICAL GEAR TOOTH BY REDUCING THE LOAD DISTRIBUTION FACTOR

Modeling and Vibration Analysis of a Drum type Washing Machine

Study of a Novel Compliant Suspension Mechanism in Low Side Type Scroll Compressor

Study Of Static And Frequency Responsible Analysis Of Hangers With Exhaust System

ELECTROMAGNETS ARRANGEMENT FOR ELECTROMAGNETIC WINDSHIELD WIPERS - PROPOSAL AND ANALYSIS

Reducing Noise and Vibration of Hydraulic Hybrid And Plug-In Hybrid Electric Vehicles

College of Mechanical & Power Engineering Of China Three Gorges University, Yichang, Hubei Province, China

Plastic Ball Bearing Design Improvement Using Finite Element Method

Evaluation of the Fatigue Life of Aluminum Bogie Structures for the Urban Maglev

Experimental investigation on vibration characteristics and frequency domain of heavy haul locomotives

Low Speed Control Enhancement for 3-phase AC Induction Machine by Using Voltage/ Frequency Technique

AXLE HOUSING AND UNITIZE BEARING PACK SET MODAL CHARACTERISATION

Experimental Investigation of Using MR Fluids in Automobiles Suspension Systems

Turbostroje 2015 Návrh spojení vysokotlaké a nízkotlaké turbíny. Turbomachinery 2015, Design of HP and LP turbine connection

Development of a low voltage Dielectric Electro-Active Polymer actuator

MODELING, TESTING AND EVALUATION OF MAGNETO-RHEOLOGICAL SHOCK ABSORBER

DESIGN OF COMPACT PERMANENT-MAGNET SYNCHRONOUS MOTORS WITH CONCENTRATED WINDINGS

Application of ABAQUS to Analyzing Shrink Fitting Process of Semi Built-up Type Marine Engine Crankshaft

Transcription:

11 th World Congress on Structural and Multidisciplinary Optimisation 07 th -12 th, June 2015, Sydney Australia Optimal design of a double coil magnetorheological fluid damper with various piston profiles Guoliang Hu 1, 2, Zheng Xie 1, Weihua Li 2, * 1 School of Mechatronic Engineering, East China Jiaotong University, Nanchang City, Jiangxi Province, China, glhu@ecjtu.edu.cn 2 School of Mechanical, Materials and Mechatronic Engineering, University of Wollongong, NSW, Australia, weihuali@uow.edu.au Abstract A magnetorheological (MR) damper is one of the most advanced devices used in a semi-active control system to mitigate unwanted vibration because the damping force can be controlled by changing the viscosity of the internal magnetorheological fluids (MRF). This study proposes a typical double coil MR damper where the damping force and dynamic range were derived from a quasi-static model based on the Bingham model of MR fluid. A finite element model was built to study the performance of this double coil MR damper by investigating seven different piston configurations, including the numbers and shapes of their chamfered ends. The objective function of an optimisation problem was proposed and then an optimisation procedure was constructed using the ANSYS parametric design language (APDL) to obtain the optimal value of a double coil MR damper. Furthermore, an experimental analysis was also carried out. These results were then compared with the optimised MR damper s simulation results, which clearly validated the simulated results. The relevant results of this analysis can easily be extended to other MR dampers. Keywords: MRF, MR damper, Double coil, Finite element analysis, Optimal design 1. Introduction The rheological properties of a magnetorheological fluid (MRF) can be continuously changed within several milliseconds by applying or removing external magnetic fields. These unique features have led to the development of many MRF-based devices such as the MR damper, MR valve, MR brake, MR clutch, and so on. Recently, some studies in literature have focused on the geometric optimisation of an MR damper, with the results showing that their performance can be improved significantly by optimising the design of the magnetic circuit that controls these systems. Zhang et al. [1] proposed the use of finite elements to improve the magnetic design of an MR damper. Kham et al. [2] used finite element software to simulate nine different configurations of the pistons for an MR damper and investigate how these configurations would affect and influence the maximum pressure drop; the results showed that a single coil piston with chamfered ends was better than the other configurations for the same magnitude of input current and piston velocity. Ferdaus et al. [] established 2D axi-symmetric and a D model of an MR damper that considered the shape of the piston, the MR fluid gap, the air gap, and the thickness of the damper s housing. All these models were simulated with different currents, different piston velocities, and different strokes. Parlak et al. [4] investigated the geometrical optimisation of an MR shock damper using the Taguchi experimental design approach by specifying four parameters (gap, flange thickness, radius of piston core, and current excitation) and by selecting the maximum dynamic range required as the target value; this analysis was carried out using analytical equations rather than experimental data. Parlak et al. [5] presented a method for optimising the design of the target damper force and maximum magnetic flux density of an MR damper; this new approach used an electromagnetic analysis of the magnetic field and a CFD analysis of MR flow together to obtain the optimal value of the design parameters. Nguyen and Choi [6] presented an optimal design of a passenger vehicle MR damper that was constrained in a specific cylindrical volume, and an advanced objective function that collectively included the damping force, the dynamic range, and an inductive time constant. Nguyen and Choi [7] also proposed two types of shear mode MR dampers for a front loader washing machine, where an optimisation methodology based on a finite element analysis integrated with an optimisation tool was used to obtain the optimal geometric dimensions of the MR dampers; the results showed that an MR damper with three coils and without a non-magnetic bobbin was the best configuration for this application. In this study a double coil MR damper was developed and prototyped, and the damping force and dynamic range were also derived. A finite element model was built to investigate the performance of the double coil MR damper by considering seven pistons with different configurations. An optimisation procedure was then constructed with the ANSYS parametric design language (APDL) to obtain the optimal parameters of the double coil MR damper. Finally, a series of dynamic experimental tests were also carried out. 1

2. Design considerations for a double coil MR damper with an annular gap Figure 1 shows a schematic for the proposed double coil MR damper under consideration. Two chambers in the cylinder are separated by a floating piston. The section with the piston head is filled with MR fluid and the accumulator that compensates for the changes in volume induced as the piston rod moves is filled with pressurised nitrogen gas. As the damper piston rod moves, the MR fluid flows through the annular gap to the other side of the piston. There is a double coil of wire inside the piston head used for winding is heat resistant and electrically insulated. When a direct current is applied to the double coil, a magnetic field occurs around the piston head. It is noted that the direction of the current applied onto the double coil can be the same or the reverse, and it can enlarge the maximum damping force or the dynamic range to some extent. Figure 2 shows the magnetic circuit of the double coil MR damper with an annular gap at both ends of the flanges and in the middle of the flange, where the flux lines are perpendicular to the flow direction and caused a field-dependent resistance to the fluid flow. The double coil MR damper is shaped to guide the magnetic flux axially through the damper core, across the length of the core flange and the gap at one end, then on through the flux return and across the gap and core flange at the opposite end. The volume of fluid through which the magnetic field passes is defined as an active volume, and the MR effects only occur within the active volume. The most effective MR dampers have a high magnetic flux density passing through a large active volume, but a lot of magnetic coils are needed to produce large magnetic fields. An optimised circuit would maintain a balance between the field produced and power required by the magnetic coils. Accumlator Floating piston Annular gap Coil Figure 1: Schematic diagram of a double coil MR damper with annular gap Damper housing h Damper core L R Magnetic circuit I 1 I 2 Rc R1 Rh wc MR flow L 1 L c L 2 L c L 1 Figure 2: Magnetic circuit of the double coil MR damper with annular gap Some of the important dimensions of double coil MR damper are also listed in Fig.2. Damper geometry is characterised by the length of the gap L, the length of the end flange L 1, the length of the middle flange L 2, the length of the coil L c, the thickness of the piston head housing R h, the width of the annular gap h, the radius of the piston head R, the internal radius of the piston core R c, the external radius of the piston core R 1, and the width of the coil W c. The total damping force F generated by the double coil MR damper consists of three components: the field-dependent force F τ due to the magnetic field, the viscous force F η due to the viscous effects, and the frictional force F f. F = Fτ + Fη + F f (1) where 2

6ηLqAp Fη = (2) π ( R h 1 + 0.5h) L1 L2 Fτ = 2c1 Apτ y1 sgn( vp) + c2 Apτ y2 sgn( vp) () h h where A p is the cross-sectional area of the piston head, and η is the plastic viscosity. q is the flow rate through the double coil MR damper, and it can be calculated from the velocity of the piston ν p. τ y1 and τ y2 are the yield stresses of the MR fluid in the end flange and middle flange, respectively. c 1 and c 2 are coefficients which depended on the flow velocity profile, and have a value ranging from a minimum value of 2.07 to a maximum value of.07. As shown in Eq. (1), the first term is called the controllable force because it varies with the applied field, whereas the sum of the latter two terms is referred to as the uncontrollable force because they generate a constant force according to the velocity of the piston. The dynamic range D is defined as the ratio of the total damping force to the uncontrollable force, and it is given by Fτ + Fη + F f Fτ D = = 1 + (6) F + F F + F η f Here, the dynamic range D was introduced to evaluate the performance of the double coil MR damper. Normally, it is better to keep the dynamic range as large as possible to maximise the effectiveness of the MR damper. The dynamic range is proportional to the shear force and is a function of the size of the gap. The width of the annular gap h is inversely proportional to the controllable force, so a small gap width will increase the range of the controllable force, but when h is less than 0.5 mm, the viscous force F η is much faster than the controllable force F τ this reduces the dynamic range. Again when the gap becomes wider both F τ and F η fall, so finding an optimal width gap that maximises the dynamic range is very important. Moreover, parameters such as the length of the flange L 1 and L 2, the radius of the piston head R, the yield stress τ y1 and τ y2, the thickness of the piston head housing R h and the internal radius of the core of the piston R c will also play an important role in searching for the right design. η f. Modelling an optimal design of a double coil MR damper using the finite element method The magnetic field in an MR damper is produced by an electromagnet. Here in the ANSYS simulation model the excitation coil is considered to be the electromagnet, and the magnetic field provided by this excitation coil is needed to energise the MR fluid. By varying the current through the excitation coil the density of the magnetic flux can be varied and the MR fluid is energised accordingly. To investigate how the shape of the piston affects the performance of a double coil MR damper, seven models with different shaped pistons were designed, as shown in Fig.. The model 1 piston was defined as having a plain end, the model 2 piston was defined as have each end of the piston chamfered, the model piston was defined as having a radius on each end, the model 4 piston was defined as having chamfers on the top, bottom, and in the middle, the model 5 piston was defined as having a radius on the top, bottom and in the middle, the model 6 piston was defined as having both ends chamfered, and the model 7 piston was defined as having a radius on every edge. (a) Model 1 (b) Model 2 (c) Model (d) Model 4 (e) Model 5 (f) Model 6 (g) Model 7 Figure : Simulation models with different piston shapes Figure 4 shows the results of the magnetic flux densities between the different models when the current varied from 0.1A to 1.0A. Here the magnetic flux densities decreased with the number of chamfers on the damper piston,

the maximum density of magnetic flux appeared in Model 1, which means that this piston had the optimal geometrical shape. Moreover, density of magnetic flux in the model with a radius was greater than the chamfered model; the reason can be found in Figure 5 which shows that the distribution of magnetic flux lines along the resistance gap in model 1 were more even that in models 6 and 7. The damper piston with chamfers had larger reluctances at its core because is the core was larger and the cross-sectional area through which the magnetic flux passed decreased, which in turn caused the magnetic flux densities in the fluid resistance gap to decrease as well. 0.6 0.4 Model 1 Model 2 Model Model 4 Model 5 Model 6 Model 7 B/T 0.2 0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 Figure 4: Magnetic flux density of different models I/A (a) Model 1 (b) Model 6 (c) Model 7 Figure 5: Distributions of magnetic flux lines of different models D 12 9 6 Model 1 Model 2 Model Model 4 Model 5 Model 6 Model 7 0 0.0 0.2 0.4 0.6 0.8 1.0 I/A Figure 6: Dynamic range of MR damper under different currents Figure 6 shows the dynamic range of the proposed MR damper under different piston configurations when the applied current varied from 0.1A to 1.0A. Here the dynamic range decreased as the number of chamfers on the pistons increased, and so too did the damping performance. Moreover, the dynamic range in the model with a radius on the edges was greater than the model with chamfered edges. When the current varied from 0.1A to 1.0A, the damping force steadily increased and then stabilised when the current was close to 1.0A due to the magnetic saturation of MRF in the resistance gap. From Figure 4 to Figure 6, model 1 with square ends had the maximum B value under different currents, while the damping force and dynamic range was also better in model 1 than in the 4

other models. So the piston with square ends was selected as the optimal geometry for the double coil MR damper. To optimise the double coil MR damper using FEM, an analysis log file to solve the damper s magnetic circuit calculate the objective function was built using the ANSYS parametric design language (APDL). In the analytical log file, the length of resistance L 1 and the radius of the damper core r 1 were used as the design variables (DVs), and initial values were assigned to them, respectively. First, starting with the initial values of the DVs, the magnetic flux density, damping force and dynamic range were calculated by executing the log file. The ANSYS optimisation tool then transformed the optimisation problem with constrained DVs to an unconstrained one via penalty functions. The search direction of DVs was assumed to be the negative of the gradient of unconstrained objective function, and a combination of a golden-section algorithm and local quadratic fitting techniques were used to calculate new DVs. If convergence occurs, the DVs at this iteration would be the optimum, but if not then subsequent iterations would be carried out. Figure 7 shows the results of the damping force between the initial design and optimal design, respectively. Here the damping force with an optimal design was greater than the initial design, although the difference between both of them improved significantly when the applied current exceeded 0.5A. 5 4 Optimal Initial F/KN 2 1 0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 I/A Figure 7: Comparison of damping force between initial design and optimal design 4. Experimental evaluation of double coil MR damper Figure 8 shows the change in the damping force under different damper displacements and applied current directions. Here, the INSTRON 8801 test machine was set at sinusoidal loading with a frequency of 0.5Hz and displacement amplitude of 5mm and 7.5mm, respectively. The applied current I 1 and I 2 in the double coils were set to 0.5A with the same direction and reverse direction, respectively. The figure shows that the damping force with the currents in a reverse direction were much greater that currents in the same direction. Moreover, the damping force increased as the displacement amplitude increased, because the increased velocity of the damper led to an increase in the viscous force in Eq.(2), so the damping force also increased. 0.6 0.4 S=5.0mm I 1 =I 2 S=7.5mm I 1 =I 2 S=5.0mm I 1 =-I 2 S=7.5mm I 1 =-I 2 0.2 F/KN 0.0-0.2-0.4-0.6-8 -6-4 -2 0 2 4 6 8 Figure 8: Damping forces under different displacements and different current directions Figure 9 shows the damping force under the applied currents changed from 0A to 1.0A. Here, the INSTRON 8801 test machine was set at a sinusoidal loading with displacement amplitude of 10mm and a frequency of 1.0Hz, while the applied currents I 1 and I 2 in the double coils were applied the reverse direction. As expected, the damping S/mm 5

force increased from 0.KN at 0A to 1.21KN at 1.0A, and the dynamic range nearly equal 4. 1.5 0A 0.25A 0.5A 0.75A 1.0A 1.0 0.5 F/KN 0.0-0.5-1.0-1.5-15 -10-5 0 5 10 15 S/mm Figure 9: Damping force under different current 5. Conclusion A double coil MR damper was proposed and the damping force and dynamic range were also derived. A finite element model was built to investigate the performance of a double coil MR damper while considering pistons with different configurations; the simulation results showed that the piston with square ends had a better damping performance. An optimisation procedure was constructed using the ANSYS APDL to obtain the optimal value of the double coil MR damper. Further, an experimental analysis was also carried out to verify the damping performance of the proposed MR damper. 6. Acknowledgements This research was financially supported by the National Natural Science Foundation of China (No. 51475165, 51165005, 11462004) and the Australian Research Council Discovery Project (No. 150100266). 7. References [1] H. H. Zhang, C. R. Liao and W. M. Chen, A magnetic design method of MR fluid dampers and FEM analysis on magnetic saturation, Journal of intelligent material systems and structures, 17(8-9), 81-818, 2006. [2] M. S. A. Khan, A. Suresh, N. S. Ramaiah, Analysis of magnetorheological fluid damper with various piston profiles. International Journal of Engineering and Advanced Technology, 2(2), 77-8, 2012. [] M. M. Ferdaus, M. M. Rashid and M. H. Hasan, Optimal design of Magneto-Rheological damper comparing different configurations by finite element analysis, Journal of mechanical science and technology, 28(9), 667-677, 2014. [4] Z. Parlak, T. Engin and V. Ari. Geometrical optimisation of vehicle shock dampers with magnetorheological fluid, International journal of vehicle design, 54(4), 71-92, 2010. [5] Z. Parlak, T. Engin and İ Çallı, Optimal design of MR damper via finite element analyses of fluid dynamic and magnetic field, Mechatronics, 22(6), 890-90, 2012. [6] Q. H. Nguyen and S. B. Choi, Optimal design of a vehicle magnetorheological damper considering the damping force and dynamic range, Smart Materials and Structures, 18(1), 01501, 2009. [7] Q. H. Nguyen, S. B. Choi and J. K. Woo, Optimal design of magnetorheological fluid-based dampers for front-loaded washing machines, Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 228(2), 294-06, 2014. 6