• <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ì)量管理
    亚洲欧美精品专区久久| 亚洲精品456在线播放app| 亚洲久久久国产精品| 91精品伊人久久大香线蕉| 看非洲黑人一级黄片| 国产白丝娇喘喷水9色精品| 蜜桃久久精品国产亚洲av| av不卡在线播放| 久久久精品94久久精品| 精品久久久噜噜| 欧美3d第一页| 最近的中文字幕免费完整| 丝袜脚勾引网站| 久久国内精品自在自线图片| 热re99久久国产66热| 久久青草综合色| 成人特级av手机在线观看| 亚洲av成人精品一二三区| 久久精品国产a三级三级三级| 一区二区av电影网| 亚洲四区av| av在线观看视频网站免费| 精品人妻熟女毛片av久久网站| 内地一区二区视频在线| 国产在视频线精品| 亚洲av电影在线观看一区二区三区| 国产真实伦视频高清在线观看| 色网站视频免费| 99热这里只有精品一区| av卡一久久| 久久久久久伊人网av| 亚洲国产成人一精品久久久| 国产 精品1| 简卡轻食公司| 2021少妇久久久久久久久久久| 18禁在线无遮挡免费观看视频| 黄色怎么调成土黄色| 久久久久久久精品精品| 免费不卡的大黄色大毛片视频在线观看| 国产午夜精品久久久久久一区二区三区| 色视频www国产| 亚洲综合精品二区| 国产黄片视频在线免费观看| 97超碰精品成人国产| 亚洲人成网站在线播| 亚洲精品第二区| 精品卡一卡二卡四卡免费| 18禁在线播放成人免费| 日日撸夜夜添| 免费高清在线观看视频在线观看| 高清黄色对白视频在线免费看 | 少妇被粗大的猛进出69影院 | 日韩av免费高清视频| 五月伊人婷婷丁香| 交换朋友夫妻互换小说| 成人免费观看视频高清| 最近中文字幕2019免费版| 中文字幕免费在线视频6| 国产亚洲精品久久久com| 国产黄频视频在线观看| 亚洲国产精品一区二区三区在线| 一个人看视频在线观看www免费| 大陆偷拍与自拍| 国产免费一级a男人的天堂| 深夜a级毛片| av免费观看日本| 一级,二级,三级黄色视频| av.在线天堂| av福利片在线观看| 色网站视频免费| 男女免费视频国产| 欧美日韩一区二区视频在线观看视频在线| 我的老师免费观看完整版| 亚洲av二区三区四区| 日本wwww免费看| 国产成人精品无人区| 天堂中文最新版在线下载| 亚洲电影在线观看av| 黄色日韩在线| 十八禁网站网址无遮挡 | 亚洲国产精品一区二区三区在线| 一级毛片黄色毛片免费观看视频| 乱人伦中国视频| 高清在线视频一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 免费黄频网站在线观看国产| 99久久人妻综合| 国产乱人偷精品视频| 高清在线视频一区二区三区| 在线观看人妻少妇| 欧美一级a爱片免费观看看| 熟妇人妻不卡中文字幕| 在线观看www视频免费| 秋霞在线观看毛片| 日本黄色片子视频| 韩国av在线不卡| 麻豆乱淫一区二区| 日韩亚洲欧美综合| 水蜜桃什么品种好| 久久人人爽人人片av| 啦啦啦啦在线视频资源| 亚洲av不卡在线观看| 三上悠亚av全集在线观看 | 91久久精品电影网| 天天操日日干夜夜撸| 伊人亚洲综合成人网| 久久免费观看电影| 久久免费观看电影| 男女边摸边吃奶| 91午夜精品亚洲一区二区三区| 在线观看av片永久免费下载| 男女啪啪激烈高潮av片| 纵有疾风起免费观看全集完整版| 成人毛片60女人毛片免费| 日本黄大片高清| 亚洲精品aⅴ在线观看| 一区二区三区精品91| 一区在线观看完整版| 久久精品久久久久久噜噜老黄| 欧美日韩精品成人综合77777| 3wmmmm亚洲av在线观看| av免费观看日本| 亚洲美女视频黄频| 亚洲欧洲国产日韩| 97超视频在线观看视频| 我的女老师完整版在线观看| 我的女老师完整版在线观看| 亚洲精品乱久久久久久| 国产av精品麻豆| 三级国产精品欧美在线观看| 国内揄拍国产精品人妻在线| 日本爱情动作片www.在线观看| 亚洲精品乱码久久久v下载方式| 肉色欧美久久久久久久蜜桃| 免费久久久久久久精品成人欧美视频 | 久久久久久久久久成人| 一级毛片 在线播放| 简卡轻食公司| 内射极品少妇av片p| 少妇丰满av| 国语对白做爰xxxⅹ性视频网站| 中文字幕久久专区| 少妇被粗大猛烈的视频| 丝袜脚勾引网站| 高清视频免费观看一区二区| 在线观看三级黄色| 如何舔出高潮| 欧美精品一区二区大全| 午夜福利视频精品| 久久毛片免费看一区二区三区| 久久久午夜欧美精品| 乱系列少妇在线播放| 少妇被粗大的猛进出69影院 | 欧美激情国产日韩精品一区| 国产男人的电影天堂91| 午夜精品国产一区二区电影| 91久久精品电影网| 男男h啪啪无遮挡| 国产色婷婷99| 精品久久久久久久久亚洲| 一级毛片我不卡| 午夜免费观看性视频| 久久国产乱子免费精品| 国产精品嫩草影院av在线观看| 国产在线男女| 亚洲av男天堂| 一级av片app| 肉色欧美久久久久久久蜜桃| 免费观看性生交大片5| 国产亚洲av片在线观看秒播厂| 亚洲激情五月婷婷啪啪| 最后的刺客免费高清国语| a级一级毛片免费在线观看| 69精品国产乱码久久久| 亚洲国产精品成人久久小说| 久久av网站| 亚洲国产精品专区欧美| 99热这里只有是精品在线观看| 黄色一级大片看看| 午夜久久久在线观看| 欧美成人精品欧美一级黄| av.在线天堂| 丰满迷人的少妇在线观看| 九九爱精品视频在线观看| 纵有疾风起免费观看全集完整版| 热re99久久国产66热| 成年女人在线观看亚洲视频| www.色视频.com| 黄色视频在线播放观看不卡| 人体艺术视频欧美日本| 在线观看美女被高潮喷水网站| 久久久久久久国产电影| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 亚洲国产精品成人久久小说| 少妇被粗大的猛进出69影院 | 三级经典国产精品| 我要看日韩黄色一级片| 欧美一级a爱片免费观看看| 如日韩欧美国产精品一区二区三区 | 一级毛片电影观看| 国产色爽女视频免费观看| 日日啪夜夜撸| 国产片特级美女逼逼视频| 欧美日本中文国产一区发布| 欧美另类一区| 亚洲精品一区蜜桃| 日本色播在线视频| 精品人妻熟女av久视频| 你懂的网址亚洲精品在线观看| 日韩强制内射视频| 欧美一级a爱片免费观看看| 一区在线观看完整版| 久久久久网色| 欧美日韩综合久久久久久| av天堂久久9| 在线观看www视频免费| 黄色日韩在线| av在线播放精品| 久久婷婷青草| 日韩三级伦理在线观看| 国产精品伦人一区二区| 男人舔奶头视频| 亚洲婷婷狠狠爱综合网| 国产一区二区三区综合在线观看 | 草草在线视频免费看| 免费黄色在线免费观看| 久久国产精品大桥未久av | 欧美最新免费一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| 七月丁香在线播放| 婷婷色av中文字幕| 在线观看av片永久免费下载| 老司机亚洲免费影院| 国产探花极品一区二区| 国国产精品蜜臀av免费| 国产亚洲av片在线观看秒播厂| 日本av免费视频播放| 亚洲四区av| 岛国毛片在线播放| 观看美女的网站| 男女无遮挡免费网站观看| 日韩精品有码人妻一区| 国产日韩欧美视频二区| 成人毛片60女人毛片免费| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 国产一区亚洲一区在线观看| 人妻系列 视频| 高清午夜精品一区二区三区| 毛片一级片免费看久久久久| 99国产精品免费福利视频| 国产精品免费大片| 丝袜喷水一区| 亚洲av成人精品一二三区| 青春草国产在线视频| 免费观看在线日韩| 18禁在线无遮挡免费观看视频| 色5月婷婷丁香| 国产亚洲欧美精品永久| 99热这里只有是精品50| 七月丁香在线播放| 亚洲精品一区蜜桃| 人人妻人人澡人人看| 久久人人爽人人片av| 制服丝袜香蕉在线| 亚洲欧美清纯卡通| 黑丝袜美女国产一区| 校园人妻丝袜中文字幕| 国产成人午夜福利电影在线观看| 久久人人爽人人爽人人片va| 日韩欧美一区视频在线观看 | videos熟女内射| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 亚洲精品aⅴ在线观看| 成人无遮挡网站| 99国产精品免费福利视频| 欧美少妇被猛烈插入视频| 久久久国产欧美日韩av| 女人精品久久久久毛片| 久久久久久久久久久丰满| 深夜a级毛片| 最近的中文字幕免费完整| 欧美日韩一区二区视频在线观看视频在线| 黄色毛片三级朝国网站 | 欧美日韩国产mv在线观看视频| 精品人妻熟女毛片av久久网站| a级毛片免费高清观看在线播放| 免费黄频网站在线观看国产| 一级毛片 在线播放| 久久久久久久国产电影| 少妇被粗大猛烈的视频| 18禁动态无遮挡网站| 自拍偷自拍亚洲精品老妇| 熟女人妻精品中文字幕| 男的添女的下面高潮视频| 人妻人人澡人人爽人人| 国产av码专区亚洲av| 亚洲欧美日韩卡通动漫| www.av在线官网国产| 18禁在线播放成人免费| 国产伦在线观看视频一区| 亚洲欧洲日产国产| 国产欧美日韩精品一区二区| 国产国拍精品亚洲av在线观看| 日本91视频免费播放| 亚洲av成人精品一二三区| 91精品国产九色| 在线观看免费日韩欧美大片 | 午夜免费男女啪啪视频观看| 一级a做视频免费观看| 国产精品麻豆人妻色哟哟久久| 亚洲欧洲国产日韩| 丝瓜视频免费看黄片| 国产真实伦视频高清在线观看| 99热全是精品| 自拍偷自拍亚洲精品老妇| 亚洲人与动物交配视频| 亚洲av中文av极速乱| 熟女人妻精品中文字幕| 一级,二级,三级黄色视频| 男人爽女人下面视频在线观看| 国产精品久久久久久久电影| 插阴视频在线观看视频| av.在线天堂| 欧美日韩国产mv在线观看视频| 性高湖久久久久久久久免费观看| av卡一久久| 国产精品一区二区三区四区免费观看| 18禁裸乳无遮挡动漫免费视频| 少妇熟女欧美另类| 久热这里只有精品99| 美女中出高潮动态图| 哪个播放器可以免费观看大片| 最近的中文字幕免费完整| 日本-黄色视频高清免费观看| 日韩不卡一区二区三区视频在线| 在线观看免费高清a一片| 国产色婷婷99| 亚洲欧美成人综合另类久久久| 建设人人有责人人尽责人人享有的| 成人亚洲精品一区在线观看| 天堂8中文在线网| 欧美精品高潮呻吟av久久| 精华霜和精华液先用哪个| 黄色视频在线播放观看不卡| 老司机影院毛片| 国产高清国产精品国产三级| 亚洲国产毛片av蜜桃av| 国产精品国产三级专区第一集| av国产久精品久网站免费入址| 97超碰精品成人国产| 成人特级av手机在线观看| 日韩欧美一区视频在线观看 | 免费黄频网站在线观看国产| 亚洲自偷自拍三级| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 国产av一区二区精品久久| 欧美少妇被猛烈插入视频| av视频免费观看在线观看| 日本猛色少妇xxxxx猛交久久| 青春草视频在线免费观看| 又粗又硬又长又爽又黄的视频| 黄色欧美视频在线观看| 高清av免费在线| 亚洲丝袜综合中文字幕| 国产视频首页在线观看| 久久久久久久久久成人| 99re6热这里在线精品视频| 丰满乱子伦码专区| 亚洲无线观看免费| 91在线精品国自产拍蜜月| 免费大片黄手机在线观看| 久久青草综合色| www.av在线官网国产| 久久 成人 亚洲| 三级经典国产精品| 亚洲欧美清纯卡通| 两个人的视频大全免费| 亚洲精品第二区| 免费人妻精品一区二区三区视频| 男人狂女人下面高潮的视频| 日韩不卡一区二区三区视频在线| 老女人水多毛片| 大码成人一级视频| 欧美变态另类bdsm刘玥| 人人妻人人添人人爽欧美一区卜| 国产精品.久久久| 观看美女的网站| 国产欧美亚洲国产| 丰满少妇做爰视频| 亚洲av日韩在线播放| a 毛片基地| 亚洲av国产av综合av卡| 成人二区视频| 99久久精品热视频| 日本欧美视频一区| 国产亚洲最大av| 人妻一区二区av| 亚洲色图综合在线观看| 久久国产精品大桥未久av | 最新的欧美精品一区二区| 久久99热6这里只有精品| 国产 一区精品| 中文字幕免费在线视频6| 菩萨蛮人人尽说江南好唐韦庄| 日本91视频免费播放| 国产亚洲91精品色在线| 97超碰精品成人国产| 观看美女的网站| 性高湖久久久久久久久免费观看| 美女cb高潮喷水在线观看| 午夜福利网站1000一区二区三区| 少妇人妻精品综合一区二区| 欧美日韩视频精品一区| 亚洲国产精品国产精品| 国产成人精品婷婷| 日韩成人伦理影院| 久久 成人 亚洲| 99视频精品全部免费 在线| 男人和女人高潮做爰伦理| 日韩一本色道免费dvd| 少妇熟女欧美另类| 日韩伦理黄色片| 新久久久久国产一级毛片| 亚洲国产欧美在线一区| 免费看不卡的av| 校园人妻丝袜中文字幕| 欧美日韩综合久久久久久| 色婷婷久久久亚洲欧美| 亚洲真实伦在线观看| 亚洲国产欧美日韩在线播放 | 国产在线男女| 久久精品久久精品一区二区三区| 九九在线视频观看精品| 简卡轻食公司| 少妇人妻一区二区三区视频| 精品人妻一区二区三区麻豆| 国产精品免费大片| 99re6热这里在线精品视频| 美女福利国产在线| 在线看a的网站| 久久人人爽人人片av| 观看av在线不卡| 亚洲av在线观看美女高潮| 亚洲国产欧美在线一区| 亚洲无线观看免费| 久久午夜福利片| 久久99一区二区三区| 久久人人爽人人片av| 最黄视频免费看| 亚洲综合色惰| 一级毛片aaaaaa免费看小| 狂野欧美激情性bbbbbb| 国产一区二区三区av在线| av线在线观看网站| 黑人高潮一二区| 新久久久久国产一级毛片| 99热国产这里只有精品6| 观看免费一级毛片| 国产精品免费大片| 精品一品国产午夜福利视频| 亚洲内射少妇av| 亚洲无线观看免费| av国产久精品久网站免费入址| 两个人的视频大全免费| 男人狂女人下面高潮的视频| 中国三级夫妇交换| 亚洲欧美成人精品一区二区| 国产成人免费无遮挡视频| 亚洲欧美精品自产自拍| 国产视频内射| 亚洲国产精品专区欧美| 少妇被粗大的猛进出69影院 | 最新中文字幕久久久久| 欧美一级a爱片免费观看看| 亚洲av电影在线观看一区二区三区| av有码第一页| 亚洲国产精品成人久久小说| 哪个播放器可以免费观看大片| 91成人精品电影| 伦理电影免费视频| 国产白丝娇喘喷水9色精品| 女性被躁到高潮视频| 亚洲精品亚洲一区二区| 亚洲国产成人一精品久久久| 啦啦啦在线观看免费高清www| 午夜福利,免费看| videossex国产| 天天躁夜夜躁狠狠久久av| 免费看日本二区| 18禁裸乳无遮挡动漫免费视频| 免费看不卡的av| 精品少妇内射三级| 久久免费观看电影| 免费黄网站久久成人精品| 大片电影免费在线观看免费| 久久久久久久久久成人| 国产亚洲精品久久久com| 免费高清在线观看视频在线观看| 亚洲精品一区蜜桃| 伦理电影免费视频| .国产精品久久| 亚洲欧美一区二区三区国产| 亚洲欧洲国产日韩| 人妻 亚洲 视频| 一二三四中文在线观看免费高清| 一级毛片aaaaaa免费看小| 91在线精品国自产拍蜜月| 永久网站在线| 嘟嘟电影网在线观看| 各种免费的搞黄视频| 色5月婷婷丁香| 人妻人人澡人人爽人人| 免费人妻精品一区二区三区视频| 国产日韩欧美亚洲二区| 搡老乐熟女国产| 欧美最新免费一区二区三区| 国产中年淑女户外野战色| 日本wwww免费看| 少妇的逼水好多| 精品少妇内射三级| 亚洲av中文av极速乱| 亚洲图色成人| 久久ye,这里只有精品| 国产男女超爽视频在线观看| 51国产日韩欧美| 日本wwww免费看| 久久毛片免费看一区二区三区| 六月丁香七月| 一级,二级,三级黄色视频| 亚洲美女搞黄在线观看| 黑人猛操日本美女一级片| 亚洲av二区三区四区| 久久久久国产精品人妻一区二区| 亚洲伊人久久精品综合| 国产又色又爽无遮挡免| 国产免费一区二区三区四区乱码| 国产极品天堂在线| 在线观看人妻少妇| 午夜免费鲁丝| 亚洲精品aⅴ在线观看| 久久久久久久国产电影| 黑人猛操日本美女一级片| .国产精品久久| 国产日韩欧美视频二区| 国产精品国产三级国产av玫瑰| 人妻一区二区av| 欧美 亚洲 国产 日韩一| 日本色播在线视频| 一级毛片aaaaaa免费看小| 少妇猛男粗大的猛烈进出视频| 日韩av不卡免费在线播放| 亚洲av在线观看美女高潮| 在线观看国产h片| av.在线天堂| 香蕉精品网在线| 国产国拍精品亚洲av在线观看| 伦理电影大哥的女人| 少妇 在线观看| av免费在线看不卡| 久久精品国产亚洲av涩爱| 在线 av 中文字幕| 国产一区二区在线观看日韩| 亚洲人成网站在线观看播放| 亚洲va在线va天堂va国产| 欧美日韩在线观看h| 久久狼人影院| 99久久综合免费| freevideosex欧美| 少妇精品久久久久久久| 亚洲情色 制服丝袜| 在线看a的网站| 内地一区二区视频在线| 久久午夜综合久久蜜桃| 亚洲欧美精品自产自拍| 99热这里只有精品一区| 久久影院123| 大话2 男鬼变身卡| 蜜桃久久精品国产亚洲av| 国产精品一区二区在线不卡| 免费观看在线日韩| 午夜精品国产一区二区电影| 狂野欧美激情性xxxx在线观看| av.在线天堂| 中文字幕亚洲精品专区| 亚洲精品国产色婷婷电影| 肉色欧美久久久久久久蜜桃| 国产欧美日韩一区二区三区在线 | 欧美xxxx性猛交bbbb| 十分钟在线观看高清视频www | 激情五月婷婷亚洲| 国产精品久久久久久久电影| 亚洲国产成人一精品久久久| 噜噜噜噜噜久久久久久91| 永久网站在线| 丝瓜视频免费看黄片| 9色porny在线观看| 日韩欧美 国产精品| 亚洲精品一区蜜桃| 校园人妻丝袜中文字幕| 在线观看av片永久免费下载| 乱码一卡2卡4卡精品| 亚洲,一卡二卡三卡| 日日啪夜夜撸| 国产免费又黄又爽又色| 王馨瑶露胸无遮挡在线观看| 亚洲熟女精品中文字幕| 少妇被粗大的猛进出69影院 | 日韩av免费高清视频|