18双孔圆形洞室围岩应力分析的交替法研究

岩石力学与工程学报

ChineseJournalofRockMechanicsandEngineering

17(5):534～543

Oct.,1998

双孔圆形洞室围岩应力分析的交替法研究*

张路青　吕爱钟

(山东矿业学院土木系　泰安　271019)

摘要　利用Schwarz交替法对双孔圆形洞室的围岩应力分析进行了研究。结果表明,完成二次迭代后,就可获得精度很高的满意结果。

关键词　Schwarz交替法,双孔圆形洞室,多余面力,二次迭代分类号　TU451

1　引言

地下工程中经常遇到两个平行洞室的开挖问题,当两洞室间距很大时,可以认为两洞室的开挖没有相互影响;当两孔洞间距较近时,孔洞的围岩应力状态就要受到较强的相互影响。两洞室的应力分析,在力学上属双连通域问题,国内外学者一直设法寻找一个行之有效的求解办法。有限单元法是求解该问题的一个常用方法,但有限单元法是一个数值方法,并且单元网格的划分对它的精度影响很大。Schwarz交替法是一种解析方法,它将双连通域问题归结为一系列的单连通域问题来求解,并且在数学上已证明了Schwarz交替法的收敛性。所以,只要该法的收敛速度快,它就可以作为解决双连通域问题的有效工具。

从60年代起,国外学者就开始使用Schwarz交替法解决一些简单的双连通域问题,文[1～4]分别对两圆洞室、两椭圆洞室问题进行了计算,但都只完成了一次迭代。迭代一次能否满足工程所需的精度,继续迭代后的精度又可提高多少?这些问题一直没有作出深入的研究。本文以双圆洞室为例得出二次迭代的结果,发现即使当两孔洞相距很近时,也可以获得满意的精度,这将为利用Schwarz交替法求解两任意形状洞室的应力分析打下基础。

2　交替法求解围岩应力的基本过程

Schwarz交替法的求解过程是这样的:首先在岩石中挖出第1个孔(图1),则可利用柯西积分法获得应力解。此时第2个孔虽未开挖,但其周边的应力可以求出。再挖第2个孔,即在此孔周边作用相应的平衡外力(以后称反面力)使孔边外载为零。如果设想此时第

1997年7月11日收到来稿。

*山东省自然科学基金资助项目(Q96A05110)。

作者张路青简介:男,25岁,硕士,1995年毕业于山东矿业学院土木系矿建专业,现于中国科学院地质研究所攻读博士学位,主要从事水文地质与工程地质方面的工作。[5]

第17卷　第5期张路青等.双孔圆形洞室围岩应力分析的交替法研究

・535・

1个孔未挖,则也是一个单连通域问题,可以求出在第2个孔周边的反面力作用下第1个孔周边的应力。如果这样算出的第1个孔周边面力为零分布,这两次计算结果叠加就相当于两孔都存在时的精确解答;如果不为零,则要在第1个孔周边加上一组平衡面力,使第1个孔周边的合面力为零。第一组平衡面力又引起第2个孔(设想第2个孔未挖)周边处的应力,这样反复做下来,直到两孔边的面力均为零,把每次计算结果叠加起来,就得原问题的解。实际上,由于迭代次数有限,只能使一个孔边的面力为零,使两个孔边的面力均为零是不可能的。由于每次每个孔上所施加的都是平衡面力,所以随着迭代的进行,第1次所要施加的平衡面力(以后称不等于零的面力为多余面力)将越来越小,当小到一定程度时,就可以得到精度很高的解析解。

本文定义每加1次反面力求解单连通域问题为完成一次迭代。

3　围岩应力分析的基本公式

当地下洞室埋置深度与孔径相比很大时,可以不考虑重力梯度的影响,把重力作用化为无限远处作用有P1,P2的外载来求解。为使求解的迭代结果具有一般性,决定两孔相对位置的C矢量可随意指定,如图1所示。图1中:　　P1——水平原岩地应力;

P2——垂直原岩地应力;Z1——X1O1Y1坐标系下坐标;Z2——X2O2Y2坐标系下坐标;

C——孔2圆心与孔1圆心相对位置。

由于利用复变函数解法,所以Z1,Z2,C都用复数表示。

在迭代过程中,需要进行坐标系变换。本文算例属坐标系平移Z1=Z2+C,坐标系平移前后两复应力函数之间变换关系[5]为

(1)　图1　双圆洞室计算示意图

′′″󰀁2(Z2)=󰀁1(Z2+C)+C󰀁1(Z2+C)(2)　Fig.1　Calculatingdiagramfortwo

由于常数不影响应力结果,分别对式(1),(2)左　circularholes右积分并舍去积分常数,可以得到变换关系为

󰀁2(Z2)=󰀁1(Z2+C)

′

󰀁2(Z2)=󰀁1(Z2+C)+C󰀁1(Z2+C)

′

󰀁′2(Z2)=󰀁1(Z2+C)

(3)(4)(5)(6)

类似地,可以得出从X2O2Y2坐标系平移至X1O1Y1坐标系时两复应力函数的变换关系为

󰀁1(Z1)=󰀁2(Z1-C)

󰀁1(Z1)=󰀁2(Z1-C)-C󰀁2(Z1-C)

标系下的两复应力函数。

有了上述关系后,就可以利用Schwarz交替法对双孔圆形洞室问题进行迭代求解。′

式中:󰀁1(Z1),󰀁1(Z1)为X1O1Y1坐标系下的两复应力函数;󰀁2(Z2),󰀁2(Z2)为X2O2Y2坐

・536・

岩石力学与工程学报1998年

第一次迭代过程为:首先求出P1,P2作用下只存在孔1的解,此时两解析函数为󰀁11

(Z1),󰀁11(Z1);󰀁11(Z1)与󰀁11(Z1)在X2O2Y2坐标系下为󰀁12(Z2),󰀁12(Z2),由于孔1的开挖而在孔2周边产生的多余面力为f12(t2)=󰀁12(t2)+t2󰀁′12(t2)+󰀁12(t2);为满足孔2的孔边应力边界条件必须加上反面力-f12(t2),进而解出在-f12(t2)作用下只存在孔2的解󰀁22(Z2),󰀁22(Z2),利用坐标变换Z2=Z1-C得到󰀁22(Z2),󰀁22(Z2)在X1O1Y1坐标系下的结果分别为󰀁21(Z1),󰀁21(Z1)。至此完成完整的一次迭代,得到如下结果:

AZ1Br1　　　󰀁11(Z1)=-(7)

2Z1

Ar21Br41

󰀁11(Z1)=BZ1-Z1-3(8)Z1r223r222()3()2

Br1Br22CC󰀁21(Z1)=--+222

Z1-Cr2(C)(Z1-C+r2)2

Z1-C+

CC

Ar21Br41BCr21

+-(9)222

r2r2r232

(Z1-C+C)(Z1-C+C)Z1-C+C

2

Br22Br21Br21r2

　　　󰀁21(Z1)=-+--Z1-Cr22(Z1-C)C2Z1-C+C

r22

󰀁′21(Z1)-C󰀁′21(Z1)(10)

Z1-CP1+P2P1-P2

式中:A=-,B=(本文应力以压为正,以拉为负)。第一次迭代后总的22

应力结果为

󰀁1(Z1)=󰀁11(Z1)+󰀁21(Z1)󰀁1(Z1)=󰀁11(Z1)+󰀁21(Z1)

(11)(12)

2

　　第二次迭代在第一次迭代的基础上继续进行,由󰀁1(Z1),󰀁21(Z1)在孔1周边产生的多余面力为f21(t1)=󰀁21(t1)+t1󰀁′21(t1)+󰀁21(t1)。为使孔1周边满足应力边界条件,在孔1周边作用反面力-f21(t1),进而求出反面力-f21(t1)作用下只有孔1的解󰀁31(Z1),󰀁31(Z1),则第二次迭代后总的应力结果为

󰀁1(Z1)=󰀁11(Z1)+󰀁21(Z1)+󰀁31(Z1)󰀁1(Z1)=󰀁11(Z1)+󰀁21(Z1)+󰀁31(Z1)

21(Z1)是经过复杂的推导而来的,本文只列出结果:式中的󰀁31(Z1),󰀁

(13)(14)

　　󰀁31(Z1)=

r21

-CZ1

2Ar2

2

-

Br1C2r21r222(-C+)Z1Cr21

-CZ1r223()2Cr22

Br1

2

Br1r2r2

++(2+C)󰀁2

r12r1

+CC(-C)-C

Z1Z1

r222()3C2222

-r21r223(-C+)Z1C+r212(Z1-C)Br2+第17卷　第5期

2

Ar21r2

2

张路青等.双孔圆形洞室围岩应力分析的交替法研究

2

3Br41r2

4

2

2BCr21r2

3

・537・

(

r1

-CZ1

2r22

r212

+C)(Z1-C)

+(

2r22

r1

-CZ1

r212

+C)(Z1-C)

-(

2r22

r1

-CZ1

r212

+C)(Z1-C)-

3Br1r222C(CC-r22)

24

r21Cr21C32

(3(2)2)26CC-r2CC-r22Br1r2

+-223󰀁r21Cr21C2C(CC-r2)(Z1-(Z1-2)2)

CC-r2CC-r2

r21Cr21Cr21C43

(4(6()2)2)CC-r2CC-r2CC-r22

2+2+2-r1Cr1Cr1C32

(Z1-)(Z1-)(Z1-)

CC-r22CC-r22CC-r22Br22

2C

22

r21Cr21Cr13r1232(3(()3()222)2)Ar1r2CC-r2CC-r2CC2+2-2+2+22r12r1r1Cr1C(CC-r)2

(Z1-)Z1-(Z1-)Z1-CCCC-r22CC-r2242

12

22

22223

Z31(r1-CZ1)CZ41-r1Z122

CZ1--CZ13Brr2BCr1r2

r21Cr21C4+3423

(Z1-(CC-r)2)(CC-r2)(Z1-2)

CC-r2CC-r2

22r22r23

23()()22Br1CCBr2Ar1

󰀁31(Z1)=22-22+2--C2r1r2r1r22r1r22

(-C+)(-C+)-C+CZ1CZ1CZ1r21

-CZ1

Br41BCr21r21r21′′2+2-Z1󰀁21(0)-Z1󰀁31(Z1)r2r2

(2+C)3(2+C)2r1r1

-C-CZ1Z1

(15)

(16)

第二次迭代完成之后,孔1周边的应力边界条件(零面力)精确满足,孔2周边存在多余面力。有了式(11)～(14),可根据式(17),(18)求出围岩内任一点的应力分量。

󰀁r+󰀁󰀁=4Re[󰀁1(Z1)]

″′2i󰀁

󰀁󰀁-󰀁r+2i󰀁r󰀁=2[Z1󰀁1(Z1)+󰀁1(Z1)]e

′

(17)(18)

式中:󰀁r,󰀁󰀁,󰀁r󰀁分别为径向应力,环向应力,剪切应力;󰀁的含义见图1。

4　收敛精度的讨论

为方便精度分析起见,取r1=r2=1.5m,并使两孔洞圆心在同一水平线上。又因为孔洞近距离情况下是我们关心的内容,取孔距为4.0m,此时两孔内侧两点距离为1.0m。本文

从以下3个方面对一次迭代和二次迭代的应力计算结果进行精度分析。4.1　由孔边多余面力来衡量迭代结果的精度

孔边多余面力包括孔边多余的󰀁r和󰀁r󰀁。由前面分析可知,孔边存在的多余面力使迭代结果与精确解不一致。要使孔边的多余面力绝对为零面力是不可能的,也是没有必要的。・538・

岩石力学与工程学报1998年

由以下计算结果可知:当孔边多余面力足够小时,计算得到的围岩应力就接近精确解。

表1,2给出了在P1=P2=20MPa时孔边不同点的应力值,表1给出的是一次迭代在孔1周边有多余面力的结果,表2给出的是二次迭代在孔2周边有多余面力的结果。由于对称性,表中只给出了两孔上半部分孔边各点的计算结果,并由于孔边位置分别以X1,X2轴正向逆时针算起,所以把表1中的孔边位置0°～180°分别与表2的孔边位置180°～0°相对应。

表1　一次迭代在孔1周边各点的应力状态

　　　　　Table1　Thestresscomponentsontheboundaryofhole1forthefirstiteration　　　MPa

孔边位置/(°)

应力

0

󰀁r󰀁󰀁

󰀁r󰀁

-1.68　48.30　0.00　

300.76　43.09　4.42　

603.46　38.31　1.82　

902.49　38.73　-0.78　

1200.91　40.18　-1.34　

150-0.12　41.16　-0.86　

180-0.45　41.49　0.00　

表2　二次迭代在孔2周边各点的应力状态

　　　　　Table2　Thestresscomponentsontheboundaryofhole2fortheseconditeration　　　MPa

孔边位置/(°)

应力

180

󰀁r󰀁󰀁󰀁r󰀁

2.95　70.27　0.00　

1500.66　41.62　0.48　

1200.22　26.98　-0.53　

900.55　32.92　-0.47　

600.　40.22　-0.25　

300.65　44.63　-0.10　

00.　46.05　0.00　

由于本算例具有几何对称性和荷载对称性,故可将第一次迭代与第二次迭代所得到的

孔边应力直接进行精度比较。本文根据孔边各点面力绝对值的均值来判定迭代的收敛效果;多余面力绝对值的均值反映了与实际零面力的差距,这个值越小说明迭代结果越接近

r󰀁精确解。把多余面力绝对值的均值分别记为󰀁r和󰀁。由表1,表2可知,第一次迭代:󰀁r1=

1.41MPa,󰀁r󰀁1=1.32MPa;第二次迭代:󰀁r2=0.90MPa,󰀁r󰀁2=0.26MPa,由此可得:󰀁r2/󰀁r1=0.,󰀁r󰀁2/󰀁r󰀁1=0.20。可见第二次迭代要比第一次迭代更接近实际的孔边应力边界条件

(󰀁r=󰀁r󰀁=0)。由表2可以看出,对第二次迭代,多余面力与外荷载P1=P2=20MPa相比已经很小,这间接地说明了二次迭代计算结果已足够精确。4.2　多余面力对围岩应力状态的影响程度

如果多余面力的变化对围岩应力状态的影响不大,这时才可以有充分的理由说二次迭代的结果已接近于精确解。由于孔边的󰀁󰀁是人们最关心的,在迭代过程中变化也是最大的,所以本文只考虑多余面力对孔边󰀁󰀁的影响程度。

一次迭代和二次迭代孔边多余面力绝对值的均值分别记为󰀁r1,󰀁r󰀁1和󰀁r2,󰀁r󰀁2,则两次迭代的变化量分别为󰀁󰀁r2-󰀁r1󰀁和󰀁󰀁r󰀁2-󰀁r󰀁1󰀁。多余面力的变化必定会引起孔边切向应力󰀁󰀁的变化。若两次迭代孔边切向应力的均值分别记为󰀁󰀁1和󰀁󰀁2,则下面两个指标:

第17卷　第5期张路青等.双孔圆形洞室围岩应力分析的交替法研究

・539・

F󰀁r=

󰀁󰀁󰀁2-󰀁󰀁1󰀁󰀁󰀁󰀁2-󰀁󰀁1󰀁,　F󰀁r󰀁=

󰀁󰀁r2-󰀁r1󰀁󰀁󰀁r󰀁2-󰀁r󰀁1󰀁

可以反映多余面力的变化对孔边切向应力󰀁󰀁变化的影响程度。前者表示孔边󰀁r单位变化量引起孔边󰀁󰀁的变化量;后者表示孔边󰀁r󰀁单位变化量引起孔边󰀁󰀁的变化量。对于本算例,F󰀁r=3.20,F󰀁r󰀁=1.54。由此可以计算出二次迭代后2个多余面力所引起的󰀁󰀁变化量分别为2.88MPa,0.40MPa,与最大切向应力(70.27MPa)相比已足够精确。

既然多余面力的变化对变化显著的孔边󰀁󰀁影响很小,那么它对稍远处的围岩应力影响将会更小。由前面分析得出二次迭代结果的多余面力已很接近零面力的结论,此时,我们有充分的理由说:二次迭代已有很高的精度,围岩应力接近精确解。4.3　利用对称性来检验迭代结果的精度

如图2所示,两孔半径r1=r2=1.5m,孔心在同一水平线上,且孔心相距4.0m,无限远处作用荷载P1=P2=20MPa。由于该问题的几何边界,荷载条件都关于Ⅰ-Ⅰ轴对称,所以两孔最终的应力状态也应关于Ⅰ-Ⅰ轴对称,但由于多余面力的存在会使所得结果并非关于Ⅰ-Ⅰ轴对称。孔边多余面力越接近零面力,Ⅰ-Ⅰ轴两侧的非对称性就会越来越弱,而对称性就会越来越强。于是,可以利用对称性的好坏来检验迭代结果的精度。

由于应力集中通常发生在孔边的应力分量󰀁󰀁,离孔越远,孔的存在对应力状态的影响就会越小。也就是说,多余面力对孔边的󰀁󰀁影响最大,孔边󰀁󰀁的对称性也最差,因此本文把对称性最差的孔边󰀁󰀁作为考察对象。如图2所示,孔边相同编号的点是关于Ⅰ-Ⅰ轴对称的点,表3,4分别给出了一次迭代和二次迭代时对称点的应力分量󰀁󰀁。

从表3和表4所列出的󰀁󰀁󰀁1-󰀁󰀁2󰀁大小可以看出,二次迭代结果要比一次迭代结果在孔边的对称性提高数倍,对称性的提高直接说明了二次迭代精度的提高。本文对其他几何边界、荷载条件下的情况进行了计算,都得到了二次迭代精度大大高于一次迭代的结论。

表3　一次迭代时对称点的应力分量󰀁󰀁

　　　　　　Table3　Thetangentstress󰀁󰀁ofsymmetricalpointsforthefirstiteration　　　　　MPa

编号

孔号

1

孔1孔2󰀁󰀁󰀁1-󰀁󰀁2󰀁

48.39　68.80　20.41　

239.82　28.23　11.59　

338.73　32.57　6.16　

440.75　42.70　1.95　

541.49　45.95　4.46　

0.75　42.70　1.95　

738.73　37.57　6.16　

839.82　28.23　11.59　

表4　二次迭代时对称点的应力分量󰀁󰀁

　　　　　　Table4　Thetangentstress󰀁󰀁ofsymmetricalpointsfortheseconditeration　　　　MPa

编号

孔号

1

孔1孔2󰀁󰀁󰀁1-󰀁󰀁2󰀁

62.23　70.27　8.04　

236.41　29.　6.52　

334.62　32.92　1.70　

443.02　42.82　0.20　

546.03　46.05　0.02　

3.02　42.82　0.20　

734.62　32.92　1.70　

836.41　29.　6.52　

・540・

岩石力学与工程学报1998年

图2　孔和外载关于Ⅰ-Ⅰ轴对称时的孔边对称点(45°间隔)

Fig.2　Symmetricalpoints(interval45°)ontheboundaryoftwoholeswiththesymmetricalⅠaxisofholesandloadⅠ-

图3　P1=10MPa,P2=20MPa,r1=r2=1.5m,C=(5.0,0.0)情形下二次迭代围岩应力状态Fig.3　ThestressconditionofsurroundingrockwithP1=10MPa,P2=20MPa,r1=r2=1.5m,C=(5.0,0.0)

以上结论是在两孔孔心位置矢量C=(4.0,0.0)的情形下得到的,随两孔孔心距离的加大,精度将会更高。

5　不同条件下的围岩应力状态

利用二次迭代的两复应力函数,可以计算得到在各种几何条件、荷载条件下两圆形洞室的围岩应力。从计算结果看出:由于孔的存在使孔边部分点有不同程度的应力集中,随着离孔距离的增加,围岩应力状态很快趋近于原岩应力;对几何对称、荷载对称下的两孔洞室,计算得到的围岩应力状态满足关于Ⅰ-Ⅰ轴对称的精度要求(图3)。

随着荷载、孔径、孔心相对位置的改变,孔边应力集中的位置和大小也要发生改变。表5,6列出的是两孔洞室在不同荷载下孔边对应点󰀁󰀁的变化。由表中数据可以发现,随侧

表5　r1=r2=1.5m,C=(5.0,0.0)时两孔洞室孔边应力分量󰀁󰀁

Table5　Thetangentstress󰀁󰀁ontheboundaryoftwoholeswith

　　　　　　　　　　　　　　r1=r2=1.5m,C=(5.0,0.0)　　　　　　　　　　　　　　MPa

荷载

对应点孔10°

90°180°270°孔2180°90°0°270°P1=5P2=20孔160.15-4.3857.74-4.38孔259.23-4.0957.51-4.09P1=10P2=20孔157.618.9853.218.98孔257.708.8753.108.87P1=20P2=20孔152.5335.7044.1735.70孔254.9034.7844.2934.78P1=20P2=10孔121.1944.5713.0444.57孔224.53.3013.3343.30P1=20P2=5孔15.5149.00-2.5349.00孔29.3947.56-2.1547.56第17卷　第5期张路青等.双孔圆形洞室围岩应力分析的交替法研究

・541・

表6　r1=1.5m,r2=2.0m,C=(5.0,0.0)时两孔洞室孔边应力分量󰀁󰀁

Table6　Thetangentstress󰀁󰀁ontheboundaryoftwoholeswith

　　　　　　　　　　　　　　r1=1.5,r2=2.0m,C=(5.0,0.0)　　　　　　　　　　　MPa

荷载

对应点孔10°90°180°270°

孔2180°90°0°270°

P1=5P2=20孔168.78-4.1260.26-4.12

孔261.73-4.9757.56-4.97

P1=10P2=20孔166.578.0855.998.08

孔261.398.6152.958.61

P1=20P2=20孔162.1532.4847.4632.48

孔260.7135.73.7235.76

P1=20P2=10孔126.6540.15.2040.

孔229.6845.0312.6345.03

P1=20P2=5孔18.9044.72-0.9344.72

孔214.19.66-2.9149.66

图4　P1=5MPa,P2=20MPa,C=(5.0,0.0)时两孔间拉应力区Fig.4　ThetensilezonewithP1=5MPa,P2=20MPa,C=(5.0,0.0)

图5　P1=20MPa,P2=5MPa,C=(5.0,0.0)时两孔间拉应力区Fig.5　ThetensilezonewithP1=20MPa,P2=5MPa,C=(5.0,0.0)

图6　P1=5MPa,P2=20MPa,C=(3.0,4.0)时两孔间拉应力区Fig.6　ThetensilezonewithP1=5MPa,P2=20MPa,C=(3.0,4.0)

图7　P1=20MPa,P2=5MPa,C=(3.0,4.0)时两孔间拉应力区Fig.7　ThetensilezonewithP1=20MPa,P2=5MPa,C=(3.0,4.0)

・542・

岩石力学与工程学报1998年

应力系数(P1/P2)的加大,最大应力集中的位置由两孔边内侧两点转向上下顶底部位;在等径孔情况下,两孔边对应点󰀁󰀁的微小差别是孔2周边的多余面力所引起(表5),在不等径孔情况下,两孔边对应点󰀁󰀁的不同主要是由于孔径的不同所引起的(表6)。

从表5,6还可发现,在侧应力系数很小的情况下,两孔的上下顶底4个区域出现拉应力区,最大拉应力在上下4个孔边顶点上;在侧应力系数很大的情况下,两孔外侧两个区域出现拉应力区,最大拉应力发生在外侧2个孔边顶点上。对于孔心不在同一水平线或竖直线上时,拉应力出现的位置和最大拉应力值有所不同。图4～7表示出了r1=r2=1.5m的两孔洞室在不同荷载和孔心位置时拉应力区出现的位置和最大拉应力值;因最大拉应力值只发生在孔边,所以图中表示󰀁󰀁最大值的位置以X1,X2轴正向逆时针转角来计算。

6　结论

(1)当两孔洞室相距很近时,由一次迭代求出的围岩应力不能满足精度要求,其结果与精确解相差甚远。

(2)二次迭代结果的精度远远高于一次迭代结果的精度,即使在两孔相距很近的情况下,二次迭代仍然有很高的精度。

(3)由于岩石的抗拉强度小,为了维护围岩,必须知道拉应力出现的位置和大小,利用本文所给出的二次迭代结果可以准确地找到各种条件下拉应力出现的位置和大小,为工程的开挖和支护提供依据。

(4)两孔洞室的空间位置不同,所产生的拉应力区和最大拉应力值也不同。由图4～7可见,当孔心布置不在同一水平线或竖直线上时(图6,7),孔边拉应力区范围和最大拉应力值要比同一水平线或竖直线上两孔洞室(图4,5)大得多。

(5)本文利用Schwarz交替法对双圆洞室围岩应力分析进行了研究,为求解两个任意形状洞室的应力分析奠定了基础。

参

考

文

献

1　SalernoVL,MahoneyJB.Stresssolutionforaninfinteplatecontainingtwoarbirarycircularholesunderequalbiax-ialstresses.JournalofEngineeringforIndustry,Trans,ASME,SeriesB,1968,90:419～426

2　UkadgaonkerVG.Stressanalysisofaplatecontainingtwocircularholeshavingtangentialstresses.AIAAJournal,

1980,18:125～128

3　KimTJ,UkadgaonkerVG.Planestressanalysisoftworigidcircularinclusions.AIAAJournal,1971,9:2294～

2296

4　UkadgaonkerVG.Interactioneffectoftwoarbitranilyorientedcracks—partI.Int.JournalofFracture,1991,51:

219～230

5　陈子荫.围岩力学分析中的解析方法.北京:煤炭工业出版社,1994

第17卷　第5期张路青等.双孔圆形洞室围岩应力分析的交替法研究

・543・

STUDYOFALTERNATIONMETHODFORSTRESSANALYSISONSURROUNDINGROCKSOFTWOCIRCULARHOLES

ZhangLuqing　LuAizhong

(ShandongInstituteofMiningandTechnology,　Tai′an　271019)

Abstract　StressanalysisismadeforthesurroundingrocksoftwocircularholesusingSchwarz′salternatingmethods.Bycalculating,itisfoundthatthesatisfyingresultswith

highprecisioncanbeacquiredaftertheseconditeration.

Keywords　Schwarz′salternatingmethods,twocircularholes,theresidualsurfaceforce,theseconditeration