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

    基于EDEM的茶葉揉捻機(jī)參數(shù)優(yōu)化及試驗(yàn)研究

    2020-06-15 10:55:38李兵李為寧柏宣丙黃劍虹
    茶葉科學(xué) 2020年3期
    關(guān)鍵詞:條率傾角茶葉

    李兵,李為寧,柏宣丙,黃劍虹

    基于EDEM的茶葉揉捻機(jī)參數(shù)優(yōu)化及試驗(yàn)研究

    李兵1,2,李為寧1,柏宣丙1,黃劍虹3

    1. 安徽農(nóng)業(yè)大學(xué)工學(xué)院,安徽 合肥 230036;2. 安徽農(nóng)業(yè)大學(xué)茶樹生物學(xué)與資源利用國(guó)家重點(diǎn)實(shí)驗(yàn)室,安徽 合肥 230036;3. 浙江上洋機(jī)械股份有限公司,浙江 衢州 324000

    為解決目前茶葉在揉捻過(guò)程中揉捻質(zhì)量不穩(wěn)定的問(wèn)題,以6CR-40型茶葉揉捻機(jī)為研究對(duì)象,基于Solidworks對(duì)茶葉揉捻機(jī)進(jìn)行三維建模,運(yùn)用離散元仿真軟件EDEM對(duì)茶葉揉捻機(jī)的揉捻過(guò)程進(jìn)行數(shù)值模擬,探求各試驗(yàn)因素對(duì)茶葉揉捻機(jī)性能指標(biāo)的影響規(guī)律。運(yùn)用二次正交旋轉(zhuǎn)試驗(yàn),通過(guò)Design-Expert軟件進(jìn)行優(yōu)化求解,獲得揉捻質(zhì)量最佳的結(jié)構(gòu)參數(shù)組合。結(jié)果表明,當(dāng)揉桶轉(zhuǎn)速為42?r·min-1、棱骨高度為10?mm、揉盤傾角為3.8°時(shí),茶葉的成條率為88.55%、碎茶率為1.83%。茶葉揉捻機(jī)具有較好的揉捻質(zhì)量。驗(yàn)證試驗(yàn)結(jié)果與仿真優(yōu)化結(jié)果基本一致。

    茶葉揉捻機(jī);EDEM;仿真分析;二次正交旋轉(zhuǎn)試驗(yàn);響應(yīng)面法

    揉捻是茶葉加工重要工序之一,茶葉加工過(guò)程中,除了白茶、綠茶和黃茶中的少類茶無(wú)需揉捻外,一般均需揉捻。殺青后的茶葉經(jīng)過(guò)揉捻工序使得茶葉體積縮小,卷緊條索,增加了茶葉容重,并且由于適當(dāng)破壞葉組織,茶汁逐漸溢出于茶葉外表面,加速多酚類化合物的酶促氧化,促進(jìn)了茶葉色香味的形成[1]?,F(xiàn)代茶葉加工生產(chǎn)過(guò)程中揉捻環(huán)節(jié)已經(jīng)實(shí)現(xiàn)機(jī)械化,主要依靠盤式茶葉揉捻機(jī)來(lái)完成,為提高茶葉揉捻機(jī)的工作性能,相關(guān)專家對(duì)茶葉揉捻機(jī)的揉捻性能進(jìn)行了大量研究。陳世輝[3]以茶葉揉捻機(jī)的外揉盤凹傾角為研究對(duì)象,通過(guò)有限元分析獲得揉捻特性較好時(shí)外揉盤的結(jié)構(gòu)參數(shù);岳鵬翔等[4]運(yùn)用機(jī)械優(yōu)化設(shè)計(jì)方法,建立優(yōu)化設(shè)計(jì)數(shù)學(xué)模型,對(duì)茶葉揉捻機(jī)的棱骨安裝數(shù)量、棱骨安裝角度等結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化設(shè)計(jì);周昊等[5]通過(guò)對(duì)棱骨的有限元?jiǎng)恿W(xué)分析,研究了不同棱骨結(jié)構(gòu)參數(shù)對(duì)揉捻葉的影響;瞿廷怡等[6]通過(guò)對(duì)揉捻機(jī)支臂優(yōu)化設(shè)計(jì)以減輕其重量,運(yùn)用電機(jī)進(jìn)行自動(dòng)加壓,實(shí)現(xiàn)揉捻自動(dòng)化。上述研究對(duì)揉捻機(jī)的參數(shù)優(yōu)化提供了新的方法,但并未能對(duì)揉捻過(guò)程中茶葉顆粒的運(yùn)動(dòng)學(xué)參數(shù)進(jìn)行量化分析。近年來(lái),隨著計(jì)算機(jī)技術(shù)不斷發(fā)展,離散元法及其虛擬仿真在農(nóng)業(yè)工程領(lǐng)域有著越來(lái)越多的應(yīng)用,為研究物料顆粒—機(jī)械零部件間相互作用的變化規(guī)律,優(yōu)化機(jī)械結(jié)構(gòu)提供了新的技術(shù)手段[7-10]。

    本文以6CR-40型茶葉揉捻機(jī)為研究載體,通過(guò)Solidworks建立茶葉揉捻機(jī)三維模型,運(yùn)用EDEM對(duì)揉捻過(guò)程中茶葉運(yùn)動(dòng)學(xué)參數(shù)進(jìn)行量化分析,得到不同結(jié)構(gòu)參數(shù)對(duì)茶葉揉捻質(zhì)量的影響規(guī)律,以成條率和碎茶率作為性能評(píng)價(jià)指標(biāo),進(jìn)行二次正交旋轉(zhuǎn)組合試驗(yàn)并借助Design-Expert繪制影響茶葉揉捻機(jī)性能的響應(yīng)曲面圖,旨在獲得茶葉揉捻機(jī)的最優(yōu)參數(shù)組合,為該型號(hào)揉捻機(jī)的優(yōu)化設(shè)計(jì)提供參考依據(jù)。

    1 茶葉揉捻機(jī)的基本結(jié)構(gòu)與工作原理

    茶葉揉捻機(jī)由電動(dòng)機(jī)、小皮帶輪、大皮帶輪、減速器、曲柄、立柱、加壓臂、加壓蓋、揉桶、揉盤、三角架、支撐機(jī)架、出茶門推桿、加壓手輪以及手柄等組成[11]。茶葉揉捻機(jī)的相關(guān)結(jié)構(gòu)設(shè)計(jì)參數(shù)如表1所示。揉盤是一個(gè)由四周向中心位置逐漸下凹的盤型結(jié)構(gòu)件,通常由不銹鋼板做成,并裝有若干不銹鋼弧形棱骨,用以增強(qiáng)茶葉揉捻的揉搓力。揉捻過(guò)程中壓力調(diào)整通過(guò)絲桿和螺母的運(yùn)動(dòng)控制加壓蓋的升降來(lái)實(shí)現(xiàn),加壓蓋的上、下浮動(dòng)由彈簧控制。揉盤的中心位置是用于開啟和關(guān)閉茶葉出料的出茶門,出茶門采用擺動(dòng)式結(jié)構(gòu),由出茶門推桿、手柄等組成。當(dāng)轉(zhuǎn)動(dòng)手柄時(shí)出茶門隨之?dāng)[動(dòng)回轉(zhuǎn),若將出茶門關(guān)閉,將出茶門推桿鎖住出茶門即可。茶葉揉捻機(jī)整體結(jié)構(gòu)圖,如圖1所示。

    茶葉揉捻機(jī)工作時(shí),電動(dòng)機(jī)作為動(dòng)力輸出,經(jīng)皮帶傳動(dòng)及齒輪箱轉(zhuǎn)換運(yùn)動(dòng)方向,驅(qū)動(dòng)主曲柄做回轉(zhuǎn)運(yùn)動(dòng),帶動(dòng)揉桶在揉盤上做偏心回轉(zhuǎn)運(yùn)動(dòng),茶葉在揉桶中受加壓蓋加壓,并在棱骨的搓揉作用下進(jìn)行揉捻作業(yè)。揉捻過(guò)程大致可以分為3個(gè)階段,分別是空壓起條階段、加壓緊條階段以及松壓解團(tuán)階段。空壓階段是促使茶葉成條的基礎(chǔ)階段,空壓可以減少茶葉的滾動(dòng)摩擦阻力,便于片狀葉起條。加壓階段增大了茶葉顆粒間的內(nèi)摩擦力以及棱骨對(duì)茶團(tuán)的揉搓作用,促使茶葉細(xì)胞破壞,條索卷緊。松壓階段,加壓蓋逐漸上移致使葉層占有的空間增大,緊壓的茶團(tuán)被振松抖散。

    表1 茶葉揉捻機(jī)的結(jié)構(gòu)設(shè)計(jì)參數(shù)

    注:1-電動(dòng)機(jī),2-小皮帶輪,3-減速器,4-大皮帶輪,5-曲柄,6-立柱,7-加壓臂,8-加壓蓋,9-揉桶,10-三腳架,11-揉盤,12-支撐機(jī)架,13-出茶門推桿,14-加壓手輪,15-手柄

    2 茶葉揉捻運(yùn)動(dòng)規(guī)律

    茶葉在揉捻過(guò)程中的運(yùn)動(dòng)情況十分復(fù)雜,不僅有揉桶、揉盤等結(jié)構(gòu)部件對(duì)茶葉的翻滾運(yùn)動(dòng),而且還有茶葉之間的相互運(yùn)動(dòng),以及茶葉散落在揉盤上的自由運(yùn)動(dòng)。揉捻葉受到多種力的綜合作用而形成一定的運(yùn)動(dòng)規(guī)律,將揉桶內(nèi)所有茶葉看成一個(gè)整體,茶團(tuán)所受的作用力包括:揉桶側(cè)壁的推力、揉盤支持力、加壓蓋正壓力、茶葉自身重力以及茶葉顆粒間的內(nèi)摩擦力[12]。以單顆茶葉顆粒為例進(jìn)行受力分析,如圖2所示。

    1=2/····························· (1)

    F=(+3) ························· (2)

    =+4································(4)

    式中,1—揉桶側(cè)壁推力,N;3—茶葉顆粒所受擠壓力,N;—茶葉顆粒單位面積所受擠壓力,N;4—加壓蓋壓力,N;—茶葉顆粒間內(nèi)摩擦力,N;—加壓蓋壓力與茶葉重力之和,N;—茶葉自身重力,N;—茶葉顆粒的質(zhì)量,kg;—揉桶半徑,m;—茶葉顆粒的線速度,m·s-1。

    揉捻葉在揉桶內(nèi)不僅隨揉桶中心作圓周運(yùn)動(dòng),而且又在翻轉(zhuǎn)作用力下向上前方運(yùn)動(dòng),最后散落在揉盤上。每一時(shí)刻揉捻葉在揉桶內(nèi)的位置都不相同,故形成了茶葉揉捻運(yùn)動(dòng)規(guī)律,可大致分為4個(gè)不同的作用區(qū)[13]??拷啾P和棱骨的區(qū)域由于揉搓作用力較大,茶葉在揉桶推力作用下向前運(yùn)動(dòng)形成揉搓區(qū)。在揉搓區(qū)內(nèi)揉捻葉運(yùn)動(dòng)速度較快,從而將茶葉揉捻成條,體積逐漸縮?。辉谌嗤?、揉盤以及加壓蓋對(duì)茶葉作用力的交點(diǎn)周圍,由于茶葉受到的作用力最大,故形成強(qiáng)壓區(qū)。當(dāng)茶葉進(jìn)入強(qiáng)壓區(qū)時(shí),茶葉被推向揉桶中心位置匯集,相互擠壓成團(tuán)狀,并在翻轉(zhuǎn)作用力下向上翻轉(zhuǎn);由于揉桶上半部揉捻葉受到的摩擦阻力較小,故在慣性力和茶葉自身重力的作用下而向前下翻,形成下翻區(qū);揉桶繼續(xù)作回轉(zhuǎn)運(yùn)動(dòng),下翻區(qū)的揉捻葉散落在揉盤尚未進(jìn)行揉捻的區(qū)域,形成散落區(qū)。

    3 仿真過(guò)程及分析

    3.1 仿真模型建立

    為了能夠更加清晰、直觀地模擬茶葉的揉捻過(guò)程,運(yùn)用三維建模軟件Solidworks對(duì)揉桶、揉盤等關(guān)鍵結(jié)構(gòu)部件實(shí)體建模(圖3),并導(dǎo)入EDEM中進(jìn)行仿真分析。

    茶葉顆粒在揉捻過(guò)程中形狀差異明顯,按照長(zhǎng)寬高的正態(tài)分布尺寸選取建模對(duì)象,故在EDEM中擬合出與茶葉顆粒外形輪廓相吻合的多球面組合體作為顆粒模型[14],按照茶葉大小建立茶葉顆粒仿真模型,如圖4所示。根據(jù)茶葉物料物理特性,設(shè)置茶葉顆粒模型的物理屬性參數(shù)和接觸屬性參數(shù)[15-18],如表2所示。

    根據(jù)茶葉顆粒表面物理特性,選擇Hertz-Mindlin(no slip)模型作為仿真接觸模型。設(shè)置EDEM顆粒工廠以300個(gè)·s-1的速率生成茶葉顆粒模型,保證揉桶內(nèi)有充足茶葉顆粒進(jìn)行仿真模擬。

    注:F1-揉桶側(cè)壁推力,F(xiàn)2—揉盤支持力,F(xiàn)4—加壓蓋壓力,F(xiàn)f—茶葉顆粒間內(nèi)摩擦力,G—茶葉自身重力

    圖3 茶葉揉捻機(jī)仿真模型

    圖4 茶葉顆粒模型

    表2 茶葉顆粒參數(shù)

    3.2 揉桶轉(zhuǎn)速對(duì)茶葉揉捻質(zhì)量的影響

    揉桶轉(zhuǎn)速對(duì)茶葉顆粒在揉桶中的運(yùn)動(dòng)有著顯著影響。隨著揉桶轉(zhuǎn)速增大,揉捻葉的平均速度明顯增大(圖5-A)。揉桶轉(zhuǎn)速對(duì)茶葉顆粒受力影響分為兩個(gè)階段:0~3.5?s是松散茶葉揉捻成茶團(tuán)的過(guò)程,此過(guò)程揉桶轉(zhuǎn)速對(duì)茶葉顆粒受力影響不大;3.5?s后茶團(tuán)形成,揉桶轉(zhuǎn)速對(duì)茶葉顆粒受力影響較大,當(dāng)揉桶轉(zhuǎn)速適當(dāng)時(shí)(45?r·min-1),茶葉顆粒受力較為均勻,有助于茶葉揉搓成條狀。揉桶轉(zhuǎn)速過(guò)低,機(jī)器對(duì)茶葉的擠壓、揉搓作用力過(guò)小,茶葉在揉桶中不能充分翻動(dòng),導(dǎo)致所制茶葉條索粗松,茶葉顆粒受力不均勻,成條率低;揉桶轉(zhuǎn)速過(guò)高,茶葉顆粒間內(nèi)摩擦力過(guò)大,碎茶率相應(yīng)增大。當(dāng)產(chǎn)生較多數(shù)量的碎茶時(shí)內(nèi)摩擦力減小,成條率降低,同時(shí)造成機(jī)器較大振動(dòng)(圖5-B)。

    3.3 棱骨高度對(duì)茶葉揉捻質(zhì)量的影響

    理論分析及實(shí)際使用均表明,棱骨在茶葉揉捻成形過(guò)程當(dāng)中發(fā)揮著重要作用,棱骨結(jié)構(gòu)有利于揉捻過(guò)程中形成茶團(tuán)、促進(jìn)其向上翻轉(zhuǎn)并扭轉(zhuǎn)成條[19]。仿真結(jié)果表明:棱骨高度變化對(duì)茶葉顆粒的平均速度影響不大(圖6-A);但棱骨高度變化對(duì)揉捻葉的作用力影響較大,棱骨高度減小,棱骨對(duì)茶葉的摩擦阻力以及茶葉向上翻轉(zhuǎn)的作用力減小,導(dǎo)致揉捻葉成條率降低。同時(shí)棱骨不能把茶葉導(dǎo)向揉捻強(qiáng)壓區(qū),揉捻葉的跑茶率相應(yīng)增大;棱骨高度過(guò)高,棱骨對(duì)揉捻葉的摩擦阻力增大,位于揉盤底部的茶葉易揉碎,碎茶率增大,揉捻質(zhì)量下降(圖6-B)。

    圖6 不同棱骨高度下茶葉顆粒的平均速度、作用力曲線

    3.4 揉盤傾角對(duì)茶葉揉捻質(zhì)量的影響

    揉盤對(duì)揉捻茶團(tuán)的支持力是偏向揉盤中心的,它是促使茶團(tuán)在揉盤翻滾的重要因素之一。如圖7-A所示,0~0.8?s,為松散茶葉加速時(shí)間,松散茶葉在揉桶轉(zhuǎn)動(dòng)帶動(dòng)下逐漸加速;0.8~3.0?s,為茶葉成團(tuán)過(guò)程,茶葉由于茶葉間互相摩擦,粘滯阻力及內(nèi)摩擦力逐漸加大,茶葉速度下降并逐步揉捻成團(tuán);3.0~3.5?s,為茶團(tuán)加速過(guò)程,成團(tuán)的茶葉在揉桶的帶動(dòng)下速度逐漸增大到最大值。揉盤傾角對(duì)茶團(tuán)的形成影響較小,但隨著揉盤傾角變大,茶團(tuán)的最大速度也增大。如圖7-B所示,所受作用力在揉盤傾角為4°時(shí)明顯增大且達(dá)到最大值。揉盤傾角過(guò)大使得強(qiáng)壓區(qū)的揉搓時(shí)間過(guò)長(zhǎng),揉桶側(cè)邊與揉盤中心位置的揉捻葉不容易發(fā)生易位,形成局部強(qiáng)壓區(qū),從而導(dǎo)致揉捻葉碎茶率相應(yīng)增大;揉盤傾角過(guò)小使得揉捻葉不容易翻滾,茶葉不容易進(jìn)入強(qiáng)壓區(qū)揉搓成形,從而造成揉捻葉成條率降低,揉捻質(zhì)量較差。

    4 正交旋轉(zhuǎn)試驗(yàn)與結(jié)果

    4.1 試驗(yàn)因素與評(píng)價(jià)指標(biāo)

    為了研究揉捻過(guò)程中影響茶葉揉捻質(zhì)量的主要因素,得到揉捻性能最優(yōu)的結(jié)構(gòu)參數(shù)組合,在仿真模型建立及EDEM結(jié)果分析基礎(chǔ)上,進(jìn)行二次正交旋轉(zhuǎn)試驗(yàn)。影響茶葉揉捻機(jī)揉捻質(zhì)量?jī)?yōu)劣的主要因素有揉桶轉(zhuǎn)速、揉盤傾角、棱骨高度、棱骨數(shù)量等多種因素,岳鵬翔等[4]對(duì)棱骨數(shù)量?jī)?yōu)化結(jié)果為12,為了探求其他影響因素的優(yōu)化結(jié)果,選取揉桶轉(zhuǎn)速、棱骨高度以及揉盤傾角作為試驗(yàn)因素,選擇成條率和碎茶率作為揉捻機(jī)的性能評(píng)價(jià)指標(biāo),設(shè)計(jì)3因素3水平的二次正交旋轉(zhuǎn)試驗(yàn),試驗(yàn)因素編碼如表3所示,每組試驗(yàn)重復(fù)3次,取平均值,試驗(yàn)結(jié)果如表4所示。

    表3 試驗(yàn)因素編碼

    4.2 試驗(yàn)原料和設(shè)備

    試驗(yàn)原料采用6CST-60型電熱滾筒殺青機(jī)進(jìn)行殺青處理,然后將殺青葉混勻后,通過(guò)更換不同結(jié)構(gòu)參數(shù)(棱骨高度、揉盤傾角)的揉盤進(jìn)行茶葉的揉捻試驗(yàn)(圖8)。加入適量殺青葉,揉捻時(shí)間為35?min,采用常規(guī)加壓方式[20]。揉捻過(guò)程結(jié)束后對(duì)揉捻葉取樣,測(cè)定所取樣品的成條率和碎茶率。為了使試驗(yàn)結(jié)果能正確反映茶葉揉捻機(jī)的揉捻性能,按照茶葉揉捻機(jī)行業(yè)標(biāo)準(zhǔn)(JB/T 9814—2007)將成條率和碎茶率指標(biāo)的計(jì)算方法規(guī)定如下:

    成條率:取揉捻葉樣50?g,剔除老梗老葉、雜質(zhì),分出該成條而未成條葉(包括扁條、扁片)、碎茶和成條葉3種,分別稱重,成條率計(jì)算公式如下:

    式中,—成條率,%;W—成條質(zhì)量,g;∑—成條葉質(zhì)量、該成條而未成條葉質(zhì)量、碎茶質(zhì)量之和,g;—焦末碎片率,%。

    碎茶率:取恒溫箱烘干揉捻葉樣200?g(含水率為7%~9%),分兩次取樣,每次取樣100?g,在轉(zhuǎn)速為200?r·min-1碎末茶測(cè)定計(jì)上用直徑為280?cm的16目樣篩篩動(dòng)5轉(zhuǎn),稱篩下的碎茶質(zhì)量,碎茶率計(jì)算公式如下:

    式中,—碎茶率,%;s—碎茶質(zhì)量,g;—焦末碎片率,%。

    4.3 結(jié)果與分析

    運(yùn)用Design-Expert 10.0對(duì)表4中的試驗(yàn)數(shù)據(jù)進(jìn)行回歸分析,并對(duì)影響評(píng)價(jià)指標(biāo)的因素進(jìn)行顯著性分析[21],得到成條率1、碎茶率2與揉桶轉(zhuǎn)速、棱骨高度以及揉盤傾角之間的回歸方程:

    1=–14.755?75+3.143?93+4.936?88+2.654?62-0.008?875-0.007?5+0.069?375-0.032?6952-0.225?52-0.389?82,

    2=2.724-0.009?35-0.247?5-0.226?75+0.000?25+0.001+0.027?5+0.000?492-0.009?752+0.006?52。

    由表5可知,茶葉揉捻機(jī)成條率的回歸模型的<0.000?1,影響極顯著;失擬項(xiàng)=0.135?5(>0.05),影響不顯著,模型的決定系數(shù)2=0.991?0。茶葉揉捻機(jī)碎茶率的回歸模型<0.000?1,影響極顯著;失擬項(xiàng)=0.097?6(>0.05),影響不顯著,模型的決定系數(shù)2=0.981?2。即表明在一定范圍內(nèi),回歸模型與實(shí)際情況擬合度較高,無(wú)失擬因素存在,故可用成條率和碎茶率的回歸方程代替試驗(yàn)真實(shí)點(diǎn)對(duì)茶葉揉捻質(zhì)量進(jìn)行分析與預(yù)測(cè)[22]。

    圖8 不同結(jié)構(gòu)參數(shù)下茶葉揉捻機(jī)揉捻試驗(yàn)

    表4 試驗(yàn)設(shè)計(jì)與結(jié)果

    表5 回歸方程方差分析

    為了更加直觀的分析試驗(yàn)因素與評(píng)價(jià)指標(biāo)之間的關(guān)系,運(yùn)用Design-Expert 10.0軟件得到不同試驗(yàn)因素與評(píng)價(jià)指標(biāo)的響應(yīng)曲面圖[23]。

    根據(jù)回歸方程和響應(yīng)曲面可得,揉桶轉(zhuǎn)速、棱骨高度以及揉盤傾角間存在交互作用,并且對(duì)成條率的影響較為顯著。當(dāng)揉桶轉(zhuǎn)速一定時(shí),成條率隨著棱骨高度增大而增大;當(dāng)棱骨高度一定時(shí),成條率隨著揉桶轉(zhuǎn)速增大先增大后減?。划?dāng)揉桶轉(zhuǎn)速變化時(shí),成條率的變化區(qū)間較大,揉桶轉(zhuǎn)速對(duì)成條率指標(biāo)的影響更為顯著(圖9-A)。當(dāng)揉盤傾角一定時(shí),成條率隨著揉桶轉(zhuǎn)速增大先增大后減??;當(dāng)揉桶轉(zhuǎn)速一定時(shí),成條率隨著揉盤傾角增大先增大后減??;當(dāng)揉盤轉(zhuǎn)速變化時(shí),成條率的變化區(qū)間較大,揉盤轉(zhuǎn)速對(duì)成條率指標(biāo)的影響更為顯著(圖9-B)。當(dāng)棱骨高度一定時(shí),成條率隨著揉盤傾角增大先增大后減??;當(dāng)揉盤傾角一定時(shí),成條率隨著棱骨高度增大而增大。當(dāng)棱骨高度變化時(shí),成條率的變化區(qū)間較大,棱骨高度對(duì)成條率指標(biāo)的影響更為顯著(圖9-C)。對(duì)成條率指標(biāo)影響顯著性的由大至小依次為:揉桶轉(zhuǎn)速、棱骨高度、揉盤傾角。

    圖9 各因素對(duì)成條率的響應(yīng)曲面

    圖10 各因素對(duì)碎茶率的響應(yīng)曲面

    當(dāng)揉桶轉(zhuǎn)速一定時(shí),碎茶率隨著棱骨高度增大而增大。當(dāng)棱骨高度一定時(shí),碎茶率隨著揉桶轉(zhuǎn)速增大而增大;當(dāng)揉桶轉(zhuǎn)速變化時(shí),碎茶率指標(biāo)的變化區(qū)間較大,揉桶轉(zhuǎn)速對(duì)碎茶率指標(biāo)的影響更為顯著(圖10-A)。當(dāng)揉盤傾角一定時(shí),碎茶率隨著揉桶轉(zhuǎn)速的增大而增大;當(dāng)揉桶轉(zhuǎn)速一定時(shí),碎茶率隨著揉盤傾角的增大而增大;當(dāng)揉桶轉(zhuǎn)速變化時(shí),碎茶率指標(biāo)的變化區(qū)間較大,揉桶轉(zhuǎn)速對(duì)碎茶率指標(biāo)的影響更為顯著(圖10-B)。當(dāng)棱骨高度一定時(shí),碎茶率隨著揉盤傾角增大而增大;當(dāng)揉盤傾角一定時(shí),碎茶率隨著棱骨高度增大而增大;當(dāng)棱骨高度變化時(shí),碎茶率的變化區(qū)間顯著,棱骨高度對(duì)碎茶率指標(biāo)的影響更為顯著(圖10-C)。對(duì)碎茶率指標(biāo)影響顯著性由大至小依次為:揉桶轉(zhuǎn)速、棱骨高度、揉盤傾角。

    為了獲得茶葉揉捻機(jī)性能的最優(yōu)參數(shù)組合,對(duì)試驗(yàn)進(jìn)行優(yōu)化設(shè)計(jì)。結(jié)合試驗(yàn)因素的邊界條件,建立茶葉揉捻機(jī)的優(yōu)化數(shù)學(xué)模型,對(duì)成條率和碎茶率的回歸方程進(jìn)行分析。運(yùn)用Design-Expert 10.0軟件中的參數(shù)優(yōu)化模塊對(duì)優(yōu)化數(shù)學(xué)模型求解[24],得到優(yōu)化結(jié)果:揉桶轉(zhuǎn)速為42?r·min-1,棱骨高度為10?mm,揉盤傾角為3.8°時(shí),成條率指標(biāo)為88.55%,碎茶率指標(biāo)為1.83%,茶葉揉捻機(jī)的揉捻質(zhì)量最優(yōu)。

    為了進(jìn)一步驗(yàn)證參數(shù)優(yōu)化結(jié)果的正確性,在安徽農(nóng)業(yè)大學(xué)茶葉加工中心進(jìn)行茶葉揉捻機(jī)的驗(yàn)證試驗(yàn)。以6CR-40型茶葉揉捻機(jī)為驗(yàn)證試驗(yàn)樣機(jī),選取揉桶轉(zhuǎn)速為42?r·min-1,棱骨高度為10?mm,揉盤傾角為3.8°的試驗(yàn)樣機(jī)進(jìn)行驗(yàn)證試驗(yàn),加入適量的殺青葉,揉捻時(shí)間20?min,采用常規(guī)加壓方式,得到茶葉揉捻機(jī)的成條率為88.32%,碎茶率為1.87%。參數(shù)優(yōu)化結(jié)果與驗(yàn)證試驗(yàn)結(jié)果基本一致。

    5 小結(jié)

    通過(guò)Solidworks建立揉捻機(jī)三維模型,基于離散元法建立茶葉顆粒仿真模型,運(yùn)用EDEM軟件對(duì)揉捻過(guò)程進(jìn)行數(shù)值模擬,解決了茶葉加工過(guò)程中的運(yùn)動(dòng)學(xué)量化分析問(wèn)題,以成條率和碎茶率為目標(biāo)函數(shù),以揉桶轉(zhuǎn)速、棱骨高度、揉盤傾角為試驗(yàn)因素,對(duì)茶葉揉捻機(jī)相關(guān)參數(shù)進(jìn)行優(yōu)化設(shè)計(jì)。運(yùn)用二次正交旋轉(zhuǎn)組合試驗(yàn)進(jìn)行數(shù)據(jù)處理及驗(yàn)證,利用Design-Expert軟件優(yōu)化模塊對(duì)試驗(yàn)結(jié)果進(jìn)行優(yōu)化求解,正交旋轉(zhuǎn)試驗(yàn)優(yōu)化結(jié)果與驗(yàn)證試驗(yàn)結(jié)果基本一致。因此運(yùn)用離散元法對(duì)茶葉揉捻機(jī)的結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化具有可行性,并且為茶葉揉捻機(jī)的優(yōu)化設(shè)計(jì)提供了一種新方法。本文僅選取茶葉揉捻機(jī)關(guān)鍵結(jié)構(gòu)參數(shù)當(dāng)中的揉桶轉(zhuǎn)速、棱骨高度以及揉盤傾角對(duì)揉捻質(zhì)量進(jìn)行分析研究,在今后的研究中為了更加精確建模,可綜合考慮棱骨寬度、棱骨截面形狀、棱骨螺旋角、棱骨數(shù)量等其他結(jié)構(gòu)參數(shù)以及茶葉物理特性參數(shù)對(duì)茶葉揉捻質(zhì)量的影響。

    [1] 張哲, 牛智有. 茶葉加工過(guò)程中的物理特性變化規(guī)律[J]. 農(nóng)機(jī)化研究, 2012, 34(8): 116-119. Zhang Z, Niu Z Y. The changes of physical parameters in tea processing [J]. Journal of Agricultural Mechanization Research, 2012, 34(8): 116-119.

    [2] 權(quán)啟愛(ài). 茶葉揉捻機(jī)的發(fā)明及我國(guó)茶葉揉捻機(jī)研制與發(fā)展[J]. 中國(guó)茶葉, 2017, 39(11): 7-9. Quan Q A. Invention of tea twisting machine and development of tea twisting machine in China [J]. Journal of China Tea, 2017, 39(11): 7-9.

    [3] 陳世輝. 茶葉揉捻機(jī)外揉盤凹傾角對(duì)揉捻特性的影響研究[J]. 農(nóng)機(jī)化研究, 2015, 37(9): 41-44. Chen S H. Reseach on the effect of the relationship between inclination angle and rolling characteristic [J]. Journal of Agricultural Mechanization Research, 2015, 37(9): 41-44.

    [4] 岳鵬翔, 張桂銀. 茶葉揉捻機(jī)棱骨安裝的優(yōu)化參數(shù)[J]. 茶葉科學(xué), 1995, 25(1): 43-48. Yue P X, Zhang G Y. Optimized designing of the batten mounting parameters for tea rolling machines [J]. Journal of Tea Science, 1995, 25(1): 43-48.

    [5] 周昊, 楊旭東, 曹菡婷, 等. 全自動(dòng)茶葉揉捻機(jī)棱骨參數(shù)的研究[J]. 現(xiàn)代機(jī)械, 2017(2): 19-23. Zhou H, Yang X D, Cao H T, et al. Study on batten parameters of automatic tea twisting machine [J]. Modern Machinery, 2017(2): 19-23.

    [6] 瞿廷怡, 王濤, 談鑫, 等. 茶葉揉捻機(jī)支臂的優(yōu)化設(shè)計(jì)[J]. 機(jī)械工程師, 2013(8): 105-106. Qu T Y, Wang T, Tan X, et al. Optimization design of the supporting beam in tea twisting machine [J]. Mechanical Engineer, 2013(8): 105-106.

    [7] 劉凡一, 張艦, 李博, 等. 基于堆積試驗(yàn)的小麥離散元參數(shù)分析及標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2016, 32(12): 247-253. Liu F Y, Zhang J, Li B, 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, 2016, 32(12): 247-253.

    [8] 劉月琴,趙滿全, 劉飛, 等. 基于離散元的氣吸式排種器工作參數(shù)仿真優(yōu)化[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2016, 47(7): 65-72. Liu Y Q, Zhao M Q, Liu F, et al. Simulation and optimization of working parameters of air suction metering device based on discrete element [J]. Transactions of The Chinese Society of Agricultural Machinery, 2016, 47(7): 65-72.

    [9] 張強(qiáng), 張旭, 孫紹安. 基于EDEM的家用榨油機(jī)壓榨腔工作性能仿真研究[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2018, 34(24): 283-291. Zhang Q, Zhang X, Sun S A. Simulation study on working performance of pressing cavity of household oil press based on EDEM [J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(24): 283-291.

    [10] 王芳, 呂冰, 王洪明, 等. 氣吸式谷子排種裝置吸種孔的結(jié)構(gòu)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2017, 33(8): 30-36. Wang F, Lv B, Wang H M, et al. Structural design and test of seed-suction hole of air-sucking seed-metering device for millet [J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(8): 30-36.

    [11] 張祖明, 張銘耀, 江子穎, 等. 自動(dòng)化茶葉揉捻機(jī)的設(shè)計(jì)與實(shí)現(xiàn)[J]. 中國(guó)高新科技, 2018(1): 68-70. Zhang Z M, Zhang M Y, Jiang Z Y, et al. Design and implementation of automatic tea rolling machine [J]. Journal of China High Tech, 2018(1): 68-70.

    [12] 趙章風(fēng), 任超, 鄭勁松, 等. 球形茶炒制設(shè)備的做功效率與茶葉成形關(guān)系研究[J]. 茶葉科學(xué), 2018, 38(5): 527-536. Zhao Z F, Ren C, Zheng J S, et al. Study on the relationship between work efficiency of spherical tea frying equipment and tea formation [J]. Journal of Tea Science, 2018, 38(5): 527-536.

    [13] 林宏清. 綠茶揉捻機(jī)主要性能結(jié)構(gòu)設(shè)計(jì)參數(shù)的初步分析[J]. 茶葉科學(xué), 1965(3): 66-72. Lin H Q. Primary analysis of main performance and structure design parameters of green tea twisting machine [J]. Journal of Tea Science, 1965(3): 66-72.

    [14] 頓國(guó)強(qiáng), 陳海濤, 馮夷寧, 等. 基于EDEM軟件的肥料調(diào)配裝置關(guān)鍵部件參數(shù)優(yōu)化與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2016, 32(7): 36-42. Dun G Q, Chen H T, Feng Y N, et al. Parameter optimization and test of key parts of fertilizer allocation device based on EDEM software [J]. Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(7): 36-42.

    [15] Shi Y Y, Xin S, Wang X C, et al. Numerical simulation and field tests of minimum-tillage planter with straw smashing and strip laying based on EDEM software [J]. Computers and Electronics in Agriculture, 2019, 166(11): 1-9.

    [16] 虞文俊, 吳瑞梅, 李紅, 等. 基于Fluent-EDEM耦合的茶葉紅外殺青機(jī)滾筒內(nèi)流場(chǎng)數(shù)值模擬[J]. 食品與機(jī)械, 2019, 35(8): 104-109. Yu W J, Wu R M, Li H, et al. The numerical simulation on temperature field inside the radiation de-enzyme machine based on Fluent-EDEM coupling [J]. Food & Machinery, 2019, 35(8): 104-109.

    [17] 李兵, 李為寧, 柏宣丙. 基于EDEM的茶鮮葉分級(jí)機(jī)的篩分率的研究[J]. 茶葉科學(xué), 2019, 39(4): 484-494. Li B, Li W N, Bai X B. Research on screening rate of fresh tea leaves classifier based on EDEM [J]. Journal of Tea Science, 2019, 39(4): 484-494.

    [18] 朱志楠, 趙章風(fēng), 鐘江, 等. 基于多相流耦合過(guò)程數(shù)值模擬的茶鮮葉離心式連續(xù)脫水設(shè)備參數(shù)模擬與優(yōu)化[J]. 茶葉科學(xué), 2017, 37(3): 280-289.Zhu Z N, Zhao Z F, Zhong J, et al. Study on parameter optimization of centrifugal and continuous dewatering of tea leaves based on numerical simulation of multiphase flow [J]. Journal of Tea Science, 2017, 37(3): 280-289.

    [19] 岳鵬翔. 茶葉揉捻機(jī)棱骨作用研究[J]. 茶葉機(jī)械雜志, 1994(1): 27-29.Yue P X. Study on the role of the ribs of tea twisting machine [J]. Journal of Tea Machine, 1994(1): 27-29.

    [20] 張哲, 牛智有. 揉捻過(guò)程中茶葉物理特性的變化規(guī)律[J]. 湖北農(nóng)業(yè)科學(xué), 2012, 51(13): 2767-2770. Zhang Z, Niu Z Y. The variation law of physical properties of tea in the rolling process [J]. Hubei Agricultural Sciences, 2012, 51(13): 2767-2770.

    [21] 彭飛, 方芳, 王紅英, 等. 基于CFD-DEM的飼料調(diào)質(zhì)器物料運(yùn)動(dòng)模擬與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2018, 49(12): 355-363. Peng F, Fang F, Wang H Y, et al. Numerical simulation and experiment on raw material in feed conditioner based on coupled CFD-DEM [J]. Transactions of the Chinese Society of Agricultural Machinery, 2018, 49(12): 355-363.

    [22] 朱德泉, 李蘭蘭, 文世昌, 等. 滑片型孔輪式水稻精量排種器排種性能數(shù)值模擬與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2018, 34(21): 17-26. Zhu D Q, Li L L, Wen S C, et al. Numerical simulation and experiment on seeding performance of slide hole-wheel precision seed-metering device for rice [J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(21): 17-26.

    [23] 趙淑紅, 陳君執(zhí), 王加一, 等. 精量播種機(jī)V型凹槽撥輪式導(dǎo)種部件設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2018, 49(6): 146-158. Zhao S H, Chen J Z, Wang J Y, et al. Design and experiment on V-groove dialing round type guiding-seed device [J]. Transactions of the Chinese Society of Agricultural Machinery, 2018, 49(6): 146-158.

    [24] 王金武, 唐漢, 王金峰, 等. 懸掛式水田單側(cè)修筑埂機(jī)數(shù)值模擬分析與性能優(yōu)化[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2017, 48(8): 72-80. Wang J W, Tang H, Wang J F, et al. Numerical analysis and performance optimization experiment on hanging unilateral ridger for paddy field [J]. Transactions of the Chinese Society of Agricultural Machinery, 2017, 48(8): 72-80.

    Parameter Optimization and Experimental Study of Tea Twisting Machine Based on EDEM

    LI Bing1,2, LI Weining1, BAI Xuanbing1, HUANG Jianhong3

    1. School of Engineering, Anhui Agricultural University, Hefei 230036, China; 2. State Key Laboratory of Tea Plant Biology and Utilization, Anhui Agricultural University, Hefei 230036, China;3. Zhejiang Sunyoung Machinery Co., Ltd, Quzhou 324000, China

    To solve the problems caused by unstable quality in the process of tea twisting, the 6CR-40 tea twisting machine was taken as the research object in this study. The three-dimensional modeling of tea twisting machine was established based on Solidworks. The numerical simulation of the twisting process of tea twisting machine was carried out by the discrete element simulation software EDEM. The influences of various experimental factors on the performance indexes of tea twisting machine were obtained. The quadratic orthogonal rotation test was carried out and the Design-Expert was used to optimize the solution and obtain the best combination of structural parameters of twisting quality. The results showed that: when the twisting barrel rotational speed, the prismatic height, twisting plate inclination angles, forming rate of tea and breaking rate of tea were 42?r·min-1, 10?mm, 3.8°, 88.55% and 1.83% respectively, the tea twisting machine had a good quality of twisting. The results of verification test and simulation optimization were basically consistent.

    tea twisting machine, EDEM, simulation analysis, quadratic orthogonal rotation test, response surface method

    S571.1;TS272.7

    A

    1000-369X(2020)03-375-11

    2019-10-22

    2019-12-09

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2018YFD0700504、2017YFD0400801)、安徽省教育廳自然重點(diǎn)項(xiàng)目(KJ2017A133)

    李兵,男,副教授,博士,主要從事茶園及茶葉機(jī)械方面的研究,libing@ahau.edu.cn

    投稿平臺(tái):http://cykk.cbpt.cnki.net

    猜你喜歡
    條率傾角茶葉
    《茶葉通訊》簡(jiǎn)介
    茶葉通訊(2022年2期)2022-11-15 08:53:56
    地球軸傾角的改斜歸正
    激光傾角儀在CT引導(dǎo)下經(jīng)皮肺穿刺活檢中的應(yīng)用
    車輪外傾角和前束角匹配研究
    北京汽車(2021年3期)2021-07-17 08:25:42
    基于線性功效系數(shù)法的紫薯面條高效干燥技術(shù)研究
    藏族對(duì)茶葉情有獨(dú)鐘
    創(chuàng)造(2020年5期)2020-09-10 09:19:22
    復(fù)合增筋劑對(duì)面條品質(zhì)的影響參數(shù)優(yōu)化研究
    黑米黑豆面條擠壓工藝參數(shù)的優(yōu)化
    香噴噴的茶葉
    PGA對(duì)紫薯面條斷條率及烹煮品質(zhì)的影響
    国产午夜福利久久久久久| 成人午夜高清在线视频| 国内少妇人妻偷人精品xxx网站| 国产熟女xx| 国产在视频线在精品| 国产v大片淫在线免费观看| 久9热在线精品视频| 久久亚洲真实| 国内毛片毛片毛片毛片毛片| 免费高清视频大片| av国产免费在线观看| 亚洲专区中文字幕在线| 国产激情欧美一区二区| 亚洲在线观看片| svipshipincom国产片| 国产成人欧美在线观看| 国产中年淑女户外野战色| 久久久国产成人免费| 成人无遮挡网站| 成人午夜高清在线视频| 性色av乱码一区二区三区2| 国产伦人伦偷精品视频| 欧美黄色片欧美黄色片| 最新中文字幕久久久久| 国产精品1区2区在线观看.| 国产色爽女视频免费观看| 成熟少妇高潮喷水视频| 久久6这里有精品| 男女之事视频高清在线观看| 国内毛片毛片毛片毛片毛片| 亚洲av日韩精品久久久久久密| 90打野战视频偷拍视频| 亚洲激情在线av| 99久久精品热视频| av国产免费在线观看| 日韩 欧美 亚洲 中文字幕| 久久久久久人人人人人| 亚洲人成网站在线播放欧美日韩| 99久久九九国产精品国产免费| 国产伦精品一区二区三区四那| 美女被艹到高潮喷水动态| 亚洲av成人精品一区久久| 亚洲欧美日韩高清在线视频| 国产成人av激情在线播放| 两个人看的免费小视频| 亚洲不卡免费看| 夜夜夜夜夜久久久久| 黄色片一级片一级黄色片| 亚洲性夜色夜夜综合| 欧美一区二区亚洲| 一个人免费在线观看电影| 在线观看美女被高潮喷水网站 | 性色avwww在线观看| 欧美成人一区二区免费高清观看| 久久九九热精品免费| 一本一本综合久久| 久久九九热精品免费| 最近在线观看免费完整版| 精品熟女少妇八av免费久了| xxxwww97欧美| xxxwww97欧美| 国产精品久久久久久久久免 | 90打野战视频偷拍视频| 久久久久久久精品吃奶| 欧美日韩乱码在线| 午夜福利在线观看吧| 麻豆成人午夜福利视频| 成年人黄色毛片网站| 久久性视频一级片| 久9热在线精品视频| 成年女人毛片免费观看观看9| 欧洲精品卡2卡3卡4卡5卡区| 一本精品99久久精品77| 真人做人爱边吃奶动态| 无限看片的www在线观看| 国产精品99久久99久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 久久久国产精品麻豆| 日韩中文字幕欧美一区二区| 又紧又爽又黄一区二区| 日本免费一区二区三区高清不卡| 欧美又色又爽又黄视频| 精品久久久久久,| 国产蜜桃级精品一区二区三区| 欧美黑人欧美精品刺激| 亚洲熟妇中文字幕五十中出| 香蕉av资源在线| 国产蜜桃级精品一区二区三区| 久久99热这里只有精品18| 日韩欧美一区二区三区在线观看| 国产精品一及| 黄色视频,在线免费观看| 91在线精品国自产拍蜜月 | netflix在线观看网站| 少妇高潮的动态图| 久久精品91蜜桃| 国产av在哪里看| 国产伦精品一区二区三区视频9 | 国产麻豆成人av免费视频| 首页视频小说图片口味搜索| 亚洲天堂国产精品一区在线| 免费一级毛片在线播放高清视频| 欧美绝顶高潮抽搐喷水| 色老头精品视频在线观看| 久久精品91无色码中文字幕| av视频在线观看入口| 真人一进一出gif抽搐免费| 成年免费大片在线观看| 精品人妻一区二区三区麻豆 | 两个人视频免费观看高清| 男女视频在线观看网站免费| 日韩亚洲欧美综合| 精品日产1卡2卡| 看免费av毛片| 波多野结衣巨乳人妻| 亚洲av五月六月丁香网| 波多野结衣巨乳人妻| 亚洲性夜色夜夜综合| 中出人妻视频一区二区| 真人做人爱边吃奶动态| av在线蜜桃| 久久久久久人人人人人| av黄色大香蕉| 国内久久婷婷六月综合欲色啪| 日本三级黄在线观看| 99久久精品一区二区三区| 天天一区二区日本电影三级| 日本黄大片高清| 亚洲精品日韩av片在线观看 | 亚洲av第一区精品v没综合| 黄片大片在线免费观看| 国产单亲对白刺激| 真实男女啪啪啪动态图| 免费av观看视频| 午夜两性在线视频| 日本与韩国留学比较| 美女黄网站色视频| 天天添夜夜摸| 亚洲精华国产精华精| 可以在线观看的亚洲视频| 亚洲国产精品久久男人天堂| 天堂av国产一区二区熟女人妻| 国产亚洲精品一区二区www| 熟女电影av网| 久久久久亚洲av毛片大全| 99久久精品热视频| 精品久久久久久久毛片微露脸| 成人无遮挡网站| 97超级碰碰碰精品色视频在线观看| 日本三级黄在线观看| 成人精品一区二区免费| 国产精品1区2区在线观看.| 久久精品人妻少妇| 久久婷婷人人爽人人干人人爱| 国产亚洲欧美在线一区二区| 久久天躁狠狠躁夜夜2o2o| 非洲黑人性xxxx精品又粗又长| 色尼玛亚洲综合影院| 亚洲国产欧美网| 精品不卡国产一区二区三区| 老熟妇仑乱视频hdxx| 久久精品夜夜夜夜夜久久蜜豆| www.999成人在线观看| 亚洲av第一区精品v没综合| 99在线视频只有这里精品首页| 狂野欧美激情性xxxx| 成人欧美大片| 亚洲色图av天堂| 日本五十路高清| 精品乱码久久久久久99久播| 欧美日韩中文字幕国产精品一区二区三区| 老鸭窝网址在线观看| 国产欧美日韩一区二区精品| 特大巨黑吊av在线直播| 国产色婷婷99| 国内毛片毛片毛片毛片毛片| 亚洲精品成人久久久久久| 国产精品一区二区三区四区久久| 久久久久久久精品吃奶| 51国产日韩欧美| 午夜久久久久精精品| 内射极品少妇av片p| 欧美一区二区精品小视频在线| 国产精品av视频在线免费观看| 黄色视频,在线免费观看| 女警被强在线播放| 亚洲avbb在线观看| 成人精品一区二区免费| 亚洲av五月六月丁香网| 欧美av亚洲av综合av国产av| 国产一区二区在线av高清观看| 欧美成人a在线观看| 国产精品美女特级片免费视频播放器| 久久草成人影院| 在线播放无遮挡| 人妻夜夜爽99麻豆av| 午夜免费观看网址| 岛国在线免费视频观看| 美女 人体艺术 gogo| 国产又黄又爽又无遮挡在线| 免费看光身美女| 国产男靠女视频免费网站| 亚洲,欧美精品.| 欧美成人a在线观看| 亚洲欧美精品综合久久99| 午夜视频国产福利| 国产午夜福利久久久久久| 草草在线视频免费看| 乱人视频在线观看| 麻豆成人午夜福利视频| 中文字幕人妻熟人妻熟丝袜美 | 波多野结衣高清作品| 性色avwww在线观看| 成人亚洲精品av一区二区| 亚洲无线观看免费| 国产乱人视频| 国产精品免费一区二区三区在线| 观看美女的网站| 熟女人妻精品中文字幕| 亚洲最大成人手机在线| 亚洲男人的天堂狠狠| 中文资源天堂在线| 一进一出抽搐gif免费好疼| 91久久精品电影网| 国产精品一及| 免费高清视频大片| 性色av乱码一区二区三区2| 国产视频内射| 精品一区二区三区av网在线观看| 一级黄色大片毛片| 91av网一区二区| av在线天堂中文字幕| 亚洲人与动物交配视频| 国内精品久久久久精免费| 亚洲成人免费电影在线观看| 午夜影院日韩av| 久久久久久大精品| 别揉我奶头~嗯~啊~动态视频| 久久久久久久亚洲中文字幕 | 岛国在线观看网站| 人人妻,人人澡人人爽秒播| 亚洲内射少妇av| 一边摸一边抽搐一进一小说| 少妇的丰满在线观看| 91九色精品人成在线观看| 欧美+亚洲+日韩+国产| 丰满乱子伦码专区| 嫁个100分男人电影在线观看| 国内久久婷婷六月综合欲色啪| 手机成人av网站| 国产一区在线观看成人免费| 精品久久久久久久毛片微露脸| 久99久视频精品免费| 99热精品在线国产| 成年女人看的毛片在线观看| 88av欧美| 男人的好看免费观看在线视频| 国产久久久一区二区三区| 婷婷亚洲欧美| 草草在线视频免费看| 男女视频在线观看网站免费| 亚洲成a人片在线一区二区| 韩国av一区二区三区四区| 国产精品影院久久| 99国产综合亚洲精品| 别揉我奶头~嗯~啊~动态视频| 两个人看的免费小视频| 亚洲精品一卡2卡三卡4卡5卡| 日本 av在线| 亚洲精华国产精华精| 精品99又大又爽又粗少妇毛片 | 欧美性感艳星| 老司机午夜十八禁免费视频| 色综合婷婷激情| 婷婷精品国产亚洲av在线| 国产精品自产拍在线观看55亚洲| 免费在线观看成人毛片| 在线国产一区二区在线| 亚洲专区中文字幕在线| 国产一区二区亚洲精品在线观看| 国产麻豆成人av免费视频| 天堂av国产一区二区熟女人妻| av天堂在线播放| 免费看美女性在线毛片视频| 国产精品久久久久久久久免 | 亚洲第一欧美日韩一区二区三区| 无限看片的www在线观看| 国产成人系列免费观看| 日日干狠狠操夜夜爽| 色av中文字幕| 99久久九九国产精品国产免费| 99久久精品热视频| 一本久久中文字幕| 一级毛片女人18水好多| 久久国产精品影院| 国产爱豆传媒在线观看| avwww免费| 日韩欧美精品免费久久 | 国产精品一区二区免费欧美| 亚洲av五月六月丁香网| 91麻豆精品激情在线观看国产| 99riav亚洲国产免费| 欧美大码av| 国产一级毛片七仙女欲春2| 亚洲av日韩精品久久久久久密| 99在线视频只有这里精品首页| 在线播放国产精品三级| 变态另类成人亚洲欧美熟女| 美女免费视频网站| 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 伊人久久精品亚洲午夜| 热99在线观看视频| 欧美激情在线99| 国产精品 欧美亚洲| 欧美三级亚洲精品| 久久久久久国产a免费观看| 国产麻豆成人av免费视频| 9191精品国产免费久久| 亚洲真实伦在线观看| 在线观看美女被高潮喷水网站 | 舔av片在线| 国产单亲对白刺激| 国产精品一区二区三区四区免费观看 | 夜夜爽天天搞| 国产 一区 欧美 日韩| 亚洲中文字幕日韩| 亚洲av成人av| 国产亚洲精品久久久com| 国内精品久久久久久久电影| 中文字幕高清在线视频| 最好的美女福利视频网| 噜噜噜噜噜久久久久久91| 国产单亲对白刺激| 老司机在亚洲福利影院| 色综合亚洲欧美另类图片| 禁无遮挡网站| 窝窝影院91人妻| av中文乱码字幕在线| 69av精品久久久久久| 国产成人a区在线观看| 欧美一区二区精品小视频在线| 亚洲性夜色夜夜综合| 99久国产av精品| 一个人看的www免费观看视频| 97超视频在线观看视频| 日韩有码中文字幕| 特级一级黄色大片| 国产精品香港三级国产av潘金莲| 日韩av在线大香蕉| 黄色丝袜av网址大全| 国产一级毛片七仙女欲春2| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成a人片在线一区二区| 成人国产综合亚洲| 美女高潮喷水抽搐中文字幕| 久9热在线精品视频| 韩国av一区二区三区四区| 色噜噜av男人的天堂激情| 禁无遮挡网站| netflix在线观看网站| av欧美777| 免费观看的影片在线观看| 国产欧美日韩精品一区二区| 亚洲国产精品sss在线观看| 亚洲五月婷婷丁香| 免费av观看视频| 一区福利在线观看| 真人一进一出gif抽搐免费| 日韩欧美在线乱码| 精品久久久久久,| 一区二区三区免费毛片| 窝窝影院91人妻| 精品久久久久久,| 高清毛片免费观看视频网站| 看片在线看免费视频| 18禁黄网站禁片免费观看直播| 在线观看午夜福利视频| netflix在线观看网站| 在线观看一区二区三区| 男人和女人高潮做爰伦理| 免费一级毛片在线播放高清视频| 2021天堂中文幕一二区在线观| 亚洲av第一区精品v没综合| 日韩欧美在线二视频| 大型黄色视频在线免费观看| 一夜夜www| xxx96com| av天堂在线播放| 亚洲精品影视一区二区三区av| 成人欧美大片| 男人舔女人下体高潮全视频| 午夜精品久久久久久毛片777| 在线国产一区二区在线| 久久亚洲精品不卡| 亚洲国产色片| 欧美激情久久久久久爽电影| 国产成年人精品一区二区| 99久久精品国产亚洲精品| 国产老妇女一区| 嫩草影院精品99| 最近最新中文字幕大全电影3| 国产三级中文精品| 日本 欧美在线| 国产亚洲精品综合一区在线观看| 男女做爰动态图高潮gif福利片| 国产淫片久久久久久久久 | 国产亚洲av嫩草精品影院| 国产激情偷乱视频一区二区| 亚洲片人在线观看| av女优亚洲男人天堂| 伊人久久精品亚洲午夜| 日本精品一区二区三区蜜桃| 两人在一起打扑克的视频| 三级男女做爰猛烈吃奶摸视频| 极品教师在线免费播放| 国产av一区在线观看免费| 亚洲精品亚洲一区二区| 色精品久久人妻99蜜桃| 看黄色毛片网站| 九九久久精品国产亚洲av麻豆| 99精品久久久久人妻精品| 国产探花极品一区二区| 亚洲18禁久久av| 成熟少妇高潮喷水视频| 国产精品综合久久久久久久免费| 欧美在线一区亚洲| 国产伦精品一区二区三区四那| 国内揄拍国产精品人妻在线| 国产美女午夜福利| 日本成人三级电影网站| 国产av在哪里看| 免费搜索国产男女视频| 99久久99久久久精品蜜桃| 国产欧美日韩精品亚洲av| 法律面前人人平等表现在哪些方面| 在线观看av片永久免费下载| 精品久久久久久久末码| 亚洲国产欧美网| 91麻豆精品激情在线观看国产| 老司机午夜福利在线观看视频| 麻豆久久精品国产亚洲av| 欧美黑人巨大hd| 成人精品一区二区免费| 免费av毛片视频| www国产在线视频色| 999久久久精品免费观看国产| 欧美性猛交黑人性爽| 欧美国产日韩亚洲一区| 午夜福利在线在线| 黄片大片在线免费观看| av片东京热男人的天堂| 国产视频内射| 欧美三级亚洲精品| 搡老妇女老女人老熟妇| 99热这里只有是精品50| 3wmmmm亚洲av在线观看| 在线免费观看不下载黄p国产 | 男人的好看免费观看在线视频| 免费看光身美女| 九九热线精品视视频播放| 中文在线观看免费www的网站| 丁香欧美五月| 一进一出抽搐动态| 欧美日韩综合久久久久久 | 国产aⅴ精品一区二区三区波| 国产成人系列免费观看| 国产精品美女特级片免费视频播放器| 国产老妇女一区| 亚洲欧美日韩卡通动漫| 亚洲国产欧美网| 男女床上黄色一级片免费看| 亚洲精品乱码久久久v下载方式 | 99久久精品国产亚洲精品| 人妻夜夜爽99麻豆av| 高清毛片免费观看视频网站| 757午夜福利合集在线观看| 乱人视频在线观看| 99精品在免费线老司机午夜| 日本精品一区二区三区蜜桃| 法律面前人人平等表现在哪些方面| 日韩成人在线观看一区二区三区| 成人性生交大片免费视频hd| 精品不卡国产一区二区三区| 91久久精品国产一区二区成人 | 亚洲精品日韩av片在线观看 | 91九色精品人成在线观看| 女生性感内裤真人,穿戴方法视频| 中文字幕久久专区| 久久久久久久久久黄片| 欧洲精品卡2卡3卡4卡5卡区| 啦啦啦免费观看视频1| 美女高潮的动态| 亚洲精品一区av在线观看| 国产成人av教育| 国内少妇人妻偷人精品xxx网站| 午夜老司机福利剧场| 少妇高潮的动态图| 香蕉丝袜av| netflix在线观看网站| 欧美一区二区国产精品久久精品| 人妻丰满熟妇av一区二区三区| 在线免费观看的www视频| 熟女少妇亚洲综合色aaa.| 久久国产乱子伦精品免费另类| 中文字幕精品亚洲无线码一区| 国产高清视频在线播放一区| 1024手机看黄色片| 动漫黄色视频在线观看| 免费看十八禁软件| 网址你懂的国产日韩在线| 母亲3免费完整高清在线观看| 免费看光身美女| 精品免费久久久久久久清纯| 在线播放无遮挡| 综合色av麻豆| 美女被艹到高潮喷水动态| 亚洲国产欧洲综合997久久,| 久久人妻av系列| 色吧在线观看| 国产综合懂色| 黄片小视频在线播放| av黄色大香蕉| 国产精品,欧美在线| 午夜久久久久精精品| 高清日韩中文字幕在线| 精品人妻一区二区三区麻豆 | 少妇人妻一区二区三区视频| 成人鲁丝片一二三区免费| 亚洲欧美日韩无卡精品| 最新在线观看一区二区三区| 国产av在哪里看| 99久久精品国产亚洲精品| 午夜福利在线在线| aaaaa片日本免费| 亚洲精品在线观看二区| 亚洲人成伊人成综合网2020| 久久欧美精品欧美久久欧美| 国产麻豆成人av免费视频| 精品电影一区二区在线| 99精品久久久久人妻精品| 成年女人看的毛片在线观看| 亚洲成av人片免费观看| 亚洲精品在线美女| 久9热在线精品视频| ponron亚洲| 亚洲成av人片免费观看| 欧美日韩黄片免| 日本一二三区视频观看| 动漫黄色视频在线观看| 久久午夜亚洲精品久久| 国内精品美女久久久久久| 欧美成狂野欧美在线观看| 色哟哟哟哟哟哟| 精品熟女少妇八av免费久了| 亚洲 国产 在线| 99热6这里只有精品| 人人妻,人人澡人人爽秒播| av片东京热男人的天堂| 美女大奶头视频| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| www.熟女人妻精品国产| 午夜免费观看网址| 男女那种视频在线观看| 天堂动漫精品| 国产伦精品一区二区三区四那| 露出奶头的视频| 无限看片的www在线观看| 女人十人毛片免费观看3o分钟| 黄色日韩在线| 日韩高清综合在线| 日本黄色视频三级网站网址| 午夜精品在线福利| 日本免费一区二区三区高清不卡| 欧美色视频一区免费| 又粗又爽又猛毛片免费看| 亚洲成a人片在线一区二区| 亚洲欧美日韩东京热| 欧美成人性av电影在线观看| 国产精品99久久99久久久不卡| www.www免费av| 国产一区二区三区视频了| 久久性视频一级片| 日本免费一区二区三区高清不卡| 一级黄片播放器| 亚洲欧美日韩无卡精品| 亚洲久久久久久中文字幕| 精品国产美女av久久久久小说| 淫妇啪啪啪对白视频| 国产三级黄色录像| 国产亚洲精品一区二区www| 亚洲av熟女| 欧美区成人在线视频| 男人和女人高潮做爰伦理| 中文字幕av成人在线电影| 欧美色欧美亚洲另类二区| 久久精品国产自在天天线| 亚洲男人的天堂狠狠| 久久久久久久久大av| 中文字幕久久专区| 亚洲无线在线观看| 亚洲成av人片在线播放无| 狠狠狠狠99中文字幕| 三级男女做爰猛烈吃奶摸视频| 精品欧美国产一区二区三| 久久久久亚洲av毛片大全| 老司机午夜福利在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成人国产综合亚洲| 国产成人影院久久av| 一进一出抽搐gif免费好疼| 亚洲人成伊人成综合网2020| 美女免费视频网站|