Embodied Human Simulation for Quantitative Design and Analysis of Interactive Robotics

TL;DR

This work establishes embodied human simulation as a scalable paradigm for interactive robotics design by simulating the coupled human-robot system, offering a systematic way to concurrently co-optimize a robot's structural parameters and control policy.

Abstract

Physical interactive robotics, ranging from wearable devices to collaborative humanoid robots, require close coordination between mechanical design and control. However, evaluating interactive dynamics is challenging due to complex human biomechanics and motor responses. Traditional experiments rely on indirect metrics without measuring human internal states, such as muscle forces or joint loads. To address this issue, we develop a scalable simulation-based framework for the quantitative analysis of physical human-robot interaction. At its core is a full-body musculoskeletal model serving as a predictive surrogate for the human dynamical system. Driven by a reinforcement learning controller, it generates adaptive, physiologically grounded motor behaviors. We employ a sequential training pipeline where the pre-trained human motion control policy acts as a consistent evaluator, making large-scale design space exploration computationally tractable. By simulating the coupled human-robot system, the framework provides access to internal biomechanical metrics, offering a systematic way to concurrently co-optimize a robot's structural parameters and control policy. We demonstrate its capability in optimizing human-exoskeleton interactions, showing improved joint alignment and reduced contact forces. This work establishes embodied human simulation as a scalable paradigm for interactive robotics design.

Figures (7)

• Figure 1: Demonstrations of the Digital Human Embodiment with interactive simulation framework. (a) The co-optimization pipeline applied to a wearable exoskeleton, the human policy is trained and applied to robotics optimization. (b) The framework's scalability to diverse interactive robotics, illustrated by a daily collaborative task with a humanoid robot.
• Figure 2: The human-exoskeleton coupled simulation model. (a) The full system, showing the MS-Human-700 model wearing the OpenExo exoskeleton. (b) A detailed view of the human-robot physical interface, where the adduction joint is passive to allow adaptation to natural leg abduction movement. (c) Illustration of the adjustable structural parameters in the interactive system. Local cuff adjustments correspond to vertical shifts at binding sites, while global assembly module adjustments modify the axis relative to attached human body segment through rotations and translations. The interactions were modeled with compliant elastic tendons elements (gray spheres) that connect the exoskeleton shell to the underlying human segment and transmit assistive forces.
• Figure 3: Flowchart of the co-optimization loop. The optimization algorithm proposes a new set of exoskeleton parameters (control and structure). The human-robot coupled simulation is executed for several gait cycles with these parameters. The resulting motion and physiological data are used to evaluate the cost function. This cost is returned to the optimizer, which updates its internal model and proposes the next set of parameters to evaluate, iterating until convergence.
• Figure 4: Validation of the digital human's walking kinematics. The simulated joint angles (blue solid lines) are compared with the reference motion capture data (green dashed lines) for major joints. The shaded areas represent the standard deviation across 10 simulation trials with different initialization time steps, indicating high consistency.
• Figure 5: Robustness evaluation of the learned walking policy against external perturbations. (a) The agent's dynamic recovery sequence after an impulsive force is applied to the pelvis, where the pink arrow visualizes the applied disturbance force. (b) Recovery time across a range of forward (positive) and backward (negative) forces. The red line shows the time required to recover the original kinematic trajectory, demonstrating a fast return to the intended motion. The shaded area represents the standard deviation over 10 randomly initialized trials.