0 引　言

复合材料加筋/夹层板壳结构由于具有比刚度高、比强度大、耐腐蚀以及具有高度可设计性等诸多优点，已在航空航天、船舶、高速列车、建筑等工程领域中得到广泛应用。目前针对此类结构的研究广泛而深入，从理论建模分析，到特定结构优化设计，再到制造工艺探讨，各环节都有大量的文献涉及。根据应用背景的不同，复合材料加筋/夹层结构的研究又通常涉及力学、声学、材料科学、控制科学等不同学科，通常较为复杂。

复合材料加筋/夹层结构声振特性研究是建立在普通板壳结构研究基础上的，从普通板壳到复合材料加筋/夹层板壳，研究所考虑的因素逐渐增多，研究方法也在继承和发展。基于由简入繁的思路，本文首先介绍普通板壳结构的声振特性研究，再逐步过渡到复合材料加筋/夹层板壳的声振特性研究。其后介绍新材料的应用给结构声振特性研究带来了新特点。最后，由于工程结构的应用离不开一定的边界条件，探讨复合材料加筋/夹层板壳在任意边界条件下的声振特性研究。通过归纳和梳理相关文献，对复合材料加筋/夹层板壳声振特性研究以及研究所采用的各种方法进行总结，并指出相关研究的发展趋势，为相关研究者提供参考。

1 结构声振特性研究现状 1.1 研究内容

结构的声振特性是指考虑结构-流体耦合作用的声振耦合特性，包含振动特性和声学特性2个部分，其研究的内容包括耦合系统中的结构振动，结构声辐射，结构声传递。在声振耦合系统中，结构振动通常由外部激励引起，用均方振速、振动加速度级等指标表征；结构声辐射由结构振动引起的声场变化造成，特征指标有辐射声压、辐射声功率和声辐射效率等；结构声传递特性是指声波通过结构后发生反射、折射、衰减从而呈现出的不同特性，通常用传递损失（TL）表征。归纳目前研究结构声振特性的文献可以发现，当结构所受外部激励为各种入射声波时，关注点通常是结构的声传递特性^[1-7]，而当外部激励为外力激励时，关注点是结构的振声响应^[8-11]。不论结构受外力激励还是声波激励，都是外部激励作功使结构振动，结构振动又引起结构周围声场的声压变化，从而将动能转化成声能辐射出去，在这个过程中结构振动和声压变化相互影响，从而形成声振耦合系统^[12]。

1.2 研究方法

研究结构声振特性的方法有解析法、半解析法、数值法。解析法适用于简单结构如板壳结构、圆柱结构等，主要包括空间谐波法^[2,13-15]、傅里叶变化法^[16-19]、模态分析法^[11,20-21]、行波法^[22]、传递矩阵法^[23]等，其中空间谐波法和傅里叶变化法适用于无限大周期结构，二者在本质上都是无限大代数方程组截断求解问题^[24]；模态分析法能利用模态正交性简化理论分析，但要实现模态解耦对结构的边界条件有所要求，主要用于有限大简支结构在中低频的声振分析；行波法采用波动形式的解来求解结构的声振耦合控制方程，可以对结构弯曲波的传播过程进行细致的机理分析。

半解析法是解析与数值相结合的求解方法。在求解结构声振特性时，计算涉及结构和声场2个域，当面对复杂结构或复杂边界时，难以获得解析解，利用离散化求解可以简化分析，将解析解和离散解组合起来的半解析法，可以兼具灵活性和计算效率，因而在研究中广泛应用^[25-28]。

对于数值方法，有限元/边界元方法可以很好地求解复杂结构的振声响应，但限于网格数量和分析频率相关^[29]而不适用于高频计算。统计能量法^[30-32]从统计学的角度求解结构的振声响应，要求模型有足够大的振动和声模态密度，适用于中高频计算，而对模态密度不够的中低频段则存在统计学上的不确定性。用统计能量法分析结构振声响应主要取决与3个参数的估计^[29]：1）单个子系统的模态密度；2）单个子系统的内部损耗因子（阻尼）；3）各个子系统间的耦合损耗因子（耦合程度）。

2 加筋/夹层板壳声振特性研究 2.1 板壳结构

普通板壳结构声振特性研究是加筋/夹层板壳结构声振特性研究的基础，从普通板壳到加筋/夹层板壳再到复合材料加筋/夹层板壳，模型所考虑的因素逐渐增多，但模型基本特征未变，所采用的研究方法也是不断继承和发展的。以经典的矩形简支薄板的声振计算为例 ^[33]，假设矩形薄板四周为无限大刚需障板，板受外部激励 $ F(x,y) $ 作用，横向位移为 $ w(x,y) $ ，矩形板的运动方程为：

$ D{\nabla ^4}w(x,y) - {\rho _D}h{\omega ^2}w(x,y) = F(x,y) - {p_0}(x,y) 。$

(1)

其中： $ D = {E_0}{h^3}/12(1 - {\nu ^2}) $ 为板的弯曲刚度； $ {F_{}}(x,y) = {F_0}\delta (x - {x_o}) \delta (y - {y_0}){e^{j\omega t}} $ 为法向点激励； $ {p_0}(x,y) $ 为板表面声压。当考虑板四周为无限大刚性障板时，声场声压可用Rayleigh积分公式求得：

$ {p_0}(x,y) = \frac{{i{\rho _0}\omega }}{{2{\text π} }}\int_{{S_p}} {\dot \omega (x',y')} \frac{{{e^{ - i{k_0}R}}}}{R}{\rm{d}}x'{\rm{d}}y'。$

(2)

其中： $ {\ \rho _0} $ 为流体密度； $ \dot \omega (x',y') $ 为板表面振速；R为板表面点与场点之间的距离。用模态叠加法求解声振耦合方程，将板的横向位移 $ w(x,y) $ 用双三角级数形式表示： $ w(x,y) = \displaystyle\sum\limits_{m = 1}^\infty {\sum\limits_{n = 1}^\infty {{a_{mn}}\sin ({k_m}x)\sin ({k_n}y)} } $ ，将声压表达式和位移表达式代入运动方程，利用模态正交性得到模态运动方程：

$ Z_{mn}^M{a_{mn}} + i\omega \sum\limits_{m = 1}^\infty {\sum\limits_{n = 1}^\infty {{Z_{mnpq}}{a_{mn}} = {F_{mn}}} } 。$

(3)

其中： $ Z_{mn}^M = {\rho _p}h(\omega _{mn}^2 - {\omega ^2}) $ 为无流体载荷时板的机械阻； $ {\omega _{mn}} = \sqrt {D/{\rho _h}} (k_m^2 + k_n^2) $ 为矩形板在真空中的模态特征频率； $ {Z_{mnpq}} $ 为板表面的声辐射阻抗。求解式（2）可得模型的位移模态，进一步计算可求得矩形板的均方振速、辐射声功率和辐射效率等指标。

2.2 加筋/夹层板壳结构

相较于普通板壳结构，加筋/夹层板壳通常有更好的力学和声学性能，在工程应用中极为广泛^[34]，在研究其声振特性时，理论模型通常会根据实际情况进行一定程度简化。例如，对于加筋板壳，可简单地将加强筋视为作用在面板上的线力^[16]，或者等效为弹簧的横向支撑和扭簧的旋转支撑^[2,18]，亦或者考虑弯曲和扭转的同时还考虑加强筋与板的耦合作用^[1]；对于夹层板壳，若夹芯层为非连续形式，如点阵结构、格栅结构、波纹结构等离散夹芯^[35]，还需要进行层芯受力简化。例如将波纹夹芯简化为具有平移刚度和扭转刚度的弹簧^[36-37]，并假设芯层的质量时均布于上下面板，或者将夹层进行均质化等效，得到等效的弹性参数^[5,38]。

在考虑加筋/夹层板壳结构声振特性时，研究对象通常有无限大周期结构和有限尺寸结构之分^[39]，这两者的关注重点不同，前者主要着眼于结构波在加筋/夹层板上的传播特性，如Mead和Pujara等^{[13-14,40-41]}研究了声波在周期弹性支撑梁和无限大周期加筋板中的传播问题，应用空间谐波法求解了周期加筋结构的声传递特性，并讨论不同加筋参数对振声响应的影响。此类研究偏向于揭示声振特性形成机理，并对结构进行波动分析^[22]。相对地，对有限尺寸结构则主要关注结构在特定激励作用下的振声响应，有限尺寸加筋/夹层板有明确的边界，可以建立起封闭的理论模型，求解得到特定频率下振动和声辐射的封闭解^[42-44,9]。

对于呈现周期排列的加筋/夹层结构，还可以通过周期结构Floquet-Bloch理论^[45]建立声振耦合方程。例如，Zhou^[25]在研究将正交加筋板的声辐射特性时，将加筋视为周期排列的十字交叉点进行分析，但这种方法用于分布稀疏周期加筋结构时会带来较大误差。对于周期加筋/夹层结构，还可以利用等效方法计算出胞元结构的等效性质^[5,38]，典型方法有代表体单元法（RVE）^[46]和渐进均匀化法（AH）^[47-49]等。

数值方法可以很好地解决复杂结构和复杂边界结构的声振计算问题^[50-52]，数值与解析相结合的方法，可以灵活适用不同结构且具有较高的计算效率。Laurent^[53]在研究水下不规矩加筋壳的振声响应时，将受流压作用的壳和肋骨框架视为不同的子系统并分别求解导纳，通过圆周导纳方法（CAA）将各子系统变成一个耦合系统，用解析和数值相结合的方法分析了水下耐压艇体在较宽频率范围内的振声性能。Yang^[15]提出一种可以预测周期结构振声响应的模型，将行波法和有限元法相结合，用有限元法建立周期结构的单元模型，求得单元的质量矩阵和刚度矩阵后代入解析模型进行分析，极大简化了声振计算过程。文献[54-57]也采用了类似的研究思路，利用有限元与等效法相结合简化计算模型。

2.3 复合材料加筋/夹层板壳结构

本文所讨论的复合材料特指纤维增强树脂基体复合材料，其最大的特点是具有分层特性和各向异性的材料属性，在研究结构声振特性时多基于复合材料力学和层合板理论考虑这一材料特性。Yin^[19]研究了双周期加筋复合材料层合板的声辐射特性，考虑加筋板受简谐点激励作用，利用层合板理论考虑复合材料面板的各向异性特性，同时利用空间谐波法求解辐射声压，最后分析加筋对不同频段声辐射的影响。Shen^[58]研究了无限大正交加筋复合材料层合板在声波和外力激励下声辐射特性，其中复合材料层合板横向位移采用一阶剪切的分层理论考虑，并过傅里叶变换法求解得到不同激励下的结构声辐射。Shen^[17]通过相似的方法研究了复合材料加筋层合板的声传递特性，并分析了不同截面形式加筋的影响。Mejdi^[59]基于模态叠加法研究了复合材料加筋板的声振特性，使用一阶剪切变形理论考虑复合材料面板和复合材料加强筋的横向变形，并将加强筋等效为作用于面板的弯矩和剪力，利用连续关系集成在面板的阻抗表达式中，分析了3种截面形式加强筋的不同影响，以及不规则加筋板和周期性加筋板的振声响应。与有限元法/边界元法计算进行比较，结果具有良好的一致性。Ghinet^[60-61]基于分层理论提出了一种通用的复合材料夹层板壳理论模型，其中每个分层采用一阶剪切理论，考虑分层的弯曲变形、横向剪切变形以及转动惯量和各向异性铺层角影响，最后用行波法求解了无限大夹层板在散射场作用下的传声损失，并用有限元法和谱有限元法分别验证了模型的有效性。Cherif^[62]对该模型进行了详细的实验验证，对结构波数、模态密度、阻尼损耗因子、辐射效率、传声损失等指标进行比较，验证了该理论模型预测复合材料夹层板壳声振特性的准确性。

2.4 小　结

复合材料加筋/夹层板壳结构声振特性研究所考虑的因素很多，采用的研究方法是从普通板壳、加筋/夹层板壳结构声振研究的方法中继承和发展的。在研究加筋/夹层板壳结构声振特性时，通常需要进行一定的模型简化，或者使用等效参数和周期理论。在各种研究方法中，空间谐波法和傅里叶变换法适合无限大周期结构的声振分析，模态叠加法适用于有限尺寸结构的声振响应计算，行波法可以用来分析结构的频散特性，带隙特性^[22]。数值法和解析法相结合进行声振分析，既能解决结构外形复杂难以解析求解的问题，也能避免计算资源消耗过大的问题。

3 新材料应用与声振特性研究

新材料的应用给复合材料加筋/夹层结构声振研究带来了新的发展^[63]。例如，声学超材料被用来改善加筋/夹层板壳结构的声传递特性^[64]，智能材料^[65-66]被用来实现加筋/夹层结构的振动主动控制，这些材料的应用使声振计算变得更复杂，也使得加筋/夹层结构有了更大的设计空间。

声学超材料是一类具备超常物理特性的人工复合材料，具有如负等效质量密度、负等效弹性模量等特性，声学超材料的概念最早是Liu等^[67-68]在研究局部共振声学材料时提出的，通过对材料组成和结构构型进行设计，可以实现小尺寸结构对大波长声波的传播抑制。从而实现对声波的调谐^[69]，这给声振特性研究带来新的变量，增大了结构声学的设计空间。Chen^[70]研究了声学薄膜板的低频吸声能力，通过建立薄膜板的声振耦合理论模型，揭示了平面入射波作用下的声能吸收机理。不同于一般声传递特性研究，局域共振（薄膜上附着质量的共振和反共振运动）是声学超材料产生低频吸声能力的关键，Chen通过求解声振耦合的积分-微分方程，准确地对结构吸声能力进行定量测定，其结果与有限元预测结果吻合较好。

智能材料也是研究的热点之一，所谓智能材料的定义相当模糊，通常这些材料都具有一些常规材料所不具备的材料特性，且有一定的可调节性，比如压电材料^[71-72]、温控阻尼材料^[73]、形状记忆合金^[74]等。就研究结构的声振特性而言，这些材料的应用主要是实现结构振动的主动控制，进而影响声振特性。文献[71-72]研究了如何利用压电纤维主动阻尼约束层来控制夹层结构振动，其中贴敷于层合梁上的约束阻尼分为压电材料层和粘弹性阻尼层2层，通过对压电纤维层实施适当的电压控制，可以主动调节粘弹性阻尼材料层的横向剪切变形，增大对结构的阻尼约束作用，从而抑制夹层梁振动，更多地耗散结构动能。文献[73]利用高分子聚合材料在玻璃化转变区具有更高的阻尼性能这一特性，提出一种温控变阻尼技术，通过主动调节芯层温度控制其阻尼性能，从而实现对夹层结构振动的控制，实验表明温控变阻尼技术抑振效果明显。

4 边界条件与声振特性研究

实际工程结构的边界一般都是非理想状态的，即结构的边界是一种弹性状态，对于有限尺寸复合材料加筋/夹层板壳，边界对其声振特性的影响是无法回避的，因而研究任意边界条件下的声振特性具有重要意义。目前对复合材料板壳结构的声振特性研究，很多都是以简支为边界条件的，这是因为在简支边界下更容易获得解析解^[75]，还有不少文献虽然也进行任意边界条件下的研究，但实质却是不同经典边界的任意组合^[76]。

研究任意边界条件下复合材料加筋/夹层板壳结构的声振特性，既考虑复合材料的各向异性同时还要考虑结构的弹性边界，这使得声振求解变得异常困难，因而大部分文献未能考虑所有条件进行综合分析，而是退而求其次地研究结构的自由振动特性^[77-78]，或结合数值方法和实验方法来进行研究。

自由振动求解难度低，且能反映结构固有频率随质量、刚度分布变化的情况，是声振特性研究的基础。Qin^[79]研究了圆柱壳在任意边界条件下自由振动，通过设置不同刚度的支撑弹簧来模拟任意边界条件，采用Sanders壳理论计算壳体弹性应变能，并使用Rayleigh-Ritz法推导控制方程，最后分析了圆柱壳在任意条件下的自由振动特性。Berry^[80]研究了矩形薄板在任意边界条件下的振动和声辐射特性，同样是利用Rayleigh-Ritz法推导获得矩形板的振动控制方程，并利用Rayleigh积分公式进行远场辐射声压计算，最后得到结论：矩形板受到的边界约束越小，板的辐射效率就越低。文献[79,81-84]均为研究各类板壳结构在任意边界下的振动，其中既有夹层结构、加筋结构、格栅结构等不同结构形式的，也有功能梯度材料，粘弹性材料等不同材料特性，但都是通过Rayleigh-Ritz法和Hamilton原理来进行振动分析的，这样可以避开复杂结构和边界带来的复杂力学分析，较为简便地得到振动控制方程。同时，上述文献在处理板壳结构的任意边界时，多是参考结构的弯曲刚度来表示边界约束条件，这样可以保证最后得到振动控制方程可解，当板壳结构是具有各向异性的复合材料层合结构时，由于耦合刚度的存在，不便直接利用弯曲刚度的表示边界约束条件，因而不少文献中都是直接假设一个刚度参数^[85-86]来模拟边界条件。

用有限元方法特别是利用商用有限软件模拟任意边界条件是非常方便的^[87]，不少文献中都将有限元计算作为理论计算的验证，但和理论分析一样，边界的刚度约束虽然可以方便地假定为任意数值，但难以合理地模拟和还原实际结构的边界约束条件，这方面的研究需要针对性的试验作为支撑。

5 结　语

随着复合材料加筋/夹层结构被不断地应用于航空航天、船舶、汽车、高速列车等不同领域，对其声振特性的研究也逐步深入。复合材料加筋/夹层板壳结构声振问题较为复杂，需要考虑的因素较多，本文从介绍声振特性研究内容和研究方法开始，先后归纳总结了加筋/夹层板壳结构声振特性研究、考虑新材料应用的声振特性研究、任意边界条件下的声振特性研究的相关文献，通过综述上述研究文献，可以发现：

1）解析与数值相结合的研究方法有更好的应用灵活性和计算效率，是研究复合材料加筋/夹层结构声振特性的有效方法，也是重要的发展方向；

2）新材料的应用为复合材料加筋/夹层结构的声振特性研究带来了新变化，增大了结构动力学设计和声学设计的设计空间；

3）现有的分析方法可以模拟出任意边界条件而无法还原出实际边界条件，因而需要结合试验研究作为支撑。