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

    基于模態(tài)置信度準則的插秧機支撐臂模態(tài)分析與結(jié)構(gòu)優(yōu)化

    2018-10-10 07:22:50陳凱康姬江濤曾欣悅
    農(nóng)業(yè)工程學報 2018年18期
    關(guān)鍵詞:插秧機固有頻率測點

    金 鑫,陳凱康,姬江濤,龐 靖,高 頌,曾欣悅

    ?

    基于模態(tài)置信度準則的插秧機支撐臂模態(tài)分析與結(jié)構(gòu)優(yōu)化

    金 鑫1,2,陳凱康1,姬江濤1,2※,龐 靖1,高 頌1,曾欣悅1

    (1.河南科技大學農(nóng)業(yè)裝備工程學院,洛陽 471003;2. 機械裝備先進制造河南省協(xié)同創(chuàng)新中心,洛陽 471003)

    針對高速插秧機栽插機構(gòu)動力傳輸過程中的振動問題,該文以2ZG-6DK插秧機為研究對象,分析了栽插系統(tǒng)動力傳輸方式及工作原理;利用SolidWorks軟件對插秧機動力傳輸系統(tǒng)中關(guān)鍵部件支撐臂進行建模,將模型導入ANSYS Workbench并結(jié)合Lanczos Method解算方法求解模態(tài)固有頻率和振型,在此基礎(chǔ)上開展基于MAC(modal assurance criterion)準則優(yōu)化的模態(tài)試驗,驗證了有限元理論分析的準確性。為使支撐臂固有頻率避開外部激振頻率,在分析外部頻率激振特點的基礎(chǔ)上,基于ISIGHT多學科軟件平臺,采用序列二次規(guī)劃法對支撐臂結(jié)構(gòu)參數(shù)進行優(yōu)化。研究結(jié)果表明:優(yōu)化后支撐臂側(cè)壁腔體厚度5.7 mm、橫梁寬度42.0 mm、臂長497.0 mm,前4階模態(tài)頻率分別調(diào)整至135.17、204.23、483.14和702.32 Hz,均可避開插秧機汽油發(fā)動機激振頻率范圍86.67~120 Hz。優(yōu)化后1階頻率下振動幅度衰減最為明顯,振動幅值最高下降9.4%,支撐臂低頻振動特性得到明顯改善。研究結(jié)果可為插秧機的振動特性分析與減振設(shè)計提供參考。

    農(nóng)業(yè)機械;振動;優(yōu)化;模態(tài)分析;MAC準則

    0 引 言

    分插機構(gòu)支撐臂又稱為栽植鏈條箱,是水稻插秧機動力傳輸系統(tǒng)核心工作部件,其工作性能直接影響整機作業(yè)質(zhì)量和效率[1-3]。插秧機高速作業(yè)條件下受路面激勵、內(nèi)燃機振動激勵、傳動系統(tǒng)振動激勵和沖擊載荷等影響,當自身固有頻率接近動載荷頻率時易發(fā)生共振現(xiàn)象[4-7]。共振產(chǎn)生的動應(yīng)力不僅影響駕駛舒適性,而且會造成與其配套的零部件因疲勞而發(fā)生脆性斷裂,影響機械設(shè)備的性能和工作壽命[8-11]。因此,開展以減振降噪技術(shù)為核心的插秧機分插機構(gòu)支撐臂的振動研究具有重要意義,可為深入探究水稻栽插作業(yè)關(guān)鍵部件振動機理提供理論借鑒。

    目前,國內(nèi)外針對農(nóng)機裝備作業(yè)工況下振動較大等問題,在振動特性分析、駕駛座椅振動與舒適性分析、模態(tài)分析與試驗、傳感器測點優(yōu)化等方面已有相關(guān)研究。李耀明等[9-10]通過有限元計算模態(tài)分析與自由模態(tài)驗證相結(jié)合的方法求解了履帶聯(lián)合收獲機機架的固有頻率和振型,對機架進行結(jié)構(gòu)優(yōu)化與試驗后達到了減振降噪的目的。姚艷春等[11]分析車架田間振動特性,并以提高1階扭轉(zhuǎn)頻率為目標優(yōu)化車架結(jié)構(gòu),優(yōu)化后車架無故障時間由20 h提高到60 h。Wang JL等[12]利用ANSYS有限元分析軟件,獲得了移栽機3種關(guān)鍵結(jié)構(gòu)件固有頻率、阻尼比和模態(tài)振型,對3種結(jié)構(gòu)的模態(tài)參數(shù)進行比較,為移栽機整體結(jié)構(gòu)設(shè)計提供了理論依據(jù)。張華標等[13]以插秧機分插機構(gòu)支撐臂為研究對象,利用有限元分析軟件對其受力和振動模態(tài)進行分析得出其6階模態(tài)參數(shù),為分插機構(gòu)設(shè)計提供理論基礎(chǔ)。

    李復輝等[14]通過有限元分析軟件對插秧機分插機構(gòu)行星輪系進行模態(tài)計算,得出行星輪前15階模態(tài)固有頻率和振型云圖,分析證明行星輪工作過程中不會發(fā)生共振現(xiàn)象。張翔等[15]以LNC8型插秧機車架為研究對象,結(jié)合有限元分析與模態(tài)試驗結(jié)果對影響較大的4個主要因素進行正交試驗設(shè)計,優(yōu)化后試驗對象振動特性有明顯改善。

    當前,在插秧機械等農(nóng)業(yè)機械裝備振動測試與特征分析方面已取得相關(guān)研究成果[16-19],但尚未完全清楚插秧機分插機構(gòu)支撐臂作業(yè)過程中振動特性及其振動規(guī)律,尤其對在高速栽插狀態(tài)下支撐臂的振動特性、振動主頻與模態(tài)振型之間的對應(yīng)關(guān)系研究相對較少,缺乏對于農(nóng)業(yè)機械振動模態(tài)試驗傳感器布置位點的研究。

    本文在現(xiàn)有研究的基礎(chǔ)上,深入探究插秧機支撐臂振動特性,利用三維建模軟件結(jié)合有限元Lanczos Method解算方法求解出支撐臂模態(tài)參數(shù),通過模態(tài)試驗對比驗證有限元模態(tài)的準確性。在對比分析支撐臂固有頻率與外部激振頻率特點的基礎(chǔ)上,基于ISIGHT多學科軟件平臺,采用序列二次規(guī)劃法對支撐臂結(jié)構(gòu)參數(shù)進行了優(yōu)化。本文旨在為栽植機械復雜結(jié)構(gòu)件振動特性分析提供方法,以期為高速栽植機械的振動特性分析與減振設(shè)計提供參考。

    1 高速插秧機栽插系統(tǒng)組成及工作原理

    高速插秧機栽插系統(tǒng)結(jié)構(gòu)簡圖如圖1a所示,主要由插秧支撐臂、插秧旋轉(zhuǎn)箱、推秧裝置等組成。其中支撐臂通過螺栓固定在機架上,支撐臂左右兩側(cè)布置2個插秧旋轉(zhuǎn)箱,插秧旋轉(zhuǎn)箱內(nèi)部行星輪系將動力傳輸至推秧裝置,每個插秧旋轉(zhuǎn)箱體上安裝有2個推秧裝置。

    栽插系統(tǒng)作業(yè)過程中,動力通過鏈條傳動至中心軸,旋轉(zhuǎn)箱內(nèi)部的行星輪和中間輪隨旋轉(zhuǎn)箱繞中心軸作勻速圓周運動。位于插秧旋轉(zhuǎn)箱兩側(cè)的推秧裝置在行星輪系帶動下構(gòu)成相對運動,完成周期性推秧動作。在插秧機工作過程中,支撐臂是動力傳輸機構(gòu)的主要承載部件,其工作穩(wěn)定性直接影響插秧的質(zhì)量和效率。支撐臂模型如圖1b所示。

    圖1 高速插秧機栽插試驗臺及插秧機支撐臂示意圖

    2 支撐臂有限元模態(tài)分析

    2.1 有限元模型建立

    南通富來威2ZG-6DK插秧機支撐臂外形尺寸長500 mm,寬98 mm,高131 mm,主要由ZL101鋁合金經(jīng)殼型鑄造后調(diào)質(zhì)處理而成。建立支撐臂有限元模型如圖2所示。

    圖2 支撐臂有限元模型

    利用三維軟件對支撐臂進行虛擬建模,將模型導入ANSYS Workbench中。因有限元建模直接影響計算結(jié)果,因此需考慮模型的簡化、網(wǎng)格數(shù)目、單元類型、網(wǎng)格質(zhì)量等因素。首先對影響模型精度的關(guān)鍵因素進行幾何清理,然后對支撐臂有限元建模作出如下簡化:1)不考慮半徑較小倒角和孔對整體的影響;2)忽略支撐臂上的凸起、側(cè)筋等非承載構(gòu)件影響;3)焊接處材料特性按照周圍構(gòu)件特性計算。

    2.2 網(wǎng)格劃分與模型設(shè)置

    如圖2所示,在對所建立的支撐臂模型進行劃分時,定義模型材料為鋁,其材料特性設(shè)置彈性模量70 GPa,泊松比0.33,質(zhì)量密度2 680 kg/m3,采用6面體網(wǎng)格,網(wǎng)格尺寸為2 mm。為使模型接近實際,采用約束模態(tài),根據(jù)實物結(jié)構(gòu),約束位置為2個軸心位置和鉸接位置。最終得到單元總數(shù)62 208個,節(jié)點總數(shù)108 698的有限元模型。

    2.3 有限元模態(tài)分析

    模態(tài)分析用于確定結(jié)構(gòu)的模態(tài)參數(shù),即固有頻率和振型。在工程應(yīng)用中,低階模態(tài)對支撐臂動態(tài)特性影響較大,分析支撐臂低階振型對預(yù)測結(jié)構(gòu)穩(wěn)定性有重要作用[20]。高頻振動一般遠離支撐臂低階振動頻率,對部件疲勞破壞明顯低于低頻振動部分。故選取Lanczos Method解算方法,求解得到前4 階固有頻率及對應(yīng)振型云圖,如圖3所示。

    圖3 插秧機支撐臂有限元模態(tài)固有頻率及模態(tài)振型

    由圖3可知,支撐臂1階模態(tài)振型主要以沿y軸正方向彎曲為主。2階振型主要表現(xiàn)為支撐臂兩端沿y軸負方向彎曲變形。3階模態(tài)振型主要表現(xiàn)為扭轉(zhuǎn)變形。4階模態(tài)振型主要表現(xiàn)為凸型結(jié)構(gòu)塊和動力輸入鏈條箱反向彎曲變形。

    3 基于MAC矩陣優(yōu)化的模態(tài)試驗

    3.1 模態(tài)試驗系統(tǒng)組成

    模態(tài)試驗是基于外界激勵和系統(tǒng)響應(yīng)的動態(tài)測試,該測試通過系統(tǒng)輸入的激振力和輸出的響應(yīng)數(shù)據(jù)經(jīng)信號處理與參數(shù)識別來確定系統(tǒng)模態(tài)參數(shù)。在模態(tài)試驗中常用激振器或模態(tài)力錘迫使結(jié)構(gòu)產(chǎn)生強迫振動,通過分析2點之間的機械導納函數(shù),識別出結(jié)構(gòu)物的的頻率、阻尼和振型等模態(tài)參數(shù),為結(jié)構(gòu)系統(tǒng)的振動特性分析以及結(jié)構(gòu)動力特性的優(yōu)化設(shè)計提供依據(jù)。

    本文模態(tài)試驗系統(tǒng)由數(shù)據(jù)采集模塊和模態(tài)分析模塊組成。其中數(shù)據(jù)采集模塊由力錘、三軸加速度傳感器和DH5902數(shù)據(jù)采集儀組成;模態(tài)分析模塊由DH5902動態(tài)信號分析系統(tǒng)和DHMA模態(tài)分析軟件組成。圖4為插秧機支撐臂模態(tài)試驗流程圖。

    圖4 插秧機支撐臂模態(tài)試驗流程圖

    3.2 試驗設(shè)備與方法

    模態(tài)試驗采用DH5902動態(tài)信號測試分析系統(tǒng)對2ZG-6DK插秧機支撐臂振動特性進行測試,試驗所需儀器設(shè)備及參數(shù)如表1所示。

    表1 模態(tài)試驗用測試設(shè)備

    如圖5所示,本文模態(tài)試驗采用脈沖激勵法,通過同時測量輸入與輸出信號,即通過同時測量激勵和響應(yīng),推知系統(tǒng)的固有特性。試驗過程為:利用遠小于被測物體剛度的黑色彈簧將支撐臂吊起,保證支撐臂平衡懸置于空中。將三軸加速度計按要求布置在支撐臂表面,順次用導線連接傳感器及整個測試系統(tǒng)。模態(tài)力錘敲擊激振點后,數(shù)據(jù)同時由模態(tài)力錘和加速度傳感器傳輸至數(shù)據(jù)采集系統(tǒng)中,并得到系統(tǒng)響應(yīng)函數(shù)。模態(tài)分析系統(tǒng)分析整體系統(tǒng)響應(yīng)函數(shù)后得到插秧機支撐臂整體振動特性。

    圖5 支撐臂模態(tài)試驗

    3.3 基于MAC置信準則的測點優(yōu)化布置

    傳統(tǒng)模態(tài)試驗依靠經(jīng)驗選取測點或均勻布點時,存在對經(jīng)驗高度依賴且試驗效率較低的弊端[21-23]。本文所測插秧機支撐臂結(jié)構(gòu)較不規(guī)則,在實際模態(tài)試驗中,由于有限的測點數(shù)及測試噪聲引起模態(tài)向量不正交等問題易對試驗準確性產(chǎn)生影響。引入模態(tài)置信度(modal assurance criterion, MAC)矩陣對測點布置方案進行優(yōu)化設(shè)計,可在保證模態(tài)試驗精度的前提下,精簡測試點數(shù)量,從而達到提高試驗精度與效率的目的[24-28]。模態(tài)置信度準則是一種傳感器優(yōu)化配置準則,其矩陣表達式為

    以支撐臂凸型結(jié)構(gòu)塊為例,在測點優(yōu)化過程中,通過有限元分析獲得前4階模態(tài)振型矩陣,然后根據(jù)MAC矩陣相結(jié)合的方法在Matlab中進行迭代計算,得到MAC矩陣最大非對角元元素的最大值與自由度數(shù)的變化曲線如圖6所示。

    圖6 MAC矩陣最大非對角元元素的最大值與自由度數(shù)的變化曲線

    由圖6可知,當自由度為3時,該矩陣非對角元最大值達到最小,即最少在凸型結(jié)構(gòu)塊上邊緣布置3個傳感器測點方可基本滿足測試要求。凸型結(jié)構(gòu)塊上邊緣布置3個傳感器后即可識別出被測部件固有頻率,為同時獲得多階模態(tài)振型,可在MAC置信準則基礎(chǔ)上結(jié)合香農(nóng)定理適當增加測點以達到準確識別振型的要求。

    使用香農(nóng)定理進行增設(shè)測點,測得支撐臂凸型結(jié)構(gòu)塊模態(tài)的最高頻率,估計該最高頻率的半波長;在半波長的每個節(jié)點上布置1個傳感器;再在半波長上等均布2個傳感器。支撐臂凸型結(jié)構(gòu)塊的頂端邊緣長80 mm,通過有限元分析獲得第4階固有頻率下滑座的一條長邊的節(jié)點位移,再將其進行曲線擬合得到滑座變形最大的y方向上的位移波形,從波形估計支撐臂凸型結(jié)構(gòu)塊的波長為75 mm,需在2節(jié)點之間增設(shè)1個測點使得振型更好分辨。同理,對支撐臂其他關(guān)鍵結(jié)構(gòu)也通過該種方法進行測點優(yōu)化。

    優(yōu)化后測點分布如圖7所示,共有78個測點,該模型結(jié)構(gòu)能較好地定義插秧機支撐臂輪廓形狀。圖中編號為有限元分析中對應(yīng)的節(jié)點編號。圖中紫色節(jié)點為依據(jù)MAC置信準則得到的節(jié)點,綠色節(jié)點為結(jié)合香農(nóng)定理增加的節(jié)點。通過對比圖7和圖3的前4階振型,發(fā)現(xiàn)通過此種方法選擇的測點大都分布在零件變形較大的部位,因此這樣選擇的測點對模態(tài)振型的貢獻最大,且滿足振型識別的要求。

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

    為保證識別出的頻響函數(shù)具有較好的一致性,并且不遺漏模態(tài),頻響函數(shù)測試采用自由激振的方法,選用2根彈簧將支撐臂懸掛呈自由狀態(tài)。為保證系統(tǒng)的可辨識性,模態(tài)試驗激振點不能靠近節(jié)點或節(jié)線太近。本試驗選取2個激振點,豎直方向(軸方向)激勵點選取在凸型結(jié)構(gòu)塊頂端,側(cè)向(軸方向)激勵點選取在回轉(zhuǎn)箱定位塊處。

    將經(jīng)由動態(tài)信號采集分析系統(tǒng)收集的數(shù)據(jù)導入DHMA模態(tài)分析軟件進行分析,并將模態(tài)參數(shù)數(shù)據(jù)匹配于所建立的試驗?zāi)P椭校罱K得到插秧機支撐臂固有頻率和振型等模態(tài)參數(shù)。一般而言低階振動對結(jié)構(gòu)的動態(tài)影響較大,低階振型決定結(jié)構(gòu)的動態(tài)特性,故提取插秧機支撐臂前4階固有頻率和振型,結(jié)果如圖8所示。

    注:圖中數(shù)字為試驗測點標號。

    3.5 試驗與仿真結(jié)果對比分析

    將有限元模態(tài)分析與試驗?zāi)B(tài)結(jié)果進行對比分析,結(jié)果如表2所示。由表2可知,插秧機支撐臂有限元模態(tài)與試驗?zāi)B(tài)結(jié)果基本一致,前4階模態(tài)振型以彎曲為主,阻尼比均不超過0.5%,固有頻率相對誤差最大值為3.2%,有限元模態(tài)分析數(shù)據(jù)可靠。

    圖8 插秧機支撐臂試驗?zāi)B(tài)下前4階固有頻率和振型

    表2 有限元分析與試驗測試的模態(tài)結(jié)果對比

    4 支撐臂結(jié)構(gòu)優(yōu)化

    4.1 激振特點分析

    本機型采用的發(fā)動機為風冷4沖程2汽缸OHV汽油機,其曲軸曲拐呈180°布置。當插秧機正常工作時,2階慣性力是其激振力,轉(zhuǎn)速達到2 600~3 600 r/min時,激振頻率為86.67~120 Hz。插秧機回轉(zhuǎn)箱在工作狀態(tài)下做回轉(zhuǎn)運動,其轉(zhuǎn)動速度處于60~110 r/min之間,激振頻率為3~3.6 Hz。實際測得插秧機傳動系統(tǒng)主軸正常工作轉(zhuǎn)速為400~600 r/min,激振頻率為13.33~20 Hz。

    通過觀察支撐臂1階模態(tài)振型云圖,結(jié)合模態(tài)試驗中物件變形狀態(tài)可知當外部激振頻率接近101.29 Hz時,高速插秧機支撐臂易產(chǎn)生共振。發(fā)動機的激振頻率范圍在86.67~120 Hz,插秧機支撐臂1階固有頻率恰在該頻率范圍內(nèi),因此在插秧機完成高速栽插作業(yè)時,發(fā)動機所引起的支撐臂共振現(xiàn)象會對秧苗栽插作業(yè)質(zhì)量產(chǎn)生影響。若對支撐臂結(jié)構(gòu)進行減振設(shè)計,需提高支撐臂整體固有頻率,使其1階固有頻率更加遠離激振頻率。

    4.2 支撐臂減振優(yōu)化設(shè)計

    根據(jù)有限元模態(tài)分析和現(xiàn)場模態(tài)試驗結(jié)果,以支撐臂工作狀態(tài)下1階固有頻率為優(yōu)化對象,對支撐臂結(jié)構(gòu)進行改進以避免共振現(xiàn)象發(fā)生。支撐臂振動特性優(yōu)化過程為復雜的多參數(shù)、跨學科優(yōu)化設(shè)計,需調(diào)用多軟件協(xié)同設(shè)計求得最優(yōu)解[29-32]。本文通過ISIGHT優(yōu)化平臺快速集成和耦合各種仿真軟件,在ISIGHT優(yōu)化平臺中集成三維模型和ANSYS源文件,采用NLPQL(序列二次規(guī)劃法),以支撐臂1階固有頻率為優(yōu)化目標,對支撐臂側(cè)壁腔體厚度,橫梁寬度,臂長進行優(yōu)化。將系統(tǒng)離散成有限多個殼體單元,其單元剛度矩陣為:

    設(shè)支撐臂材料結(jié)構(gòu)是均勻異性材料,則單元的彈性矩陣為

    優(yōu)化過程為將系統(tǒng)固有頻率求解過程集成在多學科優(yōu)化設(shè)計ISIGHT軟件框架中,并做好輸入與輸出文件的解析。ISIGHT軟件中根據(jù)序列二次規(guī)劃法對設(shè)計變量和目標響應(yīng)進行優(yōu)化,按照設(shè)定次數(shù)進行迭代計算并把過程集成為一個可自動執(zhí)行的系統(tǒng)。此系統(tǒng)優(yōu)化過程流程圖如圖9所示。

    圖9 ISIGHT優(yōu)化設(shè)計流程圖

    設(shè)計要求在保證提高支撐臂低階固有頻率的同時,盡量避開整機各總成及部件的振動的激振頻率帶。在滿足支撐臂1階固有頻率較大和質(zhì)量較小的條件下,得到支撐臂側(cè)壁腔體厚度5.7 mm、橫梁寬度42.0 mm、臂長497.0 mm。在此參數(shù)設(shè)置下支撐臂前4階有限元模態(tài)參數(shù)如表3所示。其中支撐臂前4階固有頻率分別調(diào)整至 135.17、204.23、483.14和702.32 Hz,均避開了插秧機汽油發(fā)動機86.67~120 Hz的激振頻率范圍。

    表3 優(yōu)化后支撐臂前4階有限元模態(tài)參數(shù)

    4.3 減振設(shè)計效果驗證

    為驗證支撐臂各參數(shù)優(yōu)化的效果,對優(yōu)化前后支撐臂關(guān)鍵部位(前進方向)、(橫向)、(縱向)的振動信號進行采樣分析。選取6個關(guān)鍵部位測試優(yōu)化前后振動參數(shù),現(xiàn)場試驗如圖10所示。

    圖10 優(yōu)化前后插秧機支撐臂振動參數(shù)測試

    優(yōu)化前后支撐臂的臂長變化較小,因此不考慮因臂長變化而造成傳感器測點位置發(fā)生變化帶來的影響。測試地點為河南科技大學農(nóng)業(yè)機械作業(yè)狀態(tài)監(jiān)測試驗室,測試選用356A16三軸加速度傳感器、DHDAS5902動態(tài)信號分析系統(tǒng)、DH5902多通道信號采集儀進行數(shù)據(jù)采樣與分析,試驗時發(fā)動機轉(zhuǎn)速始終在額定轉(zhuǎn)速(2 800 r/min)以保證各部件工作參數(shù)不變。根據(jù)插秧機支撐臂的結(jié)構(gòu)特點,設(shè)置采樣頻率為2 560 Hz,時域點數(shù)為4 096,頻域點數(shù)1 600。同種參數(shù)設(shè)置下采集5組振動信號數(shù)據(jù),取數(shù)據(jù)置信度較好一組數(shù)據(jù)分析所得各測點前4階振動幅值如表4所示。

    由表4可知,優(yōu)化前支撐臂在測點3、測點4、測點5處振幅較大,優(yōu)化后1階頻率下振動幅度衰減最為明顯。其中測點3處優(yōu)化后135.17 Hz頻率下振幅下降9.4%,說明結(jié)構(gòu)優(yōu)化后支撐臂振動特性有明顯改善,尤其在低階頻率范圍內(nèi)振動得到了很好控制。此方法可推廣至插秧機整機振動特性優(yōu)化中,通過改進關(guān)鍵部件結(jié)構(gòu)參數(shù),使固有頻率遠離外部激振頻率范圍即可避免共振發(fā)生。

    表4 優(yōu)化前后支撐臂各測點振動幅值與頻率對比

    支撐臂在、、方向上的振動幅值均在測點4處達到最大值,此時振動頻率約為101 Hz,振動幅值分別為2.98、2.60、3.28 m/s2,說明測點4處振動最為顯著,結(jié)構(gòu)優(yōu)化后,其在3個方向振動幅值均有所降低,但相較于其他部位振幅仍然較高,說明在栽插過程中外部激振由動力傳動系統(tǒng)、移箱機構(gòu)等部件經(jīng)機架傳至支撐臂,引起了支撐臂明顯振動。在設(shè)計時可對測點4處額外增加阻尼塊使振動衰減,也可在傳遞振動的機架與其他部件之間增加隔振裝置以達到減振目的。

    優(yōu)化后測點6處、、方向振幅衰減均不顯著,在方向1階(135.17 Hz)最大振動幅度衰減也僅有1.82%,在部分階次下振動幅度略有升高。綜合分析實際作業(yè)工況下支撐臂測點6處位置及與其他部件連接關(guān)系,發(fā)現(xiàn)測點6處為動力輸入軸與鏈輪連接部位,其振動受動力輸入軸穩(wěn)定性影響較大。且測點6處內(nèi)部結(jié)構(gòu)復雜,優(yōu)化后其局部模態(tài)特性變化不大,亦有可能其固有頻率仍在外部激振頻率范圍內(nèi),此時應(yīng)開展對局部模態(tài)的分析與改進,以降低該部位的振動幅值。

    圖11 支撐臂應(yīng)力云圖

    對優(yōu)化后的支撐臂進行靜力學分析,支撐臂主要受力有:自身重力52.36 N、插秧旋轉(zhuǎn)箱19.32 N、推秧裝置8.65 N、傳動部件重力34.18 N。將優(yōu)化后的支撐臂重新進行三維建模與網(wǎng)格劃分,把相應(yīng)作用力施加在特定位置處,并對支撐臂動力輸入連接處和與機架連接處設(shè)置固定約束,靜力學應(yīng)力仿真結(jié)果如圖11所示。由圖11可知,支撐臂整體受力較小,參照實際作業(yè)條件施加壓力與載荷后未產(chǎn)生明顯形變與裂紋。結(jié)構(gòu)優(yōu)化后應(yīng)力最大處出現(xiàn)在回轉(zhuǎn)箱定位塊位置,應(yīng)力值為57 MPa,小于材料屈服極限211 MPa,故改進后的支撐臂強度滿足設(shè)計要求。

    5 結(jié) 論

    1)通過三維建模軟件建立2ZG-6DK插秧機支撐臂模型,在ANSYS中選取Lanczos Method解算方法求解出支撐臂前4階模態(tài)固有頻率和振型,并與通過MAC模態(tài)置信度準則優(yōu)化后的試驗?zāi)B(tài)結(jié)果對比,驗證了理論模型與模態(tài)測試方法的準確性。

    2)對比分析了支撐臂固有頻率與外部激振頻率之間的關(guān)系,其中支撐臂1階固有頻率101.29 Hz在發(fā)動機激振頻率86.67~120 Hz范圍內(nèi),說明在插秧機高速栽插作業(yè)工況下,發(fā)動機激振會引起支撐臂產(chǎn)生共振現(xiàn)象。

    3)為調(diào)整插秧機支撐臂的固有頻率以避開發(fā)動機的激振頻率范圍,在多學科優(yōu)化設(shè)計ISIGHT軟件框架中采用NLPQL方法對支撐臂結(jié)構(gòu)進行了優(yōu)化。優(yōu)化后支撐臂側(cè)壁腔體厚度5.7 mm、橫梁寬度42.0 mm、臂長497.0 mm,其1階模態(tài)固有頻率135.17 Hz避開了發(fā)動機86.67~120 Hz激振頻率范圍,可在插秧機工作時有效避免共振的產(chǎn)生。

    [1] Xue Zhipeng, Li Ming, Jia Hongguang. Modal method for dynamics analysis of cantilever type structures at large rotational deformations[J]. International Journal of Mechanical Sciences, 2015(93): 22-31.

    [2] Morvan O, Emmanuel F. Model correlation and identification of experimental reduced models in vibroacoustical modal analysis[J]. Journal of Sound and Vibration, 2015(342): 200-217.

    [3] 高志朋,徐立章,李耀明,等. 履帶式稻麥聯(lián)合收獲機田間收獲工況下振動測試與分析[J]. 農(nóng)業(yè)工程學報,2017,33(20):48-55. Gao Zhipeng, Xu Lizhang, Li Yaoming, et al. Vibration measure and analys is of crawler-type rice and wheat combine harvester in field harvesting condition[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(20): 48 55. (in Chinese with English abstract)

    [4] 徐立章,李耀明,孫朋朋,等. 履帶式全喂入水稻聯(lián)合收獲機振動測試與分析[J]. 農(nóng)業(yè)工程學報,2014,30(8):49-55. Xu Lizhang, Li Yaoming, Sun Pengpeng, et al. Vibration measurement and analysis of tracked-whole feeding rice combine harvester[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(8): 49-55. (in Chinese with English abstract)

    [5] 邱白晶,何耀杰,盛云輝,等. 噴霧機噴桿有限元模態(tài)分析與結(jié)構(gòu)優(yōu)化[J]. 農(nóng)業(yè)機械學報,2014,45(8):112-116,105. Qiu Baijing, He Yaojie, Sheng Yunhu, et al.Finite element modal analysis and structure optimization of spray boom[J].Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(8): 112-116, 105.

    [6] 李發(fā)宗,童水光,王相兵. 基于模態(tài)分析的液壓挖掘機工作裝置動態(tài)優(yōu)化設(shè)計[J]. 農(nóng)業(yè)機械學報,2014,45(4):28-36. Li Fazong,Tong Shuiguang,Wang Xiangbing.Dynamic optimization design for working device of hydraulic excavator based on modal analysis[J].Transactions of the Chinese Society for Agricultural Machinery,2014,45 (4):28-36.

    [7] 陳樹人,韓紅陽,盧強. 4LZ—2.0型聯(lián)合收獲機割臺模態(tài)分析[J]. 農(nóng)業(yè)機械學報,2012,43(Z1):90-94. Chen Shuren, Han Hongyang, Lu Qiang.Modal analysis of header for type 4LZ—2.0 combine harvester[J]. Transacions of the Chinese Society for Agricultural Machinery, 2012, 43(1): 90-94.

    [8] 趙湛,李耀明,梁振偉,等. 基于振動模態(tài)分析的籽粒檢測傳感器結(jié)構(gòu)優(yōu)化設(shè)計[J]. 農(nóng)業(yè)機械學報,2011,42(Z1):103-106. Zhao Zhan, Li Yaoming, Liang Zhenwei, et al.Structure optimization of grain detecting sensor based on vibration modal analysis[J].Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(1): 103-106.

    [9] 李耀明,孫朋朋,龐靖,等. 聯(lián)合收獲機底盤機架有限元模態(tài)分析與試驗[J]. 農(nóng)業(yè)工程學報,2013,29(3):38-46. Li Yaoming, Sun Pengpeng, Pang Jing, et al. Finite element mode analysis and experiment of combine harvester chassis[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(3):38-46. (in Chinese with English abstract)

    [10] 李耀明,李有為,徐立章,等. 聯(lián)合收獲機割臺機架結(jié)構(gòu)參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學報,2014,30(18):30-37. Li Yaoming, Li Youwei, Xu Lizhang, et al. Structural parameter optimization of combine harvester cutting bench[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(18): 30-37. (in Chinese with English abstract)

    [11] 姚艷春,杜岳峰,朱忠祥,等. 基于模態(tài)的玉米收獲機車架振動特性分析與優(yōu)化[J]. 農(nóng)業(yè)工程學報,2015,31(19):46-53. Yao Yanchun, Du Yuefeng, Zhu Zhongxiang, et al. Vibration characteristics analysis and optimization of corn combine harvester frame using modal analysis method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(19): 46-53. (in Chinese with English abstract)

    [12] Wang J L, John L, Christof D, et al. Dealing with periodical loads and harmonics in operational modal analysis using time-varying transmissibility functions[J]. Mechanical Systems and Signal Processing, 2014(49): 154-164.

    [13] 張華標,朱梅,朱德泉,等. 高速水稻插秧機分插機構(gòu)振動特性研究[J]. 安徽農(nóng)業(yè)科學,2012,(18):9906-9907,9919.

    Zhang Huabiao, Zhu Mei, Zhu Dequan, et al. Study on Vibration Characteristics of High-speed Rice Transplanting Machine[J]. Anhui Agricultural Sciences, 2012, (18): 9906-9907, 9919.

    [14] 李復輝. 高速插秧機旋轉(zhuǎn)前插式分插機構(gòu)的研究[D]. 淄博:山東理工大學,2014.

    Li Fuhui.Study on the Rotating Front Insertion and Inserting Mechanism of High Speed Rice Transplanter[D].Zibo: Shandong University of Technology, 2014.

    [15] 張翔. LNC8型插秧機車架靈敏參數(shù)優(yōu)化設(shè)計及試驗驗證[D]. 杭州:浙江理工大學,2015.

    Zhang Xiang. Optimization Design and Experimental Verification of Sensitive Parameters of LNC8 Type Transplanter Frame [D]. Hangzhou: Zhejiang University of Science and Technology, 2015.

    [16] 趙曉丹,張中業(yè),駱英. 基于內(nèi)積運算與迭代算法的密集模態(tài)阻尼識別[J]. 農(nóng)業(yè)機械學報,2011,42(4):206-210. Zhao Xiaodan, Zhang Zhongye, Luo Ying.Damping identification for closely spaced modes based on Inner product calculation and iterative algorithm[J].Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(4): 206-210.

    [17] 姚艷春,宋正河,杜岳峰,等. 玉米收獲機車架應(yīng)力及模態(tài)數(shù)值模擬焊點模型優(yōu)選[J]. 農(nóng)業(yè)工程學報,2016,32(24):50-58.Yao Yanchun, Song Zhenghe, Du Yuefeng, et al. Optimum seeking of spot weld model on numerical simulation of stress and modal analysis for corn combine harvester frame[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(24): 50-58.

    [18] 姚艷春,趙雪彥,杜岳峰,等. 考慮質(zhì)量時變的收獲機械工作模態(tài)分析與試驗[J]. 農(nóng)業(yè)工程學報,2018,34(9):83-94.Yao Yanchun, Zhao Xueyan, Du Yuefeng, et al. Operating modal analysis and test of harvester induced by mass-varying process[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(9): 83-94.

    [19] 張建偉,曹克磊,趙瑜,等. 基于流固耦合模型的U型渡槽模態(tài)分析及驗證[J]. 農(nóng)業(yè)工程學報,2016,32(18):98-104. Zhang Jianwei, Cao Kelei, Zhao Yu, et al. Modal analysis and validation of U-shaped aqueducts based on fluid solid interaction model[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(18): 98-104.

    [20] 蔡力鋼,馬仕明,趙永勝,等. 多約束狀態(tài)下重載機械式主軸有限元建模及模態(tài)分析[J]. 機械工程學報,2012, 48(3):165-172. Cai Ligang, Ma Shiming, Zhao Yongsheng, et al. Finite element modeling and modal analysis of heavy-duty mechanical spindle under multiple constraints[J]. Journal of Mechanical Engineering, 2012, 48(3): 165-172. (in Chinese with English abstract)

    [21] Kang Tae Hwan, Kaizu Yutaka. Vibration analysis during grass harvesting according to ISO vibration standards[J]. Computers and Electronics in Agriculture, 2011, 79(2): 226-235.

    [22] Reza Ebrahimi, Mohsen Esfahanian, Saeed Ziaei-Rad. Vibration modeling and modification of cutting platform in a harvest combine by means of operational modal analysis(OMA)[J]. Measurement, 2013(46): 3959-3967.

    [23] Brincker Rune, Ventura Carlos E. Introduction to operational modal analysis[M]. New York:John Wiley & Sons, 2015.

    [24] Edwin Reynders. System identification methods for (operational) modal analysis: Review and comparison[J]. Archives of Computational Methods in Engineering, 2012, 19(1): 51-124.

    [25] 于天彪,王學智,關(guān)鵬,等. 超高速磨削機床主軸系統(tǒng)模態(tài)分析[J]. 機械工程學報,2012,48(17):183-188. Yu Tianbiao, Wang Xuezhi, Guan Peng, et al. Modal analysis of spindle system on ultra-high speed grinder[J]. Journal of Mechanical Engineering, 2012, 48(17): 183-188. (in Chinese with English abstract)

    [26] 戴光昊,付金波,張海福,等. 齒輪箱有限元模態(tài)分析及試驗研究[J]. 艦船科學技術(shù),2010 (8):167-170,213.

    Dai Guangwei,Fu Jinbo,Zhang Haifu, et al. Finite element modal analysis and experimental study of gearbox[J].Ship Science and Technology,2010(8):167-170,213.

    [27] 陳鋒,李忠獻. 梁橋上移動荷載識別中的測點優(yōu)化[J]. 振動、測試與診斷,2009,29(2):218-222. Chen feng, Li Zhongxian. Optimization of measurement points in mobile load identification on beam bridge [J]. Vibration, test and diagnosis, 2009,29(2): 218-222.

    [28] 王民,牛煥煥,高相勝,等. 電火花機床主軸頭的模態(tài)分析與減振設(shè)計[J]. 振動、測試與診斷,2017,(5):905-909. Wang Min, Niu Huanhuan, Gao Xiangsheng, et al. Modal analysis and vibration damping design of spindle head of edm [J]. Vibration, test and diagnosis, 2017,(5): 905-909.

    [29] 于存貴,卓艾寶,徐合剛. 復合材料定向管剛度優(yōu)化設(shè)計[J].南京理工大學學報(自然科學版),2012,36(6):1042-1046. Yucungui, Chao Aibao, Xu Hagang. Optimum design of the stiffness of the directional tube of composite materials [J]. Journal of Nanjing University of Science and Engineering (Natural Sciences Edition), 2012, 36(6): 1042-1046.

    [30] 左曙光,蔣維旭,吳旭東,等. 基于ISIGHT的電動汽車扭桿梁懸架振動傳遞特性優(yōu)化[J]. 吉林大學學報(工學版),2015,(5):1381-1387. Zuo Shuguang, Jiang Weixu, Wu Xudong, et al. Optimization of vibration transmission characteristics of electric car torsion beam suspension based on ISIGHT [J]. Journal of jilin university (Engineering Edition), 2015,(5): 1381-1387.

    [31] 張維峰,范傳帥. 車架結(jié)構(gòu)動力修改與力響應(yīng)模擬技術(shù)研究[J]. 力學與實踐,2012,34(4):19-23. Zhang Weifeng, Fan Chuanshuai. Study of dynamic modification and FRs technology in vehicle frame design[J]. Mechanics in Engineering, 2012,34(4): 19-23. (in Chinese with English abstract)

    [32] 廖宇蘭,劉世豪,孫佑攀,等. 基于靈敏度分析的木薯收獲機車架結(jié)構(gòu)優(yōu)化設(shè)計[J]. 農(nóng)業(yè)機械學報,2013,44(12):56-62. Liao Yulan, Liu Shihao, Sun Youpan, et al. Structural optimization for rack of cassava harvester based on sensitivity analysis[J]. Transactions of Chinese Society for Agricultural Machinery, 2013,44(12): 56-62. (in Chinese with English abstract)

    Modal analysis and structure optimization of transplanter support arm based on modal assurance criterion

    Jin Xin1,2, Chen Kaikang1, Ji Jiangtao1,2※, Pang Jing1, Gao Song1, Zeng Xinyue1

    (1.,,471003,; 2.,471003,

    In recent years, with the rapid development of agricultural mechanization in China, rice transplanter has been widely used. Because the transplanting parts of the transplanter have a large vibration in the working process, it will have a great impact on the operation of the machine and the operating environment of the driver. At present, the analysis of the vibration of the support arm of the transplanter and the research on how to reduce the vibration and noise are still in initial stage. In the working process, the support arm bears a variety of dynamic loads, including the working load of various parts, engine and road surface stimulation. When the frequency of the dynamic load is close to the natural frequency of the support arm, the resonance phenomenon will occur. The mechanical resonance affects the normal operation of the transplanter, reduces the service life of the transplater, and seriously affects the comfort of the working environment. Therefore, it is of great of significance to study the vibration and noise reduction of rice transplanter. In view of the vibration and noise problems in the power transmission process of high-speed rice transplanter planting mechanism, 2ZG-6DK rice transplanter was taken as the research object in this paper, and the power transmission mode and working principle of planting system were analyzed. The 3D software was used to standardize the support arm which is the key components in the power transmission system of the transplanter, and the model was imported into the ANSYS Workbench and combined with the Lanczos Method to solve the modal natural frequency and mode shape. In order to optimize the structure of the support arm based on the finite element model, it is necessary to verify the accuracy of the finite element model of the support arm. The modal test of the supporting arm of Nantong 2ZG-6DK transplanter was carried out by using DH5902 Dynamic Data acquisition and modal analysis system, the first 4-order modal frequency and mode of the support arm of the transplanter were obtained by signal analysis and parameter identification. The accuracy of the finite element model of the support arm was verified by the comparison between the modal test results and the finite element analysis results. Based on this, MAC-based (modal assurance) was carried out and the accuracy of finite element analysis was verified. The criterion test optimized the modal test to verify the accuracy of the finite element theory analysis. In order to avoid the external excitation frequency of the support arm, based on the analysis of the external frequency excitation characteristics, and on account of the ISIGHT multidisciplinary software platform, the structural parameters of the support arm were optimized by the sequential quadratic programming method. The results showed that the thickness of the optimized support arm sidewall cavity was 5.7 mm, the beam width was 42.0 mm, the arm length was 497.0 mm, and the first 4–order modal frequencies were adjusted to 135.17, 204.23, 483.14 and 702.32 Hz respectively, which all avoided the gasser engine excitation frequency range 86.67-120 Hz. After the optimization, the vibration amplitude attenuation was the most obvious at the 1-order frequency, the vibration amplitude was decreased by 9.4%, and the low-frequency vibration characteristics of the support arm were significantly improved. The research results can provide reference for the vibration characteristics analysis and vibration reduction design of the rice transplanter.

    agricultural machinery; vibration; optimization; modal analysis; MAC criterion

    10.11975/j.issn.1002-6819.2018.18.012

    S223.91

    A

    1002-6819(2018)-18-0093-09

    2018-05-21

    2018-08-02

    國家重點研發(fā)計劃項目(2017YFD0700801);國家自然科學基金項目(51505130);河南省科技創(chuàng)新杰出人才項目(184200510017)

    金 鑫,副教授,主要從事種苗高速栽插裝備與智能化技術(shù)研究。Email:jx.771@163.com

    姬江濤,教授,博士生導師,主要從事田間作物高效生產(chǎn)關(guān)鍵技術(shù)與裝備研發(fā)工作。Email:jjt0907@163.com

    金 鑫,陳凱康,姬江濤,龐 靖,高 頌,曾欣悅. 基于模態(tài)置信度準則的插秧機支撐臂模態(tài)分析與結(jié)構(gòu)優(yōu)化[J]. 農(nóng)業(yè)工程學報,2018,34(18):93-101. doi:10.11975/j.issn.1002-6819.2018.18.012 http://www.tcsae.org

    Jin Xin ,Chen Kaikang, Ji Jiangtao, Pang Jing, Gao Song, Zeng Xinyue. Modal analysis and structure optimization of transplanter support arm based on modal assurance criterion[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 93-101. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.18.012 http://www.tcsae.org

    猜你喜歡
    插秧機固有頻率測點
    液壓支架整機靜強度試驗及等效應(yīng)力分析
    水稻插秧機發(fā)展綜述
    現(xiàn)場測定大型水輪發(fā)電機組軸系的固有頻率
    基于CATIA的汽車測點批量開發(fā)的研究與應(yīng)用
    超智能插秧機
    水稻插秧機的插秧技術(shù)
    總溫總壓測頭模態(tài)振型變化規(guī)律研究
    拱壩結(jié)構(gòu)損傷的多測點R/S分析
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    轉(zhuǎn)向系統(tǒng)固有頻率設(shè)計研究
    亚洲国产精品sss在线观看| 99国产精品一区二区三区| 久久性视频一级片| 一级a爱片免费观看的视频| 日日夜夜操网爽| 男人的好看免费观看在线视频 | 亚洲国产精品999在线| 久久久久久久久中文| av天堂久久9| 国产亚洲欧美在线一区二区| 咕卡用的链子| 免费搜索国产男女视频| 大型黄色视频在线免费观看| 最新美女视频免费是黄的| 精品久久久久久久毛片微露脸| 搞女人的毛片| 给我免费播放毛片高清在线观看| 国产精品久久久av美女十八| 久久久久精品国产欧美久久久| 成人免费观看视频高清| 欧美日韩黄片免| 桃色一区二区三区在线观看| 国产精品九九99| 麻豆成人av在线观看| 精品熟女少妇八av免费久了| 国产蜜桃级精品一区二区三区| 日本vs欧美在线观看视频| 午夜激情av网站| 国产91精品成人一区二区三区| av超薄肉色丝袜交足视频| 在线观看日韩欧美| 最新美女视频免费是黄的| 亚洲 国产 在线| 国产精品九九99| 波多野结衣巨乳人妻| 久久欧美精品欧美久久欧美| а√天堂www在线а√下载| 侵犯人妻中文字幕一二三四区| 国产精品永久免费网站| 亚洲五月婷婷丁香| 一二三四社区在线视频社区8| 亚洲成人免费电影在线观看| 午夜久久久久精精品| www.精华液| 夜夜躁狠狠躁天天躁| 精品久久久久久久人妻蜜臀av | 黄色 视频免费看| 久久久久久人人人人人| 无遮挡黄片免费观看| 国产一卡二卡三卡精品| 亚洲无线在线观看| 69av精品久久久久久| 成年人黄色毛片网站| 国产成人啪精品午夜网站| 老司机深夜福利视频在线观看| 又黄又粗又硬又大视频| 亚洲国产精品sss在线观看| 国产精品av久久久久免费| 国产精品秋霞免费鲁丝片| 热re99久久国产66热| 女人被躁到高潮嗷嗷叫费观| 成人亚洲精品一区在线观看| 性欧美人与动物交配| 欧美乱码精品一区二区三区| 国产精品久久久av美女十八| 国产精品久久视频播放| 级片在线观看| 亚洲欧美精品综合一区二区三区| 亚洲精品久久成人aⅴ小说| 亚洲av第一区精品v没综合| 女人被狂操c到高潮| 亚洲国产日韩欧美精品在线观看 | 国产蜜桃级精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 精品无人区乱码1区二区| 啦啦啦韩国在线观看视频| 悠悠久久av| 男女床上黄色一级片免费看| bbb黄色大片| 真人做人爱边吃奶动态| 男女午夜视频在线观看| 黄频高清免费视频| 午夜成年电影在线免费观看| 亚洲 欧美 日韩 在线 免费| 国产精品日韩av在线免费观看 | 99re在线观看精品视频| 久久亚洲真实| 最新美女视频免费是黄的| 亚洲三区欧美一区| 88av欧美| 精品一区二区三区四区五区乱码| 国产免费男女视频| 波多野结衣高清无吗| 亚洲性夜色夜夜综合| 首页视频小说图片口味搜索| 757午夜福利合集在线观看| 国产xxxxx性猛交| av视频在线观看入口| 91麻豆av在线| 久久香蕉精品热| 欧美激情 高清一区二区三区| 在线观看免费日韩欧美大片| 91精品国产国语对白视频| 91字幕亚洲| 国产成人精品在线电影| 一卡2卡三卡四卡精品乱码亚洲| 久久久久国产精品人妻aⅴ院| 国产一区在线观看成人免费| 桃红色精品国产亚洲av| 真人做人爱边吃奶动态| 国产不卡一卡二| 日本a在线网址| 精品日产1卡2卡| 色老头精品视频在线观看| av欧美777| 久久天躁狠狠躁夜夜2o2o| 久久人妻av系列| 亚洲一区二区三区色噜噜| 19禁男女啪啪无遮挡网站| 精品一区二区三区四区五区乱码| 亚洲色图av天堂| 他把我摸到了高潮在线观看| 后天国语完整版免费观看| 国产1区2区3区精品| 男人舔女人的私密视频| 成在线人永久免费视频| 国产97色在线日韩免费| 亚洲中文日韩欧美视频| 精品久久久久久,| 大型av网站在线播放| 亚洲avbb在线观看| 黄色 视频免费看| 免费高清在线观看日韩| svipshipincom国产片| 在线观看www视频免费| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 国产精品国产高清国产av| 欧美老熟妇乱子伦牲交| 久久这里只有精品19| 无限看片的www在线观看| 国产不卡一卡二| 18禁裸乳无遮挡免费网站照片 | 男男h啪啪无遮挡| av中文乱码字幕在线| 日韩欧美免费精品| 黄色丝袜av网址大全| 精品国产亚洲在线| 亚洲熟妇熟女久久| 97超级碰碰碰精品色视频在线观看| 欧美日韩瑟瑟在线播放| 老汉色av国产亚洲站长工具| 两人在一起打扑克的视频| 美女国产高潮福利片在线看| 女生性感内裤真人,穿戴方法视频| 桃色一区二区三区在线观看| 国产精品98久久久久久宅男小说| 久久久久久久午夜电影| 中文字幕最新亚洲高清| 老熟妇仑乱视频hdxx| 久久人妻福利社区极品人妻图片| 成人18禁高潮啪啪吃奶动态图| 午夜精品国产一区二区电影| 一级片免费观看大全| 精品一区二区三区av网在线观看| 神马国产精品三级电影在线观看 | 村上凉子中文字幕在线| 女同久久另类99精品国产91| 国产高清videossex| 精品久久久久久久人妻蜜臀av | 此物有八面人人有两片| 国产不卡一卡二| 男女床上黄色一级片免费看| 亚洲国产欧美一区二区综合| 欧美黄色淫秽网站| 一区二区三区国产精品乱码| 色哟哟哟哟哟哟| 欧美黄色片欧美黄色片| 久久亚洲真实| www.精华液| 国产精品国产高清国产av| 一级a爱片免费观看的视频| 丝袜美腿诱惑在线| 久热这里只有精品99| 免费观看人在逋| 国产人伦9x9x在线观看| 国产在线观看jvid| 国产精品美女特级片免费视频播放器 | 此物有八面人人有两片| 一区二区三区精品91| 色综合婷婷激情| av在线播放免费不卡| 久久性视频一级片| 亚洲色图av天堂| 国产一区二区在线av高清观看| 精品午夜福利视频在线观看一区| 国产精品综合久久久久久久免费 | 亚洲国产日韩欧美精品在线观看 | 真人一进一出gif抽搐免费| 女人爽到高潮嗷嗷叫在线视频| 国产主播在线观看一区二区| 精品熟女少妇八av免费久了| 中文亚洲av片在线观看爽| 国产成人精品在线电影| 黄色a级毛片大全视频| 亚洲自拍偷在线| 桃红色精品国产亚洲av| 老熟妇乱子伦视频在线观看| 免费看美女性在线毛片视频| 久久国产亚洲av麻豆专区| 在线观看一区二区三区| 操美女的视频在线观看| 亚洲欧美日韩另类电影网站| 成在线人永久免费视频| 成人欧美大片| 18禁国产床啪视频网站| 自拍欧美九色日韩亚洲蝌蚪91| 又紧又爽又黄一区二区| 国产精品,欧美在线| 国产成人精品久久二区二区91| 国内精品久久久久久久电影| 久久精品亚洲精品国产色婷小说| 美女午夜性视频免费| 丝袜人妻中文字幕| 国产在线精品亚洲第一网站| 亚洲中文字幕一区二区三区有码在线看 | 国产伦一二天堂av在线观看| 婷婷精品国产亚洲av在线| 午夜免费观看网址| 黄频高清免费视频| 色综合欧美亚洲国产小说| 国产精品乱码一区二三区的特点 | 亚洲第一欧美日韩一区二区三区| 欧美乱色亚洲激情| 国产亚洲精品第一综合不卡| 国产伦人伦偷精品视频| 大码成人一级视频| 国产精品99久久99久久久不卡| 1024香蕉在线观看| 午夜精品在线福利| 日日爽夜夜爽网站| 成人国产一区最新在线观看| 国内毛片毛片毛片毛片毛片| 中文字幕色久视频| 老熟妇乱子伦视频在线观看| 一a级毛片在线观看| 一二三四社区在线视频社区8| 国产色视频综合| 免费在线观看黄色视频的| 国产成人av激情在线播放| 中文亚洲av片在线观看爽| 男人舔女人的私密视频| 人人妻,人人澡人人爽秒播| 精品久久久久久久毛片微露脸| 亚洲熟妇熟女久久| 欧美在线一区亚洲| 亚洲第一电影网av| 欧美日韩瑟瑟在线播放| 国产又爽黄色视频| 色综合站精品国产| 久久久久久久精品吃奶| 亚洲人成电影观看| 日本免费a在线| svipshipincom国产片| 悠悠久久av| 亚洲av成人av| 欧美日韩一级在线毛片| 中文字幕最新亚洲高清| 亚洲午夜精品一区,二区,三区| 成人三级黄色视频| 欧美成人性av电影在线观看| 免费久久久久久久精品成人欧美视频| 曰老女人黄片| 美女高潮喷水抽搐中文字幕| 99在线视频只有这里精品首页| 九色国产91popny在线| 成年版毛片免费区| www.999成人在线观看| 国产精品 国内视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区中文字幕在线| 免费av毛片视频| 欧美色视频一区免费| 丝袜美足系列| 人人妻人人澡欧美一区二区 | 免费一级毛片在线播放高清视频 | 亚洲免费av在线视频| 中文字幕精品免费在线观看视频| 精品国产一区二区三区四区第35| 一级片免费观看大全| 久久人人爽av亚洲精品天堂| 999久久久国产精品视频| 18禁美女被吸乳视频| 亚洲国产精品久久男人天堂| 国产成人欧美| 俄罗斯特黄特色一大片| 精品欧美一区二区三区在线| 波多野结衣巨乳人妻| 夜夜夜夜夜久久久久| 不卡av一区二区三区| 国产私拍福利视频在线观看| 大香蕉久久成人网| 国内久久婷婷六月综合欲色啪| 一进一出抽搐gif免费好疼| 国产主播在线观看一区二区| 麻豆av在线久日| 高清在线国产一区| 精品少妇一区二区三区视频日本电影| 黄片小视频在线播放| 亚洲人成电影免费在线| 亚洲久久久国产精品| 美女 人体艺术 gogo| 又黄又粗又硬又大视频| 亚洲,欧美精品.| 三级毛片av免费| 高清在线国产一区| 国内久久婷婷六月综合欲色啪| 成年女人毛片免费观看观看9| 色哟哟哟哟哟哟| 啪啪无遮挡十八禁网站| 窝窝影院91人妻| av欧美777| 亚洲免费av在线视频| 国产99白浆流出| 香蕉丝袜av| 久久人妻av系列| 韩国精品一区二区三区| 老鸭窝网址在线观看| 亚洲 欧美一区二区三区| 日韩 欧美 亚洲 中文字幕| 女警被强在线播放| 手机成人av网站| 大陆偷拍与自拍| 欧美最黄视频在线播放免费| 国产主播在线观看一区二区| 国产国语露脸激情在线看| 亚洲色图av天堂| videosex国产| 黄色丝袜av网址大全| 人人妻人人澡欧美一区二区 | 欧美国产精品va在线观看不卡| 国产成年人精品一区二区| 日韩一卡2卡3卡4卡2021年| 日韩中文字幕欧美一区二区| 亚洲精品粉嫩美女一区| av片东京热男人的天堂| 免费无遮挡裸体视频| 亚洲熟妇熟女久久| 久久久久国产一级毛片高清牌| 又紧又爽又黄一区二区| 久久久国产精品麻豆| av天堂久久9| 免费看a级黄色片| 高清毛片免费观看视频网站| 精品卡一卡二卡四卡免费| 午夜福利一区二区在线看| 可以在线观看的亚洲视频| 美女国产高潮福利片在线看| 亚洲精品粉嫩美女一区| 色av中文字幕| 悠悠久久av| 欧美在线一区亚洲| 免费一级毛片在线播放高清视频 | 欧美色欧美亚洲另类二区 | 成人国语在线视频| 欧美av亚洲av综合av国产av| 人妻久久中文字幕网| 精品国产一区二区三区四区第35| 欧洲精品卡2卡3卡4卡5卡区| 在线永久观看黄色视频| 两个人看的免费小视频| 亚洲久久久国产精品| 夜夜躁狠狠躁天天躁| 人人妻人人爽人人添夜夜欢视频| 国产99久久九九免费精品| 日本精品一区二区三区蜜桃| 日韩大尺度精品在线看网址 | 两性午夜刺激爽爽歪歪视频在线观看 | 日韩免费av在线播放| 欧美日韩一级在线毛片| 人人妻,人人澡人人爽秒播| 日本vs欧美在线观看视频| 人人妻,人人澡人人爽秒播| 老司机午夜十八禁免费视频| 91成人精品电影| 国产区一区二久久| 亚洲国产毛片av蜜桃av| 欧美国产日韩亚洲一区| av在线播放免费不卡| 国产激情久久老熟女| 国产高清激情床上av| 午夜久久久在线观看| 久久亚洲精品不卡| 欧美大码av| 国产91精品成人一区二区三区| 最近最新中文字幕大全电影3 | 亚洲av日韩精品久久久久久密| 看黄色毛片网站| 国产精品爽爽va在线观看网站 | 少妇 在线观看| 一级,二级,三级黄色视频| 一二三四在线观看免费中文在| 在线观看午夜福利视频| 色哟哟哟哟哟哟| 久99久视频精品免费| 日韩精品青青久久久久久| 美女 人体艺术 gogo| 最好的美女福利视频网| 亚洲av电影不卡..在线观看| 18禁国产床啪视频网站| 一级毛片精品| 香蕉丝袜av| 女性被躁到高潮视频| 满18在线观看网站| 久久香蕉精品热| 亚洲色图 男人天堂 中文字幕| 久久婷婷成人综合色麻豆| 操出白浆在线播放| 亚洲美女黄片视频| 法律面前人人平等表现在哪些方面| 黄色女人牲交| 50天的宝宝边吃奶边哭怎么回事| 女性生殖器流出的白浆| 国产伦一二天堂av在线观看| 母亲3免费完整高清在线观看| 国产乱人伦免费视频| 亚洲人成电影观看| 日本免费一区二区三区高清不卡 | 日韩精品免费视频一区二区三区| av天堂久久9| 亚洲精品美女久久久久99蜜臀| 91精品国产国语对白视频| 大陆偷拍与自拍| 亚洲最大成人中文| 国产成人影院久久av| 999精品在线视频| 人人妻人人澡人人看| 91成人精品电影| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美一区二区三区在线观看| 日本在线视频免费播放| 日韩欧美一区二区三区在线观看| 中文字幕久久专区| 亚洲九九香蕉| 亚洲av五月六月丁香网| 十八禁人妻一区二区| 12—13女人毛片做爰片一| 精品一区二区三区四区五区乱码| 国产真人三级小视频在线观看| 成人三级做爰电影| 久久精品国产清高在天天线| 黄色毛片三级朝国网站| 51午夜福利影视在线观看| 少妇 在线观看| 99在线视频只有这里精品首页| 91字幕亚洲| 国产单亲对白刺激| 91麻豆精品激情在线观看国产| 欧美最黄视频在线播放免费| 制服丝袜大香蕉在线| 曰老女人黄片| 精品欧美国产一区二区三| 久久久久亚洲av毛片大全| 国产单亲对白刺激| 日韩欧美一区二区三区在线观看| 可以在线观看毛片的网站| 亚洲狠狠婷婷综合久久图片| 又黄又爽又免费观看的视频| 禁无遮挡网站| 女人高潮潮喷娇喘18禁视频| 欧美一级毛片孕妇| 久久天堂一区二区三区四区| 精品久久久久久久人妻蜜臀av | 在线国产一区二区在线| 香蕉丝袜av| 国产精品98久久久久久宅男小说| 99riav亚洲国产免费| ponron亚洲| 最新在线观看一区二区三区| 日韩欧美一区视频在线观看| 国产激情欧美一区二区| 亚洲欧美一区二区三区黑人| 亚洲人成网站在线播放欧美日韩| 热99re8久久精品国产| 国产xxxxx性猛交| 久久人妻av系列| 国产一区二区激情短视频| 午夜福利影视在线免费观看| 免费看十八禁软件| 韩国精品一区二区三区| 国产精品久久电影中文字幕| 日韩三级视频一区二区三区| 少妇熟女aⅴ在线视频| 男人的好看免费观看在线视频 | 老汉色∧v一级毛片| 两人在一起打扑克的视频| 久久精品国产亚洲av香蕉五月| 精品国内亚洲2022精品成人| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费视频日本深夜| 午夜老司机福利片| 日韩欧美一区视频在线观看| 国产成人免费无遮挡视频| 国产高清有码在线观看视频 | 怎么达到女性高潮| 18禁美女被吸乳视频| 一区二区三区激情视频| aaaaa片日本免费| 日本 av在线| 亚洲色图综合在线观看| 国产成人精品久久二区二区91| 999久久久精品免费观看国产| 极品人妻少妇av视频| av欧美777| 亚洲 欧美一区二区三区| 日本a在线网址| 欧美日韩黄片免| 久久久国产成人精品二区| 国产一卡二卡三卡精品| 在线观看舔阴道视频| 久久精品成人免费网站| 午夜免费成人在线视频| av天堂久久9| 日韩视频一区二区在线观看| 在线观看免费视频日本深夜| 性少妇av在线| 在线永久观看黄色视频| 国产亚洲精品av在线| 在线观看午夜福利视频| av中文乱码字幕在线| 国产片内射在线| 亚洲av电影在线进入| 亚洲专区国产一区二区| 老司机午夜十八禁免费视频| 午夜免费观看网址| 香蕉久久夜色| 两性夫妻黄色片| 国产野战对白在线观看| 亚洲在线自拍视频| 少妇裸体淫交视频免费看高清 | 国产99久久九九免费精品| 日日干狠狠操夜夜爽| 麻豆久久精品国产亚洲av| 国产成人欧美在线观看| 亚洲久久久国产精品| 天天一区二区日本电影三级 | 搡老岳熟女国产| 怎么达到女性高潮| cao死你这个sao货| 欧美+亚洲+日韩+国产| 99国产精品一区二区蜜桃av| 丝袜在线中文字幕| 成人亚洲精品av一区二区| 国产午夜精品久久久久久| 亚洲午夜精品一区,二区,三区| 欧美日韩福利视频一区二区| 亚洲欧美日韩无卡精品| 免费在线观看日本一区| 最近最新中文字幕大全免费视频| 国产高清视频在线播放一区| 国产激情久久老熟女| 亚洲无线在线观看| 亚洲人成网站在线播放欧美日韩| 成人国语在线视频| АⅤ资源中文在线天堂| 777久久人妻少妇嫩草av网站| 精品国产一区二区久久| 亚洲成人国产一区在线观看| 女人爽到高潮嗷嗷叫在线视频| 青草久久国产| 中文字幕人成人乱码亚洲影| 国产激情欧美一区二区| 一级,二级,三级黄色视频| 女人高潮潮喷娇喘18禁视频| 精品久久久久久久久久免费视频| 最近最新中文字幕大全免费视频| 男人舔女人的私密视频| 午夜久久久久精精品| 91av网站免费观看| 女人被躁到高潮嗷嗷叫费观| 日本在线视频免费播放| 叶爱在线成人免费视频播放| 亚洲一区二区三区色噜噜| 精品久久久久久成人av| 久久久水蜜桃国产精品网| 亚洲国产精品合色在线| 久99久视频精品免费| 免费av毛片视频| 无限看片的www在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻在线不人妻| 女人爽到高潮嗷嗷叫在线视频| 制服诱惑二区| 久久精品aⅴ一区二区三区四区| 在线观看舔阴道视频| 好男人在线观看高清免费视频 | 国产精品久久视频播放| 亚洲自拍偷在线| 中文字幕av电影在线播放| 久久精品成人免费网站| 国产在线观看jvid| 精品久久久久久久人妻蜜臀av | 欧美不卡视频在线免费观看 | 国产欧美日韩一区二区精品| 国产高清激情床上av| 美国免费a级毛片| 久久久国产成人免费| 色综合亚洲欧美另类图片| 亚洲自偷自拍图片 自拍| 国产精品美女特级片免费视频播放器 | av超薄肉色丝袜交足视频| 亚洲精品国产区一区二| 一区二区三区精品91| av天堂久久9| 88av欧美|