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

Fem: Add support for 2D geometries to HeatFlux constraint

Browse Source
This commit is contained in:
marioalexis
2025-08-29 03:17:16 -03:00
parent 7afa56652c
commit 7027caf0aa
3 changed files with 31 additions and 33 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
src/Mod/Fem/Gui/TaskFemConstraintHeatflux.cpp
View File

@@ -281,10 +281,11 @@ void TaskFemConstraintHeatflux::addToSelection()
if (!subNames.empty()) {
for (const auto& subName : subNames) {
if (subName.substr(0, 4) != "Face") {
QMessageBox::warning(this,
tr("Selection error"),
tr("Selection must only consist of faces!"));
if ((subName.substr(0, 4) != "Face") && (subName.substr(0, 4) != "Edge")) {
QMessageBox::warning(
this,
tr("Selection error"),
tr("Selection must only consist of faces! (edges in 2D models)"));
return;
}
}
@@ -345,10 +346,11 @@ void TaskFemConstraintHeatflux::removeFromSelection()
if (!subNames.empty()) {
for (const auto& subName : subNames) {
if (subName.substr(0, 4) != "Face") {
QMessageBox::warning(this,
tr("Selection error"),
tr("Selection must only consist of faces!"));
if ((subName.substr(0, 4) != "Face") && (subName.substr(0, 4) != "Edge")) {
QMessageBox::warning(
this,
tr("Selection error"),
tr("Selection must only consist of faces! (edges in 2D models)"));
return;
}
}

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

@@ -533,25 +533,24 @@ class MeshSetsGetter:
def get_constraints_heatflux_faces(self):
if not self.member.cons_heatflux:
return
# TODO: use meshtools to get the surfaces (or move to mesh tools)
# see constraint contact or constraint tie and constraint force
# heatflux_obj_face_table: see force_obj_node_load_table
# [
# ("refshape_name:elemname", face_table),
# ...,
# ("refshape_name:elemname", face_table)
# ]
if not self.femnodes_mesh:
self.femnodes_mesh = self.femmesh.Nodes
if not self.femelement_table:
self.femelement_table = meshtools.get_femelement_table(self.femmesh)
if not self.femnodes_ele_table:
self.femnodes_ele_table = meshtools.get_femnodes_ele_table(
self.femnodes_mesh, self.femelement_table
)
for femobj in self.member.cons_heatflux:
# femobj --> dict, FreeCAD document object is femobj["Object"]
heatflux_obj = femobj["Object"]
femobj["HeatFluxFaceTable"] = []
for o, elem_tup in heatflux_obj.References:
for elem in elem_tup:
ho = o.Shape.getElement(elem)
if ho.ShapeType == "Face":
elem_info = f"{o.Name}:{elem}"
face_table = self.mesh_object.FemMesh.getccxVolumesByFace(ho)
femobj["HeatFluxFaceTable"].append((elem_info, face_table))
obj = femobj["Object"]
print_obj_info(obj)
result = []
ref_data = meshtools.pair_obj_reference(obj.References)
for ref_pair in ref_data:
result.append(meshtools.get_ccx_elements(self, ref_pair))
femobj["HeatFluxFaces"] = result
# ********************************************************************************************
# ********************************************************************************************

11
src/Mod/Fem/femsolver/calculix/write_constraint_heatflux.py
View File

@@ -82,10 +82,7 @@ def write_meshdata_constraint(f, femobj, heatflux_obj, ccxwriter):
heatflux_amplitude = ""
f.write(f"*{heatflux_key_word}{heatflux_amplitude}\n")
for ref_shape in femobj["HeatFluxFaceTable"]:
elem_string = ref_shape[0]
face_table = ref_shape[1]
f.write(f"** Heat flux on face {elem_string}\n")
for i in face_table:
# OvG: Only write out the VolumeIDs linked to a particular face
f.write(f"{i[0]},{heatflux_facetype}{i[1]}{heatflux_facesubtype},{heatflux_values}\n")
for feat, surf, is_sub_el in femobj["HeatFluxFaces"]:
f.write("** {0.Name}.{1[0]}\n".format(*feat))
for face, fno in surf:
f.write(f"{face},{heatflux_facetype}{fno}{heatflux_facesubtype},{heatflux_values}\n")