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

    基于似然比檢驗(yàn)的橡膠密封件加速試驗(yàn)機(jī)理一致性判別方法*

    2019-04-26 05:33:12汪亞順譚源源
    關(guān)鍵詞:估計(jì)值對數(shù)機(jī)理

    魯 相,汪亞順,陳 循,譚源源

    (1. 國防科技大學(xué) 智能科學(xué)學(xué)院, 湖南 長沙 410073; 2. 中國空氣動力研究與發(fā)展中心 設(shè)備設(shè)計(jì)及測試技術(shù)研究所, 四川 綿陽 621000)

    隨著產(chǎn)品質(zhì)量水平的不斷提高,近代研制出的產(chǎn)品普遍具有高可靠、長壽命的特點(diǎn),例如發(fā)光二極管及橡膠密封件等元器件中位壽命可以超過10年[1]。對于這類高可靠性產(chǎn)品,現(xiàn)代工業(yè)廣泛利用加速退化試驗(yàn)預(yù)測其使用壽命和貯存壽命。在加速退化試驗(yàn)中,將產(chǎn)品置于加速應(yīng)力水平下并記錄產(chǎn)品樣本的性能退化數(shù)據(jù),利用加速模型建立應(yīng)力與退化速率之間的關(guān)系。當(dāng)各試驗(yàn)應(yīng)力水平下的失效機(jī)理相同時,就可以將加速應(yīng)力下的試驗(yàn)結(jié)果外推,進(jìn)而估計(jì)在正常應(yīng)力狀況下的產(chǎn)品壽命。各試驗(yàn)應(yīng)力水平下失效機(jī)理具有一致性是加速退化試驗(yàn)結(jié)論具備有效性的首要前提,因此需要研究加速退化試驗(yàn)機(jī)理一致性判別方法。

    利用試驗(yàn)數(shù)據(jù)統(tǒng)計(jì)分析方法或物理分析測試手段,判斷失效機(jī)理在試驗(yàn)各加速應(yīng)力水平下是否發(fā)生變化,這個過程被定義為失效機(jī)理的一致性判別。目前失效機(jī)理一致性判別方法可分為3類:概率圖檢驗(yàn)法[2-4]、基于模型參數(shù)一致性檢驗(yàn)法[5-9]、基于灰色預(yù)測模型檢驗(yàn)法[10-11]。這些一致性判別方法各有優(yōu)劣,但存在共同的缺點(diǎn):通常只對壽命服從某一特定分布的產(chǎn)品有效,應(yīng)用場合受限,因而難以推廣到其他統(tǒng)計(jì)分布,并且這些方法對判別結(jié)果的誤判風(fēng)險(xiǎn)沒有開展定量分析。因此,提出基于似然比檢驗(yàn)的失效機(jī)理一致性判別方法,該方法適用于多個統(tǒng)計(jì)分布,可以定量分析誤判風(fēng)險(xiǎn),并且模型參數(shù)估計(jì)精度較高。

    1 機(jī)理一致性判別的數(shù)學(xué)模型及參數(shù)估計(jì)

    1.1 兩類加速模型

    通常利用混合效應(yīng)模型[12-13]對高可靠產(chǎn)品退化建模,并認(rèn)為性能退化量初始值B及退化軌跡形狀參數(shù)α是固定不變的,不同個體間的退化速率存在差異。產(chǎn)品的退化速率Kij在各應(yīng)力水平Si下一般服從對數(shù)正態(tài)分布LN(μi,σ2),其中σ為對數(shù)標(biāo)準(zhǔn)差,μi為對數(shù)均值(與應(yīng)力水平有關(guān))。選擇如式(1)所示的混合效應(yīng)模型。

    (1)

    在混合效應(yīng)模型下,當(dāng)失效機(jī)理不變時,提高載荷或溫度的試驗(yàn)產(chǎn)品統(tǒng)一加速模型為:

    μi=a+b·φ(Si)

    (2)

    式中,φ(Si)為轉(zhuǎn)換應(yīng)力水平,溫度為加速應(yīng)力時φ(Si)=1/Si,載荷為加速應(yīng)力時φ(Si)=lnSi。

    當(dāng)加速應(yīng)力范圍較大導(dǎo)致失效機(jī)理變化時,式(2)中的參數(shù)a和b不再為常數(shù),式(2)在較高應(yīng)力水平下無效。比式(2)更適用的模型可以表述為:

    μi=ai+bi·φ(Si)

    (3)

    可得到如下原假設(shè)H0:μi=a+b·φ(Si)和備擇假設(shè)H1:μi=ai+bi·φ(Si)。

    從上述定義可以看出,在原假設(shè)中假定退化速率-應(yīng)力關(guān)系為對數(shù)線性,將原假設(shè)對應(yīng)模型稱為對數(shù)線性加速模型。相比之下,在備擇假設(shè)中模型參數(shù)與應(yīng)力水平有關(guān),將備擇假設(shè)對應(yīng)模型稱為非對數(shù)線性加速模型。式(2)與式(3)的表達(dá)式表明原假設(shè)對應(yīng)模型是備擇假設(shè)對應(yīng)模型的一種特例。在對數(shù)線性加速模型中參數(shù)a和b被限定為常數(shù),導(dǎo)致其獨(dú)立參數(shù)少于非對數(shù)線性加速模型。如果兩種假設(shè)模型中一種模型是另一模型的特例,則可用似然比檢驗(yàn)比較兩種模型的擬合優(yōu)度[14-15],因此似然比檢驗(yàn)可用于檢驗(yàn)上述原假設(shè)。

    1.2 似然比檢驗(yàn)的基本準(zhǔn)則

    假定有來自于分布T的隨機(jī)樣本t1,t2,…,tn,令θ表示該分布未知參數(shù)向量,則該樣本的似然函數(shù)為:

    (4)

    令L0(θ)為對數(shù)線性加速模型參數(shù)的似然函數(shù),L1(θ)為非對數(shù)線性加速模型參數(shù)的似然函數(shù)。對數(shù)線性加速模型參數(shù)的極大似然估計(jì)(Maximum Likelihood Estimation, MLE)為:

    (5)

    非對數(shù)線性加速模型參數(shù)的極大似然估計(jì)為:

    (6)

    似然比統(tǒng)計(jì)量可表示為:

    (7)

    由于包含更多獨(dú)立參數(shù),通常非對數(shù)線性模型對樣本觀測值的擬合效果比對數(shù)線性模型好,但兩類模型也可能都有較好的擬合效果。非對數(shù)線性模型參數(shù)的似然函數(shù)值不小于對數(shù)線性模型的似然函數(shù)值,似然比統(tǒng)計(jì)量滿足:

    (8)

    如果對數(shù)線性模型與非對數(shù)線性模型一樣適用,則似然比較大,反之較小。當(dāng)似然比統(tǒng)計(jì)量數(shù)值過小時似然比檢驗(yàn)法會拒絕原假設(shè)(對數(shù)線性加速模型),似然比檢驗(yàn)法的臨界域或拒絕域?yàn)椋?/p>

    W={λ≤c} 0

    (9)

    式中,臨界值c與參數(shù)λ、β有關(guān)。β是似然比檢驗(yàn)中規(guī)定的顯著性水平,β表示將失效機(jī)理相同誤判為失效機(jī)理變化的風(fēng)險(xiǎn)大小。

    1.3 數(shù)學(xué)模型及參數(shù)估計(jì)

    根據(jù)上述分析,首次應(yīng)用似然比檢驗(yàn)法辨識加速退化試驗(yàn)中加速模型的變化,并計(jì)算兩種假設(shè)模型下對數(shù)似然函數(shù)的極大值。設(shè)恒定應(yīng)力加速退化試驗(yàn)(Constant Stress Accelerated Degradation Testing, CSADT)采用q個應(yīng)力水平,在每個加速應(yīng)力水平Si(i=1,2,…,q)下測試ni個試驗(yàn)樣本的性能退化量,測試mi次,測試的時間節(jié)點(diǎn)為tik(k=1,2,…,mi),在每個應(yīng)力水平下的試驗(yàn)截止時間為ti,測得退化量數(shù)據(jù)為yijk(j=1,2,…,ni)。建模時取樣本的性能退化量對數(shù)值lnyijk統(tǒng)計(jì)分析,混合效應(yīng)模型下其條件累積分布及條件概率密度函數(shù)可以表示為:

    (10)

    (11)

    其中,Φ(·)和φ(·)表示標(biāo)準(zhǔn)正態(tài)分布的累積分布函數(shù)及概率密度函數(shù)。由式(10)及式(11)可得恒定應(yīng)力加速退化試驗(yàn)數(shù)據(jù)的似然函數(shù)為:

    (12)

    式中,f(Kij)代表退化速率Kij的概率密度函數(shù)。

    (13)

    假定采用溫度作為加速應(yīng)力,在對數(shù)線性加速模型中ai=a,bi=b,取向量δ=(lnB,α,a,b),Kij=exp(a+b/Si+cij),yijk=(yij1,…,yijmi)T,tik=(ti1,…,timi)T,協(xié)方差矩陣D=σ2,cij服從正態(tài)分布N(0,σ2)。采用近似極大似然估計(jì)法求解模型參數(shù)。

    步驟1:在固定協(xié)方差矩陣的情況下,用非線性最小二乘法獲得加速模型參數(shù)a、b、cij和形狀參數(shù)α及退化量初始值B的估計(jì)值,估計(jì)值是使式(14)最小化的解。

    (14)

    步驟2:首先將模型函數(shù)在第一步獲得的各參數(shù)估計(jì)值處作Taylor展開,根據(jù)模型函數(shù)表達(dá)式,有

    其中, *代表哈達(dá)馬積。

    然后利用第一步獲得的估計(jì)值求取協(xié)方差矩陣D及其他參數(shù)的極大似然函數(shù)表達(dá)式,取

    可得對數(shù)線性加速模型的對數(shù)似然函數(shù)為:

    (15)

    設(shè)定好數(shù)值迭代算法的初始解,利用Newton-Raphson算法[16-19]求解未知參數(shù)估計(jì)值,按式(15)計(jì)算對數(shù)似然函數(shù)的極大值。

    在非對數(shù)線性加速模型中,取向量κ=(lnB,α,μ1,…,μq),Kij=exp(μi+cij),yijk=(yij1,…,yijmi)T,tik=(ti1,…,timi)T,協(xié)方差矩陣D=[σ2]。

    按以下步驟進(jìn)行參數(shù)估計(jì)。

    步驟1:在固定協(xié)方差矩陣的情況下,用非線性最小二乘法獲得加速模型參數(shù)μi、cij和形狀參數(shù)α及退化量初始值B的估計(jì)值,估計(jì)值是使式(16)最小化的解。

    (16)

    步驟2:首先將模型函數(shù)在第一步獲得的各參數(shù)估計(jì)值處作Taylor展開,根據(jù)模型函數(shù)表達(dá)式,有

    同樣求取協(xié)方差矩陣D及其他參數(shù)的極大似然函數(shù)表達(dá)式,取

    可得非對數(shù)線性加速模型的對數(shù)似然函數(shù)為:

    (17)

    利用Newton-Raphson算法類似得到對數(shù)似然函數(shù)的極大值。

    2 失效機(jī)理一致性判別準(zhǔn)則及方法流程

    算出兩種假設(shè)下似然函數(shù)極大值后,需要確定似然比統(tǒng)計(jì)量的分布,以構(gòu)建判別準(zhǔn)則。多數(shù)情況下難以確定似然比統(tǒng)計(jì)量λ的精確抽樣分布,但可以確定大樣本場合下對數(shù)似然比統(tǒng)計(jì)量的漸進(jìn)分布,令

    (18)

    根據(jù)Wilks提出的廣義似然比統(tǒng)計(jì)量的極限分布定理[20],在原假設(shè)成立的情況下,統(tǒng)計(jì)量Λ漸進(jìn)服從自由度為v的卡方分布χ2(v)。自由度v等于在假設(shè)H0和H1下的獨(dú)立參數(shù)數(shù)目之差,有v=q-2。由式(9)可知,當(dāng)λ≤c時原假設(shè)H0被拒絕,因此當(dāng)Λ≥-2lnc時對數(shù)線性加速模型無效。利用顯著性水平β及統(tǒng)計(jì)量Λ的分布,可以得出

    (19)

    式中,等式右端表示卡方分布χ2(q-2)的1-β分位數(shù)。

    利用式(15)和式(17)計(jì)算恒定應(yīng)力加速退化試驗(yàn)中兩類模型的極大對數(shù)似然函數(shù)的差值,得到:

    ΛCSM=-2[lnL0(δ,D,σε)-lnL1(κ,D,σε)]

    (20)

    失效機(jī)理一致性的判別準(zhǔn)則是:

    1)如果0<ΛCSM<-2lnc,失效機(jī)理在所有應(yīng)力水平下相同;

    2)如果ΛCSM≥-2lnc,失效機(jī)理在某一高應(yīng)力水平下變化。

    根據(jù)上述分析,基于似然比檢驗(yàn)的失效機(jī)理一致性判別方法流程如圖1所示。

    3 仿真算例

    在發(fā)光二極管(Light Emitting Diode, LED)加速試驗(yàn)中,最高溫度下激活能可能發(fā)生變化。本例中設(shè)定退化失效機(jī)理發(fā)生變化,在378 K激活能從16.1 kJ/mol轉(zhuǎn)變?yōu)?2.2 kJ/mol。

    圖1 失效機(jī)理一致性判別方法流程Fig.1 Flow chart of identification method of failure mechanism consistency

    仿真試驗(yàn)溫度應(yīng)力水平為298 K、338 K、378 K,試驗(yàn)樣本量為6,測量間隔為336 h,每個水平下測量次數(shù)為11。設(shè)定退化模型參數(shù)分別為a1=-0.546 6,b1=-1933.25(298~338 K),a2=5.173 1,b2=-3866.49(338~378 K),B=1.051 6,α=0.65,σ=0.185,σε=0.014。利用式(1)生成仿真退化軌跡。

    利用近似極大似然估計(jì)法分析仿真數(shù)據(jù),得出對數(shù)線性模型中未知參數(shù)的極大似然估計(jì),將極大似然估計(jì)值與未知參數(shù)真值進(jìn)行比較,結(jié)果如表1所示。

    表1結(jié)果表明,采用指定的加速模型進(jìn)行建模及統(tǒng)計(jì)分析時,加速模型參數(shù)估計(jì)值a、b與真值的偏差較大,估計(jì)值近似于兩真值的均值,導(dǎo)致在298 K以下的預(yù)測壽命偏大。

    表1 對數(shù)線性模型參數(shù)估計(jì)值

    非對數(shù)線性模型中未知參數(shù)極大似然估計(jì)值與真值的比較如表2所示,可以看出,模型各參數(shù)的估計(jì)值與真值相差較小,非對數(shù)線性模型擬合效果明顯優(yōu)于對數(shù)線性模型。

    進(jìn)一步可以算得溫度水平在298~378 K下檢驗(yàn)統(tǒng)計(jì)量為:

    ΛCSM=-2lnλ=-2(lnL0-lnL1)

    =-2×(510.254 7-518.963 2)=17.417

    (21)

    在298~378 K溫度范圍內(nèi)的試驗(yàn)應(yīng)力水平數(shù)為3,給定顯著性水平β=5%,臨界值滿足-2lnc=3.84<ΛCSM。因此在298~378 K范圍內(nèi)失效機(jī)理發(fā)生變化,這一判別結(jié)果與失效機(jī)理發(fā)生變化的原始設(shè)定相符,驗(yàn)證了方法的正確性。

    表2 非對數(shù)線性模型參數(shù)估計(jì)值

    4 應(yīng)用實(shí)例

    丁腈橡膠密封件是典型的高可靠性產(chǎn)品,但易于發(fā)生熱氧老化。受熱氧老化的影響,丁腈橡膠密封件在貯存或使用期密封性能會不斷退化。隨著密封性能的退化,橡膠密封性能的高彈性復(fù)原力不斷減小,壓縮永久變形不斷增大。因此可以通過觀測壓縮永久變形cs的變化來監(jiān)測退化趨勢,當(dāng)壓縮永久變形上升到某一臨界值時,認(rèn)為橡膠密封件密封失效。試驗(yàn)采用安裝于系統(tǒng)中的橡膠密封件的壓縮率,樣本可以代表在空氣介質(zhì)存儲下系統(tǒng)內(nèi)受壓丁腈橡膠密封件的退化。試驗(yàn)總樣本量為20,選擇在5個溫度應(yīng)力水平(333 K、343 K、353 K、373 K、393 K)下進(jìn)行加速退化試驗(yàn),每個應(yīng)力水平下樣本量為4。密封件在333 K及393 K應(yīng)力水平下測得的壓縮永久變形量如圖2所示,其他應(yīng)力水平下曲線形狀類似,不再給出。根據(jù)圖1所示流程對數(shù)據(jù)進(jìn)行失效機(jī)理一致性判別。建立退化模型,將測得的壓縮永久變形cs(無量綱)轉(zhuǎn)換為性能退化數(shù)據(jù)y,轉(zhuǎn)換關(guān)系式為y=1-cs,采用混合效應(yīng)模型(1)對試樣退化規(guī)律進(jìn)行描述。

    (a) 333 K丁腈密封圈加速退化試驗(yàn)數(shù)據(jù) (a) CSADT data of nitrile seal rings at 333 K

    (b) 393 K丁腈密封圈加速退化試驗(yàn)數(shù)據(jù) (b) CSADT data of nitrile seal rings at 393 K圖2 丁腈密封圈恒定應(yīng)力加速退化試驗(yàn)數(shù)據(jù)Fig.2 CSADT data of nitrile seal rings

    利用近似極大似然估計(jì)方法,令ni=4(i=1,2,3),可求解333~353 K、333~373 K及333~393 K下模型的未知參數(shù),設(shè)定顯著性水平β=10%,進(jìn)一步算得密封圈試樣在三個不同溫度范圍的假設(shè)檢驗(yàn)結(jié)果,如表3所示。

    表3 密封圈失效機(jī)理一致性判別

    根據(jù)失效機(jī)理一致性判別方法,結(jié)合表3計(jì)算結(jié)果可知,試驗(yàn)所用丁腈橡膠密封圈在373 K以上時溫度失效機(jī)理發(fā)生變化。Wise等[21]、Le等[22]對類似產(chǎn)品開展失效物理分析,得到368 K溫度以上丁腈橡膠失效機(jī)理變化的結(jié)果,與本文所得結(jié)論基本一致,驗(yàn)證了本文方法的有效性。

    5 結(jié)論

    高可靠性產(chǎn)品在進(jìn)行加速退化試驗(yàn)時,失效機(jī)理在高應(yīng)力水平下可能發(fā)生變化,導(dǎo)致錯誤的壽命預(yù)測結(jié)果?;谒迫槐葯z驗(yàn)的加速退化試驗(yàn)機(jī)理一致性判別方法,適用范圍廣,能夠定量分析誤判風(fēng)險(xiǎn),并給出機(jī)理一致性判別結(jié)論,得到準(zhǔn)確的壽命預(yù)測結(jié)果。建立非線性加速模型與高可靠長壽命產(chǎn)品失效機(jī)理變化間的聯(lián)系,并將似然比檢驗(yàn)應(yīng)用于加速退化試驗(yàn)中對加速模型的辨識,建立了有效的失效機(jī)理一致性判別準(zhǔn)則及方法。

    通過仿真案例及應(yīng)用實(shí)例,進(jìn)一步驗(yàn)證了所提方法的有效性。多應(yīng)力加速退化試驗(yàn)下產(chǎn)品失效機(jī)理的一致性判別較為復(fù)雜,相關(guān)方法還需要進(jìn)一步研究。

    猜你喜歡
    估計(jì)值對數(shù)機(jī)理
    含有對數(shù)非線性項(xiàng)Kirchhoff方程多解的存在性
    指數(shù)與對數(shù)
    隔熱纖維材料的隔熱機(jī)理及其應(yīng)用
    指數(shù)與對數(shù)
    一道樣本的數(shù)字特征與頻率分布直方圖的交匯問題
    煤層氣吸附-解吸機(jī)理再認(rèn)識
    中國煤層氣(2019年2期)2019-08-27 00:59:30
    對數(shù)簡史
    統(tǒng)計(jì)信息
    2018年4月世界粗鋼產(chǎn)量表(續(xù))萬噸
    霧霾機(jī)理之問
    欧美极品一区二区三区四区| 在线播放无遮挡| 免费高清在线观看视频在线观看| 免费观看性生交大片5| 国产精品国产av在线观看| 国产亚洲欧美精品永久| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久国产精品人妻一区二区| 亚洲一级一片aⅴ在线观看| 国产成人精品福利久久| 不卡视频在线观看欧美| 久久6这里有精品| 亚洲国产日韩一区二区| 久久韩国三级中文字幕| .国产精品久久| 一区二区三区乱码不卡18| 高清午夜精品一区二区三区| 亚洲成人手机| 九九久久精品国产亚洲av麻豆| 高清黄色对白视频在线免费看 | 美女xxoo啪啪120秒动态图| 插阴视频在线观看视频| 亚洲欧美日韩无卡精品| 国产久久久一区二区三区| 日本av手机在线免费观看| av在线观看视频网站免费| av专区在线播放| 国产精品一及| 国产日韩欧美在线精品| 少妇的逼好多水| 亚洲精华国产精华液的使用体验| 中国国产av一级| 欧美人与善性xxx| 亚洲色图综合在线观看| 人妻制服诱惑在线中文字幕| 欧美国产精品一级二级三级 | 精品午夜福利在线看| 美女中出高潮动态图| tube8黄色片| 国产精品蜜桃在线观看| 五月天丁香电影| 在线免费十八禁| 汤姆久久久久久久影院中文字幕| 99热全是精品| 街头女战士在线观看网站| 91精品国产九色| 亚洲国产av新网站| 黄色欧美视频在线观看| 国产亚洲午夜精品一区二区久久| 午夜免费男女啪啪视频观看| 国产精品一区二区在线不卡| 欧美激情国产日韩精品一区| 国产av国产精品国产| 最近中文字幕2019免费版| 午夜老司机福利剧场| 国产免费一级a男人的天堂| 久久久久国产精品人妻一区二区| 久久99蜜桃精品久久| 中文字幕制服av| 久久久久久久久大av| 日韩视频在线欧美| 丰满迷人的少妇在线观看| 国产欧美亚洲国产| 18禁裸乳无遮挡免费网站照片| av在线蜜桃| a级毛色黄片| 亚洲av成人精品一二三区| 国产成人精品一,二区| 伦理电影免费视频| 久久久久国产网址| 男人和女人高潮做爰伦理| 亚洲精品自拍成人| 国产精品蜜桃在线观看| 日韩三级伦理在线观看| 午夜免费男女啪啪视频观看| 综合色丁香网| 国产熟女欧美一区二区| 综合色丁香网| 黄色一级大片看看| 我要看黄色一级片免费的| 婷婷色麻豆天堂久久| 蜜桃久久精品国产亚洲av| 中国三级夫妇交换| 国产精品不卡视频一区二区| 欧美xxxx黑人xx丫x性爽| 男女啪啪激烈高潮av片| 午夜精品国产一区二区电影| 啦啦啦在线观看免费高清www| 国产伦在线观看视频一区| 午夜老司机福利剧场| 自拍偷自拍亚洲精品老妇| 亚洲国产精品国产精品| 国产精品爽爽va在线观看网站| 三级经典国产精品| 51国产日韩欧美| 欧美zozozo另类| 日韩免费高清中文字幕av| 日本av手机在线免费观看| 久久97久久精品| 国产免费一区二区三区四区乱码| 中文乱码字字幕精品一区二区三区| 老司机影院成人| 亚洲,一卡二卡三卡| 少妇的逼好多水| 中国三级夫妇交换| 欧美精品人与动牲交sv欧美| 精品亚洲成a人片在线观看 | 尤物成人国产欧美一区二区三区| 久久久久人妻精品一区果冻| 亚洲第一区二区三区不卡| 国产精品久久久久久久久免| 一级毛片电影观看| 熟女电影av网| 极品教师在线视频| 中国国产av一级| 日韩av免费高清视频| 丰满迷人的少妇在线观看| 麻豆国产97在线/欧美| 久久99热这里只有精品18| 国产黄色免费在线视频| 日本黄大片高清| 欧美精品人与动牲交sv欧美| 中文在线观看免费www的网站| 美女福利国产在线 | 国产精品国产三级国产av玫瑰| 久久鲁丝午夜福利片| 热re99久久精品国产66热6| 欧美少妇被猛烈插入视频| 大话2 男鬼变身卡| 我的老师免费观看完整版| 插逼视频在线观看| 99热这里只有是精品在线观看| 2018国产大陆天天弄谢| 精品国产一区二区三区久久久樱花 | av国产久精品久网站免费入址| 日本av免费视频播放| 高清毛片免费看| 综合色丁香网| 自拍偷自拍亚洲精品老妇| 少妇的逼水好多| 五月天丁香电影| 一级毛片电影观看| 日韩欧美精品免费久久| 亚洲人与动物交配视频| 最近最新中文字幕免费大全7| 大香蕉久久网| 国产色婷婷99| 内射极品少妇av片p| 成年免费大片在线观看| 欧美另类一区| 久久久亚洲精品成人影院| 日本欧美视频一区| 亚洲四区av| 黑人高潮一二区| 日韩精品有码人妻一区| 麻豆成人av视频| 人妻 亚洲 视频| 国产熟女欧美一区二区| 99久久精品热视频| 亚洲精品国产色婷婷电影| 一区二区三区免费毛片| 草草在线视频免费看| 看免费成人av毛片| 一级av片app| 成人18禁高潮啪啪吃奶动态图 | 大话2 男鬼变身卡| 黄色欧美视频在线观看| 在线 av 中文字幕| 国产男人的电影天堂91| 人妻夜夜爽99麻豆av| 精品久久久噜噜| 欧美精品亚洲一区二区| 国产国拍精品亚洲av在线观看| 亚洲精品aⅴ在线观看| 国产伦精品一区二区三区四那| 国产精品99久久99久久久不卡 | 国产精品一区二区在线观看99| 国产黄频视频在线观看| 国产精品福利在线免费观看| 免费久久久久久久精品成人欧美视频 | 色婷婷av一区二区三区视频| 国产亚洲午夜精品一区二区久久| 一区二区三区乱码不卡18| 一二三四中文在线观看免费高清| 黄片wwwwww| h日本视频在线播放| 亚洲图色成人| 欧美xxxx黑人xx丫x性爽| 日本wwww免费看| 美女脱内裤让男人舔精品视频| 国产成人a区在线观看| 国产老妇伦熟女老妇高清| 蜜桃亚洲精品一区二区三区| 蜜桃在线观看..| 免费看av在线观看网站| 91在线精品国自产拍蜜月| 日韩伦理黄色片| 国产男女超爽视频在线观看| 在线观看一区二区三区激情| 五月天丁香电影| 狂野欧美激情性xxxx在线观看| 自拍偷自拍亚洲精品老妇| 色综合色国产| 国产 一区 欧美 日韩| 日韩亚洲欧美综合| 搡老乐熟女国产| 26uuu在线亚洲综合色| 你懂的网址亚洲精品在线观看| 秋霞伦理黄片| 国产亚洲一区二区精品| 91精品伊人久久大香线蕉| 中国三级夫妇交换| 亚洲精品乱码久久久久久按摩| 免费看光身美女| 久久久a久久爽久久v久久| 午夜激情福利司机影院| 色视频在线一区二区三区| 亚洲国产精品专区欧美| 国产精品欧美亚洲77777| 午夜福利影视在线免费观看| 国产精品无大码| 自拍偷自拍亚洲精品老妇| 少妇丰满av| 男女边摸边吃奶| 亚洲经典国产精华液单| 国产欧美日韩精品一区二区| 国产精品伦人一区二区| 99久久人妻综合| 黄色日韩在线| 草草在线视频免费看| 亚洲av中文av极速乱| 国产精品.久久久| 欧美高清性xxxxhd video| 天美传媒精品一区二区| 在线天堂最新版资源| 国产极品天堂在线| www.av在线官网国产| 欧美区成人在线视频| 最新中文字幕久久久久| 久久ye,这里只有精品| av卡一久久| 成年人午夜在线观看视频| 啦啦啦啦在线视频资源| 成年av动漫网址| 黄色欧美视频在线观看| 日韩国内少妇激情av| 天天躁夜夜躁狠狠久久av| 国产精品一区二区在线不卡| 欧美精品国产亚洲| 国产一区二区在线观看日韩| 男的添女的下面高潮视频| 黄色日韩在线| 亚洲国产精品999| 国产 一区 欧美 日韩| 蜜桃亚洲精品一区二区三区| 亚洲精品国产av成人精品| 嫩草影院新地址| 国产亚洲欧美精品永久| 国内揄拍国产精品人妻在线| 一级毛片电影观看| av女优亚洲男人天堂| 国产乱人偷精品视频| 人妻系列 视频| 国产老妇伦熟女老妇高清| 欧美成人a在线观看| 久久影院123| 精品午夜福利在线看| 久久久久精品久久久久真实原创| 天天躁夜夜躁狠狠久久av| 熟女电影av网| 少妇高潮的动态图| 91狼人影院| 日韩视频在线欧美| av在线app专区| 国产一级毛片在线| 国产 精品1| 插阴视频在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 欧美国产精品一级二级三级 | 天堂中文最新版在线下载| h日本视频在线播放| 亚洲久久久国产精品| 18+在线观看网站| 在线观看免费日韩欧美大片 | 亚洲精品久久午夜乱码| 美女主播在线视频| 亚洲欧美一区二区三区黑人 | 久久久午夜欧美精品| 亚洲欧美一区二区三区黑人 | 久久国产乱子免费精品| 成人毛片a级毛片在线播放| 99精国产麻豆久久婷婷| 日本黄色日本黄色录像| 国产在线男女| 免费人成在线观看视频色| videos熟女内射| 熟妇人妻不卡中文字幕| 国产精品蜜桃在线观看| 国产精品国产三级国产专区5o| 国产成人aa在线观看| 热99国产精品久久久久久7| 国产伦精品一区二区三区四那| 2022亚洲国产成人精品| 在线观看美女被高潮喷水网站| 日本黄色片子视频| 另类亚洲欧美激情| 精品国产乱码久久久久久小说| 男女下面进入的视频免费午夜| 久久99蜜桃精品久久| 久久99精品国语久久久| 久久久午夜欧美精品| 亚洲精品日韩在线中文字幕| 边亲边吃奶的免费视频| 日本黄大片高清| 大话2 男鬼变身卡| 精品熟女少妇av免费看| 精品熟女少妇av免费看| 欧美三级亚洲精品| 亚洲欧美精品自产自拍| av.在线天堂| 成年美女黄网站色视频大全免费 | 欧美一区二区亚洲| 国产成人a∨麻豆精品| 自拍欧美九色日韩亚洲蝌蚪91 | 18禁动态无遮挡网站| 免费看日本二区| 女的被弄到高潮叫床怎么办| 久久精品国产亚洲av涩爱| 欧美精品一区二区免费开放| 亚洲美女视频黄频| 久久99热这里只有精品18| 国产成人91sexporn| 亚洲精品乱久久久久久| 婷婷色综合大香蕉| 免费播放大片免费观看视频在线观看| 久久精品久久久久久噜噜老黄| 看非洲黑人一级黄片| 少妇丰满av| 视频区图区小说| 最近最新中文字幕免费大全7| 国产欧美另类精品又又久久亚洲欧美| 中国三级夫妇交换| 大片电影免费在线观看免费| 国产精品免费大片| 在线观看国产h片| 免费观看性生交大片5| 天天躁日日操中文字幕| 美女福利国产在线 | 最近2019中文字幕mv第一页| 久久久久性生活片| 欧美3d第一页| 高清毛片免费看| 亚洲av成人精品一区久久| 99久久精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲成色77777| av一本久久久久| 99热这里只有精品一区| 亚洲欧美日韩东京热| 亚洲美女黄色视频免费看| 有码 亚洲区| 亚洲国产色片| 成人美女网站在线观看视频| 欧美精品人与动牲交sv欧美| 免费观看a级毛片全部| 久久久久视频综合| 一级毛片电影观看| 国产成人免费无遮挡视频| 亚洲欧美清纯卡通| 最近的中文字幕免费完整| 欧美成人精品欧美一级黄| 久久久久久久亚洲中文字幕| 欧美变态另类bdsm刘玥| 久久精品国产亚洲av涩爱| 插逼视频在线观看| freevideosex欧美| 亚洲怡红院男人天堂| 爱豆传媒免费全集在线观看| 丰满乱子伦码专区| 国产精品伦人一区二区| 免费av不卡在线播放| 亚洲国产精品专区欧美| 18禁裸乳无遮挡免费网站照片| 亚洲av在线观看美女高潮| 久久久久久久久久人人人人人人| 99热网站在线观看| 午夜老司机福利剧场| 在线观看一区二区三区激情| 一个人看视频在线观看www免费| 天堂中文最新版在线下载| 91aial.com中文字幕在线观看| 亚洲国产精品专区欧美| 激情 狠狠 欧美| 欧美xxxx黑人xx丫x性爽| 在线免费观看不下载黄p国产| 免费黄网站久久成人精品| 日韩电影二区| 成人免费观看视频高清| 成人毛片a级毛片在线播放| 国产一区二区在线观看日韩| 又粗又硬又长又爽又黄的视频| av在线播放精品| 国产 精品1| 亚洲图色成人| 国产伦理片在线播放av一区| 高清午夜精品一区二区三区| 婷婷色麻豆天堂久久| 日韩大片免费观看网站| 观看美女的网站| av女优亚洲男人天堂| 蜜桃久久精品国产亚洲av| 国产成人午夜福利电影在线观看| 人体艺术视频欧美日本| 亚洲真实伦在线观看| 我的女老师完整版在线观看| 国产综合精华液| 下体分泌物呈黄色| 久久久久久人妻| 国产伦理片在线播放av一区| 欧美xxⅹ黑人| 免费播放大片免费观看视频在线观看| 欧美国产精品一级二级三级 | 国产午夜精品久久久久久一区二区三区| 各种免费的搞黄视频| 啦啦啦视频在线资源免费观看| 最近的中文字幕免费完整| 国产 一区 欧美 日韩| 看免费成人av毛片| 久久久午夜欧美精品| 欧美精品国产亚洲| 欧美成人一区二区免费高清观看| 美女高潮的动态| 联通29元200g的流量卡| 夫妻性生交免费视频一级片| 国产一区二区三区av在线| 国国产精品蜜臀av免费| 成人无遮挡网站| 亚洲欧美一区二区三区国产| 内地一区二区视频在线| 人妻一区二区av| 纵有疾风起免费观看全集完整版| 亚洲伊人久久精品综合| 国产高清有码在线观看视频| 欧美xxxx性猛交bbbb| 国产日韩欧美在线精品| 久久国产精品大桥未久av | 日韩国内少妇激情av| 这个男人来自地球电影免费观看 | 久久久久国产网址| 国产黄色免费在线视频| 国产男女内射视频| 亚洲综合精品二区| 国产久久久一区二区三区| 麻豆成人av视频| 久久久午夜欧美精品| 国产精品欧美亚洲77777| 亚洲成人中文字幕在线播放| 国产精品不卡视频一区二区| 人妻夜夜爽99麻豆av| 一本—道久久a久久精品蜜桃钙片| 久久97久久精品| 国产v大片淫在线免费观看| 亚洲无线观看免费| 天天躁夜夜躁狠狠久久av| 国产亚洲欧美精品永久| 国产亚洲精品久久久com| 成年av动漫网址| 久久久久国产网址| 免费观看a级毛片全部| 人人妻人人添人人爽欧美一区卜 | 亚洲国产成人一精品久久久| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 久久精品国产亚洲网站| av在线蜜桃| 插阴视频在线观看视频| 色综合色国产| 久久国产乱子免费精品| 成人漫画全彩无遮挡| 国产精品一二三区在线看| 日韩在线高清观看一区二区三区| 中文字幕av成人在线电影| 久久久精品94久久精品| 狂野欧美激情性bbbbbb| 人妻 亚洲 视频| 亚洲精品日韩在线中文字幕| 国产精品爽爽va在线观看网站| 久久久久国产精品人妻一区二区| 小蜜桃在线观看免费完整版高清| 22中文网久久字幕| 国国产精品蜜臀av免费| 亚洲精品成人av观看孕妇| 人人妻人人看人人澡| 免费人妻精品一区二区三区视频| 在线观看国产h片| 免费看日本二区| 黄色一级大片看看| 一级毛片久久久久久久久女| 男女无遮挡免费网站观看| 一级毛片 在线播放| 91久久精品国产一区二区成人| 国产深夜福利视频在线观看| 大陆偷拍与自拍| 国产精品人妻久久久影院| 国产精品久久久久久精品电影小说 | 日韩视频在线欧美| 亚洲国产欧美在线一区| 免费少妇av软件| 最近最新中文字幕免费大全7| 国产久久久一区二区三区| 国产伦精品一区二区三区视频9| 欧美3d第一页| 亚洲精品日韩在线中文字幕| 亚洲欧美清纯卡通| 黄色配什么色好看| 狂野欧美激情性xxxx在线观看| 观看免费一级毛片| 久久精品久久久久久噜噜老黄| 国产中年淑女户外野战色| 中文天堂在线官网| 日本爱情动作片www.在线观看| 乱码一卡2卡4卡精品| 久久精品人妻少妇| 国产av国产精品国产| 中文字幕人妻熟人妻熟丝袜美| 午夜福利高清视频| 久久精品久久精品一区二区三区| 尤物成人国产欧美一区二区三区| 永久免费av网站大全| 亚洲国产精品一区三区| 老熟女久久久| 高清黄色对白视频在线免费看 | 日本与韩国留学比较| 91精品国产国语对白视频| 国产成人a区在线观看| 汤姆久久久久久久影院中文字幕| 啦啦啦视频在线资源免费观看| 久久久久久久亚洲中文字幕| 看非洲黑人一级黄片| 直男gayav资源| 99热网站在线观看| 中文在线观看免费www的网站| 色吧在线观看| 久久久久久久亚洲中文字幕| 高清日韩中文字幕在线| 内地一区二区视频在线| 免费观看a级毛片全部| 99热这里只有是精品在线观看| 中文天堂在线官网| 国产v大片淫在线免费观看| 男人舔奶头视频| 黄色视频在线播放观看不卡| 亚洲国产高清在线一区二区三| 男男h啪啪无遮挡| 国产男女超爽视频在线观看| 国产亚洲午夜精品一区二区久久| 欧美成人精品欧美一级黄| 免费大片黄手机在线观看| 精品国产乱码久久久久久小说| 国内精品宾馆在线| 亚洲最大成人中文| 在线观看三级黄色| 色婷婷久久久亚洲欧美| 男女无遮挡免费网站观看| 国产 一区精品| 天美传媒精品一区二区| 亚洲不卡免费看| 1000部很黄的大片| 啦啦啦中文免费视频观看日本| 夜夜爽夜夜爽视频| 免费播放大片免费观看视频在线观看| a 毛片基地| 网址你懂的国产日韩在线| 国产亚洲午夜精品一区二区久久| 高清毛片免费看| 只有这里有精品99| 久久久色成人| 免费大片黄手机在线观看| 欧美国产精品一级二级三级 | 我要看黄色一级片免费的| 国产成人freesex在线| 精品久久国产蜜桃| 涩涩av久久男人的天堂| 日本色播在线视频| 久久精品夜色国产| 国产精品久久久久久久电影| 最近中文字幕高清免费大全6| 青春草视频在线免费观看| 免费播放大片免费观看视频在线观看| 少妇人妻精品综合一区二区| 九九在线视频观看精品| 在线观看美女被高潮喷水网站| 日韩视频在线欧美| www.色视频.com| 晚上一个人看的免费电影| 久久久久久九九精品二区国产| 好男人视频免费观看在线| 国产在线视频一区二区| 国产精品三级大全| 新久久久久国产一级毛片| 五月伊人婷婷丁香| 免费在线观看成人毛片| 永久免费av网站大全| 午夜福利影视在线免费观看| 精品酒店卫生间| 少妇精品久久久久久久| 美女视频免费永久观看网站| 国产高潮美女av| 黄色怎么调成土黄色| 久久毛片免费看一区二区三区| 五月玫瑰六月丁香| 乱码一卡2卡4卡精品| 在线观看av片永久免费下载| 国产精品久久久久成人av| 色综合色国产|