郭守志，宋占華，2，李法德，2，宋華魯，張寧，閆銀發(fā)，2*

1.山東農(nóng)業(yè)大學(xué)機(jī)電學(xué)院，山東泰安271018

2.山東省園藝機(jī)械與裝備重點實驗室，山東泰安271018

基于Weibull分布函數(shù)的檳榔干燥模擬

郭守志1，宋占華1，2，李法德1，2，宋華魯1，張寧1，閆銀發(fā)1，2*

1.山東農(nóng)業(yè)大學(xué)機(jī)電學(xué)院，山東泰安271018

2.山東省園藝機(jī)械與裝備重點實驗室，山東泰安271018

本文利用Weibull分布函數(shù)描述檳榔在不同干燥方式（烘箱干燥和熱風(fēng)干燥）、不同干燥溫度（60、80和100℃）下的干燥過程，并模擬和分析了其干燥動力曲線。結(jié)果表明，Weibull分布函數(shù)能夠很好地模擬檳榔在試驗條件下的干燥過程；尺度參數(shù)α隨著干燥溫度的升高而減??；尺度參數(shù)β對干燥溫度影響很小。干燥過程中的水分?jǐn)U散系數(shù)Dcal分布在0.1123×10-7～0.4105×10-7m2·s-1之間，在烘箱干燥和熱風(fēng)干燥方式下的干燥活化能分別為33.19和32.55 kJ/mol。

干燥；Weibull分布函數(shù)；檳榔；模擬

由于檳榔種植的區(qū)域性，其果實的干制為檳榔的運輸和貯藏提供了方便。我國南方地區(qū)，檳榔的干燥方式通常是自然曬干或熏干，大約需要一周左右的時間，工藝周期較長，生產(chǎn)效率低［1］。因為傳熱傳質(zhì)效率、能源消耗和產(chǎn)品品質(zhì)等是檳榔干燥的過程中的重要指標(biāo)，所以對干燥過程的有效預(yù)測、控制和優(yōu)化具有十分重要的意義。一些經(jīng)典或半經(jīng)典模型，如Henderson-Pabis，Page，Modified Page，Logarithmic模型，近似擴(kuò)散模型等，被用來擬合干燥動力曲線［2，3，4］。盡管一些模型能夠很好的描述干燥過程，但干燥過程還是被認(rèn)為處在“黑箱”階段，干燥物料和干燥條件等很難與上述模型中的各參數(shù)關(guān)聯(lián)起來［5］。

Weibull分布函數(shù)以其較好的適用性和兼容性被廣泛用于材料科學(xué)、制藥學(xué)以及熱力學(xué)等領(lǐng)域［6-9］。近年來一些研究學(xué)者把Weibull分布函數(shù)逐步應(yīng)用于干燥動力學(xué)方面，并取得一些進(jìn)展。Corzo［10］等用Weibull分布函數(shù)描述不同成熟期的芒果薄片的干燥過程；Bantle［11］等研究了物料在常壓冷凍干燥方式下，不同溫度、風(fēng)速和粒徑對Weibull函數(shù)各參數(shù)的影響。但就總體而言，Weibull分布函數(shù)各參數(shù)的物理意義和影響因素有待進(jìn)一步的深入研究。

本文主要研究了不同干燥方法和干燥溫度對檳榔干燥動力學(xué)的影響，利用Weibull分布函數(shù)建立了其干燥過程的模型，并探究和闡述了模型中各參數(shù)的影響因素以及物理意義。為檳榔的干燥加工的預(yù)測和調(diào)控提供了依據(jù)，也為Weibull分布函數(shù)在干燥領(lǐng)域的運用提供一些理論支持。

1　材料與方法

1.1試驗材料

試驗所用檳榔由湖南胖哥食品有限公司提供，在其原產(chǎn)地海南省經(jīng)過坑式干燥處理。樣品平均長度和直徑分別為53±5 mm和20±2 mm，平均質(zhì)量為8.37±1.7 g。檳榔初始含水率為36%±0.2%（由烘干法測定）。試驗前，檳榔放在4±1℃的冰箱中冷藏。

1.2試驗方法

盡量挑選尺寸相當(dāng)?shù)臋壚谱鳛閷嶒灢牧?。本文主要?yīng)用兩種干燥方式：烘箱干燥（202-2型電熱干燥箱，上海第二五金廠）和熱風(fēng)干燥（自制熱風(fēng)干燥箱，主要由三部分構(gòu)成——風(fēng)機(jī)、熱交換器和干燥室。風(fēng)機(jī)提供干空氣的流動，熱交換器由電加熱管和溫控儀共同作用來調(diào)控溫度，干燥室用來放置被干燥物料。）干燥溫度分別設(shè)定為60、80和100℃，風(fēng)速恒定為1 m/s。從干燥開始起利用電子稱（G&G JJ1000型電子稱，常熟雙杰實驗儀器廠）測量檳榔的質(zhì)量，之后每隔1 h稱重1次，直到檳榔含水率降至15%時停止試驗，每組試驗重復(fù)3次。

1.3Weibull分布函數(shù)模型

干燥過程中檳榔的水分比（MR，Moisture Ratio）可以通過式（1）計算［4，12］：

式中，M0，Mt和Me分別表示檳榔在干燥起始、任意時刻t以及達(dá)到平衡時的干基含水率，g/g。干燥速率（DR，Drying Rate）可由式（2）表示：

式中，Mt1和Mt2分別表示檳榔在干燥過程中t1和t2時刻的干基含水率，g/g。Weibull分布函數(shù)由式（3）表示［4，9，10］：

式中，尺度參數(shù)α和形狀參數(shù)β決定了檳榔在任意t時刻時的水分比MR。尺度參數(shù)α表示干燥的速率常數(shù)，min，其值大約等于干燥過程完成63%時所用的時間。形狀參數(shù)β與干燥初期物料的傳質(zhì)速率有關(guān)，β值越小，說明干燥初期的干燥速率越大［3，4，13］。

用非線性回歸分析對試驗數(shù)據(jù)進(jìn)行Weibull分布函數(shù)模型的擬合。模型擬合優(yōu)劣的指標(biāo)主要是由相關(guān)系數(shù)R2、均方根誤差RMSE以及離差平方和x2進(jìn)行評估，R2值越大、RMSE和x2值越小，說明擬合的越好，其計算公式如下［1，2，8］。模型擬合過程是由分析軟件Origin 9.0完成。

式中，MRexp和MRpre分別為試驗實測值和模型預(yù)測值；N為試驗測得的數(shù)據(jù)個數(shù)。

2　試驗結(jié)果

2.1干燥動力曲線2.1.1烘箱干燥曲線不同干燥溫度下檳榔的烘箱干燥曲線如圖1-a所示。檳榔干燥至目標(biāo)含水率所用的時間隨著干燥溫度的增大而減少。提升干燥溫度可以明顯縮短干燥時間，提高效率。干燥速率曲線如圖1-b所示，干燥過程中無恒速干燥階段，并且隨著干燥的進(jìn)行，干燥速率基本呈下降趨勢。

圖1　不同干燥溫度下檳榔的烘箱干燥曲線Fig.1 Curves of drying Areca nut with dry-chamber at different temperatures

2.1.2熱風(fēng)干燥曲線不同干燥溫度下檳榔的熱風(fēng)干燥曲線如圖2-a所示。同烘箱干燥一樣，達(dá)到目標(biāo)含水率所用的時間隨著干燥溫度的增大而降低。干燥速率曲線如圖2-b所示，從圖中明顯可以看出整個干燥過程都主要處于降速階段，并且也未出現(xiàn)恒速干燥階段。盡管如此，從圖1和圖2的結(jié)果可以看出，烘箱干燥和熱風(fēng)干燥，干燥到目標(biāo)含水率所用的時間和干燥速率，幾乎沒有顯著差異。

圖2　不同干燥溫度下的檳榔熱風(fēng)干燥曲線Fig.2 Curves of drying Areca nut with hot air at different temperatures

2.2干燥曲線的Weibull分布函數(shù)模型

2.2.1尺度參數(shù)α的物理意義以及影響因素運用Weibull分布函數(shù)去擬合不同干燥條件下的檳榔干燥曲線，擬合結(jié)果見下表1。從表中可以看出，相關(guān)系數(shù)R2在0.993～0.999區(qū)間內(nèi)，均方根誤差RMSE在3.35×10-3～1.995×10-2之間，離差平方和χ2在1.12×10-5～3.97×10-4之間。由擬合結(jié)果可知，Weibull分布函數(shù)能很好的擬合檳榔的干燥曲線，也為干燥過程的分析提供基礎(chǔ)條件。

Weibull分布函數(shù)中尺度參數(shù)α，表示干燥過程的速率常數(shù)，其值大約為干燥過程完成63%所需要的時間（以min表示）。從表1可以看出，隨著干燥溫度的提升（60，80，100℃），經(jīng)烘箱干燥的檳榔尺度參數(shù)α由593.259 min減小到164.644 min，而經(jīng)熱風(fēng)干燥的檳榔尺度參數(shù)α由565.5 min減小到162.384 min。盡管相同的干燥溫度條件下，不同的干燥方式對應(yīng)的尺度參數(shù)α也不相同。

這就說明，對同一干燥物料來說，尺度參數(shù)α與干燥方式和干燥溫度有關(guān)。在相同干燥方式下，尺度參數(shù)α隨著干燥溫度的升高而減小。

2.2.2形狀參數(shù)β的物理意義以及影響因素相關(guān)研究表明［3，8，11］，形狀參數(shù)β與干燥過程中水分的遷移機(jī)理有關(guān)，即當(dāng)形狀參數(shù)β在0.3～1時，表示物料在干燥過程中受內(nèi)部水分?jǐn)U散控制，主要表現(xiàn)為降速干燥的特點；而當(dāng)形狀參數(shù)β大于1時，干燥曲線呈現(xiàn)“Z”字形，即物料在干燥初期有延滯階段，表現(xiàn)為干燥速率在干燥前期有先升高后降低的變化趨勢。因此當(dāng)描述某種干燥狀態(tài)時，可根據(jù)形狀參數(shù)β的值和1的大小關(guān)系來判斷干燥過程的進(jìn)行［4，14］。

表1　Weibull分布函數(shù)對檳榔干燥曲線的擬合結(jié)果Table 1 Results of Areca nut drying curves modeled by Weibull distribution function

2.3檳榔水分?jǐn)U散系數(shù)的計算

干燥過程中的水分的遷移是由多種因素共同作用的結(jié)果，包括毛細(xì)管流動、液體擴(kuò)散、蒸汽擴(kuò)散等，故通常被認(rèn)為是一個很復(fù)雜的過程。為綜合考量這些影響因素，本試驗測量并計算了干燥過程中的水分有效擴(kuò)散系數(shù)Deff。這對干燥過程的描述和優(yōu)化有十分重要的意義。

干燥過程中水分有效擴(kuò)散系數(shù)的計算通常是采用Fick第二定律。但Fick第二定律有其局限性，因為它成立的前提條件是整個干燥過程要始終處于降速干燥過程，這就大大限制了其適用范圍。然后Weibull分布函數(shù)在干燥過程上的另一重要應(yīng)用，就是它可以在不考慮水分遷移特點的情況下，估算出干燥過程中的水分?jǐn)U散系數(shù)Dcal?？梢酝ㄟ^下面的公式進(jìn)行計算［6，8，15］：

式中，Dcal是干燥過程中估算的水分?jǐn)U散系數(shù)，m2/s；r是檳榔的等效半徑，其值大約為0.02 m；α是Weibull分布函數(shù)中的尺度參數(shù)，s。

計算水分?jǐn)U算系數(shù)與水分有效擴(kuò)散系數(shù)的關(guān)系可由下式表示：

確定變量是指明確實驗中的自變量、因變量和無關(guān)變量。在自然狀態(tài)下，造成某一現(xiàn)象的原因可能很多，也就存在多個變量。只有準(zhǔn)確地確定變量，才能具有正確的邏輯推理前提，讓實驗設(shè)計有意義地進(jìn)行下去。一般而言，實驗課題中都包含著自變量和因變量的關(guān)系，因此只有在明確實驗課題的基礎(chǔ)上，學(xué)生才能很容易確定變量。

式中，Deff是干燥過程中的水分有效擴(kuò)散系數(shù)，m2/s；Rg是一個與幾何尺寸有關(guān)的常數(shù)。檳榔干燥過程中的水分有效擴(kuò)散系數(shù)見表2。

表2　檳榔干燥過程中的水分有效擴(kuò)散系數(shù)Table 2 Moisture effective diffusion coefficients in dryingAreca nut

由表2可知，檳榔干燥過程中的估算水分?jǐn)U散系數(shù)Dcal在0.1123×10-7～0.4105×10-7m2/s之間，并且對于烘箱干燥和熱風(fēng)干燥兩種干燥方式，估算水分?jǐn)U散系數(shù)Dcal隨著干燥溫度的提高，都呈現(xiàn)出增大的趨勢。本試驗中檳榔的干燥過程表現(xiàn)出典型的內(nèi)部水分控制的特征，故可以直接運用Fick第二定律計算其水分有效擴(kuò)散系數(shù)。

2.4干燥活化能的計算

干燥活化能表示物料在干燥過程中除去單位的水分所需要的能量。通過干燥活化能的大小可以看出物料干燥的難易程度，并可以估算出干燥能耗［16］。檳榔在干燥過程中的干燥活化能可以根據(jù)阿倫尼烏斯公式進(jìn)行計算：

式中，D0是一個定值，為有效擴(kuò)散系數(shù)的頻率因子，m2/s；Ea為物料的干燥活化能，J/mol；R為氣體摩爾常數(shù)，其值為8.314J/（mol·K）；T為物料的干燥溫度，℃。

將式（8）代入式（9）得，

從圖3中直線的斜率可以求得，烘箱干燥和熱風(fēng)干燥方式處理下的檳榔干燥活化能分別為33.19和32.55 kJ/mol。由此可見，干燥方式對于干燥活化能有一定的影響。

圖3　水分有效擴(kuò)散系數(shù)與干燥溫度的關(guān)系Fig.3 Relation between moisture effective diffusion coefficients and drying temperatures

3　討論

兩種干燥方式（烘箱干燥和熱風(fēng)干燥）檳榔的干燥動力曲線不存在顯著差異，出現(xiàn)這種情況的原因可能是，檳榔是一種多孔性物料，大部分非結(jié)合水主要分布在其表面和孔隙內(nèi)，這部分水分能很快從物料內(nèi)部遷移到物料表面，基本在干燥初期就被除去。隨著干燥的進(jìn)行，檳榔含水率逐漸降低，并且減小得越來越慢，因為這個階段主要去除的是結(jié)合水，而結(jié)合水的遷移需要更久的時間。尤其是干燥的最后時期，基本為降速干燥階段，會花費更多的干燥時間［1］。理論上，干燥溫度較高時的干燥速率應(yīng)該大于干燥溫度較低時的干燥速率，但隨著干燥的進(jìn)行，這種差異會逐漸減弱。本試驗中的檳榔樣品在其原產(chǎn)地已經(jīng)過坑式干燥，由于樣品被放在冰箱中冷藏，才使得測量的初始含水率在36%左右，也就是說，這些水份大部分為檳榔表面的自由水會在干燥初期被迅速除去，從而進(jìn)入降速干燥階段，表現(xiàn)為內(nèi)部水分?jǐn)U散控制，并且本試驗干燥停止時的目標(biāo)含水率也相對較高。另外，加快空氣流動有利于增大傳熱系數(shù)，從而提升傳質(zhì)速率，而試驗中的氣流速度并不高，也就造成了烘箱干燥和熱風(fēng)干燥兩種不同的干燥方式達(dá)到目標(biāo)含水率所用的干燥時間無顯著差異。

從表1可知，烘箱干燥方式的形狀參數(shù)β值在不同干燥溫度下位于0.636～0.679之間，而熱風(fēng)干燥方式下的形狀參數(shù)β在0.619～0.658之間。由此可以看出，對于同一種干燥物料，形狀參數(shù)β的值與干燥方式有關(guān)，而干燥溫度對其影響較小，這跟Corzo［5］的研究結(jié)果一致。本試驗中，不論是烘箱干燥還是熱風(fēng)干燥，檳榔干燥的形狀參數(shù)β值都小于1，也就是說，檳榔的整個干燥過程受內(nèi)部水分?jǐn)U散控制，從圖1-b和圖2-b中也可以看出。由于氣流速度較低，兩種干燥方式的β值沒有顯著性差異。

4　結(jié)論

1）Weibull分布函數(shù)可以很好的描述烘箱干燥和熱風(fēng)干燥方式下的檳榔干燥過程；

2）尺度參數(shù)α與干燥溫度相關(guān)，形狀參數(shù)β與干燥方式有關(guān)。同一干燥方式下，尺度參數(shù)α隨著干燥溫度的升高而減小，而形狀參數(shù)β在不同干燥溫度下變化很??；

3）通過尺度參數(shù)α計算得出的檳榔的估算水分?jǐn)U散系數(shù)Dcal在0.1123×10-7～0.4105×10-7m2/s之間。并且可以根據(jù)阿倫尼烏斯公式計算出烘箱干燥和熱風(fēng)干燥方式下的檳榔干燥活化能分別為33.19和32.55 kJ/mol。

［1］張可喜，符新，王祝年，等.檳榔熱風(fēng)干燥工藝的研究［J］.熱帶農(nóng)業(yè)工程，2006（1）：20-22

［2］林喜娜，王相友.蘋果切片紅外輻射干燥模型建立與評價［J］.農(nóng)業(yè)機(jī)械學(xué)報.2010.41（6）：128-132

［3］張茜，肖紅偉，代建武，等.哈密瓜片氣體射流沖擊干燥特性和干燥模型［J］.農(nóng)業(yè)工程學(xué)報，2011，27（增刊1）：382-388

［4］白峻文，王吉亮，肖紅偉，等.基于Weibull分布函數(shù)的葡萄干燥過程模擬及應(yīng)用［J］.農(nóng)業(yè)工程學(xué)報，2013，29（16）：278-285

［5］Corzo O，Bracho N，Pereira A，et al.Weibull distribution for modeling air drying of coroba slices［J］.LWT-Food Science and Technology，2008，41（10）：2023-2028

［6］李耀明，田進(jìn)，趙湛，等.聯(lián)合收割機(jī)行走底盤變速箱齒輪的疲勞分析［J］.農(nóng)業(yè)工程學(xué)報，2011，27（4）：106-110

［7］盧蓉蓉，錢平，何紅艷，等.超高壓殺滅低酸性食品中耐壓菌的動力學(xué)模型［J］.農(nóng)業(yè)工程學(xué)報，2010，26（9）：350-356

［8］鄭建祥，高翔，周明鋒.基于分段Weibull分布的使用可靠性統(tǒng)計模型［J］.農(nóng)業(yè)機(jī)械學(xué)報，2008，39（6）：50-53

［9］郭劍偉，吳俊珠，王成軍，等.三種方法求解藥物溶出度Weibull分布參數(shù)的探討［J］.大理學(xué)院學(xué)報，2011，10（8）：1-4

［10］Corzo O，Bracho N，Alvarez C.Weibull model for thin-layer drying of mango slices at different maturity stages［J］. Journal of Food Processing and Preservation，2010，34（6）：993-1008

［11］Bantle M，Kolsaker K，Eikevik TM.Modification of the weibull distribution for modeling atmospheric freeze-drying of food［J］.Drying Technology，2011，29（10）：1161-1169

［12］肖紅偉，張世湘，白竣文，等.杏子的氣體射流沖擊干燥特性［J］.農(nóng)業(yè)工程學(xué)報，2010，26（7）：318-323

［13］曾目成，畢金峰，陳芹芹，等.基于Weibull分布函數(shù)對獼猴桃切片中短波紅外干燥過程模擬及應(yīng)用［J］.現(xiàn)代食品科技，2014，30（6）：146-151

［14］Machado MF，Oliveira FAR，Cunha LM.Effect of milk fat and total solids concentration on the kinetics of moisture uptake by ready-to-eat breakfast cereal［J］.International Journal of Food Science and Technology，1999，34（1）：47-57

［15］王蘊博，王寶和.Weibull函數(shù)及其在干燥動力學(xué)研究中的應(yīng)用［J］.干燥技術(shù)與設(shè)備，2011，9（3）：103-109

［16］G??ü? F，Maskan M.Air drying characteristics of solid waste（pomace）of olive oil processing［J］.Journal of Food Engineering，2006，72（4）：378-382

The Simulation for Drying Areca Nut with Weibull Distribution Function

GUO Shou-zhi1，SONG Zhan-hua1，2，LI Fa-de1，2，SONG Hua-lu1，ZHANG Ning1，YAN Yin-fa1，2*
1.College of Mechanical and Electronic Engineering/Shandong Agricultural University，Tai'an 271018，China
2.Shandong Provincial Key Laboratory of Horticultural Machinery and Equipment，Tai'an 271018，China

Application of Weibull distribution was used in this study to describe the drying process of areca nut under different drying methods（dry-chamber drying，hot air drying）and drying temperature（60，80 and 100℃），simulate and analyze the drying kinetics curves.The result indicated that Weibull distribution could provide good prediction for the moisture content.Scale parameter α decreased with the increase of drying temperature.However shape parameter β had a little impact on drying temperatures.Calculated moisture diffusion coefficient Dcalof areca nut samples was calculated，which ranged from 0.1123×10-7to 0.4105×10-7m2·s-1.According to Arrhenius equation，dry activation energy of dry-chamber drying and hot air drying could be calculated to be 33.19 and 32.55 kJ/mol respectively.The result provided a reference for the application of Weibull distribution on areca nut drying.

Drying；Weibull distribution function；areca nut；simulation

TF532+.1

A

1000-2324（2016）05-0710-06

2015-04-24

2015-04-28

郭守志（1988-），男，在讀研究生.研究方向：干燥技術(shù)及干燥機(jī)械裝備研究.E-mail：452136161@qq.com

Author for correspondence.E-mail：sd28@163.com