# Development of an Experimental Model for a Magnetorheological Damper Using Artificial Neural Networks (Levenberg-Marquardt Algorithm)

Authors: Ayush Raizada; Pravin Singru; Vishnuvardhan Krishnakumar; Varun Raj
Journal: Advances in Acoustics and Vibration
 (2016)

---

## Abstract

This paper is based on the experimental study for design and control of vibrations in automotive vehicles. The objective of this paper is to develop a model for the highly nonlinear magnetorheological (MR) damper to maximize passenger comfort in an automotive vehicle. The behavior of the MR damper is studied under different loading conditions and current values in the system. The input and output parameters of the system are used as a training data to develop a suitable model using Artificial Neural Networks. To generate the training data, a test rig similar to a quarter car model was fabricated to load the MR damper with a mechanical shaker to excite it externally. With the help of the test rig the input and output parameter data points are acquired by measuring the acceleration and force of the system at different points with the help of an impedance head and accelerometers. The model is validated by measuring the error for the testing and validation data points. The output of the model is the optimum current that is supplied to the MR damper, using a controller, to increase the passenger comfort by minimizing the amplitude of vibrations transmitted to the passenger. Besides using this model for cars, bikes, and other automotive vehicles it can also be modified by retraining the algorithm and used for civil structures to make them earthquake resistant.


---

## Body

## 1. Introduction

Isolation of the forces transmitted by external application is the most important function of a suspension system. The suspension system comprises a spring element and a dissipative element, which when placed between the object to be protected and the excitation reduces the vibration transmitted to the object.Suspension systems range from active to passive suspensions. MR damper lies in between this range and behaves like a semiactive suspension system. The damping of a passive suspension system is a property of the system and cannot be varied, whereas in an active suspension system the damping of the system can be altered by using an actuator to give it an external force. This external force helps in improving the ride quality. The shortcoming of this model is that it requires considerable amount of external power source; it is costly and is difficult to incorporate in the system due to the added mass. A variation of the active suspension system is the semiactive or the adaptive suspension system. In these systems the damping of the system is varied by controlling the current thereby changing the viscous properties of the damping elements in the suspension system. PID neural network controller is one such controller used to develop a model to predict the displacement and velocity behavior of the MR damper [1]. In comparison to active suspension semiactive suspension systems’ power consumption is considerably less.The magnetorheological (MR) dampers and electrorheological dampers are the most common examples of semiactive dampers. In the past few years, research on MR damper has improved its capabilities and reduced the gap between adaptive suspension system and truly active suspension systems. Structurally, the MR damper is similar to a simple fluid damper, except that the viscosity of the fluid in the MR damper can be changed by altering the current in the system that induces a magnetic field. The MR fluid is a non-Newtonian fluid composed of mineral oil with suspended iron nanoparticles. When there is no magnetic flux (zero current), the MR damper behaves like a normal fluid damper in which the iron nanoparticles are randomly oriented as seen in Figure1(a).Figure 1
Alignment of iron particles at (a) zero current; (b) low current; (c) high current.
(a)
(b)
(c)When a magnetic field is applied to the MR damper, the iron nanoparticles align themselves along the magnetic flux, as seen in Figures1(b) and 1(c), and form chains which makes the fluid partly semisolid. These particles help in reinforcing the damping by forming chains in the oil that obstructs the movement of the oil as the magnetic field developed is in the direction perpendicular to the movement of the oil. Hence increasing the magnetic field increases damping in the system.The behavior of the MR damper can be characterized by its highly nonlinear motion [2–6]. The force versus velocity plot forms a hysteresis loop depicting the nonlinear character of the MR damper. There exist many different parametric models like Bingham model, Bingham Body model, Lee model, Spencer model, Bouc-Wen model, and Gamota-Filisko model that portray the behavior of the MR damper. These models are difficult to implement throughout the working range of the MR damper due to the hysteresis and jump type phenomena resulting from the specific properties, like viscoplasticity, viscoelastoplasticity, and viscoelasticity, of the MR fluid. It is also difficult to find a solution for the equations for some of the defined models due to the numerical and analytical complexity of the model equation. These drawbacks prevent the models from being implemented as governing equations in real time controllers [7]. In the low speed range the Bingham model can predict the rigid plastic behavior better than the involution model when the MR damper is loaded with a sinusoidal input. When triangular loading is used on the MR damper, the involution model predicts the model better than the Bingham model [8]. The deformation in the hysteresis loop of the force-velocity and force-displacement graphs deviates from the Bouc-Wen model due to the force lag phenomenon in the MR damper. The modification of the Bouc-Wen model, the Bouc-Wen-Baber-Noori model, can to a certain extent describe the pinching hysteretic behavior [9]. The evolutionary variable equation for the modified Bouc-Wen model depends on 4 parameters A, β, γ, and n, where A, β, and γ control the shape and size of the hysteresis and the parameter n controls the smoothness of the transition from elastic to plastic region [10]. Genetic algorithm assisted inverse method and nonlinear-least square error optimization in MATLAB can be used to identify the parameters and develop a model [11, 12].


## 2. Artificial Neural Networks

Artificial Neural Network is an important tool used to build models for dynamic systems. It includes a variety of modeling tools and models. This method has considerable advantages over other modeling tools. ANNs can work easily with nonlinearly separable data, and this makes them ideal for applications such as machine condition monitoring, where the training data are sparse, and the network will have to generalize well. Several applications have demonstrated that a neural network can successfully recognize and classify different faults in a number of different condition monitoring applications [13]. A good general introduction to neural networks is provided by Haykin [14] and also Rojas [15]. Multilayer perceptron (MLP) ANNs have been used for classification purposes in this experiment. The MLP used in this case consists of one hidden layer and an output layer, the hidden layer having a logistic activation function (see (1)), while the output layer uses a linear activation function (see (2)):(1)
φ
v
=
1
1
+
exp
⁡
-
v
(2)
φ
v
=
v
,where v is the sum of the weighted outputs. Different sizes of the hidden layer were tested. The size of the output layer was set at two neurons for this particular application. Training of the ANN networks is carried out using the Levenberg-Marquardt algorithm discussed in Section 3. The network is trained using a training and validation set, while testing is carried out using a further test feature set. Performance is measured in terms of the network’s classification success on unseen data in the test set. Training is stopped when the classification performance of the validation set starts to diverge from that of the training set (i.e., the network becomes overtrained on the training set) [16].The neural network architecture consists of fully connected feedforward network consisting of 100 hidden layers as shown in Figure2. The optimization is achieved using the network growing approach where one hidden neuron is used initially and is increased until there is improvement in the performance. This optimization is carried out using the fitting tool in the neural network toolbox in Matlab.Figure 2
Neural Network Architecture.


## 3. Levenberg-Marquardt Algorithm

A mathematical description of the Levenberg-Marquardt (LM) neural network training algorithm has been presented by Hagan and Menhaj [17]. The LM algorithm was originally created [18, 19] to address the drawbacks of the Gauss-Newton (GN) method and gradient descent algorithm. The quadratic convergence properties of the GN algorithm make it very fast. However, the choice of initial weights is the key to success of this method and since these values are not readily available in real life problems, the GN method becomes slightly redundant in these cases. Since, in real-world problems, the prediction of an appropriate set of initial values is not always possible, the GN method is impractical for many applications. In comparison to the GN algorithm, the success of the gradient descent algorithm is less dependent on the initial choice of weight values. However, since the gradient descent algorithm approaches the minimum in a linear (first order) manner, its speed of convergence is normally low and, thus, does not always possess adequate convergence properties.By merging the advantages of GN and gradient descent algorithms, a hybrid algorithm (LM algorithm) was created to cater to the drawbacks of GN and gradient descent algorithms. The LM algorithm possesses quadratic convergence (approximates the GN method) when it is in the vicinity of (but not too close to) a minimum [20]. LM uses gradient descent to improve on an initial guess for its parameters and transforms to the GN method as it approaches the minimum value of the cost function. Once it approaches the minimum, it transforms back to the gradient descent algorithm to improve the accuracy. The LM algorithm is used for curve fitting [13–15] and many other optimization problems [21–24]. Due to its desirable convergence capabilities, in many optimization applications, the LM method is usually preferred over many other optimization techniques [25–29].


## 4. Experimental Setup

The experimental setup consists of three components: external actuation system, data acquisition system, and the controller. The external actuation for the system was fabricated to excite the damper with the help of the electrodynamic shaker. The system was designed to closely depict a quarter car model. The shaker transmits the road disturbances to the lower plate by means of a stringer as seen in Figure3. The lower plate is equivalent to the tire of the vehicle in the quarter car model. The lower plate is coupled with the upper plate with a LORD manufactured MR damper RD-8040-1 [30] which is a part of the suspension system. The upper plate behaves as the chassis for vehicle. The MR damper is clamped to the upper and lower plates by means of an L clamp.Figure 3
Experimental setup of MR damper.The electrodynamic shaker (Model 2110E) [31] with a load capacity of 489 N (sine-peak) is connected to the linear power amplifier (Model 2050E09) [32] which provides the shaker with the input profile that is given to the system using the Spider 81 data acquisition system [33] and the EDM software. The vibration signals, taken as inputs for genetic programming, are measured using sensors. The signal is measured on the shaker and the lower plate using accelerometers “352C34 LW 155857” of sensitivity 101.3 mV/g [34] and “352C68 SN 92017” of sensitivity 102.2 mV/g [35], respectively. To measure the force transmitted and the acceleration signals at the upper plate an impedance head “288D01 SN 3176” is used of sensitivity 98.73 mV/LBF for force and 101.3 mV/g for acceleration [36]. The MR damper is operated using an external power source. The current in the damper is varied using a potentiometer type device called Wonder Box RD-3002-03 [37]. The current is measured using a multimeter.The signal parameters from these sensors are used as inputs for the algorithm and the current that is varied using the potentiometer is used as the output. This current can be directly provided to the control system to create a feedback loop.


## 5. Experimentation

The training set for the ANN is obtained from experimentation using the setup described in Section3. In the experimental setup system there are 6 factors that are monitored to develop the training set. The experimentation is carried out by varying the frequency and current parameters. The data collected is divided into three sets; that is, 70% is training set, 15% is testing set, and the remaining 15% is validation. Using this to training set the final model is developed which best fits the data having minimum error. Table 1 provides the configuration of various experiments conducted by varying the parameters.Table 1
Different test configurations.
Serial number
Description
Frequency
Steps
Current
1
Dwell
5–7 Hz
0.1 Hz
0 A
2
Dwell
5–7 Hz
0.1 Hz
0.25 A
3
Dwell
5–7 Hz
0.1 Hz
0.50 A
4
Dwell
5–7 Hz
0.1 Hz
0.75 A
5
Dwell
5–7 Hz
0.1 Hz
1.0 AThe input parameters are recorded for 10 seconds of the response. These signatures are as shown in Figure4.Figure 4
(a) Signatures of acceleration of road profile, (b) signatures of acceleration of lower plate, (c) signatures of acceleration of upper plate, and (d) signatures of force transmitted to upper plate.
(a)
(b)
(c)
(d)In the first set of tests the current is kept constant and the frequency is increased from 5 Hz to 7 Hz in steps of 0.1 Hz. The data is collected at a sampling rate of 20.48 kHz. In the similar fashion the data set is obtained for all the current values. This data is then consolidated to form the training and testing matrix. There are 4 input parameters taken into consideration in this system:(1)
Frequency of excitation.
(2)
Acceleration of the shaker.
(3)
Relative acceleration of the upper level (chassis) with respect to the lower level (tire).
(4)
Force transmitted to chassis.The current in the damper coils is taken as the output parameter. Using these input-output parameters the algorithm is trained.


## 6. Results and Discussion

The best model selected from the models generated in the different stages is based on the performance of the model. The algorithm is trained using the Levenberg-Marquardt algorithm. The error histogram for the model obtained with 100 nodes in the hidden layer is shown in Figure5. It can be seen from the histogram that the large centre peaks indicate very small errors or the output is very close to the target values. The absence of spikes towards the end implies that the there are no misclassifications in the model [38].Figure 5
Error in the best model.The best training performance for the different iterations of the model is shown in Figure6.Figure 6
Best training performance.The training state for each iteration for the various parameters is shown in Figure7.Figure 7
Training state of various parameters.Figure8 shows the vibration signature for the force exerted on the upper mass of the fabricated structure. From the figure it is visible that there is a massive reduction in the amplitude of the force exerted due to the optimization performed.Figure 8
(a) Vibration signature without optimization; (b) signature with optimization.
(a)
(b)


## 7. Conclusion

This work is the first of its kind that develops a model to predict the current value to adjust the damper force for the given input parameters, using actual experimental data, in order to maintain the driver acceleration. This proves the versatility and success of Artificial Neural Networks in highly nonlinear and hysteretic systems like MR damper. The optimization performed, based on the Artificial Neural Network, has reduced the amplitude of force transmitted to the driver. The magnitude of the force transmitted is approximately 1/10 that at the base of the vehicle. The current parameter determined based on this optimization has caused this reduction in the force transmitted to the passenger seated in the vehicle.



---

*Source: 7027259-2016-08-18.xml*