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

    不同溫度下無后坐炮內(nèi)彈道過程多維隨機變量數(shù)值模擬

    2018-03-31 02:25:37李春雷王雨時張志彪
    彈道學報 2018年1期
    關(guān)鍵詞:包絡(luò)線火藥藥量

    李春雷,王雨時,張志彪

    (南京理工大學 機械工程學院,江蘇 南京 210094)

    無后坐炮內(nèi)彈道過程中隨機影響因素很多。由于內(nèi)彈道整個過程只持續(xù)很短的時間,所以很多影響無后坐炮內(nèi)彈道過程的隨機變量無法通過實驗獲知。

    為了得到內(nèi)彈道影響因素的隨機分布,文獻[1]基于經(jīng)典內(nèi)彈道學模型,采用蒙特卡羅方法模擬了內(nèi)彈道過程中初始參量隨機性導(dǎo)致的最大膛壓和初速的隨機波動現(xiàn)象,得到某高射炮的最大膛壓和初速的隨機分布特性。其中初始參量變化只考慮了藥量、藥厚和火藥力,并未綜合考慮彈丸發(fā)射中其他初始參量,例如彈質(zhì)量和藥室容積。文獻[2]基于內(nèi)彈道勢平衡理論研究了不同火藥靜態(tài)燃燒規(guī)律與初速或然誤差之間的內(nèi)在聯(lián)系,利用建立的初速或然誤差的隨機模擬理論,得到了應(yīng)用密閉爆發(fā)器來預(yù)估初速或然誤差的工程方法。文獻[3]在考慮火藥隨機燃燒基礎(chǔ)上,利用火炮內(nèi)彈道一維兩相流模型,隨機模擬了點火管因隨機破孔而對彈道性能造成的影響;從概率統(tǒng)計的角度對點火管破孔、最大膛壓以及最大負壓差進行了統(tǒng)計分析,統(tǒng)計分析結(jié)果基本能反映出膛內(nèi)射擊過程的某些內(nèi)在統(tǒng)計規(guī)律。文獻[4]對無后坐炮整個身管和出炮口一定距離的流場進行了數(shù)值模擬,得到的流場結(jié)果符合預(yù)期。文獻[5]為了實現(xiàn)無后坐炮在有限空間內(nèi)的安全發(fā)射,設(shè)計了一種含液態(tài)平衡體的新型裝藥結(jié)構(gòu),并進行了某口徑無后坐炮內(nèi)彈道性能試驗,建立了相應(yīng)的內(nèi)彈道模型,計算結(jié)果與試驗結(jié)果吻合較好。在此基礎(chǔ)上,分析了液態(tài)平衡體初始質(zhì)量、密度、火藥弧厚、噴管喉部直徑以及擠進壓力等參數(shù)對內(nèi)彈道性能的影響。

    綜上所述,蒙特卡洛方法解決了在實際中不可能在相同條件下做大量試驗的突出問題,是深入系統(tǒng)了解問題分布的有效方法。因此,為了精確設(shè)計身管強度、彈體強度和精確分析引信解除保險性能,在考慮藥厚、火藥力、藥量、藥室容積和彈丸質(zhì)量隨機性基礎(chǔ)上,應(yīng)用蒙特卡洛方法和無后坐炮經(jīng)典內(nèi)彈道模型,隨機模擬了不同隨機因素對彈道性能的影響。本文在不同溫度下同時考慮藥厚、火藥力、藥量、藥室容積和彈丸質(zhì)量隨機性,得到了上述初始參量隨機影響因素下的某無后坐炮內(nèi)彈道膛壓曲線包絡(luò)線,該解算結(jié)果能較好地反映出膛內(nèi)射擊過程的變化規(guī)律,可用于精細設(shè)計身管強度、彈體強度和精確分析引信解除保險性能。

    1 隨機模擬方法與無后坐炮內(nèi)彈道

    1.1 隨機模擬方法與無后坐炮內(nèi)彈道的特點

    隨機模擬方法又稱蒙特卡羅(Monte Carlo)方法,是一種通過設(shè)定隨機過程,反復(fù)生成時間序列,計算參數(shù)估計量或統(tǒng)計量,進而研究其分布特征的方法。其詳細定義見文獻[6-7]。

    與一般火炮相比,無后坐炮在射擊過程中有大量火藥氣體從噴管中流出。因氣體流出影響的內(nèi)彈道性能參數(shù)分別有流量qm、總量Qm和推力F。在一維等熵條件下,有:

    (1)

    (2)

    F=CFSjp

    (3)

    無后坐炮內(nèi)彈道方程組由1個代數(shù)方程和6個一階微分方程組成。在此只需得到膛壓曲線,因此方程組中關(guān)于速度的方程可不予考慮。并且一般情況下,該方程組不存在解析解,因此,通常采用數(shù)值方法求解。

    為了便于Matlab軟件編制,將其化為無量綱的微分方程組形式,具體形式詳見文獻[8]。對其采用Matlab軟件自身提供的求常微分方程數(shù)值解的四階、五階的ode45函數(shù)求解。

    1.2 仿真可信性說明

    根據(jù)文獻[9]提供的1965年式82 mm無后坐炮的內(nèi)彈道初始數(shù)據(jù),利用編寫好的程序?qū)ζ鋬?nèi)彈道進行仿真,結(jié)果如表1所示,表中,l為彈丸行程,p為膛壓,ε為仿真結(jié)果與文獻[9]中結(jié)果之間的相對誤差。

    表1 常溫下(15 ℃)1965式82 mm無后坐炮身管內(nèi)彈道仿真結(jié)果

    由表1可知,上述內(nèi)彈道膛壓解算結(jié)果與文獻[9]計算得到的內(nèi)彈道膛壓結(jié)算結(jié)果相對誤差很小,因此本文所用的內(nèi)彈道仿真模型、軟件和解算結(jié)果是可信的。

    2 內(nèi)彈道膛壓隨機模擬

    2.1 彈質(zhì)量隨機模擬

    首先假設(shè)無后坐炮藥室容積、發(fā)射藥藥量、發(fā)射藥藥厚和發(fā)射藥火藥力恒定不變,環(huán)境溫度即藥溫為常溫(15 ℃)。根據(jù)文獻[10],彈質(zhì)量隨機變化,服從正態(tài)分布。由文獻[11],1個計算組內(nèi)不同彈質(zhì)量偏差一般不超過1個彈質(zhì)量分級符號。根據(jù)文獻[12],一般榴彈的彈質(zhì)量符號在4個以內(nèi),故在此取4個適中符號,其中彈質(zhì)量符號和適中符號的具體定義見文獻[12];中間誤差[13]Em代表隨機變量出現(xiàn)在均值左右Em范圍內(nèi)的概率為50%,在均值左右4Em范圍內(nèi)隨機變量出現(xiàn)的概率近似為1;由此即可估算出彈丸質(zhì)量分布的標準差σmp=0.009 884mp,其中,mp為彈丸質(zhì)量。

    由Monte-Carlo方法產(chǎn)生5個彈質(zhì)量隨機值,然后循環(huán)模擬內(nèi)彈道過程5次,得到的膛壓曲線如圖1所示。

    圖1 隨機模擬彈質(zhì)量獲得的膛壓曲線

    從圖1知,在5次隨機模擬后,圖示局部放大區(qū)域曲線有交叉,因此隨機模擬彈質(zhì)量所得到的膛壓曲線包絡(luò)線不是彈質(zhì)量中最大與最小所對應(yīng)的膛壓曲線。進一步由Monte-Carlo方法產(chǎn)生105個彈質(zhì)量隨機值,然后循環(huán)模擬內(nèi)彈道過程105次,通過比較某時間(某行程)對應(yīng)的壓力值,得到105條膛壓曲線的包絡(luò)線,如圖2所示。

    圖2 隨機模擬彈質(zhì)量對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖2中,最大膛壓附近的膛壓跳動最大值為0.775 3 MPa,炮口附近的膛壓跳動最大值為5.526 3 MPa。

    2.2 發(fā)射藥藥厚隨機模擬

    首先假設(shè)無后坐炮藥室容積、彈質(zhì)量、發(fā)射藥藥量和發(fā)射藥火藥力是恒定不變的,環(huán)境溫度即藥溫為常溫(15 ℃)。再假設(shè)發(fā)射藥藥厚隨機變化,服從正態(tài)分布。由文獻[14],全彈裝配時,所有尺寸均符合正態(tài)分布,其散布中心即為公差帶中心,散布范圍6σe即為公差帶寬度W。由于正態(tài)分布只有99.73%的取值在公差帶寬度W內(nèi),故要剔除在公差帶W之外的抽樣尺寸即不合格的尺寸。根據(jù)文獻[15],某無后坐炮發(fā)射藥藥厚制造公差帶寬度為0.16 mm,因此可估算出發(fā)射藥藥厚的標準差σe約為(0.16/6) mm=0.026 67 mm。

    利用2.1的方法可知,隨機模擬藥厚所得到的膛壓曲線包絡(luò)線也并不是一組隨機發(fā)射藥藥厚數(shù)值中最大與最小所對應(yīng)的膛壓曲線。進一步由Monte-Carlo方法產(chǎn)生105個發(fā)射藥藥厚隨機值,然后循環(huán)模擬內(nèi)彈道過程105次,通過比較某時間(某行程)對應(yīng)的壓力值,得到105條膛壓曲線的包絡(luò)線,如圖3所示。

    圖3 發(fā)射藥藥厚隨機值對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖3中,最大膛壓附近的膛壓跳動最大值為5.247 3 MPa,炮口附近的膛壓跳動最大值為10.785 1 MPa。

    2.3 發(fā)射藥藥量隨機模擬

    首先假設(shè)無后坐炮藥室容積、彈質(zhì)量、發(fā)射藥藥厚和發(fā)射藥火藥力恒定不變,環(huán)境溫度即藥溫為常溫(15 ℃),發(fā)射藥藥量隨機變化,服從正態(tài)分布。全彈裝配時,100 mm和105 mm口徑彈的發(fā)射藥藥量公差為±(4~6)g,由天平精度保證,因此某無后坐炮的發(fā)射藥藥量公差估取為±(1~3) g,則發(fā)射藥藥量標準差σmw分別取為0.33 g,0.67 g,1 g。

    利用2.1的方法可知,隨機模擬發(fā)射藥藥量所得到的膛壓曲線包絡(luò)線并不是一組隨機發(fā)射藥藥量中最大與最小所對應(yīng)的膛壓曲線。進一步由Monte-Carlo方法產(chǎn)生105個發(fā)射藥藥量隨機值,然后循環(huán)模擬內(nèi)彈道過程105次,通過比較某時間(某行程)對應(yīng)的壓力值,得到105條膛壓曲線的包絡(luò)線。發(fā)射藥藥量標準差為0.33 g,0.67 g,1 g,所對應(yīng)的包絡(luò)線如圖4~圖6所示。

    圖4 發(fā)射藥藥量標準差σmw=0.33 g時隨機發(fā)射藥藥量對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖5 發(fā)射藥藥量標準差σmw=0.67 g時隨機發(fā)射藥藥量對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖6 發(fā)射藥藥量標準差σmw=1 g時隨機發(fā)射藥藥量對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    由圖4~圖6可知,最大膛壓附近的膛壓跳動最大值和炮口附近的膛壓跳動最大值,在σmw=0.33 g時分別為0.429 7 MPa,0.143 4 MPa;在σmw=0.67 g時分別為0.841 4 MPa,0.441 7 MPa;在σmw=1 g時分別為1.255 7 MPa,0.523 6 MPa。

    2.4 發(fā)射藥火藥力隨機模擬

    首先假設(shè)無后坐炮發(fā)射藥藥厚、彈質(zhì)量、發(fā)射藥藥量和藥室容積恒定不變,環(huán)境溫度即藥溫為常溫(15 ℃),發(fā)射藥火藥力隨機變化,服從正態(tài)分布。據(jù)靶場統(tǒng)計可知,一批發(fā)射藥的火藥力的變動范圍在±3%以內(nèi),因此可估計發(fā)射藥火藥力的標準差σf=1%。

    利用2.1的方法可知,隨機模擬火藥力所得到的膛壓曲線包絡(luò)線并不是一組隨機火藥力中最大與最小所對應(yīng)的膛壓曲線。進一步由Monte-Carlo方法產(chǎn)生105個發(fā)射藥火藥力隨機值,然后循環(huán)模擬內(nèi)彈道過程105次,通過比較某時間(某行程)對應(yīng)的壓力值,得到105條膛壓曲線的最小輪廓和最大輪廓,如圖7所示。

    圖7 隨機發(fā)射藥火藥力對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖7中,最大膛壓附近的膛壓跳動最大值為4.490 6 MPa,炮口附近的膛壓跳動最大值為1.210 1 MPa。

    2.5 藥室容積隨機模擬

    2.5.1 藥室容積的制造誤差

    首先假設(shè)無后坐炮發(fā)射藥藥厚、彈質(zhì)量、發(fā)射藥藥量和發(fā)射藥火藥力恒定不變,環(huán)境溫度即藥溫為常溫(15 ℃),藥室容積隨機變化,服從正態(tài)分布。無后坐炮藥室直徑和藥室長度制造公差帶寬度分別為0.36 mm和0.52 mm,2.2節(jié)已介紹全彈尺寸公差帶與尺寸標準差的關(guān)系,同理可估算出某無后坐炮藥室直徑和藥室長度的標準差分別約為0.06 mm和0.09 mm,因此可得到隨機的藥室容積。

    利用2.1的方法可知,隨機模擬藥室容積所得到的膛壓曲線包絡(luò)線并不是一組隨機藥室容積值中最大與最小所對應(yīng)的膛壓曲線。利用文獻[1]的方法產(chǎn)生105個隨機藥室直徑值和隨機藥室長度值,然后由藥室容積的計算公式計算得到105個藥室容積隨機值,然后循環(huán)模擬內(nèi)彈道過程105次,通過比較某時間(某行程)對應(yīng)的壓力值,得到105條膛壓曲線的包絡(luò)線,如圖8所示。

    圖8 隨機藥室容積對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖8中,最大膛壓附近的膛壓跳動最大值為0.537 4 MPa,炮口附近的膛壓跳動最大值為0.413 7 MPa。制造公差控制良好,藥室容積的隨機變化對無后坐炮內(nèi)彈道膛壓曲線的分布區(qū)域影響很小,可忽略不計。

    2.5.2 藥室容積磨損

    藥室容積變化包括藥室長度與藥室內(nèi)徑的變化,連續(xù)發(fā)射多發(fā)炮彈時,火藥氣體對內(nèi)膛不斷地燒蝕沖刷,使火炮的藥室容積不斷增大。根據(jù)文獻[8]給出的130 mm加農(nóng)炮內(nèi)徑磨損量與射彈數(shù)的變化關(guān)系以及藥室長度增長量與射彈數(shù)的變化關(guān)系,用回歸的方法可得內(nèi)徑磨損量Δd與射彈發(fā)數(shù)N的關(guān)系式以及藥室長度增長量ΔL與射彈發(fā)數(shù)N的關(guān)系式。

    藥室長度增長量ΔL與射彈發(fā)數(shù)N的關(guān)系式:

    ΔL=5.478×10-8N3-7.692×10-6N2+0.015 75N-0.174 8

    (4)

    內(nèi)徑磨損量Δd與射彈發(fā)數(shù)N的關(guān)系式:

    Δd=1.936×10-9N3-6.171×10-6N2+0.010 3N+0.119 9

    (5)

    把回歸結(jié)果和實驗數(shù)據(jù)繪成曲線,如圖9(a)和圖9(b)所示。從圖中可以看出:離散點與回歸曲線符合很好;此外由相關(guān)系數(shù)可知,圖9(a)的相關(guān)系數(shù)ra=0.999 8,圖9(b)的相關(guān)系數(shù)rb=0.999 1,即擬合曲線很接近真實情況。

    圖9 藥室長度增長量、內(nèi)徑磨損量與射彈數(shù)關(guān)系回歸曲線

    文獻[16]提到瑞典M3無后坐炮的設(shè)計壽命為500發(fā),但是美軍靶場測試中曾經(jīng)發(fā)射2 360發(fā)炮彈仍然完好。據(jù)此,將某無后坐炮的設(shè)計壽命取為1 000發(fā)??筛鶕?jù)式(7)和式(8)計算出130 mm加農(nóng)炮在發(fā)射1 000發(fā)射彈后的內(nèi)徑磨損量Δd=5.43 mm和藥室長度增長量ΔL=62.66 mm。在沒有無后坐炮磨損量實驗數(shù)據(jù)的背景下,現(xiàn)假設(shè)在相同射擊炮彈數(shù)下,磨損量與最大膛壓成正比。因此,可得某無后坐炮的最大設(shè)計內(nèi)徑磨損量Δdmax=0.62 mm和最大設(shè)計藥室長度增長量ΔLmax=8.00 mm,進而由藥室容積計算公式得到最大設(shè)計藥室容積磨損量。由于在105次隨機模擬下,無后坐炮藥室容積已磨損到極限,因此,可取最大設(shè)計藥室磨損量下的藥室容積進行內(nèi)彈道隨機模擬。

    2.6 分析與討論

    綜上所述,發(fā)射藥藥厚的隨機變化對膛壓曲線分布區(qū)域的影響最大,且在整個內(nèi)彈道過程中影響整個身管的膛壓,其中最大膛壓附近的膛壓跳動最大值為5.247 3 MPa;發(fā)射藥藥量的隨機變化對無后坐炮內(nèi)彈道膛壓曲線的分布區(qū)域影響較小,且發(fā)射藥藥量的隨機變化影響整個內(nèi)彈道過程的膛壓,發(fā)射藥藥量標準差σmw=1 g時,其最大膛壓附近的膛壓跳動最大值為1.255 7 MPa;隨機彈質(zhì)量對無后坐炮內(nèi)彈道膛壓曲線的分布區(qū)域影響較小,且隨機彈質(zhì)量主要影響炮口附近壓力值,其最大膛壓附近的膛壓跳動最大值為0.775 3 MPa;隨機發(fā)射藥火藥力對無后坐炮內(nèi)彈道膛壓曲線的分布區(qū)域影響較小,發(fā)射藥火藥力的變化影響整個身管的膛壓,且對最大膛壓附近的影響相對較大,其最大膛壓附近的膛壓跳動最大值為4.490 6 MPa;隨機藥室容積(不考慮藥室使用磨損)的隨機變化對無后坐炮內(nèi)彈道膛壓曲線的分布區(qū)域影響很小,可忽略不計,但是必須考慮在多次隨機模擬下藥室容積的擴張,此時應(yīng)取最大設(shè)計磨損量下的藥室容積進行解算。

    3 不同溫度下4種裝填參量隨機模擬

    由文獻[17]知,火炮在作戰(zhàn)條件下使用時,裝藥溫度受環(huán)境溫度影響很大。為了保證安全,在設(shè)計身管時,膛內(nèi)壓力必須考慮裝藥溫度的變化。目前,常用的溫度設(shè)定是:常溫為15 ℃,高溫為+50 ℃,低溫為-40 ℃。

    設(shè)裝藥溫度為θ,裝藥溫度改變量Δθ=(θ-15) ℃。當θ=50 ℃時,Δθ=35 ℃;當θ=-40 ℃時,Δθ=-55 ℃。裝藥初溫變化,最大膛壓也要變化,其變化值為

    式中:cθ為最大壓力修正系數(shù)。

    初步計算時可以采用如下的關(guān)系式:硝化棉系火藥,cθ=0.002 7cθk;硝化甘油系火藥,cθ=0.003 5cθk。其中,cθk為壓力全沖量修正系數(shù),隨裝填密度Δ及最大壓力pc的不同而變化,查閱文獻[9]可得到系數(shù)cθk的值。

    由此可計算出“溫度影響系數(shù)”:J=1+cθΔθ。

    只要在第2節(jié)所提的內(nèi)彈道的燃速方程中乘以相對應(yīng)的“溫度影響系數(shù)”,就可得相對應(yīng)溫度下的膛壓曲線。

    綜上所述,單個隨機因素對內(nèi)彈道膛壓曲線包絡(luò)線影響較大的有彈質(zhì)量 、藥量、火藥力和藥厚。

    現(xiàn)取最大設(shè)計藥室磨損量下的藥室容積,并在高溫(50 ℃)、常溫(15 ℃)和低溫(-40 ℃)下,隨機模擬各影響因素,得到膛壓曲線包絡(luò)線。

    根據(jù)文獻[1]可知,影響膛內(nèi)射擊過程的因素是一個多維的問題。為了簡單起見,假定藥溫為15 ℃,不考慮點火過程的隨機因素,考慮彈質(zhì)量、裝藥量(σmw=1 g)、發(fā)射藥藥厚和發(fā)射藥火藥力等4個裝填參量的隨機變化。

    四維正態(tài)隨機變量(Y1Y2Y3Y4)的協(xié)方差矩陣為

    (6)

    式中:Kij=E[(Yi-μi)(Yj-μj)]。

    (7)

    Σ=CCT

    (8)

    由式(8)得:

    (9)

    隨機變量Y1,Y2,Y3,Y4分別代表彈質(zhì)量、發(fā)射藥火藥力、發(fā)射藥藥量和發(fā)射藥藥厚的隨機值。根據(jù)它們已知的分布可抽樣一組隨機值,帶入內(nèi)彈道方程組,得到一條膛壓曲線,如此循環(huán),可得大量內(nèi)彈道膛壓曲線,并得到其包絡(luò)線。

    高溫、常溫和低溫下的解算結(jié)果如圖10~圖12所示。

    圖10 高溫(50 ℃)下4種隨機因素對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖11 常溫(15℃)下4種隨機因素對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    圖12 低溫(-40 ℃)下4種隨機因素對某無后坐炮內(nèi)彈道膛壓曲線分布區(qū)域的影響

    由圖10~圖12可知,最大膛壓附近的膛壓跳動最大值和炮口附近的膛壓跳動最大值,在高溫時分別為14.088 8 MPa,6.352 6 MPa;在常溫時分別為7.749 8 MPa,10.832 2 MPa;在低溫時分別為5.660 3 MPa,2.614 7 MPa。

    4 結(jié)論

    在高溫(50 ℃)、常溫(15 ℃)和低溫(-40 ℃)下,考慮藥厚、火藥力、藥量和彈質(zhì)量4個隨機因素同時對某無后坐炮內(nèi)彈道膛壓曲線輪廓的影響,隨機模擬各影響因素,得到膛壓曲線包絡(luò)線。在不同溫度下,4個隨機因素對無后坐炮內(nèi)彈道膛壓曲線的分布區(qū)域影響都較大,且影響整個身管和內(nèi)彈道的膛壓分布,同時比單個隨機因素對無后坐炮內(nèi)彈道膛壓曲線的分布區(qū)域影響更為明顯,其中高溫、低溫和常溫下最大膛壓附近的膛壓跳動最大值分別為14.088 8 MPa,5.660 3 MPa,7.749 8 MPa。因此,高溫下膛壓曲線的上包絡(luò)線可用來對某無后坐炮身管和彈體進行強度校核,低溫下膛壓曲線的下包絡(luò)線可用來分析某無后坐炮引信的解除保險性能。

    本文在考慮藥厚、火藥力、藥量和彈質(zhì)量的隨機性基礎(chǔ)上,應(yīng)用蒙特卡洛方法,結(jié)合無后坐炮內(nèi)彈道經(jīng)典模型,對不同隨機因素對彈道性能造成的影響進行了隨機模擬,計算結(jié)果能較好地反映出膛內(nèi)射擊過程的變化規(guī)律,并對某無后坐炮的系統(tǒng)設(shè)計有一定的指導(dǎo)意義。同時,在制造生產(chǎn)無后坐炮時,必須將發(fā)射藥和彈丸的制造誤差控制在令人滿意的范圍內(nèi),這樣某后坐炮的膛壓跳動會相對較小,有利于提高某無后坐炮的射擊精度和使用壽命。

    [1] 金志明. 內(nèi)彈道循環(huán)的隨機模擬[J]. 南京理工大學學報,1987(3):42-51.

    JIN Zhiming. Stochastic simulation of inner trajectory cycle[J]. Journal of Nanjing University of Science and Technology,1987(3):42-51. (in Chinese)

    [2] 鮑廷鈺,張兆鈞,邱文堅,等. 火藥不均一性引起初速散布的密閉爆發(fā)器試驗預(yù)估[J]. 南京理工大學學報,1994(5):25-32.

    BAO Tingyu,ZHANG Zhaojun,QIU Wenjian,et al. Experimental prediction of closed bursts caused by the nonuniformity of gun powder[J]. Journal of Nanjing University of Science and Technology,1994(5):25-32. (in Chinese)

    [3] 馮德成,翁春生,白橋棟,等. 隨機因素對彈道性能影響的數(shù)值仿真[J]. 火炮發(fā)射與控制學報,2005(1):1-5.

    FENG Decheng,WENG Chunsheng,BAI Qiaodong,et al. Numerical simulation of the effect of stochastic factors on trajectory performance[J]. Journal of Gun Launch & Control,2005(1):1-5. (in Chinese)

    [4] NEVES A,MENDES E,SCHILLER L,et al. Recoilless gun:analysis of internal and external ballistics[C]//The 43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. Cincinnati:AIAA,2007:5160.

    [5] 曹永杰,張向明,李志飛. 基于液態(tài)平衡體的無后坐炮內(nèi)彈道研究[J]. 彈道學報,2016,28(3):71-75.

    CAO Yongjie,ZHANG Xiangming,LI Zhifei. Research on interior ballistics of recoil gun based on liquid equilibrium[J]. Journal of Ballistics,2016,28(3):71-75. (in Chinese)

    [6] 方再根. 計算機模擬和蒙特卡洛方法[M]. 北京:北京工業(yè)學院出版社,1988.

    FANG Zaigen. Computer simulation and Monte Carlo method[M]. Beijing:Beijing University of Technology Press,1988. (in Chinese)

    [7] 張晉華,聞泉,王雨時,等. 應(yīng)用蒙特卡羅方法計算彈丸偏心距[J]. 探測與控制學報,2016,38(5):42-48.

    ZHANG Jinhua,WEN Quan,WANG Yushi,et al. Application of Monte Carlo method to calculate projectile eccentricity[J]. Journal of Detection & Control,2016,38(5):42-48. (in Chinese)

    [8] 張小兵. 槍炮內(nèi)彈道學[M]. 北京:北京理工大學出版社,2014.

    ZHANG Xiaobing. Interior ballistics of guns[M]. Beijing:Beijing Institute of Technology Press,2014. (in Chinese)

    [9] 王連榮,張佩勤. 火炮內(nèi)彈道計算手冊[M]. 北京:國防工業(yè)出版社,1987.

    WANG Lianrong,ZHANG Peiqin. Manual of artillery interior ballistics calculation[M]. Beijing:National Defense Industry Press,1987. (in Chinese)

    [10] 王曉鵬,王雨時,盧鳳生,等. 155 mm口徑火炮榴彈結(jié)構(gòu)特征數(shù)分布特性研究[J]. 探測與控制學報,2015,37(5):66-72.

    WANG Xiaopeng,WANG Yushi,LU Fengsheng,et al. Study on the characteristic number distribution characteristics of 155mm caliber artillery shrapnel[J]. Journal of Detection & Control,2015,37(5):66-72. (in Chinese)

    [11] 中國人民解放軍總裝備部. 榴彈定型試驗規(guī)程:GJB 4225—2001[S]. 北京:總裝備部軍標出版發(fā)行部,2001.

    General Armament Department of Chinese People’s Liberation Army. Grenade shaping test procedure:GJB 4225—2001[S]. Beijing:General Armament Depentment’s Military Standard Publication and Distribution Department,2001. (in Chinese)

    [12] 華恭,歐林爾. 彈丸作用和設(shè)計理論[M]. 北京:國防工業(yè)出版社,1975.

    HUA Gong,OU Liner. Projectile function and design theory[M]. Beijing:National Defense Industry Press,1975. (in Chinese)

    [13] 韓子鵬. 彈箭外彈道學[M]. 北京:北京理工大學出版社,2014.

    HAN Zipeng. Rocket and projectile exterior ballistics[M]. Beijing:Beijing Institute of Technology Press,2014. (in Chinese)

    [14] 聞泉,王雨時. 子母彈徑向質(zhì)心位置的蒙特卡羅模擬[J]. 彈箭與制導(dǎo)學報,2006,26(增刊6):403-405.

    WEN Quan,WANG Yushi. A Monte Carlo simulation of radial centroid position of submunition[J]. Journal of Projectiles,Rockets,Missiles and Guidance,2006,26(S6):403-405. (in Chinese)

    [15] 李福平,等. 《火炸藥手冊》(Ⅱ)[M]. 西安:第五機械工業(yè)部第204研究所,1981.

    LI Fuping,et al. Fire and explosives manual(Ⅱ)[M]. Xi’an:No. 204 Institute of the Fifth Machinery Industry Institute Department,1981. (in Chinese)

    [16] 漢威. 步兵班的狙擊炮——“卡爾·古斯塔夫”無后坐力炮[J]. 兵器知識,2015(8):52-57.

    HAN Wei. Infantry class sniper gun:“Carl Gustav” recoilless guns[J]. Weapon Knowledge,2015(8):52-57. (in Chinese)

    [17] 張相炎,鄭建國,楊軍榮. 火炮設(shè)計理論[J]. 北京:北京理工大學出版社,2005.

    ZHANG Xiangyan,ZHENG Jianguo,YANG Junrong. Artillery design theory[J]. Beijing:Beijing Institute of Technology Press,2005. (in Chinese)

    猜你喜歡
    包絡(luò)線火藥藥量
    平面曲線運動中的包絡(luò)線問題探討
    神奇的火藥
    大灰狼畫報(2022年2期)2022-05-23 12:57:10
    碳化對飛灰穩(wěn)定化加藥量測評的影響
    兩倍藥量
    火藥的來歷
    拋體的包絡(luò)線方程的推導(dǎo)
    “火藥弟弟”
    一種用于故障隔離的參數(shù)區(qū)間包絡(luò)線計算方法
    小型水池抗爆炸沖擊極限藥量分析
    藥量適中
    小說月刊(2015年6期)2015-12-16 14:55:45
    啦啦啦中文免费视频观看日本| 国产成人av教育| 国产亚洲精品久久久久5区| 日韩成人在线观看一区二区三区| 亚洲国产成人一精品久久久| 午夜两性在线视频| 国产精品久久久久成人av| 黄片播放在线免费| 久久久精品94久久精品| 国产单亲对白刺激| 老司机午夜十八禁免费视频| 久久久国产成人免费| 成年人黄色毛片网站| 亚洲人成电影观看| 一级片免费观看大全| cao死你这个sao货| 纵有疾风起免费观看全集完整版| videosex国产| 欧美变态另类bdsm刘玥| 国产精品98久久久久久宅男小说| 国产精品欧美亚洲77777| 777米奇影视久久| 久久影院123| 欧美久久黑人一区二区| 曰老女人黄片| 欧美日韩黄片免| 成人av一区二区三区在线看| 国产一区二区三区视频了| 日韩 欧美 亚洲 中文字幕| av视频免费观看在线观看| 久久国产精品大桥未久av| 18在线观看网站| 18禁裸乳无遮挡动漫免费视频| 人人妻人人添人人爽欧美一区卜| 国产精品免费视频内射| 国产精品麻豆人妻色哟哟久久| 亚洲成av片中文字幕在线观看| 男女免费视频国产| 女人精品久久久久毛片| 精品国产乱子伦一区二区三区| 一本大道久久a久久精品| 三上悠亚av全集在线观看| 午夜福利一区二区在线看| 大陆偷拍与自拍| 一区二区三区精品91| 一本色道久久久久久精品综合| 亚洲 欧美一区二区三区| 手机成人av网站| 亚洲av成人不卡在线观看播放网| 美女高潮喷水抽搐中文字幕| 国产福利在线免费观看视频| 极品人妻少妇av视频| 国产一区二区 视频在线| 精品一区二区三区av网在线观看 | 深夜精品福利| 无遮挡黄片免费观看| 色在线成人网| 亚洲人成77777在线视频| 国产成人精品无人区| 国产精品一区二区在线观看99| 大香蕉久久网| 久久性视频一级片| √禁漫天堂资源中文www| 色综合婷婷激情| 国产日韩欧美视频二区| 欧美午夜高清在线| cao死你这个sao货| 法律面前人人平等表现在哪些方面| 大片免费播放器 马上看| 亚洲中文字幕日韩| 国产淫语在线视频| 久久久久精品人妻al黑| 国产成人精品无人区| 免费不卡黄色视频| 王馨瑶露胸无遮挡在线观看| 亚洲全国av大片| 国产亚洲精品第一综合不卡| 在线亚洲精品国产二区图片欧美| 激情视频va一区二区三区| 亚洲熟女毛片儿| 国产亚洲av高清不卡| 少妇 在线观看| 国产精品99久久99久久久不卡| 午夜日韩欧美国产| 国产成人av激情在线播放| 亚洲国产欧美网| 黑人操中国人逼视频| 高潮久久久久久久久久久不卡| svipshipincom国产片| 亚洲精品av麻豆狂野| 国产精品熟女久久久久浪| 久久久久久人人人人人| 在线观看一区二区三区激情| 久久久精品免费免费高清| 大香蕉久久网| 视频在线观看一区二区三区| 精品国产一区二区久久| 成人影院久久| 丝袜美足系列| 国产精品99久久99久久久不卡| 韩国精品一区二区三区| 国产主播在线观看一区二区| 80岁老熟妇乱子伦牲交| 欧美老熟妇乱子伦牲交| 亚洲av美国av| 91成年电影在线观看| 亚洲性夜色夜夜综合| 性色av乱码一区二区三区2| 国产精品久久久久久人妻精品电影 | 中文字幕人妻丝袜一区二区| 国产亚洲欧美在线一区二区| 精品国产一区二区三区久久久樱花| 中文字幕人妻丝袜制服| 久久久久久久精品吃奶| 免费观看人在逋| 啪啪无遮挡十八禁网站| 久久九九热精品免费| 在线观看免费午夜福利视频| 一本—道久久a久久精品蜜桃钙片| 别揉我奶头~嗯~啊~动态视频| 欧美在线一区亚洲| 老鸭窝网址在线观看| 最新的欧美精品一区二区| 91成年电影在线观看| 免费观看人在逋| 91麻豆精品激情在线观看国产 | 欧美另类亚洲清纯唯美| 正在播放国产对白刺激| 啦啦啦在线免费观看视频4| 一区二区三区国产精品乱码| 在线 av 中文字幕| 建设人人有责人人尽责人人享有的| 纯流量卡能插随身wifi吗| av天堂在线播放| h视频一区二区三区| 久久久水蜜桃国产精品网| 51午夜福利影视在线观看| 亚洲av欧美aⅴ国产| 久久香蕉激情| 国产精品秋霞免费鲁丝片| 成人三级做爰电影| 国产精品一区二区免费欧美| 9191精品国产免费久久| 老司机在亚洲福利影院| 老熟妇仑乱视频hdxx| 美女午夜性视频免费| 亚洲男人天堂网一区| 欧美精品一区二区大全| 国产成人av教育| 免费久久久久久久精品成人欧美视频| 高潮久久久久久久久久久不卡| 王馨瑶露胸无遮挡在线观看| 色老头精品视频在线观看| 黄色视频在线播放观看不卡| 伊人久久大香线蕉亚洲五| 人妻一区二区av| 久久精品成人免费网站| 男女高潮啪啪啪动态图| 激情视频va一区二区三区| 丁香六月天网| 国产又色又爽无遮挡免费看| a级片在线免费高清观看视频| 91成人精品电影| 国产精品二区激情视频| 人人妻人人澡人人爽人人夜夜| 国产有黄有色有爽视频| 日本vs欧美在线观看视频| svipshipincom国产片| 正在播放国产对白刺激| 日日爽夜夜爽网站| 热99国产精品久久久久久7| 一区二区三区激情视频| 久久久国产欧美日韩av| 国产欧美日韩一区二区三| 国产精品自产拍在线观看55亚洲 | 国产精品久久久久久精品电影小说| 999精品在线视频| 国产成人系列免费观看| 亚洲av美国av| 黄频高清免费视频| 亚洲熟女精品中文字幕| 亚洲成国产人片在线观看| 国产一区二区 视频在线| 天天躁日日躁夜夜躁夜夜| 一边摸一边做爽爽视频免费| 免费看a级黄色片| 国产一卡二卡三卡精品| 中文字幕另类日韩欧美亚洲嫩草| 老熟女久久久| 两个人免费观看高清视频| 日韩欧美三级三区| 俄罗斯特黄特色一大片| 欧美日韩亚洲综合一区二区三区_| 十八禁人妻一区二区| 三级毛片av免费| 久久影院123| 99久久99久久久精品蜜桃| 久久人人爽av亚洲精品天堂| 欧美乱码精品一区二区三区| 欧美成人午夜精品| 无遮挡黄片免费观看| 一二三四社区在线视频社区8| 国产伦人伦偷精品视频| 1024视频免费在线观看| 精品一区二区三区四区五区乱码| 久久久久精品国产欧美久久久| 国产99久久九九免费精品| 老司机深夜福利视频在线观看| 美女视频免费永久观看网站| 美女午夜性视频免费| 人人妻人人澡人人爽人人夜夜| 人成视频在线观看免费观看| 制服诱惑二区| 亚洲国产欧美在线一区| 国产av国产精品国产| 手机成人av网站| 亚洲精品在线观看二区| 美女高潮喷水抽搐中文字幕| 国产福利在线免费观看视频| 午夜老司机福利片| 一级毛片电影观看| 在线观看舔阴道视频| 搡老熟女国产l中国老女人| 亚洲精品在线美女| 精品国产一区二区三区四区第35| 成年人免费黄色播放视频| 男女免费视频国产| 视频区图区小说| 成年动漫av网址| av视频免费观看在线观看| 国内毛片毛片毛片毛片毛片| 免费在线观看黄色视频的| 看免费av毛片| 欧美亚洲日本最大视频资源| 亚洲精品美女久久久久99蜜臀| 99riav亚洲国产免费| 亚洲精品国产色婷婷电影| 精品少妇一区二区三区视频日本电影| 国产真人三级小视频在线观看| 天天操日日干夜夜撸| 亚洲情色 制服丝袜| 在线观看免费视频日本深夜| 少妇被粗大的猛进出69影院| 久久亚洲真实| 欧美中文综合在线视频| 97人妻天天添夜夜摸| avwww免费| 夜夜爽天天搞| av天堂久久9| 少妇 在线观看| 亚洲视频免费观看视频| 在线观看免费日韩欧美大片| 搡老乐熟女国产| 国产精品亚洲一级av第二区| 日韩欧美三级三区| 香蕉国产在线看| 怎么达到女性高潮| 少妇粗大呻吟视频| 美女高潮到喷水免费观看| bbb黄色大片| 中文字幕av电影在线播放| 12—13女人毛片做爰片一| 视频区欧美日本亚洲| 极品人妻少妇av视频| 咕卡用的链子| 欧美av亚洲av综合av国产av| 天堂中文最新版在线下载| 人妻久久中文字幕网| 下体分泌物呈黄色| 欧美黑人精品巨大| 日韩人妻精品一区2区三区| 亚洲综合色网址| 精品福利观看| 1024香蕉在线观看| 十八禁人妻一区二区| 久久99热这里只频精品6学生| 不卡一级毛片| 成年人免费黄色播放视频| 黄色片一级片一级黄色片| 9191精品国产免费久久| 丁香六月天网| 久久久久久久精品吃奶| 大片电影免费在线观看免费| 在线 av 中文字幕| 新久久久久国产一级毛片| 亚洲性夜色夜夜综合| 操美女的视频在线观看| 中文字幕高清在线视频| 18禁观看日本| 91成人精品电影| 久久婷婷成人综合色麻豆| 纵有疾风起免费观看全集完整版| 在线观看一区二区三区激情| 十八禁高潮呻吟视频| 丁香欧美五月| 成年动漫av网址| 亚洲精品美女久久av网站| 精品免费久久久久久久清纯 | 国产精品自产拍在线观看55亚洲 | 天天添夜夜摸| 亚洲av国产av综合av卡| 女人高潮潮喷娇喘18禁视频| 国产精品久久久人人做人人爽| 久久久水蜜桃国产精品网| 黄频高清免费视频| 母亲3免费完整高清在线观看| 久久精品亚洲av国产电影网| 久久久久视频综合| 国产精品电影一区二区三区 | 久久国产亚洲av麻豆专区| 色婷婷久久久亚洲欧美| 日韩欧美国产一区二区入口| 欧美在线黄色| 91精品三级在线观看| 高清黄色对白视频在线免费看| 免费观看a级毛片全部| 极品人妻少妇av视频| 色视频在线一区二区三区| 欧美精品啪啪一区二区三区| 久久久久久久大尺度免费视频| 99re在线观看精品视频| 19禁男女啪啪无遮挡网站| 久久亚洲真实| 久久人妻福利社区极品人妻图片| 亚洲人成电影观看| 多毛熟女@视频| 美女福利国产在线| 少妇粗大呻吟视频| 人人妻人人爽人人添夜夜欢视频| 满18在线观看网站| 夜夜夜夜夜久久久久| 国产av国产精品国产| 亚洲熟女毛片儿| 国产日韩欧美在线精品| 大型av网站在线播放| 久久免费观看电影| 91麻豆av在线| 午夜福利影视在线免费观看| 午夜福利视频在线观看免费| 精品亚洲成a人片在线观看| 丁香六月天网| 视频区欧美日本亚洲| 99riav亚洲国产免费| 男女高潮啪啪啪动态图| 9热在线视频观看99| 手机成人av网站| 在线观看www视频免费| 热re99久久精品国产66热6| 午夜福利一区二区在线看| 自线自在国产av| 999久久久国产精品视频| 热99久久久久精品小说推荐| 五月开心婷婷网| 青草久久国产| 国产激情久久老熟女| av电影中文网址| 黑人巨大精品欧美一区二区蜜桃| 亚洲成av片中文字幕在线观看| 51午夜福利影视在线观看| 一进一出抽搐动态| 国产精品一区二区免费欧美| 亚洲精品一卡2卡三卡4卡5卡| 少妇精品久久久久久久| 人人妻人人澡人人爽人人夜夜| 如日韩欧美国产精品一区二区三区| 一级毛片电影观看| 丰满饥渴人妻一区二区三| 国产精品免费一区二区三区在线 | 成年人午夜在线观看视频| 日韩欧美国产一区二区入口| 夜夜夜夜夜久久久久| 操美女的视频在线观看| av国产精品久久久久影院| 精品国产亚洲在线| 久热这里只有精品99| cao死你这个sao货| 精品午夜福利视频在线观看一区 | 新久久久久国产一级毛片| 性高湖久久久久久久久免费观看| 国产精品偷伦视频观看了| 精品欧美一区二区三区在线| 欧美乱码精品一区二区三区| 亚洲国产成人一精品久久久| 国产91精品成人一区二区三区 | 一区二区三区国产精品乱码| 天堂8中文在线网| 欧美精品亚洲一区二区| 啦啦啦免费观看视频1| 中文字幕另类日韩欧美亚洲嫩草| 菩萨蛮人人尽说江南好唐韦庄| 女人爽到高潮嗷嗷叫在线视频| 好男人电影高清在线观看| 亚洲熟妇熟女久久| 18在线观看网站| 免费久久久久久久精品成人欧美视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产精品一区二区三区在线| 两人在一起打扑克的视频| 美国免费a级毛片| 久久久久久亚洲精品国产蜜桃av| 国精品久久久久久国模美| 国产一区二区三区视频了| www.999成人在线观看| 国产黄频视频在线观看| 一本色道久久久久久精品综合| 精品国产国语对白av| 亚洲国产成人一精品久久久| 老汉色av国产亚洲站长工具| 人人妻人人澡人人看| 免费在线观看视频国产中文字幕亚洲| 精品国产一区二区三区四区第35| 亚洲av国产av综合av卡| 丝袜美足系列| 美女午夜性视频免费| 又大又爽又粗| 免费在线观看日本一区| 男女午夜视频在线观看| 满18在线观看网站| 精品第一国产精品| 日韩制服丝袜自拍偷拍| 久久中文字幕一级| 一级片'在线观看视频| 午夜两性在线视频| 亚洲中文av在线| 一级片'在线观看视频| 不卡一级毛片| 久久久国产精品麻豆| e午夜精品久久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 美女午夜性视频免费| 一区二区三区精品91| 一级毛片精品| 丁香欧美五月| 麻豆av在线久日| 啦啦啦在线免费观看视频4| 久久久久久免费高清国产稀缺| 免费一级毛片在线播放高清视频 | 黄色a级毛片大全视频| av网站免费在线观看视频| 亚洲中文字幕日韩| 国产麻豆69| kizo精华| 王馨瑶露胸无遮挡在线观看| 亚洲色图av天堂| 视频区欧美日本亚洲| 色婷婷av一区二区三区视频| 无限看片的www在线观看| 巨乳人妻的诱惑在线观看| 欧美精品人与动牲交sv欧美| 国产精品电影一区二区三区 | 97在线人人人人妻| 久久久久国内视频| 久久人妻av系列| www.自偷自拍.com| 狠狠狠狠99中文字幕| 久久国产精品人妻蜜桃| 久久久久久久大尺度免费视频| 天天添夜夜摸| 美女扒开内裤让男人捅视频| 9热在线视频观看99| 国产无遮挡羞羞视频在线观看| 狠狠婷婷综合久久久久久88av| 啦啦啦中文免费视频观看日本| 久热这里只有精品99| 极品人妻少妇av视频| 久久久精品区二区三区| 99国产精品一区二区蜜桃av | 亚洲av片天天在线观看| 午夜福利一区二区在线看| 99热网站在线观看| 中文字幕制服av| 一级,二级,三级黄色视频| 亚洲av片天天在线观看| 精品少妇黑人巨大在线播放| 老司机福利观看| 国产男靠女视频免费网站| 国产精品成人在线| 一二三四在线观看免费中文在| 午夜福利在线免费观看网站| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久av美女十八| 最黄视频免费看| 国产日韩欧美视频二区| 狠狠婷婷综合久久久久久88av| 亚洲成人手机| 亚洲成av片中文字幕在线观看| 2018国产大陆天天弄谢| 精品人妻1区二区| 麻豆国产av国片精品| a级毛片黄视频| 中文字幕人妻丝袜制服| 精品久久久精品久久久| 亚洲精品国产精品久久久不卡| 久久ye,这里只有精品| 久久狼人影院| 亚洲国产av影院在线观看| 国产精品久久久久久精品古装| 人成视频在线观看免费观看| netflix在线观看网站| 国产在线视频一区二区| 久久久久精品国产欧美久久久| 巨乳人妻的诱惑在线观看| 无人区码免费观看不卡 | 亚洲熟女精品中文字幕| 手机成人av网站| 成在线人永久免费视频| 精品一区二区三区av网在线观看 | 亚洲av欧美aⅴ国产| 亚洲国产欧美网| 欧美激情极品国产一区二区三区| 人人妻人人添人人爽欧美一区卜| 久久精品aⅴ一区二区三区四区| 别揉我奶头~嗯~啊~动态视频| 黄网站色视频无遮挡免费观看| 欧美精品一区二区大全| 丁香六月欧美| 国产精品av久久久久免费| 首页视频小说图片口味搜索| 日韩制服丝袜自拍偷拍| 操出白浆在线播放| 动漫黄色视频在线观看| 丝袜人妻中文字幕| 9热在线视频观看99| 精品福利永久在线观看| 高清欧美精品videossex| 国产三级黄色录像| 精品少妇久久久久久888优播| 欧美黄色片欧美黄色片| 亚洲欧洲精品一区二区精品久久久| 日韩精品免费视频一区二区三区| 成人国语在线视频| 亚洲精品在线观看二区| 国产熟女午夜一区二区三区| 久久久精品区二区三区| 免费人妻精品一区二区三区视频| 欧美成人午夜精品| 久久久久久久久久久久大奶| 亚洲一区中文字幕在线| 狠狠婷婷综合久久久久久88av| 亚洲中文日韩欧美视频| 久久久久精品人妻al黑| 91大片在线观看| 两个人免费观看高清视频| 中文字幕色久视频| 免费少妇av软件| tube8黄色片| 日本vs欧美在线观看视频| 欧美 日韩 精品 国产| 日本av免费视频播放| 国产精品亚洲一级av第二区| 午夜福利在线观看吧| 国产老妇伦熟女老妇高清| 大陆偷拍与自拍| 国产成+人综合+亚洲专区| 性高湖久久久久久久久免费观看| 亚洲 欧美一区二区三区| 亚洲av美国av| 黄片播放在线免费| 考比视频在线观看| 亚洲第一欧美日韩一区二区三区 | 久久久久视频综合| 丁香欧美五月| 亚洲avbb在线观看| 国产不卡一卡二| 大码成人一级视频| 国产一区有黄有色的免费视频| 国产97色在线日韩免费| 老熟女久久久| 久久久精品国产亚洲av高清涩受| 亚洲专区中文字幕在线| 亚洲人成电影免费在线| 久久久久久久大尺度免费视频| 日韩有码中文字幕| 男人舔女人的私密视频| 久久这里只有精品19| 亚洲少妇的诱惑av| 成人国语在线视频| 国产主播在线观看一区二区| 一级毛片电影观看| 国产av一区二区精品久久| 丰满人妻熟妇乱又伦精品不卡| 考比视频在线观看| 人人妻人人澡人人看| 十八禁高潮呻吟视频| 少妇的丰满在线观看| 天天影视国产精品| 美女扒开内裤让男人捅视频| 日韩成人在线观看一区二区三区| 一本—道久久a久久精品蜜桃钙片| 精品亚洲乱码少妇综合久久| 久久久久视频综合| 欧美在线黄色| 午夜久久久在线观看| 亚洲第一青青草原| 伦理电影免费视频| 久久中文看片网| 久久人人爽av亚洲精品天堂| 19禁男女啪啪无遮挡网站| 国产99久久九九免费精品| 老汉色av国产亚洲站长工具| 亚洲欧美日韩高清在线视频 | 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲| 一级片'在线观看视频| 亚洲少妇的诱惑av| 日韩一卡2卡3卡4卡2021年| 啦啦啦免费观看视频1| 午夜激情久久久久久久| 久久亚洲精品不卡| 久久ye,这里只有精品| 99re6热这里在线精品视频| 脱女人内裤的视频| 久久国产精品大桥未久av| 国产亚洲精品一区二区www | av网站免费在线观看视频|