复杂边界条件下功能梯度圆环板平稳随机振动响应特性分析

左朋 石先杰 葛任伟 罗景润

摘要： 提出一种可用于分析复杂边界条件下功能梯度圆环板平稳随机振动响应特性的谱几何法？虚拟激励法（spec？tro geometric method？pseudo excitation method，SGM？PEM）。采用沿圆环板边界均匀分布的边界约束弹簧来模拟复杂边界条件，通过虚拟激励法将平稳随机载荷转化为虚拟简谐载荷。在一阶剪切变形理论框架下，采用以简洁三角函数为内核的谱几何法来描述圆环板结构的位移容许函数。基于 Rayleigh？Ritz 法推导了平稳随机激励作用下功能梯度圆环板的动力学分析模型。通过与有限元法结果对比分析，验证了文中构建的分析模型的有效性和准确性。分析了梯度指数、厚度参数、边界条件等因素对功能梯度圆环板平稳随机振动响应特性的影响规律。

关键词： 随机振动；功能梯度圆环板；平稳随机激励；谱几何法？虚拟激励法；复杂边界条件

中图分类号： O324；TB532 文献标志码： A 文章编号： 1004-4523（2023）03-0634-11

DOI：10.16385/j.cnki.issn.1004-4523.2023.03.005

引言

作為一种新型复合材料，功能梯度材料自被提出便受到国内外学者的广泛关注。而作为基本结构单元，由功能梯度材料所制造的圆环板在航空航天工程中被广泛应用，如卫星的隔舱壁，发动机的端盖等。近年来，功能梯度圆环板结构动力学特性分析成为研究热点，国内外研究学者针对该问题开展了大量研究，提出了一系列分析方法。？ur［1］采用拟格林函数法分析了弹性支承功能梯度圆环板的自由振动特性。Zhao 等［2］采用 Rayleigh？Ritz 法获得了任意边界条件下功能梯度厚圆环板的自由振动三维解。Jodaei 等［3］采用基于状态空间的微分求积法和人工神经网络技术，分别对不同边界条件下功能梯度圆环板的自由振动特性进行分析，并对比了两种方法获取的固有频率结果，验证了人工神经网络技术在预示功能梯度圆环板自由振动特性方面的有效性。基于三维弹性理论，Jin 等［4］采用改进傅里叶级数法研究了任意边界条件下功能梯度圆环板的自由振动特性。Malekzadeh 等［5］基于微分求积法和三维弹性理论求解了热环境下功能梯度厚圆环板自由振动解。滕兆春等［6？7］基于微分求积法分析了功能梯度圆环板面内自由振动特性。针对材料属性沿径向变化的功能梯度圆环板，吕朋等［8］采用 Rayleigh？Ritz法分析了其在热环境下的面内自由振动特性，并在研究中考虑了弹性边界条件。Dong［9］基于 Ritz 法求解了功能梯度圆环板的三维自由振动解。

通过上述文献分析可以发现，虽然国内外学者针对功能梯度圆环板的动力学特性求解开展了大量工作，但是研究成果主要集中在确定性振动分析方面。然而，航空航天工程结构在实际服役期间，通常会经受各种形式的随机载荷作用，从而使得结构产生复杂的动力学行为，甚至出现疲劳破坏现象。此外，功能梯度圆环板作为工程结构的基本单元，其在载荷作用下的振动特性影响着工程结构的动力学行为。因此，有必要开展随机载荷作用下功能梯度圆环板相应的振动响应特性研究，为实际工程结构设计提供相应的理论支撑。

近年来，学者针对随机激励下板壳结构的振动响应特性开展了相关研究。在研究中采用了有限元法［10？12］，谱有限元法［13］等方法对该问题进行求解，但这些方法普遍存在计算耗时长的问题。随后，一些学者将虚拟激励法［14］用于研究板壳结构随机振动响应特性，大大缩短了计算时间。孙悦等［15］结合有限元法和虚拟激励法对硬涂层薄板随机振动响应特性进行分析，并讨论了硬涂层的减振效果。戴新进［16］结合有限元法和虚拟激励法对平稳/非平稳随机激励作用下复合材料层合板的振动响应行为开展了研究。Chen 等采用离散解析法和虚拟激励法研究获得了各向同性矩形板［17］以及各向同性圆柱壳［18］的随机振动响应解析解。Zhou 等［19］采用改善的 Ritz 法和虚拟激励法研究了功能梯度压电矩形板的非平稳随机振动响应特性。

虽然目前针对板壳结构的随机振动响应分析取得了一定的研究成果，但大部分研究工作集中在矩形板。而且目前的分析手段多采用有限元方法与虚拟激励法相结合的方式，虽然可以降低计算耗时，但有限元法对于铆接、螺栓等结构复杂边界条件的等效处理、结构快速参数化分析以及结构平稳随机振动响应分析的能力与实际相比还存在一定差距。因此，直面实际工程需要，发展一套具备有效处理复杂边界条件及快速参数化分析能力，并同时拥有计算耗时短、收敛特性优、结果保真度高的建模分析方法，具有一定理论和工程应用价值。

谱几何法［20］具有形式简洁、收敛快速、求解精度高以及满足任意边界条件等优点，近年来被广泛用于板壳结构的振动特性研究中。Jiang 等［21］在改进的傅里叶级数方法的基础上将正弦函数作为传统傅里叶级数的辅助项，提出了谱几何法。随后文献［20，22？25］基于谱几何法分析了环板、圆环板等回转板结构，并建立了旋转结构振动特性的统一分析模型。在此基础上，Zhao 等［26］分析了圆锥？圆柱？圆球组合壳结构的自由振动特性。此后，谱几何法被应用到复合材料回转组合壳结构［27？28］的自由振动和受迫振动分析中。

综合谱几何法满足任意边界条件以及虚拟激励法处理随机载荷的优势，提出一种高精度且易于参数化研究的高效计算分析方法，用于研究任意边界条件下功能梯度圆环板结构平稳随机振动问题。该方法采用谱几何法将圆环板位移容许函数描述为形式简洁的三角函数，并通过引入正弦函数项来消除边界处位移容许函数在边界处微分不连续的问题。采用虚拟激励法将平稳随机载荷转化为确定性载荷，并结合一阶剪切变形理论建立平稳随机激励下功能梯度圆环板的动力学分析模型。通过与有限元数值模拟结果对比，验证了文中方法的有效性和准确性。最后，探讨了梯度指数、厚度、边界条件等因素对功能梯度圆环板随机振动响应特性的影响。

1 理论推导

1. 1 结构模型描述

图 1 为功能梯度圆环板模型。圆环板的参考面与其中面一致，并且正交坐标系（r，θ，z）位于功能梯度圆环板结构的参考面上。其中，r，θ 和 z 分别表示圆环板的径向、周向和厚度方向。此外，R0和 R1分别表示圆环板结构的内径和外径，径向方向的长度为：R=R1-R0，而 h 则为圆环板厚度。通过在圆环板内径和外径的边界处均匀布置具有可调刚度系数的不同类型边界约束弹簧来实现对任意边界条件的模拟。kui，kvi和 kwi（i=0，1）表示三组平移约束弹簧，Kri和 Kθ（i i=0，1）表示两组旋转约束弹簧，其中，下标“0”和“1”分别表示功能梯度圆环板结构的内径和外径边界。在后续研究中，假设功能梯度圆环板经受沿 z 方向的基础加速度平稳随机激励作用或集中点平稳随机激励作用。

文中所研究的功能梯度圆环板结构上下表面材料分别为陶瓷和金属，并且功能梯度材料属性 P 在圆环板厚度 z 方向上呈现梯度变化，其变化关系如下式所示：式中 Pc和 Pm分别表示陶瓷和金属的材料参数（包括弹性模量 E，泊松比 ν 和密度 ρ）；p 表示材料的梯度指數。

1. 2 能量方程

根据一阶剪切变形理论，垂直于结构中面的法线在变形后依然为直线且长度不变，但是不一定垂直于变形后的中面［29］，因此功能梯度圆环板结构的位移场可以描述为：

1. 3 位移容许函数

为了满足位移容许函数在整个求解域内（包括边界处）充分光滑的条件，使得构造的位移容许函数能够克服不同边界条件下其微分在边界处可能存在的不连续问题，文中基于谱几何法采用余弦函数和正弦函数辅助项来描述功能梯度圆环板径向坐标下的位移容许函数，同时考虑到圆环板结构的空间对称性，采用正余弦函数描述周向坐标下的位移函数：q （ r，θ，t ）=

在一阶剪切变形理论下，圆环板的振动位移由其二阶控制微分方程决定，为保证位移函数在整个求解域内的一阶导数连续并且二阶导数存在，在径向方向引入了两项傅里叶正弦函数辅助项。从数学上可知，傅里叶余弦函数可以用来表示任何在求解域中完全可积的连续函数，并且余弦函数可以实现逐项微分。但是当边界条件中存在弯矩等条件时，若直接采用傅里叶余弦函数表示，其一阶导数在结构边界处会出现数值间断的现象。因此，为了克服数值不连续的问题，引入正弦函数辅助项消除数值间断的现象，使得位移函数本身或其导数在整个求解域内的导数连续或存在。

1. 4 求 解

基于上述构建的能量方程，功能梯度圆环板结构的拉格朗日函数可表示为：

Ξ = T - U - V + Wξ （15）式中 ξ= f，g。将式（7）～（14）代入式（15），并采用Rayleigh？Ritz法对位移容许函数的未知级数展开系数求偏导，获得平稳随机激励下功能梯度圆环板的动力学控制方程：

MQ？+ KQ = Fξ （16）式中 M 为圆环板的整体质量矩阵；K 表示圆环板的整体刚度矩阵；Q 表示未知级数展开系数向量；Fξ为 虚 拟 激 励 力 的 载 荷 向 量 ，可 以 表 示 为 Fξ =[ ] 0，0，F Wξ，0，0 T，其中 F Wξ 为沿 z 方向的集中点载荷或基础加速度载荷向量，其表达式为：将式（16）中求得的未知级数展开系数向量 Q代入式（14）中，便可得到虚拟位移响应 u，v，w，φr 和φθ，从而可以方便地获得平稳随机激励作用下功能梯度圆环板的振动响应功率谱密度。以圆环板横向位移响应为例，位移响应自功率谱函数 Sww （ω），速度响应自功率谱函数 Sw？ w？ （ω） 以及加速度响应自功率谱函数 Sw？w？（ω） 可以分别描述为：

相应地，平稳随机激励下功能梯度圆环板的响应均方根值（RMS）可按下式计算获得［18］：式中 ωl和 ωu分别表示所施加随机激励频率范围的下限和上限；Δω 代表频率步长；G ww （ω），G w？ w？ （ω） 和G w？w？（ω） 为单边功率谱。

式（20）对于随机振动响应均方根值的求解是通过对频率范围进行离散而计算得到的，NM 即为频率范围的离散频率点数，其和频率步长有如下关系：

为了避免数值计算的奇异性，文中采用复杨氏模量的方式引入结构阻尼：E'= E （1 + iη） （22）式中 E'为复杨氏模量；η 为结构损耗因子，在下述数值分析中，结构损耗因子 η=0.1。

2 振动求解与分析

在理论模型推导基础上，本节将对平稳随机激励下功能梯度圆环板的随机振动响应特性进行分析和讨论。在接下来的研究中，假设功能梯度圆环板的默认几何尺寸为：R0=0.5 m，R1=2.0 m，h=0.1 m。圆 环 板 上 下 表 面 的 材 料 分 别 为 Si3N4 和 SUS304 ，材 料 参 数 为 ：Ec=348.43 GPa ，ν c=0.24 ，ρ c=2370 kg/m3，Em=201.04 GPa，νm=0.3262，ρm=8166 kg/m3。文中研究考虑了固支（C）、简支（SS）、剪切（SD）三种经典边界条件和 E1，E2，E3，E4，E5，E6 六种弹性边界条件。不同边界条件下的相应弹簧刚度值参考文献［31］选取，相应参数如表 1 所示。此外，根据文献［20］中相关收敛性分析，文中位移容许 函 数 在 径 向 和 周 向 上 的 截 断 数 被 设 定 为 ：M=N=14。

2. 1 有效性验证

表 2 给出了不同边界条件下功能梯度圆环板的前 10 阶固有频率（Hz）。材料的梯度指数 p=1.0，圆环板的其余参数与默认参数一致，还考虑了 C？C，SS？SS，C？E1，C？E2，E1？E4 和 E2？E5 等 6 种 不 同 的边界条件。其中，以 C？C 边界条件为例，C？C 中第一个“C”代表圆环板内径边界的约束边界条件类型为固支，第二个“C”代表圆环板外径边界的约束边界条件类型为固支。通过与有限元结果对比发现，两种计算方法获得的结果吻合良好。

接下来，对采用文中方法分析功能梯度圆环板平稳随机振动的有效性进行验证。其中，集中点激励 为 功 率 谱 密 度（Power Spectral Density，PSD）S=1 N2/Hz 的 限 带 白 噪 声 谱 ，其 频 率 范 围 为［20，2000］ Hz；基础加速度激励考虑两种类型的激励载荷，Ⅰ类是功率谱密度（PSD）S0=1 g2/Hz 的限带白噪声谱，其频率范围为［20，2000］ Hz；Ⅱ类则为航空航天工程中的特定功率谱［18］，其形式如图 2 所示，频率与基础加速度激励 PSD 之间的数学关系可以表示为：式中 Sref表示频率范围 （150，700 ） Hz 内 PSD 值为0.04 g2/Hz；Ω1和 Ω2分别表示频率范围 （ 20，150 ） Hz和（700，2000 ） Hz内 PSD 曲线斜率分别为+6 dB/oc？tave和？3 dB/octave。

同时需要指出的是，文中采用向量［r，θ］来表示响应测点的位置，其中，r 和 θ 分别表示功能梯度圆环板径向和周向的位置。在后续计算分析中，测点响应的径向位置选为圆环板内径和外径的中点，而周向位置选取 θ=0。

图3 分别给出了集中点激励、Ⅰ类基础加速度激励和Ⅱ类基础加速度激励作用下 C？C 功能梯度圆环板随机振动响应（功率谱密度）对比结果。功能梯度圆环板的几何尺寸与默认尺寸一致，材料的梯度指数 p=1.0。通过图 3 可以看出，文中基于谱几何法？虚拟激励法研究获得的随机振动响应结果与有限元结果吻合较好，说明文中方法可以有效预示平稳随机激励作用下功能梯度圆环板随机振动响应特性。此外，表 3 还给出了基于谱几何法？虚拟激励法获得的集中点激励和两类基础加速度激励下 C？C 功能梯度圆环板的随机振动响应均方根值。相应地，有限元法结果也作为参考列在表中。从表 3 中可以看出，文中方法获得的随机振动响应均方根值与有限元结果具有良好的一致性，进一步验证了文中方法在预示功能梯度圆环板随机振动响应特性方面的有效性。同时，从表 3 的计算时间对比情况可知，文中方法的计算耗时约为有限元法的 0.2%，说明文中方法具有较好的计算效率。

图 4 则给出了Ⅱ类基础加速度激励作用下功能梯度圆环板随机振动响应（加速度功率谱）的对比情况。圆环板的所有参数与表2保持一致 ，考虑了C？E1，E1？E4，E4？E6，C？E2，E2？E5，E5？E6 六种不同边界条件。从图 4 的对比情况可以看出，文中模型对不同边界条件下功能梯度圆环板随机振动加速度功率谱密度响应的计算结果与有限元法结果均吻合良好。因此，文中构建的分析模型可以有效分析不同边界条件下功能梯度圆环板的平稳随机振动响应特性。

2.2 参数化分析

在验证文中方法有效性的基础上，研究关键参数对平稳随机激励下功能梯度圆环板振动响应特性的影响规律。为了简便起见，后续算例中的外界激励载荷选定为Ⅱ类基础加速度激励载荷。

首先，图 5 给出了不同梯度指数 p 下 C？C 功能梯度圆环板的随机振动响应特性。相关几何参数、材料参数和响应点位置均与默认参数保持一致。从图5 中可以发现，随着梯度指数 p 增大，圆环板随机振动响应（位移，速度，加速度）曲线的共振峰值将逐渐向低频段移动，这是因为 p 增大将导致圆环板结构的材料属性由陶瓷成分向金属成分变化，使得结构弯曲刚度不断降低，结构质量逐渐增加，从而降低了结构的固有频率。此外，以 p=9 时圆环板随机振动位移响应和速度响应在基频处的峰值作为参考，p=0，1，3，5，7 时位移响应和速度响应在基频处的峰值分 别 为 p=9 時 的 4.49%，34.47%，66.57%，82.69%，92.50% 和 21.18%，57.62%，80.33%，90.16%，96.27%。说明随着 p 增大，随机振动位移响应和速度响应基频处的峰值逐渐增加。但是，梯度指数 p 对随机振动加速度响应在基频处的峰值影响很小，在不同梯度指数 p 下，加速度响应基频处的峰值最小值为最大值的 96.42%。

圆环板外径与内径之比对随机振动响应特性的影响如图 6 所示，其中圆环板边界条件为 C？C。在该算例中，圆环板外径 R1=2 m 保持不变，外径与内径之比 R1/R0=2，3，4，6，8。测点响应位置选为圆环板内径和外径的中点，其余参数与默认参数一致。从图 6 中可以看出，R1/R0 增加将会使得圆环板随机振动响应（位移，速度和加速度）曲线的共振峰值向低频段移动。这是由于 R1/R0 增大降低了圆环板结构的弯曲刚度。此外，以 R1/R0=8 时圆环板随机振动位移响应和速度响应在基频处 的 峰 值 作 为 参 考 ，R1/R0=2，3，4，6 时 位 移 响 应和速度响应在基频处的峰值分别为 R1/R0=8 时的1.50%，12.45%，30.45%，68.41% 和 13.21%，36.61%，56.22%，83.44%。 说 明 随 着 R1/R0 增 大 ，随机振动位移响应和速度响应基频处的峰值逐渐增 加 。 对 加 速 度 响 应 而 言 ，R1/R0 对 其 影 响 较 小 ，在不同 R1/R0下，加速度响应基频处的峰值最小值为最大值的 86.87%。

图 7 则给出了不同边界条件对平稳随机激励下功能梯度圆环板的振动响应特性。该算例考虑了C？C，SS？SS，SD？SD 三 种 经 典 边 界 条 件 和 C？E1，C？E2，C？E3，E1？E4，E2？E5，E3？E6 六种组合边界条件。圆环板相关几何参数和材料参数与默认参数保持一致。对于经典边界条件而言，当边界条件从C？C 变为 SS？SS 或者 SD？SD 时，圆环板结构的弯曲刚度逐渐减小，降低了结构的固有频率，从而使得随机振动响应曲线共振峰值向低频段移动。此外，通过对比 C？C 和 SS？SS 下的随机振动响应结果可以发现，旋转约束弹簧 Kr的刚度值对圆环板的随机振动响 应 特 性 有 较 大 影 响 。 同 时 ，SS？SS 边 界 条 件 和SD？SD 边界条件对圆环板随机振动响应特性影响基本一致。这说明径向方向的约束弹簧刚度对结构随机振动响应特性影响较小 。 C？E3，E2？E5 和E3？E6 边界条件会使得随机振动响应曲线共振峰值向低频段移动，而且会明显改变圆环板随机振动响应能量在各频率上的分布密度 ，而 C？E1，C？E2 和 E1？E4边界条件对圆环板随机振动响应特性影响基本一致。从而说明厚度方向（z 方向）的弹簧刚度值对圆环板随机振动响应特性影响较大，而周向方向的弹簧刚度值对随机振动响应特性影响较小。

最后，分析了厚度参数对平稳随机激励作用下C？C 功能梯度圆环板振动响应特性的影响，相应结果如图 8 所示。圆环板的厚度分别设定为 0.050，0.075，0.100，0.125 和 0.150 m。其余参数和响应位置均与默认参数值保持一致。从图 8 中可以看出，随着厚度 h 增加，圆环板随机振动响应（位移，速度和加速度）曲线的共振峰值将向高频段移动。这是由于随着厚度 h 增加，圆环板结构的弯曲刚度和质量均提高，但是，其对刚度矩阵的贡献要大于对质量矩阵的贡献，从而提升了结构的固有频率。此外，厚度 h 的增加会明显降低位移响应和速度响应的共振峰值。在不同的厚度 h 下，加速度响应在基频处峰值的最小值为最大值的 80.40%，表明相较于梯度指数 p 和外径与内径之比 R1/R0，厚度 h 对加速度响应的影响更大。

3 结 论

文中基于谱几何法？虚拟激励法和一阶剪切变形理论，建立了平稳随机激励下功能梯度圆环板振动分析模型。通过数值算例分析对比，验证了文中所构建模型的有效性和准确性。在此基础上，探讨了材料参数、结构几何参数、边界条件等因素对功能梯度圆环板随机振动响应特性的影响，获得结论如下：

（1）梯度指数 p 和圆环板外径与内径之比 R1/R0增加，将会降低圆环板结构的弯曲刚度，使得结构固有频率降低，从而使得随机振动响应曲线的共振峰值向低频段移动。

（2）比较不同参数下加速度响应基频处峰值最大值和最小值之间的比例，发现相较于梯度指数 p和外径与内径之比 R1/R0，厚度 h 对加速度响应能量在各频率上的分布影响最大。

（3）径向约束弹簧 ku和轴向约束弹簧 kv刚度值的改变基本不会影响圆环板随机振动响应能量在各频率上的分布密度；而改变厚度方向约束弹簧 kw和旋转约束弹簧 Kr的刚度值会对结构边界势能产生很大影响，从而改变圆环板随机振动响应能量在各频率上的分布密度。

（4）在平稳随机激励作用下，功能梯度圆环板结构的随机振动位移响应、速度响应和加速度响应能量主要集中在结构基频处。此外，随机激励的高频对位移响应和速度响应的影响较小，但对加速度响应有较大影响。

参考文献：

［1］ ？ur K K. Free vibration analysis of elastically supported functionally graded annular plate via quasi-Greens func？tion method［J］. Composite Part B： Engineering，2018，144：37-55.

［2］ Zhao J， Zhang Y K， Choe K， et al. Three-dimensional exact solution for the free vibration of thick functionally graded annular sector plates with arbitrary boundary con？ditions［J］. Composite Part B： Engineering， 2019，159：418-436.

［3］ Jodaei A， Jalal M， Yas M H. Free vibration analysis of functionally graded annular plates by state-space based differential quadrature method and comparative model？ing by ANN［J］. Composite Part B： Engineering，2012，43（2）：340-353.

［4］ Jin G Y， Su Z， Ye T G， et al. Three-dimensional free vibration analysis of functionally graded annular sector plates with general boundary conditions［J］. Composite Part B： Engineering，2015，83：352-366.

［5］ Malekzadeh P， Shahpari S A， Ziaee H R. Three-dimen？sional free vibration of thick functionally graded annular plates in thermal environment［J］. Journal of Sound and Vibration，2010，329（4）：425-442.

［6］ 滕兆春，蒲育，房曉林 . FGM 圆环板面内自由振动的DQM 求 解［J］. 北 京 理 工 大 学 学 报 ，2014，34（12）：1211-1216.

TENG Zhaochun， PU Yu， FANG Xiaolin. In-plane free vibration analysis for FGM annular plates by differ？ential quadrature method［J］. Transactions of Beijing In？stitute of Technology，2014，34（12）：1211-1216.

［7］ 滕兆春，蒲育 . 温度影响下 FGM 圆环板的面内自由振动分析［J］. 振动与冲击，2015，34（9）：210-217.TENG Zhaochun， PU Yu. In-plane free vibration of FGM annular plates considering temperature effect［J］.Journal of Vibration and Shock，2015，34（9）：210-217.

［8］ 呂朋，杜敬涛，邢雪，等 . 热环境下弹性边界约束 FGM圆环板面内振动特性分析［J］. 振动工程学报，2017，30（5）：713-723.

Lü Peng， DU Jingtao， Xing Xue， et al. Study on in plane vibration characteristics of elastically restrained FGM annular panel in thermal environment［J］. Journal of Vibration Engineering，2017，30（5）：713-723. ［9］ Dong C Y. Three-dimensional free vibration analysis of functionally graded annular plates using the Chebyshev-Ritz method［J］. Materials and Design，2008，29（8）：1518-1525.

［10］ Esmailzadeh M， Lakis A A， Thomas M， et al. Predic？tion of the response of a thin structure subjected to a tur？bulent boundary-layer-induced random pressure field［J］. Journal of Sound and Vibration，2009，328（1-2）：109-128.

［11］ Chang T P， Chang H C， Liu M F. A finite element analysis on random vibration of nonlinear shell structures［J］. Journal of Sound and Vibration，2006，291（1-2）：240-257.

［12］ Franco F， De Rosa S， Ciappi E. Numerical approxima？tions on the predictive responses of plates under stochas？tic and convective loads［J］. Journal of Fluids and Struc？tures，2013，42：296-312.

［13］ Birgersson F， Ferguson N S， Finnveden S. Application of the spectral finite element method to turbulent bound？ary layer induced vibration of plates［J］. Journal of Sound and Vibration，2003，259（4）：873-891.

［14］ Lin J H， Zhao Y， Zhang Y H. Accurate and highly effi？cient algorithms for structural stationary/non-stationary random response［J］. Computer Methods in Applied Mechanics and Engineering，2001，191（1-2）：103-111.

［15］ 孙悦，刘晓峰，孙伟 . 随机激励作用下硬涂层薄板振动有限元分析及减振预估［J］. 振动与冲击，2022，41（4）：63-69.

SUN Yue， LIU Xiaofeng， SUN Wei. Vibration analy？sis of a hard-coating thin plate under random excitation using FEM and estimation of vibration reduction effect［J］. Journal of Vibration and Shock，2022，41（4）：63-69.

［16］ 戴新进 . 复合材料结构随机振动的虚拟激励法及在航空航天领域的应用［D］. 大连：大连理工大学，2007.

DAI Xinjin. PEM based random vibration analysis of composite structures and its application in aero/astro nautical engineering［D］. Dalian： Dalian University of Technology，2007.

［17］ Chen G H， Zhou J L， Yang D X. Benchmark solutions of stationary random vibration for rectangular thin plate based on discrete analytical method［J］. Probabilistic En？gineering Mechanics，2017，50：17-24.

［18］ Chen G H， Huo H， Zhan S X， et al. Analytical stochas？tic responses of thin cylindrical shells under various sta？tionary excitations［J］. International Journal of Mechani？cal Sciences，2021，190：106048.

［19］ Zhou K， Hu Z M， Hua H X. Investigation of the non？stationary stochastic response of functionally graded piezoelectric material plates with general boundary con？ditions［J］. Applied Mathematical Modelling， 2021，96：315-335.

［20］ 石先杰 . 復杂边界条件下旋转结构统一动力学模型的构建与研究［D］. 哈尔滨：哈尔滨工程大学，2014.

SHI Xianjie. The construction and analysis on unified dynamical model of revolve structures subjected to com？plex boundary conditions［D］. Harbin： Harbin Engineer？ing University，2014.

［21］ Jiang S L， Li W L， Yang T J. A spectro-geometric method for the vibration analysis of built-up structures［A］. Proceedings of the INTER-NOISE and NOISE CON Congress and Conference Proceedings［C］. Den？ver， Colorado： Institute of Noise Control Engineering，2013：948-953.

［22］ Shi X J， Shi D Y， Li W L， et al. A unified method for free vibration analysis of circular， annular and sector plates with arbitrary boundary conditions［J］. Journal of Vibration and Control，2016，22（2）：442-456.

［23］ 史冬岩，石先杰，李文龙 . 任意边界条件下环扇形板面内振动特性分析［J］. 振动工程学报，2014，27（1）：1-8.

Shi Dongyan， Shi Xianjie， Li Wenlong. In-plane vibra？tion analysis of annular sector plates with arbitrary boundary supports［J］. Journal of Vibration Engineer？ing，2014，27（1）：1-8.

［24］ 石先杰，李春丽，史冬岩 . 环板结构面内自由振动特性分析［J］. 振动工程学报，2016，29（3）：465-471.

Shi Xianjie， Li Chunli， Shi Dongyan. Free in-plane vi？bration analysis of annular plates［J］. Journal of Vibra？tion Engineering，2016，29（3）：465-471.

［25］ Shi X J， Li C L， Wang F J， et al. A unified formulation for free transverse vibration analysis of orthotropic plates of revolution with general boundary conditions［J］. Mechanics of Advanced Materials and Structures，2018，25（2）：87-99.

［26］ Zhao Y K， Shi D Y， Meng H. A unified spectro-geo？metric-Ritz solution for free vibration analysis of conical cylindrical-spherical shell combination with arbitrary boundary conditions［J］. Archive Applied Mechanics，2017，87：961-988.

［27］ Shi X J， Zuo P， Zhong R， et al. Vibration analysis of combined functionally graded cylindrical-conical shells coupled with annular plates in thermal environment［J］.Composite Structures，2022，294：115738.

［28］ Zuo P， Shi X J， Ge R W， et al. Unified series solution for thermal vibration analysis of composite laminated joined conical-cylindrical shell with general boundary conditions［J］. Thin-Walled Structures， 2022， 178：109525.

［29］ 张红 . 回转类复合材料板结构与封闭声腔耦合系统建模方法及其特性研究［D］. 哈尔滨：哈尔滨工程大学，2018.

ZHANG Hong. The modeling method and characteris？tic analysis for rotary composite laminated plate andacoustic cavity coupled system［D］. Harbin： Harbin En？gineering University，2018.

［30］ Qatu M S. Vibration of Laminated Shells and Plates［M］. Amsterdam， The Netherlands： Academic Press，2004.

［31］ Su Z， Jin G Y， Shi S X， et al. A unified solution for vi？bration analysis of functionally graded cylindrical， coni？cal shells and annular plates with general boundary con？ditions［J］. International Journal of Mechanical Scienc？es，2014，80：62-80.