唐 治，潘一山，朱小景，崔乃鑫

（遼寧工程技術(shù)大學 礦山安全技術(shù)裝備研究院；力學與工程學院，遼寧 阜新 123000）

煤礦巷道圍巖應力大小和規(guī)律是巷道支護方式選取的重要依據(jù)之一。圍巖應力大小不僅與煤礦采深、側(cè)壓系數(shù)等有關(guān)，還與巷道斷面形狀等相關(guān)（即使相同圍巖條件下，圍巖應力分布規(guī)律和圍巖變形破壞規(guī)律也因巷道斷面形狀不同而不同）。直墻拱形斷面巷道的斷面由下部分矩形和上部分拱形組成，長期實踐證明，直墻拱形斷面巷道具有較好的穩(wěn)定性，所以服務(wù)年限較長的巷道一般均采用直墻拱形斷面巷道。了解和掌握巷道圍巖應力分布規(guī)律對合理選擇巷道支護方式具有重要的理論意義和實際應用價值。對于常規(guī)的圓形、橢圓形等巷道的圍巖應力可以采用Cauchy積分法或冪級數(shù)法方便解出［1－3］，但復雜巷道圍巖應力公式需借助 復 變 函 數(shù) 彈 性 理 論［4－6］及 映 射 函 數(shù)［7－9］。 朱 大 勇等［10］求解了矩形斷面圍巖應力彈性解，趙凱等［11］利用多角形法得出了矩形硐室圍巖應力，王潤富［12］、劉金高等［13］求解了梯形孔口的應力，湯澄波等［14］、祝江鴻［15］分別利用復變函數(shù)法求解了天幕線拱形圍巖應力和兩個表示斷面圍巖應力的解析函數(shù)通式，但前提是能給出復雜斷面映射函數(shù)。針對直墻拱形巷道圍巖應力研究較少，對其他復雜巷道研究大多也只給出復雜巷道應力解析的隱函數(shù)，未給出巷道應力分布情況，也未對巷道斷面高寬比和側(cè)壓系數(shù)對其影響規(guī)律進行深入分析。

1 模型建立

圖1 模型圖Fig.1 diagram of model

圖2 單位圓外域Fig 2 The unit circle outland

2 直墻拱形巷道圍巖應力解析

2.1 保角變換

首先將直墻拱形巷道圍巖映射為單位圓形巷道圍巖，即把圖1中的六邊形的外域部分映射為圖2的單位圓外域部分。映射后A1、A2、A3、B3、B2、B1分別a1t、a2t、a3t、b3t、b2t、b1t與 對 應，如 圖 2。則 有a1t＝ie－iκ1π，a2t＝ie－iκ2π，a3t＝ie－iκ3π，b1t＝ie－iκ1π，b2t＝ie－iκ2π，b3t＝ie－iκ3π。其中K1、K2、K3由六邊形的尺寸確定。

圖3 保角變換Fig.3 Conformal transformation

2.2 由邊界條件確定

2.3 確定φ0（ζ）、ψ0（ζ）

2.5 確定應力分量

3 直墻拱形巷道圍巖應力分析

3.1 直墻拱形巷道邊界應力分布規(guī)律

巷道邊界處，ρ＝1，ζ＝eiφ，σρ＝τρφ＝0。取a＝1 m，b＝2 m，當斷面高寬c0比分別取0.5、1、1.5時，得直墻拱形巷道周邊應力分布規(guī)律分別如圖4（a）、（b）、（c），－90°～90°對應底板中點到頂板中點圍巖邊界。

圖4 應力分布規(guī)律Fig.4 stress distribution law

從圖4（a）可得：巷道斷面高寬比為0.5情況下，側(cè)壓系數(shù)大于1時圍巖邊界環(huán)向應力有相同變化趨勢，即底板中點到直墻底部再到直墻頂部的邊界環(huán)向應力先增后減，在直墻底部底角處區(qū)域出現(xiàn)較大應力集中；直墻頂部到拱形頂板中點的邊界環(huán)向應力先減后增再減。側(cè)壓系數(shù)小于1時，圍巖邊界環(huán)向應力也有相同變化趨勢，即底板中點和頂板中點附近均出現(xiàn)拉應力，直墻底部到直墻頂部的邊界環(huán)向應力較為恒定；直墻頂部到拱形頂板中點的邊界環(huán)向應力先增后減然后變?yōu)槔瓚?。應力集中區(qū)域主要集中在直墻底部底角處和拱形頂板中點附近。

從圖4（b）可得：巷道斷面高寬比為1情況下，側(cè)壓系數(shù)大于1時，圍巖邊界環(huán)向應力有相同變化趨勢，即底板中點附近均出現(xiàn)拉應力，底板中點到直墻底部再到直墻頂部的邊界環(huán)向應力先由拉應力變?yōu)閴簯?，然后增加后減小；直墻頂部到拱形頂板中點的邊界環(huán)向應力先減后增再減。側(cè)壓系數(shù)小于1時，圍巖邊界環(huán)向應力也有相同變化趨勢，即頂板中點附近均出現(xiàn)拉應力，底板中點到直墻底部再到直墻頂部的邊界環(huán)向應力先由逐漸減小；直墻頂部到拱形頂板中點的邊界環(huán)向應力先增后減然后變?yōu)槔瓚Α袇^(qū)域主要集中在底板中部和拱形頂板中點附近，且拱形頂板應力集中系數(shù)小于底板。

從圖4（c）可得：巷道斷面高寬比為1.5情況下，巷道圍巖應力分布規(guī)律與巷道斷面高寬比為1情況基本相同，不同之處在于拱形頂板應力集中系數(shù)大于底板，圍巖應力分布比巷道斷面高寬比為1時較好。

3.2 直墻拱形巷道沿水平線的應力分布規(guī)律

取φ＝0，ζ＝ρ為直墻拱形巷道水平線位置，由x＝R（1＋c1ρ2＋c3ρ4）／ρ可將曲線坐標表示的應力分量表達式轉(zhuǎn)換為直角坐標表示。取a＝1 m，b＝2 m，斷面高寬比c0分別取0.5、1、1.5時，可得沿x軸圍巖應力分布規(guī)律，如圖5（a）、（b）、（c）。

圖5 直角坐標下巷道周邊應力分布規(guī)律Fig.5 Stress distribution rule of perimeter under rectangular coordinate

從圖5可得：1）不同巷道斷面高寬比下，直墻拱形巷道沿水平線的應力分布規(guī)律基本相同。2）側(cè)壓系數(shù)大于1時，不同巷道斷面高寬比的環(huán)向應力均隨至巷道邊界距離增大而迅速增大，在距離巷道邊界2～4 m后達到穩(wěn)定；側(cè)壓系數(shù)小于1時，不同巷道斷面高寬比的環(huán)向應力均隨至巷道邊界距離增大而先增大后減小，在距離巷道邊界1 m左右達到最大值，然后較小并在距離巷道邊界2～4 m后達到穩(wěn)定。3）側(cè)壓系數(shù)大于1時，不同巷道斷面高寬比的徑向應力均隨至巷道邊界距離增大而先減小后變?yōu)槔瓚θ缓笤黾?，在距離巷道邊界4～6 m后達到穩(wěn)定；側(cè)壓系數(shù)小于1時，不同巷道斷面高寬比的徑向應力均隨至巷道邊界距離增大而增大，在距離巷道邊界2～4 m后達到穩(wěn)定。4）直墻拱形巷道邊界3 m范圍內(nèi)出現(xiàn)了剪應力，剪應力隨至巷道邊界距離增大而迅速減小。當側(cè)壓系數(shù)大于1時，最大剪應力隨側(cè)壓系數(shù)增加而增大；當側(cè)壓系數(shù)小于1時，最大剪應力隨側(cè)壓系數(shù)增加而減小。

4 結(jié)論

1）采用保角變換，應用復變函數(shù)彈性理論推導了直墻拱形巷道圍巖應力分布的解析表達式。

2）不同巷道斷面高寬比、側(cè)壓系數(shù)下，直墻拱形巷道圍巖應力集中區(qū)域均主要集中在直墻底部底角處、拱形頂板中點附近和底板中部3個位置。巷道斷面高寬比一定情況下，側(cè)壓系數(shù)大于1時，圍巖邊界環(huán)向應力有相同變化趨勢；采用巷道斷面高寬比小于1較有利于巷道穩(wěn)定；側(cè)壓系數(shù)小于等于1時，圍巖邊界環(huán)向應力也有相同變化趨勢；采用巷道斷面高寬比大于1較有利于巷道穩(wěn)定。

3）道斷面高寬比對直墻拱形巷道沿水平線的應力分布規(guī)律影響較小。側(cè)壓系數(shù)大于1時，巷道環(huán)向應力均隨至巷道邊界距離增大而迅速增大，徑向應力均隨至巷道邊界距離增大而先減小后變?yōu)槔瓚θ缓笤黾?，最大剪應力隨側(cè)壓系數(shù)增加而增大；側(cè)壓系數(shù)小于1時，巷道環(huán)向應力均隨至巷道邊界距離增大而先增大后減小，徑向應力均隨至巷道邊界距離增大而增大，最大剪應力隨側(cè)壓系數(shù)增加而減小。

［1］王明斌，李術(shù)才，李樹忱，等.圓形隧道圍巖附加應力場的解析解答［J］.巖土力學，2006，27（Sup）：207-210.

Wang M B，Li S C，Li S C，et al.Analytical solution of subsidiary stress field for circular tunnel［J］.Rock and Soil Mechanics，2006，27（Sup）：207-210.（in Chinese）.

［2］盧文超，仲政，王旭.淺埋隧道圍巖應力場的解析解［J］.力學季刊，2003，24（1）：50-54.

Lu W C，Zhong Z，Wang X.Analytical solution for stress field in surrounding rocks of shallow tunnel［J］.Journal of Chinese Quarterly Mechanics，2003，24（1）：50-54.（in Chinese）

［3］蔡曉鴻，蔡勇斌，蔡勇平，等.二向不等圍壓和內(nèi)壓作用下橢圓形洞室的計算［J］.地下空間與工程學報，2008，4（3）：453-459.

Cai X H，Cai Y B，Cai Y P，et al.Computation of elliptic tunnel under the combined action of two-dimensional unequal adjoining rock pressure and internal pressure［J］. Chinese Journal Underground Space and Engineering，2008，4（3）：453-459.（in Chinese）

［4］Muskhelishvili N I.Some basic problems of the mathematical theory of elasticity：Fundamental equations，plane theory of elasticity，torsison and bending［M］.Groningen：P Noordhoof，1953.

［5］陳子蔭.圍巖力學分析中的解析方法［M］.北京：煤炭工業(yè)出版社，1994.

Chen Z Y.Analytic method of mechanical analysis for the surrounding rock［M］.Beijing：Coal Industry Publishing House，1994.（in Chinese）

［6］Zhang Z Z，Sun Y Z.Analytical solution for a deep tunnel with arbitrary cross section in a transversely isotropic rock mass ［J］.International Journal of Rock Mechanics and Mining Sciences，2011，48 （8） ：1359-1363.

［7］房營光，孫鈞.地面荷載下淺埋隧道圍巖的黏彈性應力和變形分析［J］.巖石力學與工程學報，1998，17（3）：239-247.

Fang Y G，Sun J.Viscoelastic stress and deformation analysis of shallow tunnels under the load on the ground surrounding［J］.Chinese Journal of Rock Mechanics and Engineering，1998，17（3）：239-247.（in Chinese）

［8］朱大勇，錢七虎，周早生，等.復雜形狀洞室映射函數(shù)的新解 法 ［J］.巖 石 力 學 與 工 程 學 報，1999，18（3）：279-282.

Zhu D Y，Qian Q H，Zhou Z S，et al.New method for calculating mapping function of opening with complex shape ［J］.Chinese Journal of Rock Mechanics and Engineering，1999，18（3）：279-282.（in Chinese）

［9］皇甫鵬鵬，伍法權(quán)，郭松峰，等.基于邊界點搜索的洞室外域映射函數(shù)求解法［J］.巖石力學，2011，32（5）：1418-1424.

Huangfu P P，Wu F Q，Guo S F，et al.A new method for calculating mapping function of external area of cavern with arbitrary shape based on searching points on boundary［J］.Rock and Soil Mechanics，2011，32 （5）：1418-1424.（in Chinese）

［10］朱大勇，錢七虎，周早生，等.復雜形狀洞室圍巖應力彈性解析分析［J］.巖石力學與工程學報，1999，18（4）：402-404.

Zhu D Y，Qian Q H，Zhou Z S，et al.Elastic solution to stresses of rock masses around openings with complex shape ［J］.Chinese Journal of Rock Mechanics and Engineering，1999，18（4）：402-404.（in Chinese）

［11］趙凱，劉長武，張國良.用彈性力學的復變函數(shù)法求解矩形硐室周邊應力［J］.采礦與安全工程學報，2007，24（3）：361-365.

Zhao K，Liu C W，Zhang G L.Solution for perimeter stresses of rocks around a rectangular chamber using the complex function of elastic mechanics［J］.Journal of Mining ＆Safety Engineering，2007，24（3）：361-365.（in Chinese）

［12］王潤富.彈性力學的復變函數(shù)計算機解［J］.河海大學學報，1991，19（2）：71-74.

Wang R F.Computer solutions for complex function in elasticity［J］.Journal of Hohai University，1991，19（2）：71-74.（in Chinese）

［13］劉金高，王潤富.馬蹄形孔口和梯形孔口的應力集中問題［J］.巖土工程學報，1995，17（5）：57-63.

Liu J G，Wang R F.The stress concentration of U-shaped and ladder-shaped holes［J］.Chinese Journal of Geotechnical Engineering，1995，17（5）：57-63.（in Chinese）

［14］湯澄波，范廣勤.高地應力區(qū)天幕線拱形硐室斷面的解法［J］.煤炭工程，1989（1）：16-18.

Tang C B，F(xiàn)an G Q.Determination of the mapping function for the exterior do-main of a non-circular opening by means of the multiplication of tree absolutely convergent series ［J］.Coal Engineering，1989（1）：16-18.（in Chinese）

［15］祝江鴻.隧洞圍巖應力復變函數(shù)分析法中的解析函數(shù)求解［J］.應用數(shù)學和力學，2013，34（4）：346-354.

Zhu J H.Rock for caverns with the complex variable theory［J］.Applied Mathematics and Mechanics，2013，34（4）：346-354.（in Chinese）

（編輯王秀玲）