kindred/create
1
1
Fork 0
Code Issues 37 Pull Requests Actions Packages Projects 9 Releases 2 Wiki Activity

FEM: mesh package, all left code, use single quotes instead of double ones

Browse Source
This commit is contained in:
Bernd Hahnebach
2019-09-05 21:40:59 +02:00
parent 12e92fef31
commit 9217a08601
1 changed files with 13 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
src/Mod/Fem/femmesh/femmesh2mesh.py
View File

@@ -31,7 +31,7 @@ import time
# import Mesh
'''
"""
load FreeCADs 3D FEM example from Start Workbench
femmesh = App.ActiveDocument.getObject("Box_Mesh").FemMesh
@@ -41,7 +41,7 @@ out_mesh = femmesh.femmesh2mesh.femmesh_2_mesh(femmesh, result)
import Mesh
Mesh.show(Mesh.Mesh(out_mesh))
'''
"""
# These dictionaries list the nodes, that define faces of an element.
# The key is the face number, used internally by FreeCAD.
# The list contains the nodes in the element for each face.
@@ -97,7 +97,7 @@ def femmesh_2_mesh(myFemMesh, myResults=None):
if myFemMesh.VolumeCount > 0:
for ele in myFemMesh.Volumes:
element_nodes = myFemMesh.getElementNodes(ele)
# print('element_node: ', element_nodes)
# print("element_node: ", element_nodes)
faceDef = face_dicts[len(element_nodes)]
for key in faceDef:
@@ -113,13 +113,13 @@ def femmesh_2_mesh(myFemMesh, myResults=None):
faceCode += (node << shifter)
# x << n: x shifted left by n bits = Multiplication
shifter += shiftBits
# print('codeList: ', codeList)
# print("codeList: ", codeList)
faceCodeDict[faceCode] = nodeList
faceCodeList.append(faceCode)
elif myFemMesh.FaceCount > 0:
for ele in myFemMesh.Faces:
element_nodes = myFemMesh.getElementNodes(ele)
# print('element_node: ', element_nodes)
# print("element_node: ", element_nodes)
faceDef = {1: [0, 1, 2]}
for key in faceDef:
@@ -135,7 +135,7 @@ def femmesh_2_mesh(myFemMesh, myResults=None):
faceCode += (node << shifter)
# x << n: x shifted left by n bits = Multiplication
shifter += shiftBits
# print('codeList: ', codeList)
# print("codeList: ", codeList)
faceCodeDict[faceCode] = nodeList
faceCodeList.append(faceCode)
@@ -150,11 +150,11 @@ def femmesh_2_mesh(myFemMesh, myResults=None):
if faceCodeList[actFaceIdx] == faceCodeList[actFaceIdx + 1]:
actFaceIdx += 2
else:
# print('found a single Face: ', faceCodeList[actFaceIdx])
# print("found a single Face: ", faceCodeList[actFaceIdx])
singleFaces.append(faceCodeList[actFaceIdx])
actFaceIdx += 1
else:
print('found a last Face: ', faceCodeList[actFaceIdx])
print("found a last Face: ", faceCodeList[actFaceIdx])
singleFaces.append(faceCodeList[actFaceIdx])
actFaceIdx += 1
@@ -170,14 +170,14 @@ def femmesh_2_mesh(myFemMesh, myResults=None):
myFemMesh.getNodeById(face_nodes[1]) + dispVec1,
myFemMesh.getNodeById(face_nodes[2]) + dispVec2]
output_mesh.extend(triangle)
# print('my triangle: ', triangle)
# print("my triangle: ", triangle)
if len(face_nodes) == 4:
dispVec3 = myResults.DisplacementVectors[myResults.NodeNumbers.index(face_nodes[3])]
triangle = [myFemMesh.getNodeById(face_nodes[2]) + dispVec2,
myFemMesh.getNodeById(face_nodes[3]) + dispVec3,
myFemMesh.getNodeById(face_nodes[0]) + dispVec0]
output_mesh.extend(triangle)
# print('my 2. triangle: ', triangle)
# print("my 2. triangle: ", triangle)
else:
for myFace in singleFaces:
@@ -186,14 +186,14 @@ def femmesh_2_mesh(myFemMesh, myResults=None):
myFemMesh.getNodeById(face_nodes[1]),
myFemMesh.getNodeById(face_nodes[2])]
output_mesh.extend(triangle)
# print('my triangle: ', triangle)
# print("my triangle: ", triangle)
if len(face_nodes) == 4:
triangle = [myFemMesh.getNodeById(face_nodes[2]),
myFemMesh.getNodeById(face_nodes[3]),
myFemMesh.getNodeById(face_nodes[0])]
output_mesh.extend(triangle)
# print('my 2. triangle: ', triangle)
# print("my 2. triangle: ", triangle)
end_time = time.clock()
print('Mesh by surface search method: ', end_time - start_time)
print("Mesh by surface search method: ", end_time - start_time)
return output_mesh