多溫區(qū)冷藏車(chē)廂溫度場(chǎng)影響因素的數(shù)值模擬

官彬彬，李 強(qiáng)，張強(qiáng)華

（浙江科技學(xué)院 機(jī)械與汽車(chē)工程學(xué)院，浙江 杭州 310023）

多溫區(qū)冷藏車(chē)廂溫度場(chǎng)影響因素的數(shù)值模擬

官彬彬，李 強(qiáng)，張強(qiáng)華

（浙江科技學(xué)院 機(jī)械與汽車(chē)工程學(xué)院，浙江 杭州 310023）

針對(duì)多溫區(qū)冷藏車(chē)車(chē)廂設(shè)計(jì)參數(shù)的優(yōu)化問(wèn)題，建立了冷量在多溫區(qū)流動(dòng)的物理模型，推導(dǎo)出冷藏車(chē)多溫區(qū)溫度變化的熱平衡動(dòng)態(tài)方程。分析風(fēng)扇出風(fēng)口溫度，總體勻流板孔隙率，勻流板高度等因素對(duì)車(chē)廂總體溫度不均勻程度以及儲(chǔ)藏可靠性的影響。以楊梅為例數(shù)值模擬出評(píng)價(jià)指標(biāo)與因素間的相互影響規(guī)律。仿真結(jié)果表明：風(fēng)扇出風(fēng)口溫度為-2.5℃，總體勻流板孔隙率為0.215，高度為40mm為本實(shí)驗(yàn)多溫區(qū)冷藏車(chē)廂的較優(yōu)設(shè)計(jì)參數(shù)。

多溫區(qū)；冷藏車(chē)廂；溫度場(chǎng)；數(shù)值模擬

1 引言

近年來(lái)，由于冷凍冷藏業(yè)的飛速發(fā)展，冷凍冷藏食品的種類(lèi)和數(shù)目日益增多，以往使用的單溫區(qū)冷藏車(chē)越來(lái)越難以滿(mǎn)足市場(chǎng)需求，多溫區(qū)冷藏車(chē)則成為了高?？蒲械难芯繜狳c(diǎn)[1-3]。同時(shí)，計(jì)算流體力學(xué)（CFD）技術(shù)在冷鏈物流行業(yè)應(yīng)用越來(lái)越廣泛，用于幫助工程技術(shù)人員剖析和解決問(wèn)題[4-6]。為有效準(zhǔn)確地控制冷藏車(chē)廂內(nèi)溫度場(chǎng)，提高儲(chǔ)藏介質(zhì)質(zhì)量可靠性，在運(yùn)輸過(guò)程中，不僅僅考慮堆放方式[7]，風(fēng)扇的送風(fēng)方式、轉(zhuǎn)速及出風(fēng)口的溫度，還與總體勻流板孔隙率、儲(chǔ)藏介質(zhì)的呼吸熱等因素對(duì)車(chē)廂溫度場(chǎng)的影響有關(guān)。冷風(fēng)在車(chē)廂內(nèi)進(jìn)行流動(dòng)循環(huán)時(shí)，孔隙率大小合適的勻流板能有效減少湍流現(xiàn)象，降低車(chē)廂內(nèi)溫度不均勻程度[8]。此外，勻流板的高度也是影響車(chē)廂內(nèi)溫度場(chǎng)的重要因素。本實(shí)驗(yàn)通過(guò)分析風(fēng)扇出風(fēng)口溫度、總體勻流板孔隙率、勻流板高度對(duì)車(chē)廂總體溫度不均勻程度以及儲(chǔ)藏可靠性的影響，為多溫區(qū)冷藏車(chē)廂的參數(shù)設(shè)計(jì)提供參考。

2 數(shù)學(xué)模型構(gòu)建

2.1 冷藏車(chē)廂的熱量流動(dòng)

圖1為冷藏車(chē)廂的熱量流動(dòng)的縱向剖視圖，F(xiàn)表示

風(fēng)扇，從左到右依次為冷藏區(qū)A、B和C，箭頭表示熱量流動(dòng)方向，其中冷藏區(qū)A、B和C的出入口處均為有可變孔隙率的勻流板，其厚度為δi；δg表示隔板厚度；Qi1表示通過(guò)車(chē)廂體傳入冷藏區(qū)的熱量；Qi2表示通過(guò)水蒸氣傳入冷藏區(qū)的熱量；Qi3表示儲(chǔ)藏介質(zhì)的呼吸熱；Qf表示風(fēng)扇產(chǎn)生的熱量；Qfti表示制冷空氣通過(guò)勻流板傳入冷藏區(qū)ABC的熱量。i=a、b、c，表示冷藏區(qū)A、B和C（下同）。

2.2 熱平衡動(dòng)態(tài)方程

圖1 冷藏車(chē)廂的熱量流動(dòng)縱向剖視圖

根據(jù)熱力學(xué)第一定律，建立多溫冷藏車(chē)降溫過(guò)程的熱平衡動(dòng)態(tài)方程。假設(shè)各冷藏區(qū)內(nèi)的空氣定壓比、空氣密度均為定值；多溫區(qū)之間無(wú)空氣泄漏[9-10]。

根據(jù)整個(gè)車(chē)廂的降溫過(guò)程模擬公式如下：

式中，Q0為制冷機(jī)組實(shí)際制冷量（W）；Ci為冷藏各區(qū)內(nèi)的空氣的定壓比熱，單位J（kg·K）；ρi為各冷藏區(qū)內(nèi)的空氣密度（kg/m3）；Vi為各冷藏區(qū)的體積（m3）；τ為時(shí)間（s）。

以單個(gè)溫區(qū)為研究對(duì)象，各個(gè)溫區(qū)的降溫過(guò)程可表示為：

式中，Vfi為各冷藏區(qū)勻流板表面空氣流速（m）；aω、為車(chē)廂外表面以及車(chē)廂內(nèi)各冷藏區(qū)的內(nèi)表面換熱系數(shù)（W）；δi為除兩溫區(qū)隔板外車(chē)廂體各傳熱表面的當(dāng)量厚度（m）；δg為兩溫區(qū)隔板傳熱表面的當(dāng)量厚度（m）；λi為車(chē)廂各傳熱表面的當(dāng)量導(dǎo)熱系數(shù)，單位W/(m·K)；tω、tk、ta、tb、tc為車(chē)廂各外表面綜合溫度、車(chē)廂空氣溫度、各冷藏區(qū)內(nèi)空氣溫度（℃）；

Faj、Fbj、Fcj為除兩溫區(qū)隔板外車(chē)廂體各個(gè)溫區(qū)車(chē)廂的總傳熱面積（m2），j表示傳熱表面的數(shù)量；Fg1、Fg2為兩溫區(qū)間隔板表面?zhèn)鳠崦娣e（m2）。

式中，βi為車(chē)廂冷藏各區(qū)的漏氣倍數(shù)；γi為車(chē)廂各冷藏區(qū)的水蒸氣凝結(jié)熱（J/kg）；φiω、φin為各冷藏區(qū)車(chē)廂外、內(nèi)的空氣相對(duì)濕度（%）；xiω、xin為各冷藏區(qū)車(chē)廂外、內(nèi)空氣含濕量（g/kg）。

式中，M為車(chē)載食品貨物的質(zhì)量（kg）；

H為車(chē)載食品貨物的呼吸熱，單位J（kg·24h）；假設(shè)各冷藏區(qū)內(nèi)運(yùn)輸介質(zhì)所占的體積相等，運(yùn)輸介質(zhì)質(zhì)量以及呼吸熱近似相等。

Ffi為各冷藏區(qū)勻流板表面面積（m2）；εi為各冷藏區(qū)勻流板孔隙率；Pe為風(fēng)扇額定功率（w）；η為風(fēng)扇工作效率（%）；Ψ為風(fēng)扇熱轉(zhuǎn)換系數(shù)，取值為1。

聯(lián)立式（1）-式（9）得：

上式中，m1、m2、n1、n2、l0、l1、l2、z0、z1、z2表示與以上參數(shù)相關(guān)的常數(shù)。

聯(lián)立式（10）-式（12）解得：

P0、P1、P2表示與m1、m2、n1、n2、l0、l1、l2、z0、z1、z2有關(guān)的常數(shù)；式（13）表示冷藏區(qū)ABC的平均溫度隨時(shí)間變化的三階常系數(shù)非齊次微分方程。

ri1、ri2、ri2為對(duì)應(yīng)的特征方程的三個(gè)根。

根據(jù)式（13）解得溫度隨時(shí)間變化為

式中C1，C2，C3為常數(shù)。

3 綜合評(píng)價(jià)指標(biāo)

3.1 總體溫度不均勻系數(shù)

總體溫度不均勻系數(shù)是指溫度場(chǎng)達(dá)到相對(duì)穩(wěn)態(tài)時(shí)各冷藏車(chē)廂內(nèi)溫度的不均勻程度，其值越小，表明車(chē)廂溫度越均勻，越適合儲(chǔ)藏。Bi表示各個(gè)區(qū)的溫度不均勻系數(shù)，各個(gè)溫區(qū)取18個(gè)測(cè)點(diǎn)；B表示車(chē)廂總體溫度不均勻系數(shù)。

Wi表示第i車(chē)廂的權(quán)重系數(shù)，根據(jù)各冷藏車(chē)廂均勻性的比重，各冷藏區(qū)分別取0.5，0.3和0.2。

3.2 模擬儲(chǔ)藏可靠性系數(shù)

模擬儲(chǔ)藏可靠性系數(shù)是指運(yùn)輸介質(zhì)儲(chǔ)藏溫度在規(guī)定的溫度范圍內(nèi)所占的比例。某運(yùn)輸介質(zhì)適宜的儲(chǔ)藏溫度為t1～t2，若要較長(zhǎng)時(shí)間儲(chǔ)藏，必須嚴(yán)格控制冷庫(kù)溫度，庫(kù)溫波動(dòng)太大導(dǎo)致運(yùn)輸介質(zhì)的營(yíng)養(yǎng)流失、變質(zhì)等問(wèn)題，一般需要冷藏的運(yùn)輸介質(zhì)主要保存在溫度為-1℃-8℃范圍的冷藏區(qū)。

Ri表示冷藏區(qū)i運(yùn)輸介質(zhì)的儲(chǔ)藏可靠性。R表示模擬運(yùn)輸介質(zhì)儲(chǔ)藏可靠性系數(shù)。R值越大，說(shuō)明該運(yùn)輸介質(zhì)儲(chǔ)藏可靠性越好。

4 仿真結(jié)果與比較

以楊梅為儲(chǔ)藏運(yùn)輸介質(zhì)為例，不考慮呼吸作用對(duì)楊梅質(zhì)量的影響，楊梅的適宜儲(chǔ)藏溫度為-0.5℃-0℃，楊梅具有規(guī)則的幾何形狀，且按一定方式放置，故可將其設(shè)置成多孔介質(zhì)，楊梅的擺放方式為風(fēng)道架空擺放，這種擺放方式能使楊梅與楊梅之間存在較大空隙，空氣流動(dòng)均勻，減少局部滯留，并減少氣漩的產(chǎn)生，能使車(chē)廂內(nèi)的溫度分布更加均勻。風(fēng)扇的送風(fēng)方式為上送下回，這種方式能減緩熱量的流動(dòng)，減少氣流的漩渦，使空氣流動(dòng)相對(duì)平穩(wěn)。多溫區(qū)的布局方式為縱向布局，這種布局方式有效地縮短了各溫區(qū)長(zhǎng)度，有利于形成良好的氣流組織，冷藏效果好，裝卸方便。

4.1 勻流板孔隙率的影響

當(dāng)風(fēng)扇出風(fēng)口溫度270.15K，風(fēng)扇轉(zhuǎn)子轉(zhuǎn)速為460rad/s，勻流板高度為40mm，研究總體勻流板孔隙率對(duì)模擬楊梅儲(chǔ)藏可靠性系數(shù)以及車(chē)廂總體溫度不均勻系數(shù)的影響。

先考慮多孔板的孔隙率的取值，為了增加進(jìn)入冷藏區(qū)C的冷風(fēng)流量，應(yīng)使冷藏區(qū)C的上下多孔板的孔隙率足夠的大，冷藏區(qū)A和B保證儲(chǔ)藏溫度均勻度的情況下儲(chǔ)藏溫度下盡可能的取低。又由于勻流板工藝的限制，孔隙率最小從0.1附近取。冷藏區(qū)A上下孔板取值0.08、0.10和0.12；冷藏區(qū)B的孔隙率取值應(yīng)比冷藏區(qū)A略大一點(diǎn)，但不能太大，否則影響進(jìn)入冷藏區(qū)C的冷風(fēng)量。冷藏區(qū)B上下孔板取值0.15、0.20和0.25。冷藏區(qū)C上下孔板取0.6。既要保證溫度波動(dòng)范圍不要太高，又要保證靜態(tài)溫度不均勻系數(shù)足夠小，為了避免出入孔板孔隙率不均衡對(duì)靜態(tài)溫度不均勻系數(shù)以及溫度不均勻系數(shù)波動(dòng)指數(shù)的影響，設(shè)置時(shí)將出入口孔板孔隙率設(shè)置為一樣。根據(jù)以上結(jié)論確定各個(gè)冷藏區(qū)勻流板孔隙率。ε表示總體勻流板孔隙率。

仿真結(jié)果如圖2所示，當(dāng)冷藏區(qū)ABC選取的孔隙率不同時(shí)，總體勻流板孔隙率的大小也不同。隨著總體勻流板孔隙率的增大，模擬楊梅儲(chǔ)藏可靠性系數(shù)總體先減小后增大，第一階段，由于進(jìn)入冷藏區(qū)的冷風(fēng)增加，楊梅降溫速度變快，導(dǎo)致楊梅處在適宜儲(chǔ)藏溫度的時(shí)間范圍

減少；第二階段，楊梅處在適宜儲(chǔ)藏溫度的時(shí)間基本不變，而整體降溫時(shí)間繼續(xù)加快。冷藏區(qū)A孔隙率為0.1，冷藏區(qū)B孔隙率為0.15，冷藏區(qū)C孔隙率為0.6，總體孔隙率為0.215附近出現(xiàn)最大值為19.02%，總體變化幅度為7.15%。

圖2 綜合評(píng)價(jià)指標(biāo)隨總體勻流板孔隙率的變化規(guī)律

模擬楊梅儲(chǔ)藏可靠性系數(shù)與總體勻流板孔隙率的擬合多項(xiàng)式為：

車(chē)廂總體溫度不均勻系數(shù)先增大后減小，變化趨勢(shì)逐漸變緩。仿真結(jié)果表明：冷藏區(qū)A孔隙率為0.08，冷藏區(qū)B孔隙率為0.15，冷藏區(qū)C孔隙率為0.6，總體孔隙率為0.205時(shí)，車(chē)廂冷藏區(qū)的溫度不均勻系數(shù)達(dá)到最小，總體變化幅度為4.89%。第一階段，由于降溫速度變快導(dǎo)致車(chē)廂總體溫度不均勻系數(shù)增大；第二階段，由于最后溫度降低使得車(chē)廂總體溫度不均勻系數(shù)減小。

車(chē)廂總體溫度不均勻系數(shù)與總體勻流板孔隙率的擬合多項(xiàng)式為：

冷藏車(chē)廂總體孔隙率對(duì)車(chē)廂冷藏區(qū)溫度不均勻系數(shù)影響的變化幅度只有4.89%，而冷藏車(chē)廂總體孔隙率對(duì)模擬楊梅儲(chǔ)藏可靠性系數(shù)影響的變化幅度有7.15%，說(shuō)明冷藏區(qū)A孔隙率為0.1，冷藏區(qū)B孔隙率為0.15，冷藏區(qū)C孔隙率為0.6，總體孔隙率為0.215相對(duì)冷藏區(qū)A孔隙率為0.08，冷藏區(qū)B孔隙率為0.2，冷藏區(qū)C孔隙率為0.6，總體孔隙率為0.205較優(yōu)。

4.2 風(fēng)扇出風(fēng)口溫度的影響

當(dāng)冷藏區(qū)A孔隙率為0.1，冷藏區(qū)B孔隙率為0.15，冷藏區(qū)C孔隙率為0.6，總體孔隙率為0.215，風(fēng)扇轉(zhuǎn)子轉(zhuǎn)速為460rad/s，勻流板高度為40mm，研究風(fēng)扇出風(fēng)口溫度對(duì)模擬楊梅儲(chǔ)藏可靠性系數(shù)以及冷藏車(chē)廂溫度不均勻系數(shù)的影響。首先確定溫度的取值范圍，溫度太低或者太高都會(huì)影響楊梅的儲(chǔ)藏。楊梅適宜的冷藏溫度為-0.5℃～0℃，考慮到要克服冷藏車(chē)箱內(nèi)的各種熱量，將風(fēng)扇出風(fēng)口溫度定為-4.5℃～-2℃，取溫度間距為0.5。

圖3 綜合評(píng)價(jià)指標(biāo)隨風(fēng)扇出風(fēng)口溫度的變化規(guī)律

從圖3可以看出，在符合楊梅儲(chǔ)藏要求的前提下，隨著風(fēng)扇出口溫度的增大，模擬楊梅儲(chǔ)藏可靠性系數(shù)先增大后減小。這是因?yàn)檫m當(dāng)增加溫度，能使冷藏車(chē)廂的最終儲(chǔ)藏溫度處于-1℃到-0.5℃之間，當(dāng)溫度大于-2.5℃時(shí)，冷藏區(qū)C的楊梅儲(chǔ)藏可靠性為0，模擬楊梅可靠性系數(shù)下降。風(fēng)扇出風(fēng)口溫度為-2.5℃時(shí)取得最大，其值為19.25%，模擬楊梅儲(chǔ)藏可靠性系數(shù)的變化幅度為69.75%。模擬楊梅儲(chǔ)藏可靠性系數(shù)與風(fēng)扇出風(fēng)口溫度的擬合多項(xiàng)式為：

在符合楊梅儲(chǔ)藏要求的前提下，隨著風(fēng)扇出風(fēng)口溫度的增大，車(chē)廂總體溫度不均勻系數(shù)總體上不斷減小，這是因?yàn)轱L(fēng)扇出風(fēng)口溫度增大，導(dǎo)致總體降溫速度減小，總體溫度均勻程度更好。仿真結(jié)果表明：在溫度為-2℃時(shí)達(dá)到最小值，變化幅度為7.38%。冷藏車(chē)廂溫度不均勻系數(shù)與風(fēng)扇出風(fēng)口溫度的擬合多項(xiàng)式為：

在符合楊梅儲(chǔ)藏要求的前提下，隨著風(fēng)扇出口溫度的增大，車(chē)廂冷藏區(qū)的溫度不均勻系數(shù)的變化幅度只有7.38%，而模擬楊梅儲(chǔ)藏可靠性的變化幅度為69.75%。因此風(fēng)扇出風(fēng)口溫度設(shè)置為-2.5℃時(shí)較優(yōu)。

4.3 勻流板高度的影響

當(dāng)冷藏區(qū)A孔隙率為0.1，冷藏區(qū)B孔隙率為0.15，冷藏區(qū)C孔隙率為0.6，總體孔隙率為0.215，風(fēng)扇轉(zhuǎn)子轉(zhuǎn)速為460Rad/s，風(fēng)扇出風(fēng)口溫度為-2.5℃時(shí)，研究勻流板高度對(duì)模擬楊梅儲(chǔ)藏可靠性系數(shù)以及車(chē)廂總體溫度不均勻系數(shù)的影響。勻流板高度等同于勻流板上孔的深度。由于勻流板制作工藝的限制，取勻流板的高度在20mm到50mm之間，間隔為5mm來(lái)分析勻流板高度的影響。

圖4 綜合評(píng)價(jià)指標(biāo)隨勻流板高度的變化規(guī)律

從圖4可以看出，在符合勻流板工藝要求的前提下，隨著勻流板高度的增大，模擬楊梅儲(chǔ)藏可靠性系數(shù)先增大后減小，在高度為40mm時(shí)達(dá)到最大，其值為19.25%，其變化幅度為56.38%。這是因?yàn)楫?dāng)勻流板高度小于40mm時(shí)，增加勻流板高度，增大了勻流板體積，增加了勻流板消耗的熱量，而風(fēng)扇出口溫度不變，冷藏區(qū)的溫度升高，因此模擬楊梅儲(chǔ)藏可靠性系數(shù)增大。當(dāng)勻流板高度大于40mm時(shí)，勻流板消耗的熱量太多導(dǎo)致冷藏區(qū)C儲(chǔ)藏非常不可靠，可靠性系數(shù)明顯下降。模擬楊梅儲(chǔ)藏可靠性系數(shù)與勻流板高度的擬合多項(xiàng)式為：

在符合勻流板工藝要求的前提下，隨著勻流板高度的增大，車(chē)廂總體溫度不均勻系數(shù)增大，仿真結(jié)果表明：在高度為20時(shí)最小，總體變化幅度為8.14%。勻流板高度增加，導(dǎo)致冷藏區(qū)最終溫度升高，使楊梅呼吸熱越大，進(jìn)而導(dǎo)致車(chē)廂總體溫度不均勻系數(shù)增大。車(chē)廂總體溫度不均勻系數(shù)與風(fēng)扇出風(fēng)口溫度的擬合多項(xiàng)式為：

在符合勻流板工藝要求的前提下，隨著勻流板高度的變化，車(chē)廂總體溫度不均勻系數(shù)的變化幅度只有8.14%，而模擬楊梅儲(chǔ)藏可靠性系數(shù)的變化幅度為56.38%。因此勻流板高度設(shè)置為40mm時(shí)最佳。

5 結(jié)論

（1）隨著總體勻流板孔隙率的增大，模擬楊梅儲(chǔ)藏可靠性系數(shù)總體先減小后增大。車(chē)廂總體溫度不均勻系數(shù)先增大后減小。

（2）在符合楊梅儲(chǔ)藏要求前提下，隨著風(fēng)扇出口溫度的增大，模擬楊梅儲(chǔ)藏可靠性系數(shù)先增大后減小。車(chē)廂總體溫度不均勻系數(shù)總體上不斷減小。

（3）在符合勻流板工藝要求的前提下，隨著勻流板高度的增大，模擬楊梅儲(chǔ)藏可靠性系數(shù)先增大后減小，總體溫度不均勻系數(shù)不斷增大。

（4）風(fēng)扇出風(fēng)口溫度對(duì)模擬食物儲(chǔ)藏可靠性系數(shù)的影響程度較大，勻流板高度對(duì)冷藏車(chē)廂溫度不均勻系數(shù)的影響程度較大，由于溫度達(dá)到了相對(duì)穩(wěn)態(tài)時(shí)，運(yùn)輸介質(zhì)完成了熱交換，各處溫度值基本相等，總體溫度均勻性較好，而且各區(qū)選取的18個(gè)測(cè)點(diǎn)均空氣流動(dòng)性較好，熱交換活動(dòng)比較充分。

（5）總體孔隙率為0.215，風(fēng)扇出風(fēng)口溫度為-2.5℃，勻流板高度為40mm為本實(shí)驗(yàn)所探究的多溫區(qū)冷藏車(chē)車(chē)廂設(shè)計(jì)參數(shù)優(yōu)化問(wèn)題的較優(yōu)設(shè)計(jì)方案。值得指出的是，在設(shè)計(jì)參數(shù)的過(guò)程中，應(yīng)優(yōu)先考慮風(fēng)扇出風(fēng)口溫度，其次是勻流板高度，最后考慮勻流板的孔隙率優(yōu)化組合。

[1]謝晶,徐倩,方恒和.多溫區(qū)冷藏車(chē)熱負(fù)荷計(jì)算的研究[J].食品與機(jī)械,2007,23(4):98-101.

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

[3]趙鑫鑫,王家敏,李麗娟,等.多溫區(qū)冷藏車(chē)回風(fēng)導(dǎo)軌對(duì)廂內(nèi)溫度場(chǎng)的影響分析[J].食品與機(jī)械,2014,30(1):149-154.

[4]Xia B,Sun D W.Application of computational fluid dynamics (CFD)in the food industry:a review[J].Computers and Electronics in Agriculture,2002,(34):5-24.

[5]郭嘉明,陸華忠,呂恩利,等.CFD在食品儲(chǔ)運(yùn)設(shè)備研發(fā)中的應(yīng)用[J].農(nóng)機(jī)化研究,2012,(8):219-222.[6]張哲,郭永剛,田津津.冷板冷藏汽車(chē)箱體內(nèi)溫度場(chǎng)的數(shù)值模擬及試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(25):18-24.

[7]彭夢(mèng)瑰,段國(guó)權(quán),丁力行,等.機(jī)冷車(chē)在不同貨物堆碼方式下溫度場(chǎng)分布的數(shù)值模擬[J].鐵道學(xué)報(bào),2007,29(5):37-42.

[8]Moureh J,Flick D.Airflow Pattern and Temperature Distribution in a Typical Refrigerated Truck Configuration Loaded with Pallets[J].International Journal of Refrigeration,2004,27 (5):464-474.

[9]楊培志.冷板冷藏車(chē)整車(chē)傳熱系數(shù)的計(jì)算分析[J].制冷與空調(diào),2004,4(6):53-55.

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

Numerical Simulation of Temperature Field Influence Factors of Multi-zone Refrigeration Lorry

Guan Binbin,Li Qiang,Zhang Qianghua
(School of Mechanical&Automobile Engineering,Zhejiang University of Science&Technology,Hangzhou 310023,China)

In this paper,in view of the optimization of the design parameters of the multi-zone refrigeration lorry,we built the physical model of the thermal dynamics across the multiple temperature zones and derived the dynamic equation of the thermal balance of the refrigeration lorry.Then we analyzed the factors influencing the even thermal distribution and storage reliability in the lorry carriage.At the end,in the case of the red bayberry,we simulated the evaluation index and the mutual influence between the factors and determined the optimal design parameters for the lorry carriage.

multiple temperature zone;refrigeration carriage;temperature field;numerical simulation

U16;F224

A

1005-152X(2016)10-0087-05

10.3969/j.issn.1005-152X.2016.10.022

2016-09-03

國(guó)家科技合作專(zhuān)項(xiàng)（2013DFA31920）；浙江省科技計(jì)劃項(xiàng)目（2015C32075）；浙江省自然科學(xué)基金項(xiàng)目（LY13E050023）

官彬彬（1993-），男，碩士研究生，研究方向：新能源汽車(chē)及機(jī)電一體化模擬；李強(qiáng)（1979-），通訊作者，男，博士，副教授，研究方向：汽車(chē)機(jī)電一體化控制技術(shù)。