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

    基于離散元法的椰糠顆粒建模及仿真參數(shù)標(biāo)定

    2021-11-08 12:41:34王剛王倩薛忠郭昌進(jìn)宋剛王槊
    熱帶作物學(xué)報(bào) 2021年9期
    關(guān)鍵詞:椰糠

    王剛 王倩 薛忠 郭昌進(jìn) 宋剛 王槊

    摘? 要:為確定椰糠在離散元模擬過(guò)程中合理的仿真參數(shù),本文基于離散元法Hertz-Mindlin(no slip)接觸模型在仿真軟件EDEM 2018中建立椰糠顆粒模型并生成顆粒工廠,通過(guò)測(cè)量比對(duì)物理試驗(yàn)和仿真試驗(yàn)椰糠堆積角的方法對(duì)其仿真參數(shù)進(jìn)行標(biāo)定。首先通過(guò)物理試驗(yàn)測(cè)得椰糠堆積角、堆積密度等本征參數(shù),利用EDEM 2018內(nèi)嵌GEMM數(shù)據(jù)庫(kù)以及相關(guān)文獻(xiàn),綜合分析得到椰糠待標(biāo)定因素的高低水平;然后進(jìn)行Plackett-Burman試驗(yàn)設(shè)計(jì)篩選出對(duì)椰糠堆積角影響顯著的因素依次為:椰糠-椰糠靜摩擦系數(shù)、椰糠-椰糠滾動(dòng)摩擦系數(shù)和椰糠剪切模量;再進(jìn)行Box-Behnken試驗(yàn)設(shè)計(jì)建立堆積角與3個(gè)顯著性影響因素的回歸模型,運(yùn)用Design expert軟件優(yōu)化功能,以堆積角物理試驗(yàn)45.69°為目標(biāo),對(duì)回歸方程尋優(yōu)求解,得到顯著性影響因素最佳組合:椰糠剪切模量為1.44 MPa,椰糠-椰糠靜摩擦系數(shù)為1.1 椰糠-椰糠滾動(dòng)摩擦系數(shù)為0.15;最后通過(guò)試驗(yàn)驗(yàn)證表明此組合參數(shù)可用于椰糠物料的離散元仿真,能為椰糠輸送、混合等機(jī)械裝備的設(shè)計(jì)研發(fā)提供理論參考。

    關(guān)鍵詞:椰糠;離散元;參數(shù)標(biāo)定;堆積角;試驗(yàn)設(shè)計(jì)

    中圖分類號(hào):S22????? 文獻(xiàn)標(biāo)識(shí)碼:A

    Modeling and Simulation Parameters Calibration of Coconut Coir Particles Based on DEM

    WANG Gang1, WANG Qian1, XUE Zhong2, GUO Changjin2, SONG Gang1, WANG Shuo1

    1. Tropical Agricultural Machinery Research Institute, Chinese Academy of Tropical Agriculture Sciences, Zhanjiang, Guangdong 524091, China; 2. Key Laboratory of Agricultural Equipment for Tropical Crops, Ministry of Agriculture and Rural Affairs, Zhanjiang, Guangdong 524091, China

    Abstract: In order to determine the reasonable simulation parameters of coconut coir in the DEM process, based on the DEM Mindlin (no slip) contact model, a coconut coir model is established in the simulation software EDEM 2018 and the particle factory is generated. The simulation parameters were calibrated by comparing the physical test and virtual test of coconut coir repose angle. Firstly, the intrinsic parameters of coconut coir such as repose angle and bulk density were measured by physical tests, and the level of factors to be calibrated were obtained by comprehensive analysis using the GEMM database embedded in EDEM 2018 and relevant literature. A Plackett–Burman test design was carried out to screen out the factors that have significant influence on the repose angle of coconut coir: static friction coefficient of coconut coir with coconut coir, rolling friction coefficient of coconut coir with coconut coir, shear modulus of coconut coir. Through Box–Behnken experimental design, a regression model of repose angle and three significant influencing factors were established, using the optimization function of design expert software, taking the physical test of stacking angle of 45.69° as the objective, the regression equation was optimized and solved, and the best combination of significant influencing factors were obtained. Coconut coir shear modulus was 1.44 MPa, static friction coefficient of coconut coir with coconut coir was 1.11, and rolling friction coefficient of coconut coir with coconut coir was 0.15. The experimental results showed that the combined parameters could be used in the discrete element simulation of coconut coir materials, and provide theoretical references for the design of coconut coir conveying, mixing and other mechanical equipment.

    Keywords: coconut coir; discrete element method; parameters calibration; repose angle; experimental design

    DOI: 10.3969/j.issn.1000-2561.2021.09.036

    椰糠是椰子加工過(guò)程中的副產(chǎn)物,是從椰子外殼纖維脫落下來(lái)的一種天然有機(jī)質(zhì)媒介,椰糠因其保溫、保濕、疏松、透氣的優(yōu)良理化性質(zhì),在設(shè)施農(nóng)業(yè)中正逐步取代草炭成為應(yīng)用最為廣泛的栽培基質(zhì)[1-4]。研究表明,單一依靠椰糠作為栽培基質(zhì)是存在缺陷的,通常需要與其他物料如有機(jī)肥等按照一定比例搭配形成混合基質(zhì),才能保證栽培基質(zhì)的碳氮比、養(yǎng)分含量、保水通氣性等性狀更適應(yīng)作物生長(zhǎng)[5-6]。

    混合機(jī)械設(shè)備是椰糠栽培基質(zhì)機(jī)械化生產(chǎn)的關(guān)鍵設(shè)備之一,其生產(chǎn)過(guò)程屬于固體混合。因?yàn)楣腆w混合區(qū)別于液體攪拌,不具有擴(kuò)散特性,須施加外力才能強(qiáng)制流動(dòng),固體顆粒與混合設(shè)備(或其相關(guān)接觸部件)的接觸關(guān)系、顆粒運(yùn)動(dòng)特性以及微觀作用機(jī)理等直接關(guān)系到混合設(shè)備的作業(yè)性能和工作效率,所以不同的固體混合設(shè)備沒(méi)有統(tǒng)一的制造標(biāo)準(zhǔn),往往根據(jù)不同的功能和使用要求單獨(dú)設(shè)計(jì),而設(shè)計(jì)專用設(shè)備,就要充分了解物料特性[7-8]。

    目前,實(shí)際應(yīng)用的椰糠混合基質(zhì)生產(chǎn)設(shè)備大多為直接套用市面普通的飼料混合設(shè)備或粉粒體混合設(shè)備,生產(chǎn)過(guò)程中存在混合不均勻或功率不符等問(wèn)題。現(xiàn)階段仍缺乏對(duì)椰糠基質(zhì)混合過(guò)程及混合性能的系統(tǒng)研究,對(duì)椰糠混合基質(zhì)的研究也主要集中在基質(zhì)不同配比對(duì)農(nóng)作物生長(zhǎng)特性的影響規(guī)律等方面,椰糠作為一種典型散體物料,科研院所對(duì)其開(kāi)展物料離散元特性方面的研究較少,而近年來(lái),針對(duì)土壤、顆粒肥料、農(nóng)業(yè)物料等離散性物料特性方面的研究相對(duì)較多,并廣泛運(yùn)用了離散元法(discrete element method,簡(jiǎn)稱DEM)進(jìn)行分析,可對(duì)椰糠物料離散元特性的分析提供重要的借鑒與參考。在土壤方面,邢潔潔等[9]通過(guò)試驗(yàn)并利用離散元建模仿真確定磚紅壤顆粒特性最佳組合,為研發(fā)適用于海南熱區(qū)磚紅壤地高性能耕作裝備關(guān)鍵觸土部件的設(shè)計(jì)優(yōu)化提供支撐;Ucgul等[10]利用離散元軟件中Hertze- Mindlin及Hysteretic Spring接觸模型,研究分析了有黏結(jié)力和無(wú)黏結(jié)力的情況下土壤顆粒的塑性形變等問(wèn)題。在顆粒肥料方面,羅帥等[11]使用離散元參數(shù)標(biāo)定方法篩選蚯蚓糞基質(zhì)顆粒參數(shù),為研究蚯蚓堆肥處理及蚓糞利用相關(guān)設(shè)備奠定基礎(chǔ);Bangura等[12]采用離散元法對(duì)肥料進(jìn)行建模通過(guò)分析排肥過(guò)程比較了螺旋槽輪和直槽輪的排肥速度與排肥均勻性。在農(nóng)業(yè)物料方面,鄭智旗[13]運(yùn)用離散元分析粉碎后的玉米秸稈在輸送裝置內(nèi)的運(yùn)動(dòng)規(guī)律,通過(guò)改變拋送板葉片傾角、旋轉(zhuǎn)軸轉(zhuǎn)速和喂入量等模擬研究參數(shù)改變對(duì)輸送性能的影響;Ghodki等[14]建立了大豆離散元模型,通過(guò)與大豆形狀、靜止角等特征進(jìn)行比較,獲得了該模型的接觸參數(shù)。因此,本文在對(duì)國(guó)內(nèi)外散體物料離散元參數(shù)標(biāo)定方法運(yùn)用進(jìn)行總結(jié)的基礎(chǔ)上,針對(duì)椰糠物料具體實(shí)際,采用物理試驗(yàn)和建模仿真相結(jié)合的手段,基于離散元法分析確定椰糠的顆粒物料特性參數(shù),為后續(xù)設(shè)計(jì)高效的椰糠栽培基質(zhì)混合專用設(shè)備、優(yōu)化混合結(jié)構(gòu)關(guān)鍵部件和混合工藝參數(shù)、提高椰糠栽培基質(zhì)生產(chǎn)效率與產(chǎn)品品質(zhì)提供理論依據(jù)。

    1? 材料與方法

    1.1? 仿真原理及模型選取

    椰糠離散元模擬中需要的物料特性參數(shù)主要分為三類:一是材料本征參數(shù),包括顆粒大小、堆積密度、剪切模量和泊松比等,這是椰糠自身的特性參數(shù),和外界無(wú)關(guān);二是材料接觸參數(shù),包括碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動(dòng)摩擦系數(shù),這是椰糠顆粒相互接觸才會(huì)發(fā)生作用的參數(shù),與接觸雙方都有關(guān)系;三是接觸模型參數(shù),跟接觸模型選取有關(guān),它描述椰糠顆粒之間接觸時(shí)顆粒的行為,一些特殊的接觸模型需要特殊的模型參數(shù)[15-17]。本文選取“Hertz-Mindlin (no slip)”顆粒接觸模型,該模型在力的計(jì)算方面高效準(zhǔn)確[18-19],其基本原理如圖1所示。顆粒之間作用力包括法向彈性力 、法向阻尼力 、切向彈性力 、切向阻尼力 、滾動(dòng)摩擦力 。

    (1)

    其中

    (2)

    (3)

    式中: 為等效彈性模量, 為法向重疊量, 為等效接觸半徑, 、 為接觸顆粒的泊松比, 、 為接觸顆粒的彈性模量, 、 為接觸顆粒的接觸半徑。

    (4)

    其中

    (5)

    (6)

    (7)

    式中: 為阻尼比, 為法向剛度, 為等效質(zhì)量, 為相對(duì)速度的法向分量, 為恢復(fù)系數(shù), 為接觸顆粒的質(zhì)量。

    (8)

    其中

    (9)

    式中: 為法向剛度, 為法向重疊量, 為等效剪切模量。

    (10)

    式中: 為相對(duì)速度的切向分量。

    (11)

    式中: 為滾動(dòng)摩擦系數(shù), 為接觸點(diǎn)到顆粒質(zhì)心的距離, 為顆粒在接觸點(diǎn)處單位角速度矢量。

    1.2? 試驗(yàn)材料與物理試驗(yàn)

    1.2.1? 試驗(yàn)材料選取? 試驗(yàn)中所用的材料采用椰糠磚充分泡發(fā)后的椰糠物料,椰糠磚尺寸為200 mm×100 mm×50 mm、質(zhì)量為(550±10)g、泡發(fā)量為12 L/kg,相關(guān)術(shù)語(yǔ)取自由中國(guó)農(nóng)業(yè)機(jī)械化協(xié)會(huì)設(shè)施農(nóng)業(yè)分會(huì)發(fā)布的團(tuán)體標(biāo)準(zhǔn)T/CAMA 12—2019《無(wú)土栽培椰糠》[20]。

    1.2.2? 含水率測(cè)量? 在待測(cè)的椰糠物料堆中采用五點(diǎn)法取樣,每個(gè)取樣點(diǎn)稱取3~5 g樣品,分別放入MA45水分分析測(cè)試儀的樣品盤中,并確保樣品在樣品盤中均勻分布。根據(jù)5次樣品測(cè)量的平均值,計(jì)算出待測(cè)椰糠的含水率為86.57%。

    1.2.3? 堆積密度測(cè)量? 稱取一定質(zhì)量的椰糠將其沿喂料漏斗自然下落到量筒中,通過(guò)用質(zhì)量除以量筒內(nèi)椰糠的體積讀數(shù)計(jì)算出待測(cè)椰糠的堆積密度,重復(fù)試驗(yàn)5次取平均值為370 kg/m3。

    1.2.4? 真實(shí)密度測(cè)量? 稱取一定質(zhì)量的椰糠裝入盛有一定體積水的100 mL量筒中,量筒增加的體積讀數(shù)就是椰糠的真實(shí)體積,用椰糠質(zhì)量除以增加的體積,重復(fù)5次試驗(yàn),計(jì)算取平均值得到椰糠的真實(shí)密度為1 074 kg/m3。

    1.2.5? 顆粒尺寸測(cè)定? 稱取100 g椰糠樣品,使用分級(jí)篩進(jìn)行顆粒分級(jí),用電子天平稱取不同粒徑范圍對(duì)應(yīng)的質(zhì)量,重復(fù)5次試驗(yàn),獲得椰糠粒徑質(zhì)量分?jǐn)?shù)分布,為仿真時(shí)椰糠顆粒建模和顆粒工廠生成提供依據(jù),結(jié)果如表1所示。

    1.2.6? 堆積角的測(cè)定? 堆積角也稱安息角或休止角,是表征顆粒物料流動(dòng)、摩擦等特性的宏觀參數(shù),與顆粒表面積、顆粒密度、顆粒摩擦系數(shù)及含水率均有關(guān)系[21-22]。當(dāng)大量顆粒物料被傾倒于水平面上堆積為錐體,錐體表面與水平面所成內(nèi)角即為堆積角,堆積角物理試驗(yàn)是顆粒物料離散元參數(shù)標(biāo)定的常用手段。采用注入法將椰糠從漏斗上方慢慢加入,漏出的物料在水平面上形成圓錐狀堆積體的傾斜角即為椰糠的堆積角,重復(fù)試驗(yàn)5次,取平均值為45.69°。

    1.3? 顆粒建模與參數(shù)選取

    1.3.1? 建立椰糠仿真模型? 由上述表1可以看出,在選取的椰糠樣本中,顆粒粒徑1.25~5 mm的質(zhì)量占總質(zhì)量的93.18%,其中粒徑1.25~2 mm的質(zhì)量分?jǐn)?shù)為60.38%,占大多數(shù)。在離散元參數(shù)標(biāo)定中,若完全按照顆粒真實(shí)尺寸進(jìn)行建模,仿真效率太低,研究表明:相比真實(shí)的顆粒尺寸與形狀,顆粒之間以及顆粒與幾何體之間的相互作用更為關(guān)鍵[23-24]。因此,在建立椰糠顆粒模型時(shí),以粒徑1.25~2 mm的顆粒為基礎(chǔ),考慮縮短數(shù)值模擬時(shí)間和簡(jiǎn)化模型因素,結(jié)合椰糠實(shí)際外觀形態(tài)創(chuàng)建其離散元顆粒模型,如圖2所示。

    1.3.2? 標(biāo)定參數(shù)范圍確定? 椰糠物料特性的研究目前尚不完善,其離散元仿真的參數(shù)尤其匱乏。EDEM軟件內(nèi)置顆粒材料數(shù)據(jù)庫(kù)(Generic EDEM Material Model Database)中包含了數(shù)千種具有代表性的顆粒物料,如巖石、土壤和礦石等。通過(guò)選擇仿真規(guī)模和堆積密度范圍、輸入物料堆積角,獲得椰糠顆粒之間的接觸參數(shù)參考值范圍:恢復(fù)系數(shù)為0.15~0.75,靜摩擦系數(shù)為0.20~1.16,滾動(dòng)摩擦系數(shù)為0.05~0.15。考慮到椰糠作為農(nóng)業(yè)物料及其作為栽培基質(zhì)與土壤接近,主要參照文獻(xiàn)[25–30]中各類物料參數(shù)取值或范圍,綜合考慮設(shè)定椰糠泊松比為0.1~0.5,剪切模量為1.0~10.0 MPa,椰糠與鋼材恢復(fù)系數(shù)為0.1~0.7、靜摩擦系數(shù)為0.20~0.90、滾動(dòng)摩擦系數(shù)為0.05~0.30。通過(guò)查閱文獻(xiàn)取鋼材泊松比為0. 剪切模量為7.9×1010 Pa,密度為7 850 kg/m 重力加速度為9.81 m/s2。各參數(shù)標(biāo)定試驗(yàn)因素及水平如表2所示。

    1.4? 數(shù)據(jù)處理

    采用Design Expert 10軟件設(shè)計(jì)試驗(yàn)方案并對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行方差分析和顯著性分析,因試驗(yàn)

    參數(shù)較多,故先設(shè)計(jì)Plackett-Burman試驗(yàn)篩選出對(duì)堆積角影響顯著的試驗(yàn)因素,然后根據(jù)試驗(yàn)結(jié)果對(duì)顯著性影響因素高(1)、中(0)、低(-1)3個(gè)水平設(shè)計(jì)Box-Behnken試驗(yàn)進(jìn)行響應(yīng)面分析,得到堆積角和顯著性影響因素之間的回歸模型,最后通過(guò)參數(shù)優(yōu)化尋找最優(yōu)解和試驗(yàn)驗(yàn)證完成參數(shù)標(biāo)定。

    2? 結(jié)果與分析

    2.1? Plackett-Burman 試驗(yàn)設(shè)計(jì)與結(jié)果

    Plackett-Burman設(shè)計(jì)法是一種兩水平的試驗(yàn)設(shè)計(jì)方法,通過(guò)比較各個(gè)因素兩水平的差異與整體的差異來(lái)確定因素的顯著性。以椰糠堆積角為響應(yīng)值,延用表2中8個(gè)因素高低水平,并設(shè)置3個(gè)虛擬參數(shù)用于誤差分析,共進(jìn)行12次試驗(yàn),試驗(yàn)設(shè)計(jì)及結(jié)果如表3所示。

    通過(guò)方差分析得到8個(gè)因素對(duì)堆積角的影響效果和顯著性,如表4所示。分析可得:A、B、D、E、G對(duì)椰糠顆粒的堆積角有正效應(yīng),即堆積角隨該5個(gè)因素的增大而增大;C、F、H則有負(fù)效應(yīng),即堆積角隨該3個(gè)因素的增大而減小;其中對(duì)堆積角影響極顯著(P<0.01)的因素是椰糠-椰糠靜摩擦系數(shù)(D)和椰糠-椰糠滾動(dòng)摩擦系數(shù)(E),影響顯著(P<0.05)的因素是椰糠剪切模量(B),而其余因素影響不顯著;模型的P<0.0 說(shuō)明回歸模型極顯著,預(yù)測(cè)值 =0.994 0,說(shuō)明模型能很好地反映各因素的變化情況。

    2.2? Box-Behnken試驗(yàn)設(shè)計(jì)與結(jié)果

    2.2.1? 試驗(yàn)設(shè)計(jì)與結(jié)果? Box-Behnken設(shè)計(jì)是一種響應(yīng)曲面設(shè)計(jì)類型,根據(jù)Plackett-Burman試驗(yàn)結(jié)果,取椰糠泊松比為0.3、椰糠-椰糠恢復(fù)系數(shù)為0.45、椰糠和鋼材恢復(fù)系數(shù)為0.4、椰糠和鋼材靜摩擦系數(shù)為0.55、椰糠和鋼材滾動(dòng)摩擦系數(shù)為0.18,以椰糠-椰糠靜摩擦系數(shù)(D)、椰糠-椰糠滾動(dòng)摩擦系數(shù)(E)、椰糠剪切模量(B)為試驗(yàn)因素,堆積角為響應(yīng)指標(biāo),延用表2中3個(gè)因素的高中低水平,進(jìn)行Box-Behnken試驗(yàn)設(shè)計(jì)共進(jìn)行17模擬仿真,結(jié)果如表5所示。

    2.2.2? 堆積角回歸模型分析? 應(yīng)用Design expert軟件建立堆積角與3個(gè)顯著性影響因素的二階模型:

    (12)

    該回歸方程的P=0.00 決定系數(shù) = 0.9492,校正決定系數(shù) =0.8838,信噪比為13.47,說(shuō)明該堆積角回歸方程極顯著,相關(guān)性好,可以對(duì)目標(biāo)堆積角進(jìn)行預(yù)測(cè);變異系數(shù)CV=3.43%,失擬項(xiàng)P=0.4007>0.05,表明該方程擬合程度好,對(duì)該模型試驗(yàn)結(jié)果進(jìn)行回歸方差分析,如表6所示。結(jié)果表明,椰糠-椰糠靜摩擦系數(shù)(D)、椰糠-椰糠滾動(dòng)摩擦系數(shù)(E)對(duì)堆積角影響極顯著;椰糠-椰糠靜摩擦系數(shù)的二次方項(xiàng)( )、椰糠-椰糠滾動(dòng)摩擦系數(shù)的二次方項(xiàng)( )對(duì)堆積角影響顯著。

    2.2.3? 參數(shù)優(yōu)化與試驗(yàn)驗(yàn)證? 運(yùn)用Design expert軟件優(yōu)化功能,以椰糠堆積角試驗(yàn)值45.69°為目標(biāo)值,在設(shè)置的試驗(yàn)因素水平范圍內(nèi)對(duì)堆積角回歸模型進(jìn)行尋優(yōu),所得優(yōu)化解若干組。以優(yōu)化解作為因素水平值進(jìn)行堆積角仿真試驗(yàn),并與物理試驗(yàn)堆積角進(jìn)行對(duì)比驗(yàn)證,得到試驗(yàn)形狀相似的一組最優(yōu)解,即椰糠剪切模量(B)為1.44 MPa,椰糠-椰糠靜摩擦系數(shù)為(D)為1.1 椰糠-椰糠滾動(dòng)摩擦系數(shù)(E)為0.15,其余非顯著性因素取中間水平。此優(yōu)化解下的堆積角仿真試驗(yàn)結(jié)果45.25°與物理試驗(yàn)結(jié)果45.69°對(duì)比如圖3所示,仿真結(jié)果與真實(shí)情況一致,說(shuō)明該優(yōu)化參數(shù)可為樣品椰糠與混合設(shè)備的進(jìn)一步離散元仿真提供依據(jù)。

    3? 討論

    (1)利用椰糠+有機(jī)肥制成混合基質(zhì)替代草炭,不僅可以“變廢為寶”,避免生態(tài)資源的枯竭,還可以減少環(huán)境的二次污染。采用離散元法研究機(jī)械設(shè)備與物料相互作用的研究雖然不少,在混合設(shè)備領(lǐng)域也有相關(guān)報(bào)道,但主要集中在飼料、顆粒肥料或者混凝土的研究方面,缺乏針對(duì)椰糠基質(zhì)綜合利用的混合攪拌研究?,F(xiàn)有文獻(xiàn)中針對(duì)椰糠物料特性的研究主要分布在結(jié)構(gòu)及比例、分解度、灰分含量、孔隙度、pH、礦質(zhì)營(yíng)養(yǎng)等方面,如前所述椰糠離散元模擬中所需的本征參數(shù)、接觸參數(shù)等必備特性參數(shù)缺乏數(shù)據(jù)支撐,而這些參數(shù)以試驗(yàn)方法測(cè)量往往難以獲得。本文選擇利用EDEM 2018離散元軟件建立椰糠顆粒模型,預(yù)設(shè)相關(guān)參數(shù)進(jìn)行椰糠堆積角的模擬仿真,通過(guò)與物理試驗(yàn)相結(jié)合進(jìn)行參數(shù)標(biāo)定,確定椰糠離散元模擬所需的準(zhǔn)確參數(shù),研究結(jié)果在國(guó)內(nèi)外椰糠加工設(shè)備基礎(chǔ)研究方面尚未見(jiàn)報(bào)道。

    (2)本文通過(guò)分析文獻(xiàn)以及結(jié)合前期試驗(yàn)測(cè)量的堆積角,根據(jù)EDEM 2018材料庫(kù)確定椰糠泊松比、剪切模量,椰糠-椰糠恢復(fù)系數(shù)、靜摩擦系數(shù)、滾動(dòng)摩擦系數(shù),椰糠-鋼材恢復(fù)系數(shù)、靜摩擦系數(shù)、滾動(dòng)摩擦系數(shù)這8個(gè)因素的取值范圍。然后進(jìn)行了Plackett-Burman試驗(yàn)設(shè)計(jì)篩選出對(duì)椰糠堆積角影響顯著的接觸參數(shù)和本征參數(shù)包括:椰糠-椰糠靜摩擦系數(shù)、椰糠-椰糠滾動(dòng)摩擦系數(shù)和椰糠剪切模量。通過(guò)Box-Behnken試驗(yàn)設(shè)計(jì)建立了堆積角與顯著性影響因素的回歸模型,方差分析表明了椰糠-椰糠靜摩擦系數(shù)和椰糠-椰糠滾動(dòng)摩擦系數(shù)對(duì)堆積角影響極顯著,椰糠-椰糠靜摩擦系數(shù)的二次方項(xiàng)、椰糠-椰糠滾動(dòng)摩擦系數(shù)的二次方項(xiàng)對(duì)堆積角影響顯著。以椰糠物理試驗(yàn)堆積角為目標(biāo),對(duì)回歸方程尋優(yōu)求解,得到顯著性影響因素最佳組合為椰糠剪切模量1.44 MPa,椰糠-椰糠靜摩擦系數(shù)1.1 椰糠-椰糠滾動(dòng)摩擦系數(shù)0.15。最后通過(guò)最優(yōu)組合參數(shù)試驗(yàn)驗(yàn)證,證實(shí)了標(biāo)定的參數(shù)準(zhǔn)確、可靠,建立的模型正確、可行。

    (3)離散元法可以準(zhǔn)確高效地分析顆粒的運(yùn)動(dòng)情況,并具有可視化特點(diǎn),本文研究結(jié)果對(duì)椰糠輸送、混合等機(jī)械裝備的研發(fā)提供了基礎(chǔ)數(shù)據(jù)。但農(nóng)業(yè)物料的各種特性往往與含水率有很大關(guān)系,本文取樣以椰糠磚泡發(fā)量為依據(jù),雖盡可能縮小了與實(shí)際應(yīng)用中的差距,但受不同栽培基質(zhì)配方特性的要求,需進(jìn)一步摸清不同含水率對(duì)椰糠堆積角的影響規(guī)律,從而確定不同含水率下椰糠物料的特性參數(shù),實(shí)現(xiàn)工業(yè)化生產(chǎn)。另外,上述試驗(yàn)得出:對(duì)椰糠堆積角影響顯著的是椰糠-椰糠靜摩擦系數(shù)、椰糠-椰糠滾動(dòng)摩擦系數(shù)和椰糠剪切模量,它們都是顆粒物料特性,與鋼材特性關(guān)聯(lián)不大。而椰糠混合基質(zhì)生產(chǎn)肯定至少包含2種物料,如椰糠和有機(jī)肥等,那么有必要進(jìn)一步研究椰糠-有機(jī)肥的顆粒接觸特性,這樣才能使整個(gè)離散元仿真模擬系統(tǒng)更加完善。

    參考文獻(xiàn)

    [1]? 譚海燕, 張振榮, 林峰瑩, 等. 椰糠基質(zhì)限根栽培下櫻桃番茄肥料利用率研究[J]. 湖北農(nóng)業(yè)科學(xué), 2020, 59(17): 46-48.

    [2] 邱志豪, 湯柔穎, 韓瑩琰, 等. 椰糠與蛭石混配基質(zhì)理化性狀及其對(duì)生菜生長(zhǎng)的影響[J]. 北京農(nóng)學(xué)院學(xué)報(bào), 2019, 34(4): 62-65.

    [3] 肖守華, 趙? 西, 肖真真, 等. 以椰糠為基質(zhì)的設(shè)施甜瓜無(wú)土栽培基質(zhì)配方篩選[J]. 山東農(nóng)業(yè)科學(xué), 2019, 51(1): 61-64.

    [4]? Barbara M, Sofia M, Sabrina R, et al. Coconut coir as a sustainable nursery growing media for seedling production of the ecologically diverse quercus species[J]. Forests, 2020, 11(5): 522.

    [5]? 劉? 佳, 季延海, 王寶駒, 等. 椰糠復(fù)合基質(zhì)對(duì)溫室番茄生長(zhǎng)及品質(zhì)的影響[J]. 江蘇農(nóng)業(yè)科學(xué), 2019, 47(17): 150-154.

    [6]? 王躍華, 張明科, 惠麥俠, 等. 不同椰糠配比基質(zhì)對(duì)白菜幼苗生長(zhǎng)的影響[J]. 南方農(nóng)業(yè)學(xué)報(bào), 2019, 50(12): 2749-2754.

    [7]? 李萍萍, 吳? 碩, 張西良, 等. 番茄秸稈與發(fā)酵菌種立式螺旋混合設(shè)備參數(shù)優(yōu)化[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2016, 47(11): 114-120.

    [8]? 吳? 碩. 番茄秸稈混料立式螺旋帶式混合方法及性能試驗(yàn)研究[D]. 鎮(zhèn)江: 江蘇大學(xué), 2016.

    [9]? 邢潔潔, 張? 銳, 吳? 鵬, 等. 海南熱區(qū)磚紅壤顆粒離散元仿真模型參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2020, 36(5): 158-166.

    [10]????? Ucgul M, Fielke J M, Saunders C. Three dimensional discrete element modeling DEM of tillage accounting for soil cohesion and adhesion[J]. Biosystems Engineering, 2015(129): 298-306.

    [11]???? 羅? 帥, 袁巧霞, Gouda Shaban, 等. 基于JKR粘結(jié)模型的蚯蚓糞基質(zhì)離散元法參數(shù)標(biāo)定[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2018, 49(4): 343-350.

    [12]????? Bangura K, Gong H, Deng R L, et al. Simulation analysis of fertilizer discharge process using the Discrete Element Method (DEM)[J]. PLoS One, 2020, 15(7): e0235872 .

    [13]????? 鄭智旗. 玉米秸稈撿拾粉碎溝埋還田機(jī)研究[D]. 北京: 中國(guó)農(nóng)業(yè)大學(xué), 2017.

    [14]????? Ghodki B M, Patel M, Namdeo R, et al. Calibration of discrete element model parameters: soybeans[J]. Computational Particle Mechanics, 2018, 6(1): 3-10.

    [15]????? 彭才望, 許道軍, 賀? 喜, 等. 黑水虻處理的豬糞有機(jī)肥離散元仿真模型參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2020, 36(17): 212-218.

    [16]????? 張? 珂, 于天行, 于文達(dá), 等. 基于JKR黏結(jié)模型的混凝土離散元參數(shù)標(biāo)定[J]. 混凝土, 2020(8): 46-50, 55.

    [17]????? 文恩楊, 李玉華, 牛子孺, 等. 蒜種顆粒離散元模型參數(shù)標(biāo)定[J]. 農(nóng)機(jī)化研究, 2021, 43(5): 160-167.

    [18]????? 武? 濤, 黃偉鳳, 陳學(xué)深, 等. 考慮顆粒間黏結(jié)力的黏性土壤離散元模型參數(shù)標(biāo)定[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2017, 38(3): 93-98.

    [19]???? Cundall P A, Strack O D L. A discrete numerical model for granular assembles[J]. Geotechnique, 1979, 29(1): 47-65.

    [20]????? 團(tuán)體標(biāo)準(zhǔn)T/CAMA 12—2019無(wú)土栽培椰糠[J]. 農(nóng)業(yè)工程技術(shù), 2019, 39(31): 69-73.

    [21]????? 吳? 鵬, 張喜瑞, 李? 粵, 等. 基于離散元的磚紅土壤模型堆積特性研究[J]. 農(nóng)機(jī)化研究, 2020, 42(9): 145-150.

    [22]????? 劉羽平, 張? 拓, 劉? 妤. 稻谷顆粒模型離散元接觸參數(shù)標(biāo)定與試驗(yàn)[J]. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào), 2019, 21(11): 70-76.

    [23]????? 張? 銳, 韓佃雷, 吉巧麗, 等. 離散元模擬中沙土參數(shù)標(biāo)定方法研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2017, 48(3): 49-56.

    [24]????? Thomas S, Hassanizadeh S M, Hamed A, et al. The effect of particle shape on porosity of swelling granular materials: Discrete element method and the multi-sphere approximation[J]. Powder Technology, 2020, 360(15): 1295-1304.

    [25]????? 馬彥華, 宋春東, 宣傳忠, 等. 苜蓿秸稈壓縮仿真離散元模型參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2020, 36(11): 22-30.

    [26]????? 廖宜濤, 廖慶喜, 周? 宇, 等. 飼料油菜薹期收獲莖稈破碎離散元仿真參數(shù)標(biāo)定[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2020, 51(6): 73-82.

    [27]???? 吳文江, 郭? 斌, 高占鳳, 等. 基于DEM的沙土顆粒建模及參數(shù)標(biāo)定研究[J]. 中國(guó)農(nóng)機(jī)化學(xué)報(bào), 2019, 40(8): 182-187.

    [28]????? 向? 偉, 吳明亮, 呂江南, 等. 基于堆積試驗(yàn)的黏壤土仿真物理參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2019, 35(12): 116-123.

    [29]????? Mehrdad P, Colin H, Mojtaba G, et al. Effect of particle shape on flow in discrete element method simulation of a rotary batch seed coater[J]. Powder Technology, 2016, 296: 29-36.

    [30]????? Józef H, Marek M. Parameters and contact models for DEM simulations of agricultural granular materials: A review[J]. Biosystems Engineering, 2016, 147: 206-225.

    責(zé)任編輯:崔麗虹

    猜你喜歡
    椰糠
    椰糠-黏土植生基材崩解特性試驗(yàn)1)
    商品巖棉條和椰糠條對(duì)黃瓜生長(zhǎng)、光合、產(chǎn)量和品質(zhì)的影響
    設(shè)施栽培條件下不同椰糠栽培模式對(duì)番茄生長(zhǎng)發(fā)育的影響
    設(shè)施栽培條件下不同椰糠栽培模式對(duì)番茄生長(zhǎng)發(fā)育的影響
    浸泡時(shí)長(zhǎng)對(duì)椰糠基質(zhì)理化性狀的影響
    武山縣設(shè)施黃瓜“椰糠條水肥一體化”栽培技術(shù)探討
    灌溉液不同pH值對(duì)椰糠種植番茄的影響
    新型椰糠基質(zhì)與泥炭基質(zhì)栽培小白菜效果研究
    不同配比的椰糠與珍珠巖基質(zhì)對(duì)番茄幼苗生長(zhǎng)和育苗效果的影響
    溫室禮品西瓜椰糠無(wú)土栽培試驗(yàn)總結(jié)
    亚洲经典国产精华液单| 久久午夜亚洲精品久久| 国产一级毛片七仙女欲春2| 中文在线观看免费www的网站| 在现免费观看毛片| 国产精品国产三级国产av玫瑰| 又粗又爽又猛毛片免费看| 国产精品98久久久久久宅男小说| 国产精品国产高清国产av| 亚洲国产日韩欧美精品在线观看| 乱码一卡2卡4卡精品| 干丝袜人妻中文字幕| 亚洲人与动物交配视频| 亚洲五月天丁香| 国产精品日韩av在线免费观看| 在线观看一区二区三区| 国产高清不卡午夜福利| 成人永久免费在线观看视频| 国产大屁股一区二区在线视频| 午夜激情欧美在线| 亚洲欧美日韩卡通动漫| 久久久久久久午夜电影| 日韩欧美一区二区三区在线观看| 久久人人精品亚洲av| 动漫黄色视频在线观看| 欧美精品啪啪一区二区三区| 午夜日韩欧美国产| 亚洲第一电影网av| 级片在线观看| 99久国产av精品| 三级毛片av免费| 亚洲最大成人手机在线| а√天堂www在线а√下载| 国产精品爽爽va在线观看网站| 亚洲精品粉嫩美女一区| 国产亚洲91精品色在线| 香蕉av资源在线| 午夜福利18| 人人妻人人澡欧美一区二区| 国产三级在线视频| 日韩欧美精品v在线| 国产乱人视频| 亚洲一区二区三区色噜噜| 亚洲人成网站在线播放欧美日韩| 99久久九九国产精品国产免费| 婷婷亚洲欧美| 欧美激情久久久久久爽电影| 伦精品一区二区三区| 亚洲成av人片在线播放无| 午夜免费男女啪啪视频观看 | 久久久久久国产a免费观看| 久久精品夜夜夜夜夜久久蜜豆| 久久热精品热| 久久久久久九九精品二区国产| 22中文网久久字幕| 嫩草影视91久久| 国产精品久久久久久亚洲av鲁大| 亚洲第一区二区三区不卡| 亚洲最大成人中文| 99久久九九国产精品国产免费| 尾随美女入室| 午夜精品一区二区三区免费看| 色av中文字幕| 亚洲精品456在线播放app | 熟女电影av网| 九九在线视频观看精品| 日本黄色视频三级网站网址| 精品一区二区三区人妻视频| 国产一级毛片七仙女欲春2| 亚洲一区二区三区色噜噜| 悠悠久久av| 欧美一区二区精品小视频在线| 成人国产一区最新在线观看| 色综合婷婷激情| 干丝袜人妻中文字幕| 精品人妻熟女av久视频| 日韩亚洲欧美综合| 成人性生交大片免费视频hd| 日本撒尿小便嘘嘘汇集6| aaaaa片日本免费| 97超级碰碰碰精品色视频在线观看| 欧美3d第一页| 欧美zozozo另类| 国产精品一及| 婷婷六月久久综合丁香| 久久亚洲精品不卡| 色综合亚洲欧美另类图片| 99热只有精品国产| 久久久久久久久大av| 高清日韩中文字幕在线| 国产精品免费一区二区三区在线| 国产精品免费一区二区三区在线| 国产av不卡久久| 国产高清激情床上av| 少妇丰满av| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久精品电影| 精品一区二区免费观看| 久久精品国产亚洲av香蕉五月| 久久久久久久久久黄片| 露出奶头的视频| 嫩草影院新地址| 最近中文字幕高清免费大全6 | 精品一区二区三区视频在线| 国内精品久久久久久久电影| 亚洲美女搞黄在线观看 | 亚洲在线自拍视频| 国产精品野战在线观看| 成人永久免费在线观看视频| 黄色一级大片看看| 久久香蕉精品热| 男女啪啪激烈高潮av片| 无遮挡黄片免费观看| 在线观看一区二区三区| 丰满的人妻完整版| 欧美性猛交黑人性爽| 一区二区三区四区激情视频 | 日本免费一区二区三区高清不卡| 午夜福利18| 九九久久精品国产亚洲av麻豆| 热99在线观看视频| 亚洲自偷自拍三级| 18禁在线播放成人免费| 午夜福利在线观看吧| 少妇人妻精品综合一区二区 | 午夜福利在线观看免费完整高清在 | av天堂中文字幕网| 悠悠久久av| 亚洲精品日韩av片在线观看| 精品一区二区免费观看| 精品午夜福利在线看| 国产精品一及| 久久婷婷人人爽人人干人人爱| 欧美+亚洲+日韩+国产| 国内揄拍国产精品人妻在线| 99在线人妻在线中文字幕| 日日干狠狠操夜夜爽| 99riav亚洲国产免费| 国产单亲对白刺激| 日韩一区二区视频免费看| 黄色日韩在线| 中文字幕av成人在线电影| 亚洲成人中文字幕在线播放| 国产精品av视频在线免费观看| 亚洲精品粉嫩美女一区| 少妇裸体淫交视频免费看高清| 能在线免费观看的黄片| av黄色大香蕉| 欧美+日韩+精品| 日韩一区二区视频免费看| 女人十人毛片免费观看3o分钟| 日本免费a在线| 久久久久久久午夜电影| 日本a在线网址| 少妇裸体淫交视频免费看高清| 99国产精品一区二区蜜桃av| 久久久久国内视频| 久久精品国产亚洲av天美| av在线老鸭窝| 国产午夜精品论理片| 国产午夜福利久久久久久| 男人舔奶头视频| 国产蜜桃级精品一区二区三区| 男女边吃奶边做爰视频| 亚洲精品久久国产高清桃花| 九九爱精品视频在线观看| 日韩国内少妇激情av| 悠悠久久av| 久久99热这里只有精品18| 精品不卡国产一区二区三区| 国产极品精品免费视频能看的| 最好的美女福利视频网| 国产91精品成人一区二区三区| 男女之事视频高清在线观看| 麻豆成人午夜福利视频| 99久久成人亚洲精品观看| 亚洲欧美日韩东京热| 身体一侧抽搐| 午夜免费成人在线视频| 午夜福利在线在线| 亚洲精品成人久久久久久| 十八禁网站免费在线| 夜夜爽天天搞| 中文字幕精品亚洲无线码一区| 日本一本二区三区精品| 欧美日韩瑟瑟在线播放| 久久热精品热| 亚洲av免费在线观看| 午夜亚洲福利在线播放| 两人在一起打扑克的视频| 在线国产一区二区在线| 欧美一区二区精品小视频在线| 丰满人妻一区二区三区视频av| 深爱激情五月婷婷| 非洲黑人性xxxx精品又粗又长| 真实男女啪啪啪动态图| 三级男女做爰猛烈吃奶摸视频| 草草在线视频免费看| 国产高清有码在线观看视频| 国产成人av教育| 中文字幕久久专区| 亚洲成人免费电影在线观看| 无遮挡黄片免费观看| 精品久久国产蜜桃| 国产精华一区二区三区| 久久久久久久久久久丰满 | 国产av在哪里看| 亚洲精品粉嫩美女一区| 成人无遮挡网站| 精品日产1卡2卡| 欧美日本视频| 大型黄色视频在线免费观看| 亚洲av成人av| 久久午夜福利片| 免费看日本二区| 在线观看舔阴道视频| 国产av在哪里看| 国产精品自产拍在线观看55亚洲| av女优亚洲男人天堂| 国产久久久一区二区三区| 成年版毛片免费区| 久久婷婷人人爽人人干人人爱| 男女之事视频高清在线观看| 51国产日韩欧美| 午夜视频国产福利| 91午夜精品亚洲一区二区三区 | 欧洲精品卡2卡3卡4卡5卡区| 日韩大尺度精品在线看网址| 看免费成人av毛片| 国产三级中文精品| 长腿黑丝高跟| 国内精品宾馆在线| 国产在视频线在精品| 国产乱人伦免费视频| 日本五十路高清| 久久久久久国产a免费观看| 在现免费观看毛片| 亚洲avbb在线观看| 国产精品乱码一区二三区的特点| 日本五十路高清| 久久久久久大精品| 可以在线观看的亚洲视频| 亚洲在线观看片| 欧美日韩综合久久久久久 | 欧美日韩中文字幕国产精品一区二区三区| 在现免费观看毛片| 精品久久久久久久久久免费视频| 网址你懂的国产日韩在线| 成人三级黄色视频| 日日干狠狠操夜夜爽| 国产69精品久久久久777片| 国产午夜精品论理片| 亚洲精品在线观看二区| 综合色av麻豆| 人妻少妇偷人精品九色| 日韩一区二区视频免费看| 大型黄色视频在线免费观看| АⅤ资源中文在线天堂| 一本久久中文字幕| 国产高清三级在线| 国产亚洲精品久久久久久毛片| 香蕉av资源在线| 国产精品久久久久久av不卡| 亚洲久久久久久中文字幕| 搞女人的毛片| 少妇人妻一区二区三区视频| 嫩草影院新地址| 老司机福利观看| 少妇高潮的动态图| 精品久久久久久久久亚洲 | 欧美一级a爱片免费观看看| 亚洲国产高清在线一区二区三| 亚洲一级一片aⅴ在线观看| 亚洲av免费高清在线观看| 国产 一区 欧美 日韩| 精品不卡国产一区二区三区| 午夜a级毛片| 亚洲七黄色美女视频| xxxwww97欧美| 啦啦啦观看免费观看视频高清| 日本成人三级电影网站| 亚洲综合色惰| 一级毛片久久久久久久久女| 午夜福利成人在线免费观看| 国产乱人伦免费视频| 国产高清激情床上av| 最后的刺客免费高清国语| 久久人妻av系列| 国产爱豆传媒在线观看| 国产精品日韩av在线免费观看| 两人在一起打扑克的视频| 国产探花极品一区二区| 午夜福利在线观看免费完整高清在 | АⅤ资源中文在线天堂| 久久久色成人| 免费一级毛片在线播放高清视频| 精品人妻偷拍中文字幕| 国产高清有码在线观看视频| 日韩欧美国产在线观看| 免费观看在线日韩| 欧美黑人欧美精品刺激| 18禁黄网站禁片午夜丰满| 天天一区二区日本电影三级| 久久久久久伊人网av| 日韩欧美精品免费久久| 成人综合一区亚洲| 久久热精品热| 男女之事视频高清在线观看| 如何舔出高潮| 熟女电影av网| a级毛片免费高清观看在线播放| 国产单亲对白刺激| 夜夜看夜夜爽夜夜摸| 亚洲国产欧美人成| 国产伦精品一区二区三区四那| 日韩欧美三级三区| 久久国产精品人妻蜜桃| www.色视频.com| 欧美成人免费av一区二区三区| 成人国产一区最新在线观看| 可以在线观看的亚洲视频| 999久久久精品免费观看国产| 亚洲精品久久国产高清桃花| 久久香蕉精品热| 色精品久久人妻99蜜桃| 一级黄色大片毛片| 俺也久久电影网| 亚洲av五月六月丁香网| 欧美成人a在线观看| 成人永久免费在线观看视频| 九九久久精品国产亚洲av麻豆| 日本精品一区二区三区蜜桃| 亚洲av美国av| 亚洲人成伊人成综合网2020| 国产精品99久久久久久久久| 热99在线观看视频| 国产成人aa在线观看| 亚洲经典国产精华液单| 亚洲一级一片aⅴ在线观看| www.www免费av| 午夜精品一区二区三区免费看| 国产一级毛片七仙女欲春2| 久久天躁狠狠躁夜夜2o2o| 婷婷亚洲欧美| 精品久久久久久久末码| 亚洲人成网站在线播放欧美日韩| 禁无遮挡网站| 三级国产精品欧美在线观看| 99精品久久久久人妻精品| 麻豆久久精品国产亚洲av| 给我免费播放毛片高清在线观看| 身体一侧抽搐| 可以在线观看的亚洲视频| 日韩欧美 国产精品| 国产欧美日韩一区二区精品| 国产免费一级a男人的天堂| 国产亚洲av嫩草精品影院| 亚洲中文字幕日韩| 久久午夜亚洲精品久久| 成人毛片a级毛片在线播放| 麻豆成人午夜福利视频| av.在线天堂| 99久久精品国产国产毛片| 全区人妻精品视频| 欧美色视频一区免费| 三级国产精品欧美在线观看| 可以在线观看的亚洲视频| 久久热精品热| 99热这里只有是精品在线观看| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 日韩欧美国产一区二区入口| 国产69精品久久久久777片| 久久人人精品亚洲av| 国产69精品久久久久777片| 欧美日本亚洲视频在线播放| 欧美成人免费av一区二区三区| eeuss影院久久| 亚洲av电影不卡..在线观看| 他把我摸到了高潮在线观看| 成人美女网站在线观看视频| 女生性感内裤真人,穿戴方法视频| 日日摸夜夜添夜夜添小说| 男女边吃奶边做爰视频| 中文资源天堂在线| 一区二区三区四区激情视频 | 九九爱精品视频在线观看| 男女做爰动态图高潮gif福利片| 国产精品永久免费网站| 97超视频在线观看视频| 啦啦啦观看免费观看视频高清| 九色成人免费人妻av| 看片在线看免费视频| or卡值多少钱| 欧美精品啪啪一区二区三区| 色精品久久人妻99蜜桃| 精品一区二区三区视频在线观看免费| 午夜爱爱视频在线播放| 国产三级中文精品| 中出人妻视频一区二区| 久久久成人免费电影| 国产av一区在线观看免费| 国产私拍福利视频在线观看| 成人一区二区视频在线观看| 舔av片在线| 一个人看视频在线观看www免费| 91在线精品国自产拍蜜月| 黄色一级大片看看| 久久久久九九精品影院| 老熟妇仑乱视频hdxx| 中文字幕精品亚洲无线码一区| 亚洲精品一区av在线观看| 中文字幕精品亚洲无线码一区| 人人妻人人澡欧美一区二区| 中文字幕免费在线视频6| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩高清在线视频| 九九热线精品视视频播放| 亚洲性久久影院| av视频在线观看入口| .国产精品久久| 美女cb高潮喷水在线观看| 国产女主播在线喷水免费视频网站 | 麻豆av噜噜一区二区三区| 精品人妻视频免费看| 99热这里只有精品一区| 国产精品,欧美在线| 中文字幕久久专区| 99九九线精品视频在线观看视频| 亚洲av五月六月丁香网| 国产亚洲精品久久久com| 国产精品综合久久久久久久免费| 亚洲成a人片在线一区二区| 精华霜和精华液先用哪个| 直男gayav资源| 久久中文看片网| 欧美潮喷喷水| 国产黄色小视频在线观看| 亚洲av五月六月丁香网| 亚洲最大成人手机在线| 99国产极品粉嫩在线观看| 成年女人毛片免费观看观看9| 国产主播在线观看一区二区| 在线播放国产精品三级| 日本 av在线| 国产又黄又爽又无遮挡在线| 麻豆国产97在线/欧美| 久久久久免费精品人妻一区二区| 丝袜美腿在线中文| 国产一级毛片七仙女欲春2| 狂野欧美白嫩少妇大欣赏| 亚洲美女黄片视频| 又爽又黄a免费视频| 欧美精品国产亚洲| av在线亚洲专区| 免费搜索国产男女视频| av在线蜜桃| xxxwww97欧美| 黄色女人牲交| av专区在线播放| 一级毛片久久久久久久久女| av福利片在线观看| 可以在线观看的亚洲视频| 日韩欧美精品v在线| 国产精品永久免费网站| 色精品久久人妻99蜜桃| ponron亚洲| 一级av片app| 一进一出好大好爽视频| 国产在线精品亚洲第一网站| 欧美性猛交黑人性爽| 日本熟妇午夜| 日本欧美国产在线视频| 日本与韩国留学比较| 亚洲不卡免费看| 国产欧美日韩精品亚洲av| 日日摸夜夜添夜夜添av毛片 | 亚洲一区二区三区色噜噜| 亚洲人成伊人成综合网2020| 搡女人真爽免费视频火全软件 | 色综合婷婷激情| 成人午夜高清在线视频| 999久久久精品免费观看国产| 韩国av在线不卡| 在线观看一区二区三区| 国产探花在线观看一区二区| 观看美女的网站| 国产精品亚洲一级av第二区| 中文字幕精品亚洲无线码一区| 欧美bdsm另类| 中文字幕久久专区| 搡老妇女老女人老熟妇| 精品一区二区三区视频在线观看免费| 夜夜夜夜夜久久久久| 久久精品国产亚洲av香蕉五月| 成人一区二区视频在线观看| av视频在线观看入口| 不卡一级毛片| 久久草成人影院| 88av欧美| 老司机深夜福利视频在线观看| 久久久久九九精品影院| 免费看a级黄色片| 尤物成人国产欧美一区二区三区| 亚洲无线在线观看| 国产精品1区2区在线观看.| 国产精品一区www在线观看 | 国产探花在线观看一区二区| 最近中文字幕高清免费大全6 | 最近最新免费中文字幕在线| 大型黄色视频在线免费观看| 内地一区二区视频在线| 琪琪午夜伦伦电影理论片6080| 又黄又爽又免费观看的视频| 1024手机看黄色片| 3wmmmm亚洲av在线观看| 啦啦啦韩国在线观看视频| а√天堂www在线а√下载| 日本-黄色视频高清免费观看| 成年免费大片在线观看| 免费电影在线观看免费观看| 欧美日韩国产亚洲二区| 春色校园在线视频观看| 国产精品福利在线免费观看| 国产高清三级在线| 校园人妻丝袜中文字幕| 国产精品一区二区性色av| 久久草成人影院| 十八禁网站免费在线| 亚洲欧美日韩东京热| 久久国内精品自在自线图片| 精品人妻视频免费看| 淫妇啪啪啪对白视频| 免费大片18禁| 亚洲一区二区三区色噜噜| 91精品国产九色| 久久人人爽人人爽人人片va| 欧美日韩瑟瑟在线播放| 午夜爱爱视频在线播放| 婷婷六月久久综合丁香| 亚洲精华国产精华液的使用体验 | 丝袜美腿在线中文| 精品人妻偷拍中文字幕| 日本三级黄在线观看| 综合色av麻豆| 国产视频一区二区在线看| 欧美色欧美亚洲另类二区| 久久中文看片网| 欧美人与善性xxx| 久久久色成人| 在线观看一区二区三区| 成人精品一区二区免费| 最后的刺客免费高清国语| 国产三级在线视频| 国产精品一区二区免费欧美| 在线观看美女被高潮喷水网站| 亚洲中文字幕一区二区三区有码在线看| 国产精品伦人一区二区| 欧美精品国产亚洲| 亚洲av五月六月丁香网| 精品无人区乱码1区二区| 高清日韩中文字幕在线| 搞女人的毛片| 在线免费十八禁| 久久这里只有精品中国| 亚洲国产日韩欧美精品在线观看| 国产亚洲av嫩草精品影院| 美女免费视频网站| 国产伦精品一区二区三区四那| 欧美成人a在线观看| av天堂在线播放| 99热这里只有是精品50| 在线观看美女被高潮喷水网站| 身体一侧抽搐| 亚洲美女搞黄在线观看 | 性欧美人与动物交配| 亚洲av成人av| 国语自产精品视频在线第100页| 免费在线观看成人毛片| 国产精品国产三级国产av玫瑰| 免费看光身美女| 国产美女午夜福利| 国产私拍福利视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 1024手机看黄色片| 男女做爰动态图高潮gif福利片| 亚洲国产日韩欧美精品在线观看| 国内精品宾馆在线| 日本 av在线| av在线亚洲专区| 麻豆成人午夜福利视频| 黄色女人牲交| 欧美绝顶高潮抽搐喷水| 18禁黄网站禁片午夜丰满| 99久久中文字幕三级久久日本| 两性午夜刺激爽爽歪歪视频在线观看| 乱系列少妇在线播放| 久久精品国产清高在天天线| 免费看光身美女| 在线免费观看不下载黄p国产 | 欧美日韩国产亚洲二区| 国产真实乱freesex| 香蕉av资源在线| 日韩,欧美,国产一区二区三区 | 草草在线视频免费看| 国产伦在线观看视频一区| 国产人妻一区二区三区在| 久久精品国产清高在天天线| 亚洲欧美日韩东京热| 很黄的视频免费| 亚洲av成人av| 日本成人三级电影网站| 精品午夜福利视频在线观看一区| or卡值多少钱| 悠悠久久av|