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

2023-07-16 12:32| 来源: 网络整理| 查看: 265

任务描述：对一个三维的带网格的模型，生成一个二维平面的图像描述，其中每个空格内部显示这个坐标下在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,...

【本文地址】

公司简介

联系我们