CN109271680A - 一种将1/4汽车模型等效为二自由度模型的方法 - Google Patents
一种将1/4汽车模型等效为二自由度模型的方法 Download PDFInfo
- Publication number
- CN109271680A CN109271680A CN201810989115.6A CN201810989115A CN109271680A CN 109271680 A CN109271680 A CN 109271680A CN 201810989115 A CN201810989115 A CN 201810989115A CN 109271680 A CN109271680 A CN 109271680A
- Authority
- CN
- China
- Prior art keywords
- model
- carried mass
- freedom
- freedom model
- car
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Vehicle Body Suspensions (AREA)
Abstract
本发明公开了一种将1/4汽车模型等效为二自由度模型的方法,包括以下步骤:(1)建立二自由度模型在静平衡位置的运动微分方程;(2)获取1/4汽车模型的运动参数:(3)构造二自由度模型需要辨识的参数向量;(4)将二自由度模型的运动微分方程写出矩阵形式;(5)利用矩阵方程计算方法和最小二乘法求解出二自由度模型的辨识参数。本发明将悬架阻尼系数设置为与实际更接近的拉伸阻尼系数和压缩阻尼系数,通过参数辨识方法建立的二自由度模型能够准确反映1/4汽车模型的运动,大大增强了二自由度模型的实用性和适用范围。
Description
技术领域
本发明涉及汽车悬架系统优化设计领域,特别涉及1/4汽车模型的动力学建模方法。
背景技术
随着汽车工业的发展和人们生活水平的提高,汽车产品作为代步工具进入了越来越多的普通家庭。但在中高端汽车市场中,人们对于汽车价格还不是特别敏感,但对汽车的行驶平顺性、操纵稳定性和NVH性能等汽车的各项性能的要求则越来越高。因此,各大汽车厂商和零部件公司都在着手改善汽车的各项性能,其中在通过对悬架系统的改进以同时增加汽车的行驶平顺性和操纵稳定性方面,做了很多工作。
现代汽车悬架系统主要由弹簧减振器、导向机构和橡胶弹性元件等组成,它们自身特性和连接关系决定了悬架系统的刚度和阻尼特性。目前独立悬架按结构形式可以分为:麦弗逊悬架、双横臂悬架和多连杆悬架等。麦弗逊悬架系统的结构最为简单,并且得到众多学者的深入研究,在其多体动力学模型研究方面发展得比较成熟。相比之下,双横臂悬架和多连杆悬架的结构相对复杂,在其多体动力学建模方面的研究较少。传统的二自由度模型被广泛用于悬架系统的研究,但是其模型的簧载质量和非簧载质量的运动情况能否反映实际悬架系统的运动还没有被严格证明。尤其是对于双横臂悬架和多连杆悬架,其导向机构对悬架系统的影响较大,此时,二自由度模型的准确性可能进一步下降。虽然也有学者建立了基于麦弗逊悬架或双横臂悬架的1/4汽车二维模型,但是建模和求解过程中也都进行了简化,并且模型公式或模型的求解方法很复杂,不利于控制算法的应用。而基本上所有的半主动悬架和主动悬架控制算法都是在传统的二自由度模型上推导和验证的。如果二自由度模型能够准确反映真实悬架系统的运动,将大大增强这些控制算法的实用性。
发明内容
1/4汽车模型与二自由度模型的等效原则是:利用等效模型和利用1/4汽车模型计算得到的簧载质量与非簧载质量的响应接近。首先建立1/4汽车模型,计算其簧载质量和非簧载质量的动力学响应。通过参数辨识的方法,辨识二自由度模型的参数,使得由参数采用辨识值的二自由度模型计算得到的簧载质量和非簧载质量的响应与由1/4汽车模型计算得到的动力学响应一致。
本发明目的通过如下技术方案实现:
一种将1/4汽车模型等效为二自由度模型的方法,包含以下步骤:
(1)建立二自由度模型在静平衡位置的运动微分方程;
(2)获取1/4汽车模型的运动参数:
(3)构造二自由度模型需要辨识的参数向量;
(4)将二自由度模型的运动微分方程写出矩阵形式;
(5)利用矩阵方程计算方法和最小二乘法求解出二自由度模型的辨识参数。
进一步地,步骤(1)具体包括:
簧载质量受到弹簧力fs和阻尼力fd作用,非簧载质量除了受到弹簧力fs和阻尼力fd作用外,还受到了轮胎动载荷ft的作用,根据牛顿第二定律可以写出二自由度模型的运动微分方程:
式中fs、fd和ft分别为弹簧力、阻尼力和轮胎动载荷,分别定义如下:
fs=kszsus=ks(zs-zu),
上式中,zsus为悬架相对位移、zs为簧载质量位移、为簧载质量速度、为簧载质量加速度、zu为非簧载质量位移、为非簧载质量速度、为非簧载质量加速度、zr为路面输入位移、为路面输入速度、zsus为悬架相对位移、为悬架相对速度、ms为簧载质量、mu为非簧载质量、ks为悬架刚度、cs1为拉伸阻尼系数、cs2为压缩阻尼系数。
进一步地,所述的步骤(2)具体包括:
根据实际悬架系统的硬点数据在多体动力学软件ADAMS建立1/4汽车模型,在路面位移阶跃输入,获取1/4汽车模型的簧载质量和非簧载质量的运动参数,包括:簧载质量加速度非簧载质量加速度悬架相对位移zsus、悬架相对速度和轮胎动载荷ft,构建已知参数向量:
进一步地,所述的步骤(3)具体包括:
根据二自由度模型需要辨识的参数包括簧载质量ms、非簧载质量mu、悬架刚度ks、拉伸阻尼系数cs1和压缩阻尼系数cs2,构造辨识参数向量:A=[ms,mu,ks,cs1,cs2]。
进一步地,所述的步骤(4)具体包括:
根据二自由度模型运动微分方程、已知参数向量和辨识参数向量,得到二自由度模型的运动微分方程的矩阵形式:
进一步地,所述的步骤(5)具体包括:
采用N个时间点上的等效模型的响应和1/4汽车模型的响应误差的平方和最小来识别二自由度模型的辨识参数:
令将得到:
因此,需要辨识的参数即由下式求得:
进一步地,所述步骤(2)中,采用的路面输入为位移阶跃输入,计算公式为:
单位为mm。
相比现有技术,本发明有以下积极效果:
本发明将悬架阻尼系数设置为与实际更接近的拉伸阻尼系数和压缩阻尼系数,通过参数辨识方法建立的二自由度模型能够准确反映1/4汽车模型的运动,大大增强了二自由度模型的实用性和适用范围。
附图说明
图1是1/4汽车模型;
图2是二自由度模型;
图3是路面位移阶跃输入时域图;
图4是路面位移阶跃输入频域图;
图5是1/4汽车模型在位移阶跃输入下的簧载质量加速度;
图6是1/4汽车模型在位移阶跃输入下的非簧载质量加速度;
图7是1/4汽车模型在位移阶跃输入下的悬架相对速度;
图8是1/4汽车模型在位移阶跃输入下的悬架相对位移;
图9是1/4汽车模型在位移阶跃输入下的轮胎动载荷;
图10是1/4汽车模型和二自由度模型的簧载质量加速度频域响应对比;
图11是1/4汽车模型和二自由度模型的轮跳动载荷频域响应对比。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚、明确,以下结合附图并举实施例对本发明作进一步详细描述。
1/4汽车模型与二自由度模型的等效原则是:利用等效模型和利用1/4汽车模型计算得到的簧载质量与非簧载质量的响应接近,因此,本实施提供的一种将1/4汽车模型等效为二自由度模型的方法,包括以下步骤:
(1)建立二自由度模型在静平衡位置的运动微分方程:
对图2所示的二自由度模型进行受力分析:簧载质量受到弹簧力fs和阻尼力fd作用,非簧载质量除了受到弹簧力fs和阻尼力fd作用外,还受到了轮胎动载荷ft的作用,根据牛顿第二定律可以写出二自由度模型的运动微分方程:
式中fs、fd和ft分别为弹簧力、阻尼力和轮胎动载荷,分别定义如下:
fs=kszsus=ks(zs-zu),
上式中,zsus为悬架相对位移、zs为簧载质量位移、为簧载质量速度、为簧载质量加速度、zu为非簧载质量位移、为非簧载质量速度、为非簧载质量加速度、zr为路面输入位移、为路面输入速度、zsus为悬架相对位移、为悬架相对速度、ms为簧载质量、mu为非簧载质量、ks为悬架刚度、cs1为拉伸阻尼系数、cs2为压缩阻尼系数。
(2)获取1/4汽车模型的运动参数:
根据实际悬架系统的硬点数据在多体动力学软件ADAMS建立1/4汽车模型,如图1所示,建模时采用表1所示的名义值,即单个物理元件的参数值。在图3和图4所示的路面位移阶跃输入下,获取1/4汽车模型的簧载质量和非簧载质量的运动参数,包括:簧载质量加速度非簧载质量加速度悬架相对位移zsus、悬架相对速度和轮胎动载荷ft(如图5-图9所示),构建已知参数向量:
本步骤采用的路面输入为位移阶跃输入,计算公式为:
单位为mm。
表1
(3)构造二自由度模型需要辨识的参数向量:
根据二自由度模型需要辨识的参数有簧载质量ms、非簧载质量mu、悬架刚度ks、拉伸阻尼系数cs1和压缩阻尼系数cs2,构造辨识参数向量:A=[ms,mu,ks,cs1,cs2]。
(4)将二自由度模型的运动微分方程写出矩阵形式:
根据二自由度模型运动微分方程、已知参数向量和辨识参数向量,得到二自由度模型的运动微分方程的矩阵形式:
(5)利用矩阵方程计算方法和最小二乘法求解出二自由度模型的辨识参数:
采用N个时间点上的等效模型的响应和1/4汽车模型的响应误差的平方和最小来识别二自由度模型的辨识参数:
令将得到:
因此,需要辨识的参数即由下式求得:
计算结果如表1所示,二自由度模型中的刚度和阻尼系数的数值和1/4汽车模型中相关参数的名义值有较大的差异。相对于名义值,辨识参数中的悬架弹簧刚度、拉伸阻尼系数和压缩阻尼系数的数值都减小了56%,但是二自由度模型中的质量参数的数值与名义值接近。图10和图11分别为1/4汽车模型和参数采用辨识值的二自由度模型的簧载质量加速度和轮胎动载荷的响应对比,从图中可以看出两模型的响应去向基本重合,证明的本发明方法的正确性。
虽然本发明是结合以上实施例进行描述的,但本发明并不限定于上述实施例,而只受所附权利要求的限定,本领域的普通技术人员能轻易地对其进行修改和变化,但都落入本申请所附权利要求书所限定的保护范围。
Claims (7)
1.一种将1/4汽车模型等效为二自由度模型的方法,其特征在于,包括以下步骤:
(1)建立二自由度模型在静平衡位置的运动微分方程;
(2)获取1/4汽车模型的运动参数:
(3)构造二自由度模型需要辨识的参数向量;
(4)将二自由度模型的运动微分方程写出矩阵形式;
(5)利用矩阵方程计算方法和最小二乘法求解出二自由度模型的辨识参数。
2.根据权利要求1所述的一种将1/4汽车模型等效为二自由度模型的方法,其特征在于,步骤(1)具体包括:
簧载质量受到弹簧力fs和阻尼力fd作用,非簧载质量除了受到弹簧力fs和阻尼力fd作用外,还受到了轮胎动载荷ft的作用,根据牛顿第二定律可以写出二自由度模型的运动微分方程:
式中fs、fd和ft分别为弹簧力、阻尼力和轮胎动载荷,分别定义如下:
fs=kszsus=ks(zs-zu),
上式中,zsus为悬架相对位移、zs为簧载质量位移、为簧载质量速度、为簧载质量加速度、zu为非簧载质量位移、为非簧载质量速度、为非簧载质量加速度、zr为路面输入位移、为路面输入速度、zsus为悬架相对位移、为悬架相对速度、ms为簧载质量、mu为非簧载质量、ks为悬架刚度、cs1为拉伸阻尼系数、cs2为压缩阻尼系数。
3.根据权利要求2所述的一种将1/4汽车模型等效为二自由度模型的方法,其特征在于:所述的步骤(2)具体包括:
根据实际悬架系统的硬点数据在多体动力学软件ADAMS建立1/4汽车模型,在路面位移阶跃输入,获取1/4汽车模型的簧载质量和非簧载质量的运动参数,包括:簧载质量加速度非簧载质量加速度悬架相对位移zsus、悬架相对速度和轮胎动载荷ft,构建已知参数向量:
4.根据权利要求1所述的一种将1/4汽车模型等效为二自由度模型的方法,其特征在于:所述的步骤(3)具体包括:
二自由度模型需要辨识的参数包括簧载质量ms、非簧载质量mu、悬架刚度ks、拉伸阻尼系数cs1和压缩阻尼系数cs2,构造辨识参数向量:A=[ms,mu,ks,cs1,cs2]。
5.根据权利要求2所述的一种将1/4汽车模型等效为二自由度模型的方法,其特征在于:所述的步骤(4)具体包括:
根据二自由度模型运动微分方程、已知参数向量和辨识参数向量,得到二自由度模型的运动微分方程的矩阵形式:
6.根据权利要求5所述的一种将1/4汽车模型等效为二自由度模型的方法,其特征在于:所述的步骤(5)具体包括:
采用N个时间点上的等效模型的响应和1/4汽车模型的响应误差的平方和最小来识别二自由度模型的辨识参数:
令将得到:
因此,需要辨识的参数即由下式求得:
7.根据权利要求2所述的一种将1/4汽车模型等效为二自由度模型的方法,其特征在于:所述步骤(2)中,采用的路面输入为位移阶跃输入,计算公式为:
单位为mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810989115.6A CN109271680B (zh) | 2018-08-28 | 2018-08-28 | 一种将1/4汽车模型等效为二自由度模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810989115.6A CN109271680B (zh) | 2018-08-28 | 2018-08-28 | 一种将1/4汽车模型等效为二自由度模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109271680A true CN109271680A (zh) | 2019-01-25 |
CN109271680B CN109271680B (zh) | 2023-05-23 |
Family
ID=65154499
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810989115.6A Active CN109271680B (zh) | 2018-08-28 | 2018-08-28 | 一种将1/4汽车模型等效为二自由度模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109271680B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109977520A (zh) * | 2019-03-20 | 2019-07-05 | 吉林大学 | 一种二自由度动力吸振器联合优化方法 |
CN112182764A (zh) * | 2020-09-29 | 2021-01-05 | 奇瑞汽车股份有限公司 | 车辆平顺性测试方法及装置 |
CN113887070A (zh) * | 2021-10-21 | 2022-01-04 | 浙江吉利控股集团有限公司 | 车辆平顺性分析方法、装置、设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105045102A (zh) * | 2015-06-30 | 2015-11-11 | 吉林大学 | 一种车辆侧向稳定非线性集成控制方法 |
CN105946496A (zh) * | 2016-04-27 | 2016-09-21 | 东华大学 | 一种基于天棚控制算法的悬挂控制系统 |
CN106250591A (zh) * | 2016-07-21 | 2016-12-21 | 辽宁工业大学 | 一种考虑侧倾影响的汽车行驶状态估计方法 |
CN107380161A (zh) * | 2017-07-18 | 2017-11-24 | 吉林大学 | 一种辅助驾驶员实现理想行驶轨迹的主动转向控制装置 |
CN107891723A (zh) * | 2017-11-29 | 2018-04-10 | 辽宁工业大学 | 汽车电控空气悬架的滑模控制方法及控制装置 |
-
2018
- 2018-08-28 CN CN201810989115.6A patent/CN109271680B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105045102A (zh) * | 2015-06-30 | 2015-11-11 | 吉林大学 | 一种车辆侧向稳定非线性集成控制方法 |
CN105946496A (zh) * | 2016-04-27 | 2016-09-21 | 东华大学 | 一种基于天棚控制算法的悬挂控制系统 |
CN106250591A (zh) * | 2016-07-21 | 2016-12-21 | 辽宁工业大学 | 一种考虑侧倾影响的汽车行驶状态估计方法 |
CN107380161A (zh) * | 2017-07-18 | 2017-11-24 | 吉林大学 | 一种辅助驾驶员实现理想行驶轨迹的主动转向控制装置 |
CN107891723A (zh) * | 2017-11-29 | 2018-04-10 | 辽宁工业大学 | 汽车电控空气悬架的滑模控制方法及控制装置 |
Non-Patent Citations (3)
Title |
---|
上官文斌 等: "《基于麦弗逊悬架的1 /4 汽车模型建模及其控制研究》", 《华南理工大学学报( 自然科学版)》 * |
杨鑫: "《基于M A T L A B 两自由度1 / 4 车辆悬架模型动态特性分析》", 《科技资讯》 * |
谭喜峰 等: "《二自由度汽车悬架的动态特性分析》", 《农业装备与车辆工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109977520A (zh) * | 2019-03-20 | 2019-07-05 | 吉林大学 | 一种二自由度动力吸振器联合优化方法 |
CN109977520B (zh) * | 2019-03-20 | 2020-07-17 | 吉林大学 | 一种二自由度动力吸振器联合优化方法 |
CN112182764A (zh) * | 2020-09-29 | 2021-01-05 | 奇瑞汽车股份有限公司 | 车辆平顺性测试方法及装置 |
CN113887070A (zh) * | 2021-10-21 | 2022-01-04 | 浙江吉利控股集团有限公司 | 车辆平顺性分析方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109271680B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109271680A (zh) | 一种将1/4汽车模型等效为二自由度模型的方法 | |
CN106828005B (zh) | 一种预规定瞬态性能的汽车主动悬架自适应控制方法 | |
CN100363194C (zh) | 车辆的悬架以及控制该悬架的方法 | |
CN102609578B (zh) | 一种汽车钢板弹簧多体模型的构建方法 | |
CN103076148A (zh) | 基于跌落实验双-4自由度半车模型的车辆参数识别方法 | |
CN109334377B (zh) | 磁流变半主动悬架等效替换泰勒级数lqg时滞补偿控制系统及其构造方法 | |
CN111444623A (zh) | 一种阻尼非线性商用车悬架动力学的协同优化方法及系统 | |
CN103076146B (zh) | 基于跌落试验七自由度车辆模型的车辆参数识别方法 | |
CN114417630A (zh) | 一种基于MWorks的悬架设计、仿真、分析一体化方法 | |
Mahala et al. | Mathematical models for designing vehicles for ride comfort | |
CN109145466A (zh) | 一种基于麦弗逊悬架的1/4汽车模型建模方法 | |
CN106951588B (zh) | 一种基于耦合效应的车辆isd悬架研究方法 | |
Mahajan et al. | Modeling and system identification of a quarter car suspension using Simulink | |
Zhang et al. | Modal and dynamic analysis of a vehicle with kinetic dynamic suspension system | |
Pathare et al. | Design and development of quarter car suspension test rig model and it’s simulation | |
Fallah et al. | New nonlinear model of macpherson suspension system for ride control applications | |
Kanchwala et al. | Model Building, Hardpoint Optimization & Experimental Correlation of a Single Seater EV-Toyota COMS | |
Zhao et al. | Comparative research on optimal damping matching of seat system for an off-highway dump truck | |
Gao et al. | Frequency-Domain simulation and analysis of vehicle ride comfort based on virtual proving ground | |
Balike et al. | Kineto-dynamic performance analysis of vehicle with an asymmetric suspension damper using a roll-plane model | |
Hanafi et al. | Neuro Model for Passive Suspension of a Light Car | |
Kim et al. | Ride comfort evaluation and suspension design using axiomatic design | |
CN108303066A (zh) | 一种基于悬架行程测量路面功率谱的方法 | |
Kumar et al. | PC Modeling and simulation of car suspension system | |
De Bruyne et al. | Advanced state estimator design for an active suspension |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |