本发明专利技术涉及一种基于时变公分母模型的时频域时变结构模态参数辨识方法，属于结构动力学技术领域。首先对具有时变特性的航空器或航天器结构在工作环境下测量得到的结构动力学响应信号进行时频分析，得到非参数化估计的对应时变结构的时间相关功率谱函数，然后以时变公分母模型为时变结构动力学的参数化模型，通过时频域的最小二乘方法估计出时变公分母模型的待估参数，最后利用估计出来的时变公分母模型的待估参数计算出对应时变结构的模态频率和模态阻尼比。本发明专利技术适用于航空和航天工程应用领域的时变结构的模态参数辨识，并且所需使用者的参与度较低，具有使用简单和方便的特点。
【技术实现步骤摘要】
本专利技术涉及一种，属于结构动力学
技术介绍
在现实的生产和生活中，许多工程结构表现出这样的时变特征，如列车激励中的车桥系统、飞行过程中液体燃料逐渐减少的运载火箭、气动力附加效应下的飞机、柔性可展开的几何可变航天器、旋转机械等。在国内航天领域，大型空间站、新一代运载火箭、大柔性展开式卫星等新一代的航天器已被列入我国最新的航天发展规划中，成为未来几十年中国航天器发展的主要方向。 大型空间站、新一代运载火箭、大柔度展开式卫星的结构在运行中无一例外的存在着较强的时变因素，如未来大型空间站的空间对接问题，现役和未来新一代运载火箭的燃料质量快速消耗，以及大柔度可展开式卫星的空间展开等。因此，作为时变结构动力学特性分析的重要方法和途径，时变结构模态参数辨识研究将成为未来航天器结构动力学研究的重点之一。时变结构模态参数辨识可以辨识时变结构的模态频率、模态振型和模态阻尼，这些参数具有重要的物理意义，可以为时变结构的结构设计、结构健康监测、结构故障诊断、结构振动控制等方面的应用提供有力的支持。按照采用的数学模型的差异区分，现有的时变结构模态参数辨识的方法主要有四类第一类是从传统的时不变结构模态参数辨识发展而来的基于在线递推技术的时变模态参数辨识方法。这类方法的基础是传统的时不变结构模态参数辨识方法，不同之处为在每一个时刻数据序贯地被考虑，老的数据逐渐被遗忘，新的数据不断地加进来，模态参数的估计值在每一个时刻被修正。这类方法存在两方面的缺陷第一，存在着观测数据及遗忘因子(算法)的选取问题，需要在识别精度和跟踪能力二者之间做折中，并且对于不同结构的相关选取的适应性也很难解决；第二，这类方法来自传统的模态参数辨识方法，需要结构的输入和输出两方面的响应信息，因此很难运用于如在飞航天器等只能得到输出响应信号的结构模态参数辨识。第二类是基于短时不变假设的模态参数辨识方法。这类方法将数据(结构响应)划分成一个个小的时间段，并在每一个时间段内把结构参数看成是时不变的，然后将每一段内识别值用一定的数据处理技术(如曲线拟合) 处理得到模态参数随时间变化的规律。它的特点是估计后一段时间的模态参数时没有用到前面各段的数据信息，对参数变化较快的结构为使估计精度提高必须选取很短的数据段。 此方法包括现今较为常用的基于状态空间模型的递推的随机子空间辨识法(N4SID)和时间相关自回归滑动平均模型(Time-d印endent ARMA, TARMA)方法。这类方法的时变结构模态参数辨识方法发展时间最长，发展的也最为完善。但是一些固有的问题限制了其进一步发展和应用第一，短时不变假设限制了此类方法对于快变、突变参数辨识方面的应用； 第二，此类方法需要形式固定、明确的数学模型，如状态空间模型、时间序列的自回归滑动平均模型等，因此，在辨识中模型的定阶问题十分突出，模型阶数的不确定将引入无物理意义的虚假模态，造成辨识结果不可用，模型阶次合理选取、虚假模态的判断等问题更需要进一步深入研究；第三，作为两种主流的基于短时不变假设的模态参数辨识方法——递推的随机子空间辨识法和时间相关自回归滑动平均模型各自存在着一些其它的问题基于状态空间模型的堆积子空间方法不可避免地要使用QR分解、特征值分解(EVD)或者奇异值分解 (SVD)技术，这必然会带来方法数值实现上的复杂性，对于大型工程结构，尤其对有在线以及快速辨识要求的问题，这还需要进一步进行研究；基于时间序列模型的辨识方法研究都不能回避参数跟踪算法的设计，虽然各种改进的最小二乘法、各种滤波方法不断提出，但是当相同模型使用不同跟踪算法，以及不同模型应用相同算法结果差异非常大。第三类是人工神经网络的时变模态参数辨识方法。人工神经网络已经被广泛地应用于非线性系统辨识问题，并取得良好的效果但大部分研究工作还仅局限于时不变系统，只是近几年来被推广到时变系统。将人工神经网络用于时变模态参数辨识领域研究的公开发表的文献很少，其主要集中在针对简单结构(系统)的机理性研究。对于真实的复杂结构还存在算法复杂、计算效率低和辨识精度差等问题。第四类方法是基于时频分析的非参数化时频域的时变结构模态参数辨识方法。从信号分析的角度来看，时变结构在工作环境下的结构动力学响应信号是非平稳随机信号。经典傅立叶变换经过一个世纪的发展，已成为信号处理领域最强有力的分析方法和工具，这主要是由它的正交性和鲜明的物理意义以及快速简洁的计算方法所决定的。但是，由于傅立叶变换是对时间求积，去掉了非平稳信号中的时变信号，因而要求信号是平稳的，对时变非平稳信号难以充分刻画。为了满足对突变信号、非平稳信号分析的要求，1946 年，Gabor提出了加窗傅立叶变换分析方法,亦称短时傅立叶变换(short-time Fourier transform, STFT),通过适当窗函数的选取,就可以实现一定程度上的时频分析，但是由于时间分辨率与频率分辨率要受到窗函数宽度的限制，总是不能同时到达最佳。1948年， Ville提出了著名的维格纳-威尔分布(ffigner-Ville distribution, WVD)。它作为一种能量型时频联合分布，与其他时频分布相比有许多优良性质，如真边缘性、弱支撑性、平移不变性等，是一个非常有用的非平稳信号分析工具。由于多信号的维格纳-威尔分布出现交叉项，在不少场合会限制其应用效果，所以后来研究人员在此基础上，提出了多种改进形式，如指数分布、广义指数分布、广义双线性时频分布等，其中广义双线性时频分布又称为科恩类能量型时频分布。后来在此基础上，人们又提出了科恩类时频分布等方法，这些时频分析方法在非平稳随机信号分析中得到了广泛的应用并取得了许多令人满意的结果。近十年，由于时频分析在非平稳随机信号分析方面的优势，越来越多的研究者运用时频分析来进行时变和非线性系统辨识的研究。时频分析方法对时变和非线性结构模态参数进行辨识也渐渐成为模态参数辨识研究领域的热点之一。2000年Ghanem将结构动力学控制微分方程在一系列小波基上展开，用小波系数来代替原来的物理响应，并采用了求解展开方程的方法辨识了系统的模态参数；2003年Zhang和Xu通过对一个简单的时变结构响应的Gabor变换辨识了结构的模态频率；2007年Roshan-Ghias采用解析推导的方式对一个单自由度系统和一个三自由度系统自由振动下的响应进行了 WVD和SPWVD变换，并根据变换结果估计了系统的模态频率和阻尼比。现有的基于时频分析的时频域时变结构模态参数辨识方法都是非参数化的，虽然有的方法能够很好的辨识出时变结构的模态频率，但是非参数化的方法不同程度的依赖于使用者的主观意识和经验，并且对于随机激励下的时变结构模态阻尼比辨识还没有很好的办法。
技术实现思路
本专利技术针对航空器和航天器时变结构模态参数辨识问题，提出了一种时变公分母模型的时频域时变结构模态参数辨识方法，其基本思路为首先对具有时变特性的航空器或航天器结构在工作环境下测量得到的结构动力学响应信号进行时频分析，得到非参数化估计的对应时变结构的时间相关功率谱函数，然后以时变公分母模型为时变结构动力学的参数化模型，通过时频域的最小二乘方法估计出时变公分母模型的待估参数，最后利用估计出来的时变公分母模型的待估参数计算出对应时变结构的模态本文档来自技高网...
【技术保护点】
基于时变公分母模型的时频域时变结构模态参数辨识方法，其特征在于：具体实现步骤如下：步骤1，根据被辨识航空器或航天器时变结构工作状态、感兴趣的时间范围以及被辨识航空器或航天器时变结构的主要频率范围，设定辨识所需的采样时间和采样频率，并对辨识结构的结构动力学响应信号进行采集；步骤2，从步骤1采集的响应信号中任意选取参考信号，并结合参考信号分别对各个响应信号进行时频分析，得到响应信号的时间相关功率谱函数Gk(tτ，ωf)；其中，tτ表示时间采样点，ωf表示频率采样点，下标τ＝1，2，...，Nτ，f＝1，2，...，Nf，Nτ为时间采样数，Nf为频率采样数，k＝1，2，...，NsNr，Ns为输出响应数，Nr为响应信号参考点数；步骤3，根据采样时间和采样频率以及被辨识航空器或航天器的时变特点，建立时变公分母模型如下：G^k(tτ,ωf)=Bk(tτ,ωf)A(tτ,ωf)---(1)其中，分子多项式和分母多项式分别为：其中，为时频基函数，时间多项式和频率多项式分别为i和j阶，i＝0，1，2，...，nt，j＝0，1，2，...，nω，nt为时间多项式阶数，nω为频率多项式阶数，分子多项式系数bk，i，j和公分母ai，j写为向量形式：Bk,j=[bk,0,j,bk,1,j,...,bk,nt,j]T,Aj=[a0,j,a1,j,...,ant,j]T---(4)令：βk=[Bk,0T,Bk,1T,...,Bk,nωT]T;α=[A0T,A1T,...,TnωT]T---(5)则有：θ=[β1T,...,βkT,...,βNoNiT,αT]T---(6)其中，θ为时变公分母模型中的待估参数向量；步骤4，采用最小二乘参数方法得到时变公分母模型中的待估参数向量θ；具体过程如下：最小二乘的花费函数为：lLS(θ)=Σk=1NrNsβkTαTPkQkQkTRkβkα---(7)其中，Pk=Re(ΘkHΘk),Qk=Re(ΘkHΞk)和Rk=Re(ΞkHΞk);其中Wk(tτ，ωf)为权函数；根据下式求得约束后的待估公分母参数向量α′：D′α′＝b′????????????????????(10)其中，D′=D[1:(nt+1)(nω+1)-1,1:(nt+1)(nω+1)-1],b′=-D[1:(nt+1)(nω+1)-1,(nt+1)(nω+1)],D=2Σk=1NrNs(Rk-QkHPk-1Qk);令α=α′1,并通过下式计算出待估分子参数向量βk：βk=Pk-1Qkα---(11)从而得到时变公分母模型中的待估参数向量θ；步骤5，根据被辨识航空器或航天器模态参数的应用需要，给定需要计算模态参数的时刻点tτ′，并利用步骤4得到的公分母参数向量α以及时变公分母模型，计算出给定时刻tτ′下辨识得到的模态频率fr和模态阻尼比ξr：fr=Im(λr)2π,ζr=Re(λr)|λr|---(12)其中，Im和Re分别为取括号中值的虚部和实部，λr为以α为系数的tτ′时刻下的时变公分母模型中分母多项式A(tτ′，ω)的第r个根，ω为被辨识航空器或航天器的频率变量。FDA00002330041600012.jpg,FDA00002330041600013.jpg,FDA00002330041600014.jpg,FDA00002330041600025.jpg,FDA00002330041600026.jpg...
【技术特征摘要】
【专利技术属性】
技术研发人员：周思达，刘莉，董威利，杨武，马志赛，贺媛媛，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。