CN110737964A - 一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法 - Google Patents

一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法 Download PDF

Info

Publication number
CN110737964A
CN110737964A CN201910658287.XA CN201910658287A CN110737964A CN 110737964 A CN110737964 A CN 110737964A CN 201910658287 A CN201910658287 A CN 201910658287A CN 110737964 A CN110737964 A CN 110737964A
Authority
CN
China
Prior art keywords
rigid
model
flexible body
flexible
establishing
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.)
Pending
Application number
CN201910658287.XA
Other languages
English (en)
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.)
Nanjing Golden Dragon Bus Co Ltd
Original Assignee
Nanjing Golden Dragon Bus Co Ltd
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 Nanjing Golden Dragon Bus Co Ltd filed Critical Nanjing Golden Dragon Bus Co Ltd
Priority to CN201910658287.XA priority Critical patent/CN110737964A/zh
Publication of CN110737964A publication Critical patent/CN110737964A/zh
Pending legal-status Critical Current

Links

Images

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

  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

本发明公开了一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法,本发明运用刚柔耦合原理,运用柔性车架模型替代多刚体***中的刚性车架,建立刚柔耦合整车动力学模型,该模型更接近实际物理模型;同时应用UA轮胎模型,并建立各随机振动路况的计算模型,得出全***的动态响应,提高计算的精确度,该方法可以缩短研发周期,降低研发费用,为整车平顺性研究提供了理论支持和分析方法。

Description

一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法
技术领域:
本发明涉及汽车动力学分析领域,具体涉及一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法。
背景技术:
汽车行驶时,由于路面不平度引起的车辆行驶过程中的冲击和振动,会对构件材料的机械强度和疲劳强度产生影响,导致敏感的电气、电子和机械装置失灵,固定和活动结构的机械或结构损伤及旋转部件的过分磨损。如果人体受到强烈激发,并且其频率接近某一器官的共振频率,由于振动放大,常导致这个器官受损伤,损害人员的健康。
在车辆***的行驶平顺性研究中,车架的柔性变形不应被忽视,将车架作为刚性体的传统处理方法不能满足研究的需要。运用多刚体***动力学理论和柔性多体***动力学理论,采用计算机仿真技术,建立基于柔性车架的刚柔耦合整车模型,同时构建虚拟路面,引入轮胎模型,研究车辆在不同行驶速度下的行驶平顺性,有助于车辆设计参数的改进,避免样车的重复制造,有利于节约研发经费,缩短研发周期。
发明内容:
本发明提供一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法,建立一种运用多刚体***动力学理论和柔性多体***动力学理论,建立基于柔性车架的刚柔耦合整车模型,同时构建随机路面模型,引入轮胎模型,进行整车在不同路面和不同车速下的平顺性分析的方法。
本发明所采用的技术方案有:一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法,具体步骤为:
步骤一:在Pro/E软件中建立整车三维实体模型,并对其进行合理简化;
步骤二:利用Pro/E计算各主要部件的质量分布状况,得出质心位置、质量及转动惯量的质量特征,为整车动力学建模提供数据来源;
步骤三:通过ANSYS软件对车架进行柔性体建模时,由于车架多由板壳件组成,通过ANSYS软件对车架进行建模时,将单元类型选为SHELL63,通过由底向上的建模方法来建立其的有限元模型,并进行车架前十阶模态分析;
步骤四:在ANSYS中建立车架柔性体模型,并建立刚性连接区域后转为mnf格式导入ADAMS,替代多刚体***中的刚性车架,与悬架、车体构件建立相互约束关系,建立整车刚柔耦合动力学模型;
步骤五:轮胎模型的选取,轮胎采用UA轮胎模型,该轮胎模型考虑了非稳态效果,同时包含了轮胎的纵向、侧向松弛效应;
步骤六:建立随机路面模型,平滑后的功率谱密度函数用最小二乘法将空间频率为 0.011m-1到2.83m-1范围内的平滑数据,用一条直线进行拟合,路面功率谱密度的拟合表达式为:在MATLAB/Simulink模块中建立路面功率谱密度计算模型,得出随机路面激励,并将该随机路面激励编辑为rdf文件导入到ADAMS仿真环境中;
步骤七:在ADAMS软件中装配整车刚柔耦合动力学模型,并设置初始及边界条件,在ADAMS软件中对车辆进行刚柔耦合行驶动力学分析,得到驾驶员座椅处垂直加速度;
步骤八:计算驾驶员座椅处垂直总加权加速度的均方根值aw,用其评价人体主观舒适度,计算公式为:aw=[(1.4axw)2+(1.4ayw)2+azw 2]12,完成车辆行驶平顺性分析。
进一步地,柔性体中任意一点P的绝对坐标表示为:
Figure RE-GDA0002307825340000022
式中:rp——P在绝对坐标系中的位置;
R——柔性体坐标系在绝对坐标系中的位置矢量;
A——柔性体坐标系的方向矩阵;
Figure RE-GDA0002307825340000023
——P点相对于柔性体坐标系的位移矢量;
Figure RE-GDA0002307825340000024
——P点在柔性体坐标系中的方向矩阵;
Figure RE-GDA0002307825340000025
——柔性体弹性变形产生的位移矢量;
柔性体的动能表达式为:
Figure RE-GDA0002307825340000031
式中:M——柔性体的质量矩阵;
M表示为:
Figure RE-GDA0002307825340000032
最后运用拉格朗日乘子法建立柔性体的运动微分方程:
Figure RE-GDA0002307825340000033
式中:K——柔性体的模态刚度矩阵;
M——柔性体的质量矩阵;
D——模态阻尼矩阵;
fg——广义重力;
Kξ,
Figure RE-GDA0002307825340000034
——物体内部由于弹性变形和阻尼引起的广义力;
Q——对应于外力的广义力;
λ——对应于约束的拉格朗日乘子;
使用拉格朗日乘子法建立第n个柔性体的方程形式为:
Figure RE-GDA0002307825340000035
式中:K为动能的表达式,Qn广义力,λ为拉式乘子;
对方程(1)进一步推导得到:
Figure RE-GDA0002307825340000036
Figure RE-GDA0002307825340000037
Figure RE-GDA0002307825340000041
式中
Figure RE-GDA0002307825340000042
为速度二项式,其中包含哥氏力和陀螺力引起的力学效应,将上式带入(1)中,得
Figure RE-GDA0002307825340000043
将方程(2)与***的约束方程联立,构成多体***刚柔耦合动力学方程。
本发明具有如下有益效果:本发明运用刚柔耦合原理,运用柔性车架模型替代多刚体***中的刚性车架,建立刚柔耦合整车动力学模型,该模型更接近实际物理模型;同时应用UA轮胎模型,并建立各随机振动路况的计算模型,得出全***的动态响应,提高计算的精确度,该方法可以缩短研发周期,降低研发费用,为整车平顺性研究提供了理论支持和分析方法。
附图说明:
图1为本发明基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法的流程图。
图2为matlab随机路面功率谱密度计算模型。
图3为B级随机路面激励。
图4为20km/h时的垂向加速度。
图5为40km/h时的垂向加速度。
图6为60km/h时的垂向加速度。
具体实施方式:
下面结合附图对本发明作进一步的说明。
本发明运用刚柔耦合基本理论:在研究刚柔耦合多体***动力学问题时,刚体采用多刚体***动力学建立动力学方程,由于柔性体有无限多自由度,一般采用离散化技术将柔性体离散为有限多个节点,通过节点自由度对柔性体的变形进行描述。
柔性体中任意一点P的绝对坐标可以表示为:
式中:rp——P在绝对坐标系中的位置;
R——柔性体坐标系在绝对坐标系中的位置矢量;
A——柔性体坐标系的方向矩阵;
Figure RE-GDA0002307825340000051
——P点相对于柔性体坐标系的位移矢量;
Figure RE-GDA0002307825340000052
——P点在柔性体坐标系中的方向矩阵;
Figure RE-GDA0002307825340000053
——柔性体弹性变形产生的位移矢量。
柔性体的动能表达式为:
Figure RE-GDA0002307825340000054
式中:M——柔性体的质量矩阵。
M表示为:
Figure RE-GDA0002307825340000055
最后运用拉格朗日乘子法建立柔性体的运动微分方程:
Figure RE-GDA0002307825340000056
式中:K——柔性体的模态刚度矩阵;
M——柔性体的质量矩阵;
D——模态阻尼矩阵;
fg——广义重力;
Kξ,
Figure RE-GDA0002307825340000057
——物体内部由于弹性变形和阻尼引起的广义力;
Q——对应于外力的广义力;
λ——对应于约束的拉格朗日乘子。
使用拉格朗日乘子法建立第n个柔性体的方程形式为:
Figure RE-GDA0002307825340000058
式中:K为动能的表达式,Qn广义力,λ为拉式乘子。
对方程(1)进一步推导可以得到:
Figure RE-GDA0002307825340000061
Figure RE-GDA0002307825340000062
Figure RE-GDA0002307825340000063
式中为速度二项式,其中包含哥氏力和陀螺力引起的力学效应。将上式带入(1)中,得
Figure RE-GDA0002307825340000065
将方程(2)与***的约束方程联立,构成多体***刚柔耦合动力学方程。
基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法具体步骤为:
步骤一:在Pro/E软件中建立整车三维实体模型,并对其进行合理简化;
步骤二:利用Pro/E计算各主要部件的质量分布状况,得出质心位置、质量及转动惯量等质量特征,为整车动力学建模提供数据来源;
步骤三:通过ANSYS软件对车架进行柔性体建模时,由于车架多由板壳件组成,通过ANSYS软件对车架进行建模时,将单元类型选为SHELL63,通过由底向上的建模方法来建立其的有限元模型,并进行车架前十阶模态分析;
步骤四:在ANSYS中建立车架柔性体模型,并建立刚性连接区域后转为mnf格式导入ADAMS,替代多刚体***中的刚性车架,与悬架、车体等构件建立相互约束关系,建立整车刚柔耦合动力学模型;
步骤五:轮胎模型的选取,轮胎采用UA轮胎模型,该轮胎模型考虑了非稳态效果,同时包含了轮胎的纵向、侧向松弛效应,在只需要有限几个参数的情况下,有非常好的精度;
UA轮胎模型参量
Figure RE-GDA0002307825340000071
步骤六:建立随机路面模型,按照GB/T 7031-2005/ISO 8608:1995《机械振动—道路路面谱测量数据报告》,平滑后的功率谱密度函数可以用最小二乘法将空间频率为0.011m-1到2.83m-1范围内的平滑数据,用一条直线进行拟合,路面功率谱密度的拟合表达式为:
Figure RE-GDA0002307825340000072
在MATLAB/Simulink模块中建立路面功率谱密度计算模型,得出随机路面激励,并将该随机路面激励编辑为rdf文件导入到ADAMS仿真环境中。
步骤七:在ADAMS软件中装配整车刚柔耦合动力学模型,并设置初始及边界条件,分别以20km/h,40km/h,60km/h行驶,在ADAMS软件中对车辆进行刚柔耦合行驶动力学分析,得到驾驶员座椅处垂直加速度;
步骤八:计算驾驶员座椅处垂直总加权加速度的均方根值aw,用其评价人体主观舒适度,计算公式为:aw=[(1.4axw)2+(1.4ayw)2+azw 2]12。将计算出的aw值对应ISO 2631 —1《人体承受全身振动的评价》中加速度均方根值与人体主观感受对照表,完成车辆行驶平顺性分析。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下还可以作出若干改进,这些改进也应视为本发明的保护范围。

Claims (2)

1.一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法,其特征在于:具体步骤为:
步骤一:在Pro/E软件中建立整车三维实体模型,并对其进行合理简化;
步骤二:利用Pro/E计算各主要部件的质量分布状况,得出质心位置、质量及转动惯量的质量特征,为整车动力学建模提供数据来源;
步骤三:通过ANSYS软件对车架进行柔性体建模时,由于车架多由板壳件组成,通过ANSYS软件对车架进行建模时,将单元类型选为SHELL63,通过由底向上的建模方法来建立其的有限元模型,并进行车架前十阶模态分析;
步骤四:在ANSYS中建立车架柔性体模型,并建立刚性连接区域后转为mnf格式导入ADAMS,替代多刚体***中的刚性车架,与悬架、车体构件建立相互约束关系,建立整车刚柔耦合动力学模型;
步骤五:轮胎模型的选取,轮胎采用UA轮胎模型,该轮胎模型考虑了非稳态效果,同时包含了轮胎的纵向、侧向松弛效应;
步骤六:建立随机路面模型,平滑后的功率谱密度函数用最小二乘法将空间频率为0.011m-1到2.83m-1范围内的平滑数据,用一条直线进行拟合,路面功率谱密度的拟合表达式为:
Figure RE-FDA0002307825330000011
在MATLAB/Simulink模块中建立路面功率谱密度计算模型,得出随机路面激励,并将该随机路面激励编辑为rdf文件导入到ADAMS仿真环境中;
步骤七:在ADAMS软件中装配整车刚柔耦合动力学模型,并设置初始及边界条件,在ADAMS软件中对车辆进行刚柔耦合行驶动力学分析,得到驾驶员座椅处垂直加速度;
步骤八:计算驾驶员座椅处垂直总加权加速度的均方根值aw,用其评价人体主观舒适度,计算公式为:aw=[(1.4axw)2+(1.4ayw)2+azw 2]1/2,完成车辆行驶平顺性分析。
2.如权利要求1所述的基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法,其特征在于:柔性体中任意一点P的绝对坐标表示为:
Figure RE-FDA0002307825330000012
式中:rp——P在绝对坐标系中的位置;
R——柔性体坐标系在绝对坐标系中的位置矢量;
A——柔性体坐标系的方向矩阵;
Figure RE-FDA0002307825330000021
——P点相对于柔性体坐标系的位移矢量;
Figure RE-FDA0002307825330000022
——P点在柔性体坐标系中的方向矩阵;
Figure RE-FDA0002307825330000023
——柔性体弹性变形产生的位移矢量;
柔性体的动能表达式为:
式中:M——柔性体的质量矩阵;
M表示为:
最后运用拉格朗日乘子法建立柔性体的运动微分方程:
Figure RE-FDA0002307825330000026
式中:K——柔性体的模态刚度矩阵;
M——柔性体的质量矩阵;
D——模态阻尼矩阵;
fg——广义重力;
Kξ,
Figure RE-FDA0002307825330000027
——物体内部由于弹性变形和阻尼引起的广义力;
Q——对应于外力的广义力;
λ——对应于约束的拉格朗日乘子;
使用拉格朗日乘子法建立第n个柔性体的方程形式为:
Figure RE-FDA0002307825330000028
式中:K为动能的表达式,Qn广义力,λ为拉式乘子;
对方程(1)进一步推导得到:
Figure RE-FDA0002307825330000031
Figure RE-FDA0002307825330000032
Figure RE-FDA0002307825330000033
式中
Figure RE-FDA0002307825330000034
为速度二项式,其中包含哥氏力和陀螺力引起的力学效应,将上式带入(1)中,得
Figure RE-FDA0002307825330000035
将方程(2)与***的约束方程联立,构成多体***刚柔耦合动力学方程。
CN201910658287.XA 2019-07-21 2019-07-21 一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法 Pending CN110737964A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910658287.XA CN110737964A (zh) 2019-07-21 2019-07-21 一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910658287.XA CN110737964A (zh) 2019-07-21 2019-07-21 一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法

Publications (1)

Publication Number Publication Date
CN110737964A true CN110737964A (zh) 2020-01-31

Family

ID=69267399

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910658287.XA Pending CN110737964A (zh) 2019-07-21 2019-07-21 一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法

Country Status (1)

Country Link
CN (1) CN110737964A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111581716A (zh) * 2020-04-07 2020-08-25 北京科技大学 一种环道加速加载装置动力学性能的仿真方法
CN111597639A (zh) * 2020-05-19 2020-08-28 北京理工大学 一种履带与越野路面耦合分析方法
CN111767653A (zh) * 2020-07-01 2020-10-13 中国人民解放***箭军工程大学 一种基于解析胎侧刚度的重载轮胎建模方法
CN112035953A (zh) * 2020-08-24 2020-12-04 山东大学 基于刚柔耦合的电动三轮车车架轻量化设计方法及***
CN113722814A (zh) * 2021-07-22 2021-11-30 江铃汽车股份有限公司 一种基于虚拟道路测试的车辆平顺性分析方法
CN113792441A (zh) * 2021-09-28 2021-12-14 徐工集团工程机械股份有限公司道路机械分公司 一种压路机振动***数字模型的建立方法
CN117828313A (zh) * 2024-03-05 2024-04-05 中车工业研究院(青岛)有限公司 一种确定车体模态的方法、装置、列车及介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107103131A (zh) * 2017-04-20 2017-08-29 华中科技大学 一种基于Simulink平台的刚柔耦合动力学建模方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107103131A (zh) * 2017-04-20 2017-08-29 华中科技大学 一种基于Simulink平台的刚柔耦合动力学建模方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
张立军;张宇;赵亮;: "基于悬架刚柔耦合模型的汽车平顺性", no. 08 *
范国栋;吴光强;郑红;: "整车刚柔耦合建模及随机激励下的平顺性分析", no. 04 *
陈杰平 等: "基于Matlab/Simulink的随机路面建模与不平度仿真", 农业机械学报, vol. 41, no. 3, pages 11 - 15 *
黄泽好: "摩托车人-机刚柔耦合***动态特性研究", 中国优秀硕博士学位论文全文数据库(博士) 工程科技II辑, pages 035 - 6 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111581716A (zh) * 2020-04-07 2020-08-25 北京科技大学 一种环道加速加载装置动力学性能的仿真方法
CN111581716B (zh) * 2020-04-07 2022-05-06 北京科技大学 一种环道加速加载装置动力学性能的仿真方法
CN111597639A (zh) * 2020-05-19 2020-08-28 北京理工大学 一种履带与越野路面耦合分析方法
CN111597639B (zh) * 2020-05-19 2021-07-27 北京理工大学 一种履带与越野路面耦合分析方法
CN111767653A (zh) * 2020-07-01 2020-10-13 中国人民解放***箭军工程大学 一种基于解析胎侧刚度的重载轮胎建模方法
CN111767653B (zh) * 2020-07-01 2022-08-23 中国人民解放***箭军工程大学 一种基于解析胎侧刚度的重载轮胎建模方法
CN112035953A (zh) * 2020-08-24 2020-12-04 山东大学 基于刚柔耦合的电动三轮车车架轻量化设计方法及***
CN112035953B (zh) * 2020-08-24 2023-11-17 山东大学 基于刚柔耦合的电动三轮车车架轻量化设计方法及***
CN113722814A (zh) * 2021-07-22 2021-11-30 江铃汽车股份有限公司 一种基于虚拟道路测试的车辆平顺性分析方法
CN113792441A (zh) * 2021-09-28 2021-12-14 徐工集团工程机械股份有限公司道路机械分公司 一种压路机振动***数字模型的建立方法
CN117828313A (zh) * 2024-03-05 2024-04-05 中车工业研究院(青岛)有限公司 一种确定车体模态的方法、装置、列车及介质
CN117828313B (zh) * 2024-03-05 2024-05-28 中车工业研究院(青岛)有限公司 一种确定车体模态的方法、装置、列车及介质

Similar Documents

Publication Publication Date Title
CN110737964A (zh) 一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法
CN104965963B (zh) 一种刚柔耦合模型的参数化建模方法
Metaxas et al. Dynamic deformation of solid primitives with constraints
CN111581859B (zh) 一种悬架耦合非线性商用车的平顺性建模分析方法及***
CN105069826A (zh) 弹性物体变形运动的建模方法
CN111079323B (zh) 基于人-车-路耦合振动模型的动力响应预测方法及***
Docquier et al. Multiphysics modelling of multibody systems: application to car semi-active suspensions
WO2024078532A1 (zh) 复杂装备关键位置振动特性参数验证方法
Leonelli et al. A multibody motorcycle model with rigid-ring tyres: formulation and validation
Shiiba et al. Development of driving simulator with full vehicle model of multibody dynamics
Lee et al. The development of a sliding joint for very flexible multibody dynamics using absolute nodal coordinate formulation
CN117131703A (zh) 基于飞机-跑道相互作用的飞机滑跑动力响应获取方法
JP4612989B2 (ja) メカトロニクスシステムのシミュレーション方法
JP2003156413A (ja) 車体振動のシミュレーション方法
Karnopp et al. Engineering applications of dynamics
Vidoni et al. Comparison of model order reduction techniques for flexible multibody dynamics using an equivalent rigid-link system approach
Yarmohammadisatri et al. A family base optimization of a developed nonlinear vehicle suspension model using gray family design algorithm
CN111524237B (zh) 一种用于汽车模拟器中的汽车与地形交互仿真方法
Naidu et al. Ride analysis of three wheeled vehicle using MATLAB/Simulink
Urruzola et al. Integration of multibody systems in mechatronic simulation environments of any kind
Wang Interactive virtual prototyping of a mechanical system considering the environment effect. Part 1: Modeling dynamics
Schmitt Real-time simulation of flexible multibody systems in vehicle dynamics
CN105718615B (zh) 汽车底盘悬架弹簧的运动包络的生成方法及***
Pir et al. Integrated model-based approach for the optimisation of vehicle ride comfort and handling characteristics: integrated vehicle dynamics model design
Benoit et al. Real-time physics-based simulation of mechanisms and systems

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