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

    用戶關(guān)聯(lián)的驅(qū)動(dòng)橋試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范研究

    2022-08-05 06:23:36鄒喜紅茍林林蔣明聰袁冬梅
    中國(guó)機(jī)械工程 2022年14期
    關(guān)鍵詞:試驗(yàn)場(chǎng)路況里程

    鄒喜紅 凌 龍 陳 靜 王 超 茍林林 蔣明聰 袁冬梅

    重慶理工大學(xué)汽車(chē)零部件先進(jìn)制造技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,重慶,400054

    0 引言

    在商用車(chē)領(lǐng)域中,發(fā)動(dòng)機(jī)、驅(qū)動(dòng)橋和駕駛室一起構(gòu)成了商用車(chē)的核心三部件,其重要性不言而喻[1],因此需在新車(chē)上市前在試驗(yàn)場(chǎng)進(jìn)行道路試驗(yàn),作為考核和驗(yàn)證動(dòng)力傳動(dòng)系統(tǒng)耐久性與可靠性的最終環(huán)節(jié)[2]。國(guó)際上先進(jìn)的車(chē)企如奔馳、寶馬、豐田、通用、福特等都建立了自己的汽車(chē)試驗(yàn)場(chǎng)和試驗(yàn)規(guī)范并開(kāi)展了整車(chē)及其關(guān)鍵零部件的耐久性評(píng)價(jià)。而我國(guó)許多車(chē)企現(xiàn)有的試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范大多是參考國(guó)外的試驗(yàn)場(chǎng)規(guī)范,依賴過(guò)去的經(jīng)驗(yàn)和習(xí)慣,以能滿足最差的路況而建立的,難以真正代表我國(guó)的用戶道路,以至于汽車(chē)零部件在試驗(yàn)場(chǎng)的失效模式與用戶實(shí)際使用情況存在一定的差異。

    近年來(lái),為研究試驗(yàn)場(chǎng)與用戶道路載荷之間的當(dāng)量關(guān)系,已經(jīng)形成了相對(duì)固定的工作流程[3]。文獻(xiàn)[4-6]針對(duì)汽車(chē)承載件將車(chē)輪軸頭加速度及懸架應(yīng)變信號(hào)作為關(guān)聯(lián)通道,分別比較了試驗(yàn)場(chǎng)與用戶載荷循環(huán)的差異并構(gòu)建了關(guān)聯(lián)模型,從而制定出與用戶關(guān)聯(lián)的試驗(yàn)場(chǎng)規(guī)范。而引起傳動(dòng)系統(tǒng)零部件損傷的直接因素是力矩,與加速度等信號(hào)無(wú)必然聯(lián)系,依據(jù)這類(lèi)信號(hào)的關(guān)聯(lián)只考慮到汽車(chē)承載件,而忽略了對(duì)動(dòng)力傳動(dòng)系統(tǒng)零部件的耐久性考核。為此,文獻(xiàn)[7-8]以兩半軸輸出端扭矩為研究對(duì)象,結(jié)合擋位及轉(zhuǎn)速信號(hào),依據(jù)市場(chǎng)調(diào)研結(jié)果,采集用戶道路載荷譜,編制臺(tái)架試驗(yàn)所需程序載荷譜,從而完成對(duì)動(dòng)力傳動(dòng)系統(tǒng)的疲勞耐久性考核。文獻(xiàn)[9]利用不同工況、不同擋位下兩傳動(dòng)軸扭矩的偽損傷值來(lái)構(gòu)造目標(biāo)矩陣及優(yōu)化函數(shù),采用粒子群優(yōu)化算法對(duì)動(dòng)力傳動(dòng)系統(tǒng)的兩試驗(yàn)場(chǎng)關(guān)聯(lián)進(jìn)行計(jì)算求解,最終得到關(guān)聯(lián)兩個(gè)試驗(yàn)場(chǎng)的汽車(chē)動(dòng)力傳動(dòng)系統(tǒng)耐久性規(guī)范。

    總結(jié)目前現(xiàn)有文獻(xiàn)可知,大部分文獻(xiàn)的研究熱點(diǎn)集中在整車(chē)承載件試驗(yàn)場(chǎng)耐久性的用戶關(guān)聯(lián),對(duì)傳動(dòng)系統(tǒng)的耐久性考核更多的是利用用戶或試驗(yàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)等效到臺(tái)架試驗(yàn)上,且對(duì)于驅(qū)動(dòng)橋而言,只考慮了傳動(dòng)軸載荷單一損傷,并沒(méi)有進(jìn)一步關(guān)聯(lián)齒輪等旋轉(zhuǎn)件。其次,用戶與試驗(yàn)場(chǎng)關(guān)聯(lián)模型多為超定方程組,利用最小二乘法進(jìn)行求解時(shí),由于解的分散性較高,常常會(huì)出現(xiàn)零解的情況,因此不利于試驗(yàn)規(guī)范的制定[10]。

    綜上,本文以采集某商用車(chē)用戶和試驗(yàn)場(chǎng)載荷為基礎(chǔ),計(jì)算軸載荷與齒載荷偽損傷值,建立“用戶-試驗(yàn)場(chǎng)”等損傷關(guān)聯(lián)模型。結(jié)合帶精英策略的非支配排序遺傳算法(NSGA-Ⅱ)求解多目標(biāo)優(yōu)化模型并選取出最優(yōu)解,從而實(shí)現(xiàn)載荷損傷和分布的等效,并制定出驅(qū)動(dòng)橋試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范。

    1 用戶關(guān)聯(lián)技術(shù)

    通過(guò)對(duì)以用戶關(guān)聯(lián)技術(shù)為代表的道路載荷大數(shù)據(jù)進(jìn)行采集和統(tǒng)計(jì)分析,可以建立與特定用戶群體的駕駛習(xí)慣、路面特征等相關(guān)的且對(duì)應(yīng)于一定設(shè)計(jì)里程的道路載荷分布模型。這一分布模型對(duì)于整車(chē)正向耐久性工程而言是戰(zhàn)略性的頂層輸入[11]。獲得這一分布模型后,如何將其等效到強(qiáng)化路面上制定出試驗(yàn)場(chǎng)強(qiáng)化路面試驗(yàn)規(guī)范,從而惠及幾代同類(lèi)車(chē)型的開(kāi)發(fā),是接踵而來(lái)的一個(gè)重要問(wèn)題。圖1所示為試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范的試驗(yàn)流程。首先通過(guò)特定用戶群體路面調(diào)查和試驗(yàn)場(chǎng)工況選取來(lái)采集實(shí)測(cè)載荷;其次以載荷統(tǒng)計(jì)結(jié)果為依據(jù)進(jìn)行載荷合成及外推,從而獲得全壽命周期內(nèi)的用戶載荷分布;然后根據(jù)Miner線性疲勞損傷原理,計(jì)算用戶和試驗(yàn)場(chǎng)載荷的偽損傷值,以用戶載荷損傷為目標(biāo)損傷矩陣T、以試驗(yàn)場(chǎng)載荷損傷為系數(shù)損傷矩陣D,建立等損傷關(guān)聯(lián)模型;最后,通過(guò)優(yōu)化算法求解出各試驗(yàn)場(chǎng)工況最優(yōu)循環(huán)次數(shù)β,并從載荷損傷及分布兩個(gè)角度驗(yàn)證模型解的有效性,從而完成對(duì)關(guān)聯(lián)用戶的試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范的制定。

    圖1 用戶關(guān)聯(lián)的試驗(yàn)場(chǎng)耐久性試驗(yàn)流程

    在試驗(yàn)場(chǎng)與用戶等效轉(zhuǎn)化的過(guò)程中,選取載荷損傷特征D作為用戶載荷與試驗(yàn)場(chǎng)載荷一致性的評(píng)判指標(biāo)。由于這一損傷特征(確切來(lái)說(shuō)是偽損傷)的計(jì)算簡(jiǎn)單,適用于載荷的對(duì)比和等效關(guān)聯(lián),因此在整車(chē)及零部件耐久性工程中獲得了廣泛的應(yīng)用[12],計(jì)算表達(dá)式如下:

    SrN=C

    (1)

    (2)

    式中,S為廣義應(yīng)力幅值;N為被試件在幅值S下的理論疲勞壽命;C為材料常數(shù);r為疲勞強(qiáng)度指數(shù),與材料特性有關(guān),對(duì)于汽車(chē)零部件而言一般取值r=5;dk為給定k級(jí)應(yīng)力幅值Sk時(shí)的損傷值;nk、Nk分別為給定k級(jí)應(yīng)力幅值Sk時(shí),隨機(jī)載荷在該應(yīng)力幅值作用下的頻次和理論疲勞壽命。

    由于在驅(qū)動(dòng)橋傳動(dòng)系統(tǒng)中,除了承受路面激勵(lì)產(chǎn)生的沖擊外,還主要受到變速器輸入扭矩的作用而產(chǎn)生疲勞破壞,因此,以驅(qū)動(dòng)橋左右兩半軸扭矩為測(cè)量對(duì)象,將兩半軸扭矩?cái)?shù)據(jù)按擋位劃分為12個(gè)通道,從而構(gòu)建出試驗(yàn)場(chǎng)載荷系數(shù)矩陣與用戶載荷目標(biāo)矩陣的等損傷關(guān)聯(lián)模型,其表達(dá)式為

    Dβ=T

    (3)

    β=(β1,β2,…,βj)TT=(T1,T2,…,Ti)T

    式中,Dij為兩半軸第i通道在試驗(yàn)場(chǎng)路況j單次循環(huán)下的系數(shù)損傷值;βj為試驗(yàn)場(chǎng)路況j的循環(huán)次數(shù);Ti為用戶使用條件下第i通道的目標(biāo)損傷值。

    2 道路載荷采集與處理

    2.1 傳感器的安裝

    由于本文研究的內(nèi)容是驅(qū)動(dòng)橋傳動(dòng)系統(tǒng),對(duì)驅(qū)動(dòng)橋半軸兩端扭矩、轉(zhuǎn)速及擋位感興趣,因此在商用車(chē)后橋兩車(chē)輪端安裝六分力傳感器,如圖2所示。擋位信號(hào)可以利用CAN信號(hào)讀取。安裝GPS傳感器以便了解車(chē)速及采集線路。除此之外,為區(qū)分出后續(xù)采集路況和工況,利用邏輯開(kāi)關(guān)進(jìn)行判別。

    圖2 六分力傳感器布置示意圖

    2.2 用戶道路載荷采集

    調(diào)查發(fā)現(xiàn),該商用車(chē)用戶群體通過(guò)高速公路、城市道路、一般道路和惡劣道路的里程比例分別為40%、30%、20%、10%。其中,鑒于對(duì)國(guó)道、省道和縣道難以準(zhǔn)確區(qū)分的情況,將三種路況合并為同一類(lèi),記為一般道路。惡劣道路是指用戶通過(guò)坑洼地段以及山路、碎石路等路況差路面的統(tǒng)稱。由于市場(chǎng)調(diào)查和用戶采集所得到的數(shù)據(jù)有限,無(wú)法做到無(wú)限逼近用戶真實(shí)的使用情況[13-14],因此,本文利用多組駕駛員在典型用戶路面下實(shí)測(cè)的載荷進(jìn)行分割、外推和疊加來(lái)等效模擬用戶200 000 km總體樣本。表1所示為采集用戶特征路面信息。

    表1 路面特征參數(shù)

    2.3 試驗(yàn)場(chǎng)道路載荷采集

    按照現(xiàn)有的試車(chē)操作規(guī)范,依次在我國(guó)某大型試驗(yàn)場(chǎng)中的各個(gè)特征路況進(jìn)行兩半軸扭矩單次循環(huán)采集,如圖3所示。該試驗(yàn)場(chǎng)主要工況包括強(qiáng)化路工況、動(dòng)力工況及高速工況等,為了保證后續(xù)試驗(yàn)的連貫性,將由各種路面組成的路況作為試驗(yàn)規(guī)范工況之一,共劃分為10個(gè)子工況,下文簡(jiǎn)稱F1~F10(依次對(duì)應(yīng)j=1,2,…,10)。采集時(shí)考慮到不同用戶駕駛習(xí)慣的差異性,分別由3位駕駛員各采集3次。表2所示為試驗(yàn)場(chǎng)采集工況信息[12],其中l(wèi)j為路況j的單次里程。

    圖3 試驗(yàn)場(chǎng)各路況分布圖

    表2 各路況特征信息

    2.4 旋轉(zhuǎn)雨流計(jì)數(shù)

    為構(gòu)建等損傷關(guān)聯(lián)模型,首先需要進(jìn)行載荷的統(tǒng)計(jì)計(jì)數(shù),以求得載荷與頻次之間的對(duì)應(yīng)關(guān)系。由于驅(qū)動(dòng)橋旋轉(zhuǎn)件主要有半軸及齒輪等零部件,因此考慮利用旋轉(zhuǎn)雨流計(jì)數(shù)方法得到驅(qū)動(dòng)橋的軸載荷和齒載荷。軸載荷代表半軸扭矩雨流計(jì)數(shù)結(jié)果,是考核半軸臺(tái)架試驗(yàn)的重要數(shù)據(jù)來(lái)源;根據(jù)半軸旋轉(zhuǎn)一周、單齒嚙合一次的規(guī)律,勾勒出單齒所受脈沖扭矩時(shí)間歷程,再進(jìn)行常規(guī)的雨流計(jì)數(shù),結(jié)果記為齒載荷,如圖4所示。由于考核部件為驅(qū)動(dòng)橋總成,故最終將軸載荷和齒載荷進(jìn)行疊加得到載荷與頻次的結(jié)果記為T(mén)otal,如圖5所示。

    圖4 齒載荷構(gòu)建示意圖

    3 用戶載荷分布模型預(yù)測(cè)

    3.1 非參數(shù)二維核密度估計(jì)

    用戶載荷分布模型的構(gòu)建是車(chē)輛耐久性工程中比較稀缺的一個(gè)頂層輸入,一旦模型的構(gòu)建出現(xiàn)偏差,將會(huì)影響全部的下游工作。為此,需要對(duì)實(shí)際采集一定里程的載荷估計(jì)其概率密度分布函數(shù),從而合理地預(yù)測(cè)全壽命周期內(nèi)的用戶載荷分布情況。

    由于道路載荷的高度隨機(jī)性及環(huán)境復(fù)雜性,通常難以用特定的概率密度函數(shù)對(duì)載荷分布進(jìn)行擬合,而非參數(shù)估計(jì)方法并不需要假設(shè)載荷服從某一特定分布。本文基于核密度估計(jì)原理[15],結(jié)合From-To矩陣形成遲滯回環(huán)的計(jì)數(shù)特點(diǎn),預(yù)測(cè)該遲滯回環(huán)在全壽命周期內(nèi)可能出現(xiàn)的頻次。

    將前文旋轉(zhuǎn)雨流計(jì)數(shù)結(jié)果Total載荷的From-To矩陣進(jìn)行核密度估計(jì),概率密度函數(shù)記為f(x,y),其中x、y分別為遲滯回環(huán)的起始變量(即圖5中的坐標(biāo)軸From)和終止變量(即圖5中的坐標(biāo)軸To);選取高斯核作為核密度估計(jì)中的核函數(shù)K:

    (a)軸齒載

    (4)

    核函數(shù)K的初始核帶寬參數(shù)h的表達(dá)式為

    h=2.4σa-1/5

    (5)

    式中,σ為雨流矩陣載荷的標(biāo)準(zhǔn)差;a為雨流矩陣樣本容量。

    根據(jù)圖5中載荷分布的稀疏程度,引入計(jì)算自適應(yīng)因子從而獲得求解自適應(yīng)帶寬的方法,對(duì)雨流矩陣進(jìn)行自適應(yīng)核密度估計(jì),并計(jì)算自適應(yīng)因子:

    (6)

    式中,(xl,yl)為第l個(gè)數(shù)據(jù)點(diǎn);p(xl,yl)表示點(diǎn)(xl,yl)出現(xiàn)的概率,其值為點(diǎn)(xl,yl)處頻次與累積頻次之比。

    確定自適應(yīng)帶寬及核密度函數(shù)K(x,y)后,便可獲得雨流矩陣核密度估計(jì)的概率密度函數(shù)數(shù)學(xué)模型,其表達(dá)式為

    (7)

    3.2 旋轉(zhuǎn)雨流矩陣的隨機(jī)模擬外推

    用戶載荷的核密度估計(jì)為連續(xù)函數(shù),但為了預(yù)測(cè)各擋位下的隨機(jī)載荷,需將載荷的概率密度按載荷分級(jí)離散到載荷外推極值范圍內(nèi)的區(qū)域,并形成一個(gè)離散概率分布。運(yùn)用蒙特卡羅方法隨機(jī)放置載荷遲滯循環(huán),從而獲取全壽命周期內(nèi)的用戶載荷分布模型[16]。

    用戶道路載荷采集信息如表2所示,根據(jù)各工況所占比例可求出高速、城市、一般及惡劣道路四種工況達(dá)到目標(biāo)里程的外推系數(shù)分別為K1=1427、K2=1302、K3=610、K4=580。

    基于上述理論依據(jù),以惡劣道路工況下左半軸1擋的載荷分布為例,計(jì)算可得σ=3.0,a= 175,初始帶寬參數(shù)h=1.1,在LMS Tecware中設(shè)置外推系數(shù)K4,圖6所示為最終預(yù)測(cè)用戶載荷分布模型,圖7所示為三種載荷極值A(chǔ)作用下的頻次。

    圖6 用戶Total載荷From-To矩陣核密度估計(jì)

    圖7 全壽命周期內(nèi)用戶載荷分布曲線

    4 試驗(yàn)場(chǎng)試驗(yàn)規(guī)范制定

    4.1 約束條件設(shè)置

    求解試驗(yàn)場(chǎng)路況循環(huán)次數(shù)時(shí),為實(shí)現(xiàn)采集過(guò)程中各路況的連貫性,需要設(shè)置約束條件,從而保證解集的合理性。

    4.1.1極值范圍

    為保證試驗(yàn)在試驗(yàn)場(chǎng)各路況下均能開(kāi)展,設(shè)定各工況循環(huán)次數(shù)下限大于零,循環(huán)次數(shù)上限為試驗(yàn)場(chǎng)總里程除以每個(gè)路況的單次循環(huán)里程,可表示為

    (8)

    式中,L為試驗(yàn)場(chǎng)總里程,取L=50 000 km;lj為路況j的單次里程,單位為km,其數(shù)值大小見(jiàn)表2。

    4.1.2等式約束

    強(qiáng)化路工況包含強(qiáng)化環(huán)道路F1~F5和評(píng)價(jià)路F6,其中強(qiáng)化環(huán)道路各路況相互連接,評(píng)價(jià)路相對(duì)獨(dú)立,因此,根據(jù)試驗(yàn)場(chǎng)道路車(chē)輛測(cè)試路徑及各環(huán)道的權(quán)重可確定5個(gè)環(huán)道循環(huán)次數(shù)之間的相互比例(即各個(gè)環(huán)道試驗(yàn)循環(huán)次數(shù)之間的比值)為1∶2∶1∶1∶1,矩陣方程如下:

    (9)

    4.1.3不等式約束

    為防止試驗(yàn)場(chǎng)道路耐久性試驗(yàn)的總試驗(yàn)里程過(guò)長(zhǎng),根據(jù)經(jīng)驗(yàn)將試驗(yàn)總里程控制在50 000 km內(nèi),強(qiáng)化路工況里程控制在6000~8000 km范圍內(nèi),從而達(dá)到加速試驗(yàn)的目的,即

    (10)

    4.2 基于最小二乘法的求解

    根據(jù)圖7獲得的Total載荷與頻次關(guān)系,基于用戶關(guān)聯(lián)技術(shù),得到的方程組如式(3)所示。由于系數(shù)損傷矩陣D中i>j,D為列滿秩,因此方程組沒(méi)有精確解,此時(shí)稱為超定方程組[17]。為獲得適合的循環(huán)次數(shù)β,通常構(gòu)建多元線性回歸模型進(jìn)行求解,將式(3)所示的矩陣方程改寫(xiě)成多項(xiàng)式形式:

    Ti=Dijβ1+Dijβ2+…+Dijβj

    (11)

    其中,第i(i=1,2,…,12)通道下的用戶載荷目標(biāo)損傷值Ti見(jiàn)圖8;第i通道下,路況j(j=1,2,…,10)的試驗(yàn)場(chǎng)載荷系數(shù)損傷值Dij見(jiàn)圖9。

    圖8 用戶載荷目標(biāo)損傷

    圖9 試驗(yàn)場(chǎng)載荷系數(shù)損傷

    (12)

    由最小二乘原理可知:

    (13)

    對(duì)Q求其一階偏導(dǎo)并令各項(xiàng)為零,可得

    (14)

    進(jìn)一步轉(zhuǎn)化式(3)所示的矩陣方程可得

    DTDβ=DTT

    (15)

    則方程組解為

    β=(DTD)-1DTT

    (16)

    通過(guò)上述求解,在MATLAB中編程,設(shè)置極值、等式、不等式約束條件,獲取各路況最小二乘解(即循環(huán)次數(shù)),如圖10所示。由于約束條件的存在,使得該解集中強(qiáng)化路部分工況出現(xiàn)零次循環(huán)的情況,考慮到試驗(yàn)的完整性及解的集中性,該解集并不理想,需要進(jìn)一步進(jìn)行優(yōu)化。

    圖10 各工況循環(huán)次數(shù)

    4.3 基于NSGA-Ⅱ算法的求解

    4.3.1多目標(biāo)優(yōu)化模型構(gòu)建

    基于最小二乘法原理求其局部最優(yōu)解時(shí),由于解的分散性高,易出現(xiàn)零解的情況,不利于試驗(yàn)規(guī)范的制定,因此,為便于驅(qū)動(dòng)橋試驗(yàn)場(chǎng)規(guī)范的求解,本文提出將等損傷關(guān)聯(lián)方程轉(zhuǎn)化為多目標(biāo)優(yōu)化問(wèn)題的求解,其數(shù)學(xué)表達(dá)式如下:

    fi(x:D,T)=

    |(Di1,Di2,…,Di10)·(x1,x2,…,x10)T-(T1,T2,…,Ti)T|

    (17)

    s.t.blj≤xj≤buj

    Aeq·x=beq

    Ajxj≤bj

    式中:i=1,2,…,12;xj為利用多目標(biāo)優(yōu)化算法下求解出的路況j的循環(huán)次數(shù),j=1,2,…,10;fi(x)表示第i通道下的目標(biāo)函數(shù);bu、bl分別為上下限矩陣,分別表示極值約束中各路況的最高循環(huán)次數(shù)和最低循環(huán)次數(shù);Aeq、beq分別為等式約束中試驗(yàn)場(chǎng)各工況的權(quán)重系數(shù)矩陣和零矩陣;A為試驗(yàn)場(chǎng)各路況單次循環(huán)里程矩陣;b為試驗(yàn)場(chǎng)總里程和強(qiáng)化路總里程限定矩陣。

    帶精英策略的非支配排序遺傳算法(NSGA-Ⅱ)能夠通過(guò)一次運(yùn)算找到使各目標(biāo)函數(shù)達(dá)到比較理想的Pareto最優(yōu)解xj(也稱為非劣解[18]),具有運(yùn)行速度快、解集收斂性好的優(yōu)點(diǎn)。NSGA-Ⅱ算法的原理如圖11所示,其中G為種群迭代次數(shù)。

    圖11 NSGA-Ⅱ原理圖

    將目標(biāo)函數(shù)及各約束條件建立完成后,在MATLAB中編程,設(shè)置變量為10、最優(yōu)前端個(gè)體系數(shù)為0.3,種群規(guī)模的大小決定著收斂速度的快慢,一般在50~200之間,本文根據(jù)多次運(yùn)算結(jié)果最終選取種群規(guī)模為200。迭代上限設(shè)為500,基于NSGA-Ⅱ算法,迭代134次后,得到收斂的77組Pareto最優(yōu)解集,每組非劣解集中有10個(gè)解,表示試驗(yàn)場(chǎng)各工況下的循環(huán)次數(shù),如表3所示。

    表3 非劣解集

    由解的分布可知,環(huán)道F1~F5對(duì)應(yīng)循環(huán)次數(shù)x1~x5的權(quán)重關(guān)系滿足設(shè)置的約束條件。此外,解集被限制在[bl,bu]范圍內(nèi),將各工況循環(huán)次數(shù)與各路況單程里程相乘并累加得出的總試驗(yàn)里程均小于50 000 km,強(qiáng)化路工況里程在6000~8000 km以內(nèi)。綜上可得,解的情況滿足設(shè)計(jì)要求。

    4.3.2最優(yōu)解的選取

    分別將77組非劣解xmj(m=1,2,…,77)替換βj,可得到12個(gè)通道的關(guān)聯(lián)損傷值Tmi,將12個(gè)通道的關(guān)聯(lián)損傷值Tmi與目標(biāo)損傷值Ti之間的偽損傷比記為相對(duì)損傷:

    (18)

    式中,Pmi為第i通道下,第m組非劣解的相對(duì)損傷;Tmi為第i個(gè)通道下,第m組非劣解的關(guān)聯(lián)損傷。

    Pmi值越接近1,表明等效效果越好,最終得到的77組相對(duì)損傷如表4所示。

    表4 相對(duì)損傷

    由表4可知,相對(duì)損傷大部分都成正相關(guān),即各個(gè)通道的相對(duì)損傷值之間相互影響,一個(gè)值的增大或減小都會(huì)引起全局的變化,求解時(shí)不能同時(shí)保證全部的解都靠近1,因此選取相對(duì)損傷平均值接近1,標(biāo)準(zhǔn)差最小,即將波動(dòng)最小的一組非劣解作為最優(yōu)解[19]。通過(guò)多次運(yùn)算,選取平均值μ=1、標(biāo)準(zhǔn)差σ=0.31的第6組作為目標(biāo)解集。

    4.4 試驗(yàn)工況組合

    考慮到試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范行駛路徑規(guī)劃的合理性,需要在基于NSGA-Ⅱ算得的試驗(yàn)場(chǎng)各工況下循環(huán)次數(shù)的基礎(chǔ)上,根據(jù)工程經(jīng)驗(yàn)適當(dāng)進(jìn)行優(yōu)化。由第6組解的特性可知,取各工況循環(huán)次數(shù)為100的整數(shù)倍,從而方便試驗(yàn)場(chǎng)耐久性試驗(yàn)的進(jìn)行。所制定的驅(qū)動(dòng)橋試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范如表5所示。

    表5 試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范

    4.4.1相對(duì)損傷驗(yàn)證

    一般認(rèn)為關(guān)聯(lián)損傷通道滿足相對(duì)損傷P在0.5~2 范圍內(nèi)即視為可靠,因此將基于NSGA-Ⅱ優(yōu)化后的試驗(yàn)規(guī)范代入式(18)中,得到各通道下的相對(duì)損傷值,并計(jì)算其平均值和標(biāo)準(zhǔn)差,結(jié)果如圖12所示,可以看出,相對(duì)損傷值分布在1附近,其范圍穩(wěn)定在0.6~1.5之間,均值為0.96,標(biāo)準(zhǔn)差為0.28。但由于用戶與試驗(yàn)場(chǎng)工況的差異性以及目標(biāo)函數(shù)的復(fù)雜性,相對(duì)損傷值完全等于1是較難實(shí)現(xiàn)的。

    圖12 各通道相對(duì)損傷分布

    4.4.2載荷分布驗(yàn)證

    根據(jù)新的試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范,可以得到試驗(yàn)場(chǎng)等效于用戶路面200 000 km的載荷幅值-累積頻次曲線,通過(guò)與用戶載荷幅值-累積頻次曲線進(jìn)行對(duì)比,可以分析出在該試驗(yàn)規(guī)范下載荷的分布情況是否合理。以兩半軸1擋載荷分布為例,對(duì)比結(jié)果如圖13所示,可以看出,用戶與試驗(yàn)場(chǎng)的載荷分布趨勢(shì)較為一致。結(jié)合圖12可知,相對(duì)損傷值大于1時(shí)說(shuō)明試驗(yàn)場(chǎng)路面載荷損傷大于用戶路面載荷損傷,反之,相對(duì)損傷值小于1時(shí)則說(shuō)明試驗(yàn)場(chǎng)路面載荷損傷小于用戶路面載荷損傷。由圖13還可知,用戶載荷累積頻次多于試驗(yàn)場(chǎng)載荷累積頻次,依據(jù)式(2)所示的損傷計(jì)算原理,累積頻次越多,損傷值越大,即用戶路面載荷損傷略大于試驗(yàn)場(chǎng)路面載荷損傷,相對(duì)損傷小于1。結(jié)合圖13中用戶1擋載荷累積頻次略多于試驗(yàn)場(chǎng)1擋載荷累積頻次,則對(duì)應(yīng)圖12中1擋關(guān)聯(lián)通道1和7的相對(duì)損傷小于1,說(shuō)明關(guān)聯(lián)計(jì)算結(jié)果較理想。

    (a)左半軸1擋

    綜上所述,新的驅(qū)動(dòng)橋試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范在相對(duì)損傷及載荷分布兩方面均能較好地復(fù)現(xiàn)用戶實(shí)際使用情況的效果,滿足驅(qū)動(dòng)橋試驗(yàn)場(chǎng)耐久性試驗(yàn)的要求。

    5 強(qiáng)化系數(shù)計(jì)算

    試驗(yàn)場(chǎng)路面強(qiáng)化系數(shù)Ks是指在相同的失效模式下,驅(qū)動(dòng)橋在用戶路面上行駛的里程與試驗(yàn)場(chǎng)試驗(yàn)里程之間的比值,即

    (19)

    式中,Le為用戶路面行駛里程,km;Lp為試驗(yàn)場(chǎng)試驗(yàn)里程,km。

    計(jì)算可得到強(qiáng)化系數(shù)Ks=4.34,即當(dāng)用戶路面行駛200 000 km時(shí),在試驗(yàn)場(chǎng)僅完成46 133.3 km就能達(dá)到考核驅(qū)動(dòng)橋耐久性的要求。

    6 當(dāng)量關(guān)系計(jì)算

    試驗(yàn)場(chǎng)路況種類(lèi)多達(dá)幾十種,每種路況下零部件的載荷損傷各有不同。為比較各工況對(duì)驅(qū)動(dòng)橋零部件的損傷大小,需計(jì)算單位里程下各工況的載荷損傷與總損傷之間的當(dāng)量關(guān)系,從而為臺(tái)架試驗(yàn)加載及有限元分析的載荷輸入提供參考。

    強(qiáng)化路工況包含的路面種類(lèi)較為惡劣,如比利時(shí)鋪裝路、卵石路、坑洼路等;動(dòng)力工況對(duì)動(dòng)力系統(tǒng)要求較高,包含加速、制動(dòng)、爬坡等工況;高速工況主要考核驅(qū)動(dòng)橋部件保持高速運(yùn)轉(zhuǎn)的能力。根據(jù)表5規(guī)定的試驗(yàn)工況,比較單位里程下各工況載荷損傷與總損傷之間的當(dāng)量關(guān)系。如圖14所示,動(dòng)力工況對(duì)驅(qū)動(dòng)橋零部件的損傷最大,高速工況次之,強(qiáng)化路工況最小。研究結(jié)果表明:路面惡劣程度(即路面激勵(lì)對(duì)驅(qū)動(dòng)橋旋轉(zhuǎn)件的疲勞損傷影響)不大,使驅(qū)動(dòng)橋動(dòng)力傳動(dòng)系統(tǒng)產(chǎn)生疲勞破壞最根本的原因在于變速器輸入載荷幅值的作用頻次對(duì)疲勞耐久性能的影響,而與加速度、應(yīng)變等信號(hào)無(wú)直接關(guān)系,從側(cè)面反映出將不同擋位下的扭矩?fù)p傷作為用戶-試驗(yàn)場(chǎng)關(guān)聯(lián)特征的正確性。

    圖14 單位里程下各工況損傷占總損傷百分比

    7 結(jié)論

    (1)基于旋轉(zhuǎn)雨流計(jì)數(shù)和核密度估計(jì)方法,預(yù)測(cè)了驅(qū)動(dòng)橋軸載荷和齒載荷分布模型,為更全面地反映驅(qū)動(dòng)橋的各零部件受載特性以及建立“用戶-試驗(yàn)場(chǎng)”等損傷關(guān)聯(lián)方程奠定了基礎(chǔ)。

    (2)遺傳算法得到的解避免了最小二乘法出現(xiàn)零解的情況,其分散性較低,各通道相對(duì)損傷比值穩(wěn)定在0.6~1.5范圍內(nèi),與用戶損傷吻合。且由各種特征路面組成的試驗(yàn)規(guī)范試驗(yàn)里程為46 133.3 km時(shí)相當(dāng)于用戶行駛200 000 km,極大地縮短了試驗(yàn)時(shí)間并降低了成本。

    (3)驅(qū)動(dòng)橋傳動(dòng)系統(tǒng)產(chǎn)生損傷的主要因素來(lái)源于變速器輸入扭矩作用下的彎曲疲勞破壞,路面種類(lèi)對(duì)驅(qū)動(dòng)橋傳動(dòng)系統(tǒng)疲勞耐久性能的影響較小。所制定的試驗(yàn)場(chǎng)耐久性試驗(yàn)規(guī)范符合我國(guó)用戶實(shí)際行駛情況,本文的思路和方法為驅(qū)動(dòng)橋傳動(dòng)系統(tǒng)乃至整車(chē)試驗(yàn)場(chǎng)可靠性和耐久性試驗(yàn)規(guī)范的制定提供了參考和依據(jù)。

    猜你喜歡
    試驗(yàn)場(chǎng)路況里程
    高速公路路況信息系統(tǒng)
    汽車(chē)試驗(yàn)場(chǎng)噪音試驗(yàn)路吸聲系數(shù)研究
    石油瀝青(2022年3期)2022-08-26 09:13:12
    關(guān)于國(guó)內(nèi)主要汽車(chē)試驗(yàn)場(chǎng)運(yùn)營(yíng)情況的調(diào)研報(bào)告
    近海海洋儀器試驗(yàn)場(chǎng)微波通信系統(tǒng)構(gòu)建與測(cè)試
    世界上最大的海上無(wú)人系統(tǒng)試驗(yàn)場(chǎng)開(kāi)放
    從路況報(bào)道看廣播“類(lèi)型化”新聞的要素構(gòu)成
    騰勢(shì)400 用在上海市區(qū)的來(lái)回穿梭克服里程焦慮
    車(chē)迷(2017年12期)2018-01-18 02:16:12
    幸福合力 開(kāi)啟幸福里程
    幸福合力 開(kāi)啟幸福里程
    算里程
    香蕉av资源在线| 亚洲18禁久久av| 18禁观看日本| 国产精品1区2区在线观看.| 少妇熟女aⅴ在线视频| 国产精品一区二区精品视频观看| 国产精品99久久99久久久不卡| 夜夜看夜夜爽夜夜摸| 亚洲成人久久性| 日本五十路高清| 午夜福利欧美成人| 国产主播在线观看一区二区| 黄片小视频在线播放| 国产亚洲欧美在线一区二区| 嫩草影院精品99| 国产成人精品久久二区二区91| 国产成人aa在线观看| 怎么达到女性高潮| 亚洲中文字幕日韩| 亚洲,欧美精品.| 69av精品久久久久久| 国产精品亚洲av一区麻豆| 丁香六月欧美| 午夜福利在线在线| 香蕉丝袜av| 在线看三级毛片| 国产1区2区3区精品| 日韩中文字幕欧美一区二区| 亚洲精品美女久久久久99蜜臀| 欧美成人一区二区免费高清观看 | 色尼玛亚洲综合影院| 少妇熟女aⅴ在线视频| 在线视频色国产色| 国产成人精品无人区| 亚洲第一欧美日韩一区二区三区| 精品福利观看| 亚洲激情在线av| 亚洲精华国产精华精| 美女扒开内裤让男人捅视频| 午夜两性在线视频| 在线免费观看的www视频| 午夜福利在线观看吧| 亚洲精品国产一区二区精华液| 亚洲欧美日韩东京热| 国产亚洲精品第一综合不卡| 51午夜福利影视在线观看| 桃色一区二区三区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 我的老师免费观看完整版| 精品乱码久久久久久99久播| 麻豆成人午夜福利视频| 在线国产一区二区在线| 亚洲国产精品sss在线观看| av在线天堂中文字幕| 亚洲欧美激情综合另类| 黑人操中国人逼视频| 日韩欧美国产一区二区入口| 极品教师在线免费播放| 欧美黄色淫秽网站| 欧美大码av| 一本大道久久a久久精品| 国产野战对白在线观看| 国产高清有码在线观看视频 | 成人av一区二区三区在线看| 国产精品自产拍在线观看55亚洲| 免费看日本二区| 国产免费男女视频| 两个人免费观看高清视频| 啪啪无遮挡十八禁网站| 天堂影院成人在线观看| 麻豆久久精品国产亚洲av| 长腿黑丝高跟| 免费在线观看成人毛片| www日本在线高清视频| 亚洲自拍偷在线| 亚洲性夜色夜夜综合| 91国产中文字幕| 一本精品99久久精品77| 日本免费一区二区三区高清不卡| 国内毛片毛片毛片毛片毛片| 国产v大片淫在线免费观看| 超碰成人久久| 一级a爱片免费观看的视频| 免费观看人在逋| 97超级碰碰碰精品色视频在线观看| 一本综合久久免费| av在线播放免费不卡| 黄频高清免费视频| 国产人伦9x9x在线观看| 少妇熟女aⅴ在线视频| 久久精品人妻少妇| 亚洲精品美女久久久久99蜜臀| 午夜亚洲福利在线播放| 精品久久蜜臀av无| 亚洲成人免费电影在线观看| 一个人免费在线观看的高清视频| 99热这里只有精品一区 | 人妻夜夜爽99麻豆av| 床上黄色一级片| bbb黄色大片| 久久精品国产综合久久久| 亚洲乱码一区二区免费版| 精品国产乱码久久久久久男人| 欧美黄色淫秽网站| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女| 他把我摸到了高潮在线观看| 国产亚洲欧美在线一区二区| 一本一本综合久久| svipshipincom国产片| 国产片内射在线| 最新在线观看一区二区三区| bbb黄色大片| 久久草成人影院| 久久久久久久久免费视频了| 欧美精品啪啪一区二区三区| 午夜福利高清视频| 久久婷婷人人爽人人干人人爱| 国产精品自产拍在线观看55亚洲| 久久香蕉激情| √禁漫天堂资源中文www| 国产亚洲精品综合一区在线观看 | 亚洲人成77777在线视频| tocl精华| 久久午夜综合久久蜜桃| 成人国产综合亚洲| 99国产极品粉嫩在线观看| 中文字幕人妻丝袜一区二区| 99久久无色码亚洲精品果冻| 久久精品国产99精品国产亚洲性色| 黄色毛片三级朝国网站| 夜夜看夜夜爽夜夜摸| 日本 av在线| 这个男人来自地球电影免费观看| 一区二区三区国产精品乱码| 久久精品aⅴ一区二区三区四区| 桃红色精品国产亚洲av| 久久久久免费精品人妻一区二区| 久久天堂一区二区三区四区| 亚洲人成伊人成综合网2020| 成人国产综合亚洲| av有码第一页| 99热这里只有精品一区 | 亚洲精品在线观看二区| 亚洲国产欧美一区二区综合| 国产1区2区3区精品| 国产成人系列免费观看| 1024手机看黄色片| 亚洲国产精品合色在线| 欧美黑人欧美精品刺激| 精品一区二区三区四区五区乱码| 亚洲国产精品sss在线观看| 国产高清激情床上av| 亚洲av成人av| 99riav亚洲国产免费| 国产亚洲欧美98| 免费一级毛片在线播放高清视频| 日本黄色视频三级网站网址| 国产真人三级小视频在线观看| 亚洲无线在线观看| 麻豆成人午夜福利视频| 两个人看的免费小视频| 亚洲精品一卡2卡三卡4卡5卡| 巨乳人妻的诱惑在线观看| videosex国产| 久久午夜亚洲精品久久| 九色国产91popny在线| 变态另类成人亚洲欧美熟女| 18禁黄网站禁片免费观看直播| 欧美日韩乱码在线| 欧美最黄视频在线播放免费| 久久精品aⅴ一区二区三区四区| 麻豆国产av国片精品| 91九色精品人成在线观看| 男女做爰动态图高潮gif福利片| 欧美色视频一区免费| 国产午夜精品论理片| 日韩欧美一区二区三区在线观看| av片东京热男人的天堂| 日本黄大片高清| 听说在线观看完整版免费高清| 国产成人欧美在线观看| 777久久人妻少妇嫩草av网站| 午夜a级毛片| www日本黄色视频网| 亚洲精品美女久久久久99蜜臀| 久久精品成人免费网站| 久久午夜亚洲精品久久| 日本三级黄在线观看| 老司机在亚洲福利影院| 亚洲电影在线观看av| av在线播放免费不卡| 丝袜美腿诱惑在线| x7x7x7水蜜桃| 精品高清国产在线一区| 老司机靠b影院| 国产精品一区二区三区四区久久| 又大又爽又粗| 久久伊人香网站| 黄色成人免费大全| 国产精品av视频在线免费观看| 久久精品国产亚洲av高清一级| 久久久久国产一级毛片高清牌| 欧美绝顶高潮抽搐喷水| 国产精品久久久av美女十八| 国内毛片毛片毛片毛片毛片| 一级毛片女人18水好多| 国产高清激情床上av| 国产成人aa在线观看| 免费无遮挡裸体视频| 变态另类丝袜制服| 中文字幕最新亚洲高清| 精华霜和精华液先用哪个| 日本一二三区视频观看| 精品福利观看| 亚洲美女视频黄频| 99在线视频只有这里精品首页| 久久精品亚洲精品国产色婷小说| 五月伊人婷婷丁香| 91国产中文字幕| 日韩国内少妇激情av| 亚洲激情在线av| 最新在线观看一区二区三区| 精品国产超薄肉色丝袜足j| av在线天堂中文字幕| 欧美久久黑人一区二区| av欧美777| 无人区码免费观看不卡| 中文亚洲av片在线观看爽| 黄色视频,在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 2021天堂中文幕一二区在线观| 男人的好看免费观看在线视频 | 国产欧美日韩一区二区精品| 国内精品一区二区在线观看| 成在线人永久免费视频| 窝窝影院91人妻| 国产高清激情床上av| x7x7x7水蜜桃| 国产爱豆传媒在线观看 | 少妇粗大呻吟视频| 亚洲精品久久国产高清桃花| 国产熟女xx| 精品乱码久久久久久99久播| 日本在线视频免费播放| 久久人人精品亚洲av| 成人手机av| 一本大道久久a久久精品| АⅤ资源中文在线天堂| 91大片在线观看| 成人三级黄色视频| 国产日本99.免费观看| 动漫黄色视频在线观看| 久久久久国内视频| 激情在线观看视频在线高清| 琪琪午夜伦伦电影理论片6080| 亚洲欧洲精品一区二区精品久久久| 亚洲精品在线观看二区| bbb黄色大片| 中文亚洲av片在线观看爽| 久久精品国产亚洲av高清一级| 亚洲欧美日韩高清在线视频| 国产精品综合久久久久久久免费| 神马国产精品三级电影在线观看 | 老熟妇乱子伦视频在线观看| 国产精品,欧美在线| 女警被强在线播放| 99精品在免费线老司机午夜| 婷婷精品国产亚洲av在线| 成在线人永久免费视频| 精品久久蜜臀av无| 在线观看66精品国产| 99在线人妻在线中文字幕| 2021天堂中文幕一二区在线观| 91国产中文字幕| 99久久99久久久精品蜜桃| 欧美在线黄色| 一边摸一边抽搐一进一小说| 黄色毛片三级朝国网站| 又粗又爽又猛毛片免费看| 日本a在线网址| 久久国产精品人妻蜜桃| 一区二区三区激情视频| 国产一区二区在线av高清观看| 日本撒尿小便嘘嘘汇集6| 色老头精品视频在线观看| 午夜免费观看网址| 亚洲成人久久爱视频| 亚洲专区字幕在线| 丰满的人妻完整版| 欧美人与性动交α欧美精品济南到| 国产熟女xx| 99久久精品热视频| 在线观看日韩欧美| 搡老熟女国产l中国老女人| 亚洲美女视频黄频| 国产真实乱freesex| 久久久久免费精品人妻一区二区| 美女大奶头视频| 国产精华一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 午夜福利18| 国内毛片毛片毛片毛片毛片| 特大巨黑吊av在线直播| 99国产精品99久久久久| 久久人人精品亚洲av| 日本 欧美在线| 淫妇啪啪啪对白视频| 久久国产乱子伦精品免费另类| 久久久久国内视频| 50天的宝宝边吃奶边哭怎么回事| 国产野战对白在线观看| 午夜亚洲福利在线播放| 成年免费大片在线观看| 亚洲欧美精品综合一区二区三区| 麻豆国产97在线/欧美 | 日本熟妇午夜| 日本 av在线| 日韩有码中文字幕| 亚洲最大成人中文| 88av欧美| 首页视频小说图片口味搜索| 欧美成人午夜精品| 国产精品日韩av在线免费观看| 一本久久中文字幕| 男男h啪啪无遮挡| 婷婷六月久久综合丁香| 日韩有码中文字幕| 亚洲午夜精品一区,二区,三区| 国产高清激情床上av| 真人一进一出gif抽搐免费| 韩国av一区二区三区四区| 久久这里只有精品中国| 日本 欧美在线| 亚洲 欧美一区二区三区| 国产不卡一卡二| 精品无人区乱码1区二区| 全区人妻精品视频| 人人妻,人人澡人人爽秒播| 白带黄色成豆腐渣| 午夜免费成人在线视频| 亚洲性夜色夜夜综合| 亚洲国产中文字幕在线视频| 日韩欧美 国产精品| 一级片免费观看大全| 韩国av一区二区三区四区| 亚洲中文av在线| 伦理电影免费视频| 每晚都被弄得嗷嗷叫到高潮| 久久精品91无色码中文字幕| 国产片内射在线| 麻豆久久精品国产亚洲av| 99久久无色码亚洲精品果冻| 国产午夜精品论理片| 丰满人妻一区二区三区视频av | 黄色丝袜av网址大全| 女同久久另类99精品国产91| 免费看十八禁软件| 中文亚洲av片在线观看爽| 久久天躁狠狠躁夜夜2o2o| 久久久久国产精品人妻aⅴ院| 香蕉丝袜av| 国产午夜精品论理片| 成在线人永久免费视频| 一本大道久久a久久精品| 久久中文字幕一级| 一本大道久久a久久精品| 波多野结衣高清作品| 成熟少妇高潮喷水视频| 亚洲精品在线美女| 久久婷婷成人综合色麻豆| av中文乱码字幕在线| 国产精品久久久久久人妻精品电影| 国产又黄又爽又无遮挡在线| 午夜福利高清视频| 桃红色精品国产亚洲av| 好看av亚洲va欧美ⅴa在| 久久精品夜夜夜夜夜久久蜜豆 | 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 麻豆成人av在线观看| 黄色视频,在线免费观看| 成人18禁高潮啪啪吃奶动态图| 熟妇人妻久久中文字幕3abv| 国产精品免费一区二区三区在线| 国产v大片淫在线免费观看| 久久欧美精品欧美久久欧美| 禁无遮挡网站| 丁香六月欧美| 亚洲精品久久成人aⅴ小说| 久久这里只有精品19| 在线观看免费视频日本深夜| 舔av片在线| 亚洲欧美日韩高清在线视频| 色噜噜av男人的天堂激情| 老司机午夜十八禁免费视频| 亚洲国产高清在线一区二区三| 18禁美女被吸乳视频| 国产99久久九九免费精品| 亚洲美女视频黄频| 两人在一起打扑克的视频| 精品一区二区三区视频在线观看免费| 日日干狠狠操夜夜爽| 中文字幕高清在线视频| 99热只有精品国产| 精品欧美一区二区三区在线| ponron亚洲| 日韩大码丰满熟妇| 亚洲最大成人中文| 国产一区二区三区在线臀色熟女| 国产精品一区二区免费欧美| 婷婷六月久久综合丁香| 中出人妻视频一区二区| 日韩精品青青久久久久久| 男人的好看免费观看在线视频 | 成在线人永久免费视频| 天天一区二区日本电影三级| 国模一区二区三区四区视频 | 深夜精品福利| 真人做人爱边吃奶动态| 丝袜美腿诱惑在线| 国语自产精品视频在线第100页| 亚洲精品久久国产高清桃花| 亚洲成人久久性| 中文字幕精品亚洲无线码一区| 国产亚洲精品第一综合不卡| 欧美中文综合在线视频| 在线免费观看的www视频| 国产99久久九九免费精品| 90打野战视频偷拍视频| www日本在线高清视频| 最近最新免费中文字幕在线| 成人特级黄色片久久久久久久| 午夜精品久久久久久毛片777| 波多野结衣高清无吗| 国产区一区二久久| 男女下面进入的视频免费午夜| 中文字幕最新亚洲高清| 19禁男女啪啪无遮挡网站| 国产亚洲精品av在线| 欧美成人性av电影在线观看| 在线a可以看的网站| 久久性视频一级片| 午夜福利在线在线| 国产激情久久老熟女| 中亚洲国语对白在线视频| 亚洲自偷自拍图片 自拍| 老鸭窝网址在线观看| 女同久久另类99精品国产91| 麻豆国产av国片精品| 人人妻,人人澡人人爽秒播| 精品午夜福利视频在线观看一区| 精品国产美女av久久久久小说| 久久久精品大字幕| 露出奶头的视频| www日本在线高清视频| 在线观看午夜福利视频| 男人舔女人的私密视频| 特大巨黑吊av在线直播| 美女高潮喷水抽搐中文字幕| 久久久久久久久中文| 日日干狠狠操夜夜爽| 亚洲人成网站在线播放欧美日韩| 国产真实乱freesex| 村上凉子中文字幕在线| 伦理电影免费视频| 久久久久久久午夜电影| 色老头精品视频在线观看| 桃红色精品国产亚洲av| 久久久久亚洲av毛片大全| 精品久久蜜臀av无| 日韩欧美三级三区| 韩国av一区二区三区四区| 国产精品久久久久久人妻精品电影| 成人特级黄色片久久久久久久| 黄色片一级片一级黄色片| 亚洲精品一区av在线观看| 欧美中文综合在线视频| 床上黄色一级片| 国产免费男女视频| 国产精华一区二区三区| 人人妻人人澡欧美一区二区| 国产精品 欧美亚洲| 男女那种视频在线观看| 舔av片在线| 国语自产精品视频在线第100页| 后天国语完整版免费观看| 男女午夜视频在线观看| 国产人伦9x9x在线观看| or卡值多少钱| 三级男女做爰猛烈吃奶摸视频| 99精品久久久久人妻精品| 国产在线观看jvid| 亚洲在线自拍视频| 日韩欧美免费精品| 日本黄色视频三级网站网址| 国内精品久久久久久久电影| 国产精品av视频在线免费观看| 脱女人内裤的视频| 两个人视频免费观看高清| 全区人妻精品视频| 桃色一区二区三区在线观看| 91麻豆精品激情在线观看国产| 中文字幕高清在线视频| 一级片免费观看大全| 老熟妇仑乱视频hdxx| 欧美成人免费av一区二区三区| 国产精品一及| 韩国av一区二区三区四区| 欧美成人一区二区免费高清观看 | 精品久久久久久久毛片微露脸| 亚洲午夜精品一区,二区,三区| 1024手机看黄色片| 啦啦啦免费观看视频1| 99热这里只有是精品50| 欧美日韩乱码在线| 欧美三级亚洲精品| 国产成人欧美在线观看| 日本精品一区二区三区蜜桃| 又黄又粗又硬又大视频| 久久午夜亚洲精品久久| 美女免费视频网站| 国产成人av激情在线播放| 欧美日韩乱码在线| 日本黄大片高清| 可以免费在线观看a视频的电影网站| 国产91精品成人一区二区三区| 亚洲成av人片在线播放无| 日本 欧美在线| 精品一区二区三区视频在线观看免费| 制服诱惑二区| 亚洲一区中文字幕在线| 成人av一区二区三区在线看| 久9热在线精品视频| 日韩欧美国产一区二区入口| 女警被强在线播放| 亚洲成人久久性| 成人av在线播放网站| 女同久久另类99精品国产91| 国产三级黄色录像| 免费观看精品视频网站| 色综合亚洲欧美另类图片| 亚洲人与动物交配视频| 日韩中文字幕欧美一区二区| 国产精品 国内视频| 日本熟妇午夜| 亚洲av熟女| 欧美不卡视频在线免费观看 | 淫妇啪啪啪对白视频| 亚洲精品美女久久av网站| 黄色毛片三级朝国网站| 亚洲狠狠婷婷综合久久图片| 国产午夜精品论理片| 亚洲七黄色美女视频| 手机成人av网站| 在线观看午夜福利视频| 国产99白浆流出| 国产麻豆成人av免费视频| 国产三级在线视频| 熟女少妇亚洲综合色aaa.| 丰满的人妻完整版| 国产成人av激情在线播放| 亚洲中文av在线| 国产亚洲精品一区二区www| 亚洲av第一区精品v没综合| av免费在线观看网站| 国产亚洲av嫩草精品影院| 国内少妇人妻偷人精品xxx网站 | 欧美大码av| 午夜亚洲福利在线播放| 悠悠久久av| 日本三级黄在线观看| 成在线人永久免费视频| 在线观看www视频免费| 人妻久久中文字幕网| av有码第一页| 国产视频一区二区在线看| 国产亚洲欧美98| 一二三四社区在线视频社区8| 欧美不卡视频在线免费观看 | 久久精品国产亚洲av香蕉五月| 国产亚洲精品综合一区在线观看 | 最好的美女福利视频网| 精品久久久久久久久久免费视频| 欧美zozozo另类| 亚洲av成人不卡在线观看播放网| 成人18禁高潮啪啪吃奶动态图| 精品免费久久久久久久清纯| 午夜福利高清视频| 亚洲欧洲精品一区二区精品久久久| 91av网站免费观看| 一区二区三区高清视频在线| 亚洲国产中文字幕在线视频| av在线天堂中文字幕| 麻豆久久精品国产亚洲av| 日韩欧美在线二视频| 中文字幕熟女人妻在线| 精品国产超薄肉色丝袜足j| 校园春色视频在线观看| 久久久精品欧美日韩精品| 999久久久精品免费观看国产| 精品久久久久久成人av| 一个人免费在线观看的高清视频| 欧美三级亚洲精品| 精品久久久久久成人av| 国产精品精品国产色婷婷| 色哟哟哟哟哟哟| svipshipincom国产片| 午夜亚洲福利在线播放| 日韩欧美 国产精品| 亚洲欧洲精品一区二区精品久久久| 国产黄色小视频在线观看| 又粗又爽又猛毛片免费看| 久久久久久久久免费视频了| 啪啪无遮挡十八禁网站|