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

    砰擊載荷下輕質(zhì)波紋夾芯夾層板動力響應(yīng)特性分析

    2016-05-03 01:24:24程遠勝
    船舶力學(xué) 2016年10期
    關(guān)鍵詞:芯層中心點輕質(zhì)

    田 媛,劉 均,汪 浩,程遠勝

    (華中科技大學(xué)船舶與海洋工程學(xué)院,湖北武漢430074)

    砰擊載荷下輕質(zhì)波紋夾芯夾層板動力響應(yīng)特性分析

    田 媛,劉 均,汪 浩,程遠勝

    (華中科技大學(xué)船舶與海洋工程學(xué)院,湖北武漢430074)

    文章對輕質(zhì)波紋夾芯夾層板(Light Weight Corrugated-Core Sandwich Plates,LWCCSP)在不同入水速度下(1-6 m/s)的流-固耦合非線性動力學(xué)行為進行了分析。建立了氣—液—固三相數(shù)值模型,通過顯式動力求解獲得了輕質(zhì)波紋夾芯夾層板砰擊壓力的分布特點及結(jié)構(gòu)變形規(guī)律,并與同等質(zhì)量的加筋板在流固砰擊下的非線性力學(xué)行為進行了對比,并研究了輕質(zhì)波紋夾芯夾層板主要設(shè)計參數(shù)對其砰擊響應(yīng)的影響。研究結(jié)果表明,輕質(zhì)波紋夾芯夾層板較同等質(zhì)量的加筋板表現(xiàn)出更好的抗砰擊性能;下面板厚度、芯層厚度的增加在一定范圍內(nèi)可以有效提高輕質(zhì)波紋夾芯夾層板的抗砰擊性能。

    輕質(zhì)波紋夾芯夾層板;入水砰擊;動態(tài)響應(yīng);LS-DYNA

    0 引 言

    砰擊是一種典型的流-固沖擊現(xiàn)象,其廣泛存在于自然界中,船舶航行于惡劣的海洋環(huán)境中時會受到來自流體環(huán)境的砰擊載荷作用。當(dāng)船體結(jié)構(gòu)受到過大的砰擊載荷作用時會造成結(jié)構(gòu)的塑性變形甚至失效[1]。因此,各國軍、民用船舶規(guī)范對砰擊載荷下的結(jié)構(gòu)設(shè)計均做出了專門要求[2],各國學(xué)者對結(jié)構(gòu)入水產(chǎn)生的砰擊問題亦開展了大量的理論、數(shù)值分析和實驗研究工作。Ochi和Motter[3]基于沖量理論,根據(jù)實船的測量數(shù)據(jù)提出了預(yù)測砰擊壓力的經(jīng)驗公式。隨著計算機技術(shù)的發(fā)展,有限元仿真技術(shù)開始應(yīng)用于分析砰擊壓力問題。Yang[5]等運用FLUENT軟件,模擬分析了二維楔形體的入水砰擊問題。駱寒冰等[6-7]針對二維和三維楔形體的入水砰擊進行了數(shù)值模擬和分析。在現(xiàn)有的文獻中,對于楔形體的砰擊入水問題研究較多,而對于平底結(jié)構(gòu)砰擊入水問題(如雙體船濕甲板、巡邏艇外底板)的關(guān)注較少。同時,大多數(shù)研究僅局限于傳統(tǒng)結(jié)構(gòu)形式,而對于新興金屬夾層結(jié)構(gòu)的入水砰擊響應(yīng)研究還不多見。金屬夾層結(jié)構(gòu)具有重量輕、比強高、易于加工等優(yōu)點[8],其在生產(chǎn)工藝、靜動態(tài)力學(xué)特性等方面與傳統(tǒng)型式結(jié)構(gòu)有著較大區(qū)別,對其進行“多功能化”設(shè)計后,金屬夾層結(jié)構(gòu)在聲學(xué)、疲勞、抗沖擊等方面具有傳統(tǒng)加筋板結(jié)構(gòu)無法比擬優(yōu)勢[9-10]。近年來金屬夾層結(jié)構(gòu)在船舶與海洋工程領(lǐng)域的應(yīng)用得到了廣泛的重視[11-12]。輕質(zhì)波紋夾芯夾層板是金屬夾層結(jié)構(gòu)的一種,其由兩塊金屬面板和V型折疊式夾芯通過激光焊接連接而成,同樣具有金屬夾層結(jié)構(gòu)的優(yōu)良性能。

    本文以輕質(zhì)波紋夾芯夾層板(LWCCSP)為研究對象,利用非線性有限元程序LS-DYNA對其在砰擊載荷下的流—固耦合砰擊動力學(xué)行為進行分析,并與同質(zhì)量下的傳統(tǒng)加筋板抗砰擊性能進行對比。在此基礎(chǔ)上分析下面板厚度、芯層厚度等結(jié)構(gòu)參數(shù)對其砰擊性能的影響,旨在為指導(dǎo)輕質(zhì)波紋夾芯夾層板在雙體船等實際工程上的設(shè)計與應(yīng)用提供參考。

    1 結(jié)構(gòu)幾何模型與特征參數(shù)

    1.1 幾何模型

    圖1為LWCCSP結(jié)構(gòu)示意圖。LWCCSP由上、下面板及波紋夾芯組成。輕質(zhì)波紋夾芯夾層板長為a,寬為b,上面板厚度與下面板厚度分別為tt和tb,夾芯厚度為tc,芯層高度H,夾芯層單元間距c。本文所研究的LWCCSP結(jié)構(gòu)在長度方向包含6個相同的胞元。

    圖1 LWCCSP結(jié)構(gòu)示意圖Fig.1 A light weight corrugated-core sandwich plate

    圖2為對比加筋板結(jié)構(gòu)示意圖。其中加筋為9#球扁鋼,均勻分布于船寬方向。加筋板面板厚度為6 mm,總質(zhì)量為56.4 kg。

    LWCCSP基本參數(shù)尺寸見表1。其長度、寬度以及質(zhì)量與對比對象加筋板相同。本文將計算該尺寸下LWCCSP與加筋板兩者在1-6 m/s入水速度下的砰擊響應(yīng)特性。

    圖2 加筋板結(jié)構(gòu)示意圖Fig.2 The parameters of a stiffened plate

    表1 LWCCSP參數(shù)尺寸Tab.1 Size of LWCCSP

    2 數(shù)值計算模型

    2.1 有限元模型

    在研究平底結(jié)構(gòu)體入水時通常不能忽略空氣層的影響,位于結(jié)構(gòu)底部和水表面之間的空氣在結(jié)構(gòu)的入水過程中起到了緩沖的作用,忽略空氣的存在在理論上甚至?xí)o出奇異的解。因此需要建立一個包含氣一水一固三種介質(zhì)耦合作用的模型才能夠更加準(zhǔn)確地反映出結(jié)構(gòu)入水的真實過程。

    圖3 結(jié)構(gòu)砰擊入水模型及有限元模型Fig.3 Finite element model of structure,air and water

    圖3為LWCCSP砰擊響應(yīng)特性計算用氣—液—固三相有限元模型。流體域包括兩部分:上部分為空氣,下部分為水域。結(jié)構(gòu)初始時刻位于空氣域中,其底面距靜水面100 mm。

    圖3右側(cè)為加筋板與LWCCSP的結(jié)構(gòu)有限元模型。在實際船舶中,板格四周邊界為對稱邊界條件,同時四周均為強構(gòu)件,因此可以簡化為四邊固支。而在有限元計算中,將板格四邊固支無法模擬砰擊運動,因此會在兩種結(jié)構(gòu)四周加有等高等質(zhì)量的剛性圍板,以模擬其固支邊界。此時板格只有垂直于水面的自由度,其余自由度被四周剛形體所限制,與實際板格的邊界條件是相符的。

    加筋板與輕質(zhì)波紋夾芯夾層板網(wǎng)格劃分如圖3所示。對于流體域,采用非等分線段劃分網(wǎng)格,在結(jié)構(gòu)關(guān)注點之外的網(wǎng)格劃分得比較稀疏,設(shè)置Spacing Ratio為5;而將結(jié)構(gòu)附近的流體網(wǎng)格劃分得較密一些,采用等分網(wǎng)格,網(wǎng)格單元大小與結(jié)構(gòu)單元大小相當(dāng),以保證計算精度。

    2.2 材料模型

    文中結(jié)構(gòu)(LWCCSP與加筋板)采用Lagrange殼單元(Shell163)模擬(9號球扁鋼腹板用殼單元shell163模擬,球頭用梁單元beam161模擬),使用MAT_PLASTIC_KINEMATIC材料模型。結(jié)構(gòu)材料為某船用鋼材,其材料參數(shù):彈性模量E=210 GPa,泊松比ν=0.3,密度ρ=7 830 kg/m3,靜態(tài)屈服極限σ0= 317.6 MPa,本文采用Cowper-Symonds本構(gòu)方程描述材料在砰擊載荷下的應(yīng)變率效應(yīng):

    表2 空氣的線性多項式Polynomial狀態(tài)方程中各參數(shù)值Tab.2 The polynomial equation of state of the air parameters

    水和空氣均采用*MAT_NULL材料模式??諝獾臓顟B(tài)方程通過Linear-polynomial定義,水的壓力和體積的關(guān)系通過Gruneisen狀態(tài)方程給出,具體參數(shù)見表2和表3[13]。流體采用Euler/ALE多物質(zhì)單元算法模擬。

    2.3 流—固耦合與邊界條件

    對于LWCCSP結(jié)構(gòu),設(shè)置下面板與空氣和水域耦合;對于加筋板,設(shè)置面板與空氣和水域耦合,定義關(guān)鍵字*CONSTRAINED_LAGRANGE_IN_SOLID來實現(xiàn)結(jié)構(gòu)與流體間的耦合。在模擬結(jié)構(gòu)砰擊入水時,流體域大小必須滿足一定條件以避免邊界對結(jié)構(gòu)動力響應(yīng)的影響。本文設(shè)置流體域長度、寬度及深度應(yīng)分別是結(jié)構(gòu)長度、寬度及入水深度的4倍。同時在邊界上使用無反射邊界條件。

    表3 水的線性多項式Gruneisen狀態(tài)方程中各參數(shù)值Tab.3 The Gruneisen equation of state of the water parameters

    3 瞬態(tài)動力結(jié)果分析

    3.1 結(jié)構(gòu)砰擊入水時間歷程

    結(jié)構(gòu)的砰擊入水是一個復(fù)雜的瞬態(tài)流固耦合過程。圖4為LWCCSP以v=5 m/s恒定速度向水中下落時不同時刻液面的變化情況。結(jié)構(gòu)入水時,由于水的沖擊作用,一些水質(zhì)點脫離水面與空氣混合,形成一種新的氣水混合物,稱為空氣墊。當(dāng)結(jié)構(gòu)下落至接近空氣與水的交界面,結(jié)構(gòu)和空氣墊共同對水域產(chǎn)生擠壓作用,使得水面發(fā)生變形,某些水質(zhì)點獲得較高的運動速度而脫離原來的水面,在結(jié)構(gòu)的邊界處積累到一定量時形成射流。隨著結(jié)構(gòu)下落深度的增加,液面分離射流越來越明顯。

    圖4 LWCCSP以5 m/s速度入水砰擊過程Fig.4 Process of water entry of LWCCSP at v=5 m/s

    3.2 砰擊壓力

    結(jié)構(gòu)在砰擊入水過程中,砰擊壓力最大處一般發(fā)生在結(jié)構(gòu)的耦合面中心區(qū)域,因此本文選取兩種結(jié)構(gòu)的耦合面中心點砰擊歷程曲線進行分析。圖5為加筋板中心與LWCCSP下面板中心處在5 m/s入水速度下砰擊壓力時間歷程曲線。入水砰擊過程中加筋板與LWCCSP的耦合面受到的砰擊壓力隨時間變化的規(guī)律是一致的,即:結(jié)構(gòu)以恒定速度從空氣垂直落入水面,在空氣中運動時,壓力載荷幾乎為零。當(dāng)結(jié)構(gòu)下降到與水面相距極近時,結(jié)構(gòu)受到的壓力載荷迅速上升,在極短的時間(約5 ms)內(nèi)壓力達到峰值,又在極短的時間(約5 ms)內(nèi)衰減下去。加筋板受到的最大砰擊壓力Pmax=0.95 MPa,對應(yīng)的時間T1=24.1 ms,波紋夾層板受到的最大砰擊壓力Pmax=0.62 MPa,對應(yīng)的時間T2=23.2 ms。

    Chuang[14]在一系列砰擊試驗中觀察到剛體平板的砰擊壓力分布是不均勻的,其峰值自中心向邊緣逐漸減小。對于加筋板與LWCCSP,耦合面砰擊壓力的分布也是不均勻的,如圖6、7所示。圖6為加筋板在速度v=5 m/s,時間t=T1(加筋板中心點所受砰擊壓力最大時)時結(jié)構(gòu)的砰擊壓力分布曲線。圖7為波紋夾層板在速度v-=5 m/s,時間t=T2(波紋夾層板中心點所受砰擊壓力最大時)時結(jié)構(gòu)的砰擊壓力分布曲線。

    從圖6、7可以總結(jié)出如下規(guī)律:

    (1)兩種結(jié)構(gòu)受到的最大砰擊壓力的分布特點為中部的砰擊壓力較高,靠近邊界處的壓力較低。

    (2)兩種結(jié)構(gòu)砰擊壓力分布從中心區(qū)域至邊界呈減小趨勢,壓力峰值并不位于結(jié)構(gòu)幾何中心點。在落水速度v=5 m/s情況下,對于加筋板,結(jié)構(gòu)所受砰擊壓力最大處坐標(biāo)為(0,0.2),而波紋夾層板所受砰擊壓力最大處坐標(biāo)為(0.1,0.2)。這可能是由于結(jié)構(gòu)加筋或者夾層板芯層與面板交接線的存在,使得結(jié)構(gòu)局部剛度發(fā)生變化,從而引起了壓力分布的變化。

    經(jīng)過計算,兩種結(jié)構(gòu)在不同的入水速度下受到最大砰擊壓力的位置有所不同,但都分布于一定的區(qū)域內(nèi)。在1-6 m/s入水速度下,加筋板受到最大砰擊壓力的位置都位于以加筋板中心點為圓心,半徑r=0.2 m的圓形區(qū)域內(nèi)。而對于輕質(zhì)波紋夾層板,受到最大砰擊壓力的位置都位于x=±0.1(-0.2≤y≤0.2)的兩條直線上。

    圖5 加筋板中心點與LWCCSP下面板中心點入水砰擊歷程曲線Fig.5 Pressure-time history of the center of the stiffened plate and LWCCSP

    圖6 加筋板砰擊壓力分布Fig.6 Impact pressure distribution of the stiffened plate

    圖7 輕質(zhì)波紋夾芯夾層板砰擊壓力分布Fig.7 Impact pressure distribution of LWCCSP

    加筋板與LWCCSP耦合面受到砰擊壓力峰值隨速度的變化規(guī)律見圖8。在1-6 m/s入水速度下,兩種結(jié)構(gòu)壓力峰值隨入水速度的增加呈近似線性增長。從圖8中可以看出,在1 m/s至4 m/s入水速度下,輕質(zhì)波紋夾芯夾層板下面板所受砰擊壓力峰值與加筋板相差無幾,且隨著速度的增大(5-6 m/ s),輕質(zhì)波紋夾芯夾層板壓力峰值隨速度的增加幅度減小,使得壓力峰值小于響應(yīng)加筋板的壓力峰值。而輕質(zhì)波紋夾芯夾層板與加筋板耦合面平均壓力相差卻很小,輕質(zhì)波紋夾芯夾層板平均壓力甚至稍大,這可能由于輕質(zhì)波紋夾芯夾層板胞元沿板長分布均勻且密集,從而使整個耦合面砰擊壓力分布較加筋板均勻。

    3.3 結(jié)構(gòu)變形

    圖8 結(jié)構(gòu)砰擊壓力峰值隨速度的變化Fig.8 Impact pressure-velocity curve of structures

    圖9為兩種結(jié)構(gòu)在入水速度為5 m/s,受到砰擊壓力最大時觸水面的位移云圖。圖10、11為兩種結(jié)構(gòu)在入水速度為5 m/s時從中心點沿X,Y方向的變形曲線。

    圖9 加筋板與LWCCSP下面板位移云圖Fig.9 Displacement contours of the stiffened plate and bottom plate of the LWCCSP

    從圖10、11可以總結(jié)出如下規(guī)律:

    (1)加筋板變形的峰值位于幾何的中心點,輕質(zhì)波紋夾層板變形的峰值則有兩處,分別位于y=0 m,x=±0.1 m這兩點。

    (2)加筋板從中心點至邊界變形呈減小趨勢。在x=±0.2 m加筋處,局部剛度較大,因此變形較小。

    圖10 加筋板變形分布Fig.10 Deformation distribution of the stiffened plate

    圖11 輕質(zhì)波紋夾芯夾層板變形分布Fig.11 Deformation distribution of LWCCSP

    (3)輕質(zhì)波紋夾層板的變形自中心點沿X方向呈波浪式減小。在芯層與下面板相連處,局部剛度較大,變形減小,出現(xiàn)局部極小值;夾層板自中心點沿Y方向變形逐漸減小。

    圖12為加筋板與LWCCSP在入水砰擊過程中產(chǎn)生的最大變形隨入水速度的變化。在入水速度為1-6 m/s情況下,兩種結(jié)構(gòu)最大變形隨入水速度呈近似線性增長。在相同的入水速度下,輕質(zhì)波紋夾芯夾層板下面板所產(chǎn)生的最大變形比加筋板所產(chǎn)生的最大變形要小,且兩種結(jié)構(gòu)中心點變形的差距隨著速度的增加而增大。當(dāng)入水速度增至6 m/s時,輕質(zhì)波紋夾芯夾層板中心點變形比加筋板中心點變形減小了將近44%,表現(xiàn)出了優(yōu)越的抗砰擊響應(yīng)性能。

    圖12 加筋板與LWCCSP最大變形對比Fig.12 Max deflection comparison of the stiffened plate and LWCCSP

    4 結(jié)構(gòu)參數(shù)分析

    4.1 下面板厚度對抗砰擊性能的影響

    改變下面板厚度,分別取tb=1.5、2.5、3.5、4.5和5.5 mm,其余尺寸不變,研究觸水面板厚度對LWCCSP抗砰擊性能的影響。計算結(jié)果如圖13、14所示。

    圖13 結(jié)構(gòu)砰擊壓力隨下面板厚度的變化Fig.13 Impact pressure with different bottom plate thickness of structure

    圖14 下面板中心點變形隨下面板厚度的變化Fig.14 Deformation with different bottom plate thickness of structure

    取輕質(zhì)波紋夾層板觸水面板中心點為研究對象,分析結(jié)構(gòu)的砰擊壓力與變形隨下面板厚度的變化規(guī)律。由圖13可以看出,隨著下面板厚度的增加,耦合面平均壓力與中心點壓力都有上升的趨勢。這主要由于下面板厚度的增加增大了結(jié)構(gòu)的剛度,使結(jié)構(gòu)變硬,從而所受砰擊壓力隨之增加。圖14為不同下面板厚度對應(yīng)的結(jié)構(gòu)變形量。下面板厚度從1.5 mm增加到2.5 mm時,結(jié)構(gòu)整體變形降低了14%左右。繼續(xù)增大下面板厚度,結(jié)構(gòu)整體變形沒有明顯的改善,反而有增大的趨勢,這主要是由于隨著下面板厚度的增加,在結(jié)構(gòu)剛度增大的同時,砰擊壓力也隨之增大。綜上可知,下面板厚度的增加使輕質(zhì)波紋夾芯夾層板的砰擊壓力增加,而變形有下降趨勢。當(dāng)下面板厚度增加到一定數(shù)值,變形不再減小。且下面板厚度過大會導(dǎo)致結(jié)構(gòu)重量增加,喪失了輕質(zhì)波紋夾芯夾層板結(jié)構(gòu)的優(yōu)良性能。因此存在最優(yōu)下面板厚度使輕質(zhì)波紋夾芯夾層板抗砰擊性能最佳。

    4.2 芯層厚度對抗砰擊性能的影響

    改變夾芯層面厚度tc,分別取tc=1.5 mm、2.5 mm、3.5 mm、4.5 mm和5.5 mm,其余尺寸不變,研究夾芯層厚度對LWCCSP抗砰擊性能的影響。計算結(jié)果如圖15、16所示。

    圖15 砰擊壓力峰值隨芯層厚度的變化Fig.15 Impact pressure with different core thickness of structure

    圖16 下面板中心點變形隨芯層厚度的變化Fig.16 Deformation with different core thickness of structure

    由圖15可知,當(dāng)夾層板芯層的厚度在1.5~5.5mm范圍內(nèi),芯層厚度的增加對砰擊壓力的影響不大。而對結(jié)構(gòu)變形的減小是有一定好處的。從圖16可以看出,當(dāng)芯層厚度從1.5 mm增加至5.5 mm時,夾層板的變形減小了近50%??梢缘贸鼋Y(jié)論:在一定范圍內(nèi),輕質(zhì)波紋夾芯夾層板所受的砰擊壓力對芯層厚度不敏感,但增加芯層厚度可有效減小結(jié)構(gòu)的變形。

    5 結(jié) 論

    本文對輕質(zhì)波紋夾芯夾層板的砰擊響應(yīng)性能進行了數(shù)值仿真分析研究,并將其與同質(zhì)量的加筋板進行了對比。探討了輕質(zhì)波紋夾芯夾層板下面板厚度、芯層厚度對輕質(zhì)波紋夾芯夾層板入水砰擊壓力和結(jié)構(gòu)響應(yīng)的影響,得出以下結(jié)論:

    (1)在1-6 m/s入水速度下,輕質(zhì)波紋夾芯夾層板與加筋板所受砰擊壓力與入水速度呈近似線性關(guān)系。入水速度較小時(1-4 m/s),輕質(zhì)波紋夾芯夾層板并沒有明顯的抗砰擊優(yōu)勢,然而隨著入水速度的增加(6m/s),中心點砰擊壓力輕質(zhì)波紋夾芯夾層板較加筋板低將近15%。兩種結(jié)構(gòu)觸水面平均壓力相差不大。在相同的入水速度下,兩種結(jié)構(gòu)在入水過程中產(chǎn)生的砰擊壓力峰值自中心點至邊界呈波浪式減小趨勢。兩種結(jié)構(gòu)在不同的入水速度下受到最大砰擊壓力的位置有所不同,但都分布于一定的區(qū)域內(nèi)。在1-6 m/s入水速度下,加筋板受到最大砰擊壓力的位置都位于以加筋板中心點為圓心,半徑r=0.2 m的圓形區(qū)域內(nèi)。而對于輕質(zhì)波紋夾層板,受到最大砰擊壓力的位置都位于x=±0.1(-0.2≤y≤0.2)的兩條直線上。

    (2)在相同的入水速度下,兩種結(jié)構(gòu)在入水過程中產(chǎn)生的最大變形自中心點至邊界呈減小趨勢,局部(加筋處或夾層板芯層與下面板連接處)由于剛度較大,變形明顯減?。惠p質(zhì)波紋夾芯夾層板變形較加筋板明顯降低,且隨著入水速度的增加,這種趨勢更加顯著。當(dāng)速度增加至6 m/s,輕質(zhì)波紋夾芯夾層板最大變形比加筋板減小約44%。從結(jié)構(gòu)變形分析可以看出,輕質(zhì)波紋夾芯夾層板抗砰擊性能明顯優(yōu)于傳統(tǒng)加筋板。

    (3)在一定范圍內(nèi),輕質(zhì)波紋夾芯夾層板下面板厚度及芯層厚度的增加增大了整個輕質(zhì)波紋夾芯夾層板的剛度,使得輕質(zhì)波紋夾芯夾層板所受砰擊壓力增加,但較大程度上減小了輕質(zhì)波紋夾芯夾層板的整體變形,且存在最佳下面板厚度和芯層厚度使其砰擊性能最佳。在一定范圍內(nèi),輕質(zhì)波紋夾芯夾層板所受的砰擊壓力對芯層厚度不敏感,但增加芯層厚度可有效減小結(jié)構(gòu)的變形。

    [1]Paik J K,Shin Y S.Structural damage and strength criteria for ship stiffened panels under impact pressure actions arising from sloshing,slamming and green water loading[J].Ships and Offshore Structures,2006,1(3):249-256.

    [2]American Bureau of Shipping.ABS-Guide for building and classing high speed naval craft[S].2007.

    [3]Ochi M K,Motter L E.Prediction of slamming characteristics and hull responses for ship design[J].Transaction of the Society of Naval Architects and Marine Engineers,1973,81:144-176.

    [4]Stavovy A B,Chuang S L.Analytical determination of slamming pressures for high speed vessels in waves[J].Journal of Ship Research,1976,20:190-198.

    [5]Yang S H,Lee H H,Park T H,et al.Experimental and numerical study on the water entry of symmetric wedges and a stern section of modern containership[C]//PRADS.Houston,USA,2007:518-526.

    [6]駱寒冰,張 智,王 輝.三維加筋板楔形體結(jié)構(gòu)水彈性砰擊的ALE數(shù)值模擬研究[J].船舶力學(xué),2013,17(1-2):84-91. Luo H B,Zhang Z,Wang H.Hydroelastic simulation of one 3D free-drop wedge with stiffened panels using ALE method [J].Journal of Ship Mechanics,2013,17(1-2):84-91.

    [7]駱寒冰,等.基于顯式有限元方法的二維楔形剛體入水砰擊載荷并行計算預(yù)報[J].船舶力學(xué),2012,16(8):907-914. Luo H B,et al.Parallel computing simulation of water entry of a 2D rigid wedge using an explicit finite element method [J].Journal of Ship Mechanics,2012,16(8):907-914.

    [8]Garescì F.Static and dynamic analysis of bonded sandwich plates[J].International Journal of Adhesion and Adhesives, 2012,33:7-14.

    [9]Zhu S Q,Chai G B.Impact of aluminum,CFRP laminates,fibre-metal laminates and sandwich panels[J].Conference Proceedings of the Society for Experimental Mechanics Series,2013,7:199-205.

    [10]Sun Y G,Gao L.Structural responses of all-composite improved-pyramidal truss sandwich cores[J].Materials and Design,2013,43:50-58.

    [11]Sohel K M A,Richard Liew J Y,Yan J B,et al.Behavior of steel-concrete-steel sandwich structures with lightweight cement composite and novel shear connectors[J].Composite Structures,2012,94(12):3500-3509.

    [12]Bella G Di,Calabrese L,Borsellino C.Mechanical characterisation of a glass/polyester sandwich structure for marine applications[J].Materials and Design,2012,42:486-494.

    [13]張曉波.船底結(jié)構(gòu)砰擊時的流固耦合數(shù)值模擬[D].大連:大連理工大學(xué),2007.

    [14]Chuang S L.Investigation of impact of rigid and elastic bodies with water[R].1970,NSRDC-TR-3248.

    Dynamic response of Light Weight Corrugated-Core Sandwich Plates subjected to slamming impact

    TIAN Yuan,LIU Jun,WANG Hao,CHENG Yuan-sheng
    (School of Naval Architecture and Ocean Engineering,Huazhong University of Science and Technology,Wuhan 430074,China)

    The fluid-solid interaction dynamic responses for Light Weight Corrugated-Core Sandwich Plates (LWCCSP)under different water-impact velocities(1-6 m/s)are investigated.The FEM computational model of multi-physics(air-water-solid)is built.The characteristics of slamming pressure and structure deformation are obtained based on the explicit dynamics analysis and are compared with those of stiffened plates with the same mass.Meanwhile the influences of main parameters of LWCCSP are researched.The results show that anti-slamming properties of LWCCSP are much better than those of stiffened plates with the same mass.Adding the sheet thickness of bottom plate and sandwich core can effectively improve the anti-slamming properties of LWCCSP to some extent.

    Light Weight Corrugated-Core Sandwich Plates;slamming loads;dynamic response; LS-DYNA

    TM344.1

    :A

    10.3969/j.issn.1007-7294.2016.10.010

    1007-7294(2016)10-1299-10

    2016-04-08

    國家自然科學(xué)基金(51079058);國家自然科學(xué)基金(51279065)

    田 媛(1988-),女,碩士研究生,E-mail:tianyuan@hust.edu.cn;劉 均(1981-),男,博士,講師;汪 浩(1984-),男,博士;程遠勝(1962-),男,教授,博士生導(dǎo)師。

    猜你喜歡
    芯層中心點輕質(zhì)
    怎一個“輕質(zhì)”了得
    一種耐高溫能力強的橡膠輸送帶
    輕質(zhì)高強堇青石多孔陶瓷的制備與表征
    關(guān)于重芳烴輕質(zhì)化與分離的若干思考
    Scratch 3.9更新了什么?
    電腦報(2020年12期)2020-06-30 19:56:42
    如何設(shè)置造型中心點?
    電腦報(2019年4期)2019-09-10 07:22:44
    空中爆炸載荷下梯度波紋夾層板抗爆性能仿真研究
    車用新型結(jié)構(gòu)的夾層構(gòu)件
    汽車文摘(2017年4期)2017-04-27 02:05:07
    漢字藝術(shù)結(jié)構(gòu)解析(二)中心點處筆畫應(yīng)緊奏
    輕質(zhì)材料彈射模型飛機
    国产精品1区2区在线观看.| aaaaa片日本免费| 亚洲男人的天堂狠狠| 国产视频内射| 精品久久久久久久久久久久久 | 在线看三级毛片| 精品午夜福利视频在线观看一区| 亚洲 欧美一区二区三区| 好看av亚洲va欧美ⅴa在| 精品不卡国产一区二区三区| 精品国产一区二区三区四区第35| 久久天堂一区二区三区四区| 国产av不卡久久| 日日干狠狠操夜夜爽| 很黄的视频免费| 伦理电影免费视频| 亚洲专区国产一区二区| 亚洲国产毛片av蜜桃av| 免费观看人在逋| 亚洲成av人片免费观看| 久久久久久久久免费视频了| 精品国产超薄肉色丝袜足j| 色在线成人网| 叶爱在线成人免费视频播放| aaaaa片日本免费| 一级毛片精品| 国内少妇人妻偷人精品xxx网站 | 久久天堂一区二区三区四区| tocl精华| 亚洲人成网站高清观看| 午夜两性在线视频| 男人舔女人下体高潮全视频| 日日干狠狠操夜夜爽| 99久久精品国产亚洲精品| 一级毛片女人18水好多| 777久久人妻少妇嫩草av网站| 男女午夜视频在线观看| 亚洲精品美女久久av网站| 1024视频免费在线观看| 成熟少妇高潮喷水视频| 熟妇人妻久久中文字幕3abv| 久久久久亚洲av毛片大全| 亚洲激情在线av| 满18在线观看网站| 成人亚洲精品av一区二区| 91九色精品人成在线观看| 免费在线观看影片大全网站| 一区福利在线观看| 最近最新中文字幕大全免费视频| 亚洲第一欧美日韩一区二区三区| 一级毛片精品| 亚洲国产欧美日韩在线播放| 久久国产精品男人的天堂亚洲| av有码第一页| 国产熟女午夜一区二区三区| 精品国内亚洲2022精品成人| 老熟妇仑乱视频hdxx| 99精品欧美一区二区三区四区| 国产精品 国内视频| 老司机福利观看| 成人永久免费在线观看视频| 美女大奶头视频| 欧美日韩亚洲综合一区二区三区_| 日本一区二区免费在线视频| 人人妻人人看人人澡| 色哟哟哟哟哟哟| 欧美黄色片欧美黄色片| 欧美成狂野欧美在线观看| 亚洲欧洲精品一区二区精品久久久| 色老头精品视频在线观看| 午夜福利18| 欧美性猛交╳xxx乱大交人| 欧美日韩一级在线毛片| 亚洲国产欧美网| 男人舔女人的私密视频| 成在线人永久免费视频| 少妇的丰满在线观看| 国产伦一二天堂av在线观看| 久久国产精品男人的天堂亚洲| 国产伦人伦偷精品视频| 色综合欧美亚洲国产小说| 欧美绝顶高潮抽搐喷水| 午夜福利欧美成人| 国产成人精品无人区| 可以免费在线观看a视频的电影网站| av电影中文网址| 精品久久久久久,| 免费观看人在逋| 狂野欧美激情性xxxx| 免费在线观看视频国产中文字幕亚洲| 中文字幕人成人乱码亚洲影| 人成视频在线观看免费观看| 欧美黄色淫秽网站| 又黄又爽又免费观看的视频| 国产主播在线观看一区二区| cao死你这个sao货| 欧美国产日韩亚洲一区| 色婷婷久久久亚洲欧美| 国产av一区二区精品久久| 欧美zozozo另类| 国产精品亚洲av一区麻豆| 欧美乱色亚洲激情| 欧美av亚洲av综合av国产av| 久久久久久久精品吃奶| 欧美中文综合在线视频| 在线观看免费视频日本深夜| 国产av不卡久久| 黄色片一级片一级黄色片| 亚洲 欧美一区二区三区| 18禁观看日本| 中文字幕另类日韩欧美亚洲嫩草| 老司机福利观看| 国产精品98久久久久久宅男小说| 国产成年人精品一区二区| 国产熟女午夜一区二区三区| 欧美乱码精品一区二区三区| 91成人精品电影| 欧美日韩乱码在线| 国产久久久一区二区三区| 黄色成人免费大全| 天堂动漫精品| 日韩一卡2卡3卡4卡2021年| 无限看片的www在线观看| 国产精品久久久av美女十八| 麻豆av在线久日| 国产人伦9x9x在线观看| 午夜视频精品福利| av在线天堂中文字幕| 特大巨黑吊av在线直播 | 亚洲色图 男人天堂 中文字幕| 啦啦啦观看免费观看视频高清| 麻豆久久精品国产亚洲av| 变态另类丝袜制服| 97超级碰碰碰精品色视频在线观看| 久久精品亚洲精品国产色婷小说| 午夜福利一区二区在线看| 一本久久中文字幕| 9191精品国产免费久久| 久久人人精品亚洲av| 18禁美女被吸乳视频| www国产在线视频色| 老司机在亚洲福利影院| 免费在线观看完整版高清| 国产精华一区二区三区| 久久久久久久久中文| 成人三级黄色视频| 18禁裸乳无遮挡免费网站照片 | 久久亚洲精品不卡| 国产色视频综合| 精品久久久久久久人妻蜜臀av| 国产av在哪里看| 国产高清有码在线观看视频 | 看片在线看免费视频| 亚洲精品在线观看二区| 久久久久久久精品吃奶| 久久青草综合色| 日韩欧美国产一区二区入口| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 搡老岳熟女国产| 亚洲精华国产精华精| 淫妇啪啪啪对白视频| 亚洲三区欧美一区| 在线天堂中文资源库| 久99久视频精品免费| 黄色毛片三级朝国网站| 在线视频色国产色| 午夜免费激情av| 国产激情欧美一区二区| 亚洲国产精品999在线| 欧美成人一区二区免费高清观看 | www.熟女人妻精品国产| 日韩中文字幕欧美一区二区| 久久久精品欧美日韩精品| 精品欧美一区二区三区在线| 欧美乱色亚洲激情| 亚洲成av片中文字幕在线观看| 高潮久久久久久久久久久不卡| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av| 亚洲精品一区av在线观看| av在线天堂中文字幕| √禁漫天堂资源中文www| 精品一区二区三区视频在线观看免费| 免费高清在线观看日韩| www.熟女人妻精品国产| 亚洲美女黄片视频| 久久久久国内视频| 美国免费a级毛片| 国产一卡二卡三卡精品| 最新在线观看一区二区三区| 精品国产亚洲在线| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产欧美日韩在线播放| 母亲3免费完整高清在线观看| 我的亚洲天堂| www日本黄色视频网| 国产精华一区二区三区| 曰老女人黄片| 色综合站精品国产| 九色国产91popny在线| 日本五十路高清| 免费看十八禁软件| 国产熟女xx| 香蕉av资源在线| 一进一出抽搐动态| 熟女电影av网| 69av精品久久久久久| 精品第一国产精品| 国产又色又爽无遮挡免费看| 人人妻人人澡人人看| 免费在线观看视频国产中文字幕亚洲| 欧美激情极品国产一区二区三区| 国产精品 欧美亚洲| 18禁观看日本| 桃红色精品国产亚洲av| 欧美精品啪啪一区二区三区| 美女高潮到喷水免费观看| 欧美日韩亚洲国产一区二区在线观看| 久久久久久九九精品二区国产 | 亚洲国产精品久久男人天堂| 午夜影院日韩av| 国产成人精品无人区| 一级毛片精品| 日韩欧美一区二区三区在线观看| 成年免费大片在线观看| 亚洲国产欧美日韩在线播放| 99riav亚洲国产免费| 亚洲专区字幕在线| 国产精品一区二区三区四区久久 | 男人舔女人的私密视频| 两性夫妻黄色片| 日韩国内少妇激情av| 国产精品二区激情视频| www.精华液| 精品电影一区二区在线| 免费看十八禁软件| 两性午夜刺激爽爽歪歪视频在线观看 | 99久久99久久久精品蜜桃| 欧美黑人精品巨大| 欧美丝袜亚洲另类 | 不卡av一区二区三区| 长腿黑丝高跟| 日韩欧美免费精品| 神马国产精品三级电影在线观看 | 亚洲成人久久性| 日韩三级视频一区二区三区| 久久人妻av系列| 欧美在线黄色| 日本撒尿小便嘘嘘汇集6| 成人特级黄色片久久久久久久| 性色av乱码一区二区三区2| 夜夜躁狠狠躁天天躁| 亚洲国产看品久久| 久久久久久亚洲精品国产蜜桃av| 久久香蕉精品热| 男女之事视频高清在线观看| 看片在线看免费视频| 国产91精品成人一区二区三区| 国产成人欧美在线观看| 97碰自拍视频| 麻豆成人午夜福利视频| 一边摸一边做爽爽视频免费| 满18在线观看网站| 久久久久免费精品人妻一区二区 | 国产精品 欧美亚洲| 国产三级在线视频| 国产亚洲av嫩草精品影院| 激情在线观看视频在线高清| 狠狠狠狠99中文字幕| 一级毛片高清免费大全| 亚洲精品美女久久av网站| 丰满人妻熟妇乱又伦精品不卡| 亚洲真实伦在线观看| 99在线视频只有这里精品首页| 亚洲av五月六月丁香网| 极品教师在线免费播放| 免费在线观看日本一区| 欧美黄色片欧美黄色片| 最近最新中文字幕大全电影3 | 国产97色在线日韩免费| 久久婷婷成人综合色麻豆| 免费搜索国产男女视频| 曰老女人黄片| 欧美日韩亚洲综合一区二区三区_| 99精品在免费线老司机午夜| 天堂动漫精品| 亚洲第一电影网av| 亚洲片人在线观看| www国产在线视频色| 波多野结衣高清无吗| 给我免费播放毛片高清在线观看| 在线观看免费视频日本深夜| 在线观看午夜福利视频| 欧美不卡视频在线免费观看 | 久久草成人影院| 日韩欧美国产一区二区入口| 日本一区二区免费在线视频| 日韩免费av在线播放| 黄片大片在线免费观看| 日韩中文字幕欧美一区二区| 黄片大片在线免费观看| 中文亚洲av片在线观看爽| 99久久无色码亚洲精品果冻| 亚洲国产精品合色在线| 日日干狠狠操夜夜爽| 国产亚洲精品一区二区www| 国产午夜精品久久久久久| 亚洲欧洲精品一区二区精品久久久| 国产真人三级小视频在线观看| 久久这里只有精品19| 999精品在线视频| 一本久久中文字幕| 久久精品91蜜桃| a级毛片在线看网站| 亚洲自偷自拍图片 自拍| 女性生殖器流出的白浆| 国产精品野战在线观看| videosex国产| 9191精品国产免费久久| av在线天堂中文字幕| 亚洲中文日韩欧美视频| 欧美日韩中文字幕国产精品一区二区三区| 99精品在免费线老司机午夜| 久久香蕉国产精品| 亚洲精品av麻豆狂野| 波多野结衣高清作品| 天堂影院成人在线观看| 黄片小视频在线播放| 亚洲国产精品久久男人天堂| 国产精品香港三级国产av潘金莲| 一卡2卡三卡四卡精品乱码亚洲| 99热6这里只有精品| 老司机靠b影院| 嫁个100分男人电影在线观看| 老熟妇乱子伦视频在线观看| 国产高清有码在线观看视频 | 动漫黄色视频在线观看| 久久午夜亚洲精品久久| 亚洲精品中文字幕在线视频| 别揉我奶头~嗯~啊~动态视频| 亚洲中文字幕一区二区三区有码在线看 | 丁香欧美五月| 色综合欧美亚洲国产小说| 一级黄色大片毛片| 国产欧美日韩精品亚洲av| 成人18禁在线播放| 18禁国产床啪视频网站| 国产aⅴ精品一区二区三区波| 国产主播在线观看一区二区| 在线播放国产精品三级| 国产男靠女视频免费网站| 午夜福利在线观看吧| 女性被躁到高潮视频| 搡老妇女老女人老熟妇| 国产免费男女视频| 极品教师在线免费播放| 欧美性长视频在线观看| 国产成人精品无人区| 日韩精品中文字幕看吧| 亚洲黑人精品在线| 国语自产精品视频在线第100页| 亚洲一码二码三码区别大吗| 日韩精品中文字幕看吧| 中亚洲国语对白在线视频| 国语自产精品视频在线第100页| 亚洲黑人精品在线| 国产亚洲av高清不卡| 淫秽高清视频在线观看| 国产精品亚洲美女久久久| 日韩欧美免费精品| 亚洲一区中文字幕在线| 1024视频免费在线观看| 国产午夜福利久久久久久| 麻豆av在线久日| 精品少妇一区二区三区视频日本电影| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 给我免费播放毛片高清在线观看| 一区二区三区高清视频在线| bbb黄色大片| 久久欧美精品欧美久久欧美| av片东京热男人的天堂| 亚洲全国av大片| 国产极品粉嫩免费观看在线| av天堂在线播放| 中文字幕精品亚洲无线码一区 | 19禁男女啪啪无遮挡网站| 免费在线观看日本一区| 90打野战视频偷拍视频| 国产成+人综合+亚洲专区| 国内毛片毛片毛片毛片毛片| 99国产精品一区二区三区| 国产成人精品久久二区二区91| 白带黄色成豆腐渣| 亚洲国产精品成人综合色| 又黄又爽又免费观看的视频| 亚洲精品久久国产高清桃花| 精品一区二区三区视频在线观看免费| 操出白浆在线播放| 午夜免费成人在线视频| 91麻豆av在线| 精品国产美女av久久久久小说| 最近最新中文字幕大全电影3 | 伦理电影免费视频| 十八禁网站免费在线| 日日爽夜夜爽网站| 香蕉丝袜av| 天天一区二区日本电影三级| 国产精品 国内视频| 久久香蕉精品热| 欧美成狂野欧美在线观看| 12—13女人毛片做爰片一| 国产精品美女特级片免费视频播放器 | 国产成人av激情在线播放| 亚洲国产欧洲综合997久久, | 黄频高清免费视频| 亚洲欧美精品综合一区二区三区| 国产色视频综合| 国产爱豆传媒在线观看 | 国产伦在线观看视频一区| 男女之事视频高清在线观看| 一二三四在线观看免费中文在| 久久精品成人免费网站| 身体一侧抽搐| 99久久精品国产亚洲精品| 97人妻精品一区二区三区麻豆 | 免费看a级黄色片| 亚洲欧美精品综合久久99| 日韩欧美在线二视频| 国产又色又爽无遮挡免费看| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看完整版高清| 日日爽夜夜爽网站| 日本免费a在线| 亚洲色图av天堂| 男人的好看免费观看在线视频 | 视频区欧美日本亚洲| 国产伦一二天堂av在线观看| av超薄肉色丝袜交足视频| 亚洲成人免费电影在线观看| 国产麻豆成人av免费视频| 亚洲成人精品中文字幕电影| 国产av一区在线观看免费| 十八禁网站免费在线| 色播在线永久视频| 成人午夜高清在线视频 | 亚洲人成伊人成综合网2020| 中文字幕人成人乱码亚洲影| 精品久久久久久久久久久久久 | 欧美大码av| 国产精品,欧美在线| 免费在线观看日本一区| 精品国产乱码久久久久久男人| 午夜激情av网站| 国产成人影院久久av| 午夜视频精品福利| 成年女人毛片免费观看观看9| 男人舔女人下体高潮全视频| 黄色a级毛片大全视频| 久久久久久国产a免费观看| 又黄又粗又硬又大视频| 国产成人精品久久二区二区91| 国产视频一区二区在线看| 麻豆成人午夜福利视频| 成人三级黄色视频| 欧美大码av| 成人免费观看视频高清| 日本熟妇午夜| 嫁个100分男人电影在线观看| 亚洲专区中文字幕在线| 久久久久国产一级毛片高清牌| 香蕉av资源在线| 国产精品野战在线观看| 成人午夜高清在线视频 | 久久久久久亚洲精品国产蜜桃av| 国产av一区二区精品久久| 久久青草综合色| 精品国产乱码久久久久久男人| 亚洲精品一卡2卡三卡4卡5卡| 久久久国产成人免费| 亚洲真实伦在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产黄a三级三级三级人| 熟妇人妻久久中文字幕3abv| 亚洲五月婷婷丁香| 神马国产精品三级电影在线观看 | 又紧又爽又黄一区二区| 久久久国产成人免费| 在线视频色国产色| 久久热在线av| 亚洲性夜色夜夜综合| 少妇熟女aⅴ在线视频| a在线观看视频网站| 此物有八面人人有两片| 欧美日韩亚洲国产一区二区在线观看| 国产欧美日韩精品亚洲av| 欧美乱色亚洲激情| 91国产中文字幕| 日韩精品青青久久久久久| 99热这里只有精品一区 | 91老司机精品| 99国产综合亚洲精品| 又紧又爽又黄一区二区| 精品一区二区三区四区五区乱码| 久久婷婷人人爽人人干人人爱| 成人一区二区视频在线观看| 亚洲七黄色美女视频| 999久久久精品免费观看国产| 中文字幕精品亚洲无线码一区 | 可以免费在线观看a视频的电影网站| 18禁国产床啪视频网站| 国产熟女xx| 99riav亚洲国产免费| 99国产极品粉嫩在线观看| 国产精品久久电影中文字幕| 啦啦啦 在线观看视频| 免费在线观看视频国产中文字幕亚洲| 给我免费播放毛片高清在线观看| 久久国产乱子伦精品免费另类| 天堂影院成人在线观看| 黑丝袜美女国产一区| av免费在线观看网站| 在线观看免费日韩欧美大片| 亚洲一区二区三区色噜噜| 一本一本综合久久| 亚洲 欧美 日韩 在线 免费| a在线观看视频网站| 亚洲 欧美 日韩 在线 免费| 国产真实乱freesex| 欧美成人性av电影在线观看| 不卡一级毛片| 少妇的丰满在线观看| 男女床上黄色一级片免费看| 亚洲va日本ⅴa欧美va伊人久久| www.www免费av| 亚洲三区欧美一区| 真人一进一出gif抽搐免费| 久久久久久人人人人人| 热99re8久久精品国产| av欧美777| 久久热在线av| 日本a在线网址| 久久久久久免费高清国产稀缺| 成年免费大片在线观看| 国产成年人精品一区二区| 亚洲va日本ⅴa欧美va伊人久久| 久久午夜亚洲精品久久| 精华霜和精华液先用哪个| 人人澡人人妻人| 亚洲欧美一区二区三区黑人| 一区二区三区高清视频在线| 日韩欧美在线二视频| a级毛片a级免费在线| 亚洲人成网站在线播放欧美日韩| 国内精品久久久久久久电影| 欧美午夜高清在线| 国产亚洲精品一区二区www| 欧美乱妇无乱码| www日本黄色视频网| 日日爽夜夜爽网站| 久久天躁狠狠躁夜夜2o2o| 此物有八面人人有两片| 91大片在线观看| 性色av乱码一区二区三区2| 老司机靠b影院| 亚洲国产中文字幕在线视频| 国产亚洲欧美在线一区二区| av免费在线观看网站| 欧美色欧美亚洲另类二区| 欧美午夜高清在线| 国产免费男女视频| 日韩有码中文字幕| 午夜免费观看网址| 亚洲一区二区三区不卡视频| 一本一本综合久久| 日本 av在线| 亚洲中文日韩欧美视频| 久久久久九九精品影院| 18禁黄网站禁片免费观看直播| 国产激情偷乱视频一区二区| 色精品久久人妻99蜜桃| 久久久国产欧美日韩av| 精品国内亚洲2022精品成人| 欧美黑人巨大hd| 男女午夜视频在线观看| 中文字幕久久专区| 最近最新免费中文字幕在线| 色在线成人网| 国产精品野战在线观看| 18禁裸乳无遮挡免费网站照片 | 免费在线观看亚洲国产| 深夜精品福利| 欧美日韩精品网址| 免费在线观看亚洲国产| 波多野结衣高清作品| 婷婷精品国产亚洲av| 成人国语在线视频| 99在线人妻在线中文字幕| 久久国产精品影院| 在线av久久热| 亚洲精品美女久久av网站| 久热爱精品视频在线9| 成年免费大片在线观看| 黄色视频不卡| 黑人巨大精品欧美一区二区mp4| 久久草成人影院| 久久久国产精品麻豆| 动漫黄色视频在线观看| www.精华液| 欧美日韩乱码在线| 亚洲 国产 在线| 人人澡人人妻人| 中亚洲国语对白在线视频| 亚洲国产毛片av蜜桃av|