CN110737964A - 一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法 - Google Patents
一种基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法 Download PDFInfo
- 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
Links
- 230000008878 coupling Effects 0.000 title claims abstract description 30
- 238000010168 coupling process Methods 0.000 title claims abstract description 30
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 30
- 238000004458 analytical method Methods 0.000 title claims abstract description 20
- 230000005284 excitation Effects 0.000 title claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000001133 acceleration Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 8
- 125000000174 L-prolyl group Chemical group [H]N1C([H])([H])C([H])([H])C([H])([H])[C@@]1([H])C(*)=O 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000005489 elastic deformation Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 239000000725 suspension Substances 0.000 claims description 3
- 241000252229 Carassius auratus Species 0.000 claims description 2
- 238000012827 research and development Methods 0.000 abstract description 6
- 238000011160 research Methods 0.000 abstract description 2
- 230000004044 response Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Images
Classifications
-
- 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
- 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的绝对坐标表示为:
式中:rp——P在绝对坐标系中的位置;
R——柔性体坐标系在绝对坐标系中的位置矢量;
A——柔性体坐标系的方向矩阵;
柔性体的动能表达式为:
式中:M——柔性体的质量矩阵;
M表示为:
最后运用拉格朗日乘子法建立柔性体的运动微分方程:
式中:K——柔性体的模态刚度矩阵;
M——柔性体的质量矩阵;
D——模态阻尼矩阵;
fg——广义重力;
Q——对应于外力的广义力;
λ——对应于约束的拉格朗日乘子;
使用拉格朗日乘子法建立第n个柔性体的方程形式为:
式中:K为动能的表达式,Qn广义力,λ为拉式乘子;
对方程(1)进一步推导得到:
令
将方程(2)与***的约束方程联立,构成多体***刚柔耦合动力学方程。
本发明具有如下有益效果:本发明运用刚柔耦合原理,运用柔性车架模型替代多刚体***中的刚性车架,建立刚柔耦合整车动力学模型,该模型更接近实际物理模型;同时应用UA轮胎模型,并建立各随机振动路况的计算模型,得出全***的动态响应,提高计算的精确度,该方法可以缩短研发周期,降低研发费用,为整车平顺性研究提供了理论支持和分析方法。
附图说明:
图1为本发明基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法的流程图。
图2为matlab随机路面功率谱密度计算模型。
图3为B级随机路面激励。
图4为20km/h时的垂向加速度。
图5为40km/h时的垂向加速度。
图6为60km/h时的垂向加速度。
具体实施方式:
下面结合附图对本发明作进一步的说明。
本发明运用刚柔耦合基本理论:在研究刚柔耦合多体***动力学问题时,刚体采用多刚体***动力学建立动力学方程,由于柔性体有无限多自由度,一般采用离散化技术将柔性体离散为有限多个节点,通过节点自由度对柔性体的变形进行描述。
柔性体中任意一点P的绝对坐标可以表示为:
式中:rp——P在绝对坐标系中的位置;
R——柔性体坐标系在绝对坐标系中的位置矢量;
A——柔性体坐标系的方向矩阵;
柔性体的动能表达式为:
式中:M——柔性体的质量矩阵。
M表示为:
最后运用拉格朗日乘子法建立柔性体的运动微分方程:
式中:K——柔性体的模态刚度矩阵;
M——柔性体的质量矩阵;
D——模态阻尼矩阵;
fg——广义重力;
Q——对应于外力的广义力;
λ——对应于约束的拉格朗日乘子。
使用拉格朗日乘子法建立第n个柔性体的方程形式为:
式中:K为动能的表达式,Qn广义力,λ为拉式乘子。
对方程(1)进一步推导可以得到:
令
式中为速度二项式,其中包含哥氏力和陀螺力引起的力学效应。将上式带入(1)中,得
将方程(2)与***的约束方程联立,构成多体***刚柔耦合动力学方程。
基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法具体步骤为:
步骤一:在Pro/E软件中建立整车三维实体模型,并对其进行合理简化;
步骤二:利用Pro/E计算各主要部件的质量分布状况,得出质心位置、质量及转动惯量等质量特征,为整车动力学建模提供数据来源;
步骤三:通过ANSYS软件对车架进行柔性体建模时,由于车架多由板壳件组成,通过ANSYS软件对车架进行建模时,将单元类型选为SHELL63,通过由底向上的建模方法来建立其的有限元模型,并进行车架前十阶模态分析;
步骤四:在ANSYS中建立车架柔性体模型,并建立刚性连接区域后转为mnf格式导入ADAMS,替代多刚体***中的刚性车架,与悬架、车体等构件建立相互约束关系,建立整车刚柔耦合动力学模型;
步骤五:轮胎模型的选取,轮胎采用UA轮胎模型,该轮胎模型考虑了非稳态效果,同时包含了轮胎的纵向、侧向松弛效应,在只需要有限几个参数的情况下,有非常好的精度;
UA轮胎模型参量
步骤六:建立随机路面模型,按照GB/T 7031-2005/ISO 8608:1995《机械振动—道路路面谱测量数据报告》,平滑后的功率谱密度函数可以用最小二乘法将空间频率为0.011m-1到2.83m-1范围内的平滑数据,用一条直线进行拟合,路面功率谱密度的拟合表达式为:在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范围内的平滑数据,用一条直线进行拟合,路面功率谱密度的拟合表达式为:在MATLAB/Simulink模块中建立路面功率谱密度计算模型,得出随机路面激励,并将该随机路面激励编辑为rdf文件导入到ADAMS仿真环境中;
步骤七:在ADAMS软件中装配整车刚柔耦合动力学模型,并设置初始及边界条件,在ADAMS软件中对车辆进行刚柔耦合行驶动力学分析,得到驾驶员座椅处垂直加速度;
步骤八:计算驾驶员座椅处垂直总加权加速度的均方根值aw,用其评价人体主观舒适度,计算公式为:aw=[(1.4axw)2+(1.4ayw)2+azw 2]1/2,完成车辆行驶平顺性分析。
2.如权利要求1所述的基于随机路面激励的车辆刚柔耦合行驶平顺性分析方法,其特征在于:柔性体中任意一点P的绝对坐标表示为:
式中:rp——P在绝对坐标系中的位置;
R——柔性体坐标系在绝对坐标系中的位置矢量;
A——柔性体坐标系的方向矩阵;
柔性体的动能表达式为:
式中:M——柔性体的质量矩阵;
M表示为:
最后运用拉格朗日乘子法建立柔性体的运动微分方程:
式中:K——柔性体的模态刚度矩阵;
M——柔性体的质量矩阵;
D——模态阻尼矩阵;
fg——广义重力;
Q——对应于外力的广义力;
λ——对应于约束的拉格朗日乘子;
使用拉格朗日乘子法建立第n个柔性体的方程形式为:
式中:K为动能的表达式,Qn广义力,λ为拉式乘子;
对方程(1)进一步推导得到:
令
将方程(2)与***的约束方程联立,构成多体***刚柔耦合动力学方程。
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)
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107103131A (zh) * | 2017-04-20 | 2017-08-29 | 华中科技大学 | 一种基于Simulink平台的刚柔耦合动力学建模方法 |
-
2019
- 2019-07-21 CN CN201910658287.XA patent/CN110737964A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107103131A (zh) * | 2017-04-20 | 2017-08-29 | 华中科技大学 | 一种基于Simulink平台的刚柔耦合动力学建模方法 |
Non-Patent Citations (4)
Title |
---|
张立军;张宇;赵亮;: "基于悬架刚柔耦合模型的汽车平顺性", no. 08 * |
范国栋;吴光强;郑红;: "整车刚柔耦合建模及随机激励下的平顺性分析", no. 04 * |
陈杰平 等: "基于Matlab/Simulink的随机路面建模与不平度仿真", 农业机械学报, vol. 41, no. 3, pages 11 - 15 * |
黄泽好: "摩托车人-机刚柔耦合***动态特性研究", 中国优秀硕博士学位论文全文数据库(博士) 工程科技II辑, pages 035 - 6 * |
Cited By (12)
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 |