Smart Mobility Boosting Using High-Fidelity Magnetorheological Fluid Modeling for Adaptive Damper Control Development

Summary

This protocol details the development of a temperature-compensated magnetorheological (MR) damper control system that combines high-fidelity magnetorheological fluid modeling, optimized fluid preparation, and adaptive thermal compensation control algorithms. The methodology can be applied to automotive suspensions to significantly enhance the comfort of electric vehicles.

Abstract

With the rise of high-performance electric vehicles (EVs) which require advanced suspension systems capable of delivering precision, high force output, and rapid adaptability, magnetorheological (MR) dampers are critical for achieving rapid-response suspension control. However, their temperature sensitivity limits reliability under extreme conditions. This protocol presents a systematic approach to address this challenge. A high-performance MR fluid is synthesized using carbonyl iron particles dispersed in a thermally stable carrier fluid with anti-wear agents and antioxidants. We propose an Exponential Linear Mixed Analysis (ELMA) model and its parameter identification method, which can be considered as a superior alternative to the bi-plastic Bingham model. The ELMA framework is extended to MR dampers, with temperature compensation algorithms improving current tracking accuracy by 3.98% and force tracking accuracy by 7.75% (peak: 19.92%). Joint CarSim/Simulink simulations demonstrate that temperature-compensated Sky-hook and Mixed SH-ADD algorithms reduce vertical acceleration variance by 11.97% and peak pitch rate by 41.78% on Class D roads. This protocol bridges MR fluid physics to adaptive damper control, offering a replicable workflow for enhancing EV suspension systems in extreme thermal environments.

Introduction

The rapid development of high-performance electric vehicles (EVs) in recent years has brought more convenience to human society as well as new challenges and opportunities for conventional vehicles. For example, they require more responsive and precise powertrains¹ and advanced suspension systems capable of delivering precision, high-force output, and quick adaptability. Magnetorheological (MR) dampers are ideally suited to meet these requirements, utilizing magnetorheological fluids - a smart material that rapidly and reversibly changes its rheological properties (primarily viscosity and yield stress) in the presence of a magnetic field^2,3,4,5,6. This controllable change enables the magnetoresistive damper to precisely regulate the damping force. However, exploiting this potential requires accurate modeling, especially since maintaining consistent, predictable performance under the typical operating conditions of automotive applications (e.g., the temperature range can be from -30 °C to +80 °C) remains a major challenge.

Existing research on magnetorheological dampers (MRDs) exhibits notable limitations in systematically characterizing dynamic behaviors through magnetorheological fluid (MRF) rheological data analysis^{4,6,7,8,9,10,11,12,13,14,15,16,17,18,19}. While substantial parametric modeling efforts have been documented, these predominantly serve control algorithm development rather than practical engineering design applications^17,18,19. Furthermore, current modeling paradigms frequently employ oversimplified assumptions that restrict their generalizability across diverse operational conditions and device configurations^{20,21,22,23,24,25,26,27}.

Hence, we present an integrated approach aimed at overcoming these practical limitations through three key innovations for demanding applications: optimizing the MR fluid formulation to improve thermal stability over typical automotive temperature ranges, proposing a new Exponential Linear Mixing Analysis (ELMA) model and parameter identification methodology capable of capturing behaviors under different conditions, and developing temperature compensation algorithms for robust temperature compensation algorithms for real-time MRD control. The goal is to provide a methodology that enables reliable MRD operation, specifically addressing thermal variations. The approach provides advantages through unified fluid stability optimization, advanced modeling for improved accuracy, and temperature-compensated control. Rheological experiments validate the accuracy of the ELMA model, while vehicle dynamics simulations demonstrate the effectiveness of the temperature compensation strategy in improving suspension performance at different temperatures, thus enhancing the viability of MR dampers for applications such as electric vehicles.

Protocol

1. Synthesis of custom low-density MR fluids

1. Obtain the following materials: base oil PAO162 and additives: Organic bentonite, propylene carbonate, surfactant, and carbonyl iron powder.
2. To prepare the mix, follow the steps described below.
   1. Place PAO162 (for MRF1) in a 500 mL stainless steel vessel. Immerse the vessel in a 50 °C water bath and equilibrate for 10 min.
   2. Add organic bentonite to the base oil and stir at 600 rpm for 20 min using a mechanical stirrer (e.g., IKA RW 20). Add propylene carbonate (PC) and stir at 600 rpm for 10 min. Add surfactant T154A and stir at 600 rpm for 5 min.
   3. Gradually add CIP-SQ to the mixture. Increase stirring speed to 1,500 rpm and mix for 1 h. Monitor temperature to maintain 50 ± 2 °C throughout mixing. The experimental setup is shown in Figure 1.

2. ELMA model development

1. For ELMA model formulation, define the ELMA equation expressed as:
   
   ​where:  = Shear stress (Pa),  = Dynamic yield stress (Pa),  = Plastic viscosity (Pa·s),  = ELMA stress, quantifying nonlinear deviation magnitude (Pa),  = Exponential decay constant (s), inversely related to critical shear rate ().
2. Define key model features. At high shear rates (), simplify ELMA to the Bingham model (). At low shear rates, define the nonlinear term () that captures deviations from linear behavior. Define static yield stress as .
3. For data preparation and weighting, use the shear rate testing protocol.
   1. Perform logarithmic shear rate sweeps () using a rheometer. Acquire 20 data points for each order of magnitude.
   2. Use a weighting strategy and extend the shear rate list by adding an extrapolated value().
   3. Calculate inter-shear-rate distances using the following equation:
      
   4. Normalize weights using the following equation:
      
      where  is the average distance.
4. Use the following parameter identification algorithm (Figure 2).
   1. For linearization of parameter space, rewrite the ELMA equation in matrix form as follows:
      
      where:  = Column vector of measured shear stress,  = The matrix generated by mapping the shear rate data by a nonlinear operator,  = Parameter matrix .
   2. Use the weighted least squares solution to solve for  using the following equation:
      
      where  is a diagonal matrix of weights .
   3. Perform nonlinear optimization for b. Initialize b with a small value (e.g., ). Use MATLAB's lsqnonlin solver to minimize the residual norm:
      
      ​Iterate until convergence (tolerance: ) or maximum iterations (e.g., 1000).
5. Use the following validation metrics.
   1. For error quantification, use sum of squared errors (SSE) equation, Root Mean Square Error (RMSE) equation, and . R-squared (R2) equation as given below.
      
      
      
6. Implement the workflow as described below.
   1. For data acquisition, perform rheological tests under controlled temperatures (e.g., 40 °C, 70 °C) and magnetic fields (0-1.07 T). Export shear stress () and shear rate () data to MATLAB.
   2. For weighting and linearization, apply the weighing scheme described in step 2.3. Construct  matrix using initial . Perform parameter optimization by solving for  using weighted least squares and optimize b using nonlinear least squares.
   3. For validation, compare model predictions with experimental data. Calculate RMSE to assess accuracy.

                          

3. Temperature-compensated damper control implementation

1. Temperature-dependent ELMA model development
   1. For experimental characterization of thermal effects, establish a damper testing setup. Mount a custom MR damper (stroke: 75 mm, piston diameter: 50 mm) on an electro-hydraulic servo machine. Attach a thermocouple to the damper surface for real-time temperature monitoring. Enclose the damper in a temperature-controlled chamber (-30 °C to 80 °C).
   2. Perform thermal cycling. For each temperature (-30 °C, -20 °C, -10 °C,…,80 °C), apply sinusoidal excitation at amplitude: 40 mm and frequency: 2.085 Hz. Measure the damping force at excitation currents: 0 A, 1 A, 2 A, 3 A. Allow 60 min for temperature stabilization before data acquisition. The experimental setup is shown in Figure 3.
   3. For parameterization of the ELMA model with current dependency, use the following equation:
      
      where:  (dynamic yield stress);  (post-yield damping coefficient);   (nonlinear deviation magnitude);  (exponential decay constant).
   4. For quadratic fitting, fit parameters , etc, to experimental data using least squares optimization.
2. Temperature compensation factor design
   1. Introduce a temperature-dependent exponential correction factor,  as follows
      
      Where  = Post-yield damping coefficient at 25 °C,  = Empirical constants.
   2. To determine the modified post-yield damping coefficient, adjust the damping coefficient for temperature effects and update the ELMA model to include thermal correction using the following two equations, respectively.
      
      
      ​The expanded form of the equations is
      
3. Inverse model for real-time current control (Figure 4).
   1. To make a simplified Bingham-like model, neglect the nonlinear deviation term for real-time control as follows.
      
      The modified expanded form is
      
   2. For current calculation, solve the quadratic equation for :
      
      ​with coefficients:
      
      
      
   3. For current clamping, apply constraints to .
4. Simulation and Validation
   1. For control surface generation, perform current tracking with inputs: Desired force (), piston velocity (), temperature () and output: Current command (). Perform force tracking with inputs: Current command (), piston velocity (), temperature () and output: MRD force ().
   2. For validation metrics, carry out current tracking accuracy and calculate improvement relative to the uncompensated model (e.g., 3.98% average at 80 °C). Carry out force tracking accuracy and calculate improvement relative to the uncompensated model (e.g., 7.75% average at 80 °C).
   3. Assess vehicle suspension indicator temperature compensation effect. Use typical quarter suspension (Figure 5) for the simulation of the correlation. The Sky-hook group with compensation control shows a reduction in energy from 1.2504 to 1.1007 on Class D roads, and the Mixed SH-ADD group with compensation control shows a decrease in energy from 1.0565 to 0.9428 on Class D roads. After compensation, the Sky-hook group's absolute value of the peak pitch rate decreased from 6.2775 to 4.0287, and the Mixed SH-ADD group decreased from 6.2776 to 3.6551.
5. Implementation workflow
   1. For thermal characterization, conduct damper tests across temperatures (-30 °C to 80 °C). Export force-velocity data for parameter fitting.
   2. For model parameterization, fit ELMA model parameters (,etc.) using MATLAB's lsqcurvefit. Determine  constants () via nonlinear regression.
   3. For inverse model deployment, implement the model in Simulink with saturation blocks for current clamping.
   4. For validation, check tracking accuracy (Figure 6 and Figure 7). Compare the compensated and uncompensated vehicle suspension system responses (Figure 8 and Figure 9).

Results

The temperature compensation strategy presented in this protocol focuses primarily on magnetorheological dampers, and the compensation control is determined from comprehensive tests at the damper level over a range of current, velocity, and temperature dependencies. However, we also did the magnetorheological fluid-related temperature experiments to obtain the image shown in Figure 10, when the temperature is higher than 0 °C, we can see that the change rule ...

Discussion

The protocol provides a structured approach to address the temperature sensitivity of magnetorheological (MR) dampers through three interrelated innovations: optimized magnetorheological fluid synthesis, exponential linear mixing analysis (ELMA) modeling, and adaptive temperature compensation algorithms. Key steps include gradient-based homogenization of the magnetorheological fluid with sequential addition of additives at a controlled temperature (50 ± 2 °C) to minimize particle deposition. The original parame...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
air compressor	Dynair	DA5002CS
Carbonyl Iron Powder	BASF	CIP-SQ
Electro-hydraulic servo fatigue testing machine	DOCER	PWS-1000
MATLAB	Mathworks	R2022b
Organic bentonite	ELEMENTIS	-
Poly-alpha-olefin	INEOS	PAO162
Propylene carbonate	MACKLIN	PC
Rheometer	Anton Paar	MCR702
Surfactant	Jinzhou Xinxing	T154A
temperature control platform	Peltier	H-PTD200