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

    入水角度對球體高速入水空泡特性影響研究

    2023-04-06 05:50:28李利劍張敏弟王占瑩黃彪劉思華譚鎮(zhèn)坤
    裝備環(huán)境工程 2023年3期
    關(guān)鍵詞:空泡球體射流

    李利劍,張敏弟,王占瑩,黃彪,劉思華,譚鎮(zhèn)坤

    (1.北京理工大學(xué) 機械與車輛學(xué)院,北京 100081;2.北京宇航系統(tǒng)工程研究所,北京 100076)

    入水過程是指物體以一定初速度從空氣穿過氣水交界面進入水中的過程,也就是從物體觸及水面直至進入水中且在水中穩(wěn)定運動到水底[1-2]的整個過程。入水現(xiàn)象在自然學(xué)科、工業(yè)生產(chǎn)、仿生技術(shù)、航空航天[3]以及海戰(zhàn)裝備和海洋開發(fā)[4-5]等領(lǐng)域都有廣泛的應(yīng)用。隨著人們對能源需求的日益增長,世界各國開始重視起海洋資源開發(fā),在海洋開發(fā)過程中,都會涉及到入水問題。如船錨入水、海洋鉆井采油平臺、水下生產(chǎn)設(shè)備吊裝入水等[6]。由于入水運動涉及固、液、氣的三相耦合作用,是一個具有強瞬態(tài)和非定常特性的過程,同時還伴隨著湍動、相變、可壓縮、介質(zhì)突變等復(fù)雜流動現(xiàn)象[7-10],極易造成物體運動過程中結(jié)構(gòu)的破壞及內(nèi)部組件損壞、失靈和運動失穩(wěn)等故障。因此,研究物體入水過程中空泡演化及其流體動力特性對于優(yōu)化水下生產(chǎn)設(shè)備的設(shè)計,提升海洋開發(fā)機械的入水品質(zhì),以及推動水動力學(xué)基礎(chǔ)理論的發(fā)展都具有非常重要的意義。

    早在19 世紀末期,Worthington 等[11-12]采用單閃照相機對物體入水過程進行了初次探索,獲取了第一張入水運動的照片,首次系統(tǒng)地描述了入水過程中產(chǎn)生的噴濺、生成的空泡和空泡閉合等流動現(xiàn)象,定性分析了入水空泡的演化過程及影響空泡生成的部分因素。之后40 年對入水現(xiàn)象的試驗研究都是基于Worthington 的研究方法進行模仿或者改進,而且他對入水現(xiàn)象的部分定義一直沿用至今。對于入水問題的研究主要集中在入水空泡發(fā)展過程研究、入水沖擊與水動力研究以及射彈入水穩(wěn)定性研究3 個方面。

    對于物體入水空泡演化方面,May 等[13-15]展開了大量具有軍事應(yīng)用背景的入水機理試驗研究,開展了球體幾何外形、傅汝德數(shù)等因素對球體高速入水過程中空泡形態(tài)及球體所受阻力的影響研究。Truscott 等[16]開展了密度和表面潤濕性對球體垂直入水時空泡的發(fā)展規(guī)律及流體動力特性影響的研究,并且基于圖像處理技術(shù)獲取了入水過程中球體的位移、速度等運動參數(shù),分析了流體動力特性。張孝石等[17]通過水洞試驗對通氣空泡的潰滅特性進行了研究,結(jié)果表明,在航行體停止通氣后,空泡經(jīng)歷了空泡凹陷、空泡斷裂、空泡脫落和潰滅4 個階段。空泡脫落時,在脫落空泡區(qū)域形成射流,航行體表面承受沖擊波壓力,隨后沖擊波壓力迅速衰減。李達欽等[18-19]采用試驗觀測和數(shù)值計算相結(jié)合的方法,對不同密度比的疏水性球體垂直入水空泡形態(tài)及水下流體動力特性進行了研究,此外還分析了旋轉(zhuǎn)球體入水過程中空泡的演化歷程,認識到角速度對空泡的噴濺、形態(tài)有顯著影響,但對空泡夾斷深度影響不大。

    對于物體入水沖擊載荷與彈道穩(wěn)定性方面,El Malki Alaou 等[20]同樣結(jié)合試驗觀測和數(shù)值模擬對不同錐角回轉(zhuǎn)體入水過程開展了研究,重點關(guān)注了不同頭部錐角回轉(zhuǎn)體入水時的沖擊載荷。黃鴻鑫等[21]通過數(shù)值模擬方法,從空泡形態(tài)、速度衰減以及俯仰角變化探討了射彈頭部形狀和質(zhì)心位置對高速入水穩(wěn)定性的影響。胡青青[22]探討了彈體高速時頭型和長細比對入水彈道穩(wěn)定性的影響,并分析了其對入水過程的影響。Moxnes 等[23]通過試驗分析和數(shù)值模擬計算總結(jié)了不同類型的射彈在水中發(fā)生跳彈的規(guī)律。陳誠等[24]開展了入水角為20°時帶圓盤空化器超空泡航行器傾斜入水試驗,針對入水沖擊載荷進行了研究。侯宇等[25]通過試驗研究了射彈小角度高速入水過程中的彈道特性,分析了彈體不同側(cè)滑角入水沖擊過程的彈道軌跡、噴濺演變和水下壓力波傳播特征,發(fā)現(xiàn)在高速入水過程中,增大射彈側(cè)滑角,會降低射彈入水彈道的穩(wěn)定性。周可等[26]采用數(shù)值模擬的方法對跨介質(zhì)航行器高速斜入水跳彈現(xiàn)象進行了研究,結(jié)果表明,航行器高速入水時,頭尾兩側(cè)受力不均,形成了向上的偏力矩,改變了航行器的俯仰角,航行體入水角度越小,越容易發(fā)生跳彈現(xiàn)象。2021 年,田北晨等[27]針對入水參數(shù)對跨介質(zhì)飛行器靜水情況下近水面處滑跳運動過程進行了數(shù)值模擬,結(jié)果表明,跨介質(zhì)飛行器在近水面處滑跳運動可分為穩(wěn)定滑跳階段及水面滑行階段。穩(wěn)定滑跳階段飛行器的運動軌跡呈類正弦規(guī)律波動,逐次入水砰擊瞬間,飛行器均會受到較大的砰擊載荷??缃橘|(zhì)飛行器近水面滑跳運動過程中,入水速度越大,砰擊載荷越大,穩(wěn)定滑跳次數(shù)越多,穩(wěn)定滑跳距離越遠。

    入水過程對于反潛武器能否成功入水并準(zhǔn)確打擊目標(biāo)起著至關(guān)重要的作用。目前對于射彈高速入水穩(wěn)定性的研究主要是從入水沖擊載荷和彈體頭型等方面進行,針對入水角度對彈道穩(wěn)定性的影響也開展了一定的研究。但是,在研究入水角度對彈道穩(wěn)定性影響的過程中,從入水空泡演化的角度進行的研究還不充分,而入水空泡的演化對反潛武器入水瞬間的流體動力特性和水下彈道的穩(wěn)定性具有重大的影響。因此,本文針對球體,結(jié)合試驗觀測與數(shù)值模擬討論了高速運動的球體以不同入水角度入水時空泡的演化過程,分析明確入水角度對空泡流動特性的影響。

    1 試驗系統(tǒng)及數(shù)值計算方法

    1.1 試驗系統(tǒng)

    基于高速攝像機搭建的球體入水觀測平臺如圖1所示,平臺由發(fā)射系統(tǒng)、水箱、圖像采集和處理系統(tǒng)組成。其中,發(fā)射系統(tǒng)主要由角度調(diào)節(jié)裝置和發(fā)射裝置組成。角度調(diào)節(jié)裝置由支撐板、支架上的定位孔和支撐軸組成,試驗中通過將支撐板固定到不同的定位孔實現(xiàn)球體以不同角度入水。發(fā)射裝置由空氣壓縮機、電磁閥、安全閥和發(fā)射管組成,試驗中打開安全閥,通過壓氣機提供高壓氣體膨脹做功提供驅(qū)動力,在電磁閥開啟瞬間,球體在高壓氣體的驅(qū)動下,經(jīng)由發(fā)射管中加速發(fā)射。通過調(diào)節(jié)壓氣機內(nèi)高壓氣體壓力和流量控制球體的速度。高壓氣缸內(nèi)氣體壓力與速度的關(guān)系[28]為:

    圖1 球體入水觀測平臺 Fig.1 Water entry observation platform of sphere

    式中:P1為壓氣機內(nèi)壓力;V為壓氣機容積;k為氣體絕熱指數(shù);L為發(fā)射管長度;A為發(fā)射管面積;P為大氣壓力;Mdan為球體質(zhì)量;φ為次要功系數(shù),用于能量損耗的修正,一般為1.02~1.1;v為球體運動速度。

    試驗水箱尺寸為0.6 m×0.6 m×0.6 m,四周及底面均為10 mm 厚的有機玻璃,在水箱底部墊有厚度為0.1 m 的緩沖層,防止球體高速入水時與水箱底部發(fā)生撞擊而損傷水箱。試驗球體直徑為D=5 mm,密度為7 850 kg/m3。圖像采集和處理系統(tǒng)由高速攝像機、光源和計算機組成。試驗中,將高速攝像機放置在試驗水箱正前方,采集頻率設(shè)為4 000 fps。采用多通道同步觸發(fā)器同步觸發(fā)電磁閥和高速攝像機,實現(xiàn)球體入水過程的實時記錄和實時顯示。光源布置為陰影法布局,記錄入水空泡演變過程形成陰影,進而分析提取空泡輪廓和流場流動結(jié)構(gòu)。

    1.2 數(shù)值計算方法

    在入水過程中,存在氣、液兩相界面的運動和變形。本文假定氣、液兩相為均相流動,采用VOF 方法通過容積跟蹤法對相界面進行追蹤??刂品匠虨椋?/p>

    式中:ρm為混合介質(zhì)的密度;μm為混合動力黏性系數(shù);p為靜壓強;g為重力加速度;μt為湍流黏性系數(shù);ui、uj分別為i、j方向的速度分量。對于ρm和μm,其表達式為:

    式中:α表示蒸汽體積率,表征水、汽兩相之間的相變;ρ為介質(zhì)的密度;下標(biāo)v、g、l 分別表示蒸汽、空氣和水。

    湍流模型選擇由Launder 和Spalding 提出的k-ε標(biāo)準(zhǔn)湍流模型,它把渦黏系數(shù)和湍動能k及湍動能耗散率ε聯(lián)系在一起,控制方程為:

    式中:Pt為湍動能生成項;σk、σε、Cε1、Cε2為常數(shù),根據(jù)經(jīng)驗[29],取值分別為1.0、1.3、1.44、1.92。

    球體在高速入水時會發(fā)生空化現(xiàn)象,基于質(zhì)量傳輸?shù)目栈P停浣M分傳輸方程為:

    式中:為單位體積質(zhì)量傳遞速率;+為單位體積液體質(zhì)量傳遞速率;為單位體積蒸汽質(zhì)量傳遞速率。

    Zwart-Gerber-Βelamri 空化模型結(jié)合泡間兩相流動理論,重點考慮空穴生長和潰滅時氣泡體積變化的影響,適用于空化的非定常過程。在該模型中,單位體積內(nèi)的相間傳輸速率為:

    式中:αnuc為汽核體積分數(shù);RΒ為汽泡半徑;pv為汽化壓強;Fvap和Fcond分別為蒸發(fā)和凝結(jié)經(jīng)驗系數(shù)。

    本文采用的計算域以及邊界條件設(shè)置如圖2 所示。設(shè)置球體直徑D=5 mm,與試驗所用球體直徑一致,計算域長、寬、高分別為143.4D、120D、60D,其中空氣域高43.4D,水域高100D,球體底部距離自由液面3D。球體周圍及自由液面處網(wǎng)格進行加密,分離求解算法采用SIMPL 算法,壓力項采用PRESTO格式,空間與時間離散采用一階迎風(fēng)格式。

    圖2 計算域網(wǎng)格及邊界條件設(shè)置 Fig.2 Settings of computational domain mesh and boundary condition

    為了驗證所采用的數(shù)值計算模型,對球體以30 m/s 垂直入水過程進行了數(shù)值模擬,并與試驗結(jié)果進行對比。數(shù)值模擬與試驗所得球體入水深度隨時間變化曲線如圖3 所示,數(shù)值模擬得到的球體位置與試驗測得的位置誤差在5%以內(nèi)。試驗觀測和數(shù)值模擬得到的空泡閉合輪廓如圖4 所示。由圖4 可見。數(shù)值模擬提取的空泡輪廓與試驗采集得到結(jié)果基本完全吻合。對比圖3 和圖4 結(jié)果可知,本文構(gòu)建數(shù)值模擬方法可靠可信,能夠用于后續(xù)研究。

    圖3 入水深度隨時間變化曲線 Fig.3 Curve of water entry depth over time

    圖4 試驗與數(shù)值模擬空泡閉合形態(tài)對比 Fig.4 Experimental and numerical simulation of cavity closure patterns

    2 結(jié)果與分析

    2.1 球體垂直入水空泡演化研究

    2.1.1 演化過程

    基于試驗觀測,本文根據(jù)球體入水過程中入水空泡的形態(tài)特征將入水過程分為入水沖擊、空泡形成、空泡閉合和空泡潰滅 4 個階段。速度v=30 m/s時球體入水空泡演化過程如圖5 所示。在t=0 ms 時刻,球體撞擊水面,與水體產(chǎn)生動量交換,使得水體排向周圍,在與球體的接觸面處發(fā)生流動分離,形成噴濺。隨著球體入水深度的增加,噴濺水膜的厚度也逐漸增大,進一步發(fā)展成對稱的水冠。同時,球體在運動過程中擠壓周圍水體,在水中形成空泡,空氣流入空泡內(nèi)部,空泡逐漸被拉長。在空泡閉合階段,隨著動能逐漸耗盡,液膜在重力、壓差阻力以及表面張力的作用下逐漸收縮閉合,從而切斷了外界空氣與空泡的聯(lián)系。當(dāng)球體繼續(xù)向下運動時,空泡受到球體的作用,脫離自由液面,進一步發(fā)生面閉合,并產(chǎn)生較強的回射流??张蓍]合時產(chǎn)生的回射流會沖擊空泡壁面,導(dǎo)致壁面不穩(wěn)定,同時因空泡尾部動能耗盡而不斷收縮,導(dǎo)致空泡夾斷潰滅。然而,空泡潰滅夾斷時也會產(chǎn)生回射流,進一步引起空泡的潰滅。

    圖5 球體入水空泡演化過程 Fig.5 Evolution process of sphere’s water entry cavity

    2.1.2 演化特性分析

    在入水沖擊階段,球體撞擊自由液面,與周圍流體接觸,并發(fā)生動量交換,周圍流體獲得動量,形成噴濺。入水沖擊階段和空泡形成階段的流場圖像和壓力分布云圖如圖6 和圖7 所示,球體入水時噴濺的受力分析如圖8 所示。由圖6 壓力場可知,球體與水面接觸部分壓力明顯高于周圍流場壓力,而在球體未接觸水面部分形成了對稱的低壓區(qū)域。這是由于高速運動的球體擾動周圍空氣,使周圍空氣流速增大,壓力降低。根據(jù)壓力分布云圖可知,水體在接觸球體位置壓力最大,高達0.638 MPa,離球體越遠,壓力越小。隨著入水深度的增大,水體與球體發(fā)生流動分離,形成水冠,由于球體周圍壓力較低,空氣流入空泡內(nèi)部,球體開始攜帶空泡向下運動,且入水空泡始終與大氣相連。由圖7 壓力與矢量分布可知,空泡內(nèi)部出現(xiàn)明顯的低壓區(qū),氣體高速流入空泡內(nèi)部,同時球體頭部的高壓區(qū)域也明顯減小,此時球體頭部的壓力為0.357 MPa,比入水沖擊時減小44.0%。圖8 中,Uwater為液膜頂端運動速度,對應(yīng)加速度為as,液膜所承受的主要受力包括自身重力Fgs、由于氣流高速運動所產(chǎn)生的壓差力FΔp、流動阻力Fd以及阻礙液膜運動的表面張力Fσ。取一液膜微元,質(zhì)量記為ms,可建立受力平衡方程如下:

    圖6 入水沖擊成階段試驗圖像與壓力分布云圖 Fig.6 Experimental image and pressure distribution pattern during water entry impact

    圖7 空泡形成階試驗圖像與壓力分布云圖 Fig.7 Experimental image and pressure distribution pattern during cavity formation

    圖8 噴濺受力分析 Fig.8 Analysis of splash force

    式中:流動阻力Fd與速度平方成正比,壓差力FΔp的定義為:

    式中:ΔAs為液膜微元的截面積;ps為噴濺內(nèi)壓強;p∞為環(huán)境壓強,p∞=101 325 Pa;ρg為氣體密度;Uair為當(dāng)?shù)貧饬魉俣?;dVc入水空泡體積增量成;Ds為空泡直徑。由以上分析可知,噴濺的形狀與運動與空氣流速Uair有關(guān),空氣流速Uair越大,壓差力FΔp越大,向內(nèi)收縮越快。

    由圖9 可知,在空泡閉合過程中,空泡內(nèi)部壓力較低,空泡周圍噴濺區(qū)域的外側(cè)壓力較高,在內(nèi)外壓差的作用下,空泡不斷向中間收縮,此時向空泡內(nèi)部運動的氣體的最大速度也降低到2 m/s,遠遠小于空 泡形成階段的氣體流速。進一步分析此時流場受力分布情況,并繪制了空泡形成階段出現(xiàn)噴濺時的受力示意圖(圖10)。圖10 中,Uwater為液膜頂端運動速度,對應(yīng)的加速度為ac,液膜所承受的主要受力包括重力Fgc、氣流高速運動所產(chǎn)生的壓差力FΔp、流動阻力Fd以及阻礙液膜運動的表面張力Fσ。取一液膜微元,質(zhì)量記為mc,可建立如下受力平衡方程:

    圖9 空泡閉合過程試驗圖像與壓力、速度矢量分布 Fig.9 Experimental image and distribution of pressure and velocity vectors during cavity closure

    圖10 空泡閉合階段受力分析 Fig.10 Force analysis during cavity closure

    式(14)中,壓差力FΔp的計算方法與空泡形成階段相同。由以上分析可知,空氣的流速在空泡閉合階段起著非常重要的作用,流入空泡內(nèi)部空氣的流速越大,壓差越大,空泡閉合越快。此外,空泡閉合還受到液膜速度的影響,液膜速度越快,受到的流動阻力Fd越大。液膜速度、空氣流速都與球體的運動有關(guān),因此球體的運動參數(shù)(如角度、速度)對空泡的閉合起著至關(guān)重要的作用。

    空泡潰滅階段的試驗圖像與壓力、速度分布云圖如圖11 所示。隨著球體繼續(xù)向下運動,空泡尾部水體向四周擴張的速度小于0.5 m/s,空泡尾部壁面動能耗盡,發(fā)生潰滅。由圖11 可見,球體頭部壓力仍明顯高于其他區(qū)域的壓力,但由于球體運動速度減小,球體頭部的壓力也發(fā)生了減小,此時球體頭部壓力為0.18 MPa,減小到入水沖擊時壓力的36%,空泡形成階段壓力的71.8%。另外,在空泡的收縮處壓力較高,而收縮段下方壓力明顯低于其他區(qū)域。進一步對收縮段流場進行受力分析,收縮段空泡形態(tài)特征如圖12 所示。在空泡脫離自由液面后,隨球體向下運動的空泡的不同區(qū)域會同時存在收縮和膨脹2 種狀態(tài)。造成這種現(xiàn)象的原因主要是,空泡閉合后不再有空氣進入空泡內(nèi)部,隨著球體向下運動,空泡在拉扯力的作用下會趨于細長,同時靜壓強隨著入水深度逐漸增大,空泡會在某一動能小于壓能的位置處迅速收縮,并與靜壓強平衡,此時收縮處空泡直徑變小,氣體只有流入其他空泡區(qū)域,且導(dǎo)致體積增加。

    圖11 空泡潰滅時的試驗圖像和局部壓力、速度分布云圖 Fig.11 Experimental image and local pressure and velocity distribution pattern during cavity collapse

    圖12 空泡收縮段形態(tài)特征 Fig.12 Morphological characteristics of cavity contraction segment

    2.2 入水角度對入水空泡演化的影響

    速度為30 m/s 的球體以不同角度入水的空泡演化過程如圖13 所示。圖13a 中,t=0 ms 時,球體先撞擊液面形成噴濺。t=0.5 ms 時,噴濺繼續(xù)擴展,形成了對稱的水冠,水冠在重力、表面張力以及內(nèi)外壓差的作用下開始收縮,直至噴濺發(fā)生閉合,阻斷了外部空氣的進入(t=0.75~1.25 ms)。隨著球體繼續(xù)向下運動,入水空泡也隨著球體向下運動,已經(jīng)切斷了和空氣連接的空泡受到向下拉扯的作用力,當(dāng)作用力大于空泡所受到的表面張力時,空泡就會從自由液面脫落,發(fā)生面閉合,并產(chǎn)生比較強的回射流(t=2.25 ms)?;厣淞鞯臎_擊會導(dǎo)致空泡壁面不穩(wěn)定,同時在靜水壓的作用下,空泡尾部開始收縮潰滅,此時空泡形狀呈梭形(t=2.25~4 ms)。在空泡潰滅夾斷時,也會產(chǎn)生回射流,從而進一步地引起入水空泡的潰滅(t=5~ 6 ms)。球體在斜射入水時,空泡演變過程與垂直入水相似,均存在入水沖擊、空泡形成、空泡閉合和空泡潰滅幾個階段,且都在空泡閉合時發(fā)生面閉合。但由于斜射入水時球體水平方向分速度的存在,導(dǎo)致其迎流面的噴濺較高、背流面的噴濺較低,進而影響了水冠的形態(tài)(如圖13 b、c 所示)。其次,在入水空泡發(fā)生面閉合時,與垂直入水也有所區(qū)別,入水角度為60°時,空泡在t=3.25 ms 發(fā)生了面閉合,閉合時產(chǎn)生了微弱回射流。球體入水角度為30°時,入水空泡在t=9.5 ms 時發(fā)生了面閉合,此時空泡閉合沒有產(chǎn)生回射流。隨著入水角度的減小,空泡閉合時間逐漸增加,且產(chǎn)生的回射流在逐漸減弱。在入水角度為30°時,空泡面閉合不再產(chǎn)生回射流。入水角度為60°和30°球體第一次潰滅夾斷的時間分別為t=7.5 ms、t=36.5 ms,潰滅時間隨著入水角度的減小而逐漸延遲。在空泡第一次潰滅時,都是以小尺度氣泡團的形式潰滅。球體入水角度為60°的入水空泡與垂直入水時的入水空泡一樣,第一次潰滅都產(chǎn)生了回射流,但以30°入水時,入水空泡第一次潰滅并未產(chǎn)生回射流。此時,入水角度為30°的入水空泡也與之前不同,空泡發(fā)生了抖動,不再呈現(xiàn)出梭形(圖 13c 中t=36.5 ms)。入水角度為30°的入水空泡在第二次潰滅時與第一次不同,第二次潰滅不再以小尺度氣泡團的方式潰滅,而是斷裂成多個尺度更大的氣腔從空泡上脫離。

    圖13 球體以不同入水角度入水空泡演化過程 Fig.13 Evolution process of sphere’s water entry cavity at different water entry angles

    由于球體的運動參數(shù)對空泡的閉合起著至關(guān)重要的作用,本文開展了球體以不同角度入水過程中空泡演化過程研究。球體入水位移、速度和加速度的時間變化曲線如圖14—16 所示。這些曲線是依據(jù)對多組連續(xù)試驗圖像進行圖像處理后,確定球體位置,進而結(jié)合記錄頻率,提取球體運動速度和加速度繪制而成的。對比圖14 中各條曲線可知,球體斜射入水時都發(fā)生了不同程度的偏離,并且偏離程度隨著入水角度增大而減小,入水角度為30°時最大,入水角度90°時,幾乎沒有發(fā)生偏離。圖15 中球體的速度總體上呈現(xiàn)減小的趨勢,但在空泡閉合以及潰滅時出現(xiàn)了不同程度的增加,并且增加幅度隨著入水角度的增大逐漸減小。在空泡潰滅階段,球體以90°和60°入水產(chǎn)生的空泡在潰滅夾斷時都產(chǎn)生了回射流,導(dǎo)致球體速度增加。空泡夾斷潰滅時,球體90°入水形成空泡的回射流強度更大,速度增大也更多。入水角度為90°時球體速度相比潰滅前增加了2.08 m/s,入水角度為60°時,球體速度在夾斷潰滅時相比潰滅前增加了1.403 m/s。

    圖14 球體以不同角度入水后坐標(biāo)變化曲線 Fig.14 Coordinates of sphere after entering water at different angles

    圖15 球體以不同角度入水后速度隨時間變化曲線 Fig.15 Velocity of sphere over time after entering water at different angles

    對比圖16 中各條曲線可知,所有工況下球體加速度的絕對值都在入水的瞬間達到最大,且隨著角度的增大,加速度峰值也在增大。不同入水角度下球體 入水5 ms 內(nèi)加速度的波動范圍見表1。由表1 可知,隨著入水角度的增大,加速度波動越來越劇烈。加速度最大絕對值(即球體入水沖擊時刻加速度的絕對值)隨角度的變化曲線如圖17 所示。由圖17 可知,30°時加速度是最大絕對值為7 039 m/s2,60°時加速度的最大絕對值為14 215 m/s2,90°時為15 166 m/s2。對比各值可知,加速度最大絕對值隨入水角度的增大而增大,但其增長速度隨角度的增大而減小。

    圖16 球體以不同角度入水后加速度隨時間變化曲線 Fig.16 Acceleration curve of sphere after entering water at different angles

    圖17 球體加速度最大絕對值隨入水角度變化曲線 Fig.17 Maximum absolute acceleration curve of sphere at different water entry angles

    表1 不同入水角度球體入水5 ms 后的加速度波動的范圍 Tab.1 Range of acceleration fluctuations of sphere 5 ms after entering water at different angles

    2.2.1 入水角度對入水沖擊和空泡形成的影響

    球體以不同角度入水過程中,流場圖像與壓力分布云圖如圖18 所示。對比圖18 中各圖可知,隨著入水角度的減小,球體頭部沾濕處壓力逐漸減小,沾濕區(qū)域的最大壓力也隨著入水角度的減小而減小。球體以90°入水時,頭部壓力為0.638 MPa;60°時球體頭部壓力減小到0.51 MPa,相比90°時減小了20.1%;球體30°入水時,球體頭部壓力僅為0.443 MPa,相比90°時減小30.6%。球體未沾濕處的低壓區(qū)出現(xiàn)明顯變化,球體以90°入水時,低壓區(qū)在球體周圍對稱分布,兩側(cè)低壓區(qū)面積相同。球體斜射入水時,低壓區(qū)分布不再對稱。球體入水角度為60°和30°時,球體迎流面低壓區(qū)域的面積比背流面小,而迎流面區(qū)域的壓力高于背流面壓力。

    圖18 入水沖擊階段試驗圖像和局部壓力分布云圖 Fig.18 Experimental image and local pressure distribution pattern during water entry impact

    球體以不同角度入水時,空泡形成階段的流場圖像與壓力分布云圖如圖19 所示。對比圖19 中各幅圖像可知,球體以不同角度入水時,在空泡形成階段,球體頭部沾濕區(qū)域都承受極高壓力,隨著入水角度的減小,球體頭部沾濕處壓力逐漸減小。球體以90°入水時,頭部壓力為0.357 MPa,相比入水沖擊時減小44.0%;入水角度為 60°時,球體頭部壓力減小到0.349 MPa,相比入水沖擊時減小31.6%;入水角度為30°的球體,頭部壓力僅為0.323 MPa,相比入水沖擊時減小27.1%??张輧?nèi)部的壓力低于水面上方的空氣壓力,氣流高速流入空泡內(nèi)部。球體以90°入水 時,空泡內(nèi)的高速氣流對稱分布,隨著入水角度的變化,高速氣流的分布不再均勻。30°入水時,空泡內(nèi)部高速氣流僅分布在球體的左下側(cè)。同時,隨著入水角度的減小,流入空泡內(nèi)氣體的流速反而增加。球體以90°入水時,流入空泡內(nèi)氣體的最高流速為51 m/s;60°時氣體的最高流速為73 m/s,相比90°入水增加43.1%;30°時氣體的最高流速為132 m/s,相比90°入水增加158.8%。

    2.2.2 入水角度對空泡閉合和潰滅的影響

    球體以不同角度入水時,空泡閉合階段的流場圖像與壓力分布云圖如圖20 所示。對比圖20 中各圖可 知,入水角度為90°和60°時,當(dāng)空泡處于閉合過程中,液膜向中間收縮,流入空泡內(nèi)的空氣不斷減少,在閉合處形成高壓區(qū)域;球體以30°入水形成空泡后,在空泡閉合時,水體也向中間收縮,但由于流入空泡內(nèi)氣體的流速較低,未在閉合處形成高壓區(qū)。

    球體以不同角度入水時,潰滅階段過程中的流場圖像與壓力分布云圖如圖21 所示。對比圖21 中各圖可知,球體以90°和60°入水時,入水空泡在尾部以小尺度氣泡團的形式潰滅,同時在潰滅夾斷處產(chǎn)生了回射流。當(dāng)球體以30°入水時,空泡潰滅方式截然不同,沒有出現(xiàn)小尺度的氣泡團,而是斷裂成較大的氣腔,夾斷處也沒有回射流。由壓力分布云圖可知,球體入水角度為90°,當(dāng)空泡潰滅時,收縮處的壓力高于空泡內(nèi)部其他區(qū)域壓力,潰滅產(chǎn)生的回射流沖擊到球體表面,氣腔內(nèi)部沖擊區(qū)域的壓力升高;球體入水角度為60°,入水空泡在潰滅時產(chǎn)生的回射流沖向空泡壁,導(dǎo)致氣腔壁處壓力升高;球體入水角度為30°,空泡潰滅時沒有產(chǎn)生回射流,在空泡斷裂處的壓力與周圍水體的壓力基本一致。

    圖21 空泡潰滅階段試驗圖像和局部壓力分布云圖 Fig.21 Experimental image and local pressure distribution pattern during cavity collapse

    進一步對球體入水后的空泡進行分析,分別繪制了球體直徑和最大截面面積隨時間的變化曲線,如圖22 和圖23 所示。對比圖22 中的各條曲線可知,球體以不同角度入水過程中,在空泡閉合時,泡徑達到最大值,且空泡直徑隨著入水角度的增大反而減小。球體90°入水后,空泡閉合時的直徑為2.2D;60°空泡閉合時的直徑為3.03D,是90°空泡直徑的1.38 倍;30°空泡閉合直徑為5.32D,是90°空泡的2.42 倍。對比圖23 中的各條曲線可知,球體入水角度為30°時,空泡分別在9、22、35 ms 出現(xiàn)了閉合、一次夾斷潰滅和二次潰滅。在0~9 ms、9~22 ms、22~3 ms 時,空泡面積不斷增加,但增加幅度各不相同。在9、22、35 ms 時,區(qū)域空泡面積達到極大值,35 ms 之后,殘余空泡不斷收縮潰滅。當(dāng)球體入水角度為60°和90°時,空泡面積不斷波動,與入水角度30°時空泡面積的變化趨勢相似。在空泡閉合后,入水角度為60°和90°時,空泡面積遠小于30°時的空泡面積,且隨著入水角度的增大反而減小。球體在水下運動的受力分析如圖24 所示??张蓍]合后,球體在水下運動時,受到的力主要有自身的重力mg、附加質(zhì)量力Fa、浮力Fb、與運動速度平方成正比且與速度方向相反的流動阻力Fd以及空泡壁面的表面張力Fσ。受力平衡方程為:

    圖22 空泡直徑隨時間變化曲線 Fig.22 Cavity diameter curve over time

    圖23 空泡面積隨時間變化曲線 Fig.23 Curve of cavity area over time

    圖24 球體水下運動受力分析 Fig.24 Force analysis of sphere in underwater motion

    式中:m為球體質(zhì)量;a為球體運動加速度。Fd、Fa、Fb、Fσ合力方向與運動方向相反,被稱為總流體阻力F(t)。由于球體在水下運動時,表面張力Fσ對球體運動的影響非常微弱,可以忽略不計。根據(jù)以往學(xué)者的研究[30-31],附加質(zhì)量力Fa可以定義為:

    式中:Cm為附加質(zhì)量系數(shù);ρl為液體密度;V為球體體積;Vc為空泡體積;a為球體加速度。式(15)可以表示為:

    因此,總流動阻力系數(shù)為:

    由以上分析可知,球體在水下攜帶空泡運動的總流體阻力系數(shù)與空泡的體積有關(guān),空泡體積越大,阻力系數(shù)越大。又由于空泡的體積與空泡的最大截面積呈正相關(guān),因此空泡最大截面積越大,球體水下運動時的總流動阻力系數(shù)越大。由圖23 可知,在空泡閉合后,入水角度越小,空泡最大截面積越大,總流動阻力系數(shù)Cd越大。

    3 結(jié)論

    本文結(jié)合高速觀測技術(shù)與數(shù)值模擬開展了入水角度對球體入水空泡演化過程影響研究,探討了球體入水空泡的演化特性,分析了入水角度對入水空泡的影響。主要結(jié)論如下:

    1)依據(jù)球體入水過程中空泡的形態(tài)特征,將整個過程分為入水沖擊、空泡形成、空泡閉合和空泡潰滅4 個階段。在入水沖擊階段,球體頭部沖擊水面,生成對稱的水冠;在空泡形成階段,球體運動擠壓水體,在水中形成空泡,空泡在球體的運動下不斷被拉長;在空泡閉合階段,水冠向中間收縮,直至發(fā)生面閉合,斷開空泡與外界空氣的連接,并且空泡在球體拉扯下脫離水面;空泡潰滅階段,入水空泡尾部以小尺度氣泡團的形式泄氣潰滅,同時在潰滅夾斷處產(chǎn)生回射流,回射流沖擊空泡壁面,造成空泡進一步潰滅。

    2)球體入水過程中,空泡不斷發(fā)展,同時整個流場壓力也在不斷發(fā)生明顯變化。在入水沖擊階段,球體沖擊水面,水體與球體接觸區(qū)域壓力急劇升高,球體受到極大的沖擊力導(dǎo)致其速度迅速減小。在空泡形成階段,空泡內(nèi)部壓力低于水面氣壓,氣體連續(xù)流入空泡內(nèi)部,同時球體頭部區(qū)域水體的壓力比入水沖擊時更低。在空泡閉合階段,噴濺內(nèi)、外側(cè)壓力存在壓差,內(nèi)部較低,外側(cè)較高,在內(nèi)外壓差作用下,噴濺向中間收縮,并實現(xiàn)閉合,且進入空泡的空氣流速越高,壓力差越大,空泡閉合也越快。在空泡潰滅時,夾斷潰滅產(chǎn)生的回射流導(dǎo)致球體速度出現(xiàn)波動。

    3)球體入水過程中,入水角度對空泡形態(tài)以及球體阻力系數(shù)的影響較大,且球體的運動速度和所受沖擊載荷因入水角度不同存在較大差異。當(dāng)球體斜射入水時,產(chǎn)生的水冠不對稱。在空泡閉合后,空泡面積隨入水角度的增大而減小,拖帶著空泡的球體在運動時,總流體阻力系數(shù)隨入水角度的增大而不斷減小。在夾斷潰滅階段,空泡夾斷產(chǎn)生的回射流強度隨著入水角度的增大而增大。在入水沖擊時,隨入水角度增大,球體頭部壓力和其加速度最大絕對值不斷增大。

    猜你喜歡
    空泡球體射流
    深海逃逸艙射流注水均壓過程仿真分析
    低壓天然氣泄漏射流擴散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    水下航行體雙空泡相互作用數(shù)值模擬研究
    計算機生成均值隨機點推理三、四維球體公式和表面積公式
    消費電子(2020年5期)2020-12-28 06:58:27
    廣告創(chuàng)意新方法——球體思維兩極法
    Optimization of rice wine fermentation process based on the simultaneous saccharification and fermentation kinetic model☆
    基于LPV的超空泡航行體H∞抗飽和控制
    基于CFD的對轉(zhuǎn)槳無空泡噪聲的仿真預(yù)報
    船海工程(2015年4期)2016-01-05 15:53:28
    射流齒形噴嘴射流流場與氣動聲學(xué)分析
    SPH在水下高速物體空泡發(fā)展模擬中的應(yīng)用
    計算物理(2014年1期)2014-03-11 17:00:22
    国产熟女欧美一区二区| 国产av在哪里看| 日韩中字成人| 日韩成人av中文字幕在线观看| 国产精品熟女久久久久浪| 狠狠精品人妻久久久久久综合| 欧美日本视频| 国产黄频视频在线观看| 最近视频中文字幕2019在线8| 精品久久久久久久末码| 天堂俺去俺来也www色官网 | 免费不卡的大黄色大毛片视频在线观看 | 国语对白做爰xxxⅹ性视频网站| 亚洲av不卡在线观看| 国产伦在线观看视频一区| 亚洲精品亚洲一区二区| 国产成人福利小说| 五月天丁香电影| av在线老鸭窝| 你懂的网址亚洲精品在线观看| 草草在线视频免费看| 69人妻影院| 亚洲一区高清亚洲精品| 亚洲最大成人中文| 美女大奶头视频| 超碰av人人做人人爽久久| 非洲黑人性xxxx精品又粗又长| 成年女人在线观看亚洲视频 | 午夜激情久久久久久久| 男人爽女人下面视频在线观看| av在线老鸭窝| 成人毛片60女人毛片免费| 九九久久精品国产亚洲av麻豆| 日本av手机在线免费观看| 国产亚洲一区二区精品| 国产精品99久久久久久久久| 精华霜和精华液先用哪个| 午夜亚洲福利在线播放| 亚洲av不卡在线观看| 久久久久国产网址| 亚洲成人中文字幕在线播放| 超碰av人人做人人爽久久| 青春草国产在线视频| 国产成人午夜福利电影在线观看| 婷婷色麻豆天堂久久| 97在线视频观看| 国产精品嫩草影院av在线观看| 在线观看人妻少妇| 国产午夜福利久久久久久| 美女黄网站色视频| 禁无遮挡网站| 精品一区二区三卡| 成人毛片60女人毛片免费| 久久人人爽人人爽人人片va| 亚洲av免费高清在线观看| av在线老鸭窝| 亚洲精品日韩在线中文字幕| 在线观看美女被高潮喷水网站| 日本黄色片子视频| 欧美激情久久久久久爽电影| 少妇被粗大猛烈的视频| 18禁裸乳无遮挡免费网站照片| 亚洲天堂国产精品一区在线| 午夜福利网站1000一区二区三区| or卡值多少钱| 欧美日韩亚洲高清精品| 黄色日韩在线| 美女国产视频在线观看| 简卡轻食公司| 亚洲无线观看免费| 又大又黄又爽视频免费| 日日摸夜夜添夜夜爱| 久久久久久久国产电影| 黄色一级大片看看| 亚洲一级一片aⅴ在线观看| 最近手机中文字幕大全| 国产成人精品一,二区| 大话2 男鬼变身卡| 我要看日韩黄色一级片| 久久久久久久国产电影| 亚洲av日韩在线播放| 国产男人的电影天堂91| 久久精品国产亚洲网站| 人妻制服诱惑在线中文字幕| 日韩大片免费观看网站| 能在线免费观看的黄片| 亚洲,欧美,日韩| 亚洲乱码一区二区免费版| 少妇裸体淫交视频免费看高清| 亚洲精品456在线播放app| 天天一区二区日本电影三级| 精品国产三级普通话版| 日韩亚洲欧美综合| 日韩在线高清观看一区二区三区| 久久精品久久久久久噜噜老黄| 六月丁香七月| 一级毛片aaaaaa免费看小| 久久久亚洲精品成人影院| 搡老乐熟女国产| 超碰97精品在线观看| 人妻制服诱惑在线中文字幕| 精品一区二区免费观看| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 午夜福利在线在线| 女人十人毛片免费观看3o分钟| 内地一区二区视频在线| 永久免费av网站大全| 国产精品日韩av在线免费观看| 久久99热6这里只有精品| 不卡视频在线观看欧美| 99久久中文字幕三级久久日本| 插阴视频在线观看视频| 波多野结衣巨乳人妻| 久久99精品国语久久久| 美女内射精品一级片tv| a级一级毛片免费在线观看| 免费看a级黄色片| 中国国产av一级| 如何舔出高潮| 欧美xxⅹ黑人| 亚洲av一区综合| 熟妇人妻不卡中文字幕| 边亲边吃奶的免费视频| 日韩成人av中文字幕在线观看| 免费黄网站久久成人精品| 欧美三级亚洲精品| 青春草亚洲视频在线观看| 国产精品久久久久久久久免| 国产老妇女一区| 国产精品一区二区三区四区久久| 国产中年淑女户外野战色| 麻豆成人av视频| 久久精品夜夜夜夜夜久久蜜豆| 国模一区二区三区四区视频| 黄色一级大片看看| 午夜激情久久久久久久| 国产精品熟女久久久久浪| 成人亚洲精品av一区二区| 色综合色国产| 2018国产大陆天天弄谢| 精品久久久久久成人av| 91狼人影院| 亚洲精品日韩在线中文字幕| 久久97久久精品| 久99久视频精品免费| 国产淫语在线视频| 久久精品国产自在天天线| 一级毛片我不卡| 联通29元200g的流量卡| 男人舔女人下体高潮全视频| 亚洲四区av| 成人午夜精彩视频在线观看| 亚洲av成人精品一二三区| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 日本午夜av视频| av国产久精品久网站免费入址| 国产色爽女视频免费观看| 男的添女的下面高潮视频| 中国国产av一级| 亚洲精品乱码久久久v下载方式| 国产有黄有色有爽视频| 亚洲av电影在线观看一区二区三区 | 国产精品伦人一区二区| 日本黄大片高清| 日本猛色少妇xxxxx猛交久久| 亚洲一区高清亚洲精品| 我的老师免费观看完整版| 欧美性感艳星| 日韩欧美国产在线观看| 国产精品无大码| 美女内射精品一级片tv| 国产伦理片在线播放av一区| 亚洲av.av天堂| 国产精品美女特级片免费视频播放器| 波野结衣二区三区在线| 2018国产大陆天天弄谢| 亚洲图色成人| 街头女战士在线观看网站| 大片免费播放器 马上看| 青春草视频在线免费观看| 街头女战士在线观看网站| 国模一区二区三区四区视频| 久久久久久久久久成人| 久久久国产一区二区| 一级av片app| 国产有黄有色有爽视频| av在线播放精品| 国产精品麻豆人妻色哟哟久久 | 亚洲美女搞黄在线观看| 久久久午夜欧美精品| 国产成人a∨麻豆精品| 国产日韩欧美在线精品| 男女下面进入的视频免费午夜| 久久久色成人| 麻豆成人av视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 青春草国产在线视频| 在线a可以看的网站| 久久久精品94久久精品| 丝瓜视频免费看黄片| 成人毛片60女人毛片免费| 天美传媒精品一区二区| 国产成人精品一,二区| 亚洲自拍偷在线| 亚洲综合精品二区| 免费大片黄手机在线观看| 天堂√8在线中文| 又爽又黄无遮挡网站| 你懂的网址亚洲精品在线观看| 中国国产av一级| 天堂影院成人在线观看| 亚洲精品成人久久久久久| 禁无遮挡网站| 人人妻人人澡人人爽人人夜夜 | 十八禁国产超污无遮挡网站| 免费看a级黄色片| 午夜免费激情av| 中国美白少妇内射xxxbb| 99热网站在线观看| 中文字幕亚洲精品专区| 自拍偷自拍亚洲精品老妇| 一本久久精品| 麻豆成人av视频| 国产欧美另类精品又又久久亚洲欧美| 亚洲美女搞黄在线观看| 国产 一区精品| 中文乱码字字幕精品一区二区三区 | 床上黄色一级片| 亚洲久久久久久中文字幕| av网站免费在线观看视频 | 国产亚洲精品av在线| 久久久久久伊人网av| a级毛色黄片| 亚洲电影在线观看av| 一级毛片久久久久久久久女| 国产激情偷乱视频一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av黄色大香蕉| 国产一区二区在线观看日韩| 在线免费十八禁| 国产成人精品福利久久| 青青草视频在线视频观看| 国产成人a区在线观看| 国产一区二区三区综合在线观看 | 亚洲精品自拍成人| 大香蕉97超碰在线| 久久鲁丝午夜福利片| 2018国产大陆天天弄谢| 蜜桃久久精品国产亚洲av| 男女啪啪激烈高潮av片| 精品一区二区三区视频在线| 国产精品久久久久久久电影| 男女视频在线观看网站免费| 人妻制服诱惑在线中文字幕| av一本久久久久| 久久久a久久爽久久v久久| 91aial.com中文字幕在线观看| 2021少妇久久久久久久久久久| 搞女人的毛片| 国国产精品蜜臀av免费| 观看免费一级毛片| 精品国产三级普通话版| 亚洲一区高清亚洲精品| 汤姆久久久久久久影院中文字幕 | 国产成人午夜福利电影在线观看| 欧美日韩精品成人综合77777| 国产淫片久久久久久久久| 人人妻人人澡人人爽人人夜夜 | 久久这里有精品视频免费| 亚洲综合色惰| 免费黄网站久久成人精品| 免费看a级黄色片| 成年人午夜在线观看视频 | av在线天堂中文字幕| 日本黄大片高清| 一区二区三区四区激情视频| 亚洲性久久影院| 麻豆国产97在线/欧美| 亚洲天堂国产精品一区在线| 国产伦精品一区二区三区视频9| 成人鲁丝片一二三区免费| 国产永久视频网站| 久久久色成人| 国产成人精品久久久久久| 少妇猛男粗大的猛烈进出视频 | 欧美精品一区二区大全| 毛片女人毛片| 国产精品久久久久久久电影| 日韩欧美精品免费久久| 三级男女做爰猛烈吃奶摸视频| 一个人免费在线观看电影| 久久久久久久久中文| 亚洲av免费在线观看| 如何舔出高潮| 天天躁夜夜躁狠狠久久av| 18禁裸乳无遮挡免费网站照片| 九草在线视频观看| av在线播放精品| 中国美白少妇内射xxxbb| 床上黄色一级片| 国产综合懂色| 我的女老师完整版在线观看| 日日啪夜夜爽| 亚洲国产精品专区欧美| 国产久久久一区二区三区| 欧美日韩在线观看h| 熟女人妻精品中文字幕| 寂寞人妻少妇视频99o| 久久午夜福利片| 黄色配什么色好看| 69人妻影院| 亚洲无线观看免费| 国产精品女同一区二区软件| 亚洲欧美成人精品一区二区| 亚洲精品国产成人久久av| 日韩电影二区| 日韩av免费高清视频| 国产大屁股一区二区在线视频| 亚洲av成人精品一二三区| 国产v大片淫在线免费观看| 国产极品天堂在线| 秋霞在线观看毛片| av在线观看视频网站免费| 熟女人妻精品中文字幕| 男女边摸边吃奶| 国产片特级美女逼逼视频| 全区人妻精品视频| 久久国内精品自在自线图片| 精品熟女少妇av免费看| 午夜福利视频1000在线观看| 精品熟女少妇av免费看| 亚洲va在线va天堂va国产| 免费av不卡在线播放| 国产免费视频播放在线视频 | 大片免费播放器 马上看| 日本熟妇午夜| 一二三四中文在线观看免费高清| 伦理电影大哥的女人| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 日韩av在线大香蕉| 黄色日韩在线| 一区二区三区乱码不卡18| 国产黄片美女视频| 精品久久久久久久末码| 精品国产一区二区三区久久久樱花 | 身体一侧抽搐| 午夜爱爱视频在线播放| 在线观看av片永久免费下载| 少妇熟女欧美另类| 亚洲精品久久久久久婷婷小说| 日韩 亚洲 欧美在线| 91久久精品电影网| 日韩 亚洲 欧美在线| 国产永久视频网站| 99热全是精品| 国产黄频视频在线观看| 久久99热6这里只有精品| 亚洲精品久久久久久婷婷小说| 在线免费观看不下载黄p国产| 久久久亚洲精品成人影院| 乱系列少妇在线播放| 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 日本-黄色视频高清免费观看| 国语对白做爰xxxⅹ性视频网站| 在线a可以看的网站| av国产久精品久网站免费入址| 欧美最新免费一区二区三区| 一本久久精品| 成人毛片60女人毛片免费| 国产老妇伦熟女老妇高清| 国产精品一区二区性色av| 亚洲精品日本国产第一区| 1000部很黄的大片| 免费黄色在线免费观看| 国产午夜精品久久久久久一区二区三区| 三级国产精品片| 只有这里有精品99| 欧美成人精品欧美一级黄| 精华霜和精华液先用哪个| 久久久久久久大尺度免费视频| 亚洲成人中文字幕在线播放| 高清毛片免费看| 菩萨蛮人人尽说江南好唐韦庄| 别揉我奶头 嗯啊视频| 色综合亚洲欧美另类图片| 国产在线一区二区三区精| 啦啦啦韩国在线观看视频| 精品久久久久久久久亚洲| 91精品国产九色| 婷婷色综合大香蕉| 欧美高清性xxxxhd video| 日韩制服骚丝袜av| 老师上课跳d突然被开到最大视频| 在线观看一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 丰满人妻一区二区三区视频av| 国产 一区 欧美 日韩| 亚洲精品乱久久久久久| 联通29元200g的流量卡| 搡老乐熟女国产| 麻豆成人av视频| av在线播放精品| 久久久久久久久久久免费av| 男插女下体视频免费在线播放| 日本-黄色视频高清免费观看| 久久久久网色| 蜜臀久久99精品久久宅男| 成人鲁丝片一二三区免费| 婷婷色综合大香蕉| 天堂影院成人在线观看| 国产av不卡久久| 免费人成在线观看视频色| 久久久精品94久久精品| 搞女人的毛片| 一区二区三区免费毛片| 国产免费一级a男人的天堂| 中文乱码字字幕精品一区二区三区 | 毛片一级片免费看久久久久| 精品国产一区二区三区久久久樱花 | 日韩欧美 国产精品| 午夜爱爱视频在线播放| 国产午夜精品久久久久久一区二区三区| 国产v大片淫在线免费观看| 成人国产麻豆网| 老师上课跳d突然被开到最大视频| 亚洲18禁久久av| 成人午夜精彩视频在线观看| 校园人妻丝袜中文字幕| 色网站视频免费| 一本一本综合久久| 男女啪啪激烈高潮av片| 久热久热在线精品观看| 色综合亚洲欧美另类图片| 亚洲人成网站高清观看| 97在线视频观看| 麻豆成人av视频| 色播亚洲综合网| 夜夜看夜夜爽夜夜摸| 国产人妻一区二区三区在| 午夜福利在线在线| 亚洲人成网站在线播| 国产精品一及| 国产在线男女| av在线播放精品| 欧美日韩综合久久久久久| 五月天丁香电影| 久久久午夜欧美精品| 晚上一个人看的免费电影| 麻豆国产97在线/欧美| 亚洲天堂国产精品一区在线| 国产一区二区亚洲精品在线观看| 亚洲在久久综合| 日日干狠狠操夜夜爽| 欧美日韩综合久久久久久| 五月天丁香电影| 2022亚洲国产成人精品| 狠狠精品人妻久久久久久综合| 免费电影在线观看免费观看| 美女被艹到高潮喷水动态| 乱人视频在线观看| h日本视频在线播放| 男人和女人高潮做爰伦理| 人体艺术视频欧美日本| 欧美xxxx性猛交bbbb| 美女大奶头视频| 精品亚洲乱码少妇综合久久| 亚洲欧洲日产国产| 日韩国内少妇激情av| 能在线免费看毛片的网站| 国产高潮美女av| 成年免费大片在线观看| 99热网站在线观看| 波多野结衣巨乳人妻| av在线播放精品| 日本免费在线观看一区| 日韩一区二区三区影片| 九九久久精品国产亚洲av麻豆| 大又大粗又爽又黄少妇毛片口| 欧美性猛交╳xxx乱大交人| 看十八女毛片水多多多| 男人爽女人下面视频在线观看| 国产精品福利在线免费观看| 直男gayav资源| 国产国拍精品亚洲av在线观看| 久久亚洲国产成人精品v| 午夜福利视频精品| 中文天堂在线官网| 久久久久久久大尺度免费视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产精品专区欧美| 午夜视频国产福利| 在线观看免费高清a一片| 欧美三级亚洲精品| 日本黄色片子视频| 日日撸夜夜添| 亚洲精品一二三| 国产免费视频播放在线视频 | 国产单亲对白刺激| 爱豆传媒免费全集在线观看| 中国美白少妇内射xxxbb| 国产精品久久视频播放| 亚洲激情五月婷婷啪啪| 蜜桃久久精品国产亚洲av| 久久人人爽人人爽人人片va| 伊人久久精品亚洲午夜| 高清欧美精品videossex| 久久99精品国语久久久| 日韩欧美国产在线观看| 最近视频中文字幕2019在线8| 成人亚洲精品一区在线观看 | av一本久久久久| 女人久久www免费人成看片| 成年人午夜在线观看视频 | 非洲黑人性xxxx精品又粗又长| 最新中文字幕久久久久| 日本一本二区三区精品| 一级二级三级毛片免费看| 精品国产一区二区三区久久久樱花 | 亚洲熟女精品中文字幕| 国产综合懂色| 免费看美女性在线毛片视频| 国产精品久久久久久久电影| 亚洲一区高清亚洲精品| 国产欧美另类精品又又久久亚洲欧美| 麻豆成人av视频| 日本一二三区视频观看| 欧美日韩在线观看h| 又黄又爽又刺激的免费视频.| 嫩草影院精品99| 免费av不卡在线播放| 97精品久久久久久久久久精品| 久久精品夜夜夜夜夜久久蜜豆| 可以在线观看毛片的网站| 一二三四中文在线观看免费高清| 91久久精品电影网| 99热这里只有精品一区| 亚洲国产欧美人成| 人人妻人人澡人人爽人人夜夜 | 国产免费福利视频在线观看| 蜜桃亚洲精品一区二区三区| 国产精品一区二区三区四区久久| 久久鲁丝午夜福利片| 成年人午夜在线观看视频 | 青春草国产在线视频| 日日撸夜夜添| 国产伦在线观看视频一区| 日本一二三区视频观看| 亚洲国产精品sss在线观看| av在线老鸭窝| 免费观看性生交大片5| 麻豆久久精品国产亚洲av| 亚洲一区高清亚洲精品| 亚洲av福利一区| 黄色配什么色好看| 乱系列少妇在线播放| 亚洲美女视频黄频| 看黄色毛片网站| 好男人视频免费观看在线| 男人舔奶头视频| 中文在线观看免费www的网站| av在线天堂中文字幕| 欧美三级亚洲精品| 极品教师在线视频| 床上黄色一级片| h日本视频在线播放| 一级黄片播放器| 国产 亚洲一区二区三区 | 国产精品人妻久久久久久| 亚洲va在线va天堂va国产| 欧美日韩在线观看h| 精品一区二区三卡| 男女视频在线观看网站免费| 汤姆久久久久久久影院中文字幕 | 亚洲av成人精品一二三区| 午夜免费男女啪啪视频观看| 午夜老司机福利剧场| av播播在线观看一区| 成人综合一区亚洲| 久久韩国三级中文字幕| 2022亚洲国产成人精品| 国产精品麻豆人妻色哟哟久久 | 亚洲av日韩在线播放| 狠狠精品人妻久久久久久综合| av一本久久久久| 欧美日韩一区二区视频在线观看视频在线 | 国内精品宾馆在线| 日韩伦理黄色片| 国产又色又爽无遮挡免| 高清欧美精品videossex| 3wmmmm亚洲av在线观看| 久久久久久久久中文| 午夜福利视频1000在线观看| eeuss影院久久| 精品久久久久久久人妻蜜臀av| 一级毛片我不卡| 97人妻精品一区二区三区麻豆| a级一级毛片免费在线观看| 一区二区三区乱码不卡18| 91av网一区二区| 国产伦精品一区二区三区四那| 午夜老司机福利剧场| 免费观看a级毛片全部| 久久久精品免费免费高清| 一本久久精品| 精品久久久久久久久亚洲| 国产午夜精品一二区理论片| 一级毛片我不卡| 国产一区二区在线观看日韩| kizo精华| 久久久久网色| 亚洲综合色惰|