CN107748815A - 一种随机噪声环境下基于对偶模态方程的动响应分析方法 - Google Patents
一种随机噪声环境下基于对偶模态方程的动响应分析方法 Download PDFInfo
- Publication number
- CN107748815A CN107748815A CN201710958872.2A CN201710958872A CN107748815A CN 107748815 A CN107748815 A CN 107748815A CN 201710958872 A CN201710958872 A CN 201710958872A CN 107748815 A CN107748815 A CN 107748815A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mtr
- subsystem
- 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.)
- Granted
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Aviation & Aerospace Engineering (AREA)
- Operations Research (AREA)
- Automation & Control Theory (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- External Artificial Organs (AREA)
Abstract
本发明提出一种随机噪声环境下基于对偶模态方程的动响应分析方法,包括如下步骤:(1)将声固耦合系统中的结构和声腔划分成不同的子系统;(2)计算结构子系统和声腔子系统的模态;(3)计算相邻子系统中模态间的耦合参数;(4)建立耦合系统的对偶模态方程;(5)通过前置处理,获得随机载荷作用下,子系统模态上受到的广义力载荷的互功率谱;(6)计算对偶模态方程,获得所有模态的参与因子的互功率谱;(7)通过模态叠加,计算系统随机声固耦合响应。本发明提供的随机动响应分析方法,是一种基于对偶模态方程的随机噪声环境下动响应分析方法,该方法把系统划分成连续耦合的子系统,并用有限频带内的子系统模态描述系统的随机振动,该方法的分析效率高于传统有限元法。
Description
技术领域
本发明涉及声固耦合响应预示技术领域,尤其是一种随机噪声环境下基于对偶模态方程的动响应分析方法。
背景技术
随着航天器向高飞行速度发展,其在任务周期内面临严峻的随机噪声等环境,这可能造成结构失效或精密仪器、仪表失灵。因此,在航天器的设计过程中,需考虑随机噪声的影响。可采用试验方法、理论方法和数值方法预示系统在随机噪声激励下的动响应。其中,试验方法能得到可靠的结果,但开展试验分析的成本较高,设计周期长;理论方法只适用于简单系统,难以解决复杂系统的动响应预示问题;数值方法对复杂系统有良好的适用性,是试验分析的有效辅助手段。对偶模态方程理论用虚构的界面将系统划分成耦合的子系统,并基于有限元计算子系统的模态,而不是整个耦合系统的模态,因此,对偶模态方程方法比传统的有限元法具有更高的分析效率。
在基于对偶模态方程理论预示系统在随机噪声激励下的动响应时,需截取有限频率范围内的子系统模态参与响应预示,所选模态过少会引起误差,所选模态过多会造成计算资源浪费。因此,需要有一个准则界定模态截断的频率范围,以合理地基于对偶模态方程预示系统在随机噪声激励下的声固耦合响应。
发明内容
发明目的:为解决现有动响应分析技术中存在的技术问题,提供一个准则界定模态截断的频率范围,以合理地基于对偶模态方程预示系统在随机噪声激励下的声固耦合响应,本发明提出一种随机噪声环境下基于对偶模态方程的动响应分析方法。
技术方案:为实现上述技术效果,本发明提出的技术方案为:
一种随机噪声环境下基于对偶模态方程的动响应分析方法,该方法包括步骤:
(1)将声固耦合系统中的结构和声腔划分成在耦合界面上连续耦合的子系统,且耦合界面上相邻的两个子系统分别为声腔子系统和结构子系统;
(2)设置截断频率大于等于分析频率上限的1.25倍,截取结构子系统和声腔子系统中固有频率小于截断频率的模态;
(3)基于有限元法计算截取的各模态的模态参数,模态参数包括:模态质量、阻尼损耗系数和模态振型;
(4)根据各模态参数计算相邻子系统中截取出的模态之间的耦合参数;
(5)根据各子系统的模态参数和相邻子系统间的耦合参数建立相互耦合的两个相邻子系统的对偶模态方程为:
式中,ω为角频率,i表示虚数的虚部;Mm为结构子系统第m阶位移模态的模态质量;ωm为结构子系统第m阶位移模态的固有频率;ηm为结构子系统第m阶位移模态的阻尼损耗系数;φm(ω)为结构子系统第m阶位移模态的参与因子,Wmp为结构子系统第m阶位移模态与声腔子系统第p阶声压模态之间的耦合参数, 为声腔子系统第p阶声压模态的参与因子,Fm(ω)为结构子系统第m阶位移模态上受到的广义力载荷;
Mn为声腔子系统第n阶声压模态的模态质量;ωn为声腔子系统第n阶声压模态的固有频率;ηn为声腔子系统第n阶声压模态的阻尼损耗系数;为声腔子系统第n阶声压模态的参与因子;Wqn为结构子系统第q阶位移模态与声腔子系统第n阶声压模态之间的耦合参数;φq(ω)为结构子系统第q阶位移模态的参与因子;Fn(ω)为声腔子系统第n阶声压模态上受到的广义力载荷;
(5)将对偶模态方程转换为分块矩阵形式:
其中,
式中,上标“-1”表示矩阵的逆矩阵,上标“T”表示矩阵的转置;Hij为传递函数矩阵,i=1,2,j=1,2;矩阵元素Hij(k,l)表示当第j个子系统中第l阶模态上作用单位广义力时,第i个子系统中第k阶模态参与因子;各传递函数矩阵的计算公式为:
W(m,n)=Wmn
式中,diag()表示对角矩阵,括号内为对角矩阵元素;W(m,n)表示矩阵W中第m行第n列的元素,即结构子系统第m阶位移模态与声腔子系统第n阶声压模态之间的耦合参数Wmn;
(6)计算所述声固耦合系统中只有结构受到噪声激励时,分块矩阵S11和S22满足以下形式:
式中,为结构子系统模态载荷互功率谱矩阵,其第k行第l列元素为Skl(ω),Skl(ω)表示只有结构子系统受到随机噪声激励时,结构子系统第k阶位移模态上受到的广义力载荷和结构子系统第l阶位移模态上受到的广义力载荷之间的互谱,Skl(ω)的计算公式为:
式中,Ap为面压载荷作用面,为结构子系统第k阶位移模态的振型,为结构子系统第l阶位移模态的振型,Spp(s1,s2,ω)为面压载荷的功率谱,s1和s2为面压载荷作用面Ap上的空间位置;
(7)计算各结构子系统的位移响应和各声腔子系统的声压响应,其中,结构子系统的位移响应的计算公式为:
Sw(s,ω)表示第w个结构子系统在位置s处、角频率ω下的位移响应;
声腔子系统的声压响应的计算公式为:
式中,Sp(s,ω)表示第p个声腔子系统在位置s处、角频率ω下的声压响应。
进一步的,所述耦合参数的计算公式为:
式中,Wmn为结构子系统第m阶位移模态与声腔子系统第n阶声压模态之间的耦合参数,为结构子系统第m阶位移模态的振型,为声腔子系统第n阶声压模态的振型,Ac为结构子系统与声腔子系统之间的耦合界面,s为空间位置。
有益效果:与现有技术相比,本发明具有以下优势:
本发明是一种优于传统有限元法的随机噪声激励下的动响应预示方法,该方法可有效提高随机噪声激励下结构的动响应预示效率,缩短设计周期,节约设计成本。
附图说明
图1是本发明的逻辑流程框图;
图2是一个平板/声腔耦合系统的有限元模型;
图3是随机噪声激励下加筋板面板内各响应点处的加速度响应功率谱;
图4是随机噪声激励下声腔中各响应点处的声压响应功率谱。
具体实施方式
下面结合附图对本发明作更进一步的说明。
图1是本发明的逻辑流程框图,包括步骤:
步骤(1)将声固耦合系统中的结构和声腔划分成不同的子系统;声固耦合系统为结构与声腔耦合系统,结构振动与声场声压脉动之间存在交互作用;对子系统在耦合界面上的边界条件进行了近似,结构子系统在耦合界面上的边界条件被近似为自由状态,声腔子系统在耦合界面上的边界条件被近似为固定边界。
步骤(2)计算结构子系统和声腔子系统中固有频率低于1.25倍分析频带上限的模态;具体基于有限元法计算了结构子系统和声腔子系统的模态参数。
步骤(3)计算相邻子系统中固有频率低于1.25倍分析频带上限的模态间的耦合参数;具体由下式计算得到:
其中Wmn为结构子系统第m阶位移模态与声腔子系统第n阶声压模态之间的耦合参数,为结构子系统第m阶位移模态的振型,为声腔子系统第n阶声压模态的振型,Ac为结构子系统与声腔子系统之间的耦合界面,s为空间位置。
步骤(4)建立相邻耦合子系统的对偶模态方程:
式中,ω为角频率,i虚数的虚部;Mm为结构子系统第m阶位移模态的模态质量,ωm为结构子系统第m阶位移模态的固有频率;ηm为结构子系统第m阶位移模态的阻尼损耗系数,为结构子系统第m阶位移模态的参与因子,Wmp为结构子系统第m阶位移模态与声腔子系统第p阶声压模态之间的耦合参数,φp(ω)为声腔子系统第p阶声压模态的参与因子,Fm(ω)为结构子系统第m阶位移模态上受到的广义力载荷;Mn为声腔子系统第n阶声压模态的模态质量;ωn为声腔子系统第n阶声压模态的固有频率;ηn为声腔子系统第n阶声压模态的阻尼损耗系数;φn(ω)为声腔子系统第n阶声压模态的参与因子,Wqn为结构子系统第q阶位移模态与声腔子系统第n阶声压模态之间的耦合参数,φq(ω)为结构子系统第q阶位移模态的参与因子,Fn(ω)为声腔子系统第n阶声压模态上受到的广义力载荷;
步骤(5)通过前置处理,获得随机载荷作用下,子系统模态上受到的广义力载荷的互功率谱,具体为:
结构子系统受到随机噪声激励时,结构子系统第k阶位移模态上受到的广义力载荷和结构子系统第l阶位移模态上受到的广义力载荷之间的互谱为:
其中Ap为面压载荷作用面,为结构子系统第k阶位移模态的振型,为结构子系统第l阶位移模态的振型,Spp(s1,s2,ω)为面压载荷的功率谱,s1和s2为空间位置。
步骤(6)计算对偶模态方程,获得所有模态的参与因子的互功率谱;步骤为:
将系统对偶模态方程写成了分块矩阵的形式,并基于下式计算子系统模态参与因子的互功率谱:
其中上标“H”表示共轭转置,
其中Hij为传递函数矩阵(i=1,2;j=1,2),矩阵元素Hij(k,l)的含义为:当第j个子系统中第l阶模态上作用单位广义力时,第i个子系统中第k阶模态参与因子。传递函数矩阵可由下式获得:
其中上标“-1”表示矩阵的逆矩阵,上标“T”表示矩阵的转置,
W(m,n)=Wmn (10)
其中diag()表示对角矩阵,括号内为对角矩阵元素。只有结构受到噪声激励时,子系统模态参与因子的互功率谱为:
其中为结构子系统模态载荷互功率谱矩阵,可基于式(3)计算得到其第m行第n列元素Smn(ω)。
步骤(7)通过模态叠加,计算系统随机声固耦合响应;具体通过下式计算结构子系统的位移响应:
其中
通过下式计算声腔子系统的声压响应:
其中
下面以一个平板/声腔耦合模型为例,具体说明本发明的技术效果,该平板/声腔耦合模型如图2所示。平板的边界条件为:四条边上简支;平板的参数由表1给出:
表1平板的参数取值
声腔的边界条件为:除与平板耦合的面,其余各面为固定边界;声腔的参数由表2给出:
表2声腔的参数取值
在平板面板的外表面施加单位随机噪声载荷,该随机噪声载荷的功率谱Spp(s1,s2,ω)=1。经过上述各步骤,得到加筋板面板上坐标为(0.3m,0.1m)的响应点处的加速度响应功率谱如图3所示,以及声腔中坐标为(0.3m,0.1m,0m)的响应点处的声压响应功率谱如图4所示。
图3和图4中的参考值由有限元直接法计算得到。在对偶模态方程方法分析过程中,选取了2.5kHz以内的平板模态和声腔模态参与响应预示。图3和图4中结果显示,本发明所提供的动响应分析方法,基于对偶模态方程能准确预示随机噪声激励下系统的动响应,有效地解决随机噪声激励下的动响应预示问题,提高分析的效率。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (2)
1.一种随机噪声环境下基于对偶模态方程的动响应分析方法,其特征在于该方法包括步骤:
(1)将声固耦合系统中的结构和声腔划分成在耦合界面上连续耦合的子系统,且耦合界面上相邻的两个子系统分别为声腔子系统和结构子系统;
(2)设置截断频率大于等于分析频率上限的1.25倍,截取结构子系统和声腔子系统中固有频率小于截断频率的模态;
(3)基于有限元法计算截取的各模态的模态参数,模态参数包括:模态质量、阻尼损耗系数和模态振型;
(4)根据各模态参数计算相邻子系统中截取出的模态之间的耦合参数;
(5)根据各子系统的模态参数和相邻子系统间的耦合参数建立相互耦合的两个相邻子系统的对偶模态方程为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&omega;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>i&omega;&eta;</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>&phi;</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>i</mi>
<mi>&omega;</mi>
<munder>
<mi>&Sigma;</mi>
<mi>p</mi>
</munder>
<msub>
<mi>W</mi>
<mrow>
<mi>m</mi>
<mi>p</mi>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>m</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>&infin;</mi>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&omega;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>i&omega;&eta;</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
<mi>&omega;</mi>
<munder>
<mi>&Sigma;</mi>
<mi>p</mi>
</munder>
<msub>
<mi>W</mi>
<mrow>
<mi>q</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>&infin;</mi>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,ω为角频率,i表示虚数的虚部;Mm为结构子系统第m阶位移模态的模态质量;ωm为结构子系统第m阶位移模态的固有频率;ηm为结构子系统第m阶位移模态的阻尼损耗系数;φm(ω)为结构子系统第m阶位移模态的参与因子,Wmp为结构子系统第m阶位移模态与声腔子系统第p阶声压模态之间的耦合参数, 为声腔子系统第p阶声压模态的参与因子,Fm(ω)为结构子系统第m阶位移模态上受到的广义力载荷;
Mn为声腔子系统第n阶声压模态的模态质量;ωn为声腔子系统第n阶声压模态的固有频率;ηn为声腔子系统第n阶声压模态的阻尼损耗系数;为声腔子系统第n阶声压模态的参与因子;Wqn为结构子系统第q阶位移模态与声腔子系统第n阶声压模态之间的耦合参数;φq(ω)为结构子系统第q阶位移模态的参与因子;Fn(ω)为声腔子系统第n阶声压模态上受到的广义力载荷;
(5)将对偶模态方程转换为分块矩阵形式:
<mrow>
<msub>
<mi>S</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>X</mi>
<mn>1</mn>
<mi>H</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>H</mi>
<mrow>
<mn>1</mn>
<mi>F</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>F</mi>
<mi>F</mi>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mrow>
<mn>1</mn>
<mi>F</mi>
</mrow>
<mi>H</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>S</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<msub>
<mi>Y</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>Y</mi>
<mn>2</mn>
<mi>H</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>H</mi>
<mrow>
<mn>2</mn>
<mi>F</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>F</mi>
<mi>F</mi>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mrow>
<mn>2</mn>
<mi>F</mi>
</mrow>
<mi>H</mi>
</msubsup>
</mrow>
其中,
<mrow>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mi>m</mi>
</msub>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mi>n</mi>
</msub>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mn>1</mn>
<mi>F</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>H</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>H</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>H</mi>
<mrow>
<mn>2</mn>
<mi>F</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>H</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>H</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>F</mi>
<mi>F</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mn>1</mn>
<mi>H</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>2</mn>
<mi>H</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,上标“-1”表示矩阵的逆矩阵,上标“T”表示矩阵的转置;Hij为传递函数矩阵,i=1,2,j=1,2;矩阵元素Hij(k,l)表示当第j个子系统中第l阶模态上作用单位广义力时,第i个子系统中第k阶模态参与因子;各传递函数矩阵的计算公式为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>H</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>H</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>H</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>H</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<mrow>
<mi>j</mi>
<mi>&omega;</mi>
<mi>W</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>j&omega;W</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<msub>
<mi>M</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&omega;</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>i&omega;&eta;</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&omega;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>i&omega;&eta;</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
W(m,n)=Wmn
式中,diag()表示对角矩阵,括号内为对角矩阵元素;W(m,n)表示矩阵W中第m行第n列的元素,即结构子系统第m阶位移模态与声腔子系统第n阶声压模态之间的耦合参数Wmn;
(6)计算所述声固耦合系统中只有结构受到噪声激励时,分块矩阵S11和S22满足以下形式:
<mrow>
<msub>
<mi>S</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mn>11</mn>
</msub>
<msub>
<mi>S</mi>
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mn>11</mn>
<mi>H</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>S</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mn>21</mn>
</msub>
<msub>
<mi>S</mi>
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mn>21</mn>
<mi>H</mi>
</msubsup>
</mrow>
式中,SF1F1为结构子系统模态载荷互功率谱矩阵,其第k行第l列元素为Skl(ω),Skl(ω)表示只有结构子系统受到随机噪声激励时,结构子系统第k阶位移模态上受到的广义力载荷和结构子系统第l阶位移模态上受到的广义力载荷之间的互谱,Skl(ω)的计算公式为:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>k</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mo>&Integral;</mo>
<msub>
<mi>A</mi>
<mi>p</mi>
</msub>
</msub>
<msub>
<mo>&Integral;</mo>
<msub>
<mi>A</mi>
<mi>p</mi>
</msub>
</msub>
<msub>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>p</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>ds</mi>
<mn>1</mn>
</msub>
<msub>
<mi>ds</mi>
<mn>2</mn>
</msub>
</mrow>
式中,Ap为面压载荷作用面,为结构子系统第k阶位移模态的振型,为结构子系统第l阶位移模态的振型,Spp(s1,s2,ω)为面压载荷的功率谱,s1和s2为面压载荷作用面Ap上的空间位置;
(7)计算各结构子系统的位移响应和各声腔子系统的声压响应,其中,结构子系统的位移响应的计算公式为:
<mrow>
<msub>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Sw(s,ω)表示第w个结构子系统在位置s处、角频率ω下的位移响应;
声腔子系统的声压响应的计算公式为:
<mrow>
<msub>
<mi>S</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
<msub>
<mi>H</mi>
<mn>21</mn>
</msub>
<msub>
<mi>S</mi>
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mn>21</mn>
<mi>T</mi>
</msubsup>
<msubsup>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
<mi>T</mi>
</msubsup>
</mrow>
<mrow>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,Sp(s,ω)表示第p个声腔子系统在位置s处、角频率ω下的声压响应。
2.根据权利要求1所述的一种随机噪声环境下基于对偶模态方程的动响应分析方法,其特征在于,所述耦合参数的计算公式为:
<mrow>
<msub>
<mi>W</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mo>&Integral;</mo>
<msub>
<mi>A</mi>
<mi>c</mi>
</msub>
</msub>
<msub>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>p</mi>
<mo>~</mo>
</mover>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
式中,Wmn为结构子系统第m阶位移模态与声腔子系统第n阶声压模态之间的耦合参数,为结构子系统第m阶位移模态的振型,为声腔子系统第n阶声压模态的振型,Ac为结构子系统与声腔子系统之间的耦合界面,s为空间位置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710958872.2A CN107748815B (zh) | 2017-10-16 | 2017-10-16 | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 |
US16/647,516 US20200226309A1 (en) | 2017-10-16 | 2018-04-18 | Dynamic response analysis method based on dual-mode equation in random noise environment |
PCT/CN2018/083484 WO2019076015A1 (zh) | 2017-10-16 | 2018-04-18 | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710958872.2A CN107748815B (zh) | 2017-10-16 | 2017-10-16 | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107748815A true CN107748815A (zh) | 2018-03-02 |
CN107748815B CN107748815B (zh) | 2018-08-21 |
Family
ID=61253807
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710958872.2A Active CN107748815B (zh) | 2017-10-16 | 2017-10-16 | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20200226309A1 (zh) |
CN (1) | CN107748815B (zh) |
WO (1) | WO2019076015A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491595A (zh) * | 2018-03-07 | 2018-09-04 | 东南大学 | 一种声-固耦合结构的高频局部响应预示方法 |
CN109145369A (zh) * | 2018-07-11 | 2019-01-04 | 东南大学 | 一种计及非共振传输的中高频局部动响应预示方法 |
WO2019076015A1 (zh) * | 2017-10-16 | 2019-04-25 | 东南大学 | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 |
CN111008417A (zh) * | 2019-11-12 | 2020-04-14 | 中铁大桥勘测设计院集团有限公司 | 一种大长细比连续结构的风致振动分析方法 |
CN112528411A (zh) * | 2020-12-10 | 2021-03-19 | 中国运载火箭技术研究院 | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 |
CN113836773A (zh) * | 2021-09-29 | 2021-12-24 | 九江学院 | 声固耦合系统结构振动响应预示方法、系统及可存储介质 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111709168A (zh) * | 2020-05-29 | 2020-09-25 | 西安交通大学 | 一种基于声固耦合的壳体结构低频声辐射预报方法 |
CN112287575B (zh) * | 2020-10-09 | 2023-11-10 | 航天东方红卫星有限公司 | 一种确定小卫星光学相机环境试验随机振动功率谱的方法 |
CN112883327A (zh) * | 2021-01-12 | 2021-06-01 | 中国海洋大学 | 海上风电结构动力响应高阶精度计算方法 |
CN113343527A (zh) * | 2021-06-09 | 2021-09-03 | 长城汽车股份有限公司 | 车内中频噪声分析方法及介质 |
CN113987854A (zh) * | 2021-09-22 | 2022-01-28 | 东南大学 | 基于降阶模态能量聚类分析的统计能量分析子系统识别方法 |
CN114296353A (zh) * | 2021-12-30 | 2022-04-08 | 长光卫星技术有限公司 | 一种带双轴sada卫星的模态频率计算方法 |
CN115455685B (zh) * | 2022-09-05 | 2023-12-19 | 大连理工大学 | 一种二维紊流风激励下长索结构抖振响应的计算方法 |
CN116186838A (zh) * | 2022-12-27 | 2023-05-30 | 武汉理工大学 | 基于谐小波的结构随机响应分析方法、设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5163015A (en) * | 1989-06-30 | 1992-11-10 | Mitsubishi Denki Kabushiki Kaisha | Apparatus for and method of analyzing coupling characteristics |
US6090147A (en) * | 1997-12-05 | 2000-07-18 | Vibro-Acoustics Sciences, Inc. | Computer program media, method and system for vibration and acoustic analysis of complex structural-acoustic systems |
CN104112070A (zh) * | 2014-07-11 | 2014-10-22 | 长沙理工大学 | 用于弹性边界浅拱发生内共振时动力响应的求解方法 |
CN107133422A (zh) * | 2017-06-01 | 2017-09-05 | 东南大学 | 一种基于对偶模态方程的确定性声固耦合响应预示方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107748815B (zh) * | 2017-10-16 | 2018-08-21 | 东南大学 | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 |
-
2017
- 2017-10-16 CN CN201710958872.2A patent/CN107748815B/zh active Active
-
2018
- 2018-04-18 US US16/647,516 patent/US20200226309A1/en not_active Abandoned
- 2018-04-18 WO PCT/CN2018/083484 patent/WO2019076015A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5163015A (en) * | 1989-06-30 | 1992-11-10 | Mitsubishi Denki Kabushiki Kaisha | Apparatus for and method of analyzing coupling characteristics |
US6090147A (en) * | 1997-12-05 | 2000-07-18 | Vibro-Acoustics Sciences, Inc. | Computer program media, method and system for vibration and acoustic analysis of complex structural-acoustic systems |
CN104112070A (zh) * | 2014-07-11 | 2014-10-22 | 长沙理工大学 | 用于弹性边界浅拱发生内共振时动力响应的求解方法 |
CN107133422A (zh) * | 2017-06-01 | 2017-09-05 | 东南大学 | 一种基于对偶模态方程的确定性声固耦合响应预示方法 |
Non-Patent Citations (3)
Title |
---|
廖涛等: "噪声激励下典型复合材料铆接结构的声固耦合分析", 《固体火箭技术》 * |
张鹏等: "基于模态叠加法的面压载荷与基础激励等效方法", 《东南大学学报(自然科学版)》 * |
李彦斌等: "随机基础激励下承力筒—蒙皮结构的声-固耦合分析", 《宇航学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019076015A1 (zh) * | 2017-10-16 | 2019-04-25 | 东南大学 | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 |
CN108491595A (zh) * | 2018-03-07 | 2018-09-04 | 东南大学 | 一种声-固耦合结构的高频局部响应预示方法 |
CN108491595B (zh) * | 2018-03-07 | 2019-03-29 | 东南大学 | 一种声-固耦合结构的高频局部响应预示方法 |
CN109145369A (zh) * | 2018-07-11 | 2019-01-04 | 东南大学 | 一种计及非共振传输的中高频局部动响应预示方法 |
CN111008417A (zh) * | 2019-11-12 | 2020-04-14 | 中铁大桥勘测设计院集团有限公司 | 一种大长细比连续结构的风致振动分析方法 |
CN111008417B (zh) * | 2019-11-12 | 2022-11-25 | 中铁大桥勘测设计院集团有限公司 | 一种大长细比连续结构的风致振动分析方法 |
CN112528411A (zh) * | 2020-12-10 | 2021-03-19 | 中国运载火箭技术研究院 | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 |
CN112528411B (zh) * | 2020-12-10 | 2022-09-06 | 中国运载火箭技术研究院 | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 |
CN113836773A (zh) * | 2021-09-29 | 2021-12-24 | 九江学院 | 声固耦合系统结构振动响应预示方法、系统及可存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107748815B (zh) | 2018-08-21 |
WO2019076015A1 (zh) | 2019-04-25 |
US20200226309A1 (en) | 2020-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107748815A (zh) | 一种随机噪声环境下基于对偶模态方程的动响应分析方法 | |
CN107133422B (zh) | 一种基于对偶模态方程的确定性声固耦合响应预示方法 | |
CN105843073B (zh) | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 | |
CN107357992B (zh) | 基于聚类分析的复合材料结构有限元模型修正方法 | |
CN107368630B (zh) | 一种获取耦合损耗因子的数值方法 | |
Wang et al. | A coupled smoothed finite element method (S-FEM) for structural-acoustic analysis of shells | |
Yan et al. | Analysis and experimental validation of the middle-frequency vibro-acoustic coupling property for aircraft structural model based on the wave coupling hybrid FE-SEA method | |
CN110162826B (zh) | 薄壁结构热气动弹性动响应分析方法 | |
CN107368634A (zh) | 一种复杂结构耦合损耗因子的预示方法 | |
CN105183958B (zh) | 一种复合材料层合结构三维振动分析方法 | |
CN103942381B (zh) | 用于飞机铝合金结构性能预测的状态近场动力学方法 | |
CN107782547A (zh) | 一种传动带的模态分析方法及装置 | |
CN113466335A (zh) | 一种超声导波蜂窝夹芯结构脱粘层诊断方法 | |
CN107256295B (zh) | 一种计及非共振传输的中频动响应预示方法 | |
CN108984879A (zh) | 一种串联多自由度系统的位移频率响应计算方法 | |
CN109145369B (zh) | 一种计及非共振传输的中高频局部动响应预示方法 | |
CN105184060A (zh) | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 | |
CN106021630B (zh) | 一种适用于新型空间飞行器的结构/阻尼一体化设计方法 | |
CN107220444B (zh) | 一种计及非共振传输的中频动响应预示简化分析方法 | |
CN106295159A (zh) | 一种基于自相关函数的结构风振响应高效频域估计方法 | |
CN109946029B (zh) | 一种涡轮增压器的多频疲劳试验方法 | |
CN112434385B (zh) | 一种蜂窝夹层板有限元建模方法 | |
CN108491595B (zh) | 一种声-固耦合结构的高频局部响应预示方法 | |
CN112149323B (zh) | 一种卫星噪声预示的方法及装置 | |
CN114218672A (zh) | 跨超声速大迎角配平翼的抖振响应获取方法及相关装置 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |