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

    基于方差分析的航班離地地速影響因素研究

    2021-11-17 00:41:30陳文瑛副教授傅佳萍
    安全 2021年10期
    關(guān)鍵詞:離地因變量分組

    陳文瑛副教授 傅佳萍 齊 凱,3

    (1.首都經(jīng)濟(jì)貿(mào)易大學(xué) 管理工程學(xué)院,北京 100070;2.城市群系統(tǒng)演化與發(fā)展的決策模擬研究北京市重點(diǎn)實(shí)驗(yàn)室,北京 100071;3.山東航空股份有限公司,山東 濟(jì)南 250000)

    0 引言

    民航飛行安全關(guān)系到旅客生命和財(cái)產(chǎn)的安危,因此對(duì)飛行安全影響因素進(jìn)行研究是非常有必要的。起飛階段是飛行過(guò)程的關(guān)鍵階段之一,飛機(jī)在該階段容易發(fā)生沖出/偏出跑道事故,對(duì)飛機(jī)和機(jī)上人員的生命構(gòu)成重大威脅。據(jù)統(tǒng)計(jì),在一次飛行中,起飛階段、最后進(jìn)近和著陸階段時(shí)間約只占總飛行時(shí)間的6%,但有近50%的事故,33%的致命事故和22%的死亡人數(shù)發(fā)生在此階段,而飛機(jī)沖偏出跑道是起降階段事故發(fā)生率較高的事故之一。由此國(guó)內(nèi)外學(xué)者對(duì)飛機(jī)沖偏出跑道的原因進(jìn)行分析研究,如:國(guó)內(nèi)學(xué)者孔祥駿、楊開(kāi)等分別利用歷史數(shù)據(jù)從“人—機(jī)—環(huán)”角度,發(fā)現(xiàn)導(dǎo)致沖偏出跑道事故主要原因是飛行員操作不當(dāng)、氣象條件復(fù)雜、機(jī)組管理等,為研究沖偏出跑道事故風(fēng)險(xiǎn)因子間的關(guān)聯(lián)性奠定基礎(chǔ);張曉全等在國(guó)內(nèi)首次應(yīng)用尋因分析法對(duì)沖偏出跑道事故進(jìn)行危險(xiǎn)源辨識(shí),發(fā)現(xiàn)事故危險(xiǎn)源包括天氣原因、機(jī)組資源等方面問(wèn)題;趙寧寧等利用事故樹(shù)和貝葉斯網(wǎng)絡(luò)模型對(duì)飛機(jī)偏沖出跑道事故的風(fēng)險(xiǎn)進(jìn)行分析,發(fā)現(xiàn)“飛行員操作有誤”發(fā)生概率較高,其次分別為“飛行員判斷有誤”“飛機(jī)故障”“飛機(jī)相關(guān)系統(tǒng)突然失效”“跑道摩擦系數(shù)較小”和“跑道道面濕滑”;王浩寧等基于系統(tǒng)理論的事故模型,發(fā)現(xiàn)在飛機(jī)沖偏出跑道事故中機(jī)組人員可能存在不安全控制行為。國(guó)外學(xué)者M(jìn)arisa Jenkins研究指出風(fēng)、道面情況和起降距離是飛機(jī)沖偏出跑道事故的主要原因;Distefano等創(chuàng)新性地采用多重對(duì)應(yīng)分析法對(duì)沖偏出跑道事故進(jìn)行分析,發(fā)現(xiàn)起飛階段發(fā)生該類事故最常見(jiàn)的原因是飛機(jī)系統(tǒng)故障;Pramono等分析事故報(bào)告,第一次調(diào)查印尼航空事故的因素,指出天氣、機(jī)組資源管理、跑道條件技術(shù)故障等是導(dǎo)致沖偏出跑道的原因。大多文獻(xiàn)通過(guò)歷史事故數(shù)據(jù)或飛機(jī)構(gòu)造分析探究發(fā)生沖偏出跑道事故原因,研究發(fā)現(xiàn)機(jī)組資源管理對(duì)飛機(jī)沖偏出跑道有影響。輪胎爆胎是機(jī)組資源管理失效的一種情況,而導(dǎo)致輪胎爆胎原因很多,其中起飛時(shí)離地地速過(guò)大是導(dǎo)致輪胎爆胎的重要誘因之一。

    目前,較少有研究基于飛行數(shù)據(jù)分析該類事故的風(fēng)險(xiǎn)因素,本文根據(jù)專家意見(jiàn)及快速存取記錄器所采集的機(jī)型為B738的民航飛機(jī)飛行數(shù)據(jù)信息,應(yīng)用方差分析方法,從機(jī)重、溫度和人的操作3個(gè)角度入手,識(shí)別飛機(jī)起飛離地地速的影響因素,研究這些因素對(duì)飛機(jī)離地地速的影響程度,可為降低飛機(jī)沖偏出跑道事故風(fēng)險(xiǎn)提供科學(xué)依據(jù)。

    1 研究方法

    方差分析用于分析研究不同來(lái)源的變化因子對(duì)總體變化的貢獻(xiàn)能力,進(jìn)而確定變化因子對(duì)研究結(jié)果的影響大小。本文應(yīng)用這種方法考察不同風(fēng)險(xiǎn)因子對(duì)飛機(jī)離地地速是否有影響及影響的程度。

    方差分析是用于檢驗(yàn)不同水平下某因變量的均值變化是否顯著的一種定量分析方法,包括單因素方差分析、多因素方差分析等,利用此方法可以分析出某一因素是否為因變量的影響因素。因數(shù)據(jù)的限制,本文僅通過(guò)單因素方差分析來(lái)識(shí)別導(dǎo)致飛機(jī)離地地速

    v

    過(guò)大的風(fēng)險(xiǎn)因子。單因素方差分析,也稱一維方差分析,是用來(lái)研究一個(gè)控制變量的不同水平是否對(duì)觀測(cè)變量產(chǎn)生顯著影響的方法,換句話說(shuō),就是檢驗(yàn)3個(gè)或3個(gè)以上獨(dú)立樣本觀察值的各組平均數(shù)是否相等。其對(duì)數(shù)據(jù)的要求為:觀測(cè)變量為數(shù)值型變量,控制變量為分類型變量,分類型變量的每一種具體表現(xiàn)形式稱為控制變量的水平。除此之外,在使用方差分析前,各組觀測(cè)值必須滿足正態(tài)性、獨(dú)立性及方差齊性3個(gè)條件。對(duì)于單因素方差分析結(jié)果,可通過(guò)關(guān)聯(lián)強(qiáng)度指數(shù)

    ω

    衡量自變量的不同因子水平對(duì)因變量的貢獻(xiàn)程度,關(guān)聯(lián)強(qiáng)度計(jì)算公式:

    (1)

    式中:

    ω

    —關(guān)聯(lián)強(qiáng)度指標(biāo);

    k

    —自變量因子水平;

    MSE—組內(nèi)均方值;

    SSA—組間平方和;

    SSE—組內(nèi)平方和;

    SST—總離差平方和,SST=SSA+SSE。

    關(guān)聯(lián)強(qiáng)度指數(shù)

    ω

    高低判斷自變量與因變量關(guān)聯(lián)程度的標(biāo)準(zhǔn):關(guān)聯(lián)強(qiáng)度指數(shù)

    ω

    在6%以下,屬于微弱關(guān)系;關(guān)聯(lián)強(qiáng)度指數(shù)

    ω

    位于6%~16%之間,屬于中度關(guān)系;關(guān)聯(lián)強(qiáng)度指數(shù)

    ω

    大于16%,屬于強(qiáng)度關(guān)系。

    本文采用SPSS軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,從而得到方差分析結(jié)果,以分析機(jī)重、溫度和人的操作對(duì)飛機(jī)離地地速的影響。

    2 數(shù)據(jù)來(lái)源及預(yù)處理

    輪胎爆炸可能導(dǎo)致沖偏出跑道事故,而飛機(jī)離地地速

    v

    過(guò)大是飛機(jī)爆胎的重要誘因之一,本文主要探究導(dǎo)致輪胎爆胎因素間的關(guān)系,即溫度

    T

    、機(jī)重

    G

    和人為操作對(duì)飛機(jī)離地地速

    ν

    的影響,如圖1。在8位航空安全管理和飛行品質(zhì)管理方面專家的指導(dǎo)下,選定起飛過(guò)程的抬輪時(shí)間差

    t

    、離地時(shí)仰角

    θ

    、發(fā)動(dòng)機(jī)低壓轉(zhuǎn)子最大轉(zhuǎn)速百分比

    N

    、最大仰角增加率

    α

    、平均仰角增加率

    β

    作為表征人為操作的參數(shù)。由此,提取2019年7月246架次飛機(jī)離地瞬間的246條數(shù)據(jù)。數(shù)據(jù)包括

    v

    G

    、

    T

    、

    t

    、

    θ

    N

    、

    α

    β

    等,由于篇幅有限,僅列出部分?jǐn)?shù)據(jù), 見(jiàn)表1。由于采集的是單機(jī)場(chǎng)的數(shù)據(jù),所以未考慮標(biāo)準(zhǔn)氣壓對(duì)

    ν

    的影響。

    t

    是實(shí)際抬輪時(shí)刻與理論抬輪時(shí)刻的差值。

    圖1 可能導(dǎo)致飛機(jī)爆胎進(jìn)而沖偏出跑道因素間的邏輯影響關(guān)系

    表1 某航空公司2019年7月部分飛行紀(jì)錄

    本文將

    ν

    作為因變量,

    G

    、

    T

    、

    t

    、

    θ

    、

    N

    α

    、

    β

    作為自變量,觀察數(shù)據(jù)發(fā)現(xiàn),因變量和自變量都為數(shù)值變量。為把

    T

    等自變量數(shù)據(jù)轉(zhuǎn)化為分類變量,需要對(duì)自變量數(shù)據(jù)進(jìn)行分組。數(shù)據(jù)分組的方法有單變量值分組和組距分組2種。單變量值分組是把每一個(gè)變量值作為一組,這種分組通常只適合離散變量,且在變量值較少的情況下使用。在連續(xù)變量或變量值較多的情況下,通常采用組距分組,即將全部變量依次劃分為若干個(gè)區(qū)間,并將一個(gè)區(qū)間的變量值作為一組。根據(jù)數(shù)據(jù)特征,這里對(duì)

    T

    進(jìn)行單變量值分組,對(duì)

    G

    t

    等進(jìn)行組距分組。同一溫度下的觀察值劃為一組,共有11個(gè)不同的溫度,即11個(gè)組;對(duì)于

    G

    、

    t

    等自變量數(shù)據(jù),這里通過(guò)斯特杰斯經(jīng)驗(yàn)公式,對(duì)數(shù)據(jù)進(jìn)行分組,經(jīng)驗(yàn)公式:

    (2)

    式中:

    C

    —組距;

    n

    —總次數(shù);

    Max—最大值;

    Min—最小值。

    通過(guò)式(2)確定分組所需組距和組數(shù),再對(duì)數(shù)據(jù)進(jìn)行分組。由于對(duì)數(shù)據(jù)單因素方差分析時(shí),需要每組數(shù)據(jù)有3個(gè)或3個(gè)以上的觀察值,所以得到的分組結(jié)果若不滿足該條件,則應(yīng)去掉該組數(shù)據(jù)后重新分組。例如,對(duì)

    G

    的數(shù)據(jù)進(jìn)行分組:首先,將最大值72 974、最小值62 378以及總觀測(cè)次數(shù)246代入式(2)中進(jìn)行計(jì)算,得到組距為1 187.23,組數(shù)10組,之后發(fā)現(xiàn)位于第一組的數(shù)據(jù)只有一條,分組結(jié)果不滿足單因素方差分析的條件,原因是該條數(shù)據(jù)與其他數(shù)據(jù)間的差值較大。通過(guò)對(duì)數(shù)據(jù)進(jìn)行刪減后重新分組,將

    G

    的變量值共分為9組且每組都包含3條及3條以上的數(shù)據(jù);其他數(shù)據(jù)的分組過(guò)程類似,具體分組結(jié)果,見(jiàn)表2。

    表2 自變量數(shù)據(jù)分組情況

    3 結(jié)果分析

    3.1 正態(tài)性及方差齊性檢驗(yàn)

    方差分析需要對(duì)因變量的觀察值進(jìn)行正態(tài)及方差齊性檢驗(yàn)。首先,本文采用分位數(shù)—分位數(shù)圖(Quantile-quantile Plot,Q-Q)來(lái)判斷是否滿足正態(tài)分布。Q-Q圖的縱坐標(biāo)為預(yù)期的正態(tài)值,橫坐標(biāo)為樣本點(diǎn)的實(shí)際取值。如果樣本集合的概率分布與正態(tài)分布十分相似,那么Q-Q圖逼近一條直線。對(duì)因變量觀察值(即

    ν

    的變量值)進(jìn)行正態(tài)檢驗(yàn),如圖2。由圖2可見(jiàn),

    ν

    數(shù)據(jù)擬合后近似為一條直線,說(shuō)明數(shù)據(jù)近似服從正態(tài)分布。

    圖2 ν的正態(tài)Q-Q圖

    然后,分別對(duì)在

    G

    T

    、

    t

    θ

    、

    N

    、

    α

    、

    β

    下的

    v

    的觀察值進(jìn)行方差齊性檢驗(yàn),從而判斷數(shù)據(jù)是否滿足方差齊性,見(jiàn)表3。

    表3 自變量觀察值的方差齊性檢驗(yàn)結(jié)果

    由表3可知,在方差齊性分析時(shí),計(jì)算所得

    P

    值都大于0.05,說(shuō)明滿足方差齊性條件,即可進(jìn)行單因素方差分析。

    3.2 結(jié)果分析

    分別對(duì)

    G

    T

    、

    t

    θ

    、

    N

    α

    、

    β

    等因素對(duì)

    ν

    的影響進(jìn)行單因素方差分析,見(jiàn)表4。由表4結(jié)果可得,

    G

    、

    t

    、

    N

    所得

    P

    值都小于0.01,說(shuō)明這些因素對(duì)

    v

    的影響極其顯著;

    θ

    P

    值小于0.05,說(shuō)明

    θ

    對(duì)

    v

    的影響顯著;而

    T

    、

    α

    、

    β

    P

    值大于0.05,說(shuō)明這3個(gè)因素對(duì)

    v

    的影響不顯著。在反映人為操作的5個(gè)因素中,有3個(gè)因素對(duì)

    v

    的影響顯著,說(shuō)明人為操作對(duì)飛機(jī)離地地速的影響顯著。進(jìn)一步分析結(jié)果可以發(fā)現(xiàn),

    v

    的總變異量由

    G

    、

    t

    、

    θ

    、

    N

    解釋的百分比分別為22.44%、6.05%、3.88%、10.11%(

    ω

    =0.224 4、0.060 5、0.038 8、0.101 1)。其中,

    θ

    v

    的關(guān)聯(lián)強(qiáng)度指數(shù)為3.88%,小于6%,說(shuō)明其與

    v

    的關(guān)聯(lián)強(qiáng)度屬于微弱關(guān)系;

    t

    、

    N

    v

    的關(guān)聯(lián)強(qiáng)度指數(shù)分別為6.05%、10.11%,位于6%~16%之間,說(shuō)明這2個(gè)因素與

    v

    的關(guān)聯(lián)強(qiáng)度屬于中度關(guān)系;

    G

    v

    的關(guān)聯(lián)強(qiáng)度指數(shù)為22.44%,大于16%,說(shuō)明

    G

    v

    的關(guān)聯(lián)強(qiáng)度屬于強(qiáng)度關(guān)系。

    表4 不同自變量下因變量的單因素方差分析結(jié)果

    為進(jìn)一步了解

    v

    的觀察值均值間的差異,分別繪制不同

    G

    t

    、

    θ

    、

    N

    水平下觀察值的均值圖,如圖3。圖3可知,飛機(jī)越重,

    v

    越大;當(dāng)抬輪時(shí)間差較大時(shí),

    v

    偏大,說(shuō)明相較理論抬輪時(shí)間,若實(shí)際抬輪時(shí)間越晚,更可能導(dǎo)致

    v

    過(guò)大;當(dāng)

    θ

    偏小時(shí),

    v

    偏大;

    N

    偏大時(shí),

    v

    偏大。

    圖3 飛機(jī)重量與人為操作因素的飛機(jī)離地地速v均值圖

    可見(jiàn),經(jīng)過(guò)單因素方差分析,最終發(fā)現(xiàn)飛機(jī)重量

    G

    、人為操作(如帶桿和加油門)對(duì)民航飛機(jī)離地地速有顯著性影響。其中,在表征人為操作的因素中,

    t

    在起飛過(guò)程中,抬輪時(shí)間晚于理論值時(shí),會(huì)導(dǎo)致

    v

    偏大;若離地時(shí)仰角偏小時(shí),容易導(dǎo)致

    v

    偏大;

    N

    值偏大時(shí),

    v

    也偏大。上述結(jié)果說(shuō)明飛機(jī)重量、人為操作是飛機(jī)輪胎爆胎的間接誘因,在趙安家等對(duì)輪胎構(gòu)造以及常見(jiàn)故障分析中有所體現(xiàn)。

    4 結(jié)論

    本文分別對(duì)

    G

    T

    、

    t

    θ

    、

    N

    、

    α

    、

    β

    下的飛機(jī)離地地速觀測(cè)值進(jìn)行單因素方差分析。其中,

    t

    θ

    、

    N

    α

    、

    β

    表征人為操作,通過(guò)分析這些因素與

    v

    的關(guān)系,反映人為操作對(duì)

    v

    的影響,最終得到如下結(jié)論:(1)運(yùn)用單因素方差分析,在0.05顯著性水平上,分析得到有顯著性差異的影響因素為

    G

    t

    、

    θ

    、

    N

    ,即機(jī)重和人為操作(帶桿、加油門)對(duì)

    v

    影響顯著。(2)對(duì)單因素方差分析結(jié)果表現(xiàn)為顯著性的因素進(jìn)一步分析發(fā)現(xiàn):在不考慮其他條件(如發(fā)動(dòng)機(jī)推力、控制系統(tǒng)響應(yīng)等)下,飛機(jī)越重,

    v

    越大;在表征人為操作的因素中,若在起飛過(guò)程中,抬輪時(shí)間晚于理論值時(shí),會(huì)導(dǎo)致

    v

    偏大,若離地時(shí)仰角偏小時(shí),也容易導(dǎo)致

    v

    偏大,而就

    N

    而言,其值偏大時(shí),

    v

    偏大。

    飛機(jī)起飛過(guò)程實(shí)際上是一個(gè)動(dòng)態(tài)、多因素影響的優(yōu)化問(wèn)題,由于數(shù)據(jù)限制,本文僅通過(guò)單因素方差分析簡(jiǎn)單地得到機(jī)重和人為操作(帶桿、加油門)對(duì)離地地速有顯著性影響,希望對(duì)民航航班機(jī)組的實(shí)際操縱具有指導(dǎo)意義。

    猜你喜歡
    離地因變量分組
    調(diào)整有限因變量混合模型在藥物經(jīng)濟(jì)學(xué)健康效用量表映射中的運(yùn)用
    BJ40 刀鋒英雄版
    汽車觀察(2021年4期)2021-05-10 04:08:08
    適應(yīng)性回歸分析(Ⅳ)
    ——與非適應(yīng)性回歸分析的比較
    分組搭配
    怎么分組
    偏最小二乘回歸方法
    淺析豎直平面內(nèi)圓周運(yùn)動(dòng)中的軌道“離地”臨界條件
    物理教師(2017年6期)2017-06-27 03:23:48
    分組
    關(guān)于物體在豎直平面內(nèi)做圓周運(yùn)動(dòng)導(dǎo)致底座離地的極值問(wèn)題分析お
    空中婚禮
    日本午夜av视频| 日韩 亚洲 欧美在线| 99久久人妻综合| 熟女av电影| av网站免费在线观看视频| 99热国产这里只有精品6| 国产精品女同一区二区软件| a级毛色黄片| 国产高清三级在线| 日本黄大片高清| 亚洲欧洲国产日韩| 91久久精品国产一区二区三区| 久久久久网色| 久久精品熟女亚洲av麻豆精品| 狠狠婷婷综合久久久久久88av| 99热国产这里只有精品6| 精品少妇久久久久久888优播| 免费黄色在线免费观看| 亚洲少妇的诱惑av| 狠狠婷婷综合久久久久久88av| 亚洲国产精品专区欧美| 午夜福利影视在线免费观看| 免费少妇av软件| 免费女性裸体啪啪无遮挡网站| 亚洲伊人久久精品综合| 亚洲精品国产av成人精品| 韩国高清视频一区二区三区| 国产乱人偷精品视频| 看免费成人av毛片| 日本色播在线视频| av一本久久久久| 天美传媒精品一区二区| 丝袜脚勾引网站| 一二三四在线观看免费中文在 | 国产一区二区在线观看日韩| 在线观看免费高清a一片| 国产一级毛片在线| 看十八女毛片水多多多| 国产精品 国内视频| 久久99精品国语久久久| 一个人免费看片子| 国产成人一区二区在线| 午夜福利在线观看免费完整高清在| 亚洲一码二码三码区别大吗| 国产黄频视频在线观看| 日本av免费视频播放| 99久久中文字幕三级久久日本| av天堂久久9| 91精品伊人久久大香线蕉| 国产av国产精品国产| 成人影院久久| 欧美性感艳星| 人人妻人人添人人爽欧美一区卜| 一级毛片我不卡| 99久久精品国产国产毛片| 一级黄片播放器| 日韩欧美一区视频在线观看| av网站免费在线观看视频| 免费久久久久久久精品成人欧美视频 | 搡老乐熟女国产| 亚洲欧美日韩另类电影网站| a级片在线免费高清观看视频| 男女边摸边吃奶| 男女高潮啪啪啪动态图| 婷婷色av中文字幕| 久久国产精品大桥未久av| 少妇人妻久久综合中文| 老司机影院成人| 亚洲第一av免费看| 永久免费av网站大全| 欧美精品国产亚洲| 欧美日韩av久久| 两个人免费观看高清视频| 22中文网久久字幕| 日韩成人av中文字幕在线观看| 国产无遮挡羞羞视频在线观看| 乱人伦中国视频| 久久久久国产网址| 国产成人91sexporn| 99九九在线精品视频| av在线app专区| 国产有黄有色有爽视频| 九色成人免费人妻av| 美女国产高潮福利片在线看| 亚洲,欧美,日韩| 亚洲天堂av无毛| 啦啦啦中文免费视频观看日本| 亚洲国产成人一精品久久久| 熟女电影av网| 久久久精品94久久精品| 啦啦啦啦在线视频资源| 国产欧美亚洲国产| 各种免费的搞黄视频| 热99国产精品久久久久久7| 在线观看一区二区三区激情| 婷婷色麻豆天堂久久| 国产精品一区二区在线不卡| 纯流量卡能插随身wifi吗| av有码第一页| 黑人高潮一二区| 亚洲国产色片| 99re6热这里在线精品视频| 69精品国产乱码久久久| 夜夜骑夜夜射夜夜干| 日韩视频在线欧美| 日韩人妻精品一区2区三区| 91久久精品国产一区二区三区| 日韩一区二区三区影片| 日韩,欧美,国产一区二区三区| 最新中文字幕久久久久| 欧美日韩国产mv在线观看视频| 久久久a久久爽久久v久久| 永久免费av网站大全| 国产片内射在线| 啦啦啦啦在线视频资源| 国产白丝娇喘喷水9色精品| 美女国产视频在线观看| 欧美日韩av久久| 日韩成人伦理影院| 亚洲国产av影院在线观看| 亚洲欧美清纯卡通| 宅男免费午夜| 亚洲综合色惰| 狂野欧美激情性bbbbbb| 在线精品无人区一区二区三| 少妇的丰满在线观看| 丝袜脚勾引网站| 亚洲,欧美精品.| 一区在线观看完整版| 亚洲婷婷狠狠爱综合网| 国产午夜精品一二区理论片| 日本av免费视频播放| av卡一久久| 免费黄频网站在线观看国产| 18在线观看网站| 亚洲 欧美一区二区三区| 老熟女久久久| 精品人妻偷拍中文字幕| 免费大片18禁| 精品亚洲乱码少妇综合久久| 婷婷色麻豆天堂久久| 久久97久久精品| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看| 91久久精品国产一区二区三区| 我要看黄色一级片免费的| 亚洲精品一区蜜桃| 母亲3免费完整高清在线观看 | 午夜福利视频精品| h视频一区二区三区| 日本av免费视频播放| 王馨瑶露胸无遮挡在线观看| 国产精品国产av在线观看| 欧美日韩视频精品一区| 亚洲四区av| 亚洲精华国产精华液的使用体验| 亚洲伊人久久精品综合| 水蜜桃什么品种好| 波野结衣二区三区在线| 搡老乐熟女国产| 亚洲美女搞黄在线观看| freevideosex欧美| 又黄又粗又硬又大视频| 亚洲精品一二三| 亚洲美女视频黄频| 人妻系列 视频| 国产白丝娇喘喷水9色精品| 中文精品一卡2卡3卡4更新| a级毛色黄片| 欧美最新免费一区二区三区| 久久人人爽av亚洲精品天堂| av黄色大香蕉| 国产黄频视频在线观看| 精品一区二区三区视频在线| 看十八女毛片水多多多| 久久午夜综合久久蜜桃| 国产xxxxx性猛交| 国产精品偷伦视频观看了| 亚洲高清免费不卡视频| 亚洲精品一二三| 日韩制服骚丝袜av| 麻豆精品久久久久久蜜桃| 成人免费观看视频高清| 久久ye,这里只有精品| 内地一区二区视频在线| 欧美亚洲 丝袜 人妻 在线| 美女内射精品一级片tv| a级片在线免费高清观看视频| 亚洲欧洲国产日韩| 亚洲人与动物交配视频| 香蕉丝袜av| 国产亚洲最大av| 老熟女久久久| 成人黄色视频免费在线看| 国语对白做爰xxxⅹ性视频网站| 热99久久久久精品小说推荐| 亚洲一级一片aⅴ在线观看| 黄网站色视频无遮挡免费观看| 日韩熟女老妇一区二区性免费视频| 高清av免费在线| 欧美亚洲日本最大视频资源| 久久人妻熟女aⅴ| 黄片播放在线免费| 成人午夜精彩视频在线观看| 国产精品秋霞免费鲁丝片| 久久久久精品性色| 美女脱内裤让男人舔精品视频| 国产高清国产精品国产三级| 中文字幕亚洲精品专区| 久久久久久久久久成人| 一个人免费看片子| 日韩制服丝袜自拍偷拍| 少妇高潮的动态图| 免费观看a级毛片全部| 免费久久久久久久精品成人欧美视频 | 综合色丁香网| 在线观看www视频免费| 欧美激情 高清一区二区三区| 男女边摸边吃奶| 国产精品一区二区在线不卡| 人妻少妇偷人精品九色| 中文精品一卡2卡3卡4更新| 最近手机中文字幕大全| 97在线人人人人妻| 熟女电影av网| 亚洲一级一片aⅴ在线观看| 国产亚洲精品久久久com| 亚洲欧美色中文字幕在线| 亚洲欧美精品自产自拍| 亚洲内射少妇av| 少妇熟女欧美另类| 日本黄色日本黄色录像| 国产黄频视频在线观看| 美女xxoo啪啪120秒动态图| 女人被躁到高潮嗷嗷叫费观| 母亲3免费完整高清在线观看 | 欧美精品亚洲一区二区| 国产白丝娇喘喷水9色精品| 国产片内射在线| 乱码一卡2卡4卡精品| 蜜臀久久99精品久久宅男| 一级毛片我不卡| 亚洲久久久国产精品| 中文精品一卡2卡3卡4更新| 狠狠精品人妻久久久久久综合| 男的添女的下面高潮视频| 99久久综合免费| 蜜桃在线观看..| 18禁裸乳无遮挡动漫免费视频| 亚洲精品456在线播放app| 亚洲av欧美aⅴ国产| videossex国产| 免费在线观看完整版高清| 精品久久蜜臀av无| 内地一区二区视频在线| 99香蕉大伊视频| 51国产日韩欧美| 插逼视频在线观看| 久久韩国三级中文字幕| 国产探花极品一区二区| 国产精品一区二区在线不卡| 80岁老熟妇乱子伦牲交| 亚洲av福利一区| 国产深夜福利视频在线观看| 成人手机av| 女性被躁到高潮视频| 国产精品一区二区在线观看99| 黑人欧美特级aaaaaa片| 久久毛片免费看一区二区三区| 各种免费的搞黄视频| 国产爽快片一区二区三区| 亚洲av中文av极速乱| 久久国内精品自在自线图片| 波野结衣二区三区在线| 91精品伊人久久大香线蕉| 国产毛片在线视频| 狂野欧美激情性xxxx在线观看| 亚洲第一av免费看| 爱豆传媒免费全集在线观看| 老女人水多毛片| 美女主播在线视频| 婷婷成人精品国产| 1024视频免费在线观看| 熟女人妻精品中文字幕| 成年美女黄网站色视频大全免费| 亚洲国产最新在线播放| 在线天堂中文资源库| 久久午夜综合久久蜜桃| 日韩在线高清观看一区二区三区| 精品一区在线观看国产| 少妇的逼水好多| 日韩三级伦理在线观看| 少妇精品久久久久久久| 亚洲av男天堂| 十八禁高潮呻吟视频| 婷婷色综合大香蕉| 麻豆乱淫一区二区| 亚洲一区二区三区欧美精品| 侵犯人妻中文字幕一二三四区| 免费观看在线日韩| 久久国内精品自在自线图片| 成年美女黄网站色视频大全免费| 国产精品久久久久成人av| 伦理电影大哥的女人| 观看美女的网站| 18+在线观看网站| 日本猛色少妇xxxxx猛交久久| videos熟女内射| 成人国产麻豆网| 飞空精品影院首页| av免费观看日本| 99热这里只有是精品在线观看| 一级片免费观看大全| 夜夜骑夜夜射夜夜干| 国产男女超爽视频在线观看| 精品国产国语对白av| 女性生殖器流出的白浆| 精品少妇内射三级| 免费高清在线观看日韩| 午夜老司机福利剧场| 日韩熟女老妇一区二区性免费视频| 黄色怎么调成土黄色| 久久亚洲国产成人精品v| 免费av中文字幕在线| 美女中出高潮动态图| 最近最新中文字幕免费大全7| 自线自在国产av| 99香蕉大伊视频| 蜜桃国产av成人99| 日韩中字成人| 中文字幕制服av| 哪个播放器可以免费观看大片| 建设人人有责人人尽责人人享有的| 天天躁夜夜躁狠狠久久av| a级毛片黄视频| 日本-黄色视频高清免费观看| 亚洲婷婷狠狠爱综合网| 最近中文字幕2019免费版| 丝袜美足系列| 精品国产露脸久久av麻豆| 9191精品国产免费久久| 全区人妻精品视频| 国产精品熟女久久久久浪| 99国产综合亚洲精品| 午夜91福利影院| 亚洲av电影在线进入| 毛片一级片免费看久久久久| 亚洲成人av在线免费| 久久韩国三级中文字幕| 亚洲成人av在线免费| 好男人视频免费观看在线| 亚洲成色77777| 一二三四中文在线观看免费高清| 人人妻人人澡人人爽人人夜夜| 啦啦啦中文免费视频观看日本| 亚洲精品美女久久久久99蜜臀 | 国产熟女欧美一区二区| 一级毛片我不卡| 色网站视频免费| 成人无遮挡网站| 成人亚洲欧美一区二区av| 热re99久久国产66热| 成人亚洲精品一区在线观看| 一区二区三区乱码不卡18| 精品国产露脸久久av麻豆| 性色avwww在线观看| 人体艺术视频欧美日本| 一本大道久久a久久精品| 波野结衣二区三区在线| 夜夜骑夜夜射夜夜干| av.在线天堂| 欧美 亚洲 国产 日韩一| 你懂的网址亚洲精品在线观看| 婷婷色综合www| 你懂的网址亚洲精品在线观看| 亚洲欧美清纯卡通| 边亲边吃奶的免费视频| 亚洲欧美日韩另类电影网站| 91aial.com中文字幕在线观看| 亚洲av综合色区一区| 日本av免费视频播放| 黄色 视频免费看| 婷婷色综合大香蕉| 99久久人妻综合| 亚洲av男天堂| 男人添女人高潮全过程视频| 亚洲精品国产色婷婷电影| 亚洲欧美成人精品一区二区| 校园人妻丝袜中文字幕| 久久久久久久久久久免费av| av电影中文网址| 免费女性裸体啪啪无遮挡网站| 日本欧美国产在线视频| 伊人亚洲综合成人网| 在线免费观看不下载黄p国产| 超色免费av| 免费看av在线观看网站| 亚洲欧美成人综合另类久久久| 黑人巨大精品欧美一区二区蜜桃 | 最近的中文字幕免费完整| 日韩av不卡免费在线播放| 韩国av在线不卡| 大香蕉久久成人网| 天天躁夜夜躁狠狠躁躁| 国产精品国产三级国产av玫瑰| 国产xxxxx性猛交| 国产精品三级大全| 国产成人一区二区在线| 少妇猛男粗大的猛烈进出视频| 国产男女超爽视频在线观看| www日本在线高清视频| 国产女主播在线喷水免费视频网站| 9热在线视频观看99| 婷婷色麻豆天堂久久| 亚洲av综合色区一区| 国产极品粉嫩免费观看在线| 国产高清三级在线| www.熟女人妻精品国产 | 午夜91福利影院| 日本-黄色视频高清免费观看| 精品一区二区三区视频在线| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久久精品古装| 免费在线观看黄色视频的| 亚洲少妇的诱惑av| 一级a做视频免费观看| a级片在线免费高清观看视频| 黄色怎么调成土黄色| 成人亚洲精品一区在线观看| 免费日韩欧美在线观看| 亚洲一区二区三区欧美精品| 女人被躁到高潮嗷嗷叫费观| 色吧在线观看| 考比视频在线观看| 最近中文字幕高清免费大全6| 看十八女毛片水多多多| 在现免费观看毛片| av在线观看视频网站免费| 欧美成人精品欧美一级黄| 国产麻豆69| av在线老鸭窝| 欧美少妇被猛烈插入视频| 高清视频免费观看一区二区| 精品少妇内射三级| 午夜福利网站1000一区二区三区| 黄网站色视频无遮挡免费观看| 中文字幕人妻熟女乱码| 97在线人人人人妻| 欧美国产精品一级二级三级| 99re6热这里在线精品视频| av卡一久久| 国产成人精品久久久久久| 少妇人妻精品综合一区二区| 看非洲黑人一级黄片| 日本91视频免费播放| 2018国产大陆天天弄谢| 在线精品无人区一区二区三| 最后的刺客免费高清国语| 国产女主播在线喷水免费视频网站| 一边摸一边做爽爽视频免费| 欧美bdsm另类| 久久久久久久亚洲中文字幕| 亚洲精品国产色婷婷电影| 久久久久国产网址| 亚洲av免费高清在线观看| 婷婷色麻豆天堂久久| 午夜精品国产一区二区电影| 亚洲欧美一区二区三区黑人 | 亚洲精品久久成人aⅴ小说| 黑人欧美特级aaaaaa片| 免费av不卡在线播放| 美女国产高潮福利片在线看| 亚洲色图 男人天堂 中文字幕 | 日本色播在线视频| 少妇精品久久久久久久| 在线观看国产h片| 欧美激情极品国产一区二区三区 | 汤姆久久久久久久影院中文字幕| 国产爽快片一区二区三区| 国产一区亚洲一区在线观看| 熟女av电影| 亚洲婷婷狠狠爱综合网| 国产精品 国内视频| 欧美 日韩 精品 国产| 我的女老师完整版在线观看| kizo精华| 99九九在线精品视频| 精品人妻偷拍中文字幕| 18禁观看日本| 欧美精品国产亚洲| 国产在线视频一区二区| 热99国产精品久久久久久7| 免费高清在线观看视频在线观看| 亚洲精品一区蜜桃| 91久久精品国产一区二区三区| 久久鲁丝午夜福利片| 视频中文字幕在线观看| 亚洲精品一二三| 久久韩国三级中文字幕| 亚洲av日韩在线播放| 妹子高潮喷水视频| 桃花免费在线播放| 18禁国产床啪视频网站| 天美传媒精品一区二区| 成年人免费黄色播放视频| 在线观看免费高清a一片| 国产成人精品无人区| 亚洲人与动物交配视频| 蜜臀久久99精品久久宅男| 麻豆精品久久久久久蜜桃| 国产一区二区在线观看av| 久久久精品免费免费高清| 看免费av毛片| 国产免费一级a男人的天堂| 婷婷成人精品国产| 亚洲av福利一区| 久久精品国产亚洲av涩爱| 高清欧美精品videossex| 大香蕉久久成人网| 免费黄频网站在线观看国产| 黑人欧美特级aaaaaa片| 肉色欧美久久久久久久蜜桃| 另类精品久久| 国产欧美另类精品又又久久亚洲欧美| 深夜精品福利| 日本欧美国产在线视频| 日本色播在线视频| 日韩电影二区| 日日摸夜夜添夜夜爱| 亚洲精品aⅴ在线观看| 久久久久久人人人人人| av又黄又爽大尺度在线免费看| 伦精品一区二区三区| 国产精品久久久久成人av| 国产精品无大码| 丝袜喷水一区| 免费大片18禁| 日本黄大片高清| 国产精品久久久久久久久免| 国产成人精品婷婷| 欧美激情极品国产一区二区三区 | 免费观看a级毛片全部| 久久国产精品男人的天堂亚洲 | 美女脱内裤让男人舔精品视频| 欧美成人午夜精品| 亚洲av男天堂| 亚洲,欧美,日韩| 热99国产精品久久久久久7| 日本色播在线视频| 久久婷婷青草| 好男人视频免费观看在线| 2022亚洲国产成人精品| 久久久国产欧美日韩av| 国产日韩一区二区三区精品不卡| 99国产综合亚洲精品| 久久99热6这里只有精品| 我的女老师完整版在线观看| 欧美+日韩+精品| av不卡在线播放| 欧美日韩精品成人综合77777| 亚洲国产看品久久| 国产精品久久久久久av不卡| 国产黄色免费在线视频| 丝袜喷水一区| av国产精品久久久久影院| 精品福利永久在线观看| 天天躁夜夜躁狠狠躁躁| 久久ye,这里只有精品| 国产有黄有色有爽视频| 国产一级毛片在线| 一级a做视频免费观看| 亚洲av电影在线进入| 天堂俺去俺来也www色官网| 亚洲精品乱久久久久久| 亚洲精品美女久久av网站| 欧美日韩av久久| 极品人妻少妇av视频| 免费在线观看完整版高清| 免费日韩欧美在线观看| 亚洲欧洲国产日韩| 宅男免费午夜| 色婷婷久久久亚洲欧美| 日韩不卡一区二区三区视频在线| 亚洲欧美成人精品一区二区| av不卡在线播放| 最新中文字幕久久久久| 亚洲欧美成人精品一区二区| 国内精品宾馆在线| 日韩中字成人| 国产精品人妻久久久久久| 99久久精品国产国产毛片| 精品国产乱码久久久久久小说| 国产av一区二区精品久久| 午夜福利影视在线免费观看| 69精品国产乱码久久久| 久久毛片免费看一区二区三区| 又黄又粗又硬又大视频| 日韩伦理黄色片| 国产不卡av网站在线观看| 久久国内精品自在自线图片| 99热国产这里只有精品6| 18禁国产床啪视频网站| 老司机影院成人| 精品久久久久久电影网| 午夜福利,免费看| 下体分泌物呈黄色| 国产成人免费观看mmmm| 日韩,欧美,国产一区二区三区| 少妇人妻精品综合一区二区| 日韩中文字幕视频在线看片| 亚洲中文av在线| a级毛色黄片| 一边亲一边摸免费视频| 亚洲av电影在线进入| 美女视频免费永久观看网站|