顏建春，魏 海，謝煥雄，游兆延

·農(nóng)產(chǎn)品加工工程·

筒狀固定床花生通風(fēng)干燥性能指標(biāo)模擬與分析

顏建春，魏 海，謝煥雄※，游兆延

（1. 農(nóng)業(yè)農(nóng)村部南京農(nóng)業(yè)機(jī)械化研究所，南京 210014；2. 中美花生生產(chǎn)工程技術(shù)聯(lián)合實(shí)驗(yàn)室，南京 210014）

為了解筒狀固定床花生干燥機(jī)理、作業(yè)性能，確定合理的結(jié)構(gòu)和通風(fēng)參數(shù)，根據(jù)干燥過(guò)程花生莢果和介質(zhì)空氣間的熱質(zhì)傳遞關(guān)系，以PDE模型為理論基礎(chǔ)，建立了適用于筒狀固定床花生通風(fēng)干燥計(jì)算機(jī)模擬的離散模型，該模型可計(jì)算花生干燥實(shí)時(shí)狀態(tài)及批次干燥耗時(shí)、不均勻度、生產(chǎn)率、單位質(zhì)量能耗等干燥指標(biāo)。經(jīng)試驗(yàn)驗(yàn)證，模型模擬結(jié)果與試驗(yàn)結(jié)果基本相符，料層花生平均含水率和溫度模擬值和試驗(yàn)值的相關(guān)系數(shù)＞0.975，模型模擬可用于筒狀床花生干燥過(guò)程分析。在此基礎(chǔ)上，分析了單位面積通風(fēng)量、筒狀固定床外徑、內(nèi)徑變化對(duì)上述指標(biāo)的影響。結(jié)果表明：受介質(zhì)空氣溫度降低和相對(duì)濕度增加影響，內(nèi)層物料干燥起始時(shí)間和干燥速率存在一定的滯后性，但單位面積通風(fēng)量沿通風(fēng)方向逐漸增大的特性對(duì)內(nèi)層物料的干燥滯后有較好改善；隨著單位面積通風(fēng)量增加，干燥不均勻度明顯降低，生產(chǎn)率亦有顯著提高，但單位質(zhì)量能耗增幅較大；筒狀床外徑增加或內(nèi)徑減小都可增加生產(chǎn)率，降低單位質(zhì)量干燥能耗，但干燥不均勻問(wèn)題很難解決。為進(jìn)一步確定最優(yōu)的結(jié)構(gòu)和通風(fēng)參數(shù)，采用均勻設(shè)計(jì)模擬試驗(yàn)和綜合加權(quán)評(píng)分法，得出筒狀固定床外徑2.75 m，內(nèi)徑0.935 m，外進(jìn)風(fēng)面單位面積通風(fēng)量0.36 m3/(m2·s)時(shí)干燥效果最優(yōu)，此時(shí)干燥時(shí)耗39.2 h、生產(chǎn)率86.55 kg/h、單位質(zhì)量能耗5.87×106J/kg、干燥不均勻度1.54%。該研究可為筒狀固定床干燥設(shè)備設(shè)計(jì)優(yōu)化提供技術(shù)支撐。

農(nóng)產(chǎn)品；干燥；花生；固定床；模擬

0 引 言

花生是優(yōu)勢(shì)優(yōu)質(zhì)油料作物和重要食品蛋白資源，中國(guó)的花生產(chǎn)量和種植面積常年穩(wěn)居世界前列[1-2]。近年來(lái)受農(nóng)業(yè)人口轉(zhuǎn)移和土地流轉(zhuǎn)政策影響，花生機(jī)械化收獲水平逐年提高，收獲時(shí)間趨于集中，短期內(nèi)易堆積大量鮮果，而農(nóng)村原有的晾曬設(shè)施和勞動(dòng)力已難以滿足急劇膨脹的干燥需求，遭逢連續(xù)陰雨天氣時(shí)，霉變腐爛常有發(fā)生，黃曲霉毒素污染較嚴(yán)重[2-3]，亟需一種經(jīng)濟(jì)實(shí)用的花生干燥機(jī)。以美國(guó)為代表的花生規(guī)?；N植的發(fā)達(dá)國(guó)家，已有發(fā)展成熟的花生干燥工藝與裝備，但其干燥環(huán)節(jié)與收獲環(huán)節(jié)銜接緊密[3-8]，裝載干燥時(shí)的花生（帶殼）含水率一般20%左右[8-12]，作業(yè)原理為固定床單向通風(fēng)干燥，采用液化石油氣燃燒供熱[13-14]。但在中國(guó)農(nóng)村，為搶農(nóng)時(shí)趕種下茬作物，很少有時(shí)間進(jìn)行長(zhǎng)時(shí)間田間晾曬，收獲時(shí)以鮮花生為主，干燥降水幅度大，干燥效果欠佳，兼之中國(guó)能源價(jià)格昂貴，設(shè)備投入和使用成本高，很難在中國(guó)農(nóng)村推廣。少量花生種植大戶采用常規(guī)固定床干燥設(shè)備作為補(bǔ)充或應(yīng)急干燥手段，同批次物料干燥不均勻性嚴(yán)重，且經(jīng)濟(jì)成本高昂[15-18]，限制了花生機(jī)械化干燥的推廣應(yīng)用。

作者所在研究團(tuán)隊(duì)針對(duì)中國(guó)花生收獲方式、物料特性，提出了基于筒狀結(jié)構(gòu)固定床的花生通風(fēng)干燥技術(shù)，介質(zhì)空氣穿過(guò)外圓柱面沿半徑方向向內(nèi)穿過(guò)花生料層過(guò)程中，空氣溫度降低和濕度升高導(dǎo)致的干燥能力降低，與單位面積通風(fēng)量沿著通風(fēng)方向逐漸增加而引起的干燥能力加強(qiáng)相互抵消，達(dá)到均勻干燥的目的。論文采用計(jì)算機(jī)模擬的方法，分析了筒狀固定床不同結(jié)構(gòu)和通風(fēng)參數(shù)下的作業(yè)性能，為花生產(chǎn)地機(jī)械化干燥裝備研發(fā)提供理論支撐。

1 原理與模型

1.1 筒狀固定床通風(fēng)干燥原理

筒狀固定床通風(fēng)干燥裝置主要包括進(jìn)風(fēng)室、載料室、出風(fēng)室等，進(jìn)風(fēng)室與載料室及進(jìn)風(fēng)室與出風(fēng)室之間通過(guò)沖孔板隔開(kāi)，如圖1所示。干燥作業(yè)時(shí)，花生物料裝滿環(huán)狀筒倉(cāng)，物料底部承料板和頂部蓋板均為實(shí)心板。環(huán)狀筒倉(cāng)外側(cè)為進(jìn)風(fēng)室，內(nèi)側(cè)為出風(fēng)室，介質(zhì)空氣從進(jìn)風(fēng)室底部風(fēng)道進(jìn)入進(jìn)風(fēng)室，充分混合后，穿過(guò)外沖孔板進(jìn)入環(huán)狀物料層，介質(zhì)空氣在從外層向內(nèi)層物料穿透過(guò)程中溫度逐漸降低，相對(duì)濕度逐漸增加，吸濕能力逐漸降低，但通風(fēng)面積逐漸減少，單位面積通風(fēng)量逐漸增加。而單位面積通風(fēng)量的增加，可增加介質(zhì)空氣對(duì)物料的吸濕能力和干燥速度。因此，設(shè)置合理的結(jié)構(gòu)和通風(fēng)參數(shù)，使溫度降低濕度升高導(dǎo)致的干燥能力降低和單位面積通風(fēng)量增加導(dǎo)致干燥能力加強(qiáng)盡可能平衡，從而達(dá)到整個(gè)花生料層均勻干燥的目的。

注：箭頭表示空氣流動(dòng)方向；φa為筒狀固定床外徑，m；φb為筒狀固定床內(nèi)徑，m；H為筒倉(cāng)高度，m。

1.2 花生干燥模型

1.2.1 花生干燥理論模型

花生莢果總體由果仁和果殼組成，收獲后的新鮮花生含水率高，此時(shí)果仁緊貼果殼內(nèi)壁，干燥過(guò)程中含水率逐漸降低，果殼和果仁之間產(chǎn)生間隙并逐漸增大，阻礙了熱量的向內(nèi)傳遞和水分的向外遷移。此外，果仁種皮和果殼內(nèi)壁質(zhì)構(gòu)緊密，對(duì)水分向外遷移也造成了較大阻礙。因此花生莢果傳熱傳質(zhì)復(fù)雜，很難針對(duì)各組分具體特征形成單個(gè)花生莢果干燥模型，現(xiàn)有研究基本將花生莢果看作一個(gè)類圓柱的實(shí)體，采用半經(jīng)驗(yàn)方程和試驗(yàn)歸納的方法，得出薄層干燥方程[3,9,18]。

本文針對(duì)筒狀固定床花生物料，根據(jù)前人利用偏微分方程組干燥模型（PDE，partial differential equation）模擬固定床物料干燥過(guò)程的研究經(jīng)驗(yàn)，為簡(jiǎn)化計(jì)算，對(duì)固定床花生干燥過(guò)程中傳熱、傳質(zhì)作如下假設(shè)：忽略單個(gè)花生莢果內(nèi)部的溫度梯度；忽略花生莢果間的熱傳導(dǎo)；器壁絕熱，并忽略其熱容量；在足夠短的時(shí)間間隔內(nèi)濕空氣和花生莢果比熱視為常量[19-23]。

根據(jù)筒狀床料層花生莢果和介質(zhì)空氣間的熱質(zhì)傳遞基本理論，可得空氣質(zhì)平衡方程、空氣熱平衡方程、花生莢果熱平衡方程，如式（1）～（3）。

式中為薄層花生與筒狀固定床外徑的距離，m；為干燥時(shí)間，h；為薄層花生平均干基含水率（無(wú)量綱）；為薄層花生溫度，℃；為介質(zhì)空氣穿過(guò)薄層花生后的絕對(duì)濕度（無(wú)量綱）；為介質(zhì)空氣穿過(guò)薄層花生后的溫度，℃；a為單位送風(fēng)面積介質(zhì)空氣通風(fēng)量，m3/(m2·s)；pb為花生莢果體積密度，kg/m3；a為干空氣密度，kg/m3；為床層花生空隙率（無(wú)量綱）；a為干空氣比熱，J/(kg·℃)；v為水蒸氣比熱，J/(kg·℃)；w為液態(tài)水比熱，J/(kg·℃)；p為花生床層比熱，J/(kg·℃)；T為花生莢果-空氣傳熱系數(shù)，J/(m2·℃)；fg為花生莢果中水分的汽化熱，J/kg；s為床層花生比表面積，m2/m3。

根據(jù)已有的薄層花生干燥研究基礎(chǔ)[3]，花生莢果質(zhì)平衡半經(jīng)驗(yàn)方程，如式（4）～（7）。

其中：

式（4）～（7）適用溫度范圍25℃≤≤50℃，適用單位面積通風(fēng)量范圍0.2 m3/(m2·s)≤a≤1.2 m3/(m2·s)；為干燥常數(shù)；M為花生平衡干基含水率（無(wú)量綱）。

由于氣流穿過(guò)環(huán)狀料層過(guò)程中總通風(fēng)量不變，而通風(fēng)面積沿氣流方向逐漸減少，則有單位面積通風(fēng)量與料層位置關(guān)系

式中0為筒狀床外柱面單位面積通風(fēng)量，m3/(m2·s)；a為筒狀固定床外徑，m。

1.2.2 偏微分方程組離散化

將筒狀固定床看成若干個(gè)沿直徑方向的管狀薄層組成，每一個(gè)薄層作為一個(gè)節(jié)點(diǎn)；將烘干時(shí)間按時(shí)序分成若干等分，每個(gè)等分作為一個(gè)節(jié)點(diǎn)。當(dāng)從外向內(nèi)通風(fēng)時(shí)，模擬從筒狀固定床最外層向內(nèi)開(kāi)始，依次計(jì)算不同位置節(jié)點(diǎn)處薄層花生含水率、溫度和介質(zhì)空氣溫濕度，將計(jì)算結(jié)果作為下一節(jié)點(diǎn)的輸入值，再去計(jì)算下一節(jié)點(diǎn)處的各變量值。直到床層花生平均含水率達(dá)到貯藏要求為止。干燥過(guò)程中物料層位置節(jié)點(diǎn)、時(shí)間節(jié)點(diǎn)劃分如圖2。

注：i為料層厚度方向上的節(jié)點(diǎn)數(shù)，表示第i個(gè)薄層；j時(shí)間維度上的節(jié)點(diǎn)數(shù)，表示第j個(gè)時(shí)間等分；m為厚度方向的最大節(jié)點(diǎn)數(shù)，n為時(shí)間維度的最大節(jié)點(diǎn)數(shù)；t為干燥時(shí)間，h；x為薄層花生與筒狀固定床外徑的距離，m。

將筒狀固定床花生沿半徑方向均勻分成若干單元，每個(gè)單元厚度為Δ，=·Δ，總厚為·Δ。同時(shí)，將干燥時(shí)間劃分成若干段，每段時(shí)長(zhǎng)為Δ，=·Δ，料床平均含水率值達(dá)到烘干終止含水率時(shí)終止模擬。采用有限差分法，對(duì)上述理論模型進(jìn)行離散化，整理后得

式中eq為當(dāng)前料層在當(dāng)前溫濕度條件下上一個(gè)時(shí)間節(jié)點(diǎn)對(duì)應(yīng)的當(dāng)量干燥時(shí)間，h。其中：

通過(guò)式（12）、（9）可計(jì)算得=(+1)·Δ時(shí)刻第個(gè)厚度單元的含水率M,j+1和介質(zhì)空氣穿過(guò)該厚度單元后的絕對(duì)濕度y+1,j+1。式（10）和式（11）可看作是上述時(shí)間節(jié)點(diǎn)和位置節(jié)點(diǎn)的物料溫度θ,j+1和介質(zhì)空氣穿過(guò)該厚度單元后的溫度T+1,j+1構(gòu)成的2個(gè)二元二次方程組，如此可計(jì)算出θ,j+1和T+1,j+1的值。

1.2.3 當(dāng)量干燥時(shí)間

干燥模擬中將筒狀固定床沿半徑向內(nèi)方向劃分了若干層，介質(zhì)空氣通過(guò)第2～層物料前均經(jīng)過(guò)前層物料的吸熱降溫，屬于變溫干燥范疇，求解含水率時(shí)不能直接將干燥時(shí)間代入式（12），需求解該層物料當(dāng)前干燥溫度和單位面積通風(fēng)量下的上一時(shí)間節(jié)點(diǎn)含水率值對(duì)應(yīng)的當(dāng)量干燥時(shí)間。

在采用MATLAB軟件對(duì)整個(gè)干燥過(guò)程進(jìn)行模擬分析時(shí)，為實(shí)現(xiàn)參數(shù)化編程采用MATLAB 自帶的優(yōu)化函數(shù)“fminsearch”求解當(dāng)量干燥時(shí)間。具體方法：令函數(shù)，其中eq為當(dāng)量干燥時(shí)間，h；M值在上一個(gè)循環(huán)中已解出；(eq)為將eq代入式（4）后的求解值。求解過(guò)程中當(dāng)函數(shù)達(dá)到最小值（逼近0）時(shí)的eq值即為M層物料在=·Δ時(shí)刻含水率值對(duì)應(yīng)的當(dāng)量干燥時(shí)間。

1.2.3 初始與邊界條件

假設(shè)干燥前，花生莢果初始干基含水率為0（無(wú)量綱）；花生莢果初始溫度為0，℃。令設(shè)定的介質(zhì)空氣通風(fēng)溫度為s，℃；絕對(duì)濕度為s(無(wú)量綱)。則有：

初始條件，

邊界條件，

1.2.4 其他輔助方程

干燥模擬中模擬精度主要取決于模型中系數(shù)和參數(shù)的取值，尤其是對(duì)流換熱系數(shù)、汽化熱、平衡含水率等重要物理參數(shù)，在花生干燥脫水過(guò)程中，伴隨莢果體積收縮和內(nèi)部水分降低不停變化，模擬計(jì)算過(guò)程中不能將之視為常量，不然將影響模擬精度。國(guó)內(nèi)外一些學(xué)者通過(guò)試驗(yàn)對(duì)這些參數(shù)進(jìn)行了詳細(xì)研究，得出了經(jīng)驗(yàn)或半經(jīng)驗(yàn)公式。

Aydin等[24]對(duì)不同水分下花生莢果體積密度進(jìn)行了測(cè)量，通過(guò)非線性回歸分析，得到了體積密度隨莢果干基含水率變化的經(jīng)驗(yàn)公式

式中db為花生莢果干基含水率（無(wú)量綱）；pb為花生莢果體積密度，kg/m3。令db=0，干質(zhì)花生體積密度216.06 kg/m3。

Wright等[25]對(duì)床層花生比熱與含水率的關(guān)系進(jìn)行了研究，得出了半經(jīng)驗(yàn)公式

式中p為床層花生比熱容，kJ/(kg·℃)。令db=0，可得干質(zhì)花生莢果比熱1.687kJ/(kg·℃)。

胡志超等[3]簡(jiǎn)化了花生莢果與介質(zhì)空氣的對(duì)流傳熱系數(shù)計(jì)算公式

式中A為介質(zhì)空氣密度（包含濕分），kg/m3；A為介質(zhì)空氣速度，m/s。

Correa等[26]對(duì)花生莢果、果殼、果仁內(nèi)壁水分汽化熱進(jìn)行研究，可通過(guò)下式對(duì)花生莢果內(nèi)部水分汽化熱求解：

式中fg*為花生莢果中水的汽化熱，kJ/kg；e為花生莢果平衡干基含水率（無(wú)量綱）。

Chen[27]對(duì)各溫度、濕度條件下的花生平衡含水率（干基）與所處的環(huán)境溫度、相對(duì)濕度進(jìn)行了研究，得出花生莢果平衡含水率半經(jīng)驗(yàn)公式

式中RH為介質(zhì)空氣相對(duì)濕度，%

2 評(píng)價(jià)指標(biāo)與模擬仿真

2.1 評(píng)價(jià)指標(biāo)

2.1.1 干燥時(shí)耗

干燥時(shí)耗為仿真程序達(dá)到停止條件的時(shí)間，即筒狀固定床環(huán)狀料層平均含水率達(dá)到終止含水率的時(shí)間。

式中co為批次干燥總消耗的時(shí)間，h；final為時(shí)間節(jié)點(diǎn)數(shù)的最后一次迭代次數(shù)。

2.1.2 批次干燥生產(chǎn)率

以干燥結(jié)束時(shí)花生總質(zhì)量與總耗時(shí)之比記為批次干燥生產(chǎn)率。

其中干燥結(jié)束時(shí)花生總質(zhì)量可通過(guò)下式計(jì)算：

式中bf為干燥結(jié)束時(shí)花生莢果的體積密度，kg/m3；final為干燥結(jié)束時(shí)花生莢果的總質(zhì)量，kg；為筒倉(cāng)高度，m。

則批次干燥生產(chǎn)率

式中為批次干燥生產(chǎn)率，kg/h。

2.1.3 單位質(zhì)量能耗

實(shí)際作業(yè)中干燥能耗主要為熱耗，風(fēng)機(jī)供風(fēng)能耗很小，可忽略。則根據(jù)進(jìn)入筒狀固定床的介質(zhì)空氣溫度、絕對(duì)濕度及風(fēng)量通過(guò)式（27）可得批次干燥總能耗。

式中0為進(jìn)入料層前的介質(zhì)空氣溫度，℃；ev為環(huán)境溫度，℃；0為介質(zhì)空氣進(jìn)入料層前的絕對(duì)濕度（無(wú)量綱）；0為介質(zhì)空氣穿過(guò)筒狀固定床外柱面的單位面積體積流量，m3/(m2·s)；co為批次干燥總能耗，J。

則單位質(zhì)量干花生平均干燥能耗：

式中pum為單位質(zhì)量干花生平均干燥能耗，J/kg。

2.1.4 干燥不均勻度

干燥不均勻度根據(jù)模擬結(jié)果獲得，取干燥結(jié)束時(shí)料層最大濕基含水率和最濕基小含水率之差，即

式中un為干燥不均勻度，%；MW(:,co)為co時(shí)刻所有料層濕基含水率的集合，%

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

為解決求解最優(yōu)干燥作業(yè)參數(shù)時(shí)，不同指標(biāo)在數(shù)學(xué)試驗(yàn)因素考查范圍內(nèi)的變化趨勢(shì)不一致，不便于數(shù)據(jù)分析的問(wèn)題，選用綜合加權(quán)評(píng)分法將多個(gè)性能指標(biāo)的試驗(yàn)結(jié)果轉(zhuǎn)化為一個(gè)單指標(biāo)的試驗(yàn)結(jié)果，利用單指標(biāo)試驗(yàn)結(jié)果進(jìn)行優(yōu)化分析。

1）確定各項(xiàng)指標(biāo)的權(quán)值

干燥時(shí)耗、生產(chǎn)率、單位質(zhì)量能耗、干燥不均勻度是考核干燥參數(shù)的重要指標(biāo)，其中干燥不均勻度是評(píng)價(jià)批次干燥作業(yè)品質(zhì)的重要指標(biāo)，能耗是反映干燥成本的重要指標(biāo)，干燥時(shí)耗和生產(chǎn)率反映了干燥能力與作業(yè)效率，根據(jù)各指標(biāo)的重要性，結(jié)合生產(chǎn)經(jīng)驗(yàn)設(shè)定時(shí)耗、生產(chǎn)率、單位質(zhì)量能耗、干燥不均勻度的權(quán)重1、2、3、4分別為0.2、0.2、0.4、0.2。

2）統(tǒng)一各項(xiàng)指標(biāo)的變化趨勢(shì)

為了保證綜合加權(quán)平均值越大越好，應(yīng)將變化趨勢(shì)越小越好的指標(biāo)值轉(zhuǎn)化為越大越好，為此在其值前加負(fù)號(hào)。時(shí)耗、單位質(zhì)量能耗、不均勻度3個(gè)指標(biāo)值越小越好，生產(chǎn)率指標(biāo)值越大越好：

=1,2,3,…,16（30）

式中1,j為第號(hào)模擬仿真干燥耗時(shí)指標(biāo)值；2,j為第號(hào)模擬仿真生產(chǎn)率指標(biāo)值；3,j為第號(hào)模擬仿真單位質(zhì)量能耗指標(biāo)值；4,j為第號(hào)模擬仿真干燥不均勻度指標(biāo)值。

3）統(tǒng)一各指標(biāo)數(shù)量級(jí)和量綱

為消除各指標(biāo)數(shù)量級(jí)和量綱對(duì)其加權(quán)評(píng)分值的影響，使各指標(biāo)的加權(quán)評(píng)分值具有可比性，需統(tǒng)一各指標(biāo)的數(shù)量級(jí)和量綱。由式（31）可得到各指標(biāo)的數(shù)量級(jí)、無(wú)量綱的評(píng)分值。

4）計(jì)算綜合加權(quán)評(píng)分

把各項(xiàng)指標(biāo)的加權(quán)評(píng)分值相加即為“綜合加權(quán)評(píng)分值”。

2.2 模擬仿真

采用MATLAB軟件對(duì)筒狀固定床花生干燥過(guò)程進(jìn)行模擬仿真，根據(jù)章節(jié)1.2和2.1的計(jì)算公式編寫(xiě)仿真程序，計(jì)算eq、、M、θ、y、T、av、co、pum、un等值，仿真程序流程圖見(jiàn)圖3。根據(jù)薄層干燥方程的適用性，設(shè)定薄層單元厚度Δ=0.03 m；同時(shí)設(shè)定時(shí)間單元Δ=0.1 h。根據(jù)行業(yè)相關(guān)標(biāo)準(zhǔn)[28]，花生莢果貯藏含水率一般8%（濕基），取之為終止含水率，即當(dāng)料層平均含水率達(dá)到8%（濕基）時(shí)仿真程序停止模擬。

根據(jù)已有經(jīng)驗(yàn)和前人研究基礎(chǔ)[4-15,29-31]，固定床花生干燥中，通風(fēng)溫度一般不宜超過(guò)38 ℃，否則會(huì)影響干燥后花生品質(zhì)，干燥模擬設(shè)定通風(fēng)溫度38 ℃；同時(shí)為簡(jiǎn)化計(jì)算，參考主產(chǎn)區(qū)花生收獲期的平均環(huán)境溫濕度及花生含水率，設(shè)定干燥模擬環(huán)境溫度25 ℃，環(huán)境相對(duì)濕度50%，初始濕基含水率40%（干基含水率0.667）。通過(guò)改變不同的單位面積通風(fēng)量和筒狀固定床外徑、內(nèi)徑對(duì)其作業(yè)性能進(jìn)行分析。

注：M0為花生莢果初始干基含水率（無(wú)量綱）；θ0為花生莢果初始溫度，℃；Ts為介質(zhì)空氣通風(fēng)溫度，℃；ys為介質(zhì)空氣絕對(duì)濕度 (無(wú)量綱)；teq為當(dāng)量干燥時(shí)間，h；m為厚度方向的最大節(jié)點(diǎn)數(shù)；Mi,j為第i個(gè)位置節(jié)點(diǎn)和第j個(gè)時(shí)間節(jié)點(diǎn)的花生干基含水率（無(wú)量綱）；θi,j為第i個(gè)位置節(jié)點(diǎn)和第j個(gè)時(shí)間節(jié)點(diǎn)的花生溫度，℃；yi,j為第i個(gè)位置節(jié)點(diǎn)和第j個(gè)時(shí)間節(jié)點(diǎn)的介質(zhì)空氣絕對(duì)濕度（無(wú)量綱）；Ti,j為第i個(gè)位置節(jié)點(diǎn)和第j個(gè)時(shí)間節(jié)點(diǎn)的介質(zhì)空氣溫度，℃；Mf為花生干燥終止含水率，%；Mav為干燥結(jié)束時(shí)料層花生的平均含水率，%。

3 干燥模型驗(yàn)證試驗(yàn)

3.1 試驗(yàn)設(shè)備與材料

為了驗(yàn)證干燥模擬結(jié)果的準(zhǔn)確性，搭建了小型筒狀固定床花生干燥試驗(yàn)臺(tái)（如圖4）進(jìn)行試驗(yàn)。

圖4 小型筒狀固定床干燥試驗(yàn)臺(tái)

該設(shè)備主要由送風(fēng)機(jī)、電加熱模塊、送風(fēng)道、勻風(fēng)機(jī)構(gòu)、筒狀床通風(fēng)干燥倉(cāng)、控制系統(tǒng)等組成。其中筒狀床料床外徑1.0 m，內(nèi)徑0.3 m，載料高度1.5 m；筒狀床外風(fēng)道直徑1.3 m，風(fēng)道內(nèi)壁粘貼橡膠海綿作為保溫材料，干燥筒載料倉(cāng)與內(nèi)、外風(fēng)道通過(guò)沖孔板隔開(kāi)，沖孔板開(kāi)孔率為22.7%。通風(fēng)溫度通過(guò)PID溫控器在環(huán)境溫度～（環(huán)境溫度+30 ℃）范圍可無(wú)級(jí)調(diào)節(jié)，總通風(fēng)量在2 000～3 560m3/h范圍通過(guò)變頻器可無(wú)級(jí)調(diào)節(jié)。試驗(yàn)選用新收獲的花生莢果，滿載，產(chǎn)地為河南正陽(yáng)縣，烘干前平均初始含水率41.2%（濕基）。烘干過(guò)程環(huán)境溫度平均值約24.1 ℃，環(huán)境相對(duì)濕度平均值約43.9%。設(shè)定通風(fēng)溫度38 ℃，通風(fēng)量3 390 m3/h，即外柱面單位面積通風(fēng)量0.2 m3/(m2s)。

3.2 試驗(yàn)過(guò)程與方法

花生莢果具有顆粒大流動(dòng)性差的特點(diǎn)，干燥過(guò)程中很難采用即插式取樣器從料床中抽取少量花生莢果以獲取筒狀床外、中、內(nèi)層含水率情況。本文根據(jù)筒狀固定床沿相同半徑尺度的圓周方向上干燥狀態(tài)一致的特點(diǎn)，為了避免集中取樣導(dǎo)致的料層坍塌影響試驗(yàn)結(jié)果，將筒狀床載料區(qū)域沿圓周方向等分8個(gè)扇區(qū)，每個(gè)扇區(qū)沿半徑方向?qū)⒘蠈臃譃橥?、中、?nèi)3個(gè)測(cè)試區(qū)域，每個(gè)測(cè)試區(qū)域料層厚度相同，合計(jì)24個(gè)測(cè)試區(qū)域，如圖5所示。試驗(yàn)開(kāi)始前在每個(gè)測(cè)試區(qū)域中間位置放置已裝填好花生莢果的自制取樣筒，試驗(yàn)開(kāi)始后每隔5 h從外到內(nèi)依次取出其中一個(gè)扇區(qū)的3個(gè)取樣筒，測(cè)量外、中、內(nèi)3層物料的含水率，其中每個(gè)取樣筒取3份物料采用烘箱法[32]測(cè)含水率并取均值作為該層物料的含水率值。上述采用的取樣筒為底部和壁面密集開(kāi)孔的圓柱形薄壁容器，介質(zhì)空氣可以從底部和側(cè)壁自由進(jìn)出容器，取樣筒直徑50 mm，高80 cm。

同時(shí)為了測(cè)量干燥過(guò)程外、中、內(nèi)層物料溫度變化情況，在第一扇區(qū)外、中、內(nèi)3個(gè)測(cè)試區(qū)域中心位置分別距離載料面30、75、120 cm各放置一個(gè)溫度傳感器，共3×3=9個(gè)。每個(gè)測(cè)試區(qū)域沿高度方向布置的3個(gè)溫度傳感器作為一組，取平均值作為該測(cè)試區(qū)域的物料層溫度值。采用的溫度傳感器為DS18B20數(shù)字溫度傳感器（測(cè)量精度±0.5 ℃），通過(guò)SM1200B-160型10通道溫度采集模塊（上海搜博實(shí)業(yè)有限公司）與計(jì)算機(jī)連接，通過(guò)上位機(jī)軟件實(shí)時(shí)獲取各測(cè)試點(diǎn)的溫度讀數(shù)。

圖5 筒狀床料層測(cè)試區(qū)域劃分

相關(guān)系數(shù)用以反映變量之間相關(guān)關(guān)系密切程度的統(tǒng)計(jì)指標(biāo)，文中采用該指標(biāo)衡量試驗(yàn)值和模型模擬值之間的差異，值越接近1則說(shuō)明模型模擬對(duì)實(shí)際干燥過(guò)程的描述越準(zhǔn)確。

4 結(jié)果與分析

4.1 模擬計(jì)算與試驗(yàn)結(jié)果

根據(jù)驗(yàn)證試驗(yàn)采用的通風(fēng)參數(shù)及試驗(yàn)時(shí)的環(huán)境條件和物料初始條件，利用模型計(jì)算料層花生烘干過(guò)程中含水率分布及變化情況，將料層花生沿筒狀床半徑方向等分為外、中、內(nèi)3層，計(jì)算每層花生平均含水率和溫度隨時(shí)間的變化情況，并將計(jì)算結(jié)果與試驗(yàn)測(cè)得值進(jìn)行比較，如圖6所示。外、中、內(nèi)3層花生含水率模擬值與試驗(yàn)測(cè)得值的相關(guān)系數(shù)分別為0.997、0.995、0.998；外、中、內(nèi)3層花生溫度模擬值與試驗(yàn)測(cè)得值的相關(guān)系數(shù)分別為0.995、0.979、0.991。試驗(yàn)對(duì)比表明，各料層花生平均含水率和溫度模擬效果良好，模型模擬可以較好地預(yù)測(cè)實(shí)際烘干生產(chǎn)中花生含水率和溫度變化情況。

4.2 干燥過(guò)程物料狀態(tài)

為了解干燥過(guò)程料層花生溫度和含水率變化情況，采用MATLAB軟件模擬了筒狀固定床外徑2.5 m，內(nèi)徑1 m，外柱面單位面積通風(fēng)量0.4 m3/(m2·s)時(shí)，沿半徑向內(nèi)方向物料層花生溫度和含水率（濕基）變化情況，結(jié)果如圖7所示。

由圖7a可知，介質(zhì)空氣由外向內(nèi)穿透物料層的過(guò)程中，與花生物料逐層發(fā)生熱質(zhì)交換，空氣溫度逐漸降低，內(nèi)層物料總體呈加熱滯后狀態(tài)（起步晚，數(shù)值低）。在干燥起步階段，最內(nèi)層物料需干燥開(kāi)始3.38 h后逐漸從初始溫度開(kāi)始升溫，其余各層初始升溫滯后時(shí)間總體符合以下線性關(guān)系

式中x為薄層花生與筒狀固定床外徑的距離，m；tTR為起始升溫滯后時(shí)間，h。

4.2.1 干燥過(guò)程物料溫度變化

在物料層快速升溫階段，外層物料快速升溫至接近干燥設(shè)定溫度，升溫時(shí)間短，升溫幅值高；而內(nèi)層物料（位置靠近0.75 m處）升溫時(shí)間稍長(zhǎng)，升溫幅值低，之后逐漸緩慢接近干燥設(shè)定溫度。該現(xiàn)象越靠近筒狀固定床內(nèi)圈越明顯，這是熱量逐步從外向內(nèi)傳遞的具體體現(xiàn)，必須等較外層的物料加熱升溫至較高的溫度，介質(zhì)空氣穿過(guò)該層物料時(shí)才能損失較少的熱量和溫度，才有足夠的熱量加熱內(nèi)層物料。截止干燥終止內(nèi)外層物料溫差約0.88 ℃。

4.2.2 干燥過(guò)程物料含水率變化

由圖7b可知，在干燥降水階段，料層花生含水率降低速率隨著干燥時(shí)間的增加逐漸降低，內(nèi)層物料較外層物料干燥降水有一定的滯后性，且干燥速率略低于外層物料。

介質(zhì)空氣由外向內(nèi)穿透物料層的過(guò)程中，不僅干燥溫度在逐漸降低，相對(duì)濕度也在逐漸增加，相應(yīng)的單位體積介質(zhì)空氣吸濕能力在逐漸降低，但由于筒狀固定床的結(jié)構(gòu)優(yōu)勢(shì)，越靠近筒狀固定床內(nèi)層的物料通風(fēng)面積越小，在總通風(fēng)量不變的情況下，單位面積通風(fēng)量越大，總體吸濕排濕能力越強(qiáng)，抵消部分因介質(zhì)空氣溫度降低和相對(duì)濕度增加而造成的干燥速率的降低。但在干燥開(kāi)始階段，物料層處于低溫高水分狀態(tài)，介質(zhì)空氣在尚未穿透筒狀固定床時(shí)溫度就已經(jīng)降低至接近物料層初始溫度，且相對(duì)濕度已接近飽和或達(dá)到與物料層構(gòu)成的平衡相對(duì)濕度，不再具備吸濕能力，因此內(nèi)層物料相對(duì)于外層物料含水率開(kāi)始降低的時(shí)間有一定的滯后，滯后時(shí)間與料層位置符合以下線性關(guān)系

式中MR為含水率下降滯后時(shí)間，h。

令上式中=0.75可得最內(nèi)層物料含水率降低滯后時(shí)間約為2.04 h，比章節(jié)4.2.1中提到的最內(nèi)層物料受熱升溫滯后時(shí)間提前了1.34 h，這主要是由于干燥至2.04 h后介質(zhì)空氣穿過(guò)最內(nèi)層物料時(shí)的相對(duì)濕度已低于此時(shí)與花生物料構(gòu)成的平衡相對(duì)濕度，具備一定的吸濕排濕能力。截止干燥終止內(nèi)外層物料含水率差值約0.87%。

注：筒狀固定床外徑2.5 m，內(nèi)徑1 m，外柱面單位面積通風(fēng)量0.4 m3·m-2·s-1。

4.3 通風(fēng)和結(jié)構(gòu)參數(shù)對(duì)干燥指標(biāo)的影響

筒狀固定床通風(fēng)參數(shù)主要包括通風(fēng)溫度和外柱面單位面積通風(fēng)量，結(jié)構(gòu)參數(shù)主要包括外徑、內(nèi)徑及載料高度。章節(jié)2.2已提及鑒于已有研究基礎(chǔ)和收獲后花生含水率概況，模擬仿真中通風(fēng)溫度取38 ℃，花生物料初始含水率取40%（濕基）。此外，根據(jù)式（1）～（7）及（25）～（29）可知，載料高度與生產(chǎn)率成正比關(guān)系，與干燥時(shí)耗、單位質(zhì)量能耗、干燥不均勻度無(wú)關(guān)，而外柱面單位面積通風(fēng)量、外徑、內(nèi)徑與干燥時(shí)耗、生產(chǎn)率、單位質(zhì)量能耗、干燥不均勻度4個(gè)干燥指標(biāo)之間均存在復(fù)雜的非線性關(guān)系。因此，為了解它們之間的影響規(guī)律，需逐一分析。

4.3.1 通風(fēng)量對(duì)干燥指標(biāo)的影響

為了解不同通風(fēng)量對(duì)干燥指標(biāo)的影響，設(shè)定筒狀固定床外徑2.5 m，內(nèi)徑1 m，載料高度3 m，采用MATLAB軟件模擬了外柱面單位面積通風(fēng)量0.1～0.5 m3/(m2·s)范圍時(shí)的干燥過(guò)程，計(jì)算了干燥時(shí)耗、不均勻度、生產(chǎn)率、單位質(zhì)量能耗4個(gè)指標(biāo)值，并繪制了該4個(gè)指標(biāo)值隨單位面積通風(fēng)量變化曲線，如圖8所示。

注：筒狀固定床外徑2.5 m，內(nèi)徑1 m，載料高度3 m。

Note: The outer diameter of the barrel-shaped fixed bed is 2.5 m, the internal diameter is 1 m, and the loading height is 3 m.

圖8 不同通風(fēng)量下的各項(xiàng)干燥指標(biāo)

Fig.8 Drying indicators under different ventilation volume

由圖8可知，干燥時(shí)耗和不均勻度隨著單位面積通風(fēng)量的增加逐漸降低，其中0.1～0.3 m3/(m2·s)范圍降低幅度較大，為17.7 h和3.36%，分別為總降幅的75.3%、83.6%。生產(chǎn)率和單位質(zhì)量能耗隨通風(fēng)量的增大而增加，其中生產(chǎn)率在風(fēng)量為0.1～0.3 m3/(m2·s)范圍增幅稍高，約20 kg/h，為總增幅的64.7%。單位質(zhì)量能耗和單位面積通風(fēng)量總體呈近似線性增加的關(guān)系，單位面積通風(fēng)量0.5 m3/(m2·s)時(shí)的單位質(zhì)量能耗是0.1 m3/(m2·s)時(shí)的2.99倍，增幅約5.56×106J/kg。因此從節(jié)本增效、限制干燥不均勻度等因素綜合考慮，參考行業(yè)標(biāo)準(zhǔn)NYT 2785[28]（批次干燥降水幅度≥20%，干燥不均勻度≤3%），筒狀固定床外徑2.5 m內(nèi)徑1 m時(shí)，設(shè)定外柱面單位面積通風(fēng)量0.3 m3/(m2·s)為宜，此時(shí)干燥時(shí)耗40.5 h，不均勻度1.35%，生產(chǎn)率65.8 kg/h，單位質(zhì)量能耗5.85×106J/kg。

4.3.2 筒狀固定床外徑對(duì)干燥指標(biāo)的影響

采用MATLAB軟件模擬了筒狀固定床內(nèi)徑1 m，外柱面單位面積通風(fēng)量0.3 m3/(m2·s)，載料高度3 m時(shí)，外徑2～3 m范圍時(shí)的干燥過(guò)程，計(jì)算了干燥時(shí)耗、不均勻度、生產(chǎn)率、單位質(zhì)量能耗4個(gè)指標(biāo)值，并繪制了該4個(gè)指標(biāo)值隨外徑變化曲線，如圖9所示。

注：筒狀固定床內(nèi)徑1 m，外柱面單位面積通風(fēng)量0.3 m3·m-2·s-1，載料高度3 m。

由圖9可知，干燥時(shí)耗、不均勻度、生產(chǎn)率均隨外徑的增大而增加，單位質(zhì)量能耗隨外徑的增加而降低。這是由于隨著外徑的增加，料層厚度、筒狀床外柱面面積和裝載量相應(yīng)增加。料層厚度增加導(dǎo)致干燥時(shí)間的延長(zhǎng)和干燥不均勻度的增高，而裝載量隨著外徑的2次方增長(zhǎng)（增長(zhǎng)率平均166.7%/m），遠(yuǎn)高于干燥時(shí)耗增長(zhǎng)率（平均8.72%/m），因此生產(chǎn)率隨著外徑增大而增加（增長(zhǎng)率平均145.2%/m）。同理，外柱面通風(fēng)面積和干燥總時(shí)耗增加而導(dǎo)致的干燥總能耗增長(zhǎng)率（平均63.1%/m）小于裝載量的增長(zhǎng)率，因此均攤至單位質(zhì)量的能耗成本隨外徑增加而降低（降幅率平均38.78%/m）。因此從節(jié)本增效、限制干燥不均勻度等因素綜合考慮，參考行業(yè)標(biāo)準(zhǔn)NYT 2785[28]（批次干燥降水幅度≥20%，干燥不均勻度≤3%），筒狀固定床內(nèi)徑1m，外柱面單位面積通風(fēng)量0.3 m3/(m2·s)時(shí)，設(shè)定外徑2.6 m為宜，此時(shí)干燥時(shí)耗40.9 h，不均勻度1.49%，生產(chǎn)率71.4 kg/h，單位質(zhì)量能耗5.60×106J/kg。

4.3.3 筒狀固定床內(nèi)徑對(duì)干燥指標(biāo)的影響

采用MATLAB軟件模擬了筒狀固定床外徑2.5 m，裝載高度3 m，外柱面單位面積通風(fēng)量0.3 m3/(m2·s)時(shí)，內(nèi)徑0.6～1.4 m范圍時(shí)的干燥過(guò)程，計(jì)算了干燥時(shí)耗、不均勻度、生產(chǎn)率、單位質(zhì)量能耗4個(gè)指標(biāo)值，并繪制了該4個(gè)指標(biāo)隨外徑變化曲線，如圖10所示。

注：筒狀固定床外徑2.5 m，裝載高度3 m，外柱面單位面積通風(fēng)量0.3 m3·m-2·s-1。

隨著內(nèi)徑增加，料層厚度和裝載量相應(yīng)降低，干燥時(shí)耗、不均勻度逐漸降低，單位質(zhì)量能耗逐漸升高；此外，裝載量降低比率大于干燥時(shí)耗降低比率，因此生產(chǎn)率均逐漸降低。其中，干燥時(shí)耗降低速率較穩(wěn)定；不均勻度降低速率逐漸減小，在0.6～0.9 m內(nèi)徑范圍降幅較明顯，占總降幅的75.4%；生產(chǎn)率降低速率逐漸增大，在1.0～1.4 m內(nèi)徑范圍降幅相對(duì)明顯，占總降幅的63.6%；單位質(zhì)量能耗增加速率逐漸增大，1.0～1.4 m內(nèi)徑范圍增幅相對(duì)明顯，占總增幅的69.9%。因此從節(jié)本增效、限制干燥不均勻度等因素綜合考慮，參考行業(yè)標(biāo)準(zhǔn)NYT 2785[28]（批次干燥降水幅度≥20%，干燥不均勻度≤3%），筒狀固定床外徑2.5 m，外柱面單位面積通風(fēng)量0.3 m3/(m2·s)時(shí)，設(shè)定內(nèi)徑0.9 m為宜，此時(shí)干燥時(shí)耗40.8 h，不均勻度1.42%，生產(chǎn)率67.64 kg/h，單位質(zhì)量能耗5.69×106J/kg。

4.4 結(jié)構(gòu)與通風(fēng)參數(shù)優(yōu)化

4.4.1 均勻設(shè)計(jì)模擬試驗(yàn)與結(jié)果

為獲得最優(yōu)的結(jié)構(gòu)和通風(fēng)參數(shù)，參考上文通風(fēng)量、外徑、內(nèi)徑對(duì)4個(gè)干燥指標(biāo)影響的單因素模擬試驗(yàn)分析結(jié)果，進(jìn)一步縮小單因素模擬試驗(yàn)中各因素的取值范圍，選取單位面積通風(fēng)量0.24～0.36 m3/(m2·s)，外徑2.45～2.75 m，內(nèi)徑0.75～1.05 m；采用U13*(134)均勻設(shè)計(jì)表，開(kāi)展了以干燥時(shí)耗、生產(chǎn)率、單位質(zhì)量能耗、不均勻度為指標(biāo)的三因素13水平均勻設(shè)計(jì)模擬試驗(yàn)，在此基礎(chǔ)上計(jì)算了綜合加權(quán)評(píng)分值。模擬試驗(yàn)安排和結(jié)果見(jiàn)表1。

根據(jù)模擬試驗(yàn)結(jié)果，采用MATLAB軟件“STEPWISE”命令進(jìn)行三元二次逐步回歸分析，剔除不顯著項(xiàng)，得到綜合加權(quán)評(píng)分值與單位面積通風(fēng)量、外徑、內(nèi)徑的回歸模型，如式（35）。決定系數(shù)2、剩余標(biāo)準(zhǔn)差RMSE、值、值等模型質(zhì)量評(píng)價(jià)指標(biāo)見(jiàn)表2，各模型顯著項(xiàng)值見(jiàn)表3。

式中*為綜合加權(quán)評(píng)分值；1為單位面積通風(fēng)量，m3/(m2·s)；2為外徑，m；3為內(nèi)徑，m。

表1 U13*（134）均勻設(shè)計(jì)試驗(yàn)安排和模擬計(jì)算結(jié)果

表2 回歸模型質(zhì)量評(píng)價(jià)

表3 回歸模型各顯著項(xiàng)P值

4.4.2 最優(yōu)通風(fēng)參數(shù)的求解

為便于綜合加權(quán)評(píng)分*最大值的求解，以*相反數(shù)的最小值為目標(biāo)函數(shù)，以單位面積通風(fēng)量、外經(jīng)、內(nèi)徑3個(gè)影響因素邊界范圍為約束條件，則有：

目標(biāo)函數(shù)：

約束條件：

采用MATLAB軟件“fmincon”函數(shù)，求解上述多元線性約束優(yōu)化方程?？傻茫?=0.36 m3/(m2·s)，2=2.75 m，3=0.935 m時(shí)值最小?72.428，即獲得最大綜合加權(quán)評(píng)分值72.428，優(yōu)于表1中的所有試驗(yàn)因素組合。受各評(píng)價(jià)指標(biāo)權(quán)重分配的影響，1和2均是約束條件內(nèi)的最大值。將1、23值作為輸入值，模擬仿真筒狀固定床干燥可得干燥耗時(shí)39.2 h、生產(chǎn)率86.55 kg/h、單位質(zhì)量能耗5.87×106J/kg、干燥不均勻度1.54%。

在小型臺(tái)架試驗(yàn)驗(yàn)證、結(jié)構(gòu)與工藝參數(shù)優(yōu)化分析基礎(chǔ)上，研發(fā)適用于實(shí)際生產(chǎn)中的大、中型筒狀固定床花生干燥機(jī)械設(shè)備將是今后進(jìn)一步研究工作。

5 結(jié) 論

1）模型模擬結(jié)果與試驗(yàn)結(jié)果基本相符，外、中、內(nèi)各層花生平均含水率和溫度模擬值和試驗(yàn)測(cè)得值的相關(guān)系數(shù)均達(dá)0.975以上，模型模擬可以準(zhǔn)確描述花生實(shí)際烘干過(guò)程中含水率和溫度變化情況。

2）介質(zhì)空氣在沿半徑方向從外向內(nèi)穿過(guò)筒狀固定床對(duì)花生物料進(jìn)行干燥時(shí)，由于溫度的降低和相對(duì)濕度的增加，內(nèi)層物料在干燥起始時(shí)間和干燥速率方面均存在一定的滯后，但單位面積通風(fēng)量沿通風(fēng)方向逐漸增大，對(duì)內(nèi)層物料的干燥滯后問(wèn)題有著較好的改善。

3）隨著單位面積通風(fēng)量的增加，筒狀固定床干燥不均勻度有明顯降低，生產(chǎn)率亦有顯著提高，但單位質(zhì)量能耗增幅較大；筒狀床外徑增加或內(nèi)徑減小都可增加生產(chǎn)率，降低單位質(zhì)量干燥能耗，但干燥均勻性問(wèn)題很難解決。

4）當(dāng)外進(jìn)風(fēng)面單位面積通風(fēng)量0.36 m3/(m2·s)、筒狀固定床外徑2.75 m、內(nèi)徑0.935 m時(shí)，干燥綜合評(píng)分指標(biāo)最優(yōu)。此時(shí)綜合加權(quán)評(píng)分值72.428，干燥耗時(shí)39.2 h、生產(chǎn)率86.55 kg/h、單位質(zhì)量能耗5.87×106J/kg、干燥不均勻度1.54%。

[1]高連興，陳中玉，Charles Chen，等. 美國(guó)花生收獲機(jī)械化技術(shù)演變歷程及對(duì)中國(guó)的啟示[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(12)：1－9. Gao Lianxing, Chen Zhongyu, Charles Chen, et al. Development course of peanut harvest mechanization technology of the United States and enlightenment to China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 1－9. (in Chinese with English abstract)

[2]王海鷗，胡志超，陳守江，等. 收獲時(shí)期及干燥方式對(duì)花生品質(zhì)的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(22)：292－300. Wang Haiou, Hu Zhichao, Chen Shoujiang, et al. Effects of different harvesting dates and drying methods on peanut quality[J]. Transactions of the Chinese Society of AgriculturalEngineering (Transactions of the CSAE), 2017, 33(22): 292－300. (in Chinese with English abstract)

[3]胡志超. 花生生產(chǎn)機(jī)械化關(guān)鍵技術(shù)[M]. 鎮(zhèn)江：江蘇大學(xué)出版社，2017.

[4]Cundiff J S, Baker K D. Curing Quality Peanuts in Virginal[M]. Virgina: Virgina Tech. 2009.

[5]Jordan D, Brandenburg R, Brown A, et al. 2018 Peanut Information[M]. North Carolina State: North Carolina Cooperative Extension Service, 2018.

[6]Butts C L, Davidson J I, Lamb M C, et al. Estimating drying time for a stock peanut curing decision support system[J]. American Society of Agricultural Engineers, 2004. 47(3): 25－932.

[7]Palacios T R, Potes L B, Montenegro R A, et al. Peanut drying kinetics: Determination of the effective diffusivity for in-shell and shelled peanuts by applying a short-time analytical model of measured data[J]. 14th International Drying Symposium, 2004, 8(B): 1448－1455.

[8]Butts C L, Williams E J, Sanders T H. Algorithms for automated temperature controls to cure peanuts[J]. Postharvest Biology and Technology, 2002, 24(6): 309－316.

[9]Yang C Y, Fon D S, Lin T T. Simulation and validation of thin layer models for peanut drying[J]. Drying Technology, 2007, 25(9): 1515－1526.

[10]Krzyzanowski F C, West S H, Barros J D. Drying peanut seed using air ambient temperature allow relative humility[J]. Revista Brasileira de Sementes. 2006, 28(3): 1－5.

[11]Nakai V K, Rocha L O, Goncalez E, et al. Distribution of fungi and aflatoxins in a stored peanut variety[J]. Food Chemistry, 2008, 106(1): 285－290.

[12]Yu J, Ahmedna M, Goktepe I. Peanut protein concentrate: Production and functional properties as affected by processing[J]. Food Chemistry, 2007, 103(1): 121－129.

[13]Butts C L, Dorner J W, Brown S L, et al. Aerating farmer stock peanut storage in the southeastern US[J]. Transactions of the ASABE, 2006, 49(2): 457－465.

[14]Butts C L, Williams E J. Measuring airflow distribution in peanut drying trailers[J]. Applied Engineering in Agriculture, 2004, 20(3): 335－339.

[15]顏建春，謝煥雄，魏海，等. 5H-1.5A型花生換向通風(fēng)干燥機(jī)研制[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(10)：9－18.

Yan Jianchun, Xie Huanxiong, Wei Hai, et al. Development of 5H-1.5A peanut reversing ventilation dryer[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(10): 9－18. (in Chinese with English abstract)

[16]顏建春，吳努，胡志超，等. 花生干燥技術(shù)概況與發(fā)展[J].中國(guó)農(nóng)機(jī)化，2012，32(2)：10－13. Yan Jianchun, Wu Nu, Hu Zhichao, et al. Overview and development of peanut drying technology[J]. Chinese Agricultural Mechanization, 2012, 32(2): 10－13. (in Chinese with English abstract)

[17]劉麗，王強(qiáng)，劉紅芝. 花生干燥貯藏方法的應(yīng)用及研究現(xiàn)狀[J]. 農(nóng)產(chǎn)品加工，2011(8)：49－52. Liu Li, Wang Qiang, Liu Hongzhi. Application and its present on method for drying storage of peanut[J]. Farm Products Processing, 2011(8): 49－52. (in Chinese with English abstract)

[18]王安建，高帥平，田廣瑞，等. 花生熱泵干燥特性及動(dòng)力學(xué)模型[J]. 農(nóng)產(chǎn)品加工，2015(5)：49－52.Wang Anjian, Gao Shuaiping, Tian Guangrui, et al. Hecatomb drying characteristics and dynamics model of peanut[J]. Farm Products Processing, 2015(5): 49－52. (in Chinese with English abstract)

[19]潘永康，王喜忠，劉相東. 現(xiàn)代干燥技術(shù)[M]. 北京：化學(xué)工業(yè)出版社，2006.

[20]Zarea D, Chen G N. Evaluation of a simulation model in predicting the drying parameters for deep-bed paddy drying[J]. Computers and Electronics in Agriculture, 2009, 68(1): 78－87.

[21]Aregba A W, Nadeau J P. Comparison of two non-equilibrium models for static grain deep-bed drying by numerical simulations[J]. Journal of Food Engineering, 2007, 78(4): 1174－1187.

[22]Mujumdar A S. Handbook of industrial drying[M]. USA: CRC Press Inc, 2014.

[23]顏建春，謝煥雄，胡志超，等. 固定床上下?lián)Q向通風(fēng)小麥干燥模擬與工藝優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31(22)：292－300.

Yan Jianchun, Xie huanxiong, Hu Zhichao, et al. Simulation and process optimization of upward and downward reversing ventilating drying by fixed bed[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CASE), 2015, 31(22): 292－300. (in Chinese with English abstract)

[24]Aydin C. Some engineering properties of peanut and kernel[J]. Journal of Food Engineering, 2009, 79: 810－816.

[25]Wright M E, Porterfield J G. Specific heat of Spanish peanuts[J]. Transaction of the ASAE, 1970, 13(4): 508－510.

[26]Correa P C, Goneli A L, Jaren C，et al. Sorption isotherms and isosteric heat of peanut pods, kernels and hulls[J]. Food Science and Technology International, 2007,13(3): 231－238.

[27]Chen C C. A rapid method to determine the sorption isotherms of peanuts[J]. Journal of Agriculture Engineering Research, 2000, 75: 401－408.

[28]花生熱風(fēng)干燥技術(shù)規(guī)范：NY/T 2785-2015 [S]. 北京：中國(guó)標(biāo)準(zhǔn)出版社，2015.

[29]Baker K D, Cundiff J S, Wright F S. Peanut quality improvement through controlled curing[J]. Peanut Science, 1993, 20(12): 12－16.

[30]Butts C L. Comparison of peanut dryer control strategies[J]. Peanut Science, 1996, 23(14): 86－90.

[31]Wright F S, Desk S H, Cundiff J S. Storing peanut in trailer-sized containers[J]. Peanut Science, 1996, 23(1): 43－45.

Performance index simulation and analysis of peanut ventilation drying in barrel-shaped fixed bed

Yan Jianchun, Wei Hai, Xie Huanxiong※, You Zhaoyan

(1.,210014; 2.,210014)

Mechanical drying become gradually used to reduce post-harvest losses of peanuts because the conventional sun-drying facilities cannot meet the requirements of peanuts drying in China. A barrel-shaped fixed bed drying, a kind of economical and practical dry method, which can be effectively suitable for the specific conditions in Chinese rural production.However,there are few reported on the parameters of structure and ventilation associated with the barrel-shaped fixed bed.

In order to obtain the drying mechanism and operating performance, a set of partial differential equation models were proposed for peanut drying in the barrel-shaped fixed bed, indicating the heat- and mass transfer relationships between peanut pods and medium air during the processing of mechanical drying. Then, a fully discretized physical model was established for the computer simulation.This model could also be used to calculate the moisture content and temperature of the peanut material layers, as well as the humidity and temperature of air medium in the real-time state during the peanut drying, specifically including the time consumption, unevenness distribution, productivity and energy consumption in the per unit mass during operation simulation for each running of peanut drying. The simulation results were in good agreement with the experimental ones. Between test- and simulation values, the correlation coefficients of the average moisture content for the outer, middle and inner of peanut material layers were 0.997, 0.995 and 0.998, respectively, whereas the correlation coefficients of the temperature for the outer, middle and inner of peanut material layers were 0.995, 0.979 and 0.991, respectively. The changing mechanism was therefore established to analyze the moisture and temperature distribution on the peanut-pods bed during drying process. The ventilation volumes per unit area of the outer cylindrical surface, external diameter and inner diameter of the barrel-shaped fixed bed were also analyzed according the previous four drying indexes(drying time consumption, productivity, energy consumption per unit mass, drying unevenness). The calculated results showed that the starting time of drying and drying rate of the inner materials layers were lagged due to the decrease in the temperature of air medium, while the increase in the relative humidity, compared with the outer layers in the peanut drying.However, the drying delay of the inner materials was well relieved as the increase of the ventilation volume per unit area along the ventilating direction, whereas the drying evenness of the barrel-shaped fixed bed decreased significantly as the increase of the ventilation volume per unit area of the outer cylindrical surface, but the productivity and the energy consumption per unit mass increased significantly. The productivity of peanuts would be increased, while the energy consumption per unit mass of the batch drying would be reduced as the increase in the external diameter or decrease in the inner diameter of the barrel-shaped fixed bed, but the distribution of the peanut moisture contents after drying operation was still remained unclear during this time. Therefore, the uniform design simulation test and synthetical weighted mark method were used to obtain the optimal structure- and ventilation parameters. The simulated results demonstrated that the optimal drying condition was achieved, where the external diameter of the barrel-shaped fixed bed was 2.75m, the internal diameter was 0.935m, and the ventilation volume per unit area of the outer cylindrical surface was 0.36 m3/(m2s). At this time, the drying time consumption was 39.2 h, the productivity was 86.55 kg/h, the energy consumption per unit mass was 5.87× 106J/kg, and the drying unevenness was 1.54%. These findings can therefore provide systematically technical support to the design and optimization of the barrel-shaped fixed bed for peanuts drying.

agriculture products; drying; peanut; fixed box; simulation

顏建春，魏海，謝煥雄，游兆延. 筒狀固定床花生通風(fēng)干燥性能指標(biāo)模擬與分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(1)：292－302.doi：10.11975/j.issn.1002-6819.2020.01.035 http://www.tcsae.org

Yan Jianchun, Wei Hai, Xie Huanxiong, You Zhaoyan. Performance index simulation and analysis of peanut ventilation drying in barrel-shaped fixed bed[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(1): 292－302. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.01.035 http://www.tcsae.org

2019-09-04

2019-12-30

江蘇省農(nóng)機(jī)三新工程（NJ2018-10）；中國(guó)農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程（農(nóng)產(chǎn)品分級(jí)與貯藏裝備創(chuàng)新團(tuán)隊(duì)）

顏建春，助理研究員，主要從事農(nóng)產(chǎn)品加工技術(shù)與裝備研究。Email：yan.jc@Foxmail.com

謝煥雄，研究員，主要從事農(nóng)產(chǎn)品加工技術(shù)與裝備研究。Email：nfzhongzi@163.com

10.11975/j.issn.1002-6819.2020.01.035

S226.6

A

1002-6819(2020)-01-0292-11