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

    45°大傾角傾斜柔性圓柱渦激振動不相關原則實驗驗證

    2017-04-21 00:50:49徐萬海馬燁璇
    振動與沖擊 2017年7期
    關鍵詞:順流橫流渦激

    徐萬海, 馬燁璇, 杜 杰, 羅 浩

    (天津大學 水利工程仿真與安全國家重點實驗室, 天津 300072)

    45°大傾角傾斜柔性圓柱渦激振動不相關原則實驗驗證

    徐萬海, 馬燁璇, 杜 杰, 羅 浩

    (天津大學 水利工程仿真與安全國家重點實驗室, 天津 300072)

    渦激振動是誘發(fā)海洋管道與浮式平臺系泊纜線等細長柔性圓柱結構疲勞損傷的重要因素,已有研究大多關注結構軸向與來流垂直的情況,實際的海洋工程中,圓柱結構軸向與來流并不完全垂直,存在一定傾斜角度。針對這種復雜的情況,傾斜圓柱渦激振動的不相關原則被提出,即假定傾斜圓柱渦激振動與來流速度在結構軸向的垂直方向投影分量引起的垂直圓柱情況等價。然而,較大傾角條件下不相關原則的正確與否仍未知。為了論證大傾角傾斜柔性圓柱渦激振動不相關原則,開展了長徑比為350,傾角為45°的傾斜柔性圓柱渦激振動實驗,采用應變傳感器測量結構應變信息,運用模態(tài)分析法對實驗數(shù)據(jù)進行處理。通過大傾角傾斜圓柱與垂直圓柱的實驗結果對比,發(fā)現(xiàn)在控制模態(tài)轉化區(qū)域,傾斜圓柱更容易被激發(fā)高階模態(tài),同時傾斜圓柱與垂直圓柱在測點處的應變和位移差異較大,并均出現(xiàn)高頻特征。最終得出了較大傾角狀態(tài)下傾斜柔性圓柱渦激振動不相關原則不成立的結論。

    傾斜柔性圓柱;模態(tài)分析法;渦激振動;不相關原則

    圓柱結構置于一定來流速度流場中,其兩側會發(fā)生交替的旋渦脫落,誘發(fā)“渦激振動(Vortex-induced vibrations,VIV)”現(xiàn)象。近年來,渦激振動問題已得到人們廣泛關注[1-11]。

    然而,已有研究更多關注來流與圓柱結構軸向垂直的特殊情形[12-14]。實際工程中,特別是海洋工程中的懸鏈線立管和海洋浮式結構系泊纜線等結構軸向與來流并不垂直,存在著不同的傾斜角度。針對這種復雜的情況,傾斜圓柱渦激振動不相關原則(Independence Principle, IP)被提出(見圖1),即假定來流速度為U,傾斜角度為a的圓柱渦激振動與來流速度為Ucos(a)的垂直圓柱情況等價[15-16]。

    圖1 傾斜圓柱渦激振動不相關原則示意圖Fig.1 Independence principle of inclined cylinder

    IP原則的提出開辟了傾斜圓柱渦激振動研究的新途徑,然而其合理性一直存在爭議。FRANZINI等[17]實驗觀測了傾斜剛性圓柱渦激振動特性,傾角分別為0°、10°、20°、30°、45°,研究發(fā)現(xiàn):傾角<20°時,IP原則正確,傾角>20°,圓柱結構渦激振動幅值隨傾角的增大而變小。ZHAO 等[18]采用Petrov-Galerkin 有限元法數(shù)值求解N-S方程,分析了傾角變化范圍為0°~60°的傾斜固定圓柱渦激振動特性,發(fā)現(xiàn)傾角<30 °時,斯托羅哈數(shù)和平均阻力系數(shù)與IP原則得到的結果一致,但升力系數(shù)不滿足IP原則。JAIN等[19]實驗研究了傾角范圍0°~75°的傾斜剛性圓柱渦激振動,雷諾數(shù)的范圍為500~4 000,實驗結果表明傾角為20°時,IP原則成立,傾角為45°和60°時IP原則不再適用。國內亦開展了傾斜圓柱流激振動問題的部分研究工作。李壽英等[20]采用CFD 軟件CFX5.5 對直圓柱和斜角30°傾斜圓柱繞流問題進行了數(shù)值模擬。許常悅等[21]采用大渦模擬方法數(shù)值研究了斜角為60°的偏斜圓柱跨聲速繞流。杜曉慶等[22]同樣采用了大渦模擬方法,研究了斜置圓柱在展向剪切流作用下的氣動性能。

    近期,傾斜柔性結構渦激振動的研究逐步得以重視,BOURGUET等[23]采用DNS方法數(shù)值模擬了長徑比為50,傾角為60°的傾斜柔性圓柱渦激振動,計算雷諾數(shù)為500,重點論證了IP原則是否適用于柔性結構,結果表明,較低的軸向力(τ=T/ρD2U=13.5,τ為無量綱軸向力,ρ外部流體密度,D圓柱結構外徑,U為來流速度)不相關原則不成立,較高軸向力(τ=124)不相關原則成立。

    雖然傾斜圓柱渦激振動不相關原則已得到了初步研究,但是,至今仍缺乏足夠的研究證明不相關原則的成立條件,特別是針對于大傾角柔性傾斜圓柱結構渦激動研究還存在諸多不足。本文采用模型實驗的方法,觀測大傾角傾斜柔性圓柱結構的渦激振動特性,主要目的是對傾斜圓柱渦激振動的不相關原則正確性進行實驗驗證。

    1 實驗方案

    實驗在天津大學水利工程仿真與安全國家重點實驗室內長137 m、寬7.0 m、深3.3 m的拖曳水池(見圖2(a))中進行。實驗中拖車速度間隔為0.05 m/s,范圍為0.05~1.0 m/s,實驗雷諾數(shù)最大可達16 000,圓柱實驗模型兩端通過萬向節(jié)與鋼架相連,軸向力為450 N,鋼架通過螺栓固定于拖車之上。測量兩種工況,圓柱模型傾角為0° 和45°,實驗模型布置如圖2(b)所示。圓柱模型長為5.6 m,外徑為0.016 m,長徑比350,內芯材質為銅管,外表面包裹硅膠管,具體結構參數(shù)參見表1。圓柱模型被均等分為8份,沿軸線方向7個截面的橫流向和順流向粘貼應變片,采集結構振動信息。

    圖2 實驗示意圖Fig.2 Schematic of the experimental set-up

    參數(shù)數(shù)值長度L/m5.60外徑D/m0.016彎曲剛度EI/Nm217.45軸向力T/N450單位長度結構質量ms(kg·m-1)0.3821質量比4ms/(πρD2)1.90長徑比L/D350

    實驗過程中,采用拖車拖動實驗模型,模擬來流均勻作用。待拖車加速平穩(wěn)后,開始測量采集順流向和橫流向的應變信號,采集時間50 s,采樣頻率100 Hz,滿足奈奎斯特定律。為了驗證應變采集系統(tǒng)的正確性以及兩端萬向節(jié)提供的鉸支邊界條件是否合理,測量了傾斜柔性圓柱模型在靜水中的固有頻率 (見表2),由

    表2 固有頻率

    于僅測量7個截面的應變信息,所以最多能得到前7階固有頻率,對比測量數(shù)值與理論計算式(1)[24]獲得的理論值差異情況。

    (1)

    式中,m為單位長度質量,包括結構質量ms與附加質量ma。對比結果表明固有頻率測量值與理論值吻合的十分理想,進一步證明了實驗設計與信號采集系統(tǒng)的可信度。

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

    應變數(shù)據(jù)采集過程中會受到一些低頻干擾和高頻噪音信號影響,需對測得的應變信號進行濾波處理。實驗中拖車運動產(chǎn)生的干擾信號的頻率不到1.0 Hz,遠低于實驗中柔性圓柱結構渦激振動的響應頻率;另外,圓柱模型渦激振動的最高響應頻率(包括倍頻)<40.0 Hz。在后續(xù)實驗數(shù)據(jù)處理之前,均先采用帶通濾波的方法消除<1.0 Hz的低頻拖車運動干擾信號和>40.0 Hz的高頻噪音信號。

    模態(tài)法是柔性圓柱結構渦激振動實驗數(shù)據(jù)處理常用的一種方法[25-28],根據(jù)不同離散點獲得的應變信息,運用模態(tài)法求得圓柱結構任一點的位移信息。分別對橫流向(Cross-flow,CF)和順流向(In-line,IL)的位移進行模態(tài)分解。為簡化起見,僅以CF位移為例。IL位移信息的確定與CF位移類似,在此不再贅述。

    橫流向位移為

    (2)

    式中,y(z,t)為橫流向的位移,z為軸向坐標,t為時間;wn(t)為權重系數(shù);n為圓柱結構的模態(tài)階次;N為確定結構位移所需的模態(tài)數(shù)目;φn(z)為模態(tài)函數(shù),可表示為

    (3)

    曲率與應變具有如下關系

    (4)

    式中,ε為測量獲得的應變;R結構半徑;L結構長度。圓柱結構的振動位移可采用N個模態(tài)疊加來表示

    ΩW=Θ

    (5)

    其中,

    (6)

    W=[w1(t),w2(t),…,wN(t)]T

    (7)

    (8)

    式中,M為應變片的測點數(shù)目,本文為7。實驗中激發(fā)的橫流向最高模態(tài)為4階,小于應變片的數(shù)目7,采用最小二乘法獲得模態(tài)權重系數(shù)

    W=[ΩTΩ]-1ΩTΘ

    (9)

    將式(9)中的權重系數(shù)結果代入式(2),即可獲得結構位移信息。

    3 不相關原則的實驗驗證

    結構響應幅值是反映渦激振動特性一個重要參數(shù),在確定不相關原則是否成立時,首先對比傾斜圓柱與垂直圓柱的橫流向與順流向最大響應幅值,如圖3和圖4所示。圖中橫軸為約化速度,定義Vr=Ucos(a)/f1D,其中U為拖車的速度,a為傾斜角度,垂直時為0°,傾斜時為45°,f1為結構固有頻率,選取理論公式(1)獲得的計算值。觀察圖3可發(fā)現(xiàn),較低的約化速度,傾斜圓柱與垂直圓柱獲得的橫流向響應幅值比較接近,隨著約化速度增加,傾斜圓柱獲得的橫向位移幅值略小于垂直圓柱。對比圖4中順流向響應幅值結果,發(fā)現(xiàn)較低約化速度時,傾斜圓柱位移幅值較大,在約化速度近似為11.5~16.0時,垂直圓柱順流向位移幅值更大。順流向渦激振動較高模態(tài)會被激發(fā),因此順流向最大均方根位移幅值隨約化速度變化十分復雜。本文開展的實驗中,拖車的最大速度為1.0 m/s,采用Ucos 45°

    圖3 橫流向最大均方根位移幅值隨約化速度變化Fig.3 Max RMS CF displacement versus reduced velocity

    圖4 順流向最大均方根位移幅值隨約化速度變化Fig.4 Max RMS IL displacement versus reduced velocity

    表示垂直其軸向的來流速度,最終導致了在更高約化速度時,只有垂直圓柱結構渦激振動結果,缺少傾斜圓柱結果,圖3和圖4,以及后續(xù)圖5和圖6均會出現(xiàn)類似情況。

    圖5 橫流向控制模態(tài)隨流速變化Fig.5 CF dominate mode number versus reduced velocity

    圖6 順流向控制模態(tài)隨流速變化Fig. 6 IL dominate mode number versus reduced velocity

    控制模態(tài)是描述柔性圓柱結構渦激振動特性的另一個重要指標。圖5和圖6分別針對橫流向與順流向情形,對比傾斜圓柱與垂直圓柱控制模態(tài)差異。

    在圖5中傾斜圓柱激發(fā)最高橫流向渦激振動的控制模態(tài)為3階,垂直圓柱為4階,Vr<7.0時,獲得的控制模態(tài)結果幾乎一致,僅有一個特殊情況,約化速度Vr=2.65時,2階模態(tài)為傾斜圓柱的主要控制模態(tài),而垂直圓柱控制模態(tài)卻為1階,這也印證了傾斜柔性圓柱渦激振動的復雜性。約化速度介于7.0~8.75以及17.7附近時,恰巧是不同控制模態(tài)轉化區(qū)域,相比于垂直圓柱,傾斜圓柱更容易被激發(fā)更高階模態(tài)。圖6中傾斜圓柱順流向渦激振動最高5階模態(tài)被激發(fā),垂直圓柱為6階,與橫流向結果類似,在不同模態(tài)轉化區(qū)域,傾斜圓柱更容易發(fā)生高階的順流向渦激振動。

    分別針對具體的來流速度情況,進行對比分析。選取傾斜圓柱較低來流速度U=0.35 m/s,中等來流速度0.75 m/s和較高來流速度1.0 m/s三個實驗工況。根據(jù)不相關原則,對應的垂直圓柱來流速度分別為0.2 m/s,0.5 m/s和0.7 m/s。圖7和圖8給繪出了測點處應變隨時間變化的圖像。第一列為7個測點的應變信息,第二列為頻譜分析結果??梢园l(fā)現(xiàn)應變結果不滿足不相關原則,在距離零點較近的三個測點,垂直圓柱獲得的順流向與橫流向應變更大,距離零點最遠的三個測點,傾斜圓柱應變更大,中點處差異不明顯。通過頻譜分析結果可知,傾斜圓柱和垂直圓柱渦激振動均會出現(xiàn)高頻現(xiàn)象,但是傾斜圓柱應變振動控制頻率略低于垂直圓柱。

    圖7 測點處橫流向應變時間歷程曲線和頻譜分析結果實線為垂直圓柱,U=0.25 m/s,虛線為傾斜圓柱U=0.35 m/s)Fig.7 CF strain history and its spectra at measured positions (Solid line expressed vertical cylinder, U=0.25 m/s,dash line expressed inclined cylinder, U=0.35 m/s)

    圖8 測點處順流向應變時間歷程曲線和頻譜分析結果(實線為垂直圓柱,U=0.25 m/s,虛線為傾斜圓柱U=0.35 m/s)Fig.8 IL strain history and its spectra at measured positions(Solid line expressed inclined cylinder, U=0.25 m/s,dash line expressed inclined cylinder, U=0.35 m/s)

    圖9和圖10給出了測點處橫流向和順流向位移時間歷程曲線及頻譜分析結果,傾斜圓柱的外部來流速度為0.75 m/s,垂直圓柱為0.5 m/s,由圖9可以得到垂直和傾斜兩種情況,測點處橫流向位移響應幅值和頻率近似相同,較好的滿足不相關原則。然而,圖10中的順流向位移結果對比卻差異十分明顯,在兩端最外側4個測點及最中間測點處,垂直圓柱渦激振動順流向幅值大于傾斜圓柱,其他測點兩種情況位移幅值差距不大。振動控制頻率傾斜圓柱略大于垂直圓柱,亦可以觀測到高頻的成分。

    圖9 測點處橫流向位移時間歷程曲線和頻譜分析結果(實線為垂直圓柱,U=0.5 m/s,虛線為傾斜圓柱U=0.75 m/s)Fig.9 CF displacement history and its spectra at measured positions (Solid line expressed inclined cylinder, U=0.5 m/s,dash line expressed inclined cylinder, U=0.75 m/s)

    圖10 測點處順流向位移時間歷程曲線和頻譜分析結果(實線為垂直圓柱,U=0.5 m/s,虛線為傾斜圓柱U=0.75 m/s)Fig.10 IL displacement history and its spectra at measured positions(Solid line expressed inclined cylinder, U=0.5 m/s, dash line expressed inclined cylinder, U=0.75 m/s)

    最后繪出斜圓柱的外部流速為1.0 m/s,垂直圓柱為0.7 m/s時,橫流向與順流向位移均方根值的空間分布情況,如圖11和圖12所示,垂直圓柱橫流向渦激振動的控制模態(tài)為2階,近似空間對稱,而傾斜圓柱的控制模態(tài)為3階,不存在空間對稱。同時最大響應幅值垂直圓柱的結果略大,再次驗證了圖3和圖5中的結論。垂直圓柱的順流向渦激振動控制模態(tài)為4階,仍然具有空間對稱形式,傾斜圓柱的控制模態(tài)為5階,不具備空間對稱的特征。驗證了圖4和圖6結論的同時,也證明了不相關原則的不合理性。

    圖11 橫流向位移均方根空間分布(實線為垂直圓柱,U=0.7 m/s,虛線為傾斜圓柱U=1.0 m/s)Fig.11 Spatial distribution of CF RMS displacement(Solid line expressed inclined cylinder, U=0.7 m/s,dash line expressed inclined cylinder, U=1.0 m/s)

    圖12 順流向位移均方根空間分布(實線為垂直圓柱,U=0.7 m/s,虛線為傾斜圓柱U=1.0 m/s)Fig.12 Spatial distribution of IL RMS displacement(Solid line expressed inclined cylinder, U=0.7 m/s, dash line expressed inclined cylinder, U=1.0 m/s)

    依據(jù)響應幅值,控制模態(tài),應變和位移時間歷程變化,頻譜分析,均方根位移的空間分布等結果,發(fā)現(xiàn)傾角為45°時,柔性傾斜圓柱渦激振動不相關原則不再適用。按照BOURGUET等軸向力無量綱化的方法(τ=T/ρD2U),本文實驗的無量綱軸向力范圍為1.76×103≤τ≤3.52×104,符合BOURGUET等定義的較大軸向力條件。然而本文的實驗結果與BOURGUET等近期的數(shù)值模擬結果相矛盾,主要原因是本文實驗圓柱模型長徑比更大,為350,遠大于BOURGUET等的50,同時本文實驗雷諾數(shù)最大可達16 000,而BOURGUET等的計算雷諾數(shù)僅為500,長徑比和雷諾數(shù)是影響柔性結構渦激振動的重要參數(shù),兩個參數(shù)不同引起了本文的實驗結果與BOURGUET等的數(shù)值結果出現(xiàn)了不一致。

    4 結 論

    本文針對大傾斜圓柱渦激振動開展了實驗研究,采用了大長徑比和小質量比圓柱實驗模型,實驗雷諾數(shù)最高可達16 000,通過實驗觀測,并與垂直圓柱渦激振動實驗結果對比,發(fā)現(xiàn):

    (1) 較低約化速度時,傾斜圓柱與垂直圓柱獲得的橫流向與順流向渦激振動響應幅值均比較接近,但當約化速度進一步增加,傾斜圓柱獲得的位移幅值與垂直圓柱之間差異變大。

    (2) 在不同控制模態(tài)轉化區(qū)域,傾斜圓柱更容易被激發(fā)高階模態(tài),由于順流向渦激振動本身激發(fā)的模態(tài)較高,傾斜圓柱與垂直圓柱獲得的控制模態(tài)數(shù)值差異更加顯著。

    (3) 傾斜圓柱與垂直圓柱在測點處的應變和位移差異較大,兩種情況均出現(xiàn)了高頻現(xiàn)象。

    (4) 均方根位移的空間分布結果差異十分明顯,垂直圓柱的均方根位移具有近似空間對稱性,而傾斜圓柱的渦激振動更加復雜,對稱性不復存在。

    實驗結果進一步證明了在傾角為45°時,柔性圓柱結構不相關原則不成立。至于不相關原則的成立條件,以及與雷諾數(shù)、質量比、長徑比、傾角等關系仍有待進一步實驗和理論探討。

    [ 1 ] SARPKAYA T. Vortex induced oscillations: a selective review[J]. Journal of Applied Mechanics, 1979,46 (2): 241-258.

    [ 2 ] BEARMAN P W. Vortex shedding from oscillating bluff bodies[J]. Annual Review of Fluid Mechanics, 1984, 16: 195-222.

    [ 3 ] VANDIVER J K. Dimensionless parameters important to the prediction of vortex induced vibration of long, flexible cylinders in ocean currents[J]. Journal of Fluids and Structures, 1993, 7(5):423-455.

    [ 4 ] SARPKAYA T. A critical review of the intrinsic nature of vortex-induced vibrations[J]. Journal of Fluids and Structures, 2004, 19(4): 389-447.

    [ 5 ] FACCHINETTI M L, LANGRE E D, BIOLLEY F. Coupling of structure and wake oscillators in vortex-induced vibrations[J]. Journal of Fluids and Structures, 2004,19 (2): 123-140.

    [ 6 ] GABBAI R D, BENAROYA H. An overview of modeling and experiments of vortex-induced vibration of circular cylinders[J]. Journal of Sound and Vibration, 2005, 282: 575-616.

    [ 7 ] VIOLETTE R, LANGRE E D, SZYDLOWSKI J. Computation of vortex-induced vibrations of long structures using a wake oscillator model: comparison with DNS and experiments[J]. Computers and Structures, 2007, 85: 1134-1141.

    [ 8 ] WILLIAMSON C H K, GOVARDHAN R. A brief review of recent results in vortex-induced vibrations[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2008, 96(6/7): 713-735.

    [ 9 ] 宋芳,林黎明,林國燦.圓柱渦激振動的結構-尾流振子耦合模型研究[J]. 力學學報,2010,42(3):357-365. SONG Fang, LIN Liming, LIN Guocan. The study of vortex-induced vibrations by computation using coupling model of structure and wake oscillator[J]. Chinese Journal of Theoretical and Applied Mechanics, 2010, 42(3):357-365.

    [10] WU Xiaodong, GE Fei , HONG Youshi. A review of recent studies on vortex-induced vibrations of long slender cylinders[J]. Journal of Fluids and Structures, 2012, 28(1):292-308.

    [11] XUE Hongxiang, TANG Wenyong, QU Xue. Prediction and analysis of fatigue damage due to cross-flow and in-line VIV for marine risers in non-uniform current[J]. Ocean Engineering,2014,83(2):52-62.

    [12] 黃維平,劉娟,唐世振.考慮流固耦合的大柔性圓柱體渦激振動非線性時域模型[J]. 振動與沖擊,2012,31 (9):140-143. HUANG Weiping, LIU Juan, TANG Shizhen. Nonlinear model of vortex induced vibration of flexible cylinder in consideration of fluid-structure interaction[J].Journal of Vibration and Shock,2012,31(9): 140-143.

    [13] 唐友剛,樊娟娟,張杰,等.高雷諾數(shù)下圓柱順流向和橫向渦激振動分析[J].振動與沖擊,2013,32(13):88-92. TANG Yougang, FAN Juanjuan, ZHANG Jie, et al. In line and transverse vortex-induced vibration analysis for a circular cylinder under high Reynolds number [J]. Journal of Vibration and Shock,2013,32(13) : 88-92.

    [14] 高云,任鐵,付世曉,等.柔性立管渦激振動響應特性試驗研究[J]. 振動與沖擊,2015,34(17):6-11. GAO Yun, REN Tie, FU Shixiao, et al. Tests for response characteristics of VIV of a flexible riser [J].Journal of Vibration and Shock, 2015, 34(17): 6-11.

    [15] HANSON A R. Vortex-shedding from yawed cylinders[J]. AIAA Journal, 1966, 4(4):738-740.

    [16] VAN ATTA C W. Experiments on vortex shedding from yawed circular cylinders [J]. AIAA Journal, 1968,6(5):931-933.

    [17] FRANZINI G R, GON CALVES R T, MENEGHINI J R, et al. One and two degrees-of-freedom vortex-induced vibration experiments with yawed cylinders[J].Journal of Fluids and Structures, 2013, 42(4):401-420.

    [18] ZHAO M, CHENG L, ZHOU T. Direct numerical simulation of three-dimensional flow past a yawed circular cylinder of infinite length[J]. Journal of Fluids and Structures, 2009, 25(5):831-847.

    [19] JAIN A, MODARRES-SADEGHI Y. Vortex-induced vibrations of a flexibly mounted inclined cylinder [J]. Journal of Fluids and Structures, 2013,43(6): 28-40.

    [20] 李壽英,顧明. 斜、直圓柱繞流的CFD 模擬[J]. 空氣動力學學報,2005,23(2):222-227. LI Shouying, GU Ming. Numerical simulation for flow around perpendicular and oblique circular cylinders[J]. Acta Aerodynamic Sinica,2005,23(2):222-227.

    [21] 許常悅,孫智,王從磊.偏斜來流對圓柱跨聲速繞流的影響[J].應用數(shù)學和力學,2014,35(10):1135-1142. XU Changyue, SUN Zhi, WANG Conglei. Effects of tilted free-stream on the transonic flow past a circular cylinder[J]. Applied Mathematics and Mechanics, 2014, 35(10):1135-1142.

    [22] 杜曉慶,費陳杰,況中華,等. 展向剪切流作用下斜置圓柱氣動特性研究[J].振動與沖擊.2014, 33(21):31-37. DU Xiaoqing,F(xiàn)EI Chenjie,KUANG Zhonghua,et al. Aerodynamic characteristics of an inclined circular cylinder in a span-wise linear shear flow[J]. Journal of Vibration and Shock, 2014,33(21):31-37.

    [23] BOURGUETA R, KARNIADAKIS G E, TRIANTAFYLLOU M S. On the validity of the independence principle applied to the vortex-induced vibrations of a flexible cylinder inclined at 60°[J]. Journal of Fluids and Structures, 2015, 53:58-69.

    [24] SONG Ji’ning, LU Lin, TENG Bin, et al. Laboratory tests of vortex-induced vibrations of a long flexible riser pipe subjected to uniform flow [J]. Ocean Engineering, 2011, 38(11):1308-1322.

    [25] LIE H, KAASEN K E. Modal analysis of measurements from a large-scale VIV model test of a riser in linearly sheared flow[J]. Journal of Fluids and Structures, 2006, 22(4): 557-575.

    [26] CHAPLIN J R, BEARMAN P W, HUERA HUARTE F J, et al. Laboratory measurements of vortex-induced vibrations of a vertical tension riser in a stepped current[J]. Journal of Fluids and Structures, 2005, 21(1): 3-24.

    [27] TRIM A D, BRAATEN H, LIE H, et al. Experimental investigation of vortex-induced vibration of long marine risers[J]. Journal of Fluids and Structures, 2005, 21(3): 335-361.

    [28] 張建僑,宋吉寧,呂林,等.柔性立管渦激振動實驗的數(shù)據(jù)分析[J]. 中國海洋平臺,2009,24(4):26-32. ZHANG Jianqiao, SONG Ji’ning, Lü Lin, et al. Analysis on the vortex-induced vibrations experimental results of flexible risers[J]. China Offshore Platform, 2009,24(4):26-32.

    Test verification for independence principle applied in vortex-inducedvibrations of a flexible cylinder inclined at 45°

    XU Wanhai, MA Yexuan, DU Jie, LUO Hao

    (State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China)

    Vortex-induced vibrations (VIV) can cause fatigue damage of marine pipeline and mooring lines. The current VIV studies focus on a canonical problem that a circular cylinder freely to oscillates within a flow perpendicular to the body axis. In practical applications, cylinder structures are often inclined with respect to the direction of the oncoming flow. Aiming at this complex situation, the independence principle (IP) that the VIV of an inclined circular cylinder is assumed to be equivalent to the VIV of a vertical circular cylinder excited by the projection component of oncoming flow velocity in the direction perpendicnlar to the cylinder’s axis is proposed. However, more attentions need to be paid to the validity of the independence principle. Here, laboratory tests were conducted for vortex-induced vibrations (VIV) of an inclined flexible cylinder with a ratio of length to diameter of 350 and an incline angle of 45° in order to check the validity of the independence principle (IP). Strain gages were adopted to measure the cylinder’s dynamic response, a modal analysis approach was used to analyze the measured data. The test results of a vertical cylinder and the inclined cylinder were compared. It was shown that the higher order modes are easy to be excited for the inclined cylinder in the controlled mode transition region; the strain and displacement at measured points of both cylinders are obviously different and the higher order modes’ features of VIV responses are observed. All the results indicated that the use of IP may not be correct to describe the VIV of the flexible cylinder inclined at 45°.

    inclined flexible cylinder; modal analysis; vortex-induced vibrations (VIV); independence principle (IP)

    國家自然科學基金(51379144;51479135;51679167);國家自然科學基金創(chuàng)新群體項目科學基金(51621092);國家重點基礎研究計劃-973計劃(2014CB046801);天津市青年科學基金(15JCQNJC07700)

    2015-11-06 修改稿收到日期: 2016-02-03

    徐萬海 男,副教授,碩士生導師,1981年生 E-mail:xuwanhai@tju.edu.cn

    TV312

    A

    10.13465/j.cnki.jvs.2017.07.027

    猜你喜歡
    順流橫流渦激
    不同間距比下串聯(lián)圓柱渦激振動數(shù)值模擬研究
    橫流熱源塔換熱性能研究
    煤氣與熱力(2021年3期)2021-06-09 06:16:20
    渦激振動發(fā)電裝置及其關鍵技術
    舟山老塘山港區(qū)船舶引航與節(jié)能環(huán)保探討
    中國水運(2020年6期)2020-11-17 09:41:54
    順流與逆流
    人物(2019年9期)2019-02-13 03:00:50
    盤球立管結構抑制渦激振動的數(shù)值分析方法研究
    電子制作(2018年14期)2018-08-21 01:38:42
    基于橫流風扇技術的直升機反扭驗證
    長期股權投資逆流交易和順流交易淺議
    柔性圓管在渦激振動下的模態(tài)響應分析
    脊下橫流對PEMFC性能影響的數(shù)值分析
    電源技術(2015年11期)2015-08-22 08:50:26
    一二三四在线观看免费中文在| 国产精品久久电影中文字幕| 在线观看免费视频日本深夜| 国产精品女同一区二区软件 | 久久精品aⅴ一区二区三区四区| 亚洲片人在线观看| 97超视频在线观看视频| 在线观看午夜福利视频| 亚洲avbb在线观看| 亚洲va日本ⅴa欧美va伊人久久| 波多野结衣高清无吗| 亚洲中文av在线| 男人的好看免费观看在线视频| 久久草成人影院| 99久久99久久久精品蜜桃| 全区人妻精品视频| 欧美午夜高清在线| 久久九九热精品免费| 欧美另类亚洲清纯唯美| 婷婷精品国产亚洲av| 午夜福利18| 久久久久久大精品| 午夜激情欧美在线| 91九色精品人成在线观看| 日韩人妻高清精品专区| 欧美一区二区国产精品久久精品| 非洲黑人性xxxx精品又粗又长| 人妻丰满熟妇av一区二区三区| 国产精品98久久久久久宅男小说| 精品国产美女av久久久久小说| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲欧美98| 色综合婷婷激情| 伊人久久大香线蕉亚洲五| 老汉色av国产亚洲站长工具| 国产精品美女特级片免费视频播放器 | 午夜日韩欧美国产| www.自偷自拍.com| 免费看美女性在线毛片视频| 欧美黄色淫秽网站| 欧美乱色亚洲激情| 国产不卡一卡二| 18禁美女被吸乳视频| 久久这里只有精品19| 看片在线看免费视频| 国产高潮美女av| 亚洲,欧美精品.| 九九久久精品国产亚洲av麻豆 | 国产精品永久免费网站| 亚洲熟女毛片儿| 日本 av在线| 亚洲欧美日韩卡通动漫| 欧美成人免费av一区二区三区| 国产亚洲av高清不卡| 日韩成人在线观看一区二区三区| 国产伦人伦偷精品视频| 欧美中文综合在线视频| 国产精品久久久久久久电影 | 久久久久久大精品| 欧美中文日本在线观看视频| 亚洲专区字幕在线| 男女下面进入的视频免费午夜| а√天堂www在线а√下载| 又爽又黄无遮挡网站| 在线十欧美十亚洲十日本专区| 欧美三级亚洲精品| 国产爱豆传媒在线观看| 亚洲国产精品sss在线观看| 免费av毛片视频| 免费看美女性在线毛片视频| 中国美女看黄片| 国产成人av教育| 99久久综合精品五月天人人| 真实男女啪啪啪动态图| 欧美另类亚洲清纯唯美| 变态另类丝袜制服| 美女黄网站色视频| 亚洲国产精品sss在线观看| 亚洲精品色激情综合| 91av网站免费观看| 国产精品 国内视频| 美女 人体艺术 gogo| 国内久久婷婷六月综合欲色啪| 99国产精品99久久久久| 全区人妻精品视频| 99re在线观看精品视频| 久久精品国产亚洲av香蕉五月| 九色成人免费人妻av| 伊人久久大香线蕉亚洲五| 99热这里只有是精品50| 99久久99久久久精品蜜桃| 国产视频一区二区在线看| 最好的美女福利视频网| 国产精品免费一区二区三区在线| 俄罗斯特黄特色一大片| 一本一本综合久久| 99热这里只有精品一区 | 亚洲欧美精品综合一区二区三区| 亚洲av成人不卡在线观看播放网| 国产精品久久久久久精品电影| 听说在线观看完整版免费高清| 午夜福利视频1000在线观看| 久久精品国产99精品国产亚洲性色| 欧美一级a爱片免费观看看| 亚洲色图av天堂| 女同久久另类99精品国产91| 午夜久久久久精精品| 日韩欧美国产在线观看| 淫妇啪啪啪对白视频| 日韩欧美在线二视频| 欧美黄色淫秽网站| 日日夜夜操网爽| 搞女人的毛片| 麻豆成人av在线观看| 看片在线看免费视频| 午夜激情欧美在线| 国产 一区 欧美 日韩| 少妇熟女aⅴ在线视频| www.自偷自拍.com| 丰满人妻一区二区三区视频av | 亚洲欧美日韩无卡精品| 99在线人妻在线中文字幕| 18美女黄网站色大片免费观看| 国产成人欧美在线观看| 一级毛片精品| 欧美不卡视频在线免费观看| 成人高潮视频无遮挡免费网站| 一二三四社区在线视频社区8| 亚洲国产欧美网| 三级男女做爰猛烈吃奶摸视频| 国语自产精品视频在线第100页| 亚洲熟女毛片儿| 国产黄片美女视频| 国内精品美女久久久久久| 1000部很黄的大片| 亚洲 欧美一区二区三区| h日本视频在线播放| 男女午夜视频在线观看| 九色成人免费人妻av| 亚洲五月天丁香| 亚洲一区高清亚洲精品| 精华霜和精华液先用哪个| 国产精品 国内视频| 欧美色欧美亚洲另类二区| 国内精品久久久久精免费| 黄色日韩在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产午夜精品久久久久久| 啪啪无遮挡十八禁网站| 欧美中文日本在线观看视频| 国产伦精品一区二区三区四那| 九九在线视频观看精品| 午夜a级毛片| 搡老妇女老女人老熟妇| 国产成人aa在线观看| 亚洲欧洲精品一区二区精品久久久| 又黄又粗又硬又大视频| 最新美女视频免费是黄的| 国产激情欧美一区二区| 国语自产精品视频在线第100页| 亚洲美女黄片视频| 色综合站精品国产| 日韩高清综合在线| 亚洲av电影不卡..在线观看| 色老头精品视频在线观看| 亚洲美女黄片视频| 亚洲中文av在线| 精品久久久久久久末码| 亚洲国产色片| 欧美大码av| 亚洲欧美日韩东京热| 欧美+亚洲+日韩+国产| 九色成人免费人妻av| 在线观看免费午夜福利视频| 白带黄色成豆腐渣| 国产亚洲欧美在线一区二区| 色av中文字幕| 国内少妇人妻偷人精品xxx网站 | 亚洲av中文字字幕乱码综合| 在线永久观看黄色视频| 久99久视频精品免费| 搞女人的毛片| 国产午夜精品久久久久久| 老司机午夜福利在线观看视频| 一二三四在线观看免费中文在| 99在线人妻在线中文字幕| 国产成人欧美在线观看| 国产一级毛片七仙女欲春2| 国产一区二区三区在线臀色熟女| 真实男女啪啪啪动态图| 精华霜和精华液先用哪个| 国产亚洲av嫩草精品影院| 欧美日本视频| 88av欧美| 首页视频小说图片口味搜索| 999精品在线视频| 免费大片18禁| 可以在线观看毛片的网站| 久久久国产欧美日韩av| 麻豆国产av国片精品| 日韩国内少妇激情av| 久久中文看片网| 日本撒尿小便嘘嘘汇集6| 国产精品国产高清国产av| 国产蜜桃级精品一区二区三区| 九色成人免费人妻av| 日韩人妻高清精品专区| 99久久精品热视频| 亚洲精品中文字幕一二三四区| 久久婷婷人人爽人人干人人爱| 99国产精品一区二区蜜桃av| 最近最新免费中文字幕在线| 99久国产av精品| 看免费av毛片| 国产精品一区二区精品视频观看| 午夜福利18| 天天添夜夜摸| 中文字幕精品亚洲无线码一区| 国产探花在线观看一区二区| 999久久久国产精品视频| 国产一区二区三区在线臀色熟女| 精品一区二区三区视频在线 | 亚洲精品在线美女| 日韩av在线大香蕉| 91九色精品人成在线观看| 黄色女人牲交| 成年女人毛片免费观看观看9| 国产一级毛片七仙女欲春2| 日本免费a在线| 国产精品一区二区免费欧美| avwww免费| 国产成人一区二区三区免费视频网站| 变态另类成人亚洲欧美熟女| 欧美日韩乱码在线| 亚洲 欧美一区二区三区| 午夜a级毛片| 在线免费观看不下载黄p国产 | 可以在线观看的亚洲视频| 亚洲午夜精品一区,二区,三区| 九色成人免费人妻av| 国产精品久久电影中文字幕| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 757午夜福利合集在线观看| 亚洲精品美女久久av网站| 精品99又大又爽又粗少妇毛片 | 色综合欧美亚洲国产小说| 悠悠久久av| 欧美黑人巨大hd| 亚洲av成人一区二区三| 哪里可以看免费的av片| 制服人妻中文乱码| 热99在线观看视频| 国产v大片淫在线免费观看| 国产精品综合久久久久久久免费| 99riav亚洲国产免费| 一个人看视频在线观看www免费 | 99久久国产精品久久久| 中文亚洲av片在线观看爽| 日本与韩国留学比较| 国产免费男女视频| 麻豆国产97在线/欧美| 成人18禁在线播放| 一级a爱片免费观看的视频| 青草久久国产| 一二三四在线观看免费中文在| 最近视频中文字幕2019在线8| 国产在线精品亚洲第一网站| 黄色女人牲交| 久99久视频精品免费| 变态另类成人亚洲欧美熟女| 色播亚洲综合网| 91九色精品人成在线观看| 五月玫瑰六月丁香| 欧美极品一区二区三区四区| 2021天堂中文幕一二区在线观| 国产欧美日韩一区二区精品| 亚洲国产高清在线一区二区三| 亚洲七黄色美女视频| 一级毛片高清免费大全| 久久精品影院6| 精品久久久久久久久久久久久| 国产精品日韩av在线免费观看| 亚洲黑人精品在线| 天堂影院成人在线观看| 亚洲精品在线美女| 国产欧美日韩一区二区三| 搡老岳熟女国产| a级毛片a级免费在线| 久久久精品欧美日韩精品| 丰满人妻熟妇乱又伦精品不卡| 又粗又爽又猛毛片免费看| 色播亚洲综合网| 两个人看的免费小视频| 亚洲专区中文字幕在线| 国产成人精品无人区| 国内久久婷婷六月综合欲色啪| 亚洲欧美精品综合一区二区三区| 久久久久亚洲av毛片大全| 国产黄色小视频在线观看| 亚洲欧洲精品一区二区精品久久久| 午夜视频精品福利| 日本 欧美在线| netflix在线观看网站| 国产免费av片在线观看野外av| 亚洲欧美日韩高清专用| 搡老熟女国产l中国老女人| 久久国产乱子伦精品免费另类| 国产免费av片在线观看野外av| 色噜噜av男人的天堂激情| 久久久久亚洲av毛片大全| 成人国产一区最新在线观看| 最近视频中文字幕2019在线8| x7x7x7水蜜桃| 国产亚洲精品一区二区www| 啦啦啦免费观看视频1| 久久精品影院6| 亚洲精品中文字幕一二三四区| 色视频www国产| 免费看日本二区| 免费在线观看视频国产中文字幕亚洲| 精品国产亚洲在线| 香蕉久久夜色| 欧美日本亚洲视频在线播放| 亚洲九九香蕉| 欧美成人免费av一区二区三区| av福利片在线观看| 十八禁网站免费在线| 欧美丝袜亚洲另类 | 99国产精品一区二区蜜桃av| 国产精品av视频在线免费观看| 一本综合久久免费| 成年版毛片免费区| 国产精品国产高清国产av| 国内少妇人妻偷人精品xxx网站 | 一级作爱视频免费观看| 国产美女午夜福利| 成人18禁在线播放| 91麻豆av在线| 激情在线观看视频在线高清| 精品一区二区三区av网在线观看| 在线国产一区二区在线| h日本视频在线播放| 免费无遮挡裸体视频| 黑人巨大精品欧美一区二区mp4| 女同久久另类99精品国产91| 国产极品精品免费视频能看的| av在线蜜桃| 国产精品综合久久久久久久免费| 女人高潮潮喷娇喘18禁视频| 国产一区二区在线观看日韩 | 日本黄色片子视频| 欧美另类亚洲清纯唯美| 欧美中文综合在线视频| 国产高清三级在线| 一区二区三区国产精品乱码| 亚洲欧美日韩东京热| 成人特级黄色片久久久久久久| 欧美中文综合在线视频| 亚洲av熟女| 波多野结衣高清作品| 一二三四在线观看免费中文在| 亚洲电影在线观看av| 亚洲乱码一区二区免费版| 亚洲,欧美精品.| 最近最新中文字幕大全免费视频| 免费在线观看成人毛片| 两个人看的免费小视频| 欧美不卡视频在线免费观看| 一区二区三区高清视频在线| 免费av毛片视频| a级毛片a级免费在线| 黄频高清免费视频| 亚洲国产欧洲综合997久久,| 亚洲成人久久性| 午夜两性在线视频| 一进一出抽搐gif免费好疼| 国产av在哪里看| 少妇丰满av| 国产精品一及| 日本免费a在线| 亚洲无线在线观看| 一进一出抽搐gif免费好疼| 小说图片视频综合网站| 美女大奶头视频| 国产精品女同一区二区软件 | 欧美日韩中文字幕国产精品一区二区三区| 精品欧美国产一区二区三| 国产又黄又爽又无遮挡在线| 国产免费男女视频| 看免费av毛片| 久久这里只有精品19| 成人18禁在线播放| 国产精品一区二区精品视频观看| 日韩成人在线观看一区二区三区| 久久精品综合一区二区三区| 18禁黄网站禁片午夜丰满| 人人妻人人看人人澡| 色在线成人网| 亚洲精华国产精华精| 久久伊人香网站| 中文亚洲av片在线观看爽| 日本精品一区二区三区蜜桃| 国产成人aa在线观看| 最近最新中文字幕大全电影3| 国产毛片a区久久久久| 亚洲自拍偷在线| 一个人免费在线观看电影 | 国产精品电影一区二区三区| 国产高潮美女av| 丁香欧美五月| 亚洲18禁久久av| 久久久久国内视频| 日本三级黄在线观看| 国产探花在线观看一区二区| 亚洲av五月六月丁香网| 嫩草影视91久久| 欧美av亚洲av综合av国产av| 日本免费a在线| 99精品久久久久人妻精品| 精品午夜福利视频在线观看一区| 成熟少妇高潮喷水视频| 一个人看视频在线观看www免费 | 国产爱豆传媒在线观看| 一个人观看的视频www高清免费观看 | 中亚洲国语对白在线视频| 日韩欧美在线乱码| 99在线人妻在线中文字幕| 日韩欧美国产在线观看| 精华霜和精华液先用哪个| 一进一出抽搐动态| 久久精品91蜜桃| 观看美女的网站| 久久热在线av| 搡老岳熟女国产| 很黄的视频免费| 欧美日韩精品网址| 国产精品自产拍在线观看55亚洲| 蜜桃久久精品国产亚洲av| 女人高潮潮喷娇喘18禁视频| tocl精华| 亚洲专区国产一区二区| 色综合站精品国产| 日本在线视频免费播放| 啦啦啦免费观看视频1| 午夜精品久久久久久毛片777| 日本熟妇午夜| 美女扒开内裤让男人捅视频| 国产乱人视频| 色精品久久人妻99蜜桃| 真实男女啪啪啪动态图| 国产伦精品一区二区三区四那| 嫩草影视91久久| 又大又爽又粗| 在线观看舔阴道视频| 亚洲熟女毛片儿| 好男人电影高清在线观看| 亚洲欧美一区二区三区黑人| 成人亚洲精品av一区二区| 国产午夜福利久久久久久| 亚洲精品美女久久av网站| svipshipincom国产片| 97碰自拍视频| 亚洲五月婷婷丁香| 欧美乱色亚洲激情| 亚洲无线在线观看| 日本一本二区三区精品| 18禁黄网站禁片免费观看直播| 欧美日韩福利视频一区二区| 99视频精品全部免费 在线 | 国产v大片淫在线免费观看| 此物有八面人人有两片| 久9热在线精品视频| 在线a可以看的网站| 18禁美女被吸乳视频| 麻豆一二三区av精品| 久久中文看片网| 国产成人欧美在线观看| 18禁黄网站禁片免费观看直播| 久久伊人香网站| www.999成人在线观看| netflix在线观看网站| 国产成人精品无人区| 亚洲国产高清在线一区二区三| 熟女人妻精品中文字幕| 欧美在线黄色| 久久久精品大字幕| 欧美性猛交╳xxx乱大交人| 欧美不卡视频在线免费观看| 欧美色视频一区免费| 国产av在哪里看| 亚洲国产中文字幕在线视频| www.999成人在线观看| 搡老熟女国产l中国老女人| 免费搜索国产男女视频| 国内精品美女久久久久久| 国产男靠女视频免费网站| 久久久久性生活片| 男女下面进入的视频免费午夜| 久久久久久久久免费视频了| 最近视频中文字幕2019在线8| 免费无遮挡裸体视频| 麻豆国产av国片精品| 免费看a级黄色片| www.精华液| 长腿黑丝高跟| 又黄又粗又硬又大视频| 丝袜人妻中文字幕| 大型黄色视频在线免费观看| 人妻久久中文字幕网| 99热只有精品国产| 国产综合懂色| 丰满的人妻完整版| 十八禁人妻一区二区| 婷婷丁香在线五月| 大型黄色视频在线免费观看| 午夜福利视频1000在线观看| 91九色精品人成在线观看| 亚洲国产欧洲综合997久久,| 高清在线国产一区| av欧美777| 欧美黑人巨大hd| 国产视频内射| 国产一区二区三区视频了| 午夜亚洲福利在线播放| 国产精品久久电影中文字幕| 看黄色毛片网站| 日本熟妇午夜| av中文乱码字幕在线| 精品欧美国产一区二区三| 欧美日韩精品网址| 在线观看美女被高潮喷水网站 | 一本一本综合久久| 国产精品1区2区在线观看.| 国产私拍福利视频在线观看| 免费人成视频x8x8入口观看| 超碰成人久久| 久久久久久大精品| 麻豆一二三区av精品| 噜噜噜噜噜久久久久久91| 亚洲无线观看免费| 国产精品精品国产色婷婷| 亚洲成av人片在线播放无| 午夜免费观看网址| 亚洲国产精品成人综合色| 免费无遮挡裸体视频| 国产精品久久久久久人妻精品电影| 成年版毛片免费区| 桃红色精品国产亚洲av| 免费在线观看影片大全网站| 国产成人一区二区三区免费视频网站| 久久欧美精品欧美久久欧美| 午夜激情福利司机影院| 日韩人妻高清精品专区| 日本黄色片子视频| 18禁黄网站禁片免费观看直播| 国产久久久一区二区三区| 免费电影在线观看免费观看| 国产精品一区二区免费欧美| 长腿黑丝高跟| 国产熟女xx| cao死你这个sao货| 亚洲欧美精品综合一区二区三区| 丰满人妻一区二区三区视频av | 亚洲色图 男人天堂 中文字幕| 男人舔奶头视频| 首页视频小说图片口味搜索| 男人舔女人下体高潮全视频| 国产高清videossex| 九色成人免费人妻av| 一本精品99久久精品77| 97超视频在线观看视频| 怎么达到女性高潮| 亚洲无线在线观看| 国产免费男女视频| 国产精品 国内视频| 成年女人看的毛片在线观看| 精品熟女少妇八av免费久了| 嫩草影视91久久| 亚洲精品456在线播放app | 国模一区二区三区四区视频 | 国产淫片久久久久久久久 | 久久久久九九精品影院| 久99久视频精品免费| 国产99白浆流出| 在线观看日韩欧美| 丁香六月欧美| www日本黄色视频网| 校园春色视频在线观看| 午夜a级毛片| 高清在线国产一区| 亚洲一区二区三区不卡视频| 久久久国产精品麻豆| 91九色精品人成在线观看| 久久久久久久精品吃奶| 成人特级黄色片久久久久久久| 欧美色视频一区免费| 免费看a级黄色片| 校园春色视频在线观看| 日韩欧美三级三区| 99久久精品一区二区三区| 一个人免费在线观看电影 | 亚洲熟妇熟女久久| 国产伦在线观看视频一区| 99久久无色码亚洲精品果冻| 国产99白浆流出| x7x7x7水蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 久久国产精品人妻蜜桃| 长腿黑丝高跟| 亚洲最大成人中文| 在线观看午夜福利视频| 国内少妇人妻偷人精品xxx网站 | 少妇人妻一区二区三区视频| 日韩欧美国产在线观看|