• <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ī)理之問
    国产aⅴ精品一区二区三区波| 免费在线观看影片大全网站| 成人一区二区视频在线观看| 中出人妻视频一区二区| 国产黄色视频一区二区在线观看 | 在线免费十八禁| 一个人观看的视频www高清免费观看| 午夜免费激情av| 欧美成人a在线观看| 丝袜美腿在线中文| 亚洲内射少妇av| 国产精品国产三级国产av玫瑰| 久久亚洲精品不卡| 国产精品福利在线免费观看| 两个人的视频大全免费| 亚洲av二区三区四区| 成人亚洲欧美一区二区av| 国产亚洲精品综合一区在线观看| 免费一级毛片在线播放高清视频| 直男gayav资源| 国产精品嫩草影院av在线观看| 国产精品久久久久久av不卡| 免费看美女性在线毛片视频| 丝袜美腿在线中文| 99在线视频只有这里精品首页| 在线看三级毛片| 在线看三级毛片| 如何舔出高潮| 国产麻豆成人av免费视频| 91久久精品国产一区二区成人| 五月玫瑰六月丁香| 国产午夜精品久久久久久一区二区三区 | 日韩av在线大香蕉| 久久久久久久午夜电影| 婷婷六月久久综合丁香| 欧美日本亚洲视频在线播放| 久久欧美精品欧美久久欧美| 一进一出抽搐gif免费好疼| 99久久久亚洲精品蜜臀av| 亚洲第一区二区三区不卡| 少妇高潮的动态图| 少妇的逼水好多| 国产探花在线观看一区二区| 亚洲国产精品国产精品| a级毛色黄片| 国产精品嫩草影院av在线观看| 在线观看午夜福利视频| 免费人成视频x8x8入口观看| 欧美性感艳星| 综合色av麻豆| 床上黄色一级片| 国产精品一及| 成人国产麻豆网| 国产午夜精品久久久久久一区二区三区 | 婷婷精品国产亚洲av| 成年版毛片免费区| 成人亚洲欧美一区二区av| 网址你懂的国产日韩在线| 看免费成人av毛片| 婷婷亚洲欧美| 日韩欧美精品免费久久| 亚洲成人久久爱视频| 国产精品美女特级片免费视频播放器| 国产av在哪里看| 精品人妻视频免费看| 欧美成人精品欧美一级黄| 免费一级毛片在线播放高清视频| 国产成人精品久久久久久| 最新中文字幕久久久久| 黄色视频,在线免费观看| 亚洲无线在线观看| 亚洲美女搞黄在线观看 | 久久精品人妻少妇| 日日摸夜夜添夜夜添av毛片| 3wmmmm亚洲av在线观看| 亚洲av五月六月丁香网| 亚洲色图av天堂| 欧美区成人在线视频| 美女被艹到高潮喷水动态| 欧美中文日本在线观看视频| 欧美性猛交╳xxx乱大交人| 91久久精品国产一区二区三区| 国产精品免费一区二区三区在线| 欧美日韩精品成人综合77777| 精品人妻偷拍中文字幕| 国产毛片a区久久久久| 国产熟女欧美一区二区| 日韩欧美在线乱码| 在线国产一区二区在线| 成人永久免费在线观看视频| av女优亚洲男人天堂| 国产精华一区二区三区| 国产精品一二三区在线看| 春色校园在线视频观看| 乱系列少妇在线播放| 美女xxoo啪啪120秒动态图| 精品久久国产蜜桃| 黄色视频,在线免费观看| 美女黄网站色视频| 婷婷六月久久综合丁香| 亚洲综合色惰| 日韩欧美一区二区三区在线观看| 久久欧美精品欧美久久欧美| 久久久欧美国产精品| 欧美一级a爱片免费观看看| 直男gayav资源| 国产av一区在线观看免费| 国产在线精品亚洲第一网站| 日韩大尺度精品在线看网址| 国内久久婷婷六月综合欲色啪| 一卡2卡三卡四卡精品乱码亚洲| 欧美色欧美亚洲另类二区| 成年女人毛片免费观看观看9| 成年女人毛片免费观看观看9| 国产又黄又爽又无遮挡在线| 亚洲人成网站在线观看播放| 成人高潮视频无遮挡免费网站| 成人综合一区亚洲| 人人妻人人看人人澡| 亚洲第一电影网av| 午夜激情福利司机影院| 国产精品女同一区二区软件| 高清日韩中文字幕在线| 亚洲成人久久性| 亚洲精华国产精华液的使用体验 | 日韩一区二区视频免费看| 干丝袜人妻中文字幕| 一本久久中文字幕| 女生性感内裤真人,穿戴方法视频| 精品一区二区免费观看| 在线观看一区二区三区| aaaaa片日本免费| 秋霞在线观看毛片| 国产精品久久久久久久久免| 美女高潮的动态| 成年女人毛片免费观看观看9| 99热这里只有是精品在线观看| 麻豆国产av国片精品| 舔av片在线| 精品国产三级普通话版| av免费在线看不卡| 看非洲黑人一级黄片| 插逼视频在线观看| 亚洲经典国产精华液单| 小说图片视频综合网站| 免费av观看视频| 免费一级毛片在线播放高清视频| 校园春色视频在线观看| 成熟少妇高潮喷水视频| 日韩亚洲欧美综合| 91在线精品国自产拍蜜月| 99热只有精品国产| 男女下面进入的视频免费午夜| 国产精品一及| 天堂动漫精品| 亚洲av美国av| 又粗又爽又猛毛片免费看| 国产老妇女一区| 真实男女啪啪啪动态图| 久久久a久久爽久久v久久| 亚洲天堂国产精品一区在线| 亚洲精品乱码久久久v下载方式| 性欧美人与动物交配| 麻豆一二三区av精品| 男女啪啪激烈高潮av片| 欧美性猛交╳xxx乱大交人| 99久国产av精品| 婷婷精品国产亚洲av| av免费在线看不卡| 色视频www国产| 亚洲最大成人中文| 少妇的逼好多水| 秋霞在线观看毛片| 女同久久另类99精品国产91| a级毛片a级免费在线| 成人二区视频| 看十八女毛片水多多多| 久久人人爽人人爽人人片va| 免费看日本二区| 亚洲av五月六月丁香网| 一级毛片电影观看 | 国国产精品蜜臀av免费| 色av中文字幕| 22中文网久久字幕| 嫩草影院入口| 成人精品一区二区免费| 日韩制服骚丝袜av| 免费在线观看成人毛片| 亚洲成人av在线免费| 国产高清视频在线观看网站| 91精品国产九色| 波野结衣二区三区在线| 亚洲av中文av极速乱| 日韩av不卡免费在线播放| 婷婷精品国产亚洲av在线| 久久亚洲精品不卡| 少妇丰满av| 日韩av不卡免费在线播放| 色在线成人网| 人人妻人人澡欧美一区二区| 中文字幕av在线有码专区| 精品人妻偷拍中文字幕| 国产一区二区激情短视频| 久久精品夜色国产| 一级黄色大片毛片| 嫩草影院精品99| 亚洲美女搞黄在线观看 | 干丝袜人妻中文字幕| 插逼视频在线观看| 女人十人毛片免费观看3o分钟| 国产色婷婷99| 午夜a级毛片| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 波多野结衣高清作品| 国产v大片淫在线免费观看| 亚洲欧美精品自产自拍| 欧美性感艳星| 99精品在免费线老司机午夜| 日韩欧美一区二区三区在线观看| 久久精品国产99精品国产亚洲性色| 国产亚洲精品久久久久久毛片| 国产精品爽爽va在线观看网站| 日本三级黄在线观看| 22中文网久久字幕| 给我免费播放毛片高清在线观看| 大香蕉久久网| 欧美性猛交╳xxx乱大交人| 国产aⅴ精品一区二区三区波| 日产精品乱码卡一卡2卡三| videossex国产| 真人做人爱边吃奶动态| 搡女人真爽免费视频火全软件 | 能在线免费观看的黄片| 黄色视频,在线免费观看| 丝袜美腿在线中文| 蜜桃久久精品国产亚洲av| av在线蜜桃| 中文字幕熟女人妻在线| aaaaa片日本免费| 伦理电影大哥的女人| 精品国产三级普通话版| 久久天躁狠狠躁夜夜2o2o| 中文字幕免费在线视频6| 成人高潮视频无遮挡免费网站| 欧美bdsm另类| 夜夜爽天天搞| 在线播放国产精品三级| av.在线天堂| 亚洲熟妇熟女久久| 级片在线观看| 成人国产麻豆网| 久久精品夜色国产| 国产成人精品久久久久久| 亚洲av第一区精品v没综合| 成人欧美大片| 国产私拍福利视频在线观看| 亚洲美女搞黄在线观看 | 国产淫片久久久久久久久| 亚洲av五月六月丁香网| 日本爱情动作片www.在线观看 | 亚洲欧美日韩卡通动漫| 国产伦一二天堂av在线观看| 淫秽高清视频在线观看| 免费人成在线观看视频色| 亚洲欧美成人综合另类久久久 | 国产精品一区www在线观看| 欧美不卡视频在线免费观看| 国产精品精品国产色婷婷| 精品一区二区三区av网在线观看| 免费av毛片视频| 性欧美人与动物交配| 成人毛片a级毛片在线播放| 国产精品久久久久久av不卡| 亚洲国产精品合色在线| 免费大片18禁| 日韩欧美在线乱码| 成人特级av手机在线观看| 一本精品99久久精品77| 亚洲成人中文字幕在线播放| 性欧美人与动物交配| 日韩人妻高清精品专区| 成人av在线播放网站| 99视频精品全部免费 在线| 日韩成人av中文字幕在线观看 | а√天堂www在线а√下载| 国产精品野战在线观看| 长腿黑丝高跟| 日韩大尺度精品在线看网址| 美女黄网站色视频| 国产大屁股一区二区在线视频| 国产 一区精品| 最近的中文字幕免费完整| 久久精品91蜜桃| 性插视频无遮挡在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲性夜色夜夜综合| 天堂网av新在线| 国产蜜桃级精品一区二区三区| 欧美成人精品欧美一级黄| 99国产极品粉嫩在线观看| 亚洲国产精品成人久久小说 | 香蕉av资源在线| 欧美成人免费av一区二区三区| 日韩精品青青久久久久久| 午夜福利视频1000在线观看| av天堂中文字幕网| videossex国产| 欧美激情久久久久久爽电影| 日日干狠狠操夜夜爽| 内射极品少妇av片p| 永久网站在线| 久久久色成人| 嫩草影视91久久| 午夜福利视频1000在线观看| 亚洲国产高清在线一区二区三| 欧美成人a在线观看| 日韩,欧美,国产一区二区三区 | 长腿黑丝高跟| 性欧美人与动物交配| 在线观看免费视频日本深夜| 乱人视频在线观看| 国产精品一区二区免费欧美| 久久久午夜欧美精品| 欧美成人a在线观看| 国产欧美日韩精品一区二区| 91av网一区二区| 亚洲五月天丁香| 一级毛片aaaaaa免费看小| 黄片wwwwww| 国产欧美日韩一区二区精品| 国产探花极品一区二区| 99热网站在线观看| 国产精品嫩草影院av在线观看| 真人做人爱边吃奶动态| 在线观看一区二区三区| 久久久精品欧美日韩精品| 久久草成人影院| 欧美+亚洲+日韩+国产| 简卡轻食公司| 91午夜精品亚洲一区二区三区| 久久久国产成人免费| 亚洲国产高清在线一区二区三| 国模一区二区三区四区视频| 变态另类成人亚洲欧美熟女| 午夜老司机福利剧场| 亚洲精品一卡2卡三卡4卡5卡| 日日干狠狠操夜夜爽| 天堂网av新在线| 久久亚洲精品不卡| 久久精品国产亚洲网站| 亚洲色图av天堂| 91久久精品国产一区二区三区| 欧美国产日韩亚洲一区| 日产精品乱码卡一卡2卡三| 久久国内精品自在自线图片| 久久久欧美国产精品| 色噜噜av男人的天堂激情| 一本久久中文字幕| av在线亚洲专区| 亚洲综合色惰| 在线看三级毛片| 变态另类成人亚洲欧美熟女| 偷拍熟女少妇极品色| 天堂网av新在线| 搡女人真爽免费视频火全软件 | 亚洲色图av天堂| 少妇猛男粗大的猛烈进出视频 | 丰满的人妻完整版| 丝袜喷水一区| 十八禁网站免费在线| 久久人妻av系列| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 日韩制服骚丝袜av| 真实男女啪啪啪动态图| 久久人人爽人人爽人人片va| 国产在线男女| 免费看美女性在线毛片视频| 内射极品少妇av片p| 搡老妇女老女人老熟妇| 色综合色国产| 六月丁香七月| 熟女电影av网| av女优亚洲男人天堂| 国产高潮美女av| 精品一区二区三区视频在线| 此物有八面人人有两片| 少妇人妻精品综合一区二区 | 日韩制服骚丝袜av| 真实男女啪啪啪动态图| 变态另类丝袜制服| 波多野结衣巨乳人妻| 九九热线精品视视频播放| 国产在线男女| 亚洲av一区综合| av在线亚洲专区| 午夜亚洲福利在线播放| 精品福利观看| 97热精品久久久久久| 一个人看的www免费观看视频| 毛片一级片免费看久久久久| 在线观看免费视频日本深夜| 1024手机看黄色片| 大型黄色视频在线免费观看| 欧美3d第一页| 搡女人真爽免费视频火全软件 | 欧美bdsm另类| 亚洲天堂国产精品一区在线| 国产亚洲91精品色在线| 国产精品人妻久久久久久| 深爱激情五月婷婷| 久久鲁丝午夜福利片| 欧美成人a在线观看| 亚洲人成网站在线播放欧美日韩| 国产欧美日韩一区二区精品| 人人妻,人人澡人人爽秒播| 免费看av在线观看网站| 国产高清激情床上av| 一本精品99久久精品77| 国产黄a三级三级三级人| 久久久欧美国产精品| 一级毛片久久久久久久久女| 国产白丝娇喘喷水9色精品| 一边摸一边抽搐一进一小说| 三级国产精品欧美在线观看| 国产精品一区www在线观看| 欧美又色又爽又黄视频| 特大巨黑吊av在线直播| 久久这里只有精品中国| av在线老鸭窝| 精品久久久久久久久av| 国产精品不卡视频一区二区| 亚洲av第一区精品v没综合| 精品福利观看| 老司机影院成人| 成人毛片a级毛片在线播放| 在线看三级毛片| 18禁黄网站禁片免费观看直播| 成年av动漫网址| 亚洲一级一片aⅴ在线观看| 亚洲av成人精品一区久久| АⅤ资源中文在线天堂| 91在线观看av| 国产高清有码在线观看视频| 一级黄色大片毛片| 午夜爱爱视频在线播放| 亚洲国产精品sss在线观看| 国产久久久一区二区三区| 久久精品国产清高在天天线| 天堂av国产一区二区熟女人妻| 嫩草影院精品99| 又黄又爽又免费观看的视频| 一本久久中文字幕| 91午夜精品亚洲一区二区三区| 成人特级黄色片久久久久久久| 成人无遮挡网站| 亚洲中文字幕日韩| 1024手机看黄色片| 日本与韩国留学比较| 人妻久久中文字幕网| 热99在线观看视频| 天堂网av新在线| or卡值多少钱| 大香蕉久久网| 午夜精品在线福利| 日日啪夜夜撸| 亚洲不卡免费看| 波多野结衣高清作品| 国产亚洲精品综合一区在线观看| 亚洲国产精品成人综合色| 在线观看免费视频日本深夜| 最近2019中文字幕mv第一页| 日韩中字成人| 十八禁国产超污无遮挡网站| 久久精品综合一区二区三区| 黄色配什么色好看| 午夜免费激情av| 日韩av不卡免费在线播放| 中出人妻视频一区二区| 日韩三级伦理在线观看| 日韩欧美在线乱码| 欧美成人a在线观看| 我的老师免费观看完整版| 免费在线观看成人毛片| 日韩制服骚丝袜av| a级一级毛片免费在线观看| or卡值多少钱| 亚洲专区国产一区二区| 午夜爱爱视频在线播放| 3wmmmm亚洲av在线观看| 国内精品久久久久精免费| 久久久久久国产a免费观看| 不卡视频在线观看欧美| 亚洲欧美清纯卡通| 久久久午夜欧美精品| 欧美高清性xxxxhd video| 久久人人爽人人爽人人片va| 两性午夜刺激爽爽歪歪视频在线观看| 男人的好看免费观看在线视频| 亚洲激情五月婷婷啪啪| 高清毛片免费看| 变态另类丝袜制服| 亚洲国产精品合色在线| 麻豆国产av国片精品| 国产精品无大码| 最近最新中文字幕大全电影3| 国产精品爽爽va在线观看网站| 免费无遮挡裸体视频| 亚洲欧美日韩高清专用| 亚洲成人中文字幕在线播放| 久久久精品大字幕| 综合色丁香网| 在线天堂最新版资源| 久久精品国产亚洲网站| 国产精品久久久久久精品电影| 午夜视频国产福利| 嫩草影院入口| 嫩草影视91久久| 国产久久久一区二区三区| 有码 亚洲区| 久久久久国产精品人妻aⅴ院| 久久久色成人| www日本黄色视频网| 久久人妻av系列| 欧美性感艳星| 欧美三级亚洲精品| 99热全是精品| 一进一出好大好爽视频| 免费一级毛片在线播放高清视频| 夜夜爽天天搞| 波多野结衣高清无吗| 亚洲图色成人| 成人特级av手机在线观看| 人人妻人人澡人人爽人人夜夜 | 精品熟女少妇av免费看| 国产单亲对白刺激| 久久久久久久久久黄片| 国产午夜精品论理片| 久久久久国内视频| 国产av在哪里看| 国产综合懂色| 国产精品精品国产色婷婷| 插逼视频在线观看| 干丝袜人妻中文字幕| 亚洲自偷自拍三级| 尾随美女入室| 男人的好看免费观看在线视频| 亚洲美女黄片视频| av黄色大香蕉| av在线天堂中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 两个人的视频大全免费| 亚洲欧美精品综合久久99| 99在线人妻在线中文字幕| 国产精品国产高清国产av| 久久久国产成人免费| 最新中文字幕久久久久| av专区在线播放| 成人午夜高清在线视频| 国产乱人视频| 天堂影院成人在线观看| 日本免费a在线| av.在线天堂| 国产v大片淫在线免费观看| 九色成人免费人妻av| 日韩精品中文字幕看吧| 日韩 亚洲 欧美在线| 成人三级黄色视频| 亚洲欧美精品综合久久99| 一进一出好大好爽视频| 美女被艹到高潮喷水动态| а√天堂www在线а√下载| 老熟妇仑乱视频hdxx| 直男gayav资源| 欧美性感艳星| 在线a可以看的网站| 久久九九热精品免费| 亚洲精品亚洲一区二区| 欧美成人精品欧美一级黄| 精品久久久久久久久av| 国产午夜精品论理片| av卡一久久| 悠悠久久av| 亚洲美女搞黄在线观看 | 国产精品国产高清国产av| 51国产日韩欧美| 日日撸夜夜添| 不卡视频在线观看欧美| 可以在线观看的亚洲视频| 国产高清激情床上av| 日韩成人伦理影院| 精品欧美国产一区二区三| 久久亚洲精品不卡| 久久亚洲国产成人精品v| 久久久色成人| 婷婷六月久久综合丁香| 国产伦一二天堂av在线观看| 国产淫片久久久久久久久| 男女做爰动态图高潮gif福利片| 亚洲国产日韩欧美精品在线观看| 欧美激情在线99| 国内精品美女久久久久久| 免费在线观看成人毛片| 亚洲欧美成人综合另类久久久 | 人妻制服诱惑在线中文字幕| 高清毛片免费观看视频网站| 老熟妇仑乱视频hdxx| 99久国产av精品| 深夜精品福利| 国产亚洲91精品色在线| 国产伦一二天堂av在线观看| а√天堂www在线а√下载| 高清午夜精品一区二区三区 | 欧美高清性xxxxhd video|