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

    多次水下爆炸作用下鋼板動(dòng)態(tài)響應(yīng)數(shù)值模擬

    2020-01-10 01:55:24張斐張春輝張磊王志軍周春桂
    中國(guó)艦船研究 2019年6期
    關(guān)鍵詞:固支薄率撓度

    張斐,張春輝,張磊,王志軍,周春桂

    1 海軍研究院,北京100161

    2 中北大學(xué)機(jī)電工程學(xué)院,山西太原030051

    0 引 言

    隨著現(xiàn)代科技的高速發(fā)展,反艦武器已從單一作戰(zhàn)發(fā)展為集群作戰(zhàn),艦船舷側(cè)不可避免地會(huì)受到戰(zhàn)斗部的多次攻擊,其爆破能力會(huì)造成艦船不同程度的破損,因此艦船水下抗爆性能的研究受到了前所未有的重視[1-2]。艦船舷側(cè)受到反艦武器攻擊時(shí),其材料與結(jié)構(gòu)會(huì)發(fā)生應(yīng)變率效應(yīng)及塑性變形等非線性問題[3]。由于艦船是由大量板、梁、桿等拼裝而成的大型鋼結(jié)構(gòu)[4],直接研究艦船結(jié)構(gòu)與材料的抗爆性能具有一定困難,通常選取艦船舷側(cè)的典型結(jié)構(gòu)作為研究對(duì)象。非接觸爆炸載荷作用下艦船舷側(cè)動(dòng)態(tài)響應(yīng)的研究通常選取固支背空?qǐng)A鋼板為研究對(duì)象[5]。國(guó)外學(xué)者對(duì)水下非接觸爆炸載荷作用下鋼板的動(dòng)態(tài)塑性變形進(jìn)行了大量試驗(yàn)。美國(guó)于1990 年制定了美軍標(biāo)MIL-STD-2149A(SH)[6],確定了船用鋼母材及焊縫在快速加載條件下抗斷裂性能的評(píng)估方法。2010 年,韓國(guó)Park 等[7]設(shè)計(jì)了水下爆炸鼓脹試驗(yàn)裝置,測(cè)試了船用鋼及焊縫在多次水下爆炸加載作用下的塑性變形。2015 年,印度Kumar 等[8]研究了高強(qiáng)度低合金鋼在水下的抗沖擊性能。

    國(guó)內(nèi)于1997 年制定了空爆下測(cè)試鋼板塑性變形的試驗(yàn)標(biāo)準(zhǔn),但關(guān)于測(cè)試多次水下爆炸載荷作用下鋼板塑性變形的公開報(bào)道尚不多見。目前國(guó)內(nèi)學(xué)者的研究主要集中于一次水下爆炸載荷作用下鋼板的動(dòng)態(tài)響應(yīng)。吳成等[9]研究了水中爆炸產(chǎn)生的強(qiáng)沖擊載荷作用下固支方板的塑性動(dòng)力響應(yīng)過程,推導(dǎo)了固支方板在沖擊載荷作用下變形的最終撓度的解析解。諶勇等[10]分析了簡(jiǎn)支剛塑性圓板受水下爆炸載荷時(shí)的塑性動(dòng)力響應(yīng)。任鵬等[11]研究了水下近爆載荷作用下艦艇水下結(jié)構(gòu)的動(dòng)態(tài)變形及失效毀傷模式。牟金磊等[12]借助有限元軟件MSC.Dytran 研究了固支方板在水下爆炸載荷作用下的動(dòng)態(tài)響應(yīng)。通過上述研究發(fā)現(xiàn),國(guó)外開展的多次水下爆炸鼓脹試驗(yàn)大多為工程項(xiàng)目,國(guó)內(nèi)公開的關(guān)于多次水下爆炸載荷作用下鋼板動(dòng)態(tài)響應(yīng)的試驗(yàn)與仿真研究較少??紤]到艦船的實(shí)際服役環(huán)境以及來自反潛武器集群作戰(zhàn)的威脅,艦船結(jié)構(gòu)不可避免地會(huì)遭受多次水下爆炸沖擊,因此,研究多次水下爆炸載荷作用下鋼板的動(dòng)態(tài)響應(yīng),可為艦船水下抗爆與防護(hù)結(jié)構(gòu)設(shè)計(jì)提供一定的依據(jù),具有重要的意義。鑒此,本文首先利用有限元軟件Autodyn 模擬背空固支45#鋼板在多次水下爆炸載荷作用下的動(dòng)態(tài)響應(yīng);然后采用自行設(shè)計(jì)的水下爆炸鼓脹試驗(yàn)裝置開展背空固支45#鋼板的多次水下爆炸加載試驗(yàn),并對(duì)比分析試驗(yàn)與數(shù)值仿真結(jié)果。

    1 數(shù)值仿真計(jì)算

    1.1 計(jì)算工況

    水下爆炸數(shù)值仿真計(jì)算中炸藥為球形TNT 裝藥,藥量分別為10,20,30 和40 g,裝藥距目標(biāo)鋼板中心400 mm(即爆距為400 mm)。相同炸藥量下分別對(duì)直徑600 mm、厚度5 mm 的45#鋼板施加連續(xù)的3 次水下爆炸載荷。表1 給出了4 種藥量下的仿真計(jì)算工況,工況編號(hào)Cm-n 中,C 為“Case”的首字母,m 為藥量,n 為爆炸加載次數(shù)。

    表1 計(jì)算工況Table 1 Calculation condition

    1.2 有限元模型

    由于球形裝藥、圓形鋼板均具有對(duì)稱性,因此利用Autodyn 建立了如圖1 所示的二維1/2 數(shù)值仿真模型,主要包括球形裝藥TNT、背空鋼板、邊界擋板、水域、空氣。其中,邊界擋板與背空鋼板的材料均選擇45#鋼,其本構(gòu)選擇考慮應(yīng)變率效應(yīng)的J-C 本構(gòu)模型,材料參數(shù)如表2 所示。TNT 選用JWL 狀態(tài)方程描述;空氣選用理想氣體狀態(tài)方程描述;水選用Shock 狀態(tài)方程描述,強(qiáng)度模型為Hydro,其參數(shù)均取自Autodyn 材料庫(kù)。

    圖1 首次水下爆炸加載的有限元模型Fig.1 Finite element model of the first underwater explosion loading

    表2 45#鋼參數(shù)Table 2 Parameters of 45#steel plate

    數(shù)值仿真中涉及到固體、氣體、液體等多物質(zhì)耦合,因此選擇Autodyn 軟件中自帶的流固耦合算法。根據(jù)黃興中[13]的研究,沖擊波在邊界處的反射會(huì)影響數(shù)值仿真的精度,水域的尺寸至少應(yīng)為炸藥尺寸的50 倍以上。本研究中炸藥的最大直徑為36 mm,為了盡可能減小邊界條件對(duì)計(jì)算結(jié)果的影響,結(jié)合前期試算結(jié)果,最終水域大小選擇為4 000 mm×4 000 mm,同時(shí)在水域的邊界處添加流出邊界條件。水域采用中心網(wǎng)格加密的方法,其中中心網(wǎng)格大小為2 mm×2 mm,其余網(wǎng)格大小為5 mm×5 mm。由于研究對(duì)象為固支背空?qǐng)A鋼板,因此在直徑為600 mm 的45#鋼板背向添加邊界擋板,擋板與鋼板圍成的區(qū)域內(nèi)賦予空氣,同時(shí)在45#鋼板周向與擋板添加如圖2 所示的固支邊界條件。從圖2 中可以看出,鋼板的實(shí)際抗爆面是直徑為380 mm 的圓形區(qū)域。固支背空?qǐng)A鋼板與邊界擋板的網(wǎng)格尺寸均為1 mm×1 mm。

    圖2 鋼板與邊界擋板的邊界條件Fig.2 Boundary conditions of the steel plate and the boundary plate

    1.3 多次水下爆炸加載的數(shù)值方法

    Autodyn 模擬多次水下爆炸加載下鋼板動(dòng)態(tài)響應(yīng)的數(shù)值方法的核心是將上一次數(shù)值計(jì)算的結(jié)果作為下一次數(shù)值計(jì)算的初始條件。實(shí)際試驗(yàn)過程中,2 次連續(xù)的水下爆炸加載過程中有一定的間隔時(shí)間,第2 次施加水下爆炸載荷時(shí)水域基本恢復(fù)初始狀態(tài)。因此數(shù)值仿真時(shí)只將第1 次計(jì)算得到的鋼板結(jié)果輸出并保存,然后重新建立相同的水域并讀取保存的結(jié)果后開始下一次計(jì)算。圖3是爆距為400 mm、炸藥量為30 g 時(shí),5 mm 厚的45#鋼板第2 次水下爆炸加載時(shí)的有限元模型,其中在鋼板、邊界擋板上添加固支邊界條件、水域的邊界處設(shè)置流出邊界條件;材料的本構(gòu)模型及參數(shù)與1.2 節(jié)中相同。

    圖3 二次水下爆炸加載的有限元模型Fig.3 Finite element model of the second underwater explosion loading

    2 數(shù)值計(jì)算結(jié)果與分析

    2.1 載荷特性分析

    本研究中背空固支圓板實(shí)際受爆面的直徑為380 mm,爆距為400 mm。根據(jù)薛貴?。?4]的研究結(jié)果可知,背空固支圓板受均布載荷作用。峰值壓力是水下爆炸沖擊波的重要表征參數(shù),Cole 經(jīng)驗(yàn)公式可以計(jì)算TNT 球形裝藥在水中爆炸產(chǎn)生的沖擊波峰值壓力Pm。

    式中:Pm為沖擊波峰值壓力,MPa;W 為TNT 藥量,kg;H 為爆距,m;R 為TNT 藥包半徑,m。

    首先通過Autodyn 軟件數(shù)值計(jì)算多次水下爆炸過程中沖擊波的峰值壓力,然后與Cole 經(jīng)驗(yàn)公式計(jì)算值進(jìn)行比較,從而驗(yàn)證數(shù)值模擬方法的可行性。5 個(gè)考核點(diǎn)距離爆心分別為100,200,300,400 和500 mm。C30-1,C30-2 工況中數(shù)值模擬的5 個(gè)考核點(diǎn)處沖擊波壓力時(shí)間歷程曲線如圖4 所示。從圖4 中可以看出,在炸藥量與爆距均相同的情況下,通過Autodyn 軟件模擬的2 次水下爆炸加載過程中沖擊波壓力結(jié)果具有良好的一致性。C30-1 工況中數(shù)值模擬與經(jīng)驗(yàn)公式計(jì)算的水下沖擊波峰值壓力結(jié)果對(duì)比如表3 所示。從表3 可以看出,數(shù)值模擬結(jié)果與Cole 經(jīng)驗(yàn)公式計(jì)算結(jié)果的最小誤差僅為6.03%,最大誤差不超過19.57%。這表明本研究中的數(shù)值仿真模型及方法可以有效地模擬多次水下爆炸加載過程中爆炸流場(chǎng)的載荷特性。

    圖4 數(shù)值仿真中沖擊波壓力時(shí)程曲線Fig.4 Time history curves of shock wave pressure in numerical simulation

    表3 沖擊波峰值壓力計(jì)算結(jié)果Table 3 Results of the shock wave peak pressure

    2.2 鋼板變形模式

    對(duì)上述12 種工況進(jìn)行數(shù)值仿真計(jì)算,表4 是12 種工況中鋼板塑性變形的測(cè)試結(jié)果,厚度減薄率是指鋼板加載前后同一位置厚度的變化程度,其計(jì)算公式為

    式中:εi為第i 次爆炸后的厚度減薄率,%;t0為爆炸前原始的板厚,mm;ti為第i 次爆炸后的板厚,mm。

    表4 數(shù)值仿真中鋼板的變形結(jié)果Table 4 Deformation results of the steel plate in numerical simulation

    從表4 中可以看出,當(dāng)藥量相同時(shí),45#鋼板的中心撓度隨著水下爆炸加載次數(shù)的增加而逐漸增大,但增加量卻逐漸減少。Murr 等[15]的研究結(jié)果表明,鋼板在重復(fù)沖擊載荷作用下會(huì)出現(xiàn)剩余強(qiáng)度加強(qiáng)現(xiàn)象,導(dǎo)致鋼板內(nèi)流動(dòng)應(yīng)力增加。因此,鋼板在連續(xù)多次水下爆炸加載過程中彎曲變形逐漸減小,即鋼板的中心撓度增加量逐漸減小。進(jìn)一步研究發(fā)現(xiàn),鋼板在連續(xù)多次水下爆炸加載過程中,第2 次加載后鋼板中心撓度增加量約為首次加載后中心撓度的1/3。

    圖5 是C10-1,C20-1,C30-1,C40-1 工況下鋼板的變形形貌。從圖5 中可以看出:4 種工況下鋼板受沖擊區(qū)域均出現(xiàn)了類球冠塑性變形,且變形程度隨著裝藥量的增加而逐漸增大。結(jié)合圖5與表4 中C10-1,C20-1,C30-1,C40-1 工況下鋼板中心撓度的變化情況,可以得到作用在鋼板上的峰值壓力與中心撓度的關(guān)系,如圖6 所示。從圖6 中可以看出,二者呈線性關(guān)系,滿足公式:

    式中,D為鋼板中心撓度,mm。

    圖5 不同裝藥量下鋼板的變形形貌Fig.5 Deformation of the steel plate morphology under different charge

    圖6 峰值壓力與中心撓度的關(guān)系Fig.6 Relationship between peak pressure and deflection in the center of the plate

    下面以藥量為30 g、爆距為400 mm 的工況為例說明多次水下爆炸作用下鋼板撓度與厚度減薄率的變化規(guī)律。測(cè)量C30-1,C30-2,C30-3 工況下鋼板受沖擊區(qū)域半徑不同位置處的厚度值,代入公式(1)計(jì)算得到不同位置處的厚度減薄率,如圖7(a)所示。圖7(b)是對(duì)應(yīng)位置處的撓度,其中x 為測(cè)量位置處的半徑,R 為鋼板實(shí)際受沖擊部分的半徑。從圖7 可以看出,3 種工況下鋼板的撓度與厚度減薄率變化趨勢(shì)相似,中心位置的撓度與厚度減薄率均為最大,隨著離中心位置的距離的增加,其撓度與厚度減薄率呈現(xiàn)減小趨勢(shì)。在炸藥量一定的情況下,撓度增加量隨著加載次數(shù)的增加而逐漸減小,但厚度減薄率增加量則無明顯減小趨勢(shì)。由體積不變?cè)砜芍状渭虞d時(shí)鋼板中心位置以彎曲變形為主,隨著加載次數(shù)的增加,鋼板的雙向拉伸變形程度逐漸增大。

    圖7 多次加載下鋼板撓度與厚度減薄率的變化情況Fig.7 Variation of deflection and thickness reduction of steel sheets subjected to multiple underwater explosions

    通過對(duì)比C10-3 與C30-1 工況發(fā)現(xiàn),多次水下爆炸加載與單次水下爆炸加載作用下鋼板塑性變形并不相同。2 種工況下鋼板受沖擊區(qū)域半徑不同位置處的厚度減薄率與撓度對(duì)比分別如圖8(a)與圖8(b)所示。從圖8 中可以看出,相同的炸藥量下,1 次水下加載與均勻分為3 次連續(xù)加載下鋼板的變形程度明顯不同。單次加載后鋼板中心位置的撓度為18.528 mm,而連續(xù)3 次均勻加載后鋼板中心位置的撓度為14.776 mm,二者相差20.25%;隨著距中心位置距離的增加,二者差值逐漸減小,邊界處的撓度幾乎相同。與連續(xù)3 次均勻加載后鋼板的變形結(jié)果相比,單次加載后鋼板的厚度減薄率明顯較大。Kumar 等[8]的研究結(jié)果表明,相較于傳統(tǒng)載荷,鋼板在沖擊載荷的作用下延展性降低,出現(xiàn)明顯的硬化和強(qiáng)化。因此,連續(xù)重復(fù)水下爆炸載荷作用下鋼板的動(dòng)態(tài)塑性變形能力降低,表現(xiàn)為撓度與厚度減薄率均出現(xiàn)一定程度的降低。

    圖8 多次加載與單次加載下鋼板撓度與厚度減薄率的變化情況Fig.8 Variation of deflection and thickness reduction of steel sheets under multiple loading and single loading

    3 試驗(yàn)驗(yàn)證

    3.1 試驗(yàn)裝置及實(shí)施

    下面以藥量30 g、爆距400 mm 的工況為例驗(yàn)證數(shù)值仿真結(jié)果的準(zhǔn)確性。根據(jù)水下爆炸鼓脹試驗(yàn)特點(diǎn),設(shè)計(jì)了如圖9 所示的水下爆炸鼓脹試驗(yàn)裝置,主要由底座、密封圈、壓緊螺栓、壓鐵、固定桿等組成,其中,底座最大直徑600 mm、高300 mm,壓鐵面板厚30 mm。試驗(yàn)前將密封圈、試驗(yàn)板、壓鐵依次置于試驗(yàn)裝置底座上方,采用20 個(gè)高強(qiáng)螺栓將上述器材依次固定,保證試驗(yàn)過程中試驗(yàn)板處于背空固支狀態(tài)。從圖9 中可以看出固定后的試驗(yàn)板實(shí)際抗爆面為一直徑380 mm 的圓形區(qū)域。球形TNT 裝藥通過螺紋桿置于試驗(yàn)板中心的正上方,通過調(diào)節(jié)裝藥與試驗(yàn)板的距離(即爆距)來模擬不同爆炸載荷作用下鋼板的動(dòng)態(tài)響應(yīng)。水下爆炸試驗(yàn)水域?yàn)? m×9 m×5 m,試驗(yàn)時(shí)裝藥距水面3 m。

    將5 mm 厚(實(shí)測(cè)厚度為4.78 mm,鋼板平整度良好)的45#鋼板按試驗(yàn)要求置于水下爆炸鼓脹試驗(yàn)裝置,選取TNT 當(dāng)量為30 g,爆距為400 mm,對(duì)其施加3 次連續(xù)相同的水下爆炸載荷。距爆源1.24 m 處放置量程為35 MPa 的水壓傳感器來測(cè)量沖擊波在水中的傳遞情況。

    3.2 仿真結(jié)果與試驗(yàn)結(jié)果對(duì)比分析

    圖9 試驗(yàn)裝置示意圖及照片F(xiàn)ig.9 Schematic and photograph of the test device

    試驗(yàn)中通過水壓傳感器測(cè)得距爆源1.24 m 處的水壓時(shí)間歷程曲線以及數(shù)值仿真中C30-1 工況中同樣位置處的水壓時(shí)間歷程曲線如圖10 所示。通過經(jīng)驗(yàn)公式(1)計(jì)算得到距爆源1.24 m 處的峰值壓力為10.97 MPa;試驗(yàn)測(cè)得峰值壓力為11.34 MPa;數(shù)值仿真計(jì)算得到相同位置的峰值壓力為10.12 MPa,其與試驗(yàn)值相差10.76%,與經(jīng)驗(yàn)公式計(jì)算值相差7.75%。三者具有較好的一致性,因此本研究中的數(shù)值仿真模型及方法可以較為準(zhǔn)確地模擬多次水下爆炸試驗(yàn)中爆炸流場(chǎng)的載荷特性。

    圖10 試驗(yàn)與仿真的沖擊波壓力時(shí)程曲線對(duì)比Fig.10 Comparison of shock wave pressure time history curves between test and simulation

    試驗(yàn)與仿真中45#鋼背空?qǐng)A板經(jīng)3 次連續(xù)水下爆炸加載后的變形形貌分別如圖11 和圖12 所示,從圖中可以看出,試驗(yàn)與仿真中鋼板經(jīng)3 次水下爆炸加載后的塑性變形形貌呈類球冠形。試驗(yàn)中鋼板經(jīng)3 次水下爆炸加載后的中心撓度為27.2 mm,而仿真結(jié)果為25.28 mm,二者誤差僅為7.06%。試驗(yàn)與仿真中鋼板受沖擊區(qū)域半徑不同位置處的厚度減薄率與撓度對(duì)比分別如圖13(a)與圖13(b)所示。從圖13(a)中可以看出,與仿真結(jié)果相比,試驗(yàn)中鋼板邊界處的厚度減薄率偏大,而其他位置的厚度減薄率相差較小。從圖13(b)中可以發(fā)現(xiàn),試驗(yàn)中接近鋼板中心位置處的撓度較仿真結(jié)果偏大,遠(yuǎn)離中心位置處的試驗(yàn)與仿真結(jié)果一致性良好。這是由于仿真中邊界為固支狀態(tài),邊界處的變形模式主要為彎曲變形;而試驗(yàn)中鋼板邊界處雖有螺栓固定,但不可避免會(huì)存在一定的滑移現(xiàn)象,導(dǎo)致鋼板拉伸變形增加,因此試驗(yàn)中鋼板中心位置撓度與邊界處厚度減薄率相對(duì)偏大。整體來看,試驗(yàn)與仿真中鋼板的塑性變形結(jié)果具有較好的一致性,因此采用本研究中的數(shù)值仿真模型及方法可以較為準(zhǔn)確地模擬多次水下爆炸載荷作用下鋼板的動(dòng)態(tài)響應(yīng)。

    圖11 試驗(yàn)中3 次加載后鋼板變形形貌Fig.11 Deformation morphology of steel plate after three loading

    圖12 仿真中鋼板變形形貌Fig.12 Deformation morphology of steel plate in simulation

    圖13 試驗(yàn)與仿真中3 次加載后鋼板撓度與厚度減薄率對(duì)比Fig.13 Comparison of deflection and thickness reduction ratio of steel plate under multiple loading in test and simulation

    4 結(jié) 論

    本文采用數(shù)值仿真與試驗(yàn)相結(jié)合的方法研究了多次水下爆炸作用下5 mm 厚的45#鋼固支圓板的動(dòng)態(tài)響應(yīng),對(duì)比分析了試驗(yàn)與仿真結(jié)果,主要結(jié)論如下:

    1)多次水下爆炸載荷作用下,背空固支45#鋼板變形形貌呈類球冠形,變形模式以彎曲、雙向拉伸變形為主,且雙向拉伸變形程度隨著加載次數(shù)的增加而逐漸增大;鋼板塑性變形量隨著裝藥量的增加而增加,且中心撓度與作用在鋼板上的沖擊波峰值壓力成正比。

    2)多次水下爆炸載荷作用下,背空固支45#鋼板撓度與厚度減薄率從中心位置到邊界處呈減小趨勢(shì);鋼板中心撓度增量隨加載次數(shù)的增加而減小,第2 次加載后鋼板撓度增加量約為首次加載后撓度的1/3。

    3)同等藥量下,背空固支45#鋼板在1 次水下爆炸加載與均勻分為連續(xù)3 次水下爆炸加載作用下其塑性變形程度不同,與3 次加載后的結(jié)果相比,1 次爆炸加載下鋼板的撓度與厚度減薄率均偏大,中心撓度相差20.25%。

    4)所選用的數(shù)值仿真模型及方法可以較為準(zhǔn)確地模擬多次水下爆炸載荷作用下背空固支鋼板的動(dòng)態(tài)響應(yīng)。

    猜你喜歡
    固支薄率撓度
    固支結(jié)構(gòu)對(duì)壓阻式壓力傳感器的輸出影響研究
    基于有限元分析的某波紋管減薄率研究
    TA1板單點(diǎn)漸進(jìn)成形壁厚變化規(guī)律的數(shù)值模擬研究*
    高鋼級(jí)管道延性斷裂過程中壁厚減薄率研究*
    焊管(2021年10期)2021-11-04 03:37:24
    基于復(fù)合材料的山體滑坡固支模型及配套傳感器布置方法申請(qǐng)?zhí)? 202010514582.0
    傳感器世界(2021年2期)2021-03-27 19:36:46
    Spontaneous multivessel coronary artery spasm diagnosed with intravascular ultrasound imaging:A case report
    一角點(diǎn)支撐對(duì)面兩邊固支正交各向異性矩形薄板彎曲問題的辛疊加解
    約束分布對(duì)混合邊界薄板固有頻率影響的研究★
    山西建筑(2018年6期)2018-03-22 07:14:15
    天然氣高壓管道彎管壁厚減薄率的計(jì)算公式分析
    化工管理(2015年16期)2015-10-19 08:20:28
    懸高測(cè)量在橋梁撓度快速檢測(cè)中的應(yīng)用
    美女xxoo啪啪120秒动态图| 免费看不卡的av| 丝袜脚勾引网站| 色哟哟·www| 视频在线观看一区二区三区| 久久久欧美国产精品| 国产深夜福利视频在线观看| 免费黄色在线免费观看| 久久青草综合色| 日本猛色少妇xxxxx猛交久久| 欧美日韩精品成人综合77777| 亚洲精品乱码久久久久久按摩| 少妇被粗大的猛进出69影院 | 国产有黄有色有爽视频| 啦啦啦啦在线视频资源| 亚洲精品色激情综合| 国产探花极品一区二区| 少妇人妻精品综合一区二区| 国产免费视频播放在线视频| www.色视频.com| 日韩人妻高清精品专区| 你懂的网址亚洲精品在线观看| 麻豆成人av视频| 天天影视国产精品| 乱人伦中国视频| 99热这里只有精品一区| 久久午夜福利片| 草草在线视频免费看| 尾随美女入室| av黄色大香蕉| 亚洲欧美中文字幕日韩二区| 啦啦啦中文免费视频观看日本| 我的女老师完整版在线观看| 午夜日本视频在线| 久久女婷五月综合色啪小说| 99久久精品一区二区三区| 高清在线视频一区二区三区| 国产69精品久久久久777片| 中文乱码字字幕精品一区二区三区| 日本爱情动作片www.在线观看| 精品一区在线观看国产| 黑人高潮一二区| 国产极品粉嫩免费观看在线 | 少妇人妻 视频| 日本-黄色视频高清免费观看| 女性生殖器流出的白浆| 国产高清三级在线| 亚洲国产精品国产精品| 如何舔出高潮| 亚洲精品中文字幕在线视频| 国产精品久久久久久精品古装| 91精品伊人久久大香线蕉| 一本色道久久久久久精品综合| 亚洲色图 男人天堂 中文字幕 | 免费黄色在线免费观看| 99九九在线精品视频| 啦啦啦视频在线资源免费观看| 人人澡人人妻人| 青春草亚洲视频在线观看| 少妇的逼好多水| 亚洲国产av新网站| 亚洲美女视频黄频| 99热网站在线观看| 久久精品国产自在天天线| 菩萨蛮人人尽说江南好唐韦庄| 美女脱内裤让男人舔精品视频| 自线自在国产av| 热re99久久国产66热| 亚州av有码| 亚洲精品亚洲一区二区| 午夜免费观看性视频| av女优亚洲男人天堂| 简卡轻食公司| 亚洲成人手机| 免费观看无遮挡的男女| 久久97久久精品| 国产av国产精品国产| 亚洲,一卡二卡三卡| 欧美日韩视频高清一区二区三区二| 激情五月婷婷亚洲| 99热网站在线观看| 街头女战士在线观看网站| 免费高清在线观看日韩| 国产伦理片在线播放av一区| 乱码一卡2卡4卡精品| 亚洲国产欧美日韩在线播放| 中文乱码字字幕精品一区二区三区| 午夜福利视频精品| 99九九线精品视频在线观看视频| 国产精品三级大全| 777米奇影视久久| 国产精品国产三级国产av玫瑰| 最新中文字幕久久久久| 蜜臀久久99精品久久宅男| 永久网站在线| 国产免费福利视频在线观看| 国产精品一区二区在线观看99| 亚洲精品一二三| 免费看不卡的av| 人妻 亚洲 视频| 午夜视频国产福利| 久久av网站| 日韩成人伦理影院| 亚洲精品美女久久av网站| 一级毛片aaaaaa免费看小| 高清毛片免费看| 三上悠亚av全集在线观看| 久久精品夜色国产| 日本黄色片子视频| 国产乱来视频区| 黄色视频在线播放观看不卡| 国产精品一区www在线观看| 久久久久国产精品人妻一区二区| 三级国产精品片| 欧美亚洲日本最大视频资源| 狠狠精品人妻久久久久久综合| 中文字幕亚洲精品专区| 18禁动态无遮挡网站| 欧美成人精品欧美一级黄| 国产精品秋霞免费鲁丝片| 秋霞伦理黄片| 热re99久久精品国产66热6| 免费观看的影片在线观看| √禁漫天堂资源中文www| 免费高清在线观看视频在线观看| 久久国产精品大桥未久av| 蜜桃在线观看..| 午夜免费男女啪啪视频观看| 免费人成在线观看视频色| 久久精品熟女亚洲av麻豆精品| 男人添女人高潮全过程视频| 国产av一区二区精品久久| 另类亚洲欧美激情| 欧美+日韩+精品| 国产成人精品在线电影| 久久久久国产网址| 国产又色又爽无遮挡免| 99久久人妻综合| 乱人伦中国视频| 久久精品夜色国产| 欧美日韩av久久| 美女国产高潮福利片在线看| 好男人视频免费观看在线| 国模一区二区三区四区视频| av电影中文网址| 亚洲av免费高清在线观看| 亚洲精品日本国产第一区| 在现免费观看毛片| 亚洲国产色片| 亚洲四区av| 国产精品一区www在线观看| 精品人妻熟女av久视频| 成年人午夜在线观看视频| 精品国产国语对白av| 日韩欧美精品免费久久| 能在线免费看毛片的网站| 精品一区二区三区视频在线| 国产精品99久久99久久久不卡 | 99久久精品国产国产毛片| 国产av国产精品国产| 国产亚洲午夜精品一区二区久久| 我的女老师完整版在线观看| 中国三级夫妇交换| 中文精品一卡2卡3卡4更新| 久久午夜综合久久蜜桃| 久久人人爽av亚洲精品天堂| 天美传媒精品一区二区| 亚洲国产最新在线播放| 国产成人aa在线观看| 亚洲精品456在线播放app| 成年av动漫网址| 亚洲丝袜综合中文字幕| 亚洲在久久综合| 国产av精品麻豆| 男女边摸边吃奶| 国产精品.久久久| 精品人妻熟女毛片av久久网站| av在线老鸭窝| 少妇 在线观看| 一边亲一边摸免费视频| 国产欧美另类精品又又久久亚洲欧美| 日本av免费视频播放| 男女边摸边吃奶| 欧美精品一区二区大全| 最黄视频免费看| 精品久久蜜臀av无| 国产高清国产精品国产三级| 热re99久久国产66热| 下体分泌物呈黄色| 精品99又大又爽又粗少妇毛片| 夜夜爽夜夜爽视频| 99re6热这里在线精品视频| 在线天堂最新版资源| 国产免费福利视频在线观看| 中文天堂在线官网| 国产精品 国内视频| 亚洲欧美一区二区三区黑人 | 亚洲精品色激情综合| 五月玫瑰六月丁香| 日韩强制内射视频| 啦啦啦视频在线资源免费观看| 亚洲精品国产av成人精品| 最近2019中文字幕mv第一页| 欧美 日韩 精品 国产| 午夜免费观看性视频| 亚洲精品美女久久av网站| 国产有黄有色有爽视频| av免费观看日本| 爱豆传媒免费全集在线观看| 免费av中文字幕在线| 啦啦啦中文免费视频观看日本| 国语对白做爰xxxⅹ性视频网站| 在线 av 中文字幕| 最近的中文字幕免费完整| 熟妇人妻不卡中文字幕| 亚洲av成人精品一区久久| 99re6热这里在线精品视频| 免费观看av网站的网址| 国产 一区精品| 啦啦啦在线观看免费高清www| 免费观看性生交大片5| 99久久综合免费| 免费观看的影片在线观看| 韩国高清视频一区二区三区| 乱码一卡2卡4卡精品| 在线播放无遮挡| 久久人妻熟女aⅴ| 18在线观看网站| 亚洲伊人久久精品综合| av一本久久久久| 丰满少妇做爰视频| 美女视频免费永久观看网站| 久久久久久久亚洲中文字幕| 一个人看视频在线观看www免费| 一级毛片黄色毛片免费观看视频| 在线观看免费日韩欧美大片 | 最近的中文字幕免费完整| 插阴视频在线观看视频| 亚洲人与动物交配视频| 美女国产高潮福利片在线看| 美女国产高潮福利片在线看| 国产成人aa在线观看| 国产精品一区www在线观看| av女优亚洲男人天堂| 欧美xxxx性猛交bbbb| 久久久久精品性色| 熟妇人妻不卡中文字幕| 国产日韩欧美亚洲二区| 日韩在线高清观看一区二区三区| 国产av精品麻豆| 国产精品一国产av| 日韩av免费高清视频| 丰满少妇做爰视频| 女的被弄到高潮叫床怎么办| 日韩在线高清观看一区二区三区| 一级毛片黄色毛片免费观看视频| 久久免费观看电影| 街头女战士在线观看网站| 日韩中字成人| 自拍欧美九色日韩亚洲蝌蚪91| 久久ye,这里只有精品| av国产久精品久网站免费入址| 久久久久久久大尺度免费视频| 国产伦理片在线播放av一区| 天堂中文最新版在线下载| 一级,二级,三级黄色视频| 成人国产麻豆网| 激情五月婷婷亚洲| 99九九在线精品视频| 国产视频内射| 国产色爽女视频免费观看| 一个人看视频在线观看www免费| 天天操日日干夜夜撸| 草草在线视频免费看| 在线观看美女被高潮喷水网站| 久久久久精品性色| av在线app专区| 人妻制服诱惑在线中文字幕| 伊人亚洲综合成人网| 桃花免费在线播放| 久久99热这里只频精品6学生| 日韩亚洲欧美综合| 十八禁高潮呻吟视频| 蜜桃在线观看..| 天堂中文最新版在线下载| 日本vs欧美在线观看视频| 国产免费视频播放在线视频| 成人二区视频| 99九九在线精品视频| 亚洲精品中文字幕在线视频| 国产成人午夜福利电影在线观看| 在线看a的网站| 欧美日韩精品成人综合77777| 国产精品三级大全| 人妻人人澡人人爽人人| 亚洲国产成人一精品久久久| 18+在线观看网站| 91久久精品电影网| 亚洲av男天堂| 一级二级三级毛片免费看| 午夜激情av网站| 男人添女人高潮全过程视频| 日本-黄色视频高清免费观看| 女人精品久久久久毛片| 精品少妇内射三级| 欧美成人精品欧美一级黄| 91精品一卡2卡3卡4卡| 91午夜精品亚洲一区二区三区| 午夜激情久久久久久久| 99久国产av精品国产电影| 精品一区二区三区视频在线| 久久精品久久久久久久性| 国产深夜福利视频在线观看| 国产成人一区二区在线| .国产精品久久| 久久99蜜桃精品久久| 精品久久久久久电影网| 国国产精品蜜臀av免费| 免费观看性生交大片5| 国产在视频线精品| 香蕉精品网在线| 亚洲精品一二三| 男的添女的下面高潮视频| 精品久久蜜臀av无| 老司机亚洲免费影院| 欧美成人精品欧美一级黄| 日韩欧美一区视频在线观看| 三级国产精品欧美在线观看| 国产一区亚洲一区在线观看| 亚洲国产最新在线播放| 精品亚洲成国产av| 久久女婷五月综合色啪小说| 汤姆久久久久久久影院中文字幕| 各种免费的搞黄视频| 丝袜美足系列| 亚洲不卡免费看| 人妻制服诱惑在线中文字幕| 亚洲精品日韩av片在线观看| 777米奇影视久久| 日本猛色少妇xxxxx猛交久久| 国模一区二区三区四区视频| 亚洲国产最新在线播放| 国产欧美日韩一区二区三区在线 | 最近最新中文字幕免费大全7| 99久久综合免费| 久久精品国产亚洲av天美| 少妇高潮的动态图| 久久久久精品久久久久真实原创| av在线播放精品| 免费黄网站久久成人精品| 国产欧美亚洲国产| 黄片无遮挡物在线观看| 啦啦啦视频在线资源免费观看| 婷婷成人精品国产| 色5月婷婷丁香| 久久99一区二区三区| 日韩不卡一区二区三区视频在线| 大片免费播放器 马上看| 成年av动漫网址| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看| 欧美 亚洲 国产 日韩一| 青春草视频在线免费观看| 精品人妻偷拍中文字幕| 黄片播放在线免费| 久久久久国产精品人妻一区二区| 亚洲av在线观看美女高潮| videosex国产| 亚洲av男天堂| 亚洲国产最新在线播放| 老女人水多毛片| 免费观看av网站的网址| 欧美日韩精品成人综合77777| 国产黄片视频在线免费观看| 熟妇人妻不卡中文字幕| 99国产综合亚洲精品| 国产欧美另类精品又又久久亚洲欧美| 日韩亚洲欧美综合| 国产色爽女视频免费观看| 三级国产精品欧美在线观看| 欧美另类一区| 免费看av在线观看网站| 国产精品.久久久| 一级毛片 在线播放| 人妻少妇偷人精品九色| 新久久久久国产一级毛片| 国产精品一区二区在线观看99| 在线观看三级黄色| 国产日韩欧美视频二区| 少妇精品久久久久久久| 国产在线免费精品| 国产精品国产三级专区第一集| 又粗又硬又长又爽又黄的视频| 精品国产国语对白av| 18禁动态无遮挡网站| 色5月婷婷丁香| 男女边摸边吃奶| 2018国产大陆天天弄谢| 一本久久精品| 另类亚洲欧美激情| 国产精品一国产av| av卡一久久| 久久久久久久亚洲中文字幕| 91午夜精品亚洲一区二区三区| 综合色丁香网| 十八禁高潮呻吟视频| 91久久精品电影网| 制服人妻中文乱码| 精品99又大又爽又粗少妇毛片| 国产熟女欧美一区二区| 国产精品三级大全| 女性被躁到高潮视频| 一个人免费看片子| 亚洲激情五月婷婷啪啪| 男女高潮啪啪啪动态图| 人妻系列 视频| 最新中文字幕久久久久| 丰满迷人的少妇在线观看| 欧美日韩亚洲高清精品| 女性生殖器流出的白浆| 日韩大片免费观看网站| 熟女人妻精品中文字幕| 日日摸夜夜添夜夜添av毛片| 婷婷色综合大香蕉| 一本大道久久a久久精品| 国产午夜精品一二区理论片| 日韩中字成人| 国产精品久久久久成人av| 免费看不卡的av| 99热国产这里只有精品6| 国内精品宾馆在线| 人人妻人人爽人人添夜夜欢视频| 免费播放大片免费观看视频在线观看| av不卡在线播放| 精品人妻在线不人妻| 日韩精品有码人妻一区| 国产精品久久久久成人av| 在线播放无遮挡| 欧美bdsm另类| 国产亚洲精品久久久com| 国产毛片在线视频| 一本—道久久a久久精品蜜桃钙片| 最近2019中文字幕mv第一页| 看十八女毛片水多多多| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 久久人妻熟女aⅴ| 人人妻人人添人人爽欧美一区卜| 一个人免费看片子| av免费观看日本| 免费观看的影片在线观看| 久久久欧美国产精品| 人人妻人人澡人人看| 亚洲欧美中文字幕日韩二区| 亚洲精品美女久久av网站| 日韩大片免费观看网站| 欧美丝袜亚洲另类| 免费观看在线日韩| 在线观看美女被高潮喷水网站| videos熟女内射| 国产极品天堂在线| 午夜日本视频在线| 国产精品秋霞免费鲁丝片| 有码 亚洲区| 国产日韩欧美视频二区| 精品少妇黑人巨大在线播放| 国产欧美日韩综合在线一区二区| 丰满少妇做爰视频| 久久狼人影院| 亚洲美女黄色视频免费看| 啦啦啦视频在线资源免费观看| 老女人水多毛片| 大香蕉久久成人网| 最后的刺客免费高清国语| 日韩一区二区视频免费看| 久久精品久久久久久噜噜老黄| 简卡轻食公司| 五月开心婷婷网| 亚洲内射少妇av| 亚洲国产精品国产精品| 亚洲精品视频女| 性色avwww在线观看| kizo精华| 国产精品国产三级国产av玫瑰| 亚洲精品中文字幕在线视频| 午夜福利在线观看免费完整高清在| 久久久久久久精品精品| 亚洲欧美色中文字幕在线| 男的添女的下面高潮视频| 国产高清有码在线观看视频| 美女中出高潮动态图| 国产国语露脸激情在线看| 久久久久国产精品人妻一区二区| 亚洲成人av在线免费| 精品午夜福利在线看| 国产日韩欧美视频二区| 狠狠婷婷综合久久久久久88av| 黄片播放在线免费| 黄色欧美视频在线观看| 看非洲黑人一级黄片| 最近2019中文字幕mv第一页| 国产有黄有色有爽视频| av国产精品久久久久影院| 午夜视频国产福利| 免费不卡的大黄色大毛片视频在线观看| 成人二区视频| 亚洲精品乱码久久久久久按摩| 一本久久精品| 国产在线免费精品| 9色porny在线观看| 中文字幕久久专区| 亚洲欧洲国产日韩| 亚洲精品成人av观看孕妇| 夜夜看夜夜爽夜夜摸| av在线播放精品| 国产精品人妻久久久久久| 日韩大片免费观看网站| 国产成人aa在线观看| 波野结衣二区三区在线| 国产片内射在线| av在线app专区| 国产熟女午夜一区二区三区 | 蜜桃在线观看..| 爱豆传媒免费全集在线观看| 国产精品久久久久久久久免| 制服诱惑二区| 欧美亚洲日本最大视频资源| 欧美亚洲 丝袜 人妻 在线| 国产 精品1| 国产av精品麻豆| 一区二区三区四区激情视频| 亚洲欧美成人精品一区二区| 欧美精品国产亚洲| 日日啪夜夜爽| 亚洲天堂av无毛| 亚洲精品日本国产第一区| 女人久久www免费人成看片| 黄片播放在线免费| 精品一品国产午夜福利视频| 99久久中文字幕三级久久日本| 精品99又大又爽又粗少妇毛片| 日本与韩国留学比较| 国产精品一区www在线观看| 久久久久久久久久人人人人人人| 国产日韩一区二区三区精品不卡 | 国产成人精品在线电影| 免费观看性生交大片5| 伦精品一区二区三区| 久久久久人妻精品一区果冻| 国产色婷婷99| 精品久久蜜臀av无| 久久久国产欧美日韩av| 精品一区二区免费观看| 亚洲国产成人一精品久久久| 国产精品无大码| 天堂俺去俺来也www色官网| 亚洲激情五月婷婷啪啪| 九色亚洲精品在线播放| 亚洲精品亚洲一区二区| 婷婷色综合www| 看十八女毛片水多多多| 久久精品久久精品一区二区三区| 热99久久久久精品小说推荐| 99久久中文字幕三级久久日本| 自线自在国产av| av播播在线观看一区| 国产 精品1| 成人毛片60女人毛片免费| 一区二区三区免费毛片| 日本黄色日本黄色录像| 久久午夜综合久久蜜桃| 一级毛片aaaaaa免费看小| 日日撸夜夜添| 狠狠精品人妻久久久久久综合| 国产黄频视频在线观看| 久久人人爽人人片av| 汤姆久久久久久久影院中文字幕| 国产精品一区www在线观看| 永久免费av网站大全| 中文乱码字字幕精品一区二区三区| 国产熟女欧美一区二区| 人妻人人澡人人爽人人| 99国产综合亚洲精品| 性色avwww在线观看| 中文字幕人妻丝袜制服| 最近手机中文字幕大全| 免费观看在线日韩| 久久国内精品自在自线图片| 午夜福利在线观看免费完整高清在| 免费观看的影片在线观看| 欧美日韩一区二区视频在线观看视频在线| 午夜91福利影院| 狠狠婷婷综合久久久久久88av| 一本久久精品| 亚洲欧美精品自产自拍| 一本一本综合久久| 亚洲av日韩在线播放| 女性被躁到高潮视频| 亚洲精品aⅴ在线观看| 大香蕉97超碰在线| 亚洲国产精品一区二区三区在线| 亚洲av免费高清在线观看| 我的老师免费观看完整版| 国产精品久久久久久精品电影小说| 日韩人妻高清精品专区| 狂野欧美激情性bbbbbb| 亚洲精品日韩在线中文字幕| 人妻系列 视频| 制服诱惑二区| 三上悠亚av全集在线观看| 亚洲欧美成人综合另类久久久| 国产深夜福利视频在线观看| 少妇人妻久久综合中文| 免费黄频网站在线观看国产| 欧美精品国产亚洲|