• <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
    精品酒店卫生间| 日韩成人伦理影院| 亚洲精品日本国产第一区| 777米奇影视久久| 欧美高清成人免费视频www| 18禁在线播放成人免费| 一区二区三区高清视频在线| 黄色欧美视频在线观看| 日韩欧美国产在线观看| 大香蕉97超碰在线| 综合色丁香网| 天堂√8在线中文| 3wmmmm亚洲av在线观看| av在线蜜桃| a级一级毛片免费在线观看| 成人亚洲精品av一区二区| 日韩亚洲欧美综合| 国产黄片视频在线免费观看| 精品久久久久久久久亚洲| 国产精品1区2区在线观看.| av专区在线播放| 高清毛片免费看| 免费黄频网站在线观看国产| 日韩欧美一区视频在线观看 | 小蜜桃在线观看免费完整版高清| 国产一区二区三区综合在线观看 | 国产成人免费观看mmmm| 久久久久精品久久久久真实原创| av免费在线看不卡| 国产精品国产三级国产av玫瑰| 少妇人妻一区二区三区视频| 中文欧美无线码| 欧美+日韩+精品| 国产色婷婷99| 国产高清三级在线| 大香蕉97超碰在线| av在线播放精品| 免费看a级黄色片| 91aial.com中文字幕在线观看| 欧美精品一区二区大全| 人妻夜夜爽99麻豆av| 一个人看的www免费观看视频| 十八禁网站网址无遮挡 | 最近2019中文字幕mv第一页| 两个人视频免费观看高清| 亚洲18禁久久av| 91aial.com中文字幕在线观看| 久久久久久久久久人人人人人人| 国产在线男女| 内地一区二区视频在线| 大香蕉97超碰在线| 视频中文字幕在线观看| 高清日韩中文字幕在线| 成人毛片a级毛片在线播放| 日日摸夜夜添夜夜添av毛片| 精品亚洲乱码少妇综合久久| 亚洲精品日本国产第一区| 91久久精品国产一区二区成人| 亚州av有码| 免费高清在线观看视频在线观看| 岛国毛片在线播放| 久久久精品免费免费高清| 久久久久久久国产电影| 国产在视频线精品| 亚洲av一区综合| av免费在线看不卡| 街头女战士在线观看网站| 边亲边吃奶的免费视频| 黄色日韩在线| 国产精品麻豆人妻色哟哟久久 | 国产精品.久久久| 精品久久久噜噜| 免费看av在线观看网站| 国产精品久久久久久精品电影小说 | 精品一区二区免费观看| 80岁老熟妇乱子伦牲交| 久久草成人影院| 日韩 亚洲 欧美在线| 亚洲高清免费不卡视频| a级毛片免费高清观看在线播放| 免费播放大片免费观看视频在线观看| 一区二区三区高清视频在线| 国产成人精品福利久久| 在现免费观看毛片| 久久综合国产亚洲精品| 在线观看一区二区三区| 国产永久视频网站| 少妇高潮的动态图| 亚洲三级黄色毛片| 蜜桃久久精品国产亚洲av| 男人爽女人下面视频在线观看| 啦啦啦中文免费视频观看日本| 久久久久久久久久黄片| 日韩欧美精品v在线| 国产成人精品久久久久久| 天天一区二区日本电影三级| 久久这里有精品视频免费| 国产精品蜜桃在线观看| 久99久视频精品免费| 国精品久久久久久国模美| 日本猛色少妇xxxxx猛交久久| 少妇被粗大猛烈的视频| 在线观看av片永久免费下载| 伦理电影大哥的女人| 欧美潮喷喷水| 天堂影院成人在线观看| 日本黄色片子视频| 中文乱码字字幕精品一区二区三区 | 美女xxoo啪啪120秒动态图| 国产成年人精品一区二区| 亚洲av成人精品一二三区| 亚洲精品,欧美精品| 纵有疾风起免费观看全集完整版 | 国产乱人视频| 黄色日韩在线| 久久鲁丝午夜福利片| 美女被艹到高潮喷水动态| 亚洲人成网站高清观看| 高清毛片免费看| 国产精品蜜桃在线观看| 婷婷色av中文字幕| 日本欧美国产在线视频| 免费观看无遮挡的男女| 午夜激情久久久久久久| 免费看美女性在线毛片视频| 超碰av人人做人人爽久久| 18禁在线播放成人免费| 观看美女的网站| 一个人观看的视频www高清免费观看| av国产久精品久网站免费入址| 一级毛片 在线播放| 亚洲成人精品中文字幕电影| 少妇被粗大猛烈的视频| 26uuu在线亚洲综合色| 身体一侧抽搐| 边亲边吃奶的免费视频| 91aial.com中文字幕在线观看| 欧美日韩精品成人综合77777| 日韩大片免费观看网站| 搡女人真爽免费视频火全软件| 男人舔女人下体高潮全视频| 亚洲精品,欧美精品| 国产黄频视频在线观看| 久久久久精品久久久久真实原创| 国产伦精品一区二区三区视频9| 久久久久久伊人网av| 日韩av免费高清视频| 国产亚洲91精品色在线| 日韩强制内射视频| 69人妻影院| 夜夜看夜夜爽夜夜摸| 你懂的网址亚洲精品在线观看| 免费av观看视频| 日韩在线高清观看一区二区三区| 少妇的逼水好多| 人妻夜夜爽99麻豆av| 只有这里有精品99| av网站免费在线观看视频 | 在线天堂最新版资源| 成年版毛片免费区| 小蜜桃在线观看免费完整版高清| 欧美+日韩+精品| 精品国产一区二区三区久久久樱花 | 成人av在线播放网站| av女优亚洲男人天堂| 亚洲欧美中文字幕日韩二区| 国产成人aa在线观看| 国产精品一区www在线观看| 三级国产精品片| 亚洲成人一二三区av| 久久久久久久国产电影| 观看免费一级毛片| 五月玫瑰六月丁香| 街头女战士在线观看网站| 成人国产麻豆网| 又粗又硬又长又爽又黄的视频| or卡值多少钱| 国产av国产精品国产| 亚洲欧美一区二区三区国产| 99热这里只有是精品在线观看| 国产一级毛片在线| 蜜桃亚洲精品一区二区三区| 观看免费一级毛片| 国产成人午夜福利电影在线观看| 久久鲁丝午夜福利片| 伊人久久精品亚洲午夜| 夫妻性生交免费视频一级片| 九九在线视频观看精品| 边亲边吃奶的免费视频| 国产三级在线视频| 亚洲熟女精品中文字幕| av免费在线看不卡| 国产视频首页在线观看| 一个人免费在线观看电影| 大片免费播放器 马上看| 久久久久九九精品影院| 久99久视频精品免费| 如何舔出高潮| 精品久久久噜噜| 久久久久国产网址| 在线天堂最新版资源| 亚洲av国产av综合av卡| 26uuu在线亚洲综合色| 国内揄拍国产精品人妻在线| 亚洲精品日韩av片在线观看| 成人高潮视频无遮挡免费网站| 久久精品久久久久久久性| 女人十人毛片免费观看3o分钟| 99热网站在线观看| 午夜福利高清视频| 99re6热这里在线精品视频| 亚洲不卡免费看| 能在线免费观看的黄片| 亚洲aⅴ乱码一区二区在线播放| 一个人免费在线观看电影| 成人漫画全彩无遮挡| 久久97久久精品| 高清av免费在线| 亚洲av.av天堂| 亚洲精品成人久久久久久| 亚洲经典国产精华液单| 欧美高清成人免费视频www| 国产成人精品一,二区| 午夜久久久久精精品| 国产视频内射| 亚洲国产精品成人久久小说| 国产精品精品国产色婷婷| 舔av片在线| 亚洲精品久久午夜乱码| 日韩欧美精品免费久久| 老司机影院成人| 国产精品久久久久久久久免| 午夜精品国产一区二区电影 | 久久99精品国语久久久| 日本午夜av视频| av免费在线看不卡| 少妇的逼好多水| 超碰av人人做人人爽久久| 日韩av在线大香蕉| 男女视频在线观看网站免费| 婷婷色综合www| 国产视频内射| 精品久久久久久成人av| 又爽又黄a免费视频| 免费观看精品视频网站| 欧美高清成人免费视频www| 亚洲最大成人中文| 日韩欧美精品v在线| 男女国产视频网站| 亚洲精品国产av成人精品| 狂野欧美白嫩少妇大欣赏| 又爽又黄无遮挡网站| 蜜桃久久精品国产亚洲av| 欧美97在线视频| kizo精华| 黄色配什么色好看| 午夜福利在线观看免费完整高清在| 2022亚洲国产成人精品| 少妇丰满av| 麻豆国产97在线/欧美| 国产乱人视频| 黄片无遮挡物在线观看| 精品一区二区免费观看| 国产成人福利小说| 欧美zozozo另类| 亚洲国产最新在线播放| 成人性生交大片免费视频hd| 婷婷色av中文字幕| 丝袜喷水一区| 伦精品一区二区三区| 插逼视频在线观看| 麻豆久久精品国产亚洲av| 91aial.com中文字幕在线观看| 99热网站在线观看| 青春草亚洲视频在线观看| 六月丁香七月| 国产老妇伦熟女老妇高清| 在线免费观看不下载黄p国产| 欧美日韩亚洲高清精品| 久久精品熟女亚洲av麻豆精品 | h日本视频在线播放| av网站免费在线观看视频 | 久久久久久国产a免费观看| 校园人妻丝袜中文字幕| 精品久久久精品久久久| 久久精品夜色国产| 成人欧美大片| 人妻制服诱惑在线中文字幕| 一级毛片黄色毛片免费观看视频| 国产av在哪里看| 在线免费观看的www视频| 毛片一级片免费看久久久久| 亚洲18禁久久av| 国产黄a三级三级三级人| 啦啦啦中文免费视频观看日本| 禁无遮挡网站| 简卡轻食公司| 丝瓜视频免费看黄片| 黄色欧美视频在线观看| 亚洲丝袜综合中文字幕| 国产淫语在线视频| 黄色日韩在线| 久久久午夜欧美精品| 久久久久久伊人网av| 视频中文字幕在线观看| 精品欧美国产一区二区三| 日本三级黄在线观看| 欧美日本视频| 精品国产一区二区三区久久久樱花 | 国产精品不卡视频一区二区| 麻豆久久精品国产亚洲av| 日本午夜av视频| 久久久久久久午夜电影| 欧美xxxx黑人xx丫x性爽| 成人鲁丝片一二三区免费| 精品国内亚洲2022精品成人| 免费大片黄手机在线观看| 午夜福利网站1000一区二区三区| 久久这里只有精品中国| 两个人视频免费观看高清| 中文乱码字字幕精品一区二区三区 | 国产色爽女视频免费观看| 国内精品一区二区在线观看| 国产一区二区三区av在线| 亚洲内射少妇av| 国语对白做爰xxxⅹ性视频网站| 欧美丝袜亚洲另类| 亚洲不卡免费看| 欧美+日韩+精品| 国产黄色免费在线视频| 午夜激情久久久久久久| 日韩欧美精品免费久久| 亚洲国产成人一精品久久久| 中文字幕av在线有码专区| 身体一侧抽搐| av免费在线看不卡| 午夜亚洲福利在线播放| 亚洲精品第二区| 成人美女网站在线观看视频| www.av在线官网国产| 免费看a级黄色片| 国产男人的电影天堂91| 91精品伊人久久大香线蕉| 成人二区视频| 天堂av国产一区二区熟女人妻| 全区人妻精品视频| 美女大奶头视频| 久久综合国产亚洲精品| 国产精品久久久久久精品电影小说 | 国产大屁股一区二区在线视频| 亚洲av二区三区四区| 男人舔女人下体高潮全视频| 久久久国产一区二区| 成人午夜高清在线视频| 欧美3d第一页| 久久亚洲国产成人精品v| freevideosex欧美| 色哟哟·www| 久久精品夜色国产| 亚洲三级黄色毛片| 国产免费福利视频在线观看| 黑人高潮一二区| 能在线免费看毛片的网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一区有黄有色的免费视频 | 日韩,欧美,国产一区二区三区| 人人妻人人澡欧美一区二区| 国产精品三级大全| 亚洲精品国产av成人精品| 久久久精品免费免费高清| 男人舔女人下体高潮全视频| 久久6这里有精品| 嘟嘟电影网在线观看| 亚洲在久久综合| av国产免费在线观看| 亚洲高清免费不卡视频| 亚洲美女搞黄在线观看| 一个人看视频在线观看www免费| 日产精品乱码卡一卡2卡三| 久久久久久国产a免费观看| 狂野欧美白嫩少妇大欣赏| 天堂av国产一区二区熟女人妻| 永久免费av网站大全| 亚洲国产欧美在线一区| 日本一二三区视频观看| 一级毛片aaaaaa免费看小| 国产免费一级a男人的天堂| 丰满人妻一区二区三区视频av| 在现免费观看毛片| 一级毛片aaaaaa免费看小| 在线观看av片永久免费下载| 午夜亚洲福利在线播放| 超碰97精品在线观看| 亚洲国产日韩欧美精品在线观看| 国产永久视频网站| 一级毛片电影观看| 麻豆av噜噜一区二区三区| 少妇高潮的动态图| 日韩精品有码人妻一区| 国产成人a∨麻豆精品| 熟妇人妻不卡中文字幕| 国产精品1区2区在线观看.| 国产 亚洲一区二区三区 | 国产黄片视频在线免费观看| 免费看av在线观看网站| 熟妇人妻久久中文字幕3abv| av在线观看视频网站免费| 天美传媒精品一区二区| 国产亚洲一区二区精品| 91精品一卡2卡3卡4卡| 亚洲精品成人久久久久久| 一二三四中文在线观看免费高清| 真实男女啪啪啪动态图| 亚洲欧美精品自产自拍| 最近最新中文字幕免费大全7| 亚洲综合精品二区| 免费看不卡的av| 国产亚洲5aaaaa淫片| 两个人视频免费观看高清| 亚洲久久久久久中文字幕| 欧美区成人在线视频| 亚洲一级一片aⅴ在线观看| 亚洲精品日本国产第一区| 尾随美女入室| 亚洲精品成人久久久久久| 男女视频在线观看网站免费| 熟妇人妻不卡中文字幕| 久久久久久九九精品二区国产| 床上黄色一级片| 少妇人妻精品综合一区二区| 欧美激情国产日韩精品一区| 久久国产乱子免费精品| 婷婷色综合www| 可以在线观看毛片的网站| 亚洲av国产av综合av卡| 视频中文字幕在线观看| 亚洲精品国产av蜜桃| 欧美最新免费一区二区三区| 18禁在线播放成人免费| 老师上课跳d突然被开到最大视频| 国产熟女欧美一区二区| 九色成人免费人妻av| 最近的中文字幕免费完整| 伦精品一区二区三区| 日韩伦理黄色片| 亚洲婷婷狠狠爱综合网| av在线老鸭窝| 2022亚洲国产成人精品| 啦啦啦中文免费视频观看日本| 高清在线视频一区二区三区| 日韩av在线免费看完整版不卡| a级毛片免费高清观看在线播放| 久久精品夜色国产| 看免费成人av毛片| 最近最新中文字幕免费大全7| 成人性生交大片免费视频hd| 免费看光身美女| 七月丁香在线播放| 日韩精品青青久久久久久| 国产精品av视频在线免费观看| 在线免费观看不下载黄p国产| 精品熟女少妇av免费看| 男女视频在线观看网站免费| 亚洲三级黄色毛片| 日韩欧美精品v在线| 97热精品久久久久久| 国产午夜精品论理片| 日韩大片免费观看网站| 亚洲人与动物交配视频| 天天一区二区日本电影三级| 99久久中文字幕三级久久日本| 少妇人妻精品综合一区二区| 亚洲欧美日韩东京热| 国产亚洲91精品色在线| 干丝袜人妻中文字幕| 亚洲av日韩在线播放| 深爱激情五月婷婷| 非洲黑人性xxxx精品又粗又长| 日本熟妇午夜| 国产av在哪里看| 淫秽高清视频在线观看| 欧美高清性xxxxhd video| 永久网站在线| 午夜亚洲福利在线播放| 一级爰片在线观看| av在线观看视频网站免费| 97热精品久久久久久| 18+在线观看网站| 26uuu在线亚洲综合色| 2021少妇久久久久久久久久久| 精品酒店卫生间| 国产国拍精品亚洲av在线观看| 特大巨黑吊av在线直播| 免费高清在线观看视频在线观看| 午夜老司机福利剧场| av福利片在线观看| 亚洲aⅴ乱码一区二区在线播放| 搞女人的毛片| 国产亚洲最大av| 国产淫语在线视频| 国产欧美另类精品又又久久亚洲欧美| 天堂影院成人在线观看| 26uuu在线亚洲综合色| 极品教师在线视频| 久久99蜜桃精品久久| 亚洲成人久久爱视频| 成年女人看的毛片在线观看| 久久久久九九精品影院| 看十八女毛片水多多多| 少妇裸体淫交视频免费看高清| 成人美女网站在线观看视频| 亚洲综合色惰| 丝瓜视频免费看黄片| 成人国产麻豆网| 中文字幕制服av| 日本午夜av视频| 免费看av在线观看网站| 亚洲精品视频女| 国产不卡一卡二| 久久久久精品性色| 亚洲欧美一区二区三区黑人 | 男女那种视频在线观看| 亚洲国产av新网站| 国产v大片淫在线免费观看| 免费观看av网站的网址| 欧美xxxx黑人xx丫x性爽| 看十八女毛片水多多多| 两个人视频免费观看高清| 可以在线观看毛片的网站| 国产欧美另类精品又又久久亚洲欧美| 国产久久久一区二区三区| av免费观看日本| 干丝袜人妻中文字幕| 亚洲精品亚洲一区二区| 亚洲欧美中文字幕日韩二区| 身体一侧抽搐| 午夜精品在线福利| av.在线天堂| 成人无遮挡网站| 亚洲成人av在线免费| 久久6这里有精品| 国产精品一区二区三区四区免费观看| 听说在线观看完整版免费高清| 国产永久视频网站| 欧美一区二区亚洲| 尾随美女入室| 夫妻午夜视频| 精品久久久久久久人妻蜜臀av| 久久久久久久久久黄片| 简卡轻食公司| 老司机影院毛片| 日本爱情动作片www.在线观看| 又黄又爽又刺激的免费视频.| 只有这里有精品99| 亚洲成人一二三区av| 精品人妻偷拍中文字幕| 乱系列少妇在线播放| 天天躁夜夜躁狠狠久久av| 国产精品一区www在线观看| 国产午夜福利久久久久久| 午夜福利视频1000在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩卡通动漫| 看免费成人av毛片| 免费观看性生交大片5| 日本熟妇午夜| 69人妻影院| 精品少妇黑人巨大在线播放| 青青草视频在线视频观看| 久久亚洲国产成人精品v| 我要看日韩黄色一级片| 国产午夜精品论理片| 我要看日韩黄色一级片| 免费av观看视频| 国产高清国产精品国产三级 | 我的女老师完整版在线观看| 九九爱精品视频在线观看| 欧美日韩在线观看h| 亚洲人与动物交配视频| av网站免费在线观看视频 | 伊人久久国产一区二区| 免费播放大片免费观看视频在线观看| 精品亚洲乱码少妇综合久久| 又大又黄又爽视频免费| 亚洲久久久久久中文字幕| 国内揄拍国产精品人妻在线| 欧美 日韩 精品 国产| 男女下面进入的视频免费午夜| 久久久久久久国产电影| 深夜a级毛片| 极品教师在线视频| 男人舔女人下体高潮全视频| 亚洲人成网站在线播| 精华霜和精华液先用哪个| 人妻系列 视频| 久久99精品国语久久久| 蜜桃亚洲精品一区二区三区| 91av网一区二区| 成人亚洲精品av一区二区| 日韩精品有码人妻一区| 高清在线视频一区二区三区| 精品一区二区免费观看| 亚洲av在线观看美女高潮| 我的女老师完整版在线观看| 极品少妇高潮喷水抽搐| 精品一区二区免费观看| 久久精品国产鲁丝片午夜精品| 在线免费观看的www视频| 99久久精品一区二区三区| 亚洲av电影在线观看一区二区三区 | 丰满少妇做爰视频| 寂寞人妻少妇视频99o| 国内精品宾馆在线| 身体一侧抽搐| 最近的中文字幕免费完整|