CN104965945A - 电动汽车减振方法 - Google Patents

电动汽车减振方法 Download PDF

Info

Publication number
CN104965945A
CN104965945A CN201510364651.3A CN201510364651A CN104965945A CN 104965945 A CN104965945 A CN 104965945A CN 201510364651 A CN201510364651 A CN 201510364651A CN 104965945 A CN104965945 A CN 104965945A
Authority
CN
China
Prior art keywords
electric vehicle
damper
objective function
spring
parameters
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.)
Granted
Application number
CN201510364651.3A
Other languages
English (en)
Other versions
CN104965945B (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.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
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 Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201510364651.3A priority Critical patent/CN104965945B/zh
Publication of CN104965945A publication Critical patent/CN104965945A/zh
Application granted granted Critical
Publication of CN104965945B publication Critical patent/CN104965945B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Vehicle Body Suspensions (AREA)
  • Arrangement Or Mounting Of Propulsion Units For Vehicles (AREA)

Abstract

本发明提供一种电动汽车减振方法,特征是:电动汽车的前、后电池组与车架之间均通过弹簧和阻尼器连接,弹簧和阻尼器参数的优化利用Matlab编程实现,步骤如下:1)给出电动汽车半车模型的质量矩阵、阻尼矩阵和刚度矩阵;2)构造目标函数,求取目标函数对弹簧和阻尼器参数的偏导函数表达式;3)给出弹簧和阻尼器参数的初值,计算此时目标函数s的值,记为sa;4)给出搜索步长,计算迭代后弹簧和阻尼器的参数值以及此时目标函数s的值,记为sb;5)若sa≥sb,则令sa=sb,返回步骤4)继续循环,否则结束循环,得到弹簧和阻尼器的优化参数值。本发明用优化选择的弹簧和阻尼器连接车架和电池组,为电动汽车车身减振,减振效果优良。

Description

电动汽车减振方法
技术领域
本发明涉及一种电动汽车减振方法,更确切地说是一种利用动力吸振器思想提高车身减振效果的电动汽车方法。
背景技术
近年来,电动汽车技术快速发展,需求量和产量都日益增加,人们对电动汽车的乘坐舒适性要求也越来越高,这就要求提高电动汽车的减振效果。为提高电动汽车的减振效果,目前已有的技术方案主要是:设置电动汽车驱动电机的悬置装置或是电磁悬架来优化电动汽车的减振效果,提高电动汽车的乘坐舒适性;但是,由于驱动电机要与传动系统相连,故用其来吸收电动汽车车身的振动有一定的限制,而电磁悬架要增加很大的质量,这就影响了电动汽车动力性等性能。可见已有的电动汽车减振方法或是减振效果不够理想,或是虽有较好的减振效果但降低了其他性能,因此电动汽车减振方法还有待于进一步改进。
动力吸振器是较为常用的减振手段之一,但是利用动力吸振器来减振需要一定的附加质量,而为汽车减振需要的附加质量是非常大的。在电动汽车上增加很大的附加质量会使得电动汽车的动力性等其它性能变差,也与汽车轻量化的发展趋势背道而驰。
发明内容
本发明的目的是提供一种能克服上述缺陷、在不过多增加汽车质量的前提下提高车身减振效果的电动汽车减振方法。其技术方案为:
一种电动汽车减振方法,其特征在于:电动汽车的前端电池组的底部与电动汽车车架之间通过前端弹簧和前端阻尼器连接,后端电池组的底部与电动汽车车架之间通过后端弹簧和后端阻尼器连接,前端弹簧、后端弹簧、前端阻尼器和后端阻尼器的参数的优化利用Matlab编程实现,具体实施步骤如下:
1)给出电动汽车半车模型的质量矩阵M、阻尼矩阵C和刚度矩阵K,具体如下:
M = m t f 0 0 0 0 0 0 m t r 0 0 0 0 0 0 m c 0 0 0 0 0 0 I c 0 0 0 0 0 0 m v f 0 0 0 0 0 0 m v r
其中mtf和mtr分别是电动汽车两个前轮和两个后轮的质量,mc是电动汽车车身质量,Ic是电动汽车车身对于质心的转动惯量,mvf和mvr分别是前端电池组和后端电池组的质量;
C = c s f 0 - c s f - 1 f c s f 0 0 0 c s r - c s r - 1 r c s r 0 0 - c s f - c s r ( c s f + c s r + c v f + c v r ) / 2 0 - c v f - c v r - 1 f c s f - 1 r c s r 0 ( 1 f ( c s f + c v f ) + 1 r ( c s r + c v r ) ) / 2 - 1 f c v f - 1 r c v r 0 0 - c v f - 1 f c v f c v f 0 0 0 - c v r - 1 r c v r 0 c v r
其中csf和csr分别是电动汽车前、后悬架的阻尼器系数,lf和lr分别是电动汽车前桥、后桥与质心的水平距离,以上四个参数均为电动汽车的固定参数,cvf是待优化的前端阻尼器的参数,cvr是待优化的后端阻尼器的参数;
K = k t f - k s f 0 - k s f - 1 f c s f 0 0 0 k t r - k s r - k s r - 1 r c s r 0 0 - k s f - k s r ( k s f + k s r + k v f + k v r ) / 2 0 - k v f - k v r - 1 f k s f - 1 r k s r 0 ( 1 f ( k s f + k v f ) + 1 r ( k s r + k v r ) ) / 2 - 1 f k v f - 1 r k v r 0 0 - k v f - 1 f k v f k v f 0 0 0 - k v r - 1 r k v r 0 k v r
其中ktf和ktr分别是电动汽车前、后轮胎的等效刚度,ksf和ksr分别是是电动汽车前、后悬架的刚度,kvf是待优化的前端弹簧的参数,kvr是待优化的后端弹簧的参数;
2)构造目标函数,求取目标函数对于步骤1)中的四个待优化减振参数的偏导函数表达式,具体的指令如下:
w=6.5;
B=K+i*w*C-w*w*M;
F=[1 1 0 0 0 0]';
X=inv(B)*F;
s=0.7*X(3)+0.3*X(4);          %目标函数%
s1=diff(s,kvf);               %目标函数对于kvf的偏导函数%
s2=diff(s,cvf);               %目标函数对于cvf的偏导函数%
s3=diff(s,kvr);               %目标函数对于kvr的偏导函数%
s4=diff(s,cvr);               %目标函数对于cvr的偏导函数%
3)给出步骤1)中四个待优化的参数kvf、cvf、kvr和cvr的初值,计算此时目标函数s的值,记为sa,具体指令为:
sa=abs(eval(s));
4)给出迭代算法的搜索步长h,计算迭代后的四个变量kvf、cvf、kvr和cvr的值,并计算此时目标函数s的值,记为sb,具体指令为:
kvf=kvf+abs(s1)*h;
cvf=cvf+abs(s2)*h;
kvr=kvr+abs(s3)*h;
cvr=cvr+abs(s4)*h;
sb=abs(eval(s));
5)比较sa和sb值的大小:若sa≥sb,则令sa=sb,返回步骤4)继续循环;若sa<sb,结束循环,此时步骤4)中kvf、cvf、kvr和cvr的值即为所求弹簧和阻尼器的参数。
本发明与现有技术相比,其优点在于:
(1)利用动力吸振器思想,用弹簧和阻尼器连接电动汽车车架和电池组,以电池组为调谐质量为电动汽车车身减振,在没有过多增加电动汽车上质量的前提下,提高了车身的减振效果;
(2)利用最优化方法,设计出最优的弹簧刚度参数和阻尼器参数,使得电动汽车车身减振效果达到最理想的状态;
附图说明
图1为电动汽车结构示意图。
图2为本发明电池组悬置前、后车身位移的时域响应曲线图;
图3为本发明电池组悬置前、后车身加速度的时域响应曲线图;
图中:1、车架  2、前端电池组  3、前端弹簧  4、前端阻尼器  5、后端电池组  6、后端弹簧  7、后端阻尼器   A、电池组悬置前车身位移的时域响应曲线   B、电池组悬置后车身位移的时域响应曲线  C、电池组悬置前车身加速度的时域响应曲线  D、电池组悬置后车身加速度的时域响应曲线。
具体实施方式
下面结合实施例及附图对本发明做进一步说明:电动汽车车身质量mc=640Kg,车身对于质心的转动惯量Ic=1200Kg·m2,质心与前桥距离lf=1.3m,质心与后桥距离lr=1.5m,前悬架刚度ksf=34000N/m,前悬架阻尼系数csf=3000N·s/m,后悬架刚度ksr=28000N/m,后悬架阻尼系数csr=3000N·s/m,前轮质量mtf=80Kg,后轮质量mtr=80Kg,前轮等效刚度ktf=350000N/m,后轮等效刚度ktr=350000N/m,前端电池组质量mvf=340Kg,后端电池组质量mvr=120Kg。前端电池组2的底部与车架1之间通过前端弹簧3和前端阻尼器4连接,后端电池组5的底部与车架1之间通过后端弹簧6和后端阻尼器7连接,其中前端电池组2的质心与前桥垂直方向对齐,后端电池组5的质心与后桥垂直方向对齐,前端弹簧3、后端弹簧6、前端阻尼器4和后端阻尼器7参数的优化利用Matlab编程实现,具体实施步骤如下:
1)给出电动汽车半车模型的质量矩阵M、阻尼矩阵C和刚度矩阵K,具体如下:
M = 80 0 0 0 0 0 0 80 0 0 0 0 0 0 640 0 0 0 0 0 0 1200 0 0 0 0 0 0 240 0 0 0 0 0 0 120
C = 3000 0 - 3000 - 3900 0 0 0 3000 - 3000 - 4500 0 0 - 3000 - 3000 ( c v f + c v r ) / 2 + 3000 0 - c v f - c v r - 3900 - 4500 0 ( 1.3 c v f + 1.5 c v r ) - 1.3 v v f - 0.5 c v r 0 0 - c v f - 1.3 c v f c v f 0 0 0 - c v r - 1.5 c v r 0 c v r
其中cvf是待优化的前端阻尼器4的参数,cvr是待优化的后端阻尼器7的参数;
K = 316000 0 - 34000 - 44200 0 0 0 322000 - 28000 - 42000 0 0 - 3400 - 28000 ( k v f + k v r ) / 2 + 31000 0 - k v f - k v r - 44200 - 42000 0 ( 1.3 k v f + 1.5 k v r ) / 2 + 43100 - 1.3 k v f - 0.5 k v r 0 0 - k v f - 1.3 k v f k v f 0 0 0 - k v r - 1.5 k v r 0 k v r
其中kvf是待优化的前端弹簧3的参数,kvr是待优化的后端弹簧6的参数。
2)构造目标函数,求取目标函数对于步骤1)中的四个待优化减振参数的偏导函数表达式,具体的指令如下:
w=6.5;
B=K+i*w*C-w*w*M;
F=[1 1 0 0 0 0]';
X=inv(B)*F;
s=0.7*X(3)+0.3*X(4);         %目标函数%
s1=diff(s,kvf);              %目标函数对于kvf的偏导函数%
s2=diff(s,cvf);               %目标函数对于cvf的偏导函数%
s3=diff(s,kvr);               %目标函数对于kvr的偏导函数%
s4=diff(s,cvr);               %目标函数对于cvr的偏导函数%
3)给出四个待优化的电池悬置参数kvf、cvf、kvr、cvr的初值,具体为kvf=8000N/m,cvf=700N·s/m,kvr=5000N/m,cvr=400N·s/m,计算此时目标函数s的值,记为sa,具体指令为:
sa=abs(eval(s));
4)给出迭代算法的搜索步长h=106,计算迭代后的四个变量kvf、cvf、kvr、cvr的值,并计算此时目标函数s的值,记为sb,具体指令为:
kvf=kvf+abs(s1)*h;
cvf=cvf+abs(s2)*h;
kvr=kvr+abs(s3)*h;
cvr=cvr+abs(s4)*h;
sb=abs(eval(s));
5)比较sa和sb值的大小:若sa≥sb,则令sa=sb,返回步骤4)继续循环;若sa<sb,结束循环,输出此时4)中kvf、cvf、kvr、cvr的值即为所求弹簧和阻尼器的参数。
在MATLAB中编写程序得到优化结果是:kvf=12067N/m,cvf=1523N·s/m,kvr=8077N/m,cvr=1569N·s/m,适当圆整后四个参数值分别取为:kvf=12000N/m,cvf=1500N·s/m,kvr=8000N/m,cvr=1500N·s/m。
得出上述4个参数后,根据参数选定前端弹簧3、前端阻尼器4、后端弹簧6、后端阻尼器7,对应安装在前端电池组2、后端电池组5与车架1之间,比较电池组悬置前、后车身位移的时域响应曲线图以及车身加速度的时域响应曲线图,可以看出利用本发明所述的方法,电动汽车车身的减振效果有明显的提高。

Claims (1)

1.一种电动汽车减振方法,其特征在于:电动汽车的前端电池组(2)的底部与车架(1)之间通过前端弹簧(3)和前端阻尼器(4)连接,后端电池组(5)的底部与车架(1)之间通过后端弹簧(6)和后端阻尼器(7)连接,前端弹簧(3)、后端弹簧(6)、前端阻尼器(4)和后端阻尼器(7)的参数的优化利用Matlab编程实现,具体实施步骤如下:
1)给出电动汽车半车模型的质量矩阵M、阻尼矩阵C和刚度矩阵K,具体如下:
M = m t f 0 0 0 0 0 0 m t r 0 0 0 0 0 0 m c 0 0 0 0 0 0 I c 0 0 0 0 0 0 m v f 0 0 0 0 0 0 m v r
其中mtf和mtr分别是电动汽车两个前轮和两个后轮的质量,mc是电动汽车车身质量,Ic是电动汽车车身对于质心的转动惯量,mvf和mvr分别是前端电池组和后端电池组的质量;
C = c s f 0 - c s f - l f c s f 0 0 0 c s r - c s r - l r c s r 0 0 - c s f - c s r ( c s f + c s r + c v f + c v r ) / 2 0 - c v f - c v r - l f c s f - l r c s r 0 ( l f ( c s f + c v f ) + l r ( c s r + c v r ) ) / 2 - l f c v f - l r c v r 0 0 - c v f - l f c v f c v f 0 0 0 - c v r - l r c v r 0 c v r
其中csf和csr分别是电动汽车前、后悬架的阻尼器系数,lf和lr分别是电动汽车前桥、后桥与质心的水平距离,以上四个参数均为电动汽车的固定参数,cvf是待优化的前端阻尼器(4)的参数,cvr是待优化的后端阻尼器(7)的参数;
K = k t f - k s f 0 - k s f - l f k s f 0 0 0 k t r - k s r - k s r - l r k s r 0 0 - k s f - k s r ( k s f + k s r + k v f + k v r ) / 2 0 - k v f - k v r - l f c s f - l r c s r 0 ( l f ( k s f + k v f ) + l r ( k s r + k v r ) ) / 2 - l f k v f - l r k v r 0 0 - k v f - l f k v f k v f 0 0 0 - k v r - l r k v r 0 k v r
其中ktf和ktr分别是电动汽车前、后轮胎的等效刚度,ksf和ksr分别是是电动汽车前、后悬架的刚度,kvf是待优化的前端弹簧(3)的参数,kvr是待优化的后端弹簧(6)的参数;
2)构造目标函数,求取目标函数对于步骤1)中的四个待优化减振参数的偏导函数表达式,具体的指令如下:
w=6.5;
B=K+i*w*C-w*w*M;
F=[1 1 0 0 0 0]';
X=inv(B)*F;
s=0.7*X(3)+0.3*X(4); %目标函数%
s1=diff(s,kvf);       %目标函数对于kvf的偏导函数%
s2=diff(s,cvf);       %目标函数对于cvf的偏导函数%
s3=diff(s,kvr);       %目标函数对于kvr的偏导函数%
s4=diff(s,cvr);       %目标函数对于cvr的偏导函数%
3)给出步骤1)中四个待优化的参数kvf、cvf、kvr和cvr的初值,计算此时目标函数s的值,记为sa,具体指令为:
sa=abs(eval(s));
4)给出迭代算法的搜索步长h,计算迭代后的四个变量kvf、cvf、kvr和cvr的值,并计算此时目标函数s的值,记为sb,具体指令为:
kvf=kvf+abs(s1)*h;
cvf=cvf+abs(s2)*h;
kvr=kvr+abs(s3)*h;
cvr=cvr+abs(s4)*h;
sb=abs(eval(s));
5)比较sa和sb值的大小:若sa≥sb,则令sa=sb,返回步骤4)继续循环;若sa<sb,结束循环,此时步骤4)中kvf、cvf、kvr和cvr的值即为所求弹簧和阻尼器的参数。
CN201510364651.3A 2015-06-29 2015-06-29 电动汽车减振方法 Expired - Fee Related CN104965945B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510364651.3A CN104965945B (zh) 2015-06-29 2015-06-29 电动汽车减振方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510364651.3A CN104965945B (zh) 2015-06-29 2015-06-29 电动汽车减振方法

Publications (2)

Publication Number Publication Date
CN104965945A true CN104965945A (zh) 2015-10-07
CN104965945B CN104965945B (zh) 2018-03-27

Family

ID=54219982

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510364651.3A Expired - Fee Related CN104965945B (zh) 2015-06-29 2015-06-29 电动汽车减振方法

Country Status (1)

Country Link
CN (1) CN104965945B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112659875A (zh) * 2020-12-28 2021-04-16 奇瑞汽车股份有限公司 一种提升整车平顺性的电池悬置和设计方法
CN115230421A (zh) * 2022-07-22 2022-10-25 山东理工大学 基于车辆主动悬架振动抑制及能量收集的时滞控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076148A (zh) * 2013-01-04 2013-05-01 常州万安汽车部件科技有限公司 基于跌落实验双-4自由度半车模型的车辆参数识别方法
CN104385873A (zh) * 2014-09-24 2015-03-04 湖南大学 一种汽车悬架系统的多目标优化方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076148A (zh) * 2013-01-04 2013-05-01 常州万安汽车部件科技有限公司 基于跌落实验双-4自由度半车模型的车辆参数识别方法
CN104385873A (zh) * 2014-09-24 2015-03-04 湖南大学 一种汽车悬架系统的多目标优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
金明 等: "基于遗传算法的某微型纯电动汽车平顺性优化", 《计算机仿真》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112659875A (zh) * 2020-12-28 2021-04-16 奇瑞汽车股份有限公司 一种提升整车平顺性的电池悬置和设计方法
CN115230421A (zh) * 2022-07-22 2022-10-25 山东理工大学 基于车辆主动悬架振动抑制及能量收集的时滞控制方法
CN115230421B (zh) * 2022-07-22 2024-08-20 山东理工大学 基于车辆主动悬架振动抑制及能量收集的时滞控制方法

Also Published As

Publication number Publication date
CN104965945B (zh) 2018-03-27

Similar Documents

Publication Publication Date Title
CN103121475B (zh) 驾驶室悬置最优阻尼比的设计方法
CN101327722B (zh) 应用惯性蓄能器的车辆悬架
CN105069261B (zh) 低速轨道车辆二系横向减振器最优阻尼系数的设计方法
CN103593506B (zh) 一种两级串联式isd悬架参数优化方法
CN106427555B (zh) 一种用于电动汽车车轮全方位转向的转换装置及可全方位转向的电动汽车
CN105159094A (zh) 汽车主动悬架lqg控制器最优控制力的设计方法
CN105172511A (zh) 一种汽车悬架减振器控制系统及方法
CN103112508B (zh) 卡车驾驶室减振器最佳速度特性的设计方法
Chen et al. Approaches to diminish large unsprung mass negative effects of wheel side drive electric vehicles
CN102092256B (zh) 减小单纵臂悬架轮边电驱动系统等效簧下质量结构及方法
CN204641300U (zh) 导出式多承点独立悬架
CN104965945B (zh) 电动汽车减振方法
CN201610056U (zh) 一种车辆独立悬架
Purushotham Comparative simulation studies on MacPherson suspension system
CN202727902U (zh) 电动汽车使用的前后四轮独立悬架装置
CN203391592U (zh) 一种汽车用负刚度悬架机构及具有该机构的车辆
CN204263868U (zh) 卧式独立悬挂装置
CN102092255B (zh) 减小单斜臂悬架轮边电驱动系统等效簧下质量结构及方法
CN105302944B (zh) 基于广义耗散能的车辆被动悬架最佳阻尼比的计算方法
CN104608583A (zh) 一种汽车自抗扰悬挂系统与控制方法
Vaishnav et al. Mathematical modelling and comparison of two degree of freedom suspension system of quarter car
CN208381190U (zh) 一种曲式复合空气弹簧减震装置
CN103486174B (zh) 一种减震弹簧
Gagliardi et al. Performance analysis of linear and rotary energy harvesting shock absorber systems
CN115570923B (zh) 一种轮毂电机驱动车辆的减振装置

Legal Events

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

Granted publication date: 20180327