morphopt 4.0.4__tar.gz → 4.0.5__tar.gz

{morphopt-4.0.4 → morphopt-4.0.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: morphopt
-Version: 4.0.4
+Version: 4.0.5
 Summary: A morphing optimization framework using PyTorch and torchfea.
 Requires-Python: >=3.12
 Description-Content-Type: text/markdown
@@ -17,7 +17,7 @@ Requires-Dist: cadquery-ocp>=8.0.1.1.0
 Requires-Dist: tabulate>=0.9.0
 Requires-Dist: imageio>=2.37.3
 Requires-Dist: bspmap>=1.0.3
-Requires-Dist: torchfea>=1.1.2
+Requires-Dist: torchfea>=1.1.3
 Requires-Dist: trimesh>=4.11.5
 Requires-Dist: networkx>=3.6.1
 Requires-Dist: scipy>=1.17.1

{morphopt-4.0.4 → morphopt-4.0.5}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "morphopt"
-version = "4.0.4"
+version = "4.0.5"
 description = "A morphing optimization framework using PyTorch and torchfea."
 readme = "README.md"
 requires-python = ">=3.12"
@@ -18,7 +18,7 @@ dependencies = [
     "tabulate>=0.9.0",
     "imageio>=2.37.3",
     "bspmap>=1.0.3",
-    "torchfea>=1.1.2",
+    "torchfea>=1.1.3",
     "trimesh>=4.11.5",
     "networkx>=3.6.1",
     "scipy>=1.17.1",

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/constants.py

@@ -181,13 +181,6 @@ class _Constants:
                     class TOOL_ID:
                         DEFAULT = ""
 
-                class FRAGMENT:
-                    class OBJECT_ID:
-                        DEFAULT = ""
-
-                    class TOOL_ID:
-                        DEFAULT = ""
-
                 class TRANSLATE:
                     class INPUT_ID:
                         DEFAULT = ""

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/optcore/modelparams/parts/box_part.py

@@ -80,10 +80,10 @@ class BoxPart(BasePart):
                 elems_index=torch.arange(len(cells), device=part.nodes.device),
                 elems=torch.as_tensor(cells, device=part.nodes.device),
             ),
-            name="design",
+            name="C3D8",
         )
         if self.mesh_order == 2:
-            part.convert_linear_to_quadratic_elements(["design"], ["design"])
+            part.convert_linear_to_quadratic_elements(["C3D8"], ["C3D8"])
 
         node_coordinates = part.nodes.detach().cpu().numpy()
         bounds = np.array(

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/application/editing.py

@@ -152,7 +152,20 @@ def set_node_field(node: Node, key: str, value, problem: ProblemDefinition | Non
         if isinstance(node, FEAComponentNode) and key in {*INSTANCE_REFERENCE_FIELDS, *REFERENCE_POINT_FIELDS}:
             return problem.set_component_reference(node, key, value)
         if key == "part_name" and isinstance(node, MaterialNode):
-            return problem.set_material_part(node, value)
+            previous_part = node.part
+            if not problem.set_material_part(node, value):
+                return False
+            if node.is_updatable and node.part is not previous_part:
+                from ..visualization.geometry import PreviewGeometry
+                from ..visualization.scene import MessageItem
+
+                scene = PreviewGeometry().build(problem, "part", node.part)
+                bounds = (
+                    None if any(isinstance(item, MessageItem) for item in scene.items)
+                    else scene.bounds
+                )
+                node.bounding_box = list(bounds) if bounds is not None else [0.0] * 6
+            return True
     node.set_field(key, value)
     return True
 
@@ -266,13 +279,26 @@ class CustomClassSession:
         self._replace_uses(old, name)
         return True
 
+    def dependent_classes(self, node: CustomClassNode) -> list[CustomClassNode]:
+        """Classes whose inheritance chain includes the selected custom class."""
+        names = {node.name}
+        while True:
+            children = {custom.name for custom in self.problem.custom_classes
+                        if custom.base_class in names}
+            if children <= names:
+                break
+            names.update(children)
+        return [custom for custom in self.problem.custom_classes
+                if custom is not node and custom.name in names]
+
     def remove(
         self, node: CustomClassNode | MethodOverrideNode | CustomFunctionNode
     ) -> CustomClassNode | None:
         from ..model.problem import CustomClassNode
         if isinstance(node, CustomClassNode):
-            self.problem.custom_classes.remove(node)
-            self._replace_uses(node.name, "")
+            for custom in [node, *self.dependent_classes(node)]:
+                self.problem.custom_classes.remove(custom)
+                self._replace_uses(custom.name, "")
             return None
         owner = next(custom for custom in self.problem.custom_classes if node in custom.members)
         owner.remove_child(node)

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/application/fields.py

@@ -119,11 +119,15 @@ class SnippetInputs:
         self.problem = problem
 
     def values(self, source: str) -> list:
+        catalog = PartSetCatalog(self.problem)
         providers = {
             "load_steps": lambda: list(self.problem.steps.step_values) if self.problem.steps else [],
             "instances": self.problem.instance_names,
             "reference_points": self.problem.reference_point_names,
-            "elements": self.problem.element_names,
+            "elements": lambda: list(dict.fromkeys(
+                name for part in self.problem.parts()
+                for name in catalog.for_part(part).elements
+            )),
             "force_loads": lambda: [
                 node.name
                 for node in self.problem.fea_components()

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/application/part_sets.py

@@ -91,8 +91,11 @@ class PartSetCatalog:
                              *metadata.get("node_sets", {}))
             )
         return PartSets(
-            nodes=tuple(names) + end_nodes + tuple(f"{name}_nodes" for name in faces) + tuple(node_sets),
-            surfaces=tuple(names) + tuple(faces) + tuple(surface_sets),
+            nodes=("all", "surface_all") + tuple(names) + end_nodes + tuple(f"{name}_nodes" for name in faces) + tuple(node_sets),
+            surfaces=tuple(dict.fromkeys([
+                "extern", part.exterior_surface or "extern", *names, *faces, *surface_sets,
+            ])),
+            elements=("C3D4", "C3D6") if isinstance(part, OffsetBoundaryPartNode) else ("C3D4",),
         )
 
     @staticmethod
@@ -100,7 +103,7 @@ class PartSetCatalog:
         return PartSets(
             nodes=("all",) + tuple(f"nodes_{axis}{side}" for axis in "xyz" for side in ("min", "max")),
             surfaces=(part.exterior_surface or "extern",) + tuple(f"surface_{axis}{side}" for axis in "xyz" for side in ("min", "max")),
-            elements=("design",),
+            elements=("C3D8",),
         )
 
     @staticmethod
@@ -109,6 +112,8 @@ class PartSetCatalog:
         return PartSets(
             nodes=("all", "surface_all") + tuple(f"{name}_nodes" for name in names),
             surfaces=tuple(dict.fromkeys([part.exterior_surface or "extern", *names])),
+            # Quadratic conversion preserves the element block's name.
+            elements=("C3D4",),
         )
 
     @staticmethod

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/application/quickstart.py

@@ -578,7 +578,7 @@ def _topology(options, geometry, loads, materials):
     loads.add_component(FEAComponentNode.create("ConcentratedForce", name="tip_force", rp_name="RP_tip"))
     materials.add_material(MaterialNode.create(
         "SIMP_BSPFieldMaterials", name="density", part_name="body",
-        elementname="design", material_model="NeoHookeanLnJ", mu=options.mu,
+        elementname="C3D8", material_model="NeoHookeanLnJ", mu=options.mu,
         kappa=options.kappa, mumax=options.mu, kappamax=options.kappa,
         bounding_box=[0., options.total_length, 0., options.height, 0., options.depth],
         degree=CONSTANT.UI.QUICKSTART.GENERATION.SIMP_DEGREE,
@@ -604,7 +604,7 @@ def _topology(options, geometry, loads, materials):
                 {"type": "VolFrac", "params": {
                     "volfrac_min": options.volume_fraction, "volfrac_max": options.volume_fraction,
                     "penalty": CONSTANT.UI.QUICKSTART.GENERATION.SIMP_VOLUME_FRACTION_PENALTY,
-                    "elementname": "design",
+                    "elementname": "C3D8",
                 }},
             ],
         }],
@@ -693,8 +693,6 @@ class ObjectiveDraft:
         return len(self.objectives) - 1
 
     def remove(self, index: int) -> None:
-        if len(self.objectives) <= CONSTANT.UI.QUICKSTART.LIMITS.OBJECTIVE_MIN_COUNT:
-            raise ValueError(localized_message("ui.errors.quickstart.keep_at_least_one_objective_term"))
         self.objectives.pop(index)
 
     def set_method(self, method: str) -> None:

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/application/updaters.py

@@ -135,6 +135,18 @@ class UpdaterItems:
         ) if custom.name not in taken)
         return choices
 
+    def default_element_name(self) -> str:
+        """Use the material's block, or the sole block of its owning Part."""
+        from .part_sets import PartSetCatalog
+
+        material = self.problem.config_material_node(self.config)
+        if material is None:
+            return ""
+        if material.elementname:
+            return material.elementname
+        elements = PartSetCatalog(self.problem).for_part(material.part).elements
+        return elements[0] if len(elements) == 1 else ""
+
     def add(self, item_type: str, custom_class: str = "") -> bool:
         item = symmetry_item(item_type, self.surfaces) if self.category == "geometry_symmetry" else updater_item(
             item_type,
@@ -143,6 +155,8 @@ class UpdaterItems:
         )
         if item is None:
             return False
+        if self.group == "materials" and item_type == "VolFrac":
+            item["params"]["elementname"] = self.default_element_name()
         self.config.setdefault(self.key, []).append(item)
         return True
 

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/codegen/generator.py

@@ -562,8 +562,8 @@ def generate_base_source(
 
     Code-launched runs never open an observer window; results are watched by
     the MorphOpt UI's in-window observer (disk polling per iteration).
-    ``for_class_catalog`` emits inert objective methods so editor class discovery
-    also works while the user is still defining objective and metric terms.
+    ``for_class_catalog`` emits inert objective methods and updater classes so
+    editor class discovery also works while terms are still being defined.
     """
     if asset_paths is None:
         asset_paths, _ = collect_python_assets(problem)
@@ -672,7 +672,7 @@ def generate_base_source(
     a("")
 
     # -------- Updater -----------------------------------------------------
-    _emit_updater(a, problem, template, upd_device)
+    _emit_updater(a, problem, template, upd_device, for_class_catalog=for_class_catalog)
 
     # -------- main --------------------------------------------------------
     a("")
@@ -800,7 +800,7 @@ def _emit_part_classes(
         a("                def define_instance(self) -> None:")
         instances = part.instances
         if not instances:
-            a("                    super().define_instance()")
+            a("                    pass")
         for instance in instances:
             a("")
             pose = (
@@ -1155,6 +1155,8 @@ def _emit_updater(
     problem: ProblemDefinition,
     template: SchemeTemplate,
     device: str,
+    *,
+    for_class_catalog: bool = False,
 ) -> None:
     from .custom_classes import updater_class_name
 
@@ -1203,10 +1205,10 @@ def _emit_updater(
             if target.type == "OffsetBoundaryPart"
             else GEOMETRY_UPDATER_BASE
         )
-        _emit_nested_updater(a, class_name, base, config, problem)
+        _emit_nested_updater(a, class_name, base, config, problem, for_class_catalog=for_class_catalog)
         a("")
     for class_name, config in zip(material_classes, mat_cfgs):
-        _emit_nested_updater(a, class_name, MATERIAL_UPDATER_BASE, config, problem)
+        _emit_nested_updater(a, class_name, MATERIAL_UPDATER_BASE, config, problem, for_class_catalog=for_class_catalog)
         a("")
 
 
@@ -1272,6 +1274,8 @@ def _emit_nested_updater(
     base: str,
     cfg,
     problem: ProblemDefinition,
+    *,
+    for_class_catalog: bool = False,
 ) -> None:
     """Emit one sub-updater class.
 
@@ -1299,6 +1303,12 @@ def _emit_nested_updater(
         f"max_step_length={_render_value(cfg.get('max_step_length', defaults.MAX_STEP_LENGTH), problem)})"
     )
 
+    if for_class_catalog:
+        a("")
+        a("            def define_objective(self) -> None:")
+        a("                pass")
+        return
+
     body: list[str] = []
     # NOTE: the model stores the lists under ``objective_functions`` / ``constraints``
     # while the catalogue (UPDATER_CATALOG) is keyed ``objectives`` / ``constraints``.

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/locales/en-US.json

@@ -1,4 +1,6 @@
 {
+  "ui.component.delete_dependencies_title": "Delete bound objects",
+  "ui.component.delete_dependencies_message": "Deleting “{component}” will also delete these bound objects:\n{dependencies}\nContinue?",
   "ui.schema.fields.control_point_spacing.label": "Target control point spacing",
   "ui.schema.fields.design_faces.label": "Design faces",
   "ui.assets.brep_metadata": "BREP topology and sets",
@@ -762,7 +764,6 @@
   "ui.schema.catalog.cad_operation.fuse.label": "Fuse",
   "ui.schema.catalog.cad_operation.cut.label": "Cut",
   "ui.schema.catalog.cad_operation.intersect.label": "Intersect",
-  "ui.schema.catalog.cad_operation.fragment.label": "Fragment",
   "ui.schema.catalog.cad_operation.translate.label": "Translate",
   "ui.schema.catalog.cad_operation.rotate.label": "Rotate",
   "ui.schema.catalog.cad_operation.dilate.label": "Scale",

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/locales/zh-CN.json

@@ -1,4 +1,6 @@
 {
+  "ui.component.delete_dependencies_title": "删除绑定对象",
+  "ui.component.delete_dependencies_message": "删除“{component}”将同时删除以下绑定对象：\n{dependencies}\n是否继续？",
   "ui.schema.fields.control_point_spacing.label": "控制点目标间距",
   "ui.schema.fields.design_faces.label": "参与更新的面",
   "ui.assets.brep_metadata": "BREP 拓扑与集合",
@@ -762,7 +764,6 @@
   "ui.schema.catalog.cad_operation.fuse.label": "融合",
   "ui.schema.catalog.cad_operation.cut.label": "切除",
   "ui.schema.catalog.cad_operation.intersect.label": "相交",
-  "ui.schema.catalog.cad_operation.fragment.label": "碎片分割",
   "ui.schema.catalog.cad_operation.translate.label": "平移",
   "ui.schema.catalog.cad_operation.rotate.label": "旋转",
   "ui.schema.catalog.cad_operation.dilate.label": "缩放",

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/model/base_nodes.py

@@ -427,7 +427,7 @@ class PartNode(TypedNode):
         super().__init__(
             kind=kind, name=name, params=data, children=children
         )
-        if not self.instances:
+        if not self.instances and children is None:
             from .node_types.structural import InstanceNode
 
             self.add_instance(InstanceNode(name=f"{self.assembly_name or 'part'}-1"))

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/model/configuration.py

@@ -84,21 +84,29 @@ def materials_updater_config(
     max_step_iter: int = 50,
     max_step_length: float = CONSTANT.KERNEL.SIMP_UPDATE.MAX_STEP_LENGTH,
     if_update: bool | list[bool] = True,
-    objective_functions: tuple = (),
-    constraints: tuple = (),
+    objective_functions: tuple | None = None,
+    constraints: tuple | None = None,
 ) -> dict:
     """Config of one material sub-optimizer.
 
     ``material_name`` is the material definition it updates; an empty name
     means "the only design-carrying material in the model".
+    New configs include sensitivity, variable bounds and volume fraction.
+    Explicitly empty collections remain empty.
     """
     return {
         "material_name": str(material_name or ""),
         "max_step_iter": int(max_step_iter),
         "max_step_length": float(max_step_length),
         "if_update": if_update,
-        "objective_functions": list(objective_functions),
-        "constraints": list(constraints),
+        "objective_functions": list(objective_functions) if objective_functions is not None else [
+            updater_objective("Sensitivity"),
+        ],
+        "constraints": list(constraints) if constraints is not None else [
+            updater_constraint("MinValue"),
+            updater_constraint("MaxValue"),
+            updater_constraint("VolFrac"),
+        ],
     }
 
 

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/model/fields.py

@@ -29,6 +29,7 @@ class FieldSpec:
     ints: bool = False
     read_only: bool = False
     hidden: bool = False
+    editable: bool = True
 
     def new_default(self) -> Any:
         return _mutable_value(self.default)

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/model/node_types/__init__.py

@@ -6,7 +6,6 @@ from .cad import (
     CutOperationNode,
     CylinderOperationNode,
     DilateOperationNode,
-    FragmentOperationNode,
     FuseOperationNode,
     IntersectOperationNode,
     MirrorOperationNode,
@@ -127,7 +126,6 @@ __all__ = (
     "FairnessPenaltyNode",
     "MeshSurfaceNode",
     "BREPSurfaceNode",
-    "FragmentOperationNode",
     "FuseOperationNode",
     "GeometryNode",
     "HelperFunctionNode",

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/model/node_types/cad.py

@@ -418,33 +418,6 @@ class IntersectOperationNode(CADOperationNode):
     )
 
 
-class FragmentOperationNode(CADOperationNode):
-    """UI definition of Fragment."""
-
-    type_name: ClassVar[str] = "Fragment"
-    backend_class: ClassVar[str] = "Fragment"
-    label_key: ClassVar[str] = "ui.schema.catalog.cad_operation.fragment.label"
-    group: ClassVar[str] = "boolean"
-    operation_kind: ClassVar[str] = "fragment"
-    Defaults = CONSTANT.UI.SCHEMA_FIELDS.CAD_OPERATIONS.FRAGMENT
-    object_id: str
-    tool_id: str
-    fields: ClassVar[tuple[FieldSpec, ...]] = (
-        FieldSpec(
-            "object_id",
-            "ui.schema.fields.object_id.label.object_id",
-            "combo",
-            Defaults.OBJECT_ID.DEFAULT,
-        ),
-        FieldSpec(
-            "tool_id",
-            "ui.schema.fields.tool_id.label.tool_id",
-            "combo",
-            Defaults.TOOL_ID.DEFAULT,
-        ),
-    )
-
-
 class TranslateOperationNode(CADOperationNode):
     """UI definition of Translate."""
 

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/model/node_types/parts.py

@@ -87,6 +87,7 @@ class BoundaryPartNode(PartNode):
             Defaults.MESH_ORDER.DEFAULT,
             doc_key="ui.schema.fields.mesh_order.help.element_order_of_the_generated_mesh",
             choices=Defaults.MESH_ORDER.CHOICES,
+            editable=False,
         ),
         FieldSpec(
             "exterior_surface",
@@ -210,6 +211,7 @@ class BoxPartNode(PartNode):
             Defaults.MESH_ORDER.DEFAULT,
             doc_key="ui.schema.fields.mesh_order.help.text",
             choices=Defaults.MESH_ORDER.CHOICES,
+            editable=False,
         ),
         FieldSpec(
             "exterior_surface",
@@ -308,6 +310,7 @@ class CADPartNode(PartNode):
             Defaults.MESH_ORDER.DEFAULT,
             doc_key="ui.schema.fields.mesh_order.help.text",
             choices=Defaults.MESH_ORDER.CHOICES,
+            editable=False,
         ),
         FieldSpec(
             "exterior_surface",

{morphopt-4.0.4 → morphopt-4.0.5}/src/morphopt/ui/model/problem.py

@@ -8,6 +8,7 @@ from __future__ import annotations
 
 from collections.abc import Iterator
 from copy import deepcopy
+import re
 
 from ...constants import CONSTANT
 from ..i18n import localized_message
@@ -715,6 +716,26 @@ class ProblemDefinition:
                 )
             )
         position = owner.surfaces.index(surface)
+        for component in self.fea_components():
+            for field, attr in INSTANCE_REFERENCE_FIELDS.items():
+                if getattr(component, attr) not in owner.instances:
+                    continue
+                suffix = field.removeprefix("instance_name")
+                for key in ("surface_name" + suffix, "set_nodes_name" + suffix):
+                    if not component.has_field(key):
+                        continue
+                    name = component.get_field(key) or ""
+                    if re.fullmatch(r"face_\d+(_nodes)?", name):
+                        component.set_field(key, "")
+                        continue
+                    match = re.fullmatch(r"surface_(\d+)_(.+)", name)
+                    if match is None:
+                        continue
+                    index = int(match[1])
+                    if index == position:
+                        component.set_field(key, "")
+                    elif index > position:
+                        component.set_field(key, f"surface_{index - 1}_{match[2]}")
         self._remove_surface_state(owner, position)
         owner.remove_child(surface)
         self.align_surface_dependent_state()
@@ -799,7 +820,7 @@ class ProblemDefinition:
         return added
 
     def remove_instance(self, instance: InstanceNode) -> None:
-        """Remove one Instance, keeping a Part's first identity Instance."""
+        """Remove an Instance and clear its optional references."""
         owner = self.owner_of_instance(instance)
         if owner is None:
             raise ValueError(
@@ -807,17 +828,20 @@ class ProblemDefinition:
                     "ui.errors.problem.the_instance_is_not_part_of_this_problem"
                 )
             )
-        if len(owner.instances) <= 1:
-            raise ValueError(
-                localized_message(
-                    "ui.errors.problem.a_part_must_keep_at_least_one_instance"
-                )
-            )
         owner.remove_child(instance)
         for component in self.fea_components():
+            if (component.has_field("obj_name") and component.obj_name == instance.name
+                    and component.obj_type in ("auto", "instance")):
+                component.obj_name = ""
             for field, attr in INSTANCE_REFERENCE_FIELDS.items():
                 if getattr(component, attr) is instance:
                     component.set_reference(field, None)
+                    suffix = field.removeprefix("instance_name")
+                    for key in ("surface_name" + suffix, "set_nodes_name" + suffix,
+                                "element_name" + suffix):
+                        if component.has_field(key):
+                            component.set_field(key, "")
+        self._clear_objective_reference("instance_name", instance.name)
 
     def move_instance(self, instance: InstanceNode, offset: int) -> bool:
         """Move an Instance in its Part's declaration order."""
@@ -850,8 +874,8 @@ class ProblemDefinition:
             )
         return part.add_cad_operation(operation)
 
-    def remove_cad_operation(self, operation: CADOperationNode) -> None:
-        """Remove one CAD feature from its owning Part."""
+    def cad_operation_dependents(self, operation: CADOperationNode) -> list[CADOperationNode]:
+        """Return all features transitively dependent on this operation."""
         owner = self.owner_of_cad_operation(operation)
         if owner is None:
             raise ValueError(
@@ -859,22 +883,34 @@ class ProblemDefinition:
                     "ui.errors.problem.the_cad_operation_is_not_part_of_this_problem"
                 )
             )
-        dependents = [
-            item
-            for item in owner.operations
-            if item is not operation
-            and operation.operation_id in (item.object_id, item.tool_id, item.input_id)
-        ]
-        if dependents:
-            names = ", ".join(item.operation_id for item in dependents)
-            raise ValueError(
-                localized_message(
-                    "ui.errors.problem.cannot_remove_value_it_is_used_by_value",
-                    value_0=operation.operation_id,
-                    value_1=names,
-                )
-            )
-        owner.remove_child(operation)
+        removed = {operation.operation_id}
+        while True:
+            dependent_ids = {
+                item.operation_id for item in owner.operations
+                if removed.intersection((item.object_id, item.tool_id, item.input_id))
+            }
+            if dependent_ids <= removed:
+                break
+            removed.update(dependent_ids)
+        return [item for item in owner.operations
+                if item is not operation and item.operation_id in removed]
+
+    def remove_cad_operation(self, operation: CADOperationNode) -> None:
+        """Remove a CAD feature and its dependent features after UI confirmation."""
+        dependents = self.cad_operation_dependents(operation)
+        owner = self.owner_of_cad_operation(operation)
+        for item in [operation, *dependents]:
+            owner.remove_child(item)
+        for component in self.fea_components():
+            for field, attr in INSTANCE_REFERENCE_FIELDS.items():
+                if getattr(component, attr) not in owner.instances:
+                    continue
+                suffix = field.removeprefix("instance_name")
+                for key in ("surface_name" + suffix, "set_nodes_name" + suffix):
+                    if component.has_field(key) and re.fullmatch(
+                        r"face_\d+(_nodes)?", component.get_field(key) or ""
+                    ):
+                        component.set_field(key, "")
 
     def move_cad_operation(self, operation: CADOperationNode, offset: int) -> bool:
         """Move one CAD feature without disturbing other Part children."""
@@ -1007,17 +1043,21 @@ class ProblemDefinition:
             if material.part is part:
                 material.part = None
                 material._part_name = ""
+                material.elementname = ""
+                if self.updater is not None:
+                    for config in self.updater.materials:
+                        if config.get("target") is material:
+                            for item in config.get("constraints", []):
+                                if "elementname" in item.get("params", {}):
+                                    item["params"]["elementname"] = ""
         if self.updater is not None:
             self.updater.geometry[:] = [
                 config
                 for config in self.updater.geometry
                 if config.get("target") is not part
             ]
-        for instance in part.instances:
-            for load in self.fea_components():
-                for field, attr in INSTANCE_REFERENCE_FIELDS.items():
-                    if getattr(load, attr) is instance:
-                        load.set_reference(field, None)
+        for instance in list(part.instances):
+            self.remove_instance(instance)
         geometry.remove_child(part)
         self.align_surface_dependent_state()
         self._refresh_surface_flags()
@@ -1088,12 +1128,6 @@ class ProblemDefinition:
                     "ui.errors.problem.material_is_not_part_of_this_problem"
                 )
             )
-        if len(section.assignments) <= 1:
-            raise ValueError(
-                localized_message(
-                    "ui.errors.problem.a_problem_must_keep_at_least_one_material"
-                )
-            )
         section.remove_child(material)
         if self.updater is not None:
             self.updater.materials[:] = [
@@ -1155,9 +1189,33 @@ class ProblemDefinition:
                 item for item in self.objective.jacobian_needed if item is not component
             ]
         for load in self.fea_components():
+            if load.has_field("obj_name") and load.obj_name == component.name:
+                load.obj_name = ""
             for field, attr in REFERENCE_POINT_FIELDS.items():
                 if getattr(load, attr) is component:
                     load.set_reference(field, None)
+        if component.type == "ReferencePoint":
+            self._clear_objective_reference("reference_point_name", component.name)
+        else:
+            for key in ("force_load_name", "moment_load_name"):
+                self._clear_objective_reference(key, component.name)
+            if self.objective is not None:
+                for term in self.objective.objective_terms:
+                    if term.get("type") != "workspace":
+                        continue
+                    for entry in term.get("inputs", []):
+                        entry["parameters"] = [pair for pair in entry["parameters"]
+                                               if pair[0] != component.name]
+                    term["inputs"] = [entry for entry in term.get("inputs", [])
+                                      if entry["parameters"]]
+
+    def _clear_objective_reference(self, key: str, name: str) -> None:
+        """Clear explicit objective and metric references to a deleted object."""
+        if self.objective is None:
+            return
+        for term in [*self.objective.objective_terms, *self.objective.metrics_terms]:
+            if term.get(key) == name:
+                term[key] = ""
 
     def rename_definition(
         self,