一种微振动信号处理方法
【专利摘要】一种微振动信号处理方法，对采集的微振动信号进行处理得到扰振力和误差分析。首先将标定盘固装在六分量八传感器扰振力测试台上，用标准力锤敲击标定盘上16个标定点，每次敲击时，采集力锤和8传感器电压时域信号以及六分量扰振力测试台频响函数。通过标定过程得到标定矩阵和标定误差。其次，在六分量八传感器扰振力测试台上的活动部件对测试台有微小扰动时，通过数采系统采集时域数据，通过快速傅里叶变换得到频域结果，通过标定矩阵得到扰振力频域信号，通过反傅里叶变换得到扰振力时域数据。
【专利说明】一种微振动信号处理方法
【技术领域】
[0001]本发明涉及一种微振动信号处理方法，可用于对航天器微小扰动载荷在8个传感器上的振动信号进行动态测量，通过相关理论对微振动信号处理得到扰振力。
【背景技术】
[0002]目前的航天器大多都属于大型柔性展开式机构，且带有大量的光学元件，它们对指向精度和稳定度均提出了很高的要求。另外，在现代航天器姿态控制系统中，反作用轮、单框架力矩陀螺和太阳翼驱动机构等是其控制系统中的重要元件，它们在提供必要的控制动力的同时，也会引起一些有害振动(为简单起见，下面将上述三种系统统称为扰动源)。这些扰动主要由飞轮不平衡、轴承扰动、电机扰动、电机驱动误差等引起的，其中飞轮不平衡是导致飞轮振动的最主要原因，这些扰动力和扰动力矩会降低体太空中精密性仪器的性能指标，因此测量和分析航天器有效载荷扰动的动态特性，对于分析并消除扰动从而提高航天器的姿态控制精度和加强航天器的安全设计有着非常重要的工程意义。
[0003]由于航天器扰动源的扰动很小，个别有效载荷如动量轮在空间三个方向只能产生几十毫牛顿甚至几毫牛顿的微弱扰动，要想在具有相对强烈干扰背景噪音的地面实验室中测量此类扰动十分困难，而其对应传感器的精度要求非常高。
[0004]目前，国内外尚未见有关此类微小振动测量系统微振动信号处理的文献报导。
【发明内容】
[0005]本发明要解决的技术问题是:克服现有技术的不足，提供一种微振动信号处理方法，通过标定得到试验测试系统的频响函数，有频响函数得到试验测试系统的标定矩阵，在航天器活动部件对该测试系统有微小扰动时，通过标定矩阵得到等效到试验测试台某一点的扰振力。
[0006]本发明要解决其技术问题所采用的技术方案是:
[0007]—种微振动信号处理方法，所述的电压信号测量是在八传感器六分量刚性测力台上进行，步骤如下:
[0008](I)将标定盘安装在八传感器六分量刚性测力台上，用螺栓将其拧紧；标定盘是一个正方形铝盘，所述标定盘四个侧面上各自均布有3个敲击点，上表面上布有4个敲击点，下表面与八传感器六分量刚性测力台连接；
[0009](2)用力锤依次敲击标定盘上的16个敲击点，并采集力锤的时域信号、敲击每个点时测力台的输出电压信号以及频响函数；
[0010](3)计算锤击力O分量误差、非O分量误差以及标定矩阵；
[0011](4)拆下标定盘；
[0012](5)将微振动信号源安装在八传感器六分量刚性台上，用螺栓拧紧；
[0013](6)扰动源通电开始工作；
[0014](7)通过八传感器六分量测力台采集扰动源工作时传感器的电压时域信号；[0015](8)将电压时域信号进行傅里叶变换得到电压的频域信号；
[0016](9)将电压频域信号与标定矩阵相乘得到扰动源扰动力的频域信号；
[0017](10)将扰动力的频域信号反傅里叶变换可以得到扰动力的时域信号，即完成对微振动信号的处理。
[0018]所述八传感器六分量刚性测力台包括八个单向压电传感器，其中四个压电传感器沿水平方向布置，保证能够测试沿水平方向的两个扰动力和绕垂直方向的扰动力矩；另外四个压电传感器沿垂直方向布置，保证能够测试沿垂直方向的扰动力和绕水平方向的两个扰动力矩。
[0019]步骤(3)中所述计算锤击力O分量误差、非O分量误差以及标定矩阵具体为:
[0020]通过公式e;。=月。/<计算锤击力O分量误差<，通过公式e; = I巧~F；/F；计算非
O分量误差^行定矩阵为W6X8(?)=C6XnF' —(ωΚΥωΚ?^ωΚΥω)]-1，其中，各锤击点对中心标定点的实测等效载荷分量包含O项和非O项，分别用Ff，Fi1表示，对应的反算值用月。，月1表示，i=l~6,
[0021]矩阵下标中的η表示试验中加载的次数；C表示实际加载的载荷矩阵与等效载荷矩阵之间的转换矩阵，维数为6Xn;F'表示实际加载的载荷，是一个对角阵，矩阵中的F' π等于第ii次加载的载荷值；T为八个力传感器的响应信号，维数为8Χη。
[0022]本发明与现有技术相比具有以下优点:
[0023]( I)本发明测量误差检 验采用瞬态检验方法。稳态检验使用激振器采用扫频方式，激振器上有标准力的传感器，振动测试台的输出力和传感器测得的力比较从而得到测量误差。稳态检验对所关心的频率区间检验精度较高，缺点是稳态检验时间长，并且由于激振器体积和质量比较大，安装难度高从而引入了检验力之外的力的分量，工程应用中困难很大，而瞬态检验工程实现简单，精度满足要求。
[0024](2)本发明应用模态分析理论中的载荷识别技术即频响函数矩阵求逆法来得到电压信号与力信号之间的标定系数，实现整体标定，从而求出扰动载荷。这是一种整体标定的方法，这就减少了由于多个传感器在组成整体使用时其标定误差的累积问题；同时，每次标定时无需重新拆装测试系统，很大程度上降低了工作量，使用方便。
【专利附图】
【附图说明】
[0025]图1为本发明标定点位置示意图；
[0026]图2为标定力锤示意图；
[0027]图3为本发明方法流程图。
【具体实施方式】
[0028]如图3所示，本发明提供了一种微振动信号处理方法，所述的电压信号测量是在八传感器六分量刚性测力台上进行，其特征在于步骤如下:
[0029](I)将标定盘安装在八传感器六分量刚性测力台上，用螺栓将其拧紧；标定盘是一个正方形铝盘，如图1所示，标定盘四个侧面上各自均布有3个敲击点，上表面上布有4个敲击点，下表面与八传感器六分量刚性测力台连接；八传感器六分量刚性测力台是由八个单向压电传感器组成，其中四个压电传感器沿水平方向布置，保证能够测试沿水平方向的两个扰动力和绕垂直方向的扰动力矩；另外四个压电传感器沿垂直方向布置，保证能够测试沿垂直方向的扰动力和绕水平方向的两个扰动力矩。
[0030](2)用力锤依次敲击标定盘上的16个敲击点，并采集力锤的时域信号、敲击每个点时测力台的输出电压信号以及频响函数；力锤的示意图如图2所示，
[0031](3)计算锤击力O分量误差、非O分量误差以及标定矩阵；
[0032](4)拆下标定盘；
[0033](5)将微振动信号源安装在八传感器六分量刚性台上，用螺栓拧紧；
[0034](6)扰动源通电开始工作；
[0035](7)通过八传感器六分量测力台采集扰动源工作时传感器的电压时域信号；
[0036](8)将电压时域信号进行傅里叶变换得到电压的频域信号；
[0037](9)将电压频域信号与标定矩阵相乘得到扰动源扰动力的频域信号；
[0038](10)将扰动力的频域信号反傅里叶变换可以得到扰动力的时域信号，即完成对微振动信号的处理。
[0039]首先应该求得相应的标定矩阵W(co)，整体标定矩阵实际上是八个压电力传感器的响应信号与作用在等效中心点的载荷之间的频响函数矩阵。在本方法中，等效中心为八分量传感器装置的负载盘上表 面的中心点，在此中心点直接施加三个平动力和三个力矩有困难，本文在假设负载盘为刚体的基础上，安装一个刚度很高的标定台，并选择如图所示的16个加载点。根据空间力系平衡理论可知，利用选择的16个加载点能等效出作用在负载盘几何中心的三个平动力和三个力矩；图1中O点为负载盘上表面的中心点，表1列出了图中各参数的值。
[0040]表1标定点位置参数
[0041]
【权利要求】
1.一种微振动信号处理方法，所述的电压信号测量是在八传感器六分量刚性测力台上进行，其特征在于步骤如下: (1)将标定盘安装在八传感器六分量刚性测力台上，用螺栓将其拧紧；标定盘是一个正方形铝盘，所述标定盘四个侧面上各自均布有3个敲击点，上表面上布有4个敲击点，下表面与八传感器六分量刚性测力台连接； (2)用力锤依次敲击标定盘上的16个敲击点，并采集力锤的时域信号、敲击每个点时测力台的输出电压信号以及频响函数； (3)计算锤击力O分量误差、非O分量误差以及标定矩阵； (4)拆下标定盘； (5)将微振动信号源安装在八传感器六分量刚性台上，用螺栓拧紧； (6)扰动源通电开始工作； (7)通过八传感器六分量测力台采集扰动源工作时传感器的电压时域信号； (8)将电压时域信号进行傅里叶变换得到电压的频域信号； (9)将电压频域信号与标 定矩阵相乘得到扰动源扰动力的频域信号； (10)将扰动力的频域信号反傅里叶变换可以得到扰动力的时域信号，即完成对微振动信号的处理。
2.根据权利要求1所述的一种微振动信号处理方法，其特征在于:所述八传感器六分量刚性测力台包括八个单向压电传感器，其中四个压电传感器沿水平方向布置，保证能够测试沿水平方向的两个扰动力和绕垂直方向的扰动力矩；另外四个压电传感器沿垂直方向布置，保证能够测试沿垂直方向的扰动力和绕水平方向的两个扰动力矩。
3.根据权利要求1所述的一种微振动信号处理方法，其特征在于:步骤(3)中所述计算锤击力O分量误差、非O分量误差以及标定矩阵具体为: 通过公式< =Fi01F11计算锤击力O分量误差4，通过公式e; -1竚-Fl\lFl计算非O分量误差 < 标定矩阵为W6X8(?)=C6XnF'.(ωΚΥωΚ?^οΟτΥω)]-1，其中，各锤击点对中心标定点的实测等效载荷分量包含O项和非O项，分别用/^，<表示，对应的反算值用F F10 ,月1表示，i=l~6, 矩阵下标中的η表示试验中加载的次数；C表示实际加载的载荷矩阵与等效载荷矩阵之间的转换矩阵，维数为6Xn;F'表示实际加载的载荷，是一个对角阵，矩阵中的€等于第ii次加载的载荷值；T为八个力传感器的响应信号，维数为8Χη。
【文档编号】G01H11/06GK103499387SQ201310446695
【公开日】2014年1月8日 申请日期:2013年9月26日 优先权日:2013年9月26日
【发明者】杨文涛, 程伟, 李果, 董家俊, 王成伦, 王云峰, 刘凤晶, 赵煜, 王光远, 沈中 申请人:北京空间飞行器总体设计部, 北京航空航天大学