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

    彈丸在炮口的狀態(tài)參數(shù)對地面密集度影響研究

    2020-06-18 03:33:38錢林方陳光宋王明明
    兵工學(xué)報 2020年5期
    關(guān)鍵詞:狀態(tài)參數(shù)炮口密集度

    錢林方, 陳光宋, 王明明

    (1.西北機電工程研究所, 陜西 咸陽 712099; 2.南京理工大學(xué) 機械工程學(xué)院, 江蘇 南京 210094)

    0 引言

    火炮發(fā)射無動力飛行彈丸的地面密集度是衡量中大口徑火炮性能優(yōu)劣的一個重要指標(biāo),影響密集度的主要因素有彈丸在炮口的狀態(tài)參數(shù)、彈丸參數(shù)(物理和幾何)、氣象條件、氣動載荷等隨機誤差。彈丸在炮口的狀態(tài)參數(shù)由火炮發(fā)射過程賦予,理清彈丸在炮口的狀態(tài)參數(shù)誤差對地面密集度影響,為火炮總體設(shè)計提供依據(jù),對提高火炮射擊精度具有重要意義。

    火炮射擊精度問題一直是火炮設(shè)計和工程應(yīng)用研究的熱點和難點。眾多的國內(nèi)外學(xué)者在射擊精度方面發(fā)表了大量研究成果。例如:韓子鵬等[1]通過建立不同復(fù)雜程度的外彈道方程,系統(tǒng)地分析了影響彈丸落點分布的因素,并給出了炮口狀態(tài)參數(shù)對射擊精度影響的物理解釋;郭錫福[2]對遠程火炮射擊精度的計算分析模型、影響因素、試驗技術(shù)等進行了深入研究;王寶元[3]討論了影響火炮射擊密集度的因素,并歸納了火炮武器型號研制和提高射擊密集度所采取的技術(shù)措施;芮筱亭等[4]基于最大熵法對武器系統(tǒng)的射擊密集度進行了分析,結(jié)果顯示,彈丸落點的分布為具有偏態(tài)性質(zhì)的類正態(tài)分布;張領(lǐng)科等[5]計算分析了點火具點火延遲時間、底排裝置工作時間和底排藥劑部分脫落質(zhì)量偏差對射程散布影響;曹寧等[6]利用逐步回歸法研究了彈丸初速、質(zhì)量和偏心距等因素對火炮射擊密集度影響;張海燕等[7]基于拉丁超立方試驗設(shè)計方法和極差分析法進行了射擊密集度靈敏度計算,分析了結(jié)構(gòu)參數(shù)對射擊密集度的影響規(guī)律;王瑞林等[8]研究了高、低射擊頻率對射擊精度的影響;王麗群等[9]提出一種面向指標(biāo)要求的隨機因素參數(shù)區(qū)間計算方法,在給定密集度指標(biāo)的前提下,給出了最終的優(yōu)選參數(shù)方案;李建中等[10]分析了影響彈丸起始擾動的重要參數(shù),重點研究了彈丸質(zhì)量偏心對大口徑火炮射擊精度影響;Qian等[11]采用基于統(tǒng)計信息的區(qū)間方法分析了彈炮耦合系統(tǒng)的不確定性傳播問題。在國際上,國際彈道會議以及相關(guān)的其他會議和學(xué)術(shù)論文等對火炮射擊精度相關(guān)的理論和技術(shù)進行了探討,例如:Dursun[12]探討了炮架剛度、彈丸初速和制造公差等對彈丸散布影響;Rabbath等[13]通過比較樣本標(biāo)準差置信區(qū)間的假設(shè)檢驗方法來評估火炮的射擊精度;Khalil等[14]利用基于工程算法的火箭彈、導(dǎo)彈、炮彈快速設(shè)計平臺PRODAS軟件建立了M107型155 mm彈丸6自由度的發(fā)射動力學(xué)模型,分析了彈丸質(zhì)量、初速、炮口擺角與偏角對射擊精度影響。國內(nèi)外眾多的研究和取得的成果對認識和深入研究火炮射擊精度提供了重要的幫助,也為進一步提高火炮射擊精度奠定了良好的基礎(chǔ)。

    近年來,不確定性問題逐漸成為國內(nèi)外學(xué)者研究的熱點,開展了深入廣泛的研究,發(fā)展了多種不確定性量化和分析方法,例如概率方法[15-16]、非概率方法[17-18]、混合方法[19-21]等,這些方法為解決工程中的具體問題提供了有效途徑。最近,稀疏網(wǎng)格數(shù)值積分(SGNI)方法被成功引入到不確定性傳播分析中,并取得了較好的效果[22-23],它是基于一種特殊的張量積法則,通過低階積分點的嵌入式配置方法,在多變量空間構(gòu)建高階的積分點,形成積分精度較高的SGNI點,有效避免傳統(tǒng)高維積分法因維度較大而引起的“維度災(zāi)難”問題。

    本文基于6自由度剛體外彈道方程,考慮無動力飛行彈丸參數(shù)和炮口狀態(tài)參數(shù)誤差影響,采用SGNI獲得彈丸落點的統(tǒng)計特性參數(shù),在此基礎(chǔ)上利用最大熵估計方法,獲得彈丸落點散布的統(tǒng)計特性。最后,通過與蒙特卡洛方法對比驗證了本文方法的有效性,并利用本文方法計算分析了某155 mm加榴炮的地面密集度。

    1 剛體外彈道方程

    將彈丸視為空間運動的剛體,則彈丸有6個自由度。定義慣性坐標(biāo)系為Oixiyizi,彈體坐標(biāo)系為Opxpypzp,其中Op為彈丸的幾何中心點,Opxp沿彈軸方向,Opypzp平面垂直于Opxp,則彈丸的運動可用Opxpypzp相對Oixiyizi的3個平動位移量u和3個角度位移量ψ描述,相應(yīng)速度量分別為v和ω,加速度分別為a和ε.

    為推導(dǎo)動力學(xué)方程方便,令

    (1)

    彈丸運動關(guān)系可表示為

    (2)

    彈丸運動的虛功率方程可表示為

    (3)

    由外彈道理論可知,作用在彈丸上的合外力f包含了阻力fx、升力fy、馬格努斯力fz、哥氏力fk和重力fg,則合力f為

    f=fx+fy+fz+fk+fg.

    (4)

    作用在彈丸上的合外力矩t包含了靜力矩tz、赤道阻尼力矩tzz、極阻尼力矩txz、馬格努斯力矩ty,則合外力矩t為

    t=tz+tzz+txz+ty.

    (5)

    (2)式代入(3)式,根據(jù)變分原理可得彈丸動力學(xué)方程為

    (6)

    F=DTQ.

    (7)

    (6)式中左端兩項均為慣性力,其中:第1項為主慣性力,與加速度有關(guān);第2項為附加慣性力,包括離心力、哥氏力等。

    2 彈丸落點統(tǒng)計參數(shù)的求解

    2.1 基于SGNI的彈丸落點統(tǒng)計矩求解

    考慮上述參數(shù)影響,無動力飛行的彈丸外彈道方程可表示為

    (8)

    式中:t為彈丸飛行時間。

    隨機初始條件為彈丸在炮口的狀態(tài)參數(shù)

    (9)

    為分析參數(shù)誤差對彈丸落點分布的影響,需要根據(jù)已知隨機參數(shù)ξ的聯(lián)合概率密度函數(shù)fξ(ξ)求解方程(8)式獲得彈丸落點的概率分布。記彈丸落點的縱向和橫向位置分別為x=g(ξ)和z=h(ξ),彈丸落點中心矩即是對地面密集度的衡量。以縱向落點分布為例,落點分布的e階原點矩μe與中心矩ce分別為

    (10)

    (11)

    式中:Ωξ為隨機參數(shù)ξ的空間。由前4階統(tǒng)計矩轉(zhuǎn)化而來的均值μ、標(biāo)準差σ、偏度τ和峰度κ分別為

    (12)

    (13)

    對于上述積分點,相應(yīng)的積分權(quán)重可根據(jù)Smolyak方法[25]計算:

    (14)

    (15)

    ωi(ξ)=(ξ-ξ1)(ξ-ξ2)…(ξ-ξi).

    (16)

    2.2 基于最大熵原理的彈丸落點分布估計

    最大熵原理[26]可以描述為在給定部分信息推斷概率分布時,所有可能的分布中存在一個使信息熵取最大值的分布,該分布是最小偏見的。文獻[27]對常見分布的最大熵估計進行了推導(dǎo)和證明,并給出了一般分布的最大熵估計求解方法??紤]到彈丸落點服從類正態(tài)分布[4]的特性,可利用該方法對彈丸落點分布進行估計。

    以縱向落點x=g(ξ)分布為例,定義彈丸落點參數(shù)x的信息熵S為

    (17)

    式中:fx(x)為彈丸縱向落點的概率密度函數(shù);Ωx為積分空間。

    將落點分布的前r階中心矩(本文中取r=4)作為信息熵約束條件,則可建立如下求解fx(x)的優(yōu)化模型:

    (18)

    (18)式為等式約束優(yōu)化問題,可采用拉格朗日方法進行變換,定義如下形式的拉格朗日函數(shù):

    (19)

    式中:λe(e=0,1,…,r)為拉格朗日乘子。令(19)式對fx(x)求導(dǎo)等于0,則滿足最大熵原理的fx(x)[27]可表示為

    (20)

    將表達式(20)式代入約束方程可得如下非線性方程:

    (21)

    (21)式中待定系數(shù)λe與方程數(shù)相同,理論上存在一組最優(yōu)解,通常將(21)式的非線性方程轉(zhuǎn)化為無約束優(yōu)化問題來求解,定義誤差函數(shù)最小,即

    (22)

    式中:λ={λ0,λ1,…,λr}T.

    為確保上述優(yōu)化問題能夠得到較為準確的解,本文利用遺傳算法和序列二次規(guī)劃的方法進行求解,以此來滿足全局搜索和局部精細求解的目的。先通過遺傳算法求解(22)式,并將計算得到的結(jié)果作為初始值,利用序列二次規(guī)劃方法計算獲得最終的解。計算落點統(tǒng)計矩和概率分布的流程如圖1所示。

    3 算例分析

    以某155 mm加榴炮為研究對象,分析最大射程的彈丸落點散布,最大射程角為51°.

    圖1 計算流程Fig.1 Flowchart of computation

    3.1 有效性驗證

    由于難以獲得彈丸落點分布統(tǒng)計特性的解析解,為了驗證本文方法的有效性,本文將蒙特卡洛方法的計算結(jié)果作為參考解。采用稀疏網(wǎng)格積分的精度水平取k=2,根據(jù)前文的推導(dǎo),n=11個參數(shù)的總積分點數(shù)N=2n2+2n+1=265. 計算得到的彈丸落點縱向和橫向概率密度函數(shù)(PDF)如圖2和圖3所示,落點參數(shù)的統(tǒng)計值如表2所示。從計算結(jié)果可看出,本文方法和蒙特卡洛方法獲得的落點統(tǒng)計分布特性吻合得很好,從而驗證了本文方法的有效性。

    表1 彈丸不確定性參數(shù)

    圖2 落點縱向分布Fig.2 Longitudinal dispersion of projectiles

    3.2 彈丸在炮口的單狀態(tài)參數(shù)影響分析

    本節(jié)考慮7個炮口狀態(tài)參數(shù)單獨作用對彈丸落點密集度影響。參考文獻[3]和型號研制工程經(jīng)驗,將參數(shù)范圍進行適當(dāng)放大,給定炮口狀態(tài)參數(shù)的均值和均方差的范圍如表3所示;彈丸參數(shù)的統(tǒng)計特性與表1所示相同,在具體分析中也可根據(jù)不同型號火炮的特點確定對應(yīng)的參數(shù)范圍。采用稀疏網(wǎng)格積分的精度水平取k=2,根據(jù)前文的推導(dǎo),單變量的積分點數(shù)N=2n+1=7. 炮口狀態(tài)參數(shù)誤差與地面密集度的影響呈近似二次關(guān)系,為此,建立如下所示的炮口單狀態(tài)參數(shù)均方差和落點均方差二次形式的參數(shù)化模型:

    圖3 落點橫向分布Fig.3 Transverse dispersion of projectiles

    表2 落點統(tǒng)計特性

    Tab.2 Statistical parameters of projectiles dispersion

    方法方向均值/m均方差/m偏度峰度計算次數(shù)本文縱向30185148.8-0.05062.7656265方法30182146.8-0.06903.02998000蒙特卡縱向30183147.2-0.01542.988110000洛方法30182147.8-0.02412.975815000本文橫向1214.751.7-0.01012.9960265方法1214.051.80.03772.99108000蒙特卡橫向1214.151.70.03773.013810000洛方法1214.251.80.00283.003415000

    表3 彈丸在炮口的狀態(tài)參數(shù)

    (23)

    (24)

    (25)

    根據(jù)炮口單狀態(tài)參數(shù)影響的數(shù)據(jù),通過擬合得到各個模型系數(shù)如表4和表5所示,其中R2為擬合模型的可決系數(shù)。將(23)式、(24)式分別繪制成曲線得各炮口狀態(tài)參數(shù)均方差對彈丸落點縱向和橫向地面密集度Axi和Azi的影響規(guī)律,如圖4~圖10所示。

    表4 模型(23)式的系數(shù)

    表5 模型(24)式的系數(shù)

    圖4 初速均方差對地面密集度影響Fig.4 Effect of σv0 on ground dispersion

    圖5 高低擺角均方差對地面密集度影響Fig.5 Effect of σφv on ground dispersion

    圖6 方向擺角均方差對地面密集度影響Fig.6 Effect of σφa on ground dispersion

    圖7 高低擺角速度均方差對地面密集度影響Fig.7 Effect of on ground dispersion

    圖8 方向擺角速度均方差對地面密集度影響Fig.8 Effect of on ground dispersion

    圖9 高低偏角均方差對地面密集度影響Fig.9 Effect of σψv on ground dispersion

    圖10 方向偏角均方差對地面密集度影響Fig.10 Effect of σψa on ground dispersion

    3.3 彈丸在炮口的多參數(shù)變量影響分析

    彈丸在炮口的多狀態(tài)參數(shù)影響的綜合模型可表示為

    (26)

    (27)

    式中:σxi0和σzi0為0階綜合模型,表達式為

    (28)

    其含義是僅考慮與表1中給定的彈丸參數(shù)統(tǒng)計特性時彈丸落點的均方差。

    為驗證該模型的準確性,在表3給定的炮口均方差范圍內(nèi)隨機抽樣5組,將(26)式、(27)式計算得到的結(jié)果和稀疏網(wǎng)格法得到的結(jié)果進行對比驗證,結(jié)果如表6所示。從表6中可看出,(26)式和(27)式的計算結(jié)果和稀疏網(wǎng)格的計算結(jié)果吻合較好,從而驗證了(26)式和(27)式的準確性和可行性。

    表6 彈丸在炮口的狀態(tài)參數(shù)均方差樣本

    (26)式、(27)式對火炮射擊精度設(shè)計具有指導(dǎo)作用,當(dāng)給定彈丸最大射程密集度指標(biāo)xi和zi時,(26)式和(27)式轉(zhuǎn)換成由7個變量構(gòu)成的空間超曲面方程fxi()=xi和fzi()=zi,火炮總體密集度設(shè)計的任務(wù)是根據(jù)給定的最大射程密集度要求xi和zi,設(shè)計一組火炮關(guān)鍵參數(shù),確保彈丸在炮口的狀態(tài)參數(shù)誤差i的綜合作用滿足條件fxi(i)≤xi和fzi(i)≤zi.

    以最大射程xi,max=30 km,縱向地面密集度Axi=1/300為例,由密集度計算公式反求彈丸縱向落點的均方差

    (29)

    另外,將表4中的模型系數(shù)帶入(28)式得到僅考慮表1中給定的彈丸參數(shù)分布特性時彈丸縱向落點的均方差

    (30)

    式中:σxi0占(29)式均方差σxi的36%,此時,彈丸在炮口的狀態(tài)參數(shù)分布導(dǎo)致的彈丸縱向落點均方差占64%.

    4 結(jié)論

    本文考慮了彈丸在炮口的狀態(tài)參數(shù)誤差影響,采用SGNI方法獲得了落點散布的統(tǒng)計特性參數(shù),并利用最大熵估計方法獲得了落點散布的概率密度函數(shù);重點分析了彈丸在炮口的狀態(tài)參數(shù)誤差對地面密集度影響規(guī)律,其結(jié)果可為火炮總體設(shè)計提供理論指導(dǎo)。根據(jù)理論分析和數(shù)值仿真結(jié)果,得到如下主要結(jié)論:

    1) 通過與蒙特卡洛方法計算的統(tǒng)計矩和概率密度對比,驗證了本文方法求解彈丸落點統(tǒng)計矩及其概率密度的正確性和有效性。

    2) 通過炮口單狀態(tài)參數(shù)對地面密集度影響分析,除了初速誤差外,炮口其他6個參數(shù)的誤差對地面密集度均有影響,而且其影響近似呈二次規(guī)律。

    3) 通過炮口多狀態(tài)參數(shù)對地面密集度影響分析,建立了彈丸在炮口的狀態(tài)參數(shù)均方差與彈丸落點均方差之間的綜合映射模型,該模型可為火炮射擊密集度的總體設(shè)計提供理論指導(dǎo)。

    4) 若以射程30 km,縱向地面密集度1/300為指標(biāo)要求,則彈丸參數(shù)ξp對縱向地面密集度的貢獻約為36%,彈丸在炮口的狀態(tài)參數(shù)ξ0對縱向地面密集度的貢獻約為64%.

    猜你喜歡
    狀態(tài)參數(shù)炮口密集度
    基于CKF的大型拖拉機狀態(tài)參數(shù)估計研究
    炮口制退器研究現(xiàn)狀及發(fā)展趨勢
    基于HHT算法的分布式電力系統(tǒng)狀態(tài)自動監(jiān)測方法
    堆石料的三維應(yīng)力分數(shù)階本構(gòu)模型
    某大口徑火炮系列殺爆彈地面密集度影響因素回歸分析
    The validity analysis of ground simulation test for non-ablative thermal protection materials
    帶彈序的彈幕武器立靶密集度測試
    累積海冰密集度及其在認識北極海冰快速變化的作用
    炮口結(jié)構(gòu)參數(shù)對炮口性能的影響分析
    炮口制退器效率計算方法分析
    av片东京热男人的天堂| 久久这里只有精品中国| 久久久久久大精品| 国产真实乱freesex| 黄色丝袜av网址大全| www日本在线高清视频| 亚洲 国产 在线| 五月玫瑰六月丁香| 色播亚洲综合网| 国产精品亚洲av一区麻豆| 99热精品在线国产| 亚洲专区中文字幕在线| 欧美一区二区国产精品久久精品| 最近在线观看免费完整版| 久久久久精品国产欧美久久久| 日韩国内少妇激情av| 日本五十路高清| 国产三级黄色录像| x7x7x7水蜜桃| 国产av麻豆久久久久久久| 动漫黄色视频在线观看| 长腿黑丝高跟| 亚洲国产色片| 亚洲欧美日韩高清在线视频| 欧美一区二区精品小视频在线| 日韩欧美三级三区| 久久久国产成人免费| 国产高清三级在线| 色哟哟哟哟哟哟| 在线免费观看的www视频| 久久天躁狠狠躁夜夜2o2o| 色av中文字幕| 国产精品一区二区三区四区免费观看 | 国产视频一区二区在线看| 午夜福利欧美成人| 久久精品91无色码中文字幕| 九色国产91popny在线| 黄色日韩在线| 国产欧美日韩一区二区三| 免费搜索国产男女视频| 久久久久精品国产欧美久久久| 国产精品1区2区在线观看.| 五月伊人婷婷丁香| 久久精品91无色码中文字幕| 久久久久久久久免费视频了| 午夜视频精品福利| 极品教师在线免费播放| 亚洲欧美精品综合久久99| 无限看片的www在线观看| 亚洲成人精品中文字幕电影| 免费搜索国产男女视频| 亚洲人成网站高清观看| 一区二区三区高清视频在线| a级毛片a级免费在线| 女警被强在线播放| 欧美最黄视频在线播放免费| 日韩精品青青久久久久久| 色综合站精品国产| 欧美丝袜亚洲另类 | 狠狠狠狠99中文字幕| 成人鲁丝片一二三区免费| 国产熟女xx| 免费看a级黄色片| 中文字幕熟女人妻在线| 老司机午夜福利在线观看视频| 国产午夜精品论理片| 亚洲真实伦在线观看| 国产一区二区三区在线臀色熟女| 在线看三级毛片| 亚洲午夜理论影院| 精华霜和精华液先用哪个| 精品免费久久久久久久清纯| 日本三级黄在线观看| 变态另类成人亚洲欧美熟女| 国产蜜桃级精品一区二区三区| 啪啪无遮挡十八禁网站| www.熟女人妻精品国产| 每晚都被弄得嗷嗷叫到高潮| 久久久久久国产a免费观看| 黄色片一级片一级黄色片| 天堂av国产一区二区熟女人妻| 国产爱豆传媒在线观看| 精品国产乱码久久久久久男人| e午夜精品久久久久久久| 中文资源天堂在线| 欧美在线黄色| 精品久久久久久久末码| 中文字幕久久专区| 国产男靠女视频免费网站| 国产精品精品国产色婷婷| av片东京热男人的天堂| 国语自产精品视频在线第100页| 久久这里只有精品中国| 国产真实乱freesex| 波多野结衣巨乳人妻| 久久精品人妻少妇| 欧美激情在线99| 欧洲精品卡2卡3卡4卡5卡区| 三级国产精品欧美在线观看 | 久久精品亚洲精品国产色婷小说| 亚洲人成网站在线播放欧美日韩| 88av欧美| 国产成人一区二区三区免费视频网站| 中文资源天堂在线| 日韩欧美一区二区三区在线观看| 久久久国产成人免费| 亚洲中文日韩欧美视频| 婷婷精品国产亚洲av在线| 叶爱在线成人免费视频播放| 精品久久久久久成人av| 午夜福利在线观看免费完整高清在 | 成年版毛片免费区| 国产精品久久久人人做人人爽| 男人的好看免费观看在线视频| 亚洲人成网站高清观看| 村上凉子中文字幕在线| 男女视频在线观看网站免费| 观看免费一级毛片| 欧美日本亚洲视频在线播放| 中国美女看黄片| 夜夜躁狠狠躁天天躁| 日韩欧美 国产精品| 欧美一区二区国产精品久久精品| 法律面前人人平等表现在哪些方面| 小说图片视频综合网站| 99在线人妻在线中文字幕| 黄色片一级片一级黄色片| 成人永久免费在线观看视频| 久久精品国产清高在天天线| 免费在线观看亚洲国产| 岛国在线观看网站| 美女高潮喷水抽搐中文字幕| 亚洲人成电影免费在线| 不卡av一区二区三区| 日韩欧美一区二区三区在线观看| 精品无人区乱码1区二区| 欧美乱色亚洲激情| 国产伦精品一区二区三区视频9 | 亚洲欧美精品综合久久99| 一级黄色大片毛片| 亚洲无线在线观看| 国产精品久久电影中文字幕| 午夜福利高清视频| 国产伦人伦偷精品视频| 一级毛片女人18水好多| av国产免费在线观看| 欧美xxxx黑人xx丫x性爽| 人妻丰满熟妇av一区二区三区| 999久久久国产精品视频| or卡值多少钱| 亚洲avbb在线观看| 国产亚洲精品久久久久久毛片| 国产欧美日韩精品亚洲av| 国产成人一区二区三区免费视频网站| 又紧又爽又黄一区二区| 国产精品爽爽va在线观看网站| 国产精品亚洲av一区麻豆| 婷婷精品国产亚洲av| 岛国视频午夜一区免费看| 国产麻豆成人av免费视频| 国产午夜精品论理片| 九九久久精品国产亚洲av麻豆 | 国产乱人伦免费视频| 欧美黄色片欧美黄色片| 婷婷丁香在线五月| 12—13女人毛片做爰片一| 成在线人永久免费视频| 最新美女视频免费是黄的| 欧美+亚洲+日韩+国产| 99国产精品一区二区蜜桃av| 99久久精品一区二区三区| 色在线成人网| 黄片大片在线免费观看| 少妇熟女aⅴ在线视频| 午夜免费观看网址| 精品久久蜜臀av无| 999精品在线视频| 嫩草影院精品99| 亚洲avbb在线观看| 99久久久亚洲精品蜜臀av| 国产精品精品国产色婷婷| svipshipincom国产片| 国内精品美女久久久久久| 国产精品自产拍在线观看55亚洲| 精品久久久久久成人av| 精品人妻1区二区| 99久久无色码亚洲精品果冻| 99精品在免费线老司机午夜| 精品国产三级普通话版| 99热6这里只有精品| 曰老女人黄片| 90打野战视频偷拍视频| 法律面前人人平等表现在哪些方面| 亚洲成人久久性| 日韩欧美免费精品| 免费在线观看日本一区| 久久这里只有精品中国| bbb黄色大片| 亚洲色图av天堂| 女生性感内裤真人,穿戴方法视频| 久久久久国产一级毛片高清牌| 亚洲中文字幕日韩| 欧美在线一区亚洲| 看黄色毛片网站| 男女做爰动态图高潮gif福利片| 精品一区二区三区视频在线 | 美女午夜性视频免费| 毛片女人毛片| 久久久久国产精品人妻aⅴ院| 亚洲av美国av| 精品日产1卡2卡| 成人av一区二区三区在线看| 丝袜人妻中文字幕| 国产伦一二天堂av在线观看| 99在线人妻在线中文字幕| 亚洲av电影不卡..在线观看| 久久久久亚洲av毛片大全| a级毛片在线看网站| 国产欧美日韩精品一区二区| 久久久成人免费电影| 18禁黄网站禁片午夜丰满| 国产精品一区二区精品视频观看| 午夜亚洲福利在线播放| 国产激情偷乱视频一区二区| 国产私拍福利视频在线观看| 午夜福利成人在线免费观看| 此物有八面人人有两片| 两性夫妻黄色片| 不卡av一区二区三区| 变态另类丝袜制服| 国产综合懂色| 99re在线观看精品视频| 可以在线观看的亚洲视频| 网址你懂的国产日韩在线| 90打野战视频偷拍视频| 国产在线精品亚洲第一网站| 美女午夜性视频免费| 可以在线观看的亚洲视频| 国产精品99久久99久久久不卡| 国产三级黄色录像| 99热这里只有是精品50| 久久这里只有精品中国| 又爽又黄无遮挡网站| 日本黄大片高清| 亚洲国产精品成人综合色| 久久久久久久久中文| 一进一出好大好爽视频| 老司机福利观看| 成人av一区二区三区在线看| 毛片女人毛片| 小说图片视频综合网站| 亚洲第一电影网av| 搡老岳熟女国产| 一级作爱视频免费观看| 麻豆成人av在线观看| 在线观看美女被高潮喷水网站 | 真人一进一出gif抽搐免费| 日本熟妇午夜| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清在线视频| 色吧在线观看| 国产v大片淫在线免费观看| 我的老师免费观看完整版| 国产成人精品无人区| 亚洲中文日韩欧美视频| 国产成人福利小说| 久久久久精品国产欧美久久久| 色在线成人网| 99久久精品国产亚洲精品| 精品国产超薄肉色丝袜足j| 亚洲av熟女| 一级毛片高清免费大全| 亚洲av成人av| 99riav亚洲国产免费| 精品国内亚洲2022精品成人| 99热只有精品国产| 日韩欧美国产在线观看| 99re在线观看精品视频| 久久草成人影院| 午夜视频精品福利| 久9热在线精品视频| 免费av不卡在线播放| 性色av乱码一区二区三区2| 国产激情偷乱视频一区二区| 国产成人福利小说| 亚洲 欧美一区二区三区| 草草在线视频免费看| 国产日本99.免费观看| 国产极品精品免费视频能看的| 久久久久久久精品吃奶| 国产精品99久久久久久久久| 亚洲狠狠婷婷综合久久图片| 99精品在免费线老司机午夜| 哪里可以看免费的av片| 精品一区二区三区av网在线观看| 日韩欧美三级三区| xxx96com| 亚洲激情在线av| 人妻久久中文字幕网| 波多野结衣高清无吗| 欧美乱妇无乱码| 日本成人三级电影网站| 一边摸一边抽搐一进一小说| 婷婷亚洲欧美| 色播亚洲综合网| 亚洲自偷自拍图片 自拍| 成年女人永久免费观看视频| 可以在线观看毛片的网站| 美女黄网站色视频| 久久热在线av| 宅男免费午夜| xxx96com| 日本精品一区二区三区蜜桃| 午夜精品一区二区三区免费看| 午夜免费观看网址| 制服人妻中文乱码| 国产欧美日韩一区二区三| 校园春色视频在线观看| 国产真人三级小视频在线观看| e午夜精品久久久久久久| 男人舔女人的私密视频| 真人一进一出gif抽搐免费| 日韩中文字幕欧美一区二区| 欧美成人一区二区免费高清观看 | 久久午夜亚洲精品久久| 国产一区二区三区视频了| 久久久久久久精品吃奶| 亚洲色图av天堂| 岛国在线免费视频观看| 欧美乱色亚洲激情| 性色av乱码一区二区三区2| 91九色精品人成在线观看| 国产午夜精品久久久久久| 精品一区二区三区av网在线观看| 亚洲熟妇中文字幕五十中出| 久久久精品大字幕| 免费高清视频大片| 一级毛片女人18水好多| 久久香蕉国产精品| 波多野结衣巨乳人妻| 九色成人免费人妻av| 午夜精品一区二区三区免费看| 国产精品免费一区二区三区在线| 国产在线精品亚洲第一网站| 亚洲av成人av| 日日摸夜夜添夜夜添小说| 亚洲精品美女久久av网站| 久久热在线av| 五月玫瑰六月丁香| 亚洲在线观看片| 欧美黄色片欧美黄色片| 怎么达到女性高潮| 亚洲国产欧美人成| 看免费av毛片| 国产三级黄色录像| netflix在线观看网站| 久久久久久大精品| 亚洲一区二区三区不卡视频| 亚洲av电影在线进入| 天堂av国产一区二区熟女人妻| 国产成人福利小说| 色综合婷婷激情| 露出奶头的视频| 看片在线看免费视频| 女生性感内裤真人,穿戴方法视频| 高潮久久久久久久久久久不卡| 亚洲av美国av| 国产激情欧美一区二区| 成年女人永久免费观看视频| 亚洲成av人片在线播放无| 嫩草影院精品99| 韩国av一区二区三区四区| a级毛片a级免费在线| 日韩欧美在线乱码| 国产精品av久久久久免费| 亚洲av成人一区二区三| 午夜亚洲福利在线播放| 成人国产综合亚洲| 在线观看午夜福利视频| 国产午夜福利久久久久久| 午夜精品在线福利| 一区二区三区高清视频在线| 亚洲人成伊人成综合网2020| 日韩欧美精品v在线| 成年女人毛片免费观看观看9| 久久久国产欧美日韩av| 国产v大片淫在线免费观看| 美女午夜性视频免费| 熟女少妇亚洲综合色aaa.| 亚洲国产看品久久| 伦理电影免费视频| 97超视频在线观看视频| 美女被艹到高潮喷水动态| 麻豆国产97在线/欧美| 精品久久久久久久人妻蜜臀av| 久久久国产精品麻豆| 欧美高清成人免费视频www| 欧美性猛交╳xxx乱大交人| 韩国av一区二区三区四区| 在线观看美女被高潮喷水网站 | 校园春色视频在线观看| 久久天堂一区二区三区四区| 日韩大尺度精品在线看网址| 一个人免费在线观看的高清视频| 欧美日韩国产亚洲二区| 日日干狠狠操夜夜爽| 国产乱人伦免费视频| 丝袜人妻中文字幕| 12—13女人毛片做爰片一| 亚洲人与动物交配视频| 亚洲自拍偷在线| 亚洲精品456在线播放app | 亚洲人成网站高清观看| 日韩三级视频一区二区三区| 中文资源天堂在线| 亚洲aⅴ乱码一区二区在线播放| 一个人看的www免费观看视频| 18禁裸乳无遮挡免费网站照片| 日本黄色视频三级网站网址| 欧美日韩福利视频一区二区| 午夜福利18| 成人亚洲精品av一区二区| 午夜a级毛片| 中文在线观看免费www的网站| 黄片小视频在线播放| 国产探花在线观看一区二区| 精品久久久久久久毛片微露脸| 搞女人的毛片| 国产精品 国内视频| 日本成人三级电影网站| 午夜精品久久久久久毛片777| 久久精品人妻少妇| 91麻豆av在线| 91麻豆精品激情在线观看国产| 精品一区二区三区av网在线观看| 精品久久久久久久毛片微露脸| 亚洲自偷自拍图片 自拍| 91久久精品国产一区二区成人 | 国产成人av激情在线播放| 女生性感内裤真人,穿戴方法视频| netflix在线观看网站| 日本a在线网址| 老司机在亚洲福利影院| 一本一本综合久久| 波多野结衣高清无吗| 两个人的视频大全免费| 非洲黑人性xxxx精品又粗又长| 少妇人妻一区二区三区视频| 日本三级黄在线观看| 欧美黄色片欧美黄色片| 人人妻人人看人人澡| 91av网站免费观看| 一级黄色大片毛片| 操出白浆在线播放| 精华霜和精华液先用哪个| 久久久久精品国产欧美久久久| 一二三四社区在线视频社区8| 午夜日韩欧美国产| 两个人视频免费观看高清| 精品不卡国产一区二区三区| 亚洲成人久久爱视频| 无遮挡黄片免费观看| 亚洲avbb在线观看| 免费高清视频大片| 亚洲专区字幕在线| 精品久久久久久久末码| 一进一出抽搐gif免费好疼| 两人在一起打扑克的视频| 国产精品香港三级国产av潘金莲| 淫妇啪啪啪对白视频| 亚洲无线在线观看| 白带黄色成豆腐渣| 精品一区二区三区视频在线 | 51午夜福利影视在线观看| 日本成人三级电影网站| 黄色片一级片一级黄色片| 18禁黄网站禁片免费观看直播| 中文字幕最新亚洲高清| 国产av不卡久久| 久久精品亚洲精品国产色婷小说| 国模一区二区三区四区视频 | 欧美成人免费av一区二区三区| 色吧在线观看| 国产亚洲精品综合一区在线观看| 99久久无色码亚洲精品果冻| 日韩欧美 国产精品| 国产成人欧美在线观看| 久久中文看片网| 天堂av国产一区二区熟女人妻| 日本a在线网址| 日本黄色视频三级网站网址| 97超级碰碰碰精品色视频在线观看| 热99在线观看视频| 99久久精品一区二区三区| 一夜夜www| 日本黄色视频三级网站网址| 桃色一区二区三区在线观看| 1024香蕉在线观看| 91在线观看av| 亚洲精品久久国产高清桃花| 好男人电影高清在线观看| 亚洲欧美精品综合久久99| 欧美黄色片欧美黄色片| 欧美日韩国产亚洲二区| 午夜免费激情av| 婷婷精品国产亚洲av在线| 午夜成年电影在线免费观看| 一二三四社区在线视频社区8| 久久精品亚洲精品国产色婷小说| 高清毛片免费观看视频网站| 亚洲精品一区av在线观看| 国产av一区在线观看免费| 欧美一区二区国产精品久久精品| 老司机在亚洲福利影院| 成人av在线播放网站| 美女午夜性视频免费| 怎么达到女性高潮| 两个人视频免费观看高清| 亚洲一区二区三区色噜噜| 国产精品亚洲av一区麻豆| 国产精品免费一区二区三区在线| 韩国av一区二区三区四区| 在线观看午夜福利视频| 日本五十路高清| 亚洲欧洲精品一区二区精品久久久| 久久久国产精品麻豆| 婷婷精品国产亚洲av| 国产麻豆成人av免费视频| 久久亚洲精品不卡| 丁香六月欧美| cao死你这个sao货| 男人舔女人的私密视频| 色尼玛亚洲综合影院| 久久久久久大精品| 女同久久另类99精品国产91| 桃色一区二区三区在线观看| 欧美黄色片欧美黄色片| 超碰成人久久| 亚洲国产色片| 日本 av在线| 一级毛片精品| 亚洲五月天丁香| 日韩欧美国产一区二区入口| 国产高清videossex| 久久久精品欧美日韩精品| 九九在线视频观看精品| 欧美黄色淫秽网站| 国产精品一区二区三区四区久久| 成人高潮视频无遮挡免费网站| 成年人黄色毛片网站| 精品人妻1区二区| 午夜成年电影在线免费观看| 露出奶头的视频| 亚洲人成网站在线播放欧美日韩| 国产69精品久久久久777片 | 国产午夜精品论理片| 青草久久国产| 天天躁狠狠躁夜夜躁狠狠躁| 看免费av毛片| 99久久无色码亚洲精品果冻| 日韩欧美国产一区二区入口| 中文在线观看免费www的网站| 国产精品女同一区二区软件 | 麻豆国产av国片精品| 色噜噜av男人的天堂激情| 可以在线观看毛片的网站| 亚洲国产欧美一区二区综合| 午夜免费成人在线视频| 久久久久久久久免费视频了| 最近最新中文字幕大全免费视频| 悠悠久久av| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区三区在线臀色熟女| 国产视频一区二区在线看| 成人鲁丝片一二三区免费| 免费高清视频大片| 国产乱人伦免费视频| 亚洲欧美激情综合另类| 亚洲av美国av| 国产精品国产高清国产av| 午夜激情欧美在线| 精品乱码久久久久久99久播| 亚洲最大成人中文| 99久久国产精品久久久| 黄片小视频在线播放| 观看免费一级毛片| 狠狠狠狠99中文字幕| 久久精品人妻少妇| xxx96com| а√天堂www在线а√下载| 99视频精品全部免费 在线 | 亚洲av免费在线观看| 岛国在线观看网站| x7x7x7水蜜桃| 国产高清视频在线观看网站| 岛国在线观看网站| 99热6这里只有精品| 亚洲精品色激情综合| 老熟妇乱子伦视频在线观看| aaaaa片日本免费| 黄片小视频在线播放| 欧美激情久久久久久爽电影| 亚洲人成伊人成综合网2020| 99精品欧美一区二区三区四区| 曰老女人黄片| 夜夜夜夜夜久久久久| 久久久水蜜桃国产精品网| 在线观看66精品国产| 狂野欧美激情性xxxx| 色综合站精品国产| 日本黄色视频三级网站网址| 欧美日韩亚洲国产一区二区在线观看| 中文字幕熟女人妻在线| 人妻夜夜爽99麻豆av| 伦理电影免费视频|