FLAC渗流命令流共11页.docx

FLAC渗流命令流共11页.docx

《FLAC渗流命令流共11页.docx》由会员分享，可在线阅读，更多相关《FLAC渗流命令流共11页.docx（17页珍藏版）》请在冰豆网上搜索。

FLAC渗流命令流共11页

inipp0j=0

6q5D%E1m8G/f仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMfixppj=0

这样设置后是透水边界了

4|0Y（r.E%[;q:

[.l*HSimWe仿真论坛

底边界不透水,你可以不用管它了,默认为不透水边界

要练说，得练听。

听是说的前提，听得准确，才有条件正确模仿，才能不断地掌握高一级水平的语言。

我在教学中，注意听说结合，训练幼儿听的能力，课堂上，我特别重视教师的语言，我对幼儿说话，注意声音清楚，高低起伏，抑扬有致，富有吸引力，这样能引起幼儿的注意。

当我发现有的幼儿不专心听别人发言时，就随时表扬那些静听的幼儿，或是让他重复别人说过的内容，抓住教育时机，要求他们专心听，用心记。

平时我还通过各种趣味活动，培养幼儿边听边记，边听边想，边听边说的能力，如听词对词，听词句说意思，听句子辩正误，听故事讲述故事，听谜语猜谜底，听智力故事，动脑筋，出主意，听儿歌上句，接儿歌下句等，这样幼儿学得生动活泼，轻松愉快，既训练了听的能力，强化了记忆，又发展了思维，为说打下了基础。

我现在在三维中模拟基坑开挖，x（0~42）y（0~30）z（-40~0）,初始水位在－3米处，我的初始水位命令这样写对不对啊？

%F;d-]/F$C8vsimweinipp-3e4grad00-1e4rangex042y030z-40-3

&O/M6E9N4B,D6PSimWe仿真论坛range-3.1-2.9

$W8T5U/g4w-PSimWe仿真论坛;;

（B9]&U9b5F仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent在开挖中需要降水，降至－7米，需开挖的区域为x（0~8）y（0~30）

"~%N"z7F-o%`"b2a4P

fixpprangex842yo30z-3.1-2.9

L;_6m2v;J:

p7G#j9kSimWe仿真论坛fixpp0rangex08y030z-7.1-6.9plcontourpp

可以显示孔隙水压力场，一般就认为就是渗流场吧！

宋以后，京师所设小学馆和武学堂中的教师称谓皆称之为“教谕”。

至元明清之县学一律循之不变。

明朝入选翰林院的进士之师称“教习”。

到清末，学堂兴起，各科教师仍沿用“教习”一称。

其实“教谕”在明清时还有学官一意，即主管县一级的教育生员。

而相应府和州掌管教育生员者则谓“教授”和“学正”。

“教授”“学正”和“教谕”的副手一律称“训导”。

于民间，特别是汉代以后，对于在“校”或“学”中传授经学者也称为“经师”。

在一些特定的讲学场合，比如书院、皇室，也称教师为“院长、西席、讲席”等。

setflowon即选择是否要进行渗流分析，其前提是要setfluidon方可以！

off就是相反了。

但是，setfluidon以后不一定进行setflowon的计算，这点我的理解（有待证实）是，仅有setfluidon的话，就是把里面的pp当作是静水压力来处理的，而当setflowon打开以后，就可以进行真正意义上的渗流分析了！

再一个，setflowon的一个前提必须设定perm值（渗透系数），否则程序无法继续进行！

而setfluidon就没有那么严格，只要是configfluid以及设定模型以及fmod以后即可以开始计算了。

"f%w&Z）k2u9y7Q仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMsetmechon就表示把力学计算这部分打开，在默认的情况下，是打开的。

在流固耦合有一节中讲到，当只要使系统产生水压力的分布，而不需要由于机械变形（mechanicaldeformation）而产生的pp的时候，就必须要setmechoff，也就是说这个时候只有fluid的计算，而没有mechanical的计算了。

课本、报刊杂志中的成语、名言警句等俯首皆是,但学生写作文运用到文章中的甚少,即使运用也很难做到恰如其分。

为什么?

还是没有彻底“记死”的缘故。

要解决这个问题,方法很简单,每天花3-5分钟左右的时间记一条成语、一则名言警句即可。

可以写在后黑板的“积累专栏”上每日一换,可以在每天课前的3分钟让学生轮流讲解,也可让学生个人搜集,每天往笔记本上抄写,教师定期检查等等。

这样,一年就可记300多条成语、300多则名言警句,日积月累,终究会成为一笔不小的财富。

这些成语典故“贮藏”在学生脑中,自然会出口成章,写作时便会随心所欲地“提取”出来,使文章增色添辉。

帮助的第ftd138.pdf

bulkmodulus,K      390MPa  体积模量（土质比较硬）

shearmodulus,G      280MPa剪切模量

soildrydensity,ρd      1200kg/m3  土的密度

waterdensity,ρw      1000kg/m3  水密度

walldensity,ρwal      1500kg/m3  档土墙的密度

permeability,k      10−12m2/Pa-s  渗透系数

porosity,n      0.3  孔隙率

fluidbulkmodulus,Kf      2.0GPa  流体的体积模量

;-----以下是命令流

configfluid  ;设置流体算法

genzonebrickp11200p20120p30012size121212rat111

groupsoil;设置土组

groupexcavaterangex04y04z05;设置开挖范围组

groupwal1rangex45y05z07;设置档土墙1组

groupwal2rangex04y45z07;设置档土墙2组

groupwallrangegroupwal1anygroupwal2any;档土墙1和2合并为一个档土墙组

;---fluidflowmodel---  流体模型

modelfl_iso;各项同性流体

propperm1e-12poro0.3;设置渗透系数和？

？

比or？

？

率

inifdensity1e3;初始化水密度

inifmod2.0e9ftens-1e-3;初始化水的体积模量及？

？

modelfl_nullrangegroupexcavate  ;开挖部分没有水流

modelfl_nullrangegroupwall;档土墙不渗水

jinipp0grad001e4;初始化孔隙水压力

fixpprangez-0.10.1;约束应该是整个模型的底部吧？

fixpprangex-0.14.1y-0.14.1z4.95.1;开挖部分也给约束住了，而其他的地方不管，为什么要约束开挖部分？

modelelas;弹塑性

propbul3.9e6shea2.8e6  ;设置体积模量和剪切模量，这两个必须

modelnullrangegroupexcavate;力学模型的开挖

inidensity1.2e3;初始化密度，应该是土的密度

inidensity1.5e3rangegroupwall;初始化墙（档土墙）的密度

fixxrangex-.1.1;施加约束，据说差分法约束0的时候，必须取一个-0.1<0<0.1,范围

fixxrangex11.912.1;施加约束，或者叫做边界条件

fixyrangey-.1.1;施加约束

fixyrangey11.912.1

fixzrangez11.912.1

iniszz0grad00-1.5e4;z方向的应力大一些

inisxx0grad00-1.2e4;水平方向的，小一些

inisyy0grad00-1.2e4

applynstress0grad00-1.2e4rangex0.04.0y3.94.1z0.05.0

applynstress0grad00-1.2e4rangex3.94.1y0.04.0z0.05.0

applynstress-7.5e4rangex0.04.0y0.04.0z4.95.1

;---setting---  设置

setgravity0010;设置重力加速度

solveforce1;checkinitialequilibrium  求解平衡

;---excavation---  开挖

setfluidoff  ;关闭水

applynstress0grad00-1.e4rangex0.04.0y3.94.1z0.05.0

applynstress0grad00-1.e4rangex3.94.1y0.04.0z0.05.０

applynstress-5.e4rangex0.04.0y0.04.0z4.95.1

solve;force1;求解

saveexc1.sav;保存

;---drainage---排水

applyremovenstress;撤掉刚才的压力

defrelaxsetup  ;定义一个释放函数的参数

step0=step

end

relaxsetup;调用

defrelax;定义一个释放函数

rstep=step-step0

ifrstep

relax=1.0-（float（rstep）/float（ncyc））

else

relax=0.0

endif

end

setncyc=1000；先赋值，然后调用

applynstress0grad00-1e4hisrelaxrangex0.04.0y3.94.1z0.05.0

applynstress0 grad00-1e4hisrelaxrangex3.94.1y0.04.0z0.05.0

applynstress-5e4hisrelaxrangex0.04.0y0.04.0z4.95.1

cycncyc

solve

;---percolation---渗透

setfluidon;设置水算法

Fluentfixpp0rangex-0.14.1y-0.14.1z4.95.1;施加孔隙水压力

cyc9000;循环

saveexc3.sav;保存

你是想做渗流分析

'H*R&j8A%m+n%@5[|Simwe|仿真|设计|有限元|虚拟仪器还是静力分析

0o.C4Q!

?

7?

9E.jsimwe渗流分析的话

用configfluid做流固耦合分析

把手册中几个例子搞明白应该收获很大

（G$K2G%|8\#u）i仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM提醒生成初始应力场时

水平应力＝k0×竖向有效应力＋水应力

-y$|"Q3y9[|Simwe|仿真|设计|有限元|虚拟仪器如果仅仅是静力分析

只要设定稳定水位

手动干湿密度

#n/w&A9C1n$G#d.usimwe就行

与Interface相关的命令总结（望续）

最近用到这个了，因为初学尝尽苦头，以前整理的贴出来，望有别的还可以跟帖，simwe9j2g2K0`$^.p;j9@

仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM7v6Q8E*O1M%R4~

1）interface后面可以使用range关键字，用来限定其范围；rangeid用来设定界面元或节点的开始和最后的标识号码；

2l5M:

Y!

T$asimwe2）interface命令用来定义两个或多个子网格之间的界面，这些界面可以滑动、分离。

simwe2|5[-];];s8z7sb1V

3）Historyinterfacenkeywword:

对界面n取样其界面变量：

关键子包括：

（f.V-ly5y:

\!

A*@!

?

仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent    nstress:

法向应力SimWe仿真论坛4x6P7n5~&F

    ndisplacement:

法向位移

:

W'V4~$T8u9oi4}仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM    sstress:

剪切应力仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent#v:

F.z$b1M（t5B1H

    sdisplacement:

剪切位移仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent6x"S/S0_*M;j3X7M

4）interface命令的标准格式如下：

interfaceikeyword：

界面是子网格之间相互作用的平面。

界面用来模拟体与体之间节理、断层和摩擦界面的影响。

界面可以附着于一个子网格上，也可以放置在空间内的任何位置。

当子网格接触界面时，产生剪切力和法向力。

允许沿着界面发生滑动或分离。

界面由连接在节点顶点的三角形元组成。

每一个界面元有一个主动面和一个被动面，由界面的法向矢量的方向来确定。

界面节点可以直接附着到一个网格面上或者放置在空间内的任何位置。

下面的关键字用来产生一个界面。

仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent"@7m6P%a7A!

_9B

elementp1%A4U6b+l!

e2D1L2Y

        p2|Simwe|仿真|设计|有限元|虚拟仪器v"M;f!

E#n+}K.i

        p3

:

o1d6?

U$N9C2j2n*b这个命令产生一个三角形的界面元。

p1,p2,p3是三角形的三个顶点。

这个单元可以用空间内的三个顶点的坐标来产生，或者连接到已经存在的界面节点i上。

即使产生的界面元的位置对应于面的位置，它也没有附着到网格面上。

这个单元固定在空间位置。

单元的主动面由其所组成的边的走向来定义，从p1到p2到p3，当按顺时针方向排列时，主动面朝上。

/}$I3q;\&~）csimweface:

在设定范围内的所有网格面上产生界面元。

如果在选择的面上已经存在有界面元，会发生错误。

2S-b7g4c1X#@+d|Simwe|仿真|设计|有限元|虚拟仪器maxedgevalue：

value是界面上所有界面元的最大边长。

界面被子划分直到所有单元的边长比设定的值小。

默认情况下，对每个四边形的单元面产生两个三角形单元。

|Simwe|仿真|设计|有限元|虚拟仪器9_3g!

[&N#S4v

nodexyz：

在位置（x,y,z）上产生一个界面节点。

可以用选择的关键字id来设定ID号码，否则自动决定ID。

如果在选择的位置上已经存在一个节点，会报告有错误发生。

产生的节点被固定在空间中。

仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM）M7K%]'p"a.Z{7w

nstressx:

对选择的范围内所有界面号码i节点设置x方向的法向应力增量；法向应力增量被增加到由界面贯入计算的法向应力。

6I*[（K:

w$W*d&s,?

;Psimwesstressxyz：

对选择范围内的所有界面号码i节点设置增量的剪切应力xyz。

6p5_"u*t3I1V+D仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM可以用下面的关键字删除界面：

b.n!

m1[（E;i:

Y|Simwe|仿真|设计|有限元|虚拟仪器delete:

在选择范围内的与该界面相关的所有界面元被删除。

（Z）R6M,_3p;d0O&f3A|Simwe|仿真|设计|有限元|虚拟仪器用Coulomb滑动和/或拉伸分离来表征界面。

用property关键字设定界面的性质，在所有的情况下，必须设定法向刚度kn和剪切刚度ks，对这些性质必须使用一致的单位设置，如果不设置性质为0。

+I9c#v9L"B2j+J4t9O7O:

P仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMpropertykeywordvalue...  属性关键字包括：

&m4F&j-W;G9G:

V5a仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentcohesion:

粘结力[应力]simwe!

V6h%];g4x*f2d5c

dialation:

膨胀角[角度]!

_+L5|5P2O3R7M0s

friction:

摩擦角[角度]

0l%X3Y6I$q:

E-v1`%P.u仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentkn:

法向刚度[应力/位移]3{%~6`$Q1e7a

ks:

切向刚度[应力/位移],N:

~1|+y$Y6n（c,w（R"|

tension:

抗拉强度[应力]SimWe仿真论坛3r7@5h-z9c（j"i

如果沿着界面发生滑动，那么粘结力的值不改变；如果超过抗拉强度，那么抗拉强度就设置为0。

|Simwe|仿真|设计|有限元|虚拟仪器+j1c5U3L*vz3@

*Y#K0_4P#h3^7u3H9[仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent5）plotbcontour/contourgpnumon（默认为off）绘制界面节点IDs的网格点数目。

1T（y9T6I1x'u"]'L/h9L/o|Simwe|仿真|设计|有限元|虚拟仪器6）plotinterface下面选择的关键字绘制界面变量的等值线。

simwe,h/?

*R8i!

L"D+V*Q+O;_

location:

仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM3f6C8C+@,?

+o$Y8C

penetration:

界面节点的相互贯入量仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM4F%KH&y!

q:

l8[

ndisplacement:

相对法向位移（与贯入相同）0}n&d;n,B1V"^

nstress:

在界面元上的法向应力

-\6j7J4s"S+Q1R'{-Tsimwesdisplacement:

相对剪切位移|Simwe|仿真|设计|有限元|虚拟仪器K5g$d2l）q5^._3|

sstress:

在界面元上的剪切应力

-t#\$f6X-T1l:

J）M5x'p0h可以使用下面的开关值：

color,activiteon/off,alias,faston/off,gpnumon/off,idon/off,intervalv,maximumv,minimumv,nullon/off,outlineon/off,reverseon/off,shadeon/off

/I"S;~5[3f-t.lsimwe7）printinterface:

界面数据总结？

？

？

望续此处

4n-f&W*x-q（]8c）Y|Simwe|仿真|设计|有限元|虚拟仪器|Simwe|仿真|设计|有限元|虚拟仪器+j&n7_%y（Z:

V

[本帖最后由benjackxu于2019-9-1320:

06编辑]

与Interface相关的命令总结（望续）

最近用到这个了，因为初学尝尽苦头，以前整理的贴出来，望有别的还可以跟帖，simwe9j2g2K0`$^.p;j9@

仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM7v6Q8E*O1M%R4~

1）interface后面可以使用range关键字，用来限定其范围；rangeid用来设定界面元或节点的开始和最后的标识号码；

2l5M:

Y!

T$asimwe2）interface命令用来定义两个或多个子网格之间的界面，这些界面可以滑动、分离。

simwe2|5[-];];s8z7sb1V

3）Historyinterfacenkeywword:

对界面n取样其界面变量：

关键子包括：

（f.V-ly5y:

\!

A*@!

?

仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSO