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

    傳遞路徑分析中計(jì)算子系統(tǒng)頻響函數(shù)的方法

    2018-09-29 11:14廖旭暉李舜酩孟浩東
    振動(dòng)工程學(xué)報(bào) 2018年4期
    關(guān)鍵詞:子系統(tǒng)

    廖旭暉 李舜酩 孟浩東

    摘要: 傳遞路徑分析是在振動(dòng)噪聲控制領(lǐng)域被廣泛應(yīng)用的一種有效方法。傳遞路徑分析中將振動(dòng)系統(tǒng)分成主動(dòng)部分、被動(dòng)部分以及連接主、被動(dòng)部分的若干傳遞路徑。在傳遞路徑分析中需要對(duì)被動(dòng)部分的頻響函數(shù)進(jìn)行測(cè)量。傳統(tǒng)的傳遞路徑分析需要先拆除子結(jié)構(gòu)然后再測(cè)量頻響函數(shù),測(cè)試過程十分繁瑣。提出了一種全新的方法來計(jì)算子系統(tǒng)的頻響函數(shù),直接由整個(gè)系統(tǒng)的頻響函數(shù)矩陣推導(dǎo)得到子系統(tǒng)的頻響函數(shù)矩陣的計(jì)算公式。該方法不需要對(duì)子系統(tǒng)進(jìn)行物理解耦,大大縮短了測(cè)量子系統(tǒng)頻響函數(shù)所需要的時(shí)間。數(shù)值算例和實(shí)驗(yàn)均驗(yàn)證了該方法的正確性和有效性。

    關(guān)鍵詞: 子系統(tǒng); 傳遞路徑分析; 頻響函數(shù); 解耦

    中圖分類號(hào): O321; TB123文獻(xiàn)標(biāo)志碼:A文章編號(hào)1004-4523(2018)04-0681-07

    DOI:10.16385/j.cnki.issn.1004-4523.2018.04.016

    引言

    傳遞路徑分析(Transfer Path Analysis, TPA)在車輛的NVH(Noise, Vibration & Harshness)研究、船舶的減振降噪、隔振系統(tǒng)設(shè)計(jì)等振動(dòng)噪聲控制領(lǐng)域應(yīng)用十分廣泛,是一種對(duì)系統(tǒng)進(jìn)行減振降噪十分有效的工程手段。TPA的一個(gè)關(guān)鍵步驟是測(cè)量子系統(tǒng)的頻響函數(shù)。在傳統(tǒng)的TPA中,需要先將子結(jié)構(gòu)之間的連接拆除然后再對(duì)子結(jié)構(gòu)的頻響函數(shù)進(jìn)行測(cè)量。在實(shí)際應(yīng)用中將子系統(tǒng)物理解耦需要耗費(fèi)較長(zhǎng)的時(shí)間,并且物理解耦后子系統(tǒng)的邊界條件通常會(huì)發(fā)生改變,因此使得傳遞路徑分析方法的應(yīng)用受到了較大的限制。為了提高傳遞路徑分析方法的應(yīng)用效率,研究者們提出了許多改進(jìn)的方法,如工況傳遞路徑分析[1-2]、基于直接傳遞率的方法[3-4]、基于組件的傳遞路徑分析[5-6]等。這些方法都規(guī)避了對(duì)子系統(tǒng)頻響函數(shù)的直接測(cè)量,雖然提高了測(cè)試效率,但是由于并不能將響應(yīng)嚴(yán)格地表示成頻響函數(shù)和工況力相乘的形式,因此從本質(zhì)上來說,這些方法都已經(jīng)和傳統(tǒng)的傳遞路徑分析相背離了[7]。傳統(tǒng)的傳遞路徑分析由于有著明確的物理意義和準(zhǔn)確的數(shù)學(xué)表達(dá),因此仍然被廣泛應(yīng)用于工程實(shí)踐或者用來作為其他改進(jìn)方法的參照標(biāo)準(zhǔn)。

    相對(duì)而言,由于不需要解耦,所以獲取整個(gè)系統(tǒng)的頻響函數(shù)更加方便。從哲學(xué)上來說,整體大于部分之和,因此,子系統(tǒng)的動(dòng)力學(xué)特性一定是被包含在整個(gè)系統(tǒng)的動(dòng)力學(xué)特性之中的。如何從作為整體信息的系統(tǒng)頻響函數(shù)矩陣中“提取出”反應(yīng)局部信息的子系統(tǒng)頻響函數(shù)矩陣是非常有意義的一項(xiàng)課題。

    之前的研究者們已經(jīng)做過很多這方面的嘗試,并且取得了許多的研究成果,其中研究較多的是從動(dòng)態(tài)子結(jié)構(gòu)[8-9]的角度來獲取子系統(tǒng)的頻響函數(shù)。Okubo和Miyazaki[10-11]首先提出了根據(jù)整個(gè)系統(tǒng)和已知子系統(tǒng)動(dòng)力學(xué)特性來獲取未知子系統(tǒng)動(dòng)力學(xué)特性的方法。Gontier和Bensaibi[12]提出了一種時(shí)域的方法來對(duì)子系統(tǒng)進(jìn)行解耦。Silva等[13-14]研究了子系統(tǒng)之間具有不同連接特性時(shí)的解耦方法以及識(shí)別連接特性的方法。Kalling等[15]采用狀態(tài)空間模型識(shí)別方法來研究子系統(tǒng)解耦的問題。近年來,一種新的基于逆子結(jié)構(gòu)方法的解耦技術(shù)得到了重視,即將子系統(tǒng)看成是整個(gè)系統(tǒng)和一個(gè)虛擬系統(tǒng)相耦合而成[16],從而實(shí)現(xiàn)解耦。上述這些解耦方法都假設(shè)子系統(tǒng)之間的連接是剛性的。另外,這些子系統(tǒng)解耦方法大多需要已知整個(gè)系統(tǒng)以及剩余部分子系統(tǒng)的動(dòng)力學(xué)特性,才能實(shí)現(xiàn)對(duì)未知子系統(tǒng)的解耦。

    Keersmaekers[17]則另辟蹊徑,提出了所謂的連接保留解耦(Link-Preserving Decoupling, LPD)方法,首先列出整個(gè)系統(tǒng)的動(dòng)剛度矩陣,根據(jù)耦合連接處耦合剛度對(duì)于兩個(gè)子系統(tǒng)完全相同的特性,通過復(fù)雜的數(shù)學(xué)推導(dǎo),得出了子系統(tǒng)頻響函數(shù)矩陣的公式。該方法不需要對(duì)耦合系統(tǒng)進(jìn)行物理解耦,也無需知道相耦合的另一子系統(tǒng)的動(dòng)力學(xué)特性,根據(jù)系統(tǒng)的頻響函數(shù)矩陣即可得到解耦子系統(tǒng)的頻響函數(shù)矩陣。

    值得注意的是,LPD方法與本文所得出的結(jié)論有相似之處,但是該方法并未將內(nèi)部自由度(非耦合處的自由度)考慮進(jìn)來,而傳遞路徑分析中,所感興趣的響應(yīng)點(diǎn)通常是子系統(tǒng)的內(nèi)部自由度,因此,該方法有一定的局限性。

    本文首先給出傳統(tǒng)的TPA模型,在線性假設(shè)和彈性假設(shè)的基礎(chǔ)上,可以將路徑簡(jiǎn)化成線性彈簧和阻尼的組合。然后通過系統(tǒng)的頻響函數(shù)推導(dǎo)出解耦子系統(tǒng)的頻響函數(shù)。最后,通過數(shù)值算例和實(shí)驗(yàn)驗(yàn)證本文所提出方法的正確性。

    1TPA基本理論及耦合系統(tǒng)模型

    TPA理論將振動(dòng)系統(tǒng)分成主動(dòng)部分A、傳遞路徑和被動(dòng)部分B。激勵(lì)力作用在主動(dòng)部分上,通過若干條傳遞路徑將振動(dòng)傳遞到被動(dòng)部分上。典型的TPA模型如圖1所示。

    被動(dòng)部分上的目標(biāo)點(diǎn)的振動(dòng)響應(yīng)可以看成是所有路徑對(duì)該點(diǎn)響應(yīng)的貢獻(xiàn)之和。TPA的基本公式可以表示成yi=∑jHijFj(1)式中yi表示目標(biāo)點(diǎn)i的響應(yīng),F(xiàn)j表示第j條傳遞路徑上傳遞的力,Hij表示力Fj到響應(yīng)yi的頻響函數(shù)。從式(1)中可以發(fā)現(xiàn)運(yùn)用傳遞路徑分析方法必須測(cè)量被動(dòng)部分的頻響函數(shù)。為了簡(jiǎn)化分析,對(duì)圖1所示系統(tǒng)作如下假設(shè):(1)系統(tǒng)為線性的;(2)連接子系統(tǒng)A和B的路徑具有足夠大的彈性,確保路徑傳遞力時(shí)有充分的變形量;(3)子系統(tǒng)A和B之間相互耦合的自由度是一一對(duì)應(yīng)的,即子系統(tǒng)之間不存在多個(gè)自由度對(duì)一個(gè)自由度的耦合。

    假設(shè)(1)保證了可以通過頻響函數(shù)矩陣來描述系統(tǒng),本文的結(jié)論是在線性假設(shè)的前提下得到的。在許多情況下,系統(tǒng)可以做近似線性假設(shè)。假設(shè)(2)保證了對(duì)頻響函數(shù)矩陣進(jìn)行求逆運(yùn)算時(shí),不會(huì)出現(xiàn)奇異性。假設(shè)(3)使得模型更加簡(jiǎn)化,便于分析,且該假設(shè)在大多數(shù)情況下是能夠滿足的。通過上述幾點(diǎn)假設(shè),可以將各傳遞路徑簡(jiǎn)化成線性彈簧和阻尼的組合。上述的傳遞路徑模型實(shí)際上就是兩個(gè)子系統(tǒng)通過線性彈簧和阻尼相互耦合起來的振動(dòng)系統(tǒng)。

    2子系統(tǒng)頻響函數(shù)矩陣計(jì)算方法

    假設(shè)在向量XAC和XBC中各自自由度排列順序是一致的,即XAC的第n個(gè)元素和XBC中第n個(gè)元素正好代表一對(duì)耦合自由度,XAC-XBC則表示由各條路徑上的變形量所構(gòu)成的一組向量。式(8)的物理意義非常清楚,就是將子系統(tǒng)B上的的響應(yīng)表達(dá)成三部分之和。第一部分是作用于子系統(tǒng)B上的外力所產(chǎn)生的對(duì)響應(yīng)的貢獻(xiàn),第二部分是由于連接子系統(tǒng)A和B的路徑的變形所導(dǎo)致的對(duì)響應(yīng)的貢獻(xiàn),第三部分則是由于子系統(tǒng)A內(nèi)部自由度上的位移所導(dǎo)致的對(duì)響應(yīng)的貢獻(xiàn)。

    該方法也適用于計(jì)算子系統(tǒng)A的頻響函數(shù)矩陣??梢钥闯?,公式(10)中存在矩陣求逆運(yùn)算。如果子系統(tǒng)之間的路徑不具有足夠大的彈性,那么HAcAi和HBcAc以及HAcAc和HBcAc將會(huì)十分接近,這樣將導(dǎo)致矩陣求逆時(shí)出現(xiàn)奇異性問題。上一節(jié)中所做出的假設(shè)(2)就是為了保證矩陣求逆時(shí)不會(huì)出現(xiàn)奇異性問題。公式(10)給出了一般情況下計(jì)算子系統(tǒng)頻響函數(shù)的計(jì)算公式,在該公式中包括了所有自由度的信息。實(shí)際上,對(duì)于計(jì)算子系統(tǒng)頻響函數(shù),所有的耦合自由度信息的確是必不可少的,而內(nèi)部自由度信息則并不是必須知道的。如果不考慮某個(gè)子系統(tǒng)的內(nèi)部自由度,則將式(10)中與該內(nèi)部自由度有關(guān)的部分用空矩陣替代即可。比如,如果子系統(tǒng)A和B上的內(nèi)部自由度均不考慮,即XAi和XBi均為空向量。這種情況下,式(10)中所有下角標(biāo)中含有i的項(xiàng)將消失。式(10)將退化成HB=HBB-HBA(HAA-HBA)-1HAB-HBB(11)這與文獻(xiàn)[17]中的結(jié)論完全一致,說明文獻(xiàn)[17]的結(jié)論是本文結(jié)論的一個(gè)特例,即在不考慮任何內(nèi)部自由度的情況下所得出的結(jié)論。而在實(shí)際情況下,通常是要考慮內(nèi)部自由度的。比如,對(duì)由發(fā)動(dòng)機(jī)傳遞至車身某處的振動(dòng)進(jìn)行傳遞路徑分析,作者所感興趣的這個(gè)響應(yīng)點(diǎn)就是一個(gè)內(nèi)部自由度。因此,本文所提出的方法更具有普適性。

    3數(shù)值算例

    為驗(yàn)證上一節(jié)所提出的方法的正確性,建立如圖2所示的8自由度振動(dòng)系統(tǒng)模型。該模型中包含兩個(gè)子系統(tǒng),即主動(dòng)部分子系統(tǒng)A和被動(dòng)部分子系統(tǒng)B。兩個(gè)子系統(tǒng)之間通過3條由線性彈簧和阻尼組成的路徑相連接。子系統(tǒng)A所產(chǎn)生的振動(dòng)通過三條路徑傳遞至子系統(tǒng)B上。模型參數(shù)如表1所示。

    自由度類型所包含的自由度主動(dòng)部分內(nèi)部自由度Aim1主動(dòng)部分耦合自由度Acm2, m3, m4被動(dòng)部分耦合自由度Bcm5, m6, m7被動(dòng)部分內(nèi)部自由度Bim8

    根據(jù)表2所做的自由度類型劃分,可以將式 (13)所計(jì)算得到的系統(tǒng)頻響函數(shù)矩陣表示成如式(3)的分塊矩陣的形式。然后再根據(jù)所得出的分塊運(yùn)用式(10)進(jìn)行計(jì)算,可以求出子系統(tǒng)B的頻響函數(shù)矩陣。將該計(jì)算結(jié)果與直接計(jì)算子系統(tǒng)B的頻響函數(shù)矩陣所得到的結(jié)果進(jìn)行比較,如圖3所示,發(fā)現(xiàn)兩種方法計(jì)算結(jié)果完全一致,證明了所提出的方法的正確性。

    4實(shí)驗(yàn)驗(yàn)證

    為了進(jìn)一步驗(yàn)證本文所提出方法的正確性及其在實(shí)際應(yīng)用中的有效性,通過如圖6(a)所示的實(shí)驗(yàn)進(jìn)行了驗(yàn)證。

    實(shí)驗(yàn)對(duì)象為由兩塊鋼板通過4個(gè)彈簧相互連接而成的系統(tǒng)。下面的鋼板通過4個(gè)彈簧與臺(tái)面相連接。上面的鋼板為子系統(tǒng)A,下面的鋼板為子系統(tǒng)B??偣部紤]9個(gè)自由度,分別是兩塊鋼板與彈簧的連接處共8個(gè)點(diǎn)以及位于下面一塊鋼板上的目標(biāo)響應(yīng)點(diǎn)9(如圖6(b)所示)。首先,測(cè)量整個(gè)系統(tǒng)的9×9的頻響函數(shù)矩陣,然后根據(jù)本文所提出的方法求取下面鋼板的在解耦狀態(tài)下的5×5頻響函數(shù)矩陣,從而實(shí)現(xiàn)子系統(tǒng)的解耦。

    選取5×5頻響函數(shù)矩陣中與TPA有關(guān)的四個(gè)頻響函數(shù)進(jìn)行分析,將其與直接測(cè)量的結(jié)果進(jìn)行比較,如圖7所示。

    從圖7中可以看出,本文方法所得出的結(jié)果與直接測(cè)量的結(jié)果比較接近,進(jìn)一步證明了本文方法的正確性。測(cè)量誤差主要來自于頻響函數(shù)測(cè)量時(shí)傳感器布置及力錘敲擊時(shí)的位置偏差。

    5結(jié)論

    (1) 提出了一種計(jì)算子系統(tǒng)頻響函數(shù)矩陣的新方法。該方法不需要拆除子系統(tǒng),只需要知道整體的頻響函數(shù)矩陣即可推導(dǎo)出子系統(tǒng)的頻響函數(shù)矩陣,大大縮短了測(cè)量子系統(tǒng)頻響函數(shù)所需要的時(shí)間。

    (2) 傳遞路徑分析中被動(dòng)部分子系統(tǒng)的頻響函數(shù)的測(cè)量,可以用本文所提出的解耦方法進(jìn)行計(jì)算。

    (3) 用所提出的方法計(jì)算子系統(tǒng)頻響函數(shù)時(shí),耦合點(diǎn)的信息必須是已知的,即所測(cè)量的整體系統(tǒng)頻響函數(shù)矩陣中必須包括所有的耦合點(diǎn)在內(nèi)。在條件允許的情況下,運(yùn)用非接觸式測(cè)量會(huì)獲得更好的結(jié)果。

    參考文獻(xiàn):

    [1]Klerk D D, Ossipov A. Operational transfer path analysis: theory, guidelines and tire noise application[J]. Mechanical Systems & Signal Processing, 2010, 24(7):1950—1962.

    [2]Gajdatsy P, Janssens K, Desmet W, et al. Application of the transmissibility concept in transfer path analysis[J]. Mechanical Systems & Signal Processing, 2010, 24(7):1963—1976.

    [3]Guasch O, García C, Jové J, et al. Experimental validation of the direct transmissibility approach to classical transfer path analysis on a mechanical setup[J]. Mechanical Systems & Signal Processing, 2013, 37(1-2):353—369.

    [4]Wang Z, Zhu P. Response prediction for modified mechanical systems based on in-situ frequency response functions: Theoretical and numerical studies[J]. Journal of Sound & Vibration, 2017, 400:417—441.

    [5]Lennstrm D, Olsson M, Wullens F, et al. Validation of the blocked force method for various boundary conditions for automotive source characterization[J]. Applied Acoustics, 2016, 102:108—119.

    [6]Klerk D D, Rixen D J. Component transfer path analysis method with compensation for test bench dynamics[J]. Mechanical Systems & Signal Processing, 2010, 24(6):1693—1710.

    [7]Seijs M V V D, Klerk D D, Rixen D J. General framework for transfer path analysis: History, theory and classification of techniques[J]. Mechanical Systems & Signal Processing, 2016, (68-69): 217—244.

    [8]Klerk D D, Rixen D J, Voormeeren S N. General framework for dynamic substructuring: history, review and classification of techniques[J]. AIAA Journal, 2012, 46(5):1169—1181.

    [9]D′Ambrogio W, Sestieri A. A unified approach to substructuring and structural modification problems[J]. Shock & Vibration, 2014, 11(3):295—309.

    [10]Okubo N, Miyazaki M. Development of uncoupling technique and its application[C]. Proceedings of 4th International Modal Analysis Conference, Florida, USA, 1984:1194—1200.

    [11]Okubo N, Miyazaki M. Development of uncoupling technique to extract component′s dynamics from assembly structure[J]. Journal of the Japan Society of Precision Engineering, 1987, 53:711—716.

    [12]Gontier C, Bensaibi M. Time domain identification of a substructure from in situanalysis of the whole structure[J]. Mechanical Systems & Signal Processing, 1995, 9(4):379—396.

    [13]Silva J M M, Maia N M M, Ribeiro P M L. Dynamic modeling: application of uncoupling techniques[C]. Proceedings of the 14th International Modal Analysis Conference, Michigan, USA, 1996:1165—1172.

    [14]Maia N M M, Silva J M M, Ribeiro A M R, et al. On the dynamic characterization of joints using uncoupling techniques[C]. Proceedings of the 16th International Modal Analysis Conference, California, USA, 1998:1132—1138.

    [15]Kalling P, Abrahamsson T, Mckelvey T. Subsystem state-space model identification and its sensitivity to test variability[C]. International Conference on Noise and Vibration Engineering, Leuven, Belgium, 2004: 2729—2744.

    [16]D′Ambrogio W, Fregolent A. Inverse dynamic substructuring using the direct hybrid assembly in the frequency domain[J]. Mechanical Systems & Signal Processing, 2014, 45(2):360—377.

    [17]Keersmaekers L, Mertens L, Penne R, et al. Decoupling of mechanical systems based on in-situ frequency response functions: the link-preserving, decoupling method[J]. Mechanical Systems & Signal Processing, 2015, (58-59): 340—354.An approach to calculate subsystems′ frequency response

    Abstract: Transfer path analysis(TPA) is a widely used and effective method in the field of vibration and noise control. In TPA, vibrational systems are divided into the active part, passive part and some transfer paths through which the active and passive parts are connected. In TPA, it is necessary to measure the frequency response functions(FRFS) of the passive part. In classical TPA, the substructure needs to be disassembled firstly and then the frequency response function is measured. Therefore, the test process is very cumbersome. In this paper, a novel method is proposed to calculate the FRFs of the subsystem. The formulation of the FRF matrix of the subsystem is derived from the FRF matrix of the whole system directly. Obviously, the proposed method does not require the physical decoupling of subsystems. Consequently, the time required to measure the frequency response function of the subsystem is greatly shorten. The correctness and effectiveness of this method are validated by a numerical case and an experimental case.

    Key words: subsystem; transfer path analysis; frequency response function; decoupling

    猜你喜歡
    子系統(tǒng)
    網(wǎng)絡(luò)空間供應(yīng)鏈中入侵檢測(cè)及防御子系統(tǒng)的投資機(jī)制研究
    網(wǎng)絡(luò)空間供應(yīng)鏈中入侵檢測(cè)及防御子系統(tǒng)的投資機(jī)制研究
    高速公路機(jī)電系統(tǒng)組成及子系統(tǒng)研究
    網(wǎng)絡(luò)綜合布線系統(tǒng)設(shè)計(jì)初探
    團(tuán)隊(duì)與知識(shí)管理的關(guān)系研究
    高職院校學(xué)生管理信息系統(tǒng)功能初探
    淺談中職學(xué)校新型模塊化機(jī)房建設(shè)
    開發(fā)實(shí)驗(yàn)室管理系統(tǒng)信息發(fā)布子系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)
    西藏農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)協(xié)調(diào)度動(dòng)態(tài)分析
    IPv6環(huán)境下動(dòng)態(tài)域名系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)
    又黄又爽又免费观看的视频| 久久草成人影院| 久久久精品欧美日韩精品| 最好的美女福利视频网| 九色亚洲精品在线播放| √禁漫天堂资源中文www| 精品久久蜜臀av无| 国产精品 国内视频| 51午夜福利影视在线观看| 国产精品1区2区在线观看.| 午夜成年电影在线免费观看| 久久午夜亚洲精品久久| 亚洲精品一二三| 国产精品免费视频内射| 午夜福利影视在线免费观看| 亚洲欧美激情综合另类| av网站免费在线观看视频| 日韩欧美三级三区| 色哟哟哟哟哟哟| 男女高潮啪啪啪动态图| 国产精品美女特级片免费视频播放器 | 国产精品野战在线观看 | 在线观看免费视频日本深夜| 淫妇啪啪啪对白视频| 校园春色视频在线观看| 欧美中文日本在线观看视频| 国产亚洲欧美在线一区二区| 黄色女人牲交| 手机成人av网站| 亚洲色图综合在线观看| www.熟女人妻精品国产| 校园春色视频在线观看| 免费在线观看黄色视频的| 高清毛片免费观看视频网站 | 久久中文字幕人妻熟女| av欧美777| 妹子高潮喷水视频| 亚洲情色 制服丝袜| 成熟少妇高潮喷水视频| 国产成人一区二区三区免费视频网站| 午夜免费激情av| 男女下面插进去视频免费观看| 久热这里只有精品99| 99久久人妻综合| 身体一侧抽搐| 精品福利观看| 国产精品国产高清国产av| 国产精品1区2区在线观看.| 亚洲午夜理论影院| 50天的宝宝边吃奶边哭怎么回事| 成人国语在线视频| 丝袜人妻中文字幕| 精品免费久久久久久久清纯| 久久久久久久久免费视频了| 日本一区二区免费在线视频| 亚洲色图 男人天堂 中文字幕| 精品国产一区二区久久| 日本撒尿小便嘘嘘汇集6| 视频在线观看一区二区三区| 少妇裸体淫交视频免费看高清 | 一级毛片精品| 亚洲成国产人片在线观看| 久久香蕉激情| 最近最新中文字幕大全免费视频| 免费高清视频大片| 亚洲欧美日韩无卡精品| av天堂在线播放| 国产精品美女特级片免费视频播放器 | 亚洲精品一卡2卡三卡4卡5卡| xxxhd国产人妻xxx| 欧洲精品卡2卡3卡4卡5卡区| x7x7x7水蜜桃| 亚洲五月婷婷丁香| 国产精品自产拍在线观看55亚洲| 黄色女人牲交| 97超级碰碰碰精品色视频在线观看| 国产单亲对白刺激| 高清毛片免费观看视频网站 | 久久香蕉国产精品| 精品国产美女av久久久久小说| 韩国精品一区二区三区| 久久婷婷成人综合色麻豆| 国产精品免费视频内射| 欧美精品啪啪一区二区三区| 99国产精品一区二区三区| 91九色精品人成在线观看| 男男h啪啪无遮挡| www日本在线高清视频| 超碰成人久久| 国产乱人伦免费视频| 9热在线视频观看99| 日韩中文字幕欧美一区二区| 女人被躁到高潮嗷嗷叫费观| 在线观看66精品国产| 中文字幕精品免费在线观看视频| 国产成+人综合+亚洲专区| 精品卡一卡二卡四卡免费| 可以在线观看毛片的网站| 精品第一国产精品| 久久伊人香网站| 91麻豆av在线| 操美女的视频在线观看| 日韩有码中文字幕| 男人操女人黄网站| 一个人观看的视频www高清免费观看 | av免费在线观看网站| 亚洲一区二区三区不卡视频| 美女福利国产在线| 人妻丰满熟妇av一区二区三区| 国产成人精品久久二区二区免费| 亚洲国产欧美网| 成人手机av| 18美女黄网站色大片免费观看| 国产成人精品久久二区二区91| 久9热在线精品视频| 国产精品 国内视频| 视频在线观看一区二区三区| 国产精品日韩av在线免费观看 | 国产亚洲av高清不卡| 欧美性长视频在线观看| 久久精品亚洲熟妇少妇任你| 在线十欧美十亚洲十日本专区| 国产av一区在线观看免费| 91大片在线观看| 成熟少妇高潮喷水视频| 国内久久婷婷六月综合欲色啪| xxx96com| 成熟少妇高潮喷水视频| 97超级碰碰碰精品色视频在线观看| 三级毛片av免费| 精品久久久久久电影网| 日本一区二区免费在线视频| 黄色成人免费大全| av天堂在线播放| 亚洲成人国产一区在线观看| 久久国产精品男人的天堂亚洲| 亚洲美女黄片视频| 国产野战对白在线观看| 久久久久亚洲av毛片大全| 精品电影一区二区在线| 少妇粗大呻吟视频| 国产无遮挡羞羞视频在线观看| 亚洲欧美激情在线| 美女国产高潮福利片在线看| 又黄又爽又免费观看的视频| 岛国在线观看网站| 欧美激情久久久久久爽电影 | 色尼玛亚洲综合影院| 成熟少妇高潮喷水视频| 中文字幕最新亚洲高清| 久久久久久久久免费视频了| 色综合婷婷激情| 搡老熟女国产l中国老女人| 日韩免费高清中文字幕av| 热99国产精品久久久久久7| 日韩精品中文字幕看吧| 韩国av一区二区三区四区| 一本大道久久a久久精品| 国产精品久久久久成人av| 免费日韩欧美在线观看| 精品欧美一区二区三区在线| 性色av乱码一区二区三区2| 欧美黄色片欧美黄色片| 久久国产乱子伦精品免费另类| 91成人精品电影| 久久精品aⅴ一区二区三区四区| 深夜精品福利| 男女下面插进去视频免费观看| 母亲3免费完整高清在线观看| tocl精华| 精品久久久久久电影网| 久久久国产一区二区| 黄色怎么调成土黄色| 99久久99久久久精品蜜桃| 曰老女人黄片| 亚洲专区国产一区二区| 亚洲精品在线美女| 男人舔女人下体高潮全视频| 最好的美女福利视频网| 精品人妻在线不人妻| 一二三四社区在线视频社区8| 嫩草影院精品99| 国产有黄有色有爽视频| 精品久久久久久久毛片微露脸| 天天添夜夜摸| 日韩大码丰满熟妇| 亚洲九九香蕉| 久久99一区二区三区| 日本免费a在线| 午夜福利在线观看吧| 免费看十八禁软件| 欧美日韩av久久| 国产97色在线日韩免费| 欧美黄色淫秽网站| 成人18禁在线播放| 成人18禁高潮啪啪吃奶动态图| 人人妻人人添人人爽欧美一区卜| 真人做人爱边吃奶动态| 美国免费a级毛片| 欧美日韩精品网址| 国产一区二区三区视频了| 免费看a级黄色片| 在线永久观看黄色视频| 一二三四社区在线视频社区8| 精品国产乱码久久久久久男人| 日本免费a在线| av网站免费在线观看视频| 69av精品久久久久久| 真人一进一出gif抽搐免费| 久久久国产成人精品二区 | 18美女黄网站色大片免费观看| 十八禁人妻一区二区| 女人精品久久久久毛片| 人妻丰满熟妇av一区二区三区| 亚洲自偷自拍图片 自拍| 欧美精品亚洲一区二区| 男女午夜视频在线观看| 国产精品自产拍在线观看55亚洲| 精品无人区乱码1区二区| 丰满迷人的少妇在线观看| 曰老女人黄片| av天堂在线播放| 久久中文字幕人妻熟女| 国产国语露脸激情在线看| 国产精品电影一区二区三区| 黄色 视频免费看| 美女大奶头视频| 女生性感内裤真人,穿戴方法视频| 一区在线观看完整版| 级片在线观看| 精品少妇一区二区三区视频日本电影| 亚洲国产精品sss在线观看 | 侵犯人妻中文字幕一二三四区| 国产又爽黄色视频| 男女之事视频高清在线观看| 日本撒尿小便嘘嘘汇集6| 岛国在线观看网站| 免费av毛片视频| 在线观看午夜福利视频| 法律面前人人平等表现在哪些方面| 欧美不卡视频在线免费观看 | 在线观看免费视频网站a站| 男女下面插进去视频免费观看| www.精华液| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久人人做人人爽| 老鸭窝网址在线观看| 国产熟女午夜一区二区三区| 日韩av在线大香蕉| 女人精品久久久久毛片| 成人永久免费在线观看视频| 在线观看一区二区三区| 久久久久九九精品影院| 国产1区2区3区精品| 亚洲国产精品999在线| 电影成人av| 丰满迷人的少妇在线观看| 天堂俺去俺来也www色官网| 亚洲av成人av| 精品国产一区二区久久| 黑人操中国人逼视频| 精品国产一区二区三区四区第35| 久久久久久久久久久久大奶| 50天的宝宝边吃奶边哭怎么回事| 97碰自拍视频| 欧美在线一区亚洲| 黄片播放在线免费| 国产高清视频在线播放一区| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 麻豆久久精品国产亚洲av | www.熟女人妻精品国产| 亚洲三区欧美一区| 国产一区二区三区在线臀色熟女 | 日本 av在线| 亚洲色图 男人天堂 中文字幕| 欧美日韩国产mv在线观看视频| 亚洲精品久久成人aⅴ小说| 欧美人与性动交α欧美精品济南到| 女人爽到高潮嗷嗷叫在线视频| 久久精品aⅴ一区二区三区四区| 久久九九热精品免费| 丝袜美足系列| 国产高清激情床上av| 在线观看日韩欧美| 曰老女人黄片| 久久精品人人爽人人爽视色| 熟女少妇亚洲综合色aaa.| 免费在线观看影片大全网站| 黄频高清免费视频| 黄色丝袜av网址大全| 亚洲人成网站在线播放欧美日韩| 91成年电影在线观看| 黑人猛操日本美女一级片| 午夜免费成人在线视频| 中亚洲国语对白在线视频| 久久人妻av系列| 99久久人妻综合| 免费人成视频x8x8入口观看| 亚洲在线自拍视频| 少妇粗大呻吟视频| 国产伦人伦偷精品视频| 成人亚洲精品av一区二区 | 国产精品影院久久| 中文字幕人妻丝袜制服| 制服人妻中文乱码| 村上凉子中文字幕在线| 日韩人妻精品一区2区三区| 9191精品国产免费久久| 久久精品aⅴ一区二区三区四区| 亚洲一卡2卡3卡4卡5卡精品中文| 又黄又爽又免费观看的视频| 精品人妻在线不人妻| 视频在线观看一区二区三区| 搡老熟女国产l中国老女人| 精品国产美女av久久久久小说| 国产极品粉嫩免费观看在线| 女同久久另类99精品国产91| 黑丝袜美女国产一区| 淫妇啪啪啪对白视频| 久久久国产欧美日韩av| 国产成人精品无人区| 黑人欧美特级aaaaaa片| 精品电影一区二区在线| 久久人妻福利社区极品人妻图片| 国产高清国产精品国产三级| 九色亚洲精品在线播放| 亚洲第一av免费看| 久久精品人人爽人人爽视色| 久久久国产成人免费| 国产精品一区二区在线不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 97超级碰碰碰精品色视频在线观看| 精品久久久久久久久久免费视频 | 国产真人三级小视频在线观看| 俄罗斯特黄特色一大片| 国产成人av教育| 亚洲精品成人av观看孕妇| 大型黄色视频在线免费观看| 深夜精品福利| 桃红色精品国产亚洲av| 国产人伦9x9x在线观看| 麻豆久久精品国产亚洲av | 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 欧美亚洲日本最大视频资源| 国产精品二区激情视频| 国产精品秋霞免费鲁丝片| 成人永久免费在线观看视频| 亚洲成av片中文字幕在线观看| 精品一区二区三区四区五区乱码| 美国免费a级毛片| 国产av在哪里看| 久久久久久久精品吃奶| 男女高潮啪啪啪动态图| 亚洲成人免费电影在线观看| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区| 9热在线视频观看99| 在线观看一区二区三区| 男人的好看免费观看在线视频 | 一级作爱视频免费观看| 亚洲精品国产色婷婷电影| 在线观看66精品国产| 久久天堂一区二区三区四区| 一区二区三区精品91| 国产精品电影一区二区三区| ponron亚洲| 大型av网站在线播放| 91九色精品人成在线观看| 亚洲第一青青草原| 熟女少妇亚洲综合色aaa.| 午夜免费成人在线视频| 又大又爽又粗| 亚洲一区中文字幕在线| 99精国产麻豆久久婷婷| 国产99久久九九免费精品| 欧美黑人精品巨大| 国产精品秋霞免费鲁丝片| 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 免费av中文字幕在线| 亚洲人成网站在线播放欧美日韩| 欧美激情 高清一区二区三区| 老鸭窝网址在线观看| 中文字幕最新亚洲高清| 精品少妇一区二区三区视频日本电影| 午夜免费观看网址| 国产亚洲欧美98| 9热在线视频观看99| 天天添夜夜摸| 精品一区二区三区av网在线观看| 波多野结衣av一区二区av| 嫁个100分男人电影在线观看| 看片在线看免费视频| 午夜福利一区二区在线看| 国产极品粉嫩免费观看在线| 亚洲欧美日韩另类电影网站| 国产精品av久久久久免费| 大码成人一级视频| 手机成人av网站| 国产精品久久久人人做人人爽| 18禁黄网站禁片午夜丰满| 亚洲av成人一区二区三| 乱人伦中国视频| 交换朋友夫妻互换小说| 99精品在免费线老司机午夜| 亚洲美女黄片视频| 免费看十八禁软件| 夜夜躁狠狠躁天天躁| 久久人人97超碰香蕉20202| 亚洲性夜色夜夜综合| 老司机午夜十八禁免费视频| 午夜老司机福利片| 国产成人啪精品午夜网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产欧美一区二区综合| av电影中文网址| 黄色女人牲交| 国产精品秋霞免费鲁丝片| 老熟妇乱子伦视频在线观看| 午夜精品在线福利| 国产单亲对白刺激| 欧美日韩精品网址| 国产午夜精品久久久久久| 午夜久久久在线观看| 老熟妇仑乱视频hdxx| 69av精品久久久久久| 99精品在免费线老司机午夜| av天堂在线播放| 欧美色视频一区免费| 亚洲欧美激情在线| 一个人免费在线观看的高清视频| 中出人妻视频一区二区| 亚洲国产欧美日韩在线播放| 午夜免费观看网址| avwww免费| 中文字幕人妻熟女乱码| 国产精品永久免费网站| 久久精品亚洲精品国产色婷小说| 黄片小视频在线播放| 国产成人啪精品午夜网站| 成人三级黄色视频| 国产麻豆69| 欧美亚洲日本最大视频资源| 日韩高清综合在线| 曰老女人黄片| 欧美日韩乱码在线| 午夜日韩欧美国产| 啦啦啦 在线观看视频| 中文字幕最新亚洲高清| 波多野结衣av一区二区av| 久久久国产成人精品二区 | 国产区一区二久久| 日韩精品免费视频一区二区三区| 精品一区二区三区av网在线观看| 嫩草影院精品99| 日韩精品中文字幕看吧| 在线观看免费日韩欧美大片| www.自偷自拍.com| 国产亚洲欧美精品永久| 黄片小视频在线播放| 黄色毛片三级朝国网站| 少妇裸体淫交视频免费看高清 | 精品国产乱码久久久久久男人| 老熟妇乱子伦视频在线观看| 久久精品亚洲av国产电影网| 波多野结衣高清无吗| 亚洲欧美精品综合久久99| 丰满的人妻完整版| 在线观看免费高清a一片| 中文字幕另类日韩欧美亚洲嫩草| 亚洲五月色婷婷综合| 午夜精品久久久久久毛片777| 婷婷精品国产亚洲av在线| 老司机福利观看| 在线免费观看的www视频| 国产成人免费无遮挡视频| 男女之事视频高清在线观看| 夜夜看夜夜爽夜夜摸 | 国产无遮挡羞羞视频在线观看| 国产麻豆69| 男女下面进入的视频免费午夜 | 69精品国产乱码久久久| 亚洲欧美日韩无卡精品| 国产日韩一区二区三区精品不卡| 久久久久久久午夜电影 | 青草久久国产| 国产三级在线视频| 伦理电影免费视频| 国产成人精品在线电影| 国产精品98久久久久久宅男小说| 午夜精品在线福利| 亚洲国产毛片av蜜桃av| 法律面前人人平等表现在哪些方面| 岛国视频午夜一区免费看| netflix在线观看网站| 老鸭窝网址在线观看| 丝袜在线中文字幕| 中出人妻视频一区二区| 中文字幕人妻丝袜一区二区| 午夜福利影视在线免费观看| 1024视频免费在线观看| 国产深夜福利视频在线观看| 午夜久久久在线观看| 在线天堂中文资源库| 国产精品乱码一区二三区的特点 | 成熟少妇高潮喷水视频| 免费高清视频大片| 高清欧美精品videossex| 久久久久九九精品影院| 最好的美女福利视频网| av天堂久久9| 精品久久久久久久毛片微露脸| 97人妻天天添夜夜摸| 欧美日韩福利视频一区二区| avwww免费| 国产激情欧美一区二区| 满18在线观看网站| 老司机靠b影院| 欧美一级毛片孕妇| 欧美日本中文国产一区发布| 亚洲一区中文字幕在线| 成年人黄色毛片网站| 老熟妇仑乱视频hdxx| 欧美日韩福利视频一区二区| 国产高清视频在线播放一区| 99在线视频只有这里精品首页| 久久久久久人人人人人| 国产色视频综合| 一区二区三区激情视频| 18美女黄网站色大片免费观看| 久久国产精品影院| 久久天躁狠狠躁夜夜2o2o| 一二三四在线观看免费中文在| 1024视频免费在线观看| 中文亚洲av片在线观看爽| 黄片大片在线免费观看| 999久久久精品免费观看国产| 欧美日韩视频精品一区| 亚洲国产欧美日韩在线播放| 精品熟女少妇八av免费久了| 好男人电影高清在线观看| 国产精品秋霞免费鲁丝片| 久久人人精品亚洲av| 亚洲第一欧美日韩一区二区三区| 巨乳人妻的诱惑在线观看| 人妻久久中文字幕网| av视频免费观看在线观看| 国产日韩一区二区三区精品不卡| 免费高清视频大片| 男男h啪啪无遮挡| 欧美日韩亚洲综合一区二区三区_| 亚洲专区字幕在线| 国产黄色免费在线视频| 国产男靠女视频免费网站| 丰满迷人的少妇在线观看| 国产亚洲精品久久久久5区| 精品国产亚洲在线| 精品日产1卡2卡| 美女高潮到喷水免费观看| 亚洲欧美精品综合一区二区三区| 日本撒尿小便嘘嘘汇集6| 在线免费观看的www视频| 国产精品久久久久成人av| 久9热在线精品视频| av欧美777| 色综合婷婷激情| 午夜免费激情av| 美女大奶头视频| 一级毛片精品| 国产黄色免费在线视频| 亚洲欧美日韩无卡精品| 免费在线观看日本一区| 国产精品 国内视频| 国产男靠女视频免费网站| 淫妇啪啪啪对白视频| 久久久久久久精品吃奶| www.自偷自拍.com| 久久人妻av系列| 狠狠狠狠99中文字幕| 大型黄色视频在线免费观看| 搡老乐熟女国产| 欧美日韩乱码在线| 久久久久久人人人人人| 国产片内射在线| 少妇裸体淫交视频免费看高清 | 午夜福利在线观看吧| 涩涩av久久男人的天堂| 免费在线观看日本一区| 一级,二级,三级黄色视频| 亚洲成人精品中文字幕电影 | 亚洲成人国产一区在线观看| 国产免费男女视频| 免费在线观看黄色视频的| 免费av毛片视频| 欧美 亚洲 国产 日韩一| 免费在线观看黄色视频的| 1024香蕉在线观看| www日本在线高清视频| 日本一区二区免费在线视频| a级片在线免费高清观看视频| 久久精品亚洲熟妇少妇任你| 国产成人av教育| 美女高潮喷水抽搐中文字幕| 手机成人av网站| 亚洲成a人片在线一区二区| 天堂√8在线中文| 变态另类成人亚洲欧美熟女 | 又黄又爽又免费观看的视频| 久久这里只有精品19| 中文字幕人妻丝袜一区二区| 九色亚洲精品在线播放| 国产成人欧美在线观看| 亚洲中文av在线|