• <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
    嫩草影院入口| 精品免费久久久久久久清纯| 91在线观看av| 美女扒开内裤让男人捅视频| 久久久久九九精品影院| 国产亚洲欧美在线一区二区| 色在线成人网| 国语自产精品视频在线第100页| 亚洲av成人精品一区久久| x7x7x7水蜜桃| 男女之事视频高清在线观看| 岛国视频午夜一区免费看| 手机成人av网站| 麻豆一二三区av精品| 成年免费大片在线观看| a级毛片在线看网站| 亚洲18禁久久av| 老司机在亚洲福利影院| 啦啦啦韩国在线观看视频| 熟女少妇亚洲综合色aaa.| 好男人在线观看高清免费视频| 亚洲aⅴ乱码一区二区在线播放| 天堂√8在线中文| 久久精品国产亚洲av香蕉五月| 97人妻精品一区二区三区麻豆| ponron亚洲| 国内精品美女久久久久久| 国产亚洲av高清不卡| 欧美大码av| 欧美zozozo另类| 88av欧美| 亚洲国产看品久久| 无遮挡黄片免费观看| 婷婷亚洲欧美| 欧美乱码精品一区二区三区| 亚洲男人的天堂狠狠| 后天国语完整版免费观看| 91久久精品国产一区二区成人 | 夜夜爽天天搞| 欧美日韩中文字幕国产精品一区二区三区| 久久天堂一区二区三区四区| 亚洲 国产 在线| 变态另类成人亚洲欧美熟女| 国产精品99久久久久久久久| 亚洲天堂国产精品一区在线| 在线免费观看不下载黄p国产 | 淫秽高清视频在线观看| 久久天堂一区二区三区四区| 婷婷六月久久综合丁香| 亚洲午夜精品一区,二区,三区| 国产伦一二天堂av在线观看| 99精品久久久久人妻精品| 18美女黄网站色大片免费观看| 国内毛片毛片毛片毛片毛片| 97超视频在线观看视频| av国产免费在线观看| 欧美+亚洲+日韩+国产| 小说图片视频综合网站| 国内精品一区二区在线观看| 久久久久免费精品人妻一区二区| 午夜精品一区二区三区免费看| 亚洲av电影不卡..在线观看| 欧美日韩福利视频一区二区| 一本一本综合久久| 欧美黑人巨大hd| 日日干狠狠操夜夜爽| 色综合站精品国产| 男女床上黄色一级片免费看| 麻豆成人av在线观看| 欧美绝顶高潮抽搐喷水| 美女被艹到高潮喷水动态| 级片在线观看| 精品久久蜜臀av无| 亚洲欧美激情综合另类| 男插女下体视频免费在线播放| 亚洲精品美女久久久久99蜜臀| 午夜a级毛片| 亚洲自拍偷在线| 久久久久久九九精品二区国产| 色哟哟哟哟哟哟| 亚洲国产高清在线一区二区三| 国产精品一区二区精品视频观看| 波多野结衣巨乳人妻| 成人欧美大片| 熟女少妇亚洲综合色aaa.| 午夜福利成人在线免费观看| 亚洲成人久久爱视频| 熟女人妻精品中文字幕| 真人一进一出gif抽搐免费| 亚洲五月婷婷丁香| 亚洲欧美日韩卡通动漫| 琪琪午夜伦伦电影理论片6080| 波多野结衣高清作品| 精品国产美女av久久久久小说| 999久久久国产精品视频| 啪啪无遮挡十八禁网站| 国产免费男女视频| 国内精品一区二区在线观看| 日韩精品中文字幕看吧| 亚洲av五月六月丁香网| 国产亚洲精品久久久久久毛片| 国产亚洲av高清不卡| 亚洲中文字幕一区二区三区有码在线看 | 色在线成人网| 日日夜夜操网爽| 亚洲精华国产精华精| 亚洲熟妇熟女久久| 五月玫瑰六月丁香| 五月伊人婷婷丁香| 亚洲片人在线观看| 亚洲av日韩精品久久久久久密| 女人高潮潮喷娇喘18禁视频| 岛国在线观看网站| 精品一区二区三区av网在线观看| 精品久久久久久久久久免费视频| 一级毛片女人18水好多| 啦啦啦观看免费观看视频高清| 级片在线观看| 日韩大尺度精品在线看网址| 国产激情偷乱视频一区二区| 一区二区三区国产精品乱码| 久久久久免费精品人妻一区二区| 精品福利观看| 国产高潮美女av| 国产精品久久久av美女十八| 亚洲美女黄片视频| 亚洲人成网站在线播放欧美日韩| 国产午夜福利久久久久久| 悠悠久久av| 最好的美女福利视频网| 波多野结衣高清无吗| 1000部很黄的大片| 欧美一区二区精品小视频在线| 1000部很黄的大片| 日韩精品青青久久久久久| 成人av一区二区三区在线看| 国产男靠女视频免费网站| 最新中文字幕久久久久 | 精品国产三级普通话版| 日本黄色视频三级网站网址| 久久久久国产精品人妻aⅴ院| 久久精品亚洲精品国产色婷小说| 少妇熟女aⅴ在线视频| 999久久久国产精品视频| 国产精品一及| 日本a在线网址| 成年人黄色毛片网站| 色综合欧美亚洲国产小说| 真实男女啪啪啪动态图| 中文字幕久久专区| 天堂网av新在线| 一级毛片精品| 国产伦一二天堂av在线观看| h日本视频在线播放| 两个人的视频大全免费| 欧美又色又爽又黄视频| 一区福利在线观看| 日韩欧美三级三区| 成年女人永久免费观看视频| 制服丝袜大香蕉在线| 欧美激情在线99| 女警被强在线播放| 欧美另类亚洲清纯唯美| 脱女人内裤的视频| 国产1区2区3区精品| 亚洲av片天天在线观看| 一二三四在线观看免费中文在| 黑人巨大精品欧美一区二区mp4| 999久久久国产精品视频| 免费av毛片视频| 国产精品av久久久久免费| 在线观看美女被高潮喷水网站 | 欧美黄色淫秽网站| 色尼玛亚洲综合影院| 成人一区二区视频在线观看| 欧美日韩国产亚洲二区| 成年女人看的毛片在线观看| 国产精品久久久久久亚洲av鲁大| av天堂在线播放| 黄色女人牲交| 九九热线精品视视频播放| 亚洲黑人精品在线| 亚洲精华国产精华精| 久久午夜亚洲精品久久| 欧美中文日本在线观看视频| 国产高清videossex| 真人做人爱边吃奶动态| 免费在线观看视频国产中文字幕亚洲| h日本视频在线播放| 国产三级中文精品| 婷婷六月久久综合丁香| 亚洲无线在线观看| 给我免费播放毛片高清在线观看| 老司机午夜十八禁免费视频| 99久国产av精品| 91老司机精品| 亚洲一区二区三区不卡视频| 国产视频一区二区在线看| 老熟妇仑乱视频hdxx| 美女午夜性视频免费| 亚洲乱码一区二区免费版| 免费看a级黄色片| 欧美午夜高清在线| 听说在线观看完整版免费高清| 国产成人精品久久二区二区91| 后天国语完整版免费观看| 欧美3d第一页| 国产又黄又爽又无遮挡在线| 91九色精品人成在线观看| 国产精品综合久久久久久久免费| 亚洲乱码一区二区免费版| 国产精品一及| 欧美成人性av电影在线观看| 国产真实乱freesex| 成人午夜高清在线视频| 色播亚洲综合网| 麻豆国产av国片精品| 校园春色视频在线观看| 国产精品电影一区二区三区| 国产精品99久久99久久久不卡| 午夜福利高清视频| 婷婷六月久久综合丁香| 九九热线精品视视频播放| 一本久久中文字幕| 中出人妻视频一区二区| 国产麻豆成人av免费视频| 日韩 欧美 亚洲 中文字幕| 日本五十路高清| 成人特级av手机在线观看| 国产一区二区激情短视频| 国产伦精品一区二区三区视频9 | 国产视频一区二区在线看| 欧美最黄视频在线播放免费| 在线看三级毛片| 丰满人妻熟妇乱又伦精品不卡| 成人一区二区视频在线观看| 在线a可以看的网站| 久久久久久人人人人人| av欧美777| 亚洲,欧美精品.| 成人无遮挡网站| 亚洲人成伊人成综合网2020| 亚洲成人中文字幕在线播放| av欧美777| 亚洲av电影在线进入| 国产成人影院久久av| 老司机午夜十八禁免费视频| 国产aⅴ精品一区二区三区波| 天堂av国产一区二区熟女人妻| 国产亚洲精品久久久久久毛片| 人人妻人人澡欧美一区二区| 国产精品,欧美在线| 欧美成人性av电影在线观看| 国内精品久久久久久久电影| 级片在线观看| 一个人看视频在线观看www免费 | 日韩中文字幕欧美一区二区| 免费看a级黄色片| 久久这里只有精品中国| 久久久精品欧美日韩精品| 国产麻豆成人av免费视频| 一本精品99久久精品77| 久久精品亚洲精品国产色婷小说| 一进一出抽搐gif免费好疼| 婷婷六月久久综合丁香| 香蕉丝袜av| 成人av在线播放网站| 真实男女啪啪啪动态图| 免费av不卡在线播放| 1024香蕉在线观看| 91在线精品国自产拍蜜月 | 欧美xxxx黑人xx丫x性爽| 欧美中文日本在线观看视频| 天天一区二区日本电影三级| 久久99热这里只有精品18| 日本免费a在线| 校园春色视频在线观看| 在线观看免费视频日本深夜| 一级a爱片免费观看的视频| 国产精品av久久久久免费| 欧美黄色淫秽网站| 真人做人爱边吃奶动态| 97碰自拍视频| 国产男靠女视频免费网站| 欧美日韩国产亚洲二区| 久久久久久人人人人人| 午夜免费观看网址| 夜夜爽天天搞| 亚洲精品乱码久久久v下载方式 | 黑人操中国人逼视频| netflix在线观看网站| 久久精品91蜜桃| 免费av毛片视频| 欧美在线一区亚洲| 精品午夜福利视频在线观看一区| 美女被艹到高潮喷水动态| 一个人看视频在线观看www免费 | 成人三级黄色视频| 91在线精品国自产拍蜜月 | 久久久国产成人免费| 亚洲人成伊人成综合网2020| 最近最新中文字幕大全电影3| 亚洲人成网站高清观看| 国产aⅴ精品一区二区三区波| 国产视频一区二区在线看| 五月玫瑰六月丁香| 久久久色成人| 99热精品在线国产| 久久久久久久久免费视频了| 日本黄色视频三级网站网址| 最好的美女福利视频网| 最近视频中文字幕2019在线8| 两性夫妻黄色片| 可以在线观看的亚洲视频| 小蜜桃在线观看免费完整版高清| 99国产极品粉嫩在线观看| 亚洲国产精品sss在线观看| 美女高潮喷水抽搐中文字幕| 狠狠狠狠99中文字幕| 人妻久久中文字幕网| 免费看光身美女| 午夜福利欧美成人| 欧美黑人巨大hd| www.精华液| 欧美绝顶高潮抽搐喷水| 男人舔奶头视频| 欧美日韩中文字幕国产精品一区二区三区| 好男人在线观看高清免费视频| 亚洲 欧美一区二区三区| 黑人巨大精品欧美一区二区mp4| 热99在线观看视频| 日日干狠狠操夜夜爽| 老司机午夜十八禁免费视频| 看黄色毛片网站| 久久久久亚洲av毛片大全| 最近最新中文字幕大全免费视频| 精品无人区乱码1区二区| 特级一级黄色大片| 无遮挡黄片免费观看| 俄罗斯特黄特色一大片| 日本黄色片子视频| 91九色精品人成在线观看| 91av网一区二区| 欧美成人一区二区免费高清观看 | 国产三级中文精品| 综合色av麻豆| 非洲黑人性xxxx精品又粗又长| 国产激情偷乱视频一区二区| 黄色 视频免费看| 国产亚洲欧美在线一区二区| 国产69精品久久久久777片 | 听说在线观看完整版免费高清| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 一a级毛片在线观看| a级毛片在线看网站| 最近最新中文字幕大全电影3| 日本黄大片高清| 欧美日韩中文字幕国产精品一区二区三区| 最新在线观看一区二区三区| 又黄又爽又免费观看的视频| 久久天堂一区二区三区四区| 最近在线观看免费完整版| 9191精品国产免费久久| 亚洲精品在线观看二区| 色吧在线观看| 国产精华一区二区三区| 天堂影院成人在线观看| 国产高清视频在线播放一区| 叶爱在线成人免费视频播放| 成人鲁丝片一二三区免费| 亚洲国产精品合色在线| 少妇人妻一区二区三区视频| 精品免费久久久久久久清纯| 美女扒开内裤让男人捅视频| 国产精品久久久久久久电影 | 岛国视频午夜一区免费看| 真人一进一出gif抽搐免费| 国产熟女xx| 天天躁狠狠躁夜夜躁狠狠躁| 制服丝袜大香蕉在线| 欧美中文综合在线视频| 久久精品国产综合久久久| 成人三级黄色视频| www.www免费av| 两个人看的免费小视频| 男人舔奶头视频| 99久久无色码亚洲精品果冻| 亚洲人成网站高清观看| 日本免费一区二区三区高清不卡| 少妇裸体淫交视频免费看高清| 老司机午夜福利在线观看视频| 欧美日韩综合久久久久久 | 一级黄色大片毛片| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 国产亚洲精品av在线| 精品国产三级普通话版| 男人和女人高潮做爰伦理| 国产1区2区3区精品| 日韩欧美精品v在线| 国产成人精品无人区| 国内久久婷婷六月综合欲色啪| 国内精品美女久久久久久| 国产午夜精品久久久久久| 女人被狂操c到高潮| 女同久久另类99精品国产91| 黄色视频,在线免费观看| 国产v大片淫在线免费观看| 三级男女做爰猛烈吃奶摸视频| 免费高清视频大片| 女警被强在线播放| 搞女人的毛片| 久久天堂一区二区三区四区| 国产成人av教育| 亚洲欧美精品综合一区二区三区| 变态另类丝袜制服| 国产伦人伦偷精品视频| 亚洲精品美女久久久久99蜜臀| 国产伦精品一区二区三区视频9 | 日本成人三级电影网站| 天天添夜夜摸| 国内精品久久久久精免费| 一进一出好大好爽视频| 成人亚洲精品av一区二区| 日本黄大片高清| 日本一二三区视频观看| 91九色精品人成在线观看| 91麻豆精品激情在线观看国产| 欧美3d第一页| 午夜免费成人在线视频| 熟女人妻精品中文字幕| 久久香蕉精品热| 观看美女的网站| 国产亚洲精品综合一区在线观看| 亚洲熟妇熟女久久| 怎么达到女性高潮| 国产高清激情床上av| 女人被狂操c到高潮| 国产精品影院久久| 午夜精品久久久久久毛片777| 亚洲专区国产一区二区| 亚洲成av人片免费观看| a在线观看视频网站| 一级毛片女人18水好多| 亚洲欧美日韩无卡精品| 可以在线观看的亚洲视频| 熟女电影av网| 免费观看人在逋| 国产亚洲精品av在线| 99在线人妻在线中文字幕| 人妻丰满熟妇av一区二区三区| 久久天堂一区二区三区四区| 国产精品久久电影中文字幕| 少妇人妻一区二区三区视频| 亚洲成a人片在线一区二区| 久久久久亚洲av毛片大全| 午夜精品久久久久久毛片777| 久久精品91蜜桃| 亚洲精品久久国产高清桃花| 欧美日韩瑟瑟在线播放| 舔av片在线| 小蜜桃在线观看免费完整版高清| 国内毛片毛片毛片毛片毛片| 成人鲁丝片一二三区免费| 大型黄色视频在线免费观看| 久久亚洲真实| 国产午夜精品论理片| 国产午夜精品久久久久久| 国产精品女同一区二区软件 | 精品一区二区三区av网在线观看| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 免费在线观看日本一区| 精品国产三级普通话版| 中文字幕精品亚洲无线码一区| 免费av不卡在线播放| 国内毛片毛片毛片毛片毛片| 国产探花在线观看一区二区| 久久精品aⅴ一区二区三区四区| 18禁黄网站禁片午夜丰满| 亚洲熟女毛片儿| 亚洲aⅴ乱码一区二区在线播放| 久久久久久久久中文| 国产精品免费一区二区三区在线| 夜夜夜夜夜久久久久| 12—13女人毛片做爰片一| 亚洲精品久久国产高清桃花| 又黄又粗又硬又大视频| 欧美黑人欧美精品刺激| 法律面前人人平等表现在哪些方面| 日韩免费av在线播放| 这个男人来自地球电影免费观看| 老熟妇乱子伦视频在线观看| 色综合站精品国产| 精品久久久久久久久久免费视频| 一区福利在线观看| 99久国产av精品| 欧美+亚洲+日韩+国产| 亚洲国产色片| 1000部很黄的大片| 精品欧美国产一区二区三| 国产乱人伦免费视频| 一个人观看的视频www高清免费观看 | av在线蜜桃| 一级作爱视频免费观看| 亚洲自偷自拍图片 自拍| 亚洲欧美一区二区三区黑人| 天天一区二区日本电影三级| 久久久久久久精品吃奶| 国产精品一区二区精品视频观看| 成人无遮挡网站| 欧美日韩综合久久久久久 | 在线国产一区二区在线| 无限看片的www在线观看| 亚洲精品中文字幕一二三四区| 成人高潮视频无遮挡免费网站| 亚洲国产欧美一区二区综合| 亚洲成人免费电影在线观看| 成人无遮挡网站| 色哟哟哟哟哟哟| 狂野欧美白嫩少妇大欣赏| 1000部很黄的大片| 亚洲人与动物交配视频| 女人高潮潮喷娇喘18禁视频| www国产在线视频色| 男人舔女人下体高潮全视频| 高清在线国产一区| 欧美黄色片欧美黄色片| 欧美国产日韩亚洲一区| 精华霜和精华液先用哪个| 欧美3d第一页| 97人妻精品一区二区三区麻豆| 极品教师在线免费播放| 亚洲欧美精品综合久久99| 国产成人系列免费观看| 日韩欧美在线二视频| 日韩欧美三级三区| 国产精品野战在线观看| 免费看美女性在线毛片视频| 免费无遮挡裸体视频| 99riav亚洲国产免费| 夜夜看夜夜爽夜夜摸| 精品熟女少妇八av免费久了| 午夜福利18| 日本在线视频免费播放| 国产精品免费一区二区三区在线| 99热这里只有是精品50| 亚洲av日韩精品久久久久久密| 小蜜桃在线观看免费完整版高清| 亚洲狠狠婷婷综合久久图片| 老司机深夜福利视频在线观看| 国内毛片毛片毛片毛片毛片| 欧美绝顶高潮抽搐喷水| av视频在线观看入口| 黑人巨大精品欧美一区二区mp4| 欧美色视频一区免费| 男女午夜视频在线观看| 精品午夜福利视频在线观看一区| 草草在线视频免费看| 精品国产乱子伦一区二区三区| 久久久久久大精品| 亚洲欧美日韩东京热| 两个人视频免费观看高清| 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 99热这里只有精品一区 | 一二三四在线观看免费中文在| 国产精品女同一区二区软件 | 免费在线观看成人毛片| 1000部很黄的大片| 一本一本综合久久| 精品国产亚洲在线| 麻豆av在线久日| 黄色成人免费大全| 亚洲精华国产精华精| 国产毛片a区久久久久| 免费大片18禁| 韩国av一区二区三区四区| 亚洲国产欧美人成| 亚洲在线自拍视频| 免费av毛片视频| 国产97色在线日韩免费| 午夜福利高清视频| 午夜免费激情av| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区三区在线臀色熟女| 一级毛片精品| 亚洲av美国av| 国产综合懂色| 国产精品久久久人人做人人爽| 成人三级做爰电影| 男女午夜视频在线观看| 欧美色视频一区免费| 国产99白浆流出| 亚洲国产欧美一区二区综合| 久久精品国产综合久久久| 免费高清视频大片| 欧美日韩乱码在线| 欧美黑人巨大hd| 国产高清三级在线| 国产精品香港三级国产av潘金莲| 又爽又黄无遮挡网站| 美女午夜性视频免费| 亚洲在线自拍视频| 免费大片18禁| 夜夜看夜夜爽夜夜摸| 啪啪无遮挡十八禁网站| 99久久久亚洲精品蜜臀av| 又爽又黄无遮挡网站| 88av欧美| 久久久久免费精品人妻一区二区| 母亲3免费完整高清在线观看| 亚洲一区二区三区色噜噜| 99热6这里只有精品|