999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

2-PRC-PR(C)并聯機構減振平臺及其動力學特性

2017-05-25 03:46:19劉乃軍牛軍川
中南大學學報(自然科學版) 2017年4期
關鍵詞:模態影響系統

劉乃軍,牛軍川

?

2-PRC-PR(C)并聯機構減振平臺及其動力學特性

劉乃軍,牛軍川

(山東大學機械工程學院,高效潔凈機械制造教育部重點實驗室,山東濟南,250061)

為滿足工程中機械設備對多維振動和沖擊的隔離要求,設計一種新型2-PRC-PR(C)并聯機構減振平臺。該平臺是在3-PRC并聯機構的基礎上,通過把一支鏈和上平臺連接處的圓柱副改為一復合運動副,使該機構實現三自由度3-PRC并聯機構與四自由度2-PRC-PR^C并聯機構之間的切換,并在主動副處輔以彈簧阻尼,使該機構能夠實現三維隔振與四維隔振。對該機構進行運動學和動力學分析,在此基礎上將系統參量分類,并對參數較多類別采用正交試驗法,分析各參數對系統固有頻率的影響規律和貢獻量。研究結果表明:多維減振平臺中有關參量對系統固有頻率的影響具有不同于單自由度減振裝置的特性。

并聯機構;減振平臺;固有頻率;正交試驗法

精密儀器儀表、車載、船載、航空航天等機械設備在工作時都會受到大量的多維振動和沖擊,成為影響設備性能和使用壽命的主要因素[1?3]。傳統的隔振技術已不能滿足當前生產實際的需求,主要表現在2個方面。首先,傳統的隔振方法主要是采用彈簧、柔性橡膠等彈性材料或多層常規單自由度彈性阻尼裝置來進行隔振。由于彈性材料的性能不穩定,且對多維減振效果不明顯,而多層常規單自由度彈性阻尼裝置結構十分復雜,且僅能滿足少維隔振的需求;其次,隨著科技的發展,眾多機械設備的復雜化程度越來越高,同一機械設備的工況及工作條件存在多樣性,其所受振動和沖擊的維數會隨之改變,這對隔振裝置提出了新的要求。為解決多維隔振問題,近年來眾多學者將具有剛度大、結構穩定、承載能力強、運動慣量小等優點的并聯機構應用于振動控制領域。馬履中等[4?6]提出了可采用并聯機構作為多維減振平臺的基體,并列舉了部分適于做多維減振平臺的多自由度機型;YIN等[7]以3-RRC并聯機構作為減振主體,對多維減振系統的模態進行了分析;LI等[8]基于雙層并聯機構得到了一種隔振頻帶較寬的減振平臺;GUO等[9]基于4-CPS/RPS并聯機構,得到了一種船載主動隔振系統。近年來,諸多學者將具有結構簡單、動態特性好等特點的3-PRC并聯機構應用于多維減振裝置,在一定程度上解決了多維振動問題。馬履中等[10]論述了3-PRC并聯機構應用于空間三平動方向三維減振裝置的可行性;楊峰等[11]討論了3-PRC并聯機構的可調頻減振特性;李勇[12]則對3-PRC等少自由度并聯機構的減振進行了探討。然而,上述文獻所研究的減振裝置的減振維數都不可變更,往往不能滿足實際生產中振動和沖擊的維數隨時變化場合的隔振要求。在解決機械設備受到可變維數振動和沖擊的問題方面,牛軍川等[13]研究了一種3-CRR及其變胞減振機構,為變維數的振動隔離提供了一種思路。為很好地對受到多維且可變維數振動和沖擊的機械設備進行隔振,尤其是對如行駛在崎嶇路面上的車輛的空間三平動和一轉動的設備進行隔振,本文在3-PRC機構的基礎上給出了一種2-PRC-PR(C)并聯機構。通過調節機構中復合運動副(C),該機構可實現三自由度3-PRC并聯機構與四自由度2-PRC-PR^C并聯機構之間的切換。將2-PRC-PR(C)并聯機構作為減振平臺的主體,再在相應的主動副處輔以彈簧阻尼,能夠同時滿足三維隔振與四維隔振的要求。在此基礎上研究其運動學、動力學特性?;诓⒙摍C構設計的減振平臺與常規單自由度彈性減振裝置的區別在于影響系統固有頻率的參數眾多,分析困難。本文作者將系統參數分類,并將所含參數較多類別采用正交試驗法,分析減振平臺各參數對系統固有頻率的影響特性及貢獻量。

1 并聯機構多維減振平臺模型

2-PRC-PR(C)并聯機構是在3-PRC并聯機構的基礎上經過修改而來,即把3-PRC機構3點處的1個C副(圓柱副)修改為一復合運動副(C),如圖1所示。圖2所示為該裝置的等效簡化構型和復合運動副(C)的結構圖。可見:該機構由動平臺(上平臺)、靜平臺(下平臺)及(123),(123)和(123)3條支鏈組成。其中A和C這2條支鏈結構相同,均由1個移動副(P副)、1個與P副垂直的轉動副(R副)和1個與R副平行的圓柱副(C副)串聯而成,B支路由1個移動副(P副)、1個與P副垂直的轉動副(R副)和1個復合運動副(C)串聯而成。該復合運動副實際上是把原來的C副變成了1個移動方向不變的P副,然后在水平方向又固結了1個與之正交的十字軸,如圖2(b)所示,該十字軸3分別又與滑塊2和連桿4形成了2個轉動副。并聯機構通常存在多種構型。為使機構構型緊湊,本文中取2-PRC-PR(C)并聯機構構型如圖2(a)所示,即3個連桿3,3和3都向上平臺內側方向傾斜。

容易理解,當把復合運動副(C)中的十字軸3和滑塊2的R副固定時,則2-PRC-PR(C)并聯機構退化為3-PRC并聯機構;當復合運動副(C)中的十字軸3和滑塊2的R副存在相對轉動時,2-PRC-PR(C)并聯機構轉化為2-PRC-PR^C并聯機構,且對2-PRC-PR^C并聯機構的分析可等效轉化到3-PRC并聯機構。因此,本文僅對2-PRC-PR^C并聯機構進行分析。

圖1 2-PRC-PR(C)聯機構原理圖

(a) 構型;(b) 復合運動副

2 減振平臺運動學分析

在并聯機構上建立笛卡兒坐標系如圖1所示,根據基于約束螺旋理論的自由度分析方法[14],易得2-PRC-PR(C)并聯機構轉化所得到的2-PRC-PR^C并聯機構,其自由度為4;2-PRC-PR(C)并聯機構轉化為3-PRC并聯機構,自由度為3。

2.1 正、逆運動學分析

根據幾何條件,可得如下位姿約束方程。

(2)

(3)

若給定該機構上平臺點位置的4個運動學輸出x,yz,可得該機構的運動學逆解:

(5)

(6)

若給定該機構的4個運動學輸入1,1,1和1,由式~得該機構的運動學正解。

2.2 速度分析

為分析機構的動力學特性,這里首先對速度進行分析。由式~可得:

(9)

(10)

將式~對時間求導并整理可得

(12)

(14)

(15)

(17)

(18)

若矩陣非奇異,由式可得

3 減振機構動力學模型及其固有頻率

為方便構建基于2-PRC-PR^C并聯機構的隔振系統動力學模型,本文采用改進后的拉格朗日方程建立系統的動力學方程:

式中,=1,2,…,;s為廣義坐標;Γ為沿廣義坐標s方向作用的廣義力;為機構動能;為機構勢能;為系統的阻尼能量耗散系數。

這里只考慮上平臺質量,忽略各連桿質量和系統阻尼,可計算出機構的動能、勢能,將其代入式,可得2-PRC-PR^C機構的動力學方程為

(22)

mI分別為上平臺質量和繞軸的轉動慣量;k為系統的彈簧剛度。

系統自由振動時,令=0,由式可得

可以得到系統的特征方程為

(25)

4 系統各參數對固有頻率的影響特性

對所設計并聯機構減振平臺初始參數取為:1=1.050 m,1=1=0.146 m,2=2=2=0.6 m,3=3=3=0.5 m,4=0.2 mk(=1,2,3)=400 kN/m,m=50 kg,44 000 N?m/rad,I=0.627 kg/m2。得到該減振平臺處于2-PRC-PR^C和3-PRC這2種不同機構下的固有頻率如表1所示。進行仿真分析得減振平臺處于2-PRC-PR^C機構形位下的四階模態振型如圖3所示。

由表1及圖3可得:該減振平臺的前3階固有頻率分別沿,和軸移動方向,減振平臺處于2-PRC-PR^C機構狀態時存在的第4階固有頻率沿軸的轉動方向;減振平臺處于2種不同機構形式下的第1階、第3階固有頻率完全相同,第2階固有頻率略微存在差異,差值小于0.2 Hz??梢姡簻p振平臺處于3-PRC機構這一形位下的3階模態信息,可通過略去減振平臺處于2-PRC-PR^C機構形位下的第4階模態等效得到。

系統固有頻率與系統彈簧剛度、慣性參數、減振平臺各幾何尺寸等有關。參數眾多,規律復雜,因此,理清各參數對系統固有頻率的影響規律對減振平臺設計尤為重要。鑒于減振平臺3類參數(彈簧剛度參數、慣性參數、減振平臺各幾何尺寸參數)之間的獨立性,下面依次分類討論其對系統固有頻率的影響特性。

表1 2-PRC-PR^C和3-PRC固有頻率

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

4.1 彈簧剛度對固有頻率的影響

采用控制變量法,系統彈簧剛度1,2,3和4四因素中某單一因素在(101~1015N/m)范圍內變化,而其余3因素保持為初始參數,分別得到4因素1,2,3,4對系統4個固有頻率的影響特性,如圖4~7 所示。

由圖4~7可得如下規律:

1) 彈簧剛度在較小取值范圍(101~106N/m)內變化,減振系統4個固有頻率隨彈簧剛度1,2,3和4的增大而增大,該規律類似于單自由度減振系統。

2) 彈簧剛度在較大范圍(106~1012N/m)內變化,系統固有頻率4隨彈簧剛度的增大而增大,而系統固有頻率1,2,3都相應出現一段“停滯期”,即固有頻率隨彈簧剛度的增大保持不變。彈簧剛度在極大范圍(1012~1015N/m)內變化,系統固有頻率4隨彈簧剛度1,2,3和4的增大而增大,而1,2和3都表現出一些不規律現象。可見由于減振系統中各彈簧剛度對系統剛度矩陣的強耦合作用,使某一彈簧剛度在大范圍內變化時導致了系統固有頻率產生了不同于單自由度減振系統的特性。

為綜合分析系統彈簧剛度1,2,3和4分別對4個固有頻率1,2,3和4的影響程度,本文采用正交試驗法[15?16]進行分析。以系統4個彈簧剛度作為分析因素,綜合考慮每個變量取為3個水平,因素水平取值如表2所示,構造正交試驗表9(34)如表3所示。

將上述各數據組合代入式,對所得計算結果采用直觀分析方法,可得第因素極差:

式中:=1,2,…,;X為第因素第水平下的第個數據結果;和分別為正交試驗中的因素數和水平數;為各水平下的數據個數。

由式可得系統固有頻率1,2,3和4分別在4因素1,2,3和4下的相應極差,如圖8所示。

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

表2 系統彈簧剛度因素水平表

表3 L9(34)正交實驗表

圖8 彈簧剛度與系統固有頻率的極差關系

由極差容易判斷出系統固有頻率f(=1,2,3,4)受彈簧剛度1,2,3和4因素的影響規律,也易判斷某彈簧剛度k(=1,2,3,4)分別對系統固有頻率1,2,3和4影響的關系。由圖8可得:系統固有頻率4,3,2和1受彈簧剛度影響程度依次遞減;彈簧剛度4,1,2和3對固有頻率4影響程度依次遞減,彈 簧剛度4,1,3和2對固有頻率3影響程度依次遞減,彈簧剛度4對固有頻率2影響較大,1,2和3對固有頻率2影響較小;彈簧剛度4對固有頻率1影響較大,1,2和3對固有頻率1影響較小。

4.2 慣性參數對固有頻率的影響

上平臺質量m及其繞軸的轉動慣量I對系統固有頻率的影響特性分別如圖9及圖10所示。從圖9和圖10可以看出:減振平臺固有頻率1,2,3和4都隨上平臺的質量及繞軸的轉動慣量增加而減小,此性質類似于單自由度減振系統。

4.3 幾何尺寸對固有頻率的影響

如圖1所示,本文所設計減振平臺的幾何尺寸包括減振平臺主動輸入1,1,1,1以及減振平臺各桿件長度2,2,2,3,3,3,4和1。由式可得出,1,2和2對系統固有頻率無影響,故只需討論幾何尺寸1,1,2,3,3,3,4和1對系統4個固有頻率的影響關系。其中1和1對系統固有頻率的影響特性如圖11和圖12所示。

采用正交試驗法,以上述8個參量作為分析因素,綜合考慮每個變量取2個水平,因素水平取值如表4所示,構造正交試驗表12(28)如表5所示。

將各數據組合代入式,并將所得各結果代入式,可得系統固有頻率1,2,3,4分別在8因素1,1,2,3,3,3,4和1下的相應極差,如圖13所示。

由圖13可以看出:1) 系統固有頻率4受幾何尺寸變化影響較大,且4,1,3,3,3,1,2和1對其影響效果依次遞減;2)1,1,3,1,3,4,2和3亦對固有頻率3產生影響,作用依次遞減;3) 固有頻率1和2受幾何尺寸的影響較小。

表4 幾何尺寸因素水平表

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

(a) 第1階模態振型;(b) 第2階模態振型;(c) 第3階模態振型;(d) 第4階模態振型

表5 L12(28)正交實驗表

圖13 幾何尺寸與系統固有頻率的極差關系

5 結論

1) 基于3-PRC并聯機構,引入了1個復合運動副(C),設計了一種2-PRC-PR(C)并聯機構,該機構可通過改變復合運動副(C)的運動規律,轉化為三自由度3-PRC并聯機構或四自由度2-PRC-PR^C并聯機構?;谠摬⒙摍C構,在4個主動副處安裝彈簧阻尼,得到了基于2-PRC-PR(C)的多維減振平臺,可以滿足三維隔振與四維隔振工況要求。

2) 對該減振平臺進行了運動學、動力學分析。得到減振平臺處于3-PRC機構這一形位下的3階模態信息,可通過略去減振平臺處于2-PRC-PR^C機構形位下的第4階模態等效得到。

3) 影響減振平臺固有頻率的參數較多。鑒于部分參數之間的獨立性,本文將系統參數分為彈簧剛度,慣性參數以及機構幾何尺寸3類。對參數較多類別采用正交試驗方法,分析了各系統參數對系統固有頻率的影響特性。得到了各參數對各方向固有頻率的影響規律和貢獻量以及有關參量對系統固有頻率的影響不同于單自由度減振裝置的特性。本文為分析各參量對多維減振平臺固有頻率影響特性提供了一種新的思路,也為后期工作中設計具有特定固有頻率組合的多維減振平臺提供了依據。

[1] 牛軍川, 宋孔杰. 多激勵多支承全柔性隔振系統的傳遞特性分析[J]. 機械工程學報, 2011, 47(7): 59?64. NIU Junchuan, SONG Kongjie. Transmission characteristics of fully flexible isolation systems subjected to multi-excitations and supported by multi-mounts[J]. Journal of Mechanical Engineering, 2011, 47(7): 59?64.

[2] RAM Y M, MOTTERSHEAD J E. Multiple-input active vibration control by partial pole placement using the method of receptances[J]. Mechanical Systems and Signal Processing, 2013, 40(2): 727?735.

[3] ANTONIO Z, GIOVANNI F B, MARCO A B. Active vibration control of a composite sandwich plate[J]. Composite Structures, 2015, 128: 100?114.

[4] 馬履中, 楊啟志, 尹小琴, 等. 仿橡膠多自由度彈性阻尼減振裝置研究[J]. 機械工程學報, 2004, 40(1): 21?24. MA Luzhong, YANG Qizhi, YI Xiaoqing, et al. Analysis of combined elastic damping devices on multiple DOF parallel mechanism[J]. Journal of Mechanical Engineering, 2004, 40(1): 21?24.

[5] 許子紅, 馬履中, 劉慶波, 等. 并聯機構在多維減振平臺中的應用研究[J]. 機械設計與制造, 2007(11): 124?126. XU Zihong, MA Lüzhong, LIU Qingbo, et al. Application study of parallel mechanism for vibration reduction[J]. Machinery Design and Manufacture, 2007(11): 124?126.

[6] 朱偉, 馬履中, 吳偉光, 等. 基于三平移并聯機構的三維減振平臺建模與仿真[J]. 農業機械學報, 2008, 39(1): 142?146. ZHU Wei, MA Lüzhong, WU Weiguang, et al. Modeling and simulation analysis on multi-dimensional damping platform base on three-translation parallel Mechanism[J]. Transactions of the Chinese Society for Agricultural Machinery, 2008, 39(1): 142?146.

[7] YIN Xiaoqin, YANG Qizhi, XIE Jun, et al.Vibration analysis of multi-DIM vibration absorber based on parallel mechanism[C]//Second International Conference on Mechanic Automation and Control Engineering (MACE). Baotou, China, 2011: 7881?7884.

[8] LI Bing, ZHAO Wei, DENG Zongquan. Modeling and analysis of a multi-dimensional vibration isolator based on the parallel mechanism[J]. Journal of Manufacturing Systems, 2012, 31(1): 50?58.

[9] GUO Jun, LI Guotao, LI Bing, et al. A ship active vibration isolation system based on a novel 5-DOF parallel mechanism[C]//Proceeding of the IEEE International Conference on Information and Automation. Hailar, China, 2014: 800?805.

[10] 馬履中, 尹小琴, 楊啟志, 等. 多維減振平臺主體機構的分析研究[J]. 江蘇大學學報(自然科學版), 2004, 25(4): 281?284. MA Lüzhong, YIN Xiaoqin, YANG Qizhi, et al. Study on main mechanism of multi-DIM vibration absorber[J]. Journal of Jiangsu University (Natural Science Edition), 2004, 25(4): 281?284.

[11] 楊峰, 牛軍川, 李勇. 可調頻并聯減振機構的動力學分析與應用研究[J]. 現代振動與噪聲技術, 2012, 10: 112?115. YANG Feng, NIU Junchuan, LI Yong. Application and dynamics of tunable parallel mechanism for vibration reduction[J]. Modern Vibration and Noise Techniques, 2012, 10: 112?115.

[12] 李勇. 基于少自由度的并聯機構減振平臺設計[D]. 濟南: 山東大學機械工程學院, 2014: 11?23. LI Yong. Study on multi-dimensional vibration isolation platform based on limited-dof parallel mechanism[D]. Jinan: Shandong University. School of Mechanical Engineering, 2014: 11?23.

[13] 牛軍川, 楊峰, 李勇. 基于變胞并聯機構的多維隔振研究[J]. 振動與沖擊, 2014, 33(7): 206?209. NIU Junchuan, YANG Feng, LI Yong. Study on multi-dimensional vibration isolation based on metamorphic parallel mechanism[J]. Journal of Vibration and Shock, 2014, 33(7): 206?209.

[14] LI Qinchuan, HUANG Zhen. Mobility analysis of lower-mobility parallel manipulators based on screw theory[C]// IEEE International Conference on Robotics and Automation. Taipei, 2003: 1179?1184.

[15] FRANEK L, JIANG Xiaoyi. Orthogonal design of experiments for parameter learning in image segmentation[J]. Signal Processing, 2013, 93(6): 1694?1704.

[16] JI Lijun, SI Yunfeng, LIU Hongfei, et al. Application of orthogonal experimental design in synthesis of mesoporous bioactive glass[J]. Microporous and Mesoporous Materials, 2014, 184(1): 122?126.

(編輯 楊幼平)

2-PRC-PR(C) parallel mechanism platform for vibration isolation and its dynamics

LIU Naijun, NIU Junchuan

(Shandong University Key Laboratory of High-efficiency andClean Mechanical Manufacture,Ministry of Education, School of Mechanical Engineering, Shandong University, Jinan 250061, China)

In order to meet the requirement to reduce the multi-dimensional vibration and shock from equipment and machines, a vibration isolation platform was proposed based on a 2-PRC-PR(C) parallel mechanism. This mechanism came from a 3-PRC parallel mechanism by replacing one of the cylindrical pairs at the moving platform with a combined pair. Thanked to the combined pair, the proposed platform could be switched between 3-PRC parallel mechanism with 3 degrees of freedom and 2-PRC-PR^C parallel mechanism with 4 degrees of freedom. 3 or 4 dimensional vibration isolation could be realized when spring damping were put on the active motion pairs. The kinematics and dynamics of the system were studied, and then the influence principle and contribution of the platform’s parameters on its natural frequency were analyzed by orthogonal experiment method. The results show that the relevant parameters effects on the natural frequency of the system are different from the single degree of freedom vibration reduction device.

parallel mechanism; vibration isolation platform; natural frequency; orthogonal experiment method

TH112

A

1672?7207(2017)04?0925?11

10.11817/j.issn.1672?7207.2017.04.011

2015?05?06;

2015?06?16

國家自然科學基金資助項目(51275275,51675306);山東省優秀中青年科學家科研獎勵基金資助項目(BS2010ZZ006) (Projects(51275275, 51675306) supported by the National Natural Science Foundation of China; Project(BS2010ZZ006) supported by Research Award Fund for Outstanding Young Scientists of Shandong Province)

牛軍川,博士,教授,從事機械系統動力學、振動與噪聲控制研究;E-mail:niujc@sdu.edu.cn

猜你喜歡
模態影響系統
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
是什么影響了滑動摩擦力的大小
哪些顧慮影響擔當?
當代陜西(2021年2期)2021-03-29 07:41:24
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
擴鏈劑聯用對PETG擴鏈反應與流變性能的影響
中國塑料(2016年3期)2016-06-15 20:30:00
國內多模態教學研究回顧與展望
基于HHT和Prony算法的電力系統低頻振蕩模態識別
由單個模態構造對稱簡支梁的抗彎剛度
計算物理(2014年2期)2014-03-11 17:01:39
主站蜘蛛池模板: 亚洲成人免费看| 国产成人精品免费视频大全五级 | 亚洲香蕉久久| 国产网站免费| 日本成人一区| 波多野结衣一二三| 精品视频91| 欧美视频在线播放观看免费福利资源| 欧美日韩免费观看| 午夜福利视频一区| 国产99热| 亚洲欧美日韩另类在线一| 国产亚洲男人的天堂在线观看| 97一区二区在线播放| 婷婷六月在线| 久久免费视频6| 91小视频在线观看| 看国产毛片| 午夜天堂视频| 一级爆乳无码av| 伊人AV天堂| 女人天堂av免费| 亚洲伦理一区二区| 成人福利在线免费观看| 97久久精品人人做人人爽| 97国产在线观看| 亚洲天堂.com| 夜夜操国产| 久久精品国产在热久久2019| 中文字幕 91| 国产一级毛片yw| 国产呦精品一区二区三区下载| 亚洲成人免费看| 精品国产成人三级在线观看| 国产老女人精品免费视频| 亚洲人在线| 免费AV在线播放观看18禁强制| 最新痴汉在线无码AV| 色久综合在线| 国产三级国产精品国产普男人 | 色国产视频| 一本一本大道香蕉久在线播放| 91精品免费久久久| 青青草原国产免费av观看| 情侣午夜国产在线一区无码| 国产欧美成人不卡视频| 久久精品无码一区二区国产区| 久久久久国产精品免费免费不卡| 5388国产亚洲欧美在线观看| 在线看免费无码av天堂的| 91美女视频在线观看| 免费国产高清精品一区在线| 日本道中文字幕久久一区| 人妻精品久久无码区| 啦啦啦网站在线观看a毛片| 99精品福利视频| 一级毛片网| 亚洲人网站| 不卡无码h在线观看| 国产高清精品在线91| 精品国产乱码久久久久久一区二区| 精品久久久久久久久久久| 40岁成熟女人牲交片免费| 成人a免费α片在线视频网站| 一级一级一片免费| 久久久久人妻一区精品| a亚洲视频| 婷婷综合在线观看丁香| 国产大全韩国亚洲一区二区三区| 日韩无码黄色网站| 免费a级毛片视频| 毛片网站免费在线观看| 亚洲V日韩V无码一区二区| 一本色道久久88综合日韩精品| 夜夜操天天摸| 无码专区国产精品第一页| 日本a∨在线观看| 一级全黄毛片| 亚洲一区无码在线| 国产精品大尺度尺度视频| 欧美成人免费午夜全| 亚洲国产高清精品线久久|