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

    筒倉靜態(tài)儲糧的邊界壓力及倉壁摩擦力試驗(yàn)研究

    2018-08-10 07:08:48李東橋陳家豪段君峰
    農(nóng)業(yè)工程學(xué)報 2018年13期
    關(guān)鍵詞:深度

    韓 陽,李東橋,陳家豪,靜 行,段君峰

    ?

    筒倉靜態(tài)儲糧的邊界壓力及倉壁摩擦力試驗(yàn)研究

    韓 陽,李東橋,陳家豪,靜 行,段君峰

    (河南工業(yè)大學(xué)土木建筑學(xué)院,鄭州 450001)

    為了研究筒倉散裝糧堆的邊界壓力和倉壁摩擦力的分布規(guī)律,研制了模型筒倉試驗(yàn)裝置,基于倉體的微縫分離設(shè)計(jì),實(shí)現(xiàn)各分離倉體受力的獨(dú)立測量。以小麥為例,通過實(shí)測,發(fā)現(xiàn)不同裝糧高度下,糧堆底部壓力沿徑向呈現(xiàn)不均勻分布特征,其不均勻分布程度隨裝糧高度逐漸增加;當(dāng)裝糧高度大于筒倉直徑后,倉壁側(cè)壓力開始逐漸小于Janssen公式計(jì)算結(jié)果;而倉壁摩擦力在整個糧堆深度范圍內(nèi)均小于Janssen公式計(jì)算結(jié)果。試驗(yàn)表明,倉壁實(shí)測摩擦力與側(cè)壓力之比小于小麥與倉壁的摩擦系數(shù),且隨糧堆深度的增加不斷變化,表明靜態(tài)儲糧下儲料與倉壁邊界之間尚未達(dá)到極限平衡狀態(tài);側(cè)壓力系數(shù)接近主動態(tài),且小于主動土壓力系數(shù)。研究結(jié)果可為散體物料壓力理論提供參考。

    筒倉;試驗(yàn);壓力;側(cè)壓力;豎向壓力;摩擦力

    0 引 言

    糧食等散體物料具有復(fù)雜的力學(xué)特性。由于可以承受靜態(tài)剪切力,導(dǎo)致其空間靜應(yīng)力σ不能簡單地通過求解力平衡方程?=ρg(假定密度恒定,g為重力加速度)得到,還需要確定物料自身的本構(gòu)關(guān)系和問題的邊界條件。目前,在筒倉糧食儲料壓力理論分析中,小麥等糧食物料還沒有適宜的本構(gòu)模型,而儲料壓力問題的邊界條件也僅是簡單定義倉壁邊界和糧堆自由面兩處,限制了儲料壓力問題的求解精度和適用范圍。因此,僅以目前的條件難以推導(dǎo)普遍適用的糧食空間壓力的計(jì)算方法。而試驗(yàn)研究又受到填料方式、倉體變形、填料時內(nèi)置傳感器偏移、糧堆內(nèi)部應(yīng)力及倉壁摩擦力難以觀測等問題的困擾,有效且可重復(fù)的試驗(yàn)方案還很少見。

    目前,許多研究都是在德國學(xué)者Janssen[1]提出的連續(xù)介質(zhì)模型的基礎(chǔ)上進(jìn)行的,這個模型主要基于3個假設(shè):1)儲料的水平壓力σ與豎向壓力σ的比值(側(cè)壓力系數(shù))保持不變,這個是對筒倉內(nèi)散料應(yīng)力狀態(tài)的猜想;2)倉壁摩擦力τ與倉壁側(cè)壓力σ的比值恒等于儲料與倉壁的摩擦系數(shù)0,即假定任意深度處儲料與倉壁間均達(dá)到極限平衡狀態(tài),這個可以理解為對倉壁邊界條件的猜想。這里,和0的取值受儲料自身的物理參數(shù)、倉壁材質(zhì)及儲料充填過程等因素影響,而非由筒倉靜態(tài)系統(tǒng)所確定的常數(shù)。如側(cè)壓力系數(shù),Janssen公式中將其采用為朗肯[2]主動土壓力系數(shù),即與內(nèi)摩擦角間存在=tan2(45°?/2)的關(guān)系。但目前的一些研究表明,在某些試驗(yàn)條件[3]下甚至?xí)霈F(xiàn)大于1這種偏離朗肯主動態(tài)的現(xiàn)象。而摩擦系數(shù)0則被定義為與儲料和倉壁材質(zhì)有關(guān)的常數(shù),目前試驗(yàn)[4-7]測得的0都是散體物料與固體交界處的靜摩擦系數(shù)(有些是滑動摩擦系數(shù)),即假定儲料與倉壁間存在相對滑移的趨勢。試驗(yàn)研究[8]發(fā)現(xiàn),靜態(tài)儲糧下儲料與倉壁間并未達(dá)到極限平衡狀態(tài),實(shí)測的倉壁摩擦力τ與側(cè)壓力σ之比在整個糧堆范圍內(nèi)均小于0,并且這個比值τ/σ沿糧堆深度不斷變化。本研究將這種顆粒物質(zhì)與其容器邊界之間摩擦不充分的現(xiàn)象稱為“倉壁的摩擦效應(yīng)”,即倉壁摩擦力τ與側(cè)壓力σ的比值τ/σ是一個與糧堆深度有關(guān)的函數(shù),而不是一個與無關(guān)的常數(shù)。

    Janssen的第三點(diǎn)假設(shè)為糧堆任意水平面上豎向壓力σ均勻分布,它的大小僅與有關(guān)。由于糧堆的內(nèi)部應(yīng)力難以觀測,對于這個假設(shè)的合理性尚沒有直接的試驗(yàn)驗(yàn)證。目前關(guān)于Janssen模型的研究[9-10]通常是間接測量糧堆底部邊界處的儲料壓力,進(jìn)而推測豎向壓力σ與糧堆深度的關(guān)系?,F(xiàn)階段試驗(yàn)測量的多為堆底的平均應(yīng)力,由于受難以固定傳感器位置、填料時傳感器受擾動等因素的影響,底部壓力沿徑向的分布情況尚缺乏行之有效的試驗(yàn)方法。

    為更準(zhǔn)確研究筒倉散裝糧堆邊界壓力及倉壁摩擦力的分布規(guī)律,驗(yàn)證Janssen模型,本文設(shè)計(jì)了新型的筒倉模型試驗(yàn)倉,對不同裝糧高度下糧堆底部的豎向壓力沿徑向分布情況、倉壁側(cè)壓力和摩擦力與糧堆深度的關(guān)系進(jìn)行了試驗(yàn)研究,為散體物料壓力理論和倉儲結(jié)構(gòu)研究提供有效的試驗(yàn)方法。

    1 筒倉模型倉試驗(yàn)

    1.1 模型倉

    模型倉的設(shè)計(jì)思想主要基于倉體的微縫分離設(shè)計(jì)和各分離倉體受力的獨(dú)立測量。倉體主要由3段相互分離的筒壁以及1個可自由升降的倉底板組成。倉底板與筒壁之間以及各部分筒壁之間都設(shè)有微縫,這樣既能夠使它們相互分離,又不會使糧食顆粒流出,如圖1所示。

    注:黑色圓形為力傳感器安放位置。尺寸單位為mm。

    模型倉裝糧高徑比可為2~4,直徑約為一般糧食筒倉的1/20,可以對深倉和淺倉狀態(tài)下儲料壓力開展試驗(yàn)研究。

    筒壁用有機(jī)玻璃制成,總高度2 102 mm,內(nèi)徑480 mm,壁厚10 mm,標(biāo)記有刻度,便于觀察和測量裝糧高度。上筒高1 000 mm,固定在鋼架上,使上筒壁受到的摩擦力不向下傳遞;中筒高100 mm,由沿環(huán)向均布的3個力傳感器支撐,測量中筒受到的摩擦力,并在內(nèi)壁沿環(huán)向鑲嵌3個土壓力盒,測量中筒所受側(cè)壓力;下筒高1 000 mm,由沿環(huán)向均布的3個力傳感器支撐,測量下筒壁受到的摩擦力。

    倉底板用10 mm厚鋼板制作,直徑478 mm,分為4個同心圓盤(圖1c),各圓盤間留有小于1 mm的縫隙,分別由力傳感器支承(中心環(huán)1個,其他環(huán)各3個),可分別測量倉底不同部位承受的糧食壓力,進(jìn)而可分析倉底壓力沿徑向的分布規(guī)律;倉底板可升降,最大行程為1 000 mm。通過控制底板的升降,可使中筒位于糧堆的不同高度,從而可測量物料對倉壁任意高程的側(cè)壓力和摩擦力。

    試驗(yàn)選用力傳感器均為福建銳馬電器制造有限公司生產(chǎn)的RM-S200型,準(zhǔn)確度等級C3;土壓力盒為溧陽超源儀器廠生產(chǎn)的BW型,準(zhǔn)確度級別0.2%。編制了數(shù)據(jù)自動采集和分析系統(tǒng),可采集倉底板各圓盤所受豎向壓力、下筒所受摩擦力和中筒側(cè)壓力、摩擦力,系統(tǒng)界面如圖1d所示。

    1.2 試驗(yàn)方案

    Janssen猜想給出的一個重要結(jié)論為:

    式中為儲料重力密度(N/m3),為筒倉半徑(m)。根據(jù)Janssen的假設(shè),糧堆深度(m)處儲料豎向壓力σ(kPa)均勻分布,且側(cè)壓力σ(kPa)與豎向壓力σ(kPa)的比值恒定為側(cè)壓力系數(shù),倉壁摩擦力(kPa)等于靜摩擦系數(shù)0與側(cè)壓力σ(kPa)的乘積。

    為驗(yàn)證上述Janssen猜想及假設(shè),結(jié)合模型倉試驗(yàn)裝置的特點(diǎn)及功能,制定了以下試驗(yàn)方案:

    1)測試倉底板糧食豎向壓力σ分布和倉壁摩擦力:倉底板置于模型倉底部固定不動,每次裝糧約100 mm高,并對進(jìn)糧質(zhì)量精確稱量,裝糧至1 000 mm高。利用倉底板各自圓盤下設(shè)置的力傳感器測試倉底板不同部位糧食壓力,得到不同裝糧高度下倉底壓力σ的變化和沿徑向的分布情況;利用支撐下筒壁的力傳感器測量不同裝糧高度下倉壁摩擦力的變化趨勢。

    2)測試糧食不同深度處對倉壁的側(cè)壓力σ和摩擦力:倉底板初始位置置于下筒頂部,一次性裝糧1 000 mm高并稱總量,倉底板每次下降50 mm,靜置一段時間至各傳感器數(shù)值穩(wěn)定后,利用鑲嵌在中筒內(nèi)壁的土壓力盒和支撐中筒的力傳感器分別測量糧堆不同深度處儲料對倉壁的側(cè)壓力σ和摩擦力。

    1.3 小麥試樣

    儲料選用河南滎陽小麥,根據(jù)《糧油檢驗(yàn)容重測定》GB/T 5498-2013[11]、三軸壓縮試驗(yàn)[12]和直剪試驗(yàn)[8],測得小麥試樣參數(shù)如表1所示。

    表1 小麥試樣參數(shù)Table 1 Wheat sample parameters

    注:外摩擦角為小麥與模型倉倉壁材質(zhì)(有機(jī)玻璃)間的外摩擦角。

    Note: External friction angle is the external friction angle between the wheat and the material (PMMA) of the silo wall.

    1.4 試驗(yàn)數(shù)據(jù)處理方法

    根據(jù)方案1)進(jìn)行重復(fù)試驗(yàn)4次,方案2)重復(fù)試驗(yàn)3次,數(shù)據(jù)均取各傳感器均值,使用Excel2003進(jìn)行數(shù)據(jù)處理及圖形繪制;Minitab進(jìn)行誤差分析及數(shù)值擬合,其中誤差棒為相同裝糧高度(或糧堆深度)下,各次試驗(yàn)均值±標(biāo)準(zhǔn)差。

    2 試驗(yàn)結(jié)果

    2.1 糧堆底部豎向壓力

    倉內(nèi)糧食總重力應(yīng)等于模型倉底板各圓盤承受壓力總值與筒壁所受總摩擦力值之和,可以此作為試驗(yàn)裝置設(shè)計(jì)是否完善、傳感器是否準(zhǔn)確的依據(jù)。方案1)試驗(yàn)中裝糧任意高程下糧食重力與力傳感器總值間誤差皆小于5%,為可接受范圍內(nèi),一定程度上證明了數(shù)據(jù)的準(zhǔn)確,表明試驗(yàn)有較好的可重復(fù)性。

    取各次試驗(yàn)結(jié)果平均值,得到了不同裝糧高度下倉底壓力分布圖,如圖2所示。

    圖2 不同裝糧高度下糧堆底部豎向壓力

    由圖2可以看出,糧堆底部的儲料豎向壓力呈現(xiàn)出不均勻分布特征。隨著裝糧高度的增加,靠近倉壁處壓力與倉中心處的壓力差距逐漸明顯。在小麥堆載高度小于0.4 m,即高徑比小于1.0時,倉底各圓環(huán)壓力值與小麥堆載高度近乎呈線性增加關(guān)系;當(dāng)糧堆高度超過0.5 m后,隨著小麥堆載高度的增加,各圓環(huán)壓力值繼續(xù)增大,但增加速度減慢。倉底壓力曲線拐點(diǎn)(0.5 m)的出現(xiàn),反映了淺倉與深倉倉底壓力的區(qū)別。

    2.2 倉壁側(cè)壓力

    根據(jù)方案2)進(jìn)行試驗(yàn),記錄了倉壁側(cè)壓力與糧堆深度之間的關(guān)系,并與一些經(jīng)典公式[13]的結(jié)果進(jìn)行對比,如圖3所示。

    圖3 倉壁側(cè)壓力與糧堆深度的關(guān)系

    由于Reimbert公式是半經(jīng)驗(yàn)公式,部分內(nèi)容為試驗(yàn)曲線擬合得到,且其試驗(yàn)選取物料為砂,若不經(jīng)任何修改,它的適用范圍有限,故圖3中,其結(jié)果與試驗(yàn)結(jié)果有較大偏差。Janssen公式結(jié)果在糧堆深度較小時與試驗(yàn)結(jié)果幾乎一致;在糧堆深度大于筒倉直徑(0.5 m)時則略小于試驗(yàn)結(jié)果。這個結(jié)果與一些文獻(xiàn)[14-16]報道的相仿,即Janssen公式中得到的側(cè)壓力比實(shí)際值偏小,所以許多國家的規(guī)范[17-18]都在計(jì)算中乘以水平壓力修正系數(shù)C以保證結(jié)構(gòu)的安全。

    2.3 倉壁摩擦力

    倉壁摩擦力沿倉壁的分布可以直接由隔離的中筒測量得到,也可以由每次底板下降過程中由下筒壁測得的摩擦力的增量得到,如圖4a所示。從圖4a中可見二者基本一致,表明測試儀器工作正常。試驗(yàn)得到倉壁摩擦力與糧堆深度的關(guān)系如圖4b所示。可以看到,試驗(yàn)結(jié)果與Janssen公式和Reimbert公式計(jì)算結(jié)果都有一定偏差。Janssen公式得到的摩擦力在整個糧堆深度范圍內(nèi)均大于試驗(yàn)結(jié)果。當(dāng)糧堆深度小于0.6 m時,Reimbert公式計(jì)算得到的摩擦力小于試驗(yàn)結(jié)果;這之后,則大于試驗(yàn)結(jié)果。

    圖4 倉壁摩擦力與糧堆深度的關(guān)系

    3 分析與討論

    3.1 堆底豎向壓力分析

    筒倉儲料的靜態(tài)壓力有一種被稱為“糧倉效應(yīng)”[19-20]的奇異現(xiàn)象,即當(dāng)糧堆達(dá)到某一高度后,堆底壓力不再隨裝糧高度的增加而增加,而是趨于飽和。造成這種現(xiàn)象的原因?yàn)閮α蟽?nèi)部的剪切力、儲料與倉壁間的摩擦力等在筒倉內(nèi)形成拱效應(yīng),使儲料的一部分重力轉(zhuǎn)移到倉壁上。圖5a和圖5b分別繪制了倉底平均壓力與裝糧高度的關(guān)系及倉底不同區(qū)域豎向壓力與裝糧高度的關(guān)系。

    圖5 倉底壓力與裝糧高度的關(guān)系

    從圖5a中可以看到,在裝糧高度較低時,糧堆底部的平均壓力接近線性增長;當(dāng)糧堆高度達(dá)到倉徑(0.5 m)后,增加趨勢逐漸變緩。圖5b中,倉底不同區(qū)域處儲料壓力也呈現(xiàn)出不同程度的“糧倉效應(yīng)”。當(dāng)裝糧高度小于0.2 m時,各圓環(huán)區(qū)域的壓力曲線近乎重合;裝糧高度大于0.3 m后,各環(huán)壓力值開始出現(xiàn)區(qū)別,相同裝糧高度下,越接近倉壁邊界的圓環(huán)所承受的儲料壓力越小。三環(huán)壓力(即倉壁邊界處儲料豎向壓力)曲線與其余各環(huán)曲線區(qū)別較為明顯,它很早地結(jié)束了線性變化階段,進(jìn)入趨于飽和的非線性變化階段。而其余各環(huán)的變化趨勢與它們距倉壁的距離有關(guān):離倉壁越遠(yuǎn),越晚結(jié)束線性變化。這點(diǎn)體現(xiàn)出倉壁摩擦的影響,即相同裝糧高度下,倉壁摩擦對儲料的影響是逐漸減弱的,距倉壁越遠(yuǎn),儲料間豎向切應(yīng)力越小,這是造成糧堆底部豎向壓力不均勻分布的主要原因。

    3.2 倉壁摩擦效應(yīng)

    Janssen模型中假定倉壁摩擦力τ與側(cè)壓力σ的比值0恒定,其模型中隔離的儲料薄片的豎向平衡方程中以0·σ表示摩擦力τ。本試驗(yàn)中測得的倉壁摩擦力τ明顯小于Janssen公式中0·σ的結(jié)果,且τ/σ沿糧堆深度發(fā)生變化。本文以μ表示試驗(yàn)測得的摩擦力與側(cè)壓力的比值τ/σ,對圖表進(jìn)行無量綱處理,以糧堆深度與模型倉內(nèi)徑的比值/作為軸,μ與/的關(guān)系如圖6所示。

    注:μwg表示試驗(yàn)測得的摩擦力與側(cè)壓力的比值τw/σr,z/D表示深度比,下同。

    在圖6中,當(dāng)/小于0.3時,μ隨著/的增加而線性增大;當(dāng)/大于0.3而小于0.8時,μ非線性增大;其后到糧堆底面,μ又近似線性減小。μ在整個糧堆深度范圍內(nèi)均小于小麥與模型倉倉壁材料間的靜摩擦系數(shù)0。這表明靜態(tài)儲糧下,儲料與倉壁邊界間并未達(dá)到極限平衡狀態(tài),即儲料與倉壁間的摩擦并沒有充分“動員起來”,所以Janssen模型中假定的倉壁邊界條件(τ/σ恒定為靜摩擦系數(shù))與實(shí)際情況有一定偏離。本文將這種散體物料與其容器邊界摩擦不充分的現(xiàn)象稱為“倉壁的摩擦效應(yīng)”。有些研究[21-22]通過“上拉倉壁”或“倉底板下沉”等方法來使試驗(yàn)倉的邊界條件更接近Janssen模型,試驗(yàn)得到的結(jié)果也與Janssen吻合較好。但實(shí)際工程中,大多的糧倉幾乎都是直接填料的,很難進(jìn)行這樣的后續(xù)處理。所以如果能確定靜態(tài)儲糧下τ/σ的值(或表達(dá)式),并以此改進(jìn)Janssen的倉壁邊界條件,會為倉儲研究提供更準(zhǔn)確的儲料荷載。

    3.3 側(cè)壓力系數(shù)

    根據(jù)不同國家規(guī)范及學(xué)者們的研究[23-25],側(cè)壓力系數(shù)有許多取值方法。如中國規(guī)范將其選取為主動土壓力系數(shù)、美國規(guī)范采用靜止土壓力系數(shù),而歐洲規(guī)范[26]則對靜止土壓力系數(shù)進(jìn)行了修正。由于幾種方法對值的計(jì)算結(jié)果偏差較大,學(xué)者們對其的選取問題仍有爭議[27-30]。受測量方法的限制,目前的試驗(yàn)多以倉壁靜摩擦系數(shù)0來估算值的下限,因?yàn)槎呖偸邱詈显谝黄鸪霈F(xiàn),以此驗(yàn)證Janssen猜想。但是如3.2節(jié)所述,顆粒材料與固體邊界間摩擦力的分布規(guī)律本身就是一個十分復(fù)雜的問題,所以對側(cè)壓力系數(shù)的選取問題仍需多方面的驗(yàn)證。根據(jù)模型倉試驗(yàn)數(shù)據(jù),本文略去邊界效應(yīng)影響,以裝糧高度為的倉底壓力近似表示糧堆深度處的豎向壓力,以儲料在倉壁處的側(cè)壓力與同一水平面的豎向壓力之比作為倉壁處側(cè)壓力系數(shù)值,如圖7所示。

    圖7 側(cè)壓力系數(shù)K與z/D的關(guān)系

    在圖7中,當(dāng)/小于1.0時,隨著/的增加而非線性增大;當(dāng)/大于1.0后,值增速變緩,在糧堆底面,值又略有減小。在整個糧堆深度范圍內(nèi),試驗(yàn)所得的側(cè)壓力系數(shù)接近主動態(tài),但小于主動土壓力系數(shù)tan2(45°-/2),這點(diǎn)同樣與Janssen猜想有所偏離。

    3.4 對修正Janssen方程的建議

    糧倉系統(tǒng)的倉壁摩擦及內(nèi)部顆粒物質(zhì)的應(yīng)力狀態(tài)是兩個十分復(fù)雜的問題,因?yàn)樵囼?yàn)偶然誤差、儲料密度隨裝糧高度變化等問題的存在,目前還很難解釋它們的規(guī)律。就我們所做的有限元模擬[8]來說,兩者的分布應(yīng)與顆粒尺寸、筒倉直徑、裝糧高度及深度比等因素有關(guān)。

    以試驗(yàn)測得的μ與側(cè)壓力系數(shù)相乘,得到μK與深度比/的關(guān)系,如圖8a所示。圖8a中,在糧堆深度比/小于1時,μK非線性增大,之后增加趨勢趨于平緩;在/大于1.6后,又有所減小。雖然μK并非恒定,但可以看到在深度比大于1后,它開始圍繞某一點(diǎn)波動,其上限為0.17,下限為0.15。取深度比1~2之間的試驗(yàn)結(jié)果進(jìn)行數(shù)值擬合,發(fā)現(xiàn)本文模型倉試驗(yàn)的μK在0.159附近漲落。將μK擬合值代入Janssen模型,并將其結(jié)果與試驗(yàn)結(jié)果及Janssen公式計(jì)算結(jié)果對比,如圖8b所示??梢钥吹剑瑪M合結(jié)果與試驗(yàn)數(shù)據(jù)較為接近,Janssen公式計(jì)算結(jié)果結(jié)果則明顯小于試驗(yàn)數(shù)據(jù)。結(jié)合3.2節(jié)及3.3節(jié)試驗(yàn)結(jié)果可知,Janssen對0和的選取有所偏大。

    目前人們對Janssen模型的修正更多考慮的是如何選取適宜的壓力修正系數(shù)或值的放大系數(shù),通過增大計(jì)算值來保證結(jié)構(gòu)的安全,而對μ和的分布規(guī)律的研究還比較少見。如果能進(jìn)一步精確兩者的分布函數(shù)或不同協(xié)議下的近似值,以此改進(jìn)Janssen模型,也是一種較為可行的方法。

    圖8 試驗(yàn)數(shù)據(jù)與擬合結(jié)果的對比

    4 結(jié) 論

    通過自行研制的模型筒倉試驗(yàn)裝置,以小麥為例,對Janssen猜想及Janssen模型的邊界條件進(jìn)行了驗(yàn)證,討論了糧堆底部豎向壓力、倉壁側(cè)壓力和摩擦力的分布規(guī)律,得出以下結(jié)論:

    1)試驗(yàn)表明,倉底的豎向壓力是不均勻分布的,并隨著裝糧高度的增加,不均勻分布的程度逐漸明顯,靠近倉壁處壓力與倉中心處的壓力差距逐漸增大;倉底邊界沿徑向不同區(qū)域處的儲料壓力隨糧堆高度的增加,均會出現(xiàn)一定程度的趨于飽和的現(xiàn)象,其變化趨勢與觀測點(diǎn)至倉壁距離有關(guān)。分析表明,倉壁摩擦力及儲料間的切應(yīng)力的共同作用是造成這些現(xiàn)象的主要因素。

    2)試驗(yàn)表明,小麥深度小于筒倉直徑時,倉壁側(cè)壓力與小麥堆深度呈線性增加關(guān)系,與Janssen公式計(jì)算值基本一致;當(dāng)深度超過筒倉直徑后,側(cè)壓力值與小麥堆深度呈非線性增大關(guān)系,增加速度減緩,表明糧堆深度與倉徑比對倉壁側(cè)壓力的影響,試驗(yàn)數(shù)據(jù)略大于Janssen公式計(jì)算值。

    3)試驗(yàn)所得的倉壁摩擦力小于Janssen猜想提出的的靜摩擦系數(shù)與側(cè)壓力的乘積0·σ,并且倉壁摩擦力與倉壁側(cè)壓力的比值μ沿糧堆深度變化,在糧堆深度范圍內(nèi)均小于靜摩擦系數(shù)0,說明靜態(tài)儲糧下儲料與倉壁邊界之間尚未達(dá)到極限平衡狀態(tài);試驗(yàn)所得的側(cè)壓力系數(shù)與μ的變化趨勢相仿,同樣在糧堆深度范圍內(nèi)不斷變化,且小于主動土壓力系數(shù)tan2(45°?/2)。本文試驗(yàn)中它們的乘積μK在深度比1~2之間圍繞0.159波動,Janssen模型對二者的取值均有所高估,擬合結(jié)果表明進(jìn)一步精確μ和的選取問題,是一種修正Janssen方程的可行方法。

    試驗(yàn)結(jié)果表明,該模型試驗(yàn)倉運(yùn)行方便,每次試驗(yàn)糧食質(zhì)量與傳感器測試誤差皆在5%以內(nèi),試驗(yàn)可重復(fù)性較好。該模型倉和相關(guān)試驗(yàn)方案可為散體物料壓力理論的試驗(yàn)研究提供簡便可行的技術(shù)支持。

    [1] Janssen HA. Experiments about pressures of grain in silos[J]. Z. des Vereines Deutscher Ingenieure, 1895, 39(35): 1045-1049.

    [2] Rankine W J M. On the stability of loose earth[J]. Proceedings of the Royal Society, 1857, 147: 9-27.

    [3] 李湘群,蔣亦民,彭政,等. Rankine被動應(yīng)力狀態(tài)糧倉的Janssen行為[J]. 山東大學(xué)學(xué)報:理學(xué)版,2010,45(9):101-104.

    Li Xiangqun, Jiang Yimin, Peng Zheng, et al. Janssen behavior of silos with Rankine passive stress[J]. Journal of Shandong University: Natural Science, 2010, 45(9): 101-104. (in Chinese with English abstract)

    [4] 安蓉蓉. 糧食的內(nèi)摩擦角、彈性模量及體變模量的試驗(yàn)研究[D]. 南京:南京財經(jīng)大學(xué),2010.

    An Rongrong. Experiment Study on Angle of Internal Friction, Elastic Modulus and Bulk Strain Modulus of Grain[D]. Nanjing: Nanjing University of Finance and Economics, 2010. (in Chinese with English abstract)

    [5] 彭政,王璐珠,蔣亦民. 顆粒物質(zhì)與固體交界面靜摩擦系數(shù)的測量與分析[J]. 山東大學(xué)學(xué)報:理學(xué)版,2011,46(1):42-45.

    Peng Zheng, Wang Luzhu, Jiang Yimin. Measurement and analysis of static friction coefficient on a granular-solid interface[J]. Journal of Shandong University: Natural Science, 2011, 46(1): 42-45. (in Chinese with English abstract)

    [6] 許啟鏗,金立兵,王錄民,等. 糧食力學(xué)參數(shù)的試驗(yàn)研究[J]. 河南工業(yè)大學(xué)學(xué)報:自然科學(xué)版,2010,31(1):18-21.

    Xu Qikeng, Jin Libing, Wang Lumin, et al. Research on the determination of mechanical parameter of grains[J]. Journal of Henan University of Technology: Natural Science Edition, 2010, 31(1): 18-21. (in Chinese with English abstract)

    [7] 程緒鐸,陸琳琳,石翠霞. 小麥摩擦特性的試驗(yàn)研究[J]. 中國糧油學(xué)報,2012,27(4):15-19.

    Cheng Xuduo, Lu Linlin, Shi Cuixia. The experimental research on friction properties of wheat[J]. Journal of the Chinese Cereals and Oils Association, 2012, 27(4): 15-19. (in Chinese with English abstract)

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

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

    [9] 蔣亦民,鄭鶴鵬. 一種顆粒底部壓力不趨向飽和的糧倉系統(tǒng)[J]. 物理學(xué)報,2008, 57(11):7360-7366.

    Jiang Yimin, Zheng Hepeng. A silo of which the pressure on the bottom of granular matter does not go to saturation[J]. Acta Physica Sinica, 2008, 57(11): 7360-7366. (in Chinese with English abstract)

    [10] 李智峰,彭政,蔣亦民. 糧倉內(nèi)顆粒壓力的測量: Janssen行為及其偏差[J]. 物理學(xué)報,2014,63(10):296-303.

    Li Zhifeng, Peng Zheng, Jiang Yimin. Measurements of granular pressure in silo: Janssenbehaviour and deviation[J]. Acta Physica Sinica, 2014, 63(10): 296-303. (in Chinese with English abstract)

    [11] GB/T 5498-2013. 糧油檢驗(yàn)容重測定[S].

    [12] 許啟鏗,陳家豪,王錄民. 小麥力學(xué)參數(shù)的三軸壓縮試驗(yàn)研究[J]. 河南工業(yè)大學(xué)學(xué)報:自然科學(xué)版,2015,36(5):101-105.

    Xu Qikeng, Chen Jiahao, Wang Lumin. Mechanical parameters of wheat in triaxial compression tests[J]. Journal of Henan University of Technology: Natural Science Edition, 2015, 36(5): 101-105. (in Chinese with English abstract)

    [13] Reimbert M L, Reimbert A M. Silos: Theory and Practice[M]. Clausthal: Trans Tech Publictions, 1976.

    [14] 陳長冰. 筒倉內(nèi)散體側(cè)壓力沿倉壁分布研究[D]. 合肥:合肥工業(yè)大學(xué),2006.

    Chen Changbing. Research on Distribution of Lateral Bulk-Solid Pressures on Silos’ walls[D]. Hefei: Hefei University of Technology, 2006. (in Chinese with English abstract)

    [15] 劉定華,王建華,楊建斌. 鋼筋混凝土筒倉側(cè)壓力的試驗(yàn)研究[J]. 西安建筑科技大學(xué)學(xué)報,1995(1):8-12.

    Liu Dinghua, Wang Jianhua, Yang Jianbin. An experimental study of the lateral pressure on the wall of reinforced concrete silo[J]. Journal of Xian University of Architecture and Technology, 1995(1): 8-12. (in Chinese with English abstract)

    [16] 劉定華,郝際平. 鋼筋混凝土筒倉倉壁側(cè)壓力的研究[J]. 建筑結(jié)構(gòu)學(xué)報,1995(5):57-63.

    Liu Dinghua, Hao Jiping. Research on lateral pressure upon wall of reinforced concrete silo[J]. Journal of Building Structures, 1995(5): 57-63. (in Chinese with English abstract)

    [17] GB50077-2003. 鋼筋混凝土筒倉設(shè)計(jì)規(guī)范[S].

    [18] ACI Committee 313-97. Standard practice for design and construction of concrete silos and stacking tubes for storing granular materials[S].

    [19] 陸坤權(quán),劉寄星. 顆粒物質(zhì)(上)[J]. 物理,2004,33(9):629-635. Lu Kunquan, Liu Jixing. Static and dynamic properties of granular matter(I)[J]. Physics, 2004, 33(9): 629-635. (in Chinese with English abstract)

    [20] 王晶. 不同維度下顆粒鏈的糧倉效應(yīng)[D]. 北京:北京理工大學(xué),2016.

    Wang Jing. Bottom Stress of Static Packing of Granular Chains in 2D and 3D Silos[D]. Beijing: Beijing Institute of Technology, 2016. (in Chinese with English abstract)

    [21] Vanel L, Clément E. Pressure screening and fluctuations at the bottom of a granular column[J]. The European Physical Journal B-Condensed Matter and Complex Systems, 1999, 11(3): 525-533.

    [22] Vanel L, Claudin P, Bouchaud J P, et al. Stresses in silos: Comparison between theoretical models and new experiments[J]. Physical Review Letters, 1999, 84(7): 1439-1442.

    [23] Jaky J. The coefficient of earth pressure at rest[J]. Journal for Society of Hungarian Architects and Engineers, 1944(5): 355-358.

    [24] 張家康,黃文萃,姜濤,等. 筒倉貯料側(cè)壓力系數(shù)研究[J]. 建筑結(jié)構(gòu)學(xué)報,1999,20(1):71-74.

    Zhang Jiakang, Huang Wencui, Jiang Tao, et al. Study on lateral pressure ratio of stored material in silo[J]. Journal of Building Structures, 1999, 20(1): 71-74. (in Chinese with English abstract)

    [25] 曾丁,黃文彬. 筒倉壁壓的有限元分析[J]. 農(nóng)業(yè)工程學(xué)報,1998,14(2):44-48.

    Zeng Ding. Huang Wenbin. Finite element analysis of silo wall pressure[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1998, 14(2): 44-48. (in Chinese with English abstract)

    [26] BS EN1991-4:2006. Eurocode 1-Actions on structures-Part 4: Silos and tanks[S].

    [27] 馮云田,華云龍. 適用于柔性圓筒倉的修正的Janssen公式[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報,1996(4):107-111.

    Feng Yuntian, Hua Yunlong. Modified Janssen formula for flexible circular bins[J]. Journal of China Agricultural University, 1996(4): 107-111. (in Chinese with English abstract)

    [28] 曾丁,郝保紅,黃文彬. 筒倉靜態(tài)壁壓的有限元分析[J]. 中國粉體技術(shù),2000(5):6-10.

    Zeng Ding, Hao Baohong, Huang Wenbin. Finite element analysis of static wall pressure of silo[J]. China Powder Science and Technology, 2000(5): 6-10. (in Chinese with English abstract)

    [29] 陳長冰. 考慮倉壁彈性的糧食鋼筒倉側(cè)壓力系數(shù)修正[J]. 中國農(nóng)機(jī)化學(xué)報,2006(3):82-83.

    Chen Changbing. Lateral pressure ratio of food steel silos based on considering the wall elasticity[J]. Journal of Chinese Agricultural Mechanization, 2006(3): 82-83. (in Chinese with English abstract)

    [30] 孫珊珊,趙均海,張常光,等. 深倉和淺倉儲料側(cè)壓力新解[J]. 廣西大學(xué)學(xué)報:自然科學(xué)版,2018,43(1):168-177.

    Sun Shanshan, Zhao Junhai, Zhang Changguang, et al. New solution for lateral pressure of silos and bunkers[J]. Journal of Guangxi University: Natural Science Edition, 2018, 43(1): 168-177. (in Chinese with English abstract)

    Experimental study on boundary pressure and wall friction under static grain storage in silo

    Han Yang, Li Dongqiao, Chen Jiahao, Jing Hang, Duan Junfeng

    (,450001,)

    The authors participated in the design of the silo model test facility, which can be used to measure the storage force of the grain heap, such as the bottom pressure, the side wall pressure and the friction force of the bulk grain heap. The self-designed silo model was 480 mm in diameter and 2 102 mm high. The silo wall, constructed of organic glass, was composed of 3 separate parts. There is a 1-mm thick gap between each part of the walls to prevent the outflow of grain. The 3 parts of the wall were named Top Wall (l 000 mm high, fixed on the steel frame), Middle Wall (100 mm high) and Upper Wall (1 000 mm high). The silo bottom, which can be raised and lowered, was divided into 4 parts (1 concentric circle and 3 concentric rings). The Top Wall, the Middle Wall and each parts of silo bottom were supported on three force sensors to measure the wall friction as well as the vertical pressure on the floor. Three earth pressure cells were embedded in the middle line of the Middle Wall to measure the wall lateral pressure. Two experimental schemes were carried out. No.1 scheme is that: Fill the silo model with wheat to a height of 100 mm at each time until the height of the wheat reached 1 000 mm. Record the data obtained from each of the force sensors supported under each part of the silo bottom as well as the force sensors supported under the Upper Wall of each grain height. This scheme is designed to measure the radial distribution of bottom pressure under different grain heights. The No.2 scheme is that: Fill the silo model with wheat to a height of 1 000 mm and then raise the flat bottom to the lower edge of the Middle Wall. Drop the bottom by 50 mm at each time. Meanwhile, record the data obtained from the earth pressure cells embedded in the Middle Wall and the data obtained from the force sensors supported under the Middle Wall at difficult depths of wheat. The distribution of the storage force at the boundary was measured in the experiment. The rationality of the model test data is verified by compared with the classic theory. Combined with the experimental data, it is found that the bottom pressure of silo grain heap presents the uneven distribution characteristics of “l(fā)arge in middle and small on sides”, and the uneven distribution degree is gradually obvious with the increase of grain loading height. The lateral pressure of silo wall is slightly larger than that of Janssen formula when the depth is more than 0.5 m. The friction of the side wall is larger than that of Janssen formula along the whole depth of wheat.Andthe value of the friction coefficient between the grain and the silo wall is smaller than external friction coefficient measured by directly sheared tests along the whole depth of wheat. It proved that the value of wall friction calculated by Janssen formula is larger than what it really is and the state between the wheat and silo wall haven’t reach the limit equilibrium. The lateral pressure coefficient measured by experiment is closed to Rankine's main dynamic and is less than active earth pressure coefficient. The experiment results show that the model silo test device has good repeatability, and it can provide a simple and feasible technical support for the theory of storage pressure.

    silo; experiment; pressure; lateral pressure; vertical pressure; friction

    韓 陽,李東橋,陳家豪,靜 行,段君峰. 筒倉靜態(tài)儲糧的邊界壓力及倉壁摩擦力試驗(yàn)研究[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(13):296-302. doi:10.11975/j.issn.1002-6819.2018.13.036 http://www.tcsae.org

    Han Yang, Li Dongqiao, Chen Jiahao, Jing Hang, Duan Junfeng. Experimental study on boundary pressure and wall friction under static grain storage in silo[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(13): 296-302. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.13.036 http://www.tcsae.org

    2018-04-04

    2018-05-10

    國家自然科學(xué)基金(51608176)

    韓 陽,教授,博士生導(dǎo)師。研究方向?yàn)閭}儲結(jié)構(gòu)及防災(zāi)減災(zāi)工程。Email:hanyangh@126.com

    10.11975/j.issn.1002-6819.2018.13.036

    TS210

    A

    1002-6819(2018)-13-0296-07

    猜你喜歡
    深度
    深度理解不等關(guān)系
    四增四減 深度推進(jìn)
    深度理解一元一次方程
    深度觀察
    深度觀察
    深度觀察
    深度觀察
    芻議深度報道的深度與“文”度
    新聞傳播(2016年10期)2016-09-26 12:14:59
    提升深度報道量與質(zhì)
    新聞傳播(2015年10期)2015-07-18 11:05:40
    微小提議 深度思考
    麻豆国产av国片精品| 午夜免费激情av| 久久国内精品自在自线图片| 中文在线观看免费www的网站| 国产av麻豆久久久久久久| 中文亚洲av片在线观看爽| 国产精品无大码| 国产av不卡久久| 69人妻影院| 一区二区三区免费毛片| 欧美国产日韩亚洲一区| 亚洲最大成人中文| 亚洲成人中文字幕在线播放| 俄罗斯特黄特色一大片| 久久韩国三级中文字幕| 亚洲人成网站在线播| 欧美激情国产日韩精品一区| 日本色播在线视频| 日日干狠狠操夜夜爽| 乱系列少妇在线播放| 国产成人a区在线观看| 变态另类成人亚洲欧美熟女| 久久久久久久久大av| 69av精品久久久久久| 啦啦啦啦在线视频资源| 日日摸夜夜添夜夜添av毛片| 久久人人精品亚洲av| 精品久久国产蜜桃| 午夜爱爱视频在线播放| 成人精品一区二区免费| 成人av一区二区三区在线看| 天堂动漫精品| 国产精品一区www在线观看| 久久久久久久亚洲中文字幕| 三级国产精品欧美在线观看| 精品久久久久久久久av| 日本撒尿小便嘘嘘汇集6| 日韩欧美在线乱码| av国产免费在线观看| 91久久精品国产一区二区三区| 午夜免费激情av| 国产大屁股一区二区在线视频| 国产精品一区二区三区四区免费观看 | 成人特级黄色片久久久久久久| 日本五十路高清| 欧美另类亚洲清纯唯美| 国产黄a三级三级三级人| av天堂在线播放| 午夜视频国产福利| 波多野结衣巨乳人妻| 成人欧美大片| 国产伦在线观看视频一区| 国产高清不卡午夜福利| 久久热精品热| 少妇人妻精品综合一区二区 | 日本-黄色视频高清免费观看| 白带黄色成豆腐渣| 国产精品乱码一区二三区的特点| 国产一级毛片七仙女欲春2| 国产精品无大码| 97超视频在线观看视频| 嫩草影院新地址| 精品欧美国产一区二区三| 桃色一区二区三区在线观看| 成年av动漫网址| 日本三级黄在线观看| 亚洲欧美精品综合久久99| 亚洲av电影不卡..在线观看| 国产综合懂色| 老司机福利观看| 国产单亲对白刺激| 一进一出抽搐动态| 成年女人毛片免费观看观看9| 免费看光身美女| 亚洲精品乱码久久久v下载方式| 又爽又黄无遮挡网站| 国产精品三级大全| 亚洲18禁久久av| 97超视频在线观看视频| 搡老妇女老女人老熟妇| 69av精品久久久久久| 少妇的逼好多水| 九九在线视频观看精品| av在线观看视频网站免费| 久久精品人妻少妇| 免费看光身美女| 亚洲无线观看免费| 中文字幕久久专区| 国产在线男女| 午夜精品国产一区二区电影 | 又粗又爽又猛毛片免费看| 亚洲精品亚洲一区二区| 成人国产麻豆网| 亚州av有码| 欧美一区二区亚洲| 高清毛片免费观看视频网站| 麻豆精品久久久久久蜜桃| 别揉我奶头 嗯啊视频| 淫妇啪啪啪对白视频| 国产免费男女视频| 波多野结衣高清无吗| 亚洲最大成人av| 国产真实乱freesex| 99久久成人亚洲精品观看| 久久韩国三级中文字幕| 精品久久国产蜜桃| 中国美女看黄片| 亚洲成人久久爱视频| 特大巨黑吊av在线直播| 色5月婷婷丁香| 免费观看在线日韩| 丰满的人妻完整版| 亚洲高清免费不卡视频| 亚洲精品国产av成人精品 | 国产精品人妻久久久久久| 亚洲最大成人中文| 看非洲黑人一级黄片| 女人十人毛片免费观看3o分钟| 欧美中文日本在线观看视频| 国产精品不卡视频一区二区| 一个人看的www免费观看视频| 午夜影院日韩av| 精品久久国产蜜桃| 亚洲av免费高清在线观看| 久久人人爽人人片av| 亚洲av成人精品一区久久| 亚洲专区国产一区二区| 精品久久久久久久久亚洲| 99九九线精品视频在线观看视频| 免费观看的影片在线观看| 一级黄片播放器| 99久国产av精品| 午夜精品一区二区三区免费看| 久久久久久久久中文| 亚洲av中文字字幕乱码综合| 少妇熟女aⅴ在线视频| 成人av一区二区三区在线看| 免费av毛片视频| 欧美性猛交黑人性爽| 搡老妇女老女人老熟妇| 欧美xxxx黑人xx丫x性爽| 成年免费大片在线观看| 久久久a久久爽久久v久久| 亚洲国产精品国产精品| 天堂网av新在线| 人妻制服诱惑在线中文字幕| 亚洲一区二区三区色噜噜| 欧美bdsm另类| 有码 亚洲区| 午夜福利在线观看吧| 亚洲成a人片在线一区二区| 性欧美人与动物交配| 中出人妻视频一区二区| 十八禁国产超污无遮挡网站| 一个人看的www免费观看视频| 亚洲欧美日韩无卡精品| 少妇裸体淫交视频免费看高清| 性欧美人与动物交配| 小蜜桃在线观看免费完整版高清| 日本五十路高清| 丰满乱子伦码专区| 午夜福利在线观看吧| aaaaa片日本免费| 麻豆乱淫一区二区| 99久国产av精品| 亚洲美女搞黄在线观看 | 悠悠久久av| 国内揄拍国产精品人妻在线| 久久这里只有精品中国| 国产精品,欧美在线| 久久精品综合一区二区三区| 麻豆国产97在线/欧美| 欧美一区二区国产精品久久精品| 国产精品女同一区二区软件| 国产精品一区www在线观看| 久久久精品大字幕| 亚洲专区国产一区二区| 91麻豆精品激情在线观看国产| 我要搜黄色片| 三级毛片av免费| 成人鲁丝片一二三区免费| 亚洲精品456在线播放app| 国内精品久久久久精免费| 欧美最黄视频在线播放免费| 成人亚洲欧美一区二区av| www日本黄色视频网| 久久久久久大精品| 午夜精品在线福利| 欧美极品一区二区三区四区| 精品少妇黑人巨大在线播放 | 在线播放国产精品三级| 国内久久婷婷六月综合欲色啪| 国内精品一区二区在线观看| 深夜精品福利| 蜜桃久久精品国产亚洲av| 一区二区三区高清视频在线| av卡一久久| 内地一区二区视频在线| 国产成人freesex在线 | 亚洲精品色激情综合| 特级一级黄色大片| 欧美丝袜亚洲另类| 久久久久国内视频| 色av中文字幕| 欧美国产日韩亚洲一区| 免费观看的影片在线观看| 国产精品美女特级片免费视频播放器| 亚洲成人久久爱视频| 久久人妻av系列| 国产精品1区2区在线观看.| 亚洲成人精品中文字幕电影| АⅤ资源中文在线天堂| 久久国内精品自在自线图片| 国产91av在线免费观看| 日本熟妇午夜| 最好的美女福利视频网| 床上黄色一级片| 午夜福利视频1000在线观看| av黄色大香蕉| 国产精品无大码| 亚洲av免费高清在线观看| 国产探花极品一区二区| 亚洲精品粉嫩美女一区| 天堂影院成人在线观看| 97人妻精品一区二区三区麻豆| 99久久精品国产国产毛片| 欧美激情国产日韩精品一区| 一级毛片电影观看 | 精品久久国产蜜桃| 国产高清三级在线| 国产一区二区在线观看日韩| 久久99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 校园人妻丝袜中文字幕| 亚洲av中文字字幕乱码综合| 国产高清视频在线观看网站| 午夜精品国产一区二区电影 | 五月伊人婷婷丁香| 欧美色欧美亚洲另类二区| 欧美最新免费一区二区三区| 少妇裸体淫交视频免费看高清| 欧美一级a爱片免费观看看| 成人三级黄色视频| 狂野欧美白嫩少妇大欣赏| 小说图片视频综合网站| 欧美性感艳星| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 欧美中文日本在线观看视频| 欧美又色又爽又黄视频| 老司机影院成人| 日韩,欧美,国产一区二区三区 | 丰满的人妻完整版| 亚洲av电影不卡..在线观看| 深夜a级毛片| 午夜爱爱视频在线播放| 午夜福利在线观看免费完整高清在 | 亚洲精品粉嫩美女一区| 国产成人a区在线观看| 黄片wwwwww| 一进一出抽搐动态| 成人av一区二区三区在线看| 99久国产av精品国产电影| 亚洲av五月六月丁香网| 欧美日韩乱码在线| 97超碰精品成人国产| 久久中文看片网| 精品日产1卡2卡| 日韩av不卡免费在线播放| 一本久久中文字幕| 网址你懂的国产日韩在线| or卡值多少钱| a级一级毛片免费在线观看| 日韩欧美一区二区三区在线观看| 欧美成人免费av一区二区三区| 欧美+亚洲+日韩+国产| 国产亚洲av嫩草精品影院| 少妇熟女aⅴ在线视频| or卡值多少钱| 小蜜桃在线观看免费完整版高清| 有码 亚洲区| 老师上课跳d突然被开到最大视频| 成人三级黄色视频| 女人被狂操c到高潮| 中文字幕av成人在线电影| 日韩欧美精品免费久久| 成年版毛片免费区| а√天堂www在线а√下载| 国产欧美日韩精品一区二区| 国产一区亚洲一区在线观看| 久久久a久久爽久久v久久| 最新中文字幕久久久久| 成年女人毛片免费观看观看9| 嫩草影视91久久| 国产精品嫩草影院av在线观看| 成年版毛片免费区| а√天堂www在线а√下载| 97超碰精品成人国产| 长腿黑丝高跟| 精品久久久久久久久久免费视频| 精品乱码久久久久久99久播| 久久久久精品国产欧美久久久| www日本黄色视频网| 97超级碰碰碰精品色视频在线观看| 麻豆久久精品国产亚洲av| 亚洲五月天丁香| 搞女人的毛片| 日本免费a在线| 国产精品人妻久久久久久| 九色成人免费人妻av| 一本精品99久久精品77| 菩萨蛮人人尽说江南好唐韦庄 | 日日干狠狠操夜夜爽| 欧美日本视频| 18禁在线播放成人免费| 国产午夜精品论理片| 丰满乱子伦码专区| 看黄色毛片网站| 国产成人福利小说| 99久久精品一区二区三区| 免费av毛片视频| 午夜福利成人在线免费观看| 精品人妻视频免费看| 亚洲va在线va天堂va国产| 欧美潮喷喷水| 日产精品乱码卡一卡2卡三| 国产成年人精品一区二区| 国产黄色小视频在线观看| 国产中年淑女户外野战色| 人妻少妇偷人精品九色| 午夜福利成人在线免费观看| 伊人久久精品亚洲午夜| 亚洲精品一区av在线观看| 久久精品综合一区二区三区| 久久6这里有精品| 国产综合懂色| 成人三级黄色视频| 小蜜桃在线观看免费完整版高清| 哪里可以看免费的av片| 人妻少妇偷人精品九色| 变态另类成人亚洲欧美熟女| 国产精品女同一区二区软件| 亚洲在线观看片| 深夜a级毛片| 日韩 亚洲 欧美在线| 国产精品久久久久久久电影| 成人漫画全彩无遮挡| 观看免费一级毛片| 亚洲国产精品sss在线观看| 1024手机看黄色片| 国产精品亚洲美女久久久| 亚洲av成人av| 直男gayav资源| 搡老岳熟女国产| 国产片特级美女逼逼视频| 麻豆一二三区av精品| 亚洲自拍偷在线| 久久午夜亚洲精品久久| 国产伦精品一区二区三区四那| 亚洲成人中文字幕在线播放| 午夜激情欧美在线| 插逼视频在线观看| 舔av片在线| 麻豆av噜噜一区二区三区| 欧美绝顶高潮抽搐喷水| 国内揄拍国产精品人妻在线| 精品免费久久久久久久清纯| 最新在线观看一区二区三区| 99视频精品全部免费 在线| 久久精品影院6| 日本三级黄在线观看| 国内少妇人妻偷人精品xxx网站| 欧美xxxx黑人xx丫x性爽| 精品久久久久久久人妻蜜臀av| 欧美日韩精品成人综合77777| 最近手机中文字幕大全| 一级av片app| 伦理电影大哥的女人| 免费不卡的大黄色大毛片视频在线观看 | 国产一区亚洲一区在线观看| 乱系列少妇在线播放| 18+在线观看网站| 男女做爰动态图高潮gif福利片| 国产亚洲av嫩草精品影院| 欧美激情国产日韩精品一区| 高清日韩中文字幕在线| 简卡轻食公司| 老熟妇乱子伦视频在线观看| 波多野结衣高清无吗| 毛片一级片免费看久久久久| 婷婷亚洲欧美| 18禁在线无遮挡免费观看视频 | 久99久视频精品免费| 你懂的网址亚洲精品在线观看 | 人人妻人人澡欧美一区二区| 丰满乱子伦码专区| 国产黄色小视频在线观看| 亚洲在线观看片| 亚洲av.av天堂| 免费看av在线观看网站| 12—13女人毛片做爰片一| 1024手机看黄色片| 亚洲欧美精品自产自拍| 男女视频在线观看网站免费| 久久久久国产网址| 在线观看美女被高潮喷水网站| 亚洲三级黄色毛片| 美女黄网站色视频| 99热这里只有是精品在线观看| 久久欧美精品欧美久久欧美| 午夜老司机福利剧场| 老熟妇乱子伦视频在线观看| 久久久精品欧美日韩精品| 免费av观看视频| 麻豆国产97在线/欧美| 十八禁国产超污无遮挡网站| 亚洲综合色惰| 亚洲不卡免费看| 3wmmmm亚洲av在线观看| 秋霞在线观看毛片| 亚洲欧美精品自产自拍| 成人午夜高清在线视频| 欧美潮喷喷水| eeuss影院久久| 高清毛片免费观看视频网站| 亚洲精华国产精华液的使用体验 | av中文乱码字幕在线| 午夜福利视频1000在线观看| 日韩在线高清观看一区二区三区| 少妇被粗大猛烈的视频| 久久久久性生活片| 成人亚洲欧美一区二区av| 国产一区二区三区在线臀色熟女| 老女人水多毛片| 日韩欧美精品v在线| 国产片特级美女逼逼视频| av天堂中文字幕网| 亚洲国产欧美人成| 国产成人a区在线观看| 久久人人精品亚洲av| 日韩亚洲欧美综合| 亚洲中文字幕日韩| 97热精品久久久久久| 亚洲性久久影院| 国国产精品蜜臀av免费| 亚洲无线观看免费| 美女内射精品一级片tv| 一区二区三区高清视频在线| 欧美3d第一页| 久久久久久伊人网av| 青春草视频在线免费观看| 97在线视频观看| 欧美激情国产日韩精品一区| 欧美成人一区二区免费高清观看| 不卡一级毛片| 色视频www国产| 欧美日韩精品成人综合77777| 在线播放国产精品三级| 男女边吃奶边做爰视频| 伊人久久精品亚洲午夜| 午夜免费激情av| 天天躁夜夜躁狠狠久久av| 秋霞在线观看毛片| 久久亚洲精品不卡| 99热精品在线国产| 搡老妇女老女人老熟妇| 你懂的网址亚洲精品在线观看 | 91狼人影院| 久久精品综合一区二区三区| 伊人久久精品亚洲午夜| 高清毛片免费观看视频网站| 91午夜精品亚洲一区二区三区| 一级av片app| 色综合亚洲欧美另类图片| 亚洲婷婷狠狠爱综合网| 不卡一级毛片| 秋霞在线观看毛片| 精品国产三级普通话版| 啦啦啦观看免费观看视频高清| 中文字幕熟女人妻在线| 亚洲国产精品合色在线| 国内精品宾馆在线| 欧美日本亚洲视频在线播放| 蜜桃久久精品国产亚洲av| 亚洲国产欧洲综合997久久,| 国产午夜精品论理片| 亚洲婷婷狠狠爱综合网| 欧美一区二区精品小视频在线| 免费看av在线观看网站| 国产私拍福利视频在线观看| 午夜免费男女啪啪视频观看 | av.在线天堂| 国产美女午夜福利| 18禁黄网站禁片免费观看直播| 国产乱人偷精品视频| 久久婷婷人人爽人人干人人爱| 国产v大片淫在线免费观看| 精品国产三级普通话版| 国产精品电影一区二区三区| 狠狠狠狠99中文字幕| 日本黄色视频三级网站网址| 在线看三级毛片| 无遮挡黄片免费观看| 最新在线观看一区二区三区| 成人三级黄色视频| 亚洲人成网站高清观看| 偷拍熟女少妇极品色| 欧美日韩国产亚洲二区| 少妇裸体淫交视频免费看高清| 久久鲁丝午夜福利片| 国产亚洲精品久久久com| av中文乱码字幕在线| 久久久久久久久久久丰满| 哪里可以看免费的av片| 人妻制服诱惑在线中文字幕| 在现免费观看毛片| 国产 一区精品| 国产高清三级在线| 免费大片18禁| 国产午夜福利久久久久久| 中文字幕精品亚洲无线码一区| 日韩中字成人| 淫妇啪啪啪对白视频| 欧美一区二区亚洲| 亚洲国产高清在线一区二区三| 日本撒尿小便嘘嘘汇集6| 久久久成人免费电影| 成人精品一区二区免费| 免费人成视频x8x8入口观看| ponron亚洲| 晚上一个人看的免费电影| 国产精品一区二区三区四区久久| 一个人看视频在线观看www免费| 一边摸一边抽搐一进一小说| 此物有八面人人有两片| 尤物成人国产欧美一区二区三区| 国产精品日韩av在线免费观看| 亚洲欧美日韩东京热| 日韩亚洲欧美综合| 国产一区二区在线观看日韩| 国产欧美日韩精品一区二区| 卡戴珊不雅视频在线播放| 国产免费一级a男人的天堂| 亚洲欧美精品综合久久99| 久久人人爽人人爽人人片va| 日韩精品中文字幕看吧| 麻豆国产av国片精品| eeuss影院久久| 可以在线观看毛片的网站| 卡戴珊不雅视频在线播放| 天天躁夜夜躁狠狠久久av| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区激情短视频| 网址你懂的国产日韩在线| 热99re8久久精品国产| 一区福利在线观看| 国产精品野战在线观看| 亚洲婷婷狠狠爱综合网| 亚洲真实伦在线观看| 久久人妻av系列| 欧美日韩乱码在线| 国产免费一级a男人的天堂| 一卡2卡三卡四卡精品乱码亚洲| 天堂动漫精品| 搞女人的毛片| 亚洲国产精品成人久久小说 | 卡戴珊不雅视频在线播放| 国产乱人视频| 亚洲不卡免费看| or卡值多少钱| 亚洲电影在线观看av| 中出人妻视频一区二区| 精品久久国产蜜桃| av福利片在线观看| 秋霞在线观看毛片| 日韩欧美三级三区| 亚洲,欧美,日韩| 久久九九热精品免费| www日本黄色视频网| a级一级毛片免费在线观看| 国产精品久久视频播放| 欧美日本亚洲视频在线播放| 亚洲在线观看片| 日韩av在线大香蕉| 人人妻,人人澡人人爽秒播| 在线播放国产精品三级| 免费大片18禁| 一区二区三区高清视频在线| 别揉我奶头 嗯啊视频| 99九九线精品视频在线观看视频| 国产成人精品久久久久久| 国产v大片淫在线免费观看| 男女边吃奶边做爰视频| 九九在线视频观看精品| 成年女人毛片免费观看观看9| 精品久久久久久久久av| 国产高清不卡午夜福利| 国产单亲对白刺激| 国产一区二区激情短视频| 国产视频一区二区在线看| 成年女人毛片免费观看观看9| 女生性感内裤真人,穿戴方法视频| 一级黄色大片毛片| 桃色一区二区三区在线观看| 国产男人的电影天堂91| 亚洲av二区三区四区| 久久人人爽人人片av| 精品午夜福利在线看| 国产白丝娇喘喷水9色精品| 亚洲成a人片在线一区二区| 久久人人爽人人爽人人片va| 国产精品1区2区在线观看.| 亚洲一区高清亚洲精品| 97超视频在线观看视频| 婷婷亚洲欧美|