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,...

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！