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

新型雙臂菱形壓電柔性機構(gòu)理論設(shè)計與建模

2017-11-15 02:32:49曹軍義凌明祥曾明華
中國機械工程 2017年21期
關(guān)鍵詞:有限元模型

蔣 州 曹軍義 凌明祥,2 曾明華 林 京

1.西安交通大學機械工程學院,西安,7100492.中國工程物理研究院總體工程研究所,綿陽,621900

新型雙臂菱形壓電柔性機構(gòu)理論設(shè)計與建模

蔣 州1曹軍義1凌明祥1,2曾明華1林 京1

1.西安交通大學機械工程學院,西安,7100492.中國工程物理研究院總體工程研究所,綿陽,621900

針對目前傳統(tǒng)菱形位移放大機構(gòu)的輸出剛度、固有頻率和位移放大比無法同時提高的問題,提出了一種雙臂復合菱形柔性機構(gòu),并介紹了其理論設(shè)計與建模方法。基于歐拉-伯努利梁理論和卡氏第二定理,推導出雙臂復合菱形柔性機構(gòu)的位移放大比和剛度解析模型,利用拉格朗日方程建立了該柔性機構(gòu)的固有頻率解析模型,并通過商業(yè)有限元軟件驗證了解析模型的準確性。根據(jù)理論計算結(jié)果,通過電火花切割工藝加工出復合柔性機構(gòu),并與相同尺寸參數(shù)的傳統(tǒng)菱形位移放大機構(gòu)進行實驗比較研究。實驗結(jié)果表明,壓電雙層臂柔性機構(gòu)的固有頻率為1330 Hz,位移放大比為4.2,同時提高了固有頻率和放大比。此外,所建立的位移放大比和固有頻率力學解析模型可以為新型壓電柔性機構(gòu)的優(yōu)化設(shè)計提供理論指導。

壓電執(zhí)行器;微位移放大;柔性鉸鏈;三角形位移放大機構(gòu)

0 引言

近年來,隨著微電子[1]、光學微裝配以及超精密加工[2]等領(lǐng)域的迅猛發(fā)展,人們對精密定位技術(shù)提出了越來越高的要求。壓電陶瓷具有分辨率高、驅(qū)動力大、響應(yīng)速度快等優(yōu)點,成為精密定位平臺中使用最廣泛的驅(qū)動器件[3]。此外,壓電陶瓷在振動主動控制、精密閥、精密泵等領(lǐng)域也得到了廣泛的關(guān)注和應(yīng)用[4]。但是,壓電陶瓷最大的不足是輸出位移小,正常工作狀態(tài)下輸出位移僅是其自身尺寸的0.1%~0.2%。為實現(xiàn)亞毫米級的微位移輸出,工程實踐中常采用柔性位移放大機構(gòu)來滿足實際需求[5]。

微位移放大機構(gòu)的常用放大原理有杠桿放大、壓曲放大、液壓放大、三角形放大等[6]。張建瓴等[7]基于差式杠桿原理設(shè)計了微位移放大機構(gòu),并對該機構(gòu)放大倍數(shù)和靜態(tài)剛度進行理論建模。杠桿放大機構(gòu)的缺點是尺寸較大,輸出剛度和固有頻率較低。壓曲放大機構(gòu)利用薄板或較薄殼體受外力而產(chǎn)生的彎曲變形來進行位移放大,放大倍數(shù)與薄板的曲率以及薄板間的夾角有關(guān),放大倍數(shù)較小[8]。液壓放大機構(gòu)利用液體等體積原理,通過驅(qū)動活塞運動實現(xiàn)位移放大,主要應(yīng)用于精密伺服閥[9]。三角形放大機構(gòu)利用三角形三邊之間的幾何關(guān)系,通過改變長直角邊的位移來實現(xiàn)短直角邊位移的放大。在眾多位移放大機構(gòu)中,三角形位移放大機構(gòu)因其放大倍數(shù)大、結(jié)構(gòu)緊湊、固有頻率高、位移輸入輸出之間具有良好的線性關(guān)系等優(yōu)點,在科學研究和工程實際中得到廣泛應(yīng)用。

三角形位移放大機構(gòu)以菱形[3]、橋式(bridge-type)[10]、Moonie型[11]、Cymbal型[12]等為代表。目前,傳統(tǒng)的菱形位移放大機構(gòu)的輸出剛度、固有頻率和位移放大比三者之間無法同時達到最優(yōu),設(shè)計時一般是在這些指標中進行折中處理。為此,本文提出了一種雙層菱形柔性機構(gòu),其特點是結(jié)構(gòu)緊湊,在保持位移放大比的同時提高了機構(gòu)的輸出剛度和固有頻率,并采用卡氏第二定理和拉格朗日方程進行靜力學和動力學解析建模和理論分析,最后加工出實驗樣機,并進行實驗驗證。

1 新型雙臂菱形壓電柔性機構(gòu)

傳統(tǒng)的菱形位移放大機構(gòu)由于結(jié)構(gòu)緊湊,位移放大倍數(shù)高,抗干擾能力強等優(yōu)點得到了工程界的認可。但受限于菱形機構(gòu)放大比與固有頻率無法同時提高的缺陷,目前菱形放大機構(gòu)主要應(yīng)用于輸出剛度較小、頻響較低的工況。本文提出了一種具有雙臂的復合菱形位移放大機構(gòu),在保持較大位移放大比的同時,可以提高系統(tǒng)輸出剛度和固有頻率,其具體結(jié)構(gòu)如圖1所示。

(a)雙臂菱形柔性機構(gòu)實物圖

(b)雙臂菱形柔性機構(gòu)原理圖圖1 雙臂菱形柔性機構(gòu)Fig.1 Rhombus-type compliant mechanism with double flexible beams

整個機構(gòu)與傳統(tǒng)的菱形放大機構(gòu)構(gòu)型相似,上下、左右對稱布局,但是放大機構(gòu)的柔性臂采用雙層結(jié)構(gòu)。壓電陶瓷通過一定預緊力安裝在位移放大機構(gòu)水平對角線之間,在一定的輸入電壓下,壓電陶瓷輸出驅(qū)動力和水平位移,在其驅(qū)動作用下,菱形雙臂發(fā)生形變而使位移放大機構(gòu)產(chǎn)生豎直方向的輸出位移。合理優(yōu)化機構(gòu)參數(shù),可得到被放大的豎直方向位移輸出。機構(gòu)運動過程中保持菱形構(gòu)型不變,與傳統(tǒng)同參數(shù)的菱形放大機構(gòu)相比,雙臂菱形放大機構(gòu)在運動過程中幾何關(guān)系不變,大角度位移放大比幾乎不變;由單臂變?yōu)殡p臂結(jié)構(gòu),提高了機構(gòu)的輸出剛度和固有頻率。

2 新型柔性機構(gòu)的解析模型

2.1位移放大比模型

如圖2所示,假設(shè)壓電陶瓷輸出力和負載分別為FPZT和FN,其他幾何參數(shù)定義見圖1。當壓電陶瓷在正電壓作用下,位移放大機構(gòu)水平方向伸長2Δx時,其豎直方向?qū)⒖s減2Δy。由于機構(gòu)具有對稱性,取柔性位移放大機構(gòu)的1/4進行力學分析,則有如下關(guān)系式成立:FPZT=4fx;FN=4fy。如圖2所示,柔性雙臂機構(gòu)由柔性臂AB和柔性臂CD組成,因柔性臂AB與柔性臂CD的力學邊界條件以及幾何參數(shù)完全相同,為簡化力學分析過程,以柔性臂AB為例,對雙臂菱形位移放大機構(gòu)進行受力分析。考慮柔性臂的固定邊界條件:A、B兩點轉(zhuǎn)動角度相等,A點作水平移動,B點作豎直運動。由于受到A、B兩點的約束,柔性臂AB將受到2Mr力矩的作用,根據(jù)力平衡方程和力矩平衡方程,有

fAx=fBx=fx=FPZT/4

(1)

fAy=fBy=fy=FN/4

(2)

2Mr=fxlasinθ+fylacosθ

(3)

圖2 雙臂菱形機構(gòu)受力分析Fig.2 Force analysis of rhombus-type displacement amplifier

對簡化后的柔性臂AB進行受力分析。沿柔性臂方向,AB受到拉力的作用;垂直于柔性臂方向,AB受到一對力偶的作用,且A、B兩點的約束還導致AB存在力矩2Mr的作用,這將使柔性臂在運動過程中發(fā)生彎曲,進而輸出放大的位移。根據(jù)力和力矩平衡方程,沿AB長度方向、距A點距離為u的力和力矩分別為

(4)

基于歐拉-伯努利梁理論,柔性臂的總變形能為拉伸應(yīng)變能Vε1與彎曲應(yīng)變能Vε2之和,即

Vε=Vε1+Vε2

(5)

由卡氏第二定理,機構(gòu)沿x方向的位移

(6)

機構(gòu)沿y方向的位移

(7)

其中,EA和EI分別為柔性臂AB的抗拉剛度和抗彎剛度。基于式(6)和式(7),位移放大機構(gòu)空載時(fy=0)的位移放大比

(8)

位移放大機構(gòu)的輸入剛度

(9)

2.2固有頻率解析模型

與實際應(yīng)用工況一致,建立固有頻率解析模型時,位移放大機構(gòu)的底部邊界條件是固定的,各柔性環(huán)節(jié)的質(zhì)量和坐標選取如圖3所示。

圖3 雙臂菱形位移放大機構(gòu)及各部分質(zhì)量分布Fig.3 Mass distribution of the proposed mechanism

僅考慮柔性位移放大機構(gòu)豎直方向的輸出位移一個自由度,取廣義坐標q=2Δy。由前文力學分析和理論建模結(jié)果可知,水平輸入位移Δx、柔性臂轉(zhuǎn)動角度Δθ、柔性臂變形伸長量Δl與豎直輸出位移坐標q之間分別滿足如下關(guān)系:

(10)

(11)

式中,Kl、Kθ分別為柔性臂AB、CD的抗拉剛度和抗彎剛度。

考慮到柔性臂的抗拉剛度和抗彎剛度,則系統(tǒng)的彈性勢能

(12)

系統(tǒng)的動能

(13)

令拉格朗日函數(shù)L=T-V,對于自由振動,將式(12)和式(13)代入拉格朗日方程:

(14)

忽略推導過程,根據(jù)拉格朗日方程,可以得到自由振動的動力學方程:

(15)

雙臂菱形柔性機構(gòu)的固有頻率

(16)

3 有限元計算分析

3.1位移放大比仿真

為了驗證本文力學解析模型的準確性,保持雙臂菱形位移放大機構(gòu)的參數(shù)相同,將解析模型計算結(jié)果與有限元計算結(jié)果進行比較,并將有限元計算結(jié)果作為比較的基準。位移放大機構(gòu)具體幾何參數(shù)見表1,有限元計算的網(wǎng)格劃分如圖4所示。

進行位移放大比仿真計算時,完整約束圖4模型底部,柔性機構(gòu)水平對角線內(nèi)部兩側(cè)施加等大反向的驅(qū)動力;計算位移放大比仿真結(jié)果時,柔性機構(gòu)的位移放大比即仿真結(jié)果的豎直輸出位移與水平輸入位移的比值。計算得到的雙臂菱形位移放大比的部分結(jié)果見表2。

圖4 有限元模型Fig.4 Finite element model

參量數(shù)值參量數(shù)值la(mm)16.8d(mm)10h(mm)1.0E(GPa)210t(mm)1.5ρ(kg/m3)7850

表2 位移放大比計算結(jié)果

由表2可以看出,本文基于卡氏第二定理推導出的解析模型在角度θ較大時與有限元結(jié)果接近,在角度θ較小時還存在一定的誤差。存在計算誤差的主要原因是本文的雙層柔性臂是中空結(jié)構(gòu),在計算其彎矩時未考慮中空部分的影響。實際工程應(yīng)用中應(yīng)乘以工程系數(shù)進行調(diào)整,這樣,解析模型計算結(jié)果將與有限元計算結(jié)果更加接近。

3.2固有頻率仿真

工作頻寬和輸出位移放大倍數(shù)是評價位移放大機構(gòu)性能的兩個重要參數(shù),建立固有頻率與機構(gòu)幾何參數(shù)之間準確的解析力學模型是結(jié)構(gòu)優(yōu)化的前提。為了驗證本文建立的雙臂菱形位移放大機構(gòu)固有頻率解析模型的準確性,基于ANASYS軟件對同一模型進行了有限元模態(tài)分析,有限元計算的網(wǎng)格劃分如圖4所示,具體結(jié)構(gòu)參數(shù)見表1,其中,柔性機構(gòu)的角度θ為10°。

雙臂菱形位移放大機構(gòu)的前6階模態(tài)如圖5所示,其中,第1階模態(tài)描述了雙臂菱形位移放大機構(gòu)的平面外搖擺變形,頻率約1179 Hz;第2階模態(tài)為機構(gòu)繞Z軸扭轉(zhuǎn)變形,頻率約1257 Hz;第4階模態(tài)為機構(gòu)繞上下定點的扭轉(zhuǎn)變形,頻率約3806 Hz;第3階模態(tài)為位移放大機構(gòu)輸出位移方向的變形振型,頻率為1411.4 Hz,理論計算出該階固有頻率為1292.91 Hz,兩者的相對誤差為9.2%,誤差在工程許可范圍內(nèi)。

圖5 菱形位移放大機構(gòu)的前6階振型Fig.5 The first six vibration modes of the proposed mechanism

4 實驗驗證

按表1中的參數(shù)加工出雙臂菱形柔性機構(gòu)和傳統(tǒng)菱形柔性機構(gòu),其中,柔性機構(gòu)的角度θ為10°。雙臂菱形柔性機構(gòu)在傳統(tǒng)菱形柔性機構(gòu)的基礎(chǔ)上改進,除增加雙臂結(jié)構(gòu)外,其他尺寸參數(shù)與傳統(tǒng)菱形機構(gòu)相同。搭建實驗測試系統(tǒng)進行比較與驗證(圖6)。實驗中,由信號發(fā)生器(AFG 300)產(chǎn)生激勵信號,經(jīng)功率放大器放大后輸入壓電陶瓷,在壓電陶瓷激勵下,柔性機構(gòu)輸出微位移,該微位移由激光位移傳感器進行測量,位移傳感器所測位移和功率放大器輸出信號由示波器(Agilent,探頭內(nèi)阻為10 MΩ)進行同步采集。壓電陶瓷功率放大器型號為PI-E-482,電壓放大倍數(shù)為100;所用傳感器型號為HL-G103-A,分辨力為0.5 μm。

圖6 實驗測試系統(tǒng)Fig.6 Experimental set up

對壓電陶瓷施加正弦線性掃頻激勵,頻率范圍為1~2000 Hz,掃頻時間100 s,輸入電壓幅值20 V。實際輸入輸出曲線如圖7所示,可以看出在掃頻激勵67 s后,放大機構(gòu)出現(xiàn)共振現(xiàn)象。

圖7 掃頻輸入輸出信號Fig.7 Input and output sweep-frequency signal

放大機構(gòu)豎直方向的位移輸出由激光位移傳感器測得,分析放大機構(gòu)的輸入輸出得到系統(tǒng)的伯德圖(圖8)。實測雙臂菱形位移放大機構(gòu)固有頻率為1330 Hz,由有限元仿真分析對應(yīng)該階固有頻率為1411.4 Hz,誤差為5.77%,仿真與實驗誤差較小,表明仿真模型比較準確,可信度比較高。由圖8可以看出,本文提出的雙臂菱形柔性機構(gòu)與傳統(tǒng)的菱形機構(gòu)相比,固有頻率有了一定的提高;雙臂菱形柔性機構(gòu)的固有頻率為1330 Hz,傳統(tǒng)菱形機構(gòu)固有頻率為1140 Hz,機構(gòu)固有頻率提高了16.7%,輸出剛度增大了一倍,而位移放大比幾乎不變,充分表明了本文所設(shè)計機構(gòu)的優(yōu)點。

圖8 位移放大機構(gòu)掃頻分析伯德圖Fig.8 Bode diagram of different amplification mechanisms

5 結(jié)語

實驗測試結(jié)果表明,該雙臂菱形柔性機構(gòu)的固有頻率比傳統(tǒng)機構(gòu)提高16.7%,輸出剛度增大一倍,位移放大比為4.2,與傳統(tǒng)機構(gòu)接近。此外,基于本文建模方法所建立的位移放大比和固有頻率等力學解析模型可以準確預測位移放大機構(gòu)的靜動力學響應(yīng),揭示雙臂菱形柔性位移放大機構(gòu)各幾何參數(shù)對機構(gòu)靜、動力學行為的影響,可為機構(gòu)優(yōu)化設(shè)計和控制設(shè)計提供較為準確的依據(jù)。

[1] LING Mingxiang, CAO Junyi, JIANG Zhou, et al. Theoretical Modeling of Attenuated Displacement Amplification for Multistage Compliant Mechanism and Its Application [J]. Sensors and Actuators A: Physical,2016,15(22):15-22.

[2] 黃衛(wèi)清,史小慶,王寅,等. 菱形壓電微位移放大機構(gòu)的設(shè)計[J].光學精密工程,2015,23(3):803-809.

HUANG Weiqing, SHI Xiaoqing, WANG Yin, et al. Design of Diamond Piezoelectric Micro-displacement Amplification Mechanism [J]. Optics and Precision Engineering,2015,23(3):803-809 .

[3] 凌明祥,曹軍義,曾明華,等. 壓電柔性機構(gòu)位移放大比精確建模[J].智能材料與結(jié)構(gòu),2016,25:075022.

LING Mingxiang, CAO Junyi, ZENG Minghua, et al. Enhanced Mathematical Modeling of the Displacement Amplification Ratio for Piezoelectric Compliant Mechanisms [J]. Smart Materials and Structures,2016,25:075022.

[4] 凌明祥,劉謙,曹軍義,等.壓電位移放大機構(gòu)的力學解析模型及有限元分析[J].光學精密工程,2016,24(4):812-818.

LING Mingxiang, LIU Qian, CAO Junyi, et al.Analytical Model and Finite Element Analysis of Piezoelectric Displacement Amplification Mechanism[J]. Optics and Precision Engineering,2016,4(4):812-818 .

[5] 田延嶺,張大衛(wèi),閆兵,等.二自由度微定位平臺的研制[J].光學精密工程,2006,14(1):94-99.

TIAN Yanling, ZHANG Dawei, YAN Bing, et al. Development of a 2-DOF Micro-positioning Table [J]. Optics and Precision Engineering,2006,14(1):94-99.

[6] 林超,俞松松,程凱,等.微/納米定位平臺的動態(tài)特性分析與試驗[J].浙江大學學報(工學版),2012,46(8):1375- 1381.

LIN Chao, YU Songsong, CHENG Kai, et al. Dynamic Analysis and Testing of Micro/nano-positioning Platform [J].Journal of Zhejiang University(Engineering Science),2012,46(8):1375-1381.

[7] 張建瓴,陳萬銀,可欣榮,等.一種微位移放大機構(gòu)的設(shè)計與仿真[J].機械設(shè)計,2009,26(12):9-12.

ZHANG Jianling, CHEN Wanyin, KE Xinrong,et al.Design and Simulation of a Kind of Micro-displacement Amplification of a Kind of Micro-displacement Amplification Mechanism [J]. Journal of Machine Design,2009,26(12):9-12.

[8] 王興松,王湘江,毛燕,等.基于超磁致伸縮材料的折彎型壓曲放大機構(gòu)設(shè)計、分析與控制[J].機械工程學報,2007,43(11):27-33.

WANG Xingsong, WANG Xiangjiang, MAO Yan, et al. Design, Analysis and Control of Novel Pressing and Bending Magnifying Mechanism Driven by Giant Magnetostrictive Material[J]. Chinese Journal of Mechanical Engineering,2007,43(11):27-33.

[9] 俞軍濤,焦宗夏,吳帥.基于液壓微位移放大結(jié)構(gòu)的新型壓電陶瓷直接驅(qū)動閥設(shè)計及仿真[J].機械工程學報,2013,49(2):151-158.

YU Juntao, JIAO Zongxia, WU Shuai. Design and Simulation Study on New Servo Valve Direct Driven by Piezoelectric Actuator Using Hydraulic Amplification[J]. Journal of Mechanical Engineering,2013,49(2):151-158.

[10] 馬立,謝煒,劉波,等.柔性鉸鏈微定位平臺的設(shè)計[J].光學精密工程,2014,22(2):338-345.

MA Li, XIE Wei, LIU Bo, et al. Design of Micro-positioning Stage with Flexure Hinge[J]. Optics and Precision Engineering,2014,22(2):338-345.

[11] FREDERIC L, ZAFFIR C, CRAIG A R. A Simplified Geometrically Nonlinear Approach to the Analysis of the Moonie Actuator[J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control,1995,59(1):21-27.

[12] AYDIN D, KENJI U. Composite Piezoelectric Transducer with Truncated Conical Endcaps “Cymbal”[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control,1997,44(3):597-605.

TheoreticalDesignandModelingofaNovelPiezo-drivenRhombicFlexureMechanismwithDoubleBeams

JIANG Zhou1CAO Junyi1LING Mingxiang1,2ZENG Minghua1LIN Jing1

1.School of Mechanical Engineering,Xi’an Jiaotong University,Xi’an,710049 2.Institute of System Engineering,China Academy of Engineering Physics,Mianyang,Sichuan,621900

Design and modeling of a new diamond micro-displacement mechanism with double flexible beams was presented to overcome the issue that the output stiffness, natural frequency and displacement amplification ratio of the traditional diamond displacement amplification mechanism might not be improved at the same time. The static analysis model of displacement amplification ratio and input/output stiffness of the mechanism were derived based on Castigliano’s second theorem. Analytical model of natural frequency was also deduced by further employing the Lagrange equation .The finite element analysis was used to verify the accuracy of the analytical models. A prototype of the proposed compliant mechanism was machined by means of electric spark cutting processes and the performance measurement was performed for comparative analysis with the traditional mechanism with the same dimensions. Experimental results show that the natural frequency of the composite compliant mechanism is as 1330 Hz and displacement amplification ratio is as 4.2. In addition, the modeling method and the corresponding theoretical formula proposed herein may provide a useful reference for the optimal design and development of similar displacement amplification mechanisms.

piezoelectric actuator; displacement amplification; flexure hinge; triangle displacement amplification mechanism

TH113.2;TN384

10.3969/j.issn.1004-132X.2017.21.007

2016-12-05

國家自然科學基金資助項目(51575426,51421004)

(編輯陳勇)

蔣州,男,1993年生。西安交通大學機械工程學院碩士研究生。主要研究方向為壓電智能結(jié)構(gòu)建模與振動控制。E-mail:jiangzhou_xy@163.com。曹軍義,男,1977年生。西安交通大學機械工程學院副教授、博士研究生導師。凌明祥,男,1986年生。西安交通大學機械工程學院博士研究生,中國工程物理研究院總體工程研究所工程師。曾明華,男,1993年生。西安交通大學機械工程學院碩士研究生。林京,男, 1971年生。西安交通大學機械工程學院教授、博士研究生導師。

猜你喜歡
有限元模型
一半模型
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
新型有機玻璃在站臺門的應(yīng)用及有限元分析
基于有限元的深孔鏜削仿真及分析
基于有限元模型對踝模擬扭傷機制的探討
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
磨削淬硬殘余應(yīng)力的有限元分析
基于SolidWorks的吸嘴支撐臂有限元分析
主站蜘蛛池模板: 日本在线国产| 亚洲成年网站在线观看| 亚洲精品高清视频| 精品综合久久久久久97| 玖玖精品在线| 久草青青在线视频| 欧美日韩国产精品va| 久久先锋资源| 亚洲精品无码专区在线观看| 午夜视频www| 国产精品成人AⅤ在线一二三四| 色AV色 综合网站| 国产原创演绎剧情有字幕的| 无码AV日韩一二三区| 美女视频黄频a免费高清不卡| 欧美特级AAAAAA视频免费观看| 亚洲国产成熟视频在线多多| 五月天久久婷婷| 欧美一级色视频| a级毛片免费看| 欧美国产精品拍自| 高清欧美性猛交XXXX黑人猛交 | 婷婷中文在线| 国产性生交xxxxx免费| 午夜国产理论| 这里只有精品在线播放| 欧美日韩理论| 91麻豆久久久| 久爱午夜精品免费视频| 538国产在线| 亚洲无码熟妇人妻AV在线| 亚洲精品无码AV电影在线播放| 57pao国产成视频免费播放| 成人韩免费网站| 精品99在线观看| 久久这里只精品国产99热8| 亚亚洲乱码一二三四区| 国产电话自拍伊人| av在线手机播放| 91视频免费观看网站| 国产乱人免费视频| 国产一区二区精品福利| 超清无码一区二区三区| 天堂av综合网| 一本大道视频精品人妻 | 国产一区在线视频观看| 亚洲精品成人片在线观看| 国产乱肥老妇精品视频| 好吊妞欧美视频免费| 999精品色在线观看| 色欲不卡无码一区二区| 久久综合五月| 国产午夜精品一区二区三区软件| 亚洲第一成年免费网站| 免费毛片视频| 久久精品国产精品一区二区| 乱码国产乱码精品精在线播放| 999福利激情视频| 亚洲成年人片| 99re66精品视频在线观看| 狠狠做深爱婷婷久久一区| 四虎精品国产AV二区| 91综合色区亚洲熟妇p| 欧美日韩一区二区三区在线视频| 波多野结衣中文字幕一区二区| 99精品在线看| 日本一区二区三区精品AⅤ| 亚洲欧美日韩成人在线| 日本午夜影院| 亚洲人成电影在线播放| 免费在线一区| 国产精品永久在线| 欧美成人免费午夜全| 欧美在线伊人| 99热这里只有免费国产精品| 91视频国产高清| 无码人妻热线精品视频| 91久久国产热精品免费| 凹凸精品免费精品视频| 精品国产成人av免费| 亚洲综合九九| 国产成人三级在线观看视频|