• <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一区二区三区| 99热6这里只有精品| 久99久视频精品免费| 久久久精品国产亚洲av高清涩受| 99久久精品国产亚洲精品| 欧美黄色淫秽网站| 国产人伦9x9x在线观看| 欧美另类亚洲清纯唯美| 欧美成人免费av一区二区三区| 国产精品永久免费网站| 黄片播放在线免费| 亚洲av中文字字幕乱码综合 | 午夜两性在线视频| 色哟哟哟哟哟哟| 嫩草影视91久久| 午夜福利欧美成人| 岛国视频午夜一区免费看| 亚洲熟妇熟女久久| 亚洲狠狠婷婷综合久久图片| 老汉色av国产亚洲站长工具| 国产成年人精品一区二区| 啦啦啦 在线观看视频| 淫秽高清视频在线观看| 亚洲精品中文字幕一二三四区| 最近最新免费中文字幕在线| 久久九九热精品免费| 国产真实乱freesex| 国内精品久久久久久久电影| 亚洲精品国产一区二区精华液| www国产在线视频色| 在线播放国产精品三级| 脱女人内裤的视频| 一本综合久久免费| 777久久人妻少妇嫩草av网站| 999精品在线视频| 最好的美女福利视频网| 日韩欧美一区视频在线观看| 日本精品一区二区三区蜜桃| 国产精品爽爽va在线观看网站 | 两人在一起打扑克的视频| 99热这里只有精品一区 | 午夜影院日韩av| 黄色女人牲交| 国产一区二区激情短视频| 18禁美女被吸乳视频| 精品免费久久久久久久清纯| 国产精品久久视频播放| 真人一进一出gif抽搐免费| 人人妻人人澡人人看| 久久久久国产一级毛片高清牌| 精品久久久久久久久久免费视频| 国产在线精品亚洲第一网站| 亚洲欧美日韩高清在线视频| 50天的宝宝边吃奶边哭怎么回事| 国产精品亚洲美女久久久| 美国免费a级毛片| 欧美日韩精品网址| 欧美日韩黄片免| 亚洲熟女毛片儿| 久久香蕉精品热| 国产精品综合久久久久久久免费| 国产一区二区在线av高清观看| 国产精品久久久人人做人人爽| 欧美激情极品国产一区二区三区| 国内少妇人妻偷人精品xxx网站 | 日韩 欧美 亚洲 中文字幕| 亚洲国产看品久久| 亚洲久久久国产精品| 听说在线观看完整版免费高清| 成人午夜高清在线视频 | 不卡一级毛片| 亚洲美女黄片视频| 在线免费观看的www视频| 国产精品精品国产色婷婷| 久久精品91无色码中文字幕| 国产麻豆成人av免费视频| 日韩欧美在线二视频| 一级毛片高清免费大全| 岛国在线观看网站| 亚洲成人久久性| 午夜福利18| 久久中文看片网| 亚洲av第一区精品v没综合| 午夜精品在线福利| 亚洲欧美激情综合另类| www.熟女人妻精品国产| 在线国产一区二区在线| 国产精品一区二区精品视频观看| 天天躁夜夜躁狠狠躁躁| 黄色女人牲交| 精品国产乱码久久久久久男人| 大香蕉久久成人网| av片东京热男人的天堂| 狂野欧美激情性xxxx| 亚洲欧洲精品一区二区精品久久久| 99精品在免费线老司机午夜| 波多野结衣巨乳人妻| 国产精品99久久99久久久不卡| 动漫黄色视频在线观看| 日韩成人在线观看一区二区三区| 在线观看一区二区三区| 亚洲久久久国产精品| 免费女性裸体啪啪无遮挡网站| 国产亚洲av嫩草精品影院| 亚洲欧美激情综合另类| 国产精品 欧美亚洲| 丰满人妻熟妇乱又伦精品不卡| 免费一级毛片在线播放高清视频| 激情在线观看视频在线高清| 久久久久久九九精品二区国产 | 人人妻人人澡欧美一区二区| 国产精品香港三级国产av潘金莲| 精品久久久久久,| 国产av不卡久久| 国产精品久久视频播放| 97超级碰碰碰精品色视频在线观看| 一本精品99久久精品77| 国产熟女xx| 欧美成人午夜精品| 妹子高潮喷水视频| 久热这里只有精品99| 免费在线观看完整版高清| 亚洲成av片中文字幕在线观看| 免费高清在线观看日韩| 99久久无色码亚洲精品果冻| 两性夫妻黄色片| 国产区一区二久久| 午夜久久久久精精品| 无限看片的www在线观看| 欧美一区二区精品小视频在线| 国产精品久久久久久精品电影 | 久久国产亚洲av麻豆专区| 久久国产精品影院| 色在线成人网| 亚洲av电影在线进入| 91国产中文字幕| 久久久久精品国产欧美久久久| 757午夜福利合集在线观看| 亚洲专区中文字幕在线| 亚洲精品中文字幕一二三四区| www.熟女人妻精品国产| 老汉色∧v一级毛片| 一个人免费在线观看的高清视频| 麻豆国产av国片精品| 两性夫妻黄色片| tocl精华| 国产单亲对白刺激| 黄色毛片三级朝国网站| 美女国产高潮福利片在线看| 亚洲一区二区三区色噜噜| 美女免费视频网站| 成人免费观看视频高清| 午夜免费激情av| 亚洲av电影不卡..在线观看| 一本综合久久免费| 亚洲中文字幕日韩| 欧美av亚洲av综合av国产av| 亚洲精品一区av在线观看| 色婷婷久久久亚洲欧美| 波多野结衣高清无吗| 国产精品一区二区三区四区久久 | 十分钟在线观看高清视频www| 一级毛片精品| 精品一区二区三区视频在线观看免费| 一区福利在线观看| 日韩欧美免费精品| 国产精品美女特级片免费视频播放器 | 久久青草综合色| 一级黄色大片毛片| 99re在线观看精品视频| 很黄的视频免费| 高清在线国产一区| 中文字幕高清在线视频| 国产99白浆流出| 欧美午夜高清在线| 免费女性裸体啪啪无遮挡网站| 中文字幕精品亚洲无线码一区 | 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 午夜视频精品福利| 免费在线观看完整版高清| 夜夜爽天天搞| 色av中文字幕| 日韩国内少妇激情av| 此物有八面人人有两片| 长腿黑丝高跟| www.www免费av| 国产精品久久电影中文字幕| 韩国av一区二区三区四区| 老司机靠b影院| 亚洲男人天堂网一区| 最近最新免费中文字幕在线| 欧美丝袜亚洲另类 | 久久久久久九九精品二区国产 | 欧美日韩中文字幕国产精品一区二区三区| 国产熟女午夜一区二区三区| 两个人看的免费小视频| 午夜激情av网站| 国产不卡一卡二| 亚洲黑人精品在线| 亚洲欧美精品综合久久99| 亚洲国产欧美网| 日韩精品青青久久久久久| 黄色视频不卡| 午夜精品在线福利| 久久热在线av| 日本黄色视频三级网站网址| 香蕉av资源在线| a在线观看视频网站| 精品久久久久久久久久久久久 | avwww免费| 国产亚洲精品久久久久5区| 精品熟女少妇八av免费久了| 黑人欧美特级aaaaaa片| 欧美不卡视频在线免费观看 | 午夜a级毛片| 久久人妻福利社区极品人妻图片| 亚洲人成伊人成综合网2020| 中文字幕高清在线视频| 久久欧美精品欧美久久欧美| 淫秽高清视频在线观看| 免费在线观看成人毛片| 亚洲中文日韩欧美视频| 伊人久久大香线蕉亚洲五| 国产亚洲欧美在线一区二区| 人人妻人人澡人人看| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 露出奶头的视频| 免费在线观看亚洲国产| 国产成+人综合+亚洲专区| 午夜福利18| 久久精品国产亚洲av高清一级| 特大巨黑吊av在线直播 | а√天堂www在线а√下载| 婷婷丁香在线五月| 少妇熟女aⅴ在线视频| 美女午夜性视频免费| 欧美日韩亚洲综合一区二区三区_| 久99久视频精品免费| 人人妻人人看人人澡| 久久久国产欧美日韩av| 免费电影在线观看免费观看| 一边摸一边抽搐一进一小说| 成人三级做爰电影| 777久久人妻少妇嫩草av网站| 丰满的人妻完整版| 青草久久国产| 国产精品电影一区二区三区| 日韩精品青青久久久久久| 国产亚洲精品久久久久5区| 成人亚洲精品av一区二区| 亚洲人成网站在线播放欧美日韩| 俄罗斯特黄特色一大片| 成人欧美大片| 男女床上黄色一级片免费看| 成熟少妇高潮喷水视频| 国产一区二区三区视频了| 亚洲精品国产精品久久久不卡| 香蕉国产在线看| 一级毛片女人18水好多| 亚洲专区中文字幕在线| 99国产综合亚洲精品| 99久久综合精品五月天人人| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产清高在天天线| 俺也久久电影网| 窝窝影院91人妻| 亚洲真实伦在线观看| 一夜夜www| √禁漫天堂资源中文www| 日韩中文字幕欧美一区二区| 女生性感内裤真人,穿戴方法视频| 中文字幕精品亚洲无线码一区 | 国产成人系列免费观看| 女警被强在线播放| 一区二区三区高清视频在线| 国产av又大| 黄色成人免费大全| 母亲3免费完整高清在线观看| 99久久99久久久精品蜜桃| 中文字幕人妻丝袜一区二区| 法律面前人人平等表现在哪些方面| 婷婷精品国产亚洲av在线| 麻豆国产av国片精品| 久久国产精品影院| 啪啪无遮挡十八禁网站| 亚洲专区字幕在线| 一个人免费在线观看的高清视频| 亚洲精品av麻豆狂野| 国产一卡二卡三卡精品| 亚洲精品国产一区二区精华液| 女人高潮潮喷娇喘18禁视频| 无人区码免费观看不卡| 欧美亚洲日本最大视频资源| 久久九九热精品免费| 十八禁人妻一区二区| 岛国在线观看网站| 最近最新中文字幕大全免费视频| 在线观看舔阴道视频| 久久精品影院6| 天堂影院成人在线观看| 国产蜜桃级精品一区二区三区| 中出人妻视频一区二区| 美女免费视频网站| 亚洲第一欧美日韩一区二区三区| 最好的美女福利视频网| 色精品久久人妻99蜜桃| 怎么达到女性高潮| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品av久久久久免费| 亚洲成国产人片在线观看| 国产精品一区二区精品视频观看| 亚洲久久久国产精品| 日韩一卡2卡3卡4卡2021年| 欧美日韩瑟瑟在线播放| 香蕉av资源在线| 在线观看免费午夜福利视频| 久久久久免费精品人妻一区二区 | 一夜夜www| 黑人欧美特级aaaaaa片| 欧美+亚洲+日韩+国产| 亚洲欧美一区二区三区黑人| 黄色丝袜av网址大全| 国产精品亚洲一级av第二区| 国产午夜精品久久久久久| 在线av久久热| 亚洲五月色婷婷综合| 成人免费观看视频高清| 午夜久久久在线观看| 亚洲成人免费电影在线观看| 亚洲无线在线观看| 美女扒开内裤让男人捅视频| 久久久久久免费高清国产稀缺| 亚洲狠狠婷婷综合久久图片| 国产精品 欧美亚洲| 亚洲全国av大片| 色精品久久人妻99蜜桃| 成人一区二区视频在线观看| 在线视频色国产色| 青草久久国产| 亚洲一区二区三区不卡视频| 国产一卡二卡三卡精品| 亚洲成人久久爱视频| 亚洲国产精品合色在线| 免费高清视频大片| 婷婷丁香在线五月| 日韩欧美 国产精品| 国产亚洲欧美在线一区二区| 成人亚洲精品av一区二区| 老汉色∧v一级毛片| 男女之事视频高清在线观看| 精品国产乱子伦一区二区三区| 午夜影院日韩av| 亚洲五月婷婷丁香| 无限看片的www在线观看| 国产成年人精品一区二区| 美女 人体艺术 gogo| 亚洲中文日韩欧美视频| 一a级毛片在线观看| 最近在线观看免费完整版| 午夜日韩欧美国产| 男女下面进入的视频免费午夜 | 国产1区2区3区精品| 美女午夜性视频免费| 麻豆成人午夜福利视频| 曰老女人黄片| 女生性感内裤真人,穿戴方法视频| 校园春色视频在线观看| 91字幕亚洲| 久久久久国内视频| 9191精品国产免费久久| 宅男免费午夜| 久久精品国产99精品国产亚洲性色| 欧美日韩黄片免| 欧美成人性av电影在线观看| 亚洲 欧美一区二区三区| 午夜福利欧美成人| 18禁黄网站禁片午夜丰满| 欧美精品亚洲一区二区| 色综合站精品国产| 国产久久久一区二区三区| 老司机靠b影院| 51午夜福利影视在线观看| 国产成人av激情在线播放| 搡老岳熟女国产| 精品熟女少妇八av免费久了| 不卡av一区二区三区| 男女床上黄色一级片免费看| 可以免费在线观看a视频的电影网站| 一区二区三区激情视频| 中文字幕另类日韩欧美亚洲嫩草| 男女那种视频在线观看| 久久国产精品男人的天堂亚洲| 真人一进一出gif抽搐免费| 美女国产高潮福利片在线看| 亚洲无线在线观看| 九色国产91popny在线| 亚洲第一av免费看| 欧美日韩福利视频一区二区| av免费在线观看网站| 免费高清在线观看日韩| 久久国产精品人妻蜜桃| 在线视频色国产色| 欧美日韩一级在线毛片| 欧美精品亚洲一区二区| 叶爱在线成人免费视频播放| 免费在线观看日本一区| 国产蜜桃级精品一区二区三区| 99久久精品国产亚洲精品| 不卡一级毛片| 国产精品 国内视频| 国产成+人综合+亚洲专区| 久久久水蜜桃国产精品网| 成人18禁在线播放| 亚洲欧美日韩无卡精品| 欧美日韩福利视频一区二区| 在线观看免费视频日本深夜| 国产黄色小视频在线观看| 成人一区二区视频在线观看| 国产精品久久久久久精品电影 | 欧美日本亚洲视频在线播放| 午夜精品在线福利| 国产精品一区二区精品视频观看| 嫩草影院精品99| 久久青草综合色| 午夜免费观看网址| 成人亚洲精品av一区二区| 国产精品久久久久久亚洲av鲁大| 亚洲成av人片免费观看| 免费无遮挡裸体视频| 久久中文字幕一级| 丁香六月欧美| 美女高潮喷水抽搐中文字幕| 日韩欧美在线二视频| 日韩成人在线观看一区二区三区| 99久久国产精品久久久| 色哟哟哟哟哟哟| 日韩高清综合在线| 亚洲色图av天堂| 黄频高清免费视频| 色播在线永久视频| 国产av不卡久久| 最新美女视频免费是黄的| 久久精品夜夜夜夜夜久久蜜豆 | 国产欧美日韩一区二区精品| 日本成人三级电影网站| 制服丝袜大香蕉在线| 亚洲av片天天在线观看| 91麻豆av在线| 搞女人的毛片| 亚洲av电影在线进入| 久久青草综合色| 亚洲,欧美精品.| 法律面前人人平等表现在哪些方面| 亚洲自拍偷在线| 99热6这里只有精品| 一区二区三区高清视频在线| 国产精华一区二区三区| 婷婷精品国产亚洲av| 18禁美女被吸乳视频| 91在线观看av| 视频区欧美日本亚洲| 欧美最黄视频在线播放免费| 免费电影在线观看免费观看| 国产精品 欧美亚洲| 精品国产一区二区三区四区第35| 久久久久久亚洲精品国产蜜桃av| 精品电影一区二区在线| 91国产中文字幕| 一边摸一边做爽爽视频免费| 亚洲国产欧美一区二区综合| 少妇熟女aⅴ在线视频| 精品免费久久久久久久清纯| 亚洲专区字幕在线| 丝袜美腿诱惑在线| 亚洲国产日韩欧美精品在线观看 | 欧美成狂野欧美在线观看| 亚洲午夜精品一区,二区,三区| 国产一区二区激情短视频| 精品一区二区三区av网在线观看| 天天添夜夜摸| 少妇 在线观看| 免费在线观看黄色视频的| 美女大奶头视频| 国产高清videossex| 亚洲中文av在线| www日本在线高清视频| 丝袜人妻中文字幕| 久久久久亚洲av毛片大全| 亚洲国产中文字幕在线视频| 99在线人妻在线中文字幕| 无限看片的www在线观看| 这个男人来自地球电影免费观看| 97碰自拍视频| 国产激情久久老熟女| АⅤ资源中文在线天堂| 国产精品乱码一区二三区的特点| 午夜精品在线福利| 亚洲欧美日韩无卡精品| 欧美久久黑人一区二区| 久久人妻福利社区极品人妻图片| 日韩精品中文字幕看吧| 日本 欧美在线| 欧美一区二区精品小视频在线| 老司机午夜福利在线观看视频| 91麻豆精品激情在线观看国产| 99国产精品99久久久久| 一级片免费观看大全| 黄色 视频免费看| 性欧美人与动物交配| 1024香蕉在线观看| 在线观看舔阴道视频| 一边摸一边抽搐一进一小说| 神马国产精品三级电影在线观看 | 国产乱人伦免费视频| 午夜精品久久久久久毛片777| 97人妻精品一区二区三区麻豆 | av在线播放免费不卡| 免费在线观看日本一区| 国产成人精品久久二区二区91| 色哟哟哟哟哟哟| 成年女人毛片免费观看观看9| 亚洲成人国产一区在线观看| 中文资源天堂在线| 又大又爽又粗| 国产亚洲精品第一综合不卡| xxx96com| 男女之事视频高清在线观看| 999精品在线视频| 99久久综合精品五月天人人| 一二三四社区在线视频社区8| 一本综合久久免费| 香蕉国产在线看| 美女大奶头视频| 午夜免费鲁丝| 大型av网站在线播放| 国产高清视频在线播放一区| 91麻豆av在线| 亚洲精品国产精品久久久不卡| 亚洲午夜精品一区,二区,三区| 在线观看日韩欧美| 亚洲片人在线观看| 精品久久久久久久人妻蜜臀av| 色播亚洲综合网| 亚洲 国产 在线| 亚洲第一电影网av| 国产一区二区三区在线臀色熟女| 日本成人三级电影网站| 桃色一区二区三区在线观看| 88av欧美| 岛国在线观看网站| 国产国语露脸激情在线看| 国产主播在线观看一区二区| 亚洲国产欧美网| 18禁裸乳无遮挡免费网站照片 | 久久性视频一级片| 亚洲成人精品中文字幕电影| 啦啦啦观看免费观看视频高清| 久久精品国产亚洲av高清一级| 国产高清视频在线播放一区| 国产精品 欧美亚洲| 日韩视频一区二区在线观看| 亚洲成人免费电影在线观看| www国产在线视频色| 深夜精品福利| 日韩精品中文字幕看吧| 人人妻人人澡人人看| 一进一出好大好爽视频| 国产伦一二天堂av在线观看| 99在线视频只有这里精品首页| 日韩精品免费视频一区二区三区| 黄色毛片三级朝国网站| 欧美成人一区二区免费高清观看 | 每晚都被弄得嗷嗷叫到高潮| 午夜福利成人在线免费观看| 亚洲无线在线观看| 日本黄色视频三级网站网址| 少妇的丰满在线观看| 天堂动漫精品| 精品福利观看| 夜夜夜夜夜久久久久| 亚洲熟妇中文字幕五十中出| 老司机在亚洲福利影院| 中文字幕精品亚洲无线码一区 | 欧美成人免费av一区二区三区| 国产三级黄色录像| 麻豆成人av在线观看| 亚洲欧美激情综合另类| 国产高清激情床上av| netflix在线观看网站| 国产亚洲av嫩草精品影院| 午夜影院日韩av| 啦啦啦韩国在线观看视频| 久久天躁狠狠躁夜夜2o2o| 亚洲自拍偷在线| 免费观看精品视频网站| 最近最新免费中文字幕在线| av中文乱码字幕在线| 中文字幕最新亚洲高清| 久久人妻福利社区极品人妻图片| 12—13女人毛片做爰片一| 欧美中文综合在线视频| 久久久水蜜桃国产精品网| 侵犯人妻中文字幕一二三四区| e午夜精品久久久久久久| 国产av不卡久久| 麻豆成人av在线观看| 精品欧美一区二区三区在线| 免费高清视频大片| 欧美精品亚洲一区二区| 亚洲成国产人片在线观看| 19禁男女啪啪无遮挡网站| 少妇 在线观看|