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

    木勺生產(chǎn)線機(jī)械臂的動(dòng)靜特性及結(jié)構(gòu)優(yōu)化

    2024-04-02 08:34:00翟慶磊曹文豪張慶功王錦添
    林產(chǎn)工業(yè) 2024年3期
    關(guān)鍵詞:型板固有頻率生產(chǎn)線

    趙 輝 翟慶磊 曹文豪 張慶功 王錦添 李 明

    (東北林業(yè)大學(xué)機(jī)電工程學(xué)院,黑龍江 哈爾濱 150040)

    隨著2023 年《商務(wù)領(lǐng)域經(jīng)營(yíng)者使用、報(bào)告一次性塑料制品管理辦法》的頒布,國(guó)家對(duì)于塑料制品的限制越來(lái)越嚴(yán),同時(shí)鼓勵(lì)餐飲行業(yè)使用木制品餐具代替一次性塑料制品,木制品餐具企業(yè)的訂單數(shù)量不斷增加。然而,目前中小型企業(yè)木勺生產(chǎn)線仍采用人工上下料的方式,工人勞動(dòng)強(qiáng)度大,且工作效率較低。此外,由于木質(zhì)餐具需要在高溫下沖壓制作,這使得人工上下料的作業(yè)方式非常危險(xiǎn)。作為自動(dòng)化設(shè)備,龍門式機(jī)械臂被廣泛用于工廠車間,代替人工完成重復(fù)、枯燥且頻繁的工作[1-4]。

    本文擬根據(jù)工廠實(shí)際使用的沖壓機(jī)(圖1),以現(xiàn)有的龍門式機(jī)械臂為基礎(chǔ)(圖2),通過對(duì)其結(jié)構(gòu)進(jìn)行改進(jìn),設(shè)計(jì)可用于木勺生產(chǎn)線的機(jī)械臂。為確保機(jī)械臂順利實(shí)現(xiàn)上下料功能,要求其在上下料過程中左右晃動(dòng)幅度保持在2 mm以內(nèi)。鑒于此,本文使用建模軟件對(duì)初步設(shè)計(jì)的生產(chǎn)線機(jī)械臂進(jìn)行三維建模處理,并運(yùn)用建模軟件與有限元分析軟件的協(xié)同仿真技術(shù)對(duì)機(jī)械臂進(jìn)行靜態(tài)性能分析,確定其最不穩(wěn)定的工況,并選取此工況進(jìn)行動(dòng)態(tài)性能分析,運(yùn)用靈敏度分析機(jī)械臂的設(shè)計(jì)變量,篩選對(duì)其動(dòng)靜性能影響較大的設(shè)計(jì)變量。采用拉丁超立方抽樣對(duì)選取的設(shè)計(jì)變量進(jìn)行采樣,分別運(yùn)用Kriging模型與NPR模型構(gòu)建近似模型,并采用MOGA法(多目標(biāo)遺傳算法)對(duì)設(shè)計(jì)的目標(biāo)函數(shù)進(jìn)行求解,實(shí)現(xiàn)對(duì)生產(chǎn)線機(jī)械臂的優(yōu)化。

    圖1 沖壓機(jī)結(jié)構(gòu)Fig.1 Stamping machine structure

    圖2 生產(chǎn)線機(jī)械臂結(jié)構(gòu)Fig.2 Structure of the robotic arm on the production line

    1 靜態(tài)分析

    為節(jié)約計(jì)算的時(shí)間,忽略生產(chǎn)線機(jī)械臂模型中的倒角、螺紋等不重要細(xì)節(jié),同時(shí)忽略對(duì)生產(chǎn)線機(jī)械臂動(dòng)靜性能影響不大,且不在本研究結(jié)構(gòu)范圍內(nèi)的驅(qū)動(dòng)電機(jī)、真空吸盤等部件[5-10]。簡(jiǎn)化后的生產(chǎn)線機(jī)械臂結(jié)構(gòu)如圖3 所示。

    圖3 生產(chǎn)線機(jī)械臂簡(jiǎn)化結(jié)構(gòu)Fig.3 Simplified structure of the production line robotic arm

    龍門架采用的材料為6061 鋁合金,地腳采用的材料為HT200 灰鑄鐵,L型板與T型板采用的材料為1Cr13 不銹鋼,水平與豎直移動(dòng)機(jī)構(gòu)采用的為45 號(hào)鋼。材料的具體參數(shù)如表1 所示。

    表1 材料具體參數(shù)Tab.1 Material specific parameters

    根據(jù)生產(chǎn)線機(jī)械臂結(jié)構(gòu)的對(duì)稱性,可對(duì)其工況進(jìn)行適當(dāng)簡(jiǎn)化[3],可分為四種極限工況:

    1)工況一:豎直移動(dòng)機(jī)構(gòu)處于水平移動(dòng)機(jī)構(gòu)邊緣,且自身處于極限上工況;

    2)工況二:豎直移動(dòng)機(jī)構(gòu)處于水平移動(dòng)機(jī)構(gòu)邊緣,且自身處于極限下工況;

    3)工況三:豎直移動(dòng)機(jī)構(gòu)處于水平移動(dòng)機(jī)構(gòu)中心,且自身處于極限上工況;

    4)工況四:豎直移動(dòng)機(jī)構(gòu)處于水平移動(dòng)機(jī)構(gòu)中心,且自身處于極限下工況。

    生產(chǎn)線機(jī)械臂在進(jìn)行上下料作業(yè)時(shí),所受到的外部載荷主要來(lái)源于木勺毛坯與真空吸盤,真空吸盤吸滿木勺毛坯時(shí)的質(zhì)量為11 kg。進(jìn)行網(wǎng)格劃分,同時(shí)以Force方式在T型板對(duì)應(yīng)位置添加所受的力,并且在龍門架地腳位置添加對(duì)應(yīng)的固定支撐。各個(gè)極限工況對(duì)應(yīng)的應(yīng)力如表2 所示。

    表2 各個(gè)極限工況對(duì)應(yīng)的應(yīng)力Tab.2 Corresponding stresses for each extreme working condition

    由表2 可知,當(dāng)豎直移動(dòng)機(jī)構(gòu)處于水平移動(dòng)機(jī)構(gòu)中心,且自身處于極限上工況時(shí),其應(yīng)力值最大,此時(shí)的等效應(yīng)力圖與總變形圖分別如圖4 和圖5 所示。

    圖4 等效應(yīng)力圖Fig.4 Equivalent stress diagram

    圖5 總變形圖Fig.5 Overall deformation diagram

    由圖4 和圖5 可知,生產(chǎn)線機(jī)械臂極限工況下最大應(yīng)力發(fā)生在L型板底部,數(shù)值為1.938 8 MPa,最大變形發(fā)生在T型板,數(shù)值為0.043 3 mm。在提升安全系數(shù)的情況下,生產(chǎn)線機(jī)械臂需要滿足的剛度與強(qiáng)度分別如下:

    式中:fmax為生產(chǎn)線機(jī)械臂可以承受的最大靜變形量,mm;[f]為所選材料的許用撓度,mm;L為T型板的長(zhǎng)度,數(shù)值為196 mm。

    式中:σmax為生產(chǎn)線機(jī)械臂在受到約束情況時(shí)的最大等效應(yīng)力,MPa;[σ]為L(zhǎng)型板所選材料1Cr13 不銹鋼的屈服強(qiáng)度,數(shù)值為245 MPa;1.5 為所選的安全系數(shù)。

    根據(jù)公式(1)和(2)可知,生產(chǎn)線機(jī)械臂在極限工況下仍然符合使用的要求,且安全余量很大。同時(shí),該工況為理論上的最危險(xiǎn)工況,其余工況的最大應(yīng)力值與最大變形均小于該工況。因此,選取該工況下的模型進(jìn)行下一步的動(dòng)態(tài)分析研究。

    2 動(dòng)態(tài)分析

    2.1 模態(tài)分析

    為避免初步設(shè)計(jì)的生產(chǎn)線機(jī)械臂在上下料過程中出現(xiàn)共振問題,擬采用模態(tài)分析法[11-14]。

    在有限元分析軟件中,對(duì)簡(jiǎn)化后的生產(chǎn)線機(jī)械臂進(jìn)行網(wǎng)格劃分,同時(shí),在龍門架地腳位置添加對(duì)應(yīng)的固定支撐[15-16],將分析設(shè)置中的最大模態(tài)階數(shù)調(diào)整成為10階[17],求解后得出的前10階固有頻率如表3所示。

    表3 前10 階固有頻率Tab.3 First 10 natural frequencies

    在振動(dòng)過程中,由于外力變化速度太快,模型來(lái)不及變形或者變形速度達(dá)不到外力變化的速度,導(dǎo)致在固有頻率較低時(shí)模型變形程度較大。因此,本文選取前6階模態(tài)進(jìn)行重點(diǎn)研究[18]。

    2.2 諧響應(yīng)分析

    為了更加精確地分析在施加外部載荷后,生產(chǎn)線機(jī)械臂的振動(dòng)情況,需要對(duì)生產(chǎn)線機(jī)械臂進(jìn)行諧響應(yīng)分析[18]。

    設(shè)置諧響應(yīng)分析的最大頻率為45 Hz,并且以Force方式在T型板對(duì)應(yīng)位置添加所受的外部載荷[18]。圖6為T型板在Y方向的振幅頻率響應(yīng)曲線。

    圖6 T型板在Y方向的振幅頻率響應(yīng)曲線Fig.6 Amplitude frequency response curve of T-shaped plate in the Y direction

    從圖6可知,在1~45 Hz頻率的影響下,生產(chǎn)線機(jī)械臂的T型板在Y方向的位移最大為4.155 5 mm,對(duì)應(yīng)的頻率為34.244 Hz。然而,生產(chǎn)線機(jī)械臂在上下料過程中所允許出現(xiàn)的最大位移為2 mm,因此需要進(jìn)行優(yōu)化設(shè)計(jì)。

    3 靈敏度分析與優(yōu)化

    為實(shí)現(xiàn)對(duì)生產(chǎn)線機(jī)械臂的優(yōu)化,選擇生產(chǎn)線機(jī)械臂的立柱、橫梁、縱軸、L型板以及T型板進(jìn)行研究。同一結(jié)構(gòu)的參數(shù)分布如圖7 所示,立柱、橫梁與縱軸對(duì)應(yīng)的厚度參數(shù)位置如圖8 所示,L型板和T型板的參數(shù)位置分別如圖9 和10 所示。

    圖7 同一結(jié)構(gòu)的參數(shù)分布Fig.7 Parameter distribution of the same structure

    圖8 立柱、橫梁與縱軸對(duì)應(yīng)的參數(shù)分布Fig.8 Parameter distribution corresponding to columns,beams, and longitudinal axis

    圖9 L型板的參數(shù)分布Fig.9 Parameter distribution of L-shaped plate

    圖10 T型板的參數(shù)分布Fig.10 Parameter distribution of T-shaped plate

    利用有限元分析軟件的Parameters Correlation模塊,對(duì)各個(gè)初始設(shè)計(jì)變量進(jìn)行計(jì)算,得出各個(gè)初始設(shè)計(jì)變量對(duì)于第四階固有頻率P19、生產(chǎn)線機(jī)械臂質(zhì)量P20 以及Y方向最大振幅P39 的靈敏度,分別如表4、5、6 所示。

    表4 變量對(duì)第四階固有頻率的靈敏度Tab.4 Sensitivity of variables to the fourth order natural frequency

    表5 變量對(duì)質(zhì)量的靈敏度Tab.5 Sensitivity of variables to mass

    表6 變量對(duì)Y方向最大振幅的靈敏度Tab.6 Sensitivity of variables to the maximum amplitude in the Y direction

    當(dāng)?shù)谒碾A固有頻率與初始設(shè)計(jì)變量成正比時(shí),靈敏度為正值;當(dāng)?shù)谒碾A固有頻率與初始設(shè)計(jì)變量成反比時(shí),靈敏度為負(fù)值。質(zhì)量、Y方向最大振幅與初始設(shè)計(jì)變量的關(guān)系亦是如此。因此,為了盡可能優(yōu)化生產(chǎn)線機(jī)械臂的性能,同時(shí)提高優(yōu)化效率,選取靈敏度絕對(duì)值較大的初始設(shè)計(jì)變量進(jìn)行優(yōu)化[19-21]。

    最終,選取DS_D1、DS_D2、DS_D3、DS_D4、DS_D5、DS_D6、DS_D7、DS_D8、DS_D9、DS_D10、DS_D11、DS_D13、DS_D14、DS_D16、DS_D17、DS_D18作為優(yōu)化的對(duì)象。

    在對(duì)生產(chǎn)線機(jī)械臂進(jìn)行優(yōu)化時(shí),為確保樣本分布均勻,并快速實(shí)現(xiàn)收斂,本文運(yùn)用多維分層的拉丁超立方抽樣,具體抽選樣本點(diǎn)的過程如圖11 所示。

    圖11 拉丁超立方抽樣過程Fig.11 Latin hypercube sampling process

    選取生產(chǎn)線機(jī)械臂設(shè)計(jì)變量作為輸入,以第四階固有頻率、質(zhì)量以及Y方向最大振幅作為輸出,選取樣本數(shù)量為140,隨機(jī)生成器種子數(shù)量為10,通過求解,獲得抽取的樣本以及對(duì)應(yīng)的輸出。

    4 建立近似模型

    在實(shí)際計(jì)算過程中,使用生產(chǎn)線機(jī)械臂直接進(jìn)行有限元優(yōu)化不僅需要反復(fù)調(diào)用模型,效率低,而且容易出現(xiàn)錯(cuò)誤。因此,為節(jié)約計(jì)算時(shí)間,提高優(yōu)化效率,采用近似模型代替有限元分析。

    考慮到生產(chǎn)線機(jī)械臂存在明顯的非線性,運(yùn)用Kriging模型以及NPR模型進(jìn)行有限元分析。Kriging模型存在明顯的平滑效應(yīng),可以為生產(chǎn)線機(jī)械臂進(jìn)行無(wú)偏差預(yù)測(cè)[22],其表達(dá)式為:

    式中:GKriging(D)為是近似模型;gT(D)為全局模擬的多項(xiàng)式計(jì)算函數(shù);β為回歸系數(shù);Z(D)為是局部偏差。

    NPR模型在使用時(shí),不需要提前規(guī)定模型的形式,可以得到較寬的非線性預(yù)測(cè),而且在進(jìn)行總體估計(jì)過程中,不依賴分布形式,具有偏差小、置信區(qū)間精確等特點(diǎn)[23],其表達(dá)式為:

    式中:Y為輸出值;Pi為固定設(shè)計(jì)點(diǎn);i為抽取樣本的點(diǎn)數(shù),i=1,…,m;m為總樣本的點(diǎn)數(shù);f(Di)是關(guān)于優(yōu)化設(shè)計(jì)變量的均值函數(shù);σ(Di)為關(guān)于優(yōu)化設(shè)計(jì)變量的方差函數(shù);εi為隨機(jī)誤差;Di是第i個(gè)設(shè)計(jì)變量。

    將拉丁超立方抽樣得到的樣本數(shù)據(jù)分別運(yùn)用Kriging模型和NPR模型構(gòu)建生產(chǎn)線機(jī)械臂的近似模型。由于輸入數(shù)據(jù)及輸出數(shù)據(jù)較多,本文僅對(duì)比Kriging模型和NPR模型中的兩組,其響應(yīng)面如圖12、13、14、15 所示。

    圖12 運(yùn)用Kriging模型D1、D9、第四階固有頻率的響應(yīng)面Fig.12 Response surface of Kriging model D1, D9, and fourth order natural frequencies

    圖13 運(yùn)用Kriging模型D11、D13、最大振幅的響應(yīng)面Fig.13 Response surface of maximum amplitude using Kriging models D11, D13

    圖14 運(yùn)用NPR模型D1、D9、第四階固有頻率的響應(yīng)面Fig.14 Response surface of NPR model D1, D9, and fourth order natural frequencies

    圖15 運(yùn)用NPR模型D11、D13、最大振幅的響應(yīng)面Fig.15 Response surface using NPR models D11, D13,and maximum amplitude

    從圖12~15 中可以看出,運(yùn)用Kriging模型與NPR模型構(gòu)建的生產(chǎn)線機(jī)械臂的近似模型具有一定區(qū)別,運(yùn)用NPR模型構(gòu)建的近似模型的響應(yīng)面相對(duì)比較光滑,運(yùn)用Kriging模型構(gòu)建的近似模型的響應(yīng)面邊緣相對(duì)比較明顯,可為生產(chǎn)線機(jī)械臂做出較好的預(yù)測(cè)。

    為進(jìn)一步檢驗(yàn)運(yùn)用Kriging模型與NPR模型構(gòu)建的生產(chǎn)線機(jī)械臂近似模型能否反映實(shí)際模型,需要對(duì)其進(jìn)行詳細(xì)評(píng)估,評(píng)估參數(shù)包括測(cè)定系數(shù)、均方根誤差、相對(duì)最大絕對(duì)誤差等。測(cè)定系數(shù)表示符號(hào)為R2,最佳值為1;均方根誤差表示符號(hào)為RMSE,最佳值為0;相對(duì)最大絕對(duì)誤差表示符號(hào)為RMAE,最佳值為0。測(cè)定系數(shù)、均方根誤差以及相對(duì)最大絕對(duì)誤差計(jì)算公式如下:

    式中:yi為第i個(gè)測(cè)試點(diǎn)的仿真值;為第i個(gè)測(cè)試點(diǎn)的近似值;為測(cè)試點(diǎn)仿真值的平均值;m為總樣本的點(diǎn)數(shù);k為測(cè)試的點(diǎn)數(shù);S為標(biāo)準(zhǔn)差。

    通過計(jì)算,Kriging模型與NPR模型構(gòu)建的生產(chǎn)線機(jī)械臂近似模型的測(cè)定系數(shù)、均方根誤差、相對(duì)最大絕對(duì)誤差如表7 所示。

    表7 Kriging模型與NPR模型對(duì)比Tab.7 Comparison between Kriging model and NPR model

    由表7 可知,Kriging模型擬合精度高于NPR模型,尤其是NPR模型Y方向最大振幅的均方根誤差偏離最佳值較遠(yuǎn)。因此,選擇Kriging模型構(gòu)建的近似模型。

    5 結(jié)構(gòu)優(yōu)化

    為方便分析,需要提前構(gòu)建生產(chǎn)線機(jī)械臂設(shè)計(jì)變量與響應(yīng)值之間的數(shù)學(xué)模型。將生產(chǎn)線機(jī)械臂設(shè)計(jì)變量作為變量,將生產(chǎn)線機(jī)械臂的第四階固有頻率P19以及質(zhì)量P20作為約束,將Y方向最大振幅P39最小作為目標(biāo)[24-29],構(gòu)建好的數(shù)學(xué)模型為:

    式中:Di為設(shè)計(jì)變量,Da為設(shè)計(jì)變量下限,Db為設(shè)計(jì)變量上限。

    運(yùn)用MOGA法進(jìn)行求解優(yōu)化,設(shè)置最初生成10 000個(gè)樣本,每次迭代生成3 200個(gè)樣本,并在最多20次迭代中找到3個(gè)候選樣本,通過求解篩選出3個(gè)候選點(diǎn),如表8所示。候選點(diǎn)1的優(yōu)先級(jí)高于候選點(diǎn)2與候選點(diǎn)3,因此選取候選點(diǎn)1進(jìn)行后續(xù)優(yōu)化[30-34]。

    表8 候選點(diǎn)數(shù)值表Tab.8 Table of candidate point values

    根據(jù)得到的候選點(diǎn)1,綜合考慮生產(chǎn)線機(jī)械臂設(shè)計(jì)變量與響應(yīng)值之間的關(guān)系,按照靈敏度正負(fù)關(guān)系,進(jìn)行人工調(diào)整,最終得出滿意的數(shù)值。生產(chǎn)線機(jī)械臂經(jīng)過調(diào)整后的尺寸如表9所示。

    表9 調(diào)整后的尺寸Tab.9 Adjusted dimensions

    按照人工調(diào)整得到的數(shù)據(jù),重新對(duì)生產(chǎn)線機(jī)械臂進(jìn)行建模,并對(duì)其進(jìn)行靜態(tài)和動(dòng)態(tài)性能分析,優(yōu)化后的數(shù)據(jù)如表10所示。

    表10 調(diào)整后的模型參數(shù)Tab.10 Adjusted model parameters

    將優(yōu)化后生產(chǎn)線機(jī)械臂的動(dòng)靜性能與優(yōu)化前生產(chǎn)線機(jī)械臂的動(dòng)靜性能進(jìn)行對(duì)比,可以發(fā)現(xiàn),優(yōu)化后生產(chǎn)線機(jī)械臂的最大等效應(yīng)力為2.124 6 MPa,低于其在受到約束情況時(shí)的最大等效應(yīng)力163.33 MPa,滿足靜強(qiáng)度要求;最大變形為0.023 547 mm,低于生產(chǎn)線機(jī)械臂可以承受的最大靜變形量0.194 mm,滿足靜剛度要求。此外,優(yōu)化后生產(chǎn)線機(jī)械臂的第四階模態(tài)固有頻率提升了5.15%,并且T型板在Y方向的振幅降低了73.20%,滿足生產(chǎn)線機(jī)械臂上下料過程中左右晃動(dòng)幅度小于等于2 mm的要求。同時(shí),優(yōu)化后生產(chǎn)線機(jī)械臂的質(zhì)量也沒有顯著增加。

    6 結(jié)論

    本文運(yùn)用建模軟件與有限元分析軟件的協(xié)同仿真技術(shù),研究了生產(chǎn)線機(jī)械臂的動(dòng)靜性能,通過Parameters Correlation模塊分析生產(chǎn)線機(jī)械臂設(shè)計(jì)變量對(duì)整體質(zhì)量、第四階模態(tài)固有頻率及Y方向最大振幅的靈敏度,篩選出對(duì)生產(chǎn)線機(jī)械臂影響較為明顯的設(shè)計(jì)變量。

    運(yùn)用拉丁超立方抽樣法得到篩選出的設(shè)計(jì)變量及對(duì)應(yīng)輸出的樣本,通過Kriging模型與NPR模型分別構(gòu)建生產(chǎn)線機(jī)械臂的近似模型,通過分析響應(yīng)面及擬合精度得出,Kriging模型構(gòu)建的近似模型優(yōu)于NPR模型構(gòu)建的近似模型。運(yùn)用Kriging模型構(gòu)建的近似模型對(duì)設(shè)計(jì)變量進(jìn)行了優(yōu)化,并對(duì)優(yōu)化后的數(shù)據(jù)進(jìn)行人工調(diào)整,再次建模并進(jìn)行靜態(tài)和動(dòng)態(tài)性能分析,結(jié)果表明,優(yōu)化后的生產(chǎn)線機(jī)械臂性能顯著提高,且其晃動(dòng)的幅度也滿足上下料精度要求。

    猜你喜歡
    型板固有頻率生產(chǎn)線
    方便小米粥億級(jí)生產(chǎn)線投入運(yùn)行
    現(xiàn)場(chǎng)測(cè)定大型水輪發(fā)電機(jī)組軸系的固有頻率
    對(duì)插式型膠擠出口型板設(shè)計(jì)
    16000t鍛造壓力機(jī)生產(chǎn)線將交付
    全鋼子午線輪胎胎面擠出強(qiáng)制式口型板的設(shè)計(jì)技術(shù)
    橡膠科技(2018年4期)2018-07-21 03:31:28
    基于RFID技術(shù)的口型板管理方案設(shè)計(jì)及應(yīng)用
    智能制造(2017年9期)2018-01-26 05:29:36
    兩種大砂箱多觸頭靜壓造型機(jī)
    生長(zhǎng)在生產(chǎn)線上
    Hazelett生產(chǎn)線熔煉工藝探討
    總溫總壓測(cè)頭模態(tài)振型變化規(guī)律研究
    永久网站在线| 中文天堂在线官网| 日韩欧美三级三区| 亚洲综合色惰| 亚洲自偷自拍三级| 2018国产大陆天天弄谢| 夫妻午夜视频| 久久综合国产亚洲精品| 性插视频无遮挡在线免费观看| av卡一久久| 国产精品国产三级国产专区5o| 成人一区二区视频在线观看| 国产免费福利视频在线观看| 精品久久久久久久久av| 一级黄片播放器| 啦啦啦韩国在线观看视频| 在线天堂最新版资源| 国产伦精品一区二区三区视频9| 日本免费a在线| 51国产日韩欧美| 观看免费一级毛片| 亚洲av日韩在线播放| 我的女老师完整版在线观看| 老司机影院成人| 日本一本二区三区精品| 国产一级毛片七仙女欲春2| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区四区激情视频| 亚洲av不卡在线观看| 欧美+日韩+精品| 亚洲精品,欧美精品| 搡老乐熟女国产| 少妇人妻精品综合一区二区| 国产成人一区二区在线| 亚洲精品久久久久久婷婷小说| 久久人人爽人人爽人人片va| 少妇人妻精品综合一区二区| 内射极品少妇av片p| 亚洲乱码一区二区免费版| 成人特级av手机在线观看| 午夜爱爱视频在线播放| 97超碰精品成人国产| 男女边摸边吃奶| 国产精品日韩av在线免费观看| 男人爽女人下面视频在线观看| 天堂网av新在线| 永久网站在线| 久久久精品欧美日韩精品| 最近的中文字幕免费完整| 日韩一区二区视频免费看| 我的老师免费观看完整版| 亚洲av福利一区| 精品欧美国产一区二区三| 国产免费又黄又爽又色| 久久久午夜欧美精品| 男女国产视频网站| 亚洲精品456在线播放app| 亚洲性久久影院| 亚洲精品,欧美精品| 偷拍熟女少妇极品色| 建设人人有责人人尽责人人享有的 | 国内精品一区二区在线观看| 日本wwww免费看| 国产成人a∨麻豆精品| 中文字幕免费在线视频6| 黄色配什么色好看| 热99在线观看视频| 国产精品嫩草影院av在线观看| 丝袜喷水一区| 在线观看av片永久免费下载| 搡老妇女老女人老熟妇| av黄色大香蕉| 一级毛片黄色毛片免费观看视频| 亚洲欧洲日产国产| 日韩 亚洲 欧美在线| 99久久精品一区二区三区| 日韩av不卡免费在线播放| 男女国产视频网站| 韩国高清视频一区二区三区| 成人亚洲精品一区在线观看 | 看十八女毛片水多多多| 亚洲18禁久久av| 成年人午夜在线观看视频 | 亚洲国产精品专区欧美| 在现免费观看毛片| 亚洲欧美成人精品一区二区| 亚洲成人精品中文字幕电影| 日韩欧美精品免费久久| 国产午夜精品一二区理论片| 国产黄频视频在线观看| 国产黄片视频在线免费观看| 精品亚洲乱码少妇综合久久| 久久久久久九九精品二区国产| 成人特级av手机在线观看| 人妻系列 视频| 精品少妇黑人巨大在线播放| 亚洲av.av天堂| 久久国内精品自在自线图片| kizo精华| 久久99热这里只有精品18| 一本久久精品| 18+在线观看网站| 街头女战士在线观看网站| 欧美日韩视频高清一区二区三区二| 久热久热在线精品观看| 亚洲电影在线观看av| 日韩av免费高清视频| 国产中年淑女户外野战色| 久久亚洲国产成人精品v| 精品人妻视频免费看| 久久6这里有精品| 日本三级黄在线观看| 99视频精品全部免费 在线| 亚州av有码| 亚洲欧美清纯卡通| 人妻系列 视频| 亚洲精品乱码久久久v下载方式| 国产伦精品一区二区三区四那| 白带黄色成豆腐渣| 国产成人精品久久久久久| 欧美bdsm另类| 国产高清国产精品国产三级 | 日韩欧美一区视频在线观看 | 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区三区| 最近中文字幕高清免费大全6| 啦啦啦啦在线视频资源| 欧美xxxx性猛交bbbb| 婷婷色av中文字幕| 国产欧美日韩精品一区二区| 欧美最新免费一区二区三区| 欧美日韩在线观看h| 一级a做视频免费观看| 秋霞在线观看毛片| 久久鲁丝午夜福利片| 亚洲av.av天堂| av女优亚洲男人天堂| 久久久久久久久中文| 国产一区二区三区av在线| av卡一久久| 精品久久久噜噜| 亚洲av国产av综合av卡| 丰满少妇做爰视频| 国产精品不卡视频一区二区| 日韩制服骚丝袜av| 51国产日韩欧美| 伊人久久精品亚洲午夜| 国产一区亚洲一区在线观看| 国产精品伦人一区二区| 99视频精品全部免费 在线| 日本猛色少妇xxxxx猛交久久| 久久精品久久久久久久性| 少妇高潮的动态图| 久久久精品94久久精品| 日日啪夜夜爽| 黄片无遮挡物在线观看| 亚洲国产最新在线播放| 亚洲性久久影院| www.色视频.com| 乱系列少妇在线播放| 不卡视频在线观看欧美| 中文精品一卡2卡3卡4更新| 免费黄色在线免费观看| av天堂中文字幕网| 天堂√8在线中文| 在线 av 中文字幕| 天堂网av新在线| 熟女电影av网| 男人舔女人下体高潮全视频| 一个人免费在线观看电影| 日韩一区二区三区影片| 日韩视频在线欧美| 人体艺术视频欧美日本| 免费播放大片免费观看视频在线观看| 国产黄色小视频在线观看| 非洲黑人性xxxx精品又粗又长| 国产成人精品一,二区| 成人亚洲欧美一区二区av| 又大又黄又爽视频免费| 日本免费a在线| 青春草国产在线视频| av专区在线播放| 国产熟女欧美一区二区| 又大又黄又爽视频免费| 午夜亚洲福利在线播放| 又爽又黄无遮挡网站| 在线天堂最新版资源| 久久久精品欧美日韩精品| 别揉我奶头 嗯啊视频| av网站免费在线观看视频 | 可以在线观看毛片的网站| 深夜a级毛片| 99久久精品国产国产毛片| 国产伦理片在线播放av一区| 久久人人爽人人爽人人片va| 久久久久久久久久久丰满| 日本-黄色视频高清免费观看| 激情 狠狠 欧美| 直男gayav资源| 男女啪啪激烈高潮av片| 欧美日韩国产mv在线观看视频 | 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 国产视频内射| 日本黄大片高清| 午夜亚洲福利在线播放| 亚洲av成人精品一区久久| 国产精品福利在线免费观看| 国产人妻一区二区三区在| 日日摸夜夜添夜夜添av毛片| 2021天堂中文幕一二区在线观| 听说在线观看完整版免费高清| 欧美丝袜亚洲另类| 亚洲第一区二区三区不卡| 国产精品女同一区二区软件| 亚洲丝袜综合中文字幕| 日韩欧美三级三区| 亚洲不卡免费看| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久av不卡| 国产成人一区二区在线| 久久久精品94久久精品| 精品久久久噜噜| eeuss影院久久| 亚洲欧美中文字幕日韩二区| 国产有黄有色有爽视频| 日日干狠狠操夜夜爽| 成年女人在线观看亚洲视频 | videossex国产| 中国美白少妇内射xxxbb| 97超视频在线观看视频| 国产v大片淫在线免费观看| 成人毛片60女人毛片免费| 国产高清有码在线观看视频| 久久久久免费精品人妻一区二区| 国产亚洲91精品色在线| 亚洲人成网站在线播| 91午夜精品亚洲一区二区三区| 亚洲婷婷狠狠爱综合网| www.av在线官网国产| 国产亚洲午夜精品一区二区久久 | 欧美日韩综合久久久久久| 日韩强制内射视频| 免费观看精品视频网站| 最近最新中文字幕免费大全7| 亚洲一区高清亚洲精品| av网站免费在线观看视频 | 久久鲁丝午夜福利片| 国国产精品蜜臀av免费| 成人美女网站在线观看视频| 国产成人91sexporn| 久久99热6这里只有精品| 大话2 男鬼变身卡| 一级毛片我不卡| 国产探花极品一区二区| 又粗又硬又长又爽又黄的视频| 成人亚洲精品av一区二区| 亚洲av国产av综合av卡| 国产人妻一区二区三区在| 三级男女做爰猛烈吃奶摸视频| 免费看美女性在线毛片视频| 国产午夜精品一二区理论片| 床上黄色一级片| 亚洲国产精品国产精品| 在线观看一区二区三区| 亚洲成人一二三区av| 男的添女的下面高潮视频| 午夜老司机福利剧场| 午夜免费观看性视频| 久久99热这里只有精品18| 一级二级三级毛片免费看| 一级av片app| 九草在线视频观看| 日日啪夜夜撸| 美女xxoo啪啪120秒动态图| 麻豆精品久久久久久蜜桃| 高清视频免费观看一区二区 | 99热全是精品| 欧美激情国产日韩精品一区| 久久久久精品久久久久真实原创| 一区二区三区乱码不卡18| 精品一区二区三卡| 日韩av不卡免费在线播放| 视频中文字幕在线观看| 亚洲精品国产av蜜桃| 国产高清不卡午夜福利| xxx大片免费视频| 国产精品美女特级片免费视频播放器| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 高清欧美精品videossex| 免费黄频网站在线观看国产| 男女啪啪激烈高潮av片| 青青草视频在线视频观看| 欧美成人午夜免费资源| 国产精品一区二区在线观看99 | 97热精品久久久久久| 国产精品麻豆人妻色哟哟久久 | 国产成人aa在线观看| 97超碰精品成人国产| 97人妻精品一区二区三区麻豆| 亚洲精品影视一区二区三区av| 亚洲人成网站高清观看| 婷婷色麻豆天堂久久| 成人午夜高清在线视频| 国产综合懂色| 欧美精品一区二区大全| 我要看日韩黄色一级片| 午夜福利成人在线免费观看| 欧美最新免费一区二区三区| 男女视频在线观看网站免费| 国产成人精品福利久久| 国产免费视频播放在线视频 | 男插女下体视频免费在线播放| 亚洲色图av天堂| 嫩草影院新地址| 国产精品不卡视频一区二区| 听说在线观看完整版免费高清| 国产激情偷乱视频一区二区| 国产成人福利小说| 秋霞伦理黄片| 色综合站精品国产| 国产精品人妻久久久久久| 一级a做视频免费观看| 天堂网av新在线| 在线免费观看的www视频| 99久久人妻综合| 久久精品久久久久久噜噜老黄| 国产单亲对白刺激| 人人妻人人澡欧美一区二区| 天堂俺去俺来也www色官网 | 亚洲欧美一区二区三区国产| 网址你懂的国产日韩在线| a级毛片免费高清观看在线播放| 亚洲欧美精品自产自拍| 国产一区二区亚洲精品在线观看| 久久精品国产亚洲网站| 国产精品.久久久| 成人国产麻豆网| 男插女下体视频免费在线播放| 国产伦精品一区二区三区四那| 免费看a级黄色片| 欧美日韩国产mv在线观看视频 | 午夜日本视频在线| 97热精品久久久久久| 久久精品国产亚洲av涩爱| 亚洲欧美成人综合另类久久久| 成人av在线播放网站| 欧美不卡视频在线免费观看| 亚洲熟妇中文字幕五十中出| 国产精品久久久久久久久免| 欧美三级亚洲精品| 国产av国产精品国产| 丝袜喷水一区| 2021天堂中文幕一二区在线观| 男女视频在线观看网站免费| 99热这里只有是精品50| 两个人的视频大全免费| 99九九线精品视频在线观看视频| 又爽又黄a免费视频| 久久久久久九九精品二区国产| 亚洲欧美日韩卡通动漫| 黄色欧美视频在线观看| 亚洲国产精品sss在线观看| 国内精品宾馆在线| 免费看光身美女| 日韩av免费高清视频| 中文字幕制服av| 欧美精品一区二区大全| 欧美zozozo另类| 国产激情偷乱视频一区二区| 日韩精品有码人妻一区| 免费黄色在线免费观看| 亚洲成人av在线免费| 国产精品久久久久久av不卡| 精品国产露脸久久av麻豆 | 内射极品少妇av片p| 日韩成人av中文字幕在线观看| 搞女人的毛片| 少妇猛男粗大的猛烈进出视频 | 色综合亚洲欧美另类图片| 国产乱人偷精品视频| 久久久成人免费电影| 亚洲图色成人| 99久久中文字幕三级久久日本| 中文字幕人妻熟人妻熟丝袜美| 国产成人免费观看mmmm| 免费观看精品视频网站| 淫秽高清视频在线观看| 麻豆成人av视频| 视频中文字幕在线观看| 网址你懂的国产日韩在线| 久久午夜福利片| 在线天堂最新版资源| 亚洲av免费高清在线观看| 两个人的视频大全免费| 免费在线观看成人毛片| 五月伊人婷婷丁香| 精品人妻一区二区三区麻豆| av又黄又爽大尺度在线免费看| 青春草视频在线免费观看| 亚洲欧洲日产国产| 国产伦在线观看视频一区| 青青草视频在线视频观看| 免费看光身美女| 91久久精品国产一区二区成人| 国产精品一区二区在线观看99 | 看十八女毛片水多多多| 午夜免费激情av| 成年女人看的毛片在线观看| 国产黄色视频一区二区在线观看| 日韩欧美精品v在线| 欧美日韩亚洲高清精品| 国产成人免费观看mmmm| av专区在线播放| or卡值多少钱| 岛国毛片在线播放| 婷婷色综合www| 边亲边吃奶的免费视频| 韩国高清视频一区二区三区| 免费播放大片免费观看视频在线观看| 观看免费一级毛片| 男女边吃奶边做爰视频| 成人一区二区视频在线观看| 亚洲一级一片aⅴ在线观看| 亚洲精品456在线播放app| 国产精品综合久久久久久久免费| 成人午夜高清在线视频| 97超碰精品成人国产| 欧美精品国产亚洲| 黑人高潮一二区| 色哟哟·www| 国产乱来视频区| 国产高潮美女av| 男插女下体视频免费在线播放| 伦理电影大哥的女人| 18禁在线播放成人免费| 免费人成在线观看视频色| 高清日韩中文字幕在线| 晚上一个人看的免费电影| 亚洲精品久久久久久婷婷小说| 在线 av 中文字幕| 国产高清不卡午夜福利| 久久精品久久精品一区二区三区| 精品一区二区三卡| 国产 亚洲一区二区三区 | 日韩亚洲欧美综合| 波野结衣二区三区在线| 日本爱情动作片www.在线观看| 又爽又黄无遮挡网站| 欧美xxxx性猛交bbbb| 亚洲国产高清在线一区二区三| 美女内射精品一级片tv| 久久久久久久久久成人| 国产精品爽爽va在线观看网站| 国产探花在线观看一区二区| 一级毛片 在线播放| 熟妇人妻不卡中文字幕| 美女大奶头视频| 麻豆av噜噜一区二区三区| 汤姆久久久久久久影院中文字幕 | 国产在线一区二区三区精| 少妇丰满av| 亚洲精品456在线播放app| 2021天堂中文幕一二区在线观| 亚洲人与动物交配视频| 亚洲成人精品中文字幕电影| 中文字幕人妻熟人妻熟丝袜美| 国产黄a三级三级三级人| av在线亚洲专区| 大片免费播放器 马上看| 国产免费一级a男人的天堂| 99热这里只有是精品在线观看| 国内揄拍国产精品人妻在线| 大又大粗又爽又黄少妇毛片口| 偷拍熟女少妇极品色| 免费av不卡在线播放| 久久久久精品性色| 美女主播在线视频| 午夜精品一区二区三区免费看| 噜噜噜噜噜久久久久久91| 伦理电影大哥的女人| 一本一本综合久久| 亚洲精品,欧美精品| 内地一区二区视频在线| 有码 亚洲区| 日本熟妇午夜| 午夜精品在线福利| 国产午夜精品久久久久久一区二区三区| 精品久久久久久成人av| 床上黄色一级片| 久久精品国产亚洲av涩爱| 插阴视频在线观看视频| 日韩精品青青久久久久久| 亚洲国产成人一精品久久久| 国产成人精品婷婷| 欧美激情国产日韩精品一区| 美女cb高潮喷水在线观看| 国产精品国产三级专区第一集| 成人亚洲精品av一区二区| av又黄又爽大尺度在线免费看| 日日摸夜夜添夜夜添av毛片| 一级毛片我不卡| 亚洲av中文av极速乱| 久久久久久久国产电影| 一级a做视频免费观看| 久久鲁丝午夜福利片| 日本欧美国产在线视频| 国产一区二区三区av在线| 观看免费一级毛片| 免费看美女性在线毛片视频| 肉色欧美久久久久久久蜜桃 | 国产精品99久久久久久久久| 婷婷色综合www| 自拍偷自拍亚洲精品老妇| 天堂影院成人在线观看| 国产精品精品国产色婷婷| 色网站视频免费| 韩国av在线不卡| 久久久成人免费电影| 国产午夜精品一二区理论片| 久久午夜福利片| 日韩精品青青久久久久久| 久久精品国产亚洲av天美| 久久久久久久午夜电影| 国产精品1区2区在线观看.| 51国产日韩欧美| 亚洲怡红院男人天堂| 日韩一区二区视频免费看| 日韩av在线大香蕉| 精品人妻偷拍中文字幕| 日本三级黄在线观看| 男女国产视频网站| 欧美日韩综合久久久久久| 韩国av在线不卡| 国产精品一及| 国产伦在线观看视频一区| 美女cb高潮喷水在线观看| 欧美日韩在线观看h| 亚洲在线观看片| 久久精品久久精品一区二区三区| 亚洲精品aⅴ在线观看| 狠狠精品人妻久久久久久综合| 久久久国产一区二区| 黄色配什么色好看| 亚洲欧美一区二区三区国产| 菩萨蛮人人尽说江南好唐韦庄| 免费看不卡的av| 国语对白做爰xxxⅹ性视频网站| 人人妻人人澡欧美一区二区| 久久久久久久亚洲中文字幕| 搡老妇女老女人老熟妇| 边亲边吃奶的免费视频| 欧美成人午夜免费资源| 嫩草影院精品99| 日本免费在线观看一区| 91精品国产九色| 精品久久久久久久久av| 免费av不卡在线播放| a级毛色黄片| 亚洲自偷自拍三级| 黄色一级大片看看| 国产麻豆成人av免费视频| 99久久精品热视频| 日韩成人伦理影院| 99久久精品国产国产毛片| 日本色播在线视频| 国产伦精品一区二区三区四那| 精品久久久久久久末码| 2021天堂中文幕一二区在线观| 一级爰片在线观看| 成人午夜高清在线视频| 日韩精品青青久久久久久| 国产大屁股一区二区在线视频| 日韩成人av中文字幕在线观看| 国产女主播在线喷水免费视频网站 | 大话2 男鬼变身卡| 亚洲人成网站在线播| 狂野欧美激情性xxxx在线观看| 高清午夜精品一区二区三区| 国产v大片淫在线免费观看| 有码 亚洲区| 成人亚洲欧美一区二区av| 干丝袜人妻中文字幕| 亚洲内射少妇av| 亚洲成人av在线免费| 亚洲精品中文字幕在线视频 | 老师上课跳d突然被开到最大视频| 一级毛片久久久久久久久女| 国产女主播在线喷水免费视频网站 | freevideosex欧美| 99热这里只有是精品50| 色播亚洲综合网| 国产黄色免费在线视频| 国产一区有黄有色的免费视频 | 国产伦一二天堂av在线观看| 亚洲国产色片| 99久久精品热视频| 欧美激情国产日韩精品一区| 精品人妻熟女av久视频| 亚洲av男天堂| 久久鲁丝午夜福利片| 一级黄片播放器| xxx大片免费视频| 亚洲经典国产精华液单| 免费看av在线观看网站| 久久久久久国产a免费观看| 一级毛片黄色毛片免费观看视频| 久久久久久久大尺度免费视频| 国产精品爽爽va在线观看网站| 亚洲欧美精品自产自拍| 日本三级黄在线观看| 青青草视频在线视频观看| 午夜福利高清视频| 十八禁国产超污无遮挡网站|