suzanne.soy/FreeCAD
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

FEM: ccxInpWriter, copy Nodes dict

Browse Source
This commit is contained in:
fandaL 2016-02-21 16:27:09 +01:00 committed by Bernd Hahnebach
parent 816b54bd25
commit c1bf26bbd7
1 changed files with 16 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
src/Mod/Fem/ccxInpWriter.py
View File

@ -835,14 +835,15 @@ class inp_writer:
# [ (nodeID, length), ... , (nodeID, length) ] some nodes will have more than one entry # [ (nodeID, length), ... , (nodeID, length) ] some nodes will have more than one entry
node_length_table = [] node_length_table = []
mesh_edge_length = 0 mesh_edge_length = 0
mesh_nodes = self.mesh_object.FemMesh.Nodes
# print(len(edge_table)) # print(len(edge_table))
for me in edge_table: for me in edge_table:
if len(edge_table[me]) == 2: # 2 node mesh edge if len(edge_table[me]) == 2: # 2 node mesh edge
# end_node_length = mesh_edge_length / 2 # end_node_length = mesh_edge_length / 2
# ______ # ______
# P1 P2 # P1 P2
P1 = self.mesh_object.FemMesh.Nodes[edge_table[me][0]] P1 = mesh_nodes[edge_table[me][0]]
P2 = self.mesh_object.FemMesh.Nodes[edge_table[me][1]] P2 = mesh_nodes[edge_table[me][1]]
edge_vec = P2 - P1 edge_vec = P2 - P1
mesh_edge_length = edge_vec.Length mesh_edge_length = edge_vec.Length
# print(mesh_edge_length) # print(mesh_edge_length)
@ -855,9 +856,9 @@ class inp_writer:
# middle_node_length = mesh_face_area * 2 / 3 # middle_node_length = mesh_face_area * 2 / 3
# _______ _______ # _______ _______
# P1 P3 P2 # P1 P3 P2
P1 = self.mesh_object.FemMesh.Nodes[edge_table[me][0]] P1 = mesh_nodes[edge_table[me][0]]
P2 = self.mesh_object.FemMesh.Nodes[edge_table[me][1]] P2 = mesh_nodes[edge_table[me][1]]
P3 = self.mesh_object.FemMesh.Nodes[edge_table[me][2]] P3 = mesh_nodes[edge_table[me][2]]
edge_vec1 = P3 - P1 edge_vec1 = P3 - P1
edge_vec2 = P2 - P3 edge_vec2 = P2 - P3
mesh_edge_length = edge_vec1.Length + edge_vec2.Length mesh_edge_length = edge_vec1.Length + edge_vec2.Length
@ -879,6 +880,7 @@ class inp_writer:
# [ (nodeID,Area), ... , (nodeID,Area) ] some nodes will have more than one entry # [ (nodeID,Area), ... , (nodeID,Area) ] some nodes will have more than one entry
node_area_table = [] node_area_table = []
mesh_face_area = 0 mesh_face_area = 0
mesh_nodes = self.mesh_object.FemMesh.Nodes
for mf in face_table: for mf in face_table:
if len(face_table[mf]) == 3: # 3 node mesh face triangle if len(face_table[mf]) == 3: # 3 node mesh face triangle
# corner_node_area = mesh_face_area / 3.0 # corner_node_area = mesh_face_area / 3.0
@ -887,9 +889,9 @@ class inp_writer:
# / \ # / \
# /____\ # /____\
# P1 P2 # P1 P2
P1 = self.mesh_object.FemMesh.Nodes[face_table[mf][0]] P1 = mesh_nodes[face_table[mf][0]]
P2 = self.mesh_object.FemMesh.Nodes[face_table[mf][1]] P2 = mesh_nodes[face_table[mf][1]]
P3 = self.mesh_object.FemMesh.Nodes[face_table[mf][2]] P3 = mesh_nodes[face_table[mf][2]]
mesh_face_area = getTriangleArea(P1, P2, P3) mesh_face_area = getTriangleArea(P1, P2, P3)
corner_node_area = mesh_face_area / 3.0 corner_node_area = mesh_face_area / 3.0
@ -913,12 +915,12 @@ class inp_writer:
# /t1 \ /t2 \ # /t1 \ /t2 \
# /_____\/_____\ # /_____\/_____\
# P1 P4 P2 # P1 P4 P2
P1 = self.mesh_object.FemMesh.Nodes[face_table[mf][0]] P1 = mesh_nodes[face_table[mf][0]]
P2 = self.mesh_object.FemMesh.Nodes[face_table[mf][1]] P2 = mesh_nodes[face_table[mf][1]]
P3 = self.mesh_object.FemMesh.Nodes[face_table[mf][2]] P3 = mesh_nodes[face_table[mf][2]]
P4 = self.mesh_object.FemMesh.Nodes[face_table[mf][3]] P4 = mesh_nodes[face_table[mf][3]]
P5 = self.mesh_object.FemMesh.Nodes[face_table[mf][4]] P5 = mesh_nodes[face_table[mf][4]]
P6 = self.mesh_object.FemMesh.Nodes[face_table[mf][5]] P6 = mesh_nodes[face_table[mf][5]]
mesh_face_t1_area = getTriangleArea(P1, P4, P6) mesh_face_t1_area = getTriangleArea(P1, P4, P6)
mesh_face_t2_area = getTriangleArea(P2, P5, P4) mesh_face_t2_area = getTriangleArea(P2, P5, P4)