CN108875243A - 一种基于能量灵敏度分析的声振模型修正方法 - Google Patents
一种基于能量灵敏度分析的声振模型修正方法 Download PDFInfo
- Publication number
- CN108875243A CN108875243A CN201810697611.4A CN201810697611A CN108875243A CN 108875243 A CN108875243 A CN 108875243A CN 201810697611 A CN201810697611 A CN 201810697611A CN 108875243 A CN108875243 A CN 108875243A
- Authority
- CN
- China
- Prior art keywords
- subsystem
- energy
- parameter
- model
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
一种基于能量灵敏度分析的声振模型修正方法,通过分析和试验结果残差建立目标函数,以能量灵敏度矩阵确定模型参数的变化量,并通过迭代求解参数修正值使得目标函数达到极小,本发明根据能量灵敏度值对统计能量分析模型中耦合损耗因子或内损耗因子等待修正参数进行科学选择,并可灵活选取实验测试子系统响应与缩减子系统分析值之间的残差目标函数,采用数值稳定性高的正则化法迭代求解模型参数修正量,同时提供了一套系统完整的误差灵敏度分析流程,用于估计修正参数的误差区间。
Description
技术领域
本发明涉及一种基于能量灵敏度分析的声振模型修正方法,属于飞行器设计技术领域。
背景技术
航天器发射过程中,结构会因火箭喷气噪声和气动噪声激励产生恶劣声振环境,会对航天器结构及其组件造成严重的破坏。而且随着航天器结构轻质化、柔性化发展,也增加了结构对声振环境,尤其是高频激励的敏感性。为提高航天器的研制效率,保证结构及其设备的安全可靠,对创建结构声振响应分析模型提出了更高的要求。
有限元分析和统计能量分析目前航天器结构声振响应分析应用较广泛的两种方法。一般来说有限元分析主要适用于低频区结构动力学分析,但是在高频区域,由于结构振动波长短,计算量大,模态密集,且高频模态对结构细节敏感,表现出随机特性,更适合采用基于统计理论的统计能量分析方法(SEA)。SEA是目前被公认解决复杂系统宽带高频声振问题的有力的工具,并已经在航天领域的结构分析中得到了广泛应用。
采用统计能量分析法进行结构响应预示,首先需要将工程结构理想化并建立合理的统计能量分析模型。其关键的环节是将复杂结构系统划分为满足统计能量分析假设的多个子系统,并合理地估计各子系统的模态密度和内损耗因子以及子系统间的耦合损耗因子等SEA参数。经验表明,而且在估计SEA参数时引起的误差也会对最后振动响应的预示结果造成极大的误差。目前统计能量分析参数的理论计算公式大多只是适用于一些理想结构。对于复杂的结构系统,通常需要借助实验方法甚至根据经验确定比较合理的统计能量分析模型。近年来发展起来的以功率输入法为代表的实验统计能量分析方法丰富和方便了统计能量分析的建模与分析,但是由于实验过程耗时步骤繁琐以及计算的数值稳定性和精确性差,极大地限制了其工程的应用。
发明内容
本发明解决的技术问题是:针对目前现有技术中航天器结构的声振环境预示及高频动力学分析时模型参数估计困难或不准确的问题,提出一种基于能量灵敏度分析的声振模型修正方法。
本发明解决上述技术问题是通过如下技术方案予以实现的:
一种基于能量灵敏度分析的声振模型修正方法,具体步骤如下:
(1)开始飞行器声振试验并记录试验中飞行器各子系统能量的测量值,同时建立飞行器声振分析模型并确定所有模型参数,并模拟计算出飞行器各子系统能量的预示值;
(2)根据步骤(1)所得数据获取飞行器各子系统的能量残差,并计算各子系统能量残差的加权范数;
(3)计算各子系统模型参数的能量灵敏度,根据计算结果选取待修正的模型参数;
(4)采用迭代方式计算步骤(3)所得待修正参数的修正值;
(5)对步骤(4)所得修正值进行病态问题判断,若存在病态问题,进行正则化法修正并进入步骤(6),若不存在,则直接进入步骤(6);
(6)根据所得修正值重复进行步骤(2),判断计算所得子系统残差的加权范数是否满足收敛条件,若满足收敛条件则将该修正值作为最终修正值,若不满足收敛条件,则返回步骤(4)继续进行修正值迭代计算。
所述步骤(2)中,各子系统能量残差加权范数的计算方法如下:
其中,能量残差函数 表示第m个子系统能量的测量值,Em为第m个子系统能量的预示值,为待修正的模型参数组成的列向量,W为加权矩阵,其中T表示矩阵的转置。
所述步骤(3)中,所述能量灵敏度的计算公式为:
式中,为第m个子系统对应的模型参数的能量灵敏度,为第m个子系统能量关于第r个模型参数的能量灵敏度矩阵偏导数,r=1,2,3,…,β,其中β为第m个子系统中模型参数总个数。
所述第m个系统能量矩阵关于对应的第r个模型参数的能量灵敏度矩阵的计算方法为:
H为损耗因子矩阵,为第r个模型参数、Es为第s个模型子系统能量预示值,Πm、Πs分别为Em、Es对应子系统的输入功率,Hmm、Hms、Hss、Hsm分别为损耗因子矩阵H第m行第s列分割后由左上角沿顺时针方向所得四个分块矩阵;
其中,H矩阵的组成如下
式中,ηij为子系统i到子系统j的耦合损耗因子参数,其中i≤j,ηi为子系统i的内损耗因子参数,ni为子系统i的模态密度参数。
所述步骤(3)中,选取待修正的参数模型的依据为:
对于任一子系统,从步骤(3)所得该子系统全部模型参数的能量灵敏度由大到小排列顺序中选取前N个能量灵敏度值对应的模型参数作为待修正的模型参数,其中N为飞行器声振试验中该子系统包含的测量装置的数量。
所述步骤(4)中,所述待修正参数的修正值计算公式为:
式中,i表示迭代步数,i=1,2,3,…,为步骤(3)选取的待修正参数,为修正后参数,为关于模型参数的能量灵敏度,为关于模型参数的子系统能量残差。
所述步骤(5)中,进行正则化法修正的计算公式为:
式中,λ为正则化修正参数,WP表示参数增量的加权矩阵,其中:
WP=tr(D)D-1/tr(D-1)
式中,矩阵D是由矩阵对角元素组成的对角矩阵,tr()表示矩阵的迹。
所述步骤(6)中,所述模型参数修正值收敛条件判断方法为:
式中,为模型参数最终修正值,为第i+1个子系统的子系统能量残差,为第i+1个子系统的子系统残差的加权范数;
其中,当时,判定收敛,此时得到模型参数最终修正值;当时,返回步骤(4)继续进行待修正参数迭代计算;其中,ε为正常数,取值范围为1×10-2~1×10-6。
本发明与现有技术相比的优点在于:
(1)本发明提供的一种基于能量灵敏度分析的声振模型修正方法,针对航天器结构在发射过程中会因气动噪声激励产生恶劣声振且现有声振环境预示及高频动力学分析模型在估计SEA参数时容易对振动响应的预示结果容易造成极大的误差的问题,通过对所有待修正参数进行能量灵敏度计算,并选择满足条件的模型参数进行修正,减小了模型分析过程中的计算误差,可根据灵敏度值对统计能量分析模型中耦合损耗因子或内损耗因子等待修正参数进行科学选择,避免传统方法的盲目性,提高模型修正效率;
(2)本发明克服了传统试验方法对对激振和测量子系统的完备性限制,便于工程复杂航天器结构声振分析模型修正的应用,同时可灵活选取实验测试子系统响应与对应分析值之间的残差为目标函数,采用迭代方法求解模型参数修正量,物理意义更明确,同时基于正则化法解决了灵敏度矩阵的病态问题,算法的数值稳定性高;
(3)本发明提供了一套系统完整的误差灵敏度分析流程,用于估计修正参数的误差区间,可操作范围更大,同时简化了计算步骤及试验流程。
附图说明
图1为发明提供的修正方法流程图;
具体实施方式
一种基于能量灵敏度分析的声振模型修正方法,如图1所示,具体步骤如下:
(1)建立声振分析模型,开始飞行器声振能量分析试验并记录加速度测量数据,同时将所得加速度测量数据转化为各子系统输入功率及能量值;
(2)根据步骤(1)所得试验处理数据及子系统能量预示值获取与试验测量值残差,并计算子系统残差的加权范数,其中:
其中,能量残差函数 表示第m个子系统能量的试验测量值,Em表示由分析模型计算的第m个测量子系统能量预示值,表示待修正的模型参数组成的列向量,W为贡献加权矩阵,可以用来控制各个测量子系统能量残差对目标函数的贡献,合理的模型参数修正值应使得能量残差构成的目标函数极小。加权矩阵W根据实际需求选择,取使得子系统能量残差的相对误差向量范数最小;
(3)根据声振分析模型初始参数进行能量灵敏度分析,通过统计能量分析的稳态功率流平衡方程对所有待修正的模型参数的偏导数获取子系统分析模型的能量灵敏度矩阵,根据计算结果选取待修正的模型参数,其中:
所述步骤(3)中,所述能量灵敏度矩阵的计算公式为:
式中,为各子系统对应的模型参数的能量灵敏度,为第m个子系统能量矩阵关于第r个模型参数的能量灵敏度矩阵偏导数,r=1,2,3,…,β,β为第m个子系统的模型参数个数;
其中,所述第m个系统能量矩阵关于第r个模型参数的能量灵敏度矩阵的计算方法为:
H为损耗因子矩阵,为第r个模型参数、Es为第s个模型子系统能量预示值,Πm、Πs分别为Em、Es对应子系统的输入功率,Hmm、Hms、Hsm、Hss分别为损耗因子矩阵H第m行s列分割后由左上角沿顺时针方向所得四个分块矩阵;
其中,H矩阵的组成如下
式中,ηij为子系统i到子系统j的耦合损耗因子参数,其中i≤j,ηi为子系统i的内损耗因子参数,ni为子系统i的模态密度参数;
同时,选取待修正的参数模型的依据为:
对于任一子系统,从步骤(3)所得该子系统全部模型参数的能量灵敏度由大到小排列顺序中选取前N个能量灵敏度值对应的模型参数作为待修正的模型参数,其中N为飞行器声振试验中该子系统包含的测量装置的数量;
(4)根据修正模型计算公式迭代计算步骤(3)所得待修正参数的修正值;
所述步骤(4)中,所述待修正参数修正值迭代计算公式为:
式中,i表示迭代步数,i=1,2,3,…,为步骤(3)选取的待修正参数,为修正后参数,为关于模型参数的能量灵敏度,为关于模型参数的子系统能量残差;
(5)进行能量灵敏度矩阵病态问题判断及模型参数修正值收敛条件判断,若条件数大于等于100,则视为存在病态问题,需要进行正则化法修正,若条件数小于100,则为良态,不存在病态问题,进入步骤(6)进行修正值收敛条件判断;
(6)若步骤(2)计算所得子系统残差的加权范数满足收敛条件,则迭代所得修正值作为最终修正值,若不满足收敛条件,则返回步骤(4)继续进行修正值迭代计算,其中:
此时,判断能量灵敏度矩阵是否存在病态问题,若存在,进行正则化法修正,若不存在,继续进行迭代计算;
所述待修正参数修正值的正则化法迭代计算公式为:
式中,λ为正则化参数,由残差范数和修正参数增量范数在双对数坐标系下的L曲线得到,WP表示参数增量的加权矩阵,加权矩阵的计算方法如下:
WP=tr(D)D-1/tr(D-1)
式中,矩阵D是由矩阵对角元素组成的对角矩阵,tr()表示矩阵的迹;
所述模型参数修正值收敛条件判断方法为:
式中,表示模型参数最终修正值,为第i+1个子系统的子系统能量残差,为第i+1个子系统的子系统残差的加权范数;
当时,判定收敛,此时得到模型参数最终修正值;当时,则返回步骤(4)继续进行待修正参数迭代计算;其中,ε为比较小的正数,对于子系统能量残差的相对误差向量范数,取值范围一般为1×10-2~1×10-6;
(7)根据未修正模型参数、步骤(1)所得的各子系统输入功率及能量对应的误差灵敏度分析公式估计最终修正参数的误差区间;
所述未修正模型参数的误差灵敏度分析公式为:
式中,为α个待修正模型参数列向量,为(β-α)个未修正模型参数,E为子系统能量列向量,上标“T”表示矩阵的转置,表示矩阵的Moore-Penrose广义逆。
所述各子系统输入功率的误差灵敏度分析公式为:
式中,Π为子系统的输入功率列向量。
所述各子系统能量的误差灵敏度分析公式为:
根据未修正模型参数、步骤(1)所得的各子系统输入功率及能量的方差及误差灵敏度进行修正参数的误差区间估计,计算公式如下:
式中,j=1,2,…,α,为修正参数估计误差的标准方差,为修正参数估计误差的标准方差,为子系统能量测量误差的标准方差,为子系统能量测量误差的标准方差,m为测量子系统数。
Claims (8)
1.一种基于能量灵敏度分析的声振模型修正方法,其特征在于步骤如下:
(1)开始飞行器声振试验并记录试验中飞行器各子系统能量的测量值,同时建立飞行器声振分析模型并确定所有模型参数,并模拟计算出飞行器各子系统能量的预示值;
(2)根据步骤(1)所得数据获取飞行器各子系统的能量残差,并计算各子系统能量残差的加权范数;
(3)计算各子系统模型参数的能量灵敏度,根据计算结果选取待修正的模型参数;
(4)采用迭代方式计算步骤(3)所得待修正参数的修正值;
(5)对步骤(4)所得修正值进行病态问题判断,若存在病态问题,进行正则化法修正并进入步骤(6),若不存在,则直接进入步骤(6);
(6)根据所得修正值重复进行步骤(2),判断计算所得子系统残差的加权范数是否满足收敛条件,若满足收敛条件则将该修正值作为最终修正值,若不满足收敛条件,则返回步骤(4)继续进行修正值迭代计算。
2.根据权利要求1所述的一种基于能量灵敏度分析的声振模型修正方法,其特征在于:所述步骤(2)中,各子系统能量残差加权范数的计算方法如下:
其中,能量残差函数 表示第m个子系统能量的测量值,Em为第m个子系统能量的预示值,为待修正的模型参数组成的列向量,W为加权矩阵,其中T表示矩阵的转置。
3.根据权利要求2所述的一种基于能量灵敏度分析的声振模型修正方法,其特征在于:所述步骤(3)中,所述能量灵敏度的计算公式为:
式中,为第m个子系统对应的模型参数的能量灵敏度,为第m个子系统能量关于第r个模型参数的能量灵敏度矩阵偏导数,r=1,2,3,…,β,其中β为第m个子系统中模型参数总个数。
4.根据权利要求3所述的一种基于能量灵敏度分析的声振模型修正方法,其特征在于:所述第m个系统能量矩阵关于对应的第r个模型参数的能量灵敏度矩阵的计算方法为:
H为损耗因子矩阵,为第r个模型参数、Es为第s个模型子系统能量预示值,Πm、Πs分别为Em、Es对应子系统的输入功率,Hmm、Hms、Hss、Hsm分别为损耗因子矩阵H第m行第s列分割后由左上角沿顺时针方向所得四个分块矩阵;
其中,H矩阵的组成如下
式中,ηij为子系统i到子系统j的耦合损耗因子参数,其中i≤j,ηi为子系统i的内损耗因子参数,ni为子系统i的模态密度参数。
5.根据权利要求3所述的一种基于能量灵敏度分析的声振模型修正方法,其特征在于:所述步骤(3)中,选取待修正的参数模型的依据为:
对于任一子系统,从步骤(3)所得该子系统全部模型参数的能量灵敏度由大到小排列顺序中选取前N个能量灵敏度值对应的模型参数作为待修正的模型参数,其中N为飞行器声振试验中该子系统包含的测量装置的数量。
6.根据权利要求5所述的一种基于能量灵敏度分析的声振模型修正方法,其特征在于:所述步骤(4)中,所述待修正参数的修正值计算公式为:
式中,i表示迭代步数,i=1,2,3,…,为步骤(3)选取的待修正参数,为修正后参数,为关于模型参数的能量灵敏度,为关于模型参数的子系统能量残差。
7.根据权利要求6所述的一种基于能量灵敏度分析的声振模型修正方法,其特征在于:所述步骤(5)中,进行正则化法修正的计算公式为:
式中,λ为正则化修正参数,WP表示参数增量的加权矩阵,其中:
WP=tr(D)D-1/tr(D-1)
式中,矩阵D是由矩阵对角元素组成的对角矩阵,tr()表示矩阵的迹。
8.根据权利要求1所述的一种基于能量灵敏度分析的声振模型修正方法,其特征在于:所述步骤(6)中,所述模型参数修正值收敛条件判断方法为:
式中,为模型参数最终修正值,为第i+1个子系统的子系统能量残差,为第i+1个子系统的子系统残差的加权范数;
其中,当时,判定收敛,此时得到模型参数最终修正值;当时,返回步骤(4)继续进行待修正参数迭代计算;其中,ε为正常数,取值范围为1×10-2~1×10-6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810697611.4A CN108875243A (zh) | 2018-06-29 | 2018-06-29 | 一种基于能量灵敏度分析的声振模型修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810697611.4A CN108875243A (zh) | 2018-06-29 | 2018-06-29 | 一种基于能量灵敏度分析的声振模型修正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108875243A true CN108875243A (zh) | 2018-11-23 |
Family
ID=64297298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810697611.4A Pending CN108875243A (zh) | 2018-06-29 | 2018-06-29 | 一种基于能量灵敏度分析的声振模型修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108875243A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102411673A (zh) * | 2011-08-18 | 2012-04-11 | 西北工业大学 | 一种计算声振系统中高频动力学响应的方法 |
CN105184047A (zh) * | 2015-07-17 | 2015-12-23 | 西北工业大学 | 一种结构高频机械噪声在线预报方法 |
CN105760608A (zh) * | 2016-02-23 | 2016-07-13 | 北京空间飞行器总体设计部 | 一种整流罩声场填充系数的计算方法 |
CN105912854A (zh) * | 2016-04-11 | 2016-08-31 | 西安交通大学 | 一种机械结构动载荷识别的稀疏表征方法 |
-
2018
- 2018-06-29 CN CN201810697611.4A patent/CN108875243A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102411673A (zh) * | 2011-08-18 | 2012-04-11 | 西北工业大学 | 一种计算声振系统中高频动力学响应的方法 |
CN105184047A (zh) * | 2015-07-17 | 2015-12-23 | 西北工业大学 | 一种结构高频机械噪声在线预报方法 |
CN105760608A (zh) * | 2016-02-23 | 2016-07-13 | 北京空间飞行器总体设计部 | 一种整流罩声场填充系数的计算方法 |
CN105912854A (zh) * | 2016-04-11 | 2016-08-31 | 西安交通大学 | 一种机械结构动载荷识别的稀疏表征方法 |
Non-Patent Citations (2)
Title |
---|
向树红,王大钧,于丹 等: ""统计能量分析参数的综合确定方法及应用"", 《振动工程学报》 * |
张红亮: ""基于灵敏度分析的统计能量分析参数修正"", 《统计能量分析参数辨识及模型修正方法研究》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105843073B (zh) | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 | |
CN107085633B (zh) | 基于支持向量机的多点振动响应频域预测的装置及方法 | |
Chang et al. | A family of explicit algorithms for general pseudodynamic testing | |
CN115577436B (zh) | 一种求解不确定结构风致振动响应的组合深度学习方法 | |
Mayes et al. | Design studies for the transmission simulator method of experimental dynamic substructuring | |
CN109558635A (zh) | 一种基于单元模态应变能灵敏度的结构区间不确定性损伤识别方法 | |
CN110631792A (zh) | 基于卷积神经网络的抗震混合试验模型更新方法 | |
CN115221801B (zh) | 基于动态近似建模的飞行器不确定性传播分析方法和装置 | |
Timme et al. | Transonic aeroelastic instability searches using sampling and aerodynamic model hierarchy | |
CN101399708B (zh) | 一种建立网络性能模型的方法和设备 | |
CN101916242B (zh) | 非传统阻尼矩阵识别方法 | |
CN114139363A (zh) | 一种空间变异非平稳地震动时程的模拟方法 | |
Vendl et al. | Projection-based model order reduction for steady aerodynamics | |
CN109598371B (zh) | 一种飞行器载电子设备的温度预测方法及系统 | |
Arras et al. | Finite-element structural updating using frequency response functions | |
Bo et al. | A new method for wind-induced response analysis of long span roofs | |
CN108875243A (zh) | 一种基于能量灵敏度分析的声振模型修正方法 | |
Ayyub et al. | Methodology for the Development of Reliability‐Based Design Criteria for Surface Ship Structures | |
Morlier et al. | An EGO-like optimization framework for sensor placement optimization in modal analysis | |
CN109960776B (zh) | 一种用于水力走时和水力信号衰减反演计算的改进方法 | |
Liu et al. | Sensor placement for multi-input multi-output dynamic identification | |
Hajjout et al. | An h-adaptive element-free Galerkin meshless method using a posteriori error estimator | |
CN111428342A (zh) | 一种基于频域谱分解的随机动载荷识别方法 | |
Blais et al. | Goal-Oriented Adaptive Sampling Procedure for Projection-Based Reduced-Order Models for Aerodynamic Flows | |
Mayes | Model Correlation and Calibration. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181123 |
|
RJ01 | Rejection of invention patent application after publication |