Developing exoskeleton controllers that generalize across diverse locomotor conditions typically requires extensive motion-capture data and biomechanical labeling, limiting scalability beyond instrumented laboratory settings. Here, we present a physics-based neuromusculoskeletal learning framework that trains a hip-exoskeleton control policy entirely in simulation, without motion-capture demonstrations, and deploys it on hardware via policy distillation. A reinforcement learning teacher policy is trained using a muscle-synergy action prior over a wide range of walking speeds and slopes through a two-stage curriculum, enabling direct comparison between assisted and no-exoskeleton conditions. In simulation, exoskeleton assistance reduces mean muscle activation by up to 3.4% and mean positive joint power by up to 7.0% on level ground and ramp ascent, with benefits increasing systematically with walking speed. On hardware, the assistance profiles learned in simulation are preserved across matched speed-slope conditions (r: 0.82, RMSE: 0.03 Nm/kg), providing quantitative evidence of sim-to-real transfer without additional hardware tuning. These results demonstrate that physics-based neuromusculoskeletal simulation can serve as a practical and scalable foundation for exoskeleton controller development, substantially reducing experimental burden during the design phase.
Figures & tables
Fig. 1: Simulation-to-real workflow for learning and deploying hip-exoskeleton control policy. Left (Neuromusculoskeletal simulation): a privileged teacher policy is trained in a predictive human simulation across randomized target speeds ( 0.7 – 1.5 m/s) and slopes ( −5∘ to +5∘ ) using a two-stage curriculum: Curriculum 1 learns stable locomotion without exoskeleton actuation (no-exoskeleton), then Curriculum 2 introduces bilateral hip exoskeleton actuation (exoskeleton-assisted), enabling the policy to learn assistive torques. The teacher outputs reduced-dimensional actions via a muscle-synergy prior, comprising lower-limb muscle excitations ( esyn→elimb ), direct trunk-muscle excitations ( etrunk ), and raw bilateral hip-exoskeleton torques ( τraw ), which are smoothed by a first-order low-pass filter (LPF) to yield applied torques ( τfilt ). Simulation state and reward are fed back to optimize the teacher. Right (Real world): the teacher is distilled into a temporal convolutional network student policy that maps a short history of unilateral thigh inertial measurement unit (IMU) gyroscope signals ( θ˙thigh ) to unilateral hip-torque commands for onboard control; the resulting τraw is passed through the onboard LPF and applied to the physical exoskeleton. The dashed red arrows denote the distillation link from simulated rollouts to the deployable IMU-only policy.
Fig. 2: Human-exoskeleton musculoskeletal model used in simulation. Frontal and lateral view of the H2190 full-body musculoskeletal model, comprising 21 degrees of freedom and 90 Hill-type musculotendon actuators. Bilateral hip exoskeleton actuators apply flexion-extension torques about the left and right hip joints.
Parameter
Value
Algorithm
Soft Actor–Critic (SAC)
Implementation
Stable-Baselines3
Policy
Multilayer perceptron
Actor/Critic sizes
[512, 512, 256]
Parallel envs
30
Stage 1 steps
0 – 5×107
TABLE I: SAC hyperparameters and curriculum schedule.
Fig. 3: Comparison of thigh gyroscope signals between simulation and hardware during level-ground walking at 1.2 m/s. Simulation signals were computed as femur angular velocity expressed in the femur local frame; hardware signals were recorded from an IMU mounted on the distal thigh bar of the exoskeleton. The rotational axes were defined with respect to an upright standing posture: x corresponds to the anteroposterior axis, y to the vertical axis, and z to the mediolateral axis. Solid lines show stride-averaged means and shaded regions indicate ±1 SD across strides. The mediolateral component (gyro z ) exhibited the strongest agreement ( r=0.55 ) and was therefore selected as the sole input modality for student policy training.
Fig. 4: Teacher-student agreement in hip torque prediction. Representative time series of right hip assistance torque produced by the privileged teacher policy (black) and the distilled student controller (teal). The student is a temporal convolutional network that predicts hip torque from a short history of femur mediolateral gyroscope measurements expressed in the local segment frame. The example trial was collected on level ground while commanding a sequence of target speeds ( 0.7 , 1.1 , and 1.5 m/s), each maintained for 5 s. The first 0.95 s of the trial contains no student output, as the input history window has not yet been fully populated. For this trial, the student closely tracked the teacher output ( R2=0.93 ).
Fig. 5: Robotic hip exoskeleton hardware and onboard software architecture. (a) Robotic hip exoskeleton designed to apply bilateral hip flexion-extension torques during locomotion. (b) Onboard software architecture of the machine learning co-processor. Input sensor signals are logged via a dedicated I/O process; a parallel inference process generates real-time hip torque commands using the pre-trained student policy.
Fig. 6: Joint-level biomechanical comparison between simulation and human walking. Sagittal-plane hip, knee, and ankle joint angles (left column of each slope) and net joint moments normalized by body mass (right column of each slope) are shown over the gait cycle for ramp descent ( −5∘ ), level ground ( 0∘ ), and ramp ascent ( +5∘ ). Colored curves denote simulation results across walking speeds (0.7–1.5 m/s), while the gray band indicates the mean ± standard deviation of open-source human experimental data [ 36 ] . All simulation traces correspond to the curriculum stage 2 no-exoskeleton condition, where exoskeleton torques were clamped to zero.
Fig. 7: Learned assistance torque profile compared with conventional controllers. Representative hip assistance torque profiles over the gait cycle for a spline-based controller (Spline Control, red; 1.25 m/s) [ 37 ] , a biological-torque-based controller (Bio torque, green; 1.2 m/s) [ 15 ] , and the learned simulation policy (Simulation, blue; 1.2 m/s). Each profile was normalized by its peak hip flexion torque magnitude for visualization and peak-timing comparison. Markers indicate the timing of the hip extension peak (minimum, triangles) and hip flexion peak (maximum, circles).
Fig. 8: Peak-timing comparison between the biological hip joint moment and the learned hip-exoskeleton assistance profile. (a) Representative gait-cycle waveforms during level-ground walking ( 0∘ , 1.2 m/s), showing the sagittal-plane biological hip joint moment (black) and exoskeleton assistive torque (blue), with markers indicating the extension and flexion peaks. (b) Distribution of extension-peak and flexion-peak timings (percentage of the gait cycle) for the hip joint moment and exoskeleton torque across trials.
Fig. 9: Effect of exoskeleton assistance on muscle activation and positive joint power across conditions. Summary metrics are shown as a function of walking speed for three slope conditions (ramp descent, level ground, and ramp ascent), related to no-exoskeleton conditions. The top row reports mean muscle activation averaged across all muscles over a 15 s walking segment; the bottom row reports mean positive joint power averaged across bilateral hip, knee, and ankle joints in the sagittal plane.
Fig. 10: Comparison of hip assistance torque profiles between simulation and hardware. Gait-cycle-normalized right hip assistance torque is shown for ramp descent (RD), level ground (LG), and ramp ascent (RA). Blue curves denote simulated torque profiles across walking speeds from 0.7 to 1.5 m/s. Red curves show torque profiles generated by the distilled student controller on the physical device, evaluated at discrete speeds from 0.7 to 1.5 m/s in 0.2 m/s increments. Color intensity increases with walking speed.