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

    發(fā)動機泄氣制動數(shù)值模擬及其性能優(yōu)化研究

    2014-08-08 13:44:56龔金科潘海杰余明果江杰田應華
    湖南大學學報·自然科學版 2014年5期
    關鍵詞:數(shù)值模擬發(fā)動機

    龔金科+潘海杰+余明果+江杰+田應華

    文章編號:16742974(2014)05000107

    收稿日期:20131127

    基金項目:國家自然科學基金資助項目(51276056); 湖南大學汽車車身先進設計制造國家重點實驗室自主項目(61075002)

    作者簡介:龔金科(1954-),男,湖南臨澧人,湖南大學教授,博士生導師

    通訊聯(lián)系人,E-mail:gongjinke@126.com;380188297@qq.com

    摘 要:為優(yōu)化發(fā)動機泄氣制動性能,用UG軟件建立了三維模型,應用動態(tài)網(wǎng)格技術對模型進行了網(wǎng)格劃分并對氣門和活塞的運動規(guī)律進行設置和定義.用Fluent軟件對泄氣制動工作過程進行三維瞬態(tài)數(shù)值模擬,分析了泄氣制動過程活塞壓力變化情況,并根據(jù)活塞瞬時受力得到各狀態(tài)下發(fā)動機平均制動功率,得出排氣背壓、排氣門開度、發(fā)動機轉速對制動性能的影響.通過優(yōu)化工作參數(shù),提高了泄氣制動性能,最后進行實驗驗證,結果表明最大誤差為5.2%,最小誤差僅為4.1%,表明模擬結果有較高精度和準確性.

    關鍵詞:發(fā)動機;動態(tài)網(wǎng)格技術;數(shù)值模擬;泄氣制動;制動性能

    中圖分類號:U461.3 文獻標識碼:A

    Research onthe Numerical Simulation and Performance 

    Optimization of Engine Bleeding Brake

    

    GONG Jinke,PAN Haijie,YU Mingguo,JIANG Jie,TIAN Yinghua

    (State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan Univ,

    Changsha, Hunan 410082,China)

    Abstract: In order to optimize the performance of engine bleeding brake, a threedimensional model was built with UG. The model was divided into grid by the dynamic mesh technique. Besides, the motion rule of the intake valve and piston was set and defined. By applying the threedimensional transient numerical simulation in the process of engine brake working by Fluent, the changes of pition pressure were analyzed. According to the piston's instantaneous force, the engine average braking power in different conditions was calculated, and then the influence of exhaust valve lash, engine speed and exhaust back pressure on engine braking performance was also obtained. By optimizing working parameters, the performance of engine bleeding brake was enhanced, and the simulation results were verified with experiments. The results have shown that the maximum error is 5.2%, and the smallest error is only 4.1%. It is shown that the simulation results are of high precision and accuracy.

    Key words:engine; dynamic mesh technique; numerical simulation; bleeding brake; braking performance

    

    為限制車速,車輛往往要頻繁使用制動器,造成制動鼓和制動片過熱,從而使制動效能降低,甚至會造成制動能力完全消失.所以,在現(xiàn)代商用車特別是大型商用車上增加輔助制動裝置已成為趨勢甚至成為硬性規(guī)定[1-3].泄氣制動作為一種重要的輔助制動系統(tǒng),在國外已有較多應用,大大增加了行車安全性.國內對泄氣制動的研究、生產(chǎn)和應用還處于起步階段[4-7].

    國內外對泄氣制動研究主要是通過數(shù)學軟件對泄氣制動的模型編寫程序進行模擬計算,且多為對缸內壓力及流場的分析,對泄氣制動性能缺乏系統(tǒng)的研究[8-11].本文運用動態(tài)網(wǎng)格技術對泄氣制動工作過程進行了三維瞬態(tài)數(shù)值模擬,這比傳統(tǒng)的數(shù)值計算更精確,不但顯示了活塞壓力的變化情況,而且得出了各狀態(tài)下平均制動功率的大小.通過模擬計算直觀系統(tǒng)地對泄氣制動性能進行了研究,為開發(fā)性能優(yōu)良的泄氣式輔助制動器提供參考依據(jù).

    1 發(fā)動機泄氣制動過程

    泄氣制動原理:通過在發(fā)動機配氣系統(tǒng)中加裝一套控制機構,使排氣門在發(fā)動機泄氣制動過程中始終保持一定的開度,

    氣體充量在壓縮和膨脹沖程時通過排氣門從氣缸排出,減少了缸內充量在膨脹沖程中對活塞的做功,從而增加了發(fā)動機的制動功率[11].圖1為發(fā)動機泄氣制動工作示意圖.

    圖1 泄氣式輔助制動工作過程示意圖

    Fig.1 Schematic diagram of bleeding brake work process

    2 發(fā)動機泄氣式輔助制動數(shù)學模型

    在泄氣制動過程中,缸內停止供油,不發(fā)生燃燒現(xiàn)象,缸內工質成分在一整個循環(huán)內基本保持不變.取氣缸為開放的熱力系統(tǒng),系統(tǒng)通過氣門及氣缸壁與外界進行物質與能量交換,當處于泄氣制動時,根據(jù)熱力學第一定律,有以下基本方程[12].

    質量守恒方程(連續(xù)性方程).

    單位時間內流通微元體中質量的增加,等于同一時間間隔內流入該微元體的凈質量:

    dm=∑jdmj (1)

    在泄氣制動過程中,質量守恒方程可表示為: 

    dmdφ=dmAdφ+dmEdφ+dmSdφ(2)

    式中mA為從進氣門流進氣缸的質量;mE為通過排氣門流出氣缸的質量;mS為通過活塞環(huán)間隙竄入曲軸箱內的質量.

    能量守恒方程:

    dU=dW+∑idQi+∑jhjdmj(3)

    式中U為系統(tǒng)內氣體內能;W為系統(tǒng)的有效功;Qi為通過系統(tǒng)邊界交換的能量;h為比焓;m為系統(tǒng)氣體質量;hjdmj為進入或帶出系統(tǒng)的能量.

    在泄氣制動過程中缸內系統(tǒng)與外界不斷發(fā)生物質與能量交換,能量守恒方程可具體表示如下[11].

    進氣及氣門疊開階段;壓縮、膨脹及排氣階段:

    dTdφ=1mcv[∑idQidφ-pdVdφ+(hA-u)dmAdφ+

    (hE-u)dmEdφ+(hS-u)dmSdφ](4)

    dTdφ=1mcv[∑idQidφ-pdVdφ+(hE-u)dmEdφ+

    (hS-u)dmSdφ](5)

    式中φ為曲軸轉角;cv為工質定容比熱;mE,mA 分別為流入、流出系統(tǒng)的質量;mS為通過活塞環(huán)間隙流出氣缸的質量;hE為進氣門前的比焓;hS為活塞環(huán)間隙工質的比焓;u為缸內氣體的比內能.

    泄氣制動時外部約束方程如下[11].

    氣缸瞬時容積為:

    V=VK2[2ε-1+1-cosφ+

    1λ(1-1-λ2sin2φ)](6)

    氣缸容積隨著曲軸轉角變化率為:

    dVdφ=Vh2(sin φ+λsin φcos φ1-λ2sin 2φ)(7)

    式中Vh為氣缸的工作容積,m3;ε為壓縮比;λ為曲柄連桿比.

    單位曲軸轉角的換熱量為[11]:

    ∑idQidφ=1ω∑3i=1αg?Ai(T-Twi)(8)

    式中ω為發(fā)動機角速度,ω=6n,n為轉速;αg為瞬時平均換熱系數(shù);A為換熱面積;T為氣缸內工質瞬時溫度;Tw為壁面的平均溫度;i=1為氣缸蓋,i=2為活塞,i=3為氣缸套.

    補充方程[12] (理想氣體):

    p=ρRT(9)

    由上述基本方程和約束方程聯(lián)立,可建立發(fā)動機泄氣制動工作的數(shù)學模型,由此可以求出氣缸內工質壓力隨曲軸轉角變化的規(guī)律.發(fā)動機泄氣制動工作時單缸循環(huán)的指示功為:

    W=∮pdV=∮FdS(10)

    3三維模型數(shù)值模擬計算

    3.1 三維數(shù)值模擬的流程

    三維數(shù)值模擬的主要流程如圖2所示.

    圖2 三維數(shù)值模擬流程圖

    Fig.2 Flow chart of threedimensional

    numerical simulation

    

    根據(jù)實物建立三維模型,對模型進行網(wǎng)格劃分,設置有關參數(shù)后進行三維數(shù)值模擬計算,用origin軟件進行數(shù)據(jù)處理,分析模擬結果,然后對參數(shù)進行優(yōu)化,最后進行實驗驗證,制定最佳方案.

    3.2 泄氣制動三維模型建立及動態(tài)網(wǎng)格劃分

    首先,根據(jù)某發(fā)動機的參數(shù)建立了氣缸簡化模型.在網(wǎng)格劃分中,應用動態(tài)網(wǎng)格技術對氣門和活塞的運動規(guī)律進行設置和定義,采用四面體網(wǎng)格作為過渡性網(wǎng)格,進排氣道部分采用四面體網(wǎng)格劃分,并對氣門座及氣閥附近區(qū)域網(wǎng)格進行加密;燃燒室采用穩(wěn)定性好和精度較高的六面體網(wǎng)格劃分,此區(qū)域為動網(wǎng)格區(qū).網(wǎng)格劃分如圖3所示.

    圖3 氣缸與氣道部分的網(wǎng)格

    Fig.3 Grid of cylinder and airway

    圖4為轉速為2 100 r/min、背壓為0.3 MPa、開度為0.7 mm、曲軸轉角為270°時活塞受力云圖.

    圖4 活塞受力云圖

    Fig.4 Piston stress nephogram

    3.3 泄氣制動性能評價方法

    本文通過模擬得到一個循環(huán)中各曲軸轉角下活塞瞬時壓力的大小,由活塞瞬時受力隨活塞位移的變化關系得到活塞平均制動功率,比較各狀態(tài)下活塞消耗功率大小可以反映出制動性能的好壞.圖5為轉速為2 100 r/min、背壓為0.3 MPa、開度為0.7 mm時活塞受力隨活塞位移變化關系圖.圖中F1為壓縮、膨脹過程中活塞做功大?。ù诉^程做負功),F2為進氣、排氣過程中活塞做功大小(此過程做負功).因此,一個循環(huán)中活塞做功大小為F1與F2的面積和.

    曲軸轉角n/(°)

    圖5 活塞所受壓力隨活塞位移變化關系

    Fig.5Change relationship ofpressure on the 

    piston and the piston displacement

    4 發(fā)動機泄氣制動模擬結果分析

    應用三維數(shù)值模擬對發(fā)動機泄氣制動工作過程進行模擬研究,定量分析相關參數(shù)對活塞受力及制動性能的影響,最終實現(xiàn)制動性能的優(yōu)化.

    4.1 泄氣制動過程活塞壓力研究

    圖6為轉速為2 100 r/min、背壓為0.3 MPa時,不同排氣門開度下活塞受力的比較.

    曲軸轉角n/(°)

    圖6 不同開度下活塞受力比較

    Fig.6 Piston stress comparison under 

    different opening

    從圖6可以看出:開度越大,活塞受力峰值越小.開度越大,壓力達到峰值時泄漏的氣體越多,缸內氣壓越低,活塞受力峰值越低.在不同開度下,活塞受力相差很大,開度在0.5 mm與1.0 mm時活塞壓力分別為53 786.25 N和37 068.69 N,受力相差達16 717.56 N.另外,壓力峰值出現(xiàn)在壓縮上止點前,并且隨著開度的增加,缸內壓力峰值位置逐漸遠離壓縮上止點,因為開度越大,泄漏的氣體越多,在壓縮上止點較遠時發(fā)動機進排氣效果達到平衡.

    圖7為背壓為0.3 MPa、開度為0.7 mm時,不同轉速下活塞受力的比較.從圖7可以看出:轉速越高,活塞壓力峰值越大.轉速越高,壓力達到峰值所用的時間越短,泄漏的氣體越少,缸內氣壓越高,活塞壓力峰值越大.在不同轉速下,活塞受力值相差很大,轉速在1 800 r/min與2 300 r/min時活塞壓力分別為40 325.64 N和51 932.03 N,受力相差12 123.5 N.并且隨著轉速的升高,相應缸內最大壓力位置越靠近壓縮上止點,這是因為轉速越高,單位曲軸轉角所經(jīng)歷的時間越短,從缸內溢出的空氣越少,壓力平衡出現(xiàn)得越晚.

    曲軸轉角n/(°)

    圖7 不同轉速下活塞受力比較

    Fig.7 Piston stress comparison under different speed

    圖8為轉速為2 100 r/min、開度為0.7 mm時,不同背壓下活塞受力的比較.從圖8可以看出:背壓越大,活塞壓力峰值越大.背壓越大,壓力達到峰值時排出的氣體越少,缸內氣壓越大,活塞所受壓力峰值越大.背壓在0.15 MPa與0.4 MPa時的活塞壓力分別為36 653.23 N和53 257.56 N,受力相差16 604.33 N,并且隨著背壓的增大,缸內最大壓力位置越接近壓縮上止點,因為背壓越大,排氣越困難,排氣阻力大,壓力平衡出現(xiàn)得越晚.

    曲軸轉角n/(°)

    圖8不同背壓下活塞受力比較

    Fig.8 Piston stress comparison under 

    different back pressure

    

    4.2 各因素對泄氣制動性能的影響分析

    在研究轉速、開度、背壓對泄氣制動性能的影響時,本文選取了如下參數(shù)點:轉速分別為1 800 r/min,1 900 r/min,2 000 r/min,2 100 r/min,2 200 r/min和2 300 r/min;開度分別為0.5 mm,0.6 mm,0.7 mm,0.8 mm,0.9 mm和1.0 mm;背壓分別為0.15 MPa,0.2 MPa,0.25 MPa,0.3 MPa,0.35 MPa和0.4 MPa.

    圖9為排氣門開度0.7 mm時,不同發(fā)動機轉速在不同排氣背壓下的制動功率.從圖9可以看出:背壓越大,發(fā)動機制動性能越好.背壓增大,壓縮氣體越困難,活塞上行阻力越大,由活塞運動產(chǎn)生的制動功率越大.然而排氣背壓達到一定程度時,會引起排氣門的反跳,造成排氣門和排氣機構的損壞,它限制了平均背壓的數(shù)值.

    排氣背壓p/MPa

    圖9 不同轉速在不同背壓下的活塞制動功率

    Fig.9Pistons braking power of different speed 

    under different back pressure

    

    圖10為排氣門開度0.7 mm時,不同排氣背壓在不同發(fā)動機轉速下的活塞制動功率.從圖10可以看出:當背壓與開度一定時,轉速越高,制動性能越好.轉速增加時,氣缸壓力峰值大,壓縮過程充量對活塞做功多,并且單位時間內的工作循環(huán)次數(shù)增多,因此制動功率越大.

    發(fā)動機轉速n/(r?min-1)

    圖10 不同背壓在不同轉速下的活塞制動功率

    Fig.10 Pistons braking power of different back 

    pressure under different speed

    

    圖11為不同轉速在不同排氣門開度下的制動功率.從圖11可以看出:轉速一定時,制動功率隨開度的增加先增加后減少.

    排氣門開度/mm

    圖11 不同轉速在不同開度下的活塞制動功率

    Fig.11 Pistons braking power of different speed 

    under different opening

    

    圖12為排氣背壓為0.3 MPa時,不同開度在不同轉速下的制動功率.從圖12可以看出:對于每一個轉速,存在一個對應的開度使得制動功率最佳.通過數(shù)值模擬的方法可以找到各個轉速下的最佳開度.對于本文模擬的發(fā)動機:背壓為0.3 MPa、轉速為1 800 r/min時,最佳開度為0.6 mm;轉速為1 900 r/min和2 000 r/min時,最佳開度為0.7 mm;轉速為2 100 r/min時,最佳開度為0.8 mm;轉速為2 200 r/min和2 300 r/min時最佳開度為0.9 mm.分析趨勢可知:轉速高,對應的最佳開度越大.轉速高,氣缸壓力峰值大,但排氣時間短,留在氣缸的充量多,膨脹過程對活塞做功多;轉速低則反之.開度大,排氣徹底,但峰值低,消耗功??;開度小則反之.

    發(fā)動機轉速n/(r?min-1)

    圖12 不同開度在不同轉速下的活塞制動功率

    Fig.12 Pistons braking power of different opening 

    under different speed

    5 泄氣制動性能優(yōu)化

    對于本文研究的發(fā)動機,排氣門開度為1.0 mm時,壓縮過程功率消耗較少,制動功率較??;當開度為0.5 mm時,活塞最大受力值增加,增加了壓縮過程中的制動功率,在膨脹沖程中,缸內氣體對活塞做功也增加了;當開度為0 mm時,發(fā)動機處于倒拖的極限狀態(tài),此時發(fā)動機的制動功率很小,因此,對于每一個轉速,都存在一個對應的最佳開度使得制動功率最大.表1為不同轉速在最大、最小及對應的最佳開度下的制動功率.從表1可以看出,在各轉速下開度從最大值或最小值逐漸調為最佳值時,不同程度的增加了制動功率.

    表1 不同轉速在最大、最小、最佳開度時的制動功率

    Tab.1 Braking power of different opening under maximum, 

    minimum, best openingkW 

    開度/mm轉速/(r?min-1)

    1 800

    1 900

    2 000

    2 100

    2 200

    2 300

    0.5

    201.0

    210.6

    222.0

    231

    238.8

    244.8

    1.0

    174.6

    194.4

    212.4

    231.5

    252.0

    270.6

    最佳開度

    205.6

    217.2

    228.6

    242.4

    258.6

    277.2

    6 實驗驗證

    圖13為實驗裝置示意圖.發(fā)動機制動能力試驗中,采用電力測功機拖動發(fā)動機,在各個狀態(tài)參數(shù)下運轉,測量不同工況下的拖動功率,并通過壓電式壓力傳感器將缸內壓力等非電量信號轉化為電量信號,通過信號放大器將電量信號放大,同時通過儀器采集曲軸上止點的脈沖信號,將這些信號經(jīng)過A/D轉換器,將模擬信號轉化為數(shù)字信號,送入計算機進行數(shù)據(jù)處理并輸出數(shù)據(jù).制動臺架實驗如圖14所示,試驗發(fā)動機的主要參數(shù)如表2所示.

    圖13 發(fā)動機制動過程測量系統(tǒng)

    Fig.13 Measurement system of engine 

    bleeding brake process

    圖14發(fā)動機泄氣制動臺架

    Fig.14 Bench of engine bleeding brake

    

    表2 發(fā)動機主要參數(shù)

    Tab.2 Main parameters of the engine

    項目

    參數(shù)

    布置形式

    直列、6缸、水冷、四沖程

    氣缸直徑/mm

    126

    活塞行程/mm

    166

    排量/L

    12.4

    進氣方式

    雙級渦輪增壓

    壓縮比

    17

    氣門數(shù)

    4

    額定功率轉速/(r?min-1)

    2 100

    最大輸出功率/kW

    350

    圖15為轉速2 100 r/min、背壓0.3 MPa、開度0.7 mm時活塞受力的實驗值與模擬值的比較.從圖15可以看出:雖然模擬結果和實驗結果之間存在一定誤差,但是兩者的變化趨勢基本一致,說明在選取計算模型和設定邊界條件時基本反映了泄氣制動系統(tǒng)的實際工作情況.

    曲軸轉角n/(°)

    圖15 活塞壓力實驗值與模擬值的比較

    Fig.15 Piston pressure comparison between

    simulated value and experimental value

    

    表3為排氣背壓為0.3 MPa時,各轉速對應的最佳排氣門開度的模擬值與實驗值的比較.從表3可以看出:除轉速為2 000 r/min外,其他轉速下模擬找到的最佳開度值與實驗找到的最佳開度值相同.

    表3 各轉速下最佳開度模擬值與實驗值對比

    Tab.3 Best opening comparison between simulated value 

    and experimental value under different speed

    轉速

    /(r?min-1)

    最佳開度模擬值

    /mm

    最佳開度實驗值

    /mm

    1 800

    0.6

    0.6

    1 900

    0.7

    0.7

    2 000

    0.7

    0.8

    2 100

    0.8

    0.8

    2 200

    0.9

    0.9

    2 300

    0.9

    0.9

    圖16為背壓為0.3 MPa、開度為0.7 mm時不同轉速下單缸制動功率的模擬值與實驗測得的指示功率的比較.從圖16可以看出:數(shù)值模擬結果與實驗結果誤差十分微小.最大誤差為5.2%,最小誤差僅為4.1%,說明瞬態(tài)數(shù)值模擬具有較高的精確性.

    發(fā)動機轉速n/(r?min-1)

    圖16制動功率模擬值與實驗值的比較

    Fig.16Braking power comparison between simulated 

    value and experimental value

    7 結 論

    對泄氣制動過程進行三維數(shù)值模擬計算,比傳統(tǒng)的數(shù)值計算更加精確,分析結果表明:

    1)活塞受力受發(fā)動機轉速、排氣門開度及排氣背壓的影響.轉速越高,活塞受力峰值越大;背壓越大,活塞受力峰值越大;開度越大,活塞受力峰值越小.活塞受力峰值出現(xiàn)在壓縮上止點前,并且隨轉速的升高、背壓的增大、開度的降低越靠近上止點.

    2)轉速越高,泄氣制動性能越好;背壓越大,制動性能越好,然而背壓太大,會引起排氣門反跳,它限制了平均排氣背壓的數(shù)值;對于每一個轉速,都有一個對應的最佳排氣門開度,使得制動性能最佳.并且轉速越高,對應的最佳排氣門開度越大.由本文的方法可以找到最佳排氣門開度.

    3)對于本文研究的發(fā)動機,當排氣門開度優(yōu)化到最佳開度時,不同程度地提高了泄氣制動性能.實驗結果顯示,實驗值與模擬值誤差很小,表明模擬計算方法合理,結果可信.

    參考文獻

    [1] GB7258-2012 機動車運行安全技術條件[S].北京:中國標 準出版社,2012.

    GB7258-2012 Safety technical condition of vehicle operation[S].Beijing:China Standards Press,2012.(In Chinese)

    [2] 王勇,張雄,劉翰東.汽車輔助制動系統(tǒng)及試驗方法研究[J].客車技術與研究,2013(4):30-33.

    WANG Yong, ZHANG Xiong, LIU Handong. Study on automotive auxiliary braking system and test method[J].Bus & Coach Technology and Research,2013(4):30-33. (In Chinese)

    [3] 夏基勝.發(fā)動機輔助制動試驗研究[J].車用發(fā)動機,2010,8(4):70-73.

    XIA Jisheng. Test study on assistant braking performance of engin[J]. Vehicle Engine,2010,8(4):70-73. (In Chinese)

    [4] 厲江,孔德星,梁文芝,等.發(fā)動機緩速器整車性能計算及試驗驗證[J].合肥學院學報,2012,22(3):75-78.

    LI Jiang, KONG Dexing, LIANG Wenzhi,et al. The vehicle performances calculations and experimentation of engine retarder[J].Journal of Hefei University,2012,22(3):75-78. (In Chinese)

    [5] 周科.客車輔助制動性能試驗研究[J].客車技術與研究,2013(4):41-49.

    ZHOU Ke.Test study on assistant braking performanceof bus[J].Bus & Coach Technology and Research,2013(4):41-49. (In Chinese)

    [6] SONG Bin,LV Jiangang,ZHAO Guangjun,et al.Numerical simulation on gasliquid twophase flow in fluid couplings during braking[J].IEEE ICCSIT,2010(4):709-713.

    [7] 董穎.發(fā)動機制動的工作機理及性能分析研究[D].鎮(zhèn)江:江蘇大學汽車與交通工程學院,2007.

    DONG Ying.Working process and performance research of engine retarder[D].Zhenjiang:School of Automobile and Traffic Engineering, Jiangsu University,2007. (In Chinese)

    [8] DRUZHININAM, STEFANOPULU A G. Speed control experiment for commercial heavy vehicles with coordinated friction and engine compression brakes[J].Proceeding of the American Conference, 2002(3):2546-2551.

    [9] KURNIAWAN Wendy Hardyono, ABDULLAH Shahrir.Numerical analysis of the combustion process in a fourstroke compressed natural gas engine with direct injection system[J].Journal of Mechanical Science and Technology,2008,22(10):1937-1945.

    [10]TANNER Franz X, SRINIVASAN Seshasai. CFDbased optimization of fuel injection strategies in a diesel engine using an adaptive gradient method[J].Applied Mathematical Modelling,2009,33(3):1366-1385.

    [11]杭瑾.發(fā)動機輔助制動三維模擬研究[D].鎮(zhèn)江:江蘇大學汽車與交通工程學院,2010.

    HANG Jin. Study of engine auxiliary braking through 3d numerical simulation[D]. Zhenjiang: School of Automobile and Traffic Engineering, Jiangsu University,2010. (In Chinese)

    [12]于勇.FLUENT入門與進階教程[M].北京:北京理工大學出版社,2008.

    YU Yong.Tutorial of FLUENT introduction and advanced[M].Beijing:Beijing Institute of Technology Press,2008. (In Chinese)

    表3 各轉速下最佳開度模擬值與實驗值對比

    Tab.3 Best opening comparison between simulated value 

    and experimental value under different speed

    轉速

    /(r?min-1)

    最佳開度模擬值

    /mm

    最佳開度實驗值

    /mm

    1 800

    0.6

    0.6

    1 900

    0.7

    0.7

    2 000

    0.7

    0.8

    2 100

    0.8

    0.8

    2 200

    0.9

    0.9

    2 300

    0.9

    0.9

    圖16為背壓為0.3 MPa、開度為0.7 mm時不同轉速下單缸制動功率的模擬值與實驗測得的指示功率的比較.從圖16可以看出:數(shù)值模擬結果與實驗結果誤差十分微小.最大誤差為5.2%,最小誤差僅為4.1%,說明瞬態(tài)數(shù)值模擬具有較高的精確性.

    發(fā)動機轉速n/(r?min-1)

    圖16制動功率模擬值與實驗值的比較

    Fig.16Braking power comparison between simulated 

    value and experimental value

    7 結 論

    對泄氣制動過程進行三維數(shù)值模擬計算,比傳統(tǒng)的數(shù)值計算更加精確,分析結果表明:

    1)活塞受力受發(fā)動機轉速、排氣門開度及排氣背壓的影響.轉速越高,活塞受力峰值越大;背壓越大,活塞受力峰值越大;開度越大,活塞受力峰值越小.活塞受力峰值出現(xiàn)在壓縮上止點前,并且隨轉速的升高、背壓的增大、開度的降低越靠近上止點.

    2)轉速越高,泄氣制動性能越好;背壓越大,制動性能越好,然而背壓太大,會引起排氣門反跳,它限制了平均排氣背壓的數(shù)值;對于每一個轉速,都有一個對應的最佳排氣門開度,使得制動性能最佳.并且轉速越高,對應的最佳排氣門開度越大.由本文的方法可以找到最佳排氣門開度.

    3)對于本文研究的發(fā)動機,當排氣門開度優(yōu)化到最佳開度時,不同程度地提高了泄氣制動性能.實驗結果顯示,實驗值與模擬值誤差很小,表明模擬計算方法合理,結果可信.

    參考文獻

    [1] GB7258-2012 機動車運行安全技術條件[S].北京:中國標 準出版社,2012.

    GB7258-2012 Safety technical condition of vehicle operation[S].Beijing:China Standards Press,2012.(In Chinese)

    [2] 王勇,張雄,劉翰東.汽車輔助制動系統(tǒng)及試驗方法研究[J].客車技術與研究,2013(4):30-33.

    WANG Yong, ZHANG Xiong, LIU Handong. Study on automotive auxiliary braking system and test method[J].Bus & Coach Technology and Research,2013(4):30-33. (In Chinese)

    [3] 夏基勝.發(fā)動機輔助制動試驗研究[J].車用發(fā)動機,2010,8(4):70-73.

    XIA Jisheng. Test study on assistant braking performance of engin[J]. Vehicle Engine,2010,8(4):70-73. (In Chinese)

    [4] 厲江,孔德星,梁文芝,等.發(fā)動機緩速器整車性能計算及試驗驗證[J].合肥學院學報,2012,22(3):75-78.

    LI Jiang, KONG Dexing, LIANG Wenzhi,et al. The vehicle performances calculations and experimentation of engine retarder[J].Journal of Hefei University,2012,22(3):75-78. (In Chinese)

    [5] 周科.客車輔助制動性能試驗研究[J].客車技術與研究,2013(4):41-49.

    ZHOU Ke.Test study on assistant braking performanceof bus[J].Bus & Coach Technology and Research,2013(4):41-49. (In Chinese)

    [6] SONG Bin,LV Jiangang,ZHAO Guangjun,et al.Numerical simulation on gasliquid twophase flow in fluid couplings during braking[J].IEEE ICCSIT,2010(4):709-713.

    [7] 董穎.發(fā)動機制動的工作機理及性能分析研究[D].鎮(zhèn)江:江蘇大學汽車與交通工程學院,2007.

    DONG Ying.Working process and performance research of engine retarder[D].Zhenjiang:School of Automobile and Traffic Engineering, Jiangsu University,2007. (In Chinese)

    [8] DRUZHININAM, STEFANOPULU A G. Speed control experiment for commercial heavy vehicles with coordinated friction and engine compression brakes[J].Proceeding of the American Conference, 2002(3):2546-2551.

    [9] KURNIAWAN Wendy Hardyono, ABDULLAH Shahrir.Numerical analysis of the combustion process in a fourstroke compressed natural gas engine with direct injection system[J].Journal of Mechanical Science and Technology,2008,22(10):1937-1945.

    [10]TANNER Franz X, SRINIVASAN Seshasai. CFDbased optimization of fuel injection strategies in a diesel engine using an adaptive gradient method[J].Applied Mathematical Modelling,2009,33(3):1366-1385.

    [11]杭瑾.發(fā)動機輔助制動三維模擬研究[D].鎮(zhèn)江:江蘇大學汽車與交通工程學院,2010.

    HANG Jin. Study of engine auxiliary braking through 3d numerical simulation[D]. Zhenjiang: School of Automobile and Traffic Engineering, Jiangsu University,2010. (In Chinese)

    [12]于勇.FLUENT入門與進階教程[M].北京:北京理工大學出版社,2008.

    YU Yong.Tutorial of FLUENT introduction and advanced[M].Beijing:Beijing Institute of Technology Press,2008. (In Chinese)

    表3 各轉速下最佳開度模擬值與實驗值對比

    Tab.3 Best opening comparison between simulated value 

    and experimental value under different speed

    轉速

    /(r?min-1)

    最佳開度模擬值

    /mm

    最佳開度實驗值

    /mm

    1 800

    0.6

    0.6

    1 900

    0.7

    0.7

    2 000

    0.7

    0.8

    2 100

    0.8

    0.8

    2 200

    0.9

    0.9

    2 300

    0.9

    0.9

    圖16為背壓為0.3 MPa、開度為0.7 mm時不同轉速下單缸制動功率的模擬值與實驗測得的指示功率的比較.從圖16可以看出:數(shù)值模擬結果與實驗結果誤差十分微小.最大誤差為5.2%,最小誤差僅為4.1%,說明瞬態(tài)數(shù)值模擬具有較高的精確性.

    發(fā)動機轉速n/(r?min-1)

    圖16制動功率模擬值與實驗值的比較

    Fig.16Braking power comparison between simulated 

    value and experimental value

    7 結 論

    對泄氣制動過程進行三維數(shù)值模擬計算,比傳統(tǒng)的數(shù)值計算更加精確,分析結果表明:

    1)活塞受力受發(fā)動機轉速、排氣門開度及排氣背壓的影響.轉速越高,活塞受力峰值越大;背壓越大,活塞受力峰值越大;開度越大,活塞受力峰值越小.活塞受力峰值出現(xiàn)在壓縮上止點前,并且隨轉速的升高、背壓的增大、開度的降低越靠近上止點.

    2)轉速越高,泄氣制動性能越好;背壓越大,制動性能越好,然而背壓太大,會引起排氣門反跳,它限制了平均排氣背壓的數(shù)值;對于每一個轉速,都有一個對應的最佳排氣門開度,使得制動性能最佳.并且轉速越高,對應的最佳排氣門開度越大.由本文的方法可以找到最佳排氣門開度.

    3)對于本文研究的發(fā)動機,當排氣門開度優(yōu)化到最佳開度時,不同程度地提高了泄氣制動性能.實驗結果顯示,實驗值與模擬值誤差很小,表明模擬計算方法合理,結果可信.

    參考文獻

    [1] GB7258-2012 機動車運行安全技術條件[S].北京:中國標 準出版社,2012.

    GB7258-2012 Safety technical condition of vehicle operation[S].Beijing:China Standards Press,2012.(In Chinese)

    [2] 王勇,張雄,劉翰東.汽車輔助制動系統(tǒng)及試驗方法研究[J].客車技術與研究,2013(4):30-33.

    WANG Yong, ZHANG Xiong, LIU Handong. Study on automotive auxiliary braking system and test method[J].Bus & Coach Technology and Research,2013(4):30-33. (In Chinese)

    [3] 夏基勝.發(fā)動機輔助制動試驗研究[J].車用發(fā)動機,2010,8(4):70-73.

    XIA Jisheng. Test study on assistant braking performance of engin[J]. Vehicle Engine,2010,8(4):70-73. (In Chinese)

    [4] 厲江,孔德星,梁文芝,等.發(fā)動機緩速器整車性能計算及試驗驗證[J].合肥學院學報,2012,22(3):75-78.

    LI Jiang, KONG Dexing, LIANG Wenzhi,et al. The vehicle performances calculations and experimentation of engine retarder[J].Journal of Hefei University,2012,22(3):75-78. (In Chinese)

    [5] 周科.客車輔助制動性能試驗研究[J].客車技術與研究,2013(4):41-49.

    ZHOU Ke.Test study on assistant braking performanceof bus[J].Bus & Coach Technology and Research,2013(4):41-49. (In Chinese)

    [6] SONG Bin,LV Jiangang,ZHAO Guangjun,et al.Numerical simulation on gasliquid twophase flow in fluid couplings during braking[J].IEEE ICCSIT,2010(4):709-713.

    [7] 董穎.發(fā)動機制動的工作機理及性能分析研究[D].鎮(zhèn)江:江蘇大學汽車與交通工程學院,2007.

    DONG Ying.Working process and performance research of engine retarder[D].Zhenjiang:School of Automobile and Traffic Engineering, Jiangsu University,2007. (In Chinese)

    [8] DRUZHININAM, STEFANOPULU A G. Speed control experiment for commercial heavy vehicles with coordinated friction and engine compression brakes[J].Proceeding of the American Conference, 2002(3):2546-2551.

    [9] KURNIAWAN Wendy Hardyono, ABDULLAH Shahrir.Numerical analysis of the combustion process in a fourstroke compressed natural gas engine with direct injection system[J].Journal of Mechanical Science and Technology,2008,22(10):1937-1945.

    [10]TANNER Franz X, SRINIVASAN Seshasai. CFDbased optimization of fuel injection strategies in a diesel engine using an adaptive gradient method[J].Applied Mathematical Modelling,2009,33(3):1366-1385.

    [11]杭瑾.發(fā)動機輔助制動三維模擬研究[D].鎮(zhèn)江:江蘇大學汽車與交通工程學院,2010.

    HANG Jin. Study of engine auxiliary braking through 3d numerical simulation[D]. Zhenjiang: School of Automobile and Traffic Engineering, Jiangsu University,2010. (In Chinese)

    [12]于勇.FLUENT入門與進階教程[M].北京:北京理工大學出版社,2008.

    YU Yong.Tutorial of FLUENT introduction and advanced[M].Beijing:Beijing Institute of Technology Press,2008. (In Chinese)

    猜你喜歡
    數(shù)值模擬發(fā)動機
    2015款寶馬525Li行駛中發(fā)動機熄火
    2012年奔馳S600發(fā)動機故障燈偶爾點亮
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    一種基于液壓緩沖的減震管卡設計與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    新一代MTU2000發(fā)動機系列
    發(fā)動機的怠速停止技術i-stop
    亚洲专区字幕在线| 成人手机av| 一级片'在线观看视频| 国产免费视频播放在线视频| 建设人人有责人人尽责人人享有的| 国产精品麻豆人妻色哟哟久久| 亚洲,欧美精品.| 在线看a的网站| 欧美日韩一级在线毛片| 精品高清国产在线一区| 嫩草影视91久久| 国产一区二区三区综合在线观看| 午夜精品久久久久久毛片777| 大片电影免费在线观看免费| 成年版毛片免费区| 久久亚洲精品不卡| 欧美在线一区亚洲| 一级片'在线观看视频| 高清欧美精品videossex| 精品亚洲乱码少妇综合久久| 亚洲av日韩精品久久久久久密| 制服诱惑二区| 丁香六月欧美| 国产欧美日韩一区二区三区在线| 欧美日韩中文字幕国产精品一区二区三区 | 黄色片一级片一级黄色片| 高潮久久久久久久久久久不卡| 亚洲专区国产一区二区| 久久av网站| 亚洲av国产av综合av卡| 久久久久久久久免费视频了| 婷婷成人精品国产| 成人免费观看视频高清| 国产区一区二久久| 亚洲中文日韩欧美视频| 国产日韩欧美亚洲二区| 满18在线观看网站| 黄色丝袜av网址大全| 老司机在亚洲福利影院| 亚洲视频免费观看视频| 桃红色精品国产亚洲av| 脱女人内裤的视频| 自线自在国产av| 成人免费观看视频高清| 亚洲免费av在线视频| 精品少妇内射三级| 亚洲成人免费av在线播放| 免费高清在线观看日韩| 成人国产av品久久久| 色94色欧美一区二区| 日本一区二区免费在线视频| 在线观看舔阴道视频| 乱人伦中国视频| 国产精品一区二区在线观看99| 国产成人av教育| 热re99久久国产66热| 亚洲专区中文字幕在线| 国内毛片毛片毛片毛片毛片| 午夜激情久久久久久久| 亚洲精品成人av观看孕妇| 日韩 欧美 亚洲 中文字幕| 亚洲av日韩在线播放| 欧美一级毛片孕妇| 变态另类成人亚洲欧美熟女 | 黑人欧美特级aaaaaa片| 午夜视频精品福利| 色94色欧美一区二区| 菩萨蛮人人尽说江南好唐韦庄| 精品少妇内射三级| 久久久国产欧美日韩av| 老鸭窝网址在线观看| 国产精品一区二区在线不卡| 999精品在线视频| 久久久水蜜桃国产精品网| 如日韩欧美国产精品一区二区三区| 久久热在线av| 国产高清视频在线播放一区| 亚洲精品美女久久av网站| 久9热在线精品视频| 午夜免费鲁丝| 久久久久久久大尺度免费视频| 精品少妇黑人巨大在线播放| 欧美人与性动交α欧美精品济南到| 国产99久久九九免费精品| 悠悠久久av| 国产亚洲欧美在线一区二区| 欧美成狂野欧美在线观看| 脱女人内裤的视频| 亚洲国产毛片av蜜桃av| 久久亚洲真实| 高潮久久久久久久久久久不卡| 国产高清视频在线播放一区| 亚洲av第一区精品v没综合| 777米奇影视久久| 亚洲午夜精品一区,二区,三区| 午夜福利欧美成人| 亚洲熟女精品中文字幕| 午夜福利影视在线免费观看| 他把我摸到了高潮在线观看 | 亚洲国产成人一精品久久久| 国产av精品麻豆| 一本色道久久久久久精品综合| 成人永久免费在线观看视频 | 成人18禁高潮啪啪吃奶动态图| 一二三四在线观看免费中文在| 国产不卡av网站在线观看| 免费人妻精品一区二区三区视频| 亚洲av片天天在线观看| 国产主播在线观看一区二区| 精品国产乱码久久久久久小说| 男男h啪啪无遮挡| 女人高潮潮喷娇喘18禁视频| 黄频高清免费视频| 国产精品欧美亚洲77777| 亚洲五月色婷婷综合| 久久av网站| 国产欧美日韩综合在线一区二区| 久久久久网色| 欧美日本中文国产一区发布| 中文字幕色久视频| 在线观看免费视频日本深夜| 日韩熟女老妇一区二区性免费视频| 国产成人免费观看mmmm| 国产亚洲av高清不卡| 国产精品二区激情视频| 国产单亲对白刺激| 精品一区二区三区视频在线观看免费 | 欧美乱码精品一区二区三区| 亚洲视频免费观看视频| 一进一出抽搐动态| 黄色视频在线播放观看不卡| 国产欧美亚洲国产| 两人在一起打扑克的视频| 热99国产精品久久久久久7| 亚洲精品久久成人aⅴ小说| 嫁个100分男人电影在线观看| 国产精品电影一区二区三区 | 少妇裸体淫交视频免费看高清 | 午夜福利欧美成人| 国产精品一区二区在线观看99| 一区二区av电影网| 欧美激情久久久久久爽电影 | 精品卡一卡二卡四卡免费| tube8黄色片| 一级,二级,三级黄色视频| 亚洲欧美日韩高清在线视频 | 国产男女内射视频| 999久久久国产精品视频| 99国产综合亚洲精品| 少妇的丰满在线观看| 亚洲欧美精品综合一区二区三区| 美女高潮到喷水免费观看| 桃花免费在线播放| 精品一区二区三区视频在线观看免费 | 老司机亚洲免费影院| 久久久精品免费免费高清| 亚洲天堂av无毛| 9191精品国产免费久久| 国产不卡av网站在线观看| 在线观看免费午夜福利视频| 国产精品久久久久久精品电影小说| 国产欧美日韩一区二区三| 免费在线观看黄色视频的| 国产亚洲精品一区二区www | 欧美人与性动交α欧美精品济南到| 在线观看舔阴道视频| 人妻 亚洲 视频| 欧美日韩视频精品一区| 欧美国产精品va在线观看不卡| 国产精品免费一区二区三区在线 | 一级片'在线观看视频| 天堂俺去俺来也www色官网| 久久久国产精品麻豆| 婷婷丁香在线五月| 天天躁夜夜躁狠狠躁躁| 狠狠狠狠99中文字幕| 搡老乐熟女国产| 黑人巨大精品欧美一区二区蜜桃| 男女边摸边吃奶| 日日爽夜夜爽网站| 久久毛片免费看一区二区三区| 久9热在线精品视频| 不卡av一区二区三区| 国产精品免费一区二区三区在线 | 黑人猛操日本美女一级片| 老司机靠b影院| 性色av乱码一区二区三区2| 亚洲精品国产色婷婷电影| 大型黄色视频在线免费观看| 国产成人影院久久av| 热re99久久国产66热| 国产主播在线观看一区二区| av片东京热男人的天堂| 侵犯人妻中文字幕一二三四区| 亚洲精品国产精品久久久不卡| 欧美成人免费av一区二区三区 | 国产野战对白在线观看| 亚洲色图 男人天堂 中文字幕| 又黄又粗又硬又大视频| 亚洲精品国产精品久久久不卡| 老司机福利观看| 另类精品久久| 国产伦人伦偷精品视频| 一区二区三区国产精品乱码| 日本精品一区二区三区蜜桃| 大片电影免费在线观看免费| 中文字幕制服av| 高清黄色对白视频在线免费看| 久久99热这里只频精品6学生| 在线观看一区二区三区激情| 人人妻人人澡人人看| 国产精品98久久久久久宅男小说| 一个人免费在线观看的高清视频| 高潮久久久久久久久久久不卡| 亚洲欧美色中文字幕在线| 两性夫妻黄色片| 超碰成人久久| 日韩熟女老妇一区二区性免费视频| 国产精品免费大片| 精品人妻1区二区| 这个男人来自地球电影免费观看| 久久精品亚洲熟妇少妇任你| 国产成人精品在线电影| 高清在线国产一区| 丁香六月天网| 一进一出抽搐动态| 蜜桃国产av成人99| 中国美女看黄片| 国产成人精品久久二区二区免费| www.999成人在线观看| 中文字幕色久视频| 在线亚洲精品国产二区图片欧美| 欧美精品人与动牲交sv欧美| 欧美人与性动交α欧美软件| 亚洲av日韩精品久久久久久密| 黑人猛操日本美女一级片| 精品亚洲乱码少妇综合久久| 亚洲熟女精品中文字幕| 国产精品98久久久久久宅男小说| 51午夜福利影视在线观看| 啦啦啦中文免费视频观看日本| 一区福利在线观看| av又黄又爽大尺度在线免费看| 成人av一区二区三区在线看| 一区二区三区精品91| 又黄又粗又硬又大视频| 久久国产精品大桥未久av| 搡老乐熟女国产| 亚洲欧美色中文字幕在线| av福利片在线| 蜜桃国产av成人99| 在线观看免费日韩欧美大片| 在线观看66精品国产| 极品教师在线免费播放| 久久青草综合色| av欧美777| 菩萨蛮人人尽说江南好唐韦庄| 精品国产超薄肉色丝袜足j| 一夜夜www| 一级毛片精品| 91精品三级在线观看| 国产激情久久老熟女| 乱人伦中国视频| 制服诱惑二区| 精品国产亚洲在线| 最近最新免费中文字幕在线| 757午夜福利合集在线观看| 成人亚洲精品一区在线观看| 国产人伦9x9x在线观看| 国产淫语在线视频| 精品久久蜜臀av无| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩一级在线毛片| 黄色 视频免费看| 亚洲精华国产精华精| 夜夜骑夜夜射夜夜干| 欧美激情极品国产一区二区三区| 黄色怎么调成土黄色| 日韩精品免费视频一区二区三区| 色婷婷av一区二区三区视频| 好男人电影高清在线观看| 超色免费av| 国产精品国产高清国产av | e午夜精品久久久久久久| 蜜桃国产av成人99| 如日韩欧美国产精品一区二区三区| 成年版毛片免费区| cao死你这个sao货| 国产精品熟女久久久久浪| videosex国产| 亚洲情色 制服丝袜| 男女无遮挡免费网站观看| 国产极品粉嫩免费观看在线| 一区二区三区激情视频| 嫁个100分男人电影在线观看| 国产精品国产高清国产av | 久久久久国产一级毛片高清牌| 自线自在国产av| 香蕉国产在线看| 亚洲专区国产一区二区| 亚洲精华国产精华精| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久久毛片微露脸| 久久青草综合色| 国产91精品成人一区二区三区 | 亚洲精品美女久久av网站| 欧美大码av| 久久精品aⅴ一区二区三区四区| 狂野欧美激情性xxxx| 国产深夜福利视频在线观看| 免费在线观看黄色视频的| 妹子高潮喷水视频| 这个男人来自地球电影免费观看| 精品少妇一区二区三区视频日本电影| 亚洲国产欧美一区二区综合| 97人妻天天添夜夜摸| 免费观看av网站的网址| 欧美乱码精品一区二区三区| 亚洲av电影在线进入| 啦啦啦中文免费视频观看日本| 久久av网站| svipshipincom国产片| 在线观看人妻少妇| 国产一区二区在线观看av| 精品国产超薄肉色丝袜足j| 久久精品成人免费网站| 国产精品免费视频内射| 久久久久国内视频| 欧美国产精品va在线观看不卡| 亚洲九九香蕉| 女人爽到高潮嗷嗷叫在线视频| 国产在线免费精品| 脱女人内裤的视频| 91精品三级在线观看| 久久热在线av| 91精品三级在线观看| 亚洲欧洲精品一区二区精品久久久| 久久精品熟女亚洲av麻豆精品| 日本vs欧美在线观看视频| 中文字幕制服av| 麻豆国产av国片精品| 激情视频va一区二区三区| 丝袜人妻中文字幕| 中文亚洲av片在线观看爽 | 黑人猛操日本美女一级片| 国产真人三级小视频在线观看| 亚洲成a人片在线一区二区| 丝瓜视频免费看黄片| 麻豆国产av国片精品| 激情视频va一区二区三区| 又大又爽又粗| 精品一区二区三区视频在线观看免费 | 老司机靠b影院| 国产色视频综合| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 丰满饥渴人妻一区二区三| 一边摸一边抽搐一进一小说 | 国产91精品成人一区二区三区 | 国产免费现黄频在线看| 啦啦啦在线免费观看视频4| 欧美中文综合在线视频| 中国美女看黄片| 十八禁网站免费在线| 在线看a的网站| 日韩欧美免费精品| 精品国产乱码久久久久久小说| 多毛熟女@视频| 国产精品免费视频内射| 国产不卡一卡二| 最近最新免费中文字幕在线| 亚洲欧洲精品一区二区精品久久久| 精品福利观看| 精品人妻在线不人妻| cao死你这个sao货| 国产精品欧美亚洲77777| 国产老妇伦熟女老妇高清| 欧美亚洲日本最大视频资源| 国产精品久久久久久精品古装| 国产午夜精品久久久久久| 久久久精品国产亚洲av高清涩受| 色播在线永久视频| 丝瓜视频免费看黄片| 超色免费av| 丰满人妻熟妇乱又伦精品不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 老熟女久久久| 国产野战对白在线观看| 在线永久观看黄色视频| 国产欧美日韩一区二区三| 国产精品久久久久久人妻精品电影 | 首页视频小说图片口味搜索| 亚洲中文日韩欧美视频| 伊人久久大香线蕉亚洲五| 欧美激情 高清一区二区三区| 母亲3免费完整高清在线观看| 首页视频小说图片口味搜索| 久久精品国产综合久久久| 欧美亚洲日本最大视频资源| 久久久欧美国产精品| 国产精品久久久久久人妻精品电影 | 亚洲avbb在线观看| 老司机亚洲免费影院| 国产男靠女视频免费网站| 黑丝袜美女国产一区| 无人区码免费观看不卡 | 国产精品久久久人人做人人爽| 满18在线观看网站| 大片电影免费在线观看免费| 两性午夜刺激爽爽歪歪视频在线观看 | 大片电影免费在线观看免费| 99国产精品一区二区三区| 精品亚洲成a人片在线观看| 五月开心婷婷网| 午夜精品久久久久久毛片777| 国产成人av激情在线播放| 天堂俺去俺来也www色官网| 18禁国产床啪视频网站| 亚洲国产欧美一区二区综合| av网站在线播放免费| 亚洲自偷自拍图片 自拍| 在线观看人妻少妇| 久久中文字幕一级| 亚洲第一欧美日韩一区二区三区 | 精品一区二区三区视频在线观看免费 | 一区二区三区激情视频| 另类亚洲欧美激情| 亚洲精品在线观看二区| 欧美日韩视频精品一区| 亚洲综合色网址| 777久久人妻少妇嫩草av网站| 巨乳人妻的诱惑在线观看| 黄频高清免费视频| 久久久久久亚洲精品国产蜜桃av| 亚洲av片天天在线观看| 大型黄色视频在线免费观看| 久久精品人人爽人人爽视色| 亚洲国产av新网站| 亚洲午夜理论影院| videos熟女内射| 亚洲精品中文字幕一二三四区 | 亚洲精品成人av观看孕妇| 国产主播在线观看一区二区| 亚洲一区二区三区欧美精品| 精品一区二区三区av网在线观看 | 天天影视国产精品| 国产精品久久久久久人妻精品电影 | 黄片播放在线免费| 中文字幕制服av| 精品久久久久久久毛片微露脸| 一进一出抽搐动态| 高潮久久久久久久久久久不卡| 动漫黄色视频在线观看| 一本一本久久a久久精品综合妖精| 久久久久国内视频| 国产精品欧美亚洲77777| 国产精品一区二区精品视频观看| 男女无遮挡免费网站观看| 中文字幕人妻丝袜制服| 99久久精品国产亚洲精品| 无限看片的www在线观看| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 日韩成人在线观看一区二区三区| 夫妻午夜视频| 亚洲专区中文字幕在线| 国产精品久久久久久人妻精品电影 | 熟女少妇亚洲综合色aaa.| 中文字幕制服av| 1024视频免费在线观看| 黑人操中国人逼视频| 亚洲第一青青草原| 自拍欧美九色日韩亚洲蝌蚪91| 动漫黄色视频在线观看| 亚洲一码二码三码区别大吗| 人人妻人人添人人爽欧美一区卜| av欧美777| 人人妻人人澡人人看| 国产亚洲av高清不卡| 18禁美女被吸乳视频| 一区二区三区国产精品乱码| 日韩欧美一区二区三区在线观看 | 亚洲av日韩精品久久久久久密| 成人18禁在线播放| 精品亚洲成国产av| 国产色视频综合| 午夜激情av网站| 亚洲中文av在线| 亚洲情色 制服丝袜| av视频免费观看在线观看| 久9热在线精品视频| 午夜激情av网站| 久久人妻福利社区极品人妻图片| 亚洲情色 制服丝袜| h视频一区二区三区| 99热网站在线观看| 老司机靠b影院| 嫁个100分男人电影在线观看| 亚洲一区中文字幕在线| 露出奶头的视频| 伦理电影免费视频| 在线观看免费视频网站a站| 色婷婷av一区二区三区视频| 在线观看www视频免费| 国产精品久久久久久人妻精品电影 | 999精品在线视频| 日本av免费视频播放| 亚洲欧美日韩另类电影网站| 久久中文字幕一级| 王馨瑶露胸无遮挡在线观看| 成年人免费黄色播放视频| 国产精品一区二区在线不卡| 精品一区二区三卡| 最新在线观看一区二区三区| 国产深夜福利视频在线观看| 亚洲熟女精品中文字幕| 女同久久另类99精品国产91| 亚洲精品中文字幕一二三四区 | 岛国毛片在线播放| 亚洲情色 制服丝袜| 亚洲视频免费观看视频| 国产亚洲午夜精品一区二区久久| 国产日韩欧美在线精品| 亚洲国产欧美在线一区| 丁香欧美五月| 少妇的丰满在线观看| 无遮挡黄片免费观看| 久久久精品国产亚洲av高清涩受| 自线自在国产av| 国产成人av教育| 欧美+亚洲+日韩+国产| 午夜老司机福利片| 乱人伦中国视频| 亚洲国产av影院在线观看| 欧美+亚洲+日韩+国产| 日本wwww免费看| 久久人人97超碰香蕉20202| 亚洲人成电影免费在线| 久久久久久久大尺度免费视频| 久久久国产成人免费| 亚洲国产欧美日韩在线播放| 一区二区三区精品91| 不卡一级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 妹子高潮喷水视频| 久久九九热精品免费| 飞空精品影院首页| 国产熟女午夜一区二区三区| 婷婷成人精品国产| 欧美人与性动交α欧美精品济南到| 999精品在线视频| 99久久国产精品久久久| 久久人人97超碰香蕉20202| 国产成人一区二区三区免费视频网站| 国产在线一区二区三区精| 菩萨蛮人人尽说江南好唐韦庄| av线在线观看网站| 操美女的视频在线观看| 精品亚洲成国产av| 91麻豆av在线| 两个人免费观看高清视频| 搡老乐熟女国产| 欧美黑人欧美精品刺激| 久久久久久久久免费视频了| 国产三级黄色录像| 一级毛片精品| 视频区欧美日本亚洲| 亚洲专区字幕在线| 精品福利永久在线观看| 在线观看66精品国产| 大香蕉久久成人网| 欧美亚洲日本最大视频资源| av福利片在线| 十八禁高潮呻吟视频| 一本色道久久久久久精品综合| 国产一区有黄有色的免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 少妇被粗大的猛进出69影院| 亚洲精品国产色婷婷电影| 午夜福利,免费看| 免费在线观看日本一区| 色94色欧美一区二区| 天天躁日日躁夜夜躁夜夜| 狠狠婷婷综合久久久久久88av| 成年版毛片免费区| 最新的欧美精品一区二区| 超碰97精品在线观看| 热99re8久久精品国产| 亚洲欧美色中文字幕在线| 国产成人啪精品午夜网站| 1024香蕉在线观看| 国产1区2区3区精品| 午夜福利在线免费观看网站| 99国产精品一区二区三区| 日韩人妻精品一区2区三区| 精品国产超薄肉色丝袜足j| 国产黄频视频在线观看| 少妇粗大呻吟视频| 成人国产av品久久久| 啪啪无遮挡十八禁网站| 女同久久另类99精品国产91| 69精品国产乱码久久久| 极品少妇高潮喷水抽搐| 大码成人一级视频| 热99久久久久精品小说推荐| 黄色a级毛片大全视频| 久久精品人人爽人人爽视色| 两性午夜刺激爽爽歪歪视频在线观看 | 国产色视频综合| 999久久久国产精品视频| 日韩一卡2卡3卡4卡2021年| 在线观看一区二区三区激情| 久久久精品94久久精品| 91精品三级在线观看| 99九九在线精品视频|