CQDC-Warp: GPU-Accelerated Differentiable Quasi-Dynamic Simulation for Contact-Rich Manipulation
summary
The gist
The gist: CQDC-Warp presents a GPU-accelerated differentiable quasi-dynamic simulator for contact-rich manipulation that supports batched Model Predictive Control (MPC) by allowing intended contacts
In short
CQDC-Warp is a GPU-accelerated simulator for contact-rich manipulation that uses differentiable dynamics to support batched Model Predictive Control (MPC). It allows planners to specify intended robot contacts beforehand, which then guide the optimization process. The simulation determines the actual contact geometry and friction, enabling efficient trajectory generation for complex tasks.
Key concepts
- CQDC-Warp
- A GPU-accelerated simulator designed for batched Model Predictive Control (MPC). It enables contact-guided trajectory optimization by letting intended contacts guide planning while simulated dynamics model the resulting contact geometry and friction.
- Compliant-Contact Formulation
- This improvement replaces hard contacts with a formulation that penalizes contact deformation while keeping friction constraints. It minimizes an objective function including a term for contact softness, allowing the simulation to handle softer physical interactions.
- Contact-Solve Implicit Gradients
- A method used to compute gradients necessary for optimization. It derives the displacement Jacobian by applying the implicit function theorem to the stationarity equation of the contact dynamics, providing gradients even when contacts are solved implicitly.
Terminology used across episodes
This episode discusses
- CQDC-Warp: GPU-Accelerated Differentiable Quasi-Dynamic Simulation for Contact-Rich Manipulation · Paper Radio
- Dojo: A Differentiable Physics Engine for Robotics
- Highly-Efficient Differentiable Simulation for Robotics · Paper Radio
- Solving Rubik's Cube with a Robot Hand
The paper
CQDC-Warp: GPU-Accelerated Differentiable Quasi-Dynamic Simulation for Contact-Rich Manipulation · Read on arXiv
Yanchen Ding, Jiayi Chen, Jianfeng Xu, Wei Zhao, Li Yi, He Wang
Peking University · Galbot Institute of Technology and Research Laboratory of the Chinese Academy of Sciences (implied by context, though only Peking University and Tsinghua University are explicitly listed as affiliations)
Efficient synthesis of contact-rich manipulation trajectories is essential for downstream robot learning. Gradient-based planning offers a promising approach, but jointly optimizing motion and contact decisions remains challenging. Meanwhile, many existing planners run on CPUs with limited throughput. We present CQDC-Warp, a GPU-accelerated differentiable quasi-dynamic simulator for batched model predictive control (MPC). We use it to study contact-guided trajectory optimization, in which intended robot--object contacts are specified before planning and encouraged through soft costs, while the simulated dynamics determine contact geometry and frictional response without a prescribed contact-mode schedule. Across 17, 408 trials on the 17-object Allegro in-hand rotation benchmark, our framework achieves 79.8% success, compared with 73.7% for CFree and 65.2% for IMPACT, and reduces mean trajectory-generation time by factors of 27.0 and 50.2, respectively. On a single NVIDIA RTX 4090 with 2, 048 parallel worlds, forward and backward throughput reach 2.06 times10 5 and 2.25 times10 5 world-steps per second, respectively. Further experiments cover manipulation with LEAP and Shadow hands, environmental contact exploitation, bimanual manipulation, whole-body control, and real-world trajectory replay on a Sharpa hand.
Transcript
Introduction to the show: ident: Robotics Radio. Generated commentary on the latest robotics and control papers.
Rosa: Today's paper: "CQDC-Warp: GPU-Accelerated Differentiable Quasi-Dynamic Simulation for Contact-Rich Manipulation".
Dev: The gist:
Rosa: First, who's behind it and why it matters.
Title and authors: Rosa: So we're looking at this paper, "CQDC-Warp: GPU-Accelerated Differentiable Quasi-Dynamic Simulation for Contact-Rich Manipulation," and the main idea is that they built a simulator that lets you do batched Model Predictive Control, or MPC, which is a way to plan robot movements while thinking about contacts simultaneously.
Dev: That's right. It’s essentially a GPU-accelerated simulator designed so that you can tell the planner what contacts *should* happen before planning starts and let the simulation figure out exactly what the contact geometry and friction response will actually be.
Taro: What I find interesting is that it lets you guide trajectory optimization through those intended contacts specified via soft costs, rather than having to fix a specific contact mode schedule beforehand.
Rosa: Exactly. They extend the Convex QuasiDynamic Contact model by introducing this compliant-contact formulation which regularizes the contact solve and revises the inertial model to keep coupling between actuated and unactuated coordinates so you don't get those issues where joint motion induces base motion.
Dev: That revised metric allows them to do an analytical reduction, expressing the optimal control input as a function of delta q through cone projection, which they then substitute into the formulation. It also gives them implicit derivatives of the contact dynamics and sampled surrogate derivatives of contact geometry for gradient computation.
Taro: When I look at those gradients, it sounds like they use witness-pair surrogates to estimate things like the contact midpoint c, unit normal n, and signed gap g together using shared distance-based softmax weights. That seems like a neat way to handle the geometry part.
Rosa: It does seem clever because it's how they compute those geometry gradients, which is one of the two main ways they generate gradients for the planning process. They also have contact-solve implicit gradients derived by applying the implicit function theorem to that stationarity equation (five) <ref:2610.21349#pg2>.
Dev: And when you put this into practice, they implement CQDC-Warp using NVIDIA Warp for GPU-parallel computation and automatic differentiation, which lets them augment components like forward kinematics and collision detection so they can generate adjoints.
Title and authors: Taro: So the system is set up to handle these complex contact decisions jointly with the motion optimization in a batched MPC framework where the intended contacts are specified before planning and encouraged through soft costs in the objective function.
Rosa: That's the core of it. The planner is trying to minimize that cost function, which includes positional costs, regularization for velocity and acceleration, and these contact costs that measure how far the signed gap gk from its target g star deviates.
Dev: And what I like about this approach is that it guides contact realization without having to prescribe the witness point or the frictional mode directly in the dynamics. It lets the simulated dynamics handle those details based on what's happening during planning.
Taro: But Rosa, when we look at how they've improved upon previous work, they introduced two specific modeling enhancements over CQDC. The first is this compliant-contact formulation which replaces the hard-contact model with something that penalizes contact deformation while keeping the friction cone constraint intact and setting a contact softness parameter rho.
Rosa: And the second improvement is that revised inertial model we talked about, which keeps coupling between actuated and unactuated coordinates by including inertia at both coordinates through the full generalized mass matrix, fixing those discrepancies where joint motion causes base motion in older CQDC models.
Dev: It sounds like they are really pushing for a more robust simulation environment that handles the physics of contact better while keeping the optimization loop efficient. This is a big deal for real-time control systems.
Taro: And on the experimental side, across seventeen thousand four hundred eight trials on the seventeen-object Allegro in-hand rotation benchmark, this framework achieved a seventy-nine point eight percent success rate in in-hand rotation compared to seventy-three point seven percent for CFree and sixty-five point two percent for IMPACT <ref:2610.21349#pg1,across 17,408 trials on the 17-object Allegro in-hand rotation>.
Title and authors: Rosa: And it's not just about the success rate; they also showed that the framework significantly reduces mean trajectory-generation time, achieving speedups of factors of twenty-seven point zero and fifty point two on CFree and IMPACT, respectively, which makes the planning loop much faster than existing methods.
Dev: That throughput number is pretty impressive for a GPU simulation, reaching two point zero six times ten to the fifth world-steps per second on a single NVIDIA RTX four thousand ninety with two thousand forty-eight parallel worlds.
Taro: It also supports real-world trajectory replay on a Sharpa hand, which shows they're thinking about how this can move toward sim-to-real deployment in practice.
Rosa: So to wrap up on CQDC-Warp, the main implication is that you can now efficiently synthesize contact-rich manipulation trajectories using GPU batched MPC where intended contacts guide planning while the simulated dynamics determine the actual contact geometry and frictional response.
Dev: It means we can study complex contact decisions jointly with motion optimization in a way that was challenging before because of CPU limitations, and we get these huge speedups.
Taro: For me, it's about making sure the system actually recovers gracefully when the world misbehaves during these contact-rich tasks.
Rosa: I think the main thing to watch is that the planner still relies on suitable initializations and intended contacts specified before planning, which means we still need further evaluation on long-horizon tasks with frequent contact transitions.
Dev: Also, they noted that objective weights still require tuning across different tasks and robot embodiments even though parameters are shared within each task family.
Taro: And I'd add that the paper points out that the use of CQDC-Warp for large-scale manipulation data generation and its impact on downstream policy learning still needs to be evaluated.
Rosa: So that's where we leave it for this paper, "CQDC-Warp: GPU-Accelerated Differentiable Quasi-Dynamic Simulation for Contact-Rich Manipulation." It’s a solid step in making planning faster and more flexible for contact manipulation tasks.
Dev: It is certainly a strong tool for generating data and testing control policies under these kinds of constraints.
Taro: I'm looking forward to seeing how this simulation helps us test those failure modes in the real world.
The paper's summary: Rosa: So, CQDC-Warp is this GPU simulator that lets you plan robot movements while also planning what contacts should happen, using those intended contacts to guide the whole thing through soft costs in your objective function.
Dev: Right, it’s like you tell the planner "Link A needs to touch Object B," and then the simulation handles figuring out exactly where they touch and how much friction there is based on its own dynamics.
Taro: What I find really cool is that they didn't just use a hard contact model. They introduced this compliant-contact formulation which penalizes contact deformation while still keeping the friction cone constraint active, using a parameter rho to set how soft the contact is.
Rosa: And they also fixed a big problem with the inertial model; they made sure there’s coupling between actuated and unactuated coordinates by including inertia at both places through that full generalized mass matrix.
Dev: That coupling fix is important because in older models, you could get this weird effect where just moving a joint causes the base of the robot to move due to inertia. This revised metric lets them do an analytical reduction, expressing the optimal control input as a function of delta q through cone projection.
Taro: And they give you implicit derivatives for those contact dynamics and sampled surrogate derivatives for contact geometry so you can compute gradients for the optimization loop.
Rosa: They use witness-pair surrogates to estimate things like the contact midpoint c, unit normal n, and signed gap g together using some shared distance-based softmax weights for the geometry gradients.
Dev: That’s how they get those geometry derivatives, which is one of the two main ways they compute gradients for the planning process. They also have implicit function theorem stuff for the contact solve implicit gradients.
Taro: When you put this into practice, they use NVIDIA Warp for GPU-parallel computation and automatic differentiation to build things like forward kinematics and collision detection so they can generate those adjoints.
Rosa: So, the big picture is that this lets you do batched MPC where the intended contacts guide planning but the simulation handles the actual contact realization based on physics.
Dev: The performance numbers are pretty compelling too. they’re showing speedups of factors of twenty-seven and fifty-two on CFree and IMPACT for trajectory generation time.
Taro: That throughput number is significant, reaching two point zero six times ten to the fifth world-steps per second on a single NVIDIA RTX four thousand ninety with two thousand forty-eight parallel worlds.
Rosa: It supports real-world trajectory replay on a Sharpa hand too, which shows they’re thinking about how this moves toward sim-to-real deployment in practice.
Dev: The main things to watch are the limitations they flagged. First, the planner still needs suitable initializations and intended contacts specified before planning, especially for long horizon tasks with lots of contact transitions.
Taro: And second, you still have to tune those objective weights across different tasks and robot embodiments even though they’re shared within a family.
Rosa: Plus, they need to show what happens when the world misbehaves during these contact-rich tasks because that’s where we need more testing.
Dev: So this framework is a strong tool for generating data and testing control policies under these kinds of constraints, but it's not a complete solution yet.
Taro: That’s where we leave it for now on CQDC-Warp, but the next thing to look at is how this kind of simulation helps us test those failure modes in the real world.
The paper's improvements: Taro: So, we've seen how CQDC-Warp uses those intended contacts to guide planning while the simulation figures out the physics of what happens at contact, and now we’re looking at what they actually improved over their previous work.
Rosa: They introduced two main modeling enhancements over the original CQDC. First, they have this compliant-contact formulation which replaces the hard-contact model with something that penalizes contact deformation while still keeping the friction cone constraint active using a parameter rho for softness.
Dev: That softening parameter rho is key because it lets you tune how much compliance you want in your simulation, which is useful for testing different contact dynamics.
Taro: And the second improvement was that revised inertial model we talked about earlier; they made sure there’s coupling between actuated and unactuated coordinates by including inertia at both places through that full generalized mass matrix.
Rosa: That coupling fix is important because it solves those old discrepancies where just moving a joint causes the base of the robot to move due to inertia. It makes the whole system more physically realistic during motion planning.
Dev: So, for a controls engineer, that means when we're running real-time loops, we’re less likely to hit these weird inertial coupling failure modes when simulating complex contact tasks.
Taro: I think it really helps with graceful recovery because if the simulation is more accurate about how the robot moves under load, the resulting policy should be more robust when things go wrong in reality.
Rosa: And for someone running trials outside the lab, this means we can test longer durations because the dynamics are more consistent across different contact situations.
Dev: The way they handled those gradients through geometry surrogates and implicit contact-solve gradients is also a big part of making it differentiable enough for the MPC framework to actually work efficiently.
Taro: That’s how they get that gradient information needed for the optimization, which is crucial because it lets you optimize these contact decisions in real-time.
Rosa: The implication here is that we can use this simulator to train policies on a much wider range of robot platforms and complex tasks without having to painstakingly model every single physical detail from scratch.
Dev: It means the simulation becomes a much better testing ground for the actual control loops, which is what we care about when building a safety case.
Taro: And since they showed this works across LEAP, Shadow hands, and even whole-body control robots, it suggests these contact-guided methods are becoming versatile enough to be applied broadly.
Rosa: We need to remember the caveats too; the planner still needs good initializations and intended contacts specified beforehand, which means we still need more work on long horizon tasks with frequent contact transitions.
Dev: And they also noted that objective weights still require tuning for different tasks and robot bodies, even if you share parameters across a family.
Taro: So it’s powerful for generating data, but the planner isn't completely hands-off yet; you still have to provide some guidance upfront.
Conclusion: Rosa: So, to wrap up, CQDC-Warp is this GPU simulator that lets you plan robot movements while also planning what contacts should happen, using those intended contacts to guide the whole thing through soft costs in your objective function.
Dev: Right, it’s like you tell the planner "Link A needs to touch Object B," and then the simulation handles figuring out exactly where they touch and how much friction there is based on its own dynamics.
Taro: What I find really cool is that they didn't just use a hard contact model. They introduced this compliant-contact formulation which penalizes contact deformation while still keeping the friction cone constraint active using a parameter rho for softness.
Rosa: And they also fixed a big problem with the inertial model; they made sure there’s coupling between actuated and unactuated coordinates by including inertia at both places through that full generalized mass matrix.
Dev: That coupling fix is important because in older models, you could get this weird effect where just moving a joint causes the base of the robot to move due to inertia. This revised metric lets them do an analytical reduction, expressing the optimal control input as a function of delta q through cone projection.
Taro: And they give you implicit derivatives for those contact dynamics and sampled surrogate derivatives for contact geometry so you can compute gradients for the optimization loop.
Rosa: They use witness-pair surrogates to estimate things like the contact midpoint c, unit normal n, and signed gap g together using some shared distance-based softmax weights for the geometry gradients.
Dev: That’s how they get those geometry derivatives, which is one of the two main ways they compute gradients for the planning process. They also have implicit function theorem stuff for the contact solve implicit gradients.
Taro: When you put this into practice, they use NVIDIA Warp for GPU-parallel computation and automatic differentiation to build things like forward kinematics and collision detection so they can generate those adjoints.
Rosa: So, the big picture is that this lets you do batched MPC where the intended contacts guide planning but the simulation handles the actual contact realization based on physics.
Dev: The performance numbers are pretty compelling too. they’re showing speedups of factors of twenty-seven and fifty-two on CFree and IMPACT for trajectory generation time.
Taro: That throughput number is significant, reaching two point zero six times ten to the fifth world-steps per second on a single NVIDIA RTX four thousand ninety with two thousand forty-eight parallel worlds.
Rosa: It supports real-world trajectory replay on a Sharpa hand too, which shows they’re thinking about how this moves toward sim-to-real deployment in practice.
Dev: The main things to watch are the limitations they flagged. First, the planner still needs suitable initializations and intended contacts specified before planning, especially for long horizon tasks with lots of contact transitions.
Taro: And second, you still have to tune those objective weights across different tasks and robot embodiments even though they’re shared within a family.
Rosa: Plus, they need to show what happens when the world misbehaves during these contact-rich tasks because that’s where we need more testing.
Dev: So this framework is a strong tool for generating data and testing control policies under these kinds of constraints, but it's not a complete solution yet.
Taro: That’s where we leave it for now on CQDC-Warp, but the next thing to look at is how this kind of simulation helps us test those failure modes in the real world.
More episodes
- The day in conversation for 2026-10-19
- Daily Summary for 2026-10-19
- 2610.20797-Navigation by Skylight Polarization via Rayleigh Model Optimization on a Road Vehicle
- 2610.20815-A Response to the Reply Presented in "A Fair Comparison of Sliding-Mode and Immersion-and-Invariance Observers"
- 2610.20807-Robotic-supported Data Loss Detection in Android Applications
- 2610.20767-Multi-Dimensional Featurizers Yield Fine-Grained Steering and Data Shaping in Robot Foundation Models
- 2610.20895-BubbleCBF: Constructive Safety for Second-Order Systems and Application to Autonomous Robot Navigation in Unknown Environments
- 2610.20873-Race to Recovery: Designing Advanced Safety Systems for Autonomous Driving via Learning Autonomous Racing
- 2610.20998-MoCaP: Mobile Code-as-Policies for Zero-Shot Mobile Manipulation from Passive Human Videos
- 2610.20970-When a robot can trust reversing wind to find an odour source