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

    基于蒙特卡洛的非包容轉(zhuǎn)子爆破安全性評估

    2012-06-05 08:09:08王京婭孫有朝曾海軍徐慶宏張燕軍
    關(guān)鍵詞:災(zāi)難性概率安全性

    王京婭 孫有朝 陸 中 曾海軍 徐慶宏 張燕軍

    (1. 南京航空航天大學(xué),江蘇 南京 210016;2. 上海飛機(jī)設(shè)計(jì)研究院,上海 200232)

    非包容轉(zhuǎn)子爆破是運(yùn)輸類飛機(jī)典型的特殊風(fēng)險(xiǎn),對飛機(jī)的安全性有重要影響,雖然其發(fā)生概率很小,但一旦發(fā)生往往會造成巨大損失,嚴(yán)重威脅飛行安全。為了將非包容轉(zhuǎn)子爆破對飛機(jī)造成的危害減至最小,適航標(biāo)準(zhǔn)對相關(guān)設(shè)計(jì)提出了明確的技術(shù)要求[1-4],并在相應(yīng)咨詢通告中并給出了較為完善的安全性分析與評估方法。

    目前,非包容轉(zhuǎn)子爆破安全性評估主要依據(jù)美國聯(lián)邦航空局咨詢通告AC 20–128ADesign Considerations for Minimizing Hazards Caused by Uncontained Turbine Engine and Auxiliary Power Unit Rotor Failure,其中包括降低非包容轉(zhuǎn)子爆破危害的設(shè)計(jì)措施、非包容轉(zhuǎn)子爆破安全性分析方法以及可接受的剩余風(fēng)險(xiǎn)值[5]。根據(jù)AC 20–128A中給出的方法,進(jìn)行非包容轉(zhuǎn)子爆破安全性分析評估可操作性不強(qiáng),且工作量非常大,借助計(jì)算機(jī)分析評估,可大大提高其分析和評估效率,節(jié)約分析成本,縮短分析周期,是非包容轉(zhuǎn)子爆破安全性分析評估的必然途徑。本文基于蒙特卡洛思想給出非包容轉(zhuǎn)子爆破安全性仿真分析的方法。

    1 問題描述

    非包容轉(zhuǎn)子爆破安全性評估主要通過識別和評估飛機(jī)設(shè)計(jì)針對非包容事件存在的潛在危險(xiǎn),尤其足以導(dǎo)致機(jī)毀人亡的災(zāi)難性危險(xiǎn),計(jì)算非包容轉(zhuǎn)子爆破導(dǎo)致飛機(jī)出現(xiàn)災(zāi)難性危險(xiǎn)的概率。非包容轉(zhuǎn)子爆破安全性定量評估的目標(biāo)為:對采取了非包容設(shè)計(jì)準(zhǔn)則和危險(xiǎn)控制方法的飛機(jī)系統(tǒng),將特定非包容爆破模式導(dǎo)致飛機(jī)出現(xiàn)災(zāi)難性危險(xiǎn)的概率控制在允許范圍內(nèi)。非包容轉(zhuǎn)子爆破安全性評估是保證非包容轉(zhuǎn)子爆破安全性要求,并對其進(jìn)行驗(yàn)證的有效分析方法。

    根據(jù)AC 20–128A中給出的理論方法,本文基于蒙特卡洛方法對非包容轉(zhuǎn)子爆破進(jìn)行仿真分析。蒙特卡洛方法,或稱計(jì)算機(jī)隨機(jī)模擬方法,是一種基于“隨機(jī)數(shù)”的計(jì)算方法。其基本原理為:當(dāng)所要求的問題是某種事件的概率,或者是某個(gè)隨機(jī)變量的預(yù)期值時(shí),可以通過某種“試驗(yàn)”方法,得到這種事件的出現(xiàn)頻率,或者這個(gè)隨機(jī)變數(shù)的平均值,并用它們作為問題的解[6]。本文中,需要求出非包容轉(zhuǎn)子爆破導(dǎo)致飛機(jī)出現(xiàn)災(zāi)難性故障狀態(tài)的危險(xiǎn)概率,即轉(zhuǎn)子碎片飛出機(jī)匣打到飛機(jī)機(jī)身及其部件這一事件會導(dǎo)致飛機(jī)發(fā)生災(zāi)難性事故的概率,根據(jù)特定非包容爆破模式下轉(zhuǎn)子碎片的釋放規(guī)律,對轉(zhuǎn)子碎片飛出這一非包容事件進(jìn)行大量試驗(yàn),對每一次試驗(yàn)的結(jié)果進(jìn)行分析統(tǒng)計(jì),求出飛機(jī)災(zāi)難性事故發(fā)生的概率,當(dāng)試驗(yàn)迭代次數(shù)足夠大時(shí),可以認(rèn)為所有試驗(yàn)結(jié)果的期望值即為飛機(jī)出現(xiàn)災(zāi)難性故障狀態(tài)的危險(xiǎn)概率。

    2 假設(shè)與變量定義

    2.1 分析假設(shè)

    由于非包容問題的復(fù)雜性,非包容轉(zhuǎn)子爆破安全性設(shè)計(jì)、分析和評估均建立在諸多假設(shè)基礎(chǔ)之上。本文參考相關(guān)文獻(xiàn)總結(jié)出以下若干假設(shè)[5]。

    飛機(jī)發(fā)生發(fā)動機(jī)/APU非包容轉(zhuǎn)子爆破的概率為1.0。

    每臺發(fā)動機(jī)發(fā)生非包容轉(zhuǎn)子爆破的概率都是相同的。

    以每臺發(fā)動機(jī)的每一級轉(zhuǎn)子為單位進(jìn)行分析,每一級轉(zhuǎn)子發(fā)生爆破的概率都是相同的。

    不會有兩級、兩級以上轉(zhuǎn)子或一臺以上發(fā)動機(jī)同時(shí)獨(dú)立發(fā)生非包容轉(zhuǎn)子爆破。

    在最大飛散角范圍內(nèi),轉(zhuǎn)子碎片在所有方向上的飛出概率服從均勻分布。

    非包容碎片具有無限的能量,可切斷其掃掠路徑上的所有管路、電纜、鋼索和未加保護(hù)的結(jié)構(gòu)或構(gòu)件,并且不會從其原始軌跡變向,但防護(hù)盾板和其它發(fā)動機(jī)具有抵擋住最大能量碎片的能力。

    2.2 蒙特卡洛仿真變量定義

    用蒙特卡洛方法求解問題時(shí),首先要確定統(tǒng)計(jì)獨(dú)立的隨機(jī)變量。針對非包容轉(zhuǎn)子爆破,單個(gè)非包容轉(zhuǎn)子碎片的任意一條掃掠軌跡,可用在機(jī)身軸向剖面的碎片軌跡飛散角x和在機(jī)身徑向剖面的碎片軌跡平動角y來唯一確定,如圖1所示。

    圖1 碎片軌跡角

    根據(jù)假設(shè)“在最大飛散角范圍內(nèi),轉(zhuǎn)子碎片在所有方向上飛出的概率是均勻分布的”,定義兩個(gè)隨機(jī)變量X和Y,分別代表飛散角變量和平動角變量,兩者均為連續(xù)性隨機(jī)變量。

    2.2.1 軌跡飛散角變量X

    單級轉(zhuǎn)子上的特定類型非包容碎片在轉(zhuǎn)子軸向影響區(qū)域內(nèi)的起始飛散角為s1k、終止飛散角為s2k,軌跡飛散角變量X在區(qū)間[s1k,s2k]上服從均勻分布,X~U(s1k,s2k),X的概率密度函數(shù)

    式中:

    k——碎片類型。

    2.2.2 軌跡平動角變量Y

    碎片軌跡在發(fā)動機(jī)徑向的角度是隨機(jī)的,軌跡平動角變量Y在區(qū)間[0,360]上服從均勻分布,Y~U(0,360),Y的概率密度函數(shù)為

    3 基于蒙特卡洛仿真的安全性評估

    3.1 非包容轉(zhuǎn)子碎片生產(chǎn)與釋放

    轉(zhuǎn)子爆破一般考慮發(fā)動機(jī)的風(fēng)扇葉片、低壓壓氣機(jī)(LPC)、高壓壓氣機(jī)(HPC)、高壓渦輪(HPT)以及低壓渦輪(LPT)的轉(zhuǎn)子發(fā)生爆破的情況。

    轉(zhuǎn)子爆破時(shí),轉(zhuǎn)子碎片從機(jī)匣中沿一定軌跡飛出,基于蒙特卡洛仿真的原理對這一事件進(jìn)行模擬仿真。根據(jù)公式(1)和(2),按兩個(gè)隨機(jī)變量X和Y的分布情況對其進(jìn)行大量的隨機(jī)抽樣,用來模擬特定類型非包容碎片爆破的情況,一次抽樣結(jié)果代表一次非包容事件,用X和Y兩個(gè)隨機(jī)變量來描述每一次抽樣試驗(yàn)情況。

    通過碰撞檢測確定在此軌跡下碎片所能打到的所有飛機(jī)部件、設(shè)備,作為一次碰撞記錄,進(jìn)行N次隨機(jī)抽樣,得到特定碎片類型的N條碰撞記錄。

    3.2 危險(xiǎn)觸發(fā)分析

    用計(jì)算機(jī)進(jìn)行非包容轉(zhuǎn)子爆破安全性評估時(shí),不僅需要判斷每次非包容事件中碎片飛出后會打到哪些部件,還要確定這些被打中的部件會對飛機(jī)持續(xù)安全飛行與著陸造成的后果。

    本文采用被打中部件與最小割集對比的方法對每一次仿真的非包容事件進(jìn)行分析。根據(jù)飛機(jī)功能危險(xiǎn)分析與故障樹分析的結(jié)果,得到所有災(zāi)難性頂事件的對應(yīng)最小割集,最小割集是故障樹分析中能夠引起頂事件發(fā)生的最低數(shù)量的底事件的組合,這些底事件必須全部發(fā)生才會導(dǎo)致頂事件發(fā)生[7]。

    危險(xiǎn)觸發(fā)分析通過分析每次非包容事件是否會觸發(fā)災(zāi)難性最小割集,從而判斷是否會導(dǎo)致災(zāi)難性頂事件的發(fā)生。將最小割集中的各故障模式轉(zhuǎn)換成飛機(jī)中相應(yīng)的具體部件、設(shè)備,進(jìn)而與每次非包容事件的碰撞記錄作對比,如果某個(gè)最小割集中涉及的所有部件完全包含于該碰撞記錄,則說明該次非包容事件會觸發(fā)這個(gè)最小割集。

    通過危險(xiǎn)觸發(fā)分析,確定每次非包容事件是否會觸發(fā)最小割集,會觸發(fā)哪些最小割集,進(jìn)而確定會導(dǎo)致哪些災(zāi)難性頂事件的發(fā)生,該分析結(jié)果作為下一步分析的輸入信息。

    3.3 蒙特卡洛仿真結(jié)果處理方法

    對于每種非包容爆破模式,需要求出其災(zāi)難性危險(xiǎn)概率,判斷是否在允許范圍內(nèi)。本節(jié)對蒙特卡洛仿真分析的結(jié)果進(jìn)行處理,考慮不同風(fēng)險(xiǎn)因子對災(zāi)難性危險(xiǎn)概率的影響,針對某種非包容爆破模式,給出單級轉(zhuǎn)子及整機(jī)的災(zāi)難性危險(xiǎn)概率計(jì)算模型。

    3.3.1 風(fēng)險(xiǎn)因子

    進(jìn)行非包容轉(zhuǎn)子爆破安全性評估時(shí),必需考慮各個(gè)飛行階段的具體情況,因?yàn)槟承╋w機(jī)功能的喪失或系統(tǒng)的損壞,只在特定的飛行階段才會導(dǎo)致災(zāi)難性的后果,同時(shí)在不同飛行階段內(nèi)發(fā)生非包容轉(zhuǎn)子爆破的概率是不同的。飛機(jī)不同的飛行階段由代碼Fj表示,代表第j個(gè)飛行階段。非包容轉(zhuǎn)子爆破必須考慮以下兩種基本風(fēng)險(xiǎn)因子。

    3.3.1.1 轉(zhuǎn)子爆破風(fēng)險(xiǎn)因子

    非包容轉(zhuǎn)子爆破在不同飛行階段Fj的發(fā)生比例(概率)稱為轉(zhuǎn)子爆破風(fēng)險(xiǎn)因子,用Rj來表示。起飛過程中最可能發(fā)生非包容性轉(zhuǎn)子損壞,此時(shí)發(fā)動機(jī)轉(zhuǎn)子承受最大應(yīng)力。表1為發(fā)動機(jī)在各飛行階段發(fā)生非包容轉(zhuǎn)子爆破的比例[5]。

    表1 轉(zhuǎn)子爆破風(fēng)險(xiǎn)因子

    3.3.1.2 頂事件風(fēng)險(xiǎn)因子

    災(zāi)難性頂事件Ti發(fā)生,在第j個(gè)飛行階段會導(dǎo)致飛機(jī)災(zāi)難性事故的概率,稱為頂事件風(fēng)險(xiǎn)因子,用Dij來表示。一個(gè)災(zāi)難性頂事件對應(yīng)一組頂事件風(fēng)險(xiǎn)因子。

    頂事件風(fēng)險(xiǎn)因子根據(jù)具體災(zāi)難性頂事件在特定飛行階段的情況,考慮對飛機(jī)災(zāi)難性事故發(fā)生的影響因素,利用保守的工程經(jīng)驗(yàn)進(jìn)行評估。以“發(fā)動機(jī)推力完全喪失”為例,在不同飛行階段內(nèi),該災(zāi)難性頂事件被觸發(fā)會導(dǎo)致飛機(jī)發(fā)生災(zāi)難性事故的概率,可由該危險(xiǎn)的頂事件風(fēng)險(xiǎn)因子,即“推力完全喪失風(fēng)險(xiǎn)因子”(表2)描述[5]。

    表2 推力完全喪失風(fēng)險(xiǎn)因子

    3.3.2 單次非包容事件災(zāi)難性概率計(jì)算模型

    針對某種非包容爆破模式,以每一次非包容事件為單位進(jìn)行分析,采用第3.2節(jié)中危險(xiǎn)觸發(fā)的分析結(jié)果,根據(jù)被觸發(fā)的災(zāi)難性頂事件的具體情況,計(jì)算該次非包容事件會導(dǎo)致飛機(jī)災(zāi)難性事件的概率。安全性評估過程總共進(jìn)行N次非包容事件仿真,第n次非包容事件導(dǎo)致飛機(jī)災(zāi)難性事故的概率記為pn。下面分3種典型情況進(jìn)行考慮。

    3.3.2.1 非包容事件不會觸發(fā)任何頂事件

    這種情況下,非包容事件不會導(dǎo)致飛機(jī)發(fā)生災(zāi)難性事故,災(zāi)難性概率為:3.3.2.2 非包容事件觸發(fā)1個(gè)頂事件

    觸發(fā)的頂事件記為T1,頂事件風(fēng)險(xiǎn)因子為D1j,如表3所示。

    表3 頂事件T1風(fēng)險(xiǎn)因子

    針對該次非包容事件,導(dǎo)致飛機(jī)災(zāi)難性事故的概率為:

    3.3.2.3 非包容事件觸發(fā)2個(gè)頂事件

    一次非包容事件觸發(fā)了2個(gè)災(zāi)難性頂事件,求這2個(gè)災(zāi)難性頂事件已經(jīng)發(fā)生的情況下,飛機(jī)發(fā)生災(zāi)難性事故的概率。

    觸發(fā)的頂事件記為T1、T2,T1頂事件風(fēng)險(xiǎn)因子為D1j,T2頂事件風(fēng)險(xiǎn)因子為D2j,如表4所示。

    表4 頂事件T1、T2風(fēng)險(xiǎn)因子

    首先分析第j個(gè)飛行階段,飛行階段代碼為Fj,該次非包容事件在此飛行階段中導(dǎo)致飛機(jī)發(fā)生災(zāi)難性事故的概率,可由Fj階段的轉(zhuǎn)子爆破風(fēng)險(xiǎn)因子乘以Fj階段頂事件T1、T2同時(shí)發(fā)生會導(dǎo)致飛機(jī)災(zāi)難性事故的概率[8],即

    根據(jù)公式(5)推導(dǎo)結(jié)果,考慮所有飛行階段,針對該次非包容事件,飛機(jī)發(fā)生災(zāi)難性事故的概率為

    綜合考慮上述3種典型情況,可以將第2種情況看作是第3種情況的特殊情況,綜合公式(3)、(4)、(6)將其推廣到一般情況,第n次非包容事件會導(dǎo)致飛機(jī)發(fā)生災(zāi)難性事故的概率pn為

    式中:

    d——該次非包容事件導(dǎo)致的災(zāi)難性頂事件發(fā)生的個(gè)數(shù)。

    3.3.3 轉(zhuǎn)子爆破災(zāi)難性危險(xiǎn)定量評估模型

    3.3.3.1 單級轉(zhuǎn)子災(zāi)難性危險(xiǎn)概率

    根據(jù)蒙特卡洛仿真原理,當(dāng)?shù)螖?shù)N足夠大時(shí),每次仿真結(jié)果的期望值可以代表需求值。針對非包容轉(zhuǎn)子爆破,單級轉(zhuǎn)子某種碎片爆破模式下,飛機(jī)發(fā)生災(zāi)難性事故的概率P可以用每次非包容事件導(dǎo)致的災(zāi)難性故障概率的期望值來表示,即

    式中:

    N——仿真迭代次數(shù);

    pn——第n次非包容事件會導(dǎo)致飛機(jī)發(fā)生災(zāi)難性事故的概率。

    3.3.3.2 整機(jī)災(zāi)難性危險(xiǎn)概率

    整機(jī)災(zāi)難性危險(xiǎn)概率即非包容轉(zhuǎn)子爆破導(dǎo)致飛機(jī)出現(xiàn)災(zāi)難性故障狀態(tài)的概率。依據(jù)假設(shè),取所有發(fā)動機(jī)每一級轉(zhuǎn)子發(fā)生非包容轉(zhuǎn)子爆破而觸發(fā)災(zāi)難性危險(xiǎn)的概率的均值,作為整機(jī)災(zāi)難性危險(xiǎn)平均概率。

    基于上述分析,以每臺發(fā)動機(jī)為單位進(jìn)行計(jì)算,先求出每臺發(fā)動機(jī)上所有轉(zhuǎn)子危險(xiǎn)概率的均值,然后求出所有發(fā)動機(jī)危險(xiǎn)概率的均值,該危險(xiǎn)平均概率即整機(jī)災(zāi)難性危險(xiǎn)平均概率:

    式中:

    e——發(fā)動機(jī)編號;

    r——轉(zhuǎn)子級號;

    E——發(fā)動機(jī)個(gè)數(shù);

    R——轉(zhuǎn)子級數(shù);

    P——單級轉(zhuǎn)子災(zāi)難性危險(xiǎn)概率;Pz——整機(jī)災(zāi)難性危險(xiǎn)平均概率。

    4 結(jié)論

    本文提出了基于蒙特卡洛仿真原理的非包容轉(zhuǎn)子爆破安全性評估方法,給出了仿真原理與方法的實(shí)施步驟。

    本文給出了轉(zhuǎn)子爆破災(zāi)難性危險(xiǎn)定量評估模型。

    [1]FAR-25Part25- Airworthiness Standards: Transport Category Airplanes[S]. 2003.

    [2]FAR-23Part23-Airworthiness Standards: Normal,Utility, Acrobatic, and Commuter Category Airplanes[S].

    [3]CCAR-25-R4 運(yùn)輸類飛機(jī)適航標(biāo)準(zhǔn)[S].

    [4]CCAR-23-R3 正常類、實(shí)用類、特技類和通勤類飛機(jī)適航規(guī)定[S].

    [5]Advisory Circular AC 20-128ADesign Considerations for Minimizing Hazards Caused by Uncontained Turbine Engine and Auxiliary Power Unit Rotor Failure[R]. America: Federal Aviation Administration, 1997

    [6]陳舜麟. 計(jì)算材料科學(xué)[M]. 北京: 教材出版中心, 2005.

    [7]SAE-ARP4761Guidelines and methods for the safety assessment process on airborne systems and equipments[S]. America: The Engineering Society For Advancing Mobility Land Sea Air and Space, 1996.

    [8]魏宗舒等. 概率論與數(shù)理統(tǒng)計(jì)教程[M]. 北京: 高等教育出版社, 1983.

    猜你喜歡
    災(zāi)難性概率安全性
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    新染料可提高電動汽車安全性
    科學(xué)家說地球變暖速度比以前認(rèn)為的要快,而避免災(zāi)難性后果的窗口期即將結(jié)束
    英語文摘(2021年10期)2021-11-22 08:02:26
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    概率與統(tǒng)計(jì)(一)
    概率與統(tǒng)計(jì)(二)
    某既有隔震建筑檢測與安全性鑒定
    ApplePay橫空出世 安全性遭受質(zhì)疑 拿什么保護(hù)你,我的蘋果支付?
    綜合發(fā)力 化解災(zāi)難性醫(yī)療風(fēng)險(xiǎn)
    中國媒體在災(zāi)難性新聞報(bào)道中的改進(jìn)與不足
    新聞傳播(2015年13期)2015-07-18 11:00:41
    每晚都被弄得嗷嗷叫到高潮| 国产精品av久久久久免费| 国产成人欧美| 久久精品国产综合久久久| xxxwww97欧美| 搡老妇女老女人老熟妇| 伊人久久大香线蕉亚洲五| 欧美不卡视频在线免费观看 | 国产蜜桃级精品一区二区三区| 老司机深夜福利视频在线观看| 美女扒开内裤让男人捅视频| 中亚洲国语对白在线视频| 18禁国产床啪视频网站| 亚洲人成网站在线播放欧美日韩| 欧洲精品卡2卡3卡4卡5卡区| 亚洲黑人精品在线| 天天躁狠狠躁夜夜躁狠狠躁| 熟妇人妻久久中文字幕3abv| 国产97色在线日韩免费| 欧美激情极品国产一区二区三区| ponron亚洲| 国产熟女xx| 一边摸一边抽搐一进一小说| 亚洲片人在线观看| 亚洲五月色婷婷综合| 亚洲 欧美一区二区三区| 亚洲成av人片免费观看| АⅤ资源中文在线天堂| 国产精品久久久人人做人人爽| 亚洲激情在线av| 男女午夜视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 好看av亚洲va欧美ⅴa在| 老司机福利观看| 欧美日韩福利视频一区二区| 啪啪无遮挡十八禁网站| 欧美黑人欧美精品刺激| 国产在线观看jvid| 999久久久精品免费观看国产| 久久久久久久久免费视频了| 国产不卡一卡二| 国产在线精品亚洲第一网站| 欧美日韩一级在线毛片| 免费女性裸体啪啪无遮挡网站| 岛国视频午夜一区免费看| 成人国产一区最新在线观看| 欧美精品啪啪一区二区三区| 91老司机精品| 亚洲第一欧美日韩一区二区三区| 一本久久中文字幕| 国产三级黄色录像| 18禁国产床啪视频网站| 国产不卡一卡二| 亚洲成人久久性| 黄片播放在线免费| 国产精品一区二区三区四区久久 | 久久久久久久久免费视频了| 国产一区二区在线av高清观看| 亚洲精品一卡2卡三卡4卡5卡| 制服人妻中文乱码| 国产成人一区二区三区免费视频网站| 黑丝袜美女国产一区| 草草在线视频免费看| 亚洲欧美日韩无卡精品| 999精品在线视频| 老司机在亚洲福利影院| 天堂影院成人在线观看| 99热6这里只有精品| 久久精品91无色码中文字幕| 长腿黑丝高跟| 国产精品久久久av美女十八| 在线观看66精品国产| 日韩国内少妇激情av| 看免费av毛片| 国产av一区在线观看免费| 人人妻人人澡人人看| 1024视频免费在线观看| 制服丝袜大香蕉在线| 99精品欧美一区二区三区四区| 2021天堂中文幕一二区在线观 | 久久婷婷人人爽人人干人人爱| 免费在线观看影片大全网站| 久久精品91蜜桃| 国产精品国产高清国产av| 韩国av一区二区三区四区| 午夜久久久久精精品| 国产国语露脸激情在线看| 午夜免费激情av| 免费看十八禁软件| 亚洲欧美一区二区三区黑人| 欧美午夜高清在线| 一a级毛片在线观看| 啦啦啦免费观看视频1| 99国产极品粉嫩在线观看| 麻豆久久精品国产亚洲av| 一级黄色大片毛片| 日韩精品中文字幕看吧| 欧美国产精品va在线观看不卡| 在线观看舔阴道视频| 中文字幕av电影在线播放| 黑人欧美特级aaaaaa片| 国产精品 欧美亚洲| 免费看日本二区| 久久这里只有精品19| 精品熟女少妇八av免费久了| 亚洲欧美一区二区三区黑人| 日日干狠狠操夜夜爽| 欧美性猛交╳xxx乱大交人| 视频区欧美日本亚洲| 在线免费观看的www视频| 两个人免费观看高清视频| 黄色女人牲交| 99热只有精品国产| 听说在线观看完整版免费高清| 国产精品九九99| 日韩欧美三级三区| 日韩欧美一区视频在线观看| 亚洲精品中文字幕在线视频| 国产精品99久久99久久久不卡| 精品一区二区三区av网在线观看| 成人18禁在线播放| 99久久综合精品五月天人人| 狂野欧美激情性xxxx| 免费在线观看成人毛片| 久久久国产成人免费| 十八禁人妻一区二区| 18禁美女被吸乳视频| 午夜激情av网站| 夜夜夜夜夜久久久久| 中文字幕人妻熟女乱码| 香蕉国产在线看| 亚洲黑人精品在线| 免费搜索国产男女视频| 狠狠狠狠99中文字幕| 色播亚洲综合网| 国产一区二区三区在线臀色熟女| 国产一卡二卡三卡精品| 亚洲国产精品999在线| 男人舔奶头视频| 一区二区三区激情视频| 不卡av一区二区三区| 成人永久免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 91麻豆精品激情在线观看国产| 国产真人三级小视频在线观看| 久久精品人妻少妇| 女人高潮潮喷娇喘18禁视频| 精品电影一区二区在线| 妹子高潮喷水视频| 欧美日韩黄片免| 一二三四在线观看免费中文在| 成人国产综合亚洲| 757午夜福利合集在线观看| 久久精品亚洲精品国产色婷小说| 久久精品亚洲精品国产色婷小说| 中文字幕精品免费在线观看视频| 欧美激情久久久久久爽电影| 国产亚洲精品av在线| 亚洲在线自拍视频| 人人澡人人妻人| 国产人伦9x9x在线观看| 高清在线国产一区| 国内久久婷婷六月综合欲色啪| 国产av一区在线观看免费| 国产精品一区二区免费欧美| 中文字幕av电影在线播放| 成人一区二区视频在线观看| 国产免费男女视频| 又黄又爽又免费观看的视频| 又黄又粗又硬又大视频| 亚洲一区中文字幕在线| 日韩国内少妇激情av| 757午夜福利合集在线观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩瑟瑟在线播放| 精品高清国产在线一区| 男男h啪啪无遮挡| 午夜免费鲁丝| 观看免费一级毛片| 正在播放国产对白刺激| 特大巨黑吊av在线直播 | 老汉色∧v一级毛片| 在线观看www视频免费| 日本黄色视频三级网站网址| 少妇粗大呻吟视频| 少妇粗大呻吟视频| 午夜福利在线在线| 欧美日本亚洲视频在线播放| 在线观看一区二区三区| 午夜福利一区二区在线看| 国产精品一区二区三区四区久久 | 51午夜福利影视在线观看| 久久久久久久精品吃奶| 国产欧美日韩一区二区精品| 女警被强在线播放| 十分钟在线观看高清视频www| 亚洲精品国产一区二区精华液| 午夜久久久在线观看| 最新美女视频免费是黄的| 日韩欧美国产在线观看| 啦啦啦韩国在线观看视频| 久久久久免费精品人妻一区二区 | 久久精品91蜜桃| 亚洲 国产 在线| 精品人妻1区二区| 日韩精品免费视频一区二区三区| 岛国视频午夜一区免费看| 哪里可以看免费的av片| 精品久久久久久久末码| 亚洲人成伊人成综合网2020| 熟女少妇亚洲综合色aaa.| www.熟女人妻精品国产| 可以在线观看毛片的网站| 国产伦在线观看视频一区| 在线观看66精品国产| 身体一侧抽搐| 人人妻人人澡欧美一区二区| 在线观看一区二区三区| 1024香蕉在线观看| 国产一区二区在线av高清观看| 淫秽高清视频在线观看| 色精品久久人妻99蜜桃| 欧美成人免费av一区二区三区| 黄色a级毛片大全视频| av电影中文网址| 国产一区二区激情短视频| 99热只有精品国产| 叶爱在线成人免费视频播放| a级毛片在线看网站| 麻豆国产av国片精品| 国产成人欧美在线观看| 黄色丝袜av网址大全| 91九色精品人成在线观看| 亚洲精品中文字幕在线视频| 亚洲国产精品999在线| 亚洲久久久国产精品| 欧美成人性av电影在线观看| 亚洲五月婷婷丁香| 18禁黄网站禁片午夜丰满| 久久久久国产一级毛片高清牌| 怎么达到女性高潮| 成年版毛片免费区| 国内少妇人妻偷人精品xxx网站 | 高清毛片免费观看视频网站| 麻豆国产av国片精品| 特大巨黑吊av在线直播 | 一区福利在线观看| АⅤ资源中文在线天堂| 99国产精品一区二区蜜桃av| 18禁黄网站禁片午夜丰满| 午夜成年电影在线免费观看| 久久国产亚洲av麻豆专区| 色老头精品视频在线观看| 国产高清视频在线播放一区| 美女大奶头视频| 黄片小视频在线播放| 大型av网站在线播放| 国产成人av教育| 九色国产91popny在线| 在线国产一区二区在线| 一边摸一边做爽爽视频免费| 人人妻,人人澡人人爽秒播| 国产精品,欧美在线| 免费看美女性在线毛片视频| tocl精华| 亚洲男人的天堂狠狠| 亚洲精品av麻豆狂野| 侵犯人妻中文字幕一二三四区| 黄色 视频免费看| 无遮挡黄片免费观看| 久久久久久久久久黄片| 久久中文字幕一级| 一本一本综合久久| 日本三级黄在线观看| 美国免费a级毛片| 妹子高潮喷水视频| 日本 av在线| 精品国产乱码久久久久久男人| 国产熟女午夜一区二区三区| 久久伊人香网站| 中文字幕精品免费在线观看视频| 国产高清视频在线播放一区| 国产主播在线观看一区二区| 午夜福利免费观看在线| 成人手机av| 久久久久国产一级毛片高清牌| 国产精品野战在线观看| 欧美激情久久久久久爽电影| 久久香蕉精品热| 亚洲三区欧美一区| 天堂影院成人在线观看| 日韩免费av在线播放| 长腿黑丝高跟| 精品第一国产精品| 亚洲 欧美一区二区三区| 在线十欧美十亚洲十日本专区| 999久久久国产精品视频| 18禁国产床啪视频网站| 51午夜福利影视在线观看| 在线观看日韩欧美| 成人18禁高潮啪啪吃奶动态图| 亚洲精品一卡2卡三卡4卡5卡| 搡老岳熟女国产| 亚洲狠狠婷婷综合久久图片| 中文字幕人妻熟女乱码| 人人澡人人妻人| 欧美精品亚洲一区二区| 村上凉子中文字幕在线| 国产午夜精品久久久久久| 后天国语完整版免费观看| 国产成人啪精品午夜网站| 亚洲最大成人中文| 婷婷精品国产亚洲av在线| 国产精品日韩av在线免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| www.熟女人妻精品国产| 成人国产综合亚洲| 欧美一区二区精品小视频在线| 欧美日韩中文字幕国产精品一区二区三区| a在线观看视频网站| 欧美成人免费av一区二区三区| 久久精品亚洲精品国产色婷小说| 国产精品影院久久| 精品乱码久久久久久99久播| а√天堂www在线а√下载| 在线av久久热| 中文资源天堂在线| www.自偷自拍.com| 久久香蕉国产精品| 午夜精品久久久久久毛片777| 中文字幕最新亚洲高清| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 满18在线观看网站| 十分钟在线观看高清视频www| 亚洲av成人一区二区三| 日韩欧美 国产精品| 亚洲成av片中文字幕在线观看| 国产色视频综合| 国产片内射在线| 777久久人妻少妇嫩草av网站| 女人高潮潮喷娇喘18禁视频| 亚洲成av片中文字幕在线观看| 久久久久国内视频| www.熟女人妻精品国产| 久久久久国产精品人妻aⅴ院| 亚洲国产毛片av蜜桃av| 淫秽高清视频在线观看| 欧美日韩一级在线毛片| 窝窝影院91人妻| 国产精品综合久久久久久久免费| 天堂动漫精品| 午夜精品在线福利| 精品国产乱子伦一区二区三区| 白带黄色成豆腐渣| 丁香欧美五月| 丝袜美腿诱惑在线| e午夜精品久久久久久久| 91在线观看av| 91成人精品电影| 制服丝袜大香蕉在线| 欧美中文日本在线观看视频| 搡老熟女国产l中国老女人| 麻豆久久精品国产亚洲av| 黄色视频不卡| 精品欧美国产一区二区三| 久久国产精品人妻蜜桃| 女同久久另类99精品国产91| 欧美午夜高清在线| 这个男人来自地球电影免费观看| 九色国产91popny在线| 久久婷婷成人综合色麻豆| 久久久久久久久久黄片| 久9热在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品久久久久久毛片| 久久天堂一区二区三区四区| 午夜福利在线观看吧| 高清在线国产一区| 国产精品亚洲一级av第二区| 亚洲,欧美精品.| 又大又爽又粗| 国产高清videossex| 侵犯人妻中文字幕一二三四区| 91麻豆精品激情在线观看国产| 视频在线观看一区二区三区| 久久久久精品国产欧美久久久| 黄色成人免费大全| 国产一卡二卡三卡精品| 国产精华一区二区三区| 啦啦啦韩国在线观看视频| 悠悠久久av| 国产精品自产拍在线观看55亚洲| 国产极品粉嫩免费观看在线| 窝窝影院91人妻| www.自偷自拍.com| 露出奶头的视频| av片东京热男人的天堂| 美国免费a级毛片| 在线观看一区二区三区| 免费人成视频x8x8入口观看| 国产亚洲欧美精品永久| 老司机深夜福利视频在线观看| 精品国产一区二区三区四区第35| 男女视频在线观看网站免费 | 俺也久久电影网| 男女那种视频在线观看| 亚洲国产高清在线一区二区三 | 日本a在线网址| 十分钟在线观看高清视频www| 日本成人三级电影网站| 变态另类丝袜制服| 99久久国产精品久久久| 女人爽到高潮嗷嗷叫在线视频| 日韩 欧美 亚洲 中文字幕| 免费看十八禁软件| 日日爽夜夜爽网站| 久久久国产欧美日韩av| 亚洲七黄色美女视频| 特大巨黑吊av在线直播 | 俄罗斯特黄特色一大片| 国产精品久久电影中文字幕| 中文亚洲av片在线观看爽| 国产欧美日韩一区二区三| 精品久久久久久久末码| 高潮久久久久久久久久久不卡| 日韩免费av在线播放| 91av网站免费观看| 黄色丝袜av网址大全| 欧美三级亚洲精品| 国产91精品成人一区二区三区| 一级作爱视频免费观看| 窝窝影院91人妻| 悠悠久久av| 一区福利在线观看| 亚洲av熟女| 精品欧美一区二区三区在线| 精品国内亚洲2022精品成人| 99国产极品粉嫩在线观看| 免费在线观看影片大全网站| 亚洲国产高清在线一区二区三 | 操出白浆在线播放| www日本在线高清视频| 国产精品二区激情视频| 免费看日本二区| 中文亚洲av片在线观看爽| 免费观看人在逋| 黄片播放在线免费| 国产在线观看jvid| 国产色视频综合| 不卡av一区二区三区| АⅤ资源中文在线天堂| 熟女少妇亚洲综合色aaa.| 身体一侧抽搐| 亚洲人成网站在线播放欧美日韩| 欧美日韩中文字幕国产精品一区二区三区| av片东京热男人的天堂| 熟女电影av网| 制服丝袜大香蕉在线| 亚洲在线自拍视频| 精华霜和精华液先用哪个| 1024手机看黄色片| 亚洲最大成人中文| 亚洲精品粉嫩美女一区| √禁漫天堂资源中文www| 日本免费一区二区三区高清不卡| 国产精品 国内视频| 99久久99久久久精品蜜桃| 国产一卡二卡三卡精品| 久久久久久亚洲精品国产蜜桃av| 怎么达到女性高潮| 97人妻精品一区二区三区麻豆 | 人妻久久中文字幕网| 亚洲国产欧美日韩在线播放| 精品久久久久久久毛片微露脸| 观看免费一级毛片| 亚洲自拍偷在线| 桃红色精品国产亚洲av| 欧美日韩瑟瑟在线播放| 日本精品一区二区三区蜜桃| 亚洲成人久久爱视频| 激情在线观看视频在线高清| 欧美黑人欧美精品刺激| 久久草成人影院| 国产亚洲精品久久久久久毛片| 中文字幕人妻丝袜一区二区| 精品久久蜜臀av无| 国产一级毛片七仙女欲春2 | 正在播放国产对白刺激| 国产精品久久久久久亚洲av鲁大| а√天堂www在线а√下载| 久久伊人香网站| 国产爱豆传媒在线观看 | 国产亚洲欧美在线一区二区| 亚洲免费av在线视频| 婷婷精品国产亚洲av| 久久久久久大精品| 法律面前人人平等表现在哪些方面| 久久精品人妻少妇| 最近最新免费中文字幕在线| 色综合婷婷激情| 亚洲激情在线av| 日本撒尿小便嘘嘘汇集6| 精品国产亚洲在线| av在线天堂中文字幕| 国产高清videossex| 国产精品久久久av美女十八| 亚洲色图 男人天堂 中文字幕| 久久精品亚洲精品国产色婷小说| 亚洲精品色激情综合| 成人18禁在线播放| 午夜影院日韩av| 午夜福利在线观看吧| 国产熟女xx| 天堂动漫精品| 国产精品永久免费网站| 97人妻精品一区二区三区麻豆 | 日韩欧美国产在线观看| 日本五十路高清| 别揉我奶头~嗯~啊~动态视频| 亚洲五月色婷婷综合| 99国产精品一区二区蜜桃av| 精品国产国语对白av| 99精品欧美一区二区三区四区| √禁漫天堂资源中文www| 色在线成人网| x7x7x7水蜜桃| 日韩免费av在线播放| 精品欧美国产一区二区三| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品成人综合色| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩中文字幕国产精品一区二区三区| 午夜久久久久精精品| 久热这里只有精品99| 99久久综合精品五月天人人| 2021天堂中文幕一二区在线观 | 又紧又爽又黄一区二区| av视频在线观看入口| 搡老熟女国产l中国老女人| 一二三四社区在线视频社区8| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一av免费看| 两人在一起打扑克的视频| 欧美+亚洲+日韩+国产| 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区视频在线观看免费| 午夜福利在线观看吧| 国产乱人伦免费视频| 成人午夜高清在线视频 | 成人免费观看视频高清| 啦啦啦观看免费观看视频高清| 国产精品免费一区二区三区在线| svipshipincom国产片| 中文字幕人成人乱码亚洲影| 免费看日本二区| 欧美av亚洲av综合av国产av| 最新在线观看一区二区三区| 男女午夜视频在线观看| 免费在线观看成人毛片| 此物有八面人人有两片| 国产麻豆成人av免费视频| 欧美一区二区精品小视频在线| 麻豆成人午夜福利视频| 身体一侧抽搐| 在线观看日韩欧美| 精品不卡国产一区二区三区| 午夜福利在线观看吧| 国产精品一区二区精品视频观看| 一夜夜www| 欧美绝顶高潮抽搐喷水| 日本免费a在线| 国产亚洲精品综合一区在线观看 | 久久精品国产清高在天天线| 久久草成人影院| 亚洲中文字幕日韩| 欧美成人性av电影在线观看| 亚洲 欧美一区二区三区| 国产又爽黄色视频| 热99re8久久精品国产| 伊人久久大香线蕉亚洲五| 久久香蕉激情| 999久久久国产精品视频| 可以在线观看的亚洲视频| 一区二区三区高清视频在线| 久9热在线精品视频| 亚洲av成人不卡在线观看播放网| 日韩av在线大香蕉| 变态另类成人亚洲欧美熟女| 久久精品91无色码中文字幕| 国产精品免费视频内射| 一个人观看的视频www高清免费观看 | 精品第一国产精品| 日韩欧美三级三区| 日韩大尺度精品在线看网址| 美女高潮喷水抽搐中文字幕| 我的亚洲天堂| 国产一区二区在线av高清观看| 一级毛片精品| or卡值多少钱| 亚洲自偷自拍图片 自拍| 满18在线观看网站| 欧美日韩精品网址| 亚洲精华国产精华精| 国产精品99久久99久久久不卡| 亚洲精华国产精华精| 久久精品91蜜桃| 一本一本综合久久| 伊人久久大香线蕉亚洲五| 欧美色欧美亚洲另类二区| 一级a爱片免费观看的视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品,欧美在线| 国产精品 欧美亚洲| 成年免费大片在线观看|