999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

水下夾芯復合空腔結構聲學特性試驗研究

2014-08-11 14:50:07忠,
振動與沖擊 2014年8期
關鍵詞:結構模型

羅 忠, 周 欣

(1. 海軍工程大學 科研部, 武漢 430033;2. 海裝駐上海地區軍事代表局, 上海 201206)

水下夾芯復合空腔結構聲學特性試驗研究

羅 忠1, 周 欣2

(1. 海軍工程大學 科研部, 武漢 430033;2. 海裝駐上海地區軍事代表局, 上海 201206)

引入空腔改善夾芯復合結構低頻段吸聲性能。建立了三種典型水下聲隱身結構聲學系統模型,在數值分析吸聲性能影響因素基礎上,運用脈沖聲管試驗方法,分析了空氣背襯“硬”邊界條件和水背襯“軟”邊界條件下,有無空腔試樣、不同空腔間距、空腔形狀和深度對諧振吸聲峰值及峰值頻率等的影響。試驗結果表明,空腔顯著改善了夾芯復合結構低頻段吸聲性能;空氣背襯下,空腔間距越小,首階諧振吸聲峰值越大,峰值頻率越低,水背襯條件下則相反;空腔體積越大,首階諧振吸聲峰值越小,峰值頻率越低。

夾芯復合結構;空腔;吸聲系數;諧振;試驗

隨著潛艇巡航深度的增大和聲納探測技術的不斷發展,對潛艇的聲隱身技術提出了新的挑戰,水下聲隱身結構正朝著耐壓、低頻和寬頻段吸收的方向發展。夾芯復合聲隱身結構是一種新型水下聲隱身結構,由比強度高、比剛度大、阻尼性能好、透聲性能優異、成型工藝簡單的纖維增強復合材料表層和吸聲高分子芯材組成,可以代替原有“鋼結構+吸聲覆蓋”的吸聲結構,實現水下目標承載和聲隱身的一體化設計[1-3]。在單一均質材料中引入諧振空腔,能夠較好地改善低頻段吸聲性能,帶周期陣列空腔消聲覆蓋層已得到了廣泛應用[2-9]。傳統的對周期空腔陣列消聲覆蓋層聲學性能分析方法有傳遞矩陣法[3-4]、有限元法[2,4-7]、等效參數法[8-9]等,但由于夾芯復合空腔結構中有多種吸聲機理并存,數值方法僅能從定性角度對空腔深度、直徑、間距、位置、形狀等幾何參數的影響規律進行分析,駐波聲管和消聲水池試驗成本較大,所以本文采用脈沖聲管試驗方法,結合多種吸聲機理的影響,分析夾芯復合空腔結構幾何參數對吸聲性能的影響,并獲得了較好的低頻段吸聲性能。應該說,這些研究工作主要是新型水下夾芯復合空腔結構的初步探討,或者是在模型試驗基礎上,對理論研究規律的驗證,還需要根據探索出的吸聲機理對夾芯復合空腔結構進行深入的優化設計。

1 理論分析

圖1 典型聲隱身夾芯復合空腔結構Fig.1 The typical acoustic stealth sandwich composite structure with cavum

典型夾芯復合空腔結構見圖1,由玻璃鋼表層、吸聲材料和空腔等組成。水下聲隱身結構聲學系統簡化模型見圖2。圖2(a)模型Ⅰ與單殼體水下目標聲學環境相似,能夠采用有限元方法分析空腔諧振規律,也可以采用脈沖聲管進行試驗驗證。圖2(b)模型Ⅱ為無限水背襯條件,與聲隱身實心舵和穩定翼等凸出部位的實際聲學環境相似,能夠采用有限元方法分析空腔諧振規律;消聲水池和駐泊聲管的背襯均為尖劈吸聲材料,與模型Ⅱ聲學系統相同,但試驗成本較大;脈沖聲管試驗成本較小,但由于聲管直徑較小,難以布設尖劈吸聲材料,模型Ⅱ很難進行脈沖聲管試驗驗證。圖2(c)模型Ⅲ為耐壓殼體和空氣背襯條件,與水下目標船體的實際聲學環境最相近,也易于進行聲管試驗分析;但于吸聲結構與耐壓殼體之間存在較厚的中間水層,透過吸聲材料的聲波會和空氣背襯反射的聲波疊加,在中間水層形成諧振聲波,對理論分析帶空腔吸聲結構的諧振頻率和帶寬時,難以區分聲波疊加后的諧振和空腔引起的諧振, 一般很少采用模型Ⅲ進行帶空腔

吸聲結構的諧振規律分析。因此,在理論分析水下夾芯復合空腔結構聲學特性時,采用模型Ⅰ和模型Ⅱ;在進行脈沖聲管試驗分析時,采用模型Ⅰ和模型Ⅲ。

圖2 水下聲隱身結構聲學系統簡化模型Fig.2 The simple acoustic model of underwater stealth structure

帶空腔吸聲結構的吸聲機理非常復雜,涉及材料粘彈性損耗、空腔諧振、波型轉化、散射等多種因素,關于帶空腔吸聲結構吸聲性能的計算方法較多,有分層介質理論[3-4]、多次散射逼近法[8]、等效參數反演法[9]、基于Block-Floquet周期邊界條件的有限元法[4-6,10-11]等方法。在夾芯復合結構中引入空腔,其諧振吸聲機理由兩部分組成[2]:① 空腔側壁沿結構法線方向或徑向的振動,而且這兩種振動的頻率特性曲線最大值略有偏移,使吸聲結構在聲波入射下得到較寬的諧振吸聲頻帶;② 整個夾芯復合吸聲結構在聲波作用下整體的彎曲振動,又形成另外一個諧振吸聲系統。本文采用經過算例和試驗驗證的有限元計算方法[2],設計夾芯復合空腔結構,并通過脈沖聲管試驗對其聲學性能進行驗證,檢驗不同背襯條件下夾芯復合空腔結構的吸聲特性。

圖3 模型Ⅰ中夾芯復合空腔結構聲學特性頻響Fig.3 The frequency response of acoustic properties for sandwich composite structure with cavum in model Ⅰ

圖4 模型Ⅱ中夾芯復合空腔結構聲學特性頻響Fig.4 The frequency response of acoustic properties for sandwich composite structure with cavum in model Ⅱ

夾芯復合空腔吸聲結構內含10 mm直徑圓柱型空腔,間隔20 mm,正方形網格排列,空腔腔體位于背聲面芯材底部。模型Ⅰ和模型Ⅱ,有無圓柱空腔、空腔間距、空腔形狀對聲學特性頻響見圖3和圖4。可以看出,① 帶空腔的吸聲層聲學特性在3.5 kHz、4.3 kHz、8.2 kHz多個頻帶較寬的諧振吸聲峰值,空腔顯著改善了低頻吸聲特性;② 模型Ⅰ為空氣背襯的“硬”邊界條件,隨著空腔間距的減小, 首階諧振吸聲峰值增大, 峰值頻率往低頻移動約0.2 kHz,這是由于空腔間距減小,單位體積內空腔體積比增大,在全反射邊界條件下模型整體彎曲剛度降低引起的;圓柱形空腔體積最大、剛度最小,圓錐形空腔體積最小、剛度最大,隨著空腔體積增大,模型Ⅰ首階諧振吸聲峰值和帶寬減小,峰值頻率往低頻移動了0.6 kHz,這與文獻[4]研究規律一致。③模型Ⅱ為水背襯的“軟”邊界條件,隨著空腔間距的減小,首階諧振吸聲峰值和帶寬略微減小,峰值頻率往高頻移動了0.7 kHz,與模型Ⅰ規律相反,且變化更敏感,這是由于間距減小、單個空腔剛度增大引起的;模型Ⅱ隨著空腔體積增大,首階諧振吸聲峰值和帶寬減小,峰值頻率往低頻移動了0.3 kHz,與模型Ⅰ規律和影響幅度基本相同。

2 試驗模型

夾芯復合空腔結構聲學性能試驗采用脈沖聲管測試系統進行,由于聲學測試系統采用的是圓柱型試樣,所以無法對空腔的間距進行精確的控制,只能通過調節單個試樣中空腔的個數來控制間距。夾芯復合空腔試樣見表1,試樣直徑均為56.3 mm,上、下表層玻璃鋼厚度均為4 mm,吸聲芯層均采用聚氨酯改性環氧基微珠增強吸聲體,厚度為25 mm。試驗采用聲學系統模型Ⅰ和模型Ⅲ,模型Ⅲ中試樣位于聲管中部,距管口鋼背襯1.8 m。

表1 夾芯復合空腔結構試樣主要參數

夾芯復合空腔結構試驗模型空腔布置如圖5所示,試驗模具和未粘貼玻璃鋼試樣如圖6所示。

圖5 夾芯復合空腔結構試樣腔體分布Fig.5 The distributing of cavum of experimental model of sandwich composite structure

圖6 夾芯復合空腔結構制作模具和試樣Fig.6 The experimental mould and specimen of sandwich composite structure with cavum

采用脈沖法數字測量裝置,由位于脈沖管一端的換能器向管中發射脈沖聲波,聲波經試樣或另一端背襯反射,由同一換能器接收,測量頻率間隔為200 Hz。在測量頻率點上,應用信號采集器先測量并記錄帶試樣的反射信號相對應的電信號,然后測量并記錄與標準反射體的反射電信號,經DFT處理后得到它們的幅度和相位,最后計算反射系數的模和相位,進而計算出吸聲系數。

為驗證試驗結果的穩定性,5#樣品和10#樣品結構參數相同,按模型Ⅰ和模型Ⅲ條件測試吸聲系數見圖7。可以看出:① 模型Ⅰ條件下,兩個樣品在2-6 kHz頻段吸聲系數吻合較好,并均在6 kHz附近出現諧振峰值,在高頻段吸聲系數有一定差異,這是由于試樣為逐個澆注,工藝及環境條件對基體高頻吸聲性能產生一定影響;② 模型Ⅲ條件下,兩個試樣在2 kHz附近均出現較高的諧振峰值,且吻合較好;③ 兩種模型條件下的吸聲系數測試結果較穩定,可以采用此試驗方法分析夾芯復合空腔結構諧振特性。

3 試驗結果分析

3.1 模型Ⅰ試驗結果

模型Ⅰ條件下,均質層狀夾芯復合試樣與帶空腔試樣吸聲系數測試結果比較見圖8。可以看出:含14個圓柱空腔試樣和含14個錐臺空腔試樣吸聲系數頻譜上均出現了明顯的諧振現象,且空腔體積越大,諧振峰值往低頻移動,與理論分析規律一致。當空腔體積足夠大、空腔個數足夠多時,在吸聲系數頻響曲線上能看出完整的諧振吸聲峰值和帶寬,如錐臺型空腔的諧振峰值在5 kHz附近;圓柱型空腔在6 kHz附近出現一個較小的諧振吸聲峰值,在11 kHz附近出現較高的諧振吸聲峰值,且8-12 kHz頻段內只有半個帶寬,可見空腔有效的改善了均質層狀夾芯復合結構低頻吸聲不足。由于圓柱形試樣空腔間距及材料參數原因,造成了峰值位置與數值計算存在一定差異。

圖7 脈沖聲管試驗吸聲系數穩定性比較Fig.7Thestabilityofexperimentalabsorptioncoefficient圖8 模型Ⅰ無空腔與有空腔比較Fig.8ThecontrastofhomosphereandwithcavumSandwichcompositestructureinmodelⅠ

圓柱型、圓錐型、錐臺型空腔間距對吸聲系數頻響分別見圖9、圖10和圖11。可以看出:① 都于圓柱形空腔,空腔間距在低頻段對吸聲性能影響不大,但在5 kHz以上頻段,隨著空腔間距減小,諧振吸聲峰值增大,峰值頻率稍往低頻移動。② 對于圓錐型空腔試樣,隨著空腔間距減小,6 kHz附近諧振峰值變化不大,在11 kHz附近諧振吸聲峰值稍增大,峰值頻率稍往低頻移動,分析這是由于圓錐型空腔體積較小,且圓錐空腔在聲波作用下只有壁面的彎曲振動模式,三角形壁面相對彎曲剛度又較高,彎曲振動變形較小,不能有效的吸收聲能。③ 空腔間距對錐臺型空腔試樣影響顯著,隨著間距減小,吸聲峰值增大,峰值頻率往低頻移動,吸聲帶寬減小。④試驗結果與理論分析規律一致,夾芯復合空腔結構的諧振吸聲性能,除與空腔間間距有關外,還與空腔體積及芯材的材料參數關系較大。

圖12和圖13為模型Ⅰ條件下空腔形狀對吸聲性能頻響。可以看出:當試樣中空腔數為10個時,吸聲系數在6 kHz和11 kHz附近出現兩個諧振吸聲峰值,隨著空腔體積的增大,峰值也增大,峰值頻率稍往低頻移動,吸聲帶寬也增大;圓臺型空腔體積最大,當試樣中空腔數為14個時,能夠看到出現完整的諧振吸聲峰值。

圖9 模型Ⅰ圓柱型腔空腔間距影響Fig.9TheaffectofabsorptioncoefficientfordifferentcolumncavumintervalinmodelⅠ圖10模型Ⅰ圓錐型腔空腔間距影響Fig.10TheaffectofabsorptioncoefficientfordifferenttapercavumintervalinmodelⅠ圖11 模型Ⅰ錐臺型腔空腔間距影響Fig.11TheaffectofabsorptioncoefficientfordifferentconecavumintervalinmodelⅠ

圖12 模型Ⅰ空腔形狀影響(10個空腔)Fig.12TheaffectofabsorptioncoefficientfordifferentcavumshapeinmodelⅠ(10cavum)圖13 模型Ⅰ空腔形狀影響(14個空腔)Fig.13TheaffectofabsorptioncoefficientfordifferentcavumshapeinmodelⅠ(14cavum)圖14 模型Ⅰ圓柱腔空腔深度影響Fig.14TheaffectofabsorptioncoefficientfordifferentcolumncavumdepthinmodelⅠ

圖14為空氣背襯下14個圓柱型空腔的腔體深度對吸聲系數頻響,由圖可見:腔體深度越大,2-12 kHz頻段吸聲性能越優,且諧振吸聲峰值增大,峰值頻率往低頻移動,但諧振帶寬卻減小,這是因為空腔高度增大,側壁的彎曲剛度減小,因此側壁的彎曲振動頻率降低。

3.2 模型Ⅲ試驗結果

模型Ⅲ條件下,均質層狀夾芯復合試樣與帶空腔的夾芯復合試樣吸聲系數試驗結果見如圖15,空腔的引入有效改善了夾芯復合結構的吸聲性能,圓柱型和錐臺型空腔試樣吸聲系數頻響曲線上均出現了諧振吸聲峰值,隨著空腔體積的增大,峰值位置往低頻移動。

圓柱型、圓錐型、錐臺型空腔間距對吸聲系數頻響分別見圖16、圖17和圖18。可以看出,空腔間距減小,諧振吸聲峰值變密,在大部分頻段內試樣吸聲系數均得到提高,這是由于空腔數增多,諧振模態增多而引起的。同時,隨著空腔間距減小,圓錐型空腔試樣首階諧振峰值增大,峰值頻率往高頻移動,諧振帶寬也增大;而體積最大的錐臺型空腔試樣首階諧振峰值減小,且在高頻段吸聲系數較大。

圖15 模型Ⅲ無空腔與有空腔比較Fig.15ThecontrastofhomosphereandwithcavumSandwichcompositestructureinmodelⅢ圖16 模型Ⅲ圓柱型腔空腔間距影響Fig.16TheaffectofabsorptioncoefficientfordifferentcolumncavumintervalinmodelⅢ圖17 模型Ⅲ圓錐型腔空腔間距影響Fig.17TheaffectofabsorptioncoefficientfordifferenttapercavumintervalinmodelⅢ

圖18 模型Ⅲ錐臺型腔空腔間距影響Fig.18TheaffectofabsorptioncoefficientfordifferentconecavumintervalinmodelⅢ圖19 模型Ⅲ空腔形狀影響(10個空腔)Fig.19TheaffectofabsorptioncoefficientfordifferentcavumshapeinmodelⅢ(10cavum)圖20 模型Ⅲ空腔形狀影響(14個空腔)Fig.20TheaffectofabsorptioncoefficientfordifferentcavumshapeinmodelⅢ(14cavum)

圖19和圖20為空腔形狀對吸聲系數試驗頻響。可以看出:與空氣背襯下結果相反,隨著空腔體積的增大,吸聲系數越低,模型Ⅲ條件水背襯下透射也越大,諧振吸聲效率越低,進入吸聲芯層的聲能沒有得到有效的吸收;同時隨著空腔體積的增大,首階諧振峰值減小,峰值頻率往低頻移動,諧振帶寬也減小,與理論分析規律一致。

圖21 模型Ⅲ圓柱腔空腔深度影響Fig. 21 The affect of absorption coefficient for different column cavum depth in model Ⅲ

圖21為圓柱型空腔腔體深度對吸聲系數試驗頻響。由圖可見:隨著腔體深度的增大,首階諧振吸聲峰值頻率往低頻移動,由于2kHz以下頻段吸聲系數未能測出,所以看不出吸聲峰值是否增大,但可以明顯看出次階諧振吸聲峰值增大,在3-9 kHz頻段,吸聲性能得到改善。

4 結 論

(1) 空腔引入有效改善了夾芯復合結構低頻段吸聲性能,夾芯復合空腔結構聲學性能試驗結果與理論分析規律一致。

(2) 空氣背襯“硬”邊界條件下,隨著空腔間距減小,首階諧振吸聲峰值增大,峰值頻率往低頻移動;水背襯“軟”邊界條件下則相反。

(3) 隨著空腔體積增大,首階諧振吸聲峰值和帶寬減小,峰值頻率往低頻移動。

[ 1 ] John H, Di S C. Composite structure[P]. Britain patent: PCT/GB1998 /01377, 1997.

[ 2 ] 羅忠,楊坤. 水下夾芯復合空腔結構聲學特性計算方法研究[J]. 振動與沖擊, 2013, 32(8):147-152. LUO Zhong, YANG Kun. Acoustic property computing method for an underwater sandwich structure with cavum[J]. Journal of Vibration and Shock, 2013, 32(8):147-152.

[ 3 ] 羅忠,朱錫,簡林安,等. 斜入射下水中隱身夾芯復合材料附體結構聲學設計[J]. 振動與沖擊, 2009, 28(5):49-54. LUO Zhong, ZHU Xi, JIAN Lin-an,et al. Acoustic stealth design on the underwater sandwich composite added structure for oblique incidence[J]. Journal of Vibration and Shock, 2009, 28(5): 49-54.

[ 4 ] 王曼. 水聲吸聲覆蓋層理論與實驗研究[D]. 哈爾濱: 哈爾濱工程大學, 2004.

[ 5 ] Gaunaurd G. Comments on absorption mechanisms for water borne sound in Alberich anechoic layers[J]. Ultrasonic, 1985, 23(2): 90-91.

[ 6 ] Lee Y Y, Lee E W M,Ng C F. Sound absorption of a finite flexible micro-perforated panel backed by an cavity[J]. Journal of Sound and Vibration, 2005, 287(1-2): 227-243.

[ 7 ] Sakagami K, Morimoto M, Koike W. A numerical study of double-leaf microperforated panel absorbers[J]. Applied Acoustics, 2006, 67(7): 609-619.

[ 8 ] 何世平, 湯渭霖,何琳,等. 變截面圓柱形空腔覆蓋層吸聲系數的二維近似解[J]. 船舶力學, 2006, 10(1): 120-127. HE Shi-ping, TANG Wei-lin,HE Lin,et al. Analysis of acoustic characteristics of anechoic coating containing varying sectional cylindrical cavity[J]. Journal of Ship Mechanics, 2006, 10(1): 120-127.

[ 9 ] 朱蓓麗, 任克明. 等效參數法研究帶圓柱通道橡膠體的聲學性能[J]. 上海交通大學學報,1997,31(7):20-25. ZHU Bei-li, REN Ke-ming. Analysis of acoustics properties of rubber with column cavity with equivalent parameters[J]. Journal of Shanghai Jiaotong University, 1997, 31(7): 20-25.

[10] Mead D J. A general theory of harmonic wave propagation in linear periodic systems with multiple coupling[J]. Journal of Sound and Vibration, 1973, 27(2): 235-260.

[11] Hladky-Hennion A C, Decarpigny J N. Analysis of the scattering of a plane acoustic wave by a doubly periodic structure using the finite element method: application to Alberich anechoic coatings[J]. Journal of Acoustic Society of America, 1991, 90(6):3356-3367.

Tests for acoustic-stealth characteristics of underwater sandwich composite structures with cava

LUO Zhong1,ZHOU Xin2

(1. Department of Researching, Naval Univ. of Engineering, Wuhan 430033, China;2. Naval Deputation Department in Shanghai, Shanghai 201206, China)

Introducing cava can improve the sound-absorption ability of a sandwich composite structure within a low-frequency range. The acoustic models for 3 kinds of typical underwater acoustic-stealth structures were established. The affect factors of their sound-absorption ability were analyzed numerically. Then, using the test method of pulse sound tube, the effects of samples with cava or without cava, different distances between cava, shape and depth of cava on the peak value of sound-absorption resonance and its frequency of the structures were analyzed under hard boundary condition and air background as well as soft boundary condition and water background. The results showed that the cava obviously improve the sound-absorption ability of sandwich composite structures within a low-frequency range; under air background, the smaller the distance between cava, the larger the peak value of the 1storder sound-absorption resonance and the lower the corresponding frequency; under water background, the situation is the opposite; the bigger the volume of cava, the smaller the peak value of the 1storder sound-absorption resonance and the lower the corresponding frequency.

sandwich composite structure; cavum; absorption coefficient; resonance; experiment

國防重點基金資助項目

2013-07-16 修改稿收到日期:2013-11-28

羅忠 男,講師,博士,1982年生

TN911.7

A

10.13465/j.cnki.jvs.2014.08.022

猜你喜歡
結構模型
一半模型
《形而上學》△卷的結構和位置
哲學評論(2021年2期)2021-08-22 01:53:34
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
論結構
中華詩詞(2019年7期)2019-11-25 01:43:04
新型平衡塊結構的應用
模具制造(2019年3期)2019-06-06 02:10:54
論《日出》的結構
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉換方法初步研究
創新治理結構促進中小企業持續成長
現代企業(2015年9期)2015-02-28 18:56:50
主站蜘蛛池模板: 亚洲成人播放| 精品国产免费观看| 亚洲看片网| 91麻豆精品国产高清在线| 国产日韩丝袜一二三区| 国产女人水多毛片18| 69综合网| 久久情精品国产品免费| 久久夜夜视频| 成人国内精品久久久久影院| 精品天海翼一区二区| a级毛片免费在线观看| 日韩黄色在线| 亚洲欧洲自拍拍偷午夜色| 亚洲国产理论片在线播放| 国产 在线视频无码| 国产超碰在线观看| 欧美伊人色综合久久天天| 日韩国产黄色网站| 国产AV毛片| 老司国产精品视频91| 久久综合伊人 六十路| 深爱婷婷激情网| 在线欧美日韩| 国产屁屁影院| P尤物久久99国产综合精品| 亚洲色成人www在线观看| 试看120秒男女啪啪免费| 2021亚洲精品不卡a| 久久性妇女精品免费| 二级特黄绝大片免费视频大片| 国产欧美日韩综合在线第一| 99久久这里只精品麻豆| 香蕉视频在线观看www| 欧美一区精品| 狠狠综合久久久久综| 日韩无码白| 中文字幕免费视频| 国产精品lululu在线观看| 亚洲视频免| 亚洲成aⅴ人在线观看| 99re在线观看视频| 国产黄在线观看| 国产精品第| 97视频精品全国免费观看| 国产男女免费视频| 麻豆国产在线观看一区二区| 亚洲欧美日韩色图| 亚洲aaa视频| 老司机午夜精品视频你懂的| 黄色成年视频| av无码一区二区三区在线| 久久国产亚洲偷自| 精品欧美一区二区三区久久久| 色精品视频| 亚洲不卡av中文在线| 午夜在线不卡| 亚洲成人手机在线| 成人av专区精品无码国产| 亚洲日韩精品综合在线一区二区| 无码精品福利一区二区三区| 国产真实乱人视频| 国产尹人香蕉综合在线电影 | 国产91蝌蚪窝| 中文字幕中文字字幕码一二区| 国产一区二区色淫影院| 亚洲香蕉伊综合在人在线| 日韩经典精品无码一区二区| 亚洲综合九九| 日本精品一在线观看视频| 色综合综合网| 婷婷激情亚洲| 午夜不卡福利| 午夜影院a级片| 久久精品无码专区免费| 国产农村1级毛片| 亚洲a级在线观看| 国产成人精品男人的天堂| 亚洲天堂日韩av电影| 国产一区二区免费播放| 国产九九精品视频| 成人精品午夜福利在线播放|