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

    不同振動特性參數(shù)對杏樹振動響應的影響

    2019-02-25 01:51:32楊會民散鋆龍陳毅飛王學農(nóng)牛長河侯書林
    農(nóng)業(yè)工程學報 2019年2期
    關鍵詞:有限元振動分析

    楊會民,散鋆龍,陳毅飛,王學農(nóng),牛長河,侯書林

    ?

    不同振動特性參數(shù)對杏樹振動響應的影響

    楊會民1,散鋆龍2,陳毅飛1,王學農(nóng)1※,牛長河1,侯書林2

    (1. 新疆農(nóng)業(yè)科學院農(nóng)業(yè)機械化研究所,烏魯木齊 830091;2. 中國農(nóng)業(yè)大學工學院,北京 100083)

    為研究振動時間、振動頻率和振動激勵點振幅等不同振動特性參數(shù)對杏樹振動的影響,該文利用ANSYS軟件對杏樹進行了有限元建模分析;通過三因素三水平試驗分析不同振動特性參數(shù)對杏樹振動檢測點的影響,利用Design Expert軟件進行優(yōu)化分析,并進行實驗室驗證試驗。杏樹自由模態(tài)振動響應分析表明最佳杏樹振動采收響應頻率范圍為0~20 Hz;諧振動響應分析可知在最佳頻率范圍內,杏樹振動激勵點振幅為5、10和15 mm時,同一頻率下,隨著激勵振幅的增大,相同位置加速度增大,但振動曲線整體變化規(guī)律和趨勢一致。試驗分析可知,各因素影響檢測點1和2加速度的強弱順序一致:激勵點振幅>振動頻率>振動時間;各因素影響檢測點3加速度的強弱順序為:振動時間>激勵點振幅>振動頻率;建立3個檢測點的響應方程,由下至上3個檢測點的回歸方程決定系數(shù)分別為0.906 7、0.879 3和0.973 3;多目標參數(shù)優(yōu)化結果為:振動時間7.207 s,振動頻率15 Hz,激振點振幅10 mm,通過驗證試驗可知由下至上各檢測點加速度為10.4、10.2和9.3,與優(yōu)化值相近。該研究可為杏振動采收機械參數(shù)設計提供參考。

    振動;收獲;機械化;杏;振動響應試驗

    0 引 言

    杏(Lam.)屬薔薇科()杏屬(Mill.)落葉喬木,是世界重要的經(jīng)濟果樹之一[1]。中國杏樹栽培歷史悠久,根據(jù)遺傳學研究表明,新疆伊犁河北岸野生普通杏是世界栽培杏的起源種群[2]。截止2016年底,新疆維吾爾自治區(qū)杏樹種植面積12.487萬hm2,產(chǎn)量115.42萬t[3]。近年來新疆杏產(chǎn)業(yè)化發(fā)展迅速,但其收獲方式仍停留在利用長梯或長桿輔助的傳統(tǒng)人工采摘,作業(yè)風險較高,同時增大了收獲成本,收獲環(huán)節(jié)成為制約新疆林果產(chǎn)業(yè)發(fā)展的首要問題。振動采收是目前林果采收的主要手段之一,利用機械化采收方式可有效降低勞動成本,提高林果采收效率。

    國外已形成了一套完整的林果振動采收體系,林果振動采收理論完善,機具成熟,其中不同振動采收參數(shù)對林果收獲的影響是主要的研究方向之一[4-11]。由于果樹自身生長特性與材料特性,不同品種果樹的振動響應不同,采收機械振動頻率、振幅、振動時間以及振動激勵點位置等振動特性參數(shù)的選擇和設定將直接影響林果采收效率與質量[11-12]。早在1957年美國就已展開往復式振動采收機工作頻率和振幅對李子采收、果實損傷以及果樹搖動所需動力的研究[13-14]。Aristizábal等[15]研究使用圓周和多方向振動采收咖啡過程中,振動頻率、振幅、振動時間以及振動激勵位置對成熟果實采收率的影響。Polat等[16]為獲得巴丹木最佳振動采收參數(shù),利用慣性式懸臂振動機進行振動采收,研究相同振動時間下,不同振動頻率和振幅參數(shù)組合對果實脫落率的影響。

    國內開展了大量林果振動采收機構設計、振動傳遞特性以及果樹響應狀態(tài)等方面的基礎研究,但對振動特性參數(shù)對林果采收影響研究較淺[17-26]。范雷剛等[27]研究了海棠果樹受不同振動特性參數(shù)組合時,不同位置加速度變化。杜小強等[28]對所設計的三維激振林果采收機構進行優(yōu)化設計,并對系統(tǒng)振動響應進行了分析。王東等[29]利用有限元進行果樹響應分析,研究不同樹形果樹振動響應狀態(tài)。

    為研究不同振動采收參數(shù)對林果機械采收的影響,指導杏林果振動采收機械設計,優(yōu)化杏振動采收機械性能,本文利用ANSYS建立杏樹有限元模型,進行杏樹自由模態(tài)與諧響應仿真分析,確定杏振動采收參數(shù)范圍。在實驗室內進行振動試驗,分析振動試驗結果,研究不同振動特性參數(shù)組合對振動采收的影響,得到杏最佳振動參數(shù)組合,并進行實驗室驗證試驗,以期為杏振動采收機械參數(shù)設計提供設計及理論依據(jù)。

    1 杏振動采收機理

    根據(jù)4GZG-25型牽引式偏心式林果采收樣機工作方式與振動特性,建立偏心塊對稱放置的振動采收機-果樹模型,如圖1所示[30-31]。在振動采收機與杏樹加持位置水平面上,以杏樹夾持位置中心為坐標系原點,方向為系統(tǒng)水平運動方向,方向為與水平垂直運動方向。杏樹動力學特性由加持位置杏樹樹干等效彈性系數(shù)(常數(shù))與等效阻尼系數(shù)(常數(shù))表示。

    注:m為偏心塊質量,kg;ω為角速度,rad·s-1;t為時間,s;L1為振動箱與杏樹間距離,m; L2為偏心塊旋轉軸中心與x軸間距離,m;r為偏心塊轉動中心與質心的距離,mm;k為杏樹樹干等效彈性系數(shù),常數(shù);c為杏樹樹干等效阻尼系數(shù),常數(shù)。

    對稱放置的2個偏心塊質量相等,轉動方向相反,結合杏樹-振動機動力學模型,由牛頓第三定律可知,由于偏心塊做圓周運動,采收機振動過程中偏心機構整體在方向產(chǎn)生的力相互疊加,在方向上的力相互抵消,即方向受力為零。由此可的偏心塊產(chǎn)生的慣性力為

    式中F為偏心機構方向上產(chǎn)生的力,N;F為偏心機構方向上產(chǎn)生的力,N;為偏心機構在方向上的位移,mm。

    對系統(tǒng)整體進行受力分析,并簡化可得杏樹-振動采收機振動微分方程為

    2 杏樹振動響應試驗材料與方法

    由于所設計的4GZG-25型牽引式偏心式林果采收樣機正弦振動激勵輸出特性,為優(yōu)化樣機工作參數(shù)提供參考,研究杏樹受迫振動后的響應狀態(tài),獲得最佳的振動響應參數(shù);通過有限元建模與實驗室振動響應試驗相結合對杏樹振動響應進行研究,對優(yōu)化后的參數(shù)進行實驗室驗證試驗。杏樹振動響應試驗路線如圖2所示。

    圖2 杏樹振動響應試驗路線

    根據(jù)牛頓第二定律,振動過程中杏果實在質量一定的情況下,加速度越大,果實由于慣性產(chǎn)生的脫落力越大。因此本文以杏樹檢測點處加速度為主要評價指標進行試驗,加速度單位為(為重力加速度,取9.8 m/s2)。

    2.1 試驗材料

    以新疆主要杏品種庫麥提杏樹為主要研究對象,樹形結構為自然開心型[29]。2018年6月底杏成熟期,將新疆農(nóng)業(yè)科學院輪臺國家果樹資源圃內樹齡為3 a的果樹,移栽到新疆農(nóng)業(yè)科學院農(nóng)業(yè)機械化研究所物料特性實驗室內,進行測繪與基本參數(shù)測量。

    2.2 試驗儀器與軟件

    試驗儀器:蘇州試驗儀器設備公司DC-300振動試驗臺,正弦推力2.94 kN,工作頻率5~5 000 Hz,最大位移40 mm;SA-3振動平臺功率放大器,額定輸出功率15 kVA;YMC壓電式IEPE型加速度傳感器,靈敏度50 mV/ms2,測試范圍±10,頻率范圍0~4 000 Hz。美國NI公司NI DAQ 9 178四通道信號采集器,NI 9 234振動信號采集卡。上海一恒科學儀器有限公司BPG-9 140A型精密鼓風干燥箱,功率2 050 W,溫度范圍RT+10~25 ℃。美特斯工業(yè)系統(tǒng)(中國)有限公司CMT 6 103系列電子萬能試驗機,最大應力1 kN。

    試驗軟件:ANSYS Workbench 17.0有限元分析軟件;SolidWorks 2017三維建模軟件;Labview振動測試虛擬軟件;Design expert 10.0.3試驗分析軟件。

    2.3 試驗方法

    利用SolidWorks 2017軟件對已測繪的杏樹進行三維建模,使用ANSYS Workbench 17.0軟件對杏樹進行有限元分析,有限元模型如圖3所示。分析杏樹不附加任何振動激勵狀態(tài)下的自由模態(tài),結合模態(tài)云圖仿真結果中杏樹響應的一致性,獲得最佳的振動頻率范圍;研究杏樹不同振幅下受正弦振動激勵時沿樹干3個檢測位置的諧響應狀態(tài),每個檢測位置間距離為400 mm,以加速度大小表示檢測位置的響應狀態(tài)。

    圖3 杏樹有限元模型

    根據(jù)有限元分析相應結果,利用DC-300振動實驗平臺對杏樹施加正弦振動激勵,研究振動時間、振動頻率和加持點振幅等振動采收主要因素對杏樹不同位置加速度的影響。杏樹從果園中砍伐后立即移入實驗室,如圖4所示。去掉細小樹枝與樹葉,利用自制的連接器將杏樹與振動平臺連接,杏樹根部固定。沿杏樹樹干安裝加速度傳感器,傳感器安裝位置與有限元分析時果樹3個檢測點位置相同,通過Labview振動檢測軟件得到NI DAQ 9 178振動信號采集器所采集的加速度數(shù)據(jù)。采用Design Expert 10.0.3 軟件進行試驗方案設計,分析杏樹不同位置在不同振動激勵參數(shù)下的響應狀態(tài),得到最佳的振動激勵參數(shù),并進行實驗室驗證試驗。

    圖4 杏樹振動檢測平臺

    3 杏樹有限元建模與分析

    根據(jù)國家標準GB/T 1 931-2009木材含水率測定方法和GB/T 1 936.2-2009木材抗彎彈性模量測定方法測得庫麥提杏樹木材含水率為41.01%,彈性模量為425.35 MPa,密度為0.856 g/cm3,泊松比取0.330,由此可得杏樹木材剪切模量為159.9 MPa。將所測得的杏樹木材基本參數(shù)代入ANSYS Workbench建立的杏樹模型中,對杏樹進行有限元分析。

    3.1 杏樹自由振動模態(tài)分析

    利用ANSYS Workbench中Modal模塊對杏樹進行1~50階模態(tài)分析,得到反應杏樹模型自由模態(tài)下階次、響應頻率和杏樹整體形變關系的3D曲面圖,如圖5所示。

    圖5 杏樹自由模態(tài)響應分析

    由圖5可知,杏樹在15、26、30和40階時具有典型自由模態(tài)響應,對應的頻率分別為5.5、10.5、15.1和29.3 Hz,果樹最大形變分別為286.5、261.9、267.2和273.5 mm。

    分析杏樹典型自由模態(tài)云圖6可知,第15和26階典型自由模態(tài)中杏樹枝干末端響應性較高;第30階果樹整體響應一致性高,杏樹末端形變響應最大;第40階29.3 Hz時果樹部分枝條變形嚴重,樹形整體結構易被破壞。由此可知:杏樹振動采收過程中最佳頻率范圍可控制在5~15 Hz之間。

    圖6 杏樹典型模態(tài)響應云圖

    3.2 杏樹諧振動響應分析

    通過對杏樹模型進行諧響應分析,確定杏樹枝干在固定幅值的正弦位移作用下的穩(wěn)定響應,采用完全法進行求解。根據(jù)果樹振動采收激振點振幅特性[31],夾持點振幅設定為5、10和15 mm進行杏樹諧振動響應分析。由自由模態(tài)仿真結果設定頻率范圍為0~20 Hz,進行諧振動分析可得到檢測點1、2、3位置加速度值變化曲線,如圖7。

    圖7 不同振幅激勵下檢測點X軸方向加速度變化曲線

    由圖7可知,激勵點振幅為5、10和15 mm時,在0~20 Hz振動頻率區(qū)間內,檢測點1在4.67、7.33、15.33和18.66 Hz時具有典型加速度峰值;檢測點2在6.34和17.33 Hz時具有典型加速度峰值;檢測點3在7.01和17.3 Hz時具有典型加速度峰值。檢測點位置由下至上,振幅為5 mm時,最大加速度分別為7.52、6.57和4.93;振幅為10 mm時,最大加速度分別為14.89、12.96和9.87;振幅為15 mm時,最大加速度分別為22.36、19.71和14.8。

    由此可知,在0~20 Hz內,振幅為5、10和15 mm時,同一頻率下,隨著激勵振幅的增大,相同位置加速度增大,但振動曲線整體變化規(guī)律和趨勢一致。

    4 杏樹實驗室振動試驗

    為了驗證有限元分析結果,獲得最優(yōu)振動采收參數(shù)組合,指導振動采收樣機參數(shù)設計,根據(jù)樣機工作原理與ANSYS有限元仿真結果,在2018年6月底在新疆農(nóng)業(yè)科學院農(nóng)業(yè)機械化研究所物料特性實驗室中,按照試驗路線圖2,利用振動試驗平臺對剛移栽入實驗室的成熟期杏樹施加振動激勵,進行振動試驗,試驗數(shù)據(jù)由加速度傳感器與NI 9 234振動信號采集卡采集。

    4.1 評價指標

    選擇振動時間(s)、激勵點振幅(mm)和振動頻率(Hz)3個試驗因素,在實驗室進行振動試驗,每組試驗重復3次。樹枝加速度是果實振動脫落的主要評價指標,本文選擇3個主要檢測點所測得的加速度作為評價指標。

    4.2 振動試驗設計

    進行三因素三水平組合試驗,根據(jù)杏樹有限元自由模態(tài)響應與諧振動響應分析結果,試驗因素與水平編碼如表1所示。

    利用Design Expert 10.0.3 軟件進行試驗方案設計,共進行17組試驗,每組重復3次,試驗結果取平均值,試驗數(shù)據(jù)如表2。

    表1 試驗因素和水平編碼

    表2 試驗設計方案及響應值結果

    注:1、2、3為圖3中3個檢測點的加速度。

    Note:1,2, and3are the accelerations of the three detection points in Fig.3.

    4.3 試驗結果分析與回歸模型建立

    通過Design Expert 10.0.3軟件對試驗結果進行回歸分析,并進行多元回歸擬合,得到3個檢測點加速度1、2和3,3個試驗指標的回歸方程,并驗證其顯著性,試驗分析結果如表3所示。

    4.3.1 檢測點1加速度的顯著性分析

    由表3可知,整體模型極其顯著(<0.01),各因素影響檢測點1加速度的強弱順序為;其中模型回歸項和顯著,因素不顯著(>0.05)可忽略,擬合后得到回歸方程決定系數(shù)2為0.906 7,擬合程度高,方程為

    式中為激勵點振幅,mm;為振動頻率,Hz。

    4.3.2 檢測點2加速度的顯著性分析

    由表3可知,整體模型極其顯著(<0.01),各因素影響檢測點2加速度的強弱順序為>>,與檢測點1一致,其中模型回歸項和顯著,因素不顯著(>0.05)可忽略,擬合后得到回歸方程決定系數(shù)2為0.879 3,方程為

    4.3.3檢測點3加速度的顯著性分析

    由表3可知,整體模型極其顯著(<0.01),整體模型極其顯著,各因素影響檢測點3加速度的強弱順序為>>。其中除模型回歸項、與2的檢測點3加速度值大于0.05不顯著可忽略,其他回歸項均為顯著。由于檢測點3位于整個果樹末端,檢測點位置、振動過程中位移、果樹樹枝粗細以及果果枝形狀等原因會造成因素對不同檢測點的影響強弱順序不同。回歸方程決定系數(shù)2為0.973 3,擬合程度極高,方程為

    式中為振動時間,s。

    表3 檢測點加速度回歸模型方差F檢驗與分析

    注:<0.01,為極其顯著;0.01<<0.05為顯著。

    Note:<0.01, extremely significant; 0.01<<0.05, significant.

    4.4 各項振動特性參數(shù)優(yōu)化

    振動過程果實脫落的響應加速度越大,果實所受脫落力越大,果實更易脫落,即3個檢測點加速度響應均為最大值時滿足優(yōu)化條件。根據(jù)上述有限元分析結果、振動試驗與項目組前期的試驗結果[30-31],為避免振動過程造成果樹損傷,3個檢測點加速度優(yōu)化范圍為0~20;根據(jù)檢測點由下至上所處的位置,設定3個檢測點值的權重,檢測點1、2設置權重為1,檢測點3接近果實生長位置權重為2。

    根據(jù)上述參數(shù)設定,由Design Expert10.0.3軟件中Optimization模塊(響應優(yōu)化)對檢測點進行多目標優(yōu)化。

    目標編碼變量優(yōu)化函數(shù)為

    對目標函數(shù)進行優(yōu)化得到最優(yōu)振動采收參數(shù)組合為:振動時間7.207 s,振動頻率15 Hz,激振點振幅10 mm。

    4.5 實驗室驗證試驗

    為驗證最終優(yōu)化參數(shù),按優(yōu)化后得到參數(shù)調整振動實驗平臺設置參數(shù),在實驗室中進行振動試驗,試驗重復3次取平均值,優(yōu)化值與實際測量值之間分析驗證結果如表4所示,實驗值與優(yōu)化參數(shù)相近,偏差較小。

    表4 優(yōu)化值與實際值對比結果

    5 結 論

    1)對建立的杏樹模型進行ANSYS有限元分析,通過杏樹自由模態(tài)振動響應分析可知最佳杏樹振動采收響應頻率范圍為0~20 Hz;通過諧振動響應分析可知在0~20 Hz內,杏樹振動激勵點振幅為5、10和15 mm時,同一頻率下,隨著激勵振幅的增大,相同位置加速度增大,但振動曲線整體變化規(guī)律和趨勢一致。

    2)研究杏樹振動特性參數(shù)(振動時間、激勵點振幅和振動頻率)對不同檢測點振動加速度的影響規(guī)律,進行三因素三水平試驗分析,各參數(shù)對檢測點1和2的影響強弱順序一致為>>;對檢測點3的影響強弱順序為>>。

    3)建立3個檢測點響應方程,由下至上3個檢測點的回歸方程決定系數(shù)分別為0.906 7、0.879 3和0.973 3,為杏樹振動特性分析提供參考。

    4)利用Design Expert10.0.3軟件響應優(yōu)化模塊對檢測點響應方程進行優(yōu)化,獲得最優(yōu)振動采收參數(shù)組合為:振動時間7.207 s,振動頻率15 Hz,激振點振幅10 mm,經(jīng)試驗驗證可知由下至上3個檢測點加速度分別為10.4、10.2和9.3,試驗值與優(yōu)化參數(shù)相符,該研究可為杏振動采收機械參數(shù)設計提供設計及理論依據(jù)。

    [1] 包文泉,烏云塔娜,王淋,等. 野生杏和栽培杏的遺傳多樣性和遺傳結構分析[J]. 植物遺傳資源學報,2017,18(2):201-209.

    Bao Wenquan, Wuyun Tana, Wang Lin, et al. Genetic diversity and population structure of the wild apricot and cultivation apricot[J]. Journal of Plant Genetic Resources, 2017, 18(2): 201-209. (in Chinese with English abstract)

    [2] 廖康. 新疆野生果樹資源研究[M]. 烏魯木齊:新疆人民出版社,2013.

    [3] 中華人民共和國國家統(tǒng)計局. 中國統(tǒng)計年鑒[M]. 北京:中國統(tǒng)計出版社,2017.

    [4] BlancoRoldán G L, GilRibes J A, Kouraba K, et al. Effects of trunk shaker duration and repetitions on removal efficiency for the harvesting of oil olives[J]. Applied Engineering in Agriculture, 2009, 25(3): 329-334.

    [5] Pezzi F, Caprara C. Mechanical grape harvesting: Investigation of the transmission of vibrations[J]. Biosystems Engineering, 2009, 103(3): 281-286.

    [6] Torregrosa A, Ortí E, Martín B, et al. Mechanical harvesting of oranges and mandarins in Spain[J]. Biosystems Engineering, 2009, 104(1): 18-24.

    [7] Whitney J D, Smerage G H, Block W A. Dynamic analysis of a trunk shaker-post system[J]. Transactions of the ASAE, 1990, 33(4): 1066-1070.

    [8] Horvath E, Sitkei G. Damping properties of plum trees shaken at their trunks[J]. Transactions of the ASAE, 2005, 48(1): 19-25.

    [9] Láng Z. A one degree of freedom damped fruit tree model[J]. Agricultural Engineering International Cigr Journal, 2008, 51(3): 823-829.

    [10] Horvath E, Sitkei G. Energy consumption of selected tree shakers under different operational conditions[J]. Journal of Agricultural Engineering Research, 2001, 80(2): 191-199.

    [11] Torregrosa A, Martin B, Garcia Brunton J, et al. Mechanical harvesting of processed peaches[J]. Applied Engineering in Agriculture, 2008, 24(6): 723-729.

    [12] Castro-García S, Blanco-Roldán G L, Gil-Ribes J A. Vibrational and operational parameters in mechanical cone harvesting of stone pine (Pinus pinea L.)[J]. Biosystems Engineering, 2012, 112(4): 352-358.

    [13] Castro-García S, Blanco-Roldán G L, Gil-Ribes J A, et al. Dynamic analysis of olive trees in intensive orchards under forced vibration[J]. Trees, 2008, 22(6): 795-802.

    [14] Adrian P A, Fridley R B. Mechanical fruit tree shaking[J]. California Agriculture, 1958, 12: 3-15.

    [15] Aristizábal I D, Oliveros C E, Alvarez F. Mechanical harvest of coffee applying circular and multidirectional vibrations.[J]. Transactions of the ASAE, 2003, 46(2): 205-210.

    [16] Polat R, Guner M, Dursun E, et al. Mechanical harvesting of almond with an inertia type limb shaker[J]. Asian Journal of Plant Sciences, 2007, 6(3): 528-532.

    [17] Du X, Chen D, Zhang Q, et al. Dynamic responses of sweet cherry trees under vibratory excitations[J]. Biosystems Engineering, 2012, 111(3): 305-314.

    [18] Wu C, He L, Du X, et al. 3D reconstruction of Chinese hickory tree for dynamics analysis[J]. Biosystems Engineering, 2014, 119(1): 69-79.

    [19] Du X Q, Wu C Y, He L Y, et al. Dynamic characteristics of dwarf Chinese hickory trees under impact excitations for mechanical fruit harvesting[J]. International Journal of Agricultural & Biological Engineering, 2015, 8(1):17-25.

    [20] 杜小強,倪柯楠,潘珂,等. 可調振幅單向拽振式林果采收機構參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學報,2014,30(16):25-32.

    Du Xiaoqiang, Ni Ke’nan, Pan Ke, et al. Parameter optimization of stroke-adjustable and monodirectional pulling fruit harvester[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE) , 2014, 30(16): 25-32. (in Chinese with English abstract)

    [21] 杜小強,倪柯楠,武傳宇. 基于外旋輪線軌跡的果品振動采收機構研究[J]. 農(nóng)業(yè)機械學報,2016,47(3):59-66.

    Du Xiaoqiang, Ni Ke’nan, Wu Chuanyu. Vibratory harvesting mechanism for tree fruit based on epitrochoid[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3): 59-66. (in Chinese with English abstract).

    [22] 杜小強,李黨偉,賀磊盈,等. 基于電子果實技術的機械振動采收過程果實運動分析[J]. 農(nóng)業(yè)工程學報,2017,33(17):58-64.

    Du Xiaoqiang, Li Dangwei, He Leiying, et al. Fruit motion analysis in process of mechanical vibration harvesting based on electronic fruit technique[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(17): 58-64. (in Chinese with English abstract)

    [23] 喬園園,牛長河,孟詳金,等. 牽引式林果振動采收機的設計與田間試驗[J]. 新疆農(nóng)業(yè)科學,2015,52(3):528-534.

    [24] 王長勤,許林云,周宏平,等. 偏心式林果振動采收機的研制與試驗[J]. 農(nóng)業(yè)工程學報,2012,28(16):10-16.

    Wang Changqin, Xu Linyun, Zhou Hongping, et al. Development and experiment of eccentric-type vibratory harvester for forest-fruits[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(16): 10-16. (in Chinese with English abstract)

    [25] 蔡菲,王春耀,王學農(nóng),等. 基于高速攝像技術的振動落果慣性力研究[J]. 西北農(nóng)林科技大學學報:自然科學版,2013,41(4):208-212.

    Cai Fei, Wang Chunyao, Wang Xuenong, et al. Inertia force of fruits abscised by vibration based on high-speed video camera technology[J]. Journal of Northwest A&F University: Nat. Sci. Ed, 2013, 41(4): 208-212. (in Chinese with English abstract)

    [26] 付威,何榮,曲金麗,等. 自走式矮化密植紅棗收獲機的設計[J]. 農(nóng)機化研究,2014,36(4):106-109.

    [27] 范雷剛,王春耀,劉夢霞,等. 振動參數(shù)對果樹采收影響的試驗研究[J]. 農(nóng)機化研究,2016,38(10):165-168.

    [28] 杜小強,李松濤,賀磊盈,等. 三維激振果品采收機構優(yōu)化設計與試驗[J]. 農(nóng)業(yè)工程學報,2017,33(16):48-55.

    Du Xiaoqiang, Li Songtao, He Leiying, et al. Optimal design and experiment on vibratory fruit harvesting mechanism with three-dimensional excitation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(16): 48-55. (in Chinese with English abstract)

    [29] 王冬,陳度,王書茂,等. 基于有限元方法的整形果樹振動收獲機理分析[J]. 農(nóng)業(yè)工程學報,2017,33(增刊1):56-62.

    Wang Dong, Chen Du, Wang Shumao, et al. Analysis on vibratory harvesting mechanism for trained fruit tree based on finite element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(Supp.1): 56-62. (in Chinese with English abstract)

    [30] 散鋆龍,楊會民,王學農(nóng),等. 振動方式和頻率對杏樹振動采收響應的影響[J]. 農(nóng)業(yè)工程學報,2018,34(8):10-17.

    San Yunlong, Yang Huimin, Wang Xuenong, et al. Effects of vibration mode and frequency on vibration harvesting of apricot trees[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(8): 10-17. (in Chinese with English abstract)

    [31] 散鋆龍,楊會民,王學農(nóng),等. 振動收獲過程中杏果實脫落的動態(tài)響應分析[J]. 農(nóng)業(yè)工程學報,2018,34(18):68-75.

    San Yunlong,Yang Huimin,Wang Xuenong, et al. Dynamic response analysis of apricot fruit dropping during vibration harvesting[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 68-75. (in Chinese with English abstract)

    Influence of different vibration characteristic parameters on vibration response of apricot trees

    Yang Huimin1, San Yunlong2, Chen Yifei1, Wang Xuenong1※, Niu Changhe1, Hou Shulin2

    (1.830091,; 2.100083)

    In this paper, the combination of finite element modeling and experiment was used to study the effects of different vibration parameters such as vibration time (), amplitude () of the vibration excitation point on apricot tree vibration and vibration frequency (), and optimize the performance of forest fruit vibration harvesting prototype. The effects of different vibration parameters on the vibration detection points of apricot trees were analyzed by vibration response test. The response analysis of apricot tree free modal vibration using ANSYS Workbench software showed that apricot trees had typical free modal responses at 15, 26, 30 and 40 orders, the corresponding frequencies were 5.5, 10.5, 15.1, and 29.3 Hz, respectively; and the maximum displacement deformation of fruit trees was 286.5, 261.9, 267.2, and 273.5 mm, respectively. The typical free modal cloud image analysis showed that the typical free modalities of the 15th and 26th orders were higher in the terminal end responsiveness of the apricot branches. The 30th order fruit tree had a high overall response consistency, and the apricot tree had the largest deformation response at the end; at the 40th order of 29.3 Hz, some branches of the fruit trees were severely deformed, and the overall structure of the tree was easily destroyed. The optimal response frequency of apricot tree vibration harvesting ranged from 0 to 20 Hz. The harmonic response analysis showed that in the optimal frequency range, the acceleration of the same position increased with the increase of the excitation amplitude at the same frequency, but the overall variation of the vibration curve was consistent with the trend; When the amplitude of the vibration excitation point of the apricot tree was 5, 10 and 15 mm, at the same frequency, as the excitation amplitude increased, the acceleration at the same position increased, but the overall variation of the vibration curve was consistent with the trend. The three-factor and three-level vibration response tests were conducted to study the effects of vibration time, vibration frequency and excitation point amplitude on the acceleration of three different detection points. The multivariate regression analysis of variance showed that the accelerations1and2of detection points 1 and 2 were less than 0.000 1, the overall model was highly significant (<0.01), and the model regression termsandwere significant. The acceleration3of the detection point 3 was 0.000 1, and the overall model was highly significant (<0.01), except that the acceleration3value of the detection point 3 of the model regression terms,and2was not significant (>0.05), and other regression terms were significant (<0.05). The factors which affected the acceleration of the detection points 1 and 2 were the same as the amplitude of the excitation point, vibration frequency and vibration time. The order of the magnitude of the acceleration affecting the detection point 3 was the vibration time, the amplitude of the excitation point, and the vibration frequency. Through the response equations of the three detection points, the coefficient of determination2of the regression equations from the bottom to the top three detection points were 0.906 7, 0.879 3 and 0.973 3, respectively. Using the Design-Expert 10.0.3 software to optimize the detection point response equation, the optimal vibration recovery parameter combination was that the vibration time was 7.207 s, the vibration frequency was 15 Hz, and the amplitude of the excitation point was 10 mm. The verification test showed that the acceleration from the bottom to the top of each detection point was 10.4, 10.2 and 9.3, which was similar to the optimized value. These conclusions can provide design and theoretical basis for the design of mechanical parameters of apricot vibration harvesting.

    vibrations; harvesting; mechanization; apricot; vibration response test

    10.11975/j.issn.1002-6819.2019.02.002

    S225.93

    A

    1002-6819(2019)-02-0010-07

    2018-09-10

    2019-01-09

    國家自然科學基金項目—新疆杏振動脫落動力學特性及振動采收機理研究(51465059)

    楊會民,助理研究員,主要從事農(nóng)業(yè)工程研究。Email:yhm_shz@163.com

    王學農(nóng),研究員,主要從事農(nóng)業(yè)工程與裝備研究。Email:xjwxn2010@ sina.com

    楊會民,散鋆龍,陳毅飛,王學農(nóng),牛長河,侯書林. 不同振動特性參數(shù)對杏樹振動響應的影響[J]. 農(nóng)業(yè)工程學報,2019,35(2):10-16. doi:10.11975/j.issn.1002-6819.2019.02.002 http://www.tcsae.org

    Yang Huimin, San Yunlong, Chen Yifei, Wang Xuenong, Niu Changhe, Hou Shulin. Influence of different vibration characteristic parameters on vibration response of apricot trees[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(2): 10-16. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.02.002 http://www.tcsae.org

    猜你喜歡
    有限元振動分析
    振動的思考
    科學大眾(2023年17期)2023-10-26 07:39:14
    隱蔽失效適航要求符合性驗證分析
    振動與頻率
    天天愛科學(2020年6期)2020-09-10 07:22:44
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    中立型Emden-Fowler微分方程的振動性
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    磨削淬硬殘余應力的有限元分析
    UF6振動激發(fā)態(tài)分子的振動-振動馳豫
    計算物理(2014年2期)2014-03-11 17:01:44
    基于SolidWorks的吸嘴支撐臂有限元分析
    箱形孔軋制的有限元模擬
    上海金屬(2013年4期)2013-12-20 07:57:18
    悠悠久久av| 丰满的人妻完整版| 亚洲国产精品sss在线观看| 日韩成人在线观看一区二区三区| 久久久久亚洲av毛片大全| 成人特级黄色片久久久久久久| 午夜两性在线视频| 老熟妇仑乱视频hdxx| 在线视频色国产色| 俄罗斯特黄特色一大片| 国产成人av教育| 中文字幕久久专区| 黑人巨大精品欧美一区二区mp4| 啦啦啦免费观看视频1| 亚洲人成77777在线视频| av福利片在线| 人妻丰满熟妇av一区二区三区| 变态另类成人亚洲欧美熟女| 男女下面进入的视频免费午夜| 成在线人永久免费视频| 亚洲人与动物交配视频| 怎么达到女性高潮| 在线观看www视频免费| 他把我摸到了高潮在线观看| 免费看美女性在线毛片视频| 99久久国产精品久久久| 淫妇啪啪啪对白视频| 小说图片视频综合网站| 久久婷婷人人爽人人干人人爱| 亚洲精品中文字幕一二三四区| 99久久精品国产亚洲精品| 日日爽夜夜爽网站| 少妇粗大呻吟视频| 亚洲熟妇熟女久久| 久9热在线精品视频| 两个人看的免费小视频| 午夜两性在线视频| 免费在线观看完整版高清| 最近最新免费中文字幕在线| 国产精品久久电影中文字幕| www日本在线高清视频| 国产野战对白在线观看| 老汉色av国产亚洲站长工具| 男女之事视频高清在线观看| 亚洲一码二码三码区别大吗| 狠狠狠狠99中文字幕| 91国产中文字幕| 首页视频小说图片口味搜索| 两个人看的免费小视频| 又黄又粗又硬又大视频| 亚洲中文字幕一区二区三区有码在线看 | 欧美极品一区二区三区四区| 中文字幕久久专区| 成人国产综合亚洲| 天天躁夜夜躁狠狠躁躁| 国产精品av视频在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美在线一区亚洲| 狂野欧美白嫩少妇大欣赏| 高清在线国产一区| 色综合婷婷激情| 一级黄色大片毛片| 亚洲专区中文字幕在线| 丁香欧美五月| 长腿黑丝高跟| 久久人妻福利社区极品人妻图片| 最近最新免费中文字幕在线| 免费看十八禁软件| 欧美日韩福利视频一区二区| 精品无人区乱码1区二区| 最新在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 国内精品久久久久精免费| 香蕉久久夜色| 国产视频一区二区在线看| 999久久久国产精品视频| 国产一区二区在线观看日韩 | 老司机午夜福利在线观看视频| 一本精品99久久精品77| 午夜日韩欧美国产| 麻豆成人av在线观看| 国产精品久久久久久亚洲av鲁大| 中文字幕人成人乱码亚洲影| 男女之事视频高清在线观看| 亚洲精品国产一区二区精华液| 国语自产精品视频在线第100页| 久久这里只有精品中国| 日韩高清综合在线| 91老司机精品| 12—13女人毛片做爰片一| 床上黄色一级片| 欧美日韩一级在线毛片| 在线国产一区二区在线| 欧美日本亚洲视频在线播放| 国产成人精品无人区| 日本在线视频免费播放| 午夜两性在线视频| 黄色视频不卡| aaaaa片日本免费| 超碰成人久久| 日本 av在线| 麻豆av在线久日| 成人18禁高潮啪啪吃奶动态图| 精品国产亚洲在线| 久久久水蜜桃国产精品网| 无人区码免费观看不卡| 女人爽到高潮嗷嗷叫在线视频| 午夜福利18| 一级黄色大片毛片| 少妇被粗大的猛进出69影院| 欧美成狂野欧美在线观看| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看 | 亚洲欧美日韩东京热| 一a级毛片在线观看| 啪啪无遮挡十八禁网站| 久久人妻av系列| 身体一侧抽搐| 精品久久久久久久毛片微露脸| 成年免费大片在线观看| 手机成人av网站| 99久久精品热视频| 午夜福利视频1000在线观看| 国产激情偷乱视频一区二区| 亚洲成人精品中文字幕电影| 亚洲五月天丁香| 制服人妻中文乱码| 18禁观看日本| 国产黄色小视频在线观看| √禁漫天堂资源中文www| 国产精品亚洲一级av第二区| 国产av不卡久久| 日本撒尿小便嘘嘘汇集6| 亚洲av中文字字幕乱码综合| 亚洲av五月六月丁香网| 色哟哟哟哟哟哟| 啦啦啦免费观看视频1| 色综合站精品国产| 母亲3免费完整高清在线观看| 别揉我奶头~嗯~啊~动态视频| 国产成人一区二区三区免费视频网站| 国产视频内射| 黄色a级毛片大全视频| 91九色精品人成在线观看| 动漫黄色视频在线观看| 美女午夜性视频免费| av有码第一页| 国产精华一区二区三区| 特大巨黑吊av在线直播| 国产亚洲精品综合一区在线观看 | 久久久久久久精品吃奶| 国产成人欧美在线观看| 看片在线看免费视频| 老汉色∧v一级毛片| 久久国产乱子伦精品免费另类| 男女做爰动态图高潮gif福利片| 丝袜美腿诱惑在线| 久久久精品国产亚洲av高清涩受| 午夜精品一区二区三区免费看| 国产亚洲av高清不卡| 中文字幕人成人乱码亚洲影| 国产伦人伦偷精品视频| 很黄的视频免费| 国产激情欧美一区二区| 亚洲专区国产一区二区| 亚洲欧美日韩无卡精品| 久久久国产精品麻豆| 免费搜索国产男女视频| 无人区码免费观看不卡| 国产av麻豆久久久久久久| 亚洲国产精品合色在线| 午夜福利高清视频| 欧美黑人精品巨大| 99国产综合亚洲精品| 精品久久久久久久末码| 别揉我奶头~嗯~啊~动态视频| 91麻豆av在线| 好看av亚洲va欧美ⅴa在| 成人av一区二区三区在线看| 午夜免费观看网址| a级毛片在线看网站| 欧美丝袜亚洲另类 | 久久午夜亚洲精品久久| 蜜桃久久精品国产亚洲av| 小说图片视频综合网站| 舔av片在线| 国内久久婷婷六月综合欲色啪| 老司机靠b影院| 国产黄色小视频在线观看| 九色国产91popny在线| av天堂在线播放| 亚洲国产欧美一区二区综合| 真人做人爱边吃奶动态| 欧美黑人精品巨大| 久久久精品国产亚洲av高清涩受| 在线a可以看的网站| 琪琪午夜伦伦电影理论片6080| 欧美日韩亚洲国产一区二区在线观看| 黄色 视频免费看| 少妇裸体淫交视频免费看高清 | 国产真人三级小视频在线观看| 高清毛片免费观看视频网站| 亚洲av片天天在线观看| 99久久精品热视频| 在线看三级毛片| 女同久久另类99精品国产91| 舔av片在线| 中文字幕最新亚洲高清| 国产av在哪里看| 欧美成人一区二区免费高清观看 | 97碰自拍视频| 久久天躁狠狠躁夜夜2o2o| 亚洲中文日韩欧美视频| 欧美午夜高清在线| 亚洲一区高清亚洲精品| 亚洲精品在线美女| 大型av网站在线播放| 中文字幕久久专区| 久久精品夜夜夜夜夜久久蜜豆 | 麻豆成人av在线观看| 亚洲狠狠婷婷综合久久图片| 精品第一国产精品| 少妇人妻一区二区三区视频| 欧美人与性动交α欧美精品济南到| 亚洲欧美日韩东京热| 亚洲熟妇熟女久久| 欧美国产日韩亚洲一区| 亚洲 国产 在线| 色精品久久人妻99蜜桃| 欧美乱码精品一区二区三区| 亚洲男人天堂网一区| 亚洲第一欧美日韩一区二区三区| 老司机靠b影院| 欧美zozozo另类| 亚洲精品av麻豆狂野| 国产不卡一卡二| 别揉我奶头~嗯~啊~动态视频| 手机成人av网站| videosex国产| 国内久久婷婷六月综合欲色啪| av天堂在线播放| 男人的好看免费观看在线视频 | 香蕉久久夜色| 天堂动漫精品| 法律面前人人平等表现在哪些方面| 此物有八面人人有两片| 俺也久久电影网| 特级一级黄色大片| 国产精品日韩av在线免费观看| 日韩欧美在线二视频| 97碰自拍视频| 后天国语完整版免费观看| 曰老女人黄片| 天天添夜夜摸| 精品国产乱子伦一区二区三区| 搡老岳熟女国产| 女警被强在线播放| 亚洲人成网站在线播放欧美日韩| 精品免费久久久久久久清纯| 午夜日韩欧美国产| 免费高清视频大片| 伦理电影免费视频| 午夜福利欧美成人| 很黄的视频免费| 妹子高潮喷水视频| 啪啪无遮挡十八禁网站| 特级一级黄色大片| 90打野战视频偷拍视频| 欧美一级a爱片免费观看看 | 亚洲欧美日韩高清在线视频| 99久久精品热视频| 国产97色在线日韩免费| 国产视频一区二区在线看| 国产欧美日韩精品亚洲av| 欧美乱妇无乱码| 九色国产91popny在线| 亚洲 国产 在线| 一夜夜www| 国产片内射在线| 嫩草影视91久久| 老鸭窝网址在线观看| 免费在线观看日本一区| 此物有八面人人有两片| 亚洲av熟女| 国产片内射在线| 欧美三级亚洲精品| 变态另类丝袜制服| 美女 人体艺术 gogo| 不卡av一区二区三区| 在线永久观看黄色视频| 午夜老司机福利片| 又黄又粗又硬又大视频| 在线永久观看黄色视频| 欧美日韩福利视频一区二区| 亚洲最大成人中文| 在线永久观看黄色视频| 18禁裸乳无遮挡免费网站照片| 国产成人av激情在线播放| 精品电影一区二区在线| 丝袜美腿诱惑在线| 老司机午夜福利在线观看视频| 黑人巨大精品欧美一区二区mp4| 天堂√8在线中文| 女警被强在线播放| netflix在线观看网站| 久久精品aⅴ一区二区三区四区| 日本 欧美在线| 精品久久久久久久久久免费视频| 国内精品一区二区在线观看| 欧美午夜高清在线| 亚洲成人久久爱视频| 男女那种视频在线观看| 亚洲欧美日韩无卡精品| 中国美女看黄片| 50天的宝宝边吃奶边哭怎么回事| 国产爱豆传媒在线观看 | 亚洲av成人一区二区三| 亚洲 国产 在线| 人成视频在线观看免费观看| 成人18禁高潮啪啪吃奶动态图| 久久久久久九九精品二区国产 | 欧美日本视频| 日韩欧美精品v在线| 久久久久久久精品吃奶| 久久久久久久久免费视频了| 老司机靠b影院| 欧美黄色片欧美黄色片| 精品日产1卡2卡| 久久国产精品影院| 亚洲av成人av| 成熟少妇高潮喷水视频| 亚洲熟女毛片儿| 欧美黑人欧美精品刺激| 精品熟女少妇八av免费久了| 熟妇人妻久久中文字幕3abv| 精品国产美女av久久久久小说| 黄片大片在线免费观看| 9191精品国产免费久久| 丁香欧美五月| 久久久久国内视频| 欧美日本视频| 国产私拍福利视频在线观看| а√天堂www在线а√下载| 午夜激情福利司机影院| 国产免费男女视频| 少妇被粗大的猛进出69影院| 青草久久国产| 日日摸夜夜添夜夜添小说| 美女大奶头视频| 亚洲欧美一区二区三区黑人| 正在播放国产对白刺激| 欧美日韩黄片免| 天堂√8在线中文| av片东京热男人的天堂| 岛国在线观看网站| 国产免费av片在线观看野外av| 91国产中文字幕| 亚洲专区字幕在线| 亚洲激情在线av| 国产熟女午夜一区二区三区| 又黄又粗又硬又大视频| av在线播放免费不卡| 亚洲成av人片免费观看| 看片在线看免费视频| 午夜影院日韩av| 在线观看免费日韩欧美大片| 色老头精品视频在线观看| 日本精品一区二区三区蜜桃| 国产精品亚洲av一区麻豆| 不卡av一区二区三区| 免费在线观看完整版高清| 久久这里只有精品19| 久久久久国产精品人妻aⅴ院| 国产精品,欧美在线| 久久久久性生活片| 国产精品98久久久久久宅男小说| 免费观看人在逋| 日韩大尺度精品在线看网址| 在线视频色国产色| 欧美成狂野欧美在线观看| 在线十欧美十亚洲十日本专区| www.999成人在线观看| 老司机福利观看| 麻豆久久精品国产亚洲av| 成人午夜高清在线视频| 精品国内亚洲2022精品成人| avwww免费| 国产精品美女特级片免费视频播放器 | 黄片大片在线免费观看| 亚洲人成网站高清观看| 小说图片视频综合网站| 舔av片在线| 亚洲在线自拍视频| 日本五十路高清| 欧美在线一区亚洲| 丁香六月欧美| 18禁黄网站禁片免费观看直播| 午夜福利在线观看吧| 一级毛片高清免费大全| 99re在线观看精品视频| 精品欧美国产一区二区三| 757午夜福利合集在线观看| 亚洲熟妇中文字幕五十中出| 成年女人毛片免费观看观看9| 欧美日韩黄片免| 性色av乱码一区二区三区2| 两个人看的免费小视频| 成人国产一区最新在线观看| 亚洲五月婷婷丁香| 特级一级黄色大片| 久久亚洲精品不卡| bbb黄色大片| 久久精品亚洲精品国产色婷小说| 国产一区二区在线av高清观看| 国产av一区在线观看免费| 亚洲精品一区av在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品合色在线| 天堂√8在线中文| www国产在线视频色| 少妇粗大呻吟视频| 1024视频免费在线观看| 在线观看舔阴道视频| 亚洲av美国av| 色综合欧美亚洲国产小说| 精品电影一区二区在线| 丝袜人妻中文字幕| 欧美高清成人免费视频www| 他把我摸到了高潮在线观看| 久久久久国产一级毛片高清牌| 悠悠久久av| 18禁美女被吸乳视频| 国产精品一区二区免费欧美| 精品人妻1区二区| 好男人在线观看高清免费视频| 在线a可以看的网站| 香蕉丝袜av| 亚洲中文日韩欧美视频| 黄色视频不卡| 亚洲国产中文字幕在线视频| 国产三级中文精品| 在线观看免费日韩欧美大片| 欧美乱色亚洲激情| 亚洲精品色激情综合| 亚洲欧美精品综合久久99| 国产成人系列免费观看| 亚洲成a人片在线一区二区| 欧美精品亚洲一区二区| 麻豆国产97在线/欧美 | 搞女人的毛片| 黄色视频不卡| 啦啦啦观看免费观看视频高清| 在线永久观看黄色视频| 国产单亲对白刺激| 老熟妇仑乱视频hdxx| 99riav亚洲国产免费| 欧美性长视频在线观看| 最近最新中文字幕大全免费视频| 在线观看www视频免费| 在线观看美女被高潮喷水网站 | 精品国内亚洲2022精品成人| 亚洲国产看品久久| 老司机午夜十八禁免费视频| 成人av一区二区三区在线看| 毛片女人毛片| 午夜两性在线视频| 亚洲精华国产精华精| 日韩精品中文字幕看吧| 欧美日韩亚洲综合一区二区三区_| 中文字幕高清在线视频| 亚洲精品粉嫩美女一区| 精品一区二区三区av网在线观看| 1024香蕉在线观看| 久久伊人香网站| 在线观看一区二区三区| www国产在线视频色| 久久久国产欧美日韩av| 国产成人av教育| 久久精品人妻少妇| bbb黄色大片| 国内少妇人妻偷人精品xxx网站 | 久久久国产成人精品二区| 一级a爱片免费观看的视频| 午夜久久久久精精品| 国产黄a三级三级三级人| 中文资源天堂在线| 人成视频在线观看免费观看| 老汉色∧v一级毛片| 亚洲免费av在线视频| 国产免费av片在线观看野外av| 国产精品电影一区二区三区| 亚洲七黄色美女视频| 亚洲国产高清在线一区二区三| 香蕉国产在线看| 黄色 视频免费看| 精品欧美国产一区二区三| 欧美不卡视频在线免费观看 | 99精品在免费线老司机午夜| 18禁国产床啪视频网站| 久久热在线av| 十八禁网站免费在线| 精品国内亚洲2022精品成人| 久久香蕉国产精品| 亚洲欧美日韩高清在线视频| 精品乱码久久久久久99久播| 色av中文字幕| av天堂在线播放| 国产av又大| 日本黄色视频三级网站网址| 制服丝袜大香蕉在线| 色综合婷婷激情| 久久久久久久午夜电影| 日日摸夜夜添夜夜添小说| 国产午夜精品久久久久久| 日日夜夜操网爽| 久久香蕉国产精品| 看片在线看免费视频| 露出奶头的视频| 亚洲最大成人中文| 一边摸一边抽搐一进一小说| 日本黄大片高清| 成人av在线播放网站| 一区二区三区高清视频在线| av福利片在线| 久久精品人妻少妇| 精品无人区乱码1区二区| 免费看a级黄色片| 女人被狂操c到高潮| 欧美在线黄色| 国产成人啪精品午夜网站| 久久这里只有精品中国| www.熟女人妻精品国产| 国产精品,欧美在线| 亚洲在线自拍视频| 19禁男女啪啪无遮挡网站| 真人做人爱边吃奶动态| 老鸭窝网址在线观看| 91字幕亚洲| 久久伊人香网站| 90打野战视频偷拍视频| 精品不卡国产一区二区三区| av片东京热男人的天堂| 欧美黄色片欧美黄色片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产探花在线观看一区二区| 国产成人精品无人区| ponron亚洲| 国产久久久一区二区三区| 琪琪午夜伦伦电影理论片6080| 午夜激情av网站| 欧美成狂野欧美在线观看| 一区二区三区高清视频在线| 国产精品久久久人人做人人爽| 亚洲五月天丁香| 91国产中文字幕| 日日爽夜夜爽网站| 亚洲国产欧美人成| 国产亚洲精品av在线| 亚洲精品中文字幕一二三四区| 老熟妇乱子伦视频在线观看| av超薄肉色丝袜交足视频| 国产精品免费视频内射| 欧美乱妇无乱码| 亚洲 国产 在线| 久久久久久亚洲精品国产蜜桃av| av国产免费在线观看| 99热这里只有精品一区 | 激情在线观看视频在线高清| 好男人在线观看高清免费视频| 国产不卡一卡二| 成人精品一区二区免费| 日韩大码丰满熟妇| 村上凉子中文字幕在线| 怎么达到女性高潮| 啪啪无遮挡十八禁网站| av天堂在线播放| 国产在线精品亚洲第一网站| а√天堂www在线а√下载| 国产亚洲av高清不卡| 久久久国产成人免费| 99国产极品粉嫩在线观看| 桃色一区二区三区在线观看| 亚洲欧美日韩无卡精品| 亚洲avbb在线观看| 美女黄网站色视频| 观看免费一级毛片| 亚洲精品美女久久久久99蜜臀| 男女那种视频在线观看| 床上黄色一级片| 欧美高清成人免费视频www| 亚洲七黄色美女视频| 久久香蕉精品热| 国产成人aa在线观看| 在线视频色国产色| 高清在线国产一区| 日韩三级视频一区二区三区| 亚洲成av人片在线播放无| 国产成人精品久久二区二区91| 国产视频一区二区在线看| 在线观看免费午夜福利视频| 欧美日本亚洲视频在线播放| 日韩欧美在线二视频| 老汉色av国产亚洲站长工具| 搡老熟女国产l中国老女人| 亚洲一码二码三码区别大吗| 成人特级黄色片久久久久久久| 国产成人av教育| 99re在线观看精品视频| 久久热在线av| 午夜福利免费观看在线| 激情在线观看视频在线高清| 给我免费播放毛片高清在线观看| 又粗又爽又猛毛片免费看| 制服人妻中文乱码| 欧美zozozo另类|