杜小翠 嚴曉婕 程緒鐸 高夢瑤 馮家暢

(南京財經(jīng)大學(xué)食品科學(xué)與工程學(xué)院；江蘇省現(xiàn)代糧食流通與安全協(xié)同創(chuàng)新中心，南京 210046)

帶錐斗筒倉中稻谷的密度與應(yīng)力分布模型

杜小翠 嚴曉婕 程緒鐸 高夢瑤 馮家暢

(南京財經(jīng)大學(xué)食品科學(xué)與工程學(xué)院；江蘇省現(xiàn)代糧食流通與安全協(xié)同創(chuàng)新中心，南京 210046)

采用LHT-1糧食回彈模量儀測定稻谷堆的壓縮密度，建立帶錐斗筒倉中稻谷堆的密度、應(yīng)力與糧層深度關(guān)系的微分方程組，用數(shù)值方法計算帶錐斗筒倉中稻谷密度、應(yīng)力與糧層深度關(guān)系，由積分法計算出筒倉中稻谷的儲藏總質(zhì)量。試驗結(jié)果表明，淮稻5號(含水率為10.38%～18.30% w.b.)的密度隨豎直應(yīng)力(0.495～245.892 kPa)增大而增大(582.772～696.593 kg/m3)。模型計算結(jié)果表明，在帶錐斗筒倉的筒體部分，稻谷堆密度隨著糧層深度的增加而增大；到錐斗部分，稻谷堆密度隨著糧層深度的增加而逐漸減小。在帶錐斗筒倉的筒體部分，稻谷堆的豎直應(yīng)力隨著糧層深度的增加而增大；在錐斗部分，稻谷堆的豎直應(yīng)力則隨著糧層深度的增加而減小。在帶錐斗筒倉中的筒體部分，稻谷堆的側(cè)向應(yīng)力隨著糧層深度的增大而增大；在筒體與錐斗結(jié)合處，稻谷堆的側(cè)向應(yīng)力突然增加；到了錐斗部分，稻谷堆的側(cè)向應(yīng)力隨著糧層深度的增大先稍增大再逐漸減小。

帶錐斗筒倉 稻谷 密度 豎直應(yīng)力 側(cè)向應(yīng)力

稻谷是我國單產(chǎn)最高、總產(chǎn)最多的糧食作物，在國民經(jīng)濟和國家糧食安全中具有重要的地位和作用[1]。筒倉是儲糧倉型中機械化程度最高的一種倉型，被廣泛應(yīng)用于稻谷的儲藏[2-3]。稻谷儲藏在筒倉中，會受到自重、內(nèi)摩擦力以及倉壁的支持力，稻谷堆的內(nèi)部由于這些力的作用產(chǎn)生應(yīng)力，從而產(chǎn)生彈性和塑性形變，稻谷堆的體積縮小，密度則增大，隨著糧層深度的增加，稻谷的密度增大[4-5]。準確地計算出筒倉中密度的分布就能準確的計算筒倉中稻谷的總質(zhì)量。目前我國糧食儲藏質(zhì)量檢查方法是體積密度法[6]。體積密度法是將稻谷堆的體積與平均密度相乘來計算稻谷質(zhì)量，稻谷堆的平均密度為常規(guī)測定的密度乘以修正系數(shù)，修正系數(shù)是憑經(jīng)驗給出的，這個方法計算的糧倉中稻谷質(zhì)量誤差大。儲藏在筒倉中的稻谷壓應(yīng)力隨著糧層深度的增加而增大。筒倉中稻谷堆的壓應(yīng)力隨糧層深度的分布值是糧倉結(jié)構(gòu)設(shè)計所需要的重要參數(shù)。

Mahmoud[7]使用一個非線性雙曲線型本構(gòu)方程，建立了有限元模型，估算了由彈性波紋板構(gòu)成的圓筒倉倉壁靜態(tài)壓力。Zhang等[8]利用Lade模型計算了糧倉的靜態(tài)倉壁壓力，并研究了溫度對倉壁側(cè)壓力的影響。Ruckenbrod等[9]和Ragneau等[10]使用不同的物料本構(gòu)方程通過有限元方法估算了圓筒倉中物料的靜態(tài)壓力和卸料動壓力。Goode等[11]通過有限元方法估算了薄壁矩形截面筒倉的內(nèi)小麥的靜壓力。上述對筒倉中物料壓應(yīng)力的研究都是將物料密度視為常量，沒有考慮密度的變化對應(yīng)力分布的影響。Loewer等[12]觀測了小麥的堆密度隨壓力和水分的變化值。Thompson 等[13]試驗測定了小麥的堆密度隨壓力和水分的變化值并給出了預(yù)測密度的經(jīng)驗方程。Moya等[14]使用固結(jié)儀測定幾種谷物的孔隙率與密度與豎直壓力的關(guān)系。這些研究只是在小容器中測定密度與壓力的關(guān)系。陳德民等[15]利用壓力傳感器網(wǎng)絡(luò)進行糧食儲藏質(zhì)量在線監(jiān)測，給出了壓力傳感器外圍電路的硬件設(shè)計，并根據(jù)我國糧倉特點對傳感器進行了布置，利用VC++6.0開發(fā)了系統(tǒng)管理軟件。在筒壁和底部安裝力傳感器，由倉底與倉壁的壓力分布推出糧倉儲糧總質(zhì)量，這個方法本質(zhì)上是由倉壁與倉底壓力分布估算平均密度，再算出總質(zhì)量。這個方法測算的誤差大、成本大，難以推廣。到目前為止，筒倉中糧堆的密度分布和裝糧質(zhì)量的理論研究成果鮮見報導(dǎo)，也鮮有有效的測量方法。

本試驗采用LHT-1糧食回彈模量儀來測定稻谷堆在不同壓應(yīng)力下的壓縮密度，用試驗數(shù)據(jù)擬合壓縮密度與壓應(yīng)力的關(guān)系模型，建立筒倉中稻谷堆的密度、壓應(yīng)力與糧層深度關(guān)系的微分方程組，用數(shù)值方法計算得到筒倉中稻谷密度隨糧層深度變化關(guān)系，從而由積分法計算出筒倉中稻谷的儲藏總質(zhì)量。同時，用數(shù)值方法計算得到筒倉中隨糧層深度變化的稻谷堆的豎直應(yīng)力、側(cè)向應(yīng)力的分布值。

1 稻谷堆密度與最大主應(yīng)力關(guān)系的試驗測定

稻谷儲藏在筒倉的圓筒中，處于主動應(yīng)力狀態(tài)，最大主應(yīng)力是豎直方向；稻谷儲藏在筒倉的錐斗中，處于被動應(yīng)力狀態(tài)，最大主應(yīng)力垂直錐斗壁。密度與最大主應(yīng)力的關(guān)系(ρ=F(pmix))通過試驗來確定。

1.1 材料與方法

1.1.1 試驗材料

淮稻5號稻谷，原始水分為15.23%(濕基)，2014年產(chǎn)于江蘇淮安。2014年11至12月進行壓縮試驗。

試驗時，將其水分調(diào)節(jié)為10.38%、12.38%、13.77%、15.23%、16.99%、18.30% w.b.。稻谷籽粒的最大直徑、中直徑、最小直徑的平均值分別為7.5、3.6、2.4 mm，對應(yīng)的標準差分別為0.2、0.1、0.1 mm。(隨機取樣100粒稻谷粒，用游標卡尺測得其3個粒徑，精確到0.1 mm，最后求其平均值)

1.1.2 試驗儀器

LHT-1型糧食回彈模量測定儀：南京土壤儀器廠有限公司；HG202-2(2A/2AD)電熱干燥箱：南京盈鑫實驗儀器有限公司；AL204型分析天平、JSFD-粉碎機：上海嘉定糧油儀器有限公司。

1.1.3 試驗原理

將稻谷裝入LHT-1糧食回彈模量儀的圓筒中，頂部加載荷，圓筒中的稻谷處于主動應(yīng)力狀態(tài)，稻谷在圓筒中產(chǎn)生限側(cè)變形壓縮[16]，這種應(yīng)力狀態(tài)與變形狀態(tài)與筒倉中的稻谷應(yīng)力狀態(tài)與變形狀態(tài)相似。因此，可通過對回彈模量儀的圓筒中的糧堆施加不同的頂部壓力來模擬筒倉中不同深度的稻谷的應(yīng)力狀態(tài)與應(yīng)變狀態(tài)。筒壁對稻谷產(chǎn)生向上的摩擦力，裝樣筒內(nèi)各糧層所受到的豎直應(yīng)力與側(cè)向應(yīng)力是不同的，糧層深度增加，豎直應(yīng)力與側(cè)向應(yīng)力減小。所以，本節(jié)用圓筒中稻谷的平均豎直應(yīng)力與側(cè)向應(yīng)力表示筒倉中某一糧層的豎直應(yīng)力與側(cè)向應(yīng)力。

采用微元法對裝樣筒內(nèi)糧堆樣品進行受力分析，假設(shè)裝樣筒內(nèi)半徑為Rc，樣品高度為Hc，取距離糧面深度為y，高度為dy的微元體進行受力分析[5]，微元體在豎直方向的受力分析見圖1：

圖1 回彈模量儀裝樣筒中糧堆微元層受力圖

稻谷在裝樣筒內(nèi)受壓時處于主動應(yīng)力狀態(tài)，根據(jù)Rankine理論[17]和Janssen理論[18]，側(cè)向應(yīng)力與豎直應(yīng)力的比為：

(1)

式中：ph1為圓筒中糧堆所受側(cè)向應(yīng)力/kPa；pv1為圓筒中糧堆所受豎直應(yīng)力/kPa；φ為稻谷內(nèi)摩擦角/°。

糧堆微元體在豎直方向上受力平衡，可得平衡方程為：

(2)

式中：ρ為稻谷密度/kg/m3；μc為稻谷與裝樣筒筒壁的摩擦系數(shù)；Ac為裝樣筒的橫截面積/m2。

由式(2)，結(jié)合式(1)推導(dǎo)出下面的微分方程：

(3)

對式(3)進行積分，得：

(4)

(5)

式中：p0為糧堆頂部壓應(yīng)力/kPa；Hc為壓縮后圓筒中糧堆高度/m。

在裝入稻谷的圓筒頂部逐級加載，測定每次加載后稻谷樣品的高度，計算出對應(yīng)的稻谷密度，并根據(jù)式(5)計算出對應(yīng)的豎直應(yīng)力，即可得到稻谷樣品密度與其所受豎直應(yīng)力(最大主應(yīng)力)的關(guān)系：

ρ=F(pmax)

(6)

1.1.4 試驗方法

將回彈模量儀放在實驗室中，室內(nèi)溫度設(shè)置為20 ℃。將樣品勻速倒入裝樣筒中，并將其表面鋪平，放上傳壓板，保證傳壓板上表面與裝樣筒上端齊平。保持橫梁杠桿垂直，轉(zhuǎn)動平衡錘調(diào)整杠桿至水平以上，用M16螺母固定平衡錘。旋轉(zhuǎn)傳壓螺釘與傳壓板接觸，調(diào)整0～30 mm位移傳感器的觸頭位置，調(diào)零百分表。根據(jù)試驗要求，對試樣進行加載。隨著試樣的下沉，杠桿向下傾斜，為防止杠桿傾斜影響加載精度，調(diào)節(jié)調(diào)平手輪，使杠桿處于水平位置。試驗進行3 d后記錄樣品高度，計算樣品壓縮后的體積，結(jié)束試驗，倒出裝樣筒內(nèi)的稻谷樣品，稱重并記錄稻谷質(zhì)量。稻谷堆的密度即為稻谷的質(zhì)量除以壓縮后的體積[5]。

1.2 結(jié)果與分析

1.2.1 平均豎直應(yīng)力

淮稻5號稻谷的內(nèi)摩擦角和摩擦系數(shù)隨含水率的變化不顯著，不同壓應(yīng)力條件下內(nèi)摩擦角和摩擦系數(shù)不同。由直剪儀測得的試驗結(jié)果見表1。

使用回彈模量儀通過逐級加載的方式對稻谷堆進行壓縮試驗，加載的頂部豎直應(yīng)力分別為50、100、150、200、250、300 kPa。根據(jù)表1 中的參數(shù)與試驗測定的樣品壓縮的頂部壓力與樣品高度，根據(jù)式(5)計算得到稻谷堆所受的平均豎直應(yīng)力，結(jié)果見表2。

表1 不同壓應(yīng)力下稻谷的內(nèi)摩擦角和摩擦系數(shù)

1.2.2 稻谷堆的壓縮密度

不同含水率、不同頂部壓應(yīng)力下的稻谷堆壓縮后測定的密度見表3。

1.2.3 稻谷堆壓縮密度與平均豎直應(yīng)力的關(guān)系

根據(jù)表2和表3作圖，可以得到稻谷堆壓縮密度與平均豎直應(yīng)力的關(guān)系曲線，如圖2所示。

圖2 不同含水率下的稻谷堆密度與平均豎直應(yīng)力的關(guān)系曲線

由圖2可以得到，在同一含水率下，隨著平均豎直應(yīng)力的增大，稻谷堆的壓縮密度增大。根據(jù)圖2可擬合出稻谷堆密度與平均豎直應(yīng)力的關(guān)系方程為ρ=-apv2+bpv+c(ρ=F(pmax))，其中方程系數(shù)a、b、c隨含水率的變化而變化，結(jié)果見表4。

表2 不同含水率、不同頂部壓力下的稻谷堆所受的平均豎直應(yīng)力/kPa

表3 不同含水率、不同頂部壓力下的稻谷堆的壓縮密度/kg/m3

表4 平均豎直應(yīng)力與稻谷堆密度的擬合方程系數(shù)

2 帶錐斗筒倉內(nèi)稻谷密度、壓應(yīng)力與糧層深度關(guān)系的微分方程組

設(shè)帶錐斗筒倉的直徑為D，圓筒高為H，錐斗高為h，錐斗壁與水平面夾角為α，稻谷的內(nèi)摩擦角為φ，稻谷儲藏在倉中，倉墻與稻谷間的摩擦系數(shù)為μ，設(shè)水平壓應(yīng)力與豎直應(yīng)力大小不隨半徑的大小和方向而變化。帶錐斗的筒倉和稻谷表面下深度y處稻谷的微元層受力示意圖如圖3所示。

對于筒倉的筒體部分，其微元層的受力如圖3b所示。由于微元層在豎直方向上合力為零，則得到：

(pv+dpv)A+μphCdy-pvA-ρgAdy=0

(7)

式中：pv為稻谷豎直應(yīng)力/kPa；ph為側(cè)應(yīng)壓力/kPa；ρ為稻谷的密度/kg/m3；g為重力加速度/m/s-2；A為圓筒的橫截面面積/m2；C為圓筒的周長/m；且

(8)

式中：R為水力半徑/m；R=D/4。

圖3 筒倉的筒體中微元體受力圖

稻谷靜態(tài)儲藏在筒倉中，稻谷處于主動應(yīng)力狀態(tài)，根據(jù)Janssen理論[18]，稻谷最小主應(yīng)力(即筒壁側(cè)壓力ph)和最大主應(yīng)力(即豎直應(yīng)力pv)之比k為：

(9)

結(jié)合式(6)～式(8)得到筒體部分最大主應(yīng)力pv的微分方程：

(10)

而筒體部分最小主應(yīng)力ph為：

ph=kpv

(11)

對于筒倉的錐斗部分，其截面圖和微元層受力如圖4所示。在同一水平面上的錐斗部分微元層，其糧堆的最大主應(yīng)力的大小與方向都是變化的。在錐斗壁處稻谷處于被動應(yīng)力狀態(tài)，此處最大主應(yīng)力方向垂直于錐斗壁。筒體底層中心的最大主應(yīng)力與錐斗頂面中心最大主應(yīng)力是相等的。錐斗內(nèi)的任一水平微元層，中心豎直應(yīng)力(pvc)是該處的最大主應(yīng)力，邊緣豎直應(yīng)力并不是該處的最大主應(yīng)力，也不是該處的最小主應(yīng)力，在錐斗壁處(微元層邊緣處)取一微元體如圖4b、圖4c 所示，微元體受到的最大和最小壓應(yīng)力分別為phmax和phmin，phmax為錐斗壁側(cè)壓力，且微元體的最小壓應(yīng)力和最大壓應(yīng)力之比為k。

圖4 筒倉錐斗部分微元體受力圖

根據(jù)Rankine壓力理論[22]得到水平稻谷薄層邊緣處豎直應(yīng)力pve：

(12)

令

(13)

同一水平微元層的豎直應(yīng)力取平均值phv：

phv=(pve+pvc)/2=(λ+1)phmax/2

(14)

作用于微元層上豎直方向的合力等于0：

phvAa+ρgV-(phv+dphv)Ab-phmaxScosα=0

(15)

式中：Aa為上底面的面積/m-2；Ab為下底面的面積/m-2；S為側(cè)表面積/m-2；V為微元體的體積/m-2。

設(shè)水平面上稻谷薄層的半徑為r，則：

Aa=πr2

(16)

Ab=π(r-dycotα)2=π(r2-2rcotαdy)

(17)

(18)

(19)

結(jié)合式(13)～式(18)式并化簡得到錐斗部分最大主應(yīng)力(錐斗壁側(cè)壓應(yīng)力)方程：

(20)

結(jié)合式(10)、式(11)、式(6)得到筒體部分壓應(yīng)力、密度與糧層深度關(guān)系的方程組：

(21)

結(jié)合式(13)、式(14)、式(20)、式(6)得到錐斗部分壓應(yīng)力、密度與糧層深度關(guān)系的方程組：

y∈(H,H+h]

(22)

式中：phv為錐斗中糧堆豎直應(yīng)力/kPa；phmax為錐斗壁側(cè)壓應(yīng)力/kPa。

方程組(21)、方程組(22)的連續(xù)條件為：

pv(H)=phmax(H)

(23)

方程組(21)、方程組(22)的邊界條件為：

(24)

3 帶錐斗筒倉中稻谷堆密度、壓應(yīng)力與糧層深度的關(guān)系

3.1 淮稻5號稻谷在帶錐斗筒倉中的密度與糧層深度的關(guān)系

設(shè)筒倉的參數(shù)條件如表5所示。儲藏稻谷的含水率為12.0% w.b.，其內(nèi)摩擦角為40°，與倉壁的摩擦系數(shù)為0.35，稻谷的表面密度為600 kg/m3。稻谷堆頂部糧面是水平，并選擇糧堆的每一薄層為0.1 m。數(shù)值法求解微分方程組(21)和方程組(22)，從而計算得到不同糧層深度下稻谷堆的密度，如圖5所示。

表5 筒倉的參數(shù)

圖5 稻谷堆密度與糧層深度的關(guān)系

從圖5可以看出：1)1號、2號、3號筒倉中，在筒體部分，隨著糧層深度的增加，稻谷堆密度先快速增大，然后再慢緩增大；到錐斗部分時，隨著糧層深度的增加，稻谷堆密度開始逐漸減小。4號淺圓倉由于倉底沒有帶錐斗，其稻谷堆密度一直隨糧層深度的增加而增大。2)在同樣的糧層深度下，倉內(nèi)稻谷堆的密度隨著筒倉半徑的增大而增大。

3.2 帶錐斗筒倉中糧堆豎直應(yīng)力、側(cè)向應(yīng)力與糧層深度的關(guān)系

筒倉的倉型、大小見表5，稻谷堆的其他參數(shù)同3.1。數(shù)值計算可得到不同糧層深度下稻谷堆的豎直應(yīng)力與側(cè)向應(yīng)力。

3.2.1 豎直應(yīng)力與糧層深度的關(guān)系

計算方程組(21)和方程組(22)給出筒倉中糧堆的豎直應(yīng)力與糧層深度的關(guān)系如圖6 所示。

圖6 稻谷堆所受豎直應(yīng)力與糧層深度的關(guān)系

從圖6可以看出：1)在筒體部分(1、2、3號筒倉)，稻谷堆的豎直應(yīng)力都隨糧層深度的增加而增大，到錐斗部分(1、2、3號筒倉)時，稻谷堆豎直應(yīng)力隨糧層深度的增加而減??；4號淺圓倉底部沒有錐斗，稻谷堆的豎直應(yīng)力一直隨著糧層深度的增加而增大。2)在同樣的糧層深度下，倉內(nèi)稻谷堆的豎直應(yīng)力隨著筒倉半徑的增大而增大。

3.2.2 側(cè)向應(yīng)力與糧層深度的關(guān)系

計算方程組(21)和方程組(22)給出的筒倉中糧堆的側(cè)向應(yīng)力與糧層深度的關(guān)系如圖7 所示。

圖7 稻谷堆所受側(cè)向應(yīng)力與糧層深度的關(guān)系

從圖7可以看出：1)在筒體部分(1、2、3號帶錐斗筒倉中)的稻谷堆的側(cè)向應(yīng)力隨著糧層深度的增加而增大；在筒體與錐斗結(jié)合處，稻谷堆的側(cè)向應(yīng)力突然增加；到了錐斗部分，隨著糧層深度的增加，稻谷堆所受側(cè)向應(yīng)力先稍增大再逐漸減小。4號淺圓倉底部沒有錐斗，所以稻谷堆的側(cè)向應(yīng)力一直隨著糧層深度的增加而增大。2)在同樣的糧層深度下，稻谷堆的側(cè)向應(yīng)力隨著筒倉半徑的增大而增大。

4 基于密度變化機理的筒倉中稻谷儲藏總質(zhì)量

4.1 稻谷儲藏總質(zhì)量積分模型

給定筒倉及稻谷堆的各參數(shù)值，采用本模型計算出帶錐斗筒倉中各糧層的密度，稻谷總質(zhì)量的積分計算為式(25)。

(25)

式中，ρ1為筒體部分糧層深度y處稻谷密度/kg/m3；S1為筒體部分糧層深度y處橫截面積/m2；ρ2為錐斗部分糧層深度y處稻谷密度/kg/m3；S2為錐斗部分糧層深度y處橫截面積/m2。

4.2 模型預(yù)測結(jié)果

國家糧食局科學(xué)研究院小湯山試驗基地1、3、5、7號筒倉內(nèi)徑為3 m，高度為10 m，倉內(nèi)儲藏糧種為稻谷。筒倉與稻谷的參數(shù)及實際賬面稻谷總質(zhì)量見表6。運用式(25)計算得到各個筒倉中儲藏的稻谷總質(zhì)量，計算結(jié)果見表6。

表6 圓倉與稻谷的參數(shù)值及模型預(yù)測的筒倉中稻谷總質(zhì)量

從表6可看出，模型預(yù)測的筒倉中稻谷總質(zhì)量與實際賬面值幾乎一致，最大誤差只有0.55%。

5 模型預(yù)測與體積密度法計算結(jié)果比較

選定表6的參數(shù)，本模型計算得到的4個筒倉中稻谷總質(zhì)量與常用的體積密度法計算得到的4個筒倉中稻谷總質(zhì)量見表7。

表7 本模型、常用的體積密度法得到的稻谷儲藏總質(zhì)量與實際賬面數(shù)的比較

注：體積密度法的密度修正系數(shù)為1.015。

由表7可以看出，1、3、5、7號倉中，本模型計算得到的儲藏總質(zhì)量比常用的體積密度法誤差小，且更接近實際賬面數(shù)，體現(xiàn)出了本模型的可靠性。

6 結(jié)論

6.1 本試驗?zāi)Ｐ皖A(yù)測給出：在筒倉的筒體部分，稻谷堆密度隨著糧層深度的增加而增大，在錐斗部分，隨著糧層深度的增加而減小。由本模型預(yù)測的密度而計算得到的國家糧食局科學(xué)研究院小湯山試驗基地1、3、5、7號筒倉稻谷儲藏總質(zhì)量與實際賬面數(shù)的誤差為0.38，常用的體積密度法計算值與實際賬面數(shù)的誤差為1.06，這體現(xiàn)出了本模型的可靠性。

6.2 本試驗?zāi)Ｐ皖A(yù)測給出：在筒倉的筒體部分，稻谷堆豎直應(yīng)力隨糧層深度的增加而減?。辉阱F斗部分，稻谷堆豎直應(yīng)力隨著糧層深度的增加而減小。

6.3 本試驗?zāi)Ｐ皖A(yù)測給出：在筒倉的筒體部分，稻谷堆側(cè)向應(yīng)力隨著糧層深度的增加而增大；在筒體與錐斗結(jié)合處，稻谷堆側(cè)向應(yīng)力突然增加；到了錐斗部分，隨著糧層深度的增加，稻谷堆側(cè)向應(yīng)力先稍增大再逐漸減小。

[1]王瑞元，朱永義，謝健，等.我國稻谷加工業(yè)現(xiàn)狀與展望[J].糧食與飼料工業(yè)，2011(3)：1-5 Wang R Y，Zhu Y Y，Xie J，et al.Present situation and prospect of rice processing industry in China[J].Cereal and Feed Industry，2011(3)：1-5

[2]張承光.現(xiàn)代化儲糧倉型的選擇[J].糧油倉儲科技通訊，2002(2)：42-44 Zhang C G.The choice of modern storage silos[J].Grain and Oil Storage Science and Technology Communication，2002(2)：42-44

[3]曾丁，黃文彬，華云龍.筒倉壁壓的有限元分析[J].農(nóng)業(yè)工程學(xué)報，1998，14(2)：44-48 Zeng D，Huang W B，Hua Y L.Finite element analysis of silo wall pressure[J].Transactions of the CSAE，1998，14(2)：44-48

[4]石翠霞.筒倉內(nèi)小麥、玉米堆的壓縮特性、倉壁壓應(yīng)力及儲糧總重量的研究[D].南京：南京財經(jīng)大學(xué)，2011 Shi C X.The study on compression characteristics of wheat piles and maize piles，the stress on the wall and the total weight of the grain storage in silos[D].Nanjing：Nanjing University of Finance and Economics，2011

[5]嚴曉婕.稻谷堆與籽粒的壓縮特性及其在筒倉中密度與應(yīng)力分布的研究[D].南京：南京財經(jīng)大學(xué)，2014 Yan X J.The compression properties of paddy heap and paddy grain and study on density and stress distribution of paddy in silos[D].Nanjing：Nanjing University of Finance and Economics，2014

[6]任正曉.糧食庫存檢查實務(wù)[M].北京：中國商業(yè)出版社，2007：15-17 Ren Z X.Grain inventory examination practice[M].Beijing：China Commercial Publishing House，2007：15-17

[7]Mahmoud A A.Loading on shallow cylindrical flexible grain bins[J].Journal of Powder Bulk Solids Technology，1981(5)：12-19

[8]Zhang Q.Static and thermal loads in full size grain bins predicted by vectorized FEM on supercomputer[J].Transactionsof the ASAE，1989，32(2)：685-691

[9]Ruckenbrod C，Eibl J.Numerical results to discharge processes in silos[J]. Reliable flow of particulate Solids. 1993，503-516

[10]Ragneau E，Ooi J Y，Rotter J M. Finite element models for specific applications[J]. Brown C J，Nielsen J，editors. Fundamentals of theory behavior and design.1998,595-608

[11]Goodey R J，Brownb C.J，RotterJM.Predicted patterns of filling pressures in thin-walled square silos[J].Engineering Structures，2006，28：109-119

[12]Loewer J O，Ross I J，Kratzer D D，et al.and Walker，J.N.Properties of ground shelled corn as related to forces in bulk storage structures[J].Transactions of the ASAE，1977，20(1)：155-156

[13]Thompson S A，Ross I J.Compressibility and frictional coefficients of wheat[J].Transactions of ASAE，1983，26(4)：1171-1176

[14]Moya M，Aguado P J，Ayuga F.Mechanical properties of some granular agricultural materials used in silo design[J].International Agrophysics，2013，27：181-193

[15]陳德民，張元.基于壓力傳感器網(wǎng)絡(luò)的糧倉儲糧數(shù)量監(jiān)測系統(tǒng)[J].微計算機信息，2009，25(13)：142-144 Chen D M，Zhang Y.The system of monitoring the number of grain storage based on the pressure sensor network[J].Microcomputer Information，2009，25(13)：142-144

[16]楊平.土力學(xué)[M].北京：機械工業(yè)出版社，2005：98-107 Yang P.Soil mechanics[M].Beijing：Machinery Industry Press，2005：98-107

[17]Rankine W J M.On the stability of loose earth.[J].Philosophical Transactions of the Royal Society Londen，1857，147：1-9

[18]Janssen H A.Versuche über getreidedruck in silozellen[J].Zeitschriff des Vereines Deutscher Ingenieure，1895，39：1045-1049.

The Model of Density and Pressure Distribution of Paddy in the Hopper-Bottom Silo

Du Xiaocui Yan Xiaojie Cheng Xuduo Gao Mengyao Feng Jiachang

(College of Food Science and Engineering;Nanjing University of Finance and Economics;Collaborative Innovation Center for Modern Grain Circulation and Safety Jiangsu Province，Nanjing 210046)

The LHT-1 rebound modulus tester was used to determine the compression density of paddy.The differential equations of density，stress and grain layer depth in hopper-bottom silo were built.Numerical method was used to obtain the relationship of density，pressure and depth and then the total weight of paddy was calculated by integration.Experimental results show that as the average vertical stress of paddy(00.495～245.892 kPa)，the compression density of Huaidao No.5 paddy increased from 582.772～696.593 kg/m3in the moisture range of 10.38%～18.30% w.b.Model calculation results showed that in a silo with a cone hopper，density of paddy in cylinder of silo increased with the paddy depth，the density of paddy in cone hopper decreased with the increase of depth.The vertical stress acting on the paddy in cylinder part of silo increased with depth of paddy，and then decreased with the increase of depth in cone hopper.The lateral stress acting on the paddy in cylinder of silo increased with depth，then increased sharply in the junction of cylinder and cone hopper，the lateral stress of paddy in the cone hopper first increased slightly，and then decreased with the increase of depth of paddy heap.

hopper-bottom silo，paddy，density，vertical stress，lateral pressure

國家自然科學(xué)基金(31371856)，公益性行業(yè)(糧食)科研專項(201313001)

杜小翠，女，1991年出生，碩士，糧食儲藏工程

程緒鐸，男，1957年出生，教授，糧食儲藏工程

TS210.4

A

1003-0174(2017)05-0102-08