下载文档原格式
gen zone brick\wedge\pyramid\tetrahed size 6 8 8\\获得模型plotcreate Trench\\创建块Trenchadd surface yellow\\给块加色add axes black\\建立坐标系show\\显示create Trench2add surface yellowadd axes blackset rotation 17.495 9.847 28.481\\旋转set center 3 4 4set dist 26.12set magnification 0.8\\于原始模型的比例show\\显示q\\退出model mohr\\应用摩尔库仑准则print mem\\调入内存prop bulk = 1e8 shear = 0.3e8 fric = 35bulk\\给体弹性模量、剪切弹性模数、内摩摩擦角赋值prop coh = 1e10 tens = 1e10\\给粘聚力、抗张强度赋值set grav 0, 0, -9.81\\抗张强度设定重力加速度ini dens = 1000\\密度fix 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.1\\设定xyz的边界范围hist n = 5\\设定每计算五布显示一次hist unbal\\显示最大失稳压力记录状况hist gp zdisp 4,4,8\\显示z方向的位移在点448记录状况set mech force=50\\限定失稳压力小于50时停止计算solve\\计算plotplot create GravV\\创建gravv体plot set plane dip=90 dd=0 origin=3,4,0\\建立面与xy面垂直并垂直于y轴并以3,4,0点为起始建立gravv体plot set rot 15 0 20 ; this could be achieved interactively\\旋转plot set center 2.5 4.2 4.0 ; and is shown here simply to illustrate \\中心plot add bound behindplot add bcont szz planeplot add axesplot showqsave trench.sav\\保存plot printplot print view\\显示视图prop coh=1e3 tens=1e3\\给整个块体的粘聚力和抗拉强度赋值model null range x=2,4 y=2,6 z=5,10\\确定洞的尺寸set large\\设定最大合理拉伸力使其破坏ini xdis=0 ydis=0 zdis=0\\为计算方便与结果无关step 2000\\设置时间段步骤plotcreate DispContcopy GravV DispCont settingsadd cont disp plane behindadd axesshowprint gp position\\显示ID号及节点数和其坐标print zone\\显示块数及质心坐标1.FLAC3D是有限元程序吗？答：不是！是有限差分法。
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的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
《FLACFLAC3D基础与工程实例》阅读札记目录一、FLACFLAC3D软件概述 (2)1. 软件背景与简介 (3)1.1 FLACFLAC3D的发展历程 (4)1.2 软件的应用领域及特点 (5)2. 软件安装与运行环境 (6)2.1 系统要求 (7)2.2 安装步骤 (8)2.3 运行环境配置 (10)二、FLACFLAC3D基础知识 (11)1. 基本概念与术语 (13)1.1 有限元分析原理 (14)1.2 离散元法简介 (14)1.3 FLACFLAC3D中的相关术语解释 (15)2. 软件操作界面及功能模块 (17)2.1 操作界面介绍 (18)2.2 主要功能模块说明 (20)2.3 菜单功能详解 (20)三、工程实例分析 (22)1. 地质工程实例 (23)1.1 工程背景及问题定义 (25)1.2 模型建立与参数设置 (26)1.3 结果分析与讨论 (27)2. 土木工程实例 (29)2.1 工程概况与建模目的 (30)2.2 建模过程及计算步骤 (31)2.3 结果展示与工程应用 (32)四、FLACFLAC3D应用技巧与注意事项 (33)1. 建模技巧与优化方法 (34)1.1 建模策略及优化思路 (35)1.2 网格划分与模型简化技巧 (36)1.3 参数设置与模型验证方法 (38)2. 数据分析与处理方法 (40)2.1 数据采集与整理方法 (41)2.2 结果分析与图表展示技巧 (42)一、FLACFLAC3D软件概述3D是一种广泛使用的岩土力学与有限元分析软件。
它是一套专门用来分析连续介质中的物理力学现象的强大工具，主要应用于土木、矿山、隧道等领域，能针对各种复杂的工程问题进行数值建模和模拟分析。
3D以其高效、灵活的数值分析能力，为工程师提供了强大的技术支持。
其主要特点包括：多功能：3D能够模拟多种物理过程，包括应力分析、稳定性分析、流体流动分析等，适用于多种工程场景。
F L A C D基础知识介绍 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序以前的为DOS版本，版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
flac3D蠕变基础知识 蠕 变 模 型 将flac3d的蠕变分析option进行了简单的翻译，目的是为了搞清楚蠕变过程中系统时间是如何跟真实时间对应的。
1. 简介 Flac3d可以模拟材料的蠕变特性，即时间依赖性，flac3d2.1提供6种蠕变模型： 1. 经典粘弹型模型 model viscous 2. model burger 3. model power 4. model wipp 5. model cvisc 6. powe蠕变模型结合M-C模型产生cpow蠕变模型（model cpow） 7. 然后WIPP蠕变模型结合D-P模型产生Pwipp蠕变模型（model pwipp）； 8 model cwipp 以上模型越往下越复杂，第一个模型使用经典的maxwell蠕变公式，第二个模型使用经典的burger蠕变公式，第三个模型主要用于采矿及地下工程，第四个模型一般用于核废料地下隔离的热力学分析，第五个模型是第二个模型的M-C扩展，第六个模型是第三个模型的M-C扩展，第七个模型是第四个模型的D-P扩展，第八个模型也是第四个模型的一种变化形式，只是包含了压硬和剪缩行为。
2. flac3d解流变问题 2.1简介 流变模型和flac3d其他模型最大的不同在于模拟过程中时间概念的不同，对于蠕变，求解时间和时间步代表着真实的时间，而一般模型的静力分析中，时间步是一个人为数量，仅仅作为计算从迭代到稳态的一种手段来使用。
2.2 flac3d的蠕变时间步长 对于蠕变等时间依赖性问题，flac3d容许用户自定义一个时间步长，这个时间步长的默认值为零，那么材料对于粘弹性模型表现为线弹性，对于粘塑性模型表现为弹塑性。（命令set creep off也可以用来停止蠕变计算。）这可以用来在系统达到平衡后再开始新的蠕变计算。蠕变公式中包含时间，所以计算中时间步长对程序响应有影响。 虽然用户可以对时间步进行设置，但并不是任意的。 蠕变过程由偏应力状态控制，从数值计算的精度来讲，最大蠕变时间步长可以表示成材料粘性常数和剪切模量的比值： For the power law ----------省略。 For the WIPP law -----------省略 For the cvisc model, 上面方程应该写成： tmax = min ( ηK/GK，ηM/GM) 上标K和M分别代表Kelvin和Maxwell。 蠕变压缩的时间限制包括系统体积反应，并且估计为粘性和体积模量的比值。粘性可以表示为σ和体积蠕变压缩速率的比值。 建议利用FLAC3D作蠕变分析开始时所采用的蠕变时间步，比根据上式算得的时间tmax小两到三个数量级。通过调用SET creep dt auto on ，可以利用自动时间步自动调整。作为一项规则，时间步的最大值（SET creep maxdt ）不能超过tmax。 用来计算tmax的应力σ大小，可由蠕变开始之前的初始应力状态决定。同样，σ作为von Mises不变量，可以用FISH函数计算。 涉及体积变化响应的蠕变分析，其最大时间步长可以表示成材料粘性常数和体积模量的比值，这里粘性常数就是平均应力 和蠕变体应变率 的比值。 一般flac3d推荐使用的初始蠕变时间步长比最大时间步长（由上述公式计算得到的）约小2到3个数量级。如果使用set creep dt auto on命令，那么程序将自动调整蠕变的时间步长，同样应当记住通过命令（set creep maxdt）设置的最大蠕变时间步不能超过 。
2.3自动调整蠕变时间步长 用户可以设置蠕变时间步为一个常数值，也可以使用set creep dt auto on命令自动调节。如果时间步长自动变化，那么当最大不平衡力超过某一阀值时，它就会减小；当最大不平衡力小于某一水平时它就会增大。系统将该阀值定义为最大不平衡力和平均节点力的比值。 对于即将求解的问题，先只考虑弹性效应，通过观察接近初始平衡状态时的不平衡力就可以大概知道整个问题的不平衡力大小。 在有些算例中，又会尽量避免对时间步长的连续调整，这样的话，在某一次时间步长调整之后，可以定义一个“延迟阶段”（比如100步），在这一“延迟阶段”内不再发生进一步的时间步长调整。一般来说，为了适应开挖之类的瞬变，时间步长一开始是一个较小值，然后在模拟过程中增大，如果再发生瞬变，最好人工调小时间步长，然后让它自动增大。 要解决问题的典型的不平衡力准则，可以通过观察只有弹性作用的初始阶段接近平衡时的不平衡力来决定。很多情况下通过逐渐增大或减小时间步可以达到很好的性能。（例如默认比值lmul = 1.01和umul = 0.90).一些情况下，最好避免对时间步进行连续调整，这可能产生“干扰”。为了达到这个目的，时间步发生改变后应有一个“等待时间”（比如：100步），让系统暂停，这段时间中时间步不做进一步调整。通常时间步开始与小的数值来适应短暂的过程比如开挖，然后随模拟进程增加。如果增加一个短暂过程，则需要手工减小时间步然后在让其自动增加。命令SET creep用来设置时间步以及时间步自动调整所需要的参数。
2.4蠕变模型的指令输入 2.4.1 flac3d命令 Config creep 这个命令必须使用，以便开始蠕变分析。 History dt 时间步长 Crtime 蠕变时间 Model viscous 经典粘弹性模型 Print creep 显示蠕变模型的参数信息 Zone prop 显示分配给单元的材料属性 Property bulk（弹性体积模量，K）Shear（弹性剪切模量，G）Viscosity（动力粘性常数， ） Set creep age t 使用solve求解蠕变问题，t是蠕变时间的限值。 Set creep dt ** Set creep mindt=** maxdt=** Set creep Dt Auto on/ Auto off T就是蠕变时间步长，如果不进行赋值，则默认为mindt。Auto on auto off是自动 调整时间步长开关，默认为关；设置为开时，要通过关键字lfob，ufob，lmul，umul和latency来控制时间步长。 Lfob v 不平衡力比率低于v值，则时间步长自动增大。默认为 Ufob v 不平衡力比率超过v值，则时间步长自动减小。默认为5.0* Lmul v 不平衡力比率低于lfob，则时间步长自动乘以v值，lmul必须大于1，默认为 1.01 Umul v 不平衡力比率超过ufob，则时间步长自动乘以v值，lmul必须小于1，默认为 0.90 Latency v 延长阶段的时间步数 Maxdt v 最大蠕变时间步长，默认没有限制，但是设置的时候最好不要超过 Mindt v 最小蠕变时间步长，默认为零 On 蠕变时间开关，一旦确认蠕变分析（config creep）则默认为开 Off Time t 蠕变开始的时间，默认为t Solve age t config creep后，t为蠕变计算的时间限值 FISH中可以用的变量 crdt --------creep timestep crtime-------- creep time
2.4.2 property 经典粘弹性（Maxwell介质）— MODEL viscous (1) bulk 弹性体积模量, K (2) shear 弹性抗剪模量, G (3) viscosity 动态粘滞度, η Burger’s Model — MODEL burger (1) bulk 弹性体积模量, K (2) kshear 开尔文抗剪模量, GK (3) kviscosity 开尔文粘性系数, ηK (4) mshear 麦斯韦尔剪切模量, GM (5) mviscosity 麦斯韦尔粘性系数, ηM Burger-Creep 粘塑性模型 — MODEL cvisc (1) bulk 弹性体积模量, K (2) cohesion 内聚力, c (3) dilation 膨胀角, ψ (4) friction 内摩擦角, φ (5) kshear 开尔文抗剪模量, GK
(6)kviscosity 开尔文粘性系数, ηK (7) mshear 麦斯韦尔剪切模量, GM
(8)tension 极限抗拉强度, σt (9) mviscosity 麦斯韦尔粘性系数, ηM
3. 蠕变的例子－FLAC3D处理二维问题（采用蠕变模型） new gen zone brick size 8 1 22 p0 0 0 -6 p1 4.2 0 -6 p2 0 1 -6 p3 0 0 5 pause ；建立模型 inter 1 face range x 4.1 4.3 y 0 1 z -6 5；加3个接触面 inter 2 face range x 0 4.2 y 0 1 z -5.9 -6.1 inter 3 face range x -.1 .1 y 0 1 z 5 -6 pause； gen zone brick size 30 1 22 p0 4.2 0 -6 p1 50 0 -6 p2 4.2 1 -6 p3 4.2 0 5 ratio 1.05 1 1 pause； gen zone brick size 1 20 10 p0 0,0 0 p1 0 1 0 p2 -20 0 0 p3 0 0 5 ratio 1 1.1 1 pause； gen zone brick size 12 20 1 p0 0 0 0 p1 0 0 -6 p2 -20 0 0 p3 0 1 0 ratio 1 1.1 1 pause； gen zone brick size 30 20 1 p0 4.2 0 -6 p1 50 0 -6 p2 4.2 0 -55 p3 4.2 1 -6 ratio 1.05 1.1 1 pause； gen zone brick size 8 20 1 p0 0 0 -6 p1 4.2 0 -6 p2 0 0 -55 p3 0 1 -6 ratio 1 1.1 1 pause； gen zone brick size 20 20 1 p0 0 0 -6 p1 0 0 -55 p2 -20 0 -6 p3 0 1 -6 ratio 1.1 1.1 1 pause； group soil ；定义组 group exca1 range x -20 0 y 0 1 z 3 5 group exca2 range x -20 0 y 0 1 z 0 3 group mixingpiles range x 0 4.2 z -6 5 ；边界约束 fix x range x -19.9 -20.1 fix x range x 50.1 49.9 fix y fix z range z -55.1 -54.9 ============================================= config creep；启动蠕变模块 def burcal bkk=bme/(3*(1-2*b1)) ‘bkk——k(体积模量);bme——E（弹性模量）；b1——泊松比； bk1=bke/(2*(1+b1)) ‘bke—— bk2=bme/(2*(1+b1)) ‘bk2——G（剪切模量）；bme——E（弹性模量）；b1——泊松比；kshear 开尔文抗剪模量, GK bvis1=vis1 ‘vis——kviscosity 开尔文粘性系数, ηK bvis2=vis2 end set bme=1.35714 set bke=3.25714 set b1=0.35 set vis1=1.35714 set vis2=19.2857 burcal model moh range group mixingpiles；定义本构 model burger range group soil any group exca1 any group exca2 any prop b b1 ksh bke msh bme kvis vis1 mvis vis2 range group soil any group exca1 any group exca2 any prop s 1e9 b 2e9 c 2e10 range group mixingpiles；赋材料属性
