We present CRISP (Contact-RIch Simulation Platform), a high-fidelity physics engine tailored for complex multi-contact simulations such as tight-tolerance robotic manipulation. Achieving high physical fidelity in robotic simulation requires both expressive modeling of geometry and contact interactions, as well as accurate numerical resolution via robust collision detection and contact solvers. However, existing simulators often either rely on limited support for geometric representations and simplified modeling of contact interactions, or employ numerical resolution methods whose accuracy or robustness is inherently constrained. Accordingly, we develop a new simulator that supports diverse geometric representations with accurate optimization-based collision detection, and combines contact modeling with robust augmented Lagrangian-based contact solvers. This integration enables efficient and consistent detection of contact information across complex geometries while accurately resolving multi-contact constraints without problematic relaxations, which is essential for simulating contact-intensive and sharp interactions. We validate the physical fidelity of our simulator against state-of-the-art platforms and further demonstrate its capabilities through complex robotic manipulation scenarios. CRISP is publicly available at https://github.com/INRoL/crisp.
Figures & tables
Physics engine
Geometry 1
Formulation
Solver
ODE [ 40 ]
Mesh
LCP
Direct/PGS
Bullet [ 9 ]
Mesh
LCP
Direct/PGS
DART [ 23 ]
Mesh
LCP
Direct/PGS
MuJoCo [ 44 ]
Mesh 3 /SDF
CCP
Newton
Drake [ 43 ]
Mesh
CCP
Newton
RaiSim [ 19 ]
Mesh
- 2
Bisection
TABLE I : Comparison of existing physics engines and CRISP.
Fig. 2 : High-level architecture of CRISP. The model builder defines the low-level structure model_t , from which data_t is instantiated. At each step, collision detection and constraint solving are performed, and the resulting state is visualized.
Fig. 4 : Visualizations of DSF and TDSF geometries. It is possible to specify the degree of approximation by choosing the smoothness parameter, as described in Appendix -A .
Fig. 5 : Visualization of assembly scenarios. Initial configurations and examples of detected contact points. Left: peg insertion scenario. Right: bolt-nut assembly scenario.
Peg Insertion (Sec. VI-A )
Tol. ( \SIUnitSymbolMicrom )
CRISP (CANAL/SubADMM)
MuJoCo
Isaac Sim
DSF–TDSF
SDF–SDF
SDF–SDF
Mesh–SDF
200
43.649
80.565
147.61
143.89
756.31
47.701
100
25.646
124.46
107.95
90.769
1349.9
22.551
50
37.844
43.503
79.334
106.09
1328.6
-
Bolt-Nut Assembly (Sec. VI-B )
TABLE II : Average penetration depth across assembly scenarios. Units are in \SIUnitSymbolMicrom , and ommitted entries indicate failed assemblies.
Fig. 6 : Top-heavy stacking scenario. Left: Visualization of the scenario. Right: Height of the top block over time for mass ratios of 10 (top) and 100 (bottom).
Fig. 7 : Oblique sliding scenario. Left: Trajectories for the sliding case in CRISP (green) and MuJoCo (blue). Right: Height of the block over time for sliding case (top) and sticking case at the critical friction coefficient (bottom).
Fig. 8 : Joint constraint scenario. Joint angle trajectories over time (left) and a magnified view near the final time step (right).
Fig. 9 : Results of gear insertion and driving simulation. Number of contacts and solver residuals as a function of time.
Plane
Sphere
Capsule
Box
Cylinder
Convex
Mesh
SDF
DSF
TDSF
Plane
–
SF
GD
SF
Sphere
Analytic
BF
Direct
RTR
Capsule
Search
×
Box
MPR
FW
Cylinder
GD
Convex
MPR
TABLE III : Collision function table in CRISP. Each cell indicates the collision detection algorithm used for the corresponding geometry pair (–: cannot be defined, × : not supported). Cell colors indicate a qualitative assessment of algorithmic performance (green: highly reliable, blue: recommended, yellow: neutral, red: use with caution).
Physics simulation for contact-rich robotics is often bottlenecked by contact resolution: mainstream engines enforce non-penetration and Coulomb friction via complementarity constraints or constrained optimization, requiring per-step iterative solves whose cost grows superlinearly with contact density. We present ComFree-Sim, a GPU-parallelized analytical contact physics engine built on complementarity-free contact modeling. ComFree-Sim computes contact impulses in closed form via an impedance-style prediction--correction update in the dual cone of Coulomb friction. Contact computation decouples across contact pairs and becomes separable across cone facets, mapping naturally to GPU kernels and yielding near-linear runtime scaling with the number of contacts. We further extend the formulation to a unified 6D contact model capturing tangential, torsional, and rolling friction, and introduce a practical dual-cone impedance heuristic. ComFree-Sim is implemented in Warp and exposed through a MuJoCo-compatible interface as a drop-in backend alternative to MuJoCo Warp (MJWarp). Experiments benchmark penetration, friction behaviors, stability, and simulation runtime scaling against MJWarp, demonstrating near-linear scaling and 2--3 times higher throughput in dense contact scenes with comparable physical fidelity. We deploy ComFree-Sim in real-time MPC for in-hand dexterous manipulation on a real-world multi-fingered LEAP hand and in dynamics-aware motion retargeting, demonstrating that low-latency simulation yields higher closed-loop success rates and enables practical high-frequency control in contact-rich tasks.
Chetan Borse, Zhixian Xie, Wei-Cheng Huang +1
Intelligent Robotics and Interactive Systems (IRIS) Lab, Arizona State University · Siebel School of Computing and Data Science, University of Illinois at Urbana-Champaign
It has been difficult to accurately model the physics of the interaction between a robot and the objects it manipulates, including contact mechanics, coupled contact dynamics, deformations, stiffness and friction. We present RoboSimStudio (RSS), an interactive toolkit that makes choices of solver, contact model, and physical parameters explicit and editable, allowing both humans and coding agents to quickly create and iteratively refine simulations based on observed physical behavior. We implement RSS with the open Newton physics simulation system and describe experiments with human-guided real2sim using (1) visual evaluation, (2) force-feedback and (3) agentic Real2Sim2Real. Code, video and additional information can be found on the project website: https://robosimstudio.github.io
We propose an optimization-based framework for robust contact-rich manipulation. Recent contact-implicit methods enable online hybrid planning across contact modes, allowing closed-loop manipulation for a given target state and contact location sequence of the robot and object. However, most existing approaches lack the ability to autonomously reason and generate diverse contact location sequences and manipulation trajectories, i.e., active contact location selection, which limits their applicability to relatively simple tasks. Active contact location selection is challenging due to complementarity in contact dynamics and the sparse gradients, making the design of a unified framework for contact selection and planning difficult. To address these challenges, we introduce Simultaneous Contact Selection and Planning (SCSP), a cascaded optimization framework comprising Contact Selection Optimization (CSO) and Contact Planning Optimization (CPO). CSO leverages a surrogate contact model and discrete-continuous optimization to efficiently resolve the nonsmoothness and coupling in contact selection, enabling online global searching of optimal contact locations. CPO performs prior-guided contact planning by evaluating the reference contact locations produced by CSO and generating corresponding manipulation trajectories in real time for redundant manipulators. Extensive simulations and real-world experiments demonstrate that SCSP produces diverse manipulation behaviors and robust control under inaccurate dynamics and perceptual noise. We further validate the generalization of the framework on challenging manipulation tasks. Project website: https://sites.google.com/view/scsp-robot.