• <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)耦合響應(yīng)數(shù)值分析

    2016-04-24 03:18:20林少芬陳清林
    艦船科學(xué)技術(shù) 2016年10期
    關(guān)鍵詞:雙體主應(yīng)力波浪

    李 堯,林少芬,陳清林

    (1. 廣東海洋大學(xué) 工程學(xué)院,廣東 湛江 524000;2. 集美大學(xué) 輪機(jī)工程學(xué)院,福建 廈門 361021)

    頂浪、斜浪中復(fù)合材料雙體艇艇體結(jié)構(gòu)耦合響應(yīng)數(shù)值分析

    李 堯1,林少芬2,陳清林2

    (1. 廣東海洋大學(xué) 工程學(xué)院,廣東 湛江 524000;2. 集美大學(xué) 輪機(jī)工程學(xué)院,福建 廈門 361021)

    采用流固耦合(Fluid-Structure Interaction,F(xiàn)SI)的艇體波浪載荷和結(jié)構(gòu)響應(yīng)的數(shù)值分析方法,對頂浪、斜浪中復(fù)合材料雙體艇結(jié)構(gòu)的動態(tài)響應(yīng)進(jìn)行研究。分別建立了完整的復(fù)合材料艇體有限元模型以及流場模型,基于數(shù)值水池造波技術(shù),通過計算獲得了頂浪、斜浪中復(fù)合材料艇體結(jié)構(gòu)的時域動態(tài)響應(yīng)結(jié)果,選取高應(yīng)力梯度區(qū)域,通過網(wǎng)格加密重構(gòu)同時獲得了復(fù)合材料的層間應(yīng)力。選取具有代表性的前 10 大等效應(yīng)力與內(nèi)外面板最值主應(yīng)力,在將 FSI 與傳統(tǒng)基于經(jīng)驗公式的有限元法(Finite Elementmethod,F(xiàn)EM)的結(jié)果對比中發(fā)現(xiàn),F(xiàn)SI 中拱、中垂的計算結(jié)果更接近于 FEm彎扭組合工況,而采用《鋼制雙體船直接計算指南》計算復(fù)合材料雙體艇時,所用經(jīng)驗公式的頂浪航行波浪載荷計算值偏小。

    耦合響應(yīng);復(fù)合材料雙體艇;數(shù)值水池技術(shù);頂浪;斜浪;主應(yīng)力

    0 引 言

    流固耦合法(FSI)主要研究兩相介質(zhì)間的相互作用,其在航空航天、水利、建筑、海洋、生物等領(lǐng)域都有重要應(yīng)用[1]。在船舶領(lǐng)域,F(xiàn)SI 被作為一種結(jié)構(gòu)強(qiáng)度的計算手段,該法主要基于數(shù)值水池技術(shù)[2],該技術(shù)已被證明能夠細(xì)致描述周圍流場,最終獲得與試驗數(shù)據(jù)較吻合的結(jié)果[2]。

    復(fù)合材料因其在減重、耐腐蝕性、強(qiáng)度方面的卓越性能被認(rèn)為在造船領(lǐng)域有巨大的發(fā)展?jié)摿Γ诖邦I(lǐng)域復(fù)合材料結(jié)構(gòu)多以夾層板三明治結(jié)構(gòu)呈現(xiàn),主要包括輕質(zhì)芯材、面板以及加強(qiáng)結(jié)構(gòu)。復(fù)合材料雙體艇作為一種新艇型,具有重量輕、穩(wěn)性好、甲板面積大等優(yōu)點(diǎn),但其夾層材料、復(fù)雜的結(jié)構(gòu)卻給強(qiáng)度計算帶來了挑戰(zhàn)。目前,國內(nèi)針對其強(qiáng)度計算主要采用規(guī)范公式法,F(xiàn)Em作為計算商船強(qiáng)度的最有效方法[3],在復(fù)合材料多體船艇計算領(lǐng)域還未完全普及,鮮有采用FSI 計算的例子。國外應(yīng)用 FSI 計算復(fù)合材料雙體艇的例子[4]對結(jié)構(gòu)進(jìn)行了大量簡化并且僅計算了頂浪航行工況,未考慮更易造成結(jié)構(gòu)破壞的斜浪航行工況。鑒于目前復(fù)合材料雙體艇的 FEm主要參照鋼質(zhì)雙體船計算指南[5],基于數(shù)值水池的 FSI 計算結(jié)果可以一定程度上對同工況下的 FEm進(jìn)行檢驗。

    本文采用 FSI 對 1 艘 52 ft 復(fù)合材料雙體艇進(jìn)行頂浪、斜浪工況下結(jié)構(gòu)的時域響應(yīng)計算,將流體域模型與有限元詳細(xì)模型分開建立,借助 WORKBENCH 耦合求解器進(jìn)行求解。將結(jié)果與同工況下的 FEm對比,重點(diǎn)關(guān)注前 10 大等效應(yīng)力與主要結(jié)構(gòu)最大、小主應(yīng)力,所得結(jié)論為進(jìn)一步豐富和細(xì)化針對復(fù)合材質(zhì)雙體船的結(jié)構(gòu)數(shù)值計算提供了重要參考。

    1 材料參數(shù)及結(jié)構(gòu)模型

    主要參數(shù)包括面板/芯材材料特性,厚度以及鋪層設(shè)置,面板材料必須要有足夠的強(qiáng)度承受波浪外載荷,芯材則具有較低的密度以達(dá)到減輕整個游艇重量的目的,文中的艇體面板采用強(qiáng)度較高的 E-玻璃纖維正交布、玻璃纖維短切氈、樹脂交替鋪設(shè)的結(jié)構(gòu)形式,芯材為密度 100 kg/m3PVC 泡沫。主要力學(xué)性能如表1 所示,面板數(shù)據(jù)為實測數(shù)據(jù)(試驗樣板見圖1)。

    艇體的有限元材料模型由上下面板、芯材的“三明治”形式模擬,泡沫芯材看成各向同性材料,面板層定義為各向異性材料。主尺度為:艇長 15.83m,設(shè)計水線長 15.3m,總型寬 6.048m,型深 1.596m,片體水線寬 2.5m,計算吃水 0.78m,排水體積 16.6 t,片體中心距 4.45m。由于計算模型為薄殼結(jié)構(gòu),并且具有高柔中間層,較低的橫向剪切剛度,在 WORKBENCH中通過內(nèi)插程序的方法選用 SHELL181 單元,劃分網(wǎng)格后的模型如圖2 所示。通過定義局部坐標(biāo)系的方式,校正材料鋪層方向使其符合實際情況。

    表1 材料的力學(xué)參數(shù)列表Tab. 1 Properties of FRP panels and PVC Foam

    圖1 試驗樣板圖Fig. 1 Standard panels for testing

    2 粘性流場數(shù)值模擬

    圖2 網(wǎng)格模型Fig. 2 Meshmodel

    2.1 數(shù)值造波與消波

    目前的數(shù)值造波方法主要有仿物理造波法和純數(shù)值造波法兩大類,由于仿物理造波需要運(yùn)用動邊界和動網(wǎng)格特性,大大增加了計算時間[6]。本文采用純數(shù)值造波法中的邊界造波法,所造波浪基于線性波浪理論[7],該理論假定波浪沿 ox 軸正向傳播,模擬斜浪航行時,假定波浪的傳播方向與游艇的航行方向 ox 軸正向成 β 角,波浪速度勢經(jīng)坐標(biāo)轉(zhuǎn)換后得到各方向的速度表達(dá)式

    式中:ξa為波幅;k 為波數(shù);H 為水深;ω0為波浪的自然頻率。系數(shù) C0表示為:

    本文采用人工阻尼消波的方法來消除所造波浪,其主要方法是模擬真實波浪水池的消波層,在消波區(qū)域的動量方程中直接添加阻尼項,通過控制外部給定的阻尼大小來實現(xiàn)消波,消波區(qū)位于流場尾域一個波長的范圍。具體方法為加入 UDF 消波程序,給定計算工況下的波高、波長、波數(shù)等初始條件,完成對數(shù)值水池的模擬。

    2.2 CFD 模型

    通過對非耦合區(qū)域的簡化建立流場分析模型,導(dǎo)入 ICEM-CFD 中劃分網(wǎng)格并設(shè)置邊界條件。本文將計算域定義為:縱向船首向前 1 倍船長,船尾向后 2 倍船長;橫向左右舷側(cè)各向兩邊延伸 6 倍片體型寬;豎向從水線處向下延伸 10 倍吃水,從水線處向上延伸 5倍吃水。采用四面體網(wǎng)格為主的混合網(wǎng)格,總網(wǎng)格數(shù)為 1 889 795。邊界條件設(shè)置如圖3 所示,入口為速度入口,出口為自由出口,兩側(cè)及底面為滑移壁面,上表面為對稱邊界。

    圖3 流場計算域示意圖Fig. 3 Fluid domain for hydrodynamic analysis

    3 耦合計算與結(jié)構(gòu)響應(yīng)分析

    3.1 載荷與邊界條件

    根據(jù)實際航行情況,本文選取頂浪(中拱與中垂)、斜浪工況為計算工況,設(shè)計航速 13 kn。其中頂浪航行工況波長為艇水線長,斜浪工況波浪參數(shù)依據(jù)圖4 發(fā)生片體不同步扭轉(zhuǎn)選取,此時波長 8.55m,浪向角 74°,波高取航行區(qū)域有義波高的均值 0.8m。為了模擬自由漂浮狀態(tài),采用慣性釋放的邊界條件。

    3.2 計算結(jié)果與分析

    流固耦合計算的 3 種工況結(jié)果如圖5 所示,中拱、中垂的結(jié)果顯示,頂浪航行時除了產(chǎn)生總縱彎曲外,艇身還產(chǎn)生了小幅橫向彎曲,主要原因是興波向后傳播過程中在兩片體中間產(chǎn)生了回流、相互干擾造成片體兩側(cè)水壓力分布不均。斜浪工況下左右片體應(yīng)力分布并不對稱,主要高應(yīng)力區(qū)分布在右側(cè)片體首部。此時,游艇兩片體產(chǎn)生了不對稱扭轉(zhuǎn),但由于波浪的抨擊、上涌、破碎等作用,由此產(chǎn)生的艇體變形并不具有反對稱性。

    圖4 斜浪航行工況簡圖Fig. 4 Sketch of oblique waves

    圖5 航行等效應(yīng)力云圖Fig. 5 Equivalent stress nephogramfor hogging

    沿夾層板厚度纖維 1 方向的層間應(yīng)力曲線如圖6所示,面板的主應(yīng)力顯著高于泡沫,說明夾層板結(jié)構(gòu)面板主要承受拉壓彎矩,而泡沫非承力結(jié)構(gòu),主要起支撐、提高夾層板剛度以抵抗變形作用。圖7 為該艇在不同航速下(僅列出 2 種)的艇身結(jié)構(gòu)最大、小主應(yīng)力絕對值的時域響應(yīng)圖,選取航行時間區(qū)間 6~12 s,從應(yīng)力曲線的形狀看,該波動呈現(xiàn)出了一定的周期性,該周期性緊隨波浪周期變化但不吻合;此外,監(jiān)測結(jié)果發(fā)現(xiàn),最大最小正應(yīng)力的發(fā)生位置并不一致,其曲線的變化趨勢并不吻合;航速越大,最值主應(yīng)力也越大。

    3.3 與有限元法結(jié)果對比

    文中作為對照的 FEm是參考 CCS 對于雙體船的直接計算指南[5]完成,計算流程符合規(guī)范要求。為更好地考查 FSI 與 FEm的結(jié)果,選取全艇具有代表性的前 10 大應(yīng)力為對象比較 3 種工況的計算結(jié)果(見圖8)。由圖可知,中拱狀態(tài)下的前 10 大等效應(yīng)力中,F(xiàn)EM前兩大應(yīng)力明顯大于 FSI,二者后 8 大應(yīng)力吻合較好,F(xiàn)SI 應(yīng)力的波動幅度小于 FEM;中垂工況下 FSI 計算的前 10 大應(yīng)力值總體大于 FEM,斜浪航行工況有限元法前 10 大等效應(yīng)力值比 FSI 高出一倍,應(yīng)力曲線走勢相近。

    圖6 工況纖維1方向?qū)娱g應(yīng)力曲線Fig. 6 Normal stress distribution in 1 direction for Hogging

    圖7 全艇最大、小正應(yīng)力的時域變化曲線圖Fig. 7 Dynamic response results of the whole ship'smax/Min. principle stress

    圖8 FSI 與 FEA 前 10 大等效應(yīng)力曲線圖Fig. 8 Top 10 von-mises calculated by FSI/FEA for Hogging

    3 種工況下主要結(jié)構(gòu)最大、小主應(yīng)力如圖9和圖10所示,F(xiàn)SI 計算的 3 種工況最大、小主應(yīng)力絕對值在15mPa 以下,數(shù)值上與 FEm計算的中拱、中垂工況接近,斜浪工況 FEm各主應(yīng)力則顯著大于 FSI。中拱狀態(tài)下 FEM、FSI 計算的最大、小主應(yīng)力各結(jié)構(gòu)變化趨勢相似,根據(jù)主應(yīng)力與變形形式的對應(yīng)關(guān)系可知,中拱狀態(tài)游艇的拉伸、壓縮變形分布較均勻;中垂?fàn)顟B(tài)下全艇拉伸變形占主導(dǎo);斜浪工況下 FEM、FSI 總體拉壓變形均有分布,但注意到該工況主甲板、主艇體內(nèi)外面板存在應(yīng)力差,特別是直接受力的主艇體結(jié)構(gòu),內(nèi)外面板應(yīng)力差較大,可能造成夾層板發(fā)生層間破壞。

    中拱狀態(tài)下 2 種方法在主甲板內(nèi)外面板、主艇體外面板的主應(yīng)力極值基本吻合,曲線的走勢相近,計算的全艇最值主應(yīng)力均位于艙壁骨材;中垂?fàn)顟B(tài)下FSI 在主甲板、主艇體的主應(yīng)力均超過 FEM,若考慮FEm的結(jié)果應(yīng)具有一定的安全余量,則通過經(jīng)驗公式計算的中垂外載荷偏小;斜浪工況下 FEm各結(jié)構(gòu)最值主應(yīng)力均顯著高于 FSI,曲線的基本走勢相近,注意到 FEm主艇體最大主應(yīng)力顯著高于主甲板和艙壁骨材,這是由 FEm的計算載荷集中施加在主艇體底部所致??紤]到連接橋結(jié)構(gòu)較片體更加脆弱,而斜浪工況造成的片體不同步扭轉(zhuǎn)對連接橋結(jié)構(gòu)威脅最大,主應(yīng)力應(yīng)具有一定的安全裕度。

    3)艇身結(jié)構(gòu)最值主應(yīng)力隨著航速的升高而變大,其時域變化值呈現(xiàn)一定的周期性,并且最值主應(yīng)力的發(fā)生位置也在不斷變化,但變化周期并不與波浪周期完全吻合。

    4)前 10 大等效應(yīng)力與最值主應(yīng)力曲線顯示,中拱工況 FEm與 FSI 吻合較好,此時艇體拉壓變形分布較均勻;中垂工況 FEm略小于 FSI,全艇拉伸變形略占主導(dǎo);斜浪工況 FEm顯著大于 FSI,主應(yīng)力曲線走勢相比后者更陡??紤]到玻璃鋼存在剛度小、建造工藝不確定等問題,在許用應(yīng)力取 0.3 倍試驗強(qiáng)度的基礎(chǔ)上仍然需要留有一定的安全裕度?!朵撡|(zhì)雙體船直接計算指南》應(yīng)用于文中的復(fù)合材料艇時,據(jù)其計算的中拱、中垂外載荷值偏小,若引入等強(qiáng)度的設(shè)計思想,可以取斜浪工況下 FEm的平均安全系數(shù)作為修正系數(shù)。

    圖9 FSI 與 FEA 最大主應(yīng)力曲線圖Fig. 9 Max. principle stress of FSI/FEA

    圖10 FSI 與 FEA 最小主應(yīng)力曲線圖Fig. 10 Min. principle stress of FSI/FEA

    4 結(jié) 語

    本文采用流固耦合法對復(fù)合材料雙體艇進(jìn)行迎浪和斜浪航行工況下的結(jié)構(gòu)響應(yīng)計算,結(jié)果表明:

    1)FSI 頂浪工況下艇身結(jié)構(gòu)除了產(chǎn)生總縱彎曲外,還有小幅橫向彎曲,F(xiàn)SI 的結(jié)果更接近于 FEm組合工況。

    2)夾層板面板主要承力,而泡沫芯材主要作用為提高剛度抵抗變形。

    [1]錢若軍, 董石麟, 袁行飛. 流固耦合理論研究進(jìn)展[J]. 空間結(jié)構(gòu), 2008, 14(1): 3-15. QIAN Ruo-jun, DONG Shi-lin, YUAN Xing-fei. Advances in research on fluid-structure interaction theory[J]. Spatial Structures, 2008, 14(1): 3-15.

    [2]趙艷, 朱仁慶, 劉珍. 三維數(shù)值波浪水池的構(gòu)建和粘性影響的研究[J]. 艦船科學(xué)技術(shù), 2014, 36(5): 42-48. ZHAO Yan, ZHU Ren-qing, LIU Zhen. Simulation of 3D numerical wave tank and viscosity research[J]. Ship Science and Technology, 2014, 36(5): 42-48.

    [3]方昭昭, 趙丙乾, 朱仁傳. 頂浪中船舶運(yùn)動的數(shù)值模擬與波浪增阻計算[J]. 中國造船, 2014, 55(2): 8-17. FANG Zhao-zhao, ZHAO Bing-qian, ZHU Ren-chuan. Numerical simulation of shipmotion and calculation of added resistance in heading waves[J]. Shipbuilding of China, 2014, 55(2): 8-17.

    [4]牟宗寶, 林焰, 于雁云. 玻璃鋼漁船帽形骨材頂角變化對應(yīng)力集中的影響[J]. 船舶工程, 2012, 34(5): 17-19.mU Zong-bao, LIN Yan, YU Yan-yun. Affect of angle change on stress concentration of cap-shape reinforcing rib in FRP fishing boat[J]. Ship Engineering, 2012, 34(5): 17-19.

    [5]OJEDA R, PRUSTYbG, SALASm. Finite element investigation on the static response ofacomposite catamaran under slamming loads[J]. Ocean Engineering, 2004, 31(7): 901-929.

    [6]中國船級社. 海上高速船入級與建造規(guī)范[S]. 北京: 人民交通出版社, 2012.

    [7]李宏偉. 數(shù)值水池造波方法研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2009.

    [8]方昭昭, 朱仁傳, 繆國平. 數(shù)值波浪水池中航行船舶輻射問題的數(shù)值模擬[J]. 水動力學(xué)研究與進(jìn)展, 2011, 26(1): 65-72. FANG Zhao-zhao, ZHU Ren-chuan,mIAO Guo-ping. Numerical simulation on radiation problems ofmoving vessels in numerical wave tank[J]. Chinese Journal of Hydrodynamics, 2011, 26(1): 65-72.

    Numerical analysis on fluid-solid coupling response of the composite catamaran structures in head and oblique waves

    LI Yao1, LIN Shao-fen2, CHEN Qing-lin2
    (1. School of Engineering, Guangdong Ocean University, Zhanjiang 524000, China; 2. School ofmarine Engineering, Jimei University, Xiamen 361021, China)

    Fluid-structure interactionmethod (FSI) has been used, researching on the dynamic response of the catamaran composite structures in head and oblique waves. The whole 3D compositemodel of catamaran andanumerical wave tank are constructed, based on numerical wave tank technology, dynamic response results of composite catamaran structures in head and oblique waves are carried on. Areas with high stress gradient is chosen and developed with refinedmesh., interlaminar stresses distribution are then determined. Furthermore, top 10 von-mises andmax/Min. principal stress of the whole yacht are selected as comparative objects, comparing the results of the Fluid-Structure Interactionmethod (FSI) with finite elementmethod (FEM) ones, it finds that the results fromhogging, sagging conditions of FSI aremuch closer toacombination of bending and torsion conditions of FEM. It is also found that, using the direct calculationmethod for the steel catamaran to calculate the composite one in the paper under the head wave will lead toarelatively smaller calculation result.

    coupling response;composite catamaran;numerical wave tank technology;head wave;oblique wave;principal stress

    U661.43

    :A

    1672 - 7619(2016)10 - 0041 - 05

    10.3404/j.issn.1672-7619.2016.010.008

    2016 - 02 - 29;

    2016 - 04 - 05

    交通運(yùn)輸部應(yīng)用基礎(chǔ)研究資助項目(2014329815100);福建省高校產(chǎn)學(xué)研重大資助項目(2014H6020);廣東省創(chuàng)新訓(xùn)練資助項目(CXXL2014077)

    李堯(1989 - ),男,講師,研究方向為船舶與海洋結(jié)構(gòu)物強(qiáng)度計算及流固耦合。

    猜你喜歡
    雙體主應(yīng)力波浪
    波浪谷和波浪巖
    雙體清污船綜合減搖自適應(yīng)控制方法
    波浪谷隨想
    雙體小行星系統(tǒng)平衡態(tài)與穩(wěn)定性研究
    基于xflow軟件的雙體游艇阻力準(zhǔn)確性分析
    復(fù)合斷層對地應(yīng)力的影響研究
    去看神奇波浪谷
    深部沿空巷道圍巖主應(yīng)力差演化規(guī)律與控制
    考慮中主應(yīng)力后對隧道圍巖穩(wěn)定性的影響
    波浪中并靠兩船相對運(yùn)動的短時預(yù)報
    中國航海(2014年1期)2014-05-09 07:54:24
    亚洲少妇的诱惑av| 人妻 亚洲 视频| 中文乱码字字幕精品一区二区三区| av一本久久久久| 久久久久视频综合| 欧美精品一区二区大全| 久久av网站| 久久久久久人妻| 丰满少妇做爰视频| 中文字幕av电影在线播放| 国产一区二区三区av在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久热这里只有精品99| 麻豆精品久久久久久蜜桃| 男男h啪啪无遮挡| 美女主播在线视频| 夫妻性生交免费视频一级片| 少妇精品久久久久久久| 亚洲色图 男人天堂 中文字幕 | 国产在线免费精品| 日韩免费高清中文字幕av| 又大又黄又爽视频免费| 日本黄色日本黄色录像| 亚洲成人一二三区av| 午夜福利视频精品| av片东京热男人的天堂| 中文字幕人妻丝袜制服| 亚洲欧美成人精品一区二区| 欧美丝袜亚洲另类| 亚洲精品日本国产第一区| 免费观看性生交大片5| 中文精品一卡2卡3卡4更新| 婷婷色综合大香蕉| 母亲3免费完整高清在线观看 | 国产精品久久久久久精品古装| 国产1区2区3区精品| 国产精品国产三级国产av玫瑰| 一级片免费观看大全| 亚洲精品美女久久av网站| 久久久久久人人人人人| 国产精品久久久久成人av| 午夜激情久久久久久久| 国产永久视频网站| 在线精品无人区一区二区三| 久久精品国产鲁丝片午夜精品| 日韩 亚洲 欧美在线| 国产精品电影一区二区三区 | 久久精品亚洲精品国产色婷小说| 夫妻午夜视频| 69精品国产乱码久久久| 大陆偷拍与自拍| 亚洲精品久久成人aⅴ小说| 国产亚洲一区二区精品| 丁香六月欧美| 国产亚洲欧美精品永久| 久99久视频精品免费| 精品免费久久久久久久清纯 | 成人黄色视频免费在线看| 国产精品一区二区在线观看99| 欧美黑人精品巨大| 丝瓜视频免费看黄片| 中文亚洲av片在线观看爽 | 91成年电影在线观看| 91精品三级在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品中文字幕在线视频| 免费在线观看视频国产中文字幕亚洲| 欧美最黄视频在线播放免费 | 久久亚洲真实| 成人亚洲精品一区在线观看| 黄片小视频在线播放| 亚洲国产毛片av蜜桃av| 精品电影一区二区在线| 又紧又爽又黄一区二区| 亚洲一区二区三区不卡视频| 免费在线观看日本一区| 成人影院久久| 亚洲综合色网址| 99久久国产精品久久久| 国产色视频综合| 日韩人妻精品一区2区三区| 精品国产乱码久久久久久男人| 黄色视频不卡| 亚洲精品一二三| avwww免费| 久久国产精品男人的天堂亚洲| 天天躁夜夜躁狠狠躁躁| 91大片在线观看| 69精品国产乱码久久久| 老鸭窝网址在线观看| 交换朋友夫妻互换小说| 午夜亚洲福利在线播放| 亚洲性夜色夜夜综合| 久久久久精品人妻al黑| 欧美黑人欧美精品刺激| 80岁老熟妇乱子伦牲交| 国产蜜桃级精品一区二区三区 | 黄片大片在线免费观看| 1024视频免费在线观看| 啪啪无遮挡十八禁网站| 国产伦人伦偷精品视频| 两性夫妻黄色片| 久久精品成人免费网站| 黄频高清免费视频| 午夜福利在线免费观看网站| 在线天堂中文资源库| 精品乱码久久久久久99久播| av视频免费观看在线观看| tube8黄色片| 久热这里只有精品99| av线在线观看网站| 超碰97精品在线观看| 久久香蕉精品热| 精品少妇一区二区三区视频日本电影| 国产精品99久久99久久久不卡| 欧美乱码精品一区二区三区| 亚洲 欧美一区二区三区| 中文字幕最新亚洲高清| 亚洲国产中文字幕在线视频| 久久草成人影院| 欧美一级毛片孕妇| 成熟少妇高潮喷水视频| 免费不卡黄色视频| 99热网站在线观看| 极品人妻少妇av视频| 美女 人体艺术 gogo| 很黄的视频免费| 久久这里只有精品19| 免费不卡黄色视频| 欧美丝袜亚洲另类 | 男女床上黄色一级片免费看| 欧美午夜高清在线| 国产精品国产av在线观看| 久久精品国产亚洲av香蕉五月 | 自拍欧美九色日韩亚洲蝌蚪91| 亚洲综合色网址| 自拍欧美九色日韩亚洲蝌蚪91| 校园春色视频在线观看| 午夜成年电影在线免费观看| 久久久久精品国产欧美久久久| 成人18禁在线播放| 国内久久婷婷六月综合欲色啪| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产中文字幕在线视频| 国产精品欧美亚洲77777| 国产精品欧美亚洲77777| 欧美精品av麻豆av| 最近最新中文字幕大全电影3 | 免费不卡黄色视频| 18禁裸乳无遮挡动漫免费视频| 精品午夜福利视频在线观看一区| 亚洲国产中文字幕在线视频| 丰满人妻熟妇乱又伦精品不卡| 91大片在线观看| 精品一区二区三区av网在线观看| 91老司机精品| 日日夜夜操网爽| 日日爽夜夜爽网站| 欧美乱码精品一区二区三区| 99国产精品一区二区三区| 美女国产高潮福利片在线看| 又紧又爽又黄一区二区| 超碰97精品在线观看| 成人国产一区最新在线观看| 亚洲精品在线观看二区| 无遮挡黄片免费观看| 久久国产乱子伦精品免费另类| 久久狼人影院| 国产精品国产av在线观看| 天堂动漫精品| 色94色欧美一区二区| 久久国产精品大桥未久av| 日韩欧美免费精品| 51午夜福利影视在线观看| 亚洲黑人精品在线| 欧美成狂野欧美在线观看| 国产国语露脸激情在线看| 久久久久久亚洲精品国产蜜桃av| 日本撒尿小便嘘嘘汇集6| 国产成人欧美在线观看 | 乱人伦中国视频| 在线视频色国产色| 国产精品国产av在线观看| 国产欧美亚洲国产| 99精国产麻豆久久婷婷| 免费女性裸体啪啪无遮挡网站| 交换朋友夫妻互换小说| 成人手机av| 人人妻人人澡人人看| 婷婷丁香在线五月| 国产人伦9x9x在线观看| 精品少妇久久久久久888优播| 亚洲中文日韩欧美视频| 亚洲专区字幕在线| 女同久久另类99精品国产91| 国产一区二区三区视频了| 欧美精品啪啪一区二区三区| 99国产精品一区二区三区| 欧美成人午夜精品| 男人舔女人的私密视频| 国产真人三级小视频在线观看| 国产精品久久久人人做人人爽| 国产欧美亚洲国产| 18禁美女被吸乳视频| 视频在线观看一区二区三区| 欧美日韩一级在线毛片| 国产成人免费观看mmmm| 久久久久精品国产欧美久久久| 精品卡一卡二卡四卡免费| 午夜91福利影院| 国产在视频线精品| 激情在线观看视频在线高清 | 高清毛片免费观看视频网站 | av超薄肉色丝袜交足视频| 亚洲七黄色美女视频| 91国产中文字幕| 老司机午夜十八禁免费视频| 欧美国产精品一级二级三级| 精品国产乱子伦一区二区三区| 日韩免费高清中文字幕av| 一边摸一边做爽爽视频免费| 乱人伦中国视频| 在线永久观看黄色视频| 色94色欧美一区二区| av网站免费在线观看视频| 99精品久久久久人妻精品| 色综合欧美亚洲国产小说| 亚洲情色 制服丝袜| 高清欧美精品videossex| 天堂俺去俺来也www色官网| 国产日韩欧美亚洲二区| 国产精品久久久人人做人人爽| 成人免费观看视频高清| 日韩人妻精品一区2区三区| 热99re8久久精品国产| 老司机福利观看| 怎么达到女性高潮| 亚洲自偷自拍图片 自拍| 国产一区二区三区视频了| 在线观看一区二区三区激情| 黄片小视频在线播放| 日韩欧美免费精品| 欧美日韩视频精品一区| 91麻豆av在线| 久久青草综合色| 国产精品乱码一区二三区的特点 | 桃红色精品国产亚洲av| 首页视频小说图片口味搜索| 91老司机精品| 亚洲一区中文字幕在线| 日韩免费av在线播放| 久久精品国产亚洲av香蕉五月 | 日日夜夜操网爽| av免费在线观看网站| 自线自在国产av| 久久精品国产清高在天天线| 中文字幕另类日韩欧美亚洲嫩草| 国产av又大| 午夜久久久在线观看| 色综合婷婷激情| 亚洲欧美日韩另类电影网站| 国产亚洲精品第一综合不卡| 免费在线观看完整版高清| 久久精品国产亚洲av高清一级| 亚洲五月婷婷丁香| 久久久久久亚洲精品国产蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品合色在线| 久久久久久人人人人人| 国产精品国产av在线观看| 精品国产乱码久久久久久男人| 国产极品粉嫩免费观看在线| 久久ye,这里只有精品| 亚洲va日本ⅴa欧美va伊人久久| 国产欧美日韩一区二区三区在线| 国产蜜桃级精品一区二区三区 | 美女福利国产在线| 久久久精品区二区三区| 精品视频人人做人人爽| 成在线人永久免费视频| 欧美日韩视频精品一区| 丰满人妻熟妇乱又伦精品不卡| 超碰成人久久| 99久久人妻综合| 国产一区二区三区综合在线观看| 中亚洲国语对白在线视频| 国产成人欧美在线观看 | 婷婷精品国产亚洲av在线 | 国内久久婷婷六月综合欲色啪| 少妇的丰满在线观看| 日日摸夜夜添夜夜添小说| 日日爽夜夜爽网站| 水蜜桃什么品种好| 午夜福利视频在线观看免费| 天堂√8在线中文| 黄片大片在线免费观看| 亚洲综合色网址| 一本大道久久a久久精品| 纯流量卡能插随身wifi吗| 99热只有精品国产| 最近最新免费中文字幕在线| 婷婷丁香在线五月| 亚洲欧美一区二区三区久久| 飞空精品影院首页| 国产不卡av网站在线观看| 日本a在线网址| 夜夜夜夜夜久久久久| 日日夜夜操网爽| 国产精品久久久久久精品古装| 欧美黑人精品巨大| 国产成人精品久久二区二区91| 日韩一卡2卡3卡4卡2021年| 啦啦啦视频在线资源免费观看| 久久香蕉精品热| 丰满的人妻完整版| 日韩人妻精品一区2区三区| 欧美日韩亚洲国产一区二区在线观看 | a级毛片在线看网站| 91麻豆av在线| 中文字幕精品免费在线观看视频| 婷婷精品国产亚洲av在线 | 国精品久久久久久国模美| 亚洲中文字幕日韩| 国产欧美日韩一区二区三| 亚洲精品自拍成人| 国产成人啪精品午夜网站| 欧美激情久久久久久爽电影 | 中出人妻视频一区二区| 午夜福利影视在线免费观看| 国产精品久久久av美女十八| 久久中文字幕人妻熟女| 亚洲人成伊人成综合网2020| 日日夜夜操网爽| 12—13女人毛片做爰片一| 深夜精品福利| 桃红色精品国产亚洲av| 久久影院123| 亚洲熟女毛片儿| 在线永久观看黄色视频| 国产成人精品久久二区二区91| 中文字幕色久视频| 超碰成人久久| 亚洲第一av免费看| 满18在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 一区二区三区国产精品乱码| 天天躁夜夜躁狠狠躁躁| 成人免费观看视频高清| bbb黄色大片| 丝袜美足系列| 精品电影一区二区在线| 美女高潮到喷水免费观看| 香蕉国产在线看| 成人18禁高潮啪啪吃奶动态图| 亚洲中文字幕日韩| 高清视频免费观看一区二区| 欧美成狂野欧美在线观看| 亚洲欧美日韩高清在线视频| 悠悠久久av| 天堂√8在线中文| 女警被强在线播放| 国产1区2区3区精品| 99re在线观看精品视频| 国产视频一区二区在线看| 涩涩av久久男人的天堂| 亚洲欧美日韩另类电影网站| 日韩欧美三级三区| 久久久久国产一级毛片高清牌| 一级毛片女人18水好多| 欧美老熟妇乱子伦牲交| 色综合欧美亚洲国产小说| 天天躁夜夜躁狠狠躁躁| 亚洲精品自拍成人| 精品久久久精品久久久| 日韩 欧美 亚洲 中文字幕| 久久久久久亚洲精品国产蜜桃av| 999久久久国产精品视频| 桃红色精品国产亚洲av| 亚洲一区二区三区不卡视频| 高清视频免费观看一区二区| 大码成人一级视频| 午夜精品国产一区二区电影| cao死你这个sao货| 丁香六月欧美| 国产av一区二区精品久久| 日本一区二区免费在线视频| 欧美日韩精品网址| 久久国产精品大桥未久av| 久久ye,这里只有精品| 国产精品影院久久| 久99久视频精品免费| 老司机在亚洲福利影院| 精品人妻在线不人妻| 欧美日韩福利视频一区二区| 欧美人与性动交α欧美软件| 两个人免费观看高清视频| 免费观看a级毛片全部| 99久久99久久久精品蜜桃| 国产av一区二区精品久久| 欧美成人免费av一区二区三区 | 久久国产精品人妻蜜桃| 国产av一区二区精品久久| e午夜精品久久久久久久| 欧美精品亚洲一区二区| 欧美日韩亚洲综合一区二区三区_| 色尼玛亚洲综合影院| 久久这里只有精品19| 又黄又粗又硬又大视频| 欧美av亚洲av综合av国产av| 午夜福利在线观看吧| 50天的宝宝边吃奶边哭怎么回事| 韩国精品一区二区三区| 最近最新中文字幕大全免费视频| 欧美精品啪啪一区二区三区| 少妇被粗大的猛进出69影院| 欧美人与性动交α欧美精品济南到| 国产成人啪精品午夜网站| 脱女人内裤的视频| 老汉色∧v一级毛片| 在线观看免费视频网站a站| 桃红色精品国产亚洲av| 99国产精品99久久久久| 日韩成人在线观看一区二区三区| 久久精品91无色码中文字幕| 黄色怎么调成土黄色| 亚洲美女黄片视频| 国产黄色免费在线视频| 啦啦啦视频在线资源免费观看| 国产成人系列免费观看| 嫁个100分男人电影在线观看| 中文字幕高清在线视频| 久久精品熟女亚洲av麻豆精品| 18在线观看网站| 脱女人内裤的视频| 欧美精品高潮呻吟av久久| 在线观看免费日韩欧美大片| 亚洲午夜理论影院| 欧美日韩成人在线一区二区| 亚洲情色 制服丝袜| 新久久久久国产一级毛片| 麻豆av在线久日| 精品少妇久久久久久888优播| 一进一出抽搐gif免费好疼 | 日本精品一区二区三区蜜桃| 国产区一区二久久| 侵犯人妻中文字幕一二三四区| 久久久久久久午夜电影 | 欧美日韩亚洲国产一区二区在线观看 | 一级毛片高清免费大全| 亚洲精品久久成人aⅴ小说| 中文字幕人妻熟女乱码| 一边摸一边抽搐一进一出视频| 99在线人妻在线中文字幕 | 91成人精品电影| 大香蕉久久成人网| 无限看片的www在线观看| 99精品欧美一区二区三区四区| 成人国语在线视频| 国产一区二区三区视频了| 日本a在线网址| 亚洲在线自拍视频| 无遮挡黄片免费观看| 久久九九热精品免费| 久久精品亚洲熟妇少妇任你| 捣出白浆h1v1| 亚洲第一欧美日韩一区二区三区| 美女视频免费永久观看网站| 高清黄色对白视频在线免费看| 国产激情欧美一区二区| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 又黄又粗又硬又大视频| 亚洲熟妇中文字幕五十中出 | 精品午夜福利视频在线观看一区| 久久久久久久国产电影| 欧美成狂野欧美在线观看| 免费日韩欧美在线观看| 日韩三级视频一区二区三区| 免费av中文字幕在线| 国产深夜福利视频在线观看| 一级,二级,三级黄色视频| 精品久久久精品久久久| 99re6热这里在线精品视频| 国产亚洲欧美精品永久| 欧美+亚洲+日韩+国产| 黄片大片在线免费观看| 欧美精品啪啪一区二区三区| 女性生殖器流出的白浆| 国产aⅴ精品一区二区三区波| 欧美av亚洲av综合av国产av| 国产成人av教育| 亚洲精品国产区一区二| 啪啪无遮挡十八禁网站| 看黄色毛片网站| 99精品欧美一区二区三区四区| 国产亚洲欧美精品永久| 自线自在国产av| 成年人黄色毛片网站| 亚洲色图av天堂| 极品教师在线免费播放| 黄色a级毛片大全视频| 久久亚洲精品不卡| 久久久国产成人精品二区 | 久久香蕉精品热| cao死你这个sao货| 精品久久久精品久久久| 亚洲精品美女久久av网站| 国产人伦9x9x在线观看| 国产一区在线观看成人免费| 免费在线观看黄色视频的| 曰老女人黄片| 黑丝袜美女国产一区| 99精品在免费线老司机午夜| 久久精品国产清高在天天线| 婷婷丁香在线五月| 精品亚洲成a人片在线观看| 制服诱惑二区| 国产视频一区二区在线看| 亚洲一码二码三码区别大吗| 韩国av一区二区三区四区| 久久香蕉国产精品| 人人澡人人妻人| 他把我摸到了高潮在线观看| 大片电影免费在线观看免费| 午夜影院日韩av| 久久久精品国产亚洲av高清涩受| 久久婷婷成人综合色麻豆| 一级,二级,三级黄色视频| 757午夜福利合集在线观看| 99国产综合亚洲精品| 一级毛片高清免费大全| 老汉色av国产亚洲站长工具| 搡老乐熟女国产| 久久精品熟女亚洲av麻豆精品| 久久精品国产99精品国产亚洲性色 | 18禁观看日本| 成人黄色视频免费在线看| 黑人巨大精品欧美一区二区mp4| 免费人成视频x8x8入口观看| 欧美人与性动交α欧美软件| 黑丝袜美女国产一区| 成人精品一区二区免费| 欧美日本中文国产一区发布| xxx96com| 黄片小视频在线播放| 一二三四在线观看免费中文在| 国产伦人伦偷精品视频| 好看av亚洲va欧美ⅴa在| 露出奶头的视频| 中文字幕av电影在线播放| 久久久久精品人妻al黑| 夜夜躁狠狠躁天天躁| 18禁国产床啪视频网站| 国产人伦9x9x在线观看| 日韩欧美一区二区三区在线观看 | 在线十欧美十亚洲十日本专区| 国产有黄有色有爽视频| 搡老熟女国产l中国老女人| 久久久国产精品麻豆| 99国产精品99久久久久| 亚洲中文av在线| 在线观看免费高清a一片| 国产精品.久久久| 成人精品一区二区免费| videosex国产| 色综合婷婷激情| 欧美日韩中文字幕国产精品一区二区三区 | 国产99久久九九免费精品| 亚洲男人天堂网一区| 人人澡人人妻人| 91精品三级在线观看| 免费观看精品视频网站| 人人妻人人爽人人添夜夜欢视频| 久久久久久久久免费视频了| 91精品三级在线观看| 久久青草综合色| 国产三级黄色录像| 亚洲中文日韩欧美视频| 国产男女超爽视频在线观看| 丝袜美足系列| 99精国产麻豆久久婷婷| 无限看片的www在线观看| 男女之事视频高清在线观看| 老司机亚洲免费影院| 国产精品久久久久久人妻精品电影| 午夜激情av网站| 久久久久国内视频| 在线观看www视频免费| 成年人黄色毛片网站| 12—13女人毛片做爰片一| 久99久视频精品免费| 99久久人妻综合| 在线av久久热| 午夜福利欧美成人| 色尼玛亚洲综合影院| 欧美老熟妇乱子伦牲交| 人人澡人人妻人| 岛国毛片在线播放| 真人做人爱边吃奶动态| 国产精品自产拍在线观看55亚洲 | 90打野战视频偷拍视频| 亚洲av成人一区二区三| 精品国产乱码久久久久久男人| 欧美国产精品va在线观看不卡| av网站免费在线观看视频| 一边摸一边做爽爽视频免费| 交换朋友夫妻互换小说| 欧美乱色亚洲激情| 51午夜福利影视在线观看| 欧美日韩av久久| 大香蕉久久成人网| 老熟妇仑乱视频hdxx| 两人在一起打扑克的视频| 国产成人系列免费观看| 色尼玛亚洲综合影院|