-
Notifications
You must be signed in to change notification settings - Fork 0
/
INP_Circle_returnverticesandpolygons.py
52 lines (41 loc) · 1.87 KB
/
INP_Circle_returnverticesandpolygons.py
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
#coding=utf8
#################################################################
### ###
### Created by Bastien MARTY, 2017 ###
### ###
### École polytechnique, Palaiseau, France ###
### ###
#################################################################
### Read Abaqus File from a CIRCLE and return the vertices and polygons
def INP_Circle_returnverticesandpolygons(
mesh_filename,
verbose=0):
points = []
polygons = []
mesh_file = open(mesh_filename, "r")
context = ""
for line in mesh_file:
if line.startswith("**") : continue
line = line.strip("\n ,")
#mypy.my_print(verbose-1, "line =", line)
if (context != ""):
if line.startswith("*"):
context = ""
else:
if (context == "reading nodes"):
splitted_line = line.split(",")
points.append([float(coord) for coord in splitted_line[1:4]])
elif (context == "reading elems"):
splitted_line = line.split(",")
polygons.append([int(coord)-2 for coord in splitted_line[1:5]])
if line.startswith("1, "):
context = "reading nodes"
elif line.startswith("*ELEMENT"):
if (("TYPE=CPS4" in line) or ("type=CPS4" in line) or ("TYPE=CPE4" in line) or ("type=CPE4" in line)):
context = "reading elems"
else:
print("Warning: elements not read: "+line+".")
#we destroy the z coordinate that equals 0 evrywhere
for i in range(len(points)):
points[i].remove(points[i][-1])
return [points,polygons]