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

    串列橋墩繞流對行進船舶受力的影響研究

    2018-10-11 12:27:00牛國杰劉曉平葉玉康李安斌
    交通科學與工程 2018年3期
    關(guān)鍵詞:船尾船體橋墩

    牛國杰, 劉曉平, 葉玉康, 李安斌

    (長沙理工大學 水利工程學院,湖南 長沙 410114)

    船舶駛經(jīng)橋墩時,橋墩周圍出現(xiàn)的紊動流場而對船舶水動力造成干擾,引起船舶產(chǎn)生橫漂、旋轉(zhuǎn)運動,誘發(fā)船舶撞橋海損事故。如何保障橋區(qū)船舶通航安全已成為該領(lǐng)域?qū)W者們研究的重點。一些學者[1-3]研究了橋墩周圍的紊動流場,給出橋墩與航道邊線安全距離,并用于指導橋區(qū)航線規(guī)劃。然而,在船舶與橋墩交匯運動過程中,船舶與橋墩兩者間存在著水動力的相互作用,僅僅考慮橋墩周圍的紊流場分布難以體現(xiàn)船舶運動狀態(tài)以及船舶與橋墩交匯過程的安全性。Kijima[4]應用細長體理論,研究了2艘同型船舶在橋墩旁追越過程中兩船間水動力的相互作用規(guī)律。徐言民[5]等人視橋墩為絕對靜止的船舶,采用船間水動力干擾通用模型計算船舶上、下行經(jīng)過橋墩時的船舶與橋墩之間相互干擾的水動力。張晨曦[6-7]等人基于CFD軟件中的FLUENT軟件,模擬在橋墩鄰域直線航行船舶的粘性流場,分析了作用于船體上的橫向力和轉(zhuǎn)艏力矩歷時變化規(guī)律。乾東岳[8-10]等人采用物理模型試驗,對船舶經(jīng)過橋墩過程的受力情況進行了測量。目前,這些研究集中于單橋墩和船舶沿直線經(jīng)過橋墩情況,但無法真實地反映船、橋交匯過程中船舶的航行姿態(tài)、水動力狀態(tài)。因此,作者擬在物理模型試驗的基礎(chǔ)上,考慮船舶與橋墩繞流的耦合作用,對船舶與橋墩間的水動力相互作用進行研究。

    1 物理模型試驗

    1.1 試驗概況

    本物理模型試驗在長沙理工大學水利實驗中心的概化水槽(35 m×3.5 m×1.2 m)中進行,選用1∶50的正態(tài)比尺,原型橋墩直徑為5.0 m,船舶500 噸級駁船(45.0 m×10.8 m×1.6 m)。定義:T為船體舷側(cè)外板與橋墩壁面間距;L為兩橋墩中心間距;D為橋墩直徑;U為水流速度。采用靜水相對實驗,將扭矩傳感器安裝于船舶模型的重心上,并通過一根剛性支桿將固定支架與扭矩傳感器相連,橋墩模型固定在滑動平臺車下,平臺車在三相電磁調(diào)速電機的帶動下沿水槽頂圓形鋼管軌道勻速滑動,同時,采用扭矩傳感器測量船舶模型的艏搖力矩沿程變化情況。該模型的平面布置如圖1所示。

    圖1 模型平面布置示意Fig. 1 The layout of the model plane

    1.2 試驗結(jié)果分析

    船舶艏搖力矩歷時曲線如圖2所示。在圖2中,橫坐標為船艏與橋墩中心的縱向距離,橋墩中心處x=0 m,負值表示船舶位于橋墩上游,模型中船長為0.9 m。圖2(a)給出了在水流速度U=0.283 m/s時,單橋墩下,船與橋橫向間距不同時的艏搖力矩歷時曲線。圖2(b)給出了以船與橋的橫向間距為0.5D、兩橋墩中心間距不同時的艏搖力矩歷時曲線。

    圖2 船舶艏搖力矩歷時曲線Fig. 2 Ship's rolling moment curve

    從圖2(a)中可以看出,船舶下行經(jīng)過橋墩過程中,艏搖力矩呈正峰值→負峰值→正峰值的規(guī)律演進,且存在一定的波動。船舶以不同的船、橋橫向間距下行,艏搖力矩演進規(guī)律具有一致性。且隨著船、橋橫向間距的增加,艏搖力矩峰值減小。從圖2(b)中可以看出,對比單橋墩與串列橋墩工況,第一個艏搖力矩正峰值形成的位置相同,均為船艏到達上游橋墩時的位置。在橋墩間距為2.0D和3.0D時,船舶艏搖力矩演進規(guī)律與單橋墩工況類似,經(jīng)歷著正峰值→負峰值→正峰值的變化過程,但負峰值與第二個正峰值出現(xiàn)的時刻有所差異。當橋墩間距為4.0D和5.0D時,艏搖力矩開始呈現(xiàn)復雜的波動過程,在第一個正峰值和最后一個正峰值之間還存在多個小幅波動。將單橋墩與不同橋墩間距的串列橋墩艏搖力矩峰值的大小進行對比發(fā)現(xiàn),串列橋墩的第一個艏搖力矩正峰值均較單橋墩的大。隨著橋墩間距的增加,最大艏搖力矩負峰值和最后艏搖力矩正峰值減小。

    2 數(shù)學模型的建立及驗證

    本研究應用CFD中的FLUENT軟件,采用粘性流求解方法,模擬兩船舶會遇、追越、靠泊及經(jīng)過橋墩和岸壁時的非定常粘性流場,捕捉船舶行進過程中船舶與其他障礙物周圍的復雜流場細節(jié),對船、橋交匯運動進行二維數(shù)值模擬研究,揭示其中的水動力學機理。

    2.1 模型尺寸及邊界條件

    船舶、橋墩數(shù)學模型的參數(shù)和幾何比尺與其物理模型的一致。經(jīng)過數(shù)值計算比選后,確定計算域為80D×20D。船舶在距進口10D處釋放,保證船體與橋墩有較遠的初始距離,船體中軸線位于模型水平中心線上。上游橋墩中心距進口40D,與物理模型一致。采用速度進口和壓力出口邊界條件,將兩側(cè)壁面設(shè)定為對稱邊界,船體和橋墩表面設(shè)置為無滑移壁面邊界。計算域范圍如圖3所示。

    圖3 計算域的尺寸和邊界條件Fig. 3 Computational domain size and boundary conditions

    2.2 控制方程

    RNGk-ε模型可處理帶旋流及流線有較大畸變的流動[11]。本研究選取RNGk-ε模型進行串列橋墩繞流模擬和行進船舶周圍流場模擬。因此,在進行數(shù)值計算時,除滿足連續(xù)性方程和N-S方程外,還需滿足RNGk-ε模型的湍流動能k和湍流動能耗散率ε方程。

    2.3 數(shù)值方法

    2.3.1 網(wǎng)格劃分

    把橋區(qū)船舶二維平面運動簡化為1個三自由度問題(即沿X軸方向的縱向移動、橫漂及船體繞Z軸的轉(zhuǎn)動)。運動船體位移的尺寸遠大于網(wǎng)格的尺寸,易致使網(wǎng)格畸變率過大而導致計算不收斂,因此選擇彈性光順模型與局部重構(gòu)模型進行網(wǎng)格更新。為減少動網(wǎng)格更新時間,節(jié)約計算時間,采用多區(qū)域計算模型,如圖4所示。只將船寬外一定寬度范圍劃分為動網(wǎng)格區(qū)域,除船體邊界層網(wǎng)格使用四邊形網(wǎng)格外,均采用三角形非結(jié)構(gòu)網(wǎng)格,其他區(qū)域設(shè)定為靜止網(wǎng)格區(qū)域,采用四邊形結(jié)構(gòu)化網(wǎng)格,動網(wǎng)格區(qū)域與靜止網(wǎng)格區(qū)域之間采用交界面連接,實現(xiàn)數(shù)據(jù)傳遞。

    圖4 動網(wǎng)格區(qū)域劃分示意Fig. 4 Division diagram of the dynamic grid area

    2.3.2 用戶自定義函數(shù)UDF的使用

    為模擬船舶受橋墩紊流影響后的受力和運動情況,需采用動網(wǎng)格,將力轉(zhuǎn)化為速度和位移。本試驗采用用戶自定義函數(shù)UDF中的動網(wǎng)格宏“DEFINE_CG_MOTION”定義船舶在每一個時間步上的線速度及角速度來指定船舶運動,并通過“Compute_Force_And_Moment”宏獲取每1個時間步船體受力,以此計算下一時間步船體重心速度、角速度及重心位置,可實現(xiàn)水流與船舶運動的耦合。

    2.4 模型驗證

    2.4.1 繞流流場特性驗證

    串列橋墩繞流流場(L/D=4.0)如圖5所示。將PIV(粒子圖像測速法)物理模型試驗觀測的墩后尾流流場(如圖5(a)所示)與數(shù)學模型計算得到的墩后尾流流場(如圖5(b)所示)進行對比發(fā)現(xiàn),數(shù)學模型計算的結(jié)果與物理模型試驗的結(jié)果吻合得較好。在L/D=4.0間距下,兩者的上、下游橋墩均有渦體脫落,渦體位置和大小分布規(guī)律相似,上游圓柱分離剪切層在下游圓柱重新附著位置相同。表明:流場的數(shù)學模型計算結(jié)果能很好地與PIV物理模型試驗結(jié)果相吻合,得到合理的數(shù)值解。因此,該數(shù)學模型的網(wǎng)格和方法能夠滿足數(shù)值計算要求。

    2.4.2 繞流水力特征參數(shù)驗證

    圓柱繞流的重要參數(shù)有雷諾數(shù)Re、阻力系數(shù)CD、升力系數(shù)CL及斯特勞哈爾數(shù)St。本研究中均取Re=3.16×104,將計算得到的各參數(shù)與已有研究成果進行比較(見表1)。從表1中可以看出,本計算結(jié)果與其他研究數(shù)據(jù)[12-13]吻合得較好。

    圖5 串列橋墩繞流流場(L/D=4.0)Fig. 5 Flow field diagram of bridge piers (L/D=4.0)

    數(shù)據(jù)來源雷諾數(shù)上游橋墩阻力系數(shù)斯特勞哈爾數(shù)下游橋墩阻力系數(shù)斯特勞哈爾數(shù)Igarashi[12]3.55×1041.250.1850.420.185Kitagawa[13]2.20×1041.220.1890.490.189本研究3.16×1041.200.1860.460.186

    2.4.3 船舶受力特性驗證

    建立與物理模型試驗相對應的計算工況,對數(shù)學模型船舶的受力進行驗證。計算工況為:T=0.5D,U=0.283 m/s,L/D=3.0及單橋墩情況。單橋墩和串列橋墩物理模型測量值與數(shù)值模擬計算值的結(jié)果對比如圖6所示。從圖6中可以看出,兩工況的演進規(guī)律吻合得較好。由于物理模型力矩傳感器的數(shù)據(jù)采樣頻率限制和橋墩周圍自由液面變化的影響,艏搖力矩的測量值與計算值存在著一定的偏差,但數(shù)值模擬結(jié)果總體能夠反映船、橋交匯過程中船體所受艏搖力矩的變化規(guī)律。

    3 數(shù)值計算結(jié)果分析

    3.1 計算工況

    考慮汛期急流最不利工況,同時,根據(jù)學者橋區(qū)船舶航行安全航速的研究成果和海事部門的相關(guān)規(guī)定,取水流速度2.0 m/s,船舶速度2.5 m/s,船、橋橫向間距T=0.5D,分別計算單橋墩和兩橋墩中心間距(L=2.0D,3.0D,4.0D,5.0D,9.0D和13.0D)時船舶經(jīng)過橋墩沿程所受的艏搖力矩。

    圖6 艏搖力矩的測量值與計算值的對比Fig. 6 Comparison between the measurement value and the calculated value of the shaking moment

    3.2 船舶周圍壓力場分析

    船、橋交匯過程不同時刻的壓力云圖如圖7所示。暖色云圖為正壓,冷色云圖為負壓。本研究選擇幾個典型位置進行分析。如圖7(a)所示,船體接近橋墩過程中,船體右前端伴流斜向橋墩運動。同時,由于船體邊界的靠近,抑制了墩前駐點周圍流線的發(fā)展,船與橋間水體不斷受到擠壓,加之墩前方推船流的共同作用,使得船體右舷前端的壓力大于左舷的,船艏受排斥,遠離橋墩,航向改變。在船艏即將到達橋墩時,船與橋之間的流道最窄,排斥力最大,艏搖力矩達到一個正峰值。當船艏進入橋墩尾流區(qū)、船尾仍位于橋墩上游時,橋墩上游船體邊界依然擠壓船與橋之間的水體,如圖7(b)所示,橋墩上游船體右舷側(cè)受正壓作用,推離橋墩。受船艏部分伴流的帶動,讓水流出船與橋之間的流道后,加速向墩后擴散,加劇了尾渦的發(fā)展,較無船舶時的墩后負壓區(qū)增加。墩下游船身右舷側(cè)在負壓作用下,船艏吸向橋墩側(cè)。由于船體右舷側(cè)船艏受吸引,船尾受排斥,對船體形成順時針力矩,艏搖力矩表現(xiàn)為負值。

    船尾遠離橋墩過程中,部分尾渦撞擊船體后或吸附于船體,或與周圍渦體發(fā)生更強烈的混參。同時,在船尾伴流的帶動下,加速了船尾與橋墩間流道水流,因而墩后方負壓區(qū)域發(fā)展更加顯著。船尾右舷側(cè)處于負壓區(qū)內(nèi),形成艏搖力矩正峰值。

    圖7 L/D=5.0船橋交會不同時刻的壓力云圖Fig. 7 Pressure cloud diagram at different times of bridge crossing when L/D=5.0

    3.3 船舶受力分析

    不同橋墩間距下,船舶經(jīng)過橋墩沿程所受艏搖力矩如圖8所示。從圖8中可以看出,經(jīng)過單橋墩及不同橋墩間距下的串列橋墩船舶艏搖力矩變化過程具有一定差異:在橋墩中心間距L/D=2.0時,艏搖力矩曲線演進規(guī)律與單橋墩時的較為一致,均以正峰值→負峰值→正峰值規(guī)律演進,峰值出現(xiàn)相對位置也相同;當L/D≥3.0后,力矩曲線呈多個正負峰值交替出現(xiàn),其中,第一個正峰值及最后一個正峰值出現(xiàn)的相對位置相同。從圖7中可以看出,由于沿船體右舷交替分布著正、負壓力區(qū),且隨著船體的下行,船體各位置上的壓力區(qū)大小和強度實時變化,因此,介于第一正峰值與最后正峰值之間,不同的船橋橫向間距L/D,會使艏搖力矩呈現(xiàn)不同的波動情況。

    艏搖力矩峰值與不同橋墩間距L/D的關(guān)系 曲線如圖9所示。從圖9中可以看出,串列橋墩的第一個艏搖力矩正峰值均大于單橋墩的。隨著橋墩間距的增加,峰值逐漸減小。表明:隨著橋墩間距的增加,下游橋墩對上游橋墩的影響逐漸減弱,趨于2個孤立的單墩狀態(tài);串列橋墩最后一個艏搖力矩正峰值隨橋墩間距的增加而減小。在2.0≤L/D≤3.0時,串列橋墩的艏搖力矩峰值大于單橋墩的;而當4.0≤L/D≤5.0時,串列橋墩的艏搖力矩峰值小于單橋墩的。因此,當串列橋墩間距L/D≤3.0時,上、下游橋墩的相互作用會對船舶安全航行帶來不利影響。特別是在L/D=2.0時,串列橋墩第一個艏搖力矩正峰值達到了單橋墩的1.87倍。當L/D=3.0時,串列橋墩艏搖力矩負峰值相比單橋墩的增加了45%。隨著間距的繼續(xù)增加,串列橋墩對船舶水動力的干擾逐漸趨于單橋墩的。但在串列橋墩間距達5.0D時,部分力矩峰值比單橋墩的更低。

    圖8 不同L/D下的船舶艏搖力矩Fig. 8 Ship’s shaking moments at different L/D

    圖9 艏搖力矩峰值與橋墩間距關(guān)系曲線Fig. 9 The relationship between the peak moment of shaking torque and pier spacing

    研究兩橋墩間距達到何值時,開始類似于單橋墩情況是非常必要的。通過計算發(fā)現(xiàn),兩橋墩間距L/D≥9.0后,船舶通過上游橋墩過程中,船舶艏搖力矩曲線與單橋墩的十分貼合,如圖10所示。表明:下游橋墩對上游橋墩的水動力不構(gòu)成干擾。對于L/D=9.0工況,在船尾離開上游橋墩時,達到串列橋墩的艏搖力矩正峰值(M=1.778 N·m)卻比單橋墩的大得多,約為單橋墩的1.6倍。這是由于此處船長恰好等于橋墩間距,此時船尾將離開上游橋墩而船艏恰好到達下游橋墩,形成的力矩峰值為上游橋墩對船舶的最后正峰值與下游橋墩對船舶第一個正峰值的疊加。觀察L/D=13.0的串列橋墩艏搖力矩曲線發(fā)現(xiàn),船舶經(jīng)過下游橋墩相當于復演了經(jīng)過上游橋墩的過程,2個過程形成的艏搖力矩曲線一致且與單橋墩吻合度較高。表明:橋墩間距達13.0D后,串列橋墩對船舶水動力的干擾等同于2個孤立的單橋墩。

    4 結(jié)論

    本研究分析了不同橋墩間距下船、橋交匯過程中船舶沿程所受艏搖力矩演進規(guī)律,得出的結(jié)論為:

    圖10 L/D=9.0,13.0和單橋墩下的艏搖力矩曲線Fig. 10 Single-bridge pier swaying moment curve at L/D=9.0,13.0,respectively

    1) 物理模型試驗和數(shù)值計算結(jié)果均表明:駛經(jīng)單橋墩的船舶艏搖力矩曲線呈正峰值、負峰值、正峰值的規(guī)律演進。船體邊界靠近橋墩,擠壓船艏與橋墩之間的水體以及墩前正向橫流(推船流)作用是艏搖力矩第一個正峰值形成的原因。船尾即將離開橋墩時,船尾附近水流吸入船尾與橋墩間的過流通道向下游擴散,水體加速流動,負壓增強,吸引船尾形成第二個艏搖力矩正峰值。2個正峰值之間的負峰值的出現(xiàn)則與船體相對于橋墩特定位置,渦體沿船身分布強弱相關(guān)。

    2) 駛經(jīng)串列橋墩的船舶艏搖力矩演進規(guī)律與單橋墩的存在著差異。當兩橋墩中心間距較小、為2.0D時,其艏搖力矩演進規(guī)律與單橋墩的一致。兩橋墩中心間距達3.0D時,其艏搖力矩歷時曲線在第一個正峰值與最后正峰值之間出現(xiàn)多個波峰波谷波動現(xiàn)象。第一個正峰值與最后正峰值形成的原因與單橋墩的相同,中間的小幅波動與船舶相對于兩橋墩的位置有關(guān)。兩橋墩間距達9.0D時,下游橋墩對上游橋墩的水動力干擾消失。在兩橋墩間距大于13.0D時,船舶駛經(jīng)上、下游橋墩,船舶艏搖力矩相當于連續(xù)經(jīng)過2個孤立的單橋墩。

    猜你喜歡
    船尾船體橋墩
    船體行駛過程中的壓力監(jiān)測方法
    午睡
    船尾外形對超聲速彈丸減阻特性的影響研究
    橋墩加固方案分析研究
    生命之船
    軌道交通預制拼裝橋墩設(shè)計和施工
    上海公路(2018年1期)2018-06-26 08:37:36
    中國船尾舵開創(chuàng)航行新時代
    中國船檢(2018年3期)2018-03-29 06:38:41
    焊接殘余應力對船體結(jié)構(gòu)疲勞強度的影響分析
    焊接(2015年9期)2015-07-18 11:03:51
    赴美軍“仁慈”號醫(yī)院船駐船體會
    水下爆炸氣泡作用下船體總縱強度估算方法
    激情视频va一区二区三区| 亚洲av成人精品一二三区| 人妻一区二区av| 国产精品99久久99久久久不卡 | 女性被躁到高潮视频| 制服人妻中文乱码| av免费在线看不卡| av福利片在线| 人成视频在线观看免费观看| 日本-黄色视频高清免费观看| 国产毛片在线视频| 国产日韩欧美在线精品| 成人无遮挡网站| 欧美精品国产亚洲| 91午夜精品亚洲一区二区三区| h视频一区二区三区| 97人妻天天添夜夜摸| 欧美少妇被猛烈插入视频| 香蕉国产在线看| 九草在线视频观看| 一个人免费看片子| 免费看av在线观看网站| 卡戴珊不雅视频在线播放| 国产视频首页在线观看| 哪个播放器可以免费观看大片| 国产在线视频一区二区| 男人舔女人的私密视频| 天美传媒精品一区二区| 国产精品欧美亚洲77777| 精品午夜福利在线看| 一本色道久久久久久精品综合| 日韩伦理黄色片| 99香蕉大伊视频| 宅男免费午夜| 国产老妇伦熟女老妇高清| 日本vs欧美在线观看视频| 日本-黄色视频高清免费观看| 免费大片黄手机在线观看| 精品亚洲乱码少妇综合久久| 少妇人妻精品综合一区二区| 9191精品国产免费久久| 欧美性感艳星| 各种免费的搞黄视频| 久久久久久久久久人人人人人人| 国产毛片在线视频| 日韩制服丝袜自拍偷拍| 日本欧美国产在线视频| 免费人妻精品一区二区三区视频| 国产白丝娇喘喷水9色精品| 国产精品 国内视频| 日韩伦理黄色片| 国产精品久久久av美女十八| 女性被躁到高潮视频| 毛片一级片免费看久久久久| 最新中文字幕久久久久| 国产av码专区亚洲av| 亚洲成av片中文字幕在线观看 | 26uuu在线亚洲综合色| 国产成人精品福利久久| 欧美激情国产日韩精品一区| 亚洲精品一二三| 80岁老熟妇乱子伦牲交| 巨乳人妻的诱惑在线观看| 亚洲国产成人一精品久久久| 久久人人爽人人片av| 国产精品一区www在线观看| 黄色视频在线播放观看不卡| 日韩精品有码人妻一区| 久久精品国产自在天天线| 免费看不卡的av| 国产日韩欧美视频二区| 少妇人妻精品综合一区二区| 亚洲精品色激情综合| 人体艺术视频欧美日本| 欧美老熟妇乱子伦牲交| 日韩精品免费视频一区二区三区 | av视频免费观看在线观看| 免费人成在线观看视频色| 青春草视频在线免费观看| 亚洲国产毛片av蜜桃av| 精品99又大又爽又粗少妇毛片| 久热久热在线精品观看| 日韩一区二区三区影片| 久久精品夜色国产| 成人漫画全彩无遮挡| 黑人高潮一二区| 免费观看av网站的网址| 久久国产精品男人的天堂亚洲 | kizo精华| 国产乱人偷精品视频| 欧美成人午夜精品| 久久午夜福利片| 国产精品女同一区二区软件| 久久人人爽人人爽人人片va| 草草在线视频免费看| av播播在线观看一区| 精品熟女少妇av免费看| 插逼视频在线观看| 免费观看av网站的网址| 日本免费在线观看一区| 精品一区二区三区四区五区乱码 | av免费在线看不卡| 桃花免费在线播放| 久久人人爽人人爽人人片va| 精品一区二区免费观看| 亚洲成人av在线免费| 国产一区有黄有色的免费视频| 婷婷色av中文字幕| 超碰97精品在线观看| 热re99久久精品国产66热6| 国内精品宾馆在线| 亚洲婷婷狠狠爱综合网| 99九九在线精品视频| 最近中文字幕高清免费大全6| 伦理电影免费视频| 亚洲精品乱久久久久久| 日本色播在线视频| 日日撸夜夜添| 美女主播在线视频| 欧美日韩视频精品一区| 赤兔流量卡办理| 麻豆乱淫一区二区| 五月天丁香电影| 久久久久视频综合| 国产免费一级a男人的天堂| 十八禁网站网址无遮挡| 捣出白浆h1v1| 宅男免费午夜| 久久久精品免费免费高清| 国产又爽黄色视频| videos熟女内射| 国产一区有黄有色的免费视频| 免费少妇av软件| 国产午夜精品一二区理论片| 欧美激情国产日韩精品一区| av有码第一页| av在线老鸭窝| 日韩精品免费视频一区二区三区 | 国产欧美日韩一区二区三区在线| 少妇熟女欧美另类| av电影中文网址| 最近中文字幕高清免费大全6| 国产精品一区二区在线不卡| av又黄又爽大尺度在线免费看| 亚洲av成人精品一二三区| 一区在线观看完整版| 国产黄色免费在线视频| 男的添女的下面高潮视频| 香蕉丝袜av| 最近中文字幕2019免费版| 黄色配什么色好看| 亚洲精品一二三| a级片在线免费高清观看视频| 精品一品国产午夜福利视频| 在线观看免费视频网站a站| 一级爰片在线观看| 国产日韩欧美亚洲二区| 国产亚洲av片在线观看秒播厂| 成年动漫av网址| 国产欧美亚洲国产| 亚洲国产精品一区二区三区在线| 国产精品国产三级专区第一集| 免费在线观看黄色视频的| 精品一区二区免费观看| 一边亲一边摸免费视频| 一本色道久久久久久精品综合| 免费高清在线观看视频在线观看| 久久久久久久久久久久大奶| 国产精品 国内视频| 两性夫妻黄色片 | 男女啪啪激烈高潮av片| freevideosex欧美| 26uuu在线亚洲综合色| 爱豆传媒免费全集在线观看| 伦理电影免费视频| 91精品国产国语对白视频| 午夜老司机福利剧场| 青春草国产在线视频| 老熟女久久久| 午夜免费男女啪啪视频观看| 国产一区二区三区综合在线观看 | 国产精品偷伦视频观看了| 欧美日韩视频精品一区| 一个人免费看片子| 久久久a久久爽久久v久久| 国产乱来视频区| 男女国产视频网站| 最近中文字幕高清免费大全6| 欧美少妇被猛烈插入视频| 日韩大片免费观看网站| 草草在线视频免费看| 亚洲美女黄色视频免费看| 久久久久久久久久久久大奶| 大香蕉久久成人网| 搡女人真爽免费视频火全软件| 中国三级夫妇交换| 91久久精品国产一区二区三区| 2018国产大陆天天弄谢| 在线观看美女被高潮喷水网站| 99九九在线精品视频| 亚洲五月色婷婷综合| 狠狠精品人妻久久久久久综合| 丁香六月天网| 亚洲第一av免费看| 欧美精品高潮呻吟av久久| 狠狠精品人妻久久久久久综合| 天天躁夜夜躁狠狠躁躁| 女人被躁到高潮嗷嗷叫费观| 久久久久久人人人人人| 中文字幕亚洲精品专区| 日韩制服丝袜自拍偷拍| 日韩电影二区| 一区二区日韩欧美中文字幕 | 看十八女毛片水多多多| 肉色欧美久久久久久久蜜桃| 一级,二级,三级黄色视频| 最新中文字幕久久久久| 国产国拍精品亚洲av在线观看| 青青草视频在线视频观看| 久久精品国产综合久久久 | 久久av网站| 精品国产一区二区三区久久久樱花| 丝袜人妻中文字幕| av片东京热男人的天堂| 国产亚洲最大av| 日韩av免费高清视频| 久久精品久久精品一区二区三区| 99久国产av精品国产电影| 久久久久久人人人人人| 日韩欧美一区视频在线观看| 久久久国产精品麻豆| 亚洲经典国产精华液单| 国产成人精品无人区| 国产极品粉嫩免费观看在线| 一本—道久久a久久精品蜜桃钙片| 亚洲国产精品一区三区| 天美传媒精品一区二区| 观看美女的网站| 精品卡一卡二卡四卡免费| 国产精品久久久久成人av| 欧美日韩av久久| 青春草国产在线视频| 女的被弄到高潮叫床怎么办| 日韩视频在线欧美| 人人妻人人爽人人添夜夜欢视频| 亚洲人成77777在线视频| 久久久亚洲精品成人影院| 亚洲欧洲国产日韩| 爱豆传媒免费全集在线观看| 免费大片18禁| 美女内射精品一级片tv| 老司机影院成人| 少妇人妻精品综合一区二区| 欧美 日韩 精品 国产| 男人添女人高潮全过程视频| 国产色爽女视频免费观看| 中国美白少妇内射xxxbb| av.在线天堂| 国产熟女午夜一区二区三区| 国产成人精品婷婷| 卡戴珊不雅视频在线播放| 国产精品一区www在线观看| av在线app专区| 亚洲国产日韩一区二区| 精品人妻偷拍中文字幕| 日韩制服骚丝袜av| 国产av国产精品国产| 成人漫画全彩无遮挡| 又黄又粗又硬又大视频| 婷婷色综合www| 天天躁夜夜躁狠狠躁躁| 三上悠亚av全集在线观看| 久久亚洲国产成人精品v| 少妇的逼水好多| 免费观看在线日韩| 免费观看av网站的网址| 成人手机av| 哪个播放器可以免费观看大片| 国产激情久久老熟女| 婷婷色综合www| 九九爱精品视频在线观看| 日韩成人伦理影院| 男女下面插进去视频免费观看 | 中文乱码字字幕精品一区二区三区| 日韩制服丝袜自拍偷拍| 18禁国产床啪视频网站| 国产成人精品无人区| 国产精品三级大全| 18禁在线无遮挡免费观看视频| 欧美人与性动交α欧美软件 | 永久网站在线| 大香蕉久久成人网| 日韩大片免费观看网站| 亚洲四区av| 婷婷成人精品国产| 女人久久www免费人成看片| 这个男人来自地球电影免费观看 | 又粗又硬又长又爽又黄的视频| 久久青草综合色| 久久久久精品性色| 熟妇人妻不卡中文字幕| 久久这里有精品视频免费| 国产白丝娇喘喷水9色精品| 18在线观看网站| 777米奇影视久久| 久久 成人 亚洲| 狂野欧美激情性bbbbbb| 亚洲精品一区蜜桃| 国产午夜精品一二区理论片| 国产免费福利视频在线观看| 亚洲av成人精品一二三区| 国产精品一区二区在线不卡| 国产在视频线精品| 亚洲少妇的诱惑av| 女人久久www免费人成看片| 国产熟女午夜一区二区三区| 成人漫画全彩无遮挡| 韩国av在线不卡| 欧美变态另类bdsm刘玥| 边亲边吃奶的免费视频| 亚洲精品色激情综合| 高清欧美精品videossex| 久热久热在线精品观看| 亚洲精品久久成人aⅴ小说| 色94色欧美一区二区| 国产精品偷伦视频观看了| av一本久久久久| 国产精品熟女久久久久浪| 亚洲第一av免费看| 成人亚洲欧美一区二区av| 波多野结衣一区麻豆| 伦理电影免费视频| 免费黄色在线免费观看| 美女中出高潮动态图| 永久网站在线| 亚洲综合色网址| 国语对白做爰xxxⅹ性视频网站| 99久久人妻综合| 亚洲经典国产精华液单| 精品国产乱码久久久久久小说| 亚洲,一卡二卡三卡| 999精品在线视频| 九草在线视频观看| 欧美日韩一区二区视频在线观看视频在线| 日本与韩国留学比较| 女人精品久久久久毛片| 丰满迷人的少妇在线观看| 国产精品熟女久久久久浪| 久久久精品区二区三区| 国产成人精品一,二区| 久久久久久久久久成人| 成年美女黄网站色视频大全免费| 色婷婷av一区二区三区视频| 日本av手机在线免费观看| 国产精品一二三区在线看| 色5月婷婷丁香| 日本av免费视频播放| 久久综合国产亚洲精品| 国产精品久久久久久精品古装| 亚洲精华国产精华液的使用体验| 国产福利在线免费观看视频| 亚洲欧洲日产国产| 亚洲五月色婷婷综合| 久久久精品免费免费高清| 熟女av电影| 久久久亚洲精品成人影院| 日日爽夜夜爽网站| 最近2019中文字幕mv第一页| 老司机亚洲免费影院| 乱码一卡2卡4卡精品| av在线老鸭窝| 两个人免费观看高清视频| 国产免费福利视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久久精品久久精品一区二区三区| 国产日韩欧美视频二区| 精品一品国产午夜福利视频| 永久网站在线| 少妇的丰满在线观看| 国语对白做爰xxxⅹ性视频网站| 九九爱精品视频在线观看| 国产欧美亚洲国产| 成人黄色视频免费在线看| 观看美女的网站| 亚洲成人一二三区av| 久久av网站| 两性夫妻黄色片 | 毛片一级片免费看久久久久| 纯流量卡能插随身wifi吗| 久久久久国产网址| 一区在线观看完整版| av天堂久久9| 狂野欧美激情性bbbbbb| 精品亚洲成国产av| 大陆偷拍与自拍| 女人精品久久久久毛片| 老司机亚洲免费影院| 免费看不卡的av| 国产精品一国产av| 91精品国产国语对白视频| 精品亚洲成国产av| 久久免费观看电影| 欧美精品国产亚洲| 国内精品宾馆在线| 久久久a久久爽久久v久久| 99国产综合亚洲精品| 各种免费的搞黄视频| www日本在线高清视频| 少妇的逼水好多| 久久精品国产自在天天线| 国产成人a∨麻豆精品| 伊人久久国产一区二区| 九九在线视频观看精品| 精品久久蜜臀av无| 国产乱来视频区| 亚洲人成77777在线视频| xxx大片免费视频| 国产乱来视频区| 90打野战视频偷拍视频| 亚洲av男天堂| 97人妻天天添夜夜摸| 亚洲国产毛片av蜜桃av| 纵有疾风起免费观看全集完整版| 国产高清国产精品国产三级| 亚洲精品成人av观看孕妇| 精品少妇黑人巨大在线播放| 亚洲国产精品专区欧美| 在现免费观看毛片| 国产精品国产av在线观看| 最近最新中文字幕大全免费视频 | 肉色欧美久久久久久久蜜桃| 一级毛片 在线播放| 毛片一级片免费看久久久久| 精品亚洲成国产av| 欧美最新免费一区二区三区| 九色成人免费人妻av| 热99国产精品久久久久久7| 国产伦理片在线播放av一区| 日本色播在线视频| 一本久久精品| 热99久久久久精品小说推荐| 中文字幕人妻丝袜制服| 18禁裸乳无遮挡动漫免费视频| 国产精品99久久99久久久不卡 | 在线观看免费日韩欧美大片| 欧美少妇被猛烈插入视频| 少妇 在线观看| 午夜福利视频在线观看免费| 免费女性裸体啪啪无遮挡网站| 色哟哟·www| 天天躁夜夜躁狠狠躁躁| 91午夜精品亚洲一区二区三区| 视频中文字幕在线观看| 亚洲国产成人一精品久久久| 中文字幕另类日韩欧美亚洲嫩草| 久久综合国产亚洲精品| 亚洲精品中文字幕在线视频| 亚洲成人av在线免费| 久久久久久伊人网av| 一本大道久久a久久精品| 国产淫语在线视频| 超色免费av| 在线观看免费高清a一片| 亚洲,一卡二卡三卡| 欧美成人精品欧美一级黄| 美女脱内裤让男人舔精品视频| 激情五月婷婷亚洲| 亚洲美女搞黄在线观看| 交换朋友夫妻互换小说| 亚洲,一卡二卡三卡| 国产淫语在线视频| 欧美bdsm另类| a级毛片黄视频| 在线精品无人区一区二区三| 国产国拍精品亚洲av在线观看| 999精品在线视频| 肉色欧美久久久久久久蜜桃| 国产精品一区二区在线不卡| 高清毛片免费看| 成人毛片60女人毛片免费| 日韩成人av中文字幕在线观看| 国产亚洲精品第一综合不卡 | 久久精品久久精品一区二区三区| 91成人精品电影| 欧美精品人与动牲交sv欧美| 综合色丁香网| 亚洲第一av免费看| 哪个播放器可以免费观看大片| 99久久精品国产国产毛片| 国产毛片在线视频| 久久久久视频综合| 欧美成人精品欧美一级黄| 日本黄大片高清| 亚洲,一卡二卡三卡| h视频一区二区三区| 亚洲欧洲国产日韩| 最近的中文字幕免费完整| 日韩熟女老妇一区二区性免费视频| 夫妻午夜视频| 亚洲美女搞黄在线观看| 婷婷色麻豆天堂久久| 99久久综合免费| 国产日韩欧美亚洲二区| 伊人亚洲综合成人网| 国精品久久久久久国模美| 国产乱人偷精品视频| 日本与韩国留学比较| 大片免费播放器 马上看| 日本欧美视频一区| 精品一区二区免费观看| 中文乱码字字幕精品一区二区三区| 热99久久久久精品小说推荐| 在线观看三级黄色| 亚洲精品第二区| 免费观看av网站的网址| 9色porny在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲av国产av综合av卡| 我要看黄色一级片免费的| 免费大片黄手机在线观看| 国国产精品蜜臀av免费| 韩国精品一区二区三区 | 黄色怎么调成土黄色| 国产淫语在线视频| 少妇人妻久久综合中文| 丝袜在线中文字幕| 国产永久视频网站| 一本久久精品| 捣出白浆h1v1| 两个人看的免费小视频| 亚洲国产看品久久| 国产男人的电影天堂91| 自线自在国产av| 国产黄色免费在线视频| 2021少妇久久久久久久久久久| 国产免费现黄频在线看| 日本猛色少妇xxxxx猛交久久| 国产欧美日韩一区二区三区在线| 寂寞人妻少妇视频99o| 欧美日韩视频高清一区二区三区二| 涩涩av久久男人的天堂| 亚洲精品,欧美精品| 中文字幕精品免费在线观看视频 | 久久久久久伊人网av| 欧美3d第一页| 欧美精品一区二区大全| 99九九在线精品视频| 成人午夜精彩视频在线观看| 美女主播在线视频| 亚洲国产精品成人久久小说| 国产一区二区三区综合在线观看 | 成年动漫av网址| 美女国产视频在线观看| 欧美成人精品欧美一级黄| 亚洲一区二区三区欧美精品| www日本在线高清视频| 免费大片黄手机在线观看| 黑人猛操日本美女一级片| 丁香六月天网| 少妇被粗大猛烈的视频| 国产精品一二三区在线看| 纵有疾风起免费观看全集完整版| 美女福利国产在线| 日韩成人伦理影院| 丰满乱子伦码专区| 一级爰片在线观看| 婷婷色综合大香蕉| 黄网站色视频无遮挡免费观看| 免费观看在线日韩| 久久午夜福利片| 色94色欧美一区二区| 校园人妻丝袜中文字幕| 国产毛片在线视频| 久久热在线av| 赤兔流量卡办理| √禁漫天堂资源中文www| 久久婷婷青草| 久久人人爽人人爽人人片va| 在线观看一区二区三区激情| 久久久久久伊人网av| 丝瓜视频免费看黄片| 搡女人真爽免费视频火全软件| 在线亚洲精品国产二区图片欧美| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人久久小说| 99视频精品全部免费 在线| av在线播放精品| kizo精华| 91久久精品国产一区二区三区| 天天影视国产精品| 久久这里只有精品19| 高清毛片免费看| 五月伊人婷婷丁香| 亚洲精品中文字幕在线视频| 欧美亚洲 丝袜 人妻 在线| 亚洲国产精品成人久久小说| 免费观看在线日韩| 99久久精品国产国产毛片| 狠狠婷婷综合久久久久久88av| 咕卡用的链子| 一级毛片黄色毛片免费观看视频| 一区二区三区乱码不卡18| 最近中文字幕2019免费版| 亚洲国产最新在线播放| 黑丝袜美女国产一区| 咕卡用的链子| 熟女电影av网| 久久国内精品自在自线图片| 你懂的网址亚洲精品在线观看| 我要看黄色一级片免费的| 日本-黄色视频高清免费观看| 高清欧美精品videossex| 免费观看在线日韩| av女优亚洲男人天堂| 一级,二级,三级黄色视频|