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

    比例施肥泵吸肥活塞結(jié)構(gòu)優(yōu)化與試驗(yàn)

    2022-02-15 01:10:34易中懿
    關(guān)鍵詞:流槽肥液吸液

    湯 攀,任 妮,易中懿,李 紅

    比例施肥泵吸肥活塞結(jié)構(gòu)優(yōu)化與試驗(yàn)

    湯 攀1,任 妮2,易中懿2,李 紅1

    (1. 江蘇大學(xué)流體機(jī)械工程技術(shù)研究中心,鎮(zhèn)江 212013;2. 江蘇省農(nóng)業(yè)科學(xué)院信息中心,南京 210014)

    為了提高比例施肥泵的注肥精度,該研究分析了吸液活塞的工作原理,采用二次回歸正交組合試驗(yàn)對(duì)關(guān)鍵結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化,以吸液活塞下端直徑、泄流槽寬度以及泄流槽深度為變量,以注入流量為響應(yīng)指標(biāo),建立多元回歸模型,并通過試驗(yàn)進(jìn)行驗(yàn)證。結(jié)果表明:在不同壓差和設(shè)定肥液注入比例下,比例施肥泵的實(shí)際肥液注入比例均低于設(shè)定肥液注入比例。壓差在0.15 MPa以下時(shí),隨著設(shè)定肥液注入比例的升高,實(shí)際肥液注入比例與設(shè)定肥液注入比例的偏差減小,采用較高的設(shè)定肥液注入比例有利于提高注肥精度。吸液活塞下端直徑、泄流槽寬度、泄流槽深度對(duì)注入流量都有顯著影響(<0.01)。注入流量隨著吸液活塞下端直徑和泄流槽深度的增大而先升高后降低,隨泄流槽寬度的增大而增大。對(duì)注入流量的影響順序從大到小依次為泄流槽寬度、泄流槽深度、吸液活塞下端直徑。優(yōu)化后的吸液活塞下端直徑為16.6 mm、泄流槽寬度為5.5 mm和泄流槽深度為3.7 mm,工作壓差為0.05、0.10和0.15 MPa時(shí)的注肥精度分別提高了3.33、1.67和7.29個(gè)百分點(diǎn)。研究結(jié)果可為比例施肥泵的優(yōu)化設(shè)計(jì)及實(shí)際應(yīng)用提供理論支持。

    肥料;試驗(yàn);水肥一體化;比例施肥泵;參數(shù)優(yōu)化;回歸方程

    0 引 言

    長(zhǎng)期以來,肥料的大量與不合理施用對(duì)自然生態(tài)環(huán)境造成了嚴(yán)重破壞,如土壤肥力破壞、土壤污染、面源污染等,嚴(yán)重影響了農(nóng)業(yè)的可持續(xù)發(fā)展[1-2]。因此,減少化肥使用量、合理施肥、提高化肥利用率已成為農(nóng)業(yè)可持續(xù)發(fā)展和保障糧食安全的重要手段[3-5]。水肥一體化是精確施肥與精確灌溉相結(jié)合的產(chǎn)物,在灌溉技術(shù)中占有重要地位[6-8]。施肥設(shè)備是水肥一體化系統(tǒng)的關(guān)鍵,其性能的優(yōu)劣直接影響灌溉與施肥的質(zhì)量。比例施肥泵是一種先進(jìn)的水肥一體化施肥裝備,其通過水壓驅(qū)動(dòng)內(nèi)部吸液活塞的運(yùn)動(dòng)來向管網(wǎng)中定量添加肥液,與其他施肥設(shè)備相比,比例施肥泵的施肥精度高,且肥液注入比例可在一定范圍內(nèi)進(jìn)行調(diào)節(jié)[9-10]。

    法國(guó)Dosatron公司早在1974就開發(fā)出了比例施肥泵,此后經(jīng)過不斷發(fā)展,目前已開發(fā)出了系列產(chǎn)品,可以滿足不同應(yīng)用場(chǎng)合的需求。國(guó)內(nèi)在施肥泵方面的研究較晚,李百軍等[11]從原理及試驗(yàn)上初步研制過水動(dòng)比例施肥泵,但沒有形成產(chǎn)品。王建東等[12]制造出施肥泵樣品并進(jìn)行了試驗(yàn)測(cè)試,雖然試驗(yàn)結(jié)果表明該施肥泵可以基本滿足微灌施肥的需要,但與國(guó)外產(chǎn)品相比仍有可改進(jìn)空間。針對(duì)水力性能研究,韓啟彪等[13]對(duì)一些典型產(chǎn)品進(jìn)行了性能測(cè)試,根據(jù)試驗(yàn)結(jié)果給出了典型產(chǎn)品的工作壓差和流量關(guān)系以及工作壓差的合理控制范圍等較重要參數(shù)。楊大森等[14-15]對(duì)3種比例施肥泵的水力性能進(jìn)行了試驗(yàn),提出了比例施肥泵驅(qū)動(dòng)腔和吸肥腔的容積效率計(jì)算公式。Tang等[16-17]研究了比例施肥泵注入液體的黏度對(duì)進(jìn)口流量和注入流量的綜合影響,并建立了進(jìn)口流量計(jì)算模型和注入流量計(jì)算模型。目前已有研究主要側(cè)重于單個(gè)運(yùn)行及結(jié)構(gòu)參數(shù)對(duì)某一水力性能的簡(jiǎn)單外特性試驗(yàn)研究,缺乏總體運(yùn)行與結(jié)構(gòu)參數(shù)對(duì)綜合水力性能及其影響機(jī)理的深入研究。隨著計(jì)算流體動(dòng)力學(xué)(Computational Fluid Dynamics, CFD)的發(fā)展,數(shù)值模擬已成為目前研究流體機(jī)械內(nèi)部流動(dòng)規(guī)律的主要手段之一。湯攀等[18]在比例施肥泵運(yùn)行機(jī)理分析的基礎(chǔ)上,采用流固耦合動(dòng)網(wǎng)格技術(shù)進(jìn)行數(shù)值模擬,分析其內(nèi)部流動(dòng)以及活塞受力。此外,王睿等[19]研究了施肥泵施肥比例與肥水比對(duì)過濾器堵塞的影響,提出了網(wǎng)式和疊片式過濾器的適宜肥液濃度。吳錫凱[20]指出為了提高施肥精度,在運(yùn)行時(shí)不要使用較大壓差與較小施肥比例配合,當(dāng)肥液注入比例設(shè)定在2%~4%時(shí),實(shí)際所需的施肥時(shí)間與設(shè)定時(shí)間相比增加了近1/4,這表明實(shí)際注入流量要小于理論注入流量,即比例施肥泵的注肥精度仍有待提高。

    目前國(guó)內(nèi)比例施肥泵的研發(fā)和制造與國(guó)外相比尚有較大差距,依然缺乏自主的設(shè)計(jì)方法和優(yōu)秀產(chǎn)品。吸肥活塞作為比例施肥泵實(shí)現(xiàn)肥液抽吸的關(guān)鍵水力部件,其結(jié)構(gòu)參數(shù)對(duì)比例施肥泵的注肥精度具有重要影響?;诖耍疚耐ㄟ^旋轉(zhuǎn)優(yōu)化試驗(yàn)設(shè)計(jì)方法對(duì)吸肥活塞的關(guān)鍵結(jié)構(gòu)參數(shù)進(jìn)行研究,建立各參數(shù)與性能指標(biāo)的關(guān)系式,分析各因素的影響規(guī)律,以期為比例施肥泵的性能研究及關(guān)鍵水力部件優(yōu)化設(shè)計(jì)提供參考。

    1 吸液活塞工作原理

    比例施肥泵的主要結(jié)構(gòu)及工作原理已有詳細(xì)介紹[15-16],如圖1所示,吸液活塞在往復(fù)運(yùn)動(dòng)過程中首先需要將肥液從外部吸入到吸液活塞腔,然后將其排出到驅(qū)動(dòng)活塞腔與驅(qū)動(dòng)液體混合進(jìn)而排出泵體。吸液活塞在吸液過程中需保持密封狀態(tài),在排液過程中需起到將吸肥腔與驅(qū)動(dòng)腔連通的作用。

    如圖2所示為吸液活塞二維結(jié)構(gòu)圖,活塞頭兩側(cè)開有泄流槽,由于單獨(dú)活塞頭無法完成密封與連通功能的切換,因此在吸液活塞頭中間加入滑動(dòng)密封圈來控制泄流槽的連通與關(guān)閉,其主要結(jié)構(gòu)參數(shù)為:吸液活塞下端直徑15.5 mm、泄流槽寬度4.5 mm及泄流槽深度3 mm。

    注:d為吸液活塞下端直徑,mm;b1為泄流槽寬度,mm;b2為泄流槽深度,mm。

    如圖3a所示,當(dāng)吸液活塞向上運(yùn)動(dòng)時(shí),由于密封圈與吸液活塞腔壁面存在摩擦,使密封圈與吸液活塞頭的下端緊貼,從而阻斷泄流槽與吸液活塞腔的連通,使吸液活塞腔在活塞向上運(yùn)動(dòng)過程中形成負(fù)壓將外部液體吸進(jìn)活塞腔。如圖3b所示,當(dāng)吸液活塞向下運(yùn)動(dòng)時(shí),由于密封圈與吸液活塞腔壁面存在摩擦,使密封圈與吸液活塞頭的上端緊貼,從而通過泄流槽將吸液活塞腔與驅(qū)動(dòng)活塞腔進(jìn)行連通。吸液活塞向下運(yùn)動(dòng)的同時(shí)將活塞腔里面的液體擠壓,通過泄流槽排進(jìn)驅(qū)動(dòng)活塞腔與工作液體混合。

    注:箭頭表示液體流動(dòng)方向。

    2 試驗(yàn)方案

    為了研究比例施肥泵效率和實(shí)際注入流量與運(yùn)行參數(shù)之間的關(guān)系,搭建了比例施肥泵性能測(cè)試試驗(yàn)臺(tái),如圖4所示,通過調(diào)節(jié)比例施肥泵前后的閥門來控制其工作壓差。試驗(yàn)在江蘇大學(xué)流體機(jī)械工程技術(shù)研究中心的室內(nèi)噴灌試驗(yàn)室進(jìn)行,采用意大利Dosatron公司的D25RE2比例施肥泵,其進(jìn)出口直徑19 mm,最大承受壓力為0.60 MPa,設(shè)定肥液注入比例0.2%~2.0%。測(cè)試系統(tǒng)主要由水泵、測(cè)量裝置和控制裝置組成。系統(tǒng)管路為鋼管,主管路直徑50 mm,旁路管直徑25 mm。離心泵流量10 m3/h,揚(yáng)程70 m。主管道流量通過電磁流量計(jì)測(cè)量,精確度為±0.3%。比例施肥泵吸液速度由精度為0.1 g的精密電子天平測(cè)量。在比例施肥泵前后以及主管上安裝3個(gè)精度為±0.4%的壓力表,試驗(yàn)參照國(guó)家標(biāo)準(zhǔn)GB/T 19792-2012《農(nóng)業(yè)灌溉裝備水動(dòng)化肥—農(nóng)藥注入泵》[21]進(jìn)行。

    試驗(yàn)主要測(cè)試比例施肥泵在不同壓差和不同設(shè)定肥液注入比例下的注入流量和進(jìn)口流量。根據(jù)比例施肥泵的最大承受壓力,試驗(yàn)共設(shè)置8個(gè)水平的壓差,分別為0.05、0.10、0.15、0.20、0.25、0.30、0.35、0.40 MPa,試驗(yàn)過程中將比例施肥泵的進(jìn)口壓力保持在0.60 MPa,然后調(diào)節(jié)出口壓力,使比例施肥泵所設(shè)定的不同壓差下工作。設(shè)定肥液注入比例共設(shè)置13個(gè)水平,為了更加精確地測(cè)量肥液注入比例較小時(shí)的水力性能,當(dāng)設(shè)定肥液注入比例在0.8%以下時(shí),間隔為0.1%,當(dāng)設(shè)定肥液注入比例大于0.8%時(shí),間隔為0.2%。比例施肥泵的進(jìn)口流量直接通過電磁流量計(jì)獲得,由于注入流量相對(duì)進(jìn)口流量較小,為了提高測(cè)試精度,通過如下公式計(jì)算:

    1.電機(jī) 2.水泵 3.壓力表1 4.閥門1 5.電磁流量計(jì) 6.閥門2 7.壓力表2 8.壓力表3 9.閥門3 10.比例施肥泵 11.儲(chǔ)液罐 12.精密電子天平 13.水池1 14.水池2

    式中為注入流量,kg/h;1為儲(chǔ)液罐的初始質(zhì)量,kg;2為測(cè)試結(jié)束時(shí)儲(chǔ)液罐的質(zhì)量,kg;為測(cè)試時(shí)間,h。在本研究中,每次測(cè)試的時(shí)間為0.5 h。試驗(yàn)過程中,改變測(cè)試工作壓差和設(shè)定肥液注入比例時(shí),需待比例施肥泵運(yùn)行穩(wěn)定后開始測(cè)量,每個(gè)測(cè)試點(diǎn)3次重復(fù),取均值為最終試驗(yàn)結(jié)果。

    3 結(jié)果與分析

    表1為比例施肥泵在不同設(shè)定肥液注入比例和壓差下的注入流量。理論上,隨著壓差的增大,活塞往復(fù)運(yùn)動(dòng)頻率升高,吸液活塞吸液速度會(huì)增加,然而從表1可以看出,隨著壓差的增大,比例施肥泵在不同設(shè)定肥液注入比例下的注入流量都呈現(xiàn)先增大后趨于平穩(wěn)或減小的趨勢(shì),表明壓差的增大不能有效增加比例施肥泵的吸液速度,即活塞運(yùn)動(dòng)頻率過高會(huì)降低比例施肥泵的性能,其工作壓差應(yīng)控制在0.20 MPa以內(nèi)。表2為比例施肥泵在不同設(shè)定肥液注入比例和壓差下的進(jìn)口流量,從表2可以看出,比例施肥泵的進(jìn)口流量隨著壓差的增大而增大;在同一壓差下,隨著設(shè)定肥液注入比例的升高有較小幅度的降低,可認(rèn)為設(shè)定肥液注入比例對(duì)進(jìn)口流量的影響較小,這與吳錫凱[20]的研究相吻合。

    表1 不同設(shè)定肥液注入比例和壓差下的注入流量

    注:s為設(shè)定肥液注入比例,%;表中數(shù)據(jù)為均值±標(biāo)準(zhǔn)差,下同。

    Note:sis the set injection ratio, %; data in the table are mean ± standard deviation, the same as below.

    表2 不同設(shè)定肥液注入比例和壓差下的進(jìn)口流量

    實(shí)際肥液注入比例為實(shí)際注入流量與進(jìn)口流量的比值,即:

    式中a為實(shí)際肥液注入比例,%;as為實(shí)際注入流量,kg/h;in為進(jìn)口流量,kg/h。

    通過公式(2)計(jì)算可得如圖5所示的比例施肥泵在不同壓差下的實(shí)際肥液注入比例和設(shè)定肥液注入比例的對(duì)比,從圖中可以看出,隨著壓差的增大,實(shí)際肥液注入比例與設(shè)定肥液注入比例的偏差越大;壓差在0.15 MPa以下時(shí),隨著設(shè)定肥液注入比例的升高,實(shí)際肥液注入比例與設(shè)定肥液注入比例的偏差越小,表明采用較高的設(shè)定肥液注入比例有利于提高比例施肥泵的注入精度。Tang等[18]在研究注入液體黏度對(duì)進(jìn)口流量和注入流量影響時(shí)也同樣指出比例施肥泵的最大工作壓差不應(yīng)過高。在不同的壓差和設(shè)定肥液注入比例下,比例施肥泵的實(shí)際肥液注入比例都低于設(shè)定肥液注入比例,一方面,這是由于比例施肥泵存在一定的容積損失,無法按照設(shè)定的肥液注入比例進(jìn)行抽吸;另一方面,吸液活塞是實(shí)現(xiàn)肥液抽吸的關(guān)鍵水力部件,其關(guān)鍵結(jié)構(gòu)參數(shù)對(duì)比例施肥泵的注入流量具有重要的影響。

    圖5 不同壓差下的實(shí)際肥液注入比例和設(shè)定肥液注入比例

    4 參數(shù)優(yōu)化

    根據(jù)吸液活塞工作原理可知吸液活塞的結(jié)構(gòu)尺寸對(duì)液體流動(dòng)產(chǎn)生主要影響的有吸液活塞下端直徑、泄流槽寬度1以及泄流槽深度2。因此,選取這3個(gè)主要參數(shù)探索吸肥活塞流道結(jié)構(gòu)尺寸對(duì)比例施肥泵施肥精度的影響規(guī)律。

    比例施肥泵的注入精度是最能體現(xiàn)其水力性能優(yōu)越的重要指標(biāo),當(dāng)比例施肥泵的注入精度較低,即實(shí)際注入流量與理論注入流量差別較大時(shí),會(huì)造成施肥時(shí)間過長(zhǎng)或施藥濃度達(dá)不到目標(biāo)設(shè)定值,進(jìn)而影響施肥或施藥效果。根據(jù)圖5,比例施肥泵的實(shí)際肥液注入比例都低于設(shè)定肥液注入比例,即不同設(shè)定肥液注入比例下的實(shí)際注入流量都小于理論注入流量。比例施肥泵在同一設(shè)定肥液注入比例及壓差下工作時(shí),注入流量越大,則注入精度越高,故選擇注入流量作為本文的試驗(yàn)指標(biāo)。根據(jù)表1和圖5,設(shè)定肥液注入比例越低,實(shí)際注入流量與理論注入流量的偏差越大,故試驗(yàn)中將比例施肥泵的工作參數(shù)定為設(shè)定肥液注入比例0.2%,工作壓差0.05 MPa,研究吸液活塞流道結(jié)構(gòu)參數(shù)對(duì)注入精度的影響規(guī)律。

    4.1 二次回歸正交組合試驗(yàn)

    試驗(yàn)采用二次回歸正交組合設(shè)計(jì)[22-24],根據(jù)前期試驗(yàn)確定比例施肥泵能正常工作時(shí)各參數(shù)的取值范圍為:吸液活塞下端直徑(12~18 mm)、泄流槽寬度(4~6 mm)和泄流槽深度(2~4 mm)。由于因素?cái)?shù)=3,則二水平試驗(yàn)次數(shù)c=23=8,取零水平試驗(yàn)次數(shù)0=1,根據(jù)二次回歸正交試驗(yàn)設(shè)計(jì)中星號(hào)臂長(zhǎng)度的計(jì)算公式

    計(jì)算得到=1.215。因素的上下星號(hào)臂水平為因素的上下限,零水平為上下限的算術(shù)平均值,變化間距Δ的計(jì)算公式為

    式中x為上星號(hào)臂水平,0為零水平。

    因素的上下水平計(jì)算公式分別為

    根據(jù)試驗(yàn)設(shè)計(jì)得到試驗(yàn)因素水平編碼表見表3。

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

    根據(jù)因素水平編碼表的參數(shù)值及表4中的不同參數(shù)組合首先在Pro/E中完成三維造型,然后采用快速成型技術(shù)(Rapid Prototyping Manufacturing, RPM)加工,材料為ABS樹脂,活塞頭試件如圖6所示。試驗(yàn)測(cè)試結(jié)果見表4。

    圖6 吸液活塞頭試件

    將試驗(yàn)方案及結(jié)果導(dǎo)入Design-Expert 8.0.5.0,對(duì)試驗(yàn)數(shù)據(jù)處理和響應(yīng)面分析,采用回歸分析法建立三元二次回歸方程

    =0.084 25+0.094 167+0.263 751?0.022 752+

    0.008 3331+0.013 3332+0.052 512?

    0.005 8062?0.049 7512?0.069 7522(7)

    表4 二次回歸組合設(shè)計(jì)試驗(yàn)方案及試驗(yàn)結(jié)果

    注:為注入流量,kg·h-1。

    Note:is injection flowrate, kg·h-1.

    對(duì)試驗(yàn)結(jié)果進(jìn)行方差分析和顯著性檢驗(yàn),如表5所示。從表5中可以看出,注入流量模型的值小于0.01,注入流量回歸模型極顯著。模型決定系數(shù)2=0.957 4,表明模型能夠反映出95.74%的響應(yīng)值變化,失擬項(xiàng)不顯著,說明試驗(yàn)誤差小,試驗(yàn)各因素對(duì)注入流量影響的大小順序依次為泄流槽寬度、泄流槽深度、吸液活塞下端直徑;3個(gè)因素的交互作用對(duì)注入流量影響的大小順序依次為,泄流槽寬度×泄流槽深度、吸液活塞下端直徑×泄流槽深度、吸液活塞下端直徑×泄流槽寬度。顯著性分析結(jié)果表明所選試驗(yàn)因素都對(duì)注入流量有較為顯著的影響,且其各因素之間的相互交互作用同樣對(duì)注入流量具有顯著的影響,這表明本文所選擇的吸液活塞流道結(jié)構(gòu)參數(shù)都為關(guān)鍵水力結(jié)構(gòu)參數(shù)。

    表5 試驗(yàn)指標(biāo)方差分析表

    注:SS為離差平方和;df為自由度;MS為均方;*代表差異顯著(<0.05);**代表差異極顯著(<0.01)。

    Note: SS is sum of squares; df is degree of freedom; MS is mean squares; * is significant (<0.05); ** is very significant (<0.01).

    圖7為吸液活塞下端直徑與泄流槽寬度交互作用對(duì)注入流量的影響,從圖7可以看出,在同一泄流槽寬度下,注入流量隨著吸液活塞下端直徑的增大而先升高后降低,表明吸液活塞下端直徑在確定范圍內(nèi)存在一個(gè)最優(yōu)值。從吸液活塞結(jié)構(gòu)上分析,吸液活塞下端直徑過大或過小都不利于提高比例施肥泵的注入精度,這是由于當(dāng)吸液活塞下端直徑過大時(shí)會(huì)使吸液活塞與吸液活塞腔的間隙變小,增大被吸液體在排液過程中的阻力,從而使注入流量變?。幌喾?,當(dāng)吸液活塞下端直徑過小時(shí),吸液活塞與吸液活塞腔的間隙變大,造成吸液活塞下端對(duì)滑動(dòng)密封墊圈的支撐變小,導(dǎo)致吸液活塞在向上運(yùn)動(dòng)吸取液體時(shí)的密封性變差。吸液活塞下端直徑小于14 mm時(shí),注入流量同樣隨著泄流槽寬度的增大而先升高后降低,當(dāng)吸液活塞下端直徑在14~17 mm時(shí),注入流量隨著泄流槽寬度的增大而增大。

    圖7 吸液活塞下端直徑與泄流槽寬度交互作用對(duì)注入流量的影響

    圖8為吸液活塞下端直徑與泄流槽深度交互作用對(duì)注入流量的影響,從圖8可以看出,在設(shè)定范圍內(nèi),注入流量隨吸液活塞下端直徑與泄流槽深度的增大而先升高后降低。泄流槽深度過大或過小都會(huì)造成液體流經(jīng)吸液活塞的流道過大或過小,流道尺寸較大一方面會(huì)降低流體從吸液腔到混合器的流動(dòng)阻力,另一方面又會(huì)造成吸液活塞換向運(yùn)動(dòng)時(shí)(即由排液過程變?yōu)槲哼^程)殘留在泄流槽中的液體增多,從而影響注入流量。

    圖8 吸液活塞下端直徑與泄流槽深度交互作用對(duì)注入流量的影響

    圖9所示為泄流槽寬度與泄流槽深度交互作用對(duì)注入流量的影響,從圖中可以看出,此交互作用影響與吸液活塞下端直徑與泄流槽寬度交互作用對(duì)注入流量的影響類似。當(dāng)泄流槽深度小于2.5 mm時(shí),注入流量隨泄流槽寬度的增大基本保持不變,此時(shí)泄流槽寬度對(duì)注入流量的影響較??;當(dāng)泄流槽深度大于2.5 mm時(shí),注入流量隨泄流槽寬度的增大而升高。

    圖9 泄流槽寬度與泄流槽深度交互作用對(duì)注入流量的影響

    4.2 優(yōu)化求解及試驗(yàn)驗(yàn)證

    模型共有3個(gè)設(shè)計(jì)變量,是一個(gè)非線性數(shù)學(xué)規(guī)劃問題[25-26]。

    目標(biāo)函數(shù)如下

    根據(jù)式(8)目標(biāo)函數(shù)與約束條件進(jìn)行參數(shù)求解,得到如表6所示的吸液活塞流道結(jié)構(gòu)參數(shù)優(yōu)化結(jié)果及優(yōu)化前后的注入流量。從表6中可以看出,吸液活塞流道結(jié)構(gòu)參數(shù)優(yōu)化后的注入流量得到了較大幅度的提高,從優(yōu)化前的1.26 kg/h提高到了1.61 kg/h,提高幅度為27.78%。

    表6 優(yōu)化前后參數(shù)及性能對(duì)比

    對(duì)優(yōu)化結(jié)果進(jìn)行圓整,取吸液活塞下端直徑為16.6 mm、泄流槽寬度為5.5 mm和泄流槽深度為3.7 mm,并進(jìn)行驗(yàn)證試驗(yàn),測(cè)試優(yōu)化后比例施肥泵分別在壓差為0.05、0.10和0.15 MPa時(shí)不同設(shè)定肥液注入比例下的注入流量,并與原型比例施肥泵注入流量進(jìn)行對(duì)比,結(jié)果如表7所示。從表7中可以發(fā)現(xiàn),優(yōu)化后的注入流量在不同設(shè)定肥液注入比例和工作壓差下都較原型比例施肥泵有一定提高。特別是當(dāng)設(shè)定肥液注入比例較低時(shí),優(yōu)化后的實(shí)際注入流量更接近理論注入流量(進(jìn)口流量與設(shè)定肥液注入比例的乘積),同時(shí)實(shí)際注入流量與理論注入流量的偏差進(jìn)一步降低,在一定程度上提高了比例施肥泵的注入精度,從而有效提高了比例施肥泵的水力性能。當(dāng)工作壓差為0.15 MPa時(shí),優(yōu)化后的注入流量較優(yōu)化之前有一定的提高,且在較高設(shè)定肥液注入比例下的實(shí)際注入流量與理論注入流量的差距進(jìn)一步縮?。▔翰顬?.15 MPa和肥液注入比例為2.0%時(shí)的偏差從優(yōu)化之前的9.31%降低到優(yōu)化之后的4.24%),從而提高了比例施肥泵的正常工作壓差范圍,有利于比例施肥泵的實(shí)際應(yīng)用,因?yàn)楦呔鹊淖⒎时壤梢杂行Ы档退室惑w化系統(tǒng)肥液濃度的變化,從而提高施肥均勻性[27-29]。

    表7 不同壓差下比例施肥泵吸液活塞流道優(yōu)化前后的注入流量

    注:T.V為理論值;M.V為原型值;D.V1為原型值與理論值的偏差;O.V為優(yōu)化值;D.V2為優(yōu)化值與理論值的偏差。

    Note: T.V is the theoretical value; M.V is the prototype value; D.V1is the deviation between the prototype value and the theoretical value; O.V is the optimal value; D.V2is the deviation between the optimal value and the theoretical value.

    為進(jìn)一步分析吸肥活塞優(yōu)化前后比例施肥泵整體注肥精度的變化,將不同壓差下吸肥活塞流道結(jié)構(gòu)參數(shù)優(yōu)化前后比例施肥泵的實(shí)際肥液注入比例與對(duì)應(yīng)的設(shè)定肥液注入比例進(jìn)行對(duì)比,同時(shí)定義同一壓差下實(shí)際肥液注入比例與理論肥液注入比例比值的算術(shù)平均值為注肥精度,表示為

    式中為注肥精度,%;ti為比例施肥泵實(shí)際肥液注入比例,%;si為對(duì)應(yīng)的設(shè)定肥液注入比例,%;為同一壓差下的樣本數(shù)量;為同一壓差下的樣本編號(hào),1≤≤。通過計(jì)算得到如表8所示的對(duì)比結(jié)果,從表8中可以發(fā)現(xiàn),優(yōu)化后的實(shí)際肥液注入比例與設(shè)定肥液注入比例更接近,二者的偏差比優(yōu)化前的小。比例施肥泵在工作壓差為0.05、0.10和0.15 MPa時(shí)的注肥精度分別從82.74%、84.76%和65.29%提高到了86.07%、86.43%和72.58%,分別提高了3.33、1.67和7.29個(gè)百分點(diǎn)。

    表8 不同壓差下比例施肥泵吸液活塞流道優(yōu)化前后的注肥精度

    注:tm為原型比例施肥泵實(shí)際肥液注入比例,%;to為吸液活塞流道結(jié)構(gòu)參數(shù)優(yōu)化后比例施肥泵實(shí)際肥液注入比例,%;為注肥精度,%。

    Note:tmis the actual fertilizer injection rate of the prototype proportional fertilization pump, %;tois the actual fertilizer injection rate of the proportional fertilization pump after the structural parameters of the suction piston were optimized, %;is the fertilizer injection accuracy, %.

    5 結(jié) 論

    本研究通過二次回歸正交組合試驗(yàn)設(shè)計(jì)對(duì)比例施肥泵吸液活塞關(guān)鍵結(jié)構(gòu)參數(shù)進(jìn)行了優(yōu)化研究,得到的主要結(jié)論如下:

    1)在不同壓差和設(shè)定肥液注入比例下,比例施肥泵的實(shí)際肥液注入比例都低于設(shè)定肥液注入比例。壓差在0.15 MPa以下時(shí),隨著設(shè)定肥液注入比例的升高,實(shí)際肥液注入比例與設(shè)定肥液注入比例的偏差越小,采用較高的設(shè)定肥液注入比例有利于提高注肥精度。

    2)吸液活塞下端直徑、泄流槽寬度以及泄流槽深度都對(duì)注入流量有較為顯著的影響(<0.01);各因素在選取的水平范圍內(nèi)對(duì)注入流量的影響大小順序依次為泄流槽寬度、泄流槽深度、吸液活塞下端直徑。

    3)在同一泄流槽寬度下,注入流量隨著吸液活塞下端直徑的增大而先升高后降低。當(dāng)吸液活塞下端直徑小于14 mm時(shí),注入流量隨著泄流槽寬度的增大而先升高后降低,當(dāng)吸液活塞下端直徑在14~17 mm時(shí),注入流量隨著泄流槽寬度的增大而增大。注入流量隨著吸液活塞下端直徑與泄流槽深度的增大而先升高后降低。當(dāng)泄流槽深度小于2.5 mm時(shí),注入流量隨泄流槽寬度的增大基本保持不變;當(dāng)泄流槽深度大于2.5 mm時(shí),注入流量隨泄流槽寬度的增大而升高。

    4)優(yōu)化后的吸液活塞下端直徑為16.6 mm、泄流槽寬度為5.5 mm、泄流槽深度為3.7 mm。比例施肥泵在工作壓差為0.05、0.10和0.15 MPa時(shí)的注肥精度較優(yōu)化前分別提高了3.33、1.67和7.29個(gè)百分點(diǎn)。

    [1] 劉欽普,孫景榮,濮勵(lì)杰. 中國(guó)及歐美主要國(guó)家化肥施用強(qiáng)度與綜合效率比較研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(14):9-16.

    Liu Qinpu, Sun Jingrong, Pu Lijie. Comparative study on fertilization intensity and integrated efficiency in China and Euro-American major countries[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 9-16. (in Chinese with English abstract)

    [2] 王晶晶,卓越,張顥暉,等. 柱塞式注肥泵的設(shè)計(jì)試驗(yàn)及精準(zhǔn)施肥應(yīng)用[J]. 排灌機(jī)械工程學(xué)報(bào),2022,40(1):22-29.

    Wang Jingjing, Zhuo Yue, Zhang Haohui, et al. Design and test of piston injection pump and its application in precise fertigation[J]. Journal of Drainage and Irrigation Machinery Engineering, 2022, 40(1): 22-29. (in Chinese with English abstract)

    [3] Chiaregato C G, Souza C F, Faez R. The fertilizer release into water and soil as the biodegradation process in the sustainable material enhancing the fertilizer efficiency[J]. Environmental Technology & Innovation, 2021, 22: 101417.

    [4] Xu M, Liu M J, Liu F, et al. A safe, high fertilizer-efficiency and economical approach based on a low-volume spraying UAV loaded with chelated-zinc fertilizer to produce zinc-biofortified rice grains[J]. Journal of Cleaner Production, 2021, 323: 129188.

    [5] Zhang M, Yao Y L, Tian Y H, et al. Increasing yield and N use efficiency with organic fertilizer in Chinese intensive rice cropping systems[J]. Field Crops Research, 2018, 227: 102-109.

    [6] 劉俊萍,朱興業(yè),袁壽其,等. 中國(guó)農(nóng)業(yè)節(jié)水噴微灌裝備研究進(jìn)展及發(fā)展趨勢(shì)[J]. 排灌機(jī)械工程學(xué)報(bào),2022,40(1):87-96.

    Liu Junping, Zhu Xingye, Yuan Shouqi, et al. Research and development trend of agricultural water-saving sprinkler and micro-irrigation equipment in China[J]. Journal of Drainage and Irrigation Machinery Engineering, 2022, 40(1): 87-96. (in Chinese with English abstract)

    [7] 朱德蘭,阮漢鋮,吳普特,等. 水肥一體機(jī)肥液電導(dǎo)率遠(yuǎn)程模糊PID控制策略[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2022,53(1): 186-191.

    Zhu Delan, Ruan Hancheng, Wu Pute, et al. Strategy on remote fuzzy PID control for fertilizer liquid conductivity of water fertilizer integrated machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2022, 53(1): 186-191. (in Chinese with English abstract)

    [8] 王振華,陳瀟潔,呂德生,等. 水肥耦合對(duì)加氣滴灌加工番茄產(chǎn)量及品質(zhì)的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(19):66-75.

    Wang Zhenhua, Chen Xiaojie, Lyu Desheng, et al. Effects of water and fertilizer coupling on the yield and quality of processing tomato under aerated drip irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(19): 66-75. (in Chinese with English abstract)

    [9] 李紅,湯攀,陳超,等. 中國(guó)水肥一體化施肥設(shè)備研究現(xiàn)狀與發(fā)展趨勢(shì)[J]. 排灌機(jī)械工程學(xué)報(bào),2021,39(2):200-209.

    Li Hong, Tang Pan, Chen Chao, et al. Research status and development trend of fertilization equipment used in fertigation in China[J]. Journal of Drainage and Irrigation Machinery Engineering, 2021, 39(2): 200-209. (in Chinese with English abstract)

    [10] Tang P, Li H, Issaka Z, et al. Effect of manifold layout and fertilizer solution concentration on fertilization and flushing times and uniformity of drip irrigation systems[J]. Agricultural Water Management, 2018, 200: 71-79.

    [11] 李百軍,王曉寧. 水動(dòng)施肥裝置的設(shè)計(jì)與試驗(yàn)[J]. 江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版),2002,23(2):9-12.

    Li Baijun, Wang Xiaoning. Design and experiment on hydrodynamic fertilizer injection unit[J]. Journal of Jiangsu University (Natural Science Edition), 2002, 23(2): 9-12. (in Chinese with English abstract)

    [12] 王建東,龔時(shí)宏,徐茂云,等. 微灌用水動(dòng)活塞式施肥泵研制[J]. 農(nóng)業(yè)工程學(xué)報(bào),2006,22(6):100-103.

    Wang Jiandong, Gong Shihong, Xu Maoyun, et al. Research and development of liquamatic piston fertilizer pump for micro-irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2006, 22(6): 100-103. (in Chinese with English abstract)

    [13] 韓啟彪,吳文勇,劉洪祿,等. 三種水力驅(qū)動(dòng)比例式施肥泵吸肥性能試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2010,26(2):43-47.

    Han Qibiao, Wu Wenyong, Liu Honglu, et al. Experiment on fertilizer suction performance of three hydraulic driven pumps[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(2): 43-47. (in Chinese with English abstract)

    [14] 楊大森. 水動(dòng)比例施肥泵運(yùn)動(dòng)分析及性能優(yōu)化[D]. 鎮(zhèn)江:江蘇大學(xué),2016.

    Yang Dasen. Water-powered Proportional Dosing Pump Movement Analysis and Performance Optimization[D]. Zhenjiang: Jiangsu University, 2016. (in Chinese with English abstract)

    [15] 駱志文,水動(dòng)比例注入泵性能參數(shù)研究[D]. 鎮(zhèn)江:江蘇大學(xué),2017.

    Luo Zhiwen. Research on Performance Parameters of Hydrodynamic Proportional Injection Pump[D]. Zhenjiang: Jiangsu University, 2017. (in Chinese with English abstract)

    [16] Tang P, Li H, Issaka Z, et al. Methodology to investigate the hydraulic characteristics of a water-powered piston type proportional injector used for agricultural chemigation[J]. Applied Engineering in Agriculture, 2018, 34(3): 545-553.

    [17] Tang P, Chen C, Li H. Investigation of hydraulic performance based on response surface methodology for an agricultural chemigation proportional injector[J]. Water, 2020, 12(11): 3155.

    [18] 湯攀,李紅,駱志文,等. 比例施肥泵驅(qū)動(dòng)活塞受力分析及內(nèi)部流動(dòng)模擬與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(23):93-100.

    Tang Pan, Li Hong, Luo Zhiwen, et al. Force analysis of drive piston and simulation and experiment of internal flow for proportional fertilizer pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(23): 93-100. (in Chinese with English abstract)

    [19] 王睿,王文娥,胡笑濤,等. 微灌用施肥泵施肥比例與肥水比對(duì)過濾器堵塞的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(23):117-122.

    Wang Rui, Wang Wene, Hu Xiaotao, et al. Impact of fertilizer proportion and fertilizer-water ratio on clogging of filter by fertilizer pump in microirrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(23): 117-122. (in Chinese with English abstract)

    [20] 吳錫凱. 滴灌系統(tǒng)水力驅(qū)動(dòng)式比例施肥裝置性能影響因素研究[D]. 楊凌:西北農(nóng)林科技大學(xué),2019.

    Wu Xikai. Study on Influence Factors of Performance of Hydraulic Driven Proportional Fertilizer Device in Drip Irrigation System[D]. Yangling: Northwest A&F University, 2019. (in Chinese with English abstract)

    [21] 中國(guó)農(nóng)業(yè)機(jī)械化科學(xué)研究院. 農(nóng)業(yè)灌溉設(shè)備水動(dòng)化肥-農(nóng)藥注入泵:GB/T 19792-2012[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,2012.

    [22] Li H C, Gao F. Investigation on optimising agricultural cultivator openers using quadratic orthogonal rotation regression[J]. Acta Agriculturae Scandinavica, Section B-Soil & Plant Science, 2021, 71(9): 970-979.

    [23] Liu J J, Zhao T, Liu K, et al. Optimization of structure parameters in a coal pyrolysis filtration system based on CFD and quadratic regression orthogonal combination and a genetic algorithm[J]. Engineering Applications of Computational Fluid Mechanics, 2021, 15(1): 815-829.

    [24] 郭嘉明,吳旭東,林詩濤,等. 基于多參數(shù)耦合的蓄冷溫控箱冷板對(duì)流換熱參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(19):228-235.

    Guo Jiaming, Wu Xudong, Lin Shitao, et al. Parameter optimization on convective heat transfer of cold plate for cold storage temperature control box based on multi-parameter coupling[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(19): 228-235. (in Chinese with English abstract)

    [25] 張德勝,祁炳,趙睿杰,等. 海水淡化能量回收透平水力模型優(yōu)化設(shè)計(jì)[J]. 排灌機(jī)械工程學(xué)報(bào),2021,39(7):649-654.

    Zhang Desheng, Qi Bing, Zhao Ruijie, et al. Optimization design of hydraulic model for seawater desalination energy recovery turbine[J]. Journal of Drainage and Irrigation Machinery Engineering, 2021, 39(7): 649-654. (in Chinese with English abstract)

    [26] 諸永定,肖千豪,吳靈輝,等. 基于回歸Kriging代理模型的吸油煙機(jī)用多翼離心風(fēng)機(jī)優(yōu)化[J]. 排灌機(jī)械工程學(xué)報(bào),2022,40(1):62-67.

    Zhu Yongding, Xiao Qianhao, Wu Linghui, et al. Optimization of multi-blade centrifugal fan for range hood based on regression Kriging surrogate model[J]. Journal of Drainage and Irrigation Machinery Engineering, 2022, 40(1): 62-67. (in Chinese with English abstract)

    [27] Fan J L, Wu L F, Zhang F C, et al. Evaluation of drip fertigation uniformity affected by injector type, pressure difference and lateral layout[J]. Irrigation and Drainage, 2017, 66(4): 520-529.

    [28] 李紅,張乾坤,湯攀,等. 閥門調(diào)節(jié)式比例施肥泵性能分析與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(18):34-41.

    Li Hong, Zhang Qiankun, Tang Pan, et al. Performance analysis and test of valve-regulated proportional fertilization pumps[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(18): 34-41. (in Chinese with English abstract)

    [29] Li J S, Meng Y B, Li B. Field evaluation of fertigation uniformity as affected by injector type and manufacturing variability of emitters[J]. Irrigation Science, 2007, 25(2): 117-125.

    Optimization and experiment of the suction piston structure for proportional fertilization pump

    Tang Pan1, Ren Ni2, Yi Zhongyi2, Li Hong1

    (1.212013,;2.210014)

    Fertilization equipment is one of the most important equipment for the fertigation system in recent year. The performance of fertigation pump directly dominates the quality of irrigation and fertilization. Among them, the proportional fertilization pump can be expected to serve as an advanced fertilization equipment. Specifically, the water pressure can be used as the power source in the most areas without electricity. The fertilizer suction piston is one of the vital hydraulic components to promote the accuracy of fertilizer injection. Taking the suction piston as the study object, the quadratic regression orthogonal combination test was carried out to optimize the major structural parameters, according to the operating principle. The injection flow rate was utilized as the response index. A multiple regression model was established for the response index and factors. The variables included the diameter of the lower end of the suction piston, the width of the discharge groove, and the depth of the discharge groove. The results demonstrated that the actual injection ratio of the proportional fertilization pump was lower than the predicted injection ratio under various differential pressures and injection ratios. The difference between the actual and predicted injection ratio was reduced with an increase in the injection ratio, at the differential pressure of less than 0.15 MPa. It infers that a larger injection ratio was preferred to boost the fertilizer injection accuracy. The greatest factors were determined in the injection flow within the factor level range, including the width of the discharge groove, the depth of the discharge groove, and the diameter of the lower end of the suction piston(<0.01). The injection flow rate was significantly influenced by each of the three parameters individually, as well as significantly by their interactions. Once the discharge groove was the same width, the injection flow rate first rose and then fell, as the diameter of the lower end of the suction piston increased. The injection flow rate also increased initially before decreasing, as the width of the discharge groove expanded, when the lower end of the diameter in the suction piston was less than 14 mm. However, the injection flow rate rose, as the width of the discharge groove increased, when the diameter of the lower end of the suction piston was between 14 and 17 mm. The injection flow rate first increased and subsequently declined with the increase in the diameter of the lower end of the suction piston and the depth of the discharge groove. The injection flow rate essentially remained the constant, as the width of discharge groove increased, when the depth of the discharge groove was less than 2.5 mm. Only a little impact was found in the width of the discharge groove on the injection flow at this time. The injection flow rate increased with the discharge groove width, when the depth of the groove was larger than 2.5 mm. An optimal combination was achieved: The diameter of the lower end of the suction piston was 16.6 mm, the width of the discharge groove was 5.5 mm, and the depth of the discharge groove was 3.7 mm. The optimum structure parameters were also obtained for the suction piston flow channel. The structural parameters of the suction piston were optimized for the proportional fertilization pump under various injection ratios and differential pressures. The injection flow rate was higher than that of the prototype proportional fertilization pump, which further decreased the discrepancy between the actual and theoretical injection flow rate. The injection precision increased by 3.33, 1.67, and 7.29 percentage points at the differential pressure of 0.05, 0.10, and 0.15 MPa, respectively. Anyway, the injection precision of proportional fertilization pump was improved significantly to extend the typical working differential pressure range, indicating the better hydraulic performance of the device.

    fertilizer; experiment; fertigation; proportional fertilization pump; parameter optimization; regress equation

    10.11975/j.issn.1002-6819.2022.21.005

    S277.9+4

    A

    1002-6819(2022)-21-0033-09

    湯攀,任妮,易中懿,等. 比例施肥泵吸肥活塞結(jié)構(gòu)優(yōu)化與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(21):33-41.doi:10.11975/j.issn.1002-6819.2022.21.005 http://www.tcsae.org

    Tang Pan, Ren Ni, Yi Zhongyi, et al. Optimization and experiment of the suction piston structure for proportional fertilization pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(21): 33-41. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.21.005 http://www.tcsae.org

    2022-06-10

    2022-10-27

    江蘇省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(BE2021379);江蘇省現(xiàn)代農(nóng)機(jī)裝備與技術(shù)示范推廣項(xiàng)目(NJ2021-24)

    湯攀,博士,副研究員,研究方向?yàn)檗r(nóng)業(yè)灌溉及水肥一體化設(shè)備。Email:tangpan19@163.com

    猜你喜歡
    流槽肥液吸液
    混合肥液電導(dǎo)率預(yù)測(cè)模型構(gòu)建研究
    透析濃縮液吸管蓋聯(lián)合清洗槽改良在血透機(jī)吸液棒中的應(yīng)用
    肥液濃度對(duì)雙點(diǎn)源涌泉根灌土壤入滲特征及水分運(yùn)移的影響
    基于正交試驗(yàn)法下對(duì)噴灌均勻性影響因子的綜合分析
    新型梯級(jí)式浮沫抑除裝置在鹽水精制過程的應(yīng)用
    應(yīng)用電導(dǎo)法實(shí)時(shí)測(cè)定磷酸氫二銨和氯化鉀混合肥液濃度
    高爐出鐵場(chǎng)擺動(dòng)流槽優(yōu)化改造
    山西冶金(2020年4期)2020-09-17 18:10:48
    鋰在熱管反應(yīng)器吸液芯上毛細(xì)作用的理論分析
    預(yù)埋式電加熱流槽開發(fā)與測(cè)試
    Numerical simulation for solid-liquid phase change of metal sodium in combined wick
    天堂av国产一区二区熟女人妻| 小说图片视频综合网站| 悠悠久久av| 看非洲黑人一级黄片| 99在线人妻在线中文字幕| 欧美成人a在线观看| 99久久无色码亚洲精品果冻| 亚洲av.av天堂| 精品欧美国产一区二区三| 国产高清有码在线观看视频| 日本精品一区二区三区蜜桃| 日本精品一区二区三区蜜桃| 超碰av人人做人人爽久久| 精品久久久久久久久av| 国产高清激情床上av| 国产一区二区在线观看日韩| 精品国产三级普通话版| 亚洲乱码一区二区免费版| av天堂中文字幕网| 中文资源天堂在线| a级一级毛片免费在线观看| 国产 一区 欧美 日韩| 久久精品国产亚洲av天美| 51国产日韩欧美| 在线免费十八禁| 亚洲国产色片| 欧美xxxx性猛交bbbb| 别揉我奶头~嗯~啊~动态视频| 插阴视频在线观看视频| 九色成人免费人妻av| 日本与韩国留学比较| 国产精品久久久久久亚洲av鲁大| 精品不卡国产一区二区三区| 看免费成人av毛片| 欧美激情国产日韩精品一区| 赤兔流量卡办理| 99精品在免费线老司机午夜| 一级毛片aaaaaa免费看小| 国产色爽女视频免费观看| 亚洲三级黄色毛片| 成人综合一区亚洲| 日韩欧美 国产精品| 黑人高潮一二区| 欧美zozozo另类| 色视频www国产| 小蜜桃在线观看免费完整版高清| 99国产极品粉嫩在线观看| 亚洲国产欧美人成| 少妇熟女欧美另类| 国产精品99久久久久久久久| eeuss影院久久| 校园人妻丝袜中文字幕| 蜜桃久久精品国产亚洲av| 天天一区二区日本电影三级| 国产一区二区在线观看日韩| 免费人成视频x8x8入口观看| 久久人妻av系列| 国产视频一区二区在线看| 日韩欧美一区二区三区在线观看| 欧美xxxx黑人xx丫x性爽| 女人被狂操c到高潮| 欧美日韩国产亚洲二区| 少妇被粗大猛烈的视频| 日韩欧美精品免费久久| 国产真实乱freesex| 我要看日韩黄色一级片| 国产高清有码在线观看视频| 成人午夜高清在线视频| 99久久无色码亚洲精品果冻| 国产成人a∨麻豆精品| 嫩草影院新地址| 中国美女看黄片| 99久久无色码亚洲精品果冻| 综合色丁香网| 搡老妇女老女人老熟妇| 69av精品久久久久久| 久久午夜福利片| 丰满的人妻完整版| 久久久欧美国产精品| 天堂网av新在线| 精品国内亚洲2022精品成人| 最近最新中文字幕大全电影3| 99热精品在线国产| 亚洲专区国产一区二区| 99久久九九国产精品国产免费| 欧美区成人在线视频| 国内精品美女久久久久久| 91久久精品电影网| 亚洲真实伦在线观看| 男人舔女人下体高潮全视频| 日韩精品有码人妻一区| 日本精品一区二区三区蜜桃| 久久久久性生活片| 成年版毛片免费区| 精品午夜福利在线看| 给我免费播放毛片高清在线观看| 国产精品人妻久久久久久| 亚洲精品在线观看二区| av在线播放精品| 校园春色视频在线观看| 成人特级黄色片久久久久久久| 亚州av有码| 欧美色视频一区免费| 欧美潮喷喷水| 日本五十路高清| 国产蜜桃级精品一区二区三区| 亚洲综合色惰| 日日摸夜夜添夜夜爱| 午夜激情欧美在线| 精品日产1卡2卡| 十八禁国产超污无遮挡网站| 精品久久久久久久久久免费视频| 日韩在线高清观看一区二区三区| 国产精品一区二区三区四区免费观看 | 亚洲av电影不卡..在线观看| 别揉我奶头 嗯啊视频| 精品一区二区三区av网在线观看| 天堂av国产一区二区熟女人妻| 日本成人三级电影网站| 久久久精品94久久精品| 麻豆av噜噜一区二区三区| 国产探花极品一区二区| 精品乱码久久久久久99久播| 99热这里只有是精品在线观看| 国产探花极品一区二区| 五月玫瑰六月丁香| 最近的中文字幕免费完整| 欧美+日韩+精品| 亚洲精品日韩av片在线观看| 插阴视频在线观看视频| 国产亚洲精品综合一区在线观看| 日韩中字成人| 18禁裸乳无遮挡免费网站照片| 成人美女网站在线观看视频| 国产一区二区三区av在线 | 99久国产av精品| 变态另类成人亚洲欧美熟女| 精品午夜福利在线看| 中文字幕精品亚洲无线码一区| 国产精品99久久久久久久久| 日韩三级伦理在线观看| 麻豆精品久久久久久蜜桃| 人妻制服诱惑在线中文字幕| 久久人人爽人人爽人人片va| 中出人妻视频一区二区| 1000部很黄的大片| 免费搜索国产男女视频| 国产精品电影一区二区三区| 欧美日韩精品成人综合77777| 亚洲va在线va天堂va国产| av专区在线播放| 寂寞人妻少妇视频99o| 麻豆久久精品国产亚洲av| 麻豆乱淫一区二区| 内地一区二区视频在线| 欧美3d第一页| 色av中文字幕| 亚洲在线观看片| 熟妇人妻久久中文字幕3abv| 国产黄片美女视频| 高清午夜精品一区二区三区 | 国产黄片美女视频| 日韩av不卡免费在线播放| 国产私拍福利视频在线观看| 俺也久久电影网| 99热精品在线国产| 色噜噜av男人的天堂激情| 99久久精品国产国产毛片| 国产激情偷乱视频一区二区| 精品久久久久久久久亚洲| 亚洲欧美清纯卡通| 欧美+日韩+精品| 亚洲av.av天堂| 国产美女午夜福利| 网址你懂的国产日韩在线| 国产精品人妻久久久影院| 美女 人体艺术 gogo| 成人高潮视频无遮挡免费网站| 色吧在线观看| 国产v大片淫在线免费观看| 成人亚洲欧美一区二区av| 高清日韩中文字幕在线| 狠狠狠狠99中文字幕| 菩萨蛮人人尽说江南好唐韦庄 | 91久久精品国产一区二区成人| 午夜激情欧美在线| 国产精品一区二区免费欧美| 亚洲,欧美,日韩| 不卡一级毛片| 欧美日韩在线观看h| 变态另类丝袜制服| 六月丁香七月| 欧美xxxx性猛交bbbb| 精品国内亚洲2022精品成人| 免费无遮挡裸体视频| 日本与韩国留学比较| 亚州av有码| 一边摸一边抽搐一进一小说| 欧洲精品卡2卡3卡4卡5卡区| 伦理电影大哥的女人| 欧美bdsm另类| 日日摸夜夜添夜夜爱| 熟女人妻精品中文字幕| 在现免费观看毛片| 免费黄网站久久成人精品| 变态另类成人亚洲欧美熟女| 成人国产麻豆网| 99riav亚洲国产免费| 亚洲国产精品国产精品| 国产黄片美女视频| av国产免费在线观看| 偷拍熟女少妇极品色| 非洲黑人性xxxx精品又粗又长| 99久久精品一区二区三区| 精品熟女少妇av免费看| 国产成人91sexporn| 99久久无色码亚洲精品果冻| 欧美成人精品欧美一级黄| 免费人成在线观看视频色| 免费观看在线日韩| 久久婷婷人人爽人人干人人爱| 成人漫画全彩无遮挡| 午夜爱爱视频在线播放| 国产91av在线免费观看| 大香蕉久久网| 国产精品av视频在线免费观看| 日韩高清综合在线| 亚洲欧美日韩高清在线视频| 91午夜精品亚洲一区二区三区| 亚洲av一区综合| 亚洲av电影不卡..在线观看| 欧美最新免费一区二区三区| 精品乱码久久久久久99久播| 1024手机看黄色片| 亚洲国产日韩欧美精品在线观看| 亚洲精品日韩av片在线观看| 欧美zozozo另类| 亚洲精品456在线播放app| 可以在线观看毛片的网站| 91午夜精品亚洲一区二区三区| 噜噜噜噜噜久久久久久91| 美女 人体艺术 gogo| 日日摸夜夜添夜夜添小说| 人妻夜夜爽99麻豆av| 亚洲一级一片aⅴ在线观看| 在线观看午夜福利视频| 亚洲天堂国产精品一区在线| 白带黄色成豆腐渣| 国产蜜桃级精品一区二区三区| 午夜免费激情av| 晚上一个人看的免费电影| 搡老熟女国产l中国老女人| 国产真实伦视频高清在线观看| 日韩欧美免费精品| 欧美色欧美亚洲另类二区| 97超碰精品成人国产| 亚洲精品亚洲一区二区| 国产欧美日韩精品亚洲av| 久久久久国内视频| 老师上课跳d突然被开到最大视频| 亚洲色图av天堂| 日韩制服骚丝袜av| 中国美女看黄片| 亚洲久久久久久中文字幕| 欧美最新免费一区二区三区| 一a级毛片在线观看| 天美传媒精品一区二区| 亚洲av不卡在线观看| 伦理电影大哥的女人| 国产一区二区在线观看日韩| 99久久无色码亚洲精品果冻| 蜜桃久久精品国产亚洲av| av福利片在线观看| 国产精品国产三级国产av玫瑰| 欧美激情在线99| 精品99又大又爽又粗少妇毛片| 搞女人的毛片| 久久人人精品亚洲av| a级毛片a级免费在线| 久久精品国产99精品国产亚洲性色| 国产探花极品一区二区| 菩萨蛮人人尽说江南好唐韦庄 | 老司机影院成人| 免费人成在线观看视频色| 免费观看在线日韩| 日韩三级伦理在线观看| 性色avwww在线观看| 欧美国产日韩亚洲一区| 国产视频一区二区在线看| 波多野结衣高清作品| 久久国内精品自在自线图片| 99热这里只有是精品50| 亚洲国产高清在线一区二区三| 哪里可以看免费的av片| 亚洲婷婷狠狠爱综合网| 日韩欧美国产在线观看| 三级经典国产精品| 看非洲黑人一级黄片| 日日摸夜夜添夜夜爱| 乱码一卡2卡4卡精品| 亚洲欧美精品自产自拍| 国产国拍精品亚洲av在线观看| 日韩亚洲欧美综合| 午夜视频国产福利| 国产高潮美女av| 久久久久久大精品| 性欧美人与动物交配| 99久久九九国产精品国产免费| 色噜噜av男人的天堂激情| 六月丁香七月| ponron亚洲| 男女边吃奶边做爰视频| 可以在线观看的亚洲视频| 天堂av国产一区二区熟女人妻| 日本与韩国留学比较| 性插视频无遮挡在线免费观看| 在线播放国产精品三级| 黄色日韩在线| 国产精品国产三级国产av玫瑰| 日韩亚洲欧美综合| 在线国产一区二区在线| 色在线成人网| 99国产精品一区二区蜜桃av| 亚洲图色成人| 亚洲国产高清在线一区二区三| 国产亚洲91精品色在线| 人妻少妇偷人精品九色| 联通29元200g的流量卡| av天堂在线播放| а√天堂www在线а√下载| 女的被弄到高潮叫床怎么办| 三级经典国产精品| 丰满人妻一区二区三区视频av| 国产亚洲精品久久久com| 中文字幕久久专区| 国产男靠女视频免费网站| 天美传媒精品一区二区| 国产av一区在线观看免费| 人妻夜夜爽99麻豆av| 有码 亚洲区| a级毛色黄片| 久久久久国产精品人妻aⅴ院| 一进一出好大好爽视频| 国产片特级美女逼逼视频| 亚洲欧美成人精品一区二区| 91av网一区二区| videossex国产| 人妻丰满熟妇av一区二区三区| 亚洲精品日韩在线中文字幕 | 九色成人免费人妻av| 国产成人aa在线观看| 亚洲欧美成人精品一区二区| 97人妻精品一区二区三区麻豆| 亚洲va在线va天堂va国产| av在线老鸭窝| 蜜臀久久99精品久久宅男| 亚洲av第一区精品v没综合| 一本精品99久久精品77| 久久精品91蜜桃| 亚洲精品国产成人久久av| 成人二区视频| 三级男女做爰猛烈吃奶摸视频| 国产高清视频在线观看网站| 乱系列少妇在线播放| 欧美激情久久久久久爽电影| 久久久成人免费电影| 好男人在线观看高清免费视频| 欧美高清成人免费视频www| 婷婷色综合大香蕉| 高清毛片免费观看视频网站| 亚洲av中文字字幕乱码综合| 搡老岳熟女国产| 免费搜索国产男女视频| 久久久久久久亚洲中文字幕| 久久久午夜欧美精品| 最近中文字幕高清免费大全6| 国产乱人偷精品视频| 日韩精品青青久久久久久| 一区二区三区四区激情视频 | 久久精品国产99精品国产亚洲性色| 波野结衣二区三区在线| 亚洲自偷自拍三级| 桃色一区二区三区在线观看| 国产高清不卡午夜福利| 久久精品影院6| 如何舔出高潮| 亚洲专区国产一区二区| 亚洲精品粉嫩美女一区| 午夜精品国产一区二区电影 | 我要看日韩黄色一级片| 如何舔出高潮| 天天躁夜夜躁狠狠久久av| 一区二区三区四区激情视频 | 国产精品一区二区免费欧美| 午夜福利在线在线| 夜夜爽天天搞| 中文字幕久久专区| 国产真实伦视频高清在线观看| av卡一久久| 午夜视频国产福利| 亚洲欧美成人综合另类久久久 | av.在线天堂| 乱人视频在线观看| 高清毛片免费看| 日韩欧美 国产精品| 超碰av人人做人人爽久久| av在线老鸭窝| 午夜福利18| 狠狠狠狠99中文字幕| 久久久久精品国产欧美久久久| 一进一出抽搐gif免费好疼| 免费观看人在逋| 国产色婷婷99| 欧美bdsm另类| 简卡轻食公司| 黄色一级大片看看| 亚洲熟妇中文字幕五十中出| 又爽又黄无遮挡网站| 99视频精品全部免费 在线| 国产精品日韩av在线免费观看| 91精品国产九色| 老熟妇仑乱视频hdxx| 亚洲av熟女| 亚洲av中文字字幕乱码综合| 99热网站在线观看| 国产精品一二三区在线看| 全区人妻精品视频| 日韩欧美免费精品| 精品一区二区三区av网在线观看| 97人妻精品一区二区三区麻豆| 99久久精品国产国产毛片| 欧美日本视频| 高清日韩中文字幕在线| 99热全是精品| 欧美日韩一区二区视频在线观看视频在线 | or卡值多少钱| 亚洲美女搞黄在线观看 | 能在线免费观看的黄片| 在线观看免费视频日本深夜| 天堂√8在线中文| 欧美成人免费av一区二区三区| 身体一侧抽搐| 欧美激情在线99| 国产精品三级大全| 久久精品国产亚洲av天美| 日本色播在线视频| 亚洲av免费高清在线观看| 搞女人的毛片| 九色成人免费人妻av| 亚洲美女搞黄在线观看 | 在现免费观看毛片| 日韩,欧美,国产一区二区三区 | 成人永久免费在线观看视频| 亚洲自偷自拍三级| 91精品国产九色| 老熟妇仑乱视频hdxx| videossex国产| 欧美绝顶高潮抽搐喷水| 国产av不卡久久| 亚洲精品久久国产高清桃花| 国产高清三级在线| av国产免费在线观看| 欧美色视频一区免费| 99久久九九国产精品国产免费| 亚洲内射少妇av| 又爽又黄a免费视频| 最好的美女福利视频网| 午夜福利视频1000在线观看| 久久亚洲国产成人精品v| 精品乱码久久久久久99久播| 亚洲一区高清亚洲精品| 国产一级毛片七仙女欲春2| 免费人成视频x8x8入口观看| 免费av观看视频| 国产精品人妻久久久影院| 午夜福利成人在线免费观看| 午夜福利在线在线| 99在线视频只有这里精品首页| av天堂中文字幕网| 黄色视频,在线免费观看| 国产在线精品亚洲第一网站| 国产真实伦视频高清在线观看| 毛片一级片免费看久久久久| 大又大粗又爽又黄少妇毛片口| 亚洲综合色惰| 中文字幕免费在线视频6| 久久久午夜欧美精品| 国产亚洲精品久久久com| 国产高清不卡午夜福利| 两个人视频免费观看高清| 中文资源天堂在线| 国国产精品蜜臀av免费| 日韩大尺度精品在线看网址| 久久中文看片网| 久久久久精品国产欧美久久久| 九九久久精品国产亚洲av麻豆| 欧美成人精品欧美一级黄| 免费看a级黄色片| 最近2019中文字幕mv第一页| 老熟妇仑乱视频hdxx| 一本一本综合久久| 小说图片视频综合网站| 欧美成人a在线观看| 黄色一级大片看看| 一夜夜www| 国产精品亚洲美女久久久| 97在线视频观看| 联通29元200g的流量卡| 欧美精品国产亚洲| 精品人妻一区二区三区麻豆 | 欧美激情在线99| 欧美绝顶高潮抽搐喷水| 日韩欧美免费精品| 97超碰精品成人国产| 亚洲无线在线观看| 免费一级毛片在线播放高清视频| 日韩三级伦理在线观看| 中文字幕人妻熟人妻熟丝袜美| 看非洲黑人一级黄片| 成人鲁丝片一二三区免费| 日本五十路高清| 非洲黑人性xxxx精品又粗又长| 国产高潮美女av| 欧美色视频一区免费| 夜夜看夜夜爽夜夜摸| 欧美成人一区二区免费高清观看| 婷婷亚洲欧美| 亚洲aⅴ乱码一区二区在线播放| 少妇熟女aⅴ在线视频| 69人妻影院| 国产成人一区二区在线| 乱人视频在线观看| 亚洲国产日韩欧美精品在线观看| 综合色av麻豆| 日韩在线高清观看一区二区三区| 欧美一区二区国产精品久久精品| 日韩一区二区视频免费看| 精品乱码久久久久久99久播| 中文在线观看免费www的网站| 丰满人妻一区二区三区视频av| 又爽又黄a免费视频| 亚洲综合色惰| 亚洲不卡免费看| 成人国产麻豆网| 国产精品久久电影中文字幕| 一a级毛片在线观看| 午夜精品国产一区二区电影 | 日本一二三区视频观看| 亚洲电影在线观看av| 日韩亚洲欧美综合| 成人av一区二区三区在线看| 一进一出抽搐动态| 日韩人妻高清精品专区| 国产中年淑女户外野战色| 男人舔奶头视频| 秋霞在线观看毛片| 蜜臀久久99精品久久宅男| 国产69精品久久久久777片| 日韩欧美精品v在线| 如何舔出高潮| 热99在线观看视频| 亚洲第一电影网av| 国产精品一及| 亚洲欧美日韩高清在线视频| 久久精品综合一区二区三区| 亚洲中文日韩欧美视频| 免费看美女性在线毛片视频| 久久精品国产亚洲av涩爱 | 国产精品久久久久久av不卡| 亚洲专区国产一区二区| 午夜福利视频1000在线观看| 日本-黄色视频高清免费观看| 日本欧美国产在线视频| 亚洲中文字幕一区二区三区有码在线看| 国内精品宾馆在线| 永久网站在线| 能在线免费观看的黄片| 波多野结衣高清无吗| 欧美3d第一页| 日韩中字成人| 精品久久久久久成人av| 久久久久久伊人网av| 99久久精品一区二区三区| 美女xxoo啪啪120秒动态图| 国语自产精品视频在线第100页| 最近在线观看免费完整版| 村上凉子中文字幕在线| 别揉我奶头 嗯啊视频| 亚洲熟妇熟女久久| 精品久久久噜噜| 精品久久国产蜜桃| 国产亚洲精品av在线| 亚洲av第一区精品v没综合| h日本视频在线播放| 男人的好看免费观看在线视频| 永久网站在线| 午夜激情福利司机影院| 午夜影院日韩av| 亚洲av电影不卡..在线观看| 国产精品久久久久久精品电影| 久久人妻av系列| 亚洲18禁久久av| 欧美高清成人免费视频www| 日本a在线网址| 黄色一级大片看看| 一边摸一边抽搐一进一小说| av在线老鸭窝| 国产国拍精品亚洲av在线观看| 天堂av国产一区二区熟女人妻| 最近视频中文字幕2019在线8| 久久鲁丝午夜福利片| 在线观看午夜福利视频| 一区二区三区高清视频在线| 99热全是精品| 黄片wwwwww|