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

    基于DEM的3種勻盤結(jié)構(gòu)的玉米種子丸?;瘮?shù)值模擬

    2024-11-27 00:00:00李紅軍孫文杰張弛張成俊杜瑋陳偉
    中國粉體技術(shù) 2024年6期
    關(guān)鍵詞:數(shù)值模擬

    摘要:【目的】解決種子丸粒包衣質(zhì)量不高、種粉間混合均勻度差的問題?!痉椒ā糠治?種勻盤結(jié)構(gòu)的丸粒化機(jī)制,模擬玉米種子丸粒化過程;基于離散元軟件EDEM建立玉米種子丸?;姆抡婺P?,將離散系數(shù)作為評價指標(biāo)研究3種勻盤結(jié)構(gòu)的轉(zhuǎn)速、振動頻率、振動幅度對種粉間混合均勻度的影響;通過Design-Expert軟件建立玉米種子離散系數(shù)的二階回歸方程進(jìn)行方差分析。【結(jié)果】3種勻盤結(jié)構(gòu)的最優(yōu)工藝參數(shù)組合分別為:勻盤Ⅰ的轉(zhuǎn)速為50 r/min,振動頻率為10 Hz,振動幅度為6 mm;勻盤Ⅱ的轉(zhuǎn)速為60 r/min、振動頻率為15 Hz,振動幅度為6 mm;勻盤Ⅲ的轉(zhuǎn)速為70 r/min,振動頻率為20 Hz,振動幅度為8 mm?!窘Y(jié)論】當(dāng)勻盤Ⅰ的轉(zhuǎn)速為56. 3 r/min、振動頻率為11. 5 Hz,振動幅度為5. 8 mm時,玉米種子丸?;Ч顑?yōu),為玉米種子包衣技術(shù)提供理論依據(jù)自我評價。

    關(guān)鍵詞:丸?;浑x散元模型;包衣;數(shù)值模擬

    中圖分類號:TB44文獻(xiàn)標(biāo)志碼:A

    引用格式:

    李紅軍,孫文杰,張弛,等. 基于DEM的3種勻盤結(jié)構(gòu)的玉米種子丸?;瘮?shù)值模擬[J]. 中國粉體技術(shù),2024,30(6):149-161.

    LI Hongjun,SUN Wenjie,ZHANG Chi,et al. on pelletization of corn seeds with three uniform plate structures based on DEM[J]. China Powder Science and Technology,2024,30(6):149?161.

    20世紀(jì)初期,美國學(xué)者首次提出種子丸?;夹g(shù),經(jīng)過長期發(fā)展研究,到了90年代,西方國家的種子丸?;夹g(shù)已經(jīng)發(fā)展至基本成熟階段,丸粒化的合格率已高達(dá)90%以上[1-2]。我國對造粒技術(shù)的研究起步較晚,直至90年代才研制出各種類型的包衣機(jī),例如圓盤式、傾斜旋轉(zhuǎn)式、滾筒式,這些成功應(yīng)用于芝麻種子、煙草種子、蘿卜種子等,但是在發(fā)展過程中還存在許多技術(shù)問題,比如造粒合格率低、造粒后多種子率高、造粒包衣質(zhì)量差、造粒工藝比較落后等[3-4]。

    許多研究者針對不同種子丸?;c不同的包衣機(jī)進(jìn)行了研究與試驗。Barekatain等[5]提出對于小粒徑種子,丸?;蟛粌H可以滿足機(jī)械化精準(zhǔn)播種,還可以提高生產(chǎn)效率,節(jié)約生產(chǎn)成本。Sikhao等[6]通過研究花菜、萵筍、香菜等蔬菜種子的丸?;瘷C(jī)制,制訂了提高種子丸?;募庸すに?。劉俏等[7]在EDEM數(shù)值模擬的基礎(chǔ)上設(shè)計了種藥混合裝置,提高大豆種子丸?;陌沦|(zhì)量。Tamilselvi等[8]研究蘿卜種子的丸?;?,測試了蘿卜種子丸?;茡p率、效率及其裂解時間。Pasha等[9]通過EDEM模擬玉米種子在包衣裝置中的混合過程,不僅測試包衣的均勻性,還模擬玉米種子的包衣外層模型。李永祥等[10]通過采用EDEM中的Hertz-Mindlin with JKR接觸模型,模擬小麥粉螺旋給料裝置中的運動規(guī)律,并優(yōu)化給料裝置結(jié)構(gòu)。胡志超等[11]采用正交試驗方法,研究了花生玉米種子丸?;顑?yōu)的工藝參數(shù)。邵志威等[12]運用EDEM軟件,通過模擬牧草種子在振動下的運動規(guī)律,從而設(shè)計了牧草種子丸?;聶C(jī)模型。由于丸?;N子的粒徑大小不同,形狀不規(guī)則,在進(jìn)行種子丸?;倪^程中,種子與種子、種子與粉料、種粉與壁面之間的相互碰撞運動難以預(yù)測[13],因此一般的試驗裝置無法得出種子、粉料、壁面三者的混合運動規(guī)律,從而導(dǎo)致種子丸?;录夹g(shù)不成熟,無法解決種子包衣的質(zhì)量不高、合格率低、破碎率高等問題。

    本研究中在已有的丸?;芯炕A(chǔ)上設(shè)計丸粒機(jī)的勻盤結(jié)構(gòu),改變勻盤結(jié)構(gòu)與種粉之間的接觸面積的,以期達(dá)到改善混合效果的目的。采用正交實驗方法優(yōu)化玉米種子的丸?;に嚺c參數(shù),并基于統(tǒng)計學(xué)方法研究包衣機(jī)的勻盤結(jié)構(gòu)對顆粒包衣的影響規(guī)律,為解決包衣機(jī)的結(jié)構(gòu)與工藝設(shè)計問題提供參考。

    1理論基礎(chǔ)與包衣機(jī)制分析

    1.1DEM理論基礎(chǔ)

    根據(jù)牛頓第二定律,顆粒群的速度與位置隨著和時間的變化而改變,因此每個顆粒的控制方程為

    式中:mi為顆粒i的質(zhì)量;vi為顆粒i的移動速度;t為時間;Fij(n)顆粒i與顆粒j間的法向接觸力;Fij(t)為顆

    粒i與顆粒j間切向接觸力;g為重力加速度;Ii為顆粒i的轉(zhuǎn)動慣量;ω i為轉(zhuǎn)動角速度;Ri為顆粒i的半徑;μt為滾動摩擦系數(shù)。

    EDEM軟件仿真中最重要的環(huán)節(jié)是選擇合適的接觸模型,常用的接觸模型有Hertz-Mindlin、Hertz-Mindlin with JKR、Hysteretic Spring等。Johnson-Kendall-Roberts(JKR)接觸模型適用于描述顆粒間的黏結(jié)行為[14-15]。為此采用了JKR法向彈性接觸力FJKR來計算顆粒間的凝聚力,

    式中:FJKR為JKR法向彈性接觸力;E*為等效彈性模量;R*為等效接觸半徑;α為切向重疊量;γ為表面張力。等效彈性模量E*與等效接觸半徑R*定義為:

    (4)

    (5)

    式中:ν 1為顆粒1的泊松比;E 1為顆粒1的彈性模量;ν2為顆粒2的泊松比;E2為顆粒2的彈性模量;R 1為顆粒1的接觸半徑;R2為顆粒2的接觸半徑。

    1.2包衣過程理論分析

    包衣過程中,玉米種子與粉料落到勻盤后受到離心力作用,沿著包衣鍋的壁面作離心運動。在丸?;跗冢瑢⒂衩追N子與粉料的物理特性進(jìn)行簡化為顆粒,以便于包衣過程理論分析。勻盤以勻速進(jìn)行旋轉(zhuǎn),顆粒到達(dá)勻盤上后以相同的速度運動,忽略顆粒之間的作用力以及顆粒落到勻盤時的反彈。在丸?;衅?,顆粒在勻盤邊緣繞中心軸作圓周運動。選取其中的單個顆粒作為研究對象,并作為質(zhì)點,其受力分析如圖1所示。

    圖中:O為勻盤中心;E為種粉位置;ω為勻盤轉(zhuǎn)速,單位為r/min;G為顆粒重力,單位為N;Fr為顆粒的離心力,單位為N;Fn為勻盤對顆粒的支撐力,單位為N;Ff1為勻盤對顆粒的摩擦力,單位為N;FC為科里奧利力,單位為N。

    在玉米種子包衣過程前期,顆粒從勻盤中心向包衣鍋壁運動,則顆粒的離心力應(yīng)大于顆粒與勻盤的摩擦力,顆??朔c勻盤的摩擦向邊緣運動,假設(shè)顆粒運動到離勻盤中心軸的距離為r。對顆粒所受離心運動的動力學(xué)分析,可得

    F=Fr-Ff1-μ 1 FC=ma,(6)

    式中:F為顆粒所受合力;μ 1為顆粒與勻盤的摩擦系數(shù);a為顆粒運動的加速度;t為顆粒運動的時間。整理可得

    當(dāng)顆粒運動到勻盤邊緣時,即r=R,所需時間為tR,此時的顆粒的速度為

    式中,λ 1、λ2是特征方程λ2+2μ 1 λ2=0的根,則

    顆粒在勻盤Ⅰ的運動中,顆粒運動到離勻盤中心軸的距離r越遠(yuǎn),受到的離心作用力越大,與勻盤的接觸面積越大,如此反復(fù)地作圓周運動形成顆?;旌线_(dá)到丸粒化包衣效果。

    種粉在勻盤Ⅱ邊緣傾角、勻盤Ⅲ邊緣時,種粉的受力分析如圖2、3所示。

    圖中:R 1為勻盤的頂圓半徑,單位為mm;R2為勻盤的底圓半徑,單位為mm;h為勻盤的高度,單

    位為mm;β為勻盤邊緣傾角,單位為(°),種粉與勻盤之間的摩擦力分別為Ff1、Ff2,單位為N。如圖2、3受力分析可知

    mw2 r cos β-Ff1-mgsin β=0,F(xiàn)N-mgcosβ-mw2 r sin β=0。

    種粉在包衣鍋內(nèi)運動必須滿足勻盤對顆粒的摩擦力Ff1不小于顆粒與勻盤之間的摩擦系數(shù)μ 1與勻盤對顆粒的支撐力乘積FN,則勻盤轉(zhuǎn)速必須滿足

    (11)

    式(11)為勻盤轉(zhuǎn)速的臨界條件,當(dāng)勻盤轉(zhuǎn)速達(dá)到一定速度時,顆粒在受勻盤的離心力作用下沿著勻盤邊緣向上運動,顆粒到達(dá)包衣鍋側(cè)壁后受重力的作用落到勻盤上,如此反復(fù)運動以完成丸?;逻^程。玉米種子丸粒化包衣過程應(yīng)確保勻盤轉(zhuǎn)速大于式(11)的臨界值,否則會產(chǎn)生顆粒堆積。勻盤轉(zhuǎn)速的臨界條件與勻盤結(jié)構(gòu)有關(guān),勻盤的直徑、邊緣傾角等都會影響包衣的效果。

    2仿真模型建立與參數(shù)設(shè)置

    2.1顆粒快速填充

    對于EDEM軟件中復(fù)雜的顆粒模型,一般使用2種方法進(jìn)行顆粒建模,一是顆粒填充,二是API顆粒替換。顆粒填充更符合玉米種子外形,因此顆粒填充方法比較簡單,先將繪制好的種子三維模型先導(dǎo)入EDEM軟件,再在EDEM軟件里添加多個顆粒以填充滿種子模型即可。為了減少仿真計算時間,在保證玉米填充模型接近真實的前提下,要求填充的玉米顆粒盡量較少。如圖4所示,填充玉米種子顆粒數(shù)為33。

    2.2包衣鍋物理模型

    通過在三維繪圖軟件SolidWorks 2020中建立模型,簡化仿真對象,進(jìn)行全局參數(shù)設(shè)置,然后導(dǎo)入EDEM軟件中,具體模型如圖5、6所示。工作原理為:首先從包衣鍋鍋口倒入玉米種子,隨后倒入一定數(shù)量的包衣粉料。通過勻盤的轉(zhuǎn)動以及上下的振動,玉米種子與粉料進(jìn)行旋轉(zhuǎn)振動混合,由于濕潤的玉米種子和粉料都具有黏性,會在玉米種子表面形成包衣層,混合一定時間后包衣層達(dá)到穩(wěn)定,玉米種子丸?;^程結(jié)束。顆粒與包衣鍋的數(shù)據(jù)如表2、3所示。

    3離散元仿真

    3.1仿真過程分析

    本研究中主要是對用于玉米包衣種子的3種勻盤的形狀進(jìn)行研究,得出最有利于玉米種子包衣的勻盤結(jié)構(gòu)。根據(jù)顆??s放理論將離散元顆粒粒徑放大后進(jìn)行仿真,并不影響仿真結(jié)果[16]。為了減少仿真計算時間,炭與土壤的顆粒采用球形顆粒進(jìn)行仿真。建立粒徑為2.2 mm的土壤顆粒與2 mm炭顆粒,進(jìn)行仿真試驗。

    開展EDEM仿真試驗時,將玉米種子、土壤、炭的物理接觸參數(shù)、模型等導(dǎo)入,并設(shè)置勻盤的運動參數(shù)。在包衣鍋上方建立polygon平面并設(shè)置為顆粒工廠,玉米種子生成方式采用static,其他為dynamic方式,生成速率為10000 s-1。玉米種子、土、炭的顆粒數(shù)量分別為100、20000、10000。時間步長取Fixed Time Step的15%,為6.65×10-7,仿真總時間為10 s,每隔0.01 s記錄1次仿真數(shù)據(jù),網(wǎng)格尺寸為最小半徑的3倍。

    圖7—9為3種類型勻盤混合過程的效果圖。圖7(a)、8(a)、9(a)分別為丸?;跏茧A段,表明種粉自由落體到勻盤、受勻盤離心力的作用及種粉間的相互碰撞向包衣鍋壁面作圓周擴(kuò)散運動;圖7(b)、8(b)、9(b)所示為丸?;闹衅陔A段,分別表明種粉受勻盤的離心力作圓周運動、受勻盤振動作上下往復(fù)運動以及顆粒間的相互碰撞且沿著包衣鍋徑向作擴(kuò)散運動;圖7(c)、8(c)、9(c)所示為仿真的后期階段,分別表明顆粒自由落體到勻盤和受勻盤旋轉(zhuǎn)的離心力作用沿壁面作圓周運動。結(jié)合勻盤物理結(jié)構(gòu)和仿真結(jié)果分析,探究勻盤結(jié)構(gòu)與種粉運動對于丸?;伦鳂I(yè)質(zhì)量的影響規(guī)律,為小粒徑的丸粒化包衣機(jī)結(jié)構(gòu)設(shè)計和丸?;に噮?shù)優(yōu)化提供參考依據(jù)。

    3.2種粉間混合均勻度評價參數(shù)

    玉米種子包衣質(zhì)量均勻性是評價丸粒機(jī)丸?;鳂I(yè)效果的一個重要的指標(biāo)[17]。對于離散顆粒的混合均勻性,一般采用離散系數(shù)作為評價指標(biāo),本質(zhì)是通過標(biāo)準(zhǔn)偏差來反映混合顆粒的離散程度[18]。離散系數(shù)按下式計算。

    ×100%,(12)

    式中:CV為離散系數(shù);S為玉米樣本中包衣粉料含量的標(biāo)準(zhǔn)差;X為玉米樣本中包衣含量的平均值。通過離散系數(shù)CV計算混合物均勻度

    M=1-Cv。(13)

    在玉米種子的包衣的過程中,必須保證玉米種子不被破壞,將包衣的穩(wěn)定性系數(shù)作為參考值,進(jìn)行穩(wěn)定性離散系數(shù)的計算,得出穩(wěn)定性系數(shù)平均值。

    (14)

    標(biāo)準(zhǔn)差為

    (15)

    變異系數(shù)為

    ×100%,(16)

    式中:Xi為各測點的指標(biāo)值;X-為各測點指標(biāo)的平均值;N為測點數(shù);G為各測點指標(biāo)變異系數(shù)。

    為了獲得玉米種子丸粒化的離散系數(shù),在EDEM軟件后處理模塊中的Grid Bin Group劃分5×5×5網(wǎng)格,共計125個網(wǎng)格。

    4仿真結(jié)果與分析

    4.1轉(zhuǎn)速對玉米種子丸粒化的影響

    3種類型勻盤主要是通過改變顆粒與包衣鍋的接觸面積,對顆粒在包衣鍋內(nèi)的運動軌跡產(chǎn)生影響,從而影響玉米種子丸?;?。通過離散元EDEM軟件導(dǎo)入上述3種類型的勻盤模型,并設(shè)置勻盤轉(zhuǎn)速分別為30、40、50、60、70 r/min。圖10、11分別是玉米種子離散系數(shù)與玉米種子的配位數(shù)。其中配位數(shù)是指顆粒周圍固體顆粒的數(shù)量,也就是說一個顆粒周圍與其他相鄰的顆粒數(shù)目,配位數(shù)越大,顆粒的包衣效果越好。

    由圖10可知,勻盤的玉米種子離散系數(shù)在整個試驗轉(zhuǎn)速范圍內(nèi)呈現(xiàn)出不同的變化趨勢。勻盤Ⅰ在30~50 r/min時降幅度最大,50~70 r/min時離散系數(shù)增大;勻盤Ⅱ玉米種子離散系數(shù)在30~60 r/min時逐漸下降,超過60 r/min時離散系數(shù)變大;勻盤Ⅲ玉米種子離散系數(shù)在整個試驗轉(zhuǎn)速范圍內(nèi)逐步下降并趨于穩(wěn)定。圖11所示為3種類型勻盤在離散系數(shù)最小時玉米種子配位數(shù),勻盤Ⅰ與勻盤Ⅱ的玉米種子配位數(shù)比較接近,而勻盤Ⅲ玉米種子配位數(shù)明顯低于勻盤Ⅰ和勻盤Ⅱ的。

    4.2振動頻率對玉米種子丸?;挠绊?/p>

    為了研究振動頻率對玉米種子丸?;挠绊?,設(shè)置3種類型的勻盤轉(zhuǎn)速,分別50、60、70 r/min,在仿真開始時分別添加振動頻率0、5、10、15、20 Hz,振動幅度為6 mm的勻盤上下移動,分析振動頻率對玉米種子丸粒化的影響。如圖12、13分別為玉米種子離散系數(shù)與玉米種子配位數(shù)。

    由圖12可以看出,勻盤Ⅰ受振動頻率影響最為明顯,在振動頻率為0~10 Hz時,離散系數(shù)明顯減小,而超過10 Hz時,離散系數(shù)開始增大;勻盤Ⅱ在0~15 Hz時,離散系數(shù)逐步減小,之后開始增大,勻盤Ⅲ在0~20 Hz的范圍內(nèi),離散系數(shù)緩慢減小,并趨于平穩(wěn)的趨勢。從仿真過程中,添加振動頻率和振動幅度后有少許顆粒從包衣鍋溢出,因此選擇合適的振動頻率和振動幅度對玉米種子丸?;休^大的影響。由圖13可知,勻盤Ⅰ與勻盤Ⅱ的玉米種子配位數(shù)基本上相同,而勻盤Ⅲ的玉米種子配位數(shù)明顯小于勻盤Ⅰ和勻盤Ⅱ的,但較于圖7玉米種子配位數(shù)有明顯的增加。

    4.3振動幅度對玉米種子丸?;挠绊?/p>

    由上述振動頻率對玉米種子丸粒化的影響可知,振動頻率對玉米種子丸粒化有較大的影響。仿真中添加振動幅度為6 mm的丸粒化過程,有不少的顆粒從包衣鍋甩出,因此選擇合適的勻盤振動幅度對粉料利用率有很大的影響,同時也會增加玉米種子丸?;旌闲Ч?。在仿真過程中,分別添加振動幅度2、4、6、8 mm,3種勻盤的轉(zhuǎn)速分別為50、60、70 r/min,振動頻率為10、15、20 Hz。如圖14、15分別為玉米種子離散系數(shù)隨振幅的變化與玉米種子配位數(shù)隨時間的變化。

    由圖14所示,勻盤Ⅰ、Ⅱ在0~6 mm受勻盤振動幅度的影響,玉米種子離散系數(shù)下降最為明顯,在6~8 mm勻盤Ⅰ和勻盤Ⅱ種子離散系數(shù)上升,勻盤Ⅲ離散系數(shù)隨著振動幅度的上升逐漸減小。如圖15所示,3種類型的勻盤玉米種子配位數(shù)基本相同,說明勻盤Ⅲ在受轉(zhuǎn)速、振動頻率和振動幅度的情況可以達(dá)到勻盤Ⅰ和勻盤Ⅱ的丸?;笥衩追N子配位數(shù)。

    通過包衣勻盤的優(yōu)化后,玉米種子丸?;Ч辛孙@著的提升,勻盤Ⅰ在轉(zhuǎn)速為50 r/min、振動頻率為10 Hz、振動幅度為6 mm;勻盤Ⅱ的轉(zhuǎn)速為60 r/min、振動頻率為15 Hz、振動幅度為6 mm;勻盤Ⅲ的轉(zhuǎn)速為70 r/min、振動頻率為20 Hz、振動幅度為8 mm時它們的混合效果最優(yōu)化。

    4.4正交實驗

    通過觀察玉米種子丸粒化過程,可以得出,在不同的物理參數(shù)以及物理模型下,在初始階段,丸?;N粉之間的離散系數(shù)較大,在丸?;笃陔x散系數(shù)減少并趨于穩(wěn)定。玉米種子丸?;^程在一些結(jié)構(gòu)、參數(shù)的情況下會出現(xiàn)種子與粉料間的相互滲透作用較弱和種粉無法充分混合的問題,且種粉之間的接觸少、碰撞力過大、碰撞過少、種粉間容易發(fā)生黏結(jié)和團(tuán)聚等原因都會導(dǎo)致離散系數(shù)較大。通過上述的仿真過程及結(jié)果可知,在本文中,當(dāng)勻盤Ⅰ的轉(zhuǎn)速50 r/min、振動頻率為10 Hz、振動幅度為6 mm時離散系數(shù)最小,種粉混合效果最好。為進(jìn)一步優(yōu)化勻盤Ⅰ工況參數(shù),采用Box-Behnken原理設(shè)計試驗,將勻盤轉(zhuǎn)速、勻盤振動頻率、勻盤振動幅度設(shè)置為試驗因素,表4為試驗因素水平表。仿真試驗方案和仿真結(jié)果如表5所示。

    為了達(dá)到最好的玉米種子丸?;Ч?,必須減小離散系數(shù),通過模型優(yōu)化和Design-Expert 10軟件對回歸模型進(jìn)行優(yōu)化分析,離散系數(shù)結(jié)果分析如表6所示。模型的p值小于0.01,失擬項不明顯,說明該模型極為顯著,回歸方程與實際結(jié)果的擬合精度高,得到最佳的玉米種子丸?;瘏?shù)為使用勻盤Ⅰ,轉(zhuǎn)速為56.3 r/min,振動頻率為11.5 Hz,振動幅度為5.8 mm,離散系數(shù)為16.85%。

    5結(jié)論

    本文通過EDEM軟件建立玉米種子和包衣機(jī)的離散仿真模型,以包衣勻盤為研究對象,以離散系數(shù)作為評價指標(biāo)研究3種勻盤類型對均勻度的影響。

    1)使用Hertz-Mindlin with JKR模型進(jìn)行種粉丸?;抡?,通過種粉在3種勻盤結(jié)構(gòu)的受力分析,勻盤轉(zhuǎn)速臨界條件與勻盤結(jié)構(gòu)有關(guān),勻盤的直徑、邊緣傾角等都會影響包衣的效果。

    2)勻盤轉(zhuǎn)速、振動頻率、振動幅度對種粉間混合均勻度有影響,但是振動頻率、幅度過小對提升混合均勻效果不明顯,過大會使種粉顆粒從包衣鍋甩出。

    3)通過上述仿真分析得出,勻盤Ⅰ的轉(zhuǎn)速為56.3 r/min,振動頻率為11.5 Hz,振動幅度為5.8 mm,玉米種子丸?;Ч顬閮?yōu)。該研究結(jié)果可為玉米種子包衣提供技術(shù)依據(jù)與理論參考。

    利益沖突聲明(Conflict of Interests)

    所有作者聲明不存在利益沖突。

    All authors disclose no relevant conflict of interests.

    作者貢獻(xiàn)(Authors’Contributions)

    李紅軍、孫文杰進(jìn)行了方案設(shè)計、內(nèi)容總結(jié)與論文撰寫,孫文杰、陳偉進(jìn)行了實驗研究,張弛、張成俊、杜瑋參與了論文的審核與修改。所有作者均閱讀并同意了最終稿件的提交。

    The scheme design,content summary,and paper writing were carried out by LI Hongjun and SUN Wenjie. The experiment was conducted by SUN Wenjie and CHEN Wei. The manuscript was reviewed and revised by ZHANGChi,ZHANGChengjun,andDU Wei. Allauthors have read the final version of the paper and consented to its submission.

    參考文獻(xiàn)(References)

    [1]UTOMO HS,WENEFRIDA I,LIPSCOMBE SD. Use of an airplane and simulated aerial planting to evaluate seed coating,seed pelletization,and seeding rate on smooth cordgrass vegetation[J]. Ecological Restoration,2016,34(11):10-14.

    [2]DALSASSO SF,PIZARROA,MANFREDAS. Metrics for the quantification of seeding characteristics to enhance image velocimetry performance in rivers[J]. Remote Sensing,2020,12(11):1789.

    [3]ZENG D,LUO X,TU R. Application of bioactive coatings based on chitosan for soybean seed protection[J]. International Journal of Carbohydrate Chemistry,2012,11(12):22-43.

    [4]GE JT,GAN DF. MENGS C. Research status of seed coating and the need to implement good agricultural practices[J]. Seed,2016,35(2):45-49.

    [5]BAREKATAIN MR,WU SB,TOGHYANI M,et al. Effects of grinding and pelleting condition on efficiency of full-fat canola seed for replacing supplemental oil in broiler chicken diets[J]. Animal Feed Science and Technology,2015,207:140-149.

    [6]SIKHAO P,TAYLOR AG,MARINO ET. Development of seed agglomeration technology using lettuce and tomato as model vegetable crop seeds[J]. Scientia horticulturae,2015,184:85-92.

    [7]劉俏. 大豆種子包衣機(jī)種藥混合裝置的優(yōu)化設(shè)計與試驗研究[D]. 哈爾濱:東北農(nóng)業(yè)大學(xué),2020.

    LIUQ. Optimal design and experimental studyofseeddrugmixing device of soybean seed coating machine[D]. Harbin:Northeast Agricultural University,2020.

    [8]TAMILSELVIP,MANOHARJD. A study on physical properties of pelleted carrot seeds[J]. Advances in Life Sciences,2016,5(4):1220-1224.

    [9]PASHA M,HARE C,GHADIRI M,et al. Inter-particle coating variability in arotary batch seed coater[J]. Chemical Eng-ineering Researchamp;Design:Transactions of the Institution of Chemical Engineers,2017,120:92-101.

    [10]李永祥,李飛翔,徐雪萌. 基于顆??s放的小麥粉離散元參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(16):320-327.

    LI YX,LI FX,XU XM. Calibration of wheat flour discrete element parameters based on particle scaling[J]. Transactions of the Chinese Society of Agricultural Engineering,2019,35(16):320-327.

    [11]胡志超,田立佳,王海鷗. 種子丸?;O(shè)備的研制及丸化工藝的研究[J]. 農(nóng)機(jī)化研究,2008(3):95-97,101.

    HU ZC,TIAN LJ,WANG HO. Development of seed pelletization equipment and research on pelletization process[J]. Agricultural Mechanization Research,2008(3):95-97,101.

    [12]邵志威,陳智,侯占峰. 振動力場作用下牧草種子丸粒化包衣機(jī)試驗研究[J]. 農(nóng)機(jī)化研究,2018,40(8):124-128. SHAO ZW,CHEN Z,HOU ZF. Experimental study of forage seed pellet granulation coating machine under the action of vibration force field[J]. Agricultural Mechanization Research,2018,40(8):124-128.

    [13]趙永志,張憲旗,劉延雷. 滾筒內(nèi)非等粒徑二元顆粒體系增混機(jī)理研究[J]. 物理學(xué)報,2009,58(12):8386-8393. ZHAO YZ,ZHANG XQ,LIU YL. Study on the mixing mechanism of non-equal particle size binary particle system in the drum[J]. Acta Physica Sinica,2009,58(12):8386-8393.

    [14]TANG TQ,HE YR,REN AX,et al. Experimental study and DEM numerical simulation of dry wet particle flow behaviors in aspouted bed[J]. Industrialamp;Engineering Chemistry Research,2019,58(33):15353-15367.

    [15]MUGURUMA Y,TANAKA T,TSUJI Y. Numerical simulation of particulate flow with liquid bridge between particles[J]. Powder Technology,2000,109(1/2/3):49-57.

    [16]韓偉,王紹宗,張倩,等. 基于JKR接觸模型的微米級顆粒離散元參數(shù)標(biāo)定[J]. 中國粉體技術(shù),2021,27(6):60-69. HANW,WANGS Z,ZHANGQ,et al. Micron-scale particle discrete element parameter calibration based on the JKR contact model[J]. Chinese Science and Powder Technology,2021,27(6):60-69.

    [17]劉強(qiáng),許令峰,劉賢喜. 基于離散元的胡蘿卜機(jī)械拔取動態(tài)行為研究[J]. 中國農(nóng)機(jī)化學(xué)報,2018,39(6):61-65.

    LIU Q,XU LF,LIU XX. Dynamic behavior of carrot mechanical extraction based on discrete elements[J]. Chinese Journal of Agricultural Mechanization,2018,39(6):61-65.

    [18]DAI N,HOU Z,QIU Y. Numerical simulation and experiment of vibration pelletizer based on EDEM[J]. Inmateh-Agricultural Engineering,2021,63(1):469-478.

    Numerical simulation on pelletization of corn seeds with three uniform plate structures based on DEM

    LI Hongjun,SUN Wenjie,ZHANG Chi,ZHANG Chengjun,DU Wei,CHEN Wei

    School of Mechanical Engineering and Automation,Hubei Engineering Research Center of Industrial Detonator Intelligent Assembly,

    Wuhan Textile University,Wuhan 430200,China

    Abstract

    ObjectiveTo solve the problems of low coating quality and poor mixing uniformity of seed powder in the pelletization process of corn seeds.

    MethodsBased on the existing pelletization research,three platestructures were designed to achieve an optimal mixing effi?ciency by changing the contact area between the plate structure and seed powder. The study analyzed the pelletization mecha?nism of the three structures and simulated their pelletization process of corn seeds. A particle-filling method was adopted to make the model closer to the shape of actual corn seeds. To reduce simulation time,the number of corn seed filling particles wasreduced to 33. A simulation model of the corn seed pelletization process was established using the discrete element software(DEM)software EDEM,and the Hertz-Mindlin with Johnson-Kendall-Roberts(JKR)contact model was used to describe the bonding behavior between particles. The mixed motion of the three plate structures at the early,middle,and late stages of thesimulation was compared and analyzed. The effects of rotational speed,vibration frequency,and vibration amplitude on seed and powder mixing uniformity were studied using the dispersion coefficient as the evaluation index. The coating quality was evaluated based on the coordination number of the seeds,which represented the number of surrounding particles. A higher coor?dination number indicated better coating quality. The quadratic regression equation for corn seed dispersion coefficient was established using Design-Expert software for variance analysis. Rotational speed,vibration frequency,and vibration amplitude of uniform plate Iwere set at 50 r/min,10 Hz and 6 mm respectively. Using these test factors,17 sets of simulation tests were performed,and the resulting discrete coefficients were analyzed.

    Results and DiscussionThe results showed that the critical conditions for plate speed were related to the structure,diameter,and edge angle of the plate,all of which affected the coating quality. At the initial stage of pelletization,the seed powder wassubjected to centrifugal force and collided with the wall of the coating pan. At the middle stage,the seed powder underwent cir?cularmotion driven by centrifugal force,reciprocating up-and-down motion due to plate vibration,particle collisions,and radialdiffusion along the pan. At the final stage,particles experienced free fall and circular movement along the wall under centrifugalforce. The mixing uniformity was influenced by rotational speed,vibration frequency,and vibration amplitude. If either of themwas too low,the mixing effect would be poor. If the speed was too quick,the effect would be worse. Also,excessive vibrationfrequency and amplitude would cause seed powder to be thrown out of the pan. In some cases,the interaction between seeds andpowder was weak due to their insufficient contact and inappropriate collision force,leading to incomplete and high dispersioncoefficient. The optimal process parameters for rotational speed,vibration frequency,and vibration amplitude of the three platestructures were:(1)PlateI:50 r/min,10 Hz,and 6 mm;(2)Plate II:60 r/min,15 Hz,and 6 mm;(3)Plate III:70 r/min,20 Hz,and 8 mm. Using Design-Expert software,variance analysis on the quadratic regression equation showed that platespeed and vibration frequency significantly impacted mixing uniformity,while vibration amplitude had no significant effect.

    ConclusionThe study highlights that adjusting the contact area between the plate structure and seed powder significantly affects mixing uniformity. The particle-filling method tailored for irregular seed shapes reduces the simulation time while accurately approximating real particle behavior. The quadratic regression equation for corn seed dispersion coefficient provides agood fit with acorrected determination coefficient of 0. 95,indicating strong reliability. The lowest dispersion coefficient and the best mixing effect are achieved with plate Irotational speed of 56. 3 r/min,vibration frequency of 11. 5 Hz,and vibration amplitude of 5. 8 mm. The research results can provide references for the design of small particle granulators and the optimization of pellet?ization process parameters.

    Keywords:pelletization;discrete elements model;coating numerical simulatiovb

    (責(zé)任編輯:趙雁)

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預(yù)處理和加工研究報告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    亚洲欧美成人精品一区二区| 99热全是精品| 如何舔出高潮| 国产精品不卡视频一区二区| 婷婷精品国产亚洲av在线| 久久久久久久久大av| 美女xxoo啪啪120秒动态图| 99在线视频只有这里精品首页| 内地一区二区视频在线| 国产白丝娇喘喷水9色精品| 无遮挡黄片免费观看| 国产精品日韩av在线免费观看| 嫩草影院入口| 国产大屁股一区二区在线视频| 国产伦在线观看视频一区| 99热只有精品国产| 亚洲乱码一区二区免费版| 久久久久久久久大av| 亚洲最大成人av| 女人十人毛片免费观看3o分钟| 午夜福利视频1000在线观看| 可以在线观看的亚洲视频| 久久精品国产亚洲网站| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久噜噜| 久久久久免费精品人妻一区二区| 又粗又爽又猛毛片免费看| 九色成人免费人妻av| 亚洲av免费高清在线观看| 天堂影院成人在线观看| 精品不卡国产一区二区三区| 亚洲美女视频黄频| 久久韩国三级中文字幕| 久久久久精品国产欧美久久久| 白带黄色成豆腐渣| 我的女老师完整版在线观看| 99久国产av精品国产电影| 亚洲av成人精品一区久久| 亚洲人成网站在线播放欧美日韩| 免费电影在线观看免费观看| 18禁裸乳无遮挡免费网站照片| 国产乱人偷精品视频| 婷婷亚洲欧美| 日韩欧美三级三区| 一进一出抽搐动态| 日韩欧美在线乱码| 不卡一级毛片| 亚洲自拍偷在线| 久久久精品欧美日韩精品| 少妇人妻精品综合一区二区 | 欧美高清成人免费视频www| a级一级毛片免费在线观看| 日日干狠狠操夜夜爽| 精品人妻熟女av久视频| 久久久a久久爽久久v久久| 97超碰精品成人国产| 三级男女做爰猛烈吃奶摸视频| 成人毛片a级毛片在线播放| 老熟妇仑乱视频hdxx| 国产一区二区三区av在线 | 看免费成人av毛片| 日本与韩国留学比较| 色哟哟·www| 美女cb高潮喷水在线观看| 尤物成人国产欧美一区二区三区| 国产免费男女视频| 一进一出抽搐gif免费好疼| 久久精品夜夜夜夜夜久久蜜豆| av在线天堂中文字幕| 午夜爱爱视频在线播放| 中文字幕久久专区| 精品一区二区三区视频在线观看免费| 亚洲真实伦在线观看| 国产成年人精品一区二区| 欧美在线一区亚洲| 一级黄片播放器| 亚洲国产精品国产精品| 国产亚洲精品av在线| 久久热精品热| 国产高清激情床上av| 精品久久国产蜜桃| 欧美一区二区亚洲| 亚洲欧美中文字幕日韩二区| 精品久久久久久成人av| 99热精品在线国产| av在线亚洲专区| 国内精品一区二区在线观看| 永久网站在线| 欧美又色又爽又黄视频| 久久精品夜夜夜夜夜久久蜜豆| 变态另类丝袜制服| 欧美日韩乱码在线| 久久久久九九精品影院| 麻豆av噜噜一区二区三区| 国产在视频线在精品| 在线国产一区二区在线| 伦理电影大哥的女人| 午夜福利高清视频| 日韩欧美一区二区三区在线观看| 听说在线观看完整版免费高清| 男插女下体视频免费在线播放| 精品久久久久久久久久免费视频| 国产精品久久久久久亚洲av鲁大| 99国产精品一区二区蜜桃av| 久久久久久久亚洲中文字幕| 人人妻,人人澡人人爽秒播| 国产精品女同一区二区软件| 国产91av在线免费观看| 午夜福利在线观看免费完整高清在 | 少妇高潮的动态图| 免费高清视频大片| 国产69精品久久久久777片| 一个人看视频在线观看www免费| 18禁黄网站禁片免费观看直播| 卡戴珊不雅视频在线播放| 亚洲第一电影网av| 亚洲激情五月婷婷啪啪| 少妇猛男粗大的猛烈进出视频 | 成年女人看的毛片在线观看| 一级毛片电影观看 | 岛国在线免费视频观看| 午夜久久久久精精品| 国内精品美女久久久久久| 成人特级av手机在线观看| 久久久国产成人精品二区| 大香蕉久久网| 国产黄a三级三级三级人| 熟女人妻精品中文字幕| 婷婷亚洲欧美| 亚洲熟妇熟女久久| 国产成人a∨麻豆精品| 极品教师在线视频| 色哟哟哟哟哟哟| 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 大型黄色视频在线免费观看| 一个人观看的视频www高清免费观看| 天堂影院成人在线观看| 国产三级中文精品| 亚洲精品久久国产高清桃花| 亚洲欧美成人精品一区二区| 国产精品久久久久久精品电影| 毛片女人毛片| 91在线观看av| 精品乱码久久久久久99久播| 久久久国产成人精品二区| 波野结衣二区三区在线| 狂野欧美激情性xxxx在线观看| 成年女人毛片免费观看观看9| h日本视频在线播放| 亚洲18禁久久av| 一级黄片播放器| 色综合色国产| 一级毛片电影观看 | 国产成人a∨麻豆精品| 国产成人a区在线观看| 此物有八面人人有两片| 69av精品久久久久久| 又黄又爽又刺激的免费视频.| 国内少妇人妻偷人精品xxx网站| 午夜亚洲福利在线播放| 成人美女网站在线观看视频| av国产免费在线观看| 欧美zozozo另类| 国产精品久久久久久久久免| 免费大片18禁| 日韩欧美免费精品| 国产淫片久久久久久久久| 国产在线精品亚洲第一网站| 欧美人与善性xxx| 两个人的视频大全免费| 高清日韩中文字幕在线| 搡老妇女老女人老熟妇| 最近的中文字幕免费完整| 久久精品国产亚洲av天美| 97超碰精品成人国产| 日韩欧美 国产精品| 草草在线视频免费看| 久久综合国产亚洲精品| 可以在线观看的亚洲视频| 久久人妻av系列| 久久鲁丝午夜福利片| 人人妻,人人澡人人爽秒播| 十八禁网站免费在线| 国产成人a∨麻豆精品| 精品久久久久久成人av| 日韩,欧美,国产一区二区三区 | 在线免费十八禁| 中文字幕久久专区| ponron亚洲| 亚洲最大成人av| 亚洲一区高清亚洲精品| 国产探花极品一区二区| 深夜a级毛片| 午夜亚洲福利在线播放| 人妻少妇偷人精品九色| 啦啦啦韩国在线观看视频| 国产白丝娇喘喷水9色精品| eeuss影院久久| 超碰av人人做人人爽久久| 热99re8久久精品国产| 成人综合一区亚洲| 变态另类成人亚洲欧美熟女| 婷婷六月久久综合丁香| av卡一久久| 国产91av在线免费观看| 伊人久久精品亚洲午夜| 日韩一本色道免费dvd| 国产精品一区二区三区四区久久| 免费观看的影片在线观看| 国产色爽女视频免费观看| 国产熟女欧美一区二区| 国产成年人精品一区二区| 精品欧美国产一区二区三| 亚洲性久久影院| 99久久中文字幕三级久久日本| 桃色一区二区三区在线观看| 你懂的网址亚洲精品在线观看 | 床上黄色一级片| 精品99又大又爽又粗少妇毛片| 欧美一区二区精品小视频在线| 欧美性猛交╳xxx乱大交人| 国产精品伦人一区二区| 黄色配什么色好看| 久久久久久久久大av| 高清毛片免费观看视频网站| 国产精品电影一区二区三区| 久久欧美精品欧美久久欧美| av在线播放精品| 级片在线观看| 国产免费一级a男人的天堂| 国产av不卡久久| 国产精品美女特级片免费视频播放器| 在线天堂最新版资源| 国产黄色小视频在线观看| 精品久久久久久久久久免费视频| 国产不卡一卡二| 一进一出抽搐动态| 在线a可以看的网站| 国国产精品蜜臀av免费| 黄色欧美视频在线观看| 日日撸夜夜添| 国产精品久久电影中文字幕| 神马国产精品三级电影在线观看| 99在线视频只有这里精品首页| 国产精品一区二区三区四区免费观看 | 成年免费大片在线观看| 色哟哟哟哟哟哟| 一个人观看的视频www高清免费观看| 看黄色毛片网站| 99热网站在线观看| 99热这里只有精品一区| 五月伊人婷婷丁香| 欧美3d第一页| 一进一出抽搐gif免费好疼| 亚洲五月天丁香| 亚洲性久久影院| 99在线人妻在线中文字幕| 亚洲图色成人| 日本与韩国留学比较| 最近在线观看免费完整版| 色综合色国产| 亚洲av熟女| 在线国产一区二区在线| 国国产精品蜜臀av免费| 成人一区二区视频在线观看| 亚洲精品色激情综合| 午夜福利成人在线免费观看| 国产 一区精品| 国产高清视频在线观看网站| 一本精品99久久精品77| 搡女人真爽免费视频火全软件 | 成人毛片a级毛片在线播放| 欧美成人一区二区免费高清观看| 三级经典国产精品| 嫩草影院新地址| 成人精品一区二区免费| 亚洲熟妇熟女久久| 日韩大尺度精品在线看网址| 久久久成人免费电影| 美女高潮的动态| 国产爱豆传媒在线观看| 久久久国产成人免费| 黄色日韩在线| 免费人成视频x8x8入口观看| 国产探花极品一区二区| 亚洲av一区综合| 真人做人爱边吃奶动态| 国产国拍精品亚洲av在线观看| 99久久成人亚洲精品观看| 婷婷色综合大香蕉| 久久精品国产亚洲网站| av在线亚洲专区| 日日撸夜夜添| 久久久国产成人免费| 久久久久久久久久成人| 热99在线观看视频| 老司机影院成人| 一个人看视频在线观看www免费| 菩萨蛮人人尽说江南好唐韦庄 | 欧美日韩综合久久久久久| 精品久久久久久久久av| 人妻久久中文字幕网| 亚洲欧美精品自产自拍| 久久久精品94久久精品| 国产69精品久久久久777片| 亚洲人成网站在线播放欧美日韩| 九九热线精品视视频播放| 中出人妻视频一区二区| 国产熟女欧美一区二区| 免费看美女性在线毛片视频| 亚洲不卡免费看| 亚洲成人久久性| 淫秽高清视频在线观看| 亚洲精品456在线播放app| 欧美最新免费一区二区三区| 长腿黑丝高跟| 免费看日本二区| 成人亚洲欧美一区二区av| 久久人人爽人人爽人人片va| 卡戴珊不雅视频在线播放| 男人狂女人下面高潮的视频| 两性午夜刺激爽爽歪歪视频在线观看| 成人av一区二区三区在线看| 中国国产av一级| 国产精品永久免费网站| 不卡一级毛片| or卡值多少钱| 成人亚洲精品av一区二区| 小说图片视频综合网站| 色综合站精品国产| 老司机福利观看| 午夜福利在线观看免费完整高清在 | 午夜福利在线观看免费完整高清在 | 国内少妇人妻偷人精品xxx网站| 性色avwww在线观看| 色综合色国产| 免费人成视频x8x8入口观看| 色综合色国产| 久久韩国三级中文字幕| 精品福利观看| 国产精品久久久久久av不卡| 色综合色国产| 国产精品久久久久久av不卡| 一个人看的www免费观看视频| 校园人妻丝袜中文字幕| 麻豆一二三区av精品| 深夜精品福利| 精品一区二区三区av网在线观看| 精品一区二区免费观看| 精品久久久久久久久久久久久| 少妇丰满av| 免费看日本二区| 日本与韩国留学比较| 国产中年淑女户外野战色| 成人毛片a级毛片在线播放| 久久午夜亚洲精品久久| 日本黄色视频三级网站网址| 国产精华一区二区三区| 夜夜看夜夜爽夜夜摸| 97热精品久久久久久| 亚洲经典国产精华液单| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院| 国产精品女同一区二区软件| 国产免费男女视频| av在线播放精品| 不卡视频在线观看欧美| 最近中文字幕高清免费大全6| 亚洲无线观看免费| 在线观看美女被高潮喷水网站| 在线国产一区二区在线| 欧美性猛交╳xxx乱大交人| 亚洲人成网站在线观看播放| 嫩草影院精品99| 国产精品人妻久久久久久| 国产黄a三级三级三级人| 中文字幕久久专区| 亚洲内射少妇av| 成人综合一区亚洲| 色哟哟·www| 欧美潮喷喷水| 国产精品日韩av在线免费观看| 国产精品电影一区二区三区| 久久久久久久午夜电影| 免费人成在线观看视频色| АⅤ资源中文在线天堂| a级一级毛片免费在线观看| 色哟哟·www| 国产精品野战在线观看| 久久久久久大精品| 国产精品人妻久久久久久| 亚洲精品国产成人久久av| 又黄又爽又刺激的免费视频.| 草草在线视频免费看| 校园春色视频在线观看| 99热网站在线观看| 亚洲在线观看片| 日韩欧美一区二区三区在线观看| 国产aⅴ精品一区二区三区波| 99久久中文字幕三级久久日本| 又黄又爽又刺激的免费视频.| 波多野结衣高清作品| 欧美最黄视频在线播放免费| 久久久久久大精品| 国产成人freesex在线 | 天天躁夜夜躁狠狠久久av| 最近2019中文字幕mv第一页| 国产精品福利在线免费观看| 精品久久久久久成人av| 99热只有精品国产| 99久久精品一区二区三区| 淫秽高清视频在线观看| 国产成人freesex在线 | 日韩精品青青久久久久久| a级毛色黄片| 91av网一区二区| 内地一区二区视频在线| 亚洲一级一片aⅴ在线观看| 亚洲在线自拍视频| 91久久精品国产一区二区三区| 欧美日韩国产亚洲二区| 亚洲国产精品成人综合色| 国语自产精品视频在线第100页| 日韩高清综合在线| 一进一出抽搐gif免费好疼| 国产中年淑女户外野战色| 日韩欧美精品v在线| 欧美另类亚洲清纯唯美| 久久久久久久久久黄片| 真人做人爱边吃奶动态| 日本爱情动作片www.在线观看 | 国产真实乱freesex| 亚洲av中文av极速乱| 久久久国产成人精品二区| 久久精品夜夜夜夜夜久久蜜豆| 国产大屁股一区二区在线视频| 久久久久久久久久黄片| 在线国产一区二区在线| 欧美zozozo另类| 亚洲人成网站在线播放欧美日韩| 精品人妻熟女av久视频| 国产真实伦视频高清在线观看| 成熟少妇高潮喷水视频| 最新在线观看一区二区三区| 久久这里只有精品中国| 国产成人a区在线观看| 色综合站精品国产| 国内精品一区二区在线观看| 午夜福利视频1000在线观看| 18禁裸乳无遮挡免费网站照片| 国产aⅴ精品一区二区三区波| 成人二区视频| 国产精品一区二区三区四区免费观看 | 国产精品美女特级片免费视频播放器| 欧美又色又爽又黄视频| 亚洲欧美日韩高清专用| 自拍偷自拍亚洲精品老妇| 欧美极品一区二区三区四区| 久久精品国产清高在天天线| 午夜精品一区二区三区免费看| 久久久久久大精品| 99久久九九国产精品国产免费| 春色校园在线视频观看| 中文字幕av成人在线电影| 久久久国产成人精品二区| 久久婷婷人人爽人人干人人爱| 又黄又爽又刺激的免费视频.| 久久精品综合一区二区三区| 看免费成人av毛片| 亚洲最大成人手机在线| 免费看a级黄色片| 高清毛片免费看| 亚州av有码| 欧美高清成人免费视频www| 国产高清有码在线观看视频| 亚洲中文字幕日韩| 秋霞在线观看毛片| 男女做爰动态图高潮gif福利片| 亚洲人成网站高清观看| 欧美性猛交╳xxx乱大交人| 看十八女毛片水多多多| 国产真实乱freesex| 日韩欧美精品v在线| 精品人妻一区二区三区麻豆 | 国产综合懂色| 91在线精品国自产拍蜜月| 欧美区成人在线视频| 国产精品美女特级片免费视频播放器| 18禁黄网站禁片免费观看直播| 国产 一区 欧美 日韩| 久久久久久大精品| 一个人看的www免费观看视频| 亚洲欧美日韩高清专用| 亚洲国产精品sss在线观看| 中文在线观看免费www的网站| 一级毛片电影观看 | 亚洲图色成人| 国产精品日韩av在线免费观看| 人妻制服诱惑在线中文字幕| 国产精品嫩草影院av在线观看| 日韩精品青青久久久久久| 99热这里只有精品一区| 国产一区二区在线av高清观看| 亚洲最大成人中文| 国产欧美日韩一区二区精品| 亚洲成人中文字幕在线播放| 精品久久久噜噜| 免费观看人在逋| 国产视频一区二区在线看| 国产中年淑女户外野战色| 国产高清激情床上av| 国产av麻豆久久久久久久| 亚洲欧美日韩卡通动漫| 99久久久亚洲精品蜜臀av| 国产精品伦人一区二区| 精品午夜福利在线看| 国产伦在线观看视频一区| 插阴视频在线观看视频| 国产av在哪里看| 中文字幕人妻熟人妻熟丝袜美| 亚洲最大成人av| av天堂中文字幕网| 午夜激情福利司机影院| www.色视频.com| 赤兔流量卡办理| 美女内射精品一级片tv| 欧美又色又爽又黄视频| 一本精品99久久精品77| 18禁在线播放成人免费| 插逼视频在线观看| 欧美色欧美亚洲另类二区| 91麻豆精品激情在线观看国产| 亚洲中文日韩欧美视频| 亚洲真实伦在线观看| 欧美另类亚洲清纯唯美| 国国产精品蜜臀av免费| 国产黄片美女视频| 最近的中文字幕免费完整| 国产伦精品一区二区三区视频9| 亚洲丝袜综合中文字幕| 少妇猛男粗大的猛烈进出视频 | 一级毛片电影观看 | 久久久久久久亚洲中文字幕| 日本免费a在线| 色哟哟哟哟哟哟| 女的被弄到高潮叫床怎么办| 色哟哟·www| 黄色日韩在线| 国产日本99.免费观看| 色吧在线观看| 欧美绝顶高潮抽搐喷水| 日本欧美国产在线视频| 91av网一区二区| 免费人成在线观看视频色| 少妇猛男粗大的猛烈进出视频 | 美女xxoo啪啪120秒动态图| 你懂的网址亚洲精品在线观看 | 一级黄片播放器| 高清日韩中文字幕在线| 伦精品一区二区三区| 给我免费播放毛片高清在线观看| 日韩欧美精品v在线| 啦啦啦观看免费观看视频高清| 一区二区三区四区激情视频 | 噜噜噜噜噜久久久久久91| 亚洲人成网站在线播放欧美日韩| 成人鲁丝片一二三区免费| 国产一区亚洲一区在线观看| 亚洲国产精品合色在线| 国产成人aa在线观看| 丰满的人妻完整版| 尾随美女入室| 日韩亚洲欧美综合| 国产高清激情床上av| 欧美日本视频| 禁无遮挡网站| 国产不卡一卡二| 人人妻人人看人人澡| 日本一二三区视频观看| 亚洲精品成人久久久久久| 99久久精品一区二区三区| 国模一区二区三区四区视频| 无遮挡黄片免费观看| 欧美最黄视频在线播放免费| 性插视频无遮挡在线免费观看| 亚洲av中文av极速乱| 真实男女啪啪啪动态图| 女人十人毛片免费观看3o分钟| 精品熟女少妇av免费看| 在线看三级毛片| 国内揄拍国产精品人妻在线| 国内精品久久久久精免费| a级毛片a级免费在线| 亚洲av五月六月丁香网| 日本欧美国产在线视频| av天堂在线播放| 久久精品国产亚洲av天美| 成人一区二区视频在线观看| 男女啪啪激烈高潮av片| 男人狂女人下面高潮的视频| 久久久久久久久久黄片| 欧美高清成人免费视频www| 99国产精品一区二区蜜桃av| 欧美区成人在线视频| 最近2019中文字幕mv第一页| av国产免费在线观看| 午夜精品一区二区三区免费看| 午夜亚洲福利在线播放| 亚洲第一区二区三区不卡| 桃色一区二区三区在线观看| 日韩制服骚丝袜av| 欧美一级a爱片免费观看看| 激情 狠狠 欧美| 51国产日韩欧美|