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

    試驗與仿真相結(jié)合的發(fā)動機活塞熱負荷分析

    2014-09-18 13:46王毅楊靖鄧幫林張云飛孫承劉凱敏
    關(guān)鍵詞:溫度場活塞

    王毅+楊靖+鄧幫林+張云飛+孫承+劉凱敏

    收稿日期:20131025

    基金項目:國家高技術(shù)研究發(fā)展計劃(863計劃)資助項目(2008AA11A114);湖南省研究生科研創(chuàng)新基金資助項目(CX20108125)

    作者簡介:王毅(1986-),男,湖北利川人,湖南大學(xué)博士研究生

    通訊聯(lián)系人,Email:yangjing10@vip.com

    摘要:借助硬度塞溫度測試、數(shù)值分析手段,結(jié)合對改進前活塞的溫度場和熱應(yīng)力的計算,對發(fā)動機改進后活塞的熱負荷狀況進行評估.通過計算活塞傳熱邊界條件,利用有限元分析軟件,對活塞進行溫度場計算,然后利用試驗實測數(shù)據(jù)對其進行標(biāo)定,在計算誤差小于5%的情況下,得到改進前與改進后活塞在A工況和B工況下的溫度場分布,最后進行了對應(yīng)工況的熱應(yīng)力計算.結(jié)果表明在相同工況下發(fā)動機改進后活塞最高溫度均比改進前活塞低,最大熱應(yīng)力比改進前活塞小.在改進前活塞滿足熱負荷要求的前提下,改進后活塞滿足熱負荷要求.

    關(guān)鍵詞:活塞;硬度塞法;溫度場;熱應(yīng)力

    中圖分類號:U464.134 文獻標(biāo)識碼:A

    Thermal Load Analysis of Engine Piston Combined

    Experiment with Simulation

    WANG Yi1,2,YANG Jing1,2,DENG Banglin1,2,ZHANG Yunfei1,2,SUN Cheng1,2,LIU Kaimin1,2

    (1.State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan Univ, Changsha,Hunan410082, China;

    2.Research Center for Advanced Powertrain Technology, Hunan Univ, Changsha, Hunan 410082, China)

    Abstract: In order to investigate the thermal load distribution of modified piston, it is necessary to analyze temperature field and thermal stress combined simulation with experiments. The FE software was used to simulate the temperature field of the piston after the boundary conditions of thermal transfer were calculated, and then, the results were verified according to experiment data, and the error of calculation was less than 5%. Finally, temperature field and thermal stress were derived under different operating conditions including conditions A and B. The results of investigation show that the maximum temperature and maximum value of thermal stress of modified piston decrease compared with the original piston. The modified piston meets the demand of thermal load based on the original piston meeting the demand of thermal load.

    Key words:pistons;plug of hardness method;temperature field;thermal stress

    活塞是發(fā)動機的主要受熱零件,工作時,處于高溫、高壓、高負荷的惡劣環(huán)境,經(jīng)受周期性交變機械負荷和熱負荷,容易發(fā)生故障[1].熱負荷是造成活塞失效的主要因素之一,隨著發(fā)動機強化指標(biāo)的不斷提高,它的影響作用也越來越突出[2-5].為準(zhǔn)確計算活塞在各個工況下的應(yīng)力大小,全面考察活塞的強度,對活塞進行疲勞分析、敲缸分析以及潤滑磨損分析,確定活塞的溫度場分布是進行上述分析的前提[6].因此,活塞溫度場計算的精度顯得尤為重要.目前,活塞溫度場測試試驗主要是硬度塞法,硬度塞安裝方便,同時能夠保證一定的精度.謝琰等人研究了發(fā)動機活塞的溫度場[4],但研究工況比較單一.夏飛等人利用有限元方法研究了活塞溫度場及熱應(yīng)力[5],但缺乏試驗的驗證,其研究結(jié)果缺乏可信度.由于發(fā)動機生產(chǎn)工藝發(fā)生變化,需要對發(fā)動機活塞進行再設(shè)計以滿足現(xiàn)有的工藝要求.活塞熱負荷是否滿足要求是活塞再設(shè)計成功與否的重要評價指標(biāo).本文通過發(fā)動機多次試驗表明,改進前活塞滿足熱負荷要求.計算活塞傳熱邊界條件,利用ABAQUS有限元分析軟件,對活塞進行溫度場計算,然后利用試驗實測數(shù)據(jù)對其進行標(biāo)定,在計算誤差小于5%的狀況下,得到改進前與改進后活塞在A工況和B工況下的溫度場分布,最后進行了對應(yīng)工況的熱應(yīng)力計算,評估改進后活塞熱負荷是否滿足工程應(yīng)用要求.

    1活塞溫度測試

    1.1硬度塞回火曲線的標(biāo)定

    硬度塞法測試溫度的原理是利用某種合金經(jīng)過淬火后產(chǎn)生永久性硬度變化,在不同的溫度下進行回火,其表面硬度也將隨之變化的現(xiàn)象.對硬度塞材料而言,其回火溫度和硬度變化最好呈直線關(guān)系或近似線性關(guān)系,同時要求材料在同一回火溫度下硬度值穩(wěn)定.因此,本次試驗中選用材料為GCr15滾珠軸承鋼.

    硬度塞的淬火處理是在氣體保護電爐中進行,硬度測量采用維氏硬度計.淬火后抽取10只硬度塞進行硬度測量,在每只硬度塞上取3個測量點.測量數(shù)據(jù)表明,淬火后的硬度塞硬度均勻性很好,硬度值偏差在±5 HV內(nèi),滿足硬度塞材料要求,如圖1所示.

    硬度塞達到所需硬度要求后,進行硬度塞回火試驗,整理回火數(shù)據(jù),并利用最小二乘法擬合得到硬度溫度曲線,如圖2所示,誤差分析見表1.

    1.2活塞溫度測點布置

    活塞主要幾何參數(shù)如表2所示.在滿足活塞潤滑性能的前提下,改進前后活塞結(jié)構(gòu)對比如圖3所示,改進后活塞油環(huán)槽高度減小,第3氣環(huán)槽高度加大,頂部厚度增加.為了準(zhǔn)確預(yù)測活塞溫度場分布,本次試驗選用20個測點.根據(jù)活塞結(jié)構(gòu),布置硬度塞測點,改進前后活塞測點位置相同,如圖4所示.

    1.3發(fā)動機試驗

    本次試驗分4組進行,即改進前后活塞分別在最大扭矩工況(工況A)和最大功率工況(工況B)下進行試驗,試驗條件見表3.試驗時最大限度地滿足外界條件的一致性,保證發(fā)動機前后試驗均在相同狀態(tài)下運行.

    1.4活塞測點溫度及數(shù)據(jù)分析

    改進前后活塞在兩種不同工況下的測點溫度分別如圖5和圖6所示.從圖中可以看出,在A工況下,改進后活塞大部分溫度測點比改進前相對應(yīng)測點低,B工況亦如此.同時,A工況下活塞測點溫度較之B工況下測點溫度波動大,這與發(fā)動機運行狀況相關(guān),溫度波動大表明A工況運行工況不穩(wěn)定.

    2.2活塞材料屬性

    活塞材料為鍛鋁2A80,其密度為2.77 g/cm3,彈性模量為72 GPa,泊松比為0.33,抗拉強度為375 MPa;其他隨溫度變化屬性如表5所示.

    2.3活塞傳熱邊界條件的確定

    2.3.1活塞頂部傳熱邊界條件

    發(fā)動機工作過程中,燃料燃燒釋放的熱能轉(zhuǎn)化為機械能,該過程中活塞頂部直接與高溫燃氣接觸,熱量通過活塞頂面?zhèn)鞯交钊^部,然后通過冷卻油腔和活塞環(huán)將熱量傳給其他冷卻介質(zhì)[7].根據(jù)周期瞬態(tài)溫度波動理論[5, 8],活塞頂面的溫度沿活塞頂法線方向迅速衰減,而這個溫度的波動只發(fā)生在活塞頂面1~2 mm的表層,不對活塞的溫度場產(chǎn)生較大的影響,所以在某一特定工況下,一般將活塞溫度場近似為穩(wěn)定的溫度場.

    活塞頂面的傳熱邊界條件包括傳熱系數(shù)和燃氣溫度.此次試驗測試了發(fā)動機的氣缸壓力,已知燃氣壓力及氣缸容積后,根據(jù)式(1)計算氣缸瞬時溫度,B工況氣缸壓力曲線如圖8所示.

    Tg=PgVmR.(1)

    瞬時放熱系數(shù)采用艾歇伯格(Eickelberg)公式[9]計算,即

    αg=7.83umPgTg.(2)

    式中:um為活塞平均速度;Pg為氣體瞬時壓力(MPa);Tg為氣體瞬時溫度(K).

    計算氣缸內(nèi)瞬時壓力和溫度后,根據(jù)式(3)和式(4)計算每個循環(huán)的平均對流放熱系數(shù)和平均溫度[6].

    αgm=∫7200 αg720dθ,(3)

    Tgm=∫7200αgTg720dθαgm. (4)

    氣缸工作容積隨曲軸轉(zhuǎn)角的變化關(guān)系為[8]:

    V=

    Vs22εc-1+1-cosφ+1λs1-1-λ2ssin2φ. (5)

    式中:Vs為氣缸總?cè)莘e;εc為發(fā)動機壓縮比;λs為發(fā)動機曲柄連桿比.

    為了更精確地分析活塞頂面?zhèn)鳠徇吔?,對活塞頂部進行分區(qū).根據(jù)活塞試驗所得式(6),(7)來進行活塞頂面?zhèn)鳠嵯禂?shù)分區(qū)[6],根據(jù)本文研究發(fā)動機結(jié)構(gòu)參數(shù),計算中N=38.8 mm是式(6)和式(7)的分界點.

    αr=2αm1+e0.1N1.5e0.1r25.41.5,0

    αr=2αm1+e0.1N1.5e0.12N-r25.41.5,N

    2.3.2活塞內(nèi)腔傳熱邊界條件

    活塞內(nèi)腔換熱系數(shù)跟活塞內(nèi)腔溫度與曲軸箱內(nèi)部氣流的流動狀況相關(guān).活塞內(nèi)腔表面的換熱系數(shù)為:

    αoil=(T1-T2)kT1-Toilσ.(8)

    式中:T1,T2和Toil分別為活塞頂部溫度、活塞內(nèi)腔底部溫度和曲軸箱氣體溫度;k,σ分別為活塞的導(dǎo)熱系數(shù)、活塞頂厚度;αoil為活塞內(nèi)腔表面與曲軸箱氣體的換熱系數(shù).

    2.3.3活塞側(cè)面?zhèn)鳠徇吔鐥l件

    活塞火力岸、環(huán)槽和裙部換熱系數(shù)比較難確定,目前一般采用經(jīng)驗公式來確定.傳熱模型如圖9所示,影響這些區(qū)域的傳熱系數(shù)的因素主要有:氣膜、油膜、活塞環(huán)、缸套的厚度,還有油膜、活塞環(huán)、缸套的換熱系數(shù)和缸套與冷卻水之間的換熱系數(shù),式(9)~(14)為活塞側(cè)面區(qū)域的換熱系數(shù)求解公式.

    火力岸換熱系數(shù)為:

    1αk=n1λ3+bλ2+1αw. (9)

    第1氣環(huán)槽下緣換熱系數(shù)為:

    1αk=0.95n1λ3+bλ2+1αw. (10)

    第2氣環(huán)槽下緣換熱系數(shù)為:

    1αk=dλ3+bλ2+n2λ2+m22λ1+m22λ0. (11)

    油環(huán)上下緣表面換熱系數(shù)為:

    1αk=dλ3+bλ2+c0λ0+n0λ0+1αw. (12)

    環(huán)槽內(nèi)側(cè)表面換熱系數(shù)為:

    1αk=dλ3+bλ2+l0λ0+n0λ0+1αw. (13)

    活塞裙部表面換熱系數(shù)為:

    1αk=Δ2λ0+bλ2+1αw. (14)

    式中:λ1為燃氣的導(dǎo)熱系數(shù);λ2為缸套的導(dǎo)熱系數(shù);λ3為活塞環(huán)的導(dǎo)熱系數(shù);λ0為機油的導(dǎo)熱系數(shù);αw為缸套和水腔之間的換熱系數(shù);Δ2為活塞裙部與缸套的間隙;b為缸套厚度;l0為油環(huán)厚度;d為環(huán)傳熱中心間距;n1,n2分別為油環(huán)與第3氣環(huán)的環(huán)外間隙.活塞兩個工況下各個區(qū)域的對流換熱邊界條件見表6.

    3計算結(jié)果及分析

    3.1溫度場計算結(jié)果及分析

    不同工況下改進前與改進后活塞與試驗對比如圖10~圖13所示,對比結(jié)果表明計算值與試驗值吻合較好,誤差均小于5%,有效驗證了此次計算的精度達到工程應(yīng)用要求.

    改進前與改進后活塞在2種不同工況下溫度場計算結(jié)果如圖14~圖17所示.從結(jié)果可以看出,兩種狀態(tài)下,改進前最高溫度為298.8 ℃,改進后活塞最高溫度為278.5 ℃,改進后最高溫度比改進前有所降低,改進后活塞頂部區(qū)域溫度比改進前低.改進后油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前活塞略高.這主要是因為改進后活塞油槽結(jié)構(gòu)所引起的,油環(huán)槽高度降低,活塞頭部熱量更好的向活塞裙部傳遞,活塞裙部的熱量則通過第三氣環(huán)槽的活塞環(huán)散熱,減少了熱量的累積,同時減小了溫度梯度.

    活塞內(nèi)腔頂部路徑溫度分布如圖18和圖19所示,改進后活塞第3氣環(huán)槽加寬使得活塞散熱更好,因此,A狀態(tài)下,改進前活塞溫度梯度比改進后大,表明沿此路徑,改進前活塞熱應(yīng)力比改進后活塞大.B狀態(tài)下所表現(xiàn)出來的此種現(xiàn)象更加明顯.

    3.2活塞熱應(yīng)力分析

    根據(jù)溫度場分布,結(jié)合有限元分析軟件,計算出活塞的熱應(yīng)力,其結(jié)果對評價活塞熱強度具有一定的指導(dǎo)意義.活塞在不同工況下熱應(yīng)力計算結(jié)果如圖20~圖23所示,結(jié)果表明各處的熱應(yīng)力均在50 MPa以下,滿足材料強度要求;其最大熱應(yīng)力主要集中在活塞內(nèi)腔頂部、活塞油環(huán)槽下邊緣、活塞減重孔邊緣處.改進后的活塞熱應(yīng)力均有明顯下降,尤其是活塞頂部的內(nèi)表面處.

    4結(jié)論

    1)在相同工況下,改進后活塞最高溫度比改進前最高溫度低,改進后活塞頂部區(qū)域溫度比改進前低,油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前略高.

    2)在狀態(tài)A和狀態(tài)B下,改進后活塞熱應(yīng)力最大值比改進前活塞熱應(yīng)力均有減小.

    3)通過試驗與仿真相結(jié)合分析活塞溫度場及熱應(yīng)力分布狀況是解決活塞熱負荷問題的有效途徑,結(jié)構(gòu)改進后溫度分布和熱應(yīng)力均得到有效改善.

    參考文獻

    [1]楊進華.微型HCCI自由活塞發(fā)動機工作過程實驗及數(shù)值模擬研究[D].鎮(zhèn)江:江蘇大學(xué)能源與動力工程學(xué)院,2012.

    YANG Jinhua.Study on experiments and numerical simulation of working process of micro HCCI freepiston engine[D].Zhenjiang: School of Energy and Power Engineering,Jiangsu University, 2012. (In Chinese)

    [2]CERIT M,AYHAN V,PARLAK A.Thermal analysis of a partially ceramic coated piston: effect on cold start HC emission in a spark ignition engine[J]. Applied Thermal Engineering,2010, 31(2011):336-341.

    [3]FORMOSA F.Coupled thermodynamicdynamic semianalytical model of free piston stirling engines[J].Energy Conversion and Management, 2010,52(2011):2098-2109.

    [4]謝琰,席明智,劉曉麗.柴油機縮口四角ω燃燒室活塞溫度場試驗研究[J].內(nèi)燃機,2010,2(1):38-42.

    XIE Yan, XI Mingzhi, LIU Xiaoli. Experimental research to necking ω combustion chamberpiston temperature field[J].Internal Combustion Engine,2010,2(1):38-42. (In Chinese)

    [5]夏飛,蘇鐵熊,馮耀楠,等.基于質(zhì)量限制的全鋼活塞溫度場的三維數(shù)值模擬[J].拖拉機與農(nóng)用運輸車,2010,37(1):42-45.

    XIA Fei,SU Tiexiong, FENG Yaonan,et al.3D numerical simulation for temperature field of steel piston with quality restrictions[J].Tractor & Farm Transporter,2010,37(1):42-45.(In Chinese)

    [6]張桂昌.基于熱機耦合的柴油機活塞系統(tǒng)敲擊噪聲與潤滑研究及優(yōu)化設(shè)計[D].天津:天津大學(xué)機械工程學(xué)院,2012.

    ZHANG Guichang.Research on impact noise and lubrication of diesel engine piston assembly considering thermalmechanical coupling loads and optimization design[D].Tianjin:School of Mechanical Engineering,Tianjing University, 2012. (In Chinese)

    [7]仲杰.活塞噴油振蕩冷卻的穩(wěn)、瞬態(tài)模擬計算及活塞溫度場分析[D]. 濟南:山東大學(xué)能源與動力工程學(xué)院,2012.

    ZHONG Jie.Steady and transient numerical simulation of oilinject shaking cooling and temperature distribution in piston[D].Jinan:School of Energy and Power Engineering, Shangdong University, 2012. (In Chinese)

    [8]周龍保.內(nèi)燃機學(xué)[M].北京:機械工業(yè)出版社,2005:43-44.

    ZHOU Longbao.Fundamentals of engine[M]. Beijing: Machinery Industry Press, 2005:43-44.(In Chinese)

    [9]王國強.船舶柴油機活塞和缸套溫度場有限元分析[D].大連:大連海事大學(xué)輪機工程學(xué)院,2012.

    WANG Guoqiang.FEA for the thermal field of the piston and cylinder liner in marine diesel engine[D].Dalian:College of Marine, Dalian Maritime University, 2012. (In Chinese)

    3計算結(jié)果及分析

    3.1溫度場計算結(jié)果及分析

    不同工況下改進前與改進后活塞與試驗對比如圖10~圖13所示,對比結(jié)果表明計算值與試驗值吻合較好,誤差均小于5%,有效驗證了此次計算的精度達到工程應(yīng)用要求.

    改進前與改進后活塞在2種不同工況下溫度場計算結(jié)果如圖14~圖17所示.從結(jié)果可以看出,兩種狀態(tài)下,改進前最高溫度為298.8 ℃,改進后活塞最高溫度為278.5 ℃,改進后最高溫度比改進前有所降低,改進后活塞頂部區(qū)域溫度比改進前低.改進后油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前活塞略高.這主要是因為改進后活塞油槽結(jié)構(gòu)所引起的,油環(huán)槽高度降低,活塞頭部熱量更好的向活塞裙部傳遞,活塞裙部的熱量則通過第三氣環(huán)槽的活塞環(huán)散熱,減少了熱量的累積,同時減小了溫度梯度.

    活塞內(nèi)腔頂部路徑溫度分布如圖18和圖19所示,改進后活塞第3氣環(huán)槽加寬使得活塞散熱更好,因此,A狀態(tài)下,改進前活塞溫度梯度比改進后大,表明沿此路徑,改進前活塞熱應(yīng)力比改進后活塞大.B狀態(tài)下所表現(xiàn)出來的此種現(xiàn)象更加明顯.

    3.2活塞熱應(yīng)力分析

    根據(jù)溫度場分布,結(jié)合有限元分析軟件,計算出活塞的熱應(yīng)力,其結(jié)果對評價活塞熱強度具有一定的指導(dǎo)意義.活塞在不同工況下熱應(yīng)力計算結(jié)果如圖20~圖23所示,結(jié)果表明各處的熱應(yīng)力均在50 MPa以下,滿足材料強度要求;其最大熱應(yīng)力主要集中在活塞內(nèi)腔頂部、活塞油環(huán)槽下邊緣、活塞減重孔邊緣處.改進后的活塞熱應(yīng)力均有明顯下降,尤其是活塞頂部的內(nèi)表面處.

    4結(jié)論

    1)在相同工況下,改進后活塞最高溫度比改進前最高溫度低,改進后活塞頂部區(qū)域溫度比改進前低,油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前略高.

    2)在狀態(tài)A和狀態(tài)B下,改進后活塞熱應(yīng)力最大值比改進前活塞熱應(yīng)力均有減小.

    3)通過試驗與仿真相結(jié)合分析活塞溫度場及熱應(yīng)力分布狀況是解決活塞熱負荷問題的有效途徑,結(jié)構(gòu)改進后溫度分布和熱應(yīng)力均得到有效改善.

    參考文獻

    [1]楊進華.微型HCCI自由活塞發(fā)動機工作過程實驗及數(shù)值模擬研究[D].鎮(zhèn)江:江蘇大學(xué)能源與動力工程學(xué)院,2012.

    YANG Jinhua.Study on experiments and numerical simulation of working process of micro HCCI freepiston engine[D].Zhenjiang: School of Energy and Power Engineering,Jiangsu University, 2012. (In Chinese)

    [2]CERIT M,AYHAN V,PARLAK A.Thermal analysis of a partially ceramic coated piston: effect on cold start HC emission in a spark ignition engine[J]. Applied Thermal Engineering,2010, 31(2011):336-341.

    [3]FORMOSA F.Coupled thermodynamicdynamic semianalytical model of free piston stirling engines[J].Energy Conversion and Management, 2010,52(2011):2098-2109.

    [4]謝琰,席明智,劉曉麗.柴油機縮口四角ω燃燒室活塞溫度場試驗研究[J].內(nèi)燃機,2010,2(1):38-42.

    XIE Yan, XI Mingzhi, LIU Xiaoli. Experimental research to necking ω combustion chamberpiston temperature field[J].Internal Combustion Engine,2010,2(1):38-42. (In Chinese)

    [5]夏飛,蘇鐵熊,馮耀楠,等.基于質(zhì)量限制的全鋼活塞溫度場的三維數(shù)值模擬[J].拖拉機與農(nóng)用運輸車,2010,37(1):42-45.

    XIA Fei,SU Tiexiong, FENG Yaonan,et al.3D numerical simulation for temperature field of steel piston with quality restrictions[J].Tractor & Farm Transporter,2010,37(1):42-45.(In Chinese)

    [6]張桂昌.基于熱機耦合的柴油機活塞系統(tǒng)敲擊噪聲與潤滑研究及優(yōu)化設(shè)計[D].天津:天津大學(xué)機械工程學(xué)院,2012.

    ZHANG Guichang.Research on impact noise and lubrication of diesel engine piston assembly considering thermalmechanical coupling loads and optimization design[D].Tianjin:School of Mechanical Engineering,Tianjing University, 2012. (In Chinese)

    [7]仲杰.活塞噴油振蕩冷卻的穩(wěn)、瞬態(tài)模擬計算及活塞溫度場分析[D]. 濟南:山東大學(xué)能源與動力工程學(xué)院,2012.

    ZHONG Jie.Steady and transient numerical simulation of oilinject shaking cooling and temperature distribution in piston[D].Jinan:School of Energy and Power Engineering, Shangdong University, 2012. (In Chinese)

    [8]周龍保.內(nèi)燃機學(xué)[M].北京:機械工業(yè)出版社,2005:43-44.

    ZHOU Longbao.Fundamentals of engine[M]. Beijing: Machinery Industry Press, 2005:43-44.(In Chinese)

    [9]王國強.船舶柴油機活塞和缸套溫度場有限元分析[D].大連:大連海事大學(xué)輪機工程學(xué)院,2012.

    WANG Guoqiang.FEA for the thermal field of the piston and cylinder liner in marine diesel engine[D].Dalian:College of Marine, Dalian Maritime University, 2012. (In Chinese)

    3計算結(jié)果及分析

    3.1溫度場計算結(jié)果及分析

    不同工況下改進前與改進后活塞與試驗對比如圖10~圖13所示,對比結(jié)果表明計算值與試驗值吻合較好,誤差均小于5%,有效驗證了此次計算的精度達到工程應(yīng)用要求.

    改進前與改進后活塞在2種不同工況下溫度場計算結(jié)果如圖14~圖17所示.從結(jié)果可以看出,兩種狀態(tài)下,改進前最高溫度為298.8 ℃,改進后活塞最高溫度為278.5 ℃,改進后最高溫度比改進前有所降低,改進后活塞頂部區(qū)域溫度比改進前低.改進后油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前活塞略高.這主要是因為改進后活塞油槽結(jié)構(gòu)所引起的,油環(huán)槽高度降低,活塞頭部熱量更好的向活塞裙部傳遞,活塞裙部的熱量則通過第三氣環(huán)槽的活塞環(huán)散熱,減少了熱量的累積,同時減小了溫度梯度.

    活塞內(nèi)腔頂部路徑溫度分布如圖18和圖19所示,改進后活塞第3氣環(huán)槽加寬使得活塞散熱更好,因此,A狀態(tài)下,改進前活塞溫度梯度比改進后大,表明沿此路徑,改進前活塞熱應(yīng)力比改進后活塞大.B狀態(tài)下所表現(xiàn)出來的此種現(xiàn)象更加明顯.

    3.2活塞熱應(yīng)力分析

    根據(jù)溫度場分布,結(jié)合有限元分析軟件,計算出活塞的熱應(yīng)力,其結(jié)果對評價活塞熱強度具有一定的指導(dǎo)意義.活塞在不同工況下熱應(yīng)力計算結(jié)果如圖20~圖23所示,結(jié)果表明各處的熱應(yīng)力均在50 MPa以下,滿足材料強度要求;其最大熱應(yīng)力主要集中在活塞內(nèi)腔頂部、活塞油環(huán)槽下邊緣、活塞減重孔邊緣處.改進后的活塞熱應(yīng)力均有明顯下降,尤其是活塞頂部的內(nèi)表面處.

    4結(jié)論

    1)在相同工況下,改進后活塞最高溫度比改進前最高溫度低,改進后活塞頂部區(qū)域溫度比改進前低,油環(huán)槽溫度與活塞裙部區(qū)域溫度改進后比改進前略高.

    2)在狀態(tài)A和狀態(tài)B下,改進后活塞熱應(yīng)力最大值比改進前活塞熱應(yīng)力均有減小.

    3)通過試驗與仿真相結(jié)合分析活塞溫度場及熱應(yīng)力分布狀況是解決活塞熱負荷問題的有效途徑,結(jié)構(gòu)改進后溫度分布和熱應(yīng)力均得到有效改善.

    參考文獻

    [1]楊進華.微型HCCI自由活塞發(fā)動機工作過程實驗及數(shù)值模擬研究[D].鎮(zhèn)江:江蘇大學(xué)能源與動力工程學(xué)院,2012.

    YANG Jinhua.Study on experiments and numerical simulation of working process of micro HCCI freepiston engine[D].Zhenjiang: School of Energy and Power Engineering,Jiangsu University, 2012. (In Chinese)

    [2]CERIT M,AYHAN V,PARLAK A.Thermal analysis of a partially ceramic coated piston: effect on cold start HC emission in a spark ignition engine[J]. Applied Thermal Engineering,2010, 31(2011):336-341.

    [3]FORMOSA F.Coupled thermodynamicdynamic semianalytical model of free piston stirling engines[J].Energy Conversion and Management, 2010,52(2011):2098-2109.

    [4]謝琰,席明智,劉曉麗.柴油機縮口四角ω燃燒室活塞溫度場試驗研究[J].內(nèi)燃機,2010,2(1):38-42.

    XIE Yan, XI Mingzhi, LIU Xiaoli. Experimental research to necking ω combustion chamberpiston temperature field[J].Internal Combustion Engine,2010,2(1):38-42. (In Chinese)

    [5]夏飛,蘇鐵熊,馮耀楠,等.基于質(zhì)量限制的全鋼活塞溫度場的三維數(shù)值模擬[J].拖拉機與農(nóng)用運輸車,2010,37(1):42-45.

    XIA Fei,SU Tiexiong, FENG Yaonan,et al.3D numerical simulation for temperature field of steel piston with quality restrictions[J].Tractor & Farm Transporter,2010,37(1):42-45.(In Chinese)

    [6]張桂昌.基于熱機耦合的柴油機活塞系統(tǒng)敲擊噪聲與潤滑研究及優(yōu)化設(shè)計[D].天津:天津大學(xué)機械工程學(xué)院,2012.

    ZHANG Guichang.Research on impact noise and lubrication of diesel engine piston assembly considering thermalmechanical coupling loads and optimization design[D].Tianjin:School of Mechanical Engineering,Tianjing University, 2012. (In Chinese)

    [7]仲杰.活塞噴油振蕩冷卻的穩(wěn)、瞬態(tài)模擬計算及活塞溫度場分析[D]. 濟南:山東大學(xué)能源與動力工程學(xué)院,2012.

    ZHONG Jie.Steady and transient numerical simulation of oilinject shaking cooling and temperature distribution in piston[D].Jinan:School of Energy and Power Engineering, Shangdong University, 2012. (In Chinese)

    [8]周龍保.內(nèi)燃機學(xué)[M].北京:機械工業(yè)出版社,2005:43-44.

    ZHOU Longbao.Fundamentals of engine[M]. Beijing: Machinery Industry Press, 2005:43-44.(In Chinese)

    [9]王國強.船舶柴油機活塞和缸套溫度場有限元分析[D].大連:大連海事大學(xué)輪機工程學(xué)院,2012.

    WANG Guoqiang.FEA for the thermal field of the piston and cylinder liner in marine diesel engine[D].Dalian:College of Marine, Dalian Maritime University, 2012. (In Chinese)

    猜你喜歡
    溫度場活塞
    基于并聯(lián)雙離心送風(fēng)機內(nèi)外流場的數(shù)值模擬研究
    能源樁群溫度場分布特征數(shù)值仿真研究
    中班科學(xué)區(qū)游戲活動設(shè)計:好玩的針筒
    淺談發(fā)動機活塞冷卻的專利情況分析
    基于多物理場的T10鋼熱力耦合分析
    基于FloTHERM的抗惡劣環(huán)境計算機熱仿真
    復(fù)合式空調(diào)系統(tǒng)的數(shù)值模擬研究
    活塞常見損傷的檢查與修理
    活塞銷的連接配合與故障分析
    如何識別鋁活塞上的各種字標(biāo)
    看片在线看免费视频| 国产精品久久久久久亚洲av鲁大| 国产免费av片在线观看野外av| 亚洲无线在线观看| 国产真实乱freesex| 在线播放国产精品三级| av天堂中文字幕网| 国产不卡一卡二| 欧美一区二区亚洲| 老司机午夜十八禁免费视频| 9191精品国产免费久久| 久久欧美精品欧美久久欧美| 亚洲性夜色夜夜综合| 久99久视频精品免费| 久久人妻av系列| 亚洲av二区三区四区| 无遮挡黄片免费观看| 成人特级av手机在线观看| 毛片女人毛片| 中文字幕久久专区| 国产爱豆传媒在线观看| 丰满乱子伦码专区| 国产日本99.免费观看| 成人欧美大片| 欧美一区二区精品小视频在线| 国产一区二区在线观看日韩| 色哟哟·www| 女人十人毛片免费观看3o分钟| 嫩草影院入口| 18禁黄网站禁片午夜丰满| 午夜精品一区二区三区免费看| 在线a可以看的网站| 一区福利在线观看| 成人无遮挡网站| 内地一区二区视频在线| 亚洲人成网站在线播放欧美日韩| 九九久久精品国产亚洲av麻豆| 日日干狠狠操夜夜爽| 色视频www国产| 看黄色毛片网站| 91久久精品电影网| 嫩草影院精品99| 欧美xxxx性猛交bbbb| 免费人成在线观看视频色| 成年女人看的毛片在线观看| 日韩欧美一区二区三区在线观看| 欧美黑人巨大hd| 色av中文字幕| 国产白丝娇喘喷水9色精品| av国产免费在线观看| 我要搜黄色片| 特大巨黑吊av在线直播| 三级男女做爰猛烈吃奶摸视频| 亚洲美女视频黄频| 毛片女人毛片| 中文字幕高清在线视频| 久久国产乱子伦精品免费另类| 国产精品精品国产色婷婷| x7x7x7水蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 久久国产乱子伦精品免费另类| 一级毛片久久久久久久久女| 精品日产1卡2卡| 亚洲va日本ⅴa欧美va伊人久久| 久久国产乱子免费精品| 欧美xxxx黑人xx丫x性爽| 亚洲第一电影网av| av在线蜜桃| 亚洲av成人av| 亚洲人成伊人成综合网2020| 欧美激情久久久久久爽电影| 我的女老师完整版在线观看| 国产亚洲av嫩草精品影院| 男女那种视频在线观看| 身体一侧抽搐| 亚洲性夜色夜夜综合| 国产69精品久久久久777片| 欧美高清成人免费视频www| 国产精品久久久久久亚洲av鲁大| 欧美激情在线99| 精品日产1卡2卡| 激情在线观看视频在线高清| 黄色女人牲交| 能在线免费观看的黄片| 欧美日韩瑟瑟在线播放| 免费电影在线观看免费观看| 欧美xxxx性猛交bbbb| 色综合婷婷激情| 99国产极品粉嫩在线观看| 国产一区二区三区在线臀色熟女| 色综合欧美亚洲国产小说| 欧美激情久久久久久爽电影| 亚洲18禁久久av| 亚洲欧美日韩高清专用| 午夜日韩欧美国产| 午夜a级毛片| 性插视频无遮挡在线免费观看| 欧美国产日韩亚洲一区| 成人av一区二区三区在线看| 亚洲熟妇中文字幕五十中出| 久久久久性生活片| 国产精品国产高清国产av| 精品国产三级普通话版| 亚洲欧美日韩高清专用| 欧美+日韩+精品| 一级a爱片免费观看的视频| 久久6这里有精品| 亚洲avbb在线观看| 亚洲中文字幕一区二区三区有码在线看| 深夜精品福利| 久久久久性生活片| 婷婷精品国产亚洲av| 波多野结衣高清无吗| 国内揄拍国产精品人妻在线| 九九久久精品国产亚洲av麻豆| 嫩草影院精品99| 蜜桃亚洲精品一区二区三区| 一个人看的www免费观看视频| 免费av毛片视频| 男人的好看免费观看在线视频| 欧美一级a爱片免费观看看| 成人av一区二区三区在线看| а√天堂www在线а√下载| 成熟少妇高潮喷水视频| 动漫黄色视频在线观看| 国产一区二区亚洲精品在线观看| 国产单亲对白刺激| 在线a可以看的网站| 十八禁网站免费在线| 日韩有码中文字幕| 一个人看视频在线观看www免费| 久久这里只有精品中国| 成年人黄色毛片网站| 国产av麻豆久久久久久久| 久久欧美精品欧美久久欧美| 国产精品久久视频播放| 怎么达到女性高潮| 中文字幕人成人乱码亚洲影| 欧美xxxx黑人xx丫x性爽| 亚洲欧美日韩高清在线视频| 日韩成人在线观看一区二区三区| 高清毛片免费观看视频网站| 免费人成在线观看视频色| 国产av在哪里看| 69人妻影院| 国产 一区 欧美 日韩| 91久久精品电影网| 国产伦人伦偷精品视频| 又黄又爽又刺激的免费视频.| 每晚都被弄得嗷嗷叫到高潮| 亚洲狠狠婷婷综合久久图片| 在线播放国产精品三级| 精品人妻1区二区| 91在线精品国自产拍蜜月| 直男gayav资源| 欧美日韩黄片免| 久久久久久大精品| 成年女人毛片免费观看观看9| 久久国产精品影院| 亚洲最大成人av| 婷婷丁香在线五月| 日韩免费av在线播放| 性色av乱码一区二区三区2| 久久精品国产亚洲av天美| 99视频精品全部免费 在线| 变态另类丝袜制服| 日本三级黄在线观看| 亚洲第一区二区三区不卡| 久久久精品欧美日韩精品| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区激情视频| 国内毛片毛片毛片毛片毛片| 又黄又爽又刺激的免费视频.| 午夜福利在线在线| 99riav亚洲国产免费| 日日摸夜夜添夜夜添小说| 村上凉子中文字幕在线| 亚洲成av人片免费观看| 757午夜福利合集在线观看| 国产精品免费一区二区三区在线| 亚洲一区高清亚洲精品| 亚洲一区二区三区色噜噜| 亚洲av电影在线进入| 亚洲av电影不卡..在线观看| 久久久久久九九精品二区国产| 超碰av人人做人人爽久久| 欧美另类亚洲清纯唯美| 成人av一区二区三区在线看| 亚洲五月天丁香| 国产精品野战在线观看| 国产免费av片在线观看野外av| 欧美日韩福利视频一区二区| 少妇熟女aⅴ在线视频| 精品乱码久久久久久99久播| 亚洲第一欧美日韩一区二区三区| 在线天堂最新版资源| 我要看日韩黄色一级片| 国产色爽女视频免费观看| 亚洲av成人不卡在线观看播放网| 青草久久国产| 久久精品国产亚洲av涩爱 | 美女高潮的动态| 日本成人三级电影网站| 哪里可以看免费的av片| 精品人妻一区二区三区麻豆 | 国产伦精品一区二区三区四那| 丁香欧美五月| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区三区四区免费观看 | 国产精品乱码一区二三区的特点| 少妇的逼好多水| 国产精品亚洲av一区麻豆| 在线观看免费视频日本深夜| 亚洲最大成人av| 级片在线观看| av天堂在线播放| 国产精品久久视频播放| 久久久久免费精品人妻一区二区| 欧美精品国产亚洲| 天堂影院成人在线观看| 国产av在哪里看| 嫁个100分男人电影在线观看| 亚洲成av人片在线播放无| 我的女老师完整版在线观看| 成年免费大片在线观看| 亚洲中文字幕一区二区三区有码在线看| 一个人免费在线观看电影| 国产 一区 欧美 日韩| 亚洲自拍偷在线| 在现免费观看毛片| 欧美绝顶高潮抽搐喷水| 中文字幕熟女人妻在线| 国产综合懂色| 村上凉子中文字幕在线| 非洲黑人性xxxx精品又粗又长| 99视频精品全部免费 在线| 偷拍熟女少妇极品色| 淫妇啪啪啪对白视频| 亚洲国产精品999在线| 在现免费观看毛片| 国产视频一区二区在线看| 亚洲性夜色夜夜综合| 男女那种视频在线观看| www.色视频.com| 亚洲一区高清亚洲精品| 国产精品三级大全| 丰满人妻熟妇乱又伦精品不卡| av天堂在线播放| 亚洲最大成人手机在线| 婷婷精品国产亚洲av| 听说在线观看完整版免费高清| 国产又黄又爽又无遮挡在线| 国产伦精品一区二区三区视频9| 狂野欧美白嫩少妇大欣赏| 九九在线视频观看精品| 男女之事视频高清在线观看| 人人妻人人看人人澡| 高清日韩中文字幕在线| 亚洲人成网站在线播| 又粗又爽又猛毛片免费看| 久久午夜福利片| 搡老熟女国产l中国老女人| 亚洲综合色惰| 1024手机看黄色片| 日日摸夜夜添夜夜添小说| 国产精品国产高清国产av| 亚洲av成人不卡在线观看播放网| 国产一区二区三区在线臀色熟女| 国产精华一区二区三区| 色综合欧美亚洲国产小说| 日韩欧美在线乱码| 男人舔奶头视频| 国产一区二区亚洲精品在线观看| 国产淫片久久久久久久久 | 麻豆成人av在线观看| 国产私拍福利视频在线观看| 日韩欧美精品免费久久 | 男女床上黄色一级片免费看| 日本黄色片子视频| 成年版毛片免费区| 性色av乱码一区二区三区2| 看片在线看免费视频| 国产精品,欧美在线| 久久久国产成人精品二区| 国产乱人视频| 精品一区二区三区视频在线| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 日韩国内少妇激情av| a级毛片a级免费在线| 十八禁国产超污无遮挡网站| 日本 av在线| 桃色一区二区三区在线观看| 97碰自拍视频| 18禁黄网站禁片午夜丰满| 动漫黄色视频在线观看| 一本久久中文字幕| av天堂中文字幕网| 中亚洲国语对白在线视频| 欧美乱色亚洲激情| av在线天堂中文字幕| 免费一级毛片在线播放高清视频| 在线免费观看的www视频| 亚洲美女搞黄在线观看 | 国内少妇人妻偷人精品xxx网站| 中国美女看黄片| 国产精品乱码一区二三区的特点| 午夜福利免费观看在线| 国产精品野战在线观看| 亚洲国产欧美人成| 日日摸夜夜添夜夜添小说| 日本精品一区二区三区蜜桃| 淫妇啪啪啪对白视频| 久久久久久久久大av| 国产精品久久久久久久久免 | 男人和女人高潮做爰伦理| 在线观看66精品国产| 亚洲七黄色美女视频| 赤兔流量卡办理| АⅤ资源中文在线天堂| 午夜精品一区二区三区免费看| 高清日韩中文字幕在线| 老女人水多毛片| 国产蜜桃级精品一区二区三区| 人妻久久中文字幕网| 国产野战对白在线观看| 97超级碰碰碰精品色视频在线观看| 成人av一区二区三区在线看| 能在线免费观看的黄片| 又粗又爽又猛毛片免费看| 国产淫片久久久久久久久 | 亚洲av第一区精品v没综合| 99久久精品一区二区三区| 国产欧美日韩精品亚洲av| 十八禁国产超污无遮挡网站| 日本精品一区二区三区蜜桃| 国产色婷婷99| 九色成人免费人妻av| 精品久久久久久久久久免费视频| 深夜精品福利| 99热精品在线国产| 亚洲国产欧美人成| 色综合婷婷激情| 久久精品综合一区二区三区| 天堂av国产一区二区熟女人妻| 男女床上黄色一级片免费看| 国产探花在线观看一区二区| 国产一区二区亚洲精品在线观看| 亚洲中文日韩欧美视频| 久久精品综合一区二区三区| 国产亚洲精品综合一区在线观看| 噜噜噜噜噜久久久久久91| 国产av麻豆久久久久久久| 99在线人妻在线中文字幕| 中文字幕久久专区| 国产亚洲欧美在线一区二区| 9191精品国产免费久久| 亚洲av一区综合| 久久久久免费精品人妻一区二区| 日本免费a在线| 国产亚洲精品av在线| 亚洲av电影在线进入| 在线国产一区二区在线| av在线天堂中文字幕| 91九色精品人成在线观看| 亚洲性夜色夜夜综合| 日日摸夜夜添夜夜添av毛片 | 99久久精品国产亚洲精品| 此物有八面人人有两片| 国产69精品久久久久777片| 久久人人精品亚洲av| 日本黄色片子视频| 日韩中字成人| 欧美高清性xxxxhd video| 国产主播在线观看一区二区| 色尼玛亚洲综合影院| 色综合亚洲欧美另类图片| 日本在线视频免费播放| 97超视频在线观看视频| 欧美一级a爱片免费观看看| 欧美精品国产亚洲| 波多野结衣巨乳人妻| 久久久久久国产a免费观看| 免费一级毛片在线播放高清视频| 舔av片在线| 好看av亚洲va欧美ⅴa在| 深爱激情五月婷婷| 国产精品国产高清国产av| 夜夜爽天天搞| 欧美日韩综合久久久久久 | 中文字幕精品亚洲无线码一区| 欧美+亚洲+日韩+国产| 免费在线观看成人毛片| 午夜免费成人在线视频| 美女高潮的动态| 亚洲精品日韩av片在线观看| 久久伊人香网站| 久久99热这里只有精品18| 国内久久婷婷六月综合欲色啪| 久久婷婷人人爽人人干人人爱| 久久精品国产亚洲av天美| 久久久久久久久大av| 九色国产91popny在线| 天堂动漫精品| 99久久精品一区二区三区| 淫秽高清视频在线观看| 欧美日韩国产亚洲二区| 69人妻影院| 网址你懂的国产日韩在线| 天天躁日日操中文字幕| 中文亚洲av片在线观看爽| 免费观看的影片在线观看| 久久精品91蜜桃| 久久久国产成人免费| ponron亚洲| 国产成人av教育| 婷婷精品国产亚洲av在线| 两人在一起打扑克的视频| 无遮挡黄片免费观看| 亚洲av不卡在线观看| 少妇丰满av| 精品久久久久久久久久久久久| 欧美午夜高清在线| 欧美不卡视频在线免费观看| 99riav亚洲国产免费| 国产一区二区三区视频了| 久99久视频精品免费| 精品人妻1区二区| 老司机午夜福利在线观看视频| av国产免费在线观看| 久久精品国产亚洲av天美| 国产精品免费一区二区三区在线| 美女被艹到高潮喷水动态| 亚洲狠狠婷婷综合久久图片| 国产毛片a区久久久久| 久久久久久大精品| 日韩欧美国产在线观看| 欧美日本亚洲视频在线播放| 国产欧美日韩一区二区精品| 又紧又爽又黄一区二区| 露出奶头的视频| 成人美女网站在线观看视频| 男人舔女人下体高潮全视频| 非洲黑人性xxxx精品又粗又长| 我的老师免费观看完整版| 在线观看午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲性夜色夜夜综合| 少妇的逼水好多| 国产欧美日韩精品亚洲av| 国产aⅴ精品一区二区三区波| 亚洲国产精品sss在线观看| 一进一出抽搐动态| 亚洲中文字幕日韩| 中亚洲国语对白在线视频| 人妻制服诱惑在线中文字幕| 99国产综合亚洲精品| 国产在线男女| 国产一区二区在线av高清观看| 真实男女啪啪啪动态图| 国产主播在线观看一区二区| 亚洲av不卡在线观看| 亚洲国产精品999在线| 夜夜夜夜夜久久久久| av国产免费在线观看| 中文字幕人妻熟人妻熟丝袜美| 一二三四社区在线视频社区8| 在线播放无遮挡| 日韩av在线大香蕉| 亚洲精品在线美女| x7x7x7水蜜桃| 草草在线视频免费看| 老女人水多毛片| 欧美日韩中文字幕国产精品一区二区三区| 一本综合久久免费| 最后的刺客免费高清国语| av在线蜜桃| 美女高潮的动态| 久久久成人免费电影| 高潮久久久久久久久久久不卡| 乱人视频在线观看| 三级男女做爰猛烈吃奶摸视频| 麻豆成人午夜福利视频| 丁香欧美五月| 日本一本二区三区精品| 欧美性感艳星| 变态另类成人亚洲欧美熟女| 18禁黄网站禁片午夜丰满| 一级a爱片免费观看的视频| 国产大屁股一区二区在线视频| 有码 亚洲区| 国产精品亚洲美女久久久| 午夜福利在线观看吧| 午夜福利高清视频| 亚洲第一欧美日韩一区二区三区| 亚洲av美国av| 亚洲午夜理论影院| 18美女黄网站色大片免费观看| 99久久久亚洲精品蜜臀av| 久久99热6这里只有精品| 精品一区二区免费观看| 中文字幕高清在线视频| 性插视频无遮挡在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 蜜桃亚洲精品一区二区三区| 国产成年人精品一区二区| 搡老熟女国产l中国老女人| 一区二区三区免费毛片| 此物有八面人人有两片| 久久精品国产亚洲av天美| 麻豆国产97在线/欧美| 老司机深夜福利视频在线观看| 亚洲五月婷婷丁香| 亚洲美女黄片视频| 国产高清三级在线| 我要搜黄色片| 永久网站在线| 成人鲁丝片一二三区免费| 色综合亚洲欧美另类图片| 看片在线看免费视频| 免费人成视频x8x8入口观看| 国内少妇人妻偷人精品xxx网站| 国产单亲对白刺激| 天堂动漫精品| 在线十欧美十亚洲十日本专区| 亚洲av免费高清在线观看| 人人妻人人澡欧美一区二区| 国产 一区 欧美 日韩| 欧美性猛交╳xxx乱大交人| 亚洲精品一卡2卡三卡4卡5卡| 国产美女午夜福利| 久久婷婷人人爽人人干人人爱| 午夜视频国产福利| 男女视频在线观看网站免费| 成人国产一区最新在线观看| 色尼玛亚洲综合影院| 九九热线精品视视频播放| 激情在线观看视频在线高清| 国产美女午夜福利| 黄色女人牲交| 国产亚洲精品av在线| 俺也久久电影网| 午夜影院日韩av| 国产高清三级在线| 日本a在线网址| 国产一级毛片七仙女欲春2| 免费看日本二区| 欧美黄色片欧美黄色片| 国产伦在线观看视频一区| www.www免费av| 国产精品影院久久| 最近最新中文字幕大全电影3| 欧美性猛交黑人性爽| 看十八女毛片水多多多| 非洲黑人性xxxx精品又粗又长| 亚洲综合色惰| 亚洲一区二区三区不卡视频| 桃红色精品国产亚洲av| 亚洲美女黄片视频| 人妻丰满熟妇av一区二区三区| av视频在线观看入口| 久久热精品热| 黄色丝袜av网址大全| 18禁在线播放成人免费| 国产一区二区三区视频了| 人妻久久中文字幕网| 久久精品综合一区二区三区| 精品不卡国产一区二区三区| 首页视频小说图片口味搜索| 欧美黑人欧美精品刺激| 欧美绝顶高潮抽搐喷水| 简卡轻食公司| 色综合亚洲欧美另类图片| 亚洲无线观看免费| 青草久久国产| 色噜噜av男人的天堂激情| 赤兔流量卡办理| a在线观看视频网站| 深夜a级毛片| 窝窝影院91人妻| 久久午夜福利片| 国产午夜精品久久久久久一区二区三区 | 久久久久久久久久黄片| 欧美+日韩+精品| 午夜精品久久久久久毛片777| 国产主播在线观看一区二区| 在线观看舔阴道视频| 久久久色成人| 热99在线观看视频| 亚洲无线观看免费| 麻豆成人午夜福利视频| 国产蜜桃级精品一区二区三区| 91狼人影院| 99精品久久久久人妻精品| 亚洲片人在线观看| 亚洲国产日韩欧美精品在线观看| 午夜福利视频1000在线观看| 亚洲成人中文字幕在线播放| 成人精品一区二区免费| 日本一二三区视频观看| 亚洲精品成人久久久久久| 国产av麻豆久久久久久久| 日韩欧美精品免费久久 | 又黄又爽又刺激的免费视频.| 岛国在线免费视频观看| 国产乱人视频| 尤物成人国产欧美一区二区三区| 国产黄色小视频在线观看| 国内少妇人妻偷人精品xxx网站| 一本精品99久久精品77| 一区二区三区四区激情视频 | 国产成人啪精品午夜网站| 欧美性感艳星| 国产爱豆传媒在线观看| 九色国产91popny在线|