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

    基于貝葉斯統(tǒng)計(jì)的汽車電子器件壽命分析

    2023-10-21 07:08:22劉岳歐陽英圖葉雯
    時(shí)代汽車 2023年18期

    劉岳 歐陽英圖 葉雯

    摘 要:加速壽命試驗(yàn)是利用加速應(yīng)力水平下的壽命特征去外推評(píng)估正常應(yīng)力水平下的壽命特征的試驗(yàn)技術(shù)。針對(duì)在汽車智能化發(fā)展中起著重要作用的汽車電子器件,本文將服從對(duì)數(shù)正態(tài)分布的樣本在步進(jìn)應(yīng)力加速壽命試驗(yàn)中的統(tǒng)計(jì)分析結(jié)合貝葉斯統(tǒng)計(jì)推斷,得出了基于貝葉斯統(tǒng)計(jì)的加速壽命試驗(yàn)統(tǒng)計(jì)分析方法,從而實(shí)現(xiàn)了對(duì)汽車電子器件的壽命估計(jì)。本文通過有機(jī)發(fā)光二極管的案例應(yīng)用,驗(yàn)證了貝葉斯對(duì)數(shù)正態(tài)分布加速壽命試驗(yàn)統(tǒng)計(jì)分析方法的可行性和有效性,對(duì)相關(guān)理論研究和工程試驗(yàn)的開展具有借鑒意義。

    關(guān)鍵詞:汽車電子產(chǎn)品 壽命分析 加速壽命試驗(yàn) 貝葉斯統(tǒng)計(jì)推斷

    1 緒論

    隨著汽車電氣化乃至智能化的發(fā)展,汽車電子器件在車身各關(guān)鍵設(shè)備上的應(yīng)用日漸廣泛[1]。汽車電子器件的工作狀態(tài)、功能、壽命與汽車的正常行駛息息相關(guān),若出現(xiàn)問題,輕則造成財(cái)產(chǎn)損失,重則造成人員傷亡。因此,對(duì)汽車電子器件進(jìn)行壽命分析,具有重大的實(shí)際意義。

    1.1 加速壽命試驗(yàn)與加速模型

    為了快速地暴露產(chǎn)品的薄弱環(huán)節(jié),在較高應(yīng)力下以更短的試驗(yàn)時(shí)間推斷正常應(yīng)力下的壽命特征,常采取加速壽命試驗(yàn)(Life Accelerated Testing,ALT)。即在失效機(jī)理不變的基礎(chǔ)上,通過加速模型,利用加速應(yīng)力水平下的壽命特征去外推評(píng)估正常應(yīng)力水平下的壽命特征的試驗(yàn)技術(shù)。加速壽命試驗(yàn)方法因其可縮短試驗(yàn)時(shí)間、提高試驗(yàn)效率、降低試驗(yàn)成本等優(yōu)勢(shì)已經(jīng)被廣泛應(yīng)用于各類工程實(shí)際問題之中[2]。

    為了能夠利用ALT中搜集到的產(chǎn)品壽命信息外推產(chǎn)品在正常應(yīng)力條件下的壽命特征,必須建立產(chǎn)品壽命特征與加速應(yīng)力水平之間的關(guān)系,即加速模型。常用的加速模型分為物理模型和統(tǒng)計(jì)模型,具體有阿倫尼斯模型、艾琳模型、廣義艾琳模型、沖蝕磨損模型、逆冪律模型、Coffin-Manson模型、Norris-Landzberg模型等[3]。

    ALT的統(tǒng)計(jì)分析是通過估計(jì)壽命分布函數(shù)的參數(shù)和確定加速模型的參數(shù),從而外推評(píng)估正常應(yīng)力水平S0下的壽命特征。

    1.2 貝葉斯理論

    在工程和實(shí)際試驗(yàn)中,對(duì)于待估計(jì)參數(shù)常常會(huì)有一定的現(xiàn)有經(jīng)驗(yàn)和信息,為了利用好這一部分信息,同時(shí)通過新的數(shù)據(jù)對(duì)已有信息進(jìn)行更新,則常用貝葉斯統(tǒng)計(jì)方法[4]進(jìn)行統(tǒng)計(jì)推斷。

    p(θ|y)稱為后驗(yàn)密度函數(shù);p(θ)稱為先驗(yàn)密度函數(shù);m(y)是數(shù)據(jù)的邊沿密度函數(shù);f(y|θ)是數(shù)據(jù)的抽樣密度函數(shù)。

    由于汽車為批量生產(chǎn)的產(chǎn)品,因此其電子器件也具有相當(dāng)多的歷史信息,故采用基于貝葉斯統(tǒng)計(jì)ALT分析,能夠更準(zhǔn)確地評(píng)估汽車電子器件壽命,并對(duì)產(chǎn)品已有信息進(jìn)行更新。

    2 基于貝葉斯統(tǒng)計(jì)的ALT分析方法

    2.1 基本假設(shè)

    以對(duì)數(shù)正態(tài)分布場(chǎng)合下的步進(jìn)應(yīng)力加速壽命試驗(yàn)(Step-Stress Accelerated Life Testing, SSALT)為例,試驗(yàn)數(shù)據(jù)的統(tǒng)計(jì)分析是在以下四個(gè)假定下進(jìn)行的[2]。

    假設(shè)1:在正常應(yīng)力水平S0和k個(gè)加速應(yīng)力水平S1

    式中,?為對(duì)數(shù)均值,σ2為對(duì)數(shù)方差。

    其分布函數(shù)為

    假設(shè)2:產(chǎn)品在加速應(yīng)力水平S1

    假設(shè)3:產(chǎn)品的某壽命特征ξ與施加的應(yīng)力水平Si有如下的對(duì)數(shù)線性關(guān)系

    在對(duì)數(shù)正態(tài)分布中我們用中位壽命e?來表征壽命特征,所以常用的加速模型例如阿倫尼斯模型、逆冪律模型、指數(shù)率模型的加速模型表達(dá)式為

    假設(shè)4:產(chǎn)品的殘余壽命僅依賴于已累積失效的部分和當(dāng)時(shí)的應(yīng)力水平,而與累積方式無關(guān)(Nelson累積失效假定)。

    2.2 假設(shè)檢驗(yàn)

    假設(shè)檢驗(yàn)是指根據(jù)樣本提供的信息來推斷總體特性的某些假設(shè)是否成立,步驟如下:

    1.提出原假設(shè)H0:即需要檢驗(yàn)的假設(shè)內(nèi)容;

    2.引進(jìn)統(tǒng)計(jì)量:根據(jù)H0的內(nèi)容選取合適的統(tǒng)計(jì)量;

    3.確定統(tǒng)計(jì)量的精確分布或漸近分布;

    4.根據(jù)觀測(cè)到的樣本值算出統(tǒng)計(jì)量的值;

    5.確定適當(dāng)?shù)娘@著性水平α;

    6.根據(jù)統(tǒng)計(jì)量分布,由顯著性水平確定臨界點(diǎn)P{統(tǒng)計(jì)量計(jì)算值>Kα}=α

    現(xiàn)針對(duì)基本假設(shè)1和2做分布檢驗(yàn)和參數(shù)檢驗(yàn)。

    2.2.1 壽命分布的假設(shè)檢驗(yàn)

    分布檢驗(yàn)針對(duì)基本假設(shè)1,使用正態(tài)分布偏峰度檢驗(yàn)[2]進(jìn)行對(duì)數(shù)正態(tài)分布分布檢驗(yàn)。

    偏度描述分布密度函數(shù)的對(duì)稱程度,分布密度越對(duì)稱,偏度越??;峰度描述分布密度函數(shù)的陡峭程度,分布密度越陡峭,峰度越大。偏度和峰度分別定義為

    式中,?3為三階中心矩,?4為四階中心矩,σ2為方差。以正態(tài)分布為例,由于正態(tài)分布N(?,σ2)的偏度為0,峰度為3,因此可以通過樣本偏度和峰度是否接近0和3來判斷數(shù)據(jù)是否服從正態(tài)分布,具體判斷方法如下:

    1)偏度檢驗(yàn)

    檢驗(yàn)原假設(shè)H0:Cs=0,即原假設(shè)認(rèn)為分布密度函數(shù)是對(duì)稱的,若偏度Cs估計(jì)的絕對(duì)值超過它的1-α分位數(shù),則在顯著水平α下拒絕原假設(shè),檢驗(yàn)統(tǒng)計(jì)量|Cs|的1-α分位數(shù)可查表[5]得到,偏度檢驗(yàn)適用于樣本量n≥8。

    2)峰度檢驗(yàn)

    3)檢驗(yàn)原假設(shè)H0:Ce=3,即原假設(shè)認(rèn)為分布密度函數(shù)是對(duì)稱的,若Ce大于3,則說明峰度過度,這時(shí)備擇假設(shè)H1:Ce>3,若Ce小于3,則說明峰度不足,這時(shí)備擇假設(shè)H1:Ce<3,峰度檢驗(yàn)適用于樣本量n≥8。其中統(tǒng)計(jì)量Ce的臨界值可查[5]得到。

    2.2.2 形狀參數(shù)相等的假設(shè)檢驗(yàn)

    形狀參數(shù)檢驗(yàn)?zāi)康氖谴_認(rèn)加速應(yīng)力是否改變了產(chǎn)品的失效機(jī)理,該檢驗(yàn)針對(duì)基本假設(shè)2。

    對(duì)于對(duì)數(shù)正態(tài)分布,可以把形狀參數(shù)的檢驗(yàn)問題轉(zhuǎn)化成方差一致性的檢驗(yàn)問題[6]。對(duì)于假設(shè)檢驗(yàn)

    根據(jù)試驗(yàn)數(shù)據(jù)用最好線性無偏估計(jì)(BLUE)計(jì)算σi的估計(jì)值。構(gòu)造巴特利特檢驗(yàn)統(tǒng)計(jì)量,即

    式中,lri,ni為的方差系數(shù)。則在H0成立下,B2/C近似的服從自由度為K-1的χ2分布,給定顯著性水平?,若,拒絕H0;否則接受H0。

    2.3 特征數(shù)據(jù)的統(tǒng)計(jì)分析

    2.3.1 數(shù)據(jù)折算

    在應(yīng)力水平Si下,產(chǎn)品工作ti時(shí)間內(nèi)累積失效的概率Fi(ti)相當(dāng)于此種產(chǎn)品在Sj下工作tj時(shí)間內(nèi)累積失效的概率Fj(tj)[2]。由假設(shè)2,可得在應(yīng)力水平Si下的工作時(shí)間τi折算到應(yīng)力水平Sj下的折算時(shí)間τij為

    記a1=0,。則產(chǎn)品在Si下的真實(shí)壽命數(shù)據(jù)為

    其中b為假設(shè)3中加速方程:中的待估參數(shù);φ(Si)為已知函數(shù),記為φi。即在每一個(gè)應(yīng)力水平Si下,ri個(gè)失效產(chǎn)品的折算壽命為

    式中:;,除外,它們都是b的函數(shù)。

    2.3.2 貝葉斯統(tǒng)計(jì)推斷

    聯(lián)立壽命分布和加速模型:

    得下式:

    可知,待估參數(shù)為a、b、σ。

    1)確定先驗(yàn)分布

    進(jìn)行貝葉斯統(tǒng)計(jì)推斷的第一步是為待估參數(shù)確定先驗(yàn)分布。根據(jù)待估參數(shù)的先驗(yàn)信息的多少,可分為無信息先驗(yàn)分布和有信息先驗(yàn)分布。

    a)無信息先驗(yàn)分布

    若試驗(yàn)前有關(guān)待估參數(shù)a、b、σ的信息相對(duì)較少,可以假設(shè)其先驗(yàn)分布為均勻分布,例如:a~U[x1, y1],b~U[x2, y2],σ~U[x3,y3];xi,yi的取值可由現(xiàn)有信息、經(jīng)驗(yàn)決定。確定無信息先驗(yàn)分布的方法還有Jeffrey法則[7]等。

    b)有信息先驗(yàn)分布

    若試驗(yàn)前有比較充足的關(guān)于待估參數(shù)的先驗(yàn)信息時(shí),可根據(jù)先驗(yàn)信息確定先驗(yàn)分布。常見的先驗(yàn)信息主要來源于:物理/化學(xué)理論、仿真分析、工程試驗(yàn)和質(zhì)量控制測(cè)試結(jié)果、一些通用的可靠性數(shù)據(jù)、相似產(chǎn)品信息以及專家經(jīng)驗(yàn)。

    確定有信息先驗(yàn)分布的方法[8]主要有:直方圖法、選定先驗(yàn)密度分布函數(shù)形式再估計(jì)其超參數(shù)、變分讀法與定分度法以及多層先驗(yàn)等。

    2)觀測(cè)數(shù)據(jù)先驗(yàn)分布的更新

    先驗(yàn)分布p(θ)代表了試驗(yàn)前對(duì)待估參數(shù)的初始估計(jì)。在得到數(shù)據(jù)后,可以利用“新的信息”對(duì)先驗(yàn)分布進(jìn)行更新,得到后驗(yàn)分布。根據(jù)貝葉斯定理,后驗(yàn)分布是由似然函數(shù)和先驗(yàn)分布相乘得到的。

    后驗(yàn)分布似然函數(shù)×先驗(yàn)分布

    在本文中,似然函數(shù)即

    具體計(jì)算是由馬爾科夫鏈蒙特卡羅(MCMC)算法得出[9]:通過模擬的方法直接從后驗(yàn)分布中生成參數(shù)向量的仿真樣本,基于樣本進(jìn)行后續(xù)統(tǒng)計(jì)推斷。

    在本文算例中,具體計(jì)算通過WinBUGS軟件實(shí)現(xiàn)。在WinBUGS中輸入三個(gè)參數(shù)的先驗(yàn)分布、似然函數(shù),以及試驗(yàn)失效數(shù)據(jù),運(yùn)行后能夠給出參數(shù)的仿真樣本,即一系列服從對(duì)應(yīng)的后驗(yàn)分布的a、b、σ的值。

    3)分位壽命的計(jì)算

    第三步則是通過估計(jì)參數(shù)計(jì)算出分位壽命tp=f(a,b,p),可分為點(diǎn)估計(jì)和區(qū)間估計(jì)。

    a)點(diǎn)估計(jì)

    主要有三種方法:一是極大后驗(yàn)密度估計(jì)(maximum a posteriori estimate, MAP),即后驗(yàn)密度函數(shù)取極大值處相應(yīng)的待估參數(shù)值;二是取參數(shù)的后驗(yàn)均值;三是取參數(shù)的后驗(yàn)中位數(shù)。通過點(diǎn)估計(jì)得到a,b的值后,即可計(jì)算出分位壽命tp。

    b)區(qū)間估計(jì)

    區(qū)間估計(jì)時(shí)可以通過a, b的樣本值,兩兩配對(duì)各算出一個(gè)分位壽命——tp=f(a,b,p);然后對(duì)其取置信區(qū)間即可(例如,取置信度為α的置信區(qū)間,則除去tp取值中前α/2和后α/2的取值,所得的[tpmin, tpmax]即為置信區(qū)間)。

    3 案例分析

    本文借鑒文獻(xiàn)[10]中的案例和試驗(yàn)數(shù)據(jù)進(jìn)行貝葉斯對(duì)數(shù)正態(tài)步進(jìn)ALT統(tǒng)計(jì)分析。

    試驗(yàn)對(duì)象采用汽車電子器件——有機(jī)發(fā)光二極管(Organic Light-Emitting Diode, OLED)。對(duì)于OLED,影響其壽命的主要因素是電流,因而選擇電流為加速應(yīng)力進(jìn)行分析。根據(jù)試驗(yàn)設(shè)計(jì)要求,在對(duì)OLED電子產(chǎn)品進(jìn)行試驗(yàn)階段,施加的應(yīng)力,即電流,應(yīng)該介于最大應(yīng)力和最小應(yīng)力中,且最大應(yīng)力不應(yīng)該超過極限應(yīng)力,并且所施加的應(yīng)力水平由目前的制造技術(shù)可以達(dá)到??紤]到試驗(yàn)時(shí)長(zhǎng)與成本,采用SSALT。此外,隨著電流的增大,不能夠?qū)LED電子產(chǎn)品引發(fā)新的失效機(jī)理。OLED電子產(chǎn)品的失效判據(jù)為:當(dāng)OLED的亮度降低到原亮度的50%以下時(shí),認(rèn)為這個(gè)樣本失效,并記錄相應(yīng)的失效時(shí)間。

    3.1 基本假設(shè)

    進(jìn)行SSALT試驗(yàn)數(shù)據(jù)的統(tǒng)計(jì)分析的基本假設(shè)同2.1節(jié)。其中,在對(duì)數(shù)正態(tài)分布中用中位壽命來表征壽命特征,本文中的應(yīng)力為電流,所以選擇逆冪律模型作為加速模型:

    3.2 數(shù)據(jù)折算處理

    1.原始數(shù)據(jù)

    試驗(yàn)中一共投入了32個(gè)樣本進(jìn)行試驗(yàn),假定正常的工作條件中,應(yīng)力I0=3.20mA。在SSALT中選擇了四個(gè)應(yīng)力水平,分別為:I1=9.64mA,I2=12.36mA,I3=17.09mA,I4=22.58mA,在每個(gè)應(yīng)力水平下收集八個(gè)失效數(shù)據(jù),示意圖如圖1,試驗(yàn)數(shù)據(jù)記錄如表1。

    對(duì)數(shù)據(jù)進(jìn)行折算:

    其中,τji為加速因子,且,β=-1.8154,,,。將所有的數(shù)據(jù)折算到應(yīng)力為I1的水平。

    3.3 基于貝葉斯的SSALT統(tǒng)計(jì)分析

    將壽命分布的模型、先驗(yàn)信息以及折算數(shù)據(jù)導(dǎo)入到winbugs中進(jìn)行求解,求解的結(jié)果如圖2所示。

    通過winbugs的求解,得到參數(shù)a和b的貝葉斯估計(jì),即后驗(yàn)樣本值,代入到壽命分布函數(shù)中即可計(jì)算在正常應(yīng)力下的各種壽命指標(biāo)的估計(jì)。

    顯然,0.5分位壽命是參數(shù)a和b的函數(shù),即t0.5=(a,b);每一組a,b可以計(jì)算一個(gè)p=0.5的分位壽命,取均值,得到0.5分位壽命t0.5=24100h。

    先對(duì)參數(shù)a和b進(jìn)行點(diǎn)估計(jì),常用的方法有眾數(shù)/平均數(shù)和中位數(shù)。

    眾數(shù):求得a和b的眾數(shù),帶入到壽命分布中,得到相關(guān)的分位壽命,t0.5=20647h,這個(gè)結(jié)果和分位壽命均值的結(jié)果有所差異,通過分析a和b的概率密度函數(shù)(如圖3)可知原因?yàn)椋簠?shù)a和b通過數(shù)據(jù)的迭代會(huì)逐漸向著最終的分布形式過渡,但是由于目前的試驗(yàn)數(shù)據(jù)較少,概率密度函數(shù)都還具有均勻分布和正態(tài)分布的特點(diǎn),因而處于過渡階段。

    平均數(shù):求得a和b的平均數(shù),帶入到壽命分布中,得到相關(guān)的分位壽命t0.5=16400h;

    中位數(shù):求得a和b的中位數(shù),帶入到壽命分布中,得到相關(guān)的分位壽命t0.5=16200h;

    0.5分位點(diǎn)壽命的置信區(qū)間:對(duì)于winbugs生成的參數(shù)a和b的后驗(yàn)樣本,帶入壽命分布函數(shù)中求解一組0.5分位點(diǎn)的壽命值t0.5,對(duì)這些t0.5進(jìn)行分布擬合,可以得到如圖4的概率密度函數(shù)。

    相應(yīng)地,根據(jù)置信區(qū)間的計(jì)算方法,可以得到置信度為0.9的t0.5的置信區(qū)間為[3639.1041, 76083.1355]。

    4 總結(jié)

    本文針對(duì)汽車電子器件,介紹了基于貝葉斯統(tǒng)計(jì)推斷的對(duì)數(shù)正態(tài)步進(jìn)應(yīng)力加速壽命試驗(yàn)統(tǒng)計(jì)分析方法。利用已有歷史數(shù)據(jù)的信息,同時(shí)通過新的數(shù)據(jù)對(duì)已有信息進(jìn)行更新;結(jié)合加速壽命試驗(yàn)統(tǒng)計(jì)分析方法,對(duì)服從對(duì)數(shù)正態(tài)的汽車電子器件壽命進(jìn)行了點(diǎn)估計(jì)及區(qū)間估計(jì)。通過案例應(yīng)用,驗(yàn)證了該方法的可行性和有效性。同時(shí),基于貝葉斯統(tǒng)計(jì)推斷的對(duì)數(shù)正態(tài)步進(jìn)應(yīng)力加速壽命試驗(yàn)統(tǒng)計(jì)分析方法也可以進(jìn)一步推廣應(yīng)用于各類加速壽命試驗(yàn)統(tǒng)計(jì)分析中,對(duì)工程應(yīng)用具有指導(dǎo)意義。

    項(xiàng)目資助:

    湖南省教育廳科學(xué)研究項(xiàng)目,項(xiàng)目編號(hào)21C1659

    湖南省自然基金科教聯(lián)合項(xiàng)目,項(xiàng)目編號(hào)2022JJ60113。

    參考文獻(xiàn):

    [1]李莉. 基于售后數(shù)據(jù)的汽車電子元器件壽命分布估計(jì)研究[J]. 甘肅科技, 2021, 37(6):2.

    [2]趙宇. 可靠性數(shù)據(jù)分析[M]. 國(guó)防工業(yè)出版社, 2011.

    [3]Nelson W. Accelerated Testing: Statistical Models, Test Plans, and Data Analysis[M].? 2008.

    [4]王婧, 鮑貴. 貝葉斯統(tǒng)計(jì)與傳統(tǒng)統(tǒng)計(jì)方法的比較[J]. 統(tǒng)計(jì)與決策, 2021(037-001).

    [5]李良巧. 可靠性工程師手冊(cè)[M]. 中國(guó)人民大學(xué)出版社, 2012.

    [6]程琮,范華, Levene方差齊性檢驗(yàn). 中國(guó)衛(wèi)生統(tǒng)計(jì), 2005. 22(6): p. 408-408.

    [7]Tiao, G. C., and G. E. P. Box. Some Comments on "Bayes" Estimators. American Statistician 27.1. 1973. pp:12-14.

    [8]柴慧敏, 趙昀瑤, 方敏. 利用先驗(yàn)正態(tài)分布的貝葉斯網(wǎng)絡(luò)參數(shù)學(xué)習(xí)[J]. 系統(tǒng)工程與電子技術(shù), 2018, 40(10):6.

    [9]Hamada, Michael S., et al. Bayesian Reliability. Bayesian reliability. Springer, 2008.

    [10]Zhang J, Cheng G, Xiao C, et al. Accelerated life test of white OLED based on lognormal distribution[J]. Indian Journal of Pure and Applied Physics. 2014, 52(10):671-677.

    在线观看www视频免费| 久久久久久久国产电影| 免费人成在线观看视频色| 国产在线视频一区二区| 亚洲精品视频女| 精品午夜福利在线看| av在线播放精品| 日韩一区二区三区影片| 80岁老熟妇乱子伦牲交| 日本黄大片高清| 久久久久久伊人网av| 18禁在线播放成人免费| 免费看av在线观看网站| 麻豆成人av视频| 男女啪啪激烈高潮av片| 一边亲一边摸免费视频| 性色avwww在线观看| 欧美日韩亚洲高清精品| 日本黄色日本黄色录像| 久久国产乱子免费精品| 一级毛片久久久久久久久女| 超碰97精品在线观看| 在线观看免费高清a一片| 成人国产av品久久久| 国产女主播在线喷水免费视频网站| 国产精品99久久久久久久久| 欧美成人午夜免费资源| 桃花免费在线播放| 永久免费av网站大全| 午夜福利在线观看免费完整高清在| 欧美人与善性xxx| 亚洲情色 制服丝袜| 久久韩国三级中文字幕| 国产视频内射| 啦啦啦中文免费视频观看日本| 男女啪啪激烈高潮av片| 人人妻人人澡人人爽人人夜夜| 97超碰精品成人国产| 老女人水多毛片| 国产精品一二三区在线看| 夜夜骑夜夜射夜夜干| av卡一久久| 亚洲,一卡二卡三卡| 精品久久久久久久久av| 亚洲av福利一区| 观看美女的网站| 人妻系列 视频| kizo精华| 国产精品偷伦视频观看了| 国产亚洲最大av| 国产精品一区www在线观看| 国产免费又黄又爽又色| 亚洲欧美一区二区三区黑人 | 亚洲精品成人av观看孕妇| 赤兔流量卡办理| 纵有疾风起免费观看全集完整版| 免费av中文字幕在线| a 毛片基地| 一级毛片aaaaaa免费看小| 交换朋友夫妻互换小说| 成人18禁高潮啪啪吃奶动态图 | 国产精品久久久久久精品古装| 久久久精品免费免费高清| 22中文网久久字幕| 男女边吃奶边做爰视频| 搡老乐熟女国产| 在线观看三级黄色| 天天操日日干夜夜撸| 毛片一级片免费看久久久久| 久久精品熟女亚洲av麻豆精品| 日日摸夜夜添夜夜添av毛片| 亚洲四区av| 亚洲综合精品二区| 一级毛片电影观看| 久久午夜综合久久蜜桃| 亚洲精品国产成人久久av| 十分钟在线观看高清视频www | 成人免费观看视频高清| 九九在线视频观看精品| 国产精品一区二区性色av| 亚洲精品,欧美精品| 午夜激情久久久久久久| 亚洲,一卡二卡三卡| 中文字幕av电影在线播放| 一级毛片我不卡| 最近的中文字幕免费完整| 亚洲精品国产av蜜桃| 亚洲精品一区蜜桃| 大又大粗又爽又黄少妇毛片口| 精品99又大又爽又粗少妇毛片| 又爽又黄a免费视频| 国产在线男女| a级一级毛片免费在线观看| 熟女av电影| 亚洲成色77777| 久久久国产精品麻豆| 全区人妻精品视频| 午夜日本视频在线| 日本-黄色视频高清免费观看| 亚洲av中文av极速乱| 日本与韩国留学比较| 成人国产麻豆网| 中国三级夫妇交换| 国产精品一区二区在线不卡| 亚洲精品国产av成人精品| 香蕉精品网在线| 久久亚洲国产成人精品v| 国产精品成人在线| 久久青草综合色| 久久久久久久国产电影| 成年美女黄网站色视频大全免费 | av播播在线观看一区| 久久ye,这里只有精品| 韩国av在线不卡| 国产成人午夜福利电影在线观看| 国产欧美日韩综合在线一区二区 | 高清视频免费观看一区二区| 美女国产视频在线观看| 国产精品嫩草影院av在线观看| 精品视频人人做人人爽| 在线观看免费日韩欧美大片 | 丰满饥渴人妻一区二区三| 国产伦理片在线播放av一区| 亚洲国产av新网站| 亚洲精品国产av成人精品| 国产色婷婷99| 国内少妇人妻偷人精品xxx网站| 狂野欧美激情性xxxx在线观看| 另类精品久久| 精品酒店卫生间| 精品一区二区免费观看| 亚洲精品日本国产第一区| 蜜桃在线观看..| 在线观看免费日韩欧美大片 | 久久久久久久国产电影| 久久人妻熟女aⅴ| 国产亚洲午夜精品一区二区久久| 亚洲国产欧美日韩在线播放 | 精品人妻熟女av久视频| 免费大片黄手机在线观看| 桃花免费在线播放| 亚洲国产色片| 成人免费观看视频高清| 成人特级av手机在线观看| 免费看日本二区| 波野结衣二区三区在线| 九色成人免费人妻av| 亚洲人与动物交配视频| 男人添女人高潮全过程视频| 国产成人freesex在线| 80岁老熟妇乱子伦牲交| av.在线天堂| 少妇丰满av| 丰满迷人的少妇在线观看| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久久丰满| 永久网站在线| 色视频在线一区二区三区| 亚洲人成网站在线播| 亚洲欧美一区二区三区黑人 | 99热这里只有是精品在线观看| 欧美+日韩+精品| 最新中文字幕久久久久| 亚洲精品国产av成人精品| h视频一区二区三区| 欧美亚洲 丝袜 人妻 在线| 亚洲色图综合在线观看| 观看免费一级毛片| 深夜a级毛片| 一级毛片我不卡| 国产av一区二区精品久久| 五月伊人婷婷丁香| av不卡在线播放| 日日摸夜夜添夜夜添av毛片| 日韩中文字幕视频在线看片| 久久久久久久大尺度免费视频| 人体艺术视频欧美日本| 精品一区二区免费观看| 亚洲美女视频黄频| 亚洲精品国产av蜜桃| 午夜激情福利司机影院| 丰满饥渴人妻一区二区三| 亚洲第一av免费看| 欧美精品高潮呻吟av久久| 国产淫语在线视频| 国产精品一区二区在线观看99| 一区二区三区精品91| 一本色道久久久久久精品综合| 校园人妻丝袜中文字幕| 久久精品国产亚洲网站| tube8黄色片| 久久鲁丝午夜福利片| 寂寞人妻少妇视频99o| 亚洲av成人精品一二三区| 久久久精品94久久精品| 老司机亚洲免费影院| 黄色视频在线播放观看不卡| 大又大粗又爽又黄少妇毛片口| 精品人妻偷拍中文字幕| 视频区图区小说| 国产极品粉嫩免费观看在线 | 美女内射精品一级片tv| 国产伦理片在线播放av一区| 五月玫瑰六月丁香| 中文字幕亚洲精品专区| 成人毛片60女人毛片免费| 天美传媒精品一区二区| 精品视频人人做人人爽| 久久久久久久久大av| 午夜激情福利司机影院| 国产免费福利视频在线观看| 久久韩国三级中文字幕| 欧美高清成人免费视频www| 亚洲国产av新网站| 亚洲精品456在线播放app| 欧美xxⅹ黑人| av国产精品久久久久影院| 成人美女网站在线观看视频| 精品国产一区二区久久| 国产无遮挡羞羞视频在线观看| 夜夜看夜夜爽夜夜摸| 久久午夜综合久久蜜桃| 中文在线观看免费www的网站| av在线观看视频网站免费| 亚洲va在线va天堂va国产| 日日摸夜夜添夜夜添av毛片| 少妇的逼好多水| 亚洲精品国产av蜜桃| 中文天堂在线官网| 男男h啪啪无遮挡| 国产又色又爽无遮挡免| 国产男人的电影天堂91| 亚洲美女搞黄在线观看| 久久国产乱子免费精品| 韩国av在线不卡| 在线天堂最新版资源| 免费大片18禁| 国产高清三级在线| 少妇被粗大猛烈的视频| 色视频www国产| 亚洲欧美精品自产自拍| 男男h啪啪无遮挡| 国产又色又爽无遮挡免| 2018国产大陆天天弄谢| 女性被躁到高潮视频| 亚洲精品一二三| 丰满人妻一区二区三区视频av| 亚洲精品中文字幕在线视频 | 国产免费一区二区三区四区乱码| 国产成人freesex在线| 亚洲欧美中文字幕日韩二区| 亚洲久久久国产精品| xxx大片免费视频| 亚洲精品国产av蜜桃| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人久久小说| 亚洲成人av在线免费| 久久人人爽av亚洲精品天堂| 少妇的逼水好多| 十分钟在线观看高清视频www | 美女大奶头黄色视频| 一区二区三区免费毛片| 国产成人精品福利久久| 色视频在线一区二区三区| 免费人成在线观看视频色| 观看av在线不卡| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 国产探花极品一区二区| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人精品一区二区| 91午夜精品亚洲一区二区三区| 亚洲欧洲国产日韩| 久久午夜福利片| 国产精品久久久久久精品古装| 国产又色又爽无遮挡免| 日韩av不卡免费在线播放| 亚州av有码| 一本一本综合久久| 精品一品国产午夜福利视频| 午夜福利网站1000一区二区三区| xxx大片免费视频| 麻豆成人午夜福利视频| 精品一品国产午夜福利视频| 日韩av在线免费看完整版不卡| 成人漫画全彩无遮挡| 在线免费观看不下载黄p国产| 亚洲性久久影院| 国产精品国产三级专区第一集| 美女xxoo啪啪120秒动态图| 色视频www国产| 久久鲁丝午夜福利片| 国产成人免费观看mmmm| 国产成人aa在线观看| 在线天堂最新版资源| 2018国产大陆天天弄谢| 久久精品国产亚洲网站| 亚洲精品自拍成人| 少妇 在线观看| 午夜福利在线观看免费完整高清在| 91久久精品国产一区二区成人| 亚洲无线观看免费| 国产亚洲最大av| 伊人久久国产一区二区| 亚洲欧美成人综合另类久久久| 三级国产精品片| 成人毛片a级毛片在线播放| 亚洲国产最新在线播放| 久久久久久久久久人人人人人人| 日日爽夜夜爽网站| 久久久久国产网址| 国产高清有码在线观看视频| 中国美白少妇内射xxxbb| 欧美精品国产亚洲| 午夜影院在线不卡| 十分钟在线观看高清视频www | 伊人久久精品亚洲午夜| 人人妻人人添人人爽欧美一区卜| 麻豆成人av视频| 精品一区在线观看国产| 久久毛片免费看一区二区三区| 国产一区亚洲一区在线观看| 99热这里只有是精品50| 国产一区亚洲一区在线观看| 欧美精品高潮呻吟av久久| 精品久久久久久电影网| 欧美成人午夜免费资源| 又大又黄又爽视频免费| 亚洲精品国产av蜜桃| 国产真实伦视频高清在线观看| 99热全是精品| 亚洲精品,欧美精品| 欧美xxxx性猛交bbbb| 永久免费av网站大全| 春色校园在线视频观看| 欧美三级亚洲精品| 99久国产av精品国产电影| 大陆偷拍与自拍| 一区二区三区四区激情视频| 精品人妻熟女毛片av久久网站| 亚洲av电影在线观看一区二区三区| 日韩欧美一区视频在线观看 | 久久久欧美国产精品| 亚洲内射少妇av| 精品久久久噜噜| av免费在线看不卡| 又爽又黄a免费视频| 草草在线视频免费看| 亚洲av不卡在线观看| 日本色播在线视频| 天天躁夜夜躁狠狠久久av| 国产精品熟女久久久久浪| 色吧在线观看| 在线观看一区二区三区激情| 国产亚洲91精品色在线| 午夜免费观看性视频| 精品久久久久久久久av| 国产在线男女| 男女国产视频网站| 久久亚洲国产成人精品v| 亚洲欧洲精品一区二区精品久久久 | 久久6这里有精品| 内射极品少妇av片p| 久久人人爽av亚洲精品天堂| 内射极品少妇av片p| 欧美日韩视频高清一区二区三区二| 亚洲人成网站在线观看播放| 久久影院123| 99久久精品一区二区三区| 又粗又硬又长又爽又黄的视频| 国产免费又黄又爽又色| 色吧在线观看| 国产免费又黄又爽又色| 国产亚洲午夜精品一区二区久久| 国产成人午夜福利电影在线观看| 18禁在线播放成人免费| 国产在线视频一区二区| 最近中文字幕高清免费大全6| 久久精品国产亚洲av涩爱| 久久久久久久精品精品| 欧美激情极品国产一区二区三区 | 一边亲一边摸免费视频| av视频免费观看在线观看| 久久久精品94久久精品| a级片在线免费高清观看视频| 成人特级av手机在线观看| 少妇的逼水好多| 国产高清国产精品国产三级| 在线精品无人区一区二区三| 国产成人免费无遮挡视频| 女人久久www免费人成看片| 91精品国产国语对白视频| h日本视频在线播放| 国产精品久久久久久av不卡| 精品国产一区二区三区久久久樱花| 一级毛片我不卡| 一级毛片黄色毛片免费观看视频| 亚洲精品日本国产第一区| 亚洲精品aⅴ在线观看| 黄色视频在线播放观看不卡| 免费高清在线观看视频在线观看| 亚洲精品456在线播放app| 精品久久久久久久久亚洲| 啦啦啦啦在线视频资源| av福利片在线| 亚洲不卡免费看| 国内少妇人妻偷人精品xxx网站| 国产精品99久久99久久久不卡 | 午夜久久久在线观看| 免费人成在线观看视频色| 2018国产大陆天天弄谢| 日韩在线高清观看一区二区三区| 日韩成人av中文字幕在线观看| 国产一区二区三区综合在线观看 | 色视频www国产| 日韩欧美 国产精品| 少妇人妻精品综合一区二区| 亚洲熟女精品中文字幕| 国产在线一区二区三区精| 一本—道久久a久久精品蜜桃钙片| 久久国产精品男人的天堂亚洲 | 视频区图区小说| 国产精品久久久久成人av| 国产一区二区三区综合在线观看 | 日韩人妻高清精品专区| 欧美3d第一页| 九九久久精品国产亚洲av麻豆| 欧美成人午夜免费资源| 国产成人91sexporn| 亚洲av成人精品一区久久| 日本猛色少妇xxxxx猛交久久| 亚洲人成网站在线观看播放| 亚洲高清免费不卡视频| 夜夜看夜夜爽夜夜摸| 老司机亚洲免费影院| 久久免费观看电影| 人妻少妇偷人精品九色| 国产熟女欧美一区二区| 亚洲欧美日韩卡通动漫| 岛国毛片在线播放| 一本大道久久a久久精品| 永久免费av网站大全| 少妇 在线观看| 91aial.com中文字幕在线观看| 中文字幕亚洲精品专区| 丰满人妻一区二区三区视频av| 日韩电影二区| 天堂8中文在线网| 性高湖久久久久久久久免费观看| 国产亚洲精品久久久com| 国内少妇人妻偷人精品xxx网站| 三级国产精品片| 精品一品国产午夜福利视频| 精品久久国产蜜桃| 久久久久久久久久成人| 久久午夜福利片| 永久网站在线| 午夜激情福利司机影院| 精品久久久久久久久亚洲| av黄色大香蕉| 国产精品欧美亚洲77777| 国产成人一区二区在线| 亚洲成人手机| 亚洲欧美一区二区三区黑人 | 国产真实伦视频高清在线观看| 久久久久久人妻| 免费不卡的大黄色大毛片视频在线观看| 在线观看免费视频网站a站| 亚洲国产精品一区三区| 国产精品一区二区在线观看99| 欧美 亚洲 国产 日韩一| 91久久精品电影网| 亚洲精品中文字幕在线视频 | 国产乱来视频区| 久久精品国产自在天天线| 老女人水多毛片| 免费av中文字幕在线| 亚洲在久久综合| 国模一区二区三区四区视频| 久久热精品热| 三上悠亚av全集在线观看 | 久久这里有精品视频免费| 一级二级三级毛片免费看| 久久久久久久久久久免费av| 麻豆成人午夜福利视频| 久久国产乱子免费精品| 久久久国产欧美日韩av| 丰满乱子伦码专区| 日韩,欧美,国产一区二区三区| 人妻人人澡人人爽人人| 亚洲精品日韩在线中文字幕| 国产爽快片一区二区三区| 国产一区亚洲一区在线观看| 亚洲国产精品999| 我要看日韩黄色一级片| 久久久欧美国产精品| 精品一区在线观看国产| 80岁老熟妇乱子伦牲交| 午夜久久久在线观看| 大又大粗又爽又黄少妇毛片口| 久久久久久久久大av| 亚洲av在线观看美女高潮| 色婷婷av一区二区三区视频| 老女人水多毛片| 美女国产视频在线观看| 精品亚洲成a人片在线观看| 亚洲国产欧美在线一区| 男女无遮挡免费网站观看| 黑人猛操日本美女一级片| 日韩一区二区视频免费看| 久久久精品94久久精品| 熟女av电影| 日韩av免费高清视频| 国产高清三级在线| 日韩中文字幕视频在线看片| 一级爰片在线观看| 十分钟在线观看高清视频www | 女人久久www免费人成看片| 日本欧美国产在线视频| 中文欧美无线码| 日本91视频免费播放| 国产高清国产精品国产三级| 夫妻性生交免费视频一级片| 97精品久久久久久久久久精品| 又大又黄又爽视频免费| 日韩欧美一区视频在线观看 | 成人亚洲欧美一区二区av| 水蜜桃什么品种好| 少妇人妻 视频| 日韩欧美一区视频在线观看 | 国产av码专区亚洲av| 国产精品国产三级国产专区5o| 少妇人妻久久综合中文| 免费看不卡的av| 一区二区三区免费毛片| 亚洲精品日韩在线中文字幕| 国产乱人偷精品视频| 狂野欧美激情性xxxx在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲欧美一区二区三区国产| 搡女人真爽免费视频火全软件| 人妻 亚洲 视频| 日本爱情动作片www.在线观看| 一级毛片黄色毛片免费观看视频| 性色av一级| 桃花免费在线播放| 欧美老熟妇乱子伦牲交| 能在线免费看毛片的网站| 国产黄色免费在线视频| 精品亚洲成a人片在线观看| 我要看日韩黄色一级片| 最近中文字幕2019免费版| a级毛色黄片| 99久久中文字幕三级久久日本| 国产伦理片在线播放av一区| 最近手机中文字幕大全| 国产黄频视频在线观看| 久久久a久久爽久久v久久| 国产综合精华液| 制服丝袜香蕉在线| 久热久热在线精品观看| 国产男女内射视频| 成人免费观看视频高清| 少妇人妻精品综合一区二区| 久久久久久伊人网av| 久久av网站| 99久国产av精品国产电影| 99精国产麻豆久久婷婷| 国产亚洲一区二区精品| 777米奇影视久久| 十八禁网站网址无遮挡 | 少妇人妻一区二区三区视频| 男人添女人高潮全过程视频| 午夜福利视频精品| 2022亚洲国产成人精品| 亚洲国产毛片av蜜桃av| 国产av国产精品国产| 一级,二级,三级黄色视频| 中文字幕人妻熟人妻熟丝袜美| 在线观看美女被高潮喷水网站| 国产av精品麻豆| 婷婷色av中文字幕| 午夜久久久在线观看| 国产高清不卡午夜福利| 亚洲精品中文字幕在线视频 | 亚洲国产毛片av蜜桃av| 国产高清有码在线观看视频| 亚洲激情五月婷婷啪啪| 久久久a久久爽久久v久久| av在线播放精品| 亚洲一区二区三区欧美精品| 亚洲国产毛片av蜜桃av| 久久 成人 亚洲| 免费高清在线观看视频在线观看| 亚洲av成人精品一区久久| 日日摸夜夜添夜夜爱| 天天躁夜夜躁狠狠久久av| 香蕉精品网在线| 高清av免费在线| 色网站视频免费| 秋霞在线观看毛片| 五月玫瑰六月丁香| 亚洲经典国产精华液单| 热re99久久精品国产66热6| 免费久久久久久久精品成人欧美视频 | 国产免费视频播放在线视频| 一级毛片黄色毛片免费观看视频| 国产一级毛片在线| 久久久精品94久久精品| 久久人妻熟女aⅴ| 自线自在国产av| 人人澡人人妻人| 天天躁夜夜躁狠狠久久av| 久久久国产欧美日韩av| 国产成人免费无遮挡视频|