def get_element_properties(self, eid):
"""
Returns area, centroid, normal
"""
nodes = self.get_element_nodes(eid)
(area, centroid, normal) = triangle_area_centroid_normal(nodes)
return (area, centroid, normal)
def get_element_properties(self, eid):
"""
Returns area, centroid, normal
"""
nodes = self.get_element_nodes(eid)
(area, centroid, normal) = triangle_area_centroid_normal(nodes)
return (area, centroid, normal)
def prepare_centroid_area_normals(self):
for eid, element in enumerate(self.elements):
n1, n2, n3 = element
n1 = self.nodes[n1 - 1]
n2 = self.nodes[n2 - 1]
n3 = self.nodes[n3 - 1]
area, centroid, normal = triangle_area_centroid_normal([n1, n2, n3])
eidi = eid + 1
if centroid is not None:
if len(centroid) != 3:
msg = "eid=%s centroid=%s n1=%s n2=%s n3=%s" % (eidi, centroid, n1, n2, n3)
raise RuntimeError(msg)
self.areas[eidi] = area
self.centroids[eidi] = centroid
self.normals[eidi] = normal
def prepare_centroid_area_normals(self):
for eid, element in enumerate(self.elements):
n1, n2, n3 = element
n1 = self.nodes[n1 - 1]
n2 = self.nodes[n2 - 1]
n3 = self.nodes[n3 - 1]
area, centroid, normal = triangle_area_centroid_normal([n1, n2, n3])
eidi = eid + 1
if centroid is not None:
if len(centroid) != 3:
msg = "eid=%s centroid=%s n1=%s n2=%s n3=%s" % (eidi, centroid, n1, n2, n3)
raise RuntimeError(msg)
self.areas[eidi] = area
self.centroids[eidi] = centroid
self.normals[eidi] = normal
