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

    土體的并聯(lián)理想彈塑性彈簧損傷模型

    2010-10-18 09:25:48肖紅宇
    資源環(huán)境與工程 2010年5期
    關(guān)鍵詞:彈塑性本構(gòu)塑性

    肖紅宇,黃 英

    (1.湖南省水利水電勘測(cè)設(shè)計(jì)研究總院,湖南長(zhǎng)沙 410007;2.昆明理工大學(xué)電力工程學(xué)院,云南昆明 650051)

    土體的并聯(lián)理想彈塑性彈簧損傷模型

    肖紅宇1,黃 英2

    (1.湖南省水利水電勘測(cè)設(shè)計(jì)研究總院,湖南長(zhǎng)沙 410007;2.昆明理工大學(xué)電力工程學(xué)院,云南昆明 650051)

    通過(guò)土顆粒之間的連接其在土體受力過(guò)程中的變形特性,根據(jù)彈塑性理論建立土體的并聯(lián)彈簧模型。結(jié)合土體在受力過(guò)程中處于塑性流動(dòng)的連接數(shù)目和土體中連接總數(shù)目的比值來(lái)建立土體的損傷變量,進(jìn)而推導(dǎo)出土體受力時(shí)的應(yīng)力應(yīng)變關(guān)系。實(shí)例表明,所推導(dǎo)的土體本構(gòu)關(guān)系是合理的。

    本構(gòu)關(guān)系;損傷變量;加載;卸載

    0 引言

    材料的本構(gòu)關(guān)系是反映材料力學(xué)性質(zhì)的數(shù)學(xué)表達(dá)式,巖土材料是一個(gè)復(fù)雜的系統(tǒng),因此其力學(xué)性狀相當(dāng)復(fù)雜,不僅是因?yàn)橥馏w的組成相當(dāng)復(fù)雜,而且影響因素很多,如溫度、濕度、時(shí)間等[1]都對(duì)其有很大的影響。從1773年庫(kù)侖(coulomb)提出土質(zhì)破壞條件(后來(lái)推廣為莫爾—庫(kù)侖準(zhǔn)則)以來(lái),各國(guó)學(xué)者紛紛提出許多巖土的材料的破壞準(zhǔn)則[2],每種破壞準(zhǔn)則都有一定的適用范圍,但沒(méi)有一種可以通用的破壞準(zhǔn)則。在現(xiàn)代的巖土的本構(gòu)關(guān)系中,提出屈服面、塑性勢(shì)、內(nèi)變量等概念,都是以晶體滑移為依據(jù)的,并不一定符合土體的微觀變形機(jī)理,現(xiàn)在的巖土力學(xué)的研究應(yīng)該從微觀和宏觀相結(jié)合的道路上進(jìn)行,在搞清巖土材料的微觀變形機(jī)理的基礎(chǔ)上進(jìn)行[3]。

    損傷力學(xué)力學(xué)是從1958年前蘇聯(lián)學(xué)者Kachanov在金屬儒變的研究中提出的,主要說(shuō)明材料在外部作用下,材料將形成大量的微觀裂紋,這些微觀缺陷通過(guò)形核、擴(kuò)展、匯合造成材料劣化直到破壞[4]。經(jīng)過(guò)幾十年的發(fā)展,損傷力學(xué)已經(jīng)發(fā)展到適用階段。中國(guó)學(xué)者沈珠江在1988年把損傷力學(xué)引入土力學(xué)[5],隨后何開(kāi)勝[6],趙錫宏[7]等在這方面作了研究。本文從土顆粒微觀結(jié)構(gòu)入手,通過(guò)談塑性模型建立了土體的并聯(lián)彈簧模型,從而建立土體受力時(shí)的損傷變量來(lái)探討土體的本構(gòu)關(guān)系,在土體力學(xué)性質(zhì)研究中具有一定的理論價(jià)值。

    1 模型的建立

    1.1 模型假設(shè)

    土體是由土顆粒通過(guò)顆粒之間的連接組成的多孔材料,土體在受力過(guò)程中,土顆粒之間的連接要發(fā)生變形,一方面土顆粒之間的連接發(fā)生彈性變形,當(dāng)應(yīng)力撤消以后,土顆粒又回到其原來(lái)的位置,另一方面,土顆粒之間的連接發(fā)生不可以恢復(fù)的塑性變形,土顆粒在土體中的相對(duì)位置發(fā)生改變,這樣土體產(chǎn)生了塑性變形。由于土體中顆粒的變形一般較小,相對(duì)土體的變形而言,其變形可以忽略不計(jì),因此土顆粒之間的連接的力學(xué)性質(zhì)直接決定了土體的力學(xué)性質(zhì)。為研究的方便,對(duì)土體顆粒之間的連接做如下假設(shè)。

    人們經(jīng)常有多種元件來(lái)描述材料的特性,常用的有彈簧、滑片、粘壺。沈珠江老先生在研究土體損傷力學(xué)的時(shí)候提出了膠結(jié)桿元件用來(lái)描述土體的脆性損傷特性[8],本文將通過(guò)理想彈塑性彈簧元件來(lái)建立土體的彈塑性模型,為研究方便,對(duì)此模型做如下假設(shè):

    (1)土體可以看作由N0個(gè)彈性系數(shù)相同的彈簧并聯(lián)組成的系統(tǒng),如圖1,N0是一個(gè)很大的數(shù)目,每一個(gè)彈簧的彈性系數(shù)為e,且有土體在初始狀態(tài)時(shí)彈性模量為E0=N0e;

    (2)每一個(gè)彈簧是理想彈塑性體,當(dāng)某一彈簧應(yīng)變超過(guò)其極限應(yīng)變值時(shí),彈簧發(fā)生塑性流動(dòng),且其拉伸極限應(yīng)變和壓縮極限應(yīng)變絕對(duì)值相等,如圖2;

    (3)每個(gè)彈簧的極限應(yīng)變值不相等,某一應(yīng)變?chǔ)畔掳l(fā)生塑性流動(dòng)彈簧的數(shù)目服從平均分布,其分布密度函數(shù)為f(ε)=1/εmax,其中εmax為所有彈簧中最大的極限應(yīng)變值,也就是土體的極限應(yīng)變值。

    圖1 土體并聯(lián)彈簧模型Fig.1 The parallel springsmodel of soilmass

    圖2 理想彈塑性模型Fig.2 The ideal elastic-plastic model

    1.2 土體的本構(gòu)關(guān)系

    1.2.1 損傷理論

    損傷力學(xué)認(rèn)為,材料的破壞過(guò)程是材料中微裂慢慢演化過(guò)程,同時(shí)也是一個(gè)從完全未破壞到完全破壞的過(guò)程,描述這一變化過(guò)程的關(guān)系的參量用損傷變量D來(lái)表示,當(dāng)材料處于完全未破壞狀態(tài)時(shí),D=0,當(dāng)材料處于完全破壞狀態(tài)時(shí),D=1。由損傷變量定義可以知道,材料在某一狀態(tài)下力學(xué)特征可以表示為[9]:

    式中:S0表示在完全未損傷狀態(tài)下材料的某一力學(xué)特征參數(shù)值;Sd表示完全損傷狀態(tài)下某一力學(xué)特征參數(shù)值。損傷變量的確定有多種方法,由以上的彈塑性假設(shè)可知,土體在某一應(yīng)變狀態(tài)下時(shí),部分彈簧發(fā)生了塑性流動(dòng),當(dāng)外荷載增加時(shí),這部分彈簧不再承受所增加的應(yīng)力。這樣,我們可以定義損傷變量為處于塑性流動(dòng)狀態(tài)彈簧的數(shù)目和土體中總彈簧的數(shù)目的比值,即有:

    N(ε)為土體處于塑性流動(dòng)的連接的數(shù)量,這樣只要可以確定某一應(yīng)力狀態(tài)下處于塑性流動(dòng)狀態(tài)的彈簧數(shù)目就可以確定土體的損傷變量值。

    1.2.2 損傷狀態(tài)下土體的本構(gòu)關(guān)系

    沈珠江老先生最先把損傷力學(xué)引入土力學(xué),認(rèn)為土體的破壞也是一個(gè)從完全未破壞的到完全破壞的過(guò)程。故和金屬材料一樣,土體的力學(xué)特征參量值也是隨損傷變量演化而變化的,同樣具有如式①的關(guān)系。在研究土體彈性模量時(shí),由模型假設(shè)2可以知道,土體處于完全破壞狀態(tài)下時(shí)不再承受增加的應(yīng)力,即完全損傷狀態(tài)下土體的彈性模量Ed為0,故在某一損傷狀態(tài)下土體的彈性模量具有如下關(guān)系:

    式中:E(ε)為某一應(yīng)力狀態(tài)下土體的彈性模量;E0為完全未破壞的狀態(tài)下土體的彈性模量。

    這樣土體的應(yīng)力應(yīng)變?cè)隽筷P(guān)系可以表示為:

    由于土體受力過(guò)程是一個(gè)損傷演變的過(guò)程,我們認(rèn)為土體在固結(jié)完成未加載的時(shí)候的彈性模量為土體完全未破壞狀態(tài)下的彈性模量,即認(rèn)為土體加載時(shí)的初始彈性模量為土體完全未損傷狀態(tài)下的彈性模量。

    圖3 土體加卸載示意圖Fig.3 The schematic diagram of soilmass on load and unload

    1.2.3 土體加載時(shí)的本構(gòu)關(guān)系

    土體在加載過(guò)程中,如圖3中OA段,當(dāng)加載到某一應(yīng)變狀態(tài)ε(ε<εmax)下,極限應(yīng)變小于ε的彈簧都處于塑性流動(dòng)的狀態(tài),由模型假設(shè)條件3,可以知道,此時(shí)處于塑性流動(dòng)的彈簧數(shù)量為εf(ε)N0,這樣由式②可以建立土體加載的損傷變量為:

    其中:D1為土體加載過(guò)程中的損傷變量;εmax為土體的極限應(yīng)變值。由此可以看出,土體在加載的過(guò)程中某一狀態(tài)下,其損傷變量實(shí)際可以用該狀態(tài)下的應(yīng)變值和土體的極限應(yīng)變值的比值來(lái)決定,這樣只要已知土體的某一狀態(tài)應(yīng)變值和土體的極限應(yīng)變值就可以確定土體的損傷變量。

    將式⑤代入式④可以得到極限應(yīng)變范圍內(nèi)土體加載時(shí)的應(yīng)力應(yīng)變?cè)隽筷P(guān)系式為:

    對(duì)上式積分整理可得:

    其中C1為積分常數(shù),初始加載條件:σ=0,ε=0,把其代入上式可以得到積分常數(shù)C1=εmax2。

    由于ε<εmax,則對(duì)式⑦兩邊開(kāi)方,整理可得極限應(yīng)變范圍內(nèi)土體加載時(shí)的應(yīng)力應(yīng)變關(guān)系表達(dá)式為:

    當(dāng)土體應(yīng)變達(dá)到土體的極限應(yīng)變時(shí),則所有彈簧都發(fā)生塑性流動(dòng),即此時(shí)土體的應(yīng)變趨向于無(wú)窮,而土體的應(yīng)力值保持不變,設(shè)此時(shí)土體所受的應(yīng)力值為σmax,由式⑧可以解得綜合式⑧,可以得到加載時(shí)土體的本構(gòu)關(guān)系:

    有式⑨可以知,只要確定土體的最大應(yīng)變值和土體的初始彈性模量,就可以確定土體的加載關(guān)系。

    1.2.4 土體卸載時(shí)的本構(gòu)關(guān)系

    設(shè)土體在某一應(yīng)變狀態(tài)ε1下卸載(如圖3中曲線BFC,B點(diǎn)為卸載點(diǎn)),當(dāng)處于某一應(yīng)力狀態(tài)ε時(shí),其卸載應(yīng)變?yōu)棣?ε=ε-ε1,此時(shí)部分彈簧由于卸載應(yīng)變超過(guò)了其負(fù)應(yīng)變方向的應(yīng)變極限而發(fā)生反方向塑性流動(dòng)。由假設(shè)3可知,此時(shí)處于塑性流動(dòng)的彈簧數(shù)目由于在卸載過(guò)程中,0<ε<ε1,所以有|Δ ε|=ε1-ε,這樣由式②可以確定土體卸載時(shí)的損傷變量為:

    由此可以看出,卸載時(shí)的損傷變量實(shí)際可以由卸載應(yīng)變值的1/2和最大應(yīng)變值之間的比值來(lái)確定,這樣只要確定了卸載時(shí)的初始應(yīng)變值和土體的極限應(yīng)變值,則可以確定卸載時(shí)某一應(yīng)變狀態(tài)下的損傷變量值。

    由模型假設(shè)2可知土體在卸載時(shí)初始彈性模量和土體加載時(shí)初試彈性模量相等,故卸載時(shí)的初始彈性模量也為Ei,把式⑨代入式④就可以得到土體加載時(shí)的應(yīng)力應(yīng)變?cè)隽筷P(guān)系為:

    對(duì)上式積分可以得到:

    C2為積分常數(shù),設(shè)卸載點(diǎn)ε1所對(duì)應(yīng)的應(yīng)力值為σ1,則卸載的初始條件為σ=σ1,ε=ε1,將其代入式 ?可以得到積分常數(shù)C2:

    同時(shí)由式⑦可以得到在卸載時(shí)ε1對(duì)應(yīng)的應(yīng)力值為:

    把σ1和C1代入式 ?就可以得到:

    對(duì)式 ?變形,把C2代入可以得到土體卸載時(shí)的應(yīng)力應(yīng)變關(guān)系:

    由上式可知當(dāng)應(yīng)力值為0時(shí),此關(guān)系曲線在σ-ε坐標(biāo)系中與ε軸的交點(diǎn)值為,即土體的塑性變形值:

    土體在卸載完成時(shí),土體的應(yīng)變并不等于0,即土體發(fā)生了不可恢復(fù)的變形ε0。由土體的彈塑性理論可知,土體的并不是理想的彈性體,由以上的計(jì)算結(jié)論可知,在卸載時(shí),可以認(rèn)為ε0為土體加載過(guò)程中產(chǎn)生的塑性變形。

    1.2.5 再加載本構(gòu)關(guān)系

    土體卸載后對(duì)其重新加載時(shí),如圖3中CDB段,在沒(méi)有超過(guò)卸載時(shí)最大應(yīng)變時(shí),即ε<ε1時(shí),在卸載時(shí)沒(méi)有發(fā)生塑性流動(dòng)的彈簧在重新加載時(shí)同樣處在彈性范圍之內(nèi),不會(huì)發(fā)生塑性流動(dòng);而卸載時(shí)在負(fù)應(yīng)變方向發(fā)生塑性流動(dòng)的彈簧在再加載的過(guò)程中同樣發(fā)生正應(yīng)變方向塑性流動(dòng),同樣由假設(shè)3可知,土體中在應(yīng)變?chǔ)艩顟B(tài)下,發(fā)生塑性流動(dòng)的彈簧數(shù)目為這樣由式②可以得到土體在再加載時(shí)的損傷變量:

    其中D3為再加載時(shí)土體的損傷變量。此時(shí)的損傷變量實(shí)際上是土體在再加載時(shí)應(yīng)變值的1/2和土體的極限應(yīng)變值之間的比值,這樣只要確定了再加載時(shí)初始狀態(tài)時(shí)的應(yīng)變值和土體的極限應(yīng)變值,就可以確定某一應(yīng)變狀態(tài)下土體的損傷變量。

    同理,由假設(shè)3可知,在再加載時(shí)土體的初始彈性模量也和土體加載時(shí)的彈性模量相等,同樣也為E0,這樣通過(guò)式④可以得到土體再加載時(shí)的應(yīng)力應(yīng)變?cè)隽筷P(guān)系:

    對(duì)上式積分可得:

    C3為積分常數(shù),由式 ?和再加載時(shí)的初始邊界條件σ=0,ε=ε0,可以得到:C3=4εmax2。

    對(duì)式 ?變形可以得到再加載時(shí)應(yīng)力應(yīng)變的關(guān)系為:

    當(dāng)應(yīng)變值超過(guò)卸載時(shí)的最大應(yīng)變時(shí),應(yīng)力應(yīng)變將具有和加載同樣的關(guān)系,結(jié)合式 ?、⑨,這樣可以得到土體的再加載時(shí)的應(yīng)力應(yīng)變關(guān)系為:

    2 實(shí)例驗(yàn)證

    本文引用文獻(xiàn)[1]的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證可以得到,土體的初始彈性模量可以通過(guò)某一實(shí)測(cè)值代如式⑦得到,由實(shí)驗(yàn)資料數(shù)據(jù)代入式⑦,可以得到本土體的初始彈性模量為23.7 MPa;土體的最大應(yīng)變?nèi)⊥馏w實(shí)測(cè)曲線接近塑性流動(dòng)時(shí)的實(shí)驗(yàn)值,由實(shí)驗(yàn)資料可以取土破壞是最大應(yīng)變?yōu)?.471%,本文采用的卸載應(yīng)變?yōu)?.5%,這樣可以得到理論曲線和實(shí)驗(yàn)曲線關(guān)系如圖4。

    圖4 應(yīng)力應(yīng)變關(guān)系圖Fig.4 The graph of stress-strain relation

    由圖可以看出,預(yù)測(cè)曲線和實(shí)驗(yàn)曲線相差很小,基本上可以反映土體的本構(gòu)關(guān)系。同時(shí)土體卸載和回彈曲線組成回滯環(huán),這和巖土的實(shí)驗(yàn)加卸載特性是一致的,同時(shí)土體在回彈和重新加載時(shí)彈性模量和初始彈性模量接近。

    3 結(jié)論

    (1)土體在加卸載表現(xiàn)彈塑性性質(zhì)可以看作有由許多理想彈塑性彈簧共同作用的結(jié)果,通過(guò)理想彈塑性彈簧建立土體彈塑性損傷模型基本上可以反映土體加卸載時(shí)的本構(gòu)關(guān)系;

    (2)由本文土體模型反映土體的回彈和再加載時(shí)組成一個(gè)閉和環(huán),可以反映土體在加卸載時(shí)產(chǎn)生回滯環(huán)的特性;

    (3)在加載過(guò)程中,土體的彈性模量隨著土體的應(yīng)變的增加而減小,當(dāng)達(dá)到一定的應(yīng)變時(shí),土體的彈性模量將趨向于0;

    (4)土體中連接并不一定服從平均分布,以后的研究可以從找出土體最適應(yīng)的分布形式,從而建立更加合適土體變形特性的本構(gòu)關(guān)系。

    [1] 李廣信.高等土力學(xué)[M].北京:清華大學(xué)出版社,2004.

    [2] 鄭穎人,沈珠江,龔曉南.巖土塑性力學(xué)原理[M].北京:中國(guó)建筑工業(yè)出版社,2002.

    [3] 沈珠江.砂土本構(gòu)理論的檢討和重建[M]∥第六界全國(guó)巖土力學(xué)數(shù)值分析與解析方法討論會(huì)論文集.廣州:廣東科學(xué)技術(shù)出版社,1998.

    [4] 余壽文,馮西橋.損傷力學(xué)[M].北京:清華大學(xué)出版社,1996.

    [5] 沈珠江,章為民.損傷力學(xué)在土力學(xué)中的應(yīng)用[M]∥第三屆全國(guó)巖土力學(xué)數(shù)值分析與解析方法討論會(huì)論文集.武漢:武漢測(cè)繪科技大學(xué)出版社,1988.

    [6] 何開(kāi)勝.結(jié)構(gòu)性粘土的微觀變形機(jī)理和彈粘塑損傷模型研究[D].南京:南京水利科學(xué)研究院,2001.

    [7] 趙錫宏,孫紅,羅冠威.損傷土力學(xué)[M].上海:同濟(jì)大學(xué)出版社,2000.

    [8] 沈珠江.土體變形特性的損傷力學(xué)模擬[M]∥第五屆巖土力學(xué)數(shù)值分析與解析方法討論會(huì)論文集.武漢:武漢測(cè)繪科技大學(xué)出版社,1994.

    [9] 沈珠江.結(jié)構(gòu)性粘土的彈塑性損傷模型[J].巖土工程學(xué)報(bào),1993,15(3):21-28.

    (責(zé)任編輯:胡立智)

    Parallel Ideal Elastic-plastic Spring Damage Model of SoilMass

    XIAO Hongyu1,HUANG Ying2
    (1.Hunan Province Hydraulic and PowerDesign Institute,Changsha,Hunan410007;2.Kunm ing University of Science and Technology,School of Electric Power Engineering,Kunm ing,Yunnan650051)

    Based on the defor mation property in the process of soil mass in stress because soils was closed with each other,and the elastic-plastic theory,this paper established the parallel spring model of soil mass.And also built up the damage variable of soil based on the ratio be tween the connection number in flexible flow and the total connection number when soilmasswas in stress,then deduced the stress-strain relation of the soilmass in stress.The instance indicated that the stress-strain relation of soilmass deduced in the paperwas appropriate.

    constitutive relation;damage variable;load;unload

    TU43

    A

    1671-1211(2010)05-0566-05

    2010-07-15;改回日期:2010-09-06

    國(guó)家自然科學(xué)基金資助項(xiàng)目,項(xiàng)目編號(hào):50868009。

    肖紅宇(1979-),男,工程師,勘察技術(shù)與工程專業(yè),從事巖土工程勘察、巖土工程設(shè)計(jì)等方面的工作。E-mail:xhongyy@126.com

    猜你喜歡
    彈塑性本構(gòu)塑性
    基于應(yīng)變梯度的微尺度金屬塑性行為研究
    硬脆材料的塑性域加工
    矮塔斜拉橋彈塑性地震響應(yīng)分析
    鈹材料塑性域加工可行性研究
    離心SC柱混凝土本構(gòu)模型比較研究
    鋸齒形結(jié)構(gòu)面剪切流變及非線性本構(gòu)模型分析
    彈塑性分析在超高層結(jié)構(gòu)設(shè)計(jì)中的應(yīng)用研究
    江西建材(2018年4期)2018-04-10 12:36:52
    一種新型超固結(jié)土三維本構(gòu)模型
    石英玻璃的熱輔助高效塑性域干磨削
    動(dòng)載荷作用下冪硬化彈塑性彎曲裂紋塑性區(qū)
    亚洲精品aⅴ在线观看| 亚洲经典国产精华液单| 亚洲国产毛片av蜜桃av| 啦啦啦啦在线视频资源| av电影中文网址| 久久青草综合色| 亚洲av.av天堂| 亚洲精品乱码久久久久久按摩| videosex国产| 蜜桃在线观看..| 丝袜喷水一区| 简卡轻食公司| 欧美性感艳星| 亚洲精品色激情综合| 一区二区三区精品91| 久久久久久久久久久免费av| 考比视频在线观看| 国产av一区二区精品久久| av专区在线播放| 天堂俺去俺来也www色官网| 老司机亚洲免费影院| 综合色丁香网| 日韩av免费高清视频| 免费播放大片免费观看视频在线观看| 人妻夜夜爽99麻豆av| 岛国毛片在线播放| 国产成人一区二区在线| 国产淫语在线视频| 亚洲欧美清纯卡通| 亚洲精品色激情综合| 精品久久久久久久久亚洲| 成年美女黄网站色视频大全免费 | 亚洲精品自拍成人| 国产有黄有色有爽视频| 欧美变态另类bdsm刘玥| 丝袜在线中文字幕| 最近的中文字幕免费完整| 毛片一级片免费看久久久久| 搡老乐熟女国产| 国产亚洲午夜精品一区二区久久| 精品一区二区免费观看| 成年美女黄网站色视频大全免费 | www.色视频.com| 中文天堂在线官网| 国产白丝娇喘喷水9色精品| 日韩一区二区视频免费看| 国产高清三级在线| 亚洲图色成人| 飞空精品影院首页| 亚洲人与动物交配视频| 美女内射精品一级片tv| 观看美女的网站| 国产片内射在线| 亚洲欧洲日产国产| 大片电影免费在线观看免费| 亚州av有码| 日韩制服骚丝袜av| 国产精品偷伦视频观看了| 欧美精品一区二区大全| 精品国产露脸久久av麻豆| 免费大片黄手机在线观看| 精品一品国产午夜福利视频| 老熟女久久久| 婷婷成人精品国产| 国产精品99久久99久久久不卡 | xxx大片免费视频| 最黄视频免费看| 69精品国产乱码久久久| 亚洲精品久久久久久婷婷小说| 日韩av在线免费看完整版不卡| 永久免费av网站大全| 人成视频在线观看免费观看| 国产永久视频网站| 日韩av不卡免费在线播放| 国产免费一级a男人的天堂| 中文字幕人妻丝袜制服| 91久久精品电影网| 纵有疾风起免费观看全集完整版| 日本黄大片高清| 亚洲av日韩在线播放| 午夜福利网站1000一区二区三区| 制服丝袜香蕉在线| 日韩一本色道免费dvd| 交换朋友夫妻互换小说| 成人毛片a级毛片在线播放| 久久久久久久久久久久大奶| 满18在线观看网站| 天堂中文最新版在线下载| 永久免费av网站大全| 日韩人妻高清精品专区| 欧美亚洲 丝袜 人妻 在线| 国产精品国产三级专区第一集| 亚洲精品日韩在线中文字幕| 男女啪啪激烈高潮av片| 欧美精品人与动牲交sv欧美| 亚洲精品亚洲一区二区| 美女大奶头黄色视频| 亚洲精品亚洲一区二区| 久久免费观看电影| 超色免费av| 五月伊人婷婷丁香| 日本午夜av视频| 99热国产这里只有精品6| 七月丁香在线播放| 2018国产大陆天天弄谢| 九色成人免费人妻av| 2018国产大陆天天弄谢| 大片电影免费在线观看免费| 国产国拍精品亚洲av在线观看| 男女国产视频网站| av国产久精品久网站免费入址| 男男h啪啪无遮挡| 91久久精品国产一区二区成人| 中文字幕av电影在线播放| 丝袜美足系列| av国产精品久久久久影院| 久久99热这里只频精品6学生| 少妇猛男粗大的猛烈进出视频| 99久久精品国产国产毛片| 夜夜骑夜夜射夜夜干| 欧美精品亚洲一区二区| 韩国av在线不卡| 中文字幕制服av| 韩国av在线不卡| 少妇熟女欧美另类| 国产日韩欧美视频二区| 日本黄色片子视频| 午夜精品国产一区二区电影| 蜜桃久久精品国产亚洲av| 乱人伦中国视频| av在线app专区| 亚洲欧美色中文字幕在线| 在线观看免费视频网站a站| 97超碰精品成人国产| 国产成人freesex在线| 波野结衣二区三区在线| 高清午夜精品一区二区三区| 国产成人精品久久久久久| 亚洲一区二区三区欧美精品| 亚洲精品日韩在线中文字幕| 久久精品久久久久久噜噜老黄| 少妇被粗大猛烈的视频| 色婷婷av一区二区三区视频| 大码成人一级视频| 国产片内射在线| 精品卡一卡二卡四卡免费| av女优亚洲男人天堂| 三级国产精品片| 日本黄色片子视频| 婷婷色综合www| 久久久亚洲精品成人影院| 中文欧美无线码| 亚洲,一卡二卡三卡| 久久久欧美国产精品| 欧美97在线视频| 亚洲无线观看免费| 九草在线视频观看| 大话2 男鬼变身卡| 欧美日韩在线观看h| 一本—道久久a久久精品蜜桃钙片| 男女高潮啪啪啪动态图| 91精品国产国语对白视频| 中文乱码字字幕精品一区二区三区| 纯流量卡能插随身wifi吗| 国产淫语在线视频| 精品酒店卫生间| 国产永久视频网站| 亚洲欧洲国产日韩| 午夜久久久在线观看| 午夜福利视频精品| 精品人妻熟女av久视频| 国产欧美日韩一区二区三区在线 | av福利片在线| 又粗又硬又长又爽又黄的视频| 99九九线精品视频在线观看视频| 插阴视频在线观看视频| 久久人人爽人人片av| 亚洲av成人精品一区久久| 色5月婷婷丁香| 亚洲成人av在线免费| 婷婷色麻豆天堂久久| 少妇丰满av| 亚洲综合色网址| 国产av国产精品国产| 狂野欧美激情性bbbbbb| 高清毛片免费看| 午夜免费男女啪啪视频观看| 免费看光身美女| 久久99一区二区三区| 日韩中字成人| 欧美精品国产亚洲| 热re99久久国产66热| 国产 精品1| 精品亚洲成国产av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 22中文网久久字幕| 久久精品国产自在天天线| 国产精品免费大片| 国模一区二区三区四区视频| 国产乱来视频区| 精品亚洲成国产av| 91精品三级在线观看| 亚洲精品aⅴ在线观看| 国产欧美日韩综合在线一区二区| 91久久精品国产一区二区三区| 夜夜骑夜夜射夜夜干| 婷婷成人精品国产| 青春草视频在线免费观看| 亚洲成色77777| 中国三级夫妇交换| 国产精品久久久久成人av| 又黄又爽又刺激的免费视频.| 高清午夜精品一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国模一区二区三区四区视频| 亚洲精品456在线播放app| 亚洲人成网站在线观看播放| 精品酒店卫生间| av在线观看视频网站免费| 777米奇影视久久| 免费观看性生交大片5| 国产老妇伦熟女老妇高清| 久久久久久久久久久丰满| 国产毛片在线视频| 91午夜精品亚洲一区二区三区| 各种免费的搞黄视频| 美女内射精品一级片tv| 啦啦啦啦在线视频资源| 秋霞伦理黄片| 全区人妻精品视频| 啦啦啦啦在线视频资源| 纯流量卡能插随身wifi吗| 成人综合一区亚洲| 精品国产一区二区三区久久久樱花| 街头女战士在线观看网站| 97在线视频观看| 婷婷色综合www| 亚洲美女搞黄在线观看| 色5月婷婷丁香| 午夜福利在线观看免费完整高清在| 日本av免费视频播放| 精品一区在线观看国产| 亚洲国产色片| 久热这里只有精品99| 国产熟女欧美一区二区| 一本色道久久久久久精品综合| 综合色丁香网| 国产亚洲一区二区精品| kizo精华| 亚洲av日韩在线播放| 精品久久久久久电影网| 中国国产av一级| a 毛片基地| 日韩精品免费视频一区二区三区 | 国产乱来视频区| 免费人成在线观看视频色| av.在线天堂| 免费不卡的大黄色大毛片视频在线观看| 午夜福利视频精品| 中文字幕亚洲精品专区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 夫妻午夜视频| 人妻人人澡人人爽人人| 熟妇人妻不卡中文字幕| 欧美日本中文国产一区发布| 婷婷色麻豆天堂久久| 一级a做视频免费观看| 国内精品宾馆在线| 亚洲不卡免费看| 女性生殖器流出的白浆| 一区在线观看完整版| 18禁在线播放成人免费| 下体分泌物呈黄色| 国产亚洲av片在线观看秒播厂| 国产成人精品在线电影| 国产国语露脸激情在线看| 久久久久久久精品精品| 婷婷色av中文字幕| 久久久久久久久久成人| 在线观看免费日韩欧美大片 | 国产成人精品一,二区| 成人亚洲欧美一区二区av| 日韩精品有码人妻一区| 啦啦啦视频在线资源免费观看| 欧美丝袜亚洲另类| 亚洲无线观看免费| 中文欧美无线码| 天堂8中文在线网| 最近2019中文字幕mv第一页| 日本欧美视频一区| 国产免费现黄频在线看| 3wmmmm亚洲av在线观看| 最近中文字幕高清免费大全6| 丝袜在线中文字幕| 99国产精品免费福利视频| a级毛片在线看网站| 中文精品一卡2卡3卡4更新| 一本一本综合久久| 在线观看免费视频网站a站| 国产免费一级a男人的天堂| 大片电影免费在线观看免费| 久久久久网色| 国产午夜精品一二区理论片| 亚洲av中文av极速乱| 男女国产视频网站| 水蜜桃什么品种好| 亚洲人成网站在线观看播放| 免费人妻精品一区二区三区视频| 成人18禁高潮啪啪吃奶动态图 | 蜜桃久久精品国产亚洲av| 久久精品人人爽人人爽视色| 成人黄色视频免费在线看| 国产片内射在线| 久久99热6这里只有精品| 熟女av电影| 欧美日韩视频精品一区| 最近中文字幕高清免费大全6| 国产欧美日韩一区二区三区在线 | av天堂久久9| 免费看av在线观看网站| 亚洲,欧美,日韩| 一区二区av电影网| 亚洲国产精品成人久久小说| 一个人看视频在线观看www免费| 国产精品人妻久久久影院| 亚洲五月色婷婷综合| 岛国毛片在线播放| 日韩制服骚丝袜av| 国产精品熟女久久久久浪| 80岁老熟妇乱子伦牲交| 久久99一区二区三区| 极品人妻少妇av视频| 午夜福利影视在线免费观看| 少妇精品久久久久久久| 天堂中文最新版在线下载| 久久人人爽人人爽人人片va| 少妇 在线观看| 波野结衣二区三区在线| 午夜久久久在线观看| 久久久精品区二区三区| 久久久久久久久久人人人人人人| 亚洲精品国产色婷婷电影| 成年av动漫网址| 99精国产麻豆久久婷婷| 少妇人妻久久综合中文| 午夜老司机福利剧场| 乱码一卡2卡4卡精品| 一级毛片黄色毛片免费观看视频| 一本—道久久a久久精品蜜桃钙片| 热99国产精品久久久久久7| 黄色欧美视频在线观看| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图 | 久久久久久久久久久丰满| 在线观看一区二区三区激情| 国产有黄有色有爽视频| xxx大片免费视频| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av涩爱| 人人妻人人澡人人看| videossex国产| 久久久久久久亚洲中文字幕| 国产成人精品久久久久久| 99久久中文字幕三级久久日本| 免费日韩欧美在线观看| 少妇的逼水好多| 哪个播放器可以免费观看大片| 18禁在线无遮挡免费观看视频| 日本91视频免费播放| 国产精品久久久久久久电影| 国产男女内射视频| 国产乱来视频区| 亚洲中文av在线| 久久久久视频综合| 国产 精品1| 99re6热这里在线精品视频| 亚洲精品,欧美精品| 国产成人免费无遮挡视频| 日本91视频免费播放| 婷婷色麻豆天堂久久| 亚洲精品av麻豆狂野| 新久久久久国产一级毛片| 22中文网久久字幕| 亚洲情色 制服丝袜| 女的被弄到高潮叫床怎么办| 国产av精品麻豆| 国产精品欧美亚洲77777| 超碰97精品在线观看| 五月开心婷婷网| 国产欧美另类精品又又久久亚洲欧美| 黄色视频在线播放观看不卡| 国产伦理片在线播放av一区| 国产一区二区在线观看av| 亚洲av免费高清在线观看| 国产精品国产av在线观看| 人人妻人人澡人人看| a级毛片黄视频| 不卡视频在线观看欧美| 一级毛片 在线播放| 少妇 在线观看| 国产一区二区三区综合在线观看 | 观看av在线不卡| 春色校园在线视频观看| 91精品一卡2卡3卡4卡| 亚洲av男天堂| 亚洲精品第二区| 伦理电影大哥的女人| 中文精品一卡2卡3卡4更新| 国产精品国产三级国产专区5o| 亚洲精品美女久久av网站| 最近最新中文字幕免费大全7| 少妇 在线观看| 久久久久久久久久久免费av| 久久久国产一区二区| 国产欧美日韩一区二区三区在线 | 好男人视频免费观看在线| 美女中出高潮动态图| 欧美日韩在线观看h| 国产老妇伦熟女老妇高清| 十八禁网站网址无遮挡| 国产成人免费观看mmmm| 韩国av在线不卡| 大香蕉久久网| 国产视频首页在线观看| 亚洲少妇的诱惑av| 超色免费av| 免费观看的影片在线观看| 国产精品99久久久久久久久| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久精品古装| 成人毛片60女人毛片免费| 久久99热这里只频精品6学生| 免费人成在线观看视频色| 日韩不卡一区二区三区视频在线| 人妻少妇偷人精品九色| videosex国产| 中国国产av一级| 久久久久久久亚洲中文字幕| 国产 一区精品| 简卡轻食公司| 久久国产精品大桥未久av| 精品一区在线观看国产| 久久久国产欧美日韩av| 免费观看a级毛片全部| 狂野欧美白嫩少妇大欣赏| 99视频精品全部免费 在线| 美女xxoo啪啪120秒动态图| 亚洲国产毛片av蜜桃av| 国产一区二区在线观看日韩| 国产 精品1| 人体艺术视频欧美日本| 黑人欧美特级aaaaaa片| 日本91视频免费播放| 精品国产乱码久久久久久小说| 51国产日韩欧美| 国产精品久久久久久久电影| 啦啦啦中文免费视频观看日本| 亚洲精品久久成人aⅴ小说 | 嫩草影院入口| 日韩成人伦理影院| 爱豆传媒免费全集在线观看| 亚洲国产成人一精品久久久| 天天操日日干夜夜撸| 最近中文字幕2019免费版| 日本wwww免费看| 51国产日韩欧美| 欧美三级亚洲精品| 美女福利国产在线| 日日摸夜夜添夜夜添av毛片| 国产成人精品婷婷| 乱人伦中国视频| 久久精品国产自在天天线| 国产成人aa在线观看| 你懂的网址亚洲精品在线观看| av国产精品久久久久影院| 另类亚洲欧美激情| 亚洲综合色网址| 亚洲av成人精品一区久久| 国产精品国产三级专区第一集| 2021少妇久久久久久久久久久| 人妻制服诱惑在线中文字幕| 制服人妻中文乱码| 日韩视频在线欧美| 色吧在线观看| 飞空精品影院首页| 日本免费在线观看一区| 青春草亚洲视频在线观看| 亚洲精品一区蜜桃| 国产 一区精品| 熟女人妻精品中文字幕| 亚洲av二区三区四区| 国国产精品蜜臀av免费| 中文字幕免费在线视频6| 这个男人来自地球电影免费观看 | 老司机影院成人| 看免费成人av毛片| 欧美国产精品一级二级三级| 乱码一卡2卡4卡精品| 国产 一区精品| 最近的中文字幕免费完整| 高清午夜精品一区二区三区| 亚洲色图综合在线观看| 亚洲精华国产精华液的使用体验| 18禁裸乳无遮挡动漫免费视频| 久热久热在线精品观看| av在线观看视频网站免费| 中文乱码字字幕精品一区二区三区| 国产男女超爽视频在线观看| 黄色怎么调成土黄色| 成年女人在线观看亚洲视频| 七月丁香在线播放| 国产亚洲一区二区精品| 男女啪啪激烈高潮av片| 欧美 亚洲 国产 日韩一| 亚洲av二区三区四区| 免费大片18禁| 日韩一区二区视频免费看| 久久久久视频综合| 国产在线一区二区三区精| 精品国产乱码久久久久久小说| 热re99久久精品国产66热6| 亚洲精品亚洲一区二区| 亚洲欧美中文字幕日韩二区| 亚洲精品aⅴ在线观看| 亚洲人与动物交配视频| 日韩一本色道免费dvd| 亚洲欧美一区二区三区国产| 18在线观看网站| 性色av一级| 亚洲精品日韩在线中文字幕| av女优亚洲男人天堂| 国产av码专区亚洲av| 卡戴珊不雅视频在线播放| 中国美白少妇内射xxxbb| 欧美少妇被猛烈插入视频| 青春草视频在线免费观看| 一本—道久久a久久精品蜜桃钙片| 内地一区二区视频在线| 欧美日韩精品成人综合77777| 一级黄片播放器| 男女边摸边吃奶| 久久久久国产精品人妻一区二区| 91久久精品电影网| 99久国产av精品国产电影| 中文字幕人妻丝袜制服| 国产精品久久久久久av不卡| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 久久亚洲国产成人精品v| 黑人猛操日本美女一级片| 男人添女人高潮全过程视频| 亚洲av福利一区| 老女人水多毛片| 国产欧美另类精品又又久久亚洲欧美| 少妇精品久久久久久久| 国产精品国产三级国产av玫瑰| 国产乱人偷精品视频| 国产一区二区在线观看av| 内地一区二区视频在线| 又大又黄又爽视频免费| 蜜桃久久精品国产亚洲av| av福利片在线| 婷婷色av中文字幕| av免费观看日本| 久久久久久久大尺度免费视频| 亚洲色图 男人天堂 中文字幕 | 黄色毛片三级朝国网站| 国产在线视频一区二区| 午夜福利视频精品| 色网站视频免费| 日本欧美国产在线视频| 午夜福利网站1000一区二区三区| 人妻少妇偷人精品九色| 国产精品 国内视频| 99热全是精品| 亚洲精品久久久久久婷婷小说| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 国产无遮挡羞羞视频在线观看| 成人黄色视频免费在线看| 三上悠亚av全集在线观看| 免费播放大片免费观看视频在线观看| 亚洲精品乱码久久久v下载方式| 国产伦精品一区二区三区视频9| 极品少妇高潮喷水抽搐| 精品久久蜜臀av无| 一区二区日韩欧美中文字幕 | 只有这里有精品99| 免费看av在线观看网站| 婷婷色综合大香蕉| 亚洲精品一区蜜桃| 国产欧美另类精品又又久久亚洲欧美| 精品一区二区三卡| 欧美少妇被猛烈插入视频| 尾随美女入室| 男男h啪啪无遮挡| 在线天堂最新版资源| 2021少妇久久久久久久久久久| 人成视频在线观看免费观看| 国产精品国产三级国产专区5o| 免费观看a级毛片全部| 蜜桃久久精品国产亚洲av| 桃花免费在线播放| 日韩精品免费视频一区二区三区 | 国产一区二区在线观看日韩| 亚洲av在线观看美女高潮| av又黄又爽大尺度在线免费看| av视频免费观看在线观看| 观看美女的网站| 欧美日韩综合久久久久久| 国产片特级美女逼逼视频| 汤姆久久久久久久影院中文字幕| 亚洲精品456在线播放app| 亚洲欧美成人综合另类久久久| 亚洲欧美色中文字幕在线|