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

        ?

        空間桁架結(jié)構(gòu)的優(yōu)化配置及振動(dòng)控制研究

        2016-11-17 10:13:03許建國張靜靜董龍雷
        計(jì)算機(jī)測量與控制 2016年3期
        關(guān)鍵詞:作動(dòng)器桿件壓電

        許建國,邵 康,張靜靜,董龍雷

        (1.南京東大能源工程設(shè)計(jì)院,南京 211100; 2.西安交通大學(xué) 航天航空學(xué)院,西安 710049)

        ?

        空間桁架結(jié)構(gòu)的優(yōu)化配置及振動(dòng)控制研究

        許建國1,邵 康2,張靜靜2,董龍雷2

        (1.南京東大能源工程設(shè)計(jì)院,南京 211100; 2.西安交通大學(xué) 航天航空學(xué)院,西安 710049)

        為抑制空間柔性桁架結(jié)構(gòu)的低頻振動(dòng),采用壓電桿件進(jìn)行優(yōu)化配置實(shí)現(xiàn)桁架結(jié)構(gòu)的振動(dòng)主動(dòng)控制;建立了空間桁架結(jié)構(gòu)主動(dòng)壓電桿件的機(jī)電耦合模型,利用ANSYS前處理功能編制了壓電桁架的機(jī)電耦合有限元程序;將可控性度量指標(biāo)與逐步消減法相結(jié)合,實(shí)現(xiàn)了空間桁架結(jié)構(gòu)主動(dòng)桿件的優(yōu)化配置;對結(jié)構(gòu)進(jìn)行初始位移擾動(dòng)、正弦激勵(lì)以及隨機(jī)激勵(lì),并采用最優(yōu)模態(tài)控制算法進(jìn)行振動(dòng)抑制仿真分析,對上述方法進(jìn)行驗(yàn)證且建立振動(dòng)控制評價(jià)指標(biāo)進(jìn)行評價(jià);結(jié)果表明將可控性度量指標(biāo)與逐步消減法相結(jié)合的方法可有效抑制空間柔性桁架結(jié)構(gòu)的振動(dòng)。

        空間桁架;振動(dòng)控制;優(yōu)化配置;模態(tài)控制

        0 引言

        近年來,大型化、低剛度、柔性化是各衛(wèi)星天線結(jié)構(gòu)發(fā)展的重要趨勢,然而該類柔性結(jié)構(gòu)具有模態(tài)阻尼小等特性,當(dāng)在太空運(yùn)行受到外界激勵(lì)時(shí),若不采取措施對其振動(dòng)進(jìn)行抑制,則振動(dòng)將會(huì)持續(xù)很久,造成定向精度、姿態(tài)穩(wěn)定等問題,嚴(yán)重時(shí)可能導(dǎo)致結(jié)構(gòu)受到破壞,降低結(jié)構(gòu)使用壽命[1-2]。

        由于壓電陶瓷材料具有剛度高、頻帶寬、高轉(zhuǎn)換率、質(zhì)量輕等特點(diǎn),在桁架結(jié)構(gòu)中常被用于制作傳感器和作動(dòng)器。鄭凱等人設(shè)計(jì)了采用壓電堆作動(dòng)器的壓電主動(dòng)構(gòu)件[3];曹玉巖等人根據(jù)有限元理論和Hamilton原理建立了平面智能桁架結(jié)構(gòu)動(dòng)力學(xué)模型,并根據(jù)線性二次型最優(yōu)控制理論建立了最優(yōu)振動(dòng)控制模型[4];陳文英等人基于Lyapunov綜合法設(shè)計(jì)了一種穩(wěn)定的直接型自適應(yīng)模糊主動(dòng)振動(dòng)控制器并在一平面16桿智能桁架結(jié)構(gòu)進(jìn)行了仿真驗(yàn)證[5];Marinaki 等人研究了應(yīng)用PSO 算法優(yōu)化的模糊控制器進(jìn)行梁的振動(dòng)控制[6];Sun 研究了撓性航天器的模糊預(yù)測振動(dòng)[7];Zhu 等人研究了采用滑模模糊進(jìn)行撓性航天器的振動(dòng)和姿態(tài)機(jī)動(dòng)控制[8];所有這些研究都取得了較好的控制結(jié)果,然而關(guān)于空間桁架結(jié)構(gòu)振動(dòng)的研究很少;Armaghan等人建立了桁架結(jié)構(gòu)三維梁模型[9];李東旭等人采用遺傳算法得到了作動(dòng)器和傳感器在大型智能桁架系統(tǒng)中的最優(yōu)位置并利用T字型桁架結(jié)構(gòu)進(jìn)行了驗(yàn)證性試驗(yàn),然而遺傳算法中各參數(shù)需主觀確定且計(jì)算速度較慢[10]。

        針對作動(dòng)器、傳感器位置優(yōu)化問題,目前已經(jīng)提出了多種作動(dòng)器傳感器優(yōu)化配置準(zhǔn)則,主要有可控度可觀測準(zhǔn)則、系統(tǒng)能量準(zhǔn)則、系統(tǒng)響應(yīng)準(zhǔn)則、可靠性準(zhǔn)則、控制溢出觀測溢出準(zhǔn)則等??煽囟葴?zhǔn)則具有明確的物理意義,不需要與控制算法聯(lián)立,求解簡單有效,獲得廣泛的應(yīng)用。劉福強(qiáng)[11]利用逐步消減法對作動(dòng)器傳感器配置做了優(yōu)化,提高了求解效率,Sylvaine L從基于能量的角度出發(fā),給出了3種不同的準(zhǔn)則來配置作動(dòng)器輸入能量最小準(zhǔn)則、結(jié)構(gòu)總能量最小準(zhǔn)則,可控性格蘭姆矩陣特征值最大準(zhǔn)則,并且指出,在阻尼較小的情況下,它們是等價(jià)的,并給出了可控度綜合指標(biāo)[12]。

        本文以空間柔性桁架結(jié)構(gòu)為研究對象,首先建立了空間桁架結(jié)構(gòu)主動(dòng)壓電桿件的機(jī)電耦合模型,并利用ANSYS、Matlab編制了空間桁架結(jié)構(gòu)動(dòng)力學(xué)分析程序;其次將綜合可控度指標(biāo)與逐步消減法結(jié)合實(shí)現(xiàn)主動(dòng)壓電桿件的優(yōu)化配置;最后采用最優(yōu)模態(tài)控制算法進(jìn)行了空間桁架結(jié)構(gòu)的控制仿真驗(yàn)證,并進(jìn)行振動(dòng)控制效果評價(jià)。

        1 空間桁架動(dòng)力學(xué)模型

        1.1 壓電桿件作動(dòng)理論

        桁架結(jié)構(gòu)主動(dòng)桿件采用多層式壓電陶瓷驅(qū)動(dòng)器,如圖1所示,φ為壓電薄片兩端電勢差,t為壓電薄片厚度,L為壓電桿件長度,N1、N2分別為桿件兩端軸向力。

        圖1 壓電桿件

        若只考慮桿軸方向的壓電效應(yīng),圖1所示壓電桿件的力平衡方程和電荷平衡方程為

        (1)

        其中:e33為軸向壓電系數(shù),c33為軸向彈性系數(shù),ε33為軸向介電系數(shù)。

        每個(gè)壓電薄片均可視為平行板電容器,若沒有泄漏電流,電荷量{Q}與施加于壓電薄片外電壓{V}之間關(guān)系為

        (2)

        采用有限元方法[14],對公式(2)進(jìn)行坐標(biāo)轉(zhuǎn)換,用Guran方法消去電勢差,得到壓電作動(dòng)器的有限元?jiǎng)恿W(xué)方程為

        (3)

        壓電單元端部需要連接普通桿件或預(yù)壓彈簧,此時(shí)壓電主動(dòng)桿件需考慮耦合剛度項(xiàng)和電荷載向量,力平衡方程如式(4)所示:

        (4)

        1.2 空間桁架結(jié)構(gòu)模型及動(dòng)力學(xué)分析

        本文研究的空間桁架簡化結(jié)構(gòu)如圖2所示,該結(jié)構(gòu)共由12個(gè)正六面體柔性桁架單元組成,共52個(gè)單元,161根桿件。桁架單元如圖3所示,每個(gè)單元均由空心球體、彈簧和空心桿組成,共8個(gè)節(jié)點(diǎn)、18根空心桿,其中包括12個(gè)邊長桿及6個(gè)對角桿??招那蝮w和空心桿材料均為不銹鋼。所有節(jié)點(diǎn)均為彈性平移鉸,節(jié)點(diǎn)內(nèi)部連接情況如圖4所示。

        圖2 空間桁架結(jié)構(gòu)

        圖3 桁架單元模型

        圖4 桁架節(jié)點(diǎn)內(nèi)部結(jié)構(gòu)

        由于桿件節(jié)點(diǎn)較多,先利用ANSYS命令流生成桁架結(jié)構(gòu)有限元模型,利用APDL語言編寫提取出單元編號及對應(yīng)節(jié)點(diǎn)編號和位置信息儲(chǔ)存到數(shù)組,將數(shù)組導(dǎo)入MATALAB編寫程序。

        桿端彈簧按照與桿件串聯(lián)處理,凸出的六面體認(rèn)為固定在航天器上,空心球按集中質(zhì)量處理,對結(jié)構(gòu)進(jìn)行模態(tài)分析,分析結(jié)果如表1所示。

        表1 空間桁架模態(tài)分析結(jié)果

        結(jié)果表明,前5階模態(tài)主要為X及Y方向的振動(dòng),因此在對結(jié)構(gòu)桁架進(jìn)行振動(dòng)控制時(shí)主要對X方向及Y方向進(jìn)行控制。

        2 主動(dòng)桿件優(yōu)化配置

        隨著大型航天結(jié)構(gòu)控制技術(shù)的發(fā)展,主動(dòng)桿件的位置優(yōu)化是一個(gè)重要問題。通過主動(dòng)桿件位置的選擇,可改善控制器的性能。合適的優(yōu)化位置不但影響控制效果,且影響控制系統(tǒng)能量的消耗。本文基于可控性度量指標(biāo)和逐步消減法進(jìn)行主動(dòng)桿件的優(yōu)化配置。

        將式(3)表達(dá)為狀態(tài)方程形式,離散化后如式(5)所示:

        (5)

        系統(tǒng)的可控性格蘭姆矩陣Wc如式(6)所示:

        (6)

        定義作動(dòng)器系統(tǒng)的可控性度量指標(biāo)如式(7)所示[12,13]:

        (7)

        指標(biāo)值越大,系統(tǒng)可控度越高。逐步消減法[12]是從剩余的作動(dòng)器的候選位置中去掉一個(gè)對目標(biāo)函數(shù)貢獻(xiàn)小的候選位置,直至目標(biāo)數(shù)目的候選位置,當(dāng)s個(gè)可選位置上配置有作動(dòng)器時(shí),矩陣Pc的維數(shù)為2ncs×2nc,則式(6)表達(dá)如下

        (8)

        pi為矩陣Pc的第i行數(shù)據(jù)。刪除一個(gè)作動(dòng)器的可選位置時(shí)有:

        (9)

        pj為Pc中與這個(gè)作動(dòng)器可選位置相關(guān)聯(lián)的2nc行數(shù)據(jù)。

        優(yōu)化μ1,μ2,μ3過程中,刪除μ1(Wcs-1),μ2(Wcs-1),μ3(Wcs-1)最大值相應(yīng)位置的桿件,以保證剩余桿件的配置可控度最大。逐步消減法主動(dòng)桿件優(yōu)化配置結(jié)果如表2所示。

        表2 逐步消減法主動(dòng)桿件優(yōu)化配置

        結(jié)果表明3種指標(biāo)的優(yōu)化結(jié)果接近,在控制第二階模態(tài)時(shí)桿件配置略有差別且高階模態(tài)配置時(shí)基本包括了低階模態(tài)配置的桿件。

        3 控制仿真與分析

        為驗(yàn)證主動(dòng)桿件優(yōu)化配置的有效性,對上述模型在主動(dòng)桿件優(yōu)化配置及主動(dòng)桿件選配位置兩種情況進(jìn)行振動(dòng)控制仿真計(jì)算。對式(3)進(jìn)行模態(tài)坐標(biāo)轉(zhuǎn)換,結(jié)果如下:

        (10)

        本文基于模態(tài)濾波技術(shù)[15],利用最優(yōu)模態(tài)控制算法實(shí)現(xiàn)控制仿真分析。動(dòng)力學(xué)方程解耦后,每階模態(tài)的控制利用LQR方法,其中權(quán)矩陣Qi表征了對第i階模態(tài)的勢能、動(dòng)能的控制要求,Ri反映了對主動(dòng)桿件的控制能量要求。Qi按下式選?。?/p>

        (11)

        ki,mi為相應(yīng)i階相應(yīng)的模態(tài)剛度,模態(tài)質(zhì)量,ai,bi為對應(yīng)的權(quán)系數(shù),通過調(diào)整權(quán)系數(shù)可以實(shí)現(xiàn)不同的控制要求。

        3.1 仿真計(jì)算

        分別在桁架節(jié)點(diǎn)1,2,37,38,13,14,25,26,23,24,35,36 處的X,Y,Z三個(gè)方向施加0.1mm位移,節(jié)點(diǎn)11,12,47,48處的X,Y,Z方向施加-0.1mm位移。

        根據(jù)優(yōu)化配置結(jié)果選取6,7,17,113號桿件作為主動(dòng)桿件。節(jié)點(diǎn)1位移響應(yīng)如圖5所示。為檢驗(yàn)優(yōu)化配置的效果,選配39,95,117,146號桿件為主動(dòng)桿件,取相同的控制權(quán)系數(shù),此時(shí)節(jié)點(diǎn)1位移響應(yīng)如圖6所示,兩種配置時(shí)的控制電壓分別如圖7,圖8所示。

        圖5 優(yōu)化配置時(shí)1節(jié)點(diǎn)位移響應(yīng)

        圖6 選配位置時(shí)1節(jié)點(diǎn)位移響應(yīng)

        圖7 優(yōu)化位置控制電壓

        圖8 選配位置控制電壓

        圖5表明,XY向振動(dòng)得到明顯抑制,但Z向振動(dòng)控制效果不明顯,這是因?yàn)榍八碾A沒有Z向的主振動(dòng),結(jié)構(gòu)自身衰減較快。圖6表明,XY向振動(dòng)得到抑制,但Z向振動(dòng)出現(xiàn)較大過調(diào)。圖7表明,優(yōu)化位置控制電壓峰值在400V左右,而圖8表明,選配位置控制電壓峰值接近20 000V,因此優(yōu)化位置控制電壓分布更合理。

        圖5~圖8結(jié)果表明,優(yōu)化配置時(shí)的控制效果明顯優(yōu)于選配位置。

        3.2 控制仿真效果評價(jià)

        根據(jù)式(12)確定控制仿真過程中結(jié)構(gòu)反應(yīng)的峰值和均方值的評價(jià)指標(biāo):

        (12)

        仿真過程中作動(dòng)器評價(jià)指標(biāo)如公式(13)所示:

        (13)

        其中:s為作動(dòng)器個(gè)數(shù),W(t)=1 000V為擊穿電壓。s和W(t)可表征電壓峰值和控制能量。

        對初始位移擾動(dòng),位移衰減時(shí)間指標(biāo)具體定義為:

        (14)

        其中:tki,ti分別定義為第i節(jié)點(diǎn)施加控制和沒有控制時(shí)位移響應(yīng)衰減到結(jié)構(gòu)無控時(shí)該點(diǎn)位移響應(yīng)最大值5%所需的時(shí)間,為結(jié)構(gòu)節(jié)點(diǎn)總數(shù)。Jt可綜合評價(jià)結(jié)構(gòu)振動(dòng)衰減快慢。

        相同初始條件或相同激勵(lì)信號下,初始位移擾動(dòng)的控制評價(jià)結(jié)果見表3。

        按照定義的評價(jià)指標(biāo),對正弦激勵(lì),隨機(jī)激勵(lì)的振動(dòng)控制進(jìn)行驗(yàn)證。選配124號桿件為激勵(lì)桿件,同時(shí)配置6,7,17,113號桿件作為控制桿件進(jìn)行正弦激勵(lì)及隨機(jī)激勵(lì)的控制仿真效果評價(jià)。正弦激勵(lì)電壓信號假定為:

        隨機(jī)激勵(lì)電壓信號為0-200V平穩(wěn)隨機(jī)信號,評價(jià)結(jié)果見表4。

        從表3,表4可以清楚地比較控制效果,可以發(fā)現(xiàn)優(yōu)化位置的各項(xiàng)指標(biāo)值明顯優(yōu)于非優(yōu)化位置,振動(dòng)抑制效果明顯,控制能耗低,整體結(jié)構(gòu)的振動(dòng)控制效果良好。評價(jià)結(jié)果表明本文采用的桿件優(yōu)化措施和控制算法對大型桁架結(jié)構(gòu)的整體振動(dòng)控制是非常有效的。對于其他的不同邊界條件,可以按照同樣的方法給出評價(jià),也可以對不同的控制算法進(jìn)行控制效果評價(jià)。

        表3 初始位移擾動(dòng)振動(dòng)控制評價(jià)

        (1. X向初始位移擾動(dòng)為節(jié)點(diǎn)1X向施加0.1mm初始位移,節(jié)點(diǎn)24X向施加-0.1mm初始位移; 2. X、Y向初始位移擾動(dòng)為節(jié)點(diǎn)1X、Y向施加0.1mm初始位移,節(jié)點(diǎn)24X、Y向施-0.1mm初始位移)

        表4 正弦、隨機(jī)振動(dòng)控制評價(jià)

        4 結(jié)論

        本文研究了空間柔性桁架結(jié)構(gòu)的振動(dòng)抑制問題,完成了結(jié)構(gòu)動(dòng)力學(xué)建模、主動(dòng)桿件的優(yōu)化配置以及振動(dòng)主動(dòng)控制仿真,主要結(jié)論如下:

        1)建立了考慮剛性連接和彈性連接的壓電桿件的機(jī)電耦合模型,并采用ANSYS、Matlab編制了空間柔性桁架結(jié)構(gòu)動(dòng)力學(xué)分析程序,為動(dòng)力學(xué)分析和控制研究提供了方法和研究基礎(chǔ)。

        2)提出了可控性度量指標(biāo)與逐步消減法相結(jié)合的優(yōu)化方法,實(shí)現(xiàn)了空間桁架結(jié)構(gòu)主動(dòng)桿件的優(yōu)化配置。

        3)采用模態(tài)最優(yōu)控制算法實(shí)現(xiàn)了空間桁架結(jié)構(gòu)的振動(dòng)控制仿真,表明優(yōu)化配置的振動(dòng)抑制效果優(yōu)于非優(yōu)化配置,并且整體結(jié)構(gòu)振動(dòng)控制效果良好。

        [1] 黃文虎,王心清,張景繪,等. 航天柔性結(jié)構(gòu)振動(dòng)控制的若干新進(jìn)展[J]. 力學(xué)進(jìn)展,1997,27(1):5-18.

        [2] 王其政,劉 斌,宋文濱.航天事故與動(dòng)力學(xué)環(huán)境預(yù)示和控制技術(shù)研究述評[J].環(huán)境技術(shù),1995,4(1):1-6.

        [3] 鄭 凱,陳 杰,楊義勇. 壓電自適應(yīng)桁架結(jié)構(gòu)智能振動(dòng)控制[J]. 控制理論與應(yīng)用,2010,27(7):943-947.

        [4] 曹玉巖,付世欣,王鳴浩. 壓電智能桁架結(jié)構(gòu)的建模與最優(yōu)振動(dòng)控制[J]. 壓電與聲光,2014,36(4):523-527.

        [5] 陳文英,褚福磊,閻紹澤. 智能桁架結(jié)構(gòu)自適應(yīng)模糊主動(dòng)振動(dòng)控制[J]. 清華大學(xué)學(xué)報(bào),2008,48(5):816-819.

        [6] Marinaki M, Marinakis Y, Stavroulakis G E. Fuzzy control optimized by PSO for vibration suppression of beams[J]. Control Eng Pract, 2010, 18(6): 618-629.

        [7] Sun D Q. Hierarchical fuzzy generalized predictive control and its application in vibration suppression of flexible spacecraft[J]. J Astronautics, 2009, 30(Compendex): 1537-1543.

        [8] Zhu L, Song W, Hu Q. Active vibration suppression and attitude maneuvers of flexible spacecraft via fuzzy sliding control[A]. Proceedings of the Chinese Control and Decision Conference[C]. Yantai, 2008.

        [9] Armaghan Salehian, T.Michael Seigler, Daniel.J.Inman. Control of the continuum model of a large flexible space structure[A]. Proceedings of IMECE2006[C]. IMECE2006-13362.

        [10] 李東旭,劉 望,蔣建平,等. 空間智能桁架的傳感器作動(dòng)器位置優(yōu)化和分散化自適應(yīng)模糊振動(dòng)控制[J]. 中國科學(xué),2011,41(5):602-610.

        [11] 劉福強(qiáng),張令彌.作動(dòng)器與傳感器優(yōu)化配置的逐步消減法[J].宇航學(xué)報(bào),2000(3),64-69.

        [12] Xing G Q. Bainum P M. Actuator Placement Using Degree of Controllability for Discrete-Time Systems[J]. J.dyn.syst.meas. Control, 1992, 114(3), 508-516.

        [13] Sylvaine L, Hisham A K, Yvan B, et al. Piezoelectric actuators and sensors location for active cont rol of flexible structures[J]. IEEE Transaction On Instrument and Measurement, 2001, 50(6), 1577-1582.

        [14] 趙國忠,顧元憲.壓電智能桁架的振動(dòng)控制和結(jié)構(gòu)的一體化優(yōu)化設(shè)計(jì)[J].航空學(xué)報(bào),2003,24(4),332-335.

        [15] 歐進(jìn)萍.結(jié)構(gòu)振動(dòng)控制—主動(dòng)、半主動(dòng)和智能控制[M].北京:科學(xué)出版社,2003.

        Research on Optimal Allocation and Vibration Control of Space Truss Structure

        Xu Jianguo1, Shao Kang2, Zhang Jingjing2, Dong Longlei2

        (1.Nanjing SEU Energy Engineering Institute, Nanjing 211100,China;2.School of Aerospace, Xi’an Jiaotong University, Xi’an 710049,China)

        Large space truss structure is widely used in space crafts. The vibration of this kind of structure may lead to some serious problems. It is necessary to execute vibration control for large space truss structure. In this paper, combine the integrated controllability index and stepwise subtraction method to optimize the position of active rods for large truss structure to solve the vibration control problem. The electro-mechanical coupled equations of the active rod are established. A space truss structure is taken as an example and verification simulations are carried out. The results evaluated by some indexes used in Benchmark questions of vibration control show that the method mentioned in this paper can control the vibration effectively.

        space truss; vibration control; optimal placement; modal control

        2015-09-12;

        2015-10-21。

        許建國(1983-),男,江蘇姜堰人,碩士,工程師,主要從事結(jié)構(gòu)動(dòng)力學(xué)及振動(dòng)控制方向的研究。

        董龍雷(1973-),男,山東人,博士,副教授,主要從事綜合力學(xué)環(huán)境仿真、試驗(yàn)與控制以及振動(dòng)控制方向的研究。

        1671-4598(2016)03-0067-04

        10.16526/j.cnki.11-4762/tp.2016.03.019

        TB123

        A

        猜你喜歡
        作動(dòng)器桿件壓電
        基于臨時(shí)支撐結(jié)構(gòu)的桿件初彎曲對其軸壓性能的影響
        四川建筑(2021年1期)2021-03-31 01:01:46
        基于火箭發(fā)射起豎的機(jī)電伺服作動(dòng)器質(zhì)量特性研究
        塔式起重機(jī)拼裝式超長附著桿設(shè)計(jì)與應(yīng)用
        《壓電與聲光》征稿啟事
        壓電與聲光(2019年1期)2019-02-22 09:46:06
        新型壓電疊堆泵設(shè)計(jì)及仿真
        作動(dòng)器防擺動(dòng)控制結(jié)構(gòu)設(shè)計(jì)
        基于混合潤滑理論的航空作動(dòng)器密封性能分析
        KD379:便攜折疊式衣架
        某網(wǎng)架桿件彎曲的原因分析及處理
        航空作動(dòng)器的VL密封特性分忻
        天码av无码一区二区三区四区 | 美女胸又www又黄的网站| 国产成人一区二区三区高清| 亚洲不卡高清av在线| 精品一二三四区中文字幕| 狠狠色噜噜狠狠狠狠7777米奇| 亚洲综合欧美在线| 一本久久综合亚洲鲁鲁五月夫| 日本免费在线一区二区三区| 国产精品人妻一码二码| 欧美亚洲综合另类| 伊人影院在线观看不卡| av免费不卡一区二区| 成年女人色毛片| 国产成人精品三级91在线影院 | 97超碰中文字幕久久| 国产成人精品免费久久久久| 亚洲男人的天堂网站| 亚洲国产一区二区三区在线视频| 日本成人中文字幕亚洲一区 | 欧美性猛交xxxx富婆| 女性自慰网站免费看ww| 国产av大片久久中文字幕| 色综合天天综合欧美综合| 无码任你躁久久久久久| 国产欧美亚洲另类第一页| 成人免费av高清在线| 99国产精品自在自在久久| 久久精品中文字幕第23页| 国产目拍亚洲精品二区| 色翁荡息又大又硬又粗视频| 中国a级毛片免费观看| 国产日韩欧美视频成人| 水蜜桃男女视频在线观看网站| 亚洲熟妇少妇任你躁在线观看无码 | 久草青青91在线播放| 人妻少妇精品中文字幕av蜜桃| 国产91在线精品福利| 五月婷婷开心六月激情| 正在播放国产多p交换视频| 精品午夜久久网成年网|