亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        壓電俘能器渦激振動俘能的建模與實驗研究

        2016-12-21 02:07:48宋汝君單小彪李晉哲謝濤
        西安交通大學(xué)學(xué)報 2016年2期
        關(guān)鍵詞:俘能器渦激輸出功率

        宋汝君,單小彪,李晉哲,謝濤

        (1.哈爾濱工業(yè)大學(xué)機電工程學(xué)院, 150001, 哈爾濱;2.東北林業(yè)大學(xué)機電工程學(xué)院, 150040, 哈爾濱)

        ?

        壓電俘能器渦激振動俘能的建模與實驗研究

        宋汝君1,單小彪1,李晉哲2,謝濤1

        (1.哈爾濱工業(yè)大學(xué)機電工程學(xué)院, 150001, 哈爾濱;2.東北林業(yè)大學(xué)機電工程學(xué)院, 150040, 哈爾濱)

        針對微機電系統(tǒng)和傳感器等低能耗電子產(chǎn)品的持續(xù)供能問題,提出了一種渦激振動式壓電俘能器。該俘能器由壓電懸臂梁和末端圓柱體組成,結(jié)構(gòu)簡單,可在較低水流流速下產(chǎn)生渦激共振,得到較大的能量輸出。通過數(shù)學(xué)建模和實驗測試的方法,研究了水流速度和外接電阻對壓電俘能器振動和俘能的影響規(guī)律。實驗結(jié)果表明:壓電俘能器的振動頻率隨流速的增大而增大,振動幅值在渦激共振時最大,輸出功率受流速和外接電阻兩者影響,較小外接電阻適合較高流速,較大電阻適合較低流速,壓電俘能器在渦激共振處可獲得最大的能量輸出,當(dāng)外接電阻為0.5 MΩ、流速為0.41 m/s時,實驗測試得到了8.3 μW的最大輸出功率。數(shù)值分析結(jié)果與實驗測試結(jié)果吻合較好,驗證了數(shù)學(xué)模型的正確性。

        水流;渦激振動;壓電俘能器;俘能

        隨著集成電路、MEMS和微傳感器等低能耗電子產(chǎn)品的發(fā)展,電池供能的方式存在質(zhì)量大、體積大、存儲有限、需要定期更換等諸多弊端,因此亟需一種可持續(xù)的能源供給方式。同時,環(huán)境振動能分布廣泛,有水動能、風(fēng)動能、潮汐能、機械振動能等,其中水動能是一種清潔的可再生能源。壓電俘能器是一種結(jié)構(gòu)簡單、環(huán)境友好、使用壽命長久和無電磁干擾的振動能轉(zhuǎn)換裝置,可轉(zhuǎn)換水動能為電能,能夠滿足低能耗電子產(chǎn)品的供能需求[1]。

        既往的研究中,壓電俘能器的研究重點集中在俘能器本身結(jié)構(gòu)和俘能性能上,忽略了壓電俘能器與環(huán)境振動的相互耦合影響[2-3]。水流激勵壓電俘能同時考慮了環(huán)境振動與俘能器本身的耦合影響,同時涉及流體場、固體場和電場,是一個多物理場相互耦合的能量轉(zhuǎn)換技術(shù),具有較高的學(xué)術(shù)研究價值,目前已經(jīng)得到了國內(nèi)外研究者的重視。本課題組也從俘能器本身構(gòu)型的研究[4-6]逐漸轉(zhuǎn)向壓電俘能器耦合流體環(huán)境的研究上[7-8]。壓電俘能器的流致振動俘能技術(shù)按照流激振動形式可分為尾流馳振式[9-10]、渦激振動式[11-12]、拍振式[13-14]和馳振式[15-17]等4種。尾流馳振式壓電俘能器可工作于較低流速下,具有較寬的振動頻帶,但激勵響應(yīng)較小。渦激振動式壓電俘能器多以尾端圓柱與壓電振子構(gòu)成,可在較低流速下(如低速水流)產(chǎn)生渦激共振使壓電俘能器的振動幅值和能量輸出最大;且由于“鎖定”現(xiàn)象的存在,渦激共振區(qū)域增大,工作頻帶增寬。馳振式和拍振式壓電俘能器在流體激勵下具有相同的振動特征,即振幅發(fā)散式增幅振動,它存在一個臨界速度,當(dāng)來流流速低于臨界速度時,壓電俘能器靜止,當(dāng)速度大于臨界速度時,其振動幅值隨速度的增大而逐漸增大,因此拍振式和馳振式壓電俘能器更適合于較高流速環(huán)境下(如風(fēng)流)。

        綜上發(fā)現(xiàn),渦激振動式壓電俘能器適合于水流激勵下振動俘能,然而現(xiàn)有的研究主要圍繞空氣流進(jìn)行[11-12]。相比于空氣,水流一般流速較小,但具有較大的密度和黏度。因此,流體激勵力、激勵頻率、流體附加質(zhì)量和附加黏度對壓電俘能器的振動俘能的影響截然不同。鑒于此,本文提出了一種由壓電懸臂梁和末端圓柱體組成的渦激振動式壓電俘能器,并建立了壓電俘能器振動響應(yīng)和能量輸出的數(shù)學(xué)模型,分別分析了水流流速和外接電阻等參數(shù)對俘能器振動和發(fā)電能力的影響規(guī)律,并設(shè)計壓電俘能器樣機,通過實驗研究了壓電俘能器在水流環(huán)境下的功率輸出。

        1 結(jié)構(gòu)及數(shù)學(xué)模型

        圓柱繞流是自然界較為常見的現(xiàn)象,當(dāng)流體流過任何非流線型物體時,都會在物體兩側(cè)交替產(chǎn)生脫落的旋渦,交替脫落的旋渦會在垂直于來流方向產(chǎn)生周期性變化作用力。周期性作用力可激勵彈性固支的壓電俘能器振動并發(fā)電,圖1給出了渦激振動壓電俘能系統(tǒng)示意圖和尺寸圖。

        (a)渦激振動壓電俘能系統(tǒng)三維示意圖

        (b)壓電俘能器尺寸圖圖1 渦激振動壓電俘能系統(tǒng)示意圖和尺寸圖

        由圖1可知,壓電俘能器由壓電懸臂梁自由端與圓柱連接組成。壓電梁長度為L,寬度為b,分別由厚度為hs的支撐層和厚度為hp的壓電層組成。圓柱體直徑為D,長度為Lc,外接電阻為R。

        假設(shè)壓電懸臂梁為線性變形并僅發(fā)生在z方向上,壓電懸臂梁運動的離散控制方程表示為

        dδ(x-L)/dx]=Ftδ(x-L)

        (1)

        根據(jù)Facchinetti等的研究,圓柱渦激振動的渦激力[18]可表示為

        Ft=0.5CLρDU2Lc-

        0.5CDρDULc(?w(x,t)/?t|x=L+

        0.5D?2w(x,t)/?x?t|x=L)

        (2)

        式中:CL、CD分別是脈動升力和阻力系數(shù),CL=Clq(t)/2、CD=2,其中Cl=0.3為靜升力系數(shù),q(t)用來描述圓柱渦街脫落的尾跡運動,其運動方程為

        (A/D)(?w2(x,t)/?t2|x=L+

        (0.5D?3w(x,t)/?x?t2|x=L)

        (3)

        根據(jù)Facchinetti的研究[18],λ=1,ε=0.3,A=12,ωf是渦脫頻率,ωf=2πStU/D,其中St為斯特勞哈爾數(shù),在雷諾數(shù)為300

        根據(jù)圖1b,俘能器可等效為圖2所示的電路。

        圖2 俘能系統(tǒng)的等效耦合電路

        根據(jù)圖2,等效耦合電路方程可寫成

        (4)

        式中:Cp為壓電俘能器的等效電容。

        俘能系統(tǒng)的輸出功率表示為

        (5)

        式中:T為輸出電壓的交變周期。

        通過聯(lián)立式(1)~式(5),即可求得俘能系統(tǒng)的輸出功率。為了求解耦合方程組式(1)~式(4),根據(jù)分離變量法,將位移函數(shù)ω(x,t)表示為

        φr(x)ηr(t)

        (6)

        式中:φr(x)為r階的主振型函數(shù);ηr(t)為廣義坐標(biāo)。

        壓電俘能器的邊界條件表示為

        w(0,t)=0, ?w(x,t)/?x|x=0=0

        (7)

        (8)

        式中:M=M0+Mf,M0為圓柱質(zhì)量,Mf為流體對圓柱和壓電梁的附加質(zhì)量;It為相對于中心的轉(zhuǎn)動慣量,且It=M(D/2)2/2。流體的附加質(zhì)量[18-19]可表示為

        Mf=ρfπD2Lc/4+(33/140)πρfb2(L/4)

        (9)

        懸臂梁的各階主振型[2]為

        φr(x)=Ar[cos(λrx/L)-cosh(λrx/L)+

        γr(sin(λrx/L)-sinh(λrx/L)]

        (10)

        式中

        (11)

        Ar為振型系數(shù),滿足各階主振型正則化,具體見文獻(xiàn)[2]。

        ωr是系統(tǒng)的無阻尼固有頻率,可表示為

        (12)

        式中:λr由式(13)獲得

        cosλrcoshλr+(λrM/mL)(cosλrsinhλr-

        sinλrcoshλr)=0

        (13)

        將式(2)代入式(1),并將式(6)代入式(1)、式(3)和式(4),方程兩邊同時乘以振型函數(shù)φr(x)并沿整個懸臂壓電梁長度L積分,然后利用正則振型的正交性簡化方程,最后得到對應(yīng)各階振型的耦合振動方程組

        (14)

        (15)

        (16)

        式中:χr=θφ′r(L);ζr=csωr/(2K)+ca/(2mωr)。

        為了求解式(14)~式(16),引入狀態(tài)變量

        (17)

        因此,可將式(14)~式(16)寫成

        (18)

        式中:

        (19)

        χrX2n+3+φr(L)(ClρDU2LcX2n+1)/2

        (20)

        2 非線性數(shù)值計算分析

        壓電片選用型號為MFC-M8514-P2(產(chǎn)自Smart Material Corp?)的壓電纖維復(fù)合材料,其壓電片表面有良好的絕緣作用,可置于水流中進(jìn)行渦激振動俘能工作。壓電俘能器的支撐層選用純鋁材質(zhì),圓柱采用尼龍材料制作,并通過自由振動衰減法測量壓電俘能系統(tǒng)的一階阻尼比為ζ1=0.086。

        表1 壓電俘能器的結(jié)構(gòu)與材料參數(shù)

        根據(jù)表1所給出的參數(shù),通過MATLAB?軟件中的ode45求解器求解式(18),根據(jù)數(shù)值計算結(jié)果分析了水流流速和外接電阻對壓電俘能器的振動和俘能的影響規(guī)律,具體如圖3~圖6所示。圖3給出了在外接電阻為0.5 MΩ,水流速度分別為0.3、0.5和0.7 m/s時,壓電俘能器的振動和輸出電壓情況。由圖3可知,壓電俘能器在水流的激勵下,從靜止?fàn)顟B(tài)經(jīng)過幾個周期的調(diào)整后變?yōu)榉€(wěn)定的周期性振動;壓電俘能器在不同的流速下振動幅值不同,流速為0.5 m/s時的振動幅值大于流速為0.3和0.7 m/s時的振動幅值,流速為0.7 m/s時的振動頻率大于流速為0.3和0.5 m/s時的振動頻率,振動頻率隨流速的增大而增大,如圖 3a和圖3b所示;由圖3c和圖3d可知,輸出電壓與壓電俘能器的振動情況有相同的變化趨勢,電壓變化頻率等于振動頻率。

        圖4和圖5給出了在外接電阻為0.05、0.1、0.5和1 MΩ時,壓電俘能器振動幅值和輸出功率隨流速的變化情況。由圖4可知:外接電阻對壓電俘能器的振動幅值影響較小,流速是振動幅值的主要影響因素;振動幅值隨流速的增大先增大后減小,在流速為0.41 m/s處達(dá)到最大值,此時處于渦激共振區(qū)域內(nèi);振動幅值達(dá)到最大值后,隨著流速的增大,振動幅值緩慢降低,頻帶較寬,有利于壓電俘能器的振動俘能。由圖5可知,不同外接電阻下,壓電俘能器的輸出功率隨來流速度的變化趨勢不同:當(dāng)外接電阻較小時(例如0.05和0.1 MΩ),輸出功率隨流速的增大先增大,然后趨于平穩(wěn)變化;當(dāng)外接電阻較大時(例如0.5和1 MΩ),輸出功率隨流速的增大先增大后減小,分別在0.43 m/s處達(dá)到最大值。由此可知:較大的外接電阻更適合較低流速環(huán)境,而較小的外接電阻則更適合較高流速環(huán)境。

        圖4 不同外接電阻下壓電俘能器尾端振動幅值隨流速的變化曲線

        圖5 不同外接電阻下壓電俘能器輸出功率隨流速的變化曲線

        為了進(jìn)一步闡明外接電阻、來流速度和輸出功率三者之間的關(guān)系,圖6給出了在不同的流速情況下輸出功率隨外接電阻的變化曲線圖。由圖6可知,在不同的流速下,輸出功率都隨外接電阻的增大而先增大后減小,存在一個電阻值使壓電俘能器的輸出功率最大,稱之為最優(yōu)電阻值,且最優(yōu)電阻值隨流速的增大而逐漸減小,這是因為隨著流速的增大,壓電俘能器的振動頻率增大,相應(yīng)的其最優(yōu)電阻值逐漸減小,也是在較低流速環(huán)境中采用較大的外接電阻,在較高流速情況下采用較小的外接電阻可獲得更多的能量輸出的原因。

        圖6 不同流速下輸出功率隨外接電阻的變化曲線

        3 實驗研究

        根據(jù)表1中的壓電俘能器結(jié)構(gòu)參數(shù)設(shè)計并制作了實驗樣機,搭建的渦激振動俘能系統(tǒng)實驗平臺如圖7所示。壓電梁懸臂梁在來流前方固定,為了降低夾具對來流的干擾,采用流線型夾具加持。實驗臺采用水循環(huán)系統(tǒng),主要由水槽、離心水泵、進(jìn)水管和出水管組成。實驗段采用明渠水槽,長度為1 500 mm,寬度為400 mm。實驗中,水流深度保持恒定,即水流截面積保持不變,通過調(diào)節(jié)實驗段的傾斜度和流體流量綜合控制水流流速。實驗段水槽的傾斜度由千斤頂調(diào)節(jié);實驗段水流流量通過進(jìn)水管和出水管組合控制??烧{(diào)電阻器調(diào)節(jié)外接電阻值大小,并通過示波器顯示壓電俘能器輸出電壓值。

        圖7 渦激振動俘能實驗系統(tǒng)圖

        利用圖7所示的渦激振動俘能實驗系統(tǒng),測試了壓電俘能器在流速0.1~0.6 m/s范圍內(nèi),不同外接電阻下的輸出功率。圖8給出了實驗測得壓電俘能器在外接電阻為0.05、0.1、0.5和1 MΩ時,輸出功率隨流速變化的實驗結(jié)果。

        圖8 不同外接電阻下實驗測得壓電俘能器輸出功率隨流速的變化曲線

        由圖8可見,當(dāng)外接電阻為0.05和0.1 MΩ時,輸出功率隨流速的增大先增大,在流速為0.4 m/s之后隨流速的增大而平穩(wěn)變化;當(dāng)外接電阻為0.5和1 MΩ時,輸出功率隨流速的增大先增大后減小,在渦激共振點0.41 m/s處達(dá)到最大值,分別為8.3和6.1 μW。對比圖5發(fā)現(xiàn),數(shù)值分析結(jié)果與實驗測試結(jié)果吻合較好,驗證了渦激振動壓電俘能系統(tǒng)數(shù)學(xué)建模的正確性。

        4 結(jié) 論

        本文提出了一種基于水流環(huán)境下的渦激振動式壓電俘能器,并通過數(shù)學(xué)建模和實驗測試,研究了水流流速和外接電阻對壓電俘能器的振動特性和輸出功率的影響規(guī)律。研究表明:壓電俘能器在渦激共振時振動幅值最大,振動頻率隨流速的增大而增大;輸出功率受外接電阻和流速兩者的綜合影響,較低的流速下采用較大的外接電阻和較高的流速下采用較小的電阻可獲得更大的功率輸出;壓電俘能器在渦激共振時能量輸出最大,當(dāng)外接電阻為0.5 MΩ、流速為0.41 m/s時,實驗測試得到了8.3 μW的最大輸出功率。數(shù)學(xué)模型數(shù)值計算結(jié)果和實驗測試結(jié)果吻合較好,驗證了數(shù)學(xué)建模的正確性。

        [1] 袁江波, 謝濤, 單小彪, 等. 壓電俘能技術(shù)研究現(xiàn)狀綜述 [J]. 振動與沖擊, 2009, 28(10): 36-42. YUAN Jiangbo, XIE Tao, SHAN Xiaobiao, et al. A review of current situation for piezoelectric energy harvesting [J]. Journal of Vibration and Shock, 2009, 28(10): 36-42.

        [2] ERTURK A, INMAN D J. An experimentally validated bimorph cantilever model for piezoelectric energy harvesting from base excitations [J]. Smart Materials and Structures, 2009, 18(2): 025009.

        [3] ERTURK A, INMAN D J. A distributed parameter electromechanical model for cantilevered piezoelectric energy harvesters [J]. Journal of Vibration and Acoustics, 2008, 130(4): 041002.

        [4] 謝濤, 袁江波, 單小彪, 等. 多懸臂梁壓電振子頻率分析及發(fā)電實驗研究 [J]. 西安交通大學(xué)學(xué)報, 2010, 44(2): 98-101. XIE Tao, YUAN Jiangbo, SHAN Xiaobiao, et al. Frequency analysis and electricity generated by multiple piezoelectric cantilevers in energy harvesting [J]. Journal of Xi’an Jiaotong University, 2010, 44(2): 98-101.

        [5] 王紅艷, 謝濤, 單小彪, 等. 壓電懸臂梁俘能器多模態(tài)俘能效果研究 [J]. 西安交通大學(xué)學(xué)報, 2010, 44(6): 114-119. WANG Hongyan, XIE Tao, SHAN Xiaobiao, et al. Multi-mode vibration with piezoelectric cantilever for energy harvesting [J]. Journal of Xi’an Jiaotong University, 2010, 44(6): 114-119.

        [6] WANG Hongyan, TANG Lihua, SHAN Xiaobiao, et al. Modeling and performance evaluation of a piezoelectric energy harvester with segmented electrodes [J]. Smart Structures and Systems, 2014, 14(2): 247-266.

        [7] SONG Rujun, SHAN Xiaobiao, LV Fengchi, et al. A study of vortex-induced energy harvesting from water using Pzt piezoelectric cantilever with cylindrical extension [J]. Ceramics International, 2015, 41(S1): S768-S773.

        [8] SHAN Xiaobiao, SONG Rujun, LIU Bo, et al. Novel energy harvesting: a macro fiber composite piezoelectric energy harvester in the water vortex [J]. Ceramics International, 2015, 41(S1): S763-S767.

        [9] AKAYDIN H D, ELVIN N, ANDREOPOULOS Y. Energy harvesting from highly unsteady fluid flows using piezoelectric materials [J]. Journal of Intelligent Material Systems and Structures, 2010, 21(13): 1263-1278.

        [10]AKAYDIN H D, ELVIN N, ANDREOPOULOS Y. Wake of a cylinder: a paradigm for energy harvesting with piezoelectric materials [J]. Experiments in Fluids, 2010, 49(1): 291-304.

        [11]DAI H L, ABDELKEFI A, WANG L. Piezoelectric energy harvesting from concurrent vortex-induced vibrations and base excitations [J]. Nonlinear Dynamics, 2014, 77(3): 967-981.

        [12]GAO Xiaotong, SHIH W H, SHIH W Y. Flow energy harvesting using piezoelectric cantilever with cylindrical extension [J]. IEEE Transactions on Industrial Electronics, 2013, 60(3): 1116-1118.

        [13]MICHELIN S, DOARE O. Energy harvesting efficiency of piezoelectric flags in axial flows [J]. Journal of Fluid Mechanics, 2013, 714(1): 489-504.

        [15]ABDELKEFI A, YAN Z, HAJJ M R. Modeling and nonlinear analysis of piezoelectric energy harvesting from transverse galloping [J]. Smart Materials and Structures, 2013, 22(2): 025016.

        [16]SIROHI J, MAHADIK R. Piezoelectric wind energy harvester for low-power sensors [J]. Journal of Intelligent Material Systems and Structures, 2011, 22(18): 2215-2228.

        [17]ZHAO Liya, TANG Lihua, YANG Yiaowen. Enhanced piezoelectric galloping energy harvesting using 2 degree-of-freedom cut-out cantilever with magnetic interaction [J]. Japanese Journal of Applied Physics, 2014, 53(6): 060302.

        [18]FACCHINETTI M L, DE LANGRE E, BIOLLEY F. Coupling of structure and wake oscillators in vortex-induced vibrations [J]. Journal of Fluids and Structures, 2004, 19(2): 123-140.

        [19]BASAK S, RAMAN A, GARIMELLA S V. Hydrodynamic loading of microcantilevers vibrating in viscous fluids [J]. Journal of Applied Physics, 2006, 99(11): 114906.

        (編輯 武紅江)

        Modeling and Experimental Study of Piezoelectric Energy Harvester under Vortex-Induced Vibration

        SONG Rujun1,SHAN Xiaobiao1,LI Jinzhe2,XIE Tao1

        (1. School of Mechatronics Engineering, Harbin Institute of Technology, Harbin 150001, China;2. College of Mechanical and Electrical Engineering, Northeast Forestry University, Harbin 150040, China)

        A vortex-induced piezoelectric energy harvester (PEH) is proposed to meet the requirement of continuous energy supply for low energy-consumed electronic devices such as MEMS and sensors. The PEH is composed of a piezoelectric cantilever beam and a cylinder and has a simple mechanical structure. The PEH generates vortex-induced resonance at relatively low water velocities and obtains larger output powers. Effects of both the water velocity and the resistance on hydrodynamic response and energy harvesting ability of the PEH are investigated through mathematical modeling and experimental study. It is found that the vibration frequency increases as the velocity increases and the maximum vibration amplitude is found when the vortex-induced resonance appears. The Output power depends on both the water velocity and the resistance. More output power can be achieved through configurations of a smaller resistance with a higher velocity or a larger resistance with a lower velocity. The maximum output power is obtained with vortex-induced resonance. The maximum output power of 8.3 μW is obtained in experimental study with 0.5 MΩ and 0.41 m/s. Furthermore, the results of numerical study are consistent with those of the experimental study well, and the validity of the mathematical model is verified.

        water flow; vortex-induced vibration; piezoelectric energy harvester; energy harvesting

        2015-08-04。

        宋汝君(1988—),男,博士生;謝濤(通信作者),男,教授。 基金項目:國家自然科學(xué)基金資助項目(50875057);中央高?;究蒲袠I(yè)務(wù)費專項資金資助項目(HIT.NSRIF.2014059;HIT.KISTP.201412)。

        時間:2016-01-07

        10.7652/xjtuxb201602010

        TN 384;TP211

        A

        0253-987X(2016)02-0055-06

        網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20160107.1232.004.html

        猜你喜歡
        俘能器渦激輸出功率
        汽車懸架的振動俘能與汽車平順性仿真研究
        計算機仿真(2023年7期)2023-09-04 14:36:06
        不同間距比下串聯(lián)圓柱渦激振動數(shù)值模擬研究
        壓電疊堆—蜂鳴片復(fù)合俘能器設(shè)計*
        渦激振動發(fā)電裝置及其關(guān)鍵技術(shù)
        一種質(zhì)量塊-彈簧自參數(shù)共振壓電俘能器研究
        壓電與聲光(2020年5期)2020-10-28 05:05:26
        一種壓電式微型俘能器
        盤球立管結(jié)構(gòu)抑制渦激振動的數(shù)值分析方法研究
        電子制作(2018年14期)2018-08-21 01:38:42
        柔性圓管在渦激振動下的模態(tài)響應(yīng)分析
        適用于智能電網(wǎng)的任意波形輸出功率源
        基于雙層BP神經(jīng)網(wǎng)絡(luò)的光伏電站輸出功率預(yù)測
        男女男在线精品免费观看| 欧洲极品少妇| 北条麻妃毛片在线视频| 无码人妻丰满熟妇区免费| h视频在线免费观看视频| 亚洲av无码一区东京热久久| 8ⅹ8x擦拨擦拨成人免费视频| 国产欧美日韩在线观看一区二区三区| 一二区视频免费在线观看| 国产偷国产偷亚洲高清视频| 亚洲一区 日韩精品 中文字幕| 日韩欧美区| 国产一区二区三区av香蕉| 亚洲av区,一区二区三区色婷婷| 三年片大全在线观看免费观看大全| 人人妻人人澡人人爽精品欧美 | 中文字幕人乱码中文字幕乱码在线| 伊人情人色综合网站| 久久99热久久99精品| 日本免费一区二区三区在线看| 亚洲一区二区三区毛片| 插插射啊爱视频日a级| 四虎国产精品免费久久| 欧美日韩精品一区二区三区高清视频| 一区=区三区国产视频| 国产av无码专区亚洲av男同| 国产成人午夜精华液| 国产高清在线精品一区αpp| 日本一级三级在线观看| 一本加勒比hezyo无码专区| 欧美性videos高清精品| 超级碰碰人妻中文字幕| 97青草超碰久久国内精品91| 国模无码一区二区三区不卡| 香蕉成人啪国产精品视频综合网| 亚洲成人色黄网站久久| 亚洲精品人成中文毛片| 99久久精品费精品国产一区二区| 亚洲成人av一区二区三区| 久久成人精品国产免费网站| 人人鲁人人莫人人爱精品|