• <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)壓縮數(shù)值模擬

    2023-10-31 13:48:16巫緒濤王寶珍程長(zhǎng)征馮學(xué)凱
    關(guān)鍵詞:數(shù)目力學(xué)性能孔隙

    郭 煜, 巫緒濤, 王寶珍, 程長(zhǎng)征, 馮學(xué)凱

    (合肥工業(yè)大學(xué) 土木與水利工程學(xué)院,安徽 合肥 230009)

    混凝土作為工程中常用的建筑材料,具有良好的力學(xué)性能,廣泛地應(yīng)用于建筑、地下防護(hù)、水利等領(lǐng)域。然而混凝土在澆筑、振蕩、脫模過(guò)程中會(huì)出現(xiàn)氣體未排干凈、漏振、灌注不暢通等現(xiàn)象,從而產(chǎn)生孔隙、裂紋等初始缺陷。深入研究孔隙對(duì)混凝土力學(xué)性能的影響,對(duì)提高其在工程應(yīng)用中的可靠性極為重要。

    在含孔隙混凝土的靜態(tài)力學(xué)性能研究方面,文獻(xiàn)[1]建立孔隙的二維介觀(guān)數(shù)值模型,將孔隙形狀、尺寸等多項(xiàng)特征進(jìn)行隨機(jī)化處理,研究孔隙對(duì)混凝土斷裂模式和承載能力的影響規(guī)律;文獻(xiàn)[2]在混凝土試樣制備過(guò)程中加入云母片來(lái)模擬缺陷,并進(jìn)行單軸壓縮試驗(yàn),探討不同尺寸和傾角的初始缺陷對(duì)混凝土變形和破壞的影響。上述研究均發(fā)現(xiàn),混凝土的初始缺陷對(duì)其靜態(tài)力學(xué)性能有顯著的不利影響。

    對(duì)于一些重要的民用和軍用防護(hù)設(shè)施,混凝土結(jié)構(gòu)在設(shè)計(jì)時(shí)還必須考慮承受爆炸、沖擊等強(qiáng)度較大的動(dòng)荷載。相關(guān)研究發(fā)現(xiàn),動(dòng)荷載作用下混凝土的力學(xué)性能和靜荷載下存在較大差異[3-4]。在含孔隙混凝土的動(dòng)態(tài)力學(xué)性能研究方面,文獻(xiàn)[5]建立不同孔隙率二維混凝土簡(jiǎn)支梁有限元模型,研究其在動(dòng)荷載作用下的彎拉破壞模式和宏觀(guān)力學(xué)性能;文獻(xiàn)[3]在混凝土加載方向上預(yù)制貫穿孔洞,進(jìn)行靜、動(dòng)態(tài)壓縮試驗(yàn),研究孔洞大小、位置及應(yīng)變率對(duì)混凝土抗壓強(qiáng)度的影響;文獻(xiàn)[6]采用數(shù)值方法,對(duì)砂漿、砂漿與骨料界面處的部分單元進(jìn)行缺陷化處理,建立含細(xì)觀(guān)缺陷的混凝土計(jì)算模型,研究試樣在3種沖擊荷載作用下應(yīng)力和應(yīng)變的分布規(guī)律。

    在混凝土試樣制備中,預(yù)制孔隙的同時(shí)也可能產(chǎn)生新的缺陷,且動(dòng)荷載作用下試樣內(nèi)部結(jié)構(gòu)的變形和破壞特征不易獲取。而在數(shù)值方法的研究中,由于建立的動(dòng)態(tài)有限元模型較復(fù)雜,且計(jì)算耗時(shí)較長(zhǎng),目前對(duì)混凝土缺陷的研究主要是針對(duì)貫穿型孔隙或在二維狀態(tài)下進(jìn)行數(shù)值模擬,對(duì)含內(nèi)部孔隙缺陷的三維混凝土結(jié)構(gòu)的研究較少。

    本文構(gòu)建含不同尺寸、位置、數(shù)目球形空洞的混凝土三維計(jì)算模型,模擬采用分離式霍普金森壓桿(split Hopkinson pressure bar,SHPB)進(jìn)行混凝土沖擊壓縮試驗(yàn)的過(guò)程,分析不同沖擊荷載作用下,多種孔隙特征參數(shù)對(duì)混凝土破壞形態(tài)及動(dòng)態(tài)力學(xué)性能的量化影響。

    1 有限元模型

    混凝土試樣模型的直徑為98 mm,厚度為50 mm。整體有限元模型中,入射桿、透射桿的直徑均為100 mm,長(zhǎng)度分別為3 000、2 000 mm。入射桿、透射桿、混凝土基體和孔隙均采用Solid164單元?jiǎng)澐帧?/p>

    近年來(lái),混凝土細(xì)觀(guān)骨料模型的建模方法已經(jīng)較成熟,其中背景網(wǎng)格法[7-10]具有易參數(shù)化、網(wǎng)格精度高的優(yōu)點(diǎn)。本文采用背景網(wǎng)格,在混凝土試樣模型中生成球形孔隙,算法步驟如下所述。

    1) 將試樣進(jìn)行背景網(wǎng)格劃分,并賦予混凝土的材料參數(shù)。

    2) 確定球形孔隙的形心坐標(biāo)(xi,yi,zi)和半徑R,滿(mǎn)足的條件為:

    (1)

    其中,Rs、Ls分別為試樣的半徑和厚度。試樣形心坐標(biāo)為(0, 0, 0)。

    3) 進(jìn)行干涉判斷,判斷第i個(gè)球形孔隙與前面生成的所有球形孔隙是否發(fā)生重疊:若不重疊,則生成球體;若發(fā)生重疊,則返回步驟2);對(duì)于單孔隙混凝土,該步驟省略。

    4) 獲取所有單元形心坐標(biāo)并存放數(shù)組,若某單元的形心位于球體投影范圍內(nèi),即單元形心與球體形心距離小于球體半徑,則將該材料參數(shù)更改為孔隙材料參數(shù);反之,仍為混凝土材料參數(shù)。

    用背景網(wǎng)格法建立的含球形孔隙混凝土試樣SHPB試驗(yàn)計(jì)算模型如圖1所示。

    圖1中:z方向?yàn)闂U和試樣的軸向;試樣徑向劃分單元數(shù)為100,軸向劃分單元數(shù)為50,單元數(shù)總計(jì)375 000;入射桿、透射桿徑向劃分單元數(shù)為10,軸向劃分單元數(shù)分別為300、200,單元數(shù)分別為57 600、38 400。

    本文重點(diǎn)研究孔隙大小、孔隙數(shù)目、孔隙位置、應(yīng)變率等因素對(duì)混凝土力學(xué)性能影響的一般規(guī)律,因此混凝土基體采用文獻(xiàn)[11]提出的HJC本構(gòu)原始參數(shù),見(jiàn)表1所列。表1中:ρ、G、fc、T分別為材料的密度、剪切模量、靜態(tài)抗壓強(qiáng)度、抗拉強(qiáng)度;A、B、C、N為強(qiáng)度參數(shù);D1、D2、εf min為損傷參數(shù);pc、μc分別為彈性狀態(tài)下的極限壓力及相應(yīng)的體積應(yīng)變,pl、μl分別為塑性狀態(tài)下的極限壓力及相應(yīng)的體積應(yīng)變;k1、k2、k3為壓力參數(shù);Smax為歸一化的最大強(qiáng)度。

    表1 混凝土基體的HJC本構(gòu)參數(shù)

    為了模擬出混凝土動(dòng)態(tài)破壞效果,需要加入失效準(zhǔn)則[4,12],且失效準(zhǔn)則采用主應(yīng)變失效才能得出與試驗(yàn)相近的破壞形態(tài),因此本文采用最大主應(yīng)變失效準(zhǔn)則。在A(yíng)NSYS/LS-DYNA有限元軟件中,可通過(guò)關(guān)鍵字“Mat-add-erosion”添加失效準(zhǔn)則,當(dāng)試樣單元的最大主應(yīng)變?chǔ)?大于所定義的失效主應(yīng)變?chǔ)舖ax時(shí),該單元就被刪除。對(duì)于失效主應(yīng)變的取值,文獻(xiàn)[13]指出,侵蝕發(fā)生在完全受損的材料中,因此本文設(shè)置最大主應(yīng)變失效閾值為0.02,確保所刪單元均為損傷單元。

    入射桿、透射桿材料采用線(xiàn)彈性模型,其密度為7 800 kg/m3,彈性模量為200 GPa,泊松比為0.3。球形孔隙采用線(xiàn)彈性模型[6],并通過(guò)單元失效使其在計(jì)算開(kāi)始時(shí)失效。桿與試樣端面采用面面自動(dòng)接觸,接觸剛度罰因子[14]設(shè)置為2.0。

    直接將速度時(shí)程入射波加載到入射桿自由端面,為減小入射波波形彌散效應(yīng),并模擬不同應(yīng)變率試驗(yàn),入射波波形選擇3種半正弦波,波長(zhǎng)歷時(shí)400 μs,幅值分別為4、6、8 m/s。

    2 數(shù)值模擬結(jié)果分析

    2.1 模型有效性驗(yàn)證

    背景網(wǎng)格法生成的球形孔隙為近似球體,用很小的彈性模量代替空材料;為驗(yàn)證其有效性,利用Truegrid軟件直接建立含球形空洞的混凝土模型。入射波、反射波和透射波作用下混凝土應(yīng)變分別為εi、εr、εt,相同狀態(tài)下2種方法得到的波形如圖2所示。由圖2可知,2種方法的計(jì)算結(jié)果基本一致。背景網(wǎng)格法的有效性由此得到驗(yàn)證,而其建模過(guò)程方便快捷,可大幅提高計(jì)算效率。

    圖2 2種方法計(jì)算波形對(duì)比

    由于所有計(jì)算工況均能在入射波上升沿持續(xù)時(shí)間內(nèi)達(dá)到應(yīng)力平衡,可根據(jù)模擬得到的反射波εr(t)、透射波εt(t)計(jì)算試樣的應(yīng)力和應(yīng)變時(shí)程,在共同時(shí)間坐標(biāo)下可得試樣的應(yīng)力-應(yīng)變關(guān)系為:

    (2)

    其中:A0、c0、E0分別為壓桿橫截面面積、波速和彈性模量;As為試樣橫截面面積。

    2.2 單孔隙混凝土數(shù)值模擬結(jié)果

    首先對(duì)含單孔隙混凝土的SHPB試驗(yàn)進(jìn)行模擬,試樣y=0直徑剖面下孔隙尺寸及位置示意如圖3所示。

    在不同入射波作用下,對(duì)試樣應(yīng)變率分別為28、45、54 s-1時(shí),4個(gè)孔隙尺寸和9個(gè)孔隙位置進(jìn)行計(jì)算,得到的混凝土試樣動(dòng)態(tài)抗壓強(qiáng)度見(jiàn)表2所列。表2中:N表示完整試樣;C表示試樣孔隙位于試樣中心;H表示孔隙球心關(guān)于試樣中心x方向左偏置,偏置量Hx分別為15、30 mm;V表示孔隙球心關(guān)于試樣中心z方向上下偏置,偏置量Vz分別為25、-25 mm(正值為靠入射桿側(cè));R1、R2、R3、R4分別表示孔隙半徑r為10、15、18、20 mm;φ為孔隙率;σ1、σ2、σ3分別為試樣在3個(gè)應(yīng)變率下的強(qiáng)度。

    表2 不同應(yīng)變率作用下混凝土動(dòng)態(tài)抗壓強(qiáng)度

    部分單孔隙試樣直徑剖面單元隨加載時(shí)間變化的應(yīng)力云圖如圖4所示。從圖4a~圖4c可以看出: 800 μs時(shí)C-R1試樣孔隙周邊出現(xiàn)劇烈的應(yīng)力集中,其中沿加載方向的孔周上下產(chǎn)生較大的拉應(yīng)力,垂直于加載方向的孔周左右產(chǎn)生較大的壓應(yīng)力;840 μs時(shí)在拉、壓應(yīng)力共同作用下,孔周出現(xiàn)明顯的“X”型剪切帶,試樣局部發(fā)生破壞;890 μs時(shí)端面和孔周之間形成貫通的“X”型斷裂路徑,試樣壓縮剪切失效[15-16]。從圖4d~圖4f可以看出, 890 μs時(shí),C-R3、V2-R2和H1-R2試樣圍繞孔隙均產(chǎn)生“X”型剪切帶和斷裂路徑,但影響范圍存在顯著差異。

    圖4 單孔隙混凝土SHPB試驗(yàn)的應(yīng)力云圖

    1) 孔隙尺寸對(duì)混凝土力學(xué)性能的影響。由表2可知:隨著應(yīng)變率增加,所有試樣的動(dòng)態(tài)抗壓強(qiáng)度顯著增大;隨孔隙尺寸增大,試樣的動(dòng)態(tài)抗壓強(qiáng)度減小。

    (3)

    圖5 單孔隙混凝土孔隙尺寸對(duì)動(dòng)態(tài)抗壓強(qiáng)度的影響擬合曲線(xiàn)

    2) 孔隙位置對(duì)混凝土力學(xué)性能的影響。根據(jù)表2數(shù)據(jù)可以得到:

    在應(yīng)變率和孔隙率不變的情況下,H系列試樣強(qiáng)度相較于C系列無(wú)明顯變化,內(nèi)部缺陷位置對(duì)試樣動(dòng)態(tài)抗壓強(qiáng)度的影響較小;當(dāng)缺陷向試樣端面靠近時(shí)(V1系列和V2系列),試樣強(qiáng)度稍微降低,即接近試樣端面的缺陷會(huì)導(dǎo)致該處提前破壞,試樣與壓桿接觸處的波阻抗降低,阻止應(yīng)力波透過(guò)試樣。

    2.3 多孔隙混凝土數(shù)值模擬

    對(duì)含隨機(jī)分布的多孔隙混凝土SHPB試驗(yàn)進(jìn)行模擬,在混凝土試樣內(nèi)部生成若干個(gè)互不干涉的球形孔隙。典型含多孔隙混凝土破壞時(shí)直徑剖面單元的應(yīng)力云圖如圖6所示。試樣從各孔隙附近單元開(kāi)始破壞并向外擴(kuò)展,各個(gè)孔隙的破壞路徑在多方向產(chǎn)生交匯,從而最終導(dǎo)致試樣碎成許多大塊,失去承載能力。

    圖6 典型多孔隙混凝土直徑剖面單元應(yīng)力云圖

    不同孔隙率φ、不同孔隙數(shù)目n下試樣在不同應(yīng)變率作用下的應(yīng)力-應(yīng)變曲線(xiàn)如圖7所示。圖7中,fdc′為含多個(gè)孔隙的混凝土動(dòng)態(tài)抗壓強(qiáng)度。從圖7可以看出,所有試樣的動(dòng)態(tài)抗壓強(qiáng)度都會(huì)隨著應(yīng)變率的增加而提高,這表明多孔隙試樣也具有應(yīng)變率效應(yīng)。

    圖7 不同孔隙率和孔隙數(shù)目下混凝土的應(yīng)力-應(yīng)變曲線(xiàn)

    需要指出,多孔隙混凝土孔隙率φ、孔隙數(shù)目n及孔隙半徑r之間的關(guān)系為:

    (4)

    其中,V為試樣體積。當(dāng)φ一定時(shí),n與r3成反比,

    φ=1.11%時(shí),n對(duì)強(qiáng)度影響不顯著,此時(shí)可按單孔隙確定試樣強(qiáng)度;φ=8.89%時(shí),隨著n增大,試樣強(qiáng)度呈先降低后收斂的趨勢(shì)。文獻(xiàn)[18]對(duì)含孔隙巖石進(jìn)行數(shù)值模擬發(fā)現(xiàn),孔隙率一定時(shí),隨著孔隙數(shù)目增多,試樣強(qiáng)度降低,這是由于模型中隨機(jī)分布的孔隙較多,產(chǎn)生的失效更有可能聚結(jié)形成持續(xù)的破壞路徑,導(dǎo)致構(gòu)件喪失承載能力。文獻(xiàn)[19]在對(duì)含10%泡孔的混凝土試樣進(jìn)行SHPB試驗(yàn)研究中發(fā)現(xiàn),當(dāng)泡孔體積分?jǐn)?shù)為10%時(shí),不同泡孔直徑的混凝土強(qiáng)度基本一致,也有強(qiáng)度收斂趨勢(shì)。

    為保證觀(guān)測(cè)結(jié)果的準(zhǔn)確性,加入φ為3.75%、6.48%時(shí),對(duì)不同孔隙數(shù)目的混凝土動(dòng)態(tài)抗壓強(qiáng)度進(jìn)行分析,4種孔隙率下混凝土在3種應(yīng)變率下的動(dòng)態(tài)抗壓強(qiáng)度隨孔隙數(shù)目變化如圖8所示。

    圖8 孔隙數(shù)目和應(yīng)變率對(duì)動(dòng)態(tài)抗壓強(qiáng)度的影響

    根據(jù)上述規(guī)律,混凝土動(dòng)態(tài)抗壓強(qiáng)度經(jīng)驗(yàn)方程可表述為:

    g2(φ)=exp(-b2φ),

    (5)

    其中:a2、b2、c2為無(wú)量綱參數(shù),由圖8中的數(shù)據(jù)點(diǎn)擬合得到a2=0.75,b2=0.03,c2=24.37。圖8中,擬合曲線(xiàn)綜合相關(guān)系數(shù)R2=0.984 3,表明擬合曲線(xiàn)與數(shù)值模擬結(jié)果吻合較好。

    根據(jù)以上研究可知,當(dāng)應(yīng)變率和孔隙率一定時(shí),隨著孔隙數(shù)目增加,混凝土動(dòng)態(tài)抗壓強(qiáng)度將趨向于一個(gè)穩(wěn)定值。

    為研究不同孔隙率和應(yīng)變率下該值的變化規(guī)律,取n=40,將各個(gè)孔隙率和應(yīng)變率下的強(qiáng)度代入式(3),擬合得到a1=1.11,b1=0.69,c1=0.03,擬合效果如圖9所示。

    圖9 孔隙率和應(yīng)變率對(duì)動(dòng)態(tài)抗壓強(qiáng)度的影響擬合曲線(xiàn)

    圖9中,擬合曲線(xiàn)綜合相關(guān)系數(shù)R2=0.988 0,說(shuō)明當(dāng)孔隙數(shù)目較多時(shí),式(3)也可以很好地預(yù)測(cè)混凝土動(dòng)態(tài)抗壓強(qiáng)度隨應(yīng)變率和孔隙率的變化規(guī)律。

    3 結(jié) 論

    1) 無(wú)論孔隙的尺寸、位置、數(shù)目如何,相較于完整試樣,孔隙的存在總會(huì)使試樣強(qiáng)度有所減小;隨著應(yīng)變率增大,同一試樣的強(qiáng)度均增大,含孔隙混凝土強(qiáng)度具有顯著的應(yīng)變率效應(yīng)。

    2) 單孔隙混凝土動(dòng)態(tài)抗壓強(qiáng)度隨著孔隙尺寸增大而減小,內(nèi)部孔隙位置的改變對(duì)強(qiáng)度的影響較小。本文得到的經(jīng)驗(yàn)公式(式(3))能較好反映應(yīng)變率、孔隙率對(duì)試樣強(qiáng)度的影響規(guī)律。

    3) 多孔隙混凝土孔隙率較低時(shí),孔隙數(shù)目對(duì)強(qiáng)度的影響較小,可以按同體積單孔隙試樣考慮;而孔隙率較大時(shí),強(qiáng)度隨著缺陷數(shù)目增加呈先減小后收斂的指數(shù)衰減趨勢(shì)。本文提出的經(jīng)驗(yàn)公式(式(5))可較好地描述這種變化規(guī)律,而混凝土動(dòng)態(tài)抗壓強(qiáng)度隨缺陷數(shù)目增加的最終收斂值可采用式(3)進(jìn)行估算。

    猜你喜歡
    數(shù)目力學(xué)性能孔隙
    有機(jī)物“同分異構(gòu)體”數(shù)目的判斷方法
    Pr對(duì)20MnSi力學(xué)性能的影響
    云南化工(2021年11期)2022-01-12 06:06:14
    二氧化碳在高嶺石孔隙中吸附的分子模擬
    Mn-Si對(duì)ZG1Cr11Ni2WMoV鋼力學(xué)性能的影響
    山東冶金(2019年3期)2019-07-10 00:54:00
    Preparation of bimodal grain size 7075 aviation aluminum alloys and the ir corrosion properties
    《哲對(duì)寧諾爾》方劑數(shù)目統(tǒng)計(jì)研究
    頁(yè)巖孔隙結(jié)構(gòu)研究進(jìn)展
    牧場(chǎng)里的馬
    INCONEL625+X65復(fù)合管的焊接組織與力學(xué)性能
    焊接(2015年9期)2015-07-18 11:03:53
    蘇北盆地金湖凹陷碳酸鹽巖孔隙類(lèi)型及孔隙結(jié)構(gòu)特征
    av不卡在线播放| 亚洲国产精品999| 欧美精品人与动牲交sv欧美| 亚洲全国av大片| www.熟女人妻精品国产| 精品福利永久在线观看| 欧美97在线视频| 五月天丁香电影| 国产国语露脸激情在线看| 亚洲精品中文字幕在线视频| 久久久欧美国产精品| 久久久久久人人人人人| 人人妻人人添人人爽欧美一区卜| 无遮挡黄片免费观看| 黄色怎么调成土黄色| 夜夜夜夜夜久久久久| 国产在线免费精品| 国产精品麻豆人妻色哟哟久久| 久久精品国产综合久久久| 女人精品久久久久毛片| 91av网站免费观看| 久久精品国产亚洲av香蕉五月 | 免费在线观看影片大全网站| 青青草视频在线视频观看| 久9热在线精品视频| 午夜福利在线免费观看网站| 午夜福利视频在线观看免费| 国产精品1区2区在线观看. | 免费人妻精品一区二区三区视频| 免费在线观看完整版高清| 亚洲三区欧美一区| 国产日韩一区二区三区精品不卡| 热99re8久久精品国产| 日本a在线网址| av视频免费观看在线观看| 老熟妇仑乱视频hdxx| 国产成人精品久久二区二区免费| 国产高清videossex| 欧美日韩一级在线毛片| 男女免费视频国产| 亚洲久久久国产精品| 在线亚洲精品国产二区图片欧美| 午夜激情av网站| 91精品国产国语对白视频| 狂野欧美激情性xxxx| 99香蕉大伊视频| 99香蕉大伊视频| 母亲3免费完整高清在线观看| 精品人妻熟女毛片av久久网站| 欧美97在线视频| 丝瓜视频免费看黄片| 一边摸一边做爽爽视频免费| 99九九在线精品视频| 中国国产av一级| 一本—道久久a久久精品蜜桃钙片| 免费看十八禁软件| 中文字幕另类日韩欧美亚洲嫩草| 日韩三级视频一区二区三区| 国产免费av片在线观看野外av| 女人高潮潮喷娇喘18禁视频| 国产一区二区 视频在线| 桃花免费在线播放| 欧美日韩成人在线一区二区| 国产精品一区二区在线观看99| 国产av一区二区精品久久| 亚洲专区中文字幕在线| 欧美性长视频在线观看| 亚洲欧美日韩高清在线视频 | 久久毛片免费看一区二区三区| 久久国产精品影院| 日日爽夜夜爽网站| 国产黄色免费在线视频| 久热这里只有精品99| 亚洲精品一卡2卡三卡4卡5卡 | 久久综合国产亚洲精品| 99久久国产精品久久久| 中亚洲国语对白在线视频| 大香蕉久久网| 精品视频人人做人人爽| 国产一级毛片在线| 不卡一级毛片| 一本久久精品| 久久综合国产亚洲精品| 欧美国产精品一级二级三级| 在线观看免费午夜福利视频| 永久免费av网站大全| 亚洲人成77777在线视频| 亚洲欧美一区二区三区黑人| 满18在线观看网站| 飞空精品影院首页| 国产成人欧美| 亚洲va日本ⅴa欧美va伊人久久 | 精品卡一卡二卡四卡免费| 天堂俺去俺来也www色官网| 伊人久久大香线蕉亚洲五| 日韩一卡2卡3卡4卡2021年| 亚洲熟女精品中文字幕| 99热国产这里只有精品6| 午夜激情久久久久久久| 国产1区2区3区精品| 手机成人av网站| 纯流量卡能插随身wifi吗| 日韩视频在线欧美| 精品少妇黑人巨大在线播放| 亚洲欧美精品自产自拍| 亚洲熟女精品中文字幕| 免费久久久久久久精品成人欧美视频| 亚洲色图综合在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 久久av网站| 欧美激情高清一区二区三区| 久久久国产成人免费| 高清在线国产一区| 搡老熟女国产l中国老女人| 成年av动漫网址| 成年av动漫网址| 精品福利永久在线观看| 亚洲第一av免费看| 三上悠亚av全集在线观看| 亚洲专区国产一区二区| 国产成+人综合+亚洲专区| 亚洲国产欧美网| 黄片播放在线免费| 丰满人妻熟妇乱又伦精品不卡| 亚洲av男天堂| 日韩一卡2卡3卡4卡2021年| 久久精品久久久久久噜噜老黄| 99热国产这里只有精品6| 嫁个100分男人电影在线观看| 成人亚洲精品一区在线观看| 亚洲精品美女久久久久99蜜臀| avwww免费| 国产xxxxx性猛交| 99久久人妻综合| 日韩欧美国产一区二区入口| 男人添女人高潮全过程视频| 久久中文看片网| 久久毛片免费看一区二区三区| av在线app专区| 亚洲全国av大片| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产欧美日韩在线播放| 性色av一级| 久久久久久人人人人人| 国产在线一区二区三区精| 三上悠亚av全集在线观看| 国产深夜福利视频在线观看| 成人国产一区最新在线观看| 极品少妇高潮喷水抽搐| 丝袜喷水一区| 精品视频人人做人人爽| 久久毛片免费看一区二区三区| 80岁老熟妇乱子伦牲交| 91麻豆精品激情在线观看国产 | 国产成人欧美| 亚洲专区中文字幕在线| 伊人久久大香线蕉亚洲五| 人人妻,人人澡人人爽秒播| 久久人人97超碰香蕉20202| 韩国精品一区二区三区| 97在线人人人人妻| 亚洲精品中文字幕在线视频| 欧美av亚洲av综合av国产av| videos熟女内射| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区有黄有色的免费视频| 亚洲国产看品久久| 精品一区二区三卡| 久久久久久人人人人人| 丝袜在线中文字幕| 999精品在线视频| 另类精品久久| 国产精品99久久99久久久不卡| 欧美性长视频在线观看| 一本久久精品| 国产在视频线精品| 热re99久久国产66热| 国产成人免费无遮挡视频| 国产av国产精品国产| 精品高清国产在线一区| 国产成人精品久久二区二区免费| 国产1区2区3区精品| 精品国产国语对白av| 国产老妇伦熟女老妇高清| 岛国在线观看网站| 无限看片的www在线观看| 人妻人人澡人人爽人人| 啦啦啦中文免费视频观看日本| 国产1区2区3区精品| 精品国产国语对白av| 一级毛片精品| 黑人欧美特级aaaaaa片| 人成视频在线观看免费观看| 日本a在线网址| a级毛片在线看网站| 亚洲欧美成人综合另类久久久| 国产成人欧美在线观看 | 9191精品国产免费久久| 91精品伊人久久大香线蕉| 国产亚洲午夜精品一区二区久久| 男女国产视频网站| 天堂8中文在线网| 黄色毛片三级朝国网站| 99九九在线精品视频| 欧美 亚洲 国产 日韩一| 黄色视频在线播放观看不卡| 亚洲三区欧美一区| 成年av动漫网址| 精品亚洲乱码少妇综合久久| 国产亚洲av片在线观看秒播厂| 十八禁高潮呻吟视频| 精品国产一区二区三区四区第35| 一区二区日韩欧美中文字幕| 久久亚洲精品不卡| 一区福利在线观看| 91精品三级在线观看| 国产三级黄色录像| 欧美激情高清一区二区三区| 十八禁人妻一区二区| 久久久国产成人免费| 久久久久视频综合| 国产精品免费大片| 久久人人爽av亚洲精品天堂| 精品少妇内射三级| 欧美性长视频在线观看| 97人妻天天添夜夜摸| 国产视频一区二区在线看| 老汉色∧v一级毛片| 婷婷成人精品国产| 国产在线一区二区三区精| 美女中出高潮动态图| 国产在线视频一区二区| 国产亚洲精品久久久久5区| 69精品国产乱码久久久| 中文精品一卡2卡3卡4更新| 日本猛色少妇xxxxx猛交久久| 搡老乐熟女国产| 我的亚洲天堂| 咕卡用的链子| 欧美乱码精品一区二区三区| 无限看片的www在线观看| 热99国产精品久久久久久7| 欧美国产精品va在线观看不卡| 精品福利观看| 老熟妇仑乱视频hdxx| 亚洲少妇的诱惑av| 午夜福利乱码中文字幕| 国产在线免费精品| 国产亚洲精品久久久久5区| 一级片免费观看大全| 国产精品偷伦视频观看了| 国产国语露脸激情在线看| 国产男人的电影天堂91| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩高清在线视频 | 午夜福利在线免费观看网站| 精品一区二区三区四区五区乱码| 在线观看免费日韩欧美大片| 淫妇啪啪啪对白视频 | 日韩制服骚丝袜av| 一进一出抽搐动态| 首页视频小说图片口味搜索| 精品第一国产精品| 无限看片的www在线观看| 国产亚洲av高清不卡| 乱人伦中国视频| 91成人精品电影| 黑人巨大精品欧美一区二区mp4| 久热爱精品视频在线9| 涩涩av久久男人的天堂| avwww免费| 国产精品麻豆人妻色哟哟久久| 精品视频人人做人人爽| 国产野战对白在线观看| 蜜桃在线观看..| 国内毛片毛片毛片毛片毛片| cao死你这个sao货| 日本一区二区免费在线视频| 黑人猛操日本美女一级片| 午夜91福利影院| 久久久久久免费高清国产稀缺| 久久国产精品男人的天堂亚洲| 在线天堂中文资源库| 久久精品人人爽人人爽视色| 欧美日韩成人在线一区二区| 色婷婷av一区二区三区视频| 国产高清国产精品国产三级| 无限看片的www在线观看| 老熟妇仑乱视频hdxx| 他把我摸到了高潮在线观看 | 女人高潮潮喷娇喘18禁视频| 欧美激情 高清一区二区三区| 啪啪无遮挡十八禁网站| 久久久精品国产亚洲av高清涩受| 一级毛片电影观看| 美女福利国产在线| 性色av一级| 精品人妻在线不人妻| 久久中文字幕一级| 99热国产这里只有精品6| 丰满迷人的少妇在线观看| 色老头精品视频在线观看| 国产深夜福利视频在线观看| 亚洲,欧美精品.| 亚洲av欧美aⅴ国产| 国产国语露脸激情在线看| 国产高清videossex| 亚洲专区字幕在线| 亚洲熟女毛片儿| 免费黄频网站在线观看国产| 一二三四社区在线视频社区8| 国产精品 国内视频| 极品人妻少妇av视频| xxxhd国产人妻xxx| 国产精品一区二区精品视频观看| 精品亚洲成a人片在线观看| 19禁男女啪啪无遮挡网站| 欧美日韩成人在线一区二区| 亚洲综合色网址| 超碰成人久久| 亚洲av成人不卡在线观看播放网 | 爱豆传媒免费全集在线观看| 美女福利国产在线| 亚洲精品在线美女| 美女高潮喷水抽搐中文字幕| 1024视频免费在线观看| 色综合欧美亚洲国产小说| 亚洲专区字幕在线| 国产精品影院久久| 国产精品.久久久| 国产色视频综合| 亚洲精品国产色婷婷电影| 不卡av一区二区三区| 久久久国产欧美日韩av| 最近中文字幕2019免费版| 久久久精品国产亚洲av高清涩受| 久久99热这里只频精品6学生| 性色av一级| 日韩欧美一区二区三区在线观看 | 国产一区二区在线观看av| 国产一区二区 视频在线| 欧美日韩精品网址| 日韩视频在线欧美| 永久免费av网站大全| 免费不卡黄色视频| 热re99久久国产66热| 亚洲全国av大片| 亚洲av日韩精品久久久久久密| 又大又爽又粗| 中文字幕最新亚洲高清| 水蜜桃什么品种好| 亚洲精品久久成人aⅴ小说| cao死你这个sao货| 狂野欧美激情性xxxx| 中文字幕精品免费在线观看视频| 十八禁高潮呻吟视频| 午夜福利免费观看在线| 日本撒尿小便嘘嘘汇集6| 大片免费播放器 马上看| 丰满饥渴人妻一区二区三| 亚洲第一青青草原| 人成视频在线观看免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 女性被躁到高潮视频| 在线观看一区二区三区激情| 亚洲精品中文字幕一二三四区 | 天天操日日干夜夜撸| 一级片'在线观看视频| 老司机亚洲免费影院| 精品亚洲成a人片在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩国产mv在线观看视频| 中文字幕制服av| av有码第一页| 人人妻人人澡人人看| 91九色精品人成在线观看| 亚洲 国产 在线| 黄色视频,在线免费观看| 精品久久蜜臀av无| 精品福利观看| 亚洲成人国产一区在线观看| 两性夫妻黄色片| 啦啦啦在线免费观看视频4| 午夜视频精品福利| 亚洲精品国产精品久久久不卡| 18禁黄网站禁片午夜丰满| 午夜日韩欧美国产| 日本五十路高清| 永久免费av网站大全| 国产亚洲午夜精品一区二区久久| 9热在线视频观看99| 国产主播在线观看一区二区| 视频在线观看一区二区三区| 正在播放国产对白刺激| 欧美亚洲日本最大视频资源| 性高湖久久久久久久久免费观看| 亚洲国产av影院在线观看| 美国免费a级毛片| 国产真人三级小视频在线观看| 久久免费观看电影| 亚洲国产精品一区三区| 老司机靠b影院| 亚洲av电影在线进入| 成年女人毛片免费观看观看9 | 极品少妇高潮喷水抽搐| 精品免费久久久久久久清纯 | 亚洲精品久久午夜乱码| 欧美亚洲日本最大视频资源| 亚洲精品美女久久av网站| av福利片在线| 精品少妇内射三级| 亚洲欧美色中文字幕在线| 欧美在线一区亚洲| 欧美另类一区| 美女视频免费永久观看网站| 久久精品人人爽人人爽视色| 爱豆传媒免费全集在线观看| 黄色视频在线播放观看不卡| 精品乱码久久久久久99久播| 最黄视频免费看| 国产成+人综合+亚洲专区| 人妻 亚洲 视频| tocl精华| 91精品伊人久久大香线蕉| 丝瓜视频免费看黄片| 欧美日韩亚洲综合一区二区三区_| 女人精品久久久久毛片| 亚洲一区中文字幕在线| kizo精华| 久久精品久久久久久噜噜老黄| 午夜福利影视在线免费观看| 国产淫语在线视频| 久久精品亚洲av国产电影网| 真人做人爱边吃奶动态| 91成人精品电影| 亚洲精华国产精华精| 我的亚洲天堂| 国产亚洲精品第一综合不卡| 最近最新中文字幕大全免费视频| 欧美中文综合在线视频| 999久久久精品免费观看国产| 一个人免费在线观看的高清视频 | 一级毛片精品| 久久国产精品影院| 国产真人三级小视频在线观看| 91国产中文字幕| 亚洲欧美清纯卡通| 精品乱码久久久久久99久播| 日本wwww免费看| 黄色视频在线播放观看不卡| 亚洲精品粉嫩美女一区| 欧美少妇被猛烈插入视频| 侵犯人妻中文字幕一二三四区| 成人免费观看视频高清| 777米奇影视久久| 欧美成人午夜精品| 亚洲九九香蕉| 欧美精品av麻豆av| av天堂在线播放| 老鸭窝网址在线观看| 日日爽夜夜爽网站| av欧美777| 久久久精品国产亚洲av高清涩受| 丝袜脚勾引网站| 日韩熟女老妇一区二区性免费视频| 一级片'在线观看视频| 少妇粗大呻吟视频| 亚洲九九香蕉| 亚洲欧美日韩高清在线视频 | 国产一级毛片在线| 亚洲九九香蕉| 久久综合国产亚洲精品| 午夜福利视频在线观看免费| 后天国语完整版免费观看| 亚洲专区字幕在线| 在线精品无人区一区二区三| 一本色道久久久久久精品综合| 热re99久久精品国产66热6| 黑丝袜美女国产一区| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产a三级三级三级| 亚洲精品久久久久久婷婷小说| 中文字幕高清在线视频| 久久国产精品男人的天堂亚洲| 亚洲精品自拍成人| 黄色怎么调成土黄色| 日韩精品免费视频一区二区三区| 欧美精品人与动牲交sv欧美| 国产成人啪精品午夜网站| 自线自在国产av| 免费高清在线观看日韩| 亚洲国产av影院在线观看| 少妇 在线观看| 香蕉丝袜av| 国产av国产精品国产| 考比视频在线观看| 色婷婷av一区二区三区视频| 亚洲熟女毛片儿| 51午夜福利影视在线观看| 久久人人爽av亚洲精品天堂| 美女午夜性视频免费| 国产极品粉嫩免费观看在线| 少妇 在线观看| 51午夜福利影视在线观看| 国产成人免费无遮挡视频| 久久久精品免费免费高清| 一进一出抽搐动态| 建设人人有责人人尽责人人享有的| 日本a在线网址| 亚洲国产欧美网| 99国产精品免费福利视频| 黑人猛操日本美女一级片| 天天躁狠狠躁夜夜躁狠狠躁| 国产国语露脸激情在线看| 丰满饥渴人妻一区二区三| 纵有疾风起免费观看全集完整版| 精品国产国语对白av| 精品人妻熟女毛片av久久网站| 伦理电影免费视频| 搡老岳熟女国产| 国产精品秋霞免费鲁丝片| 国产视频一区二区在线看| 交换朋友夫妻互换小说| 91精品国产国语对白视频| 亚洲精品一二三| 成年女人毛片免费观看观看9 | 日韩熟女老妇一区二区性免费视频| 成人18禁高潮啪啪吃奶动态图| 亚洲第一av免费看| 午夜老司机福利片| 热re99久久精品国产66热6| 肉色欧美久久久久久久蜜桃| 国产男女内射视频| 国产精品国产三级国产专区5o| 久久狼人影院| 黄色毛片三级朝国网站| 久久 成人 亚洲| 人人妻人人添人人爽欧美一区卜| av一本久久久久| 国产野战对白在线观看| 日本一区二区免费在线视频| 激情视频va一区二区三区| a级毛片在线看网站| 色老头精品视频在线观看| 欧美日韩亚洲综合一区二区三区_| 黄色毛片三级朝国网站| 久久av网站| 久久影院123| 考比视频在线观看| 国产精品免费大片| 巨乳人妻的诱惑在线观看| 一本色道久久久久久精品综合| 亚洲精品一区蜜桃| 热99re8久久精品国产| 久久久久久人人人人人| 麻豆国产av国片精品| 在线观看www视频免费| 这个男人来自地球电影免费观看| 午夜激情av网站| 99国产精品一区二区蜜桃av | 中文字幕另类日韩欧美亚洲嫩草| 国产成人免费无遮挡视频| 亚洲精品美女久久久久99蜜臀| 久久热在线av| 久久亚洲国产成人精品v| 久久午夜综合久久蜜桃| 欧美黑人精品巨大| 欧美另类一区| 日韩欧美一区二区三区在线观看 | 水蜜桃什么品种好| 天天操日日干夜夜撸| 曰老女人黄片| 啪啪无遮挡十八禁网站| 曰老女人黄片| 国产精品久久久久久人妻精品电影 | 欧美日韩视频精品一区| 国产精品一二三区在线看| 久久性视频一级片| 日韩一区二区三区影片| 免费高清在线观看视频在线观看| 手机成人av网站| 一个人免费在线观看的高清视频 | 国产精品欧美亚洲77777| 午夜福利视频在线观看免费| 成人国语在线视频| 亚洲欧美色中文字幕在线| 老司机午夜十八禁免费视频| 丝袜喷水一区| 午夜老司机福利片| 巨乳人妻的诱惑在线观看| 女人精品久久久久毛片| 一个人免费看片子| 亚洲国产欧美在线一区| 国产又爽黄色视频| 91麻豆av在线| 80岁老熟妇乱子伦牲交| 久久午夜综合久久蜜桃| 久久av网站| av国产精品久久久久影院| 国产精品九九99| 看免费av毛片| 视频区欧美日本亚洲| 亚洲av成人一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 在线十欧美十亚洲十日本专区| 高清黄色对白视频在线免费看| 国产老妇伦熟女老妇高清| 日韩人妻精品一区2区三区| 日韩中文字幕欧美一区二区| 国产色视频综合| 999久久久国产精品视频| 99热网站在线观看| 一个人免费看片子| 午夜久久久在线观看| 日本91视频免费播放|