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

    柱塞式能量回收馬達(dá)配流副錐度角的分析與試驗(yàn)

    2022-02-06 00:47:44李少年包尚令周致元
    農(nóng)業(yè)工程學(xué)報 2022年20期
    關(guān)鍵詞:配流角為錐度

    李少年,梁 濤,包尚令,王 煜,周致元,李 曦

    柱塞式能量回收馬達(dá)配流副錐度角的分析與試驗(yàn)

    李少年1,梁 濤1,包尚令1,王 煜2,周致元1,李 曦1

    (1. 蘭州理工大學(xué)能源與動力工程學(xué)院,蘭州 730050;2. 燕山大學(xué)機(jī)械工程學(xué)院,秦皇島 066099)

    柱塞式能量回收馬達(dá)是將液壓馬達(dá)與發(fā)電機(jī)一體化的新一代液壓能量回收裝置,缸體-配流軸組成的配流副是其關(guān)鍵摩擦副之一,配流副配合面錐度角的選擇對馬達(dá)的配流、承載和摩擦磨損特性有重要影響。該研究采用理論分析、數(shù)值模擬和試驗(yàn)測試的方法,探討柱塞式能量回收馬達(dá)配流副錐度角的最優(yōu)值選擇。首先根據(jù)配流副結(jié)構(gòu)與尺寸,明確錐度角范圍,然后以36°、39°、42°和45°共4個配流副錐度角為對象。分別從流場仿真、弱流固耦合和摩擦磨損試驗(yàn)3個方面,評價各錐度角配流副的柱塞腔油液壓力與壓力脈動、配流副部件應(yīng)力與變形、配流副摩擦磨損等性能。結(jié)果發(fā)現(xiàn)配流副錐度角為42°和45°時,位于配流副上死點(diǎn)的柱塞腔內(nèi)油液壓力和壓力波動較小,壓力分別為4.66、4.62 MPa、壓力波動幅度分別為3.307和3.246 MPa;在柱塞腔與高壓油孔接通階段,柱塞腔油液壓力波動幅度分別為0.324、0.322 MPa;兩種錐度角下的配流軸最大等效應(yīng)力皆遠(yuǎn)小于其屈服強(qiáng)度;錐度角為42°缸體的最大等效應(yīng)力占屈服強(qiáng)度比例較45°錐度角大0.74個百分點(diǎn),最大變形量大0.251m;兩種錐度角的配流副沒有強(qiáng)度失效的風(fēng)險,雖然有微量彈性變形,但對配流副的正常工作影響極小。相較于45°錐度角,42°錐度角摩擦副的平均摩擦系數(shù)小0.012,且波動小、穩(wěn)定性好;上、下試件的磨損率分別小1.966×10-6和7.601×10-6mm3/(N·mm)。所以42°錐度角有利于能量回收馬達(dá)配流副的穩(wěn)定工作及高效運(yùn)轉(zhuǎn)。研究結(jié)果可為柱塞式能量回收馬達(dá)的設(shè)計提供參考。

    馬達(dá);流場;仿真;應(yīng)力;變形;配流副;錐度角;摩擦磨損

    0 引 言

    隨著能源危機(jī)日趨嚴(yán)重,節(jié)能液壓的關(guān)注度日趨提高[1-2]。液壓節(jié)能主要包括液壓系統(tǒng)的優(yōu)化及匹配、液壓元件設(shè)計以及節(jié)流減壓閥等元件損耗能量的回收再利用[3-4]。其中能量回收再利用是應(yīng)用能量回收裝置從液壓回路或者元件中回收能量來提高能源利用率。能量回收裝置的研制除了要滿足較高的能量回收率外,還應(yīng)滿足液壓系統(tǒng)小型化、集成化以及系列化的設(shè)計要求[5]。

    本課題組研制的柱塞式能量回收馬達(dá)將液壓馬達(dá)與發(fā)電機(jī)一體化,可將液壓泵出口溢流流量及特殊工況下的流體能量進(jìn)行收集,通過發(fā)電裝置將壓力能轉(zhuǎn)化為電能,并儲存再利用。該能量回收馬達(dá)具有結(jié)構(gòu)緊湊、體積小、效率高、無外泄漏、低速穩(wěn)定性好、轉(zhuǎn)動慣量低等優(yōu)點(diǎn)。

    與軸向柱塞泵類似,柱塞式能量回收馬達(dá)中存在3對主要摩擦副,即缸體-配流軸摩擦副、缸體-柱塞摩擦副、滾輪-滑道摩擦副。其中缸體-配流軸摩擦副承擔(dān)馬達(dá)內(nèi)部油液的供給與分配、缸體承載等任務(wù),對馬達(dá)的工作性能影響較大。缸體-配流軸摩擦面采用錐面配合形式,如果錐度角設(shè)計不合適,會引起高低壓油液在密閉容腔內(nèi)交換時發(fā)生倒灌、壓力沖擊、噪聲等現(xiàn)象[6-8]。而且配流副在極端工況下的摩擦失效會嚴(yán)重影響能量回收馬達(dá)的壽命與可靠性,最終降低其能量回收率。為了使配流副具有良好的配流特性和摩擦特性,需要對配流副錐度角開展深入研究。

    目前國內(nèi)外學(xué)者對能量回收馬達(dá)配流副及配合面錐度角的研究取得了很多成果,如王廣達(dá)[9]分析了海水淡化泵-馬達(dá)增壓能量回收集成裝置配流副的泄漏途徑,建立了雙排式配流副的水膜泄漏數(shù)學(xué)模型,計算出配流副泄漏量與馬達(dá)流量脈動特性,得到工況參數(shù)及配流盤結(jié)構(gòu)參數(shù)對摩擦轉(zhuǎn)矩的影響。喬豐立等[10-11]最早提出了錐形配流副結(jié)構(gòu)的軸向柱塞泵,并得到該軸向柱塞泵的工作壓力、徑向間隙和錐形配流副摩擦特性等參數(shù)與錐度角的關(guān)系。江浩[12]分析了錐形柱塞式軸向液壓馬達(dá)的壓力、轉(zhuǎn)速的影響因素,與原非錐形配流結(jié)構(gòu)液壓馬達(dá)相比,其額定壓力提升了11%、最高轉(zhuǎn)速提高了34.3%、軸向尺寸縮小了40%。荊崇波等[13-14]針對應(yīng)用于球塞泵中的錐形配流方式,發(fā)現(xiàn)配流軸相對缸體的偏心率以及配流間隙經(jīng)過波動后能夠在較短時間內(nèi)建立新的平衡狀態(tài)。鄒云飛等[15]應(yīng)用彈性流體動力潤滑理論,建立了錐形配流副彈流潤滑模型,提出了配流間隙和偏心率求解的數(shù)值計算方法,得到配流間隙隨工作壓力增大呈線性減小,而偏心率在工作壓力范圍內(nèi)均略小于零的結(jié)論。韋曉磊等[16]分析了含長、短兩種不同阻尼孔構(gòu)造的錐形配流副功率損失和泄漏量,得到含短阻尼孔構(gòu)造的配流副泄漏量和功率損失小,且速度和壓力梯度也低,配流副穩(wěn)定性較好的結(jié)論。Zhou等[17]采用錐形流場的雷諾方程建立錐形配流副的潤滑模型,得到包括軸偏心和分配間隙高度在內(nèi)的靜態(tài)和動態(tài)潤滑性能,并建立試驗(yàn)臺驗(yàn)證模型,該模型預(yù)測工作壓力越高,分配間隙和偏心率越大,轉(zhuǎn)速越高,分布間隙和偏心率越小。

    綜上可知,已有錐形配流副的研究主要針對柱塞腔內(nèi)的壓力、流量、間隙和偏心率開展,配流副配合面錐度角的選取對馬達(dá)的配流特性、摩擦副的摩擦特性有較大影響,確定合適的錐度角至關(guān)重要。本文基于Fluent流場仿真、弱流固耦合和摩擦磨損試驗(yàn)方法,對柱塞腔壓力、應(yīng)力與變形和摩擦磨損幾個方面進(jìn)行評價,得到配流副配合面最優(yōu)錐度角,以達(dá)到提高配流副的配流特性與摩擦特性,提高柱塞式能量回收馬達(dá)效率,延長工作壽命。

    1 柱塞式能量回收馬達(dá)

    柱塞式能量回收馬達(dá)主要由配流軸、缸體、滾輪柱塞組、主滑道、輔助滑道、傳動環(huán)、缸體卡套、傳動銷塊、定子(線圈)、轉(zhuǎn)子(永磁體、密封帶)、高低壓油接頭等組成,如圖1所示。

    1.高壓油接頭 2.低壓油接頭 3.滾輪柱塞組 4.缸體 5.配流軸 6.轉(zhuǎn)子(永磁體、密封帶) 7.定子(線圈) 8.散熱殼體 9.低壓油口 10.接線盒 11.主滑道 12.輔助滑道 13.中心彈簧 14.缸體卡套 15.傳動銷塊 16.傳動環(huán) 17.端面散熱板

    柱塞式能量回收馬達(dá)采用發(fā)電機(jī)轉(zhuǎn)子內(nèi)嵌軸向柱塞式液壓馬達(dá),缸體直接驅(qū)動永磁同步交流發(fā)電機(jī)。工作時,外部高壓油液從高壓油接頭進(jìn)入配流軸高壓油液入口,并通過配流軸上高壓油孔流入柱塞腔。柱塞腔內(nèi)油液的高壓作用使柱塞滾輪緊貼主滑道,借助主滑道的反作用力使?jié)L輪柱塞組帶動缸體旋轉(zhuǎn)。旋轉(zhuǎn)的缸體通過傳動環(huán)和傳動銷塊驅(qū)動電機(jī)轉(zhuǎn)子做旋轉(zhuǎn)運(yùn)動,鑲嵌在電機(jī)定子上的線圈因永磁體切割磁力線而產(chǎn)生感應(yīng)電流,從而實(shí)現(xiàn)發(fā)電功能。同時油液從配流軸內(nèi)部的低壓出口通過端面散熱板中間油孔和散熱殼體中的流道在低壓油口聚集流回油箱,將電機(jī)鐵心線圈產(chǎn)生的熱量帶走,起到冷卻作用。因此在安全溢流閥回路并聯(lián)該裝置,可將液壓泵溢流能量轉(zhuǎn)換為電能,儲存到超級電容或者蓄電池,實(shí)現(xiàn)能量的二次利用。

    該柱塞式能量回收馬達(dá)的額定壓力為16 MPa、轉(zhuǎn)速300 r/min。為使其輸出較大的扭矩,在縮短柱塞行程的同時優(yōu)化柱塞部分的受力,內(nèi)部采用雙排柱塞軸向?qū)ΨQ結(jié)構(gòu),可大幅降低噪聲與振動,同時提高電動機(jī)轉(zhuǎn)子轉(zhuǎn)動的穩(wěn)定性。配流副采用錐面配合面,即配流軸與缸體配合面為錐形結(jié)構(gòu)。配流軸呈對稱分布,其表面為外圓錐面,配流軸內(nèi)部設(shè)有高壓油液入口、低壓油液出口、高壓油孔及低壓油孔,缸體為內(nèi)圓錐面,缸體與配流軸構(gòu)成配流副,如圖2所示。馬達(dá)以軸配流的方式工作,滾輪在1轉(zhuǎn)有4個周期的內(nèi)曲線滑道上運(yùn)動,借助滑道曲線形式的轉(zhuǎn)換,調(diào)整馬達(dá)排量,可滿足多種工況及使用場合。

    2 配流副錐度角的仿真分析

    參考相關(guān)文獻(xiàn)[11-13],同時考慮柱塞式能量回收馬達(dá)配流副結(jié)構(gòu)特征及尺寸的限制,確定配流副錐度角的取值范圍為36°~45°,本文取36°、39°、42°和45°共4個錐度角進(jìn)行研究,通過Fluent軟件進(jìn)行流場仿真,分析不同錐度角配流副的馬達(dá)柱塞腔內(nèi)部油液壓力波動及配流副的泄漏量變化。

    2.1 配流副流體域抽取與網(wǎng)格劃分

    利用Solidworks軟件的布爾運(yùn)算抽取配流副的流體域。以42°錐度角配流副為例,得到配流副流體域模型,流體域主要包括配流軸流體域、油膜和柱塞腔流體域。

    網(wǎng)格數(shù)量對計算結(jié)果有很大影響,為了減小該值所帶來的誤差,劃分了5種網(wǎng)格數(shù)量模型,對柱塞腔壓力波動值進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證[18],結(jié)果如表1所示。

    表1 配流副網(wǎng)格無關(guān)性驗(yàn)證

    由表1可知,柱塞腔壓力波動幅值隨著網(wǎng)格數(shù)量的增加逐漸縮小,并趨于穩(wěn)定??紤]到計算的準(zhǔn)確性和經(jīng)濟(jì)性,選擇網(wǎng)格數(shù)為1 488 672的第4種網(wǎng)格方案進(jìn)行數(shù)值計算,如圖3所示,其中油膜厚度為20m。油膜網(wǎng)格先畫線網(wǎng)格,再由線網(wǎng)格拉伸構(gòu)成面網(wǎng)格,最后由面網(wǎng)格旋轉(zhuǎn)構(gòu)成體網(wǎng)格[19-20]。

    圖3 錐度角為42°時配流副流體域網(wǎng)格劃分結(jié)果

    2.2 配流副流場流態(tài)確定

    配流副流態(tài)判斷方法是計算出流體流動的雷諾數(shù),將其同臨界雷諾數(shù)作比較,當(dāng)雷諾數(shù)高于臨界雷諾數(shù)時為湍流,否則為層流[21]。雷諾數(shù)計算式如下:

    柱塞式能量回收馬達(dá)工作時,柱塞運(yùn)動速度為

    =(3)

    式中為柱塞最大運(yùn)動行程,m;為缸體的轉(zhuǎn)角,rad;為缸體的旋轉(zhuǎn)角速度,rad/s;為時間,s。

    將配流軸高壓油液入口直徑31 mm和油液流速11 m/s代入式(1),求得配流軸流體域最大雷諾數(shù)R為7 413,大于臨界雷諾數(shù)2 320[21],因此配流軸流體域的流態(tài)為湍流。將柱塞最大運(yùn)動行程60 mm和缸體角速度31.4 rad/s代入式(2)~(3)計算出柱塞速度等同為腔內(nèi)油液流速,由式(1)求得柱塞腔內(nèi)的流體最大雷諾數(shù)R為2 540,大于臨界雷諾數(shù),因此柱塞腔流體域的流態(tài)為湍流。油膜內(nèi)部流動以剪切流為主,默認(rèn)為層流。

    2.3 物理模型與邊界條件

    為了與真實(shí)湍流狀態(tài)保持一致,使計算結(jié)果更準(zhǔn)確,這里選取湍流模型中的Realizable-模型[22-23]。流場仿真采用46號液壓油作為工作介質(zhì),其相關(guān)物性參數(shù)如表 2所示。

    表2 工作介質(zhì)物性參數(shù)

    對求解器進(jìn)行設(shè)定,選取壓力修正算法的Pressure-Based。入口邊界條件設(shè)定壓力入口,根據(jù)能量回收馬達(dá)的額定工作參數(shù),壓力大小為16 MPa,出口邊界為壓力出口,壓力大小為0 MPa。求解控制參數(shù)選用SIMPLEC格式,并對松弛因子進(jìn)行相應(yīng)調(diào)整,便于收斂。由于滾輪運(yùn)動滑道為4周期的內(nèi)曲線形式,故1個周期為90°,為了提高計算效率,本文采用1個周期進(jìn)行分析,設(shè)置求解時間步長0.000 05 s,總步長設(shè)置為1 000。

    2.4 仿真結(jié)果與分析

    2.4.1 配流副流場仿真結(jié)果分析

    這里以錐度角為42°的配流副為例進(jìn)行流場仿真結(jié)果分析,仿真得到的配流副壓力云圖如圖4a所示。為了表述方便,對柱塞標(biāo)上序號1~9。以處在上死點(diǎn)位置(高壓起點(diǎn)位置)1號柱塞為例進(jìn)行分析,其他柱塞做與1號柱塞同樣的周期運(yùn)動。由圖4a可知,3、5、7、9號柱塞為高壓柱塞腔,由于油膜的存在,此時只有1、9號柱塞間壓力等值線相交,其余柱塞壓力等值線均相互獨(dú)立,不相交。等值線相交意味著存在油液從高壓油腔通過油膜向低壓油腔泄漏,會增大馬達(dá)泄漏量。1號柱塞為上死點(diǎn)位置,其左側(cè)為高壓油孔,右側(cè)為低壓油孔。1號柱塞腔內(nèi)油液壓力為4.66 MPa,其值小于進(jìn)油口壓力,這是由于1號柱塞脫離高壓油孔到上死點(diǎn)位置的過程中,柱塞腔閉死容積仍逐漸增大所導(dǎo)致的。當(dāng)缸體繼續(xù)旋轉(zhuǎn)時,柱塞腔閉死容積開始減小,油液在柱塞的作用下從出油口流出,出油口壓力為0.054 MPa。

    圖4b為42°錐度角配流副的速度矢量圖。由圖4b可以看出,油液最大流速發(fā)生在1號柱塞腔內(nèi),為42.9 m/s,由于柱塞腔脫離高壓油孔后,其閉死容積逐漸增大,造成柱塞腔內(nèi)的油液壓力發(fā)生變化,形成較大壓差,使油液流動加劇,故形成最大油液流速。從局部放大圖可以看出,1、9號柱塞油膜內(nèi)、外側(cè)均有油液泄漏,柱塞腔部分油液有進(jìn)有出,說明有油液回流情況,易產(chǎn)生壓力沖擊和噪聲等現(xiàn)象。

    2.4.2 不同錐度角配流副的流場仿真結(jié)果分析

    采用同樣的方法,對36°、39°和45°錐度角的配流副進(jìn)行流場仿真分析,不同錐度角配流副的壓力云圖和速度矢量圖如圖5所示。

    注:1~9為柱塞序號,1號柱塞位于某一滑道的上死點(diǎn)位置。下同。

    圖5 不同錐度角配流副在任一時刻的壓力和速度矢量圖

    由圖5a可看出,不同錐度角下1號柱塞腔內(nèi)的壓力均小于進(jìn)油口壓力。當(dāng)配流副錐度角為36°時,1、9號柱塞間壓力等值線相交,且壓力值均未達(dá)到進(jìn)油口壓力,3、5、7號柱塞相鄰間壓力等值線相交,1號柱塞腔內(nèi)油液壓力為9.25 MPa;當(dāng)配流副錐度角為39°時,除了1、3號柱塞不相交,1、3、5、7、9號柱塞相鄰兩者間均相交,1號柱塞腔內(nèi)油液壓力為4.82 MPa;當(dāng)配流副錐度角為45°時,只有1、9號柱塞間壓力等值線相交,其余柱塞壓力等值線均不相交,1號柱塞腔內(nèi)油液壓力為4.62 MPa。4種錐度角下,1號柱塞腔油液壓力在相鄰兩錐度角間的變化幅度分別為4.43、0.15和0.04 MPa,即當(dāng)配流副錐度角在42°~45°之間時,1號柱塞腔內(nèi)油液壓力變化較小,其余角度時變化較大。

    由圖5b可知油液最大流速均發(fā)生在1號柱塞腔內(nèi),最大達(dá)到52 m/s。隨錐度角的增大,最大流速的變化趨勢同柱塞腔內(nèi)壓力變化趨勢一致。造成該現(xiàn)象的原因是柱塞腔脫離高壓油孔后,其閉死容積逐漸增大,導(dǎo)致柱塞腔內(nèi)的油液壓力發(fā)生變化,形成壓力差,使油液流動加劇。

    2.4.3 不同錐度角配流副柱塞腔內(nèi)油液瞬時壓力變化

    圖6為不同錐度角配流副單個柱塞腔內(nèi)油液瞬時壓力變化曲線圖,選取1號柱塞為例進(jìn)行分析,其他柱塞與1號柱塞類似。

    由圖6可知,在0.01 s左右,1號柱塞腔內(nèi)油液壓力發(fā)生了突變,即產(chǎn)生了壓力波動。這是由于柱塞腔脫離高壓油孔后并未達(dá)到下死點(diǎn)位置(低壓起點(diǎn)位置),隨缸體的轉(zhuǎn)動,柱塞腔仍然被壓縮,因此柱塞腔內(nèi)閉死容積繼續(xù)減小,造成柱塞腔內(nèi)壓力突然增大。當(dāng)柱塞腔到達(dá)下死點(diǎn)位置且缸體持續(xù)轉(zhuǎn)動時,由于柱塞腔還沒有與高壓油孔連通,此時柱塞被拉伸,造成柱塞腔閉死容積增大,使柱塞腔內(nèi)壓力突然減小。直到缸體旋轉(zhuǎn)到柱塞腔與高壓油孔連通時,高壓油液注入柱塞腔,使柱塞腔壓力迅速增大到進(jìn)油口壓力。約0.035 s時,1號柱塞腔內(nèi)油液壓力再次發(fā)生突變。這是由于柱塞腔脫離高壓油孔后還沒有達(dá)到上死點(diǎn)位置,隨缸體的轉(zhuǎn)動,柱塞持續(xù)被拉伸,因此柱塞腔閉死容積繼續(xù)增大,造成柱塞腔壓力突然降低。當(dāng)柱塞腔到達(dá)上死點(diǎn)位置缸體持續(xù)轉(zhuǎn)動時,由于柱塞腔還沒有與低壓油孔連通,此時柱塞被壓縮,因此柱塞腔閉死容積開始減小,造成柱塞腔內(nèi)油液壓力突然增大。直到柱塞腔與低壓油孔連通時,油液流出柱塞腔,油液壓力迅速降低到出油口壓力。

    圖6 不同錐度角配流副單個柱塞腔內(nèi)油液瞬時壓力變化曲線

    由圖6可以看出,在0.01、0.035 s,即下死點(diǎn)和上死點(diǎn)位置附近,當(dāng)配流副錐度角為36°時,壓力波動最大,壓力波動幅度最大為7.077 MPa,壓力變化滯后于其他錐度角,且伴有負(fù)超調(diào)發(fā)生,易形成空化現(xiàn)象。當(dāng)錐度角為39°、42°和45°時,壓力波動幅度分別為4.355、3.307和3.246 MPa。由此可知,隨著配流副錐度角的增大,壓力波動幅度不斷減小。當(dāng)錐度角為42°和45°時,壓力波動幅度接近且數(shù)值較小。

    柱塞腔與高壓油孔接通時,4種錐度角下柱塞腔油液壓力波動幅度分別為0.365、0.417、0.324和0.322 MPa。配流副錐度角42°和45°時,柱塞腔內(nèi)油液壓力波動幅度接近且數(shù)值較小。

    2.4.4 不同錐度角配流副油膜處泄漏流量分析

    由上文分析可知,油膜兩側(cè)明顯存在油液泄漏現(xiàn)象,因此對不同錐度角配流副油膜處泄漏流量進(jìn)行分析,得到油膜處泄漏量變化曲線如圖7所示。

    由圖7可知,不同錐度角配流副所在油膜處的泄漏流量均較小,在(0.5~1.9)×10-3L/min范圍變化,滿足馬達(dá)泄漏量小的要求。若泄漏量很大,則會引起馬達(dá)容積效率下降[24]。配流副錐度角為36°時,泄漏量變化幅度較大,從5×10-4L/min增大到1.8×10-3L/min。而錐度角為39°、42°和45°時,泄漏量最大值分別為1.75×10-3、1.8×10-3和1.9×10-3L/min,所以配流副泄漏量峰值隨著配流軸錐度角的增大而增大,但泄漏量峰值增量均較小,對馬達(dá)容積效率影響較小。

    圖7 不同錐度角配流副的泄漏流量變化曲線

    3 配流副部件的應(yīng)力與形變分析

    由仿真結(jié)果可知,當(dāng)馬達(dá)入口壓力為16 MPa,配流副錐度角為42°和45°時,配流副壓力分布與配流特性較好,而且泄漏量較小,因此這里選擇錐度角為42°與45°的配流副進(jìn)行弱流固耦合仿真,分析錐度角對配流副部件應(yīng)力和變形的影響。

    3.1 數(shù)值計算方案

    3.1.1 弱流固耦合仿真模型

    利用Workbench平臺,將流場仿真計算結(jié)果導(dǎo)入Workbench線性靜力學(xué)模塊,搭建弱流固耦合仿真模型[25-26]。并利用Workbench自帶網(wǎng)格劃分工具對配流副劃分網(wǎng)格,得到配流副固體域網(wǎng)格模型(圖8),最低單元網(wǎng)格質(zhì)量為0.38。

    圖8 配流副固體域網(wǎng)格模型

    3.1.2 邊界條件

    摩擦副通常由一對軟硬材料搭配而成,摩擦副材料要求具有摩擦系數(shù)小、穩(wěn)定性好、耐磨等特點(diǎn),且缸體材料硬度應(yīng)略低于配流軸材料硬度[27-28]。通過調(diào)研和多次試驗(yàn),確定以4Cr13作為缸體材料,以20CrMnTi作為配流軸材料組成的摩擦副摩擦性能最好,材料屬性如表3所示。

    表3 配流副材料性能參數(shù)

    弱流固耦合仿真通過Fluent流場仿真數(shù)據(jù)導(dǎo)出date和case文件,并將其導(dǎo)入Workbench靜力學(xué)模塊,流場邊界條件不變。結(jié)構(gòu)位移邊界條件為:配流軸的進(jìn)、出油口采用固定約束,缸體兩端面采用彈性約束,其他邊界條件均為自由面,同時導(dǎo)入配流副流場仿真的流體壓力,使液體壁面與固體界面形成數(shù)據(jù)交互,進(jìn)行有限元計算。

    3.2 配流副部件的應(yīng)力分布

    通過弱流固耦合仿真得到的配流副部件等效應(yīng)力云圖如圖9所示??梢钥闯?,2種錐度角的配流副中,缸體最大等效應(yīng)力皆位于缸體錐面上的柱塞腔端部。當(dāng)錐度角為42°時,缸體最大等效應(yīng)力為114.63 MPa,占屈服強(qiáng)度(500 MPa)的22.92%;當(dāng)錐度角為45°時,缸體最大應(yīng)力105.89 MPa,占屈服強(qiáng)度的21.18%。兩者相比,45°錐度角時缸體的應(yīng)力狀況略優(yōu)于42°(大0.74個百分點(diǎn)),但兩者最大等效應(yīng)力皆遠(yuǎn)小于屈服強(qiáng)度,因此在配流副正常工作時,2種錐度角缸體應(yīng)力情況良好。配流軸最大等效應(yīng)力皆位于高壓油孔端口處,原因是此處為高壓油液通道的交匯區(qū),受力較大,且槽孔邊緣處易發(fā)生應(yīng)力集中。當(dāng)錐度角為42°時,配流軸的最大等效應(yīng)力62.508 MPa,占屈服強(qiáng)度(835 MPa)的7.49%;當(dāng)錐度角為45°時,配流軸最大等效應(yīng)力62.756 MPa,占屈服強(qiáng)度的7.52%,2種錐度角下配流軸最大等效應(yīng)力皆遠(yuǎn)小于屈服強(qiáng)度,應(yīng)力情況良好。

    圖9 配流副部件等效應(yīng)力

    3.3 配流副部件的變形分析

    圖10為配流副部件變形云圖,可以看出,2種錐度角的配流副中,缸體的最大變形位于缸體緊鄰配流軸的邊緣部位。當(dāng)錐度角為42°時,缸體最大變形量為12.069m;當(dāng)錐度角為45°時,缸體最大變形量11.818m。兩者變形量皆較?。▋H相差0.251m),屬于彈性變形,實(shí)際工作中并無風(fēng)險。同理,配流軸最大變形位于高壓油孔端口位置,由于此處為高壓油液通道的交匯區(qū),且結(jié)構(gòu)尺寸較小,易形成應(yīng)力集中所導(dǎo)致。當(dāng)錐度角為42°時,配流軸最大變形量為4.288 4m;當(dāng)錐度角為45°時,配流軸最大變形量3.471 8m,兩者皆為微小彈性變形,不影響配流副的正常工作。

    圖10 配流副部件變形

    通過弱流固耦合仿真結(jié)果可知,從應(yīng)力和變形角度看,45°錐度角配流副略優(yōu)于42°,兩者受力性能均良好,最大等效應(yīng)力均遠(yuǎn)小于其屈服強(qiáng)度。因此該結(jié)構(gòu)在正常工作時,2種錐度角下的配流副都沒有強(qiáng)度失效的風(fēng)險,雖然均有微量變形,但皆屬于微量彈性變形,對配流副的正常工作幾乎沒有影響。

    4 摩擦磨損試驗(yàn)

    4.1 試件制備

    按照缸體與配流軸間的相互運(yùn)動關(guān)系和摩擦面形狀,將摩擦副按內(nèi)徑等比例縮小1.3倍制作試驗(yàn)零件。上試件使用缸體材料,下試件使用配流軸材料,實(shí)物如圖 11所示。上試件外徑50 mm,內(nèi)徑28 mm,高度12 mm,下試件外徑54 mm,內(nèi)徑38 mm,高度17 mm。試件錐度角取42°和45°,每次試驗(yàn)重復(fù)3次,一共制作3組相同試樣。上、下試樣形位公差和表面質(zhì)量均與實(shí)際產(chǎn)品一致,同時表面需要進(jìn)行淬火、氮化處理,同時進(jìn)行冷處理和時效處理,從而獲得穩(wěn)定且硬度較高的金相組織。

    圖11 試件實(shí)物圖

    4.2 試驗(yàn)方法

    摩擦磨損試驗(yàn)系統(tǒng)主要包括立式萬能摩擦磨損試驗(yàn)機(jī)、小型液壓站以及數(shù)據(jù)采集計算機(jī),如圖12所示。

    圖12 摩擦磨損試驗(yàn)系統(tǒng)

    試驗(yàn)時,上、下試件分別安裝在夾具中,上試件由電機(jī)驅(qū)動旋轉(zhuǎn)。下試件夾具帶有油盒,油盒安裝在試驗(yàn)機(jī)下端的連接軸上,連接軸上安裝有力矩傳感器,用來檢測摩擦試件之間的摩擦力矩。采集到的數(shù)據(jù)及時傳輸?shù)綌?shù)據(jù)采集計算機(jī)并作處理。系統(tǒng)中的小型液壓站利用油液循環(huán)流動將摩擦磨損試驗(yàn)過程中產(chǎn)生的熱帶走,使摩擦副溫度處于21℃左右,消除溫度對試件摩擦特性的影響[29-30]。

    試驗(yàn)前后,將加工好的上試件、下試件放入超聲波清洗機(jī)清洗,用烘干機(jī)將試件表面水分烘干,并靜置2 h。根據(jù)實(shí)際工況與摩擦磨損試驗(yàn)機(jī)特性參數(shù),確定施加載荷為200 N,上試件轉(zhuǎn)速為 300 r/min,試驗(yàn)時間為60 min。在試驗(yàn)前后分別采用精度為0.1 mg 的微精電子天平測量磨損前后的試件質(zhì)量,結(jié)果取多次測量的平均值,得到磨損量。

    4.3 評價指標(biāo)

    4.3.1 摩擦系數(shù)

    摩擦系數(shù)是衡量摩擦性能的重要指標(biāo),在計算不同摩擦副材料對的摩擦系數(shù)時,采用控制變量法,保證材料不變的情況下,僅改變試件的錐度角。摩擦系數(shù)計算式如下[31]:

    式中為摩擦系數(shù);為摩擦力矩,N·m;為摩擦力矩對應(yīng)的摩擦半徑,m;為加載力,N;為摩擦副錐度角,(°)。

    4.3.2 磨損率

    磨損率是反映材料耐磨性能的主要參數(shù),磨損率越小,材料耐磨性越好。磨損率w計算式如下[32]:

    4.4 試驗(yàn)結(jié)果與分析

    4.4.1 不同錐度角摩擦副的摩擦系數(shù)變化

    在3組重復(fù)試驗(yàn)中,選擇摩擦系數(shù)最為穩(wěn)定的一組試驗(yàn)數(shù)據(jù),得到的摩擦系數(shù)變化曲線如圖13所示。

    圖13 摩擦系數(shù)變化曲線

    由圖13可以看出,當(dāng)錐度角為45°時,在整個摩擦試驗(yàn)過程中,摩擦系數(shù)波動劇烈,幅值變化較大。這表明45°錐度角摩擦副的摩擦穩(wěn)定性較差。當(dāng)錐度角為42°時,整體摩擦系數(shù)變化曲線低于45°錐度角,而且波動較小,幅值變化小,只是在試驗(yàn)開始的前7 min有輕微的波動,此時處于磨合階段,隨后趨于穩(wěn)定。分別對2組試件60 min內(nèi)的摩擦系數(shù)取平均值,求得42°與45°錐度角試件的平均摩擦系數(shù)分別為0.155和0.167。42°錐度角摩擦副的摩擦系數(shù)?。ㄏ噍^于45°錐度角,平均摩擦系數(shù)小0.012),穩(wěn)定性好,有利于配流副穩(wěn)定工作,并延長其使用壽命。

    4.4.2 不同錐度角摩擦副的磨損率

    試驗(yàn)得到42°錐度角上、下試件磨損量分別為1.766 7、2.9 mg;45°錐度角上、下試件磨損量分別為2.166 7、4.233 3 mg。通過式(5)計算出42°錐度角上、下試件磨損率分別為1.224 4×10-5、1.995 0×10-5mm3/(N·mm),45°錐度角上、下試件磨損率分別為1.421 0×10-5、2.755 1×10-5mm3/(N·mm)。

    可以看出,2種錐度角時,下試件的磨損率較上試件分別大0.770 6、1.334 1 mm3/(N·mm),說明下試件磨損相對嚴(yán)重一些,這是因?yàn)樯显嚰牟牧嫌捕却笥谙略嚰?2°錐度角上、下試件的磨損率較45°角分別小1.96×10-6、7.601×10-6mm3/(N·mm)。若柱塞式能量回收馬達(dá)配流副應(yīng)用45°錐度角,配流副磨損率大,易引起摩擦力和噪聲增大,以及因內(nèi)部油液泄漏導(dǎo)致馬達(dá)容積效率減小等問題。因此,針對該摩擦副材料,42°錐度角配流副的耐磨性優(yōu)于45°錐度角,工作效率更高和使用壽命更長。

    5 結(jié) 論

    分別使用流場仿真、弱流固耦合仿真和摩擦磨損試驗(yàn)方法,對36°、39°、42°和45°錐度角的能量回收馬達(dá)配流副開展仿真與試驗(yàn)研究,主要結(jié)論如下:

    1)當(dāng)配流副錐度角為42°和45°時,位于配流副上死點(diǎn)柱塞腔內(nèi)油液壓力分別為4.66、4.62 MPa。在下死點(diǎn)和上死點(diǎn)位置附近區(qū)域,柱塞腔油液壓力波動幅度分別為3.307和3.246 MPa。柱塞腔與高壓油孔接通階段,柱塞腔油液壓力波動幅度分別為0.324、0.322 MPa,柱塞腔內(nèi)油液過渡平穩(wěn),配流副泄漏量小,馬達(dá)的容積效率高。

    2)采用42°與45°錐度角的配流副正常工作時,配流副部件受力狀況良好。當(dāng)錐度角為42°時,缸體和配流軸的最大等效應(yīng)力分別占屈服強(qiáng)度的22.92%、7.49%,最大變形分別為12.069、4.288 4m。當(dāng)錐度角為45°時,缸體和配流軸的最大等效應(yīng)力分別占屈服強(qiáng)度的21.18%、7.52%,最大變形分別為11.818、3.471 8m。二者最大等效應(yīng)力皆遠(yuǎn)小于其屈服強(qiáng)度,沒有強(qiáng)度失效的風(fēng)險,雖然有微量彈性變形,但對配流副的正常工作影響極小。

    3)針對試驗(yàn)用摩擦副材料,相較于45°錐度角,42°錐度角摩擦副平均摩擦系數(shù)(0.155)小,且波動小、穩(wěn)定性好。上、下試件的磨損率較45°錐度角分別小1.96×10-6、7.601×10-6mm3/(N·mm),有利于能量回收馬達(dá)配流副的穩(wěn)定工作及高效運(yùn)行。

    [1] 李文華,劉嬌,柴博. 節(jié)能液壓泵模糊PID控制系統(tǒng)研究與仿真[J]. 控制工程,2017,24(7):1347-1351.

    Li Wenhua, Liu Jiao, Chai Bo. Research and simulation of energy efficient fuzzy PID control system for hydraulic pump stations[J]. Control Engineering of China, 2017, 24(7): 1347-1351. (in Chinese with English abstract)

    [2] Hu D M, Xia C L, Xu B, et al. Computational simulation and experimental research on the speed control of the closed-circuit energy-saving hydraulic elevator system[C]. International Symposium on Fluid Power Transmission & Control, Qinhuangdao, China, 2007. 6, 574-579.

    [3] 孟志明,賀元成,康帥帥,等. 液壓節(jié)能技術(shù)的現(xiàn)狀與發(fā)展[J]. 機(jī)械工程師,2014(1):114-117.

    Meng Zhiming, He Yuancheng, Kang Shuaishuai, et al. Status and development of hydraulic energy-saving technology[J]. Mechanical Engineer, 2014(1): 114-117. (in Chinese with English abstract)

    [4] 葛志偉. 探討工程機(jī)械液壓節(jié)能技術(shù)現(xiàn)狀與發(fā)展[J]. 液壓氣動與密封,2021,41(10):70-73.

    Ge Zhiwei. The status and development of energy conservation of hydraulic system for construction machinery[J]. Hydraulics Pneumatics & Seals, 2021, 41(10): 70-73. (in Chinese with English abstract)

    [5] 王煜. 柱塞式能量回收裝置的研制[D]. 蘭州:蘭州理工大學(xué),2021.

    Wang Yu. Development of Plunger Recovery Device[D]. Lanzhou: Lanzhou University of Technology, 2021. (in Chinese with English abstract)

    [6] 冀宏. 液壓氣壓傳動與控制[M]. 武漢:華中科技大學(xué)出版社,2009.

    [7] 許耀銘. 油膜理論與液壓泵和馬達(dá)的摩擦副設(shè)計[M]. 北京:機(jī)械工業(yè)出版社,1987.

    [8] Totten E G. Handbook of Hydraulic Fluid Technology[M]. Tarrytown, New York: Union Carbide Corporation, 2000.

    [9] 王廣達(dá). 海水淡化泵—馬達(dá)增壓能量回收集成裝置研究及配流副優(yōu)化[D]. 秦皇島:燕山大學(xué),2021.

    Wang Guangda. Research on Desalination Pump-motor Pressure Boosting and Energy Recovery Intergrated Machine and Optimization of Port Plate Pairs[D]. Qinhuangdao: Yanshan University, 2021. (in Chinese with English abstract)

    [10] 喬豐立,盧堃,李桂花. 錐形軸配流副水介質(zhì)軸向柱塞泵的研究[J]. 液壓與氣動,2007,31(10):68-70.

    Qiao Fengli, Lu Kun, Li Guihua. Study of conical spindle distributing flow pairs for water axial piston pump[J]. Chinese Hydraulics & Pneumatics, 2007, 31(10): 68-70. (in Chinese with English abstract)

    [11] 喬豐立,盧堃,李桂花,等. 水介質(zhì)液壓泵錐形軸配流副軸承間隙的計算[J]. 蘭州理工大學(xué)學(xué)報,2007,33(4):59-62.

    Qiao Fengli, Lu Kun, Li Guihua, et al. Calculation of bearing clearance of conical axle flow-distribution pairs in hydraulic pumps with water medium[J]. Journal of Lanzhou University of Technology, 2007, 33(4): 59-62. (in Chinese with English abstract)

    [12] 江浩. 一種錐形柱塞結(jié)構(gòu)的斜軸式定量馬達(dá)的研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué),2017.

    Jiang Hao. Research on Inclined Shaft Quantitative Motor with a Conical Plunger Structure[D]. Harrbin: Harbin Institute of Technology, 2017. (in Chinese with English abstract)

    [13] 荊崇波,鄒云飛,苑士華,等. 彈性變形對錐形配流副潤滑特性的影響[J]. 潤滑與密封,2009,34(8):28-32.

    Jing Chongbo, Zou Yunfei, Yuan Shihua, et al. The influence of elastic deformation lubrication characteristics of conical port plate pair of axial piston pump[J]. Lubrication Engineering, 2009, 34(8): 28-32. (in Chinese with English abstract)

    [14] 荊崇波,周俊杰,劉建豪,等. 球活塞式液壓泵錐形配流副自適應(yīng)潤滑特性[J]. 哈爾濱工業(yè)大學(xué)學(xué)報,2019,51(1):162-169.

    Jing Chongbo, Zhou Junjie, Liu Jianhao, et al. Self-adapting lubrication characteristics of conical coupling of ball piston hydraulic pump[J]. Journal of Harbin Institute of Technology, 2019, 51(1): 162-169. (in Chinese with English abstract)

    [15] 鄒云飛,荊崇波,劉晉川. 錐形配流副彈流潤滑模型與數(shù)值計算方法[J]. 農(nóng)業(yè)工程學(xué)報,2010,26(12):136-141.

    Zou Yunfei, Jing Chongbo, Liu Jinchuan. Lubricating model of conical spindle distributing flow pair and arithmetic of model[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(12): 136-141. (in Chinese with English abstract)

    [16] 韋曉磊,鄒云飛,荊崇波,等. 錐形配流副靜壓支承的節(jié)流阻尼效應(yīng)研究[J]. 液壓氣動與密封,2013,33(4):40-42.

    Wei Xiaolei, Zou Yunfei, Jing Chongbo, et al. Research of restriction orifice effect in hydraulic static bearing for conical spindle distributing flow pair[J]. Hydraulics Pneumatics & Seals, 2013, 33(4): 40-42. (in Chinese with English abstract)

    [17] Zhou J J, Jing C B, Liu J H. On the steady posture of conical spindle distribution used in ball piston pump[J]. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2020, 234(4): 490-499.

    [18] 申正精,楚武利,董瑋. 顆粒參數(shù)對螺旋離心泵流場及過流部件磨損特性的影響[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(6):58-66.

    Shen Zhengjing, Chu Wuli, Dong Wei. Effect of particle parameters on flow field and erosion wear characteristics of flow passage components in screw centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(6): 58-66. (in Chinese with English abstract)

    [19] 于勇,張俊明,姜連田. FLUENT入門與進(jìn)階教程[M]. 北京:北京理工大學(xué)出版社,2013.

    [20] 艾青林,周華,楊華勇. 軸向柱塞泵配流副潤滑特性的數(shù)值分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2004,35(6):78-81, 85.

    Ai Qinglin, Zhou Hua, Yang Huayong. Numerical analysis of lubrication property for port plate pair of axial piston pump[J]. Transactions of the Chinese Society for Agricultural Machinery, 2004, 35(6): 78-81, 85. (in Chinese with English abstract)

    [21] 羅惕乾. 流體力學(xué)(第四版)[M]. 北京:機(jī)械工業(yè)出版社,2017.

    [22] 王福軍. 計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用[M]. 北京:清華大學(xué)出版社,2004.

    [23] 郭鵬程,孫龍剛,羅興锜. 混流式水輪機(jī)葉道渦流動特性研究[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(20):43-51.

    Guo Pengcheng, Sun Longgang, Luo Xingqi. Flow characteristic investigation into inter-blade vortex for Francis turbine[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(20): 43-51. (in Chinese with English abstract)

    [24] 聞德生,商旭東,顧攀,等. 雙定子擺動液壓馬達(dá)泄漏與容積效率分析及密封改進(jìn)[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(12):74-81.

    Wen Desheng, Shang Xudong, Gu Pan, et al. Analysis of leakage and volumetric efficiency and seal improvement for double-stator swing hydraulic motor[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 74-81. (in Chinese with English abstract)

    [25] 劉厚林,徐歡,吳賢芳,等. 流固耦合作用對離心泵內(nèi)外特性的影響[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(13):82-87.

    Liu Houlin, Xu Huan, Wu Xianfang, et al. Effect of fluid-structure interaction on internal and external characteristics of centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(13): 82-87. (in Chinese with English abstract)

    [26] 王洋,王洪玉,張翔,等. 基于流固耦合理論的離心泵沖壓焊接葉輪強(qiáng)度分析[J]. 農(nóng)業(yè)工程學(xué)報,2011,27(3):131-136.

    Wang Yang, Wang Hongyu, Zhang Xiang, et al. Strength analysis on the stamping and welding impeller in centrifugal pump based on fluid-structure interaction theorem[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(3): 131-136. (in Chinese with English abstract)

    [27] 周華,焦素娟,李家鑫. 純水柱塞泵摩擦副材料摩擦磨損篩選試驗(yàn)[J]. 液壓與氣動,2002,26(7):40-41.

    Zhou Hua, Jiao Sujuan, Li Jiaxin. Material selection of friction pairs for water hydraulic piston pump[J]. Chinese Hydraulics & Pneumatics, 2002, 26(7): 40-41. (in Chinese with English abstract)

    [28] 李少年,包尚令,楊攀,等. 軸向柱塞泵配流盤摩擦副材料的磨損實(shí)驗(yàn)研究[J]. 液壓氣動與密封,2021,41(7):1-3,8.

    Li Shaonian, Bao Shangling, Yang Pan, et al. Experimental study on friction pairs wear of axial piston pump port plate[J]. Hydraulics Pneumatics & Seals, 2021, 41(7): 1-3, 8. (in Chinese with English abstract)

    [29] 李少年,張磊,常露丹,等. 固體顆粒對高壓葉片泵配流副油膜特性影響的數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(19):38-49.

    Li Shaonian, Zhang Lei, Chang Ludan, et al. Numerical simulation of effects of solid particles on oil film characteristics of port pair of high-pressure vane pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(19): 38-49. (in Chinese with English abstract)

    [30] 劉浩. 含磨粒油液子母葉片泵摩擦副摩擦特性研究[D]. 蘭州:蘭州理工大學(xué),2020.

    Liu Hao. Friction Characteristics of Friction Pairs of the Intra-vane Type Pump with Abrasive Grains in the Oil[D]. Lanzhou: Lanzhou University of Technology, 2020. (in Chinese with English abstract)

    [31] Li S N, Liu H, Li Y, et al. Tribological behavior of friction pairs in a high-pressure vane pump under variable working conditions[J]. Tribology Online, 2021, 16(3): 159-169.

    [32] Zhang Z H, Nie S L, Liao W J, et al. Tribological behaviors of carbon fiber reinforced polyetheretherketone sliding against silicon carbide ceramic under seawater lubrication[J]. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2014, 228(1): 1421 -1432.

    Analysis and tests of the taper angle of flow distribution pair in a piston type energy recovery motor

    Li Shaonian1, Liang Tao1, Bao Shangling1, Wang Yu2, Zhou Zhiyuan1, Li Xi1

    (1.,,730050,;2.,,066099,)

    A piston-type energy recovery motor is integrated the hydraulic motor and generator into a new generation of hydraulic energy recovery device. Among them, the flow distribution (friction) pair is composed of the cylinder block and flow distribution shaft. The taper angle of the friction surface on the flow distribution pair can pose an important impact on the oil flow distribution, load-bearing, friction and wear characteristics of the motor. In this study, an optimal combination of the taper angle was explored using the theoretical analysis, numerical simulation, and experimental test. Firstly, the available range of a taper angle was determined, according to the structure and size of the flow distribution pair in the motor. Four taper angles of 36°, 39°, 42°, and 45°were then selected as the research objects. Secondly, the fluid domain models were established to grid the four taper distribution pairs. The pressure contours and velocity vector diagrams were obtained for the four taper angles, in order to analyze the oil pressure fluctuation in the piston chamber and the leakage in the oil film. The solid domain grid models of distribution pairs were achieved at the taper angles of 42°and 45°and then to evaluate the stress and deformation of the cylinder block and the flow distribution shaft. Finally, the friction test pieces were made according to the motion and friction surface between the cylinder and the flow distribution shaft for the friction test. Specifically, the materials of the upper and lower test pieces were matched with those of the cylinder block and flow distribution shaft. The friction coefficient and wear amount of pairs were obtained at the taper angles of 42°and 45°during friction and wear experiments. The result showed that the oil pressures were 4.66 and 4.62 MPa in the piston chamber located in the high-pressure dead point at the taper angles of 42°and 45°, respectively, whereas, the amplitudes of pressure fluctuation were 3.307 and 3.246 MPa, respectively. Furthermore, the amplitudes of oil pressure fluctuation in the piston chamber were 0.324 and 0.322 MPa, respectively, in the process of connection between the piston chamber and the high pressure oil hole. There were the smaller amplitudes in the three-group data among the four taper angles. As such, the higher volume efficiency of the motor was achieved in the smaller leakage of the flow distribution pair, where the oil was transited smoothly at the taper angles of 42°and 45°in the piston chamber. The maximum equivalent stress of the flow distribution shaft was much less than the yield strength at the taper angle of 42°and 45°. There was 0.74 percentage points larger proportion of the maximum equivalent stress to the yield strength of cylinder block at the taper angle of 42°, compared with the taper angle of 45°. The maximum deformation of cylinder block was 0.251m. There was no failure risk of strength invalidation for the components of flow distribution pair at the two taper angles. Only a slightly impact was observed on the normal operation of the flow distribution pair, although the slight elastic deformation existed in this case. The average friction coefficient was less than 0.012 at the taper angle of 42°, indicating the small fluctuation and excellent stability. More importantly, the wear rates of the upper and lower pieces were 1.966×10-6and 7.601×10-6mm3/(N·mm) less at the taper angle of 42°than those of 45°, respectively. Therefore, the taper angle of 42° was conducive to the stable and efficient operation of the energy recovery motor. The finding can provide a strong reference to design the flow distribution pair of piston-type energy recovery motor.

    motors; flow field; simulation; stresses; deformation; flow distribution pair; taper angle; friction and wear

    10.11975/j.issn.1002-6819.2022.20.003

    TH137

    A

    1002-6819(2022)-20-0020-10

    李少年,梁濤,包尚令,等. 柱塞式能量回收馬達(dá)配流副錐度角的分析與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(20):20-29.doi:10.11975/j.issn.1002-6819.2022.20.003 http://www.tcsae.org

    Li Shaonian, Liang Tao, Bao Shangling, et al. Analysis and tests of the taper angle of flow distribution pair in a piston type energy recovery motor[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(20): 20-29. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.20.003 http://www.tcsae.org

    2022-07-27

    2022-09-27

    國家自然科學(xué)基金項(xiàng)目(52165006)

    李少年,博士,副教授,博士生導(dǎo)師,研究方向?yàn)榱黧w傳動與流體測控技術(shù)。Email:lsn19@163.com

    猜你喜歡
    配流角為錐度
    重載柱塞泵球面配流副承載特性研究*
    微觀織構(gòu)配流副熱-流-固耦合潤滑特性
    次級線圈繞線錐度對LVDT靜態(tài)特性的影響
    高速鋼電火花線切割錐度加工試驗(yàn)研究
    等效錐度曲線非線性特性及影響研究
    無鍵錐度液壓聯(lián)接在大功率高轉(zhuǎn)速偶合器中的應(yīng)用
    電力與能源(2017年6期)2017-05-14 06:19:48
    鐵路編組站動態(tài)配流分層模型
    立體幾何測試卷(B卷)
    一種以雙環(huán)鉑為有效成分的藥物檢測方法
    利用解的性質(zhì)解題
    午夜免费男女啪啪视频观看| 午夜日韩欧美国产| 中文字幕高清在线视频| 三上悠亚av全集在线观看| 欧美老熟妇乱子伦牲交| 免费高清在线观看视频在线观看| 免费久久久久久久精品成人欧美视频| 亚洲欧美中文字幕日韩二区| 欧美97在线视频| 美女脱内裤让男人舔精品视频| 免费人妻精品一区二区三区视频| 夫妻午夜视频| 亚洲熟女毛片儿| av欧美777| 久久精品国产亚洲av高清一级| 欧美日韩一级在线毛片| 在线观看免费高清a一片| 久久精品久久久久久久性| 另类亚洲欧美激情| 国产精品 国内视频| 久久久久国产一级毛片高清牌| av国产久精品久网站免费入址| 纯流量卡能插随身wifi吗| 熟女av电影| 亚洲激情五月婷婷啪啪| 午夜久久久在线观看| 老司机亚洲免费影院| 极品少妇高潮喷水抽搐| 天堂中文最新版在线下载| 操美女的视频在线观看| 欧美 亚洲 国产 日韩一| 好男人电影高清在线观看| 亚洲国产日韩一区二区| bbb黄色大片| 中文字幕av电影在线播放| 激情视频va一区二区三区| 蜜桃在线观看..| 男女边摸边吃奶| 男女边吃奶边做爰视频| 建设人人有责人人尽责人人享有的| 少妇被粗大的猛进出69影院| 久久九九热精品免费| 国产精品一区二区在线不卡| 999精品在线视频| 亚洲男人天堂网一区| 日韩大码丰满熟妇| 亚洲 欧美一区二区三区| 超碰97精品在线观看| 人人妻人人爽人人添夜夜欢视频| 制服诱惑二区| 91精品伊人久久大香线蕉| 亚洲国产成人一精品久久久| 亚洲人成77777在线视频| 亚洲国产av新网站| 99久久99久久久精品蜜桃| 亚洲av成人不卡在线观看播放网 | 成人亚洲精品一区在线观看| 欧美精品一区二区免费开放| 老鸭窝网址在线观看| 狠狠精品人妻久久久久久综合| 一区二区日韩欧美中文字幕| 久久精品人人爽人人爽视色| 一级毛片黄色毛片免费观看视频| 女人久久www免费人成看片| 欧美人与善性xxx| 欧美日韩综合久久久久久| 精品福利永久在线观看| 亚洲五月婷婷丁香| 久久这里只有精品19| 欧美日韩亚洲国产一区二区在线观看 | 飞空精品影院首页| 精品亚洲成a人片在线观看| 久久精品熟女亚洲av麻豆精品| 性高湖久久久久久久久免费观看| 国产欧美日韩综合在线一区二区| 男女高潮啪啪啪动态图| 亚洲精品久久午夜乱码| 亚洲色图综合在线观看| 黑人猛操日本美女一级片| 国产1区2区3区精品| 99国产精品一区二区三区| 性色av一级| 国产日韩欧美亚洲二区| 两个人免费观看高清视频| 黄色一级大片看看| 满18在线观看网站| 亚洲天堂av无毛| 欧美日韩视频精品一区| e午夜精品久久久久久久| 中文字幕最新亚洲高清| 美国免费a级毛片| 丝袜在线中文字幕| 韩国高清视频一区二区三区| 亚洲伊人久久精品综合| 亚洲天堂av无毛| 大码成人一级视频| 精品少妇久久久久久888优播| 欧美 亚洲 国产 日韩一| 免费看av在线观看网站| 久久国产精品男人的天堂亚洲| 国产一区二区三区综合在线观看| 国产成人精品久久久久久| 少妇人妻 视频| 久久ye,这里只有精品| 美国免费a级毛片| 国产成人欧美| 亚洲精品久久成人aⅴ小说| 最近中文字幕2019免费版| 久久精品亚洲av国产电影网| 午夜福利视频精品| 日本欧美视频一区| 国产色视频综合| 91精品国产国语对白视频| 男女午夜视频在线观看| 亚洲一码二码三码区别大吗| 黄色怎么调成土黄色| 伊人久久大香线蕉亚洲五| 精品国产乱码久久久久久男人| 亚洲国产欧美在线一区| 妹子高潮喷水视频| 久久女婷五月综合色啪小说| 国产亚洲午夜精品一区二区久久| 飞空精品影院首页| 大话2 男鬼变身卡| 丰满迷人的少妇在线观看| 欧美亚洲日本最大视频资源| 肉色欧美久久久久久久蜜桃| 丰满迷人的少妇在线观看| 亚洲国产成人一精品久久久| 丝袜美腿诱惑在线| 大片电影免费在线观看免费| 一级a爱视频在线免费观看| 亚洲,欧美,日韩| 国产国语露脸激情在线看| 一本综合久久免费| 国产不卡av网站在线观看| 成人黄色视频免费在线看| 午夜福利视频精品| 午夜福利一区二区在线看| 久久性视频一级片| 一本综合久久免费| 国产又色又爽无遮挡免| 亚洲国产av新网站| 水蜜桃什么品种好| 亚洲,欧美,日韩| 菩萨蛮人人尽说江南好唐韦庄| 我要看黄色一级片免费的| 久热这里只有精品99| 久久性视频一级片| 欧美人与性动交α欧美软件| 亚洲精品美女久久久久99蜜臀 | 亚洲精品在线美女| 美女脱内裤让男人舔精品视频| 2018国产大陆天天弄谢| 悠悠久久av| 女性生殖器流出的白浆| 狂野欧美激情性xxxx| 一级黄片播放器| 亚洲精品美女久久av网站| 久久久久久亚洲精品国产蜜桃av| 国产精品一二三区在线看| 成人亚洲精品一区在线观看| 搡老乐熟女国产| 18禁观看日本| 日韩av不卡免费在线播放| 高潮久久久久久久久久久不卡| 亚洲精品一区蜜桃| 1024视频免费在线观看| 天天躁日日躁夜夜躁夜夜| 亚洲精品久久成人aⅴ小说| 日本色播在线视频| 男人爽女人下面视频在线观看| 亚洲av美国av| 国产亚洲欧美精品永久| 免费在线观看影片大全网站 | 色视频在线一区二区三区| 精品人妻1区二区| 99国产精品一区二区蜜桃av | 国产高清不卡午夜福利| 久久久精品94久久精品| 涩涩av久久男人的天堂| 69精品国产乱码久久久| 国产黄色视频一区二区在线观看| 欧美日韩综合久久久久久| 久久久精品94久久精品| 波多野结衣av一区二区av| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 巨乳人妻的诱惑在线观看| 国产精品偷伦视频观看了| 王馨瑶露胸无遮挡在线观看| 日韩欧美一区视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 日韩 亚洲 欧美在线| 欧美日本中文国产一区发布| 美女扒开内裤让男人捅视频| 亚洲av在线观看美女高潮| 国产精品 国内视频| 啦啦啦视频在线资源免费观看| 日本欧美视频一区| 菩萨蛮人人尽说江南好唐韦庄| 国产黄色视频一区二区在线观看| 欧美日韩综合久久久久久| 久久ye,这里只有精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲五月婷婷丁香| a 毛片基地| 19禁男女啪啪无遮挡网站| 国产一区二区在线观看av| 国产老妇伦熟女老妇高清| 欧美黑人欧美精品刺激| 男女高潮啪啪啪动态图| 免费少妇av软件| 午夜av观看不卡| 在现免费观看毛片| 久久精品人人爽人人爽视色| 青春草亚洲视频在线观看| 国产爽快片一区二区三区| 欧美精品一区二区大全| 大片免费播放器 马上看| 日韩欧美一区视频在线观看| 欧美日韩福利视频一区二区| 日韩免费高清中文字幕av| 久久精品国产a三级三级三级| 欧美精品人与动牲交sv欧美| 一级毛片女人18水好多 | 一级片'在线观看视频| 侵犯人妻中文字幕一二三四区| 大片免费播放器 马上看| 午夜福利乱码中文字幕| 中文字幕色久视频| 亚洲av在线观看美女高潮| 欧美亚洲 丝袜 人妻 在线| 99久久综合免费| 国产一区有黄有色的免费视频| 久久人妻福利社区极品人妻图片 | 又黄又粗又硬又大视频| 麻豆乱淫一区二区| 成人午夜精彩视频在线观看| 母亲3免费完整高清在线观看| 日韩大片免费观看网站| 久久精品国产综合久久久| 国产日韩欧美亚洲二区| 黄频高清免费视频| 男女边吃奶边做爰视频| 亚洲中文字幕日韩| 十八禁网站网址无遮挡| 一区在线观看完整版| 久久久久久久久免费视频了| 麻豆乱淫一区二区| 啦啦啦 在线观看视频| 视频区图区小说| av电影中文网址| 亚洲国产毛片av蜜桃av| 亚洲av日韩在线播放| √禁漫天堂资源中文www| 一区二区三区激情视频| 成人亚洲精品一区在线观看| 黄色a级毛片大全视频| 日韩 欧美 亚洲 中文字幕| 国产成人av教育| 国产精品国产三级国产专区5o| 午夜福利在线免费观看网站| 校园人妻丝袜中文字幕| 久久久精品免费免费高清| 精品一区二区三区四区五区乱码 | 欧美精品一区二区大全| 国产免费一区二区三区四区乱码| 亚洲第一青青草原| 免费av中文字幕在线| av在线app专区| 男人操女人黄网站| 性色av一级| 成人亚洲欧美一区二区av| 色综合欧美亚洲国产小说| 久久综合国产亚洲精品| 热re99久久国产66热| 精品亚洲成国产av| 日韩av在线免费看完整版不卡| xxxhd国产人妻xxx| 亚洲成人免费av在线播放| 人成视频在线观看免费观看| av在线播放精品| 久久精品国产a三级三级三级| av网站免费在线观看视频| 亚洲情色 制服丝袜| 新久久久久国产一级毛片| 电影成人av| 国产成人a∨麻豆精品| 男人添女人高潮全过程视频| 汤姆久久久久久久影院中文字幕| 99国产精品99久久久久| 人妻 亚洲 视频| 色精品久久人妻99蜜桃| xxxhd国产人妻xxx| tube8黄色片| a级毛片在线看网站| av国产精品久久久久影院| 国产xxxxx性猛交| 久9热在线精品视频| 亚洲自偷自拍图片 自拍| 黄片小视频在线播放| 中文字幕人妻熟女乱码| 亚洲国产精品成人久久小说| 亚洲成人国产一区在线观看 | 欧美日韩黄片免| 精品免费久久久久久久清纯 | 国产亚洲精品第一综合不卡| 免费人妻精品一区二区三区视频| 亚洲精品国产av成人精品| 国产av一区二区精品久久| 国产一区二区在线观看av| 韩国高清视频一区二区三区| 老司机影院毛片| 欧美少妇被猛烈插入视频| 国产一级毛片在线| 国产福利在线免费观看视频| 一区二区三区四区激情视频| 精品一品国产午夜福利视频| 亚洲欧洲国产日韩| 国产片特级美女逼逼视频| www.999成人在线观看| 国产麻豆69| 亚洲午夜精品一区,二区,三区| www.999成人在线观看| 人人妻人人澡人人看| 又黄又粗又硬又大视频| av视频免费观看在线观看| 青草久久国产| 99久久综合免费| av欧美777| 久久性视频一级片| 99久久综合免费| 久久亚洲精品不卡| 99热国产这里只有精品6| 一本综合久久免费| 另类亚洲欧美激情| 免费观看a级毛片全部| 亚洲自偷自拍图片 自拍| a级片在线免费高清观看视频| 久久这里只有精品19| 一边摸一边抽搐一进一出视频| 日本黄色日本黄色录像| 色婷婷av一区二区三区视频| 一区二区三区激情视频| 高清黄色对白视频在线免费看| 91老司机精品| 一区二区日韩欧美中文字幕| 男女国产视频网站| 1024香蕉在线观看| 黄片播放在线免费| 91精品三级在线观看| 黄片播放在线免费| 亚洲精品美女久久av网站| 我的亚洲天堂| 大香蕉久久网| 色播在线永久视频| 99精国产麻豆久久婷婷| 欧美人与性动交α欧美软件| 亚洲一码二码三码区别大吗| 日韩 亚洲 欧美在线| 一级片免费观看大全| 国产精品一区二区精品视频观看| 久久精品成人免费网站| 亚洲精品在线美女| 99香蕉大伊视频| 在线亚洲精品国产二区图片欧美| 国产成人一区二区三区免费视频网站 | 精品一品国产午夜福利视频| 一级毛片电影观看| 国产精品国产av在线观看| 国产精品三级大全| 1024香蕉在线观看| xxxhd国产人妻xxx| 国产免费现黄频在线看| 一区二区三区乱码不卡18| 国产男女超爽视频在线观看| 亚洲欧美清纯卡通| 欧美在线一区亚洲| 欧美成人午夜精品| 亚洲美女黄色视频免费看| 亚洲中文字幕日韩| 婷婷色综合大香蕉| 亚洲国产精品国产精品| 中文字幕色久视频| 日日夜夜操网爽| 男女无遮挡免费网站观看| 色综合欧美亚洲国产小说| 丝袜喷水一区| 男人操女人黄网站| 国产精品免费视频内射| 精品久久久精品久久久| 男女午夜视频在线观看| 视频在线观看一区二区三区| 久久精品国产综合久久久| 1024香蕉在线观看| av又黄又爽大尺度在线免费看| 欧美日韩亚洲高清精品| 欧美精品人与动牲交sv欧美| 亚洲国产最新在线播放| 亚洲第一青青草原| 免费久久久久久久精品成人欧美视频| 老熟女久久久| 久久九九热精品免费| 亚洲午夜精品一区,二区,三区| 人人妻人人澡人人爽人人夜夜| 久久久久精品国产欧美久久久 | 美国免费a级毛片| 看免费成人av毛片| 久久影院123| 国产精品麻豆人妻色哟哟久久| 免费观看a级毛片全部| 久久久精品区二区三区| 久久久亚洲精品成人影院| 高清欧美精品videossex| 岛国毛片在线播放| xxxhd国产人妻xxx| 免费高清在线观看日韩| 国产老妇伦熟女老妇高清| a级毛片在线看网站| 国产日韩欧美视频二区| 国产熟女午夜一区二区三区| 亚洲伊人色综图| 一本色道久久久久久精品综合| 国产又爽黄色视频| 搡老乐熟女国产| 十八禁网站网址无遮挡| 国产1区2区3区精品| 亚洲九九香蕉| 曰老女人黄片| 亚洲美女黄色视频免费看| 999久久久国产精品视频| 日韩制服骚丝袜av| 日韩 欧美 亚洲 中文字幕| av一本久久久久| 在线精品无人区一区二区三| 99九九在线精品视频| 大话2 男鬼变身卡| 久久人人爽人人片av| 久久九九热精品免费| 国产爽快片一区二区三区| 男女边摸边吃奶| 亚洲精品成人av观看孕妇| 亚洲天堂av无毛| 亚洲av在线观看美女高潮| 大码成人一级视频| 亚洲成人国产一区在线观看 | 香蕉丝袜av| 亚洲欧美清纯卡通| 伊人亚洲综合成人网| 久久久精品区二区三区| 精品国产一区二区久久| 午夜福利视频在线观看免费| 男女国产视频网站| 日韩,欧美,国产一区二区三区| 精品高清国产在线一区| 男女下面插进去视频免费观看| 国产日韩欧美视频二区| 国产成人一区二区三区免费视频网站 | 视频区欧美日本亚洲| 午夜91福利影院| 男女下面插进去视频免费观看| 日韩一区二区三区影片| 香蕉国产在线看| 999久久久国产精品视频| 久久久久视频综合| 欧美老熟妇乱子伦牲交| 丝袜美足系列| 久久久久久久久久久久大奶| 飞空精品影院首页| 久久久久视频综合| 亚洲,欧美精品.| 久久 成人 亚洲| 肉色欧美久久久久久久蜜桃| 国产日韩欧美视频二区| 日韩制服骚丝袜av| 日韩欧美一区视频在线观看| 久久精品久久精品一区二区三区| 狠狠婷婷综合久久久久久88av| 777久久人妻少妇嫩草av网站| 两性夫妻黄色片| 黑人巨大精品欧美一区二区蜜桃| 国产精品免费大片| 欧美黑人欧美精品刺激| 一区二区三区乱码不卡18| 久久久亚洲精品成人影院| 晚上一个人看的免费电影| 各种免费的搞黄视频| 国产欧美日韩一区二区三 | 不卡av一区二区三区| 男女下面插进去视频免费观看| 啦啦啦啦在线视频资源| 亚洲免费av在线视频| 午夜福利一区二区在线看| 黄色怎么调成土黄色| 人体艺术视频欧美日本| 国产精品免费视频内射| 亚洲欧洲国产日韩| 午夜免费观看性视频| 一区二区三区四区激情视频| 国产福利在线免费观看视频| 老司机靠b影院| 一本—道久久a久久精品蜜桃钙片| 少妇被粗大的猛进出69影院| 啦啦啦视频在线资源免费观看| 亚洲七黄色美女视频| 亚洲精品在线美女| 一本色道久久久久久精品综合| 亚洲精品一二三| 国产免费一区二区三区四区乱码| 日韩av在线免费看完整版不卡| 成年av动漫网址| 日本猛色少妇xxxxx猛交久久| 91精品三级在线观看| 我要看黄色一级片免费的| 亚洲五月婷婷丁香| 一个人免费看片子| 天堂8中文在线网| 国产亚洲精品久久久久5区| 成年人午夜在线观看视频| 1024视频免费在线观看| 久久国产精品影院| 男人操女人黄网站| av天堂在线播放| 午夜激情久久久久久久| 久久精品熟女亚洲av麻豆精品| 搡老乐熟女国产| 亚洲,欧美,日韩| 高潮久久久久久久久久久不卡| 国产欧美日韩精品亚洲av| 久久综合国产亚洲精品| 亚洲精品美女久久久久99蜜臀 | 久久鲁丝午夜福利片| 免费在线观看日本一区| www.av在线官网国产| 国产免费视频播放在线视频| 国产免费现黄频在线看| 91九色精品人成在线观看| 性色av乱码一区二区三区2| 亚洲午夜精品一区,二区,三区| 久久久亚洲精品成人影院| 韩国高清视频一区二区三区| 亚洲天堂av无毛| 亚洲 国产 在线| 一本—道久久a久久精品蜜桃钙片| 久久国产精品男人的天堂亚洲| 熟女少妇亚洲综合色aaa.| 国产不卡av网站在线观看| 欧美日韩视频高清一区二区三区二| 在线天堂中文资源库| 精品少妇内射三级| 欧美日韩亚洲综合一区二区三区_| 手机成人av网站| 亚洲欧洲国产日韩| 亚洲视频免费观看视频| 国产在视频线精品| 日本av免费视频播放| 久久国产精品男人的天堂亚洲| 啦啦啦中文免费视频观看日本| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 国产精品一区二区在线不卡| 狠狠婷婷综合久久久久久88av| 国产一级毛片在线| 日本猛色少妇xxxxx猛交久久| 久久精品人人爽人人爽视色| 麻豆乱淫一区二区| 亚洲精品乱久久久久久| 国产欧美日韩精品亚洲av| 曰老女人黄片| 亚洲国产av影院在线观看| 国产精品国产av在线观看| 一区二区三区乱码不卡18| 日本av免费视频播放| 青春草视频在线免费观看| 在线观看免费高清a一片| 18在线观看网站| 欧美精品高潮呻吟av久久| 一区二区三区四区激情视频| 日韩av不卡免费在线播放| 涩涩av久久男人的天堂| 精品少妇久久久久久888优播| 亚洲av欧美aⅴ国产| 精品亚洲乱码少妇综合久久| 人人澡人人妻人| 国产成人精品久久二区二区91| 精品人妻熟女毛片av久久网站| 2021少妇久久久久久久久久久| 一区在线观看完整版| 精品视频人人做人人爽| 啦啦啦在线免费观看视频4| 免费观看人在逋| 少妇 在线观看| 美女主播在线视频| 在线观看人妻少妇| 欧美日韩一级在线毛片| 极品少妇高潮喷水抽搐| 自拍欧美九色日韩亚洲蝌蚪91| 热99久久久久精品小说推荐| 亚洲 欧美一区二区三区| 国产在线免费精品| 美国免费a级毛片| 麻豆av在线久日| 中文字幕制服av| 国产成人免费观看mmmm| 在线观看免费视频网站a站| 91老司机精品| 99国产综合亚洲精品| 美国免费a级毛片| 亚洲av国产av综合av卡| 天天躁夜夜躁狠狠躁躁| 精品欧美一区二区三区在线| 欧美老熟妇乱子伦牲交| 精品少妇内射三级| 天天操日日干夜夜撸|