flac3d一个多层土的例子
FLAC3D模拟实例循环开挖与支护
FLAC3D模拟实例循环开挖与支护nres ini.savset geometry=0.001ini ydis0ini xdis0ini zdis0ini yvel0ini xvel0ini zvel0m mprop bulk 4.0e9shear 2.5e9fri32coh 2.0e6& range grou diban-shayan;prop bulk 1.8e9shear 1.2e9fri25coh 1.0e6& range grou diban-niyan any grou hangdao any;prop bulk 1.2e9shear0.8e9fri22coh0.8e6& range grou diban-gentuyan;prop bulk 1.9e9shear 1.3e9fri24coh 1.0e6& range grou diban-tniyan;prop bulk0.7e9shear0.8e9fri21coh0.7e6& range grou mc any grou gzm any;prop bulk 3.0e9shear 2.5e9fri30coh 1.8e6& range grou dingban-fenshayan;prop bulk 1.5e9shear 1.2e9fri25coh 1.1e6& range grou dingban-niyan;prop bulk 3.5e9shear 2.5e9fri34coh 1.4e6& range grou dingban-shayan;添加接触面gen separate gzminterface1wrap mc gzminterface2wrap dingban-fenshayan gzm interface1prop kn20e9ks10e9tens1e9 interface2prop kn20e9ks10e9tens1e9set mech ratio=5e-4def excav_mcloop n(excav_p,excav_p_z+cut_liang);每次开挖量cut_0=excav_pcut_1=excav_p+cut_liang;开挖commandm null range grou gzm z cut_0cut_1step100end_commandn=excav_p+cut_Liangexcav_p=excav_p+cut_Liang;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
FLAC3D上机题目AND实例分析
FLAC 3D
数值模拟上机题
计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为：
容重r ＝2500 kg/m 3, 弹性模量E ＝1×108 Pa ，泊松比μ＝0.3， 抗拉强度σt ＝0.8×106 Pa ，内聚力C ＝4.2×104 Pa ，摩擦角φ＝17°
试用FLAC 3D 软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC 3D 软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。
附 换算公式： 1 kN/m 3= 100 kg/m 3 剪切弹性模量：)
1(2μ+=
E
G
体积弹性模量：)
21(3μ-=
E
K
图1 倾角为30°的边坡(单位：m)
图2 倾角为45°的边坡(单位：m)
图3 倾角为60°的边坡(单位：m)
图4 边坡开挖算例分析 eg.
以45°边坡为例：网格剖分图如下：
实习报告要求：
1、不同坡角边坡位移场和应力场特征；
2、边坡开挖前后位移场和应力场特征；
3、运用强度折减法求解边坡稳定性系数。
实例分析：
eg.
以45°边坡为例，网格剖分图如下：
图1 网格剖分图
图2 速度矢量图
图3 速度等值线图
图4 位移等值线图
最终计算边坡的稳定性系数为：Fs＝1.15625。
flac3d5.0结构单元教程
history add id=10 sel node ydisp id=7
history add id=30 sel beamsel moment mz end2 cid=1 ; moment, right of SEL-1
history add id=31 sel beamsel moment mz end1 cid=2 ; moment, left of SEL-2
• 壳单元(shell)
➢ FLAC3D中的结构单元是岩土工程中实际结 构的一种“抽象”，即采用简单的单元形式 来模拟复杂的结构体。
➢ 结构单元由结构节点(node)和结构构件 (SELs)构成。
➢ 结构单元中的节点(node)可以与周围的实体 网格(zone)或其它结构节点建立连接(link)，
• 土工格栅(geogrid)
; Print out beam responses.
list sel beam force
list sel beam moment
list sel node disp range id=7
return
21
.
4.1、简支梁（beam单元）承受两个相等集中载 荷
Beam_concent_loads_Example 1.3
1. 结构单元的类型 2. 结构单元的建模方法 3. 结构单元的参数取值 4. 结构单元实例分析 5. 关于link
3
.
1、结构单元的类型
FLAC3D中包含六种形式的结构单元，可以分成两类：
➢ 线型结构单元：
• 梁单元(beam) • 锚索单元(cable) • 桩单元(pile)
➢ 壳型结构单元：
8
2、结构单元的建模方法—线型. 结构单 元
FLAC3D快速入门及简单实例
FLAC3D快速⼊门及简单实例FLAC3D快速⼊门及简单实例李佳宇编LJY指南针教程前⾔FLAC及FLAC3D是由国际著名学者、英国皇家⼯程院院⼠、离散元的发明⼈Peter Cundall博⼠在70年代中期开始研究的，主要⾯对岩⼟⼯程的通⽤软件系统，⽬前已经在全球70多个国家得到⼴泛应⽤，在岩⼟⼯程学术界和⼯业界赢得了⼴泛的赞誉。
前国际岩⽯⼒学会主席 C.Fairhurst(1994)对FLAC程序的评价是：“现在它是国际上⼴泛应⽤的可靠程序。
”我从研⼆(2010年)开始接触FLAC3D，最初的原因是导师要求每⼀个⼈⾄少学会⼀个数值计算软件，⽽他嘴⾥每天念叨最多的就是FLAC，⾃⼰当时对数值计算⼀⽆所知，便答应⽼师要学会FLAC3D。
第⼀次打开软件界⾯，我⼼⾥就凉了⼤半截，⾯对着⼀个操作界⾯跟记事本⽆异的所谓“功能强⼤”的岩⼟⼯程专业软件，半点兴趣也提不起来。
年底，从项⽬⼯地回到学校准备论⽂开题，⽼师对我的开题报告⾮常不满意，当着全教研室师⽣的⾯，劈头盖脸⼤批⼀顿，第⼆天⼜找谈话。
在巨⼤的压⼒和强烈的⾃尊⼼驱使下，我硬着头⽪开始啃FLAC3D，⼀个半⽉之后，终于有了初步的计算结果，对⽼师有个交代，我也能回家过年了。
前⾯这⼀段过程可能是⼤多数FLAC3D初学者的必经阶段，或者是即将开始软件学习的⼈惧怕的事情。
毫⽆疑问，FLAC3D极其不友好的界⾯是阻碍初学者前进的很⼤障碍，当然还包括它是⼀个全英⽂的软件。
但是当你费尽周折的⾛进FLAC3D的世界，你就会发现它独特的魅⼒，⽐如简洁的界⾯，快捷的命令流操作，⾼效的计算⽅法，不易报错等等。
另外⼀个拿不上台⾯的优点就是它⾮常⼩巧，包括Manual在内⼀共才⼏⼗兆⼤⼩，⽽且已经被破解成绿⾊版，只要把它和命令流装进U盘，你就可以随便找⼀个⾝边功能最强⼤的电脑开始计算了，如果你有过ANSYS、ABAQUS等⼤型软件痛苦的安装经历，你便能毕业之后，本以为不⽤再接触数值计算，但⼯作需要使得我⼜⼀次开始与理解“绿⾊版”的含义，当然还请⼤家尊重知识产权，⽀持正版。
flac3d工程实例
地下建筑结构
4.分步开挖支护计算模型
对基坑放坡开挖BC段,可以采取对开挖土体和基坑围岩土体分 别建模,建模中用到基本单元的八节点定义方式。然后用null定义 开挖土体部分的模块来模拟,如图4。
本工程实例中基坑坡率较大,近乎直立开挖,也可直接采取每 步开挖支护中,再多次分步开挖得到阶梯状开挖面的方式近似模拟 坡面开挖。经模拟计算验证,该建模方式对基坑开挖支护数值模拟 计算结果影响很小,可忽略不计。
本工程沉降水准观测,对基坑周边边坡土体和建筑物沉降监测共 布设10+8个沉降观测点。其中,在基坑周边土体共布设10个水平位移 和沉降观测点,东面两栋距基坑较近的被测建筑物的首层柱上共为3个 沉降观测点。其中对周边建筑物进行13次沉降观测。观测周期与工程 进度密切联系,基坑开挖期间,土体扰动对周围环境影响较大,沉降速 率较大,故保证每隔2天一次观测,时刻注意环境动态,后期施工过程中, 根据实际情况相对减小监测频率。
2、边界条件
边界条件分为位移边界条件和受力边界条件利用FLAC3D中fish语言 程序命令:fix(位移方向)range(位移约束面),可在网格模型中设定位移 边界条件。 模型除x=0面及基坑顶面为自由面不设置位移边界条件外,其余面均 采用法向约束。坑底边界静止不动,假设为固定铰支,限制三个方向的位 移;模型x轴及y轴方向两侧避均施加边界约束条件,限制该临空方向的位 移;基坑顶面为设为自由面,不加约束。地下水己进行处理,可不考虑。 利用FLAC3D中fish语言程序命令:apply szz(荷载大小)range(荷载 范围),可在网格模型中定义受力边界条件。本工程CD段基坑周边己有建 筑物荷载,荷载大小为90KN。
图19第四步开挖支护x方向位移等值线图
图20第五步开挖支护x方向位移等值线图
FLAC3D基本原理及简单实例
V 1 2v ( 1 2 3) V E
• 如果地下岩体处于静水应力状态，则以 1 2 3 rH 带入上式可 得
V (1 2v) rH 3rH V E K
FLAC3D基础知识
FLAC 的计算循环图
3D
二．FLAC3D中材料的本构模型
1、FLAC3D内置材料本构模型
8 6 3 3
柱体网格
块体外围渐变放射网格 六面体隧道外围渐变放 射网格 柱形交叉隧道网格 棱椎体网格 退化块体网格 柱形壳体网格
cylinder
radbrick radtunnel cylint pyramid dbrick cshell
6
15 14 14 5 7 10
3
4 4 5 3 3 4
FLAC3D基础知识
f t 3 t
式中， 是摩擦角，C是粘聚力， t 是张拉强度，且有：
N
3
张拉强度不超过 值，最大值由下式给定：
1 sin 1 sin
t max
c tan
2.2 FLAC3D常用材料本构模型
Mohr-Coulomb模型
流动法则
Mohr-Coulomb破坏准则
FLAC 基本操作及简单应用
3D
一 FLAC3D基础知识
二 FLAC3D中材料的本构模型
三 FLAC3D的应用实例
FLAC3D基础知识
FLAC 分析的基本组成部分
3D
FLAC3D基础知识
FLAC 基本形状网格的基本特征
名称
六面块体网格 楔形网格
3D
关键词
brick wedge
控制点个数 单元划分的方向个数
SURFER.DAT
FLAC3D学习资料
c / kPa
/( )
MPa
1500
10
15
8.0
0.33
粘土
1800
20
20
4.0
0.33
二、分析问题的过程
建立网格
初始条件
边界条件 初始应力平衡
前处理
外荷载 求解
4
后处理
二、分析问题的过程
建立网格
FLAC3D使用命令generate zone 生成实 体网格，它实际上是通过调用网格库中的基本 形状网格，对其进行匹配、连接，最终得到用 户期望的几何形状的实体网格模型。
5
二、分析问题的过程
使用generate zone 生成基本形状网格的常用关键词
关键词 add dimension edge fill group P0 ~P16 ratio size
用途 用于以P0为原点的局部坐标系建模 定义内部区域的尺寸 定义网格边长 定义网格内部填充区域 定义某一范围内的网格组名 建立各种形状网格的控制点 定义相邻网格单元尺寸大小比率 定义网格在各坐标方向上的单元数目
7
二、分析问题的过程
边界条件
边界条件的定义通常用到apply、fix、free这三个命令
8
二、分析问题的过程
边界条件
3、如果是压应力，那么赋值的时候应该是负值。
9
二、分析问题的过程
边界条件
10
二、分析问题的过程
初始应力平衡
FLAC3D中，在模型发生变化以前必须使模型先 达到一个初始力的平衡状态。如果是简单的模型，可 以直接定义模型的边界条件和初始条件，模型可能自 然地已经达到了初始的平衡状态。但是，在大多数情 况下，模型需要在给定边界条件和初始条件的前提下， 进行一些计算步的计算使模型达到初始平衡状态，尤 其在模型比较复杂或者在模型中有多种材料的时候， 这时就需要用step或者solve命令。
FLAC3D基础介绍--FLAC3D在岩土工程中的应用-文档资料
与隐式有限元程序相比的常用模型
岩土力学通用模型(边坡稳定性分 析，地下开挖)
破坏后研究(失稳过程，立柱屈服， 顶板崩落)
松散沉积地层中的开挖
层状材料破坏后研究
粘土 岩石
16/74
主要内容
FLAC3D软件简介 FLAC3D的基本原理 FLAC3D的前后处理 流-固耦合分析 接触单元与应用 完全非线性的动力分析 自定义本构模型的基本方法 结构单元及应用
GeoHohai
21/74
前后处理功能的优点
多种zone类型 后处理快捷、方便、丰富 计算过程中的hist变量动态显示 FISH可进行参数化模型设计 单元状态的可编程 计算暂停时的后处理与可保存
GeoHohai
22/74
前后处理功能的缺点
复杂模型的建模功能不强
可以编程导入其他软件形成的网格(比如：Ansys、 Adina、GeoCAD)
PROP biot_c 0 (or INI fmod 0)
GeoHohai
34/74
无渗流计算——孔压的力学响应
不排水短期响应 两种分析方法：干法和湿法
干法：Ku=K+a2M 两种破坏形式
WATER或INI获得常孔压，不排水的c,φ (孔压改变较小) φ=0,c=cu (M>>K+4/3G)
GeoHohai
命令栏
19/74
菜单驱动(Plot)
GeoHohai
20/74
Case-2 一个最简单的例子
gen zon bri size 3 3 3 ;建立网格
model elas
;材料参数
prop bulk 3e8 shear 1e8
flac3d基坑开挖支护例子
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;七区（？？侧帮 喷层？？？）
p4 6,20,10.865 p5 2.1,20,12 p6 6,0,11.02 p7 6,20,11.02 &
size 8 20 1 ratio 1 1 1 group n3
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;四区
;part1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;一区
gen zone radcyl p0 0,0,13.4 p1 6,0,12.41 p2 0,20,13.4 p3 0,0,16.7 &
p4 6,20,12.41 p5 0,20,16.7 p6 6,0,15.13 p7 6,20,15.13 &
p4 15,20,7.36 p5 6,20,11.02 p6 15,0,8.6 p7 15,20,8.6 &
size 12 20 8 ratio 1.05 1 1 group n3
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十一区
p4 6,20,9.935 p5 1.5,20,12 p6 6,0,10.865 p7 6,20,10.865 &
size 8 20 6 ratio 1 1 1 group n3
FLAC3D实例分析教程
03
FLAC3D求解过程
求解器选择及参数的求解器，如稀疏矩阵求解器
或迭代求解器。
设置求解参数
02
根据问题特点，设置合适的求解参数，如收敛准则、迭代次数、
时间步长等。
定义材料属性
03
输入材料的物理和力学参数，如密度、弹性模量、泊松比等。
求解过程监控与调整
材料属性赋值
将定义好的材料属性赋值 给几何模型中的相应区域。
边界条件与初始条件设置
定义边界条件
根据实际问题定义模型的边界条件，如固定边界、 自由边界、位移边界等。
设置初始条件
设置模型的初始应力、初始位移等初始条件。
施加荷载
根据实际问题施加荷载，如重力、压力、温度荷 载等。
网格划分与优化
网格类型选择
操作界面及基本设置
操作界面
FLAC3D的操作界面包括菜单栏、工具栏、模型视图、属 性视图和输出窗口等部分，方便用户进行建模、分析和后 处理等操作。
基本设置 在使用FLAC3D前，需要进行一些基本设置，如选择单位 制、设置材料参数、定义边界条件等。这些设置将直接影 响模拟结果的准确性和可靠性。
建模流程
地震工程
FLAC3D可用于模拟地震波在岩土体中的传 播过程，分析地震对建筑物和岩土工程设 施的影响。
地下工程
软件可模拟地下洞室、巷道、地铁等地下 工程的开挖和支护过程，分析围岩稳定性 和支护结构受力情况。
基坑工程
软件可用于分析基坑开挖过程中的土压力 分布、支护结构变形和稳定性等问题。
隧道工程
FLAC3D可应用于隧道工程的开挖、支护和 衬砌设计，评估隧道施工对周围岩体的影 响。
FLAC3D的建模流程包括建立几何模型、划分网格、定义 材料属性、施加边界条件和初始条件等步骤。用户可根据 实际需求选择合适的建模方法和工具。
Flac3d实例分析教程
建筑面积计算规则一、计算建筑面积的范围1.单层建筑物不论其高度如何，均按一层计算建筑面积。
其建筑面积按建筑物外墙勒脚以上结构的外围水平面积计算。
单层建筑物内设有部分楼层者，首层建筑面积已包括在单层建筑物内，二层及二层以上应计算建筑面积。
高低联跨的单层建筑物，需分别计算建筑面积时，应以结构外边线为界分别计算。
2.多层建筑物建筑面积，按各层建筑面积之和计算，其首层建筑面积按外墙勒脚以上结构的外围水平面积计算，二层及二层以上按外墙结构的外围水平面积计算。
3.同一建筑物如结构、层数不同时，应分别计算建筑面积。
4.地下室、半地下室、地下车间、仓库、商店、车站、地下指挥部等及相应的出入口建筑面积，按其上口外墙（不包括采光井、防潮层及其保护墙）外围水平面积计算。
5.建于坡地的建筑物利用吊脚空间设置架空层和深基础地下架空层设计加以利用时，其层高超过2.2m，按围护结构外围水平面积计算建筑面积。
6.穿过建筑物的通道，建筑物内的门厅、大厅，不论其高度如何均按一层建筑面积计算。
门厅、大厅内设有回廊时，按其自然层的水平投影面积计算建筑面积。
7.室内楼梯间、电梯井、提物井、垃圾道、管道井等均建筑物的自然层计算建筑面积。
8.书库、立体仓库设有结构层的，按结构层计算建筑面积，没有结构层的，按承重书架层或货架层计算建筑面积。
9.有围护结构的舞台灯光控制室，按其围护结构外围水平面积乘以层数计算建筑面积。
10.建筑物内设备管道层、技术层、贮藏室其层高超过2.2m时，应计算建筑面积。
11.有柱的雨蓬、车棚、货棚、站台等、按柱外围水平面积计算建筑面积；独立柱的雨蓬、单排柱的车棚、货棚、站台等，按其顶盖水平投影面积的一半计算建筑面积。
12.屋面上部有围护结构的楼梯间、水箱间、电梯机房等，按围护结构水平面积计算建筑面积。
13.建筑物外有围护结构的门斗、眺望间、观望电梯间、阳台、橱窗、挑廊、走廊等，按其围护结构外围水平面积计算建筑面积。
14.建筑物外有柱和顶盖走廊、檐廊，按柱外围水平面积计算建筑面积；有盖无柱的走廊、檐廊按其顶盖投影面积一半计算建筑面积。
FLAC3D条带开采案例（代码）
FLAC3D条带开采案例（代码）ngen zone brick p0 0 0 0 p1 2500 0 0 p2 0 2000 0 p3 0 0 56 size 100 80 2 group 煤层底板;gen zone brick p0 0 0 56 p1 2500 0 56 p2 0 2000 56 p3 0 0 60 size 100 80 1 group 煤层;gen zone brick p0 0 0 60 p1 2500 0 60 p2 0 2000 60 p3 0 0 76 size 100 80 1 group 直接顶;顶板gen zone brick p0 0 0 76 p1 2500 0 76 p2 0 2000 76 p3 0 0 100 size 100 80 1 group 老顶;顶板gen zone brick p0 0 0 100 p1 2500 0 100 p2 0 2000 100 p3 0 0 280 size 50 40 4 group 中粒、粉砂岩互层;顶板gen zone brick p0 0 0 280 p1 2500 0 280 p2 0 2000 280 p3 0 0 480 size 50 40 4 group 泥岩、砂泥岩互层;顶板gen zone brick p0 0 0 480 p1 2500 0 480 p2 0 2000 480 p3 0 0 640 size 50 40 4 group 表土层;顶板modelmohrprop bulk 6.94e9 shear 5.21e9 coh 24.6e7 fric 38 tens 3e7 range group 煤层底板prop bulk 6.25e9 shear 0.35e9 coh 1e7 fric 25 tens 1e7 range group 煤层prop bulk 1.33e9 shear 0.44e9 coh 1.5e7 fric 32 tens 1e7 range group 直接顶prop bulk 7.44e9 shear 5.12e9 coh 8.56e7 fric 37 tens 3e7 range group 老顶prop bulk 8.26e9 shear 4.72e9 coh 6.25e7 fric 33 tens 2.6e7 range group 中粒、粉砂岩互层prop bulk 8.83e9 shear 1.89e9 coh 1.5e7 fric 33 tens 1e7 range group 泥岩、砂泥岩互层prop bulk 0.017e9 shear 0.0036e9 coh 0.1e7 fric 20 tens 1e7 range group 表土层initial density 2.6e3 range group 煤层底板initial density 1.4e3 range group 煤层initial density 2.4e3 range group 直接顶initial density 2.5e3 range group 老顶initial density 2.64e3 range group 中粒、粉砂岩互层initial density 2.6e3 range group 泥岩、砂泥岩互层initial density 1.8e3 range group 表土层fix z range z -0.1 0.1fix x range x -0.1 0.1fix x range x 2499.9 2500.1fix y range y -0.1 0.1fix y range y 1999.9 2000.1setgrav 0 0 -9.8histunbalplotcontzdissolvesave 试验平衡.savprop bulk 6.94e9 shear 5.21e9 coh 24.6e6 fric 38 tens 3e6 range group 煤层底板prop bulk 6.25e9 shear 0.35e9 coh 1e6 fric 25 tens 0.1e6 range group 煤层prop bulk 1.33e9 shear 0.44e9 coh 1.5e6 fric 32 tens 1e6 range group 直接顶prop bulk 7.44e9 shear 5.12e9 coh 8.56e6 fric 37 tens 3e6 range group 老顶prop bulk 8.26e9 shear 4.72e9 coh 6.25e6 fric 33 tens 2.6e6 range group 中粒、粉砂岩互层prop bulk 8.83e9 shear 1.89e9 coh 1.5e6 fric 33 tens 1e6 range group 泥岩、砂泥岩互层prop bulk 0.017e9 shear 0.0036e9 coh 0.01e6 fric 20 tens 0e6 range group 表土层ini state 0initialxdisp 0 ydisp 0 zdisp 0inixvel 0.0 yvel 0.0 zvel 0.0group zone 开挖range x 850 1650 y 650 780 z 56 60 group zone 开挖range x 850 1650 y 910 1040 z 56 60 group zone 开挖range x 850 1650 y 1170 1300 z 56 60 m n range group 开挖plotcontzdissolvesave 试验.sav。
复杂地质体FLAC_3D_模型快速生成_柴红保 (1)
收稿日期：2009-05-12基金项目：高等学校博士学科点专项科研基金项目(20060533071)；国家自然科学基金项目(50774093)作者简介：柴红保(1977-)，男，河南南阳人，讲师，研究方向：岩土边坡.E-mail:chbmoon@摘要：直接利用FLAC 3D 软件建立复杂三维地质模型比较困难，而使用AutoCAD 软件可以建立复杂三维模型，且方便、快捷和精确.将AutoCAD 软件建立的复杂三维模型导入到ANSYS 软件中，使用ANSYS 软件对模型进行网格划分，使用ANSYS 软件与FLAC 3D 软件的接口程序，将ANSYS 软件已经划分过网格的三维模型的单元数据转换为FLAC 3D 单元的数据，在FLAC 3D 软件中导入三维模型，从而实现在FLAC 3D 软件中快速建立复杂三维地质模型.这为利用FLAC 3D 软件对复杂三维地质模型进行分析提供了一种便利的手段.该方法借助于AutoCAD 软件的优点，使得在FLAC 3D 中精确快速建立复杂地质体三维模型更快捷.用此法建立贵州瓮福磷矿的某一复杂三维边坡FLAC 3D 模型，对该模型进行了稳定性分析.经分析该边坡的安全系数为3.03，边坡处于稳定状态.借助于该法可使利用FLAC 3D 进行复杂三维模型的分析更快捷.关键词：AutoCAD ；ANSYS ；三维模型；FLAC 3D 中图分类号：TU457文献标识码：A文章编号：1674-5876（2009）02-0025-04复杂地质体FLAC 3D 模型快速生成柴红保1，龙雪鸣2(1.煤矿安全开采技术湖南省重点实验室，湖南湘潭411201；2.湖南科技大学生命科学学院，湖南湘潭411201)FLAC 3D 软件被广泛应用于岩土工程问题的数值分析.该软件对内存的要求不是很高，计算速度很快，在解决岩土工程问题上具有许多优越性，已逐渐成为工程技术人员理想的三维数值模拟工具.但它在建模上有些不足之处[1-4]，尤其是对于复杂多介质、多边界的地质体；为解决这个问题，工程人员往往通过一些简化模型来描述岩土模型，而实际上，不同的岩性具有不同的力学性质，过度的简化使数值模拟结果的可靠性受到质疑.另一方面，三维地学模拟表现出较好的三维建模能力，它可以精确地表示出各种不同岩性的岩层、地质结构体在空间上的分布.然而，地学模拟相对独立于数值模拟，其功能一般局限在可视化和定性判断上.因此，如何把三维地质建模方面的优势和数值模拟分析结合起来成为必要.许多学者在这一方面做了有意义的工作[1-5].AutoCAD 软件是当今最流行的计算机制图软件之一，在建筑、机械、矿山等各个领域都被广泛应用，广大科研人员都能够熟练使用该软件进行绘图，该软件在数据充足的情况下可以建立精确的三维模型.AutoCAD 的图形文件可以输出为SAT 格式，该格式的文件能够在ANSYS 软件使用，并可进一步进行单元划分，经过网格划分的模型，通过ANSYS 与FLAC 3D 软件的接口程序，将ANSYS 的模型转化为FLAC 3D 模型.通过这种方式能够在FLAC 3D 中实现复杂三位地质模型快速精确地建立，从而使利用FLAC 3D 进行复杂地形地貌的对象的分析研究更加便利.1A ut oC A D 中三维模型的建立AutoCAD 软件可以绘制各种图形，并且提供了丰富的对图形进行编辑和修改的功能，使用AutoCAD 软件能够建立复杂的三维模型.AutoCAD 软件在三维绘图方面有很丰富的功能，提供很多种可以绘制三维图形的方式.例如，可以将平面拉身为三维实体图形；或直接通过三维坐标点绘制三维实体图形.并且，AutoCAD 软件提供了对三维实体进行各种布尔运算的功能，这将使使用者能够借助于这些功能建立各种矿业工程研究Mineral Engineering Research第24卷第2期2009年6月Vol.24No.2Jun.200925复杂的为模型.对于建立有复杂地形地貌的三维图形，可以通过输入复杂地形地貌的关键坐标，使用绘制三维曲面的命令绘制出复杂地形地貌的曲面图形，在此基础上建立复杂三维地质模型.图1是通过首先建立复杂地形的断面图，再通过使用拉伸面为实体的命令将断面拉伸为三维实体模型.图2是图1中三维模型中的矿体三维模型.该三维模型包含两层矿体，分别为A 、B 矿体.这种建立三维地质模型的方式对于较小距离的模型是非常实用的.对于比较复杂的三维地质模型模型，可以通过首先建立AutoCAD 的脚本文件(SCR 文件)，再在AutoCAD 中执行该脚本文件的方式快速的自动生成三维地质模型.在AutoCAD 软件中，将上面建立的三维地质模型输出为SAT 文件格式，这为在ANSYS 中顺利的倒入模型做好了准备工作.2A N SY S 中模型的导入导出ANSYS 软件提供了导入如IGES 、CIF 和SAT 等多种格式文件的功能，这些文件格式都是国际一些常用图形软件支持的一种通用的文件格式，这些文件格式方便了不同软件的数据交流共享.将前面在AutoCAD 软件中建立的三维地质模型的SAT 格式文件，通过ANSYS 菜单中的FILE-IM PORT-SAT 菜单项的操作导入到ANSYS 中.在ANSYS 软件中导入模型后，定义几种材料类型，此处材料类型可以随便定义，所采用的参数可以随便设置，该步骤只是为了在ANSYS 软件中划分网格的时候能够使模型型中不同性质材料部分具有不同的材料类型，这样在模型转换为FLAC 3D 模型时，方便在FLAC 3D 中区分不同的单元.在ANSYS 软件中定义好要进行网格划分时所有使用的单元类型，这里要划分的是三维地质模型，所以要定义的单元类型也需要时三维实体单元，此处选择的为SOLID45单元.在对模型进行网格划分前，可以对所要划分的网格单元的精度进行控制，可以使用划分网格的工具对要划分的单元的面积、单元的边长等方面进行设置，通过这些设置能够使所划分的网格的精度能够满足我们所要分析的问题的精度要求.如果网格划分的质量不高，在进行计算的时候也会出现由于网格畸变引起的计算上的错误，所以，应该根据所研究对象的具体情况对网格的尺寸进行比较好的控制，以保证计算的精度.在完成上面的步骤后，就可以进行网格划分.图1导入到ANSYS 中划分的网格如图3所示.3FL A C 3D 模型的自动生成由于FLAC 3D 能够使用ANSYS 软件所划分的单元建立模型，通过在ANSYS 软件中运行FLAC 3D 与ANSYS 软件的接口程序，将在ANSYS 软件中划分的网格的节点坐标和单元信息转换为FLAC 3D 的单元数据文件(格式为FLAC 3D 文件)，在FLAC 3D 中导入该模型数据文件，所建立的模型见图4.4工程应用贵州瓮福磷矿，处于山岭地区，地形复杂，有比较复杂的地质构造，岩层数量较多，该矿有一较大向斜，对于这样一个地表形态复杂和地质构复杂的磷矿，在图3三维地质模型的ANSYS 网格图Fig.3The ANSYS mesh of the three dimensional modelA 矿体B 矿体图2矿体三维模型Fig.2The three-dimensional model of the mine body图1三维模型Fig.1The three-dimensional model26由地下转露天的开采过程中会形成大量的边坡，对这边坡的稳定性进行分析是非常有必要的.下面对该矿的某一个边坡进行建模分析.4.1边坡模型建立边坡模型的建立经过以下步骤：(1)绘制边坡等高线(图5(a))；(2)将边坡等高线放样为三维光滑曲面(图5(b))；(3)建立完整三维边坡模型(图5(c))；(4)输出为SAT 文件，在ANSYS 中进行网格划分(图5(d))；(5)通过接口程序将ANSYS 模型转换为FLAC 3D 模型(图5(e))；4.2边坡的稳定性分析边坡体为白云岩，该岩石的力学参数：粘结力为1.569MPa ，密度为2.677kg/m 3，剪张角为0°，摩擦系数为21.63°，抗拉强度为0.36MPa ，体积模量为315MPa ，剪切模量为950M Pa.对图5中所建立的边坡模型进行稳定性分析，经计算，该边坡的安全系数为3.03，该边坡处于稳定状态，边坡的位移云图见图6.5结论1)鉴于FLAC 3D 强大的岩土工程数值计算能力和建立三维模型的不足，提出了基于AutoCAD 软件平台建立模型，导入ANSYS 软件进行网格划分，通过接口程序将ANSYS 软件生成的单元信息转换为FLAC 3D 复图4FLAC 3D 模型Fig.4The FLAC 3D model120m80m40m(a)边坡等高线(b)等高线放样为三维光滑曲面(c)边坡的AutoCAD 三维模型(d)边坡的ANSYS 三维模型(e)边坡的FLAC 3D 三维模型图5建立FLAC 3D 模型Fig.5Creating FLAC 3D model图6边坡的位移云图Fig.6The contour of displacement340m300m200m100m27杂地质模型的方法.该方法发挥了AutoCAD软件的简单易用的优点，解决了FLAC3D建立复杂三维地质模型的不足，实现在FLAC3D中复杂地质模型的快速精确建模.该方法方便快捷，大大减少使用FLAC3D软件进行建模所需要的时间和精力.2)将提出的转换方法应用于贵州翁福磷矿区的模型建立中，表明：提出的转换方法以及相应的接口程序能够实现AutoCAD模型到FLAC3D模型的转换，使FLAC3D建模快捷化，从而验证了转换方法的可行性和有效性.参考文献：[1]廖秋林,曾钱帮,刘彤,等.基于ANSYS平台复杂地质体FLAC3D模型的自动生成[J].岩石力学与工程学报,2005,24(6):1010-1013.LIAO Qiulin,ZENG Qianbang,LIU Tong,et al.Automatic model generation of complex geologic body with FLAC3D based on ANSYS platform[J].Chinese Journal of Rock Mechanics and Engineering, 2005,24(6):1010-1013.(in Chinese)[2]胡斌,张倬元,黄润秋,等.FLAC3D前处理程序的开发及仿真效果检验[J].岩石力学与工程学报,2002,21(9):1387-1391.HU Bin,ZHANG Zhuoyuan,HUANG Runqiu,et al.Development of pre processing package for FLAC3D and verification of its simulating effects [J].Chinese Journal of Rock Mechanics and Engineering,2002,21(9):1387-1391.(in Chinese)[3]徐能雄,武雄,汪小刚,等.基于三维地质建模的复杂构造岩体六面体网格剖分方法[J].岩土工程学报,2006,28(8):957-961.XU Nengxiong,WU Xiong,WANG Xiaogang,et al.Approach to automatic hexahedron mesh generation for rock mass with complex structure based on3D geological modeling[J].Chinese Journal of Geotechnical Engineering,2006,28(8):957-961.(in Chinese)[4]林航,曹平,李江腾,等.基于SURPAC的FLAC3D三维模型自动构建[J].中国矿业大学学报,2008,37(2):339-342.LIN Hang,CAO Ping,LI Jiangten,et al.Automatic generation of FLAC3D model based on SURPAC[J].Journal of China University of Mining and Technology,2008,37(2):339-342.(in Chinese)[5]罗周全,吴亚斌,刘晓明,等.基于SURPAC的复杂地质体FLAC3D模型生成技术[J].岩土力学,2008,29(5):1334-1338.LUO Zhouquan,WU Yabin,LIU Xiaoming,et al.FLAC3D modeling for complex geologic body based on SURPAC[J].Rock and Soil Mechanics, 2008,29(5):1334-1338.(in Chinese)Fast generation of the complex geology body FLAC3D modelCHAI Hongbao1,LONG Xueming2(1.Hunan Provincial Key Laboratory of Safe M ining Techniques of Coal M ines,Xiangtan411201,China;2.School of Life Science,Hunan University of Science and Technology,Xiangtan411201,China)Abstract：It is difficulty to create three dimensional numerical model in FLAC3D.But it is very easy to do this with AutoCAD software,and it is convenient and accurate.The model created by AUTOCAD is exported to ANSYS,then it is meshed.The information of the element can be translated into the format which can be recognized by FLAC3D with the interface between ANSYS and FLAC3D.The model can be import into FLAC3D.Then this is a fast method to create three dimensional numerical model in FLAC3D.One model of one slope of the WENGFU phoshate is created by this method.The stability of the slope is analyzed with this model.The slope is steady.The safety factor is3.03with this situation.The efficiency of numerical simulation can be greatly enhanced with FLAC3D software in simulating complex three dimensional model.Key words：AutoCAD;ANSYS;three dimensional model;FLAC3D28。
flac3d模拟基坑开挖
3.5 计算模型计算模型见图3-3~图3-5，X轴为水平方向，Y轴为竖直方向。
本模型采用实体单元模拟土体、桩、筏板，其中素混凝土桩长5m，筏板厚2m，筏板嵌入土层0.4m。
模型中共有12730个网格点，12542个实体单元。
图3-3 计算模型图图3-4 开挖完后模型图图3-5 筏板、桩、空洞模型图3.5 模拟计算工况计算过程先对计算域进行初始应力场平衡计算，然后模拟计算地基开挖过程，最后模拟地基土的加固，并施加竖向荷载。
计算分析地基中存在空洞时上层土层开挖后产生的卸荷回弹，以及采用筏板及置换桩加固并施加压力后土层的沉降量4 计算结果及分析为便于分析空洞部位的位移应力，对模型中的4个空洞进行编号，见图4-1。
计算结果中竖向位移向上为正，向下为负；应力以拉为正，压为负。
图4-1 空洞示意图4.1 地基中不存在空洞上层土层开挖后的竖向位移分布见图4-2，由图可知，地基开挖完后在开挖面产生较大的反弹，最大回弹位移为17.8cm。
在空洞附近，回弹量比相同高程土层要大，且最大回弹位移均发生在空洞上表面，4个空洞四周的回弹位移极值见表4-1，其中1#空洞虽然埋深较深，但由于其尺寸较大，其最大回弹量与埋深较浅的2#空洞、3#空洞接近，4#空洞则由于埋深较深，且尺寸较小，其最大回弹量也相当较小。
表4-1 地基开挖后空洞四周位移极值统计图4-2 地基开挖完后竖向位移分布云图采用混凝土桩加固，并在筏板上施加荷载后地基位移变化量分布见图4-3。
由图可知，地基加固后并施加荷载后地基土产生了一定的沉降量，在场地中央的最大沉降量为3.8cm。
空洞上表面的沉降量比相同高程的土层大，下表面的沉降量则比相同高程的土层小，空洞最大沉降量均发生在上表面，最小沉降量均发生在下表面，空洞四周的位移极值统计见表4-2，1#空洞尽管其尺寸相对较大，但由于其位于场地边缘，且埋深较深，施加荷载后位移相对较小；尺寸及埋深接近的2#、3#空洞沉降量基本一致；4#空洞虽然埋深较深且尺寸较小，但由于其更接近作用力中心，故产生的沉降量与埋深较浅的2#、3#空洞基本一致。
陕北某煤矿工业场地黄土高边坡FLAC3D模拟分析
陕北某煤矿工业场地黄土高边坡FLAC3D 模拟分析作者：李明柴卓刘晓玲王念秦来源：《城市建设理论研究》2013年第01期摘要：通过对陕北某煤矿工业场地黄土高边坡进行详细勘察和大量取样分析，获得一系列有用的边坡岩土工程参数。
在分析黄土高边坡可能破坏类型及相应机理的基础上，在ANSYS 中建立坡体的数值模拟模型，经过转换将模型导入FLAC3D进行后处理，对开挖后黄土高边坡的塑性变形进行了模拟分析，得到了一些有意义的结果。
关键字：黄土高边坡变形FLAC 3D模拟分析中图分类号：X752文献标识码：A文章编号：1 工程概况陕北某煤矿工业场地由于建设需要开挖斜坡形成黄土高陡边坡。
然而，由于黄土自身特点，黄土高边坡出现变形、破坏，危机工业场地及其建筑物的安全，为此，需要对其进行加固、整治。
现场调查、勘查表明，研究区黄土高边坡，沿场地宽220~250m，顺坡长约300m，最大高度约90m，坡向345°，从上到下主要由马兰黄土（Q3eol）和离石黄土（Q2eol）组成，离石黄土中夹古土壤层，坡体范围内未见基岩出露。
工业场地设计文件中，该边坡将被开挖成10级台阶状，而煤矿主斜井井口位于开挖后边坡的坡脚，其稳定性直接影响煤矿的正常开采，故有必要对该开挖后的边坡进行稳定性评价。
2 边坡可能的破坏类型及机理分析黄土具有独特的物质结构、力学性质及其所处的特殊的地理环境，使得黄土边坡常出现剥落、冲刷、落水洞、坍塌（滑塌）、崩塌、滑坡等灾害类型[1]。
经过现场详细调查、勘察，研究边坡的主要变形破坏类型为剥落、拉裂滑塌、坡面冲刷，甚至会发生整体的滑动破坏。
剥落、拉裂滑塌、坡面冲刷均为坡面变形破坏，影响范围有限，危害性不大。
此处仅讨论坡体整体滑动——滑坡的破坏机理。
研究边坡可能形成滑坡的原因主要是坡体结构，由于边坡中有多层古土壤分布，边坡开挖后，若遇强降雨，雨水入渗坡体，会使坡体强度降低，坡体容重增加，并形成动、静水压力，使得坡体稳定性大幅度降低，最后在外部振动或其他诱因的作用下发生整体破坏。
FLAC3D实例分析教程
FLAC3D实例分析教程假设我们要分析一个简单的边坡稳定性问题。
下面是具体的步骤：1.建立几何模型：首先，我们需要建立一个几何模型，包括边坡的形状和岩土层的属性。
在FLAC3D中，我们可以通过在网格上定义顶点和连线来创建边坡的形状。
然后，我们可以设置每个区域的岩土层属性，如密度、强度和摩擦角等。
确保模型的几何和岩土层属性与实际情况相符。
2.设定边界条件：接下来，我们需要设定边界条件，即模拟中的约束和加载条件。
在边坡稳定性问题中，我们可以设定边坡底部的约束条件，如水平位移和垂直位移。
此外，我们还可以为边坡施加水平和垂直方向的荷载，模拟边坡于不同加载条件下的行为。
3.运行模拟：在完成模型和边界条件的设置后，我们可以开始运行模拟。
FLAC3D使用多线程计算，能够利用多核处理器的能力来进行快速计算。
我们可以选择设置时间步长和计算精度等参数。
模拟运行完毕后，FLAC3D将输出边坡在不同加载条件下的应力、位移和变形等结果。
4.结果分析：最后，我们需要对模拟结果进行分析和解释。
FLAC3D 提供了丰富的结果显示和分析功能。
我们可以通过绘制曲线图、生成动画和查看计算网格等方式来可视化和分析结果。
根据模拟结果，我们可以评估边坡的稳定性，并提出针对性的建议和改进方案。
在实际应用中，我们还可以使用FLAC3D的其他高级功能来进一步分析和优化边坡设计。
例如，我们可以引入土体的非线性行为模型，模拟地下水流和渗流等复杂的工程问题。
此外，FLAC3D还支持参数化建模和优化分析，可以帮助工程师迅速评估不同方案的可行性和性能。
总结起来，FLAC3D是一个强大的岩土工程分析软件，可以用于解决各种实际问题。
通过学习和应用FLAC3D的基本使用方法和分析技巧，工程师可以更好地理解和评估岩土工程问题，为工程设计和施工提供有力支持。
FLAC3D中饱和土体中的开挖
3．饱和土体中的开挖 3.1 问题陈述在不透水的基础饱和土体层中设计开挖。
土壤层1.2m 厚，潜水面水平是一个常数，与土壤顶部表面相一致。
开挖的洞穴将有一个8×8的正方形截面。
深5m 。
为本工程做准备，开挖地址已经被1m 厚的垂直封闭墙，该封闭墙深入到开挖洞穴下方2m 。
开挖以后，安装水泵用来降低开挖洞穴下部水平面。
问题是估计开挖洞穴底部由于开挖和排水引起的总体底鼓。
问题是根据对称性，分析时在三维空间取1/4部分考虑。
在土壤表面所在平面，开挖体中心为原点，定义直角坐标系统，z 轴方向指向下方。
在这个例子中，边界x =12m 和y =12m 作为对称面。
图3-1为本题的草图。
该土壤作为弹性材料。
土壤和水具有如下属性：弹性模量, K 390 MPa 剪切模量, G 280 MPa 土壤干密度, ρd 1200 kg/m 3水密度, ρw 1000 kg/m 3 墙密度, ρwall 1500 kg/m 3 渗透系数, k 10□12 m 2/Pa-s泊松比, n 0.3 流体体积模量, K f 2.0 GPa重力加速度g 约为10m/s 2。
在平衡状态的初始状态，各向同性应力为：'''4.0zzyy xx σσσ==3.2 建模过程FLAC 3D 模型尺寸为12m ×12m ×12m ，网格是尺寸为1m ×1m ×1m 的12×12×12的总体模型。
A fluid null mode l is assigned to the zones within the excavation and wall volumes , and a mechanical null model is assigned to the zones within the excavation. (See Figure 3.2 for a plot of the grid with the excavation removed.) 注意的是，在这个问题中，土壤饱和密度w d n ρρ+值与墙密度（wall density ）相等。
flac3d
flac3dFLAC3D是一款强大的数值分析软件包，广泛应用于地质工程、岩石力学和土木工程领域。
它具有极强的计算能力和模拟功能，能够模拟和分析各种地质和工程问题，如地下隧道工程、岩土体开挖、地质灾害等。
本文将介绍FLAC3D的主要功能和应用领域，并探讨其在地质和工程领域中的重要性和优势。
首先，FLAC3D具有灵活而强大的建模功能。
用户可以利用软件提供的模型构建工具，以直观的方式创建复杂的地质模型。
该软件支持各种几何和物理参数，包括岩体的几何形状、强度参数、应力状态等。
用户可以根据具体需要对模型进行灵活的调整和修改，以满足各种工程实际问题的需求。
其次，FLAC3D具有高度准确的数值计算能力。
该软件采用离散元方法（DEM）进行数值计算，能够将复杂的地质和工程问题转化为离散粒子之间的相互作用。
离散元模型能够准确地模拟材料的物理行为，如断裂、变形、变形和变形耦合等。
通过数值模拟，用户可以预测和分析地质和工程系统的行为，为工程决策提供科学依据。
第三，FLAC3D提供了广泛的分析功能。
该软件提供了丰富的分析工具，包括应力分析、位移分析、变形分析、振动分析、渗流分析等。
用户可以根据具体需要选择不同的分析方法，以获得所需的工程和地质参数。
此外，软件还提供了可视化工具，用于直观地展示模拟结果，帮助用户更好地理解模型的行为和特性。
第四，FLAC3D的广泛应用领域。
该软件在地质工程、岩石力学和土木工程等领域具有广泛的应用。
举几个例子来说明：在地下隧道工程中，FLAC3D可以模拟地下岩层的变形和稳定性，为隧道设计和施工提供指导；在岩土体开挖中，软件可模拟岩土体的塑性变形和破坏机制，评估开挖对周围环境的影响；在地质灾害预测和评估方面，FLAC3D可以模拟地震、滑坡、地表沉降等自然灾害过程，提供科学的预警和风险评估。
最后，FLAC3D的优势和重要性。
FLAC3D作为一款成熟且可靠的数值分析软件，其重要性不言而喻。
它不仅可以提供准确的工程和地质参数，还可以帮助分析人员更好地理解系统的行为和特性。
(2021年整理)用flac3d模拟基坑开挖
用flac3d模拟基坑开挖(推荐完整)编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（用flac3d模拟基坑开挖(推荐完整)）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为用flac3d模拟基坑开挖(推荐完整)的全部内容。
用flac3d模拟基坑开挖（推荐完整)编辑整理：张嬗雒老师尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布到文库,发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是我们任然希望用flac3d模拟基坑开挖（推荐完整) 这篇文档能够给您的工作和学习带来便利。
同时我们也真诚的希望收到您的建议和反馈到下面的留言区，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请下载收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为 <用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 1gen 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.1gen 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、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
interface 8 face range cylin end1 (0 0 -27.46) end2 (0 0 -28.60) radius 0.51 &
p4 (0.5 0 27.24) p5 (0 0.5 27.24) &
size 3 2 6
gen zone cyl p0 (0 0 27.24) p1 (0.5 0 27.24) p2 (0 0 21.04) p3 (0 0.5 27.24) &
p4 (0.5 0 21.04) p5 (0 0.5 21.04) &
new
set logfile cz1000
set log on
;1# 桩长28.5米,桩直径1.0米
;-------------------------------------------------
; pile in rock for example
;-------------------------------------------------
p4 (8 0 -27.46) p5 (0 8 -27.46) p6 (8 8 -23.66) p7 (8 8 -27.46) &
p8 (0.5 0 -23.66) p9 (0 0.5 -23.66) p10 (0.5 0 -27.46) p11 (0 0.5 -27.46) &
size 3 5 6 15 ratio 1 1 1 1.15 ;强风化花岗岩
size 3 6 6
gen zone cyl p0 (0 0 9.84) p1 (0.5 0 9.84) p2 (0 0 5.84) p3 (0 0.5 9.84) &
p4 (0.5 0 5.84) p5 (0 0.5 5.84) &
size 3 5 6
gen zone cyl p0 (0 0 5.84) p1 (0.5 0 5.84) p2 (0 0 2.04) p3 (0 0.5 5.84) &
macro finesand 'bulk 5.667e6 shear 3.400e6 coh 0 fric 30' ;细砂
macro silt 'bulk 3.100e6 shear 1.033e6 coh 5500 fric 8.0' ;淤泥质土
p8 (0.5 0 -2.26) p9 (0 0.5 -2.26) p10 (0.5 0 -8.46) p11 (0 0.5 -8.46) &
size 3 7 6 15 ratio 1 1 1 1.15 ;淤泥质土
group siltl range z -8.46 -2.26
gen zone radcyl p0 (0 0 -8.46) p1 (8 0 -8.46) p2 (0 0 -9.76) p3 (0 8 -8.46) &
group cohs range z -19.66 -15.36
gen zone radcyl p0 (0 0 -19.66) p1 (8 0 -19.66) p2 (0 0 -23.66) p3 (0 8 -19.66) &
p4 (8 0 -23.66) p5 (0 8 -23.66) p6 (8 8 -19.66) p7 (8 8 -23.66) &
group pile range z 1 29.60
plo sur
pause
ini z add -29.50 range z 1 29.60
; material property
;------------------------------------------------------------------------------------
p4 (0.5 0 14.14) p5 (0 0.5 14.14) &
size 3 7 6
gen zone cyl p0 (0 0 14.14) p1 (0.5 0 14.14) p2 (0 0 9.84) p3 (0 0.5 14.14) &
p4 (0.5 0 9.84) p5 (0 0.5 9.84) &
gen zone radcyl p0 (0 0 0) p1 (8 0 0) p2 (0 0 -2.26) p3 (0 8 0) &
p4 (8 0 -2.26) p5 (0 8 -2.26) p6 (8 8 0) p7 (8 8 -2.26) &
p8 (0.5 0 0) p9 (0 0.5 0) p10 (0.5 0 -2.26) p11 (0 0.5 -2.26) &
p4 (8 0 -36.50) p5 (0 8 -36.50) p6 (8 8 -28.50) p7 (8 8 -36.50) &
p8 (0.5 0 -28.50) p9 (0 0.5 -28.50) p10 (0.5 0 -36.50) p11 (0 0.5 -36.50) &
size 3 10 6 15 ratio 1 1.15 1 1.15 fill ;中风化砾岩,此部分为无桩部分的岩体
p8 (0.5 0 -19.66) p9 (0 0.5 -19.66) p10 (0.5 0 -23.66) p11 (0 0.5 -23.66) &
size 3 5 6 15 ratio 1 1 1 1.15 ;全风化花岗岩
group quanh range z -23.66 -19.66
gen zone radcyl p0 (0 0 -23.66) p1 (8 0 -23.66) p2 (0 0 -27.46) p3 (0 8 -23.66) &
group fengs range z -9.76 -8.46
gen zone radcyl p0 (0 0 -9.76) p1 (8 0 -9.76) p2 (0 0 -15.36) p3 (0 8 -9.76) &
p4 (8 0 -15.36) p5 (0 8 -15.36) p6 (8 8 -9.76) p7 (8 8 -15.36) &
p4 (0.5 0 2.04) p5 (0 0.5 2.04) &
size 3 5 6
gen zone cyl p0 (0 0 2.04) p1 (0.5 0 2.04) p2 (0 0 1.00) p3 (0 0.5 2.04) &
p4 (0.5 0 1.00) p5 (0 0.5 1.00) &
interface 5 face range cylin end1 (0 0 -15.36) end2 (0 0 -19.66) radius 0.51
interface 6 face range cylin end1 (0 0 -19.66) end2 (0 0 -23.66) radius 0.51
cylin end1 (0 0 -27.46) end2 (0 0 -28.60) radius 0.49 not
interface 9 face range cylin end1 (0 0 -28.40) end2 (0 0 -28.60) radius 0.51
;
gen zone cyl p0 (0 0 29.50) p1 (0.5 0 29.50) p2 (0 0 27.24) p3 (0 0.5 29.50) &
size 3 7 6
gen zone cyl p0 (0 0 21.04) p1 (0.5 0 21.04) p2 (0 0 19.74) p3 (0 0.5 21.04) &
p4 (0.5 0 19.74) p5 (0 0.5 19.74) &
size 3 2 6
gen zone cyl p0 (0 0 19.74) p1 (0.5 0 19.74) p2 (0 0 14.14) p3 (0 0.5 19.74) &
group zhongh range z -36.50 -27.46
gen zone reflect dd 270 dip 90
;
interface 1 face range cylin end1 (0 0 0) end2 (0 0 -2.26) radius 0.51
interface 2 face range cylin end1 (0 0 -2.26) end2 (0 0 -8.46) radius 0.51
size 3 2 6
gen zone cyl p0 (0 0 29.60) p1 (0.5 0 29.60) p2 (0 0 29.50) p3 (0 0.5 29.60) &
p4 (0.5 0 29.50) p5 (0 0.5 29.50) &
size 3 1 6
gen zone reflect dd 270 dip 90 range z 1 29.60
size 3 2 6 15 ratio 1 1 1 1.15 ;细砂
group fiБайду номын сангаасes range z -2.26 0
gen zone radcyl p0 (0 0 -2.26) p1 (8 0 -2.26) p2 (0 0 -8.46) p3 (0 8 -2.26) &
p4 (8 0 -8.46) p5 (0 8 -8.46) p6 (8 8 -2.26) p7 (8 8 -8.46) &