m.elements[5] = Tri3(conn = (6,7,10))
m.elements[6] = Quad4(conn = (5,6,10,9))
m = m.extrude(translation = [0, 0, 1], layers = 1)
def transformation(x, y, z):
theta = np.pi / 2. * z
r = x.copy()
x = r * np.cos(theta)
z = r * np.sin(theta)
return x, y, z
m = m.transform(transformation)
mapping = m.overlapping_nodes()
m = m.merge_nodes(mapping)
m2 = m.copy()
"""
m2.elements.clear()
counter = 1
for label in m.elements.keys():
el = m.elements[label]
simplices = el.simplex_decomposition()
for simplex in simplices:
m2.elements[cunter] = simplex
counter +=1
"""
m2 = m2.clean_connectivity()
writeInp(m2, path = "clean_connectivity.inp", mapping = {"Tri3": "CAX3", "Quad4": "CAX4", "Hexa8": "C3D8", "Prism6":"C3D6", "Tetra4": "C3D4"})
writeMsh(m2, path = "clean_connectivity.msh")
m = m.extrude(translation=[0, 0, 1], layers=1)
def transformation(x, y, z):
theta = np.pi / 2. * z
r = x.copy()
x = r * np.cos(theta)
z = r * np.sin(theta)
return x, y, z
m = m.transform(transformation)
mapping = m.overlapping_nodes()
m = m.merge_nodes(mapping)
m2 = m.copy()
"""
m2.elements.clear()
counter = 1
for label in m.elements.keys():
el = m.elements[label]
simplices = el.simplex_decomposition()
for simplex in simplices:
m2.elements[cunter] = simplex
counter +=1
"""
m2 = m2.clean_connectivity()
writeInp(m2,
path="clean_connectivity.inp",
mapping={
