陳慶發(fā)，鐘瓊英，王 拓

(廣西大學(xué) 資源與冶金學(xué)院，廣西 南寧530004)

對巖石地下工程穩(wěn)定性進行監(jiān)測與預(yù)報，是保證工程設(shè)計、施工科學(xué)合理和安全生產(chǎn)的重要措施［1］。收斂觀測是一種常用的圍巖位移監(jiān)測方法，這種方法簡單易行，通過測量開挖巷道若干點的相對位移的變化，從而掌握圍巖位移的變形發(fā)展動態(tài)，進而判斷施工的安全程度以及檢驗圍巖是否穩(wěn)定，支護是否合理［2-3］。

在巷道圍巖收斂觀測研究領(lǐng)域，盛建龍［4］采用建立一種神經(jīng)網(wǎng)絡(luò)模型的方法對采場巷道進行收斂預(yù)測;許夢國等［5］通過對監(jiān)測結(jié)果的分析，確定影響巷道收斂的因素，對監(jiān)測結(jié)果中出現(xiàn)的一些特殊情況進行分析;付國彬［6］開展了考慮巖石破裂后強度衰減和體積膨脹條件下巷道圍巖位移及收斂觀測技術(shù)的研究。這些研究大多是對收斂觀測結(jié)果進行分析，未探討觀測點的布置問題。趙同彬［7］曾基于遺傳算法對巷道位移反分析進行了研究，其研究結(jié)果對于觀測點布置起到了一定的借鑒作用。本文采用復(fù)變函數(shù)的方法探討均質(zhì)巖體巷道三點收斂觀測法的測點準(zhǔn)確布置問題。

1 均質(zhì)巖體巷道收斂觀測點布置的局限性

文獻［2］指出，均質(zhì)巖體巷道圍巖收斂觀測一般分為周邊收斂量測和拱頂下沉量測，其測點與測線的布置方案如圖1 所示。

對于均質(zhì)巖體而言，由于巷道形狀的對稱性，在巷道頂?shù)装逯行狞c處的布點位置無疑是正確的，但巷道兩幫或圓弧段上的觀測點位置，大多由工程師根據(jù)個人經(jīng)驗確定，具有較大的隨機性。這也是很多大型工程無法及時、準(zhǔn)確地預(yù)測出片幫破壞現(xiàn)象的一個重要原因。

2 觀測點布置的復(fù)變函數(shù)解析

圖1 收斂觀測測點測線布置圖Fig.1 Arrangement of points and lines about convergence observation

2.1 復(fù)變函數(shù)方法

以復(fù)數(shù)作為自變量的函數(shù)就叫做復(fù)變函數(shù)，而與之相關(guān)的理論就是復(fù)變函數(shù)理論［8-9］。復(fù)變函數(shù)的優(yōu)點在于能夠把復(fù)雜的幾何構(gòu)型單連通域映射為較為簡單的幾何構(gòu)型單連通域，從而在這個區(qū)域內(nèi)利用解析函數(shù)的理論和邊界條件確定待求的解析函數(shù)。

20 世紀(jì)以來，復(fù)變函數(shù)理論被廣泛應(yīng)用于彈性與彈塑性力學(xué)領(lǐng)域，Verruijt［10］采用復(fù)變函數(shù)方法分別研究不考慮巖土體自重情形下無襯砌淺埋隧道變形導(dǎo)致的地表變形問題，Strack 等［11］研究了考慮巖土體自重情況下由于隧道構(gòu)造浮力效應(yīng)所導(dǎo)致的地表變形問題，蔚立元等［12］利用復(fù)變函數(shù)開展了含圓孔半平面體的彈性分析，王立忠等［13］利用復(fù)變函數(shù)分析了盾構(gòu)隧道施工引起的地基變形，呂愛鐘等［14］應(yīng)用復(fù)變函數(shù)方法對地下隧洞進行了詳細的力學(xué)分析。

2.2 收斂觀測點形變特征

根據(jù)文獻［2］的分析，巷道兩幫或圓弧段上收斂觀測點應(yīng)當(dāng)具有豎直方向位移為零而水平方向位移最大的形變特征，確定巷道兩幫或圓弧段上收斂觀測點位置即是尋找滿足上述形變的特征點。

為避免重力作用帶來的不必要干擾，理論分析過程中將圍巖重力分量忽略，僅考慮原巖地應(yīng)力條件下的巖體開挖。

2.3 收斂觀測點在復(fù)數(shù)平面上的位置

式中，Ck(k=0，1，2，3，…，n)取決于巷道斷面形狀，文獻［15］給出了Ck的確定方法。

設(shè)Aj為巷道周邊任意一點，其坐標(biāo)為(rj，αj)，在單位圓上相應(yīng)的映射點為Bj，Bj的極坐標(biāo)為(1，βj)。

巷道斷面內(nèi)水平方向收斂觀測點是巷道兩幫豎直位移為零的點，即u=0。同時，所需找的點在巷道的邊界上，巷道周邊對應(yīng)于ξ 平面上ρ=1 的邊界，則可知EE=0，F(xiàn)F=0。

這樣，式(1)可簡化為:

即:

將AA、PP、UU、SS 分別代入式(4)，可得:

根據(jù)地下隧洞力學(xué)分析的復(fù)變函數(shù)方法，在解決實際問題時，只要n ≥4 計算得到巷道邊界位移與實際位移相差較小，且n 取4 以及取4 以上相差不大，因此，在保證精確的情況下選取n=4 可保證計算方便。

在n 取4 時，式(4)可以化為:

運用三角函數(shù)二倍角以及三倍角公式可將式(6)化為:

將cosθ 看作一整體，可用解一元四次方程的費拉里方法來計算cosθ，這里只需求其實數(shù)根。

令a=8D，b=4C，c=2B-8D，d=A-3，e=D-B，則式(7)求解過程如下:

選取方程的解中值在(-1，1)范圍內(nèi)的值為cosθ 的值。

在ξ 平面cosθ 以及sinθ 求出之后，且巷道的邊界上ρ=1，則收斂觀測點在ξ 平面上的位置就很明確了，其極坐標(biāo)可表示為(1，θ)。

2.4 收斂點在直角平面坐標(biāo)系的位置

尋找到實際巷道上的點，需將ξ 復(fù)數(shù)平面上的點轉(zhuǎn)化到z 直角平面上。根據(jù)復(fù)變函數(shù)的分析方法，其轉(zhuǎn)換公式為:

式中，ξ=cosθ+isinθ。

將公式(8)展開，可得:

在z 平面的x 坐標(biāo)(豎直方向坐標(biāo))即是式(9)右邊的虛部，即:

在z 平面的y 坐標(biāo)(水平方向的坐標(biāo))即是式(9)右邊的實部，即:

所尋找的巷道收斂觀測點位于巷道的兩幫及或圓弧段上，在極坐標(biāo)系中，θ 值取值區(qū)間為［30°，150°］，K 值取值區(qū)間為［-0.87 ～+0.87］，省去K 的高次方，簡化后的x，y 值如下:

其次，根據(jù)復(fù)變函數(shù)分析原理，確定出直角坐標(biāo)系的原點位置，如圖2 所示。圖中h1=h2，x 軸為巷道的中心對稱軸。

求出巷道兩幫或圓弧段收斂點的豎軸坐標(biāo)后，對應(yīng)于巷道兩幫或圓弧段的點即為所求收斂點，即圖2 中的M1、M2兩點。

圖2 直角坐標(biāo)系原點與觀測點Fig.2 The origin of rectangular coordinate system and observation points

以巷道底板為起點，換算成觀測點高度，即為:

式中，H 為巷道高度。

3 討 論

①依據(jù)平面應(yīng)變問題的假設(shè)條件，巷道的軸向必須為主應(yīng)力方向，對于軸向不是主應(yīng)力方向，需進行坐標(biāo)換算后再進行求解。

②巷道復(fù)變函數(shù)豎直位移分量的公式中n 取4，誤差精度已滿足日常巖石工程收斂觀測技術(shù)要求，對于其他形狀的巷道，如需要更加精確計算，還需對誤差進一步分析研究。

③工程中采用本文理論研究成果，還需根據(jù)地應(yīng)力參數(shù)和巷道形狀參數(shù)，對觀測點坐標(biāo)值進行編程求解出具體位置。

4 結(jié) 語

本文采用復(fù)變函數(shù)的方法，通過對均質(zhì)巖體巷道三點收斂觀測點布置，巷道兩幫與圓弧段豎直方向的位移進行解析，確定出豎向位移為零而水平位移最大的形變特征點，從而解決了三點收斂觀測法測點準(zhǔn)確布置的理論依據(jù)問題。

［1］ 楊強，劉耀儒，常強，等.結(jié)構(gòu)變形穩(wěn)定與控制理論及在巖土工程中的應(yīng)用［J］.工程力學(xué)，2010，27(12):61-87.

［2］ 李造鼎.巖體測試技術(shù)［M］.北京:冶金工業(yè)出版社，1993.

［3］ 羅立強，王衛(wèi)軍，屈延嗣，等.深井巖巷分步聯(lián)合支護技術(shù)應(yīng)用［J］.廣西大學(xué)學(xué)報:自然科學(xué)版，2013，38(1):228-234.

［4］ 盛建龍.基于神經(jīng)網(wǎng)絡(luò)模型的采場巷道收斂預(yù)測研究［J］.武漢科技大學(xué)學(xué)報:自然科學(xué)版，2005，28(2):172-174.

［5］ 許夢國，伍佑倫，葉義成.無底柱分段崩落法回采巷道收斂的影響因素分析［J］.礦業(yè)研究與開發(fā)，2001，21(1):12-14.

［6］ 付國彬.巷道圍巖破裂范圍與位移的新研究［J］.煤炭學(xué)報，1995，2(3):304-310.

［7］ 趙同彬.基于遺傳算法的巷道位移反分析研究［J］.巖土力學(xué)，2004，25(增1):107-109.

［8］ 鐘玉泉.復(fù)變函數(shù)論［M］.北京:高等教育出版社，2004.

［9］ 陸明方，羅學(xué)富.彈性理論基礎(chǔ)(下冊)［M］.2 版.北京:清華大學(xué)出版社，2001.

［10］VERRUIJT A.Complex variable solution for a deforming circular tunnel in an elastic half-plane［J］.International Journal for Numerical and Analytical Methods in Geomechanics，1997，21(2):77-89.

［11］STRACK O E，VERRUIJT A.Complex variable solution for a deforming buoyant tunnel in a heavy elastic half-plane［J］.International Journal for Numerical and Analytical Methods in Geomechanics，2002，26(12):1235-1252.

［12］蔚立元，韓立軍，陳曉鵬，等.含圓孔半平面體的彈性分析及其工程應(yīng)用［J］.工程力學(xué)，2013，30(7):167-172.

［13］王立忠，呂學(xué)金.復(fù)變函數(shù)分析盾構(gòu)隧道施工引起的地基變形［J］.巖土工程學(xué)報，2007，29(3):319-327.

［14］呂愛鐘，張路青.地下隧洞力學(xué)分析的復(fù)變函數(shù)方法［M］.北京:科學(xué)出版社，2007.

［15］呂愛鐘，蔣斌松，尤春安.位移反分析有限元網(wǎng)格劃分范圍的研究［J］.土木工程學(xué)報，1999，32(1):26-30.