CN100358739C - 一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法 - Google Patents
一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法 Download PDFInfo
- Publication number
- CN100358739C CN100358739C CNB2004100401179A CN200410040117A CN100358739C CN 100358739 C CN100358739 C CN 100358739C CN B2004100401179 A CNB2004100401179 A CN B2004100401179A CN 200410040117 A CN200410040117 A CN 200410040117A CN 100358739 C CN100358739 C CN 100358739C
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- sigma
- vector
- max
- 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
Links
Images
Landscapes
- Vehicle Body Suspensions (AREA)
Abstract
本发明涉及一种基于小波分解的汽车半主动悬架阻尼控制方法,它是利用加速度传感器在线提取汽车行驶中车身的振动信号,运用小波变换的方法将振动信号分解为由高到低的若干频段,通过能量统计的方法确定当前振动信号主要分布在哪几个频段,根据阻尼系数与车轴、车身振动响应曲线确定通过可控阻尼减振器调节悬架阻尼系数,从而使汽车获得良好的平顺性和操纵稳定性。本方法运用小波变换对车身振动信号进行多分辨率分析,高频信号部分分辨率低,低频信号分辨率高,符合汽车振动信号的特点。由于算法只测量簧载质量的加速度,不需要预测路面,兼顾悬架低频振动和高频振动,能够方便地应用到汽车半主动悬架上,实现实时控制。
Description
技术领域
本发明属于汽车减振控制技术领域,具体涉及一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法。
背景技术
自1973年美国的D.A.Crosby和D.C.Karnopp提出半主动悬架的概念以来,各国学者陆续开展了基于可控阻尼减振器的半主动悬架控制方法研究,其中在实际车辆上得到运用的是天棚阻尼控制,它设想将减振器设置在簧载质量和惯性坐标-天棚之间,直接控制簧载质量的绝对运动速度。天棚阻尼控制算法简单,研究表明只要合理选择参数,能够很好的抑制簧载质量的垂直振动,提高乘坐舒适性。但是天棚阻尼需要同时测量簧载质量的绝对速度和簧载质量与非簧载质量的相对速度,而簧载质量绝对速度的获取相当困难,天棚阻尼还要求执行机构的时滞较小,否则容易导致“开关”控制的误动作,并且它以牺牲轮胎着地性为代价,在某些工况下是不可接受的。
最优控制、模糊控制、自适应控制、预见控制等控制理论也被应用于汽车半主动减振控制,但是这些方法常常是建立在各种假设的基础之上的,与实际工况有较大的差异;并且大多计算复杂,收敛慢,实时性差,因此主要停留在实验室仿真,很难应用于实车控制。
发明内容
发明目的:本发明提出一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法,运用小波变换的方法来分解振动信号,其目的是同时降低簧载质量和非簧载质量的垂直加速度,提高汽车的乘座舒适性和操纵稳定性。
本发明的技术解决方案如下:
本发明方法的总体方案是:利用加速度传感器在线提取汽车行驶中车身的振动信号,运用小波变换的方法将振动信号分解为由高到低的若干频段,通过能量统计的方法确定当前振动信号主要分布在哪几个频段,根据阻尼系数与车轴、车身振动响应曲线确定通过可控阻尼减振器调节悬架阻尼系数,从而使汽车获得良好的平顺性和操纵稳定性。
本发明的关键是运用小波变换的方法来分解振动信号,结合图2,其步骤如下:
1、设:x(n)表示滤波器在n时刻的输入,X0=[x(n),x(n-1),…x(1)]T表示由x(n)形成的输入向量;DWT表示振动信号的离散小波变换分解;Xj表示信号序列X0的第j级分解后的低频信号序列;Dj表示信号序列X0的第j级高频信号序列;σj是第j级细节信号序列的均方根值,σj+1是第j级逼近信号序列的均方根值,均方根值向量:P=[σ1,σ2,...,σj+1];wj是第j级细节信号序列对应的输出加权系数,wj+1是第j级逼近信号序列对应的输出加权系数,输出权向量:W=[w1,w2,…,wj+1]T
控制器输出量F的计算:
F=kPW (1)
式中k是输出增益。
2、振动信号的小波分解,结合图3:
设:G和H分别为小波低通滤波器和小波高通滤波器,Xj和Dj是第j级分解的逼近信号(低频)和细节信号(高频)变换公式如下:
其中h(n)和g(n)表示这滤波器H和G的冲击响应序列,并且h(n)和g(n)有如下的关系:
g(n)=(-1)1-nh(1-n) (4)
为了提高算法的实时性,选着具有对称性的小波滤波器h(n),使h(n)=h(-n)
设l=k-2m,k=2m+l,若冲击响应序列h(n)的长度为2t+1,
-(t-1)/2≤l≤(t-1)/2,所以,式(1)可以写成:
把式(5)展开成双边的形式:
因为h(l)=h(-l),所以:
同理,式(3)写成如下形式
因此利用滤波器H和滤波器G的关系以及滤波器H的对称性,可以对信号分解计算公式进行简化,减少部分计算量,提高算法的实时性。
3、输出权向量的调整:
输出权向量经离线整定后,输出权向量W一般不需要调整,当控制效果欠佳时,可以通过输出权向量调节器对进行调整。调整算法如下:
设第m次控制器输出为Fm,相应的输出权向量和均方根值向量为Wm=[wm,1,wm,2,...,wm,j+1]和Pm=[Pm,1,Pm,2,...,Pm,j+1]。
构造误差向量:Em=Pm-Pm-1=[em,1,em,2,...,em,j+1];
输出权重调节向量:ΔW=[Δwm,1,Δwm,2,...,Δwm,j+1];
式中ε为控制阈值,r根据所在频段取0或1,衰减系数R∈(0,0.618)。
Wm+1=Wm+ΔWm (10)
4.通过小波分解,振动信号就可以分解为由高到低的若干频段,通过能量统计的方法确定当前振动信号主要分布在哪几个频段,具体步骤如下:
(1)对小波分解得到的各级振动信号进行能量统计,对于j级分解后的信号向量Di=[di(1),...,di(n-1),di(n)]T,(i=1,2,...,j,j+1;Dj+1=Xi),求各向量的均方根值
(2)建立均方根值向量:P=[σ1,σ2,...,σj+1]中元素排序,选取最大的两个元素σm和σi
(3)根据小波分解信号与频率的对应关系,确定σm和σi所在频段的频率范围。如果原始信号的最高频率为fmax,可以计算出7级小波分解后各层所占的具体频带为:
X7:(0~0.0078125)×fmax
D7:(0.0078125~0.015625)×fmax
D6:(0.015625~0.03125)×fmax
D5:(0.03125~0.0625)×fmax
D4:(0.0625~0.125)×fmax
D3:(0.125~0.25)×fmax
D2:(0.25~0.5)×fmax
D1:(0.5~1)×fmax
本发明的优点:
本控制方法运用小波变换的方法将振动信号分解为由高到低的若干频段,对车身振动信号进行多分辨率分析,高频信号部分分辨率低,低频信号分辨高,符合汽车振动信号的特点。由于算法只需要利用加速度传感器在线提取汽车行驶中车身的振动信号,无须预测路面,可极大地降低***成本,提高***可靠性。通过能量统计的方法确定当前振动信号主要分布在哪几个频段,根据阻尼系数与车轴、车身振动响应曲线确定通过可控阻尼减振器调节悬架阻尼系数,兼顾悬架低频振动和高频振动,从而使汽车获得良好的平顺性和操纵稳定性。本控制方法实现容易,能够方便地应用到现有汽车悬架上,实现实时控制。
附图说明
图1是采用小波频域控制方法的汽车磁流变半主动悬架***示意图;
图2是小波频域控制算法的流程框图;
图3是信号多级分解的流程图;
图4悬架响应的频段划分,其中图4A是簧载质量加速度传递特性,图4B是轮胎动载荷传递特性;
图5振动信号的多级分解;
图6车厢地板垂直振动控制效果(40km/h);
图7车轴垂直振动控制效果(40km/h)。
具体实施方式
如图1所示,将加速度传感器1安装于车身质心处,采集车身垂直振动加速度信号,作为小波频域控制器6的输入,控制器对输入信号进行小波分解后,参见图2和图3,通过能量统计的方法确定当前振动信号主要分布在哪几个频段,根据阻尼系数与车轴、车身振动响应曲线确定通过可控阻尼减振器调节悬架阻尼系数,输出控制信号调节可控阻尼的执行器(如电/磁流变阻尼器)2-5的阻尼特性,完成一个控制循环。
图5给出了汽车以40km/h速度行驶时,车厢地板垂直振动的加速度信号的7级分解结果,其中X0是传感器采集到的原始信号,X7是第7级逼近信号(低频部分),D1~D7是1~7级细节信号(高频部分),从图中可以看到,各层小波变换如同带通或低通滤波器,通过小波变换,各个频段的信号被清晰地展现出来。原始信号的最高频率为50Hz,可以计算出各层所占的具体频带为:
X7:0~0.39Hz
D7:0.39~0.78Hz
D6:0.78~1.56Hz
D5:1.56~3.13Hz
D4:3.125~6.25Hz
D3:6.25~12.5Hz
D2:12.5~25Hz
D1:25~50Hz
通过信号能量统计,各级振动信号的均方根值为:P=[σ1,σ2,...,σj+1]=[0.002,0.005,0.004,0.01,0.09,0.23,0.16,0.03]。由此可知振动能量集中在D6:0.78~1.56Hz和D7:0.39~0.78Hz这两个频段。因此输出阻尼力应该取较小值,设输出增益k=1000;
根据图4所示的悬架响应的频段划分可知,在[0~1.1Hz]频段应该为低阻尼状态,[1.1~9Hz]频段应该为高阻尼状态,9~100Hz频段应该为中等阻尼状态,因此可以设置输出权向量为W=[0,0,0.5,1,1,0.8,0.5,0.5]。
输出阻尼力F=k*P*W=380N
图6显示了汽车以40km/h速度在B级路面上的行驶时的控制效果。图中黑线为被动悬架的簧载质量-车厢地板垂直加速度功率谱,红线为采用基于小波分解的汽车半主动悬架阻尼力控制方法后测得车厢地板垂直加速度功率谱。比较得知,本专利所采用的控制方法能有效降低人体敏感的4-12.5Hz的垂直加速度,提高汽车乘座舒适性。图7显示了非簧载质量-车轴的垂直加速度功率谱,从图中可以看出,通过控制,车轴的垂直振动也得到了较好抑制,在提高汽车乘坐舒适性的同时也提高了操纵稳定性。
Claims (2)
1、一种基于小波分解的汽车半主动悬架阻尼控制的方法,其步骤如下:
(1)利用加速度传感器在线提取汽车行驶中车身的振动信号;
(2)运用小波变换的方法将振动信号分解为由高到低的若干频段;
(3)通过能量统计的方法确定当前振动信号主要分布在哪几个频段;
(4)根据阻尼系数与车轴、车身振动响应曲线确定通过可控阻尼减振器调节悬架阻尼系数,从而使汽车获得良好的平顺性和操纵稳定性;
其中:运用小波变换的方法分解振动信号的步骤如下:
(1)设:x(n)表示滤波器在n时刻的输入,X0=[x(n),x(n-1),…x(1)]T表示由x(n)形成的输入向量;DWT表示振动信号的离散小波变换分解;Xj表示信号序列X0的第j级分解后的低频信号序列;Dj表示信号序列X0的第j级高频信号序列;σj是第j级高频信号序列的均方根值,σj+1是第j级低频信号序列的均方根值,均方根值向量:P=[σ1,σ2,...,σj+1];wj是第j级细节信号序列对应的输出加权系数,wj+1是第j级逼近信号序列对应的输出加权系数,输出权向量:W=[w1,w2,…,wj+1]T
控制器输出量F的计算:
F=kPW (1)
式中k是输出增益;
(2)振动信号的小波分解:
以G和H分别代表小波低通滤波器和小波高通滤波器,Xj和Dj是第j级分解的低频逼近信号和高频细节信号变换公式如下:
其中h(n)和g(n)表示这滤波器H和G的冲击响应序列,并且h(n)和g(n)有如下的关系:
g(n)=(-1)1-nh(1-n) (4)
为了提高算法的实时性,选着具有对称性的小波滤波器h(n),使h(n)=h(-n)
设l=k-2m,k=2m+l,若冲击响应序列h(n)的长度为2t+1,-(t-1)/2≤l≤(t-1)/2,所以,式(1)可以写成:
把式(5)展开成双边的形式:
因为h(l)=h(-l),所以:
同理,式(3)写成如下形式
(3)输出权向量的调节:
输出权向量经离线整定后,输出权向量W一般不需要调整,当控制效果欠佳时,可以通过输出权向量调节器对其进行调节,调节方法如下:
设第m次控制器输出为Fm,相应的输出权向量和均方根值向量为Wm=[wm,1,wm,2 ,...,wm,j+1]和Pm=[Pm,1,Pm,2,...,Pm,j+1];
构造误差向量:Em=Pm-Pm-1=[em,1,em,2,...,em,j+1];
输出权重调节向量:ΔW=[Δwm,1,Δwm,2,...,Δwm,j+1];
式中ε为控制阈值,r根据所在频段取0或1,衰减系数R∈(0,0.618)
Wm+1=Wm+ΔWm (10);
所述步骤(3)通过能量统计的方法确定当前振动信号主要分布在哪几个频段的具体步骤如下:
A、对于j级分解后的信号向量Di=[di(1),…,di(n-1),di(n)]T,(i=1,2,...,j,j+1;Dj+1=Xi),求各向量的均方根值
(i=1,2,...,j,j+1);
B、均方根值向量:P=[σ1,σ2,...,σj+1]中元素排序,选取最大的两个元素σm和σl;
C、根据小波分解信号与频率的对应关系,确定σm和σl所在频段的频率范围,如果原信号的最高频率为fmax,计算出7级小波分解后各层所占的具体频带为:
X7:(0~0.0078125)×fmax
D7:(0.0078125~0.015625)×fmax
D6:(0.015625~0.03125)×fmax
D5:(0.03125~0.0625)×fmax
D4:(0.0625~0.125)×fmax
D3:(0.125~0.25)×fmax
D2:(0.25~0.5)×fmax
D1:(0.5~1)×fmax。
2.根据权利要求1所述的方法,其特征在于:根据阻尼系数与车轴、车身振动响应曲线确定通过可控阻尼减振器调节悬架阻尼系数具体是,根据不同阻尼系数下簧载质量加速度相对于路面速度输入的传递函数幅频特性曲线的交点和轮胎动载荷对于路面速度输入的传递函数幅频特性曲线的交点,将0.1-100Hz的频率划为由低到高的三个频段A、B、C,在频段A使用大阻尼系数,在频段B使用小阻尼系数,对于频段C,选用中等阻尼系数,使簧载质量加速度和轮胎动载荷的均方根值都保持在一个较低的水平,A、B、C三个频率点的值应结合具体车型测试得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100401179A CN100358739C (zh) | 2004-06-29 | 2004-06-29 | 一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100401179A CN100358739C (zh) | 2004-06-29 | 2004-06-29 | 一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1593961A CN1593961A (zh) | 2005-03-16 |
CN100358739C true CN100358739C (zh) | 2008-01-02 |
Family
ID=34664478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100401179A Expired - Fee Related CN100358739C (zh) | 2004-06-29 | 2004-06-29 | 一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100358739C (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101269618B (zh) * | 2008-04-25 | 2010-06-02 | 江苏大学 | 一种电子控制空气悬架三档可调减振器阻尼值的控制方法 |
TWI630128B (zh) * | 2017-06-06 | 2018-07-21 | 樹德科技大學 | 可自動調整之避震裝置及其方法 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITMI20061403A1 (it) * | 2006-07-19 | 2008-01-20 | Milano Politecnico | Metodo ed apparato per controllare una sospensione semiattiva |
US8914193B2 (en) * | 2009-02-23 | 2014-12-16 | Toyota Jidosha Kabushiki Kaisha | Damping force control apparatus |
KR101360050B1 (ko) * | 2012-07-30 | 2014-02-10 | 기아자동차주식회사 | 차량의 경사각 측정 장치 및 방법 |
DE102013217870B4 (de) * | 2013-09-06 | 2022-10-06 | Volkswagen Aktiengesellschaft | Verfahren und Vorrichtung zum Betreiben eines Dämpfungssystems für ein Kraftfahrzeug |
CN104972860B (zh) * | 2015-01-23 | 2017-04-05 | 江苏大学 | 一种用于车辆主动悬架的全息最优滑模控制器 |
CN105138783B (zh) * | 2015-09-06 | 2018-01-16 | 山东理工大学 | 高速轨道车辆车体端部横向减振器阻尼系数的设计方法 |
CN105974821B (zh) * | 2016-05-16 | 2019-01-18 | 萨克斯汽车零部件***(上海)有限公司 | 基于阻尼多模式切换减振器的车辆半主动悬架混杂控制方法 |
CN106202653B (zh) * | 2016-06-28 | 2019-11-05 | 广州汽车集团股份有限公司 | 一种车身阻尼分布优化方法和*** |
CN107575529B (zh) * | 2017-10-16 | 2019-06-11 | 中车株洲电力机车有限公司 | 一种轨道交通车辆及其自适应减振方法和*** |
JP6939660B2 (ja) * | 2018-03-13 | 2021-09-22 | トヨタ自動車株式会社 | 車両走行制御システム |
CN113183711B (zh) * | 2021-06-04 | 2022-03-15 | 合肥工业大学 | 一种振动与冲击融合控制方法 |
CN113665311B (zh) * | 2021-07-16 | 2024-02-20 | 中国北方车辆研究所 | 一种基于频域分析的减振器控制方法及*** |
CN117445780B (zh) * | 2023-12-26 | 2024-02-23 | 常熟理工学院 | 一种变刚度变阻尼汽车座椅智能控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1121566A (zh) * | 1994-07-28 | 1996-05-01 | 华中理工大学 | 电流变流体阻尼装置振动传递率的半主动控制方法 |
CN1374207A (zh) * | 2001-03-08 | 2002-10-16 | 北京理工大学 | 半主动控制式车辆动力总成磁流变体液压悬置 |
-
2004
- 2004-06-29 CN CNB2004100401179A patent/CN100358739C/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1121566A (zh) * | 1994-07-28 | 1996-05-01 | 华中理工大学 | 电流变流体阻尼装置振动传递率的半主动控制方法 |
CN1374207A (zh) * | 2001-03-08 | 2002-10-16 | 北京理工大学 | 半主动控制式车辆动力总成磁流变体液压悬置 |
Non-Patent Citations (3)
Title |
---|
基于小波理论的复杂机械振动信号降噪分析. 杨文平,陈国定,石博强.北京科技大学学报,第24卷第4期. 2002 * |
汽车半主动悬架的多目标优化设计. 陈大跃.机械工程学报,第36卷第6期. 2000 * |
车辆悬架振动控制***研究的进展. 张庙康,胡海岩.振动、测试与诊断,第17卷第1期. 1997 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101269618B (zh) * | 2008-04-25 | 2010-06-02 | 江苏大学 | 一种电子控制空气悬架三档可调减振器阻尼值的控制方法 |
TWI630128B (zh) * | 2017-06-06 | 2018-07-21 | 樹德科技大學 | 可自動調整之避震裝置及其方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1593961A (zh) | 2005-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100358739C (zh) | 一种基于小波分解的汽车磁流变半主动悬架阻尼控制的方法 | |
CN100484788C (zh) | 汽车磁流变半主动悬架分姿态协调仿人智能控制方法 | |
CN106985627B (zh) | 一种车辆路面识别***以及悬架模式切换方法 | |
CN101367324B (zh) | 一种空气悬架车辆高度传感器的路面等级预测方法 | |
Sugahara et al. | Suppression of vertical bending and rigid-body-mode vibration in railway vehicle car body by primary and secondary suspension control: results of simulations and running tests using Shinkansen vehicle | |
JPS60248417A (ja) | アクテイブサスペンシヨン装置 | |
KR19980075986A (ko) | 차량용 현가장치의 제어방법 | |
CN110641502B (zh) | 列车转向架悬挂***动态调节方法 | |
CN102189909A (zh) | 一种车辆半主动悬架天棚阻尼频率滤波控制策略 | |
KR19980031243A (ko) | 차량의 반 능동 전자 제어 현가 장치 및 방법 | |
CN103754081B (zh) | 车辆非线性悬架***的最优模糊复合控制方法 | |
CN110341414B (zh) | 一种连续线性天棚控制下悬架自适应最优控制***及方法 | |
CN114683795A (zh) | 一种路面自适应半主动悬架控制方法和*** | |
Shalabi et al. | Control of automotive air-spring suspension system using Z-number based fuzzy system | |
CN113635726B (zh) | 整车半主动悬架***集成控制方法和*** | |
JP3128006B2 (ja) | 車両用サスペンションの制御方法及び装置 | |
Sugahara et al. | Suppressing vertical vibration in railway vehicles through air spring damping control | |
CN113665311B (zh) | 一种基于频域分析的减振器控制方法及*** | |
JP3347169B2 (ja) | 車両用サスペンションの制御装置 | |
Abouel-seoud | Vibration Behaviour Control of a Fabricated One-Passenger Electric Vehicle with Either Mechanical or Air Suspension. | |
Zhu et al. | Optimal damping matching of semi-active air suspension based on vehicle comfort and handling | |
Montazeri-Gh et al. | Vehicle ride evaluation based on a time-domain variable speed driving pattern | |
Soliman et al. | Effect of suspension spring stiffness on vehicle dynamics | |
Li et al. | A hierarchical controller for the vibration of an automotive suspension system via magnetorheological dampers | |
Tudon-Martinez et al. | Efficiency of on-off semiactive suspensions in a pick-up truck |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080102 Termination date: 20100629 |