首页 / 专利库 / 保养与维修 / 发动机试验台架 / 应用改变窗长相位差校正法的发动机激励力测量方法

应用改变窗长相位差校正法的发动机激励测量方法

阅读:351发布:2020-08-27

专利汇可以提供应用改变窗长相位差校正法的发动机激励测量方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了应用改变窗长 相位 差校正法的 发动机 激励 力 测量方法,该方法是对发动机缸体表面或悬置点上的振动 加速 度 信号 离散 采样 之后,利用在 信号处理 中加不同时间长度的汉宁窗(Hanning)而产生的 相位差 ,求出发动机激励 频率 区间的信号频率误差用求出振动位移 频谱 中发动机真实激励频率f,以及幅值和相位信息,从而精确获得激励力测量公式中振动位移矢量值,结合其它已知的发动机参数,计算出激励力F(f)。本发明的提出可以解决发动机激励力测量中因缺乏准确相位信息,而需要求解复杂非线性方程的问题。,下面是应用改变窗长相位差校正法的发动机激励测量方法专利的具体信息内容。

1.应用改变窗长相位差校正法的发动机激励测量方法,其特征在于包括如下步骤:
(1)坐标系建立与相关参数收集:利用发动机转动惯量实验台,测试出发动机质量m,发动机质心O,以O为原点建立坐标系O-XYZ,Y轴正向指向曲轴自由端,Z轴正向竖直向上,按右手定则确定X轴正向;再利用发动机转动惯量实验台测试出发动机绕坐标系X轴的转动惯量Jx、Y轴的转动惯量Jy、Z轴的转动惯量Jz、发动机对X轴和Y轴的惯性积Jxy、对Y轴和Z轴的惯性Jyz、对Z轴和X轴的惯性Jzx,组成发动机质量矩阵M;
(2)获取悬置复刚度矩阵:将发动机通过悬置支撑发动机试验台架;确定发动机悬置个数h,h=3或h=3;在发动机的坐标系O-XYZ中,按照X方向由大到小的顺序给悬置标号,定义为第一悬置、第二悬置、…第h悬置,用弹性体测试系统对各悬置的复刚度特性进行测试,获得第一悬置的复刚度 第二悬置的复刚度 …第h悬置的复刚度 组成* * * *
悬置复刚度矩阵K =[K1 K2 …Kh ];
(3)安装传感器:在发动机缸体表面或者悬置点的不同测试点上,安装L个三向加速度传感器,L≥3,每个测试点安装一个三向加速度传感器;传感器连接数据采集器,数据采集器连接便携式计算机;在发动机的坐标系O-XYZ中,按照X方向由大到小的顺序确定三向加速度传感器的坐标分别为[x1 y1 z1]、[x2 y2 z2]、…、[xL yL zL];
(4)在发动机试验台架上,设定发动机以正常工作转速w运行,待工况稳定后,在10s~
30s内的采集时间,通过数据采集器采集和同步记录各测试点的X轴正向、Y轴正向和Z轴正向的振动加速度信号X1n、X2n…XLn;Y1n、Y2n…YLn;Z1n、Z2n…ZLn;其中X1n为第一测试点X轴正向采集的具有n个点的振动加速度信号序列,XLn为第L测试点X轴正向采集的n个点的振动加速度信号序列,Y1n为第一测试点Y轴正向采集的n个点的振动加速度信号序列,XLn第L测试点Y轴正向采集的n个点的振动加速度信号序列;Z1n为第一测试点Z轴正向采集的n个点的振动加速度信号序列,ZLn第L测试点Z轴正向采集的n个点的振动加速度信号序列;采样频率fs,采样点数N表示离散序列包含的数据点个数,其中n=0、1、…、N-1,要p
求fs=N=2,P=9或10,频率分辨率Δf=fs/N=1;
(5)根据公式 测量发动机激励力,其中:
激励力F(f)是指发动机在转
速w第η谐次激励频率f下,发动机X轴正向的激励力Fx(f)、绕X轴的激励转矩Mx(f)、Y轴方向的激励力Fy(f)、绕Y轴的激励转矩My(f)、Z轴方向的激励力Fz(f)、绕Z轴的激励转矩Mz(f);
* * * *
K =[K1 K2 …Kh]
表示矩阵 的转置;
第η谐次激励频率f和其对应位移矩阵 通过如下方法确定:
选择发动机激励力分析的谐次η,η=0.5、1、1.5、2、2.5或3;指定分析发动机转速w的第η谐次后,进行以下运算:
对第一测试点X轴正向采集的具有n个点的振动加速度信号序列X1n,加时间长度T=
1的汉宁窗w(t)0≤t≤T,按照公式(1-1)进行N点快速傅里叶变换,得到幅值谱X1(f);
对X1n加时间长度 汉宁窗w(t)0≤t≤T,按照公式(1-2)进行N/2点快速傅里叶变换(FFT),得到幅值谱X2(f);
根据发动机转速w,确定校正频率区间[f1,f2],其中 取整数; 取
整数;
在校正频率区间[f1,f2]内,X1(f)的最高峰值对应频率fi,fi实部为R1(fi),虚部为I1(fi),相位为Φ1(fi),幅值为y1(fi);X2(f)的最高峰值频率fj,相位为Φ2(fj);设校正频率区间[f1,f2]内最高峰值真实频率为f0,f0对应真实相位 真实幅值A0(f0);设fi的频率误差为 fi与真实频率f0的差值 代人公式(1-3),求出频率误差
计算真实频率f0
计算真实相位
汉宁窗的谱函数 根据公式(1-6)计算出校正后的幅值
A0(f0),组成真实频率f0和对应的加速度矢量
对序列X2n…XLn;Y1n、Y2n…YLn;Z1n、Z2n…ZLn与振动加速度信号序列X1n进行同样的处理,分别得到发动机转速w第η谐次的激励频率f2x…fLx;f1y、f2y…fLy;f1z、f2z…fLz;
对于L个传感器测试点中的任意传感器H(1≤H≤L),测得它的X方向的发动机激
励频率为fHx,对应加速度矢量为 测得Y方向的发动机激励频率为fHy,对
应加速度矢量为 测得Z方向的发动机激励频率为fHz,对应加速度矢量为
传感器H的X方向位移为dXH,Y方向位移为dYH,Z方向位移为dZH,由此获
得传感器H的位移矩阵:
L个传感器的激励频率 以及
对应的位移矩阵
2.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法,其特征在于:所述的曲轴自由端为曲轴连接皮带轮的一端。
3.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法,其特征在于:所述的发动机与支撑它的悬置接触点定义为悬置点。
4.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法,其特征在于:所述的正常工作转速w为750r/min至5500r/min;所述工况稳定是指发动机转速波动在20r/min以内。
5.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法,其特征在于:所述的发动机转动惯量实验台采用三线摆测试法实验台。

说明书全文

应用改变窗长相位差校正法的发动机激励测量方法

技术领域

[0001] 本方面涉及动力机械和信号处理领域,具体是一种精确识别发动机振动信号的频率、幅值、相位方法,应用于发动机激励力识别。

背景技术

[0002] 目前发动机的激励力主要是通过刚体动力学模型计算得到,在计算中考虑的主要因素是气缸内气体爆炸力、旋转惯性力和往复惯性力所产生的激励。这种方法的主要缺点是需要获得精确的参数,而有些参数的精确获得是很困难的,如气体爆炸引起的侧倾力矩等。另一方面,在实际情况中,发动机总是与许多附件相连的,而理论计算公式通常并没有考虑发动机与这些附件的耦合,这必然会造成计算结果与工程实际有一定的偏差,而且有些情况下偏差会很大。
[0003] 在实际的工程应用中,通过实验的方法获取发动机激励力为广大学者所研究,这是一个典型的求逆问题,通常是采用频响函数矩阵求逆法。但是在发动机的振动控制问题中,需要的是发动机质心处的等效力和等效力矩,由于无法直接在质心处施加激励,也无法直接获得质心处的响应来应用互易性法,因此无法测得传递函数,导致频响函数矩阵求逆法在发动机激励力测量中无法直接使用。另一方面,对发动机质心处的激励力测量来说,频响函数矩阵求逆法相对来说太复杂。
[0004] 国外的学者J.S.TAO,G.R.LIU and K.Y.LAM等提出了一种测量发动机质心处激励力的方法(J.S.TAO,G.R.LIU and K.Y.LAM于2001在期刊《Journal of Sound and Vibration》上发表的文章《Excitation force identification of an engine with velocity data at mounting points》),指出通过提取测量得到的发动机各悬置点处的三向振动速度谱的幅值和相位信息,可精确的重建发动机的激励力,激励力测量公式但是该方法在提取相位信息时,通常无法获得绝对的相位信息,因此提出了只利用相位差来获得激励力的方法,而这种处理引入了与测量悬置点数相同的相位变量,使得问题的求解变成了求解一个复杂的非线性超定方程,增加了求解的复杂性。
发明内容
[0005] 本发明的目的是针对现有识别方法的缺陷,提出一种应用改变窗长相位差校正法的发动机激励力测量方法,解决发动机激励力测量中因缺乏准确相位信息,而需要求解复杂非线性方程的问题。
[0006] 本发明目的通过如下技术方案实现:
[0007] 一种应用改变窗长相位差校正法的发动机激励力测量方法,包括如下步骤:
[0008] (1)坐标系建立与相关参数收集:利用发动机转动惯量实验台,测试出发动机质量m,发动机质心O,以O为原点建立坐标系O-XYZ,Y轴正向指向曲轴自由端,Z轴正向竖直向上,按右手定则确定X轴正向;再利用发动机转动惯量实验台测试出发动机绕坐标系X轴的转动惯量Jx、Y轴的转动惯量Jy、Z轴的转动惯量Jz、发动机对X轴和Y轴的惯性积Jxy、对Y轴和Z轴的惯性Jyz、对Z轴和X轴的惯性Jzx组成发动机质量矩阵M;
[0009] (2)获取悬置复刚度矩阵:将发动机通过悬置支撑发动机试验台架;确定发动机悬置个数h,h=3或h=3;在发动机的坐标系O-XYZ中,按照X方向由大到小的顺序给悬置标号,定义为第一悬置、第二悬置、…第h悬置,用弹性体测试系统对各悬置的复刚度特性进行测试,获得第一悬置的复刚度 第二悬置的复刚度 …第h悬置的复刚度组成悬置复刚度矩阵K*=[K*1 K*2 …K*h ];
[0010] (3)安装传感器:在发动机缸体表面或者悬置点的不同测试点上,安装L个三向加速度传感器,L≥3,每个测试点安装一个三向加速度传感器;传感器连接数据采集器,数据采集器连接便携式计算机;在发动机的坐标系O-XYZ中,按照X方向由大到小的顺序确定三向加速度传感器的坐标分别为[x1 y1 z1]、[x2 y2 z2]、…、[xL yL zL];
[0011] (4)在发动机试验台架上,设定发动机以正常工作转速w运行,待工况稳定后,在10s~30s内的采集时间,通过数据采集器采集和同步记录各测试点的X轴正向、Y轴正向和Z轴正向的振动加速度信号X1n、X2n…XLn;Y1n、Y2n…YLn;Z1n、Z2n…ZLn;其中X1n为第一测试点X轴正向采集的具有n个点的振动加速度信号序列,XLn为第L测试点X轴正向采集的n个点的振动加速度信号序列,Y1n为第一测试点Y轴正向采集的n个点的振动加速度信号序列,XLn第L测试点Y轴正向采集的n个点的振动加速度信号序列;Z1n为第一测试点Z轴正向采集的n个点的振动加速度信号序列,ZLn第L测试点Z轴正向采集的n个点的振动加速度信号序列;采样频率fs,采样点数N表示离散序列包含的数据点个数,其中n=0、1、…、N-1,要求fs=N=2p,P=9或10,频率分辨率Δf=fs/N=1;
[0012] (6)根据公式 测量发动机激励力,其中:激励力F(f)是指发动机在转速
w第η谐次激励频率f下,发动机X轴正向的激励力Fx(f)、绕X轴的激励转矩Mx(f)、Y轴方向的激励力Fy(f)、绕Y轴的激励转矩My(f)、Z轴方向的激励力Fz(f)、绕Z轴的激励转矩Mz(f);
[0013] K*=[K*1 K*2…K*h]
[0014]
[0015]
[0016]
[0017]
[0018] 表示矩阵 的转置;
[0019] 第η谐次激励频率f和其对应位移矩阵 通过如下方法确定:
[0020] 选择发动机激励力分析的谐次η,η=0.5、1、1.5、2、2.5或3;指定分析发动机转速w的第η谐次后,进行以下运算:
[0021] 对第一测试点X轴正向采集的具有n个点的振动加速度信号序列X1n,加时间长度T=1的汉宁窗(Hanning)w(t)0≤t≤T,按照公式(1-1)进行N点快速傅里叶变换(FFT),得到幅值谱X1(f);对X1n加时间长度 汉宁窗(Hanning)w(t)0≤t≤T,按照公式(1-2)进行N/2点快速傅里叶变换(FFT),得到幅值谱X2(f);
[0022]
[0023]
[0024] 根据发动机转速w,确定校正频率区间[f1,f2],其中 取整数;取整数;
[0025] 在校正频率区间[f1,f2]内,X1(f)的最高峰值对应频率fi,fi实部为R1(fi),虚部为I1(fi),相位为Φ1(fi),幅值为y1(fi);X2(f)的最高峰值频率fj,相位为Φ2(fj);设校正频率区间[f1,f2]内最高峰值真实频率为f0,f0对应真实相位 真实幅值A0(f0);设fi的频率误差为 fi与真实频率f0的差值 代人公式(1-3),求出频率误差
[0026]
[0027] 计 算 真 实 频 率f0
[0028] 计 算 真 实 相 位
[0029] 汉宁窗(Hanning)的谱函数 根据 公式(1-6)计算出校正后的幅值A0(f0),组成真实频率f0和对应的加速度矢量
[0030] 对序列X2n…XLn;Y1n、Y2n…YLn;Z1n、Z2n…ZLn与振动加速度信号序列X1n进行同样的处理,分别得到发动机转速w第η谐次的激励频率f2x…fLx;f1y、f2y…fLy;f1z、f2z…fLz;
[0031] 对于L个传感器测试点中的任意传感器H(1≤H≤L),测得它的X方向的发动机激励频率为fHx,对应加速度矢量为 测得Y方向的发动机激励频率为fHy,对应加速度矢量为 测得Z方向的发动机激励频率为fHz,对应加速度矢量为
传感器H的X方向位移为dXH,Y方向位移为dYH,Z方向位移为dZH,
[0032] 由此获得传感器H的位移矩阵:
[0033]
[0034] L个传感器的激励频率以及对应的位移矩阵
[0035] 为进一步实现本发明目的,所述的曲轴自由端为曲轴连接皮带轮的一端。
[0036] 所述的发动机与支撑它的悬置接触点定义为悬置点。
[0037] 所述的正常工作转速w优选为750r/min至5500r/min;所述工况稳定是指发动机转速波动在20r/min以内。
[0038] 所述的发动机转动惯量实验台优选采用三线摆测试法实验台。
[0039] 本发明对发动机不同谐次η的激励力测试,只需通过变换η的取值,η=0.5、1、1.5、2、2.5、3;计算出不同的频率细化区间[f1,f2],重新计算出该谐次的激励频率f和对应位移矩阵 求解出发动机转速w的第η谐次激励力F(f)。
[0040] 相对于现有技术,本发明具有如下优点:
[0041] 本发明通过对发动机缸体表面或悬置点上的振动位移信号离散采样之后,加不同时间长度的汉宁窗(Hanning)而产生的相位差,利用该相位差校正各响应点振动信号的精确频率、幅值和相位信息,从而提高激励力的测试精度。本发明结合J.S.TAO,G.R.LIU and *K.Y.LAM提出的激励力测量公式 中K、M、求解方法,
提出应用改变窗长相位差校正法精确提取位移矩阵 中的频率、幅值和相位信息,改进了利用相位差方法来计算 的方法,避免了因利用相位差方法,将问题的求解变成了求解一个复杂的非线性超定方程,提高了求解的 的精度和效率,并提高了F(f)求解的简便性。
具体实施方式:
[0042] 下面结合实施例对本发明作进一步的描述,需要说明的是,实施例并不构成对本发明要求保护范围的限定。
[0043] 实施例1
[0044] (1)利用三线摆式发动机转动惯量实验台测试出:某1.8L直列四缸发动机质量m,发动机质心O,以质心O为原点建立坐标系O-XYZ;然后再测试出发动机绕坐标系X轴的转动惯量Jx、Y轴的转动惯量Jy、Z轴的转动惯量Jz、发动机对X轴和Y轴的惯性积Jxy、对Y轴和Z轴的惯性Jyz、对Z轴和X轴的惯性Jzx,组成发动机质量矩阵M;
[0045]
[0046] (2)确定发动机悬置个数h,h=3;在发动机的坐标系O-XYZ中,按照X方向由大到小的顺序给悬置标号,定义为悬置1、悬置2、…悬置h,采用德国MTS831型弹性体测试系统测试出各悬置的复刚度,组成悬置复刚度矩阵K*=[k*1 K*2 K*3]=[(75000+7500j)N/m(70000+7000j)N/m(80000+8000j)N/m];
[0047] (3)把1.8L直列四缸发动机安装在德国FEV发动机试验台架上,采用CPB三向加速度传感器,个数L=3,布置在悬置点处。通过三坐标测量仪测出3个传感器测试点的坐标,组成测试点坐标矩阵
[0048]
[0049] (4)CPB三向加速度传感器输入是德国米勒贝姆(BBM)测试设备的数据采集系统,数据采集信号输入便携式计算机,通过便携式计算机上安装的PAK测试分析软件,实时观测和记录信号;
[0050] (5)在发动机试验台架上,设定发动机以某正常工作转速w=1815rpm运行,待工况稳定后(即转速波动在20r/min以内),开始采集数据,CPB三向加速度传感器将信号输入德国米勒贝姆(BBM)测试设备的数据采集系统,通过便携式计算机上安装的PAK测试分析软件,实时观测和记录信号。采集时间长度为30s;采样频率fs=512,采样点数N=512,则频率分辨率Δf=fs/N=1,通过数据采集器采集和同步记录各测试点的X轴正向、Y轴正向和Z轴正向的振动加速度信号X1n、X2n、X3n;Y1n、Y2n、Y3n;Z1n、Z2n、Z3n;其中n=0、1、…、N-1;对振动加速度信号X1n、X2n、X3n;Y1n、Y2n、Y3n;Z1n、Z2n、Z3n进行N=512点快速傅立叶变换
[0051] 根据发动机转速w=1815rpm,确定频率区间[f1,f2],发动机转速w的0.5谐次激励频率f的细化区间为f1=10HZ和f2=20HZ;1谐次的细化区间为f1=25HZ和f2=35HZ;1.5谐次的细化区间为f1=40HZ和f2=50HZ;2谐次的细化区间为f1=55HZ和f2=65HZ;2.5谐次的细化区间为f1=71HZ和f2=81HZ;3谐次的细化区间为f1=86HZ和f2=96HZ;对每个频率区间[f1,f2],使用改变窗长相位差法,具体步骤按照公式(1-1)至(1-7)进行。
[0052] 3)由于篇幅关系,这里只列举出发动机转速w的1谐次和2谐次的激励频率f,发动机转速w的1谐次频率区间[f1,f2]=[25,35],使用改变窗长相位差法,获得3个传感器对应各自测试点的激励频率f与加速度矢量 如表1所示;发动机转速w的2谐次频率区间[f1,f2]=[55,65],使用改变窗长相位差法,获得3个传感器对应各自测试点的激励频率f与加速度矢量 如表2所示;
[0053] 表1发动机转速w的1谐次频率与加速度矢量
[0054]
[0055]
[0056] 表2发动机转速w的2谐次频率与加速度矢量
[0057]
[0058] 由此获得:当 时,
[0059] 位移矩阵 其中:
[0060]
[0061]
[0062]
[0063] 当 时,
[0064] 位移矩阵 其中:
[0065]
[0066]
[0067]
[0068] (7)根据发动机激励力测量的基本公式 其中: 代人(1)中的质量矩阵M,
*
(2)中的悬置复刚度矩阵K,(3)中的测点坐标矩阵 当求发动机转速w的1谐次激励力时,代人f=30.25HZ和 求
出发动机发动机转速w的1谐次激励力如表3所示;当求发动机转速w的2谐次激励力时,代人f=60.5HZ和 求出发动
机发动机转速w的2谐次激励力如表4所示;
[0069] 表3发动机转速w的1谐次激励力
[0070]
[0071] 表4发动机转速w的2谐次激励力
[0072]
[0073]
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈