• <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.

    99久久精品国产亚洲精品| 精品久久久久久久末码| 热99re8久久精品国产| 午夜久久久久精精品| 欧美精品啪啪一区二区三区| 婷婷精品国产亚洲av在线| 欧美午夜高清在线| 日本一二三区视频观看| 欧美3d第一页| 久久伊人香网站| 少妇的逼水好多| 国产在线精品亚洲第一网站| 欧美精品啪啪一区二区三区| 99久久99久久久精品蜜桃| 亚洲无线观看免费| 国产精品98久久久久久宅男小说| 动漫黄色视频在线观看| 亚洲一区二区三区色噜噜| 五月伊人婷婷丁香| 亚洲成人中文字幕在线播放| 国产欧美日韩一区二区三| avwww免费| 五月伊人婷婷丁香| 搡女人真爽免费视频火全软件 | 我的女老师完整版在线观看| 欧美黑人欧美精品刺激| 国内揄拍国产精品人妻在线| 热99re8久久精品国产| 一卡2卡三卡四卡精品乱码亚洲| 国产在线精品亚洲第一网站| 老司机深夜福利视频在线观看| 九九在线视频观看精品| 2021天堂中文幕一二区在线观| 日韩欧美三级三区| 69av精品久久久久久| 成人av在线播放网站| 中出人妻视频一区二区| 成人精品一区二区免费| 99久久九九国产精品国产免费| 1000部很黄的大片| 在线观看一区二区三区| 国产成人福利小说| 一级毛片久久久久久久久女| 久久精品国产清高在天天线| 美女免费视频网站| 亚州av有码| 超碰av人人做人人爽久久| 欧美成狂野欧美在线观看| 老司机午夜福利在线观看视频| av黄色大香蕉| 在线播放无遮挡| 国产精品1区2区在线观看.| 日本免费一区二区三区高清不卡| 久久精品国产99精品国产亚洲性色| 亚州av有码| 欧美黄色片欧美黄色片| 久久久久久久午夜电影| 日日干狠狠操夜夜爽| 免费观看的影片在线观看| 51国产日韩欧美| 国产精品久久电影中文字幕| 99久久成人亚洲精品观看| 国产精品98久久久久久宅男小说| 五月玫瑰六月丁香| 99精品久久久久人妻精品| 免费在线观看日本一区| 国产精品影院久久| 女人十人毛片免费观看3o分钟| 久久草成人影院| 内射极品少妇av片p| 久久人人精品亚洲av| 黄色视频,在线免费观看| 在线看三级毛片| 在线观看舔阴道视频| 精品人妻1区二区| 黄色视频,在线免费观看| 最新在线观看一区二区三区| 精品福利观看| 长腿黑丝高跟| 国产亚洲欧美98| 一本综合久久免费| 午夜a级毛片| 亚洲精品影视一区二区三区av| 国产精品永久免费网站| 亚洲 欧美 日韩 在线 免费| 欧美精品国产亚洲| 日本撒尿小便嘘嘘汇集6| 国产午夜精品论理片| 无遮挡黄片免费观看| 搡老妇女老女人老熟妇| 男人狂女人下面高潮的视频| 在线a可以看的网站| 丰满的人妻完整版| 日韩有码中文字幕| 老司机深夜福利视频在线观看| 婷婷色综合大香蕉| 色精品久久人妻99蜜桃| 国产探花在线观看一区二区| 此物有八面人人有两片| 我要搜黄色片| 午夜久久久久精精品| 精品一区二区三区视频在线观看免费| 九九热线精品视视频播放| 欧美色欧美亚洲另类二区| 亚洲,欧美,日韩| 97人妻精品一区二区三区麻豆| 欧美成人性av电影在线观看| 十八禁国产超污无遮挡网站| 97碰自拍视频| 色精品久久人妻99蜜桃| 免费观看精品视频网站| 亚洲欧美日韩无卡精品| 性色avwww在线观看| 国内精品久久久久精免费| 91字幕亚洲| 免费大片18禁| 国产综合懂色| 村上凉子中文字幕在线| 成人午夜高清在线视频| 怎么达到女性高潮| 一级a爱片免费观看的视频| 99久久无色码亚洲精品果冻| 欧美潮喷喷水| 亚洲精品成人久久久久久| 国产精品久久久久久人妻精品电影| 高清在线国产一区| 亚洲精品在线美女| 欧美性感艳星| 老女人水多毛片| 人妻夜夜爽99麻豆av| 淫秽高清视频在线观看| 又爽又黄a免费视频| 欧美三级亚洲精品| 嫩草影院新地址| 亚洲不卡免费看| 精品国产亚洲在线| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| av在线蜜桃| 久久久久久大精品| www.熟女人妻精品国产| 人人妻,人人澡人人爽秒播| 亚洲av第一区精品v没综合| 欧美国产日韩亚洲一区| 中文字幕av成人在线电影| 在线观看美女被高潮喷水网站 | 少妇高潮的动态图| 九色成人免费人妻av| 99久久精品国产亚洲精品| 老熟妇乱子伦视频在线观看| 色综合站精品国产| 国产主播在线观看一区二区| 日日干狠狠操夜夜爽| 精品人妻一区二区三区麻豆 | 亚洲av中文字字幕乱码综合| 少妇丰满av| 青草久久国产| 亚洲欧美日韩东京热| 国产免费av片在线观看野外av| 欧美成狂野欧美在线观看| 亚洲国产日韩欧美精品在线观看| 一进一出好大好爽视频| 18禁黄网站禁片午夜丰满| 日本熟妇午夜| 亚洲精品粉嫩美女一区| 亚洲avbb在线观看| 91字幕亚洲| 狠狠狠狠99中文字幕| 国产aⅴ精品一区二区三区波| 国产爱豆传媒在线观看| 亚洲精华国产精华精| 99精品在免费线老司机午夜| 热99在线观看视频| 欧美xxxx性猛交bbbb| 亚洲va日本ⅴa欧美va伊人久久| 亚洲真实伦在线观看| 久久6这里有精品| 久久久久久久亚洲中文字幕 | 国产精品永久免费网站| 亚洲美女黄片视频| 一本一本综合久久| 国产白丝娇喘喷水9色精品| 欧美精品国产亚洲| 久久精品国产亚洲av天美| 麻豆久久精品国产亚洲av| 免费观看人在逋| 黄色丝袜av网址大全| 日韩人妻高清精品专区| 成人无遮挡网站| 99久久精品一区二区三区| 丁香欧美五月| 国产精品美女特级片免费视频播放器| 又紧又爽又黄一区二区| 有码 亚洲区| 婷婷丁香在线五月| 大型黄色视频在线免费观看| 18禁裸乳无遮挡免费网站照片| 国产探花在线观看一区二区| 草草在线视频免费看| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩东京热| 亚洲熟妇中文字幕五十中出| 亚洲avbb在线观看| 欧美日本亚洲视频在线播放| 成人无遮挡网站| 在线播放无遮挡| 高清日韩中文字幕在线| 午夜福利成人在线免费观看| 人人妻,人人澡人人爽秒播| 欧美xxxx黑人xx丫x性爽| 99热精品在线国产| 我的女老师完整版在线观看| 中文字幕免费在线视频6| 国产精品一区二区性色av| 一个人免费在线观看电影| 亚洲欧美日韩高清在线视频| 久久亚洲精品不卡| 色视频www国产| av福利片在线观看| а√天堂www在线а√下载| 亚洲三级黄色毛片| 国产免费av片在线观看野外av| 国内精品久久久久久久电影| 亚洲一区二区三区色噜噜| 国产高清有码在线观看视频| 99久久99久久久精品蜜桃| 国产私拍福利视频在线观看| 精品99又大又爽又粗少妇毛片 | 九九热线精品视视频播放| 在线播放国产精品三级| 精品久久久久久久久亚洲 | 成人美女网站在线观看视频| 午夜福利高清视频| 欧美一区二区国产精品久久精品| 日本免费a在线| 好男人在线观看高清免费视频| 波多野结衣巨乳人妻| 日韩成人在线观看一区二区三区| 免费观看人在逋| 精品午夜福利视频在线观看一区| 国产白丝娇喘喷水9色精品| 亚洲欧美精品综合久久99| 午夜精品在线福利| 亚洲欧美日韩高清在线视频| 国产精品久久久久久精品电影| 国产黄色小视频在线观看| 久久久国产成人精品二区| 国产探花极品一区二区| 国产亚洲精品久久久久久毛片| 久久久久久久久大av| 午夜亚洲福利在线播放| 欧美极品一区二区三区四区| 在线观看66精品国产| 欧美黄色片欧美黄色片| 精品熟女少妇八av免费久了| 小蜜桃在线观看免费完整版高清| 免费看a级黄色片| 一夜夜www| 色综合欧美亚洲国产小说| 白带黄色成豆腐渣| 久久久久久久久大av| 一级作爱视频免费观看| 搡老熟女国产l中国老女人| 午夜福利欧美成人| 国模一区二区三区四区视频| 久久久久久久亚洲中文字幕 | 精品欧美国产一区二区三| 中文字幕av在线有码专区| 最好的美女福利视频网| 51午夜福利影视在线观看| 国产精品久久电影中文字幕| 永久网站在线| 成人特级av手机在线观看| 中出人妻视频一区二区| 亚洲人成网站在线播放欧美日韩| 88av欧美| 99热6这里只有精品| 日本在线视频免费播放| av女优亚洲男人天堂| 18+在线观看网站| 午夜福利欧美成人| 国产欧美日韩一区二区精品| 日韩欧美在线乱码| 久久国产乱子免费精品| 欧美成人性av电影在线观看| 色综合欧美亚洲国产小说| 如何舔出高潮| 国产精品永久免费网站| 757午夜福利合集在线观看| 国产三级在线视频| 在线观看美女被高潮喷水网站 | aaaaa片日本免费| 一个人看视频在线观看www免费| 亚洲精品色激情综合| 国产日本99.免费观看| 国产精品一区二区性色av| 舔av片在线| 99久国产av精品| 夜夜看夜夜爽夜夜摸| 亚洲av日韩精品久久久久久密| 日本撒尿小便嘘嘘汇集6| 国内久久婷婷六月综合欲色啪| 亚洲午夜理论影院| 亚洲精品在线美女| 五月伊人婷婷丁香| 欧美日韩中文字幕国产精品一区二区三区| 日本 欧美在线| 亚洲男人的天堂狠狠| 精品人妻偷拍中文字幕| 国产精品久久久久久人妻精品电影| 国产高清三级在线| 精品人妻偷拍中文字幕| 欧美zozozo另类| 老司机福利观看| 久久久久久久久久成人| 久久亚洲精品不卡| 日本一二三区视频观看| 中文字幕高清在线视频| 午夜福利18| 97热精品久久久久久| 亚洲欧美日韩高清在线视频| 午夜影院日韩av| 亚洲五月天丁香| av在线老鸭窝| 国产精品久久久久久亚洲av鲁大| 国产黄色小视频在线观看| 乱码一卡2卡4卡精品| 精品人妻1区二区| 一区二区三区免费毛片| 午夜免费男女啪啪视频观看 | 国产久久久一区二区三区| 18禁黄网站禁片免费观看直播| 国产视频内射| 村上凉子中文字幕在线| 一级av片app| 最好的美女福利视频网| 女人被狂操c到高潮| 婷婷精品国产亚洲av在线| 51国产日韩欧美| 午夜福利高清视频| 美女 人体艺术 gogo| 国产免费一级a男人的天堂| 99riav亚洲国产免费| 一个人看的www免费观看视频| 99国产精品一区二区蜜桃av| 国产精品国产高清国产av| а√天堂www在线а√下载| 国产精品99久久久久久久久| 亚洲av第一区精品v没综合| 欧美激情国产日韩精品一区| 天美传媒精品一区二区| 欧美激情久久久久久爽电影| 丝袜美腿在线中文| 国产白丝娇喘喷水9色精品| 嫩草影院新地址| 欧美一级a爱片免费观看看| www日本黄色视频网| 久久天躁狠狠躁夜夜2o2o| 精品99又大又爽又粗少妇毛片 | 国产乱人视频| 欧美成人一区二区免费高清观看| 91狼人影院| 美女黄网站色视频| 男人狂女人下面高潮的视频| 一区二区三区高清视频在线| 日本a在线网址| 色综合婷婷激情| 国产成年人精品一区二区| 午夜福利成人在线免费观看| 亚洲真实伦在线观看| 精品一区二区三区视频在线观看免费| 国产一区二区激情短视频| 老司机福利观看| 日韩人妻高清精品专区| 亚洲精品久久国产高清桃花| 脱女人内裤的视频| 久久婷婷人人爽人人干人人爱| 精品午夜福利视频在线观看一区| 日韩精品青青久久久久久| 桃红色精品国产亚洲av| 免费无遮挡裸体视频| 久久久久久大精品| 怎么达到女性高潮| 香蕉av资源在线| 国内精品久久久久久久电影| 不卡一级毛片| 免费大片18禁| 国产 一区 欧美 日韩| 国产野战对白在线观看| 欧美成狂野欧美在线观看| 波多野结衣巨乳人妻| 男女做爰动态图高潮gif福利片| 一进一出抽搐动态| 婷婷精品国产亚洲av| 亚洲av日韩精品久久久久久密| 国内久久婷婷六月综合欲色啪| 自拍偷自拍亚洲精品老妇| 精品欧美国产一区二区三| 亚洲av五月六月丁香网| 亚洲最大成人av| 久久久久久久亚洲中文字幕 | aaaaa片日本免费| 国产一区二区三区视频了| 丝袜美腿在线中文| 91九色精品人成在线观看| 久久久久久久久久黄片| 国产成人啪精品午夜网站| 波多野结衣高清无吗| 男人和女人高潮做爰伦理| 网址你懂的国产日韩在线| 一个人免费在线观看电影| 欧美高清成人免费视频www| 丰满人妻一区二区三区视频av| 亚洲天堂国产精品一区在线| 欧美黑人欧美精品刺激| 国产精品野战在线观看| 欧美性感艳星| 国产乱人伦免费视频| 99riav亚洲国产免费| 国产av麻豆久久久久久久| 深夜精品福利| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线美女| 亚洲欧美日韩卡通动漫| 人妻丰满熟妇av一区二区三区| 日本五十路高清| 非洲黑人性xxxx精品又粗又长| 免费在线观看影片大全网站| 毛片一级片免费看久久久久 | 久久久色成人| 国产精品一区二区三区四区免费观看 | 最近最新免费中文字幕在线| 久久久久国产精品人妻aⅴ院| 免费看光身美女| 看十八女毛片水多多多| 日本精品一区二区三区蜜桃| 在线看三级毛片| 淫妇啪啪啪对白视频| 久久久久久久久中文| 国产野战对白在线观看| 欧美成人a在线观看| 内射极品少妇av片p| 国产亚洲欧美在线一区二区| 女人被狂操c到高潮| 国内少妇人妻偷人精品xxx网站| 97热精品久久久久久| 免费看光身美女| 最近视频中文字幕2019在线8| 精品一区二区三区av网在线观看| av欧美777| 久久草成人影院| 亚洲精品在线观看二区| av女优亚洲男人天堂| 亚洲国产高清在线一区二区三| 亚洲性夜色夜夜综合| av专区在线播放| 亚洲18禁久久av| 国产在视频线在精品| 最新中文字幕久久久久| 久久久国产成人精品二区| 男人舔奶头视频| 精品一区二区三区视频在线| 嫩草影视91久久| 极品教师在线视频| 高清日韩中文字幕在线| 午夜视频国产福利| 91在线精品国自产拍蜜月| 国产高清激情床上av| 99久久无色码亚洲精品果冻| 性色avwww在线观看| 亚洲,欧美,日韩| 麻豆国产97在线/欧美| 中文在线观看免费www的网站| 色综合亚洲欧美另类图片| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av香蕉五月| 老司机福利观看| 波多野结衣高清无吗| 久99久视频精品免费| 日韩人妻高清精品专区| 别揉我奶头~嗯~啊~动态视频| 一区二区三区激情视频| 亚洲精品影视一区二区三区av| 欧美激情久久久久久爽电影| 色播亚洲综合网| 国语自产精品视频在线第100页| 51午夜福利影视在线观看| 熟妇人妻久久中文字幕3abv| 亚洲av一区综合| 搡老熟女国产l中国老女人| 久久精品综合一区二区三区| 国产成人a区在线观看| 日本 av在线| 99热6这里只有精品| 免费在线观看影片大全网站| 午夜老司机福利剧场| 黄色配什么色好看| 国产伦人伦偷精品视频| 欧美高清成人免费视频www| 色综合站精品国产| 国产伦精品一区二区三区四那| 亚洲国产欧洲综合997久久,| 亚洲人成网站高清观看| 亚州av有码| 亚洲在线自拍视频| 精品午夜福利在线看| 老女人水多毛片| 97超视频在线观看视频| 又黄又爽又刺激的免费视频.| 97超视频在线观看视频| 国产色爽女视频免费观看| 亚洲激情在线av| 99久久成人亚洲精品观看| 草草在线视频免费看| 国产精品久久电影中文字幕| 久久精品久久久久久噜噜老黄 | 亚洲欧美清纯卡通| 国产男靠女视频免费网站| 国产av一区在线观看免费| 琪琪午夜伦伦电影理论片6080| 黄色一级大片看看| 亚洲av免费高清在线观看| 国产午夜精品久久久久久一区二区三区 | 久久精品国产亚洲av香蕉五月| 可以在线观看毛片的网站| 69人妻影院| 亚洲一区二区三区色噜噜| 国产成人欧美在线观看| 深夜精品福利| 伊人久久精品亚洲午夜| 久久久色成人| 3wmmmm亚洲av在线观看| 窝窝影院91人妻| 69av精品久久久久久| 久久久久久九九精品二区国产| 欧美极品一区二区三区四区| 色哟哟·www| 人妻制服诱惑在线中文字幕| 在现免费观看毛片| 成人三级黄色视频| 日韩大尺度精品在线看网址| 男女之事视频高清在线观看| 国产免费一级a男人的天堂| 91麻豆av在线| 有码 亚洲区| 国产黄片美女视频| 热99re8久久精品国产| 国产蜜桃级精品一区二区三区| 人人妻人人看人人澡| 久久久久国内视频| 国产野战对白在线观看| 国产成人a区在线观看| 亚洲av第一区精品v没综合| 1000部很黄的大片| 一卡2卡三卡四卡精品乱码亚洲| 国产国拍精品亚洲av在线观看| 又黄又爽又免费观看的视频| 亚洲成人久久爱视频| 久久精品国产亚洲av香蕉五月| 色综合婷婷激情| 国产成人影院久久av| 久久国产乱子伦精品免费另类| 亚洲在线观看片| av视频在线观看入口| 精品人妻偷拍中文字幕| 99久久成人亚洲精品观看| 欧美性猛交╳xxx乱大交人| 成人特级黄色片久久久久久久| 一级av片app| 免费高清视频大片| 黄色配什么色好看| 午夜激情欧美在线| 国产精品一及| 久久久久久久久久成人| 精品一区二区三区视频在线| 久久亚洲精品不卡| 午夜日韩欧美国产| 桃红色精品国产亚洲av| 婷婷色综合大香蕉| 成人欧美大片| 在线观看午夜福利视频| 深夜a级毛片| 美女被艹到高潮喷水动态| 久久天躁狠狠躁夜夜2o2o| 最好的美女福利视频网| 国内精品久久久久精免费| 啪啪无遮挡十八禁网站| 国产av在哪里看| 99热这里只有是精品50| 一区二区三区四区激情视频 | 精品久久久久久成人av| 国产精品爽爽va在线观看网站| 夜夜看夜夜爽夜夜摸| 婷婷丁香在线五月| 国产一区二区激情短视频| 精品不卡国产一区二区三区| 国产精品伦人一区二区| 琪琪午夜伦伦电影理论片6080| 国产亚洲精品久久久久久毛片| 亚洲第一欧美日韩一区二区三区| 久久久久久国产a免费观看| 亚洲人成网站在线播放欧美日韩| 中出人妻视频一区二区| 亚洲av日韩精品久久久久久密| 99国产精品一区二区三区| 欧美成人a在线观看| 亚洲中文字幕一区二区三区有码在线看| 精品一区二区三区视频在线| 精品人妻偷拍中文字幕| 在线免费观看不下载黄p国产 | 久久久成人免费电影| 国产精品电影一区二区三区| 欧美成人性av电影在线观看| 色哟哟哟哟哟哟| 蜜桃亚洲精品一区二区三区|