收稿日期：2022-09-06""" 修回日期：2023-04-09

基金項目：國家自然科學(xué)基金資助項目（No.52074143；51874176；51874166）；煤炭科學(xué)研究院科技發(fā)展基金資助項目（No.2022CX-Ⅱ-04）；遼寧工程技術(shù)大學(xué)學(xué)科創(chuàng)新團隊資助項目（No.LNTU20TD-11）

通信作者：李宏艷。E-mail：lhylhb@163.com

引用格式：

徐海鐸，李宏艷，梁冰，等. 壓剪作用下單一閉合裂紋起裂擴展研究［J］.應(yīng)用力學(xué)學(xué)報，2024，41（3）：524-535.

XU Haiduo， LI Hongyan， LIANG Bing， et al. Study on initiation and propagation of a single closed crack under compression-shear stress［J］.Chinese journal of applied mechanics，2024，41（3）：524-535.

文章編號：1000-4939（2024）03-0524-12

摘" 要：構(gòu)建了壓剪作用下單一閉合裂紋尖端應(yīng)力場方程，應(yīng)用最大周向應(yīng)力準則分析了翼裂紋開裂角，側(cè)壓系數(shù)k對應(yīng)力場和翼裂紋起裂擴展的影響，單軸壓縮下裂紋傾角α對翼裂紋起裂應(yīng)力的影響。通過單軸壓縮試驗及三軸文獻數(shù)據(jù)對比驗證了理論解的正確性。研究表明：開裂臨界相對尺寸λ→0或裂紋傾角α→45°，T應(yīng)力對翼裂紋開裂角的影響作用可以忽略；周向拉應(yīng)力隨側(cè)壓系數(shù)k的升高而降低，翼裂紋起裂所需軸向應(yīng)力增大；圍壓對翼裂紋的擴展存在抑制作用，且裂紋面有效剪應(yīng)力大于0時抑制作用更強；試驗觀察到預(yù)制裂紋尖端出現(xiàn)翼裂紋和少量的反翼裂紋和共面剪切裂紋，翼裂紋的開裂角在預(yù)測值附近上下波動；單軸壓縮，翼裂紋起裂應(yīng)力隨裂紋傾角的增大而增大；三軸壓縮，翼裂紋起裂應(yīng)力隨側(cè)壓系數(shù)的增大而增大，與模型預(yù)測的變化關(guān)系具有一致性。

關(guān)鍵詞：T應(yīng)力；有效剪應(yīng)力；翼裂紋；側(cè)壓系數(shù)；周向應(yīng)力；單軸壓縮試驗

中圖分類號：TU452" 文獻標志碼：A

DOI：10.11776/j.issn.1000-4939.2024.03.005

Study on initiation and propagation of a single closed crackunder compression-shear stress

XU Haiduo1，2，3，LI Hongyan1，2，LIANG Bing1，SUN Zhongxue2，3，4，HE Shi2，4，MO Yunlong2

（1.School of Mechanics and Engineering，Liaoning Technical University，123000 Fuxin，China;

2.China Coal ResearchInstitute，100013 Beijing，China;

3.State Key Laboratory of Coal Mining and Clean Utilization，100013 Beijing，China;4.School of Energy and Mining Engineering，China University of Mining and Technology（Beijing），100083 Beijing，China）

Abstract：

Field equation of single and closed flaw tip under compression-shear stress is constructed in this paper.The initiation angle of wing crack，the influence of side pressure coefficient k on stress field and initiation and propagation of wing crack and the influence of flaw angle α on initiation stress of wing crack under uniaxial compression are analyzed by using the maximum tangential stress criterion.The correctness of the theoretical solution is verified by uniaxial compression test and triaxial compression literature data comparison.The research shows

that the effect of T-stress on the initiation angle of wing crack can be ignored when the critical relative size of cracking λ approaches 0 or the flaw angle α approaches 45°.Hoop tensile stress decreases with the increase of side pressure coefficient k，and axial stress required for the initiation of wing crack increases.The confining pressure has an inhibitory effect on the propagation of wing cracks and the inhibitory effect is stronger when effective shear stress of the flaw surface is greater than 0.Wing cracks and a small number of anti-wing cracks and coplanar shear cracks are observed at the tip of the pre-existing flaw.The initiation angle of wing cracks are fluctuated around the predicted value.The initiation stress of wing crack increases with the increase of flaw angle under uniaxial compression，and it also increases with the increase of side pressure coefficient under triaxial compression，and they are consistent with the changing relationship predicted by the model in this paper.

Key words：T-stress; effective shear stress; wing crack; side pressure coefficient; hoop stress; uniaxial compression test

巖石內(nèi)部的微裂紋控制著巖石的斷裂破壞過程，巖石的強度及變形性質(zhì)都將隨著裂紋的孕育、起裂及擴展而發(fā)生變化。地下巖體通常是受壓的且在長年的地質(zhì)活動作用下會產(chǎn)生不同尺度的缺陷如裂隙、斷層等，這些含有缺陷的巖體在不同的受壓狀態(tài)下會產(chǎn)生不同類型的裂紋。因此研究含有裂紋的巖石在壓剪作用下的起裂、擴展，對于掌握地下圍巖破壞機理及預(yù)測圍巖工程質(zhì)量有著重要意義。

在預(yù)制裂紋試驗方面。BRACE等［1］對板狀玻璃預(yù)制單裂紋進行單軸壓縮試驗，結(jié)果僅觀察到翼裂紋并未看到其它裂紋，預(yù)制裂紋的最終擴展平行于最大主應(yīng)力方向。LAJTAI［2］使用類巖石試件預(yù)制單裂紋進行單軸壓縮試驗，試驗后觀察到翼裂紋、反翼裂紋和次生裂紋，其中翼裂紋最先產(chǎn)生。PETIT等［3］對PMMA試件和砂巖試件進行單、三軸壓縮試驗，單軸壓縮均有翼裂紋出現(xiàn)，而三軸壓縮很少出現(xiàn)翼裂紋。李銀平等［4］在大理石中預(yù)制單、雙裂紋，通過單軸壓縮試驗得到翼裂紋和反翼裂紋。李廷春等［5］使用類巖石材料預(yù)制單裂紋進行三軸壓縮試驗，試驗得到翼裂紋擴展緩慢，次生裂紋擴展明顯，預(yù)制裂隙的擴展受圍壓影響且擴展過程艱難。劉偉韜等［6］使用砂巖試件預(yù)制45°單裂紋進行單、三軸壓縮試驗，觀察到單軸壓縮以翼裂紋和反翼裂紋為主，三軸壓縮以反翼裂紋和次生裂紋為主。

針對上述現(xiàn)象，國內(nèi)外學(xué)者進行了相應(yīng)的理論分析。GRIFFITH［7］對材料斷裂現(xiàn)象進行了研究，從能量的觀點出發(fā)提出裂紋失穩(wěn)擴展條件，并建立了準脆性材料的斷裂強度與裂紋尺寸之間的關(guān)系。MCCLINTOCK［8］和HOEK等［9］考慮了裂紋面受壓閉合產(chǎn)生摩擦對GRIFFITH理論進行了修正。MUSKHELISHVILLI［10］認為對于彈性力學(xué)平面問題的應(yīng)力場可以用兩個復(fù)變量函數(shù)Φ（z）和Ω（z）來表示，該復(fù)變量函數(shù)通常稱為Kolossoff-Muskhelishvili應(yīng)力函數(shù)，這為研究裂紋尖端的應(yīng)力函數(shù)及應(yīng)力分布提供了理論依據(jù)。

WILLIAMS［11］對裂紋尖端應(yīng)力場進行級數(shù)展開得到奇異項、常數(shù)項和高階項。IRWIN［12］引入應(yīng)力強度因子和T應(yīng)力概念用以描述奇異項和常數(shù)項。WILLIAMS［13］和YUKIO等［14］預(yù)制裂紋PMMA試件進行單軸拉伸試驗，發(fā)現(xiàn)在最大周向應(yīng)力準則中引入T應(yīng)力可以有效減少預(yù)測值和試驗值之間差距。AYATOLLAHI等［15］將T應(yīng)力引入最大周向應(yīng)力準則，預(yù)測了Ⅰ型裂紋的破裂，表明T應(yīng)力對裂紋的起裂擴展有顯著影響。趙艷華等［16］用考慮T應(yīng)力的最大周向應(yīng)力準則分析了Ⅰ-Ⅱ復(fù)合型裂紋擴展，結(jié)果表明T應(yīng)力對裂紋擴展影響不可忽略，尤其對于Ⅱ型斷裂影響明顯。師俊平等［17］根據(jù)Muskhelishvilli理論結(jié)合裂紋面應(yīng)力邊界條件及Riemann-Hilbert解法分別給出了在偽集中力下不同類型裂紋面基本解，認為裂紋面狀態(tài)不影響KⅠ而影響KⅡ。劉紅巖［18］、唐世斌等［19］認為壓縮時應(yīng)考慮裂紋面間的摩擦作用，并提出當裂紋面有效剪應(yīng)力等于0時，裂紋尖端奇異性消失；當裂紋面有效剪應(yīng)力大于0時，裂紋尖端只存在KⅡ奇異項。LI等［20］利用Muskhelishvilli理論，從遠場應(yīng)力和裂紋面應(yīng)力邊界條件出發(fā)，推導(dǎo)了壓剪作用下含閉合裂紋平板的應(yīng)力函數(shù)，得到了單軸壓縮下考慮T應(yīng)力的裂紋尖端應(yīng)力分布。

由前人的研究工作可知：常規(guī)三軸壓縮試驗由于圍壓的作用很難在裂紋尖端產(chǎn)生翼裂紋，理論并不能很好地解釋現(xiàn)象；T應(yīng)力對裂紋的起裂擴展影響顯著，可以提高預(yù)測精度；裂紋面有效剪應(yīng)力取值對于應(yīng)力函數(shù)及場方程的影響還鮮有討論?；谝陨?點，本研究通過偽力法及希爾伯特邊值解法，給出了三軸壓縮下的Γ函數(shù)及相應(yīng)的裂紋尖端應(yīng)力函數(shù)和考慮T應(yīng)力的場方程，探討了裂紋面有效剪應(yīng)力大于0和等于0兩種情況下的翼裂紋開裂角、圍壓對周向應(yīng)力變化趨勢的影響、預(yù)制裂紋傾角對翼裂紋起裂應(yīng)力的影響，通過預(yù)制裂紋巖樣單軸壓縮試驗驗證翼裂紋開裂角及起裂強度，通過文獻數(shù)據(jù)與本研究預(yù)測值對比，驗證三軸壓縮下翼裂紋起裂強度，以期對后續(xù)含裂紋巖石壓剪破壞的理論分析有一定的參考價值。

1" 理論模型

1.1" 應(yīng)力強度因子的計算

含裂紋圓柱試件在軸壓和圍壓作用下裂紋面會發(fā)生閉合，圍壓均布作用于表面，簡化為平面模型，按照彈性力學(xué)平面問題求解，如圖1所示。設(shè)預(yù)制裂紋長為2a，邊緣受到均布軸向應(yīng)力σ和圍壓kσ，k為側(cè)壓系數(shù)且0≤k≤1，預(yù)制裂紋與水平方向夾角為α。

規(guī)定拉應(yīng)力為正，壓應(yīng)力為負，試件所受遠場應(yīng)力為

σx=-σ（sin2α+kcos2α）σ

SymboleB@y=-σ（cos2α+ksin2α）

τSymboleB@xy=-σ（1-k）sinαcosα

（1）

由偽力法并考慮裂紋面摩擦，可得遠場應(yīng)力在裂紋面上的作用力為

σn=σ（cos2α+ksin2α）（2a）

τe=0，τs≤τfτs-τf=σ1-ksinαcosα-

fσcos2α+ksin2α，τs＞τf（2b）

式中：σn、τe、τs和τf分別為作用在裂紋面上的壓應(yīng)力、有效剪應(yīng)力、切向偽力和摩擦力；f為裂紋面摩擦系數(shù)。

裂紋上下表面所受合力為壓應(yīng)力和有效剪應(yīng)力，所以裂紋面應(yīng)力邊界條件為

σ+y=-σn，τ+xy=-τe

σ-y=-σn，τ-xy=-τe（3）

式中：σy、τxy右上角的正負號表示裂紋上下表面，“+”代表裂紋上表面，“-”代表裂紋下表面。

根據(jù)Muskhelishvilli理論，裂紋尖端應(yīng)力場方程可由平面內(nèi)的應(yīng)力函數(shù)Φ（z）和Ω（z）表示，二者與裂紋附近應(yīng)力關(guān)系為［10］

σx+σy=4ReΦ（z）

σy-iτxy=Φ（z）+Ω（z）+（z-z）Φ′（z）

（4）

式中：z=x+iy為復(fù)變量；z為z的共軛。

當z分別從裂紋上下表面趨近于裂紋表面上一點t時，應(yīng)力函數(shù)滿足

Φ-（t）+Ω+（t）=σ+y-iτ+xy

Φ+（t）+Ω-（t）=σ-y-iτ-xy（5）

由式（3）、（5）得

［Φ（t）+Ω（t）］++［Φ（t）+Ω（t）］-=

-2σn+2iτe=Q

［Φ（t）+Ω（t）］+-［Φ（t）+Ω（t）］-=0（6）

其中，Q為復(fù)常數(shù)。

利用Riemann-Hilbert邊值問題解法［17］，可求通解

Φ（z）=4k0-Qz4" z2-a2+Q4-12?！洇福▃）=4k0-Qz4" z2-a2+Q4+12?！洌?）

式中

k0=Γ+12?！?/p>

（8）

對于常規(guī)三軸壓縮，Γ函數(shù)為

Γ=-14σ1+kΓ′=-12σ1-k［cos（2α）+isin（2α）］（9）

將式（2）、（6）、（8）、（9）代入式（7），得到常規(guī)三軸壓縮下裂紋尖端應(yīng)力函數(shù)。

當τs＜τf時，有

Φ（z）=-σ（1+k）4+iτs2z" z2-a2-1Ω（z）=-σ43-kcos2α+（3k-1）sin2α+

iτs2z" z2-a2+1（10）

同理，當τs＞τf時，有

Φ（z）=-σ1+k4+iτf2z" z2-a2-1

Ω（z）=-σ43-kcos2α+3k-1sin2α+

iτf2z" z2-a2-1+τs

（11）

以式（10）為例進行驗算，τs＜τf，τe=0，當z→∞時，由式（4）和（10）可得

σx+σy=4limz→SymboleB@

ReΦ（z）=" -σ1+k=σSymboleB@x+σSymboleB@y

σy-iτxy=Φ（z）+Ωz-+

（z-z）Φ′（z）=σSymboleB@y-iτSymboleB@xy（12）

滿足遠場應(yīng)力邊界條件。當z→t時，由式（5）和（10）得

limz→t-Φ（z）+limz→t+Ω（z）=" -σcos2α+ksin2α=σ+y-iτ+xylimz→t+Φ（z）+limz→t-Ω（z）=

-σcos2α+ksin2α=σ-y-iτ-xy（13）

滿足裂紋面應(yīng)力邊界條件。

應(yīng)力強度因子由下式計算［12］。

K=KⅠ-iKⅡ=2" 2πl(wèi)imz→±a" z±aΦ（z）

（14）

其中，KⅠ和KⅡ分別為Ⅰ、Ⅱ型裂紋尖端應(yīng)力強度因子。

將式（10）、（11）代入式（14）得裂紋面有效剪應(yīng)力等于0時，有

KⅠ=0KⅡ=-τsπa

（15）

同理裂紋面有效剪應(yīng)力大于0時，有

KⅠ=0KⅡ=-τfπa（16）

1.2" 裂尖應(yīng)力場的理論解

將式（10）、（11）兩種情況下的應(yīng)力函數(shù)分別帶入式（4）并極坐標化得σx、σy、τxy的理論解。裂紋面有效剪應(yīng)力等于0時，有

σx=σ1-ksinαcosα" πa" 2πrsinθ2·

2+cosθ2cos3θ2-

σsin2α+kcos2α+o（r）

σy=-σ1-ksinαcosαπa2πrsinθ2cosθ2·

cos3θ2-σ（cos2α+ksin2α

）+o（r）

τxy=-σ1-ksinαcosαπa2πrcosθ2·

1-sinθ2cos3θ2+o（r）

（17）

同理裂紋面有效剪應(yīng)力大于0時，有

σx=fσcos2α+ksin2α" πa" 2πrsinθ2·

2+cosθ2cos3θ2-σsin2α+kcos2α+orσy=-fσcos2α+ksin2απa" 2πrsinθ2cosθ2·

cos3θ2-σcos2α+ksin2α+o（r）τxy=-fσcos2α+ksin2α" πa" 2πrcosθ2·

1-sinθ2cos3θ2-σ1-ksinαcosα+fσcos2α+ksin2α+o（r）（18）

經(jīng)典斷裂力學(xué)認為，靠近裂紋尖端r→0，奇異項趨近于無窮大，從而忽略方程中的常數(shù)項和高階項，上兩式退化為純Ⅱ型剪切應(yīng)力場。此時裂尖應(yīng)力趨近于無窮，意味著裂紋體稍稍承受載荷便會在裂紋尖端產(chǎn)生無窮大的應(yīng)力集中，這表明經(jīng)典理論在裂尖處并不適用。目前研究表明裂紋體受力會在尖端產(chǎn)生塑性區(qū)，從而削弱應(yīng)力集中程度，用開裂過程區(qū)尺寸rc描述塑性區(qū)的大小。對于均質(zhì)脆性材料如玻璃、PMMA等，其開裂過程區(qū)尺寸rc很小，采用經(jīng)典的斷裂判據(jù)分析裂紋擴展問題誤差較??；對于巖石或混凝土材料，其由多種礦物集合而成，是典型的非均質(zhì)各向異性材料，不同裂尖距離包含的礦物類別差別較大，所以巖石材料的開裂過程區(qū)尺寸不能像玻璃等材料小到理想的范圍。因此對于巖石材料，開裂過程區(qū)尺寸rc較大與預(yù)制裂紋半長a差距減小，奇異項與常數(shù)項差距減小，常數(shù)項的占比升高，常數(shù)項的影響不可忽略。rc雖然與a的差距減小，但是rc依然遠小于a，高階項的影響忽略不計。用T應(yīng)力項Tx、Ty、Txy代替其中的常數(shù)項，并略去高階項得圖2所示裂尖應(yīng)力場，則

σx=-KⅡ" 2πrsinθ22+cosθ2cos3θ2+Txσy=KⅡ" 2πrsinθ2cosθ2cos3θ2+Tyτxy=KⅡ" 2πrcosθ21-sinθ2sin3θ2+Txy（19）

式（19）中裂紋面有效剪應(yīng)力等于0時，有

KⅡ=-σ1-ksinαcosα" πa=-τs" πaTx=-σsin2α+kcos2αTy=-σcos2α+ksin2αTxy=0（20）

同理裂紋面有效剪應(yīng)力大于0時，有

KⅡ=-fσ（cos2α+ksin2α）" πa=-τf" πaTx=-σ（sin2α+kcos2α）Ty=-σ（cos2α+ksin2α）Txy=-σ（1-k）sinαcosα+fσ（cos2α+ksin2α）（21）

2" 考慮T應(yīng)力的斷裂準則

最大周向應(yīng)力準則由ERDOGAN等［21］提出，該準則的基本假定是：①裂紋沿最大周向應(yīng)力σθmax的方向開裂；②當此方向的周向應(yīng)力達到臨界值時，裂紋失穩(wěn)擴展。這一準則經(jīng)常被用來預(yù)測翼裂紋的開裂，根據(jù)BOBET［22］、王飛陽等［23-24］的研究，含裂紋的巖石或類巖石材料受壓過程中會產(chǎn)生3種類型裂紋：翼裂紋、剪裂紋和拉剪復(fù)合裂紋。翼裂紋的產(chǎn)生是張拉效應(yīng)的結(jié)果，屬于巖石劈裂破壞；剪裂紋包括反翼裂紋和共面剪裂紋，二者是剪切作用的結(jié)果，斷口粗糙且細觀顆粒不平整，屬于典型的摩擦滑移破壞。翼裂紋和剪裂紋常出現(xiàn)在單、三軸試驗中。拉剪復(fù)合裂紋兼具翼裂紋和剪裂紋的特點，表現(xiàn)為前期翼裂后期裂紋偏折出現(xiàn)剪切特征。拉剪復(fù)合裂紋常出現(xiàn)在混凝土偏心加載試驗中。裂紋擴展示意圖如圖3所示。本研究僅討論翼裂紋的起裂擴展情況，根據(jù)最大周向應(yīng)力準則，有

σθθ=0， 2σθθ2lt;0

（22）

將式（19）坐標變換并代入式（22），得

σθ=-3KⅡ2" 2πrcsinθcosθ2+Txsin2θ+Tycos2θ-Txysin（2θ）σθθ=-3KⅡ4" 2πrccosθ23cosθ-1+Tx-Tysin（2θ）-2Txycos（2θ）=0

2σθθ2=3KⅡ16" 2πrcsinθ2+9sin3θ2+

2Tx-Tycos2θ+4Txysin2θlt;0

（23）

式中，rc為開裂過程區(qū)半徑，用來描述裂紋尖端塑性區(qū)大小。為便于分析，令A(yù)=Tx" πaKⅡ、B=Ty" πaKⅡ、C=Txy" πaKⅡ、λ=" 2rca，其中λ為裂紋尖端開裂臨界相對尺寸，表征裂尖塑性區(qū)邊界到裂尖的相對距離。

對式（23）進行無量綱處理，得

- "πaσθKⅡ=32λsinθcosθ2-Asin2θ-" Bcos2θ+Csin（2θ）

34λcosθ23cosθ-1-A-Bsin（2θ）+" 2Ccos（2θ）=03λ16sinθ2+9sin3θ2+2（A-B）cos（2θ）+" 4Csin（2θ）gt;0（24）

3" 開裂角及周向應(yīng)力分析

3.1" 裂紋面有效剪應(yīng)力等于0

裂紋面有效剪應(yīng)力等于0，由式（20）知Txy=0，C=0，代入式（24）中第2式可得

3cosθ23cosθ-1=4λA-Bsin（2θ）（25）

當cos（θ/2）≠sin（2θ）≠0時，λ=0或A=B（α=45°），解得cosθ=1/3，θ=70.5°，即不考慮T應(yīng)力時純Ⅱ型裂紋的開裂角。而實際上λ≠0，也就是說當λ→0或者α→45°時，T應(yīng)力對開裂角的影響作用可忽略。實際上從式（24）也可看出，當λ→0，開裂角方程中的奇異項趨近于無窮，從而T應(yīng)力的影響忽略不計，這與經(jīng)典斷裂理論相同，當α→45°，A→B，T應(yīng)力項趨近于0，T應(yīng)力對開裂角方程的影響也很??；當cos（θ/2）=sin（2θ）=0，θ=π為其中一個根，代入式（24）第1式可得θ=π時均一化周向應(yīng)力為壓應(yīng)力，不會產(chǎn)生翼裂紋。

由式（2）知在側(cè)壓系數(shù)k固定的情況下，裂紋面有效剪應(yīng)力的取值只與裂紋傾角α和裂紋面摩擦系數(shù)f有關(guān)。為了控制變量，本研究通過改變裂紋面摩擦系數(shù)f來改變裂紋面有效剪應(yīng)力的值并進行分析。裂紋面有效剪應(yīng)力等于0，A、B需滿足一定的條件，而A、B的取值又與k、α有關(guān)，經(jīng)計算，α=30°，f=0.6，λ=0.2，k從0增至0.5，裂紋面有效剪應(yīng)力均為0。圖4是α=30°，f=0.6，λ=0.2，均一化周向應(yīng)力隨θ的變化關(guān)系，k從0增大至0.5，圍壓增大，裂紋面壓應(yīng)力顯著增大而切向偽力顯著減小，T應(yīng)力顯著增大，因此周向拉應(yīng)力降低而周向壓應(yīng)力增大。從圖中可以看出圍壓增大，均一化峰值周向應(yīng)力逐漸降低，因此需要增加軸壓才能達到翼裂紋起裂強度，從而延緩了翼裂紋的起裂。

on homogenized hoop stress

圖5是側(cè)壓系數(shù)k=0和k=0.4均一化周向應(yīng)力隨θ的變化關(guān)系。從圖中可以看出，首先側(cè)壓系數(shù)k=0的峰值周向應(yīng)力高于k=0.4，其次側(cè)壓系數(shù)k=0（單軸加載）時，隨著λ的增加，在裂紋尖端塑性區(qū)內(nèi)峰值周向應(yīng)力均為拉應(yīng)力，軸壓增大周向拉應(yīng)力增大，達到翼裂紋的起裂強度從而產(chǎn)生翼裂紋。側(cè)壓系數(shù)k=0.4（三軸加載）時，隨著λ的增加，在裂紋尖端塑性區(qū)內(nèi)峰值周向應(yīng)力由拉應(yīng)力逐漸減小并轉(zhuǎn)變?yōu)閴簯?yīng)力，增加軸壓只能在小的裂尖范圍內(nèi)出現(xiàn)翼裂紋，翼裂紋在峰值周向應(yīng)力轉(zhuǎn)為壓應(yīng)力時不再擴展，施加圍壓將會抑制翼裂紋的延展。

3.2" 裂紋面有效剪應(yīng)力大于0

裂紋面有效剪應(yīng)力大于0，由式（21）知A、B、C均不為0，式（24）第2式寫成

3cosθ23cosθ-1+8Cλcos（2θ）

=4λA-Bsin（2θ）（26）

由式可知，λ=0，方程有兩個根，cos（θ/2）=0和cosθ=1/3。cos（θ/2）=0，θ=π為一根，代入式（24）第1式可得均一化周向應(yīng)力為壓應(yīng)力，不會產(chǎn)生翼裂紋。cosθ=1/3，θ=70.5°，易知滿足式（24）第3式，代入第1式為拉應(yīng)力。而λ≠0，當λ→0，θ→70.5°。

裂紋面有效剪應(yīng)力大于0，A、B、C需滿足一定的條件，取f=0.2，α=30°，λ=0.2，此時側(cè)壓系數(shù)k從0增至0.5裂紋面有效剪應(yīng)力均大于0。圖6所示為均一化周向應(yīng)力隨θ變化關(guān)系，與裂紋面有效剪應(yīng)力等于0的情況相比，周向壓應(yīng)力并非隨側(cè)壓系數(shù)k整體變化，而是一部分隨k的增大在減小，另一部分隨k的增大而增大，但是峰值周向拉應(yīng)力依然隨k的增大而減小，圍壓同樣會抑制翼裂紋的產(chǎn)生。

k on homogenized hoop stress

圖7是側(cè)壓系數(shù)k=0和k=0.4的均一化周向應(yīng)力隨θ的變化關(guān)系，可以看出首先側(cè)壓系數(shù)k由0增大到0.4，峰值周向應(yīng)力顯著減少，需要更大的外載才能達到開裂臨界值，圍壓抑制翼裂紋的產(chǎn)生，這與裂紋面有效剪應(yīng)力等于0的情況相同，其次側(cè)壓系數(shù)k=0.4時，隨著開裂相對臨界尺寸λ增大，與有效剪應(yīng)力等于0相比，在裂尖塑性區(qū)內(nèi)峰值周向應(yīng)力更快地由拉應(yīng)力轉(zhuǎn)為壓應(yīng)力，翼裂紋擴展的范圍變小，圍壓對翼裂紋擴展的抑制作用增強。

3.3" 裂紋傾角對周向應(yīng)力的影響

由前面的分析可知圍壓會抑制翼裂紋的起裂擴展，所以本研究僅分析了單軸壓縮k=0的情況。圖8是f=0.2，0.6，λ=0.2得到的均一化峰值周向應(yīng)力隨裂紋傾角的變化關(guān)系。均一化峰值周向應(yīng)力前期下降較慢后期下降較快，前期預(yù)制裂紋傾角較小，裂紋面切向偽力始終等于摩擦力，裂紋面有效剪應(yīng)力等于0，根據(jù)式（20）裂紋面切向偽力緩慢增大，Tx緩慢增大，Ty緩慢減小，由于Ty在裂尖應(yīng)力場中占比較大，所以周向應(yīng)力整體呈現(xiàn)緩慢下降趨勢；當預(yù)制裂紋傾角增大，tanα逐漸增大到大于裂紋面摩擦系數(shù)f時，裂紋面切向偽力大于摩擦力，裂紋面有效剪應(yīng)力大于0，根據(jù)式（21）裂紋面摩擦力迅速減小，Tx緩慢增大，Ty迅速減小，Txy變化不明顯，導(dǎo)致周向應(yīng)力整體下降較快。前期裂紋面有效剪應(yīng)力等于0，裂紋面切向偽力始終等于摩擦力，裂紋面摩擦系數(shù)f不影響T應(yīng)力的大小，故f對周向應(yīng)力沒有影響，曲線存在重合部分；后期裂紋面有效剪應(yīng)力大于0，相同裂紋傾角，f增大，裂紋面摩擦力和Txy增大，導(dǎo)致周向應(yīng)力升高?？傊?，隨著預(yù)制裂紋傾角的增大，均一化峰值周向應(yīng)力降低，相應(yīng)的翼裂紋起裂所需軸向應(yīng)力將增大。

記峰值周向應(yīng)力對應(yīng)的極角為θm，θm即為翼裂紋開裂角，代入式（24）第1式可得θm對應(yīng)的周向應(yīng)力σθm，根據(jù)最大周向應(yīng)力準則，若σθm達到了開裂臨界強度σθc則裂紋開裂，即

σθc=-3K*Ⅱc2" 2πrcsinθmcosθm2+Txsin2θm+

Tycos2θm-Txysin2θm=σθm

（27）

由此可得考慮T應(yīng)力的斷裂韌度為

K*Ⅱc=2Txsin2θm+Tycos2θm-Txysin2θm-σθc3sinθmcosθm2·

2πrc

（28）

將式（20）、（21）代入式（28）可得兩種情況下Ⅱ型裂紋的Ⅰ型斷裂韌度，二者顯然是不同的。

4" 試驗及驗證

4.1" 制作試樣

為了能夠清晰地得到翼裂紋及其擴展，設(shè)計了含閉合裂紋砂巖單軸壓縮試驗。巖芯取自內(nèi)蒙古巴彥高勒煤礦1214m深部巷道頂板砂巖，并在煤炭科學(xué)研究總院煤炭資源高效開采與潔凈利用國家重點實驗室進行二次加工，經(jīng)一系列工序制成直徑50mm，高度100mm的標準試件，砂巖參數(shù)如表1所示。試件中心使用高壓水刀切割一條穿透裂紋，裂紋長20mm，寬約1mm。為模擬閉合裂紋，在裂紋內(nèi)部緊緊插入一厚度為0.8mm的白鋼片。通過砂巖摩擦試驗測得白鋼片與砂巖之間的摩擦系數(shù)為0.21，該試驗利用自鎖現(xiàn)象和摩擦角原理［25］。將未加工的方形砂巖放置在一端可以轉(zhuǎn)動的白鋼板上，緩慢抬起白鋼板直到砂巖滑動，記錄此時斜面的傾角β，可由臨界滑動傾角β計算二者之間摩擦系數(shù)。單軸壓縮試驗試件分為5組，每組3個試件，裂紋傾角分別為0°、30°、45°、60°、90°，加工好的試件如圖9所示。

4.2" 試驗過程及斷裂形態(tài)

單軸壓縮試驗在煤炭科學(xué)研究總院國家重點實驗室TAW-2000電液伺服試驗機上進行。加載方式采用位移加載，加載速度為0.1mm/min，試驗過程中電腦記錄負荷、應(yīng)力、應(yīng)變和位移等數(shù)據(jù)，試驗后測量翼裂紋的開裂角。圖10為試件斷裂模式圖，斷裂主裂紋為翼裂紋、反翼裂紋和共面剪切裂紋。本試驗中翼裂紋產(chǎn)生居多，另外兩種裂紋出現(xiàn)較少，在預(yù)制裂紋尖端附近和試件表面還出現(xiàn)細小的次生裂紋，試件破壞表現(xiàn)為主裂紋和次生裂紋組合的斷裂帶。宋彥琦等［26］研究表明，當軸向載荷達到一定程度，預(yù)制裂紋尖端開始出現(xiàn)由于應(yīng)力集中導(dǎo)致的帶狀損傷區(qū)域，并逐漸凝聚形成裂紋雛形。0°試件預(yù)制裂紋中心有一圓孔，受壓過程中在圓孔頂部出現(xiàn)應(yīng)力集中產(chǎn)生表面裂紋，裂紋左、右端均有翼裂紋產(chǎn)生，翼裂紋開裂角分別為78°和90°。30°試件在裂紋端部出現(xiàn)2條翼裂紋和1條反翼裂紋，翼裂紋開裂角分別為89°和51°，翼裂紋以彎曲路徑擴展，擴展?jié)u近線朝向單軸載荷方向。45°試件出現(xiàn)2條翼裂紋，翼裂紋開裂角分別為42°和44°，上部翼裂紋在壓剪作用下在表面產(chǎn)生次生裂紋。60°試件產(chǎn)生2條翼裂紋和1條共面剪切裂紋，翼裂紋開裂角分別為35°和42°，在試件底部產(chǎn)生表面裂紋。90°試件預(yù)制裂紋上端出現(xiàn)翼裂紋、反翼裂紋以及橫向劈裂裂紋，下端出現(xiàn)劈裂裂紋和表面次生裂紋，翼裂紋開裂角約為3°和2°。主裂紋貫通形成完整斷裂面或巖塊局部破碎導(dǎo)致試件最終失去承載能力。根據(jù)郭奇峰［27］、HADI等［28］的研究，巖石材料單軸壓縮過程中翼裂紋出現(xiàn)較早，加載初期便會在裂紋尖端出現(xiàn)翼型擴展；而反翼裂紋和共面剪裂紋出現(xiàn)則較晚，在加載至峰值強度左右才會出現(xiàn)。

4.3" 強度特征及預(yù)測對比

選取僅產(chǎn)生翼裂紋試件的單軸壓縮數(shù)據(jù)進行處理，分別得到不同傾角單軸壓縮下的應(yīng)力-應(yīng)變曲線（圖11）。曲線經(jīng)過壓密、彈性、屈服、塑性破壞階段，有些試件破壞過快未記錄到殘余應(yīng)力階段。根據(jù)LI等［20］的研究，應(yīng)力-應(yīng)變曲線中出現(xiàn)明顯的應(yīng)力降或應(yīng)變軟化意味著裂紋的擴展，即初次下降點所對應(yīng)的應(yīng)力為試件初始起裂應(yīng)力。本研究僅選取產(chǎn)生翼裂紋的試件，意味著曲線中出現(xiàn)的初次應(yīng)力降或者應(yīng)變屈服軟化為翼裂紋的初始起裂，后續(xù)的每一次應(yīng)力降代表其他翼裂紋的起裂或原有翼裂紋的擴展貫通。

with different inclination angles

統(tǒng)計各傾角試件應(yīng)力-應(yīng)變曲線中初次應(yīng)力降點或屈服點的應(yīng)力值，預(yù)制0°、30°、45°、60°和90°裂紋翼裂紋起裂應(yīng)力分別為峰值應(yīng)力的52%、78%、80%、88%、90%，預(yù)制裂紋傾角越大翼裂紋起裂應(yīng)力越接近峰值強度。對不同預(yù)制裂紋傾角起裂應(yīng)力進行擬合，擬合直線如圖12所示。預(yù)制裂紋傾角越大，翼裂紋的起裂應(yīng)力越大，與本研究的預(yù)測趨勢具有一致性。

對裂紋體施加相同的軸向壓力，裂尖周向應(yīng)力越大起裂應(yīng)力越小，周向應(yīng)力與起裂應(yīng)力負相關(guān)。收集文獻［29］和文獻［6］中的三軸壓縮試驗應(yīng)力數(shù)據(jù)。為了方便說明將計算所得均一化起裂應(yīng)力放大400倍，得到圖13所示對比圖。由圖可得，側(cè)壓系數(shù)較小，起裂應(yīng)力隨k幾乎線性增加，符合文獻［6］的變化趨勢；側(cè)壓系數(shù)較大，起裂應(yīng)力隨k增長變快，類似指數(shù)形式，符合文獻［29］的變化趨勢。后期隨著k的增加，翼裂紋起裂應(yīng)力將超過試件抗壓強度，從另一方面說明在較大的圍壓下，翼裂紋不會產(chǎn)生。

對于巖石或混凝土等脆性材質(zhì)， SCHIMIDT［30］根據(jù)最大周向應(yīng)力準則給出計算開裂過程區(qū)半徑rc的公式，即

rc=12πKⅠcσt2

（29）

式中：KⅠc為Ⅰ型斷裂韌度；σt為抗拉強度。根據(jù)表1可求得本次試驗所用砂巖的斷裂過程區(qū)半徑rc約為2.52mm，根據(jù)式（24）可求得砂巖翼裂紋開裂角。預(yù)測值和試驗值如圖14所示，可以看出翼裂紋開裂角隨預(yù)制裂紋傾角的增大而降低，且試驗值基本在預(yù)測值處上下波動。根據(jù)擬合曲線，裂紋傾角在0°附近預(yù)測值偏大，90°附近預(yù)測值偏小，中間角度預(yù)測值與試驗值吻合較好。

5" 結(jié)" 論

本研究通過偽力法和復(fù)變理論給出三軸壓縮的Γ函數(shù)及相應(yīng)的裂紋尖端應(yīng)力函數(shù)和考慮T應(yīng)力的場方程，根據(jù)最大周向應(yīng)力準則分析了裂紋面有效剪應(yīng)力等于0和大于0兩種情況下翼裂紋開裂角、側(cè)壓系數(shù)對周向應(yīng)力的影響和單軸壓縮預(yù)制裂紋傾角對翼裂紋起裂應(yīng)力的影響。通過預(yù)制裂紋砂巖單軸壓縮試驗得到翼裂紋開裂角和起裂應(yīng)力，將三軸壓縮的預(yù)測值與文獻值進行對比。分析理論和試驗結(jié)果可得如下結(jié)論。

1）裂紋面有效剪應(yīng)力等于0，開裂臨界相對尺寸λ→0或者預(yù)制裂紋傾角α→45°，T應(yīng)力對開裂角的影響作用可以忽略；裂紋面有效剪應(yīng)力大于0，僅當開裂臨界相對尺寸λ→0，T應(yīng)力對開裂角的影響較小，這種情況與忽略T應(yīng)力的經(jīng)典斷裂理論相同。

2）分析了側(cè)壓系數(shù)對周向應(yīng)力的影響，結(jié)果表明周向拉應(yīng)力隨側(cè)壓系數(shù)的增大而減小，翼裂紋開裂所需的載荷隨側(cè)壓系數(shù)的增大而增大。

3）三軸壓縮時，隨著開裂相對臨界尺寸λ的增大，峰值周向應(yīng)力由拉轉(zhuǎn)為壓，從而翼裂紋在擴展的過程中受到抑制。裂紋面有效剪應(yīng)力大于0時峰值周向應(yīng)力由拉應(yīng)力轉(zhuǎn)為壓應(yīng)力會更早地發(fā)生，圍壓對翼裂紋擴展的抑制作用增強。

4）單軸壓縮試驗得到的翼裂紋開裂角與本研究預(yù)測值較吻合，翼裂紋起裂應(yīng)力隨裂紋傾角的增大而增大。通過文獻數(shù)據(jù)對比，三軸壓縮下翼裂紋起裂應(yīng)力隨側(cè)壓系數(shù)的增大而增大，與本研究的預(yù)測趨勢具有一致性。

參考文獻：

［1］" BRACE W F，BOMBOLAKIS E G.A note on brittle crack growth in compression［J］.Journal of geophysical research，1963，68（12）：3709-3713.

［2］" LAJTAI E Z.A theoretical and experimental evaluation of the Griffith theory of brittle fracture［J］.Tectonophysics，1971，11（2）：129-156.

［3］" PETIT J P，BARQUINS M.Can natural faults propagate under mode II conditions？［J］.Tectonics，1988，7（6）：1243-1256.

［4］" 李銀平，王元漢，陳龍珠，等.含預(yù)制裂紋大理巖的壓剪試驗分析［J］.巖土工程學(xué)報，2004，26（1）：120-124.

LI Yinping，WANG Yuanhan，CHEN Longzhu，et al.Experimental research on pre-existing cracks in marble under compression［J］.Chinese journal of geotechnical engineering，2004，26（1）：120-124（in Chinese）.

［5］" 李廷春，呂海波.三軸壓縮載荷作用下單裂隙擴展的CT實時掃描試驗［J］.巖石力學(xué)與工程學(xué)報，2010，29（2）：289-296.

LI Tingchun，L Haibo.CT real-time scanning tests on single crack propagation under triaxial compression［J］.Chinese journal of rock mechanics and engineering，2010，29（2）：289-296（in Chinese）.

［6］" 劉偉韜，申建軍.含單裂紋真實巖石試件斷裂模式的力學(xué)試驗研究［J］.巖石力學(xué)與工程學(xué)報，2016，35（6）：1182-1189.

LIU Weitao，SHEN Jianjun.Experimental study of propagation mode of crack in real rock specimens with a single crack［J］.Chinese journal of rock mechanics and engineering，2016，35（6）：1182-1189（in Chinese）.

［7］" GRIFFITH A A.The phenomena of rupture and flow in solids［J］.Philosophical Transactions of the Royal Society of London.series a，containing papers of a mathematical or physical character，1921，221（593）：163-198.

［8］" MCCLINTOCK F A，WALSH J B.Friction on Griffith crack in rocks under pressure［R］.1962：1015-1021.

［9］" HOEK E，BIENIAWSKI Z T.Brittle fracture propagation in rock under compression［J］.International journal of fracture mechanics，1965，1（3）：137-155.

［10］MUSKHELISHVILI N I.Some basic problems of the mathematical theory of elasticity：fundamental equations，plane theory of elasticity，torsion and bending［M］.Holland：Noordhoff Press，1953.

［11］WILLIAMS M L.On the stress distribution at the base of a stationary crack［J］.Journal of applied mechanics，1957，24（1）：109-114.

［12］IRWIN G R.Analysis of stresses and strains near the end of a crack traversing a plate［J］.Journal of applied mechanics，1957，24（3）：361-364.

［13］WILLIAMS J G，EWING P D.Fracture under complex stress — the angled crack problem［J］.International journal of fracture mechanics，1972，8（4）：441-446.

［14］YUKIO U，KAZUO I，TETSUYA Y，et al.Characteristics of brittle fracture under general combined modes including those under bi-axial tensile loads［J］.Engineering fracture mechanics，1983，18（6）：1131-1158.

［15］AYATOLLAHI M R，PAVIER M J，SMITH D J.Mode I cracks subjected to large T-stresses［J］.International journal of fracture，2002，117（2）：159-174.

［16］趙艷華，陳晉，張華.T應(yīng)力對Ⅰ-Ⅱ復(fù)合型裂紋擴展的影響［J］.工程力學(xué)，2010，27（4）：5-12.

ZHAO Yanhua，CHEN Jin，ZHANG Hua.Influence of T-stress on crack propagation for Ⅰ-Ⅱ mixed mode loading［J］.Engineering mechanics，2010，27（4）：5-12（in Chinese）.

［17］師俊平，解敏，王靜.無限大平面中斜裂紋的壓剪斷裂分析［J］.工程力學(xué)，2006，23（12）：59-62.

SHI Junping，XIE Min，WANG Jing.Compression and shear fracture analysis of slanting cracks in infinite planes［J］.Engineering mechanics，2006，23（12）：59-62（in Chinese）.

［18］劉紅巖.考慮T應(yīng)力的巖石壓剪裂紋起裂機理［J］.巖土工程學(xué)報，2019，41（7）：1296-1302.

LIU Hongyan.Initiation mechanism of cracks of rock in compression and shear considering T-stress［J］.Chinese journal of geotechnical engineering，2019，41（7）：1296-1302（in Chinese）.

［19］唐世斌，黃潤秋，唐春安.T應(yīng)力對巖石裂紋擴展路徑及起裂強度的影響研究［J］.巖土力學(xué)，2016，37（6）：1521-1529，1549.

TANG Shibin，HUANG Runqiu，TANG Chunan.Effect of T-stress on crack growth path in rock and fracture strength［J］.Rock and soil mechanics，2016，37（6）：1521-1529，1549（in Chinese）.

［20］LI Q，YANG Q，JIA JC，et al.Experimental research on crack propagation and failure in rock-type materials under compression［J］.Electronic journal of geotechnical engineering，2008，13（1）：1-13.

［21］ERDOGAN F，SIH G C.On the crack extension in plates under plane loading and transverse shear［J］.Journal of basic engineering，1963，85（4）：519-525.

［22］BOBET A.The initiation of secondary cracks in compression［J］.Engineering fracture mechanics，2000，66（2）：187-219.

［23］王飛陽，黃宏偉，張東明，等.帶裂縫服役混凝土結(jié)構(gòu)力學(xué)性能的多尺度模擬方法［J］.建筑結(jié)構(gòu)學(xué)報，2019，40（12）：155-162.

WANG Feiyang，HUANG Hongwei，ZHANG Dongming，et al.Multi scale simulation method of mechanical behaviors of existing concrete structure with crack［J］.Journal of building structures，2019，40（12）：155-162（in Chinese）.

［24］王飛陽，黃宏偉.盾構(gòu)隧道襯砌結(jié)構(gòu)裂縫演化規(guī)律及其簡化模擬方法［J］.巖石力學(xué)與工程學(xué)報，2020，39（增刊1）：2902-2910.

WANG Feiyang，HUANG Hongwei.Crack evolution law of shield tunnel lining and its simplified simulation method［J］.Chinese journal of rock mechanics and engineering，2020，39（S1）：2902-2910（in Chinese）.

［25］郭少華.巖石類材料壓縮斷裂的實驗與理論研究［D］.長沙：中南大學(xué)，2003.

［26］宋彥琦，李名，王曉，等.基于高速攝影的單預(yù)制裂紋大理巖加卸載試驗［J］.中國礦業(yè)大學(xué)學(xué)報，2014，43（5）：773-781.

SONG Yanqi，LI Ming，WANG Xiao，et al.Observation of marble with single pre-existing crack using high-speed photography under loading and unloading conditions［J］.Journal of China University of Mining amp; Technology，2014，43（5）：773-781（in Chinese）.

［27］郭奇峰，武旭，蔡美峰，等.預(yù)制裂隙花崗巖的裂紋起裂機理試驗研究［J］.煤炭學(xué)報，2019，44（增刊2）：476-483.

GUO Qifeng，WU Xu，CAI Meifeng，et al.Crack initiation mechanism of pre-existing cracked granite［J］.Journal of China Coal Society，2019，44（S2）：476-483（in Chinese）.

［28］HAERI H，SHAHRIAR K，MARJI M F，et al.Cracks coalescence mechanism and cracks propagation paths in rock-like specimens containing pre-existing random cracks under compression［J］.Journal of Central South University，2014，21（6）：2404-2414.

［29］肖桃李，李新平，賈善坡.深部單裂隙巖體結(jié)構(gòu)面效應(yīng)的三軸試驗研究與力學(xué)分析［J］.巖石力學(xué)與工程學(xué)報，2012，31（8）：1666-1673.

XIAO Taoli，LI Xinping，JIA Shanpo.Triaxial test research and mechanical analysis based on structure surface effect of deep rock mass with single fissure［J］.Chinese journal of rock mechanics and engineering，2012，31（8）：1666-1673（in Chinese）.

［30］SCHMIDT R A.A microcrack model and its significance to hydraulic fracturing and fracture toughness testing［C］//The 21st U.S.Symposium on Rock Mechanics （USRMS）.Richardson，Texas，USA：OnePetro，1980：ARMA-80-0581.

（編輯" 李坤璐）