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

    儲料豎向壓力對糧倉中小麥糧堆濕熱傳遞的影響

    2020-06-04 01:20:54陳桂香劉超賽蔣敏敏陳家豪王海濤張宏偉
    農業(yè)工程學報 2020年8期
    關鍵詞:高溫區(qū)糧堆儲藏

    陳桂香,劉超賽,蔣敏敏,陳家豪,王海濤,張宏偉

    儲料豎向壓力對糧倉中小麥糧堆濕熱傳遞的影響

    陳桂香,劉超賽,蔣敏敏,陳家豪,王海濤,張宏偉

    (河南工業(yè)大學土木工程學院,鄭州 450001)

    糧倉中存在壓力場、溫度場和濕度場等多物理場,為了得出各物理因子共同影響下的糧堆內濕熱傳遞規(guī)律,該研究利用自行研制的糧堆多場耦合試驗裝置,針對倉內小麥糧堆單元體,研究在高溫邊界38.5 ℃、低溫邊界5.2 ℃,初始糧溫25.8 ℃,豎向壓力分別為50、100、150 kPa條件下小麥糧堆濕熱傳遞情況。試驗結果表明:豎向壓力增加,糧堆孔隙率減小,熱量通過糧食籽粒間傳導增加,傳遞速率加快,豎向壓力從50 kPa 增大至150 kPa,糧溫較入倉時下降約0.5~1.3 ℃,溫度梯度變化率達8.7%,不同壓力下糧堆高溫區(qū)面積隨儲藏時間呈冪函數減小。糧堆內濕空氣在邊界處累積至峰值時會有部分濕空氣向糧堆內遷移。糧堆中部與靠近低溫邊界溫差大于6.3 ℃時,糧堆內濕空氣擴散加快,糧堆中部平均相對濕度下降速率隨豎向壓力增加而加快。研究結果可為散裝糧堆多場耦合研究提供理論支持。

    溫度;相對濕度;小麥糧堆;豎向壓力;多場耦合

    0 引 言

    中國小麥年產量超1.3億t,儲備周期通常長達3 a,保障糧食儲藏安全關系國計民生和國家安全。糧堆內存在壓力場、溫度場和濕度場等多物理場,儲糧過程中,季節(jié)更替、太陽輻射等外界環(huán)境因素,以及糧食籽粒呼吸、微生物代謝等作用均會引起糧堆內溫濕度的分布發(fā)生改變[1]。糧堆內豎向壓力隨糧堆深度的增加而增大,糧堆壓力場的變化直接影響糧堆內孔隙率的分布,進而影響熱量和濕分的遷移,部分區(qū)域溫度和濕分適宜時,極易引發(fā)糧食霉變與蟲害,危及儲糧安全[2]。因此,研究儲料豎向壓力對糧堆濕熱傳遞的影響對于安全儲糧極為重要。

    針對糧堆壓力場、溫度場和濕度場已有大量研究。在壓力場方面,陳桂香等[3-4]、王錄民等[5]、陳家豪等[6-7]分別對散裝糧堆空間壓力場進行模擬和實測研究,提出筒倉儲料靜態(tài)壓力場的計算方法,得出糧堆高度及倉壁摩擦力對倉底壓力的影響規(guī)律。Mounfield等[8]通過離散元法,得出豎直方向上貯料結拱時的應力分布。研究者還提出糧堆深度與孔隙率的關系,唐福元等[9]、張小正[10]分別建立了筒倉中稻谷、小麥孔隙率與糧堆深度的關系,糧堆深度40 m時,稻谷和小麥孔隙率分別為0.56和0.33,分別較頂層糧堆減小7.9%和15.6%。王娟[11]建立小麥糧堆孔隙率與豎向壓力的關系,當豎向壓力為150 kPa時,糧堆孔隙率為0.39,其樣品較小且加載邊界不能較好模擬糧堆內籽粒與籽粒之間的柔性擠壓。糧堆孔隙率隨壓力和糧堆深度變化較大,在分析大型倉儲結構儲糧安全時,應考慮堆高對濕熱傳遞的影響。

    關于糧堆內溫濕度場的分布狀態(tài)已有相當多的研究成果。針對糧堆內的濕熱傳遞問題,Jia等[12]建立和小麥糧堆密度相關的二維非線性傳熱有限元模型,模擬小型儲糧倉在內置熱源情況下小麥糧堆7 h內溫度分布,得出離倉中心距離增加,糧堆溫度變化減小。Thorpe[13]通過有限元法,建立太陽輻射對糧堆儲藏過程中熱量傳遞影響的數學模型。陳桂香等[14]、張忠杰等[15-16]分別模擬了平房倉和淺圓倉內準靜態(tài)狀態(tài)時溫度場的變化趨勢。Jian等[17]對散裝小麥糧堆的溫濕度分布進行長期監(jiān)測,分析不同季節(jié)糧堆內濕熱傳遞規(guī)律,得出距倉底1.6 m高處溫度梯度可達32.4℃/m,糧堆內部水分變化小于表層糧堆水分變化。

    在糧堆多場耦合研究方面,Hammami等[18]基于熱濕平衡原理,建立圓形筒倉糧堆內溫濕度變化模型。Gaston等[19]通過有限元法,得出季節(jié)變化時倉儲小麥的溫度分布和水分遷移規(guī)律。王遠成等[20-21]通過理論分析與數值模擬研究靜態(tài)儲藏時的濕熱耦合傳遞過程。王振華[22]取糧堆中單元體為研究對象,考慮靜態(tài)儲藏時溫度和含水率變化及氣體自然對流情況,建立與糧堆孔隙率有關的靜態(tài)儲藏模型。尹君[23]基于多場耦合理論提出小麥糧堆受季節(jié)氣溫變化影響的糧堆內溫度場、濕度場和微氣流場多場耦合規(guī)律。王小萌等[24-25]研究糧堆在內置熱源和糧堆霉變引起自發(fā)熱時糧堆內部溫濕度場的變化規(guī)律。

    目前對于糧堆多場耦合方面的研究大多集中于糧堆壓力場、溫度場的單獨研究和糧堆內熱濕耦合的研究,對糧堆內熱濕耦合的研究主要集中于數值模擬,多數研究缺少相關試驗的驗證。糧堆壓力通過改變孔隙率和導熱系數影響糧堆溫、濕度場分布規(guī)律的研究較少,因此難以準確反映糧堆內多物理因子共同作用的儲糧環(huán)境。

    本文針對糧堆壓力通過影響孔隙率和導熱系數進而影響糧堆內的濕熱傳遞問題,以散裝小麥糧堆單元體為研究對象,自行設計了多場耦合試驗平臺,考慮糧倉內不同深度豎向壓力對糧堆濕熱傳遞的影響,較好模擬糧堆溫濕度邊界條件。通過小麥糧堆濕熱傳遞試驗,模擬倉內糧堆單元的壓力和溫濕度邊界條件,研究豎向壓力改變對糧堆孔隙率和導熱系數等的影響,得出豎向壓力改變時,糧堆內溫度梯度和糧堆內相對濕度場的變化規(guī)律及糧溫升高區(qū)面積與儲藏時間之間的數學關系,以期為倉內糧堆多場耦合理論和儲糧安全研究提供理論支持。

    1 材料與方法

    1.1 試驗材料

    本試驗以河南產小麥為研究對象,初始含水率8.56%,為避免初始含水率不均對試驗造成的影響,模擬低水分夏糧入倉工況,試驗前對小麥樣品調質,并在25.8 ℃恒溫室內靜置10 d[24],使樣品溫濕度均一,調質后的小麥含水率為10.89%。

    1.2 試驗儀器

    101型電熱鼓風干燥箱:北京永光明醫(yī)療儀器有限公司;2638A-60-DC Fluke全能型數據采集器:美國FLUKE公司;散裝糧堆多場耦合實驗裝置:河南工業(yè)大學;伺服柔性加載控制系統(tǒng):建湖惠通機械有限公司;HX-08控溫系統(tǒng):江蘇天翎儀器有限公司;PT100高精度鉑電阻探頭:OMEGA公司;HT75濕度探頭:SENSIRION公司。

    1.3 試驗方法

    前期研究表明,筒倉內豎向壓力自糧堆頂面隨糧層深度增加呈線性增大,糧堆孔隙率隨之改變,進而引起糧堆內熱量和濕分的復雜傳遞,不同糧層深度溫濕度分布存在較大差別[7,23]。為了研究大型倉儲結構糧堆內豎向壓力對濕熱傳遞影響的多物理因子相互作用問題,本文以倉內某一糧堆單元為研究對象,控制糧堆單元體壓力、溫度和濕度邊界條件,研究豎向壓力對糧堆內濕熱傳遞的影響。

    糧堆單元體左右兩側施加均勻的溫度邊界,可模擬糧堆0~100 ℃溫度范圍,溫度波動±0.05 ℃。非溫控邊界鋪設剛度較大的納米板保溫墊,保證糧堆單元體的保溫性能和熱量在垂直于軸的平面單向傳遞,消除非溫控邊界對傳遞的影響。單元體頂部為與倉儲環(huán)境相近似的柔性橡膠氣囊,氣囊四周側壁涂薄層凡士林,消除氣囊與側壁之間的摩擦效應,能夠對單元均勻施加壓力[26],較好地模擬實倉糧堆力學性能。調整糧堆內溫濕度測點呈網格狀布置,如圖1所示,溫、濕度傳感器分3層布置,沿高度方向分別為=0.1、0.3、0.5 m,沿糧堆溫度梯度方向中垂面上測點分別距高溫邊界0.05、0.30、0.55 m,中心測點(b4)位于糧堆單元體的幾何中心,便于較好的研究濕分、熱量的單向傳遞[27],每個測點均布置溫度傳感器和濕度傳感器。傳感器精度:溫度傳感器±0.15 ℃,濕度傳感器2%。

    注:“”代表糧堆中垂面測點。各測點位置分別布置溫度傳感器和濕度傳感器。A、B、C為濕度傳感器,a、b、c為溫度傳感器。

    式中~x為糧堆內不同位置溫度值(或相對濕度值),℃(%);P()為溫度(或相對濕度)隨測點位置變化的函數,℃(%);()為溫度(或相對濕度)誤差隨測點位置變化的函數,℃(%);將所得到的云圖進行灰度和二值化處理,計算糧堆內高溫區(qū)面積。

    1.4 試驗方案

    試驗中將初始含水率為10.89%的小麥樣品均勻裝入糧堆單元內,并按照圖1布置溫濕度傳感器,糧堆單元總高度0.6 m,樣品總質量143.93 kg。初始糧溫25.8 ℃,空氣相對濕度61.5%。

    根據大型糧倉內糧堆的壓力范圍,選取試驗豎向壓力分別為50、100、150 kPa共3個等級,采用伺服柔性加載控制系統(tǒng)以25 kPa逐級加壓,每組試驗初始條件相同。由于糧堆導熱系數小,較小溫度差引起的濕熱遷移緩慢[25],為便于觀察豎向壓力通過改變孔隙率進而影響糧堆濕熱傳遞的規(guī)律,根據Jian等[17]和尹君[23]關于糧堆內溫度梯度和不同糧層深度濕熱傳遞規(guī)律研究成果,取高溫邊界溫度值38.5 ℃,低溫邊界溫度值5.2 ℃,目的是利用溫差引發(fā)糧堆內濕熱傳遞[24]。糧堆單元體置于恒溫恒濕環(huán)境中,設定室內溫度25.8 ℃、空氣相對濕度61.5%,模擬華北地區(qū)秋季氣候環(huán)境??紤]到糧食籽粒排列具有隨機性以及試驗數據采集的離散性,本研究相同條件重復3組試驗,取試驗平均值進行研究。

    2 結果與分析

    2.1 溫度變化

    2.1.1 糧堆單元的溫度變化

    不同豎向壓力下沿糧堆溫度梯度方向中垂面上9個測點的溫度變化曲線如圖2所示。

    注:a2~a6, b2~b6, c2~c6含義見圖1。

    由圖2可知,相同壓力下,儲藏0~5 h,因側壁溫度邊界與糧堆間持續(xù)的熱量交換,在控溫邊界與糧堆間的較大溫差作用下,靠近溫度邊界糧食糧溫基本呈直線趨勢變化,儲藏5~43.5 h,靠近控溫邊界測點溫度變化斜率逐漸減小,表明由于熱量不斷傳遞,控溫邊界與靠近控溫邊界測點間溫差減小,糧堆溫度變化減慢。靜態(tài)儲藏72 h,豎向壓力為50 kPa時,中間截面(平行于,=300 mm,以下簡稱“中面”)測點a4、b4、c4平均溫度變化速率最小,為0.007 ℃/h,100 kPa時中面測點平均溫度變化速率為0.010 ℃/h,150 kPa時平均溫度變化速率最大,為0.018 ℃/h。表明豎向壓力增大,糧堆孔隙率減小,通過糧食籽粒間熱傳導增加,熱量傳遞速率加快。

    2.1.2 豎向壓力對熱量傳遞影響

    圖3為不同豎向壓力下糧堆中面測點(a4、b4、c4)的平均溫度變化。由圖可知,不同豎向壓力下,測點平均溫度變化趨勢基本一致,相同時刻溫度值存在差別。高溫邊界38.5 ℃、低溫邊界5.2 ℃,初始糧溫25.8 ℃條件下,隨著儲藏時間的變化,受糧堆溫度與高、低溫邊界間溫差的影響,中面糧堆平均糧溫在0~5 h時有一定升高趨勢,5 h后逐漸降低,豎向壓力為50 kPa(糧堆深度約6 m)時,平均糧溫最高,靜態(tài)儲藏72 h,糧溫較入倉時降低0.5 ℃;豎向壓力為150 kPa(糧堆深度約19 m)時糧溫最低,靜態(tài)儲藏72 h,糧溫較入倉時降低1.3 ℃。圖3結果表明,堆高較大時,糧堆孔隙率隨豎向壓力的增加而減小,糧食籽粒間接觸面積增加,通過糧食籽粒間的熱傳導增加,小麥堆導熱系數為0.129 6 W/(m·K),是空氣導熱系數的49倍[22],通過孔隙間微氣流對流傳熱減少,在低溫邊界與糧堆溫差作用下,熱量自糧堆內向近低溫邊界傳遞加快,中面平均糧溫更快趨向于平衡溫度。

    圖3 不同壓力下中面測點平均溫度變化

    溫差是影響糧堆內熱量傳遞的主要因素[25]。本文研究了不同豎向壓力下糧堆中面和靠近低溫邊界測點溫差隨時間的變化趨勢。因靠近低溫邊界測點離低溫邊界距離近,溫度下降較快,糧堆中面與近低溫邊界測點間溫差隨儲藏時間的變化逐漸增大。兩截面溫差變化趨勢如圖4所示,各壓力下溫差隨時間變化規(guī)律基本一致,相同時刻不同壓力糧堆中面與近低溫邊界測點間溫差存在差別:儲藏72 h時,50 kPa時溫度差最大,為10.3 ℃,100 kPa時,溫差為9.2 ℃,150 kPa時溫度差最小,為7.3 ℃,受豎向壓力影響,豎向壓力從50 kPa增大至150 kPa,溫差變化率達29%。結果表明,倉壁受氣候影響溫度降低時,糧堆熱量以熱傳導和自然對流的方式向糧堆邊界傳遞,糧堆深度增加,豎向壓力增大,糧堆孔隙率減小,糧食籽粒間接觸面積增加,熱傳導速率加快,自儲藏初始階段糧堆與低溫邊界間溫差更快趨向于平衡。

    圖4 不同壓力下中面和靠近低溫邊界測點間溫差

    2.1.3 高溫區(qū)面積變化

    長期儲糧過程中,由于糧食籽粒及微生物呼吸作用等,糧堆內易形成發(fā)熱點,導致糧堆局部溫度偏高,糧食長時間處于高溫區(qū)域環(huán)境易形成霉變,導致籽?;盍ο陆?,發(fā)芽率降低,甚至失去種用價值。糧堆單元體入倉糧溫為25.8 ℃,定義溫度高于25.8 ℃的區(qū)域為高溫區(qū)。圖5為不同壓力下中垂面糧堆溫度場云圖。由圖5可知,在控溫邊界與糧堆間較大的溫差作用下,熱量自控溫邊界向糧堆內部傳遞,儲藏0~5 h,在高溫邊界與糧堆間的溫差作用下,高溫區(qū)面積逐漸增大,儲藏5 h時高溫區(qū)面積達到最大,各壓力狀態(tài)下糧堆單元體溫度分布較為均勻,受豎向壓力與倉壁摩擦力影響,近高溫邊界糧堆孔隙率隨壓力增大而減小,通過糧食籽粒間熱傳導增加,孔隙間微氣流對流傳熱受阻,近控溫邊界溫度分布存在差異。因低溫邊界與糧堆間較大的溫差作用,糧堆內熱量不斷向低溫邊界傳遞,高溫區(qū)面積在靜態(tài)儲藏5 h后逐漸減小。靜態(tài)儲藏43.5 h時,由于糧堆單元試驗條件限制,難以絕對隔熱,高溫邊界上下兩端存在熱量散失,中心處較上下兩端溫度高。壓力從50 kPa增大至150 kPa時,糧堆單元體高溫區(qū)所占比例減小。

    在高溫邊界38.5 ℃、低溫邊界5.2 ℃,初始糧溫25.8 ℃條件下進行了糧堆高溫區(qū)面積比例與儲藏時間關系的研究。圖6為儲藏5 h后不同壓力下高溫區(qū)面積隨儲藏時間的變化圖,結果表明,各壓力下高溫區(qū)面積所占比例隨儲藏時間呈冪函數關系減小,表示為

    式中為高溫區(qū)面積所占比例,%;T為儲藏時間,h;a為反映糧堆高溫區(qū)面積隨儲藏時間降低的指數。對于指數a,豎向壓力為50 kPa時最大,為0.047(R2=0.97),隨著豎向壓力的增大,指數a的值越大,100 kPa時,指數a為0.052(R2=0.97),150 kPa時,為0.084(R2=0.95),R2均不小于0.95,表明所選數學模型在豎向壓力50~150 kPa范圍內適合小麥糧堆在該試驗條件下高溫區(qū)面積的擬合。

    圖6 不同壓力下高溫區(qū)面積隨儲藏時間變化

    由圖6可知,靜態(tài)儲藏72 h,豎向壓力為50 kPa時,高溫區(qū)面積所占比例從52.5%降至45.7%。豎向壓力增大到150 kPa時,高溫區(qū)面積所占比例從49.3%降至41.3%。豎向壓力增加,高溫區(qū)面積所占比例減小,由于壓力改變了糧堆孔隙率,進而影響了熱量傳遞速率,隨著儲藏時間的變化,高溫區(qū)面積變化受壓力影響顯著。

    2.1.4 溫度梯度變化

    溫度梯度是造成微氣流形成的主要原因[3],不僅可以使熱量發(fā)生傳遞,還會促使?jié)穹謹U散和再分配[28]。圖7為儲藏期間不同壓力下糧堆靠近高溫邊界測點與靠近低溫邊界測點之間溫度梯度變化曲線。儲藏72 h,豎向壓力為50 kPa時,糧堆單元溫度梯度最大,為34.6 ℃/m,100 kPa時溫度梯度為34.1 ℃/m,150 kPa時,糧堆溫度梯度最小,為31.6 ℃/m。恒定溫度場環(huán)境中,豎向壓力從50 kPa增大到150 kPa時,溫度梯度變化率達8.7%。大型倉儲結構中,糧層深度增加,糧堆豎向壓力增大,糧堆孔隙率減小,糧堆對流傳熱受阻,通過籽粒間熱傳導增加,糧堆豎向壓力影響糧堆內溫度場分布,近倉底區(qū)域糧堆內豎向壓力較大,熱量通過籽粒間熱傳導增加,糧堆溫度梯度減小。

    圖7 不同壓力下小麥糧堆溫度梯度

    2.2 相對濕度變化

    2.2.1 糧堆相對濕度變化

    溫差是引發(fā)糧堆內濕分遷移的重要因素,長期儲糧過程中,受外界環(huán)境和糧堆生態(tài)子系統(tǒng)的影響,在溫差作用下,糧堆內濕度重新分布[22]。圖8為不同壓力下中垂面糧堆相對濕度場云圖,由圖可知,小麥糧堆相對濕度分布在高、低溫邊界溫差作用下發(fā)生變化。高度≥0.3 m中上層糧堆,在糧堆與高溫邊界間較大的溫差作用下,小麥發(fā)生解吸,孔隙間水汽分壓升高,糧堆內濕空氣由高溫邊界向低溫邊界均勻遷移。高度<0.3 m中下層糧堆,受倉底不能絕對隔熱條件限制,近高溫邊界與近低溫邊界測點間溫差相比高度=0.3 m較小,所形成的微氣流作用減弱,近高溫邊界糧堆濕空氣隨微氣流向低溫邊界移動減少,近高溫邊界下端相對濕度較高。

    圖9為不同壓力下中面與近低溫邊界平均相對濕度變化。豎向壓力從50 kPa增大至150 kPa時,在本研究試驗條件下,儲藏72 h過程中,糧堆單元體中面相對濕度均逐漸降低,近低溫邊界均先升高后降低,并最終保持相對濕度不變。因兩控溫邊界與糧堆溫差的不同,初始糧溫與低溫邊界間溫差?1,大于糧堆與高溫邊界間溫差?2,所形成微氣流作用較強,因此糧堆內濕空氣隨微氣流由糧堆內部向低溫邊界遷移速率1大于由高溫邊界向糧堆內遷移速率2,圖1中A4、B4和C4測點的平均相對濕度不斷下降,形成中部低濕區(qū),如圖8所示。圖9b結果表明糧堆內濕空氣在靠近低溫邊界聚集,靠近低溫邊界平均相對濕度升高??拷蜏剡吔缂Z堆濕空氣累積至峰值,因濃度較高,部分濕空氣向糧堆內遷移,近低溫邊界相對濕度下降,并最終保持不變。

    圖8 不同壓力下中垂面糧堆相對濕度場云圖

    圖9 不同壓力下小麥糧堆平均相對濕度變化

    Fig 9 Changes of average relative humidity of wheat bulk under different pressures

    2.2.2 壓力對相對濕度分布的影響

    溫差是糧堆內濕熱遷移的動力源[25],壓力通過影響孔隙率,是影響糧堆濕熱遷移的重要因素。圖10為不同壓力下糧堆中面平均相對濕度隨靠近低溫邊界糧堆與中面溫度差的變化關系。由圖可知,隨著溫差的增大,中面糧堆平均相對濕度逐漸降低,因糧堆導熱系數小,糧堆內熱量傳遞較慢,溫差較小時,相對濕度隨溫差變化基本一致;溫差大于6.3 ℃時,所形成的微氣流作用增強,糧堆內濕空氣隨微氣流向近低溫邊界遷移,中面糧堆相對濕度迅速下降。壓力為150 kPa時,糧堆孔隙率較小,熱量通過籽粒間熱傳導增加,糧堆與低溫邊界間溫差更快趨向于平衡,平均相對濕度隨溫差變化最快。豎向壓力50 kPa時,糧堆孔隙率較大,糧堆內部熱量向低溫邊界遷移變慢,糧堆與低溫邊界間溫差相比較大,平均相對濕度隨溫差變化較慢。

    圖10 不同壓力下中面糧堆平均相對濕度隨溫差變化

    3 結 論

    本文針對大型倉儲結構儲糧環(huán)境中豎向壓力對濕熱傳遞影響問題,通過自行研制的多場耦合試驗平臺,以倉內小麥糧堆單元為研究對象,研究了高溫邊界38.5 ℃、低溫邊界5.2 ℃,初始糧溫25.8 ℃條件下,豎向壓力從50 kPa增大至150 kPa時糧堆內的熱量和濕分遷移規(guī)律,得出以下結論:

    1)豎向壓力增加,熱量通過籽粒間熱傳導作用增強,通過孔隙間微氣流對流傳熱減少,熱量傳遞速率加快,糧堆內溫度較入倉時下降約0.5~1.3 ℃,溫度梯度變化率達8.7%。50 kPa時糧堆中部與靠近低溫邊界間溫度差最大,為10.3 ℃,150 kPa時最小,為7.3 ℃。

    2)不同壓力下高溫區(qū)面積隨儲藏時間均呈冪函數關系減小,且冪函數的指數隨豎向壓力的增加而增大。

    3)糧堆內溫差大于6.3 ℃時,糧堆內濕空氣隨微氣流流動加快,遷移速率隨豎向壓力增加而加快。糧堆內濕空氣在邊界處聚集,累積至峰值時少量濕空氣會向糧堆內遷移。

    [1] 吳子丹,趙會義,曹陽,等. 糧食儲藏生態(tài)系統(tǒng)的仿真技術應用研究進展[J]. 糧油食品科技,2014,22(1):1-6.

    Wu Zidan,Zhao Huiyi,Cao Yang, et al. Research progress on the applications of simulation technology in grain storage ecosystem[J]. Science and Technology of Cereals, Oils and Foods, 2014, 22(1): 1-6. (in Chinese with English abstract)

    [2] Multon J L. Spoilage Mechanisms of Grains and Seeds in the Post-harvest Ecosystem, the Resulting Losses and Strategies for the Defense of Stocks[M]. New York: Lavoisier Publishing, 1988.

    [3] 陳桂香,劉超賽,鄭德乾,等. 倉型尺寸對大堆高散裝糧堆空間應力場影響的數值模擬研究[J]. 河南工業(yè)大學學報:自然科學版,2018,39(5):90-95.

    Chen Guixiang, Liu Chaosai, Zheng Deqian, et al. Numerical simulation of warehouse size effect on the stress field of bulk wheat under large pile height[J]. Journal of Henan University of Technology: Natural Science Edition, 2018, 39(5): 90-95. (in Chinese with English abstract)

    [4] 張達,鄭德乾,陳桂香,等. 高大平房倉散裝糧堆壓力場FLAC3D數值模擬[J]. 河南工業(yè)大學學報:自然科學版,2017,38(6):98-103.

    Zhang Da, Zheng Deqian, Chen Guixiang, et al. Numerical simulation on pressure field of bulk grain pile in large warehouse by FLAC3D[J]. Journal of Henan University of Technology: Natural Science Edition, 2017, 38(6): 98-103. (in Chinese with English abstract)

    [5] 王錄民,劉永超,許啟鏗,等. 散裝糧堆底部壓力實測研究[J]. 河南工業(yè)大學學報:自然科學版,2013,34(4):1-4.

    Wang Lumin, Liu Yongchao, Xu Qikeng, et al. Research of site measurement of base pressure of bulk grain pile[J]. Journal of Henan University of Technology: Natural Science Edition, 2013, 34(4): 1-4. (in Chinese with English abstract)

    [6] 陳家豪,陳桂香,王忠旭,等. 基于糧堆有效摩擦系數分布規(guī)律的平底筒倉儲料靜態(tài)壓力場計算方法[J]. 中國農業(yè)大學學報,2019,24(11):152-161.

    Chen Jiahao, Chen Guixiang, Wang Zhongxu, et al. Calculation method for static pressure field of bulk material in silo based on distribution law of effective friction coefficient in grain pile[J]. Journal of China Agricultural University, 2019, 24(11): 152-161. (in Chinese with English abstract)

    [7] 陳家豪. 散裝糧堆空間壓力場理論與試驗研究[D]. 鄭州:河南工業(yè)大學,2017.

    Chen Jiahao. The Theoretical and Experimental Study on Spatial Pressure Field of Granular Grain Pile[D]. Zhengzhou: Henan University of Technology, 2017. (in Chinese with English abstract)

    [8] Mounfield C C, Edwards S F. A theoretical model for the stress distribution in granular matter.Ⅱ.Forces in pipes[J]. Physica A, 1996, 226: 12-24.

    [9] 唐福元,許倩,程緒鐸,等. 筒倉中稻谷的空隙率分布研究[J]. 中國糧油學報,2017,32(12):110-116.

    Tang Fuyuan, Xu Qian, Cheng Xuduo, et al. Research on porosity distribution of paddy in silos[J]. Journal of the Chinese Cereals and Oils Association, 2017, 32(12): 110-116. (in Chinese with English abstract)

    [10] 張小正. 倉儲糧堆孔隙結構參數研究[D]. 北京:中國農業(yè)大學,2010.

    Zhang Xiaozheng. Study on the Pore Parameter of Grain Heap[D]. Beijing: China Agricultural University, 2010. (in Chinese with English abstract)

    [11] 王娟. 糧堆孔隙度的試驗研究[D]. 鄭州:河南工業(yè)大學,2016.

    Wang Juan. Test Study on the Porosity of Grain Mass[D]. Zhengzhou: Henan University of Technology, 2016. (in Chinese with English abstract)

    [12] Jia Canchun, Sun Dawen, Cao Chongwen. Mathematical simulation of temperature fields in a stored grain bin due to internal heat generation[J]. Journal of Food Engineering, 2000, 43: 227-233.

    [13] Thorpe G R. Modelling ecosystems in ventilated conical bottomed farm grain silos[J]. Ecological Modelling, 1997, 94(2/3): 255-286.

    [14] 陳桂香,岳龍飛,王振清,等. 鋼筋混凝土地下糧倉準靜態(tài)溫度場數值模擬[J]. 中國糧油學報,2014,29(3):79-83.

    Chen Guixiang, Yue Longfei, Wang Zhenqing, et al. Quasistatic temperature field simulation of reinforced concrete undergroung warehouse grain[J]. Journal of the Chinese Cereals and Oils Association, 2014, 29(3): 79-83. (in Chinese with English abstract)

    [15] 張忠杰,李瓊,楊德勇,等. 準靜態(tài)倉儲糧堆溫度場的CFD模擬[J]. 中國糧油學報,2010,25(4):46-50.

    Zhang Zhongjie, Li Qiong, Yang Deyong, et al. Temperature field simulation of stored bulk grain by CFD method[J]. Journal of the Chinese Cereals and Oils Association, 2010, 25(4): 46-50. (in Chinese with English abstract)

    [16] 尹君,吳子丹,張忠杰,等. 不同倉型的糧堆溫度場重現及對比分析[J]. 農業(yè)工程學報,2015,31(1):281-287.

    Yin Jun, Wu Zidan, Zhang Zhongjie, et al. Comparison and analysis of temperature field reappearance in stored grain of different warehouses[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(1): 281-287. (in Chinese with English abstract)

    [17] Jian Fuji, Jayas D S, White N D. Temperature fluctuations and moisture migration in wheat stored for 15 months in a metal silo in Canana[J]. Journal of Stored Products Research, 2009, 45(2): 82-90.

    [18] Hammami F, Mabrouk S, Mami A. Modelling and simulation of heat exchange and moisture content in a cereal storage silo[J]. Mathematical and Computer Modelling of Dynamical Systems, 2016, 22(3): 1-14.

    [19] Gaston A R, Abalone R E, Bartosik. Mathematical modelling of heat and moisture transfer of wheat stored in plastic bags[J]. Biosystems Engineering, 2009, 104(1): 72-85.

    [20] 王遠成,白忠權,張中濤,等.倉儲糧堆內熱濕耦合傳遞的數值模擬研究[J]. 中國糧油學報,2015,30(11):97-102.

    Wang Yuancheng, Bai Zhongquan, Zhang Zhongtao, et al. Numerical study on natural convection heat and moisture transfer in stored grain in a silo[J]. Journal of the Chinese Cereals and Oils Association, 2015, 30(11): 97-102. (in Chinese with English abstract)

    [21] 亓偉,王遠成,白忠權,等. 密閉圓筒倉內儲糧自然對流及熱濕耦合傳遞的研究[J]. 中國糧油學報,2015,30(10):83-89.

    Qi Wei, Wang Yuancheng, Bai Zhongquan, et al. Research of closed silo stored grain in the natural convection and the coupled heat and moisture transfer[J]. Journal of the Chinese Cereals and Oils Association, 2015, 30(10): 83-89. (in Chinese with English abstract)

    [22] 王振華. 倉儲糧堆濕熱傳遞過程的數值模擬與試驗研究[D]. 北京:中國農業(yè)大學,2014.

    Wang Zhenhua. Numerical Simulation and Experimental Study on Heat and Mass Transfer of Stored Grain Pile[D]. Beijing: China Agricultural University. (in Chinese with English abstract)

    [23] 尹君. 小麥糧堆多場耦合模型及結露預測研究[D]. 長春:吉林大學,2015.

    Yin Jun. Research on Multi-fields Coupling Model of Wheat Grain and Condensation Prediction[D]. Changchun: Jilin University, 2015. (in Chinese with English abstract)

    [24] 王小萌,吳文福,尹君,等. 基于溫濕度場云圖的小麥糧堆霉變與溫濕度耦合分析[J]. 農業(yè)工程學報,2018,34(10):260-266.

    Wang Xiaomeng, Wu Wenfu, Yin Jun, et al. Analysis of wheat bulk mould and temperature-humidity coupling based on temperature and humidity field cloud map[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(10): 260-266. (in Chinese with English abstract)

    [25] 王小萌,吳文福,尹君,等. 玉米糧堆霉變發(fā)熱過程中的溫濕度場變化規(guī)律研究[J]. 農業(yè)工程學報,2019,35(3):268-273.

    Wang Xiaomeng, Wu Wenfu, Yin Jun, et al. Research on temperature and humidity field change during corn bulk microbiological heating[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(3): 268-273. (in Chinese with English abstract)

    [26] 劉超賽. 散裝小麥糧堆熱-濕-力多場耦合規(guī)律研究[D]. 鄭州:河南工業(yè)大學,2019.

    Liu Chaosai. Study on the Coupling Law of Heat-Wet-Force in Wheat Grain Bulk[D]. Zhengzhou: Henan University of Technology, 2019. (in Chinese with English abstract)

    [27] Rusinek R, Kobylka R. Experimental study and discrete element method modeling of temperature distributions in rapeseed stored in a model bin[J]. Journal of Stored Products Research, 2014, 59: 254-259.

    [28] Khankari K K, Jayas D S, Smith E A, et al. Mathematical model for natural convection moisture migration in stored grain[J]. Transactions of the ASAE, 1995, 38(6): 1077-1087.

    Effects of vertical pressure on moisture and heat transfer of wheat bulk in a granary

    Chen Guixiang, Liu Chaosai, Jiang Minmin, Chen Jiahao, Wang Haitao, Zhang Hongwei

    (,450001,)

    Multiple physical fields, such as pressure, temperature and humidity, can pose grains deterioration in the wheat bulk during storage. It is necessary to explore coupled moisture and heat transfer in wheat bulk under multi-field conditions for the safety of grain storage. In most previous studies, the significant findings are: 1) the porosity decreases with the increase of grain depth; 2) the contact area between grains increases with increasing vertical pressure; 3) the thermal conductivity of wheat bulk is 49 times that of air; and 4) the thermal conductivity between grains decreases with the decrease of porosity. However, an accurate prediction of grain moisture and temperature during storage is still lacking, in order to develop efficient strategies of ventilation conditions. This paper aims to find the influence of vertical pressure on the moisture and heat transfer in the wheat bulk, taking a wheat bulk unit in the silo as the research object, and thereby an experimental study was performed on a multi-field coupling test device. Temperature and humidity of the wheat bulk were measured under three vertical pressure conditions of 50, 100 and 150 kPa, according to the pressure range in large-scale wheat bulk. The temperature was set as high temperature boundary of 38.5℃, low temperature boundary of 5.2℃, and initial grain temperature of 25.8℃. The experimental results show that the grain temperature reduced by 0.5℃ in the vertical pressure of 50 kPa, compared with the initial grain temperature, whereas, it reduced by 1.3℃ when the vertical pressure was 150 kPa. In a constant temperature field, the temperature gradient of wheat bulk decreased from 34.6℃/m to 31.6℃/m, while the rate of change reached 8.7% when the vertical pressure increased from 50 to 150 kPa. The area of high temperature decreased as a power function with storage period of wheat bulk under different vertical pressures, whereas, the power function index increased with the increase of vertical pressures. The temperature difference can cause the wet air in grain bulk to migrate from the higher temperature region to the lower temperature region. In the middle and upper grain bulk with the height greater than 0.3 m, under the action of a large temperature difference between the grain bulk and the high temperature boundary, the wet air in the grain bulk uniformly migrated from the high temperature boundary to the low temperature boundary with the micro airflow. In the middle and lower grain bulk with the height less than 0.3 m, the temperature difference between the two temperature control boundaries was small, and the micro-airflow effect was weak, due to the bottom of the unit cannot be absolutely insulated. The wet air near the high-temperature boundary moved to the low temperature, while the relative humidity was high below near high temperature boundary. A small amount of wet air in the grain bulk can migrate from the low temperature boundary to the wheat bulk after the relative humidity near the low temperature boundary region reached its peak value. When the temperature difference between the middle part of wheat bulk and the boundary near the low temperature was greater than 6.3℃, the migration rate of wet air in grain bulk increased with micro airflow, and the decrease rate of average relative humidity in the middle of wheat bulk increased with the increase of vertical pressure. The findings can provide theoretical support for the research on the evolution of grain temperature and moisture in the multi-field coupling of stored grain bulk under different ventilation conditions.

    temperature; relative humidity; thewheat bulk; vertical pressure; multi-field coupling

    陳桂香,劉超賽,蔣敏敏,等. 儲料豎向壓力對糧倉中小麥糧堆濕熱傳遞的影響[J]. 農業(yè)工程學報,2020,36(8):246-253.doi:10.11975/j.issn.1002-6819.2020.08.030 http://www.tcsae.org

    Chen Guixiang, Liu Chaosai, Jiang Minmin, et al. Effects of vertical pressure on moisture and heat transfer of wheat bulk in a granary[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(8): 246-253. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.08.030 http://www.tcsae.org

    2019-12-19

    2020-04-08

    國家糧食公益性行業(yè)科研專項(201513001);國家自然科學基金項目(51408197);河南省科技廳自然科學項目(182102210387)

    陳桂香,教授,博士,主要從事安全儲糧技術等方面研究。Email:chen-guixiang@163.com

    10.11975/j.issn.1002-6819.2020.08.030

    TU111.4

    A

    1002-6819(2020)-08-0246-08

    猜你喜歡
    高溫區(qū)糧堆儲藏
    雙溫區(qū)脈管制冷機的制冷量主動調控策略
    基于GIS和RS的南京市地表溫度動態(tài)變化研究
    WNS型燃氣鍋爐高溫區(qū)煙管管端裂紋成因及預防措施
    列車制動盤表面溫度分布演變的研究
    機械(2021年8期)2021-09-04 07:24:42
    儲糧壓力對玉米糧堆溫度場影響的實驗與模擬研究
    西藥藥品的儲藏探討
    溫濕度對稻谷糧堆結露的影響及實倉結露預警
    便于儲藏的美味
    科學啟蒙(2016年11期)2016-12-09 19:06:45
    地下庫糧堆充氮技術研究
    現代食品(2016年14期)2016-04-28 08:10:33
    “蘊藏”和“儲藏”
    他把我摸到了高潮在线观看| 亚洲专区国产一区二区| 法律面前人人平等表现在哪些方面| 无人区码免费观看不卡| 国产亚洲av嫩草精品影院| 精品免费久久久久久久清纯| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av香蕉五月| 激情在线观看视频在线高清| 亚洲av成人一区二区三| 99热这里只有是精品50| 日本三级黄在线观看| 久久久久国产一级毛片高清牌| 国产成人aa在线观看| 成人无遮挡网站| 男女视频在线观看网站免费| 黄色视频,在线免费观看| 国产精品自产拍在线观看55亚洲| 长腿黑丝高跟| 深夜精品福利| 日本 av在线| 欧美日本视频| 校园春色视频在线观看| 精品久久蜜臀av无| 欧美黄色淫秽网站| 99热6这里只有精品| 两性夫妻黄色片| 麻豆av在线久日| 亚洲国产精品成人综合色| 国内揄拍国产精品人妻在线| 午夜影院日韩av| 亚洲 国产 在线| 成人av一区二区三区在线看| 亚洲欧美精品综合一区二区三区| 欧美乱色亚洲激情| 脱女人内裤的视频| 国产av一区在线观看免费| 国产成人av教育| 免费搜索国产男女视频| 色播亚洲综合网| 狠狠狠狠99中文字幕| 免费观看人在逋| 国产精品av久久久久免费| 久久久久久大精品| 亚洲中文av在线| 国产成年人精品一区二区| 久久久久久久久久黄片| 97碰自拍视频| 在线国产一区二区在线| 亚洲男人的天堂狠狠| 午夜激情福利司机影院| 中文字幕精品亚洲无线码一区| 国产成人aa在线观看| 国产精华一区二区三区| 久久中文字幕人妻熟女| 国内毛片毛片毛片毛片毛片| 91av网站免费观看| 18禁美女被吸乳视频| 亚洲欧美一区二区三区黑人| 搡老妇女老女人老熟妇| 一夜夜www| 国产91精品成人一区二区三区| 色综合亚洲欧美另类图片| 精品国产美女av久久久久小说| 久久久成人免费电影| 国产精品免费一区二区三区在线| 国产主播在线观看一区二区| 人人妻,人人澡人人爽秒播| 午夜福利高清视频| 此物有八面人人有两片| 神马国产精品三级电影在线观看| 亚洲 国产 在线| 欧美又色又爽又黄视频| 欧美日韩精品网址| 日韩欧美在线乱码| 国产在线精品亚洲第一网站| www.精华液| 美女扒开内裤让男人捅视频| 国产69精品久久久久777片 | 美女扒开内裤让男人捅视频| 成人av在线播放网站| 午夜免费激情av| 最近最新中文字幕大全免费视频| 国产一区二区激情短视频| 久久久精品欧美日韩精品| 亚洲精品美女久久久久99蜜臀| 又粗又爽又猛毛片免费看| 天天一区二区日本电影三级| 亚洲av片天天在线观看| 99国产精品99久久久久| 免费观看精品视频网站| 两个人看的免费小视频| 99热这里只有是精品50| 99热这里只有是精品50| 久久热在线av| 蜜桃久久精品国产亚洲av| 观看免费一级毛片| 一个人看的www免费观看视频| 此物有八面人人有两片| 1024香蕉在线观看| 69av精品久久久久久| 亚洲色图 男人天堂 中文字幕| 精品99又大又爽又粗少妇毛片 | 国产视频一区二区在线看| 成人高潮视频无遮挡免费网站| 国产精品野战在线观看| xxx96com| 亚洲av片天天在线观看| 国产高清videossex| 国产成人av激情在线播放| 欧美一级毛片孕妇| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一电影网av| 在线播放国产精品三级| 99热这里只有精品一区 | 少妇熟女aⅴ在线视频| 美女被艹到高潮喷水动态| 国内精品久久久久精免费| 怎么达到女性高潮| 黑人操中国人逼视频| 久久精品亚洲精品国产色婷小说| 少妇熟女aⅴ在线视频| 无限看片的www在线观看| 国产精品99久久久久久久久| 亚洲国产精品合色在线| 99视频精品全部免费 在线 | 亚洲精品美女久久av网站| 五月玫瑰六月丁香| 国产欧美日韩精品亚洲av| 亚洲国产欧美人成| 日本 av在线| 国产蜜桃级精品一区二区三区| 久久99热这里只有精品18| 国产精品,欧美在线| av欧美777| 国产又色又爽无遮挡免费看| 日本与韩国留学比较| 久久人人精品亚洲av| 欧美在线黄色| www.自偷自拍.com| 此物有八面人人有两片| 久久欧美精品欧美久久欧美| 国产av在哪里看| 国内揄拍国产精品人妻在线| 淫秽高清视频在线观看| 亚洲成a人片在线一区二区| 天堂av国产一区二区熟女人妻| 精品人妻1区二区| 亚洲精品美女久久av网站| 亚洲欧美日韩高清专用| 看黄色毛片网站| www.自偷自拍.com| 91在线精品国自产拍蜜月 | 国产欧美日韩一区二区精品| 国产激情偷乱视频一区二区| 黑人操中国人逼视频| 中文字幕av在线有码专区| 国产精品电影一区二区三区| 真人做人爱边吃奶动态| 亚洲最大成人中文| 免费一级毛片在线播放高清视频| 少妇熟女aⅴ在线视频| 我的老师免费观看完整版| www.自偷自拍.com| 国产真人三级小视频在线观看| 97超视频在线观看视频| 免费高清视频大片| 日韩欧美免费精品| 国产精品久久视频播放| 日韩有码中文字幕| 免费看日本二区| 夜夜躁狠狠躁天天躁| 亚洲专区字幕在线| 日本成人三级电影网站| 男人舔奶头视频| 日本 av在线| АⅤ资源中文在线天堂| 久久久精品欧美日韩精品| 91久久精品国产一区二区成人 | 免费一级毛片在线播放高清视频| 国产免费男女视频| 九九在线视频观看精品| 久久欧美精品欧美久久欧美| 波多野结衣高清无吗| 亚洲成人中文字幕在线播放| 国产精品av视频在线免费观看| 天堂√8在线中文| 久久99热这里只有精品18| 日韩欧美一区二区三区在线观看| 99久久精品热视频| 天堂√8在线中文| 国产真人三级小视频在线观看| 中文字幕人妻丝袜一区二区| 丝袜人妻中文字幕| 中出人妻视频一区二区| 一区二区三区国产精品乱码| 国产高清视频在线播放一区| 亚洲精品色激情综合| 两个人视频免费观看高清| aaaaa片日本免费| 99国产精品99久久久久| 999久久久国产精品视频| 床上黄色一级片| 麻豆成人午夜福利视频| 国内精品久久久久精免费| 九九在线视频观看精品| 国产三级在线视频| 国产精品日韩av在线免费观看| 午夜免费激情av| 国产又黄又爽又无遮挡在线| 一级黄色大片毛片| 亚洲美女视频黄频| cao死你这个sao货| 91麻豆av在线| 真人做人爱边吃奶动态| 黄色成人免费大全| 美女午夜性视频免费| 两个人视频免费观看高清| 91麻豆av在线| 亚洲自偷自拍图片 自拍| av国产免费在线观看| 露出奶头的视频| 久久精品91无色码中文字幕| 18禁美女被吸乳视频| 亚洲av熟女| 别揉我奶头~嗯~啊~动态视频| 蜜桃久久精品国产亚洲av| 一区二区三区国产精品乱码| 99视频精品全部免费 在线 | 中国美女看黄片| 国产欧美日韩一区二区精品| 美女扒开内裤让男人捅视频| 看黄色毛片网站| 成年女人看的毛片在线观看| 国产三级在线视频| 亚洲国产欧洲综合997久久,| 99国产精品一区二区蜜桃av| 久久午夜亚洲精品久久| 国产激情久久老熟女| 人人妻人人看人人澡| 精品国产超薄肉色丝袜足j| 午夜福利欧美成人| 中文字幕人妻丝袜一区二区| 亚洲欧美一区二区三区黑人| 女警被强在线播放| 精品久久久久久久久久久久久| 美女被艹到高潮喷水动态| 麻豆成人av在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品免费一区二区三区在线| 最近最新中文字幕大全免费视频| 亚洲人与动物交配视频| 男女那种视频在线观看| 国产精品国产高清国产av| 国产极品精品免费视频能看的| 久久精品人妻少妇| 日韩精品青青久久久久久| 搡老岳熟女国产| 成人三级黄色视频| 国产三级黄色录像| 在线观看午夜福利视频| 欧美一级毛片孕妇| 特级一级黄色大片| 此物有八面人人有两片| 国产真人三级小视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 床上黄色一级片| 日韩高清综合在线| 成人鲁丝片一二三区免费| 国产亚洲精品av在线| 不卡一级毛片| 一个人看的www免费观看视频| 91字幕亚洲| 五月玫瑰六月丁香| a级毛片在线看网站| 精品日产1卡2卡| 丁香六月欧美| 亚洲av成人不卡在线观看播放网| 亚洲av第一区精品v没综合| 成熟少妇高潮喷水视频| 波多野结衣巨乳人妻| 夜夜看夜夜爽夜夜摸| 校园春色视频在线观看| 久久久久国产一级毛片高清牌| 午夜久久久久精精品| 国产伦精品一区二区三区四那| 国产乱人伦免费视频| 99国产极品粉嫩在线观看| 亚洲成人中文字幕在线播放| 91av网一区二区| 亚洲熟妇中文字幕五十中出| 长腿黑丝高跟| 嫩草影视91久久| 欧美最黄视频在线播放免费| 日本成人三级电影网站| 国产精品综合久久久久久久免费| 在线免费观看的www视频| 亚洲精品中文字幕一二三四区| 亚洲成人久久爱视频| 夜夜躁狠狠躁天天躁| 成人特级黄色片久久久久久久| 少妇丰满av| 欧美性猛交黑人性爽| a在线观看视频网站| 欧美乱妇无乱码| 欧美色欧美亚洲另类二区| 亚洲av片天天在线观看| 精品电影一区二区在线| 国产成年人精品一区二区| 精品久久久久久久末码| 成年女人毛片免费观看观看9| 国产熟女xx| 久久精品人妻少妇| av国产免费在线观看| 变态另类丝袜制服| 国产高清有码在线观看视频| 免费看十八禁软件| 亚洲欧洲精品一区二区精品久久久| 无限看片的www在线观看| 日本免费一区二区三区高清不卡| 美女高潮喷水抽搐中文字幕| 中国美女看黄片| 一本精品99久久精品77| 国语自产精品视频在线第100页| 少妇丰满av| 日韩av在线大香蕉| www.精华液| 又黄又爽又免费观看的视频| 一级a爱片免费观看的视频| 久久久精品欧美日韩精品| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添小说| 身体一侧抽搐| 他把我摸到了高潮在线观看| 久久天躁狠狠躁夜夜2o2o| 国产激情久久老熟女| 激情在线观看视频在线高清| 天天躁狠狠躁夜夜躁狠狠躁| 欧美性猛交╳xxx乱大交人| 国产成人aa在线观看| 国产伦精品一区二区三区四那| 校园春色视频在线观看| 亚洲第一电影网av| 小说图片视频综合网站| 国产精品久久久久久精品电影| 久久亚洲真实| 精品熟女少妇八av免费久了| 中文字幕久久专区| 国产高清视频在线播放一区| 亚洲中文av在线| 国产97色在线日韩免费| 欧美日韩亚洲国产一区二区在线观看| 欧美高清成人免费视频www| 免费看a级黄色片| 极品教师在线免费播放| 久久久国产成人精品二区| 精品一区二区三区av网在线观看| 小说图片视频综合网站| 黄片小视频在线播放| 毛片女人毛片| 国产av一区在线观看免费| 精品国产超薄肉色丝袜足j| 亚洲狠狠婷婷综合久久图片| 搡老妇女老女人老熟妇| 桃红色精品国产亚洲av| 久久精品国产综合久久久| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美日韩高清专用| 天天一区二区日本电影三级| 亚洲av成人av| 成人亚洲精品av一区二区| 国产极品精品免费视频能看的| 亚洲狠狠婷婷综合久久图片| 午夜久久久久精精品| 欧美日本亚洲视频在线播放| 美女cb高潮喷水在线观看 | 免费在线观看影片大全网站| 男女之事视频高清在线观看| 精品一区二区三区av网在线观看| 五月伊人婷婷丁香| 别揉我奶头~嗯~啊~动态视频| 精品福利观看| 听说在线观看完整版免费高清| 真实男女啪啪啪动态图| 搡老岳熟女国产| 国产美女午夜福利| av天堂中文字幕网| 最近视频中文字幕2019在线8| 国产高清videossex| 午夜免费成人在线视频| 毛片女人毛片| 又黄又爽又免费观看的视频| 午夜福利免费观看在线| 中文字幕最新亚洲高清| 午夜福利18| 国产欧美日韩一区二区精品| 亚洲自偷自拍图片 自拍| 亚洲激情在线av| 国产视频内射| 亚洲中文字幕日韩| 精品国产乱码久久久久久男人| 91老司机精品| 少妇丰满av| 两人在一起打扑克的视频| 一级黄色大片毛片| 天堂√8在线中文| 国产美女午夜福利| 欧美日韩国产亚洲二区| 香蕉av资源在线| 精品午夜福利视频在线观看一区| 久久久久久人人人人人| 亚洲aⅴ乱码一区二区在线播放| 看免费av毛片| 久久久成人免费电影| 色av中文字幕| 日韩精品青青久久久久久| 一本久久中文字幕| 岛国在线观看网站| 噜噜噜噜噜久久久久久91| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| 久久久久性生活片| 女同久久另类99精品国产91| 天天添夜夜摸| 国产高清videossex| 日韩精品青青久久久久久| 免费搜索国产男女视频| 一级毛片精品| 久久这里只有精品中国| 国产成人av教育| 亚洲电影在线观看av| 国产黄a三级三级三级人| 丰满人妻熟妇乱又伦精品不卡| 久久国产乱子伦精品免费另类| 日韩有码中文字幕| 成年女人看的毛片在线观看| 久久亚洲真实| 国产伦一二天堂av在线观看| 国产欧美日韩精品亚洲av| 女人被狂操c到高潮| 好男人电影高清在线观看| 男插女下体视频免费在线播放| 亚洲无线观看免费| 国产91精品成人一区二区三区| 国产69精品久久久久777片 | 色播亚洲综合网| 人妻丰满熟妇av一区二区三区| 亚洲国产看品久久| 国产成人av教育| 悠悠久久av| 色播亚洲综合网| 搞女人的毛片| 叶爱在线成人免费视频播放| 久久久久久久精品吃奶| 91老司机精品| 色吧在线观看| 成人亚洲精品av一区二区| 亚洲精品一区av在线观看| 久久久国产成人免费| 国产精品98久久久久久宅男小说| 两个人视频免费观看高清| 在线观看舔阴道视频| 日韩三级视频一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 我的老师免费观看完整版| 老鸭窝网址在线观看| 精品久久久久久久久久免费视频| 男女做爰动态图高潮gif福利片| 蜜桃久久精品国产亚洲av| 精品乱码久久久久久99久播| 色综合婷婷激情| 香蕉久久夜色| 亚洲精品中文字幕一二三四区| 日本成人三级电影网站| 国产麻豆成人av免费视频| 男人舔奶头视频| 很黄的视频免费| ponron亚洲| 给我免费播放毛片高清在线观看| 一级黄色大片毛片| a级毛片a级免费在线| 精品国产乱子伦一区二区三区| 久久亚洲精品不卡| 欧美中文日本在线观看视频| 好男人电影高清在线观看| 亚洲五月婷婷丁香| 精品福利观看| 在线观看美女被高潮喷水网站 | 18禁黄网站禁片免费观看直播| xxxwww97欧美| 日韩欧美国产一区二区入口| 高清毛片免费观看视频网站| 亚洲国产精品合色在线| 麻豆国产97在线/欧美| 老汉色∧v一级毛片| 中亚洲国语对白在线视频| 99国产精品一区二区三区| 露出奶头的视频| 99久久99久久久精品蜜桃| 国产一区二区三区在线臀色熟女| 亚洲 欧美一区二区三区| 九色成人免费人妻av| 怎么达到女性高潮| www.999成人在线观看| 国产成人av激情在线播放| 中文字幕人妻丝袜一区二区| xxx96com| 欧美黑人欧美精品刺激| 中文在线观看免费www的网站| 久久天堂一区二区三区四区| 午夜久久久久精精品| 1024香蕉在线观看| 中文字幕最新亚洲高清| 91麻豆精品激情在线观看国产| 级片在线观看| 宅男免费午夜| 真实男女啪啪啪动态图| 五月玫瑰六月丁香| 国内精品久久久久精免费| 美女高潮喷水抽搐中文字幕| 国产成人av激情在线播放| 欧美黑人欧美精品刺激| 国产在线精品亚洲第一网站| 在线免费观看的www视频| 真实男女啪啪啪动态图| 99精品欧美一区二区三区四区| 天堂网av新在线| 男人舔女人下体高潮全视频| 成熟少妇高潮喷水视频| 在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 少妇丰满av| 国内久久婷婷六月综合欲色啪| 他把我摸到了高潮在线观看| 中文在线观看免费www的网站| 国产私拍福利视频在线观看| 中文字幕最新亚洲高清| 日韩欧美三级三区| 欧美最黄视频在线播放免费| 成年女人永久免费观看视频| 成人精品一区二区免费| 国产精品自产拍在线观看55亚洲| 亚洲成av人片在线播放无| 我要搜黄色片| 激情在线观看视频在线高清| 一边摸一边抽搐一进一小说| 夜夜看夜夜爽夜夜摸| 欧美色视频一区免费| 日韩免费av在线播放| 很黄的视频免费| av国产免费在线观看| 亚洲国产精品久久男人天堂| 一进一出好大好爽视频| 99热这里只有是精品50| 久久久水蜜桃国产精品网| 波多野结衣巨乳人妻| 极品教师在线免费播放| 天天一区二区日本电影三级| 一个人观看的视频www高清免费观看 | 青草久久国产| 一区二区三区高清视频在线| 午夜激情欧美在线| 亚洲美女黄片视频| 亚洲欧美精品综合久久99| 草草在线视频免费看| 伦理电影免费视频| 国产91精品成人一区二区三区| 精品一区二区三区av网在线观看| 久久久水蜜桃国产精品网| 小蜜桃在线观看免费完整版高清| 亚洲无线在线观看| 麻豆一二三区av精品| 黄色视频,在线免费观看| 不卡av一区二区三区| 亚洲精华国产精华精| 午夜福利在线观看免费完整高清在 | 欧美成人免费av一区二区三区| 婷婷六月久久综合丁香| 天堂网av新在线| or卡值多少钱| 亚洲自偷自拍图片 自拍| a级毛片在线看网站| 国产午夜福利久久久久久| 亚洲人与动物交配视频| 亚洲美女视频黄频| 午夜精品一区二区三区免费看| 久久久久久国产a免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 黄片大片在线免费观看| 国产爱豆传媒在线观看| 国产69精品久久久久777片 | www.精华液| 深夜精品福利| 色综合婷婷激情| 日韩av在线大香蕉| 久久精品91蜜桃| 午夜激情欧美在线| 国产一区在线观看成人免费| 国产精品 欧美亚洲| av欧美777| 久久热在线av| 天天添夜夜摸| 国产精品1区2区在线观看.| tocl精华| 2021天堂中文幕一二区在线观| 精品久久久久久成人av| 免费电影在线观看免费观看| 国产精品电影一区二区三区| 日韩欧美免费精品| 九九在线视频观看精品| 可以在线观看的亚洲视频| 亚洲自拍偷在线| www.熟女人妻精品国产| 亚洲精品美女久久av网站| 一本综合久久免费|