流固耦合命令求解流固耦合問題

2016-09-20 by:CAE仿真在線來源:互聯網

使用ANSYS計算結構在水中的模態時, FLUID29,FLUID30單元分別用來模擬二維和三維流體部分,相應的結構模型則利用PLANE42單元和SOL ID45等單元來構造,其中,PLANE42和SOL ID45分別是用來構造二維和三維結構模型的單元。FLUID30是流體聲單元,主要用于模擬流體介質及流固耦合問題。該單元有8 個節點,每個節點上有4 個自由度,分別是XYZ上3個方向位移自由度和1個壓力自由度,為各向同性材料。輸入材料屬性時,需要輸入流體的材料密度(作為DENS 輸入)及流體聲速(作為SONC輸入),流體粘性產生的損耗效應忽略不計。FLUID29是FLUID30單元在二維上的簡化,少了一個Z向的位移。SOLID45單元用于構造三維實體結構。單元通過8 個節點來定義,每個節點有3 個沿著XYZ方向平移的自由度。PLANE42是SOLID45單元在二維上的簡化。

在利用ANSYS建模分析時,流場域單元屬性分為2種,由KEYOPT(2)(指定流體和結構分界面處結構是否存在) 控制,在流固耦合交界面上的單元KEYOPT(2) = 0 ,表示分界面處有結構,其他流體單元KEYOPT(2)=1,表示分界面處無結構。流體-結構分界面通過面載荷標志出來,指定FSI label可以把分界面處的結構運動和流體壓力耦合起來,分界面標志在分界面處的流體單元標出。

數值分析的步驟

1) 建立流體單元的實體模型。建立流體模型,需要確定流體域的范圍,可以把無限邊界流體簡化成流體區域的半徑為固體結構半徑的10倍。

2) 標記流固耦合界面。選取流體單元中流固交界面上的節點,執行FSI 命令,流固耦合交界面的處理:流體與固體是兩個獨立的實體,在劃分單元時在兩者交界面上的單元網格要劃分一致,這樣在交界面上的同一位置一般就有兩個重合的節點,一個節點屬于流體單元,一個節點屬于固體單元,這兩個重合節點在交界面的位移強制保持一致。

3) 建立固體結構實體模型。建立固體結構模型,定義單元屬性,采用映射方式進行網格的劃分。
4) 施加約束條件。由于流體區域的尺寸遠大于固體結構尺寸,故可以不考慮流體液面的重力的影響,將流體邊界處的單元節點上施加壓力(PRES) 為零的約束。因為選擇的算例為懸臂結構,在固體結構底部加全約束。

5) 選擇求解算法,進行求解。定義分析類型為模態分析,設定提取頻率階數和提取模態的方法。因為耦合問題的剛度矩陣,質量矩陣都不對稱,需要采用非對稱矩陣法(UNSYMMETRIC)求解。

6) 查看結果。進入后處理模塊,查看結構模型的頻率及振型。

以半浸沒與水中的橋墩模態問題為背景,并假設:

1. 橋墩為實心等截面的實體,實際橋墩模型應該是空心殼體,截面尺寸也非常復雜,因而需要分塊劃分單元。

2. 不考慮橋墩中的鋼筋,簡化為由單一的混凝土構成橋墩,密度1.265e4 kg/m3,彈性模量3.15e10 N/m2,泊松比0.167。

3. 不考慮水面重力影響,忽略自由表面波,水的密度1e3 kg/m3,音速 1.44e3 m/s。

4. 壩體和水的交界面上,固體節點和液體節點位移耦合在一起,在單元中液體和固體粘在一起。

劃分網格后的模型如下:

計算后的模態結果如下,為了方便對比,左側給出了空氣中的模態。

實體橋墩的固有頻率	橫向振動	縱向振動	扭轉振動
第一階	第二階	第三階	第一階	第一階
空氣中頻率	0.5545	3.357	8.97	8.63	5.27
水中頻率	0.5511	3.212	8.56	6.32	5.27
比例系數	0.994	0.957	0.954	0.732	1