基于三參數(shù)強(qiáng)度準(zhǔn)則的高強(qiáng)混凝土井壁力學(xué)特性分析

李雪梅，榮傳新，程　樺

(安徽理工大學(xué)土木建筑學(xué)院， 安徽淮南232001)

?

基于三參數(shù)強(qiáng)度準(zhǔn)則的高強(qiáng)混凝土井壁力學(xué)特性分析

李雪梅，榮傳新，程樺

(安徽理工大學(xué)土木建筑學(xué)院， 安徽淮南232001)

摘要：為了分析高強(qiáng)混凝土井壁的力學(xué)特性，采用三參數(shù)強(qiáng)度準(zhǔn)則推導(dǎo)出了高強(qiáng)混凝土立井井壁極限承載力理論解以及彈性區(qū)和塑性區(qū)應(yīng)力與半徑和荷載之間的解析表達(dá)式，并對計算結(jié)果進(jìn)行了實驗的驗證。計算分析表明：在彈性區(qū)，高強(qiáng)混凝土井壁徑向壓應(yīng)力σr隨半徑r的增大而增大，環(huán)向壓應(yīng)力σ(θ)隨半徑r的增大而減少；在塑性區(qū)，徑向壓應(yīng)力σr和環(huán)向壓應(yīng)力σ(θ)均隨半徑r的增大而增大。當(dāng)井壁內(nèi)半徑為4.0 m、厚度為1.0 m、混凝土強(qiáng)度等級為C60時，井壁極限承載力為22.87 MPa，井壁厚度每增加0.1 m，井壁極限承載力增加2.8 MPa左右。同時，混凝土井壁的環(huán)向壓應(yīng)力σ(θ)達(dá)到164.38 MPa，是混凝土立方體單軸抗壓強(qiáng)度的2.7倍左右，表明考慮了混凝土多軸強(qiáng)度影響的三參數(shù)強(qiáng)度準(zhǔn)則更適用于高強(qiáng)混凝土井壁力學(xué)特性分析，為高強(qiáng)混凝土井壁結(jié)構(gòu)的設(shè)計計算提供了的理論參考。

關(guān)鍵詞：三參數(shù)強(qiáng)度準(zhǔn)則；高強(qiáng)混凝土井壁；極限承載力；井壁應(yīng)力

因淺部煤炭資源的日益枯竭，在山東和安徽等產(chǎn)煤大省的煤炭開發(fā)已逐漸轉(zhuǎn)向深部，新井建設(shè)時煤礦立井井筒需要穿過的沖積層越來越厚，如在建的山東巨野煤田萬福礦和即將建設(shè)安徽展溝煤田口孜西礦井穿過的第三和第四系表土層達(dá)600～800 m。在如此深厚的沖積層中建井必須大幅度提高立井井壁極限承載力，以抵抗作用在其上巨大的水土壓力，為此，在井壁結(jié)構(gòu)設(shè)計中采用了C60～C80的高強(qiáng)混凝土。關(guān)于混凝土井壁結(jié)構(gòu)的研究已有較多的報道[1-4]，但關(guān)于高強(qiáng)混凝土井壁結(jié)構(gòu)的力學(xué)特性研究多以實驗研究為主[5-8]，并根據(jù)井壁模型實驗結(jié)果得到了井壁極限承載力的實驗回歸公式，同時采用Mohr-Coulomb準(zhǔn)則分析了井壁結(jié)構(gòu)的力學(xué)特性[9]。在壓力不太大的情況下，采用Mohr-Coulomb準(zhǔn)則對其進(jìn)行分析是比較符合實際，但當(dāng)壓力較大的情況下則采用三參數(shù)強(qiáng)度準(zhǔn)則[10-11]能更好地反映高強(qiáng)混凝土井壁的力學(xué)特性，因為該準(zhǔn)則考慮了混凝土在多軸應(yīng)力下的力學(xué)性能[12]。深厚沖積層凍結(jié)法施工的煤礦立井井壁一般采用雙層井壁結(jié)構(gòu)型式，凍結(jié)井內(nèi)層井壁中的混凝土一般處于雙向或三向應(yīng)力狀態(tài)[13]。本文針對深凍結(jié)井內(nèi)層井壁受力特點，采用了符合混凝土強(qiáng)度特性的三參數(shù)強(qiáng)度準(zhǔn)則[14]，通過分析計算，推導(dǎo)出高強(qiáng)混凝土立井井壁結(jié)構(gòu)極限承載力以及井壁彈性區(qū)和塑性區(qū)應(yīng)力的解析表達(dá)式，旨在為井壁的結(jié)構(gòu)設(shè)計提供更多的理論參考。

1混凝土三參數(shù)強(qiáng)度準(zhǔn)則

三參數(shù)強(qiáng)度準(zhǔn)則的破壞面為外凸的光滑曲面，其子午線為拋物線，在負(fù)靜水壓力軸的方向是開口的，描述準(zhǔn)則的三個力學(xué)特征值為單軸抗拉強(qiáng)度σl、單軸抗壓強(qiáng)度σy、雙軸抗壓強(qiáng)度σ2y[10]，三參數(shù)準(zhǔn)則通常以ρ、ξ、θ表示其應(yīng)力組合，其中：

(1)

式中：

I1=σ1+σ2+σ3=3σm，

J3=s1s2s3=(σ1-σm)(σ2-σm)(σ3-σm)。

三參數(shù)強(qiáng)度準(zhǔn)則的屈服面如圖1所示，破壞面上任一點M(ρ,ξ,θ)滿足：

M(ρ,ξ,θ)=ρ2+A1(θ)ρ+A2(θ)ξ+A3(θ)=0，

(2)

式中：

其中：

Ky=σy/σl，

K2y=σ2y/σl。

圖1　三參數(shù)強(qiáng)度準(zhǔn)則

2高強(qiáng)混凝土井壁結(jié)構(gòu)力學(xué)特性分析

圖2　高強(qiáng)混凝土井壁的力學(xué)模型Fig.2　Mechanical model of high-strength concrete shaft lining

混凝土井壁在外壓作用下，其應(yīng)力分量σr，σθ,σz為主應(yīng)力，且有σr>σz>σθ，σr<0，σθ<0。井筒的徑向和軸向應(yīng)力均為壓應(yīng)力。設(shè)煤礦立井井筒的內(nèi)半徑為a，外半徑為b，在外壓作用下井壁處于彈塑性狀態(tài)，彈塑性分界半徑為rp，如圖2所示。

主應(yīng)力符號規(guī)定，拉為“+”、壓為“-”，且σ1>σ2>σ3。

2.1三參數(shù)準(zhǔn)則應(yīng)用于井壁極限承載力分析中的公式推導(dǎo)

煤礦立井井壁可以作為平面軸對稱問題進(jìn)行分析，徑向應(yīng)力和切向應(yīng)力僅是r的函數(shù)，與θ無關(guān)，σr=σr(r)，σθ=σθ(r)。假設(shè)中間主應(yīng)力為σz=β(σr+σθ)/2，0≤β≤1，β為中間主應(yīng)力系數(shù)，在井壁塑性區(qū)，取β=1 ，則σ1=σr，σ3=σθ，σ2=(σr+σθ)/2，從而有：

(3)

(4)

cos3θ=0，

(5)

由式(3)和式(4)可得：

(6)

(7)

將式(5)代入式(2)，得：

ρ2+C1ρ+C2ξ+C3=0，

(8)

因此，求解混凝土井壁受外壓p作用的問題時應(yīng)滿足平衡條件式(7)、強(qiáng)度準(zhǔn)則式(8)以及邊界條件：

當(dāng)p很小時，整個立井井筒處于線彈性階段，各應(yīng)力分量為：

代入式(3)和式(4)，可得：

(9)

將式(9)代入強(qiáng)度條件式(8)，可得：

(10)

由此可知，當(dāng)外壓p逐漸增大時，井筒內(nèi)壁(r=a)最先滿足強(qiáng)度條件(或稱進(jìn)入塑性)，相應(yīng)的外壓為彈性極限壓力，記為pe。把r=a代入式(10)，可得：

(11)

當(dāng)p>pe時，在井壁內(nèi)緣出現(xiàn)塑性區(qū)(a≤r≤rp)，并且隨著荷載p的增加，塑性區(qū)逐漸向外擴(kuò)展。若彈塑性分界半徑為rp，由于軸對稱，在交界面上只能有徑向壓力q，即σr|r=rp=-q。在塑性區(qū)內(nèi)，應(yīng)力分量應(yīng)滿足平衡方程與強(qiáng)度準(zhǔn)則。對于混凝土，有Ky?1和K2y?1，因此，有C2>0，由式(8)可得：

(12)

將式(12)代入式(7)，有：

(13)

將式(13)積分，并利用邊界條件σr|r=a=0，整理可得：

(14)

把式(11)和r=a代入式(9)第一式，可得：

(15)

把C1、C2和C3的值代入式(15)得到ρ(a)，同時聯(lián)立式(14)、式(12)和式(6)可確定塑性區(qū)內(nèi)的應(yīng)力分量為：

(16)

式中，ρ(r)可根據(jù)式(14)在代入各參數(shù)后利用MATLAB求得，其余參數(shù)同前。

當(dāng)混凝土立井井筒出現(xiàn)塑性區(qū)時，在彈性區(qū)(rp≤r≤b)內(nèi)可視為在r=rp處受內(nèi)壓q作用，在r=b處受外壓p作用的彈性混凝土立井井筒，其應(yīng)力分量為：

(17)

將式(17)代入(3)和式(4)，有：

在r=rp處，有：

(18)

式(18)中q和rp均未知，可以由彈性與塑性交界處徑向應(yīng)力相等的條件確定二者之間的關(guān)系。在塑性區(qū)r=rp處求σr時，利用式(6)第一式及式(12)，有：

(19)

在彈性區(qū)的r=rp處，材料剛進(jìn)入塑性狀態(tài)，應(yīng)力分量應(yīng)滿足強(qiáng)度準(zhǔn)則，將式(18)代入(8)式，可得：

(20)

把r=rp和ρ(a)的值代入式(14)得ρ(rp)，再聯(lián)立式(19)和式(20)得井壁極限承載力pp與塑性區(qū)半徑rp的關(guān)系式為：

(21)

2.2理論計算與模型實驗結(jié)果對比分析

由文獻(xiàn)[9]中高強(qiáng)鋼筋混凝土井壁結(jié)構(gòu)實驗結(jié)果可知，把模型實驗的參數(shù)代入式(21)，此時rp=b，得到的高強(qiáng)混凝土井壁極限承載力如表1所示。由表1可知，計算結(jié)果和實驗結(jié)果相差很小，其誤差在±5.0%左右，表明三參數(shù)強(qiáng)度準(zhǔn)則不僅參數(shù)少、函數(shù)形式簡單，由此推導(dǎo)出的理論計算公式能夠很好地反映高強(qiáng)混凝土井壁結(jié)構(gòu)的力學(xué)特性。

表1　混凝土井壁塑性極限承載力理論分析與實驗結(jié)果的比較

2.3井壁極限承載力分析

以C60混凝土為例，井壁內(nèi)半徑a=4.0 m，外半徑b分別取4.8 m、4.9 m、5.0 m、5.1 m，5.2 m時，幾何尺寸的大小對井壁極限承載力的影響如表3和圖4所示。表3中，每個塑性區(qū)半徑所在列的第一個數(shù)據(jù)為每個半徑所對應(yīng)的井壁塑性極限承載力。由表3和圖4可見，不同厚度的井壁在同一塑性區(qū)半徑處的極限承載力相差很大，且井壁厚度每增加0.1 m，井壁的塑性極限承載力增加2.8 MPa左右。

表2　不同強(qiáng)度等級的混凝土所對應(yīng)的井壁極限承載力分析

圖3　強(qiáng)度等級不同時，極限承載力pp與

塑性區(qū)半徑rp的關(guān)系曲線圖

Fig.3Curve of ultimate bearing capacityppand

plastic radiusrpunder different concrete strength

圖4井壁厚度不同時，極限承載力pp與

塑性區(qū)半徑rp的關(guān)系曲線圖

Fig.4Curve of ultimate bearing capacityppand plastic

radiusrpunder different thickness of the shaft lining

表3井壁厚度不同時，不同位置所對應(yīng)的井壁極限承載力分析

Tab.3The analysis of ultimate bearing capacity of shaft liningpp, when different thickness of the shaft lining

混凝土井壁厚度/m井壁極限承載力/MParp=4.8mrp=4.9mrp=5.0mrp=5.1mrp=5.2m厚度每增加0.1m承載力變化值0.817.38—————0.919.9520.09———2.711.022.3822.7522.87——2.781.124.6625.2625.6125.75—2.881.226.8127.6328.2028.5828.692.94

2.4井壁應(yīng)力分析

取井壁內(nèi)半徑a=4 m，外半徑b=5 m，混凝土強(qiáng)度等級為C60。當(dāng)井壁的塑性區(qū)半徑rp=4.0 m時，對應(yīng)井壁的彈性極限承載力為11.16 MPa；rp=4.3 m時，井壁極限承載力為17.02 MPa；rp=4.5 m時，井壁極限承載力為19.85 MPa；rp=5.0 m時，對應(yīng)井壁的極限承載力為22.87 MPa。井壁外圍分別承受這四種荷載作用下，井壁應(yīng)力σr、σθ和半徑r的關(guān)系曲線如圖5和圖6所示。由此可見，在彈性區(qū)，徑向壓應(yīng)力σr隨半徑r的增大而增大，環(huán)向壓應(yīng)力σθ隨半徑r的增大而減少，均呈拋物線性分布；在塑性區(qū)，壓應(yīng)力σr、σθ均隨半徑r的增大而增大。當(dāng)高強(qiáng)混凝土井壁的極限承載力等于22.87 MPa時，由于考慮了混凝土多軸強(qiáng)度的影響，混凝土井壁的環(huán)向壓應(yīng)力σθ達(dá)到了164.38 MPa，是混凝土立方體單軸抗壓強(qiáng)度的2.7倍左右。

取r的值分別為4.3 m、4.5 m、4.7 m和5.0 m(記為r1)，計算并繪出不同位置處壓應(yīng)力σr、σθ和荷載Pp的關(guān)系曲線如圖7和圖8所示。結(jié)果表明，隨著荷載的增大，處于彈性區(qū)(a≤r≤r1)的井壁混凝土壓應(yīng)力σr、σθ逐漸增加，而當(dāng)荷載增加達(dá)到塑性半徑所對應(yīng)的極限荷載后，該位置的井壁混凝土徑向和環(huán)向應(yīng)力則保持不變，表明混凝土塑性區(qū)應(yīng)力僅與塑性半徑rp和混凝土本身的性能有關(guān)。當(dāng)荷載為17.02 MPa時，井壁的塑性區(qū)半徑rp=4.3 m，徑向應(yīng)力σr為-5.48 MPa，環(huán)向應(yīng)力σθ為-94.18 MPa，但隨著荷載的增大，在r=4.3 m處，井壁混凝土徑向和環(huán)向應(yīng)力均保持不變。

圖5pp取不同值時，徑向應(yīng)力σr

與半徑r的關(guān)系曲線圖

Fig.5Curve of radial stressσrand

radiusr, when differentppvalue

圖6pp取不同值時，環(huán)向應(yīng)力σθ

與半徑r的關(guān)系曲線圖

Fig.6Curve of hoop stressσθand

radiusr, when differentppvalue

圖7r取不同值時，徑向應(yīng)力σr和

荷載p的關(guān)系曲線圖

Fig.7Curve of radial stressσrand

loadp, when different r value

圖8r取不同值時，環(huán)向應(yīng)力σθ和

荷載p的關(guān)系曲線圖

Fig.8Curve of hoop stressσθand

loadp, when differentrvalue

綜上所述，高強(qiáng)混凝土立井井壁結(jié)構(gòu)的壓應(yīng)力σr、σθ在塑性區(qū)(a≤r≤rp)僅和半徑有關(guān)，在彈性區(qū)(rp≤r≤b)除了和半徑有關(guān)外還與所受荷載有關(guān)。由于三參數(shù)強(qiáng)度準(zhǔn)則考慮了混凝土多軸強(qiáng)度的影響，極大地提高了高強(qiáng)混凝土井壁的環(huán)向應(yīng)力，且實驗結(jié)果也表明三參數(shù)強(qiáng)度準(zhǔn)則用于高強(qiáng)混凝土井壁力學(xué)特性的分析是安全可靠的。

3結(jié)論

①文中基于三參數(shù)強(qiáng)度準(zhǔn)則推導(dǎo)出了高強(qiáng)混凝土井壁極限承載力理論解，以及彈性區(qū)和塑性區(qū)應(yīng)力與半徑和荷載之間的解析表達(dá)式。高強(qiáng)混凝土井壁極限承載力計算結(jié)果和實驗結(jié)果相差很小，其誤差均在±5.0%左右。由此可見，采用三參數(shù)強(qiáng)度準(zhǔn)則對高強(qiáng)混凝土井壁結(jié)構(gòu)進(jìn)行力學(xué)特性分析是可行的。

②凝土強(qiáng)度等級和幾何尺寸對井壁極限承載力的影響較大，當(dāng)井壁內(nèi)半徑為4.0 m、井壁厚度為1.0 m、混凝土強(qiáng)度等級為C60時，井壁極限承載力為22.87 MPa；井壁厚度每增加0.1 m，井壁極限承載力增加2.8 MPa左右。同時，混凝土井壁的環(huán)向壓應(yīng)力σθ達(dá)到164.38 MPa，是混凝土立方體單軸抗壓強(qiáng)度的2.7倍左右，表明三參數(shù)強(qiáng)度準(zhǔn)則更適用于高強(qiáng)混凝土井壁力學(xué)特性分析。另外，井壁混凝土強(qiáng)度等級從C60至C80每提高5 MPa，井壁極限承載力增加1.7 MPa左右。

③彈性區(qū)，高強(qiáng)混凝土井壁徑向壓應(yīng)力σr隨半徑r的增大而增大，環(huán)向壓應(yīng)力σθ隨半徑r的增大而減少，均呈拋物線性分布；在塑性區(qū)，徑向壓應(yīng)力σr和環(huán)向壓應(yīng)力σθ均隨半徑r的增大而增大。由此可見，井壁結(jié)構(gòu)的應(yīng)力σr、σθ在塑性區(qū)(a≤r≤rp)僅與半徑有關(guān)，在彈性區(qū)(rp≤r≤b)不僅與半徑有關(guān)外，而且還與外荷載有關(guān)。

參考文獻(xiàn)：

[1]姚直書，桂建剛，程樺，等.內(nèi)層鋼板高強(qiáng)鋼筋混凝土復(fù)合井壁數(shù)值模擬[J]. 廣西大學(xué)學(xué)報(自然科學(xué)版), 2010,35(1):35-38.

[2]徐志偉，邵鵬，商翔宇，等.深厚表土地壓誤差對凍結(jié)壁厚度設(shè)計的影響[J]. 中國礦業(yè)大學(xué)學(xué)報, 2014,43(4):606-611.

[3]LEE S W, JUNG J W, NAM S W, et al. The influence of seepage forces on ground reaction curve of circular opening[J]. Tunnelling and Underground Space Technology, 2006, 22(1) : 28-38.

[4]王衍森，黃家會，李金華，等 凍結(jié)井外壁高強(qiáng)混凝土的早期強(qiáng)度增長規(guī)律研究[J]. 中國礦業(yè)大學(xué)學(xué)報, 2008,37(5):595-598.

[5]榮傳新，程樺.凍結(jié)井鋼筋混凝土弧形板井壁力學(xué)特性研究[J]. 巖石力學(xué).2006, 27(2):193-198.

[6]YAO Z S, YANG J J, SUN W R.Experimental study on sliding shaft lining mechanical mechanisms underground subsidence conditions[J]. Journal of Coal Science and Engineering, 2003, 9 (1):95-99.

[7]姚直書，程樺，楊俊杰.深表土中高強(qiáng)鋼筋混凝土井壁力學(xué)性能的試驗研究[J]. 煤炭學(xué)報，2004，29(2)：167-171.

[8]姚直書，薛維培，宋海清，等.富水松軟巖層凍結(jié)法鑿井井壁結(jié)構(gòu)試驗研究[J]. 廣西大學(xué)學(xué)報(自然科學(xué)版), 2014, 39 (2):231-236.

[9]榮傳新, 王秀喜, 程樺，等.深厚沖積層高強(qiáng)鋼筋混凝土井壁力學(xué)特性研究[J]. 巖石力學(xué)與工程學(xué)報，2008，27(S1): 2841-2847.

[10]徐秉業(yè)，劉信聲.應(yīng)用彈塑性力學(xué)[M]. 北京:清華大學(xué)，1995.

[11]張遠(yuǎn)高，陸明萬.混凝土的三參數(shù)破壞準(zhǔn)則[J]. 上海交通大學(xué)學(xué)報, 1990,24(5,6):62-88.

[12]過鎮(zhèn)海，王傳志.多軸應(yīng)力下混凝土的強(qiáng)度破壞和破壞準(zhǔn)則研究[J]. 土木工程學(xué)報，1991,24(3): 1-13.

[13]YAO Z S, CHANG H, RONG C X.Research on stress and strength of high strength reinforced concrete drilling shaft lining in thick top soils[J]. Journal of China University Mining & Technology,2007,17(3):0432-0435.

[14]李朝弟.結(jié)構(gòu)的安定性與極限分析[D]. 北京: 清華大學(xué), 1991.

[15]KUPFER H, SHINOZUKA M, SCHELLER G, et al.Behavior of concrete under biaxial stress[J]. American Concrete Institute, 1969, 66(8): 656-666.

[16]中華人民共和國國家標(biāo)準(zhǔn)編寫組.GB50010-2010 混凝土結(jié)構(gòu)設(shè)計規(guī)范[S]. 北京：中國建筑工業(yè)出版社，2010.

(責(zé)任編輯唐漢民裴潤梅)

Mechanical properties of high-strength concrete shaft lining based on three-parameter strength criterion

LI Xue-mei, RONG Chuan-xin, CHENG Hua

(School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan 232001, China)

Abstract:In order to analyze the mechanical properties of high-strength concrete shaft lining, the three-parameter strength criterion was used to obtain the theoretical ultimate bearing capacity of the high-strength concrete shaft lining and the analytical expressions among stress, radius and the loads in both elastic and plastic regions. The calculated results were verified by experimental results. The results show that the radial stress σr increases while the circumferential stress σ(θ) decreases with the increase of radius r in the elastic region. The radial stress σr and the circumferential stress σ(θ) increase with the increase of radius r in the plastic zone. The ultimate bearing capacity of the shaft lining pp is 22.87 MPa when the thickness of the shaft lining is 1.0 m and the concrete strength is C60. The ultimate bearing capacity of the shaft lining increases by about 2.8 MPa when the thickness increases by 0.1 m. And the circumferential stress of the concrete shaft lining reaches 164.38 MPa, which is about 2.7 times of the uniaxial compressive strength of concrete cubes and demonstrates that the three-parameter strength criterion, considering the influence of multi-axial strength of concrete, is more applicable to analyze the mechanical properties of high-strength concrete shaft lining. It provides a theoretical reference for structure design and calculation of high strength-concrete shaft lining.

Key words:three-parameter strength criterion; high strength concrete shaft lining; ultimate bearing capacity; shaft lining stress

中圖分類號:TD262

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

文章編號：1001-7445(2016)02-0308-09

doi:10.13624/j.cnki.issn.1001-7445.2016.0308

通訊作者：榮傳新(1968—)，男，安徽六安人，安徽理工大學(xué)教授，博士生導(dǎo)師，博士；E-mail: chxrong@aust.edu.cn。

基金項目：國家自然科學(xué)基金資助項目(51374010，51474004)

收稿日期：2015-11-01；

修訂日期：2016-01-21

引文格式：李雪梅，榮傳新，程樺.基于三參數(shù)強(qiáng)度準(zhǔn)則的高強(qiáng)混凝土井壁力學(xué)特性分析[J].廣西大學(xué)學(xué)報(自然科學(xué)版)，2016，41(2)：308-316.