Manipulating deformable and fragile objects remains a fundamental challenge in robotics due to complex contact dynamics and strict requirements on object integrity. Existing approaches typically optimize either end-effector design or control strategies in isolation, limiting achievable performance. In this work, we present the first co-design framework that jointly optimizes end-effector morphology and manipulation control for deformable and fragile object manipulation. We introduce (1) a latent diffeomorphic shape parameterization enabling expressive yet tractable end-effector geometry optimization, (2) a stress-aware bi-level co-design pipeline coupling morphology and control optimization, and (3) a privileged-to-pointcloud policy distillation scheme for zero-shot real-world deployment. We evaluate our approach on challenging food manipulation tasks, including grasping and pushing jelly and scooping fillets. Simulation and real-world experiments demonstrate the effectiveness of the proposed method.
Figures & tables
Fig. 1: We jointly optimize end-effector morphology and motion-adaptive control to enable safe, gentle manipulation of deformable and fragile objects. The co-designed end-effector reshapes contact geometry and force distribution, achieving reliable grasping while preserving object integrity, where baseline designs fail and cause breakage.
Fig. 2: Overview of the proposed co-design framework for deformable and fragile object manipulation. (a) A diffeomorphic design space D parameterizes physically valid end-effector geometries. (b) A derivative-free evolutionary optimizer (CMA-ES) searches D to identify an optimal design d∗ that maximizes task performance while minimizing contact-induced damage. (c) For each sampled design di , a design-conditioned controller πdi,prv∗ is optimized using simulator’s privileged signals such as object stress. (d) Rollouts of the privileged policy on the optimized design d are used to collect pointcloud observations and train a student diffusion policy πd,pcd∗ . (e) The distilled policy is deployed zero-shot on the real robot using pointcloud observations, preserving gentle manipulation behaviors learned in simulation.
Fig. 3: Deformation of a 2D gripper contour under a stationary velocity field (SVF). The upper region is frozen (zero velocity), while the lower region deforms smoothly.
Fig. 4: Examples of diffeomorphic deformation parameter fitting. The target shape surface (a) and the fitted mesh (b) are close in terms of Chamfer distance.
cylinder
onigiri
large cube
score ↑
success ↑
stress ↓
score ↑
success ↑
stress ↓
score ↑
success ↑
stress ↓
Ours
143 ± 7
0.97 ± 0.03
9274 ± 124
64 ± 18
0.73 ± 0.08
13973 ± 490
19 ± 48
0.51 ± 0.21
12050 ± 1116
PJ
79 ± 12
0.63 ± 0.05
9366 ± 120
-33 ± 3
0.12 ± 0.04
10224 ± 332
-133 ± 17
0.18 ± 0.04
18893 ± 620
BO
106 ± 28
0.80 ± 0.11
9998 ± 1089
30 ± 66
0.59 ± 0.27
14946 ± 1658
33 ± 53
0.45 ± 0.40
13189 ± 1030
RL
-110 ± 62
0.04 ± 0.02
22015 ± 12105
-244 ± 97
0.09 ± 0.06
38107 ± 19327
-
-
-
TABLE I: Optimal design performance of the jelly grasping task
Fig. 5: End-effector design evolution for grasping a large jelly cube during a typical optimization run. High-scoring designs exhibit larger contact surfaces, which improve grasp stability and reduce contact-induced stress.
cylinder
onigiri
score J
success rate qsucc
stress σmax
score J
success rate qsucc
stress σmax
Ours
85 ± 5
0.760 ± 0.000
1236 ± 108
70 ± 16
0.627 ± 0.076
1177 ± 85
BO
63 ± 16
0.633 ± 0.064
1289 ± 46
65 ± 23
0.627 ± 0.115
1295 ± 64
RL
-39 ± 31
0.280 ± 0.100
2432 ± 522
-81 ± 33
0.287 ± 0.163
5032 ± 659
TABLE II: Optimal design performance of the jelly pushing task
score J↑
success rate qsucc↑
stress σmax↓
Ours
293 ± 97
0.779 ± 0.203
17154 ± 469
BO
282 ± 146
0.683 ± 0.322
13586 ± 2857
TABLE III: Optimal design performance of the fillet scooping task
Fig. 6: Visualization of the design-space landscape projected onto two principal components from the Mlatent -dimensional latent diffeomorphic space. Warmer colors indicate higher end-effector performance for the fillet scooping task. It suggests that thinner and longer end-effector geometries tend to perform better for scooping.
score J
success rate qsucc
stress σmax
Ours
64 ± 18
0.727 ± 0.081
13973 ± 490
Cubic
10 ± 17
0.473 ± 0.083
14894 ± 306
Spherical
13 ± 23
0.540 ± 0.106
15321 ± 386
TABLE IV: Comparison of design spaces
Task
Method
success rate ↑
object damage ↓
Grasp
Ours
9 /9
1 /9
BO
5/9
1 /9
Push
Ours
10 /12
0 /12
BO
9/12
0 /12
TABLE V: Real-World Experiments
Fig. 7: (a) Task-specific end-effector designs optimized, fabricated, and mounted on the robot for real-world evaluation. Two design optimizers Od are compared: CMA-ES (ours) and BO for jelly cylinder grasping and pushing tasks. The co-designed end-effector exhibits curved inner contact surfaces that improve form closure and reduce contact-induced damage to fragile objects. (b) & (c) The roll out of our policy and BO baseline in real world.
Fig. 8: Additional examples of diffeomorphic parameter fitting for jelly grasping and pushing. The target mesh (left, blue) and the fitted mesh (right, purple) closely match in terms of Chamfer distance. The target shapes are drawn from a subset of the N0 designs used to learn the raw diffeomorphic parameters.
Fig. 9: Subset of the diffeomorphic training dataset and corresponding parameter fitting results for the fillet scooping task.
Paper
Design Space Dimension
Action Space Dimension
Co-Design Paradigm
SERL Co-Design [ 8 ]
2
9
Two-stage: RL + GA
Meta-RL Co-Design [ 4 ]
4
16
Two-stage: RL + CMA-ES
Object Adaptation [ 18 ]
4
7
Single-stage: Dual MDP
Object Adaptation [ 18 ]
4
7
Single-stage: Dual MDP
Cross-Embodiment Hand Co-Design [ 14 ]
5
20
Two-stage with Design-Conditioned Policy
CageCoOpt [ 12 ]
6
4
Two-stage: PPO + BO/GA
TABLE VI: Comparison of co-design methods by design and action space dimensionality and paradigm.
Fig. 10: End-effector design evolution for pushing a jelly onigiri during a typical optimization run. High-scoring designs exhibit larger contact surfaces, which improve push stability.
Fig. 11: Example base shape in the LDM training dataset (leftmost) and its augmented variants generated via random rotation, translation, and scaling.
Fig. 12: Average Chamfer distance between original mesh and fitted mesh, with respect to latent design space dimension Mlatent . We choose Mlatent=15 in most of the experiments as a tradeoff.
Fig. 13: Real-world rollouts of the co-designed policies on the three DFOM tasks.
Fig. 14: Graphical user interface for curating end-effector designs used in diffeomorphic design space learning. Shapes from Thingi10K are attached to a standard finger mount at a user-specified pose.
Fig. 15: Parallel, independent simulation environments in Genesis.
Fig. 16: Illustration of synthesized pre-contact configurations pd⋆ and the corresponding motion primitive sequences for jelly grasping, pushing, and fillet scooping. Each row shows the ordered execution stages, including approach, interaction, and post-contact lifting, pushing or insertion.
Fig. 17: illustration of the two-stage evaluation strategy in control optimization (line 6-9 in Algorithm 1 ). pe,i,i∈{1,2,..,I} refers to I candidate pre-contact configurations synthesized from rigid body surrogate. Here I=5 . pe∗ refers to the best pre-contact configuration selected for execution.
Task
λprog
λsucc
λmax2.5
λmax0
λσˉ
Grasping
1×102
2×102
2×10−3
8×10−4
1×10−2
Pushing
3×102
5×102
2×10−3
8×10−4
1×10−2
Scooping
1×102
2×102
2×10−3
8×10−4
1×10−2
TABLE VII: Design objective weights used for different tasks.
Department of Cognitive Robotics, Delft University of Technology, 2628 CN Delft, The Netherlands · Institute of Robotics and Mechatronics, German Aerospace Center (DLR), 82234 Oberpfaffenhofen, Germany