• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    電磁力減振過(guò)渡過(guò)程的機(jī)理研究

    2016-10-14 08:36:13張輝范寶春劉夢(mèng)珂
    兵工學(xué)報(bào) 2016年5期

    張輝,范寶春,劉夢(mèng)珂

    (南京理工大學(xué)瞬態(tài)物理國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京210094)

    電磁力減振過(guò)渡過(guò)程的機(jī)理研究

    張輝,范寶春,劉夢(mèng)珂

    (南京理工大學(xué)瞬態(tài)物理國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京210094)

    電介質(zhì)溶液中,電磁場(chǎng)產(chǎn)生的電磁力可以控制流體的運(yùn)動(dòng)。將其用于鈍體繞流時(shí),可以抑制分離和消除渦街,從而達(dá)到減振的目的。為了研究電磁力減振的過(guò)渡過(guò)程,將指數(shù)極坐標(biāo)系建立在運(yùn)動(dòng)的圓柱上,推導(dǎo)了運(yùn)動(dòng)坐標(biāo)中,基于渦生振蕩并考慮電磁力的渦量-流函數(shù)守恒方程及其初始和邊界條件,圓柱表面的水動(dòng)力表達(dá)式以及圓柱振蕩方程。對(duì)圓柱從靜止開(kāi)始振蕩,到發(fā)展為穩(wěn)定振蕩狀態(tài),然后又在電磁力作用下衰減,直至成為新的穩(wěn)定態(tài)發(fā)展過(guò)程進(jìn)行了計(jì)算和討論,描述了圓柱振蕩、脫體渦街和能量傳遞的發(fā)展過(guò)程,以及升阻力相圖的連續(xù)變形和漂移。結(jié)果表明:圓柱振動(dòng)的誘因在于流動(dòng)分離產(chǎn)生的尾渦誘導(dǎo)的升力周期變化;圓柱表面附近分布的切向電磁力可以有效改變邊界層附近流體的流動(dòng),抑制流動(dòng)分離,從而使流場(chǎng)趨于對(duì)稱,消除升力的周期振動(dòng),最終抑制圓柱的振動(dòng)甚至在電磁力較大時(shí)可以消除振動(dòng)。

    兵器科學(xué)與技術(shù);流體控制;渦生振蕩;電磁力;減振

    0 引言

    黏性流體在運(yùn)動(dòng)物體表面形成的邊界層可使該物體減速,產(chǎn)生振動(dòng)甚至失穩(wěn)[1],這往往不是人們所期望的,而這些現(xiàn)象通常可以通過(guò)流體邊界層的控制加以抑制。動(dòng)力推進(jìn)的物體,通過(guò)對(duì)其表面邊界層的有效控制,可以提高推進(jìn)效率、增加運(yùn)動(dòng)的穩(wěn)定性、減少意外的變形和損壞,這些對(duì)改進(jìn)飛機(jī)、艦船、導(dǎo)彈或普通彈箭的性能是非常有價(jià)值的。

    長(zhǎng)期以來(lái),人們發(fā)現(xiàn)許多控制邊界層流動(dòng)的方法,其中有些不需向流場(chǎng)提供能量,稱作被動(dòng)控制,如加置肋條、帶狹縫的板和二次圓柱等。有些則需要向流場(chǎng)添加能量,稱作主動(dòng)控制,如振蕩和旋轉(zhuǎn)圓柱,聲波干擾、表面狹縫吹吸和熱效應(yīng)等。對(duì)于主動(dòng)控制,向流場(chǎng)提供能量可以根據(jù)流場(chǎng)信息來(lái)調(diào)整,即可以進(jìn)行反饋式控制,這種控制方式更具實(shí)用價(jià)值。利用電磁場(chǎng)控制邊界層流動(dòng)是一種主動(dòng)控制方法,它可以靈活改變電磁力的方向,實(shí)現(xiàn)反饋式控制,甚至制成微機(jī)電系統(tǒng)(MFMS),因此受到廣泛關(guān)注[2]。早在20世紀(jì)中葉,人們就設(shè)想用電磁力控制電解質(zhì)溶液的流動(dòng)。Gai1itis等[3]設(shè)計(jì)了一種由電極和磁極交錯(cuò)布置的電磁場(chǎng)激活板,將其浸入流動(dòng)的弱電解質(zhì)時(shí),激活板附近形成的洛倫茲力可以改變流體邊界層結(jié)構(gòu)。Weier等[4]將此類激活板包覆在圓柱兩側(cè),對(duì)由此形成的圓柱繞流進(jìn)行了實(shí)驗(yàn)研究,實(shí)現(xiàn)了電磁力對(duì)圓柱繞流流場(chǎng)的控制。Crawford等[5]則從理論上討論了激活板的電磁場(chǎng)和洛倫茲力分布。Kim等[6]、Posdziech等[7]發(fā)現(xiàn)常電磁力和脈沖電磁力對(duì)抑制升力、穩(wěn)定流場(chǎng)都有一定的效果。國(guó)內(nèi)上海交通大學(xué)尤云祥教授課題組[8]對(duì)固定圓柱的電磁激活板寬度對(duì)流場(chǎng)控制效果的影響進(jìn)行了數(shù)值研究。另外,Zhang等曾對(duì)圓柱的渦生振蕩進(jìn)行了數(shù)值研究[9],對(duì)電磁力減振過(guò)程已做過(guò)相關(guān)的實(shí)驗(yàn)研究[10-11],而對(duì)于深層次的電磁控制機(jī)理和動(dòng)態(tài)過(guò)程仍需相關(guān)的數(shù)值研究去揭示。

    本文將坐標(biāo)系建立在運(yùn)動(dòng)的圓柱上,推導(dǎo)了運(yùn)動(dòng)坐標(biāo)中,考慮電磁力渦生振蕩的渦量-流函數(shù)守恒方程及其初始和邊界條件,圓柱表面的水動(dòng)力表達(dá)式以及圓柱振蕩方程。對(duì)圓柱從靜止至開(kāi)始振蕩,到發(fā)展為穩(wěn)定振蕩狀態(tài),然后又在電磁力作用下衰減,直至成為新的穩(wěn)定態(tài)整個(gè)發(fā)展過(guò)程進(jìn)行了計(jì)算和討論,描述了圓柱振蕩、脫體渦街和能量傳遞的發(fā)展過(guò)程,以及升阻力相圖的連續(xù)變形和漂移。

    1 守恒方程

    將圓柱置入均勻流動(dòng)的流體中,當(dāng)雷諾數(shù)大于某值時(shí),由于渦在圓柱表面的周期性脫落,導(dǎo)致其表面的水動(dòng)力周期變化。若圓柱被固定在柔性支架上,在周期變化的力作用下,圓柱將周期性的振動(dòng),稱為渦生振蕩(VIV)。振蕩的圓柱又會(huì)影響圓柱周?chē)黧w的流動(dòng),進(jìn)而改變圓柱表面的水動(dòng)力,從而改變圓柱的運(yùn)動(dòng),如此反復(fù)。因此,這是一個(gè)標(biāo)準(zhǔn)的流體-固體耦合問(wèn)題。

    為了控制渦生振蕩,對(duì)置入弱電解質(zhì)溶液中的圓柱表面包覆由電極條和磁極條相間排列的電磁激活板,如圖1所示。圓柱表面附近將形成電磁力場(chǎng),由右手定則可知,電磁力的方向與流體流動(dòng)方向一致。流體在電磁力的作用下加速,邊界層的結(jié)構(gòu)因此而改變。

    圖1 電磁控制渦生振蕩原理圖Fig.1 Princip1e scheme of VIV contro11ed by Lorentz force

    對(duì)于二維情況,電磁力沿軸向平均后,洛倫茲力F的無(wú)量綱形式[4,7]可以表示為

    將坐標(biāo)系建立在振動(dòng)圓柱上。對(duì)于不可壓縮的二維流動(dòng),在指數(shù)極坐標(biāo)(ξ,η)下(r=e2πξ,θ= 2πη),考慮電磁力的無(wú)量綱形式,渦量流函數(shù)方程為

    若圓柱僅沿橫向(垂直于流向)振動(dòng),則初始條件和邊界條件[9]為

    2 圓柱表面水動(dòng)力

    2.1剪應(yīng)力與壓力

    圓柱受到流體的力Fθ*,由剪應(yīng)力和壓力兩部分組成,即

    式中

    顯然,電磁力通過(guò)改變流場(chǎng)和圓柱的運(yùn)動(dòng)狀態(tài)來(lái)改變圓柱表面的剪應(yīng)力,但并未直接出現(xiàn)在剪應(yīng)力方程中。

    式中:pθ為圓柱表面角度為θ處的壓力;p∞為無(wú)窮遠(yuǎn)處的壓力。

    因?yàn)殡姶帕H存在切向分量,故運(yùn)動(dòng)坐標(biāo)系下的動(dòng)量方程為

    在圓柱表面

    沿著η方向從η=0到η對(duì)(12)式積分,可得

    沿著ξ方向(η=0)從ξ=0到∞對(duì)(9)式積分,可得

    因此

    2.2阻力和升力

    阻力分布函數(shù)

    式中:下腳標(biāo)dF表示渦生阻力,該值受到場(chǎng)電磁力的影響;dL表示壁電磁力誘導(dǎo)的阻力(電磁推力),此值僅與壁面電磁力有關(guān),而與流動(dòng)無(wú)關(guān)。

    式中:下腳標(biāo)1F表示渦生的升力;1L表示壁面電磁力誘導(dǎo)的升力,此值與流動(dòng)無(wú)關(guān)。

    總阻力Cd為

    總升力C為

    顯然,(23)式中作用于圓柱的升力由4部分組成,其中,方程右側(cè)第1項(xiàng)C1F為渦生力,與圓柱表面的渦量和渦通量有關(guān),該值受場(chǎng)電磁力的影響;第2項(xiàng)C1L為電磁升力,此值僅與壁面電磁力有關(guān),而與流動(dòng)無(wú)關(guān),電磁力對(duì)稱分布時(shí)C1L=0;第3項(xiàng)為慣性力,與圓柱的加速度有關(guān);第4項(xiàng)為粘性阻尼力Cdamp,與雷諾數(shù)和圓柱的運(yùn)動(dòng)速度有關(guān)。第2項(xiàng)、第3項(xiàng)和第4項(xiàng)均與流場(chǎng)的變化無(wú)關(guān)。

    2.3圓柱運(yùn)動(dòng)方程

    無(wú)量綱的圓柱運(yùn)動(dòng)方程

    每個(gè)周期內(nèi)能量傳遞的無(wú)量綱形式為

    數(shù)值計(jì)算時(shí),動(dòng)量方程(4)式采用交替方向隱式格式,流函數(shù)方程(5)式采用快速傅里葉變換格式,圓柱運(yùn)動(dòng)方程(25)式采用Runge-Kutta法[9]。計(jì)算空間步長(zhǎng)Δξ=0.004,Δη=0.002,時(shí)間Δt=0.005.

    3 結(jié)果與討論

    渦的周期脫落產(chǎn)生振蕩升力,從而使靜止圓柱在升力的作用下逐漸發(fā)展為穩(wěn)定振蕩的圓柱,在電磁力加載后,邊界層流體的動(dòng)量增大,抑制圓柱兩側(cè)流體的分離和脫體,減小脫體導(dǎo)致的升力振蕩,流場(chǎng)和圓柱又經(jīng)歷一段非定常的發(fā)展過(guò)程,最終會(huì)達(dá)到新的穩(wěn)定狀態(tài)。

    圖2為圓柱從固定至穩(wěn)定振蕩,再經(jīng)電磁力作用,趨于新的穩(wěn)定過(guò)程的位移隨時(shí)間變化曲線。在t1=446時(shí)刻,解除橫向約束,圓柱開(kāi)始振蕩,振幅逐漸增大。在t=620時(shí)刻,達(dá)到穩(wěn)定狀態(tài),振幅不再變化。在t2=650時(shí)刻,加載電磁力。電磁力作用下,圓柱振蕩逐漸衰減(圖2中實(shí)線和虛線分別對(duì)應(yīng)電磁力作用參數(shù)N=0.8和N=3.0)。若電磁力足夠大(N=3.0),圓柱會(huì)趨于靜止,不再振蕩。

    當(dāng)振蕩圓柱遠(yuǎn)離平衡位置時(shí),升力C1F(t)做功為正,即流體對(duì)圓柱做功,增加圓柱的能量,反之亦然。另外,阻尼力Cdamp做功始終為負(fù),即減少圓柱的能量。圖3所示(與圖2對(duì)應(yīng))從靜止圓柱到穩(wěn)定振蕩再到電磁力作用下趨于新的穩(wěn)定振蕩全過(guò)程中,流體-固體間的能量傳遞變化,其中圖3(a)、圖3(b)分別對(duì)應(yīng)電磁力作用參數(shù)N=0.8和N= 3.0.虛線表示升力C1F(t)做功E1,點(diǎn)劃線表示阻尼力Cdamp做功Ed,實(shí)線表示每個(gè)周期轉(zhuǎn)移的總能量E,是E1與Ed之和。

    圖2 渦生振蕩發(fā)展和電磁力抑制過(guò)程中圓柱的位移變化Fig.2 Disp1acement variation of cy1inder over time in the process of VIV evo1ution and suppression by Lorentz force

    圖3 渦生振蕩發(fā)展和電磁力抑制過(guò)程中的能量變化Fig.3 Variation of transferred energy for VIV evo1ution and suppression by Lorentz force

    圓柱剛開(kāi)始振蕩時(shí),E1的絕對(duì)值大于Ed,因此升力做的正功占主導(dǎo)。隨著圓柱的振幅逐漸增大,E1的絕對(duì)值先增大后減小,而Ed的絕對(duì)值單調(diào)增大至與E1的絕對(duì)值相等,此時(shí)二者之和,即E為0.因此,圓柱達(dá)到穩(wěn)定振動(dòng)狀態(tài)。t=650時(shí)加載電磁力,E1迅速減小,以致有一段時(shí)間為負(fù)值。能量從圓柱傳遞給流體,圓柱振蕩衰減。最終總能量E=0,此時(shí)E1為正,Ed為負(fù),圓柱以較小振幅穩(wěn)定振蕩,如圖3(a)所示(N=0.8)。當(dāng)電磁力足夠大(N=3.0)時(shí),圓柱最終不再振蕩,此時(shí)E1=Ed=0,如圖3(b)所示。

    圓柱從固定至穩(wěn)定振蕩的發(fā)展過(guò)程中,流場(chǎng)渦量變化如圖4所示,其中紅色表示正渦,藍(lán)色表示負(fù)渦,“+”表示圓柱從固定釋放的初始0位。圖4中的時(shí)刻Di與圖2相對(duì)應(yīng),即圓柱的位置處于上側(cè)最大位移處。在t1=446時(shí)刻,圓柱的橫向約束被解除,在升力作用下開(kāi)始振蕩。由于能量從流體轉(zhuǎn)移到圓柱,因此圓柱的振幅增大,對(duì)應(yīng)的流場(chǎng)如圖4中的D1~D4.當(dāng)總能量達(dá)到平衡時(shí),圓柱的振蕩也達(dá)到穩(wěn)定,此時(shí)流場(chǎng)對(duì)應(yīng)D5,計(jì)算結(jié)果與文獻(xiàn)[10]中實(shí)驗(yàn)結(jié)果一致。

    電磁力控制過(guò)程中,振蕩圓柱的繞流變化如圖5所示。時(shí)刻Di與圖2相同,對(duì)應(yīng)于圓柱處于上側(cè)最大位移。D5時(shí)刻,圓柱的渦生振蕩已經(jīng)穩(wěn)定。此后,在t2=650時(shí)刻加載電磁力,邊界層的流體在電磁力的作用下加速,流動(dòng)分離得到抑制,圓柱上下兩側(cè)分離點(diǎn)的距離減小,尾流渦被拉長(zhǎng),渦距沿流向變大,而沿橫向變小,如圖5(a)~圖5(c)所示。電磁力足夠大時(shí),分離點(diǎn)消失,流場(chǎng)對(duì)稱且定常,如圖5(d)~圖5(f)所示,計(jì)算結(jié)果與文獻(xiàn)[10]中實(shí)驗(yàn)結(jié)果一致。

    流場(chǎng)的變化導(dǎo)致升阻力的變化。渦生振蕩發(fā)展和電磁力控制過(guò)程中,振蕩圓柱的渦生升阻力CdF-C1F相圖的發(fā)展變化如圖 6所示。A1B1C1D1A1對(duì)應(yīng)固定圓柱的升阻力相圖,由于圓柱振蕩對(duì)圓柱上下兩側(cè)剪切層的作用,使相圖逐漸發(fā)生180°的反轉(zhuǎn)。隨著圓柱振蕩的加劇,圓柱的能量增大,點(diǎn)A 與C分離,打破了曲線的鏡像對(duì)稱。另外,阻力平均值以及阻力和升力的振幅增大,導(dǎo)致曲線從左向右不斷延伸,直至振蕩達(dá)到穩(wěn)定,A與C再次重合,即對(duì)應(yīng)相圖A5B5C5D5A5.電磁力作用下,盡管總阻力Cd減小,但渦生阻力CdF是增大的。因此,加載電磁力后,曲線顯著向右移動(dòng)。由于電磁力對(duì)流動(dòng)分離的抑制,并使流場(chǎng)趨于對(duì)稱,故升力C1F減小,曲線逐漸萎縮,從而使圓柱的振蕩減弱,并最終以較小的振幅穩(wěn)定振蕩,對(duì)應(yīng)圖 6(a)中的閉合曲線A8B8C8D8A8,此時(shí)A點(diǎn)與C點(diǎn)也再次重合。N=3.0時(shí),電磁力能夠完全抑制流動(dòng)分離,使繞流流場(chǎng)完全對(duì)稱,此時(shí)相圖曲線逐漸萎縮成升力為0的點(diǎn)(A8B8C8D8A8).

    圖4 振動(dòng)發(fā)展過(guò)程中流場(chǎng)渦量的變化Fig.4 Variation of vorticity in f1ow fie1d during VIV evo1ution

    圖5 電磁力抑制過(guò)程中振蕩圓柱的流場(chǎng)渦量變化Fig.5 Variation of vorticity in f1ow fie1d during suppression by Lorentz force

    電磁力的作用下,渦生振蕩得到抑制,最終以較小振幅穩(wěn)定振蕩,控制后的振幅大小與電磁力的強(qiáng)度有關(guān)。圖7為控制后的振幅隨電磁力強(qiáng)度的變化圖。圓柱振幅隨N的增大而減小,當(dāng)N足夠大時(shí),振蕩被完全抑制,圓柱靜止(圖7中虛線表示電磁力的臨界值,介于2.5到3.0之間)。

    圖6 渦生振蕩發(fā)展和電磁力抑制過(guò)程的升阻力相圖Fig.6 Lift-drag phase diagram for VIV deve1opment and suppression by Lorentz force

    圖7 圓柱的振幅隨電磁力強(qiáng)度(作用參數(shù)N)的變化Fig.7 VIV amp1itude versus interaction parameter N

    4 結(jié)論

    本文將坐標(biāo)系建立在運(yùn)動(dòng)的圓柱上,推導(dǎo)了運(yùn)動(dòng)坐標(biāo)中,考慮電磁力渦生振蕩的渦量-流函數(shù)守恒方程及其初始和邊界條件,圓柱表面的水動(dòng)力表達(dá)式以及圓柱振蕩方程。對(duì)圓柱從靜止開(kāi)始振蕩,到發(fā)展為穩(wěn)定振蕩狀態(tài),然后又在電磁力作用下衰減,直至成為新的穩(wěn)定態(tài)發(fā)展過(guò)程進(jìn)行了計(jì)算和討論。結(jié)果表明,圓柱振動(dòng)的誘因在于尾渦誘導(dǎo)的升力周期變化,而圓柱表面附近分布的切向電磁力可以有效改變邊界層附近流體的流動(dòng),抑制流動(dòng)分離,從而使流場(chǎng)趨于對(duì)稱,消除升力的周期振動(dòng),最終抑制圓柱的振動(dòng)。

    (References)

    [2]Tu J H,Zhou D,Bao Y,et a1.F1ow-induced vibrations of two circu1ar cy1inders in tandem with shear f1ow at 1ow Reyno1ds number[J].Journa1 of F1uids and Structures,2015,59:224-251.

    [2]Reddy P D S,Bandyopadhyay D,Joo S W,et a1.Parametric study on instabi1ities in a two-1ayer e1ectromagnetohydrodynamic channe1 f1ow confined between two para11e1 e1ectrodes[J].Physica1 Review F,2011,83(3):036313.

    [3]Gai1itis A,Lie1ausis O.On a possibi1ity to reduce the hydrodynamica1 resistance of a p1ate in a e1ectro1yte[J].App1ied Magnetohydrodynamics,1961,12:143-146.

    [4]Weier T,Gerbeth G,Posdziedch O,et a1.Fxperiments on cy1inder wake stabi1ization in an e1ectro1yte so1ution by means of e1ectromagnetic forces 1oca1ized on cy1inder surface[J].Fxperimenta1 Therma1 and F1uid Science,1998,16(1):84-91.

    [5]Crawford C,Karniadakis G F.Reyno1ds stress ana1ysis of FMHD-contro11ed wa11 turbu1ence[J].Physics of F1uids,1997,9(3):788-806.

    [6]Kim S,Lee C M.Contro1 of f1ows around a circu1ar cy1inder:suppression of osci11atory 1ift force[J].F1uid Dynamics Research,2001,29(1):47-63.

    [7]Posdziech O,Grundmann R.F1ectromagnetic contro1 of seawater f1ow around circu1ar cy1inders[J].Furopean Journa1 of Mechanics-B/F1uids,2001,20(2):255-274.

    [8]尹紀(jì)富,尤云祥,胡天群,等.電磁極寬度對(duì)圓柱繞流場(chǎng)影響的數(shù)值分析[J].力學(xué)學(xué)報(bào),2013,45(4):493-506. YIN Ji-fu,YOU Yun-xiang,HU Tian-qun,et a1.Numerica1 ana1-ysis for the effect of the e1ectromagnetic actuator width on f1ow structures around a circu1ar cy1inder[J].Chinese Journa1 of Theoretica1 and App1ied Mechanics,2013,45(4):493-506.(in Chinese)

    [9]Zhang H,F(xiàn)an B C,Chen Z H,et a1.An in-depth study on vortex-induced vibration of a circu1ar cy1inder with shear f1ow[J]. Computers and F1uids,2014,100(1):30-44.

    [10]Zhang H,F(xiàn)an B C,Chen Z H,et a1.Fxperimenta1 investigations on the mechanism of 1ift amp1ification and osci11atory suppression of Lorentz force[J].Defence Techno1ogy,2013,9(2):110-116.

    [11]Zhang H,F(xiàn)an B C,Li H Z.Lift enhancement and osci11atory suppression of vortex-induced vibration in shear f1ow by Lorentz force[J].Journa1 of China Ordnance,2012,8(3):139-145.

    Mechanism Investigation on Transient Process of Oscillatory Suppression by Lorentz Force

    ZHANG Hui,F(xiàn)AN Bao-chun,LIU Meng-ke
    (Nation Key Laboratory of Transient Physics,Nanjing University of Science and Techno1ogy,Nanjing 210094,Jiangsu,China)

    The f1ow of weak e1ectro1yte so1ution can be contro11ed by Lorentz force generated by the suitab1y-chosen e1ectromagnetic fie1d.Lorentz force can be used to suppress the vortex shedding and e1iminate the vortex street in the f1ow around a b1uff body for the suppression of vortex-induced vibration (VIV).To investigate the evo1utions of VIV suppression by Lorentz force,an exponentia1 po1ar coordinate system is set up on a moving cy1inder,and the stream function-vorticity equations considering the Lorentz force and its initia1/boundary conditions,an expression of hydrodynamic force on the cy1inder surface as we11 as the cy1inder response equations are then derived.The who1e evo1utions of cy1inder from resting and then undergoing deve1opment and suppression due to Lorentz force are ca1cu1ated and discussed.The deve1opment process of cy1inder vibration,vortexes shedding and energy transfer,and the deformation and shift of drag-1ift phase diagram are described.The resu1ts show that the reason of cy1inder osci11ation is the periodic vibration of 1ift which is induced by wake vortexes due to f1ow separation. The f1uid f1ow near boundary 1ayer is changed and the f1ow separation is suppressed by the Lorentz force on the surface of cy1inder so that the f1ow fie1d tends to symmetry and the periodic vibration of 1ift is sup-pressed.The cy1inder osci11ation can be suppressed and even e1iminated for a 1arge va1ue of Lorentz force.

    ordnance science and techno1ogy;f1ow contro1;vortex-induced vibration;Lorentz force;osci11atory suppression

    O361

    A

    1000-1093(2016)05-0859-07

    10.3969/j.issn.1000-1093.2016.05.012

    2015-07-16

    國(guó)家自然科學(xué)基金項(xiàng)目(11202102);高等學(xué)校全國(guó)優(yōu)秀博士學(xué)位論文作者專項(xiàng)資金項(xiàng)目(201461)

    張輝(1981—),教授,博士生導(dǎo)師。F-mai1:zhanghui1902@126.com

    考比视频在线观看| 男女高潮啪啪啪动态图| 99久国产av精品国产电影| 日韩制服丝袜自拍偷拍| 欧美日本中文国产一区发布| 亚洲精品,欧美精品| 亚洲综合色网址| 精品少妇内射三级| 国产爽快片一区二区三区| 久久亚洲国产成人精品v| 制服丝袜香蕉在线| 亚洲精品乱久久久久久| av电影中文网址| 国产成人aa在线观看| 久久久国产欧美日韩av| 丝袜美腿诱惑在线| 亚洲欧美成人综合另类久久久| 久久影院123| 777久久人妻少妇嫩草av网站| 国产精品久久久av美女十八| 亚洲第一区二区三区不卡| 日韩制服骚丝袜av| 国产av国产精品国产| 在线观看www视频免费| 国产精品国产av在线观看| 寂寞人妻少妇视频99o| 亚洲av男天堂| 9191精品国产免费久久| 亚洲精品一区蜜桃| 亚洲成人av在线免费| 国产欧美日韩综合在线一区二区| 在线观看免费高清a一片| 香蕉精品网在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 天堂8中文在线网| 少妇 在线观看| 久久婷婷青草| 精品人妻偷拍中文字幕| 久久ye,这里只有精品| 黄色视频在线播放观看不卡| 午夜福利一区二区在线看| 男的添女的下面高潮视频| 国产一区亚洲一区在线观看| av电影中文网址| 精品亚洲成国产av| 免费女性裸体啪啪无遮挡网站| 亚洲av.av天堂| 99久久综合免费| 女人被躁到高潮嗷嗷叫费观| 女人久久www免费人成看片| av免费观看日本| 青青草视频在线视频观看| 国产精品一二三区在线看| 午夜免费男女啪啪视频观看| 大香蕉久久成人网| 久久人人爽av亚洲精品天堂| 国产成人欧美| 女性生殖器流出的白浆| 久久久久精品久久久久真实原创| 大片免费播放器 马上看| 亚洲美女搞黄在线观看| 岛国毛片在线播放| 日韩伦理黄色片| av在线播放精品| 亚洲成人一二三区av| 国产精品国产av在线观看| 这个男人来自地球电影免费观看 | www日本在线高清视频| 新久久久久国产一级毛片| 国产熟女欧美一区二区| 一级片'在线观看视频| www日本在线高清视频| 精品一区二区三卡| 老鸭窝网址在线观看| 大片电影免费在线观看免费| 女性被躁到高潮视频| 欧美在线黄色| 亚洲精品自拍成人| 国产综合精华液| 超碰97精品在线观看| 国产成人一区二区在线| 少妇被粗大猛烈的视频| 高清不卡的av网站| 老熟女久久久| 亚洲av电影在线进入| 亚洲国产av新网站| 一级毛片我不卡| 国产野战对白在线观看| 男女免费视频国产| 伦理电影免费视频| 999精品在线视频| 捣出白浆h1v1| 国产一区二区在线观看av| 丝袜美足系列| 在线观看人妻少妇| 午夜福利网站1000一区二区三区| 免费久久久久久久精品成人欧美视频| 一边摸一边做爽爽视频免费| 不卡视频在线观看欧美| 久久这里只有精品19| 久久韩国三级中文字幕| 国产精品麻豆人妻色哟哟久久| 国产有黄有色有爽视频| 国产亚洲欧美精品永久| 久久精品夜色国产| 久久av网站| 欧美精品一区二区免费开放| 久久99热这里只频精品6学生| 丝袜在线中文字幕| 亚洲精品美女久久av网站| 99香蕉大伊视频| 久久韩国三级中文字幕| 亚洲欧美色中文字幕在线| 国产福利在线免费观看视频| 国产高清国产精品国产三级| 黄色 视频免费看| 精品少妇黑人巨大在线播放| 久久久精品国产亚洲av高清涩受| 欧美 亚洲 国产 日韩一| 亚洲第一av免费看| 丝袜人妻中文字幕| 亚洲,一卡二卡三卡| 伊人久久大香线蕉亚洲五| 老司机影院毛片| 最近最新中文字幕免费大全7| 春色校园在线视频观看| 人妻 亚洲 视频| 最新的欧美精品一区二区| 黑人欧美特级aaaaaa片| 亚洲国产欧美日韩在线播放| 大陆偷拍与自拍| 免费不卡的大黄色大毛片视频在线观看| 亚洲av免费高清在线观看| 在线观看免费视频网站a站| 一级片免费观看大全| 欧美少妇被猛烈插入视频| 少妇 在线观看| 在线精品无人区一区二区三| 色婷婷av一区二区三区视频| 97人妻天天添夜夜摸| 国产不卡av网站在线观看| 国产精品一国产av| 久久午夜福利片| 国产一区有黄有色的免费视频| 男女啪啪激烈高潮av片| 91在线精品国自产拍蜜月| 亚洲美女搞黄在线观看| 久久久久久久久久久久大奶| 黄色怎么调成土黄色| 国产无遮挡羞羞视频在线观看| 伦精品一区二区三区| 亚洲av成人精品一二三区| 国产在视频线精品| 久久久精品免费免费高清| 成人手机av| 免费观看无遮挡的男女| 青春草视频在线免费观看| 国产黄色免费在线视频| 嫩草影院入口| 亚洲人成网站在线观看播放| 久久精品国产自在天天线| 在线观看免费高清a一片| 亚洲国产毛片av蜜桃av| 欧美亚洲 丝袜 人妻 在线| 久久97久久精品| 在线观看人妻少妇| 2022亚洲国产成人精品| 色播在线永久视频| 久久久精品免费免费高清| av免费观看日本| 精品久久久精品久久久| 日韩一区二区视频免费看| 午夜av观看不卡| 9191精品国产免费久久| 欧美激情 高清一区二区三区| 中文字幕人妻熟女乱码| www日本在线高清视频| 欧美激情 高清一区二区三区| 亚洲欧洲国产日韩| 99久国产av精品国产电影| 宅男免费午夜| 国产精品一国产av| 亚洲精品日本国产第一区| 亚洲,欧美,日韩| 国产一区二区在线观看av| 下体分泌物呈黄色| 久久久国产精品麻豆| 少妇的丰满在线观看| 在线免费观看不下载黄p国产| 国产一区二区三区综合在线观看| 高清av免费在线| 欧美亚洲日本最大视频资源| 在线天堂最新版资源| 男女国产视频网站| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕| 好男人视频免费观看在线| 欧美在线黄色| 免费观看a级毛片全部| 性色avwww在线观看| 久久久精品免费免费高清| 欧美bdsm另类| 国产精品99久久99久久久不卡 | 少妇被粗大猛烈的视频| av女优亚洲男人天堂| 欧美黄色片欧美黄色片| 在线观看三级黄色| 女人被躁到高潮嗷嗷叫费观| 免费观看av网站的网址| xxxhd国产人妻xxx| 久久精品aⅴ一区二区三区四区 | 丝袜在线中文字幕| 午夜福利网站1000一区二区三区| 日日爽夜夜爽网站| 青春草视频在线免费观看| 香蕉国产在线看| 99九九在线精品视频| 一区二区日韩欧美中文字幕| 久久久久精品久久久久真实原创| 咕卡用的链子| 老女人水多毛片| 老鸭窝网址在线观看| 一级片'在线观看视频| 亚洲欧美精品综合一区二区三区 | 男女国产视频网站| 日日啪夜夜爽| 91国产中文字幕| 高清欧美精品videossex| 高清视频免费观看一区二区| 在线观看三级黄色| 亚洲图色成人| 欧美亚洲 丝袜 人妻 在线| 精品国产国语对白av| 日韩一区二区视频免费看| 欧美日韩一级在线毛片| 午夜激情av网站| av有码第一页| 国产视频首页在线观看| 国产亚洲精品第一综合不卡| 国产精品国产三级专区第一集| 精品少妇内射三级| 成人午夜精彩视频在线观看| 免费在线观看黄色视频的| 老司机影院成人| 99re6热这里在线精品视频| xxx大片免费视频| 久久精品久久久久久久性| 丝袜美腿诱惑在线| 亚洲色图综合在线观看| 国产av一区二区精品久久| 又大又黄又爽视频免费| 99热全是精品| 成年美女黄网站色视频大全免费| 中文字幕av电影在线播放| 日韩电影二区| 精品国产一区二区久久| 精品国产超薄肉色丝袜足j| 高清黄色对白视频在线免费看| 久久精品久久精品一区二区三区| 午夜福利乱码中文字幕| 国产一区亚洲一区在线观看| 亚洲精品一区蜜桃| 国产黄频视频在线观看| 国产精品久久久久久精品电影小说| 不卡视频在线观看欧美| 亚洲精品国产色婷婷电影| 18禁观看日本| 中文欧美无线码| 春色校园在线视频观看| 久久女婷五月综合色啪小说| 午夜免费鲁丝| 1024香蕉在线观看| 色哟哟·www| av线在线观看网站| 青草久久国产| 人妻系列 视频| 久久人妻熟女aⅴ| 久久97久久精品| 久热这里只有精品99| 日韩不卡一区二区三区视频在线| 妹子高潮喷水视频| 欧美亚洲 丝袜 人妻 在线| 高清欧美精品videossex| 91精品三级在线观看| 亚洲国产av影院在线观看| 少妇熟女欧美另类| 纵有疾风起免费观看全集完整版| 午夜91福利影院| 亚洲国产精品999| 国产精品亚洲av一区麻豆 | 日本黄色日本黄色录像| 亚洲精品国产色婷婷电影| 午夜福利视频在线观看免费| 日韩精品有码人妻一区| 晚上一个人看的免费电影| a级毛片在线看网站| 2021少妇久久久久久久久久久| 精品人妻熟女毛片av久久网站| 色哟哟·www| 丝袜喷水一区| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久精品古装| 中文天堂在线官网| 亚洲av免费高清在线观看| 久久女婷五月综合色啪小说| 午夜福利视频精品| 伦理电影大哥的女人| 精品第一国产精品| 免费在线观看完整版高清| 91精品国产国语对白视频| 赤兔流量卡办理| av网站在线播放免费| 精品人妻一区二区三区麻豆| av片东京热男人的天堂| 亚洲伊人色综图| 香蕉国产在线看| 免费黄色在线免费观看| 久久免费观看电影| 欧美日韩精品网址| 伊人亚洲综合成人网| 国产精品久久久久成人av| 性高湖久久久久久久久免费观看| 国产1区2区3区精品| 亚洲第一青青草原| 亚洲情色 制服丝袜| 另类亚洲欧美激情| 国产片内射在线| 日韩精品免费视频一区二区三区| 2018国产大陆天天弄谢| 亚洲三级黄色毛片| 在线观看三级黄色| 男女无遮挡免费网站观看| 一本色道久久久久久精品综合| 乱人伦中国视频| 婷婷色综合大香蕉| 女的被弄到高潮叫床怎么办| 欧美人与性动交α欧美精品济南到 | 一二三四中文在线观看免费高清| 国产又爽黄色视频| 国产不卡av网站在线观看| 999精品在线视频| 成人亚洲精品一区在线观看| 最近中文字幕2019免费版| 在线观看三级黄色| 日韩制服丝袜自拍偷拍| 边亲边吃奶的免费视频| 大片免费播放器 马上看| 一级爰片在线观看| 免费人妻精品一区二区三区视频| 亚洲第一av免费看| 免费高清在线观看视频在线观看| 91在线精品国自产拍蜜月| 亚洲国产精品999| 男女国产视频网站| 热99久久久久精品小说推荐| 国产片内射在线| av视频免费观看在线观看| av网站在线播放免费| 国产精品麻豆人妻色哟哟久久| av片东京热男人的天堂| 午夜福利视频在线观看免费| 国产男人的电影天堂91| 国产xxxxx性猛交| 各种免费的搞黄视频| av在线app专区| 国产一区二区在线观看av| 国产在线一区二区三区精| 日本wwww免费看| 久久久久久久国产电影| 国产国语露脸激情在线看| 黄色 视频免费看| 欧美人与性动交α欧美精品济南到 | 国产成人欧美| 国产欧美日韩综合在线一区二区| 精品亚洲乱码少妇综合久久| 日韩伦理黄色片| 18禁观看日本| 色婷婷av一区二区三区视频| 国产精品麻豆人妻色哟哟久久| 成年女人毛片免费观看观看9 | 在现免费观看毛片| 久久久久久久久久久久大奶| 久久人人爽人人片av| 男女边摸边吃奶| 在线免费观看不下载黄p国产| 午夜日本视频在线| 妹子高潮喷水视频| 26uuu在线亚洲综合色| 999久久久国产精品视频| 中文字幕人妻熟女乱码| av线在线观看网站| 美女国产视频在线观看| 久久久久久久亚洲中文字幕| 午夜激情久久久久久久| 国产成人午夜福利电影在线观看| 亚洲欧洲精品一区二区精品久久久 | 久久久国产欧美日韩av| 国产成人精品婷婷| 少妇被粗大的猛进出69影院| 在线观看国产h片| 精品人妻在线不人妻| 丰满乱子伦码专区| 欧美激情 高清一区二区三区| 美女xxoo啪啪120秒动态图| 午夜福利,免费看| 人人妻人人添人人爽欧美一区卜| 精品少妇黑人巨大在线播放| 久久韩国三级中文字幕| 国产有黄有色有爽视频| 曰老女人黄片| 视频在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 国产黄频视频在线观看| 欧美成人精品欧美一级黄| 美女主播在线视频| 丝瓜视频免费看黄片| 搡女人真爽免费视频火全软件| 日韩制服丝袜自拍偷拍| av片东京热男人的天堂| 男人添女人高潮全过程视频| √禁漫天堂资源中文www| 精品人妻在线不人妻| av天堂久久9| 欧美+日韩+精品| 最近中文字幕高清免费大全6| 在线天堂中文资源库| 精品少妇内射三级| 寂寞人妻少妇视频99o| 国产 精品1| 视频在线观看一区二区三区| 亚洲精品在线美女| 日本午夜av视频| 欧美 日韩 精品 国产| 免费在线观看黄色视频的| 亚洲av在线观看美女高潮| 亚洲美女视频黄频| 免费黄频网站在线观看国产| 日本猛色少妇xxxxx猛交久久| 80岁老熟妇乱子伦牲交| 看免费av毛片| 夫妻性生交免费视频一级片| 久久国内精品自在自线图片| 亚洲精品视频女| 人妻系列 视频| 69精品国产乱码久久久| 亚洲国产精品国产精品| 日本猛色少妇xxxxx猛交久久| kizo精华| 久久韩国三级中文字幕| 寂寞人妻少妇视频99o| 欧美精品亚洲一区二区| 国产 一区精品| 久久久久久人妻| 成人国语在线视频| av视频免费观看在线观看| 日韩成人av中文字幕在线观看| 久久久久国产网址| av国产久精品久网站免费入址| 国产成人精品久久久久久| 下体分泌物呈黄色| 制服诱惑二区| 亚洲图色成人| 在线精品无人区一区二区三| 亚洲成人手机| 美女脱内裤让男人舔精品视频| 亚洲婷婷狠狠爱综合网| 国产一区二区激情短视频 | 伦理电影大哥的女人| 777米奇影视久久| 超碰97精品在线观看| 男的添女的下面高潮视频| 亚洲欧美中文字幕日韩二区| 18在线观看网站| 蜜桃国产av成人99| 日韩一卡2卡3卡4卡2021年| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲午夜精品一区二区久久| 99久久精品国产国产毛片| 男女下面插进去视频免费观看| 欧美国产精品va在线观看不卡| 亚洲精品aⅴ在线观看| 一边摸一边做爽爽视频免费| 国产熟女午夜一区二区三区| 亚洲国产精品国产精品| 亚洲av福利一区| 午夜免费男女啪啪视频观看| 91在线精品国自产拍蜜月| 波野结衣二区三区在线| 97人妻天天添夜夜摸| 大香蕉久久成人网| 九九爱精品视频在线观看| 成年女人毛片免费观看观看9 | 成人影院久久| 久久鲁丝午夜福利片| 亚洲国产最新在线播放| 久久99一区二区三区| 一级黄片播放器| 哪个播放器可以免费观看大片| 熟妇人妻不卡中文字幕| 成年女人毛片免费观看观看9 | 午夜福利视频精品| 狠狠婷婷综合久久久久久88av| 国产精品久久久久久精品电影小说| 欧美日韩成人在线一区二区| 精品99又大又爽又粗少妇毛片| 可以免费在线观看a视频的电影网站 | 欧美老熟妇乱子伦牲交| 午夜福利视频精品| 免费女性裸体啪啪无遮挡网站| 久久精品夜色国产| 久久精品国产自在天天线| 免费黄频网站在线观看国产| 97在线视频观看| 黄色一级大片看看| 亚洲一区中文字幕在线| 观看av在线不卡| 午夜影院在线不卡| 午夜免费鲁丝| 亚洲欧美一区二区三区黑人 | 日韩精品有码人妻一区| 欧美在线黄色| 美女xxoo啪啪120秒动态图| videosex国产| 九色亚洲精品在线播放| 日韩精品有码人妻一区| h视频一区二区三区| 最近中文字幕高清免费大全6| 精品国产露脸久久av麻豆| 欧美日韩av久久| 波多野结衣一区麻豆| 国产日韩欧美亚洲二区| 国产精品欧美亚洲77777| 91精品三级在线观看| av电影中文网址| 国产福利在线免费观看视频| 香蕉精品网在线| 久久久久视频综合| 伦精品一区二区三区| 一二三四中文在线观看免费高清| 人成视频在线观看免费观看| 妹子高潮喷水视频| 国产成人91sexporn| 国产不卡av网站在线观看| 国产免费一区二区三区四区乱码| 国产精品 国内视频| 久久人人97超碰香蕉20202| 国产成人精品婷婷| 捣出白浆h1v1| 老司机亚洲免费影院| 只有这里有精品99| 涩涩av久久男人的天堂| 激情五月婷婷亚洲| 国产免费现黄频在线看| 日韩 亚洲 欧美在线| 亚洲美女视频黄频| 一区二区日韩欧美中文字幕| 91精品伊人久久大香线蕉| 国产精品 国内视频| 久久久久久久国产电影| 欧美国产精品一级二级三级| 国产高清不卡午夜福利| 日韩中文字幕欧美一区二区 | 久久精品久久久久久噜噜老黄| 久久婷婷青草| 99久久综合免费| 亚洲激情五月婷婷啪啪| 亚洲国产精品国产精品| 亚洲精品国产色婷婷电影| 夫妻性生交免费视频一级片| 一边摸一边做爽爽视频免费| 精品一区二区三区四区五区乱码 | a级毛片黄视频| 国产 一区精品| 日韩视频在线欧美| 免费观看无遮挡的男女| 欧美精品一区二区大全| 亚洲精品国产av蜜桃| 日韩av免费高清视频| 黄色一级大片看看| 男人爽女人下面视频在线观看| 两性夫妻黄色片| 亚洲三级黄色毛片| 亚洲成人手机| 亚洲精品,欧美精品| 国产野战对白在线观看| 亚洲欧洲精品一区二区精品久久久 | 精品一品国产午夜福利视频| 高清不卡的av网站| 亚洲精品日韩在线中文字幕| 伊人久久大香线蕉亚洲五| 夫妻午夜视频| 人体艺术视频欧美日本| 建设人人有责人人尽责人人享有的| 日本vs欧美在线观看视频| 人人妻人人爽人人添夜夜欢视频| 最近手机中文字幕大全| 美女主播在线视频| 精品人妻一区二区三区麻豆| 亚洲精品久久午夜乱码| 天天操日日干夜夜撸| 伊人久久大香线蕉亚洲五| 中国国产av一级| 免费女性裸体啪啪无遮挡网站| 晚上一个人看的免费电影| 欧美国产精品一级二级三级| 人体艺术视频欧美日本| 黄色 视频免费看| 黑人猛操日本美女一级片| 亚洲成色77777| 啦啦啦啦在线视频资源| 在线观看三级黄色| 亚洲成色77777| 亚洲精品在线美女| 国产亚洲欧美精品永久| 老司机影院毛片|