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

    多車道寬幅懸浮隧道管體基于CFD的繞流特性與縮尺效應

    2025-02-24 00:00:00陳志威劉宏專易壯鵬
    應用力學學報 2025年1期

    摘 要:為滿足通行的需求懸浮隧道管體宜設計成大尺度的多車道寬幅截面,洋流是施工及運營階段所面臨的主要環(huán)境荷載。洋流作用下寬幅截面的繞流特性,以及縮尺試驗或數值模擬時的尺寸效應值得重點關注。以橢圓形、圓端形、八邊形3種按雙向八車道設計的寬幅截面為對象,并以圓形截面為對照,用CFD軟件fluent建立了非定常流下的二維數值模型。研究了不同形狀及不同縮尺比的寬幅截面在均勻流下的復雜繞流特性,建立了各類截面升力系數幅值、阻力系數均值和力矩系數幅值與流速的關系圖,獲取了不同縮尺比下三分力系數的時程曲線,并在此基礎上給出了穩(wěn)定周期內的渦量云圖。結果表明:

    當流速增大時,寬幅截面三分力系數變化不大,圓形截面阻力系數、升力系數則明顯減??;3種寬幅截面中,八邊形升力系數、阻力系數和力矩系數均最大,脫渦頻率最小,橢圓形則與之相反;各類寬幅截面不同縮尺比之間的三分力系數存在不同程度的差異,需考慮尺寸效應帶來的影響;尾渦脫落對結構所受三分力有一定影響,阻力系數與渦量云圖中尾渦強度、分離點有關。

    關鍵詞:懸浮隧道;多車道寬幅截面;繞流特性;三分力系數;縮尺效應;CFD數值模擬

    中圖分類號:U451" 文獻標志碼:A

    文章編號:1000-4939(2025)01-0223-13

    Flow characteristics and scale effects of multi-lane wide tube of

    submerged floating tunnel based on CFD

    CHEN Zhiwei,LIU Hongzhuan,YI Zhuangpeng

    (School of Civil Engineering,Changsha University of Science amp; Technology,410114 Changsha,China)

    Abstract:To meet the traffic demand,the tube of the submerged floating tunnel (SFT) should be designed as a large-scale multi-lane wide section.And the current is the main environmental load faced during the construction and service stages.The flow characteristics of the wide section under the action of ocean currents,as well as the scale effect in scale test or numerical simulation,deserve special attention.The two-dimensional numerical models under unsteady flow are established by using CFD software fluent,by taking the elliptical,round-end and octagonal wide sections designed according to two-way eight lanes as the objects,and taking the circular section as the contrast.The complex flow characteristics of wide sections with different shapes and different scale ratios under uniform flow are studied.The relationship between lift coefficient amplitude,drag coefficient mean-value,moment coefficient amplitude and flow velocity for various sections is established.The time history curves of the three component coefficient under different scale ratios are obtained.On this basis,the vorticity cloud chart in a stable period is given.The results show that the three component coefficients of wide section change little in the increasing course of the flow velocity,while the drag and lift coefficients of circular section decrease significantly.For the three wide sections,the octagonal section has the largest lift coefficient,drag coefficient and moment coefficient,but the smallest vortex shedding frequency,while the elliptical section is just the opposite.The three component coefficients for various wide sections and scale ratios have differences in different degrees,so the influence of scale effect should be considered.The wake shedding has a certain influence on the three component forces,and the drag coefficient is related to the wake strength and separation point in the vorticity cloud chart.

    Key words:SFT;multi-lane wide section;flow characteristic;three-component coefficient;scale effect;CFD numerical simulation

    懸浮隧道[1-2]是一種新型的懸浮于水下一定深度的超長交通結構物,利用自身重力、浮力和錨索支撐力的共同作用實現受力平衡,在超長深水海峽通道中有很好的應用前景,其流固耦合特性與結構體系的力學行為近年來受到了廣泛的關注[3-6],并成為21世紀最具競爭力的跨海交通結構[7-8]。洋流是懸浮隧道施工及運營階段面臨的主要環(huán)境荷載,也是引起結構渦激響應[9-10]的主要原因。厘清洋流作用下懸浮隧道管體的繞流特性與流體作用力是研究各類流致振動的前提,而圓柱繞流問題[11-14]是研究懸浮隧道繞流特性和尺寸效應的理論基礎。

    圓柱繞流問題有著廣泛的工程應用背景,研究者們對其進行了許多研究,方法以試驗和數值模擬為主。如:梁圣召等[15]采用CFD數值模擬與試驗結合的方法,研究了圓柱在不同長徑比和雷諾數工況下的流動特性;陳松濤等[16]利用數值模擬探究了自由液面與自由端對典型鈍體繞流的影響,結果表明自由液面與自由端的存在會顯著改變流場分布;高洋洋等[17]基于CFD模擬了不同雷諾數下傾斜圓柱繞流,結果表明隨著圓柱傾角的增大,圓柱尾流渦旋脫落受到明顯抑制;劉闖等[18]采用大渦數值模擬了高雷諾數下圓柱尾渦模式,發(fā)現高雷諾數下圓柱繞流的尾跡變化呈現出不穩(wěn)定特征;李聰洲等[19]采用改進后的延遲分離渦方法模擬了高雷諾數下圓柱與方柱的繞流形式,數值結果與試驗結果吻合良好。隨著圓柱繞流問題的研究越來越深入,將圓柱繞流的研究成果用于懸浮隧道截面設計和優(yōu)化具有重要的理論意義和工程價值。

    近年來,研究者采用基于CFD的數值模擬方法結合試驗方法,對懸浮隧道的繞流現象做了許多研究。如:王長春[20]采用數值方法模擬了不同工況下洋流與懸浮隧道管段的相互作用規(guī)律,并與管段模型的試驗結果進行了對比驗證;秦銀剛等[21]采用相似準則對懸浮隧道管體進行了節(jié)段模型試驗,研究了管體受到洋流沖擊時的力學特性;王廣地等[22]基于RNG k-ε湍流模型采用數值模擬方法研究了斷面形狀對懸浮隧道管段升力系數、阻力系數的影響規(guī)律;郭曉玲等[23]基于SST k-ω湍流模型的雷諾平均方程,利用數值模擬方法對不同雷諾數下懸浮隧道管體的渦激振動問題進行了研究;曾繁旭等[24]采用CFD軟件中Fluent建立了均勻流作用下懸浮隧道不同管體截面的運動模型,研究了結構運動行為與流場參數之間的關系。這些研究為進一步探尋懸浮隧道管體的繞流特性奠定了基礎。未來懸浮隧道為滿足通行的需求,管體宜設計成大尺度的多車道寬幅截面[25],相應寬幅截面的升力系數、阻力系數、力矩系數和尾渦脫落規(guī)律相對于已有研究涉及的圓形截面或小尺度截面將發(fā)生變化。同時,為獲取大尺寸多車道寬幅懸浮隧道截面繞流特性,試驗一般為某一特定范圍雷諾數的縮尺試驗,而1∶1試驗難以實現。不同雷諾數對應的縮尺試驗或數值模型的尺寸效應同樣值得關注。

    本研究以橢圓形、圓端形、八邊形3種按雙向八車道設計的寬幅懸浮隧道截面為對象,并以圓形截面為對照,用CFD軟件Fluent建立了非定常流下的二維數值模型。建立各類截面不同縮尺比的升力系數、阻力系數和力矩系數與流速的關系,獲取三分力系數的時程曲線與穩(wěn)定周期內的渦量云圖。在此基礎上,研究多車道寬幅懸浮隧道截面在均勻流下的繞流特性以及不同縮尺比之間的尺寸效應,以期為懸浮隧道截面設計提供參考依據。

    1 研究方法

    1.1 流域設置與參數

    建立寬幅截面管體繞流的流場區(qū)域及參數條件如下:斷面寬度設為B,計算域入口距管體中心10B,管體中心距出口25B,上邊界距管體中心10B,下邊界距管體中心10B。計算域左側設置為速度入口,右側邊界為壓力出口,上下側邊界為壁面,隧道管體表面為靜止壁面。模型示意圖見圖1。計算模型采用四邊形網格劃分,斷面邊界層網格加密,近壁面第一層網格密度y+值小于1。y+值是無量綱化的壁面距離,使第一層網格y+值小于1可以使得第一層網格劃分到黏性底層,化解了網格數量與大小的矛盾,避免了第一層網格高度過大造成的數值誤差與不穩(wěn)定。動量方程及湍流方程采用二階迎風離散格式,求解采用SIMPLEC方法[24]。為進行懸浮隧道寬幅管體繞流特性與縮尺效應的初步探討,不考慮管體本身運動的影響。描述懸浮隧道斷面繞流的重要參數雷諾數Re和Strouhal數St為

    流體經過寬幅懸浮隧道截面時與管體之間產生的相互作用力包含3個部分:與來流方向垂直的分力為升力[26],來流方向的分力為阻力,合力作用點與截面形心不一致而對截面產生的扭轉力矩。一般采用無量綱的三分力系數來描述截面在流體作用過程中的受力特征,其中升力系數、阻力系數、力矩系數分別定義為

    式中:FL、FD和FM表示斷面升力、阻力和扭矩;B為斷面寬度。三分力系數在穩(wěn)定之后,升力系數和扭矩系數一般是平均值為0的周期曲線,阻力系數一般是平均值為某個非零值的周期曲線。因此為了進行更精確地描述,用CL,AM表示升力系數幅值,CM,AM表示扭矩系數幅值,CD表示阻力系數均值。

    1.2 計算模型選擇

    Fluent軟件的多種湍流模型中,SST k-ω模型在流動領域具有更高的精度和可信度,可以很好地模擬旋流、強壓力梯度的邊界層流動、分離流等,不包括復雜非線性黏性衰減函數,可較精確地模擬計算物體邊界層,適用于外部復雜流場計算,其中湍流動能k和比耗散率ω的輸運方程為[27]

    t(ρk)+xi(ρkui)=xj(Γkkxj)+Gk-Yk(3)

    式中:xi、xj分別為平面上2個坐標方向;ui、uj分別為流體沿2個方向的速度分量;Gk為湍流動能的速度梯度;Gω為比耗散率的速度梯度;Γk和Γω分別表示k和ω的有效擴散率;Yk和Yω表示湍流耗散項;Dω為交叉擴散項。

    1.3 截面形式選擇與計算工況

    考慮橢圓形、圓端形、八邊形3種按雙向八車道設計的寬幅懸浮隧道管體截面[25],并以圓形截面作為對照,各種截面形式與尺寸如圖2所示。初始入口的流體速度U取0.5~3.0m/s范圍內的多組值,用以計算三分力系數??紤]1∶10和1∶100兩種縮尺比,截面尺寸按照幾何相似準則考慮。Froude相似準則是懸浮隧道與洋流耦合作用時非常重要的相似準則,研究者們在試驗研究時[24-25]一般采用Froude相似準則處理縮尺后的洋流速度來保證試驗流場與實際流場相似,本研究中僅研究通過Froude相似準則處理過洋流速度后的試驗流場與實際流場管體繞流特性的差別。采用Froude相似準則對速度進行處理,有

    2 數值模擬結果與分析

    2.1 模型驗證

    為確保CFD數值模擬結果的有效性,消除不同參數對計算結果的不利影響,對數值模型的計算步長Δt與網格尺寸d進行無關性驗證。以直徑1.35m圓形截面超臨界狀態(tài)(

    Re=3.5×106)的阻力系數時程曲線為對象,采用不同時間步長得到的阻力系數如圖3(a)所示,可以看到當Δt減小至5×

    10-3s時已經收斂,此時阻力系數曲線與Δt=1×10-3、5×10-4 s時非常吻合。

    保持其他參數不變,選用不同的截面邊界層第一層網格尺寸d,計算所得結果如圖3(b)所示。可知當d加密至1.5×10-5m時,阻力系數時程曲線隨d的變化已經很小,此時不同d值對應的時程曲線基本吻合,其均值是一致的。綜合考慮計算效率與計算精度的要求,多車道寬幅懸浮隧道管體的數值模擬取值選擇Δt=5×10-3s,d=1.5×10-5m。另外考慮到計算精度,選用4套不同的網格數量進行無關性驗證。網格數量為10512、85751、153042、305212。當網格數量為10512時計算結果出入較大,不予考慮。當網格數量分別為85751、153042、305212時,阻力系數均值分別為0.56、0.57、0.55。計算所得結果比較吻合上文中結果。網格數量為153042的網格模型如圖4所示。

    為進一步驗證本研究方法的數值模擬結果,選用3組Re數的圓形截面計算結果與已有文獻進行對比驗證。當Re=2×104時,阻力系數均值為1.27,與FRANKE等[28]、詹昊等[29]的計算結果1.25、1.20相差不超過5%;當Re=1.4×105時,阻力系數均值為1.20,與CANTWELL[30]、SCHEWE[31]的計算結果1.23、1.18幾乎相等;當Re=3.5×106時,阻力系數均值為0.54,與程友良等[32]、劉二朋等[33]、SCHLICHING[34]的計算結果0.56、0.55、0.56高度吻合。證明本研究計算方法的精確性及計算模型的有效性。

    2.2 升力系數分析

    由于懸浮隧道斷面的存在改變了流場的分布與特征,在懸浮隧道管體表面流動的流體速度不同,從而導致斷面上下表面產生壓強差繼而產生升力。不同懸浮隧道截面形式的受力特性存在差異,選取橢圓形、圓端形和八邊形3種寬幅截面形式,并以圓形截面作為參照,圖5給出了這些截面形式在不同縮尺比λ情況下的升力系數幅值CL,AM隨流速U1的變化圖。由圖5可知,不同λ對應的圓形截面CL,AM均大于寬幅截面CL,AM,以對比圓形截面和八邊形截面的CL,AM最大值為例,λ=0.01(U1=0.05m/s),λ=0.1(U1=0.16 m/s)與λ=1(U1=0.5m/s)時圓形截面CL,AM最大值分別為八邊形截面的7.81倍、7.06倍與1.61倍。在3種寬幅懸浮隧道截面中,不同λ值對應CL,AM的大小均為:八邊形gt;圓端形gt;橢圓形,此外,λ=0.01、0.1時3種截面形式對應的升力系數幅值接近,λ=1時3種截面形式對應的升力系數幅值則差異明顯,因此3種縮尺比之間的尺寸效應明顯,值得關注。另外,在流速U1變化過程中,3種寬幅截面升力系數的穩(wěn)定性較好,即不同λ值時

    CL,AM

    隨U1變化而變化的程度很小,相對而言圓形截面對U1變化較為敏感,如λ=0.1時圓形截面CL,AM隨流速增大存在明顯下降。因此,圓形截面

    CL,AM大于寬幅形截面;3種寬幅截面中八邊形CL,AM最大、橢圓形最小,而各類截面CL,AM尺寸效應明顯;3種寬幅截面流速變化過程中CL,AM的穩(wěn)定性較好。

    為了更進一步描述一定流速作用下管體升力系數CL隨時間變化的過程,圖6繪出了3種縮尺情況下最大流速對應的升力系數時程曲線。由圖6可知,升力系數曲線均為均值為0的周期曲線,即升力以周期力的形式豎向作用于管體,是激發(fā)管體各類運動的主要荷載來源。同時,對比λ=0.01、0.1、1 (對應流速0.3、0.95、3.0m/s)下升力系數時程曲線,可獲取各類截面的尺寸效應規(guī)律。圓形截面CL,AM隨著λ的增大而減小,λ=1時的CL,AM約為λ=0.01時的55%;橢圓形截面CL,AM在λ=0.01時最大,λ=0.1、1時約為λ=0.01時的50 %;圓端形截面CL,AM隨著λ的增大而增大,在λ=1時的最大值約為λ=0.01時最小值的2.3倍;八邊形截面CL,AM在λ=0.1時最小,且最大值(λ=1)約為其4.3倍。因此,采用縮尺法對寬幅截面進行數值模擬和模型試驗時,升力系數隨著縮尺比的不同出現一定的變化,由此產生的差別需要考慮在內。

    對圖6中的升力系數時程曲線進行傅里葉變換之后可以得到懸浮隧道截面的St數。St數是描述脫渦頻率的重要參數之一,通過傅里葉變換所得St數分別為:圓形截面在λ=0.01、0.1、1時St數分別為0.22、0.215、0.255。橢圓形截面分別為0.38、0.42、0.53,圓端形截面分別為0.28、0.27、0.26,八邊形截面分別為0.20、0.18、0.18。計算結果顯示橢圓形截面St數最大、八邊形截面最小。比較某一種截面不同縮尺比之間的St數發(fā)現,隨著λ的增大,橢圓形、圓形截面的St數增大,其中橢圓形截面λ=1時的St數約為λ=0.01時的1.39倍,圓形截面λ=1時的St數約為λ=0.01時的1.16倍,而圓端形、八邊形截面的St數則隨著λ增大有少許減少。

    2.3 阻力系數分析

    處于流場中的懸浮隧道管體前后表面的壓強差的面積分就是結構所受的阻力。圖7給出了橢圓形、圓端形和八邊形3種寬幅截面以及圓形截面在不同的縮尺比λ情況下的阻力系數均值

    CD隨流速U1的變化圖。由圖7可知,在寬幅截面中3組λ值對應CD的大小均為:八邊形gt;圓端形gt;橢圓形,3種截面形式對應的阻力系數均值在λ=0.01時較接近,λ=1時差異明顯,λ=0.1時居于二者之間,3種縮尺比之間的尺寸效應明顯。此外,在流速U1變化過程中,

    3種寬幅截面阻力系數的穩(wěn)定性較好,即不同λ值時CD

    隨U1變化而變化的程度很小。相對而言,圓形截面在3種λ下的

    CD變化明顯且尺寸效應顯著:λ=0.01時圓形截面的

    CD整體上大于3種寬幅截面,且在U1變化過程中有少許下降;λ=0.1時圓形截面的CD一直大于橢圓形、圓端形截面,但隨著U1的增大CD迅速降至小于八邊形截面,減小幅度最大達到45%;λ=1時圓形截面的CD整體上大于橢圓形、圓端形截面而小于八邊形截面,且在U1變化過程中減幅很小。

    圖8給出了3種縮尺比分別對應流速0.3、0.95、3.0m/s時各類截面的阻力系數時程曲線。

    由圖8可知阻力系數曲線為均值(對應圖7中的

    C

    D)非零的周期曲線,其振幅隨著截面形式和縮尺比大小而變。

    同時,對比λ=0.01、0.1、1不同縮尺比下阻力系數時程曲線可發(fā)現:圓形與圓端形、橢圓形截面的阻力系數均隨λ的增大而下降,3種截面在λ=1時的阻力系數均值相對于λ=0.01時的值下降幅度分別為61%、65%和54%;相對而言,八邊形截面在λ=1時阻力系數均值最大,3種縮尺比對應的CD值依次為0.68、0.65、0.84。因此,不同λ下的阻力系數之間存在明顯的差別,在進行數值研究和模型試驗時需考慮阻力系數尺寸效應帶來的影響。

    2.4 力矩系數分析

    圖9給出了3種縮尺比下橢圓形、圓端形、八邊形及圓形截面的力矩系數幅值CM,AM隨流速U1的變化圖。

    由圖9可知,各類截面在不同λ下的CM,AM隨流速改變時的變化幅度較小,即不同流速作用下某類截面在固定縮尺比下的扭轉力矩差別較小。此外,八邊形截面的CM,AM整體上明顯大于其余3種截面,其中圓形截面和橢圓形截面的力矩系數處于較小水平。對比3種λ對應的

    CM,AM

    變化情況發(fā)現,不同截面類型的CM,AM尺寸效應各異:圓形截面與橢圓形截面的CM,AM

    在不同縮尺比時其值均很小,對應的尺寸效應可忽略;圓端形截面的

    CM,AM隨著λ的增大有小幅上升;八邊形截面的CM,AM在不同λ值時差異明顯,對應的尺寸效應顯著。

    對于懸浮隧道的寬幅截面,力矩系數增大是相對于圓形截面的一個重要特點,其本質是由于“寬幅”的結構特點導致升力、阻力的合力作用點與寬幅斷面形心不一致而產生扭矩。八邊形截面由于上、下不對稱,這種趨勢更加明顯,導致CM,AM明顯大于其他類型的截面。對于懸浮隧道管體而言,為滿足寬幅截面的通行需求,宜采用對稱設置截面等措施盡量減少扭矩,以避免可能引起的結構傾覆,同時配合上下游空間錨索等結構措施減小管體扭矩的不利影響。

    圖10給出了λ=0.01、0.1、1(對應流速0.3、 0.95、3.0m/s)時各類截面的力矩系數時程曲線。

    由圖10可知,圓形、橢圓形和圓端形截面的力矩系數曲線為零均值周期曲線,八邊形截面力矩系數曲線為非零均值周期曲線。對比不同截面形式和縮尺比的力矩系數曲線可發(fā)現:圓形截面的力矩系數非常小,當λ增大時力矩系數增加,λ=1時CM,AM約為0.005;橢圓截面的力矩系數比另外2種寬幅截面(圓端形、八邊形)小,且隨著λ增大力矩系數稍有減?。粓A端形截面的力矩系數隨著λ的增大而增大,其中λ=1時CM,AM約為λ=0.01時的1.5倍;八邊形截面的力矩系數在λ=1時最大而在λ=0.1時最小,且λ=1時CM,AM約為λ=0.1時的2.2倍。整體上,八邊形截面相對于其他截面所受傾覆力矩較大,各類截面在不同λ下力矩系數均存在不同程度的差異,采用縮尺比試驗和數值模擬時存在尺寸效應。

    2.5 渦量云圖對比

    圖11給出了典型參數(λ=0.01,U1=0.3m/s)下圓形、橢圓形和圓端形截面的三分力系數全過程時程曲線,在此基礎上進行局部放大并給出了一個穩(wěn)定周期內多個時間點(P1~P9)的渦量云圖。與之對應,圖12則給出了八邊形截面在3種縮尺比下λ=0.01、0.1、1(對應流速0.3、0.95、3.0m/s)對應的三分力系數時程曲線與一個周期內的渦量云圖。從圖11~圖12中穩(wěn)定之后的三分力系數可以看出,CL,AM、C

    D、CM,AM與圖5~圖10中系數分布規(guī)律一致,圓形截面與各類寬幅截面的卡門渦街規(guī)律十分明顯,各類截面的三分力系數不同時,渦的脫落位置與尾渦脫落的強度也存在差異。一般說來,繞流阻力包括摩擦阻力和壓差阻力[35],前者與黏性系數、繞流物幾何形狀有關,繞流物體沿流向尺寸越大摩擦阻力也越大。壓差阻力由邊界層分離與前后壓強差產生,分離點越靠后繞流體前后壓強差越小,壓差阻力也越小。由于懸浮隧道管體為鈍體截面,繞流阻力主要為壓差阻力,與尾流脫渦、渦量分離息息相關。三分力系數曲線經歷一段時間后呈現穩(wěn)定的周期性變化規(guī)律,升力系數變化一個周期即產生一個渦的脫落,同時產生2個周期的阻力系數變化。對比圖11(b)、圖11(d)、圖11(f)中3類截面與圖12(b)中八邊形截面在一個周期內的渦量云圖分布,當λ=0.01時圓形截面尾渦強度最大,分離點最靠前,橢圓形截面尾渦強度最小,分離點最靠后,八邊形、圓端形截面的尾渦強度介于圓形截面和橢圓形截面之間。因此,圓形截面的阻力系數最大,橢圓截面最小,這與圓形、八邊形、圓端形、橢圓形的阻力系數均值分別為1.21、0.68、0.48、0.23(圖7a)一致。這表明尾渦脫落對結構所受三分力有一定影響,觀察尾渦分布特點可以判斷三分力系數分布情況。

    由圖12所示的八邊形截面在3種λ下的時程曲線與渦量云圖可知,不同的λ對應的八邊形截面尾渦脫落非常明顯,當λ增大時尾渦強度降低且渦旋脫落沒有統一的規(guī)律。這說明不同縮尺比例下的尾渦脫落規(guī)律存在較大的差別,當使用縮尺方法研究八邊形結構的繞流特性時存在一定的尺寸效應,會帶來結果上的差值。

    3 結 論

    基于CFD軟件Fluent對多車道寬幅懸浮隧道截面在不同縮尺比下的繞流特性進行了數值模擬,得到了三分力系數、渦量云圖的規(guī)律,所得結論如下。

    1)橢圓形、圓端形和八邊形寬幅截面三分力系數的穩(wěn)定性較好,當流速增大時三分力系數變化不大。λ=0.1時,圓形截面阻力系數、升力系數隨著流速的增大明顯減小。

    2)圓形截面升力系數幅值CL,AM大于寬幅形截面;3種寬幅截面中,不同縮尺比對應CL,AM的大小均為八邊形>圓端形>橢圓形;橢圓形截面脫渦頻率最大即St數最大,八邊形截面脫渦頻率最小即St數最小。

    3)3種寬幅截面對應的阻力系數均值CD為八邊形>圓端形>橢圓形,寬幅截面形式在不同縮尺比之間阻力系數的尺寸效應明顯。

    4)八邊形截面力矩系數幅值CM,AM大于其余3種截面,且不同λ值時差異明顯,對應的尺寸效應顯著。其原因在于“寬幅”的特點導致升力、阻力的合力與形心不一致而產生扭矩,八邊形截面由于上、下不對稱而更加明顯。

    5)各類寬幅截面在不同縮尺比之間的三分力系數存在不同程度的差異,在進行數值研究和模型試驗時需考慮尺寸效應帶來的影響。

    6)圓形截面尾渦強度最大,分離點最靠前,對應阻力系數最大;橢圓形截面尾渦強度最小,分離點最靠后,對應阻力系數最?。话诉呅?、圓端形截面的尾渦強度介于這兩類截面之間。

    參考文獻:

    [1] STLID H.When is SFT competitive?[J].Procedia engineering,2010,4:3-11.

    [2] MARTIRE G,FAGGIANO B,MAZZOLANI F M.Compared cost evaluation among traditional versus innovative strait crossing solutions[J].Procedia engineering,2010,4:293-301.

    [3] 項貽強,陳政陽,楊贏.懸浮隧道動力響應分析方法及模擬的研究進展[J].中國公路學報,2017,30(1):69-76.

    XIANG Yiqiang,CHEN Zhengyang,YANG Ying.Research development of method and simulation for analyzing dynamic response of submerged floating tunnel[J].China journal of highway and transport,2017,30(1):69-76(in Chinese).

    [4] 丁浩,程亮,李科.懸浮隧道結構動力響應研究進展與展望[J].隧道建設(中英文),2019,39(6):901-912.

    DING Hao,CHENG Liang,LI Ke.Research progress and prospect on dynamic response of SFT structures[J].Tunnel construction,2019,39(6):901-912(in Chinese).

    [5] 周泰翔,賈軍波,鄧揚,等.懸浮隧道結構體系與動力響應研究進展[J].工程力學,2022,39(4):15-28.

    ZHOU Taixiang,JIA Junbo,DENG Yang,et al.Recent developments in the structural system and dynamic response of submerged floating tunnel[J].Engineering mechanics,2022,39(4):15-28(in Chinese).

    [6] 林巍,尹海卿,劉凌鋒,等.懸浮隧道實驗工程研究[J].中國港灣建設,2021,41(3):41-49.

    LIN Wei,YIN Haiqing,LIU Lingfeng,et al.Experimental engineering research on submerged floating tunnel[J].China harbour engineering,2021,41(3):41-49(in Chinese).

    [7] ZHANG H Q,YANG Z W,LI J Z,et al.A global review for the hydrodynamic response investigation method of submerged floating tunnels[J].Ocean engineering,2021,225:108825.

    [8] AADNESEN L,JAKOBSEN B,MOKSNES J.The case for floating submerged tunnels[J].Tunnels amp; tunnelling international,1999,31(6):32-34.

    [9] 范澤旭,袁勇,何任飛,等.固定支承式懸浮隧道在洋流渦激作用下的動力響應研究[J].鐵道科學與工程學報,2020,17(3):653-659.

    FAN Zexu,YUAN Yong,HE Renfei,et al.Dynamic response analysis of submerged floating tunnel supported on columns in vortex-induced vibration[J].Journal of railway science and engineering,2020,17(3):653-659(in Chinese).

    [10]葉騰超.定常流狀態(tài)下懸浮隧道結構外型與渦激響應理論研究[D].重慶:重慶交通大學,2018.

    [11]武玉濤,任華堂,夏建新.圓柱繞流研究進展及展望[J].水運工程,2017(2):19-26.

    WU Yutao,REN Huatang,XIA Jianxin.Research progress and prospects on flow across cylinder[J].Port amp; waterway engineering,2017(2):19-26(in Chinese).

    [12]KIM S,WILSON P A,CHEN Z M.Large-eddy simulation of the turbulent near wake behind a circular cylinder:Reynolds number effect[J].Applied ocean research,2015,49:1-8.

    [13]JIANG H Y,CHENG L.Strouhal-reynolds number relationship for flow past a circular cylinder[J].Journal of fluid mechanics,2017,832:170-188.

    [14]ALJURE D E,LEHMKHUL O,RODRí GUEZ I,et al.Three dimensionality in the wake of the flow around a circular cylinder at Reynolds number 5000[J].Computers amp; fluids,2017,147:102-118.

    [15]梁圣召,周軍偉,梅蕾.不同工況下有限長圓柱繞流實驗測量及分析[J].艦船科學技術,2021,43(11):19-24.

    LIANG Shengzhao,ZHOU Junwei,MEI Lei.Experiments and numerical simulations of the three-dimensional flow around a finite-length circular cylinder[J].Ship science and technology,2021,43(11):19-24(in Chinese).

    [16]陳松濤,趙偉文,萬德成,等.帶自由液面有限長圓柱繞流數值模擬[J].中國艦船研究,2022,17(1):91-98.

    CHEN Songtao,ZHAO Weiwen,WAN Decheng,et al.Numerical simulation of flows around a finite-length cylinder with free surface[J].Chinese journal of ship research,2022,17(1):91-98(in Chinese).

    [17]高洋洋,張演明,劉彩,等.不同雷諾數下傾斜圓柱繞流三維數值模擬研究[J].海洋工程,2020,38(1):86-100.

    GAO Yangyang,ZHANG Yanming,LIU Cai,et al.Three-dimensional numerical simulation on flow past an inclined circular cylinder at different Reynolds numbers[J].The ocean engineering,2020,38(1):86-100(in Chinese).

    [18]劉闖,黃劍峰,邵晨.基于LBM-LES的三維串列多圓柱繞流尾渦結構研究[J].應用力學學報,2022,39(1):186-194.

    LIU Chuang,HUANG Jianfeng,SHAO Chen.Research on wake structure of three-dimensional tandem multi-cylinder flow based on LBM-LES[J].Chinese journal of applied mechanics,2022,39(1):186-194(in Chinese).

    [19]李聰洲,張新曙,胡曉峰,等.高雷諾數下多柱繞流特性研究[J].力學學報,2018,50(2):233-243.

    LI Congzhou,ZHANG Xinshu,HU Xiaofeng,et al.The study of flow past multiple cylinders at high Reynolds numbers[J].Chinese journal of theoretical and applied mechanics,2018,50(2):233-243(in Chinese).

    [20]王長春.水中懸浮隧道與洋流耦合作用的模型試驗[D].成都:西南交通大學,2005.

    [21]秦銀剛,周曉軍.洋流作用下懸浮隧道動力學行為試驗研究[J].公路交通科技,2009,26(12):69-72.

    QIN Yingang,ZHOU Xiaojun.Experimental study of dynamic behavior of submerged floating tunnel under current effect[J].Journal of highway and transportation research and development,2009,26(12):69-72(in Chinese).

    [22]王廣地,周曉軍,高波.均勻流作用下懸浮隧道管段水動力參數分析[J].四川大學學報(工程科學版),2009,41(2):96-102.

    WANG Guangdi,ZHOU Xiaojun,GAO Bo.Analysis of hydrodynamic coefficients of the submerged floating tunnel subjected to homogeneous flow[J].Journal of Sichuan University(engineering science edition),2009,41(2):96-102(in Chinese).

    [23]郭曉玲,金瑞佳,劉名名,等.水流作用下懸浮隧道耦合動力數值模擬分析[J].水道港口,2020,41(3):275-283.

    GUO Xiaoling,JIN Ruijia,LIU Mingming,et al.Coupled numerical dynamic analysis of flow and submerged floating tunnel[J].Journal of waterway and harbor,2020,41(3):275-283(in Chinese).

    [24]曾繁旭,陳進,林巍.懸浮隧道特征橫斷面過水特性研究4:繞流場特征[J].中國港灣建設,2020,40(2):76-81.

    ZENG Fanxu,CHEN Jin,LIN Wei.Submerged floating tunnel conceptualized cross-section hydrodynamic characteristic study 4:flow field characteristics[J].China harbour engineering,2020,40(2):76-81(in Chinese).

    [25]李勤熙.波--流耦合作用下懸浮隧道關鍵參數數值仿真與物理試驗[D].重慶:重慶交通大學,2019.

    [26]周云龍,郭婷婷.高等流體力學[M].北京:中國電力出版社,2008.

    [27]MENTER F R.Two-equation eddy-viscosity turbulence models for engineering applications[J].AIAA journal,1994,32(8):1598-1605.

    [28]FRANKE J,FRANK W.Large eddy simulation of the flow past a circular cylinder at ReD=3900[J].Journal of wind engineering and industrial aerodynamics,2002,90(10):1191-1206.

    [29]詹昊,李萬平,方秦漢,等.不同雷諾數下圓柱繞流仿真計算[J].武漢理工大學學報,2008,30(12):129-132.

    ZHAN Hao,LI Wanping,FANG Qinhan,et al.Numerical simulation of the flow around a circular cylinder at varies Reynolds number[J].Journal of Wuhan University of Technology,2008,30(12):129-132(in Chinese).

    [30]CANTWELL B J.A flying hot wire study of the turbulent near wake of a circular cylinder at a Reynolds number of 140000 show affiliations[D].Pasadena:California Institute of Technology,1976.

    [31]SCHEWE G.On the force fluctuations acting on a circular cylinder in crossflow from subcritical up to transcritical Reynolds numbers[J].Journal of fluid mechanics,1983,133:265-285.

    [32]程友良,汪輝,趙洪嵩.超臨界雷諾數下旋轉圓柱繞流的三維大渦[J].可再生能源,2017,35(7):1094-1100.

    CHENG Youliang,WANG Hui,ZHAO Hongsong.Three-dimensional large eddy simulation of flow around a rotating circular cylinder under supercritical Reynolds number[J].Renewable energy resources,2017,35(7):1094-1100(in Chinese).

    [33]劉二朋,陳威,林永水,等.不同雷諾數下二維橢圓柱繞流的數值模擬研究[J].應用力學學報,2021,38(5):2025-2031.

    LIU Erpeng,CHEN Wei,LIN Yongshui,et al.Numerical simulation of flow around a two-dimensional elliptical cylinder with different Reynolds numbers[J].Chinese journal of applied mechanics,2021,38(5):2025-2031(in Chinese).

    [34]SCHLICHING H.Boundary-layer theory[M].New York:MC Graw-Hill,1979.

    [35]王福軍.計算流體動力學分析[M].北京:清華大學出版社,2004.

    (編輯 李坤璐)

    麻豆国产av国片精品| 亚洲成国产人片在线观看| 少妇熟女aⅴ在线视频| 亚洲aⅴ乱码一区二区在线播放 | 一a级毛片在线观看| www.熟女人妻精品国产| 大码成人一级视频| 99精品在免费线老司机午夜| 免费在线观看视频国产中文字幕亚洲| 动漫黄色视频在线观看| 久久香蕉精品热| 99国产极品粉嫩在线观看| 男女下面进入的视频免费午夜 | 在线视频色国产色| 欧美黄色片欧美黄色片| 51午夜福利影视在线观看| 国产精品久久久久久人妻精品电影| 久久精品91无色码中文字幕| 亚洲人成伊人成综合网2020| 久久久国产成人精品二区| 午夜成年电影在线免费观看| 看黄色毛片网站| 女性被躁到高潮视频| 又紧又爽又黄一区二区| 亚洲精品一区av在线观看| 日韩 欧美 亚洲 中文字幕| 宅男免费午夜| 成在线人永久免费视频| 亚洲色图 男人天堂 中文字幕| 老司机午夜十八禁免费视频| 久久人妻av系列| 久久久久久久久中文| 99国产精品一区二区三区| 美女国产高潮福利片在线看| 日本免费一区二区三区高清不卡 | 色综合亚洲欧美另类图片| 精品一品国产午夜福利视频| 老汉色av国产亚洲站长工具| 亚洲七黄色美女视频| 99riav亚洲国产免费| 午夜老司机福利片| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩高清在线视频| 久久久久国产精品人妻aⅴ院| 乱人伦中国视频| 国产一区二区激情短视频| 好男人在线观看高清免费视频 | 女性被躁到高潮视频| 国产精品久久久av美女十八| 两个人视频免费观看高清| 这个男人来自地球电影免费观看| 成年女人毛片免费观看观看9| 欧美大码av| 真人做人爱边吃奶动态| 一级黄色大片毛片| 国产亚洲精品第一综合不卡| 欧美老熟妇乱子伦牲交| 免费高清视频大片| cao死你这个sao货| 精品国产一区二区久久| 亚洲av电影不卡..在线观看| 9191精品国产免费久久| 亚洲国产高清在线一区二区三 | 一级毛片精品| 我的亚洲天堂| 亚洲欧美日韩高清在线视频| 女同久久另类99精品国产91| 国产成人系列免费观看| 88av欧美| 欧美在线一区亚洲| 国产精品美女特级片免费视频播放器 | 老司机靠b影院| 国产一区二区三区视频了| 在线视频色国产色| 久久久水蜜桃国产精品网| 满18在线观看网站| 操出白浆在线播放| 国产aⅴ精品一区二区三区波| 夜夜躁狠狠躁天天躁| 成年女人毛片免费观看观看9| 一进一出抽搐gif免费好疼| 日日爽夜夜爽网站| 色播在线永久视频| 正在播放国产对白刺激| 色综合欧美亚洲国产小说| 中亚洲国语对白在线视频| 免费av毛片视频| 亚洲美女黄片视频| 一边摸一边做爽爽视频免费| 精品卡一卡二卡四卡免费| 精品久久久久久成人av| 在线观看免费午夜福利视频| 最近最新免费中文字幕在线| 搡老妇女老女人老熟妇| 久久精品影院6| 国产xxxxx性猛交| 天天躁夜夜躁狠狠躁躁| 天堂动漫精品| 欧美成人性av电影在线观看| 亚洲成a人片在线一区二区| 亚洲天堂国产精品一区在线| 曰老女人黄片| 无遮挡黄片免费观看| 久久国产精品男人的天堂亚洲| 久9热在线精品视频| 欧美一区二区精品小视频在线| 一级毛片精品| 黑丝袜美女国产一区| 又黄又爽又免费观看的视频| 免费一级毛片在线播放高清视频 | 99riav亚洲国产免费| 18美女黄网站色大片免费观看| 国产精品日韩av在线免费观看 | 精品一品国产午夜福利视频| 成在线人永久免费视频| 国产麻豆69| 亚洲在线自拍视频| 精品久久久久久久人妻蜜臀av | 欧美乱色亚洲激情| 韩国精品一区二区三区| 久久久水蜜桃国产精品网| a级毛片在线看网站| 极品人妻少妇av视频| 欧美一级a爱片免费观看看 | 欧美日韩福利视频一区二区| 欧美大码av| 亚洲av五月六月丁香网| 男人舔女人的私密视频| 欧美日本中文国产一区发布| 亚洲五月婷婷丁香| 一区在线观看完整版| 99香蕉大伊视频| 亚洲精品粉嫩美女一区| 国产av一区在线观看免费| 国内精品久久久久久久电影| 中文字幕人妻丝袜一区二区| 黄色女人牲交| 国产黄a三级三级三级人| 叶爱在线成人免费视频播放| 一级,二级,三级黄色视频| 精品一区二区三区四区五区乱码| 精品久久久久久成人av| 女人精品久久久久毛片| 久久久久久亚洲精品国产蜜桃av| 亚洲欧洲精品一区二区精品久久久| 精品国产国语对白av| 国产亚洲精品一区二区www| 老司机深夜福利视频在线观看| 亚洲最大成人中文| 国产精品一区二区精品视频观看| 欧美日韩黄片免| 久久久久亚洲av毛片大全| 日本 av在线| 国产三级黄色录像| 欧美日韩乱码在线| 一区二区三区精品91| 桃色一区二区三区在线观看| 最新美女视频免费是黄的| 99国产精品99久久久久| 无限看片的www在线观看| 纯流量卡能插随身wifi吗| 亚洲成人国产一区在线观看| 中亚洲国语对白在线视频| 国产97色在线日韩免费| 久久国产精品人妻蜜桃| 又黄又粗又硬又大视频| 亚洲电影在线观看av| 人人妻,人人澡人人爽秒播| 在线观看一区二区三区| 99热只有精品国产| 亚洲精品中文字幕一二三四区| 两性夫妻黄色片| 午夜免费鲁丝| 麻豆久久精品国产亚洲av| 精品久久久久久,| 国产精品久久久人人做人人爽| 91麻豆精品激情在线观看国产| 免费av毛片视频| 久久久久久久午夜电影| 大陆偷拍与自拍| 久久精品成人免费网站| 国产亚洲精品av在线| 男女做爰动态图高潮gif福利片 | 国产aⅴ精品一区二区三区波| 日本免费一区二区三区高清不卡 | 久久精品91蜜桃| 亚洲狠狠婷婷综合久久图片| 国产成人啪精品午夜网站| 亚洲国产高清在线一区二区三 | 1024视频免费在线观看| 精品高清国产在线一区| 免费在线观看完整版高清| 精品国产一区二区久久| 日韩精品中文字幕看吧| videosex国产| 最新美女视频免费是黄的| 亚洲伊人色综图| 一个人观看的视频www高清免费观看 | 激情视频va一区二区三区| 久久精品成人免费网站| 中文字幕人妻熟女乱码| 宅男免费午夜| 国产色视频综合| 国产成人精品久久二区二区免费| 97碰自拍视频| 精品免费久久久久久久清纯| 法律面前人人平等表现在哪些方面| 9191精品国产免费久久| 午夜福利高清视频| 在线天堂中文资源库| 国产熟女午夜一区二区三区| 国产亚洲精品综合一区在线观看 | 国产成人av激情在线播放| 国语自产精品视频在线第100页| 1024香蕉在线观看| 亚洲五月天丁香| 精品第一国产精品| 黄色 视频免费看| 91麻豆精品激情在线观看国产| 一区在线观看完整版| 国产成+人综合+亚洲专区| 99re在线观看精品视频| 亚洲精华国产精华精| 国产精品一区二区在线不卡| av免费在线观看网站| 一a级毛片在线观看| 欧美成人性av电影在线观看| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品一区二区www| 日本 av在线| 两性夫妻黄色片| 99国产精品99久久久久| 成人国产一区最新在线观看| 成人亚洲精品一区在线观看| 人人妻人人澡人人看| 99国产精品99久久久久| 一进一出抽搐gif免费好疼| 欧美老熟妇乱子伦牲交| cao死你这个sao货| 女人精品久久久久毛片| 亚洲色图综合在线观看| 妹子高潮喷水视频| 91成年电影在线观看| 久久伊人香网站| 黑人巨大精品欧美一区二区蜜桃| 麻豆av在线久日| 国产亚洲精品一区二区www| 99精品欧美一区二区三区四区| 女人被狂操c到高潮| 妹子高潮喷水视频| 久久精品91无色码中文字幕| 91九色精品人成在线观看| 欧美老熟妇乱子伦牲交| 久久伊人香网站| 国产亚洲精品第一综合不卡| av福利片在线| av免费在线观看网站| av在线天堂中文字幕| 亚洲情色 制服丝袜| 国产成人系列免费观看| 国产亚洲欧美98| 亚洲熟女毛片儿| 免费观看人在逋| 精品久久久久久久久久免费视频| 亚洲国产高清在线一区二区三 | 九色国产91popny在线| 亚洲欧洲精品一区二区精品久久久| 欧美性长视频在线观看| 国产精品影院久久| 免费高清在线观看日韩| 操美女的视频在线观看| 我的亚洲天堂| 久久精品91蜜桃| 丰满的人妻完整版| 满18在线观看网站| 亚洲第一电影网av| 国产亚洲av高清不卡| 男女午夜视频在线观看| 99国产综合亚洲精品| 成年女人毛片免费观看观看9| 日韩精品中文字幕看吧| 宅男免费午夜| 亚洲国产欧美一区二区综合| 不卡av一区二区三区| 欧美激情极品国产一区二区三区| 69精品国产乱码久久久| 高潮久久久久久久久久久不卡| 国产精品1区2区在线观看.| 婷婷精品国产亚洲av在线| 久久人妻福利社区极品人妻图片| 十分钟在线观看高清视频www| 欧美日韩瑟瑟在线播放| 欧美在线黄色| 最近最新免费中文字幕在线| 看免费av毛片| 男女床上黄色一级片免费看| 国产在线精品亚洲第一网站| 国产欧美日韩一区二区三| 老司机福利观看| 男女午夜视频在线观看| e午夜精品久久久久久久| 亚洲成人久久性| 在线观看午夜福利视频| www日本在线高清视频| 国产成人精品在线电影| 中文字幕av电影在线播放| 亚洲av美国av| 91大片在线观看| 婷婷六月久久综合丁香| 精品久久久久久成人av| 一区福利在线观看| 午夜成年电影在线免费观看| 精品人妻在线不人妻| 女人被狂操c到高潮| 亚洲av成人不卡在线观看播放网| 国产精品秋霞免费鲁丝片| 少妇粗大呻吟视频| 波多野结衣一区麻豆| 成人永久免费在线观看视频| 欧美黑人精品巨大| 日本 欧美在线| 欧美久久黑人一区二区| 国产一级毛片七仙女欲春2 | www.999成人在线观看| 日本一区二区免费在线视频| 久久精品亚洲精品国产色婷小说| 国产欧美日韩一区二区三| 成年女人毛片免费观看观看9| 亚洲伊人色综图| 99久久国产精品久久久| 黄网站色视频无遮挡免费观看| 国产精品二区激情视频| 99riav亚洲国产免费| 久久久久久久久免费视频了| 成人av一区二区三区在线看| 日本精品一区二区三区蜜桃| 精品欧美国产一区二区三| 免费一级毛片在线播放高清视频 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲专区字幕在线| 亚洲成av片中文字幕在线观看| 成人亚洲精品av一区二区| 一进一出抽搐gif免费好疼| 免费女性裸体啪啪无遮挡网站| 精品久久久久久久久久免费视频| 神马国产精品三级电影在线观看 | 精品一区二区三区av网在线观看| 国产99白浆流出| 国产一区二区三区在线臀色熟女| 欧美中文综合在线视频| 午夜精品国产一区二区电影| 搡老妇女老女人老熟妇| 亚洲精品国产一区二区精华液| 自线自在国产av| 免费女性裸体啪啪无遮挡网站| 人妻丰满熟妇av一区二区三区| 丁香六月欧美| 国产精品永久免费网站| 国产av在哪里看| 电影成人av| 精品久久久久久久人妻蜜臀av | 亚洲激情在线av| 中出人妻视频一区二区| 国产一区二区三区视频了| 久久久久久国产a免费观看| 日韩欧美在线二视频| 精品久久久久久,| 欧美色欧美亚洲另类二区 | 色综合站精品国产| 国产精品亚洲av一区麻豆| 日本黄色视频三级网站网址| 亚洲黑人精品在线| 黄色a级毛片大全视频| 日本 av在线| 国产单亲对白刺激| www.自偷自拍.com| 日韩有码中文字幕| 国产亚洲精品久久久久久毛片| 日日干狠狠操夜夜爽| 黑人欧美特级aaaaaa片| 美女免费视频网站| 女人被狂操c到高潮| 999久久久国产精品视频| 国产成人影院久久av| 99国产精品一区二区三区| 啦啦啦 在线观看视频| 久久久久精品国产欧美久久久| 亚洲电影在线观看av| 69av精品久久久久久| 精品高清国产在线一区| 午夜福利,免费看| 亚洲国产精品久久男人天堂| 村上凉子中文字幕在线| 免费在线观看黄色视频的| 不卡av一区二区三区| 国产一区二区三区在线臀色熟女| 国产成人av激情在线播放| 精品人妻在线不人妻| 成人国语在线视频| 男人舔女人的私密视频| 国产色视频综合| 国产精品一区二区精品视频观看| 国产亚洲精品久久久久5区| 怎么达到女性高潮| 亚洲视频免费观看视频| 亚洲欧美精品综合久久99| 伊人久久大香线蕉亚洲五| 久久久久九九精品影院| 成年版毛片免费区| 国产伦人伦偷精品视频| 午夜日韩欧美国产| 99在线人妻在线中文字幕| 久久午夜综合久久蜜桃| 亚洲第一av免费看| 亚洲成a人片在线一区二区| 午夜激情av网站| 99在线人妻在线中文字幕| 亚洲国产精品久久男人天堂| 亚洲国产精品999在线| 热99re8久久精品国产| 亚洲精品一卡2卡三卡4卡5卡| 久久亚洲真实| 欧美精品啪啪一区二区三区| 99国产综合亚洲精品| 一级毛片高清免费大全| 成人免费观看视频高清| 黄片播放在线免费| 亚洲视频免费观看视频| www.熟女人妻精品国产| 在线观看免费视频网站a站| 欧美一级毛片孕妇| 无遮挡黄片免费观看| 性欧美人与动物交配| 大型黄色视频在线免费观看| 精品熟女少妇八av免费久了| 搞女人的毛片| 麻豆一二三区av精品| 亚洲成人国产一区在线观看| 老司机福利观看| 涩涩av久久男人的天堂| 无人区码免费观看不卡| 亚洲中文字幕一区二区三区有码在线看 | 国产精品 国内视频| 少妇 在线观看| 亚洲五月婷婷丁香| 亚洲欧美日韩另类电影网站| 亚洲aⅴ乱码一区二区在线播放 | 午夜日韩欧美国产| 亚洲第一青青草原| 精品日产1卡2卡| 国产野战对白在线观看| 一进一出抽搐动态| 美女国产高潮福利片在线看| 黄色视频,在线免费观看| 国产aⅴ精品一区二区三区波| 丁香欧美五月| 日本vs欧美在线观看视频| 亚洲熟女毛片儿| 欧美日韩亚洲综合一区二区三区_| 精品久久久精品久久久| 成在线人永久免费视频| 国产精品1区2区在线观看.| 成年人黄色毛片网站| 国产精品久久久久久亚洲av鲁大| 热99re8久久精品国产| 国产一区在线观看成人免费| 精品第一国产精品| 黄色视频,在线免费观看| 老熟妇乱子伦视频在线观看| 好男人在线观看高清免费视频 | 在线永久观看黄色视频| 9色porny在线观看| 亚洲性夜色夜夜综合| 9热在线视频观看99| 亚洲成av人片免费观看| 国产精品久久久久久亚洲av鲁大| 免费久久久久久久精品成人欧美视频| 亚洲精品美女久久久久99蜜臀| 欧美日本视频| 午夜福利在线观看吧| 老司机靠b影院| 最新在线观看一区二区三区| 两个人免费观看高清视频| 久久青草综合色| 久久中文字幕人妻熟女| 国产精品1区2区在线观看.| 如日韩欧美国产精品一区二区三区| 国产精品影院久久| www.熟女人妻精品国产| 亚洲av日韩精品久久久久久密| 色婷婷久久久亚洲欧美| 国产亚洲av嫩草精品影院| 久久久久亚洲av毛片大全| 国产激情久久老熟女| 老汉色av国产亚洲站长工具| 精品国产国语对白av| 亚洲五月天丁香| 999久久久国产精品视频| 婷婷丁香在线五月| 免费高清视频大片| 亚洲一区二区三区不卡视频| 又大又爽又粗| 日韩精品青青久久久久久| 免费久久久久久久精品成人欧美视频| 久久久水蜜桃国产精品网| 又黄又粗又硬又大视频| 99国产精品一区二区蜜桃av| 国产精品99久久99久久久不卡| 亚洲五月天丁香| e午夜精品久久久久久久| 黄色 视频免费看| 亚洲熟妇中文字幕五十中出| 日日夜夜操网爽| 免费高清在线观看日韩| 在线十欧美十亚洲十日本专区| 亚洲精华国产精华精| 国产一区二区三区视频了| 精品熟女少妇八av免费久了| 久久久久亚洲av毛片大全| 精品久久蜜臀av无| 在线观看免费视频网站a站| 国产片内射在线| 亚洲成人国产一区在线观看| 最新美女视频免费是黄的| 91精品三级在线观看| 男男h啪啪无遮挡| 久久中文字幕人妻熟女| 成人国产综合亚洲| 母亲3免费完整高清在线观看| 激情视频va一区二区三区| 女警被强在线播放| 乱人伦中国视频| 狂野欧美激情性xxxx| 国产不卡一卡二| 高清在线国产一区| 一二三四社区在线视频社区8| 国产人伦9x9x在线观看| 一进一出好大好爽视频| 黄色a级毛片大全视频| videosex国产| 国产成+人综合+亚洲专区| 免费人成视频x8x8入口观看| 精品不卡国产一区二区三区| 国产亚洲av嫩草精品影院| 伊人久久大香线蕉亚洲五| 少妇的丰满在线观看| 午夜a级毛片| 久久午夜综合久久蜜桃| 一二三四社区在线视频社区8| 一级毛片女人18水好多| 在线天堂中文资源库| 真人做人爱边吃奶动态| 极品教师在线免费播放| 一本大道久久a久久精品| av视频在线观看入口| 成人国产综合亚洲| 久久精品成人免费网站| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 人人妻人人澡欧美一区二区 | 国产高清激情床上av| 午夜日韩欧美国产| 男人的好看免费观看在线视频 | 性欧美人与动物交配| 岛国在线观看网站| 男女下面插进去视频免费观看| av免费在线观看网站| 亚洲五月婷婷丁香| 12—13女人毛片做爰片一| 99在线视频只有这里精品首页| 伊人久久大香线蕉亚洲五| 国产97色在线日韩免费| 精品人妻1区二区| 99久久综合精品五月天人人| 搡老岳熟女国产| 正在播放国产对白刺激| 欧美日本中文国产一区发布| 波多野结衣巨乳人妻| 国产亚洲精品久久久久5区| 欧美国产精品va在线观看不卡| 麻豆av在线久日| 国产99白浆流出| 国产单亲对白刺激| 18禁国产床啪视频网站| 亚洲精品av麻豆狂野| 国产亚洲精品久久久久5区| 国产精品98久久久久久宅男小说| 一本久久中文字幕| 999久久久精品免费观看国产| 在线观看免费午夜福利视频| 97碰自拍视频| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区| 久久国产精品人妻蜜桃| 国产精品日韩av在线免费观看 | 丰满人妻熟妇乱又伦精品不卡| av在线天堂中文字幕| 国产成人欧美| 欧美日韩亚洲国产一区二区在线观看| or卡值多少钱| 不卡一级毛片| 亚洲第一欧美日韩一区二区三区| 一区福利在线观看| 色av中文字幕| 岛国视频午夜一区免费看| 脱女人内裤的视频| 久久中文看片网| 免费久久久久久久精品成人欧美视频| 国产精品久久久久久精品电影 | 一区二区三区国产精品乱码| 香蕉丝袜av| 夜夜看夜夜爽夜夜摸| 欧美日韩亚洲国产一区二区在线观看| 嫩草影院精品99| 成人av一区二区三区在线看|