应用改变窗长相位差校正法的发动机激励力测量方法-山东科威数控机床有限公司

专利名称：：应用改变窗长相位差校正法的发动机激励力测量方法
技术领域：
：本方面涉及动力机械和信号处理领域，具体是一种精确识别发动机振动信号的频率、幅值、相位方法，应用于发动机激励力识别。
背景技术：
：目前发动机的激励力主要是通过刚体动力学模型计算得到，在计算中考虑的主要因素是气缸内气体爆炸力、旋转惯性力和往复惯性力所产生的激励。这种方法的主要缺点是需要获得精确的参数，而有些参数的精确获得是很困难的，如气体爆炸引起的侧倾力矩等。另一方面，在实际情况中，发动机总是与许多附件相连的，而理论计算公式通常并没有考虑发动机与这些附件的耦合，这必然会造成计算结果与工程实际有一定的偏差，而且有些情况下偏差会很大。在实际的工程应用中，通过实验的方法获取发动机激励力为广大学者所研究，这是一个典型的求逆问题，通常是采用频响函数矩阵求逆法。但是在发动机的振动控制问题中，需要的是发动机质心处的等效力和等效力矩，由于无法直接在质心处施加激励，也无法直接获得质心处的响应来应用互易性法，因此无法测得传递函数，导致频响函数矩阵求逆法在发动机激励力测量中无法直接使用。另一方面，对发动机质心处的激励力测量来说，频响函数矩阵求逆法相对来说太复杂。国外的学者J.S.TAO，G.R.LIUandK.Y.LAM等提出了一种测量发动机质心处激励力的方法(J.S.TAO，G.R.LIUandK.Y.LAM于2001在期刊《JournalofSoundandVibration》上发表的文章《Excitationforceidentificationofanenginewithvelocitydataatmountingpoints》)，指出通过提取测量得到的发动机各悬置点处的三向振动速度谱的幅值和相位信息，可精确的重建发动机的激励力，激励力测量公式V{f)=[K'-{iTtfjM'iEEYE'sif)。但是该方法在提取相位信息时，通常无法获得绝对的相位信息，因此提出了只利用相位差来获得激励力的方法，而这种处理引入了与测量悬置点数相同的相位变量，使得问题的求解变成了求解一个复杂的非线性超定方程，增加了求解的复杂性。
发明内容本发明的目的是针对现有识别方法的缺陷，提出一种应用改变窗长相位差校正法的发动机激励力测量方法，解决发动机激励力测量中因缺乏准确相位信息，而需要求解复杂非线性方程的问题。本发明目的通过如下技术方案实现—种应用改变窗长相位差校正法的发动机激励力测量方法，包括如下步骤(1)坐标系建立与相关参数收集利用发动机转动惯量实验台，测试出发动机质量m，发动机质心0，以0为原点建立坐标系0-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；在发动机的坐标系0-XYZ中，按照X方向由大到小的顺序给悬置标号，定义为第一悬置、第二悬置、…第h悬置，用弹性体测试系统对各悬置的复刚度特性进行测试，获得第一悬置的复刚度、第二悬置的复刚度、…第h悬置的复刚度，组成悬置复刚度矩阵K*=K*2…！^]；(3)安装传感器在发动机缸体表面或者悬置点的不同测试点上，安装L个三向加速度传感器，L^3，每个测试点安装一个三向加速度传感器；传感器连接数据采集器，数据采集器连接便携式计算机；在发动机的坐标系0-XYZ中，按照X方向由大到小的顺序确定三向加速度传感器的坐标分别为hZl]、[x2y2z2]、…、[xLyLzL]；(4)在发动机试验台架上，设定发动机以正常工作转速W运行，待工况稳定后，在10s30s内的采集时间，通过数据采集器采集和同步记录各测试点的X轴正向、Y轴正向和z轴正向的振动加速度信号Xln、X2n-XLn；Yln、Y2n-YLn；Zln、Z2n…ZLn；其中Xln为第一测试点X轴正向采集的具有n个点的振动加速度信号序列，XLn为第L测试点X轴正向采集的n个点的振动加速度信号序列，Yln为第一测试点Y轴正向采集的n个点的振动加速度信号序列，第L测试点Y轴正向采集的n个点的振动加速度信号序列；Zln为第一测试点Z轴正向采集的n个点的振动加速度信号序列，第L测试点Z轴正向采集的n个点的振动加速度信号序列；采样频率fs，采样点数N表示离散序列包含的数据点个数，其中n=0、1、…、N-1，要求fs=N=2P，P=9或10，频率分辨率Af=fs/N=1；(6)根据公式F(/)=[火‘一位；^拟！^才；^⑴测量发动机激励力，其中:=Fy{f)Fz(f)Mx{f)My{f)Mz(/)j，激励力F(f)是指发动机在转速w第n谐次激励频率f下，发动机X轴正向的激励*Fx(f)、绕X轴的激励转矩Mx(f)、Y轴方向的激励力Fy(f)、绕Y轴的激励转矩My(f)、Z轴方向的激励力Fz(f)、绕Z轴的激励转矩Mz(f)；权利要求应用改变窗长相位差校正法的发动机激励力测量方法，其特征在于包括如下步骤(1)坐标系建立与相关参数收集利用发动机转动惯量实验台，测试出发动机质量m，发动机质心O，以O为原点建立坐标系OXYZ，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；在发动机的坐标系OXYZ中，按照X方向由大到小的顺序给悬置标号，定义为第一悬置、第二悬置、…第h悬置，用弹性体测试系统对各悬置的复刚度特性进行测试，获得第一悬置的复刚度第二悬置的复刚度…第h悬置的复刚度组成悬置复刚度矩阵K*＝[K*1K*2…K*h]；(3)安装传感器在发动机缸体表面或者悬置点的不同测试点上，安装L个三向加速度传感器，L≥3，每个测试点安装一个三向加速度传感器；传感器连接数据采集器，数据采集器连接便携式计算机；在发动机的坐标系OXYZ中，按照X方向由大到小的顺序确定三向加速度传感器的坐标分别为[x1y1z1]、[x2y2z2]、…、[xLyLzL]；(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、…、N1，要求fs＝N＝2p，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*＝[K*1K*2…K*h]<mrow><mi>M</mi><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mi>m</mi></mtd><mtd></mtd><mtd></mtd><mtd></mtd><mtd></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd><mi>m</mi></mtd><mtd></mtd><mtd></mtd><mtd></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd></mtd><mtd><mi>m</mi></mtd><mtd></mtd><mtd></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd></mtd><mtd></mtd><mtd><msub><mi>J</mi><mi>x</mi></msub></mtd><mtd><mo>-</mo><msub><mi>J</mi><mi>xy</mi></msub></mtd><mtd><mo>-</mo><msub><mi>J</mi><mi>zx</mi></msub></mtd></mtr><mtr><mtd></mtd><mtd></mtd><mtd></mtd><mtd><mo>-</mo><msub><mi>J</mi><mi>xy</mi></msub></mtd><mtd><msub><mi>J</mi><mi>y</mi></msub></mtd><mtd><mo>-</mo><msub><mi>J</mi><mi>yz</mi></msub></mtd></mtr><mtr><mtd></mtd><mtd></mtd><mtd></mtd><mtd><msub><mrow><mo>-</mo><mi>J</mi></mrow><mi>zx</mi></msub></mtd><mtd><mo>-</mo><msub><mi>J</mi><mi>yz</mi></msub></mtd><mtd><msub><mi>J</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced></mrow><mrow><mover><mi>E</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>z</mi><mn>1</mn></msub></mtd><mtd><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>x</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><msub><mi>y</mi><mn>1</mn></msub></mtd><mtd><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>z</mi><mn>2</mn></msub></mtd><mtd><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mrow><mo>-</mo><mi>z</mi></mrow><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>x</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><msub><mi>y</mi><mn>2</mn></msub></mtd><mtd><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd></mtr><mtr><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd></mtr><mtr><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>z</mi><mn>3</mn></msub></mtd><mtd><mo>-</mo><msub><mi>y</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mrow><mo>-</mo><mi>z</mi></mrow><mn>3</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>x</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><msub><mi>y</mi><mn>3</mn></msub></mtd><mtd><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow><mrow><mi>f</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>3</mn><mi>L</mi></mrow></mfrac><mrow><mo>(</mo><msub><mi>f</mi><mrow><mn>1</mn><mi>x</mi></mrow></msub><mo>+</mo><msub><mi>f</mi><mrow><mn>1</mn><mi>y</mi></mrow></msub><mo>+</mo><msub><mi>f</mi><mrow><mn>1</mn><mi>z</mi></mrow></msub><mo>+</mo><msub><mi>f</mi><mrow><mn>2</mn><mi>x</mi></mrow></msub><mo>+</mo><msub><mi>f</mi><mrow><mn>2</mn><mi>y</mi></mrow></msub><mo>+</mo><msub><mi>f</mi><mrow><mn>2</mn><mi>z</mi></mrow></msub><mo>+</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>+</mo><msub><mi>f</mi><mi>Lx</mi></msub><mo>+</mo><msub><mi>f</mi><mi>Ly</mi></msub><mo>+</mo><msub><mi>f</mi><mi>Lz</mi></msub><mo>)</mo></mrow></mrow><mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msup><mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>dX</mi><mn>1</mn></msub></mtd><mtd><mi>d</mi><msub><mi>Y</mi><mn>1</mn></msub></mtd><mtd><msub><mi>dZ</mi><mn>1</mn></msub></mtd><mtd><msub><mi>dX</mi><mn>2</mn></msub></mtd><mtd><msub><mi>dY</mi><mn>2</mn></msub></mtd><mtd><msub><mi>dZ</mi><mn>2</mn></msub><mo>·</mo><mo>·</mo><mo>·</mo><msub><mi>dX</mi><mi>L</mi></msub></mtd><mtd><msub><mi>dY</mi><mi>L</mi></msub></mtd><mtd><msub><mi>dZ</mi><mi>L</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>表示矩阵的转置；第η谐次激励频率f和其对应位移矩阵通过如下方法确定选择发动机激励力分析的谐次η，η＝0.5、1、1.5、2、2.5或3；指定分析发动机转速w的第η谐次后，进行以下运算对第一测试点X轴正向采集的具有n个点的振动加速度信号序列X1n，加时间长度T＝1的汉宁窗w(t)0≤t≤T，按照公式(11)进行N点快速傅里叶变换，得到幅值谱X1(f)；对X1n加时间长度汉宁窗w(t)0≤t≤T，按照公式(12)进行N/2点快速傅里叶变换(FFT)，得到幅值谱X2(f)；<mrow><msub><mi>X</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>x</mi><mi>n</mi></msub><mo>[</mo><mn>0.5</mn><mo>-</mo><mn>0.5</mn><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>πn</mi></mrow><mi>N</mi></mfrac><mo>)</mo></mrow><msup><mrow><mo>]</mo><mi>e</mi></mrow><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>N</mi></mfrac><mi>fn</mi></mrow></msup><mo>=</mo><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>jI</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>y</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>Φ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></msup></mrow><mrow><mi>f</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>X</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>x</mi><mi>n</mi></msub><mo>[</mo><mn>0.5</mn><mo>-</mo><mn>0.5</mn><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>4</mn><mi>πn</mi></mrow><mi>N</mi></mfrac><mo>)</mo></mrow><msup><mrow><mo>]</mo><mi>e</mi></mrow><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>π</mi></mrow><mi>N</mi></mfrac><mi>fn</mi></mrow></msup><mo>=</mo><msub><mi>R</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>jI</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>y</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>Φ</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></msup></mrow><mrow><mi>f</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mn>2</mn><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mfrac><mi>N</mi><mn>4</mn></mfrac><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><mo>)</mo></mrow></mrow>根据发动机转速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的差值代人公式(13)，求出频率误差<mrow><mo>&dtri;</mo><mi>f</mi><mo>=</mo><mo>-</mo><mn>2</mn><msub><mrow><mo>[</mo><mi>Φ</mi></mrow><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>Φ</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mi>π</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>f</mi><mi>i</mi></msub><mo>-</mo><msub><mi>f</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>]</mo><mo>/</mo><mi>π</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>3</mn><mo>)</mo></mrow></mrow>计算真实频率f0计算真实相位汉宁窗的谱函数根据公式(16)计算出校正后的幅值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个传感器的激励频率以及对应的位移矩阵FSA00000278397300011.tif,FSA00000278397300012.tif,FSA00000278397300013.tif,FSA00000278397300014.tif,FSA00000278397300015.tif,FSA00000278397300025.tif,FSA00000278397300026.tif,FSA00000278397300027.tif,FSA00000278397300028.tif,FSA00000278397300033.tif,FSA00000278397300034.tif,FSA00000278397300035.tif,FSA00000278397300036.tif,FSA00000278397300037.tif,FSA00000278397300038.tif,FSA000002783973000310.tif,FSA000002783973000311.tif,FSA000002783973000312.tif,FSA000002783973000313.tif,FSA000002783973000314.tif,FSA000002783973000315.tif,FSA000002783973000316.tif,FSA000002783973000317.tif,FSA000002783973000318.tif,FSA00000278397300041.tif,FSA00000278397300042.tif,FSA00000278397300043.tif2.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法，其特征在于所述的曲轴自由端为曲轴连接皮带轮的一端。3.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法，其特征在于所述的发动机与支撑它的悬置接触点定义为悬置点。4.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法，其特征在于所述的正常工作转速w为750r/min至5500r/min；所述工况稳定是指发动机转速波动在20r/min以内。5.根据权利要求1所述的应用改变窗长相位差校正法的发动机激励力测量方法，其特征在于所述的发动机转动惯量实验台采用三线摆测试法实验台。全文摘要本发明公开了应用改变窗长相位差校正法的发动机激励力测量方法，该方法是对发动机缸体表面或悬置点上的振动加速度信号离散采样之后，利用在信号处理中加不同时间长度的汉宁窗(Hanning)而产生的相位差，求出发动机激励频率区间的信号频率误差用求出振动位移频谱中发动机真实激励频率f，以及幅值和相位信息，从而精确获得激励力测量公式中振动位移矢量值，结合其它已知的发动机参数，计算出激励力F(f)。本发明的提出可以解决发动机激励力测量中因缺乏准确相位信息，而需要求解复杂非线性方程的问题。文档编号G01L5/00GK101975634SQ20101028808公开日2011年2月16日申请日期2010年9月17日优先权日2010年9月17日发明者丁康,张绍辉,李魏华,杨志坚,王永胜申请人:华南理工大学