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

FEM: Allow shell faces for tie constraint (#18325)

Browse Source 
FEM: Update meshtools.py
This commit is contained in:
FEA-eng
2024-12-09 19:21:09 +01:00
committed by GitHub
parent 74e7df9676
commit 1c35d08d87
1 changed files with 15 additions and 20 deletions
Download Patch File Download Diff File Expand all files Collapse all files

35
src/Mod/Fem/femmesh/meshtools.py
View File

@@ -1585,7 +1585,6 @@ def get_contact_obj_faces(femmesh, femelement_table, femnodes_ele_table, femobj)
# ***** tie faces ****************************************************************************
def get_tie_obj_faces(femmesh, femelement_table, femnodes_ele_table, femobj):
# see comment get_contact_obj_faces
# solid mesh is same as contact, but face mesh is not allowed for tie
# TODO get rid of duplicate code for contact and tie
slave_faces, master_faces = [], []
@@ -1620,27 +1619,23 @@ def get_tie_obj_faces(femmesh, femelement_table, femnodes_ele_table, femobj):
FreeCAD.Console.PrintLog(f"Slave: {slave_ref[0].Name}, {slave_ref}\n")
FreeCAD.Console.PrintLog(f"Master: {master_ref[0].Name}, {master_ref}\n")
if is_solid_femmesh(femmesh):
# get the nodes, sorted and duplicates removed
slaveface_nds = sorted(list(set(get_femnodes_by_refshape(femmesh, slave_ref))))
masterface_nds = sorted(list(set(get_femnodes_by_refshape(femmesh, master_ref))))
# FreeCAD.Console.PrintLog("slaveface_nds: {}\n".format(slaveface_nds))
# FreeCAD.Console.PrintLog("masterface_nds: {}\n".format(slaveface_nds))
# get the nodes, sorted and duplicates removed
slaveface_nds = sorted(list(set(get_femnodes_by_refshape(femmesh, slave_ref))))
masterface_nds = sorted(list(set(get_femnodes_by_refshape(femmesh, master_ref))))
# FreeCAD.Console.PrintLog("slaveface_nds: {}\n".format(slaveface_nds))
# FreeCAD.Console.PrintLog("masterface_nds: {}\n".format(slaveface_nds))
# fill the bit_pattern_dict and search for the faces
slave_bit_pattern_dict = get_bit_pattern_dict(
femelement_table, femnodes_ele_table, slaveface_nds
)
master_bit_pattern_dict = get_bit_pattern_dict(
femelement_table, femnodes_ele_table, masterface_nds
)
# fill the bit_pattern_dict and search for the faces
slave_bit_pattern_dict = get_bit_pattern_dict(
femelement_table, femnodes_ele_table, slaveface_nds
)
master_bit_pattern_dict = get_bit_pattern_dict(
femelement_table, femnodes_ele_table, masterface_nds
)
# get the faces ids
slave_faces = get_ccxelement_faces_from_binary_search(slave_bit_pattern_dict)
master_faces = get_ccxelement_faces_from_binary_search(master_bit_pattern_dict)
elif is_face_femmesh(femmesh):
FreeCAD.Console.PrintError("Shell mesh is not allowed for constraint tie.\n")
# get the faces ids
slave_faces = get_ccxelement_faces_from_binary_search(slave_bit_pattern_dict)
master_faces = get_ccxelement_faces_from_binary_search(master_bit_pattern_dict)
FreeCAD.Console.PrintLog(f"slave_faces: {slave_faces}\n")
FreeCAD.Console.PrintLog(f"master_faces: {master_faces}\n")