Nothing Special   »   [go: up one dir, main page]

CN116498391B - 一种地下空间围岩灾害综合预警与辅助决策方法 - Google Patents

一种地下空间围岩灾害综合预警与辅助决策方法 Download PDF

Info

Publication number
CN116498391B
CN116498391B CN202310778679.6A CN202310778679A CN116498391B CN 116498391 B CN116498391 B CN 116498391B CN 202310778679 A CN202310778679 A CN 202310778679A CN 116498391 B CN116498391 B CN 116498391B
Authority
CN
China
Prior art keywords
stress
displacement
boundary
early warning
points
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.)
Active
Application number
CN202310778679.6A
Other languages
English (en)
Other versions
CN116498391A (zh
Inventor
帅永建
谢基祥
田元勋
吕作海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sinohydro Bureau 7 Co Ltd
Original Assignee
Sinohydro Bureau 7 Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinohydro Bureau 7 Co Ltd filed Critical Sinohydro Bureau 7 Co Ltd
Priority to CN202310778679.6A priority Critical patent/CN116498391B/zh
Publication of CN116498391A publication Critical patent/CN116498391A/zh
Application granted granted Critical
Publication of CN116498391B publication Critical patent/CN116498391B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18Special adaptations of signalling or alarm devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/025Geometry of the test
    • G01N2203/0252Monoaxial, i.e. the forces being applied along a single axis of the specimen

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Geology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Computational Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明公开了一种地下空间围岩灾害综合预警与辅助决策方法,属于地下洞室灾害预警领域,所述方法包括:获取洞室的若干个边界点;将所述边界点映射到单位圆上,获取映射函数公式;基于映射函数公式,计算洞室边界上各监测点的应力和位移,所述监测点为多点位移计的布置点;对监测点岩样进行单轴压缩实验,获得单轴压缩应力应变特性曲线,利用复变函数计算应力和位移之间的线性关系,基于应力和位移之间的线性关系来获得每一特征应力对应的特征位移;根据所述特征位移建立基于岩石单轴压缩曲线的变形预警分级体系。本发明可以快速准确地对围岩安全进行预警。

Description

一种地下空间围岩灾害综合预警与辅助决策方法
技术领域
本发明属于地下洞室灾害预警领域,特别是涉及一种地下空间围岩灾害综合预警与辅助决策方法。
背景技术
地下洞室是修建于自然地质体中复杂的地下结构,因此,在施工过程中围岩稳定性受到地应力、岩石力学特性等多种未知因素的影响,导致计算分析结果与实际存在较大差异;而随着工程技术的不断进步,地下洞室的单洞断面尺寸越来越大、开挖边界越来越复杂,会因开挖卸荷、局部应力集中以及应力反复调整诱发开挖面附近的岩体变形失稳,产生片帮、破裂,严重时还会产生松弛垮塌等过度变形或变形不收敛现象。这些不确定因素和变形失稳会导致施工过程中的许多安全隐患。在高地应力环境下,硬岩灾害的反应时间短、频率高、后果严重且预警难度大,因此,根据实时位移监测信息,快速准确地对围岩安全进行预警,建立一个相对合理的变形分级预警体系与辅助决策系统,直接关系到洞室施工的安全和成功。
现有技术中公开了一种适用于岩体条件地下洞室群变形分级预警指标的确定方法(专利公告号:CN110847969B),但是该方法不能针对性的获得某监测点的分级预警指标,只能获得同一岩性和岩体质量条件下的分级预警指标,而且不能快速得到分级预警指标,需要进行开挖支护仿真模拟。
发明内容
本发明的目的在于克服现有技术的不足,提供一种地下空间围岩灾害综合预警与辅助决策方法。
本发明的目的是通过以下技术方案来实现的:一种地下空间围岩灾害综合预警与辅助决策方法,包括:
获取洞室的若干个边界点;
将所述边界点映射到单位圆上,获取映射函数公式;
基于映射函数公式,计算洞室边界上各监测点的应力和位移,所述监测点为多点位移计的布置点;
对监测点岩样进行单轴压缩实验,获得单轴压缩应力应变特性曲线,利用复变函数计算中应力和位移之间的线性关系,基于应力和位移之间的线性关系来获得每一特征应力对应的特征位移;
根据所述特征位移建立基于岩石单轴压缩曲线的变形预警分级体系。
进一步地,通过共形将所述边界点映射到单位圆上。
进一步地,所述映射函数公式为,/>,其中:
(1)
(2)
(3)
式中,表示映射前边界点的坐标,/>表示映射后边界点的坐标,/>为共形映射参数,/>、/>为实数,/>为洛朗级数展开项数,/>为开挖边界和单位圆上对应点的组数,/>为单位圆周线上各点相应于实轴的转角。
进一步地,基于映射函数公式,计算洞室边界上各监测点的应力和位移,包括:
基于复变函数的应力求解公式、复变函数的位移求解公式和映射函数公式,计算洞室边界上各监测点的应力和位移。
进一步地,复变函数的应力求解公式和复变函数的位移求解公式如下:
(4)
(5)
(6)
应力边界的求解公式为:(7)
位移边界的求解公式为:(8)
(9)
(10)
(11)
(12)
式中,、/>分别是直角坐标系下沿/>、/>轴应力分量,/>为/>平面内的剪应力,/>,/>、/>称为复应力函数,在无限平面的任一点/>都解析;/>、/>为沿/>、/>轴位移分量;/>、/>,/>;/>为弹性模量,为泊松比,/>为开挖边界上的点;/>、/>为边界上沿/>、/>轴的面力分量;/>、/>为应力边界的积分起点和终点,/>、/>为关于/>的解析函数,/>、/>为表征应力场的过程参数,假设初始应力场为无穷远处与/>轴呈一角度的均匀应力场,/>、/>为远场应力,为/>与/>轴夹角,/>为/>
通过Cauchy积分算子、留数理论和共轭理论求解出来的复应力函数、/>分别如式(13)~(14)所示:
(13)
(14)式中,/>表示/>中n取值为1时的值,/>表示/>中n取值为2时的值。
进一步地,所述特征应力包括起裂应力、扩容应力和单轴抗压强度。
进一步地,所述特征位移包括起裂位移、扩容位移和破坏位移。
进一步地,所述方法还包括:
针对不同的分级预警,设置相应的工程建议响应措施。
本发明的有益效果是:本发明提出的基于岩体强度的地下空间围岩预警临界变形确定方法,利用复变函数方法,可以快速准确地对围岩安全进行时空预警,利用共形映射可以求解每次开挖后洞室边界的映射函数,进而实现不同开挖部下同一点的预警位移,此为随时间变化的预警;而边界上不同点对应的坐标不同,通过输入不同的坐标可实现对不同点的预警,此为从空间层面上的预警。
附图说明
图1为本发明中发明的一种实施例的流程图;
图2为映射示意图;
图3为边界平面力学模型;
图4为基于岩石单轴压缩曲线的分级体系图;
图5为基于分级体系的响应工作流程图。
具体实施方式
下面将结合实施例,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
参阅图1至图5,本发明提供一种地下空间围岩灾害综合预警与辅助决策方法:
步骤S100.获取洞室的若干个边界点。
所述边界点的取点方法为:在开挖洞室的边界进行取点,取点时直线段可只取首尾两端点,弧线段则适当增加取点数量,以体现弧线段的形状,边界点的坐标为。一般的,在进行弧线段取点时,在可以提现弧线段形状的前提下尽量少取点,可以减少后续的计算量。
步骤S200.将所述边界点映射到单位圆上,映射后边界点的坐标为,获得映射函数公式/>;其中,/>为映射函数公式的项数,根据洞室边界复杂程度而确定;/>的计算公式如式(1)~式(3)。
(1)
(2)
(3)
式中,为共形映射参数,/>、/>为实数,/>为洛朗级数展开项数,/>为开挖边界和单位圆上对应点的组数,/>为单位圆周线上各点相应于实轴的转角。
在一些实施例中,通过共形将所述边界点映射到单位圆上,映射示意如图2所示。
步骤S300.基于复变函数的应力求解公式、复变函数的位移求解公式和映射函数公式,计算洞室边界上各监测点的应力和位移,监测点为实际工程中多点位移计的布置点。
具体的,复变函数的应力求解公式(公式(4)和公式(5))和复变函数的位移求解公式(公式(6))如下:
(4)
(5)
(6)
应力边界的求解公式为:(7)
位移边界的求解公式为: (8)
(9)
(10)
(11)
(12)
式中,、/>分别是直角坐标系下沿/>、/>轴应力分量,/>为/>平面内的剪应力,,/>、/>称为复应力函数,在无限平面的任一点/>都解析;/>、/>为沿/>轴位移分量;/>(平面应变)、/>(平面应力),;/>为弹性模量,/>为泊松比,/>为开挖边界上的点;/>、/>为边界上沿/>轴的面力分量;/>、/>为应力边界的积分起点和终点。/>、/>为关于/>的解析函数,/>、/>为表征应力场的过程参数,假设初始应力场为无穷远处与/>轴呈一角度的均匀应力场,如图3所示。/>、/>为远场应力,/>为/>与/>轴夹角(逆时针为正),/>为/>
通过Cauchy积分算子、留数理论和共轭理论求解出来的复应力函数、/>分别如式(13)~(14)所示。
(13)
(14)式中,/>、/>为/>的前两项,即当/>和/>时。
步骤S400.在开挖卸荷作用下,临空面法向应力降为零,周向应力增加,洞壁岩石处于单轴应力状态,故对监测点岩样进行单轴压缩实验,获得单轴压缩应力应变特性曲线,利用复变函数计算应力和位移之间的线性关系,基于应力和位移之间的线性关系来获得每一特征应力对应的特征位移。
所述特征应力包括起裂应力、扩容应力/>和单轴抗压强度/>
所述特征位移包括起裂位移、扩容位移/>和破坏位移/>
步骤S500.根据所述特征位移建立基于岩石单轴压缩曲线的变形预警分级体系。所述变形预警分级体系如表1和图4所示。
表1
在一些实施例中,针对不同的分级预警,给出相应的工程建议响应措施,如图5所示。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。

Claims (5)

1.一种地下空间围岩灾害综合预警与辅助决策方法,其特征在于,包括:
获取洞室的若干个边界点;
将所述边界点映射到单位圆上,获取映射函数公式;
基于映射函数公式,计算洞室边界上各监测点的应力和位移,所述监测点为多点位移计的布置点;
对监测点岩样进行单轴压缩实验,获得单轴压缩应力应变特性曲线,利用复变函数计算应力和位移之间的线性关系,基于应力和位移之间的线性关系来获得每一特征应力对应的特征位移;
根据所述特征位移建立基于岩石单轴压缩曲线的变形预警分级体系;
所述映射函数公式为,/>,其中:
(1)
(2)
(3)
式中,表示映射前边界点的坐标,/>表示映射后边界点的坐标,/>为共形映射参数,/>、/>为实数,/>为洛朗级数展开项数,/>为开挖边界和单位圆上对应点的组数,/>为单位圆周线上各点相应于实轴的转角;
基于映射函数公式,计算洞室边界上各监测点的应力和位移,包括:
基于复变函数的应力求解公式、复变函数的位移求解公式和映射函数公式,计算洞室边界上各监测点的应力和位移;
复变函数的应力求解公式和复变函数的位移求解公式如下:
(4)
(5)
(6)
应力边界的求解公式为:
(7)
位移边界的求解公式为:
(8)
(9)
(10)
(11)
(12)
式中,、/>分别是直角坐标系下沿/>、/>轴应力分量,/>为/>平面内的剪应力,,/>、/>称为复应力函数,在无限平面的任一点/>都解析;/>、/>为沿、/>轴位移分量;/>,/>;/>为弹性模量,/>为泊松比,/>为开挖边界上的点;/>、/>为边界上沿/>、/>轴的面力分量;/>、/>为应力边界的积分起点和终点,/>、/>为关于/>的解析函数,/>、/>为表征应力场的过程参数,假设初始应力场为无穷远处与/>轴呈一角度的均匀应力场,/>、/>为远场应力,为/>与/>轴夹角,/>为/>
通过Cauchy积分算子、留数理论和共轭理论求解出来的复应力函数、/>分别如式(13)~(14)所示:
(13)
(14)
式中,表示/>中n取值为1时的值,/>表示/>中n取值为2时的值。
2.根据权利要求1所述的一种地下空间围岩灾害综合预警与辅助决策方法,其特征在于,通过共形将所述边界点映射到单位圆上。
3.根据权利要求1所述的一种地下空间围岩灾害综合预警与辅助决策方法,其特征在于,所述特征应力包括起裂应力、扩容应力和单轴抗压强度。
4.根据权利要求1所述的一种地下空间围岩灾害综合预警与辅助决策方法,其特征在于,所述特征位移包括起裂位移、扩容位移和破坏位移。
5.根据权利要求1所述的一种地下空间围岩灾害综合预警与辅助决策方法,其特征在于,所述方法还包括:
针对不同的分级预警,设置相应的工程建议响应措施。
CN202310778679.6A 2023-06-29 2023-06-29 一种地下空间围岩灾害综合预警与辅助决策方法 Active CN116498391B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310778679.6A CN116498391B (zh) 2023-06-29 2023-06-29 一种地下空间围岩灾害综合预警与辅助决策方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310778679.6A CN116498391B (zh) 2023-06-29 2023-06-29 一种地下空间围岩灾害综合预警与辅助决策方法

Publications (2)

Publication Number Publication Date
CN116498391A CN116498391A (zh) 2023-07-28
CN116498391B true CN116498391B (zh) 2023-09-22

Family

ID=87330549

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310778679.6A Active CN116498391B (zh) 2023-06-29 2023-06-29 一种地下空间围岩灾害综合预警与辅助决策方法

Country Status (1)

Country Link
CN (1) CN116498391B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1460354A1 (ru) * 1986-03-10 1989-02-23 Всесоюзный Научно-Исследовательский Горнометаллургический Институт Цветных Металлов Рудничный аккумул тор сжатого воздуха
CN102155231A (zh) * 2011-03-18 2011-08-17 大连海事大学 一种隧道施工过程的快速反馈分析系统
JP2011256525A (ja) * 2010-06-04 2011-12-22 Ohbayashi Corp トンネル監視方法
CN102982247A (zh) * 2012-12-05 2013-03-20 浙江科技学院 一种任意开挖断面隧洞的围岩应力分析方法
CN104361165A (zh) * 2014-11-04 2015-02-18 南京航空航天大学 多孔压电材料电弹场分析方法
CN105117516A (zh) * 2015-07-28 2015-12-02 西北工业大学 一种含孔复合材料层合板孔周应力分布的计算方法
CN106194262A (zh) * 2016-08-03 2016-12-07 太原理工大学 煤矿底板奥灰岩突水在线监测智能预警系统
CN108614035A (zh) * 2018-03-23 2018-10-02 长江水利委员会长江科学院 大埋深洞室硬岩岩爆孕灾风险识别方法
CN108871946A (zh) * 2018-03-23 2018-11-23 长江水利委员会长江科学院 大埋深洞室硬岩岩爆灾变风险等级评估方法
CN109184795A (zh) * 2018-09-03 2019-01-11 煤炭科学技术研究院有限公司 一种煤矿应急救援系统、安全监测方法、装置及存储介质
CN109268028A (zh) * 2018-08-30 2019-01-25 中国水利水电第十四工程局有限公司 一种大跨度地下洞室半幅错距开挖支护方法
CN110847969A (zh) * 2019-10-29 2020-02-28 中国电建集团华东勘测设计研究院有限公司 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4495672B2 (ja) * 2003-07-02 2010-07-07 ニューブレクス株式会社 構造体監視システム
GB2407440B (en) * 2003-09-23 2006-02-22 Responsiveload Ltd Grid stabilising system
CN111368405B (zh) * 2020-02-25 2021-12-28 山东大学 隧道岩体破坏突涌水灾变模拟的近场动力学方法及系统
CN111927558B (zh) * 2020-10-13 2021-01-12 中国科学院武汉岩土力学研究所 动水软弱围岩隧道全断面掘进的安全预警方法及装置

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1460354A1 (ru) * 1986-03-10 1989-02-23 Всесоюзный Научно-Исследовательский Горнометаллургический Институт Цветных Металлов Рудничный аккумул тор сжатого воздуха
JP2011256525A (ja) * 2010-06-04 2011-12-22 Ohbayashi Corp トンネル監視方法
CN102155231A (zh) * 2011-03-18 2011-08-17 大连海事大学 一种隧道施工过程的快速反馈分析系统
CN102982247A (zh) * 2012-12-05 2013-03-20 浙江科技学院 一种任意开挖断面隧洞的围岩应力分析方法
CN104361165A (zh) * 2014-11-04 2015-02-18 南京航空航天大学 多孔压电材料电弹场分析方法
CN105117516A (zh) * 2015-07-28 2015-12-02 西北工业大学 一种含孔复合材料层合板孔周应力分布的计算方法
CN106194262A (zh) * 2016-08-03 2016-12-07 太原理工大学 煤矿底板奥灰岩突水在线监测智能预警系统
CN108614035A (zh) * 2018-03-23 2018-10-02 长江水利委员会长江科学院 大埋深洞室硬岩岩爆孕灾风险识别方法
CN108871946A (zh) * 2018-03-23 2018-11-23 长江水利委员会长江科学院 大埋深洞室硬岩岩爆灾变风险等级评估方法
CN109268028A (zh) * 2018-08-30 2019-01-25 中国水利水电第十四工程局有限公司 一种大跨度地下洞室半幅错距开挖支护方法
CN109184795A (zh) * 2018-09-03 2019-01-11 煤炭科学技术研究院有限公司 一种煤矿应急救援系统、安全监测方法、装置及存储介质
CN110847969A (zh) * 2019-10-29 2020-02-28 中国电建集团华东勘测设计研究院有限公司 一种适用于岩体条件地下洞室群变形分级预警指标的确定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
单位圆外域到任意开挖断面隧洞外域共形映射的计算方法;祝江鸿;杨建辉;施高萍;王;蔡建平;;岩土力学(第01期);全文 *

Also Published As

Publication number Publication date
CN116498391A (zh) 2023-07-28

Similar Documents

Publication Publication Date Title
Khatir et al. Fast simulations for solving fracture mechanics inverse problems using POD-RBF XIGA and Jaya algorithm
Bazilevs et al. Isogeometric fatigue damage prediction in large-scale composite structures driven by dynamic sensor data
Klerck et al. Discrete fracture in quasi-brittle materials under compressive and tensile stress states
Wang et al. Fracture prediction of welded steel connections using traditional fracture mechanics and calibrated micromechanics based models
CN109900554B (zh) 一种利用压痕法计算断裂韧性的方法
Zhang et al. Identifying accurate crack initiation and propagation thresholds in siliceous siltstone and limestone
CN105823661B (zh) 可调控裂纹大小和电导率的模拟应力腐蚀裂纹制备方法
CN111339700B (zh) 核电汽轮机叶片的疲劳损伤评估方法、装置和存储介质
CN105004789B (zh) 一种超声波测量焊接残余应力的误差修正方法
Cisilino et al. Three-dimensional boundary element analysis of fatigue crack growth in linear and non-linear fracture problems
Zhou et al. XFEM based node scheme for the frictional contact crack problem
CN109060539A (zh) 一种岩石微米尺度弹性模量及屈服强度获取方法
Yu et al. Experimental and numerical study on casing wear in a directional well under in situ stress for oil and gas drilling
Chen Advances in conservation laws and energy release rates: Theoretical treatments and applications
Hay et al. Small correction required when applying the Hertzian contact model to instrumented indentation data
Lukasiewicz et al. Calculation of strains in dents based on high resolution in-line caliper survey
CN116498391B (zh) 一种地下空间围岩灾害综合预警与辅助决策方法
Mu et al. Study on rough surfaces: A novel method for high-precision simulation and interface contact performances analysis
Zhao et al. Experimental and numerical studies on small-scale direct tension test for rock
CN108446251B (zh) 一种岩质边坡主动岩石压力的计算方法
Kulesa et al. Low cycle fatigue of steel in strain controled cyclic bending
Horton et al. The inclusion and role of micro mechanical residual stress on deformation of stainless steel type 316L at grain level
Wen et al. Evaluation method of equivalent initial flaw size and fatigue life prediction of nickel-based single crystal superalloy
Vallejos Characterization of large diameter reinforcement under large strain cyclic reversals
CN115221727B (zh) 一种基于含水率的岩体的数值仿真模型参数确定方法

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
GR01 Patent grant
GR01 Patent grant