目錄
結構模型的建立
聲腔模型的建立
結構，聲腔模型的模態(tài)分析
流固耦合模態(tài)分析
結論

摘要

本次分享，對一個簡單的模型進行流固耦合的模態(tài)分析，主要給大家演示如何使用Hypermesh與Nastran對流固耦合的結構進行模態(tài)分析，以及了解聲腔對結構模態(tài)的影響。

流固耦合法被廣泛地應用于聲學和噪音控制領域，對空腔結構（比如汽車車室）進行流固耦合模態(tài)分析，可以了解到聲腔對結構模態(tài)的影響，為研究耦合系統(tǒng)的聲學特性提供可靠的理論參考。

結構模型的建立
本次分享采用一個壁厚為0.5mm，長寬高為100 x 50 x 50mm的長方體作為分析模型，模擬車身車室，命名為structure。單元采用片體shell單元，單元邊長控制在5mm范圍。本次分析給結構模型附上一般的鋼材材料屬性：彈性模量E=206800Mpa, 泊松比NU=0.3，密度RHO=7.85e-9ton/mm3
結構模型如下圖所示：


聲腔模型的建立
對于復雜的結構模型，例如車身車室的聲腔模型，我們一般會采用Hypermesh中的一個專門劃分聲腔網(wǎng)格的小功能 "Acoustic Cavity Mesh" 進行建立，詳細教程下次再作介紹 。但是對于今天我們這個長方體小模型，我們將按以下步驟建立聲腔模型：
首先先創(chuàng)建一個新component的層，命名為cavity；
按照下圖的指示，選擇 3D后點擊 teramesh 進行實體網(wǎng)格的劃分；


進入到體網(wǎng)格劃分界面，點擊Tetranesh parameters，在紅圈內選擇User Controlled，然后把下面的Growth rate改成1.000，這個選項意義在于模型會以1.0的比例進行體網(wǎng)格的劃分，即劃分出來的體網(wǎng)格，從外到內每一層單元的大小將一致。下圖給出兩組不同Growth rate的體網(wǎng)格模型供讀者比較。


1.0

1.5

返回，然后點選Tetra mesh，點擊comps選擇建立好的structure模型，選上split quads into trias表示建立三菱錐實體單元，其余設置如圖所示，點擊mesh進行體網(wǎng)格劃分，聲腔實體網(wǎng)格模型就完成了；


聲腔網(wǎng)格模型出來后，還需要對其進行進一步設置。點擊圖示圖標，把聲腔模型單獨顯示出來全選上所有的節(jié)點，然后點擊edit，進去設置界面把CD-1勾選，界面將出現(xiàn)-1的字樣，點擊return完成設置。這一步設置把聲腔模型的節(jié)點設置為流固耦合的節(jié)點，為后續(xù)軟件計算查找流固耦合的節(jié)點做標志；



聲腔模型的材料選擇流體材料類型MAT10，設置材料密度1.2e-12ton/mm^3，聲速340000mm/s.


聲腔模型單元屬性選擇PSOLID，附上材料cavity，勾選全積分FULL，F(xiàn)CTN選擇PFLUID.


結構，聲腔模型的模態(tài)分析
為了了解聲腔模態(tài)與結構模態(tài)的相互影響，在這里我們先來分別單獨計算出聲腔和結構的自由模態(tài)。聲腔的第一階自由模態(tài)為剛體模態(tài)，聲腔內各點的聲壓幅值相同；結構自由模態(tài)前6階為6個自由度的剛體模態(tài)。要注意，在這里我們使用nastran求解器中EIGRL的方式進行模態(tài)提取，把聲腔與結構模型分開單獨導出提交計算即可。詳細的求解過程與 "白車身自由模態(tài)分析” 相同，點擊超鏈接回顧。下面直接給出聲腔與結構各自的模態(tài)計算結果，由于分析的結構幾何形狀具有對稱性，計算結果也同理具有對稱性，所以你會發(fā)現(xiàn)有些計算出來的模態(tài)頻率與振形結果相等或高度相似。



結構第7~11階自由模態(tài)位移云圖


聲腔第2~5階自由模態(tài)聲壓云圖

流固耦合模態(tài)分析
接下來，我們來討論一下聲腔與結構相互作用下的模態(tài)分析。把結構與聲腔同時導入到hypermesh中進行流固耦合模態(tài)分析的相關設置。

在做結構與空腔的流固耦合分析時，導入的聲腔模型與結構模型必須保證結構單元的結點編號與聲腔單元的結點編號獨立分開，即兩類單元沒有共用的結點編號，否則軟件識別編號出錯無法計算；
由于結構和空氣在界面的相互作用，導致質量和剛度陣中引入了對角線外的耦合項，對耦合系統(tǒng)進行模態(tài)分析時將出現(xiàn)復特征值，所以此時我們采用Nastran的復特征值(Dir. Complex Eigenvalues)的直接法求解器SOL107求解，而不是像以前單獨計算模態(tài)那樣，采用SOL103求解。通過control card > SOL 進行更改；
提取模態(tài)的方式也不再使用EIGRL，而改為EIGC；


同樣，loadstep設置中Type也更改為Complex Eigen(direct)，CMETHOD選上EIGC；

流固耦合分析要記得在control card里面選擇ACMODL，設置耦合界面參數(shù)；提交到nastran計算。

計算完后，我們使用hyperview14.0打開結果文件.op2，此時計算出來的自由模態(tài)結果中，前7階模態(tài)頻率為0（結構6個剛體模態(tài)加聲腔1個等聲壓模態(tài)），下圖給出流固耦合的第8，9階模態(tài)計算結果，動圖的上方分別是第8階流固耦合的結構模態(tài)（左邊）與聲腔模態(tài)（右邊），下方為第9階流固耦合的結構模態(tài)（左邊）與聲腔模態(tài)（右邊）結果。


結構模態(tài)（左）聲腔模態(tài)（右）
第8階模態(tài)（上） 第9階模態(tài)（右）

結構耦合前后的模態(tài)頻率結果對比如下表所示：



總結
流固耦合后，結構出現(xiàn)局部模態(tài)的地方會導致聲腔在該位置也出現(xiàn)局部模態(tài)；
對比一下耦合前后結構的模態(tài)頻率結果，我們發(fā)現(xiàn)，考慮聲腔耦合作用后，結構的模態(tài)頻率都有明顯的提高；
耦合系統(tǒng)并不是兩個系統(tǒng)簡單的疊加，二者之間的相互作用會一定程度上改變原系統(tǒng)模態(tài)的頻率和振型。