• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    溫度和粉碎粒度對不同能量飼料原料熱物理特性的影響

    2019-05-11 07:02:42孔丹丹段恩澤陳計遠王紅英
    農(nóng)業(yè)工程學報 2019年6期
    關(guān)鍵詞:麥麩米糠粉料

    孔丹丹,方 鵬,金 楠,段恩澤,陳計遠,王紅英

    ?

    溫度和粉碎粒度對不同能量飼料原料熱物理特性的影響

    孔丹丹,方 鵬,金 楠,段恩澤,陳計遠,王紅英※

    (中國農(nóng)業(yè)大學工學院,國家農(nóng)產(chǎn)品加工技術(shù)裝備研發(fā)分中心,北京 100083)

    為探究常見能量飼料原料在不同溫度、不同粉碎粒度下的熱物理特性差異,該文以4種谷物(玉米、小麥、大麥和高粱)和4種加工副產(chǎn)品(小麥麩、木薯渣、甜菜渣和米糠)原料為研究對象,分別粉碎過孔徑1.5、2.0和2.5 mm篩片,得到3種不同粉碎粒度的粉料,利用差示掃描量熱法(differential scanning calorimetry,DSC)和線熱源法分別測定了不同粉料在25~100 ℃范圍內(nèi)的比熱和導熱率,通過計算得到相應(yīng)的導溫系數(shù),分析了溫度、粉碎粒度對原料熱特性的影響以及不同原料之間的熱特性差異,并建立了熱特性參數(shù)關(guān)于溫度的回歸預測模型。結(jié)果顯示:8種原料粉碎過1.5 mm篩孔的粉料比熱、導熱率和導溫系數(shù)隨溫度的升高分別在1.580~2.671 kJ/(kg·K),0.054~0.362 W/(m·K)和6.694×10-8~23.254×10-8m2/s范圍內(nèi)變化。整體上,原料的熱特性均隨溫度的升高而呈線性或非線性上升趨勢。不同粉碎粒度的同一原料,在相同溫度下的比熱差異均不顯著(>0.05);隨著粒度的增大,粉料的導熱率和導溫系數(shù)均有逐漸下降的趨勢。4種谷物的比熱與溫度之間呈線性關(guān)系;小麥麩、木薯渣和米糠的比熱可用溫度的三次多項式表示;而甜菜渣的比熱與溫度呈二次關(guān)系。4種谷物和甜菜渣在3種粒度下的導熱率均可用溫度的三次多項式表示,而小麥麩、木薯渣和米糠的導熱率則可用溫度的二次多項式表示。小麥麩和米糠的導溫系數(shù)與溫度呈二次關(guān)系,其余6種原料的導溫系數(shù)則可用溫度的三次多項式表示。研究結(jié)果可為配合顆粒飼料配方變換所需的調(diào)質(zhì)、制粒等熱加工過程的工藝參數(shù)的調(diào)整、優(yōu)化提供理論依據(jù)。

    飼料;比熱;導熱率;導溫系數(shù);溫度;粉碎粒度

    0 引 言

    能量飼料原料是顆粒狀全價配合飼料配方中占比最高、必不可少的組成成分。它們以可溶性碳水化合物、脂肪和油脂等形式提供能量[1]。能量飼料主要包括谷物類、糠麩類、塊根塊莖類、糖蜜、油脂等[2]。其中,玉米在畜禽顆粒飼料配方中占比高達40%~70%[2],是全球消耗量最多的飼料原料(世界上約65%的玉米用作飼料[3])。大麥、小麥和高粱通過熱加工(膨化、制粒)、添加復合酶制劑等方式處理,可以降低原料中的抗營養(yǎng)因子含量,實現(xiàn)部分甚至全部替代畜禽日糧中的玉米[4-6],節(jié)約顆粒飼料的生產(chǎn)成本。甜菜渣、木薯渣等纖維含量較高的能量飼料在反芻動物的日糧中可等量替代7%~15%(日糧配比)的玉米,不僅能降低日糧成本,還能補充有效纖維[7-10]??符燁愒谛笄蓊w粒飼料配方中占比為10%~40%[2]。

    飼料企業(yè)配方設(shè)計師在作多元化、品牌化產(chǎn)品的配方設(shè)計時,會根據(jù)原料的市場價格、供應(yīng)情況等及時變換配方。而生產(chǎn)技術(shù)人員則需對調(diào)質(zhì)、制粒等熱加工過程工藝參數(shù)做出相應(yīng)調(diào)整。但在實際生產(chǎn)中,他們因缺乏對配方中不同原料熱特性的了解,無法準確有效地調(diào)整工藝參數(shù),導致更換配方后首次生產(chǎn)時會出現(xiàn)開機試制時間延長、生產(chǎn)率降低、顆粒品質(zhì)穩(wěn)定性變差等問題,嚴重影響飼料企業(yè)個性化產(chǎn)品的開發(fā)質(zhì)量及效益。因此,研究不同種類能量飼料的熱特性差異,掌握原料在熱處理過程中的物性變化規(guī)律,從而合理設(shè)計熱加工過程的工藝參數(shù),穩(wěn)定顆粒產(chǎn)品質(zhì)量變得尤為重要。

    比熱、導熱率和導溫系數(shù)作為3個最重要的熱物理特性參數(shù),廣泛應(yīng)用于涉及熱傳遞的熱加工過程的設(shè)計、優(yōu)化和建模模擬計算[11-12]。目前,已有許多學者針對糧食谷粒(面粉)[13-18]、果核種子[11,19-22]、果蔬類[12,23-27]食品加工領(lǐng)域物料的熱特性進行了許多深入的研究,而針對飼料加工領(lǐng)域原料熱特性的研究則相對缺乏。王紅英等[28-29]研究了飼用玉米、主要飼料原料不同配比的混合料的比熱??椎ささ萚30-31]報道了溫度、含水率、粉碎粒度以及乳清粉含量對仔豬配合粉料比熱的影響規(guī)律。針對不同種類能量飼料熱特性差異的系統(tǒng)研究還沒未見報道。

    本文以畜禽顆粒料常用的8種能量飼料(4種谷物:玉米、小麥、大麥和高粱;4種加工副產(chǎn)品:小麥麩、木薯渣、甜菜渣和米糠)為研究對象,粉碎過孔徑1.5、2.0和2.5 mm篩片,得到3種不同粒度的粉料,利用差示掃描量熱法和線熱源法分別測定了所有粉料在25~100 ℃范圍內(nèi)的比熱和導熱率,并計算出導溫系數(shù),建立了3個熱特性參數(shù)關(guān)于溫度的可靠預測模型;分析了原料不同種類、不同粒度之間的熱特性差異,以期為顆粒飼料配方變換所需的熱加工工藝參數(shù)的調(diào)整、優(yōu)化提供理論依據(jù)。

    1 材料與方法

    1.1 材料與儀器

    1.1.1 試驗材料

    玉米:鄭單958,源自中國農(nóng)業(yè)大學涿州試驗基地;小麥、大麥、高粱、木薯渣顆粒(泰國進口)和甜菜渣顆粒:取自北京首農(nóng)畜牧發(fā)展有限公司飼料分公司;小麥麩:北京古船食品有限公司;米糠:北京新成谷源糧油工貿(mào)有限公司。

    1.1.2 儀器設(shè)備

    15B萬能粉碎機:江陰市宏達粉體設(shè)備有限公司;DHG-9240A電熱恒溫鼓風干燥箱:上海精宏實驗設(shè)備有限公司;AL204分析天平:梅特勒-托利多儀器有限公司; OCTAGON 200標準振篩機和不銹鋼標準篩(ISO3310):英國Endecotts(恩德)公司;DSC 214 Polyma差示掃描量熱儀:德國Netzsch(耐馳)公司;KD2 Pro熱特性分析儀:美國Decagon devices公司;HH-S數(shù)顯恒溫油浴鍋:常州國華電器有限公司。

    1.2 試驗方法

    1.2.1 樣品處理

    將采集到的除米糠外的7種飼料原料,用配有孔徑為1.5、2.0和2.5 mm篩片的粉碎機粉碎,得到3種不同粒度的粉料。米糠可100%通過孔徑1.5 mm的篩片,故無需粉碎。

    由于物料含水率對熱物理特性的影響較大,本研究將所有粉料的含水率統(tǒng)一調(diào)整到12%(濕基)。若初始含水率較高,則通過30℃低溫烘干的方式將其含水率降到12%。若初始含水率較低,則通過往已知質(zhì)量的粉料中添加計算好的蒸餾水量,置于混合機中充分混合15 min后裝入密封袋,在4 ℃的冷藏柜靜置1周,期間每隔6 h晃動1 min,使其水分達到平衡[15]。處理后各原料的水分含量如表1所示。

    1.2.2 含水率及粉碎粒度測定

    玉米的含水率按照GB/T 10362-2008進行測定,小麥、大米和高粱的含水率則依照GB 5009.3-2016進行測定,其余加工副產(chǎn)品的含水率依照GB/T 6435-2014進行測定。粉料的粉碎粒度依照ANSI/ASAE S319.4-2008標準的十四層篩法進行篩分、稱量篩上物質(zhì)量,按式(1)進行計算,以幾何平均直徑表示。

    表1 經(jīng)干燥或加濕處理后能量飼料的實測含水率

    注:結(jié)果表示為平均數(shù)±標準差;<1.5/2.0/2.5 mm篩孔指飼料粉碎過孔徑為1.5/2.0/2.5 mm的篩片,下同。

    Note: Results expressed as means ± standard deviation. <1.5/2.0/2.5 mm mesh size means feed ingredients being ground passing through screen with aperture size of 1.5/2.0/2.5 mm, the same below.

    1.2.3 比熱測定

    本試驗采用DSC掃描了不同粉碎粒度的8種能量飼料在5~120 ℃范圍內(nèi)熱譜曲線,并計算其在不同溫度下的比熱值。試驗開始前,先采用C10H16、In、Sn、Bi、Zn、CsCl共6種標準物質(zhì)對儀器進行溫度校正和靈敏度校正。所有測定過程通入氮氣(流量:40 mL/min)對樣品進行保護,裝樣品的坩堝均加蓋密封,以避免測定過程樣品水分蒸發(fā)。具體測定方法同孔丹丹等[30]描述的一致。每個樣品重復測定3次,取其平均值作為最終結(jié)果。

    1.2.4 導熱率測定

    導熱率采用KD2 Pro探針測定,它是利用非穩(wěn)態(tài)熱傳導的線熱源原理制成的。其測定原理基于埋入無限大均勻介質(zhì)中的線熱源(半徑大小忽略),最初是等溫的,與其周圍的環(huán)境處于平衡狀態(tài),介質(zhì)中線熱源的熱流可以用如下的溫度傳導方程來表示:

    其中為樣品在處的溫度,℃;為距離線熱源的半徑,m;為加熱時間,s;為導溫系數(shù),m2/s。該溫度的解由Hopper和Lepper[32]給出:

    而1與2之間的溫度變化為

    因此

    即只要測出121和2的值,便可計算出導熱率k。

    將封口后的玻璃杯用平坦的重物(塑料塊)壓住,以防止升溫過程中體積膨大。而后將其置于恒溫油浴鍋中(水浴鍋在高溫時大量蒸發(fā)的水蒸氣不利于試驗測定),在預定的溫度加熱1 h左右使其平衡(加熱時間由先前的預試驗獲得)。將整個探針垂直插入物料的中部,固定住,等待5 min,當探針溫度達到與物料溫度一致時,開始測定。每種粉料每個處理溫度進行3次重復測定,取3次的平均值為導熱率測定的結(jié)果。

    1.2.5 導溫系數(shù)測定

    粉料的導溫系數(shù)采用已經(jīng)測得的導熱率、密度和比熱的值計算得到:

    其中ρ為堆積密度,kg/m3;C為比熱,J/(kg·K)。

    1.3 數(shù)據(jù)分析

    利用SPSS 22.0軟件對試驗數(shù)據(jù)進行單因素方差分析和Duncan檢驗(<0.05,方差不齊性時采用Dunnett’s T3檢驗),用OriginPro 9.1軟件繪制圖形。能量飼料熱特性與溫度的最佳關(guān)系采用SPSS 22.0軟件的線性與非線性逐步回歸來確定,剔除回歸系數(shù)不顯著的項,只保留顯著項。選取決定系數(shù)(2)最高、均方根誤差(root mean square error,RMSE)和平均相對百分誤差(mean relative percent error,)最小的為最優(yōu)模型。RMSE和分別由式(8)、式(9)計算得到。

    其中exp表示實測值,pred表示預測值,表示試驗數(shù)據(jù)的個數(shù)。

    2 結(jié)果與分析

    2.1 能量飼料粉碎粒度差異分析

    玉米、小麥等7種飼料粉碎過孔徑1.5、2.0和2.5 mm篩片后的粉碎粒度以及未經(jīng)粉碎的米糠粒度差異如表2所示。同一飼料原料,隨著粉碎篩孔的增大,粉碎粒度顯著增加(<0.05)。粉碎過同一篩孔的不同飼料原料之間,粉碎粒度差異較大。以過2.0 mm篩孔的粉料為例,7種飼料粉碎粒度在325.65~594.83m范圍內(nèi)變化,且不同原料之間粒度均表現(xiàn)出顯著的差異(<0.05)。其中,小麥麩的粒徑最大,甜菜渣顆粒次之,木薯渣顆粒最??;4種谷物中,玉米和小麥的粒徑要遠小于大麥和高粱的。這是由于不同種類的能量飼料,其顆粒大小、形狀、硬度、質(zhì)地等差異較大,導致其粉碎特性也表現(xiàn)出較大差異。

    表2 能量飼料粉碎過不同篩孔后的粒度

    注:結(jié)果表示為平均數(shù)±標準差,同列肩標(大寫字母)或同行肩標(小寫字母)不同字母表示差異顯著(<0.05),不同。

    Note: Results are expressed as mean ± standard deviation. Values not sharing a common capital/lowercase letter in the same column/row are significantly different (<0.05), the same below.

    2.2 能量飼料比熱差異分析

    8種粉碎過1.5 mm篩孔的飼料粉料的比熱隨溫度的變化如圖1所示。除米糠外,其余7種飼料的比熱均隨溫度的升高而增大,且25~100 ℃范圍內(nèi),玉米、小麥、大麥和高粱4種谷物飼料的比熱分別在1.678±0.055~2.421±0.083、1.644±0.027~2.402±0.025、1.642±0.018~2.385±0.047和1.649±0.026~2.441±0.044 kJ/(kg·K)范圍內(nèi)變化,小麥麩、木薯渣和甜菜渣3種加工副產(chǎn)品飼料的比熱分別在1.621±0.030~2.540±0.009、1.580±0.010~2.195±0.011和1.746±0.026~2.351±0.027 kJ/(kg·K)的范圍內(nèi)變化。溫度平均每升高1 ℃,7種原料的比熱值分別升高0.010、0.010、0.010、0.011、0.012、0.008和0.008 kJ/(kg·K)。與其他7種原料不同的是,米糠的比熱在25.0~33.5 ℃范圍內(nèi),呈現(xiàn)先升高后降低的趨勢,并在28.0 ℃時出現(xiàn)比熱的峰值點,為1.985 kJ/(kg·K);33.5~100 ℃范圍內(nèi),其比熱隨溫度的升高而在1.868~2.671 kJ/(kg·K)范圍內(nèi)增大。比熱峰值的出現(xiàn)可能與升溫時米糠中含量較高的油脂(≥14%)在脂肪酶、脂肪氧合酶的催化作用下在發(fā)生水解[33-34](吸熱過程)、氧化酸敗[35]有關(guān)。米糠中的油脂水解生成游離脂肪酸,其中大多為不飽和脂肪酸(含量達80%以上[36]),容易發(fā)生氧化酸敗。上述8種能量飼料的比熱值范圍與文獻中報道的低含水率的谷物[14,17]、面粉[18]、油菜籽[20]、枯茗籽[22]、仔豬配合粉料[30-31]等農(nóng)業(yè)物料相接近,但遠低于高含水率的果汁[12]、葉類蔬菜[26]和白蘿卜[27]等果蔬類物料??椎ささ萚30]的研究顯示,12%含水率時,仔豬配合粉料的比熱隨溫度的升高(25~100 ℃)而在1.679~2.503 kJ/(kg·K)范圍增大,這與本文中的玉米、小麥和高粱的比熱范圍較為接近。

    圖1 8種能量飼料粉料(1.5 mm篩孔)比熱隨溫度的變化

    以高粱和甜菜渣為例,同種原料粉碎過3種篩孔的粉料比熱隨溫度的變化規(guī)律如圖2所示??梢钥闯?,雖然粉料的粉碎粒度差異顯著,但其比熱隨溫度變化的規(guī)律均一致,且比熱值差異較小。方差分析顯示,不同粉碎粒度的7種能量飼料(米糠除外),在同一溫度下的比熱差異均不顯著(>0.05),表明,粉碎粒度對上述7種能量飼料的比熱無顯著影響。這與孔丹丹等[30]有關(guān)粉碎粒度對低含水率(10%~16%)仔豬配合粉料的比熱影響不顯著的研究結(jié)果相一致。

    圖2 不同粉碎粒度的高粱和甜菜渣比熱隨溫度的變化

    由于粉碎粒度對飼料原料的比熱影響不顯著,故本文以3種粒度的數(shù)據(jù)為基礎(chǔ),分別建立了8種飼料原料比熱關(guān)于溫度的回歸預測模型(見表3)。模型的所有系數(shù)均顯著,2≥0.991,且RMSE≤0.024 kJ/(kg·K),≤0.951%,說明模型擬合度較佳,充分反映了響應(yīng)與因子之間的關(guān)系。模型顯示:在25~100 ℃的范圍內(nèi),4種谷物的比熱與溫度均呈線性正相關(guān)關(guān)系;小麥麩、木薯渣和米糠3種加工副產(chǎn)品的比熱與溫度的關(guān)系均可以用溫度的三次多項式表示;而甜菜渣的比熱隨溫度變化遵循二次多項式的關(guān)系。

    以大麥、小麥麩和甜菜渣的實測比熱數(shù)據(jù)為例,對其相應(yīng)的線性、三次和二次模型進行驗證,結(jié)果如式(10)~(12)所示。可以看出,預測值()和實測值()呈良好的一階線性關(guān)系(2≥0.995)。斜率越接近于1,截距越接近于0,模型的預測效果越好。驗證結(jié)果表明,本文建立的模型可以較為真實地描述原料比熱與溫度的關(guān)系。

    許多研究報道了農(nóng)業(yè)物料的比熱與溫度的線性正相關(guān)關(guān)系。Kaletun?[18]的研究顯示小麥、玉米和大米3種谷物面粉的比熱表現(xiàn)出基于溫度(20~110 ℃)的線性依存關(guān)系。Obot等[27]報道了20~80 ℃范圍內(nèi),白蘿卜比熱隨溫度的升高而線性增加。Aviara等[37]測定的乳木果果核比熱與溫度(30~90 ℃)呈良好的線性關(guān)系。還有研究報道了麻風樹種仁[11]、雙低油菜籽[21]、弗羅里達側(cè)耳[25]、玉米[31]、阿月渾子果仁[38]的比熱與溫度的線性正相關(guān)關(guān)系。這些均與本文中4種谷物的比熱和溫度關(guān)系的研究結(jié)果相一致。

    表3 8種能量飼料比熱的回歸模型及統(tǒng)計信息

    注:C表示比熱,kJ·(kg–1·K–1);表示溫度,25~100 ℃范圍內(nèi),下同。

    Note:Cmeans specific heat, kJ·(kg–1·K–1);means temperature, in the range of 25-100 ℃, the same below.

    此外,還有少數(shù)研究報道了農(nóng)業(yè)物料比熱關(guān)于溫度的非線性關(guān)系。Shrestha等[19]報道了王不留行籽比熱關(guān)于溫度(25~80 ℃)的二次關(guān)系。Yu等[20]建立了雙低油菜籽的比熱在40~90 ℃范圍內(nèi)關(guān)于溫度的二次模型。Singh等[22]的研究表明,低干基含水率(1.8%~11.1%)的枯茗籽比熱可以用溫度(-70~50 ℃)的二次多項式表示。還有研究報道了魚粉[31]、琉璃苣種子[39]比熱與溫度的二次關(guān)系。這些物料比熱與溫度關(guān)系均和本文中甜菜渣顆粒的研究結(jié)果相符??椎ささ萚31]建立了乳清粉比熱關(guān)于溫度(25~110 ℃)的三次模型,這與本文中小麥麩、木薯渣和米糠的研究結(jié)果相似。

    比熱表示單位質(zhì)量的物質(zhì)每升高(降低)1 ℃所吸收(放出)的熱量,故能量飼料比熱的高低直接關(guān)系到調(diào)質(zhì)升溫過程中每噸飼料所需供給的蒸汽量的多少。目前,飼料廠加工畜禽顆粒料的常見調(diào)質(zhì)溫度范圍為65~85 ℃,所以本文基于已建立的8種原料的比熱預測模型,分別積分求解出每噸原料從室溫25 ℃升至65、75和85 ℃時理論上所需吸收的熱量(見表4)。

    整體來看,溫度越高,所需吸收的能量越多;原料從25 ℃升至75和85 ℃所需的熱量大約是25 ℃升至65 ℃的1.28倍和1.58倍。同一升溫條件下,8種原料中,木薯渣所需的熱量最低(75 ℃時,為89 755.219 kJ/t),米糠所需的熱量最高(75 ℃時,為103 191.229 kJ/t),后者所需熱量約為前者的1.15倍。4種谷物中,小麥、大麥和高粱升高同等溫度所需吸收的熱量均較玉米稍低,分別約為玉米的0.98、0.98和0.99倍。4種加工副產(chǎn)品與玉米相比,同一升溫范圍,小麥麩所需吸收的熱量與玉米相接近,而木薯渣為玉米的0.94倍;甜菜渣和米糠所需的熱量均較玉米高,約為玉米的1.04和1.08倍。因此,在用其他能量飼料部分替代玉米更換配方或變換調(diào)質(zhì)溫度時,需根據(jù)上述吸熱規(guī)律,及時調(diào)整調(diào)質(zhì)過程中蒸汽的添加量,以確保快速、精確地達到預期的調(diào)質(zhì)效果。

    表4 能量飼料升溫過程所需吸收熱量

    2.3 能量飼料導熱率差異分析

    8種過1.5 mm篩孔的能量飼料導熱率隨溫度的變化如圖3所示。顯而易見,在25~100 ℃溫度范圍內(nèi),所有飼料的導熱率均隨溫度的升高呈非線性上升的趨勢,這與許多農(nóng)業(yè)物料導熱率的性質(zhì)相似[19-20,22]。這是由于,物料中的離子和偶極子在高溫下會表現(xiàn)出較高的晶格振動,從而加快傳熱的速度[20]。而在55~100 ℃高溫時,能量飼料會發(fā)生淀粉糊化、非酶褐變反應(yīng)等,導致其流變特性發(fā)生變化,黏度增大,顆粒腫脹,孔隙率降低,導熱率急劇升高,這可能是導致曲線偏離線性行為的主要原因。

    圖3 8種能量飼料(1.5 mm篩孔)導熱率隨溫度的變化曲線

    玉米、小麥、大麥和高粱4種谷物的導熱率隨溫度的升高分別在0.089~0.299、0.078~0.336、0.088~0.276和0.082~0.288 W/(m·K)的范圍內(nèi)變化;小麥麩、木薯渣、甜菜渣和米糠4種加工副產(chǎn)物的導熱率則分別在0.054~0.190、0.092~0.362、0.085~0.263和0.064~0.202 W/(m·K)的范圍內(nèi)變化。方差分析顯示,不同能量飼料之間,導熱率差異較大,表現(xiàn)為:常溫25 ℃時,除玉米和大麥之間導熱率差異不顯著外(>0.05),其余原料之間導熱率差異均顯著(<0.05),且木薯渣導熱率最高,玉米和大麥次之,小麥麩導熱率最低,米糠次之;55 ℃時,除玉米和甜菜渣、小麥和大麥之間導熱率差異不顯著(>0.05)外,其余原料之間導熱率差異均顯著(<0.05);高溫85和100 ℃時,各原料之間,導熱率差異均顯著(<0.05),木薯渣的導熱率最高,小麥次之,而小麥麩的最低,米糠次之。在整個研究溫度范圍內(nèi),4種加工副產(chǎn)品之間導熱率差異均顯著(<0.05)。對于散粒材料而言,結(jié)構(gòu)越緊實,空隙率越少,傳熱效果越好[19-20]。已有的研究也顯示王不留行籽[19]和雙低油菜籽[20]在堆積密度下的導熱率明顯低于在粒子密度下的。小麥麩和米糠因質(zhì)地松軟,孔隙率高,其堆積密度顯著低于其余6種原料(見表5),導致其導熱率也較其余幾種原料低。8種原料中,高粱和小麥的堆積密度較高(分別為0.728和0.722 g/cm3),玉米和木薯渣的次之(分別為0.714和0.711 g/cm3)。木薯渣較高的導熱率可能源于其較高堆積密度和較小的平均粒徑。高溫時,各種物料之間顯著的導熱率差異可能源自其淀粉糊化、非酶褐變的程度不同而導致的化學組分、物質(zhì)狀態(tài)的巨大差異。

    表5 8種不同粉碎粒度的能量飼料的堆積密度

    注:粒度值參照表2,下同。

    Note: The particle sizes refer to Table 2, the same below.

    Yu等[20]的研究表明11%含水率的雙低油菜籽的導熱率隨溫度(40~90 ℃)的增加而在0.07~0.13 W/(m·K)范圍內(nèi)升高。據(jù)Shrestha等[19]建立的預測模型可得出12%含水率(濕基)的王不留行籽在25~55 ℃下的導熱率在0.060~0.071 W/(m·K)范圍內(nèi)變化。上述2種物料低溫時的導熱率與本文中小麥麩、米糠的相接近,而高溫時則明顯較8種原料的低。據(jù)Jian等[40]建立的高油含量的雙低油菜籽導熱率的預測模型可計算出其在12%含水率(濕基)時,導熱率隨溫度升高(-10~30 ℃)而在0.083~0.101 W/(m·K)范圍內(nèi)線性增大,同等溫度下,其導熱率較本文中的8種原料稍高。

    圖4a、4b和4c分別顯示了25、55和85 ℃時7種原料(除米糠)在3種粒度下的導熱率差異。從圖中可以看出,大體上,隨著粒度的增大,粉料的導熱率有逐漸下降的趨勢,且這種差異隨著溫度的升高而逐漸增大。

    常溫25 ℃時,過1.5和2.0 mm篩孔的玉米、小麥和高粱導熱率顯著高于過2.5 mm的(<0.05),大麥不同粒度間導熱率差異顯著(<0.05),而不同粒度的小麥麩、木薯渣和甜菜渣差異則不顯著(>0.05);55 ℃時,玉米不同粒度間導熱率差異顯著,過1.5和2.0 mm篩孔小麥和大麥導熱率顯著高于過2.5 mm的(<0.05),高粱、木薯渣和甜菜渣則表現(xiàn)為1.5 mm的顯著高于2.0和2.5 mm的(<0.05),小麥麩不同粒度間差異仍不顯著(>0.05);高溫85 ℃時,7種原料不同粒度之間導熱率差異均顯著(<0.05)。這主要是由于高溫時,粒度較細的粉料更容易發(fā)生糊化變性,散粒體結(jié)構(gòu)變得更為緊密,因而導熱率也較粗粉高。

    不同粉碎粒度的8種原料導熱率關(guān)于溫度的回歸模型及統(tǒng)計信息如表6所示。回歸模型的系數(shù)均顯著,2≥0.980,<0.001,RMSE≤0.008 W/(m·K),≤5.546%,表明模型擬合度良好,預測精度高。從表中可以看出,4種谷物和甜菜渣在3種粒度下的導熱率都可以用溫度的三次多項式來表示,而小麥麩、木薯渣和米糠的導熱率與溫度遵循二次多項式的關(guān)系。

    a. 25 ℃b. 55 ℃c. 85 ℃

    表6 不同粉碎粒度的能量飼料導熱率的回歸模型及統(tǒng)計信息

    注:表示導熱率,W·(m–1·K–1)。

    Note:means thermal conductivity, W·(m–1·K–1).

    以高粱(<2.0 mm)、木薯渣(<2.0 mm)和米糠的實測導熱率數(shù)據(jù)為例,對其相應(yīng)的三次和二次模型進行驗證,結(jié)果如式(13)~(15)所示??梢钥闯?,預測值()和實測值()呈良好的一階線性關(guān)系(2≥0.992),且斜率接近于1,截距約等于0,表明本文建立的模型可以對原料的導熱率與溫度的關(guān)系作較為真實的表達。

    許多文獻報道了物料導熱率與溫度的非線性關(guān)系。Singh等[22]的研究表明,枯茗籽的導熱率與溫度(-50~50 ℃)之間呈二次關(guān)系。Shrestha等[19]報道了堆積密度下王不留行籽的導熱率是關(guān)于溫度(25~55℃)的二次多項式。這些均與本文小麥麩、木薯渣和米糠導熱率與溫度的研究結(jié)果相一致。Yu等[20]的研究揭示了雙低油菜籽在堆積密度下的導熱率與溫度(40~90 ℃)之間呈三次關(guān)系。這與本研究中4種谷物和甜菜渣的導熱率變化規(guī)律相似。此外,也有研究報道了導熱率與溫度之間的線性關(guān)系。Shrivastava等[25]的研究顯示,弗羅里達側(cè)耳的導熱率與溫度(40~70 ℃)線性相關(guān)。Obot等[27]用簡易構(gòu)造的探針測定了白蘿卜的導熱率,顯示其在20~80 ℃范圍內(nèi)與溫度呈良好的線性正相關(guān)關(guān)系。Jian等[40]建立了基于貯藏溫度范圍(-10~30 ℃)用于預測高油含量的雙低油菜籽導熱率的線性回歸方程。

    2.4 能量飼料導溫系數(shù)差異分析

    8種過1.5 mm篩孔的能量飼料導溫系數(shù)隨溫度的變化如圖5所示。25~100 ℃范圍內(nèi),玉米、小麥、大麥和高粱4種谷物的導溫系數(shù)分別在7.459×10-8~17.320×10-8、6.694×10-8~19.195×10-8、7.878×10-8~17.045×10-8和6.857×10-8~16.109× 10-8m2/s范圍內(nèi)變化,小麥麩、木薯渣、甜菜渣和米糠4種加工副產(chǎn)品的導溫系數(shù)分別在9.543×10-8~20.892×10-8、8.217×10-8~23.254×10-8、7.258×10-8~16.646×10-8和8.092×10-8~18.718 ×10-8m2/s的范圍內(nèi)變化。大體上(除高粱外),原料的導溫系數(shù)呈現(xiàn)隨溫度在25~55 ℃范圍內(nèi)緩慢增加,55~100 ℃范圍內(nèi)急劇上升的趨勢;高粱的導溫系數(shù)在40~55 ℃范圍內(nèi)稍有下降(從7.556×10-8降至6.989×10-8m2/s)。導溫系數(shù)隨溫度的增加主要歸功于導熱率隨溫度的增加速率高于比熱的增加速率[20],導熱率的變化對物料的導溫系數(shù)有著顯著的影響[19]。

    圖5 8種能量飼料(1.5 mm篩孔)導溫系數(shù)隨溫度的變化曲線

    導溫系數(shù)反映了物料導熱能力與儲熱能力之間的比值關(guān)系,比熱與堆積密度的乘積稱作體積熱容量,它表示物料儲熱的能力。在55~85 ℃(調(diào)質(zhì)溫度)范圍內(nèi),小麥麩、米糠和木薯渣的導溫系數(shù)位居8種原料的前三,明顯高于4種谷物飼料,而甜菜渣的導溫系數(shù)則與玉米、小麥和大麥3種谷物相接近,高粱明顯較上述3種谷物低。小麥麩、米糠的導熱率較低,儲熱能力也較差,故導溫系數(shù)的值較高。而木薯渣較高的導溫系數(shù)則是由其較高的導熱率和較低比熱值導致的,說明其容易被加熱和冷卻。而高粱較低的導溫系數(shù)則源于其較低的導熱率和較高的體積熱容量,說明其導熱能力較弱,而儲熱能力較強,不易升溫和冷卻,與文獻報道的大豆性質(zhì)類似[17]。

    圖6顯示了溫度為70 ℃時7種原料在3種粒度下導溫系數(shù)的差異。隨著粒度的增大,7種原料的導溫系數(shù)均有下降的趨勢,小麥麩(14.683×10-8~13.175×10-8m2/s)、大麥(9.847×10-8~8.374×10-8m2/s)和小麥(9.086×10-8~7.910×10-8m2/s)下降最為明顯,而高粱(7.852×10-8~7.640×10-8m2/s)則表現(xiàn)不明顯。

    圖6 70 ℃時3種粉碎粒度的7種能量飼料導溫系數(shù)的差異

    8種過1.5 mm篩孔的原料導溫系數(shù)關(guān)于溫度的回歸模型及統(tǒng)計信息如表7所示。顯然,除小麥麩和米糠外,其余6種原料的導溫系數(shù)隨溫度的變化均遵循三次多項式的關(guān)系,而小麥麩和米糠的導溫系數(shù)與溫度則可以用更為簡單的二次關(guān)系來描述。以上回歸模型的系數(shù)均顯著,2≥0.987,<0.000 1,RMSE≤0.488×10-8m2/s,≤3.224%,表明模型擬合度良好,可以直接對8種原料導溫系數(shù)與溫度的關(guān)系作較為準確的預測。Shrestha等[19]的研究報道了王不留行籽導溫系數(shù)與溫度(25~55 ℃)的二次關(guān)系。Singh等[22]的研究表明枯茗籽的導溫系數(shù)在-50~50 ℃范圍隨溫度的升高而增大,并可用溫度的二次多項式表示。上述2種物料導溫系數(shù)關(guān)于溫度的依存關(guān)系與本文中小麥麩、米糠的相似。Yu等[20]建立的預測模型顯示雙低油菜籽在堆積密度下導溫系數(shù)與溫度(40~90 ℃)呈三次關(guān)系。這與本文中4種谷物、木薯渣和甜菜渣的研究結(jié)果相符。此外,還有研究報道了白蘿卜導溫系數(shù)與溫度的(20~80 ℃)線性正相關(guān)關(guān)系。

    表7 能量飼料原料(1.5 mm 篩孔)導溫系數(shù)的回歸模型及統(tǒng)計信息

    注:表示導溫系數(shù),×10-8·m2·s–1。

    Note:means thermal diffusivity, ×10-8·m2·s–1.

    3 結(jié) 論

    本文測定了粉碎過孔徑1.5、2.0和2.5 mm篩片的4種谷物(玉米、小麥、大麥和高粱)和4種加工副產(chǎn)品(小麥麩、木薯渣、甜菜渣和米糠)在25~100 ℃范圍內(nèi)的比熱、導熱率,并計算出相應(yīng)的導溫系數(shù);分析了不同原料、不同粒度粉料的熱特性差異。主要結(jié)論如下:

    1)粉碎過1.5 mm篩孔的4種谷物比熱隨溫度的升高分別在1.678~2.421、1.644~2.402、1.642~2.385和1.649~2.441 kJ/(kg·K)范圍內(nèi)增大,小麥麩、木薯渣和甜菜渣比熱則分別在1.621~2.540、1.580~2.195和1.746~2.351 kJ/(kg·K)的范圍內(nèi)增大。米糠在28.0 ℃時出現(xiàn)比熱的峰值點,為1.985 kJ/(kg·K);33.5~100 ℃范圍內(nèi),其比熱隨溫度的升高而在1.868~2.671 kJ/(kg·K)范圍內(nèi)增大。不同粉碎粒度的同種原料,在相同溫度下的比熱差異均不顯著(>0.05)。4種谷物的比熱與溫度之間均呈線性關(guān)系;小麥麩、木薯渣和米糠的比熱可用溫度的三次多項式表示;甜菜渣的比熱與溫度呈二次關(guān)系。

    2)同一升溫范圍,小麥、大麥、高粱和木薯渣所需吸收的熱量均較玉米稍低,分別約為玉米的0.98、0.98、0.99和0.94倍;小麥麩所需吸收的熱量與玉米相接近;甜菜渣和米糠所需的熱量較高,約為玉米的1.04和1.08倍。

    3)粉碎過1.5 mm篩孔的4種谷物導熱率隨溫度的升高分別在0.089~0.299、0.078~0.336、0.088~0.276和0.082~0.288 W/(m·K)的范圍內(nèi)增大;4種加工副產(chǎn)品的導熱率則分別在0.054~0.190、0.092~0.362、0.085~0.263和0.064~0.202 W/(m·K)的范圍內(nèi)增加。高溫(85和100 ℃)時,各原料之間導熱率差異均顯著(>0.05),木薯渣導熱率最高,小麥次之,而小麥麩導熱率最低,米糠次之。隨著粒度的增大,粉料的導熱率有逐漸下降的趨勢,且這種差異隨著溫度的升高而逐漸增大。4種谷物和甜菜渣在3種粒度下的導熱率都可以用溫度的三次多項式來表示,而小麥麩、木薯渣和米糠的導熱率與溫度遵循二次多項式的關(guān)系。

    4)粉碎過1.5 mm篩孔的4種谷物的導溫系數(shù)隨溫度的升高分別在7.459×10-8~17.320×10-8、6.694×10-8~19.195×10-8、7.878×10-8~17.045×10-8和6.857×10-8~16.109×10-8m2/s范圍內(nèi)變化,4種加工副產(chǎn)品的導溫系數(shù)分別在9.543×10-8~20.892×10-8、8.217×10-8~23.254×10-8、7.258×10-8~16.646×10-8和8.092×10-8~18.718×10-8m2/s的范圍內(nèi)增大。大體上,原料的導溫系數(shù)呈現(xiàn)隨溫度在25~55 ℃范圍內(nèi)緩慢增加,55~100 ℃范圍內(nèi)急劇上升的趨勢。調(diào)質(zhì)溫度55~85 ℃范圍內(nèi),小麥麩、米糠和木薯渣的導溫系數(shù)位居8種原料的前三,明顯高于4種谷物飼料。隨著粒度的增大,原料的導溫系數(shù)有下降的趨勢。小麥麩和米糠的導溫系數(shù)與溫度呈二次關(guān)系,其余6種原料的導溫系數(shù)則可用溫度的三次多項式表示。

    [1] Kellems R O,Church D C. 畜禽飼料與飼養(yǎng)學:第5版[M]. 姜成鋼,張輝譯. 北京:中國農(nóng)業(yè)大學出版社,2006.

    [2] 曹康,郝波. 中國現(xiàn)代飼料工程學:上卷[M]. 上海:上??茖W技術(shù)文獻出版社,2014.

    [3] 李德發(fā),宋國隆,趙麗丹. 飼料工業(yè)對玉米的數(shù)量需求和質(zhì)量要求[J]. 玉米科學,2003,11(1):83-87.

    [4] 公衍玲,徐銀龍,鄒新華. 大麥、高粱在養(yǎng)豬生產(chǎn)中的應(yīng)用[J]. 豬業(yè)科學,2016,33(4):80-82.

    [5] Torres K A A, Pizauro J M, Soares C P, et al. Effects of corn replacement by sorghum in broiler diets on performance and intestinal mucosa integrity[J]. Poultry Science, 2013, 92(6): 1564-1571.

    [6] 李建軍,雷小文,邱光忠,等. 不同能量飼料配方對肥育豬生長性能的影響研究[J]. 豬業(yè)科學,2017,34(6):88-89.

    [7] 樊懿萱,王鋒,王強,等. 發(fā)酵木薯渣替代部分玉米對湖羊生長性能、血清生化指標、屠宰性能和肉品質(zhì)的影響[J]. 草業(yè)學報,2017,26(3):91-99. Fan Yixuan, Wang Feng, Wang Qiang, et al. Effects of substitution of corn and fermented cassava residue in the diet on growth, carcass characteristics and meat quality of Hu lamb[J]. Acta Prataculturae Sinica, 2017, 26(3): 91-99. (in Chinese with English abstract)

    [8] 李景偉. 木薯渣對肉牛生產(chǎn)性能、屠宰性能、胴體品質(zhì)及血清生化指標的影響[D]. 泰安:山東農(nóng)業(yè)大學,2015.

    [9] Guo Yongqing, Xu Xiaofeng, Zou Yang, et al. Changes in feed intake, nutrient digestion, plasma metabolites, and oxidative stress parameters in dairy cows with subacute ruminal acidosis and its regulation with pelleted beet pulp[J]. Journal of Animal Science and Biotechnology, 2013, 4(1): 1-10.

    [10] 郭勇慶,劉潔,柏兆海,等. 甜菜顆粒粕在反芻動物日糧中的應(yīng)用[J]. 中國飼料,2015(22):10-15.

    [11] Sirisomboon P, Posom J. Thermal properties ofkernels[J]. Biosystems Engineering, 2012, 113(4): 402-409.

    [12] Mercali G D, Sarkis J R, Jaeschke D P, et al. Physical properties of acerola and blueberry pulps[J]. Journal of Food Engineering, 2011, 106(4): 283-289.

    [13] Gharibzahedi S M T, Ghahderijani M, Lajevardi Z S. Specific heat, thermal conductivity and thermal diffusivity of red lentil seed as a function of moisture content[J]. Journal of Food Processing & Preservation, 2014, 38(4): 1807-1811.

    [14] Kara M, Boydas M G, Kara Y A, et al. The effect of moisture content on the thermal properties of red lentil seeds[J]. Transactions of the ASABE, 2012, 55(6): 2301-2306.

    [15] Subramanian S, Viswanathan R. Thermal properties of minor millet grains and flours[J]. Biosystems Engineering, 2003, 84(3): 289-296.

    [16] 楊洲,羅錫文,李長友. 稻谷熱特性參數(shù)的試驗測定[J]. 農(nóng)業(yè)機械學報,2003,34(4):76-78. Yang Zhou, Luo Xiwen, Li Changyou. Measurement of thermal properties of rough rice[J]. Transactions of the Chinese Society for Agricultural Machinery, 2003, 34(4): 76-78. (in Chinese with English abstract)

    [17] 周祖鍔,趙世宏,曹崇文. 谷物和種子的熱特性研究[J]. 北京農(nóng)業(yè)工程大學學報,1988,8(3):31-39. Zhou Zu’e, Zhao Shihong, Cao Chongwen. Thermal properties of grain[J]. Journal of Beijing Agricultural Engineering University, 1988, 8(3): 31-39. (in Chinese with English abstract)

    [18] Kaletun? G. Prediction of specific heat of cereal flours: A quantitative empirical correlation[J]. Journal of Food Engineering, 2007, 82(4): 589-594.

    [19] Shrestha B L, Baik O D. Thermal conductivity, specific heat, and thermal diffusivity of Saponaria vaccaria seed particles[J]. Transactions of the ASABE, 2010, 53(5): 1717-1725.

    [20] Yu D U, Shrestha B L, Baik O D. Thermal conductivity, specific heat, thermal diffusivity, and emissivity of stored canola seeds with their temperature and moisture content[J]. Journal of Food Engineering, 2015, 165: 156-165.

    [21] Jian F, Jayas D S, White N D G. Specific heat, thermal diffusivity, and bulk density of genetically modified canola with high oil content at different moisture contents, temperatures, and storage times[J]. Transactions of the ASABE, 2013, 56(3): 1077-1083.

    [22] Singh K K, Goswami T K. Thermal properties of cumin seed[J]. Journal of Food Engineering, 2000, 45(4): 181-187.

    [23] Shiva K M, Durga P B, Basavaraj M. Effect of moisture content and temperature on thermal conductivity ofby line heat source method (transient analysis)[J]. International Journal of Heat and Mass Transfer, 2014, 78: 354-359.

    [24] Perussello C A, Mariani V C, Masson M L, et al. Thermophysical properties of yacon (): Experimental determination and effect of moisture content[J]. Acta Scientiarum Technology, 2015, 37(1): 167-173.

    [25] Shrivastava M, Datta A K. Determination of specific heat and thermal conductivity of mushrooms ()[J]. Journal of Food Engineering, 1999, 39(3): 255-260.

    [26] Quirion S B, Villeneuve S, Leblanc D I, et al. Thermophysical properties and thermal behavior of leafy vegetables packaged in clamshells[J]. Journal of Food Engineering, 2012, 113(1): 27-32.

    [27] Obot M S, Li C, Fang T, et al. Measurement of thermal properties of white radish () using easily constructed probes[J]. PLoS ONE, 2017, 12(3): e0171016. doi:10.1371/journal.pone.0171016

    [28] 王紅英,李倪薇,高蕊,等. 不同前處理對飼料玉米比熱的影響[J]. 農(nóng)業(yè)工程學報,2012,28(14):269-276. Wang Hongying, Li Niwei, Gao Rui, et al. Effects of different pretreatments on specific heat of forage maize[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(14): 269-276. (in Chinese with English abstract)

    [29] 王紅英,高蕊,李軍國,等. 不同原料組分的配合飼料比熱模型[J]. 農(nóng)業(yè)工程學報,2013,29(9):285-292. Wang Hongying, Gao Rui, Li Junguo, et al. Model of formula feed specific heat based on different components of feed ingredients[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(9): 285-292. (in Chinese with English abstract)

    [30] 孔丹丹,陳嘯,楊潔,等. 仔豬配合料比熱預測模型的構(gòu)建[J]. 農(nóng)業(yè)工程學報,2016,32(18):307-314. Kong Dandan, Chen Xiao, Yang Jie, et al. Establishment of specific heat prediction model for weaned piglet mash feed[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(18): 307-314. (in Chinese with English abstract)

    [31] 孔丹丹,方鵬,王紅英,等. 高含量乳清粉的仔豬配合飼料熱特性及調(diào)質(zhì)溫度控制[J]. 農(nóng)業(yè)工程學報,2017,33(16):299-307. Kong Dandan, Fang Peng, Wang Hongying, et al. Thermal properties and conditioning temperature control of formula feeds containing high content of whey powder for weanling pigs[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(16): 299-307. (in Chinese with English abstract)

    [32] Hooper F C, Lepper F R. Transient heat flow apparatus for the determination of thermal conductivities[J]. Transactions of ASHVE, 1950, 56: 309-322.

    [33] Rafe A, Sadeghian A. Stabilization of Tarom and Domesiah cultivars rice bran: Physicochemical, functional and nutritional properties[J]. Journal of Cereal Science, 2017, 74: 64-71.

    [34] Park J W, Jang E Y, Ji Y K, et al. Effects of visible light irradiation on the oxidative stability in rice bran[J]. Journal of Cereal Science, 2013, 58(1): 178-181.

    [35] 曹文明,薛斌,袁超,等. 油脂氧化酸敗研究進展[J]. 糧食與油脂,2013,26(3):1-5. Cao Wenming, Xue Bin, Yuan Chao, et al. Research progress on the oxidative rancidity of oils and fats[J]. Cereals & Oils, 2013, 26(3): 1-5. (in Chinese with English abstract)

    [36] 徐浩,張海玲,顧廣東,等. 米糠加工及綜合利用研究進展[J]. 糧油食品科技,2017,25(5):37-41. Xu Hao, Zhang Hailing, Gu Guangdong, et al. Research progress on processing and comprehensive utilization of rice bran[J]. Science and Technology of Cereals, Oils and Foods, 2017, 25(5): 37-41. (in Chinese with English abstract)

    [37] Aviara N A, Haque M A. Moisture dependence of thermal properties of sheanut kernel[J]. Journal of Food Engineering, 2001, 47(2): 109-113.

    [38] Razavi S M A, Taghizadeh M. The specific heat of pistachio nuts as affected by moisture content, temperature, and variety[J]. Journal of Food Engineering, 2007, 79(1): 158-167.

    [39] Yang W, Sokhansanj S, Tang J, et al. Determination of thermal conductivity, specific heat and thermal diffusivity of borage seeds[J]. Biosystems Engineering, 2002, 82(2): 169-176.

    [40] Jian F, Jayas D S, White N D G. Thermal conductivity, bulk density, and germination of a canola variety with high oil content under different temperatures, moisture contents, and storage periods[J]. Transactions of the ASABE, 2012, 55(5): 1837-1843.

    Effect of temperature and particle size on thermophysical properties of different energy feedstuffs material

    Kong Dandan, Fang Peng, Jin Nan, Duan Enze, Chen Jiyuan, Wang Hongying※

    (,,,100083,)

    Energy feedstuffs, the most important and indispensable components in the pelleted complete feeds, provide energy in the form of soluble carbohydrate, oil and fat. Energy feedstuffs mainly includes cereals, cereal brans, molasses, oil, tuber and tuberous root, among which corn, wheat, barley, sorghum, wheat/rice bran, pelleted cassava residue and beet pulp are widely utilized in diets for livestock and poultry. Knowledge of the thermal properties of different feedstuffs is required for designing and parameter optimization of thermal process such as conditioning and pelleting processing, especially when energy costs, pellet quality and uniformity are the main consideration. The specific heat, thermal conductivity and thermal diffusivity are the most important thermal properties. In this paper, eight energy feed ingredients (4 cereal grains: corn, wheat, barley, sorghum; 4 processed by-products:wheat bran, pelletedcassava residue, pelleted beet pulp, rice bran) were obtained from a feed mill in the suburb of Beijing, and ground in a mini-mill equipped with screens having a aperture size of 1.5, 2.0 or 2.5 mm to obtain feed meals with 3 particle sizes (except rice bran, which could pass through the 1.5 mm screen completely without being ground). The moisture content of all feed meal samples was conditioned to 12% wet basis in order to eliminate the effect of moisture content difference on thermal properties. The specific heat, thermal conductivity of feed meals at the temperature range of 25-100 ℃were determined by DSC (differential scanning calorimetry) and KD2 Pro (thermal properties analyzer), respectively. Thermal diffusivity of feed meals was calculated using measured thermal conductivity, specific heat, and bulk density. Based on the experimental data, prediction models of three thermal properties as a function of temperature were established. The effects of particle size on thermal properties of feedstuffs were analyzed. The particle size of 8 feedstuffs grinding through screens with aperture size of 1.5, 2.0 and 2.5 mm ranged from 257.61 to 511.79m, 325.65 to 594.83m and 335.28 to 671.05m, respectively, which increased significantly with the increase in mesh size. There were significant differences between the values of particle size of 8 feedstuffs through the same mesh size. The values of specific heat of 4 grain meals (passing a screen with aperture size of 1.5 mm) increased with temperature (25-100 ℃) and ranged from 1.678 to 2.421, 1.644 to 2.402, 1.642 to 2.385, and 1.649 to 2.441 kJ/(kg·K), respectively. For wheat bran, cassava residue and beet pulp, the values were found to range from 1.621 to 2.540, 1.580 to 2.195, 1.746 to 2.351 kJ/(kg·K). The specific heat of rice bran raised from 1.966 kJ/(kg·K) at 25 ℃ to 1.985 kJ/(kg·K) at peak temperature (28.0 ℃), then declined to 1.868 kJ/(kg·K) at end-set temperature (33.5 ℃), and finally increased to 2.671 kJ/(kg·K) at 100 ℃. There was no significant effect of particle size on the specific heat of all feedstuffs (>0.05). The specific heat followed linear relationships with temperature at the range of 25-100 ℃ for 4 grains, and cubic polynomial relationships forwheat bran, cassava residue and rice bran, and a second order polynomial relationship for beet pulp. The values of thermal conductivity of 4 grain meals (passing a 1.5 mm screen) increased nonlinearly with temperature (25-100 ℃) and varied from 0.089 to 0.299, 0.078 to 0.336, 0.088 to 0.276 and 0.082 to 0.288 W/(m·K), respectively. For 4 processed by-products, the thermal conductivity was found to raise from 0.054 to 0.190, 0.092 to 0.362, 0.085 to 0.263, and 0.064 to 0.202 W/(m·K), respectively. Significant differences were observed between the thermal conductivity of 8 feed mills (passing a 1.5 mm screen) at a higher temperature (85 or 100 ℃), and the thermal conductivity of cassava residue was the highest, followed by wheat, wheat bran with the lowest value. For the same feedstuffs, the thermal conductivity of feed mills, at the same temperature, decreased generally with the increase in particle size, and the higher temperature, the stronger decrease. Cubic polynomial relationships of thermal conductivity with temperature (25-100 ℃) were established for 4 grains and beet pulp with different particle sizes, and quadratic polynomial relationships were obtained for wheat bran, cassava residue and rice bran. The values of thermal diffusivity of 4 grain meals (passing a 1.5 mm screen) varied from 7.459×10-8to 17.320×10-8, 6.694×10-8to 19.195×10-8, 7.878×10-8to 17.045×10-8, and 6.857×10-8to 16.109 ×10-8m2/s at temperature range of 25-100 ℃, respectively. The values was observed to ascend from 9.543×10-8to 20.892×10-8, 8.217×10-8to 23.254×10-8, 7.258×10-8to 16.646×10-8and 8.092×10-8to 18.718×10-8m2/s for 4 processed by-products, respectively. In general, the thermal diffusivity of the 8 energy feedstuffs had a slow increase with temperature at the range of 25-55 ℃, and a sharply raise with temperature at the range of 55-100 ℃. The thermal diffusivity of all 8 feedstuffs tended to decline with increasing particle size. The thermal diffusivity displayed second order polynomial relationship with temperature at the range of 25-100 ℃ for wheat bran and rice bran, and cubic polynomial relationships for the other 6 feedstuffs.

    feed; specific heat; thermal conductivity; thermal diffusivity; temperature; particle size

    10.11975/j.issn.1002-6819.2019.06.036

    S816.8

    A

    1002-6819(2019)-06-0296-11

    2018-10-24

    2019-03-13

    國家重點研發(fā)計劃項目“畜禽養(yǎng)殖綠色安全飼料飼養(yǎng)新技術(shù)研發(fā)(2018YFD0500600)”

    孔丹丹,博士生,主要從事飼料加工工藝研究。Email:dandank@cau.edu.cn

    王紅英,教授,博士生導師,主要從事飼料加工工藝技術(shù)與設(shè)備及畜禽養(yǎng)殖技術(shù)與裝備研究。Email:hongyingw@cau.edu.cn

    中國農(nóng)業(yè)工程學會會員:王紅英(E041200500S)

    孔丹丹,方 鵬,金 楠,段恩澤,陳計遠,王紅英. 溫度和粉碎粒度對不同能量飼料原料熱物理特性的影響[J]. 農(nóng)業(yè)工程學報,2019,35(6):296-306. doi:10.11975/j.issn.1002-6819.2019.06.036 http://www.tcsae.org

    Kong Dandan, Fang Peng, Jin Nan, Duan Enze, Chen Jiyuan, Wang Hongying. Effect of temperature and particle size on thermophysical properties of different energy feedstuffs material[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(6): 296-306. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.06.036 http://www.tcsae.org

    猜你喜歡
    麥麩米糠粉料
    近凈成型粉體制備技術(shù)的研究
    廣東建材(2023年10期)2023-10-17 08:25:42
    麥麩喂豬你做對了嗎
    麥麩喂豬,你做對了嗎?
    預拌混凝土粉料入庫智能防錯系統(tǒng)研究
    四川建筑(2019年6期)2019-07-20 09:39:32
    麥麩喂豬,你做對了嗎?
    干法制粉與噴霧干燥制粉的粉料性能對比研究
    佛山陶瓷(2017年7期)2017-09-06 21:27:34
    用針線“補”
    一種粉料氣力輸送槽
    化學法米糠穩(wěn)定化技術(shù)研究
    模糊數(shù)學方法應(yīng)用于麥麩香茶研制
    食品科學(2013年17期)2013-03-11 18:26:48
    一级爰片在线观看| 热99久久久久精品小说推荐| 亚洲精品一二三| 新久久久久国产一级毛片| 另类精品久久| 一二三四在线观看免费中文在| 日本欧美国产在线视频| 91精品伊人久久大香线蕉| 成人国产av品久久久| 满18在线观看网站| 久久99一区二区三区| 又大又爽又粗| 中文字幕人妻熟女乱码| 女人精品久久久久毛片| 自拍欧美九色日韩亚洲蝌蚪91| 日韩av免费高清视频| 中国国产av一级| 亚洲婷婷狠狠爱综合网| a 毛片基地| 欧美成人精品欧美一级黄| 宅男免费午夜| 亚洲精品久久久久久婷婷小说| 一级毛片 在线播放| 久久精品亚洲熟妇少妇任你| 久久久国产欧美日韩av| 免费黄色在线免费观看| 亚洲精品国产av成人精品| 婷婷色综合www| 国产精品免费视频内射| 黄色 视频免费看| 亚洲少妇的诱惑av| 成人午夜精彩视频在线观看| 黄色一级大片看看| 黄频高清免费视频| 这个男人来自地球电影免费观看 | 亚洲国产最新在线播放| a 毛片基地| 亚洲专区中文字幕在线 | 亚洲综合色网址| 黑人猛操日本美女一级片| 久久精品亚洲熟妇少妇任你| avwww免费| 天堂中文最新版在线下载| 十分钟在线观看高清视频www| 国产1区2区3区精品| 国产精品久久久久久精品古装| 亚洲精华国产精华液的使用体验| 十八禁人妻一区二区| 多毛熟女@视频| 久久久久精品性色| 男人添女人高潮全过程视频| 母亲3免费完整高清在线观看| av网站在线播放免费| 老司机影院成人| 国产免费又黄又爽又色| 男女免费视频国产| 色综合欧美亚洲国产小说| av不卡在线播放| 精品亚洲成a人片在线观看| 国产av精品麻豆| 性少妇av在线| 黄色视频在线播放观看不卡| 中文字幕另类日韩欧美亚洲嫩草| 极品少妇高潮喷水抽搐| 国产一区二区三区av在线| 最新在线观看一区二区三区 | 精品国产一区二区三区四区第35| 九草在线视频观看| 国产伦人伦偷精品视频| 国产欧美日韩综合在线一区二区| 一级片'在线观看视频| 天堂8中文在线网| 国产成人精品久久久久久| 成年动漫av网址| 国产探花极品一区二区| 免费观看av网站的网址| 免费人妻精品一区二区三区视频| 各种免费的搞黄视频| 青青草视频在线视频观看| 国产精品熟女久久久久浪| 三上悠亚av全集在线观看| 日韩精品免费视频一区二区三区| 女人久久www免费人成看片| 亚洲国产欧美一区二区综合| 性少妇av在线| av女优亚洲男人天堂| 精品卡一卡二卡四卡免费| 一级毛片黄色毛片免费观看视频| 香蕉国产在线看| 免费av中文字幕在线| 一边摸一边抽搐一进一出视频| 高清在线视频一区二区三区| 在线精品无人区一区二区三| 国产福利在线免费观看视频| 日韩不卡一区二区三区视频在线| 久久久久人妻精品一区果冻| 韩国av在线不卡| 国产av一区二区精品久久| 国产成人精品久久久久久| 色吧在线观看| 91精品伊人久久大香线蕉| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 美女大奶头黄色视频| 色婷婷av一区二区三区视频| 精品第一国产精品| 婷婷色综合www| 亚洲成人免费av在线播放| 午夜福利免费观看在线| 亚洲一卡2卡3卡4卡5卡精品中文| 青春草亚洲视频在线观看| 好男人视频免费观看在线| 国产精品一二三区在线看| 国产精品香港三级国产av潘金莲 | 欧美黑人精品巨大| 狠狠婷婷综合久久久久久88av| 精品国产乱码久久久久久小说| 欧美中文综合在线视频| 丰满少妇做爰视频| 国产精品99久久99久久久不卡 | 日韩一本色道免费dvd| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲最大av| 亚洲自偷自拍图片 自拍| 女人被躁到高潮嗷嗷叫费观| 考比视频在线观看| 精品国产露脸久久av麻豆| 免费观看a级毛片全部| 9色porny在线观看| 欧美精品亚洲一区二区| 丁香六月天网| 看免费成人av毛片| 国产极品天堂在线| 精品一区二区三区av网在线观看 | 精品国产乱码久久久久久男人| 香蕉国产在线看| 亚洲,欧美,日韩| 啦啦啦 在线观看视频| 国产黄色免费在线视频| 亚洲精品国产色婷婷电影| 免费黄色在线免费观看| 欧美激情 高清一区二区三区| 亚洲精品久久久久久婷婷小说| 精品第一国产精品| 韩国av在线不卡| 欧美在线黄色| 香蕉丝袜av| 美女主播在线视频| 丁香六月欧美| 97在线人人人人妻| 欧美xxⅹ黑人| 自拍欧美九色日韩亚洲蝌蚪91| 精品久久久精品久久久| 999久久久国产精品视频| 亚洲成色77777| 国产精品秋霞免费鲁丝片| 男男h啪啪无遮挡| 美女视频免费永久观看网站| 国产av一区二区精品久久| 最近中文字幕2019免费版| 在线观看免费午夜福利视频| 日韩av免费高清视频| 免费日韩欧美在线观看| 欧美97在线视频| 亚洲一码二码三码区别大吗| 99久久精品国产亚洲精品| 人人澡人人妻人| 不卡av一区二区三区| 欧美成人精品欧美一级黄| 黄片小视频在线播放| 日韩av不卡免费在线播放| 中文精品一卡2卡3卡4更新| 国产免费一区二区三区四区乱码| 亚洲欧美激情在线| 国产精品久久久人人做人人爽| 制服丝袜香蕉在线| 少妇人妻久久综合中文| 午夜福利视频精品| www.自偷自拍.com| 九草在线视频观看| 国产熟女午夜一区二区三区| 午夜免费鲁丝| 超碰97精品在线观看| 国产av一区二区精品久久| 国产福利在线免费观看视频| 久久精品久久精品一区二区三区| 午夜福利视频在线观看免费| 亚洲av国产av综合av卡| 最近手机中文字幕大全| 国产伦理片在线播放av一区| av又黄又爽大尺度在线免费看| 亚洲三区欧美一区| 成人毛片60女人毛片免费| 综合色丁香网| 国产成人系列免费观看| 久久性视频一级片| 老司机在亚洲福利影院| 成人免费观看视频高清| 中国国产av一级| 高清在线视频一区二区三区| 一区在线观看完整版| 99久国产av精品国产电影| 亚洲成av片中文字幕在线观看| 国产成人免费观看mmmm| 亚洲精品中文字幕在线视频| 纵有疾风起免费观看全集完整版| 王馨瑶露胸无遮挡在线观看| 捣出白浆h1v1| 搡老乐熟女国产| 伊人久久国产一区二区| 亚洲欧美一区二区三区国产| 精品人妻熟女毛片av久久网站| 操美女的视频在线观看| 老汉色∧v一级毛片| 制服丝袜香蕉在线| xxxhd国产人妻xxx| kizo精华| 亚洲成国产人片在线观看| 一区二区三区激情视频| 亚洲国产精品一区二区三区在线| 亚洲国产欧美网| av福利片在线| 久久久久久久精品精品| 日本午夜av视频| 精品一区二区免费观看| 午夜日韩欧美国产| 老司机靠b影院| 一级毛片我不卡| 久久久久久久大尺度免费视频| 考比视频在线观看| 精品国产露脸久久av麻豆| 国产精品秋霞免费鲁丝片| 老司机影院成人| 无遮挡黄片免费观看| 精品免费久久久久久久清纯 | 精品酒店卫生间| 看非洲黑人一级黄片| 国产成人a∨麻豆精品| 人人妻,人人澡人人爽秒播 | 欧美人与善性xxx| 欧美人与性动交α欧美精品济南到| 校园人妻丝袜中文字幕| 天天躁夜夜躁狠狠躁躁| 久久人妻熟女aⅴ| √禁漫天堂资源中文www| 少妇 在线观看| 免费高清在线观看日韩| 免费观看性生交大片5| 日韩电影二区| 性少妇av在线| 自线自在国产av| 国产精品国产三级专区第一集| 我要看黄色一级片免费的| 啦啦啦 在线观看视频| 男女边吃奶边做爰视频| 精品人妻一区二区三区麻豆| 国产精品一国产av| 免费黄色在线免费观看| 亚洲av欧美aⅴ国产| 日韩欧美精品免费久久| 90打野战视频偷拍视频| 亚洲色图综合在线观看| 免费看不卡的av| 国产一区亚洲一区在线观看| av在线app专区| 99香蕉大伊视频| 精品少妇一区二区三区视频日本电影 | 日韩电影二区| 大码成人一级视频| 波多野结衣av一区二区av| av免费观看日本| 亚洲欧美一区二区三区国产| 免费观看av网站的网址| 日韩欧美精品免费久久| 午夜福利在线免费观看网站| 在线观看免费高清a一片| 亚洲,欧美精品.| 啦啦啦在线观看免费高清www| 亚洲av成人不卡在线观看播放网 | 国产成人精品无人区| 丝袜美足系列| 亚洲中文av在线| 亚洲av电影在线观看一区二区三区| 日本欧美视频一区| 亚洲欧洲国产日韩| 亚洲av综合色区一区| 日韩一本色道免费dvd| 欧美人与性动交α欧美软件| 欧美日韩综合久久久久久| 一区二区三区四区激情视频| 亚洲国产精品999| 国产精品免费大片| 青春草视频在线免费观看| 亚洲伊人久久精品综合| 99久久人妻综合| 欧美人与性动交α欧美精品济南到| 侵犯人妻中文字幕一二三四区| 欧美激情高清一区二区三区 | 国产成人精品久久二区二区91 | 九草在线视频观看| 亚洲色图综合在线观看| 啦啦啦 在线观看视频| 在线观看国产h片| 秋霞伦理黄片| 多毛熟女@视频| 巨乳人妻的诱惑在线观看| 精品一区二区免费观看| 高清在线视频一区二区三区| 少妇人妻久久综合中文| 亚洲精品国产av蜜桃| 亚洲美女黄色视频免费看| 久久久久国产精品人妻一区二区| 久久久亚洲精品成人影院| 在现免费观看毛片| 免费观看人在逋| 9191精品国产免费久久| 人人妻人人澡人人爽人人夜夜| 午夜免费鲁丝| 高清视频免费观看一区二区| 人人澡人人妻人| 国产不卡av网站在线观看| 日本猛色少妇xxxxx猛交久久| 丝袜喷水一区| 国产精品亚洲av一区麻豆 | 欧美中文综合在线视频| 久久精品aⅴ一区二区三区四区| av又黄又爽大尺度在线免费看| 极品少妇高潮喷水抽搐| 七月丁香在线播放| 男女午夜视频在线观看| 最新在线观看一区二区三区 | 久久精品aⅴ一区二区三区四区| 中文字幕色久视频| 国产精品久久久av美女十八| 在线观看三级黄色| 久久婷婷青草| 色综合欧美亚洲国产小说| 免费黄频网站在线观看国产| 中文字幕色久视频| 免费黄频网站在线观看国产| 国产精品人妻久久久影院| 日本91视频免费播放| 少妇人妻 视频| 亚洲精品国产av蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美成人精品欧美一级黄| 国产精品一国产av| 在线亚洲精品国产二区图片欧美| 观看av在线不卡| 亚洲国产欧美日韩在线播放| 久久精品久久久久久噜噜老黄| 晚上一个人看的免费电影| 如日韩欧美国产精品一区二区三区| 欧美成人午夜精品| 国产精品无大码| 妹子高潮喷水视频| 国产精品一国产av| 纯流量卡能插随身wifi吗| 制服诱惑二区| 狠狠精品人妻久久久久久综合| 亚洲一卡2卡3卡4卡5卡精品中文| 男人舔女人的私密视频| 人人妻人人澡人人看| 三上悠亚av全集在线观看| 中文字幕人妻熟女乱码| 伊人久久大香线蕉亚洲五| 99精品久久久久人妻精品| 黄色 视频免费看| 日韩熟女老妇一区二区性免费视频| 人人妻人人添人人爽欧美一区卜| 19禁男女啪啪无遮挡网站| 国产av国产精品国产| 在线天堂中文资源库| 亚洲精品久久午夜乱码| 五月开心婷婷网| 亚洲精品av麻豆狂野| 国产一区二区三区综合在线观看| 午夜福利视频精品| 丝袜美腿诱惑在线| 国产精品熟女久久久久浪| 亚洲七黄色美女视频| 老司机亚洲免费影院| 美女脱内裤让男人舔精品视频| 一级毛片我不卡| 国产有黄有色有爽视频| 成年人午夜在线观看视频| 亚洲av欧美aⅴ国产| 九色亚洲精品在线播放| 亚洲人成电影观看| 亚洲三区欧美一区| 99精品久久久久人妻精品| 国产欧美日韩一区二区三区在线| 极品人妻少妇av视频| 精品人妻一区二区三区麻豆| 777米奇影视久久| 亚洲国产看品久久| 国产视频首页在线观看| 国产精品二区激情视频| 啦啦啦视频在线资源免费观看| 国产一区二区激情短视频 | 黄片播放在线免费| 国产在线免费精品| 国产免费一区二区三区四区乱码| 天堂中文最新版在线下载| 欧美国产精品va在线观看不卡| 高清黄色对白视频在线免费看| 丝袜脚勾引网站| 日韩av在线免费看完整版不卡| 国产乱来视频区| 成人手机av| 女人被躁到高潮嗷嗷叫费观| 性少妇av在线| 国产淫语在线视频| 久久午夜综合久久蜜桃| 国产欧美亚洲国产| 美国免费a级毛片| 制服丝袜香蕉在线| 侵犯人妻中文字幕一二三四区| 只有这里有精品99| 男人舔女人的私密视频| 18禁裸乳无遮挡动漫免费视频| 欧美在线黄色| 制服人妻中文乱码| 欧美人与善性xxx| 黄色毛片三级朝国网站| 免费观看人在逋| 美女国产高潮福利片在线看| 一区二区三区精品91| 如何舔出高潮| 十分钟在线观看高清视频www| 国产精品久久久av美女十八| 激情五月婷婷亚洲| 亚洲欧美成人精品一区二区| 高清不卡的av网站| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久久精品电影小说| 1024视频免费在线观看| 欧美 日韩 精品 国产| 亚洲国产av新网站| 国产淫语在线视频| 青青草视频在线视频观看| 五月开心婷婷网| 婷婷色综合www| 国产精品久久久久久人妻精品电影 | 大香蕉久久网| 日韩免费高清中文字幕av| 欧美精品一区二区大全| 超碰97精品在线观看| 少妇人妻 视频| 丁香六月天网| www.熟女人妻精品国产| 九草在线视频观看| 亚洲四区av| 午夜日本视频在线| 亚洲国产欧美日韩在线播放| 亚洲免费av在线视频| 日韩成人av中文字幕在线观看| 欧美黄色片欧美黄色片| 少妇猛男粗大的猛烈进出视频| 免费观看人在逋| 一级a爱视频在线免费观看| 国产精品一国产av| 丁香六月天网| 欧美精品av麻豆av| 久久久国产精品麻豆| 热99国产精品久久久久久7| 亚洲av国产av综合av卡| 三上悠亚av全集在线观看| 18禁裸乳无遮挡动漫免费视频| 国产av一区二区精品久久| xxxhd国产人妻xxx| av.在线天堂| 国产精品一区二区精品视频观看| 国产熟女午夜一区二区三区| 亚洲国产毛片av蜜桃av| 成人18禁高潮啪啪吃奶动态图| 黄色视频不卡| 99久国产av精品国产电影| 国产极品天堂在线| av线在线观看网站| www.精华液| 不卡视频在线观看欧美| 免费在线观看视频国产中文字幕亚洲 | 成年动漫av网址| 91老司机精品| 老司机影院毛片| 久久久久精品性色| 观看av在线不卡| 少妇精品久久久久久久| 亚洲国产精品一区三区| 大陆偷拍与自拍| 午夜91福利影院| 黄频高清免费视频| 超碰成人久久| 成年美女黄网站色视频大全免费| 亚洲人成电影观看| a 毛片基地| 亚洲av日韩在线播放| 日韩 亚洲 欧美在线| 中文天堂在线官网| 青春草国产在线视频| 日韩一区二区视频免费看| 搡老乐熟女国产| 久久av网站| 人人妻人人添人人爽欧美一区卜| 国产在线视频一区二区| 久久久久久久精品精品| 黄色一级大片看看| 精品久久久精品久久久| 日韩 欧美 亚洲 中文字幕| 亚洲综合色网址| 午夜福利一区二区在线看| 精品一区在线观看国产| 国产伦理片在线播放av一区| 最近中文字幕高清免费大全6| 久久精品人人爽人人爽视色| 天天躁日日躁夜夜躁夜夜| 亚洲第一区二区三区不卡| a级毛片黄视频| 男女免费视频国产| 侵犯人妻中文字幕一二三四区| 日本av免费视频播放| 亚洲美女搞黄在线观看| 日本wwww免费看| 久久人人97超碰香蕉20202| 国产97色在线日韩免费| 国产又爽黄色视频| 欧美少妇被猛烈插入视频| 高清视频免费观看一区二区| 日日爽夜夜爽网站| 精品国产乱码久久久久久男人| 日韩欧美精品免费久久| 男女下面插进去视频免费观看| 亚洲欧美日韩另类电影网站| 亚洲美女黄色视频免费看| 自线自在国产av| 成人国产av品久久久| 日本爱情动作片www.在线观看| 在线天堂中文资源库| 大片免费播放器 马上看| 中文字幕人妻丝袜制服| 下体分泌物呈黄色| a级毛片在线看网站| 国产片特级美女逼逼视频| 久久久亚洲精品成人影院| 一二三四在线观看免费中文在| 国产精品 国内视频| 久热爱精品视频在线9| 国产有黄有色有爽视频| 久久久久久久精品精品| 美女高潮到喷水免费观看| 9191精品国产免费久久| 亚洲七黄色美女视频| 精品国产国语对白av| 亚洲久久久国产精品| 国产成人啪精品午夜网站| 国产午夜精品一二区理论片| 亚洲情色 制服丝袜| 国产一区二区激情短视频 | 国产精品偷伦视频观看了| 美女视频免费永久观看网站| bbb黄色大片| 新久久久久国产一级毛片| av网站在线播放免费| h视频一区二区三区| 国产日韩一区二区三区精品不卡| 国产成人一区二区在线| 国产成人免费无遮挡视频| 免费高清在线观看日韩| av不卡在线播放| 老鸭窝网址在线观看| 男人舔女人的私密视频| 亚洲精品国产色婷婷电影| 丝袜美腿诱惑在线| 欧美亚洲日本最大视频资源| 免费久久久久久久精品成人欧美视频| 色94色欧美一区二区| 一边摸一边做爽爽视频免费| 免费不卡黄色视频| 中文字幕人妻丝袜一区二区 | 欧美精品一区二区免费开放| 国产精品麻豆人妻色哟哟久久| 久久精品久久精品一区二区三区| av网站免费在线观看视频| 九九爱精品视频在线观看| 激情五月婷婷亚洲| 中文字幕人妻丝袜制服| 纯流量卡能插随身wifi吗| 伦理电影免费视频| 黄频高清免费视频| 999久久久国产精品视频| 国产亚洲最大av| 中文字幕人妻丝袜制服| 纯流量卡能插随身wifi吗| 欧美精品av麻豆av| 精品免费久久久久久久清纯 | 亚洲av日韩精品久久久久久密 | www.自偷自拍.com| 久久ye,这里只有精品| 制服丝袜香蕉在线| 嫩草影视91久久| 久久综合国产亚洲精品| 国产成人精品无人区| 精品亚洲乱码少妇综合久久| 观看美女的网站| 天天躁夜夜躁狠狠久久av| videos熟女内射| av国产精品久久久久影院| 男女边摸边吃奶| 亚洲欧美色中文字幕在线| 国产成人a∨麻豆精品| 久久精品亚洲av国产电影网| 精品视频人人做人人爽| 精品久久久精品久久久| av福利片在线|