• <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
    国产伦一二天堂av在线观看| 在线观看午夜福利视频| 美女免费视频网站| 国产色婷婷99| 国产成人aa在线观看| 国产淫片久久久久久久久 | 亚洲五月婷婷丁香| 国产99白浆流出| 校园春色视频在线观看| 69人妻影院| 亚洲美女视频黄频| 天天一区二区日本电影三级| 美女高潮的动态| 极品教师在线免费播放| 此物有八面人人有两片| 看片在线看免费视频| 十八禁人妻一区二区| 精品国产三级普通话版| 国产毛片a区久久久久| 99久久综合精品五月天人人| 国内毛片毛片毛片毛片毛片| 欧美日本亚洲视频在线播放| 国产亚洲精品av在线| 中文字幕精品亚洲无线码一区| АⅤ资源中文在线天堂| 成人特级黄色片久久久久久久| 精品人妻一区二区三区麻豆 | 欧美日本视频| 一进一出抽搐动态| 国产亚洲欧美在线一区二区| 国产乱人视频| 亚洲一区二区三区不卡视频| 非洲黑人性xxxx精品又粗又长| 亚洲国产欧美人成| 国产精品久久久久久精品电影| eeuss影院久久| 熟女电影av网| 人妻久久中文字幕网| 欧美xxxx黑人xx丫x性爽| 亚洲精品亚洲一区二区| 午夜视频国产福利| 国产麻豆成人av免费视频| 女生性感内裤真人,穿戴方法视频| 精品电影一区二区在线| 亚洲中文日韩欧美视频| 精品无人区乱码1区二区| 乱人视频在线观看| 搡老岳熟女国产| www日本在线高清视频| 在线播放国产精品三级| 成年免费大片在线观看| 男人的好看免费观看在线视频| 国产av在哪里看| 国产在视频线在精品| 日本一二三区视频观看| 深夜精品福利| 国产乱人视频| 在线十欧美十亚洲十日本专区| avwww免费| 高清在线国产一区| 一个人观看的视频www高清免费观看| av欧美777| 小说图片视频综合网站| 久久久久免费精品人妻一区二区| 国产v大片淫在线免费观看| 久久午夜亚洲精品久久| 亚洲国产精品成人综合色| 午夜亚洲福利在线播放| 欧美黑人欧美精品刺激| 国产视频一区二区在线看| 欧美一级a爱片免费观看看| 三级毛片av免费| 亚洲七黄色美女视频| 欧美大码av| 国产男靠女视频免费网站| 午夜福利成人在线免费观看| 婷婷六月久久综合丁香| 国产激情偷乱视频一区二区| 欧美zozozo另类| 久久久久久九九精品二区国产| 久久人妻av系列| 欧美黑人巨大hd| 中文字幕av成人在线电影| 女警被强在线播放| www.色视频.com| 国内精品久久久久精免费| 亚洲av一区综合| 国产精品三级大全| 亚洲精品美女久久久久99蜜臀| 国产精品久久久人人做人人爽| 午夜福利高清视频| 欧美日韩瑟瑟在线播放| 国产精品一区二区免费欧美| xxxwww97欧美| www.www免费av| 亚洲无线在线观看| avwww免费| 好男人电影高清在线观看| 久久久久国产精品人妻aⅴ院| 国产色爽女视频免费观看| 欧美黄色淫秽网站| 高清毛片免费观看视频网站| 90打野战视频偷拍视频| 日本 av在线| 男人舔奶头视频| 欧美不卡视频在线免费观看| 日日夜夜操网爽| 91麻豆av在线| 一边摸一边抽搐一进一小说| 岛国在线免费视频观看| 亚洲性夜色夜夜综合| 亚洲狠狠婷婷综合久久图片| 久久亚洲精品不卡| 国产一区二区三区视频了| 岛国在线观看网站| 国产久久久一区二区三区| 日韩有码中文字幕| 免费看十八禁软件| 男女之事视频高清在线观看| АⅤ资源中文在线天堂| 久久精品国产清高在天天线| 午夜福利18| 亚洲中文字幕一区二区三区有码在线看| 香蕉久久夜色| 婷婷亚洲欧美| 亚洲在线自拍视频| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区精品| 亚洲av熟女| 久久久国产精品麻豆| 亚洲激情在线av| 午夜影院日韩av| 波多野结衣巨乳人妻| 亚洲成a人片在线一区二区| 亚洲七黄色美女视频| 变态另类丝袜制服| 亚洲成人中文字幕在线播放| 香蕉久久夜色| av专区在线播放| 国产亚洲精品综合一区在线观看| 精品日产1卡2卡| 男人和女人高潮做爰伦理| 69人妻影院| 亚洲内射少妇av| 99久久精品热视频| 最好的美女福利视频网| 1000部很黄的大片| 成年版毛片免费区| 久久精品国产综合久久久| 人妻丰满熟妇av一区二区三区| 99热这里只有精品一区| 很黄的视频免费| 成年人黄色毛片网站| 99国产精品一区二区蜜桃av| 亚洲无线在线观看| 中文字幕人妻丝袜一区二区| 国内揄拍国产精品人妻在线| 国产精品影院久久| 999久久久精品免费观看国产| 午夜福利视频1000在线观看| 天美传媒精品一区二区| 国产精品 国内视频| 麻豆国产av国片精品| 国产精品女同一区二区软件 | 超碰av人人做人人爽久久 | 国产精品久久久久久人妻精品电影| 99riav亚洲国产免费| 女人被狂操c到高潮| 制服人妻中文乱码| 韩国av一区二区三区四区| 欧美日韩精品网址| 夜夜看夜夜爽夜夜摸| 麻豆久久精品国产亚洲av| 精品人妻1区二区| 成人国产一区最新在线观看| 国产高清视频在线播放一区| 亚洲一区二区三区色噜噜| 国产精品亚洲一级av第二区| 欧美日韩瑟瑟在线播放| 欧美+日韩+精品| 国语自产精品视频在线第100页| 国内精品美女久久久久久| 国产一区二区在线观看日韩 | 欧美成人a在线观看| 午夜激情福利司机影院| 成熟少妇高潮喷水视频| 熟女少妇亚洲综合色aaa.| 制服人妻中文乱码| 波多野结衣巨乳人妻| 又黄又爽又免费观看的视频| 国产黄片美女视频| 人妻夜夜爽99麻豆av| 日本五十路高清| 午夜福利欧美成人| 国产高清视频在线观看网站| 色吧在线观看| 日本精品一区二区三区蜜桃| 丰满乱子伦码专区| 国产精品香港三级国产av潘金莲| 国产99白浆流出| 搡老熟女国产l中国老女人| 美女高潮的动态| 深爱激情五月婷婷| 国产精品久久久久久亚洲av鲁大| 小蜜桃在线观看免费完整版高清| 成人av在线播放网站| 亚洲熟妇熟女久久| 欧美三级亚洲精品| 一本精品99久久精品77| 精品一区二区三区视频在线 | 欧美zozozo另类| 精品国产三级普通话版| 午夜福利18| 久久人人精品亚洲av| 亚洲乱码一区二区免费版| 精品久久久久久,| 搡老岳熟女国产| АⅤ资源中文在线天堂| 偷拍熟女少妇极品色| 午夜福利视频1000在线观看| 免费在线观看亚洲国产| 午夜福利成人在线免费观看| 国产色婷婷99| 99国产精品一区二区三区| 三级男女做爰猛烈吃奶摸视频| 日韩免费av在线播放| 国内精品美女久久久久久| 成人国产一区最新在线观看| 波多野结衣高清无吗| 午夜免费成人在线视频| 午夜福利欧美成人| 黄片小视频在线播放| 国产精品三级大全| 性色av乱码一区二区三区2| 欧美黑人巨大hd| 亚洲国产色片| 色综合欧美亚洲国产小说| 欧美中文日本在线观看视频| 亚洲在线观看片| 亚洲国产日韩欧美精品在线观看 | 3wmmmm亚洲av在线观看| 久久亚洲精品不卡| 高清日韩中文字幕在线| www.999成人在线观看| 少妇的逼好多水| 国产伦人伦偷精品视频| 亚洲精品粉嫩美女一区| 免费观看人在逋| 久久天躁狠狠躁夜夜2o2o| 久久草成人影院| 亚洲av成人精品一区久久| 亚洲av五月六月丁香网| 麻豆久久精品国产亚洲av| 亚洲精品在线观看二区| 大型黄色视频在线免费观看| 无遮挡黄片免费观看| 国模一区二区三区四区视频| 99精品久久久久人妻精品| 国产亚洲欧美98| 黄色丝袜av网址大全| 一区二区三区国产精品乱码| 无限看片的www在线观看| 国产毛片a区久久久久| 日韩中文字幕欧美一区二区| 精品一区二区三区人妻视频| 又粗又爽又猛毛片免费看| 国产成人系列免费观看| 色综合亚洲欧美另类图片| 长腿黑丝高跟| 最近视频中文字幕2019在线8| 久久久久久久精品吃奶| 一本久久中文字幕| 狂野欧美白嫩少妇大欣赏| 亚洲一区高清亚洲精品| 91麻豆av在线| 国产在线精品亚洲第一网站| 小蜜桃在线观看免费完整版高清| 亚洲国产精品合色在线| 亚洲电影在线观看av| 91字幕亚洲| 五月伊人婷婷丁香| 亚洲,欧美精品.| 亚洲欧美日韩卡通动漫| 国产成人影院久久av| av片东京热男人的天堂| 脱女人内裤的视频| 女同久久另类99精品国产91| 美女免费视频网站| 淫秽高清视频在线观看| 两个人视频免费观看高清| 午夜精品久久久久久毛片777| 成人鲁丝片一二三区免费| 久久久久精品国产欧美久久久| 在线播放无遮挡| 99视频精品全部免费 在线| 人人妻人人澡欧美一区二区| 美女黄网站色视频| 男女之事视频高清在线观看| 精品久久久久久久末码| 亚洲成a人片在线一区二区| 国内久久婷婷六月综合欲色啪| 国产精品野战在线观看| 国产精品国产高清国产av| 精品一区二区三区人妻视频| 亚洲国产精品sss在线观看| 岛国视频午夜一区免费看| 国产一区二区激情短视频| 一夜夜www| 久久午夜亚洲精品久久| 欧美丝袜亚洲另类 | 又黄又粗又硬又大视频| 精品国产超薄肉色丝袜足j| 亚洲欧美日韩无卡精品| 美女高潮喷水抽搐中文字幕| 一个人免费在线观看电影| 欧美色视频一区免费| 人妻久久中文字幕网| 日日摸夜夜添夜夜添小说| 色老头精品视频在线观看| 亚洲成人久久性| 国产高清激情床上av| 一个人免费在线观看的高清视频| 99久久无色码亚洲精品果冻| 亚洲人与动物交配视频| 欧美日韩亚洲国产一区二区在线观看| 中文资源天堂在线| 国语自产精品视频在线第100页| 神马国产精品三级电影在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品女同一区二区软件 | 午夜久久久久精精品| 欧美绝顶高潮抽搐喷水| 12—13女人毛片做爰片一| 一a级毛片在线观看| 极品教师在线免费播放| 最近在线观看免费完整版| 久久伊人香网站| 国产熟女xx| 两个人的视频大全免费| 久久精品国产清高在天天线| 在线观看午夜福利视频| 精品久久久久久久久久久久久| 亚洲av二区三区四区| 精品国产三级普通话版| 亚洲黑人精品在线| 日本撒尿小便嘘嘘汇集6| 丁香六月欧美| 一进一出抽搐gif免费好疼| 88av欧美| 日韩av在线大香蕉| 老司机福利观看| 国产免费男女视频| 久久精品国产自在天天线| 男女之事视频高清在线观看| 桃色一区二区三区在线观看| 无限看片的www在线观看| 日韩人妻高清精品专区| 在线免费观看的www视频| 成人午夜高清在线视频| 亚洲人成伊人成综合网2020| 久久人妻av系列| 国产黄片美女视频| 精品久久久久久,| 国产精品免费一区二区三区在线| 日日干狠狠操夜夜爽| 少妇的丰满在线观看| 搞女人的毛片| 久久精品国产自在天天线| 无限看片的www在线观看| 成人三级黄色视频| 免费看十八禁软件| 国产探花在线观看一区二区| 国产精品影院久久| 亚洲精品影视一区二区三区av| 亚洲在线自拍视频| 午夜免费成人在线视频| 女人十人毛片免费观看3o分钟| 国产成人av激情在线播放| 男女午夜视频在线观看| 欧美黑人巨大hd| 日日摸夜夜添夜夜添小说| 国产真人三级小视频在线观看| 国产成人影院久久av| netflix在线观看网站| 一区福利在线观看| 亚洲人成电影免费在线| 又黄又爽又免费观看的视频| 天堂动漫精品| 精品无人区乱码1区二区| 国产视频一区二区在线看| 中文字幕熟女人妻在线| 国产97色在线日韩免费| 亚洲一区二区三区不卡视频| 精华霜和精华液先用哪个| 日韩亚洲欧美综合| 亚洲 欧美 日韩 在线 免费| 高清毛片免费观看视频网站| 亚洲中文字幕一区二区三区有码在线看| 一本久久中文字幕| 国内精品美女久久久久久| 男人和女人高潮做爰伦理| 国产99白浆流出| 久久久色成人| 久久久久久九九精品二区国产| 麻豆国产av国片精品| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩福利视频一区二区| 国产精品影院久久| 又爽又黄无遮挡网站| 亚洲国产精品sss在线观看| 国产成人欧美在线观看| 久久精品国产自在天天线| 精品久久久久久成人av| 久久精品国产亚洲av香蕉五月| svipshipincom国产片| 黄色女人牲交| 国产又黄又爽又无遮挡在线| 老汉色av国产亚洲站长工具| 一区二区三区国产精品乱码| 国内毛片毛片毛片毛片毛片| 久久久久国产精品人妻aⅴ院| 69av精品久久久久久| 亚洲色图av天堂| 天天躁日日操中文字幕| 熟女电影av网| 亚洲国产色片| 精品人妻1区二区| 欧美黄色淫秽网站| 男女下面进入的视频免费午夜| 日韩欧美在线二视频| 黄色女人牲交| 欧美一级a爱片免费观看看| 国产成年人精品一区二区| 免费观看精品视频网站| 国产精品爽爽va在线观看网站| 国产熟女xx| 搡老妇女老女人老熟妇| 老熟妇乱子伦视频在线观看| 午夜老司机福利剧场| 成人特级黄色片久久久久久久| 日本一二三区视频观看| 久久久久久九九精品二区国产| 国产成人av教育| 午夜影院日韩av| 成年女人毛片免费观看观看9| 亚洲成人中文字幕在线播放| 法律面前人人平等表现在哪些方面| 搡女人真爽免费视频火全软件 | 久久国产精品影院| 久久性视频一级片| 一进一出抽搐gif免费好疼| 中文字幕精品亚洲无线码一区| 最好的美女福利视频网| 精品久久久久久成人av| 黄片小视频在线播放| 免费观看精品视频网站| 国产精品精品国产色婷婷| 国产亚洲精品一区二区www| 一个人免费在线观看电影| 少妇高潮的动态图| www.www免费av| 精品一区二区三区人妻视频| 在线免费观看不下载黄p国产 | 免费av不卡在线播放| 久久久久久久午夜电影| 日韩大尺度精品在线看网址| 丰满人妻一区二区三区视频av | 亚洲色图av天堂| 在线观看美女被高潮喷水网站 | 亚洲成a人片在线一区二区| 麻豆久久精品国产亚洲av| 欧美xxxx黑人xx丫x性爽| 国产成+人综合+亚洲专区| 国产精品日韩av在线免费观看| 亚洲天堂国产精品一区在线| 一区福利在线观看| 婷婷丁香在线五月| 国产av在哪里看| 欧美日韩一级在线毛片| 国产99白浆流出| 18美女黄网站色大片免费观看| 久久精品国产亚洲av涩爱 | 亚洲电影在线观看av| 99久久无色码亚洲精品果冻| 亚洲欧美日韩东京热| 最近视频中文字幕2019在线8| 成人一区二区视频在线观看| 亚洲,欧美精品.| 亚洲精品乱码久久久v下载方式 | aaaaa片日本免费| 午夜福利在线在线| 91麻豆av在线| 亚洲,欧美精品.| 国产伦精品一区二区三区四那| 噜噜噜噜噜久久久久久91| 特级一级黄色大片| www.www免费av| 久久久久久九九精品二区国产| 夜夜夜夜夜久久久久| 18禁黄网站禁片免费观看直播| 久久久久精品国产欧美久久久| 亚洲精品美女久久久久99蜜臀| 国产亚洲精品综合一区在线观看| 亚洲aⅴ乱码一区二区在线播放| 全区人妻精品视频| 99热6这里只有精品| 欧美av亚洲av综合av国产av| 久久久久精品国产欧美久久久| 国产国拍精品亚洲av在线观看 | 亚洲18禁久久av| 9191精品国产免费久久| 岛国在线免费视频观看| 动漫黄色视频在线观看| 热99re8久久精品国产| 啦啦啦免费观看视频1| 欧美大码av| 国产成人福利小说| 日本熟妇午夜| 桃红色精品国产亚洲av| 欧美又色又爽又黄视频| 99热这里只有精品一区| 日本精品一区二区三区蜜桃| 亚洲无线观看免费| 在线观看午夜福利视频| 嫩草影视91久久| 欧美一级毛片孕妇| 亚洲av免费在线观看| netflix在线观看网站| 夜夜爽天天搞| 欧美一区二区精品小视频在线| 成年女人毛片免费观看观看9| 久久人人精品亚洲av| 亚洲不卡免费看| 变态另类丝袜制服| 亚洲欧美日韩东京热| 美女cb高潮喷水在线观看| 久久欧美精品欧美久久欧美| 亚洲久久久久久中文字幕| 国产av不卡久久| 国产伦在线观看视频一区| 最后的刺客免费高清国语| 伊人久久大香线蕉亚洲五| 亚洲avbb在线观看| 亚洲熟妇熟女久久| 成人高潮视频无遮挡免费网站| 亚洲成av人片免费观看| 无遮挡黄片免费观看| 欧美大码av| 亚洲狠狠婷婷综合久久图片| 亚洲色图av天堂| 欧美日韩一级在线毛片| 国产一区二区在线观看日韩 | 欧美bdsm另类| 啦啦啦韩国在线观看视频| 亚洲av第一区精品v没综合| 午夜福利在线观看免费完整高清在 | 国产黄a三级三级三级人| 国产精品99久久久久久久久| 国产精华一区二区三区| 久久久久久人人人人人| 成人鲁丝片一二三区免费| 丰满人妻一区二区三区视频av | 欧美黑人欧美精品刺激| 国产精品99久久久久久久久| 国产亚洲精品久久久com| 国产午夜精品久久久久久一区二区三区 | 欧美最黄视频在线播放免费| 性色av乱码一区二区三区2| 动漫黄色视频在线观看| 又爽又黄无遮挡网站| 美女高潮喷水抽搐中文字幕| 欧美日韩国产亚洲二区| 青草久久国产| 午夜影院日韩av| 国产高清三级在线| 欧美一区二区亚洲| 午夜两性在线视频| 欧美丝袜亚洲另类 | 好男人电影高清在线观看| 成熟少妇高潮喷水视频| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av涩爱 | 欧美极品一区二区三区四区| 久久精品影院6| 成人国产综合亚洲| 一进一出抽搐动态| 国产免费av片在线观看野外av| 18+在线观看网站| 国产成人影院久久av| 村上凉子中文字幕在线| 亚洲在线观看片| 成人午夜高清在线视频| 精品人妻偷拍中文字幕| 亚洲第一电影网av| 听说在线观看完整版免费高清| 国语自产精品视频在线第100页| 在线免费观看不下载黄p国产 | 亚洲 国产 在线| 高清在线国产一区| 欧美性猛交╳xxx乱大交人| 精品国产亚洲在线| 99久久综合精品五月天人人| 久久精品影院6| 免费无遮挡裸体视频| 午夜精品久久久久久毛片777| 男女下面进入的视频免费午夜| 啦啦啦观看免费观看视频高清| 69av精品久久久久久| 又粗又爽又猛毛片免费看| 国产高清videossex| 淫秽高清视频在线观看| 亚洲中文字幕日韩| 欧美黄色淫秽网站| 亚洲18禁久久av| 老司机在亚洲福利影院|