CN107356387A - 一种模态试验中多传感器附加质量消除方法 - Google Patents
一种模态试验中多传感器附加质量消除方法 Download PDFInfo
- Publication number
- CN107356387A CN107356387A CN201710598555.4A CN201710598555A CN107356387A CN 107356387 A CN107356387 A CN 107356387A CN 201710598555 A CN201710598555 A CN 201710598555A CN 107356387 A CN107356387 A CN 107356387A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtr
- mtd
- msup
- 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
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种模态试验中多传感器附加质量消除方法,推导模态试验中多传感器质量影响消除公式,确定模态试验中布置的传感器个数及消除传感器质量的顺序,基于实测的频响函数根据消除公式逐一计算消除各传感器质量影响的频响函数。本发明首先基于Sherman‑Morrison公式推导传感器质量影响消除公式,根据公式实测所需的频响函数,代入公式依次逐一消除,最终实现多传感器质量消除。本发明实现了消除模态试验接触式测量方法中传感器附加质量造成的测量误差,通过对实测频响函数信号的处理消除了多传感器质量对频响函数的不利影响,具有实际工程意义。
Description
技术领域
本发明涉及一种模态试验,具体涉及一种模态试验中多传感器附加质量消除方法。
背景技术
模态试验主要有接触式与非接触式两种测量方法。接触式测量需要在结构上布置传感器(如加速度传感器)来测量结构的动态响应;非接触式测量(如激光测振仪),与被测对象之间没有因为接触产生附加影响,但测量成本高,测量环境要求苛刻,通用性较差。
因此,接触式测量方法的应用更为广泛,但针对轻质柔性结构,传感器的布置不可避免地引起附加质量,导致结构特性改变,从而造成测量误差。特别是当对轻质柔性结构布置多个传感器时,附加质量的影响更为严重,使得结构的固有特性发生变化,导致模态试验结果的准确性降低。如何在多传感器布置的情况下,消除附加质量对实测频响函数的影响,得到准确的待测结构动态特性,已成为亟待解决的实际工程问题。
发明内容
发明目的:本发明的目的在于针对现有技术的不足,提供一种模态试验中多传感器附加质量消除方法。
技术方案:本发明提供了一种模态试验中多传感器附加质量消除方法,包括以下步骤:
(1)基于Sherman-Morrison公式推导模态试验中多传感器质量影响消除公式;
(2)确定模态试验中布置的传感器个数及消除传感器质量的顺序;
(3)基于实测的频响函数根据步骤(1)的消除公式逐一计算消除各传感器质量影响的频响函数。
进一步,步骤(1)包括以下步骤:
(1.1)模态试验布置的传感器个数为n,对应的各位置节点分别为a1,a2,…an,假定n个传感器的附加质量分别为在j点激励、i点测量的加速度频响函数记为其中右上标(a1,a2,…an)表示该频响函数包含a1,a2,…an点处加速度传感器附加质量影响,而对应原结构即不包含传感器附加质量影响的结构的加速度频响函数记为Aij;
(1.2)实测过程中结构系统动态刚度矩阵为Z,亦表示为右上标(a1,a2,…an)表示该动刚度矩阵包含a1,a2,…an节点处加速度传感器附加质量的影响,传感器质量矩阵ΔM,消除传感器质量的影响后系统动态刚度矩阵为Z*,则Z*与矩阵变换量ΔZ可分别表示为:
其中,ω为角频率;
(1.3)将公式(2)分解成n项之和:
其中,1的各下标表示其所处的节点位置;
为了简化公式(3),分别令
et=[0 … 1t … 0]T (16)
Jt=[0 … ω2mat … 0] (17)
其中,ΔZt表示在节点t处的质量修正量,et与Jt表示只在节点t处有非零项,分别为1与将公式(4)、(5)和(6)代入公式(3)得:
(1.3)将公式(7)代入公式(1),最终获得消除传感器质量影响的质量矩阵:
其中,由与组合可获得消除an处传感器质量的影响后的再与组合可获得消除an-1处传感器质量的影响后依次类推,最终消除所有传感器质量对动态刚度矩阵的影响;
(1.4)根据位移频响函数α与动刚度矩阵Z的关系,可知α=Z-1,若ΔZ表示成
ΔZ={u}{v}T (20)
其中{u}、{v}分别为行向量,代入基本的Sherman-Morrison公式得到频响函数矩阵为:
其中,α*、Z*为消除传感器质量影响以后的位移频响函数和动刚度矩阵;
(1.5)公式(8)中每一项皆可表示为公式(9)的形式,故每一步消除代入公式(10),再根据位移导纳与加速度导纳的关系A=-ω2α,其中A加速度频响函数矩阵,最终可以推导获得多传感器质量消除一般性公式:
其中,为消除节点an处传感器质量影响后的加速度频响函数,每一项为对应的加速度频响函数矩阵A中第i行j列的元素。
进一步,步骤(2)对于n个传感器,按照节点an、an-1、an-2、…、a2、a1的顺序依次消除。
进一步,步骤(3)包括以下步骤:
(3.1)按照式(11)首先消除节点an处传感器质量的影响,测量并获得所需加速度频响函数和计算得到消除节点an处传感器质量影响的加速度频响函数
(3.2)再对节点an-1处的传感器质量影响进行消除,根据式(11)得到消除公式:
其中,所需已通过步骤(3.1)获得,而和分别通过公式(11)计算获得;
(3.3)重复上述操作,依次进行消除。
有益效果:本发明首先基于Sherman-Morrison公式推导多传感器质量影响消除公式,根据公式实测所需的频响函数,代入公式依次逐一消除,最终实现多传感器质量消除。
本发明实现了消除模态试验接触式测量方法中传感器附加质量造成的测量误差,通过对实测频响函数信号的处理消除了多传感器质量对频响函数的不利影响,具有实际工程意义。
附图说明
图1是本发明实际例中三自由度弹簧-阻尼-质量系统示意图;
图2是消除处理的过程说明;
图3是理论值FRFA13与最终修正值对比。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
实施例:一种模态试验中多传感器附加质量消除方法,采用一个简单的三自由度弹簧-阻尼-质量系统来验证,系统的参数分别为:质量块质量为m1=m2=m3=0.1kg,弹簧刚度系数k1=k3=500N/m,k2=800N/m,阻尼系数为c1=0.4N·m/s,c2=c3=0.2N·m/s,加速度传感器质量ma1=ma2=ma2=0.02kg,如图1所示。
具体操作如下:
3个加速度传感器附加质量分别mi(i=1,2,…3),布置位置节点分别为1,2,3。在3点激励、1点测量的频响函数记为其中右上标(1,2,3)表示该频响函数包含1,2,3点处加速度传感器附加质量影响,而对应原结构频响函数记为A13,原结构为不包含传感器附加质量影响的结构。
确定依次消除传感器质量的顺序,依次按照节点3、2、1逐一消除:
确定对实测频响函数进行消除,先消除节点3传感器质量的影响,推导出消除公式:
根据公式(13)获取计算所需的消除节点3处加速度传感器质量影响的频响函数和将需要实测的频响函数和代入到公式(13),即可获得消除3号节点处传感器质量影响修正值
再对节点2处的传感器质量影响进行消除,推导出消除公式:
为了获取频响函数计算所需和其中已通过式(13)计算获得,而和项可带入一般性公式计算获得,分析过程见图2。
重复上述操作,并再消除最后1节点处的加速度传感器质量的影响:
为了获取频响函数计算所需和其中已通过公式(14)处理获得,而可重复上述步骤获得,分析过程见图2,最终实现多传感器质量消除目的。
由图3可知,横坐标f为频率,纵坐标FRF为频响函数,分别给出了有传感器影响的频响函数及无传感器影响的频响函数,将无传感器质量影响的频响函数作为理论值A13,有传感器质量影响的作为实测值对实测值采用本专利消除方法获得修正后的曲线与理论值吻合一致,最终实现多传感器质量消除目的。
Claims (4)
1.一种模态试验中多传感器附加质量消除方法,其特征在于:包括以下步骤:
(1)基于Sherman-Morrison公式推导模态试验中多传感器质量影响消除公式;
(2)确定模态试验中布置的传感器个数及消除传感器质量的顺序;
(3)基于实测的频响函数根据步骤(1)的消除公式逐一计算消除各传感器质量影响的频响函数。
2.根据权利要求1所述的模态试验中多传感器附加质量消除方法,其特征在于:步骤(1)包括以下步骤:
(1.1)模态试验布置的传感器个数为n,对应的各位置节点分别为a1,a2,…an,假定n个传感器的附加质量分别为(i=1,2,…n),在j点激励、i点测量的加速度频响函数记为其中右上标(a1,a2,…an)表示该频响函数包含a1,a2,…an点处加速度传感器附加质量影响,而对应原结构即不包含传感器附加质量影响的结构的加速度频响函数记为Aij;
(1.2)实测过程中结构系统动态刚度矩阵为Z,亦表示为右上标(a1,a2,…an)表示该动刚度矩阵包含a1,a2,…an节点处加速度传感器附加质量的影响,传感器质量矩阵ΔM,消除传感器质量的影响后系统动态刚度矩阵为Z*,则Z*与矩阵变换量ΔZ可分别表示为:
<mrow>
<msup>
<mi>Z</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mi>Z</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>Z</mi>
<mo>=</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>Z</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ω为角频率;
(1.3)将公式(2)分解成n项之和:
<mrow>
<mi>&Delta;</mi>
<mi>Z</mi>
<mo>=</mo>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mn>1</mn>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>m</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mn>1</mn>
<mrow>
<mi>a</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>m</mi>
<mrow>
<mi>a</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,1的各下标表示其所处的节点位置;
为了简化公式(3),分别令
<mrow>
<msup>
<mi>&Delta;Z</mi>
<mi>t</mi>
</msup>
<mo>=</mo>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mn>1</mn>
<mi>t</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>m</mi>
<msub>
<mi>a</mi>
<mi>t</mi>
</msub>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
et=[0 … 1t … 0]T (5)
Jt=[0 … ω2mat … 0] (6)
其中,ΔZt表示在节点t处的质量修正量,et与Jt表示只在节点t处有非零项,分别为1与将公式(4)、(5)和(6)代入公式(3)得:
<mrow>
<mi>&Delta;</mi>
<mi>Z</mi>
<mo>=</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
</msup>
<mo>+</mo>
<mn>...</mn>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
(1.3)将公式(7)代入公式(1),最终获得消除传感器质量影响的质量矩阵:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>Z</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mi>Z</mi>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</msup>
<mo>+</mo>
<mn>...</mn>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</msup>
<mo>+</mo>
<mn>...</mn>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</msup>
<mo>+</mo>
<mn>...</mn>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</msup>
<mo>+</mo>
<mn>...</mn>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</msup>
<mo>+</mo>
<mn>...</mn>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mn>...</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>&Delta;Z</mi>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,由与组合可获得消除an处传感器质量的影响后的再与组合可获得消除an-1处传感器质量的影响后依次类推,最终消除所有传感器质量对动态刚度矩阵的影响;
(1.4)根据位移频响函数α与动刚度矩阵Z的关系,可知α=Z-1,若ΔZ表示成
ΔZ={u}{v}T (9)
其中{u}、{v}分别为行向量,代入基本的Sherman-Morrison公式得到频响函数矩阵为:
<mrow>
<msup>
<mi>&alpha;</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>*</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mi>&alpha;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>{</mo>
<mi>u</mi>
<mo>}</mo>
<mo>)</mo>
<mo>(</mo>
<msup>
<mrow>
<mo>{</mo>
<mi>v</mi>
<mo>}</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mrow>
<mo>{</mo>
<mi>v</mi>
<mo>}</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>&alpha;</mi>
<mo>{</mo>
<mi>u</mi>
<mo>}</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,α*、Z*为消除传感器质量影响以后的位移频响函数和动刚度矩阵;
(1.5)公式(8)中每一项皆可表示为公式(9)的形式,故每一步消除代入公式(10),再根据位移导纳与加速度导纳的关系A=-ω2α,其中A加速度频响函数矩阵,最终可以推导获得多传感器质量消除一般性公式:
<mrow>
<msubsup>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</msub>
<msubsup>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>A</mi>
<mrow>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>m</mi>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</msub>
<msubsup>
<mi>A</mi>
<mrow>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为消除节点an处传感器质量影响后的加速度频响函数,每一项为对应的加速度频响函数矩阵A中第i行j列的元素。
3.根据权利要求2所述的模态试验中多传感器附加质量消除方法,其特征在于:步骤(2)对于n个传感器,按照节点an、an-1、an-2、…、a2、a1的顺序依次消除。
4.根据权利要求3所述的模态试验中多传感器附加质量消除方法,其特征在于:步骤(3)包括以下步骤:
(3.1)按照式(11)首先消除节点an处传感器质量的影响,测量并获得所需加速度频响函数和计算得到消除节点an处传感器质量影响的加速度频响函数
(3.2)再对节点an-1处的传感器质量影响进行消除,根据式(11)得到消除公式:
<mrow>
<msubsup>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>n</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</msub>
<msubsup>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>A</mi>
<mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>m</mi>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</msub>
<msubsup>
<mi>A</mi>
<mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
2
其中,所需已通过步骤(3.1)获得,而和分别通过公式(11)计算获得;
(3.3)重复上述操作,依次进行消除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710598555.4A CN107356387B (zh) | 2017-07-21 | 2017-07-21 | 一种模态试验中多传感器附加质量消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710598555.4A CN107356387B (zh) | 2017-07-21 | 2017-07-21 | 一种模态试验中多传感器附加质量消除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107356387A true CN107356387A (zh) | 2017-11-17 |
CN107356387B CN107356387B (zh) | 2018-06-19 |
Family
ID=60286002
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710598555.4A Active CN107356387B (zh) | 2017-07-21 | 2017-07-21 | 一种模态试验中多传感器附加质量消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107356387B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108038343A (zh) * | 2018-01-31 | 2018-05-15 | 东南大学 | 一种保证结构特定频率不变的弹簧修正方法 |
CN108256256A (zh) * | 2018-01-31 | 2018-07-06 | 东南大学 | 一种保证阻尼系统特定模态频率不变的方法 |
CN108984480A (zh) * | 2018-06-13 | 2018-12-11 | 东南大学 | 一种多个加速度传感器附加质量影响的一步消除方法 |
CN109299512A (zh) * | 2018-08-27 | 2019-02-01 | 东南大学 | 一种基于质量影响的快速灵敏度分析方法 |
CN109344433A (zh) * | 2018-08-27 | 2019-02-15 | 东南大学 | 基于响应信号的灵敏度数值计算方法 |
CN109506696A (zh) * | 2018-12-17 | 2019-03-22 | 湖北工业大学 | 一种基于传感器质量辨识的频响函数校准方法 |
CN110569568A (zh) * | 2019-08-21 | 2019-12-13 | 湖北工业大学 | 自由模态测试中消除悬挂刚度和传感器质量影响的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106052743A (zh) * | 2016-05-17 | 2016-10-26 | 湖北工业大学 | 一种评估传感器质量对频响函数影响大小的方法 |
CN106871933A (zh) * | 2017-03-03 | 2017-06-20 | 湖北工业大学 | 一种消除频响函数中多传感器附加质量的方法 |
US20170177769A1 (en) * | 2015-12-22 | 2017-06-22 | Dassault Systemes Simulia Corp. | Effectively Solving Structural Dynamics Problems With Modal Damping In Physical Coordinates |
-
2017
- 2017-07-21 CN CN201710598555.4A patent/CN107356387B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170177769A1 (en) * | 2015-12-22 | 2017-06-22 | Dassault Systemes Simulia Corp. | Effectively Solving Structural Dynamics Problems With Modal Damping In Physical Coordinates |
CN106052743A (zh) * | 2016-05-17 | 2016-10-26 | 湖北工业大学 | 一种评估传感器质量对频响函数影响大小的方法 |
CN106871933A (zh) * | 2017-03-03 | 2017-06-20 | 湖北工业大学 | 一种消除频响函数中多传感器附加质量的方法 |
Non-Patent Citations (3)
Title |
---|
任军 等: "模态测试中力传感器附加质量辨识及消除方法研究", 《振动与冲击》 * |
尹涛: "时域响应中传感器附加质量影响消除", 《振动与冲击》 * |
董自卫 等: "消除传感器附加质量影响的机枪动力特性测试技术", 《自动测量与控制》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108038343A (zh) * | 2018-01-31 | 2018-05-15 | 东南大学 | 一种保证结构特定频率不变的弹簧修正方法 |
CN108256256A (zh) * | 2018-01-31 | 2018-07-06 | 东南大学 | 一种保证阻尼系统特定模态频率不变的方法 |
CN108984480A (zh) * | 2018-06-13 | 2018-12-11 | 东南大学 | 一种多个加速度传感器附加质量影响的一步消除方法 |
CN109299512A (zh) * | 2018-08-27 | 2019-02-01 | 东南大学 | 一种基于质量影响的快速灵敏度分析方法 |
CN109344433A (zh) * | 2018-08-27 | 2019-02-15 | 东南大学 | 基于响应信号的灵敏度数值计算方法 |
CN109299512B (zh) * | 2018-08-27 | 2019-06-18 | 东南大学 | 一种基于质量影响的快速灵敏度分析方法 |
CN109506696A (zh) * | 2018-12-17 | 2019-03-22 | 湖北工业大学 | 一种基于传感器质量辨识的频响函数校准方法 |
CN110569568A (zh) * | 2019-08-21 | 2019-12-13 | 湖北工业大学 | 自由模态测试中消除悬挂刚度和传感器质量影响的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107356387B (zh) | 2018-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107356387B (zh) | 一种模态试验中多传感器附加质量消除方法 | |
CN108592952A (zh) | 基于杆臂补偿与正反倍速率同时标定多mimu误差的方法 | |
CN103808349B (zh) | 矢量传感器的误差校正方法和装置 | |
CN108387951B (zh) | 一种基于重复线校正海空重力仪格值的新算法 | |
CN110879302B (zh) | 一种石英谐振差动式加速度计温度补偿方法 | |
CN106871933B (zh) | 一种消除频响函数中多传感器附加质量的方法 | |
CN104020036B (zh) | 一种确定横向集中载荷下环形预应力薄膜最大挠度的方法 | |
CN102889076A (zh) | 陀螺测斜仪标定方法 | |
EP3714240B1 (en) | Weighting method and storage medium thereof | |
JP3313028B2 (ja) | 張力のかかっているケーブルの曲げ剛性及び張力の測定方法 | |
CN112182697A (zh) | 一种有阻尼吊杆系统张力的高精度动测法 | |
CN104697712B (zh) | 一种回转体工件质心检验方法 | |
CN106383003A (zh) | 基于柔度识别的索结构索力的测量方法及测量系统 | |
CN110619173B (zh) | 一种利用曲率指标辨识一维杆结构刚度薄弱环节的方法 | |
CN114577318B (zh) | 一种车载称重模块及其传感方法 | |
CN108333061A (zh) | 一种测量应力松弛的系统及测量方法 | |
CN103823083B (zh) | 提高加速度计校准精度的方法与系统 | |
CN116187053A (zh) | 基于蒙特卡洛法的臂架应力测量不确定度方法及装置 | |
CN110487392A (zh) | 一种设备弹性安装时基座导纳测量的修正方法 | |
CN105842481A (zh) | 一种惯性导航系统中基于十字叉的三轴加速度计标定方法 | |
CN106289492B (zh) | 一种砝码量值分量组合检定法 | |
RU2562273C2 (ru) | Стенд для измерения массо-инерционных характеристик изделия | |
KR101326280B1 (ko) | 소형 무게중심 측정장치 및 방법 | |
CN111649903A (zh) | 一种利用动模型实验间接测量气动阻力系数的方法 | |
CN109883611A (zh) | 一种简易的力传感器标定装置和方法 |
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 |