flac3d模拟隧道开挖

new

;网格建立

;;;;;;;;;;;;;;;;;;;;;;;;;

gen zone brick p0 90 0 -30 p1 202 0 -30 p2 90 4 -30 p3 90 0 0 size 112 4 30 ratio 1 1 1 gen zone brick p0 90 0 -30 p1 90 0 0 p2 90 4 -30 p3 0 0 -30 size 30 4 25 ratio 1 1 1.1 gen zone brick p0 90 0 -30 p1 0 0 -30 p2 90 4 -30 p3 90 0 -75 size 25 4 18 ratio 1.1 1 1.1 gen zone brick p0 90 0 -30 p1 90 0 -75 p2 90 4 -30 p3 202 0 -30 size 18 4 112 ratio 1.1 1 1 gen zone brick p0 202 0 -30 p1 292 0 -30 p2 202 4 -30 p3 202 0 0 size 25 4 30 ratio 1.1 1 1 gen zone brick p0 202 0 -30 p1 202 0 -75

FLAC 3D基础知识介绍

一、概述

FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。因此，大大发护展了计算规模。FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。

FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。 FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动

1. FLAC3D的基本知识介绍

http://www.simwe.com/forum/viewthread.php?tid=209662

岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的差分法。差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。 FLAC3D（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存（64K），所以，程序求解的最大结点数仅限于2000个以内。1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。因此，大大发护展了计算规模。FLAC3D是一个三维有限差分程序，目前已发展到V2。1版本。

FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入

《地质工程数值法》 GeoStudio上机实习作业

内容：用GeoStudio的Seep/w模块模拟滑坡不同工况下的地下水渗流曲线。

要求：提供模拟结果文件，现场验收操作。 相关资料：

? 滑坡工程地质剖面图； ? 岩土参数； ? 渗流模拟工况。

1 计算模型的建立

根据2#滑坡的工程地质条件，选择该滑坡钻探剖面（图1）作为模拟剖面。其方位NW113°，滑坡前缘高程1190m，后缘高程1398m，前缘和后部较薄，中部较厚，根据钻孔揭露，滑体物质主要是灰黑色碎块石土，滑带土主要为灰黒色含角砾粉质粘土、角砾质土或含角砾泥炭质土，滑带厚度0.6～1.5m，滑床为炭质页岩和炭质灰岩的碎块石土。基于以上条件，对该堆积层滑坡的坡体地质模型进行网格剖分，共剖分为2829个网格单元，如图2所示。

1450高程(m)1400135013001250Qcol4第四系全新统崩积层碎块石土滑动面及滑动方向113°比例尺 1:1000高程(m)145014001350ZK31322.94推测潜在滑面26.6(1296.34)33.8(1289.14)碎块石土ZK41286.011300西河ZK51237.00潜在滑体碎块石土27.8(1258.21)41.71250

《地质工程数值法》 GeoStudio上机实习作业

内容：用GeoStudio的Seep/w模块模拟滑坡不同工况下的地下水渗流曲线。

要求：提供模拟结果文件，现场验收操作。 相关资料：

? 滑坡工程地质剖面图； ? 岩土参数； ? 渗流模拟工况。

1 计算模型的建立

根据2#滑坡的工程地质条件，选择该滑坡钻探剖面（图1）作为模拟剖面。其方位NW113°，滑坡前缘高程1190m，后缘高程1398m，前缘和后部较薄，中部较厚，根据钻孔揭露，滑体物质主要是灰黑色碎块石土，滑带土主要为灰黒色含角砾粉质粘土、角砾质土或含角砾泥炭质土，滑带厚度0.6～1.5m，滑床为炭质页岩和炭质灰岩的碎块石土。基于以上条件，对该堆积层滑坡的坡体地质模型进行网格剖分，共剖分为2829个网格单元，如图2所示。

1450高程(m)1400135013001250Qcol4第四系全新统崩积层碎块石土滑动面及滑动方向113°比例尺 1:1000高程(m)145014001350ZK31322.94推测潜在滑面26.6(1296.34)33.8(1289.14)碎块石土ZK41286.011300西河ZK51237.00潜在滑体碎块石土27.8(1258.21)41.71250

FLAC3D使用教程

总第115期

　　　　　　　　西部探矿工程seriesNo.115

2005年第11期WEST-CHINAEXPLORATIONENGINEERINGNov.2005

文章编号:1004—5716(2005)11—0001—04

中图分类号:TV55114+2　文献标识码:B

岩土工程

FLAC3D及其在地下空间开挖分析中的一些问题

尹清锋,白　冰

(北京交通大学土木建筑工程学院,北京100044)

摘　要:首先阐述了FLAC3D的基本原理、主要特点和适用范围,,给出了在地下空间开挖分析中处理特殊问题的经验,关键词:FLAC3D;显式有限差分法;大变形;地下空间

1　概述

FLAC3D是FastLagrangianin32mensions的简称,—四面体的的表面所围成的空间闭区域;∑———四面体的外表面;

nj———外表面的单位法向向量分量,i,j=1～3。

维显式有限差分程序AC应用软件的拓展,可以模拟土质、,可以精确的模拟屈服、塑性流动、软化直至破坏的整个过程,尤其适用于软弱介质材料的弹塑性分析、大变形分析以及施工过程模拟[1～3],并且可以在初始模型中加入诸如断裂、节理构造等地质因素[4]。

FLAC

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

本科生毕业论文（设计）

题 目 ：基于flac3D深基坑开挖模拟与支护设计

指导教师: 职称: 评 阅 人: 职称:

摘 要

随着城市化过程中不断涌现的高层建筑和超高层建筑以及城市地下空间的开发，深基坑工程越来越多，深基坑工程项目的规模和复杂性日益增大，给深基坑工程的设计和施工带来了更大的挑战。在这样的背景下，深基坑支护结构设计和变形量预测已成为岩土工程领域的重要研究课题之一。本文以武汉市万达广场深基坑工程作为研究对象，利用勘查资料和深基坑支护结构设计要求，比选合理的基坑支护方案并进行相应的计算设计。同时，本文针对深基坑工程变形量验算等难以解决的问题引用了flac3D数值模拟方法，对基坑开挖、支护结构施工进行全方位的模拟监测，将计算设计结果和模拟计算结果进行对比验算，得出比较合理的支护结构设计方案和变形量控制方案。

根据基坑实际情况和勘查资料，本文选择的围护方案为以大直径混凝土排桩、双排桩、角撑与对顶撑相结合的内支撑为主的多种联合支护方案，结合坡顶大面积卸土减载、坑内被动区加固的措施。计算部分主要设计计算大直径混凝土排桩（钻孔灌注桩）桩长、内力和配筋，而对