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

    寬型水田輪作業(yè)過(guò)程動(dòng)力學(xué)仿真

    2016-03-23 04:34:09楊紅艷
    農(nóng)機(jī)化研究 2016年4期
    關(guān)鍵詞:機(jī)理

    楊紅艷,楊 望,楊 堅(jiān),路 宏

    (廣西大學(xué) 機(jī)械工程學(xué)院,南寧 530004)

    ?

    寬型水田輪作業(yè)過(guò)程動(dòng)力學(xué)仿真

    楊紅艷,楊望,楊堅(jiān),路宏

    (廣西大學(xué) 機(jī)械工程學(xué)院,南寧530004)

    摘要:寬型水田輪是前耕后驅(qū)微耕機(jī)的行走驅(qū)動(dòng)裝置,除具有驅(qū)動(dòng)性能之外,還具有碾壓碎土和攪動(dòng)混土的作用,其設(shè)計(jì)參數(shù)對(duì)作業(yè)質(zhì)量有大的影響。為此,以1WGQ4型微耕機(jī)為對(duì)象,采用有限元法和光滑粒子流體動(dòng)力學(xué)方法(Smoothed particle hydrodynamics, SPH)相結(jié)合的方法,構(gòu)建土壤-寬輪系統(tǒng)的動(dòng)力學(xué)仿真模型,對(duì)其作業(yè)過(guò)程進(jìn)行動(dòng)力學(xué)仿真分析,在細(xì)微上研究其與土壤的作用機(jī)理。結(jié)果表明:隨著寬輪的轉(zhuǎn)動(dòng),輪葉與土壤的接觸應(yīng)力和輪葉水平推進(jìn)力均呈現(xiàn)先增大后減小的變化趨勢(shì),但接觸應(yīng)力比水平推進(jìn)力先達(dá)到最大值;滑轉(zhuǎn)率不同,輪葉出現(xiàn)水平推力為零的位置不同;輪葉折角對(duì)驅(qū)動(dòng)性能、碾壓碎土和攪動(dòng)混土作用有較大的影響。

    關(guān)鍵詞:寬型水田輪;作業(yè)過(guò)程;動(dòng)力學(xué)仿真;機(jī)理

    0引言

    南方水田廣泛應(yīng)用微型水田耕整機(jī)械(簡(jiǎn)稱微耕機(jī))進(jìn)行耕作,其水田驅(qū)動(dòng)輪(簡(jiǎn)稱水田輪)的結(jié)構(gòu)參數(shù)和運(yùn)動(dòng)參數(shù)對(duì)微耕機(jī)的驅(qū)動(dòng)性能影響較大。因此,國(guó)內(nèi)外學(xué)者通過(guò)物理試驗(yàn)和理論分析對(duì)水田輪的作業(yè)性能影響參數(shù)進(jìn)行了較多的研究,且取得了一定的成果[1-7]。

    1WGQ4型微耕機(jī)是一種前耕后驅(qū)的新型微耕機(jī),其水田輪為寬型水田輪(簡(jiǎn)稱寬輪),寬輪的寬度和旋耕刀具寬度相同;除具有驅(qū)動(dòng)性能之外,還具有碾壓碎土和攪動(dòng)混土作用,且其設(shè)計(jì)參數(shù)對(duì)作業(yè)質(zhì)量影響較大;同時(shí),其設(shè)計(jì)方法和一般的微耕機(jī)水田輪有較大的差別。目前,前耕后驅(qū)微耕機(jī)寬輪的設(shè)計(jì)方法未見有報(bào)導(dǎo),因此開展寬型水田輪設(shè)計(jì)方法研究具有重要意義。本文以1WGQ4型微耕機(jī)的寬輪為對(duì)象,采用有限元法和光滑粒子流體動(dòng)力學(xué)方法(SPH)相結(jié)合的方法[8-9],構(gòu)建土壤-寬輪系統(tǒng)的動(dòng)力學(xué)仿真模型,對(duì)其作業(yè)過(guò)程進(jìn)行動(dòng)力學(xué)分析,在微觀上研究其與土壤的作用機(jī)理,為寬輪的優(yōu)化設(shè)計(jì)提供依據(jù)。

    1結(jié)構(gòu)和工作原理

    1WGQ4型微耕機(jī)主要由機(jī)架、發(fā)動(dòng)機(jī)、減速箱、寬輪及旋耕刀具等組成,如圖1所示。

    圖1 1WGQ4型微耕機(jī)

    微耕機(jī)作業(yè)時(shí),發(fā)動(dòng)機(jī)輸出的動(dòng)力經(jīng)減速箱和傳動(dòng)機(jī)構(gòu)傳輸?shù)綄捿喓托毒?,在寬輪的?qū)動(dòng)作用下機(jī)組前進(jìn),微耕機(jī)前面的旋耕刀旋轉(zhuǎn)切削土壤,且把切下的土塊拋向擋泥板進(jìn)行破碎。機(jī)組后面的寬輪除了具有驅(qū)動(dòng)微耕機(jī)前進(jìn)之外,還具有碾壓大塊土壤使其進(jìn)一步破碎和攪動(dòng)混土作用。

    2動(dòng)力學(xué)仿真模型

    2.1模型的建立

    1WGQ4型微耕機(jī)寬輪由連接圓板、輻條、輪圈及輪葉組成,輪葉為近似長(zhǎng)方體形,如圖2所示。

    由于寬輪比土壤的剛度大得多,為了節(jié)省仿真計(jì)算時(shí)間,建模時(shí)將其整體建成剛體。同時(shí),為了便于建模,將寬輪做以下簡(jiǎn)化:

    1)省去連接圓板和輻條,僅保留輪圈和輪葉;

    2)輪葉簡(jiǎn)化成長(zhǎng)方體。

    由于寬輪的對(duì)稱性,只建其1/2模型,簡(jiǎn)化后寬輪的模型如圖3所示。模型尺寸根據(jù)樣機(jī)的實(shí)際尺寸確定:輪圈外直徑450mm、內(nèi)直徑420mm、厚5mm;輪葉長(zhǎng)120mm、寬72mm、厚6mm;輪葉傾角γ為45°。

    圖2 寬輪

    圖3 寬輪模型

    由于水田的耕作層土壤(泥腳層)和硬底層土壤的硬度及物理特性參數(shù)差別較大,因此建模時(shí)土壤分上下兩層建模,形狀為長(zhǎng)方體:上層土壤為泥腳層,其尺寸(長(zhǎng)×寬×高)為1 550mm×195mm×100mm;下層土壤為硬底層,其尺寸(長(zhǎng)×寬×高)為1 550mm×195mm×40mm。土壤外圍設(shè)有剛性墻,約束土壤。建立的土壤-寬輪系統(tǒng)模型如圖4所示。

    圖4 土壤-寬輪系統(tǒng)模型

    2.2模型材料

    寬輪選用SOLID164實(shí)體單元類型建模,定義為剛性材料,材料參數(shù):密度為1.112×104kg/m3(因建模時(shí)省去了輻條和連接圓板,故適度增大材料密度,以保持其原有的質(zhì)量),泊松比為0.27,彈性模量為0.2×1012Pa。

    土壤材料采用LS-DYNA971中MAT_FHWA_SOIL材料模型,其針對(duì)實(shí)體單元且考慮了含水率、應(yīng)變軟化、應(yīng)變率效應(yīng)、孔隙比及孔隙水壓力等的影響和單元?jiǎng)h除[10]。為了解決土壤產(chǎn)生大變形時(shí)引起網(wǎng)格畸變,造成數(shù)值模擬計(jì)算失效的問(wèn)題,土壤使用SPH單元進(jìn)行建模。土壤的具體參數(shù)如表1所示。

    表1 土壤參數(shù)

    2.3網(wǎng)格劃分和加載

    仿真模型網(wǎng)格劃分直接影響到模型的計(jì)算精度和運(yùn)算時(shí)間。網(wǎng)格劃分細(xì),計(jì)算精度高,但運(yùn)算時(shí)間長(zhǎng);網(wǎng)格劃分過(guò)粗,計(jì)算精度低,易發(fā)生穿透。故在劃分網(wǎng)格時(shí),要合理匹配寬輪和土壤的網(wǎng)格尺寸,使其接觸良好,相互作用過(guò)程不產(chǎn)生穿透現(xiàn)象,保證模型運(yùn)算的精度。

    Smart Size是ANSYS提供的一種強(qiáng)大的自動(dòng)劃分網(wǎng)格算法,有利于劃分網(wǎng)格時(shí)生成合理的網(wǎng)格單元[11]。因此,本文采用Smart Size算法對(duì)寬輪進(jìn)行自由網(wǎng)格劃分,寬輪的單元數(shù)為12 010。土壤模型為長(zhǎng)方形,形狀規(guī)則,故土壤選擇映射式網(wǎng)格劃分方式進(jìn)行網(wǎng)格劃分:上層土壤的單元數(shù)為33 033,下層土壤的單元數(shù)為15 015。

    1WGQ4型微耕機(jī)總質(zhì)量為157kg,寬輪的質(zhì)量為14kg。正常作業(yè)時(shí),人通過(guò)扶手把施加在微耕機(jī)上的豎直向下的力約為210N,水平向后的力為140N。由于微耕機(jī)左右兩個(gè)寬輪受力基本一樣,且模型只建立了1/2寬輪,則1/2寬輪承受1/4的外力。因此,建模時(shí),在寬輪的質(zhì)量中心豎直向下施加410N的力,向后施加35N的力水平力。

    由于正常作業(yè)條件下寬輪的轉(zhuǎn)速為1.8rad/s,故建模時(shí)給寬輪施加1.8rad/s的轉(zhuǎn)速。土壤和寬輪之間定義為面面自動(dòng)接觸,靜摩擦因數(shù)為0.3,動(dòng)摩擦因數(shù)為0.25[12]。

    土壤在ANSYS/LS-DYNA中建模后,生成K文件,導(dǎo)入Ls-Prepost中,將土壤的有限元網(wǎng)格轉(zhuǎn)化為SPH粒子,SPH粒子的疏密程度取決于有限元網(wǎng)格劃分的大小。建立的土壤-寬輪系統(tǒng)動(dòng)力學(xué)仿真模型如圖5所示。

    圖5 土壤-寬輪系統(tǒng)動(dòng)力學(xué)仿真模型

    3模型的驗(yàn)證

    3.1試驗(yàn)方法及設(shè)備

    本文采用物理試驗(yàn)的方法對(duì)模型精度進(jìn)行驗(yàn)證。試驗(yàn)時(shí),測(cè)量空載條件下(微耕機(jī)的旋耕刀具離地不工作)1WGQ4型微耕機(jī)寬輪的下陷深度及滑轉(zhuǎn)率后,與相同條件下的仿真試驗(yàn)結(jié)果比較驗(yàn)證模型的精度。試驗(yàn)設(shè)備:1WGQ4型微耕機(jī)(水田刀配置,傳動(dòng)比為1檔38.5,2檔23,倒檔50),紅外轉(zhuǎn)速儀,卷尺,標(biāo)桿。試驗(yàn)地:廣西大學(xué)農(nóng)場(chǎng)水稻田。

    驗(yàn)證試驗(yàn)時(shí),選微耕機(jī)1檔速度行走,用直尺測(cè)量寬輪的下陷深度;采用紅外轉(zhuǎn)速儀測(cè)量與減速箱輸入軸相連的皮帶輪的轉(zhuǎn)速,后由減速箱傳動(dòng)比(一檔為38.5)計(jì)算出寬輪的實(shí)際轉(zhuǎn)速。寬輪的滑轉(zhuǎn)率采用傳統(tǒng)方式測(cè)量,即微耕機(jī)行走一定距離L1時(shí),數(shù)出寬輪轉(zhuǎn)的圈數(shù),后由寬輪轉(zhuǎn)的圈數(shù)和半徑計(jì)算出理論上微耕機(jī)行走的距離L,則滑轉(zhuǎn)率為

    (1)

    物理試驗(yàn)重復(fù)6次。

    3.2驗(yàn)證結(jié)果及分析

    寬輪物理試驗(yàn)的平均下陷深度約為8cm,滑轉(zhuǎn)率的試驗(yàn)結(jié)果如表2所示。

    圖6是仿真試驗(yàn)時(shí)寬輪的下陷深度曲線。由圖6可知:寬輪運(yùn)行穩(wěn)定后,下陷深度在8cm附近波動(dòng),即平均值約為8cm,和物理試驗(yàn)的結(jié)果一致。由表2可知:物理試驗(yàn)的平均滑轉(zhuǎn)率為3.18%,仿真試驗(yàn)的滑轉(zhuǎn)率為3.10%,平均滑轉(zhuǎn)率相對(duì)誤差為2.5%。這表明,建立的土壤-寬輪系統(tǒng)動(dòng)力學(xué)仿真模型精度高。

    表2  滑轉(zhuǎn)率試驗(yàn)結(jié)果

    圖6 下陷深度曲線圖

    4作業(yè)過(guò)程分析

    4.1單輪葉

    圖7是在寬輪和土壤作用過(guò)程中單一輪葉與土壤作用過(guò)程的動(dòng)力學(xué)仿真截圖(滑轉(zhuǎn)率3.10%),圖8是單一輪葉的水平推力變化圖。由圖7可知:隨著輪葉平面向下移動(dòng),輪葉對(duì)土壤的壓力不斷加強(qiáng),輪葉與土壤的接觸應(yīng)力逐漸增大。當(dāng)t=0.34s時(shí),接觸應(yīng)力最大;當(dāng)t>0.34s時(shí),雖然輪葉的前端還繼續(xù)下行,但輪葉的后端已開始上行,輪葉對(duì)土壤的壓力開始減小,輪葉與土壤的接觸應(yīng)力逐漸減??;當(dāng)t>0.76s時(shí),輪葉雖然已不向下后方壓土,但輪葉開始前推土壤,故輪葉與土壤間還存在接觸應(yīng)力。由圖8可知:t在0~0.22s和0.76~0.94s區(qū)間時(shí),輪葉水平推力為負(fù)值;t在0.22~0.76s區(qū)間時(shí),為正值,且隨時(shí)間的增加呈先增大后減小的變化趨勢(shì)。其原因是:t在0~0.22s區(qū)間時(shí),雖然輪葉平面向下移動(dòng),輪葉對(duì)土壤的壓力不斷加強(qiáng),輪葉與土壤的接觸應(yīng)力逐漸增大;但輪葉傾角較小,且輪葉水平方向的運(yùn)動(dòng)速度較大,輪葉厚度方向產(chǎn)生向前推土作用,故輪葉的水平推力為負(fù)值。當(dāng)t在0.76~0.94s區(qū)間時(shí),輪葉葉面已向前運(yùn)動(dòng),葉面對(duì)土壤產(chǎn)生向前的推壓作用,故輪葉的水平推力為負(fù)值;當(dāng)t在0.22~0.5s區(qū)間時(shí),輪葉向下后方運(yùn)動(dòng),向后對(duì)土壤的壓力不斷增大,且輪葉的傾角不斷增大,故輪葉的推力不斷增大,且為正值;當(dāng)t在0.5~0.76s區(qū)間時(shí),開始雖然輪葉前端還繼續(xù)向下后方運(yùn)動(dòng),輪葉傾角也在增大,但輪葉后端已上移較多,故輪葉對(duì)土壤的向后壓力逐漸減小,輪葉推力逐漸減小。同時(shí),由圖7和圖8可知:土壤應(yīng)力(0.34s時(shí)最大)比輪葉水平推力(0.5s時(shí)最大)先達(dá)到最大值。當(dāng)t=0.76s時(shí),即當(dāng)輪葉轉(zhuǎn)到最低點(diǎn)再轉(zhuǎn)17.5°時(shí),輪葉的水平推力為零,比文獻(xiàn)[10]研究結(jié)果的大4.5°。不同滑轉(zhuǎn)率條件下的仿真分析表明:滑轉(zhuǎn)率不同,出現(xiàn)水平推力為零的位置也不同,故出現(xiàn)與文獻(xiàn)研究結(jié)果不一致的情況。

    圖7 單一輪葉與土壤作用過(guò)程的動(dòng)力學(xué)仿真截圖

    圖8 單一輪葉的水平推力變化圖

    4.2寬輪

    寬輪除具有驅(qū)動(dòng)微耕機(jī)前進(jìn)的作用之外,還具有碾壓碎土和攪動(dòng)混土作用。因此,為了在寬輪作業(yè)過(guò)程的動(dòng)力學(xué)仿真圖中便于觀察寬輪的破碎土壤和攪動(dòng)混土作用效果,把耕作層土壤分為兩層(土壤材料參數(shù)不變)建模,使其產(chǎn)生不同顏色層(兩層土壤粒子相互混合多,表示攪動(dòng)混土能力強(qiáng))。圖9是寬輪對(duì)土壤的碾壓過(guò)程仿真截圖。其中,圖9(a)是輪葉為直板型寬輪對(duì)土壤碾壓過(guò)程的仿真截圖,圖9(b)是輪葉為折角型寬輪對(duì)土壤碾壓過(guò)程的仿真截圖。

    圖9 寬輪碾壓痕跡

    由圖9(a)、(b)可知:輪葉帶折角的寬輪比輪葉為直板型的土壤碾壓破碎作用強(qiáng),且耕作層的兩層土壤相互混合的粒子數(shù)也比直板型的多,表明其攪動(dòng)混土作用也強(qiáng)。兩種寬輪動(dòng)力學(xué)仿真表明:直板型輪葉寬輪的平均滑轉(zhuǎn)率為3.28%,牽引功率為650.239W,而折角型輪葉寬輪的平均滑轉(zhuǎn)率為5.97%,牽引功率為632.148W;直板型的打滑率小,牽引效率高,驅(qū)動(dòng)性能好,表明折角因素對(duì)寬輪的驅(qū)動(dòng)性能、碾壓碎土和攪動(dòng)混土作用影響較大。因此,為了使寬輪具有良好的驅(qū)動(dòng)性能,且具有較強(qiáng)的碾壓碎土和攪動(dòng)混土作用,其輪葉折角的角度應(yīng)進(jìn)行優(yōu)化設(shè)計(jì)。

    5結(jié)論

    1)采用有限元法和光滑粒子流體動(dòng)力學(xué)方法 (SPH) 相結(jié)合的方法構(gòu)建的土壤-寬輪系統(tǒng)的動(dòng)力學(xué)仿真模型精度高,可用于寬輪作業(yè)過(guò)程分析。

    2)隨著寬輪的轉(zhuǎn)動(dòng),輪葉與土壤的接觸應(yīng)力和輪葉水平推進(jìn)力均呈現(xiàn)先增大后減小的變化趨勢(shì),但接觸應(yīng)力比水平推進(jìn)力先達(dá)到最大值;且滑轉(zhuǎn)率不同,輪葉出現(xiàn)水平推力為零的位置不同;輪葉轉(zhuǎn)到最低點(diǎn)、再轉(zhuǎn)17.5°角時(shí),輪葉的水平推力為零,且滑轉(zhuǎn)率不同時(shí),出現(xiàn)輪葉水平推力為零的位置也不同。

    3)輪葉折角因素對(duì)寬輪的驅(qū)動(dòng)性能、碾壓碎土和攪動(dòng)混土作用影響較大,為使寬輪具有良好的驅(qū)動(dòng)性能,且具有較強(qiáng)的碾壓破土和攪動(dòng)混土作用,輪葉折角的角度應(yīng)進(jìn)行優(yōu)化設(shè)計(jì)。

    參考文獻(xiàn):

    [1]Arthur L. Fajardo, Delfin C. Suministrado, Engelbert K. Peralta, et al. Force and puddling characteristics of the tilling wheel of float-assisted tillers at different lug angle and shaft speed[J]. Soil and Tillage Research,2014,140:118-125.

    [2]Wawan Hermawan, Minoru Yamazaki, Akira Oida. Theoretical analysis of soil reaction on a lug of the movable lug cage wheel[J]. Journal of Terramechanics, 2000, 37(2) : 65-86.

    [3]C Watyotha, D Gee-Clough, V.M Salokhe. Effect of circumferential angle, lug spacing and slip on lug wheel forces[J].Journal of Terramechanics, 2001,38(1) :1-14.

    [4]C Watyotha, V.M Salokhe. Development of a data acquisition system for measuring the characteristics of real time forces by cage wheels[J].Journal of Terramechanics, 2001,38(4):201-210.

    [5]Yongzhi Zhang, Chuncheng Zuo, Yang Wang, et al. Study of a bionic paddy impeller inspired by buffalo hoof[J].Journal of Bionic Engineering, 2008,5(S1):79-86.

    [6]C Watyotha, V.M Salokhe. Pull, lift and side force characteristics of cage wheels with opposing circumferential lugs[J].Soil and Tillage Research, 2001,60(3-4):123-134.

    [7]陸華忠,羅錫文.水田驅(qū)動(dòng)葉輪輪葉下土壤流動(dòng)特性與動(dòng)力性能研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2010,41(7):50-53.

    [8]張永智,左春檉,孫少明,等.水田驅(qū)動(dòng)葉輪仿生葉片機(jī)理數(shù)值模擬分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2008,39(11):176-179,202.

    [9]邵耀堅(jiān),周大軍,羅錫文.沿前刺孔行駛的水田葉輪動(dòng)力學(xué)的試驗(yàn)研究[J].農(nóng)業(yè)工程學(xué)報(bào),1995,11(1):58-64.

    [10]邵耀堅(jiān),羅錫文.水田耕整機(jī)驅(qū)動(dòng)輪的設(shè)計(jì)和試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),1992,8(2):80-87.

    [11]鐘江,蔣建東, 姜濤,等.基于光滑粒子流體動(dòng)力學(xué)仿真的板結(jié)土壤深旋耕技術(shù)[J].機(jī)械工程學(xué)報(bào),2010,46(19):63-69.

    [12]楊望, 楊堅(jiān), 賈豐云,等.種植于紅粘土的木薯塊根挖掘作業(yè)數(shù)值模擬[J].機(jī)械工程學(xué)報(bào),2013,49(9):135-143.

    [13]LEWIS B A.Manual for LS-DYNA soil material model147[R].Department of Transportation: Federal Highway Administration,2004:10-13.

    [14]胡國(guó)良,任繼文,龍銘.ANSYS13.0有限元分析實(shí)用基礎(chǔ)教程[M].北京:國(guó)防工業(yè)出版社,2012:108-110.

    [15]中國(guó)農(nóng)業(yè)機(jī)械化科學(xué)研究院.農(nóng)業(yè)機(jī)械設(shè)計(jì)手冊(cè)[K].北京:機(jī)械工業(yè)出版社,1988:171-192.

    Abstract ID:1003-188X(2016)04-0029-EA

    Dynamic Simulation of Working Process of Wide Type Paddy Wheel

    Yang Hongyan, Yang Wang, Yang Jian, Lu Hong

    (College of Mechanical Engineering, Guangxi University, Nanning 530004, China)

    Abstract:Wide type paddy-field wheel is a driving device for walking of the micro-tillage machine which is tilling at first and then driving. Except driving performance, the wheel can roll broken soil and stir mixed soil. Its design parameters has great influence on operation quality. “1WGQ4” micro-tillage machine was used as the study object. Dynamic simulation model of soil -wide wheel system was established by using the finite element method and the smooth particle hydrodynamics method. Working process of the wheel was analyzed. And mechanism of interaction between the wheel and soil was studied at meso level. The results show that with rotation of the wheel, contact stress of the soil and lug of the wheel and horizontal propelling force of the lug increase first and then decrease. But the stress reached the maximum first. With the variation of slip rate, the location where the force is zero is different. Bending angle of the wheel has great influence on driving performance, rolling broken soil and stirring mixed soil.

    Key words:wide type paddy wheel; working process; dynamic simulation; mechanism

    文章編號(hào):1003-188X(2016)04-0029-05

    中圖分類號(hào):S222.3

    文獻(xiàn)標(biāo)識(shí)碼:A

    作者簡(jiǎn)介:楊紅艷(1989-),女,安徽亳州人,碩士研究生,(E-mail) zuihongyan1020@163.com。通訊作者:楊望(1984-),男,廣西合浦人,副教授,碩士生導(dǎo)師,博士,(E-mail)yanghope@163.com。

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(51365005) ; 玉林市科技項(xiàng)目(玉市??飘a(chǎn)201306802)

    收稿日期:2015-03-21

    猜你喜歡
    機(jī)理
    揭示S型電荷轉(zhuǎn)移機(jī)理
    隔熱纖維材料的隔熱機(jī)理及其應(yīng)用
    周禮與儒學(xué)的機(jī)理
    煤層氣吸附-解吸機(jī)理再認(rèn)識(shí)
    停用鍋爐的腐蝕機(jī)理及保養(yǎng)措施
    霧霾機(jī)理之問(wèn)
    運(yùn)動(dòng)性血紅蛋白降低的機(jī)理及營(yíng)養(yǎng)防治
    60Co γ-輻照對(duì)硅橡膠GD414損傷機(jī)理的研究
    真空與低溫(2015年4期)2015-06-18 10:47:22
    傳統(tǒng)肥皂的洗滌機(jī)理
    球形ADN的吸濕機(jī)理
    女性被躁到高潮视频| 五月伊人婷婷丁香| 精品一区二区三区视频在线| 2021少妇久久久久久久久久久| 大片电影免费在线观看免费| 国产精品三级大全| 日韩三级伦理在线观看| 五月伊人婷婷丁香| 七月丁香在线播放| 国产视频首页在线观看| 国产一区亚洲一区在线观看| 国产av精品麻豆| 国产欧美日韩综合在线一区二区 | 另类精品久久| 久久午夜综合久久蜜桃| 亚洲精品aⅴ在线观看| 亚洲成色77777| 国产精品麻豆人妻色哟哟久久| 卡戴珊不雅视频在线播放| 一区二区三区四区激情视频| 青春草亚洲视频在线观看| 精品久久国产蜜桃| 日本wwww免费看| xxx大片免费视频| 久久青草综合色| 国产综合精华液| 欧美日韩av久久| 最近的中文字幕免费完整| 一级毛片我不卡| 黄色毛片三级朝国网站 | 国产毛片在线视频| 成人漫画全彩无遮挡| 色94色欧美一区二区| 黄色配什么色好看| 一级毛片黄色毛片免费观看视频| 插逼视频在线观看| 精品酒店卫生间| 亚洲精品,欧美精品| 丰满乱子伦码专区| 天堂俺去俺来也www色官网| 久久久精品免费免费高清| 我的老师免费观看完整版| 美女国产视频在线观看| 日韩亚洲欧美综合| 久热这里只有精品99| 欧美区成人在线视频| 国产成人91sexporn| 99久久综合免费| 内地一区二区视频在线| 黑丝袜美女国产一区| 成人国产麻豆网| 久久97久久精品| 极品教师在线视频| 亚洲成人av在线免费| 搡女人真爽免费视频火全软件| 午夜福利视频精品| 伊人久久国产一区二区| 欧美xxⅹ黑人| 一级a做视频免费观看| 男人爽女人下面视频在线观看| 国产成人精品一,二区| 国产av码专区亚洲av| 少妇人妻 视频| 91精品国产国语对白视频| 国产91av在线免费观看| 亚洲成人av在线免费| 亚洲av男天堂| 精品一区二区三区视频在线| 日韩成人伦理影院| 3wmmmm亚洲av在线观看| 成人无遮挡网站| 十八禁高潮呻吟视频 | 国产精品久久久久久精品电影小说| 免费少妇av软件| 亚洲第一av免费看| www.av在线官网国产| 久久精品国产自在天天线| 国产高清有码在线观看视频| 丝袜在线中文字幕| 黑丝袜美女国产一区| 女人精品久久久久毛片| av女优亚洲男人天堂| 精品视频人人做人人爽| 欧美亚洲 丝袜 人妻 在线| 丰满乱子伦码专区| 国产无遮挡羞羞视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 狂野欧美白嫩少妇大欣赏| 日韩电影二区| 亚洲精品,欧美精品| 人人澡人人妻人| 欧美 日韩 精品 国产| 成人免费观看视频高清| 国产欧美另类精品又又久久亚洲欧美| 狂野欧美激情性bbbbbb| 一个人看视频在线观看www免费| 国产精品无大码| 亚洲精品久久午夜乱码| 国产乱人偷精品视频| 精品国产一区二区三区久久久樱花| 国产成人91sexporn| 欧美 亚洲 国产 日韩一| 国产色婷婷99| 只有这里有精品99| 嫩草影院入口| 午夜影院在线不卡| 国产精品偷伦视频观看了| 少妇裸体淫交视频免费看高清| 丝袜在线中文字幕| 大陆偷拍与自拍| 成人国产av品久久久| 久久精品夜色国产| 看免费成人av毛片| 五月玫瑰六月丁香| 久久99热6这里只有精品| 国产精品99久久久久久久久| www.色视频.com| 国产在视频线精品| 久久亚洲国产成人精品v| av播播在线观看一区| 久久久久久久精品精品| 国产亚洲av片在线观看秒播厂| 日本黄大片高清| 午夜久久久在线观看| 国产中年淑女户外野战色| 麻豆成人av视频| 最黄视频免费看| 久久午夜福利片| 在线观看一区二区三区激情| 国产亚洲5aaaaa淫片| av国产精品久久久久影院| 久久精品国产鲁丝片午夜精品| 一级毛片黄色毛片免费观看视频| 18禁在线播放成人免费| 久久久久久久久久久丰满| 高清午夜精品一区二区三区| 青春草亚洲视频在线观看| 欧美日韩亚洲高清精品| 午夜福利在线观看免费完整高清在| 一边亲一边摸免费视频| 丝袜在线中文字幕| 午夜福利网站1000一区二区三区| 国产黄片视频在线免费观看| 久久久久久久亚洲中文字幕| 欧美日韩在线观看h| 精品少妇内射三级| 日韩视频在线欧美| 一区在线观看完整版| 亚洲国产最新在线播放| 免费看不卡的av| 国产一区二区在线观看日韩| 国产精品99久久99久久久不卡 | 美女大奶头黄色视频| 久久午夜综合久久蜜桃| 极品少妇高潮喷水抽搐| 国产成人精品无人区| 久久久久人妻精品一区果冻| 亚洲精品一区蜜桃| 欧美三级亚洲精品| 色视频在线一区二区三区| 黑人巨大精品欧美一区二区蜜桃 | 欧美激情极品国产一区二区三区 | 观看av在线不卡| videossex国产| 在线观看免费视频网站a站| 免费高清在线观看视频在线观看| 欧美日韩国产mv在线观看视频| 国产淫语在线视频| 人人妻人人澡人人看| 免费观看的影片在线观看| 欧美高清成人免费视频www| 欧美日韩精品成人综合77777| 日韩免费高清中文字幕av| 秋霞伦理黄片| 99久久精品热视频| 黄色视频在线播放观看不卡| 国产欧美亚洲国产| 欧美少妇被猛烈插入视频| 亚洲国产欧美日韩在线播放 | 国产黄色视频一区二区在线观看| 人人妻人人澡人人爽人人夜夜| 高清av免费在线| 插阴视频在线观看视频| 中文资源天堂在线| 丰满人妻一区二区三区视频av| 水蜜桃什么品种好| 国产黄片美女视频| 国产精品一区二区在线不卡| tube8黄色片| 亚洲av日韩在线播放| 丰满迷人的少妇在线观看| 免费黄网站久久成人精品| 在线看a的网站| 9色porny在线观看| 一个人免费看片子| 国产欧美日韩综合在线一区二区 | 久久国产精品大桥未久av | 黄色一级大片看看| 亚洲精品国产成人久久av| 亚洲精品国产av蜜桃| 亚洲欧美精品自产自拍| 啦啦啦视频在线资源免费观看| 国产一区有黄有色的免费视频| 久久99精品国语久久久| 美女内射精品一级片tv| a级一级毛片免费在线观看| 国产 一区精品| 亚洲四区av| 午夜福利,免费看| 亚洲综合色惰| 精华霜和精华液先用哪个| 一本久久精品| 一个人看视频在线观看www免费| 爱豆传媒免费全集在线观看| 亚洲成人一二三区av| 一级毛片黄色毛片免费观看视频| 日日爽夜夜爽网站| 91精品一卡2卡3卡4卡| 国产一区二区在线观看av| 欧美精品亚洲一区二区| 免费观看av网站的网址| 热re99久久国产66热| 99热这里只有是精品50| 亚洲国产av新网站| 欧美xxxx性猛交bbbb| av有码第一页| 丰满迷人的少妇在线观看| 国产色爽女视频免费观看| 久久99精品国语久久久| 丝袜脚勾引网站| 岛国毛片在线播放| 一级毛片久久久久久久久女| 国产色婷婷99| 一级毛片 在线播放| 成人二区视频| 交换朋友夫妻互换小说| 久久午夜福利片| 91精品一卡2卡3卡4卡| 我的女老师完整版在线观看| 国产精品蜜桃在线观看| 欧美亚洲 丝袜 人妻 在线| 最近手机中文字幕大全| 日韩熟女老妇一区二区性免费视频| √禁漫天堂资源中文www| 乱人伦中国视频| 最新的欧美精品一区二区| 夜夜骑夜夜射夜夜干| 人妻一区二区av| 日韩中文字幕视频在线看片| 丝瓜视频免费看黄片| 亚洲av中文av极速乱| 亚洲综合精品二区| 亚洲中文av在线| 亚洲国产欧美在线一区| 水蜜桃什么品种好| h日本视频在线播放| 亚洲av男天堂| 国产成人aa在线观看| 亚洲国产最新在线播放| 女人精品久久久久毛片| √禁漫天堂资源中文www| 简卡轻食公司| 亚洲第一av免费看| 男人舔奶头视频| 男女国产视频网站| 99热6这里只有精品| av在线老鸭窝| 日韩欧美 国产精品| 高清av免费在线| 亚洲成人一二三区av| 久久精品国产亚洲av涩爱| 精品久久久久久久久亚洲| 一级毛片电影观看| 狂野欧美激情性bbbbbb| 80岁老熟妇乱子伦牲交| 国内少妇人妻偷人精品xxx网站| 色婷婷av一区二区三区视频| 亚洲婷婷狠狠爱综合网| av在线app专区| 精品国产一区二区三区久久久樱花| 亚洲熟女精品中文字幕| 黄色怎么调成土黄色| 欧美日韩亚洲高清精品| 波野结衣二区三区在线| 人人妻人人澡人人看| 91精品伊人久久大香线蕉| 亚洲图色成人| 日本黄大片高清| 最近最新中文字幕免费大全7| 国产无遮挡羞羞视频在线观看| av网站免费在线观看视频| videossex国产| 国产成人freesex在线| 免费av中文字幕在线| 一个人免费看片子| 丰满人妻一区二区三区视频av| 午夜精品国产一区二区电影| 3wmmmm亚洲av在线观看| 欧美 日韩 精品 国产| 亚洲精品久久久久久婷婷小说| 亚洲四区av| av天堂中文字幕网| 精品少妇内射三级| 国产黄片美女视频| 毛片一级片免费看久久久久| 精品人妻熟女毛片av久久网站| av一本久久久久| 男人添女人高潮全过程视频| 亚洲自偷自拍三级| 亚洲第一av免费看| 亚洲国产欧美在线一区| 国产精品一二三区在线看| 欧美日韩综合久久久久久| 在线观看三级黄色| 亚洲精品日本国产第一区| 免费观看性生交大片5| 免费观看的影片在线观看| 人妻制服诱惑在线中文字幕| 一边亲一边摸免费视频| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 国产国拍精品亚洲av在线观看| 国产精品一二三区在线看| 久久精品国产亚洲av涩爱| 亚州av有码| 天堂俺去俺来也www色官网| 欧美激情极品国产一区二区三区 | 亚洲性久久影院| 极品人妻少妇av视频| 丝袜脚勾引网站| 国产男女超爽视频在线观看| 男女国产视频网站| 美女cb高潮喷水在线观看| 少妇被粗大的猛进出69影院 | 中文字幕制服av| 天天操日日干夜夜撸| 91精品伊人久久大香线蕉| 久久毛片免费看一区二区三区| 搡女人真爽免费视频火全软件| 亚洲精品成人av观看孕妇| 国产男女超爽视频在线观看| 嫩草影院新地址| 日韩一区二区视频免费看| 久久久久精品久久久久真实原创| 午夜激情福利司机影院| 亚洲av二区三区四区| 国产日韩一区二区三区精品不卡 | 丰满饥渴人妻一区二区三| 日韩,欧美,国产一区二区三区| 久久精品国产自在天天线| 成年人午夜在线观看视频| 成年人免费黄色播放视频 | 在线免费观看不下载黄p国产| 国产av码专区亚洲av| 久久青草综合色| 精品视频人人做人人爽| 国产成人91sexporn| 精品酒店卫生间| 日本vs欧美在线观看视频 | 纵有疾风起免费观看全集完整版| 黄色怎么调成土黄色| 亚洲人成网站在线播| 免费人妻精品一区二区三区视频| 国产精品熟女久久久久浪| 在线精品无人区一区二区三| 日韩中文字幕视频在线看片| 纯流量卡能插随身wifi吗| 免费观看在线日韩| 美女xxoo啪啪120秒动态图| 欧美丝袜亚洲另类| 永久网站在线| 最近中文字幕高清免费大全6| 成人国产麻豆网| 日韩精品有码人妻一区| 国产一级毛片在线| 欧美 日韩 精品 国产| 亚洲情色 制服丝袜| 国产一区二区在线观看av| 十分钟在线观看高清视频www | 九九在线视频观看精品| freevideosex欧美| 亚洲精品乱码久久久久久按摩| 人人妻人人澡人人看| 一级二级三级毛片免费看| 国产精品久久久久久精品电影小说| 妹子高潮喷水视频| h视频一区二区三区| 国产黄片美女视频| 久久久亚洲精品成人影院| 极品人妻少妇av视频| 久久午夜综合久久蜜桃| 亚洲欧美清纯卡通| 国产精品蜜桃在线观看| 欧美高清成人免费视频www| 三级国产精品欧美在线观看| 成人黄色视频免费在线看| 日韩av免费高清视频| 啦啦啦视频在线资源免费观看| 国产日韩欧美视频二区| 免费人成在线观看视频色| 免费观看av网站的网址| 成人国产麻豆网| 久久影院123| 亚洲综合精品二区| 丝袜在线中文字幕| 熟女av电影| a级毛片在线看网站| 色吧在线观看| 久久97久久精品| 欧美xxxx性猛交bbbb| 欧美一级a爱片免费观看看| 三级经典国产精品| 国产黄频视频在线观看| 中国国产av一级| 狂野欧美白嫩少妇大欣赏| 亚洲第一av免费看| 亚洲经典国产精华液单| 视频区图区小说| 男男h啪啪无遮挡| 在线亚洲精品国产二区图片欧美 | 精品99又大又爽又粗少妇毛片| 成年女人在线观看亚洲视频| 成人毛片60女人毛片免费| 国产成人免费观看mmmm| 久久国产亚洲av麻豆专区| 欧美丝袜亚洲另类| kizo精华| 欧美日韩在线观看h| 伦理电影免费视频| 日韩免费高清中文字幕av| 亚洲不卡免费看| 国产在线免费精品| 国产精品福利在线免费观看| 在线观看美女被高潮喷水网站| 超碰97精品在线观看| 欧美日本中文国产一区发布| 亚洲欧美成人综合另类久久久| 国产成人a∨麻豆精品| 精品亚洲成国产av| 三级国产精品欧美在线观看| 免费久久久久久久精品成人欧美视频 | 久久鲁丝午夜福利片| 日本午夜av视频| 在线观看免费高清a一片| 国产亚洲av片在线观看秒播厂| 91久久精品电影网| 午夜久久久在线观看| 国产日韩欧美视频二区| 亚洲欧美日韩东京热| 99久久精品热视频| 亚洲综合精品二区| 久热这里只有精品99| 国产精品国产三级专区第一集| 久久精品久久精品一区二区三区| 内地一区二区视频在线| 欧美最新免费一区二区三区| 777米奇影视久久| 欧美 亚洲 国产 日韩一| 久久婷婷青草| 91aial.com中文字幕在线观看| 国产成人a∨麻豆精品| 少妇精品久久久久久久| 国产亚洲5aaaaa淫片| 王馨瑶露胸无遮挡在线观看| 日韩欧美一区视频在线观看 | 日韩免费高清中文字幕av| 亚洲人成网站在线观看播放| 久久久久久久久久久久大奶| 亚洲精品国产成人久久av| 久久精品国产鲁丝片午夜精品| 亚洲精品国产成人久久av| 又黄又爽又刺激的免费视频.| 如日韩欧美国产精品一区二区三区 | 日本黄色日本黄色录像| 噜噜噜噜噜久久久久久91| 国产日韩欧美在线精品| 婷婷色麻豆天堂久久| 亚洲成人手机| 国产欧美日韩精品一区二区| 最新中文字幕久久久久| 亚洲va在线va天堂va国产| 黄色一级大片看看| 91在线精品国自产拍蜜月| 亚洲国产精品999| xxx大片免费视频| 精品亚洲成a人片在线观看| 五月伊人婷婷丁香| 99热这里只有是精品在线观看| 永久网站在线| 日韩制服骚丝袜av| 丰满人妻一区二区三区视频av| 婷婷色综合大香蕉| 永久网站在线| 国产在线视频一区二区| 久热久热在线精品观看| 久久久久精品性色| 嫩草影院新地址| 亚洲av.av天堂| 汤姆久久久久久久影院中文字幕| 人人妻人人爽人人添夜夜欢视频 | 色吧在线观看| 99热6这里只有精品| av有码第一页| 成人影院久久| 97在线人人人人妻| 国产午夜精品久久久久久一区二区三区| 99久久综合免费| 久久亚洲国产成人精品v| 极品人妻少妇av视频| 另类精品久久| 国产免费一级a男人的天堂| 啦啦啦视频在线资源免费观看| h视频一区二区三区| 99久久精品热视频| 国产精品国产三级国产专区5o| 人人妻人人澡人人爽人人夜夜| 成人影院久久| 日本91视频免费播放| av播播在线观看一区| 亚洲精品第二区| 交换朋友夫妻互换小说| 亚洲精品国产av成人精品| 在线精品无人区一区二区三| 久久久国产一区二区| 久久精品国产鲁丝片午夜精品| 亚洲av中文av极速乱| 在现免费观看毛片| 成人亚洲精品一区在线观看| 91在线精品国自产拍蜜月| 在线播放无遮挡| 自拍欧美九色日韩亚洲蝌蚪91 | 交换朋友夫妻互换小说| av卡一久久| 亚洲,欧美,日韩| 97在线视频观看| 老女人水多毛片| 亚洲精品日韩av片在线观看| 99热6这里只有精品| av网站免费在线观看视频| 人妻人人澡人人爽人人| 女人久久www免费人成看片| 日韩成人av中文字幕在线观看| 国产日韩欧美在线精品| 看非洲黑人一级黄片| 久久精品国产亚洲网站| 99视频精品全部免费 在线| 我的女老师完整版在线观看| 亚洲精品第二区| 欧美日韩一区二区视频在线观看视频在线| 亚洲综合色惰| 99热全是精品| 国产免费一级a男人的天堂| 一区二区三区乱码不卡18| 日韩,欧美,国产一区二区三区| 伦理电影免费视频| 日韩一区二区视频免费看| 七月丁香在线播放| 亚洲精品中文字幕在线视频 | 一本—道久久a久久精品蜜桃钙片| 香蕉精品网在线| 亚洲精品乱码久久久久久按摩| 麻豆成人午夜福利视频| 国产成人精品久久久久久| 一本大道久久a久久精品| 肉色欧美久久久久久久蜜桃| 日本欧美国产在线视频| 国产免费福利视频在线观看| 亚洲精品日本国产第一区| 永久网站在线| 亚洲国产色片| 亚洲一区二区三区欧美精品| 少妇人妻精品综合一区二区| 亚洲精品视频女| 啦啦啦啦在线视频资源| 91aial.com中文字幕在线观看| 亚洲精品国产av蜜桃| 亚洲综合色惰| 久久久久久久亚洲中文字幕| 爱豆传媒免费全集在线观看| 大又大粗又爽又黄少妇毛片口| 少妇的逼好多水| 午夜视频国产福利| 9色porny在线观看| 99国产精品免费福利视频| 少妇裸体淫交视频免费看高清| 精品人妻熟女毛片av久久网站| 国产av国产精品国产| 看十八女毛片水多多多| 丰满少妇做爰视频| 欧美日韩综合久久久久久| 国产精品熟女久久久久浪| 国产精品一区二区性色av| 51国产日韩欧美| 欧美日韩在线观看h| 乱人伦中国视频| 中国三级夫妇交换| 国产色婷婷99| 青春草亚洲视频在线观看| 亚洲在久久综合| 亚洲高清免费不卡视频| 精品熟女少妇av免费看| 99久国产av精品国产电影| 亚洲精品亚洲一区二区| 天堂俺去俺来也www色官网| 深夜a级毛片| 精品亚洲乱码少妇综合久久| av免费在线看不卡| 一级毛片aaaaaa免费看小| 国产精品女同一区二区软件| 久久久久久久久久成人| 少妇猛男粗大的猛烈进出视频| 你懂的网址亚洲精品在线观看| 男人和女人高潮做爰伦理| 一级,二级,三级黄色视频| 啦啦啦啦在线视频资源|