• <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| 最近最新中文字幕大全电影3 | 亚洲第一av免费看| 亚洲伊人色综图| 亚洲伊人色综图| 大型黄色视频在线免费观看| 欧美+亚洲+日韩+国产| 婷婷精品国产亚洲av在线| 村上凉子中文字幕在线| av在线播放免费不卡| 久久精品国产清高在天天线| 在线观看一区二区三区| 免费少妇av软件| 久久国产精品男人的天堂亚洲| 国产又爽黄色视频| 99热只有精品国产| 久久九九热精品免费| 国产成人欧美在线观看| 久久久久精品国产欧美久久久| 伦理电影免费视频| 久久午夜综合久久蜜桃| 国产精品野战在线观看 | 90打野战视频偷拍视频| 亚洲欧美一区二区三区黑人| 这个男人来自地球电影免费观看| 国产色视频综合| 日韩一卡2卡3卡4卡2021年| 夜夜夜夜夜久久久久| 欧美另类亚洲清纯唯美| 国产在线精品亚洲第一网站| 一区二区三区国产精品乱码| 欧美日韩国产mv在线观看视频| 一边摸一边抽搐一进一出视频| 一区二区三区激情视频| 黑人巨大精品欧美一区二区蜜桃| 国产又爽黄色视频| 国产成年人精品一区二区 | 国产成人精品久久二区二区免费| 久久久久国产精品人妻aⅴ院| 在线观看www视频免费| 99久久人妻综合| 男女做爰动态图高潮gif福利片 | 在线视频色国产色| 黄色女人牲交| 女性生殖器流出的白浆| 日日干狠狠操夜夜爽| 欧美成狂野欧美在线观看| 精品一区二区三区四区五区乱码| 美女扒开内裤让男人捅视频| 精品国产乱子伦一区二区三区| 香蕉久久夜色| 在线观看免费视频日本深夜| 伊人久久大香线蕉亚洲五| 在线观看一区二区三区| 亚洲精品粉嫩美女一区| 在线看a的网站| 真人一进一出gif抽搐免费| 久久久久久亚洲精品国产蜜桃av| 国产高清videossex| 窝窝影院91人妻| 日韩精品中文字幕看吧| 老司机靠b影院| 国产真人三级小视频在线观看| 中文字幕人妻丝袜一区二区| 欧美日韩亚洲综合一区二区三区_| 中亚洲国语对白在线视频| 国产极品粉嫩免费观看在线| 久久草成人影院| 精品久久久久久成人av| 两个人看的免费小视频| 亚洲午夜理论影院| 国产精品99久久99久久久不卡| 99久久久亚洲精品蜜臀av| 亚洲欧美日韩无卡精品| 亚洲男人天堂网一区| 99精国产麻豆久久婷婷| 欧美激情极品国产一区二区三区| 51午夜福利影视在线观看| a在线观看视频网站| 亚洲精华国产精华精| 在线观看日韩欧美| 老司机福利观看| 亚洲色图综合在线观看| 欧美激情久久久久久爽电影 | 亚洲自拍偷在线| 亚洲男人天堂网一区| 18禁美女被吸乳视频| 成人黄色视频免费在线看| 老司机午夜福利在线观看视频| 精品久久久久久久久久免费视频 | 国产又色又爽无遮挡免费看| 99在线人妻在线中文字幕| 亚洲av第一区精品v没综合| 99精品久久久久人妻精品| 久久久久国产一级毛片高清牌| 怎么达到女性高潮| 亚洲色图av天堂| 99国产极品粉嫩在线观看| 丰满的人妻完整版| 黑丝袜美女国产一区| 18禁裸乳无遮挡免费网站照片 | 男女午夜视频在线观看| 精品乱码久久久久久99久播| 男女之事视频高清在线观看| 久9热在线精品视频| 精品久久久久久成人av| 国产亚洲精品综合一区在线观看 | 亚洲男人天堂网一区| 亚洲精品中文字幕一二三四区| 欧美人与性动交α欧美精品济南到| 精品日产1卡2卡| 可以在线观看毛片的网站| 激情在线观看视频在线高清| 亚洲 国产 在线| a级片在线免费高清观看视频| 久久人妻福利社区极品人妻图片| 一区二区三区精品91| 精品欧美一区二区三区在线| 日日夜夜操网爽| 色婷婷久久久亚洲欧美| 999久久久精品免费观看国产| 亚洲国产精品999在线| 黄色视频不卡| 久久天躁狠狠躁夜夜2o2o| 十分钟在线观看高清视频www| 成人三级做爰电影| 国产1区2区3区精品| 九色亚洲精品在线播放| av中文乱码字幕在线| 在线免费观看的www视频| 久久99一区二区三区| 精品国产一区二区久久| 欧美大码av| 曰老女人黄片| 免费看a级黄色片| 国产欧美日韩精品亚洲av| 国产精品久久久久久人妻精品电影| 性少妇av在线| 亚洲色图综合在线观看| 一区二区三区激情视频| 精品一品国产午夜福利视频| 叶爱在线成人免费视频播放| 69精品国产乱码久久久| 人妻久久中文字幕网| 亚洲欧美一区二区三区久久| 成人av一区二区三区在线看| 他把我摸到了高潮在线观看| 久久精品国产亚洲av高清一级| 久久这里只有精品19| 亚洲在线自拍视频| av免费在线观看网站| 免费看a级黄色片| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片 | 色在线成人网| 久久久国产精品麻豆| 国产欧美日韩一区二区精品| 伦理电影免费视频| 老司机亚洲免费影院| 久久欧美精品欧美久久欧美| www.999成人在线观看| 亚洲五月天丁香| 中文字幕色久视频| 一区二区三区国产精品乱码| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品一二三| 精品福利永久在线观看| 在线av久久热| 91在线观看av| 少妇粗大呻吟视频| 1024香蕉在线观看| 母亲3免费完整高清在线观看| 在线av久久热| 色播在线永久视频| 老熟妇仑乱视频hdxx| 欧美黄色片欧美黄色片| 一边摸一边抽搐一进一小说| 欧美中文日本在线观看视频| 亚洲国产欧美网| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜一区二区| 大香蕉久久成人网| 成人18禁在线播放| 欧美精品亚洲一区二区| 在线观看一区二区三区激情| av在线天堂中文字幕 | 精品久久久久久久久久免费视频 | 纯流量卡能插随身wifi吗| 国产xxxxx性猛交| 69av精品久久久久久| 人人妻人人澡人人看| 国产日韩一区二区三区精品不卡| 少妇 在线观看| 日韩 欧美 亚洲 中文字幕| 性色av乱码一区二区三区2| 亚洲成国产人片在线观看| 一本综合久久免费| 婷婷精品国产亚洲av在线| 女性生殖器流出的白浆| 在线观看免费午夜福利视频| 一级毛片女人18水好多| 精品一区二区三区视频在线观看免费 | 久久久久久人人人人人| 亚洲精品美女久久av网站| 在线观看免费日韩欧美大片| 在线播放国产精品三级| 国产精品99久久99久久久不卡| 久久精品aⅴ一区二区三区四区| 亚洲国产精品一区二区三区在线| 欧美丝袜亚洲另类 | 亚洲精华国产精华精| 午夜久久久在线观看| 一边摸一边抽搐一进一出视频| 国产xxxxx性猛交| 日韩视频一区二区在线观看| 老司机午夜十八禁免费视频| 亚洲色图综合在线观看| 麻豆久久精品国产亚洲av | 亚洲色图综合在线观看| 动漫黄色视频在线观看| 国产一卡二卡三卡精品| 超色免费av| 久久午夜综合久久蜜桃| 精品国产超薄肉色丝袜足j| 国产精品免费一区二区三区在线| 嫁个100分男人电影在线观看| 久久人妻av系列| 精品国产亚洲在线| 99热只有精品国产| 午夜精品在线福利| 国产99久久九九免费精品| 少妇粗大呻吟视频| 亚洲中文av在线| 在线观看免费午夜福利视频| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩精品网址| 在线观看一区二区三区| 国产精品 国内视频| 午夜成年电影在线免费观看| 老司机在亚洲福利影院| 亚洲片人在线观看| 麻豆av在线久日| 成人18禁高潮啪啪吃奶动态图| 亚洲人成网站在线播放欧美日韩| 免费av毛片视频| 1024视频免费在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美精品综合一区二区三区| 一区福利在线观看| 亚洲国产欧美一区二区综合| 亚洲男人的天堂狠狠| 久热爱精品视频在线9| 亚洲激情在线av| 日韩人妻精品一区2区三区| 国产片内射在线| 午夜两性在线视频| 欧美日本中文国产一区发布| 欧美日韩福利视频一区二区| 国产片内射在线| 国产精品国产高清国产av| 精品国产一区二区三区四区第35| 国产成人精品无人区| 日本三级黄在线观看| 精品少妇一区二区三区视频日本电影| 一进一出好大好爽视频| 欧美成狂野欧美在线观看| 国产成人av激情在线播放| ponron亚洲| 国产一区二区在线av高清观看| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 亚洲国产毛片av蜜桃av| 欧美日韩乱码在线| 亚洲国产中文字幕在线视频| 黑人巨大精品欧美一区二区mp4| 人人妻人人爽人人添夜夜欢视频| 国产熟女午夜一区二区三区| 麻豆一二三区av精品| 99久久人妻综合| 免费女性裸体啪啪无遮挡网站| 成人三级做爰电影| 变态另类成人亚洲欧美熟女 | 色综合站精品国产| 女人被狂操c到高潮| 国产乱人伦免费视频| 欧美激情 高清一区二区三区| 免费在线观看黄色视频的| 午夜日韩欧美国产| 日韩成人在线观看一区二区三区| 亚洲黑人精品在线| 夫妻午夜视频| 精品国产超薄肉色丝袜足j| 久久久国产精品麻豆| 久久草成人影院| 久久国产精品影院| 1024香蕉在线观看| 久9热在线精品视频| 韩国av一区二区三区四区| 巨乳人妻的诱惑在线观看| 又大又爽又粗| 亚洲精品av麻豆狂野| 亚洲美女黄片视频| 欧美在线一区亚洲| 黄色a级毛片大全视频| 午夜精品国产一区二区电影| 亚洲九九香蕉| svipshipincom国产片| 久久婷婷成人综合色麻豆| 嫩草影院精品99| 亚洲男人的天堂狠狠| 国产一区二区三区视频了| 亚洲av电影在线进入| 精品久久久久久成人av| 丰满饥渴人妻一区二区三| 国产免费现黄频在线看| 免费在线观看亚洲国产| 日韩免费高清中文字幕av| 日本 av在线| 成人黄色视频免费在线看| 亚洲一区二区三区欧美精品| 怎么达到女性高潮| 亚洲成人免费av在线播放| 老鸭窝网址在线观看| 国产在线精品亚洲第一网站| 亚洲av成人不卡在线观看播放网| 亚洲av第一区精品v没综合| 正在播放国产对白刺激| 亚洲 国产 在线| 亚洲一码二码三码区别大吗| 国产精品一区二区免费欧美| 亚洲伊人色综图| 在线观看午夜福利视频| 男男h啪啪无遮挡| 80岁老熟妇乱子伦牲交| 69精品国产乱码久久久| 国产精品乱码一区二三区的特点 | 午夜福利欧美成人| 1024视频免费在线观看| 丰满人妻熟妇乱又伦精品不卡| 悠悠久久av| 亚洲av电影在线进入| 中出人妻视频一区二区| 国产深夜福利视频在线观看| 久久午夜亚洲精品久久| 久久久久国产精品人妻aⅴ院| 亚洲久久久国产精品| 成人18禁高潮啪啪吃奶动态图| 在线永久观看黄色视频| 校园春色视频在线观看| 日韩精品免费视频一区二区三区| 国产精品 欧美亚洲| 精品一品国产午夜福利视频| 国产精品永久免费网站| 亚洲成人国产一区在线观看| 精品久久久久久久久久免费视频 | 国产成人啪精品午夜网站| 久久欧美精品欧美久久欧美| 最近最新免费中文字幕在线| 脱女人内裤的视频| 免费看a级黄色片| 久久久久九九精品影院| 成人18禁高潮啪啪吃奶动态图| 欧美+亚洲+日韩+国产| 麻豆av在线久日| 看黄色毛片网站| 久久久水蜜桃国产精品网| 国产在线精品亚洲第一网站| 久久精品影院6| 看免费av毛片| videosex国产| 欧美乱码精品一区二区三区| 亚洲av第一区精品v没综合| 久久精品亚洲熟妇少妇任你| 亚洲黑人精品在线| 777久久人妻少妇嫩草av网站| 亚洲黑人精品在线| 国产亚洲av高清不卡| 18禁裸乳无遮挡免费网站照片 | 午夜免费激情av| 亚洲欧洲精品一区二区精品久久久| 纯流量卡能插随身wifi吗| 亚洲七黄色美女视频| 亚洲 欧美一区二区三区| 亚洲第一av免费看| 欧美日韩福利视频一区二区| 欧美激情久久久久久爽电影 | 丝袜美腿诱惑在线| 国产精品一区二区在线不卡| 国产三级在线视频| 午夜精品在线福利| 欧美日韩精品网址| 亚洲性夜色夜夜综合| 欧美大码av| 欧美av亚洲av综合av国产av| 欧美精品一区二区免费开放| 日韩精品中文字幕看吧| 国产免费现黄频在线看| 天堂影院成人在线观看| 亚洲自拍偷在线| 69精品国产乱码久久久| 欧美日韩精品网址| 久久精品人人爽人人爽视色| 级片在线观看| 久久国产精品男人的天堂亚洲| 久99久视频精品免费| 女同久久另类99精品国产91| 久久久久久大精品| 级片在线观看| 国产成人系列免费观看| 别揉我奶头~嗯~啊~动态视频| 国产乱人伦免费视频| 亚洲人成77777在线视频| 久久人妻av系列| 亚洲国产精品999在线| 国产极品粉嫩免费观看在线| 亚洲专区国产一区二区| 性色av乱码一区二区三区2| 日日摸夜夜添夜夜添小说| 亚洲精品一二三| 欧美一级毛片孕妇| 国产精品综合久久久久久久免费 | 亚洲熟妇熟女久久| 黄色丝袜av网址大全| 免费少妇av软件| 国产成人一区二区三区免费视频网站| 一区二区三区激情视频| 日本a在线网址| 看片在线看免费视频| 国产精品免费一区二区三区在线| 大码成人一级视频| 嫩草影视91久久| 日韩高清综合在线| 久久久久国产精品人妻aⅴ院| 成年人免费黄色播放视频| 午夜老司机福利片| 不卡一级毛片| 狠狠狠狠99中文字幕| 精品高清国产在线一区| 午夜福利,免费看| 欧美另类亚洲清纯唯美| 国产亚洲精品第一综合不卡| 在线观看日韩欧美| 亚洲 欧美 日韩 在线 免费| cao死你这个sao货| 亚洲专区字幕在线| 国产三级黄色录像| 淫妇啪啪啪对白视频| 丰满的人妻完整版| 色综合欧美亚洲国产小说| 97碰自拍视频| 精品熟女少妇八av免费久了| 精品国产亚洲在线| 亚洲精华国产精华精| 又黄又粗又硬又大视频| 亚洲av电影在线进入| 50天的宝宝边吃奶边哭怎么回事| 久热这里只有精品99| 91麻豆av在线| 国产三级在线视频| 美女 人体艺术 gogo| 97碰自拍视频| 久久这里只有精品19| 男人舔女人的私密视频| 国产麻豆69| 水蜜桃什么品种好| 超色免费av| 亚洲 国产 在线| 成年女人毛片免费观看观看9| 国产精品久久久人人做人人爽| xxxhd国产人妻xxx| 中文字幕人妻丝袜制服| 国产亚洲精品久久久久久毛片| 欧美日韩福利视频一区二区| 搡老岳熟女国产| 国产高清国产精品国产三级| 90打野战视频偷拍视频| 亚洲国产精品合色在线| 欧美成狂野欧美在线观看| 亚洲精品在线观看二区| 亚洲国产中文字幕在线视频| 99在线人妻在线中文字幕| 国产伦一二天堂av在线观看| 丝袜人妻中文字幕| 波多野结衣高清无吗| 熟女少妇亚洲综合色aaa.| 国产精品二区激情视频| 一本综合久久免费| 高清毛片免费观看视频网站 | 又大又爽又粗| 久久欧美精品欧美久久欧美| 老司机午夜福利在线观看视频| 黄色女人牲交| 一级毛片高清免费大全| 90打野战视频偷拍视频| 亚洲国产欧美网| 午夜福利在线免费观看网站| 日日摸夜夜添夜夜添小说| 长腿黑丝高跟| 免费观看精品视频网站| 亚洲久久久国产精品| 午夜成年电影在线免费观看| 国产精品久久久av美女十八| 女人精品久久久久毛片| 婷婷精品国产亚洲av在线| 很黄的视频免费| 99久久人妻综合| 精品一区二区三区四区五区乱码| 国产欧美日韩一区二区三区在线| 波多野结衣一区麻豆| 999精品在线视频| 成人黄色视频免费在线看| 99久久久亚洲精品蜜臀av| 国产欧美日韩一区二区精品| 999精品在线视频| 久9热在线精品视频| 亚洲成人久久性| 欧美日韩福利视频一区二区| 亚洲 国产 在线| 精品久久久久久电影网| 色哟哟哟哟哟哟| 亚洲色图 男人天堂 中文字幕| 欧美日韩黄片免| 一区在线观看完整版| 日本免费一区二区三区高清不卡 | 国产区一区二久久| 亚洲成人精品中文字幕电影 | 亚洲色图av天堂| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站在线播放欧美日韩| 久久久久九九精品影院| 丰满饥渴人妻一区二区三| 久久九九热精品免费| 免费在线观看影片大全网站| 欧美日韩一级在线毛片| 黑人巨大精品欧美一区二区蜜桃| 久久久精品国产亚洲av高清涩受| 三上悠亚av全集在线观看| 热99re8久久精品国产| 日韩欧美国产一区二区入口| 亚洲欧美激情综合另类| 亚洲国产精品999在线| 精品日产1卡2卡| 亚洲片人在线观看| 黄色丝袜av网址大全| 99热只有精品国产| 黄色怎么调成土黄色| 国产三级黄色录像| 色综合婷婷激情| 极品教师在线免费播放| 老司机深夜福利视频在线观看| 亚洲欧美一区二区三区黑人| 亚洲av五月六月丁香网| 亚洲欧美一区二区三区黑人| 美国免费a级毛片| 久久久久久大精品| 久久 成人 亚洲| 亚洲熟妇中文字幕五十中出 | av中文乱码字幕在线| 操出白浆在线播放| 色哟哟哟哟哟哟| 18禁裸乳无遮挡免费网站照片 | 成人国产一区最新在线观看| 久久九九热精品免费| 侵犯人妻中文字幕一二三四区| 精品福利永久在线观看| 欧美乱色亚洲激情| 大陆偷拍与自拍| 怎么达到女性高潮| 亚洲久久久国产精品| 久久人人爽av亚洲精品天堂| 国产一卡二卡三卡精品| 国产精品 欧美亚洲| 亚洲国产看品久久| 中出人妻视频一区二区| 亚洲国产精品999在线| 亚洲成国产人片在线观看| 新久久久久国产一级毛片| 日韩一卡2卡3卡4卡2021年| 日韩欧美国产一区二区入口| 日本撒尿小便嘘嘘汇集6| 精品久久久久久,| 可以免费在线观看a视频的电影网站| 午夜影院日韩av| 狂野欧美激情性xxxx| 精品久久久久久,| 欧美一级毛片孕妇| 亚洲第一av免费看| 老司机午夜福利在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 在线av久久热| 精品久久久久久,| 欧美中文日本在线观看视频| 亚洲国产毛片av蜜桃av| e午夜精品久久久久久久| 精品人妻1区二区| aaaaa片日本免费| 国产深夜福利视频在线观看| 久久久国产成人免费| 黄色 视频免费看| 亚洲精华国产精华精| 亚洲一区二区三区不卡视频| 91成年电影在线观看| 亚洲中文日韩欧美视频| 亚洲中文av在线| 99精品在免费线老司机午夜| 国产一卡二卡三卡精品| 多毛熟女@视频| 国产成年人精品一区二区 | 黄频高清免费视频| 丰满饥渴人妻一区二区三| 国产亚洲欧美98| 一进一出好大好爽视频| 亚洲av片天天在线观看| 国产成人影院久久av|