一种描述减振器滞回特性的Bouc–Wen改进模型

工程科学学报 Chinese Journal of Engineering 一种描述减振器滞回特性的BoucWenp改进模型 赵义伟刘永强杨绍普陈祖晨 An improved Bouc-Wen model for describing hysteretic characteristics of shock absorbers ZHAO Yi-wei,LIU Yong-qiang.YANG Shao-pu,CHEN Zu-chen 引用本文： 赵义伟，刘永强，杨绍普，陈祖晨.一种描述减振器滞回特性的Bouc Wenp改进模型.工程科学学报，2020,42(10)：1352-1361. doi10.13374j.issn2095-9389.2019.10.18.001 ZHAO Yi-wei,LIU Yong-qiang,YANG Shao-pu,CHEN Zu-chen.An improved BoucWen model for describing hysteretic characteristics of shock absorbers[J].Chinese Journal of Engineering,2020,42(10):1352-1361.doi:10.13374/j.issn2095- 9389.2019.10.18.001 在线阅读View online:https::/oi.org10.13374.issn2095-9389.2019.10.18.001 您可能感兴趣的其他文章 Articles you may be interested in 考虑磁滞的铁稼磁致伸缩位移传感器输出电压模型及结构设计 Output voltage model of Fe-Ga magnetostrictive displacement sensor considering hysteresis and structural design 工程科学学报.2017,398)：1232 https:/1doi.org/10.13374.issn2095-9389.2017.08.013 改进人工鱼群算法及其在时滞系统辨识中的应用 An improved artificial fish swarm algorithm and its application on system identification with a time-delay system 工程科学学报.2017,394：619 https::/1doi.org/10.13374斩.issn2095-9389.2017.04.018 一种改进的非刚性图像配准算法 An improved non-rigid image registration approach 工程科学学报.2019,41(7)：955 https:1doi.org/10.13374j.issn2095-9389.2019.07.015 膏体流变参数影响机制及计算模型 Influence mechanism and calculation model of CPB rheological parameters 工程科学学报.2017,392：190 https:/1doi.org10.13374.issn2095-9389.2017.02.004 一种改进的人工蜂群算法—粒子蜂群算法 An improved artificial bee colony algorithm:particle bee colony 工程科学学报.2018,40(7)：871 https:loi.org10.13374j.issn2095-9389.2018.07.014 基于改进CV模型的金相图像分割 Segmentation of metallographic images based on improved CV model 工程科学学报.2017,3912：1866 https:1doi.org/10.13374j.issn2095-9389.2017.12.013

一种描述减振器滞回特性的BoucWen改进模型 赵义伟 刘永强 杨绍普 陈祖晨 An improved Bouc –Wen model for describing hysteretic characteristics of shock absorbers ZHAO Yi-wei, LIU Yong-qiang, YANG Shao-pu, CHEN Zu-chen 引用本文: 赵义伟, 刘永强, 杨绍普, 陈祖晨. 一种描述减振器滞回特性的BoucWen改进模型[J]. 工程科学学报, 2020, 42(10): 1352-1361. doi: 10.13374/j.issn2095-9389.2019.10.18.001 ZHAO Yi-wei, LIU Yong-qiang, YANG Shao-pu, CHEN Zu-chen. An improved BoucWen model for describing hysteretic characteristics of shock absorbers[J]. Chinese Journal of Engineering, 2020, 42(10): 1352-1361. doi: 10.13374/j.issn2095- 9389.2019.10.18.001 在线阅读 View online: https://doi.org/10.13374/j.issn2095-9389.2019.10.18.001 您可能感兴趣的其他文章 Articles you may be interested in 考虑磁滞的铁稼磁致伸缩位移传感器输出电压模型及结构设计 Output voltage model of Fe-Ga magnetostrictive displacement sensor considering hysteresis and structural design 工程科学学报. 2017, 39(8): 1232 https://doi.org/10.13374/j.issn2095-9389.2017.08.013 改进人工鱼群算法及其在时滞系统辨识中的应用 An improved artificial fish swarm algorithm and its application on system identification with a time-delay system 工程科学学报. 2017, 39(4): 619 https://doi.org/10.13374/j.issn2095-9389.2017.04.018 一种改进的非刚性图像配准算法 An improved non-rigid image registration approach 工程科学学报. 2019, 41(7): 955 https://doi.org/10.13374/j.issn2095-9389.2019.07.015 膏体流变参数影响机制及计算模型 Influence mechanism and calculation model of CPB rheological parameters 工程科学学报. 2017, 39(2): 190 https://doi.org/10.13374/j.issn2095-9389.2017.02.004 一种改进的人工蜂群算法——粒子蜂群算法 An improved artificial bee colony algorithm: particle bee colony 工程科学学报. 2018, 40(7): 871 https://doi.org/10.13374/j.issn2095-9389.2018.07.014 基于改进CV模型的金相图像分割 Segmentation of metallographic images based on improved CV model 工程科学学报. 2017, 39(12): 1866 https://doi.org/10.13374/j.issn2095-9389.2017.12.013

工程科学学报.第42卷.第10期：1352-1361.2020年10月 Chinese Journal of Engineering,Vol.42,No.10:1352-1361,October 2020 https://doi.org/10.13374/j.issn2095-9389.2019.10.18.001;http://cje.ustb.edu.cn 一种描述减振器滞回特性的Bouc-Wen改进模型 赵义伟2)，刘永强2,3)四，杨绍普2，)，陈祖晨2) 1)石家庄铁道大学交通运输学院，石家庄0500432)石家庄铁道大学省部共建交通工程结构力学行为与系统安全国家重点实验室，石家 庄0500433)石家庄铁道大学机械工程学院.石家庄050043 ☒通信作者，E-mail:liuyq@stdu.edu.cn 摘要Bouc-Wn模型在非识别激励工况下模拟的阻尼力与实际阻尼力误差较大，对非识别激励振幅过于敏感，针对这一 问题，提出了一种描述减振器滞回特性的改进模型.首先用Mechanical testing and simulation(MTS)疲劳试验机对磁流变减振 器进行力学性能试验，获得在多种激励幅值、频率和电流作用下的阻尼力.采用阻尼力对位移的斜率与阻尼力关系来模拟滞 回环特性曲线.根据滞回曲线特点利用二次多项式函数来表征滞回环斜率与阻尼力的关系，同时，引人关于速度的指数函数 修正项.进而对改进后的Bouc-Wn模型进行参数识别.并对其进行仿真及验证.与试验得到的阻尼力进行对比，发现在非 识别激励工况下，曲线吻合效果较好.对改进前后Bouc-Wn模型模拟的阻尼力特性曲线进行对比.结果表明：改进后模型 得到的阻尼力仿真值能够较好地模拟试验得到的各种工况下阻尼力的值，且优于Bouc-Wen模型，同时Bouc-Wen模型在非 识别激励工况下模拟阻尼力精度较差这一问题得到了改善.新模型为保证车辆悬架系统在多变工况下仿真响应的准确性打 下了基础. 关键词磁流变减振器：性能试验：改进模型：参数识别：滞回特性 分类号TG156 An improved Bouc-Wen model for describing hysteretic characteristics of shock absorbers ZHAO Yi-wei2).LIU Yong-qiang,YANG Shao-pu CHEN Zu-chen2 1)School of Transportation,Shijiazhuang Tiedao University,Shijiazhuang 050043,China 2)State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures,Shijiazhuang Tiedao University,Shijiazhuang 050043,China 3)School of Mechanical Engineering,Shijiazhuang Tiedao University,Shijiazhuang 050043,China Corresponding author,E-mail:liuyq @stdu.edu.cn ABSTRACT The error between the actual damping force and the simulated damping force obtained using the Bouc-Wen model under non-identification excitation conditions is large,and the model is too sensitive to non-identification excitation amplitude and thus features poor accuracy.To solve this sensitivity problem,an improved model describing hysteretic characteristics of shock absorbers was proposed.Firstly,the mechanical properties of a magnetorheological(MR)damper were tested to obtain the damping force under various excitation amplitudes,frequencies and currents using a mechanical testing and simulation(MTS)fatigue testing machine.The smooth hysteresis loop curve was simulated based on the relationship between the slope of the hysteresis loop and the damping force.The quadratic polynomial function was used to characterize the relationship between the slope of hysteresis loop and the damping force according to the hysteresis curve characteristics.At the same time,the revision term of the exponential function for the velocity value 收稿日期：2019-10-18 基金项目：国家自然科学基金资助项目(11790282.12072208.52072249.11802184)：河北省自然科学基金资助项目(A2020210028):河北省 研究生创新资助项目(CXZZBS2019148)

一种描述减振器滞回特性的 Bouc–Wen 改进模型 赵义伟1,2)，刘永强2,3) 苣，杨绍普2,3)，陈祖晨1,2) 1) 石家庄铁道大学交通运输学院，石家庄 050043    2) 石家庄铁道大学省部共建交通工程结构力学行为与系统安全国家重点实验室，石家 庄 050043    3) 石家庄铁道大学机械工程学院，石家庄 050043 苣通信作者，E-mail：liuyq@stdu.edu.cn 摘    要    Bouc–Wen 模型在非识别激励工况下模拟的阻尼力与实际阻尼力误差较大，对非识别激励振幅过于敏感，针对这一 问题，提出了一种描述减振器滞回特性的改进模型. 首先用 Mechanical testing and simulation(MTS) 疲劳试验机对磁流变减振 器进行力学性能试验，获得在多种激励幅值、频率和电流作用下的阻尼力. 采用阻尼力对位移的斜率与阻尼力关系来模拟滞 回环特性曲线. 根据滞回曲线特点利用二次多项式函数来表征滞回环斜率与阻尼力的关系，同时，引入关于速度的指数函数 修正项，进而对改进后的 Bouc–Wen 模型进行参数识别，并对其进行仿真及验证. 与试验得到的阻尼力进行对比，发现在非 识别激励工况下，曲线吻合效果较好. 对改进前后 Bouc–Wen 模型模拟的阻尼力特性曲线进行对比，结果表明：改进后模型 得到的阻尼力仿真值能够较好地模拟试验得到的各种工况下阻尼力的值，且优于 Bouc–Wen 模型，同时 Bouc–Wen 模型在非 识别激励工况下模拟阻尼力精度较差这一问题得到了改善. 新模型为保证车辆悬架系统在多变工况下仿真响应的准确性打 下了基础. 关键词    磁流变减振器；性能试验；改进模型；参数识别；滞回特性 分类号    TG156 An  improved  Bouc–Wen  model  for  describing  hysteretic  characteristics  of  shock absorbers ZHAO Yi-wei1,2) ，LIU Yong-qiang2,3) 苣 ，YANG Shao-pu2,3) ，CHEN Zu-chen1,2) 1) School of Transportation, Shijiazhuang Tiedao University, Shijiazhuang 050043, China 2) State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures, Shijiazhuang Tiedao University, Shijiazhuang 050043, China 3) School of Mechanical Engineering, Shijiazhuang Tiedao University, Shijiazhuang 050043, China 苣 Corresponding author, E-mail: liuyq@stdu.edu.cn ABSTRACT    The error between the actual damping force and the simulated damping force obtained using the Bouc–Wen model under non-identification  excitation  conditions  is  large,  and  the  model  is  too  sensitive  to  non-identification  excitation  amplitude  and  thus features poor accuracy. To solve this sensitivity problem, an improved model describing hysteretic characteristics of shock absorbers was proposed. Firstly, the mechanical properties of a magnetorheological (MR) damper were tested to obtain the damping force under various excitation amplitudes, frequencies and currents using a mechanical testing and simulation(MTS) fatigue testing machine. The smooth hysteresis  loop  curve  was  simulated  based  on  the  relationship  between  the  slope  of  the  hysteresis  loop  and  the  damping  force.  The quadratic  polynomial  function  was  used  to  characterize  the  relationship  between  the  slope  of  hysteresis  loop  and  the  damping  force according to the hysteresis curve characteristics. At the same time, the revision term of the exponential function for the velocity value 收稿日期: 2019−10−18 基金项目: 国家自然科学基金资助项目（11790282，12072208，52072249，11802184）；河北省自然科学基金资助项目（A2020210028）；河北省 研究生创新资助项目（CXZZBS2019148） 工程科学学报，第 42 卷，第 10 期：1352−1361，2020 年 10 月 Chinese Journal of Engineering, Vol. 42, No. 10: 1352−1361, October 2020 https://doi.org/10.13374/j.issn2095-9389.2019.10.18.001; http://cje.ustb.edu.cn

赵义伟等：一种描述减振器滞回特性的Bouc-Wen改进模型 1353· was introduced,and the parameters of the established improved Simulink model were identified.The damping forces under different working conditions were obtained from the experiment,and the new model was simulated and validated.The damping forces obtained from new model and the experiment were compared,and the curves obtained from the model agree well with the experimental results under different working conditions.Meanwhile,the improved model was compared with the Bouc-Wen model based on the characteristic curves of the damping force.The results show that the improved model can better simulate the damping force values obtained from tests under different conditions,and is superior to the Bouc-Wen model.At the same time,the problem of poor accuracy of the Bouc-Wen model under non identification excitation conditions was improved.The new model lays the foundation for ensuring the accuracy of the vehicle suspension system response under various working conditions. KEY WORDS magnetorheological dampers:performance test;improved model;parameters identification;hysteresis characteristic 当前，诸多专家学者们将精力都投入到磁流 悬架模型上，进行了半主动控制分析.上述文献 变减振器的研究当中-刀.其作用大体和传统减振 对Bouc-Wen模型都进行了不同程度的研究与应 器相同，即抑制物体振动.但是其构造和传统减振 用，但是大多数都是在特定的或者识别激励幅值 器大不相同.磁流变阻尼器作为目前较为新型的 下进行研究，对非识别激励幅值下Bouc-Wen模 减振器，有着传统减振器不能比拟的优点.它的阻 型是否能够准确描述出减振器滞回特性没有进行 尼系数可以连续调节，结构简单紧凑.磁流变减振 验证研究.文献[20]针对Bouc-Wen模型提出了 器需要很小的电流或电压就可正常工作.其中，研 一种参数识别方法，并对模型在不同激励幅值下 究减振器滞回特性的建模是基础，同时也是一个 进行验证，发现Bouc-Wen模型对激励幅值较敏 关键的环节.只有根据减振器滞回特性，并用数学 感，在非识别激励幅值情况下阻尼力精度较差，且 模型描述出来，才能进行接下来的减振器仿真设 阻尼力越大此问题越明显.这种问题会影响模型 计.因为磁流变减振器阻尼力特性呈强烈非线性 在悬架等减振系统中的仿真应用，但针对此问题 并带有滞回环，用数学模型精确简洁地描述其特 的研究较少.因此，论文对Bouc-Wen模型进行改 性是比较困难的.现在常用的减振器数学模型有 进，对其进行参数识别与仿真验证，来重点解决此 Bouc-Wen模型I-⑧、Bingham模型Io-、多项式模 问题.论文的研究目的是，改进Bouc-Wen模型， 型II-l等.Bouc-Wen模型能够很好地反映MRD 使其能够在非识别激励幅值下模拟出的阻尼力更 的动态性能，同时能较好地反映低速时的滞回情 接近于实际阻尼力值，使其更具有实用性，为其应 况，且模拟出的滞回曲线较为平滑.因此论文采用 用于车辆悬架仿真中打下基础 此模型来展开研究 Bouc-Wen模型是1997年由Bouc和Wen提 1Bouc-Wen模型及其不足 出的，是应用比较广泛的一种模型，受到国内外许 多专家学者的关注.文献[15]设计了一种磁流变 Bouc-Wen模型可以模拟多种不同类型的滞回 特性曲线，能够较准确地描述磁流变减振器在低 减振器，建立其Bouc-Wen模型，并将其应用于列 速区的非线性滞回特性.其结构示意图如图1所 车半主动悬挂系统，并与被动悬挂系统进行了对 比.文献[I6]对Bouc-Wen模型不能模拟带有力 示.其数学模型描述为 滞后现象的滞回环这一不足进行研究，提出了用 之=-y-1-B+Ar (1) Bouc-Wen的改进模型BWBN来模拟滞回环，取 F=co+ko(-x0)+az (2) 得了较好的效果.文献[I7刀针对Bouc-Wen模型 Bouc-Wen 无法表述单出杆减振器中蓄能器造成的阻尼力偏 置这一现象，提出了一种改进的Bouc-Wen模型， 并验证了其正确性.文献[I8]对Bouc-Wen模型 进行改进，并用来描述钢板装配式屈曲约束支撑 (BRB)的滞回特性，并与试验数据进行了对比，验 证了其有效性.文献[19]设计了一种应用于汽车 悬架系统上的磁流变减振器，并用Bouc-Wen模 图1Bouc-Wen模型示意图 型描述了其滞回特性，最后将其应用在四分之一 Fig.1 Bouc-Wen model schematic

was introduced, and the parameters of the established improved Simulink model were identified. The damping forces under different working conditions were obtained from the experiment, and the new model was simulated and validated. The damping forces obtained from new model and the experiment were compared, and the curves obtained from the model agree well with the experimental results under  different  working  conditions.  Meanwhile,  the  improved  model  was  compared  with  the  Bouc –Wen  model  based  on  the characteristic  curves  of  the  damping  force.  The  results  show  that  the  improved  model  can  better  simulate  the  damping  force  values obtained from tests under different conditions, and is superior to the Bouc–Wen model. At the same time, the problem of poor accuracy of the Bouc–Wen model under non identification excitation conditions was improved. The new model lays the foundation for ensuring the accuracy of the vehicle suspension system response under various working conditions. KEY WORDS    magnetorheological dampers；performance test；improved model；parameters identification；hysteresis characteristic 当前，诸多专家学者们将精力都投入到磁流 变减振器的研究当中[1−7] . 其作用大体和传统减振 器相同，即抑制物体振动. 但是其构造和传统减振 器大不相同. 磁流变阻尼器作为目前较为新型的 减振器，有着传统减振器不能比拟的优点. 它的阻 尼系数可以连续调节，结构简单紧凑. 磁流变减振 器需要很小的电流或电压就可正常工作. 其中，研 究减振器滞回特性的建模是基础，同时也是一个 关键的环节. 只有根据减振器滞回特性，并用数学 模型描述出来，才能进行接下来的减振器仿真设 计. 因为磁流变减振器阻尼力特性呈强烈非线性 并带有滞回环, 用数学模型精确简洁地描述其特 性是比较困难的. 现在常用的减振器数学模型有 Bouc–Wen 模型[8−9]、Bingham 模型[10−12]、多项式模 型[13−14] 等. Bouc–Wen 模型能够很好地反映 MRD 的动态性能，同时能较好地反映低速时的滞回情 况，且模拟出的滞回曲线较为平滑. 因此论文采用 此模型来展开研究. Bouc–Wen 模型是 1997 年由 Bouc 和 Wen 提 出的，是应用比较广泛的一种模型，受到国内外许 多专家学者的关注. 文献 [15] 设计了一种磁流变 减振器，建立其 Bouc–Wen 模型，并将其应用于列 车半主动悬挂系统，并与被动悬挂系统进行了对 比. 文献 [16] 对 Bouc–Wen 模型不能模拟带有力 滞后现象的滞回环这一不足进行研究，提出了用 Bouc–Wen 的改进模型 BWBN 来模拟滞回环，取 得了较好的效果. 文献 [17] 针对 Bouc–Wen 模型 无法表述单出杆减振器中蓄能器造成的阻尼力偏 置这一现象，提出了一种改进的 Bouc–Wen 模型， 并验证了其正确性. 文献 [18] 对 Bouc–Wen 模型 进行改进，并用来描述钢板装配式屈曲约束支撑 （BRB）的滞回特性，并与试验数据进行了对比，验 证了其有效性. 文献 [19] 设计了一种应用于汽车 悬架系统上的磁流变减振器，并用 Bouc–Wen 模 型描述了其滞回特性，最后将其应用在四分之一 悬架模型上，进行了半主动控制分析. 上述文献 对 Bouc–Wen 模型都进行了不同程度的研究与应 用，但是大多数都是在特定的或者识别激励幅值 下进行研究，对非识别激励幅值下 Bouc–Wen 模 型是否能够准确描述出减振器滞回特性没有进行 验证研究. 文献 [20] 针对 Bouc–Wen 模型提出了 一种参数识别方法，并对模型在不同激励幅值下 进行验证，发现 Bouc–Wen 模型对激励幅值较敏 感，在非识别激励幅值情况下阻尼力精度较差，且 阻尼力越大此问题越明显. 这种问题会影响模型 在悬架等减振系统中的仿真应用，但针对此问题 的研究较少. 因此，论文对 Bouc–Wen 模型进行改 进，对其进行参数识别与仿真验证，来重点解决此 问题. 论文的研究目的是，改进 Bouc–Wen 模型， 使其能够在非识别激励幅值下模拟出的阻尼力更 接近于实际阻尼力值，使其更具有实用性，为其应 用于车辆悬架仿真中打下基础. 1    Bouc–Wen 模型及其不足 Bouc–Wen 模型可以模拟多种不同类型的滞回 特性曲线，能够较准确地描述磁流变减振器在低 速区的非线性滞回特性. 其结构示意图如图 1 所 示. 其数学模型描述为 z˙ = −γ|x˙|z|z| n−1 −βx˙|z| n + Ax˙ （1） F = c0 x˙ +k0 (x− x0)+αz （2） Bouc–Wen F x c0 k0 图 1    Bouc–Wen 模型示意图 Fig.1    Bouc–Wen model schematic 赵义伟等： 一种描述减振器滞回特性的 Bouc–Wen 改进模型 · 1353 ·

·1354 工程科学学报，第42卷，第10期 其中，是阻尼器的相对速度，同下文dx,ms;x是 i0)=g[u,r,sgn(i0小i (4) 阻尼器相对位移，m;A,B,y,a,xo和n是控制滞回 其中，为隔振器的相对速度，ms；u为隔振器的 环形状的参数；：和是滞回阻尼力及其微分，N; 相对位移，m;r和r分别为恢复力以及它的微分，N: co是名义阻尼系数，Nsm;k是名义刚度，Nm 1为时间，s:81和g2分别是上升段和下降段滞回曲 所有参数下同 线的曲率函数，当速度大于0时，看作上升段，速 用Simulink软件建立Bouc-Wen模型，如图2所 度小于0时，即下降段 示.文献[20]对Bouc-Wen模型进行了详细的仿 然而，磁流变阻尼器的滞回特性曲线和线圈隔 真验证，并指出了模型存在的不足.Bouc-Wen 振器的有着很大区别.本文对其进行改进，可以得到 模型在识别的激励幅值下，能够准确的描述出阻 立=g(U,sgn(x)·sgn(x)·元 (5) 尼力滞回特性，但是在非识别激励振幅下不能够 此方程整理一下可得 准确地描述出滞回特性.在仿真应用中，激励幅值 三=[g1(Uのx+(0+g2(U)x-(]·元+()+ 往往不是特定的，而是随机的连续变化的.因此， (6) [g3(Ux+(t)+g4(U)x-(t]元-（①） Bouc-Wen模型在激励振幅为连续变化的系统中， 其中，U为试验得到的阻尼力，N;g(U)为试验得到 会给系统响应带来一定的误差，从而影响仿真结果 的阻尼器滞回环的斜率与阻尼力的函数关系式， 此外，试验得到的阻尼力滞回环曲线不规则时， 将阻尼器滞回环分为四段，即共有四个函数关系 Bouc-Wen模型不能够很准确地描述出来s 式，即g(U(i=1,2,3,4),进一步可得到 g1(U元 x≥0andt≥0 82(Uのx x<0and>0 = (7) g3(U文 x>0andt≤0 g4(Uの：文x≤0andx≤0 通过观察，滞回环的斜率与阻尼力的关系可 以用许多类型的初等函数来表示，如正弦函数、多 项式函数、幂函数等，还可以用分段函数，从而模 拟出不同形状的滞回环四，例如含滑移捏拢效应、 九u 刚度退化等滞回特性.论文采用二次多项式来表 示g(U),即 g(U=m;U2+nU+p:(i=1,2,3,4) (8) 针对Bouc-Wen模型对非识别激励幅值的变 Ou 化比较敏感20这一问题，对模型继续进行修正 Force 假如识别的曲线在激励幅值为10mm的情况下， 识别出的曲线与试验曲线具有准确的吻合度，然 图2 Simulink建立的Bouc-Wen模型 而当幅值为5mm时，阻尼力衰减的过快，或者当 Fig.2 Bouc-Wen model established by Simulink 幅值为15mm时，阻尼力增大的过快.修正的原理 2新型Bouc-Wen改进模型 就是：当振幅变小时，模型能够延缓力的衰减，当 振幅变大时，能够抑制阻尼力快速增大.因此，可 为了更好地描述出阻尼力滞回特性曲线弥补 以考虑在模型阻尼力公式(2)中引入一个修正力. Bouc-Wen模型的缺点，文章将此数学模型中的 如果模型中引入的修正力为一常值，经过参数识 项用另外一种形式来表示，即用滞回环斜率与阻 别后，此常数项会分担阻尼力数值的一部分，因此 尼力的函数关系来模拟滞回特性，此方法灵活性 阻尼力在不同幅值下变化就会缓慢许多，即模型 较强，有着更为广泛的应用.Ni等P利用Dunem 对振幅的敏感度下降.但如果引入常数项会破坏 微分算子对一种线圈隔振器的非线性滞回特性建 滞回环的光滑性和连续性.因此，为了保证滞回环 立了非参数模型，采用的滞回曲线方程如下： 的光滑性和连续性，论文采用指数函数作为修正 rt)=g1(u,r)·i+(t)+g2(u,r）·i(t) (3) 项，因此公式(2)可变为 此方程也可以写成 F=co+k(-x0)+az+F (9)

x˙ A, β, γ, α, x0 z˙ 其中， 是阻尼器的相对速度，同下文 dx，m·s−1 ；x 是 阻尼器相对位移，m； 和 n 是控制滞回 环形状的参数； z 和 是滞回阻尼力及其微分，N； c0 是名义阻尼系数，N·s·m−1 ；k0 是名义刚度，N·m−1 . 所有参数下同. 用Simulink 软件建立Bouc–Wen 模型，如图2 所 示. 文献 [20] 对 Bouc–Wen 模型进行了详细的仿 真验证，并指出了模型存在的不足. Bouc–Wen 模型在识别的激励幅值下，能够准确的描述出阻 尼力滞回特性，但是在非识别激励振幅下不能够 准确地描述出滞回特性. 在仿真应用中，激励幅值 往往不是特定的，而是随机的连续变化的. 因此， Bouc–Wen 模型在激励振幅为连续变化的系统中， 会给系统响应带来一定的误差，从而影响仿真结果. 此外，试验得到的阻尼力滞回环曲线不规则时， Bouc–Wen 模型不能够很准确地描述出来[16] . 2    新型 Bouc–Wen 改进模型 z˙ 为了更好地描述出阻尼力滞回特性曲线弥补 Bouc–Wen 模型的缺点，文章将此数学模型中的 项用另外一种形式来表示，即用滞回环斜率与阻 尼力的函数关系来模拟滞回特性，此方法灵活性 较强，有着更为广泛的应用. Ni 等[21] 利用 Dunem 微分算子对一种线圈隔振器的非线性滞回特性建 立了非参数模型，采用的滞回曲线方程如下： r˙(t) = g1(u,r) · u˙+(t)+g2(u,r) · u˙_(t) （3） 此方程也可以写成 r˙(t) = g[u,r,sgn( ˙u)]· u˙ （4） u˙ r˙ 其中， 为隔振器的相对速度，m·s−1 ；u 为隔振器的 相对位移，m；r 和 分别为恢复力以及它的微分，N； t 为时间，s；g1 和 g2 分别是上升段和下降段滞回曲 线的曲率函数，当速度大于 0 时，看作上升段，速 度小于 0 时，即下降段. 然而，磁流变阻尼器的滞回特性曲线和线圈隔 振器的有着很大区别. 本文对其进行改进，可以得到 z˙ = g(U,sgn(x))·sgn( ˙x)· x˙ （5） 此方程整理一下可得 z˙ = [g1(U)x+(t)+g2(U)x−(t)]· x˙+(t)+ [g3(U)x+(t)+g4(U)x−(t)]· x˙−(t) （6） 其中，U 为试验得到的阻尼力，N； g(U) 为试验得到 的阻尼器滞回环的斜率与阻尼力的函数关系式， 将阻尼器滞回环分为四段，即共有四个函数关系 式，即 gi (U)(i=1,2,3,4)，进一步可得到 z˙ =    g1(U)· x x ˙ ⩾ 0 and ˙x ⩾ 0 g2(U)· x x ˙ < 0 and ˙x > 0 g3(U)· x x ˙ > 0 and ˙x ⩽ 0 g4(U)· x x ˙ ⩽ 0 and ˙x ⩽ 0 （7） g(U) 通过观察，滞回环的斜率与阻尼力的关系可 以用许多类型的初等函数来表示，如正弦函数、多 项式函数、幂函数等，还可以用分段函数，从而模 拟出不同形状的滞回环[22] ，例如含滑移捏拢效应、 刚度退化等滞回特性. 论文采用二次多项式来表 示 ，即 gi(U) = miU 2 +niU + pi (i = 1, 2, 3, 4) （8）δF 针对 Bouc–Wen 模型对非识别激励幅值的变 化比较敏感[20] 这一问题，对模型继续进行修正. 假如识别的曲线在激励幅值为 10 mm 的情况下， 识别出的曲线与试验曲线具有准确的吻合度，然 而当幅值为 5 mm 时，阻尼力衰减的过快，或者当 幅值为 15 mm 时，阻尼力增大的过快. 修正的原理 就是：当振幅变小时，模型能够延缓力的衰减，当 振幅变大时，能够抑制阻尼力快速增大. 因此，可 以考虑在模型阻尼力公式（2）中引入一个修正力 . 如果模型中引入的修正力为一常值，经过参数识 别后，此常数项会分担阻尼力数值的一部分，因此 阻尼力在不同幅值下变化就会缓慢许多，即模型 对振幅的敏感度下降. 但如果引入常数项会破坏 滞回环的光滑性和连续性. 因此，为了保证滞回环 的光滑性和连续性，论文采用指数函数作为修正 项，因此公式（2）可变为 F = c0 +k(x− x0)+αz+δF （9） dx γ γ n z α β c0 k0 x0 β dx z n Out Out Out + + + Force Out dx A dx A x c0 k0 dx x a x0 z n −s 1 图 2    Simulink 建立的 Bouc–Wen 模型 Fig.2    Bouc–Wen model established by Simulink · 1354 · 工程科学学报，第 42 卷，第 10 期

赵义伟等：一种描述减振器滞回特性的Bouc-Wen改进模型 ·1355· 其中，dr=b·sgn()e-;b为待优化参数 3000 改进后的新模型具有17个待识别参数.接下 2000 来对修正后的模型进行参数识别 1000 0A0.5A 3MRD性能试验及模型参数识别 0 1.5A1.0A 3.1MRD特性试验分析 -1000 采用材料疲劳性能测试试验台MTS793对磁 -2000 流变阻尼器进行力学特性试验.如图3所示.对阻 2.5A2.0A -3000 0 0.5 尼器采用正弦信号激励进行加载.按不同频率和 1.0 1.5 2.0 Time/s 幅值分为8组：（0.5Hz,5mm),(1.0Hz,5mm), 图4不同电流下的阻尼力曲线图 (1.5Hz,5mm),(0.5Hz,10mm),(1.0Hz,10mm), Fig.4 Damping force curves at different currents (1.5Hz,10mm),(0.5Hz,15mm),(0.5Hz,20mm). 搜索法的优点是可以快速地得到全局的最优解 每组工况下，给阻尼器以不同的电流信号.该阻尼 尽管模式搜索法的计算效率高速度快，但是它的 器最大电流允许值为2.5A,因此电流可分为六档： 缺点在于过度地依赖初值，单独使用的时候效果 0,0.5,1.0,1.5,2.0,2.5A.最后，将采集得到的数据 不是很好.利用matlab表述为 进行处理，并画出振幅和频率分别为10mm、0.5Hz [x,fval]=patternsearch(@objfun,xx); 且不同电流情况下的力-位移曲线，如图4所示 其中，fval代表目标函数值；xx代表被识别参数初 始值 综合分析可知，传统遗传算法容易陷入局部 最优解，而模式搜索法对初值具有依赖性，因此采 用这种联合方法既考虑了遗传算法易得到局部最 优解也考虑了模式搜索法的初始值依赖性，充分 利用二者的优势，快速得到参数识别所需要的全 局最优解，识别流程图如图5所示 Read test data Set the range of initial identification parameter 图3试验装置及减振器 Fig.3 Testing device and shock absorber 3.2模型参数识别方法 Genetic algorithm(GA)parameter identification 新模型包含有17个未知参数，需要对模型中 的各个参数进行识别.论文采用的是GA-PS识别 No Satisfy the stopping condition 算法23-2 Holland教授在1962年提出了遗传算法 (Genetic algorithm).遗传算法优点是可以解决其 Yes Pattern search(PS)parameter identification 他的一些标准优化算法无法解决的问题，比如目 标函数为不连续、不可微等等.遗传算法的全局 搜索能力较强，能在较为宽广的范围内快速识别 No Satisfy the stopping condition 到需要的解.但是遗传算法的缺点是容易过早收 Yes 敛从而陷入局部最优解.利用matlab表述为 End [x,~]=ga(@fitnessfun,c,[],[].[],[],LB,UB.[],OPS): 其中，c表示被识别参数的个数；LB和UB分别表 图5基于GA和PS的参数识别过程 示参数的下界和上界；OPS表示其他选项设置. Fig.5 Parameter identification process based on GA and PS 模式搜索法(Pattern search)是一种直接搜索算 3.3识别结果及验证 法，其目标函数不一定可微也可能是非连续.模式 采用幅值10mm频率0.5Hz激励的试验数据

δF = b ·sgn( ˙x)· e 其中， −|x˙| ；b 为待优化参数. 改进后的新模型具有 17 个待识别参数. 接下 来对修正后的模型进行参数识别. 3    MRD 性能试验及模型参数识别 3.1    MRD 特性试验分析 采用材料疲劳性能测试试验台 MTS793 对磁 流变阻尼器进行力学特性试验. 如图 3 所示. 对阻 尼器采用正弦信号激励进行加载. 按不同频率和 幅值分为 8 组：（ 0.5 Hz， 5 mm），（ 1.0 Hz， 5 mm） ， （1.5 Hz，5 mm），（0.5 Hz，10 mm），（1.0 Hz，10 mm）， （1.5 Hz，10 mm），（0.5 Hz，15 mm），（0.5 Hz，20 mm）. 每组工况下，给阻尼器以不同的电流信号. 该阻尼 器最大电流允许值为 2.5 A，因此电流可分为六档： 0，0.5，1.0，1.5，2.0，2.5 A. 最后，将采集得到的数据 进行处理，并画出振幅和频率分别为 10 mm、0.5 Hz 且不同电流情况下的力–位移曲线，如图 4 所示. 3.2    模型参数识别方法 新模型包含有 17 个未知参数，需要对模型中 的各个参数进行识别. 论文采用的是 GA–PS 识别 算法[23−25] . Holland 教授在 1962 年提出了遗传算法 （Genetic algorithm）. 遗传算法优点是可以解决其 他的一些标准优化算法无法解决的问题，比如目 标函数为不连续、不可微等等. 遗传算法的全局 搜索能力较强，能在较为宽广的范围内快速识别 到需要的解. 但是遗传算法的缺点是容易过早收 敛从而陷入局部最优解. 利用 matlab 表述为 [x,∼]= ga(@fitnessfun,c,[],[],[],[],LB,UB,[],OPS); 其中，c 表示被识别参数的个数；LB 和 UB 分别表 示参数的下界和上界；OPS 表示其他选项设置. 模式搜索法（Pattern search）是一种直接搜索算 法，其目标函数不一定可微也可能是非连续. 模式 搜索法的优点是可以快速地得到全局的最优解. 尽管模式搜索法的计算效率高速度快，但是它的 缺点在于过度地依赖初值，单独使用的时候效果 不是很好. 利用 matlab 表述为 [x,fval]= patternsearch(@objfun, xx) ; 其中，fval 代表目标函数值；xx 代表被识别参数初 始值. 综合分析可知，传统遗传算法容易陷入局部 最优解，而模式搜索法对初值具有依赖性，因此采 用这种联合方法既考虑了遗传算法易得到局部最 优解也考虑了模式搜索法的初始值依赖性，充分 利用二者的优势，快速得到参数识别所需要的全 局最优解，识别流程图如图 5 所示. 3.3    识别结果及验证 采用幅值 10 mm 频率 0.5 Hz 激励的试验数据 图 3    试验装置及减振器 Fig.3    Testing device and shock absorber 0 A 0.5 A 1.5 A 2.5 A 2.0 A 1.0 A 0 0.5 1.0 1.5 2.0 −3000 −2000 −1000 0 1000 Damping force/N 2000 3000 Time/s 图 4    不同电流下的阻尼力曲线图 Fig.4    Damping force curves at different currents Yes Yes End Read test data Set the range of initial identification parameter Genetic algorithm (GA) parameter identification Pattern search (PS) parameter identification No Satisfy the stopping condition Satisfy the stopping condition No 图 5    基于 GA 和 PS 的参数识别过程 Fig.5    Parameter identification process based on GA and PS 赵义伟等： 一种描述减振器滞回特性的 Bouc–Wen 改进模型 · 1355 ·