CN105956376B - 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 - Google Patents
一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 Download PDFInfo
- Publication number
- CN105956376B CN105956376B CN201610261281.5A CN201610261281A CN105956376B CN 105956376 B CN105956376 B CN 105956376B CN 201610261281 A CN201610261281 A CN 201610261281A CN 105956376 B CN105956376 B CN 105956376B
- Authority
- CN
- China
- Prior art keywords
- model scope
- displacement error
- model
- formula
- tunnel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明提供了一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法,包括步骤:建立不同模型范围下的数值计算模型,并进行隧洞开挖数值模拟;计算位移误差并提取整体位移误差,采用最小二乘法分析不同方向的模型范围对整体位移误差的影响,探讨整体位移误差与模型范围的函数关系;基于整体位移误差和模型范围函数关系,限定某一整体位移误差允许值并以模型范围最小为目标确定模型范围最优值。本发明能同时考虑不同方向的模型尺寸影响,从而确定最优的模型范围,可广泛应用于深埋隧洞施工开挖数值模拟中,为深埋隧洞数值模型的建立提供指导。
Description
技术领域
本发明涉及深埋隧洞数值模拟分析领域,具体地说是一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法。
背景技术
数值模拟是指依靠计算机,通过数值计算方法,对工程问题和物理问题乃至自然界各类问题进行研究。在深埋隧洞数值模拟分析中,由于左右和下部范围是无限的,其范围取值已是数值模拟计算中不可回避的问题。过大的模型范围将导致庞大的计算量,影响计算效率,而过小的模型范围又将产生较大的边界效应,影响数值模拟的精度,因此截取合适范围的岩体为研究对象,建立数值计算模型具有重要的指导意义和实践价值。
深埋隧洞模型范围的选取分为纵向和横向两个部分。目前纵向模型范围的取值尚缺乏研究,人们一般根据经验取5~20倍洞径;而对于横向模型范围的选取,人们通常将模型横截面内的上部、下部和左右模型尺寸等同起来,不加以区别,然后根据经验取3~10倍的洞径。
因此现有的数值模型范围选取存在问题:(1)模型范围的选取大多根据经验,取值不统一,取值原则不明确;(2)模型范围选定以后,无法确定选定的模型范围对计算结果的影响程度;(3)横向模型范围选取中,未对各个方向的模型尺寸加以区分。
发明内容
本发明的目的是克服现有技术存在的不足,提供一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法,该方法可评价数值模型的计算精度,同时可以考虑模型各个方向尺寸影响。
本发明一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法,包括如下步骤:
步骤1,建立数值计算模型,在不同的工况下进行隧洞开挖数值模拟,所述工况包括不同的岩性和不同的地应力;对于纵向模型范围的选取,建立横向模型范围为定值、纵向模型范围足够大的模型,纵向模型范围分为近端边界到监测断面的距离H0和远端边界到监测断面的距离H4;对于横向模型范围的选取,建立纵向模型范围为定值、横向模型范围不同的模型,横向模型范围包括上部范围H1、下部范围H2和左右范围H3;
步骤2,基于步骤1中的数值模拟结果,计算位移误差并提取整体位移误差,采用最小二乘法分析不同方向的模型范围对整体位移误差的影响,计算整体位移误差与模型范围的函数关系;位移数据包括隧洞顶部、隧洞腰部和隧洞底部三个部位的径向位移;对于纵向模型范围的确定,整体位移误差指监测断面上的最大位移误差;对于横向模型范围的确定,整体位移误差指监测断面上的平均位移误差;
步骤3,基于步骤2计算出的整体位移误差和模型范围的函数关系,限定某一整体位移误差允许值并以模型范围最小为目标确定模型范围最优值,对于纵向模型范围的确定,模型范围最小指纵向模型尺寸最小;对于横向模型范围的确定,模型范围最小指横截面面积最小。
1、纵向模型范围的确定:
监测断面上各监测点的位移误差wi采用公式(1)
式中,si为监测断面上某监测点的位移,包括顶部、腰部和底部三个部位;sa为模型取足够大时,中部范围内监测点位移平均值,即没有边界效应影响的稳定值,同样包括顶部、腰部和底部三个部位;
纵向模型范围选取时,整体位移误差取监测断面的最大位移误差wm采用公式(2)
wm=max(wA,wB,wC) (2)
式中,wm为整体位移误差,作为拟合分析的基础数据;wA、wB和wC分别为隧洞顶部、腰部和底部监测点的位移误差;
纵向模型范围与监测断面的整体位移误差的函数关系采用公式(3)
wm=a(H0/D)b+c(H4/D)d+e (3)
式中,H0、H4分别为近端边界与监测断面的距离以及远端边界与监测断面的距离;D为隧洞的开挖洞径;a、b、c、d、e为用最小二乘法拟合的常数。
最优纵向模型范围的选取:
式中,Hop为最优纵向模型范围,相当于纵向模型尺寸最小;[w]为某一既定整体位移误差允许值。
2、横向模型范围的确定:
横向模型范围选取时,整体位移误差取平均位移误差采用公式(5)
式中,为平均位移误差,作为拟合分析的基础数据;wA、wB和wC分别为顶部、腰部和底部监测点的位移误差。
横向模型范围与监测断面平均位移误差的函数关系采用公式(6)
式中,H1、H2、H3分别表示模型的上部范围、下部范围和左右范围;
D为隧洞的开挖洞径;a、b、c、d、j、f、g为用最小二乘法拟合的常数。
最优横向模型范围的确定为:
式中,Aop为最优横向模型范围,相当于模型的横截面面积最小;[w]为某一既定整体位移误差允许值。
与现有技术相比,本发明的有益效果在于:
1、采用本发明方法,可以在保证计算精度的前提下,很方便的确定最优的计算模型范围,提高计算效率,解决了现阶段依靠经验进行模型范围选取的问题,避免了人为选取的主观性;
2、引入数学方法,可以定量的评价采用某计算模型时,模拟结果的位移误差,提出了一种全新的数值模拟精度评价方法;
3、本发明将位移误差与各个方向的模型范围相关联,可以考虑不同方向的模型范围对计算结果的影响,更符合实际情况。
附图说明
图1为本发明方法的流程图。
图2为本发明的模型示意图。
图中,1为监测断面,2为开挖隧洞,3为近端边界,4为远端边界,5为远端边界到监测断面的距离H4,6为近端边界到监测断面的距离H0,7为上部模型范围H1,8为下部模型范围H2,9为左右模型范围H3。
具体实施方式
下面结合附图对本发明进行进一步说明。
如图1、图2所示,一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法,具体步骤如下:
步骤1,建立数值计算模型,在不同的工况下进行隧洞开挖数值模拟,具体工况包括不同的岩性和不同的地应力。对于纵向模型范围的选取,建立横向模型范围一定,纵向模型范围足够大的模型,纵向模型范围分为近端边界到监测断面的距离H0和远端边界到监测断面的距离H4;对于横向模型范围的选取,保证纵向模型范围一定,改变横向范围建立不同的模型,横向模型范围包括上部范围H1、下部范围H2和左右范围H3。
步骤2,基于步骤1中的数值模拟结果,计算位移误差并提取整体位移误差,进而引入最小二乘法分析不同方向的模型范围对整体位移误差的影响,探讨整体位移误差与模型范围的函数关系,位移数据包括隧洞顶部、隧洞腰部和隧洞底部三个部位的径向位移。对于纵向模型范围的确定,整体位移误差指监测断面上的最大位移误差;对于横向模型范围的确定,整体位移误差指监测断面上的平均位移误差。
步骤3,基于步骤2计算出的整体位移误差和模型范围的函数关系,限定某一整体位移误差允许值并以模型范围最小为目标确定模型范围最优值,对于纵向模型范围的确定,模型范围最小指纵向模型尺寸最小;对于横向模型范围的确定,模型范围最小是指横截面面积最小。
针对纵向和横向模型范围采用不同的方法进行确定:
(1)纵向模型范围的确定:
监测断面上各监测点的位移误差wi采用公式(1)
式(1)中:
si为监测断面上某监测点的位移,包括顶部、腰部和底部三个部位;
sa为模型取足够大时,中部范围内监测点位移平均值,即没有边界效应影响的稳定值,同样包括顶部、腰部和底部三个部位。
纵向模型范围选取时,整体位移误差取监测断面的最大位移误差wm采用公式(2)
wm=max(wA,wB,wC) (2)
式(2)中:
wm为整体位移误差,作为拟合分析的基础数据;
wA、wB和wC分别为隧洞顶部、腰部和底部监测点的位移误差。
纵向模型范围与监测断面的整体位移误差的函数关系采用公式(3)
wm=a(H0/D)b+c(H4/D)d+e (3)
式(3)中:
H0、H4分别为近端边界与监测断面的距离以及远端边界与监测断面的距离;
D为隧洞的开挖洞径;
a、b、c、d、e为用最小二乘法拟合的常数。
最优纵向模型范围的选取:
式(4)中:
Hop为最优纵向模型范围,相当于纵向模型尺寸最小;
[w]为某一既定整体位移误差允许值。
(2)横向模型范围的确定:
横向模型范围选取时,整体位移误差取平均位移误差采用公式(5)
式(5)中:
为平均位移误差,作为拟合分析的基础数据;
wA、wB和wC分别为顶部、腰部和底部监测点的位移误差。
横向模型范围与监测断面平均位移误差的函数关系采用公式(6)
式(6)中:
H1、H2、H3分别表示模型的上部范围、下部范围和左右范围;
D为隧洞的开挖洞径;
a、b、c、d、j、f、g为用最小二乘法拟合的常数。
最优横向模型范围的确定为:
式(7)中:
Aop为最优横向模型范围,相当于模型的横截面面积最小;
[w]为某一既定整体位移误差允许值。
Claims (3)
1.一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法,其特征在于包括如下步骤:
步骤(1),建立数值计算模型,在不同的工况下进行隧洞开挖数值模拟,所述工况包括不同的岩性和不同的地应力;对于纵向模型范围的选取,建立横向模型范围为定值、纵向模型范围足够大的模型,纵向模型范围分为近端边界到监测断面的距离H0和远端边界到监测断面的距离H4;对于横向模型范围的选取,建立纵向模型范围为定值、横向模型范围不同的模型,横向模型范围包括上部范围H1、下部范围H2和左右范围H3;
步骤(2),基于步骤(1)中的数值模拟结果,计算位移误差并提取整体位移误差,采用最小二乘法分析不同方向的模型范围对整体位移误差的影响,计算整体位移误差与模型范围的函数关系;位移数据包括隧洞顶部、隧洞腰部和隧洞底部三个部位的径向位移;对于纵向模型范围的确定,整体位移误差指监测断面上的最大位移误差;对于横向模型范围的确定,整体位移误差指监测断面上的平均位移误差;
步骤(3),基于步骤(2)计算出的整体位移误差和模型范围的函数关系,限定某一整体位移误差允许值并以模型范围最小为目标确定模型范围最优值,对于纵向模型范围的确定,模型范围最小指纵向模型尺寸最小;对于横向模型范围的确定,模型范围最小指横截面面积最小。
2.如权利要求1所述的模型范围取值方法,其特征在于所述纵向模型范围的确定方法如下:
监测断面上各监测点的位移误差wi采用公式(1)
式中,si为监测断面上某监测点的位移,包括顶部、腰部和底部三个部位;sa为模型取足够大时,中部范围内监测点位移平均值,即没有边界效应影响的稳定值,同样包括顶部、腰部和底部三个部位;
纵向模型范围选取时,整体位移误差取监测断面的最大位移误差wm采用公式(2)
wm=max(wA,wB,wC) (2)
式中,wm为整体位移误差,作为拟合分析的基础数据;wA、wB和wC分别为隧洞顶部、腰部和底部监测点的位移误差;
纵向模型范围与监测断面的整体位移误差的函数关系采用公式(3)
wm=a(H0/D)b+c(H4/D)d+e (3)
式中,H0、H4分别为近端边界与监测断面的距离以及远端边界与监测断面的距离;D为隧洞的开挖洞径;a、b、c、d、e为用最小二乘法拟合的常数;
最优纵向模型范围的选取:
式中,Hop为最优纵向模型范围,相当于纵向模型尺寸最小;[w]为某一既定整体位移误差允许值。
3.如权利要求1所述的模型范围取值方法,其特征在于所述横向模型范围的确定方法如下:
横向模型范围选取时,整体位移误差取平均位移误差采用公式(5)
式中,为平均位移误差,作为拟合分析的基础数据;wA、wB和wC分别为顶部、腰部和底部监测点的位移误差;
横向模型范围与监测断面平均位移误差的函数关系采用公式(6)
式中,H1、H2、H3分别表示模型的上部范围、下部范围和左右范围;
D为隧洞的开挖洞径;a、b、c、d、j、f、g为用最小二乘法拟合的常数;
最优横向模型范围的确定为:
式中,Aop为最优横向模型范围,相当于模型的横截面面积最小;[w]为某一既定整体位移误差允许值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610261281.5A CN105956376B (zh) | 2016-04-25 | 2016-04-25 | 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610261281.5A CN105956376B (zh) | 2016-04-25 | 2016-04-25 | 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105956376A CN105956376A (zh) | 2016-09-21 |
CN105956376B true CN105956376B (zh) | 2018-07-27 |
Family
ID=56915731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610261281.5A Expired - Fee Related CN105956376B (zh) | 2016-04-25 | 2016-04-25 | 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105956376B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106682318B (zh) * | 2016-12-29 | 2019-12-10 | 武汉大学 | 一种基于隧洞围岩纵向变形曲线的初次支护时机计算方法 |
CN106919745B (zh) * | 2017-02-23 | 2019-12-10 | 武汉大学 | 隧洞管片衬砌内力分析方法 |
CN113987632B (zh) * | 2021-10-19 | 2023-07-21 | 武汉大学 | 预测场地液化横向位移的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102982247A (zh) * | 2012-12-05 | 2013-03-20 | 浙江科技学院 | 一种任意开挖断面隧洞的围岩应力分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100981322B1 (ko) * | 2007-10-09 | 2010-09-10 | 재단법인서울대학교산학협력재단 | 터널굴착기의 최적 설계모델에 의한 터널굴착기 설계방법 |
-
2016
- 2016-04-25 CN CN201610261281.5A patent/CN105956376B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102982247A (zh) * | 2012-12-05 | 2013-03-20 | 浙江科技学院 | 一种任意开挖断面隧洞的围岩应力分析方法 |
Non-Patent Citations (2)
Title |
---|
深部实验隧洞围岩脆性破坏及数值模拟;张传庆 等;《岩石力学与工程学报》;20101031;第29卷(第10期);第2063-2068页 * |
隧道开挖面附近约束损失分析;唐雄俊 等;《华中科技大学学报》;20091231;第26卷(第4期);第49-66页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105956376A (zh) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105956376B (zh) | 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 | |
CN103411607B (zh) | 行人步长估计及航位推算方法 | |
CN105353620B (zh) | 一种基于地表水深信息的地面灌溉控制方法 | |
CN104881547B (zh) | 一种用于定向井井眼轨迹的误差分析方法 | |
CN106709207A (zh) | 一种考虑粗糙表面微凸体相互作用影响的确定受载结合部法向接触刚度的方法 | |
CN104750884B (zh) | 基于多因素非线性回归的页岩油气富集指数定量评价方法 | |
WO2014059282A3 (en) | Attitude reference for tieback/overlap processing | |
CN109948272A (zh) | 基于井间连通性的调堵动态预测方法和系统 | |
CN109405817A (zh) | 一种圆形立柱垂直度检测方法 | |
CN106446349A (zh) | 一种基于revit软件的土方原地摊销的计算方法 | |
CN102322256A (zh) | 水平井着陆导向方法 | |
CN104265279B (zh) | 断层条件下随钻测井曲线预测方法 | |
GB2575418A (en) | System and methods for evaluating a formation using pixelated solutions of formation data | |
CN103759701A (zh) | 基于安卓平台的手机智能测距方法 | |
FR3046479B1 (fr) | Extrapolation de la permeabilite effective d'un reseau de fractures discretes | |
CN105005646B (zh) | 一种边坡稳定安全分析全息法 | |
CN105953794A (zh) | 一种基于mems传感器的计步导航方法 | |
CN106484930A (zh) | 用于确定非均质油藏复杂结构井产量的方法及系统 | |
CN104504236A (zh) | 泵功图计量稳定性新方法 | |
CN108062631A (zh) | 城市内涝风险评估方法、装置及终端 | |
RU2015134392A (ru) | Способ моделирования подземного объема | |
CN107945261A (zh) | 一种地质钻孔测斜计算方法 | |
CN106703842A (zh) | 盾构隧道衬砌结构渗透系数计算方法 | |
CN112943198B (zh) | 一种深层页岩复杂构造地层非均匀应力场计算方法 | |
CN107883858A (zh) | 大型深坑回填方量的测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180727 |