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

    基于傳遞率函數(shù)的運(yùn)行模態(tài)分析方法*

    2015-03-13 02:30:38張永年夏遵平
    振動、測試與診斷 2015年5期
    關(guān)鍵詞:機(jī)翼阻尼模態(tài)

    張永年, 王 彤, 夏遵平

    (南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國家重點(diǎn)實驗室 南京,210016)

    ?

    基于傳遞率函數(shù)的運(yùn)行模態(tài)分析方法*

    張永年, 王 彤, 夏遵平

    (南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國家重點(diǎn)實驗室 南京,210016)

    傳統(tǒng)的運(yùn)行模態(tài)分析方法在推導(dǎo)過程中多假設(shè)激勵為白噪聲,造成在應(yīng)用上有一定的局限性。提出了一種基于傳遞率函數(shù)的運(yùn)行模態(tài)分析方法,無需采用白噪聲假設(shè),利用兩種不同載荷情況下的傳遞率函數(shù)構(gòu)造有理函數(shù),通過有理分式Forsythe正交多項式法對其進(jìn)行擬合,得到模態(tài)頻率、阻尼和振型參數(shù)。最后,采用機(jī)翼模型仿真算例和懸臂梁實驗,驗證了在非白噪聲激勵情況下該方法的有效性與可靠性。

    模態(tài)參數(shù)識別; 運(yùn)行模態(tài)分析; 傳遞率函數(shù); 有理函數(shù); 正交多項式

    引 言

    隨著科學(xué)技術(shù)的發(fā)展,在航空航天、車輛、機(jī)械、土木等工程領(lǐng)域,結(jié)構(gòu)動力學(xué)問題日益突出。模態(tài)分析技術(shù)已成為解決此類問題不可或缺的手段,其中,傳統(tǒng)的實驗?zāi)B(tài)分析方法(experimental modal analysis,簡稱EMA)是在實驗室中完成,實驗結(jié)果準(zhǔn)確,相關(guān)的實驗方法也較為完善[1]。運(yùn)行模態(tài)分析(operational modal analysis,簡稱OMA)是只從振動響應(yīng)信號中提取模態(tài)參數(shù)的一項技術(shù)[2]。這與EMA相比,OMA不需要人工激勵,測試成本低且能在系統(tǒng)正常運(yùn)行情況下進(jìn)行測試,更能反映系統(tǒng)真實的動力學(xué)性能。

    近年來發(fā)展了許多運(yùn)行模態(tài)識別方法。其中,時域方法占大多數(shù)并且應(yīng)用更為廣泛,代表性的有隨機(jī)減量(random decrement technique,簡稱RDT)法[3]、基于ARMA模型的時序分析法[4]、自然激勵識別技術(shù)(natural excitation technique,簡稱NExT)[5]和隨機(jī)子空間(stochastic subspace identification,簡稱SSI)法[6]等。與時域方法相比,在頻域方面運(yùn)行模態(tài)識別方法較少。早期主要采用基于響應(yīng)功率譜密度函數(shù)進(jìn)行模態(tài)參數(shù)識別的峰值選取法(peak-picking,簡稱PP)。1999年,Brinker等[7]提出了頻域分解(frequency domain decomposition,簡稱FDD)法,該方法在工程領(lǐng)域應(yīng)用較多,但對于阻尼的識別精度不高。王彤等[8-9]提出頻域空間域分解法,該方法通過增強(qiáng)功率譜得到準(zhǔn)確的頻率和阻尼。這些方法都基于環(huán)境激勵為白噪聲的假設(shè)。在實際的工程應(yīng)用中,系統(tǒng)所受的環(huán)境激勵是豐富多樣的,并非是理想的白噪聲類型。若系統(tǒng)所受的環(huán)境激勵不是白噪聲,采用傳統(tǒng)的運(yùn)行模態(tài)分析技術(shù)進(jìn)行測試分析會帶來誤差,甚至得到錯誤的結(jié)果?;趥鬟f率函數(shù)的運(yùn)行模態(tài)分析方法對于環(huán)境激勵力的性質(zhì)要求較低,激勵力可以是白噪聲、脈沖激勵和諧波激勵等。由于它無需“白噪聲激勵”這一前提,因此一方面擴(kuò)大了運(yùn)行模態(tài)分析技術(shù)的使用范圍,另一方面也提高了運(yùn)行模態(tài)參數(shù)識別的可靠度與精度。因此,基于傳遞率的運(yùn)行模態(tài)分析方法[10]的研究具有較大的工程意義。

    筆者采用頻響函數(shù)的Hv計算來估計傳遞率函數(shù)。利用有理分式正交多項式法[11]擬合構(gòu)造的有理函數(shù),得到系統(tǒng)的模態(tài)頻率與阻尼,并通過對傳遞率函數(shù)矩陣奇異值分解得到結(jié)構(gòu)的振型。通過機(jī)翼模型仿真算例和懸臂梁實驗檢驗了該算法在非白噪聲激勵下的有效性和可靠性。

    1 理論背景

    與頻響函數(shù)的定義類似(運(yùn)動響應(yīng)/激勵力),傳遞率函數(shù)的定義為(運(yùn)動響應(yīng)/運(yùn)動輸入)。在拉普拉斯域中,輸出自由度i與參考自由度j的傳遞率函數(shù)Tij定義為

    Tij(s)=Xi(s)/Xj(s)

    (1)

    φir/φjr

    (2)

    其中:φir和φjr為系統(tǒng)第i和j自由度的第r階振型。

    φir/φjr

    (3)

    由式(2)和式(3)可得

    (4)

    (5)

    其倒數(shù)為

    (6)

    (7)

    其中:p,q分別為分子分母多項式的階數(shù);ak和bk為多項式分子分母多項式系數(shù)向量。

    為減少方程求解過程中的病態(tài)性并解耦系統(tǒng)矩陣,采用Forsythe正交多項式有理分式形式

    (8)

    通過最小二乘擬合獲得系統(tǒng)的模態(tài)頻率與阻尼。

    系統(tǒng)振型可通過對構(gòu)造的傳遞率函數(shù)矩陣在系統(tǒng)模態(tài)頻率處奇異值分解得到[12]。

    (9)

    對式(9)奇異值分解取奇異值向量第一列,即得到結(jié)構(gòu)的振型。

    2 仿真算例

    彎扭二自由度機(jī)翼模型被廣泛應(yīng)用于飛機(jī)飛行顫振分析的仿真。該模型仿真真實機(jī)翼的截面,模型平動和轉(zhuǎn)動可仿真真實飛機(jī)機(jī)翼的彎曲和扭轉(zhuǎn)模態(tài)。機(jī)翼物理模型和簡化的力學(xué)模型如圖1,2所示。

    圖1 機(jī)翼物理模型Fig.1 The physical model of the wing

    圖2 機(jī)翼力學(xué)模型Fig.2 The mechanical model of the wing

    機(jī)翼模型物理參數(shù)設(shè)置如表1所示,對模型添加比例阻尼C=0.000 4K。

    表1 機(jī)翼的物理參數(shù)

    為驗證該算法對非白噪聲激勵情況下模態(tài)參數(shù)的識別效果,激勵信號分別采用兩種典型的非白噪聲信號仿真得到響應(yīng)信號:情況1,粉紅噪聲(Pink noise);情況2,間歇脈沖激勵。通過改變激勵點(diǎn),得到兩種載荷情況下的加速度響應(yīng)數(shù)據(jù)。

    仿真采樣頻率為64 Hz,數(shù)據(jù)長度為6 400點(diǎn)。仿真得到兩種信號激勵情況下的響應(yīng)信號分別如圖3和圖4所示。對應(yīng)的激勵信號自功率譜(power spectrum density,簡稱PSD)如圖5和圖6所示,為非白噪聲信號。

    選取2號點(diǎn)為參考點(diǎn),根據(jù)兩種載荷下的加速度響應(yīng)數(shù)據(jù)構(gòu)造有理函數(shù),情況1和情況2下有理函數(shù)及其擬合如圖7和圖8所示。

    圖3 模型的仿真時間響應(yīng)(情況1)

    圖4 模型的仿真時間響應(yīng)(情況2)

    圖5 激勵信號功率譜(情況1)

    圖6 激勵信號功率譜(情況2)

    圖7 有理函數(shù)及其擬合(情況1)

    圖8 有理函數(shù)及其擬合(情況2)

    圖9和圖10為采用響應(yīng)譜和傳遞率函數(shù)矩陣奇異值分解得到的模態(tài)指示曲線。該仿真算例中數(shù)據(jù)長度為6 400點(diǎn),由于激勵信號非白噪聲,且采集得到的響應(yīng)信號數(shù)據(jù)量較少,采用傳統(tǒng)方法得到的模態(tài)指示曲線毛刺較多,不夠清晰。

    識別得到情況1和情況2下機(jī)翼的模態(tài)參數(shù)與其理論值對比如表2和表3所示。

    圖9 模態(tài)指示曲線(情況1)Fig.9 The curve of mode indicator (case 1)

    圖10 模態(tài)指示曲線(情況2)Fig.10 The curve of mode indicator (case 2)

    理論頻率/Hz識別頻率/Hz頻率誤差/%理論阻尼/%識別阻尼/%阻尼誤差/%8.978.970.001.131.12-0.8814.5214.51-0.071.821.830.55

    表3 結(jié)構(gòu)模態(tài)參數(shù)及誤差分析(情況2)

    由上表可看出,該算法對非白噪聲激勵情況下也適用,且具有較高的識別精度。另外該算法與傳統(tǒng)的方法相比所需數(shù)據(jù)量少,提高了識別效率。經(jīng)測試,在Inter(R) Core(TM) i3-2310M處理器下,參數(shù)識別耗時0.69 s,識別速度較快。

    3 實 驗

    采用基于傳遞率函數(shù)的運(yùn)行模態(tài)分析方法對懸臂梁進(jìn)行模態(tài)參數(shù)識別。實驗裝置如圖11所示,在一根鐵質(zhì)的懸臂梁上等間距布置4個ICP加速度傳感器,激勵采用螺絲釘人工連續(xù)敲擊(激勵為脈沖激勵),并且采用分別敲擊1號和3號點(diǎn)得到兩組不同載荷情況下的響應(yīng)信號。

    采樣頻率為1 280 Hz,采樣時間設(shè)置為12.8 s,采樣點(diǎn)數(shù)為16 384點(diǎn)。首先敲擊1號點(diǎn),采集得到第一種載荷情況下的輸出信號,該激勵情況下2號點(diǎn)的輸出信號,如圖12所示。然后再敲擊3號點(diǎn),采集得到第二種載荷情況下的輸出信號。

    因為梁的自由端4號點(diǎn)振幅較大,相對信噪比較高,故選取該點(diǎn)作為參考點(diǎn)。選取3號點(diǎn)分析,構(gòu)造得到基于傳遞率函數(shù)的有理函數(shù)。采用有理分式正交多項式法對構(gòu)造得到的有理多項式函數(shù)進(jìn)行擬合,擬合結(jié)果如圖13所示。

    為驗證該算法識別結(jié)果的精度,作為對比,采用力錘激勵對懸臂梁做EMA實驗。由于采用力錘激勵,人工激勵輸入可測,且測試信噪比高,識別結(jié)果精確可靠,故可視EMA識別結(jié)果為結(jié)構(gòu)真實參數(shù)?;趥鬟f率函數(shù)(transmissibility functions,簡稱TF)法和頻域空間域分解(frequency spatial domain decomposition,簡稱FSDD)法識別結(jié)果與EMA頻率和阻尼參數(shù)識別結(jié)果對比如表4和表5所示。

    圖11 實驗裝置圖

    圖12 2號點(diǎn)的響應(yīng)

    圖13 有理函數(shù)擬合

    表4 頻率參數(shù)及誤差分析

    表5 阻尼參數(shù)及誤差分析

    由表4,5可看出,基于傳遞率函數(shù)法識別得到的頻率和阻尼參數(shù)與EMA識別得到的結(jié)果相對誤差較小,頻率誤差1%以內(nèi),阻尼誤差10%以內(nèi),識別結(jié)果可靠。該算法識別得參數(shù)共耗時0.739 s,具有較高的識別效率。

    在非白噪聲和短數(shù)據(jù)情況下,采用基于傳遞率函數(shù)法仍能得到較為理想的有理函數(shù),能清晰的指示出模態(tài)頻率。另外,兩種方法得到的模態(tài)指示曲線如圖14所示,可以看到功率譜法得到的模態(tài)指示曲線毛刺較多,模態(tài)指示不清晰。

    圖14 模態(tài)指示曲線Fig.14 The curve of mode indicator

    因此基于傳遞率函數(shù)的運(yùn)行模態(tài)識別方法在非白噪聲激勵和短數(shù)據(jù)情況下具有一定的優(yōu)越性。該算法識別得到懸臂梁的振型如圖15所示。

    圖15 懸臂梁的4階模態(tài)振型Fig.15 Four mode shapes of simply cantilever beam

    4 結(jié) 論

    1) 該方法需要至少二種不同載荷情況下的響應(yīng)信號來識別模態(tài)參數(shù)。在實際工程應(yīng)用中,由于環(huán)境激勵的隨機(jī)性,兩種不同載荷情況下的數(shù)據(jù)理論上可通過不同時間段采集的數(shù)據(jù)得到。

    2) 仿真和實驗表明,該方法即使在非白噪聲信號激勵和短數(shù)據(jù)情況下也能得到較為理想的識別結(jié)果。因此,該方法有望拓展OMA的應(yīng)用范圍。

    3) 適當(dāng)提高擬合階次能得到更好的識別結(jié)果,與此同時帶來的計算模態(tài)可通過判斷阻尼大小、查看有理函數(shù)擬合圖以及振型動畫等方法予以剔除。

    [1] 王彤, 張令彌.多頻段擬合的正交多項式方法及其Matlab工具箱[J].振動工程學(xué)報,2004,25(6):560-564.

    Wang Tong, Zhang Lingmi. A multiple frequency band fitting method using orthogonal polynomials and corresponding Matlab toolbox[J]. Journal of Vibration Engineering, 2004,25(6):560-564. (in Chinese)

    [2] 張義民,張守元,李鶴,等.運(yùn)行模態(tài)分析中諧波模態(tài)識別方法研究及應(yīng)用[J].振動、測試與診斷,2008,28(3):197-200.

    Zhang Yimin, Zhang Shouyuan, Li He, et al. Harmonic mode identification in the operational modal analysis and its application[J]. Journal of Vibration, Measurement & Diagnosis, 2008,28(3):197-200. (in Chinese)

    [3] Ibrahim S R. Random decrement technique for modal identification of structures[J]. AIAA Journal of Spacecraft and Rockets, 1977,14(11):696-700.

    [4] Bonnecase D, Provosto M. Application of a multidimensional ARMA model to modal analysis under natural excitation[C]∥Proceeding of the 8th International Modal Analysis Conference (IMAC). Kissemee, Florida, USA: Society for Experimental Mechanics, 1990:382-388.

    [5] James G H, Carne T G, Lauffer J P. The natural excitation technique (NExT) for modal parameter extraction from operating structures[J]. Journal of Analytical and Experimental Modal Analysis, 1995,10(4):260-277.

    [6] Overschee P V, Moor B D. Identification for linear systems: theory, implementation, application[M]. UK: Kluwer Academic Publishers, 1996:57-89.

    [7] Brincker R, Zhang L M, Anderson P. Modal identification from ambient response using frequency domain decomposition[C]∥Proceeding of the 18th International Modal Analysis Conference (IMAC). San Antonio, USA: Society for Experimental Mechanics, 2000:628-630.

    [8] Wang Tong, Zhang Lingmi, Tamura Yukio. An operational modal analysis method in frequency and spatial domain[J]. Earthquake Engineering and Engineering Vibration, 2005,4(2):295-300.

    [9] 王彤,張令彌.運(yùn)行模態(tài)分析的頻域空間域分解法及其應(yīng)用[J].航空學(xué)報,2006,27(1):62-66.

    Wang Tong, Zhang Lingmi. Frequency and spatial domain decomposition for operational modal analysis and its application[J]. Acta Aeronautica et Astronautica Sinica, 2006,27(1):62-66. (in Chinese)

    [10]Devriendt C, Guillaume P. The use of transmissibility measurements in output-only modal analysis[J]. Mechanical Systems and Signal Processing, 2007,21(7):2689-2696.

    [11]王彤,張令彌.有理分式正交多項式頻響函數(shù)模態(tài)參數(shù)識別[J].航空學(xué)報,2003,24(2):140-143.

    Wang Tong, Zhang Lingmi. Modal identification with frequency response function based on rational fraction orthogonal polynomials[J]. Acta Aeronautica et Astronautica Sinica, 2003,24(2):140-143. (in Chinese)

    [12]Araujo I G, Laier J E. Operational modal analysis using SVD of power spectral density transmissibility matrices[J]. Mechanical Systems and Signal Processing, 2014,46(1):129-145.

    10.16450/j.cnki.issn.1004-6801.2015.05.023

    *南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國家重點(diǎn)實驗室自主研究課題資助項目(0113G01);江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目

    2014-05-27;

    2014-09-01

    O329; TB123

    張永年,男,1990年9月生,碩士。主要研究方向為振動測試與模態(tài)參數(shù)識別。 E-mail:zyn@dytac.com.cn

    猜你喜歡
    機(jī)翼阻尼模態(tài)
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項的擴(kuò)散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    變時滯間隙非線性機(jī)翼顫振主動控制方法
    具阻尼項的Boussinesq型方程的長時間行為
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    機(jī)翼跨聲速抖振研究進(jìn)展
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    由單個模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
    計算物理(2014年2期)2014-03-11 17:01:39
    基于模糊自適應(yīng)的高超聲速機(jī)翼顫振的主動控制
    国产成人精品久久久久久| 久久精品国产自在天天线| 中文资源天堂在线| 午夜精品国产一区二区电影 | 美女大奶头视频| 国产免费男女视频| 午夜久久久久精精品| 午夜爱爱视频在线播放| 亚洲精品国产av成人精品| 91在线精品国自产拍蜜月| 99九九线精品视频在线观看视频| 七月丁香在线播放| 国产高潮美女av| 网址你懂的国产日韩在线| 国产激情偷乱视频一区二区| 大又大粗又爽又黄少妇毛片口| 国产综合懂色| 亚洲人成网站在线播| 免费搜索国产男女视频| 淫秽高清视频在线观看| 午夜福利成人在线免费观看| 麻豆精品久久久久久蜜桃| 国产69精品久久久久777片| a级毛色黄片| 嘟嘟电影网在线观看| 欧美成人免费av一区二区三区| 成人二区视频| 欧美丝袜亚洲另类| a级毛色黄片| 嘟嘟电影网在线观看| 老司机福利观看| 日本三级黄在线观看| 日韩欧美国产在线观看| 色尼玛亚洲综合影院| 国产高清三级在线| 真实男女啪啪啪动态图| 99久久人妻综合| 中文字幕人妻熟人妻熟丝袜美| 99热这里只有精品一区| 久久精品国产99精品国产亚洲性色| 免费黄色在线免费观看| 久久人人爽人人片av| 桃色一区二区三区在线观看| 久久精品国产99精品国产亚洲性色| 嘟嘟电影网在线观看| 日本三级黄在线观看| videos熟女内射| 国产亚洲av片在线观看秒播厂 | 深爱激情五月婷婷| 成年女人永久免费观看视频| 国产精品日韩av在线免费观看| 激情 狠狠 欧美| 一区二区三区四区激情视频| 国产在线男女| 高清视频免费观看一区二区 | 九草在线视频观看| 久久久久久久久久久免费av| 三级男女做爰猛烈吃奶摸视频| 国产精品1区2区在线观看.| 高清午夜精品一区二区三区| 1000部很黄的大片| 午夜福利高清视频| 亚洲精品,欧美精品| 亚洲精品成人久久久久久| 亚洲精品,欧美精品| 只有这里有精品99| 五月伊人婷婷丁香| 97人妻精品一区二区三区麻豆| 日韩欧美精品v在线| 久久99热这里只有精品18| 国产男人的电影天堂91| 狂野欧美激情性xxxx在线观看| 日韩成人av中文字幕在线观看| 精品久久久久久久末码| 禁无遮挡网站| 韩国av在线不卡| 韩国高清视频一区二区三区| 色尼玛亚洲综合影院| av在线老鸭窝| 国产精品野战在线观看| 熟妇人妻久久中文字幕3abv| 亚洲成人av在线免费| 天天躁夜夜躁狠狠久久av| 3wmmmm亚洲av在线观看| 久久久久久伊人网av| 亚洲熟妇中文字幕五十中出| 九九热线精品视视频播放| 美女xxoo啪啪120秒动态图| 嫩草影院入口| av女优亚洲男人天堂| 青春草国产在线视频| 亚洲人成网站在线播| 久久久精品大字幕| 人人妻人人看人人澡| 尾随美女入室| 一级黄色大片毛片| 成人毛片60女人毛片免费| 亚洲精品久久久久久婷婷小说 | 国产一区二区亚洲精品在线观看| 久久久久久大精品| 美女xxoo啪啪120秒动态图| 欧美xxxx黑人xx丫x性爽| 1024手机看黄色片| 亚洲欧美精品综合久久99| 老女人水多毛片| 大又大粗又爽又黄少妇毛片口| 免费观看精品视频网站| 美女被艹到高潮喷水动态| 天天一区二区日本电影三级| 亚洲中文字幕日韩| 国内精品美女久久久久久| 免费看a级黄色片| 99热网站在线观看| 激情 狠狠 欧美| av国产久精品久网站免费入址| 精品欧美国产一区二区三| 女人被狂操c到高潮| 国产高清三级在线| 看黄色毛片网站| 国产精品1区2区在线观看.| 欧美三级亚洲精品| 国产精品综合久久久久久久免费| 成人三级黄色视频| АⅤ资源中文在线天堂| 成人午夜高清在线视频| 精品一区二区三区视频在线| 啦啦啦啦在线视频资源| 十八禁国产超污无遮挡网站| 亚洲欧美成人精品一区二区| 免费av不卡在线播放| 日日啪夜夜撸| 最新中文字幕久久久久| 中文资源天堂在线| 免费av不卡在线播放| 日本五十路高清| 七月丁香在线播放| 青青草视频在线视频观看| 免费观看的影片在线观看| 国产69精品久久久久777片| 亚洲人成网站高清观看| 国产亚洲精品av在线| av天堂中文字幕网| 男人舔奶头视频| 国产一区二区在线观看日韩| 男女边吃奶边做爰视频| 丝袜美腿在线中文| 国产精品熟女久久久久浪| 久久精品国产99精品国产亚洲性色| 最近最新中文字幕大全电影3| 色5月婷婷丁香| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久久丰满| 国产日韩欧美在线精品| 欧美日韩国产亚洲二区| 熟女电影av网| 深夜a级毛片| 国产成人福利小说| 一边摸一边抽搐一进一小说| 成年免费大片在线观看| 国产亚洲av片在线观看秒播厂 | 国产一级毛片在线| 可以在线观看毛片的网站| 色综合色国产| 久久久亚洲精品成人影院| av免费在线看不卡| 草草在线视频免费看| 长腿黑丝高跟| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美色视频一区免费| 国产一级毛片在线| 久久精品熟女亚洲av麻豆精品 | 男女那种视频在线观看| 岛国毛片在线播放| 久久久午夜欧美精品| 国产精品,欧美在线| 一夜夜www| 中文字幕精品亚洲无线码一区| 国产 一区 欧美 日韩| www.av在线官网国产| 精品人妻视频免费看| 免费电影在线观看免费观看| 两个人视频免费观看高清| 欧美xxxx性猛交bbbb| 亚洲在线自拍视频| 亚洲色图av天堂| 黄片无遮挡物在线观看| 久久久久久国产a免费观看| 日本一二三区视频观看| 91精品国产九色| 国产又色又爽无遮挡免| 熟女电影av网| 亚洲欧美成人综合另类久久久 | 国产精品人妻久久久久久| 中文亚洲av片在线观看爽| 在线观看一区二区三区| 我要搜黄色片| 三级毛片av免费| 99热这里只有是精品50| 波野结衣二区三区在线| 亚洲精品日韩av片在线观看| 嫩草影院精品99| 国产人妻一区二区三区在| 久久久a久久爽久久v久久| 白带黄色成豆腐渣| 日韩欧美精品v在线| 狂野欧美激情性xxxx在线观看| 亚洲四区av| 亚洲天堂国产精品一区在线| 99久久九九国产精品国产免费| 黄片无遮挡物在线观看| 91久久精品国产一区二区三区| 午夜精品在线福利| 别揉我奶头 嗯啊视频| 少妇熟女aⅴ在线视频| 欧美又色又爽又黄视频| 99在线人妻在线中文字幕| 久久国内精品自在自线图片| 久久久久久久久久久丰满| 欧美不卡视频在线免费观看| 深夜a级毛片| 久久综合国产亚洲精品| 精品无人区乱码1区二区| 成人欧美大片| 麻豆久久精品国产亚洲av| 日韩大片免费观看网站 | 日本三级黄在线观看| www.av在线官网国产| 日日啪夜夜撸| or卡值多少钱| 国产麻豆成人av免费视频| 亚洲美女搞黄在线观看| 特级一级黄色大片| 只有这里有精品99| 中文字幕熟女人妻在线| 最近的中文字幕免费完整| 久久这里有精品视频免费| 免费观看在线日韩| 亚洲精品乱久久久久久| 日韩欧美精品免费久久| 国产精品永久免费网站| 国产成人freesex在线| 成人毛片a级毛片在线播放| 成人国产麻豆网| 亚洲国产最新在线播放| 亚洲久久久久久中文字幕| 日韩国内少妇激情av| 欧美日韩在线观看h| 久久综合国产亚洲精品| 色视频www国产| 亚洲性久久影院| 黄片无遮挡物在线观看| 免费看美女性在线毛片视频| 成年免费大片在线观看| 中国美白少妇内射xxxbb| 人人妻人人澡欧美一区二区| 国产亚洲最大av| 麻豆成人午夜福利视频| 国产亚洲av片在线观看秒播厂 | 99热精品在线国产| 一区二区三区高清视频在线| 日本猛色少妇xxxxx猛交久久| 国产片特级美女逼逼视频| 成人美女网站在线观看视频| 日本黄色视频三级网站网址| 免费不卡的大黄色大毛片视频在线观看 | 在线观看一区二区三区| 99久久无色码亚洲精品果冻| 国产久久久一区二区三区| 亚洲av二区三区四区| 婷婷六月久久综合丁香| av国产久精品久网站免费入址| 日韩欧美国产在线观看| 你懂的网址亚洲精品在线观看 | 国产在视频线在精品| 最近中文字幕2019免费版| 国产私拍福利视频在线观看| 国产精华一区二区三区| 国产色婷婷99| 三级国产精品欧美在线观看| 亚洲av成人av| 国产麻豆成人av免费视频| 天堂av国产一区二区熟女人妻| 又粗又爽又猛毛片免费看| 国产精品三级大全| 亚洲欧美一区二区三区国产| 国产又色又爽无遮挡免| 日本与韩国留学比较| 两个人视频免费观看高清| 亚洲美女搞黄在线观看| 精品国产一区二区三区久久久樱花 | 小蜜桃在线观看免费完整版高清| 亚洲精品乱码久久久久久按摩| 淫秽高清视频在线观看| 欧美3d第一页| 国产精品乱码一区二三区的特点| 看非洲黑人一级黄片| 卡戴珊不雅视频在线播放| av线在线观看网站| 你懂的网址亚洲精品在线观看 | 国产成人a区在线观看| 亚洲精品乱码久久久久久按摩| 两个人视频免费观看高清| 亚洲国产精品久久男人天堂| www.av在线官网国产| 菩萨蛮人人尽说江南好唐韦庄 | 高清毛片免费看| 亚洲成色77777| 美女黄网站色视频| 久久久国产成人免费| 99热这里只有是精品50| 成人亚洲欧美一区二区av| av免费在线看不卡| 亚洲国产精品成人久久小说| 97超碰精品成人国产| 深夜a级毛片| 亚洲高清免费不卡视频| 国产成人a区在线观看| 国内少妇人妻偷人精品xxx网站| 淫秽高清视频在线观看| 男人狂女人下面高潮的视频| 精品酒店卫生间| 美女黄网站色视频| 国产人妻一区二区三区在| 欧美日韩国产亚洲二区| 日韩精品有码人妻一区| 日韩欧美精品免费久久| 欧美潮喷喷水| 免费看a级黄色片| 毛片一级片免费看久久久久| 国产精品麻豆人妻色哟哟久久 | 婷婷色麻豆天堂久久 | av在线亚洲专区| 亚洲精品成人久久久久久| 一个人观看的视频www高清免费观看| 中文字幕亚洲精品专区| 国产探花极品一区二区| 哪个播放器可以免费观看大片| 我要搜黄色片| 超碰97精品在线观看| 国产亚洲精品av在线| 男人舔女人下体高潮全视频| 天天躁日日操中文字幕| 国产免费又黄又爽又色| 日本免费一区二区三区高清不卡| 亚洲av一区综合| 亚洲不卡免费看| 一区二区三区高清视频在线| 视频中文字幕在线观看| 别揉我奶头 嗯啊视频| 欧美精品一区二区大全| 听说在线观看完整版免费高清| 国产一级毛片在线| 全区人妻精品视频| 国产精品日韩av在线免费观看| 久久精品国产99精品国产亚洲性色| 一级黄色大片毛片| 色综合站精品国产| 韩国高清视频一区二区三区| 熟妇人妻久久中文字幕3abv| 国产精品,欧美在线| 三级毛片av免费| 大香蕉久久网| 亚洲18禁久久av| 麻豆乱淫一区二区| 国产高清不卡午夜福利| 男女边吃奶边做爰视频| 国产精品久久久久久久久免| 在线天堂最新版资源| 久久精品综合一区二区三区| 晚上一个人看的免费电影| 全区人妻精品视频| 亚洲av电影在线观看一区二区三区 | ponron亚洲| 一区二区三区四区激情视频| 欧美一级a爱片免费观看看| 久久欧美精品欧美久久欧美| 亚洲婷婷狠狠爱综合网| 国产激情偷乱视频一区二区| 欧美bdsm另类| 99在线视频只有这里精品首页| 国产伦一二天堂av在线观看| 丰满人妻一区二区三区视频av| 精品无人区乱码1区二区| 人人妻人人澡人人爽人人夜夜 | ponron亚洲| 一个人免费在线观看电影| 中文在线观看免费www的网站| ponron亚洲| 久久综合国产亚洲精品| av线在线观看网站| 2022亚洲国产成人精品| 国产精品久久视频播放| 成人综合一区亚洲| 欧美bdsm另类| 亚洲四区av| 成人高潮视频无遮挡免费网站| 丰满人妻一区二区三区视频av| 麻豆乱淫一区二区| 亚洲中文字幕日韩| 午夜亚洲福利在线播放| 色尼玛亚洲综合影院| 亚洲精品久久久久久婷婷小说 | 国产乱人偷精品视频| 2022亚洲国产成人精品| 午夜免费激情av| 特级一级黄色大片| 亚洲精品乱久久久久久| 国产精品国产三级国产专区5o | 久久久久九九精品影院| 熟妇人妻久久中文字幕3abv| 亚州av有码| 欧美丝袜亚洲另类| av.在线天堂| 91在线精品国自产拍蜜月| 久久这里有精品视频免费| 日韩制服骚丝袜av| 国产亚洲一区二区精品| 日本黄大片高清| 一个人免费在线观看电影| 精品国内亚洲2022精品成人| 精品人妻视频免费看| 午夜福利视频1000在线观看| 精品少妇黑人巨大在线播放 | 国产又黄又爽又无遮挡在线| 精品人妻偷拍中文字幕| 国产精品久久视频播放| 亚洲性久久影院| 亚洲国产欧洲综合997久久,| 精品99又大又爽又粗少妇毛片| 日韩人妻高清精品专区| 免费观看性生交大片5| 十八禁国产超污无遮挡网站| 久久精品人妻少妇| 天天躁日日操中文字幕| 国产精品人妻久久久久久| 国产成年人精品一区二区| 精品久久久久久成人av| av在线观看视频网站免费| 99久久精品国产国产毛片| 91狼人影院| 日韩欧美三级三区| 干丝袜人妻中文字幕| 国产精品综合久久久久久久免费| 亚洲婷婷狠狠爱综合网| 麻豆成人av视频| 床上黄色一级片| 国产又黄又爽又无遮挡在线| 亚洲乱码一区二区免费版| 国产精品一区www在线观看| 波野结衣二区三区在线| 少妇裸体淫交视频免费看高清| 国产色婷婷99| 国产精品av视频在线免费观看| 久久精品久久久久久久性| 亚洲电影在线观看av| av视频在线观看入口| kizo精华| 国产成人福利小说| 久久久a久久爽久久v久久| 国产精品国产三级国产av玫瑰| 中文字幕精品亚洲无线码一区| 欧美日韩一区二区视频在线观看视频在线 | 一区二区三区乱码不卡18| 国产伦一二天堂av在线观看| 又粗又硬又长又爽又黄的视频| 男人狂女人下面高潮的视频| 久久精品人妻少妇| 好男人在线观看高清免费视频| 欧美不卡视频在线免费观看| 嫩草影院精品99| 啦啦啦韩国在线观看视频| www.色视频.com| 欧美3d第一页| 纵有疾风起免费观看全集完整版 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人妻制服诱惑在线中文字幕| av又黄又爽大尺度在线免费看 | 女人久久www免费人成看片 | 秋霞伦理黄片| 日韩大片免费观看网站 | 岛国毛片在线播放| 国产一区二区亚洲精品在线观看| 国产日韩欧美在线精品| 免费观看a级毛片全部| 久久精品夜色国产| 免费观看在线日韩| 成人无遮挡网站| 欧美人与善性xxx| 精品久久久久久久久av| 亚洲天堂国产精品一区在线| 欧美bdsm另类| 久热久热在线精品观看| 婷婷色麻豆天堂久久 | 欧美三级亚洲精品| 午夜老司机福利剧场| 一个人看视频在线观看www免费| 国产私拍福利视频在线观看| 中文字幕精品亚洲无线码一区| 成人毛片a级毛片在线播放| 免费观看的影片在线观看| 精品不卡国产一区二区三区| 69av精品久久久久久| 高清视频免费观看一区二区 | av天堂中文字幕网| 一个人观看的视频www高清免费观看| 日韩成人伦理影院| 嫩草影院精品99| 日韩一区二区视频免费看| АⅤ资源中文在线天堂| 最近最新中文字幕大全电影3| 精品久久久久久电影网 | 亚洲欧洲国产日韩| 免费看日本二区| 少妇人妻一区二区三区视频| 精品人妻熟女av久视频| h日本视频在线播放| 卡戴珊不雅视频在线播放| 久久精品国产亚洲av涩爱| 亚洲最大成人av| 大又大粗又爽又黄少妇毛片口| 国产精品一区二区三区四区久久| 精品久久久久久久久久久久久| 亚洲精品成人久久久久久| 人妻系列 视频| 亚洲精品,欧美精品| 亚洲成人久久爱视频| 亚洲电影在线观看av| 性色avwww在线观看| 久久久久国产网址| 一区二区三区四区激情视频| 精品久久久久久久人妻蜜臀av| 亚洲精品乱久久久久久| 婷婷色av中文字幕| 2021天堂中文幕一二区在线观| 国产精品一区www在线观看| 久久久久久九九精品二区国产| 男人和女人高潮做爰伦理| 我的老师免费观看完整版| 一区二区三区免费毛片| 欧美区成人在线视频| 日韩av在线大香蕉| 三级经典国产精品| 久久99热这里只频精品6学生 | 久久99蜜桃精品久久| 日本黄色片子视频| 欧美日本亚洲视频在线播放| 精品无人区乱码1区二区| 日日啪夜夜撸| 免费在线观看成人毛片| 免费看a级黄色片| 久久草成人影院| 久久精品久久久久久噜噜老黄 | 免费播放大片免费观看视频在线观看 | 成年免费大片在线观看| 亚洲欧美日韩无卡精品| 91精品国产九色| 久久精品夜夜夜夜夜久久蜜豆| 欧美另类亚洲清纯唯美| 九九久久精品国产亚洲av麻豆| 国产91av在线免费观看| 亚洲av成人精品一二三区| 高清毛片免费看| 最近最新中文字幕免费大全7| 中文字幕制服av| 久久婷婷人人爽人人干人人爱| 蜜桃久久精品国产亚洲av| 国产精品久久久久久精品电影| 国产av一区在线观看免费| 国产精品久久久久久精品电影| 真实男女啪啪啪动态图| 国产av不卡久久| 国产又色又爽无遮挡免| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 午夜福利在线观看吧| 亚洲第一区二区三区不卡| 简卡轻食公司| 神马国产精品三级电影在线观看| 亚洲一区高清亚洲精品| 久久精品久久久久久噜噜老黄 | 精华霜和精华液先用哪个| 99热这里只有精品一区| 国产成人免费观看mmmm| 亚洲三级黄色毛片| 亚洲最大成人av| 国产精品99久久久久久久久| 亚洲aⅴ乱码一区二区在线播放| 国产高清不卡午夜福利| 久久久色成人| 午夜福利在线在线| 亚洲国产成人一精品久久久| 亚洲综合色惰| 啦啦啦韩国在线观看视频| 国产视频首页在线观看| 亚洲在久久综合| videossex国产| 亚洲中文字幕一区二区三区有码在线看| 中文乱码字字幕精品一区二区三区 | 国产午夜精品论理片| 少妇猛男粗大的猛烈进出视频 | 麻豆精品久久久久久蜜桃| 黄片无遮挡物在线观看| 高清毛片免费看| 亚洲av福利一区| 小蜜桃在线观看免费完整版高清| 日本免费a在线| 国内精品一区二区在线观看| 性色avwww在线观看| 22中文网久久字幕| 久久精品久久久久久久性| 亚洲欧洲国产日韩| 丝袜喷水一区| 综合色av麻豆| 99久久九九国产精品国产免费| 美女cb高潮喷水在线观看|