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

    離子型金屬聚合物復(fù)合材料的分?jǐn)?shù)階建模及二階滑??刂?/h1>
    2020-04-22 10:37:52陳靖威黃梁松李甜甜李連鵬李玉霞
    科學(xué)技術(shù)與工程 2020年7期
    關(guān)鍵詞:致動器二階控制算法

    陳靖威, 黃梁松, 李甜甜, 李連鵬, 李玉霞

    (山東科技大學(xué)山東省機(jī)器人與智能技術(shù)重點(diǎn)實(shí)驗(yàn)室,青島266590)

    離子型金屬聚合物復(fù)合材料 (IPMC)[1]是一種電驅(qū)動智能復(fù)合材料,能夠在外加電場的驅(qū)動下實(shí)現(xiàn)大幅度的彎曲形變,撤去電場后又能快速恢復(fù)至初始狀態(tài)。近年來,由于其作為致動器具有致動電壓低(1~5 V),響應(yīng)速度快,擺動位移大,可根據(jù)需要裁剪成任意形狀且柔韌性好等特點(diǎn),引起大量國內(nèi)外學(xué)者進(jìn)行了深入的研究與探索。美國NASA(National Aeronautics and Space Administration) 的 JPL(Jet Propulsion Laboratory)實(shí)驗(yàn)室用 IPMC 驅(qū)動器制造了與自動刮雨器相類似的除塵刷,用于太空儀器表面的除塵[2];美國密歇根州立大學(xué)的Tant等設(shè)計(jì)了一種用 IPMC 驅(qū)動的機(jī)器魚[3];日本Eamex公司開發(fā)了以 IPMC 作為驅(qū)動材料的商業(yè)化產(chǎn)品人工魚[4];日本的名古屋大學(xué)研制了基于 IPMC 材料的泳動機(jī)器人[5]。因此,IPMC在仿生學(xué)、生物醫(yī)學(xué)和微機(jī)電等領(lǐng)域的應(yīng)用前景十分廣闊。

    分?jǐn)?shù)階系統(tǒng),即是用分?jǐn)?shù)階概念里的導(dǎo)數(shù)和積分算子來代替?zhèn)鹘y(tǒng)整數(shù)階動力系統(tǒng)當(dāng)中的導(dǎo)數(shù)和積分,近年來受到了越來越多的關(guān)注。國外學(xué)者Sabatier指出采用分?jǐn)?shù)階模型可以有效地模擬大量的物理現(xiàn)象[6]。針對IPMC模型,Bao等證明了IPMC具有分形電極結(jié)構(gòu)[7]。并由Caponetto在2008年確定了能夠在較大的頻率范圍內(nèi)描述IPMC執(zhí)行器動態(tài)特性的分?jǐn)?shù)階傳遞函數(shù)[8]。證明了IPMC可以用具有長期時(shí)間或空間依賴現(xiàn)象特性的分?jǐn)?shù)階系統(tǒng)來表示,即建立一個(gè)分?jǐn)?shù)階IPMC模型能更好地描述它的運(yùn)動特性。

    在大多數(shù)情況下,滑模是通過向控制或觀測系統(tǒng)中注入一個(gè)非線性不連續(xù)項(xiàng)來獲得的,該項(xiàng)取決于輸出誤差。在設(shè)計(jì)不連續(xù)注入時(shí),系統(tǒng)的運(yùn)動軌跡必須保持在誤差空間的某個(gè)滑動面上,產(chǎn)生的運(yùn)動稱為滑模[9]?;?刂品椒ㄓ捎谄渚哂杏邢迺r(shí)間的收斂性和對不確定情況的強(qiáng)魯棒性等優(yōu)點(diǎn),得到了廣泛的應(yīng)用。本文根據(jù)Davila在2005年提出的二階滑??刂朴^測器[10],設(shè)計(jì)了一個(gè)基于分?jǐn)?shù)階的二階滑??刂破?,從理論方面論證了其可行性,并在仿真與實(shí)驗(yàn)兩個(gè)方面分別與分?jǐn)?shù)階PIλDμ控制方法作對比,證明所提出的控制器具有更好的控制精度。

    1 實(shí)驗(yàn)平臺的搭建

    離子聚合物是IPMC的重要組成部分,也是用來區(qū)分膜體種類的重要依據(jù),主要分為全氟磺酸離子聚合物(Nafion)和全氟碳酸離子聚合物(Flemion)兩種,通常選用美國杜邦公司生產(chǎn)的Nafion膜作為IPMC的基底膜。除了基膜外,還有上下表面的電極,通常為各種貴金屬材料,如鉑、鈀、金和銀等。值得注意的是,電極和基膜的界面特性對IPMC致動器的變形特性十分重要,通??梢詫﹄x子交換膜進(jìn)行打磨處理,然后對糙化后的基膜采用浸泡還原的方法沉積電極層,所獲得的電極與基膜相互滲透,可以極大提高電極與基膜的接觸面積。圖1為IPMC致動器結(jié)構(gòu)圖。

    圖1 IPMC致動器結(jié)構(gòu)圖Fig.1 IPMC actuator structure diagram

    實(shí)驗(yàn)平臺如圖2所示,使用LABVIEW搭建了上位機(jī)部分并通過RS485與控制器進(jìn)行串口通信,可以實(shí)現(xiàn)IPMC性能測試和輸出位移控制兩種功能。

    圖2 IPMC致動器實(shí)驗(yàn)平臺Fig.2 IPMC actuator experimental platform

    本實(shí)驗(yàn)平臺的硬件主要包括控制器芯片、電源模塊、可調(diào)電壓輸出模塊、數(shù)據(jù)轉(zhuǎn)換傳輸模塊、激光位移傳感器等。

    (1)控制器芯片選用TI公司的TMS320C28335浮點(diǎn)DSP(digital signal processing),該控制器相比于以往的DSP,多了一個(gè)MAC單元,具有運(yùn)算速度快、數(shù)據(jù)及程序存儲量大等優(yōu)勢,可滿足實(shí)驗(yàn)平臺功能需要。

    (2)電源模塊主要為開關(guān)電源直接供電220 V轉(zhuǎn)24 V,為激光位移傳感器供電24 V,為控制器芯片供電+3.3 V和+1.9 V,為IPMC致動器提供可調(diào)的輸入電壓。

    (3)可調(diào)電壓輸出模塊是為了滿足實(shí)驗(yàn)中IPMC致動器兩極的電壓不僅要有直流電壓,還要有正弦電壓。所以為了滿足需要,該部分采用了正弦波脈寬調(diào)制法(SPWM),該方法中使用的驅(qū)動電路是全橋驅(qū)動芯片HIP4082IB,直流可調(diào)降壓電源為LM2596S。

    (4)數(shù)據(jù)轉(zhuǎn)換傳輸模塊中,使用的模數(shù)轉(zhuǎn)換芯片為ADS8688,吞吐量可達(dá)500 kSPS;數(shù)字隔離器為ISO7241M,它具有四通道配置和輸出使能功能,傳輸速度最高可達(dá)150 Mbit/s。

    (5)為了實(shí)現(xiàn)IPMC致動器尖端位移的檢測,本實(shí)驗(yàn)平臺選用美國邦納公司生產(chǎn)的L-GAGE LE250激光位移傳感器作為檢測設(shè)備。它帶有顯示屏,在量程范圍(位移100~400 mm,電流0~10 V)內(nèi)可以實(shí)時(shí)顯示當(dāng)前的測量位移和輸出電流。

    使用該實(shí)驗(yàn)平臺,對IPMC薄膜施加一個(gè)電壓為2 V,頻率為10/2π Hz的正弦電壓,可得其尖端位移的曲線如圖3所示。

    圖3 驅(qū)動電壓2 V、10/2π Hz時(shí)的尖端位移曲線Fig.3 Tip displacement curve when driving voltage is 2 V and 10/2π Hz

    并且可以得到在該正弦電壓下的輸入輸出關(guān)系圖如圖4所示,可以清楚地看出IPMC具有的磁滯特性。

    圖4 IPMC的輸入輸出關(guān)系Fig.4 The input/output relationship of IPMC

    2 IPMC致動模型的建立

    2.1 分?jǐn)?shù)階微積分

    早在300多年前,隨著古典微積分概念的提出,數(shù)學(xué)家們就推導(dǎo)出了分?jǐn)?shù)階微積分,雖然僅僅是純數(shù)學(xué)理論分析和推導(dǎo),但也受到了眾多學(xué)者的青睞。直到19世紀(jì)后期,分?jǐn)?shù)階微積分理論才開始逐漸在實(shí)際工程領(lǐng)域中得到初步的應(yīng)用[11]。近年來,由于運(yùn)用傳統(tǒng)控制理論和方法無法達(dá)到期望指標(biāo),更無法滿足人們對控制精度日益提高的要求,激發(fā)了學(xué)者們將分?jǐn)?shù)階微積分理論應(yīng)用在實(shí)際工程的研究熱情以及對分?jǐn)?shù)階領(lǐng)域進(jìn)行更加深入的探索。

    (1)

    分?jǐn)?shù)階微積分普遍定義如下。

    Riemann-Liouville 分?jǐn)?shù)階積分定義:存在α>0,t≥0時(shí),f(t)被定義為

    (2)

    式(2)中:Γ表示歐拉伽馬函數(shù):

    (3)

    (4)

    Caputo 分?jǐn)?shù)階微分定義:存在α>0,t≥0時(shí),f(t)被定義為

    (5)

    當(dāng)α∈(0,1)時(shí),Riemann-Liouville與Caputo分?jǐn)?shù)階微分之間的聯(lián)系為

    (6)

    2.2 IPMC分?jǐn)?shù)階建模

    Sun在2014年的時(shí)候提出了一個(gè)基于IPMC的傳遞函數(shù)模型,分別采用了線性不可逆電動力學(xué)模型和動力學(xué)模型對離子聚合物-金屬復(fù)合材料的性能進(jìn)行了分析[12]。其中第一個(gè)模型根據(jù)Onsager方程準(zhǔn)確地預(yù)測了其靜態(tài)特性,而第二個(gè)模型能夠揭示其磁滯特性。通過這兩個(gè)模型推導(dǎo)出的IPMC動力學(xué)傳遞函數(shù)模型將能較好地?cái)M合出IPMC的實(shí)際運(yùn)動效果。

    該傳遞函數(shù)模型如下:

    (7)

    式(7)中:C、τ1、τ2、K1、K2為需要辨識的參數(shù)。

    通過上述的實(shí)驗(yàn)平臺,對IPMC采用驅(qū)動電壓為2 V,頻率為10/2π Hz的正弦信號,得到所需的輸入輸出數(shù)據(jù),隨后結(jié)合人工蜂群算法來對上述傳遞函數(shù)中的5個(gè)未知參數(shù)進(jìn)行辨識。人工蜂群算法是由土耳其學(xué)者提出的一種基于蜜蜂群體覓食行為的智能優(yōu)化算法[13]。該算法相比粒子群算法、蟻群算法等其他智能優(yōu)化算法,具有參數(shù)少、易于實(shí)現(xiàn)、魯棒性較好等優(yōu)點(diǎn),而且全局搜索和局部搜索在每次迭代過程中都會進(jìn)行,這一點(diǎn)不僅增加了找到最優(yōu)解的概率,還相對減小了陷入局部最優(yōu)的可能性。整個(gè)優(yōu)化步驟如下。

    (1)初始化ABC(artificial bee colony)算法中的各種參數(shù)、食物源的數(shù)量SN、食物源的最大循環(huán)次數(shù)和終止循環(huán)次數(shù)。

    本次人工蜂群優(yōu)化算法中,設(shè)蜜源數(shù)N=100,參數(shù)limit為100,最大迭代次數(shù)為500,程序運(yùn)行次數(shù)為100。最終得到5個(gè)待辨識的參數(shù)分別為:τ1=0.145 2,τ2=12.722 4,C=1,K1=1.275,K2=-22.246 8。即傳遞函數(shù)模型為

    (8)

    得到了IPMC整數(shù)階的傳遞函數(shù),此時(shí)再根據(jù)上述的輸入輸出數(shù)據(jù),采用結(jié)合了泰勒計(jì)數(shù)法和梯度法的Marquardt算法將其轉(zhuǎn)化成分?jǐn)?shù)階模型[14]。該算法實(shí)際上是使模型滿足:

    (9)

    (10)

    (11)

    根據(jù)圖5可以看出分?jǐn)?shù)階模型與實(shí)驗(yàn)數(shù)據(jù)的擬合度更高,即分?jǐn)?shù)階模型能夠描述出IPMC所具有的非整數(shù)階特性。

    圖5 分?jǐn)?shù)階模型、整數(shù)階模型和實(shí)驗(yàn)輸出結(jié)果比較Fig.5 Comparison of fractional order model, integer order model and experimental output results

    3 基于分?jǐn)?shù)階的二階滑模控制器設(shè)計(jì)

    考慮到一個(gè)受動態(tài)模型控制的帶擾動的一階分?jǐn)?shù)階系統(tǒng):

    CDαx(t)=ax(t)+bu(t)+φ(t)

    (12)

    式(12)中:α∈(0,1)是已知的非整數(shù)階微分;x(t)∈R是測量出的輸出變量;u(t)∈R是一個(gè)可調(diào)輸入;a∈R是一個(gè)不確定的參數(shù);b是一個(gè)已知的正常量;φ(t)是一個(gè)不確定的干擾。

    基于上述分?jǐn)?shù)階系統(tǒng),設(shè)計(jì)了一個(gè)基于分?jǐn)?shù)階的二階滑??刂破鳛?/p>

    (13)

    設(shè)xr(t)為一個(gè)足夠光滑的參考面,此時(shí),需要證明該控制器能夠保證跟蹤誤差能在有限時(shí)間歸零。誤差e可以表示為

    e(t)=x(t)-xr(t)

    (14)

    考慮跟蹤誤差變量的動態(tài)特性:

    CDαe(t)=ax(t)+bu(t)+φ(t)-CDαxr(t)

    (15)

    此時(shí)把除了可調(diào)輸入的其余三項(xiàng)合并為

    Δ (x,xr,t)=ax(t)+φ(t)-CDαxr(t)

    (16)

    很容易可以證得閉環(huán)跟蹤誤差動力學(xué)滿足這種形式:

    (17)

    定義:

    (18)

    此時(shí):

    Δ2(xr,t)=φ(t)-CDαxr(t)

    (19)

    (20)

    根據(jù)式(19)可得:

    (21)

    此時(shí),根據(jù)式(12),假設(shè)受控輸出x(t)的時(shí)間演化足夠光滑,即系統(tǒng)輸出x(t)由一個(gè)已知的常量xd和有限時(shí)間td組成,即:

    |RLD2-αx(t)|≤xd,t≥td

    (22)

    在該假設(shè)的基礎(chǔ)上,擾動φt,參考面xr(t)以及不確定的常量a有以下關(guān)系:

    (23)

    式(23)中:xa、xb、xc為已知的常量;tb、tc為有限時(shí)間。

    根據(jù)式(22)~式(23),推得:

    RLD2-αΔ2(xr,t)≤xb+xc

    (24)

    針對所提出的控制器,其整數(shù)階模型和分?jǐn)?shù)階模型在有限時(shí)間歸零條件下是等價(jià)的,所以通過形式更加簡單的整數(shù)階模型進(jìn)行李雅普諾夫定理的證明,簡化后的整數(shù)階模型為

    (25)

    轉(zhuǎn)換Lyapunov函數(shù)為

    (26)

    式(26)在任何地方都連續(xù),但是在x1=0處不可微。

    因?yàn)橛墒?25)得到的微分子集的狀態(tài)軌跡φ(t,x0)是絕對連續(xù)函數(shù),因此V[φ(t,x0)]是一個(gè)連續(xù)時(shí)間的函數(shù),但于V(x)缺少Lipschitz性質(zhì),所以無法保證V[φ(t,x0)]的絕對連續(xù)性和它的可微性是無處不在的。然而V(x)是除了在S={(x1,x2)∈R2|x1=0}處是連續(xù)可微的。很容易看出系統(tǒng)(13)的軌跡只是穿過曲面S,不能停留在上面,除非達(dá)到原點(diǎn)x=0。這意味著V[φ(t,x0)]幾乎在每一個(gè)t都是可微的,并且這些點(diǎn)導(dǎo)數(shù)都可以用鏈?zhǔn)椒▌t這種普通的計(jì)算方法得到。這表明,在運(yùn)用Lyapunov定理時(shí),可以只考慮V(x)可微的點(diǎn)。

    (27)

    注意V(x)在x=0處是連續(xù)不可微的。當(dāng)K3>0時(shí)為正定無界,即:

    (28)

    (29)

    (30)

    (31)

    4 仿真和實(shí)驗(yàn)研究

    4.1 仿真研究

    選取了目前在分?jǐn)?shù)階控制系統(tǒng)中較為常用的分?jǐn)?shù)階PIλDμ控制算法與所提出的基于分?jǐn)?shù)階的二階滑??刂扑惴ㄗ鲗Ρ取7?jǐn)?shù)階PIλDμ控制算法是分?jǐn)?shù)階微積分理論與傳統(tǒng)PID控制算法的融合,在傳統(tǒng)PID控制的三個(gè)參數(shù)KP、KI、KD的基礎(chǔ)上,又多了兩個(gè)可調(diào)的參數(shù)λ和μ,所以控制器參數(shù)的正定范圍變大,能夠更加靈活地對受控對象進(jìn)行控制,得到更好的控制效果[15]。設(shè)計(jì)PIλDμ控制器如下:

    GFOPID(s)=0.95+2.5s0.98+0.04s-1.9

    (32)

    圖7所示即為上述仿真系統(tǒng)中的單位負(fù)反饋階躍響應(yīng)。其中FOSM表示基于分?jǐn)?shù)階的二階滑模控制,F(xiàn)OPID表示分?jǐn)?shù)階PIλDμ控制。

    由圖7、表1可見,在上升時(shí)間、超調(diào)量、調(diào)節(jié)時(shí)間三個(gè)方面,基于分?jǐn)?shù)階的二階滑??刂贫家∮诜?jǐn)?shù)階PIλDμ控制。即從仿真結(jié)果來說,基于分?jǐn)?shù)階的二階滑??刂破鞯目刂菩Ч獌?yōu)于分?jǐn)?shù)階PIλDμ控制器。

    圖6 分?jǐn)?shù)階PIλDμ與基于分?jǐn)?shù)階的二階滑??刂瓶驁DFig.6 The control block diagrams of fractional PIλDμ and second order sliding model based on fractional order

    圖7 FOSM與FOPID的比較Fig.7 Comparison of FOSM and FOPID

    表1 兩種控制器的數(shù)據(jù)對比

    4.2 實(shí)驗(yàn)研究

    已經(jīng)通過第一章介紹的實(shí)驗(yàn)平臺對IPMC致動器進(jìn)行了性能測試,得到了實(shí)驗(yàn)數(shù)據(jù),下面則運(yùn)用該實(shí)驗(yàn)平臺對IPMC尖端位移進(jìn)行控制。通過上位機(jī)發(fā)送出期望位移,RS485接收到該位移并進(jìn)行相應(yīng)的控制器的計(jì)算,然后將計(jì)算得到電壓附加到IPMC致動器上,以實(shí)現(xiàn)對IPMC致動器的控制。軟件流程圖如圖8所示。

    圖8 控制模式軟件流程圖Fig.8 Software flow chart of control mode

    由文獻(xiàn)[14]可知,分?jǐn)?shù)階PIλDμ控制算法中的微分項(xiàng)對IPMC控制的影響很小,所以在本實(shí)驗(yàn)中忽略其參數(shù)KD和μ??梢酝ㄟ^上位機(jī)隨時(shí)更改分?jǐn)?shù)階PIλ控制算法中的KP、KI、λ三個(gè)參數(shù)以及基于分?jǐn)?shù)階的二階滑??刂扑惴ㄖ械摩?、K1、K3三個(gè)參數(shù)。通過對參數(shù)不斷調(diào)節(jié),以期望達(dá)到每種算法的最優(yōu)跟蹤效果。設(shè)定期望位移為幅值為3 mm,頻率為1/2π Hz的正弦曲線。兩種算法的最優(yōu)跟蹤曲線及誤差分別如圖9~圖12所示。

    圖9 分?jǐn)?shù)階PIλ控制時(shí)的輸出曲線Fig.9 The output curve under fractional PIλcontrol

    圖10 分?jǐn)?shù)階PIλ控制器的控制誤差Fig.10 The control error of fractional PIλcontroller

    圖11 基于分?jǐn)?shù)階的二階滑??刂破鞯妮敵銮€Fig.11 The output curve under second order sliding mode controller based on fractional order

    圖12 基于分?jǐn)?shù)階的二階滑??刂破鞯目刂普`差Fig.12 The control error of second order sliding mode controller based on fractional order

    由圖9~圖12可以看出,基于分?jǐn)?shù)階的二階滑模制器的跟蹤效果要好于分?jǐn)?shù)階PIλDμ控制器,且位移誤差更小。通過計(jì)算均方根誤差可得,基于分?jǐn)?shù)階的二階滑模制控制器的均方根誤差為RMSE=0.289 5,分?jǐn)?shù)階PIλDμ控制器的均方根誤差為RMSE=0.431 9。因此可以得出,所提出的控制器在針對IPMC實(shí)際的控制效果上要優(yōu)于傳統(tǒng)的分?jǐn)?shù)階PIλDμ控制器,且擁有較好的控制效果。

    5 結(jié)論

    根據(jù)IPMC特性,搭建了實(shí)驗(yàn)平臺,并為IPMC薄膜施加了驅(qū)動電壓為2 V,頻率為10/2π Hz的正弦電壓,得到了實(shí)驗(yàn)數(shù)據(jù)。隨后,選用人工蜂群算法對參考模型中的未知參數(shù)進(jìn)行尋優(yōu),得到了整數(shù)階模型;再使用Marquardt算法辨識出了分?jǐn)?shù)階模型。通過兩種模型與實(shí)驗(yàn)數(shù)據(jù)的比較可以得出,分?jǐn)?shù)階模型與實(shí)驗(yàn)數(shù)據(jù)擬合度更高,能更精確地反映出IPMC的非整數(shù)階特性。

    然后針對該IPMC的分?jǐn)?shù)階模型,設(shè)計(jì)了一個(gè)基于分?jǐn)?shù)階的二階滑??刂扑惴?,并從理論方面證明了其合理性。后將其與常用于分?jǐn)?shù)階系統(tǒng)控制的分?jǐn)?shù)階PIλDμ算法進(jìn)行對比,并分別采用仿真與實(shí)驗(yàn)兩個(gè)方面進(jìn)行。無論是仿真結(jié)果還是實(shí)驗(yàn)結(jié)果,所提出的控制算法均優(yōu)于分?jǐn)?shù)階PIλDμ算法。即證明了所提方法是可行的,可以更加精確地控制IPMC致動器。

    通過實(shí)驗(yàn)過程可知,控制器參數(shù)不僅關(guān)系著控制器本身,而且對于分?jǐn)?shù)階系統(tǒng)的影響也很大。所以,下一步應(yīng)該針對所提出控制器的參數(shù)進(jìn)行深入的研究,提供一種有效的參數(shù)選取方法,以便在任何控制情況下都能高效地設(shè)計(jì)出所需的控制器。

    猜你喜歡
    致動器二階控制算法
    一類二階迭代泛函微分方程的周期解
    一類二階中立隨機(jī)偏微分方程的吸引集和擬不變集
    二階線性微分方程的解法
    MEMS致動器結(jié)構(gòu)參數(shù)變化對致動器的影響
    一類二階中立隨機(jī)偏微分方程的吸引集和擬不變集
    基于ARM+FPGA的模塊化同步控制算法研究
    超磁致伸縮致動器的復(fù)合反饋控制及其在變橢圓銷孔精密加工中的應(yīng)用
    精密磁致伸縮致動器的動態(tài)非線性多場耦合建模
    一種優(yōu)化的基于ARM Cortex-M3電池組均衡控制算法應(yīng)用
    FE-GA-AL磁致伸縮致動器輸出特性研究

    a级毛片a级免费在线| 无限看片的www在线观看| 一进一出抽搐动态| 麻豆av在线久日| 久久中文字幕人妻熟女| 88av欧美| 亚洲国产精品sss在线观看| 曰老女人黄片| 国产伦在线观看视频一区| 国产亚洲精品av在线| 日日摸夜夜添夜夜添小说| 亚洲av电影不卡..在线观看| 日本在线视频免费播放| 国产97色在线日韩免费| 91国产中文字幕| 十八禁人妻一区二区| 宅男免费午夜| 99久久精品热视频| 黄色a级毛片大全视频| av欧美777| 亚洲黑人精品在线| 婷婷六月久久综合丁香| 在线观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲精品中文字幕在线视频| 妹子高潮喷水视频| 成人特级黄色片久久久久久久| 动漫黄色视频在线观看| 欧美性长视频在线观看| 国产成人av激情在线播放| 男人的好看免费观看在线视频 | 久久精品国产综合久久久| 亚洲av中文字字幕乱码综合| 亚洲精品av麻豆狂野| 一本精品99久久精品77| 美女扒开内裤让男人捅视频| 国产精品亚洲av一区麻豆| 一级作爱视频免费观看| av片东京热男人的天堂| 99热6这里只有精品| 久久这里只有精品19| 国产视频内射| 中文字幕精品亚洲无线码一区| 99久久精品国产亚洲精品| tocl精华| 国产真实乱freesex| 亚洲人成电影免费在线| 99热6这里只有精品| 国产视频内射| 国产成人啪精品午夜网站| 国内久久婷婷六月综合欲色啪| 午夜成年电影在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 两人在一起打扑克的视频| 国内精品久久久久久久电影| 首页视频小说图片口味搜索| 成年免费大片在线观看| 午夜老司机福利片| 免费在线观看亚洲国产| 婷婷丁香在线五月| 可以免费在线观看a视频的电影网站| 一进一出抽搐动态| 国产真实乱freesex| 午夜精品在线福利| 日本一二三区视频观看| 青草久久国产| 成人av在线播放网站| 真人一进一出gif抽搐免费| 成人av在线播放网站| 99久久精品国产亚洲精品| 老司机午夜福利在线观看视频| aaaaa片日本免费| 特大巨黑吊av在线直播| 亚洲第一电影网av| 中文字幕久久专区| 成人一区二区视频在线观看| 可以在线观看毛片的网站| 好看av亚洲va欧美ⅴa在| 亚洲国产日韩欧美精品在线观看 | 老鸭窝网址在线观看| 18禁美女被吸乳视频| 中亚洲国语对白在线视频| 久久久国产成人免费| 精品熟女少妇八av免费久了| 免费在线观看成人毛片| 别揉我奶头~嗯~啊~动态视频| 很黄的视频免费| 桃色一区二区三区在线观看| 法律面前人人平等表现在哪些方面| 日本在线视频免费播放| 国产精品电影一区二区三区| netflix在线观看网站| 一本大道久久a久久精品| 白带黄色成豆腐渣| 又紧又爽又黄一区二区| 一区二区三区高清视频在线| av天堂在线播放| 久久精品aⅴ一区二区三区四区| 丰满人妻一区二区三区视频av | 免费在线观看黄色视频的| 国产成人一区二区三区免费视频网站| 亚洲成人久久爱视频| 亚洲av电影不卡..在线观看| 亚洲欧美日韩高清在线视频| 香蕉久久夜色| 国产精品1区2区在线观看.| 亚洲精品在线美女| 麻豆久久精品国产亚洲av| 久久久久久九九精品二区国产 | 午夜福利成人在线免费观看| 一级毛片女人18水好多| 亚洲国产日韩欧美精品在线观看 | 99久久久亚洲精品蜜臀av| 国产精品98久久久久久宅男小说| 麻豆久久精品国产亚洲av| 中文字幕人成人乱码亚洲影| 美女黄网站色视频| 亚洲无线在线观看| 一个人观看的视频www高清免费观看 | 99热这里只有精品一区 | 成人av在线播放网站| 最新在线观看一区二区三区| 亚洲av成人精品一区久久| 免费一级毛片在线播放高清视频| 99国产精品一区二区三区| 久久国产精品影院| 久久久久久人人人人人| av国产免费在线观看| 亚洲,欧美精品.| 少妇裸体淫交视频免费看高清 | 亚洲av成人精品一区久久| 熟妇人妻久久中文字幕3abv| 精品国产亚洲在线| 精品久久久久久,| 美女大奶头视频| 久久精品影院6| 亚洲精华国产精华精| 非洲黑人性xxxx精品又粗又长| 99热这里只有精品一区 | 精品日产1卡2卡| 欧美另类亚洲清纯唯美| 婷婷丁香在线五月| 国产探花在线观看一区二区| av免费在线观看网站| 午夜视频精品福利| 黄色女人牲交| 午夜亚洲福利在线播放| 欧美色视频一区免费| 99久久综合精品五月天人人| 亚洲 欧美一区二区三区| 精品少妇一区二区三区视频日本电影| 成人一区二区视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 日本黄大片高清| 亚洲乱码一区二区免费版| 99久久精品热视频| 搡老妇女老女人老熟妇| 国产真实乱freesex| 制服人妻中文乱码| 无限看片的www在线观看| 亚洲全国av大片| 亚洲成人中文字幕在线播放| 国产av不卡久久| 桃红色精品国产亚洲av| 三级毛片av免费| 国产精品久久电影中文字幕| 韩国av一区二区三区四区| 男人舔女人的私密视频| 国内揄拍国产精品人妻在线| 欧美日韩亚洲综合一区二区三区_| 午夜免费激情av| 亚洲精品国产一区二区精华液| 91大片在线观看| 免费在线观看视频国产中文字幕亚洲| 国产一级毛片七仙女欲春2| 老熟妇仑乱视频hdxx| 欧美午夜高清在线| 久久久久久久精品吃奶| 最近在线观看免费完整版| 校园春色视频在线观看| 日日爽夜夜爽网站| 久久久久国产精品人妻aⅴ院| 亚洲av成人不卡在线观看播放网| 国内揄拍国产精品人妻在线| 黄色视频,在线免费观看| 国产av又大| 国产aⅴ精品一区二区三区波| 国产免费男女视频| 免费在线观看影片大全网站| 久久精品aⅴ一区二区三区四区| 免费高清视频大片| 午夜亚洲福利在线播放| 日本一区二区免费在线视频| 久久国产精品人妻蜜桃| 欧美精品啪啪一区二区三区| 国产三级中文精品| 亚洲,欧美精品.| 久久精品国产亚洲av高清一级| 日本一区二区免费在线视频| 怎么达到女性高潮| 国产亚洲av高清不卡| 亚洲成人中文字幕在线播放| 国产精品精品国产色婷婷| 久久久久久亚洲精品国产蜜桃av| 国产成人精品久久二区二区免费| 久久精品91蜜桃| 欧美黄色片欧美黄色片| 国产精品久久视频播放| 美女扒开内裤让男人捅视频| 欧美性长视频在线观看| 久久国产乱子伦精品免费另类| 久久精品影院6| 久久久国产成人免费| 黄色a级毛片大全视频| 身体一侧抽搐| 亚洲av美国av| 久久精品91无色码中文字幕| 久久草成人影院| 久久精品国产综合久久久| 此物有八面人人有两片| 国产精品一区二区三区四区久久| 免费人成视频x8x8入口观看| 国产日本99.免费观看| 99精品在免费线老司机午夜| 久久精品91无色码中文字幕| 亚洲avbb在线观看| 波多野结衣巨乳人妻| 国产高清视频在线播放一区| 亚洲精品色激情综合| 欧美成人一区二区免费高清观看 | 最新在线观看一区二区三区| 一区二区三区国产精品乱码| 日本免费一区二区三区高清不卡| 一进一出好大好爽视频| 69av精品久久久久久| 国内精品久久久久久久电影| 亚洲天堂国产精品一区在线| 久久久久免费精品人妻一区二区| 久久精品国产综合久久久| 日韩免费av在线播放| 91字幕亚洲| 少妇人妻一区二区三区视频| 亚洲精品粉嫩美女一区| bbb黄色大片| 看黄色毛片网站| 伊人久久大香线蕉亚洲五| 国产av不卡久久| 女警被强在线播放| 正在播放国产对白刺激| 亚洲欧美精品综合一区二区三区| 美女免费视频网站| 狠狠狠狠99中文字幕| 精品久久蜜臀av无| 无限看片的www在线观看| 日韩欧美精品v在线| 99国产精品一区二区蜜桃av| 成年免费大片在线观看| 成人永久免费在线观看视频| 久久久国产欧美日韩av| 九九热线精品视视频播放| 非洲黑人性xxxx精品又粗又长| 日韩欧美在线乱码| 神马国产精品三级电影在线观看 | 国产激情久久老熟女| 老熟妇乱子伦视频在线观看| 国产野战对白在线观看| 91成年电影在线观看| 久久这里只有精品19| 黄色毛片三级朝国网站| 全区人妻精品视频| 午夜a级毛片| 啪啪无遮挡十八禁网站| 成人国语在线视频| 在线观看舔阴道视频| 每晚都被弄得嗷嗷叫到高潮| 五月玫瑰六月丁香| 在线观看午夜福利视频| 少妇的丰满在线观看| 成人一区二区视频在线观看| 国产高清视频在线播放一区| 99在线人妻在线中文字幕| 欧美黑人巨大hd| 女人被狂操c到高潮| 日本熟妇午夜| 男人舔女人下体高潮全视频| 最近最新免费中文字幕在线| 午夜久久久久精精品| 黄色成人免费大全| avwww免费| 母亲3免费完整高清在线观看| 免费在线观看亚洲国产| 黑人操中国人逼视频| 亚洲乱码一区二区免费版| 性色av乱码一区二区三区2| 欧美日韩亚洲国产一区二区在线观看| 国产乱人伦免费视频| 国产精品一及| 中亚洲国语对白在线视频| 很黄的视频免费| 国产精品野战在线观看| 久久久久久九九精品二区国产 | 夜夜躁狠狠躁天天躁| 亚洲中文字幕一区二区三区有码在线看 | 国产99久久九九免费精品| 身体一侧抽搐| 两个人免费观看高清视频| 亚洲狠狠婷婷综合久久图片| 美女 人体艺术 gogo| 亚洲午夜理论影院| av中文乱码字幕在线| 国产一区二区在线av高清观看| 亚洲免费av在线视频| 国产亚洲精品av在线| 一级a爱片免费观看的视频| а√天堂www在线а√下载| 久久久久性生活片| 久久 成人 亚洲| 村上凉子中文字幕在线| 国产精品av久久久久免费| 欧美黄色淫秽网站| 婷婷精品国产亚洲av| 长腿黑丝高跟| 别揉我奶头~嗯~啊~动态视频| netflix在线观看网站| 91字幕亚洲| 变态另类丝袜制服| 在线看三级毛片| 91九色精品人成在线观看| 少妇人妻一区二区三区视频| 啪啪无遮挡十八禁网站| 精品久久久久久久久久久久久| 人妻夜夜爽99麻豆av| 久久久久久久午夜电影| www国产在线视频色| 国产精品一区二区三区四区免费观看 | 婷婷精品国产亚洲av在线| 国产伦人伦偷精品视频| 在线播放国产精品三级| 亚洲电影在线观看av| 九九热线精品视视频播放| aaaaa片日本免费| 欧美日韩瑟瑟在线播放| 久久精品aⅴ一区二区三区四区| a在线观看视频网站| 亚洲在线自拍视频| 男人舔奶头视频| 啪啪无遮挡十八禁网站| 在线视频色国产色| 欧美zozozo另类| 丰满人妻熟妇乱又伦精品不卡| 色哟哟哟哟哟哟| 亚洲人成77777在线视频| 国产单亲对白刺激| 欧美精品亚洲一区二区| 一二三四社区在线视频社区8| 99久久精品国产亚洲精品| 欧美不卡视频在线免费观看 | 精品久久久久久成人av| 在线观看www视频免费| av在线天堂中文字幕| 99热这里只有精品一区 | 久久久久久久久中文| 国产av在哪里看| 精品国产乱码久久久久久男人| √禁漫天堂资源中文www| 免费在线观看影片大全网站| 亚洲九九香蕉| 天堂影院成人在线观看| 悠悠久久av| 岛国在线免费视频观看| 一卡2卡三卡四卡精品乱码亚洲| 老汉色∧v一级毛片| 性色av乱码一区二区三区2| 99热这里只有是精品50| 日韩欧美国产一区二区入口| 久久天躁狠狠躁夜夜2o2o| 真人做人爱边吃奶动态| 成人三级黄色视频| 成人手机av| 欧美 亚洲 国产 日韩一| 日韩欧美三级三区| 国产日本99.免费观看| 欧美成人午夜精品| 听说在线观看完整版免费高清| 国内久久婷婷六月综合欲色啪| ponron亚洲| 日本 欧美在线| 淫妇啪啪啪对白视频| 国产野战对白在线观看| 欧美成人性av电影在线观看| 啦啦啦观看免费观看视频高清| 人人妻人人澡欧美一区二区| 国产精品 国内视频| 可以在线观看毛片的网站| 99国产极品粉嫩在线观看| 最近最新免费中文字幕在线| 波多野结衣巨乳人妻| 国产又色又爽无遮挡免费看| 91国产中文字幕| 国产成+人综合+亚洲专区| 久久天堂一区二区三区四区| 一进一出抽搐动态| 亚洲欧美精品综合一区二区三区| АⅤ资源中文在线天堂| 精品午夜福利视频在线观看一区| 亚洲 欧美 日韩 在线 免费| 嫩草影院精品99| 午夜久久久久精精品| 婷婷精品国产亚洲av在线| cao死你这个sao货| 一级毛片精品| 男女下面进入的视频免费午夜| 1024手机看黄色片| 国内精品一区二区在线观看| 变态另类丝袜制服| 老司机福利观看| 我的老师免费观看完整版| 亚洲精品中文字幕一二三四区| 黑人欧美特级aaaaaa片| av中文乱码字幕在线| 国产成人影院久久av| 香蕉久久夜色| 十八禁网站免费在线| 亚洲男人的天堂狠狠| 每晚都被弄得嗷嗷叫到高潮| 一个人观看的视频www高清免费观看 | 免费在线观看亚洲国产| 欧美黑人巨大hd| 亚洲人成网站高清观看| 欧美日本亚洲视频在线播放| tocl精华| 变态另类成人亚洲欧美熟女| 日本免费a在线| 可以在线观看毛片的网站| 婷婷六月久久综合丁香| 夜夜夜夜夜久久久久| 999久久久精品免费观看国产| 最近视频中文字幕2019在线8| 熟女少妇亚洲综合色aaa.| 99国产极品粉嫩在线观看| 91国产中文字幕| 欧美激情久久久久久爽电影| 淫秽高清视频在线观看| 最新美女视频免费是黄的| av国产免费在线观看| 日韩有码中文字幕| 黄色成人免费大全| 欧美日韩亚洲综合一区二区三区_| 亚洲人成网站高清观看| 欧美日韩瑟瑟在线播放| 久久这里只有精品中国| 三级毛片av免费| 亚洲专区国产一区二区| 欧美日韩乱码在线| 久久 成人 亚洲| 十八禁网站免费在线| 亚洲av电影不卡..在线观看| 两个人的视频大全免费| 极品教师在线免费播放| 亚洲人成网站高清观看| 日韩欧美精品v在线| 国产熟女午夜一区二区三区| 亚洲国产精品合色在线| 久久久久久大精品| 亚洲一区高清亚洲精品| 久久久久久九九精品二区国产 | 曰老女人黄片| 制服丝袜大香蕉在线| 精品久久久久久久人妻蜜臀av| 一进一出好大好爽视频| 成人三级做爰电影| 国产高清视频在线播放一区| 午夜福利在线观看吧| 日韩国内少妇激情av| 亚洲av日韩精品久久久久久密| 午夜久久久久精精品| 亚洲精品久久成人aⅴ小说| 免费在线观看成人毛片| 亚洲精品在线美女| cao死你这个sao货| 久久精品亚洲精品国产色婷小说| 久久99热这里只有精品18| 色av中文字幕| 日韩成人在线观看一区二区三区| 欧美在线一区亚洲| 亚洲中文日韩欧美视频| 男女床上黄色一级片免费看| 成人三级做爰电影| 2021天堂中文幕一二区在线观| 欧美乱妇无乱码| 全区人妻精品视频| 日本 欧美在线| 亚洲av五月六月丁香网| 国产精品av久久久久免费| 亚洲欧美日韩高清在线视频| 两个人免费观看高清视频| 成人亚洲精品av一区二区| 无遮挡黄片免费观看| 亚洲欧美激情综合另类| 精品电影一区二区在线| 真人一进一出gif抽搐免费| 精品少妇一区二区三区视频日本电影| 精品免费久久久久久久清纯| 搡老熟女国产l中国老女人| 女生性感内裤真人,穿戴方法视频| 黄色丝袜av网址大全| avwww免费| 美女扒开内裤让男人捅视频| 亚洲真实伦在线观看| 日韩有码中文字幕| 一区二区三区高清视频在线| 午夜福利视频1000在线观看| 婷婷精品国产亚洲av| 久久亚洲真实| 少妇熟女aⅴ在线视频| 欧美性猛交黑人性爽| 亚洲 欧美一区二区三区| 在线免费观看的www视频| 精品久久久久久久久久久久久| 1024香蕉在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产爱豆传媒在线观看 | 麻豆成人av在线观看| 国产成人影院久久av| 亚洲成人免费电影在线观看| 午夜福利在线观看吧| 亚洲欧美精品综合一区二区三区| 欧美成人一区二区免费高清观看 | 麻豆成人av在线观看| 亚洲av熟女| 淫妇啪啪啪对白视频| 亚洲精品美女久久av网站| 看免费av毛片| 久久精品影院6| 免费一级毛片在线播放高清视频| 午夜福利免费观看在线| 国内精品久久久久精免费| 美女高潮喷水抽搐中文字幕| 精品一区二区三区av网在线观看| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 久久这里只有精品中国| 伦理电影免费视频| 国产精品久久久人人做人人爽| 日韩欧美在线二视频| 亚洲精品在线观看二区| 久久久国产成人免费| 亚洲专区中文字幕在线| 好男人在线观看高清免费视频| 我的老师免费观看完整版| 成人手机av| 麻豆久久精品国产亚洲av| 99精品久久久久人妻精品| 国产精品综合久久久久久久免费| 最近在线观看免费完整版| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久久久人妻精品电影| 精品福利观看| 老鸭窝网址在线观看| 激情在线观看视频在线高清| 一级片免费观看大全| 日本成人三级电影网站| 嫁个100分男人电影在线观看| 成人欧美大片| 国产av不卡久久| 亚洲午夜理论影院| 欧美日韩乱码在线| 欧美乱色亚洲激情| 国产成年人精品一区二区| 中文字幕久久专区| 老汉色av国产亚洲站长工具| 可以在线观看毛片的网站| 国产精品爽爽va在线观看网站| 日日摸夜夜添夜夜添小说| 91成年电影在线观看| 好看av亚洲va欧美ⅴa在| 两个人视频免费观看高清| 欧美不卡视频在线免费观看 | 亚洲国产日韩欧美精品在线观看 | 麻豆成人av在线观看| 俺也久久电影网| 夜夜看夜夜爽夜夜摸| 久久久久国内视频| 身体一侧抽搐| 欧美人与性动交α欧美精品济南到| 一级黄色大片毛片| 两个人视频免费观看高清| 欧美一级a爱片免费观看看 | 女生性感内裤真人,穿戴方法视频| 99在线视频只有这里精品首页| 免费一级毛片在线播放高清视频| 亚洲av成人不卡在线观看播放网| 极品教师在线免费播放| 国产人伦9x9x在线观看| 久久婷婷成人综合色麻豆| 免费搜索国产男女视频| 岛国在线观看网站| 国产成人精品久久二区二区91| 日韩精品中文字幕看吧| 最近最新中文字幕大全免费视频| 亚洲色图av天堂| av视频在线观看入口| 韩国av一区二区三区四区| 变态另类丝袜制服| 在线观看舔阴道视频| 国内毛片毛片毛片毛片毛片| 不卡av一区二区三区| 免费av毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 男女床上黄色一级片免费看| 国产99久久九九免费精品| 国产精品精品国产色婷婷| 久久久久国产精品人妻aⅴ院| www.www免费av| 夜夜看夜夜爽夜夜摸|