应用连续细化分析傅立叶变换法的发动机激励力测量方法-山东科威数控机床有限公司

专利名称：：应用连续细化分析傅立叶变换法的发动机激励力测量方法
技术领域：
：本方面涉及动力机械和信号处理中发动机激励力的测量，特别是涉及将连续细化分析傅立叶变换信号处理方法应用到发动机振动信号的频率、幅值和相位的校正，获得发动机激励力测量的方法。
背景技术：
：目前发动机的激励力主要是通过刚体动力学模型计算得到，在计算中考虑的主要因素是气缸内气体爆炸力、旋转惯性力和往复惯性力所产生的激励。这种方法的主要缺点是需要获得精确的参数，而有些参数的精确获得是很困难的，如气体爆炸引起的侧倾力矩等。另一方面，在实际情况中，发动机总是与许多附件相连的，而理论计算公式通常并没有考虑发动机与这些附件的耦合，这必然会造成计算结果与工程实际有一定的偏差，而且有些情况下偏差会很大。在实际的工程应用中，通过实验的方法获取发动机激励力为广大学者所研究，这是一个典型的求逆问题，通常是采用频响函数矩阵求逆法。但是在发动机的振动控制问题中，需要的是发动机质心处的等效力和等效力矩，由于无法直接在质心处施加激励，也无法直接获得质心处的响应来应用互易性法，因此无法测得传递函数，导致频响函数矩阵求逆法在发动机激励力测量中无法直接使用。另一方面，对发动机质心处的激励力测量来说，频响函数矩阵求逆法相对来说太复杂。国外的学者J.S.TAO，G.R.LIUandK.Y.LAM等提出了一种测量发动机质心处激励力的方法(J.S.TAO，G.R.LIUandK.Y.LAM于2001在期刊《JournalofSoundandVibration》上发表的文章《Excitationforceidentificationofanenginewithvelocitydataatmountingpoints》)，指出通过提取测量得到的发动机各悬置点处的三向振动速度谱的幅值和相位信息，可精确的重建发动机的激励力，激励力测量公式权利要求应用连续细化分析傅立叶变换法的发动机激励力测量方法，其特征在于包括如下步骤(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><msub><mrow><mo>-</mo><mi>J</mi></mrow><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><msub><mrow><mo>-</mo><mi>J</mi></mrow><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><msub><mrow><mo>-</mo><mi>y</mi></mrow><mn>1</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>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><msub><mrow><mo>-</mo><mi>x</mi></mrow><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><msub><mrow><mo>-</mo><mi>y</mi></mrow><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><msub><mrow><mo>-</mo><mi>x</mi></mrow><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><msub><mrow><mo>-</mo><mi>y</mi></mrow><mn>3</mn></msub><mi></mi></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><msub><mrow><mo>-</mo><mi>x</mi></mrow><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><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><msub><mi>dY</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，进行采样点数N点快速傅立叶变换得到离散频谱，对于采样点数N中的任意第k点，0≤k≤N/2的整数，k点处离散频谱的实部XR(k)和虚部XI(k)为；<mrow><msub><mi>X</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>k</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><mrow><mn>1</mn><mi>n</mi></mrow></msub><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>πkn</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow></mrow>k＝0，1，2…，N/21(11)<mrow><msub><mi>X</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><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><mrow><mn>1</mn><mi>n</mi></mrow></msub><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mi>πkn</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow></mrow>k＝0，1，2…，N/21(12)将式(11)和式(12)通过傅里叶变换得到连续频谱X(f)、连续频谱实部XR(f)、连续频谱虚部XI(f)；<mrow><msub><mi>X</mi><mi>R</mi></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><mrow><mn>1</mn><mi>n</mi></mrow></msub><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>πfn</mi><mo>/</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow>0≤f＜fs/2(13)<mrow><msub><mi>X</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><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><mrow><mn>1</mn><mi>n</mi></mrow></msub><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mi>πfn</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow></mrow>0≤f＜fs/2(14)<mrow><mi>X</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><msub><mi>X</mi><mi>R</mi></msub><mn>2</mn></msup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>+</mo><msup><msub><mi>X</mi><mi>I</mi></msub><mn>2</mn></msup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></msqrt></mrow>0≤f＜fs/2(15)对连续频谱X(f)，根据第η谐次，确定傅里叶变换的频率细化区间[f1，f2]，其中取整数；取整数；对式(13)和式(14)进行q倍的间隔谱分析，q＝100*(f2f1)；确定细化后频率分辨率增量Δf1＝(f2f1)/q(16)确定细化频率序列为{f1，f1+Δf1，f1+2Δf1，…，f1+qΔf1＝f2}(17)对频率细化区间[f1，f2]内的细化频率序列，以序列中的每个频率为变量，代人公式(15)，获得对应频谱值{X(f1)，X(f1+Δf1)，X(f1+2Δf1)，…，X(f2)}，则其中最大值频谱值对应的细化频率定义为f1x，f1x为第η谐次激励频率，用式(13)和式(14)计算得到f1x点的连续频谱实部XR(f1x)和虚部XI(f1x)，再利用极坐标与直角坐标关系式(18)，计算得到f1x频率下的幅值A(f1x)和相位对序列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个传感器的第η谐次激励频率f和对应的位移矩阵<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><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><msub><mi>dY</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><mo>.</mo></mrow>FSA00000275917600011.tif,FSA00000275917600012.tif,FSA00000275917600013.tif,FSA00000275917600014.tif,FSA00000275917600015.tif,FSA00000275917600025.tif,FSA00000275917600026.tif,FSA00000275917600027.tif,FSA00000275917600034.tif,FSA00000275917600035.tif,FSA00000275917600036.tif,FSA00000275917600037.tif,FSA00000275917600038.tif,FSA00000275917600039.tif,FSA000002759176000310.tif,FSA000002759176000311.tif,FSA000002759176000312.tif2.根据权利要求1所述的应用连续细化分析傅立叶变换法的发动机激励力测量方法，其特征在于所述曲轴自由端为曲轴连接皮带轮的一端。3.根据权利要求1所述的应用连续细化分析傅立叶变换法的发动机激励力测量方法，其特征在于所述发动机与支撑它的悬置接触点定义为悬置点。4.根据权利要求1所述的应用连续细化分析傅立叶变换法的发动机激励力测量方法，其特征在于所述正常工作转速w为750r/min至5500r/min；所述工况稳定是指发动机转速波动在20r/min以内。5.根据权利要求1所述的应用连续细化分析傅立叶变换法的发动机激励力测量方法，其特征在于所述的发动机转动惯量实验台采用三线摆测试法实验台。全文摘要本发明公开了一种基于振动响应的发动机激励力测量方法，涉及动力机械和信号处理领域，包括激励力测量公式和连续细化分析傅立叶变换法，该方法对发动机缸体表面或悬置点上的振动加速度信号作全景谱，对全景谱中的发动机激励频率区间，采用傅立叶变换(FT)细化运算，精确获得振动位移频谱中发动机真实激励频率f，以及幅值和相位信息，求出激励力测量公式中的振动位移矢量精确值，结合其它已知的发动机参数，计算出激励力F(f)。本发明的方法可以解决发动机激励力测量中因缺乏准确相位信息，而需要求解复杂非线性方程的问题。文档编号G01L5/00GK101975633SQ20101028626公开日2011年2月16日申请日期2010年9月16日优先权日2010年9月16日发明者丁康,徐传燕,杨志坚,杨翠丽,王永胜申请人:华南理工大学