振动型能量采集器的非线性动力学参数辨识方法技术

本发明专利技术涉及一种振动型能量采集器的非线性动力学参数辨识方法，属于能量采集领域。解决了当采集器含有未知非线性恢复力情况下无法准确估计采集器参数的问题。基于采集器恢复力方程，利用实验数据构建恢复力三维数据点集，通过内插值法构造恢复力曲面的截面，通过升序重排数据实现刚度曲线和阻尼曲线的辨识；采用最小二乘法，利用多项式方程实现刚度或阻尼数据点的拟合。本发明专利技术由于不需要预判采集器恢复力模型，因此适合各种复杂的非线性采集器的动力学参数识别。与现有技术相比，本发明专利技术解决了恢复力模型未知情况下，振动能量采集器的动力学模型辨识和参数识别问题，为振动能量采集器的设计和优化奠定了理论基础。

【技术实现步骤摘要】

本专利技术涉及振动能量采集领域，特别是一种振动型能量采集器的非线性动力学参数辨识方法。
技术介绍
无线传感节点和监测单元等低功耗系统目前被广泛地应用在结构健康监测、故障诊断和早期故障预警系统中。从环境中采集振动能量是为其提供电源重要的途径之一，形成研究热点。因此，发展和完善振动型能量采集器系统辨识和参数识别方法，对振动型能量采集器的结构设计、优化参数和工程应用具有重要的价值和意义。目前，在振动能量型能量采集器的动力学参数辨识或识别方面，需要对恢复力模型本身有一定预判，预判的准确性很大程度上影响了参数识别的准确性；由于未知非线性恢复力的存在，这种预判很难与实际情况相符。因此，对于具有非线性恢复力的能量采集器进行参数识别是一个难点。
技术实现思路
针对现有技术存在的缺陷，本专利技术的目的是提供一种振动型能量采集器的非线性动力学参数辨识方法，解决当采集器含有未知非线性恢复力时，无法准确识别动力学参数的问题。为达到上述目的，本专利技术的构思如下：恢复力表达式为： f ( z , z · ) = - m z ·· b - m z ·· - η u - - - ( 4 ) ]]>给予采集器随机激励，在第i个采样时刻，若质量m事先称量得到，机电耦合系数η由压电材料特性得到，加速度和电压ui也已测出(位移zi和可以通过对加速度积分得到)，则每个采样时刻下的fi均可得到。构造三维数值点集则可以绘出该系统恢复力曲面。若弹性恢复力与阻尼恢复力，则有： f ( z , z · ) = f d ( z · ) + f s ( z ) - - - ( 5 ) ]]>根据公式(5)，通过在恢复力曲面上构造或z＝0的截面，即可获得刚度恢复力数据序列或阻尼恢复力数据序列。根据上述构思，本专利技术采用如下技术方案：一种振动型能量采集器的非线性动力学参数辨识方法，利用实验数据构造三维恢复力曲面，利用截面法获得采集器刚度和阻尼数据序列，包括如下步骤：步骤1：遍历点集搜寻并记录所有当时所对应的数据序列：位移(zi,zi+1)、速度和恢复力(fi,fi+1)；步骤2：通过在以上数据对中内插值，搜寻时的位移zi’和恢复力fi’数据，得到刚度曲线数据序列(zi’,fi’)；步骤3：对数据序列(zi’,fi’)按照zi’进行升序排列，绘制刚度曲线z’→f’；步骤4：采用最小二乘法，利用多项式方程对刚度曲线z→fs进行拟合，得到多项式系数，完成恢复力模型辨识和参数识别；步骤5：识别阻尼曲线时，重复步骤1～4，但需要记录所有当zi×zi+1<0时所对应的数据序列，并在zi’＝0处插值。与现有技术相比，本专利技术具有如下的优点：本方法首先通过实验测试的手段，构造恢复力数据点集；通过插值方法，完成刚度曲线和阻尼曲线的辨识工作；最后利用多项式方程，利用最小二乘法拟合得到各项系数，完成整个辨识和识别工作。该方法由于不需要预先假定恢复力模型，因此对各种复杂非线性恢复力形式均能可靠辨识和识别。附图说明图1为振动型能量采集器等效模型。图2为振动型能量采集器动力学学参数识别结果：(a)为恢复力曲面；(b)为刚度拟合曲线与实验数据散点。具体实施方式下面结合附图对本专利技术实施例作进一步说明。振动型能量采集器等效模型如图1所示，模型由质量块m，非线性弹簧kn，阻尼系数c和机电耦合系数为η的压电换能器组成。za为绝对坐标下，采集器的振动位移，zb为采集器外壳的基础激励，则相对位移z可以写作方程(1)z＝za-zb (1)系统方程可以写为： m z ·· + f ( z , z · ) + η u = - m z ·· b - - - ( 2 ) ]]> - η z · + C p u · + u / R L = 0 - - - ( 3 ) ]]>其中：为与kn和c有关的未知恢复力函数，η为机电耦合系数，Cp为压电材料等效电容，RL为外接负载电阻，u为输出电压。将公式(2)改写为恢复力函数的表达式(4)： f ( z , z · ) = - m z ·· b - m z ·· - η u - - - ( 4 ) ]]>在第i个采样时刻，若质量m事先称量得到，机电耦合系数η由压电材料特性得到，加速度和电压ui也已测出(位移zi和可以通过对加速度积分得到)，则每个采样时刻下的fi均可得到。构造三维数值点集则可以绘出该系统恢复力曲面。若弹性恢复力与阻尼恢复力，则有： f ( z , z · ) = f d ( z 本文档来自技高网...

【技术保护点】
一种振动型能量采集器的非线性动力学参数辨识方法，利用实验数据构造三维恢复力曲面，利用截面法获得采集器刚度和阻尼数据序列，其特征在于，包括如下步骤：步骤1：遍历点集搜寻并记录所有当时所对应的数据序列：位移(zi,zi+1)、速度和恢复力(fi,fi+1)；步骤2：通过在以上数据对中内插值，搜寻时的位移zi’和恢复力fi’数据，得到刚度曲线数据序列(zi’,fi’)；步骤3：对数据序列(zi’,fi’)按照zi’进行升序排列，绘制刚度曲线z’→f’；步骤4：采用最小二乘法，利用多项式方程对刚度曲线z→fs进行拟合，得到多项式系数，完成恢复力模型辨识和参数识别；步骤5：识别阻尼曲线时，重复步骤1～4，但需要记录所有当zi×zi+1<0时所对应的数据序列，并在zi’＝0处插值。

【技术特征摘要】
1.一种振动型能量采集器的非线性动力学参数辨识方法，利用实验数据构造三维恢复力曲面，利用截面法获得采集器刚度和阻尼数据序列，其特征在于，包括如下步骤：步骤1：遍历点集搜寻并记录所有当时所对应的数据序列：位移(zi,zi+1)、速度和恢复力(fi,fi+1)；步骤2：通过在以上数据对中内插值，搜寻时的位移zi’和恢复力fi’数据，得到刚度曲线数...

【专利技术属性】
技术研发人员：陈立群，袁天辰，杨俭，丁虎，
申请(专利权)人：上海大学，
类型：发明
国别省市：上海;31