王 晴 洪 明 周艷秋

1中國艦船研究設(shè)計(jì)中心，湖北 武漢 430064

2大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點(diǎn)實(shí)驗(yàn)室運(yùn)載工程與力學(xué)學(xué)部船舶工程學(xué)院，遼寧大連 116024

3中國艦船研究設(shè)計(jì)中心，上海 201108

夾層結(jié)構(gòu)模型假設(shè)對振動(dòng)和聲傳輸特性影響研究

王 晴1洪 明2周艷秋3

1中國艦船研究設(shè)計(jì)中心，湖北 武漢 430064

2大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國家重點(diǎn)實(shí)驗(yàn)室運(yùn)載工程與力學(xué)學(xué)部船舶工程學(xué)院，遼寧大連 116024

3中國艦船研究設(shè)計(jì)中心，上海 201108

不同模型假定下夾層結(jié)構(gòu)的聲振特性各不相同?；赗eissner模型和Zig－Zag模型位移場假設(shè)，分別建立了兩種復(fù)合材料夾層板有限元?jiǎng)恿卧Ｐ?，并通過推導(dǎo)的模型編制了計(jì)算程序，采用子空間迭代法求得了其固有頻率和固有振型，采用混響聲場聲傳輸理論計(jì)算了傳聲損失，并對夾層板的動(dòng)力特性和聲傳輸特性進(jìn)行了數(shù)值模擬研究，同時(shí)還考慮了夾層結(jié)構(gòu)面板與芯板的厚度比和彈性模量比不同時(shí)各種模型的適用范圍。

Reissner模型；Zig－Zag模型；Mindlin理論；振動(dòng)特性；傳聲損失

1 引言

隨著各個(gè)工程領(lǐng)域?qū)Σ牧弦蟮牟粩嗵岣撸瑥?fù)合材料夾層板也得到了廣泛的應(yīng)用。尤其是在船舶設(shè)計(jì)中，開始越來越多地應(yīng)用夾層結(jié)構(gòu)。通常，設(shè)計(jì)者們均十分關(guān)注船舶的居住性和舒適性，而這又與其結(jié)構(gòu)振動(dòng)狀況和隔聲狀況密切相關(guān)，因此，對夾層結(jié)構(gòu)的振動(dòng)特性和聲傳輸特性研究具有實(shí)際意義。

對于夾層板問題的研究，大多是以上世紀(jì)40～50年 代 建 立 起 來 的 Reissner、Hoff和∏pycakob 理論為基礎(chǔ)［1－3］，這 3 種理論對求解均勻材料面板的夾層結(jié)構(gòu)比較有效。

1）Reissner理論

該理論認(rèn)為面板很薄，假設(shè)其主要承受面內(nèi)應(yīng)力為 σx、σy和 τxy，且都沿板厚度均勻分布。 在芯層內(nèi)，σx＝σy＝τxy＝0， 芯層可抵抗橫向剪切力，并考慮了芯層橫向剪切變形的影響。當(dāng)夾層板關(guān)于中面做反對稱彎曲變形時(shí)，假設(shè)面板和芯層內(nèi)都有εz＝0。在面板和芯層內(nèi)，都不計(jì)σz對形變的影響。

該理論把面板作為只承受面力的薄膜處理，忽略了本身的抗彎剛度，而夾芯則假設(shè)其僅具有抗橫向剪切能力，這是一種最簡單的夾層板理論。

2） Hoff理論

該理論認(rèn)為表層面板很薄，只把它們看成是通常的彈性薄板。 在芯層內(nèi)，σx＝σy＝τxy＝0，考慮了芯層橫向剪切變形的影響。當(dāng)夾層板對中面作反對稱變形時(shí)，假定芯層εz＝0，忽略了σz對形變的影響。

該理論將面板視為普通薄板，即只考慮了其面內(nèi)應(yīng)力和抗彎剛度，而對于夾芯部分，則仍認(rèn)為其僅具有抗橫向剪切能力。該理論和Reissner理論都沒有考慮夾芯的其他變形能力。

夾層板在轉(zhuǎn)動(dòng)時(shí)，面板和芯層的轉(zhuǎn)角是不一致的，其位移具有齒形特征，尤其是當(dāng)厚度較大或材料較軟時(shí)，忽略轉(zhuǎn)動(dòng)項(xiàng)將會(huì)帶來較大誤差［4］。本文推導(dǎo)了忽略面板轉(zhuǎn)動(dòng)項(xiàng)的Reissner理論有限單元模型和考慮面板轉(zhuǎn)動(dòng)項(xiàng)的基于Hoff理論基礎(chǔ)上Srinivas提出的Zig－Zag模型。計(jì)算了夾層板的面板和芯板厚度比逐漸增加和彈性模量比逐漸增加時(shí)，結(jié)構(gòu)自由振動(dòng)固有頻率和混響聲場入射情況下夾層板傳聲損失的變化情況。

2 夾層板自由振動(dòng)有限元分析理論

2.1 模型示意圖、自由度及位移場假設(shè)

夾層板由上下面板和中間芯層組成。以夾層板中性面為坐標(biāo)平面，建立整體坐標(biāo)系，其形式如圖1所示。

基于Reissner模型假定和Mindlin一階剪切理論，單元每個(gè)節(jié)點(diǎn)有5個(gè)自由度，即沿3個(gè)方向的位移 u0、ν0、w0，以及繞 2 個(gè)軸的轉(zhuǎn)角 θx和 θy，夾層板單元的位移場假設(shè)為：

基于Hoff理論基礎(chǔ)上的Zig－Zag模型假定和Mindlin一階剪切理論，單元每個(gè)節(jié)點(diǎn)有9個(gè)自由度：即芯板的 5 個(gè)自由度 u0、ν0、w0、θx和 θy，以及上下面板各兩個(gè)自由度 φx、φy和 φx、φy， 夾層板單元各層的位移場假設(shè)為：

2.2 兩種模型夾層板結(jié)構(gòu)振動(dòng)有限元推導(dǎo)

對于夾層板自由振動(dòng)有限元理論，其與普通板的區(qū)別在于，由于各層的材料參數(shù)和幾何參數(shù)不同，其質(zhì)量陣和剛度陣需進(jìn)行疊加。

4節(jié)點(diǎn)夾層板的剛度陣和質(zhì)量陣的形成步驟為：按照位移插值函數(shù)和材料本構(gòu)關(guān)系，可得到形函數(shù)矩陣［N］、密度矩陣［ρ］和幾何矩陣［B］。

對于Reissner模型，其彈性矩陣可將芯板彈性 陣［Dc］、上面板彈性陣［Du］和下面板彈性陣［Dd］進(jìn)行疊加，從而得到夾層板總的彈性矩陣［5］：

進(jìn)而，能夠得到夾層板總的剛度矩陣：

對于Zig－Zag模型，上下面板的幾何矩陣［B］也不相同，需要每一層形成各自的剛度矩陣之后再進(jìn)行疊加［6］。

質(zhì)量陣的形成方法與剛度陣相同。

設(shè)在結(jié)構(gòu)上作用有激勵(lì)｛f（t）｝，結(jié)構(gòu)產(chǎn)生的響應(yīng)為｛x（t）｝，由 Hamilton 原理，在動(dòng)載荷作用下，結(jié)構(gòu)的振動(dòng)方程可以表示為：

式中，［M］、［C］、［K］分別為結(jié)構(gòu)的質(zhì)量矩陣、阻尼矩陣和剛度矩陣。求解對應(yīng)的齊次方程，可得到該結(jié)構(gòu)相應(yīng)的固有頻率 ωj和振型｛φ｝j。

2.3 兩種模型夾層板結(jié)構(gòu)傳聲損失模型

選取混響聲場模型為比較理想的單頻模型，其假設(shè)為，是在各個(gè)方向上聲壓的幅值和相位隨機(jī)。這是一種很接近真實(shí)的混響聲場的模型，具體表示為［7－8］：

因此，混響聲場的入射聲功率可以認(rèn)為是I×J列平面聲波的入射聲功率相疊加。若空氣中板的長和寬分別為a和b，以角度（θ，φ）入射的平面聲波的入射聲功率∏i為：

式中，ρ為輻射介質(zhì)的密度；c為聲音在介質(zhì)中的傳播速度。

夾層板在平面聲波作用下的傳聲損失TL可以定義為：

式中，∏為結(jié)構(gòu)在平面聲波激勵(lì)下的輻射聲功率。

計(jì)算結(jié)構(gòu)輻射聲功率的方法很多，本文采用的是基于結(jié)構(gòu)表面振速計(jì)算輻射聲功率。

式中，Re（）表示取實(shí)部；νn*（P）表示表面 νn（P）的共軛復(fù)數(shù)。

3 數(shù)值模擬

下面，將選取典型算例來驗(yàn)證模型的適用性問題。四邊剛固夾層板長2 m、寬1.5 m，芯層材料參數(shù) E2＝ 8.96 × 108Pa， μ2＝ 0.42， ρ2＝ 999 kg／m3，厚1 cm。上下表層為厚度相等的鋼板，分別為1 mm、3 mm、6 mm 和 9 mm，E1＝2.1×1011Pa，μ1＝0.3， ρ1＝7 800 kg／m3，夾層板分析網(wǎng)格為8×6 個(gè)，單元?jiǎng)澐趾凸?jié)點(diǎn)編號(hào)模型如圖2所示。

上下面板厚度為1 mm、3 mm、6 mm和9 mm的夾層板前四階固有頻率對比如圖3所示。

在混響聲場激勵(lì)下，夾層板的傳聲損失對比如圖4所示。

若夾層板邊界條件和板的長寬不變，芯層彈性 模 量 分 別 取 E2＝8.96 ×106Pa、8.96 ×107Pa、8.96 ×108Pa 和 8.96 ×109Pa， μ2＝0.42， ρ2＝ 999 kg／m3， 厚度為 1 cm。 上下表層為鋼板，E1＝ 2.1 ×1011Pa， μ1＝ 0.3， ρ1＝7 800 kg／m3， 厚度相等，分別為1 mm。單元?jiǎng)澐趾凸?jié)點(diǎn)編號(hào)不變。

芯板的彈性模量為 8.96 ×106Pa，8.96 ×107Pa，8.96×108Pa和 8.96×109Pa的夾層板前 4 階固有頻率對比如圖5所示。

在混響聲場激勵(lì)下，夾層板的傳聲損失對比如圖6所示。

4 結(jié)論

經(jīng)過對2種模型的計(jì)算分析，得出以下結(jié)論：

1）在2種模型夾層板面板厚度逐漸增加的情況下，結(jié)構(gòu)振動(dòng)的固有頻率在面板厚度不足芯板厚度的1／10時(shí)幾乎相同。隨著厚度的增加，2種模型所得結(jié)果的差距會(huì)逐漸增大。在混響聲場入射情況下，夾層板的傳聲損失整體保持不變，但傳聲損失低谷產(chǎn)生的位置由于固有頻率不同，也有差別。

2）在2種模型夾層板芯板彈性模量逐漸增加的情況下，結(jié)構(gòu)振動(dòng)的固有頻率在芯板的彈性模量不足面板的0.4%時(shí)相差很大。隨著芯板彈性模量的增加，2種模型所得到的結(jié)果趨于相等。在混響聲場入射情況下，夾層板的傳聲損失整體保持不變，但傳聲損失低谷產(chǎn)生的位置由于固有頻率不同，也有差別。

3）模型化的差異對結(jié)構(gòu)振動(dòng)的固有頻率影響很大。在工程計(jì)算中，當(dāng)面板厚度與芯板厚度相比較大或面板剛度與芯板剛度相比較大時(shí)，Reissner理論已不能滿足計(jì)算精度的要求，而應(yīng)選擇一些精度更高的單元。但在面板厚度較小或彈性模量相差不大時(shí)，Reissner理論的計(jì)算方法由于節(jié)省計(jì)算空間，也具有其自身的優(yōu)越性。

［1］白瑞祥，陳浩然，蘇長健.考慮面板和芯體剪切效應(yīng)的復(fù)合材料夾層板穩(wěn)定性的有限元分析［J］.玻璃鋼／復(fù)合材料，1999（5）：3－6.

［2］師俊平，劉協(xié)會(huì)，越居才，等.復(fù)合材料夾層板的振動(dòng)及阻尼分析［J］.應(yīng)用力學(xué)學(xué)報(bào)，1996，13（2）：132－136.

［3］BANERJEE J R，SOBEY A J.Dynamic stiffness formulation and free vibration analysis of a three－layered sandwich beam ［J］.International Journal of Solid and Structures，2005，42（8）：2181－2197.

［4］CHO Y B，AVERILL R C.First－order zig－zag sublaminate plate theory and finite element model for laminated composite and sandwich panels［J］.Composite Structure，2000，50（1）：1－15.

［5］周艷秋.艦船基層結(jié)構(gòu)動(dòng)力響應(yīng)特性分析［D］.大連：大連理工大學(xué)，2005.

［6］白瑞祥，張志峰，陳浩然.基于Zig－Zag變形假定的復(fù)合材料夾層板的自由振動(dòng) ［J］.力學(xué)季刊，2004，25（12）：528－534.

［7］RAFAELY B： Spatial－temporal correlation of a diffuse sound field ［J］.Journal Acoustical Society of America，2000，107（6）：3254－3258.

［8］JACOBSEN F，ROISIN T.The coherence of reverberant sound field ［J］.Journal Acoustical Society of America，2000，108（1）：204－210.

Effects of Model Assumption on the Characteristics of Vibration and Sound Transmission in Sandwich Structure

Wang Qing1Hong Ming2Zhou Yan－qiu3
1 China Ship Development and Design Center， Wuhan 430064，China
2 School of Naval Architecture Engineering， Faculty of Vehicle Engineering and Mechanics，State Key Laboratory of Structural Analysis for Industrial Equipment， Dalian University of Technology， Dalian 116024，China
3 China Ship Development and Design Center，Shanghai Division，Shanghai 201108，China

With variable assumed sandwich model， the vibration and acoustic characteristics are distinct.Based on the displacement field assumption of the Reissner and Zig－Zag model， the paper established two kinds of composite sandwich plate Finite Element models， and its inherent frequency and Eigen－mode shapes were required by subspace iteration method with the help of the derivation model program.The paper also calculated the acoustic transmission loss by transmission theory of diffused sound field，and studied the numerical simulation of vibration and sound transmission character，the application of sandwich plate＇s face and core with variable thickness ratio and elastic modulus ratio were considered as well.

Reissner model； Zig－Zag model； Mindlin method； vibration character； sound transmission loss

TB532

A

1673－3185（2011）03－68－05

10.3969/j.issn.1673-3185.2011.03.015

2010-03-18

國家自然科學(xué)基金資助項(xiàng)目（51079027）

王 晴（1984－），女，碩士。研究方向：復(fù)合材料結(jié)構(gòu)動(dòng)力學(xué)。E-mail：sunnywangqing＠126.com

洪 明（1959－ ） ，男，教授。 研究方向：結(jié)構(gòu)失效、振動(dòng)和噪聲機(jī)理、預(yù)極及控制研究。E-mail：mhong＠dlut.edu.cn