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

    含裂隙巖體單軸壓縮裂紋擴展機制離散元分析

    2015-02-17 07:41:50蔣明鏡申志福
    巖土力學(xué) 2015年11期
    關(guān)鍵詞:單軸力學(xué)裂隙

    蔣明鏡,張 寧,申志福,陳 賀

    (1.同濟(jì)大學(xué) 土木工程地下建筑與工程系,上海 200092;2.同濟(jì)大學(xué) 土木工程防災(zāi)國家重點實驗室,上海 200092;3.同濟(jì)大學(xué)巖土及地下工程教育部重點實驗室,上海 200092;4.云南省交通規(guī)劃設(shè)計研究院 陸地交通氣象災(zāi)害防治技術(shù)國家工程試驗室,云南 昆明 650011)

    1 引 言

    自然界中的巖石經(jīng)歷了長期的地質(zhì)構(gòu)造作用,使堅硬的巖石材料含有大量的節(jié)理、裂隙等地質(zhì)缺陷,致使巖石的強度降低且易發(fā)生脆性破裂。大量的巖土工程實踐表明,工程的失穩(wěn)破壞與裂紋缺陷的擴展和貫通有密切的關(guān)系[1]。為理清裂紋的擴展與貫通機制,很多學(xué)者以預(yù)制單條、雙條以及多條裂隙巖石為研究對象,從室內(nèi)試驗、理論計算及數(shù)值模擬等方面開展了大量卓有成效的研究工作[2-10]:將單軸壓縮破壞中的裂紋分為翼裂紋和次生裂紋;提出了雙參數(shù)復(fù)合準(zhǔn)則、拉應(yīng)變斷裂準(zhǔn)則、最小J2準(zhǔn)則和有效剪應(yīng)力準(zhǔn)則等壓剪斷裂理論。然而,理論分析很多情況下并不能合理地考慮巖石的各向異性、非均勻及不連續(xù)等特性;室內(nèi)試驗中成功預(yù)制裂隙具有較大的難度,而采用模型材料又不能完全模擬真實巖石的特性[11]。

    而通過大量的細(xì)觀試驗研究,發(fā)現(xiàn)巖石破壞主要是裂紋沿著巖石晶體顆粒邊界的擴展造成的[12]。而離散元正是一種從微觀機制出發(fā),利用非連續(xù)介質(zhì)方法,在不需要較為復(fù)雜的本構(gòu)關(guān)系的前提下,分析材料宏觀破壞過程的數(shù)值計算方法。離散元模擬巖石破損的主要原理為將巖石抽象為巖石晶體顆粒與粒間膠結(jié)物所組成,建立適合的膠結(jié)接觸模型反映巖石的在不同情況下的破壞過程。近年來,國內(nèi)外已有許多學(xué)者采用離散單元法模擬巖石裂紋的擴展過程,分析了巖石的破裂機制,并通過與室內(nèi)試驗資料對比驗證了該方法模擬巖石的可行性[13-14]。然而,通過蔣明鏡等[13-14]的微觀接觸試驗發(fā)現(xiàn),膠結(jié)物質(zhì)的破壞是法向、切向、轉(zhuǎn)動向荷載共同作用的結(jié)果,而現(xiàn)階段的離散元模型并不能完全考慮這些因素。

    本文將基于室內(nèi)微觀膠結(jié)接觸試驗的巖石微觀膠結(jié)模型嵌入離散元軟件,對預(yù)制單裂隙的Lac du Bonnet 花崗巖石進(jìn)行單軸壓縮試驗,對比分析壓縮過程中試樣中應(yīng)力分布與理論計算的差異,并對裂紋擴展判據(jù)的適用性進(jìn)行了對比研究。

    2 離散元數(shù)值模擬

    2.1 模型簡介

    由于巖石中膠結(jié)的存在,常規(guī)的離散元模型不能直接應(yīng)用于巖石的模擬,基于有厚度水泥膠結(jié)顆粒的試驗結(jié)果[15],通過理論分析推導(dǎo)出無厚度膠結(jié)顆粒的強度包線,采用C++語言與FISH 語言編寫用戶自定義的膠結(jié)接觸模型,然后,通過離散元自帶的模型添加接口即可將該模型嵌入。由于篇幅的限制,推導(dǎo)過程另文做詳細(xì)介紹,此處只給出相應(yīng)的結(jié)果。式(1)~(4)為室內(nèi)試驗總結(jié)得到的強度包線,圖1為無厚度的膠結(jié)顆粒的強度包線。

    式中:Rs為峰值剪切荷載;μb為室內(nèi)試驗得到的參數(shù);Rc、Rt分別為壓縮強度和拉伸強度;Fnb、Fnp分別為作用在膠結(jié)和顆粒上的法向力;gs、fs分別為室內(nèi)試驗得到的參數(shù);μp為顆粒的摩擦系數(shù)。

    式中:Rr為峰值扭轉(zhuǎn)荷載;βb為室內(nèi)試驗得到的參數(shù);為兩個膠結(jié)顆粒(半徑分別為R1,R2)的平均半徑,按計算;gr、fr分別為室內(nèi)試驗得到的參數(shù);βp為顆??罐D(zhuǎn)動系數(shù)。

    圖1 膠結(jié)模型強度包面Fig.1 Strength envelop of bond model

    膠結(jié)破壞后,膠結(jié)和顆粒會同時受力,定義此時為參與狀態(tài),則殘余強度包線:

    式中:Rsr為殘余剪切荷載。

    式中:Rrr為殘余扭轉(zhuǎn)荷載。

    選用蔣明鏡等[16-17]提出的考慮抗轉(zhuǎn)動的膠結(jié)模型,其力學(xué)響應(yīng)如圖2 所示。圖2(a)為膠結(jié)接觸模型的法向力學(xué)響應(yīng)。該力學(xué)響應(yīng)中,法向力隨著法向位移的增加而線性增加,當(dāng)拉力達(dá)到膠結(jié)峰值拉伸荷載Rt時,膠結(jié)被拉壞,拉力突然跌落為0。圖2(b)為膠結(jié)接觸模型的切向力學(xué)響應(yīng)。剪切力Fs隨著切向位移us的增加首先線性增加,當(dāng)達(dá)到峰值剪切荷載Rs時,膠結(jié)發(fā)生破壞,接觸點的切向力突然跌落至恒定值,該恒定值由Mohr-Coulomb 準(zhǔn)則計算得到。圖2(c)為膠結(jié)接觸模型的轉(zhuǎn)動向力學(xué)響應(yīng)。相對轉(zhuǎn)角θ 同扭矩M 呈線性關(guān)系,達(dá)到峰值扭轉(zhuǎn)荷載Rr之后,殘余扭轉(zhuǎn)荷載Rrr由顆粒材料形狀和荷載類型進(jìn)行確定。

    圖2 無厚度膠結(jié)模型的力學(xué)響應(yīng)[16-17]Fig.2 Mechanical responses of zero-thickness bond contact model

    2.2 模型參數(shù)的選取

    參照Lac du Bonnet 花崗巖特性[2],本文選擇了如圖3 所示的級配,用分層欠壓法[18]制成平面孔隙比為0.2 的試樣,然后保持一定的豎向壓力為顆粒接觸賦予接觸模型成樣。通過圖4 所示的單軸壓縮試驗和巴西試驗確定了試樣的微觀參數(shù),如表1所示,試驗得到試樣的宏觀參數(shù)為:單軸壓縮強度為205 MPa,拉伸強度為11.36 MPa,泊松比為0.256。

    圖3 DEM 巖石試樣的級配曲線Fig.3 Particle size distribution of DEM rock sample

    圖4 DEM 試驗結(jié)果Fig.4 Results of DEM simulation

    表1 巖石膠結(jié)模型中的微觀參數(shù)Table 1 DEM micro parameters of rock sample

    2.3 預(yù)制裂隙的制樣

    預(yù)制裂隙代表此處物質(zhì)的缺失,因此,先按照前節(jié)方法制成單軸試樣備用,然后采用刪除部分顆粒的方法進(jìn)行預(yù)制裂隙,此時裂隙為張開狀態(tài)。對于預(yù)制的單裂隙,裂隙中心與試樣中心重合,為保持與室內(nèi)試驗相同的預(yù)制裂隙角度,選取與水平方向的角度α 分別為0°、30°、45°、60°、75°進(jìn)行研究;預(yù)制裂紋形式和尺寸如圖5 所示。在預(yù)制裂隙形成之后,循環(huán)至平衡,隨后對巖石離散元試樣進(jìn)行單軸壓縮試驗。

    圖5 預(yù)制裂隙試樣Fig.5 Pre-cracked samples

    3 試驗結(jié)果分析

    3.1 試驗結(jié)果

    表2為各傾角預(yù)制裂隙試樣的破壞形態(tài)圖與室內(nèi)試驗的對比,詳細(xì)內(nèi)容已另文介紹[19],本文只列出結(jié)果說明模型的合理性。

    表2 DEM 模擬與室內(nèi)試驗[19]試樣破壞形態(tài)對比Table 2 Comparison of failure modes between DEM and experimental results [19]

    從表2 可以看到,DEM 模擬裂隙擴展形態(tài)與室內(nèi)試驗的相似度較高,都可以觀察到明顯的翼裂紋擴展:發(fā)啟于裂紋端部向主應(yīng)力方向發(fā)展。試驗結(jié)果證明了模型模擬巖石的適用性,下文將從試樣的應(yīng)力分布和擴展判據(jù)來論述裂隙的擴展機制。

    3.2 裂隙周圍的應(yīng)力分布

    巖石中的裂隙是由巖石中應(yīng)力引起膠結(jié)破壞導(dǎo)致的,預(yù)制裂隙的存在使巖石中應(yīng)力的分布更加復(fù)雜和不均勻。因而,研究預(yù)制裂隙周圍的應(yīng)力分布對于揭示巖石裂隙的萌生和發(fā)育具有重要意義。

    平面問題中的裂紋主要有兩種類型:Ⅰ型(張開型裂紋)、Ⅱ型(滑開型裂紋)。Westergaard[20]提出了一個復(fù)變應(yīng)力函數(shù),通過假定給出了兩種類型裂紋周圍的應(yīng)力場和位移。如圖6 所示,中心有長為2a 的貫穿裂隙,邊界作用拉應(yīng)力σx和σy致使裂紋張開,即張開型裂紋(Ⅰ型);作用剪應(yīng)力τ 使得裂紋相互剪切,即剪切型裂紋(Ⅱ型)。試樣裂隙類型可視為兩者的疊加。

    圖6 裂紋類型Fig.6 Crack modes

    根據(jù)其推導(dǎo)結(jié)論,在拉剪作用下裂紋周圍的全場應(yīng)力的表達(dá)式為

    其中:K?和KⅡ分別為Ⅰ型和Ⅱ型裂紋的應(yīng)力強度因子;r 和θ為圖6 中的幾何參數(shù)。

    對于如圖7 所示的單軸壓縮,試樣承受著上下兩側(cè)的均布壓力q,試樣的中心位置為一張開型裂紋,該裂紋與加載方向的夾角為β,裂紋的長度為2a,設(shè)置笛卡爾坐標(biāo)系(oxy)和極坐標(biāo)系的原點分別在試樣的中心和裂紋的右端點。據(jù)彈性力學(xué)坐標(biāo)旋轉(zhuǎn)公式可得oxy 坐標(biāo)系下的等效邊界應(yīng)力

    圖7 裂紋擴展判據(jù)的計算模型Fig.7 Model for crack propagation criterion

    將式(6)帶入式(5)即得裂隙巖石單軸壓縮的應(yīng)力分布:

    其中,應(yīng)力強度因子為

    需要說明的是,在前文推導(dǎo)中裂紋受力為拉剪,若裂隙尚未閉合,則式(5)依然適用。在單軸壓縮中裂紋周圍為壓剪受力。若裂紋在壓剪受力時裂隙已經(jīng)閉合,此時會有界面之間的摩擦力。在離散元模擬過程中,觀察到在試樣發(fā)生膠結(jié)破壞,即開裂之前,裂隙并未閉合,因此,式(5)依然適用。

    表3為DEM 數(shù)值試驗和理論分析得到的軸向應(yīng)變0.7εp(εp為峰值應(yīng)力對應(yīng)應(yīng)變)所對應(yīng)的應(yīng)力分布圖。由表可以看出:(1)盡管由于各向同性的假設(shè),預(yù)制裂隙試樣中理論計算應(yīng)力分布的數(shù)值與DEM 模擬得到的數(shù)值有差異,但當(dāng)預(yù)制角度小于75°時,在定性規(guī)律上是一致的,在預(yù)制裂隙的周圍出現(xiàn)了明顯的應(yīng)力集中現(xiàn)象,側(cè)向應(yīng)力σxx中拉應(yīng)力集中于裂隙的上下兩面、壓應(yīng)力集中于裂隙的左右兩端,而切應(yīng)力則在預(yù)制裂隙的上下兩面和左右兩端呈現(xiàn)相反的方向。產(chǎn)生這樣現(xiàn)象的原因應(yīng)該為軸向的壓縮使裂隙有閉合的趨勢,此時對于裂隙兩端較大的壓應(yīng)變使裂隙上下兩面產(chǎn)生拉應(yīng)變,從而使裂隙周圍產(chǎn)生拉應(yīng)力。(2)而當(dāng)預(yù)制裂隙角度為75°時,理論計算值中側(cè)向應(yīng)力σxx處于拉應(yīng)力狀態(tài),而DEM 計算結(jié)果表明,此時試樣處于拉壓應(yīng)力狀態(tài),但與軸向應(yīng)力σyy相比數(shù)值較小,說明此時裂隙對試樣整體的應(yīng)力分布影響已經(jīng)不大。(3)隨著預(yù)制裂隙角度的增加,在達(dá)到0.7εp時,軸向應(yīng)力會增加,而側(cè)向應(yīng)力的最大拉應(yīng)力卻越來越小,說明傾角較大的裂隙可以有效地減小應(yīng)力的集中,也因此預(yù)制裂隙角度較大的試樣具有更高的強度。

    3.3 裂紋擴展判據(jù)

    裂紋的擴展總是從應(yīng)力集中的端部開始的,而從上節(jié)的分析中可知,裂隙的左右兩端是壓應(yīng)力集中的地方,裂隙的上下兩面是拉應(yīng)力集中的地方。而裂隙擴展正是從拉壓應(yīng)力的交界部分開始,而由于裂隙周圍應(yīng)力較難測量,因而,許多學(xué)者開始從理論上分析裂紋的擴展方向,并為此提出了最大周向應(yīng)力準(zhǔn)則、最大拉應(yīng)變準(zhǔn)則、最小應(yīng)變能密度準(zhǔn)則、最大能量釋放準(zhǔn)則、最小J2準(zhǔn)則等。

    表3 DEM 數(shù)值試驗和理論分析得到的單裂隙試樣應(yīng)力場Table 3 Stress fields in rock samples containing single flaw obtained from DEM numerical simulations and theoretical analyses

    基于圖7 中坐標(biāo)系,巖石受壓剪荷載作用時,裂紋尖端的奇異應(yīng)力在極坐標(biāo)中應(yīng)力分量表達(dá)式為

    巖石受壓剪荷載作用時,裂紋尖端的奇異應(yīng)力場在笛卡爾坐標(biāo)系中的應(yīng)力分量表達(dá)式為

    (1)最大周向應(yīng)力準(zhǔn)則[21]

    最大周向拉應(yīng)力準(zhǔn)則認(rèn)為,裂紋沿著最大的周向應(yīng)力σθθmax擴展,所以該方向滿足下列條件:

    聯(lián)立式(8)、(9)、(10)、(12)可得到裂紋擴展角為

    (2)最大拉應(yīng)變準(zhǔn)則[22]

    最大拉應(yīng)變準(zhǔn)則認(rèn)為,環(huán)向拉應(yīng)變 εθθ達(dá)到最大值的方向為裂紋擴展方向。對于平面應(yīng)變問題,由巖石斷裂力學(xué)可知,裂紋尖端的環(huán)向拉應(yīng)變?yōu)?/p>

    式中:E為彈性模量;ν為泊松比。

    將式(8)、(9)代入式(14),推導(dǎo)得到:

    (3)最小應(yīng)變能密度因子準(zhǔn)則[23]

    最小應(yīng)變能密度準(zhǔn)則認(rèn)為,裂紋沿應(yīng)變能力最小方向發(fā)展。由斷裂力學(xué)可知,單位體積的應(yīng)變能(又稱為應(yīng)變能密度)為

    聯(lián)立式(8)、(9)、(11)、(16),求導(dǎo)得到極小值時裂紋的擴展角θ0為

    (4)最大能量釋放率準(zhǔn)則

    最大能量釋放準(zhǔn)則認(rèn)為,裂紋沿著能量釋放率最大的方向擴展,由文獻(xiàn)[24]可得到能量釋放率G的表達(dá)式為

    將式(8)、(9)代入式(18),求得G 取極大值時的裂紋擴展角滿足:

    (5)最小J2準(zhǔn)則

    最小J2準(zhǔn)則認(rèn)為,裂紋沿J2最小的方向擴展,J2的表達(dá)式為

    將式(11)代入式(20),求得J2極小值時得到裂紋的起裂角θ0滿足[25]:

    圖8為DEM 數(shù)值試驗、室內(nèi)試驗和理論分析得到的預(yù)制裂隙傾角與裂紋起裂角之間的關(guān)系。由圖可見,三者的結(jié)果,裂紋的起裂角均隨著預(yù)制裂隙傾角的增加而減小;除了最小J2準(zhǔn)則之外,其他的理論分析結(jié)果與DEM 數(shù)值試驗及室內(nèi)試驗得到的結(jié)果在定性上的一致性良好;DEM 數(shù)值試驗與室內(nèi)試驗得到的裂紋起裂角在數(shù)值上的一致性較好。另外,最大周向應(yīng)力準(zhǔn)則與最大能量釋放率準(zhǔn)則得到的裂紋起裂角與DEM 數(shù)值試驗得到的結(jié)果較為相似。以上分析表明,最大周向應(yīng)力準(zhǔn)則和最大能量釋放率準(zhǔn)則在預(yù)測裂紋起裂角上較為準(zhǔn)確。

    圖8 不同方法所得預(yù)制裂隙傾角與裂紋起裂角之間的關(guān)系Fig.8 Relationships between the pre-existing crack inclination and crack initiation inclination obtained by different methods

    5 結(jié) 論

    (1)由室內(nèi)試驗總結(jié)得到的無厚度由抗轉(zhuǎn)動膠結(jié)模型適用于模擬巖石裂隙的擴展規(guī)律。

    (2)與DEM 模擬結(jié)果相比,理論分析得到的應(yīng)力在定性上一致,而在定量上有差距。主要是由于理論分析中假設(shè)了巖石為各向同性彈性體,這不能反映真實巖石的非均勻、不連續(xù)的特性。當(dāng)預(yù)制裂隙角度較小時,可以觀察到明顯的應(yīng)力集中現(xiàn)象:在側(cè)向應(yīng)力中,裂隙左右兩端壓應(yīng)力集中,并由此導(dǎo)致了裂紋上下兩面的拉應(yīng)力集中;而當(dāng)裂隙角度較大時,裂紋周圍的應(yīng)力集中現(xiàn)象已經(jīng)不明顯,這導(dǎo)致了DEM 結(jié)果與理論分析的差別。

    (3)DEM 數(shù)值試驗得到的裂紋起裂角與室內(nèi)試驗得到的結(jié)果較一致,且與最大周應(yīng)力準(zhǔn)則和最大能量釋放率準(zhǔn)則得到的裂紋擴展角較為吻合。

    [1]陳衛(wèi)忠,李術(shù)才,朱維申,等.巖石裂紋擴展的實驗與數(shù)值分析試驗[J].巖石力學(xué)與工程學(xué)報,2003,22(1):18-23.CHEN Wei-zhong,LI Shu-cai,ZHU Wei-shen,et al.Experimental and numerical research on crack propagation in rock under compression[J].Chinese Journal of Rock Mechanics and Engineering,2003 22(1):18-23.

    [2]MARTIN C D,CHANDLER N A.The progressive fracture of Lac du Bonnet granite[J].International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts,1994,31(6):643-659

    [3]HEEKWANG L,SEOKWON J.An experimental and numerical study of fracture coalescence in pre-cracked specimen under uniaxial compression[J].International Journal of Solids and Structures,2011,48(6):979-999.

    [4]邵冬亮,李術(shù)才,李明田,等.單軸壓縮下不同裂紋制作方式的裂紋破壞分析及其聲發(fā)射特征研究[J].山東大學(xué)學(xué)報(工學(xué)版),2011,41(3):131-136.SHAO Dong-liang,LI Shu-cai,LI Ming-tian,et al.Study of failure analysis and the acoustic emission signature of 3-D crack from different origins under uniaxial compression[J].Journal of Shandong University:Engineering Science,2011,41(3):131-136.

    [5]LI Y P,CHEN L Z,WANG Y H.Experimental research on pre-cracked marble under compression[J].International Journal of Solids and Structures,2005,42(9):2505-2516.

    [6]YANG S D.Experimental study on mechanical behavior of brittle marble samples containing different flaws under uniaxial compression[J].Engineering Fracture Mechanics,2009,76(12):1833-1845.

    [7]REYES O,EINSTEIN H H.Failure mechanisms of fractured rock—A fracture coalescence model[C]//Proceedings of 7th International Congress on Rock Mechanics.Germany:Aachen,1991,333-340.

    [8]TANG C A.Numerical simulation of progressive rock failure and associated seismicity[C]//International Journal of Rock Mechanics and Mining Sciences.[S.l.]:[s.n.],1997,34(2):249-262.

    [9]TANG C A,WONG R H C,CHAU K T,et al.Analysis of crack coalescence in rock-like materials containing three flaws—part II:Numerical approach[J].International Journal of Rock Mechanics &Mining Sciences,2001,38(7):925-939.

    [10]IRWIN G R.Analysis of stress and strains near the end of a crack extension force[J].Journal of Applied Mechanics,1957,24:361-364.

    [11]李銀平,王元漢,陳龍珠,等.含預(yù)制裂紋的大理巖的壓剪試驗分析[J].巖土工程學(xué)報,2004,26(1):120-124.LI Yin-ping,WANG Yuan-han,CHEN Long-zhu,et al.Experimental research on pre-existing cracks in marble under compression[J].Chinese Journal of Geotechnical Engineering,2004,26(1):120-124.

    [12]陳芳,秦昊.細(xì)觀尺度下巖石沿晶斷裂應(yīng)力強度因子計算研究[J].巖土力學(xué),2011,32(3):941-945.CHEN Fang,QIN Hao.Study of intercrystalline factures stress intensity factor of rock in meso-scale[J].Rock and Soil Mechanics,2011,32(3):941-945.

    [13]蔣明鏡,孫渝剛,李立青.水泥膠結(jié)顆粒的微觀力學(xué)模型試驗[J].巖土力學(xué),2012,33(5):1293-1299.JIANG Ming-jing,SUN Yu-gang,LI Li-qing.Experimental investigation on micromechanical model of cementbonded particles[J].Rock and Soil Mechanics,2012,33(5):1293-1299.

    [14]蔣明鏡,周雅萍,陳賀.不同膠結(jié)厚度下粒間膠結(jié)力學(xué)特性的試驗研究[J].巖土力學(xué),2013,34(5):1264-1273.JIANG Ming-jing,ZHOU Ya-ping,CHEN He.Experimental study of mechanical behaviors of bonded granules under different bond thickness[J].Rock and Soil Mechanics,2013,34(5):1264-1273.

    [15]蔣明鏡,陳賀,劉芳.巖石微觀膠結(jié)模型及離散元數(shù)值仿真方法初探[J].巖石力學(xué)與工程學(xué)報,2013,32(1):15-23.JIANG Ming-jing,CHEN He,LIU Fang.A microscopic bond model for rock and preliminary study of numerical simulation method by distinct element method[J].Chinese Journal of Rock Mechanics and Engineering,2013,32(1):15-23.

    [16]JIANG M J,YU H S,HARRIS D.A novel discrete model for granular material incorporating rolling resistance[J].Computers and Geotechnics,2005,32(5):340-357.

    [17]JIANG M J,YU H S,HARRIS D.Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses[J].International Journal for Numerical and Analytical Methods in Geomechnics,2006,30(8):723-761.

    [18]JIANG M J,KONRAD J M,LEROUEIL S.An efficient technique for generating homogeneous specimens for DEM studies[J].Computers and Geotechnics,2003,30(7):579-597.

    [19]蔣明鏡,陳賀.預(yù)制單裂隙巖石的裂紋擴展機制的離散元數(shù)值分析[J].巖石力學(xué)與工程學(xué)報,在審.JIANG Ming-jing,CHEN He.Distinct element analysis of crack propagation in rocks containing pre-existing single flaws[J].Chinese Journal of Rock Mechanics and Engineering,under review.

    [20]WESTERGAARD H M.Bearing pressures and cracks[J].Journal of Applied Mechanics,1939,61:49-53.

    [21]PARIS P C,ERDOGAN F.A critical analysis of crackpropagation laws[J].Journal of Fluids Engineering,1963,85(4):528-533.

    [22]MAZARS J,PIJAUDIER C G.Continuum damage theory—Application to concrete[J].Journal of Engineering Mechanics(ASCE),1987,115(2):345-365.

    [23]SIH G C,BARTHELEMY B M.Mixed mode fatigue crack growth predictions[J].Engineering Fracture Mechanics,1980,13:439-451.

    [24]崔德渝,吳紹富,張行.復(fù)合裂紋的能量釋放率計算與斷裂判據(jù)[J].北京航空航天大學(xué)學(xué)報,1990,3:43-49.CUI De-yu,WU Shao-fu,ZHANG Xing.Energy release rate and failure criterion for mixed mode cracking[J].Journal of Beijing University of Aeronautics and Astronautics,1990,3:43-49.

    [25]周家文,徐衛(wèi)亞,石崇.基于破壞準(zhǔn)則的巖石壓剪斷裂判據(jù)研究[J].巖石力學(xué)與工程學(xué)報,2007,26(6):1194-1201.ZHOU Jia-wen,XU Wei-ya,SHI Chong.Investigation on compression-shear fracture criterion of rock based on failure criteria[J].Chinese Journal of Rock Mechanics and Engineering,2007,26(6):1194-1201.

    猜你喜歡
    單軸力學(xué)裂隙
    力學(xué)
    裂隙腦室綜合征的診斷治療新進(jìn)展
    弟子規(guī)·余力學(xué)文(十)
    快樂語文(2021年11期)2021-07-20 07:41:32
    弟子規(guī)·余力學(xué)文(四)
    快樂語文(2020年30期)2021-01-14 01:05:28
    單軸壓縮條件下巖石峰后第Ⅱ種類型應(yīng)力——應(yīng)變曲線的新解釋
    CFRP-鋼復(fù)合板的單軸拉伸力學(xué)性能
    裂隙燈檢查的個性化應(yīng)用(下)
    單軸應(yīng)變Si NMOS電流模型研究
    電子測試(2017年12期)2017-12-18 06:35:42
    斜單軸跟蹤式光伏組件的安裝傾角優(yōu)化設(shè)計
    力學(xué) 等
    亚洲欧洲国产日韩| 久久人人爽人人片av| 在线观看人妻少妇| 99re6热这里在线精品视频| 亚洲欧美一区二区三区黑人 | 伊人久久国产一区二区| 亚洲,欧美精品.| 国产男女超爽视频在线观看| 欧美bdsm另类| 黄色怎么调成土黄色| 国产老妇伦熟女老妇高清| 中文字幕亚洲精品专区| 国产一区二区三区综合在线观看| 国产黄频视频在线观看| 成人影院久久| 最近2019中文字幕mv第一页| 成年女人在线观看亚洲视频| 男女免费视频国产| 成年女人毛片免费观看观看9 | av片东京热男人的天堂| 国产免费视频播放在线视频| 老汉色∧v一级毛片| 亚洲欧美成人精品一区二区| 国产黄色免费在线视频| 97在线视频观看| 99久国产av精品国产电影| 日本91视频免费播放| 欧美人与善性xxx| 国产黄色免费在线视频| 国产高清国产精品国产三级| 久久久久国产精品人妻一区二区| 久久这里有精品视频免费| 亚洲国产精品一区二区三区在线| 精品一区二区三区四区五区乱码 | 国产熟女午夜一区二区三区| 国产精品一区二区在线不卡| 亚洲三区欧美一区| a级片在线免费高清观看视频| 久久久亚洲精品成人影院| 亚洲综合色网址| 中文字幕精品免费在线观看视频| 新久久久久国产一级毛片| 国产白丝娇喘喷水9色精品| 国产精品一国产av| 免费av中文字幕在线| 国产精品99久久99久久久不卡 | 成人毛片a级毛片在线播放| 久久久久视频综合| 中文字幕精品免费在线观看视频| 在线观看一区二区三区激情| 亚洲欧美精品自产自拍| 国产精品人妻久久久影院| 婷婷色综合www| 日日啪夜夜爽| 人体艺术视频欧美日本| 亚洲av.av天堂| av国产久精品久网站免费入址| 亚洲图色成人| 久久久久久久久久人人人人人人| 成人国产麻豆网| 黄色毛片三级朝国网站| 午夜精品国产一区二区电影| av有码第一页| 日韩制服骚丝袜av| 婷婷色综合www| 午夜影院在线不卡| 成人毛片60女人毛片免费| 2018国产大陆天天弄谢| 777久久人妻少妇嫩草av网站| 在线观看免费视频网站a站| 天天操日日干夜夜撸| 亚洲av福利一区| 热99国产精品久久久久久7| 久久久久久久亚洲中文字幕| 亚洲国产精品一区二区三区在线| 多毛熟女@视频| 久久精品夜色国产| 国产精品不卡视频一区二区| 国产精品一国产av| 精品少妇内射三级| 国产日韩欧美亚洲二区| 亚洲第一区二区三区不卡| 日韩中文字幕视频在线看片| 中文乱码字字幕精品一区二区三区| 精品视频人人做人人爽| 91午夜精品亚洲一区二区三区| 精品一区在线观看国产| 建设人人有责人人尽责人人享有的| 午夜福利,免费看| 国产精品成人在线| 日韩成人av中文字幕在线观看| 国产一级毛片在线| 久久女婷五月综合色啪小说| 女人精品久久久久毛片| 99久国产av精品国产电影| 日韩电影二区| 欧美日韩成人在线一区二区| 午夜免费鲁丝| 日韩制服丝袜自拍偷拍| 一级毛片黄色毛片免费观看视频| 亚洲精品乱久久久久久| 亚洲欧美中文字幕日韩二区| 蜜桃在线观看..| 欧美日韩视频精品一区| 在线观看免费高清a一片| 乱人伦中国视频| 美女脱内裤让男人舔精品视频| 亚洲成国产人片在线观看| 日本欧美国产在线视频| 另类亚洲欧美激情| 亚洲色图 男人天堂 中文字幕| 久久综合国产亚洲精品| 999久久久国产精品视频| 色94色欧美一区二区| 成人国语在线视频| 下体分泌物呈黄色| 国产伦理片在线播放av一区| 日韩伦理黄色片| 好男人视频免费观看在线| 亚洲欧美色中文字幕在线| 国产女主播在线喷水免费视频网站| 国产精品欧美亚洲77777| 麻豆乱淫一区二区| 亚洲国产日韩一区二区| 亚洲av男天堂| 午夜激情av网站| 人人妻人人爽人人添夜夜欢视频| 日本午夜av视频| 热re99久久国产66热| 亚洲国产毛片av蜜桃av| 午夜日本视频在线| av在线播放精品| 久久人人97超碰香蕉20202| 国产av码专区亚洲av| 各种免费的搞黄视频| 亚洲国产色片| 中文字幕人妻丝袜一区二区 | 亚洲国产欧美网| 巨乳人妻的诱惑在线观看| 1024视频免费在线观看| 精品少妇久久久久久888优播| 男人舔女人的私密视频| 久久人妻熟女aⅴ| 精品一区在线观看国产| 亚洲精品国产av蜜桃| 99国产综合亚洲精品| 国产亚洲欧美精品永久| 国产亚洲一区二区精品| 午夜福利在线免费观看网站| xxx大片免费视频| 国产精品免费视频内射| 美女高潮到喷水免费观看| 亚洲成国产人片在线观看| 一级a爱视频在线免费观看| 伦精品一区二区三区| 一区在线观看完整版| 日韩欧美精品免费久久| 精品人妻在线不人妻| 免费人妻精品一区二区三区视频| 亚洲国产欧美网| 一区福利在线观看| 国产 一区精品| 最近最新中文字幕大全免费视频 | 免费在线观看黄色视频的| 婷婷色综合www| kizo精华| 亚洲伊人久久精品综合| 水蜜桃什么品种好| 岛国毛片在线播放| 熟女少妇亚洲综合色aaa.| 另类精品久久| 你懂的网址亚洲精品在线观看| 久久久精品94久久精品| 中文字幕制服av| 丝袜美腿诱惑在线| 成人影院久久| 新久久久久国产一级毛片| 天天躁夜夜躁狠狠久久av| 我的亚洲天堂| 日本欧美视频一区| 国产 一区精品| 一本大道久久a久久精品| 久久精品aⅴ一区二区三区四区 | 麻豆精品久久久久久蜜桃| 大香蕉久久网| √禁漫天堂资源中文www| 中文字幕制服av| 大话2 男鬼变身卡| 日本爱情动作片www.在线观看| 国产野战对白在线观看| 国产精品国产三级专区第一集| 国产成人a∨麻豆精品| 日韩精品免费视频一区二区三区| 婷婷色综合大香蕉| 美女国产高潮福利片在线看| 九草在线视频观看| 国产精品免费大片| 晚上一个人看的免费电影| 肉色欧美久久久久久久蜜桃| 久久久久久久精品精品| tube8黄色片| 欧美日韩成人在线一区二区| 精品久久久精品久久久| 黄片小视频在线播放| 亚洲av.av天堂| 欧美成人午夜精品| 精品国产露脸久久av麻豆| 久久久久人妻精品一区果冻| 精品亚洲成a人片在线观看| 国产极品天堂在线| 精品国产国语对白av| 考比视频在线观看| 亚洲av男天堂| 999精品在线视频| 欧美日韩一区二区视频在线观看视频在线| 在线观看美女被高潮喷水网站| 在线观看人妻少妇| 欧美人与性动交α欧美精品济南到 | 下体分泌物呈黄色| 成人亚洲精品一区在线观看| 精品一区在线观看国产| 丝袜喷水一区| 日本色播在线视频| 日韩 亚洲 欧美在线| 国产白丝娇喘喷水9色精品| 丰满乱子伦码专区| 亚洲av在线观看美女高潮| 日韩一本色道免费dvd| 卡戴珊不雅视频在线播放| 国产日韩欧美在线精品| 99香蕉大伊视频| 亚洲欧美中文字幕日韩二区| 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 日本91视频免费播放| 国产 一区精品| 免费日韩欧美在线观看| 伦理电影大哥的女人| 男女边吃奶边做爰视频| 一个人免费看片子| 久久精品久久久久久噜噜老黄| 亚洲欧美精品综合一区二区三区 | 在线观看美女被高潮喷水网站| av在线app专区| 日本av免费视频播放| 曰老女人黄片| 最近中文字幕2019免费版| 亚洲av电影在线观看一区二区三区| 欧美av亚洲av综合av国产av | 日日爽夜夜爽网站| 久久久欧美国产精品| 成人国语在线视频| 久久人人97超碰香蕉20202| freevideosex欧美| 男人爽女人下面视频在线观看| 国产精品.久久久| 欧美精品av麻豆av| 天天影视国产精品| 国产麻豆69| 亚洲精品在线美女| 制服丝袜香蕉在线| 国产精品久久久久久av不卡| 欧美黄色片欧美黄色片| 午夜激情久久久久久久| 极品少妇高潮喷水抽搐| 国产 一区精品| 国产 精品1| 最近中文字幕2019免费版| 色视频在线一区二区三区| 精品一品国产午夜福利视频| 国产97色在线日韩免费| 男人舔女人的私密视频| 街头女战士在线观看网站| 免费日韩欧美在线观看| 国产av码专区亚洲av| 国产av国产精品国产| 午夜福利一区二区在线看| 99热网站在线观看| 在线亚洲精品国产二区图片欧美| 大香蕉久久成人网| 侵犯人妻中文字幕一二三四区| 亚洲成色77777| 国产一区二区 视频在线| 国产淫语在线视频| 精品久久蜜臀av无| 精品亚洲乱码少妇综合久久| 又大又黄又爽视频免费| 热re99久久精品国产66热6| 成人毛片60女人毛片免费| 丰满少妇做爰视频| 国产一区二区激情短视频 | av有码第一页| 十分钟在线观看高清视频www| 国产 精品1| 高清视频免费观看一区二区| 国产极品天堂在线| 精品少妇久久久久久888优播| av线在线观看网站| 亚洲av日韩在线播放| 亚洲欧洲精品一区二区精品久久久 | 国产麻豆69| www日本在线高清视频| 国产精品国产av在线观看| 美女主播在线视频| 精品久久久精品久久久| av视频免费观看在线观看| 亚洲天堂av无毛| 一本久久精品| 18+在线观看网站| 看非洲黑人一级黄片| 亚洲av免费高清在线观看| 1024香蕉在线观看| 日本色播在线视频| 成人影院久久| 亚洲精品在线美女| 日韩制服丝袜自拍偷拍| 在线天堂最新版资源| 国产精品免费大片| 成人二区视频| 黄片小视频在线播放| 国产极品天堂在线| 国产1区2区3区精品| 天天躁狠狠躁夜夜躁狠狠躁| 麻豆av在线久日| 国精品久久久久久国模美| 在线观看一区二区三区激情| 只有这里有精品99| 好男人视频免费观看在线| 久久精品aⅴ一区二区三区四区 | 最近中文字幕高清免费大全6| 日本91视频免费播放| 亚洲精品视频女| 国产成人a∨麻豆精品| 精品人妻在线不人妻| 久久精品久久久久久噜噜老黄| 色婷婷久久久亚洲欧美| 亚洲精品美女久久av网站| 日韩中文字幕视频在线看片| 超碰97精品在线观看| 国产成人精品久久二区二区91 | 最黄视频免费看| 最近中文字幕高清免费大全6| 午夜福利视频在线观看免费| 久久免费观看电影| 午夜日韩欧美国产| 色网站视频免费| 国产av国产精品国产| 不卡视频在线观看欧美| 91在线精品国自产拍蜜月| 十八禁网站网址无遮挡| 日韩成人av中文字幕在线观看| 午夜福利,免费看| 亚洲av电影在线观看一区二区三区| 久久毛片免费看一区二区三区| 欧美少妇被猛烈插入视频| 日韩电影二区| 欧美少妇被猛烈插入视频| 午夜免费观看性视频| 国产成人精品婷婷| 久久这里有精品视频免费| 午夜久久久在线观看| 成年女人在线观看亚洲视频| 97在线视频观看| 人妻 亚洲 视频| 熟女少妇亚洲综合色aaa.| 国产精品香港三级国产av潘金莲 | 七月丁香在线播放| 一级爰片在线观看| www.精华液| 久热这里只有精品99| 日韩av免费高清视频| 欧美97在线视频| 久久久国产一区二区| 男人添女人高潮全过程视频| 一级毛片电影观看| 国产毛片在线视频| 日本爱情动作片www.在线观看| 亚洲精品美女久久av网站| 久久97久久精品| 美女高潮到喷水免费观看| 久久韩国三级中文字幕| 国产精品av久久久久免费| 好男人视频免费观看在线| 99精国产麻豆久久婷婷| 免费在线观看视频国产中文字幕亚洲 | 久久热在线av| 老汉色av国产亚洲站长工具| av卡一久久| 色94色欧美一区二区| 精品国产国语对白av| 国产av精品麻豆| 亚洲综合精品二区| 欧美精品一区二区大全| 国产乱来视频区| 亚洲美女黄色视频免费看| 另类精品久久| 日韩不卡一区二区三区视频在线| 久久久国产欧美日韩av| 日本猛色少妇xxxxx猛交久久| 80岁老熟妇乱子伦牲交| 两个人看的免费小视频| 精品酒店卫生间| 亚洲婷婷狠狠爱综合网| 国产成人精品久久久久久| 国产精品嫩草影院av在线观看| 高清黄色对白视频在线免费看| 中文乱码字字幕精品一区二区三区| 精品福利永久在线观看| 人妻少妇偷人精品九色| 成人手机av| 18禁动态无遮挡网站| 视频在线观看一区二区三区| 亚洲av免费高清在线观看| 国产成人精品久久久久久| 老司机亚洲免费影院| av有码第一页| 久久99热这里只频精品6学生| 伦理电影大哥的女人| 亚洲国产最新在线播放| 男女免费视频国产| 久久人人爽人人片av| 国产精品一国产av| 亚洲国产精品国产精品| 精品午夜福利在线看| 精品人妻在线不人妻| 在线观看国产h片| 深夜精品福利| 91精品伊人久久大香线蕉| 精品国产一区二区三区久久久樱花| 不卡av一区二区三区| 老鸭窝网址在线观看| 国产精品一二三区在线看| 免费观看a级毛片全部| www.自偷自拍.com| xxx大片免费视频| 赤兔流量卡办理| 我的亚洲天堂| 欧美日韩综合久久久久久| 国产精品嫩草影院av在线观看| 一级,二级,三级黄色视频| 热re99久久精品国产66热6| 99热国产这里只有精品6| 国产野战对白在线观看| 男女边吃奶边做爰视频| 国产在线视频一区二区| 午夜福利影视在线免费观看| 欧美精品人与动牲交sv欧美| 999精品在线视频| 国产麻豆69| 日本欧美视频一区| 国产日韩一区二区三区精品不卡| 看十八女毛片水多多多| 久久久久网色| 中文欧美无线码| 午夜精品国产一区二区电影| 久久久久视频综合| 国产精品不卡视频一区二区| 美女福利国产在线| 国产成人精品无人区| 在线免费观看不下载黄p国产| 青草久久国产| 伊人久久大香线蕉亚洲五| 午夜免费男女啪啪视频观看| 日韩成人av中文字幕在线观看| 亚洲人成77777在线视频| 超碰成人久久| av在线观看视频网站免费| 天天躁夜夜躁狠狠躁躁| 日韩三级伦理在线观看| 亚洲,欧美精品.| 毛片一级片免费看久久久久| 高清黄色对白视频在线免费看| 日韩一区二区三区影片| 国产精品秋霞免费鲁丝片| 欧美在线黄色| 人人妻人人添人人爽欧美一区卜| 大陆偷拍与自拍| 午夜激情av网站| 一区在线观看完整版| videos熟女内射| 亚洲内射少妇av| 亚洲成av片中文字幕在线观看 | 亚洲在久久综合| av卡一久久| 欧美精品一区二区大全| 如何舔出高潮| 在线观看三级黄色| 国产成人一区二区在线| 天天躁日日躁夜夜躁夜夜| 交换朋友夫妻互换小说| 欧美 日韩 精品 国产| 国产成人免费观看mmmm| kizo精华| av国产久精品久网站免费入址| 亚洲av欧美aⅴ国产| 国产无遮挡羞羞视频在线观看| 久久久久久伊人网av| 欧美变态另类bdsm刘玥| 国产免费视频播放在线视频| 色婷婷久久久亚洲欧美| 午夜av观看不卡| 丝袜喷水一区| 久久久久久久久免费视频了| 国产熟女欧美一区二区| 国产精品一二三区在线看| 性少妇av在线| av视频免费观看在线观看| 啦啦啦视频在线资源免费观看| 亚洲国产精品一区二区三区在线| 亚洲精品中文字幕在线视频| 免费不卡的大黄色大毛片视频在线观看| 蜜桃在线观看..| 久久午夜福利片| 91午夜精品亚洲一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产精品三级大全| 亚洲中文av在线| 亚洲欧美中文字幕日韩二区| 免费观看a级毛片全部| av在线app专区| 亚洲欧洲精品一区二区精品久久久 | 18+在线观看网站| 色播在线永久视频| 午夜日本视频在线| 国产乱人偷精品视频| 欧美国产精品一级二级三级| 国产一区二区激情短视频 | 国产熟女午夜一区二区三区| 亚洲第一av免费看| 日本av手机在线免费观看| 欧美激情高清一区二区三区 | 老汉色av国产亚洲站长工具| 老司机影院成人| 亚洲色图 男人天堂 中文字幕| 一级爰片在线观看| 久久99热这里只频精品6学生| 久久精品国产亚洲av高清一级| 成人二区视频| 久久久国产一区二区| 岛国毛片在线播放| 亚洲欧美一区二区三区国产| 国产一区二区在线观看av| 少妇的逼水好多| 9191精品国产免费久久| 只有这里有精品99| 国产一区有黄有色的免费视频| 欧美日韩视频高清一区二区三区二| 免费女性裸体啪啪无遮挡网站| 超碰97精品在线观看| 国产男人的电影天堂91| 国产一区二区激情短视频 | 亚洲精品久久久久久婷婷小说| 丰满饥渴人妻一区二区三| 哪个播放器可以免费观看大片| 狠狠精品人妻久久久久久综合| 人人妻人人澡人人看| 丝袜喷水一区| 最近中文字幕高清免费大全6| 纵有疾风起免费观看全集完整版| 熟妇人妻不卡中文字幕| 久久ye,这里只有精品| 国产有黄有色有爽视频| 在线观看www视频免费| xxx大片免费视频| www.av在线官网国产| 人妻少妇偷人精品九色| 精品一区二区免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 1024香蕉在线观看| 春色校园在线视频观看| 哪个播放器可以免费观看大片| 国产精品香港三级国产av潘金莲 | xxxhd国产人妻xxx| 国产福利在线免费观看视频| 免费黄网站久久成人精品| 丰满乱子伦码专区| 国产有黄有色有爽视频| 黄色 视频免费看| 男女边吃奶边做爰视频| 99精国产麻豆久久婷婷| 婷婷成人精品国产| 午夜久久久在线观看| 老司机影院毛片| 免费看av在线观看网站| 国产激情久久老熟女| 成人亚洲欧美一区二区av| 超碰97精品在线观看| 国产精品免费大片| 色网站视频免费| 婷婷成人精品国产| √禁漫天堂资源中文www| 欧美日韩综合久久久久久| 亚洲精品中文字幕在线视频| 色哟哟·www| 叶爱在线成人免费视频播放| 国产精品久久久久久久久免| 97在线视频观看| 精品国产乱码久久久久久男人| 香蕉精品网在线| 男女边吃奶边做爰视频| 欧美激情高清一区二区三区 | 人人妻人人澡人人爽人人夜夜| 久久人人爽人人片av| av又黄又爽大尺度在线免费看| 女性生殖器流出的白浆| 天堂俺去俺来也www色官网| 2018国产大陆天天弄谢| 欧美日韩综合久久久久久| 国产精品久久久久久久久免| 男女国产视频网站| 欧美xxⅹ黑人| 老熟女久久久| xxxhd国产人妻xxx| 午夜日本视频在线| 日韩视频在线欧美|