• <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| 在线天堂最新版资源| 亚洲 国产 在线| 亚洲国产欧洲综合997久久,| 老鸭窝网址在线观看| 日韩亚洲欧美综合| 中文字幕久久专区| 精品人妻偷拍中文字幕| 亚洲熟妇中文字幕五十中出| 美女cb高潮喷水在线观看| 日本三级黄在线观看| 91九色精品人成在线观看| 日韩人妻高清精品专区| 中文字幕熟女人妻在线| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址| 久久精品夜夜夜夜夜久久蜜豆| 色综合站精品国产| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 国产免费一级a男人的天堂| 他把我摸到了高潮在线观看| 国产精华一区二区三区| 怎么达到女性高潮| 国产爱豆传媒在线观看| 亚洲无线在线观看| 最新在线观看一区二区三区| 熟妇人妻久久中文字幕3abv| 精品人妻视频免费看| 亚洲性夜色夜夜综合| 国产成人aa在线观看| 直男gayav资源| 一个人看的www免费观看视频| 美女大奶头视频| 免费大片18禁| 99热这里只有精品一区| 国产老妇女一区| 欧美午夜高清在线| 日韩高清综合在线| 欧美一级a爱片免费观看看| 国产视频一区二区在线看| 草草在线视频免费看| 嫩草影院新地址| 亚洲第一欧美日韩一区二区三区| 一个人看的www免费观看视频| avwww免费| 亚洲av第一区精品v没综合| 欧美丝袜亚洲另类 | 97碰自拍视频| 国内毛片毛片毛片毛片毛片| 成人亚洲精品av一区二区| 国内精品一区二区在线观看| 国产精品久久久久久精品电影| 老鸭窝网址在线观看| 国产久久久一区二区三区| 少妇高潮的动态图| 国产精品,欧美在线| 中国美女看黄片| 亚洲男人的天堂狠狠| 99久久九九国产精品国产免费| 一区福利在线观看| 午夜免费激情av| 精品久久久久久成人av| 久久香蕉精品热| 欧美午夜高清在线| 日本一二三区视频观看| 看片在线看免费视频| 一进一出抽搐gif免费好疼| АⅤ资源中文在线天堂| xxxwww97欧美| 国产白丝娇喘喷水9色精品| 18禁黄网站禁片免费观看直播| 国产伦人伦偷精品视频| 蜜桃亚洲精品一区二区三区| 黄色配什么色好看| 国产精品免费一区二区三区在线| 久99久视频精品免费| 午夜视频国产福利| 无遮挡黄片免费观看| 天美传媒精品一区二区| 国产av一区在线观看免费| 免费av不卡在线播放| 国产精品久久久久久精品电影| 99riav亚洲国产免费| 国产又黄又爽又无遮挡在线| ponron亚洲| 90打野战视频偷拍视频| 亚洲第一欧美日韩一区二区三区| 熟女电影av网| 国产精品日韩av在线免费观看| 无人区码免费观看不卡| 国产精品亚洲美女久久久| 亚洲精品成人久久久久久| av视频在线观看入口| 精品99又大又爽又粗少妇毛片 | 久久精品夜夜夜夜夜久久蜜豆| 色精品久久人妻99蜜桃| 免费观看精品视频网站| 成人三级黄色视频| 亚洲国产精品成人综合色| 亚洲精品乱码久久久v下载方式| 亚洲av二区三区四区| 免费电影在线观看免费观看| 免费观看的影片在线观看| 久久久精品大字幕| 欧美又色又爽又黄视频| 欧美在线一区亚洲| 夜夜躁狠狠躁天天躁| 国产高清视频在线播放一区| 亚洲狠狠婷婷综合久久图片| 女人被狂操c到高潮| 又黄又爽又免费观看的视频| 久久婷婷人人爽人人干人人爱| 十八禁人妻一区二区| 国产 一区 欧美 日韩| 动漫黄色视频在线观看| 深爱激情五月婷婷| 日韩欧美在线二视频| 国产老妇女一区| 长腿黑丝高跟| 99久国产av精品| 精品福利观看| 久久精品综合一区二区三区| 欧美极品一区二区三区四区| 国产精品亚洲一级av第二区| 亚洲狠狠婷婷综合久久图片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99国产精品一区二区蜜桃av| 国产成人a区在线观看| 成人av一区二区三区在线看| 国产精品自产拍在线观看55亚洲| 久久久国产成人精品二区| 伦理电影大哥的女人| 国产伦精品一区二区三区视频9| 成人无遮挡网站| 久9热在线精品视频| 免费在线观看日本一区| 亚洲乱码一区二区免费版| 国产高清有码在线观看视频| 成人美女网站在线观看视频| 悠悠久久av| 欧美一区二区亚洲| 精品99又大又爽又粗少妇毛片 | 久久久久久久久久黄片| 午夜激情福利司机影院| 黄色一级大片看看| 变态另类丝袜制服| 色综合欧美亚洲国产小说| 一级黄色大片毛片| 一个人免费在线观看的高清视频| 国产精品久久久久久久久免 | 搡老妇女老女人老熟妇| 国内精品久久久久精免费| 久久精品人妻少妇| 嫩草影院入口| 9191精品国产免费久久| 一进一出好大好爽视频| 欧美黄色淫秽网站| 日本黄大片高清| 国产精品三级大全| 免费看日本二区| av天堂中文字幕网| 国内毛片毛片毛片毛片毛片| 亚洲欧美日韩高清在线视频| 五月伊人婷婷丁香| 99在线视频只有这里精品首页| 国产精品人妻久久久久久| 五月伊人婷婷丁香| 在现免费观看毛片| 欧美日本视频| 一二三四社区在线视频社区8| 精品无人区乱码1区二区| 激情在线观看视频在线高清| 黄片小视频在线播放| 在线免费观看不下载黄p国产 | 男女做爰动态图高潮gif福利片| 好男人在线观看高清免费视频| 国产伦一二天堂av在线观看| 丰满人妻一区二区三区视频av| 51国产日韩欧美| 欧美精品啪啪一区二区三区| 久久精品影院6| 女人被狂操c到高潮| 99视频精品全部免费 在线| 国产精品日韩av在线免费观看| 老司机福利观看| 人妻夜夜爽99麻豆av| 免费看美女性在线毛片视频| 人妻丰满熟妇av一区二区三区| 国产三级在线视频| 好男人在线观看高清免费视频| 内射极品少妇av片p| 亚洲激情在线av| 国产成人av教育| 国产精品一区二区三区四区免费观看 | 两人在一起打扑克的视频| 亚洲色图av天堂| 国产欧美日韩精品亚洲av| 久久久国产成人精品二区| 亚洲电影在线观看av| 欧美性猛交╳xxx乱大交人| 久久久久精品国产欧美久久久| 亚洲精品亚洲一区二区| 18禁裸乳无遮挡免费网站照片| 午夜日韩欧美国产| 国产精品免费一区二区三区在线| 亚洲av电影不卡..在线观看| 91午夜精品亚洲一区二区三区 | 别揉我奶头 嗯啊视频| 精品人妻偷拍中文字幕| 日本一二三区视频观看| 精品福利观看| 午夜福利免费观看在线| 日本一本二区三区精品| 国产高清激情床上av| 能在线免费观看的黄片| 亚洲国产高清在线一区二区三| 一个人免费在线观看的高清视频| 又爽又黄a免费视频| 国产男靠女视频免费网站| 久久久久久久精品吃奶| 欧美日韩国产亚洲二区| 亚洲av成人av| 国产综合懂色| 亚洲成人中文字幕在线播放| 美女 人体艺术 gogo| 一个人看的www免费观看视频| 九九热线精品视视频播放| 又爽又黄a免费视频| 亚洲乱码一区二区免费版| 给我免费播放毛片高清在线观看| 久久精品国产亚洲av涩爱 | 亚洲av中文字字幕乱码综合| 欧美黑人欧美精品刺激| 中出人妻视频一区二区| 97碰自拍视频| 久久久国产成人免费| 成年版毛片免费区| 嫩草影院新地址| av专区在线播放| 国产主播在线观看一区二区| 久久草成人影院| 麻豆久久精品国产亚洲av| 日日干狠狠操夜夜爽| 国产在线男女| 天堂动漫精品| 精品久久久久久久久久久久久| 亚洲乱码一区二区免费版| 久9热在线精品视频| 国产黄片美女视频| а√天堂www在线а√下载| 成人一区二区视频在线观看| 久久精品国产清高在天天线| 麻豆av噜噜一区二区三区| 亚洲中文字幕日韩| 色播亚洲综合网| 国产精品不卡视频一区二区 | 国产亚洲精品综合一区在线观看| 亚洲av.av天堂| 蜜桃亚洲精品一区二区三区| 自拍偷自拍亚洲精品老妇| 欧美色视频一区免费| 国产精品一区二区三区四区免费观看 | 国产精品综合久久久久久久免费| 男插女下体视频免费在线播放| 成人永久免费在线观看视频| 国产精品自产拍在线观看55亚洲| 小蜜桃在线观看免费完整版高清| 他把我摸到了高潮在线观看| 国产久久久一区二区三区| 国产三级在线视频| 免费在线观看成人毛片| 999久久久精品免费观看国产| 欧美在线一区亚洲| 少妇熟女aⅴ在线视频| 亚洲国产欧美人成| 国产精品电影一区二区三区| 国产精品,欧美在线| 色av中文字幕| 97人妻精品一区二区三区麻豆| 可以在线观看的亚洲视频| 成人特级黄色片久久久久久久| 最近中文字幕高清免费大全6 | 激情在线观看视频在线高清| 国产精品自产拍在线观看55亚洲| 看黄色毛片网站| 99热只有精品国产| 亚洲天堂国产精品一区在线| 麻豆成人午夜福利视频| 不卡一级毛片| 亚洲五月婷婷丁香| 日本成人三级电影网站| 午夜福利在线观看吧| 欧美乱色亚洲激情| 最近视频中文字幕2019在线8| 99热这里只有是精品50| 一a级毛片在线观看| 亚洲成a人片在线一区二区| 精品一区二区三区人妻视频| 男人的好看免费观看在线视频| 人妻制服诱惑在线中文字幕| 最近中文字幕高清免费大全6 | 欧美最新免费一区二区三区 | 亚洲中文日韩欧美视频| 美女黄网站色视频| 中文在线观看免费www的网站| 亚洲精品成人久久久久久| 欧美色欧美亚洲另类二区| 国产精品综合久久久久久久免费| 久9热在线精品视频| 国产精品久久久久久久电影| www.www免费av| 国产精品一区二区三区四区免费观看 | 国产男靠女视频免费网站| 精品人妻熟女av久视频| 大型黄色视频在线免费观看| 无人区码免费观看不卡| 此物有八面人人有两片| 国产成人aa在线观看| 人妻夜夜爽99麻豆av| 男人的好看免费观看在线视频| 亚洲av二区三区四区| av视频在线观看入口| 国产伦在线观看视频一区| 成年女人永久免费观看视频| 成人欧美大片| 床上黄色一级片| 99热这里只有是精品50| 一二三四社区在线视频社区8| 老女人水多毛片| 色综合亚洲欧美另类图片| 欧美高清性xxxxhd video| 深夜a级毛片| 成年免费大片在线观看| 亚洲精品一区av在线观看| 亚洲成人久久爱视频| 成人av一区二区三区在线看| 精品人妻熟女av久视频| 国产三级中文精品| 99久久精品国产亚洲精品| 97碰自拍视频| 国产视频内射| 国产真实伦视频高清在线观看 | 国产午夜精品论理片| 色综合站精品国产| 亚洲成av人片在线播放无| 亚洲精品粉嫩美女一区| 免费一级毛片在线播放高清视频| a级毛片免费高清观看在线播放| 窝窝影院91人妻| 长腿黑丝高跟| 天堂影院成人在线观看| 一区二区三区高清视频在线| 美女高潮喷水抽搐中文字幕| xxxwww97欧美| 色在线成人网| 精品日产1卡2卡| 亚洲成人中文字幕在线播放| 亚洲人成网站在线播| 国产不卡一卡二| 婷婷丁香在线五月| 搡女人真爽免费视频火全软件 | 欧美激情在线99| 欧美精品啪啪一区二区三区| 夜夜爽天天搞| 日本成人三级电影网站| 毛片女人毛片| 精品国产亚洲在线| 99精品久久久久人妻精品| 精品一区二区三区人妻视频| 亚洲成人久久爱视频| 大型黄色视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 欧美又色又爽又黄视频| 国产成人欧美在线观看| 久久精品国产自在天天线| 精品99又大又爽又粗少妇毛片 | 18美女黄网站色大片免费观看| 国产白丝娇喘喷水9色精品| 日日干狠狠操夜夜爽| 日韩精品中文字幕看吧| 国模一区二区三区四区视频| 午夜福利在线观看免费完整高清在 | 啦啦啦韩国在线观看视频| 亚洲内射少妇av| 可以在线观看毛片的网站| 天堂√8在线中文| 啦啦啦韩国在线观看视频| 国产一区二区在线av高清观看| 99热精品在线国产| 一个人免费在线观看的高清视频| 久久久久久久久久成人| 99久久无色码亚洲精品果冻| 全区人妻精品视频| 99在线视频只有这里精品首页| 每晚都被弄得嗷嗷叫到高潮| 能在线免费观看的黄片| 国产三级在线视频| 校园春色视频在线观看| 精品一区二区免费观看| av在线老鸭窝| 琪琪午夜伦伦电影理论片6080| 欧美高清成人免费视频www| 美女xxoo啪啪120秒动态图 | 久久久色成人| 此物有八面人人有两片| 男人狂女人下面高潮的视频| 97人妻精品一区二区三区麻豆| 成人亚洲精品av一区二区| 久久精品国产亚洲av香蕉五月| 成熟少妇高潮喷水视频| 国产91精品成人一区二区三区| 国产综合懂色| 国产一区二区亚洲精品在线观看| 国产精品亚洲美女久久久| 成人午夜高清在线视频| 熟女电影av网| 国产成人aa在线观看| 亚洲专区中文字幕在线| 在现免费观看毛片| 九九在线视频观看精品| 日本 av在线| 日本免费一区二区三区高清不卡| 性插视频无遮挡在线免费观看| 国产伦精品一区二区三区四那| avwww免费| 国产淫片久久久久久久久 | 少妇高潮的动态图| 18禁在线播放成人免费| 在线播放国产精品三级| 亚洲精品日韩av片在线观看| 欧美日韩福利视频一区二区| 小说图片视频综合网站| 黄色视频,在线免费观看| 亚洲av成人精品一区久久| av国产免费在线观看| 免费在线观看亚洲国产| 色哟哟哟哟哟哟| 午夜福利在线观看吧| 欧美日韩瑟瑟在线播放| 亚洲av电影在线进入| 高潮久久久久久久久久久不卡| 久久99热这里只有精品18| 免费观看精品视频网站| 久久精品夜夜夜夜夜久久蜜豆| 91麻豆av在线| 国产精品久久久久久精品电影| 高潮久久久久久久久久久不卡| 淫秽高清视频在线观看| 亚洲欧美日韩东京热| 亚洲七黄色美女视频| 亚洲专区国产一区二区| 能在线免费观看的黄片| 美女免费视频网站| 观看免费一级毛片| 欧美不卡视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲av二区三区四区| 久久久色成人| 亚洲乱码一区二区免费版| 亚洲不卡免费看| 亚洲国产欧洲综合997久久,| 亚洲无线观看免费| 免费人成视频x8x8入口观看| 国内精品美女久久久久久| 18禁黄网站禁片午夜丰满| 美女大奶头视频| 精品人妻一区二区三区麻豆 | 成人国产一区最新在线观看| 久久精品国产99精品国产亚洲性色| 午夜a级毛片| 亚洲天堂国产精品一区在线| 久久久久久大精品| 久久久久国产精品人妻aⅴ院| 亚洲欧美清纯卡通| 精品人妻1区二区| 欧美极品一区二区三区四区| 一二三四社区在线视频社区8| 在线国产一区二区在线| 男女视频在线观看网站免费| 欧美一区二区精品小视频在线| 国产精品亚洲av一区麻豆| 婷婷色综合大香蕉| 亚洲在线自拍视频| 亚洲av一区综合| 女同久久另类99精品国产91| 国产中年淑女户外野战色| 男女下面进入的视频免费午夜| 淫妇啪啪啪对白视频| 国产三级黄色录像| 欧美不卡视频在线免费观看| 日韩亚洲欧美综合| 亚洲精品日韩av片在线观看| 精品久久久久久久久久免费视频| 国产精品自产拍在线观看55亚洲| 麻豆成人av在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品sss在线观看| 黄色丝袜av网址大全| 久久精品综合一区二区三区| 老熟妇乱子伦视频在线观看| 亚洲人成电影免费在线| 3wmmmm亚洲av在线观看| 精品国内亚洲2022精品成人| 两个人视频免费观看高清| 一级黄色大片毛片| x7x7x7水蜜桃| 他把我摸到了高潮在线观看| 国产亚洲精品综合一区在线观看| 一级毛片久久久久久久久女| 国产成+人综合+亚洲专区| 一个人免费在线观看电影| 一级黄片播放器| 色综合婷婷激情| 亚洲av电影在线进入| 国产亚洲精品综合一区在线观看| 国产精品日韩av在线免费观看| 国产精品久久久久久久久免 | 精品乱码久久久久久99久播| 99久久成人亚洲精品观看| 色噜噜av男人的天堂激情| 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久,| 51国产日韩欧美| 日韩精品青青久久久久久| 国产精品美女特级片免费视频播放器| 国产一区二区三区视频了| 成年人黄色毛片网站| 国内精品美女久久久久久| 国产亚洲av嫩草精品影院| 在线天堂最新版资源| 国产精华一区二区三区| 变态另类成人亚洲欧美熟女| 欧美成人一区二区免费高清观看| 国产乱人伦免费视频| 97热精品久久久久久| 日韩中文字幕欧美一区二区| 亚洲av免费高清在线观看| 悠悠久久av| av女优亚洲男人天堂| 好男人电影高清在线观看| 中文字幕精品亚洲无线码一区| av在线观看视频网站免费| 国产精品久久久久久久久免 | 国产亚洲av嫩草精品影院| 悠悠久久av| 少妇人妻精品综合一区二区 | 亚洲欧美精品综合久久99| 天堂动漫精品| 国产精品国产高清国产av| 简卡轻食公司| 国产免费av片在线观看野外av| 亚洲成av人片在线播放无| 97热精品久久久久久| 国产野战对白在线观看| 国产在线精品亚洲第一网站| 亚洲七黄色美女视频| 国产爱豆传媒在线观看| 99热这里只有精品一区| or卡值多少钱| 欧美三级亚洲精品| 婷婷丁香在线五月| 久久精品综合一区二区三区| 在线观看免费视频日本深夜| 亚洲一区高清亚洲精品| www.999成人在线观看| 免费看光身美女| 久99久视频精品免费| 九九久久精品国产亚洲av麻豆| 亚洲成av人片在线播放无| 国产亚洲av嫩草精品影院| 五月伊人婷婷丁香| 成人精品一区二区免费| www.www免费av| 亚洲国产精品sss在线观看| 国产精品电影一区二区三区| 国产精品久久视频播放| 色在线成人网| 可以在线观看毛片的网站| 18美女黄网站色大片免费观看| 欧美日本亚洲视频在线播放| 啦啦啦观看免费观看视频高清| 啪啪无遮挡十八禁网站| 人人妻,人人澡人人爽秒播| 99热这里只有是精品50| 日韩av在线大香蕉| 久久久久久九九精品二区国产| 在线播放国产精品三级| 观看美女的网站| 1000部很黄的大片| 高清在线国产一区| 一个人免费在线观看电影| 亚洲成人免费电影在线观看| 精品日产1卡2卡| 国产一区二区三区在线臀色熟女| 国产午夜福利久久久久久| 免费黄网站久久成人精品 | 免费av不卡在线播放| 国产69精品久久久久777片| 97超视频在线观看视频| 国产精品国产高清国产av| 在线看三级毛片| 亚洲av成人av| 久久国产精品影院| 亚洲第一欧美日韩一区二区三区| 久久久久久国产a免费观看| 日韩中文字幕欧美一区二区| 中文字幕免费在线视频6| 色哟哟·www| 亚洲电影在线观看av| 欧美高清性xxxxhd video| 我要搜黄色片|