基于功率谱密度的路面不平整度数值模拟

合集下载

路面不平度构造以及功率谱密度生成检测

路面不平度构造以及功率谱密度生成检测

路面不平度构造以及功率谱密度生成检测一构造C级路面不平度构造,采用谐波叠加法(又称为三角级数法),常见空间频率范围为:0.011<n<2.83,车速为30km/h, 采用Parseval 公式和相关定理,谐波叠加,生成C级路面不平度时域信号。

程序如下:n1=0.011;nh=2.83;n0=0.1;G0=256*10^(-6);N=2000;v=30/3.6;L=1000; %区间数车速距离fl=n1*v;fh=nh*v; %时间频域fk=fl+(fh-fl)/(2*N):(fh-fl)/N:fh-(fh-fl)/(2*N); %中心频率faik=2*pi*rand(1,N); %相位角ak=zeros(1,N); %系数for i=1:Nak(i)=sqrt(2*G0*n0^2*v*(fh-fl)/N/(fl+(i-1)*(fh-fl)/N)/(fl+i*(fh-fl)/N ));endt=0:1/fh/2:L/v; %时间采样频率为最高时间频率的两倍Lt=length(t);qk=zeros(N,Lt);for i=1:Nqk(i,:)=ak(i)*sin(2*pi*fk(i)*t+faik(i));endqt=zeros(Lt,1);qt(:,1)=sum(qk)';plot(t,qt(:,1)),grid ontitle('c级路面时域不平度q(t)')xlabel('时间t/s')ylabel('位移q/m')GC=sqrt(sum(qt.^2)/length(qt)) %路面不平度均方值L=0:0.05/2.83:1000;Ll=length(L);qll=zeros(N,Ll);nk=0.011+(2.83-0.011)/(2*N):(2.83-0.011)/N:2.83-(2.83-0.011)/(2*N);for i=1:Nqll(i,:)=ak(i)*sin(2*pi*nk(i)*L+faik(i)); endql=zeros(Ll,1);ql(:,1)=sum(qll)';figure(2);plot(L,ql(:,1)),grid onGC =0.0143路面功率谱密度:figure(3)ns=10*nh;% 空间采样频率[pxx,m]=psd(qt,512,ns);% 512为快速傅里叶变换点数loglog(m,pxx)grid on ;hold on ;n=0.011:0.001:2.83;y1=128e-6./n.^2 ;y2=512e-6./n.^2 ;loglog(n,y1) ;grid on ;loglog(n,y2) ;grid on;xlabel('空间频率/(1/m)')ylabel('路面公路车谱密度G(n)/[m^2/(1/m)]')二构造F级路面不平度,常见空间频率范围为:0.011<n<2.83,车速为100km/h。

路面不平度的模拟与汽车非线性随机振动的研究报告

路面不平度的模拟与汽车非线性随机振动的研究报告

清华大学学报自然科学版JOURNALOF TSINGHUAUNIVERSITY SCIENCEANDTECHNOLOGY1999年第39卷第8期Vol.39 No.81999路面不平度的模拟与汽车非线性随机振动的研究*金睿臣,宋健文摘预测汽车的随机振动响应对汽车的开发设计是非常重要的。

实际汽车存在许多非线性环节,需采用非线性振动模型进行研究,在这种情况下,通常采用的频域分析方法一般不再适用。

应用机械系统分析软件ADAMS建立了11自由度汽车非线性振动模型,并用由伪白噪声法生成的符合实际路面统计特性的伪随机序列来模拟路面不平度。

在此基础上,利用数值算法在时域中对汽车的非线性随机振动响应进行了计算机仿真计算研究。

结果表明,这种方法对研究汽车的非线性随机振动是有效的。

关键词汽车动力学;ADAMS软件;非线性随机振动;路面不平度分类号U 461;O 322Simulation of the road irregularity and study of nonlinear randomvibration of the automobileJIN Ruichen,SONG JianDepartment of Automotive Engineering,State Key Laboratory of Automotive Safety and Energy Conservation,Tsinghua University,Beijing 100084,China Abstract To use the simulation technique is very important to predict the random vibration of the automobile.Because there are many nonlinear factors in a real automobile,a nonlinear vibration model should be necessarily used.In this case,the frequency domain methods can not be applicable.Under the help of the mechanical system simulation program ADAMS,an 11 DOF nonlinear vibration model of the automobile was built.By means of pseudo white noise,pseudo random sequences,which can simulate the random irregularities of a road,were generated.Based on these,using numerical method,the random vibration of the automobile was studied.The results of simulation have demonstrated the validity of the method.Key words vehicle dynamics;program ADAMS;nonlinear random vibration; road irregularities汽车以一定的速度行驶时,路面的随机不平度通过轮胎、悬架等弹性、阻尼元件传递到车身上,并通过座椅将振动传递到人体。

基于ADAMS的随机路面不平度建模及参数选择

基于ADAMS的随机路面不平度建模及参数选择

2.Hefei University of Technology ,Hefei 230009 ,China ) Abstract :The vertical motions and pitch of vehicle has an effect on the comfort and safety of passengers.The vertical vibration is depends on the road unevenness profiles besides engine excitation.It shows that the road unevenness imitated by Partial Wave Adding Model and compared with national standards.Analysed numbers of harmonic have affection for RMS value of road unevenness.demonstrate that the PSD yields of the random time signal are consistent with the road grades standard.It offers reliable excitation signals for controlling research of vehicle vertical dynamic. Key words :random vibration ;road unevenness ;harmonic superposition ;Parameters Selection
基于 ADAMS 的随机路面不平度建模及参数选择 / 王俊龙 , 汪
洋 , 王吉华
设 计·研 究

基于IFFT法的路面不平度时域模拟方法

基于IFFT法的路面不平度时域模拟方法

基于IFFT法的路面不平度时域模拟方法作者:鲍家定伍建伟王瀚超莫秋云来源:《现代电子技术》2016年第20期摘要:路面不平度时域模拟方法存在模拟精度低下等问题。

考虑到IFFT法能够精确而又简单地重构道路的时域模型,是一种普适性的方法,其具有计算量小、计算简单高效等的特点,可为后续车辆动力学仿真分析提供实时的时域模型。

基于IFFT法(逆快速傅里叶变换法)进行路面不平度的时域模拟,总结了IFFT法时域模拟的五个步骤,并利用Matlab语言开发了基于该方法的时域模拟函数。

通过与标准的功率谱密度曲线比较,其吻合程度高,表明该方法是正确可靠的,同时也说明了模拟参数选取的合理性。

关键词:路面不平度; IFFT;时域模拟; Matlab;功率谱密度中图分类号: TN911.71⁃34; U461.4 文献标识码: A 文章编号: 1004⁃373X(2016)20⁃0008⁃04Abstract: Currently, there are some problems in the time⁃domain simulation of road roughness, such as low simulation precision. Therefore, IFFT (inverse fast Fourier transform)is used to do time⁃domain simulation of road roughness because it can reconstruct time⁃domain model of road accurately and simply. It has the characteristics of small computational amount,simple calculation, etc. The time⁃domain simulation of road roughness was performed on the basis of IFFT method. Five steps of the time⁃domain simulation of IFFT method were summarized. The time⁃domain simulation function based on the method was developed with Matlab language. In comparison with the standard power spectral density curve, its match degree is high, which shows that the method is right and reliable, and also illustrates the rationality of the simulation parameter.Keywords: road roughness; IFFT; time⁃domain simulation; Matlab; power spectral density随着人们对汽车平顺性的要求不断提高,研究人员逐渐认识到汽车振动系统非线性动力学分析的重要性,纷纷通过建立非线性动力学模型来提高平顺性仿真的精度和优化的质量。

路面功率谱密度换算及不平度建模理论研究

路面功率谱密度换算及不平度建模理论研究
长春 , 1 3 0 0 2 2 )
重庆 , 4 0 0 0 3 0 )
( 3 . 吉林 大 学 汽 车 仿 真 与 控 制 国家 重 点 实 验 室 摘要
为 了 明晰 路 面不 平度 空 间 域 统 计 量 的 计 算 , 以及几种重要功 率谱密度 ( p o we r s p e c t r a l d e n s i t y , 简称 P S D ) 之
的傅 里 叶逆 变 换 法 是 路 面 空 间建 模 的常 用 方 法 之


其 建模 思想 比较 成熟 , 但 尚有一 个关键 步 骤缺乏
严格 理论 分析 ] —— 功率 谱密 度离 散化 。现 有文 献 通 常都 略过这 些 问题| 9 ] 。
另外 , 在应 用傅 里 叶逆变 换法 做路 面模 拟时 , 细
频域均离散从而可以由7yi项级数和与离散傅里叶变换为空间频率分辨率从以上分析可看到功率谱密度离散化本质是将截断路面信号的幅值谱密度离散化得到延拓周期信号的傅里叶级数以此来合成路面信号离散化结果验证功率谱密度离散化的目的是得到信号的幅值谱那么发来给以验证将长度和周期均为的路面不平度连续信号的傅里叶级数近似表达成项级数和的形式可得到连续周期信号傅里叶级数与7yi的关系只存在一个系数差别这个系数将在逆变换时得以消除故可以代表信号幅值谱亦可如此分析假设信号点根据频域卷积定理信号频谱被延拓成个路面客观存在即傅里叶级数不变于是得到即为信号幅值谱采样点数频谱点对称设置根据上述分析离散功率谱密度估计式为为空间频率分辨率即采样路面总长度的倒数相位角在上均匀分布在210中可以通过傅里叶逆变换生成路面空间域随机信号离散傅里叶变换的第一个点直流分量不存在对称对于模拟路面数据如果是偶数则关于点对称若为奇数则关于8l的共轭复数将两段频谱数据组合后进行离散傅里叶逆变换即可得到路面信号功率谱密度验证正确的模拟路面除了要满足均值和标准差的条件外即均值接近于零标准差与标准值接近并不随车速

路面不平度构造以及功率谱密度生成检测

路面不平度构造以及功率谱密度生成检测

路面不平度构造以及功率谱密度生成检测一构造C级路面不平度构造,采用谐波叠加法(又称为三角级数法),常见空间频率范围为:0.011<n<2.83,车速为30km/h, 采用Parseval 公式和相关定理,谐波叠加,生成C级路面不平度时域信号。

程序如下:n1=0.011;nh=2.83;n0=0.1;G0=256*10^(-6);N=2000;v=30/3.6;L=1000; %区间数车速距离fl=n1*v;fh=nh*v; %时间频域fk=fl+(fh-fl)/(2*N):(fh-fl)/N:fh-(fh-fl)/(2*N); %中心频率faik=2*pi*rand(1,N); %相位角ak=zeros(1,N); %系数for i=1:Nak(i)=sqrt(2*G0*n0^2*v*(fh-fl)/N/(fl+(i-1)*(fh-fl)/N)/(fl+i*(fh-fl)/N ));endt=0:1/fh/2:L/v; %时间采样频率为最高时间频率的两倍Lt=length(t);qk=zeros(N,Lt);for i=1:Nqk(i,:)=ak(i)*sin(2*pi*fk(i)*t+faik(i));endqt=zeros(Lt,1);qt(:,1)=sum(qk)';plot(t,qt(:,1)),grid ontitle('c级路面时域不平度q(t)')xlabel('时间t/s')ylabel('位移q/m')GC=sqrt(sum(qt.^2)/length(qt)) %路面不平度均方值L=0:0.05/2.83:1000;Ll=length(L);qll=zeros(N,Ll);nk=0.011+(2.83-0.011)/(2*N):(2.83-0.011)/N:2.83-(2.83-0.011)/(2*N);for i=1:Nqll(i,:)=ak(i)*sin(2*pi*nk(i)*L+faik(i)); endql=zeros(Ll,1);ql(:,1)=sum(qll)';figure(2);plot(L,ql(:,1)),grid onGC =0.0143路面功率谱密度:figure(3)ns=10*nh;% 空间采样频率[pxx,m]=psd(qt,512,ns);% 512为快速傅里叶变换点数loglog(m,pxx)grid on ;hold on ;n=0.011:0.001:2.83;y1=128e-6./n.^2 ;y2=512e-6./n.^2 ;loglog(n,y1) ;grid on ;loglog(n,y2) ;grid on;xlabel('空间频率/(1/m)')ylabel('路面公路车谱密度G(n)/[m^2/(1/m)]')二构造F级路面不平度,常见空间频率范围为:0.011<n<2.83,车速为100km/h。

路面不平度的数值模拟方法研究

路面不平度的数值模拟方法研究

目录目录............................................... 错误!未定义书签。

摘要.. (3)A BSTRACT (4)1 前言 (5)1.1课题研究的意义 (5)1.2相关领域发展历史和发展趋势 (5)1.3汽车振动的建模方法 (6)1.4主要研究内容 (7)2 悬架工作过程分析 (9)2.1悬架的功用 (9)2.2悬架的组成 (9)2.3汽车悬架的分类 (11)2.4汽车悬架工作过程 (13)3 MATLAB/SIMULINK简介 (14)4路面激励的时域模型 (16)4.1路面不平度的功率谱密度 (16)4.2路面激励的时域模型 (17)5汽车悬架的五自由度模型 (20)5.1汽车模型简化的假设 (20)5.2五自由度动力学模型 (20)5.3系统动力学方程 (22)5.4路面激励时域模型 (23)5.5悬架特性五自由度仿真模型 (25)5.5某汽车振动仿真结果及分析 (27)6 悬架特性分析与评价 (32)6.1汽车行驶速度对平顺性的影响分析 (32)6.2路面不平度对平顺性的影响分析 (33)6.3悬架弹性元件刚度对平顺性的影响分析 (35)6.4悬架减振器阻尼对平顺性的影响分析 (40)7 结论 (43)总结与体会 (47)谢辞 (48)参考文献 (49)路面不平度的数值模拟方法研究摘要汽车平顺性是汽车的主要性能之一。

它的优劣直接关系到驾驶员的驾驶舒适性和乘员的舒适性。

如何提高汽车平顺性已成为汽车设计人员越来越关心的问题。

随着随机振动理论的发展,计算机仿真也大量运用到汽车性能分析当中。

本文基于汽车振动动力学分析方法,建立五自由度汽车振动模型。

该五自由度模型的五个自由度包括纵向两车轮的两个自由度、车身的两个自由度和驾驶员一个自由度。

汽车振动模型的输入有路面不平度、车速。

因此,本次采用白噪声法以白噪声作为输入建立了A、B、C、D各级路面模型。

通过Matlab/ simulink 软件建立仿真模型,以某汽车参数为基础仿真得到了汽车驾驶员加速度、悬架动挠度和车轮动载荷的仿真结果,这样得到了汽车平顺性与车速、路面不平度、悬架参数的关系。

路面不平度的统计特性-

路面不平度的统计特性-

G42 n G2*4 n Gyx n
16
第二节 路面不平度得统计特性
左、右轮迹间得互谱可以表示为
Gxy n Gxy n e jxy n
两个轮迹得相关函数为
coh
2 xy
n
Gxy n 2 Gxx nGyy n
侧倾角位移功率谱密度Gθ(n)与垂直位移功率谱密度 Gq 得n
比值与相干函数 cohxy n得 关系为
路面不平度的统计特性
第二节 路面不平度得统计特性
一、路面不平度得功率谱密度
1.路面不平度函数
➢路面相对基准平面得高度 q ,沿道路走向长度 I 得变化 q(I)称为路面不平度函数。
➢用水准仪或路面计可以得到路面不平度函数。
2
第二节 路面不平度得统计特性
2.路面不平度得功率谱密度 Gq (n)
1) Gq (n得) 拟合公式
Gq f 2πf 2 Gq f 4π2Gq n0 n02u
Gq f 2πf 4 Gq f 16π4Gq n0 n02uf 2
10
第二节 路面不平度得统计特性
Gq f Gq n0 n02
u f2
对上式得等 式两边取对数 后作图,得到位 移功率谱密度。
u Gq(n0)
Gxy n Gyx n cohxy nGq n
路面对四轮汽车输入得谱矩阵最后可以表示为
1
Gik
n
Gq
n
e j2πnL
coh n
coh
n
e
j2πnL
e j2πnL
1
coh n e j2πnL coh n
coh n coh n e j2πnL
1 e j2πnL
coh n e j2πnL
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
式( 3 ) 可写 成 :
e - 2  ̄ f . m a a At
( , ) 。 ( f ) e - Z , : , d t l
【 3 1

l 2 =
一 2 7 t n m _ i l
2At X

( 4 )
式中: x 为x Ⅲ 的离散傅里 叶变换 ( m= 0 , 1 …N 一 1 o 以上分析推导 出离散 时间信号 与离 散功率谱 密度之 间的相关 关系, 按照式( 4 ) 即可计算得出离散功率谱 密度 G  ̄ ( f 3 。 3采样 与 混淆

5 6・
科 技 论 坛
基于功率谱 密度 的路 面不平整度数值模 拟

任思润
雒 晓 东
( 西 安 市 政设 计研 究 院有 限公 司 . 陕 西 西安 7 1 0 0 6 8 )
摘 要: 路 面不平整是 引起车辆振动进而导致车辆对路 面产生 附加动载 的主要原 因。为了了解路 面不平整对车辆动力特性的影响 , 采 用基 于离散傅里叶逆 变换的方法 , 首先将路面功率谱 密度进行 离散 , 然后模 拟路 面的不平顺情况。模 拟 出的路 面不平整信 号的功率谱
r _ o。 : 、
由式( 6 ) 确定出合适的 △ f 和L 后, 即可得 出频率分辨率为 采样点总数 为 一 。根据傅里叶变换可得 :
△f

l i mX r f t ) =x ( t 1 。 此时, 存在傅里叶变换对:

( , ) = e ( ¨ e - 2  ̄ : , d t
Ⅳ ≥ 壶 ’ 正 号

( 6 )
引 起 激 励 的 位移 或加 速 度 响应 极 小 , 所 以 在 积分 计算 过 程 中 不 需 要 波叠加法 、 积分 单 位 白噪 声 法 、 滤波器整形 白噪声法和 A R MA模 型 考虑这些频率范围的影响 。 设路面有效空间频率范围的上下限分别 法等 , 这些所得到的路面不平整度与实际情况都有较大差别 。本文 为n 、 n 。 , 常用 的车速 为 v , 车辆 固有频率 为 , f 0 , 路 面不平度空 间频 采用 基于离散傅里 叶逆变换 的方 法 , 首先将路面功率谱密度进行离 率为 n , 则频率 f 与 n的关系为 f = v n 。所以: 散, 然后模拟路面的不平顺情况。 2功率谱 密度的离散化
4算 例 分 析 某 汽 车振 动 系 统 的 固有 频率 为 0 . 7 — 1 5 H z , 常 见 车 速 为 1 0 ~ 0 m / s 。现取路面对轮胎激励 的时间频 率上下限分别为 3 0 H z和 ( , ) 二, - l i — m1 0 . 5 Hz 。根据 以上分析可得空 间频 率上下限分别 为 3 m 和 0 . 0 1 m — I 。 以上是概念性的推导 , f 可为正也可为负 ,然而在实际工程 中 f 根 据 式 ( 6 ) 可得 : 为负是没有意义的 , 所 以根据 S ( f ) 偶 函数 的性质 常采用单边 功率 谱密度 G ( , ) 来表示 。 实际情况下 T是有限 的, 即: 本研 究取 A/ =0 . 1 m 、 L = 1 0 0 0 m, 于是 采样 点数 N = 1 0 0 0 0 , 频率 分辨率为 A n 0 0 0 0 01 m- 。
X r ( 列 l i m 1 f 而 ( f ) e f

^ f 去 一 o ・

者 = 1 0 0 m

G B / T 7 0 3 1 — 2 0 0 5按 照功率谱密度 把公路路 面分为 A ~ H 8个等 要把连续模 拟量转换 为离 散数字量 , 需要 对连续模拟 量 x ( 1 1 的 级 ,我 国高等级公 路基 本属 于 A、 B和 c三个等 级 ,本 研究 采用 时间历程进行采样 采样一般都是 以等间隔 取值 ,采样点 数为 MA T L A B软件对 A、 B和 c三个等级 的路 面不 平度 进行 了数 字模 N, 采样 时间为 T (T=N A t) 。则采样频率分辨率为 A t :1: 一, 拟 。模 拟 结 果 如 图 1 ~图 3 所示 。

实践 中经 常遇 到的是非周期函数 ,对于在无限区间 『 — o q + 上 的非周期 函数 x ( t ) , 可以看做是某个周期 函数 ( r ) 当 一士 D 。 时 的情
V 告 1 ,
, =
t 7 , J

况。在 I 一 , l之外按周期T 无限延拓, 构成周期函 数 而( f ) ,
l ( f ) : . J + x r ( , ) e d f ( 1 )
根据随机振动理论可知【 l l , x ( t ) ,t c [ 0 , + c y o ) 的功率谱密度 为 :

J ( 靠 ( = 0 ' l ' 2 … N )
( 8 )
式 中: n e : k A n ; 为随机相位角 , 在 【 o , 2 7 r 】 内均匀分布 。
密度 与 理 论 谱 密 度 几 乎 完 全 吻 合 。 关键词 : 不平整度 ; 功率谱密度 ; 傅里叶 变换 ; 数值模拟
1 概 述
求解上式可得到 :
随着车辆 一路面系统 动力 学分析要求 的逐步提高 和模 拟实 际 情况的需要 , 如何真实的再现路面空间 的不平顺状态成为计算 和试 对具有 N个数据的离散时域信号 x ( n A t ) ,其折叠频率或奈奎 验结果能否符合实 际运行情况 的一个重要保证 。 车辆在运行时会受 到路 面不平整的激励 , 为 了了解路面不平整对车辆动力特性 的影 响 斯特频率为 N 所以 G x ( f J 的 n取值应为 : 0 1 1 - . . 。 就必须从 时域 和频 域两方面进行分析 。 由于汽车隔振系统 的隔离作用 , 使 得车体 自身对某些路 面波形 针对路面不平 整度 的数值模 拟问题 ,目前常使用的方法有 : 谐
相关文档
最新文档