完整word版,FLAC3D中应变软化硬化模型参数取值(table)

一FLAC3D中应变软化\硬化模型参数取值（table）
例：model ss
prop den 2500 bulk 2e8 shear 1e8 co 2e6 fric 45 ten 1e6 dil 10
prop ftab 1 ctab 2 dtab 3
table 1 0，45 0. 05，42 0. 1，40
table 2 0，2e6 0. 05，1e6 0. 1，5e5
table 3 0，10 0. 05，3 0. 1，0
Table 的功能是，用来定义某一变量的值随另一变量的值而发生变化的变化规律。

例如，对于应变软化模型（strain-softening model ），可使用PROP ftable, PROP ctable 和PROP dtable等命令来设置摩擦角（friction）、黏聚力（cohesion）和剪胀角（dilation），随着累积塑性剪切应变（accumulated plastic shear strain ）而变化的规律；
或者使用PROP ttable命令，来定义抗拉强度（tensile strength ）随着累积塑性拉应变（accumulated plastic tensile strain ）而变化的规律。

Table命令的调用格式为：
TABLE n <keyword> x1 y1 <x2 y2> <x3 y3> . . .
定义1号table（ftab 1；即定义摩擦角的变化规律），
累积塑性切应变取值为0 时，摩擦角取值为45；
累积塑性切应变取值为0.05 时，摩擦角取值为42；
累积塑性切应变取值为0.1 时，摩擦角取值为40；。

Flac3D中⽂⼿册Flac3D 中⽂⼿册FLAC3D的计算模式中是否需要做孔压分析取决于是否采⽤config fluid命令。

1 ⽆渗流模式（不使⽤config fluid）即使不使⽤命令config fluid，仍然可以在节点上施加孔压。

这种模式下，孔压将保持为常量。

如果采⽤塑性本构模型的话，材料的破坏将由有效应⼒状态来控制。

节点上的孔压分布可由initial pp命令或water table命令来设定。

如果采⽤water table命令，由程序⾃动计算⽔位线以下的静⽔孔压分布。

此时，必须施加流体密度（water density）和重⼒（set gravity）。

流体密度值和⽔位位置可以⽤命令print water显⽰。

如果⽔位线是由face关键字来定义的，则可⽤命令plot water命令显⽰⽔位。

这两种情况，单元的孔压都由节点孔压值平均求出，并在本构模型计算中⽤作有效应⼒。

这种计算模式下，体积⼒中不反映流体的出现：⽤户必须根据⽔位线以上或以下相应地指定⼲密度和湿密度。

使⽤命令print gp pp和priint zone pp可分别得到节点或单元孔压。

plot contour pp命令可绘出节点孔压云图。

2 渗流模式（使⽤config fluid）如果使⽤命令config fluid，则可进⾏瞬时渗流分析，孔压改变和潜⽔⾯的改变都可能出现。

在config fluid模式下，有效应⼒计算（静态孔压分布）和⾮排⽔计算均被执⾏。

除此之外，还可进⾏全耦合分析，这种情况下，孔压改变将使固体产⽣变形，同时体积应变反过来影响孔压的变化。

如果采⽤渗流模式，单元孔压仍由节点孔压平均求出。

但这种模式，⽤户只能指定⼲密度（不论是⽔位以上还是以下），因为FLAC3D 将流体的影响考虑到了体积⼒的计算中。

采⽤渗流模式时，渗流模型必须施加到单元上，使⽤命令modelfl_isotropic模拟各向同性渗流，model fl_anisotropic模拟各向异性渗流，model fl_null模拟⾮渗透物质。

1-1栈型多it代码可参考m a n u a l中冬个章节的c o m m a n d命令及说朗，注意单住。

用p r o p赋值。

正交冬的异性暉性挨型＜鲁克•普柱格棋晏余比松棋型＞节現棋型丄应隻戌化/缺化模型修正创林棋愛他劝力爭栈型形挨型均术淙体挨型冬衙异技流体棋型均朋範导棋型1.2履熨追用说明遍布节理核世迨用于Mohr-Coulomb材蚪来期确显示力在冬个方向上的差异性。

玖线性欲化应变遍布节理栈空综合了软化应变Mohr-Coulomb棋黑和遍布节理栈空，这种栈世包含而向矩阵和遍布节理的一个玖线性斯裂点集。

改迸.的Cam-clay栈空反映了形支度和抗破坏能力对体孩雯化的彩响。

Mohr-Coulomb核型呆迨用于一般工程研克，同肘，Mohr-Coulomb的聚力和摩擦角茨数相对于地施工程材料的其它為性，灵汆易莪得。

飲化应变和遍布节理塑性栈禮实际上是Mohr-Coulomb钱世的支形，这些楼型如果在附加材蚪参数的值较离时将得出与Mohr-Coulomb核空同择的结果。

Druck-Prager橈型是一个相对于Mohr-Coulomb棋型的破坏标准的简化体.但是它一般不迨于用来福述地质工程材料的破坏情况。

它主要是用来杷FLAC3D与其它一些有Druck-Prager核型但却没冇Mohr-Coulomb侯空的数学软件作比较。

札哮撫力为零的对候诗注意,此肘Mohr-Coulomb钱型退化为Tresca核童，而Druck-Prager 棋空退化为Von Mises栈空。

Druck-Prager楼禮和Mohr-Coulomb棋型旻计算起来效率呆爲的型性侯型.而其它的复性核禮A计算时却需要灵多的存和额外的对间。

例如，題性应变不能在Mohr-Coulomb 楼型中直挂计算出来（参见附录G 儿如果需要计算覆性应变，则必常要用应雯轶化核型。

这科核世主要是用于破坏后的悄况对工程彩响重大的工程活动中，如专曲柱.开采堀掠戎凹珑研兗。

1.1模型参数代码可参考manual中各个章节的command命令及说明，注意单位。

用prop 赋值。

1.1.1各向同性弹性模型1.1.2横向同行弹性模型1.1.3正交各向异性弹性模型1.1.4德鲁克-普拉格模型1.1.5摩尔-库伦模型1.1.6多节理模型1.1.7应变硬化/软化模型1.1.8双线性应变硬化/软化多节理模型1.1.9D-Y模型1.1.10修正剑桥模型1.1.11纯动力学模型1.1.12经典粘弹性模型1.1.13粘弹性模型1.1.14二分幂律模型7 rs_2 参考应力，σ2ref8 shear 弹性剪切模量，G1.1.15蠕变模型蠕变模型材料参数（WIPP Model –MODEL mechanical wipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_wipp 常数，B4 bulk 弹性体积模量，K5 d_wipp 常数，D6 e_dot_star临界稳定状态蠕变率，7 gas_c 气体常数，R8 n_wipp 指数，n9 shear 弹性剪切模量，G10 temp 温度，T下列参数可以显示、绘图和通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率1.1.16Burger、蠕变组合材料模型Burger、蠕变组合材料模型的材料参数（Burgers-Creep Viscoplastic Model –MODEL mechanical cvisc）1 bulk 弹性体积模量，K2 cohesion 聚力，c3 density 密度，ρ4 dilation 剪胀角，Ψ1.1.17幂律模型1.1.18粘塑形模型粘塑形模型的材料参数（WIPP-Creep Viscoplastic Model –MODEL mechanical pwipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_wipp 常数，B4 bulk 弹性体积模量，K5 d_wipp 常数，D6 e_dot_star临界稳定状态蠕变率，7 gas_c 气体常数，R8 kshear 材料参数，KΦ9 n_wipp 指数，n10 kdil 材料参数，q k11 kvol 材料参数，qΦ12 shear 弹性切变模量，G13 temp 温度，T14 tension 抗拉强度，σt以下计算参数可以显示、绘图和通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率3 es_plastic 累积塑性切应变4 et_plastic 累积塑性拉应变1.1.19碎盐变形模型碎盐变形模型的材料参数（Crushed-Salt Model –MODEL mechanical cwipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_f 最终体积模量，K f4 b_wipp 常数，B5 b0 蠕变压实系数，B06 b1 蠕变压实系数，B17 b2 蠕变压实系数，B28 bulk 弹性体积模量，K9 d_f 最终密度，ρf10 d_wipp 常数，D11 e_dot_star临界稳定状态蠕变率，12 gas_c 气体常数，R13 n_wipp 指数，n14 rho 密度，ρ15 s_f 最终切变模量，G f16 shear 弹性切变模量，G17 temp 温度，T以下计算参数可以显示、绘图和通过fish访问1 frac_d 当前碎片密度，ρd2 s_g1 蠕变压实参数，G3 s_k1 蠕变压实参数，K1.1.20均质流体模型均质流体模型的材料参数1 permeability 等方向渗透性，k2 porosity 孔隙率，n（默认时，n=0.5）1.1.21各向异性流体模型1.1.22均质热导模型1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。

计算说明1、计算方法1）内力计算采用弹性支点法；2）土的水平抗力系数按 M法确定；3）主动土压力与被动土压力采用矩形分布模式；4)采用力法分析环形内支撑内力；5）采用"理正深基坑支护结构软件 FSPW 5.2"计算，计算采用单元计算与协同计算相结合，并采用FLAC-3D进行验证；6）土层参数选取表1 各土层参数取值表2、单元计算1）基坑分为4个区，安全等级为一级，基坑重要性系数为 1.1；2）荷载：施工荷载：10kPa；地面超载：4区活动荷载为25kPa，1区、2区和3区超载按10kPa考虑；水压力；基坑外侧为常水位，内侧坑底以下 0.5m。

3）基坑开挖深度：根据现场地形确定，按开挖12.50m确定；4）支撑水平刚度系数：式中α取0.8，E取28000MPa，L取7.0m，sa取1.20m，s取7.0m，经计算，kT大于800 MN/m，本计算中，取800MN/m。

5）计算过程详见附件1，其中1区选用钻孔ZK1，2区选用钻孔ZK4，3区选用钻孔ZK16，4区选用钻孔ZK5。

各区计算结果汇总如下：表2 计算结果汇总表3、协同计算1）计算方法简介协同计算采用考虑支护结构、内支撑结构及土空间整体协同作用有限元的计算方法。

有限元方程如下：（[K n]+[Kz]+[Kt]）{W)}={F}式中：[K n]－内支撑结构的刚度矩阵；[K z]－支护结构的刚度矩阵；[Kt]－开挖面以下桩侧土抗力的刚度矩阵；{W}－位移矩阵；{F}－荷载矩阵。

计算时采用如下简化计算方法：（1）将基坑周边分成几个计算区域，同一计算区域的支护信息相同，地质条件相同。

（2）将每一个桩或每单位长度的墙看成是一个超级的子结构，这一子结构包括桩墙，土，主动和被动土压力。

（3）将第三道锚索等效为弹性支承点，作为支承系统的一部份进行计算。

（4）单独求解（2)中的子结构，可采用单桩内力计算的一套方法，将刚度和荷载凝聚到与支锚的公共节点上，这是一个一维梁计算问题。

FLAC3D常见命令与使用技巧1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。

2.最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

3.怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布4.怎样看模型的边界情况？答：plo gpfix red5.怎样看模型的体力分布？答：plo fap red6.怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx,sxy, syz, sxz)7.怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8.怎样看模型有多少单元、节点？答：pri info9.怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10.怎样调用一个文件？答：File/call或者call命令10.如何施加面力？答：app nstress11.如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键12.如何进行边界约束？答：fix x ran（约束的是速度，在初始情况下约束等效于位移约束）13.如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14.如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片)15.如何保存计算结果？答：save +文件名.16.如何调用已保存的结果？答：rest +文件名；或者File / Restore17.如何暂停计算？答：Esc18.如何在程序中进行暂停，并可恢复计算？答：在命令中加入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 end2radius r2 not，其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

1.1模型参数代码可参考manual中各个章节的command命令及说明,注意单位。

用prop 赋值。

16 ttable 塑性拉应变-抗拉强度的表号下列参数可以显示、绘图与通过fish访问1 es_plastic 塑性切应变2 et_plastic 塑性拉应变3 ff_count 检测切应变反向的数4 ff_cvd 体应变,εvd经典粘弹性模型的材料参数(Classical Viscoelastic (Maxwell Substance) –MODEL mechanical viscous)1 bulk 弹性体积模量,K2 shear 弹性剪切模量,G3 viscosity 动力粘度,η粘弹性模型的材料参数(Burgers Model –MODEL mechanical burgers)1 bulk 弹性体积模量,K2 kshear Kelvin弹性剪切模量,G K3 kviscosity Kelvin动力粘度,ηK4 mkshear Maxwell切边模量,G M5 mviscosity Maxwell动力粘度,ηM二分幂律模型的材料参数(Power Law –MODEL mechanical power)1 a_1 常数,A12 a_2 常数,A23 bulk 弹性体积模量,K4 n_1 指数,n15 n_2 指数,n26 rs_1 参考应力,σ1ref7 rs_2 参考应力,σ2ref8 shear 弹性剪切模量,G蠕变模型材料参数(WIPP Model –MODEL mechanical wipp)1 act_energy 活化能,Q2 a_wipp 常数,A3 b_wipp 常数,B4 bulk 弹性体积模量,K5 d_wipp 常数,D6 e_dot_star临界稳定状态蠕变率,7 gas_c 气体常数,R8 n_wipp 指数,n9 shear 弹性剪切模量,G10 temp 温度,T下列参数可以显示、绘图与通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率Burger、蠕变组合材料模型的材料参数(Burgers-Creep Viscoplastic Model –MODEL mechanical cvisc)1 bulk 弹性体积模量,K2 cohesion 内聚力,c3 density 密度,ρ4 dilation 剪胀角,Ψ5 friction 内摩擦角,Φ6 kshear Kelvin弹性剪切模量,G K7 kviscosity Kelvin粘度,ηK8 shear 弹性剪切模量,G9 tension 抗拉强度,σt10 mviscosity Maxwell动力粘度,ηM下列计算参数可以显示、绘图与通过fish访问1 es_plastic 累积塑性切应变2 et_plastic 累积塑性拉应变幂律模型的材料参数(Power-Law Viscoplastic Model –MODEL mechanical cpower)1 a_1 常数,A12 a_2 常数,A23 bulk 弹性体积模量,K4 cohesion 内聚力,c5 dilation 剪胀角,Ψ6 friction 内摩擦角,Φ7 n_1 指数,n18 n_2 指数,n29 rs_1 参考应力,σ1ref10 rs_2 参考应力,σ2ref11 shear 弹性剪切模量,G12 tension 抗拉强度,σt粘塑形模型的材料参数(WIPP-Creep Viscoplastic Model –MODEL mechanical pwipp)1 act_energy 活化能,Q2 a_wipp 常数,A3 b_wipp 常数,B4 bulk 弹性体积模量,K5 d_wipp 常数,D6 e_dot_star临界稳定状态蠕变率,7 gas_c 气体常数,R8 kshear 材料参数,KΦ9 n_wipp 指数,n10 kdil 材料参数,q k11 kvol 材料参数,qΦ12 shear 弹性切变模量,G13 temp 温度,T14 tension 抗拉强度,σt以下计算参数可以显示、绘图与通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率3 es_plastic 累积塑性切应变4 et_plastic 累积塑性拉应变碎盐变形模型的材料参数(Crushed-Salt Model –MODEL mechanical cwipp)1 act_energy 活化能,Q2 a_wipp 常数,A3 b_f 最终体积模量,K f4 b_wipp 常数,B5 b0 蠕变压实系数,B06 b1 蠕变压实系数,B17 b2 蠕变压实系数,B28 bulk 弹性体积模量,K9 d_f 最终密度,ρf10 d_wipp 常数,D11 e_dot_star临界稳定状态蠕变率,12 gas_c 气体常数,R13 n_wipp 指数,n14 rho 密度,ρ15 s_f 最终切变模量,G f16 shear 弹性切变模量,G17 temp 温度,T以下计算参数可以显示、绘图与通过fish访问1 frac_d 当前碎片密度,ρd2 s_g1 蠕变压实参数,G3 s_k1 蠕变压实参数,K均质流体模型的材料参数1 permeability 等方向渗透性,k2 porosity 孔隙率,n(默认时,n=0、5)各向异性流体模型的材料参数1 fdd k1-k2的平面倾向2 fdip k1-k2的平面倾角1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。

1.1模型参数代码可参考manual中各个章节的command命令及说明,注意单位.用prop 赋值。

1.1.12经典粘弹性模型的材料参数(Classical Viscoelastic （Maxwell Substance) –MODEL mechanical viscous）1bulk弹性体积模量，K2shear弹性剪切模量，G3viscosity动力粘度，η1.1.13粘弹性模型粘弹性模型的材料参数（Burgers Model –MODEL mechanical burgers）1bulk弹性体积模量，K2kshear Kelvin弹性剪切模量，G K3kviscosity Kelvin动力粘度，ηK4mkshear Maxwell切边模量，G M5mviscosity Maxwell动力粘度,ηM1.1.14二分幂律模型二分幂律模型的材料参数（Power Law –MODEL mechanical power）1a_1常数,A12a_2常数,A23bulk弹性体积模量，K4n_1指数，n15n_2指数,n26rs_1参考应力,σ1ref7rs_2参考应力，σ2ref8shear弹性剪切模量，G1.1.15蠕变模型蠕变模型材料参数（WIPP Model –MODEL mechanical wipp）1act_energy活化能，Q2a_wipp常数，A3b_wipp常数,B4bulk弹性体积模量,K5d_wipp常数，D6e_dot_star临界稳定状态蠕变率,7gas_c气体常数，R8n_wipp指数,n9shear弹性剪切模量，G10temp温度，T下列参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变2e_rate累积主蠕变应变率Burger、蠕变组合材料模型的材料参数(Burgers—Creep Viscoplastic Model –MODEL mechanical cvisc)1bulk弹性体积模量,K2cohesion内聚力，c3density密度，ρ4dilation剪胀角，Ψ5friction内摩擦角,Φ6kshear Kelvin弹性剪切模量，G K7kviscosity Kelvin粘度,ηK8shear弹性剪切模量,G9tension抗拉强度，σt10mviscosity Maxwell动力粘度，ηM下列计算参数可以显示、绘图和通过fish访问1es_plastic累积塑性切应变2et_plastic累积塑性拉应变1.1.17幂律模型幂律模型的材料参数（Power-Law Viscoplastic Model –MODEL mechanical cpower）1a_1常数，A12a_2常数,A23bulk弹性体积模量，K4cohesion内聚力，c5dilation剪胀角，Ψ6friction内摩擦角，Φ7n_1指数，n18n_2指数,n29rs_1参考应力,σ1ref10rs_2参考应力，σ2ref11shear弹性剪切模量，G12tension抗拉强度，σt1.1.18粘塑形模型粘塑形模型的材料参数(WIPP—Creep Viscoplastic Model –MODEL mechanical pwipp）1act_energy活化能，Q2a_wipp常数,A3b_wipp常数，B4bulk弹性体积模量,K5d_wipp常数，D6e_dot_star临界稳定状态蠕变率，7gas_c气体常数，R8kshear材料参数，KΦ9n_wipp指数，n10kdil材料参数,q k11kvol材料参数，qΦ12shear弹性切变模量,G13temp温度，T14tension抗拉强度，σt以下计算参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变2e_rate累积主蠕变应变率3es_plastic累积塑性切应变4et_plastic累积塑性拉应变碎盐变形模型的材料参数(Crushed—Salt Model –MODEL mechanical cwipp）1act_energy活化能,Q2a_wipp常数,A3b_f最终体积模量，K f4b_wipp常数，B5b0蠕变压实系数,B06b1蠕变压实系数，B17b2蠕变压实系数，B28bulk弹性体积模量，K9d_f最终密度，ρf10d_wipp常数,D11e_dot_star临界稳定状态蠕变率，12gas_c气体常数,R13n_wipp指数，n14rho密度，ρ15s_f最终切变模量，G f16shear弹性切变模量，G17temp温度，T以下计算参数可以显示、绘图和通过fish访问1frac_d当前碎片密度，ρd2s_g1蠕变压实参数,G3s_k1蠕变压实参数，K均质流体模型的材料参数1permeability等方向渗透性，k2porosity孔隙率,n(默认时，n=0。

1.FLAC3D知识基本介绍SimWe岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。

由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的差分法。

差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。

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 对模拟塑性破坏和塑性流动采用的是“混合离散法“。

第1部分命令流按照顺序进行2-1定义一个FISH函数newdef abcabc = 25 * 3 + 5Endprint abc2-2使用一个变量newdef abchh = 25abc = hh * 3 + 5EndPrint hhPrint abc2-3对变量和函数的理解newdef abchh = 25abc = hh * 3 + 5Endset abc=0 hh=0print hhprint abcprint hhnewdef abcabc = hh * 3 + 5endset hh=25print abcset abc=0 hh=0print hhprint abcprint hh2-4获取变量的历史记录newgen zone brick size 1 2 1model mohrprop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10fix x y z range y -0.1 0.1apply yvel -1e-5 range y 1.9 2.1plot set rotation 0 0 45plot block groupdef get_adad1 = gp_near(0,2,0)ad2 = gp_near(1,2,0)ad3 = gp_near(0,2,1)ad4 = gp_near(1,2,1)endget_addef loadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) endhist loadhist gp ydis 0,2,0step 1000plot his 1 vs -22-5用FISH函数计算体积模量和剪砌模量newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25deriveprint b_modprint s_mod2-6 在FLAC输入中使用符号变量Newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25derivegen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear2-7 控制循环Newdef xxxsum = 0prod = 1loop n (1,10)sum = sum + nprod = prod * nend_loopendxxxprint sum, prodnewgen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10) model elasplot set rotation 0 0 45plot block groupdef installpnt = zone_headloop while pnt #nullz_depth = -z_zcen(pnt)y_mod = y_zero + cc * sqrt(z_depth)z_prop(pnt, ’shear’) = y_mod / (2.0*(1.0+p_ratio))z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio))pnt = zone_next(pnt)end_loopendset p_ratio=0.25 y_zero=1e7 cc=1e8install2-8 拆分命令行new ;example of a sum of many thingsdef long_sumtemp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10 long_sum = temp + v11 + v12 + v13 + v14 + v15end2-9 变量类型newdef haveoneaa = 2bb = 3.4cc = ’Have a nice day’dd = aa * bbee = cc + ’, old chap’endhaveoneprint fish2-10 IF条件语句newdef abcif xx > 0 thenabc = 1000elseabc = -1000end_ifendset xx = 10print abcset xx = 0print abc2-11 索单元自动生成newgen zone brick size 10 3 5plot set rotation 0 0 45plot block groupdef place_cablesloop n (1,5)z_d = float(n) - 0.5commandsel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7 end_commandend_loopendplace_cablesplot grid sel geom rednewgen zone brick size 10 3 5plot set rotation 15 0 60plot block groupmod mohrprop bulk 1e8 shear .3e8 fric 35prop coh 1e3 tens 1e3ini dens 1000set grav 0,0,-10fix x y z range z -.1 .1fix y range y -.1 .1fix y range y 2.9 3.1fix x range x -.1 .1fix x range x 9.9 10.1set largehist unbalsolvesave cab_str.savini xdis 0 ydis 0 zdis 0hist gp xdisp 0,1,5def place_cablesloop n (1,5)z_d = 5.5 - float(n)z_t = z_d + 0.5z_b = z_d - 0.5commandfree x range x -.1,.1 z z_b z_tsolvesel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 &gr_k 2e10 gr_coh 1e10 gr_per 1.0end_commandend_loopendplace_cablessave cab_end.savplot sketch sel cable force red2-12圆形隧道开挖模拟计算;建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 &size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderockgen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock gen zon reflect dip 90 dd 90 orig 0 0 0gen zon reflect dip 0 dd 0 ori 0 0 0gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1 gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2 gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3 gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1;绘制模型图plot block groupplot add axes red;plot set rotation 0 0 45 用于显示三维模型;设置重力set gravity 0 0 -10;给定边界条件fix z range z -12.01,-11.99fix x range x -21.01,-20.99fix x range x 20.99,21.01fix y range y -0.01 0.01fix y range y 0.99,1.01;求解自重应力场model mohrini density 1800 ;围岩的密度prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度set mech ratio=1e-4solvesave Gravsol.savplot cont zdisp outl onplot cont szz;毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0model null range group insiderock any group concretliner anyplot block groupplot add axes redset mech ratio=5e-4solvesave Kaiwsol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;模筑衬砌计算model elas range group concretliner anyplot block groupplot add axes redini density 2500 range group concretliner any ;衬砌混凝土的密度prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模set mech ratio=1e-4solvesave zhihusol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;完成计算分析。