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

FEM: constraint pressure, add support for shell meshes if they have group data

Browse Source
This commit is contained in:
Bernd Hahnebach
2017-05-26 21:41:09 +01:00
committed by Yorik van Havre
parent 34cac07b43
commit 7ed786cf86
3 changed files with 42 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/Mod/Fem/FemInputWriter.py
View File

@@ -175,6 +175,6 @@ class FemInputWriter():
pressure_faces = FemMeshTools.get_pressure_obj_faces(self.femmesh, self.femelement_table, self.femnodes_ele_table, femobj)
# print(len(pressure_faces))
femobj['PressureFaces'] = [(femobj['Object'].Name + ': face load', pressure_faces)]
# print(femobj['PressureFaces'])
print(femobj['PressureFaces'])
## @}

7
src/Mod/Fem/FemInputWriterCcx.py
View File

@@ -897,7 +897,12 @@ class FemInputWriterCcx(FemInputWriter.FemInputWriter):
for ref_shape in femobj['PressureFaces']:
f.write('** ' + ref_shape[0] + '\n')
for face, fno in ref_shape[1]:
f.write("{},P{},{}\n".format(face, fno, rev * prs_obj.Pressure))
if fno > 0: # solid mesh face
f.write("{},P{},{}\n".format(face, fno, rev * prs_obj.Pressure))
elif fno == 0: # on shell mesh face: fno == 0 --> normal of element face == face normal
f.write("{},P,{}\n".format(face, rev * prs_obj.Pressure))
elif fno == -1: # on shell mesh face: fno == -1 --> normal of element face oposite direction face normal
f.write("{},P,{}\n".format(face, -1 * rev * prs_obj.Pressure))
def write_constraints_temperature(self, f):
f.write('\n***********************************************************\n')

41
src/Mod/Fem/FemMeshTools.py
View File

@@ -373,6 +373,20 @@ def get_femnode_set_from_group_data(femmesh, fem_object):
return group_nodes
def get_femelementface_sets_from_group_data(femmesh, fem_object):
# get femfaceelements from femmesh face groupdata for reference shapes of obj.References
obj = fem_object['Object']
group_faces = None
if femmesh.GroupCount:
for g in femmesh.Groups:
grp_name = femmesh.getGroupName(g)
if grp_name.startswith(obj.Name + "_"):
if femmesh.getGroupElementType(g) == "Face":
print("Constraint: " + obj.Name + " --> " + "mesh group: " + grp_name)
group_faces = femmesh.getGroupElements(g) # == ref_shape_femelements
return group_faces
def get_femelement_sets_from_group_data(femmesh, fem_objects):
# get femelements from femmesh groupdata for reference shapes of each obj.References
count_femelements = 0
@@ -537,12 +551,27 @@ def get_pressure_obj_faces_depreciated(femmesh, femobj):
def get_pressure_obj_faces(femmesh, femelement_table, femnodes_ele_table, femobj):
# get the nodes
prs_face_node_set = get_femnodes_by_femobj_with_references(femmesh, femobj) # sorted and duplicates removed
# print('prs_face_node_set: ', prs_face_node_set)
# fill the bit_pattern_dict and search for the faces
bit_pattern_dict = get_bit_pattern_dict(femelement_table, femnodes_ele_table, prs_face_node_set)
pressure_faces = get_ccxelement_faces_from_binary_search(bit_pattern_dict)
if is_solid_femmesh(femmesh):
# get the nodes
prs_face_node_set = get_femnodes_by_femobj_with_references(femmesh, femobj) # sorted and duplicates removed
# print('prs_face_node_set: ', prs_face_node_set)
# fill the bit_pattern_dict and search for the faces
bit_pattern_dict = get_bit_pattern_dict(femelement_table, femnodes_ele_table, prs_face_node_set)
pressure_faces = get_ccxelement_faces_from_binary_search(bit_pattern_dict)
elif is_face_femmesh(femmesh):
pressure_faces = []
# normaly we should call get_femelements_by_references and the group check should be integrated there
if femmesh.GroupCount:
meshfaces = get_femelementface_sets_from_group_data(femmesh, femobj)
# print(meshfaces)
for mf in meshfaces:
# pressure_faces.append([mf, 0])
pressure_faces.append([mf, -1])
# 0 if femmeshface normal == reference face normal direction
# -1 if femmeshface normal opposite reference face normal direction
# easy on plane faces, but on a half sphere ... ?!?
else:
print("Pressure on shell mesh at the moment only supported for meshes with appropriate group data.")
return pressure_faces