CN109002600A - 一种基于逆Newmark方法的路面平整度识别方法 - Google Patents
一种基于逆Newmark方法的路面平整度识别方法 Download PDFInfo
- Publication number
- CN109002600A CN109002600A CN201810729023.4A CN201810729023A CN109002600A CN 109002600 A CN109002600 A CN 109002600A CN 201810729023 A CN201810729023 A CN 201810729023A CN 109002600 A CN109002600 A CN 109002600A
- Authority
- CN
- China
- Prior art keywords
- matrix
- surface evenness
- wheel
- inverse
- newmark
- 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
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01C—CONSTRUCTION OF, OR SURFACES FOR, ROADS, SPORTS GROUNDS, OR THE LIKE; MACHINES OR AUXILIARY TOOLS FOR CONSTRUCTION OR REPAIR
- E01C23/00—Auxiliary devices or arrangements for constructing, repairing, reconditioning, or taking-up road or like surfaces
- E01C23/01—Devices or auxiliary means for setting-out or checking the configuration of new surfacing, e.g. templates, screed or reference line supports; Applications of apparatus for measuring, indicating, or recording the surface configuration of existing surfacing, e.g. profilographs
-
- 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)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Vehicle Body Suspensions (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于逆Newmark载荷识别技术的路面平整度识别方法。本发明方法首先根据牛顿运动定律建立半车模型的运动微分方程,从而得到路面平整度的模型,然后根据Newmark逆问题分析方法进行载荷识别,最后运用复合高斯积分得到路面平整度的参数信息。本方法仅利用车辆的动态特性以及在许多传统路面平整度评估方法中容易收集的车辆振动加速度,相较于传统方法更节约时间,具有很高精度且易于实现。
Description
技术领域
本发明涉及一种路面平整度的识别方法,特别是一种基于逆Newmark载荷识别技术的路面平整度识别方法。
背景技术
路面平整度作为一种垂直输入载荷,在车辆动力学设计和车辆疲劳寿命研究中具有重要作用,并且影响车辆的安全和舒适性能。常用的评估路面平整度的方法是模拟方法和直接测量法,模拟方法主要是基于数学模型通过获得的通用谱密度函数来模拟路面的平整度,直接测量方法是通过高精度的测量仪器来识别路面平整度。
对路面平整度进行精确的评估对于车辆悬架系统动力学设计非常关键,但是获得路面平整度的精确信息往往是非常困难的,而且非常消耗时间。传统的路面平整度模拟方法较为复杂,费时较长;而测量方法需要高精度仪器,价格昂贵。
发明内容
本发明所要解决的技术问题是:
本发明运用逆Newmark方法来进行路面平整度的识别,根据半车模型的动力学参数及加速度响应,结合Newmark载荷识别方法和复合高斯积分从而识别出路面平整度,该方法能够有效节约时间且误差较小,适用于工程实际中路面平整度的识别问题。
本发明为解决上述技术问题采用以下技术方案:
一种基于逆Newmark方法的路面平整度识别方法,包括:
步骤1建立路面平整度模型:采用四自由度半车模型,得到车辆的运动微分方程如下:
其中,M是质量矩阵,K是刚度矩阵,C是阻尼矩阵,F是载荷矩阵,U是位移矩阵,是速度矩阵,是加速度矩阵,x,θ,yf,yr为半车模型的四个自由度;
mc:车辆总质量,Ic:转动惯量,mtf:前簧下质量,mtr:后簧下质量,cf:前悬架线性阻尼率,cr:后悬架线性阻尼率,kf:前悬架弹簧刚度,kr:后悬架弹簧钢度,a:前轮到质心的距离,b:后轮到质心的距离,ctf:前轮阻尼率,ctr:后轮阻尼率,ktf:前轮弹簧刚度,ktr:后轮弹簧刚度,zf:前轮的道路位移激励,zr:后轮的道路位移激励,zf的一阶导数,即前轮的道路速度激励,zr的一阶导数,即后轮的道路速度激励;
步骤2根据逆Newmark方法,求解如下t+Δt时刻的递推方程,计算载荷F:
K*U(t+Δt)=F*
其中,K*为等效刚度矩阵,U(t+Δt)是t+Δt时刻的位移矩阵,F*为总等效力,F(t+Δt)是t+Δt时刻的载荷矩阵,Δt为采样间隔,β、γ为Newmark参数,U(t)、分别是t时刻的位移、速度、加速度矩阵;
步骤3识别路面平整度,
根据识别出的载荷F,利用复合高斯积分计算路面平整度参数。
优选地,所述步骤2中计算载荷F包括:
步骤2.1根据计算初始激励F(0),其中及U(0)分别是0时刻的加速度矩阵、速度矩阵及位移矩阵,令t=0;
步骤2.2根据K*U2(t+Δt)=F2(t)且计算U2(t+Δt);
步骤2.3根据计算U(t+Δt);
步骤2.4根据U1(t+Δt)=U(t+Δt)-U2(t+Δt),计算U1(t+Δt);
步骤2.5根据K*di(t+Δt)=Fi(t+Δt),计算di(t+Δt),其中i∈{3,4},Fi(t+Δt)是假设在i自由度上作用的单位载荷,di(t+Δt)是第i自由度作用单位载荷下结构的体系位移向量;
步骤2.6根据得到λi(t+Δt),其中m是测量已知的加速度响应个数,m∈{2,3,4},(di(t+Δt))j是di(t+Δt)的第j个元素,其中j∈{1,2,3,4},u1(t+Δt)j是步骤2.4得到的U1(t+Δt)中的第j自由度的位移,λi(t+Δt)是t+Δt时刻第i自由度上的载荷幅值;
步骤2.7根据F(t+Δt)=λi(t+Δt)Fi(t+Δt),可得到t+Δt时刻的载荷矩阵F(t+Δt),令t=t+Δt,返回步骤2.2。
优选地,所述步骤3中利用复合高斯积分计算路面平整度参数包括:
根据以下方程组,计算zf和zr
优选地,对步骤2获得的载荷F进行四阶高通滤波。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
本发明提出的基于逆Newmark方法的路面平整度识别方法仅利用车辆的动态特性以及在许多传统路面平整度评估方法中容易收集的车辆振动加速度,相较于传统方法更节约时间,具有很高精度且易于实现。
附图说明
图1是车辆路面平整度模型;
图2是本发明方法的总体技术方案流程;
图3展示简谐路况下无噪声污染的识别效果;
图4展示简谐路况下受到噪声污染的识别效果;
图5展示简谐路况下经过滤波后的识别效果;
图6展示随机路况下无噪声污染的识别效果;
图7展示随机路况下受到噪声污染的识别效果;
图8展示随机路况下经过滤波后识别效果。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
Newmark方法是广泛应用于结构动力响应的数值计算方法,通过对加速度进行合理的假设,根据初始条件,通过数值积分递推得到下一时刻的响应。经过不断发展,研究出逆Newmark方法应用于载荷识别问题,该方法稳定性好,误差小。在路面平整度识别领域,尚无类似的这种方法。
本发明方法首先根据牛顿运动定律建立半车模型的运动微分方程,从而得到路面平整度的模型,然后根据Newmark逆问题分析方法进行载荷识别,最后运用复合高斯积分得到路面平整度的参数信息。具体实施步骤如下:
1.路面平整度模型的建立
如图1所示,采用四个自由度的半车模型,根据牛顿第二定律和动量定理得到车辆的运动微分方程:
其中
其中,M是质量矩阵,K是刚度矩阵,C是阻尼矩阵,F是载荷矩阵,U是位移矩阵,是速度矩阵,是加速度矩阵,x,θ,yf,yr为半车模型的四个自由度;
mc:车辆总质量,Ic:转动惯量,mtf:前簧下质量,mtr:后簧下质量,cf:前悬架线性阻尼率,cr:后悬架线性阻尼率,kf:前悬架弹簧刚度,kr:后悬架弹簧钢度,a:前轮到质心的距离,b:后轮到质心的距离,ctf:前轮阻尼率,ctr:后轮阻尼率,ktf:前轮弹簧刚度,ktr:后轮弹簧刚度,zf:前轮的道路位移激励,zr:后轮的道路位移激励,zf的一阶导数,即前轮的道路速度激励,zr的一阶导数,即后轮的道路速度激励;
2.根据逆Newmark方法识别载荷F
根据Newmark方法,t+Δt时刻的递推方程可以写成K*U(t+Δt)=F*,
其中K*为等效刚度矩阵,F*为等效力,Δt为采样间隔,β、γ为Newmark参数,U(t)、分别是位移、速度、加速度矩阵。
根据叠加原理,将U(t+Δt)分成两部分U(t+Δt)=U1(t+Δt)+U2(t+Δt)
其中:
根据Newmark方法的假设有
当所有自由度上加速度响应已知时,容易得到t+Δt时刻的载荷F(t+Δt);
当只已知部分自由度的加速度响应时,假设Fi(t+Δt)是在i自由度上作用的单位载荷,di(t+Δt)是第i自由度作用单位载荷下结构的体系位移向量,即:
F(t+Δt)=λi(t+Δt)Fi(t+Δt)
K*di(t+Δt)=Fi(t+Δt)
根据以上公式有:
根据上面方程组,可解得λi(t+Δt),于是有F(t+Δt)。
3.路面平整度的识别
在半车模型里,zf、zr是表示路面平整度的参数,根据上一步识别出的载荷F利用复合高斯积分便可识别路面平整度。
总体技术方案流程如图2所示。
根据该方法识别出路面平整度,运用四阶高通滤波方法减小误差,精度较高,并且只需要已知车辆动力学特性和测量加速度响应,节约时间,便于工程应用。
为了验证所提出的方法的有效性,引入了某种车辆模型,运用仿真技术模拟了具有简谐激励特征和随机波特征的路面平整度的识别情况。仿真结果表明:该方法能够有效应用于路面平整度的识别,效果很好。
案例一:简谐特征的路面平整度的识别
利用Matlab编程模拟路面剖面情况如正弦激励的识别效果,图3显示了加速度响应未受到噪声污染的识别结果;假设加速度响应受到20db的高斯白噪声污染时,其识别结果如图4所示,识别误差会随着时间逐渐增大,因此,噪声对识别结果影响很大,主要包括了积分误差和系统误差;为了减小加速度响应中噪声的影响,引入高通双向滤波器,如图5所示,经过滤波后误差明显减小,识别效果显著提高。
案例二:随机特征的路面平整度的识别
图6显示了测量的加速度响应未受到噪音干扰的识别效果;当加速度响应受到20db的高斯白噪声干扰时,识别误差随着时间越来越大,如图7所示;为了减小误差,使用四阶高通双向滤波器滤波,经过滤波后识别误差明显减小,如图8所示。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
本技术领域技术人员可以理解的是,本发明可以涉及用于执行本申请中所述操作中的一项或多项操作的设备。所述设备可以为所需的目的而专门设计和制造,或者也可以包括通用计算机中的已知设备,所述通用计算机有存储在其内的程序选择性地激活或重构。这样的计算机程序可以被存储在设备(例如,计算机)可读介质中或者存储在适于存储电子指令并分别耦联到总线的任何类型的介质中,所述计算机可读介质包括但不限于任何类型的盘(包括软盘、硬盘、光盘、CD-ROM、和磁光盘)、随机存储器(RAM)、只读存储器(ROM)、电可编程ROM、电可擦ROM(EPROM)、电可擦除可编程ROM(EEPROM)、闪存、磁性卡片或光线卡片。可读介质包括用于以由设备(例如,计算机)可读的形式存储或传输信息的任何机构。例如,可读介质包括随机存储器(RAM)、只读存储器(ROM)、磁盘存储介质、光学存储介质、闪存装置、以电的、光的、声的或其他的形式传播的信号(例如载波、红外信号、数字信号)等。
本技术领域技术人员可以理解的是,可以用计算机程序指令来实现这些结构图和/或框图和/或流图中的每个框以及这些结构图和/或框图和/或流图中的框的组合。可以将这些计算机程序指令提供给通用计算机、专业计算机或其他可编程数据处理方法的处理器来生成机器,从而通过计算机或其他可编程数据处理方法的处理器来执行的指令创建了用于实现结构图和/或框图和/或流图的框或多个框中指定的方法。
本技术领域技术人员可以理解的是,本发明中已经讨论过的各种操作、方法、流程中的步骤、措施、方案可以被交替、更改、组合或删除。进一步地,具有本发明中已经讨论过的各种操作、方法、流程中的其他步骤、措施、方案也可以被交替、更改、重排、分解、组合或删除。进一步地,现有技术中的具有与本发明中公开的各种操作、方法、流程中的步骤、措施、方案也可以被交替、更改、重排、分解、组合或删除。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (4)
1.一种基于逆Newmark方法的路面平整度识别方法,其特征在于,包括:
步骤1 建立路面平整度模型:采用四自由度半车模型,得到车辆的运动微分方程如下:
其中,M是质量矩阵,K是刚度矩阵,C是阻尼矩阵,F是载荷矩阵,U是位移矩阵,是速度矩阵,是加速度矩阵,x,θ,yf,yr为半车模型的四个自由度;
mc:车辆总质量,Ic:转动惯量,mtf:前簧下质量,mtr:后簧下质量,cf:前悬架线性阻尼率,cr:后悬架线性阻尼率,kf:前悬架弹簧刚度,kr:后悬架弹簧钢度,a:前轮到质心的距离,b:后轮到质心的距离,ctf:前轮阻尼率,ctr:后轮阻尼率,ktf:前轮弹簧刚度,ktr:后轮弹簧刚度,zf:前轮的道路位移激励,zr:后轮的道路位移激励,zf的一阶导数,即前轮的道路速度激励,zr的一阶导数,即后轮的道路速度激励;
步骤2 根据逆Newmark方法,求解如下t+△t时刻的递推方程,计算载荷F:
K*U(t+△t)=F*
其中,K*为等效刚度矩阵,U(t+△t)是t+△t时刻的位移矩阵,F*为总等效力,F(t+△t)是t+△t时刻的载荷矩阵,△t为采样间隔,β、γ为Newmark参数,U(t)、分别是t时刻的位移、速度、加速度矩阵;
步骤3 识别路面平整度,
根据识别出的载荷F,利用复合高斯积分计算路面平整度参数。
2.根据权利要求1所述的基于逆Newmark方法的路面平整度识别方法,其特征在于,所述步骤2中计算载荷F包括:
步骤2.1 根据计算初始激励F(0),其中及U(0)分别是0时刻的加速度矩阵、速度矩阵及位移矩阵,令t=0;
步骤2.2根据K*U2(t+△t)=F2(t)且计算U2(t+△t);
步骤2.3 根据计算U(t+△t);
步骤2.4 根据U1(t+△t)=U(t+△t)-U2(t+△t),计算U1(t+△t);
步骤2.5 根据K*di(t+△t)=Fi(t+△t),计算di(t+△t),其中i∈{3,4},Fi(t+△t)是假设在i自由度上作用的单位载荷,di(t+△t)是第i自由度作用单位载荷下结构的体系位移向量;
步骤2.6 根据得到λi(t+△t),其中m是测量已知的加速度响应个数,m∈{2,3,4},(di(t+△t))j是di(t+△t)的第j个元素,其中j∈{1,2,3,4},u1(t+△t)j是步骤2.4得到的U1(t+△t)中的第j自由度的位移,λi(t+△t)是t+△t时刻第i自由度上的载荷幅值;
步骤2.7 根据F(t+△t)=λi(t+△t)Fi(t+△t),可得到t+△t时刻的载荷矩阵F(t+△t),令t=t+△t,返回步骤2.2。
3.根据权利要求1所述的基于逆Newmark方法的路面平整度识别方法,其特征在于,所述步骤3中利用复合高斯积分计算路面平整度参数包括:
根据以下方程组,计算zf和zr
4.根据权利要求1所述的基于逆Newmark方法的路面平整度识别方法,其特征在于,还包括对步骤2获得的载荷F进行四阶高通滤波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810729023.4A CN109002600B (zh) | 2018-07-05 | 2018-07-05 | 一种基于逆Newmark方法的路面平整度识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810729023.4A CN109002600B (zh) | 2018-07-05 | 2018-07-05 | 一种基于逆Newmark方法的路面平整度识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109002600A true CN109002600A (zh) | 2018-12-14 |
CN109002600B CN109002600B (zh) | 2023-02-03 |
Family
ID=64599389
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810729023.4A Active CN109002600B (zh) | 2018-07-05 | 2018-07-05 | 一种基于逆Newmark方法的路面平整度识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109002600B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020499A (zh) * | 2019-04-19 | 2019-07-16 | 南京航空航天大学 | 一种基于逆Wilson方法的路面起伏扰动识别方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1945264A (zh) * | 2006-11-08 | 2007-04-11 | 重庆交通大学 | 由路面平整度预估沥青路面使用寿命的方法 |
-
2018
- 2018-07-05 CN CN201810729023.4A patent/CN109002600B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1945264A (zh) * | 2006-11-08 | 2007-04-11 | 重庆交通大学 | 由路面平整度预估沥青路面使用寿命的方法 |
Non-Patent Citations (1)
Title |
---|
吴庆雄等: "路面平整度PSD和IRI评价方法比较", 《交通运输工程学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020499A (zh) * | 2019-04-19 | 2019-07-16 | 南京航空航天大学 | 一种基于逆Wilson方法的路面起伏扰动识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109002600B (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Holtz-Eakin et al. | Infrastructure in a structural model of economic growth | |
Orszag et al. | Numerical computation of turbulent shear flows | |
CN107729706B (zh) | 一种非线性机械系统的动力学模型构建方法 | |
Sun | Simulation of pavement roughness and IRI based on power spectral density | |
Bogdanski et al. | Kalman and particle filtering methods for full vehicle and tyre identification | |
CN114692473B (zh) | 一种不确定车桥耦合振动系统响应的求解方法 | |
CN108009584A (zh) | 基于单源点检测的欠定盲源分离方法 | |
CN104598753A (zh) | 一种基于Brakhage v方法的桥梁移动车辆荷载识别方法 | |
CN107238500A (zh) | 汽车操纵稳定性试验快速评价系统建立方法 | |
Tota et al. | An intelligent predictive algorithm for the anti-rollover prevention of heavy vehicles for off-road applications | |
CN115130175A (zh) | 一种桥梁结构非线性随机振动和地震易损性分析方法 | |
CN106932206B (zh) | 路谱处理方法及装置 | |
CN115186377A (zh) | 一种商用车悬置系统的稳健性优化方法 | |
Medina et al. | Statistical validation of crowdsourced pavement ride quality measurements from smartphones | |
CN109002600A (zh) | 一种基于逆Newmark方法的路面平整度识别方法 | |
Zhou | Research and simulation on new active suspension control system | |
CN110096779B (zh) | 一种伺服机构动特性分析方法 | |
Monteil et al. | Robust parameter estimation of car-following models considering practical non-identifiability | |
Fraikin et al. | A fast and accurate hybrid simulation model for the large-scale testing of automated driving functions | |
Ming et al. | Application of subspace-based method in vehicle handling dynamic model identification and properties estimation | |
Kukreja et al. | Nonlinear black-box modeling of aeroelastic systems using structure detection approach: application to F/A-18 aircraft data | |
CN112784442A (zh) | 一种基于大数据的仿真分析方法及系统 | |
Chang et al. | Cause and effect chains analysis of rollover behavior with respect to chassis design | |
Anhalt et al. | Influence of the road profile accuracy on disturbance compensation in active suspension systems | |
CN110020499A (zh) | 一种基于逆Wilson方法的路面起伏扰动识别方法 |
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 |