CN105956376A - 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 - Google Patents
一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 Download PDFInfo
- Publication number
- CN105956376A CN105956376A CN201610261281.5A CN201610261281A CN105956376A CN 105956376 A CN105956376 A CN 105956376A CN 201610261281 A CN201610261281 A CN 201610261281A CN 105956376 A CN105956376 A CN 105956376A
- Authority
- CN
- China
- Prior art keywords
- model scope
- model
- displacement error
- formula
- longitudinal
- 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.)
- Granted
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 true CN105956376A (zh) | 2016-09-21 |
CN105956376B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106682318A (zh) * | 2016-12-29 | 2017-05-17 | 武汉大学 | 一种基于隧洞围岩纵向变形曲线的初次支护时机计算方法 |
CN106919745A (zh) * | 2017-02-23 | 2017-07-04 | 武汉大学 | 隧洞管片衬砌内力分析方法 |
CN113987632A (zh) * | 2021-10-19 | 2022-01-28 | 武汉大学 | 预测场地液化横向位移的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090036503A (ko) * | 2007-10-09 | 2009-04-14 | 한국건설기술연구원 | 터널굴착기의 최적 설계모델에 의한 터널굴착기 설계방법 |
CN102982247A (zh) * | 2012-12-05 | 2013-03-20 | 浙江科技学院 | 一种任意开挖断面隧洞的围岩应力分析方法 |
-
2016
- 2016-04-25 CN CN201610261281.5A patent/CN105956376B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090036503A (ko) * | 2007-10-09 | 2009-04-14 | 한국건설기술연구원 | 터널굴착기의 최적 설계모델에 의한 터널굴착기 설계방법 |
CN102982247A (zh) * | 2012-12-05 | 2013-03-20 | 浙江科技学院 | 一种任意开挖断面隧洞的围岩应力分析方法 |
Non-Patent Citations (2)
Title |
---|
唐雄俊 等: "隧道开挖面附近约束损失分析", 《华中科技大学学报》 * |
张传庆 等: "深部实验隧洞围岩脆性破坏及数值模拟", 《岩石力学与工程学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106682318A (zh) * | 2016-12-29 | 2017-05-17 | 武汉大学 | 一种基于隧洞围岩纵向变形曲线的初次支护时机计算方法 |
CN106682318B (zh) * | 2016-12-29 | 2019-12-10 | 武汉大学 | 一种基于隧洞围岩纵向变形曲线的初次支护时机计算方法 |
CN106919745A (zh) * | 2017-02-23 | 2017-07-04 | 武汉大学 | 隧洞管片衬砌内力分析方法 |
CN106919745B (zh) * | 2017-02-23 | 2019-12-10 | 武汉大学 | 隧洞管片衬砌内力分析方法 |
CN113987632A (zh) * | 2021-10-19 | 2022-01-28 | 武汉大学 | 预测场地液化横向位移的方法 |
CN113987632B (zh) * | 2021-10-19 | 2023-07-21 | 武汉大学 | 预测场地液化横向位移的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105956376B (zh) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hickey | Slope angle and slope length solutions for GIS | |
CN105740563B (zh) | 一种成熟油田二次开发之优势通道识别方法 | |
WO2016192077A1 (zh) | 一种致密气压裂水平井数值试井模型建立求解方法 | |
CN105426619B (zh) | 一种确定盾构隧道纵向弯曲刚度的分析方法 | |
CN106709207A (zh) | 一种考虑粗糙表面微凸体相互作用影响的确定受载结合部法向接触刚度的方法 | |
CN105956376A (zh) | 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 | |
CN104252576B (zh) | 一种地铁隧道岩土体等效计算参数反演方法 | |
CN105205864B (zh) | 基于多源数据的地质结构面三维模型自动建模方法和系统 | |
CN101847171A (zh) | 基于安全监控的边坡位移反分析方法 | |
CN104933268A (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
CN111382472A (zh) | 随机森林融合svm预测盾构引起近接结构变形方法及装置 | |
CN109653728B (zh) | 一种基于向量相似度的井眼轨迹钻前模拟方法 | |
CN112227292A (zh) | 基于一维水动力模型的空腹拱桥壅水智能分析方法 | |
CN106897554A (zh) | 一种基于matlab进行p‑y曲线的自动化计算方法 | |
RU2015134392A (ru) | Способ моделирования подземного объема | |
CN104504236B (zh) | 泵功图计量稳定性方法 | |
KR20110072551A (ko) | 텐서형 와점성계수를 가진 2차원 하천흐름모형을 이용하여 천수흐름을 해석하는 방법 | |
CN110847112B (zh) | 一种基于水力学模拟的河道泄洪预警方法 | |
CN116882023A (zh) | 一种地铁隧道暗挖施工上覆土层横向沉降预测方法 | |
CN102288144B (zh) | 确定含水层平均厚度的定量方法 | |
CN107832482A (zh) | 致密储层多尺度裂缝网络建模及模拟方法 | |
CN106971422A (zh) | 一种基于catia的拱坝坝肩稳定分析模型编绘方法 | |
CN107944127A (zh) | 一种井工开采地表沉陷参数拐点偏距的确定方法 | |
CN109271734B (zh) | 一种基于机器学习技术的场地地震液化势评价方法 | |
CN109002580A (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 |