CN107256301A - 一种模拟岩溶土洞形成演化的三维建模方法 - Google Patents
一种模拟岩溶土洞形成演化的三维建模方法 Download PDFInfo
- Publication number
- CN107256301A CN107256301A CN201710414291.2A CN201710414291A CN107256301A CN 107256301 A CN107256301 A CN 107256301A CN 201710414291 A CN201710414291 A CN 201710414291A CN 107256301 A CN107256301 A CN 107256301A
- Authority
- CN
- China
- Prior art keywords
- soil
- cave
- soil cave
- modeling
- model
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种模拟岩溶土洞形成演化的三维建模方法,充分利用Flac3D软件本身的建模功能和内置编程语言FISH,具体的技术方案:①建立由三个棱锥体组成的正方体;②建立块体外围渐变式网格,其内部大小与步骤1建立的正方体一致;③利用FISH函数,使网格节点成放射状分布5;④对建立模型镜像,完成建模;⑤对模型进行开挖,模拟土洞由小到大演化发育。本发明优点和效果:充分利用了FLAC3D软件本身的建模功能,运用FISH语言编制了简单实用的建模程序,大大简化了建模流程,提高了建模效率。建模精度可根据工程实际情况灵活调整,建模精细度高,实现了岩溶土洞演化发育的三维动态建模方法。
Description
技术领域
本发明涉及岩土数值模拟计算技术领域,具体是涉及一种模拟岩溶土洞形成演化的三维建模方法。
背景技术
岩溶土洞的形成发育与土,水,岩的性质结构密切相关,土洞的形成须有开口向上的溶洞、溶隙或断裂裂隙等可容纳覆盖土体的地下空间。地下水流及自重力等动力因素可导致覆盖层土体剥落破坏,与母体分离并进入地下空间中,于是就形成了近似球体的土洞。一旦破坏开始,土洞便迅速向上扩展,直到覆盖土体应力小于土的强度,土洞便停止发展,处于暂时稳定状态。当动力条件变化时,土洞会继续扩展,最后土洞顶板整体发生垮塌,形成地面塌陷,如图1所示。目前,随着计算机计算技术的飞速发展,数值模拟越来越多的被应用到岩土工程领域中,各种数值模拟软件获得了广阔的应用空间,尤其是FLAC系列软件,以其强大的计算能力成为目前最主要的岩土数值模拟软件。FLAC系列软件的计算能力非常出众,但是建模能力不足,尤其对于较复杂的岩溶土洞形成演化过程的构建。背景技术主要采取两种办法来进行解决,一是对球形土洞形状进行简化,用正方形代替。二是用二维模型建模,但上述方法存在明显的缺点和不足:(1)数值模拟计算结果是否准确的关键之一是构建的数值计算模型与实际模型是否匹配;对岩溶土洞进行不合理的简化或等效,大大的影响了建模的准确性,进而影响了数值计算的准确性。(2)用二维模型建模,不能真实反映土洞的发育演化过程。因此,目前利用FLAC软件进行建模存在数值计算模型与实际模型不匹配的问题。
发明内容
本发明目的是提供一种模拟岩溶土洞形成演化的三维建模方法,充分利用Flac3D软件本身的建模功能和内置编程语言FISH,提高三维建模的精细度和效率,解决背景技术中存在的上述问题。
本发明的技术方案是:一种模拟岩溶土洞形成演化的三维建模方法,包含如下步骤:
S1、确定建模参数,土层厚度(rad),土洞半径(len),网格节点调整参数(rad_size,ab_size):
通过现场地质勘察,在岩溶土洞区域进行原位测试,测量得到所述土层厚度、土洞半径数据;再根据网格节点的偏移量确定网格节点调整参数;
S2、建立模型,并调模型四面体节点位置:
所述建立模型的具体步骤为:
①建立由三个棱锥体组成的正方体;
②建立块体外围渐变式网格,其内部大小与步骤①建立的正方体一致;
③利用FISH函数,使网格节点成放射状分布;
④对建立模型镜像,完成建模;
S3、对模型进行开挖,模拟土洞由小到大演化发育:
将所述土洞顶板概化为半球体,按照模型的网格划分以半球体形状的层体为单位由内向外一层一层开挖。
进一步地,在上述方案中,步骤S2所述建立模型是利用FISH语言编制建模程序。
进一步地,在上述方案中,步骤S2所述建立模型运用的是FLAC3D软件进行建模的。
2、进一步地,在上述方案中,步骤S2所述利用FISH函数,使网格节点成放射状分布的FISH语言如下:
进一步地,在上述方案中,所述对模型进行开挖,是利用null命令,开挖球形区域,令半径=1,模拟土洞形成,然后继续利用null命令,开挖球形区域,令半径=2,以此类推,通过网格划分一层一层由内向外开挖,模拟土洞发展演化。进一步地,所述开挖在具体施工中的操作方法为:
1)生成土体的初始应力场,进行应力平衡计算,生成k0条件下的土体初始应力场;
2)开挖球体最内层土体,在最内层施加开挖面均布力;激活最内层土体的衬砌单元,在最内层土体施加收敛位移;
3)开挖球体第二层土体,移除最内层施加的均布力,在第二层土体中施加开挖面均布力,激活第二层土体的衬砌单元,在第二层土体施加收敛位移;
4)以后各层土体按以上述2)、3)类推,直至球体开挖至最后一层。
更进一步地,在具体施工中,所述开挖主要是根据施工仿真参数,采用钻孔爆破法;采用钻爆法反复循环进行开挖,每掘进一次,工作面按炮孔深度向前推进,完成一次循环作业。
更进一步地,所述施工仿真参数中包含岩溶土洞尺寸信息,以及施工所需循环进尺、机械设备仿真参数;利用不同的参数进行仿真计算会得到相应的仿真成果;通过修改、调整仿真参数可以得到最优的施工进度仿真成果。
本发明的有益效果是:本发明充分利用了FLAC3D软件本身的建模功能,运用FISH语言编制了简单实用的建模程序,大大简化了建模流程,提高了建模效率;建模精度可根据工程实际情况灵活调整,建模精细度高,实现了岩溶土洞演化发育的三维动态建模方法。
附图说明
图1为岩溶土洞发育示意图;
图2为岩溶土洞建模示意图;
图3为三个棱锥体组成的正方体;
图4外围渐变式网格体与正方体组合示意图;
图5使网格节点调整图;
图6土洞形成前模型示意图;
图7土洞演化发育模型示意图(土洞半径1m);
图8土洞演化发育过程模型示意图(土洞半径2m)。
具体实施方式
利用FISH语言编制简单实用的建模程序,具体实施例如下:
一种模拟岩溶土洞形成演化的三维建模方法,包含如下步骤:
(1)定义参数,拟模拟的最大土洞半径rad=4,模型高度和宽度len=5,具体FISH语言如下:
def parm
rad=4;土洞半径
len=5;模型长宽
rad_size=12
ab_size=6
end
parm
(2)运用FLAC3D软件进行建模,先建立由三个棱锥体组成的正方体,正方体的边长为rad=4,如图3。再建立块体外围渐变式网格,外部边长为len=5,其内部大小与由三个棱锥体组成的正方体一致,如图4。具体FISH语言如下:
gen zone pyramid p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 p4 radrad rad size rad_size rad_size rad_size group 1;建立棱锥体1
gen zone pyramid p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 p4 radrad rad size rad_size rad_size rad_size group 2;建立棱锥体2
gen zone pyramid p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 p4 radrad rad size rad_size rad_size rad_size group 3;建立棱锥体3
gen zone radbrick edge len size rad_size rad_size rad_size rad_sizerat 1 1 1 1.2 dim rad rad rad group 4;块体外围渐变放射网格
(3)FISH函数,使网格节点成放射状分布,如图5,具体FISH语言如下:
(4)对建立模型镜像,完成建模,如图6,具体FISH语言如下:
gen zone ref dip 90 dd 90
(5)利用null命令,开挖球形区域(半径=1),模拟土洞形成,如图7,具体FISH语言如下:
model null range annulus center 0 0 0 radius 0 1
(6)利用null命令,开挖球形区域(半径=2),模拟土洞发展演化,如图,8,具体FISH语言如下:
model null range annulus center 0 0 0 radius 1 2
其中,确定建模参数是通过现场地质勘察,在岩溶土洞区域进行原位测试,测量得到所述土层厚度、土洞半径数据;再根据网格节点的偏移量确定网格节点调整参数。
开挖在具体施工中的操作方法为:
1)生成土体的初始应力场,进行应力平衡计算,生成k0条件下的土体初始应力场;
2)开挖球体最内层土体,在最内层施加开挖面均布力;激活最内层土体的衬砌单元,在最内层土体施加收敛位移;
3)开挖球体第二层土体,移除最内层施加的均布力,在第二层土体中施加开挖面均布力,激活第二层土体的衬砌单元,在第二层土体施加收敛位移;
4)以后各层土体按以上述2)、3)类推,直至球体开挖至最后一层。
开挖主要是根据施工仿真参数,采用钻孔爆破法;采用钻爆法反复循环进行开挖,每掘进一次,工作面按炮孔深度向前推进,完成一次循环作业。施工仿真参数中包含岩溶土洞尺寸信息,以及施工所需循环进尺、机械设备仿真参数;利用不同的参数进行仿真计算会得到相应的仿真成果。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。
Claims (8)
1.一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,包括以下步骤:
S1、确定建模参数,土层厚度,土洞半径,网格节点调整参数:
通过现场地质勘察,在岩溶土洞区域进行原位测试,测量得到所述土层厚度、土洞半径数据;再根据网格节点的偏移量确定网格节点调整参数;
S2、建立模型,并调模型四面体节点位置:
所述建立模型的具体步骤为:
①建立由三个棱锥体组成的正方体;
②建立块体外围渐变式网格,其内部大小与步骤①建立的正方体一致;
③利用FISH函数,使网格节点成放射状分布;
④对建立模型镜像,完成建模;
S3、对模型进行开挖,模拟土洞由小到大演化发育:
将所述土洞顶板概化为半球体,按照模型的网格划分以半球体形状的层体为单位由内向外一层一层开挖。
2.根据权利要求1所述的一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,步骤S2所述建立模型是利用FISH语言编制建模程序。
3.根据权利要求1所述的一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,步骤S2所述建立模型运用的是FLAC3D软件进行建模的。
4.根据权利要求1所述的一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,步骤S2所述利用FISH函数,使网格节点成放射状分布的FISH语言如下:
5.根据权利要求1所述的一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,所述对模型进行开挖,是利用null命令,开挖球形区域,令半径=1,模拟土洞形成,然后继续利用null命令,开挖球形区域,令半径=2,以此类推,通过网格划分一层一层由内向外开挖,模拟土洞发展演化。
6.根据权利要求5所述的一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,所述开挖在具体施工中的操作方法为:
1)生成土体的初始应力场,进行应力平衡计算,生成k0条件下的土体初始应力场;
2)开挖球体最内层土体,在最内层施加开挖面均布力;激活最内层土体的衬砌单元,在最内层土体施加收敛位移;
3)开挖球体第二层土体,移除最内层施加的均布力,在第二层土体中施加开挖面均布力,激活第二层土体的衬砌单元,在第二层土体施加收敛位移;
4)以后各层土体按以上述2)、3)类推,直至球体开挖至最后一层。
7.根据权利要求6所述的一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,在具体施工中,所述开挖主要是根据施工仿真参数,采用钻孔爆破法;采用钻爆法反复循环进行开挖,每掘进一次,工作面按炮孔深度向前推进,完成一次循环作业。
8.根据权利要求7所述的一种模拟岩溶土洞形成演化的三维建模方法,其特征在于,所述施工仿真参数中包含岩溶土洞尺寸信息,以及施工所需循环进尺、机械设备仿真参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710414291.2A CN107256301A (zh) | 2017-06-05 | 2017-06-05 | 一种模拟岩溶土洞形成演化的三维建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710414291.2A CN107256301A (zh) | 2017-06-05 | 2017-06-05 | 一种模拟岩溶土洞形成演化的三维建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107256301A true CN107256301A (zh) | 2017-10-17 |
Family
ID=60023544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710414291.2A Pending CN107256301A (zh) | 2017-06-05 | 2017-06-05 | 一种模拟岩溶土洞形成演化的三维建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107256301A (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101739716A (zh) * | 2009-12-04 | 2010-06-16 | 东北大学 | 工程岩体三维空间结构建模与关键块识别方法 |
CN102945570A (zh) * | 2012-11-23 | 2013-02-27 | 华东师范大学 | 一种全空间三维数字地球模型的构建方法 |
-
2017
- 2017-06-05 CN CN201710414291.2A patent/CN107256301A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101739716A (zh) * | 2009-12-04 | 2010-06-16 | 东北大学 | 工程岩体三维空间结构建模与关键块识别方法 |
CN102945570A (zh) * | 2012-11-23 | 2013-02-27 | 华东师范大学 | 一种全空间三维数字地球模型的构建方法 |
Non-Patent Citations (1)
Title |
---|
贾龙 等: "岩溶土洞演化及其数值模拟分析", 《中国岩溶》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cai | Influence of stress path on tunnel excavation response–Numerical tool selection and modeling strategy | |
Geyer et al. | A short review of our current understanding of the development of ring faults during collapse caldera formation | |
CN105787206A (zh) | 一种公路隧道锚杆支护设计方法 | |
CN108109203B (zh) | 基于bim环境下的地质界线点云数据提取方法及系统 | |
Fattahi et al. | Hybrid Monte Carlo simulation and ANFIS-subtractive clustering method for reliability analysis of the excavation damaged zone in underground spaces | |
CN106227913A (zh) | 用于水电工程的堆石料爆破级配预测方法 | |
CN116975988A (zh) | 一种基于数字孪生模型的隧道施工安全实时分析方法 | |
He et al. | Overburden Damage Degree‐Based Optimization of High‐Intensity Mining Parameters and Engineering Practices in China’s Western Mining Area | |
CN112487517B (zh) | 浅埋盾构隧道围岩压力的二维数值求解方法 | |
CN107256301A (zh) | 一种模拟岩溶土洞形成演化的三维建模方法 | |
CN111046471B (zh) | 一种帷幕灌浆三维可视化模型构建方法 | |
CN113010997A (zh) | 一种岩溶发育区桩基冲击成孔施工事故的预警处治方法 | |
CN104915477B (zh) | 一种矿用液压挖掘机挖掘阻力的仿真方法 | |
CN111524225B (zh) | 一种基于bim技术的人工挖孔灌注桩施工方法 | |
Founta et al. | Numerical simulation of ground movements due to EPB tunnelling in clay | |
Zou et al. | Development and application of an intelligent evaluation and control platform for tunnel smooth blasting | |
Alsahly et al. | Three dimensional re-meshing for real time modeling of advancing process in mechanized tunneling | |
Niu et al. | An intelligent prediction method of the karst curtain grouting volume based on support vector machine | |
Alsahly et al. | Adaptive computational simulation of TBM-soil interactions during machine-driven tunnel construction in saturated soft soils | |
CN104636511A (zh) | 基于Dasarathy模型的地铁施工风险快速评估方法 | |
Campbell et al. | Integrated slope stability assessment in a complex geotechnical and hydrogeological setting | |
Nikakhtar et al. | Performance comparison of particle swarm optimization and genetic algorithm for back-analysis of soil layer geotechnical parameters | |
CN107103830A (zh) | 隧道开挖模型试验的模拟方法及所用工具 | |
CN102880774A (zh) | 基于三维离散填充体的三下压煤开采沉陷治理分析方法 | |
Liu et al. | The Comparison and Analysis between Numerical Simulation and Monitor Data of Deep Foundation Pit Excavation Deformation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171017 |
|
RJ01 | Rejection of invention patent application after publication |