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

    復(fù)合材料薄板結(jié)構(gòu)中的聲學(xué)黑洞效應(yīng)探究

    2023-01-31 13:59:58鄭鋒黃薇季宏麗裘進浩
    航空學(xué)報 2023年1期
    關(guān)鍵詞:波場鋪層波速

    鄭鋒,黃薇,季宏麗,*,裘進浩

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

    2.南京理工大學(xué) 機械工程學(xué)院,南京 210094

    聲學(xué)黑洞(Acoustic Black Hole,ABH)作為一種新型波操控和被動阻尼技術(shù),近年來得到了科研工作者的廣泛關(guān)注和研究,在實際工程中也取得了很好的運用。ABH一般通過調(diào)整結(jié)構(gòu)相關(guān)參數(shù)以改變結(jié)構(gòu)阻抗,使得結(jié)構(gòu)中彎曲波的相速度和群速度降低,從而實現(xiàn)波在某一區(qū)域高度聚集。在聚集區(qū)域粘貼少量阻尼材料就能實現(xiàn)高效的能量耗散,最終實現(xiàn)減振降噪的目的。

    ABH常見的實現(xiàn)方式是裁剪結(jié)構(gòu)厚度,基于此產(chǎn)生了多種設(shè)計形式,如一維的有錐形楔[1]、螺旋式[2]、雙葉式[3]等,二維的有圓形、環(huán)形和矩形內(nèi) 嵌 式[4]以 及 周 期 帶 狀A(yù)BH隧 穿 結(jié) 構(gòu)[5-6]。此外,ABH還可以通過改變材料特性來實現(xiàn),例如在黏彈性層上周期性嵌入ABHs以降低振動與噪聲[7],利用梁末端的厚度與模量參數(shù)梯度變化以降低反射[8]。聲學(xué)黑洞從波動角度實現(xiàn)了對結(jié)構(gòu)中彈性波傳播、能量傳遞和消耗的控制,廣泛應(yīng)用于結(jié)構(gòu)的減振降噪[9-10]、能量回收[11-12]和波操縱[13-18]等方面。對于波操縱,研究人員采用有限元法仿真[13-15]、波場ABH實驗測試[14-17]、波軌跡分析[18]等多種方法,對ABH結(jié)構(gòu)中的波傳播和聚集過程進行了深入研究。

    現(xiàn)有關(guān)于聲學(xué)黑洞的研究主要是針對均勻材料結(jié)構(gòu),例如金屬和樹脂材料,而對復(fù)合材料結(jié)構(gòu)中的ABH效應(yīng)的研究還很少,其主要原因是復(fù)合材料結(jié)構(gòu)的復(fù)雜性和各向異性給建模和仿真分析帶來了巨大的挑戰(zhàn)。碳纖維復(fù)合材料具有強度高、模量大、耐腐蝕性能好等優(yōu)點,廣泛應(yīng)用于航空、汽車和船舶等領(lǐng)域,因此研究碳纖維復(fù)合材料結(jié)構(gòu)中的ABH效應(yīng)對擴展ABH的應(yīng)用有著重要意義。碳纖維復(fù)合材料結(jié)構(gòu)中的波能量聚焦與耗散效應(yīng)主要有2個因素:不同鋪層角度碳纖維層引起波速的各向異性[19-21]和較高的損耗因子[22]。

    因此,本文選取內(nèi)嵌式的碳纖維復(fù)合材料ABH薄板結(jié)構(gòu)(CFRP-ABH),對其波傳播特性進行研究分析。通過有限元仿真計算,探究CFRP-ABH結(jié)構(gòu)中的能量聚集特性,驗證ABH效應(yīng)的存在性;同時通過對波傳播過程的量化分析,探究不同鋪層角度對能量聚集特性的影響。此外,在ABH中心區(qū)域粘貼少量阻尼材料,計算系統(tǒng)的動力學(xué)響應(yīng),并通過實驗測試來驗證CFRP-ABH的寬頻減振特性。

    1 CFRP-ABH結(jié)構(gòu)的設(shè)計

    飛機壁板是飛機的組件單元之一,一般由蒙皮、長桁、角片、框等零件組成,主要存在于飛機的機身、機翼等結(jié)構(gòu)中。針對飛機壁板結(jié)構(gòu),本文設(shè)計出內(nèi)嵌入ABH的碳纖維復(fù)合材料薄板結(jié)構(gòu),板的整體尺寸為680 mm×480 mm×5 mm,其中內(nèi)嵌的ABH結(jié)構(gòu)位于薄板幾何中心,結(jié)構(gòu)的整體及剖面示意圖如圖1所示,圖中:h和r分別表示厚度與半徑;h0和r1為中心平臺厚度與半徑;h1為 總 厚 度;r2為ABH區(qū) 域 半 徑;m為 冪 指數(shù);ε為系數(shù)。由于ABH截斷厚度的存在,CFRP-ABH結(jié)構(gòu)由3部分組成:中心平臺(均勻厚度)、ABH區(qū)域(變厚度)、其他區(qū)域(均勻厚度),其厚度滿足式(1):

    圖1 CFRP-ABH結(jié)構(gòu)示意圖Fig. 1 Schematic diagram of CFRP-ABH structure

    內(nèi)嵌式ABH設(shè)計利用結(jié)構(gòu)剛度的局部破壞實現(xiàn)振動能量在局部區(qū)域聚集,因此一定程度上會降低原有結(jié)構(gòu)的強度和可靠性。對于本文設(shè)計的CFRP-ABH,變厚度區(qū)域的纖維層會被切斷以滿足厚度變化規(guī)律。纖維層的破壞會直接影響結(jié)構(gòu)的強度,但在ABH尺寸設(shè)計合理的情況下,原結(jié)構(gòu)整體上仍保留一定的力學(xué)性能,可用于航空結(jié)構(gòu)的非主要承力構(gòu)件。

    2 CFRP-ABH結(jié)構(gòu)的波能量聚集特性

    2.1 有限元模型

    在ABAQUS有限元軟件中建立模型,探究彎曲波在CFRP-ABH結(jié)構(gòu)中的傳播和聚集特性。為了避免自由邊界和固支邊界處反射波對分析結(jié)果的影響,需要對圖1的設(shè)計模型進行適當(dāng)調(diào)整,最終得到用于時域波場分析的模型,如圖2所示。延長區(qū)域的主要作用是為了在分析的時間范圍內(nèi),波在ABH區(qū)域的傳播不受固支邊界的影響。板的左端為固支邊界,在板的右側(cè)面施加一個中心頻率為20 kHz的五波峰激勵,方向垂直于板平面,以此產(chǎn)生平行入射波,讓彎曲波主要沿長度方向傳播并經(jīng)過ABH區(qū)域。為了保證計算精度,每個波長至少有10個網(wǎng)格單元,單元類型為C3D20R。

    圖2 時域分析的結(jié)構(gòu)示意圖Fig. 2 Structure diagram of time domain analysis

    為了研究鋪層角度對波能量聚集特性的影響,設(shè)置5種常用鋪層角度的CFRP-ABH板,并與金屬均勻材料ABH鋁板(Al-ABH)的結(jié)果進行對比。CFRP-ABH板的鋪層共有27層,其中上下2個表面是正交編織布,簡化為0°層;剩余鋪層為纖維層,各個模型的纖維層角度各不相同,分別為0°單向?qū)印?0°單向?qū)?、?ˉ/90]6S、[4ˉ5ˉ/?45]6S和[0ˉ/?45/45/90]3S。為 方 便 表 示,將5個CFRP-ABH模型分別記為0/0、90/90、0/90、45/?45、0/?45/45/90,將Al-ABH記為Al。

    CFRP材料選擇T700/Epoxy,阻尼材料選擇丁基橡膠,相關(guān)材料參數(shù)(包括金屬Al)如表1所示。已有研究表明,在ABH區(qū)域粘貼一般常見的阻尼材料對波能量聚集現(xiàn)象的影響較?。?3]。因此,時域分析時板均不需要粘貼阻尼材料,而在減振特性的仿真和實驗分析中,需要在ABH區(qū)域粘貼少量阻尼材料以實現(xiàn)振動能量的高效耗散。

    2.2 不同鋪層角度下的波能量聚集特性

    提取6個模型在不同時刻下垂直于板面的位移,作為結(jié)構(gòu)的波場源數(shù)據(jù)。觀察在整個時間段的波場,在約0.5 ms左右時彎曲波發(fā)生了完全聚集。以模型Al和45/?45為例,其波場如圖3所示,CFRP-ABH板具有和傳統(tǒng)均勻Al-ABH板中類似的波能量聚集現(xiàn)象,表明了CFRP-ABH結(jié)構(gòu)也存在ABH效應(yīng)。6個模型中波最快在0.63 ms時到達左端固支邊界,因此分析選取的時間段為0~0.6 ms。

    圖3 模型Al和45/?45完全聚集時的波場Fig. 3 Wave fields of model Al and 45/?45 at complete aggregation

    為了進一步分析鋪層角度對波能量聚集特性的影響,需要對波能量聚集現(xiàn)象進行量化分析。根據(jù)T時刻ABH區(qū)域中的波幅值能量分布情況,定義2個指標(biāo):平均分布位置xˉ、聚集程度(Focusing level, Fl),分別為式(2)和式(3):

    式中:uABH*為ABH區(qū)域中超過設(shè)定閾值的節(jié)點位移(U3方向);SABH*為對應(yīng)區(qū)域的總面積;uABH為ABH區(qū)域中的節(jié)點位移;SABH為對應(yīng)區(qū)域的總面積;∑表示對滿足條件的所有節(jié)點求和;xi為節(jié)點的橫坐標(biāo)值。

    設(shè)定一個閾值,形式為某時刻ABH區(qū)域中節(jié)點位移幅值最大值的百分比。在t時刻下,統(tǒng)計ABH區(qū)域中超過閾值的波能量,作為ABH區(qū)域中的主要能量。在相同的閾值標(biāo)準(zhǔn)下,如果ABH區(qū)域的主要能量分布越集中,則波的聚集程度越大,F(xiàn)l值也越大。本文中選取30%作為閾值。例如,對于Al-ABH板,波在0.49 ms時發(fā)生了很好的聚集,如圖4所示。在30%的閾值標(biāo)準(zhǔn)下,需要統(tǒng)計的幅值能量均位于圖中紅色虛線劃定的區(qū)域內(nèi),說明設(shè)定的閾值能反映ABH區(qū)域中的波幅值能量分布情況,且在同一閾值標(biāo)準(zhǔn)下能反映出不同鋪層之間的差異。

    圖4 模型Al在0.49 ms時的波場Fig. 4 Wave field of model Al at 0.49 ms

    6個模型中彎曲波的平均分布位置xˉ和聚集程度Fl隨時間變化情況分別如圖5和圖6所示。圖5中虛線xˉ=0對應(yīng)ABH中心位置,ABH區(qū)域?qū)?yīng)的xˉ范圍為?0.14~0.14 m。彎曲波進入ABH后波長被壓縮,幅值增大,在波傳播經(jīng)過ABH區(qū)域的過程中逐漸產(chǎn)生能量聚集現(xiàn)象。xˉ和Fl值首次出現(xiàn)的時刻,反映了波在不同鋪層角度中的沿入射方向的波速差異。ABH區(qū)域中,在波不斷向前傳播并最終完全聚集的過程中,波幅值能量的平均分布位置前移,xˉ值一直變小。而在這個過程中,波的聚集程度先降低后升高。波向前傳播,而在這一過程中波長會被壓縮,波速降低,使得波的聚集程度會在一定時間內(nèi)保持較低的水平。當(dāng)波開始聚集時,其聚集程度會迅速增大直至最大,此時波發(fā)生完全聚集,xˉ值達到傳播過程中的最小值,并在一定的時間段內(nèi)保持相對穩(wěn)定。如圖5中黑色虛線框標(biāo)出,除90/90鋪層外的其他鋪層,在0.47~0.53 ms時間段內(nèi)xˉ值的變化范圍很小,說明在該時間段內(nèi)波完成了很好的聚集。定義的Fl對波的幅值能量分布較為敏感,因此在0.47~0.53 ms時間段內(nèi)Fl值變化并不平穩(wěn),而是存在一個最大值,如圖6所示,表明此時的波幅值能量聚集程度達到最大。圖7為各個模型的Fl值達到最大時的波場。

    圖5 各個模型的平均分布位置xˉFig. 5 Average location of each modelxˉ

    圖6 各個模型的聚集程度FlFig. 6 Aggregation level of each model F1

    圖7 各個模型的Fl值最大時的波場Fig. 7 Wave fields of each model with maximum Fl value

    綜合考慮xˉ和Fl,選擇Fl最大值對應(yīng)的xˉ值作為波在x方向上的聚集位置,由此得到各個模型中波的聚集時間、位置和Fl峰值,如表2所示。其中,波在0/0鋪層中波速最大,波速降低最難,波幅值能量的聚集位置離ABH中心最遠(yuǎn),但聚集程度最低;在90/90鋪層中波速最小,在經(jīng)過ABH中心之前就完全聚集,聚集位置為正值,離ABH中心最近,但聚集程度最高;其他鋪層(包括Al)中的聚集位置和聚集程度則位于0/0鋪層和90/90鋪層之間??梢园l(fā)現(xiàn),沿入射方向上的波速越小,聚集位置離入射位置更近,而波聚集程度受到聚集位置的影響,即波聚集位置處的厚度越小,波的聚集程度越高。因此,可以得出結(jié)論:在CFRP-ABH薄板結(jié)構(gòu)中,鋪層角度主要通過改變結(jié)構(gòu)中的波速來影響波在ABH區(qū)域中的聚集位置,從而影響波最終的聚集程度。0/0鋪層和90/90鋪層分別表示結(jié)構(gòu)中波速最快和最慢時的鋪層角度。

    2.3 橢圓形ABH設(shè)計

    已有學(xué)者研究了其他形狀的ABH,例如將其設(shè)計為橢圓形透鏡實現(xiàn)對非ABH區(qū)域的波操控[24]。本文主要關(guān)注ABH區(qū)域的波動現(xiàn)象,結(jié)合各向異性對波能量聚集特性的影響規(guī)律,對于CFRP-ABH結(jié)構(gòu),可以利用鋪層角度和ABH形狀設(shè)計進一步實現(xiàn)對波傳播的操控。在特殊鋪層角度設(shè)計下,波在未傳播至ABH中心之前就能夠發(fā)生聚集,這一現(xiàn)象是非常特殊的,而將ABH設(shè)計為橢圓形對這一波操控現(xiàn)象有促進作用。以90/90鋪層為例,將ABH形狀設(shè)計為橢圓形,如圖8所示,ABH區(qū)域的厚度變化仍然滿足ABH原理,即式(4)~式(6):

    圖8 橢圓CFRP-ABH的ABH截面Fig. 8 ABH profile of elliptical CFRP-ABH

    式中:長軸半徑A=280 mm,短軸半徑B=100 mm;中心平臺也是橢圓,其長短軸半徑a與b分別為20 mm和14 mm;θ表示徑向方向與x軸的夾角。

    計算橢圓CFRP-ABH結(jié)構(gòu)(Ell-90/90)中波完全聚集時的聚集位置、時間和聚集程度如表2所示,聚集時刻的波場如圖9所示。可以看出,橢圓形設(shè)計會使得波進一步提前聚集,聚集位置離入射邊界更近,但聚集程度也相應(yīng)地降低。因此,復(fù)合材料結(jié)構(gòu)的各向異性特征有利于實現(xiàn)ABH效應(yīng)的波操控目的。

    表2 各個模型中波的聚集時間、位置和Fl峰值Table 2 Aggregation time, location and Fl peak of wave in each model

    圖9 橢圓CFRP-ABH完全聚集時的波場Fig. 9 Wave field of elliptical CFRP-ABH with com?plete aggregation

    3 CFRP-ABH結(jié)構(gòu)的寬頻減振特性

    通過第2節(jié)的分析,CFRP-ABH結(jié)構(gòu)的能量聚集特性得到了驗證。而對于CFRP-ABH結(jié)構(gòu)的頻域特性,需要在ABAQUS中建立模型,計算結(jié)構(gòu)的動力學(xué)響應(yīng),分析結(jié)構(gòu)的寬頻減振性能。為了方便有限元建模計算和實驗測試,模型選用圖1所示的模型,即尺寸為680 mm×480 mm×5 mm,鋪層角度采用0/90正交。

    為了實現(xiàn)能量耗散,需要在ABH中心位置粘貼少量阻尼材料(位于ABH凹痕一側(cè)),直徑為129 mm,厚度為2.2 mm。設(shè)置相同鋪層結(jié)構(gòu)、相同厚度的均勻厚度板(CFRP-uniform)作為對比,在同一位置處粘貼相同的阻尼層。CFRP-ABH板和CFRP-uniform板均處于自由邊界狀態(tài),在板的(154,?35) mm位置處施加單位恒力,掃頻范圍為10~3 000 Hz。選用C3D20R網(wǎng)格單元對模型進行網(wǎng)格劃分,單元尺寸滿足精度要求。

    提取系統(tǒng)的原點響應(yīng)和平均振速曲線,如圖10所示,由曲線可以得出:在200~3 000 Hz的寬頻帶內(nèi),CFRP-ABH板的振動水平相對于CFRP-uniform均有非常明顯的抑制,且在大多數(shù)頻率下均有10~20 dB的降幅。而在200 Hz頻率以下,CFRP-ABH板相對沒有減振效果,這是由ABH原理決定的,計算得出從波的角度定義的起始頻率(Cut-on frequency)為310 Hz。CFRP-ABH結(jié)構(gòu)配合使用少量阻尼材料能夠?qū)崿F(xiàn)很好的減振效果,表明CFRP-ABH結(jié)構(gòu)也具有優(yōu)越的高頻減振性能和寬頻帶減振性能。

    圖10 仿真計算的動力學(xué)響應(yīng)Fig. 10 Dynamic responses of simulation calculation

    4 實驗驗證

    4.1 能量聚集特性驗證

    圖11 LDV測試系統(tǒng)Fig. 11 LDV test system

    搭建如圖11所示的基于激光超聲掃描技術(shù)的LDV(Laser Doppler Vibrometer)測試系統(tǒng),系統(tǒng)設(shè)備包含計算機、激光頭、信號采集卡(PXI-5105-1)、AE(Acoustic Emission)傳感器和實驗件等。為了驗證復(fù)合材料結(jié)構(gòu)中的ABH效應(yīng),需測試彎曲波在CFRP-ABH結(jié)構(gòu)中的傳播和聚集過程,本文中的實驗件為0/90鋪層(正交)CFRP-ABH板,其整體尺寸及ABH截面參數(shù)與第2.1節(jié)的仿真模型設(shè)置保持一致。定義包含ABH區(qū)域的二維矩形(345 mm×340 mm)掃描區(qū)域,如圖12所示(紅色虛線標(biāo)出),用激光掃描板的平整面,AE傳感器位于板中線上,距離矩形掃描區(qū)域3 mm。

    提取彎曲波在進入ABH區(qū)域過程中的波場(如t=0.23 ms,如圖13(a)所示),可以發(fā)現(xiàn):隨著彎曲波向前傳播,結(jié)構(gòu)的厚度逐漸減小,波速降低,波長被壓縮,波動幅值增大,且波傳播方向發(fā)生偏轉(zhuǎn)。為觀察CFRP-ABH結(jié)構(gòu)的能量聚集特性,提取彎曲波在完全聚集時(對應(yīng)時間t=0.5 ms)的波場,如圖13(b)所示,可以看出,彎曲波大部分聚集在ABH中心平臺靠后的位置(圖中紅色虛線標(biāo)出),與第2.2節(jié)的仿真分析結(jié)論相符。

    4.2 寬頻減振特性驗證

    搭建如圖14所示的激光測振實驗平臺,設(shè)備系統(tǒng)包括多普勒掃描式激光測振儀(Polytec Scanning Vibrometer,PSV)(PSV500)、功 率 放大器(B&K2706)、激振器(B&K4809)、力傳感器(B&K8230)和2塊CFRP薄板。用吊繩懸掛實驗件模擬自由邊界,其余設(shè)置保持和第3節(jié)的仿真條件保持一致,分別測試CFRP-ABH板和CFRP-uniform板的原點響應(yīng)和板平面的平均振速,如圖15所示。

    對比分析CFRP-uniform板和CFRP-ABH板的原點響應(yīng)和平均振速響應(yīng),可以得出:在200 Hz以上頻率范圍,CFRP-ABH板相對于CFRP-uniform板的共振峰峰值有5~18 dB的降低,且除個別頻率下峰值降低不明顯,大多數(shù)頻率下的峰值降幅都超過10 dB。在2 500 Hz左右,CFRP-ABH板的減振效果減弱,相對CFRP均勻板僅有4 dB的降低。相對其他頻段的減振性能有所降低,其原因是從波的角度分析,聲學(xué)黑洞效應(yīng)存在起始頻率,在該頻率以上的范圍內(nèi)都會起作用。但對于一個有限尺寸的、特定的聲學(xué)黑洞結(jié)構(gòu),結(jié)構(gòu)的振動特性與結(jié)構(gòu)的模態(tài)振型相關(guān),在起始頻率以上的頻率范圍內(nèi),結(jié)構(gòu)的陣型以聲學(xué)黑洞區(qū)域的局部模態(tài)為主,此時結(jié)構(gòu)表現(xiàn)出的阻尼較高,減振效果非常明顯;但在高頻范圍內(nèi)存在個別模態(tài)并非表現(xiàn)為局部模態(tài),此時結(jié)構(gòu)表現(xiàn)出的阻尼比沒有大幅提高,相應(yīng)的減振效果并不理想。在低于200 Hz的低頻范圍,ABH結(jié)構(gòu)的引入沒能實現(xiàn)減振效果。實驗得出的結(jié)論與仿真結(jié)果相符,驗證了CFRP-ABH結(jié)構(gòu)的寬頻減振性能。

    5 結(jié) 論

    本文分別對CFRP-ABH薄板進行時域分析和頻域分析,研究了結(jié)構(gòu)的波傳播特性和動力學(xué)特性,得出以下結(jié)論:

    1)通過時域分析,對結(jié)構(gòu)中的波聚集現(xiàn)象進行量化計算,驗證了CFRP-ABH結(jié)構(gòu)中的能量聚集特性;結(jié)合波場測試實驗,驗證了復(fù)合材料結(jié)構(gòu)中的ABH效應(yīng)。

    2)探究了不同鋪層對波幅值能量的分布位置和聚集程度的影響,表明復(fù)合材料結(jié)構(gòu)的各向異性特征對波操控有著重要影響,造成這些影響的本質(zhì)是波速的變化。根據(jù)各向異性對復(fù)合材料結(jié)構(gòu)中ABH效應(yīng)的影響規(guī)律,利用復(fù)合材料特有的鋪層結(jié)構(gòu)設(shè)計以及ABH形狀設(shè)計,可以實現(xiàn)波在特定區(qū)域的聚集,為實現(xiàn)波操控效果的復(fù)合材料ABH結(jié)構(gòu)設(shè)計提供了重要參考。

    3)通過頻域分析,從有限元仿真和實驗測試2方面分析了CFRP-ABH結(jié)構(gòu)的動力學(xué)特性,驗證了其優(yōu)越的寬頻減振性能,表明了ABH結(jié)構(gòu)在復(fù)合材料工程結(jié)構(gòu)的減振降噪等方面有很好的應(yīng)用前景。

    猜你喜歡
    波場鋪層波速
    基于實測波速探討地震反射波法超前預(yù)報解譯標(biāo)志
    彈性波波場分離方法對比及其在逆時偏移成像中的應(yīng)用
    衛(wèi)星天線復(fù)合材料框架的鋪層優(yōu)化設(shè)計
    復(fù)合材料軸結(jié)構(gòu)力學(xué)性能預(yù)測及鋪層方案設(shè)計
    交錯網(wǎng)格與旋轉(zhuǎn)交錯網(wǎng)格對VTI介質(zhì)波場分離的影響分析
    基于Hilbert變換的全波場分離逆時偏移成像
    吉林地區(qū)波速比分布特征及構(gòu)造意義
    CPD鋪層自動創(chuàng)建技術(shù)
    FIBERSIM軟件在工字梁鋪層設(shè)計上的應(yīng)用
    教練機(2015年2期)2015-04-03 06:51:48
    旋轉(zhuǎn)交錯網(wǎng)格VTI介質(zhì)波場模擬與波場分解
    亚洲精华国产精华精| 欧美激情在线99| 国产成+人综合+亚洲专区| 久久久国产成人精品二区| 国产激情久久老熟女| 亚洲熟妇中文字幕五十中出| 日韩欧美 国产精品| 久久午夜综合久久蜜桃| 18禁裸乳无遮挡免费网站照片| 91麻豆精品激情在线观看国产| 五月玫瑰六月丁香| 欧美+亚洲+日韩+国产| 黑人巨大精品欧美一区二区mp4| 精品日产1卡2卡| 蜜桃久久精品国产亚洲av| 最新美女视频免费是黄的| 国产亚洲精品av在线| 亚洲一区二区三区不卡视频| 日韩精品中文字幕看吧| 1000部很黄的大片| 神马国产精品三级电影在线观看| 女人高潮潮喷娇喘18禁视频| av欧美777| 久久天堂一区二区三区四区| 成年女人永久免费观看视频| 国产高清激情床上av| 99精品欧美一区二区三区四区| 国产精品av久久久久免费| 淫秽高清视频在线观看| 我的老师免费观看完整版| 女人被狂操c到高潮| 午夜福利在线观看吧| 成人性生交大片免费视频hd| 美女cb高潮喷水在线观看 | 在线观看日韩欧美| 久99久视频精品免费| 热99re8久久精品国产| aaaaa片日本免费| 成年免费大片在线观看| 欧美成人免费av一区二区三区| 欧美中文综合在线视频| 亚洲精品乱码久久久v下载方式 | 国产探花在线观看一区二区| 久久久久久久久中文| 无限看片的www在线观看| 国产精品亚洲一级av第二区| 久久欧美精品欧美久久欧美| 小说图片视频综合网站| 亚洲中文字幕日韩| av在线蜜桃| 久久久久久九九精品二区国产| 99国产极品粉嫩在线观看| 国产亚洲欧美在线一区二区| x7x7x7水蜜桃| 国产亚洲精品av在线| 精品人妻1区二区| 久久久久久久午夜电影| 首页视频小说图片口味搜索| 亚洲av五月六月丁香网| 国产高清激情床上av| 老鸭窝网址在线观看| 欧美极品一区二区三区四区| 白带黄色成豆腐渣| 午夜免费观看网址| 99在线视频只有这里精品首页| 高潮久久久久久久久久久不卡| 国产伦一二天堂av在线观看| 99精品久久久久人妻精品| 女生性感内裤真人,穿戴方法视频| 亚洲aⅴ乱码一区二区在线播放| 在线观看美女被高潮喷水网站 | 99久久成人亚洲精品观看| 日日摸夜夜添夜夜添小说| 白带黄色成豆腐渣| 在线观看66精品国产| 国产亚洲精品一区二区www| 天天一区二区日本电影三级| 老熟妇乱子伦视频在线观看| 国产亚洲精品久久久久久毛片| 久久精品91蜜桃| 好男人电影高清在线观看| 女人高潮潮喷娇喘18禁视频| 国产主播在线观看一区二区| 宅男免费午夜| 成人特级av手机在线观看| or卡值多少钱| 日日干狠狠操夜夜爽| 在线观看午夜福利视频| 美女cb高潮喷水在线观看 | 亚洲av日韩精品久久久久久密| 可以在线观看的亚洲视频| 国产伦精品一区二区三区四那| 老熟妇仑乱视频hdxx| 国产精品九九99| 无人区码免费观看不卡| 国产精品乱码一区二三区的特点| 美女大奶头视频| av视频在线观看入口| 亚洲精品乱码久久久v下载方式 | www.999成人在线观看| 又大又爽又粗| 夜夜躁狠狠躁天天躁| 亚洲一区二区三区色噜噜| 波多野结衣高清作品| 欧美色视频一区免费| 一级毛片精品| 在线视频色国产色| 成人无遮挡网站| 中文字幕最新亚洲高清| 午夜福利在线观看免费完整高清在 | 精品一区二区三区av网在线观看| 999久久久国产精品视频| 久久久国产精品麻豆| 亚洲欧美激情综合另类| 国产毛片a区久久久久| 婷婷六月久久综合丁香| 啪啪无遮挡十八禁网站| 桃色一区二区三区在线观看| 国产精品乱码一区二三区的特点| 久久精品国产亚洲av香蕉五月| 国产成人精品久久二区二区91| 18禁裸乳无遮挡免费网站照片| 男女那种视频在线观看| 亚洲专区国产一区二区| 欧美在线黄色| 99国产综合亚洲精品| 亚洲人与动物交配视频| 18禁裸乳无遮挡免费网站照片| 精华霜和精华液先用哪个| 精品国产三级普通话版| 亚洲av免费在线观看| 欧美性猛交╳xxx乱大交人| 欧美黑人巨大hd| 久久这里只有精品19| 99热6这里只有精品| 校园春色视频在线观看| 午夜激情福利司机影院| 老鸭窝网址在线观看| 久久久久久久午夜电影| 一夜夜www| 午夜免费激情av| 国产69精品久久久久777片 | 日韩免费av在线播放| 国产伦精品一区二区三区视频9 | 亚洲av美国av| 舔av片在线| 国产亚洲精品久久久com| 欧美一级毛片孕妇| 一级毛片精品| 女同久久另类99精品国产91| 人人妻人人澡欧美一区二区| 日本精品一区二区三区蜜桃| 亚洲aⅴ乱码一区二区在线播放| bbb黄色大片| 99久久精品一区二区三区| 婷婷精品国产亚洲av在线| 亚洲av熟女| 精品久久久久久久久久久久久| 人妻夜夜爽99麻豆av| 国产一区二区三区在线臀色熟女| 亚洲天堂国产精品一区在线| av中文乱码字幕在线| 免费搜索国产男女视频| 国产av麻豆久久久久久久| 久久久国产欧美日韩av| 午夜精品在线福利| 曰老女人黄片| 色综合站精品国产| 国产人伦9x9x在线观看| 男女午夜视频在线观看| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 热99在线观看视频| 国产欧美日韩一区二区精品| 亚洲无线观看免费| 岛国视频午夜一区免费看| 国内精品久久久久久久电影| av天堂中文字幕网| 国产精品日韩av在线免费观看| 久久久国产成人精品二区| 国产精品亚洲一级av第二区| 黄频高清免费视频| 一级毛片精品| 午夜福利欧美成人| 午夜激情福利司机影院| 亚洲成人免费电影在线观看| or卡值多少钱| 99热这里只有是精品50| 午夜a级毛片| 欧美色视频一区免费| 亚洲专区字幕在线| 波多野结衣高清作品| 在线看三级毛片| 一个人看的www免费观看视频| 天堂√8在线中文| 非洲黑人性xxxx精品又粗又长| 熟妇人妻久久中文字幕3abv| 日本 av在线| 日本一二三区视频观看| 给我免费播放毛片高清在线观看| 国内精品一区二区在线观看| 亚洲av成人av| 美女高潮的动态| 国产精品1区2区在线观看.| 久久精品国产亚洲av香蕉五月| 女警被强在线播放| 欧美日韩瑟瑟在线播放| 久久香蕉国产精品| 91麻豆av在线| 国产高潮美女av| 久久久久久大精品| 国产在线精品亚洲第一网站| 国产高潮美女av| 亚洲七黄色美女视频| www.熟女人妻精品国产| 亚洲欧美日韩东京热| 欧美又色又爽又黄视频| 欧美黑人巨大hd| 真实男女啪啪啪动态图| 成年女人看的毛片在线观看| 亚洲性夜色夜夜综合| 99久久成人亚洲精品观看| 国产一区二区三区视频了| 免费看美女性在线毛片视频| 人妻久久中文字幕网| 国产成人精品久久二区二区免费| 欧美黑人欧美精品刺激| 亚洲一区高清亚洲精品| 精品国产超薄肉色丝袜足j| 国产精品九九99| 法律面前人人平等表现在哪些方面| 久久香蕉精品热| 国产精品精品国产色婷婷| 成人鲁丝片一二三区免费| 999久久久精品免费观看国产| 国产亚洲欧美在线一区二区| 欧美色欧美亚洲另类二区| 亚洲熟女毛片儿| 亚洲精品美女久久av网站| 欧美丝袜亚洲另类 | 久久久久久久久久黄片| 国产成人av教育| 国产精品国产高清国产av| 国产精品久久久久久精品电影| 国产成人一区二区三区免费视频网站| 观看免费一级毛片| 亚洲中文字幕日韩| 久久久成人免费电影| 人人妻,人人澡人人爽秒播| 99精品欧美一区二区三区四区| 欧美黑人巨大hd| 久久久国产成人免费| 两人在一起打扑克的视频| 久久久久精品国产欧美久久久| 男女那种视频在线观看| 亚洲午夜理论影院| 男女午夜视频在线观看| 国产精品一及| 国产精品日韩av在线免费观看| 热99re8久久精品国产| 桃色一区二区三区在线观看| 欧美大码av| 精品不卡国产一区二区三区| www日本黄色视频网| 国产午夜精品久久久久久| 熟女人妻精品中文字幕| 欧美国产日韩亚洲一区| 黄色成人免费大全| 国产精品av久久久久免费| 国产精品免费一区二区三区在线| 亚洲成人中文字幕在线播放| 免费高清视频大片| 香蕉丝袜av| 给我免费播放毛片高清在线观看| 成熟少妇高潮喷水视频| 日韩精品中文字幕看吧| 亚洲国产精品合色在线| 我的老师免费观看完整版| 久久久久久大精品| 国产精品久久视频播放| 99久久综合精品五月天人人| www日本在线高清视频| 国产一级毛片七仙女欲春2| 草草在线视频免费看| 免费观看人在逋| 免费看日本二区| 欧美最黄视频在线播放免费| 色哟哟哟哟哟哟| 久久久久久大精品| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品久久男人天堂| 国产精品亚洲一级av第二区| 一本精品99久久精品77| 国产精品乱码一区二三区的特点| 国产一区二区三区视频了| 国产一区二区在线av高清观看| 日本免费a在线| 18禁裸乳无遮挡免费网站照片| 桃色一区二区三区在线观看| 美女被艹到高潮喷水动态| 一级黄色大片毛片| 国产精品av视频在线免费观看| 亚洲无线在线观看| 国内精品久久久久久久电影| 可以在线观看毛片的网站| 国内精品美女久久久久久| 久久精品国产亚洲av香蕉五月| 国产高清有码在线观看视频| e午夜精品久久久久久久| 无限看片的www在线观看| 日韩欧美精品v在线| 亚洲性夜色夜夜综合| 欧美丝袜亚洲另类 | 99久久久亚洲精品蜜臀av| 精品人妻1区二区| 日韩欧美三级三区| 亚洲欧美日韩高清专用| 久久久久久久久免费视频了| 我的老师免费观看完整版| 亚洲av电影在线进入| 在线看三级毛片| 国产精品野战在线观看| 亚洲av五月六月丁香网| 日韩欧美精品v在线| 欧美性猛交黑人性爽| 日本黄色片子视频| 欧美高清成人免费视频www| 伊人久久大香线蕉亚洲五| 亚洲av免费在线观看| 无人区码免费观看不卡| 99热精品在线国产| 97碰自拍视频| 久久久久久久精品吃奶| 久久精品国产99精品国产亚洲性色| 国产乱人伦免费视频| 欧美黄色片欧美黄色片| 午夜亚洲福利在线播放| 欧美黄色淫秽网站| 男女那种视频在线观看| 床上黄色一级片| av在线天堂中文字幕| 热99re8久久精品国产| 亚洲一区二区三区不卡视频| 精华霜和精华液先用哪个| 偷拍熟女少妇极品色| 欧美3d第一页| 国产亚洲精品综合一区在线观看| 亚洲精品在线美女| 长腿黑丝高跟| 久久久国产成人免费| 99久国产av精品| 成人国产综合亚洲| 婷婷六月久久综合丁香| 搡老妇女老女人老熟妇| 热99re8久久精品国产| 免费看光身美女| 国产三级中文精品| 成人特级av手机在线观看| 99久久综合精品五月天人人| 免费在线观看视频国产中文字幕亚洲| 丁香欧美五月| 免费看十八禁软件| 亚洲在线自拍视频| 免费在线观看影片大全网站| 九九热线精品视视频播放| tocl精华| 午夜日韩欧美国产| 午夜福利欧美成人| 51午夜福利影视在线观看| 久9热在线精品视频| 一个人看视频在线观看www免费 | 亚洲天堂国产精品一区在线| 亚洲av片天天在线观看| 女人高潮潮喷娇喘18禁视频| 欧美黑人巨大hd| 国产久久久一区二区三区| 亚洲国产欧美网| 午夜a级毛片| cao死你这个sao货| 桃色一区二区三区在线观看| 国产伦人伦偷精品视频| 欧美中文综合在线视频| 国产伦一二天堂av在线观看| 国产成人系列免费观看| 成年人黄色毛片网站| 窝窝影院91人妻| 欧美+亚洲+日韩+国产| 香蕉丝袜av| 亚洲国产中文字幕在线视频| 午夜视频精品福利| 精品日产1卡2卡| 变态另类成人亚洲欧美熟女| 成人特级av手机在线观看| av片东京热男人的天堂| 熟妇人妻久久中文字幕3abv| 精品人妻1区二区| 99国产精品一区二区蜜桃av| 国产精品永久免费网站| 小说图片视频综合网站| 国产一区二区三区在线臀色熟女| 亚洲av熟女| 一夜夜www| 在线播放国产精品三级| 亚洲熟女毛片儿| 国产高清视频在线观看网站| 免费大片18禁| 日韩大尺度精品在线看网址| 亚洲中文字幕一区二区三区有码在线看 | 99国产极品粉嫩在线观看| 亚洲美女视频黄频| 18美女黄网站色大片免费观看| 国产在线精品亚洲第一网站| ponron亚洲| 成人午夜高清在线视频| 蜜桃久久精品国产亚洲av| 免费一级毛片在线播放高清视频| 在线视频色国产色| 男女视频在线观看网站免费| 伦理电影免费视频| 成人国产一区最新在线观看| 日韩欧美免费精品| 美女 人体艺术 gogo| h日本视频在线播放| 亚洲18禁久久av| 特大巨黑吊av在线直播| 一级黄色大片毛片| 国产高清视频在线观看网站| 熟女少妇亚洲综合色aaa.| 久久精品人妻少妇| 色视频www国产| 精品电影一区二区在线| 国内精品久久久久久久电影| 97碰自拍视频| 特大巨黑吊av在线直播| 18禁黄网站禁片免费观看直播| www.精华液| 欧美av亚洲av综合av国产av| 国产激情久久老熟女| 天天躁狠狠躁夜夜躁狠狠躁| 欧美一级毛片孕妇| 久久久成人免费电影| 欧美极品一区二区三区四区| 香蕉av资源在线| 欧美中文综合在线视频| 国产精品一区二区三区四区久久| 国产又黄又爽又无遮挡在线| 男人舔女人的私密视频| 婷婷亚洲欧美| 国产午夜精品论理片| 婷婷六月久久综合丁香| 欧美日韩亚洲国产一区二区在线观看| 黄色日韩在线| 岛国在线免费视频观看| 激情在线观看视频在线高清| 久久这里只有精品中国| 欧美色欧美亚洲另类二区| 国产精品1区2区在线观看.| 一个人免费在线观看电影 | 神马国产精品三级电影在线观看| 欧美在线一区亚洲| 免费在线观看日本一区| 国产精品,欧美在线| 中文字幕精品亚洲无线码一区| 日本黄色视频三级网站网址| av在线天堂中文字幕| 免费在线观看日本一区| 一级a爱片免费观看的视频| 高清毛片免费观看视频网站| 国产免费男女视频| 欧美黄色片欧美黄色片| 99久久精品国产亚洲精品| 999精品在线视频| 免费观看人在逋| 亚洲av成人一区二区三| 黄片小视频在线播放| 欧洲精品卡2卡3卡4卡5卡区| 久久精品人妻少妇| 色播亚洲综合网| 好看av亚洲va欧美ⅴa在| 欧美午夜高清在线| 法律面前人人平等表现在哪些方面| 免费在线观看影片大全网站| 性色avwww在线观看| 亚洲国产日韩欧美精品在线观看 | 人妻夜夜爽99麻豆av| 久久精品91无色码中文字幕| 老司机在亚洲福利影院| 高清毛片免费观看视频网站| 动漫黄色视频在线观看| 午夜激情福利司机影院| 亚洲精品在线观看二区| 日本在线视频免费播放| 久久香蕉国产精品| 日韩欧美在线乱码| 中亚洲国语对白在线视频| 成人av一区二区三区在线看| 18禁黄网站禁片午夜丰满| 欧美在线黄色| 日本三级黄在线观看| 久久久水蜜桃国产精品网| 91老司机精品| 床上黄色一级片| 不卡av一区二区三区| 亚洲av熟女| 999精品在线视频| 欧美日韩福利视频一区二区| 亚洲成人中文字幕在线播放| 国产69精品久久久久777片 | 亚洲国产精品久久男人天堂| 久久99热这里只有精品18| 国产探花在线观看一区二区| 亚洲国产欧美人成| 国产乱人伦免费视频| 亚洲精品中文字幕一二三四区| 成年女人看的毛片在线观看| 国产欧美日韩一区二区精品| 91字幕亚洲| 一卡2卡三卡四卡精品乱码亚洲| 男人舔女人下体高潮全视频| 国产激情欧美一区二区| 噜噜噜噜噜久久久久久91| 亚洲欧美精品综合一区二区三区| 搡老妇女老女人老熟妇| 成年女人看的毛片在线观看| 欧美黄色淫秽网站| av欧美777| 免费一级毛片在线播放高清视频| 亚洲性夜色夜夜综合| 精品不卡国产一区二区三区| 久久这里只有精品19| 天天躁狠狠躁夜夜躁狠狠躁| 日本一本二区三区精品| 亚洲熟妇中文字幕五十中出| 少妇熟女aⅴ在线视频| а√天堂www在线а√下载| 国产午夜福利久久久久久| 女人被狂操c到高潮| 国产主播在线观看一区二区| 久久久国产成人精品二区| 午夜福利欧美成人| av黄色大香蕉| 亚洲精品美女久久av网站| 国产精品野战在线观看| 国产精品女同一区二区软件 | 动漫黄色视频在线观看| 成人鲁丝片一二三区免费| 一区二区三区高清视频在线| 美女扒开内裤让男人捅视频| 国产精品日韩av在线免费观看| 色综合婷婷激情| 精品免费久久久久久久清纯| 日日夜夜操网爽| 亚洲精品国产精品久久久不卡| 亚洲av成人一区二区三| 嫩草影院精品99| 日日夜夜操网爽| 老鸭窝网址在线观看| 精品国产亚洲在线| 夜夜夜夜夜久久久久| 99热只有精品国产| 日日摸夜夜添夜夜添小说| 国产精品 国内视频| 亚洲精品在线观看二区| 熟女少妇亚洲综合色aaa.| 久久久久久久午夜电影| 亚洲中文日韩欧美视频| 亚洲av第一区精品v没综合| 1000部很黄的大片| 夜夜夜夜夜久久久久| 亚洲18禁久久av| 久久久久国产精品人妻aⅴ院| www.www免费av| 午夜福利18| 日本在线视频免费播放| 99国产综合亚洲精品| 不卡一级毛片| 久久这里只有精品19| 一区福利在线观看| 成人18禁在线播放| 国产精品av视频在线免费观看| 波多野结衣高清作品| 天堂影院成人在线观看| 国产精品av视频在线免费观看| 国产在线精品亚洲第一网站| 国产精品亚洲一级av第二区| 亚洲欧美日韩卡通动漫| 成人特级黄色片久久久久久久| 天天躁日日操中文字幕| 国产精品一区二区精品视频观看| 一个人免费在线观看电影 | 久久久久久九九精品二区国产| 桃红色精品国产亚洲av| 国内精品美女久久久久久| 日本一二三区视频观看| 全区人妻精品视频| 桃红色精品国产亚洲av| 免费在线观看影片大全网站| 亚洲欧美日韩高清专用| 老汉色∧v一级毛片| 这个男人来自地球电影免费观看| 午夜视频精品福利| 一个人免费在线观看电影 | 久久精品夜夜夜夜夜久久蜜豆| 国产又黄又爽又无遮挡在线| 色噜噜av男人的天堂激情| www.精华液| 日韩 欧美 亚洲 中文字幕| 黄色丝袜av网址大全| 91av网一区二区| 久久久久精品国产欧美久久久| 国产精品久久久av美女十八| 在线播放国产精品三级| 精品国产美女av久久久久小说| 99热精品在线国产| 麻豆国产97在线/欧美|