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

    單軸壓縮鹽巖聲發(fā)射特征及損傷演化探討

    2016-05-30 02:59:02周志威劉建鋒徐楊夢迪
    長江科學(xué)院院報(bào) 2016年5期
    關(guān)鍵詞:聲發(fā)射鹽巖

    周志威,劉建鋒,鄒 航,卓 越,徐楊夢迪

    (四川大學(xué) 水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,成都 610065)

    ?

    單軸壓縮鹽巖聲發(fā)射特征及損傷演化探討

    周志威,劉建鋒,鄒航,卓越,徐楊夢迪

    (四川大學(xué) 水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,成都610065)

    摘要:為了探究鹽巖壓縮變形破壞過程中,聲發(fā)射現(xiàn)象與內(nèi)部損傷的相關(guān)關(guān)系,利用四川大學(xué)MTS815巖石力學(xué)試驗(yàn)系統(tǒng),對(duì)云南安寧鹽巖進(jìn)行單軸壓縮試驗(yàn),分析了鹽巖聲發(fā)射特征。依據(jù)損傷理論,得到了基于聲發(fā)射參數(shù)的鹽巖的損傷演化模型,并引入修正系數(shù)對(duì)其進(jìn)行理論修正,分析損傷終值對(duì)損傷本構(gòu)模型的影響。研究結(jié)果表明:鹽巖的聲發(fā)射(acoustic emission,簡稱AE)振鈴計(jì)數(shù)率和能量率在彈性極限附近達(dá)到最大值,之后呈現(xiàn)高頻低幅值的聲發(fā)射現(xiàn)象;基于聲發(fā)射的鹽巖損傷變量值在壓縮變形前期增加較快;基于AE累計(jì)振鈴計(jì)數(shù)得到的損傷值<基于AE累計(jì)能量計(jì)算得到的鹽巖損傷值;修正后的聲發(fā)射損傷理論模型曲線與試驗(yàn)曲線全過程吻合較好,能夠較好地反映鹽巖內(nèi)部損傷演化過程;振鈴計(jì)數(shù)參數(shù)模型比能量參數(shù)模型能更好地模擬鹽巖的損傷演化特征。

    關(guān)鍵詞:鹽巖;聲發(fā)射;損傷模型;損傷變量;修正系數(shù)

    1研究背景

    鹽巖是國際上公認(rèn)的石油、天然氣等能源地下儲(chǔ)備的理想介質(zhì),也是中國能源地下儲(chǔ)備的主要建庫場所。鹽巖損傷會(huì)直接危及儲(chǔ)庫的安全運(yùn)營,因此鹽巖損傷研究是鹽穴能源儲(chǔ)庫研究的重要內(nèi)容之一。謝和平等[1-4]基于循環(huán)荷載試驗(yàn)對(duì)鹽巖損傷演化特征進(jìn)行了探討,揭示了鹽巖應(yīng)力狀態(tài)對(duì)損傷發(fā)展的影響。Schulze等[5-6]根據(jù)鹽巖變形過程中滲透性的變化規(guī)律,對(duì)鹽巖損傷發(fā)展特征進(jìn)行了研究。

    聲發(fā)射(acoustic emission,簡稱AE)作為一種無損監(jiān)測技術(shù),被廣泛應(yīng)用于材料變形破壞全過程的損傷研究。國內(nèi)外利用聲發(fā)射對(duì)巖石和混凝土等材料的損傷進(jìn)行了大量的研究。Ohtsu[7-8]運(yùn)用速率過程理論將聲發(fā)射用于混凝土的損傷評(píng)估當(dāng)中,認(rèn)為聲發(fā)射活性與材料所含裂紋的多少成正相關(guān)。Heiple等[9-10]通過大量實(shí)驗(yàn)分析,認(rèn)為聲發(fā)射振鈴計(jì)數(shù)更能反映材料內(nèi)部損傷演化的過程。Tang和Xu[11]采用連續(xù)損傷力學(xué)方法,從細(xì)觀角度出發(fā),得到了單軸壓縮條件下巖石材料基于聲發(fā)射的損傷模型。國內(nèi)一些學(xué)者利用聲發(fā)射對(duì)煤等巖石的損傷進(jìn)行了研究,揭示了研究巖石的聲發(fā)射特征及基于聲發(fā)射的應(yīng)力-應(yīng)變損傷演化全過程[12-14]。

    損傷理論將橫截面破壞面積與無損狀態(tài)下橫截面積之比定義為損傷變量D。在對(duì)壓縮應(yīng)力狀態(tài)下的巖石損傷演化研究中,為獲得損傷演化模型中的材料參數(shù),通常假定加載前巖石無損傷,即損傷變量初值D=0,當(dāng)巖石達(dá)到抗壓強(qiáng)度時(shí)完全損傷,即損傷變量D=1[15]。事實(shí)上,巖石在壓縮應(yīng)力狀態(tài)下,當(dāng)壓縮應(yīng)力達(dá)到抗壓強(qiáng)度時(shí),橫截面的損傷呈單裂紋或多裂紋破壞特征,仍具有有效承載面積,橫截面并非完全失效,因此,根據(jù)損傷變量定義得到的破壞后的損傷變量D(即損傷終值)應(yīng)<1。假定D=1必然不能反映巖石的真實(shí)損傷特征,為盡可能真實(shí)地反映巖石的損傷演化特征,文獻(xiàn)[12-14]通過聲發(fā)射試驗(yàn),提出了單軸壓縮應(yīng)力下殘余應(yīng)力與峰值應(yīng)力之比的損傷變量修正方法,并對(duì)巖石破壞過程中的損傷本構(gòu)模型進(jìn)行了研究。

    大量研究表明,鹽巖變形破壞過程中的聲發(fā)射特征更為顯著[16-18],因此,基于聲發(fā)射對(duì)鹽巖的損傷研究具有較強(qiáng)的現(xiàn)實(shí)意義,而目前針對(duì)聲發(fā)射特征對(duì)鹽巖損傷演化的研究還相對(duì)較少。

    本文通過鹽巖單軸壓縮聲發(fā)射試驗(yàn),利用損傷修正系數(shù)η對(duì)損傷本構(gòu)模型進(jìn)行修正,分析損傷終值對(duì)損傷本構(gòu)模型的影響,探討利用聲發(fā)射作為描述損傷演化的參量進(jìn)行損傷研究時(shí),聲發(fā)射能量和振鈴計(jì)數(shù)哪一個(gè)更能反映鹽巖的損傷演化特征。

    2基于聲發(fā)射特征參數(shù)的損傷模型

    20世紀(jì)60年代,Kachanov[19]首次提出了連續(xù)度φ這一概念,并用連續(xù)度來描述材料內(nèi)部的逐漸損傷。在連續(xù)度φ的基礎(chǔ)上,對(duì)于單軸應(yīng)力下的各向同性損傷,Kachanov隨后又提出損傷因子(即損傷變量)D,即

    (1)

    式中:A′為損傷后截面有效面積;A為無損狀態(tài)時(shí)的橫截面面積;Ad為破壞截面面積。

    對(duì)于上述表述材料損傷的公式,由于材料的有效承載面積需通過復(fù)雜的統(tǒng)計(jì)分析才能確定,故為了簡化考慮損傷因素后的結(jié)構(gòu)計(jì)算物性方程,1972年,Lemaitre[20]提出了應(yīng)變等效性假設(shè),通過簡單的等效變換可在原有無損材料的本構(gòu)方程基礎(chǔ)上建立損傷材料的本構(gòu)方程,從而在一定程度上促進(jìn)了損傷力學(xué)的快速發(fā)展。該假設(shè)認(rèn)為,對(duì)于彈性材料,在真實(shí)應(yīng)力σ作用下,受損狀態(tài)的應(yīng)變等效于在有效應(yīng)力σ′作用下虛構(gòu)的無損狀態(tài)的應(yīng)變,那么則有

    (2)

    式中:E為變形模量;ε為應(yīng)變。

    假設(shè)斷面微缺陷面積達(dá)到A時(shí)的各累積聲發(fā)射參數(shù)值為N;斷面面積A完全破壞時(shí),即巖石發(fā)生破壞時(shí)的聲發(fā)射參數(shù)累積值為Nm,那么當(dāng)截面損傷破壞面積達(dá)到Ad時(shí),這一時(shí)刻聲發(fā)射參數(shù)累計(jì)值即為

    (3)

    聯(lián)立式(1)和式(3),可以得到損傷變量D與聲發(fā)射參數(shù)值之間的歸一化方程,即

    (4)

    將式(4)代入式(2),即可得出應(yīng)力-應(yīng)變與聲發(fā)射參數(shù)之間的相關(guān)關(guān)系方程,即

    (5)

    通常假定巖石破壞后,試件承載面積全部破壞,事實(shí)上并非如此。

    圖1 鹽巖試驗(yàn)前、后的CT掃描結(jié)果Fig.1 CT scanning images of rock saltbefore and after test

    圖1給出了鹽巖在單軸應(yīng)力下破壞后CT掃描得到的試件內(nèi)部裂紋面分布特征,雖然試件破壞后,產(chǎn)生了大量的裂紋,然而從統(tǒng)計(jì)分析學(xué)來講,試驗(yàn)結(jié)束后試件承載面并未發(fā)生完全破壞,有效面積A′確實(shí)是存在的。該圖表明破壞后的橫截面并非完全失效,故此時(shí)的損傷變量D(損傷終值)應(yīng)<1,同時(shí)也進(jìn)一步說明了單軸壓縮應(yīng)力下,巖石峰后仍具有承載能力的原因。因此,針對(duì)損傷終值<1這一實(shí)際情況,引入損傷修正系數(shù)η(η≤1,本文假定初始時(shí)鹽巖處于無損狀態(tài))對(duì)鹽巖的損傷作出修正,那么可得

    (6)

    (7)

    式(7)為修正后基于聲發(fā)射參數(shù)值的損傷本構(gòu)方程。

    圖2 標(biāo)準(zhǔn)試驗(yàn)試件Fig.2 Standard testspecimens

    3試件制備及試驗(yàn)設(shè)備

    3.1試驗(yàn)試件及制備

    本文鹽巖試件取自云

    南安寧擬建天然儲(chǔ)氣庫鹽礦,平均埋深約500 m。試驗(yàn)參照相關(guān)規(guī)定,采用車床干車法將巖樣制成尺寸為φ90 mm×180 mm的標(biāo)準(zhǔn)圓柱形試件,編號(hào)分別為1#,2#(如圖2)。

    3.2試驗(yàn)設(shè)備

    加載設(shè)備為四川大學(xué)MTS815 Flex Test GT巖石力學(xué)測試系統(tǒng)和PCI-Ⅱ聲發(fā)射(AE)三維定位系統(tǒng)(圖3)。

    圖3 MTS815巖石力學(xué)測試系統(tǒng)Fig.3 Rock mechanics testing system MTS815

    試驗(yàn)軸向荷載傳感器1 000 kN,軸向引伸計(jì)量程-4~4 mm,環(huán)向引伸計(jì)量程為-2.5~8 mm,各測量傳感器精度均達(dá)到標(biāo)定量程點(diǎn)的0.5%。聲發(fā)射監(jiān)測采用8個(gè)傳感器對(duì)稱分布于試件柱面上下兩端,每端4個(gè),傳感器頻率為100 kHz。為保證傳感器與鹽巖試件的耦合效果,在二者接觸部位涂凡士林。

    4試驗(yàn)結(jié)果及分析

    4.1聲發(fā)射特征分析

    鹽巖的應(yīng)力-應(yīng)變?nèi)^程示例曲線和鹽巖試件破壞前、后照片的對(duì)比見圖4和圖5。

    圖4 應(yīng)力-應(yīng)變?nèi)^程曲線Fig.4 Complete stress-strain curves

    圖5 鹽巖破壞后示例Fig.5 Photos of salt rock specimens after failure

    圖6 AE各特征參數(shù)曲線Fig.6 AE characteristic curves

    從圖4可看出,鹽巖的應(yīng)力-應(yīng)變?nèi)^程大致分為3個(gè)階段:①彈性變形階段(0~10 MPa),應(yīng)力增加較快而應(yīng)變?cè)鲩L較慢;②塑性變形階段,隨著應(yīng)力增加,應(yīng)變開始快速增大,加載應(yīng)力上升速度減緩,此過程持續(xù)時(shí)間較長,應(yīng)變值變化較大,塑性應(yīng)變量占整個(gè)應(yīng)變值的80%~90%;③峰后破壞階段,應(yīng)力到達(dá)峰值應(yīng)力后,試件發(fā)生破壞,應(yīng)力值下降一定幅度之后即終止試驗(yàn)。1#試件的峰值應(yīng)力(26.3 MPa)比2#試件的峰值應(yīng)力(23.73 MPa)大10%左右,而峰值應(yīng)力對(duì)應(yīng)的應(yīng)變值大概為2#試件的75%。對(duì)試驗(yàn)后的鹽巖試件進(jìn)行化學(xué)成分分析,得到1#試件的雜質(zhì)含量(約20%)高于2#試件的雜質(zhì)含量(約10%),通常情況下,雜質(zhì)含量越高,鹽巖的強(qiáng)度越大,塑性變形能力相應(yīng)越低[21]。

    由圖5可知,在單軸壓縮條件下,鹽巖試件主要表現(xiàn)為破裂面平行于軸向力方向的貫通張性破壞特征。

    圖6為鹽巖的AE特征曲線。由圖6可看出,對(duì)于鹽巖這種塑性較強(qiáng)的巖石,其聲發(fā)射總體表現(xiàn)為頻度高、強(qiáng)度低。整個(gè)加載試驗(yàn)過程,聲發(fā)射振鈴計(jì)數(shù)率均處于一個(gè)比較高的計(jì)數(shù)水平,最大值都出現(xiàn)在600~700次/s之間,整個(gè)過程并沒有太大的起伏變化,1#試件平均維持在300~400次/s,2#試件整個(gè)過程則維持在400~500次/s,因而累計(jì)振鈴計(jì)數(shù)接近于線性增長。測試鹽巖的聲發(fā)射能量率變化特征基本一致,都是先增大,后逐步減小并最終趨于低能穩(wěn)定,最大值都在28左右,因而累計(jì)能量值在加載前期(彈性階段)快速增長,占整個(gè)試驗(yàn)總值的40%以上,后期則逐漸趨于低速率緩慢增加。

    對(duì)圖6(a)和圖6(b)進(jìn)行細(xì)致分析,發(fā)現(xiàn)鹽巖的聲發(fā)射現(xiàn)象在各個(gè)階段又不大相同。加載初期AE活動(dòng)較為活躍,聲發(fā)射振鈴計(jì)數(shù)和能量隨著加載應(yīng)力增大而增多,試件處在彈性壓密階段,原生缺陷出現(xiàn)閉合、鹽晶晶粒間產(chǎn)生摩擦和滑移而產(chǎn)生聲發(fā)射信號(hào),同時(shí)鹽巖的各向異性和非均勻性也會(huì)引起聲發(fā)射現(xiàn)象。當(dāng)應(yīng)力加載到彈性極限承載力左右時(shí),能量率達(dá)到最大值,此時(shí)鹽巖內(nèi)部初始缺陷基本閉合,開始產(chǎn)生新裂紋以及初始裂紋由于應(yīng)力集中開始擴(kuò)展,試件出現(xiàn)擴(kuò)容現(xiàn)象,因而聲發(fā)射比較活躍,并出現(xiàn)高幅值的能量率[16];隨著應(yīng)力增加,應(yīng)變開始快速增大,試件進(jìn)入塑性變形階段,加載應(yīng)力上升速度減緩,微裂紋的衍生和擴(kuò)展逐漸趨于穩(wěn)定,振鈴計(jì)數(shù)率相對(duì)塑性變形前期稍微減小并最終趨于一個(gè)穩(wěn)定值,能量率顯著降低。王祖蔭[22]認(rèn)為:材料處于穩(wěn)定的位錯(cuò)運(yùn)動(dòng)時(shí),內(nèi)部處于最低能狀態(tài),因而產(chǎn)生高頻度而低(強(qiáng)度)幅值的聲發(fā)射現(xiàn)象。到了加載后期,應(yīng)力逐漸達(dá)到峰值強(qiáng)度,微裂紋逐漸擴(kuò)展并貫通,AE信號(hào)主要由微裂紋貫通及破壞面晶?;飘a(chǎn)生并有所增強(qiáng),但相對(duì)加載前期的聲發(fā)射計(jì)數(shù)頻度和能量率,特別是能量率,已大為降低。

    4.2損傷演化分析

    圖7 煤巖試驗(yàn)擬合結(jié)果[13]Fig.7 Fitting result of coal testing[13]

    文獻(xiàn)[12-13]利用聲發(fā)射研究單軸壓縮下煤損傷本構(gòu)模型時(shí),用峰后殘余應(yīng)力與峰值應(yīng)力之比,并將1-σc/σm(σc為殘余應(yīng)力,σm為峰值應(yīng)力)作為損傷修正因子,對(duì)損傷邊界作出了修正,圖7為修正后理論曲線與試驗(yàn)曲線的對(duì)比。需要指出的是,巖石壓縮試驗(yàn)時(shí),殘余強(qiáng)度的取值本身帶有一定的主觀性,導(dǎo)致試驗(yàn)曲線與理論曲線差異較大。同時(shí),由于巖石峰后繼續(xù)破壞過程中仍有聲發(fā)射信號(hào),峰后取值點(diǎn)不同,不僅導(dǎo)致?lián)p傷理論公式得到不同結(jié)果,而且也不能反映鹽巖的真實(shí)損傷特征,故本文理論公式中Nm均取值至峰值應(yīng)力點(diǎn)。

    根據(jù)文獻(xiàn)[2-3]對(duì)鹽巖損傷的研究,單軸循環(huán)荷載下,鹽巖的損傷終值為0.85~0.95,由于循環(huán)荷載下試件損傷一般較常規(guī)荷載下試件損傷更為嚴(yán)重,并結(jié)合已有鹽巖破壞后CT掃描試驗(yàn)結(jié)果,假定損傷修正系數(shù)η取值為0.80~0.85,得出基于振鈴計(jì)數(shù)和能量參數(shù)修正前、后的損傷變量D與應(yīng)變的關(guān)系如圖8。

    圖8 修正前、后損傷因子D與應(yīng)變的關(guān)系Fig.8 Relationship between damage variable D andstrain before and after modification

    由于聲發(fā)射能量主要集中于加載前期,振鈴計(jì)數(shù)在整個(gè)加載過程比較均勻分布,因而計(jì)算所得的基于振鈴計(jì)數(shù)的理論損傷值要<基于能量的損傷值。式(8)和式(9)分別是基于圖8數(shù)據(jù)擬合得到的1#和2#鹽巖單軸壓縮荷載下修正后的基于聲發(fā)射振鈴計(jì)數(shù)和能量的損傷演化方程。

    1#鹽巖:

    (8)

    2#鹽巖:

    (9)

    將其分別代入式(7)即可得單軸壓縮荷載下鹽巖一維損傷的本構(gòu)方程。

    由圖8可知,基于聲發(fā)射的鹽巖損傷演化大致分為2個(gè)階段:

    (1) 初始損傷階段,主要對(duì)應(yīng)鹽巖的彈性變形階段。損傷理論認(rèn)為,此時(shí)微裂紋和裂隙的發(fā)展較少,變形屬于可恢復(fù)變形,因此,損傷變量應(yīng)該是趨于0的。然而,從圖8可以很清楚地看出,鹽巖在這一階段的損傷變量變化幅度較大,占整個(gè)過程的20%~40%,損傷值在壓縮變形早期增加較快。究其原因:一是受內(nèi)部晶粒之間的膠結(jié)情況以及其他物質(zhì)成分影響,鹽巖的強(qiáng)度、致密程度以及內(nèi)部結(jié)構(gòu)特征與其他巖石具有一定的差異性;二是由于鹽巖在這一階段并不是完全單純的彈性體,原始缺陷的壓密閉合和新裂紋的擴(kuò)展會(huì)產(chǎn)生大量的聲發(fā)射信號(hào)并產(chǎn)生不可逆變形,因而導(dǎo)致其初始損傷變量迅速增加。

    (2) 損傷逐步穩(wěn)定演化發(fā)展階段。從塑性變形前期一直延續(xù)到峰值應(yīng)力出現(xiàn),損傷變量增長幅度逐漸減小并最終趨于穩(wěn)定增加,鹽巖的微裂紋開始擴(kuò)展,同時(shí)產(chǎn)生新的微裂紋,并緩慢發(fā)展貫通,此階段損傷是連續(xù)穩(wěn)定的。由于鹽巖的延展性強(qiáng),密集及高幅值的聲發(fā)射信號(hào)發(fā)生在屈服階段,在破壞斷裂時(shí)產(chǎn)生可察覺但幅值不高的聲發(fā)射信號(hào)[22],因而在破壞階段,損傷變量變化不顯著,即沒有明顯的損傷加速階段出現(xiàn)。

    根據(jù)基于聲發(fā)射參數(shù)的損傷模型理論公式(5)和修正公式(7),得到基于聲發(fā)射參數(shù)值的應(yīng)力-應(yīng)變修正前后的曲線,并將其與試驗(yàn)結(jié)果作對(duì)比(見圖9)。

    圖9 修正前、后理論曲線與試驗(yàn)結(jié)果對(duì)比Fig.9 Comparison between theoretical curves andtest results before and after modification

    從圖9可以看出,未經(jīng)修正的聲發(fā)射模型與試驗(yàn)情況相差較大,得到的峰值應(yīng)力值也明顯小于試驗(yàn)所測得的峰值應(yīng)力;而經(jīng)過修正后的理論曲線在整個(gè)階段都能較好地與試驗(yàn)曲線吻合,修正后的理論峰值應(yīng)力值也與試驗(yàn)值非常接近。對(duì)比圖7和圖9理論曲線與試驗(yàn)曲線的吻合度可知,本文針對(duì)損傷終值提出的修正系數(shù)η及取值是有效的??梢?,聲發(fā)射作為鹽巖變形破壞過程中的伴生現(xiàn)象,如果獲得破壞時(shí)的真實(shí)損傷,便能夠較好地用來描述鹽巖壓縮狀態(tài)下的應(yīng)力-應(yīng)變損傷演化發(fā)展特征。

    由圖9對(duì)修正后基于振鈴計(jì)數(shù)和能量的曲線作比較,發(fā)現(xiàn)相對(duì)于聲發(fā)射能量參數(shù)的理論模型,基于聲發(fā)射振鈴計(jì)數(shù)的理論模型與試驗(yàn)曲線更為吻合,表明基于振鈴計(jì)數(shù)參數(shù)的模型比基于能量參數(shù)的模型更適合用于描述鹽巖應(yīng)力-應(yīng)變損傷全過程。Heiple和Carpenter[9]認(rèn)為:振鈴計(jì)數(shù)與材料中位錯(cuò)運(yùn)動(dòng)、夾雜物、第二相粒子剝離斷裂及裂紋擴(kuò)展所釋放的應(yīng)變能呈一定的相關(guān)關(guān)系,因而在一定程度上能夠更好地反映材料的損傷演化過程。

    4.3討論

    巖石損傷量的準(zhǔn)確試驗(yàn)測試具有很大困難。為準(zhǔn)確測定損傷,任建喜、葛修潤等[23-24]利用CT對(duì)巖石損傷分布特征、變形破壞過程中的損傷閾值等進(jìn)行了分析,取得了豐富研究成果。本文僅根據(jù)現(xiàn)有的鹽巖損傷研究文獻(xiàn)及已有CT掃描結(jié)果,對(duì)損傷修正系數(shù)取值進(jìn)行假定,未能實(shí)現(xiàn)聲發(fā)射與CT掃描的同步試驗(yàn)測試。鑒于CT掃描可為準(zhǔn)確獲得巖石破壞后各橫截面損傷及整體損傷量值提供有效測試手段,若能利用CT掃描獲得真實(shí)的損傷修正系數(shù),得到的本構(gòu)模型將更能反映實(shí)際損傷演化特征。筆者所在實(shí)驗(yàn)室已引進(jìn)了加載式工業(yè)CT,對(duì)巖石變形破壞全過程的損傷演化特征進(jìn)行CT掃描與聲發(fā)射三維定位的同步實(shí)時(shí)監(jiān)測將在今后進(jìn)行系統(tǒng)深入研究。

    5結(jié)論

    (1) 單軸壓縮下,鹽巖的聲發(fā)射振鈴計(jì)數(shù)率在整個(gè)試驗(yàn)過程維持在一個(gè)較高的水平,整體變化幅度不大,能量率在彈性極限附近達(dá)到最大值,塑性變形后期至破壞階段均呈現(xiàn)高頻低幅值的聲發(fā)射現(xiàn)象。

    (2) 變形初期階段,損傷變量值增長速率較快;基于能量計(jì)算得到的損傷變量值要大于基于振鈴計(jì)數(shù)的損傷變量值;引入修正系數(shù)前的理論曲線與試驗(yàn)曲線相差較大,且峰值應(yīng)力明顯小于試驗(yàn)值,針對(duì)損傷終值,引入修正系數(shù)后的理論曲線與試驗(yàn)曲線吻合較好,能更好地模擬鹽巖應(yīng)力-應(yīng)變損傷全過程。

    (3) 基于振鈴計(jì)數(shù)參數(shù)的模型比基于能量參數(shù)的理論模型計(jì)算得到的理論曲線更接近于試驗(yàn)結(jié)果,說明采用振鈴計(jì)數(shù)這一聲發(fā)射參數(shù)能夠較好地模擬鹽巖的應(yīng)力-應(yīng)變損傷過程。

    參考文獻(xiàn):

    [1]XIE He-ping,LIU Jian-feng,JU Yang,etal. Fractal Property of Spatial Distribution of Acoustic Emissions during the Failure Process of Bedded Rock Salt[J]. International Journal of Rock Mechanics and Mining Sciences, 2011, 48(8): 1344-1351.

    [2]LIU Jian-feng,XIE He-ping,HOU Zheng-meng,etal. Damage Evolution of Rock Salt under Cyclic Loading in Uniaxial Tests[J]. Acta Geotechnica,2014, 9(1):153-160.

    [3]劉建鋒,裴建良,張茹,等. 基于多級(jí)荷載試驗(yàn)的巖石損傷模量探討[J]. 巖石力學(xué)與工程學(xué)報(bào),2012,31(增1):3145-3151. [4]郭印同,趙克烈,孫冠華,等.周期荷載下鹽巖的疲勞變形及損傷特性研究[J]. 巖土力學(xué),2011,32(5):1353-1359.

    [5]SCHULZE O, POPP T, KERN H. Development of Damage and Permeability in Deforming Rock Salt[J]. Engineering Geology, 2001, 61(2/3): 163-180.

    [6]STORMONT J C. In situ Gas Permeability Measurements to Delineate Damage in Rock Salt[J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(7): 1055-1064.[7]OHTSU M. Acoustic Emission Characteristics in Concrete and Diagnostic Application [J]. Acoustic Emission, 1987, 6(2): 99-108.[8]OHTSU M. Rate Process Analysis of Acoustic Emission Activity in Core Test of Concrete[J]. Concrete, 1992 (20):143-153 .

    [9]HEIPLE C R,CARPENTER S H.Acoustic Emission from Dislocation Motion[M]. In: Acoustic Emission, New York: Gordon and Breach Science Publishers, 1983.

    [10]WADLEY H N G, SCRUBY C B, SPEAKE J H. Acoustic Emission for Physical Examination of Metals[J]. International Metals Reviews, 1980, 25(1): 41-64.

    [11]TANG C A,XU X H. Evolution and Propagation of Material Defects and Kaiser Effect Function[J]. Journal of Seismological Research,1990,13(2):203-213.

    [12]劉保縣,黃敬林,王澤云,等.單軸壓縮煤巖損傷演化及聲發(fā)射特征研究[J].巖石力學(xué)與工程學(xué)報(bào),2009,28(增1):3234-3238.

    [13]高保彬,李回貴. 基于聲發(fā)射參量的煤樣損傷模型研究[J].防災(zāi)減災(zāi)工程學(xué)報(bào),2014,34(1):101-106.

    [14]楊永杰,王德超,郭明福,等. 基于三軸壓縮聲發(fā)射試驗(yàn)的巖石損傷特征研究[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(1):98-104.

    [15]謝和平. 巖石混凝土損傷力學(xué)[M]. 徐州:中國礦業(yè)大學(xué)出版社,1990.

    [16]姜德義,陳結(jié),任松,等. 鹽巖單軸應(yīng)變率效應(yīng)與聲發(fā)射特征試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(2):326-336.

    [17]任松,白月明,姜德義,等. 周期荷載作用下鹽巖聲發(fā)射特征試驗(yàn)研究[J].巖土力學(xué),2012,33(6):1613-1618.

    [18]李浩然,楊春和,劉玉剛,等.單軸荷載作用下鹽巖聲波與聲發(fā)射特征試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(10):2107-2116.

    [19]KACHANOV L M. On the Time to Failure under Creep Conditions [J]. Izv AN SSSR Otd Tekhn Nauk,1958, 8(2): 26-31.

    [20]LEMAITURE J. Evaluation of Dissipation and Damage in Metals Submitted to Dynamic Loading[C]∥Proceedings of the ICM-1. Kyoto, Japan, 1972:540-549.

    [21]鄭雅麗,張華賓,王芝銀,等. 含雜質(zhì)鹽巖力學(xué)特性對(duì)比試驗(yàn)研究[J].煤炭學(xué)報(bào),2012,37(1):17-20.

    [22]王祖蔭.聲發(fā)射技術(shù)基礎(chǔ)[M].濟(jì)南:山東科學(xué)技術(shù)出版社,1990.

    [23]任建喜,葛修潤.單軸壓縮巖石損傷演化細(xì)觀機(jī)理及其本構(gòu)模型研究[J].巖石力學(xué)與工程學(xué)報(bào),2001,20(4):425-431.

    [24]葛修潤,任建喜,蒲毅彬,等. 巖石疲勞損傷擴(kuò)展規(guī)律CT細(xì)觀分析初探[J].巖土工程學(xué)報(bào),2001,23(2):191-195.

    (編輯:趙衛(wèi)兵)

    Acoustic Emission Characteristics and Damage Evolution ofRock Salt under Uniaxial Compression

    ZHOU Zhi-wei , LIU Jian-feng, ZOU Hang, ZHUO Yue, XU Yang-meng-di

    (State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu610065, China)

    Abstract:In order to explore the relation between acoustic emission phenomenon and internal damage in the failure process of compression deformation of rock salt, we conducted uniaxial compression test on salt rock from Anning, Yunnan province using MTS815 rock mechanics test system of Sichuan University. The damage evolution model of rock salt was established with acoustic emission parameters based on damage theory and the damage variable D was modified by coefficient of correction. The research results show that the maximum values of ring count rate and energy rate appeared near the elastic limit and the phenomenon of high frequency and low amplitude appeared later; the damage value based on AE increased rapidly in the early stage of compressive deformation, and the damage value based on AE accumulated ring counts is smaller than that based on cumulative energy; the revised damage model curve is in agreement with the whole process of test curve and could reflect the damage evolution process of rock salt well.Simulation using ring count can reflect the damage evolution characteristics of rock salt better than using energy.

    Key words:rock salt;acoustic emission;damage model ;damage variable;coefficient of correction

    中圖分類號(hào):TU45

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1001-5485(2016)05-0063-06

    doi:10.11988/ckyyb.201501682016,33(05):63-68,75

    作者簡介:周志威(1990-),男,四川樂山人,碩士研究生,主要從事巖石力學(xué)與工程方面的研究工作,(電話)15882485585(電子信箱)154762469@qq.com。通訊作者:劉建鋒(1979-),男,河南新鄉(xiāng)人,副教授,博士,研究方向?yàn)閹r石力學(xué)與工程,(電話)13550106392(電子信箱)ljfscu@163.com。

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目(51120145001,51104101,51374148)

    收稿日期:2015-03-10;修回日期:2015-04-14

    猜你喜歡
    聲發(fā)射鹽巖
    水熱綜合作用下鈣芒硝鹽巖強(qiáng)度等參數(shù)的衰減規(guī)律研究*
    鹽巖路基工程特性研究進(jìn)展
    公路工程(2021年6期)2021-02-14 12:34:06
    鹽巖巴西劈裂損傷愈合特性實(shí)驗(yàn)研究
    鹽巖儲(chǔ)庫腔底堆積物空隙體積試驗(yàn)與計(jì)算
    聲發(fā)射技術(shù)在壓力容器中的應(yīng)用與發(fā)展
    濾波器在聲發(fā)射信號(hào)噪聲預(yù)處理中的應(yīng)用
    典型鹽巖力學(xué)特性分析
    修正的鹽巖擴(kuò)容模型及擴(kuò)容界限研究
    基于LPC1768的聲發(fā)射數(shù)據(jù)采集系統(tǒng)設(shè)計(jì)
    聲發(fā)射信號(hào)處理關(guān)鍵技術(shù)研究
    观看美女的网站| 久久精品综合一区二区三区| av中文乱码字幕在线| 99久久久亚洲精品蜜臀av| 不卡一级毛片| 波野结衣二区三区在线| 欧美日韩瑟瑟在线播放| 每晚都被弄得嗷嗷叫到高潮| 国产精品影院久久| 一进一出好大好爽视频| 国产av不卡久久| av在线天堂中文字幕| 免费搜索国产男女视频| 在线看三级毛片| 男人狂女人下面高潮的视频| 午夜久久久久精精品| 亚洲av中文字字幕乱码综合| 欧美精品啪啪一区二区三区| 欧美最新免费一区二区三区 | 国产aⅴ精品一区二区三区波| 精品久久久久久久久久免费视频| 国产伦精品一区二区三区四那| 天天一区二区日本电影三级| 国产极品精品免费视频能看的| 日韩精品青青久久久久久| 少妇被粗大猛烈的视频| 精品久久国产蜜桃| 久久天躁狠狠躁夜夜2o2o| 最近在线观看免费完整版| 欧美3d第一页| 国产色婷婷99| 中文亚洲av片在线观看爽| 欧美日本亚洲视频在线播放| 亚洲自偷自拍三级| 亚洲国产欧洲综合997久久,| 一区二区三区四区激情视频 | 免费在线观看成人毛片| 尤物成人国产欧美一区二区三区| 国产欧美日韩精品一区二区| 国产老妇女一区| 亚洲自拍偷在线| 国产av不卡久久| 亚洲av一区综合| 免费观看的影片在线观看| 国产麻豆成人av免费视频| 美女高潮的动态| 性欧美人与动物交配| 久久草成人影院| 国产久久久一区二区三区| 极品教师在线视频| 国产又黄又爽又无遮挡在线| av在线老鸭窝| 免费人成在线观看视频色| 午夜免费男女啪啪视频观看 | 男插女下体视频免费在线播放| 日韩有码中文字幕| 日韩欧美在线乱码| 亚洲欧美日韩无卡精品| 综合色av麻豆| 国产日本99.免费观看| 亚洲av中文字字幕乱码综合| 人妻丰满熟妇av一区二区三区| 欧美日本视频| 亚洲精品日韩av片在线观看| 欧美激情在线99| 首页视频小说图片口味搜索| 色5月婷婷丁香| 99国产极品粉嫩在线观看| 琪琪午夜伦伦电影理论片6080| 久久午夜亚洲精品久久| 国产探花在线观看一区二区| 欧美又色又爽又黄视频| 深爱激情五月婷婷| 欧美日韩亚洲国产一区二区在线观看| 久久久久久久久久成人| 五月伊人婷婷丁香| 中国美女看黄片| 亚洲美女搞黄在线观看 | 亚洲成av人片在线播放无| 91字幕亚洲| 国产一区二区激情短视频| 简卡轻食公司| 国产成+人综合+亚洲专区| 国产午夜精品久久久久久一区二区三区 | 久久久色成人| 国产精品野战在线观看| 国产午夜精品久久久久久一区二区三区 | 99热这里只有精品一区| 一个人观看的视频www高清免费观看| 日本 欧美在线| 久久久久国内视频| 18禁黄网站禁片午夜丰满| 亚洲人成网站在线播| 国产综合懂色| 欧美高清成人免费视频www| 欧美日韩瑟瑟在线播放| 免费在线观看影片大全网站| 自拍偷自拍亚洲精品老妇| 男女做爰动态图高潮gif福利片| 亚洲熟妇熟女久久| 91在线精品国自产拍蜜月| 久久久色成人| 精品久久久久久久久久久久久| 一本综合久久免费| 少妇人妻一区二区三区视频| 欧美成人a在线观看| 欧美国产日韩亚洲一区| 身体一侧抽搐| 成人特级av手机在线观看| 国产成人aa在线观看| 99久久99久久久精品蜜桃| 亚洲精品一区av在线观看| 国产黄色小视频在线观看| 亚洲欧美精品综合久久99| 亚洲色图av天堂| 国产高清激情床上av| 日本 欧美在线| 日韩欧美国产在线观看| 久久香蕉精品热| 一二三四社区在线视频社区8| 久久午夜亚洲精品久久| 国产在线男女| 18禁黄网站禁片免费观看直播| 国产伦人伦偷精品视频| 精品久久久久久成人av| 女同久久另类99精品国产91| 天堂√8在线中文| 一本一本综合久久| 亚洲午夜理论影院| 国内揄拍国产精品人妻在线| 久久午夜亚洲精品久久| 国产亚洲精品av在线| 亚洲国产精品999在线| 51午夜福利影视在线观看| 91麻豆av在线| 动漫黄色视频在线观看| 最近视频中文字幕2019在线8| 欧美3d第一页| 一区二区三区激情视频| 亚洲av电影在线进入| 国产69精品久久久久777片| 人妻制服诱惑在线中文字幕| 在线免费观看的www视频| 欧美日韩福利视频一区二区| 欧美精品啪啪一区二区三区| www.999成人在线观看| 国产精品美女特级片免费视频播放器| 欧美激情久久久久久爽电影| www.999成人在线观看| 午夜激情福利司机影院| 亚洲av二区三区四区| 久久亚洲精品不卡| av女优亚洲男人天堂| 国产精品人妻久久久久久| 精品久久久久久,| 国产精品久久久久久人妻精品电影| 国产大屁股一区二区在线视频| 日韩欧美 国产精品| 成人特级av手机在线观看| 亚洲精品456在线播放app | 亚洲精品一区av在线观看| а√天堂www在线а√下载| 日本在线视频免费播放| 一进一出抽搐gif免费好疼| 中文字幕熟女人妻在线| 日韩欧美国产在线观看| 在线观看av片永久免费下载| 日本免费a在线| 久久伊人香网站| 日本黄大片高清| 噜噜噜噜噜久久久久久91| 如何舔出高潮| 国产精品久久久久久久电影| 乱人视频在线观看| 亚洲不卡免费看| 永久网站在线| 欧美性感艳星| 亚洲人成网站高清观看| 2021天堂中文幕一二区在线观| 久久6这里有精品| 亚洲天堂国产精品一区在线| 免费看a级黄色片| 嫩草影院入口| 久久久久免费精品人妻一区二区| 国产精品免费一区二区三区在线| 久久精品国产99精品国产亚洲性色| 亚洲色图av天堂| 国产欧美日韩精品亚洲av| 欧美色欧美亚洲另类二区| 婷婷亚洲欧美| 欧美性猛交黑人性爽| 亚洲欧美精品综合久久99| 999久久久精品免费观看国产| 国产免费av片在线观看野外av| 简卡轻食公司| 变态另类成人亚洲欧美熟女| 18禁在线播放成人免费| 九色国产91popny在线| 好男人电影高清在线观看| 亚洲av熟女| 高潮久久久久久久久久久不卡| 亚洲精品一区av在线观看| 久久天躁狠狠躁夜夜2o2o| 丁香六月欧美| 久久精品国产清高在天天线| 久久精品国产99精品国产亚洲性色| 欧美区成人在线视频| 国产一区二区激情短视频| 精品日产1卡2卡| 亚洲第一区二区三区不卡| 免费看美女性在线毛片视频| 看免费av毛片| 亚洲综合色惰| 最新在线观看一区二区三区| 国产一区二区三区视频了| 我要看日韩黄色一级片| 国产精品一区二区免费欧美| 久久国产乱子伦精品免费另类| 中国美女看黄片| 少妇的逼水好多| 少妇丰满av| 精品人妻熟女av久视频| 女同久久另类99精品国产91| 国产黄色小视频在线观看| 国产成人a区在线观看| 国产av在哪里看| 99riav亚洲国产免费| 别揉我奶头 嗯啊视频| 国产毛片a区久久久久| 国产爱豆传媒在线观看| 欧美黄色片欧美黄色片| 天天躁日日操中文字幕| 女人被狂操c到高潮| 国产av麻豆久久久久久久| 国产一区二区三区视频了| 亚州av有码| 中国美女看黄片| 蜜桃久久精品国产亚洲av| 日韩精品中文字幕看吧| 午夜福利欧美成人| 狠狠狠狠99中文字幕| 国产成人av教育| 欧美最黄视频在线播放免费| 99视频精品全部免费 在线| 国产主播在线观看一区二区| 18禁裸乳无遮挡免费网站照片| 国产色婷婷99| 国产精品久久久久久精品电影| 一级黄片播放器| 日本黄色视频三级网站网址| 日韩av在线大香蕉| 国产精品不卡视频一区二区 | 国产伦精品一区二区三区四那| 国内揄拍国产精品人妻在线| 国产成人欧美在线观看| 亚洲中文字幕一区二区三区有码在线看| 91在线精品国自产拍蜜月| 成人一区二区视频在线观看| 国产欧美日韩一区二区三| 免费在线观看成人毛片| 老熟妇仑乱视频hdxx| 夜夜爽天天搞| 欧美激情国产日韩精品一区| 国产精品三级大全| 天堂网av新在线| 亚洲美女视频黄频| 午夜激情福利司机影院| 悠悠久久av| 性插视频无遮挡在线免费观看| 亚洲无线观看免费| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式| 人人妻,人人澡人人爽秒播| 嫁个100分男人电影在线观看| 一个人看视频在线观看www免费| 亚洲av成人av| av在线老鸭窝| 久久久久久久久大av| 久久99热6这里只有精品| 999久久久精品免费观看国产| 久久久久久久精品吃奶| 18+在线观看网站| 一进一出好大好爽视频| 日韩欧美在线二视频| av在线观看视频网站免费| av女优亚洲男人天堂| 在线天堂最新版资源| 岛国在线免费视频观看| 看片在线看免费视频| 精品一区二区免费观看| 免费黄网站久久成人精品 | 国产aⅴ精品一区二区三区波| 成人三级黄色视频| 在线看三级毛片| 老熟妇仑乱视频hdxx| 99riav亚洲国产免费| 日本熟妇午夜| 日本免费a在线| 精品久久久久久久人妻蜜臀av| 国产精品免费一区二区三区在线| 亚洲天堂国产精品一区在线| 久久亚洲精品不卡| 如何舔出高潮| 婷婷色综合大香蕉| 蜜桃亚洲精品一区二区三区| 美女大奶头视频| 午夜激情福利司机影院| 久久欧美精品欧美久久欧美| 国产成人av教育| 两个人的视频大全免费| av女优亚洲男人天堂| 精品国产亚洲在线| 亚洲av免费高清在线观看| 麻豆成人av在线观看| 69人妻影院| 国产熟女xx| 久久久久久国产a免费观看| 人妻丰满熟妇av一区二区三区| 日韩欧美一区二区三区在线观看| 一卡2卡三卡四卡精品乱码亚洲| 自拍偷自拍亚洲精品老妇| 欧美日韩福利视频一区二区| 免费在线观看影片大全网站| 好男人在线观看高清免费视频| 在线天堂最新版资源| 91久久精品国产一区二区成人| 国产午夜精品论理片| a级毛片免费高清观看在线播放| 精品一区二区三区av网在线观看| 国产精品,欧美在线| 欧洲精品卡2卡3卡4卡5卡区| 国产高清激情床上av| 成熟少妇高潮喷水视频| www日本黄色视频网| 国产亚洲欧美98| 日韩人妻高清精品专区| 无人区码免费观看不卡| 91午夜精品亚洲一区二区三区 | 成人欧美大片| 精品久久久久久久久av| 丁香六月欧美| 国产三级黄色录像| 麻豆一二三区av精品| 久久精品影院6| 久久人妻av系列| 国产免费av片在线观看野外av| 国产成人啪精品午夜网站| 少妇被粗大猛烈的视频| 日本三级黄在线观看| 男人舔女人下体高潮全视频| 久久精品国产亚洲av天美| 少妇丰满av| 精品午夜福利视频在线观看一区| www日本黄色视频网| 国产v大片淫在线免费观看| 人人妻人人澡欧美一区二区| 两人在一起打扑克的视频| 亚洲专区中文字幕在线| 一进一出抽搐动态| 在线观看av片永久免费下载| 亚洲成a人片在线一区二区| 欧美一区二区精品小视频在线| 国产亚洲精品久久久com| 波野结衣二区三区在线| 国模一区二区三区四区视频| 日本成人三级电影网站| 国产精品久久视频播放| 久久国产精品影院| 色精品久久人妻99蜜桃| 久久香蕉精品热| 中亚洲国语对白在线视频| 国产精品国产高清国产av| 一夜夜www| 色5月婷婷丁香| 我要看日韩黄色一级片| 好看av亚洲va欧美ⅴa在| 国产精品自产拍在线观看55亚洲| 91狼人影院| 偷拍熟女少妇极品色| 91九色精品人成在线观看| 淫秽高清视频在线观看| 中文字幕av在线有码专区| 国产精品三级大全| 国产精品人妻久久久久久| 国产精品爽爽va在线观看网站| 国产av麻豆久久久久久久| 床上黄色一级片| 中文字幕久久专区| 久久精品综合一区二区三区| netflix在线观看网站| 综合色av麻豆| 日韩国内少妇激情av| 男女之事视频高清在线观看| 亚洲av二区三区四区| 欧美黑人巨大hd| 99久久九九国产精品国产免费| 欧美成人a在线观看| 婷婷精品国产亚洲av| 最近最新免费中文字幕在线| 日韩中文字幕欧美一区二区| 国产人妻一区二区三区在| 亚洲成人久久性| 国产人妻一区二区三区在| 国产极品精品免费视频能看的| 91久久精品国产一区二区成人| 人人妻人人看人人澡| 色噜噜av男人的天堂激情| 久99久视频精品免费| 免费无遮挡裸体视频| 国产国拍精品亚洲av在线观看| 亚洲一区高清亚洲精品| 999久久久精品免费观看国产| 国产精品爽爽va在线观看网站| 久久伊人香网站| 99国产极品粉嫩在线观看| 亚洲精品久久国产高清桃花| 18+在线观看网站| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 国产黄片美女视频| 亚洲最大成人av| 国产 一区 欧美 日韩| 精品人妻1区二区| 久久久久国内视频| 午夜日韩欧美国产| 精品国产三级普通话版| 99久久精品国产亚洲精品| 在线a可以看的网站| 搡老妇女老女人老熟妇| 亚洲人成网站高清观看| 久久精品影院6| 看片在线看免费视频| 综合色av麻豆| 两性午夜刺激爽爽歪歪视频在线观看| 在线观看美女被高潮喷水网站 | 琪琪午夜伦伦电影理论片6080| 人人妻人人看人人澡| 免费看美女性在线毛片视频| 欧美不卡视频在线免费观看| 男人的好看免费观看在线视频| 日本黄大片高清| 亚洲人成电影免费在线| av在线观看视频网站免费| 一级av片app| 国内揄拍国产精品人妻在线| 国产色爽女视频免费观看| 亚洲自偷自拍三级| 国产av麻豆久久久久久久| 在线免费观看不下载黄p国产 | 成年女人永久免费观看视频| 99久久99久久久精品蜜桃| 在线观看美女被高潮喷水网站 | 黄色视频,在线免费观看| 欧美绝顶高潮抽搐喷水| 99riav亚洲国产免费| 国产免费一级a男人的天堂| 国产人妻一区二区三区在| 自拍偷自拍亚洲精品老妇| 婷婷丁香在线五月| 国产免费av片在线观看野外av| 亚洲人成网站高清观看| 午夜福利免费观看在线| 免费看日本二区| 日日摸夜夜添夜夜添av毛片 | 日本熟妇午夜| 国产国拍精品亚洲av在线观看| 国产精品一区二区免费欧美| 国产精品一及| 国产精品一区二区免费欧美| 黄片小视频在线播放| 日韩欧美国产在线观看| 久久99热这里只有精品18| 别揉我奶头 嗯啊视频| 一区福利在线观看| 日本撒尿小便嘘嘘汇集6| 国产av在哪里看| 三级毛片av免费| 国产爱豆传媒在线观看| 好男人在线观看高清免费视频| 国产一级毛片七仙女欲春2| 午夜免费激情av| 日韩中文字幕欧美一区二区| 久9热在线精品视频| 精品一区二区三区人妻视频| 在现免费观看毛片| 亚洲第一区二区三区不卡| 久久久色成人| 成人av一区二区三区在线看| 欧美在线一区亚洲| 脱女人内裤的视频| 亚洲18禁久久av| 观看免费一级毛片| 国产午夜精品久久久久久一区二区三区 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 天天躁日日操中文字幕| 少妇高潮的动态图| 亚洲精品影视一区二区三区av| 国产精品98久久久久久宅男小说| 免费大片18禁| 中文资源天堂在线| 偷拍熟女少妇极品色| 久久久久久久久久成人| av在线天堂中文字幕| 国产熟女xx| 日韩欧美免费精品| 久久久久免费精品人妻一区二区| 亚洲av电影在线进入| 亚洲午夜理论影院| 欧美性猛交黑人性爽| 在线国产一区二区在线| 亚洲国产欧美人成| 精品国内亚洲2022精品成人| www.999成人在线观看| 99热精品在线国产| www.999成人在线观看| 亚洲欧美日韩高清在线视频| 丰满人妻一区二区三区视频av| 国产又黄又爽又无遮挡在线| 久久99热6这里只有精品| 亚洲电影在线观看av| 真人一进一出gif抽搐免费| 国产精品久久久久久久电影| 最后的刺客免费高清国语| av天堂中文字幕网| www.www免费av| 一区二区三区激情视频| 久久午夜亚洲精品久久| 日本成人三级电影网站| 午夜福利在线观看吧| 99视频精品全部免费 在线| 国产亚洲av嫩草精品影院| 精品国产亚洲在线| 久久久色成人| 日韩欧美 国产精品| 热99re8久久精品国产| 午夜福利在线在线| 亚洲欧美激情综合另类| 少妇人妻一区二区三区视频| 狂野欧美白嫩少妇大欣赏| 日韩欧美国产一区二区入口| 老司机午夜十八禁免费视频| 波多野结衣高清无吗| 精品无人区乱码1区二区| 国内少妇人妻偷人精品xxx网站| 国产淫片久久久久久久久 | 免费电影在线观看免费观看| 男女那种视频在线观看| www.999成人在线观看| 99精品在免费线老司机午夜| 午夜福利免费观看在线| 国产熟女xx| 欧美在线一区亚洲| 变态另类成人亚洲欧美熟女| a在线观看视频网站| 1000部很黄的大片| 99精品久久久久人妻精品| 一进一出抽搐gif免费好疼| 99精品久久久久人妻精品| 欧美+亚洲+日韩+国产| 一级黄片播放器| 免费看a级黄色片| 免费看日本二区| 一区福利在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产蜜桃级精品一区二区三区| 国产在视频线在精品| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 欧美黑人欧美精品刺激| 啦啦啦观看免费观看视频高清| 亚洲国产欧美人成| 国产亚洲精品av在线| 偷拍熟女少妇极品色| 国内精品一区二区在线观看| 69av精品久久久久久| 女同久久另类99精品国产91| 国产精品久久视频播放| 婷婷精品国产亚洲av| 中亚洲国语对白在线视频| 亚洲专区国产一区二区| 欧美xxxx黑人xx丫x性爽| 午夜两性在线视频| 久9热在线精品视频| 嫩草影视91久久| 麻豆av噜噜一区二区三区| 精品久久久久久久久亚洲 | 精品人妻一区二区三区麻豆 | 色尼玛亚洲综合影院| 久久精品国产清高在天天线| 亚洲人成网站高清观看| 变态另类丝袜制服| 在线播放无遮挡| 啪啪无遮挡十八禁网站| 国产视频一区二区在线看| av女优亚洲男人天堂| 日本五十路高清| 国产v大片淫在线免费观看| 久久这里只有精品中国| 国产精品野战在线观看| 日韩欧美精品免费久久 | 观看免费一级毛片| 亚洲av熟女| 婷婷精品国产亚洲av在线| 亚洲人成网站高清观看| 亚洲黑人精品在线| av天堂在线播放| 免费观看精品视频网站| 欧美成人一区二区免费高清观看| 亚洲国产精品999在线| 国产在线男女| 国产亚洲精品久久久com| 日本 av在线| 中文资源天堂在线| 国产真实乱freesex| 色吧在线观看|