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

    楔形減阻旋耕刀設(shè)計與試驗

    2019-05-24 07:39:18郝建軍于海杰趙建國李建昌馬志凱蔡金金
    農(nóng)業(yè)工程學(xué)報 2019年8期
    關(guān)鍵詞:楔形國標(biāo)耕作

    郝建軍,于海杰,趙建國,李建昌,馬志凱,蔡金金

    ?

    楔形減阻旋耕刀設(shè)計與試驗

    郝建軍,于海杰,趙建國,李建昌,馬志凱,蔡金金

    (河北農(nóng)業(yè)大學(xué)機電工程學(xué)院,保定 071001)

    為解決旋耕整地作業(yè)阻力大、能耗高等問題,該文基于旋耕刀理論受力模型設(shè)計了一種楔形減阻旋耕刀。采用Inventor和HyperMesh軟件分別創(chuàng)建國標(biāo)旋耕刀及楔形減阻旋耕刀的三維模型和切削土壤模型,分析了楔形減阻旋耕刀的應(yīng)力強度,對比了國標(biāo)旋耕刀與楔形減阻旋耕刀的切削阻力。通過田間試驗對比了國標(biāo)旋耕刀與楔形減阻旋耕刀的扭矩、功耗與碎土率。結(jié)果表明:楔形減阻旋耕刀所受最大應(yīng)力為29.49 MPa,遠(yuǎn)小于材料的屈服強度430 MPa,在保證刀身強度的前提下,與國標(biāo)旋耕刀相比,楔形減阻旋耕刀質(zhì)量減輕8.3%;平均切削阻力較國標(biāo)旋耕刀下降10.65%。在相同工況條件下,楔形減阻旋耕刀的平均扭矩為648.916 N·m,較國標(biāo)旋耕刀下降11.35%;楔形減阻旋耕刀的平均功耗為67.3 kW,較國標(biāo)旋耕刀下降9.29%,碎土率提高4%,耐磨性能與國標(biāo)旋耕刀持平,能夠達(dá)到在降低作業(yè)功耗的同時,提高耕作質(zhì)量并保證刀具使用壽命。

    農(nóng)業(yè)機械;土壤;模型;楔形減阻旋耕刀;田間試驗

    0 引 言

    旋耕可一次完成耕地、碎土、整地等作業(yè),達(dá)到犁耙?guī)状胃氐男Ч倚蟮牡乇砥秸?、松軟,一次性使土壤達(dá)到待播狀態(tài),是目前平原地區(qū)采用的主要耕作方式[1-3]。旋耕阻力及能源消耗主要來源于旋耕刀[4]。為減小旋耕阻力,生產(chǎn)中通常采用理論減阻、振動減阻和仿生減阻等方法。理論減阻主要是通過改變旋耕刀的工作參數(shù)、折彎角和刃口角等結(jié)構(gòu),來降低切削阻力[5-7],如陳雪等采用ADAMS軟件對旋耕刀工作參數(shù)進行了優(yōu)化,極大地降低了旋耕刀阻力[8]。蓋超等采用COSMOSwork軟件分析了不同旋耕刀折彎角的切削功耗,結(jié)果表明折彎角在125°~130°時功耗最小[9]。康松林等采用LS-DYNA軟件分析了雙面刃口與單面刃口對切削阻力的影響,結(jié)果表明雙面刃口阻力小于單面刃口阻力[10];振動減阻通過凸輪機構(gòu)的擺動使旋耕刀在切削土壤的同時,通過振動使土壤疏松,從而達(dá)到降低耕作阻力的目的[11-13]。如蔣建東等在旋耕刀軸部位加裝凸輪機構(gòu),利用刀具振動來降低切削阻力,結(jié)果表明合適的震動頻率及幅值能夠?qū)崿F(xiàn)土壤耕作減阻[14];仿生減阻主要是參考蚯蚓、鼴鼠等動物表皮、爪趾等特點,對旋耕刀的材料和結(jié)構(gòu)進行優(yōu)化設(shè)計,從而達(dá)到減少耕作阻力的目的[15-17]。如郭俊等根據(jù)鼴鼠指爪結(jié)構(gòu)優(yōu)化旋耕刀,結(jié)果表明在秸稈覆蓋50%以上時,優(yōu)化的旋耕刀耕作阻力低于國標(biāo)旋耕刀[18]。上述減阻技術(shù)雖在一定程度上都能降低旋耕耕作阻力,但振動減阻需要加裝相應(yīng)裝置,在增加了制造成本的同時也增加了機具故障率;仿生減阻多為被動模仿生物的生物學(xué)特性,結(jié)構(gòu)復(fù)雜,不易加 工[19-21]。而理論減阻通過優(yōu)化改進刀具自身結(jié)構(gòu)參數(shù)實現(xiàn)減阻,不增加機械結(jié)構(gòu)且易成型加工。為此,本文針對典型的IT245型國標(biāo)旋耕刀旋耕阻力大的不足,設(shè)計一種楔形減阻旋耕刀,旨在通過改變旋耕刀表面形狀結(jié)構(gòu)以降低切削阻力。

    1 旋耕刀結(jié)構(gòu)與受力分析

    1.1 國標(biāo)旋耕刀結(jié)構(gòu)、工作過程及受力分析

    根據(jù)其彎轉(zhuǎn)向不同,旋耕刀可分為左旋刀與右旋刀2類。典型的IT245型右旋國標(biāo)旋耕刀結(jié)構(gòu)及其參數(shù)如圖1。旋耕作業(yè)時,旋耕刀側(cè)切刃與正切刃先后切入土壤,通過對土壤進行周期性擠壓、切削、破碎、拋甩,完成旋耕作業(yè)[22-23]。

    如圖2,以國標(biāo)旋耕刀切削部任一截面為研究對象,進行受力分析。假設(shè)旋耕為純切削(主要以刃口貫入與切開阻力為主,以Kostritsyn建立的切削計算模型,未考慮土壤粘附力,只考慮了土壤變形破壞后對刀具所起反力的作用平衡)[24],土壤為均質(zhì)土壤[25],并認(rèn)為不發(fā)生其他形式的土壤破壞[26],刀具勻速切入土壤。當(dāng)考慮摩擦?xí)r,依據(jù)文獻[10],土壤相對旋耕刀的滑動摩擦系數(shù)為0.2,土壤與刃口發(fā)生滑動摩擦,產(chǎn)生摩擦角,且相當(dāng)于把刃口角增加,如圖2c所示,刀刃面上修正后的正壓力1為:

    式中1為修正后刃面上的正壓力,N;為修正前刃面上正壓力,N;為土壤和滑切刃之間的動摩擦角,(°)。

    1. 刀柄 2. 刀頭 3. 主切削刃 4. 副切削刃 5. 國標(biāo)旋耕刀(右旋)6. 土壤 7. 刀座 8.刀軸

    1. Blade handle 2. Blade head 3. Main cutting edge 4.Auxiliary cutting edge 5. National standard rotary blade (right-rotary) 6. Soil 7. Blade holder 8. Blade shaft

    注:刀柄寬=30 mm;刀軸孔直徑=13 mm;刀軸孔邊距=15 mm;旋轉(zhuǎn)半徑=245 mm;切削刃寬度=12 mm;切削刃厚度=4 mm;正切面高度=50 mm;耕作寬幅=50 mm;折彎角=120°;折彎半徑=30 mm;為旋耕刀旋轉(zhuǎn)速度,r·min–1;為拖拉機前進速度,km·h–1。為坐標(biāo)原點。

    Note: Blade handle width=30 mm; Rotary blade hole diameter=13 mm; Rotary blade hole space=15 mm; Radius of gyration=245 mm; Cutting edge width=12 mm; Cutting edge thickness=4 mm; Blade height=50 mm; Tillage width=50 mm; Bending angle=120°; Bend radius=30 mm;is rotation speed of rotary blade, r·min–1;is tractor forward speed, km·h–1.is origin of coordinates.

    圖1 右旋國標(biāo)旋耕刀結(jié)構(gòu)參數(shù)與工作示意圖

    Fig.1 Structural parameters and working sketch of right-handed national standard rotary blade

    依上所述,對國標(biāo)旋耕刀進行受力分析。圖2為國標(biāo)旋耕刀切削部任一截面受力分析圖。設(shè)沿軸向右、軸向上為正方向,旋耕刀沿、方向的受力如下:

    式中為力的個數(shù);X為國標(biāo)旋耕刀切削部截面投影在方向的第個力,N;Y為國標(biāo)旋耕刀切削部截面投影在方向的第個力,N。

    注:x為橫坐標(biāo),y為縱坐標(biāo),T1為刃面剪力(即摩擦力),N;N2為刀身正面的正壓力,N;T2為刀身正面的剪力,N;N3為刀身背面的正壓力,N;T3為刀身背面的剪力,N;F貫為旋耕刀切入土壤的貫入阻力,N;F阻為旋耕刀切入土壤所受到的阻力和(即切削阻力),N;j為動摩擦角,(°)。N、N1為修正前后刃面正壓力,N

    需要說明的是,刃面剪力1在方向上的投影未作用在國標(biāo)旋耕刀上,所以公式(3)未對其進行考慮,當(dāng)考慮摩擦力時,1、2、3均為土壤相對于旋耕刀表面的滑動摩擦力,其滑動摩擦力公式為:

    式中T為刀面及刃面上的第個剪力,N;為滑動摩擦系數(shù);tan;N為刀面及刃面上的第個正壓力,N。將式(4)代入式(3),求得3為:

    將式(5)代入式(2),求得阻為:

    由公式(6)可知,切削阻力阻與刀面及刃面上的正壓力N、滑動摩擦因數(shù)和刃口角有關(guān),N值越大切削阻力越大,表明土壤越不易切削;在其他條件不變的情況下,滑動摩擦因數(shù)越大,切削阻力越大。旋耕刀刃口角過小時,土壤與滑切刃之間摩擦劇烈,土壤與刃口面直接摩擦,使得刃口面發(fā)生變形失效,如果刃口角度過大,旋耕刀的切削阻力也會隨之增加[27]。研究表明,在進行土壤切削時,旋耕刀最小阻力刃口角在20°~45°之間[28],為便于沖壓加工,國標(biāo)旋耕刀對刃口進行45°倒角形成刃口角[29]。

    1.2 楔形減阻旋耕刀設(shè)計及受力分析

    1.2.1 楔形減阻旋耕刀設(shè)計

    由上文可知,旋耕刀的切削阻力主要與刀面及刃面上的正壓力、滑動摩擦因數(shù)和刃口角有關(guān),由于滑動摩擦因數(shù)固定(滑動摩擦因數(shù)與接觸材料和相對運動速度有關(guān),耕作同一地塊時,土壤參數(shù)基本一致,刀具相對運動速度一定),楔形減阻旋耕刀設(shè)計主要通過增大刃面寬度及減小刃口角來降低切削阻力。楔形減阻旋耕刀依據(jù)楔子原理將國標(biāo)旋耕刀刀身變成一個較鋒利的刃口,即楔形刃口,如圖3所示。

    1. 國標(biāo)旋耕刀切削部截面 2. 優(yōu)化部分 3. 楔形減阻旋耕刀切削部截面

    1. Cross section of cutting part of national standardrotary blade 2. Optimized part 3. Cross section of cutting part of wedge drag reduction rotary blade

    注:為國標(biāo)旋耕刀刃面寬度,mm;為楔形減阻旋耕刀刃面寬度,mm;為國標(biāo)旋耕刀刃口角,(°);為楔形減阻旋耕刀刃口角,(°)。

    Note:is edge width of rotary blade, mm;is edge width of wedge drag reduction rotary blade, mm;is edge angle of rotary blade, (°);¢is edge angle of wedge drag reduction rotary blade, (°).

    圖3 楔形減阻旋耕刀設(shè)計原理

    Fig.3 Design principle of wedge drag reduction rotary blade

    如圖3所示,1為國標(biāo)旋耕刀切削部截面,去除優(yōu)化部分2,得到楔形減阻旋耕刀的切削部截面3。楔形減阻旋耕刀的刃面寬度大于國標(biāo)旋耕刀的刃面寬度,刃口角小于國標(biāo)旋耕刀的刃口角,刀柄寬度、刀孔直徑、刀孔距、旋轉(zhuǎn)半徑、刃口厚度、正切面寬度、耕作寬幅、折彎角均與國標(biāo)旋耕刀相同。

    如圖4所示,為楔形減阻旋耕刀(以右旋為例)結(jié)構(gòu)參數(shù)示意圖。

    耕作深度是評價旋耕質(zhì)量的指標(biāo)之一,旋耕刀的旋轉(zhuǎn)半徑及彎折角影響耕作深度。楔形減阻旋耕刀與國標(biāo)旋耕刀的旋轉(zhuǎn)半徑與彎折角相同,因此,評價耕作質(zhì)量時耕作深度不予考慮。楔形減阻旋耕刀刃口角變小,與國標(biāo)旋耕刀相比更容易切入并切碎土壤,因此,本文僅以碎土率評價耕作質(zhì)量。

    注:刃面寬度l'=31~50.09 mm。

    1.2.2 楔形減阻旋耕刀受力分析

    為明確楔形減阻旋耕刀是否滿足減阻要求,對楔形減阻旋耕刀進行受力分析,圖5所示為楔形減阻旋耕刀切削部任一截面的受力情況。

    注:o為坐標(biāo)原點;x為橫坐標(biāo);y為縱坐標(biāo);N1¢為修正后刀刃面上的正壓力,N;T1¢為刃面剪力(即摩擦力),N;N3¢為刀身背面的正壓力,N;T3'為刀身背面的剪力,N;F阻¢為楔形減阻旋耕刀切入土壤所受到的阻力和(即切削阻力),N。

    與1.2.1節(jié)國標(biāo)旋耕刀給定情況相同,其沿、方向的受力為:

    將式(4)代入式(8),求得3¢為:

    再將式(9)代入式(7),求得阻¢為:

    由圖2d、圖5可知,旋耕刀刃口厚度不發(fā)生變化,即貫大小相同,比較阻與阻¢,由式(4)可知,當(dāng)動摩擦因數(shù)一定時,TN正相關(guān),由《機械土壤動力學(xué)》可知,當(dāng)?shù)毒邌挝惠d荷一定(即刀具均勻受力)時,其所受正壓力與承載面積成正比。楔形減阻旋耕刀刀身表面所受壓力小于國標(biāo)旋耕刀,即阻¢<阻,因此楔形減阻旋耕刀能夠達(dá)到降低切削阻力的目的,且更易于切削土壤。

    2 旋耕刀切削運動仿真

    采用HyperMesh軟件分別對國標(biāo)旋耕刀與楔形減阻旋耕刀的切削土壤過程進行仿真,分析其切削過程、應(yīng)力大小及切削阻力,以驗證楔形減阻旋耕刀設(shè)計的合理性。

    2.1 有限元模型建立

    2.1.1 刀具有限元模型

    在Inventor軟件中按照GB/T5669—2008《旋耕機械刀與刀座》中旋耕刀幾何參數(shù)進行三維實體建模,同時對楔形減阻旋耕刀建模,如圖6所示,利用HyperMesh(能夠測量刀具質(zhì)量,楔形減阻旋耕刀質(zhì)量減輕50 g,較國標(biāo)旋耕刀下降8.3%)完成刀具有限元模型的網(wǎng)格劃分,采用5 mm網(wǎng)格對2種旋耕刀進行劃分,國標(biāo)旋耕刀的網(wǎng)格數(shù)量為931個,楔形減阻旋耕刀的網(wǎng)格數(shù)量為670個,參照LS-DYNA-971用戶手冊,刀具單元定義為剛性材料MATL20,2種刀具模型材料均為65 Mn,密度為7.8×10-6kg/mm3,彈性模量為2.1×108kPa,泊松比為0.3。

    a. 國標(biāo)旋耕刀a. National standard rotary bladeb. 楔形減阻旋耕刀b. Wedge drag reduction rotary blade

    2.1.2 土壤-刀具模型

    由于HyperMesh軟件中沒有自帶土壤模型,在仿真過程中先選擇HyperMesh自帶的材料模型MATL1,然后根據(jù)關(guān)鍵字手冊把土壤模型改為MAT193(MAT_drucke_ prager)塑性土壤模型,并結(jié)合華北地區(qū)土壤特性[30]輸入材料參數(shù)。土壤模型的主要參數(shù)取值如表2所示(依據(jù)后文田間試驗土壤實際狀況,秸稈均已粉碎還田,仿真試驗未考慮秸稈狀況)。

    表2 土壤參數(shù)

    圖7所示為土壤-旋耕刀初始模型,土壤模型尺寸為520 mm×200 mm×240 mm,為了縮短計算時間,采用分割法劃分網(wǎng)格,將土壤網(wǎng)格劃分為切削部分(尺寸為4 mm)和未切削部分(尺寸為8 mm)。

    2.2 邊界條件

    對刀具與土壤進行相應(yīng)的條件約束,第一步是對土壤模型的四周與底面進行全約束,保證在旋耕刀切削土壤過程中土壤能夠保持恒定不動。第二步定義旋耕刀運動,使旋耕刀能夠圍繞旋轉(zhuǎn)中心做旋轉(zhuǎn)運動且有向前行進的速度。第三步對旋耕刀與土壤接觸添加定義,選用面面接觸AUTOMATIC_SUFACE_TO_SUFACE中的Eroding侵蝕接觸,這種定義能夠真實反映旋耕刀與土壤進行接觸情況,設(shè)定滑動摩擦因數(shù)為0.2。

    2.3 刀具切削過程及應(yīng)力分析

    刀具切削土壤時,隨著載荷的增大,土壤先發(fā)生彈性變形,然后發(fā)生塑性變形,隨后進入屈服狀態(tài)并出現(xiàn)硬化現(xiàn)象。隨著刀具的切入,應(yīng)力繼續(xù)增大,土壤發(fā)生失效,進而破裂[31]。

    圖8為國標(biāo)旋耕刀在LS-DYNA中模擬仿真切削土壤過程。切削參數(shù)依據(jù)《GB/T 5668-2008 旋耕機》[32]選?。盒掇D(zhuǎn)速300 r/min(5 r/s),前進速度4.10 km/h(1.1 m/s),耕深150 mm。

    由圖8可知,0s時,旋耕刀與土壤進行接觸;400s時,旋耕刀先沿側(cè)切刃由近及遠(yuǎn)切入土壤,土壤受到剪切力的作用,開始發(fā)生變形;750~3 000s時,正切刃切入土壤,在旋耕刀的剪切與擠壓共同作用下,土壤沿著旋耕刀正面發(fā)生破壞變形。

    由圖9a可知,國標(biāo)旋耕刀工作時所受的最大應(yīng)力為35.68 MPa,其遠(yuǎn)小于刀具材料的屈服強度430 MPa,驗證了國標(biāo)旋耕刀設(shè)計合理性;由圖9b可知,楔形減阻旋耕刀工作時所受最大應(yīng)力(29.49 MPa)小于國標(biāo)旋耕刀最大應(yīng)力(35.68 MPa),滿足強度設(shè)計要求。另外,在滿足強度設(shè)計的要求下,與國標(biāo)旋耕刀相比,楔形減阻旋耕刀質(zhì)量較小。

    圖8 旋耕刀切削土壤過程

    圖9 2種旋耕刀應(yīng)力圖

    2.4 切削阻力分析

    依照上文,已對2種旋耕刀的材料參數(shù)進行定義與邊界約束,其工作參數(shù)為前進速度1.1 m/s,旋轉(zhuǎn)速度5 r/s,耕作深度150 mm,定義旋耕刀為正旋,以切削阻力為指標(biāo),在HyperMesh軟件中對2種旋耕刀的切削阻力進行分析,導(dǎo)出K文件,將K文件導(dǎo)入LS-DYNA中進行計算,對比2種旋耕刀切削阻力。

    圖10為國標(biāo)旋耕刀與楔形減阻旋耕刀切削阻力仿真結(jié)果對比。由圖10可知,旋耕刀由0s開始,切削阻力先上升,達(dá)到最大切削阻力之后再下降,直至旋耕刀切出土壤。模擬仿真的切削阻力趨勢與參考文獻[10,33]相似,證明了其研究可行性。國內(nèi)學(xué)者葛云[34]、陳雪等[8]對國標(biāo)旋耕刀切削阻力做了研究,得出國標(biāo)旋耕刀的最大切削阻力大概為450~500 N,本文中,國標(biāo)旋耕刀最大切削阻力為530 N,平均切削阻力為375.5 N,楔形減阻旋耕刀最大切削阻力為510 N,平均切削阻力為335.5 N,楔形減阻旋耕刀平均切削阻力比國標(biāo)旋耕刀平均切削阻力降低了10.65%。

    圖10 旋耕刀切削阻力仿真結(jié)果對比

    3 田間試驗

    為了進一步驗證楔形減阻旋耕刀的減阻降耗效果與耕作質(zhì)量,與國標(biāo)旋耕刀進行了田間對比試驗,以旋轉(zhuǎn)扭矩[35]、功率消耗及碎土率為指標(biāo),評價楔形減阻旋耕刀的減阻降耗效果與耕作質(zhì)量。

    3.1 試驗方法

    3.1.1 設(shè)備安裝及扭矩傳感器標(biāo)定

    如圖11所示,將試驗旋耕刀安裝在旋耕機刀軸上,扭矩傳感器在拖拉機后輸出軸上與旋耕機相連。扭矩-轉(zhuǎn)速-功率顯示儀與扭矩傳感器相連,并連接至電腦,將扭矩、功率信號通過數(shù)據(jù)線傳遞給電腦,用M400數(shù)據(jù)分析軟件進行實時記錄。

    1. 拖拉機2. 扭矩-轉(zhuǎn)速-功率顯示儀3. 電腦4. 扭矩傳感器5. 旋耕機6. 旋耕刀7. 輸出軸

    扭矩傳感器標(biāo)定結(jié)果如圖12所示。根據(jù)標(biāo)定結(jié)果,靜止時,扭矩傳感器的頻率在10 000~10 040 Hz之間波動,扭矩傳感器存在的誤差為0.4%(即靜止時,傳感器頻率波動值與初始頻率的比值),誤差在可接受范圍之內(nèi),滿足試驗要求。

    圖12 扭矩傳感器標(biāo)定結(jié)果

    3.1.2 試驗設(shè)置

    試驗于2018年9月26號在河北省定州市內(nèi)化村進行。

    試驗地為壤土,地勢平坦,雖有部分秸稈,但均已粉碎還田,對土壤旋耕功耗不造成影響,0~150 mm土層內(nèi)平均硬度為6.5 kg/cm2,平均含水率為15%。

    將楔形減阻旋耕刀和國標(biāo)旋耕刀分2次各自安裝在河北雙天機械制造有限公司生產(chǎn)的1GKN- 220A1旋耕機上,依次進行田間試驗。

    選取地勢平坦的100 m×30 m地塊,去除地頭地尾各10 m,在中間80 m的旋耕機穩(wěn)定運行段進行5次重復(fù)試驗,測試旋轉(zhuǎn)扭矩和功率消耗,取平均值。依據(jù)上文作業(yè)參數(shù)為:拖拉機前進速度取4.10 km/h,后輸出轉(zhuǎn)速540 r/min,刀輥轉(zhuǎn)速為300 r/min,耕作深度為150 mm。試驗裝備及實物圖如圖13所示。

    3.2 試驗結(jié)果與分析

    每次試驗取穩(wěn)定段測試數(shù)據(jù),每次記錄50 s數(shù)據(jù),2種旋耕刀的扭矩和功耗測試結(jié)果如圖14所示,質(zhì)地復(fù)雜,旋耕扭矩波動比較大,2種旋耕刀扭矩趨勢大致相似,符合試驗要求。

    由圖14a可知,國標(biāo)旋耕刀的平均扭矩值為732.038 N×m,楔形減阻旋耕刀的平均扭矩值為648.916 N×m,楔形減阻旋耕刀較國標(biāo)旋耕刀平均扭矩下降11.35%。由前文可知,拖拉機后輸出轉(zhuǎn)速固定,結(jié)合扭矩計算公式[35],扭矩值能夠?qū)倪M行表達(dá)且與功耗呈正線性相關(guān),楔形減阻旋耕刀在作業(yè)過程中平均扭矩優(yōu)于國標(biāo)旋耕刀,進一步驗證了楔形減阻旋耕刀節(jié)能降耗的可行性。

    由圖14b可知,楔形減阻旋耕刀功耗整體優(yōu)于國標(biāo)旋耕刀,國標(biāo)旋耕刀平均功耗為74.2 kW,楔形減阻旋耕刀平均功耗為67.3 kW,楔形減阻旋耕刀較國標(biāo)旋耕刀平均功耗降低了9.29%。

    切削阻力與旋轉(zhuǎn)扭矩都能夠表示旋耕功耗,由仿真試驗與田間試驗可知,楔形減阻旋耕刀的切削阻力與旋轉(zhuǎn)扭矩(均與功耗呈正線性相關(guān))均較國標(biāo)旋耕刀有所下降,功耗的下降雖然存在一定誤差,但都在可接受范圍之內(nèi)。

    圖13 試驗裝備及實物圖

    如圖15所示,對試驗土壤采樣并計算碎土率以評價旋耕刀耕作質(zhì)量。在已旋耕的500 mm×500 mm耕層進行土壤采樣5次,分別測試采樣土壤質(zhì)量總及土壤中直徑大于40 mm的土塊質(zhì)量>40 mm,根據(jù)公式(13)計算碎土率R。

    計算結(jié)果表明,國標(biāo)旋耕刀與楔形減阻旋耕刀的平均碎土率分別為75%和79%。楔形減阻旋耕刀碎土率較國標(biāo)旋耕刀提高了4%,耕作質(zhì)量優(yōu)于國標(biāo)旋耕刀。

    2018年9月30日,在河北保定龐口試驗田對國標(biāo)旋耕刀與楔形減阻旋耕刀的耐用性進行了田間對比試驗。試驗結(jié)果表明:經(jīng)過連續(xù)不間斷33 hm2耕作后,楔形減阻旋耕刀未發(fā)生崩刃和斷裂現(xiàn)象,且其磨損程度與國標(biāo)旋耕刀相當(dāng):國標(biāo)旋耕刀的平均磨損量為98.5 g,楔形減阻旋耕刀的平均磨損量為97.3 g,磨損量相差1.3%。按規(guī)定旋耕刀磨損量超過100 g即無法繼續(xù)使用。上述結(jié)果表明,在保證與國標(biāo)旋耕刀耐用性基本持平的情況下,楔形減阻旋耕刀能提高耕作質(zhì)量,降低作業(yè)功耗。

    圖14 2種旋耕刀的扭矩和功耗測試結(jié)果

    圖15 碎土率測算

    4 結(jié) 論

    1)基于楔子理論與刀具受力模型設(shè)計一種楔形減阻旋耕刀,分析國標(biāo)旋耕刀與楔形減阻旋耕刀受力,可知切削阻力主要與動摩擦因數(shù)、刃口角與刃口角表面正壓力相關(guān);且國標(biāo)旋耕刀切削阻力大于楔形減阻旋耕刀。

    2)以旋耕刀與楔形減阻旋耕刀為研究對象,運用Inventor軟件建立刀具三維模型,基于HyperMesh-LS- DYNA軟件對刀具進行了切削過程及切削阻力分析。結(jié)果表明:模擬動態(tài)仿真能夠客觀的反映旋耕刀切削土壤過程;國標(biāo)旋耕刀最大應(yīng)力為35.68 MPa,楔形減阻旋耕刀的最大應(yīng)力為29.49 MPa,其應(yīng)力遠(yuǎn)遠(yuǎn)小于刀具材料的屈服強度430 MPa,在質(zhì)量減輕8.3%的情況下,保證了刀身強度。國標(biāo)旋耕刀的最大切削阻力為530 N,楔形減阻旋耕刀最大切削阻力為510 N,楔形減阻旋耕刀平均切削阻力比國標(biāo)旋耕刀平均切削阻力低10.65%,驗證了楔形減阻旋耕刀設(shè)計可行性。

    3)田間試驗結(jié)果表明:在標(biāo)準(zhǔn)工況條件下,楔形減阻旋耕刀的平均扭矩較國標(biāo)旋耕刀下降11.35%,平均功耗較國標(biāo)旋耕刀下降9.29%,平均碎土率較國標(biāo)刀提高4%,楔形減阻旋耕刀與國標(biāo)旋耕刀磨損量相差1.3%,耐用性基本保持一致。進一步驗證了楔形減阻旋耕刀設(shè)計的合理性。

    [1] 車剛,張偉,萬霖,等. 基于滅茬圓盤驅(qū)動旋耕刀多功能耕整機設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(20): 34-40. Che Gang, Zhang Wei, Wan Lin, et al. Design and experiment of multifunctional tillage machine with driven bent blade by stubble ploughing disk[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2012, 28(20): 34-40. (in Chinese with English abstract)

    [2] Temesgen M, Hoogmoed W B, Rockstrom J, et al. Conservation tillage implements and systems for smallholder farms insemi-arid Ethiopia[J]. Soil and Tillage Research, 2009, 104(1): 185-191.

    [3] 熊平原,楊洲,孫志全,等. 基于離散元法的旋耕刀三向工作阻力仿真分析與試驗[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 tne Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 113-121. (in Chinese with English abstract)

    [4] 方會敏,姬長英,張慶怡,等. 基于離散元法的旋耕刀受力分析[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(21):54-59. Fang Huimin, Ji Changying, Zhang Qingyi, et al. Force analysis of rotary blade based on distinct element method[J]. Transactions of tne Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(21): 54-59. (in Chinese with English abstract)

    [5] 張居敏,周勇,夏俊芳,等. 旋耕埋草機螺旋橫刀的數(shù)學(xué)建模與參數(shù)分析[J]. 農(nóng)業(yè)工程學(xué)報,2013,29(1):18-25. Zhang Jumin, Zhou Yong, Xia Junfang, et al. Mathematical modeling and analysis of helical blade for stul rotary tiller[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(1): 18-25. (in Chinese with English abstract)

    [6] 賈洪雷,汲文峰,韓偉峰,等. 旋耕-碎茬通用刀片結(jié)構(gòu)參數(shù)優(yōu)化試驗[J]. 農(nóng)業(yè)機械學(xué)報,2009,40(7):45-50. Jia Honglei, Ji Wenfeng, Han Weifeng, et al. Optimization experiment of structure parameters rototilling and stubble breaking universal blade[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(7): 45-50. (in Chinese with English abstract)

    [7] 湯楚宙,謝方平,向衛(wèi)兵,等. 自推進耕耘機械研究現(xiàn)狀與分析[J]. 農(nóng)業(yè)機械學(xué)報,2001,32(5):112-114. Tang Chuzhou, Xie Fangping, Xiang Weibing, et al. Current situation of study on self-propelled cultivation machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2001, 32(5): 112-114. (in Chinese with English abstract)

    [8] 陳雪,張周,黃化剛,等. 基于ADAMS旋耕刀工作參數(shù)優(yōu)化及ANSYS仿真分析[J]. 現(xiàn)代農(nóng)業(yè)科技,2018(2):167-170. Chen Xue, Zhang Zhou, Huang Huagang, et al. Optimization of rotary blade working parameters based on ADAMS and ANSYS simulation analysis[J]. Modern Agricultural Sciences and Technology, 2018(2): 167-170. (in Chinese with English abstract)

    [9] 蓋超,董玉平. 基于COSMOS的還田機械旋耕刀彎折角優(yōu)化[J]. 農(nóng)機化研究,2011,33(3):30-33. Gai Chao, Dong Yuping. Optimization of bending angle of rotary machine rotary tilling blade based on COSMOS[J]. Journl of Agricultural Mechanization Research, 2011, 33(3): 30-33. (in Chinese with English abstract)

    [10] 康松林. 微耕機刀具的有限元分析及優(yōu)化[D]. 重慶:重慶理工大學(xué),2015. Kang Songlin. Finite Element Analysis and Optimization of the Tool of Rotary Tillers[D]. Chongqing: Chongqing University of Technology, 2015. (in Chinese with English abstract)

    [11] Butson M J, Mac I D. Vibratory cutting soil tank studies of draught and power requirements[J]. Journal of Terramechanics, 2004, 40 (1): 186-190.

    [12] 殷涌光,程悅蓀,李俊明. 振動式二維切削土壤減小阻力機理[J]. 農(nóng)業(yè)機械學(xué)報,1992,23(2):11-16. Yin Yongguang, Cheng Yuesun, Li Junming. The reasons of resistance reduction for vibaatory cutting soil in two directions[J]. Transactions of the Chinese Society for Agricultural Machinery, 1992, 23(2): 11-16. (in Chinese with English abstract)

    [13] 高潔,趙穎娣,Willem Hoogmoed,等. 基于ALE有限元仿真的土壤切削振動減阻[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(增刊1):33-38. Gao Jie, Zhao Yingdi, Willem Hoogmoed, et al. Numerical simulation on resistance reduction of soil vibratory tillage using ALE equation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(Supp.1): 33-38. (in Chinese with English abstract)

    [14] 蔣建東,高潔,趙穎娣,等. 土壤旋切振動減阻的有限元分析[J]. 農(nóng)業(yè)機械學(xué)報,2012,43(1):58-62. Jiang Jiandong, Gao Jie, Zhao Yingdi, et al. Finite element simulation and analysis on soil rotary tillage with external vibration excitation[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(1): 58-62. (in Chinese with English abstract)

    [15] 陳秉聰,任露泉,李安琪,等. 蚯蚓體表液收取方法的初步研究[J]. 農(nóng)業(yè)工程學(xué)報,1990,6(2):7-12. Chen Bingcong, Ren Lunquan, Li Anqi, et al. Initial study on the method of collecting the body surface fluid of erathworms[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1990, 6(2): 7-12. (in Chinese with English abstract)

    [16] 李曉鵬. 馬鈴薯耦合仿生挖掘及其減阻研究[D]. 成都:西華大學(xué),2018. Li Xiaopeng. Research on Potato Coupled-Bionic Digging and its Drag Reduction[D]. Chengdu:Xihua University. 2018. (in Chinese with English abstract)

    [17] 汲文峰,賈洪雷,佟金. 旋耕-碎茬仿生刀片田間作業(yè)性能的試驗研究[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(12):24-30. Ji Wenfeng, Jia Honglei, Tong Jin. Experiment on working performance of bionic blade for soil-rototilling and stubble-breaking[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(12): 24-30. (in Chinese with English abstract)

    [18] 郭俊,張慶怡,Muhammad Sohail Memon,等. 仿鼴鼠足趾排列的旋耕-秸稈粉碎鋸齒刀片設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(6):43-50. Guo Jun, Zhang Qingyi, Muhammad Sohail Memon, et al. Design and experiment of bionic mole’s toe arrangement serrated blade for soil-rototilling and straw-shattering[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(6): 43-50. (in Chinese with English abstract)

    [19] 薛維良. 仿生減阻雙圓盤開溝器設(shè)計與研究[D]. 長春:吉林大學(xué),2017. Xue Weiliang. Design and Research on Bionic Double Disc Opner with Drag Reducing[D]. Changchun: Jilin University, 2017. (in Chinese with English abstract)

    [20] 楊曉東,任露泉. 形體減阻類型、減阻機理與仿生[J]. 農(nóng)業(yè)機械學(xué)報,2003,18(1):130-133. Yang Xiaodong, Ren Luquan. Typers and mechanisms of shape drag reduction[J]. Transactions of the Chinese Society for Agricultural Machinery, 2003,18(1): 130-133. (in Chinese with English abstract)

    [21] 張清珠. 仿生幾何結(jié)構(gòu)表面土壤鎮(zhèn)壓輥[D]. 長春:吉林大學(xué),2014.Zhang Qingzhu. Soil Press Roller with Bionically Geometrically Structured Surfaces[D]. Changchun: Jilin University, 2014. (in Chinese with English abstract)

    [22] 丁為民,徐志剛,汪小函. 斜置旋耕刀滑切角及其方程[J]. 農(nóng)業(yè)工程學(xué)報,2002,18(3):49-53. Ding Weimin, Xu Zhigang, Wang Xiaohan. Grass sliding cutting angles and their equations of oblique rotary blades[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2002, 18(3): 49-53. (in Chinese with English abstract)

    [23] 劉永清,桑正中. 潛土逆轉(zhuǎn)旋耕刀數(shù)學(xué)模型及參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報,2000,16(4):88-91. Liu Yongqing, Sang Zhengzhong. Mathematical model of submerged reverse-rotary tiller and parameters optimization[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2000, 16(4): 88-91. (in Chinese with English abstract)

    [24] 曾德超. 機械土壤動力學(xué)[M]. 北京:北京科學(xué)技術(shù)出版社. 1995.

    [25] 趙家豐,汪偉,孫中興,等. 均質(zhì)土壤承壓下陷模型改進及驗證[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(21):60-66. Zhao Jiafeng, Wang Wei, Sun Zhongxing, et al. Improvement and verification of pressure-sinkage model in homogeneous soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(21): 60-66. (in Chinese with English abstract)

    [26] 熊平原,楊洲,孫志全,等. 旋耕刀三向工作阻力試驗及作業(yè)參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(19):51-58. Xiong Pingyuan, Yang Zhou, Sun Zhiquan, et al. Experiment on three-axis working resistances of rotary blade and working parameters optimization[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(19): 51-58. (in Chinese with English abstract)

    [27] Gill , Vanden Berg. Soil dynamics in tillage and traction[J]. Journal of Terramechanics, 1968, 5(4): 65-66.

    [28] 鄭侃,何進,李洪文,等. 基于離散元深松土壤模型的折線破土刃深松鏟研究[J]. 農(nóng)業(yè)機械學(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)

    [29] 朱留憲. 基于SPH算法的微耕機旋耕刀有限元仿真與優(yōu)化[D]. 重慶:西南大學(xué),2012. Zhu Liuxian. Finite Elements Simulation and Optimization of Rotary Blade of Mini-Tiller Based on SPH Algrithm[D]. Chongqing: Southwest University, 2012. (in Chinese with English abstract)

    [30] 盧彩云,何進,李洪文,等. 基于SPH算法的平面刀土壤切削過程模擬[J]. 農(nóng)業(yè)機械學(xué)報,2014,45(8):134-139. Lu Caiyun, He Jin, Li Hongwen, et al. Simulation of soil cutting process by plane blade based on SPH method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(8): 134-139. (in Chinese with English abstract)

    [31] 顧耀權(quán),賈洪雷,郭慧,等. 滑刀式開溝器設(shè)計與試驗[J]. 農(nóng)業(yè)機械學(xué)報,2013,44(2):38-42. Gu Yaoquan, Jia Honglei, Guo Hui, et al. Design and experiment of sliding knife furrow openner[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(2): 38-42. (in Chinese with English abstract)

    [32] 中國機械工業(yè)聯(lián)合會. 旋耕機:GB/T 5668-2008 [S]. 北京:中國標(biāo)準(zhǔn)出版社. 2009.

    [33] 劉謙文,楊有剛. 基于ANSYS/LS-DYNA的旋耕刀強度和功耗研究[J]. 中國農(nóng)機化學(xué)報,2017,38(6):16-19. Liu Qianwen, Yang Yougang. Research on strength and power consumption of rotary blade based on ANSYS/LS- DYNA[J]. Journal of Chinese Agricultural Mechanization , 2017, 38(6): 16-19. (in Chinese with English abstract)

    [34] 葛云,吳雪飛,王磊,等. 基于ANSYS微型旋耕機旋耕彎刀的應(yīng)力仿真[J]. 石河子大學(xué)學(xué)報:自然科學(xué)版,2007,25(5):627-629. Ge Yun, Wu Xuefei, Wang Lei, et al. The strss simulation for the rotary blade of self propelled rotary cultivator based on ANSYS[J]. Journal of Shihezi University: Natural Science, 2007, 25(5): 627-629. (in Chinese with English abstract)

    [35] 方會敏. 基于離散元法的秸稈—土壤—旋耕刀相互作用機理研究[D]. 南京:南京農(nóng)業(yè)大學(xué),2016. Fang Huimin. Research on the Straw-Soil-Rotary Blade Interaction Using Discrete Element[D]. Nanjing: Nanjing Agricultural University, 2016. (in Chinese with English abstract)

    Design and test of wedge drag reduction rotary blade

    Hao Jianjun, Yu Haijie, Zhao Jianguo, Li Jianchang, Ma Zhikai, Cai Jinjin

    (071001)

    In order to solve the problems of high resistance and high power consumption of rotary tillage and land preparation, the rotary blade is studied based on theoretical analysis, simulation and field experimentin this paper. First, a wedge drag reduction rotary blade is designed based on the theoretical force model of cutting tools in mechanical soil dynamics. Second , the three-dimensional model of the national rotary blade and the wedge drag reduction rotary blade is established by using the Inventor software and the simulation experiment is carried on. The cutting soil model of 2 kinds of rotary blades is established by using HyperMesh software, and the material parameters, boundary constraints and working parameters are defined. On this basis, the soil cutting process, stress and cutting resistance are analyzed. It is proved that the wedge drag reduction rotary blade design is reasonable. In order to verify whether the wedge drag reduction rotary blade meets the strength requirement, the stress of the rotary blades is analyzed, and the cutting resistance of the national rotary blade and the wedge drag reduction rotary blade are compared. Third, field comparative experiments are carried out on national rotary blade and wedge drag reducing rotating blade with the evaluation indexes of torque, power consumption and soil breakage. The results show that the cutting resistance of the wedge drag reduction rotary blade is smaller than that of the national rotary blade; the maximum stress of the national rotary blade is 35.68 MPa and the maximum stress of the wedge drag reduction rotary blade is 29.49 MPa, which is less than that of the national rotary blade and far less than the yield strength of the material 430 MPa. On the premise of ensuring the strength of the blade body, the weight of the wedge drag reduction rotary blade is reduced by about 50 g, which is about 8.3% lower than that of the national rotary blade. The maximum cutting resistance of the national rotary blade is 530 N, the average cutting resistance is 375.5 N, and the maximum cutting resistance of the wedge drag reduction rotary blade is 510 N, the average cutting resistance is 335.5 N, the average cutting resistance of the wedge drag reduction rotary blade is 10.65% lower than that of the national rotary blade, which proves the feasibility of the wedge drag reduction rotary blade designed; The average torque of the wedge drag reduction rotary blade is 11.35% lower than that of the national rotary blade. In the working process of rotary blades, the average power consumption of the national rotary blade is 74.2 kW, the average power consumption of the wedge drag reduction rotary blade is 67.3 kW, and the average power consumption of the wedge drag reduction rotary blade is 9.29% lower than that of the national rotary blade. The average breaking rate of soil of the national rotary blade is 75%, and the average breaking rate of wedge drag reduction rotary blade is 79%, which is 4% higher than that of the national rotary blade; After continuous 33 hm2tillage, rotary blades did not break, and the average wear of the national rotary blade is 98.5 g, the average wear of the wedge shaped drag reduction rotary blade is 97.3 g, the wear difference is 1.3%, and the wear degree is the same. The wedge drag reduction rotary blade can effectively solve the problems of high resistance and high power consumption of rotary blade, which is helpful to accelerate the development of arable land machinery, improve the mechanization level of cultivated land, and reduce the production cost of crops.

    agricultural machinery; soil; models; wedge drag reduction rotary blade; field experiment

    2019-01-09

    2019-03-20

    國家重點研發(fā)計劃(2017YFD0300907)

    郝建軍,博士,教授,博士生導(dǎo)師,從事農(nóng)機裝備設(shè)計與制造。Email:hjjpaper@163.com

    10.11975/j.issn.1002-6819.2019.08.007

    S222.12+9

    A

    1002-6819(2019)-08-0055-10

    郝建軍,于海杰,趙建國,李建昌,馬志凱,蔡金金.楔形減阻旋耕刀設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(8):55-64. doi:10.11975/j.issn.1002-6819.2019.08.007 http://www.tcsae.org

    Hao Jianjun, Yu Haijie, Zhao Jianguo, Li Jianchang, Ma Zhikai, Cai Jinjin. Design and test of wedge drag reduction rotary blade[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(8): 55-64. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.08.007 http://www.tcsae.org

    猜你喜歡
    楔形國標(biāo)耕作
    History of the Alphabet
    鋼絲繩楔形接頭連接失效分析與預(yù)防
    Eight Surprising Foods You’er Never Tried to Grill Before
    凝心聚力,共同迎接“新國標(biāo)”時代
    中國自行車(2018年3期)2018-04-18 07:16:33
    腹腔鏡下胃楔形切除術(shù)治療胃間質(zhì)瘤30例
    車輛碰撞預(yù)警系統(tǒng)有望進入國標(biāo)
    汽車與安全(2016年5期)2016-12-01 05:22:15
    耕作深度對紫色土坡地旋耕機耕作侵蝕的影響
    玉米保護性耕作的技術(shù)要領(lǐng)
    中小學(xué)塑膠跑道“新國標(biāo)”立項公示
    草地耕作技術(shù)在澳大利亞的應(yīng)用
    土壤與作物(2015年3期)2015-12-08 00:46:58
    国国产精品蜜臀av免费| 亚洲伊人久久精品综合| 最近2019中文字幕mv第一页| 国产亚洲5aaaaa淫片| 亚洲无线观看免费| 欧美丝袜亚洲另类| 久久久成人免费电影| 国产亚洲精品av在线| 亚洲国产色片| 2021少妇久久久久久久久久久| 免费大片18禁| 久久久久久国产a免费观看| 五月玫瑰六月丁香| 日韩欧美三级三区| 国产 亚洲一区二区三区 | 亚洲精品日韩av片在线观看| 亚洲成人精品中文字幕电影| 国产精品综合久久久久久久免费| 插逼视频在线观看| 不卡视频在线观看欧美| 国产乱人偷精品视频| 禁无遮挡网站| 久久精品综合一区二区三区| 日产精品乱码卡一卡2卡三| 白带黄色成豆腐渣| av在线老鸭窝| 国产人妻一区二区三区在| 久久综合国产亚洲精品| 亚洲怡红院男人天堂| 精品一区二区三区人妻视频| 亚洲性久久影院| 久久久a久久爽久久v久久| 99热6这里只有精品| 日韩欧美精品免费久久| 精品一区在线观看国产| 亚洲综合色惰| 97人妻精品一区二区三区麻豆| 老女人水多毛片| 好男人在线观看高清免费视频| av在线天堂中文字幕| 色综合站精品国产| 97超视频在线观看视频| 亚洲精品乱久久久久久| 国产成年人精品一区二区| 精品久久国产蜜桃| 成人特级av手机在线观看| 欧美性感艳星| 黑人高潮一二区| 国产成人精品久久久久久| 久久久久国产网址| 色5月婷婷丁香| 亚洲,欧美,日韩| 国产成人午夜福利电影在线观看| 丝瓜视频免费看黄片| 国产v大片淫在线免费观看| 亚洲成人一二三区av| 欧美成人精品欧美一级黄| 蜜桃亚洲精品一区二区三区| 亚洲成人一二三区av| 狠狠精品人妻久久久久久综合| 国产v大片淫在线免费观看| 九九爱精品视频在线观看| 别揉我奶头 嗯啊视频| av黄色大香蕉| 亚洲婷婷狠狠爱综合网| eeuss影院久久| 男女啪啪激烈高潮av片| 熟妇人妻久久中文字幕3abv| 久久99热6这里只有精品| 一级毛片久久久久久久久女| 97超视频在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 边亲边吃奶的免费视频| 永久网站在线| 男女边摸边吃奶| 国产av在哪里看| 五月伊人婷婷丁香| 亚洲怡红院男人天堂| 亚洲一级一片aⅴ在线观看| 亚洲18禁久久av| 777米奇影视久久| 在线免费观看的www视频| 欧美日韩在线观看h| 大陆偷拍与自拍| 亚洲欧美成人精品一区二区| 尤物成人国产欧美一区二区三区| 亚洲精品第二区| 日韩欧美精品免费久久| 嫩草影院新地址| av一本久久久久| 久久久久久国产a免费观看| 国产黄a三级三级三级人| 爱豆传媒免费全集在线观看| 少妇人妻一区二区三区视频| 久久久亚洲精品成人影院| av在线蜜桃| 老司机影院毛片| 白带黄色成豆腐渣| 人妻制服诱惑在线中文字幕| 国产亚洲精品av在线| 日韩欧美一区视频在线观看 | 天美传媒精品一区二区| 好男人视频免费观看在线| 精品欧美国产一区二区三| 久久久久精品性色| 嘟嘟电影网在线观看| 国产亚洲一区二区精品| 国产淫片久久久久久久久| 午夜福利视频精品| 国产av国产精品国产| 午夜福利在线观看吧| 麻豆成人午夜福利视频| 久久久久精品性色| 亚洲精品成人av观看孕妇| 久久草成人影院| 高清在线视频一区二区三区| 欧美潮喷喷水| 久久久久免费精品人妻一区二区| av天堂中文字幕网| 91午夜精品亚洲一区二区三区| 亚洲精品日韩av片在线观看| 国产成人91sexporn| 久久久久久久午夜电影| 日韩av在线大香蕉| 亚洲国产色片| 99九九线精品视频在线观看视频| 国产不卡一卡二| 日本色播在线视频| 老司机影院毛片| 国产伦理片在线播放av一区| 伊人久久精品亚洲午夜| 色尼玛亚洲综合影院| 亚洲无线观看免费| 亚洲人成网站在线观看播放| 欧美高清成人免费视频www| 日韩,欧美,国产一区二区三区| 国产精品一区二区在线观看99 | 国产午夜精品久久久久久一区二区三区| 日韩中字成人| 亚洲国产欧美在线一区| 18禁在线播放成人免费| 69av精品久久久久久| 噜噜噜噜噜久久久久久91| 精品国产三级普通话版| 国产一区二区三区av在线| 久久精品人妻少妇| 国产黄片美女视频| 免费看光身美女| 欧美+日韩+精品| 91狼人影院| 尾随美女入室| 国产精品不卡视频一区二区| kizo精华| 国产伦理片在线播放av一区| 大香蕉久久网| 免费大片黄手机在线观看| 网址你懂的国产日韩在线| 日本免费在线观看一区| 麻豆久久精品国产亚洲av| 老司机影院毛片| 三级国产精品欧美在线观看| 成人毛片a级毛片在线播放| 日韩一本色道免费dvd| 在线观看一区二区三区| 色哟哟·www| 天天躁日日操中文字幕| 日韩 亚洲 欧美在线| 亚洲美女搞黄在线观看| 赤兔流量卡办理| 啦啦啦韩国在线观看视频| av专区在线播放| 久久久精品免费免费高清| 91精品一卡2卡3卡4卡| 肉色欧美久久久久久久蜜桃 | 80岁老熟妇乱子伦牲交| 欧美激情久久久久久爽电影| av在线蜜桃| eeuss影院久久| 97在线视频观看| 久久午夜福利片| 亚洲精品第二区| 国产高清三级在线| 熟妇人妻久久中文字幕3abv| 国产女主播在线喷水免费视频网站 | av.在线天堂| 少妇人妻一区二区三区视频| 国产成人精品久久久久久| 欧美精品一区二区大全| 哪个播放器可以免费观看大片| 麻豆国产97在线/欧美| 日日干狠狠操夜夜爽| 综合色av麻豆| 三级国产精品片| av专区在线播放| 亚洲av二区三区四区| 深夜a级毛片| 在线观看免费高清a一片| 日本av手机在线免费观看| 99久久精品国产国产毛片| 亚洲一区高清亚洲精品| 日日摸夜夜添夜夜添av毛片| 日韩精品青青久久久久久| 久久精品综合一区二区三区| 一级黄片播放器| 亚洲精华国产精华液的使用体验| 亚洲乱码一区二区免费版| 国产成人精品婷婷| 免费观看无遮挡的男女| 又爽又黄无遮挡网站| 秋霞伦理黄片| 国产高清不卡午夜福利| 91aial.com中文字幕在线观看| 一边亲一边摸免费视频| 欧美成人一区二区免费高清观看| 精品亚洲乱码少妇综合久久| 亚洲精华国产精华液的使用体验| 天天躁夜夜躁狠狠久久av| 久久久久久伊人网av| 国产黄色小视频在线观看| 国产黄色视频一区二区在线观看| 真实男女啪啪啪动态图| 国产伦一二天堂av在线观看| 国产熟女欧美一区二区| 两个人的视频大全免费| 男人舔奶头视频| 亚洲熟女精品中文字幕| 97超视频在线观看视频| 亚洲av一区综合| 性插视频无遮挡在线免费观看| 男插女下体视频免费在线播放| 欧美xxxx性猛交bbbb| 亚洲国产成人一精品久久久| 观看免费一级毛片| 成人午夜精彩视频在线观看| 国产探花极品一区二区| 国产熟女欧美一区二区| 真实男女啪啪啪动态图| 中文字幕av在线有码专区| 九草在线视频观看| 在线天堂最新版资源| 男人舔奶头视频| 亚洲欧洲国产日韩| 亚洲18禁久久av| 偷拍熟女少妇极品色| 久久午夜福利片| 国产精品美女特级片免费视频播放器| 超碰97精品在线观看| 午夜久久久久精精品| 国产伦精品一区二区三区四那| 精品酒店卫生间| 久久这里只有精品中国| 卡戴珊不雅视频在线播放| 女的被弄到高潮叫床怎么办| 看十八女毛片水多多多| 亚洲精品成人久久久久久| 美女xxoo啪啪120秒动态图| 欧美丝袜亚洲另类| 青春草视频在线免费观看| 日本黄大片高清| 国产精品av视频在线免费观看| 视频中文字幕在线观看| 国产黄色视频一区二区在线观看| 成年女人在线观看亚洲视频 | 小蜜桃在线观看免费完整版高清| 国产女主播在线喷水免费视频网站 | 女人十人毛片免费观看3o分钟| 99热全是精品| 亚洲av成人精品一区久久| 夫妻性生交免费视频一级片| 婷婷六月久久综合丁香| 亚洲高清免费不卡视频| 日本三级黄在线观看| 网址你懂的国产日韩在线| 亚洲欧美成人综合另类久久久| 欧美zozozo另类| 又爽又黄无遮挡网站| 99re6热这里在线精品视频| 菩萨蛮人人尽说江南好唐韦庄| 免费观看无遮挡的男女| 日韩一区二区三区影片| 午夜精品一区二区三区免费看| 成人高潮视频无遮挡免费网站| 免费观看a级毛片全部| 国产精品一区二区三区四区久久| 嫩草影院新地址| 国产精品美女特级片免费视频播放器| 欧美97在线视频| 欧美xxⅹ黑人| 国产成年人精品一区二区| 性插视频无遮挡在线免费观看| 亚洲精品中文字幕在线视频 | 免费观看的影片在线观看| 亚洲一区高清亚洲精品| 免费观看精品视频网站| 黄色配什么色好看| 成人毛片60女人毛片免费| 国产精品无大码| 国产精品久久久久久精品电影| av又黄又爽大尺度在线免费看| 国产精品女同一区二区软件| 啦啦啦中文免费视频观看日本| 国产精品嫩草影院av在线观看| 欧美成人精品欧美一级黄| 人人妻人人看人人澡| 精品熟女少妇av免费看| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 日本黄色片子视频| 搞女人的毛片| 黄片wwwwww| 欧美成人a在线观看| 午夜久久久久精精品| 国产熟女欧美一区二区| 国产成年人精品一区二区| 亚洲国产色片| 美女被艹到高潮喷水动态| 十八禁网站网址无遮挡 | 日韩欧美一区视频在线观看 | 三级经典国产精品| 2021天堂中文幕一二区在线观| 久久久国产一区二区| 水蜜桃什么品种好| 午夜精品一区二区三区免费看| 国产精品.久久久| av免费观看日本| 日韩av在线大香蕉| 午夜福利成人在线免费观看| 看十八女毛片水多多多| 久久国产乱子免费精品| 一个人免费在线观看电影| 97精品久久久久久久久久精品| 免费人成在线观看视频色| 国产免费一级a男人的天堂| 麻豆国产97在线/欧美| 日本欧美国产在线视频| 成人亚洲精品av一区二区| 久久精品熟女亚洲av麻豆精品 | 亚洲精品自拍成人| 亚洲18禁久久av| 一级毛片我不卡| 99久久人妻综合| 亚洲成色77777| 深夜a级毛片| 日日干狠狠操夜夜爽| 亚洲图色成人| 国产精品一区二区性色av| 国产精品日韩av在线免费观看| 久久久久精品性色| 亚洲图色成人| 精品国产露脸久久av麻豆 | 久久99热这里只频精品6学生| 亚洲不卡免费看| 精品少妇黑人巨大在线播放| 22中文网久久字幕| 国产精品三级大全| 日本黄色片子视频| 日韩,欧美,国产一区二区三区| 久久这里只有精品中国| 亚洲精品国产av成人精品| 老女人水多毛片| 亚洲国产高清在线一区二区三| 国产一区二区在线观看日韩| 免费看美女性在线毛片视频| 久久精品国产亚洲网站| 80岁老熟妇乱子伦牲交| av在线观看视频网站免费| 日韩三级伦理在线观看| 成年人午夜在线观看视频 | 国产探花极品一区二区| 免费播放大片免费观看视频在线观看| 国产精品久久久久久av不卡| 欧美成人午夜免费资源| 国产精品1区2区在线观看.| 一级毛片 在线播放| av女优亚洲男人天堂| 国产亚洲精品av在线| 亚洲乱码一区二区免费版| 免费看美女性在线毛片视频| 色视频www国产| 久久久久久久久久成人| 少妇人妻一区二区三区视频| 日本免费a在线| 免费看av在线观看网站| 一个人看视频在线观看www免费| 亚洲熟女精品中文字幕| 日日干狠狠操夜夜爽| 99热这里只有精品一区| 人体艺术视频欧美日本| 久久久久久久久久黄片| 久久久久网色| 欧美激情国产日韩精品一区| 亚洲伊人久久精品综合| 亚洲欧美中文字幕日韩二区| 久久久久久久久久成人| 亚洲综合色惰| 在线免费十八禁| 内射极品少妇av片p| 菩萨蛮人人尽说江南好唐韦庄| 国内精品宾馆在线| 可以在线观看毛片的网站| 边亲边吃奶的免费视频| 久久久久久久午夜电影| 女人被狂操c到高潮| 人人妻人人澡人人爽人人夜夜 | 毛片一级片免费看久久久久| 日本一本二区三区精品| 国产爱豆传媒在线观看| or卡值多少钱| 中文字幕制服av| 你懂的网址亚洲精品在线观看| 丰满人妻一区二区三区视频av| 国产亚洲精品久久久com| 欧美潮喷喷水| 久久国产乱子免费精品| 日韩欧美精品免费久久| 少妇的逼水好多| 日韩精品有码人妻一区| 好男人在线观看高清免费视频| 国产又色又爽无遮挡免| 免费观看a级毛片全部| 一区二区三区乱码不卡18| 午夜免费男女啪啪视频观看| 久久久成人免费电影| 一级黄片播放器| 全区人妻精品视频| 嫩草影院新地址| 夜夜爽夜夜爽视频| 少妇人妻一区二区三区视频| 日本一二三区视频观看| 国产在线男女| av免费观看日本| 丰满少妇做爰视频| 成年女人在线观看亚洲视频 | 肉色欧美久久久久久久蜜桃 | 日韩欧美三级三区| 熟女人妻精品中文字幕| 久久久欧美国产精品| 免费观看a级毛片全部| 精品久久久久久久末码| 亚洲三级黄色毛片| av在线观看视频网站免费| 91久久精品国产一区二区成人| 日日啪夜夜爽| 22中文网久久字幕| 天美传媒精品一区二区| 精品久久久久久久久久久久久| 纵有疾风起免费观看全集完整版 | 国内精品宾馆在线| 亚洲欧美成人精品一区二区| 国产成人免费观看mmmm| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜 | 日韩亚洲欧美综合| 九草在线视频观看| 一级毛片电影观看| 特大巨黑吊av在线直播| 欧美成人精品欧美一级黄| 国产精品1区2区在线观看.| 亚洲内射少妇av| 内地一区二区视频在线| 日韩av不卡免费在线播放| 在线观看免费高清a一片| 赤兔流量卡办理| 亚洲综合色惰| 成人鲁丝片一二三区免费| 亚洲欧美精品自产自拍| 亚洲av一区综合| 丝袜喷水一区| 成人无遮挡网站| 看黄色毛片网站| 色综合亚洲欧美另类图片| 色视频www国产| 欧美精品一区二区大全| 丰满少妇做爰视频| 国产午夜精品论理片| 丰满人妻一区二区三区视频av| 久久人人爽人人片av| 免费av不卡在线播放| 成人特级av手机在线观看| 亚洲不卡免费看| 亚洲精品视频女| 大片免费播放器 马上看| 成人午夜精彩视频在线观看| 又大又黄又爽视频免费| 中文乱码字字幕精品一区二区三区 | 亚洲性久久影院| 久久人人爽人人爽人人片va| 国产91av在线免费观看| 亚洲激情五月婷婷啪啪| 99久久精品国产国产毛片| 成年版毛片免费区| 超碰av人人做人人爽久久| 国产欧美日韩精品一区二区| 美女国产视频在线观看| 午夜精品国产一区二区电影 | 国产午夜精品一二区理论片| 丰满人妻一区二区三区视频av| 国产91av在线免费观看| 亚洲精品,欧美精品| 久久久久网色| 波野结衣二区三区在线| 美女xxoo啪啪120秒动态图| 日韩电影二区| 日韩欧美三级三区| 久久精品国产亚洲av天美| 日本色播在线视频| 国产视频内射| av网站免费在线观看视频 | 超碰97精品在线观看| 午夜福利网站1000一区二区三区| 国产精品久久久久久av不卡| 91狼人影院| 精品久久久久久久久av| 午夜久久久久精精品| 欧美日韩在线观看h| 日韩 亚洲 欧美在线| 午夜福利在线观看免费完整高清在| 国产高清有码在线观看视频| 色5月婷婷丁香| 男女国产视频网站| 成人综合一区亚洲| 免费不卡的大黄色大毛片视频在线观看 | 免费黄网站久久成人精品| 亚洲乱码一区二区免费版| 亚洲成人久久爱视频| av在线蜜桃| 午夜精品在线福利| 国产精品久久久久久久电影| 美女内射精品一级片tv| 99久久九九国产精品国产免费| 亚洲精品久久午夜乱码| 国产又色又爽无遮挡免| 欧美精品国产亚洲| 日韩欧美国产在线观看| 欧美日韩精品成人综合77777| 97精品久久久久久久久久精品| 日本免费在线观看一区| 欧美zozozo另类| 日日摸夜夜添夜夜添av毛片| 2018国产大陆天天弄谢| 午夜久久久久精精品| 欧美日韩精品成人综合77777| 免费黄色在线免费观看| 身体一侧抽搐| 18禁裸乳无遮挡免费网站照片| 九草在线视频观看| 国产精品女同一区二区软件| 青春草视频在线免费观看| 国产精品一及| 91精品国产九色| 婷婷色麻豆天堂久久| 国产片特级美女逼逼视频| 精品国产一区二区三区久久久樱花 | 国产av不卡久久| 高清在线视频一区二区三区| 黄片wwwwww| 国产v大片淫在线免费观看| 在线观看一区二区三区| 久热久热在线精品观看| 久久精品久久久久久噜噜老黄| 极品教师在线视频| 男人舔奶头视频| 精品人妻熟女av久视频| 亚洲欧美一区二区三区黑人 | 成年女人看的毛片在线观看| 午夜福利在线观看吧| 精品酒店卫生间| 91在线精品国自产拍蜜月| 亚洲性久久影院| av在线天堂中文字幕| 18+在线观看网站| 一级毛片电影观看| 国产伦精品一区二区三区四那| 国产精品一区二区三区四区久久| 老师上课跳d突然被开到最大视频| 亚洲美女视频黄频| 精品久久久久久成人av| 看非洲黑人一级黄片| 欧美 日韩 精品 国产| 久久国产乱子免费精品| 国产精品蜜桃在线观看| 日韩精品有码人妻一区| av线在线观看网站| 日韩制服骚丝袜av| 国产乱人偷精品视频| 久久久久久久久大av| 国产大屁股一区二区在线视频| 亚洲av.av天堂| 国产精品麻豆人妻色哟哟久久 | 精品久久久久久成人av| 九九久久精品国产亚洲av麻豆| 亚洲精品一区蜜桃| 一级毛片 在线播放| 国产日韩欧美在线精品| 日韩av免费高清视频| 少妇的逼水好多| 一级毛片久久久久久久久女| 久久久久久久久久成人| 一夜夜www| 精品酒店卫生间| 精品久久久精品久久久| 肉色欧美久久久久久久蜜桃 | 高清欧美精品videossex| 国产精品一区二区性色av| 国产成人aa在线观看| av在线观看视频网站免费| 亚洲成人一二三区av| 久久精品久久久久久噜噜老黄| 欧美成人午夜免费资源| 99久久精品热视频| 国产大屁股一区二区在线视频| 亚洲国产欧美在线一区| av专区在线播放| 中国国产av一级| 搞女人的毛片|