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

    基于自適應(yīng)帶寬核密度估計(jì)的載荷外推方法研究

    2021-02-01 12:02:04牛文鐵才福友付景靜
    關(guān)鍵詞:四叉樹密度估計(jì)概率密度

    牛文鐵 才福友 付景靜

    (天津大學(xué)機(jī)構(gòu)理論與裝備設(shè)計(jì)教育部重點(diǎn)實(shí)驗(yàn)室, 天津 300350)

    0 引言

    載荷譜[1-2]是進(jìn)行零部件疲勞壽命預(yù)測(cè)和優(yōu)化設(shè)計(jì)的重要依據(jù)。目前,載荷譜已經(jīng)在航天、汽車、高鐵等領(lǐng)域得到了廣泛研究與應(yīng)用[3-5]。而我國(guó)大型玉米收獲機(jī)關(guān)鍵零部件載荷譜的編制方法研究尚處于起步階段,與國(guó)外發(fā)達(dá)國(guó)家相比還有較大的差距。國(guó)產(chǎn)玉米收獲機(jī)故障頻發(fā),穩(wěn)定性較差,其主要原因是設(shè)計(jì)水平不高,在研發(fā)設(shè)計(jì)初期大多以仿制和借鑒為主,很多零部件的結(jié)構(gòu)設(shè)計(jì)并沒(méi)有足夠的載荷譜數(shù)據(jù)作為支撐。玉米收獲機(jī)的車架是整車的核心承載部件,承受著割臺(tái)、駕駛室、發(fā)動(dòng)機(jī)、糧倉(cāng)等傳遞的力和力矩,在實(shí)際工作過(guò)程中,經(jīng)常會(huì)出現(xiàn)車架裂痕、甚至斷裂的現(xiàn)象,嚴(yán)重影響了玉米收獲機(jī)的整車性能。因此,采用合理的載荷外推方法進(jìn)行多種工況下玉米收獲機(jī)車架的載荷譜編制,對(duì)于玉米收獲機(jī)車架的結(jié)構(gòu)優(yōu)化和整車性能提高具有重要意義。

    在實(shí)際載荷測(cè)試試驗(yàn)中,由于時(shí)間、天氣、試驗(yàn)場(chǎng)地和成本等多種原因,只能進(jìn)行一定時(shí)間內(nèi)載荷時(shí)間歷程的測(cè)量,實(shí)測(cè)載荷僅反映試驗(yàn)期間的載荷時(shí)間歷程,不能反映全壽命周期下的載荷分布,因此需要對(duì)有限的載荷時(shí)間歷程進(jìn)行合理有效的載荷外推[6]。載荷外推是影響載荷譜編制準(zhǔn)確性的關(guān)鍵環(huán)節(jié)。JOHANNESSON[7]提出了基于POT參數(shù)模型的時(shí)域外推方法,通過(guò)設(shè)置峰值的閾值,估計(jì)峰值的分布函數(shù),在最大限度保留原始載荷時(shí)間序列的基礎(chǔ)上,重構(gòu)得到外推后的載荷時(shí)間歷程。楊子涵等[8]對(duì)時(shí)域外推過(guò)程中的閾值選取方法進(jìn)行了改進(jìn),提出了一種時(shí)域外推過(guò)程中閾值選取的量化方法,解決了時(shí)域外推過(guò)程中閾值選取主觀性較強(qiáng)的問(wèn)題。張英爽等[9]將經(jīng)過(guò)雨流計(jì)數(shù)法統(tǒng)計(jì)后的載荷循環(huán)均值和載荷循環(huán)幅值的分布分別用正態(tài)分布和威布爾分布進(jìn)行參數(shù)擬合,并以此進(jìn)行載荷外推。翟新婷等[10]、GENG等[11]以混合分布函數(shù)作為函數(shù)擬合模型進(jìn)行參數(shù)估計(jì),該方法的擬合效果優(yōu)于單分布的參數(shù)估計(jì)。但是,參數(shù)外推方法是用純粹的分布函數(shù)來(lái)描述載荷的分布規(guī)律,對(duì)于復(fù)雜、隨機(jī)性較大的載荷,則會(huì)產(chǎn)生較大的誤差,大大降低了外推載荷譜與實(shí)際載荷譜的等效性。隨著研究的不斷深入,基于核密度估計(jì)法的非參數(shù)估計(jì)方法被應(yīng)用于載荷外推中,該方法既可很好地保留載荷數(shù)據(jù)本身的分布規(guī)律,又能實(shí)現(xiàn)對(duì)任意載荷分布的擬合。李凡松等[12]采用基于自適應(yīng)帶寬的核密度估計(jì)方法進(jìn)行載荷外推,每一個(gè)數(shù)據(jù)點(diǎn)都有其自身對(duì)應(yīng)的帶寬,該方法相較于固定帶寬的核密度估計(jì)方法具有更好的擬合效果。

    本文針對(duì)核密度估計(jì)載荷非參數(shù)外推方法中帶寬的選擇問(wèn)題,結(jié)合改進(jìn)的四叉樹分割算法,對(duì)核密度估計(jì)的帶寬計(jì)算進(jìn)行優(yōu)化,提出一種改進(jìn)的自適應(yīng)帶寬核密度估計(jì)的載荷外推方法。以玉米收獲機(jī)車架為研究對(duì)象,通過(guò)與傳統(tǒng)的基于核密度估計(jì)的固定帶寬外推方法和自適應(yīng)帶寬外推方法進(jìn)行比較,驗(yàn)證該自適應(yīng)帶寬核密度估計(jì)的載荷外推方法的準(zhǔn)確性和合理性。

    1 核密度估計(jì)

    核密度估計(jì)是由ROSENBLATT[13]和PARZEN[14]提出的一種由實(shí)測(cè)分布函數(shù)未知的隨機(jī)變量來(lái)估計(jì)其概率密度函數(shù)的非參數(shù)估計(jì)方法。該方法不需要知道數(shù)據(jù)的先驗(yàn)分布,也不必對(duì)數(shù)據(jù)進(jìn)行任何的假設(shè),只需要確定輸入的數(shù)據(jù)變量、核函數(shù)以及帶寬就可以估計(jì)輸入數(shù)據(jù)的概率密度函數(shù)。目前,核密度估計(jì)已經(jīng)在電力、醫(yī)療、地理等領(lǐng)域得到了非常廣泛的應(yīng)用[15-18]。玉米收獲機(jī)工作時(shí)由于其結(jié)構(gòu)復(fù)雜,工況多樣,導(dǎo)致負(fù)載波動(dòng)大,實(shí)測(cè)載荷信號(hào)具有很強(qiáng)的分散性和隨機(jī)性,相較于一般的參數(shù)估計(jì)法來(lái)說(shuō),核密度估計(jì)能更好地描述載荷數(shù)據(jù)的分布規(guī)律,從而使得載荷外推更加準(zhǔn)確。

    基于雨流計(jì)數(shù)矩陣的載荷外推是關(guān)于載荷循環(huán)均值和幅值的二維問(wèn)題,二維核密度估計(jì)表達(dá)式為

    (1)

    其中

    (2)

    式中n——輸入數(shù)據(jù)點(diǎn)的個(gè)數(shù)

    hx、hy——核密度估計(jì)載荷循環(huán)幅值、均值的帶寬

    xi、yi——輸入的第i個(gè)載荷循環(huán)幅值、均值

    K(·)——核函數(shù)

    只要選取的核函數(shù)和帶寬合適,就可以無(wú)限制地去逼近任何隨機(jī)變量真實(shí)的概率密度函數(shù)。大量研究表明,當(dāng)輸入的樣本數(shù)據(jù)量足夠大時(shí),核函數(shù)的具體形式對(duì)概率密度估計(jì)的準(zhǔn)確性產(chǎn)生的影響相對(duì)較小[19],本文選擇光滑且連續(xù)、明顯單峰分布的高斯核函數(shù)。二維高斯核函數(shù)表達(dá)式為

    (3)

    與核函數(shù)相比,帶寬h的選擇對(duì)核密度估計(jì)的準(zhǔn)確性影響更大[20]。帶寬h決定了(x,y)的光滑程度,若h較大,則有較多的數(shù)據(jù)點(diǎn)影響此處的概率密度計(jì)算,(x,y)曲線在此處較光滑,但是其與實(shí)際概率密度曲線的偏差較大;若h較小,則有較少的數(shù)據(jù)點(diǎn)影響此處的概率密度計(jì)算,(x,y)曲線在此處較陡峭,但是其與實(shí)際概率密度曲線的偏差較小。因此,為了更加準(zhǔn)確地進(jìn)行核密度估計(jì),帶寬h的選擇尤為重要。為了定量化實(shí)現(xiàn)最優(yōu)帶寬計(jì)算,提出了平均積分平方誤差(MISE),MISE表達(dá)式為

    (4)

    f(x)——實(shí)測(cè)數(shù)據(jù)真實(shí)的概率密度

    E(·)——求均值函數(shù)

    拇指法則是目前應(yīng)用最多的固定最優(yōu)帶寬的計(jì)算法則,拇指法則的最優(yōu)帶寬計(jì)算公式為

    (5)

    式中d——核密度估計(jì)的維數(shù),取2

    σ——輸入二維數(shù)據(jù)樣本的標(biāo)準(zhǔn)差

    此外,還可用無(wú)偏交叉驗(yàn)證(Unbiased cross validation)[21]、插入法(Plug-in)[22]等進(jìn)行最優(yōu)帶寬的計(jì)算。上述方法計(jì)算的最優(yōu)帶寬是固定的,即每一個(gè)數(shù)據(jù)都有著同樣的帶寬,不能自動(dòng)調(diào)節(jié)和變化,然而,在核密度估計(jì)的實(shí)際應(yīng)用中,由于數(shù)據(jù)的隨機(jī)性較強(qiáng),分布不均勻,通過(guò)固定帶寬核密度估計(jì)計(jì)算得到的概率密度可能與實(shí)際的分布相差較大,所以,希望帶寬能夠隨數(shù)據(jù)的變化而變化,在數(shù)據(jù)密集的地方取得小一些,在數(shù)據(jù)稀疏的地方取得大一些。因此,自適應(yīng)帶寬核密度估計(jì)[23-24]得到了廣泛的應(yīng)用,該方法的具體計(jì)算公式為

    (x,y)=

    (6)

    其中

    λi=(g-1(xi,yi))-α

    (7)

    (8)

    式中λi——帶寬的自適應(yīng)修訂系數(shù)

    α——敏感性參數(shù),取0~1

    一般情況下當(dāng)α取0.5時(shí),核密度估計(jì)的擬合效果較好[25],相對(duì)于固定帶寬的核密度估計(jì)來(lái)說(shuō),自適應(yīng)帶寬在一定程度上提高了核密度估計(jì)的準(zhǔn)確性,但由于數(shù)據(jù)集中的每一個(gè)數(shù)據(jù)都會(huì)影響彼此的帶寬,當(dāng)輸入數(shù)據(jù)集較大時(shí),會(huì)大大降低核密度估計(jì)的計(jì)算效率。當(dāng)一定數(shù)量的數(shù)據(jù)聚集在一個(gè)很小的區(qū)域時(shí),此區(qū)域內(nèi)帶寬的差異對(duì)核密度估計(jì)的影響十分有限,因此可以根據(jù)數(shù)據(jù)的聚集程度將同一區(qū)域內(nèi)所有數(shù)據(jù)的帶寬設(shè)置成相同的,就可以避免對(duì)這些點(diǎn)進(jìn)行不必要的詳細(xì)的帶寬計(jì)算,在保證一定核密度估計(jì)準(zhǔn)確性的前提下,提高了計(jì)算效率,所以針對(duì)核密度估計(jì)最優(yōu)帶寬的選擇問(wèn)題,為了更好兼顧核密度估計(jì)的計(jì)算效率與準(zhǔn)確性,本文在上述方法的基礎(chǔ)上進(jìn)行改進(jìn),提出了一種新的自適應(yīng)帶寬計(jì)算方法。

    2 基于四叉樹算法的自適應(yīng)帶寬核密度估計(jì)載荷外推

    基于四叉樹算法的自適應(yīng)帶寬核密度估計(jì)載荷外推主要由3個(gè)步驟組成,包括:核密度估計(jì)的數(shù)據(jù)輸入、基于四叉樹算法的數(shù)據(jù)區(qū)域分割以及基于四叉樹算法的自適應(yīng)帶寬核密度估計(jì)。具體的載荷外推流程圖如圖1所示。

    2.1 核密度估計(jì)的數(shù)據(jù)輸入

    雨流計(jì)數(shù)法得到的應(yīng)變與材料的應(yīng)力-應(yīng)變遲滯回線具有很好的一致性,并且考慮了載荷循環(huán)幅值和均值2個(gè)變量,符合疲勞載荷本身固有的特性,能夠滿足疲勞壽命預(yù)估以及載荷譜編制的條件,是目前應(yīng)用最為廣泛的實(shí)測(cè)載荷時(shí)間歷程的計(jì)數(shù)統(tǒng)計(jì)方法。首先將已經(jīng)預(yù)處理的實(shí)測(cè)載荷利用該方法進(jìn)行計(jì)數(shù)統(tǒng)計(jì)處理,得到載荷循環(huán)均幅值矩陣M(Mm、Ma),其中Mm表示載荷循環(huán)均值,Ma表示載荷循環(huán)幅值;在實(shí)際工作中,小幅值載荷循環(huán)的數(shù)量比大幅值載荷循環(huán)大,但是其造成的疲勞損傷非常小[26],所以可以將小于最大載荷循環(huán)幅值10%的載荷循環(huán)只進(jìn)行簡(jiǎn)單的線性比例外推,這很大程度上減少了核密度估計(jì)的數(shù)據(jù)輸入量,從而降低核密度估計(jì)的計(jì)算量,提高了計(jì)算效率。已經(jīng)過(guò)濾掉小幅值載荷循環(huán)的載荷循環(huán)均幅值矩陣記為Mn(Mmn、Man),并以此作為核密度估計(jì)的數(shù)據(jù)輸入。

    2.2 基于四叉樹算法的數(shù)據(jù)區(qū)域分割

    四叉樹算法[27]是一種經(jīng)典的空間分割與索引技術(shù),目前在圖像分割與空間索引領(lǐng)域已經(jīng)得到了廣泛的應(yīng)用[28-29]。該算法可以通過(guò)對(duì)數(shù)據(jù)空間進(jìn)行遞歸四等分割將數(shù)據(jù)區(qū)域分割成密度不等的塊,數(shù)據(jù)越密集的地方,塊的數(shù)量越多,塊相對(duì)較??;而數(shù)據(jù)越稀疏的地方,塊的數(shù)量越少,塊相對(duì)較大,分割出的塊的數(shù)量與大小從一定程度上反映出數(shù)據(jù)的密集程度。數(shù)據(jù)密集的地方采用較小的帶寬,能夠更好地反映出該區(qū)域分布的細(xì)節(jié),與實(shí)際分布更加接近,擬合效果較好。根據(jù)四叉樹分割算法分割出不同的數(shù)據(jù)塊,計(jì)算各自的局部最優(yōu)帶寬,達(dá)到自適應(yīng)帶寬核密度估計(jì)的目的。將一個(gè)區(qū)域進(jìn)行四叉樹分割的流程如圖2所示。

    通過(guò)對(duì)區(qū)域不斷執(zhí)行上述四叉樹分割流程,直到所有區(qū)域都不滿足區(qū)域分割的條件,最后可以將原始數(shù)據(jù)區(qū)域分割成大小不盡相同的數(shù)據(jù)塊。判斷數(shù)據(jù)區(qū)域是否需要分割的閾值,數(shù)據(jù)區(qū)域內(nèi)數(shù)據(jù)點(diǎn)個(gè)數(shù)的最大值Nmax和數(shù)據(jù)區(qū)域最小寬度Lmin對(duì)分割完成后數(shù)據(jù)塊的大小和數(shù)量有著重要的影響,不同Nmax對(duì)概率密度函數(shù)的影響如圖3所示。

    由圖3可知,Nmax較大時(shí),圖3b不能很好地反映出該區(qū)域內(nèi)數(shù)據(jù)點(diǎn)分布的細(xì)節(jié)信息;Nmax較小時(shí),圖3c反映的細(xì)節(jié)信息和圖3a沒(méi)有明顯的差別,但是Nmax越小,四叉樹分割的層次就越深,計(jì)算量越大。

    Δan=(maxMan-minMan)/64

    (9)

    Δmn=(maxMmn-minMmn)/64

    (10)

    式中Δan——載荷循環(huán)幅值每一級(jí)的長(zhǎng)度

    Δmn——載荷循環(huán)均值每一級(jí)的長(zhǎng)度

    因此定義四叉樹分割數(shù)據(jù)塊的最小寬度Lmin=min(Δan,Δmn),這既能提高四叉樹分割的計(jì)算效率,也能最小限度影響核密度估計(jì)的準(zhǔn)確性。

    2.3 基于四叉樹算法的自適應(yīng)帶寬核密度估計(jì)

    不同數(shù)據(jù)塊內(nèi)數(shù)據(jù)點(diǎn)的帶寬因?yàn)閿?shù)據(jù)的差別而有所不同,但同一數(shù)據(jù)塊內(nèi)數(shù)據(jù)點(diǎn)的帶寬是相同的。當(dāng)一個(gè)數(shù)據(jù)塊內(nèi)的數(shù)據(jù)點(diǎn)密集程度比較高,彼此相差不大時(shí),每一個(gè)數(shù)據(jù)點(diǎn)貢獻(xiàn)的核密度估計(jì)分量的差別很微小,數(shù)據(jù)落入此數(shù)據(jù)塊不同位置的概率近似相同。因此,為了進(jìn)一步提高自適應(yīng)核密度的計(jì)算效率,在數(shù)據(jù)密集程度比較高的塊,不需要每一個(gè)數(shù)據(jù)點(diǎn)都進(jìn)行核密度估計(jì)的計(jì)算,只需要將該數(shù)據(jù)區(qū)域數(shù)據(jù)點(diǎn)的平均值作為核密度估計(jì)的數(shù)據(jù)輸入,再乘以該區(qū)域數(shù)據(jù)點(diǎn)的個(gè)數(shù)即可。在此將表示數(shù)據(jù)塊內(nèi)數(shù)據(jù)點(diǎn)密集程度的參數(shù)定義為

    (11)

    其中

    Lside=min(Δan,Δmn)/2

    (12)

    H=max(hx,hy)

    (13)

    當(dāng)γ>1時(shí),表示該數(shù)據(jù)塊內(nèi)的數(shù)據(jù)點(diǎn)密集程度高,可以對(duì)此數(shù)據(jù)塊內(nèi)的核密度估計(jì)計(jì)算進(jìn)行優(yōu)化。依此改進(jìn)的自適應(yīng)核密度估計(jì)公式為

    (x,y)=

    (14)

    式中N——優(yōu)化前載荷數(shù)據(jù)點(diǎn)總個(gè)數(shù)

    N0——優(yōu)化后載荷數(shù)據(jù)點(diǎn)總個(gè)數(shù)

    當(dāng)γ>1時(shí),n(xi,yi)表示此數(shù)據(jù)塊內(nèi)數(shù)據(jù)點(diǎn)的個(gè)數(shù),xi、yi分別表示此數(shù)據(jù)塊內(nèi)數(shù)據(jù)點(diǎn)的載荷循環(huán)幅值以及載荷循環(huán)均值的平均值;當(dāng)γ<1時(shí),n(xi,yi)等于1,xi、yi分別表示此數(shù)據(jù)塊內(nèi)每一個(gè)數(shù)據(jù)點(diǎn)的載荷循環(huán)幅值以及載荷循環(huán)均值。

    通過(guò)上述方法進(jìn)行自適應(yīng)帶寬核密度估計(jì)后,再結(jié)合Monte Carlo模擬算法進(jìn)行載荷外推。

    3 實(shí)例驗(yàn)證

    3.1 試驗(yàn)方案

    玉米收獲機(jī)在實(shí)際工作中,工況較多且復(fù)雜多變,由于地形、車輛載重以及操作行為不規(guī)范等原因,車架時(shí)常會(huì)出現(xiàn)裂紋,大大降低了車架的疲勞強(qiáng)度,影響整車性能,受到很大沖擊時(shí),甚至?xí)霈F(xiàn)車架整根斷裂的問(wèn)題,嚴(yán)重影響了糧食收獲。由于試驗(yàn)環(huán)境較惡劣,并且測(cè)試部位空間有限,需要選用安裝方便并且固定牢靠的傳感器。因此本試驗(yàn)采用中航工業(yè)公司BE350-4AA型應(yīng)變片進(jìn)行數(shù)據(jù)測(cè)取,并通過(guò)有線連接的方式與HBM公司的SoMat eDAQ型數(shù)據(jù)采集儀連接。該數(shù)據(jù)采集儀具有極佳的密封性能和抗震性能,能有效應(yīng)對(duì)水濺、揚(yáng)塵、顛簸等測(cè)試環(huán)境,有線連接的方式能夠確保信號(hào)傳輸?shù)姆€(wěn)定性,通過(guò)數(shù)采的計(jì)算機(jī)終端軟件將采樣頻率設(shè)置為500 Hz。應(yīng)變片的安裝如圖5所示,現(xiàn)場(chǎng)試驗(yàn)如圖6所示,數(shù)據(jù)采集系統(tǒng)組成如圖7所示。

    砂石路面工況下行駛較為平穩(wěn),載荷變化相對(duì)較小,大幅值載荷循環(huán)較少,載荷循環(huán)的分布相對(duì)集中;田間收獲工況下載荷容易受田間地面的軟硬程度以及糧倉(cāng)的使用率等因素的影響,載荷的隨機(jī)性較高,平穩(wěn)性較差,大幅值載荷循環(huán)相對(duì)多一些,載荷循環(huán)的分布也更為分散。根據(jù)通過(guò)四叉樹分割算法得到的數(shù)據(jù)塊的大小以及密集程度,能夠很好地反映出這兩種工況下數(shù)據(jù)的密集程度以及分布情況,并且這兩種工況占玉米收獲機(jī)實(shí)際工作工況中很大的比重,具有一定的代表性,為了降低結(jié)論的偶然性,本文采用2種工況下各3組樣本進(jìn)行對(duì)比驗(yàn)證,經(jīng)過(guò)預(yù)處理的試驗(yàn)數(shù)據(jù)如圖8所示。

    3.2 結(jié)果驗(yàn)證

    通過(guò)雨流計(jì)數(shù)法對(duì)上述已經(jīng)完成預(yù)處理的車架實(shí)測(cè)原始載荷進(jìn)行統(tǒng)計(jì)計(jì)數(shù)處理,得到的部分載荷循環(huán)均幅值頻次分布圖如圖9所示。

    由圖9可知,幅值較小的載荷循環(huán)占比較大,將小于最大幅值10%的載荷循環(huán)過(guò)濾掉,過(guò)濾完成的數(shù)據(jù)作為固定帶寬和自適應(yīng)帶寬核密度估計(jì)的數(shù)據(jù)輸入,本文方法的核密度估計(jì)的數(shù)據(jù)輸入在過(guò)濾的基礎(chǔ)上還需要通過(guò)改進(jìn)的四叉樹算法進(jìn)行優(yōu)化。首先將已經(jīng)篩選好的載荷循環(huán)均幅值矩陣中的數(shù)據(jù)點(diǎn)繪制成如圖10所示的散點(diǎn)圖,并在此散點(diǎn)圖上利用改進(jìn)的四叉樹算法進(jìn)行數(shù)據(jù)區(qū)域的分割操作,數(shù)據(jù)區(qū)域信息如表1所示。

    表1 數(shù)據(jù)區(qū)域信息Tab.1 Information of data area

    經(jīng)過(guò)基于四叉樹算法的數(shù)據(jù)區(qū)域分割的計(jì)算,數(shù)據(jù)區(qū)域被分割成大小不相同的數(shù)據(jù)塊,每個(gè)數(shù)據(jù)塊都有各自的帶寬,利用式(14)對(duì)數(shù)據(jù)塊的自適應(yīng)帶寬核密度估計(jì)進(jìn)行優(yōu)化,在之前濾除小幅值載荷循環(huán)的基礎(chǔ)上又大幅減少了核密度估計(jì)的輸入,進(jìn)一步提高了核密度估計(jì)的計(jì)算效率。

    為了驗(yàn)證本文提出的基于四叉樹算法的自適應(yīng)帶寬核密度估計(jì)的載荷外推方法的準(zhǔn)確性與合理性,與傳統(tǒng)的基于核密度估計(jì)的固定帶寬載荷外推和自適應(yīng)帶寬載荷外推以及實(shí)測(cè)原始載荷進(jìn)行對(duì)比論證。

    由3種帶寬選擇方法計(jì)算得到的載荷循環(huán)均幅值概率密度分布如圖11所示。將圖11與圖10進(jìn)行對(duì)比,可以看出傳統(tǒng)的固定帶寬與自適應(yīng)帶寬核密度估計(jì)的曲線較為平滑,而本文提出的自適應(yīng)帶寬核密度估計(jì)的曲線較為曲折,能夠很好地反映圖10數(shù)據(jù)的密集程度,并且數(shù)據(jù)密集點(diǎn)的分布情況劃分更加具體,細(xì)節(jié)信息表達(dá)更為清楚,能很好地?cái)M合真實(shí)概率密度。為了進(jìn)一步評(píng)價(jià)概率密度函數(shù)的擬合程度,將固定帶寬、自適應(yīng)帶寬以及本文提出的自適應(yīng)帶寬核密度估計(jì)的概率密度函數(shù)的結(jié)果分別結(jié)合Monte Carlo模擬算法進(jìn)行載荷外推。將實(shí)測(cè)原始載荷數(shù)據(jù)和上述3種方法進(jìn)行載荷外推后得到的載荷循環(huán)均幅值分別進(jìn)行分級(jí)統(tǒng)計(jì)處理,并計(jì)算實(shí)測(cè)原始數(shù)據(jù)與3種帶寬選擇外推方法的均方根誤差(RMSE),計(jì)算結(jié)果如表2所示。

    表2 載荷循環(huán)均幅值的均方根誤差Tab.2 RMSE of load cycle mean and range με

    由表2可知,相較于傳統(tǒng)的固定帶寬和自適應(yīng)帶寬核密度估計(jì),本文提出的自適應(yīng)帶寬核密度估計(jì)載荷外推方法的均方根誤差RMSE最小,表明通過(guò)本文方法計(jì)算得到的概率密度更加接近于真實(shí)的概率密度。

    將實(shí)測(cè)原始載荷數(shù)據(jù)和上述3種方法進(jìn)行載荷外推后得到的載荷循環(huán)均值和幅值數(shù)據(jù)分別進(jìn)行分級(jí)統(tǒng)計(jì)處理,結(jié)果如圖12所示。為進(jìn)一步分析頻次分布的相關(guān)性和擬合程度,計(jì)算實(shí)測(cè)原始數(shù)據(jù)與3種帶寬選擇外推方法的相關(guān)系數(shù)(R)以及均方根誤差(RMSE),結(jié)果如表3、4所示。

    表3 載荷循環(huán)均值的相關(guān)系數(shù)及均方根誤差Tab.3 R and RMSE of load cycle mean

    表4 載荷循環(huán)幅值的相關(guān)系數(shù)及均方根誤差Tab.4 R and RMSE of load cycle range

    由表3、4可知,相較于傳統(tǒng)的固定帶寬和自適應(yīng)帶寬核密度估計(jì),本文提出的自適應(yīng)帶寬核密度估計(jì)載荷外推與實(shí)測(cè)原始載荷均值和幅值的相關(guān)系數(shù)R最大,均方根誤差最小,表明通過(guò)本文提出的自適應(yīng)帶寬核密度估計(jì)進(jìn)行載荷外推得到的載荷循環(huán)均值和幅值分布與實(shí)測(cè)原始載荷的載荷循環(huán)均值和幅值分布具有很強(qiáng)的相似性,分布擬合效果好,更加接近實(shí)測(cè)原始載荷的真實(shí)分布規(guī)律。

    載荷循環(huán)幅值累積頻次曲線為后續(xù)程序載荷譜的編制提供數(shù)據(jù)支撐,影響程序載荷譜編制的準(zhǔn)確性。將實(shí)測(cè)原始數(shù)據(jù)與傳統(tǒng)的固定帶寬載荷外推、自適應(yīng)帶寬載荷外推以及本文提出的自適應(yīng)帶寬載荷外推的載荷循環(huán)幅值累積頻次曲線進(jìn)行對(duì)比驗(yàn)證,如圖13所示。由圖13可以看出,通過(guò)本文提出的帶寬選擇方法進(jìn)行載荷外推得到的曲線與實(shí)測(cè)原始載荷曲線更為接近,擬合效果更好。為了更加直觀進(jìn)行對(duì)比,更好評(píng)價(jià)載荷循環(huán)幅值累積頻次曲線的擬合效果,采用決定系數(shù)R2作為曲線擬合的檢驗(yàn)指標(biāo),計(jì)算結(jié)果如表5所示。

    表5 載荷循環(huán)幅值累積頻次曲線的決定系數(shù)Tab.5 Coefficient of determination of load cycle range cumulative frequency curve

    由表5可知,本文提出的自適應(yīng)帶寬核密度估計(jì)外推方法的R2更加接近于1,這表明,與另外兩種方法相比,本文提出的方法與實(shí)測(cè)原始載荷的載荷循環(huán)幅值累積頻次曲線的擬合度更高,更為相似,能夠?yàn)槌绦蜉d荷譜編制提供更加接近實(shí)際的數(shù)據(jù)支撐。

    4 結(jié)論

    (1)提出了一種基于四叉樹算法的自適應(yīng)帶寬核密度估計(jì)載荷外推方法,該方法利用四叉樹分割算法將載荷循環(huán)數(shù)據(jù)區(qū)域進(jìn)行分割,以分割完成后的數(shù)據(jù)塊為依據(jù)計(jì)算帶寬,并根據(jù)數(shù)據(jù)塊的密集程度對(duì)核密度估計(jì)的數(shù)據(jù)輸入進(jìn)行優(yōu)化。

    (2)以玉米收獲機(jī)車架為研究對(duì)象,在砂石路面工況和田間收獲工況下,采用本文提出的方法大幅降低了核密度估計(jì)的數(shù)據(jù)輸入,極大提高了核密度估計(jì)的計(jì)算效率。

    (3)將固定帶寬、自適應(yīng)帶寬與本文提出的基于四叉樹算法的自適應(yīng)帶寬核密度估計(jì)載荷外推方法進(jìn)行了對(duì)比驗(yàn)證。結(jié)果表明,通過(guò)本文提出的載荷外推方法計(jì)算得到的載荷循環(huán)均值和幅值分布的相關(guān)系數(shù)均更加接近于1,均方根誤差更小;載荷循環(huán)幅值累積頻次曲線的決定系數(shù)均大于0.99。與傳統(tǒng)的固定帶寬以及自適應(yīng)帶寬的核密度估計(jì)方法相比,本文提出的方法具有更高的準(zhǔn)確性,載荷外推結(jié)果更接近實(shí)際載荷的分布情況,能夠?yàn)檩d荷譜的編制、零部件結(jié)構(gòu)優(yōu)化提供參考。

    猜你喜歡
    四叉樹密度估計(jì)概率密度
    中國(guó)人均可支配收入的空間區(qū)域動(dòng)態(tài)演變與差異分析
    m-NOD樣本最近鄰密度估計(jì)的相合性
    面向魚眼圖像的人群密度估計(jì)
    基于MATLAB 的核密度估計(jì)研究
    科技視界(2021年4期)2021-04-13 06:03:56
    連續(xù)型隨機(jī)變量函數(shù)的概率密度公式
    基于WebGL的三維點(diǎn)云可視化研究
    基于四叉樹的高效梯度域圖像融合
    Hunt過(guò)程在Girsanov變換下的轉(zhuǎn)移概率密度的表示公式
    隨機(jī)變量線性組合的分布的一個(gè)算法
    隨機(jī)結(jié)構(gòu)-TMD優(yōu)化設(shè)計(jì)與概率密度演化研究
    日韩大尺度精品在线看网址| 桃红色精品国产亚洲av| 草草在线视频免费看| 人人妻,人人澡人人爽秒播| 亚洲中文字幕日韩| 亚洲无线观看免费| 最近最新中文字幕大全电影3| 久久草成人影院| x7x7x7水蜜桃| 99国产精品一区二区蜜桃av| 成年女人永久免费观看视频| 国产精品国产三级国产av玫瑰| 午夜影院日韩av| 精品人妻偷拍中文字幕| 人妻制服诱惑在线中文字幕| 日韩一区二区视频免费看| 尤物成人国产欧美一区二区三区| 嫁个100分男人电影在线观看| 天美传媒精品一区二区| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久伊人网av| eeuss影院久久| 亚洲乱码一区二区免费版| 日韩精品青青久久久久久| 国产视频一区二区在线看| 特级一级黄色大片| 真人做人爱边吃奶动态| 免费在线观看日本一区| 精品久久久久久久久久免费视频| 色尼玛亚洲综合影院| 亚洲不卡免费看| 成人亚洲精品av一区二区| 欧美色欧美亚洲另类二区| 美女 人体艺术 gogo| 最近在线观看免费完整版| 69人妻影院| 亚洲av成人av| 我的女老师完整版在线观看| 少妇被粗大猛烈的视频| 国产午夜精品久久久久久一区二区三区 | 少妇的逼水好多| 国产在线精品亚洲第一网站| 国内久久婷婷六月综合欲色啪| 成年女人毛片免费观看观看9| 久久精品国产99精品国产亚洲性色| 无人区码免费观看不卡| 可以在线观看的亚洲视频| 国产精品嫩草影院av在线观看 | 中文字幕精品亚洲无线码一区| 免费在线观看成人毛片| 久9热在线精品视频| 午夜免费激情av| 男插女下体视频免费在线播放| 亚洲精品久久国产高清桃花| 欧美成人a在线观看| 欧美黑人巨大hd| 嫩草影院新地址| 男女之事视频高清在线观看| 搡老妇女老女人老熟妇| 搡老妇女老女人老熟妇| 欧美在线一区亚洲| 日韩欧美免费精品| 18禁在线播放成人免费| 在线观看美女被高潮喷水网站| 搡老妇女老女人老熟妇| 两个人视频免费观看高清| 久久6这里有精品| 两个人视频免费观看高清| 亚洲国产精品成人综合色| 亚洲av第一区精品v没综合| 国产黄色小视频在线观看| 日韩欧美免费精品| 精品免费久久久久久久清纯| 一区二区三区免费毛片| 最近在线观看免费完整版| 亚洲国产色片| 国产真实乱freesex| 少妇人妻精品综合一区二区 | 国产又黄又爽又无遮挡在线| 午夜精品一区二区三区免费看| 少妇的逼水好多| 韩国av一区二区三区四区| 亚洲国产日韩欧美精品在线观看| 尾随美女入室| 欧美人与善性xxx| 欧美+日韩+精品| 一个人看的www免费观看视频| av天堂中文字幕网| 精品一区二区三区人妻视频| 一进一出抽搐gif免费好疼| 22中文网久久字幕| 国产乱人视频| 国产三级中文精品| 黄片wwwwww| 91麻豆精品激情在线观看国产| 免费不卡的大黄色大毛片视频在线观看 | 极品教师在线免费播放| 国产一区二区激情短视频| 俄罗斯特黄特色一大片| 亚洲黑人精品在线| 久9热在线精品视频| 极品教师在线视频| 成人亚洲精品av一区二区| 色尼玛亚洲综合影院| 精品久久久久久久久久免费视频| 国产激情偷乱视频一区二区| 亚洲,欧美,日韩| 夜夜看夜夜爽夜夜摸| av.在线天堂| 日日啪夜夜撸| 国产真实伦视频高清在线观看 | 午夜爱爱视频在线播放| 成年免费大片在线观看| 蜜桃久久精品国产亚洲av| 亚洲国产精品久久男人天堂| 麻豆国产97在线/欧美| 久久久久久久精品吃奶| 三级毛片av免费| 男女之事视频高清在线观看| 国产精品人妻久久久久久| 国产又黄又爽又无遮挡在线| 久久国产精品人妻蜜桃| 精品久久久噜噜| 国产精品99久久久久久久久| 日本 av在线| 精品人妻视频免费看| 久久精品国产清高在天天线| 床上黄色一级片| 99国产精品一区二区蜜桃av| 18禁黄网站禁片午夜丰满| 国产v大片淫在线免费观看| 九九爱精品视频在线观看| 国产综合懂色| 少妇的逼水好多| 精品福利观看| 亚洲中文字幕一区二区三区有码在线看| 91午夜精品亚洲一区二区三区 | 国产真实乱freesex| 蜜桃亚洲精品一区二区三区| 国产aⅴ精品一区二区三区波| 国产中年淑女户外野战色| 成人二区视频| 精品一区二区三区av网在线观看| 亚洲成av人片在线播放无| 国产亚洲精品久久久久久毛片| 精品人妻1区二区| 精品国产三级普通话版| 变态另类丝袜制服| 日本爱情动作片www.在线观看 | 性插视频无遮挡在线免费观看| 久久久久九九精品影院| aaaaa片日本免费| 黄色欧美视频在线观看| 国产精品亚洲美女久久久| 国内精品久久久久精免费| 伦精品一区二区三区| 欧美三级亚洲精品| 亚洲成人精品中文字幕电影| 国产精品一区二区三区四区久久| 精品免费久久久久久久清纯| 欧美色视频一区免费| 一卡2卡三卡四卡精品乱码亚洲| 国产av麻豆久久久久久久| 亚洲专区国产一区二区| 2021天堂中文幕一二区在线观| 日韩欧美三级三区| 校园人妻丝袜中文字幕| 欧美色视频一区免费| 亚洲久久久久久中文字幕| 精品久久久久久久人妻蜜臀av| 波野结衣二区三区在线| 国产三级在线视频| 亚洲精品色激情综合| 国产精品久久久久久av不卡| 久久久国产成人精品二区| 亚洲午夜理论影院| 成年版毛片免费区| 97人妻精品一区二区三区麻豆| 亚洲av.av天堂| 舔av片在线| 好男人在线观看高清免费视频| 亚洲avbb在线观看| 中国美白少妇内射xxxbb| 国产爱豆传媒在线观看| 亚洲av美国av| 欧美黑人巨大hd| 欧美3d第一页| 中文字幕av在线有码专区| 永久网站在线| 黄色女人牲交| 最近中文字幕高清免费大全6 | 精品国内亚洲2022精品成人| 亚洲无线观看免费| 国产乱人视频| 十八禁网站免费在线| a在线观看视频网站| 能在线免费观看的黄片| 精品国内亚洲2022精品成人| 国产精品一区二区性色av| 国产主播在线观看一区二区| 成人午夜高清在线视频| 少妇人妻精品综合一区二区 | 国产探花在线观看一区二区| 国产黄片美女视频| 美女大奶头视频| 淫秽高清视频在线观看| 在线观看午夜福利视频| 老司机午夜福利在线观看视频| 国产欧美日韩精品一区二区| 国产精品国产高清国产av| 一区二区三区高清视频在线| 在线观看一区二区三区| 极品教师在线免费播放| 亚洲四区av| 观看免费一级毛片| 亚洲色图av天堂| 99国产极品粉嫩在线观看| 国国产精品蜜臀av免费| 日韩欧美精品免费久久| 亚洲久久久久久中文字幕| 成人永久免费在线观看视频| 99久久久亚洲精品蜜臀av| 在线观看一区二区三区| 国产一区二区在线av高清观看| 听说在线观看完整版免费高清| 麻豆国产av国片精品| 乱码一卡2卡4卡精品| 在线a可以看的网站| 欧美bdsm另类| 男人和女人高潮做爰伦理| 亚洲人与动物交配视频| 中文字幕久久专区| 亚洲中文字幕日韩| 99热只有精品国产| 中国美女看黄片| 内地一区二区视频在线| 人人妻人人看人人澡| 国产一区二区三区视频了| 亚洲一级一片aⅴ在线观看| 久久久久国内视频| 亚洲第一电影网av| 波野结衣二区三区在线| 久久婷婷人人爽人人干人人爱| 亚洲黑人精品在线| 最新在线观看一区二区三区| 免费一级毛片在线播放高清视频| ponron亚洲| 国产国拍精品亚洲av在线观看| 久9热在线精品视频| 99热精品在线国产| 99视频精品全部免费 在线| 久久精品夜夜夜夜夜久久蜜豆| 日日干狠狠操夜夜爽| 精品久久久久久成人av| 伦理电影大哥的女人| 久久婷婷人人爽人人干人人爱| 在线观看午夜福利视频| 在线播放国产精品三级| 中文亚洲av片在线观看爽| 色吧在线观看| 中国美女看黄片| 级片在线观看| 亚洲第一区二区三区不卡| 亚洲自偷自拍三级| 简卡轻食公司| 真人做人爱边吃奶动态| 男女啪啪激烈高潮av片| 国产国拍精品亚洲av在线观看| 高清毛片免费观看视频网站| 又爽又黄a免费视频| 老熟妇乱子伦视频在线观看| 亚洲欧美日韩无卡精品| 我要搜黄色片| 夜夜夜夜夜久久久久| 在线观看舔阴道视频| 少妇人妻一区二区三区视频| 俄罗斯特黄特色一大片| 亚洲成a人片在线一区二区| 成人欧美大片| 亚洲精品亚洲一区二区| 国产一区二区三区在线臀色熟女| 韩国av在线不卡| 熟妇人妻久久中文字幕3abv| 无遮挡黄片免费观看| 我要搜黄色片| 又爽又黄无遮挡网站| 欧美日韩瑟瑟在线播放| 午夜福利欧美成人| 毛片一级片免费看久久久久 | 三级男女做爰猛烈吃奶摸视频| 九九久久精品国产亚洲av麻豆| 舔av片在线| 性欧美人与动物交配| 国产免费男女视频| 18禁黄网站禁片免费观看直播| 制服丝袜大香蕉在线| 国产单亲对白刺激| 成人高潮视频无遮挡免费网站| 黄片wwwwww| 一本久久中文字幕| 亚洲av成人精品一区久久| 一个人看视频在线观看www免费| 91av网一区二区| 久久亚洲真实| 最近最新中文字幕大全电影3| 全区人妻精品视频| 天堂动漫精品| 一进一出抽搐动态| 欧美高清成人免费视频www| 国产精品福利在线免费观看| 两个人视频免费观看高清| 久久精品国产亚洲av天美| 国产高清不卡午夜福利| 变态另类丝袜制服| 麻豆精品久久久久久蜜桃| 欧美日本视频| 久久精品国产自在天天线| 成人无遮挡网站| 久久久久久久亚洲中文字幕| 欧美三级亚洲精品| 欧美日韩中文字幕国产精品一区二区三区| av女优亚洲男人天堂| 精品一区二区三区av网在线观看| 国产午夜福利久久久久久| 哪里可以看免费的av片| 国产日本99.免费观看| 中文资源天堂在线| 亚洲自拍偷在线| 亚洲精华国产精华精| 精品99又大又爽又粗少妇毛片 | 有码 亚洲区| 色视频www国产| 国产白丝娇喘喷水9色精品| 精品一区二区免费观看| 久久国内精品自在自线图片| 精品国内亚洲2022精品成人| 久久精品人妻少妇| 欧美性猛交╳xxx乱大交人| 欧美成人免费av一区二区三区| 国产精品国产高清国产av| 久久久久久久久久黄片| 在现免费观看毛片| 中文字幕高清在线视频| 日韩欧美精品v在线| 亚洲精品久久国产高清桃花| 日日摸夜夜添夜夜添小说| av国产免费在线观看| 亚洲五月天丁香| 亚洲三级黄色毛片| 国产男靠女视频免费网站| 两个人视频免费观看高清| 天堂影院成人在线观看| 麻豆久久精品国产亚洲av| 日韩强制内射视频| 如何舔出高潮| 国产一区二区三区视频了| 91久久精品国产一区二区三区| 国产精品久久电影中文字幕| 九九在线视频观看精品| a级毛片免费高清观看在线播放| 国产单亲对白刺激| 精品不卡国产一区二区三区| 十八禁网站免费在线| xxxwww97欧美| 成年版毛片免费区| 国产大屁股一区二区在线视频| 老司机深夜福利视频在线观看| 搡女人真爽免费视频火全软件 | 有码 亚洲区| 女人被狂操c到高潮| 九九在线视频观看精品| 亚洲精品国产成人久久av| 欧美日韩亚洲国产一区二区在线观看| 熟女电影av网| a级毛片免费高清观看在线播放| 久久99热6这里只有精品| 免费无遮挡裸体视频| 国产私拍福利视频在线观看| 久久久国产成人免费| 欧美人与善性xxx| 国产精品一及| 熟女人妻精品中文字幕| 中国美白少妇内射xxxbb| 国产精品久久久久久精品电影| 乱系列少妇在线播放| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 99热这里只有是精品在线观看| 波野结衣二区三区在线| 五月玫瑰六月丁香| 成人二区视频| 日本-黄色视频高清免费观看| 99九九线精品视频在线观看视频| 亚洲avbb在线观看| 色哟哟哟哟哟哟| 欧美zozozo另类| 成人一区二区视频在线观看| 精品人妻1区二区| 久久午夜福利片| 有码 亚洲区| 51国产日韩欧美| 少妇高潮的动态图| 色综合色国产| 搡老岳熟女国产| 嫁个100分男人电影在线观看| netflix在线观看网站| 一a级毛片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产精品嫩草影院av在线观看 | 变态另类成人亚洲欧美熟女| 免费看a级黄色片| 成人特级av手机在线观看| x7x7x7水蜜桃| 精品不卡国产一区二区三区| 琪琪午夜伦伦电影理论片6080| 亚洲avbb在线观看| 亚洲av中文av极速乱 | 老司机深夜福利视频在线观看| 国产成年人精品一区二区| 欧美3d第一页| 1000部很黄的大片| 欧美激情在线99| 人人妻人人看人人澡| 一区福利在线观看| 三级毛片av免费| 99视频精品全部免费 在线| 一级毛片久久久久久久久女| 精品久久久久久久久久免费视频| 97人妻精品一区二区三区麻豆| 成人av在线播放网站| 老熟妇乱子伦视频在线观看| 女生性感内裤真人,穿戴方法视频| 直男gayav资源| 国产三级中文精品| 97热精品久久久久久| 在线观看舔阴道视频| 国产男人的电影天堂91| 男人狂女人下面高潮的视频| 我的老师免费观看完整版| 一区二区三区四区激情视频 | 嫩草影院精品99| 淫妇啪啪啪对白视频| 3wmmmm亚洲av在线观看| 很黄的视频免费| 亚洲中文字幕一区二区三区有码在线看| av视频在线观看入口| 国产亚洲精品综合一区在线观看| 国产精品,欧美在线| 午夜福利在线在线| 免费人成在线观看视频色| 深夜a级毛片| av女优亚洲男人天堂| 欧美区成人在线视频| 搡老妇女老女人老熟妇| 精品人妻一区二区三区麻豆 | 日日摸夜夜添夜夜添小说| 亚洲乱码一区二区免费版| 禁无遮挡网站| 国产精品98久久久久久宅男小说| 最近视频中文字幕2019在线8| 国产精品一区二区免费欧美| 久久精品久久久久久噜噜老黄 | 午夜久久久久精精品| 久久久午夜欧美精品| 亚洲va在线va天堂va国产| 搞女人的毛片| 九色成人免费人妻av| 成人二区视频| 日韩一本色道免费dvd| www.www免费av| 国产乱人视频| 亚洲精品影视一区二区三区av| 国产在线男女| 美女高潮喷水抽搐中文字幕| 亚洲av不卡在线观看| 人人妻人人澡欧美一区二区| 日韩亚洲欧美综合| av.在线天堂| 性插视频无遮挡在线免费观看| 人妻丰满熟妇av一区二区三区| 黄色女人牲交| 老司机福利观看| 日韩欧美国产一区二区入口| 午夜精品在线福利| 亚洲精品一卡2卡三卡4卡5卡| 欧美性感艳星| 日韩一区二区视频免费看| 99久久中文字幕三级久久日本| 女的被弄到高潮叫床怎么办 | 久久久国产成人免费| 国产视频一区二区在线看| а√天堂www在线а√下载| 国产亚洲精品综合一区在线观看| 高清日韩中文字幕在线| 亚洲精品乱码久久久v下载方式| 国产一区二区在线观看日韩| 麻豆精品久久久久久蜜桃| 又粗又爽又猛毛片免费看| 国产精品嫩草影院av在线观看 | 丰满的人妻完整版| 九九爱精品视频在线观看| 欧美日韩精品成人综合77777| 午夜激情福利司机影院| 91久久精品电影网| 亚洲精品久久国产高清桃花| 欧美精品啪啪一区二区三区| 18禁在线播放成人免费| 欧美高清成人免费视频www| 九色成人免费人妻av| 国产在线精品亚洲第一网站| 久久久成人免费电影| 亚洲久久久久久中文字幕| 美女大奶头视频| 国产精品久久久久久久久免| 人人妻人人澡欧美一区二区| 国产中年淑女户外野战色| 老熟妇仑乱视频hdxx| 久9热在线精品视频| av.在线天堂| 午夜免费激情av| 亚洲精品456在线播放app | 久久久国产成人免费| 国产极品精品免费视频能看的| 天天躁日日操中文字幕| 看免费成人av毛片| 欧美一区二区亚洲| 国产午夜精品论理片| 久久久精品欧美日韩精品| 免费一级毛片在线播放高清视频| 久久精品影院6| 亚洲欧美日韩高清在线视频| 久久中文看片网| 亚洲天堂国产精品一区在线| 他把我摸到了高潮在线观看| 国产久久久一区二区三区| 欧美bdsm另类| 婷婷丁香在线五月| 91av网一区二区| 特大巨黑吊av在线直播| 波多野结衣高清无吗| 国产麻豆成人av免费视频| 22中文网久久字幕| 可以在线观看的亚洲视频| 国产精品综合久久久久久久免费| 麻豆国产av国片精品| 日本黄大片高清| 欧美色视频一区免费| 亚洲成人免费电影在线观看| a在线观看视频网站| 熟女电影av网| 国产精华一区二区三区| 熟女电影av网| 日韩高清综合在线| 国产精品久久久久久久电影| 久久九九热精品免费| 国产极品精品免费视频能看的| 久久久久九九精品影院| 国产熟女欧美一区二区| 午夜免费激情av| 在线天堂最新版资源| 午夜精品在线福利| 22中文网久久字幕| 日本撒尿小便嘘嘘汇集6| 欧美日本亚洲视频在线播放| 在线观看午夜福利视频| 精品不卡国产一区二区三区| 18禁裸乳无遮挡免费网站照片| 亚洲电影在线观看av| 真人一进一出gif抽搐免费| 国产午夜福利久久久久久| 97热精品久久久久久| 午夜福利在线观看免费完整高清在 | 淫秽高清视频在线观看| 波多野结衣巨乳人妻| 国产精品免费一区二区三区在线| 国产精品永久免费网站| 男女做爰动态图高潮gif福利片| 人妻久久中文字幕网| 国产一区二区三区av在线 | 最新中文字幕久久久久| 久久久久久久久久久丰满 | 国产成人影院久久av| 精品一区二区三区视频在线观看免费| 在现免费观看毛片| 国产女主播在线喷水免费视频网站 | 欧美激情久久久久久爽电影| 69人妻影院| 99久久九九国产精品国产免费| 欧美xxxx黑人xx丫x性爽| 日本黄大片高清| 国产真实乱freesex| 欧美又色又爽又黄视频| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 国内久久婷婷六月综合欲色啪| 国产精品国产三级国产av玫瑰| 99精品在免费线老司机午夜| 色综合站精品国产| 成人永久免费在线观看视频| 狠狠狠狠99中文字幕| 婷婷亚洲欧美| 最后的刺客免费高清国语| 亚洲成人久久性| 欧美色欧美亚洲另类二区| 精华霜和精华液先用哪个| 精品99又大又爽又粗少妇毛片 | 内地一区二区视频在线| 我的女老师完整版在线观看| 99精品在免费线老司机午夜| 最近视频中文字幕2019在线8| 国产探花极品一区二区| 午夜福利在线观看吧| 不卡一级毛片| 亚洲欧美激情综合另类| 国产一区二区三区av在线 |