一种获取耦合损耗因子的数值方法技术

本发明专利技术提出一种获取耦合损耗因子的数值方法，该方法包括如下步骤：(1)将系统划分成连续耦合的子系统；(2)计算子系统的模态；(3)计算相邻子系统中模态间的耦合参数；(4)基于模态能量法计算模态上载荷输入功率及模态振动能量在频域内的分布；(5)计算分析频带内子系统上载荷的总输入功率及子系统的总振动能量；(6)基于功率输入法计算耦合损耗因子。本发明专利技术提供的获取统计能量法理论中的耦合损耗因子的方法，是一种基于模态能量法和功率输入法的获取耦合损耗因子的数值方法，该方法把模态能量法与功率输入法结合，可考虑非共振模态间的功率传输的影响，能够精确计算统计能量法中子系统间的耦合损耗因子。

A numerical method for obtaining the coupling loss factor

The invention provides a method for obtaining the coupling loss factor numerical method, the method comprises the following steps: (1) divide the system into subsystems for coupling; (2) calculating modal system; (3) calculating modal coupling parameters of adjacent sub system; (4) calculating modal load input power and modal vibration energy in the frequency domain distribution mode based on energy method; (5) to calculate the total vibration energy analysis of the total input power and frequency subsystem subsystem load; (6) calculating power input method based on coupling loss factor. Method to obtain the coupling loss factor in the theory of statistical energy analysis method provided by the invention is a numerical method of modal energy method and power input method to obtain the coupling loss factor based on the modal energy method combined with power input method, can consider the impact of power transmission between the non resonant mode, can accurately calculate the coupling loss factor between subsystems in the statistical energy analysis method.

【技术实现步骤摘要】
一种获取耦合损耗因子的数值方法
本专利技术涉及基于统计能量法的高频动响应预示
，具体涉及一种获取耦合损耗因子的数值方法。
技术介绍
航天器在任务周期内面临严峻的高强、宽频噪声等环境，这可能造成结构失效或精密仪器、仪表失灵。因此，在航天器的设计过程中，需考虑高强、宽带噪声的影响。可采用试验方法、理论方法和数值方法预示系统在高频噪声激励下的动响应。其中，试验方法能得到可靠的结果，但开展试验分析的成本较高，设计周期长；理论方法只适用于简单系统，难以解决复杂系统的动响应预示问题；数值方法对复杂系统有良好的适用性，是试验分析的有效辅助手段。统计能量法能有效预示噪声环境下系统的高频动响应，得到广泛应用。基于统计能量法解决工程问题的工作量集中体现在构建统计能量法基本方程，包括子系统的划分、子系统耦合关系的定义以及统计能量分析参数的获取。子系统的划分和子系统耦合关系的定义比较依赖工程人员的经验；统计能量分析参数的获取是统计能量法整个分析过程中极为重要的一环。耦合损耗因子是统计能量分析参数之一，准确获取耦合损耗因子是基于统计能量法开展高频动响应预示的必要条件。
技术实现思路
专利技术目的：针对高频动响应预示方法——统计能量法中存在的问题，本专利技术公开了获取耦合损耗因子的一种数值方法，该方法把模态能量法与功率输入法结合，可考虑非共振模态间的功率传输的影响，能够精确计算统计能量法中子系统间的耦合损耗因子。技术方案：为实现上述技术效果，本专利技术提供的技术方案为：一种获取耦合损耗因子的数值方法，包括步骤：(1)将系统划分成连续耦合的子系统，并分别用位移模态和应力/声压模态描述耦合界面上相邻的两个子系统的振动特性；(2)获取各子系统的模态参数；(3)计算相邻子系统间的模态耦合参数；(4)根据得到的各子系统的模态参数和相邻子系统间的模态耦合参数，基于模态能量法构建各模态上的功率平衡方程，并根据功率平衡方程得到任意模态上载荷输入功率与模态振动能量在频域内的分布关系；(5)向子系统施加激励，计算分析频带内子系统上载荷的总输入功率及子系统的总振动能量；(6)建立统计能量法系数矩阵,并根据分析频带内子系统上载荷的总输入功率及子系统的总振动能量计算子系统间的耦合损耗因子；统计能量法系数矩阵为：其中，PN×N＝diag(Ptk)式中，η为统计能量法系数矩阵，ηk为子系统k的内损耗因子，ηkl为子系统k对子系统l的耦合损耗因子；N为子系统的个数，ω0为分析频带Δω的中心频率，上标“-1”表示矩阵的逆矩阵；PN×N为系统输入功率矩阵，diag()表示对角矩阵，括号内为对角矩阵元素，Ptk为子系统k内所有模态受到单位激励时分析频带Δω内子系统k上载荷的总输入功率；EN×N为系统能量响应矩阵，为只有子系统k内所有模态受到单位激励时，分析频带Δω内子系统l的总振动能量。进一步的，所述相邻子系统间的模态耦合参数的计算公式为：式中，γmn为位移模态m与应力模态n之间的耦合参数，为位移模态m的第i个分量，s为空间位置，为应力模态n的应力张量，nj为应力模态n所在子系统在耦合边界上单位外法向向量的第j个分量，Sc为相邻子系统的耦合界面，Mm和Mn分别为模态m和模态n的模态质量。进一步的，所述功率平衡方程为：式中，Piinj(ω)为模态i上的载荷输入功率，ω为角频率；Pidiss(ω)为模态i的阻尼损耗功率；Ei(ω)和Ej(ω)分别为模态i和模态j的单频振动能量；αij(ω)和αji(ω)为单频耦合损耗系数；Piinj(ω)、Pidiss(ω)、αij(ω)和αji(ω)的表达式分别为：式中，Sii(ω)为模态i上载荷的自功率谱，Re(·)表示虚数的实部，Hii(ω)为频响函数，其含义为：模态i上受到单位力激励时模态i的位移响应，λi和ωi分别为模态i的模态阻尼和模态固有频率，λj和ωj分别为模态j的模态阻尼和模态固有频率；γij为模态i和模态j之间的耦合参数。进一步的，所述分析频带Δω内子系统上载荷的总输入功率及子系统的总振动能量的计算公式为：式中，Ptk表示分析频带Δω内子系统k上载荷的总输入功率，表示只有子系统k受到激励时，分析频带Δω内子系统l的总振动能量。进一步的，所述各子系统的模态参数基于有限元法计算得到。有益效果：与现有技术相比，本专利技术具有以下优势：本专利技术把模态能量法与功率输入法结合，可考虑非共振模态间的功率传输的影响，能够精确计算统计能量法中子系统间的耦合损耗因子。附图说明图1是本专利技术的原理流程框图；图2是一个平板/声腔耦合系统的示意图；图3是平板的有限元模型；图4是声腔的有限元模型；图5是在不同中心频率分析频带内的平板对声腔的耦合损耗因子；图6是在不同中心频率分析频带内的声腔对平板的耦合损耗因子；图7是采用本专利技术计算得到的平板的振动能量和平板振动能量的参考值在不同频率上的对比图；图8是采用本专利技术计算得到声腔的振动能量和声腔振动能量的参考值在不同频率上的对比图。具体实施方式下面结合附图对本专利技术作更进一步的说明。如图1所示为本专利技术的原理流程框图，主要包括以下步骤：步骤(1)将系统划分成连续耦合的子系统：用虚构的界面将系统划分成连续耦合的子系统，并对子系统在耦合界面上的边界条件进行近似，假定较“刚”子系统在耦合界面上为自由边界，假定较“柔”子系统在耦合界面上受到固定支撑。分别用位移模态和应力/声压模态描述相邻子系统的振动。步骤(2)计算子系统的模态：基于有限元法计算了子系统的模态参数。步骤(3)计算相邻子系统中模态间的耦合参数：其中，γmn为位移模态m与应力模态n之间的耦合参数，为位移模态m的第i个分量，s为空间位置，为应力模态n的应力张量，nj为应力模态n所在子系统在耦合边界上单位外法向向量的第j个分量，Sc为相邻子系统的耦合界面，Mm和Mn分别为模态m和模态n的模态质量。步骤(4)基于模态能量法计算模态上载荷输入功率及模态振动能量在频域内的分布：(4.1)基于下式计算模态上载荷输入功率：其中，ω为角频率，Sii(ω)为模态i上载荷的自功率谱，Re(·)表示虚数的实部，Hii(ω)为频响函数，其含义为：模态i上受到单位力激励时模态i的位移响应。(4.2)基于下式计算模态阻尼损耗功率：其中λi和ωi分别为模态i的模态阻尼和模态固有频率，Ei(ω)为模态i的单频振动能量。(4.3)基于下式计算模态间单频耦合损耗系数：其中γij为模态i和模态j之间的耦合参数，由式(1)计算得到，λj和ωj分别为模态j的模态阻尼和模态固有频率。将式(4)中下标i和j对调后可得到αji(ω)的表达式。步骤(4.4)基于下式计算单频模态振动能量：其中Ej(ω)为模态j的单频振动能量。步骤(5)计算分析频带内子系统上载荷的总输入功率及子系统的总振动能量：步骤(5.1)基于下式计算分析频带Δω内子系统k上载荷的总输入功率：Ptk＝∫Δω(ΣiPiinj(ω))dω(6)步骤(5.2)基于下式计算只有子系统k受到激励时，分析频带Δω内子系统l的总振动能量：步骤(6)基于功率输入法计算耦合损耗因子；步骤(6.1)统计能量法系数矩阵为其中ηk为子系统k的内损耗因子，ηkl为子系统k对子系统l的耦合损耗因子。步骤(6.2)基于下式计算统计能量法系数矩阵：其中N为子系统的个数，ω0为本文档来自技高网...

【技术保护点】
一种获取耦合损耗因子的数值方法，其特征在于，包括步骤：(1)将系统划分成连续耦合的子系统，并分别用位移模态和应力模态或位移模态和声压模态描述耦合界面上相邻的两个子系统的振动特性；(2)获取各子系统的模态参数；(3)计算相邻子系统间的模态耦合参数；(4)根据得到的各子系统的模态参数和相邻子系统间的模态耦合参数，基于模态能量法构建各模态上的功率平衡方程，并根据功率平衡方程得到任意模态上载荷输入功率与模态振动能量在频域内的分布；(5)向子系统施加激励，计算分析频带内子系统上载荷的总输入功率及子系统的总振动能量；(6)建立统计能量法系数矩阵,并根据分析频带内子系统上载荷的总输入功率及子系统的总振动能量计算子系统间的耦合损耗因子；统计能量法系数矩阵为：

【技术特征摘要】
1.一种获取耦合损耗因子的数值方法，其特征在于，包括步骤：(1)将系统划分成连续耦合的子系统，并分别用位移模态和应力模态或位移模态和声压模态描述耦合界面上相邻的两个子系统的振动特性；(2)获取各子系统的模态参数；(3)计算相邻子系统间的模态耦合参数；(4)根据得到的各子系统的模态参数和相邻子系统间的模态耦合参数，基于模态能量法构建各模态上的功率平衡方程，并根据功率平衡方程得到任意模态上载荷输入功率与模态振动能量在频域内的分布；(5)向子系统施加激励，计算分析频带内子系统上载荷的总输入功率及子系统的总振动能量；(6)建立统计能量法系数矩阵,并根据分析频带内子系统上载荷的总输入功率及子系统的总振动能量计算子系统间的耦合损耗因子；统计能量法系数矩阵为：其中，PN×N＝diag(Ptk)式中，η为统计能量法系数矩阵，ηk为子系统k的内损耗因子，ηkl为子系统k对子系统l的耦合损耗因子；N为子系统的个数，ω0为分析频带Δω的中心频率，上标“-1”表示矩阵的逆矩阵；PN×N为系统输入功率矩阵，diag()表示对角矩阵，括号内为对角矩阵元素，Ptk为子系统k内所有模态受到单位激励时分析频带Δω内子系统k上载荷的总输入功率；EN×N为系统能量响应矩阵，为只有子系统k内所有模态受到单位激励时，分析频带Δω内子系统l的总振动能量。2.根据权利要求1所述的一种获取耦合损耗因子的数值方法，其特征在于，所述相邻子系统间的模态耦合参数的计算公式为：1式中，γmn为位移模态m与应力模态n之间的耦合参数，为位移模态m的第i个分量，s为空间位置，为应力模态n的应力张量，nj为应力模态n所在子系统在耦合边界上单位外法向向量的第j个分量，Sc为相邻子系统的耦合界面，Mm和Mn分别为模态m和模态n的模态质量。3.根据权利要求2所述的一种获取耦合损耗因子的数值方法，其特征在于，所述功率平衡方程为：式中，为模态i上的载荷输入功率...

【专利技术属性】
技术研发人员：李彦斌，张鹏，费庆国，吴邵庆，杨轩，姜东，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏,32