多鈍體流致振動(dòng)是一個(gè)較為復(fù)雜的流固耦合過程,普遍存在于自然界和工程領(lǐng)域。為了減小高幅振動(dòng)時(shí)網(wǎng)格變形引起的計(jì)算誤差,基于非定常Navier-Stokes方程對(duì)二維雙圓柱和三圓柱、三維雙圓柱流致振動(dòng)進(jìn)行數(shù)值求解,采用耦合界面結(jié)合拓?fù)渚W(wǎng)格變形技術(shù),實(shí)現(xiàn)流體與多個(gè)運(yùn)動(dòng)鈍體之間的耦合計(jì)算。將數(shù)值結(jié)果與實(shí)驗(yàn)進(jìn)行比較分析,驗(yàn)證了該數(shù)值方法是處理高振幅多鈍體流致振動(dòng)的有效方法。研究結(jié)果表明上游圓柱的存在對(duì)下游圓柱流致振動(dòng)和旋渦形成產(chǎn)生明顯影響。串列雙圓柱流致振動(dòng)振幅和頻率響應(yīng)與實(shí)驗(yàn)測(cè)試趨勢(shì)一致,清晰觀察到了渦致振動(dòng)初始分支和上部分支;并且當(dāng)Re>8×10⁴時(shí),圓柱流致振動(dòng)由渦致振動(dòng)向馳振過渡。圓柱尾渦形態(tài)隨流致振動(dòng)分支切換發(fā)生變化,當(dāng)馳振發(fā)生時(shí),下游圓柱的尾渦形態(tài)受上游圓柱影響難以捕捉。隨著雙圓柱間距增大,低Re時(shí)下游圓柱受到上游圓柱的抑制作用減弱。三維多柱體流致振動(dòng)計(jì)算結(jié)果更接近實(shí)驗(yàn)值,如何提高三維數(shù)值計(jì)算速度將是下一步研究工作的重點(diǎn)。 

【文章來(lái)源】：振動(dòng)與沖擊. 2019,38(22)北大核心EICSCD

【文章頁(yè)數(shù)】：8 頁(yè)

【部分圖文】：

物理模型Fig．1Physicalmodel

，不要求相鄰界面的網(wǎng)格節(jié)點(diǎn)一一匹配。在具有旋轉(zhuǎn)部件的數(shù)值計(jì)算中要求所有界面網(wǎng)格相互匹配往往非常困難，CI界面是處理這類問題的有效途徑。對(duì)于界面網(wǎng)格不連續(xù)的區(qū)域，分別進(jìn)行網(wǎng)格劃分，然后使用一個(gè)或者多個(gè)CI界面進(jìn)行連接。與滑移網(wǎng)格不同，CI界面包含主界面和副界面，使用加權(quán)插值進(jìn)行流率計(jì)算和傳遞。為避免柱體運(yùn)動(dòng)造成網(wǎng)格拉伸或擠壓變形，保證在數(shù)值計(jì)算過程中維持較高的網(wǎng)格質(zhì)量，通過網(wǎng)格拓?fù)渥冃渭夹g(shù)實(shí)現(xiàn)運(yùn)動(dòng)子塊上下網(wǎng)格的自動(dòng)生成或塌縮，避免網(wǎng)格單元體積的劇烈變化。如圖2所示，振動(dòng)柱體所在子塊的上下側(cè)為動(dòng)態(tài)層，h為動(dòng)態(tài)層與鄰近網(wǎng)格(n)的距離。當(dāng)流致振動(dòng)發(fā)生時(shí)，子塊內(nèi)網(wǎng)格整體運(yùn)動(dòng)，過程如下:(a)子塊向上運(yùn)動(dòng)時(shí)，h變小，當(dāng)h＜hmin時(shí)，網(wǎng)格n崩塌消失;(b)子塊向下運(yùn)動(dòng)時(shí)，h增大，當(dāng)h＞hmax時(shí)，將會(huì)自動(dòng)生成網(wǎng)格n－1。其中hmin和hmax分別是根據(jù)整體網(wǎng)格質(zhì)量而設(shè)定的最小和最大網(wǎng)格厚度。拓?fù)渥冃渭夹g(shù)能自動(dòng)對(duì)網(wǎng)格進(jìn)行更新，有效解決旋轉(zhuǎn)、拉伸和擠壓造成的網(wǎng)格變形問題，有效減小網(wǎng)格變形造成的計(jì)算誤差。圖2拓?fù)渚W(wǎng)格及耦合界面Fig．2Topologicalmeshandcouplinginterface3結(jié)果與討論為了研究彈性支撐多鈍體的流致振動(dòng)響應(yīng)特性，同時(shí)驗(yàn)證本文提出的數(shù)值方法尤其是動(dòng)網(wǎng)格模型的可靠性，本文分別對(duì)二維雙圓柱和三圓柱、三維雙圓柱模型進(jìn)行了數(shù)值計(jì)算。3．1二維雙圓柱流致振動(dòng)二維雙圓柱計(jì)算區(qū)域如圖3所示，頂部和底部距離為9D，與實(shí)驗(yàn)水深相同，速度入口距上游圓柱軸心為25D，下游圓柱軸心距離壓力出口為25D，兩種不同間距的串列雙圓柱模型的軸心間距d分別為2

?hmin和hmax分別是根據(jù)整體網(wǎng)格質(zhì)量而設(shè)定的最小和最大網(wǎng)格厚度。拓?fù)渥冃渭夹g(shù)能自動(dòng)對(duì)網(wǎng)格進(jìn)行更新，有效解決旋轉(zhuǎn)、拉伸和擠壓造成的網(wǎng)格變形問題，有效減小網(wǎng)格變形造成的計(jì)算誤差。圖2拓?fù)渚W(wǎng)格及耦合界面Fig．2Topologicalmeshandcouplinginterface3結(jié)果與討論為了研究彈性支撐多鈍體的流致振動(dòng)響應(yīng)特性，同時(shí)驗(yàn)證本文提出的數(shù)值方法尤其是動(dòng)網(wǎng)格模型的可靠性，本文分別對(duì)二維雙圓柱和三圓柱、三維雙圓柱模型進(jìn)行了數(shù)值計(jì)算。3．1二維雙圓柱流致振動(dòng)二維雙圓柱計(jì)算區(qū)域如圖3所示，頂部和底部距離為9D，與實(shí)驗(yàn)水深相同，速度入口距上游圓柱軸心為25D，下游圓柱軸心距離壓力出口為25D，兩種不同間距的串列雙圓柱模型的軸心間距d分別為2D和2．5D。頂部和底面為壁面邊界。為了強(qiáng)化圓柱的流致振動(dòng)，在圓柱表面沿軸向安裝兩條附屬粗糙帶，粗糙帶覆蓋16°圓柱表面，厚度為0．847mm，數(shù)值計(jì)算中采用壁面函數(shù)進(jìn)行求解。本文計(jì)算采用結(jié)構(gòu)化網(wǎng)格，均經(jīng)過網(wǎng)格無(wú)關(guān)性驗(yàn)證，并將計(jì)算結(jié)果與實(shí)驗(yàn)進(jìn)行比較。圖3計(jì)算區(qū)域及邊界條件Fig．3Calculationregionandboundaryconditions3．1．1雙圓柱軸心距離d=2D當(dāng)串列雙圓柱軸心距離d=2D時(shí)，雙圓柱所在方形網(wǎng)格子區(qū)域之間存在共用耦合界面。當(dāng)圓柱發(fā)生流致振動(dòng)時(shí)，圓柱所在子區(qū)域整體運(yùn)動(dòng)，并在動(dòng)態(tài)層經(jīng)過區(qū)域成功實(shí)現(xiàn)網(wǎng)格自動(dòng)崩塌與生成，如圖4所示。計(jì)算雷諾數(shù)(Ｒe)范圍為3×104≤Ｒe≤10×104，對(duì)應(yīng)折減832振動(dòng)與沖擊2019年第38卷

【參考文獻(xiàn)】：
期刊論文
[1]雙自由度渦激振動(dòng)數(shù)值模擬方法研究[J]. 孫麗萍,張旭,倪問池.  振動(dòng)與沖擊. 2017(23)
[2]質(zhì)量比對(duì)剛性圓柱體渦激振動(dòng)影響的研究[J]. 陳正壽,趙宗文,張國(guó)輝,鄭武,顏盛漢.  振動(dòng)與沖擊. 2017(11)
[3]串列圓柱體尾流、尾渦耦合振動(dòng)試驗(yàn)研究[J]. 關(guān)德寶,黃維平,宋虹,田會(huì)元,耿翱翔,魏立波.  振動(dòng)與沖擊. 2014(22)
[4]串列雙圓柱繞流下游圓柱兩自由度渦致振動(dòng)研究[J]. 陳文曲,任安祿,李廣望.  力學(xué)學(xué)報(bào). 2004(06)



本文編號(hào)：3529421