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

    基于基函數(shù)分解的連續(xù)變速顫振邊界預(yù)測方法

    2014-07-09 17:56:09譚博鄭華裴承鳴
    現(xiàn)代電子技術(shù) 2014年13期

    譚博+鄭華+裴承鳴

    摘 要: 針對連續(xù)變速顫振試驗(yàn)實(shí)測信號的特點(diǎn)及在線處理需求,提出一種基于基函數(shù)展開的時(shí)變參數(shù)建模方法。通過對結(jié)構(gòu)響應(yīng)信號建模和穩(wěn)定性判據(jù)的趨勢分析,得到顫振邊界隨時(shí)間變化的預(yù)測曲線。依據(jù)顫振試驗(yàn)機(jī)理,采用白噪聲激勵(lì)多模態(tài)耦合系統(tǒng)的方式生成了相應(yīng)的仿真信號,進(jìn)而在不同信噪比下驗(yàn)證了方法的數(shù)值性能。最后通過氣彈模型低速風(fēng)洞試驗(yàn)進(jìn)一步檢驗(yàn)本文方法的工程實(shí)用性。

    關(guān)鍵詞: 連續(xù)變速顫振試驗(yàn); 基函數(shù); 時(shí)變參數(shù)建模; 顫振邊界預(yù)測

    中圖分類號: TN911.7?34; TP391.9 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號: 1004?373X(2014)13?0047?03

    Method for progression variable speed flutter boundary prediction

    based on basis function decomposition

    TAN Bo, ZHENG Hua, PEI Cheng?ming

    (Data Processing Center, Northwestern Polytechnical University, Xian 710072, China)

    Abstract: According to the characteristics and online processing requirements of the actual measuring signal in flutter test with progression variable speed (FTPVS), a time?varying parameter modeling algorithm based on basis function expansion is presented in this paper. The prediction curves of flutter boundary versus time variation were obtained on the basis of modeling of structural response signal and trend analysis of the stability criterion. According to the mechanism of flutter test, the correspon?ding simulation signal is generated with the mode that the white noise stimulates the multi?modal coupling system. The numerical performance of the method was verified at different SNRs. The practical applicability of the method was checked in aeroelastic model wind?tunnel test.

    Keywords: flutter test with progression variable speed; basis function; time?varying parameter modeling; flutter boundary prediction

    0 引 言

    連續(xù)變速顫振試驗(yàn)是近年來從飛機(jī)結(jié)構(gòu)強(qiáng)度專業(yè)角度提出的一種全新顫振試驗(yàn)概念,指由于研究對象的實(shí)際工作狀態(tài)或試驗(yàn)方法的需要,馬赫數(shù)、速度、高度等參數(shù)連續(xù)變化的一類顫振試驗(yàn)技術(shù)。與傳統(tǒng)的臺(tái)階式顫振試驗(yàn)方法相比,具有試驗(yàn)周期短、試驗(yàn)開支小、技術(shù)約束少、更符合實(shí)際使用狀態(tài)等優(yōu)勢,但同時(shí)也存在更大的試驗(yàn)風(fēng)險(xiǎn)[4]。因此,在試驗(yàn)中能否基于實(shí)測信號快速準(zhǔn)確地完成顫振邊界預(yù)測是保證該類試驗(yàn)安全有效的關(guān)鍵技術(shù)之一。

    連續(xù)變速顫振試驗(yàn)的結(jié)構(gòu)響應(yīng)信號是非平穩(wěn)隨機(jī)過程,這一特點(diǎn)使得傳統(tǒng)的顫振試驗(yàn)數(shù)據(jù)處理方法難以直接應(yīng)用。因此,本文提出一種基于基函數(shù)展開的時(shí)變參數(shù)建模方法,通過實(shí)測信號建立測試對象的時(shí)變參數(shù)模型,從而提取信號中所蘊(yùn)含的穩(wěn)定性特征量,并根據(jù)該特征量的趨勢進(jìn)行外推分析,以獲取每個(gè)時(shí)刻的顫振邊界預(yù)測值,為試驗(yàn)過程的實(shí)時(shí)監(jiān)控及試驗(yàn)安全提供技術(shù)保障。

    1 方法簡述

    在平穩(wěn)系統(tǒng)模型建立中,非時(shí)變隨機(jī)信號可以看作是線性系統(tǒng)的白噪聲激勵(lì)響應(yīng),如下式:

    [y(n)=-k=1paky(n-k)+k=0qbku(n-k)] (1)

    式中:[y(n)]為非時(shí)變隨機(jī)信號;[u(n)]為系統(tǒng)的白噪聲激勵(lì)。

    相應(yīng)的,由于在連續(xù)變速顫振試驗(yàn)中,試驗(yàn)對象的結(jié)構(gòu)模態(tài)參數(shù)會(huì)隨時(shí)間發(fā)生變化,因此,對公式(1)中的系數(shù)進(jìn)行改變,即可得到時(shí)變線性系統(tǒng)的參數(shù)模型,表示為:

    [y(n)=-k=1pak(n)y(n-k)+k=0qbk(n)u(n-k)] (2)

    公式(2)也被稱為非平穩(wěn)隨機(jī)信號的時(shí)變自回歸滑動(dòng)平均模型。若系數(shù)[bk(n)=0, k=1,2,3,…,]則模型變?yōu)榉瞧椒€(wěn)隨機(jī)信號的時(shí)變自回歸模型(TVAR)。

    類似于平穩(wěn)系統(tǒng)模型的建立,時(shí)變隨機(jī)信號可以看作非平穩(wěn)線性系統(tǒng)的白噪聲激勵(lì)響應(yīng),在實(shí)際應(yīng)用中,通常采用時(shí)變自回歸模型,由公式⑵可以得到TVAR的參數(shù)模型為:

    [y(n)=-k=1pak(n)y(n-k)+u(n)] (3)

    式中:[y(n)]是時(shí)變隨機(jī)信號;[p]是模型階數(shù);[ak(n)]是時(shí)變AR系數(shù);[u(n)]是白噪聲激勵(lì)。

    為求解時(shí)變AR系數(shù),將[ak(n)]用函數(shù)空間基函數(shù)的加權(quán)和表示為:

    [ak(n)=i=0mbkifi(n)] (4)

    式中:[bki]為加權(quán)系數(shù);[fi(n)]為選定函數(shù)空間的基函數(shù);[m]為函數(shù)空間的維數(shù)。

    相應(yīng)的TVAR模型公式變?yōu)椋?/p>

    [y(n)=-k=1pi=0mbki[fi(n)y(n-k)]+u(n)] (5)

    定義矢量[Y(N-k)]和[θ:]

    [Y(N-k)=[f0(n)y(n-k),…,fm(n)y(n-k)]] (6)

    [θ=[a10,…,a1m,…,ap0,…,qpm]] (7)

    這時(shí),TVAR模型的矩陣形式為:

    [y(n)=-[Y(N-1),…,Y(n-p)]θT+u(n)] (8)

    由此,將線性時(shí)變AR系數(shù)[ak(n)]的求解轉(zhuǎn)化成為求解線性時(shí)不變系統(tǒng)的加權(quán)系數(shù)向量[θ]的問題。解出[θ]后,結(jié)合選定函數(shù)空間的基函數(shù)[fi(n),]由公式(4)即可得到所建模型的時(shí)變AR系數(shù)。

    本文選擇傅里葉基函數(shù)對時(shí)變AR系數(shù)進(jìn)行求解,其構(gòu)造公式為:

    [fk(n)=coskπn2N,k為偶數(shù)sin(k+1)πn2N,k為奇數(shù)] (9)

    由解得的時(shí)變AR系數(shù)[ak(n),]可以提取[n]時(shí)刻的穩(wěn)定性判據(jù),如Jury判據(jù),Lyapunov判據(jù)以及阻尼比系數(shù)等。依據(jù)[n]時(shí)刻對應(yīng)的速度值,以曲線擬合的方式對顫振邊界進(jìn)行預(yù)測。

    2 仿真試驗(yàn)及結(jié)果

    根據(jù)典型結(jié)構(gòu)顫振機(jī)理及其響應(yīng)信號特征,采用白噪聲激勵(lì)時(shí)變系統(tǒng)生成的仿真信號來測試本文方法的數(shù)值性能。

    不排除一般性,典型的顫振現(xiàn)象往往在兩階結(jié)構(gòu)模態(tài)耦合時(shí)出現(xiàn),為此,模擬連續(xù)變速顫振響應(yīng)信號可由以下系統(tǒng)生成:

    [x(t)=i=1nAie-ξi(t)fi(t)2πtcos2πfi(t)1-ξ2i(t)t] (10)

    式中:[n]為模態(tài)個(gè)數(shù);[f]為頻率;[ξ]為阻尼比。在0時(shí)刻,各參數(shù)的初值分別為:[f1=]10 Hz,[f2=]20 Hz,[v=0,][ξ1=0.11,ξ2=0.1。]

    為了模擬顫振發(fā)生過程,在歸一化速度線性增加的情況下,兩階模態(tài)的頻率逐漸靠攏直至完全重合,而阻尼比則逐漸衰減至零。經(jīng)過64 s運(yùn)行后仿真系統(tǒng)發(fā)生顫振,在采樣頻率為128 Hz情況下得到長度為8 192的隨機(jī)響應(yīng)信號,此時(shí)對應(yīng)的頻率為[f1=f2=]15 Hz,相應(yīng)的阻尼比系數(shù)為[ξ1=ξ2=0,]歸一化顫振速度為[v=1。]

    在上述條件下生成的無噪聲仿真信號的時(shí)間歷程如圖1所示,相應(yīng)的聯(lián)合時(shí)頻分布如圖2所示。

    圖1 仿真信號的時(shí)間歷程

    圖2 仿真信號的聯(lián)合時(shí)頻圖

    應(yīng)用本文方法對無噪聲干擾的仿真信號進(jìn)行顫振邊界外推計(jì)算,所得結(jié)果如圖3所示。圖中,縱軸為速度,[v=1]的直線標(biāo)示了顫振邊界的真值,橫軸為時(shí)間。由圖3可以看出,在無噪聲的情況下,本文方法計(jì)算的外推值從時(shí)刻[n=40]開始快速上升在[n=50]時(shí)到達(dá)真實(shí)值[v=1,]在其后的時(shí)間內(nèi),外推值在真值附近小幅波動(dòng),也就是說可以快速得到準(zhǔn)確穩(wěn)定的顫振邊界預(yù)測值。

    為了進(jìn)一步測試本文方法在噪聲環(huán)境下的數(shù)值性能,定義信噪比:

    [SNR=10×logPsignalPnoise] (11)

    式中:[P]為信號的功率,在此定義下,計(jì)算并得到了信噪比分別為20 dB,14 dB,10 dB和8 dB的外推值曲線,如圖4所示。

    圖3 無噪聲環(huán)境的顫振邊界預(yù)測結(jié)果

    圖4 噪聲環(huán)境下顫振邊界預(yù)測結(jié)果

    由圖4可以看出,信噪比為20 dB時(shí),外推值曲線盡管隨時(shí)間變化的趨勢有少許不同,但還是在[n=50]左右達(dá)到了真實(shí)值。當(dāng)信噪比下降到14 dB時(shí),外推值趨勢與無噪聲情況類似,但最終外推值在到達(dá)[v=]1.2附近時(shí)才趨于穩(wěn)定,顫振邊界預(yù)測值的精度也有所下降。若信噪比進(jìn)一步降低,可以看到外推值的平穩(wěn)趨勢開始消失,所得的外推結(jié)果已經(jīng)不宜在試驗(yàn)中作為參考。由此可見,本文方法對噪聲相對敏感。為此,在實(shí)際應(yīng)用中應(yīng)盡可能保證測試、采集、記錄等響應(yīng)信號獲取過程各環(huán)節(jié)的質(zhì)量。

    3 試驗(yàn)數(shù)據(jù)處理

    為了驗(yàn)證所研究方法的工程實(shí)用性,先后進(jìn)行了多次物理試驗(yàn)。作為舉例,這里僅給出一組試驗(yàn)結(jié)果。

    試驗(yàn)數(shù)據(jù)取自某飛機(jī)氣彈模型低速風(fēng)洞試驗(yàn),對應(yīng)試驗(yàn)狀態(tài)由傳統(tǒng)臺(tái)階等速試驗(yàn)得到的顫振臨界速度為40.8 m/s。

    實(shí)測信號的時(shí)間歷程及對應(yīng)的試驗(yàn)風(fēng)速變化曲線如圖5所示。

    圖5 連續(xù)變速顫振試驗(yàn)實(shí)測數(shù)據(jù)

    由圖5可見,選取的實(shí)測信號對應(yīng)的風(fēng)速以近似于線性增長的方式,經(jīng)過29 s從25 m/s提升至38 m/s。采用本文方法對該段實(shí)測數(shù)據(jù)進(jìn)行處理,所求的顫振邊界預(yù)測值曲線如圖6所示。

    圖6 實(shí)測數(shù)據(jù)的顫振邊界預(yù)測結(jié)果

    由圖6可以看出,本文方法的計(jì)算結(jié)果,于20 s左右達(dá)到了顫振臨界速度[v=]40.8 m/s附近,并趨于平穩(wěn)。此時(shí)對應(yīng)的速度為[v=]35 m/s,距離顫振速度尚有較大的安全空間。這表明,本文方法在連續(xù)變速顫振試驗(yàn)中,可以對顫振邊界進(jìn)行有效預(yù)測,并且保留了較大的安全裕量。

    4 結(jié) 論

    本文研究了基于基函數(shù)展開的時(shí)變參數(shù)建模方法及其在連續(xù)變速顫振試驗(yàn)中的應(yīng)用問題。通過數(shù)值仿真,驗(yàn)證并分析了方法的精度及抗噪性能。結(jié)果表明,該方法在較高信噪比下,可以得到滿意的預(yù)測值,但其精度會(huì)隨信噪比的減小而下降。而氣彈模型風(fēng)洞試驗(yàn)的應(yīng)用結(jié)果表明,該方法完全可以滿足工程實(shí)際的物理需求。

    本文研究為豐富和擴(kuò)展連續(xù)變速顫振試驗(yàn)的在線監(jiān)控方法提供了一種新的手段,在航空、航天、導(dǎo)彈、亞太空等領(lǐng)域的結(jié)構(gòu)顫振試驗(yàn)方面有著潛在的應(yīng)用前景。

    參考文獻(xiàn)

    [1] COWAN T J, ARENA A S, GUPTA K K. Accelerating CFD?based aero elastic predictions using system identification [R/OL]. [2010?01?02]. http://www. citeseerx.ist.psu.edu.

    [2] NAMARA J J, FRIEDMANN P P. Flutter boundary identification for time?domain computational aero elasticity [J]. American Institute of Aeronautics and Astronautics. 2007, 45(7): 1546?1554.

    [3] MAHER Damien, Combined time and frequency domain approaches to the operational identification of vehicle suspension [D]. Ireland: School of Mechanical and Manufacturing Engineering, Dublin City University, 2011.

    [4] Garrick I E, Reed W H. Historical development of aircraft flutter [J]. Journal of Aircraft, 1981, 18(11): 897?912.

    [5] BAE J, KIM J, LEE I, et al. Extension of flutter prediction parameter for multimode flutter systems [J] Journal of Aircraft, 2005, 42(1): 285?288.

    [6] MATSUZAKI Y, TORII H. Flutter boundary prediction of an adaptive smart wing during process of adaptation using steady?state response data [C]// 47th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference. Newport, RI: AIAA, 2006: 2132?2139.

    [7] THURUTHIMATTAM, B J, FRIEDMANN P P, MCNAMARA J J, et al. Modeling approaches to hypersonic aerothermoelasticity with application to reusable launch vehicles [C]// 44th AIAA/ASME/ASCE/AHS Structures, Structural Dynamics and Materials Conference. Norfolk, VA: AIAA, 2003: 1967?1977.

    [8] 朝倫巴根,賈德彬.數(shù)值計(jì)算方法[M].北京:中國水利水電出版社,2006.

    [9] 王宏禹.非平穩(wěn)隨機(jī)信號分析與處理[M].北京:國防工業(yè)出版社,1999.

    [10] 林青,戴慧珺,馬文濤.基于正交基函數(shù)神經(jīng)網(wǎng)絡(luò)的圖像加密算法仿真[J].計(jì)算機(jī)仿真,2013,30(10):416?421.

    參考文獻(xiàn)

    [1] COWAN T J, ARENA A S, GUPTA K K. Accelerating CFD?based aero elastic predictions using system identification [R/OL]. [2010?01?02]. http://www. citeseerx.ist.psu.edu.

    [2] NAMARA J J, FRIEDMANN P P. Flutter boundary identification for time?domain computational aero elasticity [J]. American Institute of Aeronautics and Astronautics. 2007, 45(7): 1546?1554.

    [3] MAHER Damien, Combined time and frequency domain approaches to the operational identification of vehicle suspension [D]. Ireland: School of Mechanical and Manufacturing Engineering, Dublin City University, 2011.

    [4] Garrick I E, Reed W H. Historical development of aircraft flutter [J]. Journal of Aircraft, 1981, 18(11): 897?912.

    [5] BAE J, KIM J, LEE I, et al. Extension of flutter prediction parameter for multimode flutter systems [J] Journal of Aircraft, 2005, 42(1): 285?288.

    [6] MATSUZAKI Y, TORII H. Flutter boundary prediction of an adaptive smart wing during process of adaptation using steady?state response data [C]// 47th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference. Newport, RI: AIAA, 2006: 2132?2139.

    [7] THURUTHIMATTAM, B J, FRIEDMANN P P, MCNAMARA J J, et al. Modeling approaches to hypersonic aerothermoelasticity with application to reusable launch vehicles [C]// 44th AIAA/ASME/ASCE/AHS Structures, Structural Dynamics and Materials Conference. Norfolk, VA: AIAA, 2003: 1967?1977.

    [8] 朝倫巴根,賈德彬.數(shù)值計(jì)算方法[M].北京:中國水利水電出版社,2006.

    [9] 王宏禹.非平穩(wěn)隨機(jī)信號分析與處理[M].北京:國防工業(yè)出版社,1999.

    [10] 林青,戴慧珺,馬文濤.基于正交基函數(shù)神經(jīng)網(wǎng)絡(luò)的圖像加密算法仿真[J].計(jì)算機(jī)仿真,2013,30(10):416?421.

    參考文獻(xiàn)

    [1] COWAN T J, ARENA A S, GUPTA K K. Accelerating CFD?based aero elastic predictions using system identification [R/OL]. [2010?01?02]. http://www. citeseerx.ist.psu.edu.

    [2] NAMARA J J, FRIEDMANN P P. Flutter boundary identification for time?domain computational aero elasticity [J]. American Institute of Aeronautics and Astronautics. 2007, 45(7): 1546?1554.

    [3] MAHER Damien, Combined time and frequency domain approaches to the operational identification of vehicle suspension [D]. Ireland: School of Mechanical and Manufacturing Engineering, Dublin City University, 2011.

    [4] Garrick I E, Reed W H. Historical development of aircraft flutter [J]. Journal of Aircraft, 1981, 18(11): 897?912.

    [5] BAE J, KIM J, LEE I, et al. Extension of flutter prediction parameter for multimode flutter systems [J] Journal of Aircraft, 2005, 42(1): 285?288.

    [6] MATSUZAKI Y, TORII H. Flutter boundary prediction of an adaptive smart wing during process of adaptation using steady?state response data [C]// 47th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference. Newport, RI: AIAA, 2006: 2132?2139.

    [7] THURUTHIMATTAM, B J, FRIEDMANN P P, MCNAMARA J J, et al. Modeling approaches to hypersonic aerothermoelasticity with application to reusable launch vehicles [C]// 44th AIAA/ASME/ASCE/AHS Structures, Structural Dynamics and Materials Conference. Norfolk, VA: AIAA, 2003: 1967?1977.

    [8] 朝倫巴根,賈德彬.數(shù)值計(jì)算方法[M].北京:中國水利水電出版社,2006.

    [9] 王宏禹.非平穩(wěn)隨機(jī)信號分析與處理[M].北京:國防工業(yè)出版社,1999.

    [10] 林青,戴慧珺,馬文濤.基于正交基函數(shù)神經(jīng)網(wǎng)絡(luò)的圖像加密算法仿真[J].計(jì)算機(jī)仿真,2013,30(10):416?421.

    欧美另类一区| 波多野结衣av一区二区av| 国产精品久久久人人做人人爽| 欧美激情极品国产一区二区三区| 亚洲一区二区三区欧美精品| 亚洲综合色网址| 午夜日韩欧美国产| 12—13女人毛片做爰片一| 国产淫语在线视频| 51午夜福利影视在线观看| 亚洲伊人久久精品综合| 国产亚洲精品一区二区www | 韩国高清视频一区二区三区| 精品人妻在线不人妻| 欧美xxⅹ黑人| 中文精品一卡2卡3卡4更新| 欧美乱码精品一区二区三区| tube8黄色片| 久久女婷五月综合色啪小说| 精品熟女少妇八av免费久了| 啦啦啦视频在线资源免费观看| 热re99久久精品国产66热6| 国产精品久久久人人做人人爽| 69av精品久久久久久 | 欧美国产精品va在线观看不卡| 亚洲精品一区蜜桃| 亚洲中文av在线| 国产欧美日韩一区二区精品| 日日摸夜夜添夜夜添小说| 不卡一级毛片| 亚洲精品粉嫩美女一区| 国产麻豆69| 丁香六月欧美| 亚洲精品一区蜜桃| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 丝袜脚勾引网站| 国产xxxxx性猛交| 国产伦人伦偷精品视频| 国产精品免费大片| 国产视频一区二区在线看| 欧美变态另类bdsm刘玥| 成人手机av| 青青草视频在线视频观看| 日韩一区二区三区影片| 99re6热这里在线精品视频| 超色免费av| 搡老乐熟女国产| 少妇 在线观看| 国产主播在线观看一区二区| 在线观看免费日韩欧美大片| 国产99久久九九免费精品| 亚洲精品国产区一区二| 日韩熟女老妇一区二区性免费视频| 国产av又大| 两个人看的免费小视频| 韩国精品一区二区三区| 亚洲中文av在线| 99久久国产精品久久久| 亚洲精品国产av成人精品| 精品人妻一区二区三区麻豆| 最近中文字幕2019免费版| 夜夜夜夜夜久久久久| 午夜福利影视在线免费观看| 精品一品国产午夜福利视频| 久久精品成人免费网站| 免费av中文字幕在线| 97人妻天天添夜夜摸| 午夜福利乱码中文字幕| 波多野结衣一区麻豆| 我的亚洲天堂| 国产精品影院久久| 男女免费视频国产| 丰满迷人的少妇在线观看| 亚洲成人免费电影在线观看| 国产精品二区激情视频| 真人做人爱边吃奶动态| 国产精品1区2区在线观看. | 亚洲欧美成人综合另类久久久| 正在播放国产对白刺激| 欧美性长视频在线观看| 纯流量卡能插随身wifi吗| 日韩电影二区| 亚洲av日韩精品久久久久久密| 交换朋友夫妻互换小说| 在线观看免费午夜福利视频| 婷婷丁香在线五月| 色视频在线一区二区三区| 国产无遮挡羞羞视频在线观看| 亚洲av美国av| 亚洲伊人久久精品综合| 一本久久精品| 不卡一级毛片| 亚洲精品美女久久av网站| av网站免费在线观看视频| 久久香蕉激情| 老司机福利观看| 欧美国产精品va在线观看不卡| 欧美97在线视频| 啪啪无遮挡十八禁网站| 国产成人av激情在线播放| 女人被躁到高潮嗷嗷叫费观| 国产成人精品在线电影| 爱豆传媒免费全集在线观看| 亚洲五月色婷婷综合| 黄片小视频在线播放| 欧美日韩视频精品一区| 黄片播放在线免费| 欧美性长视频在线观看| 国产老妇伦熟女老妇高清| 日本一区二区免费在线视频| 亚洲精品乱久久久久久| 亚洲av成人一区二区三| 亚洲精品国产av蜜桃| 亚洲精品一二三| 国产成人精品在线电影| 最新的欧美精品一区二区| 欧美另类一区| 91精品国产国语对白视频| 下体分泌物呈黄色| 国产精品1区2区在线观看. | 99精品欧美一区二区三区四区| 国产免费现黄频在线看| 桃红色精品国产亚洲av| 婷婷丁香在线五月| 亚洲欧美清纯卡通| 涩涩av久久男人的天堂| 久久女婷五月综合色啪小说| 一本—道久久a久久精品蜜桃钙片| 91精品伊人久久大香线蕉| 国产亚洲一区二区精品| 他把我摸到了高潮在线观看 | 老司机影院毛片| 欧美成狂野欧美在线观看| 精品一区二区三区av网在线观看 | 精品少妇内射三级| 婷婷成人精品国产| 久久精品aⅴ一区二区三区四区| 狠狠狠狠99中文字幕| 操美女的视频在线观看| 看免费av毛片| 亚洲中文字幕日韩| 日韩欧美一区视频在线观看| 天堂中文最新版在线下载| 久久这里只有精品19| 久久国产亚洲av麻豆专区| 老司机靠b影院| 亚洲国产日韩一区二区| 51午夜福利影视在线观看| 国产男女超爽视频在线观看| 久久中文字幕一级| 久久久久久久久久久久大奶| 亚洲激情五月婷婷啪啪| 久久天躁狠狠躁夜夜2o2o| 亚洲色图 男人天堂 中文字幕| 久久久久国内视频| 日本a在线网址| 两性夫妻黄色片| 精品国产乱码久久久久久小说| 美女大奶头黄色视频| 丝袜人妻中文字幕| 人人妻人人澡人人看| 日本撒尿小便嘘嘘汇集6| e午夜精品久久久久久久| 两人在一起打扑克的视频| 美女扒开内裤让男人捅视频| 国产人伦9x9x在线观看| 制服人妻中文乱码| 国产精品欧美亚洲77777| 国产亚洲午夜精品一区二区久久| 亚洲欧美日韩高清在线视频 | 久久午夜综合久久蜜桃| 免费一级毛片在线播放高清视频 | 少妇的丰满在线观看| 岛国在线观看网站| 久久久久精品人妻al黑| 午夜久久久在线观看| 婷婷色av中文字幕| 美女大奶头黄色视频| 国产99久久九九免费精品| 精品亚洲成a人片在线观看| cao死你这个sao货| 少妇猛男粗大的猛烈进出视频| 国产亚洲精品第一综合不卡| 色精品久久人妻99蜜桃| 制服人妻中文乱码| 国产精品久久久久久精品古装| kizo精华| 中文字幕人妻丝袜制服| 中国国产av一级| 嫩草影视91久久| 天天操日日干夜夜撸| 亚洲精品久久午夜乱码| 免费少妇av软件| 香蕉国产在线看| 咕卡用的链子| 亚洲第一欧美日韩一区二区三区 | 亚洲成人免费电影在线观看| 黄频高清免费视频| 日日摸夜夜添夜夜添小说| 国产日韩欧美在线精品| 久9热在线精品视频| 老司机深夜福利视频在线观看 | 91精品国产国语对白视频| 美女中出高潮动态图| 午夜精品久久久久久毛片777| av电影中文网址| 久久中文看片网| 女人精品久久久久毛片| 韩国高清视频一区二区三区| 亚洲av欧美aⅴ国产| 窝窝影院91人妻| 国产真人三级小视频在线观看| 大香蕉久久成人网| 欧美精品啪啪一区二区三区 | 欧美中文综合在线视频| 国产日韩欧美亚洲二区| av福利片在线| 动漫黄色视频在线观看| 免费观看av网站的网址| 国产免费视频播放在线视频| 欧美老熟妇乱子伦牲交| av网站在线播放免费| 亚洲国产成人一精品久久久| 欧美日韩成人在线一区二区| 亚洲一区中文字幕在线| 天天添夜夜摸| 久久ye,这里只有精品| www.自偷自拍.com| 天堂中文最新版在线下载| www.精华液| 女人被躁到高潮嗷嗷叫费观| 99国产综合亚洲精品| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美激情在线| 日本a在线网址| 啦啦啦在线免费观看视频4| 国产真人三级小视频在线观看| tocl精华| 一级a爱视频在线免费观看| 国产成人精品在线电影| 淫妇啪啪啪对白视频 | 女警被强在线播放| 激情视频va一区二区三区| 国产一区二区三区在线臀色熟女 | 久久久国产精品麻豆| 一级毛片电影观看| 亚洲国产看品久久| 女警被强在线播放| 黄色 视频免费看| 国产福利在线免费观看视频| 午夜福利影视在线免费观看| 色94色欧美一区二区| 成人影院久久| avwww免费| 性少妇av在线| 伊人久久大香线蕉亚洲五| 国产成人影院久久av| 国产在视频线精品| 一个人免费在线观看的高清视频 | 日本欧美视频一区| 性少妇av在线| 亚洲av成人不卡在线观看播放网 | 中文字幕制服av| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人久久小说| 午夜福利免费观看在线| av国产精品久久久久影院| 久久久精品免费免费高清| e午夜精品久久久久久久| 国产一级毛片在线| 欧美亚洲 丝袜 人妻 在线| 色婷婷久久久亚洲欧美| 香蕉丝袜av| 啦啦啦在线免费观看视频4| 欧美午夜高清在线| 丰满人妻熟妇乱又伦精品不卡| 精品人妻1区二区| 亚洲色图 男人天堂 中文字幕| 欧美日韩精品网址| 99国产精品一区二区蜜桃av | 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| e午夜精品久久久久久久| 成人av一区二区三区在线看 | 狠狠精品人妻久久久久久综合| 久久久久久久久久久久大奶| 我要看黄色一级片免费的| 成人免费观看视频高清| 丁香六月欧美| 岛国毛片在线播放| 国精品久久久久久国模美| 久久亚洲国产成人精品v| 亚洲综合色网址| 亚洲少妇的诱惑av| 99久久精品国产亚洲精品| 女人被躁到高潮嗷嗷叫费观| 免费观看a级毛片全部| 黑人巨大精品欧美一区二区mp4| 男女之事视频高清在线观看| 丰满少妇做爰视频| 午夜福利在线免费观看网站| 国精品久久久久久国模美| 久久av网站| 国产精品 欧美亚洲| 天天影视国产精品| 91麻豆av在线| 美女午夜性视频免费| 最新在线观看一区二区三区| 丁香六月天网| 极品人妻少妇av视频| 999精品在线视频| 亚洲av国产av综合av卡| 亚洲欧美精品综合一区二区三区| 亚洲七黄色美女视频| 天天躁日日躁夜夜躁夜夜| 日韩熟女老妇一区二区性免费视频| 人人澡人人妻人| 久久精品国产亚洲av香蕉五月 | 精品久久久久久电影网| 日韩制服丝袜自拍偷拍| 一进一出抽搐动态| 好男人电影高清在线观看| 一区在线观看完整版| 国产精品国产三级国产专区5o| 精品福利永久在线观看| 777久久人妻少妇嫩草av网站| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美网| www.av在线官网国产| 午夜免费鲁丝| 久久中文看片网| 高清在线国产一区| 日韩欧美一区二区三区在线观看 | av超薄肉色丝袜交足视频| www日本在线高清视频| 亚洲性夜色夜夜综合| 亚洲欧洲日产国产| 自线自在国产av| 电影成人av| 啦啦啦中文免费视频观看日本| 国产亚洲av片在线观看秒播厂| 这个男人来自地球电影免费观看| 午夜福利视频精品| 五月开心婷婷网| 窝窝影院91人妻| 国产一区二区三区综合在线观看| 精品少妇久久久久久888优播| 欧美老熟妇乱子伦牲交| 青青草视频在线视频观看| 欧美激情 高清一区二区三区| 男女高潮啪啪啪动态图| 久久中文看片网| 中国国产av一级| 1024香蕉在线观看| 亚洲精品av麻豆狂野| 免费在线观看完整版高清| 久久久久精品国产欧美久久久 | 老司机午夜十八禁免费视频| 韩国精品一区二区三区| 操美女的视频在线观看| 丁香六月天网| 成年女人毛片免费观看观看9 | 亚洲专区字幕在线| 91av网站免费观看| 免费在线观看视频国产中文字幕亚洲 | 日韩中文字幕视频在线看片| 无遮挡黄片免费观看| 夜夜夜夜夜久久久久| 97在线人人人人妻| tocl精华| 91字幕亚洲| 操出白浆在线播放| av天堂在线播放| 日本撒尿小便嘘嘘汇集6| 女人高潮潮喷娇喘18禁视频| 国产精品亚洲av一区麻豆| 精品亚洲成a人片在线观看| 夜夜夜夜夜久久久久| 国产亚洲av片在线观看秒播厂| 日本vs欧美在线观看视频| 国产在线一区二区三区精| 视频区图区小说| 天天躁夜夜躁狠狠躁躁| 两人在一起打扑克的视频| 国产亚洲av片在线观看秒播厂| 国产成人免费观看mmmm| 美女午夜性视频免费| 亚洲avbb在线观看| 他把我摸到了高潮在线观看 | 久久人人爽人人片av| 丝袜人妻中文字幕| 黄色 视频免费看| 99精品久久久久人妻精品| 91字幕亚洲| 老司机午夜十八禁免费视频| 国产精品久久久av美女十八| 国产精品欧美亚洲77777| 在线永久观看黄色视频| 两个人看的免费小视频| av有码第一页| 久久精品国产综合久久久| 香蕉丝袜av| 色婷婷久久久亚洲欧美| 久久精品成人免费网站| 国产成人欧美| 欧美日韩一级在线毛片| 亚洲欧美日韩另类电影网站| 国产成人精品久久二区二区91| 美女福利国产在线| 国产精品久久久av美女十八| 别揉我奶头~嗯~啊~动态视频 | 免费女性裸体啪啪无遮挡网站| 淫妇啪啪啪对白视频 | 亚洲av欧美aⅴ国产| 日本av手机在线免费观看| 国产主播在线观看一区二区| 极品人妻少妇av视频| 黄色 视频免费看| 国产黄色免费在线视频| 免费在线观看视频国产中文字幕亚洲 | avwww免费| 最近最新免费中文字幕在线| 黑人操中国人逼视频| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 老汉色∧v一级毛片| 人成视频在线观看免费观看| 大片免费播放器 马上看| 热99久久久久精品小说推荐| 女人精品久久久久毛片| 婷婷色av中文字幕| 一级,二级,三级黄色视频| 国产亚洲精品一区二区www | www日本在线高清视频| 成人亚洲精品一区在线观看| 成人影院久久| 久久精品熟女亚洲av麻豆精品| 桃花免费在线播放| 久久亚洲精品不卡| 9色porny在线观看| 性高湖久久久久久久久免费观看| 夜夜骑夜夜射夜夜干| 欧美精品人与动牲交sv欧美| 久久久久久免费高清国产稀缺| 久9热在线精品视频| 国产在线一区二区三区精| 可以免费在线观看a视频的电影网站| 精品一品国产午夜福利视频| 人人妻人人添人人爽欧美一区卜| 亚洲人成77777在线视频| 又黄又粗又硬又大视频| 日韩人妻精品一区2区三区| 黄网站色视频无遮挡免费观看| 午夜精品久久久久久毛片777| 国产精品影院久久| 免费久久久久久久精品成人欧美视频| 热99久久久久精品小说推荐| 日日摸夜夜添夜夜添小说| 日韩欧美免费精品| 亚洲五月色婷婷综合| 大片免费播放器 马上看| 精品国产乱码久久久久久男人| 不卡av一区二区三区| 狠狠婷婷综合久久久久久88av| 中文字幕高清在线视频| 美女高潮到喷水免费观看| 精品少妇久久久久久888优播| 欧美午夜高清在线| 91精品国产国语对白视频| 久久人妻熟女aⅴ| 亚洲欧美日韩另类电影网站| bbb黄色大片| 丝袜喷水一区| 日韩中文字幕视频在线看片| 国产亚洲av片在线观看秒播厂| 国产一区二区在线观看av| 少妇粗大呻吟视频| 99精品久久久久人妻精品| 91成人精品电影| svipshipincom国产片| 日本精品一区二区三区蜜桃| 美国免费a级毛片| 色精品久久人妻99蜜桃| 亚洲专区中文字幕在线| 欧美日韩亚洲国产一区二区在线观看 | 搡老乐熟女国产| 黄色视频不卡| 在线 av 中文字幕| 后天国语完整版免费观看| 涩涩av久久男人的天堂| 搡老乐熟女国产| 久久久国产精品麻豆| 亚洲午夜精品一区,二区,三区| 国产黄频视频在线观看| 国产在线一区二区三区精| 国产激情久久老熟女| 日本黄色日本黄色录像| av不卡在线播放| 丝袜人妻中文字幕| 精品亚洲乱码少妇综合久久| 人妻久久中文字幕网| 日韩 欧美 亚洲 中文字幕| 法律面前人人平等表现在哪些方面 | 女人久久www免费人成看片| 久久99热这里只频精品6学生| 欧美激情 高清一区二区三区| 久久ye,这里只有精品| 老汉色av国产亚洲站长工具| av不卡在线播放| 国产精品偷伦视频观看了| 日韩欧美国产一区二区入口| 欧美少妇被猛烈插入视频| 日韩三级视频一区二区三区| 久久久久网色| 精品熟女少妇八av免费久了| 精品视频人人做人人爽| 日日爽夜夜爽网站| 亚洲黑人精品在线| 狠狠精品人妻久久久久久综合| 精品欧美一区二区三区在线| 十八禁网站网址无遮挡| 色婷婷久久久亚洲欧美| 美女主播在线视频| 国产精品香港三级国产av潘金莲| 久久99热这里只频精品6学生| 少妇裸体淫交视频免费看高清 | www.精华液| 一二三四社区在线视频社区8| 搡老熟女国产l中国老女人| 亚洲精华国产精华精| 999久久久国产精品视频| 日本猛色少妇xxxxx猛交久久| 亚洲,欧美精品.| 人妻人人澡人人爽人人| 亚洲人成电影观看| 无遮挡黄片免费观看| 精品一区二区三区av网在线观看 | 丝袜人妻中文字幕| 新久久久久国产一级毛片| 91大片在线观看| 国产男女超爽视频在线观看| 午夜两性在线视频| 一本—道久久a久久精品蜜桃钙片| 精品乱码久久久久久99久播| 我要看黄色一级片免费的| 香蕉国产在线看| 国产日韩欧美在线精品| 精品第一国产精品| 天堂中文最新版在线下载| 蜜桃国产av成人99| 老司机福利观看| av有码第一页| 一级黄色大片毛片| 久久久久久亚洲精品国产蜜桃av| 亚洲男人天堂网一区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产欧美日韩av| av网站在线播放免费| 在线天堂中文资源库| 精品福利观看| 极品人妻少妇av视频| 久久亚洲国产成人精品v| 中文字幕精品免费在线观看视频| 成人手机av| 巨乳人妻的诱惑在线观看| 啦啦啦在线免费观看视频4| 欧美精品人与动牲交sv欧美| 乱人伦中国视频| 久久热在线av| 日本a在线网址| 纯流量卡能插随身wifi吗| 一级毛片电影观看| 欧美日韩亚洲综合一区二区三区_| 国产日韩欧美视频二区| 日韩制服骚丝袜av| 欧美黄色片欧美黄色片| 五月天丁香电影| 操出白浆在线播放| 欧美午夜高清在线| 亚洲自偷自拍图片 自拍| 别揉我奶头~嗯~啊~动态视频 | 久久青草综合色| 永久免费av网站大全| 三上悠亚av全集在线观看| 国产片内射在线| 肉色欧美久久久久久久蜜桃| 少妇粗大呻吟视频| 精品一区在线观看国产| 免费女性裸体啪啪无遮挡网站| 免费在线观看影片大全网站| 国产精品 欧美亚洲| 叶爱在线成人免费视频播放| 热re99久久精品国产66热6| 欧美另类亚洲清纯唯美| 黑丝袜美女国产一区| 999精品在线视频| 夫妻午夜视频| 最近中文字幕2019免费版| 黄色 视频免费看| 久久久精品94久久精品| 肉色欧美久久久久久久蜜桃| 三上悠亚av全集在线观看| 国产高清国产精品国产三级| 欧美xxⅹ黑人| 丝瓜视频免费看黄片| 淫妇啪啪啪对白视频 | 亚洲精品国产av成人精品| 99香蕉大伊视频| 国产一卡二卡三卡精品| 老司机影院成人| 少妇人妻久久综合中文| 亚洲成av片中文字幕在线观看| 久久久久网色| av片东京热男人的天堂|