当前位置: 首页 > 专利查询>华侨大学专利>正文

复杂声振模拟实验环境下进行不相关多源频域载荷识别的装置和方法制造方法及图纸

技术编号:12528782 阅读:110 留言:0更新日期:2015-12-17 23:56
本发明专利技术涉及一种不相关的声振载荷联合施加和不相关多源频域载荷识别研究的试验装置和方法,以及利用该装置在复杂的声振模拟环境下进行不相关多源频域载荷识别的三种方法,分别是最小二乘广义逆法、改进的正则化方法以及多输入多输出支持向量回归机法,三种方法均能根据系统上多个测点的振动响应同时识别出多个不相关频域载荷源的大小本发明专利技术的最小二乘广义逆法不需要测定传递函数的相位;改进的正则化方法对应的每一个频率都有最佳的正则化参数;多输入多输出支持向量回归机法能够有避免过学习现象。

【技术实现步骤摘要】
一种声振载荷联合施加的试验装置及其载荷识别方法
本专利技术涉及一种声振载荷联合施加的试验装置,以及利用该装置在复杂的声振模拟环境下进行不相关多源频域载荷识别的三种方法,分别是最小二乘广义逆法、改进的正则化方法以及多输入多输出支持向量回归机法,三种方法均能根据系统上多个测点的振动响应同时识别出多个不相关频域载荷源的大小。
技术介绍
在实际工程中,如导弹在空中飞行、海洋平台等大型建筑物受风浪及交通激励作用等情况下,很难对作用于结构的外载荷进行直接测量或计算,甚至有时因载荷作用点不可到达,使这种动态载荷不可测。载荷识别是通过测量结构动态响应和系统特征来求结构所受激励的方法,载荷识别技术为那些无法直接测量载荷的结构或者系统提供了一种识别动态载荷的有效方法,而准确地确定载荷、科学地制定相应的载荷谱是可靠性试验、振动主被动控制、铁路交通和桥梁设计等重大工程设计时面临的迫切问题。载荷识别是通过测量结构的动态响应和系统特性来反求结构所受激励的方法,属于振动问题中的第二类反问题,存在不适定性。从结构响应数据中反求载荷是目前研究的热点和难点。载荷识别主要分为频域法和时域法两类。其中频域法提出较早,主要利用激励和响应之间的频响函数求逆来实现,但是矩阵求逆法在应用时通常需要求解广义逆,其经常会遇到系数矩阵的病态问题和奇异值分解问题。
技术实现思路
专利技术的目的在于克服现有技术之不足,提供一种声振载荷联合施加的试验装置,用于模拟复杂的声振环境,用于载荷识别试验研究,并提供三种不相关多源载荷的频域识别方法,分别是最小二乘广义逆法,改进的正则化方法以及多输入多输出支持向量回归机法。本专利技术解决其技术问题所采用的技术方案是:提供一种声振载荷联合施加的实验装置,包括:两端加板密闭圆柱壳装置,通过弹性橡胶绳悬挂;所述圆柱壳装置内部有一个球型噪声激励源,由声传感器记录球型声源所施加的噪声激励大小,所述圆柱壳装置外部有一悬挂式振动台激振器,由振动传感器记录悬挂式振动台输入的振动激励;所述圆柱壳装置外表面和内表面布置有18个振动加速度传感器,用于测量机械结构在球型噪声源和悬挂式振动台联合激励下的振动响应;球型噪声激励源激励与悬挂式振动台激振器激励每次加载的时候位置和方向均固定不变,且18个振动加速度传感器分布在圆柱壳内部的各个地方,能反映系统的主要振动。优选的,所述独立的球型噪声激励源激励,有3种量级激励,而且量级逐渐增大;所述独立的悬挂式振动台激振器振动激励,有5种量级激励,而且量级逐渐增大;当噪声激励和振动激励联合加载时,噪声激励和振动激励的量级两两组合,形成了15种不同的量级,从而实现了模拟复杂的声振环境。对声振实验装置分别加载15种不同量级的噪声激励和振动激励的联合激励,通过传感器分别测得振动激励的激振力,振动激励的激振加速度和声激励的激振声压,以及通过加速度传感器测得响应,并记录相应的试验结果数据,用于载荷识别试验研究。优选的,一种基于一元线性回归模型和传递函数最小二乘广义逆的不相关多源频域载荷识别方法;该方法可根据测得的多个响应点的频域振动响应的自谱和多个载荷点到响应点传递函数模的平方的最小二乘广义逆,同时识别出多个不相关频域载荷源的大小;此最小二乘广义逆方法仅需测量载荷点到响应点传递函数模的平方和多个振动响应点的自谱,不需要测量相位,且能缓解矩阵求逆所出现的病态问题。优选的,根据施加的联合激励以及测得的响应,可进行不相关多源载荷识别的理论推导;步骤A1该系统有m个载荷激励输入fi(i=1,…,m),在该联合激励下,测得线性时不变系统的n个测点输出为yj(j=1,…,n)。根据叠加原理,线性系统的每一个输出都可以由各个分立输入所引起的响应叠加而成。其输入各激励之间的互功率谱密度矩阵Sff(ω)与输出各响应间的互功率谱密度矩阵Syy(ω)的关系为:(1)式中h(u)是系统的单位脉冲响应矩阵,Cff(τ)∈Rm×m是输入的协方差函数矩阵,是系统频率特性矩阵,是系统频率特性矩阵的共轭;式(1)给出了多输入/多输出情形下输出功率谱矩阵与输入功率谱矩阵之间的关系式;它显示了输入与输出功率谱关系的简明特点,正是频域分析法的优点所在;步骤A2在实际情况中,m与n不相等,因此要求取载荷谱矩阵,须对频响函数矩阵求广义逆,则在频域中的载荷识别公式可表示为:(2)式的主要问题是用试验获得系统的复频响应函数矩阵H(ω)的工作量太大,而用有限元法来获得H(ω)又存在仿真建模与试验的误差问题;步骤A3在m个输入载荷激励都是零均值的平稳随机过程,且在互不相关的情况下,m个输入载荷激励的协方差函数矩阵Cff(τ)∈Rm×m为对角阵,即:其对应的输入功率谱矩阵Sff(ω)也为对角阵此时,输出功率谱中主对角线上的任意一元素满足:(3)式写成矩阵后的形式为:其中,|Hj,i(ω)|2是输入fi对响应yj的传递函数模的平方,是待识别的载荷源fi的自功率谱,是响应yj的自功率谱;步骤A4记(4)式可简写为:1)当n>m,(4)式为超定方程,无对应的满足(2)式的解。它的最小二乘解为:2)当n=m,(4)式为正定方程,对应的满足(4)式的解唯一,其解为:3)当n<m,(4)式为欠定方程,对应,满足(4)式的解有无穷组;为保证反演出载荷激励的精度,(4)式中应满足n>m,并将该问题转化为一个优化问题,目标是找一组m个不相关平稳载荷激励使得系统的n个测点的响应能达到为验证该方法的正确性和精度,识别出来的激励可以与实际加载的激励进行比较;但是,(4)式本身是一个多目标优化问题,目标是找一组m个不相关平稳载荷激励使得在该组载荷激励作用下,系统的n个测点的响应与误差最小。在工程实践中,该问题需要转化成单目标优化问题,才能进行求解计算。步骤A5对于(4)式,当n≥m时,在响应误差平方和最小的单目标优化准则下的解的为:步骤A6证明:在(4)式中,响应误差平方和的一半为:为了使J最小化,以为参数,求J的梯度,可以得到(10)式:为了使J最小化,使(10)式最后结果等于零,从而得到以下等式:(11)式化简后得到最后的结果(8)式,证毕;(8)式又叫(4)式的最小二乘广义逆;步骤A5中,单优化目标准则为响应误差平方和最小。优选的,一种基于改进的正则化方法的不相关多源频域载荷识别方法;该方法可根据测得的多个响应点的频域振动响应和多个载荷点到响应点传递函数模的平方的正则化方法,同时识别出多个不相关频域源的大小;此改进的正则化方法,对应的每一个频率都存在最佳的正则化参数,从而解决了条件数较大时矩阵求逆的病态问题,提高了固有频率处多个不相关频域载荷源的识别精度。优选的,采用一种基于改进的正则化方法的不相关多源频域载荷识别方法,具体步骤如下:步骤B1根据权利2中的公式(4),当n≥m时,在正则化最小二乘代价函数下的解为:在(11)式中λ(ω)≥0称为正则化参数。步骤B2证明:由于正则化最小二乘代价函数为:为了使J最小化,以为参数,求J的梯度,可以得到(14)式:为了使J最小化,使(14)式最后结果等于零,从而得到以下等式:求解得到最后的结果(12)式,证毕;步骤B3对每一个频率选择最佳的正则化参数,选择正则化参数的原则是:使得所识别的多个载荷源的最大相对误差最小,公式(16)如下:所述步骤本文档来自技高网...
<a href="http://www.xjishu.com/zhuanli/55/CN105159865.html" title="复杂声振模拟实验环境下进行不相关多源频域载荷识别的装置和方法原文来自X技术">复杂声振模拟实验环境下进行不相关多源频域载荷识别的装置和方法</a>

【技术保护点】
一种不相关的声振载荷联合施加和不相关多源频域载荷识别研究的试验装置和方法,其特征在于,包括:两端加板密闭圆柱壳装置,通过弹性橡胶绳悬挂;所述圆柱壳装置内部有一个球型噪声激励源,由声传感器记录球型声源所施加的噪声激励大小,所述圆柱壳装置外部有一悬挂式振动台激振器,由振动传感器记录悬挂式振动台输入的振动激励;所述圆柱壳装置外表面和内表面布置有18个振动加速度传感器,用于测量机械结构在球型噪声源和悬挂式振动台联合激励下的振动响应;球型噪声激励源激励与悬挂式振动台激振器激励每次加载的时候位置和方向均固定不变,且18个振动加速度传感器分布在圆柱壳内部的各个地方,能反映两端加板密闭圆柱壳装置系统的主要振动。

【技术特征摘要】
1.一种声振载荷联合施加的试验装置,其特征在于,包括:两端加板密闭圆柱壳装置,通过弹性橡胶绳悬挂;所述圆柱壳装置内部有一个球型噪声激励源,由声传感器记录球型声源所施加的噪声激励大小,所述圆柱壳装置外部有一悬挂式振动台激振器,由振动传感器记录悬挂式振动台输入的振动激励;所述圆柱壳装置外表面和内表面布置有18个振动加速度传感器,用于测量机械结构在球型噪声源和悬挂式振动台联合激励下的振动响应;球型噪声激励源激励与悬挂式振动台激振器激励每次加载的时候位置和方向均固定不变,且18个振动加速度传感器分布在圆柱壳内部的各个地方,能反映两端加板密闭圆柱壳装置系统的主要振动。2.根据权利要求1所述的一种声振载荷联合施加的试验装置,其特征在于:所述独立的球型噪声激励源激励,有3种量级激励,而且量级逐渐增大;所述独立的悬挂式振动台激振器振动激励,有5种量级激励,而且量级逐渐增大;当噪声激励和振动激励联合加载时,噪声激励和振动激励互不相关,其载荷量级两两组合,形成了15种不同的量级;对声振试验装置分别加载15种不同量级的噪声激励和振动激励的联合激励,通过传感器分别测得振动激励的激振力,振动激励的激振加速度和声激励的激振声压,以及通过18个加速度传感器测得圆柱壳装置外表面和内表面的振动响应,并记录相应的试验结果数据,从而实现了模拟复杂的声振环境,用于载荷识别试验研究。3.一种基于一元线性回归模型和传递函数最小二乘广义逆的不相关多源频域载荷识别方法,采用利用权利要求1和2所述的声振载荷联合施加的试验装置;其特征在于:该方法可根据测得的多个响应点的频域振动响应的自谱和多个载荷点到响应点传递函数模的平方的最小二乘广义逆,同时识别出多个不相关频域载荷源的大小;此最小二乘广义逆方法仅需测量载荷点到响应点传递函数模的平方和多个振动响应点的自谱,不需要测量相位,且能缓解矩阵求逆所出现的病态问题。4.根据权利要求3所述的一种基于一元线性回归模型和传递函数最小二乘广义逆的不相关多源频域载荷识别方法,其特征在于:根据施加的联合激励以及测得的响应,可进行不相关多源载荷识别的理论推导;步骤A1该系统有m个载荷激励输入fi(i=1,…,m),在该联合激励下,测得线性时不变系统的n个测点输出为yj(j=1,…,n);根据叠加原理,线性系统的每一个输出都可以由各个分立输入所引起的响应叠加而成,其输入各激励之间的互功率谱密度矩阵Sff(ω)与输出各响应间的互功率谱密度矩阵Syy(ω)的关系为:(1)式中h(u)是系统的单位脉冲响应矩阵,Cff(τ)∈Rm×m是输入的协方差函数矩阵,是系统频率特性矩阵,是系统频率特性矩阵的共轭;式(1)给出了多输入/多输出情形下输出功率谱矩阵与输入功率谱矩阵之间的关系式;它显示了输入与输出功率谱关系的简明特点,正是频域分析法的优点所在;步骤A2在实际情况中,m与n不相等,因此要求取载荷谱矩阵,须对频响函数矩阵求广义逆,则在频域中的载荷识别公式可表示为:(2)式的主要问题是用试验获得系统的复频响应函数矩阵H(ω)的工作量太大,而用有限元法来获得H(ω)又存在仿真建模与试验的误差问题;步骤A3在m个输入载荷激励都是零均值的平稳随机过程,且在互不相关的情况下,m个输入载荷激励的协方差函数矩阵Cff(τ)∈Rm×m为对角阵,即:其对应的输入功率谱矩阵Sff(ω)也为对角阵此时,输出功率谱中主对角线上的任意一元素满足:(3)式写成矩阵后的形式为:其中,|Hj,i(ω)|2是输入fi对响应yj的传递函数模的平方,是待识别的载荷源fi的自功率谱,是响应yj的自功率谱;步骤A4记(4)式可简写为:1)当n>m,(4)式为超定方程,无对应的满足(2)式的解,它的最小二乘解为:2)当n=m,(4)式为正定方程,对应的满足(4)式的解唯一,其解为:3)当n<m,(4)式为欠定方程,对应,满足(4)式的解有无穷组;为保证反演出载荷激励的精度,(4)式中应满足n>m,并将该问题转化为一个优化问题,目标是找一组m个不相关平稳载荷激励使得系统的n个测点的响应能达到为验证该方法的正确性和精度,识别出来的激励可以与实际加载的激励进行比较;但是,(4)式本身是一个多目标优化问题,目标是找一组m个不相关平稳载荷激励使得在该组载荷激励作用下,系统的n个测点的响应与误差最小,在工程实践中,该问题需要转化成单目标优化问题,才能进行求解计算;步骤A5对于(4)式,当n≥m时,在响应误差平方和最小的单目标优化准则下的解的为:步骤A6证明:在(4)式中,响应误差平方和的一半为:为了使J最小化,以为参数,求J的梯度,可以得到(10)式:为了使J最小化,使(10)式最后结果等于零,从而得到以下等式:(11)式化简后得到最后的结果(8)式,证毕;(8)式又叫(4)式的最小二乘广义逆;步骤A5中,单优化目标准则为响应误差平方和最小。5.一种基于改进的正则化方法的不相关多源频域载荷识别方法,利用权利要求1和2所述的声振载荷联合施加的试验装置,;其特征在于:该方法可根据测得的多个响应点的频域振动响应和多个载荷点到响应点传递函数模的平方的正则化方法,同时识别出多个不相关频域源的大小;此改进的正则化方法,对应的每一个频率都存在最佳的正则化参数,从而解决了条件数较大时矩阵求逆的病态问题,提高了固有频率处多个不相关频域载荷源的识别精度。6.根据权利要求5所述的一种基于改进的正则化方法的不相关多源频域载荷识别方法,其特征在于:采用一种基于改进的正则化方法的不相关多源频域载荷识别方法,具体步骤如下:步骤B1根据权利4中的公式(4),当n≥m时,在正则化最小二乘代价函数下的解为:在(12)式中λ(ω)≥0称为正则化参数;步骤B2证明:由于正则化最小二乘代价函数为:为了使J最小化,以为参数,求J的梯度,可以得到(14)式:

【专利技术属性】
技术研发人员:王成王建英官威赖雄鸣钟必能张惠臻蔡奕侨李静彭佳林应晖郑黎晓
申请(专利权)人:华侨大学
类型:发明
国别省市:福建;35

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1