cyc 500
fix pp ran x 16 40 y 0 20 z -2.1 -1.9 fix pp 0 ran x -0.1 16.1 y -0.1 20.1 z -7.1 -6.9 group excavate ran x 0 16 y 0 20 z -7 -3
model fl_null ran group excavate model fl_null ran sel group wall model fl_null ran group soil1 model null ran group excavate
sel node cid 10000 0 4 -2 sel node cid 10001 0 8 -2 sel node cid 10002 0 12 -2 sel node cid 10003 0 16 -2
sel beamsel 10000 509 sel beamsel 10001 513 sel beamsel 10002 517 sel beamsel 10003 521
sel node free xr yr zr z y ran cid 10000 509 sel node free xr yr zr z y ran cid 10001 513 sel node free xr yr zr z y ran cid 10002 517 sel node free xr yr zr z y ran cid 10003 521
sel beam prop dens 7800 emod 210e9 nu 0.25 xcarea 0.0364 xcj 0.0064 xciy 0.0032 xciz 0.0032 ydirection 0 1 0 set fluid on cyc 500
fix pp ran x 16 40 y 0 20 z -2.1 -1.9 fix pp 0 ran x -0.1 16.1 y -0.1 20.1 z -11.1 -10.9 group excavate ran x 0 16 y 0 20 z -11 -7
model fl_null ran group excavate model fl_null ran sel group wall model fl_null ran group soil1 model null ran group excavate
sel node cid 10004 0 4 -6 sel node cid 10005 0 8 -6 sel node cid 10006 0 12 -6 sel node cid 10007 0 16 -6
sel beamsel 10004 425 sel beamsel 10005 429 sel beamsel 10006 433 sel beamsel 10007 437
sel node free xr yr zr z y ran cid 10004 425 sel node free xr yr zr z y ran cid 10005 429 sel node free xr yr zr z y ran cid 10006 433 sel node free xr yr zr z y ran cid 10007 437
sel beam prop dens 7800 emod 210e9 nu 0.25 xcarea 0.0364 xcj 0.0064 xciy 0.0032 xciz 0.0032 ydirection 0 1 0 set fluid on cyc 500
fix pp ran x 16 40 y 0 20 z -2.1 -1.9 fix pp 0 ran x -0.1 16.1 y -0.1 20.1 z -17.1 -16.9 group excavate ran x 0 16 y 0 20 z -17 -11
model fl_null ran group excavate model fl_null ran sel group wall model fl_null ran group soil1 model null ran group excavate
sel node cid 10008 0 4 -10 sel node cid 10009 0 8 -10 sel node cid 10010 0 12 -10 sel node cid 10011 0 16 -10
sel beamsel 10008 341 sel beamsel 10009 345 sel beamsel 10010 349 sel beamsel 10011 353
sel node free xr yr zr z y ran cid 10008 341 sel node free xr yr zr z y ran cid 10009 345 sel node free xr yr zr z y ran cid 10010 349 sel node free xr yr zr z y ran cid 10011 353
sel beam prop dens 7800 emod 210e9 nu 0.25 xcarea 0.0364 xcj 0.0064 xciy 0.0032 xciz 0.0032 ydirection 0 1 0 set fluid on cyc 500
sel node cid 10012 0 4 -14 sel node cid 10013 0 8 -14 sel node cid 10014 0 12 -14 sel node cid 10015 0 16 -14
sel beamsel 10012 257 sel beamsel 10013 261 sel beamsel 10014 265 sel beamsel 10015 269
sel node free xr yr zr z y ran cid 10012 257 sel node free xr yr zr z y ran cid 10013 261 sel node free xr yr zr z y ran cid 10014 265 sel node free xr yr zr z y ran cid 10015 269
sel beam prop dens 7800 emod 210e9 nu 0.25 xcarea 0.0364 xcj 0.0064 xciy 0.0032 xciz 0.0032 ydirection 0 1 0 hist gp xdisp 16 10 40 ;2 hist gp xdisp 16 10 38 ;3 hist gp xdisp 16 10 34 ;4 hist gp xdisp 16 10 30 ;5 hist gp xdisp 16 10 26 ;6 hist gp xdisp 16 10 23 ;7
hist gp xdisp 17 10 40 ;8 hist gp xdisp 17 10 38 ;9 hist gp xdisp 17 10 34 ;10 hist gp xdisp 17 10 30 ;11 hist gp xdisp 17 10 26 ;12 hist gp xdisp 17 10 23 ;13
20
hist gp xdisp 18 10 40 ;14 hist gp xdisp 20 10 40 ;15 hist gp xdisp 23 10 40 ;16 hist gp xdisp 26 10 40 ;17 hist gp xdisp 29 10 40 ;18 hist gp xdisp 32 10 40 ;19 hist gp xdisp 35 10 40 ;20 hist gp xdisp 38 10 40 ;21 hist gp xdisp 40 10 40 ;22
hist gp zdisp 16 10 40 ;23 hist gp zdisp 17 10 40 ;24 hist gp zdisp 18 10 40 ;25 hist gp zdisp 20 10 40 ;26 hist gp zdisp 23 10 40 ;27 hist gp zdisp 26 10 40 ;28 hist gp zdisp 29 10 40 ;29 hist gp zdisp 32 10 40 ;30 hist gp zdisp 35 10 40 ;31 hist gp zdisp 38 10 40 ;32 hist gp zdisp 40 10 40 ;33 hist gp zdisp 0 10 23 ;34 hist gp zdisp 3 10 23 ;35 hist gp zdisp 4 10 23 ;36 hist gp zdisp 6 10 23 ;37 hist gp zdisp 8 10 23 ;38 hist gp zdisp 9 10 23 ;39 hist gp zdisp 12 10 23 ;40 hist gp zdisp 15 10 23 ;41 cyc 500 ret
\单元与shell单元共享节点了吗？
你用 interface 和 shell 来模拟地下连续墙是出于什么考虑？ 似乎有点勉强 sel node cid 10012 0 4 -14\单元与shell单元共享节点了吗？ 感谢你的回复！！ 共享了，beam单元的第二节点与壳单元节点是同一节点。 to dynamax: 感谢回复！！ 我不是用interface来模拟地下连续墙的，是用来模拟墙与土之间的接触相互作用的。我是用shell单元来模拟地下连续墙的，但是弯矩很小，并且支撑轴力很小，好象是墙体上的侧压力没有作用到支撑上。 建议作者先不要考虑water，或者简单设置水位就可以，慎用config fluid命令 我也对考虑水与不考虑水这俩种情况进行过比较，考虑水时，地表沉降小一些，不考虑水时，地表沉降大一些，但对shell单元上的弯矩和beam单元上的轴力影响不大。是不是先开挖再做shell连续墙会好点，有这样做成功的高人吗？请指点！地下连续墙用FLAC3D中的什么单元来模拟？到底如何模拟？用实体zone单元模拟行吗？我看到手册中是这样做的！ 刚开始试的时候可以把单元放少一点，我觉得太多了，运行速度太慢，而且对这么多单元来讲每个计算只有500步好象太少吧，运行后应该还不能达到稳定状态。 用了shell单元，最好不要用interface了吧？这样的话是不是直接用linear更好。 sel node free xr yr zr z y ran cid 10000 509 sel beamsel 10000 509 我想问一下，上面两行里的509是什么，怎么得到的？ 509是shell单元上的与beam对应的节点的ID号，它作为beam单元的终止节点，由于shell和beam单元共享这一节点，所以他们应该协同变形，共同受力！这一点可以查出来的！shell单元上的任意一节点的ID都可以查出来的！
geotechman
关于地下连续墙的模拟 1.地下连续墙用shell单元模拟是否合适？ 2.因为地下连续墙是在开挖以前就作好了，就应该在 model null 语句之 前用 sel shell id=1 range#### 语句先将其设好？ 3.shell单元是否只能设在边界面上或开挖的面上， 而不能先将shell单元设置在网格（土）里？ 可以先在计算区域外生成shell单元，然后通过ini node 命令将其置于实际位置。不过有时，先mode null再施加shell单元一样可以反应问题。 上述建议仅拱参考。
15．一个汇的小例子
一个汇的小例子;------------------------------------------------------ ; Excavation in a saturated soil
;------------------------------------------------------ config fluid
; --- geometrical model ---
gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1
21
group soil
group excavate range x 0 4 y 0 4 z 0 5 group wal1 range x 4 5 y 0 5 z 0 7 group wal2 range x 0 4 y 4 5 z 0 7
group wall range group wal1 any group wal2 any ; --- fluid flow model --- model fl_iso
prop perm 1e-12 poro 0.3 ini fdensity 1e3
ini fmod 2.0e9 ftens -1e-3
;model fl_null range group excavate ;model fl_null range group wall ini pp 0 grad 0 0 1e4 fix pp range z -0.1 0.1
fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ; --- mechanical model --- model elas
prop bul 3.9e6 shea 2.8e6
;model null range group excavate ini density 1.2e3
ini density 1.5e3 range group wall fix x range x -.1 .1 fix x range x 11.9 12.1 fix y range y -.1 .1 fix y range y 11.9 12.1 fix z range z 11.9 12.1
; initial total stresses ini szz 0 grad 0 0 -1.5e4 ini sxx 0 grad 0 0 -1.2e4 ini syy 0 grad 0 0 -1.2e4
set gravity 0 0 10 ; --- initial state --- plot cont pp
apply vwell -0.03e-3 range x 3.9 5.1 y 3.9 5.1 z 7.9 9.1 ;apply Pwell -0.3 GP 508 pause
solve force 1 ; check initial equilibrium
只是增加了一个vwell
还是可以看到降水曲线的，不过就是在背面， 奇怪怎么汇应力增大呢？
22
计算水，收敛起来，太慢了， 耦合，耦合就是慢，
就像人一样，一个人睡觉，很快就睡着了，
但是要2个人睡觉，那么就耦合半天，也不一定睡安稳， 还是1个人睡觉好哦。
16．用3DEC生成岩体随机节理网络
岩体内部裂隙分布在空间上具有随机性，岩体节理网络计算机模拟是近代计算机技术在岩石力学领域应用的一个重要方面，是研究岩体内大量节理裂隙随机分布规律的有效手段。 节理面的方位（或产状）可用节理面的产状要素即倾向和倾角表示或用节理面法向单位矢量的方向余弦来表示。根据现场实测统计可以获得一定地质区域具有不同产状的若干节理面，将这些节理面的法向矢量进行等面积赤平投影，可以得到用等面积网作出的节理极点图和等密度图。根据等密度图找出极高等值区即“极密区”，此“极密区”所代表的节理面即是概化优势节理面，从而确定优势产状并对节理面进行分组，最后得到具有优势产状的若干概化节理组。下图是某岩体节理倾向玫瑰花、节理极点和节理等密度图。
为便于进行三维节理网络随机模拟，首先需要生成节理随机数，这里采用直接法来生成随机数（引用他人成果）。该方法通过样本直接生成随机数，不需要建立样本分布函数，可以弥补传统方法中由于样本分散以及一些随机数不能很好地服从随机分布形式而带来误差的缺点。其具体过程如下：
(1) 将节理按产状进行分组；
23
(2) 对每一个节理参数的样本在其变化区间内进行离散化处理，即分成若干小区间，求每个小区间内该参数样本的频率（占样本总体的百分比）；
(3) 根据需要模拟区域的大小和节理间距的统计结果计算出在该区域范围内，在每一组节理的法向方向上需要生成节理的数目。由此可以计算出每个小区间内需要生成的节理随机数数目，作为目标数；
(4) 利用Mont-Carlo法生成一个在[0,1]范围内服从均匀分布的随机数，再把它换算到需要的区间范围内（与该节理参数统计样本的分布范围一致）；
(5) 设目前考察的是第i条待生成节理的某个参数的一个随机数，如果把它引入到它所属的小区间内，计算此时该区间内随机数目。如果此时的数目不超过目标数，保留该随机数，反之，生成并考察下一个随机数；
(6) 不断重复上述过程，直到生成全部随机数为止。 节理网络中的各个随机节理可以采用节理模拟程序产生，这里通过在Itasca公司基于离散元原理的3DEC软件中编写FISH函数来实现。根据节理来生成Block块，划分模型是3DEC建模的主要方法之一，因此，使用3DEC可以很方便地生成节理网络。 3DEC中生成节理的命令为：
jset dip * dd * origin *,*,*jset表示生成一条贯穿整个模型的节理，其参数dip表示节理的倾向，dd表示节理的倾角，origin表示节理面上一点的三维坐标。根据该命令并借助FISH函数生成随机数就可以很简单、快捷地生成随机节理网络。urand为FISH内部函数，其功能是生成0~1之间的均匀分布的随机数。a+(b-a)*urand则生成在区间[a,b]之间的均匀分布的随机数。每条节理的倾向和倾角由命令jset dip jdip dd jdd ori x 2 0获得，jdip和jdd分别是节理倾角dip和节理倾向dd的值，根据a+(b-a)*urand生成的随机数获得。遵照测线法，先指定一个节理初始点，下一节理面上一点的坐标根据节理初始点位置和节理间距获得，当前小区间的节理全部生成完后，进入下一小区间，直至生成整个节理网络。下图为生成的一组随机节理网络模型。
17．固结小算例
title
two dimension consolidation config fluid
24
FLAC3D中一些问题的讨论
