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

    上浮式水面發(fā)射筒彈射內(nèi)彈道特性

    2022-07-05 06:56:40曾培高姜毅楊哩娜
    兵工學(xué)報 2022年6期
    關(guān)鍵詞:發(fā)射筒彈體水面

    曾培高,姜毅,楊哩娜

    (北京理工大學(xué) 宇航學(xué)院,北京 100081)

    0 引言

    上浮式水面發(fā)射筒是一種新型導(dǎo)彈發(fā)射系統(tǒng),其工作原理是,發(fā)射筒在深水區(qū)域釋放,受浮力作用向水面運動,當(dāng)發(fā)射筒的筒蓋到達(dá)水面時,整個發(fā)射系統(tǒng)獲得一定的出水速度,然后打開筒蓋,采用燃?xì)鈴椛浞绞綄?dǎo)彈發(fā)射出筒,使導(dǎo)彈獲得一定的出筒速度,當(dāng)導(dǎo)彈運動到水面一定高度后,導(dǎo)彈上的發(fā)動機(jī)開始啟動工作。這種發(fā)射系統(tǒng)能夠充分利用發(fā)射筒在水下的浮力勢能,大大節(jié)省導(dǎo)彈的出水助推能量,因此這種發(fā)射系統(tǒng)具有結(jié)構(gòu)簡單、隱蔽性強(qiáng)等優(yōu)點,而且可以在水下進(jìn)行區(qū)域陣地布置,具有很強(qiáng)的戰(zhàn)略意義。水面發(fā)射過程中,發(fā)射筒處于浮動狀態(tài),受到水的復(fù)雜作用力,不同于地面發(fā)射時發(fā)射筒處于固定狀態(tài),水面發(fā)射的彈筒相對運動過程不同于地面發(fā)射的彈筒相對運動過程,這將直接影響發(fā)射筒彈射內(nèi)彈道性能,進(jìn)而影響導(dǎo)彈的打擊精度,甚至當(dāng)導(dǎo)彈的出筒姿態(tài)與出筒速度達(dá)不到要求時,會引起導(dǎo)彈點火失敗、姿態(tài)失穩(wěn)進(jìn)而導(dǎo)致自毀。

    隨著燃?xì)鈴椛浼夹g(shù)的廣泛應(yīng)用,不同應(yīng)用下的燃?xì)鈴椛鋬?nèi)彈道問題得到了國內(nèi)外研究人員的關(guān)注和研究,其中研究方法普遍采用理論與實驗相結(jié)合的方法。燃?xì)鈴椛涞母拍钭钤缬蒑cKinnis等提出。袁曾鳳采用經(jīng)典內(nèi)彈道理論建立了高、低壓室相似準(zhǔn)則和內(nèi)彈道方程。楊珺凡等圍繞水下彈射過程,利用Realizable-湍流模型的計算流體力學(xué)方法和嵌套動網(wǎng)格技術(shù),對水下彈射過程筒口氣泡發(fā)展進(jìn)行了仿真分析。唐垚等針對多級活塞缸式的新型燃?xì)鈴椛浒l(fā)射方式,采用動態(tài)網(wǎng)格分層技術(shù)模擬導(dǎo)彈的運動,對活塞缸的展開過程、流場變化以及內(nèi)彈道特性進(jìn)行了研究。邱海強(qiáng)等以潛射導(dǎo)彈為研究對象,利用動網(wǎng)格技術(shù)和Mixture模型對導(dǎo)彈彈射出筒后筒口氣泡的發(fā)展情況進(jìn)行了研究,探究了筒口氣泡擴(kuò)張、收縮和拉斷的一系列過程。高賢志等對自力發(fā)射的發(fā)射箱燃?xì)馀艑?dǎo)方案進(jìn)行了分析,設(shè)計了彈箱間隙直接排導(dǎo)方案。傅德彬等采用計算流體動力學(xué)(CFD)方法,分析了內(nèi)外筒間隙大小對同心筒自力發(fā)射排氣的影響,并提出優(yōu)化設(shè)計方案。陳舟等討論了不同水深對水下發(fā)射裝置的最大膛壓以及彈丸出膛速度的影響,提出了不同水深下的內(nèi)彈道計算方法。周鵬等使用Spalart-Allmaras湍流模型,采用動網(wǎng)格技術(shù),對彈丸在高壓氣體作用下的運動規(guī)律進(jìn)行了研究。胡曉磊等、李仁鳳等、胡曉磊等采用Realizable-湍流模型、動態(tài)網(wǎng)格分層技術(shù)和有限速率/渦耗散模型,分析了燃燒產(chǎn)物特性對燃?xì)鈴椛鋬?nèi)彈道的影響。趙謝等為了消除低溫燃?xì)鈴椛涞膲毫﹄p峰問題,對環(huán)形腔結(jié)構(gòu)進(jìn)行了優(yōu)化。孟艷等研究了裝配姿態(tài)偏差等對導(dǎo)彈離臺姿態(tài)的影響。徐勤超等對水下提拉缸發(fā)射輕型魚雷的過程進(jìn)行了模擬。牛鈺森等針對一種單提拉桿活塞氣缸式彈射發(fā)射裝置,研究了導(dǎo)軌排氣孔對發(fā)射箱內(nèi)流場的影響。這些研究為上浮式水面發(fā)射筒彈射內(nèi)彈道性能研究提供了良好的研究基礎(chǔ),然而上浮水面發(fā)射筒彈射過程的核心問題,即燃?xì)鈴膹椡查g隙泄漏問題以及發(fā)射筒在水面浮動問題,其相關(guān)研究較少。本文提出了一種新型的導(dǎo)彈發(fā)射系統(tǒng),即上浮式水面發(fā)射筒,對彈射內(nèi)彈道涉及發(fā)射平臺浮動與燃?xì)庑孤﹩栴}進(jìn)行了深入研究。

    上浮式水面發(fā)射筒的發(fā)射過程中,發(fā)射筒在燃?xì)夂秃K墓餐饔孟逻\動,筒內(nèi)的燃?xì)飧綦x板直接受到燃?xì)獾耐屏ψ饔茫M(jìn)而帶動彈體運動。由于加工裝配要求,發(fā)射筒內(nèi)壁與隔離板有一定的間隙,燃?xì)鈺拈g隙中泄漏,整個發(fā)射過程中間隙內(nèi)的流動會變得復(fù)雜,存在亞音速流和超音速流的轉(zhuǎn)化現(xiàn)象,使得零維內(nèi)彈道理論不能精確地預(yù)測低壓室的壓力。本文針對這類復(fù)雜流動問題,對地面固定發(fā)射筒發(fā)射與水面發(fā)射筒發(fā)射分別進(jìn)行了數(shù)值研究,并對比分析了兩種發(fā)射形式下彈射內(nèi)彈道性能與發(fā)射筒的運動軌跡等。

    本文基于Fluent軟件,采用雷諾平均的瞬態(tài)Realizable-湍流模型,求解高壓高溫氣體在筒內(nèi)壁和隔離板間隙中的流動狀態(tài);采用動態(tài)網(wǎng)格分層技術(shù),模擬發(fā)射筒和彈體的各自運動,并考慮了它們的運動對流動的影響;采用Mixture多相流模型以及非均勻初始條件和邊界條件,模擬發(fā)射筒的浮力變化。計算結(jié)果表明,由于水面發(fā)射中發(fā)射筒具有相對較高的設(shè)定出水速度,水面發(fā)射彈體的出筒速度大于地面發(fā)射的出筒速度。

    1 物理模型和計算方法

    1.1 物理模型

    上浮式水面發(fā)射筒的工作原理是,發(fā)射筒在深水區(qū)域釋放后,受浮力作用向上運動,當(dāng)筒蓋運動到水面時,彈與筒具有共同的出水速度,此時筒蓋打開,高壓室開始工作,燃?xì)忾_始推動隔離板和彈體一起向上運動。本文以高壓室開始工作的時刻作為模擬的初始時刻,此時彈筒位置示意圖如圖1所示。上浮式水面發(fā)射筒由高壓室、低壓室、發(fā)射筒、底座、隔離板和4根導(dǎo)軌等組成。監(jiān)測點為實驗和數(shù)值計算的壓力觀測點,用于監(jiān)測低壓室內(nèi)的壓力變化,監(jiān)測點位置隨彈體一起運動,用于監(jiān)測間隙附近的壓力變化。由于導(dǎo)軌的使用,隔離板與發(fā)射筒內(nèi)壁和導(dǎo)軌存在一定的間隙,其中最小間隙位于圖1中-截面處。如圖2所示,隔離板與導(dǎo)軌的間隙約占總間隙面積的30%,大于隔離板與筒內(nèi)壁的最小間隙。

    圖1 發(fā)射筒發(fā)射結(jié)構(gòu)示意圖Fig.1 Schematic diagram of gas-ejection launcher

    圖2 最小間隙示意圖Fig.2 Schematic diagram of minimum gap

    1.2 計算方法

    1.2.1 控制方程

    由于發(fā)射筒不具有軸對稱性,采用三維多組分Navier-Stokes(N-S)控制方程作為質(zhì)量、動量與能量交換以及組分?jǐn)U散的控制方程。其中:

    質(zhì)量交換方程即連續(xù)方程,形式為

    (1)

    式中:為混合氣體密度;為流動速度矢量。

    動量交換方程為帶有黏性流動的N-S方程,

    (2)

    式中:為混合氣體壓強(qiáng);為重力加速度,由于模型需要計算水對筒的作用力,液相的質(zhì)量力不能忽略,重力加速度需要保留;′為湍流引起的雷諾應(yīng)力張量,由湍流模型確定;為分子黏性引起的應(yīng)力張量,

    (3)

    為分子動力黏性系數(shù);為單位矩陣。

    組分?jǐn)U散方程為空氣- 燃?xì)怆p組分輸運方程,求解組分的質(zhì)量分?jǐn)?shù)

    (4)

    式中:為化學(xué)反應(yīng)中組分的生成率,模型不考慮化學(xué)反應(yīng),為0;為其他相向組分轉(zhuǎn)移的生成率,模型不考慮相變質(zhì)量轉(zhuǎn)移過程,因此該項為0;為組分的擴(kuò)散通量,

    (5)

    ,m為分子黏性引起的組分質(zhì)量擴(kuò)散系數(shù),,t為湍流引起的組分質(zhì)量擴(kuò)散系數(shù),由湍流模型確定,T,為組分的熱致擴(kuò)散系數(shù),為溫度。

    能量交換方程為總內(nèi)能輸運方程,由(6)式求解氣體單位質(zhì)量的總內(nèi)能,

    (6)

    式中:為氣體導(dǎo)熱系數(shù);為湍流引起的導(dǎo)熱系數(shù),由湍流模型確定;為組分的單位質(zhì)量焓;為能量源項,模型中不加入相變傳熱過程,因此該項為0。

    1.2.2 湍流模型

    發(fā)射過程中,高壓室內(nèi)的燃?xì)饨?jīng)過噴管進(jìn)入低壓室時會形成超音速射流,提高低壓室的壓力。隨著低壓室的壓力升高,低壓室內(nèi)的燃?xì)鈺耐矁?nèi)壁和隔離板間隙中流出,發(fā)生二次射流現(xiàn)象。兩次射流過程中均涉及壁面的流動分離過程。Realizable-湍流模型能夠較好地模擬剪切流、分離流以及射流等現(xiàn)象,因此模型中采用Realizable-湍流模型模擬湍流效應(yīng)。

    1.2.3 多相流模型

    在高壓室工作的開始時刻,發(fā)射筒下方完全浸在水中,發(fā)射筒口與液面齊平。由于筒具有向上的運動速度,筒的運動為跨液- 氣兩相運動,筒的外表面既受到氣相作用又受到液相作用,因此計算模型需要考慮多相流模型。筒內(nèi)的燃?xì)? 空氣混合氣體在出筒后可能會與筒口附近的液面發(fā)生作用,使得液相和氣相發(fā)生交互。Mixture模型能夠較好地處理相間交互作用,因此模型采用Mixture模型模擬發(fā)射階段的兩相運動過程。

    1.2.4 導(dǎo)彈與發(fā)射筒的運動規(guī)律

    彈射過程中,導(dǎo)彈在燃?xì)馔屏Φ淖饔孟卵刂l(fā)射筒軸線向上運動。彈體在軸線方向上的受力有燃?xì)馔屏?、重力、彈體表面空氣阻力和摩擦力。其合力的公式為

    =---

    (7)

    而在計算開始時刻,發(fā)射筒具有向上的出水速度,在水面上懸浮運動。因此,筒在軸線方向上的受力有氣體阻力、重力、水對筒的作用力和摩擦力。其合力的公式為

    =--+

    (8)

    非定常計算到第個時間步時,導(dǎo)彈沿軸線方向的速度m,和位移m,、發(fā)射筒沿軸線方向的速度c,和位移c,分別由(9)式求得,其中Δ為計算時間步長。

    (9)

    式中:m,-1為第-1個時間步彈受到的合力。

    仿真計算將(7)式~(9)式通過udf加載到Fluent軟件中,動態(tài)計算彈、筒的受力、加速度與速度,采用動態(tài)分層技術(shù)對每個時刻的網(wǎng)格進(jìn)行更新,模擬彈、筒在流場中的運動。

    1.2.5 初邊值條件與網(wǎng)格模型

    由于模型中發(fā)射筒的運動涉及外流場的流動,計算域不僅包含高壓室、低壓室等內(nèi)流場區(qū)域,而且包含由筒外水域和筒口上方氣域組成的外流場區(qū)域。發(fā)射筒的幾何模型雖然不具有軸對稱性,但呈中心對稱,因此可采用1/4模型對發(fā)射筒彈射過程進(jìn)行模擬。計算域以及邊界條件如圖3所示。計算域的外邊界除對稱面外均為壓力出口,內(nèi)流場中高壓室入口設(shè)為壓力入口,作為燃?xì)馊肟谶吔?。彈、筒運動為沿同一方向的直線運動,適合使用動態(tài)網(wǎng)格分層技術(shù)模擬其在流場中的運動過程。動態(tài)網(wǎng)格分層技術(shù)模擬不同速度的多個運動體時,需要采用interface設(shè)置。因此,計算域中,包絡(luò)彈的運動域與包絡(luò)筒的運動域的共同面拆分為一對interface(interface 1),包絡(luò)筒的運動域與周圍靜止域的共同面拆分為一對interface(interface 2)。計算域中壁面條件均為絕熱壁面條件。

    圖3 計算域及邊界條件示意圖(紅色面為壓力出口,綠色面為interface,黃色面為對稱面,藍(lán)色面為壓力入口,灰色面為壁面)Fig.3 Computational domain and boundary conditions(The red faces are the pressure outlets,the green faces are the interfaces,the yellow faces are the symmetries,the blue face is the pressure inlets,and the gray faces are the walls)

    由實驗獲得的壓力入口的總壓隨時間變化曲線如圖4所示,其中為高壓室的參考壓力,為燃?xì)獍l(fā)生劑燃燒的參考時間。由圖4可見:=0時,高壓室的堵蓋打開,高壓室的壓力約為02;=0033時,高壓室的壓力達(dá)到了第1個峰值,約為03;隨后高壓室壓力降低,在=0154時壓力曲線達(dá)到了第1個波谷,約為02;由于高壓室裝藥為增面燃燒藥,在=069時壓力曲線達(dá)到了第2個峰值,約為098。

    圖4 高壓室總壓力曲線Fig.4 Total pressure of high pressure chamber

    發(fā)射計算域的初始值分布如圖5所示。從圖5中可知,在計算開始時刻,筒口與氣液交界面齊平;高壓室的初始壓力為沖開堵蓋的壓力,氣域的初始壓力值為大氣壓值,水域的初始壓力隨水深線性變化;在初始時刻,彈、筒具有相同的出水速度,方向向上。

    圖5 計算域初始值分布云圖Fig.5 Initial contours of water volume fraction and pressure

    在建立網(wǎng)格前,對筒內(nèi)壁與隔離板間隙的網(wǎng)格進(jìn)行了無關(guān)性驗證。圖6為網(wǎng)格無關(guān)性驗證曲線。為了驗證網(wǎng)格的無關(guān)性,本文對4種不同間隙網(wǎng)格尺度的工況進(jìn)行了數(shù)值模擬,分別為0.5 mm、1 mm、2 mm和3 mm。間隙網(wǎng)格尺度為Δ,彈的出筒速度為Δ。圖6中,Δ為間隙網(wǎng)格尺度臨界值。從圖6中可以看出,在間隙網(wǎng)格尺度為2 mm即取臨界值時,出筒速度的相對誤差梯度出現(xiàn)階躍性降低,其相對誤差在0.8%內(nèi),同時該尺度下筒內(nèi)間隙的平均溫度開始趨于一致。因此,間隙網(wǎng)格尺度設(shè)定為2 mm。為了保證運動過程中interface上配對的網(wǎng)格尺度相近,計算域的網(wǎng)格數(shù)約為590萬。

    圖6 網(wǎng)格無關(guān)性驗證曲線Fig.6 Mesh size independence verification curve

    2 數(shù)值方法驗證

    為驗證本文數(shù)值方法的可靠性,對發(fā)射筒做了地面實驗。監(jiān)測點的計算值和實驗值對比如圖7示,其中為實驗參考壓力。由圖7可見:點處數(shù)值結(jié)果與實驗結(jié)果的壓力變化趨勢基本一致。數(shù)值結(jié)果的第1個壓力峰值為0610,實驗值為0624,計算結(jié)果的第2個壓力峰值為0835,實驗值為0792;計算值與實驗值在壓力峰值上的誤差在55內(nèi);實驗中彈體出筒速度為767 m/s,數(shù)值計算中彈出筒速度為820 m/s,計算值和實驗值在彈體出筒速度上的誤差在7內(nèi)。因此,從計算結(jié)果和實驗結(jié)果的對比可以看出,采用本文所使用物理模型和數(shù)值計算方法能夠有效地模擬發(fā)射筒的彈射過程。

    圖7 P0點計算壓力值與實驗結(jié)果對比Fig.7 Comparison of numerical and experimental pressure results at Point P0

    3 計算結(jié)果分析

    3.1 內(nèi)彈道特性分析

    水面彈射啟動前,彈體在發(fā)射筒內(nèi)受到浮力作用一起向上運動,彈筒具有相同的向上運動速度。筒蓋浮出水面后,筒蓋打開,并且高壓室開始工作,高壓室的壓力快速升高。高壓室內(nèi)的壓力升高至破膜壓力后,燃?xì)鉀_開高壓室的堵蓋進(jìn)入低壓室內(nèi),低壓室的壓力升高。隔離板受到低壓室的壓力作用,產(chǎn)生向上的力,當(dāng)隔離板受到的推力能克服隔離板、彈體以及適配器等的重力、彈體運動的空氣阻力和適配器與筒接觸產(chǎn)生的摩擦力時,彈體與發(fā)射筒開始做相對運動。由于隔離板與內(nèi)筒壁有間隙,燃?xì)鈱⑼ㄟ^間隙從低壓室進(jìn)入大氣環(huán)境中,在間隙中形成二次射流,并在附近形成復(fù)雜的波系結(jié)構(gòu),使低壓室的壓力產(chǎn)生較為復(fù)雜的變化。

    地面發(fā)射與水面發(fā)射監(jiān)測點、監(jiān)測點的壓力隨時間變化曲線,彈體加速度、彈體速度增量Δ(即彈體速度與出水速度的差)隨彈筒相對位移變化曲線以及彈筒相對位移隨時間變化曲線如圖8所示,發(fā)射筒的運動曲線如圖9所示,其中為破膜時刻彈筒的共同運動速度(出水速度),為彈體的參考加速度,為低壓室的參考壓力,為環(huán)境壓力,為彈體出筒行程,=1101,=1171,為筒的位移,為筒的加速度。

    由圖8(a)可知,根據(jù)低壓室的壓力升降,可以將發(fā)射過程分為4個階段,分別為段、段、段以及段,為計算初始時刻,點在=008左右,點在=024左右,點在=069左右,為彈出筒時刻。

    在段,點時刻高壓室的堵蓋沖開,高溫高壓燃?xì)獬淙氲蛪菏?,并壓縮低壓室的剩余空氣,低壓室的壓力快速傳播并升高。由圖8(d)可知,在段,水面發(fā)射與地面發(fā)射中彈、筒沒有明顯的相對位移,二者低壓室的容積幾乎沒有變化,而且二者高壓室采用同一壓力曲線,因此二者低壓室內(nèi)的流動幾乎相同,同一監(jiān)測處的壓力也幾乎一致,如圖8(a)中段所示。對比圖4可知,監(jiān)測點出現(xiàn)第1個壓力波峰的點時刻出現(xiàn)高壓室的第1壓力波峰與高壓室的第1壓力波谷之間,此時表壓值約為052,表明在低壓室的壓力到達(dá)第1個壓力波峰時刻,高壓室向低壓室補充的燃?xì)獬霈F(xiàn)了減少。由圖9(a)可知,在時刻,筒受到向上的正浮力作用,隨后低壓室壓力升高,發(fā)射筒受到低壓室的下壓作用,筒向上做減速運動。

    圖8 水面發(fā)射和地面發(fā)射內(nèi)彈道曲線Fig.8 Interior ballistic curves of launching on ground and water-surface

    圖9 筒運動曲線Fig.9 Motion curves of canister on water-surface

    在段,彈筒開始有明顯的相對運動,彈體相對運動使低壓室氣體出現(xiàn)明顯膨脹。點時刻位于高壓室的第1個壓力波谷與第2個壓力波峰之間,因此在段高壓室的壓力先降低后升高。由于高壓室的壓力先降低,低壓室內(nèi)補充的燃?xì)獬霈F(xiàn)減少,但彈筒的相對位移在加速變大,使得低壓室內(nèi)的氣體膨脹加劇,進(jìn)而低壓室壓力出現(xiàn)降低,彈體的加速度也隨之減小。在這一階段,水面發(fā)射與地面發(fā)射中低壓室壓力具有相似的下降趨勢。由于點靠近高壓室,其燃?xì)鈮嚎s效應(yīng)要強(qiáng)于燃?xì)馀蛎浶?yīng),而點靠近相對運動的隔離板處,其燃?xì)馀蛎浶?yīng)要強(qiáng)于燃?xì)鈮嚎s效應(yīng),使得點的壓力高于點的壓力。由圖8(b)可知,水面發(fā)射彈筒相對加速度大于地面發(fā)射中彈筒相對加速度(即地面彈體自身的加速度),在點時刻,水面發(fā)射中彈的加速度極小值029要小于地面發(fā)射彈的加速度極小值037,約小21,水面發(fā)射的彈筒相對位移0038 8要大于地面發(fā)射的彈筒相對位移0024 5,其低壓室約多膨脹了58,這也表明了水面發(fā)射低壓室的膨脹速度要大于地面發(fā)射的膨脹速度,水面發(fā)射低壓室的壓力降低得更快。因此,在時刻,水面發(fā)射低壓室表壓的極小值025要小于地面發(fā)射低壓室表壓的極小值029,約小14,如圖8(a)中段所示。

    段高壓室壓力一直處于上升趨勢,即高壓室向低壓室補充的燃?xì)庖恢痹谠龆?。在點,低壓室內(nèi)燃?xì)獾难a充已經(jīng)足夠抵消燃?xì)獾呐蛎涀饔门c泄漏作用,低壓室氣體再次受到補充燃?xì)獾膲嚎s作用,使得低壓室壓力再次升高,并在高壓室壓力達(dá)到第2個峰值的附近時刻,地面發(fā)射和水下發(fā)射的低壓室壓力幾乎同時到達(dá)第2個峰值,即圖8(a)中點附近。在點時刻,由于彈筒相對速度較大,隔離板附近的膨脹現(xiàn)象較為明顯,監(jiān)測點、的壓力值出現(xiàn)明顯的差異,二者的壓力差約為其壓力值的6。由圖8(b)可知,在段,地面發(fā)射與水面發(fā)射燃?xì)馀蛎洭F(xiàn)象與段近似,在時刻,地面發(fā)射彈筒相對位移為0315,而水面發(fā)射彈筒相對位移為0412,約多膨脹了30,使得水面發(fā)射低壓室的平均表壓062低于地面發(fā)射低壓室的平均表壓0744,約低16,水面發(fā)射彈體的最大加速度約為110低于地面發(fā)射的最大加速度1374,約低20。這也使得在時刻水面發(fā)射彈體的速度增量比地面發(fā)射的速度增量小25,如圖8(c)所示。

    在段,高壓室壓力開始降低,直至降至環(huán)境壓力。由于彈筒相對位移繼續(xù)增大,低壓室的氣體膨脹作用逐漸占主導(dǎo)地位,監(jiān)測點、的壓力值逐漸貼近,同時隨著彈體完全出筒,水面發(fā)射和地面發(fā)射筒內(nèi)的壓力都降至環(huán)境壓力附近,監(jiān)測點上的壓力也就趨于一致。由圖8(d)可知,水面發(fā)射彈體的出筒時間=1013要快于地面發(fā)射的出筒時間=1171,約快16。由圖8(c)中可知,水面發(fā)射彈體的出筒速度增量小于地面發(fā)射的出筒速度增量,彈體的動量增量約損失15。由于水面發(fā)射筒有初始出水速度,使得水面發(fā)射彈體的出筒速度為,大于地面發(fā)射彈體的出筒速度。由圖9(b)可知,從高壓室破膜后到彈體出筒過程中,雖然筒的運動一直受到燃?xì)獾淖璧K作用,但由于其具有出水速度,其位移一直為正,使得筒口一直處于水面上方。

    3.2 流場變化規(guī)律分析

    由于高壓室的壓力與低壓室的壓力相差較大以及低壓室的壓力和環(huán)境壓力也相差較大,復(fù)雜的波系結(jié)構(gòu)主要存在于高壓室噴管出口附近以及隔離板與筒壁間隙處及其出口附近,波系的變化可能會導(dǎo)致高壓室的出口流量以及低壓室的流出流量發(fā)生變化,進(jìn)而影響內(nèi)彈道的性能。

    高壓室的出口流量隨時間變化曲線如圖10所示,其中為參考流量。由圖10可見:由于水面發(fā)射和地面發(fā)射采用相同的高壓室條件,雖然二者低壓室的壓力存在明顯差異,但在段,高壓室壓力足夠大,噴管喉部馬赫數(shù)一直等于1,因此在段噴管流量一直處于壅塞狀態(tài),二者高壓室的出口流量幾乎一致;在段,高壓室的壓力進(jìn)入衰減階段,雖然高壓室噴管內(nèi)流動存在超音速轉(zhuǎn)向亞音速的轉(zhuǎn)變過程,但高壓室噴管內(nèi)流動轉(zhuǎn)為亞音速時,水面發(fā)射和地面發(fā)射低壓室內(nèi)的壓力幾乎趨于一致,因此整個發(fā)射階段,水面發(fā)射和地面發(fā)射高壓室向低壓室內(nèi)補充的燃?xì)饬渴菐缀跸嗤摹?/p>

    圖10 高壓室出口流量隨時間變化曲線Fig.10 Mass flux in high pressure chamber on ground and water-surface

    由于隔離板和內(nèi)筒壁存在間隙,低壓室內(nèi)的氣體受到高壓室流入的燃?xì)鈮嚎s后便會從間隙中流向外界環(huán)境中,并進(jìn)一步膨脹。低壓室的流出流量隨時間變化曲線如圖11所示。由圖11可見:低壓室的第1個流出流量峰值到達(dá)時刻與高壓室的第1個壓力峰值到達(dá)時相近,由圖8(a)可知,此時水面發(fā)射和地面發(fā)射低壓室的壓力幾乎一致,因此二者流出流量的第1個峰值也幾乎一致,約為0095,約占流入低壓室燃?xì)饬髁康?3;低壓室的第2個流出流量峰值約在高壓室的第2個壓力峰值時達(dá)到,其流量約占流入低壓室燃?xì)饬髁康?。此時水面發(fā)射低壓室的壓力低于地面發(fā)射低壓室的壓力,因此水面發(fā)射低壓室的流出流量小于地面發(fā)射的流出流量。雖然第1個峰值時低壓室的壓力低于第2個峰值時低壓室的壓力,但由于在第1個峰值時低壓室流出的氣體為常溫空氣,而在第2個峰值時低壓室流出的氣體為高溫燃?xì)?,常溫空氣的密度遠(yuǎn)遠(yuǎn)高于高溫燃?xì)獾拿芏?,使得低壓室的?個流出流量峰值大于第2個流出流量峰值。

    圖11 低壓室流出流量隨時間變化曲線Fig.11 Mass fluxes in low pressure chamber on ground and water-surface

    間隙處馬赫數(shù)1等值圖如圖12所示。由圖12可見:在=0017至=098時,最小間隙處流動一直處于超音速流動狀態(tài),即紅色等值面包絡(luò)了最小間隙,因此在這一階段,低壓室從最小間隙流出的流量也一直處于壅塞狀態(tài),壅塞時間約占總發(fā)射時間的95;由于隔離板與導(dǎo)軌的間隙要寬于最小間隙,且占總間隙面積的30左右,在整個發(fā)射過程導(dǎo)軌處間隙的流動在亞音速流動和超音速流動中反復(fù)變化。其具體變化有,在=008時導(dǎo)軌上端間隙剛好被紅色等值面完全封閉,在=017時封閉該間隙的等值面發(fā)生破碎,在=035時該間隙再次被等值面封閉,在=083時封閉該間隙的等值面再次破碎。因此,在=008前、=017至=035以及=083后導(dǎo)軌間隙處的流動均有亞音速流,其時間約占總發(fā)射時間的45,其余時間處于超音速流動狀態(tài)。由于導(dǎo)軌處間隙的流出流量比重較大,而且其處于亞音速流動狀態(tài)占時較長,因此無法直接采用零維內(nèi)彈道理論計算間隙的流出量,同時也無法采用二維等效模型精確計算出間隙的流出量,進(jìn)一步表明采用三維計算的必要性。

    圖12 間隙處馬赫數(shù)1等值面圖Fig.12 Mach=1 iso-surfaces at the gap on water-surface

    水面發(fā)射筒口附近的流動如圖13所示,左側(cè)為速度矢量圖,速度上限為20 m/s,右側(cè)為溫度云圖,溫度上限為370 K。從圖13中可知,在=0325時,燃?xì)忾_始從筒口溢出并直接向上方運動,此時筒口與氣- 液交界面處有一定的距離,氣- 液交界面處僅受常溫低速的空氣作用,該條件不足以引起相變發(fā)生。在低壓室的壓力到達(dá)最大時=069,以及隔離板快出筒時=10,從筒口流出的燃?xì)庖恢北3窒蛏线\動,雖然燃?xì)獬鐾部诤蟀l(fā)生了膨脹,但筒口遠(yuǎn)離氣- 液交界面,燃?xì)獠粫苯幼饔糜谝后w表面,氣- 液交界面處依然僅受常溫低速的空氣作用。因此,在水面發(fā)射過程中,由于在燃?xì)饬鞒鐾部跁r筒口已遠(yuǎn)離氣液交界面,燃?xì)獠粫苯优c水面發(fā)生作用,僅有低速空氣會與水面產(chǎn)生微弱的作用,這一過程對發(fā)射的影響可以忽略不計。

    圖13 筒口流動速度矢量圖與溫度云圖Fig.13 Vector diagrams and temperature contours near the canister lip on water-surface

    3.3 出水速度對內(nèi)彈道的影響

    由于發(fā)射筒受到燃?xì)獾臏p速作用,當(dāng)發(fā)射筒的出水速度不足時,發(fā)射筒的筒口將沒入水下,水將傾入發(fā)射筒內(nèi),在高溫燃?xì)庾饔孟拢畬l(fā)生相變,消耗燃?xì)獯罅康哪芰?,大大降低彈的出筒速度,甚至?xí)绊憦椀陌l(fā)射安全性,這種情況將視為發(fā)射失敗。由31節(jié)的內(nèi)彈道性能分析可知,水面發(fā)射彈體的出筒速度增量要小于地面發(fā)射的出筒速度增量,但由于水面發(fā)射中發(fā)射筒有出水速度,水面發(fā)射彈體的出筒速度要快于地面發(fā)射的出筒速度。為了使發(fā)射安全,以及彈體的出筒速度不小于地面發(fā)射的出筒速度,水面發(fā)射筒的出水速度需要滿足一定的條件,因此有必要分析出水速度對水面發(fā)射內(nèi)彈道的影響。

    (10)

    式中:為導(dǎo)彈的位移;為筒的位移;為重力和摩擦力提供的加速度。從(10)式可知,彈筒的相對位移與彈筒的出水速度無關(guān),而彈、筒受到的燃?xì)饬H與彈筒相對位移以及其相對速度有關(guān),因此燃?xì)饬? 相對位移變化曲線不會隨出水速度發(fā)生改變。

    當(dāng)考慮水對筒的作用力、彈受到的空氣阻力作用時,二者可以近似表達(dá)為

    (11)

    Δ=((-)+05(-))+
    (-)+(-)

    (12)

    式中:、分別為彈體速度、下筒的加速度。

    從(12)式可知,由于發(fā)射筒出水速度有限,當(dāng)前條件下,其不大于10 m/s,因此速度項產(chǎn)生的水對筒的作用力差值與燃?xì)饬ο啾冗h(yuǎn)遠(yuǎn)較小。發(fā)射筒的加速度主要由燃?xì)饬ω暙I(xiàn),對于不同的出水速度,發(fā)射筒的加速度幾乎相同。

    3種不同出水速度的彈體速度曲線如圖14所示,這3種出水速度分別為08、以及12,彈出筒時刻,發(fā)射筒的筒口均未沒入水下。從圖14中可以看出,不同出水速度的彈體速度曲線幾乎為同一曲線的平移曲線。這與上述理論分析一致,出水速度在一定范圍內(nèi),彈體的出筒速度增量幾乎不隨出水速度改變而改變,彈體的出筒速度為彈體出筒速度增量與出水速度的疊加。由于彈體的出筒速度增量不隨出水速度發(fā)生改變,對于不同的出水速度,彈射用時、彈體受到的最大載荷等幾乎一致,水面發(fā)射與地面發(fā)射相比較分析的主要結(jié)論幾乎不隨出水速度改變而改變。根據(jù)線性疊加原理可知,水面發(fā)射彈體出筒速度不低于地面發(fā)射的出筒速度,出水速度需大于0485,而水不進(jìn)入發(fā)射筒內(nèi),出水速度需大于075。因此,在發(fā)射安全條件下,水面發(fā)射彈體的出筒速度一定大于地面發(fā)射的出筒速度。

    圖14 不同出水速度的彈體速度曲線Fig.14 Velocity curves of projetile at different initial velocities

    4 結(jié)論

    本文提出一種新型發(fā)射方式即水面發(fā)射筒彈射發(fā)射方式,具有耗能少,隱蔽性強(qiáng)等優(yōu)點,在理論上,本文是首次對彈射內(nèi)彈道涉及發(fā)射平臺浮動與燃?xì)庑孤﹩栴}進(jìn)行了深入研究。得出主要結(jié)論如下:

    1)水面發(fā)射和地面發(fā)射具有近似的流動狀態(tài),但水面發(fā)射低壓室的膨脹速度要快于地面發(fā)射的膨脹速度,主要原因是水面發(fā)射時彈筒相對速度更大。水面發(fā)射彈體的出筒速度增量小于地面發(fā)射的出筒速度增量,但由于水面發(fā)射筒有出水速度,水面發(fā)射彈體的出筒速度反而大于地面發(fā)射的出筒速度,而且水面發(fā)射彈體受到的最大載荷僅為地面發(fā)射的80,水面發(fā)射出筒用時僅為地面發(fā)射的86。

    2)水面發(fā)射過程中,在給定出水速度下,燃?xì)饬鞒鐾部跁r筒口已遠(yuǎn)離氣液交界面,燃?xì)獠粫苯优c水面發(fā)生作用,僅有少量的低速空氣會與水面產(chǎn)生微弱的作用,這一過程對發(fā)射的影響可以忽略不計。

    3)在筒口不沒入水時,與筒受到的燃?xì)饬ο啾?,不同出水速度下水對筒的作用力相差較小,筒受到的加速度幾乎相同,使得彈體的出筒速度增量不隨出水速度改變而改變,彈體的出筒速度為彈體出筒速度增量與出水速度的疊加。因此,在安全發(fā)射條件下,水面發(fā)射與地面發(fā)射相比較分析的主要結(jié)論幾乎不隨出水速度改變而改變,水面發(fā)射彈體出筒速度一定大于地面發(fā)射的出筒速度。

    猜你喜歡
    發(fā)射筒彈體水面
    尾錐角對彈體斜侵徹過程中姿態(tài)的影響研究
    某型增強(qiáng)型復(fù)合材料發(fā)射筒強(qiáng)度特性評估
    橢圓截面彈體斜侵徹金屬靶體彈道研究*
    爆炸與沖擊(2022年2期)2022-03-17 07:28:44
    水黽是怎樣浮在水面的
    發(fā)射筒熱力耦合場下力學(xué)性能研究
    地空導(dǎo)彈發(fā)射筒焊接系統(tǒng)的校準(zhǔn)技術(shù)研究
    電子測試(2018年9期)2018-06-26 06:45:46
    降低某產(chǎn)品貯運發(fā)射筒內(nèi)氣體濕度的工藝方法研究
    STOPAQ粘彈體技術(shù)在管道施工中的應(yīng)用
    上海煤氣(2018年6期)2018-03-07 01:03:22
    創(chuàng)造足以亂真的水面反光
    爭奪水面光伏
    能源(2016年3期)2016-12-01 05:11:02
    他把我摸到了高潮在线观看| www.www免费av| 一本久久中文字幕| 国产真实乱freesex| 两人在一起打扑克的视频| 久久中文字幕人妻熟女| 最近最新免费中文字幕在线| 黄色片一级片一级黄色片| 在线免费观看不下载黄p国产 | 波多野结衣高清作品| 90打野战视频偷拍视频| 国产精品美女特级片免费视频播放器 | 又紧又爽又黄一区二区| 国产伦精品一区二区三区视频9 | 一进一出抽搐动态| 国产激情久久老熟女| 亚洲av中文字字幕乱码综合| 国产三级中文精品| 琪琪午夜伦伦电影理论片6080| 国产精品 国内视频| aaaaa片日本免费| 99国产精品99久久久久| 一区二区三区高清视频在线| 成年女人永久免费观看视频| 99精品久久久久人妻精品| 最新美女视频免费是黄的| h日本视频在线播放| 搡老岳熟女国产| 亚洲九九香蕉| 中出人妻视频一区二区| 成人国产综合亚洲| 亚洲激情在线av| 老鸭窝网址在线观看| 亚洲精品456在线播放app | 亚洲国产精品成人综合色| 精品久久蜜臀av无| 亚洲av五月六月丁香网| 丰满人妻一区二区三区视频av | 美女cb高潮喷水在线观看 | 久久久久久国产a免费观看| 午夜福利在线观看吧| 国产激情久久老熟女| 欧美不卡视频在线免费观看| 久久99热这里只有精品18| cao死你这个sao货| 久久午夜综合久久蜜桃| 夜夜看夜夜爽夜夜摸| 国产精品一及| 久久久色成人| 午夜久久久久精精品| 老司机在亚洲福利影院| 亚洲,欧美精品.| 黄片大片在线免费观看| 99精品在免费线老司机午夜| 精品电影一区二区在线| 国产精品爽爽va在线观看网站| 别揉我奶头~嗯~啊~动态视频| 午夜精品一区二区三区免费看| 国产精品99久久99久久久不卡| 国产不卡一卡二| 精品免费久久久久久久清纯| 免费看日本二区| 欧美乱色亚洲激情| 国产不卡一卡二| 欧美成人性av电影在线观看| 国产精品久久电影中文字幕| 国产精品一区二区精品视频观看| 亚洲熟妇中文字幕五十中出| 精品无人区乱码1区二区| 一级毛片高清免费大全| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 成人鲁丝片一二三区免费| 搞女人的毛片| 香蕉av资源在线| 十八禁人妻一区二区| 国产成人啪精品午夜网站| 精品福利观看| 日韩欧美一区二区三区在线观看| 欧美xxxx黑人xx丫x性爽| 欧美国产日韩亚洲一区| 亚洲一区二区三区色噜噜| 亚洲18禁久久av| 亚洲色图 男人天堂 中文字幕| 91老司机精品| 久久久久亚洲av毛片大全| 久久中文字幕一级| 国产精品久久电影中文字幕| 老汉色∧v一级毛片| 国产精品九九99| 美女cb高潮喷水在线观看 | 免费大片18禁| 久久伊人香网站| 三级男女做爰猛烈吃奶摸视频| 97人妻精品一区二区三区麻豆| 亚洲av第一区精品v没综合| 老司机午夜福利在线观看视频| 日本熟妇午夜| 久久精品国产亚洲av香蕉五月| 精品久久久久久久末码| 99热这里只有是精品50| 久久久久久人人人人人| 亚洲av第一区精品v没综合| 99热这里只有精品一区 | 国产伦人伦偷精品视频| 在线播放国产精品三级| 大型黄色视频在线免费观看| 亚洲九九香蕉| 级片在线观看| 伊人久久大香线蕉亚洲五| 最近视频中文字幕2019在线8| 国产单亲对白刺激| 亚洲 欧美 日韩 在线 免费| 淫秽高清视频在线观看| 国产成人啪精品午夜网站| 中文亚洲av片在线观看爽| 美女高潮喷水抽搐中文字幕| 国产日本99.免费观看| 亚洲自偷自拍图片 自拍| 国内揄拍国产精品人妻在线| 精品一区二区三区视频在线 | 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美一区二区综合| 窝窝影院91人妻| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看 | 精品一区二区三区视频在线观看免费| 久久精品综合一区二区三区| 久久久久国产一级毛片高清牌| 手机成人av网站| 香蕉av资源在线| 最近在线观看免费完整版| 真人做人爱边吃奶动态| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久末码| 免费在线观看日本一区| 亚洲自拍偷在线| 91av网站免费观看| 亚洲在线观看片| 神马国产精品三级电影在线观看| 国产成人精品久久二区二区91| 国产精品 欧美亚洲| 午夜两性在线视频| 精品国产乱码久久久久久男人| 亚洲av电影不卡..在线观看| 999精品在线视频| 90打野战视频偷拍视频| 成人三级做爰电影| 久久久久久人人人人人| 制服人妻中文乱码| 12—13女人毛片做爰片一| 天天添夜夜摸| 国产成人啪精品午夜网站| 亚洲中文字幕一区二区三区有码在线看 | 国产精品永久免费网站| 国产高潮美女av| 国产精品女同一区二区软件 | 国产精品,欧美在线| 久久性视频一级片| netflix在线观看网站| av国产免费在线观看| 免费搜索国产男女视频| 国产伦精品一区二区三区四那| 精品久久久久久,| 亚洲国产精品成人综合色| 久久人妻av系列| 久久精品人妻少妇| 欧美zozozo另类| 观看免费一级毛片| www.自偷自拍.com| 国产精品 国内视频| 欧美另类亚洲清纯唯美| 黄色 视频免费看| 国产亚洲欧美在线一区二区| 人妻久久中文字幕网| 久久香蕉精品热| 久久久国产精品麻豆| 老司机福利观看| 精品国产亚洲在线| 国产成人精品久久二区二区91| 夜夜躁狠狠躁天天躁| 午夜免费成人在线视频| 午夜福利视频1000在线观看| 一个人看视频在线观看www免费 | 在线观看免费午夜福利视频| 又紧又爽又黄一区二区| 男女之事视频高清在线观看| 人妻丰满熟妇av一区二区三区| 国产三级在线视频| 后天国语完整版免费观看| 91在线观看av| 日韩欧美三级三区| 欧美另类亚洲清纯唯美| 精品人妻1区二区| 欧美黑人巨大hd| 可以在线观看毛片的网站| 中文字幕人成人乱码亚洲影| 亚洲成人久久爱视频| 欧美精品啪啪一区二区三区| 九九热线精品视视频播放| 成熟少妇高潮喷水视频| 欧美又色又爽又黄视频| 桃红色精品国产亚洲av| 少妇熟女aⅴ在线视频| 亚洲av成人av| 亚洲国产精品久久男人天堂| 色在线成人网| 国产成人aa在线观看| 老司机在亚洲福利影院| 我要搜黄色片| 一级毛片高清免费大全| 五月玫瑰六月丁香| 黄色丝袜av网址大全| 午夜免费成人在线视频| 精品一区二区三区视频在线 | 最新中文字幕久久久久 | 九九在线视频观看精品| 搡老妇女老女人老熟妇| 高潮久久久久久久久久久不卡| 亚洲自拍偷在线| 一本一本综合久久| av福利片在线观看| 观看美女的网站| 淫妇啪啪啪对白视频| 精品久久久久久,| 亚洲欧美日韩高清专用| 白带黄色成豆腐渣| 俄罗斯特黄特色一大片| 黑人欧美特级aaaaaa片| 色噜噜av男人的天堂激情| 成人性生交大片免费视频hd| 白带黄色成豆腐渣| 欧美日韩乱码在线| 日日摸夜夜添夜夜添小说| 国产又色又爽无遮挡免费看| 成年免费大片在线观看| 国产一区二区三区视频了| 亚洲欧美一区二区三区黑人| 久久精品国产综合久久久| 他把我摸到了高潮在线观看| 欧美丝袜亚洲另类 | bbb黄色大片| 网址你懂的国产日韩在线| 亚洲av电影在线进入| 亚洲第一电影网av| 国产一区二区在线观看日韩 | 成人三级做爰电影| 村上凉子中文字幕在线| 久久这里只有精品19| 欧美日韩综合久久久久久 | 国产精品电影一区二区三区| 香蕉av资源在线| 免费看日本二区| 午夜福利高清视频| 欧美成狂野欧美在线观看| 99国产精品一区二区三区| 亚洲精品国产精品久久久不卡| 国产久久久一区二区三区| 97碰自拍视频| 在线观看免费午夜福利视频| 亚洲一区二区三区色噜噜| 亚洲黑人精品在线| 特级一级黄色大片| 亚洲九九香蕉| 黑人操中国人逼视频| 97超视频在线观看视频| www.熟女人妻精品国产| 久久精品夜夜夜夜夜久久蜜豆| 国产麻豆成人av免费视频| 中文亚洲av片在线观看爽| 最近最新中文字幕大全免费视频| 无限看片的www在线观看| 久久久水蜜桃国产精品网| 亚洲国产中文字幕在线视频| 精品久久久久久,| 熟女人妻精品中文字幕| 亚洲狠狠婷婷综合久久图片| 狂野欧美白嫩少妇大欣赏| 婷婷亚洲欧美| 国产精品久久久av美女十八| 高清毛片免费观看视频网站| 久久久久免费精品人妻一区二区| 中文亚洲av片在线观看爽| 免费观看人在逋| 啦啦啦韩国在线观看视频| 村上凉子中文字幕在线| 这个男人来自地球电影免费观看| 色av中文字幕| 国产高潮美女av| 黄色视频,在线免费观看| 国产精品久久电影中文字幕| 老司机在亚洲福利影院| 少妇人妻一区二区三区视频| 亚洲精华国产精华精| 午夜精品一区二区三区免费看| 国产野战对白在线观看| 久久久久久久精品吃奶| 91老司机精品| 女生性感内裤真人,穿戴方法视频| 免费在线观看成人毛片| 色综合婷婷激情| 日韩 欧美 亚洲 中文字幕| 给我免费播放毛片高清在线观看| 亚洲自偷自拍图片 自拍| 亚洲18禁久久av| 国产激情久久老熟女| 999久久久精品免费观看国产| 精品免费久久久久久久清纯| 日日夜夜操网爽| 91av网一区二区| 亚洲国产欧美一区二区综合| 日日摸夜夜添夜夜添小说| 国产精品av视频在线免费观看| 在线免费观看不下载黄p国产 | 99久久精品热视频| 在线观看免费午夜福利视频| 久久久久国产精品人妻aⅴ院| 欧美黄色片欧美黄色片| 精品无人区乱码1区二区| 国产精品亚洲一级av第二区| 免费看美女性在线毛片视频| av黄色大香蕉| 三级国产精品欧美在线观看 | 小蜜桃在线观看免费完整版高清| 1000部很黄的大片| 天天添夜夜摸| 看免费av毛片| 男女做爰动态图高潮gif福利片| 麻豆成人午夜福利视频| 国产午夜精品论理片| 18禁黄网站禁片免费观看直播| 中国美女看黄片| 亚洲av成人一区二区三| 婷婷六月久久综合丁香| 十八禁人妻一区二区| 欧美最黄视频在线播放免费| 男人舔女人的私密视频| 亚洲aⅴ乱码一区二区在线播放| 国内久久婷婷六月综合欲色啪| 好男人电影高清在线观看| 99久久无色码亚洲精品果冻| 亚洲成av人片在线播放无| 国产成年人精品一区二区| 国产v大片淫在线免费观看| 欧美日韩综合久久久久久 | 长腿黑丝高跟| 午夜激情福利司机影院| 一二三四在线观看免费中文在| 欧美日韩中文字幕国产精品一区二区三区| 人妻久久中文字幕网| 日韩中文字幕欧美一区二区| 99在线人妻在线中文字幕| 夜夜夜夜夜久久久久| 成人高潮视频无遮挡免费网站| 国产高清激情床上av| 在线观看午夜福利视频| 精品国内亚洲2022精品成人| 最好的美女福利视频网| 麻豆av在线久日| 国产探花在线观看一区二区| 欧美一区二区精品小视频在线| 久久久久久大精品| 国产日本99.免费观看| 熟妇人妻久久中文字幕3abv| 99re在线观看精品视频| 国产精品久久视频播放| 国产 一区 欧美 日韩| 一区二区三区高清视频在线| 天天添夜夜摸| 欧美国产日韩亚洲一区| 国产乱人视频| 我要搜黄色片| 真实男女啪啪啪动态图| 日本免费一区二区三区高清不卡| 国产黄色小视频在线观看| 国产综合懂色| 日韩欧美精品v在线| 一本一本综合久久| 老熟妇乱子伦视频在线观看| 欧美黄色片欧美黄色片| 欧美成狂野欧美在线观看| 亚洲一区二区三区不卡视频| 黄色女人牲交| 99久久国产精品久久久| 在线十欧美十亚洲十日本专区| av欧美777| av女优亚洲男人天堂 | 亚洲精品国产精品久久久不卡| 老司机午夜福利在线观看视频| 亚洲电影在线观看av| 俄罗斯特黄特色一大片| 国产人伦9x9x在线观看| 日韩欧美国产一区二区入口| 欧美zozozo另类| 天堂网av新在线| 免费看a级黄色片| 最近最新免费中文字幕在线| 麻豆一二三区av精品| 啦啦啦免费观看视频1| 这个男人来自地球电影免费观看| 精品久久久久久久人妻蜜臀av| 天天躁狠狠躁夜夜躁狠狠躁| 性色av乱码一区二区三区2| 亚洲成人久久爱视频| 特级一级黄色大片| 99久久精品热视频| 99久久无色码亚洲精品果冻| 亚洲 欧美一区二区三区| 午夜影院日韩av| 欧美在线一区亚洲| 精品久久久久久久久久久久久| 亚洲精品在线美女| 19禁男女啪啪无遮挡网站| 很黄的视频免费| 欧美日韩瑟瑟在线播放| 久久亚洲精品不卡| 露出奶头的视频| 成人三级做爰电影| 三级男女做爰猛烈吃奶摸视频| 19禁男女啪啪无遮挡网站| 婷婷丁香在线五月| 婷婷精品国产亚洲av| 久久久水蜜桃国产精品网| 天天躁日日操中文字幕| 久久伊人香网站| 网址你懂的国产日韩在线| 香蕉av资源在线| 成年免费大片在线观看| 在线观看免费视频日本深夜| 美女高潮喷水抽搐中文字幕| aaaaa片日本免费| 一个人看视频在线观看www免费 | 欧美高清成人免费视频www| 亚洲成人中文字幕在线播放| 亚洲成av人片免费观看| 午夜久久久久精精品| 伊人久久大香线蕉亚洲五| 久久久成人免费电影| 99在线人妻在线中文字幕| 成年免费大片在线观看| 夜夜看夜夜爽夜夜摸| 91av网站免费观看| 亚洲欧美日韩高清在线视频| 丁香欧美五月| 久久久国产成人免费| 国产亚洲精品av在线| 国产精品久久久人人做人人爽| 亚洲精华国产精华精| 精品国内亚洲2022精品成人| 香蕉丝袜av| 亚洲欧美日韩高清专用| 熟妇人妻久久中文字幕3abv| 国产欧美日韩一区二区三| 亚洲欧美激情综合另类| 小说图片视频综合网站| 偷拍熟女少妇极品色| 99热这里只有是精品50| 男人舔女人下体高潮全视频| 午夜精品久久久久久毛片777| 亚洲无线观看免费| 国产成年人精品一区二区| 色综合欧美亚洲国产小说| 午夜福利在线观看免费完整高清在 | 午夜福利高清视频| 母亲3免费完整高清在线观看| 免费观看精品视频网站| 一级a爱片免费观看的视频| 美女cb高潮喷水在线观看 | 国产精品自产拍在线观看55亚洲| 亚洲国产中文字幕在线视频| 成年女人看的毛片在线观看| 亚洲九九香蕉| 久久国产精品人妻蜜桃| 日本与韩国留学比较| 一本久久中文字幕| 好看av亚洲va欧美ⅴa在| 91麻豆精品激情在线观看国产| 午夜激情欧美在线| 麻豆成人av在线观看| 国产精品精品国产色婷婷| www.熟女人妻精品国产| 日韩精品青青久久久久久| 久久精品91无色码中文字幕| 麻豆一二三区av精品| 国产伦精品一区二区三区视频9 | 日本精品一区二区三区蜜桃| 女人被狂操c到高潮| 欧美三级亚洲精品| 精品免费久久久久久久清纯| 国内精品久久久久精免费| 美女午夜性视频免费| 在线观看免费午夜福利视频| 国产黄片美女视频| 久久香蕉国产精品| 丰满人妻熟妇乱又伦精品不卡| 久久久久亚洲av毛片大全| 高清毛片免费观看视频网站| 无限看片的www在线观看| 国产精品一区二区三区四区久久| 美女扒开内裤让男人捅视频| 又粗又爽又猛毛片免费看| 国产伦人伦偷精品视频| 成人欧美大片| 五月玫瑰六月丁香| 日本在线视频免费播放| 亚洲国产精品sss在线观看| 午夜激情福利司机影院| 国产极品精品免费视频能看的| 精品国产乱码久久久久久男人| 91老司机精品| 精品国产三级普通话版| 成年免费大片在线观看| 久久久久久九九精品二区国产| 精品一区二区三区视频在线 | av在线蜜桃| 一级作爱视频免费观看| 成人av一区二区三区在线看| 亚洲国产色片| 亚洲色图av天堂| 精品免费久久久久久久清纯| 99久久99久久久精品蜜桃| 国产成人福利小说| 亚洲avbb在线观看| 美女免费视频网站| 成人性生交大片免费视频hd| 夜夜看夜夜爽夜夜摸| 久久天躁狠狠躁夜夜2o2o| 欧美不卡视频在线免费观看| 欧美日韩乱码在线| 黑人操中国人逼视频| 国产黄a三级三级三级人| 黄色丝袜av网址大全| 亚洲国产精品999在线| 国产美女午夜福利| 免费av不卡在线播放| 中文字幕熟女人妻在线| 国产综合懂色| 久久久成人免费电影| 法律面前人人平等表现在哪些方面| 丁香欧美五月| 免费看日本二区| 97人妻精品一区二区三区麻豆| 狂野欧美激情性xxxx| 国产av在哪里看| 午夜激情福利司机影院| 999精品在线视频| 久久久久国内视频| 免费av不卡在线播放| 免费搜索国产男女视频| 成人永久免费在线观看视频| 国产亚洲欧美在线一区二区| 天堂影院成人在线观看| 一级a爱片免费观看的视频| www国产在线视频色| 欧美av亚洲av综合av国产av| 美女高潮的动态| 亚洲国产精品sss在线观看| 国产成人欧美在线观看| 国产三级中文精品| 麻豆国产av国片精品| 成人特级黄色片久久久久久久| 69av精品久久久久久| 国产毛片a区久久久久| 色播亚洲综合网| 国产免费男女视频| 成在线人永久免费视频| 精品久久久久久久人妻蜜臀av| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久久久免费视频| 深夜精品福利| 日韩大尺度精品在线看网址| 黄色日韩在线| 嫩草影院精品99| 午夜日韩欧美国产| 九色成人免费人妻av| 久久伊人香网站| 少妇丰满av| 国产亚洲欧美在线一区二区| 看片在线看免费视频| 亚洲精品国产精品久久久不卡| 99精品在免费线老司机午夜| 久久香蕉国产精品| 欧美日本视频| 观看美女的网站| 免费在线观看视频国产中文字幕亚洲| 国产精品爽爽va在线观看网站| 色综合婷婷激情| 久久中文字幕人妻熟女| 国产伦在线观看视频一区| 免费观看人在逋| 在线十欧美十亚洲十日本专区| 亚洲av熟女| 黄色日韩在线| 18禁裸乳无遮挡免费网站照片| 亚洲中文字幕日韩| 久久久色成人| 日本免费一区二区三区高清不卡| 三级毛片av免费| 国产男靠女视频免费网站| 又爽又黄无遮挡网站| 成人三级黄色视频| 蜜桃久久精品国产亚洲av| 九九久久精品国产亚洲av麻豆 | 亚洲av片天天在线观看| 中文字幕熟女人妻在线| 搡老熟女国产l中国老女人| 叶爱在线成人免费视频播放| 国产成+人综合+亚洲专区| 偷拍熟女少妇极品色| 特级一级黄色大片| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧| 国产精品 国内视频| 日本成人三级电影网站| 国产亚洲精品久久久com|