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

    基于離散元法的川產(chǎn)法半夏仿真參數(shù)標(biāo)定

    2023-07-12 03:16:22謝宇鴻王金棟吳展揚(yáng)
    現(xiàn)代農(nóng)業(yè)裝備 2023年3期
    關(guān)鍵詞:恢復(fù)系數(shù)摩擦系數(shù)螢火蟲

    謝宇鴻,王金棟,陳 燚,吳展揚(yáng)

    (西南交通大學(xué)機(jī)械工程學(xué)院,四川 成都 610031)

    0 引言

    半夏是一種典型川產(chǎn)中藥材,在鎮(zhèn)咳祛痰、抗腫瘤方面具有明顯療效,法半夏則是半夏通過(guò)生石灰、甘草汁等進(jìn)行炮制后的加工品。半夏的粉碎細(xì)度會(huì)影響其有效成分的析出[1],采用離散元法進(jìn)行半夏粉碎的仿真分析,有利于揭示半夏的粉碎機(jī)理,對(duì)粉碎設(shè)備研磨副的研制、提高半夏粉碎細(xì)度具有指導(dǎo)意義。

    離散元法(Discrete element method,DEM)是將研究對(duì)象視為由大量顆粒組成的離散系統(tǒng),通過(guò)研究顆粒間的接觸關(guān)系、受力情況和運(yùn)動(dòng)狀態(tài)等進(jìn)而獲取研究對(duì)象的宏觀和微觀物理量信息的研究方法。經(jīng)過(guò)數(shù)十年的發(fā)展,離散元法已經(jīng)廣泛用于工農(nóng)業(yè)領(lǐng)域[2-4]。半夏的粉碎加工是一個(gè)典型粉體加工過(guò)程,半夏與半夏、半夏與研磨副間具有復(fù)雜的相對(duì)運(yùn)動(dòng)和碰撞,因此適于采用離散元法進(jìn)行仿真分析。

    為確保離散元法分析的準(zhǔn)確性,需確立研究對(duì)象準(zhǔn)確的物性參數(shù)[5-8]。物性參數(shù)主要分為物體的本征參數(shù)和物體間的接觸參數(shù),本征參數(shù)(包括泊松比、剪切模量和密度等)為物體自身固有的特性參數(shù);接觸參數(shù)(包括碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動(dòng)摩擦系數(shù))為表征兩物體接觸特性的參數(shù),由于復(fù)雜的接觸特性通常無(wú)法通過(guò)直接測(cè)量獲取,因此需要對(duì)接觸參數(shù)進(jìn)行標(biāo)定。目前,眾多國(guó)內(nèi)外學(xué)者開(kāi)展了有關(guān)離散元仿真參數(shù)標(biāo)定的研究。張勇等[9]通過(guò)臺(tái)架試驗(yàn)和理論計(jì)算獲取綠豆與不同材料(鋼板、有機(jī)塑料板和橡膠)間的滑動(dòng)摩擦系數(shù)和碰撞恢復(fù)系數(shù)。于慶旭等[10]基于離散元仿真和臺(tái)架試驗(yàn),對(duì)三七種子間、三七種子與ABS 塑料板間的接觸參數(shù)進(jìn)行了標(biāo)定。劉敏等[11]通過(guò)實(shí)際試驗(yàn)獲取紫花苜宿種子的本征參數(shù),以休止角為指標(biāo)標(biāo)定了紫花苜宿種子與鋼板的接觸參數(shù)。田劍鋒等[12]通過(guò)修改滾動(dòng)摩擦系數(shù),使仿真試驗(yàn)所得的堆積角與實(shí)際堆積角接近的方法獲取黨參種子間的滾動(dòng)摩擦系數(shù)。戴念祖等[13]、吳孟宸等[14]、馬文鵬等[15]以堆積角相對(duì)誤差為指標(biāo)進(jìn)行無(wú)物料的接觸參數(shù)標(biāo)定。鹿芳媛等[16]、RORATO 等[17],XIA 等[18]運(yùn)用圖 像處理 技術(shù),以堆積角為指標(biāo)分別對(duì)水稻、砂石和濕煤粉進(jìn)行參數(shù)標(biāo)定。邢潔潔等[19]、BAHRAMI 等[20],采用離散元法對(duì)特定土壤進(jìn)行仿真參數(shù)標(biāo)定。GHODKI等[21]采用側(cè)壁坍塌試驗(yàn)對(duì)大豆的離散元仿真參數(shù)進(jìn)行標(biāo)定。YE 等[22]運(yùn)用BP 神經(jīng)網(wǎng)絡(luò)進(jìn)行顆粒的參數(shù)標(biāo)定。ROESSLER 等[23]考慮到最優(yōu)參數(shù)的不確定性,通過(guò)對(duì)比堆積試驗(yàn)、側(cè)壁坍塌試驗(yàn)和抬升試驗(yàn),發(fā)現(xiàn)不同試驗(yàn)方法和指標(biāo)結(jié)合能夠有效減小最優(yōu)參數(shù)組合的數(shù)量。綜合國(guó)內(nèi)外研究,目前對(duì)半夏的接觸參數(shù)標(biāo)定鮮有涉及,導(dǎo)致無(wú)法建立準(zhǔn)確的半夏離散元仿真分析模型,從而對(duì)半夏高效粉碎加工及有效成分提取帶來(lái)較大影響。

    本文以干燥的川產(chǎn)法半夏(以下簡(jiǎn)稱半夏)為研究對(duì)象,進(jìn)行半夏的接觸參數(shù)標(biāo)定。采用EDEM軟件建立半夏的離散元模型,并結(jié)合物理試驗(yàn)和仿真試驗(yàn)標(biāo)定半夏與粉碎設(shè)備研磨副材料(4Cr13 不銹鋼)間的接觸參數(shù);利用堆積試驗(yàn)和圖像識(shí)別方法獲取半夏堆積角及誤差,采用響應(yīng)面法(RSM)構(gòu)建誤差與接觸參數(shù)的數(shù)學(xué)模型;以誤差最小為目標(biāo),通過(guò)螢火蟲算法(FA)求解半夏間最優(yōu)接觸參數(shù)組合;采用側(cè)壁坍塌試驗(yàn)進(jìn)行參數(shù)組合準(zhǔn)確性的試驗(yàn)驗(yàn)證。

    1 材料與離散元模型

    1.1 材料的本征參數(shù)

    半夏的整體形狀為稍扁斜的半球形,下端鈍圓較光滑,頂端有凹狀莖痕,莖痕周邊密布點(diǎn)狀根痕,實(shí)物如圖1 所示。隨機(jī)選取100 粒半夏,采用分度值為0.02 mm 的游標(biāo)卡尺測(cè)量半夏的外形尺寸,統(tǒng)計(jì)得半夏的直徑為11.70~15.16 mm,高度為8.14~11.64 mm。采用檢定分度為0.1 g 的電子天平進(jìn)行稱重,將半夏置于WH9020A 型電熱恒溫干燥箱(杭州聚同電子有限公司),105 ℃烘干至恒重,計(jì)算法制旱半夏含水率6.25%~9.63%。采用精度為1 mL的量筒,通過(guò)液浸法測(cè)量半夏的密度,統(tǒng)計(jì)得半夏密度為963~1 153 kg/m3。

    圖1 半夏實(shí)物

    采用壓縮試驗(yàn)測(cè)試半夏的彈性模量、泊松比和剪切模量。將半夏切割為高度5 mm,直徑8 mm 的圓柱體試樣,采用TMS-PRO 研究型食品物性分析儀(美國(guó)FTC 公司),使用直徑為10 mm 的圓形探頭,沿試樣的軸線方向,按設(shè)置的加載速度和加載時(shí)間施加載荷,試驗(yàn)重復(fù)5 次,并記錄半夏在壓縮過(guò)程中的載荷-位移數(shù)據(jù)。根據(jù)公式(1)計(jì)算半夏的彈性模量

    式中:

    E——彈性模量,Pa;

    σ ——最大壓應(yīng)力,Pa;

    ε1——半夏的縱向線應(yīng)變。

    根據(jù)公式(2)計(jì)算半夏的泊松比

    式中:

    μ——泊松比;

    ε2——半夏的橫線線應(yīng)變。

    根據(jù)公式(3)計(jì)算半夏的剪切模量

    式中:

    G——剪切模量。

    4Cr13 不銹鋼的相關(guān)材料參數(shù)通過(guò)查表可得。結(jié)合壓縮試驗(yàn)、EDEM 材料數(shù)據(jù)庫(kù)及相關(guān)物料參數(shù)[24-27],本研究設(shè)定半夏及4Cr13 不銹鋼參數(shù)如表1 所示。

    表1 試驗(yàn)材料本征參數(shù)

    1.2 半夏的離散元模型

    由于實(shí)際物體的形狀并不規(guī)則,在離散元分析軟件EDEM 中,通常采用多球聚合模型(Multisphere method,MSM)和粘結(jié)顆粒模型(Bonded particle method,BPM)構(gòu)建離散元模型,離散元模型還原度越高,分析誤差越小。MSM 通過(guò)若干不同大小的球形顆粒堆疊,建立不同形狀的模型;BPM通過(guò)若干相同大小、互不重疊的球形顆粒,利用粘結(jié)鍵連接形成物體的離散元模型,互相粘結(jié)的球形顆粒半徑越小,所需顆粒數(shù)量越多,離散元模型越逼近物體輪廓,所需的仿真時(shí)間越長(zhǎng)。

    由于半夏的外形較為簡(jiǎn)單、規(guī)則,本研究將其簡(jiǎn)化為下端鈍圓的半球形,使用三維建模軟件NX 12.0,以外形尺寸平均值建立半夏的三維模型,并導(dǎo)出為.msh 文件格式,最終獲得半夏的輪廓網(wǎng)格模型(圖2a)。將輪廓網(wǎng)格模型導(dǎo)入EDEM 軟件,分別采用多球聚合模型和粘結(jié)顆粒模型構(gòu)建半夏的離散元模型。模型如圖2 所示,半夏的粘結(jié)顆粒模型(圖2b)由82 300 顆半徑均為0.15 mm 的球形顆粒組成,半夏的多球聚合模型(圖2c)由11 顆不同半徑的球形顆粒組成,綜合考慮仿真試驗(yàn)的精度和效率,選擇多球聚合模型所構(gòu)建的離散元模型。

    2 接觸參數(shù)標(biāo)定試驗(yàn)與結(jié)果討論

    半夏所需標(biāo)定的接觸參數(shù)有:半夏與不銹鋼的碰撞恢復(fù)系數(shù)a、半夏與不銹鋼的靜摩擦系數(shù)b、半夏與不銹鋼的滾動(dòng)摩擦系數(shù)c;半夏間的碰撞恢復(fù)系數(shù)A,半夏間的靜摩擦系數(shù)B和半夏間的滾動(dòng)摩擦系數(shù)C。由于實(shí)際半夏的外形與簡(jiǎn)化的離散元模型存在差異,為減小標(biāo)定的接觸參數(shù)與實(shí)際接觸參數(shù)的誤差,本文將采用實(shí)際試驗(yàn)與離散元仿真試驗(yàn)結(jié)合的方式對(duì)半夏與不銹鋼的接觸參數(shù)進(jìn)行標(biāo)定。

    2.1 半夏-不銹鋼接觸參數(shù)標(biāo)定

    2.1.1 半夏-不銹鋼的碰撞恢復(fù)系數(shù)標(biāo)定

    如公式(4)所示,碰撞恢復(fù)系數(shù)定義為:碰撞后兩物體相對(duì)分離速度與碰撞前相對(duì)接近速度之比,只與物體的材料有關(guān),用于表征物體碰撞后恢復(fù)原形的能力。

    假設(shè)物體1 在運(yùn)動(dòng)過(guò)程中只受重力影響,且物體2 碰撞前后速度均為0,則公式(4)可簡(jiǎn)化為

    式中:

    g——重力加速度,m/s2;

    H——物體初始下落高度,mm;

    h——物體最大彈起高度,mm。

    根據(jù)公式(5),本文設(shè)計(jì)了圖3 所示彈跳試驗(yàn)方案,規(guī)定半夏中心到不銹鋼板的距離為半夏高度,將半夏從H=105 mm 處自由下落,與不銹鋼板碰撞后彈起,采用高速相機(jī)拍攝全過(guò)程,并記錄彈起的最大高度。試驗(yàn)重復(fù)5 次取平均值,確定半夏平均最大彈起高度為33.76 mm,碰撞恢復(fù)系數(shù)為0.567。將碰撞恢復(fù)系數(shù)a輸入EDEM 軟件,接觸參數(shù)b、c、A、B、C均設(shè)置為0,進(jìn)行5 次半夏下落碰撞仿真試驗(yàn),得到平均最大彈起高度為35.12 mm,與實(shí)際高度的相對(duì)誤差為4.03%,小于5%,表明標(biāo)定的碰撞恢復(fù)系數(shù)準(zhǔn)確性高。確定半夏-不銹鋼的碰撞恢復(fù)系數(shù)為0.567。

    2.1.2 半夏-不銹鋼的靜摩擦系數(shù)標(biāo)定

    采用斜面試驗(yàn)進(jìn)行靜摩擦系數(shù)的標(biāo)定。試驗(yàn)原理如圖4 所示,將半夏置于不銹鋼板上表面,不銹鋼板繞一棱線緩慢傾斜,當(dāng)半夏開(kāi)始滑移時(shí)立即停止轉(zhuǎn)動(dòng),記錄此時(shí)平板的傾斜角度y1。仿真試驗(yàn)中,為消除其他參數(shù)的影響,設(shè)置接觸參數(shù)a=0.567,其余接觸參數(shù)c、A、B、C均設(shè)置為0;通過(guò)最陡爬坡試驗(yàn)確定靜摩擦系數(shù)的取值范圍為0.35~0.70,步長(zhǎng)取0.05,進(jìn)行8 組仿真試驗(yàn)。實(shí)際試驗(yàn)采用分度值為0.01° 的數(shù)顯角度尺測(cè)量?jī)A斜角度;仿真試驗(yàn)中采用EDEM 量角器工具測(cè)量,每組試驗(yàn)重復(fù)5 次取平均值。得到不銹鋼板的實(shí)際平均傾斜角度為34.72°,仿真試驗(yàn)結(jié)果如表2 所示。

    表2 靜摩擦系數(shù)仿真試驗(yàn)方案與結(jié)果

    圖4 斜面試驗(yàn)

    將表中數(shù)據(jù)采用Matlab 軟件進(jìn)行多項(xiàng)式擬合,擬合曲線如圖5 所示,曲線函數(shù)為

    圖5 仿真試驗(yàn)靜摩擦系數(shù)與平板傾斜角度擬合曲線

    函數(shù)決定系數(shù)R2=0.996 5,趨近1,說(shuō)明式(6)擬合精度高。將實(shí)際傾斜角度34.72°代入上式,解得b=0.649。將求解所得靜摩擦系數(shù)輸入EDEM 軟件進(jìn)行5 次重復(fù)試驗(yàn)取平均值,得到傾斜角度36.23°,與實(shí)際角度的相對(duì)誤差為4.35%,小于5%,表明標(biāo)定的靜摩擦系數(shù)準(zhǔn)確性高。通過(guò)上述過(guò)程,確定半夏-不銹鋼的靜摩擦系數(shù)為0.649。

    2.1.3 半夏-不銹鋼的滾動(dòng)摩擦系數(shù)標(biāo)定

    半夏與不銹鋼的滾動(dòng)摩擦系數(shù)試驗(yàn)如圖6 所示,將半夏從斜度為30°的斜坡頂端靜止釋放,并沿斜坡滾動(dòng)至水平面,記錄半夏在水平面上滾動(dòng)的最遠(yuǎn)距離y2。在離散元仿真試驗(yàn)中,為消除其他參數(shù)的影響,設(shè)置接觸參數(shù)a=0.567,b=0.649,其余接觸參數(shù)A、B、C均設(shè)置為0;通過(guò)最陡爬坡試驗(yàn)確定滾動(dòng)摩擦系數(shù)的取值范圍為0.06~0.13,步長(zhǎng)取0.01,進(jìn)行8 組仿真試驗(yàn)。實(shí)際試驗(yàn)采用精度為0.02 mm 的游標(biāo)卡尺進(jìn)行測(cè)量;仿真試驗(yàn)中采用EDEM 長(zhǎng)度測(cè)量工具,每組試驗(yàn)重復(fù)5 次取平均值。得到半夏的實(shí)際平均最遠(yuǎn)滾動(dòng)距離為149.40 mm,仿真試驗(yàn)結(jié)果如表3 所示。

    表3 滾動(dòng)摩擦系數(shù)仿真試驗(yàn)方案及結(jié)果

    圖6 斜坡滾動(dòng)試驗(yàn)

    將表中數(shù)據(jù)采用Matlab 軟件進(jìn)行二次多項(xiàng)式擬合,擬合曲線如圖7 所示,曲線函數(shù)為

    圖7 仿真滾動(dòng)摩擦系數(shù)與最大水平滾動(dòng)距離擬合曲線

    函數(shù)決定系數(shù)R2=0.994 1,趨近1,說(shuō)明式(7)擬合精度高。將實(shí)際最遠(yuǎn)距離149.40 mm 代入上式,解得c=0.103。將求解所得靜摩擦系數(shù)輸入EDEM軟件進(jìn)行5 次重復(fù)試驗(yàn)取均值,得到最遠(yuǎn)滾動(dòng)距離143.62 mm,與實(shí)際距離的相對(duì)誤差為3.87%,小于5%,表明標(biāo)定的滾動(dòng)摩擦系數(shù)準(zhǔn)確性高。通過(guò)上述過(guò)程,確定半夏-不銹鋼的滾動(dòng)摩擦系數(shù)為0.103。

    2.2 半夏-半夏的接觸參數(shù)標(biāo)定

    堆積角是大量顆粒自由下落堆積所形成錐體的自由表面與水平面的夾角,主要受顆粒表面形態(tài)、摩擦系數(shù)等因素影響。由于采用試驗(yàn)直接進(jìn)行半夏間的接觸參數(shù)標(biāo)定存在困難,因此本文采用堆積試驗(yàn),以仿真試驗(yàn)堆積角與實(shí)際試驗(yàn)堆積角的相對(duì)誤差為指標(biāo)間接求解半夏間的最佳接觸參數(shù)組合。

    2.2.1 堆積試驗(yàn)

    如圖8(a)所示,試驗(yàn)裝置由帶孔的有機(jī)玻璃箱和不銹鋼圓盤組成。試驗(yàn)過(guò)程中,使用擋板將種箱孔封堵,再填入種箱體積2/3 左右的半夏,迅速抽離擋板使半夏從孔中自由下落,半夏在圓盤上堆積形成錐體,并采用相機(jī)拍攝錐體正視圖。如圖8(b)所示,錐體斜面與水平面的內(nèi)角即為堆積角。重復(fù)試驗(yàn)結(jié)果表明錐體兩側(cè)的堆積角相對(duì)誤差小于5%,本文選擇錐體右側(cè)進(jìn)行堆積角測(cè)量。

    圖8 堆積試驗(yàn)

    為降低測(cè)量誤差,采用圖像識(shí)別技術(shù)進(jìn)行堆積角的測(cè)量。如圖9 所示,依次對(duì)圖像進(jìn)行灰度處理、二值化、邊緣檢測(cè)和邊緣直線擬合,并對(duì)堆積角角度進(jìn)行測(cè)量。試驗(yàn)重復(fù)5 次,得到半夏的實(shí)際堆積角平均值為39.65°。

    圖9 圖像處理

    在仿真試驗(yàn)中,將試驗(yàn)裝置的1∶1 三維模型導(dǎo)入EDEM 軟件中,半夏離散元模型選擇多球堆疊模型,顆粒接觸模型選擇Hertz-Mindlin 無(wú)滑移模型,設(shè)置接觸參數(shù)a=0.567,b=0.649,c=0.103,仿真試驗(yàn)流程與實(shí)際試驗(yàn)相同。

    2.2.2 仿真試驗(yàn)方案設(shè)計(jì)及試驗(yàn)結(jié)果

    為構(gòu)建準(zhǔn)確的半夏間接觸參數(shù)與堆積角相對(duì)誤差的多元回歸模型,本文基于三因素二次旋轉(zhuǎn)正交組合試驗(yàn)方案設(shè)計(jì)仿真堆積試驗(yàn)。

    如表4 所示,首先,依據(jù)最陡爬坡試驗(yàn)結(jié)果,分別以第3、4、5 組的碰撞恢復(fù)系數(shù)A、靜摩擦系數(shù)B和滾動(dòng)摩擦系數(shù)C為正交試驗(yàn)的-1、0 和1 水平因素。再采用Design-Expert 12 軟件進(jìn)行仿真試驗(yàn)因素編碼和正交試驗(yàn)方案設(shè)計(jì)。仿真試驗(yàn)因素編碼如表5 所示,仿真試驗(yàn)方案及結(jié)果如表6 所示,Y為仿真試驗(yàn)堆積角與實(shí)際堆積角的相對(duì)誤差。

    表4 最陡爬坡試驗(yàn)方案及結(jié)果

    表5 仿真試驗(yàn)因素編碼

    表6 試驗(yàn)設(shè)計(jì)方案及結(jié)果

    2.2.3 仿真試驗(yàn)結(jié)果分析

    對(duì)仿真試驗(yàn)結(jié)果進(jìn)行多元回歸擬合,得到接觸參數(shù)對(duì)堆積角相對(duì)誤差的回歸方程為

    回歸模型的方差分析結(jié)果如表7 所示,該模型的擬合度極其顯著(P<0.000 1)。其中A、B、C、B2和C2的P值均小于0.01,表明相關(guān)因素對(duì)堆積角相對(duì)誤差有顯著影響。失擬項(xiàng)P值為0.081 6,不顯著,表明無(wú)其他顯著因素影響試驗(yàn)指標(biāo)。方差分析結(jié)果驗(yàn)證了影響堆積角的主要因素為半夏間的碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動(dòng)摩擦系數(shù)。

    表7 回歸模型方差分析

    2.2.4 求解最優(yōu)參數(shù)組合

    以半夏間的碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動(dòng)摩擦系數(shù)為優(yōu)化對(duì)象,以堆積角相對(duì)誤差最小為優(yōu)化目標(biāo)進(jìn)行最優(yōu)參數(shù)組合研究。目標(biāo)函數(shù)和約束函數(shù)為

    螢火蟲算法(Firefly Algorithm,FA)是受自然界螢火蟲發(fā)光行為啟發(fā)的一種新型群智能優(yōu)化算法[28],該算法提出以下前提:①所有螢火蟲不分性別,每只螢火蟲會(huì)受比它更亮的螢火蟲吸引;②吸引力與發(fā)光強(qiáng)度成正比;③個(gè)體的亮度由目標(biāo)函數(shù)值決定。依據(jù)以上前提,定義螢火蟲i在r=0 時(shí)的亮度為絕對(duì)亮度,即為Ii,則螢火蟲i在螢火蟲j處的相對(duì)亮度為

    式中:

    γ——光吸收系數(shù);

    rij——螢火蟲i和j的距離。

    由于發(fā)光強(qiáng)度與吸引力成比例,則螢火蟲i對(duì)螢火蟲j的吸引力表示為

    式中:

    β0——r=0 時(shí)的吸引力。

    此時(shí)螢火蟲j受螢火蟲i的吸引,螢火蟲j的位置發(fā)生變化,位置更新公式為

    式中:

    xi、xj——螢火蟲i和j在空間中的位置;

    k——迭代次數(shù);

    α——步長(zhǎng);

    ε——高斯分布的隨機(jī)向量。

    螢火蟲算法的基本流程如圖10 所示。

    圖10 螢火蟲算法基本流程

    相較于其他優(yōu)化算法,螢火蟲算法具有設(shè)置參數(shù)少,操作簡(jiǎn)便穩(wěn)定的優(yōu)點(diǎn)。由于每只螢火蟲具有獨(dú)立的感知和搜索能力,螢火蟲算法不僅具有局部尋優(yōu)能力,也能尋找全局最優(yōu)解,因此,在求解連續(xù)的峰值函數(shù)時(shí),螢火蟲算法具有一定的優(yōu)勢(shì)?;谖灮鹣x算法的求解優(yōu)勢(shì),本文采用螢火蟲算法求解半夏間最優(yōu)參數(shù)組合,得到半夏間最優(yōu)接觸參數(shù)組合為:碰撞恢復(fù)系數(shù)為0.584,靜摩擦系數(shù)為0.293,滾動(dòng)摩擦系數(shù)為0.084,此時(shí)堆積角的最小相對(duì)誤差為0.58%。

    3 試驗(yàn)驗(yàn)證

    采用半夏的側(cè)壁坍塌試驗(yàn)對(duì)上述試驗(yàn)標(biāo)定的接觸參數(shù)組合準(zhǔn)確性進(jìn)行分析驗(yàn)證。如圖11 所示,在有機(jī)玻璃箱內(nèi)填充2/3 體積的半夏,并迅速抽離右側(cè)擋板,使半夏從容器右側(cè)自然坍塌在不銹鋼板上,待坍塌過(guò)程穩(wěn)定后,采用相機(jī)拍攝半夏堆體的正視圖。

    圖11 半夏側(cè)壁坍塌試驗(yàn)

    半夏堆體的斜面與水平面的夾角定義為休止角,休止角的測(cè)量方法與堆積角的測(cè)量方法相同,測(cè)量結(jié)果表明,實(shí)際休止角角度為27.08°。仿真試驗(yàn)中,將已標(biāo)定的6 個(gè)接觸參數(shù)輸入EDEM 軟件進(jìn)行半夏的仿真坍塌試驗(yàn),試驗(yàn)重復(fù)5 次。試驗(yàn)結(jié)果如表8所示,5 次試驗(yàn)的休止角相對(duì)誤差均小于5%。綜上,本文建立的半夏離散元模型及接觸參數(shù),在誤差允許范圍內(nèi)可為半夏的離散元法仿真試驗(yàn)提供參考。

    表8 側(cè)壁坍塌試驗(yàn)驗(yàn)證結(jié)果

    4 結(jié)論

    1)基于多球聚合模型,建立半夏離散元模型。依據(jù)碰撞恢復(fù)系數(shù)定義,設(shè)計(jì)彈跳試驗(yàn),標(biāo)定半夏-不銹鋼碰撞恢復(fù)系數(shù)為0.567。采用試驗(yàn)與仿真相結(jié)合的方法,通過(guò)斜面試驗(yàn)和斜坡滾動(dòng)試驗(yàn)分別對(duì)半夏-不銹鋼靜摩擦系數(shù)、半夏-不銹鋼滾動(dòng)摩擦系數(shù)進(jìn)行標(biāo)定,建立目標(biāo)參數(shù)與試驗(yàn)指標(biāo)的數(shù)學(xué)模型,求解得半夏-不銹鋼靜摩擦系數(shù)為0.649,半夏-不銹鋼滾動(dòng)摩擦系數(shù)為0.103。

    2)基于堆積試驗(yàn),以實(shí)際試驗(yàn)和仿真試驗(yàn)的堆積角相對(duì)誤差為指標(biāo),對(duì)半夏間的接觸參數(shù)進(jìn)行標(biāo)定。以半夏間碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動(dòng)摩擦系數(shù)為試驗(yàn)因素,設(shè)計(jì)了三因素二次旋轉(zhuǎn)正交組合仿真試驗(yàn)方案,采用圖像識(shí)別技術(shù)測(cè)量半夏堆積角。建立半夏間接觸參數(shù)與堆積角相對(duì)誤差的多元回歸模型,并通過(guò)方差分析和響應(yīng)面分析,驗(yàn)證了回歸模型擬合的顯著性(P<0.000 1)、表明了各因素及其交互作用對(duì)指標(biāo)存在顯著影響。

    3)以堆積角相對(duì)誤差最小為優(yōu)化目標(biāo),采用螢火蟲算法對(duì)半夏間接觸參數(shù)與堆積角相對(duì)誤差的多元回歸模型進(jìn)行求解,得到半夏間碰撞恢復(fù)系數(shù)為0.584,靜摩擦系數(shù)為0.293,滾動(dòng)摩擦系數(shù)為0.084。側(cè)壁坍塌試驗(yàn)驗(yàn)證結(jié)果表明,5 次試驗(yàn)休止角相對(duì)誤差均小于5%,因此該半夏離散元模型和接觸參數(shù)可用于后續(xù)半夏的離散元法粉碎仿真試驗(yàn)。

    猜你喜歡
    恢復(fù)系數(shù)摩擦系數(shù)螢火蟲
    隧道內(nèi)水泥混凝土路面微銑刨后摩擦系數(shù)衰減規(guī)律研究
    中外公路(2022年1期)2022-05-14 08:13:26
    落石法向恢復(fù)系數(shù)的多因素聯(lián)合影響研究
    利用恢復(fù)系數(shù)巧解碰撞問(wèn)題
    摩擦系數(shù)對(duì)直齒輪副振動(dòng)特性的影響
    螢火蟲
    螢火蟲
    抱抱就不哭了
    夏天的螢火蟲
    落石碰撞法向恢復(fù)系數(shù)的模型試驗(yàn)研究
    CSP生產(chǎn)線摩擦系數(shù)與軋制力模型的研究
    上海金屬(2014年3期)2014-12-19 13:09:12
    国产一区有黄有色的免费视频| 亚洲国产精品999| 熟女av电影| 狠狠婷婷综合久久久久久88av| 国产黄频视频在线观看| 观看av在线不卡| 激情五月婷婷亚洲| 99热全是精品| 午夜免费观看性视频| 你懂的网址亚洲精品在线观看| 大香蕉久久成人网| 国产日韩一区二区三区精品不卡| 一边摸一边做爽爽视频免费| av视频免费观看在线观看| 精品国产一区二区三区久久久樱花| 亚洲熟女精品中文字幕| 夫妻午夜视频| 亚洲五月色婷婷综合| 91成人精品电影| 国产成人精品福利久久| 日韩精品有码人妻一区| 国产片内射在线| 一本大道久久a久久精品| 午夜精品国产一区二区电影| 看十八女毛片水多多多| 成人国产麻豆网| 久久99热这里只频精品6学生| 麻豆乱淫一区二区| 久久亚洲国产成人精品v| 精品人妻偷拍中文字幕| 国产高清国产精品国产三级| 一边摸一边做爽爽视频免费| 久久精品国产亚洲av涩爱| 午夜日本视频在线| 亚洲欧洲日产国产| 三级国产精品片| 久久久精品免费免费高清| 亚洲国产精品一区二区三区在线| 国产亚洲午夜精品一区二区久久| 亚洲精品第二区| 午夜福利乱码中文字幕| 亚洲精品日韩在线中文字幕| 久久综合国产亚洲精品| 女人被躁到高潮嗷嗷叫费观| 久久久久久久久久久免费av| 欧美bdsm另类| 卡戴珊不雅视频在线播放| 国产精品久久久久久久久免| 国产精品.久久久| 五月伊人婷婷丁香| 黄色怎么调成土黄色| 麻豆乱淫一区二区| 九草在线视频观看| 美女国产视频在线观看| 欧美国产精品va在线观看不卡| 国产亚洲av片在线观看秒播厂| 色94色欧美一区二区| 免费播放大片免费观看视频在线观看| 黄色怎么调成土黄色| 999精品在线视频| 亚洲在久久综合| 国产黄频视频在线观看| 卡戴珊不雅视频在线播放| 香蕉精品网在线| 亚洲成av片中文字幕在线观看 | 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 五月开心婷婷网| 观看av在线不卡| 青春草亚洲视频在线观看| 狂野欧美激情性bbbbbb| 国产成人午夜福利电影在线观看| 汤姆久久久久久久影院中文字幕| 国产一区二区在线观看av| 国产高清国产精品国产三级| 91精品伊人久久大香线蕉| 亚洲精品国产一区二区精华液| 免费人妻精品一区二区三区视频| 九草在线视频观看| 在线精品无人区一区二区三| 蜜桃在线观看..| kizo精华| 成人免费观看视频高清| 久久这里有精品视频免费| 国产xxxxx性猛交| 久久久精品免费免费高清| 赤兔流量卡办理| 美女脱内裤让男人舔精品视频| 男人操女人黄网站| www.精华液| 男女高潮啪啪啪动态图| 嫩草影院入口| 久久久久久久精品精品| 91午夜精品亚洲一区二区三区| 国产精品一区二区在线观看99| 久久久精品免费免费高清| 日韩欧美一区视频在线观看| 成人午夜精彩视频在线观看| 蜜桃在线观看..| 久热这里只有精品99| 久久av网站| av视频免费观看在线观看| 久久鲁丝午夜福利片| 亚洲欧美精品自产自拍| 你懂的网址亚洲精品在线观看| 1024香蕉在线观看| 国产欧美日韩一区二区三区在线| 在线亚洲精品国产二区图片欧美| 国产精品免费视频内射| 亚洲精品在线美女| 欧美激情 高清一区二区三区| 亚洲av日韩在线播放| 精品一区二区免费观看| 亚洲国产精品一区三区| 国产精品一区二区在线不卡| 国产欧美日韩综合在线一区二区| 男女啪啪激烈高潮av片| 汤姆久久久久久久影院中文字幕| 黄频高清免费视频| 亚洲精品aⅴ在线观看| 七月丁香在线播放| 老司机亚洲免费影院| 久久婷婷青草| 青草久久国产| 高清av免费在线| 久久久久久久精品精品| 男女高潮啪啪啪动态图| 国产精品人妻久久久影院| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩视频高清一区二区三区二| 2021少妇久久久久久久久久久| av卡一久久| 久久久久网色| 婷婷色综合www| 三级国产精品片| 97精品久久久久久久久久精品| 一二三四在线观看免费中文在| 曰老女人黄片| 亚洲图色成人| tube8黄色片| 宅男免费午夜| 欧美精品一区二区大全| 国产精品偷伦视频观看了| 亚洲人成电影观看| 七月丁香在线播放| 久久久精品国产亚洲av高清涩受| 少妇人妻久久综合中文| 美女主播在线视频| 国产成人91sexporn| 91久久精品国产一区二区三区| 夜夜骑夜夜射夜夜干| av线在线观看网站| 久久亚洲国产成人精品v| 亚洲精品国产av蜜桃| 国产在视频线精品| 国产免费一区二区三区四区乱码| 精品少妇一区二区三区视频日本电影 | 婷婷成人精品国产| 美女脱内裤让男人舔精品视频| 亚洲精品国产一区二区精华液| 国产精品国产三级国产专区5o| 国产午夜精品一二区理论片| 欧美精品一区二区大全| 国产成人精品婷婷| 亚洲精品国产一区二区精华液| 少妇的丰满在线观看| 久久久久国产一级毛片高清牌| 青青草视频在线视频观看| www日本在线高清视频| 精品国产超薄肉色丝袜足j| 免费女性裸体啪啪无遮挡网站| 亚洲av日韩在线播放| 亚洲av成人精品一二三区| www.精华液| 日韩av免费高清视频| 午夜久久久在线观看| 亚洲精品国产色婷婷电影| 午夜老司机福利剧场| 99久久精品国产国产毛片| 国产精品国产三级国产专区5o| 日本欧美视频一区| 国产精品二区激情视频| 亚洲欧美精品自产自拍| 日本vs欧美在线观看视频| 在线观看免费高清a一片| 啦啦啦啦在线视频资源| 国产又爽黄色视频| 国产又爽黄色视频| 国产xxxxx性猛交| 国产成人精品福利久久| 性高湖久久久久久久久免费观看| 91午夜精品亚洲一区二区三区| 大码成人一级视频| 一区二区三区四区激情视频| 久久国产亚洲av麻豆专区| 国产综合精华液| 秋霞在线观看毛片| 国产又爽黄色视频| 亚洲av国产av综合av卡| 激情五月婷婷亚洲| 9热在线视频观看99| 久久久久久久国产电影| 乱人伦中国视频| 国产野战对白在线观看| 波野结衣二区三区在线| 精品国产一区二区久久| 亚洲婷婷狠狠爱综合网| 欧美成人精品欧美一级黄| 久久久欧美国产精品| 欧美在线黄色| 18+在线观看网站| 熟女电影av网| 精品国产乱码久久久久久小说| 最近中文字幕高清免费大全6| 欧美成人精品欧美一级黄| 大陆偷拍与自拍| 久久久国产一区二区| 欧美日韩av久久| 日日啪夜夜爽| 日韩欧美一区视频在线观看| 女人久久www免费人成看片| 亚洲av欧美aⅴ国产| 日本av免费视频播放| 丰满迷人的少妇在线观看| 在线观看免费高清a一片| 自线自在国产av| 90打野战视频偷拍视频| 人人澡人人妻人| 99久久人妻综合| 你懂的网址亚洲精品在线观看| 大码成人一级视频| 咕卡用的链子| 午夜福利影视在线免费观看| 一级黄片播放器| 99久久精品国产国产毛片| 国产亚洲一区二区精品| 极品少妇高潮喷水抽搐| 日韩精品免费视频一区二区三区| 两性夫妻黄色片| 欧美人与性动交α欧美软件| 性高湖久久久久久久久免费观看| 色视频在线一区二区三区| 国产日韩欧美亚洲二区| 亚洲av日韩在线播放| 日韩中文字幕视频在线看片| 亚洲精品av麻豆狂野| 亚洲国产欧美网| 激情视频va一区二区三区| 久久精品久久久久久久性| 久久精品aⅴ一区二区三区四区 | 看免费成人av毛片| 国产一区二区在线观看av| 麻豆精品久久久久久蜜桃| 一级毛片黄色毛片免费观看视频| 国产成人免费无遮挡视频| 午夜福利乱码中文字幕| 免费高清在线观看视频在线观看| 欧美日韩视频高清一区二区三区二| 香蕉精品网在线| 日韩视频在线欧美| 不卡视频在线观看欧美| 午夜av观看不卡| 亚洲av.av天堂| 亚洲综合精品二区| 一个人免费看片子| 黄频高清免费视频| 各种免费的搞黄视频| 大片免费播放器 马上看| 在线亚洲精品国产二区图片欧美| 国产高清不卡午夜福利| 可以免费在线观看a视频的电影网站 | 日日啪夜夜爽| 少妇的逼水好多| 制服丝袜香蕉在线| 亚洲精品乱久久久久久| 最近中文字幕高清免费大全6| 少妇熟女欧美另类| 看十八女毛片水多多多| 国产av国产精品国产| 日韩在线高清观看一区二区三区| 多毛熟女@视频| 久久久久久人妻| 精品国产露脸久久av麻豆| 又粗又硬又长又爽又黄的视频| 日韩中文字幕欧美一区二区 | 精品人妻偷拍中文字幕| 亚洲美女视频黄频| 亚洲 欧美一区二区三区| 丝瓜视频免费看黄片| 男女午夜视频在线观看| 午夜福利视频在线观看免费| 最黄视频免费看| 少妇 在线观看| 你懂的网址亚洲精品在线观看| 日日爽夜夜爽网站| 国产欧美亚洲国产| 天美传媒精品一区二区| 中文字幕人妻丝袜一区二区 | 97人妻天天添夜夜摸| 国产有黄有色有爽视频| 黄片播放在线免费| 久久精品国产鲁丝片午夜精品| 亚洲一码二码三码区别大吗| 国产亚洲欧美精品永久| 国产成人免费无遮挡视频| 亚洲五月色婷婷综合| 高清在线视频一区二区三区| 熟女av电影| 老鸭窝网址在线观看| 色吧在线观看| 下体分泌物呈黄色| 尾随美女入室| 少妇人妻 视频| 国产精品嫩草影院av在线观看| 黄色 视频免费看| 黑丝袜美女国产一区| 亚洲一级一片aⅴ在线观看| 啦啦啦在线免费观看视频4| 晚上一个人看的免费电影| 日韩制服丝袜自拍偷拍| 久久这里有精品视频免费| 亚洲精品自拍成人| 2018国产大陆天天弄谢| 亚洲精品在线美女| 国产成人精品久久久久久| 卡戴珊不雅视频在线播放| 国产伦理片在线播放av一区| 中国三级夫妇交换| 少妇熟女欧美另类| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 一区二区三区四区激情视频| 一边摸一边做爽爽视频免费| 考比视频在线观看| 国产福利在线免费观看视频| 两个人看的免费小视频| 久久久久国产精品人妻一区二区| videos熟女内射| 视频在线观看一区二区三区| 亚洲欧美一区二区三区国产| 人体艺术视频欧美日本| 中国国产av一级| 波多野结衣av一区二区av| 免费高清在线观看视频在线观看| 国产片内射在线| 亚洲欧美一区二区三区国产| 亚洲第一青青草原| 免费大片黄手机在线观看| 久久精品久久久久久噜噜老黄| 国产1区2区3区精品| 欧美成人午夜免费资源| 精品国产露脸久久av麻豆| 午夜日韩欧美国产| 精品国产一区二区三区久久久樱花| 精品少妇内射三级| 99国产综合亚洲精品| 9191精品国产免费久久| 日韩在线高清观看一区二区三区| 亚洲av.av天堂| 天堂俺去俺来也www色官网| 人人妻人人爽人人添夜夜欢视频| 免费大片黄手机在线观看| 美女主播在线视频| 日本av手机在线免费观看| 国产精品99久久99久久久不卡 | 最近的中文字幕免费完整| 久久久久国产网址| 啦啦啦视频在线资源免费观看| 午夜福利一区二区在线看| 日韩不卡一区二区三区视频在线| 美女国产高潮福利片在线看| 国产毛片在线视频| 国产在线免费精品| 日本色播在线视频| 在线观看美女被高潮喷水网站| 免费在线观看视频国产中文字幕亚洲 | 欧美国产精品va在线观看不卡| 国产精品av久久久久免费| 一二三四中文在线观看免费高清| h视频一区二区三区| 国产成人91sexporn| 亚洲五月色婷婷综合| 日本黄色日本黄色录像| 汤姆久久久久久久影院中文字幕| av国产精品久久久久影院| 狂野欧美激情性bbbbbb| 国产伦理片在线播放av一区| 麻豆av在线久日| 黑人欧美特级aaaaaa片| 日本91视频免费播放| 中国三级夫妇交换| 国产精品香港三级国产av潘金莲 | 亚洲第一青青草原| 777米奇影视久久| 久久国内精品自在自线图片| 新久久久久国产一级毛片| 性高湖久久久久久久久免费观看| 国产av码专区亚洲av| 日本猛色少妇xxxxx猛交久久| 久久精品夜色国产| 黄色配什么色好看| 精品午夜福利在线看| 女人高潮潮喷娇喘18禁视频| 大码成人一级视频| 久久毛片免费看一区二区三区| 国产男女内射视频| 精品亚洲成a人片在线观看| 国产精品偷伦视频观看了| 欧美国产精品va在线观看不卡| 一个人免费看片子| 亚洲国产日韩一区二区| 久久精品久久久久久噜噜老黄| 黄片播放在线免费| 一级爰片在线观看| freevideosex欧美| 丝袜美足系列| 赤兔流量卡办理| 国产激情久久老熟女| 另类亚洲欧美激情| 久久午夜综合久久蜜桃| 日韩电影二区| 香蕉国产在线看| a级毛片黄视频| 国产成人一区二区在线| 观看av在线不卡| 男人爽女人下面视频在线观看| 一级爰片在线观看| 亚洲熟女精品中文字幕| 亚洲精品成人av观看孕妇| 亚洲国产色片| 免费大片黄手机在线观看| 中文字幕最新亚洲高清| 国产精品二区激情视频| 久久免费观看电影| 69精品国产乱码久久久| 日韩一区二区三区影片| a级毛片在线看网站| 好男人视频免费观看在线| 午夜福利一区二区在线看| 欧美人与性动交α欧美精品济南到 | 中国三级夫妇交换| 最近2019中文字幕mv第一页| 国产有黄有色有爽视频| 美女xxoo啪啪120秒动态图| 国产又爽黄色视频| 久久久久精品久久久久真实原创| 国产精品99久久99久久久不卡 | 久久久亚洲精品成人影院| 天堂俺去俺来也www色官网| 国产熟女午夜一区二区三区| 国产乱来视频区| 视频在线观看一区二区三区| 欧美国产精品va在线观看不卡| 岛国毛片在线播放| 国产精品欧美亚洲77777| 老汉色∧v一级毛片| 精品国产露脸久久av麻豆| 人人妻人人添人人爽欧美一区卜| 男人添女人高潮全过程视频| 免费在线观看黄色视频的| 秋霞在线观看毛片| 国产成人精品久久久久久| 在线观看免费高清a一片| 精品人妻熟女毛片av久久网站| 色哟哟·www| 国产毛片在线视频| 日本av手机在线免费观看| 亚洲男人天堂网一区| 国产精品女同一区二区软件| 久久久精品免费免费高清| 精品一区二区三区四区五区乱码 | 亚洲第一av免费看| 街头女战士在线观看网站| 国产片特级美女逼逼视频| 波野结衣二区三区在线| 国产成人av激情在线播放| 日韩电影二区| 高清av免费在线| 国产伦理片在线播放av一区| av国产久精品久网站免费入址| 国产高清不卡午夜福利| 王馨瑶露胸无遮挡在线观看| 岛国毛片在线播放| 天美传媒精品一区二区| 国产免费又黄又爽又色| 在线观看免费日韩欧美大片| 一本久久精品| 深夜精品福利| 啦啦啦在线观看免费高清www| 日日啪夜夜爽| av免费观看日本| 建设人人有责人人尽责人人享有的| 一级片免费观看大全| www日本在线高清视频| 国产黄色免费在线视频| av线在线观看网站| 日韩三级伦理在线观看| 久久久久久久亚洲中文字幕| 色吧在线观看| av在线老鸭窝| 久久99精品国语久久久| 夜夜骑夜夜射夜夜干| 99re6热这里在线精品视频| 男女国产视频网站| 亚洲三级黄色毛片| 亚洲精品乱久久久久久| 国产亚洲av片在线观看秒播厂| 日韩伦理黄色片| 国产无遮挡羞羞视频在线观看| 街头女战士在线观看网站| 晚上一个人看的免费电影| 久久久久视频综合| 午夜福利,免费看| 天美传媒精品一区二区| 国产一区二区三区综合在线观看| 国产日韩一区二区三区精品不卡| 亚洲视频免费观看视频| 七月丁香在线播放| 在线免费观看不下载黄p国产| 久久女婷五月综合色啪小说| 啦啦啦中文免费视频观看日本| 国产精品国产av在线观看| 国产亚洲欧美精品永久| 精品一品国产午夜福利视频| 国产男人的电影天堂91| 美女视频免费永久观看网站| 如日韩欧美国产精品一区二区三区| 狂野欧美激情性bbbbbb| 26uuu在线亚洲综合色| 亚洲第一av免费看| 91在线精品国自产拍蜜月| 99国产综合亚洲精品| 热99国产精品久久久久久7| 男女下面插进去视频免费观看| 午夜福利影视在线免费观看| 成人免费观看视频高清| 国产亚洲最大av| 久久国产精品男人的天堂亚洲| 日本wwww免费看| 国产精品蜜桃在线观看| 这个男人来自地球电影免费观看 | 日韩欧美精品免费久久| 少妇的逼水好多| 久久热在线av| 久久久久国产网址| 午夜91福利影院| 色网站视频免费| 老司机亚洲免费影院| 人成视频在线观看免费观看| 亚洲欧美中文字幕日韩二区| 免费观看在线日韩| 一区二区三区四区激情视频| 国产一区有黄有色的免费视频| 日韩精品有码人妻一区| 久久久久久免费高清国产稀缺| 最新的欧美精品一区二区| av女优亚洲男人天堂| 国产女主播在线喷水免费视频网站| 欧美日韩av久久| 国产精品国产三级专区第一集| 久久精品久久久久久久性| 日本色播在线视频| 三上悠亚av全集在线观看| 久久久国产一区二区| 成人免费观看视频高清| 亚洲欧美精品综合一区二区三区 | 夜夜骑夜夜射夜夜干| freevideosex欧美| tube8黄色片| 亚洲 欧美一区二区三区| 一级片免费观看大全| 日日啪夜夜爽| 天堂8中文在线网| 午夜福利网站1000一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产97色在线日韩免费| 成年动漫av网址| 中文字幕最新亚洲高清| 捣出白浆h1v1| 亚洲 欧美一区二区三区| 丝袜脚勾引网站| 国产精品蜜桃在线观看| 最黄视频免费看| 熟妇人妻不卡中文字幕| 国产国语露脸激情在线看| 国产野战对白在线观看| 国产精品三级大全| 黄色毛片三级朝国网站| 久久久亚洲精品成人影院| 极品少妇高潮喷水抽搐| 天美传媒精品一区二区| 老鸭窝网址在线观看| 女人高潮潮喷娇喘18禁视频| 午夜免费鲁丝| 精品人妻偷拍中文字幕| 在线观看人妻少妇| 国精品久久久久久国模美| 国产黄色免费在线视频| 日本欧美视频一区| 综合色丁香网| 人体艺术视频欧美日本| 日韩 亚洲 欧美在线| 波多野结衣一区麻豆| 亚洲av在线观看美女高潮| 一本色道久久久久久精品综合| 亚洲三级黄色毛片| 亚洲av电影在线观看一区二区三区| 91精品三级在线观看| 午夜福利视频精品| 中文字幕人妻丝袜制服| 国产老妇伦熟女老妇高清| 婷婷色麻豆天堂久久| 在线观看www视频免费| 免费在线观看黄色视频的| 久久国产精品男人的天堂亚洲|