CN114722614B - 一种基于二维-三维协同的***动力响应模拟方法 - Google Patents
一种基于二维-三维协同的***动力响应模拟方法 Download PDFInfo
- Publication number
- CN114722614B CN114722614B CN202210404710.5A CN202210404710A CN114722614B CN 114722614 B CN114722614 B CN 114722614B CN 202210404710 A CN202210404710 A CN 202210404710A CN 114722614 B CN114722614 B CN 114722614B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- model
- blast hole
- explosive
- blasting
- 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
Links
Images
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于二维‑三维协同的***动力响应模拟方法,属于隧道或巷道钻孔***数值模拟领域。包括:根据现场钻孔***参数建立二维平面应变模型,设置边界条件和材料参数,计算出各段别炮孔壁上的***荷载压力时程曲线,并结合理论及经验公式进行验算,再将所得***荷载曲线继承施加至三维真实尺度模型的炮孔壁上,根据现场实测***振动数据,进一步调试***及岩体材料与模型参数,最后进行三维真实尺度模型***动力响应分析。本发明适用于水利、采矿、交通及市政等岩土***领域大型复杂模型的计算分析,相比传统数值模拟方法更高效稳定且结果更加真实可靠。
Description
技术领域
本发明属于隧道或巷道钻孔***数值模拟领域,具体涉及一种基于二维-三维协同的***动力响应模拟方法。
背景技术
随着我国基础设施的逐渐发展,如近接隧道,大型地下洞室等复杂工程不断涌现。
由于数值模拟技术计算精度和可靠性高,其计算结果已经成为各类工程问题分析的依据。在隧道或巷道***动力响应模拟中,关于***荷载的模拟方法主要有2类:一类是直接建立***单元,定义***参数和状态方程,并通过流固耦合算法实现***-岩石的荷载传递;另一类是将爆源荷载简化为函数,再将其施加于模型炮孔壁、等效边界或炮孔联心线上。其中,第一类方法可以较真实地模拟***的爆轰过程,但其对***及周围介质的网格划分有严格的要求,若网格尺寸划分质量较差,计算结果易失真,而若网格划分过于精细,则单元数量显著增加,需消耗大量计算时间,严重时会使计算停滞;第二类方法将***荷载进行简化,省去了***爆轰过程的模拟,虽提高了计算效率,但其计算精度及可靠性无法保证。
发明内容
为了克服上述现有技术存在的不足,本发明提供了一种基于二维-三维协同的***动力响应模拟方法。
为了实现上述目的,本发明提供如下技术方案:
一种基于二维-三维协同的***动力响应模拟方法,包括:
基于现场各段别炮孔***参数和装药结构,建立二维平面应变模型;
设置二维平面应变模型的边界条件;设置二维平面应变模型中***和空气参数及状态方程、围岩的本构模型和物理力学参数;
对二维平面应变模型进行求解,得到各段别炮孔的***压力荷载曲线;
基于圣维南原理利用实际工程形状、尺寸和研究区域围岩类别,建立三维真实尺度模型;
设置三维真实尺度模型的边界条件;设置三维真实尺度模型中围岩、初期支护和二次衬砌的本构模型和物理力学参数;
将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔;
求解施加***荷载压力曲线的三维真实尺度模型得到三维真实尺度模型二次衬砌的振动波形,利用二次衬砌的振动波形验证三维真实尺度模型及选取参数的可靠性,根据三维真实尺度模型分析***动力响应。
优选的,还包括:
对二维平面应变模型进行求解时,直至各段别炮孔的***荷载压力曲线的荷载峰值符合理论和经验公式、且三维真实尺度模型二次衬砌的振动波形与现场监测振动数据吻合,得到各段别炮孔的***压力荷载曲线;否则重复设置二维平面应变模型中***和空气参数及状态方程、围岩的本构模型和物理力学参数,对二维平面应变模型进行求解。
优选的,当炮孔的***荷载压力曲线的荷载峰值与理论计算值的相对误差小于10%时,判定炮孔的***荷载压力曲线的荷载峰值符合理论和经验公式;
其中,耦合装药时,炮孔的***荷载压力理论计算公式为:
不耦合装药时,炮孔的***荷载压力理论计算公式为:
优选的,所述将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔的步骤包括:
根据现场******段别分别建立part,将不同段别的***荷载压力时程曲线按照起爆时间依次施加至对应的part。
优选的,还包括:
对三维真实尺度模型进行求解时,得到三维数值模型二次衬砌的振动波形,直至三维真实尺度模型二次衬砌的振动波形与现场监测振动数据吻合;否则重新设置三维真实尺度模型中围岩、初期支护和二次衬砌的本构模型和物理力学参数,将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔,对三维真实尺度模型进行求解时。
优选的,
利用ANSYS/LS-DYNA、AUTODYNA和ABAQUS中任一个软件建立二维平面应变模型;
利用ANSYS/LS-DYNA、AUTODYNA或ABAQUS中任一个软件建立三维真实尺度模型。
优选的,所述二维平面应变模型为边长5m的正方形模型,对二维平面应变模型划分得到多个尺寸在10mm以内的网格。
优选的,
依据现场情况对二维平面应变模型施加非反射边界条件;
依据现场隧道实际埋深情况对三维真实尺度模型施加边界条件。
优选的,所述研究区域为包含隧道且大于等于隧道直径3倍的区域。
本发明提供的基于二维-三维协同的***动力响应模拟方法具有以下有益效果:1、本发明能求解得到较准确的孔壁***荷载压力时程曲线。2、本发明提供一种简单可行的大型复杂模型***动力响应模拟方法,能够模拟分析大型地下洞室或近接隧道***开挖过程中围岩或二衬的动力响应特性,符合现场工程情况,更具有研究价值和参考意义。
附图说明
为了更清楚地说明本发明实施例及其设计方案,下面将对本实施例所需的附图作简单地介绍。下面描述中的附图仅仅是本发明的部分实施例,对于本领域普通技术人员来说,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1的基于二维-三维协同的***动力响应模拟方法的流程图;
图2为本发明实施例1提供的二维平面应变模型图;
图3为本发明实施例1提供的***荷载压力时程曲线图;
图4为本发明实施例1提供的三维真实尺度模型;
图5为本发明实施例1提供的真实模型***荷载加载示意图;
图6为本发明实施例1提供的数值模拟与实测振动波形对比图。
具体实施方式
为了使本领域技术人员更好的理解本发明的技术方案并能予以实施,下面结合附图和具体实施例对本发明进行详细说明。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
实施例1
参阅图1,一种基于二维-三维协同的***动力响应模拟方法,包括:
基于现场各段别炮孔***参数和装药结构,建立二维平面应变模型;设置二维平面应变模型的边界条件(依据现场情况对二维平面应变模型施加非反射边界条件);设置二维平面应变模型中***和空气参数及状态方程、围岩的本构模型和物理力学参数,在本实施例中,围岩的本构模型选用弹塑性本构模型,初期支护和二次衬砌厚度分别为25cm和50cm,并选用RHT材料模型。参阅图2,图2为本发明建立的二维平面应变模型;对二维平面应变模型进行求解时,直至各段别炮孔的***荷载压力曲线的荷载峰值符合理论和经验公式且三维数值模型二次衬砌的振动波形与现场监测振动数据吻合,否则重复设置二维平面应变模型中***和空气参数及状态方程、围岩的本构模型和物理力学参数,对二维平面应变模型进行求解。
参阅图4,基于圣维南原理利用实际工程形状、尺寸和研究区域围岩类别,建立三维真实尺度模型;设置三维真实尺度模型的边界条件(依据现场隧道实际埋深情况对三维真实尺度模型施加边界条件,除已开挖处其余部分均设置为非反射边界);设置三维真实尺度模型中围岩、初期支护和二次衬砌的本构模型和物理力学参数;参阅图5,将各段别炮孔的***压力荷载曲线分别施加至对应的三维数值模型各段别炮孔,某一炮孔中,考虑爆轰传播的时间效应,P(t)表示随时间变化的炮孔***荷载压力,t=L/D表示距起爆点距离为L的***起爆滞后时间。对三维真实尺度模型进行求解时,将二维平面应变模型中提取出的***荷载时程曲线时间效应依次施加于炮孔孔壁上,求解施加***荷载压力曲线的三维真实尺度模型得到三维真实尺度模型二次衬砌的振动波形,直至三维真实尺度模型二次衬砌的振动波形与现场监测振动数据吻合,否则重新设置三维真实尺度模型中围岩、初期支护和二次衬砌的本构模型和物理力学参数,将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔,求解施加***荷载压力曲线的三维真实尺度模型得到三维真实尺度模型二次衬砌的振动波形,利用二次衬砌的振动波形验证三维真实尺度模型及选取参数的可靠性,根据三维真实尺度模型分析***动力响应。
当炮孔的***荷载压力曲线的荷载峰值与理论计算值的相对误差小于10%,判定炮孔的***荷载压力曲线的荷载峰值符合理论和经验公式。其中,耦合装药时,炮孔的***荷载压力理论计算公式为:
不耦合装药时,炮孔的***荷载压力理论计算公式为:
在本实施例中,现场炮孔孔径为50mm,药径42mm,为不耦合装药结构,计算得孔壁压力峰值为0.978GPa,二维平面应变计算结果参考图3,得孔壁压力峰值为0.96GPa,误差为1.84%,满足要求。
具体的,将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔的步骤包括:根据现场******段别分别建立part,将不同段别的***荷载压力时程曲线按照起爆时间依次施加至对应的part。
在本实施例中,利用ANSYS/LS-DYNA、AUTODYNA和ABAQUS中任一个软件建立二维平面应变模型;利用ANSYS/LS-DYNA、AUTODYNA或ABAQUS中任一个软件建立三维真实尺度模型。同时,利用LS-DYNA软件求解器对模型进行求解。
同时,二维平面应变模型为边长5m的正方形模型,对二维平面应变模型精细划分得到多个尺寸在10mm以内的网格。研究区域为包含隧道且大于等于隧道直径3倍的区域,以选取掏槽段第一段炮孔壁设置为独立部分。
为保证模型建立与所选参数的可靠性,依据现场实测结果,对模型及参数进行了反演调试,图6为某测点处数值计算波形与实测振动波形的对比验证。如图6所示,数值模拟与现场实测波形的幅值、频率及持续时间基本吻合。
以上实施例仅为本发明较佳的具体实施方式,本发明的保护范围不限于此,任何熟悉本领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换,均属于本发明的保护范围。
Claims (8)
1.一种基于二维-三维协同的***动力响应模拟方法,其特征在于,包括:
基于现场各段别炮孔***参数和装药结构,建立二维平面应变模型;
设置二维平面应变模型的边界条件;设置二维平面应变模型中***和空气参数及状态方程、围岩的本构模型和物理力学参数;
对二维平面应变模型进行求解,得到各段别炮孔的***压力荷载曲线;
基于圣维南原理利用实际工程形状、尺寸和研究区域围岩类别,建立三维真实尺度模型;
设置三维真实尺度模型的边界条件;设置三维真实尺度模型中围岩、初期支护和二次衬砌的本构模型和物理力学参数;
将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔;
求解施加***荷载压力曲线的三维真实尺度模型得到三维真实尺度模型二次衬砌的振动波形,利用二次衬砌的振动波形验证三维真实尺度模型及选取参数的可靠性,根据三维真实尺度模型分析***动力响应;
当炮孔的***荷载压力曲线的荷载峰值与理论计算值的相对误差小于10%时,判定炮孔的***荷载压力曲线的荷载峰值符合理论和经验公式;
其中,耦合装药时,炮孔的***荷载压力理论计算公式为:
不耦合装药时,炮孔的***荷载压力理论计算公式为:
2.根据权利要求1所述的基于二维-三维协同的***动力响应模拟方法,其特征在于,还包括:
对二维平面应变模型进行求解时,直至各段别炮孔的***荷载压力曲线的荷载峰值符合理论和经验公式、且三维真实尺度模型二次衬砌的振动波形与现场监测振动数据吻合,得到各段别炮孔的***压力荷载曲线;否则重复设置二维平面应变模型中***和空气参数及状态方程、围岩的本构模型和物理力学参数,对二维平面应变模型进行求解。
3.根据权利要求1所述的基于二维-三维协同的***动力响应模拟方法,其特征在于,所述将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔的步骤包括:
根据现场******段别分别建立part,将不同段别的***荷载压力时程曲线按照起爆时间依次施加至对应的part。
4.根据权利要求1所述的基于二维-三维协同的***动力响应模拟方法,其特征在于,还包括:
对三维真实尺度模型进行求解时,得到三维数值模型二次衬砌的振动波形,直至三维真实尺度模型二次衬砌的振动波形与现场监测振动数据吻合;否则重新设置三维真实尺度模型中围岩、初期支护和二次衬砌的本构模型和物理力学参数,将各段别炮孔的***荷载压力曲线分别施加至对应的三维真实尺度模型各段别炮孔,对三维真实尺度模型进行求解。
5.根据权利要求1所述的基于二维-三维协同的***动力响应模拟方法,其特征在于,
利用ANSYS/LS-DYNA、AUTODYNA和ABAQUS中任一个软件建立二维平面应变模型;
利用ANSYS/LS-DYNA、AUTODYNA或ABAQUS中任一个软件建立三维真实尺度模型。
6.根据权利要求1所述的基于二维-三维协同的***动力响应模拟方法,其特征在于,所述二维平面应变模型为边长5m的正方形模型,对二维平面应变模型划分得到多个尺寸在10mm以内的网格。
7.根据权利要求1所述的基于二维-三维协同的***动力响应模拟方法,其特征在于,
依据现场情况对二维平面应变模型施加非反射边界条件;
依据现场隧道实际埋深情况对三维真实尺度模型施加边界条件。
8.根据权利要求1所述的基于二维-三维协同的***动力响应模拟方法,其特征在于,所述研究区域为包含隧道且大于等于隧道直径3倍的区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210404710.5A CN114722614B (zh) | 2022-04-18 | 2022-04-18 | 一种基于二维-三维协同的***动力响应模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210404710.5A CN114722614B (zh) | 2022-04-18 | 2022-04-18 | 一种基于二维-三维协同的***动力响应模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114722614A CN114722614A (zh) | 2022-07-08 |
CN114722614B true CN114722614B (zh) | 2022-11-25 |
Family
ID=82243677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210404710.5A Active CN114722614B (zh) | 2022-04-18 | 2022-04-18 | 一种基于二维-三维协同的***动力响应模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114722614B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111188622B (zh) * | 2020-03-18 | 2021-07-06 | 山东大学(齐河)新材料与智能装备研究院 | 一种地下工程围岩危石垮塌模拟用平台、模拟***及方法 |
CN115859714B (zh) * | 2022-11-23 | 2023-06-23 | 长安大学 | 一种基于fem-dem联合仿真的岩石***全过程模拟方法 |
CN116796453B (zh) * | 2023-04-18 | 2024-03-08 | 江汉大学 | 一种减振孔设计方法、装置、设备及可读存储介质 |
CN116822311B (zh) * | 2023-08-30 | 2023-11-07 | 昆明理工大学 | 一种***瞬时炮孔周围应力预测方法及*** |
CN117313322B (zh) * | 2023-09-01 | 2024-06-28 | 长安大学 | 一种基于图像模拟的***块度分析方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104359776A (zh) * | 2014-11-14 | 2015-02-18 | 四川大学 | 一种***载荷下i型裂纹动态断裂韧度测试方法 |
CN109374408A (zh) * | 2018-09-28 | 2019-02-22 | 长安大学 | 一种人工充填节理岩体动力特性试验方法 |
CN109580912A (zh) * | 2018-12-19 | 2019-04-05 | 中国地质大学(武汉) | 一种研究***惯性力作用下断层边坡动力响应特征的方法 |
WO2021122576A1 (en) * | 2019-12-16 | 2021-06-24 | Sandvik Mining And Construction Oy | Mine vehicle safety control |
CN113255179A (zh) * | 2021-04-21 | 2021-08-13 | 中铁十六局集团第三工程有限公司 | 隧道***开挖围岩损伤深度计算方法、装置及存储介质 |
WO2021213280A1 (zh) * | 2020-04-21 | 2021-10-28 | 长安大学 | 一种单洞四车道公路隧道断面参数的设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111188622B (zh) * | 2020-03-18 | 2021-07-06 | 山东大学(齐河)新材料与智能装备研究院 | 一种地下工程围岩危石垮塌模拟用平台、模拟***及方法 |
-
2022
- 2022-04-18 CN CN202210404710.5A patent/CN114722614B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104359776A (zh) * | 2014-11-14 | 2015-02-18 | 四川大学 | 一种***载荷下i型裂纹动态断裂韧度测试方法 |
CN109374408A (zh) * | 2018-09-28 | 2019-02-22 | 长安大学 | 一种人工充填节理岩体动力特性试验方法 |
CN109580912A (zh) * | 2018-12-19 | 2019-04-05 | 中国地质大学(武汉) | 一种研究***惯性力作用下断层边坡动力响应特征的方法 |
WO2021122576A1 (en) * | 2019-12-16 | 2021-06-24 | Sandvik Mining And Construction Oy | Mine vehicle safety control |
WO2021213280A1 (zh) * | 2020-04-21 | 2021-10-28 | 长安大学 | 一种单洞四车道公路隧道断面参数的设计方法 |
CN113255179A (zh) * | 2021-04-21 | 2021-08-13 | 中铁十六局集团第三工程有限公司 | 隧道***开挖围岩损伤深度计算方法、装置及存储介质 |
Non-Patent Citations (4)
Title |
---|
Numerical modelling of dynamic response of underground openings under blasting based on field tests;Ulf Nyberg 等;《Ground Support 2016》;20160930;第1-12页 * |
炮孔水耦合装药***应力分布特性试验研究;宗琦 等;《实验力学》;20060930;第21卷(第3期);第393-398页 * |
炮孔空气不耦合装药***的孔壁冲击压力初探;朱子良;《煤矿***》;20040815(第3期);第5-7页 * |
***荷载作用下岩石边坡动态响应的FLAC~(3D)模拟研究;陈占军 等;《***》;20051225;第22卷(第4期);第8-13页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114722614A (zh) | 2022-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114722614B (zh) | 一种基于二维-三维协同的***动力响应模拟方法 | |
CN110390152B (zh) | 一种模拟巷道围岩裂隙演化的离散元方法 | |
Wang et al. | Propagation characteristics and prediction of blast-induced vibration on closely spaced rock tunnels | |
Silva et al. | Practical assessment of rock damage due to blasting | |
Ainalis et al. | Modelling the source of blasting for the numerical simulation of blast-induced ground vibrations: a review | |
Trivino et al. | Assessment of crack initiation and propagation in rock from explosion-induced stress waves and gas expansion by cross-hole seismometry and FEM–DEM method | |
Lu et al. | An equivalent method for blasting vibration simulation | |
CN106326546B (zh) | 一种预测核电工程岩体***损伤深度的方法 | |
Wu et al. | Investigation of the excavation damaged zone around deep TBM tunnel using a Voronoi-element based explicit numerical manifold method | |
CN105787206A (zh) | 一种公路隧道锚杆支护设计方法 | |
Chen et al. | Experimental and numerical study of granite blocks containing two side flaws and a tunnel-shaped opening | |
CN109975119B (zh) | 一种岩石双轴压缩***设计方法 | |
WO2024021907A1 (zh) | 一种隧道掘进***合理微差时间控制方法 | |
CN105589999B (zh) | 应用于地下工程围岩注浆方案的确定方法 | |
Luo et al. | Cavity 3D modeling and correlative techniques based on cavity monitoring | |
CN114297864B (zh) | 一种受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法 | |
CN104239637A (zh) | 一种离散元爆堆形态模拟方法 | |
Gao et al. | Study on numerical simulation of geometric elements of blasting funnel based on PFC5. 0 | |
CN116894291A (zh) | 基于fem-dem数值计算的不同地质隧洞围岩***振动安全判据构建方法 | |
US11920472B2 (en) | Reasonable millisecond time control method for excavation blasting of tunnel | |
Lin et al. | A data assimilation method for blast load prediction | |
Xie et al. | Large-scale field experiments on blast-induced vibration and crater in sand medium | |
CN115292975A (zh) | 护盾式tbm开挖卡机三维连续-非连续耦合数值模拟方法 | |
Basaligheh et al. | Evaluating the safety and the effect of blast loading on the shotcrete together with lattice girder support of tunnels | |
Yu et al. | Numerical analysis of rock joints in tunnel construction during blasting |
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 |