基于高保真磁流变液建模的智能移动性提升与自适应减振器控制开发

摘要

本方案详细介绍了温度补偿型磁流变（MR）阻尼器控制系统的设计方法，该系统结合了高保真度磁流变液建模、优化的流体制备工艺以及自适应热补偿控制算法。该方法可应用于汽车悬架系统，显著提升电动汽车的乘坐舒适性。

摘要

随着高性能电动汽车（EV）的发展，需要先进的悬架系统以实现高精度、高输出力和快速适应能力，磁流变（MR）减振器在实现快速响应悬架控制方面至关重要。然而，其对温度的敏感性限制了在极端条件下的可靠性。本实验方案提出了一种系统性方法以应对这一挑战。通过将羰基铁颗粒分散于含抗磨剂和抗氧化剂的热稳定载液中，合成一种高性能磁流变液。我们提出一种指数线性混合分析（ELMA）模型及其参数辨识方法，可作为双塑性宾汉模型的更优替代方案。该ELMA框架被拓展应用于磁流变减振器，并结合温度补偿算法，使电流跟踪精度提高3.98%，力输出跟踪精度提高7.75%（峰值提升达19.92%）。联合CarSim/Simulink仿真结果表明，在D级路面上，采用温度补偿的天棚控制（Sky-hook）与混合SH-ADD控制算法可使车身垂向加速度方差降低11.97%，俯仰角速度峰值降低41.78%。本方案实现了磁流变流体物理特性与自适应减振器控制之间的有效衔接，为提升电动汽车在极端热环境下的悬架系统性能提供了可复现的技术流程。

引言

近年来，高性能电动汽车（EV）的快速发展为人类社会带来了更多便利，同时也为传统车辆带来了新的挑战与机遇。例如，它们要求更灵敏、更精确的动力系统¹，以及能够提供高精度、高输出力和快速适应能力的先进悬架系统。磁流变（MR）减振器非常适合满足这些要求，其利用磁流变液——一种智能材料，在外加磁场作用下可迅速且可逆地改变其流变特性（主要是黏度和屈服应力）^2,3,4,5,6。这种可控的变化使磁流变减振器能够精确调节阻尼力。然而，充分发挥其潜力需要精确的建模，尤其是在汽车应用典型工况下（例如，温度范围可从 -30 °C 到 +80 °C）保持性能的一致性和可预测性，仍是当前面临的主要挑战。

现有针对磁流变阻尼器（MRD）的研究在通过磁流变液（MRF）流变学数据分析系统表征其动态行为方面存在显著局限性^{4,6,7,8,9,10,11,12,13,14,15,16,17,18,19}。尽管已有大量参数化建模研究报道，但这些模型主要服务于控制算法的开发，而非实际的工程设计应用^17,18,19。此外，当前的建模范式常采用过度简化的假设，限制了其在不同工况和装置结构下的普适性^{20,21,22,23,24,25,26,27}。

因此，我们提出一种综合方法，旨在通过三项关键技术突破来克服上述实际应用中的局限性：优化磁流变液配方，以提升其在典型汽车工作温度范围内的热稳定性；提出一种新的指数线性混合分析（ELMA）模型及参数辨识方法，能够准确捕捉不同工况下的系统行为；并开发用于实时磁流变阻尼器（MRD）控制的温度补偿算法。该方法的目标是提供一种能够确保MRD可靠运行的解决方案，特别针对温度变化问题。该方法通过统一的流体稳定性优化、高精度的先进建模以及具备温度补偿能力的控制策略，带来了显著优势。流变学实验验证了ELMA模型的准确性，而车辆动力学仿真则证明了温度补偿策略在不同温度下改善悬架性能的有效性，从而提升了磁流变阻尼器在电动汽车等应用中的可行性。

方案

1. 自定义低密度磁流变液的合成

1. 准备以下材料：基础油 PAO162 以及添加剂：有机膨润土、碳酸丙烯酯、表面活性剂和羰基铁粉。
2. 配制混合物的步骤如下：
   1. 将 PAO162（用于 MRF1）放入一个 500 mL 的不锈钢容器中。将容器浸入 50 °C 的水浴中，平衡 10 分钟。
   2. 向基础油中加入有机膨润土，使用机械搅拌器（例如 IKA RW 20）以 600 rpm 的速度搅拌 20 分钟。加入碳酸丙烯酯（PC），以 600 rpm 搅拌 10 分钟。再加入表面活性剂 T154A，以 600 rpm 搅拌 5 分钟。
   3. 逐步将 CIP-SQ 加入混合物中。将搅拌速度提高至 1,500 rpm，持续混合 1 小时。在整个混合过程中监测温度，保持在 50 ± 2 °C。实验装置如图 1所示。

2. ELMA 模型的建立

1. 对于 ELMA 模型的构建，定义 ELMA 方程表达式为：
   
   ​其中： = 剪切应力（Pa）， = 动态屈服应力（Pa）， = 塑性黏度（Pa·s）， = ELMA 应力，用于量化非线性偏离程度（Pa）， = 指数衰减常数（s），与临界剪切速率（）成反比。
2. 定义模型的关键特征。在高剪切速率（）下，将 ELMA 简化为宾汉模型（）。在低剪切速率下，定义非线性项（），以捕捉对线性行为的偏离。将静态屈服应力定义为 。
3. 对于数据准备和加权，使用剪切速率测试协议。
   1. 使用流变仪进行对数剪切速率扫描（）。每个数量级采集 20 个数据点。
   2. 采用加权策略，并通过添加一个外推值（）扩展剪切速率列表。
   3. 使用以下方程计算剪切速率之间的距离：
      
   4. 使用以下方程对权重进行归一化：
      
      其中 为平均距离。
4. 使用以下参数识别算法（图 2）。
   1. 为实现参数空间的线性化，将 ELMA 方程重写为矩阵形式：
      
      其中： = 测量剪切应力的列向量， = 通过非线性算子对剪切速率数据映射生成的矩阵， = 参数矩阵 。
   2. 使用加权最小二乘法求解 ，所用方程如下：
      
      其中 为权重 的对角矩阵。
   3. 对参数 b 进行非线性优化。将 b 初始化为较小值（例如 ）。使用 MATLAB 的 lsqnonlin 求解器最小化残差范数：
      
      ​迭代直至收敛（容差：）或达到最大迭代次数（例如 1000 次）。
5. 使用以下验证指标。
   1. 对于误差量化，使用平方误差和（SSE）方程、均方根误差（RMSE）方程和决定系数（R²）方程，如下所示：
      
      
      
6. 按以下方式实施工作流程。
   1. 对于数据采集，在受控温度（例如 40 °C、70 °C）和磁场（0–1.07 T）条件下进行流变测试。将剪切应力（）和剪切速率（）数据导出至 MATLAB。
   2. 对于加权和线性化，应用步骤 2.3 中描述的加权方案。使用初始 构建 矩阵。通过加权最小二乘法求解 并使用非线性最小二乘法优化 b，完成参数优化。
   3. 对于验证，将模型预测结果与实验数据进行比较。计算 RMSE 以评估准确性。

                          

3. 温度补偿型阻尼器控制的实现

1. 温度依赖性ELMA模型开发
   1. 为实验表征热效应，建立减振器测试装置。将定制的磁流变减振器（行程：75 mm，活塞直径：50 mm）安装在电液伺服试验机上。在减振器表面安装热电偶以实现温度的实时监测。将减振器置于温控箱内（温度范围：-30 °C 至 80 °C）。
   2. 进行热循环。在每个温度点（-30 °C、-20 °C、-10 °C、…、80 °C）下，施加振幅为 40 mm、频率为 2.085 Hz 的正弦激励。在激励电流分别为 0 A、1 A、2 A、3 A 的条件下，测量阻尼力。每次达到设定温度后，等待 60 分钟以确保温度稳定，然后开始数据采集。实验装置如图所示 图3.
   3. 对于具有电流依赖性的ELMA模型参数化，使用以下方程：
      
      其中：  （动态屈服应力）；  （屈服后阻尼系数）；   非线性偏差幅度  （指数衰减常数）
   4. 对于二次拟合，拟合参数 ，等等，通过最小二乘法优化拟合实验数据。
2. 温度补偿因子设计
   1. 引入一个与温度相关的指数修正因子，  如下所示
      
      何处  25 °C 下的屈服后阻尼系数  经验常数。
   2. 为确定修正后的屈服后阻尼系数，需分别利用以下两个方程对温度效应进行阻尼系数修正，并更新包含热修正项的ELMA模型。
      
      
      方程的展开形式为
      
3. 实时电流控制的逆模型（图4).
   1. 为构建简化的宾汉流体类似模型，忽略实时控制中的非线性偏差项，如下所示。
      
      修饰后的扩展形式为
      
   2. 对于电流计算，求解该二次方程 :
      
      系数如下：
      
      
      
   3. 对于电流钳制，施加约束条件以 .
4. 模拟与验证
   1. 用于对照表面生成，使用以下输入进行电流跟踪：期望力（），活塞速度（），温度（）和输出：电流指令（）。使用以下输入进行力追踪：电流指令（），活塞速度（），温度（）以及输出：MRD 力（).
   2. 针对验证指标，进行电流跟踪精度测试，并计算相对于未补偿模型的提升程度（例如，在80 °C时平均提升3.98%）。进行力跟踪精度测试，并计算相对于未补偿模型的提升程度（例如，在80 °C时平均提升7.75%）。
   3. 评估车辆悬架指示器温度补偿效果。使用典型的四分之一悬架图5用于模拟相关性分析。在D级路面上，带补偿控制的天棚阻尼（Sky-hook）组能量值从1.2504降低至1.1007，带补偿控制的混合天棚-自适应阻尼分配（Mixed SH-ADD）组能量值从1.0565降低至0.9428。补偿后，天棚阻尼组的俯仰角速度峰值绝对值从6.2775降至4.0287，混合天棚-自适应阻尼分配组则从6.2776降至3.6551。
5. 实施流程
   1. 进行热表征时，在不同温度条件下（-30 °C 至 80 °C）开展阻尼器测试。导出力-速度数据用于参数拟合。
   2. 用于模型参数化，拟合 ELMA 模型参数（,等等使用 MATLAB 的 lsqcurvefit 确定  常数（通过非线性回归。
   3. 对于逆模型部署，使用饱和模块在 Simulink 中实现模型以进行电流钳位。
   4. 为验证，检查追踪的准确性（图6 和 图7）。比较经补偿与未经补偿的车辆悬架系统的响应（图8 和 图9).

结果

本方案中提出的温度补偿策略主要针对磁流变阻尼器，补偿控制是基于在不同电流、速度和温度条件下对阻尼器进行的综合测试结果确定的。此外，我们也开展了与磁流变液相关的温度实验，获得了如图10所示的结果。当温度高于0 °C时，可以看出零磁场屈服应力和塑性黏度随温度的变化规律基本一致，且零磁场屈服应力对黏度的影响较小。磁流变液的零磁场屈服应力受多种因素影响，如触变剂网络结构和固体颗粒等，其温度效应的机理和建模较为复杂，因此本文仅研究宽温范围内温度对磁流变液塑性黏度的影响。通过本方案获得的结果验证了所提出的指数线性混合分析（ELMA）模型辨识方法及磁流变（MR）阻尼器系统温度补偿控制策略的有效性。

图2 展示了无约束参数辨识方案的流程。如图4所示，利用已验证的ELMA模型参数，可为磁流变（...

讨论

该方案通过三项相互关联的创新，系统性地解决了磁流变（MR）阻尼器的温度敏感性问题：优化的磁流变液合成、指数线性混合分析（ELMA）建模以及自适应温度补偿算法。关键步骤包括在受控温度（50 ± 2 °C）下通过逐步添加添加剂实现磁流变液的梯度均质化，以最小化颗粒沉积。通过非线性映射算子对原始参数空间进行线性化处理，降低了待优化参数空间的维度，减少了参数辨识陷入局部最优的可能性，仅需通过非线性最小二乘法确定使代价函数最小化的非线性参数 b，即可获得全局最优参数集。提出了一种指数型热校正因子用于补偿控制，实现了精确的电流跟踪与力跟踪，并进一步优化了车辆悬架的舒适性指数。

尽管本方案为建模与温度补偿提供了稳健的框架，用户在使用过程中仍可能遇到困难或希望对方法进行调整。在ELMA模型参数辨识过程中，一个常见问题是参数初值的敏感性或实验数据中的噪声干扰¹⁷。若基于最小二乘法的方法无法可靠收敛，可考虑在参数辨识前实施数据滤波技术，或根据函数特性选择合适的初始值范围。准确的温度测量也至关重要；用户应仔细校准温度传感器，并考虑传感器的...

材料

本文使用的材料清单

姓名	公司	目录编号	评论
空气压缩机	Dynair	DA5002CS
羰基铁粉	BASF	CIP-SQ
电液伺服疲劳试验机	DOCER	PWS-1000
MATLAB	Mathworks	R2022b
有机膨润土	ELEMENTIS	-
聚α-烯烃	INEOS	PAO162
碳酸丙烯酯	MACKLIN	PC
流变仪	Anton Paar	MCR702
表面活性剂	Jinzhou Xinxing	T154A
温控平台	Peltier	H-PTD200