• <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ì)算
    久久久国产成人精品二区| 人人妻人人澡人人爽人人夜夜 | 白带黄色成豆腐渣| 午夜免费激情av| 国产午夜福利久久久久久| 啦啦啦观看免费观看视频高清| 亚洲av电影在线观看一区二区三区 | 在线播放国产精品三级| 欧美日韩在线观看h| 晚上一个人看的免费电影| 久久精品久久久久久噜噜老黄 | 国产色婷婷99| 午夜日本视频在线| 婷婷色av中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产色爽女视频免费观看| 99久久无色码亚洲精品果冻| 久久久亚洲精品成人影院| 午夜福利成人在线免费观看| 女人十人毛片免费观看3o分钟| 亚洲精品乱码久久久久久按摩| 观看美女的网站| 亚洲一区高清亚洲精品| 色吧在线观看| 直男gayav资源| 噜噜噜噜噜久久久久久91| 又粗又爽又猛毛片免费看| 亚洲一区高清亚洲精品| 国产精品国产高清国产av| 在线免费观看不下载黄p国产| 国产精品麻豆人妻色哟哟久久 | 成年免费大片在线观看| 国产私拍福利视频在线观看| 人人妻人人澡人人爽人人夜夜 | 美女大奶头视频| 久久久久久久久中文| 一二三四中文在线观看免费高清| 性色avwww在线观看| 超碰97精品在线观看| 亚洲av成人精品一区久久| 日本熟妇午夜| 亚洲成人精品中文字幕电影| av在线天堂中文字幕| 好男人在线观看高清免费视频| 亚洲久久久久久中文字幕| 午夜免费男女啪啪视频观看| 亚洲精品aⅴ在线观看| 搡女人真爽免费视频火全软件| 我要看日韩黄色一级片| 国产精品久久久久久精品电影小说 | 亚洲中文字幕一区二区三区有码在线看| 成人三级黄色视频| 免费观看人在逋| 国产淫片久久久久久久久| 免费黄网站久久成人精品| 91精品一卡2卡3卡4卡| 成年女人看的毛片在线观看| 丰满人妻一区二区三区视频av| 午夜亚洲福利在线播放| a级毛色黄片| 99久久精品一区二区三区| 网址你懂的国产日韩在线| 人人妻人人澡欧美一区二区| or卡值多少钱| 欧美xxxx性猛交bbbb| 国产色爽女视频免费观看| 国产精品1区2区在线观看.| 欧美日韩在线观看h| 精品久久国产蜜桃| 国产精品一区www在线观看| 国产在视频线在精品| 国产探花极品一区二区| 欧美一区二区精品小视频在线| 国产成人免费观看mmmm| 丰满少妇做爰视频| 国产乱人偷精品视频| 少妇熟女aⅴ在线视频| 国产精品综合久久久久久久免费| 亚州av有码| av卡一久久| 亚洲精品456在线播放app| 久久久欧美国产精品| 两个人的视频大全免费| 九色成人免费人妻av| 校园人妻丝袜中文字幕| 成人亚洲精品av一区二区| 一二三四中文在线观看免费高清| 好男人视频免费观看在线| 国产精品久久久久久久久免| 日韩国内少妇激情av| 国产乱人视频| 中文字幕熟女人妻在线| 免费在线观看成人毛片| 色5月婷婷丁香| 又粗又爽又猛毛片免费看| 亚洲国产欧洲综合997久久,| 乱系列少妇在线播放| 精品国内亚洲2022精品成人| 日韩人妻高清精品专区| 日韩强制内射视频| 麻豆乱淫一区二区| 久久久久免费精品人妻一区二区| 国产真实伦视频高清在线观看| 一级毛片电影观看 | 成人亚洲欧美一区二区av| 久久精品夜夜夜夜夜久久蜜豆| 国产黄片视频在线免费观看| 日本三级黄在线观看| 一级黄色大片毛片| 日韩三级伦理在线观看| 午夜激情福利司机影院| 少妇熟女aⅴ在线视频| 五月伊人婷婷丁香| 日韩欧美精品免费久久| 99热精品在线国产| 国产精品福利在线免费观看| 国产精品久久久久久av不卡| 十八禁国产超污无遮挡网站| 男女啪啪激烈高潮av片| 综合色丁香网| 国产在线一区二区三区精 | av在线播放精品| 欧美极品一区二区三区四区| 日本免费在线观看一区| 淫秽高清视频在线观看| 中文字幕免费在线视频6| 伦理电影大哥的女人| 亚洲欧美成人综合另类久久久 | 精品人妻偷拍中文字幕| 91在线精品国自产拍蜜月| 欧美高清性xxxxhd video| 亚洲成人av在线免费| 日日撸夜夜添| 国产美女午夜福利| 亚洲最大成人中文| 我的女老师完整版在线观看| 日韩 亚洲 欧美在线| 看片在线看免费视频| 久久久久久九九精品二区国产| 欧美性猛交黑人性爽| 久久人人爽人人片av| 国产亚洲精品av在线| 日本色播在线视频| 亚洲欧美日韩卡通动漫| 又爽又黄a免费视频| 免费观看在线日韩| 国产视频首页在线观看| 18禁动态无遮挡网站| 免费无遮挡裸体视频| 国产一区二区在线观看日韩| 性色avwww在线观看| 免费观看精品视频网站| 91aial.com中文字幕在线观看| 免费观看人在逋| 国产精品蜜桃在线观看| 国产老妇女一区| 丰满人妻一区二区三区视频av| 午夜精品在线福利| 久久亚洲国产成人精品v| 亚洲成人精品中文字幕电影| 亚洲五月天丁香| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| 麻豆国产97在线/欧美| 久久人人爽人人片av| 天天一区二区日本电影三级| 亚洲av成人精品一二三区| 亚洲高清免费不卡视频| 小蜜桃在线观看免费完整版高清| 国内揄拍国产精品人妻在线| 精品欧美国产一区二区三| 国产精品野战在线观看| 欧美bdsm另类| 人妻少妇偷人精品九色| 成人高潮视频无遮挡免费网站| 成人国产麻豆网| 亚洲va在线va天堂va国产| 日韩视频在线欧美| 哪个播放器可以免费观看大片| 卡戴珊不雅视频在线播放| 别揉我奶头 嗯啊视频| 亚洲电影在线观看av| 男女国产视频网站| 亚洲国产精品国产精品| 亚洲精品久久久久久婷婷小说 | 波多野结衣巨乳人妻| 国产真实伦视频高清在线观看| 精品久久久久久久末码| 久久久久久久久中文| 少妇人妻一区二区三区视频| 午夜精品在线福利| 国产精品乱码一区二三区的特点| 又粗又硬又长又爽又黄的视频| 99热全是精品| 熟女电影av网| 一夜夜www| 国产综合懂色| 久久99蜜桃精品久久| 亚洲精品影视一区二区三区av| 国产色婷婷99| 亚洲精品成人久久久久久| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 亚洲av成人精品一二三区| 久久国内精品自在自线图片| 特级一级黄色大片| av又黄又爽大尺度在线免费看 | 国产女主播在线喷水免费视频网站 | 免费看美女性在线毛片视频| 国产成人精品一,二区| 日韩av不卡免费在线播放| av视频在线观看入口| 国产一区有黄有色的免费视频 | 久久精品久久久久久噜噜老黄 | 国产精品福利在线免费观看| 老司机影院成人| 久久久久久久久久成人| 亚洲内射少妇av| 国产精品久久久久久久久免| 中文字幕av成人在线电影| 亚洲aⅴ乱码一区二区在线播放| 一卡2卡三卡四卡精品乱码亚洲| 天堂av国产一区二区熟女人妻| 尾随美女入室| 女人十人毛片免费观看3o分钟| 国产淫片久久久久久久久| 亚洲精品456在线播放app| 欧美不卡视频在线免费观看| 久久99热这里只有精品18| 午夜福利视频1000在线观看| av专区在线播放| 国产精品一区二区性色av| 丰满乱子伦码专区| 51国产日韩欧美| 人人妻人人看人人澡| 日产精品乱码卡一卡2卡三| 免费观看性生交大片5| 成人亚洲精品av一区二区| 国产精品.久久久| 男的添女的下面高潮视频| 久久亚洲精品不卡| 黑人高潮一二区| 国产精品久久久久久av不卡| 欧美高清性xxxxhd video| 国产白丝娇喘喷水9色精品| 国产成人免费观看mmmm| 久久久久九九精品影院| 久久草成人影院| 91狼人影院| 男人舔奶头视频| 国产一区二区在线av高清观看| 欧美成人a在线观看| 国产精品久久久久久精品电影小说 | 久久6这里有精品| 免费看av在线观看网站| 啦啦啦韩国在线观看视频| 伦理电影大哥的女人| 汤姆久久久久久久影院中文字幕 | 亚洲乱码一区二区免费版| 有码 亚洲区| 国产一级毛片在线| 精品人妻偷拍中文字幕| 丰满少妇做爰视频| 黄色一级大片看看| 久久久精品94久久精品| 精品久久久噜噜| 3wmmmm亚洲av在线观看| 少妇熟女欧美另类| 两个人的视频大全免费| 亚洲高清免费不卡视频| 日韩,欧美,国产一区二区三区 | 日韩国内少妇激情av| 少妇的逼好多水| 欧美成人免费av一区二区三区| av在线蜜桃| 免费搜索国产男女视频| 国产精品99久久久久久久久| 天堂√8在线中文| 日本三级黄在线观看| 国产av在哪里看| 亚洲av免费高清在线观看| 久久鲁丝午夜福利片| 一级毛片电影观看 | 亚洲成人中文字幕在线播放| 国产真实乱freesex| 黄色配什么色好看| 国产精品爽爽va在线观看网站| 国产单亲对白刺激| av在线蜜桃| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 婷婷六月久久综合丁香| 韩国高清视频一区二区三区| 麻豆久久精品国产亚洲av| 18禁在线播放成人免费| 日韩一区二区视频免费看| 国内揄拍国产精品人妻在线| 国产乱来视频区| 日本-黄色视频高清免费观看| 免费大片18禁| 欧美97在线视频| 久久久久精品久久久久真实原创| 日本黄色片子视频| 婷婷色麻豆天堂久久 | 91午夜精品亚洲一区二区三区| 精品久久久久久久末码| 亚洲国产精品合色在线| 美女xxoo啪啪120秒动态图| 免费人成在线观看视频色| 日本wwww免费看| 亚洲国产精品成人久久小说| 人人妻人人澡人人爽人人夜夜 | 99久久无色码亚洲精品果冻| 看免费成人av毛片| 免费大片18禁| 人人妻人人看人人澡| 亚洲精品乱久久久久久| 久久久久久久久久久丰满| 2021天堂中文幕一二区在线观| 男人和女人高潮做爰伦理| 亚洲色图av天堂| 国产视频首页在线观看| 秋霞伦理黄片| 又粗又硬又长又爽又黄的视频| 亚洲国产精品久久男人天堂| 村上凉子中文字幕在线| 国产人妻一区二区三区在| 精品不卡国产一区二区三区| 日韩av在线大香蕉| 天堂av国产一区二区熟女人妻| 免费不卡的大黄色大毛片视频在线观看 | 一级av片app| 久久热精品热| 久久久色成人| 亚洲第一区二区三区不卡| 天天躁日日操中文字幕| 91狼人影院| 色哟哟·www| 精品人妻一区二区三区麻豆| 精品久久久久久电影网 | 国产成年人精品一区二区| 国产在线一区二区三区精 | 人妻系列 视频| 两个人的视频大全免费| 国产亚洲精品久久久com| 伦理电影大哥的女人| 日本免费a在线| 欧美日本视频| 中国国产av一级| 一级毛片aaaaaa免费看小| 日韩欧美三级三区| av在线老鸭窝| 色吧在线观看| 日本一本二区三区精品| av又黄又爽大尺度在线免费看 | 男人舔女人下体高潮全视频| 亚洲精品456在线播放app| 日本五十路高清| 国产成人aa在线观看| 大话2 男鬼变身卡| 久久久国产成人免费| 国产极品精品免费视频能看的| 精品久久久久久久末码| 欧美日韩综合久久久久久| 欧美激情在线99| 天天一区二区日本电影三级| 女人十人毛片免费观看3o分钟| 国产伦在线观看视频一区| or卡值多少钱| 日韩欧美在线乱码| 国产一区二区亚洲精品在线观看| 成人午夜精彩视频在线观看| 亚洲精品色激情综合| 高清在线视频一区二区三区 | 精品久久久久久久末码| 国语自产精品视频在线第100页| 免费观看人在逋| 97超碰精品成人国产| 久久亚洲精品不卡| 欧美高清性xxxxhd video| 91av网一区二区| 国产视频内射| 国产av在哪里看| 久久精品夜色国产| 日韩av不卡免费在线播放| 国产老妇伦熟女老妇高清| 亚洲电影在线观看av| 18禁在线无遮挡免费观看视频| 久久精品国产亚洲av涩爱| 国产精品久久久久久精品电影| 免费看av在线观看网站| 黄色一级大片看看| 国产探花极品一区二区| 99久久精品热视频| 国产中年淑女户外野战色| 国产精品1区2区在线观看.| 久久精品综合一区二区三区| 国产91av在线免费观看| 欧美区成人在线视频| 少妇猛男粗大的猛烈进出视频 | 国产免费视频播放在线视频 | 国产伦在线观看视频一区| 床上黄色一级片| 特大巨黑吊av在线直播| 婷婷色麻豆天堂久久 | 黄色日韩在线| 国产av不卡久久| 日本黄色片子视频| 欧美+日韩+精品| 亚洲自偷自拍三级| 女人被狂操c到高潮| 国产精品1区2区在线观看.| 免费看美女性在线毛片视频| 国产日韩欧美在线精品| 亚洲欧美日韩无卡精品| 精品久久久久久久末码| 国产精品人妻久久久久久| 日本wwww免费看| 三级毛片av免费| h日本视频在线播放| 色综合站精品国产| 乱码一卡2卡4卡精品| 欧美日韩在线观看h| 国内精品宾馆在线| 国产白丝娇喘喷水9色精品| 乱人视频在线观看| 大又大粗又爽又黄少妇毛片口| 国产伦精品一区二区三区视频9| 国产成人aa在线观看| 成年女人看的毛片在线观看| 嫩草影院入口| 亚洲性久久影院| 日韩强制内射视频| 国产伦精品一区二区三区四那| 久久久久久伊人网av| 精品久久久久久久久久久久久| 蜜臀久久99精品久久宅男| 成人av在线播放网站| 成人毛片60女人毛片免费| 亚洲欧美清纯卡通| 你懂的网址亚洲精品在线观看 | h日本视频在线播放| 卡戴珊不雅视频在线播放| 岛国毛片在线播放| 久久6这里有精品| 久久久欧美国产精品| 国产成人福利小说| 精品人妻一区二区三区麻豆| 搡女人真爽免费视频火全软件| 国产老妇女一区| 亚洲精品乱久久久久久| 亚洲人成网站在线观看播放| 综合色av麻豆| 欧美不卡视频在线免费观看| 一本一本综合久久| 亚洲精品,欧美精品| 中文字幕亚洲精品专区| 亚洲自拍偷在线| 国产精品久久久久久久久免| 我的女老师完整版在线观看| 一二三四中文在线观看免费高清| 精品欧美国产一区二区三| 99在线人妻在线中文字幕| 国产在线男女| 色哟哟·www| 久久久久久久久久成人| 国模一区二区三区四区视频| 欧美一区二区精品小视频在线| 久久99热6这里只有精品| 久久精品久久久久久噜噜老黄 | 精品人妻视频免费看| 99久久人妻综合| 亚州av有码| 国产欧美日韩精品一区二区| 国产黄片视频在线免费观看| 国产精品一区二区三区四区免费观看| 99久久无色码亚洲精品果冻| 日韩中字成人| 亚洲国产精品国产精品| 国产亚洲午夜精品一区二区久久 | 亚洲精品乱码久久久v下载方式| 99久久中文字幕三级久久日本| 免费一级毛片在线播放高清视频| 噜噜噜噜噜久久久久久91| 亚洲在久久综合| 亚洲婷婷狠狠爱综合网| 22中文网久久字幕| 久久草成人影院| 日韩大片免费观看网站 | 国产精品国产三级专区第一集| 男人狂女人下面高潮的视频| 国产色婷婷99| 日本爱情动作片www.在线观看| 国产成人精品一,二区| 国产亚洲av嫩草精品影院| 精品人妻熟女av久视频| 亚洲欧美清纯卡通| 啦啦啦啦在线视频资源| 亚洲精品国产成人久久av| 秋霞在线观看毛片| 国产欧美日韩精品一区二区| 国产私拍福利视频在线观看| 国产v大片淫在线免费观看| 色噜噜av男人的天堂激情| 非洲黑人性xxxx精品又粗又长| 日韩av在线免费看完整版不卡| 久久人人爽人人片av| 精品熟女少妇av免费看| 国产片特级美女逼逼视频| 久久精品人妻少妇| 人妻少妇偷人精品九色| 草草在线视频免费看| 黄片无遮挡物在线观看| 久久久久久久久久久免费av| 国产三级在线视频| 日韩亚洲欧美综合| 久久精品久久久久久久性| 午夜久久久久精精品| 岛国毛片在线播放| 国产极品天堂在线| 日韩av在线大香蕉| 欧美另类亚洲清纯唯美| 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 久久久精品94久久精品| 国产精品乱码一区二三区的特点| 插阴视频在线观看视频| 人体艺术视频欧美日本| 在线免费观看不下载黄p国产| 婷婷色综合大香蕉| 欧美3d第一页| 国产精品,欧美在线| 中文字幕av在线有码专区| 麻豆国产97在线/欧美| 在线播放无遮挡| 国产精品久久久久久av不卡| av.在线天堂| av卡一久久| 成年版毛片免费区| 国产精品久久久久久av不卡| 又粗又爽又猛毛片免费看| 午夜激情欧美在线| 99热这里只有精品一区| 国产高潮美女av| 一夜夜www| 精品久久久久久成人av| 午夜久久久久精精品| 精品人妻偷拍中文字幕| 好男人在线观看高清免费视频| 国产探花极品一区二区| 成人毛片a级毛片在线播放| 欧美一区二区国产精品久久精品| 一级爰片在线观看| 最新中文字幕久久久久| 91久久精品国产一区二区成人| 国产在线男女| 搡女人真爽免费视频火全软件| 如何舔出高潮| 亚洲精品自拍成人| 国产一区二区亚洲精品在线观看| 中文字幕人妻熟人妻熟丝袜美| 纵有疾风起免费观看全集完整版 | 日韩精品青青久久久久久| 男女下面进入的视频免费午夜| 男插女下体视频免费在线播放| 一卡2卡三卡四卡精品乱码亚洲| 亚洲性久久影院| 精品午夜福利在线看| 中文字幕精品亚洲无线码一区| 欧美精品一区二区大全| 成人午夜高清在线视频| 久久99热6这里只有精品| 国产成人freesex在线| eeuss影院久久| 亚洲av不卡在线观看| 免费播放大片免费观看视频在线观看 | 日本免费a在线| 免费黄色在线免费观看| 欧美激情国产日韩精品一区| 国产成人精品婷婷| 国产大屁股一区二区在线视频| 日韩精品青青久久久久久| 午夜精品在线福利| 国产成人午夜福利电影在线观看| 激情 狠狠 欧美| 欧美日韩国产亚洲二区| 国产亚洲精品久久久com| 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人久久小说| 九九久久精品国产亚洲av麻豆| 人人妻人人看人人澡| 深爱激情五月婷婷| 日韩成人伦理影院| 日本色播在线视频| 国产在视频线在精品| 一级二级三级毛片免费看| 人妻制服诱惑在线中文字幕| 午夜激情欧美在线| 亚洲精品乱久久久久久| 免费看光身美女| 成人性生交大片免费视频hd| 久久这里只有精品中国| 国产成人一区二区在线| 在线免费十八禁| 国内精品美女久久久久久| 久久精品夜色国产| 村上凉子中文字幕在线| 乱码一卡2卡4卡精品| 高清午夜精品一区二区三区| av视频在线观看入口| 国产精品综合久久久久久久免费| 亚洲图色成人| 91av网一区二区| 亚洲丝袜综合中文字幕| 欧美成人一区二区免费高清观看|