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

    加肋圓柱體水下低頻輻射聲特性研究

    2015-04-25 10:11:00祁立波鄒明松中國船舶科學(xué)研究中心江蘇無錫214082
    船舶力學(xué) 2015年7期
    關(guān)鍵詞:圓柱體聲源軸向

    祁立波,鄒明松(中國船舶科學(xué)研究中心,江蘇 無錫 214082)

    加肋圓柱體水下低頻輻射聲特性研究

    祁立波,鄒明松
    (中國船舶科學(xué)研究中心,江蘇 無錫 214082)

    為掌握不同激勵下加肋圓柱體水下低頻輻射聲頻譜特征,揭示輻射聲峰值對應(yīng)的優(yōu)勢模態(tài)及彎縱耦合現(xiàn)象規(guī)律,該文利用三維水彈性聲輻射計算方法,針對加肋圓柱體的梁模型和有限元?dú)つP?,計算分析其水下不同激勵的聲源級傳遞函數(shù),分析聲源級曲線峰值與單個模態(tài)聲源級分量的對應(yīng)關(guān)系,得出垂向激勵時,輻射聲峰值對應(yīng)于垂向彎曲模態(tài);軸向激勵時,輻射聲峰值對應(yīng)于軸向模態(tài),同時運(yùn)用殼模型計算結(jié)果發(fā)現(xiàn)梁模型所不能反映的彎縱耦合現(xiàn)象。比較了不同激勵輻射聲源級傳遞函數(shù),得出所分析頻帶內(nèi),軸向單位力激勵總聲級比垂向略大。并比較不同模型的結(jié)果,指出梁模型的適用范圍。

    輻射聲;彎縱耦合;優(yōu)勢模態(tài)

    0 引 言

    隨著加肋圓柱體在民用、軍事領(lǐng)域的應(yīng)用價值越來越大,其水下噪聲也得到越來越多的關(guān)注和研究。李冰茹等[1]以水下簡支薄殼板為研究對象,討論了殼體的縱橫比以及材質(zhì)對有限長圓柱殼體聲輻射特性的影響,但文中的研究對象簡單,與實(shí)際加肋圓柱體結(jié)構(gòu)差異較大;彭旭等[2]從薄殼理論和Helmholtz波動方程出發(fā),根據(jù)模態(tài)疊加原理分析了有限長加肋圓柱殼內(nèi)點(diǎn)力和點(diǎn)聲源作用的聲輻射特性及傳遞損失;曾革委[3]采用Donnel殼體理論研究結(jié)構(gòu)在環(huán)頻率以下的聲輻射特性,激振力為任意方向的簡諧集中力。但上述研究工作并未針對實(shí)際加肋圓柱體結(jié)構(gòu)運(yùn)行工況中低頻輻射聲特性加以研究,輻射聲峰值與對應(yīng)實(shí)際結(jié)構(gòu)優(yōu)勢模態(tài)的對應(yīng)關(guān)系也缺乏認(rèn)識,對于彎縱耦合現(xiàn)象更是從未提及。

    吳有生[4]將三維適航性理論與結(jié)構(gòu)動力學(xué)理論相結(jié)合,提出了廣義流固界面條件及圍繞彈性體的三維勢流理論,從而形成了一個適用于分析在水下運(yùn)動的任意三維彈性體,承受內(nèi)外激勵時的動響應(yīng)性能的三維線性水彈性力學(xué)理論,其中水處理為不可壓介質(zhì)。鄒明松[5]在三維線性水彈性力學(xué)理論基礎(chǔ)上,采用可壓流體的Green函數(shù),得到適用于分析彈性體結(jié)構(gòu)在水下受激勵時濕表面振動產(chǎn)生的近場點(diǎn)聲壓和遠(yuǎn)場輻射聲的三維水彈性聲輻射計算方法。同時對方法進(jìn)行了驗(yàn)證,并對非規(guī)則頻率作了處理。

    低頻輻射聲主要由推進(jìn)激勵引起,本文利用零航速三維水彈性聲輻射計算方法,針對水下加肋圓柱體梁模型和有限元?dú)つP停M(jìn)行了不同單位力激勵輻射聲傳遞函數(shù)計算分析,獲得結(jié)構(gòu)在水下輻射聲的頻譜特性及決定聲源級峰值的優(yōu)勢模態(tài),同時通過比較兩種模型的計算結(jié)果,指出梁模型的適用范圍。本文結(jié)論有利于提高減振降噪措施的針對性和有效性[6],對低噪聲水下加肋圓柱體設(shè)計有著重要指導(dǎo)意義。

    1 加肋圓柱體水下輻射聲頻域分析方法

    在加肋圓柱體作微幅振動假定的條件下,結(jié)構(gòu)相對其平衡位置的運(yùn)動可采用模態(tài)疊加的表達(dá)形式:

    式中:qr(r=1,2,…,m )為相應(yīng)于第r階干模態(tài)位移ur的主坐標(biāo)分量。

    將水處理成可壓理想聲介質(zhì)。不考慮航速時,根據(jù)線性系統(tǒng)疊加原理可知:流場速度勢Φ為各階聲波輻射速度勢φr的線性疊加。存在如下關(guān)系:

    其中:O(x,y,z),為場點(diǎn)的坐標(biāo)。

    在結(jié)構(gòu)平均濕表面上,零航速的流固耦合邊界條件是:

    式中:ur,vr,wr為相應(yīng)于第r階干模態(tài)位移,n→為加肋圓柱體濕表面單位法向量。零航速時,輻射速度勢自由液面聲學(xué)邊界條件為:

    則與自由液面邊界條件對應(yīng)的Green函數(shù)為[7]:

    式中:r1=為流體中聲波波數(shù),O(x,y,z),為場點(diǎn)P的坐標(biāo),(ξ,η,ζ)為源點(diǎn)Q的坐標(biāo)。

    聲介質(zhì)中水彈性力學(xué)運(yùn)動方程可表示為如下形式:

    式中:{q}為廣義位移向量,{Ξ }為廣義力向量,[a],[b],[c]為干結(jié)構(gòu)廣義質(zhì)量、阻尼和剛度矩陣;[A],[B],[C]為廣義流體附加質(zhì)量、附加阻尼和恢復(fù)力系數(shù)矩陣,其元素可表示為:

    水中輻射聲壓可表示為:

    2 加肋圓柱體梁模型

    根據(jù)加肋圓柱體結(jié)構(gòu)剛度、質(zhì)量分布信息,建立分段非均勻等直梁模型與對應(yīng)的船體濕表面如圖1所示。在低頻段,可以假設(shè)艦船以整體振動為主,橫剖面內(nèi)不發(fā)生變形,即橫剖面周向第一階固有頻率前。同一位置的梁單元對應(yīng)一周濕面元,振動相同,再由濕面元計算聲場。本文考慮推進(jìn)激勵的作用點(diǎn)位置,推進(jìn)激勵本身的頻譜特性不作研究。利用零航速三維水彈性聲輻射計算方法分別求解不同激勵工況[8],得到梁模型聲源級傳遞函數(shù)頻譜曲線。便于分析,本文頻率作無量綱化處理可表示為

    式中:f為頻率,R為加肋圓柱體半徑,g為重力加速度,Ω為無量綱頻率。

    圖1 等效梁與濕表面模型示意圖Fig.1 Model and wet surface

    2.1 梁模型垂向單位力激勵

    梁模型尾部垂向激勵工況如圖2所示。梁模型尾部垂向單位力激勵輻射聲源級曲線如圖3所示,從圖中可以看出曲線的峰值出現(xiàn)在1.2,1.9,2.7,3.6,4.6,5.7。分析每個峰值相關(guān)聯(lián)的優(yōu)勢模態(tài),結(jié)果如圖4所示??紤]到篇幅,只給出第2階和第4階垂向彎曲模態(tài)。

    可以得出,梁模型尾部垂向單位力激勵時,第2~7階垂向彎曲模態(tài)聲源級分量依次對應(yīng)于聲源級曲線的對應(yīng)峰值,不存在彎縱耦合現(xiàn)象。梁模型的第2~7階垂向彎曲干固有頻率分別為1.4,2.3,3.2,4.4,5.5,6.6,對應(yīng)于濕固有頻率1.2,1.9,2.7,3.6,4.6,5.7。

    圖2 梁模型尾部垂向單位力激勵示意圖Fig.2 Beam model under unit vertical force exciting on tail

    圖3 梁模型尾部垂向單位力激勵輻射聲源級曲線圖Fig.3 Source level of beam model under unit vertical force exciting on tail

    圖4 梁模型垂向單位力激勵聲源級與垂向彎曲模態(tài)聲源級分量曲線圖Fig.4 Total source level and components of vertical bending modes

    2.2 梁模型軸向單位力和水平單位力矩(1-2平面)激勵

    梁模型軸向單位力和水平單位力矩激勵工況如圖5所示。這里水平單位力矩是考慮到梁模型軸向力激勵位置不在中心時,對中心的力矩作用。梁模型軸向單位力和水平單位力矩激勵輻射聲源級頻譜曲線如圖6所示,從圖中可以看出曲線的峰值出現(xiàn)在1.2,1.9,2.5,4.3。分析每個峰值相關(guān)聯(lián)的優(yōu)勢模態(tài),結(jié)果如圖7所示。

    圖5 軸向單位力和水平單位力矩激勵示意圖Fig.5 Beam model under unit longitudinal force and unit horizontal moment exciting

    圖6 梁模型軸向單位力和水平單位力矩激勵輻射聲源級曲線圖Fig.6 Source level of beam model under unit longitudinal force and unit horizontal moment exciting

    圖7 梁模型軸向單位力和水平單位力矩激勵聲源級與垂向彎曲模態(tài)和軸向模態(tài)聲源級分量曲線圖Fig.7 Total source level and components of different modes

    由于水平單位力矩作用,梁模型第2,3階垂向彎曲模態(tài)聲源級分量分別對應(yīng)于聲源級曲線的前兩個次要峰值;其次軸向單位力作用,梁模型第1,2階軸向模態(tài)聲源級分量分別對應(yīng)于聲源級曲線的兩個主要峰值,不存在彎縱耦合現(xiàn)象。干固有頻率為1.4,2.3,2.6,4.5,對應(yīng)于濕固有頻率1.2,1.9,2.5,4.3。

    2.3 梁模型兩種激勵輻射聲源級曲線比較

    從圖8中可以看出梁模型垂向單位力激勵聲源級曲線峰值較多,但峰值都較小。由聲源級頻譜曲線計算得出:梁模型推力軸承處軸向單位力和水平單位力矩激勵工況輻射總聲壓級比尾部垂向單位力激勵大4 dB左右。

    圖8 梁模型兩種激勵輻射聲源級曲線比較圖Fig.8 Comparison of source levels of beam model under different exciting

    3 水下加肋圓柱體殼模型

    圖9 整體建模效果圖Fig.9 The finite element shell model

    圖10 殼模型尾部垂向單位力激勵示意圖Fig.10 Shell model under unit vertical force exciting on tail

    水下加肋圓柱體是個非常復(fù)雜的結(jié)構(gòu),其中包括有耐壓結(jié)構(gòu)、各種結(jié)構(gòu)增強(qiáng)構(gòu)件、設(shè)備等各種元素。在建立有限元模型時,不可能完全按照其實(shí)際結(jié)構(gòu)去建模。本文中采取質(zhì)量平均布置,模擬主要結(jié)構(gòu)增強(qiáng)構(gòu)件,再與水下加肋圓柱體整體梁模型模態(tài)分析結(jié)果比較,得到既能準(zhǔn)確表達(dá)水下加肋圓柱體低頻振動特性,又能夠應(yīng)用于水彈性噪聲計算的簡化殼模型[9],如圖9所示。這兒仍然采用與梁模型同樣的濕表面,此時濕表面與同一位置的殼結(jié)構(gòu)單元對應(yīng)。

    3.1 殼模型垂向單位力激勵

    殼模型尾部垂向單位力激勵工況與梁模型相同如圖10所示。殼模型尾部垂向單位力激勵輻射聲源級頻譜曲線如圖11所示,從圖中可以看出曲線的峰值出現(xiàn)在1.3,2.1,2.5,2.9,3.8,4.6,5.4。分析每個峰值相關(guān)聯(lián)的優(yōu)勢模態(tài),如圖12所示??紤]到篇幅,只給出第2,5階垂向彎曲模態(tài)。

    圖11 殼模型尾部垂向單位力激勵輻射聲源級曲線圖Fig.11 Source level of shell model under unit vertical force exciting on tail

    圖12 殼模型尾部垂向單位力激勵聲源級與垂向彎曲模態(tài)聲源級分量曲線圖Fig.12 Total source level and components of vertical bending modes

    圖13 殼模型尾部垂向單位力激勵聲源級與軸向模態(tài)聲源級分量曲線和對應(yīng)模態(tài)振型圖Fig.13 Source level components of longitudinal modes and the corresponding modes of shell model

    我們可以看出殼模型第2~7階垂向彎曲模態(tài)聲源級分量對應(yīng)于聲源級曲線的峰值。可得出垂向激勵時,殼模型整體垂向彎曲模態(tài)對聲源級曲線的對應(yīng)峰值有主要貢獻(xiàn),干固有頻率為1.5,2.4,3.4,4.5,5.4,6.3,對應(yīng)于濕固有頻率峰值點(diǎn)1.3,2.1,2.9,3.8,4.6,5.4??疾檩S向模態(tài)聲源級分量貢獻(xiàn)如圖13所示。

    圖13(a),(b)中可以看出殼模型第1階軸向模態(tài)聲源級分量對應(yīng)于聲源級曲線第二、三個峰值。從振型圖中,我們可以看出,由于質(zhì)量中心與剛度中心的不重合,軸向模態(tài)與垂向彎曲模態(tài)出現(xiàn)耦合現(xiàn)象,因此軸向模態(tài)引起的聲輻射在彎曲模態(tài)頻率點(diǎn)上也有較大貢獻(xiàn)。進(jìn)一步可以看出,1階軸向模態(tài)主要和3階垂向彎曲模態(tài)耦合??傻贸鰵つP偷?階軸向模態(tài)對聲源級曲線的第二、三個峰值都有主要貢獻(xiàn),第1階軸向模態(tài)干固有頻率為2.5,對應(yīng)于濕固有頻率2.4和第3階垂向彎曲模態(tài)對應(yīng)濕頻率2.1。

    圖13(c),(d)中可以看出殼模型第2階軸向模態(tài)聲源級分量無對應(yīng)峰值。從振型圖中看出,2階軸向模態(tài)并未明顯地與彎曲模態(tài)發(fā)生耦合現(xiàn)象。

    3.2 殼模型軸向單位力激勵

    上飛院在上世紀(jì)末設(shè)計并完成了“大攻角多媒體氣動數(shù)據(jù)庫”,初步驗(yàn)證了氣動數(shù)據(jù)庫在軍機(jī)和民機(jī)的研究設(shè)計過程中可以起到的重要作用。該數(shù)據(jù)庫是中國航空工業(yè)總公司第一個包含大量各類媒體動態(tài)變化的氣動研究試驗(yàn)信息,并且便于查詢管理、實(shí)用性很強(qiáng)的多媒體數(shù)據(jù)庫[3]。

    殼模型軸向單位力激勵工況如圖14所示。殼模型軸向單位力激勵輻射聲源級頻譜曲線如圖15所示,從圖中可以看出曲線有兩個主要峰值出現(xiàn)在2.4和4.7,另外兩個次要峰值出現(xiàn)在3.8和5.4。分析峰值相關(guān)聯(lián)的優(yōu)勢模態(tài),如圖16所示。

    圖14 殼模型軸向單位力激勵示意圖Fig.14 Shell model under unit longitudinal force exciting

    圖15 殼模型軸向單位力激勵輻射聲源級曲線圖Fig.15 Source level of shell model under unit longitudinal force exciting

    圖16殼模型軸向單位力激勵聲源級與軸向模態(tài)聲源級分量曲線圖Fig.16 Total source level and components of longitudinal modes

    圖16 中可以看出,軸向單位力作用時,殼模型第1,2階軸向模態(tài)聲源級分量對應(yīng)于聲源級曲線的兩個主要峰值。干固有頻率分別為2.5,4.8,對應(yīng)于濕固有頻率2.4,4.7處曲線的主要峰值。

    圖17殼模型軸向單位力激勵聲源級與垂向彎曲模態(tài)聲源級分量曲線圖Fig.17 Total source level and components of vertical bending modes

    圖17 中可以看出,軸向單位力作用時,殼模型第5,7階垂向彎曲模態(tài)聲源級分量對應(yīng)于聲源級曲線的兩個次要峰值,同樣是因?yàn)槟P偷馁|(zhì)量中心與剛度中心的不重合,軸向激勵時,模型的垂向彎曲模態(tài)也被激發(fā)??傻贸鰵つP偷?,7階整體垂向彎曲模態(tài)干固有頻率分別為4.5,6.3,對應(yīng)于濕固有頻率3.8,5.4處曲線的次要峰值。

    3.3 殼模型兩種激勵輻射聲源級曲線比較

    圖18中可以看出殼模型垂向單位力激勵聲源級曲線峰值同樣較多,但峰值較小。由聲源級頻譜曲線計算得出:殼模型軸向單位力激勵工況輻射總聲壓級比垂向單位力激勵大2 dB左右。

    圖18 殼模型兩種激勵輻射聲源級曲線比較圖Fig.18 Comparison of source levels of shell model under different exciting

    4 兩種模型比較

    梁模型與殼模型頻率比較見表1。

    表1 梁模型與殼模型干濕頻率比較Tab.1 Frequencies of beam and shell model

    從表中可以看出梁模型與殼模型濕頻率誤差大于干頻率,另外由于彎縱耦合的影響,梁模型三階彎曲濕頻率有了較大誤差。梁模型與殼模型聲源級曲線比較如圖19,20。從圖19可以看出同樣因?yàn)閺澘v耦合的影響,梁模型在三階彎曲與一階縱向濕頻率處存在較大誤差。從圖20可以看出,隨著頻率的升高,梁模型的誤差逐漸增大。

    5 結(jié) 論

    (2)水下加肋圓柱體殼模型在尾部垂向單位力激勵時,輻射聲主要由整體彎曲模態(tài)和軸向模態(tài)引起。聲源級曲線在垂向彎曲模態(tài)濕頻率點(diǎn)處存在峰值。同時在軸向模態(tài)濕頻率點(diǎn)處也存在峰值。存在彎縱耦合現(xiàn)象。

    (3)水下加肋圓柱體殼模型在軸向單位力激勵時,輻射聲主要由軸向模態(tài)引起,聲源級曲線在軸向模態(tài)濕頻率點(diǎn)處存在主要峰值,同時在相應(yīng)彎曲模態(tài)濕頻率點(diǎn)處也存在次要峰值。存在彎縱耦合現(xiàn)象。

    (4)水下加肋圓柱體軸向單位激勵輻射總聲壓級略大于垂向單位激勵。

    圖19 梁模型與殼模型尾部垂向激勵輻射聲源級曲線比較圖Fig.19 Comparison of source levels of shell model and beam model under different exciting

    圖20 梁模型與殼模型軸向激勵輻射聲源級曲線比較圖Fig.20 Comparison of source levels of shell model and beam model under different exciting

    (5)梁模型適用于無彎縱耦合模態(tài)的結(jié)構(gòu),在彎縱耦合影響頻率附近會產(chǎn)生較大誤差。隨著頻率的升高,梁模型誤差逐漸增大。

    [1]李冰茹,王宣銀,葛輝良.有限長圓柱殼體聲輻射特性影響因素研究[J].浙江大學(xué)學(xué)報(工學(xué)版),2007,41(5):785-789. Li Bingru,Wang Xuanyin,Ge huiliang.Study on influencing factors for sound radiation characteristics of finite length cylindrical shells[J].Journal of Zhejiang University(Engineering Science),2007,41(5):785-789.

    [2]彭 旭,關(guān)珊珊,駱東平,陳美霞.內(nèi)部聲激勵下加筋圓柱殼的聲輻射特性分析[J].固體力學(xué)學(xué)報,2007,28(4):355-361. Peng Xu,Guan Shanshan,Luo Dongping,Chen Meixia.Sound radiation analysis of stiffened cylindrical shells excited by interior sound excitation[J].Acta Mechanica Solida Sinica,2007,28(4):355-361.

    [3]曾革委.潛艇結(jié)構(gòu)輻射噪聲的建模、求解及其聲特性研究[D].武漢:華中科技大學(xué),2002,10. Zeng Gewei.Research of model building and solving of acoustic radiation of submarine[D].Wuhan:Huazhong University of Science and Technology,2002,10.

    [4]Wu Yousheng.Hydroelasticity of floating bodies[D].London:Brunel University,1984.

    [5]Zou Mingsong,Wu Yousheng,Ye Yonglin.Threee-dimensional hydroelasticity analysis of acoustic responses of ship structures[C].9th International Conference on Hydrodynamics,2010,10:844-851.

    [6]吳經(jīng)彪,姜 哲,朱利鋒.基于聲輻射模態(tài)的有源控制解耦[J].聲學(xué)學(xué)報,2009,34(5):453-461.Wu Jingbiao,Jiang Zhe,Zhu Lifeng.The decoupling of the active structural acoustic control based on the theory of radiation modes[J].Acta Acustica,2009,34(5):453-461.

    [7]鄒明松,吳有生,沈順根,吳文偉.考慮航速及自由液面影響的聲介質(zhì)中三維結(jié)構(gòu)水彈性力學(xué)研究[J].船舶力學(xué), 2010,14(11):1304-1311. Zou Mingsong,Wu Yousheng,Shen Shungen,Wu Wenwei.Three-dimensional hydroelasticity with forward speed and free surface in acoustic medium[J].Journal of Ship Mechanics,2010,14(11):1304-1311.

    [8]謝基榕.螺旋槳、軸系激勵艇體結(jié)構(gòu)引起水下噪聲研究[D].無錫:中國船舶科學(xué)研究中心,2003. Xie Jirong.Propeller excitated noise radiation of ship hull structure[D].Wuxi:China Ship Scientific Research Center, 2003.

    [9]曹志遠(yuǎn).板殼振動理論[M].北京:中國鐵道出版社,1989. Cao Zhiyuan.Vibration theory of plate and shell[M].Beijing:China Railway Press,1989.

    Research on acoustic radiation of stiffened cylinder in low-frequency

    QI Li-bo,ZOU Ming-song
    (China Ship Scientific Research Center,Wuxi 214082,China)

    To obtain the character of acoustic radiation of stiffened cylinder in low-frequency,and disclose the dominant modes corresponding to the peaks of acoustic radiation and the phenomenon of the coupled mode of bending and longitudinal deformation,three-dimensional sono-elastical method is applied to solve the transfer function of acoustic radiation of the beam model and shell model of stiffened cylinder under different exciting.The vertical bending modes corresponding to the peaks under vertical exciting and the longitudinal modes corresponding to the peaks under longitudinal exciting are obtained by using the relation between the source level curve and the component of a single mode.The coupled mode of both bending and longitudinal deformation of the shell model is also obtained,which could not be found of the beam model.By comparing the transfer function of source level under different exciting,the total acoustic power under longitudinal exciting is slightly higher than that under vertical exciting is obtained.The applicability of beam model is identified by comparing the results of shell model and beam model.

    acoustic radiation;coupled mode of bending and longitudinal deformation;dominant mode

    TB53

    A

    10.3969/j.issn.1007-7294.2015.07.014

    1007-7294(2015)07-0874-10

    2015-05-07

    祁立波(1985-),男,高級工程師,博士研究生,E-mail:qilibo1984@163.com;

    鄒明松(1982-),男,高級工程師,博士研究生。

    猜你喜歡
    圓柱體聲源軸向
    虛擬聲源定位的等效源近場聲全息算法
    大型立式單級引黃離心泵軸向力平衡的研究
    基于GCC-nearest時延估計的室內(nèi)聲源定位
    電子制作(2019年23期)2019-02-23 13:21:12
    荒銑加工軸向切深識別方法
    運(yùn)用內(nèi)積相關(guān)性結(jié)合迭代相減識別兩點(diǎn)聲源
    找出圓柱體
    娃娃畫報(2016年6期)2016-07-29 19:31:34
    圓柱體上的最短路徑
    微小型薄底零件的軸向車銑實(shí)驗(yàn)研究
    力-聲互易在水下聲源強(qiáng)度測量中的應(yīng)用
    特殊多線螺紋的快速軸向分線法
    河南科技(2014年7期)2014-02-27 14:11:18
    亚洲情色 制服丝袜| 欧美午夜高清在线| 在线看a的网站| 国产成人av教育| 淫妇啪啪啪对白视频| 欧美日韩精品网址| 日韩欧美一区视频在线观看| а√天堂www在线а√下载 | 老熟妇仑乱视频hdxx| 国产精品久久久人人做人人爽| а√天堂www在线а√下载 | 国产欧美日韩精品亚洲av| 丝袜人妻中文字幕| 最近最新中文字幕大全电影3 | 丝袜美腿诱惑在线| 老司机午夜福利在线观看视频| 一区二区三区国产精品乱码| 免费看a级黄色片| 国产免费现黄频在线看| 一级作爱视频免费观看| 久久影院123| 国产精品一区二区精品视频观看| 欧美午夜高清在线| 午夜精品在线福利| 看片在线看免费视频| 国产一区二区三区视频了| 国产不卡av网站在线观看| 国产精品久久久久久人妻精品电影| 欧美黄色淫秽网站| 波多野结衣一区麻豆| 国产高清国产精品国产三级| 欧美激情久久久久久爽电影 | 国产精品影院久久| 91av网站免费观看| 美女扒开内裤让男人捅视频| 久久国产精品影院| 日本精品一区二区三区蜜桃| 操出白浆在线播放| 韩国av一区二区三区四区| 黑人巨大精品欧美一区二区蜜桃| 国产免费现黄频在线看| 18禁裸乳无遮挡动漫免费视频| 精品国产国语对白av| 丝瓜视频免费看黄片| 女人久久www免费人成看片| 久久久国产成人免费| 亚洲精品久久午夜乱码| 日韩欧美国产一区二区入口| 精品国产乱码久久久久久男人| 操美女的视频在线观看| 757午夜福利合集在线观看| 一级a爱片免费观看的视频| 久久人妻熟女aⅴ| 久久久精品国产亚洲av高清涩受| 国产在视频线精品| 亚洲国产精品sss在线观看 | 久久久水蜜桃国产精品网| 一个人免费在线观看的高清视频| 操美女的视频在线观看| 午夜激情av网站| av欧美777| 啪啪无遮挡十八禁网站| 91大片在线观看| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜一区二区| 精品国产一区二区久久| 欧美人与性动交α欧美软件| 久久精品国产亚洲av香蕉五月 | 国产精品综合久久久久久久免费 | 黑人巨大精品欧美一区二区蜜桃| а√天堂www在线а√下载 | 天天添夜夜摸| 亚洲综合色网址| 757午夜福利合集在线观看| 国产又色又爽无遮挡免费看| 日韩人妻精品一区2区三区| 免费av中文字幕在线| 一二三四在线观看免费中文在| 丰满的人妻完整版| 午夜精品在线福利| 精品久久久精品久久久| 精品久久久久久,| 久久国产精品男人的天堂亚洲| 国产精品免费大片| 18禁裸乳无遮挡动漫免费视频| 欧美中文综合在线视频| 亚洲一区二区三区不卡视频| 老司机靠b影院| 免费在线观看完整版高清| 最近最新中文字幕大全免费视频| bbb黄色大片| 国产在线精品亚洲第一网站| 亚洲午夜精品一区,二区,三区| 国产在线一区二区三区精| 久久狼人影院| 1024香蕉在线观看| 久久99一区二区三区| 成人国语在线视频| 亚洲中文日韩欧美视频| 12—13女人毛片做爰片一| 1024香蕉在线观看| 大陆偷拍与自拍| 一区二区三区精品91| 美女扒开内裤让男人捅视频| 日韩免费高清中文字幕av| 91字幕亚洲| 亚洲第一青青草原| 黄色视频,在线免费观看| 国产精品美女特级片免费视频播放器 | 黄片小视频在线播放| 久久精品人人爽人人爽视色| 大香蕉久久成人网| 欧美黄色淫秽网站| 国产成人啪精品午夜网站| 真人做人爱边吃奶动态| 美国免费a级毛片| 国产精品影院久久| 国产精品九九99| 亚洲国产欧美一区二区综合| 两人在一起打扑克的视频| 亚洲成a人片在线一区二区| 黑人猛操日本美女一级片| 窝窝影院91人妻| 悠悠久久av| 精品一区二区三卡| 久久精品国产亚洲av香蕉五月 | 久久久久国产一级毛片高清牌| 伦理电影免费视频| 欧美日韩精品网址| 婷婷成人精品国产| 99国产极品粉嫩在线观看| 欧美激情极品国产一区二区三区| 欧美在线一区亚洲| 国产精品亚洲av一区麻豆| 欧美最黄视频在线播放免费 | 后天国语完整版免费观看| 成人三级做爰电影| 久久人人爽av亚洲精品天堂| 精品国产国语对白av| 亚洲色图av天堂| 12—13女人毛片做爰片一| 1024香蕉在线观看| 交换朋友夫妻互换小说| 成人国产一区最新在线观看| 亚洲精品一二三| 男女之事视频高清在线观看| 99久久国产精品久久久| xxxhd国产人妻xxx| 国产精品免费视频内射| tube8黄色片| 日韩三级视频一区二区三区| 免费看十八禁软件| 久久久国产成人精品二区 | 欧美丝袜亚洲另类 | 正在播放国产对白刺激| 搡老熟女国产l中国老女人| 少妇 在线观看| 12—13女人毛片做爰片一| 久久精品aⅴ一区二区三区四区| 亚洲精品一二三| 日韩三级视频一区二区三区| 国产午夜精品久久久久久| 久久久久久久久久久久大奶| 亚洲人成77777在线视频| 欧美色视频一区免费| 国产99久久九九免费精品| 亚洲成人免费电影在线观看| 国产又爽黄色视频| 欧美激情极品国产一区二区三区| 国产极品粉嫩免费观看在线| av天堂久久9| 不卡av一区二区三区| 丁香欧美五月| 欧美日本中文国产一区发布| 久久久久久久精品吃奶| 午夜老司机福利片| 99精品在免费线老司机午夜| 一边摸一边做爽爽视频免费| 精品久久久久久,| 满18在线观看网站| 高清在线国产一区| 亚洲精品成人av观看孕妇| 在线视频色国产色| 亚洲av成人av| 青草久久国产| √禁漫天堂资源中文www| 欧美日韩av久久| 亚洲av成人不卡在线观看播放网| 国产av一区二区精品久久| 国产伦人伦偷精品视频| 一级作爱视频免费观看| 国产精品电影一区二区三区 | 啪啪无遮挡十八禁网站| 国产麻豆69| 国产在线精品亚洲第一网站| 亚洲午夜精品一区,二区,三区| 男人舔女人的私密视频| 十八禁网站免费在线| 久久精品成人免费网站| 在线av久久热| 黑人操中国人逼视频| 亚洲五月色婷婷综合| 80岁老熟妇乱子伦牲交| 日韩视频一区二区在线观看| 欧美人与性动交α欧美精品济南到| 国精品久久久久久国模美| 欧美 日韩 精品 国产| 国产免费av片在线观看野外av| 夜夜爽天天搞| svipshipincom国产片| 9热在线视频观看99| 国产xxxxx性猛交| 欧美激情高清一区二区三区| 久99久视频精品免费| 久久人妻av系列| 国产精品久久久久成人av| 99国产极品粉嫩在线观看| 久久久国产精品麻豆| tocl精华| 我的亚洲天堂| 精品少妇一区二区三区视频日本电影| 91精品三级在线观看| 高清毛片免费观看视频网站 | 精品一区二区三区四区五区乱码| 在线观看免费高清a一片| 乱人伦中国视频| 亚洲精品美女久久久久99蜜臀| 日韩有码中文字幕| 脱女人内裤的视频| 国产亚洲一区二区精品| 两个人免费观看高清视频| 亚洲精华国产精华精| 成人18禁在线播放| 久久久久国产精品人妻aⅴ院 | 三级毛片av免费| 国产99久久九九免费精品| 亚洲欧美激情在线| 国产精品免费一区二区三区在线 | 欧美大码av| 亚洲va日本ⅴa欧美va伊人久久| 老熟女久久久| 黑丝袜美女国产一区| 性少妇av在线| 脱女人内裤的视频| 精品国产亚洲在线| 亚洲片人在线观看| 电影成人av| 精品一区二区三区av网在线观看| 91字幕亚洲| 日本黄色日本黄色录像| 桃红色精品国产亚洲av| 成熟少妇高潮喷水视频| 757午夜福利合集在线观看| 精品一区二区三卡| av网站免费在线观看视频| 国产单亲对白刺激| av不卡在线播放| 亚洲伊人色综图| 看黄色毛片网站| 午夜视频精品福利| 日韩三级视频一区二区三区| 黄色片一级片一级黄色片| 脱女人内裤的视频| 极品少妇高潮喷水抽搐| 久久久久久久久久久久大奶| 欧美中文综合在线视频| 欧美乱色亚洲激情| 精品国产一区二区三区四区第35| 免费人成视频x8x8入口观看| 正在播放国产对白刺激| 不卡一级毛片| 日韩视频一区二区在线观看| 精品午夜福利视频在线观看一区| 一级a爱视频在线免费观看| 一本一本久久a久久精品综合妖精| 一区在线观看完整版| 国产成人免费无遮挡视频| 乱人伦中国视频| 亚洲av第一区精品v没综合| 成人三级做爰电影| 亚洲三区欧美一区| 中亚洲国语对白在线视频| 精品久久久精品久久久| 天天添夜夜摸| 精品国产国语对白av| 久久久水蜜桃国产精品网| av天堂久久9| 日韩精品免费视频一区二区三区| e午夜精品久久久久久久| 女人被躁到高潮嗷嗷叫费观| 两性夫妻黄色片| 亚洲欧美精品综合一区二区三区| 亚洲国产精品一区二区三区在线| 国产精品乱码一区二三区的特点 | 欧美最黄视频在线播放免费 | 黄色视频,在线免费观看| 99久久综合精品五月天人人| 身体一侧抽搐| 国产精品av久久久久免费| 成人三级做爰电影| 亚洲全国av大片| 亚洲一区二区三区欧美精品| 王馨瑶露胸无遮挡在线观看| 在线视频色国产色| 国产av又大| 老熟妇乱子伦视频在线观看| 热99久久久久精品小说推荐| 精品熟女少妇八av免费久了| 欧美国产精品va在线观看不卡| 亚洲,欧美精品.| 国产成人av教育| 岛国在线观看网站| 久久久久久久午夜电影 | 狂野欧美激情性xxxx| 国产99白浆流出| 国产亚洲精品久久久久5区| 久久午夜综合久久蜜桃| av中文乱码字幕在线| 在线观看午夜福利视频| 91九色精品人成在线观看| 亚洲精品美女久久久久99蜜臀| 久久久国产成人精品二区 | 亚洲成人免费av在线播放| 欧美日韩一级在线毛片| 在线观看免费视频日本深夜| 99精国产麻豆久久婷婷| 99精品久久久久人妻精品| 在线观看66精品国产| 91成年电影在线观看| 中文字幕最新亚洲高清| 国产一区有黄有色的免费视频| 国产极品粉嫩免费观看在线| 欧美成人午夜精品| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 人妻 亚洲 视频| 国产成人av教育| 久热爱精品视频在线9| 国产亚洲精品久久久久5区| 精品国内亚洲2022精品成人 | 国产极品粉嫩免费观看在线| 丰满饥渴人妻一区二区三| 人成视频在线观看免费观看| 色综合婷婷激情| 亚洲视频免费观看视频| 免费不卡黄色视频| 妹子高潮喷水视频| 日韩人妻精品一区2区三区| 在线观看舔阴道视频| 成年人免费黄色播放视频| 国产精品 欧美亚洲| 丝瓜视频免费看黄片| 精品福利永久在线观看| netflix在线观看网站| 少妇被粗大的猛进出69影院| 午夜老司机福利片| 一本综合久久免费| 日韩 欧美 亚洲 中文字幕| 91国产中文字幕| 亚洲一区二区三区不卡视频| 人妻 亚洲 视频| 99热网站在线观看| av免费在线观看网站| 国产一区二区三区在线臀色熟女 | 欧美日韩亚洲国产一区二区在线观看 | 国产乱人伦免费视频| 亚洲熟妇熟女久久| av网站在线播放免费| 美女国产高潮福利片在线看| 99re6热这里在线精品视频| 国产区一区二久久| 精品福利观看| 十八禁网站免费在线| 亚洲成国产人片在线观看| 亚洲精品国产区一区二| 欧美日韩乱码在线| 高潮久久久久久久久久久不卡| 亚洲欧美色中文字幕在线| 亚洲五月婷婷丁香| 韩国av一区二区三区四区| 精品一品国产午夜福利视频| tocl精华| 一级a爱视频在线免费观看| 久久久水蜜桃国产精品网| 777米奇影视久久| 好男人电影高清在线观看| 国产有黄有色有爽视频| 在线播放国产精品三级| 免费看a级黄色片| 午夜福利在线观看吧| 午夜日韩欧美国产| avwww免费| 嫁个100分男人电影在线观看| 一本大道久久a久久精品| 国产精品二区激情视频| 国产单亲对白刺激| 精品国产超薄肉色丝袜足j| videosex国产| 国产精品一区二区在线不卡| 在线观看一区二区三区激情| 18在线观看网站| 国产精品乱码一区二三区的特点 | 一区二区三区国产精品乱码| 国产免费现黄频在线看| 国产单亲对白刺激| 悠悠久久av| а√天堂www在线а√下载 | 久久精品熟女亚洲av麻豆精品| 操出白浆在线播放| av免费在线观看网站| 老司机影院毛片| 久久久久久久久免费视频了| 国产淫语在线视频| 亚洲综合色网址| 亚洲第一av免费看| 国产精品秋霞免费鲁丝片| 国产精品98久久久久久宅男小说| 欧美日韩乱码在线| 亚洲欧美一区二区三区久久| 久久精品熟女亚洲av麻豆精品| 99久久精品国产亚洲精品| 亚洲全国av大片| 一级作爱视频免费观看| 丰满饥渴人妻一区二区三| 一级片免费观看大全| 不卡av一区二区三区| 亚洲精品国产精品久久久不卡| 久久久久精品国产欧美久久久| 91成人精品电影| 波多野结衣一区麻豆| 少妇 在线观看| a级毛片在线看网站| 日韩中文字幕欧美一区二区| 视频区图区小说| 色在线成人网| 欧洲精品卡2卡3卡4卡5卡区| 深夜精品福利| 一区二区三区国产精品乱码| 狠狠狠狠99中文字幕| 一二三四在线观看免费中文在| 青草久久国产| 午夜免费成人在线视频| 欧美另类亚洲清纯唯美| 深夜精品福利| 韩国av一区二区三区四区| 最近最新中文字幕大全免费视频| 中文亚洲av片在线观看爽 | 日韩欧美国产一区二区入口| 人妻久久中文字幕网| 1024香蕉在线观看| 在线观看免费高清a一片| 免费观看人在逋| 国产精品一区二区精品视频观看| 一区在线观看完整版| 十八禁高潮呻吟视频| 欧美日韩亚洲国产一区二区在线观看 | 精品福利永久在线观看| 不卡一级毛片| 亚洲免费av在线视频| 午夜福利欧美成人| 精品第一国产精品| 超碰97精品在线观看| 欧美日韩福利视频一区二区| 免费看a级黄色片| 国产激情久久老熟女| 深夜精品福利| 国产亚洲精品一区二区www | 999久久久国产精品视频| 最新在线观看一区二区三区| 一进一出好大好爽视频| 美女福利国产在线| 亚洲av熟女| 激情视频va一区二区三区| 国产精品美女特级片免费视频播放器 | 一本大道久久a久久精品| 国产精品 国内视频| 亚洲精品中文字幕在线视频| 国产亚洲欧美98| 久久精品国产亚洲av香蕉五月 | 亚洲av第一区精品v没综合| 女同久久另类99精品国产91| 无人区码免费观看不卡| 久久狼人影院| 脱女人内裤的视频| 亚洲av第一区精品v没综合| 亚洲熟女毛片儿| 久久久精品国产亚洲av高清涩受| 99re在线观看精品视频| 9色porny在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲精品一卡2卡三卡4卡5卡| 欧美人与性动交α欧美精品济南到| 精品亚洲成a人片在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品九九99| 色婷婷久久久亚洲欧美| tocl精华| 国产在线观看jvid| 淫妇啪啪啪对白视频| 国内毛片毛片毛片毛片毛片| 国产视频一区二区在线看| 国产精品秋霞免费鲁丝片| 国产主播在线观看一区二区| 黄色毛片三级朝国网站| 亚洲av欧美aⅴ国产| 女人被躁到高潮嗷嗷叫费观| 亚洲综合色网址| 亚洲久久久国产精品| 在线观看舔阴道视频| 中国美女看黄片| 丝瓜视频免费看黄片| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文av在线| 老熟妇乱子伦视频在线观看| 亚洲人成电影免费在线| av在线播放免费不卡| 99精品在免费线老司机午夜| 99国产精品一区二区蜜桃av | 亚洲成国产人片在线观看| 久久中文看片网| 搡老熟女国产l中国老女人| 在线观看午夜福利视频| 婷婷精品国产亚洲av在线 | 亚洲av日韩精品久久久久久密| a级毛片黄视频| 中文字幕人妻熟女乱码| av欧美777| 黄色毛片三级朝国网站| 一级片免费观看大全| 成人国产一区最新在线观看| 国内久久婷婷六月综合欲色啪| videosex国产| 久热爱精品视频在线9| 国产精品秋霞免费鲁丝片| 国产免费av片在线观看野外av| 男男h啪啪无遮挡| 色精品久久人妻99蜜桃| av线在线观看网站| 欧美精品av麻豆av| 亚洲九九香蕉| 动漫黄色视频在线观看| 中出人妻视频一区二区| 成年人免费黄色播放视频| a级毛片黄视频| 国产欧美日韩综合在线一区二区| 亚洲片人在线观看| 国产成人精品久久二区二区免费| 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产一区二区入口| 成人影院久久| 免费观看人在逋| 欧美亚洲日本最大视频资源| 久久久国产一区二区| 美女高潮到喷水免费观看| 黄色成人免费大全| 日本vs欧美在线观看视频| 久久精品国产清高在天天线| 亚洲专区字幕在线| 一区在线观看完整版| avwww免费| 丝袜美足系列| 黄片大片在线免费观看| 欧美av亚洲av综合av国产av| 精品国内亚洲2022精品成人 | 久久精品亚洲精品国产色婷小说| 如日韩欧美国产精品一区二区三区| 久久国产精品男人的天堂亚洲| 国精品久久久久久国模美| 亚洲国产欧美一区二区综合| 99国产极品粉嫩在线观看| 精品第一国产精品| 亚洲国产精品一区二区三区在线| ponron亚洲| 成人18禁在线播放| 丝袜美腿诱惑在线| 亚洲精品国产区一区二| 美女扒开内裤让男人捅视频| 亚洲人成77777在线视频| 国产亚洲一区二区精品| 成人三级做爰电影| 91成年电影在线观看| 黄网站色视频无遮挡免费观看| 亚洲第一欧美日韩一区二区三区| 19禁男女啪啪无遮挡网站| 变态另类成人亚洲欧美熟女 | 国产成人一区二区三区免费视频网站| 91麻豆精品激情在线观看国产 | 免费少妇av软件| 真人做人爱边吃奶动态| 十八禁人妻一区二区| 亚洲欧美色中文字幕在线| 免费一级毛片在线播放高清视频 | cao死你这个sao货| 久久精品国产a三级三级三级| 国产精品一区二区在线不卡| 狠狠狠狠99中文字幕| 丝瓜视频免费看黄片| 国产在线精品亚洲第一网站| cao死你这个sao货| av中文乱码字幕在线| 怎么达到女性高潮| 国产精品久久视频播放| 免费黄频网站在线观看国产| 国产激情欧美一区二区| 18在线观看网站| 99re在线观看精品视频| 亚洲午夜理论影院| 国产欧美日韩一区二区精品| 少妇被粗大的猛进出69影院| 老司机深夜福利视频在线观看| 国产极品粉嫩免费观看在线| 日韩 欧美 亚洲 中文字幕| 亚洲情色 制服丝袜| 90打野战视频偷拍视频| 日本撒尿小便嘘嘘汇集6|