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

    多弧槽球面配流副潤(rùn)滑特性分析與多目標(biāo)優(yōu)化

    2023-03-07 07:21:38葉紹干胡佳俊賴偉群卜祥建
    關(guān)鍵詞:配流缸體球面

    葉紹干 胡佳俊 侯 亮 賴偉群 卜祥建

    (廈門(mén)大學(xué)機(jī)電工程系, 廈門(mén) 361021)

    0 引言

    軸向柱塞泵作為重要的液壓泵之一,廣泛應(yīng)用于航空航天、農(nóng)業(yè)機(jī)械、車輛工程等領(lǐng)域[1-4],摩擦副的磨損問(wèn)題是其性能衰減甚至失效的重要原因,其中配流副[5]作為軸向柱塞泵最重要的摩擦副之一,潤(rùn)滑性能對(duì)軸向柱塞泵的工況等級(jí)和壽命有著直接影響。

    針對(duì)配流副潤(rùn)滑性能問(wèn)題,國(guó)內(nèi)外學(xué)者開(kāi)展了大量研究。BERGADA等[6]基于流體微分方程考慮了配流副的力學(xué)特性、泄漏量、摩擦力矩,獲得了不同位置下缸體負(fù)載及其波動(dòng)。王凱麗等[7]分析了深海柱塞泵工作參數(shù)、水深環(huán)境參數(shù)等參數(shù)變化時(shí)油膜壓力分布、溫度分布變化規(guī)律。ZHU等[8]測(cè)試分析了柱塞泵配流盤(pán)和缸體孔接觸間摩擦磨損影響。葉紹干等[9]利用多目標(biāo)遺傳算法對(duì)配流盤(pán)密封環(huán)結(jié)構(gòu)參數(shù)進(jìn)行了優(yōu)化。寧志強(qiáng)等[10]利用粒子群算法對(duì)非對(duì)稱軸向柱塞泵的三角槽寬度參數(shù)進(jìn)行了優(yōu)化。竇振華等[11]基于有限元軟件分析了配流副的摩擦磨損及熱力耦合特性。荊崇波等[12]對(duì)球活塞式液壓泵中錐形配流副潤(rùn)滑特性進(jìn)行了分析。張曉剛等[13]通過(guò)改變柱塞泵缸體結(jié)構(gòu)、柱塞數(shù)等,設(shè)計(jì)了一種雙排油內(nèi)外環(huán)并聯(lián)配流結(jié)構(gòu)的軸向柱塞泵。潘陽(yáng)等[14]為了改善流量和壓力突變?cè)肼暣蟮膯?wèn)題,對(duì)一種雙聯(lián)軸向柱塞泵配流盤(pán)結(jié)構(gòu)進(jìn)行了優(yōu)化。以上研究對(duì)于改善配流副潤(rùn)滑特性提供了一定的思路,同時(shí)對(duì)配流副結(jié)構(gòu)優(yōu)化設(shè)計(jì)也有相關(guān)指導(dǎo)。

    為進(jìn)一步改善配流副潤(rùn)滑性能,國(guó)內(nèi)外學(xué)者提出了多種優(yōu)化方案,并進(jìn)行研究與改進(jìn)。CHACON[15]考慮在配流盤(pán)表面引入規(guī)則形狀波紋,進(jìn)而調(diào)整間隙油膜厚度來(lái)改善配流副表面接觸情況。GEFFROY等[16]通過(guò)在缸體腰形槽設(shè)計(jì)壓力槽的方式提供額外支撐作用,提升了油膜承載能力。趙愿等[17]運(yùn)用類似開(kāi)槽方法來(lái)解決配流副磨損問(wèn)題,分析開(kāi)槽端面的油膜動(dòng)壓支撐作用,并結(jié)合數(shù)學(xué)模型仿真與試驗(yàn)對(duì)比分析。童哲銘等[18]基于參考點(diǎn)的非支配排序遺傳算法對(duì)離心泵葉片參數(shù)進(jìn)行了多目標(biāo)尋優(yōu)。CHEN等[19]建立了多尺寸維度下的配流副織構(gòu)優(yōu)化模型,對(duì)配流副微觀織構(gòu)形狀進(jìn)行了優(yōu)化。JI等[20]通過(guò)多重網(wǎng)格法,分析局部織構(gòu)化參數(shù)對(duì)無(wú)量綱化平均壓力的影響。ZHANG等[21]將激光表面紋理技術(shù)應(yīng)用于配流盤(pán),發(fā)現(xiàn)合理的表面織構(gòu)可以起到減小磨損和降低缸體傾覆角的效果。以上優(yōu)化方法研究表明,合理的配流副結(jié)構(gòu)對(duì)提高配流副潤(rùn)滑特性有著顯著影響,能為球面配流副結(jié)構(gòu)優(yōu)化提供理論指導(dǎo)。

    本文提出一種多弧槽球面配流副結(jié)構(gòu),并對(duì)其結(jié)構(gòu)進(jìn)行優(yōu)化。首先,對(duì)多弧槽球面配流副進(jìn)行理論建模;然后,對(duì)多弧槽球面配流副承載特性進(jìn)行仿真分析;最后,利用多目標(biāo)遺傳算法對(duì)多弧槽球面配流副的結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化。

    1 球面配流副潤(rùn)滑模型

    1.1 多弧槽球面配流盤(pán)結(jié)構(gòu)

    本文提出一種沿內(nèi)、外密封帶各布置若干個(gè)圓弧槽的球面配流盤(pán)結(jié)構(gòu)。如圖1所示,弧槽的分布域?yàn)棣誫s到φge之間的夾角區(qū)域,φgs為分布域起點(diǎn)張角,φge為分布域終點(diǎn)張角;由于弧槽數(shù)目為N,故將分布域沿周向分割成N個(gè)區(qū)域,則每個(gè)弧槽所占區(qū)域的包角為φa,有φa=(φge-φgs)/N,φg為槽區(qū)包角;Rgi、Rgo分別為內(nèi)、外密封帶槽區(qū)節(jié)線半徑;wg為槽寬;由于球面配流盤(pán)內(nèi)、外密封帶為球帶狀,因此槽區(qū)內(nèi)深度不一,如圖1中槽區(qū)截面所示,hgs、hge分別為槽區(qū)內(nèi)外側(cè)槽深,θgs、θge分別為槽區(qū)內(nèi)外側(cè)的仰角。

    圖1 球面配流副弧槽結(jié)構(gòu)示意圖Fig.1 Schematic of arc groove structure of spherical valve-plate pair

    定義槽區(qū)中心槽深為hgc,槽區(qū)內(nèi)任一球心張角θ處的槽深為

    hg=hge+Rssin(θge-θ)tan(θge-θgs)

    (1)

    式中Rs——配流盤(pán)球面半徑

    定義槽區(qū)周向占比為槽區(qū)張角φg對(duì)角度φa的比值,則槽區(qū)張角φg為

    φg=δgφa

    (2)

    (3)

    式中δg——槽區(qū)周向占比

    hr為位于槽區(qū)內(nèi)任一點(diǎn)(θ,φ)不考慮槽深的油膜厚度,而考慮槽深度的油膜厚度為

    (4)

    如圖1中槽區(qū)截面所示,由于配流盤(pán)結(jié)構(gòu)存在夾角θgr,該夾角區(qū)域的槽內(nèi)油膜深度被槽內(nèi)壁面截?cái)?,而由幾何分析可得夾角θgr為

    (5)

    假設(shè)槽深極大值為8 μm,球面半徑Rs為400 mm,因此計(jì)算夾角θgr極大值為2×10-5rad。

    由于θgr極小,因此為簡(jiǎn)化運(yùn)算,忽略該夾角區(qū)域,統(tǒng)一采用式(1)計(jì)算槽區(qū)內(nèi)膜厚。由于密封帶寬度較小,考慮便于加工以及保證弧槽具備足夠的槽寬,將內(nèi)外密封帶槽區(qū)節(jié)線半徑Rgi、Rgo分別設(shè)置在內(nèi)、外密封帶中心線上,分別為38.5、60.5 mm。

    1.2 油膜潤(rùn)滑理論模型

    為了分析球面配流副油膜壓力分布,本研究建立球面間隙下的油膜流體動(dòng)力潤(rùn)滑模型[22]

    (6)

    式中η——油液動(dòng)力粘度ρ——油液密度

    h——油膜厚度p——油膜壓力

    vr0——油膜微元沿經(jīng)線切向速度

    vw0——油液微元沿緯線切向速度

    采用有限容積法對(duì)油膜壓力控制方程進(jìn)行離散化,由交錯(cuò)網(wǎng)格法建立如圖2所示網(wǎng)格劃分模型。其中,白節(jié)點(diǎn)代表壓力分布節(jié)點(diǎn),黑節(jié)點(diǎn)代表速度分布節(jié)點(diǎn),每一個(gè)壓力節(jié)點(diǎn)與其周圍的速度節(jié)點(diǎn)之間的區(qū)域組成控制容積(虛線框部分)[22]。

    圖2 網(wǎng)格劃分及節(jié)點(diǎn)控制容積Fig.2 Discretization of mesh and control volume of mesh point

    對(duì)式(6)兩側(cè)化簡(jiǎn)后可得[22]

    aPpP=aNpN+aSpS+aEpE+aWpW+S

    (7)

    其中

    aP=aN+aS+aE+aW

    (8)

    (9)

    (10)

    (11)

    (12)

    (13)

    運(yùn)用環(huán)形三對(duì)角矩陣算法(Circular tridiagonal matrix algorithm,CTDMA)求解壓力場(chǎng)各節(jié)點(diǎn)壓力。在球坐標(biāo)系下,配流副密封帶上同一仰角所在緯線上各壓力節(jié)點(diǎn)的離散化壓力控制方程可由相鄰節(jié)點(diǎn)壓力之間關(guān)系式表述[22],即

    Ajpj=Bjpj+1+Cjpj-1+Dj(j=1,2,…,m)

    (14)

    式中pj-1、pj、pj+1——同一緯線上相鄰節(jié)點(diǎn)上的壓力

    Dj則包含其余壓力節(jié)點(diǎn)的作用。

    式(14)可進(jìn)一步簡(jiǎn)化為相鄰兩節(jié)點(diǎn)及一個(gè)固定節(jié)點(diǎn)之間的關(guān)系式[22],即

    pj=Ejpj+1+Fjpm+Gj

    (15)

    當(dāng)j=1,有

    (16)

    當(dāng)1

    Aj+1pj+1=Bj+1pj+2+Cj+1pj+Dj+1=

    Bj+1pj+2+Cj+1(Ejpj+1+Fjpm+Gj)+Dj+1

    (17)

    其中

    (18)

    當(dāng)j=m,有

    (Am-BmF1)pm=BmE1p2+BmG1+Cmpm-1+Dm

    (19)

    對(duì)式(19)做系數(shù)替換有

    J2pm=K2(E2p3+F2pm+G2)+Cmpm-1+L2

    (20)

    化簡(jiǎn)可得

    (J2-K2F2)pm=K2E2p3+K2G2+L2+Cmpm-1

    (21)

    類比式(21)有

    J3pm=K3p3+L3+Cmpm-1

    (22)

    進(jìn)而推導(dǎo)可得

    (23)

    由此類推,進(jìn)行反復(fù)運(yùn)算后可得

    (24)

    結(jié)合式(14)~(24),可解得同一緯線上各壓力節(jié)點(diǎn)壓力,進(jìn)而再由密封帶內(nèi)外邊界向內(nèi)部推進(jìn)求解,獲得配流副油膜壓力分布。將前后兩次迭代所得壓力的相對(duì)誤差作為收斂指標(biāo),即

    (25)

    式中k——迭代次數(shù)ε——收斂判據(jù)

    2 仿真

    2.1 多弧槽球面配流副厚度分布與壓力分布

    弧槽結(jié)構(gòu)主要影響配流副油膜流體動(dòng)壓特性,進(jìn)而改善球面配流副的潤(rùn)滑特性。因此,本文設(shè)計(jì)3種弧槽結(jié)構(gòu)參數(shù),建立球面配流副弧槽結(jié)構(gòu)仿真模型,分析不同球面配流副弧槽結(jié)構(gòu)下的油膜壓力分布規(guī)律。

    結(jié)構(gòu)具體參數(shù)如表1所示。選擇仿真工況為:工作壓力35 MPa,工作轉(zhuǎn)速1 500 r/min,斜盤(pán)傾角15.5°。

    表1 弧槽結(jié)構(gòu)參數(shù)Tab.1 Parameters of arc groove structures

    為分析原始結(jié)構(gòu)與弧槽結(jié)構(gòu)的油膜厚度分布和壓力分布規(guī)律,分別提取不同結(jié)構(gòu)下整體油膜厚度分布與對(duì)應(yīng)弧槽結(jié)構(gòu)下外密封帶厚度分布如圖3所示,分別提取不同結(jié)構(gòu)下整體油膜壓力與對(duì)應(yīng)弧槽結(jié)構(gòu)下外密封帶壓力如圖4所示。

    圖3 球面配流副油膜厚度分布對(duì)比Fig.3 Comparisons of thickness distributions of oil film of spherical valve-plate pair

    如圖3所示,相比原始結(jié)構(gòu),弧槽結(jié)構(gòu)在外密封帶最低膜厚分別為3.24、3.13、3.11、3.14 μm,弧槽結(jié)構(gòu)相對(duì)原始結(jié)構(gòu)下降3.1%~4.0%;如圖4所示,對(duì)比不同弧槽結(jié)構(gòu),整體油膜壓力分布最大壓力約為40.5 MPa,原始結(jié)構(gòu)以及弧槽結(jié)構(gòu)1、2、3在外密封帶下的最大壓力分別為35.0、40.7、36.6、35.3 MPa,相比原始結(jié)構(gòu)與弧槽結(jié)構(gòu),最大壓力有顯著上升,增幅最大為16.3%。對(duì)比圖4中外密封帶高壓分布區(qū)域,可見(jiàn)弧槽結(jié)構(gòu)的高壓分布區(qū)域得到延伸。

    圖4 球面配流副油膜壓力分布對(duì)比Fig.4 Comparisons of pressure distribution of oil film of spherical valve-plate pair

    綜上所述,弧槽對(duì)外密封帶油膜厚度與壓力分布特性有著顯著影響。與原結(jié)構(gòu)相比,弧槽結(jié)構(gòu)最低膜厚有所減小,但降低幅度減小,而外密封帶壓力增幅較為明顯,表明弧槽結(jié)構(gòu)可有效提升球面配流副油膜承載能力。

    2.2 弧槽結(jié)構(gòu)對(duì)油膜潤(rùn)滑特性的影響

    為了進(jìn)一步說(shuō)明弧槽結(jié)構(gòu)對(duì)球面配流副油膜潤(rùn)滑特性的影響,本研究對(duì)多弧槽球面配流副和原始球面配流副的潤(rùn)滑特性進(jìn)行對(duì)比分析。

    圖5~7給出了原始結(jié)構(gòu)與弧槽結(jié)構(gòu)1、2、3的缸體傾角、泄漏量以及摩擦轉(zhuǎn)矩隨方位角的變化規(guī)律。

    由圖5可知,相比原始結(jié)構(gòu),弧槽結(jié)構(gòu)1、3的缸體傾角有所下降,分別下降5.1%、3.2%,而弧槽結(jié)構(gòu)2相對(duì)原始結(jié)構(gòu)有輕微提高,提高0.6%。由 圖6、7可知,弧槽結(jié)構(gòu)1、2、3均在泄漏量及摩擦轉(zhuǎn)矩上出現(xiàn)明顯下降,弧槽結(jié)構(gòu)球面配流副泄漏量最大下降8.1%,摩擦轉(zhuǎn)矩最大下降5.9%。

    圖5 不同球面配流副結(jié)構(gòu)下缸體傾角變化曲線Fig.5 Variations of cylinder overturning angles under different spherical valve-plate pair structures

    圖6 不同球面配流副結(jié)構(gòu)下泄漏量變化曲線Fig.6 Variations of leakage losses under different spherical valve-plate pair structures

    圖7 不同球面配流副結(jié)構(gòu)下摩擦轉(zhuǎn)矩變化曲線Fig.7 Variations of friction torques under different spherical valve-plate pair structures

    仿真得到的原始結(jié)構(gòu)與弧槽結(jié)構(gòu)平均缸體傾角、平均泄漏量以及平均摩擦轉(zhuǎn)矩如表2所示,其雷達(dá)圖如圖8所示。

    表2 不同球面配流副結(jié)構(gòu)的平均缸體傾角、平均泄漏量和平均摩擦轉(zhuǎn)矩對(duì)比Tab.2 Comparisons of average cylinder overturning angle, average leakage loss and average friction torque of different spherical valve-plate pair structures

    圖8 不同球面配流副結(jié)構(gòu)下的特征雷達(dá)圖Fig.8 Feature radar chart under different spherical valve-plate pair structures

    如圖8所示,相對(duì)原始結(jié)構(gòu),3種弧槽結(jié)構(gòu)的特征參數(shù)均明顯向內(nèi)收縮,表明弧槽結(jié)構(gòu)可有效改善球面配流副潤(rùn)滑性能。由于3組特征參數(shù)雷達(dá)圖存在互相交叉,說(shuō)明可進(jìn)一步優(yōu)化弧槽結(jié)構(gòu)參數(shù),以獲得更優(yōu)特征值的弧槽結(jié)構(gòu)。

    3 多目標(biāo)優(yōu)化

    3.1 優(yōu)化參數(shù)設(shè)置

    本文優(yōu)化變量為弧槽結(jié)構(gòu)參數(shù),主要包括:分布域起點(diǎn)張角φgs、分布域終點(diǎn)張角φge、槽寬wg、周向槽占比δg、槽數(shù)N、中心槽深hgc,則相應(yīng)的優(yōu)化變量為

    X=(φgs,φge,wg,δg,N,hgc)

    (26)

    本文優(yōu)化目的主要為增強(qiáng)球面配流副抗傾覆能力、提高油膜潤(rùn)滑能力,因此選用缸體傾角γ、泄漏量Q、摩擦轉(zhuǎn)矩Mf為優(yōu)化指標(biāo)。

    由于球面配流副運(yùn)轉(zhuǎn)過(guò)程中,缸體傾角γ、泄漏量Q、摩擦轉(zhuǎn)矩Mf始終處于變化狀態(tài),為便于計(jì)算分析,采用單個(gè)周期內(nèi)的平均值作為目標(biāo)函數(shù),即

    (27)

    最終優(yōu)化目標(biāo)函數(shù)為

    (28)

    優(yōu)化變量取值范圍如表3所示。

    表3 優(yōu)化變量取值范圍Tab.3 Value ranges of optimizable variables

    3.2 優(yōu)化模型

    由于多目標(biāo)優(yōu)化不存在唯一最優(yōu)解,因此最終解需要由決策者從求解得到的可行域中篩選。為保證解集個(gè)體的多樣性,方便后續(xù)解集個(gè)體篩選具備

    更多的可操作空間,本研究運(yùn)用遺傳算法進(jìn)行求解,并且采用并行計(jì)算解決計(jì)算效率問(wèn)題。算法優(yōu)化參數(shù)如表4所示,求解流程如圖9所示[22]。

    圖9 多目標(biāo)遺傳算法優(yōu)化流程圖Fig.9 Optimization process of multi-objective using genetic algorithm[22]

    表4 遺傳算法優(yōu)化參數(shù)Tab.4 Optimization parameters of genetic algorithm

    多目標(biāo)遺傳算法求解流程基本為:①隨機(jī)產(chǎn)生初始種群。②計(jì)算種群個(gè)體適應(yīng)度。③根據(jù)適應(yīng)度和選擇概率進(jìn)行個(gè)體選擇、交叉和變異,產(chǎn)生下一代種群新個(gè)體。④新個(gè)體組成新一代種群進(jìn)行計(jì)算。⑤判定是否滿足終止條件,否則將下一代種群代回步驟②。⑥輸出帕累托最優(yōu)解集。⑦決策者根據(jù)偏好或其它準(zhǔn)則選擇合適的個(gè)體。

    3.3 優(yōu)化結(jié)果及選擇

    為分析不同目標(biāo)組合對(duì)優(yōu)化結(jié)果的影響,將缸體傾角、泄漏量以及摩擦轉(zhuǎn)矩3個(gè)目標(biāo)函數(shù)進(jìn)行兩兩組合,求解分別考慮兩個(gè)目標(biāo)函數(shù)時(shí)的目標(biāo)解集,獲得不同目標(biāo)組合下的帕累托解集U1、U2、U3,如 圖10 所示。

    由圖10可知,在初始種群個(gè)體數(shù)量一致的情況下,不同的目標(biāo)函數(shù)組合得到的帕累托解集個(gè)體數(shù)量存在較大差異:其中,解集U1僅得到4個(gè)個(gè)體;而另外兩種組合分別得到14、17個(gè)個(gè)體,表明缸體傾角和泄漏量?jī)蓚€(gè)目標(biāo)函數(shù)之間的沖突較小,組合得到的解集個(gè)體也較少,相對(duì)容易得到最優(yōu)解。而缸體傾角與摩擦轉(zhuǎn)矩,泄漏量與摩擦轉(zhuǎn)矩表現(xiàn)出較強(qiáng)的沖突,在對(duì)一個(gè)目標(biāo)進(jìn)行優(yōu)化的同時(shí),另一個(gè)目標(biāo)表現(xiàn)出明顯的增大趨勢(shì),相應(yīng)的解集個(gè)體也較多。

    圖10 考慮兩個(gè)目標(biāo)函數(shù)的帕累托解集Fig.10 Pareto solution set considering two objective functions

    由圖10a可知,解集U1個(gè)體的泄漏量在1.13~1.15 L/min之間,差異較小,而當(dāng)泄漏量進(jìn)一步下降時(shí),相應(yīng)個(gè)體的缸體傾角出現(xiàn)了較大跨度,表明泄漏量與缸體傾角之間同樣存在沖突。

    由于目標(biāo)函數(shù)之間存在沖突,并且實(shí)際優(yōu)化通常需要考慮3個(gè)及以上的目標(biāo)函數(shù)。因此,本研究同時(shí)考慮3個(gè)目標(biāo)函數(shù),經(jīng)過(guò)迭代求解,解得最終優(yōu)化種群的帕累托解集U如圖11所示。

    圖11 考慮3個(gè)目標(biāo)函數(shù)的帕累托解集Fig.11 Pareto solution set considering three objective functions

    考慮3個(gè)目標(biāo)函數(shù)的帕累托解集包括46個(gè)個(gè)體,需要根據(jù)實(shí)際需要或者其它選擇準(zhǔn)則對(duì)其篩選,選取一個(gè)或幾個(gè)個(gè)體作為優(yōu)化最優(yōu)解。

    為直觀分析優(yōu)化后的多弧槽球面配流副潤(rùn)滑性能的改善程度,將解集U的個(gè)體目標(biāo)值與原始結(jié)構(gòu)進(jìn)行對(duì)比,如圖12所示。

    圖12 解集U的個(gè)體目標(biāo)函數(shù)Fig.12 Individual objective functions of solution set U

    檢查個(gè)體的所有目標(biāo)函數(shù)值均不高于原始結(jié)構(gòu),并在此基礎(chǔ)上根據(jù)解集U中的所有目標(biāo)函數(shù)值進(jìn)行個(gè)體過(guò)濾,剔除存在一個(gè)或多個(gè)目標(biāo)函數(shù)值大于原始結(jié)構(gòu)的情況的個(gè)體,濾除標(biāo)準(zhǔn)為

    (29)

    式中γ0——原始球面配流副缸體傾角

    Q0——原始球面配流副泄漏量

    Mf0——原始球面配流副摩擦轉(zhuǎn)矩

    由此對(duì)解集U篩選獲得解集L,解集L個(gè)體目標(biāo)值如圖13所示。

    圖13 解集L個(gè)體目標(biāo)函數(shù)Fig.13 Individual objective functions of solution set L

    篩選后解集L中個(gè)體總數(shù)為19個(gè),解集L中的每個(gè)個(gè)體的所有目標(biāo)函數(shù)均小于原始結(jié)構(gòu),有效避免了后續(xù)選擇的解集個(gè)體在個(gè)別目標(biāo)性能上出現(xiàn)表現(xiàn)極端不佳的現(xiàn)象。

    采用線性加權(quán)篩選方法對(duì)個(gè)體目標(biāo)函數(shù)的遞減比例進(jìn)行加權(quán)獲得綜合目標(biāo)函數(shù),從而進(jìn)行二次篩選。設(shè)置3個(gè)目標(biāo)函數(shù)的權(quán)重相等均為1/3,最優(yōu)解被選擇為最佳油膜特性的結(jié)構(gòu)參數(shù),綜合目標(biāo)函數(shù)為

    (30)

    計(jì)算可得每個(gè)個(gè)體的綜合目標(biāo)函數(shù)值如圖14所示。

    圖14 綜合目標(biāo)函數(shù)值Fig.14 Comprehensive objective function values

    經(jīng)過(guò)初次篩選后,解集L保證了個(gè)體在所有目標(biāo)函數(shù)上均優(yōu)于原始結(jié)構(gòu),由圖14可知,解集L中每個(gè)個(gè)體至少在1個(gè)目標(biāo)函數(shù)上相對(duì)原始結(jié)構(gòu)取得了優(yōu)化,其中綜合性能最大優(yōu)化個(gè)體為個(gè)體15,綜合性能提升達(dá)到10.5%,最小優(yōu)化個(gè)體為個(gè)體11,綜合性能提升4.0%。

    選擇優(yōu)化性能最高且相近的前4個(gè)個(gè)體,分別為個(gè)體4、8、14及15,其目標(biāo)函數(shù)值對(duì)比雷達(dá)圖如圖15所示。由圖15可知,尋優(yōu)過(guò)程中出現(xiàn)兩個(gè)趨勢(shì):一個(gè)是沿著缸體傾覆角和泄漏量越來(lái)越小的趨勢(shì),如個(gè)體4、14及15;另一個(gè)是沿著摩擦轉(zhuǎn)矩越來(lái)越小的趨勢(shì),如個(gè)體8。這是由于目標(biāo)函數(shù)之間相互制衡,優(yōu)化算法能有效避免對(duì)單個(gè)目標(biāo)函數(shù)進(jìn)行尋優(yōu),同時(shí)缸體傾角與泄漏量之間為弱沖突;相比之下,缸體傾角、泄漏量均與摩擦轉(zhuǎn)矩之間呈現(xiàn)負(fù)相關(guān)現(xiàn)象,形成較強(qiáng)的矛盾關(guān)系。

    圖15 目標(biāo)函數(shù)值雷達(dá)圖Fig.15 Radar chart of objective function values

    因此,決策者可根據(jù)以下兩種需求進(jìn)行選擇或重新分配權(quán)重尋優(yōu):①對(duì)缸體位姿穩(wěn)定性或容積效率有較高要求,可對(duì)缸體傾角和泄漏量賦予高權(quán)重或設(shè)為主要函數(shù)。②對(duì)球面配流副機(jī)械效率損耗控制要求較高,可對(duì)摩擦轉(zhuǎn)矩賦予高權(quán)重或設(shè)為主要函數(shù)。

    綜上所述,本文提出的多弧槽球面配流副結(jié)構(gòu)可有效提高球面配流副潤(rùn)滑性能,并且可同時(shí)考慮對(duì)多個(gè)性能指標(biāo)優(yōu)化,降低缸體傾角以提高球面配流副抗傾覆能力,降低泄漏量與摩擦轉(zhuǎn)矩以提高潤(rùn)滑性能。

    4 結(jié)論

    (1)提出了一種多弧槽球面配流副結(jié)構(gòu),進(jìn)行了弧槽結(jié)構(gòu)參數(shù)設(shè)計(jì),將弧槽分別置于球面配流盤(pán)內(nèi)、外密封帶,對(duì)多弧槽球面配流副進(jìn)行了理論建模。

    (2)對(duì)弧槽結(jié)構(gòu)承載特性進(jìn)行仿真分析,相比原始結(jié)構(gòu),弧槽結(jié)構(gòu)最小膜厚下降幅度為3.1%~4.0%,弧槽結(jié)構(gòu)最大壓力顯著提高,增幅提高為16.3%。

    (3)基于遺傳算法對(duì)多弧槽球面配流盤(pán)結(jié)構(gòu)進(jìn)行多目標(biāo)參數(shù)尋優(yōu)。對(duì)比原始結(jié)構(gòu),弧槽結(jié)構(gòu)球面配流副缸體傾角最大下降5.1%,泄漏量最大下降8.1%,摩擦轉(zhuǎn)矩最大下降5.9%,最終獲得的最優(yōu)解集個(gè)體綜合目標(biāo)函數(shù)提升10.5%,有效提升了球面配流副的潤(rùn)滑性能。

    猜你喜歡
    配流缸體球面
    重載柱塞泵球面配流副承載特性研究*
    如何消除灰鑄鐵缸體麻點(diǎn)缺陷
    缸體頂面輪廓度加工能力提升分析
    鑄鐵缸體新型電機(jī)孔粗鏜刀具的應(yīng)用
    微觀織構(gòu)配流副熱-流-固耦合潤(rùn)滑特性
    球面檢測(cè)量具的開(kāi)發(fā)
    缸體鑄件清理打磨自動(dòng)化輸送線設(shè)計(jì)
    Heisenberg群上移動(dòng)球面法的應(yīng)用——一類半線性方程的Liouville型定理
    基于博弈論的集裝箱港口聯(lián)盟模型構(gòu)建與應(yīng)用
    基于博弈論的集裝箱港口聯(lián)盟模型構(gòu)建與應(yīng)用
    成人毛片a级毛片在线播放| 亚洲熟妇熟女久久| 舔av片在线| 欧美激情在线99| a在线观看视频网站| 国产真实乱freesex| 女人被狂操c到高潮| 中文字幕人妻熟人妻熟丝袜美| 免费人成在线观看视频色| 日日摸夜夜添夜夜添av毛片 | 久久久成人免费电影| 网址你懂的国产日韩在线| 久久久久久九九精品二区国产| 日本-黄色视频高清免费观看| 国产白丝娇喘喷水9色精品| ponron亚洲| 国产免费一级a男人的天堂| netflix在线观看网站| 女人被狂操c到高潮| 一级毛片久久久久久久久女| 亚洲精品日韩av片在线观看| 伦精品一区二区三区| 精品久久久久久久久av| 欧美日韩精品成人综合77777| 国产女主播在线喷水免费视频网站 | 亚洲四区av| 制服丝袜大香蕉在线| 黄色日韩在线| 国产精品久久电影中文字幕| 禁无遮挡网站| 午夜福利高清视频| 成年版毛片免费区| 久久精品国产亚洲av涩爱 | 日韩精品有码人妻一区| 国产真实伦视频高清在线观看 | 99精品久久久久人妻精品| 嫁个100分男人电影在线观看| 精品欧美国产一区二区三| 国产精品一区www在线观看 | 国产三级中文精品| 国产 一区精品| 热99在线观看视频| 我的女老师完整版在线观看| 偷拍熟女少妇极品色| 亚洲内射少妇av| 国产成年人精品一区二区| 他把我摸到了高潮在线观看| 色吧在线观看| 亚洲av免费在线观看| 亚州av有码| 美女高潮喷水抽搐中文字幕| 亚洲精品国产成人久久av| 国产高清有码在线观看视频| 春色校园在线视频观看| 国产成人aa在线观看| 国产单亲对白刺激| 国产大屁股一区二区在线视频| 99热精品在线国产| 亚洲午夜理论影院| 真人做人爱边吃奶动态| 联通29元200g的流量卡| 精品午夜福利视频在线观看一区| 性欧美人与动物交配| 久久久国产成人精品二区| 欧美成人免费av一区二区三区| 中出人妻视频一区二区| 欧洲精品卡2卡3卡4卡5卡区| 国产精品永久免费网站| 最近最新中文字幕大全电影3| 淫妇啪啪啪对白视频| 美女被艹到高潮喷水动态| 国产在线男女| 国产白丝娇喘喷水9色精品| 一区二区三区高清视频在线| 男人舔女人下体高潮全视频| 国产欧美日韩精品亚洲av| 男人的好看免费观看在线视频| 久久草成人影院| 国产亚洲精品综合一区在线观看| 91午夜精品亚洲一区二区三区 | 午夜福利在线观看吧| 成人亚洲精品av一区二区| 亚洲国产日韩欧美精品在线观看| 人妻夜夜爽99麻豆av| 久久99热这里只有精品18| 欧美性猛交黑人性爽| 一本久久中文字幕| 欧美另类亚洲清纯唯美| 看黄色毛片网站| 欧美日韩黄片免| 久久久久久久久久黄片| 春色校园在线视频观看| 亚洲avbb在线观看| 99久久精品一区二区三区| 伦精品一区二区三区| 久久精品国产亚洲av涩爱 | 亚洲性久久影院| 在线观看av片永久免费下载| 色吧在线观看| 观看美女的网站| 亚洲国产色片| 午夜视频国产福利| 免费看a级黄色片| a级毛片免费高清观看在线播放| 国模一区二区三区四区视频| 直男gayav资源| 国内精品一区二区在线观看| 国产一区二区三区av在线 | 国产主播在线观看一区二区| 超碰av人人做人人爽久久| 三级毛片av免费| 可以在线观看毛片的网站| 99精品在免费线老司机午夜| 精品一区二区三区人妻视频| 亚洲精品色激情综合| 久久天躁狠狠躁夜夜2o2o| 精品日产1卡2卡| 中文资源天堂在线| 欧美日本亚洲视频在线播放| 狠狠狠狠99中文字幕| 免费看av在线观看网站| 欧美日本视频| 欧美日本亚洲视频在线播放| 免费大片18禁| 亚洲一区高清亚洲精品| 国产免费男女视频| h日本视频在线播放| 久久人人精品亚洲av| 一个人免费在线观看电影| 成年免费大片在线观看| 欧美精品啪啪一区二区三区| 日韩一本色道免费dvd| 哪里可以看免费的av片| 精品一区二区三区视频在线| 黄色欧美视频在线观看| 国产精品无大码| 国产精品女同一区二区软件 | 神马国产精品三级电影在线观看| 一区福利在线观看| 亚洲最大成人手机在线| 九色国产91popny在线| 亚洲成人免费电影在线观看| 日韩中字成人| 真人一进一出gif抽搐免费| 亚洲在线自拍视频| 波多野结衣高清无吗| 欧美日韩乱码在线| 别揉我奶头 嗯啊视频| 亚洲av二区三区四区| 久久这里只有精品中国| 精品欧美国产一区二区三| 精品国内亚洲2022精品成人| 日韩欧美 国产精品| 亚洲中文字幕日韩| 香蕉av资源在线| 国产成人影院久久av| av在线老鸭窝| 免费在线观看影片大全网站| 成年女人永久免费观看视频| 51国产日韩欧美| 中文字幕高清在线视频| 国产午夜精品久久久久久一区二区三区 | 欧美高清性xxxxhd video| 看免费成人av毛片| 99在线人妻在线中文字幕| 中国美女看黄片| 中文字幕人妻熟人妻熟丝袜美| 91狼人影院| 婷婷六月久久综合丁香| 能在线免费观看的黄片| 亚洲午夜理论影院| 国产不卡一卡二| 欧美在线一区亚洲| 精品久久久久久久久久免费视频| 三级国产精品欧美在线观看| 九色成人免费人妻av| 久99久视频精品免费| 午夜精品一区二区三区免费看| 999久久久精品免费观看国产| 欧美精品国产亚洲| 97超级碰碰碰精品色视频在线观看| 九色国产91popny在线| 亚洲电影在线观看av| 精品人妻视频免费看| 99热网站在线观看| 琪琪午夜伦伦电影理论片6080| 国产伦精品一区二区三区视频9| 欧美+亚洲+日韩+国产| 亚洲av二区三区四区| 亚洲性久久影院| 看黄色毛片网站| 久久久久久久久久成人| 亚洲精华国产精华液的使用体验 | 99久久精品热视频| 99在线人妻在线中文字幕| 久久6这里有精品| 美女高潮的动态| 啦啦啦观看免费观看视频高清| 日韩国内少妇激情av| 午夜福利在线观看免费完整高清在 | 国产精品免费一区二区三区在线| 校园人妻丝袜中文字幕| 两人在一起打扑克的视频| 久久久久久久亚洲中文字幕| 免费人成在线观看视频色| 精品久久久久久久久av| 最近中文字幕高清免费大全6 | 午夜亚洲福利在线播放| 国产精华一区二区三区| 国产日本99.免费观看| 99久久成人亚洲精品观看| 久久久久久久久中文| 在线天堂最新版资源| 97人妻精品一区二区三区麻豆| 精品午夜福利在线看| 久久久久久久精品吃奶| 国产日本99.免费观看| 99久久精品热视频| 欧美性感艳星| 亚洲图色成人| 国产女主播在线喷水免费视频网站 | 亚洲成人久久性| 国产精品久久久久久精品电影| 97热精品久久久久久| 精品一区二区三区av网在线观看| 大型黄色视频在线免费观看| 嫩草影院新地址| 午夜爱爱视频在线播放| 天堂影院成人在线观看| 免费看a级黄色片| 91久久精品国产一区二区成人| 99热这里只有精品一区| 免费在线观看影片大全网站| 我要看日韩黄色一级片| 国产色爽女视频免费观看| 久久欧美精品欧美久久欧美| 欧美激情在线99| 国产aⅴ精品一区二区三区波| 国产精品一区二区性色av| 综合色av麻豆| 精品久久久久久成人av| 22中文网久久字幕| 91久久精品国产一区二区成人| 国产av麻豆久久久久久久| 午夜影院日韩av| 国产精品久久久久久av不卡| 成人综合一区亚洲| 一级a爱片免费观看的视频| 毛片一级片免费看久久久久 | 免费观看在线日韩| 超碰av人人做人人爽久久| av在线老鸭窝| 久久久午夜欧美精品| 久久6这里有精品| 久久久久久久亚洲中文字幕| 自拍偷自拍亚洲精品老妇| 免费人成视频x8x8入口观看| 亚洲久久久久久中文字幕| 国内精品美女久久久久久| 国产精品爽爽va在线观看网站| 成年免费大片在线观看| 亚洲人成网站在线播放欧美日韩| 成年女人永久免费观看视频| 成熟少妇高潮喷水视频| 日韩欧美 国产精品| 少妇裸体淫交视频免费看高清| 一卡2卡三卡四卡精品乱码亚洲| 国产毛片a区久久久久| 99热网站在线观看| 啦啦啦啦在线视频资源| 国产老妇女一区| 精品久久国产蜜桃| 色尼玛亚洲综合影院| 午夜福利高清视频| 欧美zozozo另类| 天堂影院成人在线观看| 97人妻精品一区二区三区麻豆| 免费观看的影片在线观看| 韩国av一区二区三区四区| 亚洲欧美日韩无卡精品| 18禁在线播放成人免费| 色视频www国产| 波多野结衣巨乳人妻| 人妻夜夜爽99麻豆av| 久久中文看片网| 69av精品久久久久久| 啦啦啦啦在线视频资源| 伊人久久精品亚洲午夜| 免费看a级黄色片| 淫秽高清视频在线观看| 国产成人aa在线观看| 日韩欧美国产在线观看| 在线国产一区二区在线| 男女啪啪激烈高潮av片| 久久99热这里只有精品18| 99国产精品一区二区蜜桃av| 91麻豆av在线| 国产精品1区2区在线观看.| 成人三级黄色视频| 欧美zozozo另类| 校园春色视频在线观看| 成年版毛片免费区| 国产免费av片在线观看野外av| 91久久精品国产一区二区成人| 国产色爽女视频免费观看| 色在线成人网| 99精品久久久久人妻精品| 少妇高潮的动态图| 国产一区二区三区在线臀色熟女| 日韩欧美国产在线观看| 久久久久久久久大av| 麻豆av噜噜一区二区三区| ponron亚洲| 我要看日韩黄色一级片| 免费在线观看成人毛片| 日日摸夜夜添夜夜添小说| 国产精品一区二区性色av| 99久久中文字幕三级久久日本| 久久久久九九精品影院| 99久久成人亚洲精品观看| 在线天堂最新版资源| 少妇丰满av| 给我免费播放毛片高清在线观看| 亚洲av成人精品一区久久| 成人国产综合亚洲| 搡老岳熟女国产| 少妇猛男粗大的猛烈进出视频 | 日本黄色片子视频| 久久久久久久久久成人| 国产主播在线观看一区二区| 天天躁日日操中文字幕| 午夜久久久久精精品| 91久久精品电影网| 看片在线看免费视频| 一进一出好大好爽视频| 精品午夜福利视频在线观看一区| 国产乱人视频| 俺也久久电影网| 日韩,欧美,国产一区二区三区 | a在线观看视频网站| 成人国产综合亚洲| a级毛片a级免费在线| 国产久久久一区二区三区| 日本色播在线视频| 国产日本99.免费观看| 日韩人妻高清精品专区| 亚洲欧美日韩高清专用| 小说图片视频综合网站| 亚洲黑人精品在线| 熟妇人妻久久中文字幕3abv| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 国产免费av片在线观看野外av| 亚洲一区高清亚洲精品| 成人毛片a级毛片在线播放| 国产精品嫩草影院av在线观看 | 此物有八面人人有两片| 精品欧美国产一区二区三| av.在线天堂| 国产高清有码在线观看视频| 亚洲成av人片在线播放无| 国产一区二区激情短视频| 国产又黄又爽又无遮挡在线| 午夜福利欧美成人| 日本撒尿小便嘘嘘汇集6| 欧美+亚洲+日韩+国产| 亚洲四区av| 国产精品1区2区在线观看.| 免费搜索国产男女视频| 欧美精品啪啪一区二区三区| 免费观看人在逋| 久久人妻av系列| 国内揄拍国产精品人妻在线| 干丝袜人妻中文字幕| 日本免费一区二区三区高清不卡| 国产真实乱freesex| 女同久久另类99精品国产91| 免费人成视频x8x8入口观看| 亚洲va日本ⅴa欧美va伊人久久| 日本免费a在线| 人人妻人人看人人澡| 国产男人的电影天堂91| 成人美女网站在线观看视频| 国产高清视频在线观看网站| 九色国产91popny在线| 亚洲精品成人久久久久久| 免费看日本二区| 国产男靠女视频免费网站| 国产精品美女特级片免费视频播放器| 嫩草影视91久久| 好男人在线观看高清免费视频| 又紧又爽又黄一区二区| 久久天躁狠狠躁夜夜2o2o| 在线a可以看的网站| а√天堂www在线а√下载| 搡老熟女国产l中国老女人| 亚洲国产欧美人成| 在线观看舔阴道视频| 久久久久久国产a免费观看| 国产色婷婷99| 欧美日韩精品成人综合77777| 两人在一起打扑克的视频| 又粗又爽又猛毛片免费看| 麻豆国产av国片精品| 91久久精品电影网| 人妻少妇偷人精品九色| 国产高清视频在线播放一区| 熟女人妻精品中文字幕| 日韩一本色道免费dvd| 精品无人区乱码1区二区| 欧美日韩瑟瑟在线播放| 久久天躁狠狠躁夜夜2o2o| 白带黄色成豆腐渣| 国产精品久久久久久av不卡| 少妇的逼水好多| 成年女人永久免费观看视频| 免费观看精品视频网站| 国产黄a三级三级三级人| 小蜜桃在线观看免费完整版高清| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 国产高清视频在线观看网站| 欧美一级a爱片免费观看看| 在线观看一区二区三区| 女的被弄到高潮叫床怎么办 | 亚洲国产欧洲综合997久久,| 久久精品国产鲁丝片午夜精品 | 很黄的视频免费| av在线蜜桃| 一本一本综合久久| 日本熟妇午夜| 色哟哟·www| 色av中文字幕| 欧美一区二区精品小视频在线| 亚洲aⅴ乱码一区二区在线播放| 精品午夜福利视频在线观看一区| 动漫黄色视频在线观看| 18禁在线播放成人免费| 九色国产91popny在线| 日本与韩国留学比较| www.www免费av| 国产精品不卡视频一区二区| 一个人看的www免费观看视频| 十八禁网站免费在线| 久久久精品大字幕| 亚洲av.av天堂| 亚洲图色成人| 日韩欧美一区二区三区在线观看| 国产综合懂色| 精品乱码久久久久久99久播| 久久久精品大字幕| 国产主播在线观看一区二区| 狂野欧美白嫩少妇大欣赏| 亚洲av中文av极速乱 | 两人在一起打扑克的视频| av福利片在线观看| 欧美精品国产亚洲| 亚洲avbb在线观看| 久久热精品热| 一个人观看的视频www高清免费观看| 欧美最黄视频在线播放免费| 悠悠久久av| 国产免费一级a男人的天堂| 深夜精品福利| 男人狂女人下面高潮的视频| 国产极品精品免费视频能看的| 午夜亚洲福利在线播放| 亚洲欧美日韩高清在线视频| 午夜a级毛片| 成人性生交大片免费视频hd| 国产精品一区二区三区四区免费观看 | 久99久视频精品免费| 一级黄色大片毛片| 深爱激情五月婷婷| 免费观看精品视频网站| 国产高清激情床上av| 免费在线观看成人毛片| 亚洲一区二区三区色噜噜| 午夜免费成人在线视频| 日韩亚洲欧美综合| 欧美精品啪啪一区二区三区| 成人特级av手机在线观看| 国产精品久久电影中文字幕| 日本黄色视频三级网站网址| 成人av在线播放网站| 欧美+亚洲+日韩+国产| aaaaa片日本免费| 自拍偷自拍亚洲精品老妇| 国产精品伦人一区二区| 搡女人真爽免费视频火全软件 | 国产精品电影一区二区三区| 免费看美女性在线毛片视频| 在线观看免费视频日本深夜| 99热6这里只有精品| 五月伊人婷婷丁香| 最近最新中文字幕大全电影3| 麻豆国产97在线/欧美| 国产伦一二天堂av在线观看| 高清毛片免费观看视频网站| 久久久午夜欧美精品| 免费黄网站久久成人精品| 国产成年人精品一区二区| 尾随美女入室| a在线观看视频网站| 此物有八面人人有两片| 日本三级黄在线观看| 女生性感内裤真人,穿戴方法视频| 日韩国内少妇激情av| 99久久精品国产国产毛片| av在线观看视频网站免费| 亚洲图色成人| 欧美性猛交黑人性爽| 亚洲av五月六月丁香网| 成人性生交大片免费视频hd| 搡老岳熟女国产| 国产精品永久免费网站| 国产欧美日韩精品亚洲av| 少妇的逼水好多| 欧美日韩乱码在线| 在线a可以看的网站| 亚洲最大成人手机在线| 美女大奶头视频| 色尼玛亚洲综合影院| 国产成年人精品一区二区| 亚洲av中文字字幕乱码综合| 天美传媒精品一区二区| 欧美最新免费一区二区三区| 久久久国产成人免费| 精品人妻偷拍中文字幕| 麻豆av噜噜一区二区三区| 国产精品,欧美在线| 99热精品在线国产| 最近视频中文字幕2019在线8| 天天一区二区日本电影三级| 我要看日韩黄色一级片| 国产乱人伦免费视频| а√天堂www在线а√下载| 麻豆成人av在线观看| 亚洲aⅴ乱码一区二区在线播放| 1000部很黄的大片| 在线看三级毛片| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 国产精品三级大全| 国产探花极品一区二区| 在线观看午夜福利视频| 三级国产精品欧美在线观看| 人妻少妇偷人精品九色| 国产 一区精品| 尤物成人国产欧美一区二区三区| 亚洲内射少妇av| 欧美人与善性xxx| 午夜福利在线观看吧| 欧美日韩黄片免| 特大巨黑吊av在线直播| 88av欧美| 又爽又黄a免费视频| 欧美一区二区国产精品久久精品| 国内精品一区二区在线观看| 午夜爱爱视频在线播放| 精华霜和精华液先用哪个| 国产色爽女视频免费观看| 观看美女的网站| 亚洲成人免费电影在线观看| 亚洲人成伊人成综合网2020| 亚洲一级一片aⅴ在线观看| 中文字幕高清在线视频| 搡老妇女老女人老熟妇| 亚洲成人久久性| 免费看日本二区| 亚洲美女视频黄频| 色综合婷婷激情| 少妇猛男粗大的猛烈进出视频 | 婷婷丁香在线五月| 久久精品国产清高在天天线| 国产欧美日韩精品一区二区| 亚洲最大成人av| 欧美最新免费一区二区三区| 能在线免费观看的黄片| 少妇被粗大猛烈的视频| 亚洲欧美日韩高清在线视频| 免费观看人在逋| 精品不卡国产一区二区三区| 精品欧美国产一区二区三| 色在线成人网| 久久人妻av系列| 一级毛片久久久久久久久女| 在线天堂最新版资源| 18禁黄网站禁片午夜丰满| 高清日韩中文字幕在线| 日韩在线高清观看一区二区三区 | 欧美潮喷喷水| 久久久精品欧美日韩精品| 久久精品91蜜桃| 啦啦啦啦在线视频资源| 久久久久久久久大av| 一个人看视频在线观看www免费| 久久国内精品自在自线图片| 成人毛片a级毛片在线播放| 亚洲avbb在线观看| 国产一区二区亚洲精品在线观看| 午夜免费男女啪啪视频观看 | 亚洲av第一区精品v没综合| 国产白丝娇喘喷水9色精品| 99热6这里只有精品| 色综合亚洲欧美另类图片| 日韩一区二区视频免费看| 成人无遮挡网站| 婷婷精品国产亚洲av在线| 蜜桃久久精品国产亚洲av| 丰满乱子伦码专区| 国产精品98久久久久久宅男小说| 在线观看美女被高潮喷水网站| 国产精品一区二区三区四区免费观看 | 国产精品久久久久久久电影|