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

    爆炸波問題中偶然不確定度的量化*

    2016-04-25 08:11:12王瑞利
    高壓物理學(xué)報(bào) 2016年6期
    關(guān)鍵詞:概率密度函數(shù)激波標(biāo)準(zhǔn)差

    梁 霄,王瑞利

    (1.北京應(yīng)用物理與計(jì)算數(shù)學(xué)研究所,北京 100094;2.山東科技大學(xué)數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院,山東青島 266590)

    1 引 言

    基于建模與模擬(Modeling and Simulation,M&S)的設(shè)計(jì)可以明顯縮短產(chǎn)品研發(fā)周期、降低其生產(chǎn)成本,在生產(chǎn)和研發(fā)過程中起著重要的作用,尤其是試驗(yàn)成本很高、風(fēng)險(xiǎn)很大和危害人身安全的研究領(lǐng)域。如何提高建模的逼真度及模擬結(jié)果的精度,從而通過物理或者數(shù)學(xué)模型模擬客觀世界及其演變過程,已成為一個(gè)亟待解決的重要問題[1-8]。

    研究者總是希望M&S方法能夠?qū)ρ芯靠腕w進(jìn)行逼真的刻畫和精確的描述,然而,實(shí)際研究過程中,誤差和不確定度的存在是不可避免的[8]。 特別是對(duì)于一些復(fù)雜的非線性系統(tǒng),初值或邊值的微小差異有時(shí)會(huì)引起模擬結(jié)果劇烈的變化,影響模型的預(yù)測(cè)效果。M&S方法的關(guān)鍵任務(wù)之一就在于掌握誤差與不確定度的規(guī)律,通過采用比較完善的M&S技術(shù),使用更加可靠的算法,將不確定度和誤差減小到最低限度。為了達(dá)到這一目的,分析M&S不確定度的來(lái)源及其性質(zhì)是非常重要的。了解不確定度對(duì)計(jì)算結(jié)果的影響,可以確定最終結(jié)果的可信度。

    不確定度量化(UQ)是一門蓬勃發(fā)展的學(xué)科。美國(guó)非常重視UQ技術(shù)的研究,近幾年一直將UQ技術(shù)作為數(shù)值模擬軟件及其應(yīng)用研究的重點(diǎn)。1991年暫停地下核試驗(yàn)后,美國(guó)為了維持武器庫(kù)的水平與能力,保證武器的可靠性,開展了先進(jìn)模擬和計(jì)算計(jì)劃,進(jìn)一步刺激了UQ的發(fā)展。另外,2001年隸屬于美國(guó)能源部的三大國(guó)家核武器安全實(shí)驗(yàn)室(SNL,LANL,LLNL)共同倡導(dǎo)將UQ方法應(yīng)用于核武器庫(kù)存安全管理,給出了認(rèn)知不確定度和偶然不確定度在反應(yīng)堆故障排除和放射性核廢料處理中的應(yīng)用案例[9-11]。本研究以爆炸波問題為例,利用非嵌入多項(xiàng)式混沌方法進(jìn)行偶然不確定度的量化,并將結(jié)果與Monte Carlo(MC)所得結(jié)果比較,以期為非嵌入多項(xiàng)式系統(tǒng)不確定度的量化提供參考。

    2 不確定度的分類和量化方法

    如圖1所示,不確定度一般分為偶然不確定度和認(rèn)知不確定度。偶然不確定度是由于事物固有的隨機(jī)性而引起的,有時(shí)也稱為可變性和隨機(jī)不確定性,是系統(tǒng)固有的特點(diǎn),即使在理論上也不能約減,如火箭助推器系統(tǒng)中質(zhì)量特性和幾何形狀的變異。目前對(duì)偶然不確定度的數(shù)學(xué)描述是采用概率論的方法,給出輸出量的概率描述。增加數(shù)據(jù)的樣本容量有助于更準(zhǔn)確地確定樣本的分布類型,但并不會(huì)減小不確定度。認(rèn)知不確定度是由于缺乏知識(shí)導(dǎo)致的不確定度,有時(shí)也稱為知識(shí)的不確定性、模型形式的不確定性或主觀不確定性。認(rèn)知不確定度分兩類:一類是認(rèn)可不確定度,即所謂“知道的不知道”,來(lái)源于對(duì)物理現(xiàn)象的理解較差,如斷裂力學(xué); 另一類是盲目不確定度,即所謂“不知道的不知道”,例如軟件中的程序代碼出錯(cuò)[8]。

    圖1 不確定度分類Fig.1 Classification of uncertainty

    不確定度量化方法是利用各種數(shù)學(xué)和統(tǒng)計(jì)學(xué)方法,對(duì)隨機(jī)因素或認(rèn)知缺陷盡可能地給出分布或區(qū)間描述,建立M&S中偶然或認(rèn)知不確定性因素與系統(tǒng)響應(yīng)量之間的關(guān)系,進(jìn)而量化影響系統(tǒng)性能的不確定性,給出系統(tǒng)響應(yīng)量性能分析的不確定度度量方法。

    在計(jì)算流體力學(xué)和計(jì)算固體力學(xué)領(lǐng)域中,不確定度量化的概念和方法發(fā)展迅速[8]。許多方法都可以應(yīng)用于偶然不確定度的量化,其中敏感度分析方法、抽樣方法、代理模型方法、隨機(jī)微分方程方法、參數(shù)率定與優(yōu)化方法已經(jīng)應(yīng)用于復(fù)雜工程中。MC方法是出現(xiàn)最早、應(yīng)用最廣泛的一種抽樣方法,并且使用直接,不需要任何假設(shè),可以得出模型在不確定意義下的精確解,通常作為基準(zhǔn)解[12]。但是MC方法的計(jì)算成本過高,還會(huì)遇到維數(shù)災(zāi)難問題,因此人們已經(jīng)開始尋找MC的替代方法,如隨機(jī)譜方法。

    隨機(jī)譜方法是使用譜逼近隨機(jī)微分方程或參量,然后將其分解為獨(dú)立的確定性分量和隨機(jī)分量來(lái)量化不確定性。其數(shù)學(xué)思想是數(shù)學(xué)模型中的每一個(gè)參量利用正交多項(xiàng)式(如Hermite多項(xiàng)式)展開成無(wú)窮級(jí)數(shù)項(xiàng),在實(shí)際應(yīng)用中取有限項(xiàng)[13-15]。根據(jù)求解系數(shù)方法的不同,將多項(xiàng)式混沌分為嵌入多項(xiàng)式混沌(Intrusive Polynomial Chaos,IPC)和非嵌入多項(xiàng)式混沌(Non-Intrusive Polynomial Chaos,NIPC)。IPC無(wú)法利用已有程序,需要對(duì)已有的數(shù)值求解程序進(jìn)行大量修改或者重新編制計(jì)算程序。

    3 非嵌入多項(xiàng)式混沌方法

    NIPC是把已有的數(shù)值求解程序作為一個(gè)黑匣子,在隨機(jī)空間里通過一定的抽樣方法獲得若干個(gè)樣本點(diǎn),將樣本點(diǎn)輸入到確定性程序求解,然后對(duì)確定性輸出結(jié)果進(jìn)行統(tǒng)計(jì)分析,以獲得相關(guān)數(shù)值求解結(jié)果的統(tǒng)計(jì)特征,進(jìn)而評(píng)估輸入?yún)?shù)或計(jì)算條件的不確定性在計(jì)算過程中傳播的影響。NIPC采用已有的數(shù)值模擬程序,不需要對(duì)控制方程進(jìn)行修改,不需要重新編寫程序。

    NIPC不確定度評(píng)估的計(jì)算流程如下。

    (1) 選取參數(shù)λ服從何種分布。

    (2) 利用相應(yīng)標(biāo)準(zhǔn)隨機(jī)變量,確定參數(shù)λ多項(xiàng)式混沌譜展開,如λ(ξ)=λ0+sξ。

    (3) 對(duì)參數(shù)構(gòu)成的隨機(jī)變量ξ進(jìn)行抽樣{θ1,θ2,…,θM+1},將其帶入?yún)?shù)λ多項(xiàng)式混沌譜展開式中,得到參數(shù)λ的抽樣{ω1,ω2,…,ωM+1},對(duì)每一次抽樣,模型參數(shù)值為確定值ωj,j=1,…,M+1。

    (4) 對(duì)每個(gè)抽樣確定后的模型參數(shù)值,利用確定性程序(例如CFD程序)計(jì)算得到相應(yīng)的數(shù)值解u(x,t),如第j次抽樣得到解為[u(x,t)]j。

    (1)式可采用眾多積分公式求出,例如:假定正交多項(xiàng)式歸一,采用Latin等概率抽樣,則

    (2)

    4 偶然不確定度量化在爆炸波問題中的應(yīng)用

    4.1 爆炸波問題的控制方程和迎風(fēng)格式

    一維Lagrange坐標(biāo)系中爆炸波問題的控制方程為[16]

    (3)

    (3)式也可以寫成非守恒形式

    (4)

    4.2 爆炸波初邊值問題的不確定度量化

    與Lax和Sod問題類似,爆炸波問題也是一種激波管問題。根據(jù)炸藥PBX9401的特性,將初始狀態(tài)等效于炸藥爆轟波CJ狀態(tài),建立一維平面模型:在無(wú)限長(zhǎng)管道中,以x=0為界,左側(cè)是炸藥,采用高壓氣體模型描述,右側(cè)是低壓氣體。初始時(shí)刻,兩側(cè)物質(zhì)均處于靜止?fàn)顟B(tài),中間用薄膜隔離。t=0時(shí)刻,薄膜發(fā)生破裂,產(chǎn)生左行稀疏波和右行激波。假設(shè)該爆炸波問題的計(jì)算區(qū)域?yàn)閇-1,1],初值條件如下:在x≤0區(qū)域,炸藥密度ρL服從期望μ為1.4 g/cm3、標(biāo)準(zhǔn)差σ為5.0×10-4g/cm3的正態(tài)分布,速度uL=-2.6 km/s,壓強(qiáng)pL=37.18 GPa,多方指數(shù)γL=3.1;在x>0區(qū)域,氣體密度ρR=15.3 g/cm3,速度uR=0,壓強(qiáng)pR=1.0 GPa,多方指數(shù)γR=1.4。

    該問題的特點(diǎn)是左、右兩側(cè)氣體的多方指數(shù)不同,屬于多介質(zhì)激波管問題,并且有解析解[17]。該問題的解析解表明,強(qiáng)間斷橫跨激波面和激波后有一個(gè)窄的密度峰,可以用于驗(yàn)證流體力學(xué)程序的計(jì)算能力和格式精度。對(duì)于左端炸藥密度,由于加工過程中炸藥顆粒凝結(jié)的不均勻性,其密度存在隨機(jī)性。假設(shè)炸藥密度的概率密度函數(shù)(PDF)為

    (5)

    圖2 炸藥密度分布直方圖Fig.2 Histogram of explosive density distribution

    確定型爆炸波方程通常使用數(shù)學(xué)模型的單一解描述。對(duì)于不確定系統(tǒng),這是不現(xiàn)實(shí)的,實(shí)際上需要執(zhí)行一系列運(yùn)算評(píng)估模型中輸入不確定度對(duì)輸出不確定度的影響,即不確定度傳播。確定型方程利用迎風(fēng)格式,計(jì)算到T=0.10時(shí)刻,得到1 000組狀態(tài)響應(yīng)量構(gòu)成的“馬尾圖”。計(jì)算結(jié)果的信息量巨大,必須對(duì)其量化以提取有用的信息。

    圖3 1 000組抽樣下爆炸波問題的計(jì)算結(jié)果Fig.3 Results of blast wave problem solved with 1 000 times sampling by NIPC method

    4.3 爆炸波問題的不確定度傳播

    采用NIPC和MC方法對(duì)爆炸波初值問題開展不確定度的量化,得到T=0.10時(shí)刻計(jì)算區(qū)域中密度、壓力、速度、能量的期望和標(biāo)準(zhǔn)差如圖4所示,其中E為期望,σ為標(biāo)準(zhǔn)差。由圖4可知,MC方法和NIPC方法吻合較好,說明NIPC應(yīng)用于不確定度量化分析是可行的。

    圖4 T=0.10時(shí)刻爆炸波問題中密度、壓強(qiáng)、速度和能量的期望與標(biāo)準(zhǔn)差Fig.4 Expectations and standard variances of density,pressure,velocity and internal energy in blast wave problem at T=0.10

    4.4 激波位置的概率密度函數(shù)

    為進(jìn)一步驗(yàn)證NIPC方法的有效性,利用NIPC方法研究初邊值輸入不確定度對(duì)爆炸波問題激波位置的影響,計(jì)算得到T=0.10時(shí)刻激波位置的概率密度函數(shù),如圖5所示。可以看出,激波位置近似服從正態(tài)分布,激波位置的期望為0.19,標(biāo)準(zhǔn)差為4.13×10-4。圖6為采用MC方法取樣106次與采用NIPC取樣103次計(jì)算所得激波位置的概率密度函數(shù),通過對(duì)比可以看出,MC方法的概率密度函數(shù)和NIPC方法得到的概率密度函數(shù)基本吻合,說明了NIPC方法的可靠性和有效性。

    圖5 激波位置的概率密度分布函數(shù)Fig.5 Probability density funtion of shock wave position by NIPC at T=0.10

    圖6 激波位置的概率密度分布函數(shù)Fig.6 Comparison of probability density funtion of shock wave position by MC and NIPC at T=0.10

    5 結(jié) 論

    給出建模中輸入初邊值偶然不確定度抽樣方法及量化結(jié)果,并利用MC和NIPC方法計(jì)算得到偶然不確定度對(duì)輸出結(jié)果期望、標(biāo)準(zhǔn)差及激波位置的影響,驗(yàn)證了NIPC方法在偶然不確定度量化方面的有效性。 此方法對(duì)復(fù)雜工程M&S中偶然不確定度量化具有重要的參考價(jià)值。

    [1] 楊振虎.假設(shè)檢驗(yàn)與CFD不確定度量化分析 [J].航空計(jì)算技術(shù),2003,33(2):5-8.

    YANG Z H.Hypothesis tests and analysis of CFD uncertainty quantification [J].Aeronautical Computer Technique,2003,33(2):5-8.

    [2] 王瑞利,江 松.多物理耦合非線性偏微分方程與數(shù)值解不確定度量化數(shù)學(xué)方法 [J].中國(guó)科學(xué)A輯,2015,45(5):723-728.

    WANG R L,JIANG S.Mathematical methods for uncertainty quantification of the nonlinear partial differential equation and numerical solution [J].Science China Series A,2015,45(5):723-728.

    [3] 劉 全,王瑞利,林 忠,等.流體力學(xué)拉氏程序收斂性及數(shù)值計(jì)算不確定度初探 [J].計(jì)算物理,2013,30(3):346-352.

    LIU Q,WANG R L,LIN Z,et al.Asymptotic convergence analysis and quantification of uncertainty in Lagrangian computation [J].Chinese Journal of Computational Physics,2013,30(3):346-352.

    [4] 王瑞利,劉 全,溫萬(wàn)治.非嵌入式多項(xiàng)式混沌法在爆轟產(chǎn)物JWL參數(shù)評(píng)估中的應(yīng)用 [J].爆炸與沖擊,2015,35(1):9-15.

    WANG R L,LIU Q,WEN W Z.Uncertainty quantification for JWL EOS parameters using non-intrusive polynomial chaos [J].Explosion and Shock Waves,2015,35(1):9-15.

    [5] 唐 煒,王玉明.復(fù)雜系統(tǒng)關(guān)鍵特征參數(shù)確定方法 [J].信息與電子工程,2011,9(1):83-86.

    TANG W,WANG Y M.A method of confirming critical characteristic parameters for complex system [J].Information and Electrical Engineering,2011,9(1):83-86.

    [6] 陳堅(jiān)強(qiáng),張益榮.基于Richardson插值法的CFD驗(yàn)證和確認(rèn)方法的研究 [J].空氣動(dòng)力學(xué)學(xué)報(bào),2012,30(2):176-183.

    CHEN J Q,ZHANG Y R.Verification and validation in CFD based on the Richardson extrapolation method [J].Acta Aerodynamical Sinica,2012,30(2):176-183.

    [7] 劉 全,王瑞利,林 忠.非嵌入多項(xiàng)式混沌方法在拉氏計(jì)算中的應(yīng)用 [J].固體力學(xué)學(xué)報(bào),2013,33(增刊):224-233.

    LIU Q,WANG R L,LIN Z.Uncertainty quantification for Lagrangian computational using non-intrusive polynomial chaos [J].Chinese Journal of Solid Mechanics,2013,33(Suppl):224-233.

    [8] 張 楠,沈泓萃,姚惠之.阻力和流場(chǎng)的CFD不確定度分析探討 [J].船舶力學(xué),2008,12(2):211-224.

    ZHANG N,SHEN H C,YAO H Z.Uncertainty analysis in CFD for resistance and flow field [J].Journal of Ship Mechanics,2008,12(2):211-224.

    [9] OBERKAMPF W,ROY C.Verification and validation in the scientific computing [M].New York:Cambridge University Press,2010.

    [10] 馬智博,鄭 淼,尹建偉,等.爆轟模擬不確定度的量化方法 [J].計(jì)算物理,2011,28(1):66-74.

    MA Z B,ZHENG M,YIN J W,et al.Quantification of uncertainties in detonation simulations [J].Chinese Journal of Computational Physics,2011,28(1):66-74.

    [11] HELTON J,JOHNSON J,SALLABERRY C,et al.Survey of sampling-based method for uncertainty and sensitivity analysis [J].Reliab Eng Syst Safe,2006,91(11):1175-1209.

    [12] HELTONl J,JOHNSON J,SALLABERRY C.Quantification of margins and uncertainties:example analyses from reactor safety and radioactive waste disposal involving the separation of aleatory and epistemic uncertainty [J].Reliab Eng Syst Safe,2008,99(1):175-195.

    [13] JANSSEN H.Monte-Carlo based uncertainty analysis:sampling efficiency and sampling convergence [J].Reliab Eng Syst Safe,2013,109:123-132.

    [14] 皮 霆,張?jiān)魄?吳景錸.基于多項(xiàng)式混沌方法的柔性多體系統(tǒng)不確定性分析 [J].中國(guó)機(jī)械工程,2011,22(19):2341-2348.

    PI T,ZHANG Y Q,WU J L.Uncertainty analysis of flexible multibody systems using the polynomial chaos methods [J].China mechanical Engineering,2011,22(19):2341-2348.

    [15] 王曉東,康 順.多項(xiàng)式混沌方法在隨機(jī)方腔流動(dòng)模擬中的應(yīng)用 [J].中國(guó)科學(xué)E輯,2011,41(6):790-798.

    WANG X D,KANG S.Application of polynomial chaos on numerical simulation of stochastic cavity flow [J].Science China Series E,2011,41(6):790-798.

    [16] CHEN Q,GOTTLIEB D,HESTHAVEN J.Uncertainty analysis for the steady state flows in a dual throat nozzle [J].J Comput Phys,2005,204:119-148.

    [17] 水鴻壽.一維流體力學(xué)差分方法 [M].北京:國(guó)防工業(yè)出版社,1998.

    SHUI H S.The difference method of one dimensional fluid mechanics [M].Beijing:National Defense Press,1998.

    猜你喜歡
    概率密度函數(shù)激波標(biāo)準(zhǔn)差
    冪分布的有效估計(jì)*
    用Pro-Kin Line平衡反饋訓(xùn)練儀對(duì)早期帕金森病患者進(jìn)行治療對(duì)其動(dòng)態(tài)平衡功能的影響
    一種基于聚類分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斜激波入射V形鈍前緣溢流口激波干擾研究
    已知f(x)如何求F(x)
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    對(duì)于平均差與標(biāo)準(zhǔn)差的數(shù)學(xué)關(guān)系和應(yīng)用價(jià)值比較研究
    基于概率密度函數(shù)的控制系統(tǒng)性能評(píng)價(jià)
    非高斯隨機(jī)分布系統(tǒng)自適應(yīng)控制算法的研究
    777米奇影视久久| 2022亚洲国产成人精品| 久久久久国产网址| 国产高清三级在线| 天堂网av新在线| 丝瓜视频免费看黄片| 国产精品久久久久久精品古装| 又爽又黄无遮挡网站| 午夜免费观看性视频| 美女脱内裤让男人舔精品视频| 日本免费在线观看一区| 精品一区二区免费观看| 一个人观看的视频www高清免费观看| 下体分泌物呈黄色| 国产精品久久久久久精品古装| 日韩在线高清观看一区二区三区| 日韩欧美 国产精品| 国产乱人视频| 噜噜噜噜噜久久久久久91| 简卡轻食公司| 日韩一区二区视频免费看| 三级男女做爰猛烈吃奶摸视频| 国产精品无大码| 日本猛色少妇xxxxx猛交久久| 高清欧美精品videossex| 国产在线一区二区三区精| 免费大片黄手机在线观看| 免费看a级黄色片| 精品久久久久久久久亚洲| 国内揄拍国产精品人妻在线| 2021少妇久久久久久久久久久| 日韩大片免费观看网站| 少妇丰满av| 亚洲第一区二区三区不卡| a级一级毛片免费在线观看| 极品少妇高潮喷水抽搐| 成年av动漫网址| 99久久精品国产国产毛片| 亚洲最大成人手机在线| 欧美激情在线99| 国产亚洲最大av| 各种免费的搞黄视频| 午夜激情久久久久久久| 久久精品综合一区二区三区| 久久久久国产精品人妻一区二区| 人人妻人人看人人澡| 久久久久久九九精品二区国产| 国产男人的电影天堂91| 国产黄色视频一区二区在线观看| 人人妻人人看人人澡| 亚洲av中文av极速乱| 久久久久久国产a免费观看| 精品一区二区三区视频在线| 黄色欧美视频在线观看| 久久精品人妻少妇| 激情 狠狠 欧美| 水蜜桃什么品种好| 国产毛片在线视频| 最新中文字幕久久久久| 久久久久九九精品影院| 男女无遮挡免费网站观看| 免费观看av网站的网址| 爱豆传媒免费全集在线观看| 午夜激情福利司机影院| 在线精品无人区一区二区三 | 黄片wwwwww| 久久97久久精品| 免费黄网站久久成人精品| 黄色一级大片看看| 精品久久久久久电影网| 久久久色成人| 国产伦精品一区二区三区四那| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区 | 久久这里有精品视频免费| 22中文网久久字幕| 欧美性感艳星| 亚洲精品乱码久久久v下载方式| 国产白丝娇喘喷水9色精品| av国产免费在线观看| 国产黄片视频在线免费观看| 成年av动漫网址| 久久久久久久久久人人人人人人| 日韩制服骚丝袜av| 五月天丁香电影| 99久久精品一区二区三区| 亚洲一区二区三区欧美精品 | h日本视频在线播放| 欧美bdsm另类| 久久国内精品自在自线图片| 亚洲av中文字字幕乱码综合| 免费看a级黄色片| 国产男人的电影天堂91| 欧美成人a在线观看| 大香蕉97超碰在线| 秋霞伦理黄片| 少妇人妻一区二区三区视频| 久久精品久久精品一区二区三区| 观看免费一级毛片| 在线观看av片永久免费下载| 校园人妻丝袜中文字幕| tube8黄色片| 永久网站在线| 精品一区在线观看国产| 极品教师在线视频| 一级毛片黄色毛片免费观看视频| 欧美国产精品一级二级三级 | 午夜爱爱视频在线播放| 麻豆成人av视频| 成人黄色视频免费在线看| 久久精品久久精品一区二区三区| av在线老鸭窝| 在线观看av片永久免费下载| 卡戴珊不雅视频在线播放| 国产久久久一区二区三区| 成人亚洲欧美一区二区av| 18禁裸乳无遮挡动漫免费视频 | 国产真实伦视频高清在线观看| 观看美女的网站| 亚洲不卡免费看| 老司机影院毛片| 又粗又硬又长又爽又黄的视频| 国产人妻一区二区三区在| 看黄色毛片网站| 少妇人妻一区二区三区视频| 99久国产av精品国产电影| 成人一区二区视频在线观看| 午夜日本视频在线| 91精品伊人久久大香线蕉| 极品少妇高潮喷水抽搐| 欧美亚洲 丝袜 人妻 在线| 亚洲,一卡二卡三卡| 久久久久久久久久久丰满| av.在线天堂| 国产一区二区亚洲精品在线观看| 又爽又黄无遮挡网站| 国产精品一区二区在线观看99| 丰满人妻一区二区三区视频av| 五月天丁香电影| 亚洲精品一二三| 97人妻精品一区二区三区麻豆| 欧美成人a在线观看| 国产日韩欧美在线精品| freevideosex欧美| 亚洲高清免费不卡视频| 80岁老熟妇乱子伦牲交| 97人妻精品一区二区三区麻豆| 欧美日韩国产mv在线观看视频 | 午夜福利高清视频| 欧美激情在线99| av在线老鸭窝| 亚洲熟女精品中文字幕| 日韩免费高清中文字幕av| 精品久久久久久久人妻蜜臀av| 国产午夜精品久久久久久一区二区三区| 天天躁日日操中文字幕| 国产高清三级在线| av.在线天堂| 成人毛片a级毛片在线播放| 国产老妇女一区| 国产精品久久久久久久电影| 五月开心婷婷网| 国产亚洲5aaaaa淫片| 免费高清在线观看视频在线观看| 亚洲国产欧美人成| 99九九线精品视频在线观看视频| 一级毛片电影观看| 99九九线精品视频在线观看视频| 色视频在线一区二区三区| 18禁在线播放成人免费| 久久精品久久精品一区二区三区| 男插女下体视频免费在线播放| 欧美日韩亚洲高清精品| 日韩电影二区| 国产精品av视频在线免费观看| 久久国内精品自在自线图片| 日本熟妇午夜| 丰满人妻一区二区三区视频av| 亚洲精品成人av观看孕妇| 久久鲁丝午夜福利片| 99热网站在线观看| 国产亚洲一区二区精品| av在线蜜桃| 深爱激情五月婷婷| 国产日韩欧美亚洲二区| 22中文网久久字幕| 一本久久精品| 国产精品国产av在线观看| eeuss影院久久| 一本久久精品| 内射极品少妇av片p| 亚洲欧美成人综合另类久久久| 男人和女人高潮做爰伦理| 熟女av电影| 久久久久精品久久久久真实原创| 51国产日韩欧美| 大香蕉久久网| 91久久精品国产一区二区三区| 亚洲av不卡在线观看| 国产免费一级a男人的天堂| 深爱激情五月婷婷| 国产欧美日韩精品一区二区| 一区二区av电影网| 建设人人有责人人尽责人人享有的 | 人妻 亚洲 视频| av国产免费在线观看| 在现免费观看毛片| 91精品一卡2卡3卡4卡| 草草在线视频免费看| 丝袜脚勾引网站| 日韩,欧美,国产一区二区三区| 一本色道久久久久久精品综合| 一级黄片播放器| 精品国产一区二区三区久久久樱花 | 国产精品一及| 久热久热在线精品观看| 天堂中文最新版在线下载 | 美女被艹到高潮喷水动态| av国产免费在线观看| 成年人午夜在线观看视频| 精品少妇黑人巨大在线播放| 国产精品三级大全| 国产综合懂色| 少妇人妻 视频| 一级黄片播放器| 99视频精品全部免费 在线| 嫩草影院精品99| 欧美高清性xxxxhd video| 男的添女的下面高潮视频| 啦啦啦中文免费视频观看日本| 欧美一级a爱片免费观看看| 国产在线男女| 91久久精品电影网| 少妇猛男粗大的猛烈进出视频 | av国产精品久久久久影院| 成年女人在线观看亚洲视频 | 欧美97在线视频| 成人黄色视频免费在线看| 国产片特级美女逼逼视频| 日本熟妇午夜| 欧美性感艳星| 久久久欧美国产精品| 亚洲怡红院男人天堂| 国产精品国产三级国产专区5o| 欧美日韩一区二区视频在线观看视频在线 | 少妇人妻精品综合一区二区| 欧美日韩一区二区视频在线观看视频在线 | 成人鲁丝片一二三区免费| 欧美三级亚洲精品| 毛片一级片免费看久久久久| 国产伦理片在线播放av一区| 男女国产视频网站| 亚洲丝袜综合中文字幕| 亚洲一级一片aⅴ在线观看| 热re99久久精品国产66热6| 成年女人看的毛片在线观看| 成人无遮挡网站| 熟女av电影| 亚洲av福利一区| 欧美极品一区二区三区四区| 国产 一区精品| 麻豆成人午夜福利视频| 欧美最新免费一区二区三区| av国产精品久久久久影院| 赤兔流量卡办理| 免费av毛片视频| 欧美3d第一页| 亚洲精品日韩在线中文字幕| 婷婷色综合www| 黄色一级大片看看| 免费观看性生交大片5| 日韩一区二区视频免费看| 一级毛片 在线播放| 国产精品蜜桃在线观看| 亚洲精品视频女| 精品久久久久久久末码| 91狼人影院| 日本猛色少妇xxxxx猛交久久| 又爽又黄a免费视频| 亚洲成人中文字幕在线播放| 麻豆精品久久久久久蜜桃| 永久免费av网站大全| 97热精品久久久久久| 欧美人与善性xxx| 97在线视频观看| 国产 一区 欧美 日韩| 成人黄色视频免费在线看| 在线免费观看不下载黄p国产| 午夜福利在线观看免费完整高清在| 国产av国产精品国产| 另类亚洲欧美激情| 国产欧美另类精品又又久久亚洲欧美| 老女人水多毛片| 婷婷色综合大香蕉| 免费黄色在线免费观看| 身体一侧抽搐| 99久久精品国产国产毛片| 网址你懂的国产日韩在线| 美女xxoo啪啪120秒动态图| 最近最新中文字幕大全电影3| 少妇丰满av| 国产精品久久久久久精品古装| 亚洲一区二区三区欧美精品 | 日韩 亚洲 欧美在线| 一级黄片播放器| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人 | 插阴视频在线观看视频| 精品久久久久久久人妻蜜臀av| 日本黄大片高清| 哪个播放器可以免费观看大片| 我要看日韩黄色一级片| 久久久久久久久久久丰满| 国产精品无大码| 最近2019中文字幕mv第一页| 久久久亚洲精品成人影院| 日韩亚洲欧美综合| 99精国产麻豆久久婷婷| 久久韩国三级中文字幕| 伦精品一区二区三区| 噜噜噜噜噜久久久久久91| 久久精品熟女亚洲av麻豆精品| 美女国产视频在线观看| 亚洲av不卡在线观看| 黄片无遮挡物在线观看| 99热全是精品| 中文字幕av成人在线电影| 男人爽女人下面视频在线观看| 成年女人看的毛片在线观看| 久久久色成人| 边亲边吃奶的免费视频| 亚洲欧美清纯卡通| 亚洲经典国产精华液单| 亚洲av一区综合| 亚洲人成网站在线观看播放| 欧美少妇被猛烈插入视频| 成人鲁丝片一二三区免费| 天堂中文最新版在线下载 | 国产白丝娇喘喷水9色精品| av又黄又爽大尺度在线免费看| 亚洲aⅴ乱码一区二区在线播放| 精品久久久噜噜| 美女主播在线视频| 亚洲色图综合在线观看| 精品久久久久久久人妻蜜臀av| 最新中文字幕久久久久| 在线 av 中文字幕| 在线精品无人区一区二区三 | 国产av不卡久久| 高清毛片免费看| 中文字幕av成人在线电影| 久久97久久精品| 成人午夜精彩视频在线观看| 精品国产三级普通话版| 国产精品国产三级国产专区5o| 国产精品一区二区三区四区免费观看| 亚洲av男天堂| 国精品久久久久久国模美| 久久久久国产精品人妻一区二区| 午夜免费男女啪啪视频观看| 色网站视频免费| 高清日韩中文字幕在线| 97热精品久久久久久| 亚洲精品一区蜜桃| 国产av码专区亚洲av| 我的女老师完整版在线观看| 2021少妇久久久久久久久久久| 国产大屁股一区二区在线视频| av网站免费在线观看视频| 欧美另类一区| 精品人妻熟女av久视频| 亚洲av二区三区四区| 亚洲精品久久午夜乱码| 天天一区二区日本电影三级| 亚洲欧美日韩卡通动漫| 成人高潮视频无遮挡免费网站| 一区二区三区四区激情视频| 女人久久www免费人成看片| 久久韩国三级中文字幕| 一个人看的www免费观看视频| 嫩草影院入口| 老司机影院成人| 国产白丝娇喘喷水9色精品| 美女内射精品一级片tv| 成年版毛片免费区| av卡一久久| 成人亚洲精品一区在线观看 | av女优亚洲男人天堂| 大话2 男鬼变身卡| 午夜爱爱视频在线播放| 亚洲一区二区三区欧美精品 | 亚洲国产精品专区欧美| 久久久久久久久久久丰满| 久久久久久久久大av| 精品人妻一区二区三区麻豆| 在线观看三级黄色| 国产色婷婷99| 国产成人精品婷婷| 99久久精品一区二区三区| 97在线视频观看| 国产精品人妻久久久久久| 美女xxoo啪啪120秒动态图| 在线精品无人区一区二区三 | 国产久久久一区二区三区| 亚洲国产精品专区欧美| av在线老鸭窝| 中文精品一卡2卡3卡4更新| 黄色欧美视频在线观看| 国产精品三级大全| 亚洲欧美成人精品一区二区| 亚洲,一卡二卡三卡| 人人妻人人澡人人爽人人夜夜| 久久午夜福利片| 人体艺术视频欧美日本| 美女cb高潮喷水在线观看| 欧美xxⅹ黑人| 久久人人爽人人片av| 91aial.com中文字幕在线观看| 女的被弄到高潮叫床怎么办| 久久99热这里只有精品18| av黄色大香蕉| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美一级a爱片免费观看看| 日本猛色少妇xxxxx猛交久久| 国产精品国产av在线观看| 青青草视频在线视频观看| 国产成人freesex在线| 日日摸夜夜添夜夜爱| 欧美区成人在线视频| 大香蕉97超碰在线| 天天躁夜夜躁狠狠久久av| 国产精品一二三区在线看| 亚洲熟女精品中文字幕| 高清毛片免费看| 爱豆传媒免费全集在线观看| 看黄色毛片网站| 天美传媒精品一区二区| 色综合色国产| 亚洲精品国产av成人精品| 色网站视频免费| 久久久久网色| 国产黄色免费在线视频| 亚洲成人一二三区av| 搡女人真爽免费视频火全软件| 丝袜喷水一区| 夫妻午夜视频| 天天一区二区日本电影三级| 插阴视频在线观看视频| 午夜福利高清视频| 亚洲成人一二三区av| 草草在线视频免费看| 又黄又爽又刺激的免费视频.| av免费在线看不卡| 亚洲人成网站在线观看播放| 久久97久久精品| 精品国产一区二区三区久久久樱花 | 亚洲欧美日韩卡通动漫| 免费观看在线日韩| 麻豆精品久久久久久蜜桃| 亚洲精品日韩在线中文字幕| 熟女人妻精品中文字幕| 成人免费观看视频高清| 日本黄大片高清| 久久韩国三级中文字幕| 波多野结衣巨乳人妻| 丰满人妻一区二区三区视频av| 亚洲久久久久久中文字幕| 亚洲熟女精品中文字幕| 国产精品一二三区在线看| 亚洲精品视频女| 人妻少妇偷人精品九色| 老司机影院毛片| 69av精品久久久久久| 亚洲,欧美,日韩| 国产亚洲最大av| 80岁老熟妇乱子伦牲交| 2018国产大陆天天弄谢| 亚洲人成网站高清观看| 99久久中文字幕三级久久日本| 直男gayav资源| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频 | 少妇人妻精品综合一区二区| 久久精品人妻少妇| 国产一级毛片在线| 免费不卡的大黄色大毛片视频在线观看| 国产伦精品一区二区三区四那| 嫩草影院新地址| 欧美日韩一区二区视频在线观看视频在线 | 欧美97在线视频| 美女xxoo啪啪120秒动态图| .国产精品久久| 精品人妻熟女av久视频| 少妇丰满av| 亚洲三级黄色毛片| 最近中文字幕高清免费大全6| 高清午夜精品一区二区三区| 99热这里只有是精品在线观看| 九九久久精品国产亚洲av麻豆| 久久这里有精品视频免费| 91在线精品国自产拍蜜月| 成年女人看的毛片在线观看| 成人鲁丝片一二三区免费| 七月丁香在线播放| 国产精品国产av在线观看| 少妇猛男粗大的猛烈进出视频 | 香蕉精品网在线| 一级毛片 在线播放| 18+在线观看网站| 91久久精品国产一区二区成人| 在线看a的网站| 天堂俺去俺来也www色官网| 草草在线视频免费看| 日韩欧美一区视频在线观看 | 欧美性感艳星| 在线天堂最新版资源| 在线免费观看不下载黄p国产| 亚洲av成人精品一二三区| 国产欧美日韩精品一区二区| 卡戴珊不雅视频在线播放| 91久久精品电影网| 成年av动漫网址| 国产亚洲午夜精品一区二区久久 | 中文乱码字字幕精品一区二区三区| 深夜a级毛片| 国产成人a区在线观看| 听说在线观看完整版免费高清| 欧美日韩综合久久久久久| 秋霞在线观看毛片| 久久99热这里只有精品18| 国产午夜精品一二区理论片| 男男h啪啪无遮挡| 亚洲av欧美aⅴ国产| 女的被弄到高潮叫床怎么办| 精品久久久久久久久av| 国产白丝娇喘喷水9色精品| 国产美女午夜福利| 久久99热这里只有精品18| 国产精品国产av在线观看| 王馨瑶露胸无遮挡在线观看| 男女啪啪激烈高潮av片| 欧美成人午夜免费资源| 久久精品综合一区二区三区| 老司机影院毛片| 日韩一区二区三区影片| 国产精品国产三级国产av玫瑰| 国产成人午夜福利电影在线观看| 最后的刺客免费高清国语| 2021天堂中文幕一二区在线观| 男女那种视频在线观看| 久久精品久久精品一区二区三区| 交换朋友夫妻互换小说| 国产久久久一区二区三区| 综合色av麻豆| 一级毛片黄色毛片免费观看视频| 男人狂女人下面高潮的视频| 国产成人精品婷婷| 亚洲色图av天堂| 久久人人爽人人片av| 最近最新中文字幕免费大全7| 欧美老熟妇乱子伦牲交| 亚洲最大成人手机在线| 少妇被粗大猛烈的视频| 欧美变态另类bdsm刘玥| 精品久久久精品久久久| 久久国产乱子免费精品| 亚洲三级黄色毛片| 欧美成人a在线观看| 国产精品国产三级专区第一集| 国产乱人偷精品视频| 激情 狠狠 欧美| 精品熟女少妇av免费看| 在线亚洲精品国产二区图片欧美 | 九草在线视频观看| 狂野欧美白嫩少妇大欣赏| 高清日韩中文字幕在线| 日韩国内少妇激情av| 精品人妻熟女av久视频| 日日撸夜夜添| 日本熟妇午夜| 久久久久久伊人网av| 伦精品一区二区三区| 亚洲国产精品成人久久小说| 日韩制服骚丝袜av| 亚洲,欧美,日韩| 亚洲欧美中文字幕日韩二区| 精品一区二区免费观看| 在线精品无人区一区二区三 | 少妇的逼水好多| 最近手机中文字幕大全| 日韩中字成人| 免费av观看视频| 在现免费观看毛片| 国产伦精品一区二区三区视频9| 亚洲国产精品999| 久久久久网色| www.色视频.com| 日韩国内少妇激情av| 久久精品夜色国产| 黑人高潮一二区| 免费大片黄手机在线观看| 国产中年淑女户外野战色| 国产精品人妻久久久影院| 久久久久久久久久久免费av| 日本黄色片子视频| 日本猛色少妇xxxxx猛交久久| 大陆偷拍与自拍| 欧美成人精品欧美一级黄| 亚洲精品国产色婷婷电影| 亚洲欧美一区二区三区黑人 | 99热这里只有精品一区| 听说在线观看完整版免费高清| 国产免费一级a男人的天堂| av黄色大香蕉| 国精品久久久久久国模美|