1、ICS 17160J 04 a亘中华人民共和国国家标准GBT 204851 3-2007IS0 1 6063-1 3：200 1振动与冲击传感器校准方法第1 3部分：激光干涉法冲击绝对校准Methods for the calibration of vibration and shock transducers-Part 1 3：Primary shock calibration by laser interferometry2007-07-02发布(IS0 1606313：2001，IDT)2007-12-01实施宰瞀髁鬻瓣警矬赞星发布中国国家标准化管理委员会及111目 次前言引言-1范围“
2、2规范性引用文件”3测量不确定度”4仪器设备要求-41 总则-42基于砧体(刚体)运动的冲击机-43基于细长棒中波传播的冲击机“44冲击机和激光干涉仪的隔振块”45激光器”46干涉仪一47示波器48带计算机接口的波形记录仪49带数据处理程序的计算机410滤波器411其他要求5环境条件6优选的加速度及脉冲持续时间7方法71测量步骤72数据采集73数据处理8校准结果报告附录A(规范性附录)校准中测量不确定度的表述附录B(资料性附录)程序介绍附录C(资料性附录)复灵敏度幅值和相移计算的替代方法参考文献GBT 2048513-2007150 16063-13：200111l11124445555566
3、66669m坨坫前 言GBT 2048513-2007ISO 16063-13：2001GBT 20485振动与冲击传感器的校准方法主要由以下几大类组成：第1类：基本概念；第2类：绝对法校准(如：激光干涉法振动绝对校准、激光干涉法冲击绝对校准等)第3类：比较法校准(如：振动比较法校准、冲击比较法校准等)，第4类：环境模拟校准(如：磁灵敏度、声灵敏度、基座应变灵敏度、横向振动灵敏度等)；本部分是GBT 20485振动与冲击传感器的校准方法第2类中的一个部分。本部分等同采用ISO 1606313：2001振动与冲击传感器校准方法第13部分：激光干涉法冲击绝对校准(英文版)。本部分等同翻译ISO 1
4、606313：2001。为使用方便，本部分作了如下编辑性修改：“国际标准的本部分”一词改为“本部分”。删除国际标准的前言。用小数点“”代替作为小数点的逗号“，”。对IsO 1606313：2001中引用的其他国际标准，有被采用为我国标准的，用我国标准代替对应的国际标准，未被采用为我国标准的直接引用国际标准。本部分的附录A为规范性附录，附录B、附录C为资料性附录。本部分由中国计量科学研究院提出。本部分由全国机械振动与冲击标准化技术委员会归口。本部分起草单位：中国航空工业第一集团公司北京长城计量测试技术研究所、中国计量科学研究院。本部分主要起草人：李新良、曾吾、于梅、孙桥。本部分为首次制定。GBT
5、 2048513-2007IS0 16063-13：2001引 言按照定义，冲击灵敏度s“由加速度计输出峰值和输入加速度峰值关系来决定。s“不是一个唯一的量，是随着冲击脉冲持续时间、冲击波形、被校传感器的带宽、配套信号适调仪的频响而变化。在频域计算某一个频率点对应的复灵敏度可以确定加速度计的线性。本部分提出了计算复灵敏度幅值s。和相移吼的数据处理程序，再计算冲击灵敏度s。“(见附录c)。本部分给出了运动的时间历程的绝对法测量运动的方法，从根本上不同于GBT 204851介绍的速度改变冲击校准法。因此，得到的冲击灵敏度不同于后者得到的冲击校准因子，但和I$O 53474(将转换为ISO 1605
6、322)所得到的校准因子一致。1范围GBT 2048513-2007IS0 16063-13：2001振动与冲击传感器校准方法第13部分：激光干涉法冲击绝对校准本部分规定了用激光干涉法测量冲击过程时间一位移信号对直线加速度计进行绝对法冲击校准的仪器设备和操作程序。该方法适用的冲击脉冲持续时间范围为005 mslO ms，峰值范围为102 n1s2105 rns2(与脉冲持续时间有关)，给出的是冲击灵敏度。2规范性引用文件下列文件中的条款通过GBT 20485的本部分的引用而成为本部分的条款。凡是注日期的引用文件，其随后所有的修改单(不包括勘误的内容)或修订版均不适用于本部分。然而，鼓励根据本部
7、分达成协议的各方研究是否可使用这些文件的最新版本。凡是不注日期的引用文件，其最新版本适用于本部分。GBT 204851振动与冲击传感器校准方法第1部分：基本概念(Iso 160631：1998，IDT)GBT 204851l一2006 振动与冲击传感器校准方法 第1l部分：激光干涉法振动绝对校准(IS0 1606311：1999，IDT)ISO 534722：1997振动与冲击传感器校准方法第22部分：加速度计共振测试通用方法3测量不确定度冲击灵敏度的测量不确定度：在设定参考加速度峰值i03 ms2，参考冲击脉冲持续时间2 ms和放大器参考增益时，为读数的l；对所有冲击峰值、脉冲持续时间，不超
8、过读数的2。上述规定适用于高精度等级传感器(例如，参考标准加速度计)的校准，这样的校准应仔细操作以确保给出的全部不确定度分量应足够小到符合相应的规定(不确定度的分量参看附录A)。特别是由传感器或冲击机的固有模态激发的频谱能量应比校准频率范围内的频谱能量要小得多，传感器固有频率的测试按照ISO 534722：1997进行。通常应避免使用第1章、第6章中给出的相对短持续时间的脉冲。本部分的所有使用者应按照附录A编制实际的不确定度分量表。注：根据GBT 204851，测量不确定度表示为测量的扩展不确定度(简称不确定度)。4仪器设备要求41 总则为了满足第1章中的测量范围及达到第3章中的不确定度要求，
9、本章给出了所需仪器设备的推荐技术指标。42基于砧体(刚体)运动的冲击机该冲击机是用一个能自由移动的锤(弹体)来撞击安装有加速度计的砧体(目标)。锤撞击砧体后，砧体能够自由直线加速，同时锤被自动锁住。为获得期望的脉冲宽度和脉冲波形，在锤和砧之间应安装钢弹簧或缓冲垫(由橡胶、纸等某种可实现脉冲波形的材料制成)，所得冲击脉冲应是一个近似半正弦、半正弦平方或高斯(Gaussian)的加速度波形。锤和砧的共振频率至少为IOT，T为脉冲持续时间。】GBT 2048513-2007ISO 1606313：2001为了避免冲击机结构共振的影响，锤和砧应尽量与结构隔离。锤和砧的安装中心线间距最大值为o2 mm。
10、砧应由对称力支撑，而非对称力支撑会导致其旋转、偏离于直线运动。安装加速度计的表面的粗糙度，用算术平均偏差Ra表示，其值小于l,um。表面平面度的要求是：在对应于任何被校准传感器的最大安装表面内，该表面应包容在距离为5,um两个平行平面之间。联接传感器用的螺纹孔与安装表面的垂直偏差应小于10 pm，即该孔的轴线允许在一个直径为10 pm的柱体区域内，并且柱面的高度与孔的深度相同。注1：当砧(或砧和锤)安装在空气轴承上时，应满足上述要求(见图1和参考文献1)。图1所示的冲击机可以产生半正弦平方加速度波形。注2。按照ISO 53474t1993进行比较冲击校准的传统冲击机(见参考文献2和3)，不可能
11、产生能用激光干涉法精确测量的运动。l冲击机(42)2弹簧单元l3空气轴承冲击锤(钢，直径30 minx长200 ram)；4垫f5空气轴承砧体(钢直径30 rnm长200 turn)I6加速度计7放大器，8数字波形记录仪(48)l9激光器(45)；10干涉仪(46)11光电接收器(46)12-隔振块1(44)I13隔振块2(44)。图1 基于砧体(刚体)运动的冲击校准测量系统实例(加速度峰值范围102 ms,z一5103 ms2)43基于细长棒中波传播的冲击机该冲击机主要由一个可运动部件构成，例如，被加速的钢球(弹体)撞击和棒连接的缓冲件(如另一个相同直径的钢球)，而棒的另一端面安装了加速度计
12、。这里棒以柔性支撑来避免冲击机结构共振的2GBT 2048513-2007ISO 16063-13：2001影响。钢球和棒应尽量安装在同一直线上以满足第3章中的不确定度要求。为达到第3章中所述不确定度，在测量期间，即数据采集的最重要时间段(最大值：l ms)，应尽量保证由加速度计安装表面引起的相对于直线运动的偏差最小。冲击机应具备数据采集的触发功能。安装加速度计的表面的粗糙度，用偏差的算术平均Ra表示，小于1 Ixm。表面平面度的要求是：在对应于任何被校准传感器的最大安装表面内，该表面应包容在距离为5 pm两个平行平面之间。连接传感器用的螺纹孔与安装表面的垂直偏差应小于10Lm，即该孔的轴线允
13、许在一个直径为10 pm的柱体区域内，并且柱面的高度与孔的深度相同。棒的尺寸设计(见参考文献4、Is)应考虑以下两个因素：校准单端安装加速度计时，激光光束必须能够射到棒的端面；保证足够的时间长度。最大冲击持续时间和可用数据采集测量时间是从主脉冲开始时刻到该脉冲在安装表面发生反射的时刻之间的时间长度(例如，如图2所示，棒长2 m时，时间间隔08 ms)。1冲击机2阀门(压缩空气)f3空气腔I4碰撞球对(tr径50 mm)，5硅橡胶；6铝管l7o形环，8棒(钛金属，直径25 nlm长2 000 mm)；9加速度计l10-放大器1l数字波形记录仪(48)I12应变片I13桥路放大器；14触发单元，1
14、5干涉仪(46)16激光器(45)t17光电接收器(46)。圈2基于细长棒中应力波传播的冲击校准测量系统实例(加速度峰值范圈103 ms2105 ms2)GBT 2048513-2007IS0 16063-13：2001图2是一个基于应力波在细长棒中传播的冲击机的实例。为了获得触发信号，在棒的柱面对称位置上贴两个应变片。带有两个钢球的冲击激励装置(图2所示)可产生所需加速度波形，用高斯函数的微分(即高斯速度脉冲)来描述(见参考文献6)。这种装置的优点是：在多次冲击校准中，具有良好的重复性；在不同加速度峰值上冲击频谱的频谱分量变化相对小(见参考文献13)。为满足不同校准条件，也可使用与图2尺寸不
15、同的棒。通常，棒内轴向位移是径向位置和频率的复杂函数，而频率取决于棒的材料特性和直径。这样，对于被校传感器可能会带来和频率有关的基座应变，或增加校准的不确定度，或二者都有。44冲击机和激光干涉仪的隔振块为了避免大地脉动引起的相对运动，以及冲击激励源支撑结构的反冲力对校准结果产生过大的影响，冲击机和干涉仪应安装在起隔振作用的同一个或两个不同的大质量块上。45激光器应使用发红光的氨氖激光器。在实验室条件下(即：大气压为100 kPa，温度23，相对湿度50)，激光波长为0632 81 pm。如果激光器具有手动或自动气压补偿，应将其设置为零或关闭。另外，还可以选择其他已知稳定波长的单频激光器。46干
16、涉仪应选用改进型迈克尔逊(Michelson)激光干涉仪，正交信号输出，带有两只频响满足所需带宽的光电接收器。所需带宽根据最大速度。，由下式得到：，。，一口。；316106 In-1 (1)改进的迈克尔逊干涉仪的结构见图3。四分之一波片将人射的线偏振光分成偏振方向相互垂直、相位差为90。的两个测量光束。在与线偏振参考光束发生干涉后，两个正交偏振光束由光学器件例如：沃尔斯顿(Wollaston)棱镜或偏振分光镜分开，并由两个光电二极管检测。1激光器，2偏振片，3分光镜4参考镜I5模拟质量块(测量反射镜)6加速度计7袄尔斯顿棱镜，8光电接收器。圈3正交输出的激光干涉仪GBT 20485131200
17、7ISO 1606313：2001改进的迈克尔逊干涉仪的两路输出幅值的偏移量应小于-L5，相对幅值偏差应小于士5，与90。名义角度的偏差应小于士5。为确保这些允差。要采用适当的方法调整两路干涉信号的偏移量、幅值大小及相角。干涉仪的测量光束和砧、棒或砧棒轴在一条直线上，以满足第3章中的测量不确定度的要求。如果传感器是单端形式的，砧、棒或砧棒端面应该抛光，以满足测量光束反射需要如果是背靠背传感器，其顶面也应该抛光，避免使用反射镜。在大加速度峰值测量时，需要带宽很宽。例如，在测量冲击峰值为105 ms2、脉冲持续时间为200zs的冲击过程(加速度波形见43)，干涉仪输出信号频率最大值为32 MHz，
18、详见参考文献6。也可以用其他合适的双光束干涉仪代替(改进)迈克尔逊干涉仪，例如，(改进)马赫一泽德(MathZehnder)干涉仪。47示波器如果48不包括示渡器功能，应该有一个DC50 MHz或上限频率更高的示渡器检查干涉仪和加速度计的输出信号波形。48带计算机接口的波形记录仪采用带计算机接口的波形记录仪，该记录仪具有模一数转换功能，并能存储干涉仪的两路正交输出信号及加速度计的输出信号。其幅值分辨力、采样率及存储容量应足够大，在校准的幅值范围内，满足第3章中规定的不确定度。通常，用于加速度计输出的幅值分辨力应不小于lo比特，用于干涉仪正交输出信号的分辨力应不小于8比特。加速度输出测量的典型幅
19、值分辨力应该不小于10比特，用一个两通道的波形记录仪记录干涉仪的输出信号，同时另一个波形记录仪(具有更高的分辨力和更低的采样率)记录加速度计的输出信号。示例；在加速度峰值为2 500 ms。、脉冲持续时间为2 ms的条件下校准加速度计采样频率为50 MHz(存储容量应1 M字节，才能满足3个通道存储需要)或更高。49带数据处理程序的计算机应使用带有数据处理程序的计算机，其程序应按照73中描述的计算过程编制。410滤波器加速度计和干涉仪输出信号均用模拟滤波器滤波。为了避免混叠和(或)抑制噪声，该滤波应具有适合的幅频和相频响应，以满足测量不确定度要求(见第3章)。数据处理程序中的数字滤波器也应满足
20、上述要求，见73。对干涉仪输出信号进行滤波时。应考虑参考文献63介绍的误差。4”其他要求在校准中为获得较小测量不确定度(如1)，加速度计和放大器应作为一个整体一起进行校准。加速度计应具有刚性结构。在计算测量不确定度时，应考虑基座应变灵敏度、横向灵敏度及加速度计放大器套组的稳定性(见附录A)。校准一个背靠背参考加速度计的灵敏度(幅值和或相移)时，应加装一个模拟质量块，其质量与使用背靠背参考加速度计通过比较法校准(见ISO 1606322)的加速度计质量相同。测量光点位于模拟质量块的顶部(外表面)或参考加速度计的上表面。如果检测模拟质量块顶部运动，质量块的顶部应具有抛光的表面，测量光点的位置应靠近
21、该表面的几何中心。如果质量块的运动不能视为刚体运动，应考虑其顶部与底部表面的相对运动。另一方面，也可以通过模拟质量块上的纵向通孔测量参考加速度计的顶端运动。用于比较法冲击校准的一些传统冲击机，被校加速度计的运动是通过背靠背参考加速度计的基准面传递的，两个加速度计固定为一个刚性结构。在这种情况下，背靠背的参考加速度计在进行冲击运动校准时检测紧靠加速度计的运动部位(如，棒的端面)。5GBT 2048513-2007ISO 16063-13：20015环境条件校准应在以下环境条件下进行a)室内温度：233；b)相对湿度：最大75。6优选的加速度殛脉冲持续时间加速度(峰值)和冲击脉冲持续时间应按以下系
22、列值优先选用：a)加速度值，ms2：100，200，500，1 000，2 000，5 000，10 000，20 000，50 000，100 000。b) 冲击脉冲持续时间，ms；005，01，02，051，2，5，10。注意：冲击校准脉冲持续时间要比厂家提供的最窄冲击脉冲持续时间要长一些，以避免损坏被校加速度计。7方法71测量步骤按照图1、图2、图3所示安装。调整激光干涉仪(见图3)，使其输出为两路相位正交的信号“，和地，允差满足46的规定。在做校准之前，应测量干扰量如交流声和噪声，其值应尽量小以满足校准不确定度要求。把干涉仪(46)调整到最佳状态、设置放大器档位后，应优选推荐(见第6章
23、)加速度峰值和脉冲持续时间，按72和73进行校准。注1：Fh数字信号处理程序和参考文献7、8所给出的特定的外差干涉技术可以产生允差满足46所述的正交信号，可实现比零差技术更小的不确定度。适当的外差技术在宽带(见46)干涉信号的光电传输中具有特殊的优势，但成本比零差技术高。注2：在测量时变加速度信号采用下列手段时就无须产生正交信号t单通道输出的迈克尔逊干涉仪一台时间间隔分析仪和从测得的位移值及对应时同来确定加速度值的特定数学方法(见GBT 204851f详见参考文献9)。72数据采集应该设置低通滤波器和高通滤波器(如使用)的截止频率，使低、高通滤波对校准结果的干扰影响在允差范围以内(见参考文献6
24、)。采样率的设置应满足奈奎斯特(Nyquist)定律，使最高频率小于采样频率的二分之一。在一个测量周期t。f岛+T。中对正交信号进行等间隔采样，开始于岛，结束于to+Tu。测量应在加速度计安装面上的冲击发生前时刻开始(推荐的时移为01T，T为脉冲持续时间)，在反射脉冲到达之前结束(推荐的时移为005T)。在totta+TM。内采样得到的(“。(“)和蚍(t。)两路测量值序列。其采样间隔Att一ti-1为常数。加速度计输出值的采样序列为u(t。)。将数据传送到计算机的存储器。73数据处理731概述图4给出了计算冲击灵敏度(加速度计输出峰值比加速度峰值，见引言)的三种数据处理方法的框图，三种方法的
25、区别在于使用或不使用DFT(离散傅里叶变换)：a)不使用DFT方法；6GBT 2048513-2007IS0 16063-13：2001b)速度值DFT方法；c)位移值DFT方法。对于方法a)和b)，大部分计算步骤都相同，不同的仅是在由速度序列计算加速度序列时一个是在时域微分，一个是在频域微分。在频域计算时。应用了DFT和IDFT(离散傅里叶逆变换)，而中间数据可用于加速度计灵敏度幅值和相位校准(见附录c)。方法c)仅在频域对位移序列s(t。)进行两次微分得到加速度序列口(t。)。方法a)用于“较长”的冲击，例如T-口。-o5 ms的正弦平方加速度波形(见42)和Tnm15 ms的高斯速度波形
26、(见43)，T为脉冲持续时间，以s(秒)为单位。n，k为加速度峰值，以ms2为单位。对于“较短”的冲击，可以用方法b)和方法c)。两种方法也可用于长持续时间的冲击，可有效地抑制高频干扰振动的影响。如果在冲击前后加速度计的速度不是零，使用特殊的数据处理过程(剪切处理)可正确进行傅里叶变换。注1：冲击脉冲持续时间的定义和关系见参考文献61。对于高斯速度冲击冲击脉冲持续时间就是速度大于0606vp,d的时间段。注2：上述关系下给出的冲击脉冲持续时间丁应考虑46所给的允差范围内的正交信号的失真而带来的影响(见参考文献6)。注3t参数T和n础的乘积约等于速度峰值，干涉仪输出瞬时频率和速度呈线性关系。如果
27、上述条件满足，可以用低通滤波器对干涉仪输出信号进行滤波得到满足第3章不确定度要求的信噪比。732不使用DFI方法计算冲击灵敏度通过下述的a)至i)数据处理方式计算出加速度计的冲击灵敏度，见图4的a)。a)利用采样得到的干涉仪输出值“。(f。)和批()。由下列公式计算调相值序列铷“(。)；一，、弛od(f。)一arctan竺兮*+n(2)式中：n一0，1，2，为了使调相值序列(“(。)连续。应选择nz项中的n为整数。注：关于整数n计算的程序见参考文献10。用下式计算位移序列SD(；)：跖()一鑫阳od(tz)下标D表示受高频噪声影响的序列值。a)不使用DFT的计算方法竹，。=鼽(，-)礼，=他(
28、)b)速度值DFT方法图4计算冲击夏敏度的数据处理框图(3)GBT 2048513-2007IS0 1606313：2001勘Cti) 亿。=他() R一钆()c)位移值DFT方法田4(续)b)选择数字低通滤波器的算法和参数，能够抑制高频噪声而不使信号发生畸变，对位移序列s。(t。)进行滤波，得到一“平滑”位移序列，记作s(t。)。注：具有平坦幅值响应的滤波器如4阶递归巴特沃斯(Butterworth)低通滤波器可实现这一目的。为了抑制典型噪声截止频率在正弦平方加速度波形(见42冲击机)时不应大于16T，在高斯速度波形(见43冲击机)时不应大于5T。然而为了保持信号本身的失真不被忽略截止频率也
29、不应太小。c)对位移一时间函数进行一次微分获得速度一时间函数，即速度序列(。o(，)。注：在t时刻用下式计算一次微分t口D(t。)一焘j(H1)一s(t1) “(4)d)选择数字低通滤波器的算法和参数。能够抑制高频噪声而不使信号发生畸变，对速度序列”o(f。)进行滤波，得到一“平滑”速度序列，记作v(t。)。注：同b)注。e)对速度一时间函数进行一次微分获得加速度一时间函数，即加速度序列a(t，)。注：在t时刻用下式计算一次微分ta(t)一去口(“)一v(t，1) (5)f)从计算所得的加速度计输入值的序列a(t，)判断最大值：maxa(t，)，作为加速度峰值d，-。g)选择数字低通滤波器的算
30、法和参数，能够抑制高频噪声而不使信号发生畸变，对采样得到的加速度计输出序列UD(。)进行滤波，得到一“平滑”输出序列，记作u(t。)。注：具有平坦幅值响应的滤波器如4阶递归巴特沃斯低通滤波器可实现这一目的。h)从滤波所得的加速度计输出值的序列“(t。)判断最大值：maxu(t)，作为加速度计输出峰值“。k。如果信号有零漂，冲击之前的最近零点和冲击之后的最近漂移零点应连接成一条直线，作为输出的最终零轴。允许有相对于输出峰值1的最大零点漂移。如果零点漂移过大。应考虑对不确定度的影响，给出零点漂移数值。i) 由步骤h)和f)得到“，k、n，-，用下式计算冲击灵敏度sS“一兰唑 (6)oHk在校准结果
31、的报告中，应包含按照附录A计算出的校准中的测量扩展不确定度。733速度值Dgr方法计算冲击炅敏度通过下述的a)至j)步骤数据处理方法计算出加速度计的冲击灵敏度，见图4的b)。a)同732的a)；b)同732的b)；c)同732的c)；d)对步骤c)得到的v(t。)进行DFT计算，得到速度的复频谱；e)速度复频谱乘以j2“，得到加速度的复频谱；f)用IDFT计算加速度序列a(t；)；8GBT 2048513-2007ISO 16063-18：2001g)同732的f)；h)同732的g)；i)同732的h)；j)同732的i)。784位移值DFI方法计算冲击灵敏度通过下述的a)至h)步骤数据处理
32、方式计算出加速度计的冲击灵敏度，见图4的c)。a)同732的a)b)同732的b)c)对步骤b)得到的s(t。)进行DFT计算，得到位移的复频谱d) 给位移复频谱乘以(j2n，)2得到加速度的复频谱；e)同733的f)；f)同732的f)；g)同732的g)；h)同732的h)；i)同732的i)。8校准结果报告在校准结果报告中，除校准方法以外，至少应对下列校准条件和特性参数给予说明a)环境条件：加速度计的温度环境空气温度。b)安装方法：安装表面的材料；安装力矩(如果加速度计是螺栓安装)，油或油脂(如果使用)；电缆固定；安装方位(垂直或水平)。c)模拟质量块(如果使用)：材料(如：钢)、尺寸(
33、长度、直径)和质量；安装力矩；为了补偿上下表面之间相对运动的影响，应给出灵敏度修正系数值。d)激光反射：反射镜(抛光的棒端面)；激光光点在反射面的位置。e)所有放大器的设置(如果可调)，例如：增益滤波器的截止频率。f)校准结果：冲击加速度幅值和冲击脉冲持续时间；冲击灵敏度值；测量的扩展不确定度，如果包含因子不等于2，给出的值。GBT 20485。13-20071SO 16063-13：2001跗录A(规范性附录)校准中测量不确定度的衰述A1计算给定加速度峰值、冲击脉冲持续时间、放大器增益设置和滤波器截止频率时的冲击灵敏度s。的扩展不确定度u。(sm)应按照以下GBT 204851中的公式，计算
34、给定加速度峰值、冲击脉冲持续时间、放大器增益设置和滤波器截止频率时的冲击灵敏度测量的相对扩展不确定度u。(s。)：U州(S“)一kuI(S，h) (A1)式中，包含因子一2。不确定度分量见表A1。瓶)=掣一瓦1裹A1(A2)引入灵敏度不标准不确定度分量不确定度来源 确定度的分量“(z)毗(y)加速度计输出电压峰值测量(波形记录仪；如：ADC的1 “(“v) “1(S“)分辨力)电压滤波对加速度计输出电压峰值的影响(频率带宽2 “(“kF) “。(S“)的限制)电压干扰对加速度计输出电压峰值的影响(如：交流3 “(kD) “3(Sm)声和噪声)横向、摇摆和弯曲加速度对加速度计输出电压峰值的4 (
35、“mT) “。(S“)影响(横向灵敏度)干涉仪正交输出信号的干扰对加速度峰值的影响5 “(8附k,ct) (Sn)(如t偏移、电压幅值偏差、与90名义角度差的偏差)干涉仪信号滤波对加速度峰值的影响(频率带宽的6 “(口眦k。F) “6(S“)限制)电压干扰对加速度峰值的影响(如：光电铡量链中的7 “(4“vD) 蜥(S“)随机噪声)运动干扰对加速度峰值的影响(如t偏移，加速度计参8 “(n。止MD) 口s(S“)考平面与干涉仪光点之间的相对运动)相位干扰对加速度峰值的影响(如t干涉仪信号的相9 “(4。mPD) 蛳(S“)位噪声)10 “(口附kRE) 其他干涉技术对加速度峰值的影响(干涉仪功
36、能) “lo(S，h)对冲击灵敏度测量的其他影响(如t在传感器或冲击ll u(S“RE) 机中的共振激励的影响重复测量中的随机影响算 “11(Sd)术平均值的实验标准偏差)一)S(一“TGBT 2048513-2007ISO 1606313：2001A2计算加速度峰值和冲击脉冲持续时间整个范围的冲击灵敏度s“。的扩展不确定度u。-(s“。)按照式(A1)计算的冲击灵敏度测量的扩展不确定度u。-(s。n)，只对校准的加速度峰值、冲击脉冲持续时间、放大器增益和滤波器截止频率的设置有效。在连续校准过程中的任何时间间隔里，应按照以下GBT 204851中的公式计算加速度峰值和冲击脉冲持续时间整个范围内
37、的冲击灵敏度s“。的相对扩展不确定度：U删(S札。)一hI(S。) (A3)式中，包含因子k=2。不确定度分量见表A2。一csm扣罐一瓦1裹A2引入灵敏度标准不确定度分量不确定度来源 不确定度的分量u(z)。(y)在参考峰值、冲击脉冲持续时间和放大器增益设1 u(S女) 置情况下，按照式(A1)计算的冲击灵敏度的不 “l(S也t)确定度参考放大器跟踪特性(对不同放大器设置的增益2 “(oT) 地(S“。)和相移的佣差)与参考放大器恒定的幅值一频率特性和线性的相3 M(oL L) (Sb。)位-期率特性的偏差与参考加速度计恒定的幅值一频率特性和线性的4 “(乱f*P) 地(S，h；)相位一频率特
38、性的偏差5 u(eL。) 参考放大器幅值线性的儡差 os(S,ht)6 “(丸山P) 参考加速度计幅值线性的偏差 “日(Slht)参考放大器增益的不稳定性和源阻抗对增益和7 u(eI) “7(S“I)相移的影响8 “(01P) 加速度计灵敏度的不稳定性(幅值和相移) “8(S“I)9 “(8EA) 环境对参考放大器增益和相移的影响 49(S“t)环境对参考加速度计灵敏度(幅值和相移)的10 u(eEP) “lq(Sh。)影响面辱GBT 2048513-20071SO 1606313：2001B1概述附录B(资料性附录)程序介绍按照GBT 204851，对承受平行于规定轴正弦激励的加速度计，其加
39、速度复灵敏度s。定义为：S一S一一 (B1)式中：s=：S加速度灵敏度的幅值；口加速度计正弦输出“的幅值(最好是加速度计放大器套组的输出电压)；d正弦加速度a的幅值；礼输出的初相角吼加速度的初相角；(礼一亿)复灵敏度的相移P，P一吼一10 (B3)由于复灵敏度的幅值和相移取决于频率，所以必须在任何一个所关心的频率下确定。在GBT 20485112006中描述了由正弦激励进行加速度计绝对校准的方法。(复)灵敏度的定义和校准基于线性系统理论。应用这个理论时，可以使用冲击激励的绝对校准作为一种替代方法，来确定复灵敏度的幅值和相移。如果加速度计的幅值线性偏离了理想值，冲击校准结果就可能偏离振动校准值。
40、可以按照附录C进行加速度计输入和输出的傅里叶谱计算，并在任何频率分量下可按照公式(B2)和公式(B3)分别得出复灵敏度的幅值和相移。在谱的某些频率点，n下，通过观察加速度峰值对复灵敏度幅值s。和相移他的影响，可以确定与加速度计(带或不带电荷放大器)理想幅值线性的偏差值。根据测量范围(见第1章)。本部分对加速度某些参数(峰值、冲击脉冲持续时间和冲击脉冲形状)测量使用冲击激励，用加速度计的输出峰值与输入峰值的比值来确定加速度计(带或不带电荷放大器)的冲击灵敏度。在这种情况下，应用公式(B2)时，用峰值“pH-和n，“代替口和d，灵敏度用s。表示：S“=“peak (B4)*“在这种情况下，灵敏度取
41、决于冲击脉冲持续时间，如果加速度计与理想的幅值线性发生偏离，则灵敏度还取决于峰值。当冲击机产生一个加速度。()时，加速度计输出信号满足关系式“()。第一个光电接收器的输出可以写为：“1()一nlcosouca()一位lcos90+仇d() (B5)式中，调制相位9M。d一+于M(幻 (B6)是光电接收器信号初相位角吼和调制项釉(#)合成。而彻()与位移s(幻成比例；：!至；堕 (B7)它是假设在位移s()和正弦相位项孙()之间投有时间延迟。下面说明位移是如何计算的，这样12GBT 2048513-2007IS0 1606313：2001速度”和加速度n可以基于定义从位移s中得出：础)dr矿(t
42、)一等 另一个光电接收器正交输出可以描述为：“2()一n2 sinq)M耐()一n2sin90+9M()式中：n2一口1。在t。岛+T。测量周期内，以恒定的采样速率对正交信号进行采样。在to岛+Tu。内采样得到的“。(#。)和“。(t。)测量值序列具有血一一f一一1，Iconst(常数)的一个采样间隔，其中，I为采样频率。对于两个正交信号，在测量过程中可使用下列公式连续计算它们的相位值蛳“()：9Moa(t1)_arct“搿+n(BB 10)其中n一0，1，2，这些计算程序，尤其是具有连续的“相位展开”(见参考文献10)的反正切函数的计算程序都是数字信号处理中的标准程序。使用由公式(B9)得到
43、的()值，离散的位移值序列(50(t，)可以通过下式计算：SD(t)一轰舳() ”(B11)如果90o，则s。(。)存在偏差。因为按照公式(B8)求导数，这个常数项不影响速度和加速度的测量结果。下角标D表示这些值受到高频噪音的干扰。这些干扰可通过离散时间的低通滤波器来抑制。要获得复灵敏度的幅值和相移(额域形式)和冲击灵敏度(时域形式)，需要进行B2和B3中叙述的运算。B2离散时间的低通滤波方法为实现低通滤波，可以根据输出值y(t。)y(t z)、输人采样值x(t；)、x(t一)、x(t一)，使用一种任意时刻f，输出值的递归计算算法：y(t，)=Aly(t1)+A2y(t，2)+Box(t。)+
44、B1x(t，1)+B2x(t，2)(B12)式中：i一0，1，N一1A，、Az、B0、B1、Bz是滤波器系数。这个二次微分方程的频率响应函数由下式确定： Hc加，一譬豢苇式中：n标称的角频率，,q=2xf，。为获得该滤波器通带和阻带时的单调频率响应，使用由下式给出的巴特沃斯近似函数：渺川2一可1 t碚an声tan研14+ 詈l l等)l式中n。为标称的截止频率。滤波器系数A。、A。、B。、B。、B。均根据这个近似函数确定。当应用相关计算公式时，依据理想的截止频率n。得到它们的值。对输出信号的递归计算会造成非线性的相移结果。为满足低通滤波器无任何相移的技术要求，使用下式对输出信号进行反方向的重复
45、计算：GBT 2048513-20071S0 16063-13：2001Y。(肛l一。)一Aly。(。一。)+A2y(肛)+B。，(-l一)+B1y(tN-z一；)+B2y(o肛=H)(B15)式中：i一0，l，N一1。则合成的频率响应函数可以表示为：H(e柏)H。(一)一I H(一)l 2 (B16)这样，在信号，。(。)和x(t，)之间应用两次滤波运算之后，就不会出现相移。B3微分的方法多种方法可用于求导。如果高频干扰的影响能够得到有效抑制，可使用本部分推荐的一种最简单的方法。为此目的，在B2介绍了离散时间的低通滤波方法。在对由公式(B10)得到的位移值进行离散低通滤波之后，形成以(s(t。)表示的“平滑”的位移值序列，根据下式可以进行一次微分计算：”。(f，)一面1*)一s(t-i
