Abaqus obj inp 文件内部读取点坐标并统计

任务描述：对一个三维的带网格的模型，生成一个二维平面的图像描述，其中每个空格内部显示这个坐标下在z方向具有多少网格，类似于三维模型的深度图
方法：有两个方法
1.从网格文件OBJ内部读取点坐标，统计（x,y）的个数
2.从Inp文件内部读取点坐标，统计（x,y）的个数
在实验过程中，发现方法1 的obj文件内部的点的信息比较杂，可能包含除了网格以外的其他的点的信息，所以放弃了这种方法。这里主要介绍方法2

inp文件读取

writenode = []
#writenode 点集
write_flag = False
input = open("new_1.inp","r")
#读取inp文件
astring = input.readlines()
#这里这种读入方式的好处是astring是list类型，其余方式可能是str类型
for i,line in enumerate(astring):
    if line =="*Node\n":
    #找到inp文件内部开始记录点信息的地方
        write_flag = True
    #写入文件标志位开始
        continue
    #跳出当前循环，执行下一次循环，这样*Node就不会写入文件内部
        print(i)
    if line =="*Element, type=C3D8R\n":
        break
    #点描述信息结束后是描述单元信息，正好这里截止，挑出for循环
    if write_flag ==True:
       writenode.append(line.replace(",",""))

点集写出

接下来可以将我们的点信息输出另一个文件：

f=open('point.txt','w') 
for i in writenode:
    f.write(i)
f.close()

重新读入生成numpy array

a = np.loadtxt('zz.txt') 
a = np.around(a)
#四舍五入，可能出现一些不准的点，防止出问题 

删除无关信息

c= np.delete(a,0,axis=1)
#删除第一列：各个点的标号[1,2,3,4,5...]
c= np.delete(c,2,axis=1)
#删除Z坐标，只需要统计xy坐标

统计xy相同的点的个数

e = []
for (x,y) in (c):
#这里如果直接统计列表是会报错的，无法生成hash表所以要先生成这样的数组的类型，后面才好统计
    e.append((x,y))

然后做一个字典，记录各个点的个数

h = {}
for i in e:
    h[i] = e.count(i)
print()

这样就可以看到各个点的坐标的个数了

{(-139.0, -18.0): 23,
 (199.0, -96.0): 23,
 (-199.0, -60.0): 23,
 (-182.0, -150.0): 23,
 (251.0, -60.0): 23,
 (-243.0, -102.0): 23,
 (208.0, 66.0): 23,...