第一作者：江文松

（-），男，安徽安庆人，副教授，博士，硕士研究生导师，现于中国计量大学任教，主要研究方向为精密测试理论。

通讯作者：张力

（—），男，辽宁本溪人，研究员，博士生导师，博士，主要从事动态测试与校准技术、先进传感技术研究。

论文引用格式

江文松，吴娅辉，李新良，等.动态力计量建模与分析方法综述［J］.计测技术，，42（2）：2-12.

动态力计量建模与分析方法综述

江文松1，吴娅辉2，李新良2，王中宇3，张力2*

（1.中国计量大学计量测试工程学院，2.航空工业北京长城计量测试技术研究所计量与校准技术重点实验室，3.北京航空航天大学仪器科学与光电工程学院）

　　摘要：针对动态力计量中大规模复杂数据分析与处理难度大的问题，本文从原理性和实用性角度，系统地对动态力的标准激励源、校准仪表建模方法、动态时频域评价和测量不确定度评定等进行了综述讨论。梳理了数据建模和评价在动态力校准、量值溯源与传递等关键计量链路中的作用。为动态力计量数据评价机制的建立提供了一定的理论基础。

　　关键词：计量；系统建模；数据评价；溯源；动态力

0　引言

　　动态计量数据作为一种典型的科研数据，是国家质量基础的核心内容之一［1-2］。在计量领域，测量数据的正确建模与分析是处理复杂系统量值溯源与传递的关键。

　　计量需要选取合适的动态发生信号作为标准量的激励源，用以激发被校准传感器、测量装置、实验装置等工业仪器仪表的动态特性［3-4］。例如：在动态力校准中，既可选择电磁激振器、机械激振器等产生标准正弦激励信号来激发被校准仪表的动态频域特性，也可以选择落锤式冲击力校准装置产生标准半正弦脉冲激励信号来激发被校准仪表的瞬态响应特性。

　　对被校准仪表建立数学模型既能帮助计量检定人员了解其工作原理，选择合适的计量方式，也可为生产厂家开展仪表参数的辨识和校准提供量值溯源的可能［5-6］。例如：在建立好压力传感器的二阶参数化模型之后，采用激波校准法实现其参数的辨识和二阶模型响应特性的量化［7］。

　　动态力校准装置通过施加标准激励和可溯源的量值测量，就能通过辨识的数学模型进行动态时域特性和频域特性评价。例如：通过时域上的上升时间、峰值时间或频域上的穿越频率可以评价被校准仪表的响应速度；通过时域上的调节时间、超调量、或频域上的相位裕量可以评价被校准仪表的稳定性；通过稳态误差或动态降落可以判断被校准仪表的抗干扰特性等［8-12］。

　　由于对计量活动认识的局限性，测量过程不可避免地存在测量误差、建模误差、数据分析误差等，产生具有一定量值水平的测量不确定度［13］。通过这些评价指标，实现动态量值的立体计量。例如：Heisenberg在年针对量子物理学实验首次提出了测不准原理及概念，年国际不确定度工作组联合各个专业国际组织和各国计量部门，制定了《测量不确定度表示指南》（《GuidetotheExpressionofUncertaintyinMeasurement》，简称《GUM》），被世界各国及国际组织沿用至今，并被广泛应用于航空航天、电力电子等现代高端制造行业［14-15］。

　　因此，激励源设计、被校准仪表建模、动态时频域评价、测量不确定度评定等数据处理方法是开展动态校准、量值溯源与传递等的关键环节。本文以动态力校准为例，分别通过讨论激励源的分类形式、工作原理及模型，分析动态计量参数的建模方法，建立其计量数据的时频域评价体系，介绍动态校准数据的测量不确定度评定方法，系统地概述动态力计量的建模和分析理论，为动态力计量数据评价机制的建立提供理论基础。

1　动态力激励信号的量值表征

　　在动态计量领域，通常采用阶跃信号、脉冲信号、正弦信号、斜坡信号、抛物线信号等特征数据作为典型的激励源。在动态力校准实践中，前三种类型的信号较为常见，属于常规激励信号，通常利用这些激励开展不同形式动态量值的测量、溯源、校准及评价等计量活动，如图1所示。

图1　动态量值计量过程的建模与分析方法

　　1.1　常规激励源

　　1.1.1　正弦激励信号

　　正弦激励信号r（t）属于周期变化的信号，在交流供电电源、动态力、动态压力、振动、动态温度等现场校准过程中较为常见，经常在动态校准过程中，用于信号的频谱分析，其数学表达式为

　　式中：A为信号的幅值；ω为角频率；t为时间，s；在频域表达上，r（t）需要进行拉普拉斯变换，L（r（t））为对r（t）进行拉普拉斯变换，得到其频域的表征形式值R（s）；s为频率，Hz。

　　在动态校准过程中，将正弦信号作为激励，通过频率特性法计算被校准仪表的频率响应特性，实现其幅频特性、稳定性等动态过程指标的分析。例如：正弦激振器为正弦力激励源，输出周期频率可变的动态力，如图2所示。

图2正弦激励源及信号发生原理

　　由于常规激励源具有激发动态特性的良好能力，研究人员近年来针对正弦激励源设计了动态力计量装置，并做了大量的研究。通过被校准力传感器、负载质量块和参考传感器来实现动态载荷的精确计量与校准［16］，动态力校准的矢量受力方程为

　　式中：m1为标准质量块质量，kg；m2为连接件的质量，kg。

　　正弦力校准装置最早由德国联邦物理技术研究院（Physikalisch-TechnischeBundesanstalt，PTB）提出，利用波形发生器、功率放大器和高精度振动台产生特定频率和幅值的正弦力。年，T.Bruns研究员和R.Kumme教授在此基础上对引入气浮技术、减震技术和加速度传感器网络等，实现横向偏摆、惯性质量、加速度分布不均匀等影响因素的干扰消除，在10Hz~2kHz的频率范围内实现了，N和10kN的正弦力，测量不确定度达0.4%［19］。法国国家实验室（Laboratoirenationaldemétrologieetd‘essais，LNE）和西班牙计量中心（CentroEspa?oldeMetrología，CEM）及北京长城计量测试技术研究所（ChangchengInstituteofMe?trologyMeasurement，CIMM）等也基于该原理设计了系列正弦力校准装置（如图3所示），并实现国际比对［20］。

图3不同型号正弦力计量装置

　　1.1.2冲击激励信号

　　冲击信号持续时间极短，其数学表达式为

　　式中：H为冲击幅值；ε为脉宽。通常当H=1且脉宽逼近时，即为常见的单位理想冲击函数δ（t），其面积为脉冲强度。据此可知δ（t）信号的拉普拉斯变换为1，具有瞬间无穷大的能量，因此这种信号无法在实际工程中获得，如图4所示。

图4　冲击信号与单位理想冲击函数

　　在计量实践中，通常把ε远小于被测传感器响应时间的窄脉冲半正弦信号近似为δ（t）函数。例如：落锤式发生装置可以近似为冲击力激励源，能够输出不同脉冲宽度和脉冲幅值的动态力，如图5所示。

图5　冲击激励源及信号发生原理图

　　冲击力校准装置由激励源和测量单元组成。各国也基于该激励原理实现冲击力计量装置的研制，如图6所示。例如：中国计量科学研究院（NationalInstituteofMetrology，China，NIM）20世纪80年代就开展了N~kN冲击力校准典型装置研制，其在0.6ms的脉冲宽度的测量不确定度为2%，此时自由度取k=3［21］；印度国家物理实验室（NationalPhysicalLaboratoryIndia，NPLI）在年联合德国PTB研制了1kN~1MN的冲击力校准装置，实现了，，，0kN等冲击力的溯源，重复性误差接近1%［22］；CIMM研制的冲击力校准装置，能够实现0.5ms脉宽、kN冲击力的溯源，重复性测量误差不大于1%［23］。

图6　不同型号冲击力计量装置

　　1.1.3　阶跃激励信号

　　阶跃信号为通过瞬间突变后保持的信号，其数学模型可表示为

　　阶跃信号的发生比较特殊，既需要具备瞬态突变特征，又需要有量值保持能力。例如：激波管的膜片破裂、脆性梁的断裂、电源突然通断、负载突变等均近似地认为会输出阶跃信号。另外，阶跃信号是传感器或测量系统时域性能分析的常用激励，因此被广泛应用于动态性能分析。在动态力校准中，脆性梁的瞬间断裂产生的负载通常为阶跃力，用来产生标准阶跃激励，如图7所示。

图7　阶跃激励源及信号发生原理

　　阶跃力校准主要以卸荷的方式产生，例如，德国PTB实现了20kN和kN的阶跃力发生，北京航空航天大学黄俊钦教授实现了kN阶跃力校准［24］；北京理工大学也利用该方法实现了kN负阶跃力校准［25］；北京长城计量测试技术研究所采用快速卸载方式在10~μs的上升时间实现了1.2MN负阶跃力校准等［18］。但相比于正弦力和冲击力校准方式，由于阶跃信号的保持难度较大，因此在计量实践中并不多见。

　　1.2　其他激励源

　　1.2.1　斜坡激励信号

　　斜坡信号是随时间t线性增长的信号，其数字表达式为

　　1.2.2　抛物线信号

　　抛物线信号随时间t以等加速度增长，因此也称为等加速度信号，其数学表达式及拉普拉斯变换形式可以表示为

图8　其它激励源

　　由于斜坡信号和抛物线信号作为激励源不仅不容易准确实现信号发生与放大，还很难激发被校准传感器或仪器的动态性能，因此在动态校准中并不常见。

2　动态力校准参数的建模与分析方法

　　2.1　动态系统建模

　　当确定了激励源之后，建立被校准仪表的模型是动态计量的关键环节。为此需要利用被校准仪表上的激励和响应信号，对其模型进行辨识。

　　针对原理已知的白箱系统，通常对其各个环节所遵循的基本物理规律建立微分方程组，通过消除中间变量和归一化处理，建立传感器输入与输出之间的微分方程（由输出量的各阶导数和输入量的各阶导数及系统结构参数构成），从而得到动态系统的数学模型。

　　令n阶线性定常系统，其数学模型可用n阶线性定常微分方程描述为

　　通过式（11）可得到G（s）的幅频特性A（ω）和相频特性φ（ω）。因此，只要确定了被校准力传感器的频率特性，就能辨识其关键参数水平。力传感器动态特性由其传递函数确定，采用上述的正弦、阶跃或脉冲等常规激励作用，即可求解传递函数模型的各参数，实现动态力测量过程的动态建模。

　　2.2时域评价方法

　　在动态计量中，被校准传感器的时域动态指标通常采用阶跃激励r（t）作用，评价被校准仪表的跟随性能及稳定性能。典型的被校准传感器的阶跃响应曲线如图9所示，其中，反映响应特性的指标包括上升时间tr、峰值时间tp、超调量σ%、调节时间ts、振荡次数、动态降落、恢复时间及稳态误差等。

图9阶跃响应曲线及性能指标

　　设某被校准传感器为二阶线性定常系统，即n=2，在阶跃激励r（t）作用下的输出c（t）的拉普拉斯变换为

　　通常根据ξ在［0，+∞）上取值的不同，特征值s1，2有可能为实数根、复数根或重根。被校准传感器动态时域性能评价指标分别如下：

　　1）上升时间

　　上升时间tr是系统输出响应从稳态值的10%上升至稳态值的90%所需时间。因此，根据式（8）可求得c（t）为

　　2）峰值时间

　　峰值时间tp是指传感器或系统的输出响应自零时刻开始到首次达到峰值所经历的时间。对于二阶系统，此时的c（t）在tp时刻的一阶导数为0，即

　　3）超调量

　　超调量σ%为传感器或系统的输出响应超出稳态值的最大偏离量占稳态值的百分比，即

　　由式（18）可知，σ%越小，被校准仪表的动态响应越平稳。

　　4）调节时间

　　调节时间ts是从零时刻开始，被校准仪表的输出响应进入并保持在某一误差带Δ之内所需的时间，通常Δ取±5%（或±2%）。

　　显然，误差带越小，对被校准仪表的计量特性要求就越高；调节时间越小，表明被校准仪表的动态响应过程越短，即响应越快速。

　　5）稳态误差

　　从定性的角度看，被校准仪表的输出量围绕稳态值摆动的次数越少、摆动的幅度越小，表明系统稳定性越好。这种摆动次数和摆动幅度通常用稳态误差ess来定量表示。稳态误差是指系统期望值r0（t）与实际稳态值r（t）之间的差值，即

　　由式（19）可知，ess越小，系统稳态精度越高。

　　6）动态降落值

　　动态降落值ΔCmax用于衡量被校准仪表的抗扰动性能，通过对稳态过程施加一个标准负阶跃激励，就能测量出被校准仪表响应的最大降落值ΔCmax，从扰动激励到最大降落值之间的过渡时间称为降落时间tm。

　　当采用脉冲信号、正弦信号等其它激励源时，还需要通过被校准传感器的输出计量数据的峰谷值、脉冲宽度、平均值等特征指标来评价其动态计量性能。

　　2.3频域评价方法

　　以典型二阶系统为例，无论是时域评价方法还是频域评价方法， 　　令s=jω，根据式（8）可知，被校准仪表的幅频特性A（ω）可表示为

图10典型二阶系统幅频特性曲线

　　另外，在进行典型高阶测试系统的频域特性评价时，谐振频率和阻尼比系数也常常被用于判断其系统的跟随性能及稳定性能。由于二者均可与相位裕量、穿越频率建立算法转换联系，因此在功能原理上类似，在此不再赘述。

3　测量不确定度评定

　　测量不确定度是一个体现测量值分散性的参数，它由基于概率论和统计理论的A类评定和基于先验知识的B类评定得到，是两类传统评定方法。

　　3.1　两类传统评定方法

　　3.1.1　A类评定方法

　　假设某被测量x，在相同测量条件下独立重复测量n次，观测值为xi（i=1，2，…，n），任意一次测量值xk的A类标准不确定度用其实验标准差s（xk）表示，因此，采用统计分析的方法求得的标准差即为A类评定方法，以贝塞尔公式法为例，标准差为

　　除此之外，标准差还可采用极差法、最大残差法、最大方差法、彼特斯法等进行求解。

　　3.1.2B类评定方法

　　该方法是基于信任度的主观概率方法。当输入量Xi的估计值xi不由重复观测得到时，估计方差或标准不确定度可根据Xi的可能变化的全部有关信息的判断来评定，这些信息来源可能包含：以前的测量数据、对有关材料和仪器特性的了解、生产厂提供的技术说明书、校准证书或其他证书提供的数据、手册给出的参考数据的不确定度、规定实验方法的国家标准或类似技术文件中给出的重复性限r或复现性限R。

　　在动态力校准过程中，误差主要来源于各个观测值。其中，冲击力校准装置主要测量误差源包括有效碰撞质量的测量误差、激光光轴与运动方向的夹角误差、激光干涉仪的测量误差、落锤的加速度分布不均匀引入的误差、落锤的横向偏摆引入的误差、二次仪表误差和重复性测量误差等；正弦力校准装置的测量误差主要来源于等效质量的测量误差、激光光轴与运动方向的夹角误差、激光干涉仪的测量误差、质量块的加速度分布不均匀引入的误差、激振器的横向偏摆引入的误差、激振器波形失真和重复性测量误差等。其中，重复性测量误差可采用A类评定方法求解标准差，其他误差来源可采用B类评定方法求解标准差。

　　3.2蒙特卡洛法

　　传统评定方法仅适用于高斯分布模型，当测量模型具有偏导、积分等复杂函数时，测量不确定度的传统评定方法很难解决。为此，ISO/IECGuide98-3系列标准规定了蒙特卡洛法评定测量不确定度的实施步骤。

　　已知测量模型为Y=f（X），抽样样本X=（x1，x2，…，xN），其概率分布函数为px（ξ），对应输出模型样本为Y=（y1，y2，…，yN），对应的概率分布函数为py（ξ）。蒙特卡洛法的具体步骤如下：

　　1）选择蒙特卡洛仿真试验次数

　　蒙特卡洛仿真试验次数M决定输出量的样本容量。M值越大，样本容量越大，越接近于输出量的真实情况。通常M的取值使包含区间达到95%及以上为宜。

　　2）确定输入样本的概率密度函数

　　蒙特卡洛法的评定结果依赖于测量模型，因此各输入量的概率密度函数是评定结果的约束。由于各输入量的概率密度函数依赖于数据本身规律，因此可以通过直接法、贝叶斯法、最大熵原理等统计分析方法来设定X的概率密度函数。

　　3）输入样本的概率密度函数抽样

　　根据上述概率密度函数，对输入样本X进行M次随机抽样，得到M×N维输入样本的伪随机数矩阵VM×N。

　　4）输出样本值的计算

　　将上述矩阵VM×N代入测量模型Y=f（X）中，即可得到对应的M组输出量WM×N。

　　5）输出评定结果

　　根据输出量Y的抽样WM×N，计算输出量Y的数学期望和标准差，其中数学期望即为Y的估计值，标准差即为标准测量不确定度；通过分布函数py（ξ）得到Y的包含区间和置信概率。

　　3.3其他评定方法

　　由于传统评定方法依据大数定律和中心极限定理，当先验信息缺乏或测量数据为小样本量时，GUM因缺乏理论基础会出现评估偏差和置信度下降的问题。为此，依据先验知识的贝叶斯理论、灰模型法及自助法等非统计分析方法被用于测量不确定度的评定上。

　　3.3.1贝叶斯理论方法

　　设测量模型的参数为θ，测量数据样本集为X=（x1，x2，…，xN），该样本的似然函数为L（x

θ），采用贝叶斯理论可计算后验分布，即

　　式中：c为灰度系数，由灰色GM（0，N）模型求得。

　　由此可以看出，贝叶斯法侧重于建立后验分布规律，灰模型法侧重于小样本数据的挖掘，这两种方法都旨在解决数据规律未知情况下的不确定度评价。但鉴于这类新方法在适用场合上的特殊性和对数据的二次加工处理，通常只会作为两类传统评定方法和MCM方法的补充。

　　3.4不确定度合成与表示方法

　　为了保证扩展不确定度是被测量可能值包含区间的半宽度，扩展不确定度U由合成标准不确定度uc乘以包含因子k得到，即

4结论

　　系统性地综述了动态力计量的建模和分析方法，探讨了激励源、参数建模、数据评价、测量不确定度评定等计量基础理论在实际计量工程案例中的应用。针对动态力校准过程中标准激励发生方式，探讨了各类激励信号源的原理、数学表示形式，简要介绍了不同量值溯源中的激励装置。针对被校准仪表的黑白箱模型，以动态力校准为例，分别介绍了基于物理环节的微分-传递函数建模方法和基于卷积、归一化、聚类、模态分解或正则化等信号处理的参数辨识方法。建立动态力校准模型的时域评价体系和频域评价体系，实现动态力校准模型的动态性、跟随性和稳定性的综合评价；系统地建立了动态力校准数据的测量不确定度评价方法体系，实现测量数据的可溯源。

长按识别

转载请注明:http://www.aideyishus.com/lkyy/844.html