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

    混凝土板裂紋擴(kuò)展的態(tài)型近場(chǎng)動(dòng)力學(xué)模擬

    2016-11-25 20:42劉一鳴黃丹秦洪遠(yuǎn)
    計(jì)算機(jī)輔助工程 2016年5期
    關(guān)鍵詞:混凝土

    劉一鳴+黃丹+秦洪遠(yuǎn)

    摘要: 構(gòu)建考慮混凝土拉壓異性和宏觀斷裂特征的混凝土類材料非局部態(tài)型近場(chǎng)動(dòng)力學(xué)本構(gòu)模型,并通過引入動(dòng)態(tài)松弛、系統(tǒng)失衡判斷和力邊界等效等算法,構(gòu)建適于分析混凝土類材料和結(jié)構(gòu)變形破壞過程的態(tài)型近場(chǎng)動(dòng)力學(xué)數(shù)值模擬體系.通過分組模擬和定量計(jì)算,分析算法的收斂性、計(jì)算精度和效率等問題;在此基礎(chǔ)上開展含不同角度中心裂紋混凝土板的破壞模擬.

    關(guān)鍵詞: 混凝土; 裂紋擴(kuò)展; 態(tài)型近場(chǎng)動(dòng)力學(xué); 非局部模型

    中圖分類號(hào): O346.1 文獻(xiàn)標(biāo)志碼: A

    0 引 言

    材料和結(jié)構(gòu)的破壞機(jī)制及其數(shù)值模擬是計(jì)算力學(xué)研究的經(jīng)典難題,也是諸多工程領(lǐng)域關(guān)注的重點(diǎn).傳統(tǒng)的有限元等數(shù)值方法由于在連續(xù)介質(zhì)框架下求解偏微分方程,所以在分析破壞問題時(shí)必須預(yù)先知道裂紋的存在與否及其位置和尺寸,計(jì)算時(shí)又需要判斷裂紋是否擴(kuò)展及擴(kuò)展路徑,并重新剖分網(wǎng)格,具有一定的復(fù)雜性.邊界元、擴(kuò)展有限元、非連續(xù)有限元等方法及虛擬裂紋閉合等措施可以較好地處理不連續(xù)問題,但由于依然在連續(xù)性假設(shè)框架下進(jìn)行不連續(xù)區(qū)域的特殊處理,模擬復(fù)雜的諸如三維群裂紋、多尺度動(dòng)力破壞等問題時(shí),依然有待進(jìn)一步研究.

    近年來,基于非局部積分思想的近場(chǎng)動(dòng)力學(xué)(PeriDynamics,PD)方法[1-2]憑借其不需要求解空間微分方程而在模擬大變形及裂紋擴(kuò)展、爆炸和沖擊破壞等強(qiáng)不連續(xù)力學(xué)問題方面的突出優(yōu)勢(shì)[3-5],成為計(jì)算力學(xué)和相關(guān)領(lǐng)域研究的熱點(diǎn).針對(duì)建筑工程中最廣泛使用的混凝土材料和結(jié)構(gòu),SILLING[6],DEMMIE[7]和KILIC等[8-9]采用近場(chǎng)動(dòng)力學(xué)方法模擬簡(jiǎn)單混凝土梁和柱的沖擊破壞、失穩(wěn)等問題,GERSTLE等[10-12]在模擬混凝土梁的基礎(chǔ)上,還進(jìn)一步分析近場(chǎng)動(dòng)力學(xué)模型參數(shù)對(duì)計(jì)算精度的影響.本課題組也曾根據(jù)近場(chǎng)動(dòng)力學(xué)思想建模分析常規(guī)混凝土構(gòu)件的拉、壓和沖擊破壞過程.[13-15]

    然而,已有的相關(guān)工作主要基于常規(guī)“鍵型”近場(chǎng)動(dòng)力學(xué)模型,不論是采用早期的單參數(shù)微彈脆性模型[6-9],還是應(yīng)用改進(jìn)后的微極鍵型近場(chǎng)動(dòng)力學(xué)模型[11-12],均無法真正滿足實(shí)際混凝土材料和結(jié)構(gòu)破壞模擬特別是定量計(jì)算分析的需要(例如,固有的泊松比限制問題[1,8],三維模擬時(shí)單參數(shù)微觀脆性模型中泊松比限定于0.25,而微極模型[11-12]的泊松比則必須小于0.25),對(duì)材料的泊松比具有一定的限制.為修正包括泊松比限制在內(nèi)的鍵型近場(chǎng)動(dòng)力學(xué)模型的缺陷,SILLING等[16-17]進(jìn)一步發(fā)展態(tài)型近場(chǎng)動(dòng)力學(xué)理論,近年來已成為關(guān)注的熱點(diǎn)[18-19],但其相關(guān)工作主要圍繞各種傳統(tǒng)本構(gòu)模型和態(tài)型近場(chǎng)動(dòng)力學(xué)模型之間的轉(zhuǎn)化來開展,對(duì)于混凝土材料和結(jié)構(gòu)破壞這一工程實(shí)際問題的態(tài)型近場(chǎng)動(dòng)力學(xué)模擬尚未見諸報(bào)道.

    本文基于態(tài)型近場(chǎng)動(dòng)力學(xué)理論,構(gòu)建考慮混凝土拉壓異性和宏觀斷裂特征的彈脆性態(tài)型近場(chǎng)動(dòng)力學(xué)本構(gòu)模型,并且通過引入動(dòng)態(tài)松弛、系統(tǒng)失衡判斷、力邊界等效等算法,構(gòu)建完整的統(tǒng)一求解混凝土材料和結(jié)構(gòu)變形破壞問題的數(shù)值體系.通過對(duì)混凝土板單軸拉伸問題的定量計(jì)算和分組模擬,分析方法的收斂性、計(jì)算精度和效率.在此基礎(chǔ)上,應(yīng)用本文的模型和方法模擬含不同角度中心裂紋混凝土板的裂紋擴(kuò)展與破壞過程.

    1 態(tài)型近場(chǎng)動(dòng)力學(xué)理論和模型

    1.1 態(tài)型近場(chǎng)動(dòng)力學(xué)理論

    將空間物質(zhì)視為由帶質(zhì)量、有代表性體積的系列物質(zhì)點(diǎn)x組成,其

    僅與其近場(chǎng)范圍H(|ξ|<δ)內(nèi)的其他物質(zhì)點(diǎn)相互作用,見圖1.圖1中:x和x′分別表示參考構(gòu)形中2個(gè)物質(zhì)點(diǎn)位置矢量;u和u′分別表示當(dāng)前構(gòu)形中2個(gè)物質(zhì)點(diǎn)位移矢量;y=x+u和y′=x′+u′分別表示當(dāng)前構(gòu)形中2個(gè)物質(zhì)點(diǎn)的坐標(biāo)矢量;ξ=x′-x表示2個(gè)物質(zhì)點(diǎn)的相對(duì)位置矢量;η=u′-u表示2個(gè)物質(zhì)點(diǎn)的相對(duì)位移矢量.

    1.2 二維常規(guī)態(tài)型近場(chǎng)動(dòng)力學(xué)本構(gòu)模型

    2 數(shù)值實(shí)現(xiàn)

    2.1 局部阻尼引入

    為便于應(yīng)用近場(chǎng)動(dòng)力學(xué)方法求解準(zhǔn)靜力問題,根據(jù)準(zhǔn)靜力結(jié)果驗(yàn)證模型的可靠型,引入經(jīng)典力學(xué)中的動(dòng)態(tài)松弛法,在運(yùn)動(dòng)方程中引入局部阻尼項(xiàng)[8,15],將物質(zhì)點(diǎn)運(yùn)動(dòng)方程轉(zhuǎn)化為

    2.2 數(shù)值離散

    對(duì)平面問題采用四邊形均勻離散(單位厚度),離散間距為Δx,則物質(zhì)點(diǎn)的代表面積為VΔx=Δx2,對(duì)于物質(zhì)點(diǎn)xi考慮近場(chǎng)范圍內(nèi)(xj-xi≤δ)的相互作用

    2.3 數(shù)值算法

    近場(chǎng)動(dòng)力學(xué)方法的顯式動(dòng)力學(xué)實(shí)現(xiàn)涉及空間離散方式、時(shí)間差分格式及數(shù)值積分方法,數(shù)值積分方法與空間離散方式相關(guān).本文采用中心差分對(duì)時(shí)間序列進(jìn)行離散,

    2.4 系統(tǒng)失衡力準(zhǔn)則

    3 數(shù)值算例

    為驗(yàn)證本文模型和算法,首先通過對(duì)混凝土板的單向拉伸變形進(jìn)行定量計(jì)算分析,進(jìn)而分別對(duì)含不同角度中心裂紋的混凝土板受拉伸時(shí)的裂紋擴(kuò)展和破壞過程進(jìn)行分析.根據(jù)實(shí)際情況,將問題簡(jiǎn)化為平面應(yīng)力問題分析.

    對(duì)粒子系統(tǒng)的力邊界條件施加問題是粒子類方法中的難點(diǎn).對(duì)于長(zhǎng)為L(zhǎng),寬為W的矩形板(見圖2),本文在近場(chǎng)動(dòng)力學(xué)模擬時(shí)在模型的加載邊界上增加寬度為δ的加載區(qū)域,將外力載荷等效為加載區(qū)物質(zhì)點(diǎn)體力施加于加載區(qū)物質(zhì)點(diǎn)上,使模型邊界上應(yīng)力等效邊界條件一致.這樣做還可以消除往常近場(chǎng)動(dòng)力學(xué)計(jì)算中常見的邊界弱化問題.

    3.1 板拉伸定量分析

    取200 mm×200 mm混凝土方板,密度ρ=2 400 kg/m3,彈性模量E=30 GPa,泊松比v=0.2,體積模量和剪切模量分別為K=16.67 GPa,G=12.5 GPa.在豎直方向施加均布拉伸載荷σ=3 MPa.模型離散情況見圖2,離散物質(zhì)點(diǎn)間距Δx=1 mm,模型共包含40 600個(gè)物質(zhì)點(diǎn),其中邊界物質(zhì)點(diǎn)600個(gè),近場(chǎng)范圍取δ取3Δx.使用控制變量法,通過改變時(shí)間步長(zhǎng)、人工阻尼因數(shù)C等計(jì)算參數(shù),觀察迭代過程中板下邊緣中點(diǎn)y方向位移的收斂變化情況.為方便觀察位移隨時(shí)間步變化情況,將收斂容差-設(shè)為0,程序根據(jù)最大迭代步終止計(jì)算.

    (1)時(shí)間步長(zhǎng)取Δt=2×10-7s,改變阻尼系數(shù)C時(shí)的位移收斂情況見圖3.y方向的位移斂于1.011 8×10-5 mm,與彈性力學(xué)解1×10-5 mm相比,誤差為1.18%.

    由此可見,當(dāng)阻尼系數(shù)相對(duì)較小時(shí),位移值在理論解附近劇烈振蕩且無法收斂;隨著阻尼系數(shù)增大,位移振蕩幅值降低,且最終收斂于理論解.當(dāng)阻尼系數(shù)達(dá)到1×108 kg/(m3·s)時(shí),位移值呈單調(diào)收斂,表現(xiàn)出類似超阻尼現(xiàn)象.阻尼系數(shù)繼續(xù)增大,位移收斂逐漸變緩.合適地選取阻尼系數(shù)可以大幅提高準(zhǔn)靜態(tài)求解的計(jì)算效率.根據(jù)本文分組比較與分析,建議將局部阻尼項(xiàng)阻尼系數(shù)取為5×107~5×108 kg/(m3·s).

    (2)取阻尼系數(shù)C=2×108 kg/(m3·s)不變,改變時(shí)間步長(zhǎng)Δt進(jìn)行分析,見圖4.隨著時(shí)間步長(zhǎng)的逐漸增大,位移的收斂也加快,其中時(shí)間步長(zhǎng)為3×10-7s時(shí)收斂速度最快,但當(dāng)時(shí)間步長(zhǎng)取Δt=5×10-7s時(shí),差分迭代出現(xiàn)不收斂.關(guān)于態(tài)型近場(chǎng)動(dòng)力學(xué)計(jì)算中動(dòng)力學(xué)迭代時(shí)間步長(zhǎng)的理論推導(dǎo)目前還需進(jìn)一步深入研究和討論.

    通過對(duì)人工阻尼系數(shù)和時(shí)間步長(zhǎng)尺寸的敏感性分析可知,在基于態(tài)型近場(chǎng)動(dòng)力學(xué)的準(zhǔn)靜力問題分析中,阻尼系數(shù)和時(shí)間步長(zhǎng)的選取對(duì)計(jì)算效率影響非常關(guān)鍵.若選取合適的阻尼系數(shù)和時(shí)間步長(zhǎng),則計(jì)算將快速收斂(如圖3和4).此外,依靠最大迭代數(shù)終止計(jì)算具有很強(qiáng)的不穩(wěn)定性且無法反映計(jì)算結(jié)果的收斂情況,因此引入收斂準(zhǔn)則顯得極為重要.

    根據(jù)式(19)的失衡力準(zhǔn)則,分別選取不同收斂容差,進(jìn)行分組計(jì)算.分別取C=1×108 kg/(m3·s)和Δt=2×10-7 s,計(jì)算結(jié)果見表1:當(dāng)收斂容差-=1×10-4時(shí),計(jì)算誤差接近1%;當(dāng)收斂容差繼續(xù)減小時(shí),計(jì)算誤差僅有細(xì)微變化,而收斂迭代步數(shù)不斷增加.綜合考慮計(jì)算的精度和效率,收斂容差一般取1×10-4即可.

    3.2 含中心裂紋板的裂紋擴(kuò)展模擬

    在圖2所示的混凝土板中心位置設(shè)置20 mm×1 mm的初始裂紋,并在豎直方向施加均布拉伸載荷,計(jì)算混凝土材料的抗壓強(qiáng)度和抗拉強(qiáng)度.對(duì)含不同角度中心裂紋的混凝土板進(jìn)行分級(jí)準(zhǔn)靜力加載,觀察裂紋開展情況,并確定其開裂載荷.時(shí)間步長(zhǎng)取為Δt=2×10-7s,局部阻尼系數(shù)取C=2×108 kg/(m3·s),收斂容差取為-=1×10-4.

    在上述平面模型中,近場(chǎng)范圍為3倍物質(zhì)點(diǎn)間距即δ=3 mm,裂紋寬度小于近場(chǎng)范圍,在形成物質(zhì)點(diǎn)對(duì)時(shí)斷開所有跨越裂紋線的物質(zhì)點(diǎn)對(duì)形成初始裂紋.

    對(duì)于典型的平面I-II復(fù)合型裂紋擴(kuò)展問題,根據(jù)脆性材料的最大拉應(yīng)力準(zhǔn)則可確定起裂角應(yīng)滿足sin θ0+(3cos θ0-1)cot β=0,其中θ0為起裂角,β為裂紋與加載方向夾角.當(dāng)β分別為30°,45°,60°和90°時(shí),理論計(jì)算所得起裂角θ0分別為60°,53.1°,43.2°和0°.

    逐級(jí)加載過程如下:在初始加載時(shí),以0.030 MPa為載荷增量進(jìn)行施加,每次確定系統(tǒng)平衡后施加下一增量步;在系統(tǒng)出現(xiàn)損傷后,退回一個(gè)增量步將載荷增量降低為0.003 MPa進(jìn)行加載,模擬裂紋的穩(wěn)定擴(kuò)展,直至最后失穩(wěn)擴(kuò)展,系統(tǒng)無法平衡,不滿足式(19).

    圖5a~5d分別顯示含45°初始裂紋混凝土板中從裂尖出現(xiàn)損傷、演化成裂紋起裂直至最后裂紋失穩(wěn)擴(kuò)展直至貫穿的全過程.根據(jù)分級(jí)加載算法可以確定,含45°初始裂紋混凝土板的最終破壞載荷為0.441 MPa,即載荷增加到0.441 MPa時(shí),系統(tǒng)無法達(dá)到平衡狀態(tài),無法承受下一載荷增量,隨著動(dòng)力學(xué)方程的繼續(xù)迭代,裂紋自然的快速失穩(wěn)擴(kuò)展.

    其他含不同初始角度裂紋混凝土板的破壞過程與含45°初始裂紋板類似,均由彈性變形、損傷累積、起裂和穩(wěn)定擴(kuò)展、失穩(wěn)擴(kuò)展等階段組成,主要區(qū)別在于最終的破壞載荷以及裂紋擴(kuò)展路徑,詳見圖6.由圖6可知,盡管初始裂紋的傾角不同,裂紋最終基本會(huì)沿著垂直于加載方向快速擴(kuò)展,與試驗(yàn)結(jié)果一致.各種情況下裂紋的起裂角也與根據(jù)最大周向拉應(yīng)力計(jì)算的理論結(jié)果較為吻合.此外,由模擬結(jié)果可知,當(dāng)初始裂紋與加載方向夾角β增大時(shí),I型裂紋所占的比例增加,板的承載能力也隨之下降,最終破壞載荷減小.這一結(jié)果也符合實(shí)際情況.

    4 結(jié) 論

    (1)考慮混凝土類材料的拉壓異性和宏觀斷裂

    特征,建立混凝土新的態(tài)型近場(chǎng)動(dòng)力學(xué)本構(gòu)模型,并通過引入動(dòng)態(tài)松弛、系統(tǒng)失衡判斷、力邊界等效等系列算法,構(gòu)建完整的適合于求解混凝土材料和結(jié)構(gòu)變形破壞過程的態(tài)型近場(chǎng)動(dòng)力學(xué)數(shù)值體系.

    (2)通過對(duì)典型混凝土板的準(zhǔn)靜態(tài)定量變形分析,驗(yàn)證本文模型和算法的收斂性、計(jì)算精度和計(jì)算效率,并討論局部阻尼、時(shí)間步長(zhǎng)和收斂容差等計(jì)算變量的選取對(duì)態(tài)型近場(chǎng)動(dòng)力學(xué)計(jì)算精度和效率的影響,給出合理取值范圍建議.

    (3)通過分級(jí)加載實(shí)現(xiàn)混凝土類材料和結(jié)構(gòu)從彈性變形、損傷累積、裂紋穩(wěn)定擴(kuò)展以及失穩(wěn)破壞全過程的連續(xù)模擬與定量分析,并成功開展含不同角度中心裂紋混凝土板的破壞過程模擬,模擬所得起裂角和裂紋擴(kuò)展路徑與試驗(yàn)結(jié)果和理論值吻合,對(duì)結(jié)構(gòu)破壞載荷的定量計(jì)算結(jié)果與變化規(guī)律符合實(shí)際情況.

    參考文獻(xiàn):

    [1] SILLING S A. Reformulation of elasticity theory for discontinuities and long-range forces[J]. Journal of the Mechanics and Physics of Solids, 2000, 48(1):175-209. DOI: 10.1016/S0022-5096(99)00029-0.

    [2] 黃丹, 章青, 喬丕忠, 等. 近場(chǎng)動(dòng)力學(xué)方法及其應(yīng)用[J]. 力學(xué)進(jìn)展, 2010, 40(4): 448-459.

    HUANG D, ZHANG Q, QIAO P Z, et al. A review on peridynamics(PD) method and its application[J]. Advance in Mechanics, 2010, 40(4): 448-459.

    [3] LAI X, REN B, FAN H, et al. Peridynamics simulations of geomaterial fragmentation by impulse loads[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2015, 39(12):1304-1330. DOI: 10.1002/nag.2356.

    [4] FAN H, BERGEL G L, LI S. A hybrid peridynamics-SPH simulation of soil fragmentation by blast loads of buried explosive[J].International Journal of Impact Engineering, 2016(87): 14-27. DOI: 10.1016/j.ijimpeng.2015.08.006.

    [5] SILLING S A, ASKARI A. Peridynamic model for fatigue cracking, SAND2014-18590[R]. Sandia Report. New Mexico, 2014.

    [6] SILLING S A, DEMMIE P, WARREN T L. Peridynamic simulation of high-rate material failure, SAND2007-3464C[R]. Proceedings of 2007 Applied Mechanics and Materials conference. Austin TX, 2007.

    [7] DEMMIE P N, SILLING S A. An approach to modeling extreme loading of structure using peridynamics[J]. Journal of Mechanics of Materials and Structures, 2007, 2(10): 1921-1945.

    [8] KILIC B. Peridynamic theory for progressive failure prediction in homogeneous and heterogeneous materials[D].Tucson: The University of Arizona, 2008.

    [9] KILIC B, MADENCI E. Structure stability and failure analysis using peridynamic theory[J]. International of Journal of Non-linear Mechanics, 2009, 44(8): 845-854.

    [10] GERSTLE W, SAU N, SILLING S A. Peridynamic modeling of plain and reinforced concrete structures[C]// Proceedings of 18th International Conference Structural Mechanics in Reactor Technology(SMiRT18). Beijing, 2005.

    [11] GERSTLE W, SAU N, AGUILERA E. Micropolar peridynamic constitutive model for concrete[C]// Proceedings of 19th International Conference Structural Mechanics in Reactor Technology (SMiRT19), Toronto,2007: 12-17.

    [12] GERSTLE W H, SAU N, SAKHAVAND N. On peridynamic computational simulation of concrete structures[J]. ACI Special Publication, 2009, 65: 245-264.

    [13] 沈峰, 章青, 黃丹, 等. 基于近場(chǎng)動(dòng)力學(xué)理論的混凝土軸拉破壞過程模擬[J]. 計(jì)算力學(xué)學(xué)報(bào), 2013, 30(S1): 79-83. DOI: 10.7511/jslx2013z017.

    SHEN F, ZHANG Q, HUANG D, et al. Damage and failure process of concrete structure under uni-axialtension based on peridynamics modeling[J]. Chinese Journal of Computational Mechanics, 2013, 30(S1): 79-83. DOI: 10.7511/jslx2013z017.

    [14] HUANG D, ZHANG Q, QIAO P Z. Damage and progressive failure of concrete structures using non-local peridynamic modeling[J].Science China: Technological Science, 2011, 54(3): 591-596. DOI: 10.1007/s11431-011-4306-3.

    [15] HUANG D, LU G, QIAO P. An improved peridynamic approach for quasi-static elastic deformation and brittle fracture analysis[J]. International Journal of Mechanical Sciences, 2015(94/95): 111-122. DOI: 10.1016/j.ijmecsci.2015.02.018.

    [16] SILLING S A, EPTON M, WECKNER O, et al. Peridynamic states and constitutive modeling[J]. Journal of Elasticity, 2007, 88(2): 151-184. DOI: 10.1007/s10659-007-9125-1.

    [17] SILLING S A. Linearized theory of peridynamic states[J]. Journal of Elasticity, 2010, 99(1): 85-111. DOI: 10.1007/s10659-009-9234-0.

    [18] LE Q V, CHAN W K, SCHWARTZ J. A two-dimensional ordinary, state-based peridynamic model for linearly elastic solids[J]. International Journal for Numerical Methods in Engineering 2014, 98(8): 547-561. DOI: 10.1002/nme.4642.

    [19] WU C T. Kinematic constraints in the state-based peridynamics with mixed local/nonlocal gradient approximations[J]. Computational Mechanics, 2014, 54(5): 1255-1267. DOI: 10.1007/s00466-014-1055-8.

    猜你喜歡
    混凝土
    淺談商品混凝土使用控制技術(shù)措施
    混凝土早期變形對(duì)其耐久性的影響
    自我修復(fù)的神奇混凝土
    哈爾濱市觀江國(guó)際小區(qū)地下室大體積混凝土底板溫度計(jì)算及裂縫控制計(jì)算書
    混凝土
    分析膨脹混凝土施工技術(shù)在建筑結(jié)構(gòu)施工中的應(yīng)用
    再生骨料混凝土的基本性能分析研究
    混凝土的質(zhì)量檢測(cè)與控制
    大體積混凝土水化熱及溫度計(jì)算
    加強(qiáng)混凝土的質(zhì)量管理
    国产成人91sexporn| 精品人妻偷拍中文字幕| 久久久a久久爽久久v久久| 亚洲精品久久国产高清桃花| 国产精品免费一区二区三区在线| 久久99热这里只有精品18| 亚洲中文字幕日韩| 国产精品综合久久久久久久免费| 成人毛片a级毛片在线播放| 国产老妇女一区| 日日啪夜夜撸| 欧美日本亚洲视频在线播放| 日日啪夜夜撸| 免费看a级黄色片| 身体一侧抽搐| 国产精品一区二区性色av| 国产免费一级a男人的天堂| 日本黄大片高清| 色哟哟哟哟哟哟| 99久久久亚洲精品蜜臀av| 亚洲中文字幕日韩| 人人妻,人人澡人人爽秒播| 看非洲黑人一级黄片| 香蕉av资源在线| 精品久久久久久久久av| 欧美成人a在线观看| 精品久久久久久久末码| 欧美xxxx性猛交bbbb| 蜜桃亚洲精品一区二区三区| 亚洲国产精品久久男人天堂| 如何舔出高潮| 欧美xxxx性猛交bbbb| 国产人妻一区二区三区在| 国产视频内射| 久久久久久久久中文| 两个人的视频大全免费| 有码 亚洲区| 在线免费十八禁| 精品久久久久久久久久免费视频| 99久久久亚洲精品蜜臀av| 精品午夜福利在线看| 美女cb高潮喷水在线观看| 国产 一区 欧美 日韩| 婷婷精品国产亚洲av| 亚洲成人中文字幕在线播放| 免费搜索国产男女视频| 国产综合懂色| av.在线天堂| 久久精品国产清高在天天线| 给我免费播放毛片高清在线观看| 啦啦啦韩国在线观看视频| 国产午夜福利久久久久久| 在线免费观看的www视频| 亚洲精品久久国产高清桃花| 联通29元200g的流量卡| 12—13女人毛片做爰片一| 黑人高潮一二区| 性色avwww在线观看| 久久人妻av系列| 在线天堂最新版资源| 久久久久久久久大av| 成人无遮挡网站| 国内久久婷婷六月综合欲色啪| 精品福利观看| 高清日韩中文字幕在线| 久99久视频精品免费| 成人欧美大片| 久久人人精品亚洲av| 一级毛片久久久久久久久女| 午夜精品国产一区二区电影 | 日本三级黄在线观看| 国产精品国产高清国产av| 婷婷色综合大香蕉| 69人妻影院| 99久国产av精品| 精品久久久久久久人妻蜜臀av| 18+在线观看网站| 国产在线精品亚洲第一网站| 男人狂女人下面高潮的视频| 亚洲人成网站高清观看| 久久精品夜色国产| 亚洲国产精品成人久久小说 | 日韩高清综合在线| 给我免费播放毛片高清在线观看| 精品人妻熟女av久视频| 日本撒尿小便嘘嘘汇集6| 老女人水多毛片| 一级毛片aaaaaa免费看小| 国产成人一区二区在线| 午夜视频国产福利| 免费看日本二区| 成人二区视频| 99热这里只有是精品在线观看| 成人亚洲精品av一区二区| 一进一出抽搐动态| 国产v大片淫在线免费观看| 别揉我奶头 嗯啊视频| 国产探花在线观看一区二区| 日韩成人伦理影院| 18+在线观看网站| 亚洲国产精品合色在线| 大香蕉久久网| 日本黄色视频三级网站网址| 精品国产三级普通话版| 成年女人永久免费观看视频| 女人十人毛片免费观看3o分钟| 1024手机看黄色片| 在线免费观看不下载黄p国产| 欧美极品一区二区三区四区| 亚洲欧美精品综合久久99| 国产精品一区二区免费欧美| 久久精品国产亚洲av涩爱 | 国产精品亚洲美女久久久| 成年版毛片免费区| 日韩av不卡免费在线播放| 亚洲一级一片aⅴ在线观看| ponron亚洲| 国语自产精品视频在线第100页| 女同久久另类99精品国产91| 久久精品91蜜桃| 看免费成人av毛片| 日日摸夜夜添夜夜添av毛片| 女同久久另类99精品国产91| 久久久久久久久久成人| 人妻制服诱惑在线中文字幕| 天天一区二区日本电影三级| 国产高清三级在线| 一级毛片aaaaaa免费看小| 99精品在免费线老司机午夜| 久久久久精品国产欧美久久久| 欧美中文日本在线观看视频| 1024手机看黄色片| or卡值多少钱| 欧美性猛交╳xxx乱大交人| 国产av不卡久久| 亚洲精品国产成人久久av| 最近视频中文字幕2019在线8| 亚洲熟妇中文字幕五十中出| 国产精品无大码| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产精品国产精品| 成年女人看的毛片在线观看| av在线天堂中文字幕| 内地一区二区视频在线| 成人毛片a级毛片在线播放| 听说在线观看完整版免费高清| 免费看日本二区| 精品不卡国产一区二区三区| 国产国拍精品亚洲av在线观看| 国产日本99.免费观看| 久久久久国产网址| 99在线人妻在线中文字幕| 97超级碰碰碰精品色视频在线观看| 成年版毛片免费区| 一级毛片电影观看 | 免费大片18禁| 国产成人a区在线观看| 91午夜精品亚洲一区二区三区| 亚洲七黄色美女视频| 丝袜喷水一区| 亚洲精品成人久久久久久| 欧美xxxx黑人xx丫x性爽| 成人综合一区亚洲| 国产精品福利在线免费观看| 久久九九热精品免费| 亚洲中文日韩欧美视频| 一进一出好大好爽视频| 免费在线观看成人毛片| 小说图片视频综合网站| 亚洲av一区综合| 国产精品永久免费网站| 免费在线观看影片大全网站| 久久综合国产亚洲精品| 日韩大尺度精品在线看网址| 日本 av在线| 亚洲欧美日韩高清专用| 亚洲色图av天堂| 国产av不卡久久| 久久精品国产99精品国产亚洲性色| 成人亚洲欧美一区二区av| 色综合站精品国产| 久久久久精品国产欧美久久久| 国产乱人视频| 日韩av在线大香蕉| 少妇的逼水好多| 色噜噜av男人的天堂激情| 老司机影院成人| 成人永久免费在线观看视频| 长腿黑丝高跟| 内地一区二区视频在线| 在线观看av片永久免费下载| 一本久久中文字幕| 伊人久久精品亚洲午夜| 亚州av有码| 日本黄色片子视频| 啦啦啦韩国在线观看视频| 全区人妻精品视频| 久久这里只有精品中国| а√天堂www在线а√下载| 美女高潮的动态| 亚洲成av人片在线播放无| 两个人的视频大全免费| 久久中文看片网| videossex国产| 成熟少妇高潮喷水视频| 欧美日韩国产亚洲二区| 免费看a级黄色片| 超碰av人人做人人爽久久| 干丝袜人妻中文字幕| 特级一级黄色大片| 成人三级黄色视频| 岛国在线免费视频观看| 免费黄网站久久成人精品| 亚洲专区国产一区二区| 草草在线视频免费看| 亚洲精品久久国产高清桃花| 久久人妻av系列| 1000部很黄的大片| 久久精品综合一区二区三区| 非洲黑人性xxxx精品又粗又长| 国语自产精品视频在线第100页| 国产精品日韩av在线免费观看| 亚洲久久久久久中文字幕| 中文字幕av成人在线电影| 日日摸夜夜添夜夜爱| 可以在线观看毛片的网站| 看非洲黑人一级黄片| 熟女人妻精品中文字幕| 亚洲av二区三区四区| 久久久久久国产a免费观看| 亚洲国产精品久久男人天堂| 国产 一区精品| 最近在线观看免费完整版| 国产精品一区二区三区四区免费观看 | 亚洲一区高清亚洲精品| 22中文网久久字幕| 给我免费播放毛片高清在线观看| 波野结衣二区三区在线| 精品免费久久久久久久清纯| 精品国内亚洲2022精品成人| 国产精品av视频在线免费观看| 亚洲人与动物交配视频| 国内精品美女久久久久久| 日日摸夜夜添夜夜爱| 少妇的逼水好多| 国产精品一及| 校园春色视频在线观看| 成人三级黄色视频| 久久99热6这里只有精品| 永久网站在线| 插阴视频在线观看视频| 伦理电影大哥的女人| 久久久久久国产a免费观看| 久久久久精品国产欧美久久久| 午夜影院日韩av| 日韩亚洲欧美综合| 97超碰精品成人国产| 中文亚洲av片在线观看爽| 九九爱精品视频在线观看| 国产精品不卡视频一区二区| 欧美区成人在线视频| 国产精品免费一区二区三区在线| 午夜福利视频1000在线观看| 久久久精品大字幕| av在线蜜桃| 麻豆国产97在线/欧美| 国产欧美日韩一区二区精品| 国产探花在线观看一区二区| 精品熟女少妇av免费看| 国产成人91sexporn| 毛片女人毛片| 老师上课跳d突然被开到最大视频| 麻豆久久精品国产亚洲av| 精品久久久久久久久av| 狂野欧美白嫩少妇大欣赏| 国产精品亚洲美女久久久| 国产真实伦视频高清在线观看| 99九九线精品视频在线观看视频| 久久久国产成人精品二区| 成人高潮视频无遮挡免费网站| 精品午夜福利在线看| 成人亚洲精品av一区二区| 日本五十路高清| 最近在线观看免费完整版| 成人亚洲欧美一区二区av| 最近中文字幕高清免费大全6| 成人精品一区二区免费| 少妇的逼好多水| 老女人水多毛片| 男女下面进入的视频免费午夜| 99热这里只有是精品在线观看| 国产高潮美女av| 最近手机中文字幕大全| 国产精品国产高清国产av| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 国产亚洲91精品色在线| 天堂影院成人在线观看| 香蕉av资源在线| 午夜精品国产一区二区电影 | 在线观看av片永久免费下载| 亚洲av免费高清在线观看| 欧美+日韩+精品| 精品国产三级普通话版| 久久草成人影院| 国产精品综合久久久久久久免费| 女人十人毛片免费观看3o分钟| 联通29元200g的流量卡| 亚洲熟妇熟女久久| 免费观看精品视频网站| 变态另类丝袜制服| av专区在线播放| 五月玫瑰六月丁香| 午夜精品一区二区三区免费看| 草草在线视频免费看| 精品日产1卡2卡| 日本五十路高清| 最近的中文字幕免费完整| 麻豆久久精品国产亚洲av| 久久久久九九精品影院| 99久久九九国产精品国产免费| 日本免费a在线| 久久久久久久久中文| 少妇的逼好多水| 亚洲在线自拍视频| 亚洲一级一片aⅴ在线观看| 亚洲美女搞黄在线观看 | 久久人妻av系列| www.色视频.com| 国产91av在线免费观看| 成人国产麻豆网| 女的被弄到高潮叫床怎么办| 久久精品国产自在天天线| 国产69精品久久久久777片| 国产三级中文精品| 在线免费观看不下载黄p国产| 欧美一区二区精品小视频在线| 午夜福利在线观看吧| 国产白丝娇喘喷水9色精品| 精品免费久久久久久久清纯| 看免费成人av毛片| 超碰av人人做人人爽久久| 日韩在线高清观看一区二区三区| 俄罗斯特黄特色一大片| 国产免费男女视频| 狠狠狠狠99中文字幕| 内射极品少妇av片p| 听说在线观看完整版免费高清| 久久久久久久久中文| 久久久久久久午夜电影| 1024手机看黄色片| 午夜日韩欧美国产| 色噜噜av男人的天堂激情| 国产真实乱freesex| 成年版毛片免费区| 国产精品爽爽va在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 一进一出抽搐gif免费好疼| 男女之事视频高清在线观看| 国产成人a区在线观看| 亚洲精品日韩av片在线观看| 国产欧美日韩一区二区精品| 欧美激情久久久久久爽电影| 国产欧美日韩精品亚洲av| 亚洲一区高清亚洲精品| 不卡一级毛片| 欧美激情在线99| 亚洲中文日韩欧美视频| 亚洲美女视频黄频| 在线免费观看不下载黄p国产| 日韩精品中文字幕看吧| 男女边吃奶边做爰视频| 极品教师在线视频| 久久精品人妻少妇| 日韩 亚洲 欧美在线| 天天躁夜夜躁狠狠久久av| 国产真实乱freesex| 午夜福利高清视频| 在线看三级毛片| 国产av在哪里看| 桃色一区二区三区在线观看| 精品人妻视频免费看| 欧美中文日本在线观看视频| 久久精品久久久久久噜噜老黄 | 最近2019中文字幕mv第一页| 波多野结衣高清无吗| 国产日本99.免费观看| 最好的美女福利视频网| 91精品国产九色| 欧美bdsm另类| 亚洲天堂国产精品一区在线| 亚洲av成人av| 国产伦精品一区二区三区四那| 欧美性猛交黑人性爽| 亚洲欧美成人精品一区二区| 亚洲七黄色美女视频| 婷婷亚洲欧美| 高清毛片免费观看视频网站| 成人av一区二区三区在线看| 99热6这里只有精品| 欧美又色又爽又黄视频| 国产v大片淫在线免费观看| 一a级毛片在线观看| 别揉我奶头 嗯啊视频| av免费在线看不卡| 99精品在免费线老司机午夜| 我要看日韩黄色一级片| 日日摸夜夜添夜夜添av毛片| 久久6这里有精品| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品久久国产高清桃花| 两个人的视频大全免费| a级毛色黄片| 亚洲精品日韩av片在线观看| 欧美国产日韩亚洲一区| 国产一区亚洲一区在线观看| 国产激情偷乱视频一区二区| 超碰av人人做人人爽久久| 成人高潮视频无遮挡免费网站| av卡一久久| 高清午夜精品一区二区三区 | 精品国内亚洲2022精品成人| 国产精品爽爽va在线观看网站| 国产伦一二天堂av在线观看| 性插视频无遮挡在线免费观看| 九九在线视频观看精品| 美女内射精品一级片tv| 欧美激情国产日韩精品一区| 97在线视频观看| 淫妇啪啪啪对白视频| 欧美又色又爽又黄视频| 男人的好看免费观看在线视频| 亚洲专区国产一区二区| 国产欧美日韩一区二区精品| av在线观看视频网站免费| 国产精品伦人一区二区| 久久精品国产亚洲av香蕉五月| av在线老鸭窝| 亚洲熟妇熟女久久| 国产伦一二天堂av在线观看| 久久天躁狠狠躁夜夜2o2o| 精品一区二区三区视频在线观看免费| 欧美xxxx性猛交bbbb| 国产精品无大码| 亚洲不卡免费看| 国产一区亚洲一区在线观看| 亚洲精华国产精华液的使用体验 | 性欧美人与动物交配| 久久精品久久久久久噜噜老黄 | 午夜福利成人在线免费观看| 内地一区二区视频在线| 国产私拍福利视频在线观看| 看片在线看免费视频| 超碰av人人做人人爽久久| 精品福利观看| 亚洲国产色片| 国产精品福利在线免费观看| 日本三级黄在线观看| 99热只有精品国产| 啦啦啦韩国在线观看视频| 亚洲性夜色夜夜综合| 婷婷亚洲欧美| 久久午夜亚洲精品久久| 亚洲精品一区av在线观看| 午夜精品国产一区二区电影 | 嫩草影院精品99| 国产精品国产高清国产av| 春色校园在线视频观看| 久久久午夜欧美精品| 欧美成人精品欧美一级黄| 国产成人91sexporn| 成人一区二区视频在线观看| 偷拍熟女少妇极品色| 看十八女毛片水多多多| 亚洲国产精品成人久久小说 | 日本a在线网址| .国产精品久久| 亚洲国产色片| 搡老岳熟女国产| 观看免费一级毛片| 综合色av麻豆| 国产精品一区二区三区四区免费观看 | 欧美日本视频| 夜夜爽天天搞| 国产爱豆传媒在线观看| 一级av片app| 精品一区二区免费观看| 麻豆一二三区av精品| 国产黄色小视频在线观看| 免费av不卡在线播放| 国产人妻一区二区三区在| 亚洲成人精品中文字幕电影| 国产 一区 欧美 日韩| 亚洲精品国产成人久久av| 婷婷亚洲欧美| 国产真实乱freesex| 国产 一区精品| 国产aⅴ精品一区二区三区波| 色哟哟·www| 久久久久国产网址| 国产精品嫩草影院av在线观看| 欧美日本亚洲视频在线播放| 一本精品99久久精品77| 狂野欧美激情性xxxx在线观看| 国产白丝娇喘喷水9色精品| 国产淫片久久久久久久久| 美女大奶头视频| 久久人妻av系列| 中国国产av一级| 变态另类成人亚洲欧美熟女| 国产成人aa在线观看| 一a级毛片在线观看| 啦啦啦观看免费观看视频高清| 亚洲av一区综合| 你懂的网址亚洲精品在线观看 | 1024手机看黄色片| 美女大奶头视频| 淫妇啪啪啪对白视频| 在线观看美女被高潮喷水网站| 九九在线视频观看精品| 国产视频内射| 精品乱码久久久久久99久播| 91在线观看av| 久久久久国内视频| 蜜臀久久99精品久久宅男| 久久精品国产亚洲av天美| 国产黄a三级三级三级人| 日韩,欧美,国产一区二区三区 | 精品久久久久久久末码| 日韩精品青青久久久久久| 99热这里只有精品一区| 久久久久久国产a免费观看| 网址你懂的国产日韩在线| 亚洲成人久久爱视频| 久久精品夜夜夜夜夜久久蜜豆| av专区在线播放| 国产成人a∨麻豆精品| 国产乱人视频| 色尼玛亚洲综合影院| 欧美性猛交黑人性爽| 国产高清不卡午夜福利| 村上凉子中文字幕在线| 91精品国产九色| 在线观看一区二区三区| 中文字幕免费在线视频6| 九九久久精品国产亚洲av麻豆| 欧美国产日韩亚洲一区| 一进一出好大好爽视频| 最近最新中文字幕大全电影3| av免费在线看不卡| 又黄又爽又刺激的免费视频.| 亚洲自偷自拍三级| 国产精品电影一区二区三区| 国产精品日韩av在线免费观看| 少妇的逼水好多| 免费看美女性在线毛片视频| 国产精品乱码一区二三区的特点| 欧美在线一区亚洲| 国产精品1区2区在线观看.| 亚洲精品久久国产高清桃花| .国产精品久久| 免费观看的影片在线观看| 成人毛片a级毛片在线播放| 99热全是精品| 国产真实伦视频高清在线观看| 国产精品久久电影中文字幕| 国产精品一及| 国产欧美日韩精品一区二区| 偷拍熟女少妇极品色| h日本视频在线播放| 亚洲精品国产av成人精品 | 美女 人体艺术 gogo| 亚洲自拍偷在线| 久久精品国产鲁丝片午夜精品| 99热只有精品国产| 少妇的逼水好多| av视频在线观看入口| 成人特级黄色片久久久久久久| 亚洲欧美中文字幕日韩二区| 人人妻,人人澡人人爽秒播| 99热这里只有精品一区| 丰满人妻一区二区三区视频av| a级一级毛片免费在线观看| 搡老妇女老女人老熟妇| 小说图片视频综合网站| 性欧美人与动物交配| 精品久久久久久久久av| 三级男女做爰猛烈吃奶摸视频| 欧美3d第一页| 99在线人妻在线中文字幕| 欧美最新免费一区二区三区| 美女高潮的动态| 少妇的逼好多水| 一级黄色大片毛片| 国产精品伦人一区二区| 国产精品日韩av在线免费观看| 欧美区成人在线视频| 极品教师在线视频| 久久久久久大精品| 国产极品精品免费视频能看的| 午夜精品国产一区二区电影 | 亚洲精品影视一区二区三区av| av在线观看视频网站免费| 黄色欧美视频在线观看| 亚洲精品影视一区二区三区av| 国产精品福利在线免费观看| 国产精品电影一区二区三区| 久久久久国内视频| 婷婷色综合大香蕉| 免费av观看视频| 好男人在线观看高清免费视频| 国产伦精品一区二区三区视频9| 国产亚洲精品av在线| 在现免费观看毛片| 日韩欧美一区二区三区在线观看|