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

    剪切來流下柔性立管渦激振動抑制裝置試驗研究

    2016-05-04 01:45:02付世曉楊家棟
    船舶力學(xué) 2016年4期
    關(guān)鍵詞:渦激螺距立管

    任 鐵,高 云,付世曉,楊家棟,趙 勇

    (1上海交通大學(xué)海洋工程國家重點實驗室,上海200240;2西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點實驗室,成都610500;3大連海事大學(xué)交通運輸裝備與海洋工程學(xué)院,遼寧大連116026)

    剪切來流下柔性立管渦激振動抑制裝置試驗研究

    任 鐵1,高 云2,付世曉1,楊家棟2,趙 勇3

    (1上海交通大學(xué)海洋工程國家重點實驗室,上海200240;2西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點實驗室,成都610500;3大連海事大學(xué)交通運輸裝備與海洋工程學(xué)院,遼寧大連116026)

    文章針對柔性立管螺旋列板抑制裝置在剪切流場中的渦激振動響應(yīng)特性進(jìn)行了試驗研究,試驗過程中通過旋轉(zhuǎn)臂架從而形成相對剪切來流。通過測試得到的應(yīng)變數(shù)據(jù),基于模態(tài)疊加法,可得到立管的位移響應(yīng)等參數(shù)。試驗中針對螺旋列板的螺距和鰭高的變化進(jìn)行了分析,系統(tǒng)地研究了不同螺旋列板狀態(tài)下立管的主導(dǎo)頻率、主導(dǎo)模態(tài)、無量綱振幅比以及疲勞損傷等參數(shù)。研究結(jié)果表明:螺旋列板可以很好地抑制立管的渦激振動響應(yīng);與螺距相比,鰭高對立管渦激振動響應(yīng)會帶來更大的影響;剪切流場中立管螺距恒定為5.0D(D為立管外徑)時,鰭高為0.15D時的立管具有最好的抑制效果。

    剪切來流;柔性立管;渦激振動;螺旋列板

    0 引 言

    立管在一定的來流下,會在立管兩側(cè)形成交替的漩渦,漩渦脫落會產(chǎn)生周期性的橫向(Cross-flow, CF)升力以及流向(In-line,IL)拖曳力。若立管為彈性支撐,周期性的升力以及拖曳力會產(chǎn)生橫向以及流向振動,稱之為渦激振動(Vortex-induced vibration,VIV)[1]。渦激振動會引起立管產(chǎn)生疲勞損傷進(jìn)而導(dǎo)致立管疲勞斷裂引起油氣泄漏,給工程帶來經(jīng)濟損失、以及對環(huán)境帶來生態(tài)污染,因此需要對立管VIV加以抑制。

    抑制立管產(chǎn)生渦激振動的方法主要有兩種[2],一是主動控制,二是被動控制。被動控制由于方便實現(xiàn),在海洋工程領(lǐng)域得到了廣泛的應(yīng)用,而螺旋列板(Helical strakes)裝置則是工程中用來抑制渦激振動普遍采用的抑制裝置之一。螺旋列板的抑制效果分為兩個方面:(1)螺旋列板會在流向方向上破壞正常的漩渦脫落;(2)螺旋列板會降低漩渦脫落在立管軸線方向上的相關(guān)性,從而降低立管的VIV響應(yīng)。理論上來說,如果增加立管上螺旋列板的高度,會減小立管的VIV響應(yīng),但是高度的增加會導(dǎo)致立管在流向的迎水面面積的加大,從而增大立管在流向方向的拖曳力。因此,找到一個合適的螺旋列板,既能抑制VIV響應(yīng)又能最大程度地控制增加的拖曳力,一直是立管海洋工程的研究熱點問題。如表1所示,過去十年的時間內(nèi),國內(nèi)外眾多學(xué)者對此進(jìn)行了研究[3-15]。

    表1 過去十年內(nèi)關(guān)于立管螺旋列板的研究Tab.1 Study on riser fitted with helical strakes in last decade

    由表1可以得到如下結(jié)論:(1)過去大多數(shù)研究立管的細(xì)長比均較小,但是隨著海洋油氣開發(fā)深度的增加,對大細(xì)長比柔性立管的研究需要進(jìn)一步地加強;(2)過去大多研究都是在水中進(jìn)行,這是由于空氣中螺旋列板的VIV抑制效果比水中更為明顯,水中的螺旋列板抑制特性更需要加以關(guān)注;(3)與試驗研究相比,數(shù)值研究要少很多,且考慮的細(xì)長比要小、研究的雷諾數(shù)范圍更窄,這是由計算機的條件所限制;(4)螺旋列板的研究鰭高范圍基本穩(wěn)定在0.10-0.30D(D為立管外徑),螺距范圍在5.0-25.0D,這個范圍是國內(nèi)外眾多學(xué)者多年研究經(jīng)驗所總結(jié)出來的螺旋列板幾何形狀的推薦取值區(qū)間。

    由表1同樣可以看出過去絕大多數(shù)研究均在均勻流下進(jìn)行,但是實際海洋工程環(huán)境中,海水的流速并不是均勻分布的,而是隨著深度的變化流速呈現(xiàn)剪切分布特性。為了更切合實際地研究真實海洋工程下的帶抑制裝置的立管VIV響應(yīng)特性,需要對剪切流場中的帶抑制裝置的立管VIV響應(yīng)特性加以研究。本文以尋找剪切流場中不帶抑制裝置以及帶抑制裝置的柔性立管VIV響應(yīng)特性為研究出發(fā)點,對立管的主導(dǎo)頻率、主導(dǎo)模態(tài)、無量綱振幅比以及疲勞損傷進(jìn)行了系統(tǒng)的研究和討論。

    1 試驗裝置介紹

    如圖1所示,試驗中剪切來流的模擬方法是將立管一端固定在垂直的固定圓柱上,另一端在旋轉(zhuǎn)臂架的帶動下進(jìn)行旋轉(zhuǎn)運動,從而在立管和流場中間形成相對的線性剪切來流。

    圖1 試驗裝置示意圖Fig.1 Sketch of the device during the experiment

    由圖1可以看出:剪流流場中最大流速分布在頂部,最小流速分布在底部為0,沿立管長度方向為線性剪切分布。立管模型的長度為6.75 m,外徑為0.03 m,試驗存在6種流速工況,分別對應(yīng)的立管頂端最大流速為0.8-2.8 m/s,流速間隔為0.4 m/s。立管的其它主要特性參數(shù)見表2所示。

    試驗中立管模型共采用88個光纖光柵應(yīng)變傳感器,分別布置于CF1,CF2,IL1和IL2四個方向(如圖2所示),CF每個方向布置19個傳感器,記為G01-G19;IL每個方向布置25個傳感器,記為F01-F25。

    表2 立管模型主要參數(shù)Tab.2 Key parameters of the riser model

    圖2 應(yīng)變片分布位置示意圖Fig.2 Sketch of the strain sensor locations

    試驗中使用的立管渦激振動抑制裝置是目前海洋工程中應(yīng)用比較廣泛的螺旋列板(Helical Strakes)抑制裝置。螺旋列板的幾何尺寸包括螺距 (pitch)、鰭高(height)以及列板個數(shù) (start number)三個。螺距即螺旋列板繞立管旋轉(zhuǎn)一周在立管軸向的長度,鰭高是螺旋列板的高度,如圖3所示。

    本文試驗研究的列板個數(shù)均取3個,螺距和鰭高采用無量綱來表示,即其與立管直徑的比值。螺距分別取為 5.0D,17.5D以及20.0D,鰭高分別取為0.10D,0.15D以及0.25D,試驗中列板的覆蓋面積為100%,試驗共6種不同的流速工況,立管狀態(tài)共有6種,因此共有36種試驗工況,如表3所示。

    圖3 帶有螺旋列板的立管Fig.3 Riser fitted with helical strakes

    表3 立管試驗工況Tab.3 Test Matrix

    2 數(shù)據(jù)分析方法

    立管發(fā)生渦激振動時,由于立管的周期振動,立管的軸向張力會發(fā)生周期性的變化,這使得測量的應(yīng)變信號包括兩部分:由初始張力產(chǎn)生的軸向應(yīng)變以及由渦激振動產(chǎn)生的軸向應(yīng)變。由預(yù)張力產(chǎn)生的應(yīng)變必須加以消除。圖2中,CF1和CF2相互對稱,因此由VIV產(chǎn)生的彎曲應(yīng)變大小相等,方向相反,而軸向力產(chǎn)生的應(yīng)變是相同的,那么CF1和CF2處的應(yīng)變可寫為:

    對上式進(jìn)行簡單的變換,便可得到橫流方向的由VIV引起的彎曲應(yīng)變?yōu)椋?/p>

    流向方向則與橫向方向情況不同,主要是由于IL方向在初始拖曳力的作用下,立管會在流向產(chǎn)生一個初始的彎曲應(yīng)變εinitial,那么測得的應(yīng)變包括三個部分:由初始張力產(chǎn)生的應(yīng)變、由初始拖曳力產(chǎn)生的應(yīng)變以及由渦激振動產(chǎn)生的軸向應(yīng)變。那么IL1和IL2處的應(yīng)變可寫為:

    為了計算(3)式,如果試驗選取的穩(wěn)定段時間足夠長,可認(rèn)為渦激振動引起的彎曲應(yīng)變的時間歷程均值為零,引入假設(shè)1,可表示如下:

    IL方向初始拖曳力引起的初始應(yīng)變的變化只可能由拖車速度發(fā)生略微的波動所導(dǎo)致,但是我們拖車速度的精確度達(dá)到了0.2%,因此可以引入假設(shè)2,假設(shè)初始應(yīng)變不隨時間發(fā)生變化,表示如下:

    假設(shè)受軸向力作用的立管做小變形的振動,則CF和IL方向的響應(yīng)均可基于模態(tài)疊加法寫為:

    其中:φi(z)是立管的第i階模態(tài)振型,pi(t)為立管的第i階模態(tài)位移權(quán)重,z為立管的位置,基于小變形假設(shè),立管的曲率可表示為立管位移響應(yīng)對空間求二次導(dǎo)數(shù):

    其中:φi″(z)為立管的第i階模態(tài)曲率,立管的曲率和彎曲應(yīng)變之間的關(guān)系,可表示為:

    式中:R為立管的外部半徑,由于(10)和(11)式,可以看出:給定了測點的應(yīng)變可求出對應(yīng)的模態(tài)權(quán)重,進(jìn)一步根據(jù)(9)式便可以求出位移響應(yīng)。本文試驗的立管模型可以簡化為兩端鉸接的索模型,因此第i階模態(tài)振型可寫成:

    依據(jù)測試得到的應(yīng)變數(shù)據(jù)便可根據(jù)(15)式計算得到應(yīng)變模態(tài)權(quán)重,并進(jìn)一步可通過(16)式得到位移模態(tài)權(quán)重,最后便可由(9)式得到結(jié)構(gòu)的位移響應(yīng)。由于試驗中原始應(yīng)變信號中含有噪聲信號,因此需要對原始應(yīng)變信號進(jìn)行濾波處理,把原始信號中頻率低于1 Hz以及高于50 Hz成分的信號進(jìn)行消除。

    3 分析與討論

    在對流速V進(jìn)行分析時采用了無量綱化,引入折合速度,定義如下:

    式中:V為流速,D為立管外徑,f1為立管在水中的一階固有頻率,可由下式計算得到:

    實際上F(t)是隨時間在預(yù)張力附近發(fā)生變化的,為數(shù)據(jù)處理方便,此處F(t)取為預(yù)張力3 000 N,為定值。m為單位長度振動系統(tǒng)質(zhì)量,包括立管結(jié)構(gòu)質(zhì)量、內(nèi)部流體質(zhì)量以及附加質(zhì)量,這里附加質(zhì)量系數(shù)取為1.0,l為立管長度,E為彈性模量,I為立管的慣性矩,這些參數(shù)均可由表2獲得,n為立管固有頻率的階數(shù)。經(jīng)過(18)式計算得到立管的一階固有頻率為2.3 Hz,已知立管的一階固有頻率,立管直徑以及流速,便可根據(jù)(17)式計算得到6種流速對應(yīng)的6種折合速度,如表4所示。

    這里為了進(jìn)一步驗證立管固有頻率計算的可靠性,采用了試驗方法對立管在水中的真實固有頻率進(jìn)行了研究。通過對水中立管中點處的應(yīng)變時間歷程做快速傅立葉變換(Fast Fourier Transform, FFT),便可得到立真實固有頻率。如圖4所示,給出了水中立管中點處的應(yīng)變時歷曲線以及應(yīng)變幅值譜,由圖可以看出經(jīng)試驗測試得到的一階固有頻率為2.41 Hz,與計算固有頻率2.3 Hz較接近。

    表4 不同流速下對應(yīng)的相對折合速度Tab.4 The relative reduced velocity versus different velocity

    圖4 衰減試驗中立管中點處的應(yīng)變時間歷程曲線以及應(yīng)變幅值譜Fig.4 Time trace and amplitude spectrum of the CF strain at midpoint during free decay test

    圖4同樣給出了立管中點處第N個以及N+M個立管的應(yīng)變幅值,便可根據(jù)下式計算得到立管在水中的阻尼比。

    經(jīng)計算得到的水中阻尼比為0.014,如表2所示。

    3.1 鰭高對立管渦激振動響應(yīng)特性的影響

    為了研究方便,這里分別以Strake/5.0D/0.10D、Strake/5.0D/0.15D以及Strake/5.0D/0.25D表示螺距為5.0D,鰭高依次為0.10D、0.15D以及0.25D時的立管。圖5給出了裸管以及不同鰭高立管狀態(tài)下CF以及IL方向的主導(dǎo)頻率,表5給出了對應(yīng)的主導(dǎo)模態(tài)。由圖5可以看出:剪切流場中,裸管的主導(dǎo)頻率不再出現(xiàn)在固定圓柱體尾流漩渦泄放頻率(fst=St×U/D)附近,且流向主導(dǎo)頻率與橫向主導(dǎo)頻率不再存在2倍關(guān)系。

    由圖5和表5可以看出:與裸管相比,帶strake立管的主導(dǎo)頻率以及主導(dǎo)模態(tài)均得到了大幅度地降低。CF方向的主導(dǎo)模態(tài)由4階降為2階,IL方向的主導(dǎo)模態(tài)由5階降為1階。與CF方向相比,IL方向模樣階數(shù)的降低程度更為明顯。

    圖5 不同鰭高立管狀態(tài)下對應(yīng)的CF和IL方向的主導(dǎo)頻率Fig.5 Dominant frequency versus reduced velocity for risers with different heights in both CF and IL directions

    表5 不同鰭高立管狀態(tài)下對應(yīng)的CF和IL方向的主導(dǎo)模態(tài)Tab.5 Dominant mode versus reduced velocity for risers with different heights in both CF and IL directions

    圖6 不同鰭高立管狀態(tài)下對應(yīng)的CF和IL方向的無量綱位移標(biāo)準(zhǔn)差最大值Fig.6 Maximum RMS value of non-dimensional displacement versus reduced velocity for risers with different heights in both CF and IL directions

    圖6給出了不同鰭高狀態(tài)下立管的CF以及IL方向的最大無量綱振幅比,可以看出裸管CF方向的最大值為0.95D,發(fā)生在Vr=29.0處;IL方向的最大值為0.33D,發(fā)生在Vr=40.6處;對于裸管,CF方向最大值大約是IL方向的3倍。與裸管相比,加上strake抑制裝置后,無論是CF還是IL方向,無量綱振幅比均得到了大幅度降低。無量綱振幅比并沒有隨著鰭高的增加而一直上升:當(dāng)鰭高由0.10D增加到0.15D時,無論是CF方向還是IL方向,無量綱振幅比均成下降趨勢;當(dāng)鰭高繼續(xù)增加,由0.15D增加到0.25D時,無量綱振幅比隨著鰭高的增加而上升。

    這里為了進(jìn)一步研究strake幾何形狀對立管VIV響應(yīng)的影響,我們進(jìn)一步對立管的疲勞損傷進(jìn)行了分析。如果立管上的應(yīng)力分布滿足Rayleigh分布,那么疲勞損傷便可表示為:

    式中:Srms是某個特定點的應(yīng)力均方根值,可以通過Srms=εrms×E計算得到,其中E為彈性模量,如表2所示,εrms為應(yīng)變均方根值,由試驗測試得到。Ty為一年內(nèi)的時間(以秒計),b和c分別為材料S-N曲線中對應(yīng)的材料常數(shù)。

    為了更為準(zhǔn)確地計算剪切流中的VIV引起的疲勞損傷,我們引入了跨零頻率fz,其表達(dá)如下:

    式中:S(f)為應(yīng)力響應(yīng)頻率譜(或應(yīng)變響應(yīng)譜,二者差一個常系數(shù)彈性模量),該譜的最大值處對應(yīng)的頻率就是我們通常所說的譜峰頻率fp,若計算均勻流,由于均勻流流場隨時間變化非常穩(wěn)定,譜能量分布比較集中,在計算疲勞損傷時,可以將譜峰頻率作為計算頻率[12];但是當(dāng)來流為復(fù)雜剪切來流時,應(yīng)力響應(yīng)頻率非常分散,通常存在若干個與峰值頻率相當(dāng)?shù)念l率,因此光憑譜峰頻率無法準(zhǔn)確判斷應(yīng)力響應(yīng)譜的平均能量分布范圍,因此在這里我們引入了跨零頻率。

    圖7 不同鰭高立管狀態(tài)下對應(yīng)的CF和IL方向的最大疲勞損傷Fig.7 Maximum fatigue damage versus reduced velocity for risers with different heights in both CF and IL directions

    計算疲勞損傷時對所有測點(CF方向19個測點,IL方向25個測點)均進(jìn)行了計算,然后取所有測點處的最大疲勞損傷。圖7給出了CF以及IL方向下不同鰭高立管狀態(tài)下的疲勞損傷。由圖可以看出:當(dāng)采用了strake后,無論是CF還是IL方向疲勞損傷均得到了大幅度地降低。對于CF方向,當(dāng)Vr小于23.2時,0.15D立管具有最小的疲勞損傷;當(dāng)Vr超過23.2時,0.25D立管具有最低疲勞損傷。對于IL方向,所有折合速度下,0.15D立管均具有最低的疲勞損傷。

    3.2 螺距對立管渦激振動響應(yīng)特性的影響

    圖8和表6分別給出了不同螺距立管狀態(tài)下主導(dǎo)頻率以及主導(dǎo)模態(tài),由圖8可看出:strake可以很好地抑制立管的主導(dǎo)頻率;表6可看出:裸管CF方向的最大主導(dǎo)模態(tài)為4階,IL方向的最大主導(dǎo)模態(tài)為5階;加上strake抑制裝置后,CF方向以及IL方向的主導(dǎo)模態(tài)分別降為1階和2階,這表明了strake可以有效地抑制立管的主導(dǎo)模態(tài)。

    圖8 不同螺距立管狀態(tài)下對應(yīng)的CF和IL方向的主導(dǎo)頻率Fig.8 Dominant frequency versus reduced velocity for risers with different pitches in both CF and IL directions

    表6 不同螺距立管狀態(tài)下對應(yīng)的CF和IL方向的主導(dǎo)模態(tài)Tab.6 Dominant mode versus reduced velocity for risers with different pitches in both CF and IL directions

    圖9和圖10分別給出了不同螺距立管狀態(tài)下CF以及IL方向的無量綱位移標(biāo)準(zhǔn)差最大值以及最大疲勞損傷。由圖9和圖10可以看出:當(dāng)鰭高一定(0.25D)時,無論是CF還是IL方向,螺距為5.0D的立管均具有最大的無量綱振幅比;螺距為17.5D和20.0D的立管的最大無量綱振幅比很接近,均小于5.0D立管。當(dāng)螺距變化時,疲勞損傷變化非常不明顯,三種螺距下,5.0D立管具有最大疲勞損傷;而17.5D與20.0D立管疲勞損傷接近,略小于5.0D立管。

    圖9 不同螺距立管狀態(tài)下對應(yīng)的CF和IL方向的無量綱位移標(biāo)準(zhǔn)差最大值Fig.9 Maximum RMS value of non-dimensional displacement versus reduced velocity for risers with different pitches in both CF and IL directions

    圖10 不同螺距立管狀態(tài)下對應(yīng)的CF和IL方向的最大疲勞損傷Fig.10 Maximum fatigue damage versus reduced velocity for risers with different pitches in both CF and IL directions

    4 結(jié) 語

    本文針對裸管以及帶有螺旋列板的柔性立管在剪切來流下的渦激振動響應(yīng)特性進(jìn)行了研究,并系統(tǒng)地研究了不同鰭高以及螺距的strake立管狀態(tài)下的響應(yīng)特性。特性參數(shù)包括應(yīng)變、位移、主導(dǎo)頻率、主導(dǎo)模態(tài)以及疲勞損傷等參數(shù),注意在計算疲勞損傷時并沒有使用通常在均勻來流計算中采用的峰值頻率,而是采用了更能反映應(yīng)力分布平均能量的跨零頻率。通過以上研究可得到如下結(jié)論:

    (1)對于剪切流場中的裸管,CF以及IL方向的主導(dǎo)頻率不再存在通常均勻流場中的2倍關(guān)系,且CF以及IL方向的主導(dǎo)頻率不再分布在斯脫哈爾漩渦泄放頻率附近。

    (2)剪切流場中,當(dāng)裸管采用了strake抑制裝置后,CF以及IL方向的主導(dǎo)頻率、主導(dǎo)模態(tài)、無量綱振幅比以及疲勞損傷均得到了大幅地降低,這說明了strake立管可以很好地抑制立管的VIV響應(yīng)。

    (3)由不同立管狀態(tài)下的疲勞損傷可以看出:當(dāng)立管螺距一定、鰭高發(fā)生變化時,立管的疲勞損傷發(fā)生較大的變化;但當(dāng)立管鰭高一定、螺距發(fā)生變化時候,立管的疲勞損傷發(fā)生很小的變化。這說明了與螺距變化相比,鰭高變化給立管疲勞損傷帶來的變化更大。

    (4)在螺距一定(5.0D)時,鰭高為0.15D的立管在三種立管狀態(tài)中(0.10D,0.15D,0.25D)具有最低的疲勞損傷;當(dāng)鰭高一定(0.25D)時,無論是無量綱振幅比還是疲勞損傷,螺距為5.0D時的立管均具有最大的疲勞損傷,而螺距為17.5D以及20.0D的立管,疲勞損傷相當(dāng),均小于5.0D立管。

    [1]Blevins R D.Flow-Induced Vibration[M].2nd edn.Krieger Publishing,Inc.,Malabar/Florida,USA,2001.

    [2]Choi H,Jeon W P,Kim J.Control of flow over a bluff body[J].Annual Review of Fluid Mechanics,2008,40:113-139.

    [3]Brankovi M,Bearman P W.Measurements of transverse forces on circular cylinders undergoing vortex-induced vibration [J].Journal of Fluids and Structures,2006,22:829-836.

    [4]Vandiver J K,Swithenbank S,Jaiswai V.The effectiveness of helical strakes in the suppression of high-mode-number VIV [C].Offshore Technology Conference,Paper 18276,2006.

    [5]Allen D W,Henning D L,Lee L.High Reynolds number flow tests of flexible cylinders with helical strakes[C].Proceedings of the OMAE 2006 25th International Conference on Offshore Mechanics and Arctic Engineering,Paper 92332,2006.

    [6]Korkischko I,Meneghini J R,Gioria R S,Jabardo P J.Casaprima,E.,Franciss,R.:An experimental investigation of the flow around straked cylinders[C].Proceedings of the OMAE 2007 26th International Conference on Offshore Mechanics and Arctic Engineering,Paper 29057,2007.

    [7]Hao Z,Zhou T,Wang X,Tan S K.Experimental studies of vortex structures in the wake of a cylinder with helical strakes [C].Proceedings of the ASME 2010 29th International Conference on Ocean,Offshore and Arctic Engineering,Paper 20181, 2010.

    [8]Korischko I,Meneghini J R.Experimental investigation of flow-induced vibration on isolated and tandem circular cylinders fitted with strakes[J].Journal of Fluids and Structures,2006,26:611-625.

    [9]Zhou T,Razali S F M,Hao Z,Cheng L.On the study of vortex-induced vibration of a cylinder with helical strakes[J].Journal of Fluids and Structures,2011,27:903-917.

    [10]Gao Y,Fu S X,Song L J,Peng T,Li R P.Experimental investigation on the suppression device of VIV of a flexible riser [C].Proceedings of the ASME 2014 33rd International Conference on Ocean,Offshore and Arctic Engineering,Paper 23427, 2014.

    [11]Gao Y,Fu S X,Cao J,Chen Y F.Experimental study on response performance of VIV of a flexible riser with helical strakes[J].China Ocean Engineering,2015,29(5):673 690.

    [12]Constantinides Y,Oakley O H.Numerical prediction of bare and straked cylinder VIV[C].Proceedings of the OMAE 2006 25th International Conference on Offshore Mechanics and Arctic Engineering,Paper 92334,2006.

    [13]Pinto A,Broglia R,Ciappi E,Mascio A D,Campana E F,Racco P.Vortex suppression efficiency of discontinuous helicoidal fins[C].Proceedings of the OMAE 2007 26th International Conference on Offshore Mechanics and Arctic Engineering,Paper 29255,2007.

    [14]Carmo B S,Gioria R S,Korkischko I,Freire C M,Meneghini J R.Two-and three-dimensional simulations of the flow around a cylinder fitted with strake[C].Proceedings of the ASME 2012 31st International Conference on Ocean,Offshore and Arctic Engineering,Paper 83603,2012.

    [15]Zhu R Q,Zheng T T,Ni Y Y.CFD analysis of the efficiency of helical strakes in suppressing vortex induced vibrations[C]. Proceedings of the ASME 2013 32nd International Conference on Ocean,Offshore and Arctic Engineering,Paper 11494, 2013.

    Experimental investigation on the suppression device of VIV of a flexible riser in sheared currents

    REN Tie1,GAO Yun2,FU Shi-xiao1,YANG Jia-dong2,ZHAO Yong3
    (1 State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China; 2 State Key Laboratory of Oil and Gas Reservoir Geology and Exploration,Southwest Petroleum University, Chengdu 610500,China;3 Transportation Equipment and Ocean Engineering College,Dalian Maritime University,Dalian 116026,China)

    Laboratory test was conducted on the flexible riser with helical strakes in order to understand its response characteristics in sheared currents.The relative sheared currents were simulated by rotating the cantilever during the experiment.Based on the modal analysis method,the displacement responses can be obtained by the measured strain.The strakes with different heights and pitches were analyzed,and the response parameters such as dominant frequency,dominant mode,non-dimensional displacement and fatigue damage were studied.The analysis results indicate that helical strakes can suppress vortex induced vibration(VIV)response effectively.Compared with strake pitch,strake height has greater influence on the VIV response.When the given pitch is 5.0D(D is the external diameter of the riser)in the sheared current,the 0.15D height always has the best suppression efficiency.

    sheared currents;flexible riser;Vortex-Induced Vibration;helical strakes

    O357

    :Adoi:10.3969/j.issn.1007-7294.2016.04.014

    1007-7294(2016)04-0497-11

    2015-12-22

    國家自然科學(xué)基金(51239007,51490674,51490675,51279101);西南石油大學(xué)科研啟航計劃項目(2015QHZ005);西南石油大學(xué)青年教師“過學(xué)術(shù)關(guān)”基金(201499010114)

    任 鐵(1985-),男,博士研究生;高 云(1985-),男,講師,碩士生導(dǎo)師,E-mail:dutgaoyun@163.com;付世曉(1976-),男,研究員,博士生導(dǎo)師,E-mail:shixiao.fu@sjtu.edu.cn。

    猜你喜歡
    渦激螺距立管
    不同間距比下串聯(lián)圓柱渦激振動數(shù)值模擬研究
    變螺距螺紋數(shù)控車削技巧
    渦激振動發(fā)電裝置及其關(guān)鍵技術(shù)
    常見高層建筑物室內(nèi)給水立管材質(zhì)解析
    盤球立管結(jié)構(gòu)抑制渦激振動的數(shù)值分析方法研究
    電子制作(2018年14期)2018-08-21 01:38:42
    常用粗牙螺紋螺距簡易計算公式及應(yīng)用
    機械工程師(2018年7期)2018-07-30 09:05:26
    深水鋼懸鏈立管J型鋪設(shè)研究
    雙源CT大螺距掃描在胸部篩查的應(yīng)用
    柔性圓管在渦激振動下的模態(tài)響應(yīng)分析
    The Power of Integration
    Beijing Review(2015年43期)2015-11-25 03:12:04
    午夜福利视频1000在线观看| 中亚洲国语对白在线视频| 老熟妇仑乱视频hdxx| 99久久99久久久精品蜜桃| 日本撒尿小便嘘嘘汇集6| 色播亚洲综合网| 91成人精品电影| 欧美精品啪啪一区二区三区| 午夜久久久久精精品| 国产国语露脸激情在线看| 国产精品久久久人人做人人爽| 不卡一级毛片| 我的亚洲天堂| 亚洲人成伊人成综合网2020| 亚洲精华国产精华精| 女性生殖器流出的白浆| 亚洲人成77777在线视频| 亚洲欧美一区二区三区黑人| 青草久久国产| 欧美丝袜亚洲另类 | 免费在线观看日本一区| 啦啦啦观看免费观看视频高清| 欧美激情久久久久久爽电影| 一二三四在线观看免费中文在| 超碰成人久久| 亚洲色图av天堂| 国产亚洲av高清不卡| 国产精品亚洲一级av第二区| 国产精品98久久久久久宅男小说| 亚洲人成网站在线播放欧美日韩| 人人妻人人澡欧美一区二区| 男人舔女人下体高潮全视频| 中文字幕人妻丝袜一区二区| 久久中文看片网| 啦啦啦 在线观看视频| 国产精品九九99| 久久婷婷人人爽人人干人人爱| 久热爱精品视频在线9| 婷婷精品国产亚洲av| 91字幕亚洲| 国产一区在线观看成人免费| 一a级毛片在线观看| 91成人精品电影| 久久久久久九九精品二区国产 | 夜夜躁狠狠躁天天躁| 91大片在线观看| 黄色a级毛片大全视频| 窝窝影院91人妻| 韩国av一区二区三区四区| 欧美绝顶高潮抽搐喷水| 又黄又粗又硬又大视频| 黄色a级毛片大全视频| 我的亚洲天堂| 麻豆一二三区av精品| 99国产精品99久久久久| 我的亚洲天堂| 我的亚洲天堂| 久久午夜综合久久蜜桃| 亚洲天堂国产精品一区在线| 欧美大码av| 日韩一卡2卡3卡4卡2021年| av片东京热男人的天堂| 久久久久国内视频| 国产午夜精品久久久久久| 男人舔奶头视频| 亚洲精品国产区一区二| 一进一出抽搐gif免费好疼| 麻豆成人av在线观看| 亚洲一区二区三区不卡视频| 成人亚洲精品av一区二区| 国产一区二区三区在线臀色熟女| 午夜福利高清视频| www.自偷自拍.com| av天堂在线播放| 久久精品aⅴ一区二区三区四区| 日本精品一区二区三区蜜桃| 成人免费观看视频高清| 久久中文看片网| 亚洲国产欧洲综合997久久, | 国产亚洲av嫩草精品影院| 校园春色视频在线观看| 9191精品国产免费久久| 亚洲午夜理论影院| 亚洲精品av麻豆狂野| 亚洲av中文字字幕乱码综合 | 国产免费av片在线观看野外av| 国产精品影院久久| 女同久久另类99精品国产91| 麻豆成人av在线观看| 久久久久免费精品人妻一区二区 | 曰老女人黄片| 亚洲av五月六月丁香网| av免费在线观看网站| 久久久久免费精品人妻一区二区 | 两个人看的免费小视频| 国产v大片淫在线免费观看| 性欧美人与动物交配| 色av中文字幕| 成人亚洲精品av一区二区| 亚洲精品久久国产高清桃花| 亚洲av熟女| 亚洲九九香蕉| 淫秽高清视频在线观看| 亚洲熟妇中文字幕五十中出| 草草在线视频免费看| 一a级毛片在线观看| 久久欧美精品欧美久久欧美| 日本黄色视频三级网站网址| 曰老女人黄片| 亚洲 欧美一区二区三区| 亚洲av片天天在线观看| 亚洲人成网站高清观看| 亚洲久久久国产精品| 99国产极品粉嫩在线观看| 麻豆成人av在线观看| 激情在线观看视频在线高清| 久久精品国产综合久久久| 亚洲 欧美一区二区三区| 久久人妻av系列| 91麻豆av在线| 久99久视频精品免费| 在线观看免费日韩欧美大片| 夜夜爽天天搞| 一夜夜www| 俄罗斯特黄特色一大片| 亚洲av第一区精品v没综合| av有码第一页| 18禁黄网站禁片免费观看直播| 人人妻人人澡欧美一区二区| 黄片小视频在线播放| 亚洲熟妇熟女久久| 一二三四社区在线视频社区8| 美女高潮喷水抽搐中文字幕| 日韩视频一区二区在线观看| 母亲3免费完整高清在线观看| 男人操女人黄网站| 免费搜索国产男女视频| 午夜亚洲福利在线播放| 国产精品精品国产色婷婷| www.999成人在线观看| 免费观看人在逋| 给我免费播放毛片高清在线观看| 99热这里只有精品一区 | 9191精品国产免费久久| 久久国产精品影院| ponron亚洲| 国产一卡二卡三卡精品| 十八禁人妻一区二区| 在线av久久热| 精品久久久久久久久久久久久 | 日韩免费av在线播放| 午夜福利在线在线| 亚洲国产毛片av蜜桃av| 99在线人妻在线中文字幕| 在线观看午夜福利视频| 女人被狂操c到高潮| 国产精品一区二区免费欧美| 亚洲精品国产精品久久久不卡| 热99re8久久精品国产| 国产精品美女特级片免费视频播放器 | 美女午夜性视频免费| 宅男免费午夜| 亚洲无线在线观看| 极品教师在线免费播放| 夜夜爽天天搞| 成人手机av| 国产国语露脸激情在线看| 夜夜夜夜夜久久久久| 欧美大码av| 欧美日本视频| 亚洲中文字幕日韩| 精品久久久久久,| 国产成人系列免费观看| 美女免费视频网站| 欧美激情 高清一区二区三区| www日本黄色视频网| 国产在线精品亚洲第一网站| 成人精品一区二区免费| 国产视频一区二区在线看| 日本免费一区二区三区高清不卡| 女人高潮潮喷娇喘18禁视频| 宅男免费午夜| 国产黄色小视频在线观看| 最近在线观看免费完整版| 99热这里只有精品一区 | 欧美一级a爱片免费观看看 | 国内久久婷婷六月综合欲色啪| 高潮久久久久久久久久久不卡| 精品福利观看| 国产av在哪里看| 日韩欧美三级三区| 亚洲av熟女| 久久香蕉激情| 成人亚洲精品av一区二区| 精品一区二区三区av网在线观看| 久久青草综合色| 天天一区二区日本电影三级| 国产午夜福利久久久久久| 男女下面进入的视频免费午夜 | 国产高清有码在线观看视频 | 久久精品人妻少妇| 一级a爱片免费观看的视频| av欧美777| 亚洲,欧美精品.| 国产三级黄色录像| 国产伦人伦偷精品视频| 亚洲五月天丁香| 人人妻人人澡人人看| 欧美成人免费av一区二区三区| 91大片在线观看| 女性被躁到高潮视频| 黄色成人免费大全| 女同久久另类99精品国产91| 麻豆久久精品国产亚洲av| 国内少妇人妻偷人精品xxx网站 | 视频在线观看一区二区三区| 久久久国产成人免费| 成年免费大片在线观看| 色av中文字幕| 亚洲人成网站高清观看| 天天一区二区日本电影三级| 欧美一级毛片孕妇| 欧美最黄视频在线播放免费| 欧美乱码精品一区二区三区| 精品不卡国产一区二区三区| 91字幕亚洲| 99热只有精品国产| 91av网站免费观看| 亚洲成a人片在线一区二区| 精品久久久久久久久久免费视频| 国产乱人伦免费视频| 999久久久国产精品视频| 久久久国产欧美日韩av| 国产主播在线观看一区二区| 天堂√8在线中文| 中国美女看黄片| 精品欧美一区二区三区在线| 9191精品国产免费久久| 侵犯人妻中文字幕一二三四区| 好男人电影高清在线观看| 他把我摸到了高潮在线观看| 69av精品久久久久久| 中文字幕人妻熟女乱码| 最新美女视频免费是黄的| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品色激情综合| 国产精品爽爽va在线观看网站 | 精品日产1卡2卡| 精品国产国语对白av| 大型黄色视频在线免费观看| 国内揄拍国产精品人妻在线 | 日韩视频一区二区在线观看| 久久精品成人免费网站| 法律面前人人平等表现在哪些方面| 怎么达到女性高潮| 99国产精品一区二区三区| 亚洲国产高清在线一区二区三 | 亚洲aⅴ乱码一区二区在线播放 | 一本久久中文字幕| 欧美丝袜亚洲另类 | 日本一本二区三区精品| 777久久人妻少妇嫩草av网站| 国产不卡一卡二| 日本a在线网址| 色综合欧美亚洲国产小说| 麻豆一二三区av精品| 搞女人的毛片| 国产亚洲精品久久久久久毛片| av中文乱码字幕在线| 中文字幕人妻丝袜一区二区| 亚洲av熟女| 久久精品国产99精品国产亚洲性色| 一区二区三区国产精品乱码| 在线观看www视频免费| 波多野结衣高清无吗| 丰满的人妻完整版| 日韩欧美免费精品| 男人操女人黄网站| 国产精品 欧美亚洲| 欧美乱码精品一区二区三区| 身体一侧抽搐| 国产欧美日韩一区二区三| 搡老妇女老女人老熟妇| 很黄的视频免费| 搡老熟女国产l中国老女人| ponron亚洲| 亚洲国产看品久久| 女性生殖器流出的白浆| 色精品久久人妻99蜜桃| 脱女人内裤的视频| 熟女电影av网| 色在线成人网| 亚洲熟妇熟女久久| 国产精品久久视频播放| 97超级碰碰碰精品色视频在线观看| 国产aⅴ精品一区二区三区波| 天天一区二区日本电影三级| 久久久久久人人人人人| 午夜福利一区二区在线看| 亚洲aⅴ乱码一区二区在线播放 | 母亲3免费完整高清在线观看| 亚洲成国产人片在线观看| 美国免费a级毛片| 男人舔奶头视频| 久久久久久国产a免费观看| 日日干狠狠操夜夜爽| 国产麻豆成人av免费视频| 91老司机精品| 一二三四在线观看免费中文在| 色综合亚洲欧美另类图片| 免费看十八禁软件| 日韩欧美国产在线观看| 成人精品一区二区免费| 99久久久亚洲精品蜜臀av| 免费高清视频大片| 1024手机看黄色片| 亚洲av中文字字幕乱码综合 | 久久精品成人免费网站| 老司机福利观看| 国产成+人综合+亚洲专区| 夜夜看夜夜爽夜夜摸| 美女大奶头视频| 午夜福利免费观看在线| 国产精品久久视频播放| 午夜福利成人在线免费观看| 国产亚洲av高清不卡| 亚洲成av片中文字幕在线观看| 国产真实乱freesex| 亚洲中文日韩欧美视频| 午夜福利在线在线| 久久中文字幕人妻熟女| 麻豆国产av国片精品| 两个人看的免费小视频| 国产精品98久久久久久宅男小说| 一a级毛片在线观看| 狠狠狠狠99中文字幕| 色尼玛亚洲综合影院| 一级a爱视频在线免费观看| 免费av毛片视频| 欧美色视频一区免费| 国产成人精品久久二区二区免费| 精品人妻1区二区| 香蕉丝袜av| 欧美中文综合在线视频| 精品国产超薄肉色丝袜足j| 91麻豆精品激情在线观看国产| 亚洲精品久久国产高清桃花| 国产激情偷乱视频一区二区| 国产野战对白在线观看| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| 中文字幕精品免费在线观看视频| 国产精品香港三级国产av潘金莲| 亚洲五月天丁香| 日韩大尺度精品在线看网址| 国产三级在线视频| 亚洲五月天丁香| 免费在线观看日本一区| 麻豆久久精品国产亚洲av| 国产高清激情床上av| 亚洲av成人不卡在线观看播放网| 搞女人的毛片| 色在线成人网| 男女之事视频高清在线观看| 欧美黄色淫秽网站| 老司机深夜福利视频在线观看| 亚洲国产毛片av蜜桃av| 午夜亚洲福利在线播放| 国产激情久久老熟女| 国内毛片毛片毛片毛片毛片| 两性夫妻黄色片| 在线国产一区二区在线| 午夜福利18| 香蕉丝袜av| 国产单亲对白刺激| 免费在线观看视频国产中文字幕亚洲| 亚洲精品一区av在线观看| av在线天堂中文字幕| 精品无人区乱码1区二区| 久久精品国产综合久久久| 国产亚洲精品第一综合不卡| 99在线人妻在线中文字幕| bbb黄色大片| 成人永久免费在线观看视频| 久久精品aⅴ一区二区三区四区| 两性夫妻黄色片| 韩国av一区二区三区四区| 欧美+亚洲+日韩+国产| 亚洲av电影在线进入| 久久久精品国产亚洲av高清涩受| 91在线观看av| 久久99热这里只有精品18| 欧美人与性动交α欧美精品济南到| 黑人操中国人逼视频| 久久天堂一区二区三区四区| 99热6这里只有精品| 757午夜福利合集在线观看| 国产精品香港三级国产av潘金莲| 国产色视频综合| 好看av亚洲va欧美ⅴa在| 亚洲第一青青草原| 热99re8久久精品国产| 亚洲一区中文字幕在线| 免费高清视频大片| 免费在线观看成人毛片| 激情在线观看视频在线高清| 色综合亚洲欧美另类图片| 亚洲精品久久国产高清桃花| 制服人妻中文乱码| 99国产极品粉嫩在线观看| 一边摸一边抽搐一进一小说| 看片在线看免费视频| 男女下面进入的视频免费午夜 | 十分钟在线观看高清视频www| 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| www.www免费av| 真人做人爱边吃奶动态| 精品久久久久久,| 中国美女看黄片| 精品一区二区三区av网在线观看| 日韩欧美国产一区二区入口| 午夜亚洲福利在线播放| 人人澡人人妻人| 美女大奶头视频| 欧美人与性动交α欧美精品济南到| 国产免费男女视频| 日韩欧美一区二区三区在线观看| 99在线视频只有这里精品首页| 757午夜福利合集在线观看| 老汉色av国产亚洲站长工具| 成人亚洲精品av一区二区| 夜夜看夜夜爽夜夜摸| 国产精品日韩av在线免费观看| 午夜成年电影在线免费观看| www.自偷自拍.com| 视频在线观看一区二区三区| 久久精品国产综合久久久| 国产成年人精品一区二区| 成年女人毛片免费观看观看9| 侵犯人妻中文字幕一二三四区| 91字幕亚洲| 国产午夜精品久久久久久| e午夜精品久久久久久久| 欧美精品亚洲一区二区| 99精品欧美一区二区三区四区| 国产成人一区二区三区免费视频网站| 在线十欧美十亚洲十日本专区| 久久99热这里只有精品18| 制服诱惑二区| 色婷婷久久久亚洲欧美| 亚洲熟妇熟女久久| 伦理电影免费视频| 亚洲自拍偷在线| 国产成人精品无人区| 99久久综合精品五月天人人| 少妇 在线观看| 国产精品日韩av在线免费观看| 正在播放国产对白刺激| 人人妻,人人澡人人爽秒播| 国产黄色小视频在线观看| 正在播放国产对白刺激| 亚洲全国av大片| 69av精品久久久久久| 人妻丰满熟妇av一区二区三区| 亚洲美女黄片视频| 高清毛片免费观看视频网站| 国产激情偷乱视频一区二区| 久久国产精品男人的天堂亚洲| 一区二区三区精品91| 91九色精品人成在线观看| 国产欧美日韩一区二区三| 日本一区二区免费在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩欧美在线二视频| 国内少妇人妻偷人精品xxx网站 | 老司机午夜福利在线观看视频| 在线国产一区二区在线| 亚洲国产精品999在线| 热99re8久久精品国产| 久久香蕉精品热| 日韩视频一区二区在线观看| 我的亚洲天堂| 亚洲成人国产一区在线观看| 老司机福利观看| 成人av一区二区三区在线看| 最近最新中文字幕大全免费视频| 啦啦啦观看免费观看视频高清| 97人妻精品一区二区三区麻豆 | 欧美另类亚洲清纯唯美| 最好的美女福利视频网| 真人一进一出gif抽搐免费| 亚洲最大成人中文| 美女高潮喷水抽搐中文字幕| 国产国语露脸激情在线看| 午夜激情av网站| 女性生殖器流出的白浆| 天堂影院成人在线观看| 亚洲欧美日韩无卡精品| 中文字幕精品亚洲无线码一区 | 午夜福利欧美成人| 在线看三级毛片| 亚洲国产精品久久男人天堂| 午夜两性在线视频| or卡值多少钱| 欧美国产日韩亚洲一区| 伊人久久大香线蕉亚洲五| 久久久久九九精品影院| 激情在线观看视频在线高清| 欧美日韩中文字幕国产精品一区二区三区| 色尼玛亚洲综合影院| 巨乳人妻的诱惑在线观看| 免费看日本二区| 波多野结衣av一区二区av| 美女 人体艺术 gogo| 丝袜在线中文字幕| 给我免费播放毛片高清在线观看| 麻豆成人av在线观看| 两个人视频免费观看高清| 亚洲av中文字字幕乱码综合 | av福利片在线| 丝袜在线中文字幕| 给我免费播放毛片高清在线观看| 国产免费男女视频| 国产精品爽爽va在线观看网站 | 香蕉久久夜色| 亚洲av日韩精品久久久久久密| 看免费av毛片| 中文字幕人妻丝袜一区二区| 久久国产精品人妻蜜桃| 一进一出好大好爽视频| 精品日产1卡2卡| 此物有八面人人有两片| 国产色视频综合| 亚洲欧美精品综合一区二区三区| 国产精品自产拍在线观看55亚洲| 看免费av毛片| 天堂√8在线中文| 高潮久久久久久久久久久不卡| 18禁国产床啪视频网站| 丰满人妻熟妇乱又伦精品不卡| 一a级毛片在线观看| 97超级碰碰碰精品色视频在线观看| 国产亚洲av高清不卡| 满18在线观看网站| 91在线观看av| 亚洲无线在线观看| 窝窝影院91人妻| 亚洲国产欧洲综合997久久, | www.自偷自拍.com| 欧美日韩一级在线毛片| 不卡一级毛片| 国产亚洲av嫩草精品影院| av免费在线观看网站| 久久久精品欧美日韩精品| 中文字幕最新亚洲高清| 国产欧美日韩一区二区三| 中文字幕最新亚洲高清| 一进一出好大好爽视频| 制服人妻中文乱码| 国产免费av片在线观看野外av| 久久久久国产一级毛片高清牌| av天堂在线播放| bbb黄色大片| 午夜视频精品福利| 制服诱惑二区| 日韩中文字幕欧美一区二区| 色精品久久人妻99蜜桃| 美女国产高潮福利片在线看| 久久中文字幕一级| 欧美日本视频| 一个人免费在线观看的高清视频| 亚洲精品av麻豆狂野| 亚洲av成人av| 淫妇啪啪啪对白视频| 日韩欧美三级三区| 最近最新中文字幕大全电影3 | 亚洲九九香蕉| 久久精品国产综合久久久| 90打野战视频偷拍视频| 亚洲av中文字字幕乱码综合 | 久久伊人香网站| 午夜a级毛片| 亚洲男人的天堂狠狠| 午夜福利欧美成人| 狂野欧美激情性xxxx| 丁香欧美五月| 国产v大片淫在线免费观看| 国产男靠女视频免费网站| 黄色片一级片一级黄色片| 亚洲全国av大片| 国产国语露脸激情在线看| 一区二区三区激情视频| 国产成人精品久久二区二区91| 亚洲专区字幕在线| 欧美日韩亚洲综合一区二区三区_| 国产亚洲精品av在线| 嫩草影院精品99| 国产熟女xx| 精品国产乱码久久久久久男人| 欧美国产精品va在线观看不卡| 久久亚洲真实| 午夜两性在线视频| 国产av又大| 自线自在国产av| 亚洲成国产人片在线观看| 中文字幕人成人乱码亚洲影| 亚洲第一电影网av| 婷婷六月久久综合丁香| 国产av一区在线观看免费| 久久天堂一区二区三区四区| 一卡2卡三卡四卡精品乱码亚洲| 国产又色又爽无遮挡免费看| 亚洲欧洲精品一区二区精品久久久| 国产精品一区二区三区四区久久 |