关于surfer绘制等直线，由于需要结合自己实际需要绘制，sufer只能绘制2D等直线，所以首先确定自己绘制等直线的面，然后flac3d导出的单元信息（坐标及应力，应变），通过简单的fortran程序找出所要求的面上的单元的，并记录其坐标位置，即应力应变值。输出格式为surfer认可的格式即可以了。
授人以鱼，不如授人以渔。我编的程序只能针对自己的问题，所以大家不要关心程序，如果，确实需要，按照上述思路实现应该很容易。
建立模型以后直接solve让程序产生应力分布!虽然时间花费上多一点，可是不用那么伤脑筋了。 我得懒人做法是：1）现把模量×10，c，Fi取很大，然后solve for balance；2）ini disp＝0，ini state ＝0，；3）回复原来得E，G，C，Fi等，
subzone和zone哪个是“最小几何单位”？二者有和区别呢？ 我们通常说的“单元element”指的是哪一个？ 用户接触到的最小几何单位应该是zone 通常说的element也是指zone
subzone和FLAC3d的 Mixed Discretization算法有关 比如说一个brick
计算时又自动划分为 2套overlay 每套overlay中各有5个tetra 即subzone
据说这样可以更好地模拟材料的塑性变形
这是flac3d为了避免求解中单元出现“沙漏”现象而采取的一种单元离散技术
Q:gen zone brick size 6 6 6 model elas
prop bulk 1e8 shear 7e7 fix x range x -0.1 0.1
table 1 0,0 100,-1e5 /我不知道table在这里什么意思,有什么作用 apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 hist zone sxx 6,0,0 step 100
我对这个table不是很清楚
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 gen zone reflect norm 1,0,0 gen zone reflect norm 0,0,1 model ss
prop bulk 1.19e10 shear 1.1e10 prop coh 2.72e5 fric 44 ten 2e5 prop ttab 1
table 1 0,2e5 2e-5,0
请问这里的table和后面的数字到底是什么意思啊? A:table相当于数已经赋值的数组
请问sel cable id=1 begin=(0.0, 0.5, 3.5) end=(8, 0.5, 3.5) nseg=8 中的nseg是什么意思? 把cable分为几段,比如nseg=8 ,就是分为八段.
gen zone reflect dip 90 dd 270 origin 0,0,0
应该是以90 和270度 还有原点组成的 平面进行对称拷贝吧？？
ini state＝0
初始化单元为弹性。一般来说，mc材料在初始化应力的时候，由于是手动设置的k0，因此局部单元可能已经进入了塑性状态，这会影响后继的计算步骤，因此需要ini state＝0
我这里编写了一个ansys to flac3d的前处理小程序，用起来比较方便。希望能和大家共享和交流。（可以自动形成flac开挖信息）
ansys-flac2002转换程序使用说明:
本程序用于将ansys 生成的数据格式转换为flac计算所需的数据格式。程序需要输入5个文件名。 1。单元信息文件，格式如下： 总单元数
单元信息。。。，例如： 10（总单元数）
1 1 1 1 0 1 8123 8713 8816 8812 8122 8826 8831 8827 2 1 1 1 0 1 8122 8826 8831 8827 8127 8687 7791 8456 3 1 1 1 0 1 8812 8816 8817 8813 8827 8831 8832 8828 4 1 1 1 0 1 8827 8831 8832 8828 8456 7791 7790 8455 5 1 1 1 0 1 8813 8817 8818 8814 8828 8832 8833 8829 6 1 1 1 0 1 8828 8832 8833 8829 8455 7790 7789 8454 7 1 1 1 0 1 8814 8818 8819 8815 8829 8833 8834 8830 8 1 1 1 0 1 8829 8833 8834 8830 8454 7789 7788 8453 9 1 1 1 0 1 8815 8819 7629 7632 8830 8834 7699 7704 10 1 1 1 0 1 8830 8834 7699 7704 8453 7788 7700 7703 2。节点信息 节点总数
节点信息。。。。例如： 20（节点总数）
1 210.750000000 228.398805249 -327.458594084 2 210.750000000 227.440685854 -299.440511342 3 210.750000000 226.118354064 -279.017920414 4 241.820265858 240.000000000 -331.304948353 5 246.293378743 240.000000000 -331.239755597 6 252.148198534 240.000000000 -331.154425346 7 258.555372995 240.000000000 -331.061044881 8 266.830116362 240.000000000 -330.940445794 9 274.869142734 240.000000000 -330.823282132 10 284.683736459 240.000000000 -330.680240712 11 295.556691364 240.000000000 -330.521774355 12 306.689058734 240.000000000 -330.359527228 13 318.078795786 240.000000000 -330.193529103 14 329.460877328 240.000000000 -330.027642551 15 340.835169125 240.000000000 -329.861869531
16 352.899895166 240.000000000 -329.686033873 17 367.026374896 240.000000000 -329.480149475 18 380.592277707 240.000000000 -329.282435127 19 390.569611461 240.000000000 -329.137021875 20 242.519404127 240.000000000 -307.394667740 3。输入开挖信息 开挖步数
第一步开挖单元数 开挖单元信息 第二步开挖单元数 开挖单元信息 。。。。。 例如
2（两步开挖）
5（第一步开挖5个单元）
12527 5 1 1 0 1 17639 17740 8430 8221 17891 17869 8457 8689 12528 5 1 1 0 1 17891 17869 8457 8689 17894 17872 8127 8126 12529 5 1 1 0 1 17638 17739 17740 17639 17890 17868 17869 17891 12530 5 1 1 0 1 17890 17868 17869 17891 17893 17871 17872 17894 12531 5 1 1 0 1 17637 17738 17739 17638 17889 17867 17868 17890 15（第二步开挖15个单元）
12532 5 1 1 0 1 17889 17867 17868 17890 17892 17870 17871 17893 12533 5 1 1 0 1 17641 17742 17738 17637 17887 17873 17867 17889 12534 5 1 1 0 1 17887 17873 17867 17889 17888 17874 17870 17892 12535 5 1 1 0 1 17643 17730 17742 17641 17885 17875 17873 17887 12536 5 1 1 0 1 17885 17875 17873 17887 17886 17876 17874 17888 12537 5 1 1 0 1 17645 17732 17730 17643 17883 17877 17875 17885 12538 5 1 1 0 1 17883 17877 17875 17885 17884 17878 17876 17886 12539 5 1 1 0 1 17647 17734 17732 17645 17881 17879 17877 17883 12540 5 1 1 0 1 17881 17879 17877 17883 17882 17880 17878 17884 12541 5 1 1 0 1 17633 17697 17734 17647 17832 17836 17879 17881 12542 5 1 1 0 1 17832 17836 17879 17881 17831 17835 17880 17882 12543 5 1 1 0 1 8383 8283 17349 17325 8410 8685 17559 17537 12544 5 1 1 0 1 8410 8685 17559 17537 8411 8207 17562 17540 12545 5 1 1 0 1 17325 17349 17348 17324 17537 17559 17558 17536 12546 5 1 1 0 1 17537 17559 17558 17536 17540 17562 17561 17539 4。输入计算结果文件名：该文件存放各种几何信息，包括材料分组信息 5。输入开挖信息文件名：该文件存放开挖分组信息。
所有材料的分组和开挖分组（group）都是自动定义的 在ansys里将要分步开挖的单元分别找出来，按使用说明输入即可对不同开挖步信息分别给出它的group名
本code最大单元和结点数量可达1000000个。转换后的文件可能是大了点，不过很好用，转换后的文件直接call就行了。
将开挖信息列出来,是为了给它们一个group name,以便再开挖计算时使用。在ansys里列开挖单元是不难
的。
第一步，在ansys中建立三维模型，并给每个单元赋予相应的材料编号 第二步，将ansys中的单元信息和结点信息导出，格式如上述例子 第三步，在ansys中找出每一步开挖的单元号信息，格式如上 第四步，执行转换程序即可
材料分组和开挖单元分组由转换程序自动完成。
该程序绝对能用，如果出现问题，检查输入文件是否正确
在ＦＬＡＣ现有的模型，包括摩尔库伦和应力应变软化模型，运用到隧道等开挖支护的时候，并不能考虑到由于岩层土体应力变化对其材料参数的影响，很多情况下，没有锚杆支护与使用锚杆支护的效果并不明显，与实际出入较大．这是模型的缺陷．如果中间调整土体岩层材料的参数后，效果要明显许多，但参数的重新选取只能是经验来决定，误差较大．另外，锚杆支护时有个初锚力，然后再逐渐加载到一定值．模拟时，可以按此分步实现，但发现效果也并不明显．如果用ANSYS模拟锚杆支护，只能对杆端起到作用，不能准确模拟整个杆的受力状况．具体可以查阅中国矿业大学出版社＜锚杆支护的数值模拟＞． 一个煤矿巷道的例子： new
;100m*100m gr 120 100 ;ROCK m m j=1,30 m m j=32,43 m m j=45,56 m m j=57,60 m m j=62,65 m m j=67,76 m m j=78,79 m m j=81,84 m m j=86,100 ;INTERFACE m n j=31 m n j=44 m n j=61 m n j=66 m n j=77 m n j=80 m n j=85 ;40m
gen 0.00, -32.50 0.0,30.00 40.00,17.00 40.00,-32.50 i=1,41 j=1,31 gen 0.00, 30.00 0.0,38.00 40.00,25.00 40.00,17.00 i=1,41 j=32,44 gen 0.00, 38.00 0.0,44.00 40.00,31.00 40.00,25.00 i=1,41 j=45,57 gen 0.00, 44.00 0.0,46.00 40.00,33.00 40.00,31.00 i=1,41 j=57,61 gen 0.00, 46.00 0.0,48.00 40.00,35.00 40.00,33.00 i=1,41 j=62,66
gen 0.00, 48.00 0.0,55.50 40.00,42.50 40.00,35.00 i=1,41 j=67,77 gen 0.00, 55.50 0.0,57.30 40.00,44.30 40.00,42.50 i=1,41 j=78,80 gen 0.00, 57.30 0.0,61.30 40.00,48.30 40.00,44.30 i=1,41 j=81,85 gen 0.00, 61.30 0.0,67.50 40.00,67.50 40.00,48.30 i=1,41 j=86,101 ;20m
gen 40.00, -32.50 40.0,17.00 60.00,10.50 60.00,-32.50 i=41,81 j=1,31 gen 40.00, 17.00 40.0,25.00 60.00,18.50 60.00,10.50 i=41,81 j=32,44
gen 40.00, 25.00 40.0,27.50 47.90,24.7675 47.90,22.4325 i=41,57 j=45,50 gen 40.00, 27.50 40.0,31.00 47.90,28.4325 47.90,24.7675 i=41,57 j=50,57
gen 47.90,22.4325 47.90,24.7675 52.1,24.7675 52.1,21.0675 i=57,65 j=45,50 gen 47.90,24.7675 47.90,28.4325 52.1,27.0675 52.1,24.7675 i=57,65 j=50,57
gen 52.1,21.0675 52.1,24.7675 60,21 60,18.5 i=65,81 j=45,50 gen 52.1,24.7675 52.1,27.0675 60,24.5 60,21 i=65,81 j=50,57
gen 40.00, 31.00 40.0,33.00 60.00,26.50 60.00,24.50 i=41,81 j=57,61 gen 40.00, 33.00 40.0,35.00 60.00,28.50 60.00,26.50 i=41,81 j=62,66 gen 40.00, 35.00 40.0,42.50 60.00,36.00 60.00,28.50 i=41,81 j=67,77 gen 40.00, 42.50 40.0,44.30 60.00,37.80 60.00,36.00 i=41,81 j=78,80 gen 40.00, 44.30 40.0,48.30 60.00,41.80 60.00,37.80 i=41,81 j=81,85 gen 40.00, 48.30 40.0,67.50 60.00,67.50 60.00,41.80 i=41,81 j=86,101 ;40m
gen 60.00, -32.50 60.0,10.50 100.00,-2.50 100.00,-32.50 i=81,121 j=1,31 gen 60.00, 10.50 60.0,18.50 100.00,5.50 100.00,-2.50 i=81,121 j=32,44 gen 60.00, 18.50 60.0,24.50 100.00,11.50 100.00,5.50 i=81,121 j=45,57 gen 60.00, 24.50 60.0,26.50 100.00,13.50 100.00,11.50 i=81,121 j=57,61 gen 60.00, 26.50 60.0,28.50 100.00,15.50 100.00,13.50 i=81,121 j=62,66 gen 60.00, 28.50 60.0,36.00 100.00,23.00 100.00,15.50 i=81,121 j=67,77 gen 60.00, 36.00 60.0,37.80 100.00,24.80 100.00,23.00 i=81,121 j=78,80 gen 60.00, 37.80 60.0,41.80 100.00,28.80 100.00,24.80 i=81,121 j=81,85 gen 60.00, 41.80 60.0,67.50 100.00,67.50 100.00,28.80 i=81,121 j=86,101
;INTERFACE
interface 1 aside from 1,31 to 121,31 bside from 1,32 to 121,32 interface 1 c= 0.2e+6 f=28 kn=588.00e+6 ks=14.70e+6 interface 2 aside from 1,44 to 121,44 bside from 1,45 to 121,45 interface 2 c= 0.2e+6 f=28 kn=588.00e+6 ks= 14.70e+6 ;interface 3 aside from 1,57 to 121,57 bside from 1,58 to 121,58 ;interface 3 c= 0.05e+6 f=20 kn=480.00e+6 ks= 9.00e+6 interface 4 aside from 1,61 to 121,61 bside from 1,62 to 121,62 interface 4 c= 0.1e+6 f=25 kn=588.00e+6 ks= 14.70e+6