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

    內(nèi)燃動(dòng)力包隔振參數(shù)靈敏度分析及優(yōu)化設(shè)計(jì)

    2022-01-07 08:53:50吳楊俊徐翠強(qiáng)陳杰賀小龍張立民
    關(guān)鍵詞:烈度靈敏度動(dòng)力

    吳楊俊,徐翠強(qiáng),陳杰,賀小龍,張立民

    (1. 西南交通大學(xué)牽引動(dòng)力國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川成都,610031;2. 中車青島四方機(jī)車車輛股份有限公司,山東青島,266111;3. 重慶文理學(xué)院智能制造工程學(xué)院,重慶,402160)

    目前,高速鐵路客運(yùn)廣泛采用電力動(dòng)車運(yùn)輸,但在鐵路電氣化水平不高的國(guó)家與地區(qū),為了提高鐵路運(yùn)輸能力,其鐵路客運(yùn)主要采用內(nèi)燃動(dòng)車運(yùn)輸。柴油發(fā)電機(jī)組為內(nèi)燃動(dòng)車的動(dòng)力源,由于車下安裝空間的限制,柴油發(fā)電機(jī)機(jī)組及其他附屬設(shè)備被安裝在1個(gè)基礎(chǔ)框架上,框架通過(guò)二級(jí)隔振器與車體相連,從而構(gòu)成了動(dòng)力包雙層隔振系統(tǒng)。動(dòng)力包工作狀態(tài)下產(chǎn)生的振動(dòng)既會(huì)影響動(dòng)力包內(nèi)部設(shè)備的振動(dòng)狀態(tài),也會(huì)通過(guò)隔振器傳遞到車體,降低乘客的乘坐舒適性,因此,動(dòng)力包隔振參數(shù)的優(yōu)化設(shè)計(jì)也是內(nèi)燃動(dòng)車組設(shè)計(jì)中不可或缺的一環(huán)。

    針對(duì)雙層隔振理論和技術(shù),許多學(xué)者從不同的角度進(jìn)行了大量研究[1?3]。目前,動(dòng)力包雙層隔振系統(tǒng)的隔振技術(shù)還處在初步探索和應(yīng)用階段,帶有源子隔振系統(tǒng)的雙層隔振系統(tǒng)隔振設(shè)計(jì)的相關(guān)研究還較少。GINA 等[4]研究了運(yùn)載火箭及其子系統(tǒng)的隔振設(shè)計(jì),但并沒有深入探討兩者間的耦合振動(dòng)特性。孫玉華等[5?7]針對(duì)內(nèi)燃動(dòng)車動(dòng)力包雙層隔振系統(tǒng),建立了動(dòng)力包雙層隔振系統(tǒng)有限元模型,將解耦率、隔振效率及振動(dòng)烈度作為評(píng)價(jià)指標(biāo),運(yùn)用枚舉法從多個(gè)設(shè)計(jì)方案中選擇出滿足工程要求的方案。此方法雖然取得了較好的效果,但設(shè)計(jì)效率較低、計(jì)算工作量較大,且不易獲得最優(yōu)參數(shù)方案。為此,時(shí)威振[8]以內(nèi)燃機(jī)車動(dòng)力總成為研究對(duì)象,分析了隔振系統(tǒng)懸掛參數(shù)對(duì)其隔振性能的影響規(guī)律,并選用基于Pareto最優(yōu)解的多目標(biāo)遺傳算法對(duì)系統(tǒng)進(jìn)行多目標(biāo)優(yōu)化設(shè)計(jì),從而提高了設(shè)計(jì)效率。陳俊等[9]研究了子系統(tǒng)對(duì)雙層隔振主系統(tǒng)固有特性的影響規(guī)律,并對(duì)子系統(tǒng)設(shè)計(jì)提出了合理建議。

    從以上研究成果可以看出,目前,針對(duì)內(nèi)燃動(dòng)力包雙層隔振系統(tǒng)設(shè)計(jì)并未形成統(tǒng)一的標(biāo)準(zhǔn)或者行業(yè)規(guī)范。由于動(dòng)力包隔振參數(shù)優(yōu)化設(shè)計(jì)涉及變量較多,導(dǎo)致優(yōu)化設(shè)計(jì)過(guò)程中計(jì)算量較大,設(shè)計(jì)效率較低。為解決上述問(wèn)題,本文將全局靈敏度分析方法引入動(dòng)力包隔振參數(shù)設(shè)計(jì)中,通過(guò)全局靈敏度方法確定對(duì)動(dòng)力包隔振性能影響較大的參數(shù),并將其作為優(yōu)化變量對(duì)動(dòng)力包隔振參數(shù)進(jìn)行優(yōu)化設(shè)計(jì),以期有效縮減設(shè)計(jì)中所需的優(yōu)化變量,提高優(yōu)化效率。

    1 內(nèi)燃動(dòng)力包系統(tǒng)數(shù)學(xué)建模

    以某型號(hào)內(nèi)燃動(dòng)力包雙層隔振系統(tǒng)(見圖1)為研究對(duì)象,建立內(nèi)燃動(dòng)力包數(shù)學(xué)模型,該模型包含1個(gè)框架(圖2(a))、1個(gè)柴油發(fā)電機(jī)組(圖1(b))和1個(gè)冷卻風(fēng)機(jī)(圖2(c))。動(dòng)力包模型含有11 個(gè)隔振器,其中包括機(jī)組與框架連接處的3個(gè)一級(jí)隔振器(布置在a1,a2和a3懸掛點(diǎn)),冷卻裝置與框架連接處的4 個(gè)一級(jí)隔振器(布置在a4,a5,a6和a7懸掛點(diǎn))以及框架與基礎(chǔ)連接處的4 個(gè)二級(jí)隔振器(布置在b1,b2,b3和b4懸掛點(diǎn))。

    圖1 動(dòng)力包結(jié)構(gòu)俯視圖Fig.1 Top view of power pack structure

    圖2 動(dòng)力包各部件平面圖Fig.2 Plane view of each part of power pack

    內(nèi)燃動(dòng)力包振動(dòng)模型的18 個(gè)自由度分別為:框架沿x,y和z方向平動(dòng)的自由度xo1,yo1和zo1;過(guò)框架質(zhì)心o1繞x,y和z方向的轉(zhuǎn)動(dòng)自由度αo1,βo1和γo1;柴油機(jī)沿x,y和z方向平動(dòng)的自由度xo1,yo2和zo2,過(guò)機(jī)組質(zhì)心o2繞x,y和z方向的轉(zhuǎn)動(dòng)自由度αo2,βo2和γo2;冷卻風(fēng)機(jī)沿x,y和z方向平動(dòng)的自由度xo3,yo3和zo3,過(guò)冷卻風(fēng)機(jī)質(zhì)心o3繞x,y和z軸的轉(zhuǎn)動(dòng)自由度αo3,βo3和γo3。動(dòng)力包物理參數(shù)見表1。

    常規(guī)的內(nèi)燃動(dòng)力包隔振系統(tǒng)設(shè)計(jì)一般將各個(gè)隔振器三向剛度作為優(yōu)化變量。為減少設(shè)計(jì)變量的數(shù)量,本文將各個(gè)隔振器垂向剛度、橫垂比(隔振器橫向剛度與垂向剛度比值)、縱垂比(隔振器縱向剛度與垂向剛度比值)作為設(shè)計(jì)參數(shù),同時(shí),動(dòng)力包各參數(shù)計(jì)算值由相關(guān)合作單位提供,見表1。

    表1 動(dòng)力包物理參數(shù)Table 1 Physical parameter of power pack

    根據(jù)振動(dòng)理論可得出柴油發(fā)動(dòng)機(jī)沿x,y和z軸平動(dòng)自由度的振動(dòng)方程為

    柴油發(fā)動(dòng)機(jī)繞x,y,z軸轉(zhuǎn)動(dòng)自由度的振動(dòng)方程為

    式中:Faix,F(xiàn)aiy和Faiz分別為ai隔振器在x,y和z方向的作用力;Fo2x,F(xiàn)o2y和Fo2z分別為作用在柴油機(jī)x,y和z軸上的激振力;Mo2x,Mo2y和Mo2z分別為作用在柴油機(jī)x,y和z軸上的激振力偶。

    冷卻風(fēng)機(jī)沿x,y和z軸平動(dòng)自由度的振動(dòng)方程為

    冷卻風(fēng)機(jī)繞x,y和z軸轉(zhuǎn)動(dòng)自由度的振動(dòng)方程為

    式中:Fo3x,F(xiàn)o3y和Fo3z分別為作用在冷卻風(fēng)機(jī)x,y和z軸上的激振力;Mo3x,Mo3y和Mo3z分別為作用在冷卻風(fēng)機(jī)x,y和z軸上的激振力偶。

    框架沿x,y和z軸平動(dòng)自由度的振動(dòng)方程為

    框架繞x,y和z軸轉(zhuǎn)動(dòng)自由度的振動(dòng)方程為

    式中:Fbix,F(xiàn)biy和Fbiz分別為bi隔振器在x,y和z方向上作用力。

    2 動(dòng)力包系統(tǒng)振動(dòng)響應(yīng)及隔振性能指標(biāo)的計(jì)算方法

    2.1 動(dòng)力包系統(tǒng)數(shù)值求解方法

    本文采用翟婉明[10]提出的新型快速顯示數(shù)值積分法求解車輛系統(tǒng)的振動(dòng)響應(yīng),假設(shè)系統(tǒng)在t=(h+ 1)Δt瞬時(shí)的振動(dòng)方程為

    式中:M為系統(tǒng)的質(zhì)量矩陣;Kh+1和Ch+1分別為第h+1 步迭代系統(tǒng)的剛度與阻尼矩陣;Ph+1為第h+1步迭代系統(tǒng)的廣義載荷矢量;為第h+1步迭代系統(tǒng)的廣義加速度矢量;為第h+1步迭代系統(tǒng)的廣義速度矢量;Xh+1為第h+1步迭代系統(tǒng)的廣義位移矢量。Δt為時(shí)間積分步長(zhǎng),φ與φ為控制參數(shù)。

    系統(tǒng)初始條件為

    根據(jù)式(7)~(9)逐次求出各迭代步的位移、速度與加速度離散值。起步時(shí)只需令φ=φ= 0,則可使本方法具有積分“自開始”的特性。

    2.2 動(dòng)力包機(jī)組振動(dòng)烈度求解方法

    機(jī)組振動(dòng)烈度反映了柴油發(fā)電機(jī)組自身振動(dòng)環(huán)境,若機(jī)組振動(dòng)烈度過(guò)大,則會(huì)嚴(yán)重影響機(jī)組工作質(zhì)量并縮短其使用壽命[8],因此,可將動(dòng)力包機(jī)組振動(dòng)烈度作為系統(tǒng)隔振性能評(píng)價(jià)指標(biāo)之一。參照TB/T 3164—2007“柴油機(jī)車車內(nèi)設(shè)備機(jī)械振動(dòng)烈度評(píng)定方法”[11]中柴油發(fā)電機(jī)組測(cè)點(diǎn)布置規(guī)定,結(jié)合該柴油發(fā)電機(jī)組結(jié)構(gòu)特點(diǎn),機(jī)組振動(dòng)烈度的測(cè)點(diǎn)位置Du(u=1~8)分布如圖3所示。

    圖3 振動(dòng)烈度考核點(diǎn)位置分布Fig.3 Location distribution of vibration intensity measurement points

    基于動(dòng)力包機(jī)組振動(dòng)烈度測(cè)點(diǎn)位置,可得機(jī)組振動(dòng)烈度Vrms計(jì)算公式:

    式中:VrxDu,VryDu和VrzDu分別為測(cè)點(diǎn)Du(u=1~8)在x,y和z方向的均方根速度。

    2.3 動(dòng)力包系統(tǒng)隔振效率求解方法

    動(dòng)力包雙層隔振系統(tǒng)設(shè)計(jì)的目的之一是最大限度地減少內(nèi)部設(shè)備激勵(lì)的傳遞,避免動(dòng)力包與車體之間的振動(dòng)耦合,提高整車的乘坐舒適性[5]。動(dòng)力包系統(tǒng)隔振效率能夠反映系統(tǒng)激勵(lì)傳遞到車體上的衰減情況,可作為動(dòng)力包隔振性能評(píng)價(jià)指標(biāo)。根據(jù)文獻(xiàn)[5]可知內(nèi)燃動(dòng)力包系統(tǒng)傳遞到基座上的當(dāng)量力F為

    式中:Frbix,F(xiàn)rbiy與Frbiz(i=1~4)分別為bi處隔振器在x,y和z方向上的力的均方根值;Nx,Ny,Nz分別為各隔振器在x,y,z這3個(gè)方向上的測(cè)點(diǎn)個(gè)數(shù)。

    由于動(dòng)力包內(nèi)部激勵(lì)既含有激振力又含有激振力矩,因此,在計(jì)算動(dòng)力包隔振效率之前,需將內(nèi)部激勵(lì)轉(zhuǎn)化為基座上4個(gè)隔振器的三向力,即動(dòng)力包內(nèi)部激勵(lì)與4個(gè)隔振器上的三向力處于平衡狀態(tài)。根據(jù)超靜定方程可獲得內(nèi)部激勵(lì)轉(zhuǎn)化到基座上4個(gè)隔振器的三向力,進(jìn)而得到轉(zhuǎn)化后的激勵(lì)當(dāng)量力F0為

    式中:F′rbix,F(xiàn)′rbiy與F′rbiz(i=1~4)分別為內(nèi)部激勵(lì)轉(zhuǎn)化到bi處隔振器的力在x,y和z方向上的均方根值;

    動(dòng)力包系統(tǒng)的隔振效率η為

    3 動(dòng)力包隔振參數(shù)全局靈敏度分析

    為研究動(dòng)力包各隔振參數(shù)對(duì)動(dòng)力包隔振性能影響的程度,進(jìn)而確定影響系統(tǒng)隔振性能的主要參數(shù),需對(duì)系統(tǒng)隔振參數(shù)進(jìn)行靈敏度分析。早期主要使用的是局部靈敏度分析法,此類方法概念明確,計(jì)算方便,但鑒于其是以微分或差分理論為基礎(chǔ),因此,系統(tǒng)結(jié)構(gòu)參數(shù)的變動(dòng)范圍不能過(guò)大。當(dāng)系統(tǒng)非線性特征較強(qiáng)或者參數(shù)變化范圍較大時(shí),局部靈敏度分析方法往往不能得出有效結(jié)果[12]。為克服局部靈敏度分析法所存在的缺陷,全局靈敏度分析法應(yīng)運(yùn)而生[13]。Sobol 法[14]是一種基于方差的全局靈敏度分析法,該方法能夠快速簡(jiǎn)便地計(jì)算出高階交叉影響項(xiàng),目前已被廣泛應(yīng)用于經(jīng)濟(jì)、環(huán)境、生物、物理、化學(xué)、控制及神經(jīng)網(wǎng)絡(luò)等領(lǐng)域的研究。

    3.1 Sobol靈敏度分析法原理

    定義單位區(qū)間I為[0,1],In為n維超立方單元體,假設(shè)函數(shù)f(g)的變量g=(g1,g2,g3,…,gn),g∈In,可將函數(shù)f(g) 看作是2n個(gè)遞增子項(xiàng)之和[14]:

    式中:f0為期望。

    假設(shè)q={g1,g2,…,gm}?{g1,g2,…,gn},同時(shí)Mq={1,2,…,m}, 其補(bǔ)集v={g1,g2,…,gn}q,Mv={1,2,…,n}Mq,q的一階和高階靈敏度表達(dá)式為式中:D為函數(shù)f(g)的總方差;Dv為v的偏方差;Dq為q的偏方差;Sv為v的一階靈敏度;Sq和S分別為參數(shù)q的一階靈敏度與高階靈敏度,其中高階靈敏度也可稱為總靈敏度,0 ≤Sq≤≤1。當(dāng)Sq== 1 時(shí),f(g) 只與q有關(guān);當(dāng)Sq== 0時(shí),f(g)與q無(wú)關(guān)。

    按式(15)直接計(jì)算靈敏度會(huì)碰到很多困難,因此,Sobol 靈敏度可運(yùn)用蒙特卡羅積分獲得,計(jì)算公式如下:

    式中:gi=(qi,vi)和g′i=(q′i,v′i)為2 組樣本數(shù)據(jù);i= 1,2,…,N;N為采樣點(diǎn)數(shù)。

    將式(16)中的Dv,D和Dq代入式(15)可得到相應(yīng)q的一階和高階靈敏度。

    參數(shù)的一階靈敏度只反映了此參數(shù)單獨(dú)變化時(shí)對(duì)目標(biāo)函數(shù)的影響程度,參數(shù)的高靈敏度不僅反映了該參數(shù)單獨(dú)變化的影響,也反映了該參數(shù)與其他所有參數(shù)的交互作用對(duì)目標(biāo)函數(shù)的影響。

    3.2 動(dòng)力包系統(tǒng)隔振參數(shù)靈敏度結(jié)果分析

    采用Sobol法分析各隔振參數(shù)對(duì)動(dòng)力包隔振性能的影響程度時(shí),需給定參數(shù)的變化區(qū)間。各參數(shù)的取值范圍見表2。

    表2 動(dòng)力包隔振參數(shù)取值范圍Table 2 Value range of vibration isolation parameter of power pack

    基于動(dòng)力包18自由度數(shù)學(xué)模型,運(yùn)用Sobol法計(jì)算得到不同轉(zhuǎn)速工況下動(dòng)力包各隔振參數(shù)對(duì)2個(gè)隔振性能指標(biāo)的一階與高階靈敏度。

    圖4和圖5所示分別為各個(gè)轉(zhuǎn)速工況下動(dòng)力包各隔振參數(shù)對(duì)機(jī)組振動(dòng)烈度的一階與高階靈敏度分布。圖中工況1 對(duì)應(yīng)轉(zhuǎn)速為1 000 r/min(空載),工況2~10 分別對(duì)應(yīng)轉(zhuǎn)速為1 000,1 100,…,1 800 r/min(負(fù)載)。

    從圖5可以看出:隨著動(dòng)力包轉(zhuǎn)速增加,nV與nH的高階靈敏度整體上呈增大趨勢(shì)。此外,從圖4和圖5可以看出,一階靈敏度與高階靈敏度所反映的情況不完全相同,例如,當(dāng)動(dòng)力包轉(zhuǎn)速工況為1 700 r/min(負(fù)載)時(shí),nV的一階靈敏度比nH的大,而nV的高階靈敏度比nH的小,這是因?yàn)橐浑A靈敏度并不能反映評(píng)價(jià)參數(shù)與其他參數(shù)的交互作用對(duì)目標(biāo)函數(shù)的影響。從圖4和圖5看出,對(duì)機(jī)組振動(dòng)烈度而言,一階與高階靈敏度較大的參數(shù)為Kzda1,Kzda2,Kzda3,Kzdb1,Kzdb2,Kzdb3,Kzdb4,nV和nH。

    圖4 隔振參數(shù)對(duì)機(jī)組振動(dòng)烈度的一階靈敏度Fig.4 The first-order sensitivity of vibration isolation parameters to unit’s vibration intensity

    圖5 隔振參數(shù)對(duì)機(jī)組振動(dòng)烈度的高階靈敏度Fig.5 High-order sensitivity of vibration isolation parameters to unit’s vibration intensity

    圖6和圖7所示分別為各個(gè)轉(zhuǎn)速工況下動(dòng)力包各隔振參數(shù)對(duì)系統(tǒng)隔振效率的一階與高階靈敏度分布情況。

    圖6 隔振參數(shù)對(duì)系統(tǒng)隔振效率的一階靈敏度Fig.6 The first-order sensitivity of vibration isolation parameters to vibration isolation efficiency

    圖7 隔振參數(shù)對(duì)系統(tǒng)隔振效率的高階靈敏度Fig.7 High-order sensitivity of vibration isolation parameters to vibration isolation efficiency

    從圖6和圖7可以看出,在1 000 r/min(空載)與1 000 r/min(負(fù)載)工況下,對(duì)系統(tǒng)隔振效率而言,一階靈敏度最大的參數(shù)為Kzdb4,而高階靈敏度最大的參數(shù)為Kzda3,可見一階靈敏度對(duì)參數(shù)影響程度的評(píng)價(jià)結(jié)果存在誤差。從圖6和圖7還可以看出,對(duì)于系統(tǒng)隔振效率而言,一階與高階敏度較大的參數(shù)為Kzda1,Kzda2,Kzda3,Kzdb1,Kzdb2,Kzdb3,Kzdb4,nV和nH。

    內(nèi)燃動(dòng)力包系統(tǒng)隔振參數(shù)全局靈敏度分析結(jié)果表明,在所有轉(zhuǎn)速工況下,冷卻風(fēng)機(jī)上4個(gè)一級(jí)隔振器垂向剛度對(duì)系統(tǒng)2個(gè)隔振性能指標(biāo)的靈敏度較小,而其他9個(gè)參數(shù)對(duì)系統(tǒng)隔振性能指標(biāo)的靈敏度較大,因此,可將這9個(gè)參數(shù)作為影響動(dòng)力包系統(tǒng)隔振性能的主要參數(shù)。

    4 動(dòng)力包隔振參數(shù)優(yōu)化設(shè)計(jì)

    4.1 優(yōu)化變量

    根據(jù)第3節(jié)動(dòng)力包隔振參數(shù)靈敏度分析結(jié)果可知,影響機(jī)組振動(dòng)烈度與系統(tǒng)隔振效率的主要參數(shù)為Kzda1,Kzda2,Kzda3,Kzdb1,Kzdb2,Kzdb3,Kzdb4,nV和nH,因此,本文將橫垂比、縱垂比、機(jī)組a1,a2和a3以及框架b1,b2,b3和b4懸掛點(diǎn)上的隔振器的垂向懸掛剛度作為優(yōu)化變量,優(yōu)化變量U可以表示為

    4.2 優(yōu)化目標(biāo)函數(shù)的建立

    將動(dòng)力包機(jī)組振動(dòng)烈度作為優(yōu)化目標(biāo)之一,基于本文2.2節(jié)中振動(dòng)烈度計(jì)算公式,可得到機(jī)組振動(dòng)烈度目標(biāo)函數(shù)J1為

    同理,動(dòng)力包系統(tǒng)隔振效率也可作為優(yōu)化目標(biāo)之一。同時(shí),為使2個(gè)目標(biāo)函數(shù)都以最小值為最優(yōu),將目標(biāo)隔振效率轉(zhuǎn)化為力傳遞率,基于本文2.3 節(jié)中隔振效率計(jì)算公式,可得到系統(tǒng)力傳遞率目標(biāo)函數(shù)J2為

    為滿足工程需求,系統(tǒng)的隔振效率需大于80%,即力傳遞率要小于20%。同時(shí)機(jī)組振動(dòng)烈度應(yīng)該為B 級(jí)或A 級(jí)(Vrms≤0.018 m/s)。為便于計(jì)算,需運(yùn)用加權(quán)系數(shù)法將多目標(biāo)函數(shù)問(wèn)題轉(zhuǎn)化為單目標(biāo)函數(shù)問(wèn)題,由于2個(gè)目標(biāo)函數(shù)存在差別,需對(duì)目標(biāo)函數(shù)進(jìn)行歸一化處理:

    式中:λ1和λ2為加權(quán)因子,且有

    如果沒強(qiáng)調(diào)以某一目標(biāo)為主要優(yōu)化目標(biāo),那么2個(gè)加權(quán)因子均取0.5。

    4.3 約束條件

    4.3.1 模態(tài)頻率匹配約束

    根據(jù)隔振理論可知,為避免隔振系統(tǒng)激勵(lì)頻率與模態(tài)頻率出現(xiàn)共振現(xiàn)象,同時(shí)保證系統(tǒng)具有良好的隔振效果,系統(tǒng)的激勵(lì)頻率與模態(tài)頻率的比值應(yīng)該大于,若工程要求激勵(lì)頻率必須小于模態(tài)頻率時(shí),則激勵(lì)頻率與模態(tài)頻率的比值應(yīng)該小于。其約束條件為

    式中:ωj為第j個(gè)激勵(lì)頻率,ωi0為系統(tǒng)第i階模態(tài)頻率。

    4.3.2 靜平衡條件約束

    為防止動(dòng)力包發(fā)生傾覆,要求裝車后機(jī)組與安裝框架不產(chǎn)生較大傾斜,因此,要求相同部件上的同級(jí)隔振器垂向靜壓縮量差不超過(guò)1 mm[15];此外,由于車下安裝空間限制,還需將一級(jí)隔振器的垂向靜壓縮量范圍控制在2~7 mm,并且二級(jí)隔振器的垂向靜壓縮量范圍控制在3~9 mm。其約束條件如下:

    式中:Zbl和Zae分別為懸掛點(diǎn)bl和ae上隔振器的垂向靜壓縮量,其中l(wèi)取值范圍為1~4,e取值范圍為1~7。

    本文采用的隔振器為橡膠隔振器,其動(dòng)靜比(動(dòng)剛度與靜剛度的比值)取值范圍一般為1.3~1.6[16],本文橡膠隔振器動(dòng)靜比取為1.5。通過(guò)上述靜平衡約束可對(duì)各隔振器的靜剛度進(jìn)行約束,進(jìn)而結(jié)合動(dòng)靜比對(duì)各隔振器的動(dòng)剛度進(jìn)行約束。

    4.4 動(dòng)力包隔振參數(shù)優(yōu)化分析

    基于本文所建內(nèi)燃動(dòng)力包優(yōu)化模型,在轉(zhuǎn)速為1 800 r/min(負(fù)載)激勵(lì)條件下,運(yùn)用遺傳算法對(duì)動(dòng)力包系統(tǒng)進(jìn)行優(yōu)化設(shè)計(jì)。

    圖8所示為各代種群中所有個(gè)體的最佳適應(yīng)度值與平均適應(yīng)度的變化曲線,種群中所有個(gè)體的最佳適應(yīng)度隨種群迭代數(shù)不斷下降,同時(shí)可以看出45代之后,最佳適應(yīng)度趨于穩(wěn)定,直到75代時(shí)滿足終止條件,其種群最佳適應(yīng)度為0.538,因此,可認(rèn)為遺傳算法優(yōu)化過(guò)程達(dá)到了收斂狀態(tài)。

    圖8 適應(yīng)度曲線Fig.8 Fitness curves

    表3 所示為內(nèi)燃動(dòng)力包參數(shù)優(yōu)化前后的取值。圖9 和圖10 分別為不同轉(zhuǎn)速工況下內(nèi)燃動(dòng)力包中柴油發(fā)電機(jī)組振動(dòng)烈度與系統(tǒng)隔振效率優(yōu)化前后結(jié)果對(duì)比。圖11 所示為系統(tǒng)各隔振性能指標(biāo)在不同轉(zhuǎn)速工況下的變化幅度。

    表3 優(yōu)化前后動(dòng)力包隔振參數(shù)Table 3 Vibration isolation parameters of power pack before and after optimization

    從圖9可以看出,參數(shù)優(yōu)化后,各工況下的機(jī)組振動(dòng)烈度明顯降低且都小于0.018 m/s。

    圖9 機(jī)組振動(dòng)烈度優(yōu)化前后對(duì)比Fig.9 Comparison of unit's vibration intensity before and after optimization

    從圖10 可見,參數(shù)優(yōu)化后,系統(tǒng)的隔振效率得到了顯著提升且各工況下的隔振效率都達(dá)到了85%以上。

    圖10 隔振效率優(yōu)化前后對(duì)比Fig.10 Comparison of vibration isolation efficiency before and after optimization

    從圖11 可知,不同工況下機(jī)組振動(dòng)烈度與系統(tǒng)隔振效率的變化幅度趨勢(shì)相近,且都在1 000 r/min(空載)工況時(shí)變化幅度最大,分別達(dá)到了12.00%與29.61%。

    圖11 隔振性能指標(biāo)變化幅度Fig.11 Change amplitudes of index of vibration isolation performance

    綜上所述,經(jīng)過(guò)優(yōu)化設(shè)計(jì)后的內(nèi)燃動(dòng)力包隔振系統(tǒng)的振動(dòng)烈度明顯下降,系統(tǒng)隔振效率明顯提升,說(shuō)明此優(yōu)化方法是可行的。

    5 動(dòng)力包隔振性能試驗(yàn)與仿真結(jié)果對(duì)比分析

    將動(dòng)力包安裝于地面臺(tái)架上(見圖12),對(duì)各工況下動(dòng)力包隔振性能進(jìn)行測(cè)試,其中各隔振器剛度為本文第4節(jié)所提優(yōu)化設(shè)計(jì)的剛度。

    圖12 動(dòng)力包臺(tái)架Fig.12 Platform of power pack

    5.1 機(jī)組振動(dòng)烈度試驗(yàn)與仿真對(duì)比分析

    圖13所示為振動(dòng)烈度測(cè)點(diǎn)D5在各轉(zhuǎn)速工況下的三向加速度振動(dòng)信號(hào)。

    圖13 D5點(diǎn)三向加速度振動(dòng)信號(hào)Fig.13 Three-directional acceleration vibration signal of measuring point D5

    實(shí)驗(yàn)測(cè)試結(jié)束后,獲取各測(cè)試點(diǎn)的加速度響應(yīng)信號(hào),并通過(guò)積分將加速度信號(hào)轉(zhuǎn)變?yōu)闇y(cè)點(diǎn)速度信號(hào),分別求出各個(gè)測(cè)點(diǎn)在測(cè)試時(shí)間段內(nèi)的三向速度,從而得到動(dòng)力包機(jī)組振動(dòng)烈度。

    圖14 所示為柴油發(fā)電機(jī)組振動(dòng)烈度實(shí)驗(yàn)值與仿真值對(duì)比。從圖14 可以看出:振動(dòng)烈度實(shí)驗(yàn)與仿真結(jié)果相近,且變化趨勢(shì)一致。由于仿真模型將動(dòng)力包隔振系統(tǒng)考慮成純剛體模型,而實(shí)驗(yàn)中的動(dòng)力包系統(tǒng)是柔性模型,柔性構(gòu)架在動(dòng)力包運(yùn)行過(guò)程中會(huì)產(chǎn)生一定的彈性振動(dòng),導(dǎo)致仿真結(jié)果與實(shí)驗(yàn)結(jié)果產(chǎn)生了一定偏差。實(shí)驗(yàn)結(jié)果表明,動(dòng)力包在1 000 r/min空載與負(fù)載2個(gè)工況下機(jī)組振動(dòng)烈度的等級(jí)達(dá)到A級(jí),其他工況下都達(dá)到B級(jí)。測(cè)試結(jié)果表明,本文設(shè)計(jì)的動(dòng)力包結(jié)構(gòu)機(jī)組振動(dòng)烈度滿足工程要求。

    圖14 機(jī)組振動(dòng)烈度實(shí)驗(yàn)值與仿真值對(duì)比Fig.14 Comparison of experiment values and simulation values of unit's vibration intensity

    5.2 動(dòng)力包隔振效率試驗(yàn)與仿真對(duì)比分析

    分別在4個(gè)二級(jí)隔振器b1,b2,b3和b4的上下2 個(gè)位置各布置1 個(gè)測(cè)點(diǎn),從而測(cè)得隔振器上下測(cè)點(diǎn)加速度時(shí)域信號(hào),再通過(guò)二次積分得到相應(yīng)測(cè)點(diǎn)的位移時(shí)域信號(hào),接著將各個(gè)二級(jí)隔振器上下測(cè)點(diǎn)的位移差與對(duì)應(yīng)隔振器的剛度相乘即可以求出傳遞到基礎(chǔ)框架上的力時(shí)域信號(hào),從而獲得本次地面臺(tái)架實(shí)驗(yàn)動(dòng)力包系統(tǒng)的隔振效率。

    圖15 所示為動(dòng)力包隔振系統(tǒng)的隔振效率實(shí)驗(yàn)值與仿真值對(duì)比。從圖15 可以看出:隔振效率仿真結(jié)果與實(shí)驗(yàn)結(jié)果偏差較小,并且變化趨勢(shì)相近;各工況下動(dòng)力包系統(tǒng)隔振效率都在85%以上;當(dāng)轉(zhuǎn)速在1 000 r/min(負(fù)載)以上時(shí),系統(tǒng)隔振效率均在90%以上。實(shí)驗(yàn)結(jié)果表明,本文設(shè)計(jì)的動(dòng)力包雙層隔振系統(tǒng)在各工況下都具有良好的隔振效果。

    圖15 機(jī)組隔振效率實(shí)驗(yàn)值與仿真值對(duì)比Fig.15 Comparison of experiment and simulation values of vibration isolation efficiency

    6 結(jié)論

    1)對(duì)機(jī)組振動(dòng)烈度與系統(tǒng)隔振效率而言,靈敏度較大的參數(shù)變量為Kzda1,Kzda2,Kzda3,Kzdb1,Kzdb2,Kzdb3,Kzdb4,nV和nH。因此,上述參數(shù)可作為動(dòng)力包系統(tǒng)隔振性能的主要影響參數(shù)。

    2)與優(yōu)化前相比,優(yōu)化后的動(dòng)力包系統(tǒng)各個(gè)工況下的隔振性能得到了顯著改善,在1 000 r/min(空載)工況下,機(jī)組振動(dòng)烈度降低了12.00%,系統(tǒng)的隔振效率提升了29.61%。

    3) 所有工況下動(dòng)力包機(jī)組的振動(dòng)烈度都為B級(jí)或A 級(jí),且系統(tǒng)隔振效率都在85%以上,滿足工程要求。同時(shí),動(dòng)力包2個(gè)隔振指標(biāo)的實(shí)驗(yàn)結(jié)果與仿真結(jié)果差異較小,且變化趨勢(shì)一致,驗(yàn)證了內(nèi)燃動(dòng)力包隔振系統(tǒng)數(shù)學(xué)模型的準(zhǔn)確性及系統(tǒng)隔振參數(shù)優(yōu)化設(shè)計(jì)方法的可行性。

    猜你喜歡
    烈度靈敏度動(dòng)力
    學(xué)習(xí)動(dòng)力不足如何自給自足
    高烈度區(qū)域深基坑基坑支護(hù)設(shè)計(jì)
    高烈度地震區(qū)非規(guī)則多跨長(zhǎng)聯(lián)連續(xù)梁抗震分析
    導(dǎo)磁環(huán)對(duì)LVDT線性度和靈敏度的影響
    胖胖一家和瘦瘦一家(10)
    地下水非穩(wěn)定流的靈敏度分析
    動(dòng)力船
    穿甲爆破彈引信對(duì)薄弱目標(biāo)的靈敏度分析
    318國(guó)道沿線蘆山地震的震害特征與烈度區(qū)劃探討
    基于多動(dòng)力總成的六點(diǎn)懸置匹配計(jì)算
    亚洲精品美女久久久久99蜜臀 | 99久久综合免费| 色综合欧美亚洲国产小说| 国产麻豆69| 桃花免费在线播放| 日韩大码丰满熟妇| 母亲3免费完整高清在线观看| 一级毛片黄色毛片免费观看视频| 亚洲国产精品一区二区三区在线| 欧美激情高清一区二区三区 | 一区二区三区乱码不卡18| 在现免费观看毛片| 老熟女久久久| 女性生殖器流出的白浆| 日韩不卡一区二区三区视频在线| 久久精品久久精品一区二区三区| 久久久久久久大尺度免费视频| 久久毛片免费看一区二区三区| 免费观看av网站的网址| 成人国产av品久久久| 精品人妻一区二区三区麻豆| 国产免费现黄频在线看| av福利片在线| 黄色毛片三级朝国网站| 亚洲成人一二三区av| 国产亚洲一区二区精品| av线在线观看网站| 精品人妻一区二区三区麻豆| 成人午夜精彩视频在线观看| 日本黄色日本黄色录像| 亚洲av电影在线进入| 亚洲av电影在线进入| 美女脱内裤让男人舔精品视频| 老司机亚洲免费影院| 在线观看一区二区三区激情| 黑人巨大精品欧美一区二区蜜桃| 亚洲成国产人片在线观看| 亚洲av国产av综合av卡| 又粗又硬又长又爽又黄的视频| 中国国产av一级| 男女高潮啪啪啪动态图| 亚洲精品av麻豆狂野| 老司机影院成人| 18禁动态无遮挡网站| 成年动漫av网址| 国产欧美亚洲国产| 婷婷色综合大香蕉| 国产免费视频播放在线视频| 国产无遮挡羞羞视频在线观看| 91成人精品电影| 伦理电影大哥的女人| 国产一区二区三区av在线| 男女高潮啪啪啪动态图| 99热全是精品| 欧美人与善性xxx| 视频在线观看一区二区三区| 男女高潮啪啪啪动态图| 老司机在亚洲福利影院| 搡老岳熟女国产| 久久久国产精品麻豆| 中文欧美无线码| 久久久久精品国产欧美久久久 | 国产av码专区亚洲av| 日韩伦理黄色片| 久久热在线av| 久久久久国产精品人妻一区二区| 精品一品国产午夜福利视频| 麻豆av在线久日| 午夜福利视频精品| 国产精品蜜桃在线观看| 久久久久久久久久久免费av| 免费不卡黄色视频| 亚洲精品国产av成人精品| 99久久人妻综合| 成年动漫av网址| 亚洲成av片中文字幕在线观看| 亚洲人成77777在线视频| 精品国产乱码久久久久久男人| 国产麻豆69| 久久精品久久精品一区二区三区| 成年av动漫网址| 亚洲精品美女久久av网站| 中文欧美无线码| 久久综合国产亚洲精品| 免费观看性生交大片5| 午夜福利视频精品| 中国国产av一级| 日韩精品免费视频一区二区三区| 一级黄片播放器| 热re99久久精品国产66热6| 久久久国产一区二区| 电影成人av| 在线天堂最新版资源| 久久久久久久久久久免费av| 久久精品国产亚洲av涩爱| 波多野结衣av一区二区av| 国产精品久久久久久精品电影小说| av福利片在线| 少妇的丰满在线观看| av一本久久久久| 久久人人爽av亚洲精品天堂| 国产精品一区二区在线不卡| 一边亲一边摸免费视频| 精品国产乱码久久久久久小说| 国产精品成人在线| 黄片无遮挡物在线观看| 亚洲av欧美aⅴ国产| 色婷婷久久久亚洲欧美| 少妇人妻精品综合一区二区| 热re99久久精品国产66热6| 韩国精品一区二区三区| 熟妇人妻不卡中文字幕| 美国免费a级毛片| 一本久久精品| 国产色婷婷99| 国产精品国产三级专区第一集| 久久热在线av| 精品亚洲乱码少妇综合久久| 满18在线观看网站| 青春草亚洲视频在线观看| 亚洲精品第二区| 亚洲精华国产精华液的使用体验| 亚洲天堂av无毛| 亚洲 欧美一区二区三区| 久久久久久久大尺度免费视频| 亚洲,欧美,日韩| 男人操女人黄网站| 亚洲欧美一区二区三区久久| 午夜福利视频在线观看免费| 亚洲av中文av极速乱| www.熟女人妻精品国产| 晚上一个人看的免费电影| 大香蕉久久网| 可以免费在线观看a视频的电影网站 | 欧美日韩视频高清一区二区三区二| 一边摸一边做爽爽视频免费| 国产成人精品福利久久| 中文字幕色久视频| 欧美成人午夜精品| 美女高潮到喷水免费观看| 精品国产乱码久久久久久小说| 99久久99久久久精品蜜桃| 无遮挡黄片免费观看| 制服丝袜香蕉在线| 在线观看国产h片| 久久精品久久久久久久性| 人妻一区二区av| 国产熟女午夜一区二区三区| 亚洲综合精品二区| 免费少妇av软件| 久久久久久久久久久久大奶| 人人妻人人添人人爽欧美一区卜| 极品少妇高潮喷水抽搐| 一级,二级,三级黄色视频| 亚洲国产精品一区三区| 亚洲伊人色综图| 99精国产麻豆久久婷婷| 中文字幕最新亚洲高清| 国产乱人偷精品视频| 肉色欧美久久久久久久蜜桃| 国产高清国产精品国产三级| 999久久久国产精品视频| 国产av精品麻豆| 高清av免费在线| 少妇 在线观看| 男女床上黄色一级片免费看| 色吧在线观看| 日韩一区二区三区影片| 夫妻性生交免费视频一级片| svipshipincom国产片| 天堂俺去俺来也www色官网| 在线天堂中文资源库| 一级毛片 在线播放| 亚洲精品日本国产第一区| 国产精品香港三级国产av潘金莲 | 夫妻午夜视频| 久久久久久久国产电影| 老鸭窝网址在线观看| 久久人人爽av亚洲精品天堂| 日本91视频免费播放| 亚洲激情五月婷婷啪啪| 日韩欧美精品免费久久| 1024香蕉在线观看| 最黄视频免费看| 亚洲综合精品二区| 人妻一区二区av| 欧美xxⅹ黑人| 亚洲男人天堂网一区| 午夜福利视频精品| 欧美日韩亚洲国产一区二区在线观看 | 亚洲第一av免费看| 精品一品国产午夜福利视频| 十八禁网站网址无遮挡| 久久国产精品大桥未久av| 中文欧美无线码| 我的亚洲天堂| 丝袜在线中文字幕| 国产毛片在线视频| 欧美日韩成人在线一区二区| svipshipincom国产片| 午夜激情av网站| 日韩视频在线欧美| 老鸭窝网址在线观看| 亚洲国产欧美在线一区| 婷婷色综合大香蕉| 久久精品久久久久久噜噜老黄| 国产激情久久老熟女| 一级毛片我不卡| 最近的中文字幕免费完整| 精品一区二区三卡| 亚洲一码二码三码区别大吗| 啦啦啦在线观看免费高清www| 中文乱码字字幕精品一区二区三区| tube8黄色片| 51午夜福利影视在线观看| 十分钟在线观看高清视频www| 婷婷色综合www| 亚洲欧美色中文字幕在线| 欧美日韩一区二区视频在线观看视频在线| 最新在线观看一区二区三区 | 99久久综合免费| 高清在线视频一区二区三区| 亚洲第一区二区三区不卡| 国产一区亚洲一区在线观看| 久久国产亚洲av麻豆专区| 一二三四在线观看免费中文在| 狂野欧美激情性xxxx| 免费不卡黄色视频| 日韩大片免费观看网站| 欧美日韩国产mv在线观看视频| 亚洲图色成人| 久久国产精品大桥未久av| 最近最新中文字幕免费大全7| 老司机深夜福利视频在线观看 | 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 亚洲成色77777| 色婷婷av一区二区三区视频| 老司机影院成人| 中文字幕人妻丝袜一区二区 | 国产日韩欧美亚洲二区| 午夜福利在线免费观看网站| 大片电影免费在线观看免费| 免费黄色在线免费观看| 两个人免费观看高清视频| 国产精品无大码| 亚洲精品久久午夜乱码| 亚洲av在线观看美女高潮| 亚洲av日韩精品久久久久久密 | 视频区图区小说| 国产精品熟女久久久久浪| 亚洲,欧美精品.| 亚洲天堂av无毛| 王馨瑶露胸无遮挡在线观看| 纯流量卡能插随身wifi吗| 精品国产露脸久久av麻豆| 少妇人妻久久综合中文| 老司机深夜福利视频在线观看 | 久久精品国产综合久久久| 自线自在国产av| 久久久久久久久久久免费av| 2018国产大陆天天弄谢| 亚洲成人一二三区av| 看免费成人av毛片| e午夜精品久久久久久久| 亚洲四区av| 男女无遮挡免费网站观看| 日本av免费视频播放| 最近中文字幕高清免费大全6| 久久久久久人人人人人| 国产男女超爽视频在线观看| 亚洲精品乱久久久久久| 亚洲精品国产av成人精品| 97在线人人人人妻| 欧美精品一区二区大全| 亚洲欧洲精品一区二区精品久久久 | 老司机深夜福利视频在线观看 | 欧美97在线视频| 亚洲精品aⅴ在线观看| 多毛熟女@视频| 久久久久久免费高清国产稀缺| 亚洲久久久国产精品| 少妇被粗大猛烈的视频| 丰满少妇做爰视频| 在线观看免费视频网站a站| 一区二区三区精品91| 日本爱情动作片www.在线观看| 乱人伦中国视频| 中文字幕人妻丝袜制服| 搡老乐熟女国产| 999久久久国产精品视频| 亚洲欧美精品自产自拍| 两个人看的免费小视频| 亚洲美女搞黄在线观看| 高清黄色对白视频在线免费看| 99久久99久久久精品蜜桃| 老汉色av国产亚洲站长工具| 午夜免费鲁丝| 午夜福利乱码中文字幕| 狂野欧美激情性xxxx| 亚洲图色成人| 久久久精品国产亚洲av高清涩受| 老汉色∧v一级毛片| 亚洲精品日韩在线中文字幕| 97精品久久久久久久久久精品| 国产免费福利视频在线观看| 最近手机中文字幕大全| 免费高清在线观看日韩| 99久久99久久久精品蜜桃| 国产欧美亚洲国产| 国精品久久久久久国模美| 午夜福利,免费看| www.精华液| 在线观看人妻少妇| 色婷婷av一区二区三区视频| 男人舔女人的私密视频| 亚洲欧洲国产日韩| 女的被弄到高潮叫床怎么办| 久久精品亚洲av国产电影网| 最近2019中文字幕mv第一页| 日韩人妻精品一区2区三区| 国产精品久久久久成人av| videosex国产| 久久韩国三级中文字幕| 19禁男女啪啪无遮挡网站| 91老司机精品| 人妻一区二区av| 中文字幕最新亚洲高清| 亚洲第一青青草原| 91老司机精品| 捣出白浆h1v1| 久久99热这里只频精品6学生| 18禁国产床啪视频网站| 国产精品欧美亚洲77777| 波多野结衣av一区二区av| 99久久精品国产亚洲精品| 久久人人97超碰香蕉20202| 国产熟女午夜一区二区三区| 免费看av在线观看网站| 一区二区三区精品91| 亚洲欧美日韩另类电影网站| 男女边摸边吃奶| 黄色视频在线播放观看不卡| 久久免费观看电影| 国产亚洲av高清不卡| 欧美精品一区二区大全| 一级a爱视频在线免费观看| 欧美在线一区亚洲| 亚洲欧美精品综合一区二区三区| 国产精品 国内视频| 热99国产精品久久久久久7| 99国产综合亚洲精品| 中文字幕高清在线视频| 老司机影院毛片| 亚洲国产欧美网| www.av在线官网国产| 国产成人av激情在线播放| a级片在线免费高清观看视频| 国产精品三级大全| 亚洲视频免费观看视频| 久久久精品94久久精品| 熟女av电影| 亚洲av电影在线观看一区二区三区| 热99久久久久精品小说推荐| 欧美久久黑人一区二区| 久久精品国产亚洲av涩爱| 色视频在线一区二区三区| 宅男免费午夜| 国产亚洲最大av| 久久精品aⅴ一区二区三区四区| 中文精品一卡2卡3卡4更新| 久久久国产欧美日韩av| 大片电影免费在线观看免费| 午夜福利一区二区在线看| 桃花免费在线播放| 亚洲av在线观看美女高潮| 国产麻豆69| 精品亚洲乱码少妇综合久久| 成人国语在线视频| 国产成人系列免费观看| 国产激情久久老熟女| 久热这里只有精品99| 亚洲综合色网址| 国产精品久久久av美女十八| 最近最新中文字幕免费大全7| 在线观看人妻少妇| 一级片'在线观看视频| 日本欧美视频一区| 一边亲一边摸免费视频| 亚洲视频免费观看视频| 午夜免费鲁丝| 女性生殖器流出的白浆| 午夜福利影视在线免费观看| 制服丝袜香蕉在线| 亚洲欧美精品综合一区二区三区| 午夜日本视频在线| 一区二区三区激情视频| 免费在线观看黄色视频的| 日韩一区二区三区影片| 观看av在线不卡| 欧美成人精品欧美一级黄| 哪个播放器可以免费观看大片| 中文字幕亚洲精品专区| 亚洲国产av新网站| 天天影视国产精品| 国产成人啪精品午夜网站| 久久毛片免费看一区二区三区| 激情五月婷婷亚洲| 久热这里只有精品99| 国产高清国产精品国产三级| 国产一区二区三区综合在线观看| 91精品伊人久久大香线蕉| 十八禁网站网址无遮挡| 好男人视频免费观看在线| 中文字幕另类日韩欧美亚洲嫩草| a 毛片基地| 伊人亚洲综合成人网| 99热全是精品| av福利片在线| 91成人精品电影| 成人三级做爰电影| 亚洲精品一二三| 久久热在线av| 久久久久久久久久久免费av| 免费在线观看视频国产中文字幕亚洲 | 久久久久久久国产电影| 亚洲免费av在线视频| 极品少妇高潮喷水抽搐| 国产男女内射视频| 亚洲男人天堂网一区| 亚洲精品av麻豆狂野| 啦啦啦 在线观看视频| 9191精品国产免费久久| 精品免费久久久久久久清纯 | 女人高潮潮喷娇喘18禁视频| 日韩伦理黄色片| 国产无遮挡羞羞视频在线观看| 十八禁高潮呻吟视频| 人人妻人人爽人人添夜夜欢视频| 中文字幕色久视频| 十八禁高潮呻吟视频| 亚洲色图 男人天堂 中文字幕| 一二三四在线观看免费中文在| 性高湖久久久久久久久免费观看| 一本色道久久久久久精品综合| 日本vs欧美在线观看视频| 国产日韩一区二区三区精品不卡| 亚洲欧美色中文字幕在线| 在线观看免费午夜福利视频| 午夜激情av网站| 黄频高清免费视频| 一区二区三区四区激情视频| 在线精品无人区一区二区三| 青春草亚洲视频在线观看| 午夜福利乱码中文字幕| 国产精品免费大片| 亚洲国产看品久久| 国产97色在线日韩免费| 大香蕉久久网| 亚洲欧洲精品一区二区精品久久久 | 成人手机av| 国产精品免费视频内射| 中文字幕人妻丝袜制服| 青春草亚洲视频在线观看| 日日啪夜夜爽| 国产精品偷伦视频观看了| 色婷婷久久久亚洲欧美| 国产男女超爽视频在线观看| 国产成人免费无遮挡视频| 国产又爽黄色视频| 亚洲美女黄色视频免费看| 欧美国产精品va在线观看不卡| 日韩视频在线欧美| 国产精品偷伦视频观看了| 国产亚洲一区二区精品| 国产伦理片在线播放av一区| 高清黄色对白视频在线免费看| 亚洲国产最新在线播放| 国产精品久久久人人做人人爽| 免费观看av网站的网址| 一二三四中文在线观看免费高清| 欧美在线黄色| 国产黄频视频在线观看| 不卡av一区二区三区| 久久鲁丝午夜福利片| 午夜久久久在线观看| 999精品在线视频| 超色免费av| 亚洲欧美一区二区三区国产| 天天添夜夜摸| 亚洲国产最新在线播放| 99精国产麻豆久久婷婷| 9热在线视频观看99| 亚洲色图 男人天堂 中文字幕| 女人高潮潮喷娇喘18禁视频| av一本久久久久| 美女高潮到喷水免费观看| 老司机亚洲免费影院| 91国产中文字幕| 熟妇人妻不卡中文字幕| 啦啦啦在线观看免费高清www| 亚洲国产精品999| 国产成人a∨麻豆精品| 久久精品国产亚洲av高清一级| 精品少妇一区二区三区视频日本电影 | 成人18禁高潮啪啪吃奶动态图| 国产免费现黄频在线看| 曰老女人黄片| 看非洲黑人一级黄片| 巨乳人妻的诱惑在线观看| 久久人人爽av亚洲精品天堂| 欧美日韩亚洲综合一区二区三区_| 国产精品久久久久久精品电影小说| 两性夫妻黄色片| 国产亚洲精品第一综合不卡| 中国国产av一级| 国产高清不卡午夜福利| 亚洲欧美一区二区三区国产| 日韩不卡一区二区三区视频在线| 男人添女人高潮全过程视频| 亚洲色图综合在线观看| 女性被躁到高潮视频| 两性夫妻黄色片| 丝袜喷水一区| 欧美日韩视频高清一区二区三区二| 欧美日韩亚洲综合一区二区三区_| 韩国av在线不卡| 亚洲av成人精品一二三区| 午夜福利乱码中文字幕| 卡戴珊不雅视频在线播放| 最近最新中文字幕大全免费视频 | 精品少妇黑人巨大在线播放| 亚洲综合精品二区| 黄片无遮挡物在线观看| 亚洲欧美色中文字幕在线| 91老司机精品| 国产免费福利视频在线观看| 国产精品久久久久久久久免| 精品久久蜜臀av无| 国产在线视频一区二区| 成年人午夜在线观看视频| 欧美日韩视频精品一区| 国产日韩欧美在线精品| 最近中文字幕高清免费大全6| 久久99一区二区三区| 免费看av在线观看网站| 在线观看三级黄色| 黄色视频不卡| 一边摸一边抽搐一进一出视频| 波多野结衣一区麻豆| 久久国产精品男人的天堂亚洲| 91老司机精品| 亚洲av日韩在线播放| 欧美亚洲日本最大视频资源| av网站在线播放免费| 精品酒店卫生间| a级毛片黄视频| 久久精品久久久久久噜噜老黄| 久久精品aⅴ一区二区三区四区| 美女福利国产在线| 伊人亚洲综合成人网| 嫩草影视91久久| 国产高清国产精品国产三级| 香蕉丝袜av| 蜜桃国产av成人99| 一边亲一边摸免费视频| 色94色欧美一区二区| 亚洲欧美精品综合一区二区三区| 国产成人91sexporn| 国产又色又爽无遮挡免| 国产在视频线精品| 哪个播放器可以免费观看大片| 久久精品熟女亚洲av麻豆精品| 亚洲,一卡二卡三卡| a级毛片在线看网站| 亚洲少妇的诱惑av| 人人澡人人妻人| 国产成人精品久久二区二区91 | 在线观看免费日韩欧美大片| 涩涩av久久男人的天堂| 校园人妻丝袜中文字幕| 男的添女的下面高潮视频| 乱人伦中国视频| 亚洲情色 制服丝袜| 免费观看av网站的网址| 久久久久久久大尺度免费视频| 欧美另类一区| 一级毛片黄色毛片免费观看视频| 悠悠久久av| 美女扒开内裤让男人捅视频| 欧美最新免费一区二区三区| 久久久精品区二区三区| 一级片'在线观看视频| 色94色欧美一区二区| 天堂俺去俺来也www色官网| 成人三级做爰电影| 国产免费现黄频在线看| 国产亚洲一区二区精品| 国产亚洲欧美精品永久| 国产精品一区二区在线观看99| 一区二区三区激情视频| 欧美日韩亚洲国产一区二区在线观看 | 男人添女人高潮全过程视频| av卡一久久| 久久精品亚洲av国产电影网| 亚洲国产看品久久| 99久久99久久久精品蜜桃| 如日韩欧美国产精品一区二区三区| 国产精品香港三级国产av潘金莲 | 亚洲少妇的诱惑av| 精品少妇久久久久久888优播| a级毛片在线看网站| 一个人免费看片子| 久久97久久精品| 亚洲人成电影观看|