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

    巴西圓盤劈裂二維及三維數(shù)值模擬研究

    2022-01-14 06:45:52李慶文趙國(guó)彥
    工程科學(xué)學(xué)報(bào) 2022年1期
    關(guān)鍵詞:圓盤圓心端面

    喬 蘭,劉 建,李慶文?,趙國(guó)彥

    1) 北京科技大學(xué)土木與資源工程學(xué)院,北京 100083 2) 中南大學(xué)資源與安全工程學(xué)院,長(zhǎng)沙 410083

    巖石、混凝土和陶瓷等脆性材料的抗拉強(qiáng)度遠(yuǎn)低于其抗壓強(qiáng)度及抗剪強(qiáng)度,巖石的宏觀剪切破壞一般源于其內(nèi)部大量張拉裂紋的相互作用及貫通[1?2],Hoek 和 Martin[3]也曾指出深部巖體的破壞主要是應(yīng)力驅(qū)動(dòng)的張拉破壞,因此巖石類脆性材料的抗拉強(qiáng)度是人們尤其關(guān)注的重要參數(shù).相對(duì)于直接拉伸試驗(yàn),以巴西劈裂試驗(yàn)為代表的抗拉強(qiáng)度間接測(cè)試方法因其操作簡(jiǎn)單方便而在國(guó)內(nèi)外得到廣泛應(yīng)用.國(guó)際巖石力學(xué)學(xué)會(huì)(ISRM)于1978年將巴西劈裂試驗(yàn)推薦為測(cè)定巖石材料抗拉強(qiáng)度的建議方法[4].與此同時(shí),自巴西劈裂試驗(yàn)方法于1943年被提出以來(lái),關(guān)于其有效性的探討一直是國(guó)內(nèi)外學(xué)者的研究熱點(diǎn)[5?7].根據(jù)Griffith強(qiáng)度理論,只有在滿足中心起裂的條件下通過(guò)巴西劈裂試驗(yàn)測(cè)得的巖石抗拉強(qiáng)度才是有效的[5],并且從破壞結(jié)果來(lái)看其破裂面至少應(yīng)表現(xiàn)為與張拉破壞對(duì)應(yīng)的單一徑向破裂面.但是,F(xiàn)airhurst[5]指出對(duì)于壓拉比較低的脆性材料采用小加載角時(shí)起裂點(diǎn)可能由于應(yīng)力集中效應(yīng)而位于巖樣端部.Erarslan等[8?9]、Komurlu 等[10?11]采用 ISRM 推薦的弧形壓模開(kāi)展巴西劈裂試驗(yàn),發(fā)現(xiàn)圓盤試樣發(fā)生多破裂壓碎破壞.Jaeger和Hoskins[12]提出采用弧形均布載荷開(kāi)展巴西劈裂試驗(yàn),通過(guò)增大加載角以減小加載端的應(yīng)力集中程度從而促進(jìn)中心起裂.Erarslan 等[8?9]認(rèn)為加載角為 20°或 30°時(shí)可以獲得可靠的抗拉強(qiáng)度.對(duì)于混凝土巴西圓盤試樣,García等[13]建議加載角應(yīng)大于 12°.Gutiérrez-Moizant等[14]則認(rèn)為開(kāi)展巴西劈裂試驗(yàn)最合適加載角為20°.Bahaaddini等[15]通過(guò)開(kāi)展數(shù)值模擬研究發(fā)現(xiàn),對(duì)于巖石類材料加載角為30°或更高時(shí)才能確保中心起裂.目前,鮮少有學(xué)者探究壓拉比對(duì)試樣起裂破壞模式的影響.此外,國(guó)內(nèi)外學(xué)者比較關(guān)注圓盤試樣受壓直徑上的起裂點(diǎn)位置及破裂發(fā)展方式,較少有學(xué)者探究巴西劈裂三維破裂面的起裂及擴(kuò)展過(guò)程.由于張拉破裂面沿受壓直徑向圓盤內(nèi)部擴(kuò)展,受限于目前的觀測(cè)手段,以高速攝影儀[16]及數(shù)字圖像相關(guān)技術(shù)[17]為代表的表面變形觀測(cè)技術(shù)便不再適用,而數(shù)值模擬方法便顯示出其獨(dú)特的優(yōu)勢(shì).喻勇[6]、Li和 Wong[7]及鄧華鋒等[18]通過(guò)開(kāi)展三維彈性分析探討圓盤內(nèi)部的應(yīng)力分布規(guī)律及起裂點(diǎn)位置,但未模擬三維破裂面的形成過(guò)程.許學(xué)良[19]運(yùn)用PFC3D開(kāi)展巴西劈裂三維模擬,以細(xì)觀裂紋集聚表征張拉破裂面的擴(kuò)展過(guò)程,模擬結(jié)果中破裂面仍然從加載端開(kāi)始起裂,并且起裂階段初始裂紋全部為張拉裂紋.鑒于此,本文借助FLAC3D進(jìn)行連續(xù)介質(zhì)彈塑性分析,一方面開(kāi)展巴西劈裂二維模擬研究,探究圓盤試樣壓拉比及加載角對(duì)試樣起裂破壞模式的影響;另一方面開(kāi)展三維模擬研究,探究圓盤試樣三維破裂面的形成及擴(kuò)展過(guò)程.

    1 巴西圓盤應(yīng)力解析解

    Hondros[20]推導(dǎo)出平面問(wèn)題中弧形均布載荷作用下巴西圓盤內(nèi)部任意一點(diǎn)應(yīng)力分布的解析解,如圖1所示,則加載直徑上的應(yīng)力分布為(以拉應(yīng)力為正):

    圖1 巴西圓盤弧形均布載荷加載示意圖Fig.1 A disk subjected to diametrically distributed loads in the Brazilian test

    式中:σr,σθ,τrθ分別為受壓直徑上的徑向應(yīng)力、周向應(yīng)力及剪應(yīng)力,P為施加的外部載荷,D為巴西圓盤直徑,R為圓盤半徑,t為圓盤厚度,2α為加載角度,r為受壓直徑上一點(diǎn)到圓心的距離.在滿足中心起裂的基礎(chǔ)上,巖石抗拉強(qiáng)度的計(jì)算公式為

    當(dāng)接觸加載角小于15°時(shí),采用近似公式計(jì)算巖石抗拉強(qiáng)度引起的誤差不大于2.3%.此外,ISRM推薦的弧形壓模與巴西圓盤接觸在試樣破壞時(shí)產(chǎn)生的加載角為10°左右[4].

    2 二維數(shù)值模擬

    二維巴西劈裂數(shù)值模擬采用平面應(yīng)變模型,數(shù)值試樣直徑為 50 mm,厚度為 0.5 mm.模型網(wǎng)格采用六面體立方體單元,單元尺寸為0.5 mm,整個(gè)模型共包含7864個(gè)單元,16130個(gè)節(jié)點(diǎn).加載過(guò)程采用位移控制,邊界條件為:固定模型y方向的變形,在模型頂端和底端施加z方向的速度邊界條件,加載速率為每步 1×10?10m.二維加載示意圖如圖2所示.采用彈塑性應(yīng)變軟化本構(gòu)模型及Mohr-Coulomb強(qiáng)度準(zhǔn)則,參考Iddefjord花崗巖[21]的力學(xué)參數(shù),數(shù)值模型參數(shù)如表1所示.當(dāng)模型單元發(fā)生屈服后,單元內(nèi)聚力及抗拉強(qiáng)度隨塑性剪應(yīng)變及塑性拉應(yīng)變累計(jì)逐漸線性軟化;當(dāng)塑性剪應(yīng)變及塑性拉應(yīng)變累計(jì)至臨界值后,單元內(nèi)聚力及抗拉強(qiáng)度弱化至殘余強(qiáng)度并保持不變.為探究圓盤試樣壓拉比及加載角對(duì)試樣起裂破壞模式的影響,分別考慮 6 種壓拉比 (λ)及加載角 (2α):λ=20、18、16、14、12和 10及 2α=4.58°、6.88°、9.18°、11.48°、13.78°和 16.10°.針對(duì)每種壓拉比,保持單軸抗壓強(qiáng)度為200 MPa不變.數(shù)值模擬結(jié)果如表2所示.

    圖2 二維巴西劈裂加載示意圖Fig.2 Loading setup of the 2D numerical Brazilian test

    表1 數(shù)值模型參數(shù)Table 1 Material properties of the numerical model

    表2 不同壓拉比及加載角巴西試樣的破壞模式Table 2 Failure modes of the numerical Brazilian disks with different compression?tension ratios (λ) and contact loading angles (2α)

    從表2橫向分析,接觸加載角越大,巴西圓盤試樣越容易發(fā)生中心起裂,因?yàn)樵龃蠹虞d角可以減小圓盤試樣頂?shù)锥说膽?yīng)力集中程度.以λ=14為例,當(dāng) 2α≤ 9.18°時(shí),巴西試樣發(fā)生端部起裂;當(dāng) 2α≥ 11.48°時(shí),巴西試樣發(fā)生中心起裂.沿表2縱向分析,壓拉比越大,巴西圓盤試樣越容易發(fā)生中心起裂.以 2α=11.48°為例,當(dāng)λ≥ 14 時(shí),試樣發(fā)生中心起裂;而當(dāng)λ≤ 12時(shí),圓盤試樣發(fā)生端部起裂.較大的壓拉比,一方面表示巖樣的抗拉強(qiáng)度相對(duì)較低,因此試樣更容易發(fā)生張拉破壞;另一方面則表示巖樣的抗壓強(qiáng)度相對(duì)較高,抵抗剪切破壞的能力強(qiáng).因此,巴西圓盤試樣的起裂破壞模式受內(nèi)部因素與外部因素共同控制,總體上,圓盤試樣壓拉比越大,外部加載角越大,試樣破壞越容易產(chǎn)生中心起裂,反之則越容易發(fā)生端部起裂.

    圖3(a)~(b)為均勻巴西圓盤試樣的加載破壞過(guò)程,包括峰值載荷(Fpeak)前及峰后階段,接觸加載角為 9.18°,壓拉比分別為λ=20、10.當(dāng)λ=20時(shí),圓盤試樣發(fā)生中心起裂,隨外載進(jìn)一步增加,徑向張拉裂紋自中心起裂點(diǎn)逐漸向兩端擴(kuò)展,并在加載點(diǎn)附近趨于停止,這與Erarslan等[8?9]通過(guò)實(shí)驗(yàn)觀察到的裂紋擴(kuò)展模式相似,陳沙等[22]也通過(guò)數(shù)值模擬獲得過(guò)類似的裂紋擴(kuò)展過(guò)程.而當(dāng)λ=10時(shí),兩條起裂裂紋從圓盤受力區(qū)域邊緣向圓盤內(nèi)部擴(kuò)展,初始起裂裂紋主要發(fā)生剪切破壞,向試樣中心擴(kuò)展后則轉(zhuǎn)變?yōu)閺埨茐模罱K形成兩條貫通的劈裂縫,十分明顯這是由于各向同性圓盤試樣及外載關(guān)于受壓直徑的對(duì)稱性引起.圖3(c)~(e)為非均勻巴西圓盤試樣的加載破壞過(guò)程,加載角仍為 9.18°,壓拉比分別為λ=20、10、8.運(yùn)用劉建等[23]提出的巖石介質(zhì)細(xì)觀物理力學(xué)參數(shù)統(tǒng)計(jì)分布模型,采用對(duì)數(shù)正太分布描述細(xì)觀單元彈性模量及單軸抗壓強(qiáng)度的空間變異性[24?25],試樣均質(zhì)度參數(shù)設(shè)為0.3.當(dāng)λ=20時(shí),數(shù)值試樣同樣產(chǎn)生中心起裂,受模型非均質(zhì)性的影響起裂點(diǎn)偏離圓盤正中心,同時(shí)在圓盤內(nèi)部受壓直徑周圍分布有離散的張拉破裂單元.這是由于一方面材料的非均質(zhì)性改變了圓盤內(nèi)部的應(yīng)力分布,引起局部應(yīng)力集中;另一方面,可能是由于軟弱介質(zhì)單元的過(guò)早破壞引起.而當(dāng)λ=10、8時(shí),與均勻條件下的模擬結(jié)果明顯不同的是,圓盤試樣發(fā)生端部起裂后最終僅形成一條劈裂縫;在劈裂縫形成之前,試樣兩端出現(xiàn)剪切破壞單元,而試樣內(nèi)部基本無(wú)破壞單元出現(xiàn),壓拉比為8的數(shù)值試樣更明顯,這與楊志鵬等[26]關(guān)于頁(yè)巖巴西劈裂的聲發(fā)射定位結(jié)果類似.值得指出的是,無(wú)論是均勻試樣還是非均勻試樣,端部起裂均由剪切破壞引起,而劈裂裂紋進(jìn)一步擴(kuò)展則由張拉破壞驅(qū)動(dòng).

    圖3 不同壓拉比時(shí)巴西圓盤試樣的破壞過(guò)程.(a)均勻試樣 λ=20;(b) 均勻試樣 λ=10;(c)非均勻試樣 λ=20;(d)非均勻試樣 λ=10;(e)非均勻試樣λ=8Fig.3 Failure processes of the numerical disks with different compression–tension ratios (λ) under the loading angle of 9.18°: (a) homogeneous disk λ =20; (b) homogeneous disk λ = 10; (c) heterogeneous disk λ = 20; (d) heterogeneous disk λ = 10; (e) heterogeneous disk λ = 8

    3 三維數(shù)值模擬

    三維巴西圓盤數(shù)值模型直徑為50 mm,厚度為25 mm;采用六面體網(wǎng)格進(jìn)行劃分,對(duì)受壓直徑附近的網(wǎng)格進(jìn)行適當(dāng)加密,模型共包含69125個(gè)單元,74360個(gè)節(jié)點(diǎn),模型網(wǎng)格如圖4所示.采用應(yīng)力邊界條件,在受壓區(qū)域內(nèi)沿模型表面施加法向正應(yīng)力,加載速率為每步500 Pa.三維數(shù)值模擬仍采用表1中的力學(xué)參數(shù);為模擬不同接觸加載角下巴西圓盤三維破裂面的起裂及擴(kuò)展過(guò)程,考慮以下 8種加載角度:2α=29.7°、25.7°、21.7°、17.9°、14.4°、10.9°、7.6°、4.4°.

    圖4 三維巴西圓盤模型網(wǎng)格示意圖Fig.4 Grid model of the 3D numerical Brazilian disk

    圖5 為加載角 2α=29.7°、25.7°、14.4°和 4.4°時(shí)巴西圓盤試樣三維破裂面的起裂及擴(kuò)展過(guò)程,由于篇幅限制,其他4種加載角下的數(shù)值模擬結(jié)果沒(méi)有給出,但破裂面發(fā)展規(guī)律基本類似.當(dāng)2α=29.7°時(shí),初始起裂點(diǎn)位于圓盤端面中心,隨外部載荷增加,破裂面以弧形邊界向試樣內(nèi)部擴(kuò)展,在圓盤厚度方向首先貫通后繼續(xù)沿受壓直徑向兩端擴(kuò)展并在加載點(diǎn)附近趨于停止.值得注意的是,貫通后的破裂面沿受壓直徑方向的擴(kuò)展邊界也并非直線,這與平面應(yīng)變假設(shè)下貫通破裂面的直線擴(kuò)展情形相矛盾.當(dāng) 2α=25.7°、14.4°時(shí),起裂點(diǎn)仍位于圓盤端面受壓直徑上,但卻偏離端面中心,同時(shí)出現(xiàn)4個(gè)對(duì)稱的起裂點(diǎn);加載角越小,起裂點(diǎn)越靠近上下加載端,并且破裂面擴(kuò)展邊界線仍為弧形.當(dāng)2α=4.4°時(shí),起裂點(diǎn)位于圓盤端面加載區(qū)域邊緣,即破裂沒(méi)有在受壓中心面上發(fā)展,這與二維模擬中端部起裂的情況類似.事實(shí)上,8種接觸加載角下圓盤試樣的起裂點(diǎn)均位于圓盤端面,而不是圓盤內(nèi)部中心.喻勇[8]及鄧華鋒等[20]通過(guò)對(duì)線荷載作用下的巴西圓盤試樣進(jìn)行三維彈性分析指出,圓盤橫截面上的應(yīng)力值沿試樣厚度方向是有變化的,越靠近兩端面,水平拉應(yīng)力越大,因此試樣破壞首先從端面開(kāi)始.

    圖5 不同加載角時(shí)巴西圓盤試樣三維破裂面的起裂及擴(kuò)展過(guò)程.(a)2α=29.7°;(b)2α=25.7°;(c)2α=14.4°;(d)2α=4.4°Fig.5 Fracture initiation and propagation processes of the 3D Brazilian disks with different contact loading angles: (a) 2α = 29.7°; (b) 2α = 25.7°; (c) 2α =14.4°; (d) 2α = 4.4°

    4 解析驗(yàn)證及應(yīng)力分析

    4.1 二維模擬結(jié)果分析

    圖6為λ=20時(shí)不同加載角下二維巴西圓盤受壓直徑上切向應(yīng)力與徑向應(yīng)力的數(shù)值模擬結(jié)果與理論計(jì)算結(jié)果的對(duì)比示意圖.由圖可知,除在加載端附近數(shù)值模擬結(jié)果與解析解存在較小誤差外,受壓直徑上大部分區(qū)域具有較好的吻合性,說(shuō)明了二維數(shù)值模擬的有效性.需要指出的是,由于二維數(shù)值模型采用正六面體單元逼近圓盤試樣進(jìn)行網(wǎng)格劃分,因此在模型頂?shù)锥藶橹本€邊界而非弧形邊界,并且施加的速度邊界條件與受壓直徑平行而非沿圓盤半徑指向圓盤中心,如圖3所示,這是造成上述誤差的原因.在加載端附近,徑向應(yīng)力與切向應(yīng)力均為壓應(yīng)力,靠近圓盤中心,兩者快速減小,并且切向應(yīng)力由壓應(yīng)力逐漸轉(zhuǎn)變?yōu)槔瓚?yīng)力.隨著接觸加載角增大,受壓直徑上張拉區(qū)域逐漸減小,加載端附近壓縮區(qū)域逐漸增大,徑向與切向壓應(yīng)力也逐漸減小,這說(shuō)明端部應(yīng)力集中程度在減弱.

    圖6 不同加載角時(shí)二維巴西圓盤受壓直徑上應(yīng)力分布的數(shù)值模擬結(jié)果與理論計(jì)算結(jié)果對(duì)比示意圖.(a) 2α=4.58°;(b)2α=6.88°;(c)2α=9.18°;(d)2α=11.48°;(e)2α=13.78°;(f)2α=16.10°Fig.6 Comparison of normalized stresses along the compressed diameter between the numerical results and Hondros’ solutions with different contact loading angles: (a) 2α = 4.58°; (b) 2α = 6.88°; (c) 2α = 9.18°; (d) 2α = 11.48°; (e) 2α = 13.78°; (f) 2α = 16.10°

    表3為不同加載角及壓拉比下二維巴西圓盤試樣抗拉強(qiáng)度的計(jì)算結(jié)果.由表3可知,對(duì)于發(fā)生中心起裂的圓盤試樣,通過(guò)公式(4)計(jì)算的抗拉強(qiáng)度與圓盤試樣的真實(shí)抗拉強(qiáng)度比較接近,相對(duì)誤差在3%以內(nèi)(括號(hào)內(nèi)數(shù)據(jù)為相對(duì)誤差);而對(duì)于發(fā)生端部起裂的圓盤試樣來(lái)說(shuō),其計(jì)算抗拉強(qiáng)度明顯低于試樣的真實(shí)抗拉強(qiáng)度,并且壓拉比及加載角越小,兩者相對(duì)誤差也越大.上述分析表明只有在滿足中心起裂的基礎(chǔ)上依據(jù)峰值荷載計(jì)算出的抗拉強(qiáng)度才是有效的,而端部起裂條件下巴西劈裂實(shí)驗(yàn)會(huì)低估巖石的抗拉強(qiáng)度.

    表3 不同加載角及壓拉比下二維巴西圓盤試樣抗拉強(qiáng)度計(jì)算結(jié)果Table 3 Calculated tensile strength of 2D numerical Brazilian disks with different contact loading angles and compression –tension ratios

    4.2 三維模擬結(jié)果分析

    圖7 為接觸加載角 2α=4.4°、14.4°、25.7°、29.7°時(shí)三維巴西圓盤試樣端面(Y=0.025)及中心剖面(Y=0.0125)受壓直徑上徑向應(yīng)力與切向應(yīng)力的分布示意圖,圖中同時(shí)給出了Hondros解析解的計(jì)算結(jié)果.由圖可知,4種加載角下中心剖面受壓直徑上徑向與切向應(yīng)力的數(shù)值模擬結(jié)果與Hondros解析解的計(jì)算結(jié)果十分吻合,而端面受壓直徑上應(yīng)力分布的數(shù)值模擬結(jié)果與Hondros解析結(jié)果存在一定誤差,這是由于三維圓盤試樣及外載關(guān)于中心剖面對(duì)稱,因此中心剖面更貼近平面應(yīng)變條件.圓盤端面受壓直徑上切向拉應(yīng)力大于中心剖面受壓直徑上的切向拉應(yīng)力,2α=25.7°、29.7°時(shí)尤其明顯.圖8為8種加載角下三維巴西圓盤中心剖面及端面受壓直徑上切向應(yīng)力分布示意圖,中心剖面受壓直徑上最大切向拉應(yīng)力始終位于剖面中心;而在端面受壓直徑上,當(dāng)加載角比較小時(shí),最大切向拉應(yīng)力靠近圓盤頂?shù)锥耍S著加載角增大,最大切向拉應(yīng)力逐漸向端面圓心移動(dòng),當(dāng)2α=29.7°時(shí),最大切向拉應(yīng)力出現(xiàn)在端面圓心處.

    圖7 不同加載角下三維巴西圓盤受壓直徑上應(yīng)力分布的數(shù)值模擬結(jié)果與理論計(jì)算結(jié)果對(duì)比示意圖.(a)2α=4.4°;(b)2α=14.4°;(c)2α=25.7°;(d)2α=29.7°Fig.7 Comparison of normalized stresses along the compressed diameter between the numerical results and Hondros’ solutions with different contact loading angles: (a) 2α = 4.4°; (b) 2α = 14.4°; (c) 2α = 25.7°; (d) 2α = 29.7°

    圖8 不同加載角下三維巴西圓盤試樣中心剖面及端面受壓直徑上切向應(yīng)力分布示意圖.(a)中心剖面(Y=0.0125)切向應(yīng)力;(b)端面(Y=0.025)切向應(yīng)力Fig.8 Normalized tangential stresses along the compressed diameter of the surface and middle section of 3D disks with different contact loading angles:(a) middle section (Y=0.0125); (b) surface (Y=0.025)

    圖9 為 2α=29.7°、25.7°、14.4°、7.6°、4.4°時(shí)巴西圓盤軸向受壓平面上(X=0)切向拉應(yīng)力分布云圖.十分明顯,圓盤內(nèi)部切向拉應(yīng)力明顯小于圓盤外部切向拉應(yīng)力,這是由于各質(zhì)點(diǎn)受到周圍質(zhì)點(diǎn)的約束程度不一樣引起的.從圓盤中心剖面(Y=0.0125)越靠近圓盤端面,質(zhì)點(diǎn)受到的水平約束也越小,水平張拉變形越大,因此張拉應(yīng)力就相對(duì)越大,所以三維圓盤試樣的最先起裂點(diǎn)出現(xiàn)在圓盤端面.而在圓盤端面上,除 2α=4.4°外,其他 7 種加載角下,起裂點(diǎn)位置就是端面受壓直徑上最大切向拉應(yīng)力出現(xiàn)的點(diǎn).本文運(yùn)用考慮張拉截?cái)嗟腗ohr-Coulomb強(qiáng)度準(zhǔn)則進(jìn)行模擬分析,當(dāng)圓盤端面受壓直徑上某一點(diǎn)的最大切向拉應(yīng)力大于試樣抗拉強(qiáng)度時(shí)就會(huì)發(fā)生張拉破壞,其本質(zhì)為最大拉應(yīng)力強(qiáng)度準(zhǔn)則.Li和Wong[7]運(yùn)用Stacey[27]張拉應(yīng)變準(zhǔn)則對(duì)三維巴西圓盤起裂點(diǎn)位置進(jìn)行分析,當(dāng)圓盤受壓直徑上某一點(diǎn)的最大張拉應(yīng)變超過(guò)某一限值時(shí)就會(huì)發(fā)生開(kāi)裂,其本質(zhì)為最大拉應(yīng)變強(qiáng)度準(zhǔn)則.喻勇[6]指出三維巴西圓盤的應(yīng)力應(yīng)變分布規(guī)律與材料泊松比有關(guān),而與材料彈性模量無(wú)關(guān).Li和Wong[7]通過(guò)數(shù)值模擬研究發(fā)現(xiàn),三維圓盤端面圓心處的拉應(yīng)變隨泊松比線性增大,而端面受壓直徑上的最大拉應(yīng)變隨泊松比呈現(xiàn)非線性增大趨勢(shì),但是其未探討受壓直徑上最大張拉應(yīng)變位置的變化規(guī)律.圖10 為泊松比μ=0.15、0.20、0.25、0.30、0.35時(shí)最大切向拉應(yīng)力點(diǎn)及最大切向拉應(yīng)變點(diǎn)位置隨接觸加載角的變化情況.由圖可知,隨泊松比減小,受壓直徑上最大切向拉應(yīng)力點(diǎn)及最大切向拉應(yīng)變點(diǎn)位置均會(huì)向端面圓心移動(dòng).上述模擬結(jié)果表明,對(duì)于三維巴西圓盤,其端面起裂點(diǎn)位置還受圓盤材料的泊松比影響,較小的泊松比會(huì)促進(jìn)中心起裂.

    圖9 不同加載角時(shí)三維巴西圓盤試樣軸向受壓平面上(X=0)切向應(yīng)力分布云圖.(a)2α=29.7°;(b)2α=25.7°;(c)2α=14.4°;(d)2α=7.6°;(e)2α=4.4°Fig.9 Contour plots of normalized tangential stresses (σxx) on the compressed middle section plane (X = 0) of the 3D disks with different contact loading angles: (a) 2α = 29.7°; (b) 2α = 25.7°; (c) 2α = 14.4°; (d) 2α = 7.6°; (e) 2α = 4.4°

    圖10 不同泊松比時(shí)最大切向拉應(yīng)力點(diǎn)及最大切向拉應(yīng)變點(diǎn)位置隨接觸加載角的變化情況.(a)最大切向拉應(yīng)力點(diǎn);(b)最大切向拉應(yīng)變點(diǎn)Fig.10 Change in the position of the maximum tangential tensile stress point and strain point with the contact loading angles under different Poisson’s ratios: (a) the maximum tangential tensile stress point; (b) the maximum tangential tensile strain point

    表4為9種加載角下三維巴西圓盤試樣抗拉強(qiáng)度計(jì)算結(jié)果,除上述8種加載角外增加2α=33.7°這一情況,表中Distance表示起裂點(diǎn)至端面圓心的距離.當(dāng) 2α=4.4°時(shí),圓盤試樣發(fā)生端部起裂,峰值載荷為5427 N,計(jì)算抗拉強(qiáng)度最小,而相對(duì)誤差最大為65.5%;隨著加載角增大,起裂點(diǎn)逐漸向圓盤端面圓心移動(dòng),峰值載荷逐漸增大,計(jì)算抗拉強(qiáng)度逐漸接近試樣的真實(shí)抗拉強(qiáng)度.但是,當(dāng)2α=29.7°及33.7°時(shí),圓盤試樣發(fā)生端面中心起裂,此時(shí)計(jì)算抗拉強(qiáng)度分別為 7.17 MPa、7.16 MPa,仍低于試樣的真實(shí)抗拉強(qiáng)度,相對(duì)誤差均為10.4%.上述計(jì)算及分析結(jié)果表明,即便在保證中心起裂的條件下,由于三維效應(yīng)巴西劈裂試驗(yàn)可能仍會(huì)低估巖石試樣的抗拉強(qiáng)度.但從工程角度而言,其結(jié)果是偏于保守的,有利于工程穩(wěn)定.

    表4 不同加載角下三維巴西圓盤試樣抗拉強(qiáng)度計(jì)算結(jié)果Table 4 Calculated tensile strength of the 3D numerical Brazilian disks with different loading angles

    5 討論

    本文運(yùn)用Mohr-Coulomb強(qiáng)度準(zhǔn)則進(jìn)行模擬,由于圓盤破裂面上發(fā)生張拉破壞,因此其實(shí)質(zhì)為最大拉應(yīng)力強(qiáng)度準(zhǔn)則.當(dāng)圓盤受壓直徑上最大切向拉應(yīng)力位于圓盤圓心時(shí),隨外載增加圓盤最終就會(huì)出現(xiàn)中心起裂.由于直接對(duì)公式(2)關(guān)于r求極值十分困難,因此采用數(shù)值枚舉法分析公式(2)在不同加載角下的最大值的取值位置,r/R、α取值間距分別為 0.001、0.1,結(jié)果當(dāng)α≤ 30°時(shí)公式(2)的最大值均在r/R=0即圓心處取得.這是二維數(shù)值模擬中試樣出現(xiàn)中心起裂的內(nèi)在原因.

    Griffith強(qiáng)度準(zhǔn)則的實(shí)質(zhì)為當(dāng)一點(diǎn)的Griffith等效拉應(yīng)力等于其抗拉強(qiáng)度時(shí)就會(huì)發(fā)生破壞.對(duì)于線荷載作用下的平面巴西圓盤,其受壓直徑上的應(yīng)力分布為

    式中:r1,r2分別為受壓直徑上一點(diǎn)到上下加載點(diǎn)處的距離,滿足r1+r2=2R.將上述兩式代入Griffith強(qiáng)度準(zhǔn)則,則Griffith等效拉應(yīng)力為

    在圓心處,σG取最小值,為 2P/πDt;而在加載點(diǎn)處,σG為無(wú)窮大.因此,根據(jù)Griffith強(qiáng)度準(zhǔn)則加載點(diǎn)處最先起裂.喻勇[6]通過(guò)數(shù)值模擬研究同樣發(fā)現(xiàn)數(shù)值試樣圓心處等效拉應(yīng)力最小,而上下加載點(diǎn)處等效拉應(yīng)力最大.如果假設(shè)圓心處首先起裂,將r1=r2=R代入上式,則

    這就是通常的抗拉強(qiáng)度的計(jì)算公式.對(duì)于弧形均布荷載作用下的平面巴西圓盤,受壓直徑上的Griffith等效拉應(yīng)力為

    同樣采用數(shù)值枚舉法分析其在不同加載角下受壓直徑上最大等效拉應(yīng)力點(diǎn)的位置變化特征.結(jié)果如圖11所示,圖中同時(shí)給出了三維巴西圓盤端面受壓直徑上最大等效拉應(yīng)力位置隨加載角度的變化情況.對(duì)于公式(9),當(dāng) 2α≥18.2°時(shí)最大等效拉應(yīng)力點(diǎn)位于圓心處.而對(duì)于三維巴西圓盤來(lái)說(shuō),以泊松比μ=0.20為例,由于三維效應(yīng),當(dāng)2α≥29.7°時(shí),最大等效拉應(yīng)力點(diǎn)才移動(dòng)到端面圓心處.當(dāng)發(fā)生中心破壞時(shí),將r/R=0代入上式,則

    圖11 最大Griffith等效拉應(yīng)力位置隨加載角度的變化情況Fig.11 Change in the position of the maximum Griffith equivalent tensile stresses with the contact loading angles

    上式即為根據(jù)Griffith強(qiáng)度準(zhǔn)則推出的弧形均布荷載作用下平面巴西圓盤的抗拉強(qiáng)度計(jì)算公式.對(duì)于近似公式,當(dāng) 2α=20°相對(duì)誤差為 4.0%,而當(dāng) 2α=30°相對(duì)誤差達(dá)到 9.0%.需要指出的是,公式(4)是根據(jù)最大拉應(yīng)力強(qiáng)度準(zhǔn)則導(dǎo)出的.隨著加載角度減小,公式(4)及(10)逐漸逼近公式(8).由于一般巖石的泊松比介于0.15~0.25,通過(guò)三維巴西圓盤端面受壓直徑上最大切向拉應(yīng)力及最大Griffith等效拉應(yīng)力位置隨加載角度的變化趨勢(shì)(圖10(a)及圖11),本文建議加載角度至少為 30°,因?yàn)榇藭r(shí)無(wú)論在最大拉應(yīng)力強(qiáng)度準(zhǔn)則還是在Griffith強(qiáng)度準(zhǔn)則下,端面初始破壞點(diǎn)均位于圓心處.

    6 結(jié)論

    (1)二維數(shù)值模擬結(jié)果表明,接觸加載角及壓拉比越大,巴西圓盤試樣越容易發(fā)生中心起裂;端部起裂由剪切破壞引起,而劈裂裂紋進(jìn)一步擴(kuò)展則由張拉破壞驅(qū)動(dòng);端部起裂條件下巴西劈裂實(shí)驗(yàn)會(huì)低估巖石的抗拉強(qiáng)度.

    (2)三維數(shù)值模擬結(jié)果表明,初始起裂點(diǎn)位于三維圓盤端面,隨加載角增大其沿受壓直徑逐漸向端面圓心移動(dòng);當(dāng)圓盤發(fā)生端面中心起裂時(shí),三維破裂面以弧形邊界向試樣內(nèi)部發(fā)散擴(kuò)展;對(duì)于三維巴西圓盤,其端面起裂點(diǎn)位置還受圓盤材料的泊松比影響,較小的泊松比會(huì)促進(jìn)中心起裂.

    (3)無(wú)論圓盤試樣發(fā)生中心起裂還是端部起裂,由于三維效應(yīng)巴西劈裂試驗(yàn)可能都會(huì)低估巖石的抗拉強(qiáng)度.

    猜你喜歡
    圓盤圓心端面
    二次曲線的一個(gè)類似圓心的性質(zhì)
    KDF3E成型機(jī)濾棒端面觸頭的原因及排除方法
    圓盤鋸刀頭的一種改進(jìn)工藝
    石材(2020年6期)2020-08-24 08:27:00
    以圓周上一點(diǎn)為圓心作圓的圖的性質(zhì)及應(yīng)用
    銅基合金襯套端面鍍鉻質(zhì)量的改善
    單位圓盤上全純映照模的精細(xì)Schwarz引理
    奇怪的大圓盤
    老舊端面磨齒機(jī)故障處理
    基于Profibus-DP的圓盤澆鑄控制系統(tǒng)的應(yīng)用
    貫通軸端面花鍵拉床的設(shè)計(jì)
    欧美区成人在线视频| 少妇人妻一区二区三区视频| 人妻制服诱惑在线中文字幕| 国产免费又黄又爽又色| 青春草国产在线视频| 亚洲高清免费不卡视频| 国产精品人妻久久久影院| av在线蜜桃| 不卡视频在线观看欧美| 大又大粗又爽又黄少妇毛片口| 国产成人aa在线观看| 亚洲av一区综合| 老司机影院成人| 国产精品福利在线免费观看| 欧美人与善性xxx| 欧美高清成人免费视频www| 国产精品国产三级国产专区5o| 成年人午夜在线观看视频| 日本爱情动作片www.在线观看| 寂寞人妻少妇视频99o| 久久久欧美国产精品| 秋霞伦理黄片| 日日啪夜夜爽| 在线免费十八禁| 午夜亚洲福利在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最近最新中文字幕免费大全7| 国产女主播在线喷水免费视频网站| 特级一级黄色大片| 欧美一级a爱片免费观看看| 午夜精品国产一区二区电影 | 少妇猛男粗大的猛烈进出视频 | 亚洲美女搞黄在线观看| 69人妻影院| 日韩av在线免费看完整版不卡| 蜜臀久久99精品久久宅男| 五月伊人婷婷丁香| 国产av不卡久久| a级毛片免费高清观看在线播放| 色视频在线一区二区三区| 国产亚洲精品久久久com| 青春草国产在线视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲人成网站在线观看播放| 成人亚洲欧美一区二区av| freevideosex欧美| 精品人妻视频免费看| 欧美亚洲 丝袜 人妻 在线| 久久久久精品久久久久真实原创| 精品久久久久久电影网| 亚洲av免费在线观看| 黑人高潮一二区| av又黄又爽大尺度在线免费看| 青春草亚洲视频在线观看| 国产在视频线精品| 99久国产av精品国产电影| 国产精品99久久99久久久不卡 | 亚洲欧美日韩无卡精品| 91午夜精品亚洲一区二区三区| 欧美亚洲 丝袜 人妻 在线| 天天躁夜夜躁狠狠久久av| 免费看光身美女| 国产探花在线观看一区二区| 欧美丝袜亚洲另类| 51国产日韩欧美| 亚洲精品国产av成人精品| 午夜亚洲福利在线播放| 有码 亚洲区| 日本与韩国留学比较| 亚洲欧美日韩东京热| a级一级毛片免费在线观看| 成人美女网站在线观看视频| 国产综合懂色| 国产亚洲一区二区精品| 热99国产精品久久久久久7| 久久ye,这里只有精品| 深夜a级毛片| 大话2 男鬼变身卡| 内射极品少妇av片p| 夜夜爽夜夜爽视频| 国产91av在线免费观看| 亚洲欧美中文字幕日韩二区| 一区二区三区四区激情视频| 亚洲av福利一区| 国产有黄有色有爽视频| 国产精品久久久久久精品古装| 建设人人有责人人尽责人人享有的 | 一级毛片电影观看| 夫妻性生交免费视频一级片| 久久久欧美国产精品| 赤兔流量卡办理| 97在线人人人人妻| 中文字幕人妻熟人妻熟丝袜美| 中文字幕免费在线视频6| 熟女人妻精品中文字幕| 最近中文字幕高清免费大全6| 日本爱情动作片www.在线观看| 欧美成人a在线观看| 少妇裸体淫交视频免费看高清| 精品久久国产蜜桃| 五月玫瑰六月丁香| 神马国产精品三级电影在线观看| 亚洲精品视频女| 91aial.com中文字幕在线观看| 日产精品乱码卡一卡2卡三| 国产精品一区二区三区四区免费观看| 亚洲av中文字字幕乱码综合| 国产欧美日韩一区二区三区在线 | 欧美精品国产亚洲| 亚洲精品成人久久久久久| 亚洲精品视频女| av.在线天堂| 国产精品国产av在线观看| 中文在线观看免费www的网站| 国产精品一区www在线观看| 91久久精品国产一区二区三区| www.色视频.com| 国产精品成人在线| 国产成人freesex在线| 男女无遮挡免费网站观看| 欧美潮喷喷水| 亚洲天堂国产精品一区在线| av线在线观看网站| 99久久人妻综合| 99久久人妻综合| 日日啪夜夜爽| 老女人水多毛片| 免费黄网站久久成人精品| 日本三级黄在线观看| 亚洲久久久久久中文字幕| 国产精品伦人一区二区| 久久精品国产亚洲网站| 又爽又黄无遮挡网站| 免费观看的影片在线观看| 在线免费观看不下载黄p国产| 亚洲最大成人手机在线| 国产亚洲精品久久久com| 春色校园在线视频观看| 久久97久久精品| 久久久精品欧美日韩精品| 人妻少妇偷人精品九色| 欧美97在线视频| 观看免费一级毛片| 欧美性感艳星| 少妇猛男粗大的猛烈进出视频 | 成人二区视频| 蜜桃亚洲精品一区二区三区| 欧美3d第一页| 在线亚洲精品国产二区图片欧美 | 极品少妇高潮喷水抽搐| 在线观看美女被高潮喷水网站| 99久久中文字幕三级久久日本| 国产亚洲5aaaaa淫片| 精品人妻偷拍中文字幕| 成人高潮视频无遮挡免费网站| 国产淫语在线视频| 国产色爽女视频免费观看| 中文字幕制服av| a级毛色黄片| 国产成人aa在线观看| 好男人在线观看高清免费视频| 六月丁香七月| 免费电影在线观看免费观看| 91久久精品电影网| 国产高潮美女av| 亚洲精品乱码久久久久久按摩| 亚洲一区二区三区欧美精品 | av在线播放精品| 国精品久久久久久国模美| 天美传媒精品一区二区| 久久精品久久精品一区二区三区| 国产av国产精品国产| 国模一区二区三区四区视频| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品成人综合色| 国产av国产精品国产| 亚洲国产最新在线播放| 亚洲最大成人手机在线| 麻豆成人午夜福利视频| 一级黄片播放器| 亚洲成人久久爱视频| 国产人妻一区二区三区在| 少妇 在线观看| 99热全是精品| 久久99热6这里只有精品| 热re99久久精品国产66热6| 免费看av在线观看网站| 特大巨黑吊av在线直播| 国产爽快片一区二区三区| 99热这里只有精品一区| 亚洲图色成人| 国产精品蜜桃在线观看| 国产 一区 欧美 日韩| 国产一区有黄有色的免费视频| 熟女av电影| 九九在线视频观看精品| 麻豆久久精品国产亚洲av| 深夜a级毛片| 99热国产这里只有精品6| 国产亚洲91精品色在线| 免费大片18禁| videossex国产| 插阴视频在线观看视频| 久久精品国产亚洲网站| 美女视频免费永久观看网站| 日韩免费高清中文字幕av| 日韩av免费高清视频| 久久热精品热| 夫妻性生交免费视频一级片| 国产欧美日韩精品一区二区| 成人鲁丝片一二三区免费| 久久午夜福利片| 九色成人免费人妻av| 亚洲自拍偷在线| av天堂中文字幕网| 成人欧美大片| 免费黄频网站在线观看国产| 久久久久国产精品人妻一区二区| 午夜亚洲福利在线播放| 中国三级夫妇交换| 只有这里有精品99| 内射极品少妇av片p| 男女那种视频在线观看| 日韩中字成人| 日本-黄色视频高清免费观看| 欧美精品人与动牲交sv欧美| 久久国产乱子免费精品| 亚州av有码| av卡一久久| 国产精品99久久99久久久不卡 | 亚洲国产日韩一区二区| 久久久午夜欧美精品| 国产精品国产三级国产av玫瑰| 亚洲自偷自拍三级| 听说在线观看完整版免费高清| 久久综合国产亚洲精品| 久久韩国三级中文字幕| 日韩亚洲欧美综合| 在线免费观看不下载黄p国产| 国产精品蜜桃在线观看| 久久精品国产鲁丝片午夜精品| 国产精品麻豆人妻色哟哟久久| 午夜精品一区二区三区免费看| 十八禁网站网址无遮挡 | 有码 亚洲区| 亚洲国产精品成人综合色| 国产精品99久久99久久久不卡 | 国产精品国产三级国产av玫瑰| 狂野欧美激情性xxxx在线观看| 少妇裸体淫交视频免费看高清| 婷婷色综合www| 国产视频内射| 青春草视频在线免费观看| 国产欧美日韩一区二区三区在线 | 亚洲精品乱码久久久v下载方式| 国产伦精品一区二区三区四那| 夜夜看夜夜爽夜夜摸| 亚洲精品aⅴ在线观看| 黄片无遮挡物在线观看| 国产国拍精品亚洲av在线观看| 色播亚洲综合网| 亚洲国产成人一精品久久久| 波多野结衣巨乳人妻| av在线亚洲专区| 亚洲av二区三区四区| 精品熟女少妇av免费看| 一区二区三区精品91| 99热国产这里只有精品6| 97超视频在线观看视频| 日韩强制内射视频| 天天躁日日操中文字幕| 在线免费十八禁| 国产视频内射| 亚洲欧美日韩另类电影网站 | 中文字幕人妻熟人妻熟丝袜美| 少妇熟女欧美另类| 欧美成人a在线观看| 欧美另类一区| 国产精品av视频在线免费观看| 国产爱豆传媒在线观看| 黄色日韩在线| 午夜免费男女啪啪视频观看| 18禁动态无遮挡网站| 少妇 在线观看| 亚洲四区av| 免费观看性生交大片5| 日本猛色少妇xxxxx猛交久久| 一二三四中文在线观看免费高清| av女优亚洲男人天堂| 亚洲精品国产av蜜桃| 波多野结衣巨乳人妻| 欧美日韩综合久久久久久| 2021天堂中文幕一二区在线观| 九九久久精品国产亚洲av麻豆| 激情五月婷婷亚洲| 99久国产av精品国产电影| 中文字幕免费在线视频6| 搞女人的毛片| 国产乱人视频| 自拍欧美九色日韩亚洲蝌蚪91 | 全区人妻精品视频| 99久久中文字幕三级久久日本| 大香蕉97超碰在线| 黄色欧美视频在线观看| 亚洲国产高清在线一区二区三| 99久久中文字幕三级久久日本| 一二三四中文在线观看免费高清| 久久99热6这里只有精品| 精品视频人人做人人爽| 国产有黄有色有爽视频| 免费播放大片免费观看视频在线观看| 日韩强制内射视频| 男人和女人高潮做爰伦理| 欧美成人午夜免费资源| 在线播放无遮挡| 超碰97精品在线观看| 久久精品熟女亚洲av麻豆精品| 99视频精品全部免费 在线| 蜜臀久久99精品久久宅男| 日韩大片免费观看网站| 美女内射精品一级片tv| 国产精品蜜桃在线观看| 伊人久久国产一区二区| 成年版毛片免费区| 中文欧美无线码| xxx大片免费视频| 欧美激情在线99| 国产v大片淫在线免费观看| 国产又色又爽无遮挡免| 交换朋友夫妻互换小说| 日本wwww免费看| 王馨瑶露胸无遮挡在线观看| 日韩电影二区| 男的添女的下面高潮视频| 日韩av免费高清视频| 欧美国产精品一级二级三级 | 久久久欧美国产精品| 亚洲真实伦在线观看| 观看美女的网站| 美女视频免费永久观看网站| 成人亚洲精品一区在线观看 | 人人妻人人澡人人爽人人夜夜| 色吧在线观看| 最近2019中文字幕mv第一页| 18+在线观看网站| 蜜臀久久99精品久久宅男| 国产精品女同一区二区软件| 欧美 日韩 精品 国产| 国产白丝娇喘喷水9色精品| 成人国产麻豆网| 欧美xxxx黑人xx丫x性爽| 久久久久久久久久久丰满| 国产有黄有色有爽视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产乱人偷精品视频| 天堂中文最新版在线下载 | 男插女下体视频免费在线播放| 国产免费又黄又爽又色| 精品国产一区二区三区久久久樱花 | 免费大片18禁| 日韩 亚洲 欧美在线| 男女无遮挡免费网站观看| 91午夜精品亚洲一区二区三区| 日韩强制内射视频| 韩国高清视频一区二区三区| 国产精品国产av在线观看| 久久99蜜桃精品久久| 丰满少妇做爰视频| 午夜精品一区二区三区免费看| 久久精品人妻少妇| 嫩草影院新地址| 老女人水多毛片| 交换朋友夫妻互换小说| 亚洲精品国产色婷婷电影| 一级av片app| av天堂中文字幕网| 欧美一区二区亚洲| 成人二区视频| 国产免费一级a男人的天堂| 在线观看国产h片| 国产91av在线免费观看| 人妻夜夜爽99麻豆av| 69人妻影院| 国产色爽女视频免费观看| 国产精品无大码| 国产亚洲av嫩草精品影院| 国产欧美亚洲国产| 精品一区在线观看国产| 精品少妇久久久久久888优播| 一级毛片我不卡| 亚洲欧美日韩东京热| 七月丁香在线播放| 插阴视频在线观看视频| 国产精品一区二区三区四区免费观看| 国产精品爽爽va在线观看网站| 欧美潮喷喷水| 少妇被粗大猛烈的视频| videossex国产| 欧美+日韩+精品| 五月玫瑰六月丁香| av女优亚洲男人天堂| 亚洲在久久综合| 国产亚洲午夜精品一区二区久久 | 国产在线一区二区三区精| av.在线天堂| 国产精品女同一区二区软件| av专区在线播放| 九九久久精品国产亚洲av麻豆| 女人被狂操c到高潮| 精品人妻偷拍中文字幕| 欧美日韩国产mv在线观看视频 | 老师上课跳d突然被开到最大视频| 中文天堂在线官网| 国产高潮美女av| 色视频在线一区二区三区| av黄色大香蕉| 只有这里有精品99| 18禁裸乳无遮挡免费网站照片| 欧美97在线视频| 精品人妻偷拍中文字幕| 欧美少妇被猛烈插入视频| 蜜桃亚洲精品一区二区三区| 欧美区成人在线视频| 校园人妻丝袜中文字幕| 国产精品麻豆人妻色哟哟久久| 亚洲欧美精品专区久久| 亚洲精品影视一区二区三区av| 国产日韩欧美在线精品| 只有这里有精品99| 性色avwww在线观看| 国产伦精品一区二区三区四那| 亚洲怡红院男人天堂| 欧美xxxx性猛交bbbb| 久久久久久久久久久免费av| 夫妻午夜视频| 久久精品国产亚洲av天美| 欧美xxⅹ黑人| 青青草视频在线视频观看| 中文在线观看免费www的网站| 午夜日本视频在线| 欧美一区二区亚洲| 夜夜爽夜夜爽视频| 99久久精品热视频| 成年免费大片在线观看| 熟女电影av网| av在线播放精品| 日韩成人av中文字幕在线观看| 欧美97在线视频| 日本午夜av视频| 麻豆成人午夜福利视频| h日本视频在线播放| 啦啦啦在线观看免费高清www| 国产一区二区三区综合在线观看 | 欧美三级亚洲精品| 丝袜美腿在线中文| 一个人看视频在线观看www免费| 久久久精品94久久精品| 人人妻人人看人人澡| 欧美少妇被猛烈插入视频| 成人综合一区亚洲| 亚洲av男天堂| 丝瓜视频免费看黄片| 免费观看a级毛片全部| a级毛色黄片| 亚洲av国产av综合av卡| 欧美xxxx性猛交bbbb| 欧美xxⅹ黑人| 国产淫片久久久久久久久| 国产免费一级a男人的天堂| 色婷婷久久久亚洲欧美| 国产成人精品久久久久久| 一区二区三区免费毛片| 国产成人freesex在线| 亚洲国产精品国产精品| 午夜福利高清视频| 日韩亚洲欧美综合| 国产精品一区二区性色av| 在线观看免费高清a一片| 九草在线视频观看| 在线观看免费高清a一片| 免费av观看视频| 成人漫画全彩无遮挡| av在线亚洲专区| 午夜免费鲁丝| 久久久久久久久久成人| 成人漫画全彩无遮挡| 国产高清三级在线| 午夜日本视频在线| 乱系列少妇在线播放| 亚洲av日韩在线播放| 最近手机中文字幕大全| 黄色配什么色好看| 亚洲av免费在线观看| 尤物成人国产欧美一区二区三区| 美女主播在线视频| 少妇 在线观看| 夫妻性生交免费视频一级片| 一个人观看的视频www高清免费观看| 欧美人与善性xxx| 偷拍熟女少妇极品色| 成人午夜精彩视频在线观看| 高清日韩中文字幕在线| 国产在线一区二区三区精| 街头女战士在线观看网站| 男插女下体视频免费在线播放| 看免费成人av毛片| 男女国产视频网站| 国产亚洲av片在线观看秒播厂| 国产真实伦视频高清在线观看| 好男人在线观看高清免费视频| 国产毛片a区久久久久| 91久久精品国产一区二区三区| av在线app专区| 亚洲精品日韩av片在线观看| 丝袜脚勾引网站| 黄色日韩在线| 18禁在线播放成人免费| 亚洲在线观看片| 精品国产三级普通话版| 青春草国产在线视频| 日本-黄色视频高清免费观看| 亚洲国产精品成人综合色| 人妻 亚洲 视频| 国产亚洲午夜精品一区二区久久 | 制服丝袜香蕉在线| 禁无遮挡网站| 全区人妻精品视频| 一区二区三区乱码不卡18| 免费av不卡在线播放| 男男h啪啪无遮挡| 色综合色国产| 性插视频无遮挡在线免费观看| 精品一区二区三卡| 国产精品一及| 高清在线视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕制服av| 日韩伦理黄色片| 狠狠精品人妻久久久久久综合| 日韩中字成人| 18禁动态无遮挡网站| 久久精品久久久久久噜噜老黄| 久热久热在线精品观看| 2021天堂中文幕一二区在线观| 国产乱来视频区| 午夜日本视频在线| 最近2019中文字幕mv第一页| 欧美精品一区二区大全| 美女cb高潮喷水在线观看| av天堂中文字幕网| 能在线免费看毛片的网站| 国产精品.久久久| 午夜亚洲福利在线播放| 在线精品无人区一区二区三 | 亚洲av免费在线观看| 精品人妻一区二区三区麻豆| 久久久久精品久久久久真实原创| 国产精品蜜桃在线观看| 久久国产乱子免费精品| 成人二区视频| 国产一区二区三区av在线| 少妇熟女欧美另类| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 国产精品.久久久| 亚洲色图综合在线观看| 久久精品国产亚洲av涩爱| 嫩草影院新地址| 国产一区二区三区综合在线观看 | 国产一区二区三区综合在线观看 | 18+在线观看网站| 国产成人精品婷婷| 最近最新中文字幕大全电影3| 日本wwww免费看| 九九爱精品视频在线观看| 黄片wwwwww| 亚洲av不卡在线观看| 亚洲自拍偷在线| 性色av一级| av在线观看视频网站免费| 久久人人爽人人片av| 神马国产精品三级电影在线观看| 乱系列少妇在线播放| 国产乱人视频| 久久人人爽人人片av| av在线观看视频网站免费| 午夜亚洲福利在线播放| 国产毛片在线视频| 九色成人免费人妻av| 午夜免费观看性视频| 久久久久久久午夜电影| 亚洲人成网站在线观看播放| 国产一区二区亚洲精品在线观看| 日韩视频在线欧美| 久久久久久国产a免费观看| 欧美激情国产日韩精品一区| 91精品伊人久久大香线蕉| 最近中文字幕2019免费版| 男人爽女人下面视频在线观看| 91久久精品电影网| 亚洲精品第二区| 午夜免费男女啪啪视频观看| 成人黄色视频免费在线看| 99久久人妻综合| 国产男人的电影天堂91| 国产成人一区二区在线| videossex国产| 交换朋友夫妻互换小说| 国产一级毛片在线| 99热网站在线观看| 黄色欧美视频在线观看| a级毛片免费高清观看在线播放| 搞女人的毛片| 能在线免费看毛片的网站| 免费观看在线日韩|