CN103522862B - 一种确定半主动悬架等效阻尼最大值的方法 - Google Patents

一种确定半主动悬架等效阻尼最大值的方法 Download PDF

Info

Publication number
CN103522862B
CN103522862B CN201310475719.6A CN201310475719A CN103522862B CN 103522862 B CN103522862 B CN 103522862B CN 201310475719 A CN201310475719 A CN 201310475719A CN 103522862 B CN103522862 B CN 103522862B
Authority
CN
China
Prior art keywords
semi
suspension
damping
maxim
active suspension
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.)
Expired - Fee Related
Application number
CN201310475719.6A
Other languages
English (en)
Other versions
CN103522862A (zh
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.)
Jiangsu University
Original Assignee
Jiangsu 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 Jiangsu University filed Critical Jiangsu University
Priority to CN201310475719.6A priority Critical patent/CN103522862B/zh
Publication of CN103522862A publication Critical patent/CN103522862A/zh
Application granted granted Critical
Publication of CN103522862B publication Critical patent/CN103522862B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开一种车辆悬架领域中确定半主动悬架等效阻尼最大值的方法,半主动悬架包括簧载质量、悬架弹簧、粘性阻尼可调减震器和控制器,采用数值仿真方法确定二次型性能指标值和实时所需的理想半主动控制力;以0.1为初值、0.1为公比的等比数列设置一系列悬架相对运动速度的下限值,通过仿真得到不同下限值对应的二次型性能指标值和阻尼可调减震器实时阻尼的标准差,根据二次型性能指标值设计要求,取满足性能要求的最大相对运动速度下限值及其对应的标准差;根据概率统计理论和标准差,以满足悬架使用概率要求确定粘性阻尼的最大值;本发明为半主动悬架的制造成本、安装方便性及使用性能之间的综合平衡提供设计依据,简单易行。

Description

一种确定半主动悬架等效阻尼最大值的方法
技术领域
本发明属于车辆悬架技术领域,尤其涉及一种确定半主动悬架等效阻尼最大值的方法。
背景技术
半主动悬架与主动悬架相比无需外界动力源,能取得比传统被动悬架好得多的减振效果。如图1所示的1/4车半主动悬架的结构,在铅垂方向上(与悬架运动相关的所有物理量的方向都是在铅垂方向上),车轮质量6与轮胎等效弹簧7组成车轮,车轮位于簧载质量1的下方,车轮与簧载质量1之间并联有悬架弹簧2与阻尼可调减震器4,车轮直接与地面相互作用而使悬架产生振动;在簧载质量1上固定设有簧载质量加速度传感器3,在车轮质量6上固定设有车轮质量加速度传感器5,簧载质量加速度传感器3与车轮质量加速度传感器5各自通过信号线连接于控制器8,阻尼可调减震器4也通过信号线连接于控制器8,控制器8可采用PID或LQG控制、模糊控制等控制策略,这几种控制策略可单独使用或组合使用。
从理论上讲,半主动悬架等效阻尼的最大值越大,阻尼可调减震器的工作能力越强,半主动悬架的使用性能越好。此外,等效阻尼的最大值也与阻尼可调减震器制造成本及安装方便性有关,等效阻尼的最大值越大,阻尼可调减震器制造成本越高,体积也越大,安装也就越不方便。因此对于半主动悬架等效阻尼最大值的确定,应在制造成本、安装方便性及使用性能三者之间做综合平衡的考虑。
现有半主动悬架减震器等效阻尼最大值的确定一般根据经验选取,并未在理论和技术手段上提供一种确定半主动悬架阻尼可调减震器等效阻尼最大值的方法。
发明内容
本发明提供一种确定半主动悬架等效阻尼最大值的方法,该方法简单精确,能提高半主动悬架的使用性能和安装方便性,同时降低制造成本。
为实现上述目的,本发明采用的技术方案是:半主动悬架包括簧载质量、悬架弹簧、阻尼可调减震器和控制器,车轮与簧载质量之间并联有悬架弹簧与粘性阻尼可调减震器,包括以下步骤:(1)采用数值仿真方法确定悬架二次型性能指标值J和实时的理想半主动控制力F idsa ;(2)以0.1为初值、0.1为公比的等比数列设置一系列悬架相对运动速度Δv的下限值v flr ,通过仿真得到不同相对运动速度下限值v flr 对应的阻尼可调减震器实时阻尼c vvd 的阻尼标准差;(3)根据半主动悬架设计对二次型性能指标值J的设计要求,取满足该要求的最大相对运动速度下限值v flr 及其对应的阻尼标准差;(4)根据概率统计理论和步骤(3)所述阻尼标准差,以满足半主动悬架的使用概率要求确定粘性阻尼的最大值。
本发明简单易行、便于操作,为半主动悬架的制造成本、安装方便性及使用性能三者之间的综合平衡提供设计依据,同样适用于其它类型半主动悬架等效阻尼最大值的确定,可为各种类型的半主动悬架的设计及相关的理论与试验研究提供基础,具有重要的实际工程意义。
附图说明
图1是背景技术中1/4车半主动悬架的结构与控制原理示意图。
图中:1.簧载质量;2.悬架弹簧;3.簧载质量加速度传感器;4.阻尼可调减震器;5.车轮质量加速度传感器;6.车轮质量;7.轮胎等效弹簧;8.控制器。
具体实施方式
本发明以传统的粘性阻尼可调减震器半主动悬架以及其LQG控制策略为例,其它类型的半主动悬架及控制策略与本发明雷同,不再赘述,具体如下:
1、采用数值仿真方法确定实时所需的理想半主动控制力,参见图1,根据牛顿第二运动定律进行动力学分析有:
(1)
式中:m 1为车轮质量;m 2为簧载质量;q为悬架系统的路面位移输入,,其中:n 0为参考空间频率;w为高斯分布白噪声,为路面不平度系数;u为汽车行驶速度;f 0为下截止频率;x 1为车轮垂直位移;为车轮速度;为车轮加速度;x 2为车身垂直位移;为车身速度;车身加速度;k 1为轮胎刚度;k 2为悬架刚度;c s 为半主动悬架的结构阻尼;F为半主动悬架的实际控制力。
选取悬架系统状态空间变量X=(q,x 1,x 2,,)T,结合式(1)得到悬架系统运动的状态方程为:
(2)
式中:
J为悬架二次型性能指标值,J值越小则悬架综合性能越好。
(3)
式中:为车身加速度;x 1-q为车轮动变形;x 2-x 1为悬架动挠度;δ 1δ 2分别为(x 1-q)2与(x 2-x 1)2的加权系数;T为车辆行驶时间。
半主动悬架的理想半主动控制力F idsa 表示如下:
(4)
其中,理想主动控制力F ida 通过LQG控制得到:
F ida =- YX (5)
( Y , S , E )=LQR( A , B , Q , R , N )(6)
LQR()是matlab软件的自带函数,利用LQR()函数求出Y,因X=(q,x 1,x 2,,)T,可由式(5)求出理想主动控制力F ida
2、按以0.1为初值、以0.1为公比的等比数列设置一系列具有显著差异的悬架相对运动速度Δv的下限值v flr ,通过仿真得到不同相对运动速度下限值v flr 对应的悬架二次型性能指标J和阻尼标准差,悬架二次型性能指标J按式(3)计算。
理论上,半主动悬架的阻尼可调减震器实时阻尼c vvd 等于理想半主动控制力F idsa 除以相对运动速度Δv,这意味着如果Δv→0,c vvd →∞。为了避免c vvd 出现无穷大,设置Δv的下限值v flr 。粘性阻尼可调减震器实时阻尼c vvd 按式(7)求取,阻尼标准差按式(8)计算:
(7)
(8)
当车轮质量m 1=350kg、簧载质量m 2=5000kg、k 1=3000000N/m,k 2=505000N/m、c s=1000N·s/m、n =0.1m-1T=20s、u=20m/s、f 0=0.011uδ 1=5.2894×104δ 2=4.4051×103,阻尼不变的原型被动悬架的阻尼值为30150N·s/m时,相对运动速度下限值v flr 分别等于0.1、0.01、0.001和0.0001m·s-1所对应的半主动悬架的二次型性能指标值J与阻尼标准差以及理想主动悬架的二次型性能指标值J统计结果如表1所列;
表1
3、以阻尼标准差与被动悬架减震器阻尼值相近为基准,根据半主动悬架设计的悬架性能要求(例如:悬架二次型性能指标值J相对理想状态的增加不超1.5%或2%等特定设计要求),取满足悬架性能要求的最大相对运动速度下限值v flr 以及该下限值v flr 对应的阻尼标准差
表1显示当v flr 等于0.001m·s-1,半主动悬架粘性阻尼可调减震器的J值基本不增加,说明当v flr 等于0.001m·s-1,半主动悬架的J值已经接近理想值。考虑到当v flr 等于0.1m·s-1时半主动悬架的J值仅比v flr 等于0.001m·s-1时半主动悬架的J值大1.25%,满足悬架的二次型悬架性能指标值J相对理想状态的增加不超过1.5%要求,且v flr 等于0.1m·s-1时的阻尼标准差与被动悬架减震器的阻尼值非常接近,且仅等于v flr 为0.001m·s-1时阻尼标准差的0.25%。
4、根据概率统计理论和上述确定的阻尼标准差,即上述确定的最大相对运动速度下限值v flr 对应的阻尼标准差,以满足半主动悬架的特定工作能力为目标确定粘性阻尼的最大值,例如以满足99%或95%等的使用概率要求确定阻尼的最大值。下表2是标准差倍数与概率统计的对应关系表。
表2
根据表2,将阻尼可调减震器实时阻尼c vvd 代替实时阻尼χ,代替σ χ ,这样,实时阻尼χ最大值超过限值χ 0λσ χ 的概率P χ 可由表2的λP χ 的关系查得,也就是c vvd 最大值超过限值χ 0λ 的概率P χ 可由表2的λP χ 的关系查得。表2中λ为某一实时阻尼χ最大值的限值χ 0相对于标准差σ χ 的倍数。
表2显示:当λ为2.00时可满足95.4%的使用概率(工作需求);当λ为2.58时可满足99%的使用概率;当λ为3.00可满足99.7%的使用概率;考虑到实际工程误差,本发明使用λ为2.58时,即当c vvd 的最大值等于c vvd 标准差的2.58倍,满足99%的工作要求,确定c vvd 的最大值等于92271N·s·m-1,此时J值等于3.6073较理想半主动悬架的J值3.5594约大1.35%,而较理想主动悬架的J值3.3538约大7.56%。

Claims (1)

1.一种确定半主动悬架等效阻尼最大值的方法,半主动悬架包括簧载质量、悬架弹簧、阻尼可调减震器和控制器,车轮与簧载质量之间并联有悬架弹簧与阻尼可调减震器,其特征是包括以下步骤:
1)采用数值仿真方法确定悬架二次型性能指标值J和实时的理想半主动控制力F idsa ;所述二次型性能指标值
所述实时的理想半主动控制力
T为车辆行驶时间,q为悬架路面位移,x 1为车轮垂直位移,为车轮速度;为车轮加速度;x 2为车身垂直位移;为车身速度;车身加速度;
x 1-q为车轮动变形;x 2-x 1为悬架动挠度;δ 1δ 2分别为(x 1-q)2与(x 2-x 1)2的加权系数,F ida 为理想主动控制力;
(2)以0.1为初值、0.1为公比的等比数列设置一系列悬架相对运动速度Δv的下限值v flr ,通过仿真得到不同相对运动速度下限值v flr 对应的阻尼可调减震器实时阻尼c vvd 的阻尼标准差;所述阻尼标准差
T为车辆行驶时间,
(3)根据半主动悬架设计对二次型性能指标值J的设计要求,取满足该要求的最大相对运动速度下限值v flr 及其对应的阻尼标准差
(4)根据概率统计理论和步骤(3)所述阻尼标准差,以满足半主动悬架的使用概率要求确定半主动悬架等效阻尼的最大值;阻尼可调减震器实时阻尼c vvd 最大值超过限值χ 0λ 的概率P χ 由限值χ 0相对于标准差的倍数λ与概率统计的对应关系表查得:
CN201310475719.6A 2013-10-14 2013-10-14 一种确定半主动悬架等效阻尼最大值的方法 Expired - Fee Related CN103522862B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310475719.6A CN103522862B (zh) 2013-10-14 2013-10-14 一种确定半主动悬架等效阻尼最大值的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310475719.6A CN103522862B (zh) 2013-10-14 2013-10-14 一种确定半主动悬架等效阻尼最大值的方法

Publications (2)

Publication Number Publication Date
CN103522862A CN103522862A (zh) 2014-01-22
CN103522862B true CN103522862B (zh) 2015-12-09

Family

ID=49925363

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310475719.6A Expired - Fee Related CN103522862B (zh) 2013-10-14 2013-10-14 一种确定半主动悬架等效阻尼最大值的方法

Country Status (1)

Country Link
CN (1) CN103522862B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105751847B (zh) * 2016-03-24 2018-02-27 江苏大学 一种车辆多模式减振器的控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102582389A (zh) * 2012-03-02 2012-07-18 江苏大学 一种基于虚拟阻尼的车辆悬架控制方法
CN102729760A (zh) * 2012-07-17 2012-10-17 山东理工大学 汽车半主动悬架系统实时最佳阻尼控制算法
CN103204043A (zh) * 2013-04-01 2013-07-17 中国人民解放军装甲兵工程学院 汽车半主动悬架系统的频域控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7321816B2 (en) * 2004-04-27 2008-01-22 Tenneco Automotive Operating Company, Inc. Model free semi-active vehicle suspension system
JP5700190B2 (ja) * 2009-08-04 2015-04-15 日立オートモティブシステムズ株式会社 サスペンション制御装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102582389A (zh) * 2012-03-02 2012-07-18 江苏大学 一种基于虚拟阻尼的车辆悬架控制方法
CN102729760A (zh) * 2012-07-17 2012-10-17 山东理工大学 汽车半主动悬架系统实时最佳阻尼控制算法
CN103204043A (zh) * 2013-04-01 2013-07-17 中国人民解放军装甲兵工程学院 汽车半主动悬架系统的频域控制方法

Also Published As

Publication number Publication date
CN103522862A (zh) 2014-01-22

Similar Documents

Publication Publication Date Title
Shen et al. Improved design of dynamic vibration absorber by using the inerter and its application in vehicle suspension
Sun et al. Saturated adaptive robust control for active suspension systems
CN102582389B (zh) 一种基于虚拟阻尼的车辆悬架控制方法
CN106647256B (zh) 一种基于h∞ pid的主动悬架控制系统及控制方法
CN103079849A (zh) 车辆的控制装置
CN109334376B (zh) 磁流变半主动悬架的减振器力学特性时变补偿控制系统及其构造方法
CN106926660B (zh) 一种基于轮边驱动电动车的电磁悬架系统及其控制方法
CN105539052A (zh) 一种以车辆平稳状态为参考的可控悬架滑模追踪控制器
CN109334377B (zh) 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制系统及其构造方法
CN110321665A (zh) 车辆悬架系统的控制方法、装置、设备、介质及车辆
Wei et al. A new method of static output-feedback H∞ controller design for 5 DOF vehicle active suspension system
CN111055650A (zh) 磁流变半主动悬架粒子群-时滞依赖h∞鲁棒控制方法
Sulaiman et al. Groundhook control of semi-active suspension for heavy vehicle
CN102975587A (zh) 基于双可控阻尼器的车辆半主动悬架及其控制方法
CN103522862B (zh) 一种确定半主动悬架等效阻尼最大值的方法
Shi et al. Feedback linearization and sliding mode control for active hydropneumatic suspension of a special-purpose vehicle
Ismail et al. A control performance of linear model and the MacPherson model for active suspension system using composite nonlinear feedback
Tran et al. Advanced passive suspension with Inerter devices and optimization design for vehicle oscillation
CN103034253A (zh) 动态物理过程振动抑制方法和系统
Yue et al. Investigation of semi-active hydro-pneumatic suspension for a heavy vehicle based on electro-hydraulic proportional valve
Junyao et al. Vibration control for vehicle active suspension based on ANFIS method
Chen et al. Stiffness-damping matching method of an ECAS system based on LQG control
JP2013241076A (ja) サスペンション制御装置
Najm et al. Mathematical Modelling and PID Controller Implementation to Control Linear and Nonlinear Quarter Car Active Suspension
Jibril Comparison of active and semi-active suspension systems using robust controller

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151209

Termination date: 20161014

CF01 Termination of patent right due to non-payment of annual fee