From cd4952dee4858fae0cb312f3cd40ad3aec9dbfb5 Mon Sep 17 00:00:00 2001 From: Przemo Firszt Date: Tue, 26 May 2015 14:20:18 +0100 Subject: [PATCH] FEM: fix white spaces in ccxInpWriter.py Signed-off-by: Przemo Firszt --- src/Mod/Fem/ccxInpWriter.py | 33 +++++++++++++++++++-------------- 1 file changed, 19 insertions(+), 14 deletions(-) diff --git a/src/Mod/Fem/ccxInpWriter.py b/src/Mod/Fem/ccxInpWriter.py index a902cd852..4b427eb24 100644 --- a/src/Mod/Fem/ccxInpWriter.py +++ b/src/Mod/Fem/ccxInpWriter.py @@ -156,11 +156,13 @@ class inp_writer: f.write(fix_obj_name + ',3\n\n') def write_constraints_force(self, f): - def getTriangleArea(P1,P2,P3): + + def getTriangleArea(P1, P2, P3): vec1 = P2 - P1 vec2 = P3 - P1 vec3 = vec1.cross(vec2) return 0.5 * vec3.Length + f.write('\n***********************************************************\n') f.write('** Node loads\n') f.write('** written by {} function\n'.format(sys._getframe().f_code.co_name)) @@ -202,14 +204,17 @@ class inp_writer: f.write('** node loads on element face: ' + o.Name + '.' + elem + '\n') volume_faces = self.mesh_object.FemMesh.getVolumesByFace(ref_face) - face_table = {} # { meshfaceID : ( nodeID, ... , nodeID ) } - for mv,mf in volume_faces: + face_table = {} # { meshfaceID : ( nodeID, ... , nodeID ) } + for mv, mf in volume_faces: face_table[mf] = self.mesh_object.FemMesh.getElementNodes(mf) # calulate the appropriate node_areas for every node of every mesh face (mf) # G. Lakshmi Narasaiah, Finite Element Analysis, p206ff - node_area_table = [] # [ (nodeID,Area), ... , (nodeID,Area) ] some nodes will have more than one entries - node_sumarea_table = {} # { nodeID : Area, ... , nodeID:Area } AreaSum for each node, one entry for each node + + # [ (nodeID,Area), ... , (nodeID,Area) ] some nodes will have more than one entry + node_area_table = [] + # { nodeID : Area, ... , nodeID:Area } AreaSum for each node, one entry for each node + node_sumarea_table = {} mesh_face_area = 0 for mf in face_table: # print ' ', mf, ' --> ', face_table[mf] @@ -224,7 +229,7 @@ class inp_writer: P2 = self.mesh_object.FemMesh.Nodes[face_table[mf][1]] P3 = self.mesh_object.FemMesh.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 node_area_table.append((face_table[mf][0], corner_node_area)) @@ -250,11 +255,11 @@ class inp_writer: P5 = self.mesh_object.FemMesh.Nodes[face_table[mf][4]] P6 = self.mesh_object.FemMesh.Nodes[face_table[mf][5]] - mesh_face_t1_area = getTriangleArea(P1,P4,P6) - mesh_face_t2_area = getTriangleArea(P2,P5,P4) - mesh_face_t3_area = getTriangleArea(P3,P6,P5) - mesh_face_t4_area = getTriangleArea(P4,P5,P6) - mesh_face_area = mesh_face_t1_area + mesh_face_t2_area + mesh_face_t3_area + mesh_face_t4_area + mesh_face_t1_area = getTriangleArea(P1, P4, P6) + mesh_face_t2_area = getTriangleArea(P2, P5, P4) + mesh_face_t3_area = getTriangleArea(P3, P6, P5) + mesh_face_t4_area = getTriangleArea(P4, P5, P6) + mesh_face_area = mesh_face_t1_area + mesh_face_t2_area + mesh_face_t3_area + mesh_face_t4_area middle_node_area = mesh_face_area / 3.0 node_area_table.append((face_table[mf][0], 0)) @@ -269,12 +274,12 @@ class inp_writer: # print n, ' --> ', A if n in node_sumarea_table: node_sumarea_table[n] = node_sumarea_table[n] + A - else: + else: node_sumarea_table[n] = A sum_node_areas = 0 for n in node_sumarea_table: - # print n, ' --> ', node_sumarea_table[n] + # print n, ' --> ', node_sumarea_table[n] sum_node_areas = sum_node_areas + node_sumarea_table[n] print ' sum_node_areas ', sum_node_areas, ' ref_face.Area: ', ref_face.Area sum_ref_face_node_area += sum_node_areas @@ -300,7 +305,7 @@ class inp_writer: # print ' sum_ref_face_area : ', sum_ref_face_area # print ' sum_ref_face_node_area * force_per_sum_ref_face_area: ', sum_ref_face_node_area * force_per_sum_ref_face_area # print ' sum_node_load: ', sum_node_load - # print ' frc_obj.Force: ', frc_obj.Force + # print ' frc_obj.Force: ', frc_obj.Force f.write('\n') def write_face_load(self, f):