CN107391820A - 隧洞开挖支护无网格法分析方法 - Google Patents
隧洞开挖支护无网格法分析方法 Download PDFInfo
- Publication number
- CN107391820A CN107391820A CN201710554841.0A CN201710554841A CN107391820A CN 107391820 A CN107391820 A CN 107391820A CN 201710554841 A CN201710554841 A CN 201710554841A CN 107391820 A CN107391820 A CN 107391820A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- mover
- omega
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
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)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明涉及隧洞技术领域,尤其是隧洞开挖支护无网格法分析方法，其方法步骤包括：(1)启动ABQUS的CAE，建立洞室的二维模型；(2)划分三角单元模型；(3)输出模型文件Job‑1.inp(4)启动MATLAB，打开MFVM主程序Main_Program.M文件；(5)读取模型文件Job‑1.inp；(6)生成节点信息；(7)使用MFVM计算围岩初始地应力；(8)计算第一开挖步应力场、位移场。本发明有益效果：无网格法的近似函数没有网格依赖性，减少了因网格畸变而引起的困难，更易与CAD软件接口，适用于求解大型科学与工程问题，无网格法的自适应性很强，无网格法的前处理只要节点位置信息，不用网格信息，容易分析复杂的三维结构，无网格计算的结果是光滑连续的，不必再进行应力光顺化处理。
Description
技术领域
本发明涉及隧洞技术领域,尤其是隧洞开挖支护无网格法分析方法。
背景技术
有限元(FEM)是一种鲁棒而成熟的方法，由于它在处理复杂几何体的各类线性和非线性问题时所表现出的通用性及灵活性，使其在工程领域得到了广泛的应用。目前大多数涉及固体及结构的实际工程问题均可利用商业化的FEM软件包进行求解。
然而，FEM有其数值方法上固有的弱点:
(I)形成FEM网格时的计算成本高；
(II)应力精度低
(III)自适应分析困难。
因此，对于上述问题有必要提出隧洞开挖支护无网格法分析方法。
发明内容
针对上述现有技术中存在的不足，本发明的目的在于提供隧洞开挖支护无网格法分析方法。
隧洞开挖支护无网格法分析方法，其方法步骤包括：(1)启动 ABQUS的CAE，建立洞室的二维模型；(2)划分三角单元模型；(3)输出模型文件Job-1.inp
(4)启动MATLAB，打开MFVM主程序Main_Program.M文件；(5)读取模型文件Job-1.inp；(6)生成节点信息；(7)使用MFVM计算围岩初始地应力；(8)计算第一开挖步应力场、位移场。
优选地，其方法步骤进一步包括：(9)锚杆支护计算，流变分析； (10)计算第二开挖步应力场、位移场；(11)开挖区域增加锚杆支护，流变分析；
(12)以文本格式输出位移、应力计算结果；(13)启动TECPLOT后处理软件，生成围岩的各分量的等值线图。
优选地，步骤(8)进一步的包括开挖边界力换算法和开挖边界应力转换法。
优选地，所述开挖边界力换算法的方法为：位于开挖边界上的节点xI受到开挖边界上的面力p及其子域内的体积力b的作用，可建立杂交型无网格有限体积法的平衡方程
其中开挖边界上的面力可由MLS拟合得到，即
N(x)为计算点x(这里是高斯点)在开挖边界上的外法线向量，Φ(x)为x的应力形函数矩阵，为各节点的初始应力场。将式(2)代入 (1)，得到
则开挖荷载为
优选地，所述开挖边界应力转换法的方法为：对于开挖边界上的节点xI，位于Ωex内的积分子域为其边界位于Ωex内的部分为作用在该点上的外载除了面力载荷p(开挖载荷)外，还有体积力b(自重或渗流力等)，由它们共同产生开挖边界上节点内力因此，可建立无网格有限体积形式的平衡方程，
则开挖载荷为
其中N是位于上的高斯点的外法线向量，是高斯积分点的初始应力场，它可以由MLS或者Shepard拟合得到，即
将上式代入式(6)，得到开挖荷载为
同时，开挖载荷又作用在剩余区域ΩR的开挖边界上，因此对位于该边界上的节点xI可建立无网格有限体积法平衡方程
其中，为xI的积分子域边界位于ΩR内的部分。
由于采用上述技术方案，本发明有益效果：无网格法的近似函数没有网格依赖性，减少了因网格畸变而引起的困难，更易与CAD软件接口，适用于求解大型科学与工程问题，无网格法的自适应性很强，无网格法的前处理只要节点位置信息，不用网格信息，容易分析复杂的三维结构，无网格计算的结果是光滑连续的，不必再进行应力光顺化处理。
附图说明
图1至2是本发明的方法流程图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
如图1并结合图2所示，隧洞开挖支护无网格法分析方法，其方法步骤包括：(1)启动ABQUS的CAE，建立洞室的二维模型；(2)划分三角单元模型；(3)输出模型文件Job-1.inp
(4)启动MATLAB，打开MFVM主程序Main_Program.M文件；(5)读取模型文件Job-1.inp；(6)生成节点信息；(7)使用MFVM计算围岩初始地应力；(8)计算第一开挖步应力场、位移场。
其方法步骤进一步包括：(9)锚杆支护计算，流变分析；(10) 计算第二开挖步应力场、位移场；(11)开挖区域增加锚杆支护，流变分析；
(12)以文本格式输出位移、应力计算结果；(13)启动TECPLOT后处理软件，生成围岩的各分量的等值线图。
其中，步骤(8)进一步的包括开挖边界力换算法和开挖边界应力转换法。
进一步的，所述开挖边界力换算法的方法为：位于开挖边界上的节点xI受到开挖边界上的面力p及其子域内的体积力b的作用，可建立杂交型无网格有限体积法的平衡方程
其中开挖边界上的面力可由MLS拟合得到，即
N(x)为计算点x(这里是高斯点)在开挖边界上的外法线向量，Φ(x)为x的应力形函数矩阵，为各节点的初始应力场。将式(2)代入 (1)，得到
则开挖荷载为
进一步的，所述开挖边界应力转换法的方法为：对于开挖边界上的节点xI，位于Ωex内的积分子域为其边界位于Ωex内的部分为作用在该点上的外载除了面力载荷p(开挖载荷)外，还有体积力b(自重或渗流力等)，由它们共同产生开挖边界上节点内力因此，可建立无网格有限体积形式的平衡方程，
则开挖载荷为
其中N是位于上的高斯点的外法线向量，是高斯积分点的初始应力场，它可以由MLS或者Shepard拟合得到，即
将上式代入式(6)，得到开挖荷载为
同时，开挖载荷又作用在剩余区域ΩR的开挖边界上，因此对位于该边界上的节点xI可建立无网格有限体积法平衡方程
其中，为xI的积分子域边界位于ΩR内的部分。
无网格法的近似函数没有网格依赖性，减少了因网格畸变而引起的困难，更易与CAD软件接口，适用于求解大型科学与工程问题，无网格法的自适应性很强，无网格法的前处理只要节点位置信息，不用网格信息，容易分析复杂的三维结构，无网格计算的结果是光滑连续的，不必再进行应力光顺化处理。
以上所述仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (5)
1.隧洞开挖支护无网格法分析方法，其特征在于：其方法步骤包括：
(1)启动ABQUS的CAE，建立洞室的二维模型；
(2)划分三角单元模型；
(3)输出模型文件Job-1.inp
(4)启动MATLAB，打开MFVM主程序Main_Program.M文件；
(5)读取模型文件Job-1.inp；
(6)生成节点信息；
(7)使用MFVM计算围岩初始地应力；
(8)计算第一开挖步应力场、位移场。
2.根据权利要求1所述的隧洞开挖支护无网格法分析方法，其特征在于：其方法步骤进一步包括：
(9)锚杆支护计算，流变分析；
(10)计算第二开挖步应力场、位移场；
(11)开挖区域增加锚杆支护，流变分析；
(12)以文本格式输出位移、应力计算结果；
(13)启动TECPLOT后处理软件，生成围岩的各分量的等值线图。
3.根据权利要求1所述的隧洞开挖支护无网格法分析方法，其特征在于：步骤(8)进一步的包括开挖边界力换算法和开挖边界应力转换法。
4.根据权利要求3所述的隧洞开挖支护无网格法分析方法，其特征在于：所述开挖边界力换算法的方法为：
位于开挖边界上的节点xI受到开挖边界上的面力p及其子域内的体积力b的作用，可建立杂交型无网格有限体积法的平衡方程
<mrow>
<mo>(</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>R</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>N</mi>
<mi>&Phi;</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>)</mo>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
</mrow>
</msub>
<mi>p</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>s</mi>
<mi>R</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>b</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中开挖边界上的面力可由MLS拟合得到，即
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>0</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
N(x)为计算点x(这里是高斯点)在开挖边界上的外法线向量，Φ(x)为x的应力形函数矩阵，为各节点的初始应力场。将式(2)代入(1)，得到
<mrow>
<mo>(</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>R</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>N</mi>
<mi>&Phi;</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>)</mo>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<mo>(</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
</mrow>
</msub>
<mi>N</mi>
<mi>&Phi;</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>)</mo>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>s</mi>
<mi>R</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>b</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
则开挖荷载为
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
</mrow>
</msub>
<mi>N</mi>
<mi>&Phi;</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>s</mi>
<mi>R</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>b</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
5.根据权利要求3所述的隧洞开挖支护无网格法分析方法，其特征在于：所述开挖边界应力转换法的方法为：
对于开挖边界上的节点xI，位于Ωex内的积分子域为其边界位于Ωex内的部分为作用在该点上的外载除了面力载荷p(开挖载荷)外，还有体积力b(自重或渗流力等)，由它们共同产生开挖边界上节点内力因此，可建立无网格有限体积形式的平衡方程，
<mrow>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>E</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>N</mi>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>~</mo>
</mover>
<mn>0</mn>
</msub>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
</mrow>
</msub>
<mi>p</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>s</mi>
<mi>E</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>b</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
则开挖载荷为
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>=</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
</mrow>
</msub>
<mi>p</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>E</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>N</mi>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>~</mo>
</mover>
<mn>0</mn>
</msub>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>s</mi>
<mi>E</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>b</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中N是位于上的高斯点的外法线向量，是高斯积分点的初始应力场，它可以由MLS或者Shepard拟合得到，即
<mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>~</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
将上式代入式(6)，得到开挖荷载为
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>E</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>N</mi>
<mi>&Phi;</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>s</mi>
<mi>E</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>b</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
同时，开挖载荷又作用在剩余区域ΩR的开挖边界上，因此对位于该边界上的节点xI可建立无网格有限体积法平衡方程
<mrow>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>R</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>N</mi>
<mi>&sigma;</mi>
<mi>d</mi>
<mi>&Gamma;</mi>
<mo>=</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<msubsup>
<mi>&Omega;</mi>
<mrow>
<mi>s</mi>
<mi>R</mi>
</mrow>
<mi>I</mi>
</msubsup>
</msub>
<mi>b</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，为xI的积分子域边界位于ΩR内的部分。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710554841.0A CN107391820A (zh) | 2017-07-10 | 2017-07-10 | 隧洞开挖支护无网格法分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710554841.0A CN107391820A (zh) | 2017-07-10 | 2017-07-10 | 隧洞开挖支护无网格法分析方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107391820A true CN107391820A (zh) | 2017-11-24 |
Family
ID=60335537
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710554841.0A Pending CN107391820A (zh) | 2017-07-10 | 2017-07-10 | 隧洞开挖支护无网格法分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107391820A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108344535A (zh) * | 2018-04-19 | 2018-07-31 | 四川大学 | 考虑钻井液压力的水平有效地应力测试方法及测试装置 |
| CN108984908A (zh) * | 2018-07-19 | 2018-12-11 | 中国建筑第八工程局有限公司 | 利用bim技术辅助人防洞室改造的施工方法 |
| CN109165404A (zh) * | 2018-07-16 | 2019-01-08 | 四川大学 | 一种扩展光滑无网格伽辽金法 |
| CN111898187A (zh) * | 2020-07-23 | 2020-11-06 | 武汉大学 | 隧洞开挖渗流模拟分析的纵向模型范围取值的确定方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101950318A (zh) * | 2010-10-02 | 2011-01-19 | 上海交通大学 | 基于混合模型的输水隧道模拟方法 |
| CN102982247A (zh) * | 2012-12-05 | 2013-03-20 | 浙江科技学院 | 一种任意开挖断面隧洞的围岩应力分析方法 |
| CN104021254A (zh) * | 2014-06-20 | 2014-09-03 | 暨南大学 | 盾构隧道开挖后应力分布获取方法 |
| CN106682318A (zh) * | 2016-12-29 | 2017-05-17 | 武汉大学 | 一种基于隧洞围岩纵向变形曲线的初次支护时机计算方法 |
-
2017
- 2017-07-10 CN CN201710554841.0A patent/CN107391820A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101950318A (zh) * | 2010-10-02 | 2011-01-19 | 上海交通大学 | 基于混合模型的输水隧道模拟方法 |
| CN102982247A (zh) * | 2012-12-05 | 2013-03-20 | 浙江科技学院 | 一种任意开挖断面隧洞的围岩应力分析方法 |
| CN104021254A (zh) * | 2014-06-20 | 2014-09-03 | 暨南大学 | 盾构隧道开挖后应力分布获取方法 |
| CN106682318A (zh) * | 2016-12-29 | 2017-05-17 | 武汉大学 | 一种基于隧洞围岩纵向变形曲线的初次支护时机计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 黄哲聪: "杂交型无网格法研究与应用", 《中国博士学位论文全文数据库》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108344535A (zh) * | 2018-04-19 | 2018-07-31 | 四川大学 | 考虑钻井液压力的水平有效地应力测试方法及测试装置 |
| CN109165404A (zh) * | 2018-07-16 | 2019-01-08 | 四川大学 | 一种扩展光滑无网格伽辽金法 |
| CN109165404B (zh) * | 2018-07-16 | 2021-04-23 | 四川大学 | 一种扩展光滑无网格伽辽金法 |
| CN108984908A (zh) * | 2018-07-19 | 2018-12-11 | 中国建筑第八工程局有限公司 | 利用bim技术辅助人防洞室改造的施工方法 |
| CN108984908B (zh) * | 2018-07-19 | 2021-09-21 | 中国建筑第八工程局有限公司 | 利用bim技术辅助人防洞室改造的施工方法 |
| CN111898187A (zh) * | 2020-07-23 | 2020-11-06 | 武汉大学 | 隧洞开挖渗流模拟分析的纵向模型范围取值的确定方法 |
| CN111898187B (zh) * | 2020-07-23 | 2022-10-11 | 武汉大学 | 隧洞开挖渗流模拟分析的纵向模型范围取值的确定方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107391820A (zh) | 隧洞开挖支护无网格法分析方法 | |
| Khoei et al. | A mesh-independent finite element formulation for modeling crack growth in saturated porous media based on an enriched-FEM technique | |
| CN104021254B (zh) | 盾构隧道开挖后应力分布获取方法 | |
| CN104598296A (zh) | 软土深基坑变形数值模拟分析方法 | |
| CN113919197B (zh) | 一种模拟非均质含水层中地下水流的新型三层网格多尺度有限元法 | |
| CN104123400A (zh) | 基于力法的Global-Local细节有限元分析方法 | |
| CN111681316A (zh) | 一种高精度河道地形插值方法 | |
| CN108763625A (zh) | 一种flac3d软件前处理的方法及系统 | |
| Cai et al. | A mixed cover meshless method for elasticity and fracture problems | |
| CN112861409A (zh) | 单桩基础承载能力计算方法、系统、存储介质及设备 | |
| CN106780723B (zh) | 一种利用bim技术实现基坑快速参数化三维建模的方法 | |
| CN113011072B (zh) | 基于midas-pfc3d的离散元复杂模型识别方法 | |
| CN112100890B (zh) | 一种压裂水平缝起裂扩展全三维数学模型的计算方法 | |
| CN110377990B (zh) | 重力作用下含空洞地层浅埋隧道围岩应力及位移求解方法 | |
| CN110704916B (zh) | 基于bim技术的大型复杂三维地质模型网格粗化方法 | |
| CN114970159A (zh) | 一种采用三阶剪切变形板理论的地连墙侧向变形计算方法 | |
| CN115146402A (zh) | 高温合金的氧化过程模拟方法及装置 | |
| CN112580118A (zh) | 批量套绘多测次多个断面地形的方法 | |
| CN106022944B (zh) | 油藏数值模拟中常规与流线模型的数据体转换方法及装置 | |
| Mao et al. | Modeling method of FLAC3D based on RHINO-KUBRIX and deformation analysis of tunnel lining structure | |
| Guo et al. | An elastic solution for twin circular tunnels’ stress in hydrostatic stress field | |
| Rui | Superconvergence of a mixed covolume method for elliptic problems | |
| CN114036785B (zh) | 一种基于洋壳速度结构的深海热液区深部构造应力场模拟方法 | |
| Shen et al. | Parametric Three-Dimensional Modeling of Underground Cable Channels on Flac3D/Python | |
| CN109241209B (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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171124 |
|
| WD01 | Invention patent application deemed withdrawn after publication |