1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61

我们都知道，自从flac3d500版本加入了自身建模功能后，可以很方便的通过描图、建块、划分、拉伸等步骤建立起适合flac3d分析的模型，而flac3d300版本软件不能直接读入500版本生成的模型，那么怎样将500建立的模型读入300版本中且保留组（group）的信息呢？
目的：将flac3d500版本建模信息完整的导入flac3d300中
原因：部分同学不会用flac3d500高级命令，但500版本建模易学；flac3d300不能直接读入500版本的建模信息
首先让我们分析一下flac3d300中模型是按怎样的格式读入读出的。
new
gen zon bri size 3 3 3
group roof range z 2 3 
group floor range z 0 1
expgrid flac3d300modelgroup.flac3d

;观察生成的带组的节点文档内容，可见组成如下：
* FLAC3D grid produced by FLAC3D
* GRIDPOINTS
G 1 0.000000000e+000 0.000000000e+000 0.000000000e+000
G 2 1.000000000e+000 0.000000000e+000 0.000000000e+000
G 3 0.000000000e+000 1.000000000e+000 0.000000000e+000
......
* ZONES
Z B8 1 1 2 3 4 5 6 7 8
Z B8 2 2 9 5 7 10 8 11 12
Z B8 3 9 13 10 11 14 12 15 16
......
* GROUPS
ZGROUP roof
 19 20 21 22 23 24 25 26 27
ZGROUP floor
 1 2 3 4 5 6 7 8 9
通过上面可以看出有gridpoint、zones和zgroup三部分组成，其中省略号是我添加上代表还有很多信息没贴上。
那flac3d500是怎样的输出格式呢？建立模型：
gen zon bri size 3 3 3
group roof range z 2 3 
group floor range z 0 1
plot zone
expgrid flac3d300modelgroup.flac3d

观察500版本输出的模型信息如下：

new
* FLAC3D grid produced by FLAC3D 5.00 Release 093
* 周四 四月 23 13:51:42 2015
* GRIDPOINTS
G 1   0.00000000000000e+00   0.00000000000000e+00   0.00000000000000e+00
G 2   1.00000000000000e+00   0.00000000000000e+00   0.00000000000000e+00
G 3   0.00000000000000e+00   1.00000000000000e+00   0.00000000000000e+00
......
* ZONES
Z B8 1 1 2 3 4 5 6 7 8
Z B8 2 2 9 5 7 10 8 11 12
Z B8 3 9 13 10 11 14 12 15 16
......
* GROUPS
ZGROUP 'floor' SLOT 1
 1 2 3 4 5 6 7 8 9
ZGROUP 'roof' SLOT 1
 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
 25 26 27

那么将flac33d500建立的模型输出为*.flac3d后，把格式改为完全相同即可。
做法就是把相应的题头改为相同即可。把组名的单引号和后面的SLOT 1去掉，比如 'floor' SLOT 1改为floor即可，余同。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

# -*- coding: cp936 -*-
"""
Created on Sun Dec 27 08:00:44 2015
@author: Zhubinglong
"""
import os
os.chdir('D:\\')       # 跳到D盘
if not os.path.exists('test1.txt'): # 看一下这个文件是否存在
    exit(-1)  
else:                       #不存在就退出
    lines = open('test1.txt').readlines()  #打开文件，读入每一行
    fp = open('test2.txt','w')  #打开你要写得文件test2.txt
    for s in lines:
    # replace是替换，write是写入
        fp.write( s.replace("'",'').replace('SLOT 1',''))    
    fp.close()  # 关闭文件
print 'end'

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

# -*- coding: cp936 -*-
"""
Created on Sun Dec 27 08:00:44 2015

@author: Administrator
"""
import os
os.chdir(os.getcwd())       # 跳到D盘
if not os.path.exists('test1.txt'): # 看一下这个文件是否存在
    exit(-1)  
else:                       #不存在就退出
    lines = open('test1.txt').readlines()  #打开文件，读入每一行
    fp = open('test2.txt','w')  #打开你要写得文件test2.txt
    for s in lines:
    # replace是替换，write是写入
        fp.write( s.replace("'",'').replace('SLOT 1',''))    
    fp.close()  # 关闭文件
print 'end'