FLAC3D中初始地应力的生成方法

flac结构单元用法本人没有做过桩锚或桩＋支撑的围护形式，不知道pile单元在靠近开挖面一侧当土体被开挖时，pile于土体的接触会是什么样子？（1）cable单元模拟锚杆（全长连接、非全长连接，预应力锚杆），土钉等（2）beam单元模拟支撑、围檩、冠梁等（3）liner单元模拟面层，地下连续墙（4）pile单元模拟钻孔灌注桩，SMW工法，超前支护（注浆钢管、微型桩），抗滑桩等。

一般的搅拌桩可用zone＋interface解决。

FLAC学习总结一个在使用flac3d的兄弟(QQ:65006196)1、得到初始应力的方法：方法1、可以先给一些材料参数很大的值，进行初始求解，在计算之前再将材料参数设为正常值，即可。

如在手册中给的第一个示例中就是这样做的。

下面是例子，These are only initial values that are used during the development of gravitational stresses within the body. In effect, we are forcing the body to behave elastically during the development of the initial in-situ stress state.* This prevents any plastic yield during the initial loading phase of the analysis.Gen zone brick size 6 8 8Mode mohrProp bulk 1e8 shear 0.3e8 fric 35Prop cohesion 1e10 tens 1e10 ;注意在此这个值给的很大。

Init dens 1000Set gravity 0 0 -10Fix x range x -0.1 0.1Fix x range x 5.9 6.1Fix y range y -0.1 0.1Fix y range y 7.9 8.1Fix z range z -0.1 0.1set mech force=50solve;---------------------- mode null ---------------------Prop coh 1e3 tens 1e3 ;改为正常值(在此例中我们故意给小值) Mode null range x 2 ,4 y 2 , 6 z 5, 10Set largeIni xdis 0 ydis 0 zdis 0 ;清零，不影响结果，为画图方便。

ini syy=-1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
ini szz=-1e8 ran z -.1 .1
fix x y z ran z -0.1 0.1
plot cont szz
solve
apply szz=-1e8 ran z -.1 .1
plot cont szz
step 30000
最后要位移场和速度场清零的哦，用step 30000而不用solve是因为计算不收敛，呵呵，但是应力场确实跟假设比较接近。
下面是传统方法的搞法，确实只生成了自重应力场
new
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 &
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 &
size 6 6 10
model elas
pro bulk 1e10 she 1e10
ini density 2500
set gravity 0 0 -10
apply sxx=-1e9 grad 0 0 6.6666666e6 range x 59.9 60.1
apply syy=-1e9 grad 0 0 8.3333333e6 range y -.1 .1
apply syy=-1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
ini sxx=-1e9 grad 0 0 1.1111111e7 range x -.1 .1
ini sxx=-1e9 grad 0 0 6.6666666e6 range x 59.9 60.1

数值计算模型 ，可以通过多次试计算 的方式确定 速度边 界条 件与应 力边 界条件 的加载 范 围；生成 的初始 地应力 场与
矿 山深部工程地应力 实测数据 基本 一致 ；解决 了前述 ２种 边界条件加 载方法存 在的 问题 。研 究成果 可为矿 山深部工
程 的数值模拟分 析提 供参 考。
关 键 词 矿 山 深 部 工 程 初 始 地 应 力场 速 度 一应 力 边 界 法 ＦＬＡＣ
ｓｉｍｕｌａｔｉｏｎ ｏｆ ｒｏｃｋ ｕｎｄｅｒｇｒｏｕｎｄ ｅｎｇｉｎｅｅｒ ｉｎｇ ｐｒｏｖｅｄ ｔｈａｔ ｔｈｅ ｗｉｄｅｌｙ ｕｓｅｄ ｍｅｔｈｏｄｓ（ｉ．ｅ．，ｖｅｌｏｃｉｔｙ ｂｏｕｎｄａｒｙ ｍｅｔｈｏｄ ａｎｄ ｒａｐｉｄ ｓｔｒｅｓｓ ｂｏｕｎｄａｒｙ ｍｅｔｈｏｄ） ｃａｎ ｇｅｎｅｒａｔｅ ｔｈｅ ｉｎｉｔｉａｌ ｓｔｒｅｓｓ ｆｉｅｌｄ ｍｏｄｅｒａｔｅｌｙ ａｇｒｅｅｄ ｗｉｔｈ ｔｈｅ ｐｒａｃｔｉｃｅ ｉｎ ＦＬＡＣ ｐｒｏｇｒａｍ ；Ｈｏｗｅｖｅｒ，
ቤተ መጻሕፍቲ ባይዱ
力与实际不符”的问题 。针对 这些问题 ，综合前述 ２种边界 条件加 载方 法的优 点 ，提 出 了矿 山深部 工程 ＦＬＡＣ ’初 始
地应力场生成 的“速度一应力边界法 ”，探讨 了模 型速度边 界条 件与应力边界条件加载范 围的确定 方式 ，通过算例 分析
验证了“速度 一应力边界法”生成 初始地应力场 的合理 性。研究 表明 ：采用 “速 度一应力 边界法 ”构 建 的矿 山深部 工程
ｆｏｒ ｄｅｅｐ ｅｎｇｉｎｅｅｒｉｎｇ，ｔｈｅ ｇｒｏｕｎｄ ｓｔｒｅｓｓ ｇｅｎｅｒａｔｅｄ ｂｙ ｔｈｅ ｖｅｌｏｃｉｔｙ ｂｏｕｎｄａｒ ｙ ｍ ｅｔｈｏｄ ｉｓ ｎｏｔ ｃｏｎｓｉｓｔｅｎｔ ｗｉｔｈ ｔｈｅ ａｃｔｕａｌ，ａｎｄ ｗｈｅｎ Ｕ— ｓｉｎｇ ｔｈｅ ｒａｐｉｄ ｓｔｒｅｓｓ ｍ ｅｔｈｏｄ，ｔｈｅ ｂｏｕｎｄａｒｙ ｃｏｎｄｉｔｉｏｎ ｍａｙ ｂｅ ａｐｐｌｉｅｄ ｄｉｆ ｉ ｃｕｌｔｌｙ ａｎｄ ｔｈｅ ｇｅｎｅｒａｔｅｄ ｇｒｏｕｎｄ ｓｔｒｅｓｓ ａｔ ｔｈｅ ｓｈａｌｌｏｗ ｌｅｖｅｌ ｄｉｓａｇｒｅｅｓ ｗｉｔｈ ｔｈｅ ａｃｔｕａ１．Ｔｈｕｓ，ｂｙ ｃｏｍ ｂｉｎｉｎｇ ｔｈｅ ａｄｖａｎｔａｇｅｓ ｏｆ ｔｈｅ ａｂｏｖｅ ｔｗｏ ｍｅｔｈｏｄｓ，ｔｈｅ ｖｅｌｏｃｉｔｙ—ｓｔｒｅｓｓ ｂｏｕｎｄａｒ ｙ ｍｅｔｈ— ｏｄ ｔｏ ｇｅｎｅｒａｔｅ ｔｈｅ ｉｎｉｔｉａｌ Ｓｔｒｅｓｓ ｆｉｅｌｄ ｏｆ ｄｅｅｐ ｅｎｇｉｎｅｅｒ ｉｎｇ ｉｓ ｐｒｏｐｏｓｅｄ ｉｎ ＦＬＡＣ ．Ｔｈｅｎ。ｔｈｅ ｌｏａｄｉｎｇ ｒａｎｇｅ ｏｆ ｔｈｅ ｖｅｌｏｃｉｔｙ ｂｏｕｎｄａ— ｒ ｙ ｃｏｎｄｉｔｉｏｎ ａｎｄ ｔｈｅ ｓｔｒｅｓｓ ｂｏｕｎｄａｒｙ ｃｏｎｄｉｔｉｏｎ ｏｆ ｔｈｅ ｍ ｏｄｅｌ ｉｓ ｄｉｓｃｕｓｓｅｄ ｔｏ ｄｅｔｅｒｍ ｉｎｅ．Ｆｉｎａ ｌｌｙ，ｔｈｅ ｒｅａｓｏｎａｂｌｅｎｅｓｓ ｏｆ ｔｈｅ ｖｅｌｏｃｉ— ｔｙ—ｓｔｒｅｓｓ ｂｏｕｎｄａｒｙ ｍｅｔｈｏｄ ｉｓ ｐｒｏｖｅｎ ｂｙ ｔｈｅ ａｎａｌｙｓｉｓ ｏｆ ａｎ ｅｘａｍｐｌｅ．Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ：ｔｈｅ ｎｕｍｅｒ ｉｃａｌ ｍ ｏｄｅｌ ｆｏｒ ｄｅｅｐ ｅｎｇｉ— ｎｅｅｒｉｎｇ ｃａｎ ｂｅ ｂｕｉｌｔ ｂｙ ｔｈｅ ｖｅｌｏｃｉｔｙ—ｓｔｒｅｓｓ ｂｏｕｎｄａｒｙ ｍｅｔｈｏｄ ａｎｄ ｉｔｓ ｌｏａｄｉｎｇ ｒａｎｇｅ ｏｆ ｖｅｌｏｃｉｔｙ ｂｏｕｎｄａｒ ｙ ａｎｄ ｓｔｒｅｓｓ ｂｏｕｎｄａｒ ｙ ｃｏｎ— ｄｉｔｉｏｎ ｍａｙ ｂｅ ｄｅｔｅｒｍｉｎｅｄ ｔｈｒｏｕｇｈ ｃｏｍｐｕｔｉｎｇ ｍａｎｙ ｔｉｍｅｓ；ｔｈｅ ｉｎｉｔｉａｌ ｓｔｒｅｓｓ ｆ ｉｅｌｄ ｇｅｎｅｒａｔｅｄ ｂｙ ｔｈｅ ｖｅｌｏｃｉｔｙ－－ｓｔｒｅｓｓ ｂｏｕｎｄａｒ ｙ ｍ ｅｔｈ·－ ｏｄ ｃｏｎｓｉｓｔ ｗｉｔｈ ｔｈｅ ｍｅａｓｕｒｅｄ ｄａｔａ；Ｔｈｅ ｐｒｏｂｌｅｍｓ ｅｘｉｓｔｉｎｇ ｉｎ ｔｈｅ ａｂｏｖｅ ｔｗｏ ｍｅｔｈｏｄｓ ａｒｅ ｓｏｌｖｅｄ．Ｔｈｅ ｒｅｓｕｌｔｓ ｃａｎ ｓｅｒｖｅ ａｓ ｒｅｆｅｒ— ｅｎｃｅ ａｎｄ ｇｕｉｄａｎｃｅ ｆｏｒ ｔｈｅ ｎｕｍｅｒｉｃａｌ ｓｉｍ ｕｌａｔｉｏｎ ａｎｄ ａｎａｌｙｓｉｓ ｏｆ ｄｅｅｐ ｍｉｎｅ ｅｎｇｉｎｅｅｒ ｉｎｇ．

《FLAC/FLAC 3D基础与工程实例》全部命令流1gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 &p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 ＆size 10 10 10 rat 1。

0 0。

9 1.1 group brick_1gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 ＆size 10 10 10 rat 1。

0 0.9 1.1 group brick_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1。

0 0。

9 1.1 group brick_3plot sur2—1newgen zon bri size 3 3 3model elasprop bulk 3e8 shear 1e8ini dens 2000fix z ran z -.1 .1fix x ran x —.1 .1fix x ran x 2。

9 3。

1fix y ran y —.1 。

1fix y ran y 2.9 3.1set grav 0 0 -10solveapp nstress -10e4 ran z 3 x 1 2 y 1 2hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve3—1；———----——----—--—-—-——-—-—-——--------工程信息;Project Record Tree export；Title:Simple test;-—-—-—-——-—---———-—---—-——--——--—计算第一步;。

.. STATE： STATE1 ...。

第7章初始地应力场的生成及应用在土木工程或采矿工程领域中，初始地应力场的存在和影响不容忽略，它既是影响岩体力学性质的重要控制因素，也是岩体所处环境条件下发生改变时引起变形和破坏的重要力源之一。

因此，要想较真实地进行工程模拟仿真，就必须保证初始地应力场的可靠性。

初始地应力场生成的主要目的是为了模拟所关注分析阶段之前岩、土体已存在的应力状态。

本章即介绍FLAC3D中初始地应力场的生成方法及应用。

本章重点：✓常用的初始地应力场生成方法✓常见工程初始地应力场的生成✓路基施工过程的模拟7.1 初始地应力场生成方法在FLAC3D中，初始应力场的生成办法较多，但通常用的是以下三种方法，即弹性求解法、改变参数的弹塑性求解法以及分阶段弹塑性求解法。

下面将以表7-1所述简单模型为例，介绍这三种生成初始地应力场的方法。

表7-1 模型尺寸、土体密度及变形参数1×1×2 1×1×2 2000 30 10 0.35 7.1.1 弹性求解法初始地应力的弹性求解法生成是指将材料的本构模型设置为弹性模型，并将体积模量与剪切模量设置为大值，然后求解生成初始地应力场。

例叙述的是采用该法生成上述简单模型的初始地应力场的过程。

例7.1弹性求解生成初始地应力场newgen zone brick size 1 1 2 model elasprop bulk 3e7 shear 1e7 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1 ;开始一个新的分析;生成网格模型;设置弹性本构模型;设置体积模量和剪切模量;固定z=0平面所有节点z向速度;固定x=0平面所有节点x向速度;固定x=1平面所有节点x向速度;固定y=0平面所有节点y向速度;固定y=1平面所有节点y向速度ini dens 2000 set grav 0 0 -10 solve;设置密度 ;设置重力加速度 ;按软件默认精度求解图 7-1为运行上述命令文件后得到的初始地应力场应力云图。

设置初始应力的弹塑性求解：
gen zon bri size 1 1 2 model mohr prop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 ini szz -40e3 grad 0 0 20e3 ran z 0 2 ini syy -20e3 grad 0 0 10e3 ran z 0 2 ini sxx -20e3 grad 0 0 10e3 ran z 0 2 set grav 0 0 -10 solve
4、边界条件及初始条件
在FLAC3D中，包含多种边界条件，边界方位 可以任意变化，边界条件可以是速度边界、应力边 界，单元内部可以给定初始应力，节点可以给定初 始位移、速度等，还可以给定地下水位以计算有效 应力等。这众多的边界条件主要通过apply或fix命 令来进行设置。而初始条件则主要通过initial命令 来执行，对所提的这两个命令必须严格区分并了解 其差异。通常我们所计算的模型均采用力学边界， 初始条件也基本是初始地应力的输入，对此两种不 同的力，其设置存在差别，同时在计算过程中，该 二者的变化情况也各不相同。
对于这两种基本的 网格，其公共面上的 关键点的对应关系更 需校核好，否则将出 现杂乱错误的网格。
对此马蹄形隧道，其公 共面处，p0 — p0,p1—p3, p2—p2,p4—p5 , p8—p9,p10 —p11
对于对称的模型也可以采 用镜像命令：
gen zone reflect norm -1 0 0 & origin 0,0,0
对于任何形状的单元体， 其建立单元模型时关键

• size 6 6 10
• •
mproodeblulelkas10e10常sh规e 法10e10
• ini den 2500
• ini sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1
• ini sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1
FLAC3D 5.0培训教程（武汉）
工程师 李振 2014.3.27-3.28 Itasca(武汉)安排
2014.3.27~ 2014.3.28
1. FLAC3D V5.0界面操作 2. FLAC3D基本操作方法vs应用流程；
initial stress 3. FLAC3D内置Fish语言的应用； 4. FLAC3D结构单元vs接触单元； 5. FLAC3D渗流模块 6. 其他
精品课件
几种形成初始应力的方法
1. 弹性求解的方法 2. 更改强度参数的弹塑性求解 3. 分阶段弹塑性求解(mohr—solve elas) 4. 存在静水压力的初始应力（水下构筑物） 5. 深埋地应力场
精品课件
1、弹性求解的方法
精品课件
2、更改强度参数的弹塑性求解
精品课件
2、更改强度参数的弹塑性求解
• ini syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1
• ini syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1精品课件
精品课件
深埋工程地应力场
精品课件
• size 6 6 10
• model elas
• pro深bu埋lk工10程e10地sh应e 力10e场10 —SB法

３ 初 始 地应 力场 反演 过程
射 。它是 一种 快 速梯 度下 降 的方 法 ， 目的是使 实 际输 出
和 预期输 出之 间 的误 差 的平 方 最 小 。 通 过 调 整 Ｂ Ｐ网 络 中的连 接权 值 、 网络 规模 （ 括 、 和 隐含 层 神 经元 包 ｍ
结合 Ｂ Ｐ神 经 网络 和有 限差 分程 序 Ｆ ＡＣ Ｄ，初 始 Ｌ ３
第 １ 含 层 隐
第 ｎ隐 含层
确选 取地 应力 的量 值 和方 向是 设计 工 作 的关键 ， 始地 初
应力 场 的 确 定 历 来 是 岩 石 力 学 的 一 个 重 要 研 究 分 支 ‘。 ‘ 。初始 地 应力 场是 随着 时 间 、 间不断 变 化 的非 空
稳定 场 ， 人们 很难 了解 其 内部 构 造 ， 能 通 过 它 与外 部 只 的联 系 （ 如 输人 一 出关 系 ） 把 握 其 信息 。 因此 初 例 输 来
凌影 ： 于 Ｆ Ａ ３ 基 Ｌ Ｃ Ｄ及 Ｂ Ｐ神 经 网 络 的厂 区初 始 地 应 力 场反 演 分 析
新 的 速 度 和 位 移
图 ２ Ｆ Ａ ３ 的计 算 循 环 示 意 ＬＣＤ
按正交 或均匀 设 计 原理 构 建 参 数组 合 表 ，并 通 过 Ａ — Ｎ
Ｓ Ｓ Ｌ Ｃ等数 值计算 软件 或其 它一 些计 算 Ｔ 具获 得 Ｙ 、Ｆ Ａ 量 测位置 的计算 地 应力 值 信 息 ，从 而建 立 神 经 网络 学 习与训 练的样本 ；② 将 计算 地 应 力值 作 为 网 络 的输 人
ｉｎＡ ａ ｓ ｆ ｏ ｔ ｕ Ｄ ｍ ｎｉ ｓ 是 由美 国 Ｉ ｓａ ａ ｎ ｌ ｉ ｏ ｎｉ ａｉ ３ ｉ ｅ ｓ ｎ ） ｙｓ Ｃ ｎ ｎ ｏ ｔ ｃ ａ

Itasca Consulting Group, Inc. Minneapolis, MN USA
12:30
FLAC3D 3.00
Step 163 Model Perspective 11:24:07 Fri Jun 10 2011 Center: X: 5.000e-001 Y: 5.000e-001 Z: 1.000e+000 Dist: 6.030e+000 Rotation: X: 20.000 Y: 0.000 Z: 40.000 Mag.: 1 Ang.: 22.500
初始地应力场的生成及应用
12:30
学习内容

常用的初始地应力场生成方法 常见工程初始地应力场的生成


路基施工过程的模拟
12:30
4.1 初始地应力场生成方法
在土木工程或采矿工程领域中，初始地应力场 的存在和影响不容忽略，它既是影响岩体力学性质 的重要控制因素，也是岩体所处环境条件下发生改 变时引起变形和破坏的重要力源之一。因此，要想 较真实地进行工程模拟仿真，必须保证初始地应力 场的可靠性。初始地应力场生成的主要目的是为了 模拟所关注分析阶段之前岩、土体已存在的应力状 态。本章介绍FLAC3D中初始地应力场的生成方法 及应用。
12:30
4.1.3 分阶段弹塑性求解法
该方法是分成两个阶段求解:首先，程序自 动将模型所有组成材料的粘聚力和抗拉强度分 别设置为较大值，进行弹性求解，直至体系达 到力平衡状态，接着将粘聚力和抗拉强度重置 为初始设定值进行弹塑性阶段的求解，直至体 系达到力平衡状态。 该方法只适合计算模型采用MC模型的情况。
12:30
n gen zone brick size 1 1 2 model mohr prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ini dens 2000 set grav 0 0 -10 solve elas plo con sz

以应力回归分析方法为基础针对flac3d程序施加荷载以及处理边界条件方面的特殊性和传统地应力场反演方法的局限性提出一种基于连续加载法的flac3d三维初始地应力场反演的改进方法给出了该方法的具体实现步骤通过简单算例分析验证了该方法的优越性并将该方法应用于双江口水电站大型地下厂房区域地应力场反演回归分析回归结果与实测值基本相符说明该方法是可行的可获得合理可靠的三维地应力场并具有易于实现的特点
别对应于 6个应力分量; L i 为相应于自变量的回归系数; n为分 项荷载模式总数。
假定有 m 个观测点, 则最小二乘残差平方和为
m6
n
E E E S残 =
( R*jk -
L
i
R
i jk
)2
( 2)jk 为 k 观测点的第 j 个应力分量的观测值。
根据最小二乘法 原理, 使 S残 为最小值的 方程式为
收稿日期: 2010-02-11 基金项目: 国家自然科学基 金资助项目 ( 50809017) ; 中国 博士后科学基金资助 项目及特别资助项目 ( 20080440812, 200902354) 。 作者简介: 苏国韶 ( 1973) ) , 男, 广西南宁人, 副教授, 博士后, 主要 研究方向为 水利水电工程。
回归结果与实测值基本 相符, 说明该方法是可行的, 可获得合理可靠的三维地应力场, 并具有易于实现的特点。
关 键 词: 地下厂房; 地应力场; 回归分析; FLAC3D
中图分类号: TV 512
文献标识码: A
do :i 10. 3969 / .j issn. 1000-1379. 2011. 02. 061
E ER
1 jk
R
2 jk
k= 1 j= 1

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的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。

调整三维网格中的多面体单元来拟合实际的结构。

单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。

FLAC3 D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。

由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。

三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。

三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。

三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。

FLAC3D各种命令笔记编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（FLAC3D各种命令笔记）的内容能够给您的工作和学习带来便利。

同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。

本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为FLAC3D各种命令笔记的全部内容。

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 剪胀角FISH语言：Degrade;∏/180Pi;∏Ngp；节点总数Nzone;单元体总数单元节点遍历：P_z = zone_headLoop while p_z ＃ null………P_z = z_next(p_z)End loopRead（ar，n）;将数组ar的前n个记录写到文件print zone state 显示塑性区域编写fish函数时，尽量采用浮点型，例如1。

关键词 ：初始地应 力场 ；ｒＥＡｄ。；分析
１ 概 述
初始地应力 场是 围岩稳定 的重要影 响因素之 一 ， 也是工程设计的重要初始条件之一 ，因此在计 算中采取的初始应力场是否合理 ，将对工程的可靠 性 和安 全性 产生 直接 影 响 …１。现 阶段 初 始地 应力 场 的模拟方法 主要有位移反 分析法和应 力 回归 分析 法 ，前 者通 过对 大地 构造 变 形 的观测 或 工程 施工 中 的位 移监 测 对工 程 区域 的应 力场 进行 反 演研 究 。后 者依据对工程区域地应力场产生条件 的规律性认识 建立三维模型 ，以实测地应力值为基础 ，在模型边 界上施加不同组 合的荷载 和约束 ，通过数 值计算 ， 使得计算值与实测值达到最优拟合 ，以求得工程 区 域初始地应力场 ，这科Ｉ方法可以很好的模拟 出构造 应力场［１－２］。本文结合某水电站工程 ，基于岩土工
摘 要 ：介绍利 用 Ｌ Ｄ进行初 始地 应力场模拟 分析 中的几种方法 ，着重描述 了生成构造应 力场所采 用的几种边界 条 件 ，借助这几种 不同的方法分别对 某水电站调 压 室初始地应 力场进行模拟 分析 ，并将 其结 果与 实测值进 行 比较 ，对每 一种 方案的 实际模拟 效果进行分析 ，为后 续的工程分析提供 可靠依 据。 图 ８幅，表 ４个。
Ｉｉ ｃ０ｓ ｓｉｎ口ｉ
３ 数值计算及结果比较
３．１ 计算模型及计算范围 根据计算要求和一般经验 ，对于地下洞室进行
数值计算时模型范围选取的一般原则是 ：引水隧洞 一 般 模 型 范 围大 于 ３—５倍 的洞 径 ，对 于 大 型 的地
ｍ ｉ＝ ｃｏｓ ｉｃｏｓ口 ｉ
ＩＺｉ＝ｓｉｎｆｌｉ 式 中 ，Ｚｉ、 ｍｉ、／／，ｉ分别 为 ｉ对 ，ｌ，，ｚ轴 的 方向余弦。口为主应力在水平面上 的投影方位角 ， 以北为 ０，顺时针旋转。 为主应力 的倾角 ，仰角 为 正 。

Flac3D相关命令_加油_⼤家⼀起学习!1. 基础实例分析命令：Gen zone ……；model ……；prop ……（材料参数）；set grav 0,0,-9.81（重⼒加速度） plot add block group red yellow 把在group 中的部分染成红⾊和黄⾊plot add axes black 坐标轴线为⿊⾊；print zone stress% K 单元应⼒结果输出 ini dens 2000 ran z a b （设置初始密度，有时不同层密度不同）；ini ……（设置初始条件）；fix ……（固定界⾯）hist unbal 监测不平衡⼒，并保留历史记录；save t1.sav 保存到⽂件t1；hist gp zdisp 4,4,8监测⽹格坐标点（4，4，8）在z ⽅向的位移，并保留历史纪录 plot set rot 20 0 30视图的旋转⾓度为（20，0，30）PLOT block group 查看计算模型 Ctrl+G 变成灰⾊图像Ctrl+Z 模型某个部分进⾏放⼤ Home 、End 拉近远离图像X ，Y ，Z 旋转 Shift+ X ，Y ，Z 反向旋转set plot jpg ；set plot quality 100 ；plot hard file 1.jpg 图像输出（格式、像素、名称） plot set magf 1.0视图的放⼤倍数为1.0；plo con szz z ⽅向应⼒云图2. ini z add -1 range group one 群one 的所有单元，在z ⽅向上向下移动1m ；然后合并命令 gen merge 1e-5 range z 0此命令是接触⾯单元合并成⼀个整体，1e-5是容差3. (基坑开挖步骤)：Step 1: create initial model state （建⽴初始模型）Step 2: excavatetrench （开挖隧道）4. group Top range group Base not 定义（群组Base 以外的为）群组Top5. plot blo gro 使得各个群组不同颜⾊显⽰6. （两个部分间设置界⾯；切割法）：gen separate Top 使两部分的接触⽹格分离为两部分；interface 1 wrap Base Top 在（Base 和Top ）这两部分之间添加接触单元；plot create view_int 显⽰，并创建标题view_int ；plot add surface 显⽰表⾯；plot add interface red 界⾯颜⾊红⾊7. （简单的定义函数及运⾏函数）new ；def setup 定义函数setup ；numy = 8定义常量numy 为8；depth = 10.0 定义depth 为10；end 结束对函数的定义；setup 运⾏函数setup8. （隧道⽣成）上部圆形放射性圆柱及下部块体单元体的建⽴，然后镜像。

flac怎样输出制定节点的应力值怎样输出制定节点的应力值我计算一个模型有几万个节点，我怎么才能把指定的节点应力输出。

指定节点的应力？到底是节点的力嘛，还是单元的应力哦？两个办法：1. 先plot block model id on，然后找到你想看的那个单元的id，假设是10，然后就plot con szz out on range id 10；2. 用fish，通过gp_head，gp_near，gp_id等等获得你要的节点，然后可以得到gridpoint的各种信息。

通过flac的fish语言，写节点信息和相关的力到文件既然，ansys的格式文件能够导入到flac文件，那么用fish肯定能够导出相关的数据，然后再看处理结果。

如题，我建议，通过flac的fish语言，写节点信息和相关的力，请帮主支持。

如果能够把相关的信息导入到文件，那么就很容易的对数据进行处理。

我先贴一个：用fish写文件的，当然也是从这个网络上找到的，我解释一下，也不算直接抄袭。

;; Initialization;初始化def ini_mesh2tecIO_READ = 0;定义读文件关闭吧IO_WRITE = 1;定义写文件IO_FISH = 0;定义一个IO_ASCII = 1;定义用asci码iarray buf(1);定义一个数组，tec_file = 'tec10.dat';打开文件的名字zone_ngp = z_numgp(zone_head);zone 的标志数endini_mesh2tec;; Write Tecplot File Headdef write_mesh_headbuf(1) = 'TITLE = "Flac3d Mesh to Tecplot Version 10" \n'buf(1) = buf(1) + 'VARIABLES = "X" "Y" "Z" "XDISP" "YDISP" "ZDISP"\n' buf(1) = buf(1) + 'ZONE T="Tecplot v10" \n' buf(1) = buf(1) + ' N=' + string(ngp) + ','buf(1) = buf(1) + ' E=' + string(nzone) + ','if zone_ngp = 4 thenbuf(1) = buf(1) + ' ZONETYPE=FETETRAHEDRON \n'elsebuf(1) = buf(1) + ' ZONETYPE=FEBrick \n'endifbuf(1) = buf(1) + ' DATAPACKING=POINT \n'buf(1) = buf(1) + ' DT=(SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE)' status = write(buf,1)end;; Write Grid Point Coordinates and Displacements;写节点信息def write_gpp_gp = gp_headloop while p_gp # null;这个循环，叫做遍历所有节点buf(1) = string(gp_xpos(p_gp))buf(1) = buf(1) + ' ' + string(gp_ypos(p_gp)) buf(1) = buf(1) + ' ' + string(gp_zpos(p_gp)) buf(1) = buf(1) + ' ' + string(gp_xdisp(p_gp)) buf(1) = buf(1) + ' ' + string(gp_ydisp(p_gp)) buf(1) = buf(1) + ' ' + string(gp_zdisp(p_gp));这就是传说中的写status = write(buf,1)p_gp = gp_next(p_gp)endloopend;; Write Zone Connectivity;写单元信息def write_zonep_z = zone_head;获得单元的最高号码loop while p_z # null;循环，或者叫做遍历buf(1) = ''if zone_ngp = 4 thenbuf(1) = buf(1) + string(gp_id(z_gp(p_z, 1))) + ' ' buf(1) = buf(1) + string(gp_id(z_gp(p_z, 2))) + ' ' buf(1) = buf(1) + string(gp_id(z_gp(p_z, 3))) + ' ' buf(1) = buf(1) + string(gp_id(z_gp(p_z, 4))) + ' ' elsebuf(1) = buf(1) + string(gp_id(z_gp(p_z, 1))) + ' ' buf(1) = buf(1) + string(gp_id(z_gp(p_z, 2))) + ' ' buf(1) = buf(1) + string(gp_id(z_gp(p_z, 5))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 3))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 4))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 7))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 8))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 6))) + ' 'endifstatus = write(buf,1)p_z = z_next(p_z)endloopend;; Main Functiondef mesh2tecstatus = closestatus = open(tec_file,IO_WRITE,IO_ASCII)if status = 0 thenwrite_mesh_headwrite_gpwrite_zonestatus = close;关闭文件ii = out('Successfully Write Data Into File ' + tec_file);显示信息elseii=out('Open File Error! Status = ' + string(status))endifendmesh2tec;调用把所有的记录文件写到文件中，名字不一样，那出来给大家共享，如果谁知道谁自动获得最大的记录的编号，那么这个程序就真正完美了，不过，可以凑合一下哦文件名是循环生成的哦，最大数是300，请使用的时候，修改这个值 loop iixx_ (1,300)上面的300，是根据自己的项目变化的哦def write_logloop iixx_ (1,300)logfile11= 'log' + string ( iixx_ ) + '.his'commandhis write iixx_ file logfile11end_commandend_loopendwrite_log自动生成的，乖乖，好多哦一个一个的看吧一下是异想天开几点：1.一个文件写n行，现在是一个文件写2行2.直接导入到excel，然后就作图3.想办法和上面的那个loop结合起来，导入到能够画等值线的程序中，如ansys，显示结果set log onpr histset log off然后算一下有多少行就知道有多少纪录了如何寻求cable中的应力值试试plot add sel cable grout可以通过直径显示应力大小，加颜色可以表示出应力方向补充一下plot add sel cable grout stress首先谢谢ying111您所说的是，从图像上确定cable的受力性状；如何确定cable中沿轴线方向任何一点的具体应力值，能否说的更清楚一些！请赐教。

FLAC3D中初始地应⼒的⽣成⽅法
《FLAC/FLAC3D基础与⼯程实例》——陈育民
FLAC3D初始地应⼒场的⽣成通常⽤三种⽅法：弹性求解法，改变参数的弹塑性求解法，分阶段的弹塑性求解法。

弹性求解法是指，弹性模型为材料本构模型，将体积模量与剪切模量设置为⼤值，求解⽣成初始地应⼒场。

其结果与公式“主应⼒（zz⽅向）=密度*g*⼟层深度”，“主应⼒（xx⽅向）=主应⼒（yy⽅向）=主应⼒（zz⽅向）*泊松⽐/1-泊松⽐”的计算结构⼀致。

此法常⽤于浅埋⼯程和地表⼯程数值模拟时的初始地应⼒场⽣成。

改变强度参数的弹塑性求解法指，求解过程中始终采⽤塑性模型，但为防⽌在计算过程中出现屈服区域，将粘聚⼒和抗拉强度设为⼤值，计算⾄平衡后，再将粘聚⼒和抗拉强度改为分析所采⽤的值计算⾄最终平衡状态。

分阶段的弹塑性求解法与前述改变强度参数的弹塑性求解法计算⽅法⼀致，⽬前，在FLAC3D中，此法只适合计算模型采⽤摩尔-——库仑模型的情况。

该求解过程分两个阶段进⾏：⾸先，程序⾃动将模型所有组成材料的粘聚⼒和抗拉强度分别设置为较⼤值，进⾏弹性求解，直⾄平衡；接着将粘聚⼒和抗拉强度重置为初始设定值进⾏塑性阶段的求解，直⾄平衡。

FLAC3D500培训教程1.引言FLAC3D500是一款基于三维快速拉格朗日法的岩土工程数值分析软件，广泛应用于岩土工程、地质工程、矿业工程等领域。

本教程旨在帮助用户了解FLAC3D500的基本操作和功能，为实际工程问题提供有效的数值模拟解决方案。

2.FLAC3D500软件安装与启动2.1软件安装请确保您的计算机满足FLAC3D500的运行要求。

然后，从官网FLAC3D500安装包，按照提示完成安装。

2.2软件启动安装完成后，在开始菜单中找到FLAC3D500，启动。

软件启动后，您将看到主界面。

3.FLAC3D500基本操作3.1创建新项目“文件”菜单，选择“新建项目”，在弹出的对话框中输入项目名称，“确定”创建新项目。

3.2导入模型“文件”菜单，选择“导入模型”，在弹出的对话框中选择模型文件（.flac3d或.f3grid），“打开”导入模型。

3.3设置模型参数在“模型”菜单中，可以设置模型的基本参数，如材料属性、边界条件、初始应力等。

3.4创建网格在“网格”菜单中，可以创建和编辑网格。

选择“创建网格”，在弹出的对话框中设置网格参数，“确定”网格。

3.5设置分析类型在“分析”菜单中，选择分析类型（如静态分析、动态分析等），并设置相应的分析参数。

3.6运行分析在“分析”菜单中，选择“开始分析”，软件将开始计算。

计算过程中，您可以在“输出”菜单中查看计算结果。

3.7结果查看与导出分析完成后，您可以在“输出”菜单中查看计算结果，如位移、应力等。

还可以将结果导出为文本、图片等格式。

4.FLAC3D500高级功能4.1参数化分析通过参数化分析，可以方便地研究不同参数对计算结果的影响。

在“分析”菜单中，选择“参数化分析”，设置参数范围和步长，“开始分析”进行计算。

4.2剖面分析剖面分析可以帮助用户更好地了解模型内部的应力、位移等分布情况。

在“分析”菜单中，选择“剖面分析”，设置剖面位置和方向，“开始分析”进行计算。