周海安，王曉明，宋鳳敏，梅玉林

（1山東理工大學(xué)機械工程學(xué)院，山東淄博255049；2大連理工大學(xué)a.機械工程學(xué)院；b.汽車工程學(xué)院，遼寧大連116024；3山東水利職業(yè)學(xué)院，山東日照276826）

諧振面力作用下的無窮大雙周期加筋板的振動響應(yīng)分析

周海安1，王曉明2a，宋鳳敏3，梅玉林2b

（1山東理工大學(xué)機械工程學(xué)院，山東淄博255049；2大連理工大學(xué)a.機械工程學(xué)院；b.汽車工程學(xué)院，遼寧大連116024；3山東水利職業(yè)學(xué)院，山東日照276826）

文章主要研究了流體負載下的無窮大雙周期加筋薄板，在周期諧振力作用下的振動響應(yīng)。假設(shè)薄板和加強筋只存在法向力的作用，通過薄板和加強筋的位移連續(xù)條件，以及格林函數(shù)及傅立葉變換公式的應(yīng)用，建立了頻域內(nèi)雙周期加筋薄板的振動位移方程。利用空間波數(shù)法，將加筋薄板的振動位移表達為波數(shù)分量的迭加形式。通過數(shù)值計算的方法對波數(shù)分量進行求解，最后通過方程的傅立葉逆變換，得到了加筋薄板的振動方程。數(shù)值計算的結(jié)果表明文中的計算方法收斂速度快。

振動響應(yīng)；加筋板；格林函數(shù)；空間波數(shù)法

0 引言

加筋板是許多工程結(jié)構(gòu)的組成部分，在工程界得到了廣泛的應(yīng)用，尤其在機械，船舶和航天領(lǐng)域，加筋板的應(yīng)用更為常見。到目前為止，許多學(xué)者對加筋板的振動特性進行過研究，并且提出了許多不同的求解方法[1-8]。然而大部分的科研工作主要集中在有限尺寸或單周期無窮大加筋板的振動研究上，而對于水下無窮大雙周期加筋薄板振動響應(yīng)的研究則要少得多。本文的主要目的是尋求一種合適的數(shù)學(xué)方法以便對雙周期加筋薄板結(jié)構(gòu)進行振動分析。但是由于加筋薄板的周期性以及流固耦合問題的復(fù)雜性，幾種重要的數(shù)學(xué)方法如Rayleigh-Ritz能量法﹑有限差分法﹑傳遞矩陣法﹑模態(tài)法等并不能很好地解決無窮大周期結(jié)構(gòu)的振動問題。由Mead和Pujara[9]提出的空間波數(shù)法在處理流固耦合的周期結(jié)構(gòu)振動問題上體現(xiàn)出了明顯的優(yōu)勢：Mace[1]在空間波數(shù)法的基礎(chǔ)上，推導(dǎo)出流體作用下雙周期加筋薄板的振動和輻射聲壓方程，并且分析了點力作用下的水中遠場輻射聲壓。Lee和Kim[10]等亦通過空間波數(shù)法對單周期加筋薄板的振動響應(yīng)及聲透射系數(shù)進行了理論分析與計算。此外，傅里葉變換的方法也是研究結(jié)構(gòu)振動和聲輻射的一種有力工具。吳文偉等[11]曾經(jīng)采用傅里葉變換的方法對單周期加筋薄板的聲輻射進行過研究。

本文對水下無窮大雙周期加筋薄板的振動響應(yīng)進行了分析，并且根據(jù)格林函數(shù)與空間波數(shù)法以及傅里葉變換法，建立了加筋薄板的振動方程，并最終通過數(shù)值計算的方法得到了方程的解。在計算結(jié)果中，本文重點分析了加強筋橫截面尺寸以及周期間距等幾個主要物理參數(shù)對薄板振動響應(yīng)帶來的影響。

1 加筋薄板的振動方程

1.1 數(shù)學(xué)模型

無窮大雙周期加筋薄板的上表面空間充滿水流體，薄板的下表面空間是空氣。在x軸方向的加強筋沿x=mLx周期排列，在y軸方向的加強筋沿y=nLy（m，n=0，±1，…，±∞）周期排列，周期加筋薄板的局部結(jié)構(gòu)示意圖如圖1所示。

圖1 雙周期加筋薄板Fig.1 Diagram of the two-dimensional beam-stiffened plate

雙周期加筋薄板受到隨時間t周期變化的諧振面力f（x，y）exp（jωt）的作用，其方程表達式是[1]

pe是諧振力的幅值，ω是圓頻率，k0是諧振力波數(shù)，θ是諧振力與z軸方向的夾角，φ是諧振力在水平面上的分量與x軸方向的夾角。

對雙周期加筋薄板進行振動響應(yīng)分析，假設(shè)加筋薄板的振動位移是w（x，y）exp（jωt），加強筋與板的相互作用力僅考慮法向力的作用，則周期加筋板的振動方程是（為了公式的簡潔，在以下所有方程中都將省略時間因子exp（jωt））

式中：G（x-ξ，y-η）是格林函數(shù)，其物理含義是：薄板由于受到作用在（ξ，η）點的單位激勵力而在（x，y）位置處的振動位移。F1（ξ，η）δ（ξ-mLx）和F2（ξ，η）δ（η-nLy）是加強筋沿直線ξ=mLx和η=nLy作用于薄板的法向力。這里假設(shè)兩組加強筋是完全相同的，利用薄板和加強筋的位移連續(xù)條件可以得到

式中：E1，I，ρb，Sb分別是加強筋的楊氏模量，慣性矩，質(zhì)量密度和橫截面面積。

1.2 振動方程的傅里葉變換

要求解方程（3），需要對方程進行傅里葉變換。對任意滿足要求的函數(shù)g（x，y），其傅里葉變換是

傅里葉逆變換是

對于滿足傅里葉積分定理的任意函數(shù)g1（x，y）和g2（x，y），根據(jù)函數(shù)卷積的傅里葉變換等于函數(shù)的傅里葉變換的乘積的性質(zhì)可知

將（4）式代入（3）式中，對（3）式進行傅里葉變換，并根據(jù)（7）式的積分性質(zhì)可得

式中：w*（α，nLy），w*（mLx，β）分別是w（x，nLy），w（mLx，y）關(guān)于x軸和y軸的單重傅里葉變換，并且

流體負載下薄板的格林函數(shù)滿足方程

式中：c是水流體中的聲速。

薄板和水流體接觸面的位移邊界條件是

對（12）式和（13）式進行傅里葉變換，求解得到

對（11）式進行傅里葉變換，結(jié)合（14）式可以得到頻域內(nèi)流固耦合的薄板阻抗

2 振動方程的求解

按照空間波數(shù)法理論[9]可知，二維周期加筋薄板的振動位移可以表達為如下形式[1]

式中：Wmn是（m，n）th諧振分量振幅，且kx，mx=kx+2mπ/Lx，ky，n=ky+2nπ/Ly。

由（16）式可以得到振動位移的周期條件

根據(jù)（17）式加筋薄板的周期條件，可知

將（18）式代入（8）式中，可得

雙周期加筋薄板在原點處的振動位移分別可以表達為

根據(jù)（16）式和（20）式可知，加筋薄板沿x=0和y=0直線上的振動位移分別是

對（21a）式和（21b）式分別進行關(guān)于y軸和x軸的單重傅立葉變換，得到：

由泊松迭加公式可知

同理可得

將（9）、（22a）、（22b）式和（23a）、（23b）式代入（19）式中，可以得到頻域內(nèi)加筋薄板關(guān)于未知量Vn和Nm的振動方程

對（24）式進行關(guān)于α的單重傅立葉逆變換，并且令變量x=0以及β=ky，n，根據(jù)（22a）式得到方程[1]

同理，對（24）式進行關(guān)于β的單重傅立葉逆變換，并令變量y=0以及α=kx，m，根據(jù)（22b）式得到方程?

對（24）式進行傅立葉逆變換，最終得到時域內(nèi)加筋薄板的振動方程

本文提出了無窮大雙周期加筋板振動響應(yīng)的格林函數(shù)法，通過比較可以發(fā)現(xiàn)，雖然本文與文獻[1]采用的方法不同，但是本文最終求解的振動位移表達式（方程（28））和文獻[1]中的表達式（方程（34））和實際物理含義是相同的，從而證明了本文方法的正確性。相對于文獻[1]本文方法的優(yōu)點在于：格林函數(shù)的方法不僅可以用來計算諧振面力的作用，而且可以計算包括點力激勵在內(nèi)的任何作用力形式；通過格林函數(shù)的應(yīng)用，得到了頻域內(nèi)流固耦合的薄板阻抗，揭示了加強筋作用力與薄板振動位移的周期關(guān)系（（8）式）；本文中將加強筋的作用力表達為薄板振動位移的分量形式，其求解方法也與文獻[1]不同。

3 計算分析

下面將利用推導(dǎo)的公式來計算薄板和加筋薄板的振動，并且分析加強筋對薄板振動的影響。計算中用到的部分物理參數(shù)如下：E=E1=2.01×1011×（1+0.02j），v=0.2，h=0.005，ρ=ρb=7 850，ρ0=1 000，c=1 500，pe=1 000，θ=φ=π/4，k0=ω/c。加強筋的橫截面面積是Sb=0.005×0.05，短邊與薄板相接。所有的單位均采用國際標(biāo)準(zhǔn)單位。本文在下面的算例分析中，通過計算發(fā)現(xiàn)當(dāng)Q=10時，已經(jīng)得到滿足收斂性要求的計算精度。

圖2是薄板在頻率為2 000 Hz諧振力作用下，角度分別為φ=π/4和φ=π/6，在-1≤x，y≤1且時間t=0的振動響應(yīng)。

圖2 薄板的振動響應(yīng)Fig.2 The vibration responses of the unstiffened plate

薄板的振動特性主要取決于它的振動相位角和幅值大小。薄板的振動相位角變化方向由φ決定。從圖中可以看出，當(dāng)φ=π/4時，振動相位沿與x軸成45°的直線方向變化；如果改變φ的大小，如圖所示，振動相位角變化方向?qū)S之改變。與φ不同，θ并不改變相位角變化方向，但是可以改變薄板的振動幅值大小。

圖3是不同周期間距的加筋薄板在頻率為2 000 Hz諧振力作用下，在-1≤x，y≤1且時間t=0的振動響應(yīng)。

圖3 加筋薄板的振動響應(yīng)Fig.3 The vibration responses of the beam-stiffened plate

對比圖2和圖3可以發(fā)現(xiàn)，相對于薄板結(jié)構(gòu)，雙周期加筋薄板的振動具有二維周期性質(zhì)，而且加筋薄板的振幅變化比較明顯。從圖3（a）中可以看出，加筋薄板的振動沿x軸和y軸方向上是完全一致的，這是由于沿兩個軸方向的周期間距相等而且壓力角φ=π/4。如果改變兩軸的周期間距或壓力角φ，如圖3（b）所示，加筋薄板沿x軸和y軸方向的振動將不再具有對稱性。

圖4給出了薄板和不同周期間距加筋薄板在諧振力頻率為50～2 500 Hz范圍內(nèi)分別在（0，0）點和（0.4，0.4）點處的振幅。周期間距分別是Lx=Ly=0.2和Lx=Ly=0.4。

從圖4中可以發(fā)現(xiàn)，薄板在（0，0）點的振幅幾乎為0，且隨頻率的增加不發(fā)生改變，表明在薄板中傳播的振動彎曲波在通過原點位置時受到阻抑作用；在（0.4，0.4）點薄板的振幅隨頻率的增加而變化，振動波的傳播可以順利通過。加強筋對薄板振動的影響在100～1 300 Hz赫茲范圍內(nèi)作用明顯，相對于薄板，加筋薄板的振幅振蕩較大。隨著頻率的進一步增加，從整體分析，加強筋的作用減弱，加筋薄板振幅的振蕩峰值減小，振幅的變化隨頻率的增加趨于平緩。當(dāng)周期間距增大后，加筋薄板在更多的頻率位置出現(xiàn)共振現(xiàn)象，表明適當(dāng)增大周期間距可以增強加筋薄板的振蕩效果，從而在較多的頻率范圍內(nèi)影響薄板的振動。而且由于加強筋周期間距的增大，使得加筋薄板的共振頻率向低頻范圍偏移，因此可以通過調(diào)節(jié)周期間距的方法改變加筋薄板的共振頻率。

圖4 在不同位置點的振動響應(yīng)Fig.4 The vibration responses at different points

圖5 不同加強筋對振幅的影響Fig.5 The influence of different beam-stiffeners on the vibration amplitudes of the stiffened plate

為了分析不同尺寸的加強筋對薄板的影響，圖5給出了加強筋橫截面面積分別為0.01×0.1和0.01× 0.15的周期加筋薄板在（0.4，0.4）點隨頻率變化的振幅曲線。

對比圖5（a）和5（b）發(fā)現(xiàn)，增加加強筋的橫截面高度對加筋薄板在200 Hz頻率處的振幅影響比較明顯，但對其它頻率振幅的影響相對較小，而且振幅曲線在400 Hz之后隨頻率的增加其變化趨勢是相似的。當(dāng)同時改變加強筋橫截面的寬度和高度尺寸后，通過對比圖4（b）和5（a）可以發(fā)現(xiàn)，周期間距為Lx=Ly=0.4的加筋薄板的振幅曲線發(fā)生變化。從總體上分析，加強筋橫截面尺寸的變化對1 000 Hz以內(nèi)振動響應(yīng)的作用較為顯著，而對較高頻率振動響應(yīng)的影響較弱。

4 結(jié)論

根據(jù)格林函數(shù)和空間波數(shù)法以及傅立葉變換的應(yīng)用，本文提出了一種求解受諧振力作用的流體加載下的無窮大雙周期加筋薄板振動方程的方法。通過將加筋薄板的振動位移展開為波數(shù)分量迭加的形式，并采用數(shù)值計算的方法求解分量的值，從而得到了加筋薄板的振動方程表達式。在計算結(jié)果中，本文詳細分析了加強筋對薄板振動響應(yīng)的作用，結(jié)果表明周期加強筋對薄板的低頻振動影響較大。

參考文獻：

[1]Mace B R.Sound radiation from fluid loaded orthogonally stiffened plates[J].Journal of Sound and Vibration,1981,71(3): 435-441.

[2]Alain Berry.A new formulation for the vibrations and sound radiation of fluid-loaded plates with elastic boundary conditions[J].Journal of the Acoustical Society of America,1994,96(2):889-901.

[3]洪明,孫盈盈.船上含敷料板與加筋結(jié)構(gòu)的固有振動特性研究[J].船舶力學(xué),2008,12(2):311-318. Hong Ming,Sun Yingying.Natural vibration characteristics of stiffened plates and plates with covering layer in ships[J]. Journal of Ship Mechanics,2008,12(2):311-318.

[4]Lin T R,Pan J.A closed form solution for the dynamic response of finite ribbed plates[J].Journal of the Acoustical Society of America,2006,119(2):917-925.

[5]黎勝,楊婧媛.水下加筋板振動聲輻射的代理模型研究[J].聲學(xué)學(xué)報,2010,35(6):659-664. Li Sheng,Yang Jingyuan.Research on surrogate models for structural vibration and acoustic radiation of underwater stiffened plates[J].Acta Acustica,2010,35(6):659-664.

[6]Maxit L.Wavenumber space and physical space responses of a periodically ribbed plate to a point drive:A discrete approach[J].Applied Acoustics,2009,70(4):563-578.

[7]Mace B R.Sound radiation from a plate reinforced by two sets of parallel stiffeners[J].Journal of Sound and Vibration, 1980,71(3):435-441.

[8]陳海龍,金葉青,龐福振.一種分析單向加筋板隔聲性能的解析—數(shù)值耦合法[J].船舶力學(xué),2013,17(8):952-959. Chen Hailong,Jin Yeqing,pang Fuzhen.Sound insulation performance analysis of single stiffened plate by a combined analytical-numerical method[J].Journal of Ship Mechanics,2013,17(8):952-959.

[9]Mead D J,Pujara K K.Space harmonic analysis of periodically supported beams:Response to convected random loading [J].Journal of Sound and Vibration,1971,14(2):525-532.

[10]Lee J H,Kim J.Analysis of sound transmission through periodically stiffened panels by space-harmonic expansion method [J].Journal of Sound and Vibration,2002,251(3):349-366.

[11]吳文偉,冷文浩,沈順根.具有等間距相同加強筋板的聲輻射[J].中國造船,1999,146(3):72-81. Wu Wenwei,Leng Wenhao,Shen Shungen.Sound radiation of plate with equidistant reinforced ribs[J].Shipbuilding of China,1999,146(3):72-81.

Investigation of the vibration response of an infinite two-dimensional rib-stiffened plate under a harmonic plane pressure

ZHOU Hai-an1,WANG Xiao-ming2a,SONG Feng-min3,MEI Yu-lin2b
(1.School of Mechanical Engineering,Shandong University of Technology,Zibo 255049,China;2a.School of Mechanical Engineering;b.School of Automotive Engineering,Dalian University of Technology,Dalian 116024,China;3.Shandong Water Polytechnic,Rizhao 276826,China)

The vibration response of an infinite fluid-loaded plate,reinforced with two sets of orthogonal rib stiffeners,is investigated.The stiffeners are assumed to exert only forces on the plate.With the displacement continuity condition of the plate and the stiffeners,the vibration response of the rib-stiffened plate to a harmonic pressure in the wavenumber space is established by using the Green function and Fourier transforms.Applying the space harmonic method,the vibration displacement of the stiffened plate is then expressed as infinite sets of space harmonic amplitudes.After having a numerical solution for these amplitudes,the Fourier inverse transform is employed to find an expression for the vibration response of the stiffened plate in the physical space.Numerical results show that this method converges rapidly.

vibration response;rib-stiffened plate;Green function;space harmonic method

U661.4

A

10.3969/j.issn.1007-7294.2015.11.008

1007-7294（2015）11-1352-08

2015-07-15

國家自然科學(xué)基金資助（51405276，50875030）

周海安（1981-)，男，博士，講師，通訊作者，E-mail：zhouhaian1981@163.com；

王曉明（1965-），男，教授，博士生導(dǎo)師，E-mail：xiaoming@dlut.edu.cn；

宋鳳敏（1978-），女，講師，E-mail：stj08@163.com；

梅玉林（1969-），女，副教授，E-mail：meiyulin@dlut.edu.cn。