FLAC3D学习

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的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。

FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发变形和移动（大变形模式）。FLAC3D采用的显式拉格朗日算法和混合-离散分区技术能够非常准确发模拟材料的塑性破坏和流动。由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。FLAC3D采用ANSI C++语言编写的。

FLAC3D有以下几个优点：

1 对模拟塑性破坏和塑性流动采用的是“混合离散法”。这种方法比有限元法中通常采用的“离散集成法”更为准确、合理。

2 即使模拟的系统是静态的，仍采用了动态运动方程，这使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍。

1

3 采用了一个“显式解”方案。因此，显式解方案对非线性的应力-应变关系的求解所花费的时间，几乎与线性本构关系相同，而隐式求解方案将会花费较长的时间求解非线性问题。面且，它没有必要存储刚度矩阵，这就意味着，采用中等容量的内存可以求解多单元结构； 模拟大变形问题几乎并不比小变形问题多消耗更多的计算时间，因为没有任何刚度矩阵要被修改。 当然，

它也存在以下几个不足之处：

1 对于线性问题的求解，FLAC3D比有限元程序运行得要慢；因此，当进行大变形非线性问题或模拟实际可能出现不稳定问题时，FLAC3D是最有效的工具。

2 用FLAC3D求解时间取决于最长的自然周期和最短的自然周期之比。但某些问题对模型是无效的。

Flac3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了10种材料模型：

1. 开挖模型

2. 3个弹性模型（各向同性，横观各向同性和正交各向同性弹性模型） 3 6个塑性模型（Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型）。

还包含了节理单元，也称为界面单元，能够模拟两种或多种材料界面不同材料性质的间断特性。节理允许发生滑动或分离，因此可以用来模拟岩体中的断层、节理或摩擦边界。

FLAC3D中的网格生成器gen，通过匹配、连接由网格生成器生成局部网格，能够方便地生成所需要的三维结构网格。还可以自动产生交叉结构网格（比如说相交的巷道），三维网格由整体坐标系x，y，z系统所确定，不同于FLAC程序是由行列方式确定。这就提供了比较灵活的产生和定义三维空间参数。

边界条件和初始条件：

定义方式与FLAC相同。在边界区域可以指定速度（位移）边界条件或应力（力）边界条件。也可以给出初始应力条件，包括重力荷载以及地下水位线。所

2

有的条件都充许指定变化梯度。

FLAC3D还包含了模拟区域地下水流动、孔隙水压力的扩散以及可形的多孔隙固体和在孔隙内粘性流动流体的相互耦合。流体被认为是服从各向同性的达西定律。流体和孔隙固体中的颗粒是可变形的，将稳态流处理为紊态流可以模拟非稳态流。同时能够考虑固定的孔隙压力和常流的边界条件，也能模拟源和井。流体模型可以与结构的力学分析独立进行。

FLAC3D做计算分析的一般步骤：

与大多数程序采用数据输入方式不同，FLAC采用的是命令驱动方式。命令字控制着程序的运行。在必要时，尤其是绘图，还可以启动FLAc用户交互式图形界面。为了建立FLAC计算模型，必须进 行以下三个方面的工作：

1. 有限差分网格 2. 本构特性与材料性质 3. 边界条件与初始条件

完成上述工作后，可以获 得模型的初始平衡状态，也就是模拟开挖前的原岩应力状态。然后，进行工程开挖或改变边界条件来进行工程的响应分析，类似于FLAC的显式有限差分程序的问题求解。与传统的隐式求解程序不同，FLAC采用一种显式的时间步来求解代数方程。进行一系列计算步后达到问题的解。

在FLAC中，达到问题所需的计算步能够通过程序或用户加以控制，但是，用户必须确定计算步是否已经达到问题的最终的解。 最后，进行结果的分析与总结，这就是用户自己的事情了。

3

2. 关于flac3d—内置fish语言精讲

FISH语言是FLAC3D程序的内置编程语言。

因为FLAC3D的最佳操作方式是命令流文件方式。/这一点与ANSYS很相似，而FISH就相当于ANSYS的APDL语言。它包括循环、判断等结构。如果你还用过其它高级语言，那么从形式上讲你也可以把它理解为子函数。FISH语言的引入极大的方便了用户进行复杂的程序建模。它不但可以嵌入命令流文件里工作而且还可以引用FLAC3D本身的任何命令。所以说它实现了对FLAC3D的完全控制引入FISH语言的作法值得其它通用软件效仿。

Fish函数增加了以下几种新的特性：

1.增加了fish变量来获取结点、单元和界面变量 2.fish提供了获取结构单元变量的途径

3.休单元和面单元性质目前可以通过单元变量名z_prop(i_z string)和界面单元名i_prop(i_z string)分别加以识别

4.fish函数可以获取单元应变和应力速率，还提供了全应变增量张量和应变速率张量

5.提供了fish绘图子程序函数够生用户定义的图形内容

FLAC3D是一个强大的软件，但是不得不承认的，它的界面没有ansys好用，商业运行不是很强的说。

注意点：

1.fish函数可以嵌套使用;

2.以save命令保存模型时，fish函数和变量也同时保存；

3.fish函数不支持缩写，这与flac3d命令不同，另外所有的fish函数或变量不区分大小写，程序同意转化为大写进行编译，当然也可以通过执行set case-sensitivity on来区分大小写；

4.变量或函数名不能以一个数字开头也不能是下列字

符： . ， * / ^ = > < # ( ) [ ] @ ; “

5.如果用命令set safe on 指定了编译安全模式，则用户调用fish函数时，函数名前必须加@；

6.如果变量不曾赋值，则系统默认为零（整形），如果赋值，其类型由值的

4

类型决定；

7.fish函数的调用方法：

? 可以出项在其他fish函数的单独行中； ? 可以出现在其他fish函数的表达式中； ? 出现在flac3d的命令行中 ? 作为命令set，print，hist的参数。 8.fish函数不支持递归调用;

9.fish函数可以重新定义，其中变量初始化为原先值； 10.数组名不得和已定义fish函数同名（和参变量不同）；

11.在flac3d命令提示下，set命令可以用已存在的参数对另一参数进行赋值，但数组类不支持。

5

3. FLAC/FLAC3D的FAQ常见问题问答

1. FLAC3D是有限元程序吗？ 答：不是！是有限差分法。 2. 最先需要掌握的命令有哪些？

答：需要掌握gen， ini， app， plo， solve等建模、初始条件、边界条件、后处理和求解的命令。

3. 怎样看模型的样子？

答：plo blo gro可以看到不同的group的颜色分布 4. 怎样看模型的边界情况？ 答：plo gpfix red sk 5. 怎样看模型的体力分布？ 答：plo fap red sk 6. 怎样看模型的云图？

答：位移：plo con dis (xdis， ydis， zdis) 应力：plo con sz (sy， sx， sxy， syz， sxz) 7. 怎样看模型的矢量图？ 答：plo dis (xdis， ydis， zdis) 8. 怎样看模型有多少单元、节点？ 答：plo info

9. 怎样输出模型的后处理图？

答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件 10. 怎样调用一个文件？ 答：File/call或者call命令 10. 如何施加面力？ 答：app nstress

11. 如何调整视图的大小、角度？

答：综合使用x， y， z， m， Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键

12. 如何进行边界约束？

6

答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束） 13. 如何知道每个单元的ID？

答：用鼠标双击单元的表面，可以知道单元的ID和坐标 14. 如何进行切片？

答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 15. 如何保存计算结果？ 答：save +文件名

16. 如何调用已保存的结果？ 答：rest +文件名；或者File / Restor 17. 如何暂停计算？ 答：Esc

18. 如何在程序中进行暂停，并可恢复计算？ 答：在命令中加入pause命令，用continue进行继续 19. 如何跳过某个计算步？

答：在计算中按空格键跳过本次计算，自动进入下一步 20. Fish是什么东西？

答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能

21. Fish是否一定要学？

答：可以不用，需要的时候查Mannual获得需要的变量就可以了 22. FLAC3D允许的命令文件格式有哪些？ 答：无所谓，只要是文本文件，什么后缀都可以 23. 如何调用一些可选模块？

答：config dyn (fluid， creep， cppudm) 24 .如何在圆柱体四周如何施加约束条件？

答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2

25. 如何能把一个PLOT的图像数据导出来以便用其他软件绘图？

7

答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

26.关于dd、dip的理解：

答：dip-direction angle， measured in the global xy-plane clockwise from the positive y-axis 这个是从y轴到x顺时针方向所转过的角度；dip dip angle， measured in the negative z-direction from the global xy-plane 这个是xy平面与z轴负方向的夹角。

dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线 ，dd表示y轴正向顺时针到那条射线的夹角。首先应该按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．

27. 如何调用FISH文件？

答：新手调用FISH文件出错的原因是：文件名为 fishcall.fis.txt(因为在记事本里这个txt的后缀是默认的）。修改方法是：将 fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存，保存时候将保存类型改为所有文件(*，*)，这时你便可在文件名栏内输入 fishcall.fis，类型也就正确了。

28. 如何理解Tension-p shear-p tension-n shear-n表示的含义？

答：在计算的循环里面，每个循环中，每个zone都依据failure criterion处于不同的状态，shear 和tension分别表示因受剪和受拉而处于塑性状态.n表示now，p表示previous，即分别表示在现在和以前的循环当中处于塑性状态.

29. 绘制相关曲线命令的一些解释

1) 在plot hist m vs n的形式里，m代表y轴，n代表x轴(不管m，n的正负); 2) \表示对其值作\，如果大家使用AUTOCAD的话就明白\的含义了，比如上面的 hist 1 gp ydisp的值是0~100，那么vs -1就变成-100~0.以此类推。

3) 为什么要这么做呢?主要是为了符合岩土工程的习惯需要，例如桩载荷试验曲线，沉降曲线等都不是画在第一象限内.

30. 如何查看剖面上任一点的位移值？

8

答：print gp disp range x() y( ) z( )或print gp disp range id plot gp disp range x() y( ) z( )或plot gp disp range id 31. fix x y z range x -0.1 0.1 ..........................a

fix x range x -0.1 0.1................................b 请问a式与b式有何区别？

a表示在x=0这个平面是固定的

b表示在x=0这个平面只x方向不能发生运动 32. initial 与 apply 有何区别？

initial初始化命令，如初始化计算体的应力状态等； apply边界条件限制命令，如施加边界的力、位移等约束等。

initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化。

33. 如何查看各个时段不平衡力的具体数值？ 采用his来记录计算，包括位移应力等 命令his unbal his gp(zone) zdis range (0 0 0) 或者id=？

导出数据命令his write n vs m begin 时步 end 时步 file filename.his n表示纪录的id m表示时步 要导出不平衡力的具体数值 his unbal

step 100000 or solve

his write 1 vs step begin 1 end 1000 file 123.his

使用上述命令就可以查看各个时步下的不平衡力的具体数值 34. 如何显示变形轮廓线的命令？ plo ske magf 10 其中10为放大系数

35． 判断模型中网格点是一个还是两个的方法

显示模型网格点的ID，如果有两个ID号，说明是两个。节理面上的网格点应该有两个ID号

36． 改变模型应力边界时FLAC与3DEC的不同之处

在3DEC中，当改变应力边界时，模型表面的应力直接被更换，FLAC中改

9

变应力边界时，表面的应力为两者应力相加。

37．对DD,DIP的理解

DD表示节理面正方向在XY面上的投影与Y轴的夹角，也等于节理面与XY面的交线与X轴的夹角；DIP表示节理面与XY面的夹角。

38．塑性状态的解释

-p 表示某一区域的应力在模型运行过程中进入过屈服状态，但现在已经推出了屈服状态；

－n 表示某一区域正处于屈服状态 39．对网格点速度的理解

flac中网格点的速度表示的是单位时步网格点发生的位移，查看网格点的速度有两种方法：

（1）绘制整个区域的速度图 （2）记录某一网格点的速度 40．flac中影响运行时间的因素

（1）单元材料，结构面之间的刚度对比很大；

（2）划分单元尺寸对比很大，尽量避免长细比大于5：1的细长单元。 41．阻尼的选择

（1）对于静力分析，局部阻尼比较有效，是默认状态的(SET MECHANIC DANM LOCAL)。

（2）组合阻尼，比局部阻尼更容易在大量网格点发生刚体运动时更好的收敛到稳定状态，使用SET MECH DAMP COMBINE来调用组合阻尼，组合阻尼对于减小动能方面不如局部阻尼有效。

42．FLAC中有多少中材料模型 内置11中本构模型

（1）零模型 用于开挖回填的区域 （2）各向同性弹性模型 （3）正交各向异性弹性模型 （4）横观各向同性弹性模型 （5）德鲁克－普拉格塑性模型

10

（6）莫尔－库仑塑性模型 一般岩土体的力学行为

（7）节理化塑性模型 用于研究材料破坏后行为（如边破稳定，屈服煤矿开采）

（8）应变硬化/软化模型－库仑塑性模型 薄板层状材料破坏后力学行为 （9）双线性应变硬化/软化节理化塑性模型 临近层状土层开挖 （10）双屈服塑性模型 （11）修正的剑桥粘土模型

43．初始化应力时为什么用了ini命令后还用apply命令

ini命令给模型所有区域应力赋予改定的值，但是它不能保证应力是平衡的。这导致了至少存在两个潜在的问题，第一应力可能违反该区域设定的非线性模型屈服准则，step命令运行后，模型会出现塑性流动。第二，网格边界的指定应力可能和假定的初始应力不相等，当step命令执行后网格点可能移动。为保证应力确实达到初始状态，必须使用ini和apply命令。

44．对不平衡力的理解

flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小，该值为百分数，且从不会减小到零，因此，模型不会达到绝对稳定的平横状态，可以根据要求的精度设定R值，模型默认的R值是1e-5。注意R值小仅表示所有网格点的力都平衡，但模型还可能发生塑性流动，通过 plot vel 命令查看仿网格点速度，一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。

45

plot add cont disp out on 时时显示云图 plot block group gpnum on 显示节点编号 plot inter red ske 看接触面的效果 del range group 1 not 不删除group信息 del range group 2 删除group信息

ini z add 数值大小 range group name 顺z方向移动

dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线 ，dd表示y轴正向顺时针到那条射线的夹角。首先应该

11

按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．如gen zone reflect dip 90 dd 90是以YZ为对称面，gen zone reflect dip 90 dd 0 是以xz为对称面

建立接触面 interface 1 face range x(范围）y(范围） z(范围）或者 interface 1 face range plane orig(x y z 坐标） norm（法线方向）

water table clear 清除水位线 pl set back white ;设置背景

set plot bitmap size (800,510) ;设置图片大小

plot set caption size 37 ; 设置图例数字大小 (default=35 范围10-50) plot set center 100 0 40 fish命令的解释（网上摘录） def aux1 定义函数

loop n(1,nptab1) 开始一个循环，从1循环到变量nptab1 rr=xtable(1,n) 从table 1 的x列中读第n个数付给变量rr zz=ytable(1,n) 从table 1 的y列中读第n个数付给变量zz loop k(1,nprof) 开始第二个循环，从1循环到变量nprof case_of k 下面是一个选择结构，利用变量k判断 alfa=0.5*(pi+.1)*float(k-2)/float(nprof-2) xtable(n+10,k)=rr*cos(alfa) ytable(n+10,k)=rr*sin(alfa)

case 1 当k＝1的时候执行下面语句

xtable(n+10,k)=rr table n+10的x列第k个值等于变量rr ytable(n+10,k)=-40.0 table n+10的y列第k个值等于变量-40 case 2 当k＝2的时候执行下面语句

xtable(n+10,k)=rr table n+10的x列第k个值等于变量rr ytable(n+10,k)=0.0 table n+10的y列第k个值等于变量0 end_case 结束选择结构 end_loop 结束第二个循环

12

end_loop 结束第一个循环 end 函数定义结束 aux1 执行什么定义的函数

fix x y z range x -0.1 0.1 表示在x=0 这个平面是固定的 fixed boundary fix x range x -0.1 0.1表示在x=0 这个平面只x方向不能发生运动roller boundary

46．在固定面上施加位移时，fix与apply不能同时使用，应该用ini和fix,才能施加上。

47．菜单 plotitems/add/structure elements/geometry/结构单元的形式，把对话框里的mag fac的值改成不同的数，就可以看到位移放大后的结构单元。

48．你可以用键盘上的上、下、左、右键在窗口中移动模型。 49.unbal

1)不平衡力在宏观上，是由于外力和内力的差值引起的，这是本质上的来源，外力来源可以是荷载、体积力、开挖引起的应力释放、温度应力、固结等。

2)在本构模型上，由于弹塑性计算中的应力“估算”，都是按照某个应力状态下的”弹性模量“下的“弹性应力”，经常发生高估实际应力状态的情况，这时候，将土体所能承受的应力（由本构模型决定）和计算所得的弹性应力的差值也可以理解为不平衡力。这些不平衡力随着，模型的塑性流动，屈服面的向外扩张，以能量的形式消耗调了，也就是说，由于系统的自适应，不平衡力是可以“消化掉”的。fem的非线性计算也是采用这套思路。

3)对于不平衡力在物理意义上的真实性，我的看法是：如果阻尼等参数符合土体的真实属性（对于流变性土体可能就容易接受一点），不平衡力的传递过程极其消散过程是可以理解为真实的过程。如果参数的选取只是为了迅速达到quasi－static状态，那么就没有真实性可言。

50.节点信息

pri gp pos输出节点信息 pri zone gp输出单元的节点编号 51.attach与merge命令

gen merge 是针对模型中的节点运行的，在模型中坐标容差在1e5以内的相

13

邻节点会合并在一起，而attach face 的作用是把容差范围内的相邻面粘在一起

顺序要小心，先gen merge后attach face，从以上两条命令作用的级别也可以看出，应该先合并点，后粘面，否则会发生系统致命错误而退出计算

14

4. 快捷键命令

1. ctrl+G：有彩色变黑白，有色曲线变成有区别的点画线 2. ctrl+R：变回原来的状态 3. ctrl+c：brings up the camera dialog 4. ctrl+l: cuttingplane

5. 双击鼠标左键，可获得单元的 id 和 position

6. 如果不晓得一个命令后面跟什么关键字，可以在命令后加 ？ 可以获得一点

help！比如键入 plot ？ 然后回车，看会有什么提示 7. 直接按delete,图形就缩小,M是图形放大

8. 如果不想按shift+m键，可以先按一下caps lock键，然后再按m键就可以实

现图形缩小的功能

9. 放大和缩小的操作可以分别通过Insert和Delete来实现\\

15

5. 边坡稳定性计算

1. FLAC3D做了个边坡模型，用step命令计算的最大位移是20cm，并且ubal

最大不平力也接近零了，最大位移也基本稳定了，应该可以认为达到平衡了。再用solve fos 计算安全系数时，fs＝1.02，可以认为在临界状态，可是现在最大位移达到50cm？

? 用solve fos计算的时候，强度会折减，位移肯定会变大。 ? zdis其实不是实际的，是不可用的位移

? slove fos 除了得到安全系数以外，其它没有什么用处。它的位移场和速

度场的结果都是不能用的，因为这是采用折减后的参数算出来的结果，肯定要比采用未折减的参数算出来的结果大一些。所以千万不要用这个结果去进行应力－应变分析，会得出错误结论的。

2. 计算边坡安全系数时候，出现

search for fos terminated prematurely fos is larger than 64.00 ? 安全系数过大 停止计算 3. 怎么用位移来判断边坡破坏了呢？

? 如果用位移判据，当位移突变时即认为边坡失稳了。

16

6. 关于结构单元之间的链接问题体会和心得

http://www.simwe.com/forum/thread-756746-1-1.html

1. 结构单元建模中node的关系处理问题心得与总结：

相同结构单元之间的节点的链接问题：

相同单元类型，节点（ID）编号相同，那么不会再产生新的node，而在该位置处公用原来单元的node。

比如在隧道掘进的过程中，要涉及到前、后开挖步过程中喷砼之间的纵向连接，那么在Flac3D 2.1以上版本中，只要所前后开挖步中产生的shell单元是同一ID，那么在相同位置处的node会被自动链接，这里的链接不仅仅是几何上的连接，也会使内力在该shell中传递（轴力，弯矩）。

（1）

gen zone brick p0 0 0 0 p1 20 0 0 p2 0 10 0 p3 0 0 10 sel shell id 1 range x 0 10 z 0 10 y -0.1 0.1 sel shell id 1 range x 10 20 z 0 10 y -0.1 0.1

相同ID，不再产生新的node，几何链接，内力会传递。 （2）

gen zone brick p0 0 0 0 p1 20 0 0 p2 0 10 0 p3 0 0 10 sel shell id 1 range x 0 10 z 0 10 y -0.1 0.1 sel shell id 2 range x 10 20 z 0 10 y -0.1 0.1

不同ID，连接处会产生新的node，几何不链接，内力不传递。 “类似的，对于cable，beam，pile等来说，这个法则同样适用。” 2. 不同结构单元之间的节点的链接问题：

以上是对于相同类型的结构单元来说的。对于不同结构单元来说，应用时分两种情况。

一种是共用节点。比如说：先生成一个shell，然后生成一个beam，其中beam的一个端点是shell单元中的一个node。该node就是共用节点。该节点同时具有shell单元和beam单元的属性，力学机理比较复杂。

另外一种是通过产生结构单元之间的link来定义结构单元之间的结构关系

17

(sel link id link_id struct1_id target node tgt_num struct2_id)。

（1）一般立柱与梁板的关系——铰接。

sel link attach xdir=rigid ydir=rigid zdir=rigid range id link_id sel link attach xrdir=free yrdir=free zrdir=free range id link_id （2）悬臂梁端部——固接。

sel link attach xdir=rigid ydir=rigid zdir=rigid range id link_id sel link attach xrdir=rigid yrdir=rigid zrdir=rigid range id link_id 或者其他转动或移动限制组合方式。

具体实例参看：http://www.simwe.com/forum/thread-726231-1-1.html& 里面有比较详尽的fish命令流。

3. 隧道初期支护中的“钢拱架”应该如何建立

http://www.simwe.com/forum/thread-726231-1-1.html 1、钢拱架一般都采用多个beamsel小段来模拟（调整形状），同一榀钢拱架只要保证组成它的beamsel都属同一id就可以了。这一部分好解决。

2、设置每个node与网格的关系。我个人认为：实际的钢拱架都与围岩密贴接触（加垫块调整），这和shell模拟的喷混凝土与围岩的接触关系有异曲同工之处，区别仅在于网格不能带动node（钢拱架）旋转。因此node与zone之间的Link就应该只保留x,y,z，而将link中的xrot,yrot,zrot删除。

3、如果钢拱架安装之前，已经有shell（喷混凝土）附在网格上。那么在生成beamsel时，两个端点node可以与shell的node设在同一位置（两个node，分属于shell和beam），然后删除node（beam）与zone的Link，最后建立beam node与shell node 的联结（仅是x,y,z,不包括xrot,yrot,zrot），这保证了beam和shell的实际联结情况——密贴但不会随之转动。

按照这个思路我试着建立了模型，效果还是不错的.我把link的fish处理程序（一个处理Link的小思路）贴出来吧，大家分享一下，可能有的地方不是很合理，欢迎斧正。

restore 111.sav& M: def parmater ;起始位置 z_=0

;榀数（纵向）

18

p=5 ;半径- Rbeam=5.850006 Rshell=6.0084(洞室半径) ;首个node的角度位置 a=02 det(a)=11.0589 end: -

Parmater

;beam半径稍小于shell（洞室）半径 def operate_link) link_id=1000% z=z_ loop m(1,p) loop n(0,13)

x_shellnode=Rshell*cos((a+n*det(a))*degrad) y_shellnode=Rshell*sin((a+n*det(a))*degrad) z_shellnode=z. _( x_beamnode=Rbeam*cos((a+n*det(a))*degrad) y_beamnode=Rbeam*sin((a+n*det(a))*degrad) z_beamnode=z

sn_pointer=nd_near(x_shellnode,y_shellnode,z_shellnode)' bn_pointer=nd_near(x_beamnode,y_beamnode,z_beamnode)

shellnode_id=nd_id(sn_pointer) beamnode_id=nd_id(bn_pointer)

command

sel set link node_tol 0.3 sel

link id link_id beamnode_id target node tgt_num shellnode_id sel link

attach xdir=rigid ydir=rigid zdir=rigid range id link_id* sel link

attach xrdir=free yrdir=free zrdir=free range id link_id0 end_command + l$ n5 e9 _+ k; V) Q+ g& B8 link_id=link_id+1 end_loop z=z+0.8 end_loop

19

end operate_link

sel set link node_tol 0.3 表示只要两个node之间距离小于0.3，就产生link x_shellnode=Rshell*cos((a+n*det(a))*degrad) 这是给出和beam的node产生link的node（shell）坐标值

sel link id link_id beamnode_id target node tgt_num shellnode_id 是在将beam的node 和shell的node link起来。手册中有这个命令行的解释。搜索Link章节即可查得。

sel link attach xdir=rigid ydir=rigid zdir=rigid range id link_id 是设置link的初始条件。手册中亦有这个命令行的解释。搜索Link章节即可查得。

你好！你提出的建立shell 与beem建立连接整体受力的方案很好。 但你帖出的“link的fish处理程序”中没有删除原来的shellnode或beemnode与zoom之间的连接（link ），而直接创建beamnode与shellnode之间的连接能行吗？一个节点已存在某种连接时不能再和其它节点创建连接吧。

产生beamnode时，我控制了beam产生的半径略小于洞室半径，这样beam不会和zone产生Link。另外，shellnode与zone之间的link不用删除。因为对于一个source node(shellnode)来说，它只能拥有一个属于它自己的Link（shellnode和zone之间的Link，对于这个Link来说，shellnode是source node），也就是说这个shellnode不能再去\其他的zone或者node。但是这个shellnode却可以被多个其他的zone或node去“link”，在上面的程序中，我设置beamnode去\。

20

7. 模拟预应力锚索

方式1、通过删除－建立link模拟。

删除掉锚索端头原来CABLE和ZONE之间自动建立的连接，然后建立钢性连接。示例命令如下：

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10 sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17 sel delete link range id 1 sel link id=100 1 target zone100

&S(} ~6?0s;_;O仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent\0g9K:f6R$D9_-K$C4j

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id sel cable id=1 pretension 60e3 range cid 1,10

SimWe仿真论坛*b:x)T7f6m,y2C方式2、通过赋三段属性模拟。即将锚索的端头、自由段、锚固段赋不同的属性来模拟预应力锚索。示例命令如下：

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 2,10

5e-y6i(i7Q:])G#Qwww.simwe.com8b2L)g9l3`.t&E N/^2p d

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17 gr_coh 10e8 gr_k 2e10 range cid 1,1 sel cable id=1 pretension 60e3 range cid 1,10

方式3、借助别的结构单元（如liner单元）来模拟托盘。删除掉锚索端头的link，然后建立新的link，新link的target为liner上的node。示例命令如下： sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent3A5I/y%k9Q0F5t9Q7]1d1z+B5swww.simwe.com9f9l:J F){.a'g8k

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM.I5?2u$U5B5Jsel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10 gr_coh 10e5 gr_k 2e7 range cid 11,17 sel liner range y=-.1, .1 x=-1,1 z=-1,1

sel liner PROP iso=( 25e9, 0.15) thick=0.1 ; concrete

(A$j9P2d\n2f|Simwe.com仿真设计有限元虚拟仪器+v9Y'Y1{'t,|sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

sel liner PROP cs_nk=8e8 cs_sk=8e8 &

cs_ncut=0.0 cs_scoh=0.0 cs_scohres=0.0 cs_sfric=0.0

$C'} ~/R\sel delete link range id 1

www.simwe.com(\\#|2T.I9u-c u

&Z$w!E(x/H L(Qsel link id=100 1 target node tgt_num 18

21

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100

|Simwe.com|仿真设计有限元虚拟仪器-O7{1S'o#z-y(c*_

sel cable id=1 pretension 60e3 range cid 1,10 总结

1、 预应力和全锚均通过删除－建立LINK模拟。

2、 预应力通过sel cable id=1 pretension 60e3 range cid 1,10方式施加 3、 预应力加在自由端。 4、 全锚预应力加全根

锚索属性讨论

最后想提一点，预应力锚索的几个属性参数非常重要emod ytension xcarea gr_coh gr_k 模拟的时候务必搞清楚几个参数的真正含义！！许多在公开刊物（包括 岩石力学与工程学报）发表的论文，里面的参数都不恰当，有一篇模拟动载对锚杆影响的论文（不方便具体指出），给出的粘结强度竟然为100e5！这是什么概念，每米锚固剂内能承受的最大剪切力为1000吨力！我用这样的锚杆支护井下巷道的话，只怕再大的采动影响，巷道还是稳稳的不动。另外还有一点，另一个参数gr_fric，锚固剂的摩擦角，指定锚固剂内所能承受的最大剪力随着围压的增加而增加，个人觉得没有给出的必要。我查过许多文献，目前还没有查到锚固剂剪切力会随围压增加而增大的论证。也许是本人水平有限，请高手不吝指教。谢谢。

再说预应力施加位置

前面说的预应力加在自由段合适，现在有疑问：实际工况中，在对锚杆进行拉拔加载预应力时，锚固剂肯定是受力的（岩体通过锚固剂给锚杆剪力以平衡拉拔力），那么在模拟的时候，预应力是不是不能只加载到自由段，也要加载到锚固段？如果锚固段也加的话，肯定不是全锚固段都加，而是加靠近自由段的一部分，那么，应该加多长合适？？锚杆锚固机理，上面的问题就必须要考虑了。 模型

很简单，就是一块长方体，不设任何应力，只算附加应力。如下：

www.simwe.com;y4s7U7J&B5Z0h7S,i _

其实如果只模拟巷道支护，不需要考虑这么细，加在自由段就行了。但如果要分析预应力

gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 & ratio 1 1 1 1.1 fill

gen zone reflect normal 1 0 0 ori 0 0 0 gen zone reflect normal 0 0 1 ori 0 0 0 mo mohr

pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000

fix x range x -25.1,-24.9 fix x range x 24.9 25.1 fix y range y 49.9 50.1 fix z range z -25.1 -24.9

22

fix z range z 24.9 25.1

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

gr_coh 1 gr_k 1 gr_per 0.0785 range cid 2,10

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e8 gr_k 2e10 range cid 1,1

sel cable id=1 pretension 60e3 range cid 1,10

step 2000

没做过锚杆的,有两个问题请教tg0215老大

虚拟仪器0n)|:M%Y;D*K3x/K'{SimWe仿真论坛4N'y+K5R.V,g.`

|Simwe.com仿真设计有限元1.对于形成锚杆的三个控制点(0,0,0),(0,29,0 ) ,(0,35,0),是不是需要与模型的节点重合? 1、不需要

2、不需要，但要注意端头node位于最靠近开挖面的zone上。FLAC3D里是这样的情况，2D里另有别论。

怎么显示link及其id号？ pl sel geom link on id on range

|Simwe.com|仿真设计有限元虚拟仪器'q*}6\\6w \\2.锚杆的起点(端头)是不是必须要位于模型表面?因为在一些不规则的面上,很难精确得 ...

gr_per是指孔径周长，也即锚固圈周长，我没有给出是因为我认为锚固剂的剪切强度跟围压没有关系，gr_per是用来计算sbond的。gr_per是用来求sbond的，如果给出了gr_per，又给出了fric，则sbond会随有效围压的增大而增大。但是现在没有文献证明锚固剂的最大抗剪强度会随有效围压增大而增大（我没有看到过，如果你有这方面资料，希望能共享一下）。实际情况中，砂浆厚度肯定影响锚固力，数值模拟中通过sbond和kbond来描述砂浆的锚固力，并不表明砂浆厚度对锚固力没影响。kbond是锚固剂刚度，sbond是锚固剂强度，注意区分。两个参数一半通过实验得到或用公式计算，具体见手册。

对《受邀谈论FLAC3D预应力锚杆模拟心得》帖子的一点补充 sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10 sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17

23

delete link range id 1 sel link id=100 1 target zone

SimWe仿真论坛+P&?2L-?2y

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent#K\&M9S.H&Wsel cable id=1 pretension 60e3 range cid 1,10100

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent%|/^5C2S;Y在这个锚杆支护例子中，如果工程只需要几个锚杆还方便处理，如果锚杆数量多了后就不好处理了。

www.simwe.com#m {#~\

建议：把sel link id=100 1 target zone这段命名中的100改为1，因为delete link range id 1已经将代号为1的连接去除了，这样就可以再生成一个代号为1的连接，同时还要修改

限元模拟计算力学航空航仿真分析有sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100

仿真分析有限元中的100为1。

|Simwe.com|仿真设计有限元虚拟仪器-k,Q#e3Q*r1C

这样修改的理由是：

1、工程中锚杆是随着开挖而施加的，那么就是说锚杆不会一次性全部加上，而做

天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent)o4J#y4`!C\仿真分析有限元模拟计算力学航空航连接的时候不可能一次全部连接。

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent:C.A8f'V!X4O/@0k*e

2、如果分开加锚杆，那么先用id＝100生成了新的link，之后加载的锚杆的link号就会 从100开始增加，这样在再此需要修改link时，就不方便寻找到link号，也就是

!L4h!H+~\delete link range id 1 sel link id=100 1 target zone两个语句中的id 1。起始号。 例如： ;开挖第一步

mo null range group 1 x -14.6 0 z 9.3 12.1

www.simwe.com$]6x!}6y#N,s:}%z8y1L|Simwe.com|仿真设计有限元虚拟仪器*d!V+b#Z7k%@-y仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM*G'}4E*~/`&[&N

3、如果把id＝100改成1 ，那么之后重新修改link时，需要修改的link号还是这个锚杆的

mo null range group 2 x -14.6 0 z 9.3 12.1

sel cable id=1 begin 0 0.5 10.6 end 12.48 0.5 7.22 nseg 10 sel cable id=1 begin 0 0.5 9.30 end 12.48 0.5 5.92 nseg 10range cid 1,20

sel delete link rang id=1 sel delete link rang id=12 sel link id=1 1 target zone sel link id=12 12 target zone

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM/L,q&O6t1e&\\/S2s\仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM8c8@,V2f-V2v.P4^仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM0E4f,E!r+Q0j5F3^(_

sel cable pro emod 200e9 xcarea 0.00126 gr_per 1 yten 4e6 gr_k 6e9 gr_c 3e9 gr_fric 30

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 1

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 12

www.simwe.com/T4h#J1}1X!X8E

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM*u!a+z!B'z'B8Q M2@7l;开挖第二步

mo null range group 2 x -14.6 0 z 8 9.3

sel cable id=1 begin 0 0.5 8 end 12.48 0.5 4.62 nseg 10

/c6|7J![&F4`sel cable pro emod 200e9 xcarea 0.00126 gr_per 1 yten 4e6 gr_k 6e9 gr_c 3e9 gr_fric 30 range cid 21,30

24

sel delete link rang id=23

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent9|+I,l$A#J+L

sel link id=23 23 target zone

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM&x!J(v0C9r J9~*~.c%S4Q www.simwe.com2X&I4V5m3f+~0D

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 23

.b\7W2w;v\

最后要说明的是：我也不能肯定我说的一定对，如果错误希望大家纠正。

模拟锚杆在flac里可以用pilesel或者cablesel，前者适用模拟锚杆加固作用例如土体中锚杆，后者适用锚杆悬吊作用，例如模拟锚杆悬吊隧道顶板楔形体。不加选择使用cablesel必然导致上面说的结果。

关于这个问题，讨论了不少了，俺也来发个言，仅供菜鸟参考，请大牛指正！ 俺认为，这个问题首先要从预应力锚杆这种工程措施开始：

3 S. g& ~+ v$ Z. N) B$ g+ W9

1、破坏形式，有三种破坏形式，a,钢筋的屈服破坏.b,锚固体（砂浆）破坏.c,锚固体和钢筋以及围岩的接触面破坏。对于这三类破坏，软件都使用理想弹塑性模型。对于A类破坏，使用屈服强度ytens和emod来描述，如果你知道锚杆不会屈服，可以设置一个较大的ytens来保证钢筋不会屈服，当然如果你使用一个很小的ytens，那么，就可能出现塑性流动。对于B类破坏，是指锚固砂浆被剪破，发生单位位移时每米锚杆能承担的力，这里一定要记住是理想弹塑性模型啊！俺的理解，这种破坏应该在锚固体内。软件中这个指标叫gr_k，只会与锚固砂浆的剪切模量和锚固体的截面几何形状有关，gr_k=2*Pi*G/10*ln(1+2t/D)，式中G为砂浆的剪切模量，t为截面砂浆厚度，D为锚杆钢筋的直径。对于C类破坏，是发生在锚固体的两个接触面上的，但你只能定义一种。首先，你必须明确，既然是接触类型的问题，就应该有C和Phi，这和边坡是类似的。这两个参数的确定，是应该按抗拔试验来确定的，将试验锚杆的最大抗拔力/锚固长度作为y轴，围岩压力*锚固体直径为x轴，截距就是C，软件中是gr_coh，倾角就是gr_fric。实际应用中，软件可以采用gr_coh等于D（钢筋拔出）或者D+2t（锚固体拔出）乘于pi再乘于max（围岩和锚固砂浆的单轴抗压强度）/2来计算；如果不考虑围压的影响，gr_fric可以设置为0。

2、预应力锚杆的布置形式，一根预应力锚杆可以由许多段的sel cable组成，只要begin和end相互连接就可以了，这些点也不一定在单元的节点上，而且可以不是一条直线，自由段的模拟俺认为可以将gr_coh、gr_fric、gr_k设置为0，但gr_per应该和锚固段一致。

3、预应力的施加，俺认为使用sel cable pretension较好，直接加在锚固段上，锚头嘛，你应该和梁连在一起，呵呵，俺针对的是边坡啊

锚杆和节点不相连计算的结果相差很大。我没试过锚杆的，但是桩的单元我试过，桩体头尾和节点相连的计算结果和一段一段建立桩单元计算的结果差非常大。但是一段一段建立桩单元计算的结果和头尾和节点相连并在桩土接触之间建立接触面单元的计算结果很相近。有没有人有兴趣讨论一下？

25

8. 其他

set large在什么情况下使用？

大变形一般指的是大应变，在应变为5-10%的时候一般称为大应变，需要采用large mode。在实际操作中，对于基坑开挖、软土变形、洞室等工程一般都要用大应变模式来求解。FLAC的大应变也不能计算非常大的应变问题，也只能相对小变性下的应变，当然还与你的网格划分有关系。变形太大时，网格畸形，会停止计算。在进行边坡fos计算时，在进行sol fos时候，若是设置了大变形，有的时候网格变形太大，会导致sol fos计算终止。

E－－弹性模量 Es－－压缩模量 Eo－－变形模量 弹性模量＝应力/弹性应变，它主要用于计算瞬时沉降；压缩模量和变形模量均＝应力/总应变，压缩模量是通过现场取原状土进行实验室有侧限压缩实验得出的，而变形模量则是通过现场的原位载荷试验得出的，它是无侧限的。弹性模量要远大于压缩模量和变形模量，而压缩模量又大于变形模量。

地堪报告中，一般给出的是土的压缩模量Es与变形模量Eo，而一般不会给出弹性模量E。

按规范的规定，在地基变形验算中要用的是压缩模量Es，但因Es是通过现场取原状土进行试验的，这对于粘性土来说很容易做到，但对于一些砂土和砾石土等粘聚力较小的土来说，取原状土是很困难的，很容易散掉，因此对砂土的砾石土通常都是通过现场载荷试验得到Eo，所以在地堪报告上，对于砂土的砾石土一般都仅给出Eo，即使给出Es，也是根据Eo换算来的，而不是试验直接得出的。理论上Es和Eo有一定的关系，但根据该关系换算误差较大，所以二者关系一般都根据地区经验进行换算。弹性模量一般可取为压缩模量的3～5倍，上海地区经验一般为2.5～3.5倍（见同济大学 杨敏教授相关论文）数值分析时可以适当加大一些。

请问如何能让模型的变形显示的更明显一些？ magfac (2D的情况下用: plot plas no_past）

26

初始应力不仅仅是重力引起的啊！它还包括土压力，水浮力，上方建筑物的力，这些信息的勘测一般都不容易获得。通常可以把模型材料设为弹性，再把两个模量设到很大solve后再把模量和材料还原可得。设置重力的原因不是为了平衡减少程序计算步骤，而是因为如果未初始化应力的话，在set grav的时候，模型会在自重应力下发生应力变形。所以，初始化应力的原因是防止模型在自重应力下发生变形（我们认为地层自重变形已经发生，当然欠固结等现象是另外一回事）。ini命令来对地层进行初始应力的施加，而APPLY是施加应力边界条件

27

本文来源：https://www.bwwdw.com/article/m7oa.html