CN113761768A - 面向整车振动抑制的磁流变阻尼器一体化优化设计方法 - Google Patents

面向整车振动抑制的磁流变阻尼器一体化优化设计方法 Download PDF

Info

Publication number
CN113761768A
CN113761768A CN202111017459.9A CN202111017459A CN113761768A CN 113761768 A CN113761768 A CN 113761768A CN 202111017459 A CN202111017459 A CN 202111017459A CN 113761768 A CN113761768 A CN 113761768A
Authority
CN
China
Prior art keywords
magneto
vehicle
whole vehicle
optimization
damper
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111017459.9A
Other languages
English (en)
Inventor
邓召学
韦鑫鑫
董轲建
蔡强
朱孙科
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Jiaotong University
Original Assignee
Chongqing Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing Jiaotong University filed Critical Chongqing Jiaotong University
Publication of CN113761768A publication Critical patent/CN113761768A/zh
Priority to NL2031692A priority Critical patent/NL2031692B1/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Vehicle Body Suspensions (AREA)

Abstract

本发明公开了一种面向整车振动抑制的磁流变阻尼器一体化优化设计方法,针对磁流变阻尼器结构与控制器之间存在的耦合等问题,以磁流变阻尼器结构参数为优化变量,以磁流变阻尼器磁感应强度为约束条件,以车辆目标性能参数的加速度均方根值以及频域幅值均取得最小值为优化目标,对磁流变阻尼器进行一体化优化设计,获得磁流变阻尼器在整车动力学模型下的全局最优结构参数。本发明能够兼顾车辆整车性能,提升车辆的平顺性和操稳性,可满足不同工况下的隔振需求,实用性强。

Description

面向整车振动抑制的磁流变阻尼器一体化优化设计方法
技术领域
本发明涉及车辆零部件设计领域,具体涉及一种面向整车振动抑制的磁流变阻尼器一体化优化设计方法。
背景技术
为满足车辆行驶时复杂工况对隔振性能的不同要求,性能卓越的磁流变阻尼器应具有较大的可控阻尼力和较小的粘性阻尼力。然而,在获取较大的可控阻尼力时,粘性阻尼力也随之增加。为了折衷性能指标之间的矛盾,部分学者针对磁流变阻尼器单体性能展开了相关研究,而磁流变阻尼器作为底盘的重要子系统,在割裂阻尼器与车轮和底盘的耦合振动关系时,会忽略结构与控制之间存在的耦合问题,孤立地研究磁流变阻尼器本体的优化设计,不能兼顾整车性能,很难实现磁流变阻尼器与整车NVH品质之间的良好匹配关系。
发明内容
有鉴于此,本发明的目的是克服现有技术中的缺陷,提供面向整车振动抑制的磁流变阻尼器一体化优化设计方法,能够兼顾车辆整车性能,提升车辆的平顺性和操稳性,可满足不同工况下的隔振需求,实用性强。
本发明的面向整车振动抑制的磁流变阻尼器一体化优化设计方法,包括:
考虑结构设计和控制器之间的耦合,以磁路磁感应强度为约束条件,以结构参数为优化变量,以目标工况下车辆目标性能参数的加速度均方根值以及频域幅值最小为优化目标,构建整车动力学一体化优化模型;并使用整车动力学一体化优化模型对磁流变阻尼器进行优化设计,获得磁流变阻尼器的优化结构参数。
进一步,所述整车动力学一体化优化模型包括7自由度的整车动力学模型;所述7自由度包括车身垂向运动、车身侧倾运动、车身俯仰运动以及4个车轮的垂向运动。
进一步,设计模糊PID控制器,并使用模糊PID控制器对磁流变阻尼器的电流进行实时控制。
进一步,对磁路进行参数化建模与仿真分析,并将得到的磁路磁感应强度作为约束条件。
进一步,基于结构/控制一体化平台,筛选出结构参数,并将所述结构参数作为优化变量。
进一步,选取表征平顺性和操稳性的评价指标,并将所述评价指标作为优化目标。
本发明的有益效果是:本发明公开的一种面向整车振动抑制的磁流变阻尼器一体化优化设计方法,通过考虑磁流变阻尼器结构和控制之间的耦合关系,构建整车动力学一体化优化模型,对磁流变阻尼器的结构参数进行优化设计,使得车辆的目标性能参数的加速度均方根值以及频域值均取得最小值,进而实现对车辆平顺性和操稳性的提升,兼顾了车辆整车性能,满足了不同工况下的隔振需求。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1为本发明的原理设计示意图;
图2为本发明的匀速行驶工况下一体化优化设计前后的时域图;
图3为本发明的匀速行驶工况下一体化优化设计前后的频域图;
图4为本发明的过减速带工况下一体化优化设计前后的时域图;
图5为本发明的过减速带工况下一体化优化设计前后的频域图。
具体实施方式
以下结合说明书附图对本发明做出进一步的说明,如图1所示:
本发明的面向整车振动抑制的磁流变阻尼器一体化优化设计方法,包括:
考虑结构设计和控制器之间的耦合,以磁路磁感应强度作为约束条件,以磁流变阻尼器的结构参数作为优化变量,以目标工况下车辆目标性能参数的加速度均方根值以及频域幅值最小为优化目标,构建整车动力学一体化优化模型;并使用整车动力学一体化优化模型对磁流变阻尼器进行优化设计,获得磁流变阻尼器的优化结构参数。其中,所述结构参数包括阻尼间隙、倾斜角度、内径尺寸以及磁芯长度;所述目标工况包括匀速行驶工况和过减速带工况;所述目标性能参数包括悬架动挠度、车身垂向加速度以及轮胎动载荷。
本实施例中,所述整车动力学一体化优化模型包括7自由度的整车动力学模型;所述7自由度包括车身垂向运动、车身侧倾运动、车身俯仰运动以及4个车轮的垂向运动;其中,可以采用动力学仿真软件搭建7自由度的整车动力学模型。
本实施例中,通过考虑结构和控制之间的耦合,设计一款实用性强的模糊PID控制器,使用模糊PID控制器控制整车动力学一体化优化模型,将所述优化模型输出的优化目标的值输入到所述控制器而形成正反馈,将所述控制器输出的电流输入到所述优化模型而形成负反馈,通过上述设计,从而形成一个系统闭环。并通过采用现有的智能优化算法来计算车辆的目标性能参数的加速度均方根值以及频域幅值;其中,现有的智能优化算法可以为非支配排序遗传算法。
本实施例中,根据实际工况设定约束条件的经验值,并将所述经验值作为磁路磁感应强度最小值;选取磁流变液及磁芯的磁饱和值作为磁路磁感应强度最大值;则将磁路磁感应强度限定在所述磁路磁感应强度最小值与所述磁路磁感应强度最大值之间。
本实施例中,选取表征平顺性和操稳性的常用评价指标,并将所述评价指标作为优化目标,通过时频分析方法获得优化目标的最小值,并将取得最小值时设置的优化变量作为磁流变阻尼器的最优结构参数。
本实施例中,基于结构/控制一体化优化平台,利用灵敏度分析方法筛选出对车辆目标性能影响较大的结构参数,研究磁流变阻尼器结构参数对车辆NVH品质的影响规律,为一体化优化变量的制定提供决策依据;其中,为实现磁流变阻尼器与整车NVH品质之间的优良匹配关系,提取车辆振动观测点目标工况下振动响应的时频特征,通过分析时频特征量与整车固有特性、外部激励等因素之间的关联性,揭示其映射关系。
本实施例中,根据如下方法确定磁路磁感应强度:
对磁路进行参数化建模,得到建模后的磁路,对建模后的磁路的电磁场进行仿真分析,得到磁路磁感应强度,并将所述磁路磁感应强度作为约束条件。其中,在有限元软件内对磁路进行APDL参数化建模。
本实施例中,根据如下方法确定结构参数:
获取初始结构参数;其中,所述初始结构参数包括线圈槽长度、阻尼间隙、倾斜角度、内径尺寸以及磁芯长度;
通过灵敏度分析,从初始结构参数中筛选出对车辆时域和频域性能影响较大的结构参数,并将所述结构参数作为优化变量。
所述有限元软件采用现有软件,如ANSYS、ABAQUS、Hypermesh等;使用动力学仿真软件进行仿真分析,所述动力学仿真软件采用现有软件,如:MATLAB/Simulink、Adams、Carsim等;
本实施例中,根据如下公式确定车身垂向运动:
Figure BDA0003240437610000041
其中,mb为车身质量,
Figure BDA0003240437610000042
为车身质心处的垂向加速度,ks1为前轴左侧悬架弹簧的刚度系数,ks2为前轴右侧悬架弹簧的刚度系数,ks3为后轴左侧悬架弹簧的刚度系数,ks4为后轴右侧悬架弹簧的刚度系数;zb1为车身在前轴左侧悬架处的垂向位移,zb2为车身在前轴右侧悬架处的垂向位移,zb3为车身在后轴左侧悬架处的垂向位移,zb4为车身在后轴右侧悬架处的垂向位移;
Figure BDA0003240437610000043
为车身在前轴左侧悬架处的垂向速度,
Figure BDA0003240437610000044
为车身在前轴右侧悬架处的垂向速度,
Figure BDA0003240437610000045
为车身在后轴左侧悬架处的垂向速度,
Figure BDA0003240437610000051
为车身在后轴右侧悬架处的垂向速度;zw1为前轴左侧车轮的垂向位移,zw2为前轴右侧车轮的垂向位移,zw3为后轴左侧车轮的垂向位移,zw4为后轴右侧车轮的垂向位移;
Figure BDA0003240437610000052
为前轴左侧车轮的垂向速度,
Figure BDA0003240437610000053
为前轴右侧车轮的垂向速度,
Figure BDA0003240437610000054
为后轴左侧车轮的垂向速度,
Figure BDA0003240437610000055
为后轴右侧车轮的垂向速度;c1为前轴左侧悬架弹簧阻尼系数,c2为前轴右侧悬架弹簧阻尼系数,c3为后轴左侧悬架弹簧阻尼系数,c4为后轴右侧悬架弹簧阻尼系数,f1为前轴左侧磁流变阻尼器提供的阻尼器,f2为前轴右侧磁流变阻尼器提供的阻尼器,f3为后轴左侧磁流变阻尼器提供的阻尼器,f4为后轴右侧磁流变阻尼器提供的阻尼器;
根据如下公式确定车身侧倾运动:
Figure BDA0003240437610000056
其中,Ir为车身绕其质心的侧倾转动惯量,
Figure BDA0003240437610000057
为车身侧倾角加速度,a、b分别为前后车轮的轮距;
根据如下公式确定车身俯仰运动:
Figure BDA0003240437610000058
其中,Ip为车身绕其质心的俯仰转动惯量,
Figure BDA0003240437610000059
为车身俯仰角加速度,Lf、Lr分别为前后轮的中心轴线到车身质心的距离;
根据如下公式确定前轴左侧车轮的垂向运动:
Figure BDA00032404376100000510
其中,m1为前左的簧下质量,
Figure BDA00032404376100000511
为前轴左侧车轮的垂向加速度,kt1为前左轮胎刚度系数;
根据如下公式确定前轴右侧车轮的垂向运动:
Figure BDA0003240437610000061
其中,m2为前右的簧下质量,
Figure BDA0003240437610000062
为前轴右侧车轮的垂向加速度,kt2为前右轮胎刚度系数;
根据如下公式确定后轴左侧车轮的垂向运动:
Figure BDA0003240437610000063
其中,m3为后左的簧下质量,
Figure BDA0003240437610000064
为后轴左侧车轮的垂向加速度,kt3为后左轮胎刚度系数;
根据如下公式确定后轴右侧车轮的垂向运动:
Figure BDA0003240437610000065
其中,m4为后右的簧下质量,
Figure BDA0003240437610000066
为后轴右侧车轮的垂向加速度,kt4为后右轮胎刚度系数;
由于车身俯仰角
Figure BDA0003240437610000067
和车身侧倾角θ均是在小范围内变化,故可以得到四角处簧上质量垂向位移坐标与车身质心坐标存在如下的位移约束关系式:
Figure BDA0003240437610000068
Figure BDA0003240437610000069
Figure BDA00032404376100000610
Figure BDA00032404376100000611
其中,zb为车身质心处的垂向位移,
Figure BDA00032404376100000612
为车身绕其质心的俯仰角位移,θ为车身绕其质心的侧倾角位移。
为了能更好地理解本发明,对使用现有方法(优化设计前)与使用本发明(优化设计后)进行对比分析如下:
图2和图3分别为车辆在B级路面60km/h匀速行驶,优化前后时域图形和频域图形对比,可以看出,使用本发明进行优化设计后的磁流变阻尼器相比于初始的磁流变阻尼器,悬架动挠度、轮胎动载荷、车身垂向加速度在时域和频域范围内均得到不同程度的改善。
图4和图5分别为车辆在B级路面10km/h过减速带行驶,优化前后时域图形和频域图形对比,也可以看出,使用本发明进行优化设计后的磁流变阻尼器相比于初始的磁流变阻尼器,悬架动挠度、轮胎动载荷、车身垂向加速度在时域和频域范围内均得到不同程度的改善。
综合上述分析结果可知,在B级路面60Km/h匀速行驶工况和车辆10Km/h过减速带工况下,使用本发明进行优化设计的磁流变阻尼器的隔振性能与初始设计的磁流变阻尼器相比较,悬架动挠度、轮胎动载荷和车身垂向加速度都有一定的提高;车辆的平顺性和操稳性都得到了相应的改善,NVH性能得到明显提升。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (6)

1.一种面向整车振动抑制的磁流变阻尼器一体化优化设计方法,其特征在于:考虑结构设计和控制器之间的耦合,以磁路磁感应强度为约束条件,以结构参数为优化变量,以目标工况下车辆目标性能参数的加速度均方根值以及频域幅值最小为优化目标,构建整车动力学一体化优化模型;并使用整车动力学一体化优化模型对磁流变阻尼器进行优化设计,获得磁流变阻尼器的优化结构参数。
2.根据权利要求1所述的面向整车振动抑制的磁流变阻尼器一体化优化设计方法,其特征在于:所述整车动力学一体化优化模型包括7自由度的整车动力学模型;所述7自由度包括车身垂向运动、车身侧倾运动、车身俯仰运动以及4个车轮的垂向运动。
3.根据权利要求1所述的面向整车振动抑制的磁流变阻尼器一体化优化设计方法,其特征在于:设计模糊PID控制器,并使用模糊PID控制器对磁流变阻尼器的电流进行实时控制。
4.根据权利要求1所述的面向整车振动抑制的磁流变阻尼器一体化优化设计方法,其特征在于:对磁路进行参数化建模与仿真分析,并将得到的磁路磁感应强度作为约束条件。
5.根据权利要求1所述的面向整车振动抑制的磁流变阻尼器一体化优化设计方法,其特征在于:基于结构/控制一体化平台,筛选出结构参数,并将所述结构参数作为优化变量。
6.根据权利要求1所述的面向整车振动抑制的磁流变阻尼器一体化优化设计方法,其特征在于:选取表征平顺性和操稳性的评价指标,并将所述评价指标作为优化目标。
CN202111017459.9A 2021-07-30 2021-08-31 面向整车振动抑制的磁流变阻尼器一体化优化设计方法 Pending CN113761768A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
NL2031692A NL2031692B1 (en) 2021-08-31 2022-04-25 Integrated optimization design method of magnetorheological damper for whole vehicle vibration suppression

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110875907 2021-07-30
CN2021108759072 2021-07-30

Publications (1)

Publication Number Publication Date
CN113761768A true CN113761768A (zh) 2021-12-07

Family

ID=78792375

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111017459.9A Pending CN113761768A (zh) 2021-07-30 2021-08-31 面向整车振动抑制的磁流变阻尼器一体化优化设计方法

Country Status (1)

Country Link
CN (1) CN113761768A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL2031692A (en) * 2021-08-31 2023-03-09 Univ Chongqing Jiaotong Integrated optimization design method of magnetorheological damper for whole vehicle vibration suppression
CN117521229A (zh) * 2023-12-27 2024-02-06 石家庄铁道大学 一种路面位移响应检测方法、系统及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112257182A (zh) * 2020-10-26 2021-01-22 重庆交通大学 面向整车振动抑制的磁流变悬置时频特征多目标优化方法
CN112560299A (zh) * 2020-10-26 2021-03-26 重庆交通大学 基于整车动力学模型的磁流变阻尼器结构参数优化方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112257182A (zh) * 2020-10-26 2021-01-22 重庆交通大学 面向整车振动抑制的磁流变悬置时频特征多目标优化方法
CN112560299A (zh) * 2020-10-26 2021-03-26 重庆交通大学 基于整车动力学模型的磁流变阻尼器结构参数优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
盛鹏程: "基于磁流变阻尼器的发动机振动模糊PID控制", 《汽车工程》, vol. 31, no. 2, 25 February 2009 (2009-02-25), pages 123 - 126 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL2031692A (en) * 2021-08-31 2023-03-09 Univ Chongqing Jiaotong Integrated optimization design method of magnetorheological damper for whole vehicle vibration suppression
CN117521229A (zh) * 2023-12-27 2024-02-06 石家庄铁道大学 一种路面位移响应检测方法、系统及存储介质
CN117521229B (zh) * 2023-12-27 2024-03-19 石家庄铁道大学 一种路面位移响应检测方法、系统及存储介质

Similar Documents

Publication Publication Date Title
Theunissen et al. Preview-based techniques for vehicle suspension control: A state-of-the-art review
Ming et al. Semi-active suspension control based on deep reinforcement learning
CN113761768A (zh) 面向整车振动抑制的磁流变阻尼器一体化优化设计方法
CN110001337B (zh) 一种基于add正实网络优化的车辆isd悬架二阶理想模型的建立方法
CN103754081B (zh) 车辆非线性悬架系统的最优模糊复合控制方法
WO2023174327A1 (zh) 一种协同控制模块、自适应巡航系统及其控制方法、交通工具
CN113591360B (zh) 基于整车动力学模型的磁流变阻尼器结构参数优化方法
Corno et al. Active control of aerodynamic surfaces for ride control in sport vehicles
CN112356633A (zh) 一种考虑时滞干扰的车辆主动悬架系统的自适应控制方法
Corno et al. Performance assessment of active aerodynamic surfaces for comfort and handling optimization in sport cars
CN108859648A (zh) 一种基于神经网络的悬架减震器阻尼控制切换加权系数确定方法
CN113569340A (zh) 面向整车振动抑制的磁流变悬置时频特征多目标优化方法
CN110001336B (zh) 一种基于add正实网络优化的车辆isd悬架一阶理想模型
JP2007233985A (ja) システムの最適制御方法
CN112906133B (zh) 一种地棚控制动惯性悬架的垂向振动负效应抑制方法
CN110901326A (zh) 一种具有状态约束与死区输入的主动悬架系统的控制方法
CN111439086B (zh) 一种基于粒子群优化的车辆isd悬架pdd控制理想模型
CN115782496B (zh) 一种基于map控制的半主动悬架系统智能进化方法
Marzbani et al. Smart flat ride tuning
Gustafsson et al. Neural network controller for semi-active suspension systems with road preview
Feng et al. Bandwidth-limited active suspension controller for an off-road vehicle based on co-simulation technology
CN111231595B (zh) 一种考虑汽车前后轴动态耦合的半主动悬架控制方法
NL2031692B1 (en) Integrated optimization design method of magnetorheological damper for whole vehicle vibration suppression
CN111444660A (zh) 新能源车用机电悬架模型构建方法、优化仿真方法及系统
CN117261521B (zh) 一种高度和阻尼自适应的半主动悬架智能控制与进化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination