• <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ì)算
    午夜福利影视在线免费观看| 欧美精品啪啪一区二区三区| 天堂动漫精品| 一级a爱视频在线免费观看| 亚洲欧美精品综合久久99| 色精品久久人妻99蜜桃| 亚洲成av片中文字幕在线观看| 级片在线观看| 亚洲片人在线观看| 波多野结衣高清无吗| 婷婷精品国产亚洲av在线| 国产精品综合久久久久久久免费 | avwww免费| 久久精品91蜜桃| 欧美日韩精品网址| 中国美女看黄片| 国产亚洲av嫩草精品影院| 母亲3免费完整高清在线观看| 在线观看午夜福利视频| 禁无遮挡网站| tocl精华| 黄片大片在线免费观看| 免费在线观看黄色视频的| 中文字幕色久视频| 久久人妻av系列| 丝袜美足系列| 亚洲最大成人中文| 悠悠久久av| 人人妻人人澡欧美一区二区 | 亚洲精品久久成人aⅴ小说| 此物有八面人人有两片| www.自偷自拍.com| 一进一出抽搐动态| 一区在线观看完整版| 亚洲aⅴ乱码一区二区在线播放 | 一区在线观看完整版| 91国产中文字幕| 18美女黄网站色大片免费观看| 制服丝袜大香蕉在线| 国产精品久久久人人做人人爽| 在线观看www视频免费| 国产精品98久久久久久宅男小说| 午夜日韩欧美国产| 亚洲国产日韩欧美精品在线观看 | 成人18禁高潮啪啪吃奶动态图| 国产成人系列免费观看| 国产成人欧美| 妹子高潮喷水视频| 欧美亚洲日本最大视频资源| 久久这里只有精品19| 久久性视频一级片| tocl精华| 久久精品aⅴ一区二区三区四区| 少妇被粗大的猛进出69影院| 久久婷婷人人爽人人干人人爱 | 国产精品国产高清国产av| 中出人妻视频一区二区| 夜夜躁狠狠躁天天躁| 亚洲精品国产一区二区精华液| 波多野结衣高清无吗| 乱人伦中国视频| 18禁观看日本| 搡老岳熟女国产| 一区二区日韩欧美中文字幕| 又紧又爽又黄一区二区| 亚洲国产欧美网| 欧美最黄视频在线播放免费| 国产成人精品久久二区二区免费| 女性被躁到高潮视频| 亚洲中文字幕日韩| 自线自在国产av| 国产av在哪里看| 久久亚洲精品不卡| 99在线人妻在线中文字幕| av在线天堂中文字幕| 麻豆av在线久日| bbb黄色大片| 免费久久久久久久精品成人欧美视频| 国产高清videossex| 精品国产乱子伦一区二区三区| 中文字幕av电影在线播放| 大型黄色视频在线免费观看| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 精品卡一卡二卡四卡免费| 此物有八面人人有两片| 高清黄色对白视频在线免费看| 人妻久久中文字幕网| 色婷婷久久久亚洲欧美| 亚洲激情在线av| 国产精品久久久久久人妻精品电影| 成人国语在线视频| 一边摸一边抽搐一进一出视频| 色在线成人网| 欧美在线一区亚洲| 9色porny在线观看| 欧美+亚洲+日韩+国产| 90打野战视频偷拍视频| 亚洲欧美精品综合久久99| 精品福利观看| 国产精品二区激情视频| 18禁黄网站禁片午夜丰满| 人妻丰满熟妇av一区二区三区| 亚洲成人久久性| 国产精品九九99| 免费在线观看黄色视频的| 亚洲国产欧美网| 色在线成人网| 亚洲黑人精品在线| 国产av一区在线观看免费| 国产乱人伦免费视频| 久久精品国产亚洲av香蕉五月| 啪啪无遮挡十八禁网站| 国产精品永久免费网站| 丝袜在线中文字幕| 欧美日韩亚洲综合一区二区三区_| 99久久国产精品久久久| 欧美激情久久久久久爽电影 | 亚洲成av人片免费观看| netflix在线观看网站| 人人妻人人澡欧美一区二区 | 精品少妇一区二区三区视频日本电影| 日韩欧美三级三区| 国产欧美日韩一区二区精品| 美国免费a级毛片| 日本欧美视频一区| 国产高清有码在线观看视频 | а√天堂www在线а√下载| 中文字幕色久视频| 又紧又爽又黄一区二区| 亚洲狠狠婷婷综合久久图片| 午夜成年电影在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | av天堂久久9| 在线观看免费视频网站a站| 久久精品国产亚洲av香蕉五月| 亚洲av美国av| 亚洲一码二码三码区别大吗| 麻豆国产av国片精品| 亚洲av电影不卡..在线观看| 99精品在免费线老司机午夜| 亚洲专区字幕在线| 亚洲伊人色综图| 欧美激情高清一区二区三区| 女性生殖器流出的白浆| 欧美日韩乱码在线| 香蕉国产在线看| av天堂久久9| 777久久人妻少妇嫩草av网站| 男女下面进入的视频免费午夜 | 非洲黑人性xxxx精品又粗又长| 久久久国产欧美日韩av| 在线永久观看黄色视频| 99国产综合亚洲精品| 亚洲五月天丁香| 亚洲aⅴ乱码一区二区在线播放 | 一夜夜www| 亚洲人成电影观看| 精品熟女少妇八av免费久了| 婷婷丁香在线五月| 成人18禁在线播放| 亚洲av成人一区二区三| 少妇 在线观看| 久久人人精品亚洲av| 少妇被粗大的猛进出69影院| 精品电影一区二区在线| 天天一区二区日本电影三级 | 麻豆av在线久日| 欧美成人免费av一区二区三区| 午夜日韩欧美国产| 国产精品99久久99久久久不卡| 午夜精品国产一区二区电影| 久久婷婷成人综合色麻豆| 黄色毛片三级朝国网站| 悠悠久久av| 日韩 欧美 亚洲 中文字幕| 欧美激情高清一区二区三区| 国产精品秋霞免费鲁丝片| 美女扒开内裤让男人捅视频| 久久人人爽av亚洲精品天堂| 国产精品免费视频内射| 欧美av亚洲av综合av国产av| 亚洲九九香蕉| 99久久精品国产亚洲精品| 久久精品国产综合久久久| 久久婷婷成人综合色麻豆| 欧美黑人欧美精品刺激| 视频在线观看一区二区三区| 亚洲avbb在线观看| 757午夜福利合集在线观看| 欧美日本亚洲视频在线播放| 天天添夜夜摸| 国产色视频综合| 法律面前人人平等表现在哪些方面| 九色亚洲精品在线播放| 欧美黑人精品巨大| 一级片免费观看大全| 最新美女视频免费是黄的| 我的亚洲天堂| 亚洲自拍偷在线| 桃红色精品国产亚洲av| 精品久久久久久久人妻蜜臀av | 变态另类成人亚洲欧美熟女 | 亚洲精品国产色婷婷电影| 午夜日韩欧美国产| 亚洲欧美精品综合一区二区三区| 在线永久观看黄色视频| 大香蕉久久成人网| www.熟女人妻精品国产| 久久国产精品影院| 男人舔女人的私密视频| av在线天堂中文字幕| 一进一出抽搐动态| 成人永久免费在线观看视频| 欧美一区二区精品小视频在线| 亚洲熟妇中文字幕五十中出| 两性夫妻黄色片| 中文字幕久久专区| 777久久人妻少妇嫩草av网站| 亚洲三区欧美一区| 少妇粗大呻吟视频| 久久精品国产亚洲av香蕉五月| 视频在线观看一区二区三区| 国产高清videossex| 999久久久精品免费观看国产| 国产精品二区激情视频| 日本一区二区免费在线视频| 色播在线永久视频| 人人妻,人人澡人人爽秒播| 亚洲精品国产精品久久久不卡| 国产av又大| 麻豆一二三区av精品| 伊人久久大香线蕉亚洲五| 国产精品av久久久久免费| 91麻豆av在线| 日本黄色视频三级网站网址| 丝袜美腿诱惑在线| 国产精品 欧美亚洲| 国产一级毛片七仙女欲春2 | 亚洲精品国产精品久久久不卡| 久久香蕉激情| 国产成+人综合+亚洲专区| 两个人免费观看高清视频| 国产亚洲av嫩草精品影院| 黄色毛片三级朝国网站| 午夜影院日韩av| 十八禁人妻一区二区| 99精品在免费线老司机午夜| 午夜影院日韩av| 亚洲五月天丁香| 久久久久久亚洲精品国产蜜桃av| 国产乱人伦免费视频| 88av欧美| 两人在一起打扑克的视频| 色老头精品视频在线观看| 午夜福利,免费看| 成人免费观看视频高清| 大陆偷拍与自拍| 久久香蕉激情| 亚洲专区字幕在线| 欧美黄色片欧美黄色片| 国产精品秋霞免费鲁丝片| 国产精品国产高清国产av| 精品乱码久久久久久99久播| 免费在线观看亚洲国产| 久久久久国产精品人妻aⅴ院| 日韩欧美国产一区二区入口| 男女下面插进去视频免费观看| 久久伊人香网站| 亚洲中文字幕日韩| 黄片小视频在线播放| 成人手机av| 国产1区2区3区精品| 午夜福利18| 一级片免费观看大全| 久久久久久久久免费视频了| 亚洲全国av大片| 男女之事视频高清在线观看| 色av中文字幕| 一本大道久久a久久精品| 亚洲精品av麻豆狂野| 亚洲九九香蕉| 久久人人精品亚洲av| 免费高清在线观看日韩| 99香蕉大伊视频| 中文字幕久久专区| 久久久久久人人人人人| 久久久国产精品麻豆| 免费在线观看黄色视频的| 人妻丰满熟妇av一区二区三区| 欧美在线黄色| 色精品久久人妻99蜜桃| 99精品在免费线老司机午夜| 亚洲国产精品sss在线观看| 夜夜爽天天搞| 久久国产精品人妻蜜桃| 亚洲全国av大片| 久久人妻福利社区极品人妻图片| 黄片播放在线免费| 老司机靠b影院| 大陆偷拍与自拍| 国产亚洲欧美精品永久| 国产精品1区2区在线观看.| 巨乳人妻的诱惑在线观看| 国产精品亚洲av一区麻豆| 少妇的丰满在线观看| 国产精品美女特级片免费视频播放器 | 国产私拍福利视频在线观看| 99国产综合亚洲精品| 亚洲第一欧美日韩一区二区三区| 久久久久久免费高清国产稀缺| 亚洲性夜色夜夜综合| 日本在线视频免费播放| 亚洲伊人色综图| 久久久久国内视频| 亚洲一区中文字幕在线| 免费在线观看完整版高清| 亚洲免费av在线视频| 精品一区二区三区视频在线观看免费| 亚洲国产欧美一区二区综合| 欧美绝顶高潮抽搐喷水| 色综合站精品国产| 成人欧美大片| 十八禁网站免费在线| 丰满人妻熟妇乱又伦精品不卡| 一级a爱片免费观看的视频| 操出白浆在线播放| 欧美成人午夜精品| 国产精品久久久久久人妻精品电影| 真人做人爱边吃奶动态| 日本精品一区二区三区蜜桃| 久久午夜综合久久蜜桃| 日本 av在线| 夜夜夜夜夜久久久久| 少妇粗大呻吟视频| 好男人在线观看高清免费视频 | 亚洲精品一区av在线观看| 真人一进一出gif抽搐免费| 国产乱人伦免费视频| 香蕉丝袜av| 可以在线观看毛片的网站| 国产私拍福利视频在线观看| 国产精品一区二区在线不卡| 日韩欧美一区二区三区在线观看| 日韩欧美在线二视频| 久久国产乱子伦精品免费另类| 在线观看免费午夜福利视频| 亚洲国产欧美日韩在线播放| 成年人黄色毛片网站| 黄色成人免费大全| 一个人观看的视频www高清免费观看 | 每晚都被弄得嗷嗷叫到高潮| 两性夫妻黄色片| 久久香蕉精品热| 电影成人av| 天天一区二区日本电影三级 | 国产乱人伦免费视频| 亚洲精品中文字幕一二三四区| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 人妻丰满熟妇av一区二区三区| 一进一出好大好爽视频| 国产av一区在线观看免费| 国产一区在线观看成人免费| 色av中文字幕| 99re在线观看精品视频| 成人三级做爰电影| 久久精品aⅴ一区二区三区四区| 中文字幕色久视频| 国产亚洲精品一区二区www| 久久香蕉国产精品| 成人精品一区二区免费| 久久久久久久午夜电影| 看片在线看免费视频| 日本精品一区二区三区蜜桃| 免费观看精品视频网站| 老汉色av国产亚洲站长工具| 一级毛片精品| 亚洲人成77777在线视频| 国产精品电影一区二区三区| 在线天堂中文资源库| 国产成人av激情在线播放| 久久精品国产清高在天天线| 亚洲天堂国产精品一区在线| 亚洲视频免费观看视频| 最新美女视频免费是黄的| 亚洲激情在线av| 亚洲avbb在线观看| 男女做爰动态图高潮gif福利片 | 亚洲第一av免费看| 精品午夜福利视频在线观看一区| 一级a爱片免费观看的视频| 久久草成人影院| 免费看a级黄色片| 国内精品久久久久久久电影| www.自偷自拍.com| 亚洲人成伊人成综合网2020| 精品人妻1区二区| 国产麻豆69| 久久欧美精品欧美久久欧美| 欧美丝袜亚洲另类 | 久久久国产成人精品二区| 久久久久亚洲av毛片大全| 黄色片一级片一级黄色片| 国产精品综合久久久久久久免费 | 国产精品免费视频内射| 黄片大片在线免费观看| 激情视频va一区二区三区| 国产成人影院久久av| 女同久久另类99精品国产91| 一区二区三区精品91| 亚洲熟女毛片儿| 高清毛片免费观看视频网站| 亚洲精品在线美女| 好男人在线观看高清免费视频 | 国产精品日韩av在线免费观看 | 久久久久久久精品吃奶| 国产乱人伦免费视频| 视频在线观看一区二区三区| 亚洲中文字幕日韩| 精品欧美国产一区二区三| 黄色 视频免费看| 国产极品粉嫩免费观看在线| 国产精品影院久久| 欧美丝袜亚洲另类 | 久久国产精品人妻蜜桃| 男女床上黄色一级片免费看| 看片在线看免费视频| 亚洲狠狠婷婷综合久久图片| 国产伦一二天堂av在线观看| а√天堂www在线а√下载| 老司机深夜福利视频在线观看| 精品免费久久久久久久清纯| av免费在线观看网站| 欧美成人一区二区免费高清观看 | 操出白浆在线播放| 久久欧美精品欧美久久欧美| 亚洲第一欧美日韩一区二区三区| 两个人免费观看高清视频| 伦理电影免费视频| 男女床上黄色一级片免费看| 精品高清国产在线一区| 99国产综合亚洲精品| 亚洲最大成人中文| 最近最新免费中文字幕在线| 中文亚洲av片在线观看爽| 一级,二级,三级黄色视频| 亚洲人成网站在线播放欧美日韩| 后天国语完整版免费观看| 久久久久久人人人人人| 最新在线观看一区二区三区| 中文字幕人成人乱码亚洲影| 正在播放国产对白刺激| 99国产精品免费福利视频| 啦啦啦韩国在线观看视频| 国产精品日韩av在线免费观看 | 又大又爽又粗| 丝袜在线中文字幕| 非洲黑人性xxxx精品又粗又长| 久久国产乱子伦精品免费另类| 精品久久久久久成人av| 欧美老熟妇乱子伦牲交| www.自偷自拍.com| 成人18禁高潮啪啪吃奶动态图| 国产精品,欧美在线| 国内精品久久久久精免费| 成人免费观看视频高清| 欧美精品啪啪一区二区三区| 一本久久中文字幕| 天天躁夜夜躁狠狠躁躁| 久久久国产成人精品二区| 国产av一区二区精品久久| 久久久久久久久免费视频了| 90打野战视频偷拍视频| 久久久久国产精品人妻aⅴ院| 麻豆av在线久日| 身体一侧抽搐| 极品人妻少妇av视频| 久久国产精品男人的天堂亚洲| 日本三级黄在线观看| 亚洲男人天堂网一区| 国产成人免费无遮挡视频| 国产av精品麻豆| 美女国产高潮福利片在线看| 成人手机av| 神马国产精品三级电影在线观看 | av超薄肉色丝袜交足视频| 亚洲电影在线观看av| 一进一出抽搐动态| av在线播放免费不卡| 久久久久久久久免费视频了| 叶爱在线成人免费视频播放| 一边摸一边做爽爽视频免费| 极品人妻少妇av视频| 久久亚洲精品不卡| 中文字幕久久专区| 黄频高清免费视频| 亚洲欧美激情在线| 精品国产一区二区久久| 啦啦啦免费观看视频1| 久久久久久久久免费视频了| 免费看美女性在线毛片视频| 中文亚洲av片在线观看爽| 久久久精品国产亚洲av高清涩受| 日韩精品免费视频一区二区三区| 黄色a级毛片大全视频| 国产精品爽爽va在线观看网站 | 亚洲精品国产区一区二| 欧美中文日本在线观看视频| 91麻豆精品激情在线观看国产| 中文字幕av电影在线播放| 国产高清有码在线观看视频 | 国产主播在线观看一区二区| 人妻丰满熟妇av一区二区三区| 91av网站免费观看| 一边摸一边做爽爽视频免费| 日韩一卡2卡3卡4卡2021年| 日本 av在线| 黄色毛片三级朝国网站| 啪啪无遮挡十八禁网站| 悠悠久久av| av片东京热男人的天堂| 国产不卡一卡二| 黄色a级毛片大全视频| 国产伦人伦偷精品视频| 国产一区二区三区视频了| 久久精品国产99精品国产亚洲性色 | 亚洲欧美一区二区三区黑人| 亚洲国产高清在线一区二区三 | www.精华液| 久久精品影院6| 搞女人的毛片| 好男人在线观看高清免费视频 | 亚洲国产精品成人综合色| 妹子高潮喷水视频| www国产在线视频色| 天堂动漫精品| 国产熟女xx| 少妇裸体淫交视频免费看高清 | 成年女人毛片免费观看观看9| 久久精品国产亚洲av香蕉五月| 91精品国产国语对白视频| 亚洲九九香蕉| 欧美性长视频在线观看| 国产精品 国内视频| 老汉色∧v一级毛片| 欧美丝袜亚洲另类 | 精品国产一区二区久久| 国产精品久久久久久精品电影 | 神马国产精品三级电影在线观看 | 老司机午夜十八禁免费视频| 91麻豆精品激情在线观看国产| 咕卡用的链子| 亚洲五月婷婷丁香| 国产野战对白在线观看| 成人18禁高潮啪啪吃奶动态图| 久久人妻福利社区极品人妻图片| 亚洲欧美日韩另类电影网站| 国产av又大| 免费在线观看影片大全网站| 亚洲国产中文字幕在线视频| 一个人免费在线观看的高清视频| 亚洲欧美精品综合一区二区三区| 极品人妻少妇av视频| 久久中文看片网| 一本综合久久免费| tocl精华| 18禁观看日本| 天堂影院成人在线观看| 极品教师在线免费播放| 亚洲中文字幕一区二区三区有码在线看 | 日本黄色视频三级网站网址| 国产97色在线日韩免费| 在线免费观看的www视频| 日韩国内少妇激情av| 欧美日韩一级在线毛片| 久久精品aⅴ一区二区三区四区| 亚洲国产欧美网| 禁无遮挡网站| 美女大奶头视频| 窝窝影院91人妻| 高清黄色对白视频在线免费看| 亚洲全国av大片| 多毛熟女@视频| av在线天堂中文字幕| 成人特级黄色片久久久久久久| 十分钟在线观看高清视频www| 免费在线观看完整版高清| 亚洲va日本ⅴa欧美va伊人久久| 日韩有码中文字幕| 婷婷六月久久综合丁香| 国产亚洲精品av在线| 嫁个100分男人电影在线观看| 亚洲 欧美 日韩 在线 免费| 少妇的丰满在线观看| 国产精品九九99| 黄色成人免费大全| 成人国语在线视频| 国产黄a三级三级三级人| 亚洲色图av天堂| 国产成人一区二区三区免费视频网站| 丰满的人妻完整版| 国产伦一二天堂av在线观看| 中文字幕人妻熟女乱码| 国产蜜桃级精品一区二区三区| 国产亚洲精品av在线| 一区在线观看完整版| av欧美777| 国产亚洲精品综合一区在线观看 | 51午夜福利影视在线观看| 一进一出抽搐gif免费好疼| 国产高清videossex| 日韩欧美免费精品| 久久婷婷成人综合色麻豆| 亚洲第一电影网av|