CN112595521B - 车辆平顺性模拟测试方法 - Google Patents
车辆平顺性模拟测试方法 Download PDFInfo
- Publication number
- CN112595521B CN112595521B CN202010852883.4A CN202010852883A CN112595521B CN 112595521 B CN112595521 B CN 112595521B CN 202010852883 A CN202010852883 A CN 202010852883A CN 112595521 B CN112595521 B CN 112595521B
- Authority
- CN
- China
- Prior art keywords
- signal
- vehicle
- ride comfort
- frequency
- output signal
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M17/00—Testing of vehicles
- G01M17/007—Wheeled or endless-tracked vehicles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M17/00—Testing of vehicles
- G01M17/007—Wheeled or endless-tracked vehicles
- G01M17/04—Suspension or damping
-
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种车辆平顺性模拟测试方法。该车辆平顺性模拟测试方法包括以下步骤:(a)确定模拟目标信号Yt;(b)求解振动试验台架的响应函数;(c)迭代:得到等速扫频输入信号;(d)将车辆放在振动试验台架上,输入所述等速扫频输入信号,进行模拟振动输入,采集测量点的加速度信号。该方法为车辆平顺性的客观测试方法,输入的等速扫频信号可准确模拟真实的随机路面输入,所得测试结果精确可靠,有利于对车辆的平顺性能进行准确把控、缩短项目开发周期、降低成本和风险,特别是可实现车辆平顺性的前期虚拟分析和后期台架试验验证,以及汽车平顺性的弹性元件及阻尼元件的调校。
Description
技术领域
本发明涉及汽车性能开发与测试领域,具体而言,涉及一种车辆平顺性模拟测试方法。
背景技术
汽车行驶时,由于路面不平整引起的车辆行驶过程中的冲击和振动,会对车辆构件材料的机械强度和疲劳强度产生影响,导致电气、电子和机械装置失灵,机械结构损伤或旋转部件的过分磨损。如果人体受到强烈激发,并且其频率接近某一器官的共振频率,由于振动放大,容易导致该器官受损伤,损害人体健康。因此,在汽车开发过程中,车辆行驶过程中的平顺性研究显得非常重要,提高平顺性有利于行车安全,保护乘员的身体健康,提高汽车综合性能。
目前汽车平顺性客观评价方法,包括以下几种:吸收功率法、总体乘坐值法(BS6841-1987)、VDI2057、ISO2631和GB/T4970等,以上是世界范围内几种主要的平顺性客观评价方法,虽然存在一些争议,但是各个标准都有自己的适用范围和地域,在汽车平顺性开发中仍没有建立起广泛采用的分析和试验方法。多数汽车企业采用平顺性主观评价的方法,即根据乘员主观感觉的乘坐舒适性来进行评价,匹配悬架系统参数完成平顺性能的提升,无法很好地在设计阶段进行平顺性能把控,精确性差,导致项目开发周期的延长、经费及风险的增加。
有鉴于此,特提出本发明。
发明内容
本发明的目的在于提供一种车辆平顺性模拟测试方法,该方法为车辆平顺性的客观测试方法,输入的等速扫频信号可准确模拟真实的随机路面输入,所得测试结果精确可靠,有利于对车辆的平顺性能进行准确把控、缩短项目开发周期、降低成本和风险,特别是可实现车辆平顺性的前期虚拟分析和后期台架试验验证,以及汽车平顺性的弹性元件及阻尼元件的调校。
为了实现本发明的上述目的,特采用以下技术方案:
本发明提供了一种车辆平顺性模拟测试方法,包括以下步骤:
(a)确定模拟目标信号Yt;
(b)求解振动试验台架的响应函数:加载随机信号X到振动试验台架,测量振动试验台架的输出信号Y,通过计算Y=[H]X,得出振动试验台架传递函数H,求逆函数X=[H]-1Y;
(c)迭代:根据模拟目标信号Yt和所述逆函数,计算第一次输入信号X0=[H]-1Yt,加载X0到振动试验台架得到输出信号Y0,计算输出信号和模拟目标信号的误差并修正输入信号X1=X0+α0[H]-1(Yt-Y0),其中α0为修正系数,再依次重复加载、计算误差和修正输入信号的步骤,直到输出信号和模拟目标信号相吻合,得到等速扫频输入信号;
(d)将车辆放在振动试验台架上,输入所述等速扫频输入信号,进行模拟振动输入,采集测量点的加速度信号。
作为进一步优选的技术方案,所述振动试验台架包括四立柱振动试验台架。
作为进一步优选的技术方案,判断输出信号和模拟目标信号是否相吻合的方法包括:设定评分阈值,当输出信号的评分值大于或等于设定的评分阈值时,判断输出信号和模拟目标信号相吻合;当输出信号的评分值小于设定的评分阈值时,判断输出信号和模拟目标信号不吻合。
作为进一步优选的技术方案,所述测量点包括试验台台面、车轮轮毂、减振器车身安装点、驾驶员脚部地板或座椅导轨中的至少一种。
作为进一步优选的技术方案,在采集测量点的加速度信号后还包括数据处理的步骤。
作为进一步优选的技术方案,所述数据处理包括:采用所述等速扫频输入信号和实际输出信号的互功率谱Gxy除以所述等速扫频输入信号的自功率谱Gxx,计算得到传递特性的幅频特性曲线。
作为进一步优选的技术方案,所述数据处理包括:采用输出信号自功率谱Gyy和频率加权函数,计算加权加速度均方根值进行车辆平顺性的评价;
与现有技术相比,本发明的有益效果为:
本发明提供的车辆平顺性模拟测试方法中首先确定所需输入的模拟目标信号,然后求解振动试验台架的响应函数,再根据模拟目标信号和所述响应函数进行数据迭代,得到等速扫频输入信号,再以此信号输入到振动试验台架上进行振动试验,采集测量点的加速度信号,完成测试过程。
该方法为车辆平顺性的客观测试方法,输入的等速扫频信号可准确模拟真实的随机路面输入,所得测试结果精确可靠,有利于对车辆的平顺性能进行准确把控、缩短项目开发周期、降低成本和风险,特别是可实现车辆平顺性的前期虚拟分析和后期台架试验验证,以及汽车平顺性的弹性元件及阻尼元件的调校。
附图说明
图1为C级路面速度和速度功率谱的示意图;
图2为等速扫频输入速度和速度功率谱的示意图;
图3a为1/4车辆二自由度振动模型图;
图3b为整车七自由度振动模型图;
图4为传递特性的幅频特性曲线图;
图5为输出信号自功率谱Gyy示意图;
图6为频率加权函数图。
具体实施方式
下面将结合实施例对本发明的实施方案进行详细描述,但是本领域技术人员将会理解,下列实施例仅用于说明本发明,而不应视为限制本发明的范围。实施例中未注明具体条件者,按照常规条件或制造商建议的条件进行。
实施例1
本实施例提供了一种车辆平顺性模拟测试方法,包括以下步骤:
(a)确定模拟目标信号Yt;
(b)求解振动试验台架的响应函数:加载随机信号X到振动试验台架,测量振动试验台架的输出信号Y,通过计算Y=[H]X,得出振动试验台架传递函数H,求逆函数X=[H]-1Y;
(c)迭代:根据模拟目标信号Yt和所述逆函数,计算第一次输入信号X0=[H]-1Yt,加载X0到振动试验台架得到输出信号Y0,计算输出信号和模拟目标信号的误差并修正输入信号X1=X0+α0[H]-1(Yt-Y0),其中α0为修正系数,再依次重复加载、计算误差和修正输入信号的步骤,直到输出信号和模拟目标信号相吻合,得到等速扫频输入信号;
(d)将车辆放在振动试验台架上,输入所述等速扫频输入信号,进行模拟振动输入,采集测量点的加速度信号。
上述车辆平顺性模拟测试方法中首先确定所需输入的模拟目标信号,然后求解振动试验台架的响应函数,再根据模拟目标信号和所述响应函数进行数据迭代,得到等速扫频输入信号,再以此信号输入到振动试验台架上进行振动试验,采集测量点的加速度信号,完成测试过程。
该方法为车辆平顺性的客观测试方法,输入的等速扫频信号可准确模拟真实的随机路面输入,所得测试结果精确可靠,有利于对车辆的平顺性能进行准确把控、缩短项目开发周期、降低成本和风险,特别是可实现车辆平顺性的前期虚拟分析和后期台架试验验证,以及汽车平顺性的弹性元件及阻尼元件的调校。
在以上测试方法中,本发明之所以创造性的采用了等速扫频输入信号来模拟真实的路面输入,是基于对路面输入和扫频输入的研究而来的。发明人经研究发现,路面的速度功率谱呈常数状态(如图1所示),这与等速扫频输入形式的功率谱(如图2所示)非常接近,由此得出以等速扫频输入信号来模拟真实的路面输入的测试思路。
进一步地,原始扫频输入信号在经振动试验台架后,所输出的信号往往与原始扫频输入信号不同,如果将原始扫频输入信号直接设定为模拟目标信号,则所输出的信号并非理想的信号,无法真实模拟出路面输入情况,因而,本发明通过采用特定的数据迭代过程,使得等速扫频输入信号在振动试验台架输出后的信号能够与模拟目标信号相吻合,从而模拟出真实的出路面输入情况,保证测试的可靠性。
需要说明的是,在进行数据迭代时,可采用模拟驱动文件进行数据迭代,提高迭代效率。
应当理解的是,模拟目标信号的确定根据路面输入来合理确定,只要能够模拟真实的路面输入即可,本发明对此不做特别限制。
在一种优选的实施方式中,所述振动试验台架包括四立柱振动试验台架。四立柱振动试验台架可利用四轮信号方向差异或时间延迟,在试验过程中模拟包括单轮跳动、俯仰运动、侧倾运动、扭曲运动、过减速带等使用工况。
可选地,等速扫频原始信号可由以下matlab程序获得:
t=0:0.01:300
Vz=200*chirp(t,0.5,300,30)
Dz=cumtrapz(t,y)
其中t为数据时间信号,采样时间0.01s,数据总时间300s,Vz为垂直速度时间信号,幅值200mm/s,0时刻初始频率0.5Hz,t1(300s)时刻频率30Hz,Dz为垂直位移时间信号,单位为mm。
可选地,目标信号文件有几个方式:物理道路载荷谱采集;虚拟道路载荷谱提取;人工合成的道路驱动文件。
可选地,振动试验台架作动器位移低频复现能力为0.5Hz以上,作动器高频限值设置为50Hz。滤波器被使用过滤0.5Hz~50Hz之外的驱动和目标响应信号和传递函数。
可选地,以1/4车辆二自由度振动模型(如图3a所示)分析悬架系统相关参数对于振动特性的影响。m2为簧上质量;m1为簧下质量;K为悬架刚度;C为减振器阻尼系数;Kt为轮胎刚度;z1为车轮垂直位移;z2为车身垂直位移;q为路面不平度。车轮与车身垂直位移坐标分别为z1、z2,坐标原点选在各自的平衡位置,其运动方程为:
经推导,车身加速度对路面速度的频率特性为:
其中,A1=jωC+K,A2=-ω2m2+jωC+K,A3=-ω2m1+jωC+K+Kt,N=A3A2-A1 2。
图3b为整车七自由度振动模型图,其相关参数对振动特性的影响与1/4车辆二自由度振动模型的类似,因此不再赘述。
在一种优选的实施方式中,判断输出信号和模拟目标信号是否相吻合的方法包括:设定评分阈值,当输出信号的评分值大于或等于设定的评分阈值时,判断输出信号和模拟目标信号相吻合;当输出信号的评分值小于设定的评分阈值时,判断输出信号和模拟目标信号不吻合。本优选实施方式提供了一种具体的判断方法,需要提前预设评分阈值,然后判断输出信号的评分值与评分阈值的大小关系,进而判定输出信号与模拟目标信号是否相吻合。
在一种优选的实施方式中,所述测量点包括试验台台面、车轮轮毂、减振器车身安装点、驾驶员脚部地板或座椅导轨中的至少一种。上述测量点包括但不限于试验台台面,车轮轮毂,减振器车身安装点,驾驶员脚部地板,座椅导轨,试验台台面和车轮轮毂,减振器车身安装点和驾驶员脚部地板,驾驶员脚部地板和座椅导轨,试验台台面、车轮轮毂和减振器车身安装点,减振器车身安装点、驾驶员脚部地板和座椅导轨,或试验台台面、车轮轮毂、减振器车身安装点、驾驶员脚部地板和座椅导轨等。以上测量点的选择较为合理,可反映出真实全面的车辆振动的情况。
在一种优选的实施方式中,在采集测量点的加速度信号后还包括数据处理的步骤。对加速度信号进行数据处理后,可以直接来确定或评价车辆的平顺性好坏。
在一种优选的实施方式中,所述数据处理包括:采用所述等速扫频输入信号和实际输出信号的互功率谱Gxy除以所述等速扫频输入信号的自功率谱Gxx,计算得到传递特性的幅频特性曲线(如图4所示)。根据该幅频特性曲线,结合车辆不同系统的振动频率和传递特性,可以来确定平顺性好坏和优化方案。
在一种优选的实施方式中,所述数据处理包括:采用输出信号自功率谱Gyy(如图5所示)和频率加权函数(如图6所示),计算加权加速度均方根值进行车辆平顺性的评价;
本优选实施方式通过特定公式计算加权加速度均方根值进行车辆平顺性的评价,准确性更高。
尽管已用具体实施例来说明和描述了本发明,然而应意识到,在不背离本发明的精神和范围的情况下可以作出许多其它的更改和修改。因此,这意味着在所附权利要求中包括属于本发明范围内的所有这些变化和修改。
Claims (7)
1.一种车辆平顺性模拟测试方法,其特征在于,包括以下步骤:
(a)确定模拟目标信号Yt;
(b)求解振动试验台架的响应函数:加载随机信号X到振动试验台架,测量振动试验台架的输出信号Y,通过计算Y=[H]X,得出振动试验台架传递函数H,求逆函数X=[H]-1Y;
(c)迭代:根据模拟目标信号Yt和所述逆函数,计算第一次输入信号X0=[H]-1Yt,加载X0到振动试验台架得到输出信号Y0,计算输出信号和模拟目标信号的误差并修正输入信号X1=X0+α0[H]t1(Yt-Y0),其中α0为修正系数,再依次重复加载、计算误差和修正输入信号的步骤,直到输出信号和模拟目标信号相吻合,得到等速扫频输入信号;
(d)将车辆放在振动试验台架上,输入所述等速扫频输入信号,进行模拟振动输入,采集测量点的加速度信号;
等速扫频原始信号可由以下matlab程序获得:
t=0:0.01:300
Vz=200*chirp(t,0.5,300,30)
Dz=cumtrapz(t,y)
其中t为数据时间信号,采样时间0.01s,数据总时间300s,Vz为垂直速度时间信号,幅值200mm/s,0时刻初始频率0.5Hz,t1(300s)时刻频率30Hz,Dz为垂直位移时间信号,单位为mm。
2.根据权利要求1所述的车辆平顺性模拟测试方法,其特征在于,所述振动试验台架包括四立柱振动试验台架。
3.根据权利要求1所述的车辆平顺性模拟测试方法,其特征在于,判断输出信号和模拟目标信号是否相吻合的方法包括:设定评分阈值,当输出信号的评分值大于或等于设定的评分阈值时,判断输出信号和模拟目标信号相吻合;当输出信号的评分值小于设定的评分阈值时,判断输出信号和模拟目标信号不吻合。
4.根据权利要求1所述的车辆平顺性模拟测试方法,其特征在于,所述测量点包括试验台台面、车轮轮毂、减振器车身安装点、驾驶员脚部地板或座椅导轨中的至少一种。
5.根据权利要求1-4任一项所述的车辆平顺性模拟测试方法,其特征在于,在采集测量点的加速度信号后还包括数据处理的步骤。
6.根据权利要求5所述的车辆平顺性模拟测试方法,其特征在于,所述数据处理包括:采用所述等速扫频输入信号和实际输出信号的互功率谱Gxy除以所述等速扫频输入信号的自功率谱Gxx,计算得到传递特性的幅频特性曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010852883.4A CN112595521B (zh) | 2020-08-24 | 2020-08-24 | 车辆平顺性模拟测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010852883.4A CN112595521B (zh) | 2020-08-24 | 2020-08-24 | 车辆平顺性模拟测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112595521A CN112595521A (zh) | 2021-04-02 |
CN112595521B true CN112595521B (zh) | 2022-12-09 |
Family
ID=75180122
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010852883.4A Active CN112595521B (zh) | 2020-08-24 | 2020-08-24 | 车辆平顺性模拟测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112595521B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113532879A (zh) * | 2021-07-15 | 2021-10-22 | 中国第一汽车股份有限公司 | 数据采集方法、装置、电子设备及存储介质 |
CN114035719B (zh) * | 2021-10-09 | 2024-04-26 | 湖南麒麟信安科技股份有限公司 | 一种远程桌面流畅度性能评价方法、系统及介质 |
CN113887070A (zh) * | 2021-10-21 | 2022-01-04 | 浙江吉利控股集团有限公司 | 车辆平顺性分析方法、装置、设备及存储介质 |
CN114646475B (zh) * | 2022-03-16 | 2024-05-03 | 中国第一汽车股份有限公司 | 四柱平顺性试验数据处理方法、装置和电子设备 |
CN115452290A (zh) * | 2022-08-31 | 2022-12-09 | 中国第一汽车股份有限公司 | 车辆的数据处理方法和目标车辆 |
CN116256134B (zh) * | 2023-01-09 | 2024-05-28 | 河南科技大学 | 一种车辆振动测试装置系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104120644A (zh) * | 2013-04-26 | 2014-10-29 | 同济大学 | 一种基于重力加速度传感器的路面平整度检测方法 |
CN104316328A (zh) * | 2014-05-29 | 2015-01-28 | 潍柴动力股份有限公司 | 一种车辆平顺性的测定方法 |
CN104344962A (zh) * | 2013-08-02 | 2015-02-11 | 重庆车辆检测研究院有限公司 | 摩托车道路模拟试验方法 |
CN109933887A (zh) * | 2019-03-11 | 2019-06-25 | 桂林电子科技大学 | 一种基于牵引车平顺性的悬架非线性特性优化方法 |
WO2020118338A1 (de) * | 2018-12-10 | 2020-06-18 | Avl List Gmbh | Prüfstand mit einem prüflingsmodell |
-
2020
- 2020-08-24 CN CN202010852883.4A patent/CN112595521B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104120644A (zh) * | 2013-04-26 | 2014-10-29 | 同济大学 | 一种基于重力加速度传感器的路面平整度检测方法 |
CN104344962A (zh) * | 2013-08-02 | 2015-02-11 | 重庆车辆检测研究院有限公司 | 摩托车道路模拟试验方法 |
CN104316328A (zh) * | 2014-05-29 | 2015-01-28 | 潍柴动力股份有限公司 | 一种车辆平顺性的测定方法 |
WO2020118338A1 (de) * | 2018-12-10 | 2020-06-18 | Avl List Gmbh | Prüfstand mit einem prüflingsmodell |
CN109933887A (zh) * | 2019-03-11 | 2019-06-25 | 桂林电子科技大学 | 一种基于牵引车平顺性的悬架非线性特性优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112595521A (zh) | 2021-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112595521B (zh) | 车辆平顺性模拟测试方法 | |
US5610330A (en) | Effective road profile control method for a spindle-coupled road simulator | |
CN104573343B (zh) | 一种沥青路面平整度舒适性现场评价方法及仿真评价方法 | |
CN104175920B (zh) | 车辆座椅悬置磁流变减振器最优控制电流的设计方法 | |
Dodds et al. | Laboratory road simulation for full vehicle testing: a review | |
CN102944422B (zh) | 汽车传动系统行驶载荷疲劳再现试验方法 | |
Konieczny | Analysis of simplifications applied in vibration damping modelling for a passive car shock absorber | |
CN113570057B (zh) | 一种基于模型训练的车辆轮心垂向位移测量方法和装置 | |
CN113408204A (zh) | 一种路基服役性能的健康诊断装置与方法 | |
CN111896275A (zh) | 一种整车道路模拟试验的测试系统及方法 | |
CN101644629A (zh) | 测量小动态振幅下转向柱的扭转动力学特性的方法 | |
Lépine et al. | A laboratory excitation technique to test road bike vibration transmission | |
CN113434954B (zh) | 一种振动式路面平整度测试车的标定方法 | |
US5750890A (en) | Method and apparatus for modelling a tire for use with a vehicle spindle-coupled simulator | |
Soria et al. | Operational modal analysis and the performance assessment of vehicle suspension systems | |
CN115270296A (zh) | 一种商用车驾驶室疲劳耐久性分析方法及系统 | |
Doumiati et al. | Adaptive control scheme for road profile estimation: application to vehicle dynamics | |
CN108731908B (zh) | 用于模拟行车过程的外界激励加载试验装置 | |
CN110399683A (zh) | 基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 | |
CN104330269B (zh) | 一种用于测试轮胎轮辋组合体阻尼比的方法 | |
CN104455157A (zh) | 车辆座椅悬置液压减振器非线性速度特性参数的获取方法 | |
CN115358130B (zh) | 一种车辆段振动噪声仿真荷载实现方法 | |
KR100941968B1 (ko) | 인공윤하중을 이용한 교량의 동적거동 해석방법 | |
CN116863905A (zh) | 针对汽车路噪特征峰值的主动降噪方法及车辆 | |
CN115524086A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |