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

    不同含水率黏重黑土與觸土部件互作的離散元仿真參數(shù)標(biāo)定

    2019-05-11 07:02:30李俊偉王虎彪毛春昱馬云海
    農(nóng)業(yè)工程學(xué)報 2019年6期
    關(guān)鍵詞:黑土摩擦系數(shù)標(biāo)定

    李俊偉,佟 金,胡 斌,王虎彪,毛春昱,5,馬云?!?/p>

    ?

    不同含水率黏重黑土與觸土部件互作的離散元仿真參數(shù)標(biāo)定

    李俊偉1,2,3,4,佟 金1,3,胡 斌2,4,王虎彪1,3,毛春昱1,3,5,馬云海1,3※

    (1. 吉林大學(xué)工程仿生教育部重點實驗室,長春 130022;2. 石河子大學(xué)機(jī)械電氣工程學(xué)院,石河子 832003;3. 吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院,長春 130022; 4. 石河子大學(xué)農(nóng)業(yè)部西北農(nóng)業(yè)裝備重點實驗室,石河子 832003; 5. 吉林工程師范技術(shù)學(xué)院,長春 130022)

    為了獲得可用于東北地區(qū)黏重黑土與觸土部件相互作用的離散元仿真模擬參數(shù),該文利用EDEM中Hertz-Mindlin with JKR Cohesion接觸模型對不同含水率的東北地區(qū)黏重黑土進(jìn)行相關(guān)參數(shù)標(biāo)定,針對含水率在10%~20%的實際作業(yè)環(huán)境,分別配置含水率為12.46%±1%和17.15%±1%的2種黏重黑土,以土壤顆粒間的滾動摩擦系數(shù)、恢復(fù)系數(shù)、表面能參數(shù)及靜摩擦系數(shù)為標(biāo)定對象,并以土壤顆粒的仿真堆積角為響應(yīng)值,基于Box-Behnken的響應(yīng)面優(yōu)化方法得到堆積角回歸模型,并對回歸模型進(jìn)行尋優(yōu),得到2種含水率的模型參數(shù)優(yōu)化解,并給出了模型參數(shù)范圍。測定了4種含水率下黏重黑土對3種觸土部件材料(65Mn、UHMW-PE和PTFE)的靜摩擦系數(shù),并以此為基礎(chǔ)分別對65Mn(典型鐵基材料)和PTFE(典型低表面能材料)板進(jìn)行斜面試驗,以含水率為17.15%±1%的黏重黑土為試驗對象,分別搭建斜面物理試驗平臺和仿真模型,以土壤顆粒與觸土部件材料之間的滾動摩擦系數(shù)、恢復(fù)系數(shù)、表面能參數(shù)及靜摩擦系數(shù)為標(biāo)定對象,以仿真得到的土球在65Mn和PTFE板上的滾動距離為響應(yīng)值,基于響應(yīng)面優(yōu)化法得到滾動距離回歸模型,以實測的滾動距離為目標(biāo)對回歸模型進(jìn)行尋優(yōu),得到黏重黑土對2類典型觸土部件材料接觸模型參數(shù)的優(yōu)化解。研究結(jié)果表明,標(biāo)定優(yōu)化后的土壤模型能夠近似代替真實的東北地區(qū)黏重黑土進(jìn)行仿真,可利用標(biāo)定后的參數(shù)進(jìn)行黏重黑土與觸土部件間的離散元仿真,可為東北黏重黑土作業(yè)條件下的農(nóng)業(yè)機(jī)械觸土部件仿生減阻設(shè)計與優(yōu)化提供基礎(chǔ)數(shù)據(jù)。

    土壤;含水率;離散元法;黏重黑土;仿真參數(shù)標(biāo)定

    0 引 言

    觸土部件與土壤接觸或在切割土壤過程中所涉及到的相互作用力學(xué)問題一直是該領(lǐng)域研究的難點與熱點。由于土壤的復(fù)雜性與多樣性,目前尚未有一種準(zhǔn)確的數(shù)學(xué)模型能替代所有類型土壤的本構(gòu)模型。對于現(xiàn)代觸土部件的設(shè)計及其優(yōu)化,通常需要對觸土部件與土壤相互作用的過程進(jìn)行仿真分析,而現(xiàn)有的有限元土壤模型由于其采用的是連續(xù)介質(zhì)的表達(dá)形式而無法模擬土壤顆粒之間運動及其對耕作部件相互作用的動態(tài)過程[1-2]。離散元素建模(discrete element modeling)是一種用于對粒狀材料的力學(xué)行為進(jìn)行建模的數(shù)值方法[3],該方法將散粒體簡化成具有質(zhì)量和形狀的顆粒集合,解決各個顆粒體之間及其與邊界間的相互作用,從而揭示顆粒與邊界之間的物理作用性質(zhì)[4]。近幾年來,離散元法在工程中得到越來越廣泛的應(yīng)用,越來越多的學(xué)者運用離散元法進(jìn)行散粒物料的流動性、固體破碎及機(jī)具-土壤相互作用方面的研究。在農(nóng)業(yè)工程領(lǐng)域,離散元法的應(yīng)用主要集中在以散粒體為主的農(nóng)業(yè)物料學(xué),土壤與農(nóng)機(jī)部件的相互作用等方面的研究[5-7]。通常在利用離散元法進(jìn)行散粒體研究之前,需對散粒體的物料特性參數(shù)進(jìn)行標(biāo)定。除了模擬土壤顆粒之外,基于離散元法的接觸模型參數(shù)標(biāo)定還被應(yīng)用于比如小麥、水稻、馬鈴薯、大豆等其他散粒體[8-13]。

    土壤與觸土部件之間相互作用的精準(zhǔn)模型對設(shè)計與優(yōu)化觸土部件具有重要作用[14-15]。然而,觸土部件在土壤中進(jìn)行耕作作業(yè)時,土壤-機(jī)具之間的相互作用受土壤性質(zhì)及土壤的動態(tài)作用直接影響[16-17],對不同類型的土壤,其影響效應(yīng)差異顯著,因此,在運用離散元法對觸土部件與土壤的相互作用進(jìn)行仿真前,需先標(biāo)定土壤顆粒與觸土部件之間的相關(guān)仿真參數(shù)。

    目前,國內(nèi)外學(xué)者主要針對EDEM 軟件中以Hertze-Mindlin等模型為主的內(nèi)嵌模型對土壤顆粒接觸參數(shù)進(jìn)行研究,Mustafa等[1,3,14-15]利用Hertze-Mindlin及Hysteretic Spring接觸模型,研究分析了土壤在有黏結(jié)力和無黏結(jié)力的情況下土壤顆粒的塑性形變等問題。丁啟朔等[18]利用Hertz-Mindlin with Bonding模型建立了黏性水稻土的深松耕作離散元模型。雖然Hertz-Mindlin with Bonding接觸模型可以用來模擬粘結(jié)顆粒,但是該模型是用一個有限大小的“膠粘劑”粘結(jié),適用于模擬混凝土和巖石等硬粒介質(zhì)。張銳等[19]提出了一種標(biāo)定沙土顆粒相互作用參數(shù)的方法,該方法使用標(biāo)準(zhǔn)球和非標(biāo)準(zhǔn)球?qū)ι惩令w粒的靜摩擦系數(shù)、恢復(fù)系數(shù)等關(guān)鍵參數(shù)進(jìn)行標(biāo)定,為沙土顆粒的離散元接觸模型參數(shù)標(biāo)定提供了有價值的參考。王憲良等[4]基于代理模型的方法,利用The Edinburgh Elasto-Plastic Cohesion Model(ECM)彈塑性接觸模型表示土壤顆粒接觸模型,根據(jù)試驗測定堆積角及剪切試驗得到的參數(shù)建立離散元土壤模型,仿真分析模型參數(shù)的敏感性,最后以敏感性參數(shù)為變量,以試驗值為目標(biāo)構(gòu)建代理模型并進(jìn)行參數(shù)優(yōu)化,ECM模型是一種將土壤應(yīng)變滯后性、內(nèi)聚力及范德華力等考慮在內(nèi)的非線性模型,適用于研究土壤、有機(jī)材料等的壓縮問題。石林榕等[20-21]整合延遲彈性模型(hysteretic spring contact model,HSCM)和線性內(nèi)聚力模型(liner cohesion model,LCM)優(yōu)勢,建立了西北旱區(qū)農(nóng)田土壤模型,對土壤顆粒之間的靜摩擦系數(shù)、滾動摩擦系數(shù)和內(nèi)聚強度等關(guān)鍵參數(shù)進(jìn)行離散元標(biāo)定。

    上述EDEM軟件中的接觸模型大多適用于黏性弱含水率較低的土壤顆粒,如沙土、旱土或者巖石等,而對于東北地區(qū)的黏重黑土,其宏觀表現(xiàn)出的黏性及塑性特征遠(yuǎn)大于土壤顆粒之間的彈性特征,土壤顆粒之間粘結(jié)力較大,上述接觸模型并不適用。

    武濤等[22]基于土壤堆積角物理試驗結(jié)果,采用考慮顆粒間粘結(jié)力的“Hertz-Mindlin with JKR Cohesion”接觸模型進(jìn)行土壤堆積角仿真試驗,對考慮粘結(jié)力的土壤顆粒之間的關(guān)鍵參數(shù)如JKR表面能參數(shù)、恢復(fù)系數(shù)、靜摩擦系數(shù)與滾動摩擦系數(shù)進(jìn)行接觸模型的參數(shù)標(biāo)定,試驗驗證了JKR模型可以用來模擬含水率較高的黏濕南方土壤。課題組在實際調(diào)研中發(fā)現(xiàn),對于東北地區(qū)的黑黏土,土壤孔隙率低,黏粒多,腐殖質(zhì)含量高,土壤性質(zhì)復(fù)雜,含水率較高,在作物收獲期間的含水率一般在10%~20%之間,特別是下雨后的土壤,其含水率更高,當(dāng)農(nóng)機(jī)耕作部件進(jìn)行作業(yè)時,觸土部件的土壤粘附現(xiàn)象嚴(yán)重,耕作阻力較大,嚴(yán)重影響生產(chǎn)效能。

    本文前期研究發(fā)現(xiàn),由于粘附能力很強,當(dāng)含水率超過20%后,黏重黑土近于流變體性質(zhì),其黏性及塑性性質(zhì)很突出,在實際試驗過程中發(fā)現(xiàn),當(dāng)含水率超過12%后,土壤顆粒之間就會發(fā)生較為嚴(yán)重的團(tuán)聚現(xiàn)象,由此可見,對于黏重黑土,其顆粒之間的粘結(jié)力是必須要重點考慮的因素。因此,本文根據(jù)東北地區(qū)的黏重黑土特性,配置2種含水率的黏重黑土,在EDEM2.7軟件平臺上,基于顆粒間粘結(jié)力的“Hertz-Mindlin with JKR Cohesion”接觸模型,采用物理試驗與仿真試驗相結(jié)合的方法,分別對含水率為12.46%±1%和17.15%±1%的黏重黑土進(jìn)行堆積角仿真試驗,基于Box-Behnken的響應(yīng)面優(yōu)化方法,對不同含水率黏重黑土顆粒之間的接觸模型仿真參數(shù)進(jìn)行標(biāo)定,對4種含水率的黏重黑土與3種觸土部件材料(典型鐵基類觸土部件材料65Mn及2種低表面能材料超高分子量聚乙烯(UHMW-PE)和聚四氟乙烯(PTFE))的靜摩擦系數(shù)進(jìn)行試驗測定,并在此基礎(chǔ)上利用斜面試驗方法,分別對較高含水率下(含水率17.15%±1%)的黏重黑土與2類典型觸土部件材料(65Mn鋼板和PTFE板)的接觸模型參數(shù)進(jìn)行標(biāo)定,以期為進(jìn)一步研究東北地區(qū)黏重黑土或者含水率較高的黏濕土壤與觸土部件之間的相互作用提供離散元仿真所需的基礎(chǔ)參數(shù),并為觸土部件減粘降阻仿生設(shè)計提供必要的仿真參數(shù)。

    1 土壤顆粒間的接觸模型參數(shù)標(biāo)定

    1.1 接觸模型的選取

    式中pullout為將2個顆粒分開所需的力,;為接觸角,(°);為顆粒半徑,mm。

    1.2 模型參數(shù)類型及標(biāo)定

    對于含水率較高的黏重黑土,顆粒之間的黏結(jié)力或者顆粒與觸土部件材料之間的粘結(jié)力是重要考慮對象,在EDEM JKR模型中用表面能來表征顆粒之間或顆粒與觸土部件材料之間的黏結(jié)力。

    1.2.1 試驗材料及其本征參數(shù)

    試驗用土壤樣品采自吉林農(nóng)業(yè)大學(xué)試驗田,土壤類型為典型東北黏重黑土,土壤堆積密度為1.07 kg/m3,為了更準(zhǔn)確地反映實際作業(yè)土壤狀況,考慮到雨后土壤含水率實測值與連續(xù)晴天下的土壤含水率差異,以2種含水率的黏重黑土為對象進(jìn)行離散元土壤顆粒接觸模型的參數(shù)標(biāo)定,因此本文采用配置的含水率為12.46%±1%和17.15%±1%的黏重黑土為對象進(jìn)行堆積角物理試驗。土壤的本征參數(shù)通過圖1所示的SLB-1型應(yīng)力應(yīng)變控制式三軸試驗儀測得,結(jié)果如表1所示。

    1.2.2 土壤堆積角物理試驗

    圖1 土壤參數(shù)測試

    圖2 土壤堆積角測試

    1.2.3 土壤堆積角仿真試驗?zāi)P偷慕?/p>

    參考文獻(xiàn)[22]的方法,在EDEM 2.7軟件中建立堆積角仿真模型。漏斗頂部直徑為150 mm,底部直徑為30 mm,高為200 mm,接料底板為400 mm×400 mm的方形板。漏斗底部與接料底板的距離為400 mm。采用隨機(jī)分布方式生成土壤顆粒,顆粒半徑設(shè)置為2~4 mm,仿真參數(shù)設(shè)置詳見表1.

    表1 土壤堆積角EDEM仿真參數(shù)

    1.3 標(biāo)定試驗設(shè)計

    根據(jù)不同含水條件的黏重黑土顆粒之間的粘結(jié)性質(zhì),擬定顆粒之間的JKR表面能參數(shù)、恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動摩擦系數(shù)作為待標(biāo)定的接觸模型參數(shù),采用響應(yīng)面優(yōu)化方法,以堆積角為試驗指標(biāo),在Design-expert軟件平臺上,運用Box-Behnken方法進(jìn)行試驗設(shè)計,參考EDEM 通用顆粒材料數(shù)據(jù)庫GEMM ( generic EDEM material model database,GEMM)中土壤顆粒的相關(guān)參數(shù),結(jié)合仿真預(yù)試驗的結(jié)果及文獻(xiàn)[22],堆積角仿真試驗因素的水平值設(shè)置如表2所示。

    表2 土壤堆積角仿真試驗因素及水平

    在EDEM2.7軟件平臺下,分別從錐面的方向和方向截取圖片,導(dǎo)入CAD軟件中標(biāo)注土壤堆積角[20],共測量4個堆積角,取平均值,仿真試驗結(jié)果如表3所示。

    表3 土壤堆積角仿真試驗設(shè)計及結(jié)果

    1.4 土壤堆積角仿真試驗結(jié)果與分析

    1.4.1 堆積角回歸模型

    按表2進(jìn)行堆積角仿真試驗,并對試驗結(jié)果進(jìn)行回歸方差分析,結(jié)果如表4所示,通過對試驗數(shù)據(jù)的多元回歸,得到東北黏重黑土的堆積角回歸模型為

    回歸方程(2)的決定系數(shù)2= 0.997 7,校正決定系數(shù)Adj-2=0.983 8。

    表4 土壤堆積角回歸模型的方差分析

    1.4.2 不同含水率的黏重黑土堆積角參數(shù)優(yōu)化與驗證

    通過Design-expert軟件的優(yōu)化模塊,以12.46%±1%含水率的黏重黑土形成的堆積角40.92°為目標(biāo)對堆積角回歸模型進(jìn)行尋優(yōu),得到多組優(yōu)化解,通過堆積角仿真試驗驗證這些優(yōu)化解,選取與物理試驗得到的土壤堆積角相近的1組優(yōu)化解,即土壤顆粒之間的JKR表面能7.46 J/m2,恢復(fù)系數(shù)0.56,靜摩擦系數(shù)0.78,滾動摩擦系數(shù)0.22,此優(yōu)化解下的堆積角仿真結(jié)果為 40.31°,與物理試驗結(jié)果的相對誤差為1.49%。圖3為堆積角仿真試驗與物理試驗堆形的對比。從圖3可以看出,對于優(yōu)化參數(shù)的仿真試驗,土壤堆形與物理試驗得到的結(jié)果具有較高的相似性。在土壤顆粒堆積角標(biāo)定的4個參數(shù)中,影響最顯著的是滾動摩擦系數(shù),其次是恢復(fù)系數(shù)與表面能參數(shù)。對于含水率12.46%±1%的黏重黑土,仿真參數(shù)在一定范圍內(nèi)取值均能獲得在物理試驗堆積角誤差范圍之內(nèi)的結(jié)果,仿真結(jié)果與試驗結(jié)果的相對誤差在2.5%以內(nèi),即滾動摩擦系數(shù)在0.2~0.24之間,JKR表面能范圍在6~8 J/m2之間,恢復(fù)系數(shù)在0.45~0.6之間,靜摩擦系數(shù)在0.65~0.9之間。同理,對于含水率17.15%±1%的黏重黑土,以44.15°的堆積角為優(yōu)化目標(biāo),最終得到土壤間JKR表面能參數(shù)為8.41 J/m2,土壤顆粒的恢復(fù)系數(shù)為0.5,靜摩擦系數(shù)系數(shù)為0.8,土壤間的滾動摩擦系數(shù)為0.23,利用該優(yōu)化參數(shù)獲得的堆積角仿真試驗結(jié)果與物理試驗結(jié)果在角度和堆形上均具有較高的相似性。同理,得到含水率17.15%±1%黏重黑土的仿真參數(shù)范圍,即滾動摩擦系數(shù)在0.22~0.24之間,JKR表面能范圍在7~9 J/m2之間,恢復(fù)系數(shù)在0.45~0.6之間,靜摩擦系數(shù)在0.65~0.9之間。

    圖3 不同含水率黏重土壤堆積角仿真與物理試驗的堆形對比 Fig.3 Stacked shape comparison between simulation and physical test of repose angle for clayey soil under different moisture content

    2 土壤與不同觸土部件材料的接觸模型參數(shù)標(biāo)定

    2.1 材料與方法

    為了給針對黏重黑土地區(qū)觸土部件的減阻性能優(yōu)化及仿生設(shè)計提供準(zhǔn)確可靠的仿真參數(shù),在1.1節(jié)研究結(jié)果的基礎(chǔ)上,針對典型觸土部件材料65Mn和2種減粘材料(超高分子量聚乙烯UHMW-PE和聚四氟乙烯PTFE),采用靜摩擦試驗測定4種含水率的黏重黑土對典型觸土部件材料的靜摩擦系數(shù),并基于靜摩擦試驗結(jié)果,采用斜面試驗標(biāo)定高含水率(含水率為17.15%±1%)的黏重黑土對2種典型觸土部件材料(65Mn和PTFE)的離散元仿真參數(shù)。

    2.1.1 土壤靜摩擦試驗

    土壤對不同材料的靜摩擦系數(shù)是土壤與觸土部件材料之間的特有屬性,不會因為土壤的形狀而改變,但是土壤的含水率對其具有較大影響[29-30]。

    為了減少仿真因素數(shù),提高仿真優(yōu)化水平與精度,土壤對不同材料的靜摩擦系數(shù)由如圖4所示的靜摩擦試驗臺測得。按前述方法配置含水率分別為0,12.46%,17.15%和23.5%的4種黏重黑土,制成1 cm×1 cm×1 cm的土塊,每次試驗前,保證土塊的含水率保持均衡,誤差控制在±1%。每種含水率的黏重黑土均重復(fù)測試20次,利用斜面儀測量每次土塊下滑瞬間的角度,通過式 (3)~(4)計算得到對應(yīng)的靜摩擦系數(shù),取平均值作為最終結(jié)果。

    注:F1為拉力,N;F2為壓力,N;G為重力,N;f為摩擦力,N;N為支撐力,N;θ為底座與斜面夾角,(°)。B為底座長度,mm;H為斜面末端高度,mm。

    2.1.2 斜面試驗

    表5 斜面試驗臺參數(shù)

    利用Desin-expert軟件中的Box-Behnken法進(jìn)行斜面試驗仿真參數(shù)標(biāo)定的試驗設(shè)計,參考EDEM通用顆粒材料數(shù)據(jù)庫GEMM ( generic EDEM material model database,GEMM)中鋼鐵及塑料的相關(guān)參數(shù),根據(jù)仿真預(yù)試驗及物理試驗的結(jié)果及文獻(xiàn)[25],擬定土壤顆粒與65Mn及PTFE板之間的斜面仿真試驗因素與水平,如表 6所示。

    注:S為土球自然滾動距離,mm。

    Fig 5 Slope test of soil and different contact materials

    表6 斜面仿真試驗因素及水平

    2.2 結(jié)果與分析

    2.2.1 靜摩擦試驗結(jié)果與分析

    圖6為靜摩擦系數(shù)試驗測試結(jié)果。通過靜摩擦試驗,測得含水率為0(干土)、12.46%±1%、17.15%±1%、23.5%±1%的黏重黑土與PTFE板之間的靜摩擦系數(shù)依次為0.447±0.075、0.332±0.034、0.522±0.11、0.633±0.16;與UHMW-PE板之間的靜摩擦系數(shù)依次為0.351±0.09、0.435±0.069、0.593±0.17、0.83±0.25;含水率為0(干土)、12.46%±1%、17.15%±1%的黏重黑土與65Mn板之間的靜摩擦系數(shù)依次為0.197±0.03、0.351±0.075、0.571±0.16。由圖6可知,對于觸土材料65Mn,隨著土壤含水率的增大,其靜摩擦系數(shù)亦增大且增大幅度明顯,而對于PTFE板與UHMW-PE板,隨著土壤含水率的增加,靜摩擦系數(shù)提高的趨勢較平緩,當(dāng)含水率超過23.5%時,土壤仍然具有下滑趨勢,摩擦系數(shù)并無顯著增大。

    圖7為含水率23.5%的黏重黑土粘附在豎直(90°)的65Mn鋼板上。如圖7所示,含水率為23.5%(小于液限值)時,試驗土壤已經(jīng)完全粘附在65Mn板上,試驗臺升角>90°后仍然不掉落,這說明在此含水率下的黏重黑土與65Mn板觸土面的粘附力大于其自身重力。而在低表面能的PTFE板和UHMW-PE板上含水率為23.5%的黏重黑土在一定升角范圍內(nèi)自由滑動,說明低表面能觸土部件具有減粘效果,這為仿生減粘觸土部件的設(shè)計提供了思路。

    圖6 不同含水率黏重黑土與不同觸土材料間的靜摩擦系數(shù)

    圖7 含水率23.5%的黏重黑土粘附在豎直65Mn板上

    2.2.2 斜面仿真試驗結(jié)果與分析

    1)土壤與65Mn板的斜面仿真試驗結(jié)果與分析

    含水率為17.15%±1%的土球在65Mn板上的斜面仿真試驗結(jié)果見表7,滾動距離回歸模型的方差分析見表8。

    對表7中的試驗數(shù)據(jù)進(jìn)行多元回歸,得到土壤顆粒在65Mn板上的滾動距離回歸模型,其回歸方程為

    表7 含水率17.15%的土壤在65Mn板上的斜面仿真試驗結(jié)果

    表8 土壤在65Mn板上滾動距離的回歸模型方差分析

    2)土壤與PTFE板的斜面仿真試驗結(jié)果分析

    土球在PTFE板上的斜面仿真試驗結(jié)果見表9,回歸模型的方差分析見表10。

    表9 含水率17.15%的土壤在PTFE板上的斜面仿真試驗結(jié)果

    對表9中的試驗數(shù)據(jù)進(jìn)行多元回歸,得到土壤顆粒在PTFE板上的滾動距離回歸模型,其回歸方程為

    表10 土壤在PTFE板上滾動距離的回歸模型方差分析

    回歸方程(6)的決定系數(shù)2= 0.999 8,校正決定系數(shù)Adj-2=0.998 9。

    對含水率為17.15%±1%的土球,測得其在PTFE板上的滾動距離為269.55±11.25 mm。利用Design-expert軟件中的優(yōu)化模塊,以滾動距離269.55 mm為目標(biāo)對回歸模型進(jìn)行尋優(yōu),選取與土球滾動距離相近的1組優(yōu)化解,逐一進(jìn)行仿真驗證,結(jié)合前文靜摩擦試驗結(jié)果,選取黏重黑土與PTFE板的靜摩擦系數(shù)為0.52,仿真結(jié)果顯示當(dāng)土壤與PTFE板之間的JKR表面能為4.08 J/m2、恢復(fù)系數(shù)為0.6、靜摩擦系數(shù)為0.52、滾動摩擦系數(shù)為0.045時,滾動距離仿真結(jié)果為269.35 mm,與物理試驗得到的平均滾動距離269.55 mm接近,相對誤差為0.07%。

    3 結(jié) 論

    1)針對東北地區(qū)高含水率的黏重黑土,基于離散元方法,采用Hertz-Mindlin with JKR Cohesion模型分別對含水率為12.46%±1%與含水率17.15%±1%的2種土壤進(jìn)行離散元仿真,通過堆積角仿真試驗標(biāo)定2種含水率下的黏重黑土顆粒之間的接觸模型仿真參數(shù),與物理試驗進(jìn)行對比,得到堆積角回歸方程,利用響應(yīng)面優(yōu)化方法得到2種含水率的黏重黑土接觸模型的參數(shù)分別為:含水率12.46%±1%下土壤顆粒之間滾動摩擦系數(shù)為0.22,JKR表面能為7.46 J/m2,恢復(fù)系數(shù)為0.56,靜摩擦系數(shù)0.78;含水率17.15%±1%的土壤顆粒JKR表面能為8.41 J/m2,恢復(fù)系數(shù)為0.5,靜摩擦系數(shù)系數(shù)為0.8,土壤之間的滾動摩擦系數(shù)0.23。由此可得,對于含水率在10%~20%的黏重黑土,在進(jìn)行相關(guān)離散元仿真試驗時,各對應(yīng)參數(shù)值可取上下限之間的值。

    3)通過斜面物理試驗與仿真試驗相結(jié)合的方法,得到含水率為17.15%±1%的黏重黑土對2類觸土部件材料的離散元仿真參數(shù),即黏重黑土與PTFE板之間的JKR表面能為4.08 J/m2,恢復(fù)系數(shù)為0.6,靜摩擦系數(shù)為0.52,滾動摩擦系數(shù)為0.045;與65Mn板之間的JKR表面能為5.5 J/m2,恢復(fù)系數(shù)為0.61,靜摩擦系數(shù)為0.57,滾動摩擦系數(shù)為0.056。

    4)通過2種含水率的黏重黑土堆積角的回歸模型方差分析以及土壤對2種觸土部件材料65Mn、PTFE的斜面試驗的回歸模型方差分析,得到東北地區(qū)黏重黑土顆粒之間(或顆粒與觸土部件材料之間)的滾動摩擦系數(shù)以及顆粒之間(或顆粒對觸土部件材料)的表面能參數(shù)對接觸模型影響均極顯著,這與物理試驗現(xiàn)象相吻合。

    本文研究結(jié)果對東北黏重黑土地區(qū)農(nóng)業(yè)生產(chǎn)具有一定指導(dǎo)意義,含水率越高,黏重黑土的粘附性越強,對觸土部件的粘附亦越嚴(yán)重,因此,對于黏重黑土地區(qū)土壤進(jìn)行耕作,可以選擇減粘材料的觸土部件或者設(shè)計仿生減粘觸土部件,有利于減少土壤粘附;對于塊莖類作物的機(jī)械收獲作業(yè),可以盡量讓土壤晾曬一些時日,以降低土壤含水率,減少土壤與觸土部件的粘附,且盡量避免雨后作業(yè)。

    [1] Mustafa Ucgul, John M Fielke, Chris Saunders. 3D DEM tillage simulation: Validation of a hysteretic spring (plastic) contact model for a sweep tool operating in a cohesionless soil[J]. Soil & Tillage Research, 2014(144): 220-227.

    [2] 張銳,李建橋,周長海,等. 推土板表面形態(tài)對土壤動態(tài)行為影響的離散元模擬[J]. 農(nóng)業(yè)工程學(xué)報,2007,23(9):13-19. Zhang Rui, Li Jianqiao, Zhou Changhai, et al. Simulation of dynamic behavior of soil ahead of the bulldo zing plates with different surface configurations by discrete element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(9): 13-19. (in Chinese with English abstract).

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

    [4] 王憲良,胡紅,王慶杰,等. 基于離散元的土壤模型參數(shù)標(biāo)定方法[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2017,48(12):78-85. Wang Xianliang, Hu Hong, Wang Qingjie, et al. Calibration method of soil contact characteristic parameters based on DEM theory[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(12): 78-85. (in Chinese with English abstract).

    [5] B Li, Y Chen, J Chen. Modeling of soil–claw interaction using the discrete element method (DEM)[J]. Soil & Tillage Research, 2016(158): 177–185.

    [6] 鄭侃,何進(jìn),李洪文,等. 基于離散元深松土壤模型的折線破土刃深松鏟研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2016,47(9):62-72. Zheng Kan, He Jin, Li Hongwen, et al. Research on polyline soil-breaking blade subsoiler based on subsoiling soil model using discrete element method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(9): 62-72. (in Chinese with English abstract).

    [7] 賀一鳴,吳明亮,向偉,等. 離散元法在農(nóng)業(yè)工程領(lǐng)域的應(yīng)用進(jìn)展[J]. 中國農(nóng)學(xué)通報,2017,33(20):133-137. He Yiming, Wu Mingliang, Xiang Wei, et al. Application progress of discrete element method in agricultural engineering[J]. Chinese Agricultural Science Bulletin, 2017, 33(20): 133-137. (in Chinese with English abstract).

    [8] 鹿芳媛,馬旭,譚穗妍,等. 水稻芽種離散元主要接觸參數(shù)仿真標(biāo)定與試驗. 農(nóng)業(yè)機(jī)械學(xué)報[J],2018,49(2):94-99. Lu Fangyuan, Ma Xu, Tan Suiyan, et al. Simulative calibration and experiment on main contact parameters of discrete elements for rice bud seeds[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(2): 94-99. (in Chinese with English abstract).

    [9] 王云霞,梁志杰,張東興,等. 基于離散元的玉米種子顆粒模型種間接觸參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(22):36-42. Wang Yunxia, Liang Zhijie, Zhang Dongxing, et al. Calibration method of contact characteristic parameters for corn seeds based on EDEM[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(22): 36-42. (in Chinese with English abstract)

    [10] 劉凡一,張艦,李博,等. 基于堆積試驗的小麥離散元參數(shù)分析及標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(12):247-253. Liu Fanyi, Zhang Jian, Li Bo, et al. Calibration of parameters of wheat required in discrete element method simulation based on repose angle of particle heap[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(12): 247-253. (in Chinese with English abstract)

    [11] 劉凡一,張艦,陳軍. 小麥籽粒振動篩分黏彈塑性接觸模型構(gòu)建及其參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(15):37-43. Liu Fanyi, Zhang Jian, Chen Jun. Construction of visco-elasto-plasticity contact model of vibratory screening and its parameters calibration for wheat[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(15): 37-43. (in Chinese with English abstract)

    [12] 彭飛,王紅英,方芳,等. 基于注入截面法的顆粒飼料離散元模型參數(shù)標(biāo)定[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2018,49(4):140-147. Peng Fei, Wang Hongying, Fang Fang, et al. Calibration of discrete element model parameters for pellet feed based on injected cection method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(4): 140-147. (in Chinese with English abstract)

    [13] 張濤,劉飛,趙滿全,等. 大豆種子與排種器接觸物理參數(shù)的測定與離散元仿真標(biāo)定[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報,2017,22(9):86-92. Zhang Tao, Liu Fei, Zhao Manquan, et al. Measurement of physical parameters of contact between soybean seed and metering device and discrete element simulation calibration[J]. Journal of China Agricultural University, 2017, 22(9): 86-92. (in Chinese with English abstract).

    [14] Mustafa Ucgul, John M Fielke, Chris Saunders. Defining the effect of sweep tillage tool cutting edge geometry on tillage forces using 3D discrete element modeling[J]. Information Processing in Agriculture, 2015(2): 130-141.

    [15] Mustafa Ucgul, John M Fielke, Chris Saunders. Three-dimensional discrete element modelling of tillage: Determination of a suitable contact model and parameters for a cohesionless soil[J]. Biosystems Engineering, 2014(121): 105-107.

    [16] 于建群,付宏,李紅,等. 離散元法及其在農(nóng)業(yè)機(jī)械工作部件研究與設(shè)計中的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報,2005,21(5):1-6. Yu Jianqun, Fu Hong, Li Hong, et al. Application of discrete element method to research and design of working parts of agricultural machines[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2005, 21(5): 1-6. (in Chinese without English abstract).

    [17] Tong Jin, Mohammad Almagzoub Mohammad, Zhang Jinbo, et al. DEM numerical simulation of abrasive wear characteristics of a bioinspired ridged surface[J]. Journal of Bionic Engineering, 2010 (7): 175-181.

    [18] 丁啟朔,任駿,Belal Eisa Adam,等. 濕粘水稻土深松過程離散元分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2017,48(3):38-48. Ding Qishuo, Ren Jun, Belal Eisa Adam, et al. DEM analysis of subsoiling process in wet clayey paddy soil[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(3): 38-48. (in Chinese with English abstract)

    [19] 張銳,韓佃雷,吉巧麗,等. 離散元模擬中沙土參數(shù)標(biāo)定方法研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2017,48(3):50-56. Zhang Rui, Han Dianlei, Ji Qiaoli, et al. Calibration methods of sandy soil parameters in simulation of discrete element method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(3): 50-56. (in Chinese with English abstract).

    [20] 石林榕,趙武云,孫偉. 基于離散元的西北旱區(qū)農(nóng)田土壤顆粒接觸模型和參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(21):181-187. Shi Linrong, Zhao Wuyun, Sun Wei. Parameter calibration of soil particles contact model of farmland soil in northwest arid region based on discrete element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(21): 181-187. (in Chinese with English abstract)

    [21] 石林榕,吳建民,趙武云,等. 基于離散單元法農(nóng)田土壤單軸壓縮模型的建立及參數(shù)校核[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報,2015,20(4):174-182. Shi Linrong, Wu Jianmin, Zhao Wuyun, et al. Establishement and parameter verification of farmland soil model in uniaxial compression based on discrete element method[J]. Journal of China Agricultural University, 2015, 20(4): 174-182. (in Chinese with English abstract)

    [22] 武濤,黃偉鳳,陳學(xué)深,等. 考慮顆粒間黏結(jié)力的黏性土壤離散元模型參數(shù)標(biāo)定[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報,2017,38(3):93-98. Wu Tao, Huang Weifeng, Chen Xueshen, et al. Calibration of discrete element model parameters for cohesive soil considering the cohesion between particles[J]. Journal of South China Agricultural University, 2017, 38(3): 93-98. (in Chinese with English abstract)

    [23] EDEM 2.5 theory reference guide[R/OL]. 2014-12-05 [2015-06-07].http://www.docin.com/p-980174717. html

    [24] 熊平原,楊洲,孫志全,等. 基于離散元法的旋耕刀三向工作阻力仿真分析與試驗[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(18):113-121. Xiong Pingyuan, Yang Zhou, Sun Zhiquan, et al. Simulation analysis and experiment for three-axis working resistances of rotary blade based on discrete element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 113-121. (in Chinese with English abstract)

    [25] 郝新敏,楊元,黃斌香. 聚四氟乙烯微孔膜及纖維[M]. 北京:化學(xué)工業(yè)出版社,2011.

    [26] 滕兵. 端面扭動摩擦界面行為研究[D]. 北京:中國礦業(yè)大學(xué),2014. Teng Bing. Study on Face-on-Face Torsional Tribological Interfacial Behaviors[D]. Beijing: China University of Mining and Technology, 2014. (in Chinese with English abstract)

    [27] Mustafa Ucgul, Chris Saunders, John M Fielke. Discrete element modelling of tillage forces and soil movement of a one-third scale mouldboard plough[J]. Biosystems Engineering, 2017(155): 44-54.

    [28] J Y Sun, Y M Wang, Y H Ma, et al. DEM simulation of bionic subsoilers (tillage depth >40 cm) with drag reduction and lower soil disturbance characteristics[J]. Advances in Engineering Software, 2018(119): 30-37.

    [29] 任露泉. 土壤粘附力學(xué)[M]. 北京:機(jī)械工業(yè)出版社. 2011.

    [30] 錢定華,張際先. 土壤對金屬材料粘附和摩擦研究狀況概述[J]. 農(nóng)業(yè)機(jī)械學(xué)報,1984(1):69-78. Qian Dinghua, Zhang Jixian. A Summary of study of adhesion and friction between soil and metals[J]. Transactions of the Chinese Society for Agricultural Machinery, 1984(1): 69-78. (in Chinese with English abstract).

    Calibration of parameters of interaction between clayey black soil with different moisture content and soil-engaging component in northeast China

    Li Junwei1,2,3,4, Tong Jin1,3, Hu Bin2,4, Wang Hubiao1,3, Mao Chunyu1,3,5, Ma Yunhai1,3※

    (1.,,,130022,; 2.,,832003,3.,,130022,4.,,832003,5.,130022,;)

    In order to obtain the interaction parameters which can be used for discrete element simulation between the soil-engaging components and clayey black soil in northeast China, the Hertz-Mindlin with JKR Cohesion contact model in EDEM2.7 was used to simulate the clayey black soil. Aiming at the actual working environment with moisture content of 10%-20%, 2 kinds of clayey black soil with moisture content of 12.46%±1% and 17.15%±1% were prepared respectively. The rolling friction coefficient, recovery coefficient, JKR surface energy parameter and the static friction coefficient between the soil particles were used as the model calibrated parameters. The soil repose angle simulated under the simulation parameters was set as the response value. The regression model of the soil repose angle was established based on the Box-Behnken response surface method, and the regression model was optimized by using the soil repose angles under 2 moisture contents obtained by physical experiments. The optimal solution of the contact model parameters of the clayey black soil particles with 2 moisture contents was obtained. The parameter range of the discrete element contact model of clayey black soil with moisture content between 10%-20% were given. The soil repose angle and stacked shape obtained by the simulation experiment with the optimized solution had a high similarity with that of the physical test. At the same time, the static friction coefficients of 3 kinds of soil contact materials such as 65Mn plate, UHMW-PE plate and PTFE plate were tested under the conditions of 4 kinds of moisture content, which were 0, 12.46%, 17.15% and 23.5% respectively. Then for the 2 typical soil contact materials such as 65Mn and PTFE, the slope physical test bench were constructed for the clayey black soil with moisture content of 17.15%±1%. The slope simulation models of 65Mn and PTFE materials were established under the EDEM2.7 software platform. Based on the critical simulation parameters such as JKR surface energy parameter, static friction coefficient, rolling friction coefficient and recovery coefficient between soil particle and soil contact material, the simulated rolling distance of the soil ball on different soil contact materials was set as the response value. The excellent fit regression model of rolling distance for the 2 kinds of typical soil contact materials were established based on the response surface optimization method of Box-Behnkende. The 2 regression models were optimized by the sliding distance obtained by the physical experiment. Then the optimized solution of the contact model parameters of different materials was obtained. The simulation results showed that when the surface energy of JKR between soil and 65Mn was 5.5 J/m2, the recovery coefficient was 0.61, the static friction coefficient was 0.57 and the rolling friction coefficient was 0.056, the rolling distance simulation result was 153.56 mm, which was close to the average rolling distance obtained from physical experiments of 155.93 mm, and the relative error was 1.52%; When the surface energy of JKR between soil and PTFE was 4.08 J/m2, the recovery coefficient was 0.6, the static friction coefficient was 0.52 and the rolling friction coefficient was 0.045, the simulation result of rolling distance was 269.35 mm, which was close to the average rolling distance of 269.55 mm obtained from physical experiment, and the relative error was 0.07%. And the optimized parameters can be used to simulate the discrete parameters between the clayey black soil and the soil-engaging components. The study provides credible basic data for the design and simulation of agricultural machinery under clayey black soil conditions.

    soils; water content; discrete element method; clayey black soil; calibration of simulation parameter

    2018-10-11

    2019-03-03

    國家重點研發(fā)計劃項目(YS2016YFNC050065)課題(2016YFD0701601);國家自然科學(xué)基金項目(51805338,51875242,51865051);吉林省科技廳項目(20170101173JC,20190302129GX)

    李俊偉,講師,博士生,主要從事農(nóng)業(yè)機(jī)械減阻耐磨仿生技術(shù)研究。Email:ljwjdxy@shzu.edu.cn

    馬云海,教授,博士生導(dǎo)師,主要從事農(nóng)業(yè)工程仿生技術(shù)研究。Email:myh@jlu.edu.cn

    10.11975/j.issn.1002-6819.2019.06.016

    S347.7;S22

    A

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

    李俊偉,佟 金,胡 斌,王虎彪,毛春昱,馬云海. 不同含水率黏重黑土與觸土部件互作的離散元仿真參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(6):130-140. doi:10.11975/j.issn.1002-6819.2019.06.016 http://www.tcsae.org

    Li Junwei, Tong Jin, Hu Bin, Wang Hubiao, Mao Chunyu, Ma Yunhai. Calibration of parameters of interaction between clayey black soil with different moisture content and soil-engaging component in northeast China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(6): 130-140. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.06.016 http://www.tcsae.org

    猜你喜歡
    黑土摩擦系數(shù)標(biāo)定
    自然資源部:加強黑土耕地保護(hù)
    基于Sentinel-2遙感影像的黑土區(qū)土壤有效磷反演
    隧道內(nèi)水泥混凝土路面微銑刨后摩擦系數(shù)衰減規(guī)律研究
    中外公路(2022年1期)2022-05-14 08:13:26
    摩擦系數(shù)對直齒輪副振動特性的影響
    使用朗仁H6 Pro標(biāo)定北汽紳寶轉(zhuǎn)向角傳感器
    基于勻速率26位置法的iIMU-FSAS光纖陀螺儀標(biāo)定
    船載高精度星敏感器安裝角的標(biāo)定
    CSP生產(chǎn)線摩擦系數(shù)與軋制力模型的研究
    上海金屬(2014年3期)2014-12-19 13:09:12
    基于Harris-張正友平面標(biāo)定法的攝像機(jī)標(biāo)定算法
    典型黑土區(qū)不同尺度觀測場地融雪徑流
    日本黄色视频三级网站网址 | 露出奶头的视频| 男女床上黄色一级片免费看| 老鸭窝网址在线观看| 精品国产一区二区久久| 精品第一国产精品| 国内久久婷婷六月综合欲色啪| a级毛片在线看网站| 国产91精品成人一区二区三区| 久久 成人 亚洲| 中文欧美无线码| 老汉色av国产亚洲站长工具| 国产亚洲av高清不卡| 老司机午夜福利在线观看视频| 欧美丝袜亚洲另类 | 国产精品自产拍在线观看55亚洲 | 欧美国产精品va在线观看不卡| 国产精品免费视频内射| 精品欧美一区二区三区在线| 热99re8久久精品国产| 大码成人一级视频| 久久精品亚洲熟妇少妇任你| 一区二区三区激情视频| www.999成人在线观看| 黄片大片在线免费观看| 国产人伦9x9x在线观看| 亚洲一区中文字幕在线| 极品少妇高潮喷水抽搐| 亚洲人成77777在线视频| 两性夫妻黄色片| 久久热在线av| 999精品在线视频| av视频免费观看在线观看| 亚洲av第一区精品v没综合| 黄色片一级片一级黄色片| 免费观看精品视频网站| 激情视频va一区二区三区| 在线永久观看黄色视频| 久久精品91无色码中文字幕| 亚洲欧洲精品一区二区精品久久久| 午夜福利影视在线免费观看| 十八禁高潮呻吟视频| 久久精品人人爽人人爽视色| 极品人妻少妇av视频| 日日爽夜夜爽网站| 日韩免费av在线播放| 国产一区二区激情短视频| 午夜精品在线福利| 黄色怎么调成土黄色| 国产亚洲精品久久久久久毛片 | 国产高清videossex| 一边摸一边抽搐一进一出视频| 国产一区有黄有色的免费视频| 成在线人永久免费视频| 久久久久久亚洲精品国产蜜桃av| 国产在视频线精品| www.自偷自拍.com| 俄罗斯特黄特色一大片| 老司机午夜福利在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 十八禁高潮呻吟视频| 国产男女超爽视频在线观看| 日本撒尿小便嘘嘘汇集6| 看片在线看免费视频| 亚洲全国av大片| 在线天堂中文资源库| 欧美成狂野欧美在线观看| 成人18禁高潮啪啪吃奶动态图| 九色亚洲精品在线播放| 80岁老熟妇乱子伦牲交| 99riav亚洲国产免费| 久久久久精品国产欧美久久久| 欧美日韩精品网址| 午夜免费观看网址| 欧美日韩视频精品一区| 黄色毛片三级朝国网站| 日韩欧美国产一区二区入口| 久久久久久久午夜电影 | 热re99久久国产66热| 亚洲专区中文字幕在线| 天天添夜夜摸| 久久久久视频综合| 久久久久国产精品人妻aⅴ院 | 91老司机精品| 欧美日韩亚洲高清精品| 大香蕉久久网| 日韩一卡2卡3卡4卡2021年| 亚洲成人手机| 成人国产一区最新在线观看| 午夜亚洲福利在线播放| 可以免费在线观看a视频的电影网站| 无限看片的www在线观看| 日本五十路高清| 一个人免费在线观看的高清视频| 久久这里只有精品19| 99国产精品99久久久久| av天堂在线播放| 一二三四社区在线视频社区8| 久久久久久久国产电影| 亚洲黑人精品在线| 午夜福利在线免费观看网站| 久久久精品免费免费高清| 久久狼人影院| 国产深夜福利视频在线观看| 无遮挡黄片免费观看| 好男人电影高清在线观看| 久久国产精品男人的天堂亚洲| 黄片小视频在线播放| 露出奶头的视频| 国产成人影院久久av| 性少妇av在线| 久久影院123| 亚洲欧美一区二区三区久久| 国产精品免费视频内射| 日韩欧美国产一区二区入口| 另类亚洲欧美激情| 成在线人永久免费视频| 久久天躁狠狠躁夜夜2o2o| aaaaa片日本免费| 国产精品久久视频播放| 国产男女内射视频| 久久久精品区二区三区| 大陆偷拍与自拍| 免费在线观看亚洲国产| 高清欧美精品videossex| 亚洲少妇的诱惑av| 国产成人精品无人区| 精品午夜福利视频在线观看一区| 老熟妇乱子伦视频在线观看| 一区福利在线观看| 色精品久久人妻99蜜桃| 久久ye,这里只有精品| 国产亚洲av高清不卡| 一a级毛片在线观看| 亚洲,欧美精品.| 国产成人系列免费观看| 欧美老熟妇乱子伦牲交| 国产淫语在线视频| 99久久人妻综合| 国产国语露脸激情在线看| 一区福利在线观看| 免费女性裸体啪啪无遮挡网站| 天堂中文最新版在线下载| 亚洲一区二区三区欧美精品| 中亚洲国语对白在线视频| 国产成人av教育| 成人18禁高潮啪啪吃奶动态图| 热re99久久国产66热| 亚洲熟女毛片儿| 十八禁高潮呻吟视频| 18禁国产床啪视频网站| videosex国产| 黄频高清免费视频| 亚洲五月天丁香| 国产99白浆流出| 欧美+亚洲+日韩+国产| 丰满的人妻完整版| 欧美av亚洲av综合av国产av| 18禁国产床啪视频网站| 1024视频免费在线观看| 黄片大片在线免费观看| 夜夜夜夜夜久久久久| av天堂久久9| 国产男女超爽视频在线观看| 91精品三级在线观看| 精品国内亚洲2022精品成人 | 亚洲国产欧美一区二区综合| 91字幕亚洲| av电影中文网址| 亚洲中文日韩欧美视频| 国产精品永久免费网站| 亚洲一区二区三区不卡视频| 亚洲全国av大片| 丰满的人妻完整版| 又大又爽又粗| 亚洲精品国产色婷婷电影| 国产精品国产高清国产av | av网站免费在线观看视频| 国产单亲对白刺激| 99国产精品99久久久久| 老熟女久久久| 久久人妻熟女aⅴ| 中出人妻视频一区二区| 国产在视频线精品| 啦啦啦免费观看视频1| 高清黄色对白视频在线免费看| 中文字幕人妻熟女乱码| 午夜成年电影在线免费观看| 亚洲午夜精品一区,二区,三区| 国产不卡一卡二| 天天添夜夜摸| 国产区一区二久久| 日韩免费高清中文字幕av| 一级a爱片免费观看的视频| 欧美日本中文国产一区发布| 高清黄色对白视频在线免费看| 母亲3免费完整高清在线观看| 国产亚洲欧美98| 日日夜夜操网爽| 久久久久久免费高清国产稀缺| 午夜成年电影在线免费观看| 搡老熟女国产l中国老女人| 国产高清视频在线播放一区| 免费在线观看完整版高清| 精品福利永久在线观看| 成人国语在线视频| 丁香六月欧美| 国产99白浆流出| 99re6热这里在线精品视频| 久久久久久久精品吃奶| 欧美大码av| 动漫黄色视频在线观看| 丁香欧美五月| 欧美黄色淫秽网站| 欧美人与性动交α欧美精品济南到| 可以免费在线观看a视频的电影网站| 我的亚洲天堂| 成人影院久久| 亚洲色图综合在线观看| 欧美日韩国产mv在线观看视频| 超碰97精品在线观看| 亚洲,欧美精品.| 丝袜人妻中文字幕| 免费av中文字幕在线| 欧美不卡视频在线免费观看 | 免费观看精品视频网站| 欧美亚洲日本最大视频资源| 久久中文看片网| 在线看a的网站| 中出人妻视频一区二区| 一边摸一边抽搐一进一出视频| 丰满饥渴人妻一区二区三| 亚洲精品美女久久久久99蜜臀| 精品无人区乱码1区二区| 成年人午夜在线观看视频| 日本黄色日本黄色录像| 免费av中文字幕在线| 国产深夜福利视频在线观看| 亚洲五月婷婷丁香| 九色亚洲精品在线播放| 久久人妻熟女aⅴ| 黄色怎么调成土黄色| 巨乳人妻的诱惑在线观看| 波多野结衣一区麻豆| 又黄又爽又免费观看的视频| 亚洲av电影在线进入| 丰满人妻熟妇乱又伦精品不卡| 久久人人97超碰香蕉20202| 久久久久久久久久久久大奶| 亚洲国产看品久久| 大型av网站在线播放| 国产在线精品亚洲第一网站| 国产精品香港三级国产av潘金莲| 亚洲精品中文字幕一二三四区| 久久香蕉激情| 国产高清国产精品国产三级| 一级毛片高清免费大全| 飞空精品影院首页| 在线永久观看黄色视频| 久久久久国内视频| 亚洲精品乱久久久久久| 久久精品国产综合久久久| 久久精品国产a三级三级三级| 亚洲av日韩精品久久久久久密| 中文字幕人妻熟女乱码| 亚洲精品久久午夜乱码| 国产亚洲一区二区精品| 国产午夜精品久久久久久| 无人区码免费观看不卡| 亚洲av成人一区二区三| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区| 高清av免费在线| 少妇 在线观看| 午夜福利视频在线观看免费| 国产av精品麻豆| 人妻久久中文字幕网| 婷婷成人精品国产| 在线观看午夜福利视频| 亚洲一区中文字幕在线| 少妇裸体淫交视频免费看高清 | 宅男免费午夜| 亚洲一区二区三区欧美精品| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三| 午夜福利欧美成人| 天天躁狠狠躁夜夜躁狠狠躁| 日本a在线网址| 国内久久婷婷六月综合欲色啪| 亚洲精品国产区一区二| 国产亚洲欧美在线一区二区| 麻豆av在线久日| 亚洲五月天丁香| 精品人妻在线不人妻| 日韩视频一区二区在线观看| 亚洲人成电影免费在线| 啦啦啦视频在线资源免费观看| 好男人电影高清在线观看| 中文字幕人妻丝袜制服| 国产亚洲欧美98| 中文字幕av电影在线播放| 久久精品亚洲熟妇少妇任你| 久99久视频精品免费| 91大片在线观看| 久久国产精品影院| 午夜免费观看网址| 后天国语完整版免费观看| 在线国产一区二区在线| videos熟女内射| 一夜夜www| 美女扒开内裤让男人捅视频| 欧美日韩av久久| 无遮挡黄片免费观看| 国产精品自产拍在线观看55亚洲 | 免费日韩欧美在线观看| 1024香蕉在线观看| 亚洲人成电影免费在线| 一本综合久久免费| 精品人妻在线不人妻| 亚洲五月婷婷丁香| 看片在线看免费视频| 涩涩av久久男人的天堂| 丝袜美腿诱惑在线| 欧美中文综合在线视频| 成人特级黄色片久久久久久久| 亚洲熟妇熟女久久| 欧美黄色淫秽网站| 又黄又爽又免费观看的视频| 久久久久久久午夜电影 | 建设人人有责人人尽责人人享有的| 亚洲av第一区精品v没综合| 日日摸夜夜添夜夜添小说| 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 亚洲片人在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产免费男女视频| 精品国产一区二区三区久久久樱花| 黄色a级毛片大全视频| 人妻 亚洲 视频| 成人三级做爰电影| 国产xxxxx性猛交| 极品人妻少妇av视频| 老汉色∧v一级毛片| 久久久久国内视频| 久久国产精品影院| 亚洲aⅴ乱码一区二区在线播放 | 99久久综合精品五月天人人| 最近最新中文字幕大全免费视频| 黄色片一级片一级黄色片| 久久国产亚洲av麻豆专区| 精品亚洲成a人片在线观看| 亚洲性夜色夜夜综合| 免费久久久久久久精品成人欧美视频| 欧美日韩瑟瑟在线播放| 在线永久观看黄色视频| 国产av又大| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 成在线人永久免费视频| 十分钟在线观看高清视频www| 欧美人与性动交α欧美软件| 一级a爱片免费观看的视频| 亚洲av第一区精品v没综合| 成熟少妇高潮喷水视频| 欧美精品亚洲一区二区| 精品一品国产午夜福利视频| 欧美激情极品国产一区二区三区| 国产欧美日韩综合在线一区二区| 一级毛片精品| 乱人伦中国视频| 51午夜福利影视在线观看| 国产亚洲欧美精品永久| 777久久人妻少妇嫩草av网站| 亚洲综合色网址| 国产精品乱码一区二三区的特点 | 精品国产一区二区三区久久久樱花| 日本vs欧美在线观看视频| 又黄又爽又免费观看的视频| 国产主播在线观看一区二区| 一区在线观看完整版| 久久ye,这里只有精品| 动漫黄色视频在线观看| 亚洲av美国av| 一区在线观看完整版| 欧美人与性动交α欧美软件| 丰满的人妻完整版| 国产97色在线日韩免费| 久久婷婷成人综合色麻豆| 男女之事视频高清在线观看| 午夜成年电影在线免费观看| 日韩免费av在线播放| 最近最新免费中文字幕在线| 色播在线永久视频| 一级a爱片免费观看的视频| 国产精品香港三级国产av潘金莲| 色94色欧美一区二区| 国产亚洲欧美98| 精品国产超薄肉色丝袜足j| 国产精品一区二区在线观看99| 欧美国产精品一级二级三级| 久久人人97超碰香蕉20202| 久久亚洲真实| 中文字幕av电影在线播放| 久久精品熟女亚洲av麻豆精品| 又大又爽又粗| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩一区二区三| 两人在一起打扑克的视频| 精品国产超薄肉色丝袜足j| 中文字幕另类日韩欧美亚洲嫩草| 久久人人爽av亚洲精品天堂| 午夜免费成人在线视频| 在线观看午夜福利视频| 91九色精品人成在线观看| 国产欧美日韩一区二区三| aaaaa片日本免费| 亚洲熟女毛片儿| 免费女性裸体啪啪无遮挡网站| 久久久久久久久免费视频了| 久久久国产欧美日韩av| 久久九九热精品免费| 色婷婷久久久亚洲欧美| 午夜福利免费观看在线| 人人澡人人妻人| 国产在线观看jvid| 老司机午夜十八禁免费视频| 丰满饥渴人妻一区二区三| 亚洲九九香蕉| 夜夜夜夜夜久久久久| 高清欧美精品videossex| 美女午夜性视频免费| 男人操女人黄网站| 亚洲一区中文字幕在线| 精品电影一区二区在线| 黄频高清免费视频| 黄色毛片三级朝国网站| 久久精品国产亚洲av香蕉五月 | 欧美日韩黄片免| 久久亚洲真实| 久久精品国产综合久久久| 亚洲久久久国产精品| 亚洲精品国产一区二区精华液| 99国产极品粉嫩在线观看| 日韩免费高清中文字幕av| 亚洲精品在线美女| 交换朋友夫妻互换小说| 在线观看舔阴道视频| 午夜福利,免费看| 国产单亲对白刺激| 久久久久久人人人人人| 亚洲五月天丁香| 香蕉久久夜色| 亚洲精品粉嫩美女一区| 久久草成人影院| 热re99久久精品国产66热6| 一夜夜www| 成人精品一区二区免费| a级片在线免费高清观看视频| 99久久综合精品五月天人人| 午夜免费成人在线视频| 亚洲色图 男人天堂 中文字幕| 国产乱人伦免费视频| 亚洲九九香蕉| 真人做人爱边吃奶动态| 18禁国产床啪视频网站| 欧美成人免费av一区二区三区 | 亚洲精品粉嫩美女一区| 欧美精品av麻豆av| 美女午夜性视频免费| 欧美日韩一级在线毛片| 欧美不卡视频在线免费观看 | 久久国产精品影院| 久久人妻福利社区极品人妻图片| 人妻丰满熟妇av一区二区三区 | 久久久久久久久久久久大奶| 成人精品一区二区免费| 亚洲,欧美精品.| 国产精品亚洲av一区麻豆| 亚洲国产欧美日韩在线播放| 久久久久久久久久久久大奶| 女人爽到高潮嗷嗷叫在线视频| 国产午夜精品久久久久久| 夫妻午夜视频| 国产成人系列免费观看| 黄色片一级片一级黄色片| 激情视频va一区二区三区| 久久天躁狠狠躁夜夜2o2o| 午夜激情av网站| 亚洲欧美一区二区三区黑人| 国产精品久久久久成人av| 国内毛片毛片毛片毛片毛片| 女人爽到高潮嗷嗷叫在线视频| 99国产综合亚洲精品| 国产亚洲精品一区二区www | 咕卡用的链子| 成在线人永久免费视频| 欧美日韩国产mv在线观看视频| ponron亚洲| 大香蕉久久成人网| 黄网站色视频无遮挡免费观看| 亚洲一码二码三码区别大吗| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 日韩欧美在线二视频 | 久久久久久亚洲精品国产蜜桃av| 久久久国产欧美日韩av| 欧美黄色淫秽网站| 无遮挡黄片免费观看| 大码成人一级视频| 搡老岳熟女国产| 婷婷丁香在线五月| 亚洲国产毛片av蜜桃av| 久久精品国产99精品国产亚洲性色 | 午夜福利乱码中文字幕| 婷婷丁香在线五月| 国产成人系列免费观看| 婷婷丁香在线五月| 亚洲少妇的诱惑av| 久久久久久亚洲精品国产蜜桃av| 亚洲伊人色综图| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 美女高潮喷水抽搐中文字幕| 久久人妻av系列| 99精品欧美一区二区三区四区| 精品国产乱子伦一区二区三区| 日本黄色日本黄色录像| 看片在线看免费视频| 国产男女超爽视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 黑丝袜美女国产一区| 又黄又粗又硬又大视频| 狠狠狠狠99中文字幕| 又黄又爽又免费观看的视频| 久久亚洲真实| 黄网站色视频无遮挡免费观看| 欧美色视频一区免费| 91字幕亚洲| 欧美黄色淫秽网站| 另类亚洲欧美激情| 欧美日韩国产mv在线观看视频| 99riav亚洲国产免费| 精品亚洲成a人片在线观看| 丰满迷人的少妇在线观看| 国产有黄有色有爽视频| 妹子高潮喷水视频| 首页视频小说图片口味搜索| www.999成人在线观看| 欧美精品啪啪一区二区三区| 欧美日韩一级在线毛片| 久热这里只有精品99| 亚洲欧洲精品一区二区精品久久久| 亚洲三区欧美一区| 99香蕉大伊视频| 精品一区二区三卡| 久久人人爽av亚洲精品天堂| 丰满饥渴人妻一区二区三| 日韩欧美免费精品| 真人做人爱边吃奶动态| 欧美在线一区亚洲| 日韩欧美国产一区二区入口| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三卡| 精品欧美一区二区三区在线| 91成人精品电影| 高清毛片免费观看视频网站 | 91精品三级在线观看| 国产亚洲欧美精品永久| 每晚都被弄得嗷嗷叫到高潮| 99精品久久久久人妻精品| 亚洲色图av天堂| 免费观看人在逋| 国产精品乱码一区二三区的特点 | 久久久久久久精品吃奶| 好看av亚洲va欧美ⅴa在| 一夜夜www| 黄色毛片三级朝国网站| 国产av精品麻豆| 50天的宝宝边吃奶边哭怎么回事| 啦啦啦 在线观看视频| 精品第一国产精品| 久久香蕉国产精品| 两个人看的免费小视频| 久久久久久久午夜电影 | 欧美精品啪啪一区二区三区| 久久久国产成人免费| 18禁观看日本| 91大片在线观看| 国产精品一区二区在线观看99| 欧美精品啪啪一区二区三区| 99久久国产精品久久久| 老熟妇仑乱视频hdxx| 久久ye,这里只有精品| 精品免费久久久久久久清纯 | a在线观看视频网站| 亚洲第一av免费看| 色综合欧美亚洲国产小说| 色94色欧美一区二区| 精品国产超薄肉色丝袜足j| 国产精品美女特级片免费视频播放器 | 久久精品国产亚洲av香蕉五月 | 成人三级做爰电影| 乱人伦中国视频| 99久久综合精品五月天人人| 老熟女久久久| 国产精品亚洲av一区麻豆| 欧美+亚洲+日韩+国产| 少妇的丰满在线观看| 国产在线一区二区三区精| 国产黄色免费在线视频| 身体一侧抽搐| 亚洲人成电影观看| 丝袜在线中文字幕|