基于FLUENT嵌套網(wǎng)格的串列雙圓柱渦激振動數(shù)值模擬

2017-01-13 by:CAE仿真在線來源:互聯(lián)網(wǎng)

摘要:管束類結(jié)構(gòu)的流固耦合振動在海洋工程、航天航空、核工程等領(lǐng)域中非常普遍。目前國內(nèi)外學(xué)者針對渦激振動的研究主要采用CFD與CSD聯(lián)合仿真的形式進行,在計算流體力學(xué)中不可避免的需要采用動網(wǎng)格模型進行求解計算,傳統(tǒng)的動網(wǎng)格內(nèi)省包括光順、層鋪與局部重構(gòu)。本文采用fluent17.0版本的嵌套網(wǎng)格方法同時結(jié)合udf進行雙圓柱渦激振動分析,通過分析表明,采用嵌套網(wǎng)格進行渦激振動分析方法可行。

1引言

管束類結(jié)構(gòu)的流固耦合振動在海洋工程、航天航空、核工程等領(lǐng)域中非常普遍。運行經(jīng)驗和科學(xué)研究結(jié)果表明,流致振動是導(dǎo)致管道破裂失效的關(guān)鍵因素,由其導(dǎo)致的管道將會影響設(shè)備運行的可靠性。因此,對流體誘發(fā)的彈性管振動進行深入研究十分必要。本文采用雙向耦合的方法結(jié)合fluent嵌套網(wǎng)格技術(shù),同時考慮流體-結(jié)構(gòu)兩個物理場之間的相互作用,聯(lián)合計算流體力學(xué)和計算結(jié)構(gòu)力學(xué)建立管束的流固耦合振動模型,研究雙彈性管的流固耦合振動。

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

2.1 CFD模型

對于二維不可壓縮的粘性牛頓流體的控制方程為N-S方程,在直角坐標(biāo)系下,連續(xù)性方程為:

動量方程為:

本文利用嵌套網(wǎng)格進行計算分析,具體網(wǎng)格如圖1-1-1所示,計算網(wǎng)格包含一套背景網(wǎng)格及兩套圓柱區(qū)域網(wǎng)格。本文計算中背景網(wǎng)格區(qū)域大小為30D*60D;圓柱直徑為D=0.0023878m,串列兩圓柱中心距為2D。計算過程中,背景網(wǎng)格左側(cè)為速度入口,右側(cè)為壓力出口,上下表面為對稱面;component網(wǎng)格與backgroud網(wǎng)格通過overset邊界條件進行連接。

2.2 流固耦合分析模型

圓柱的結(jié)構(gòu)參數(shù)取為,m=0.003575kg,圓柱直徑D=0.0023878m;密度為67.44kg/m^3,自振圓頻率為44.08,結(jié)構(gòu)剛度為69.47。

流固耦合作用時,其振動方程可寫為:

其中Ft為流體對結(jié)構(gòu)的力,采用四階龍格庫塔法進行離散求解得到每個時間步對應(yīng)的圓柱運動速度和位移,并通過得到的速度和位移來控制圓柱的運動和網(wǎng)格的更新。采用全隱式格式來處理方程可以得到函數(shù):

通過離散得到速度與質(zhì)心位置的關(guān)系式;

上式中,K1、K2、K3、K4的表達式分別為:

考慮到計算是在時域內(nèi)進行的,故給定初始條件如下:

本文考慮了串列兩圓柱兩自由度的渦激振動,對應(yīng)需要利用以上方法求解兩個圓柱x\y方向的速度與位移。

3 結(jié)果分析

圖3-1給出了串列兩圓柱不同時刻尾渦結(jié)構(gòu)圖,由圖可知:兩圓柱中心距為2D時,其尾渦結(jié)構(gòu)與單圓柱擾流尾渦結(jié)構(gòu)相當(dāng)。

圖3-2給出了串列兩圓柱運動軌跡圖,由圖可知:在雷諾數(shù)為400,中心距為2D的情況下,串列兩圓柱的運動軌跡存在一定的差別;在來流速度為0.166m/s的流速下,兩圓柱在流向與橫向發(fā)生大幅振動,圓柱1運動軌跡為“8”字形圖,如圖3-1a所示,同時阻力使圓柱往下游方向移動;上游圓柱的振幅大于下游圓柱的振幅。

圖3-3給出了串列兩圓柱的生阻力系數(shù)時程曲線,由圖可知:t=2s后,兩圓柱的組升力系數(shù)均呈現(xiàn)周期性變化,上游圓柱阻力系數(shù)大于下游圓柱阻力系數(shù);下游圓柱升力系數(shù)振動幅值大于上游圓柱升力系數(shù)振動幅值。

4 結(jié)論

通過上文利用fluent嵌套網(wǎng)格結(jié)合動網(wǎng)格模型對二維串列圓柱渦激振動的分析可知:

通過上文利用fluent嵌套網(wǎng)格結(jié)合動網(wǎng)格模型對二維串列圓柱渦激振動的分析可知:采用嵌套網(wǎng)格進行流固耦合的模擬該方法可行,考慮到嵌套網(wǎng)格自身的特點,計算過程中不涉及網(wǎng)格的重構(gòu),不涉及網(wǎng)格變形,相比其他動網(wǎng)格方法其計算速度更快,模型網(wǎng)格量越大效果越明顯;

附,udf文件