基于浸入邊界-譜元法的流體柔性體耦合運(yùn)動(dòng)研究
發(fā)布時(shí)間:2021-03-20 12:54
為了模擬大變形柔性體與黏性流體耦合運(yùn)動(dòng)的動(dòng)力學(xué)行為,該文發(fā)展了一套基于浸入邊界-譜元法的數(shù)值模擬方法。該數(shù)值方法中,流體運(yùn)動(dòng)在歐拉坐標(biāo)下采用譜元法模擬,柔性體運(yùn)動(dòng)在拉格朗日坐標(biāo)系采用有限差分法求解,兩者的耦合通過浸入邊界法實(shí)現(xiàn)。首先對(duì)低雷諾數(shù)方柱繞流后柔性懸臂梁流固耦合振動(dòng)問題進(jìn)行數(shù)值模擬,并將計(jì)算結(jié)果與其它文獻(xiàn)計(jì)算結(jié)果進(jìn)行比較,驗(yàn)證了該耦合計(jì)算方法的可靠性。然后研究了高雷諾數(shù)下方柱繞流后柔性懸臂梁流固耦合振動(dòng)情況,得到柔性結(jié)構(gòu)的耦合振動(dòng)特性和流場(chǎng)動(dòng)態(tài)分布特征。該文的數(shù)值方法放寬了由大拉伸系數(shù)引起的穩(wěn)定性限制,能更靈活地被運(yùn)用于復(fù)雜柔性邊界條件的情況,為后續(xù)對(duì)各種柔性體流體耦合的研究提供了基礎(chǔ)。
【文章來源】:水動(dòng)力學(xué)研究與進(jìn)展(A輯). 2020,35(03)北大核心CSCD
【文章頁數(shù)】:8 頁
【部分圖文】:
算例示意圖
圖4展示了柔性懸臂自由端的位移和速度信息。從模擬結(jié)果得出,由于不可壓縮黏性流經(jīng)過障礙物的不穩(wěn)定分離,旋渦從方柱突出邊緣以一定頻率脫落。方柱尾流中的渦形成非定常的阻力和升力,激勵(lì)柔性懸臂在流場(chǎng)中自由振動(dòng)。根據(jù)Liew等的研究,位移歷程被劃分為兩個(gè)階段:鎖定和跳動(dòng)。結(jié)構(gòu)位移在第一階段逐漸放大,懸臂的振蕩頻率接近其第一個(gè)本征頻率,激發(fā)鎖定(或共振)。在第二階段,結(jié)構(gòu)響應(yīng)收斂于穩(wěn)定的長(zhǎng)期振蕩運(yùn)動(dòng),形成跳動(dòng)(或調(diào)制)。在整個(gè)VIV過程中,渦旋脫落引起鎖定,從而驅(qū)動(dòng)懸臂急劇振蕩[24]。從圖4可以清晰看出這一演變過程。穩(wěn)定后的振動(dòng)頻率經(jīng)過計(jì)算可得為f=3.06,自由端振幅為d=1.22。表3展示了已有文獻(xiàn)對(duì)這個(gè)基本算例的模擬結(jié)果。對(duì)比發(fā)現(xiàn),本文結(jié)果與已有文獻(xiàn)的模擬結(jié)果非常吻合。圖5截取了結(jié)構(gòu)震蕩穩(wěn)定后的三幀典型渦量場(chǎng)?梢钥闯,在振蕩期間,流動(dòng)模式和結(jié)構(gòu)振蕩的不穩(wěn)定特征明顯地展示在不同階段中。在圖中,渦作用在懸臂梁的表面上,并同時(shí)從懸臂兩側(cè)交替脫落,由此激發(fā)懸臂梁大幅度振動(dòng)。值得注意的是,懸臂梁開始移動(dòng)的一側(cè)的渦被削弱并最終被附近的流動(dòng)抑制。因此,在這一側(cè)受到附近的高壓縮。然而,在相反的一側(cè),渦被加強(qiáng)并且在懸臂梁的下游平行脫落。這些渦是由振動(dòng)懸臂梁的尖端產(chǎn)生的速度梯度產(chǎn)的。由于立即暴露于流動(dòng),方柱正面上的壓力分布是正的,但是沿著方柱的側(cè)面檢測(cè)到低壓。沿著柔性懸臂也有低壓區(qū)域促使懸臂梁振動(dòng)。圖3 方柱周邊網(wǎng)格
方柱周邊網(wǎng)格
【參考文獻(xiàn)】:
期刊論文
[1]上游圓柱固定條件下串列三圓柱渦激振動(dòng)響應(yīng)和尾流特性[J]. 張志猛,及春寧,許棟,陳威霖,楊梟梟. 水動(dòng)力學(xué)研究與進(jìn)展(A輯). 2019(02)
[2]不同控制角下附加圓柱對(duì)圓柱渦激振動(dòng)影響[J]. 陳威霖,及春寧,許棟. 力學(xué)學(xué)報(bào). 2019(02)
[3]錯(cuò)列角度對(duì)雙圓柱渦激振動(dòng)影響的數(shù)值模擬研究[J]. 段松長(zhǎng),趙西增,葉洲騰,王凱鵬. 力學(xué)學(xué)報(bào). 2018(02)
[4]五葉側(cè)斜螺旋槳流固耦合特性的非定常數(shù)值模擬[J]. 曹峰,李巍,尤云祥. 水動(dòng)力學(xué)研究與進(jìn)展A輯. 2013(05)
博士論文
[1]譜元方法求解不可壓縮流體流動(dòng)及流動(dòng)線性穩(wěn)定性分析[D]. 梅歡.重慶大學(xué) 2012
本文編號(hào):3091033
【文章來源】:水動(dòng)力學(xué)研究與進(jìn)展(A輯). 2020,35(03)北大核心CSCD
【文章頁數(shù)】:8 頁
【部分圖文】:
算例示意圖
圖4展示了柔性懸臂自由端的位移和速度信息。從模擬結(jié)果得出,由于不可壓縮黏性流經(jīng)過障礙物的不穩(wěn)定分離,旋渦從方柱突出邊緣以一定頻率脫落。方柱尾流中的渦形成非定常的阻力和升力,激勵(lì)柔性懸臂在流場(chǎng)中自由振動(dòng)。根據(jù)Liew等的研究,位移歷程被劃分為兩個(gè)階段:鎖定和跳動(dòng)。結(jié)構(gòu)位移在第一階段逐漸放大,懸臂的振蕩頻率接近其第一個(gè)本征頻率,激發(fā)鎖定(或共振)。在第二階段,結(jié)構(gòu)響應(yīng)收斂于穩(wěn)定的長(zhǎng)期振蕩運(yùn)動(dòng),形成跳動(dòng)(或調(diào)制)。在整個(gè)VIV過程中,渦旋脫落引起鎖定,從而驅(qū)動(dòng)懸臂急劇振蕩[24]。從圖4可以清晰看出這一演變過程。穩(wěn)定后的振動(dòng)頻率經(jīng)過計(jì)算可得為f=3.06,自由端振幅為d=1.22。表3展示了已有文獻(xiàn)對(duì)這個(gè)基本算例的模擬結(jié)果。對(duì)比發(fā)現(xiàn),本文結(jié)果與已有文獻(xiàn)的模擬結(jié)果非常吻合。圖5截取了結(jié)構(gòu)震蕩穩(wěn)定后的三幀典型渦量場(chǎng)?梢钥闯,在振蕩期間,流動(dòng)模式和結(jié)構(gòu)振蕩的不穩(wěn)定特征明顯地展示在不同階段中。在圖中,渦作用在懸臂梁的表面上,并同時(shí)從懸臂兩側(cè)交替脫落,由此激發(fā)懸臂梁大幅度振動(dòng)。值得注意的是,懸臂梁開始移動(dòng)的一側(cè)的渦被削弱并最終被附近的流動(dòng)抑制。因此,在這一側(cè)受到附近的高壓縮。然而,在相反的一側(cè),渦被加強(qiáng)并且在懸臂梁的下游平行脫落。這些渦是由振動(dòng)懸臂梁的尖端產(chǎn)生的速度梯度產(chǎn)的。由于立即暴露于流動(dòng),方柱正面上的壓力分布是正的,但是沿著方柱的側(cè)面檢測(cè)到低壓。沿著柔性懸臂也有低壓區(qū)域促使懸臂梁振動(dòng)。圖3 方柱周邊網(wǎng)格
方柱周邊網(wǎng)格
【參考文獻(xiàn)】:
期刊論文
[1]上游圓柱固定條件下串列三圓柱渦激振動(dòng)響應(yīng)和尾流特性[J]. 張志猛,及春寧,許棟,陳威霖,楊梟梟. 水動(dòng)力學(xué)研究與進(jìn)展(A輯). 2019(02)
[2]不同控制角下附加圓柱對(duì)圓柱渦激振動(dòng)影響[J]. 陳威霖,及春寧,許棟. 力學(xué)學(xué)報(bào). 2019(02)
[3]錯(cuò)列角度對(duì)雙圓柱渦激振動(dòng)影響的數(shù)值模擬研究[J]. 段松長(zhǎng),趙西增,葉洲騰,王凱鵬. 力學(xué)學(xué)報(bào). 2018(02)
[4]五葉側(cè)斜螺旋槳流固耦合特性的非定常數(shù)值模擬[J]. 曹峰,李巍,尤云祥. 水動(dòng)力學(xué)研究與進(jìn)展A輯. 2013(05)
博士論文
[1]譜元方法求解不可壓縮流體流動(dòng)及流動(dòng)線性穩(wěn)定性分析[D]. 梅歡.重慶大學(xué) 2012
本文編號(hào):3091033
本文鏈接:http://sikaile.net/kejilunwen/lxlw/3091033.html
最近更新
教材專著