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

FEM: fix white spaces in ccxInpWriter.py

Browse Source 
Signed-off-by: Przemo Firszt <przemo@firszt.eu>
This commit is contained in:
Przemo Firszt 2015-05-26 14:20:18 +01:00
parent e13a8cf304
commit cd4952dee4
1 changed files with 19 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -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):