CN108733917A - 一种基于传递矩阵的载货汽车动力学建模方法 - Google Patents
一种基于传递矩阵的载货汽车动力学建模方法 Download PDFInfo
- Publication number
- CN108733917A CN108733917A CN201810475810.0A CN201810475810A CN108733917A CN 108733917 A CN108733917 A CN 108733917A CN 201810475810 A CN201810475810 A CN 201810475810A CN 108733917 A CN108733917 A CN 108733917A
- Authority
- CN
- China
- Prior art keywords
- matrix
- cargo vehicle
- transfer matrix
- transfer
- rigid body
- 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
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Vehicle Body Suspensions (AREA)
Abstract
本发明公开了一种基于传递矩阵的载货汽车动力学建模方法,通过建立载货汽车动力学模型中力学元件的传递矩阵;对其进行合理的假设,得到具有载货汽车特点的刚柔耦合动力学模型;基于传递矩阵法理论,确定载货汽车元件编号、传递方向和切割位置;模块化建立载货汽车的传递矩阵方法,离散每个模块中的元件,进行子模块传递矩阵的建立,通过拼装得到载货汽车动力学模型的总传递矩阵,建立载货汽车动力学模型;在频域内,通过路面谱输入,根据随机振动理论,对建立模型,进行固有特性和动态特性分析。本发明方法提高了建模的效率、准确性和代码的重用率,为载货汽车平顺性分析提供理论的依据。
Description
技术领域
本发明属于汽车工程领域,具体是一种基于传递矩阵的载货汽车动力学建模方法。
背景技术
随着运输行业的发展,载货汽车成为必不可少的交通工具,人们对其要求越来越高,对平顺性也提出更高要求。平顺性是货车重要的性能指标之一,该指标反映了货车的振动激烈情况,不仅关系着运输的效率及可靠性、零部件的寿命和燃油经济性等,而且严重的会危害到人体安全。
为改善货车的平顺性,首先需要进行的是对汽车系统建模研究。然而,对于汽车这一复杂的多体系统,目前常用的牛顿-欧拉法、凯恩法和拉格朗日法等的动力学建模与仿真方法,都存在着总体动力学方程建立繁琐、易出错、矩阵存储量大以及计算效率低下等问题,成为平顺性建模与仿真的技术瓶颈。
因此,针对载货汽车多体系统,人们迫切需要一种精确高效的建模方法,为载货汽车结构设计、分析和优化提供理论依据。
发明内容
本发明的目的在克服现有技术的不足,而提供一种基于传递矩阵的载货汽车动力学建模方法,该方法解决了现有的方法中无法高效、准确的建立载货汽车动力学模型的问题,易编程,涉及矩阵存储量低和代码重用率高,能准确快速的进行载货汽车动力学模型建立,并用于动态特性分析。
实现本发明目的的技术方案是:
一种基于传递矩阵的载货汽车动力学建模方法,具体包括如下步骤:
1)建立载货汽车动力学模型中力学元件的传递矩阵,并以库的形式存储,以便调用;
2)对载货汽车进行合理的假设,得到具有载货汽车特点的刚柔耦合动力学模型,并确定载货汽车动力学模型的物理参数、力学参数和外形尺寸参数;
3)基于多体系统传递矩阵法理论,确定步骤2)中载货汽车动力学模型的元件编号、传递方向和切割位置;
4)基于模块化建立载货汽车的传递矩阵方法,把载货汽车分割成前簧下质量、后簧下质量、车架总成和驾驶室及座椅系统四大子模块,并利用步骤3)和调用步骤1)的库,进行子模块传递矩阵的建立,通过拼装得到载货汽车动力学模型的总传递矩阵,建立载货汽车动力学模型;
5)在频域内,通过路面谱输入,根据随机振动理论,对步骤4)所建立的载货汽车动力学模型,进行动态特性分析,得到驾驶室座椅底部中心的加速度功率谱密度曲线。
所述(1)中,建立的力学元件传递矩阵包括集中质量块矩阵、弹簧矩阵、多输入单输出的平面刚体矩阵、欧拉伯努利梁矩阵和虚拟刚体矩阵,其传递矩阵分别为:
所述集中质量块矩阵为:
所述弹簧矩阵为:
所述多输入单输出的平面刚体矩阵为:
式中
所述欧拉伯努利梁矩阵为:
式中
所述虚拟刚体D1包括一端为梁输入,另一端为弹簧输入;并且梁的连接点作为第一输入点,另一端输出为梁的虚拟刚体,
所述虚拟刚体D2包括两根梁输入,一个弹簧输出;
所述虚拟刚体D3包括一根梁和两个弹簧输入,一根梁输出;
步骤2)中,所述的合理假设为:
(1)货车做匀速直线运动;
(2)货车不考虑来自x,z方向的外部激励,即窜动、晃动及横摆位移较小可忽视,重点分析部件的垂向运动和绕x,z的转动;
(3)静止时,货车保持在平衡处产生微小的振幅;
(4)货车驾驶室仅有三个自由度,忽略其变形影响;
(5)货车零部件的刚度和阻尼都为线性,即线性的刚柔耦合载货汽车动力学模型;
(6)车架部分仅考虑最主要受力的4根横梁与2根纵梁,仅考虑欧拉梁的扭转和弯曲变形;
(7)驾驶者与座椅的总质量,在模型中以集中质量块表示;
(8)忽略发动机整体对货车动力学仿真的影响;
(9)车架总成和驾驶室分别关于各自中性面对称;
(10)路面的不平度冲击为平稳随机高程,仅对车轮与地面的接触点上进行输入。
所述物理参数包括质量、转动惯量,力学参数包括刚度、阻尼、抗弯刚度、抗扭刚度和材料密度;所述外形尺寸包括轮胎半径、轮距、悬架位置、各梁的位置和驾驶室悬置位置。
步骤3)中,所述的载货汽车动力学模型的元件,包括虚拟刚体元件、弹簧元件、刚体元件、弹簧阻尼器元件、梁元件以及边界部分,所述虚拟刚体元件包括u1、u2、u3、u4、u5、u6、u7、u8、u9、u10、u11和u12,所述弹簧元件包括1、 2、3、4,所述刚体元件包括5、6、25,所述弹簧阻尼器7、8、9、10、21、22、 23、24、26,所述梁元件包括11~20、19、16的梁元件被分成三部分,以19p1、 19p2、19p3、16p1、16p2、16p3进行编号,所述边界部分包括路面、顶端输出边界下标编号为0。
步骤4)中,所述的模块化建立载货汽车的传递矩阵方法,是把载货汽车分割成前簧下质量、后簧下质量、车架总成和驾驶室及座椅系统四大子模块;簧下质量,其概念为非由悬架总成中承载的质量,一般包括有前后桥、弹簧阻尼器、车轮等部件,所述的四大子模块的传递矩阵分别如下:
(1)前簧下质量的输入与输出状态矢量形式为[Y Qy]和其传递矩阵如下:
UFSZ=O
式中Z=[Z5,1 T Z5,2 T Zu9,10 T Z9,5 T]T,
(2)后簧下质量的输入与输出状态矢量形式为[Y Qy]和其传递矩阵如下:
URSZ=O
式中,Z=[Z5,3 T Z5,4 T Zu3,8 T Z7,6 T]T,
(3)车架总成中梁元件的输入与输出状态矢量形式为 [Y Θx Θz Mx Mz Qy]T及其传递矩阵如下,其中符号下标分别为传递矩阵的行和列,未标出元素都是关于0的矩阵,可根据对应状态矢量与几何矩阵判断其维数:
式中A与A1为不同维数的空间姿态坐标转换矩阵,可利用欧拉定理得到,补充方程也不较细的列出已与几何方程矩阵合并,由于虚拟剪切产生切割点的状态矢量存在关系,可把第1和7列、2和8列、4和10列通过虚拟切割理论进行合并,得到车架位置的总状态矢量如下所示及传递矩阵:
(4)驾驶室及座椅状态矢量形式为[Y Qy]和其传递矩阵如下:
(5)通过上述四个模块传递矩阵搭建,得到载货汽车多体系统总状态矢量及其维数为33×38的总传递矩阵为:
Z=[Z1,0 T Z2,0 T Z3,0 T Z4,0 T Zu1,0 T Zu3,0 T Zu4,0 T … Zu6,0 T Zu7,0 T Zu8,0 T Zu9,0 T Zu11,0 TZ27,0 T]T
式中存在8个未知的状态矢量,5个已知的边界矢量,通过把多体系统整车模型边界点的代入,即去除总传递矩阵1、3、5、7、38列可求得总传递矩阵如下式:
(6)由于其它状态矢量不为零,有如下行列式等式:
通过对上式求解,得到载货汽车多体系统的固有特性。
所述步骤5)中载货汽车动力学模型为线性时不变系统,包括单输入单输出和多输入单输出端,所述不变系统的激励与响应之间关系如下:
a、系统受到随机平稳输入{x(t)}作用下,其输出响应为{y(t)};
b、假设它们自功率谱密度分别为[Sx(ω)]、[Sy(ω)],互功率谱为[Sxy(ω)]、[Syx(ω)],而激励与响应之间的互功率谱输入输出关系为:
[Sy(ω)]=[H(ω)]*[Sx(ω)][H(ω)]T
[Sxy(ω)]=[Sx(ω)][H(ω)]T
[Syx(ω)]=[H(ω)]*[Sx(ω)]
式中,[H(ω)]表示线性系统受到单位脉冲函数产生的频响函数矩阵;[H(ω)]T表示频响函数矩阵的转置;[H(ω)]*表示频响函数矩阵的共轭,可利用输入自功率谱密度矩阵与频响函数的关系计算得到系统关于随机平稳激励的振动响应。
有益效果:本发明提供的一种基于传递矩阵的载货汽车动力学建模方法,为载货汽车平顺性分析提供理论基础,同时为产品结构设计和优化带来了方便,也起着减短产品的研发周期、减少研发成本和提高设计水平的重要作用。
附图说明
图1为本发明一种基于传递矩阵的载货汽车动力学建模方法流程图;
图2为本发明实施例提供的某型号载货汽车刚柔耦合动力学模型;
图3为本发明实施例提供的某型号载货汽车刚柔耦合动力学模型主视图;
图4为本发明实施例提供的某型号载货汽车刚柔耦合动力学模型左视图;
图5为本发明实施例提供的某型号载货汽车动力学模型元件编号、传递方向和切割位置图;
图6为本发明实施例提供的某型号载货汽车固有特性零点寻根图;
图7为本发明实施例提供的某型号载货汽车在B级路面下车速60km/h的座椅加速度功率谱密度曲线图。
具体实施方式
下面结合附图和实施例对本发明做进一步阐述,但不是对本发明的限定。
实施例:
一种基于传递矩阵的载货汽车动力学建模方法,如图1所示,具体步骤如下:
1)建立载货汽车动力学模型中力学元件的传递矩阵,并利用Matlab语言以类的形式建立每个传递矩阵库,以便调用;
建立的载货汽车动力学模型中力学元件传递矩阵包括集中质量块矩阵、弹簧矩阵、多输入单输出的平面刚体矩阵、欧拉伯努利梁矩阵和虚拟刚体矩阵,其传递矩阵分别为:
所述集中质量块矩阵为:
所述弹簧矩阵为:
所述多输入单输出平面刚体矩阵为:
式中,
所述欧拉伯努利梁矩阵为:
式中,
所述虚拟刚体D1包括一端为梁输入,另一端为弹簧输入;并且梁的连接点作为第一输入点,另一端输出为梁的虚拟刚体,
所述虚拟刚体D2包括两根梁输入,一个弹簧输出;
所述虚拟刚体D3包括一根梁和两个弹簧输入,一根梁输出;
所述物理参数包括质量、转动惯量,力学参数包括刚度、阻尼、抗弯刚度、抗扭刚度和材料密度;所述外形尺寸包括轮胎半径、轮距、悬架位置、各梁的位置和驾驶室悬置位置;
2)对载货汽车进行合理的假设,得到具有载货汽车特点的刚柔耦合动力学模型,确定物理、力学和外形尺寸参数,结合图2、图3和图4,得到载货汽车参数如表1、2所示,表中和图中符号定义为前簧下质量m5、后簧下质量m6、驾驶室m25、人椅质量m27,前簧下质量关于第一个输入点轴o1z1的转动惯量后簧下质量关于第一个输入点轴o2z2的转动惯量驾驶室关于第一个输入点轴 o3x3的转动惯量驾驶室关于第一个输入点轴o3z3的转动惯量轮胎垂直刚度kti,i=1,2,3,4,悬架垂直刚度ksi,i=7,8,9,10及阻尼csi,i=7,8,9,10,驾驶室悬置垂直刚度ksi,i=21,22,23,24及阻尼csi,i=21,22,23,24,座椅垂直刚度ks26和阻尼cs26,车架上各部分梁的抗弯刚度EIz、抗扭刚度GJP和材料密度ρ,轮胎半径R;
表1载货汽车力学参数
表2载货汽车尺寸参数(mm)
所述的合理假设为:
(1)货车做匀速直线运动;
(2)货车不考虑来自x,z方向的外部激励,即窜动、晃动及横摆位移较小可忽视,重点分析部件的垂向运动和绕x,z的转动;
(3)静止时,货车保持在平衡处产生微小的振幅;
(4)货车驾驶室仅有三个自由度,忽略其变形影响;
(5)货车零部件的刚度和阻尼都为线性,即线性的刚柔耦合载货汽车动力学模型;
(6)车架部分仅考虑最主要受力的4根横梁与2根纵梁,仅考虑欧拉梁的扭转和弯曲变形;
(7)驾驶者与座椅的总质量,在模型中以集中质量块表示;
(8)忽略发动机整体对货车动力学仿真的影响;
(9)车架总成和驾驶室分别关于各自中性面对称;
(10)路面的不平度冲击为平稳随机高程,仅对车轮与地面的接触点上进行输入。
3)基于多体系统传递矩阵法理论,确定步骤2)中载货汽车动力学模型的元件编号、传递方向和切割位置,如图5所示;所述载货汽车动力学模型的元件包括虚拟刚体元件、弹簧元件、刚体元件、弹簧阻尼器元件、梁元件以及边界部分,所述虚拟刚体元件包括u1、u2、u3、u4、u5、u6、u7、u8、u9、u10、u11和u12,所述弹簧元件包括1、2、3、4,所述刚体元件包括5、6、25,所述弹簧阻尼器7、 8、9、10、21、22、23、24、26,所述梁元件包括11~20。19、16的梁元件被分成三部分,以19p1、19p2、19p3、16p1、16p2、16p3进行编号,所述边界部分包括路面、顶端输出边界下标编号为0。
4)基于模块化建立载货汽车的传递矩阵方法,把载货汽车分割成前簧下质量、后簧下质量、车架总成和驾驶室及座椅系统四大子模块,其概念为非由悬架总成中承载的质量,一般包括有前后桥、弹簧阻尼器、车轮等部件,并利用步骤3)和调用步骤1)的库,进行子模块传递矩阵的建立,通过拼装得到载货汽车动力学模型的总传递矩阵,建立载货汽车动力学模型;
具体是模块化建立前后簧下质量传递矩阵、车架总成传递矩阵和驾驶室及座椅传递矩阵,进而搭建载货汽车多体系统总传递矩阵,在主程序中调用力学元件传递矩阵,并把编号的元件与对应的力学元件传递矩阵匹配,再将步骤2)中物理、力学和尺寸参数替代力学元件传递矩阵中定义的变量,最后通过Matlab语言编写的fMin2D.m寻根算法文件,并设置参数为寻根区域[-10,10]×[-30,30],布点个数Nx0=100,Ny0=100,精度控制dx=10-6,dy=10-6,进行固有特性的求解,如图6所示,其寻根结果如表3所示;
表3载货汽车固有频率值
所述的四大子模块的传递矩阵分别如下:
(1)前簧下质量的输入与输出状态矢量形式为[Y Qy]和其传递矩阵如下:
UFSZ=O
式中Z=[Z5,1 T Z5,2 T Zu9,10 T Z9,5 T]T
(2)后簧下质量的输入与输出状态矢量形式为[Y Qy]和其传递矩阵如下:
URSZ=O
式中,Z=[Z5,3 T Z5,4 T Zu3,8 T Z7,6 T]T
(3)车架总成中梁元件的输入与输出状态矢量形式为 [Y Θx Θz Mx Mz Qy]T及其传递矩阵如下,其中符号下标分别为传递矩阵的行和列,未标出元素都是关于0的矩阵,可根据对应状态矢量与几何矩阵判断其维数:
式中A与A1为不同维数的空间姿态坐标转换矩阵,可利用欧拉定理得到,补充方程也不较细的列出已与几何方程矩阵合并,由于虚拟剪切产生切割点的状态矢量存在关系,可把第1和7列、2和8列、4和10列通过虚拟切割理论进行合并,得到车架位置的总状态矢量如下所示及传递矩阵:
(4)驾驶室及座椅状态矢量形式为[Y Qy]和其传递矩阵如下:
(5)通过上述四个模块传递矩阵搭建,得到载货汽车多体系统总状态矢量及其维数为33×38的总传递矩阵为:
Z=[Z1,0 T Z2,0 T Z3,0 T Z4,0 T Zu1,0 T Zu3,0 T Zu4,0 T… Zu6,0 T Zu7,0 T Zu8,0 T Zu9,0 T Zu11,0 TZ27,0 T]T
式中存在8个未知的状态矢量,5个已知的边界矢量,通过把多体系统整车模型边界点的代入,即去除总传递矩阵1、3、5、7、38列可求得总传递矩阵如下式:
(6)由于其它状态矢量不为零,有如下行列式等式:
通过对上式求解,得到载货汽车多体系统的固有特性。
5)在频域内,通过路面谱输入,根据随机振动理论,对步骤4)所建立的载货汽车动力学模型,进行动态特性分析,得到驾驶室座椅底部中心的加速度功率谱密度曲线,
所述的载货汽车动力学模型为线性时不变系统,包括单输入单输出和多输入单输出,所述不变系统的激励与响应之间关系如下:
a.系统受到随机平稳输入{x(t)}作用下,其输出响应为{y(t)};
b.假设它们自功率谱密度分别为[Sx(ω)]、[Sy(ω)],互功率谱为[Sxy(ω)]、 [Syx(ω)],而激励与响应之间的互功率谱输入输出关系为:
[Sy(ω)]=[H(ω)]*[Sx(ω)][H(ω)]T
[Sxy(ω)]=[Sx(ω)][H(ω)]T
[Syx(ω)]=[H(ω)]*[Sx(ω)]
式中,[H(ω)]表示线性系统受到单位脉冲函数产生的频响函数矩阵;[H(ω)]T表示频响函数矩阵的转置;[H(ω)]*表示频响函数矩阵的共轭,可利用输入自功率谱密度矩阵与频响函数的关系计算得到系统关于随机平稳激励的振动响应,得到某型号载货汽车在B级路面下车速60km/h的座椅加速度功率谱密度曲线图,如图7所示。
Claims (6)
1.一种基于传递矩阵的载货汽车动力学建模方法,其特征在于,具体包括如下步骤:
1)建立载货汽车动力学模型中力学元件的传递矩阵,并以库的形式存储,以便调用;
2)对载货汽车进行合理的假设,得到具有载货汽车特点的刚柔耦合动力学模型,并确定载货汽车动力学模型的物理参数、力学参数和外形尺寸参数;
3)基于多体系统传递矩阵法理论,确定步骤2)中载货汽车动力学模型的元件编号、传递方向和切割位置;
4)基于模块化建立载货汽车的传递矩阵方法,把载货汽车分割成前簧下质量、后簧下质量、车架总成和驾驶室及座椅系统四大子模块,并利用步骤3)和调用步骤1)的库,进行子模块传递矩阵的建立,通过拼装得到载货汽车动力学模型的总传递矩阵,建立载货汽车动力学模型;
5)在频域内,通过路面谱输入,根据随机振动理论,对步骤4)所建立的载货汽车动力学模型,进行动态特性分析,得到驾驶室座椅底部中心的加速度功率谱密度曲线。
2.根据权利要求1所述的一种基于传递矩阵的载货汽车动力学建模方法,其特征在于,步骤1)中,建立的力学元件传递矩阵包括集中质量块矩阵、弹簧矩阵、多输入单输出的平面刚体矩阵、欧拉伯努利梁矩阵和虚拟刚体矩阵,其传递矩阵分别为:
所述集中质量块矩阵为:
所述弹簧矩阵为:
所述多输入单输出的平面刚体矩阵为:
式中
所述欧拉伯努利梁矩阵为:
式中
所述虚拟刚体D1包括一端为梁输入,另一端为弹簧输入;并且梁的连接点作为第一输入点,另一端输出为梁的虚拟刚体,
所述虚拟刚体D2包括两根梁输入,一个弹簧输出;
所述虚拟刚体D3包括一根梁和两个弹簧输入,一根梁输出;
3.根据权利要求1所述的一种基于传递矩阵的载货汽车动力学建模方法,其特征在于,步骤2)中,所述的合理假设为:
(1)货车做匀速直线运动;
(2)货车不考虑来自x,z方向的外部激励,即窜动、晃动及横摆位移较小可忽视,重点分析部件的垂向运动和绕x,z的转动;
(3)静止时,货车保持在平衡处产生微小的振幅;
(4)货车驾驶室仅有三个自由度,忽略其变形影响;
(5)货车零部件的刚度和阻尼都为线性,即线性的刚柔耦合载货汽车动力学模型;
(6)车架部分仅考虑最主要受力的4根横梁与2根纵梁,仅考虑欧拉梁的扭转和弯曲变形;
(7)驾驶者与座椅的总质量,在模型中以集中质量块表示;
(8)忽略发动机整体对货车动力学仿真的影响;
(9)车架总成和驾驶室分别关于各自中性面对称;
(10)路面的不平度冲击为平稳随机高程,仅对车轮与地面的接触点上进行输入;
所述物理参数包括质量、转动惯量,力学参数包括刚度、阻尼、抗弯刚度、抗扭刚度和材料密度;所述外形尺寸包括轮胎半径、轮距、悬架位置、各梁的位置和驾驶室悬置位置。
4.根据权利要求1所述的一种基于传递矩阵的载货汽车动力学建模方法,其特征在于,步骤3)中,所述的载货汽车动力学模型的元件,包括虚拟刚体元件、弹簧元件、刚体元件、弹簧阻尼器元件、梁元件以及边界部分,所述虚拟刚体元件包括u1、u2、u3、u4、u5、u6、u7、u8、u9、u10、u11和u12,所述弹簧元件包括1、2、3、4,所述刚体元件包括5、6、25,所述弹簧阻尼器7、8、9、10、21、22、23、24、26,所述梁元件包括11~20、19、16的梁元件被分成三部分,以19p1、19p2、19p3、16p1、16p2、16p3进行编号,所述边界部分包括路面、顶端输出边界下标编号为0。
5.根据权利要求1所述的一种基于传递矩阵的载货汽车动力学建模方法,其特征在于,步骤4)中,所述的模块化建立载货汽车的传递矩阵方法,是把载货汽车分割成前簧下质量、后簧下质量、车架总成和驾驶室及座椅系统四大子模块,所述四大子模块的传递矩阵分别如下:
(1)前簧下质量的输入与输出状态矢量形式为[Y Qy]和其传递矩阵如下:
UFSZ=O
式中Z=[Z5,1 T Z5,2 T Zu9,10 T Z9,5 T]T,
(2)后簧下质量的输入与输出状态矢量形式为[Y Qy]和其传递矩阵如下:
URSZ=O
式中,Z=[Z5,3 T Z5,4 T Zu3,8 T Z7,6 T]T,
(3)车架总成中梁元件的输入与输出状态矢量形式为[Y Θx Θz Mx Mz Qy]T及其传递矩阵如下,其中符号下标分别为传递矩阵的行和列,未标出元素都是关于0的矩阵,可根据对应状态矢量与几何矩阵判断其维数:
式中A与A1为不同维数的空间姿态坐标转换矩阵,可利用欧拉定理得到,补充方程也不较细的列出已与几何方程矩阵合并,由于虚拟剪切产生切割点的状态矢量存在关系,可把第1和7列、2和8列、4和10列通过虚拟切割理论进行合并,得到车架位置的总状态矢量如下所示及传递矩阵:
(4)驾驶室及座椅状态矢量形式为[Y Qy]和其传递矩阵如下:
(5)通过上述四个模块传递矩阵搭建,得到载货汽车多体系统总状态矢量及其维数为33×38的总传递矩阵为:
Z=[Z1,0 T Z2,0 T Z3,0 T Z4,0 T Zu1,0 T Zu3,0 T Zu4,0 T …
Zu6,0 T Zu7,0 T Zu8,0 T Zu9,0 T Zu11,0 T Z27,0 T]T
式中存在8个未知的状态矢量,5个已知的边界矢量,通过把多体系统整车模型边界点的代入,即去除总传递矩阵1、3、5、7、38列可求得总传递矩阵如下式:
(6)由于其它状态矢量不为零,有如下行列式等式:
通过对上式求解,得到载货汽车多体系统的固有特性。
6.根据权利要求1所述的一种基于传递矩阵的载货汽车动力学建模方法,其特征在于,所述步骤5)中载货汽车动力学模型为线性时不变系统,包括单输入单输出和多输入单输出端,所述不变系统的激励与响应之间关系如下:
a、系统受到随机平稳输入{x(t)}作用下,其输出响应为{y(t)};
b、假设它们自功率谱密度分别为[Sx(ω)]、[Sy(ω)],互功率谱为[Sxy(ω)]、[Syx(ω)],而激励与响应之间的互功率谱输入输出关系为:
[Sy(ω)]=[H(ω)]*[Sx(ω)][H(ω)]T
[Sxy(ω)]=[Sx(ω)][H(ω)]T
[Syx(ω)]=[H(ω)]*[Sx(ω)]
式中,[H(ω)]表示线性系统受到单位脉冲函数产生的频响函数矩阵;[H(ω)]T表示频响函数矩阵的转置;[H(ω)]*表示频响函数矩阵的共轭,可利用输入自功率谱密度矩阵与频响函数的关系计算得到系统关于随机平稳激励的振动响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810475810.0A CN108733917B (zh) | 2018-05-17 | 2018-05-17 | 一种基于传递矩阵的载货汽车动力学建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810475810.0A CN108733917B (zh) | 2018-05-17 | 2018-05-17 | 一种基于传递矩阵的载货汽车动力学建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108733917A true CN108733917A (zh) | 2018-11-02 |
CN108733917B CN108733917B (zh) | 2022-04-19 |
Family
ID=63938627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810475810.0A Active CN108733917B (zh) | 2018-05-17 | 2018-05-17 | 一种基于传递矩阵的载货汽车动力学建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108733917B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110147600A (zh) * | 2019-05-13 | 2019-08-20 | 桂林电子科技大学 | 基于典型冲击工况的载货汽车动力学模型修正方法 |
CN110489801A (zh) * | 2019-07-19 | 2019-11-22 | 国网新疆电力有限公司电力科学研究院 | 利用粒子群优化算法的发电机轴系多质块参数简化方法 |
CN110837676A (zh) * | 2019-11-05 | 2020-02-25 | 扬州大学 | 一种基于多体系统传递矩阵法的舵系统振动特性预测方法 |
CN110879872A (zh) * | 2019-11-15 | 2020-03-13 | 中国科学院长春光学精密机械与物理研究所 | 一种基于传递矩阵的车辆系统传递函数求解方法 |
CN111651849A (zh) * | 2019-02-18 | 2020-09-11 | 广州汽车集团股份有限公司 | 一种汽车实时动力学模块化建模方法、装置及存储介质 |
CN113408072A (zh) * | 2021-06-25 | 2021-09-17 | 扬州大学 | 风力机柔塔系统固有振动特性快速建模及仿真方法 |
CN116151079A (zh) * | 2023-03-09 | 2023-05-23 | 南京理工大学 | 基于多体系统传递矩阵法的声学黑洞动力吸振计算方法 |
CN116151079B (zh) * | 2023-03-09 | 2024-05-28 | 南京理工大学 | 基于多体系统传递矩阵法的声学黑洞动力吸振计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100318336A1 (en) * | 2009-06-13 | 2010-12-16 | Falangas Eric T | Method of modeling dynamic characteristics of a flight vehicle |
CN103034784A (zh) * | 2012-12-15 | 2013-04-10 | 福州大学 | 基于多体系统传递矩阵的柴油机配气系统动力学计算方法 |
CN104166772A (zh) * | 2014-08-26 | 2014-11-26 | 桂林电子科技大学 | 一种载货汽车平顺性的稳定性的仿真分析方法 |
CN104573255A (zh) * | 2015-01-22 | 2015-04-29 | 天津理工大学 | 一种基于改进多体系统传递矩阵的机械臂建模与求解方法 |
-
2018
- 2018-05-17 CN CN201810475810.0A patent/CN108733917B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100318336A1 (en) * | 2009-06-13 | 2010-12-16 | Falangas Eric T | Method of modeling dynamic characteristics of a flight vehicle |
CN103034784A (zh) * | 2012-12-15 | 2013-04-10 | 福州大学 | 基于多体系统传递矩阵的柴油机配气系统动力学计算方法 |
CN104166772A (zh) * | 2014-08-26 | 2014-11-26 | 桂林电子科技大学 | 一种载货汽车平顺性的稳定性的仿真分析方法 |
CN104573255A (zh) * | 2015-01-22 | 2015-04-29 | 天津理工大学 | 一种基于改进多体系统传递矩阵的机械臂建模与求解方法 |
Non-Patent Citations (1)
Title |
---|
任武 等: "带铰弹性开环多体臂架变位姿固有频率算法", 《中南大学学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111651849A (zh) * | 2019-02-18 | 2020-09-11 | 广州汽车集团股份有限公司 | 一种汽车实时动力学模块化建模方法、装置及存储介质 |
CN110147600A (zh) * | 2019-05-13 | 2019-08-20 | 桂林电子科技大学 | 基于典型冲击工况的载货汽车动力学模型修正方法 |
CN110147600B (zh) * | 2019-05-13 | 2022-07-05 | 桂林电子科技大学 | 基于典型冲击工况的载货汽车动力学模型修正方法 |
CN110489801A (zh) * | 2019-07-19 | 2019-11-22 | 国网新疆电力有限公司电力科学研究院 | 利用粒子群优化算法的发电机轴系多质块参数简化方法 |
CN110837676A (zh) * | 2019-11-05 | 2020-02-25 | 扬州大学 | 一种基于多体系统传递矩阵法的舵系统振动特性预测方法 |
CN110879872A (zh) * | 2019-11-15 | 2020-03-13 | 中国科学院长春光学精密机械与物理研究所 | 一种基于传递矩阵的车辆系统传递函数求解方法 |
CN113408072A (zh) * | 2021-06-25 | 2021-09-17 | 扬州大学 | 风力机柔塔系统固有振动特性快速建模及仿真方法 |
CN113408072B (zh) * | 2021-06-25 | 2023-10-13 | 扬州大学 | 风力机柔塔系统固有振动特性快速建模及仿真方法 |
CN116151079A (zh) * | 2023-03-09 | 2023-05-23 | 南京理工大学 | 基于多体系统传递矩阵法的声学黑洞动力吸振计算方法 |
CN116151079B (zh) * | 2023-03-09 | 2024-05-28 | 南京理工大学 | 基于多体系统传递矩阵法的声学黑洞动力吸振计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108733917B (zh) | 2022-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108733917A (zh) | 一种基于传递矩阵的载货汽车动力学建模方法 | |
Sharma | Ride, eigenvalue and stability analysis of three-wheel vehicle using Lagrangian dynamics | |
CN103942392A (zh) | 一种基于全生命周期的汽车底盘技术参数稳健设计方法 | |
CN102975587B (zh) | 基于双可控阻尼器的车辆半主动悬架及其控制方法 | |
CN103273976A (zh) | 一种基于行驶平顺性的坦克悬挂系统的设计方法 | |
Ozcan et al. | Optimisation of Nonlinear Spring and Damper Characteristics for Vehicle Ride and Handling Improvement | |
CN107584983A (zh) | 汽车主动悬架系统的参数化控制方法 | |
Sampò | Modelling chassis flexibility in vehicle dynamics simulation | |
yin Tang et al. | Research on suspension system based on genetic algorithm and neural network control | |
Yuezhen et al. | Study on the influence of the center of gravity of fuel cell city bus on its handling characteristics | |
CN112035953A (zh) | 基于刚柔耦合的电动三轮车车架轻量化设计方法及系统 | |
Yuan et al. | Modeling and simulation of motorcycle ride comfort based on bump road | |
Deulgaonkar et al. | Design and Analysis of Self-Balancing Two Wheeler Vehicle | |
Xu et al. | Design of FSAE Formula Racing Car Frame and Finite Element Analysis | |
Ma et al. | Simulation and Test Research of Tracked Vehicle Ride Comfort | |
Al-Oqla et al. | Integrated TRIZ-CAD Design Approach for an Ecofriendly Solar Driven Vehicle for Handicap Transportation. | |
Guosheng et al. | The research of automobile suspension system performance based on the fuzzy neural network control | |
Chu et al. | A combined dynamics model for active roll stability control of ground vehicles | |
Xu | ANSYS application in modal analysis of bus body frame | |
Lu et al. | Dynamic analysis of aircraft towing slip-out system considering vertical wheel constraints | |
Paul | Design of Chassis, Impact Attenuator, Suspension and Aerodynamic Systems of a Formula SAE Car | |
O'reilly et al. | Models Of Vehicular Collision: Development And Simulation With Emphasis On Safety I: Development of a Model for a Single Vehicle | |
CN115544661A (zh) | 一种牵引车驾驶室主动悬置动力学模型的构建方法 | |
Gadge et al. | Design, development and analysis of self-balancing electric bike | |
Fan et al. | Study on the interaction of electric vehicle roll and battery pack stability |
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 |