[发明专利]一种基于对偶模态方程的确定性声固耦合响应预示方法

说明书

技术领域

本发明涉及声固耦合响应预示技术领域，具体涉及一种基于对偶模态方程的确定性声固耦合响应预示方法。

背景技术

航天器在任务周期内面临严峻的机械振动、噪声等环境，这可能造成结构失效或精密仪器、仪表失灵。因此，在航天器的设计过程中，需考虑机械振动和噪声的影响。可采用试验方法、理论方法和数值方法预示系统在确定性载荷激励下的声固耦合响应。其中，试验方法能得到可靠的结果，但开展试验分析的成本较高，设计周期长；理论方法只适用于简单系统，难以解决复杂系统声固耦合响应预示问题；数值方法对复杂系统有良好的适用性，是试验分析的有效辅助手段。对偶模态方程理论用虚构的界面将系统划分成耦合的子系统，并基于有限元计算子系统的模态，而不是整个耦合系统的模态，因此，对偶模态方程方法比传统的有限元法具有更高的分析效率。

现有对偶模态方程理论中存在模态截断问题，即需截取有限频率范围内的子系统模态参与响应预示，所选模态过少会引起误差，所选模态过多会造成计算资源浪费。因此，需要有一个准则界定模态截断的频率范围，以合理地基于对偶模态方程预示系统在确定性载荷激励下的声固耦合响应。

发明内容

发明目的：针对现有的一种对偶模态方程方法在应用中存在的问题，本发明公开了一种基于对偶模态方程的确定性声固耦合响应预示方法，该方法可有效提高确定性载荷激励下结构的声固耦合响应预示效率。

技术方案：为实现上述目的，本发明采用的技术方案为：

一种基于对偶模态方程的确定性声固耦合响应预示方法，包括以下步骤：

(1)将声固耦合系统中的结构和声腔划分成不同的子系统；

(2)计算结构子系统和声腔子系统的模态；

(3)计算相邻子系统中模态间的耦合参数；

(4)建立耦合系统的对偶模态方程；

(5)通过前置处理，获得确定性载荷作用下，子系统模态上受到的广义力载荷；

(6)计算对偶模态方程，获得所有模态的参与因子；

(7)通过模态叠加，计算系统确定性声固耦合响应。

进一步地，所述步骤(1)中的声固耦合系统为结构与声腔耦合系统，结构振动会影响声腔声压脉动，声腔声压脉动也会影响结构振动；结构子系统在与声腔子系统耦合界面上的边界条件被近似为自由状态，声腔子系统在与结构子系统耦合界面上的边界条件被近似为固定边界。

进一步地，所述步骤(2)中基于有限元法计算了结构子系统和声腔子系统的模态参数和模态振型。

进一步地，所述步骤(3)中相邻子系统中模态间的耦合参数由下式计算得到：

其中W_mn为结构子系统第m阶位移模态与声腔子系统第n阶声压模态之间的耦合参数，为结构子系统第m阶位移模态的振型，为声腔子系统第n阶声压模态的振型，S_c为结构子系统与声腔子系统之间的耦合界面，s为空间位置。

进一步地，所述步骤(4)中建立的耦合系统的对偶模态方程为：

其中M_m为结构子系统第m阶位移模态质量，M_n为声腔子系统第n阶声压模态质量，η_m和η_n分别为模态m和模态n的阻尼损耗系数，ω为角频率，φ_m(ω)为模态m的参与因子，为模态n的参与因子，F_m(ω)和F_n(ω)分别为模态m和模态n上受到的广义力载荷。

进一步地，所述步骤(5)中，结构子系统受到确定性的面压载荷P(s,ω)激励时，模态m上受到的广义力载荷由下式给出：

其中S_p为面压载荷作用面，当结构子系统受到确定性集中力载荷F₀(ω)激励时，模态m上受到的广义力载荷由下式给出：

其中s₀为集中力载荷F₀(ω)的作用位置。

进一步地，所述步骤(6)中将系统对偶模态方程写成了分块矩阵的形式：

其中上标“T”表示矩阵的转置，

W＝[W_mn] (21)

其中diag()表示对角矩阵，括号内为对角矩阵元素，W_mn为矩阵W的第m行第n列元素，其取值由式(1)计算得到，基于式(5)求得各子系统内的模态参与因子：