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

    基于空載溫度場(chǎng)模擬與試驗(yàn)的冷藏車(chē)?yán)浒宀贾梅绞絻?yōu)選

    2018-01-09 08:30:13謝如鶴唐海洋陶文博劉廣海劉康佳吳俊章
    關(guān)鍵詞:冷板冷藏車(chē)冷藏

    謝如鶴,唐海洋,陶文博,劉廣海,劉康佳,吳俊章

    ?

    基于空載溫度場(chǎng)模擬與試驗(yàn)的冷藏車(chē)?yán)浒宀贾梅绞絻?yōu)選

    謝如鶴,唐海洋,陶文博,劉廣海,劉康佳,吳俊章

    (廣州大學(xué)物流與運(yùn)輸研究所,廣州 510006)

    該文通過(guò)計(jì)算流體力學(xué)(computational fluid dynamics,CFD)對(duì)常見(jiàn)的3種不同冷板布置方式(冷板頂置、冷板側(cè)置、冷板部分頂置部分側(cè)置)下的空載冷藏車(chē)廂內(nèi)溫度場(chǎng)進(jìn)行了10 h的非穩(wěn)態(tài)數(shù)值模擬研究。結(jié)果表明:在10 h的非穩(wěn)態(tài)模擬研究中,在車(chē)廂高度方向上由上往下,空氣溫度都呈現(xiàn)逐漸上升的趨勢(shì);在車(chē)廂長(zhǎng)度方向上,靠近車(chē)廂中心的區(qū)域空氣溫度更低。壁面附近空氣溫度較高,溫度梯度較大。其中,在冷板側(cè)置的方式下,廂內(nèi)整體溫度較低,各測(cè)點(diǎn)溫度在-9.5~-7 ℃范圍內(nèi),且溫度均勻性較好;在冷板部分頂置部分側(cè)置的方式中,各測(cè)點(diǎn)溫度在起始階段溫差較大,在末段溫差有所縮小,最終各測(cè)點(diǎn)溫度為-9~-7.2 ℃,測(cè)點(diǎn)整體溫度與溫差稍好于冷板頂置的方式。同時(shí),冷板側(cè)置的方式下,車(chē)廂的可用空間較?。焕浒屙斨玫姆绞降目捎萌莘e較大,但是冷板布置于車(chē)廂頂部,充注共晶液后會(huì)導(dǎo)致車(chē)輛重心的上升,降低了冷藏車(chē)在運(yùn)行時(shí)的穩(wěn)定性;冷板部分頂置部分側(cè)置的方式可用空間及重心介于冷板側(cè)置的方式與冷板頂置的方式之間。綜上所述,考慮到冷板側(cè)置的冷板布置方式對(duì)貨物的影響,優(yōu)先選用冷板部分頂置部分側(cè)置、冷板頂置的方式。綜合考慮冷板布置方式對(duì)貨物的影響,冷藏車(chē)廂內(nèi)可用容積、冷藏車(chē)廂的重心高度等因素,研制了2種不同冷板布置方式的冷板冷藏車(chē)廂,對(duì)車(chē)內(nèi)的溫度場(chǎng)進(jìn)行空載試驗(yàn),結(jié)果表明冷板部分側(cè)置、部分頂置的方式較為合理。同時(shí)經(jīng)過(guò)10 h的試驗(yàn),實(shí)測(cè)溫度與模擬仿真溫度的平均偏差為0.9 ℃,絕對(duì)誤差在2.5 ℃以內(nèi),說(shuō)明計(jì)算模型可以用于冷藏車(chē)廂內(nèi)溫度場(chǎng)的模擬。研究結(jié)果可為優(yōu)化冷板冷藏車(chē)的冷板布置方式提供依據(jù)。

    冷藏;數(shù)值模擬;溫度分布;冷板冷藏車(chē);布置方式

    0 引 言

    冷板冷藏車(chē)?yán)眯罾鋭┠毯笏鶅?chǔ)存的相變潛熱來(lái)為冷藏車(chē)廂降溫,維持貨物儲(chǔ)運(yùn)所需溫度。冷板可以在充冷站預(yù)先凍結(jié)后移至冷藏車(chē)廂中,也可以在冷藏車(chē)使用前直接將制冷機(jī)組與外接電源接通,對(duì)安裝在冷藏車(chē)廂內(nèi)的冷板充冷。常用的蓄冷劑均為低融點(diǎn)的無(wú)機(jī)鹽混合溶液,與機(jī)械式冷藏運(yùn)輸車(chē)廂比具有如下優(yōu)點(diǎn):運(yùn)行可靠、車(chē)廂內(nèi)溫度穩(wěn)定、冷板結(jié)構(gòu)簡(jiǎn)單,易于制造、經(jīng)濟(jì)性好、無(wú)噪聲污染。在越來(lái)越倡導(dǎo)可持續(xù)化發(fā)展的今天,針對(duì)冷板冷藏車(chē)的研究顯得十分必要[1]。

    針對(duì)冷藏車(chē)的研究,在車(chē)內(nèi)氣流組織方面,Moureh等[2]利用CFD軟件和雷諾應(yīng)力湍流模型對(duì)冷藏車(chē)中氣流進(jìn)行了數(shù)值模擬。Wang等[3]模擬了帶有側(cè)壁風(fēng)道和底部導(dǎo)流槽的多溫區(qū)冷藏卡車(chē)車(chē)廂內(nèi)流場(chǎng),探究了回風(fēng)口位置和風(fēng)速對(duì)流場(chǎng)的影響。在熱工特性方面,劉廣海等[4]在熱平衡法的基礎(chǔ)上,建立了冷藏車(chē)熱狀況分析體系。謝如鶴等[5]通過(guò)因素分析,建立了鐵路冷藏車(chē)傳熱系數(shù)與隔熱層厚度的優(yōu)化模型。李錦等[6]建立了冷藏車(chē)廂內(nèi)溫度隨時(shí)間變化的降溫?cái)?shù)學(xué)模型,并對(duì)影響因素進(jìn)行了分析。在堆棧方式方面,趙春江等[7]模擬了不同邊界條件和貨物不同堆棧方式車(chē)廂內(nèi)溫度場(chǎng)分布情況。張哲等[8]采用數(shù)值模擬方法,分析了不同堆碼方式對(duì)貨物區(qū)溫度場(chǎng)及流場(chǎng)的影響。在冷板冷藏車(chē)的研究方面,王延覺(jué)等[9]運(yùn)用多孔介質(zhì)理論建立番茄在冷板冷藏車(chē)內(nèi)及貨物的數(shù)學(xué)模型模擬了2種不同堆碼方式的動(dòng)態(tài)溫度場(chǎng),并對(duì)其變化及分布的特點(diǎn)進(jìn)行分析討論。張杰等[10]建立了冷板放冷的數(shù)學(xué)模型,采用顯熱容法對(duì)冷板放冷的過(guò)程進(jìn)行了分析。王奕明等[11]介紹了冷板冷藏車(chē)運(yùn)行熱負(fù)荷的計(jì)算方法,并對(duì)運(yùn)輸冷凍和冷卻貨物時(shí)的熱負(fù)荷進(jìn)行了實(shí)例計(jì)算。陳煥新等[12]采用數(shù)值模擬方法,計(jì)算了運(yùn)輸過(guò)程中冷凍貨物的動(dòng)態(tài)溫度場(chǎng)。蔡敏等[13]通過(guò)對(duì)運(yùn)輸中貨物傳熱過(guò)程的分析,建立了車(chē)廂內(nèi)復(fù)雜耦合傳熱的-數(shù)學(xué)模型;采用數(shù)值模擬的方法,計(jì)算了3種不同堆碼方式貨物動(dòng)態(tài)溫度場(chǎng)。目前,無(wú)論是常規(guī)冷藏車(chē)還是冷板冷藏車(chē),研究重點(diǎn)都集中在氣流組織、保溫材料、貨物堆碼方式等因素對(duì)溫度場(chǎng)的影響這一塊[14-21],而沒(méi)能將冷板冷藏車(chē)相對(duì)常規(guī)冷藏車(chē)的優(yōu)勢(shì)最大化,針對(duì)于冷板布置方式這一塊的研究還顯得很匱乏。

    近年來(lái),詹耀立等[22]應(yīng)用有限元法,對(duì)冷板側(cè)置及在不同冷板凍結(jié)溫度和儲(chǔ)藏食品溫度下的車(chē)廂內(nèi)溫度場(chǎng)進(jìn)行了三維模擬計(jì)算研究,發(fā)現(xiàn)了車(chē)廂內(nèi)溫度場(chǎng)分布和變化的規(guī)律及影響因素,并提出在車(chē)門(mén)附近的頂部增加兩塊頂置式冷板的建議,說(shuō)明了單一的冷板布置方式并不能滿足對(duì)車(chē)廂內(nèi)溫度場(chǎng)分布的要求,針對(duì)不同冷板的布置方式的研究具有一定的研究?jī)r(jià)值。之后,詹耀立等[23]又應(yīng)用有限元法,對(duì)冷板側(cè)置、頂置及在不同冷板凍結(jié)溫度和儲(chǔ)藏食品溫度下的鐵路冷板車(chē)車(chē)內(nèi)溫度場(chǎng)進(jìn)行數(shù)值模擬,結(jié)果顯示冷板頂置時(shí)在相同條件下總比冷板側(cè)置時(shí)溫升高,因此,可以通過(guò)改變冷板布置方式來(lái)提升車(chē)內(nèi)溫度場(chǎng)的均勻性和穩(wěn)定性。然而,詹耀立等使用的是純模擬的方式,并沒(méi)有進(jìn)行試驗(yàn)驗(yàn)證,而且研究方式單一,僅有冷板頂置與側(cè)置2種方式,另外,其研究對(duì)象是鐵路冷板冷藏車(chē),它與冷板冷藏汽車(chē)內(nèi)部結(jié)構(gòu)不同。特別是二者頂部結(jié)構(gòu)差異明顯:鐵路冷板冷藏車(chē)的頂部為圓弧棚頂結(jié)構(gòu),而冷板冷藏汽車(chē)為平頂結(jié)構(gòu)。因此,需要針對(duì)冷板冷藏汽車(chē)進(jìn)行研究,來(lái)填補(bǔ)這一方面資料的空白。在現(xiàn)用的冷板冷藏汽車(chē)中,常見(jiàn)的冷板布置方式有冷板頂置以及冷板側(cè)置2種方式,本文在2種方式的基礎(chǔ)上進(jìn)行綜合考慮,以可用容積的大小作為依據(jù),增加了將2塊冷板側(cè)置、8塊冷板頂置的第3種冷板布置方式,以期能夠更合理(仔細(xì))的觀察到冷板布置方式的改變對(duì)溫度場(chǎng)所造成的影響,并利用CFD建立了冷板冷藏車(chē)廂的傳熱模型,通過(guò)數(shù)值模擬與試驗(yàn)驗(yàn)證的方法對(duì)這3種方式進(jìn)行了比較,希望能為改善冷板的布局方式,優(yōu)化冷板冷藏車(chē)的冷板布置方式提供依據(jù)。

    1 模型的建立

    1.1 物理模型

    本文以自行研制的冷板冷藏車(chē)作為研究對(duì)象,該冷板冷藏車(chē)通過(guò)制冷機(jī)組對(duì)冷板內(nèi)的蓄冷劑(氯化鈉15~16.5%,丙三醇12.5%,水余量)充冷,使蓄冷劑完全凝固后關(guān)閉制冷機(jī)組,再通過(guò)冷板與周?chē)諝鈸Q熱制冷,它的基本資料如下:

    1)車(chē)廂外尺寸為長(zhǎng)×寬×高=6 240 mm× 2 370 mm× 2 440 mm,車(chē)廂壁厚均為120 mm,車(chē)廂內(nèi)尺寸為長(zhǎng)×寬×高=6 000 mm×2 130 mm×2 200 mm,內(nèi)部容積28.12 m3;

    2)雙開(kāi)后門(mén),開(kāi)右側(cè)門(mén)一扇,后門(mén)和側(cè)門(mén)均使用耐低溫密封膠條進(jìn)行密封,所有門(mén)均采用了不銹鋼板鎖,并采用不銹鋼鉸鏈鏈接車(chē)廂門(mén)與車(chē)壁。

    3)單塊冷板外尺寸為長(zhǎng)×寬×高=5 600×180×68 mm,板厚均為3 mm,采用鋁板制作。

    圖1 冷板冷藏車(chē)實(shí)物圖

    冷藏車(chē)箱裝配10塊冷板,在現(xiàn)用的冷板冷藏車(chē)中,常見(jiàn)的冷板布置方式有冷板頂置以及冷板側(cè)置2種方式,本文在綜合考慮這2種方式的基礎(chǔ)上,以可用容積的大小作為依據(jù),增加了將2塊冷板側(cè)置、8塊冷板頂置的第3種冷板布置方式,以期能夠更全面的觀察到冷板布置方式的改變對(duì)溫度場(chǎng)所造成的影響:1)方式1,將10塊冷板2塊一組均勻放置于車(chē)廂的頂部;2)方式2,將頂置時(shí)外側(cè)的冷板取下,分別放置在車(chē)廂的兩側(cè);3)方式3,分別在車(chē)廂兩側(cè)放置5塊冷板。

    圖2 冷藏車(chē)結(jié)構(gòu)圖

    為了便于觀察冷藏車(chē)廂內(nèi)各處的溫度隨時(shí)間的變化情況,沿車(chē)廂長(zhǎng)度方向選定3個(gè)縱斷面,截點(diǎn)分別為=100 mm、=1 500 mm、=3 000 mm,沿車(chē)并在每個(gè)縱斷面上選出5個(gè)代表點(diǎn)作為數(shù)值模擬研究點(diǎn),選出的5個(gè)代表點(diǎn)位置如圖3所示。

    本文選取溫度不均勻系數(shù)作為車(chē)廂氣流組織的評(píng)價(jià)指標(biāo),分析各溫區(qū)在不同因素水平下的溫度場(chǎng)分布差異。不均勻系數(shù)越小說(shuō)明冷藏車(chē)廂內(nèi)氣流分布越均勻,均勻的氣流分布有利于避免出現(xiàn)氣流死角、局部溫度過(guò)高現(xiàn)象。不均勻系數(shù)計(jì)算方法較為統(tǒng)一,需要在車(chē)廂載貨區(qū)選取個(gè)檢測(cè)點(diǎn),根據(jù)試驗(yàn)結(jié)果測(cè)得各點(diǎn)的溫度值,計(jì)算其算術(shù)平均值。為了研究冷藏車(chē)溫度場(chǎng)的不均勻度對(duì)蔬菜保鮮的影響,引入絕對(duì)不均勻度,的大小代表在相同測(cè)點(diǎn)數(shù)量情況下,冷藏車(chē)溫度場(chǎng)的不均勻程度,其值越大溫度場(chǎng)的不均勻程度也越大,表示為:

    式中t為第個(gè)測(cè)點(diǎn)的溫度,℃;t為個(gè)測(cè)點(diǎn)的平均溫度,℃。

    為保障冷板冷藏車(chē)內(nèi)氣流組織的均勻性,在車(chē)廂載貨區(qū)選擇3個(gè)斷面,每個(gè)斷面布置5個(gè)測(cè)點(diǎn),其布點(diǎn)位置如圖2所示。共計(jì)15個(gè)檢測(cè)點(diǎn),各檢測(cè)點(diǎn)坐標(biāo)位置如圖3所示。

    注:×為溫度測(cè)點(diǎn)。

    冷藏車(chē)廂利用Gambit軟件進(jìn)行建模,對(duì)于冷板采用六面體網(wǎng)格劃分,并對(duì)冷板局部進(jìn)行網(wǎng)格加密,對(duì)于廂內(nèi)流場(chǎng)采用非結(jié)構(gòu)化網(wǎng)格劃分,以上3個(gè)模型的網(wǎng)格單元數(shù)量維持在25萬(wàn)個(gè)左右,最大長(zhǎng)寬比在3.5左右,網(wǎng)格質(zhì)量良好,而且節(jié)約計(jì)算時(shí)間[24-26]。

    1.2 數(shù)學(xué)模型

    本文建模不考慮運(yùn)輸過(guò)程中貨物水分蒸發(fā)等傳質(zhì)因素的影響;車(chē)廂內(nèi)空氣為低速不可壓縮流體;不考慮車(chē)廂內(nèi)的熱輻射;不考慮溫度變化對(duì)空氣物性參數(shù)的影響。經(jīng)計(jì)算瑞利數(shù)Ra=1.092×1012,屬于湍流,在此采用了Boussinesq假設(shè)[27-29]。

    綜合考慮以上的各影響因素,在直角坐標(biāo)系中,計(jì)算所需的連續(xù)性方程、動(dòng)量方程、能量方程及-方程的通用形式如下:

    對(duì)應(yīng)于不同的通用變量,式(2)中各項(xiàng)參數(shù)及?模型中的系數(shù)詳見(jiàn)文獻(xiàn)[30-31]。

    1.3 冷藏車(chē)邊界條件的設(shè)定

    影響冷藏車(chē)運(yùn)行時(shí)的熱負(fù)荷的因素有很多,具體可以分為以下幾類(lèi):車(chē)體傳熱的冷消耗、漏氣的冷消耗、太陽(yáng)輻射的冷消耗、通風(fēng)的冷消耗、貨物降溫的冷消耗、車(chē)體降溫的冷消耗、貨物呼吸的冷消耗[32]。由于研究空載情況,且無(wú)需對(duì)車(chē)廂進(jìn)行通風(fēng),在模擬時(shí)只考慮車(chē)體傳熱的冷消耗、太陽(yáng)輻射的冷消耗、漏氣的冷消耗、車(chē)體降溫的冷消耗。經(jīng)計(jì)算冷藏車(chē)箱的面熱流密度=29.28 W/m3。

    冷板中蓄冷劑的比熱容為4.74 kJ/(kg·K),融解熱為175.3 kJ/kg,融點(diǎn)為?30 ℃。蓄冷劑在冷板冷藏車(chē)運(yùn)行過(guò)程中不斷吸熱,由凝固狀態(tài)逐漸轉(zhuǎn)變?yōu)橐簯B(tài)。車(chē)廂隔熱所采用的聚氨酯密度為45 kg/m3,導(dǎo)熱率為0.026 W/(m·K),比熱容為1.8 kJ/(kg·K)。假定冷板初始溫度為?34 ℃,冷板冷藏車(chē)車(chē)廂外部環(huán)境溫度設(shè)置為35 ℃,冷藏車(chē)廂內(nèi)空氣初溫及廂體初溫均為35 ℃。

    本研究在模擬計(jì)算時(shí)采用標(biāo)準(zhǔn)-模型[33],開(kāi)啟凝固/融化模塊,采用非穩(wěn)態(tài)的求解方法計(jì)算冷藏車(chē)廂的溫度變化情況,控制微分方程的離散均采用一階迎風(fēng)格式,同時(shí)采用SIMPLE算法對(duì)計(jì)算域的速度和壓力進(jìn)行耦合。時(shí)間步長(zhǎng)設(shè)為1 s,總步數(shù)為36 000。

    2 冷板冷藏車(chē)空載數(shù)值模擬仿真

    2.1 仿真結(jié)果

    本文所研究的冷藏汽車(chē)主要是用于城市之間的配送,所以運(yùn)行時(shí)間一般不超過(guò)10 h,故本文對(duì)3種不同冷板布置形式下的冷藏車(chē)廂內(nèi)溫度場(chǎng)模擬10 h。冷藏汽車(chē)空載時(shí),主要是外界熱空氣通過(guò)與車(chē)廂外壁的對(duì)流換熱、保溫層內(nèi)的熱傳導(dǎo)以及內(nèi)部空氣的自然對(duì)流將熱量傳入[12-13],從而影響車(chē)廂內(nèi)的溫度分布。

    圖4是3種冷板布置方式的冷藏車(chē)在車(chē)廂空載條件下,經(jīng)過(guò)10 h后,不同位置截面的溫度分布云圖。

    圖4a~圖4c為空載車(chē)廂長(zhǎng)度方向(=0.1 m)距車(chē)門(mén)處0.1 m的溫度分布云圖,由結(jié)果顯示,方式1中,中間部分的溫度較均勻,從?8.5 ℃升至?8 ℃,但是在靠近壁面的區(qū)域空氣的溫度梯度較大,迅速由?8 ℃升至?6 ℃左右;方式2中,中間空氣的溫度均勻性也比較好,絕大部分空氣的溫度由上至下從?9 ℃升至?8.5 ℃,同時(shí)平均溫度要比方式1中低,靠近壁面附近溫度梯度略?。环绞?中空氣在該截面中部及上部保持在?8.5 ℃左右,在截面下方及靠近壁面位置的空氣溫度較高,最高溫度為?7.5 ℃左右。這是由于縱斷面1的位置靠近車(chē)廂端部,距離冷板較遠(yuǎn),特別是下方兩側(cè)的角落,因此這部分區(qū)域的空氣溫度較高;此外,由于車(chē)廂壁面不斷與外界換熱,導(dǎo)致靠近廂壁的空氣溫度梯度變大,溫度均勻性都較差。

    圖4d~圖4f為空載車(chē)廂長(zhǎng)度方向(=1.5 m)距車(chē)門(mén)處1.5 m處的溫度分布云圖,方式1中,空氣溫度維持在?8.5~?8 ℃的范圍內(nèi),均勻性較好,但是同樣的在靠近車(chē)廂下面的區(qū)域空氣的溫度梯度較大;方式2中,該斷面內(nèi)大部分區(qū)域的空氣溫度在?9~?8.5 ℃左右,由于車(chē)廂下面的區(qū)域空氣的溫度梯度較方式1要小一些,使得整體的均勻性要相比方式1較好,同時(shí)平均溫度要比方式1中低0.5 ℃左右;在方式3中可以看出,空氣溫度在?9.5~?8.5 ℃,雖然平均溫度要比其他方式要低,這是由于在車(chē)廂兩側(cè)都放置了冷板,能同時(shí)降低兩邊壁面進(jìn)入車(chē)廂內(nèi)的熱量,但是其均勻性要差些。3種方式中冷板附近區(qū)域的流場(chǎng)溫度均最低,較高的溫度分布在距離冷板較遠(yuǎn)的地方,中部區(qū)域以及車(chē)廂底部。冷板以及壁面附近的溫度梯度較大,這是因?yàn)榭諝馐艿嚼浒謇鋮s作用最強(qiáng),冷板與周?chē)臒峥諝獯嬖谳^大溫差,不斷進(jìn)行換熱,且由于流速相對(duì)較大,換熱劇烈,因而形成較大的溫度梯度;壁面溫度梯度較大,是外界熱空氣通過(guò)與車(chē)廂外壁、保溫車(chē)的熱傳導(dǎo)以及內(nèi)部空氣的自然對(duì)流造成的。

    圖g~i為空載車(chē)廂長(zhǎng)度方向(=3 m)中間截面處的溫度分布云圖,方式1中,截面中空氣溫度維持在?9.5~?8.5 ℃的范圍內(nèi);方式2中,大部分區(qū)域空氣溫度低于?9.5 ℃;方式3中,空氣溫度在?10.5~?9.5 ℃左右。3種方式在中間截面處的平均溫度要比前2個(gè)截面處要低一些,這說(shuō)明車(chē)廂中間區(qū)域的空氣與冷板換熱效果要更好,保溫效果也更好一些。在冷板以及壁面也都存在這較大的溫度梯度。

    縱斷面3(=3 m)作為車(chē)廂長(zhǎng)度方向的中心斷面,相比于縱斷面1(=0.1 m)、縱斷面2(=1.5 m),縱斷面3中各方式截面的共晶液溫度基本一致,但截面中空氣溫度比縱斷面1和2都有所下降,這是因?yàn)橛捎谑艿娇諝饧袄浒逯g的自然對(duì)流換熱影響,空氣在密度差的驅(qū)動(dòng)下產(chǎn)生擾動(dòng),而車(chē)廂中心位置的擾動(dòng)較為強(qiáng)烈,所以該區(qū)域的空氣溫度更低。

    注:T=t+273.15,T指開(kāi)爾文溫度,t指攝氏溫度。

    在3種冷板布置方式中,冷板側(cè)置時(shí)的車(chē)廂內(nèi)的溫度場(chǎng)相對(duì)均勻一些,這是由于冷板冷藏車(chē)內(nèi)的空氣靠自然對(duì)流進(jìn)行換熱,引起流體這種運(yùn)動(dòng)的內(nèi)在力量是溫度差。流體內(nèi)部存在著溫度差,使得各部分流體的密度不同,溫度高的流體密度小,必然上升;溫度低的流體密度大,必然下降。當(dāng)冷板全部頂置時(shí),頂部的溫度最低,密度最大,使得低溫的頂部空氣會(huì)全部向下流動(dòng),而底部離冷板距離最遠(yuǎn),因此溫度最高。頂部空氣下沉的過(guò)程中,也會(huì)出現(xiàn)比較明顯的溫度分層情況。冷板部分頂置部分側(cè)置時(shí),由于冷板也主要集中于頂部,因此車(chē)內(nèi)溫度分布也具有從上到下的分層情況,不過(guò)由于在側(cè)面各放置了一塊冷板,因此相對(duì)頂置時(shí),溫度分布相對(duì)好一些,但是側(cè)置冷板也接近頂部,所以不是很明顯。而冷板全部側(cè)置時(shí),靠近兩側(cè)的溫度低,密度大,氣體主要是在靠著車(chē)廂壁的區(qū)域向下流動(dòng),向中部空氣換熱速度較慢,所以會(huì)出現(xiàn)類(lèi)似于拋物線的溫度層界線。相對(duì)于上下分層的冷板頂置以及部分頂置部分側(cè)置的方式,在中部區(qū)域,冷板側(cè)置時(shí)車(chē)廂內(nèi)的溫度分布均勻性較好。

    盡管3種冷板布置方式下的冷板數(shù)量是一定的,但是在不同的冷板布置方式下,冷藏車(chē)廂內(nèi)的可用空間是不同的,所能裝載的貨物量也將受到影響。根據(jù)國(guó)家行業(yè)標(biāo)準(zhǔn)《易腐食品冷藏鏈操作規(guī)范禽肉類(lèi)》[34]中的規(guī)定,貨物裝載時(shí),貨物與廂壁應(yīng)留有縫隙,貨物與后門(mén)之間宜保留至少10 cm距離,天花板和貨物之間宜留出至少25 cm距離,使用固定裝置防止貨物移動(dòng)。本次研究的對(duì)象為冷板冷藏車(chē),貨物距放置于側(cè)壁頂部的冷板的距離都設(shè)定為25 cm。

    按上述條件計(jì)算出冷板頂置時(shí)的可用容積為18.37 m3,冷板部分頂置部分側(cè)置時(shí)的可用容積為18.21 m3,冷板側(cè)置時(shí)的可用容積為17.43 m3。

    表1 各冷板布置方式的參數(shù)

    2.2 試驗(yàn)驗(yàn)證

    為了驗(yàn)證模擬結(jié)果的準(zhǔn)確性,試驗(yàn)以本單位自行研制的冷板冷藏車(chē)作為測(cè)試對(duì)象。由于考慮到方式3的冷板布置方式對(duì)貨物安全的影響以及可用容積的大小,優(yōu)先選用方式1、2的冷板布置形式。

    對(duì)2種方式的空載冷板冷藏車(chē)內(nèi)的溫度場(chǎng)模擬10 h后的結(jié)果進(jìn)行了試驗(yàn)驗(yàn)證,模擬和實(shí)測(cè)溫度的平均誤差為0.9 ℃,絕對(duì)誤差在2.5 ℃以內(nèi),這是由于本文中數(shù)值模擬的是非穩(wěn)態(tài)情況下冷板冷藏車(chē)廂內(nèi)溫度分布情況,冷藏車(chē)在運(yùn)行過(guò)程中外界熱環(huán)境包括空氣溫度、太陽(yáng)輻射、風(fēng)速等影響參數(shù)是不斷變化的,而在模擬過(guò)程中未考慮廂外環(huán)境的變化。模擬值與實(shí)測(cè)值偏差不大,說(shuō)明所建立的模型可以用于模擬冷板冷藏車(chē)廂體內(nèi)的溫度場(chǎng)。

    圖5 試驗(yàn)驗(yàn)證

    3 2種冷板冷藏車(chē)的試驗(yàn)對(duì)比

    試驗(yàn)儀器和材料包括:德國(guó)德圖Testo926T型熱電偶,測(cè)溫范圍為?180~350 ℃;安捷倫34972A溫度采集器(美國(guó)安捷倫科技有限公司);德國(guó)德圖Testo174溫度記錄儀,測(cè)溫范圍為?30.0~70.0 ℃,分辨率為±0.1 ℃,最大儲(chǔ)存3900個(gè)數(shù)據(jù)。測(cè)點(diǎn)布置方式如圖3所示。通過(guò)制冷機(jī)組對(duì)冷板內(nèi)的蓄冷劑充冷,使蓄冷劑完全凝固后關(guān)閉制冷機(jī)組,通過(guò)溫度記錄儀與熱電偶記錄車(chē)廂內(nèi)的空間溫度與車(chē)廂內(nèi)壁溫度,并對(duì)該種冷板布置方式下的車(chē)廂溫度場(chǎng)的溫度以及溫度均勻性進(jìn)行評(píng)估;改變部分冷板布置位置后按照以上的試驗(yàn)內(nèi)容重新試驗(yàn),將2種工況下的溫度不均勻度進(jìn)行對(duì)比。

    3.1 方式1(冷板頂置)試驗(yàn)結(jié)果

    方式1外環(huán)境在10 h內(nèi)的平均溫度為34.6 ℃,溫度在35 ℃以上的時(shí)間約為3 h。

    圖6a,制冷機(jī)組在8:00時(shí)關(guān)閉,此時(shí)蓄冷劑溫度為?31.7 ℃,已全部凝固,此時(shí)斷面1的中心溫度為?12.2 ℃,斷面2的中心溫度為?12.7 ℃,斷面3的中心溫度為?13.8 ℃。隨著時(shí)間推移,廂內(nèi)空氣溫度繼續(xù)下降,至9:00時(shí)廂內(nèi)溫度達(dá)到最低值,然后不斷上升。至18:00時(shí),斷面1的平均溫度為?8.5 ℃,斷面2的平均溫度為?8.8 ℃,斷面3的平均溫度為?9.3 ℃。

    由圖6b可知,斷面1在開(kāi)始階段溫度均勻度較好,并保持平穩(wěn)不變,10:30之后勻速上升;斷面2、3的溫度不均勻度在制冷機(jī)組關(guān)機(jī)初始階段有略微上升,然后均勻下降,12:00之后逐漸上升;測(cè)試結(jié)束時(shí),斷面2、3較斷面1具有更好的溫度均勻性,這是由于斷面1更為靠近冷藏車(chē)廂后門(mén),此處的溫升較大。

    3.2 方式2(冷板部分頂置部分側(cè)置)試驗(yàn)結(jié)果

    圖7a中,制冷機(jī)組在8:00時(shí)關(guān)閉,此時(shí)蓄冷劑溫度為?31.8 ℃,蓄冷劑已全部凝固,此時(shí)斷面1的平均溫度為?13.2 ℃,斷面2的平均溫度為?14.3 ℃,斷面3的平均溫度為?14.7 ℃。停機(jī)后廂內(nèi)空氣溫度開(kāi)始穩(wěn)定上升。至18:00時(shí),斷面1的中心為?10.2 ℃,斷面2的中心溫度為?11 ℃,斷面3的中心溫度為?11.4 ℃。

    由圖7b可知,開(kāi)始階段(8:00~10:00之間)各斷面的溫度不均勻度大小為斷面1<斷面2<斷面3,在12:00之后溫度不均勻度大小為斷面3<斷面2<斷面1;測(cè)試結(jié)束時(shí)斷面1與斷面3的不均勻度相差0.6左右。

    圖6 冷板頂置試驗(yàn)結(jié)果

    圖7 冷板部分頂置部分側(cè)置試驗(yàn)結(jié)果

    方式1與方式2中各斷面的溫度不均勻度的變化趨勢(shì)是一致的:斷面1在開(kāi)始階段較為平穩(wěn),10:00之后勻速上升;斷面2、3在12:00之前先下降,之后逐漸上升。方式1在斷面1上各點(diǎn)的溫度不均勻度大小與方式2基本相同,方式1在斷面2、3上各點(diǎn)的溫度不均勻度比方式2稍大,因此方式2的冷板位置放置方式溫度均勻性更好。

    3.3 2種方式冷板溫度變化曲線

    由圖8可知,制冷機(jī)組停機(jī)時(shí),冷板溫度基本接近;在8:00~11:00這段時(shí)間內(nèi),2種布置方式的冷板溫度變化曲線比較一致,且在9:00~11:00中冷板的溫度變化較為平緩,這是由于蓄冷劑在這個(gè)時(shí)間段內(nèi)處于融化階段,溫度維持恒定;12:00之后方式1的冷板溫度上升速度較方式2更快,測(cè)試結(jié)束時(shí)方式1冷板溫度比方式2溫度高1.5 ℃左右。這樣會(huì)使得在充冷過(guò)程中,冷板頂置的方式將耗費(fèi)更多的電量以及更多的時(shí)間,在節(jié)能與效率方面也不如部分頂置部分側(cè)置的方式。

    圖8 2種方式冷板溫度變化曲線

    4 結(jié) 論

    1)對(duì)2種方式的空載冷板冷藏車(chē)內(nèi)的溫度場(chǎng)模擬10 h后的結(jié)果進(jìn)行了試驗(yàn)驗(yàn)證,模擬和實(shí)測(cè)溫度的平均誤差為0.9 ℃,絕對(duì)誤差在2.5 ℃以內(nèi),模擬值與實(shí)測(cè)值偏差不大,說(shuō)明所建立的模型可以用于模擬冷板冷藏車(chē)廂體內(nèi)的溫度場(chǎng)。

    2)冷板部分頂置部分側(cè)置的冷板布置方式更利于車(chē)廂內(nèi)溫度場(chǎng)的溫度均勻分布,保溫時(shí)間也更長(zhǎng),而且車(chē)輛重心要較全部頂置時(shí)更低,故冷板布置應(yīng)優(yōu)先選擇部分側(cè)置部分頂置的方式。

    3)測(cè)試結(jié)束時(shí)冷板頂置的冷板溫度要比部分頂置部分側(cè)置時(shí)溫度高1.5 ℃左右。這樣會(huì)使得在充冷過(guò)程中,冷板頂置的方式將耗費(fèi)更多的電量以及更多的時(shí)間,在節(jié)能與效率方面也不如部分頂置部分側(cè)置的方式。

    4)對(duì)于冷板冷藏車(chē)而言,僅靠溫差產(chǎn)生的自然對(duì)流難以形成均勻的溫度場(chǎng)及流場(chǎng),建議采用合理的送風(fēng)形式,增強(qiáng)空氣對(duì)流,從而達(dá)到提高冷板冷藏車(chē)?yán)洳匦实哪康摹?/p>

    5)本文中數(shù)值模擬的是非穩(wěn)態(tài)情況下冷板冷藏車(chē)廂內(nèi)溫度分布情況,冷藏車(chē)在運(yùn)行過(guò)程中外界熱環(huán)境包括空氣溫度、太陽(yáng)輻射、風(fēng)速等影響參數(shù)是不斷變化的,因此為得到更為精確的數(shù)值模擬結(jié)果,應(yīng)當(dāng)綜合考慮上述因素的影響。此外,冷藏車(chē)自身的漏氣、漏熱等因素一定程度影響了空氣與冷板的換熱,還需進(jìn)行載貨試驗(yàn)進(jìn)一步驗(yàn)證冷板冷藏車(chē)廂的蓄冷性能。

    [1] 廖明義,金立芝.蓄冷板冷藏車(chē)[J]. 制冷學(xué)報(bào),1987(3):58-59.

    Liao Mingyi, Jin Lizhi. Hold-over plate refrigerated truck[J]. Journal of refrigeration, 1987(3): 58-59. (in Chinese with English abstract)

    [2] Moureh J, Tapsoba S, Derens E, et al. Air velocity characteristics within vented palletsloaded in a refrigerated vehicle with and without air ducts[J]. InternationalJournal of Refrigeration, 2009(32): 220-234.

    [3] Wang J M, Zhao X X, Li L J, et al. Air supply system sptimization for the multi-Temperature refrigerated Trucks[J]. Digital Design and Manufacturing Technology III, 2013.

    [4] 劉廣海,謝如鶴. 冷藏車(chē)熱狀況的模擬分析與實(shí)驗(yàn)研究[J]. 武漢理工大學(xué)學(xué)報(bào):交通科學(xué)與工程版,2009(4):691-694.

    Liu Guanghai, Xie Ruhe. Simulation and experiment on the heat status of refrigerated car[J]. Journal of Wuhan University of Technology(Transportation Science & Engineering) 2009(4): 691-694. (in Chinese with English abstract)

    [5] 謝如鶴,歐陽(yáng)仲志. 鐵路冷藏車(chē)傳熱系數(shù)的優(yōu)化[J]. 鐵道學(xué)報(bào),1999,21(1):34-38.

    Xie Ruhe, Ouyang Zhongzhi. Optimum for the coefficient of heat transfer of railwayrefrigerator car[J]. Journal of the China Railway Society, 1999, 21(1): 33-37. (in Chinese with English abstract)

    [6] 李錦,謝如鶴,劉廣海,等. 冷藏車(chē)降溫?cái)?shù)學(xué)模型建立與影響因素分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(6):175-182.

    Li Jin, Xie Ruhe, Liu Guanghai, et al. Establishment of cooling mathematical model and analysis ofinfluence factors for refrigerated trucks[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(6): 175-182. (in Chinese with English abstract)

    [7] 趙春江,韓佳偉,楊信廷,等. 基于CFD的冷藏車(chē)車(chē)廂內(nèi)部溫度場(chǎng)空間分布數(shù)值模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(11):168-173.

    Zhao Chunjiang, Han Jiawei, Yang Xinting, et al. Numerical simulation of temperature field distribution in refrigerated truck based on CFD[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(11): 168-173. (in Chinese with English abstract)

    [8] 張哲,郭永剛,田津津,等. 冷板冷藏汽車(chē)箱體內(nèi)溫度場(chǎng)的數(shù)值模擬及試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(25):18-24.

    Zhang Zhe, Guo Yonggang, Tian Jinjin, et al. Numericalsimulation and experiment of temperature fielddistribution in box of cold plate refrigerated truck[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(25): 18-24. (in Chinese with English abstract)

    [9] 王延覺(jué),周麗,陳煥新,等. 番茄在冷板冷藏車(chē)運(yùn)輸過(guò)程中的溫度分布[J]. 鐵道科學(xué)與工程學(xué)報(bào),2007,4(3):87-90.

    Wang Yanjue, Zhou Li, Chen Huanxin, et al.Temperature distribution of tomatoes in hold-over platerefrigerated vehicle during transporting period[J]. Journalof Railway Science and Engineering, 2007, 4(3): 87-90. (in Chinese with English abstract)

    [10] 張杰,莫羚. 機(jī)械冷板冷藏車(chē)?yán)浒宸爬鋽?shù)值模擬[J]. 制冷與空調(diào),2006(1):24-26.

    Zhang Jie, Mo Ling. The numerically calculation on discharging period of hold-over plateof hold-over plate refrigerated vehicles[J]. Refrigerationand Air-Conditioning, 2006(1): 24-26. (in Chinese with English abstract)

    [11] 王奕明,蔡敏. 冷板冷藏車(chē)運(yùn)行熱負(fù)荷的分析與研究[J].鐵道車(chē)輛,2004,42(4):4-7.

    Wang Yiming, Cai Min. Analysis and research on heat load of cool plate refrigerator car in operation[J]. Rolling Stock, 2004, 42(4): 4-7. (in Chinese with English abstract)

    [12] 陳煥新,蔡敏,張春安. 冷板冷藏車(chē)運(yùn)輸過(guò)程中冷凍貨物的溫度分布[J]. 流體機(jī)械,2005,33(4):65-68.

    Chen Huanxin, Cai Min, Zhang Chun’an. Temperature distribution of frozen cargoes in hold-over platerefrigerated vehicle during transporting period[J]. Fluid Machinery, 2005, 33(4): 65-68. (in Chinese with English abstract)

    [13] 蔡敏,陳煥新,朱先鋒. 冷板冷藏車(chē)中貨物堆碼對(duì)其溫度場(chǎng)的影響[J]. 鐵道科學(xué)與工程學(xué)報(bào),2005,2(3):78-82.

    Cai Min, Chen Huanxin, Zhu Xianfeng. Influence of stack method on the temperature field of cargoes in hold -over plate refrigerated vehicles[J]. Journal of Railway Science and Engineering, 2005, 2(3): 78-82. (in Chinesewith English abstract)

    [14] Chourasia M K, Goswami T K. Simulation of effect of stack dimensions and stacking arrangement on cool-down characteristics of potato in a cold store by computation fluid dynamics[J]. Biosystems Engineering, 2007, 96(4): 503-515.

    [15] Gan Guohui. CFD modelling of transparent bubble cavity envelopes for energy efficient greenhouses[J]. Building and Environment, 2009, 44(12): 2486-2500.

    [16] 劉廣海,謝如鶴,鄒毅峰,等. 多溫區(qū)冷藏車(chē)氣密性能影響參數(shù)理論分析與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(1):289-296.

    Liu Guanghai, Xie Ruhe, Zou Yifeng, et al. Theoretical analysis and experiment of air tightness of multi- temperature refrigerated truck[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(1): 289-296. (in Chinese with English abstract)

    [17] 徐倩,謝晶. 多溫區(qū)冷藏車(chē)廂體內(nèi)溫度場(chǎng)和速度場(chǎng)的數(shù)值模擬及優(yōu)化設(shè)計(jì)[J]. 食品與機(jī)械,2008,24(6):88-93.

    Xu Qian, Xie Jing. Optimization design and numerical simulation on the field of fluid and temperature in the multi-temperature refrigerated trucks[J]. Food & Machinery, 2008, 24(6): 88-93. (in Chinese with English abstract)

    [18] 郭嘉明,呂恩利,陸華忠,等. 保鮮運(yùn)輸車(chē)果蔬堆碼方式對(duì)溫度場(chǎng)影響的數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(13):231-236.

    Guo Jiaming, Lü Enli, Lu Huazhong, et al. Numerical simulation on temperature field effect of stack method of garden stuff forfresh-keeping transportation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(13): 231-236. (in Chinese with English abstract)

    [19] 張東霞,呂恩利,陸華忠,等. 保鮮運(yùn)輸車(chē)溫度場(chǎng)分布特性試驗(yàn)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2008,28(11):254-259.

    Zhang Dongxia, Lü Enli, Lu Huazhong, et al. Experimental study on temperature field distribution characteristics in fresh-keeping transportation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 28(11): 254-259. (in Chinese with English abstract)

    [20] 陳天及,余克志. 裝配式冷庫(kù)溫度場(chǎng)特性的實(shí)驗(yàn)研究[J].制冷,2001,20(1):21-23.

    Chen Tianji, Yu Kezhi. The experimental study on the characteristics of temperature field inside a fabricated cold store[J]. Refrigeration, 2001, 20(1): 21-23. (in Chinese with English abstract)

    [21] 李錦,謝如鶴,劉廣海,等. 車(chē)外綜合溫度條件下典型冷藏車(chē)廂內(nèi)熱穩(wěn)定性研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2012,43(8):141-147.

    Li Jin, Xie Ruhe, Liu Guanghai, et al. Thermal stability of typical refrigerated trucks compartment under outdoor sol-airtemperature conditions[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(8): 141-147. (in Chinese with English abstract)

    [22] 詹耀立,董奇志,王淞. 冷板側(cè)置時(shí)鐵路冷板冷藏車(chē)車(chē)內(nèi)溫度場(chǎng)模擬[J]. 低溫工程,2008(3):61-64.

    Zhan Yaoli, Dong Qizhi, Wang Song. Simulation on inner temperature field of railroad cold-plate refrigerator car with cold plates laid on side[J]. Cryogenics, 2008(3): 61-64. (in Chinese with English abstract)

    [23] 詹耀立,馮國(guó)會(huì). 不同冷板布置方式的冷板車(chē)內(nèi)溫度場(chǎng)模擬[J]. 沈陽(yáng)建筑大學(xué)學(xué)報(bào),2008,24(5):859-862.

    Zhan Yaoli, Feng Guohui. The simulation of temperature

    field inside cold -plate refrigeratorcar following different setting of cold plate[J]. Journal of Shenyang Jianzhu University, 2008, 24(5): 859-862. (in Chinese with English abstract)

    [24] 謝晶,湯毅,王金鋒,等. 三維流體力學(xué)預(yù)測(cè)風(fēng)機(jī)不同布置形式對(duì)冷庫(kù)氣流的影響[J]. 食品工業(yè)科技,2011,32(11):349-351.

    Xie Jing, Tang Yi, Wang Jinfeng, et al. Computational fluid dynamics simulation of influence of different arrangements of fans to the cold store[J]. Science and Technology of Food Industry, 2011, 32(11): 349-351. (in Chinese with English abstract)

    [25] 湯毅,謝晶,王金鋒,等. CFD預(yù)測(cè)風(fēng)機(jī)擺設(shè)形式對(duì)冷庫(kù)的影響[J]. 食品與機(jī)械,2012,28(1):124-128.

    Tang Yi, Xie Jing, Wang Jinfeng, et al. CFD simulationof air flow distribution with different arrangements of fanthough a cold store[J]. Food & Machinery, 2012, 28(1): 124-128. (in Chinese with English abstract)

    [26] Foster A M, Swain M J, Barrett R, et al. Experimental verification of analytical and CFD predictions of infiltration through cold store entrances[J]. International Journal of Refrigeration, 2003, 26(8): 918-925.

    [27] 王永生,楊?lèi)?ài)玲,張華,等. 基于FLUENT 軟件的臥式冷柜內(nèi)溫度場(chǎng)的數(shù)值模擬[C]//第五屆全國(guó)制冷空調(diào)新技術(shù)研討會(huì)論文集. 中國(guó)廈門(mén),中國(guó)制冷學(xué)會(huì),2008.

    Wang Yongsheng, Yang Ailing, Zhang Hua, et al. A numerical study of the temperature distribution inside a horizontal freezer with fluent software[C]//The Fifth National Refrigeration and air conditioning new technology conference. Xiamen, Chinese Association of Refrigeration, 2008. (in Chinese with English abstract)

    [28] Rodríguez-Bermejoa J, Barreirob P, Roblaa J I, et al. Thermal study of a transport container[J]. Journal of Food Engineering, 2007, 80(2): 517-527.

    [29] Foster A M, Swain M J, Barrett R, et al. Experimental verification of analytical and CFD predictions of infiltration through cold store entrances[J]. International Journal of Refrigeration, 2003, 26(8): 918-925.

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

    [31] 陶文銓. 數(shù)值傳熱學(xué)[M]. 西安:西安交通大學(xué)出版社,2001.

    [32] 羅榮武. 冷藏集裝箱的冷消耗計(jì)算[J]. 集裝箱化,2001(11):22-23.

    [33] 韓占忠,王敬,蘭小平. 流體工程仿真計(jì)算實(shí)例與應(yīng)用.北京理工大學(xué)出版社,2004.

    [34] SB/T 10731-2012,易腐食品冷藏鏈操作規(guī)范[S].

    Optimization of cold-plate location in refrigerated vehicles based on simulation and test of no-load temperature field

    Xie Ruhe,Tang Haiyang, Tao Wenbo, Liu Guanghai, Liu Kangjia, Wu Junzhang

    (,,510006,)

    In this paper, we studied the uniformity of temperature field distribution of different cold plate arrangement, and learned the influence of the arrangement of cold plate on the distribution uniformity of temperature field in the carriage, so as to provide basis for further research. In this paper, the unsteady numerical simulation of the temperature field in the no-load refrigerated compartment with three different cold plate arrangements ( cold plate set by top, by side,and partially by top partially by side) was carried for 10h by computing CFD (computational fluid dynamics ).It was found that the air temperature showed a rising trend from top to bottom in the height direction of the carriage. In the direction of the length of the car, the regional air temperature near the center of the car was lower. The air temperature near the wall was higher, and the temperature gradient was larger. When the cold plates were set on side, the whole compartment temperature was low after 10 hours, and the temperature of each measuring point was in the range of -9.5--7 ℃. The temperature uniformity was good.When the cold plates were partially set by top and partially set on side, the temperatures of different measuring points were vastly different at the beginning, and the temperature difference was reduced at the end, with the final temperature of each measuring point in the range of -9-7.2 ℃. In this situation, the overall temperature and the temperature differences were slightly better than that with cold plate completely set on top. At the same time, when the cold plates were set on side, the available space of the compartment was small, and there was the risk that there might be collisions between the goods and the cold plate in the operation of the refrigerated vehicle, which might be a threat to the safe transportation of the goods. The available area of the cold plate was larger, but because the cold plates were arranged on the top of the vehicle, the Eutectic liquid will lead to the increase of the center of gravity of the vehicle, and the stability of the refrigerated vehicle in operation could be reduced. The available space and the height of the center of gravity of the compartment with cold plates partially set on top and partially set on side were between the two arrangements above. In view of the above, considering the effect of the arrangement of the cold plates on the goods, the cold plates partially set on top and partially set on side were selected as a priority. Considering the influence of the arrangement of the cold plate on the cargo, the available volume in refrigerated vehicles and the height of center of gravity of refrigerated vehicles, two kinds of cold plate refrigerated vehicles with different cold plates arrangements were prepared. At the same time, after 10 h, the deviation between the measured temperature and the simulated temperature was about 2.5 ℃, which showed that the calculation model can be used to simulate the temperature field in refrigerated compartment. The results provide a basis for improving the arrangement of the cold plate and optimizing the arrangement of the cold plate of the cold plate refrigerated vehicle.

    cold storage; numerical simulation; temperature distribution; cold plate; arrangement form

    10.11975/j.issn.1002-6819.2017.24.038

    U272.5

    A

    1002-6819(2017)-24-0290-09

    2017-06-03

    2017-12-11

    廣東省科技計(jì)劃項(xiàng)目(No.2016B020205004,No.2017B090907028,No.611138153066);國(guó)家科技支撐計(jì)劃農(nóng)業(yè)領(lǐng)域項(xiàng)目(新型冷藏及保鮮運(yùn)輸技術(shù)與裝備/2013BAD19B01-1)

    謝如鶴,教授,博士,博士生導(dǎo)師,主要從事冷鏈物流方面研究。Email:583385752@qq.com

    謝如鶴,唐海洋,陶文博,劉廣海,劉康佳,吳俊章.基于空載溫度場(chǎng)模擬與試驗(yàn)的冷藏車(chē)?yán)浒宀贾梅绞絻?yōu)選[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(24):290-298. doi:10.11975/j.issn.1002-6819.2017.24.038 http://www.tcsae.org

    Xie Ruhe, Tang Haiyang, Tao Wenbo, Liu Guanghai, Liu Kangjia, Wu Junzhang. Optimization of cold-plate location in refrigerated vehicles based on simulation and test of no-load temperature field[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(24): 290-298. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.24.038 http://www.tcsae.org

    猜你喜歡
    冷板冷藏車(chē)冷藏
    東風(fēng)汽車(chē)股份簽約500臺(tái)冷藏車(chē)!
    福特phve冷板檢測(cè)定位pin自動(dòng)分揀機(jī)控制系統(tǒng)的研究
    利用光伏發(fā)電制冷的冷藏車(chē)設(shè)計(jì)選型
    煤氣與熱力(2021年4期)2021-06-09 06:17:00
    某組合型冷板的結(jié)構(gòu)設(shè)計(jì)與優(yōu)化
    歐洲冷藏車(chē)主流技術(shù)介紹
    機(jī)載大長(zhǎng)寬比風(fēng)冷均溫冷板優(yōu)化設(shè)計(jì)
    食物冷藏不要超過(guò)多少天
    哪些應(yīng)該放冷藏?哪些應(yīng)該放冷凍?哪些不用放冰箱?
    媽媽寶寶(2017年2期)2017-02-21 01:21:04
    冷藏保溫車(chē)發(fā)展?jié)摿Ρ患ぐl(fā)
    再談冷藏保溫車(chē):市場(chǎng)已升溫
    男女无遮挡免费网站观看| 韩国精品一区二区三区| 亚洲国产欧美日韩在线播放| 国产亚洲av片在线观看秒播厂| 日本欧美国产在线视频| 国产熟女欧美一区二区| 国产成人a∨麻豆精品| 国精品久久久久久国模美| 久久久精品国产亚洲av高清涩受| 最黄视频免费看| 精品国产一区二区三区久久久樱花| 国产熟女欧美一区二区| 黄频高清免费视频| 国产女主播在线喷水免费视频网站| 国产女主播在线喷水免费视频网站| 亚洲国产欧美一区二区综合| 久久久久久久久久久久大奶| 久久久亚洲精品成人影院| 51午夜福利影视在线观看| 欧美激情高清一区二区三区 | 精品少妇黑人巨大在线播放| 精品一区二区三区av网在线观看 | 人人妻人人澡人人爽人人夜夜| 熟女av电影| 午夜福利网站1000一区二区三区| 国产精品 欧美亚洲| 老司机亚洲免费影院| 涩涩av久久男人的天堂| 大话2 男鬼变身卡| 免费av中文字幕在线| 老司机亚洲免费影院| 国产精品久久久av美女十八| 久久99一区二区三区| 国产成人欧美在线观看 | 精品亚洲成国产av| 精品免费久久久久久久清纯 | 免费日韩欧美在线观看| 免费少妇av软件| 国产av码专区亚洲av| 日本wwww免费看| 啦啦啦在线免费观看视频4| 好男人视频免费观看在线| 久久热在线av| 亚洲中文av在线| 国产有黄有色有爽视频| 男人操女人黄网站| 18禁裸乳无遮挡动漫免费视频| 在线观看国产h片| 老司机亚洲免费影院| 国产精品欧美亚洲77777| 人人妻人人添人人爽欧美一区卜| 欧美日韩精品网址| 国产亚洲av片在线观看秒播厂| 亚洲久久久国产精品| 亚洲av国产av综合av卡| 成人三级做爰电影| 国产成人a∨麻豆精品| 侵犯人妻中文字幕一二三四区| 丝袜在线中文字幕| 乱人伦中国视频| 午夜福利一区二区在线看| 成年人免费黄色播放视频| 中文字幕人妻丝袜制服| 亚洲精品日韩在线中文字幕| 亚洲天堂av无毛| 狂野欧美激情性xxxx| 国产精品欧美亚洲77777| 高清不卡的av网站| 在现免费观看毛片| 国产精品一二三区在线看| 九九爱精品视频在线观看| tube8黄色片| 精品人妻熟女毛片av久久网站| 亚洲国产欧美一区二区综合| 精品久久蜜臀av无| 在线免费观看不下载黄p国产| 欧美变态另类bdsm刘玥| 哪个播放器可以免费观看大片| 黄色怎么调成土黄色| 久久久久久久大尺度免费视频| 久久鲁丝午夜福利片| 欧美精品一区二区免费开放| 成年女人毛片免费观看观看9 | 国产精品一区二区在线观看99| 肉色欧美久久久久久久蜜桃| 男女免费视频国产| 国产黄频视频在线观看| 狂野欧美激情性xxxx| 天堂8中文在线网| 午夜91福利影院| 曰老女人黄片| 女性被躁到高潮视频| 国产一级毛片在线| 国产精品麻豆人妻色哟哟久久| 欧美老熟妇乱子伦牲交| 久久久久久久久久久免费av| 人体艺术视频欧美日本| 麻豆av在线久日| 国产国语露脸激情在线看| 男女边吃奶边做爰视频| 久久国产亚洲av麻豆专区| av女优亚洲男人天堂| 美女午夜性视频免费| 中文天堂在线官网| 欧美精品av麻豆av| 天天躁夜夜躁狠狠躁躁| 成人亚洲欧美一区二区av| 亚洲成国产人片在线观看| 免费日韩欧美在线观看| a级片在线免费高清观看视频| 91国产中文字幕| 黄色怎么调成土黄色| 国产爽快片一区二区三区| 国产av精品麻豆| 午夜激情久久久久久久| 国产日韩欧美视频二区| 亚洲av日韩精品久久久久久密 | 亚洲国产欧美在线一区| 在线观看www视频免费| 九草在线视频观看| 18禁裸乳无遮挡动漫免费视频| 日本爱情动作片www.在线观看| 免费黄色在线免费观看| av一本久久久久| 亚洲精品第二区| 欧美黑人欧美精品刺激| 最新在线观看一区二区三区 | 日韩中文字幕欧美一区二区 | 日韩不卡一区二区三区视频在线| 久久久久网色| 亚洲成av片中文字幕在线观看| 在线观看免费午夜福利视频| 超碰成人久久| a级片在线免费高清观看视频| 欧美日韩一级在线毛片| 超碰成人久久| 人妻 亚洲 视频| 纵有疾风起免费观看全集完整版| av又黄又爽大尺度在线免费看| 亚洲国产精品成人久久小说| 黄色一级大片看看| 精品一区二区免费观看| 免费黄色在线免费观看| 天天添夜夜摸| av一本久久久久| 晚上一个人看的免费电影| 亚洲成人av在线免费| 午夜激情久久久久久久| 纯流量卡能插随身wifi吗| 色网站视频免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产黄色视频一区二区在线观看| 岛国毛片在线播放| 熟女少妇亚洲综合色aaa.| 性色av一级| 久久久久久久久久久免费av| 国产成人啪精品午夜网站| 婷婷成人精品国产| 一级毛片黄色毛片免费观看视频| 欧美精品av麻豆av| tube8黄色片| 18在线观看网站| 欧美精品人与动牲交sv欧美| 久久久久久久久免费视频了| 91成人精品电影| 男女下面插进去视频免费观看| 国产一区二区 视频在线| 九九爱精品视频在线观看| 满18在线观看网站| 久久综合国产亚洲精品| 少妇精品久久久久久久| 国产 一区精品| 咕卡用的链子| 高清视频免费观看一区二区| 中文精品一卡2卡3卡4更新| 亚洲五月色婷婷综合| 啦啦啦 在线观看视频| 狂野欧美激情性xxxx| 国产乱来视频区| 免费高清在线观看日韩| 亚洲国产精品一区三区| 伊人久久国产一区二区| av网站免费在线观看视频| 观看美女的网站| 亚洲少妇的诱惑av| 青草久久国产| a 毛片基地| 捣出白浆h1v1| 自拍欧美九色日韩亚洲蝌蚪91| 免费看av在线观看网站| 99国产综合亚洲精品| 国产成人av激情在线播放| 亚洲av在线观看美女高潮| 中文天堂在线官网| 美国免费a级毛片| 国产不卡av网站在线观看| bbb黄色大片| 国产 精品1| 美女高潮到喷水免费观看| 熟女少妇亚洲综合色aaa.| 女的被弄到高潮叫床怎么办| 亚洲国产欧美一区二区综合| 女人爽到高潮嗷嗷叫在线视频| 成人午夜精彩视频在线观看| 操美女的视频在线观看| 日韩中文字幕欧美一区二区 | 桃花免费在线播放| 国产亚洲一区二区精品| 美女大奶头黄色视频| 国精品久久久久久国模美| 免费少妇av软件| 免费观看a级毛片全部| 欧美 亚洲 国产 日韩一| 久久亚洲国产成人精品v| 亚洲免费av在线视频| 99热全是精品| 亚洲精品aⅴ在线观看| 亚洲色图综合在线观看| 久久精品熟女亚洲av麻豆精品| 69精品国产乱码久久久| 欧美精品亚洲一区二区| 国产成人精品在线电影| 黄色怎么调成土黄色| 国产xxxxx性猛交| 日本av免费视频播放| 免费在线观看视频国产中文字幕亚洲 | 久久久精品国产亚洲av高清涩受| 久久久久久免费高清国产稀缺| 免费女性裸体啪啪无遮挡网站| 久久人人爽av亚洲精品天堂| 久久精品aⅴ一区二区三区四区| 日本一区二区免费在线视频| 亚洲欧美激情在线| 亚洲国产精品999| 69精品国产乱码久久久| 久久午夜综合久久蜜桃| 国产精品欧美亚洲77777| av网站免费在线观看视频| 国产精品免费大片| 2021少妇久久久久久久久久久| 国产欧美日韩一区二区三区在线| 777米奇影视久久| 日韩av免费高清视频| 精品国产国语对白av| 久久精品国产a三级三级三级| 91精品三级在线观看| 午夜av观看不卡| 中文字幕av电影在线播放| 国产深夜福利视频在线观看| 成人亚洲欧美一区二区av| 中国三级夫妇交换| 国产熟女欧美一区二区| 成年女人毛片免费观看观看9 | 免费观看av网站的网址| 大码成人一级视频| 亚洲欧美一区二区三区国产| 尾随美女入室| 男人操女人黄网站| 精品国产国语对白av| 久久精品熟女亚洲av麻豆精品| 国产国语露脸激情在线看| 亚洲五月色婷婷综合| 免费女性裸体啪啪无遮挡网站| 啦啦啦在线免费观看视频4| 另类亚洲欧美激情| 亚洲综合色网址| 18禁观看日本| 99精品久久久久人妻精品| 国产精品99久久99久久久不卡 | 国产 精品1| 精品国产一区二区三区久久久樱花| 日韩成人av中文字幕在线观看| 久久精品熟女亚洲av麻豆精品| 99热网站在线观看| 丝瓜视频免费看黄片| 亚洲人成电影观看| 最近的中文字幕免费完整| 男女之事视频高清在线观看 | 欧美97在线视频| 视频区图区小说| 建设人人有责人人尽责人人享有的| 久久人妻熟女aⅴ| 热re99久久国产66热| 欧美人与善性xxx| 欧美国产精品一级二级三级| 一区福利在线观看| 一级片'在线观看视频| 亚洲少妇的诱惑av| 久久久久精品久久久久真实原创| 飞空精品影院首页| 欧美精品一区二区免费开放| 性色av一级| 久久久亚洲精品成人影院| 免费在线观看视频国产中文字幕亚洲 | 久久久久久免费高清国产稀缺| 亚洲av欧美aⅴ国产| 一级毛片我不卡| 女人久久www免费人成看片| 搡老乐熟女国产| 色综合欧美亚洲国产小说| 丝袜在线中文字幕| 亚洲美女黄色视频免费看| 一级a爱视频在线免费观看| 欧美xxⅹ黑人| 免费日韩欧美在线观看| 丝袜美腿诱惑在线| 亚洲国产最新在线播放| 99国产综合亚洲精品| 国产1区2区3区精品| 宅男免费午夜| 欧美日韩亚洲国产一区二区在线观看 | 性色av一级| 男的添女的下面高潮视频| 又黄又粗又硬又大视频| 大陆偷拍与自拍| 精品少妇久久久久久888优播| videos熟女内射| 久久免费观看电影| 日韩av免费高清视频| 亚洲精品国产区一区二| 国产精品偷伦视频观看了| 九色亚洲精品在线播放| 久久精品国产亚洲av高清一级| 丰满少妇做爰视频| 在线天堂最新版资源| 亚洲婷婷狠狠爱综合网| 成人亚洲欧美一区二区av| 91精品国产国语对白视频| 精品人妻熟女毛片av久久网站| 丰满迷人的少妇在线观看| h视频一区二区三区| xxx大片免费视频| 无遮挡黄片免费观看| 色婷婷av一区二区三区视频| 国产亚洲精品第一综合不卡| 亚洲av男天堂| 高清在线视频一区二区三区| 十分钟在线观看高清视频www| 人人妻,人人澡人人爽秒播 | 99久久99久久久精品蜜桃| 不卡视频在线观看欧美| 国产99久久九九免费精品| 免费在线观看完整版高清| 成人影院久久| 深夜精品福利| 看免费av毛片| 黄色视频不卡| 嫩草影院入口| 看十八女毛片水多多多| 在线观看免费高清a一片| 在线天堂中文资源库| 国产成人av激情在线播放| 免费在线观看完整版高清| 在线观看国产h片| 欧美变态另类bdsm刘玥| 久久影院123| 国产精品国产av在线观看| 丰满乱子伦码专区| 亚洲人成网站在线观看播放| 久久ye,这里只有精品| 热99久久久久精品小说推荐| 一级毛片电影观看| 欧美黄色片欧美黄色片| 人人妻人人添人人爽欧美一区卜| 老汉色av国产亚洲站长工具| 人人妻人人添人人爽欧美一区卜| 欧美激情极品国产一区二区三区| 卡戴珊不雅视频在线播放| 精品酒店卫生间| 久久97久久精品| 大香蕉久久网| 色播在线永久视频| 国产熟女欧美一区二区| 色婷婷久久久亚洲欧美| 亚洲av欧美aⅴ国产| 老司机在亚洲福利影院| 国产精品国产三级专区第一集| 亚洲熟女毛片儿| 国产黄频视频在线观看| 久久久久久久大尺度免费视频| 欧美日韩一级在线毛片| 日韩制服丝袜自拍偷拍| av片东京热男人的天堂| 国产免费又黄又爽又色| 国产精品三级大全| 另类亚洲欧美激情| 成年女人毛片免费观看观看9 | 午夜福利一区二区在线看| 一级毛片电影观看| 欧美人与善性xxx| 看免费av毛片| 午夜福利在线免费观看网站| 我要看黄色一级片免费的| 肉色欧美久久久久久久蜜桃| 欧美日韩av久久| 1024视频免费在线观看| 婷婷色av中文字幕| 国产 一区精品| 中文字幕亚洲精品专区| 国产成人免费观看mmmm| 日韩 欧美 亚洲 中文字幕| 日韩成人av中文字幕在线观看| 国产毛片在线视频| 久久99热这里只频精品6学生| 捣出白浆h1v1| av国产久精品久网站免费入址| 精品国产露脸久久av麻豆| 亚洲综合精品二区| 亚洲av欧美aⅴ国产| 少妇 在线观看| 久久99精品国语久久久| 女人被躁到高潮嗷嗷叫费观| 色婷婷av一区二区三区视频| 精品一区二区三卡| av一本久久久久| 亚洲男人天堂网一区| 看免费av毛片| 人体艺术视频欧美日本| a 毛片基地| 丝袜人妻中文字幕| 婷婷成人精品国产| bbb黄色大片| 国产亚洲精品第一综合不卡| 免费久久久久久久精品成人欧美视频| 18禁国产床啪视频网站| 黄色毛片三级朝国网站| 欧美激情高清一区二区三区 | 男人操女人黄网站| 亚洲国产精品国产精品| 亚洲国产欧美在线一区| 宅男免费午夜| 卡戴珊不雅视频在线播放| 老司机影院毛片| 亚洲四区av| 日韩av不卡免费在线播放| 久久久精品区二区三区| 国产精品秋霞免费鲁丝片| 无遮挡黄片免费观看| 午夜久久久在线观看| 国产又色又爽无遮挡免| 亚洲色图综合在线观看| 中文字幕人妻丝袜制服| 亚洲国产精品成人久久小说| 街头女战士在线观看网站| 国产爽快片一区二区三区| 一级毛片 在线播放| 国产又爽黄色视频| 国产日韩一区二区三区精品不卡| www.av在线官网国产| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品一区二区三区在线| 亚洲四区av| 在线观看免费日韩欧美大片| 欧美日韩精品网址| 大话2 男鬼变身卡| 天美传媒精品一区二区| 好男人视频免费观看在线| 日本黄色日本黄色录像| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| 桃花免费在线播放| 亚洲成人免费av在线播放| 哪个播放器可以免费观看大片| 国产无遮挡羞羞视频在线观看| 18在线观看网站| 久热爱精品视频在线9| 国产成人精品久久久久久| 精品久久久精品久久久| videosex国产| 久久97久久精品| 人人妻人人爽人人添夜夜欢视频| 欧美少妇被猛烈插入视频| 国产xxxxx性猛交| 国产精品亚洲av一区麻豆 | a级毛片在线看网站| 国产老妇伦熟女老妇高清| 免费高清在线观看日韩| 久久国产亚洲av麻豆专区| 精品人妻熟女毛片av久久网站| 国产一区二区 视频在线| kizo精华| 欧美精品av麻豆av| 欧美变态另类bdsm刘玥| 高清av免费在线| 日韩视频在线欧美| 爱豆传媒免费全集在线观看| 天天躁夜夜躁狠狠躁躁| 最黄视频免费看| 欧美成人精品欧美一级黄| 成人黄色视频免费在线看| 欧美亚洲日本最大视频资源| 黄片无遮挡物在线观看| 久久精品久久久久久久性| 成人手机av| 国产无遮挡羞羞视频在线观看| 亚洲色图 男人天堂 中文字幕| av线在线观看网站| 久久久精品94久久精品| 久久久久久久久久久免费av| 777久久人妻少妇嫩草av网站| 国产精品一国产av| 亚洲国产毛片av蜜桃av| 国产日韩一区二区三区精品不卡| 色网站视频免费| 亚洲欧美日韩另类电影网站| 国产一区二区在线观看av| 在线观看免费视频网站a站| av国产精品久久久久影院| 免费女性裸体啪啪无遮挡网站| 国产乱来视频区| 日韩一卡2卡3卡4卡2021年| 人妻 亚洲 视频| 在线观看免费高清a一片| 精品福利永久在线观看| 久久久久久久国产电影| av网站在线播放免费| 国产探花极品一区二区| 国产精品嫩草影院av在线观看| 日本av免费视频播放| av在线老鸭窝| 久久综合国产亚洲精品| 侵犯人妻中文字幕一二三四区| 欧美人与性动交α欧美精品济南到| 丁香六月天网| 丰满迷人的少妇在线观看| 人体艺术视频欧美日本| 建设人人有责人人尽责人人享有的| 性少妇av在线| 久久亚洲国产成人精品v| 天天躁日日躁夜夜躁夜夜| 国产极品天堂在线| 97在线人人人人妻| av天堂久久9| 91av网站免费观看| 精品国产超薄肉色丝袜足j| 国产亚洲欧美98| 欧美激情 高清一区二区三区| 精品久久久精品久久久| 成人手机av| 国产成人系列免费观看| 久久青草综合色| 人成视频在线观看免费观看| 好男人电影高清在线观看| 精品久久久久久成人av| 亚洲无线在线观看| 变态另类成人亚洲欧美熟女 | 国产精品久久久人人做人人爽| 国产1区2区3区精品| 日本 欧美在线| 国产精品影院久久| xxx96com| 狠狠狠狠99中文字幕| 久久久久久人人人人人| 亚洲精品国产一区二区精华液| 亚洲中文av在线| 少妇 在线观看| 夜夜夜夜夜久久久久| 国产91精品成人一区二区三区| 午夜福利影视在线免费观看| 制服诱惑二区| 午夜福利欧美成人| 精品高清国产在线一区| 一本久久中文字幕| 在线视频色国产色| 国产99久久九九免费精品| 久久久久久久久久久久大奶| 男人舔女人的私密视频| 日韩欧美国产在线观看| 性色av乱码一区二区三区2| 亚洲一区二区三区色噜噜| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品九九99| 777久久人妻少妇嫩草av网站| 国内精品久久久久久久电影| 免费在线观看完整版高清| 在线十欧美十亚洲十日本专区| 久久香蕉激情| 国产高清激情床上av| 精品卡一卡二卡四卡免费| 禁无遮挡网站| 午夜福利欧美成人| 丝袜美足系列| 给我免费播放毛片高清在线观看| 国产又色又爽无遮挡免费看| 精品高清国产在线一区| 一区二区三区精品91| 99re在线观看精品视频| 国产成人精品久久二区二区免费| 十八禁人妻一区二区| 亚洲国产毛片av蜜桃av| 免费一级毛片在线播放高清视频 | www.自偷自拍.com| 久久热在线av| 国产成人精品久久二区二区91| 自线自在国产av| 精品第一国产精品| 高清黄色对白视频在线免费看| 亚洲av成人不卡在线观看播放网| 日韩欧美国产在线观看| 久久精品国产亚洲av高清一级| 成人国产一区最新在线观看| 久久热在线av| 人成视频在线观看免费观看| 国产男靠女视频免费网站| 免费在线观看日本一区| 18美女黄网站色大片免费观看| 波多野结衣av一区二区av| 国产成人免费无遮挡视频| 久久精品成人免费网站| 免费少妇av软件| 麻豆一二三区av精品| 这个男人来自地球电影免费观看| 亚洲第一av免费看| 国产亚洲精品综合一区在线观看 | 在线观看午夜福利视频|