劉汝寬，楊星星，肖志紅，李昌珠，黃志輝，葉紅齊

（1. 湖南省林業(yè)科學(xué)院 生物能源研究所，湖南 長(zhǎng)沙 410004；2. 中南大學(xué) a.化學(xué)與化工學(xué)院；b.機(jī)電工程學(xué)院，湖南 長(zhǎng)沙 410083）

蓖麻籽冷態(tài)壓榨制油過程中油料散體固相力學(xué)模型及其參數(shù)求解

劉汝寬1,2a，楊星星2b，肖志紅1，李昌珠1，黃志輝2b，葉紅齊2a

（1. 湖南省林業(yè)科學(xué)院 生物能源研究所，湖南 長(zhǎng)沙 410004；2. 中南大學(xué) a.化學(xué)與化工學(xué)院；b.機(jī)電工程學(xué)院，湖南 長(zhǎng)沙 410083）

以蓖麻籽為原料，采用自制的軸向—側(cè)向壓力測(cè)量試驗(yàn)裝置獲取壓榨過程中上端應(yīng)力、側(cè)向應(yīng)力和底部應(yīng)力，建立了油料固相應(yīng)力模型，并求解得到了模型參數(shù)，相對(duì)誤差為9.80%，可用于描述蓖麻籽壓榨過程中沿榨筒軸向固相顆粒應(yīng)力分布。根據(jù)油料固相應(yīng)力模型可知，在壓榨油料確定的情況下，油料散體軸向應(yīng)力主要受榨筒內(nèi)圈半徑r、油料散體高度H影響，即軸向應(yīng)力隨物料高度的增加逐漸減小，隨榨筒半徑增大而增大。同時(shí)，結(jié)合流體力學(xué)可知，半徑越大，油料流動(dòng)阻力越大，因此榨筒半徑不宜過大。

蓖麻籽；冷態(tài)壓榨；散體模型

蓖麻是一種重要的高含油特種植物油料，其油脂中90%以上為富含羥基蓖麻油[1]，有利于改性制備可再生化工產(chǎn)品[2]。蓖麻籽多采用壓榨或浸提等方式制備蓖麻油[3]，以壓榨制油較為常見。該法又有螺旋和液壓制油等操作方式[4]，前者可以獲得較低的餅粕殘油率（約8%），制油效率較高，但是油料會(huì)經(jīng)過一個(gè)高溫（＞200 ℃）壓縮過程[5]，不利于后續(xù)高蛋白餅粕[6]（蛋白含量＞40%）的再利用。因此，若采用直筒式液壓制油，一定程度上可以保證油料加工過程中的低溫操作，再采用溶劑浸提[7]的方式提取殘油（約15%），整個(gè)過程有利于餅粕的再利用，進(jìn)而能夠?qū)崿F(xiàn)蓖麻加工綜合效益的提升。

目前，針對(duì)油料直筒式壓榨的研究相對(duì)較少，主要是其工藝過程不連續(xù)，多應(yīng)用于特種油料的加工。筆者[8]已設(shè)計(jì)出連續(xù)式低溫直筒壓榨系統(tǒng)，實(shí)現(xiàn)了油料連續(xù)加工。為實(shí)現(xiàn)壓榨后餅粕的殘油率降低，應(yīng)系統(tǒng)開展直筒式液壓壓榨制油研究。針對(duì)油料壓榨制油過程的研究，早期關(guān)注的是出油率（或殘油率）與壓力及壓榨時(shí)間的關(guān)系，基本是實(shí)驗(yàn)性總結(jié)得到的關(guān)系式[9-12]，未見針對(duì)蓖麻籽壓榨過程中散體固相力學(xué)的分析。

本研究重點(diǎn)圍繞蓖麻籽散體顆粒在直筒式壓榨制油中的受力狀態(tài)，構(gòu)建其固相力學(xué)模型，并對(duì)其模型參數(shù)求解，為多工位連續(xù)式低溫液壓制油裝置的研制提供參考。

1 材料與方法

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

1.1.1 試驗(yàn)材料

蓖麻籽（湘蓖1號(hào)，自然干燥，含水量6.0%）：湖南省林業(yè)科技示范園。

1.1.2 試驗(yàn)設(shè)備

游標(biāo)卡尺（0～150 mm），上海力衡儀器儀表有限公司；液壓加載保壓試驗(yàn)機(jī)（WEW），上海和晟儀器科技有限公司；軸向—側(cè)向壓力測(cè)量試驗(yàn)裝置（套筒內(nèi)徑39 mm，見圖1），自制。

圖1 軸向—側(cè)向壓力搜集試驗(yàn)裝置Fig.1 Experimental device for axial and lateral pressure

1.2 試驗(yàn)方法

1.2.1 直筒壓榨過程的基本假設(shè)

壓榨時(shí)油料散體主要沿著軸向向下壓縮滑動(dòng)，假設(shè)壓榨過程中滿足如下條件：

① 同一高度上的油料散體固相顆粒在壓榨過程中始終保持在同一水平面；

② 忽略榨筒側(cè)面細(xì)孔對(duì)油料散體中固體顆粒的影響；

③ 忽略油料散體重力、慣性力及油料散體固相內(nèi)部徑向摩擦的影響；

④ 在微小單元厚度范圍內(nèi)，徑向應(yīng)力基本不變，并引入側(cè)壓系數(shù)[13-14]。

1.2.2 油料受力參數(shù)搜集

為同時(shí)獲取油料底部軸向應(yīng)力、油料上端軸向應(yīng)力及上端側(cè)向應(yīng)力，利用自制試驗(yàn)設(shè)備——軸向—側(cè)向壓力測(cè)量試驗(yàn)裝置進(jìn)行壓榨，讀取油料壓榨時(shí)上端壓力、側(cè)向壓力以及底部壓力。側(cè)向壓力傳感器（NS-F，量程0～500 kPa）用于測(cè)量側(cè)壓，底部壓力傳感器（HY603稱重傳感器，量程0～2.5 kN）用于測(cè)量油料底部壓力，其中油料上端壓力通過測(cè)量加裝液壓缸壓力值確定。

1.2.3 試驗(yàn)分組

為使試驗(yàn)達(dá)到不同壓力值（即不同應(yīng)變），改變每次加料高度H0，且油料壓榨完成后高度均為H1，即可達(dá)到不同應(yīng)變，從而產(chǎn)生不同應(yīng)力值。為使應(yīng)變分別能達(dá)到0.20、0.25、0.30、0.35、0.40、0.45和0.50，每次試驗(yàn)加料高度分別為63、67、71、77、83、91和100 mm，共7組試驗(yàn)。

1.2.4 直筒式壓榨操作

在套筒中加入一定量的蓖麻籽，利用液壓加載系統(tǒng)進(jìn)行加載[15]，加載速度選用v=2.2 mm/s，根據(jù)加料高度H0和壓完后油料高度H1，控制活塞桿加載位移為H0-H1，然后停止加載完成壓榨操作。直筒式壓榨模型見圖2。

圖2 直筒式壓榨模型Fig.2 Model for straight tube press

2 結(jié)果與分析

2.1 油料散體固相受力分析

2.1.1 油料散粒體在直筒內(nèi)的受力分析

蓖麻籽直筒壓榨時(shí)，榨筒與油料始終相互作用，許多大小不同、形狀各異的固體顆粒所構(gòu)成的機(jī)械混合物稱為散粒體或散體[16]，這種油料散體主要是由固相顆粒和液相油脂組成的多孔介質(zhì)[17]。蓖麻籽在活塞的軸向作用下進(jìn)行壓榨，初始階段油料處于散粒體狀態(tài)，油脂并未開始流出榨筒，可忽略孔隙流體壓力，油料散體固相顆粒承受活塞全部軸向正壓力，此時(shí)油料的壓榨受力情況與粉體成形過程中粉體顆粒相似[13]。油料不斷被壓縮，油料中固相顆粒開始向下做壓縮運(yùn)動(dòng)，由于向下運(yùn)動(dòng)的油料與榨筒壁產(chǎn)生摩擦力，消耗了一部分固相顆粒承受的軸向應(yīng)力，使其不能毫無損耗地向下傳遞，進(jìn)而影響油料壓榨效果。油料在被壓縮的過程中，除軸向壓縮外，徑向上油料散體力向側(cè)面膨脹和變形，但這種趨勢(shì)會(huì)受到榨筒筒壁對(duì)其的側(cè)向限制，由此產(chǎn)生一個(gè)側(cè)向應(yīng)力。

2.1.2 蓖麻籽直筒壓榨制油時(shí)有限元分析

油料散體受力情況如圖3所示，其中H為油料高度（即壓榨過程中榨筒內(nèi)植物油料的實(shí)時(shí)高度）；σg(H)為油料散體底部受到油盤的應(yīng)力；r為圓柱形榨筒內(nèi)圈半徑，稱為壓榨半徑；Ff為油料散體傳遞給榨筒的摩擦力，也等于油盤對(duì)榨筒的支撐力；σg為整個(gè)油料散體固相顆粒所受軸向應(yīng)力。

圖3 油料散體固相受力情況(a, 左)及其有限元分析(b, 右)Fig.3 Analysis of solid state(left) and its fi nite element (right) for oilseeds

油料散體固相顆粒壓榨過程中所受外力關(guān)于榨筒軸線軸向?qū)ΨQ，可將其轉(zhuǎn)化為平面問題。利用有限元思想，在z坐標(biāo)軸上h處取一微小單元的油料散體固相顆粒層，厚度為dh，微元受力見圖 3（b）。其中σg(h)、σg(h-dh)分別為所取油料散體固相微元層下表面和上表面所受應(yīng)力；dFf為微元所受摩擦力；σr(h)為油料散體h處對(duì)榨筒內(nèi)壁的側(cè)向應(yīng)力。

2.2 油料散體固相力學(xué)模型

2.2.1 摩擦力的計(jì)算

根據(jù)油料散體固相受力分析，油料散體與榨筒內(nèi)壁之間的摩擦力是由于油料散體有側(cè)向應(yīng)力作用，當(dāng)兩者發(fā)生相對(duì)運(yùn)動(dòng)時(shí)即產(chǎn)生摩擦力，該摩擦力大小取決于油料散體對(duì)榨筒內(nèi)壁的側(cè)向應(yīng)力、油料散體與榨筒內(nèi)壁的摩擦系數(shù)以及摩擦面面積，由于微元薄片厚度極小，可視側(cè)向應(yīng)力在軸向基本不變，對(duì)圖3中所取的油料散體微元有：

dFf= 2πrfσr(h)dh。 （1）式中：f為油料散體固相顆粒與榨筒內(nèi)壁之間的摩擦系數(shù)；r為榨筒內(nèi)圈半徑，即壓榨半徑，mm；σr(h)為油料散體固相對(duì)榨筒內(nèi)壁的側(cè)向應(yīng)力，MPa；dFf為微元所受摩擦力，N；dh為微元厚度，mm。

2.2.2 側(cè)壓系數(shù)的計(jì)算

由側(cè)壓系數(shù)定義[13]及假設(shè)條件可得：

式中 ：ξ為側(cè)壓系數(shù)。

根據(jù)文獻(xiàn)[14]中多種食用油料壓榨過程中不同種油料側(cè)壓系數(shù)經(jīng)驗(yàn)公式，可統(tǒng)一表達(dá)為：

式中：A、B為常數(shù)，其值取決于油料物理特性，不同油料取值不同。

2.2.3 軸向力的計(jì)算

z坐標(biāo)軸上任意z處固相所受軸向力等于整個(gè)油料散體所受活塞壓力與坐標(biāo)0到z之間摩擦力之差，即：

式中：Fz為z坐標(biāo)軸上z處固相所受軸向力；F0為z坐標(biāo)軸上0處軸向力，即整個(gè)油料散體所受壓力；Ff|z

0為坐標(biāo)0到z之間油料散體固相所受摩擦力。

2.2.4 固相力學(xué)模型的構(gòu)建

聯(lián)立式(1)、(2)和(3)積分，代入式(4)，即：

將式（5）積分方程求導(dǎo)轉(zhuǎn)化為微分方程，并利用伯努利方程的求解方法求解，并代入初始條件σg(0)=σg，得到其解為：

式（6）即為壓榨過程中油料散體固相顆粒所受軸向應(yīng)力沿軸向的變化力學(xué)模型。式（6）中摩擦系數(shù)f、側(cè)壓系數(shù)中油料特性常量A和B均為油料物理特性參數(shù)，非壓榨過程所能控制，因此油料散體力學(xué)模型主要受榨筒內(nèi)圈半徑r（稱為壓榨半徑）、油料散體高度H影響。

2.3 固相力學(xué)模型參數(shù)的求解及其驗(yàn)證

2.3.1 側(cè)壓模型參數(shù)的計(jì)算

式（6）中與油料性質(zhì)有關(guān)的參數(shù)A、B的值，可通過油料上端軸向應(yīng)力與油料上端處側(cè)向應(yīng)力確定。由于壓榨出油后，活塞壓力同時(shí)由孔隙中流體和油料固相共同承擔(dān)，超出固相力學(xué)模型的范圍，取應(yīng)變值在0～0.6之間。將底部壓力及加載壓力換算成相應(yīng)應(yīng)力值，結(jié)果如表1所示。

由式（2）和（3）可知，油料上端應(yīng)力與側(cè)向應(yīng)力滿足以下關(guān)系：

代入表1中油料上端應(yīng)力值及側(cè)向應(yīng)力值，利用Matlab按式（7）進(jìn)行非線性擬合，得到蓖麻籽直筒壓榨對(duì)應(yīng)的側(cè)壓模型參數(shù)：A=0.012 6，B=0.003 2。

2.3.2 模型的驗(yàn)證

榨筒半徑及壓榨高度不僅影響油料壓榨過程中固相顆粒受力情況，而且決定了壓榨機(jī)總壓榨量，有必要對(duì)油料散體固相模型進(jìn)行驗(yàn)證分析。將A、B、壓榨后油料高度H1及表1中油料上端應(yīng)力值代入式（2）～（6），得到油料底部應(yīng)力模擬值，并與試驗(yàn)值對(duì)比（如圖4所示）。

圖4 油料底部應(yīng)力試驗(yàn)值與理論值對(duì)比Fig.4 Comparison of experimental and theoretical values

圖4中油料底部應(yīng)力理論值與試驗(yàn)值平均相對(duì)誤差為9.80%，可見參數(shù)A、B和式（6）所表示的固相力學(xué)模型適用于描述蓖麻籽壓榨過程中沿榨筒軸向固相顆粒應(yīng)力分布情況。

2.4 模型的應(yīng)用分析

實(shí)際生產(chǎn)中所關(guān)注的是油料散體最下面一層油料是否達(dá)到壓榨出油中的應(yīng)力，可將油料散體力學(xué)模型式（6）中z換為油料高度H，即可得到底部油料散體所受應(yīng)力σg(H)。設(shè)整個(gè)油料散體軸向壓榨應(yīng)力σg=50 MPa，取蓖麻的物性參數(shù)A=0.012 6，B=0.003 2，農(nóng)作物與金屬鋼板的摩擦角范圍為10°～33°[16]，油料散體與榨筒內(nèi)壁間的摩擦系數(shù)f=0.25，可得油料散體力學(xué)模型與油料高度H、榨筒半徑r的關(guān)系（如圖5和圖6所示）。

結(jié)合圖5和圖6可知，當(dāng)榨筒半徑為定值時(shí)，油料散體底部軸向應(yīng)力σg(H)隨著油料高度的增加而減小；當(dāng)油料高度為定值時(shí)，軸向應(yīng)力σg(H)隨著榨筒半徑增大而增大；并且，當(dāng)榨筒半徑增加到一定值后油料散體底部軸向應(yīng)力σg(H)有著不再變化的趨勢(shì)，結(jié)合流體力學(xué)可知，半徑越大，油料流動(dòng)阻力越大，因此榨筒半徑不宜過大。

3 結(jié) 論

以蓖麻籽為原料，采用自制的軸向—側(cè)向壓力測(cè)量試驗(yàn)裝置獲取蓖麻籽壓榨過程中上端應(yīng)力、側(cè)向應(yīng)力和底部應(yīng)力，建立了壓榨過程中油料固相應(yīng)力模型，求解得到了模型參數(shù)，相對(duì)誤差為9.80%，可用于描述蓖麻籽壓榨過程中沿榨筒軸向固相顆粒應(yīng)力分布。

圖5 底部軸向應(yīng)力σg(H)與油料高度H的關(guān)系Fig.5 Relationship between bottom axial stress and height of oilseeds

圖6 底部軸向應(yīng)力σg(H)與榨筒半徑r的關(guān)系Fig.6 Relationship between the bottom axial stress and the cylinder radius

根據(jù)獲得的油料固相應(yīng)力模型可知，在壓榨油料確定的情況下，油料散體軸向應(yīng)力主要受榨筒內(nèi)圈半徑r、油料散體高度H影響，即軸向應(yīng)力隨物料高度的增加逐漸減小，隨榨筒半徑增大而增大。同時(shí)，結(jié)合流體力學(xué)可知，半徑越大，油料流動(dòng)阻力越大，因此榨筒半徑不宜過大。

[1]張 研, 喬金友. 淺析中國(guó)蓖麻產(chǎn)業(yè)化發(fā)展前景[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2009, 25(16): 316-319.

[2]劉汝寬,肖志紅,李昌珠,等. 正交法優(yōu)化植物油脂調(diào)配潤(rùn)滑防銹油[J]. 糧油加工, 2009, (2): 68-69.

[3]劉汝寬, 施亮林, 肖志紅, 等. 壓榨時(shí)間和壓力對(duì)蓖麻籽壓榨性能的影響及其數(shù)學(xué)擬合研究[J]. 中國(guó)糧油學(xué)報(bào), 2014,29(10): 51-55.

[4]劉汝寬,肖志紅,李昌珠,等. 光皮樹果實(shí)內(nèi)含物含量的測(cè)定與綜合利用[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2013, 33(12): 19-21.

[5]解士聰, 黃志輝, 李昌珠, 等. 1種螺旋榨油機(jī)榨膛內(nèi)表面溫度預(yù)估方法[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2014, 35(6): 104-107.

[6]Liu Rukuan, Xu Fanglei, Xiao Zhihong,et al.Model of Plug Flow and Equivalent Friction Coefficient in Process of Castor Beans by Single Screw Press[J]. Journal of Bioprocess Engineering and Biore fi nery, 2014, 3(3): 227-231.

[7]Zeng Fantao, Liu Rukuan, Xiao Zhihong,et al.Castor Oil Extraction Using Subcritical n-Butane under Low Temperature and Refining Process of Oil[J]. Journal of Bioprocess Engineering and Biore fi nery, 2014, 3(3): 196-201.

[8]劉汝寬, 黃志輝, 李昌珠, 等. 一種植物油料壓榨試驗(yàn)方法及裝置: 中國(guó), 201310270033.3[P]. 2013-10-30.

[9]李昌珠, 肖志紅, 劉汝寬, 等. 光皮樹果實(shí)專用型螺旋冷榨機(jī)的壓榨性能研究[J]. 湖南林業(yè)科技, 2013, 40(2): 6-8.

[10]劉汝寬,肖志紅,張愛華,等. 基于松弛試驗(yàn)的光皮樹果實(shí)冷態(tài)壓榨流變模型[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(16): 278-283.

[11]趙 揚(yáng). 濾餅微觀結(jié)構(gòu)與壓榨過濾理論的研究[D]. 杭州:浙江大學(xué), 2006.

[12]Bargale P C. Mechanical oil expression from selected oilseeds under uniaxial compression[C]. PhD Thesis, Department of Agricultural & Bioresource Engineering, University of Saskatchewan, Saskatoon, Canada, 1997.

[13]吳成義, 張麗英. 粉體成形力學(xué)原理[M]. 北京: 冶金工業(yè)出版社, 2003.

[14]黃志輝, 許方雷, 李昌珠,等. 蓖麻籽不同狀態(tài)不同壓榨方式下壓縮比與軸向應(yīng)力關(guān)系的研究[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào),2015, 36(1): 101-105.

[15]楊星星. 蓖麻籽多工位直筒式壓榨工藝及控制功能規(guī)劃研究[D]. 長(zhǎng)沙:中南大學(xué), 2014.

[16]黃松元. 散體力學(xué)[M]. 北京：機(jī)械工業(yè)出版社, 1993.

[17]鄭 曉, 宛 農(nóng), 林國(guó)祥, 等. 基于孔隙度的冷榨菜籽仁餅微觀結(jié)構(gòu)研究[J]. 中國(guó)油脂, 2005, 29(12): 14-17.

Solid model and its parameters for castor beans in cold press

LIU Ru-kuan1,2a, YANG Xing-xing2b, XIAO Zhi-hong1, LI Chang-zhu1, HUANG Zhi-hui2b, YE Hong-qi2a
(1. Institute of Bioenergy, Hunan Academy of Forestry, Changsha 410004, Hunan, China; 2a. College of Chemistry and Chemical Engineering; b. College of Mechanical and Electrical Engineering, Central South University, Changsha 410083, Hunan, China)

Castor is an important oil plant, of which the seeds are rich in castor oil usually obtained by pressing instruments. In this paper, pressures of castor bean in top, lateral and bottom were obtained using self-made axial and lateral pressure measurement test device, and solid model for oilseeds was founded. The model parameters were obtained with the relative error was 9.80%, which can be used to describe the distribution of axial stress along the cylinder in squeezing process of castor beans. Based on the model, the axial stress is mainly affected by inner radius of squeezing barrel and height of oilseeds, which means the axial stress decreases with the increase of the height of the oilseeds, and increases with the increase of the cylinder radius. According to fl uid mechanics, the larger the radius is, the greater the oil fl ow resistance, leading to a result that the cylinder radius could not too large.

Castor bean; cold press; dispersion model

S789.7；TQ641

A

1673-923X(2016)05-0133-05

10.14067/j.cnki.1673-923x.2016.05.024

2014-12-12

林業(yè)公益性行業(yè)專項(xiàng)（201304608）；國(guó)家自然科學(xué)基金項(xiàng)目（31470594）

劉汝寬，副研究員，博士 通訊作者：肖志紅，研究員；E-mail：xzhh1015@163.com

劉汝寬， 楊星星， 肖志紅，等. 蓖麻籽冷態(tài)壓榨制油過程中油料散體固相力學(xué)模型及其參數(shù)求解[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2016, 36(5): 133-137.

[本文編校：謝榮秀]