《一种任意开挖断面隧洞的围岩应力分析方法.pdf》由会员分享，可在线阅读，更多相关《一种任意开挖断面隧洞的围岩应力分析方法.pdf（18页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102982247 A(43)申请公布日 2013.03.20CN102982247A*CN102982247A*(21)申请号 201210525284.7(22)申请日 2012.12.05G06F 19/00(2006.01)(71)申请人浙江科技学院地址 310023 浙江省杭州市西湖区留和路318号(72)发明人祝江鸿 李明友 王渊 朱鹏(74)专利代理机构浙江杭州金通专利事务所有限公司 33100代理人王桂名(54) 发明名称一种任意开挖断面隧洞的围岩应力分析方法(57) 摘要本发明涉及一种隧洞围岩稳定性分析和评价领域里围岩应力分布的确定方法，包括下述步骤。
2、：（1）建立力学模型：设围岩容重为、侧压力系数为、隧道埋深为H及开挖跨度为D当隧道埋深H与开挖跨度D满足H10D时,隧道原始地应力场简化为无穷远处均布的竖向应力、水平应力，剪应力；（2）建立任意开挖断面与变量的映射函数；（3）以及表示的解析函数和(z)；（4）按下列复变函数计算而得隧洞开挖后的围岩应力；（5）对计算后的围岩应力进行分析。本发明填补了地下隧洞围岩应力解析解分析的空白，具有较高的精度和较高的工程可靠度。(51)Int.Cl.权利要求书2页 说明书8页 附图7页(19)中华人民共和国国家知识产权局(12)发明专利申请权利要求书 2 页 说明书 8 页 附图 7 页1/2页21.一种任。
3、意开挖断面隧洞的围岩应力分析方法，其特征在于包括下述步骤：（1）建立力学模型：设围岩容重为、侧压力系数为、埋深为H和开挖跨度为D，当隧道埋深H与开挖跨度D满足H10D时,隧道原始地应力场简化为无穷远处均布的竖向应力、水平应力，剪应力；（2）建立任意开挖断面与变量的映射函数：对隧洞的任意开挖断面，将任意开挖断面看成由平面上的单位圆映射而得，从而将单位圆外域映射到任意断面隧洞外域，式中：z 为 z 平面上单连通外域内的点坐标；为平面上单位圆外域内的点坐标；为映射函数，m为映射函数项数，其取值与实际开挖洞形的复杂程度及映射要求的精度有关，m1；cn为与隧道开挖断面形状相关的复常数；（3）采用及表示解。
4、析函数和(z)：其中， ，为c1的共轭；（4）隧洞开挖后的围岩应力按下列复变函数计算而得：式中x和y分别为围岩任意点的x和y方向正应力，xy为剪应力；Re和Im分别表示取复数的实部和虚部；为z的共轭；（5）对计算后的围岩应力进行分析。2.根据权利要求1所述的任意开挖断面隧洞的围岩应力分析方法，其特征在于：所述权 利 要 求 书CN 102982247 A2/2页3的开挖断面形状为圆形，式(1)中m=1, z=()=R,R为实际隧洞开挖断面半径。3.根据权利要求1所述的任意开挖断面隧洞的围岩应力分析方法，其特征在于：式(1)中cn=An+iBn，An、Bn按三角插值法求解而得。权 利 要 求 书。
5、CN 102982247 A1/8页4一种任意开挖断面隧洞的围岩应力分析方法技术领域0001 本发明涉及一种隧洞围岩稳定性分析和评价领域里围岩应力分布的确定方法，它用来在隧洞开挖后的围岩稳定性分析与评价过程中快捷方便高效地确定围岩应力的分布和开展围岩应力分析。背景技术0002 通常的，地下隧道围岩应力分析的方法有以下四种。一种是数值计算法，需要对地下隧洞围岩进行弹塑性力学分析，确定具体边界条件下的围岩位移场或应力场，然后根据求得的结果进行围岩应力分析和稳定性评价，以确定地下隧洞结构设计的可靠性和施工方案的合理性。具体计算确定方法有有限元法、边界元法和离散元法等。这个方法的缺点是：针对相类似的工。
6、程分析结果往往差异较大，且计算结果与实际情况之间也存在着一定的差距，这与计算时材料模型选取、计算参数取值、计算模式确定及对问题的简化程度等方面是有关的，因此对于数值模拟当前采取的是“定量分析，定性应用”一种理念；其次是效率低，特别是当边界条件发生变化时（荷载及材料属性等），需要针对每一种边界条件重新建模计算应力场或位移场。对于工程问题，每次有限元法或边界元法计算的工作量都较大，所以总效率很低；第二种方法是模型试验，该方法能把工程结构与围岩作为统一体考虑，较好地反映了岩体特性且能模拟复杂工程结构与地质环境。只要能满足一定的相似关系，不必建立复杂的本构关系或进行严密的计算分析，可直接通过测试得出结。
7、果，省去了数学、力学上的麻烦。但模型试验也存在尺寸效应、试验难度大以及费用高等不足。第三种方法是现场监控量测，通过该方法可动态地掌握围岩应力、变形及支护结构的工作状态。通过该方法获取的具体位移和应力信息就工程本身实践来讲意义是重大的，但无法对同类问题做规律性研究。第四种方法是解析法，该方法采用数学力学的计算取得围岩应力和位移的闭合解。利用解析方法讨论隧道围岩稳定性的优势在于所获得的是精确解，对计算所涉及到的各参数讨论比较方便并容易得到规律性的认识。然而，受制于数学计算方法研究的缺失，单洞隧道围岩应力场和位移场解析解的研究成果更多地集中在圆形隧道方面，而对于任意开挖断面的隧洞围岩应力分析成果尚未。
8、成熟。0003 现有技术存在的缺点是：(1)运用数值计算方法进行任意开挖断面隧洞围岩应力分析时，针对相类似的工程分析结果往往差异较大，且计算结果与实际情况之间也存在着一定的差距，对实际工程分析时的效率也较低；当前对该方法的使用是基于一种“定量分析，定性应用”的理念。（2）模型试验存在尺寸效应、试验难度大以及费用高等不足。(3) 通过现场监控量测仅能获取具体工程的位移和应力信息，由于实际工程的边界条件千差万别，该方法对同类问题无法做规律性研究。（4）受制于数学计算方法研究的缺失，任意开挖断面隧洞的围岩应力解析解的研究成果主要集中在圆形和椭圆形隧洞方面，对其它任意开挖断面的研究成果不成熟，且未给出。
9、解析解显式，通用性差，利用现有成果无法对任意开挖断面的隧洞围岩应力进行弹塑性解析分析。说 明 书CN 102982247 A2/8页5发明内容0004 为克服已有的地下隧道任意开挖断面围岩应力分析方法的不足，弥补由于数学力学方法研究的缺失而造成的任意开挖断面隧洞围岩应力解析分析的不成熟，填补任意开挖断面围岩应力尚未有解析解通式的空白，本发明提供一种独特的地下隧洞围岩应力解析分析方法，能高精度和高效率地解决任意开挖断面围岩应力的求解问题，快捷地进行围岩应力分析，使进行任意开挖断面隧洞的围岩应力弹塑性解析分析成为可能。0005 本发明解决其技术问题所采用的技术方案是：0006 一种任意开挖断面隧洞。
10、的围岩应力分析方法，包括下述步骤：0007 （1）建立力学模型：0008 设围岩容重为、侧压力系数为、隧道埋深为H及开挖跨度为D，当隧道埋深H与开挖跨度D满足H10D时,隧道原始地应力场简化为无穷远处均布的竖向应力、水平应力，剪应力；0009 （2）建立任意开挖断面与变量的映射函数：0010 对隧洞的任意开挖断面，将任意开挖断面看成由平面上的单位圆映射而得，从而将单位圆外域映射到任意断面隧洞外域，0011 0012 式中：z 为 z 平面上单连通外域内的点坐标； 为平面上单位圆外域内的点坐标；为映射函数，m为映射函数项数，其取值与实际开挖洞形的复杂程度及映射要求的精度有关，m1，cn是与开挖断。
11、面形状相关的复常数；0013 （3）以及表示的解析函数和(z)：0014 0015 0016 其中， ，为c1的共轭；0017 （4）按下列复变函数计算而得隧洞开挖后的围岩应力：0018 0019 说 明 书CN 102982247 A3/8页60020 0021 式中x和y分别为围岩任意点的x和y方向正应力，xy为剪应力；Re和Im分别表示取复数的实部和虚部；为z的共轭；0022 （5）对计算后的围岩应力进行分析。0023 进一步，所述的开挖断面形状为圆形，式(1)中m=1,z=()=R,R为实际隧洞开挖断面半径。0024 式(1)中，映射函数 可采用当前本领域技术人员所熟知的数字方法求解，。
12、如朱大勇、钱七虎等在岩石力学与工程学报 第18卷第3期所发表的“复杂形状洞室映射函数的新解法” ，或者皇甫鹏鹏、伍法权等在岩土力学第32卷第5期所发表的“基于边界点搜索的洞室外域映射函数求解法”，或者按通用的三角插值法求解。0025 进一步，式(1)中cn=An+iBn，An、Bn按三角插值法求解而得。0026 本发明步骤（1）力学模型的建立基于以下理论：对地下隧洞围岩的弹性应力分析，若满足隧洞围岩均质、各向同性和线弹性的假设条件，隧道长度远大于开挖跨度，则可采用平面应变的方法，取隧洞的任一开挖断面作为其代表研究。研究表明，当隧洞符合埋深条件，即隧道埋深H与开挖跨度D满足H10D时，忽略隧洞影。
13、响范围内岩石（1.5-2.5倍的D）的自重，原问题的误差不超过10%。于是，原岩应力仅考虑自重应力场时，隧洞原岩应力可简化为无穷远处均布的竖向应力和水平应力，对隧道围岩的弹性应力分析可按平面弹性洞口问题考虑。 0027 本发明的技术方案要点是：0028 首先是建立其进行地下任意开挖断面隧道围岩应力分析的力学模型；其次是建立起隧洞围岩应力及边界条件方程在实数域与复数域之间的联系方程；第三是通过三角插值理论解决联系方程中的映射函数的求解问题；第四是将任意开挖断面隧洞围岩应力边界条件方程转化成积分方程，通过对积分方程中各被积函数在讨论域内的解析性研究获知每一被积函数的解析性；第五是利用留数理论直接求。
14、得边界条件方程中的两个解析函数；第六是将获得的解析函数代回隧洞围岩应力在实数域与复数域之间的联系方程中，得到实数域内的围岩应力求解式。0029 当具体地下隧洞的边界条件确定后，利用上述方法可以方便快捷的进行任意开挖断面隧洞的围岩应力计算和分析，具体分析时仅需将具体边界条件输入便可获得围岩应力计算结果。本方法的实际执行效率，可以比现有的其它方法提高几十倍乃至上百倍；特别是对于边界条件有微小变化时，利用现有的其它方法，每一微小变化就需建模和计算一次，计算周期较长，且所获得的计算结果由于精度问题往往不理想。利用本方法可以在数分钟至几小时内完成具体的工程分析，且分析结果符合工程精度要求。0030 相对。
15、于现有技术，本发明的有益效果主要表现在：0031 (1)、填补了地下隧洞围岩应力解析解分析的空白。现有的地下隧洞围岩应力的解说 明 书CN 102982247 A4/8页7析法分析，仅在圆形和椭圆形隧洞方面较成熟；对于其它隧洞开挖洞形的围岩应力分析，现有的解析方法计算复杂、通用性差，且无法利用它进行后续的围岩稳定性弹塑性分析。本方法不仅能快速计算获得任意开挖断面隧洞的围岩应力并展开围岩应力分析，且给出了弹性解的显式，通用性强，运用该成果使围岩稳定性的弹塑性解析解分析成为可能，填补了该领域的空白。0032 (2)、具有较高的精度。只要映射函数是符合精度要求的，将其代入两个解析函数通式进行隧洞围岩。
16、应力求解时，无需对映射函数项数作截项处理，且所获得的应力较其它解析近似解更为精确。0033 （3）较高的工程可靠度。由于所获得的围岩应力是解析解，即在具体隧洞围岩应力边界积分方程求解过程中，被积函数的奇点个数和类型在讨论域内一旦明确，则可应用留数理论直接求得每一积分项的积分，无需采用新的级数来代替被积函数，因此其在数学和力学上属于闭合解，具有较高的精度；边界条件的任何微小调整均体现在围岩应力值的变化上，在此基础上获得的围岩应力分析结果较其他方法更为可靠。0034 （4）具有较强的适用性。实际地下隧洞工程中，由于工程功能性需要，实际开挖断面形状是多种多样的。利用本方法可对任意开挖断面隧洞的围岩应。
17、力展开分析，而无需借助其他手段，具有较强的适用性。附图说明0035 图1是本发明地下隧洞任意开挖断面的围岩应力分析方法的流程图。0036 图2是任意开挖断面隧洞围岩应力的力学模型。0037 图3是任意开挖断面隧洞围岩应力的计算模型。0038 图4是单位圆任意断面外域映射关系图。0039 图5是隧道圆形开挖断面的围岩应力分析力学模型。0040 图6是圆形隧道环向应力分布图。0041 图7是圆形隧道=0 ,45和90三个方向上围岩应力变化规律图。0042 图8是高速公路两车道隧道断面图。0043 图9是高速公路两车道隧道孔边环向应力分布图。0044 图10是高速公路两车道隧道围岩应力分析图。004。
18、5 图11是高速公路两车道隧道分析点围岩应力分析图。具体实施方式0046 下面结合附图对本发明作进一步描述，但本发明的保护范围并不限于此。0047 实施例10048 参照图1、2、3、4和5，利用一种地下隧道任意开挖断面的围岩应力分析方法，对圆形隧道围岩应力展开分析，该方法主要包括以下步骤：0049 （1）计算条件和参数。0050 有一R=3m的圆形隧道，力学模型见图5，计算模型参见图2。由图5知： H=630m20R0。围岩容重为=23KN/m3，侧压力系数为0.54。相关参数计算如下：说 明 书CN 102982247 A5/8页80051 0052 0053 （2）、映射函数求解：005。
19、4 由于研究对象为圆形隧洞，由图4知映射函数为：0055 z=()=R (1)0056 其中R为实际隧洞开挖断面半径，从图5的初始条件易知R为3m。0057 （3）解析函数获取：0058 将式(1)代入式（2）和（3），并转化成z变量表示：0059 0060 0061 （4）围岩应力求解。0062 将(2)、(3)式以及计算参数代入下式便可求得圆形隧洞围岩应力：0063 0064 0065 0066 （5）围岩应力分析。0067 洞周环向应力分布0068 图6为洞周环向应力分布图。由图可知，洞周只存在环向应力而r,r为零，圆形隧道洞周环向应力均为压应力，呈中心对称关系。应力大小分布与角度有关，。
20、最大压应力在=0及=180位置，达到35.68MPa，为边界竖向应力的2.46倍，最小压应力在在=90及=270的位置,数值为 8.96MPa，与边界水平应力较接近。由此看出，对于断面形状为圆形的隧道，洞周环向应力分布较均匀，有利于保证围岩稳定性。0069 图7分别为=0 ,45和90三个方向上围岩各分析点的x、y和xy应力分布图，其中分析点是按单位圆外上述三个方向上的点映射到洞室外的实际点考虑，因此，图7中的横坐标是按映射点到单位圆洞壁的距离表示，r为单位圆的半径。由图7（a）和说 明 书CN 102982247 A6/8页9（c）知，在=0 和90处围岩各点的xy均为0，即隧道的水平和铅直。
21、面为主应力面，在该面上只有正应力没有剪应力，其余方向截面都有剪应力，这是由计算模型和计算条件决定的。计算模型中仅考虑自重应力场，将埋深范围内的岩石自重和侧压力简化为无穷边界处的作用力；且不考虑重力梯度。因此计算结果和计算模型是相吻合的。从图7（a）可以看出，x从洞周开始随着离洞周的距离加大应力值不断增加，在2倍洞径范围内增加幅度最大，随后增加值减缓，当距离洞周3r以上时，x趋近于7.82MPa；y的变化刚好相反，随着距离洞周越远，应力值从洞周的最大值逐渐趋近于14.49MPa，且在一倍洞周范围内应力下降幅度最大。图7（b）和（c）也可以看出同样的规律，这说明隧道开挖引起了洞周附近围岩应力场发生。
22、较大变化，随着离洞周距离的增加，围岩受到的扰动越来越小。以上所反应的规律和隧道工程实践中所得规律较一致，说明了应用复变函数法求解隧道围岩应力场是可行和可信的。0070 实施例20071 参照图1、2、3、4和8，利用一种地下隧道任意开挖断面的围岩应力分析方法，对以两车道公路隧道标准轮廓线为基础的隧道进行开挖后的隧道围岩应力分析，该方法主要包括以下步骤：0072 （1）计算条件和参数。0073 计算对象为埋深120m的两车道公路隧道，围岩级别，重度为23kN/m3，侧压力系数为0.54。公路隧道设计规范(JTG D702004)中的两车道公路隧道标准内轮廓断面为图8（a）；隧道实际开挖断面一般距。
23、标准内轮廓线50-500mm左右，视围岩等级、埋深情况及支护设计等条件确定。本次计算考虑开挖宽度为距标准轮廓线300mm，具体断面尺寸见图8（b）。其它相关计算所需参数计算如下：0074 0075 0076 （2）建立任意开挖断面与变量的映射函数：0077 映射函数采用式（1），将单位圆外域映射到图8（b）的隧道实际开挖断面外域，函数项数m取32项，映射断面与实际开挖断面的边界线误差控制在5mm以内，系数cn=An+iBn，An、Bn采用三角插值法求解，表1为映射函数系数的计算结果。0078 表1高速公路两车道隧道开挖断面的映射函数系数说 明 书CN 102982247 A7/8页100079。
24、 0080 （3）解析函数获取。0081 将上述映射函数代入下式便可求得两解析函数。0082 0083 0084 其中；0085 （4）围岩应力求解。0086 参照实施例1中的围岩应力求解公式进行围岩应力计算。0087 （5）围岩应力分析。0088 孔边环向应力0089 开挖后洞周不考虑支护，因此隧道的孔边边界面力为零，即只存在环向应力，而径向应力和剪应力均为零。图9为隧道开挖后的孔边环向应力分布。由图知环向应力均为压应力，且呈轴对称分布；孔边不同曲率的弧段应力分布差异较大，拱底的曲率半径最大环向应力最小，拱底和拱腰的结合段曲率半径最小应力值最大。另外，不同弧段曲率越接近意味着顺接越光滑，则洞周环向应力突变值越小，如果两弧段曲率半径差异较大，则存在应力突变，这从拱底和R=1.6m小弧段的结合部可以看出，这也诠释了实际工程中总是避免出现角点的原因。0090 围岩应力分布0091 隧道开挖后洞周围岩产生次生应力场，随距孔边距离的增加，次生应力逐渐趋近于原岩应力。因此，有必要对隧道的拱顶、拱肩、拱腰、拱脚和拱底等五个特征部位附近的围岩应力展开分析。由于所有的应力计算都是在将单位圆外域映射到实际洞室外域的基础上获得的，因此按单位圆外5个相应方向上的点映射到洞室外的实际点考虑，对于数据的提说 明 书CN 102982247 A10。