pith. sign in

arxiv: 2605.13192 · v1 · pith:ABTDNUPZnew · submitted 2026-05-13 · 💻 cs.RO

Dynamics Computation of Soft-Rigid Hybrid-Link System and Its Application to Motion Analysis of an Athlete Wearing Sport Prosthesis

Pith reviewed 2026-05-14 18:16 UTC · model grok-4.3

classification 💻 cs.RO
keywords systemforcemotionhybrid-linkprosthesisanalysisdynamicsestimation
0
0 comments X

The pith

A soft-rigid hybrid-link model enables inverse kinematics and dynamics computation for athletes wearing flexible prostheses, with 12% error reported on ground-reaction-force estimation in a human-subject trial.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The work combines two modeling traditions. Human bodies are usually treated as rigid links connected at joints, while flexible prostheses bend and store energy like soft rods. The authors merge both into one hybrid system so that the deformation of the prosthesis and the forces it exchanges with the remaining limb can be calculated together. From standard motion-capture markers they first solve an inverse-kinematics problem to recover the full pose of the hybrid chain. They then run inverse dynamics to obtain joint torques and ground-reaction forces. In a single-subject experiment the computed ground-reaction force matched measured force-plate data within roughly 12 percent. The same framework also supplies estimates of muscle forces that remain after amputation and that interact with the bending prosthesis.

Core claim

Through a human subject experiment, we show that the inverse dynamics achieved approximately 12% error on the ground reaction force estimation. Furthermore, we provide the muscle force estimation considering muscle amputation and interaction force with the prosthesis leg deformation.

Load-bearing premise

The hybrid-link formulation accurately captures the interaction force between the rigid human skeleton and the flexible prosthesis without requiring additional parameters or post-hoc tuning that would invalidate the inverse-dynamics results.

Figures

Figures reproduced from arXiv: 2605.13192 by Ko Yamamoto, Sunghee Kim, Taiki Ishigaki, Yuta Shimane.

Figure 1
Figure 1. Figure 1: Overview for the motion analysis framework based on the hybrid-link system. From a motion capture system, we reconstruct human motion and prosthesis deformation by inverse kinematics. Then, inverse dynamics calculates the joint torques, inner force in the prosthesis and the ground reaction force. From the obtained joint torques and interaction force from the prosthesis, we estimate muscle tendon forces by … view at source ↗
Figure 2
Figure 2. Figure 2: Schematic view of the PCS model H(s, t) ∈ SE(3) as H(s, t) =  R(s, t) p(s, t) 0 T 1  (4) where R(s, t) ∈ SO(3) is the rotational matrix and p(s, t) ∈ R 3 is the position vector. Hereinafter, t or s is omitted for simplicity when its meaning is obvious. Given the configuration curve H(s), a six-dimensional displacement ξ(s) ∈ R 6 is defined with the operator [ξ(s)×] as [ξ(s)×] := H(s) −1 ∂H ∂s . (5) In th… view at source ↗
Figure 3
Figure 3. Figure 3: Conceptual illustration of a hybrid-link system that consists of rigid-links and soft-segments rep￾resented by the PCS model. We consider that the base-link is not fixed to the environment and that there are multiple contacts where f C,i denote a contact force. generalized velocity vector ψ as follows: q := n H0, [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schemetics for the transtibial amputation surgery by myodesis focused on muscles 2.4. Muscle Force Estimation 2.4.1. Muscle Force Optimization From the joint torques estimated by the ID, we can estimate muscle forces that actuate each joint. Because a human joint is redundantly actuated by multiple muscles, muscle forces f are usually estimated by an optimization given as follows [2]: min f 1 2 ∥τ R − J T … view at source ↗
Figure 5
Figure 5. Figure 5: Schematic of the side view of the MRI scanning for leg with transtibial amputation operated by myodesis. Gastrocnemius and Soleus [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Arrangement of the transected muscles based on the MRI image of the subject. • Anterior compartment musculature was removed. • Deep posterior musculature (Gastrocnemius and Soleus) was transected and se￾cured to the tibia via small holes. If the transtibial amputation is operated by myodesis, the side view of the MRI scan￾ning of the amputated leg can be illustrated as [PITH_FULL_IMAGE:figures/full_fig_p0… view at source ↗
Figure 7
Figure 7. Figure 7: Hybrid-link model that integrates human skeleton and a leaf-spring prosthesis leg [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Layout of optical markers for motion capture measurement. than the amputation stump. Also, the end of the muscles are secured to the front of the tibia. Therefore, these muscles has the shape as wrapping the tibia stump [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Walking motion reconstructed by the inverse kinematics of the hybrid-link system. 1st Force Plate Right (Prosthesis) 3rd Force Plate Right (Prosthesis) 2nd Force Plate Left (Foot) [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: A subject was asked to walk or run with a single stride on each force plate in the experiment The number of segments was determined based on the previous report [23], where the change ratio of strain along the central axis coordinate was large. We assume that each segment has only angular strain with 3 DOFs ignoring linear strain while qs,0 is calculated so that each segment has 6 DOFs for a precise model… view at source ↗
Figure 11
Figure 11. Figure 11: Setting of contact points in the inverse dynamics [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Reconstructed walking motion and estimated ground reaction force. for the prosthesis, as shown in [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Reconstructed running motion and estimated ground reaction force. 0 0.5 1 1.5 2 2.5 time(s) -500 0 500 1000 1500 F x(N) 0 0.5 1 1.5 2 2.5 time(s) -500 0 500 1000 1500 F y(N) Contact force for 1st force plate 0 0.5 1 1.5 2 2.5 time(s) -500 0 500 1000 1500 F z(N) Force plate ID Result 0 0.5 1 1.5 2 2.5 time(s) -500 0 500 1000 1500 F x(N) 0 0.5 1 1.5 2 2.5 time(s) -500 0 500 1000 1500 F y(N) Contact force fo… view at source ↗
Figure 14
Figure 14. Figure 14: Estimated contact force and its ground truth along xyz axes in the walking motion. W2 = O. This would degrade the precision of the estimation; however, we excluded the EMG signal factor from the estimation for a pure comparison as a trial. During the experiment, the subject wore tights to protect the connecting part be￾tween the socket and amputated leg. If the air enters inside the tights, the subject co… view at source ↗
Figure 15
Figure 15. Figure 15: Estimated contact force and its ground truth along xyz axes in the running motion [PITH_FULL_IMAGE:figures/full_fig_p015_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Visualization of estimated muscle activities, in which the red part indicates higher activity, and the green part indicates lower activity. is the most affected joint from the transtibial muscle amputation. If the change of transected muscle tensions affects the knee joint torque, the action of these muscles will be altered together. Therefore, we consider that the evaluation of these two muscles is still… view at source ↗
Figure 17
Figure 17. Figure 17: Location of the muscle and EMG sensors for the experiment 0 0.2 0.4 0.6 0.8 1 Time [s] -500 -400 -300 -200 -100 0 Muscle Tension [N ] Biceps Femoris Caput Longum (Right) emg msid Single Phase (Right) Single Phase (Left) Single Phase (Right) [PITH_FULL_IMAGE:figures/full_fig_p016_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Muscle tension comparison between the estimations by the optimization and the EMG signal: Biceps Femoris Caput Longum 4. Discussions 4.1. Ground Reaction Force Estimation Tables 1 and 2 show the root mean square error (RMSE) and related root mean square error (rRMSE) between the estimated and ground truth values, which are defined as follows [24]: RMSE = sPn i=1(fi − ˆfi) 2 n , and (38) rRMSE = RMSE × 100… view at source ↗
Figure 19
Figure 19. Figure 19: Muscle tension comparison between the estimations by the optimization and the EMG signal: Semitendinosus [PITH_FULL_IMAGE:figures/full_fig_p017_19.png] view at source ↗
Figure 11
Figure 11. Figure 11: (2) Running motion result: As seen in [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗
read the original abstract

This paper presents a motion analysis framework for an athlete wearing sport-specific flexible prosthesis based on the soft-rigid hybrid-link system. Such a motion analysis is a challenging problem because we need to consider the interaction force between the rigid human skeleton system and a flexible prosthesis. However, most of human musculoskeletal models are based on the computation framework of a rigid-body multi-link system. Recently in soft robotics research field, fast and efficient modeling methods were developed for a flexible rod deformation, which allows us to build a hybrid-link system that integrates rigid-link and soft-bodies in a unified formulation. We apply inverse kinematics of the hybrid-link system to motion reconstruction from a motion captured data, and also present the estimation of the joint torques and ground reaction force by inverse dynamics. Through a human subject experiment, we show that the inverse dynamics achieved approximately 12% error on the ground reaction force estimation. Furthermore, we provide the muscle force estimation considering muscle amputation and interaction force with the prosthesis leg deformation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. This paper presents a framework for motion analysis of an athlete with a sport prosthesis using a soft-rigid hybrid-link system. It integrates rigid-body dynamics with soft-body modeling to account for interaction forces, applies inverse kinematics to motion capture data, and uses inverse dynamics to estimate joint torques, ground reaction forces (GRF), and muscle forces considering amputation and prosthesis deformation. A human subject experiment reports an approximately 12% error in GRF estimation.

Significance. Should the hybrid model prove robust, the work offers a novel approach to analyzing sports movements in prosthetic users by unifying soft robotics techniques with human dynamics. It enables estimation of internal forces and muscle activations that rigid models cannot capture, potentially aiding in performance optimization and injury prevention for amputee athletes. The experimental application to real motion data is a positive step toward practical use.

major comments (3)
  1. [Human subject experiment (results)] The GRF estimation achieves approximately 12% error, but the abstract and results provide no error bars, no description of data exclusion rules, no baseline comparison to rigid models, and no count of free parameters fitted (such as prosthesis stiffness and damping coefficients). This limits assessment of the result's reliability and the hybrid formulation's specific role.
  2. [Hybrid-link formulation (methods)] The interaction force computation between the rigid skeleton and flexible prosthesis depends on stiffness and damping parameters. The paper should specify their identification method and confirm they were not adjusted post-hoc to fit the force-plate data used in validation, as this could introduce circularity in the 12% error claim.
  3. [Inverse dynamics application] While GRF error is reported, there is no independent validation of the soft-body deformation predictions, such as comparison to measured prosthesis bending or interface pressures. The integrated GRF metric alone does not isolate the accuracy of the hybrid coupling terms.
minor comments (2)
  1. [Introduction] The background on soft robotics methods for flexible rod deformation could include more specific citations to recent efficient modeling techniques.
  2. [Figures] Ensure that diagrams of the hybrid-link system clearly label rigid and soft segments and the interaction points.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major point below, indicating where revisions will be made to improve clarity and rigor.

read point-by-point responses
  1. Referee: [Human subject experiment (results)] The GRF estimation achieves approximately 12% error, but the abstract and results provide no error bars, no description of data exclusion rules, no baseline comparison to rigid models, and no count of free parameters fitted (such as prosthesis stiffness and damping coefficients). This limits assessment of the result's reliability and the hybrid formulation's specific role.

    Authors: We agree that these details are important for evaluating reliability. In the revised manuscript we will add error bars computed from the temporal variation of the GRF estimation error across the movement cycle, describe the data processing pipeline and confirm that no trials were excluded, include a direct comparison against a standard rigid-body inverse-dynamics baseline, and explicitly state the number of fitted parameters (two stiffness and two damping coefficients identified in separate calibration tests). revision: yes

  2. Referee: [Hybrid-link formulation (methods)] The interaction force computation between the rigid skeleton and flexible prosthesis depends on stiffness and damping parameters. The paper should specify their identification method and confirm they were not adjusted post-hoc to fit the force-plate data used in validation, as this could introduce circularity in the 12% error claim.

    Authors: The stiffness and damping coefficients were obtained from independent material characterization experiments performed on the prosthesis prior to the motion-capture session and were held fixed during inverse-dynamics computation. They were not tuned to the force-plate recordings. We will expand the Methods section with a full description of the calibration procedure to eliminate any ambiguity regarding circularity. revision: yes

  3. Referee: [Inverse dynamics application] While GRF error is reported, there is no independent validation of the soft-body deformation predictions, such as comparison to measured prosthesis bending or interface pressures. The integrated GRF metric alone does not isolate the accuracy of the hybrid coupling terms.

    Authors: We acknowledge that the present validation relies solely on the integrated GRF metric and does not include direct measurements of prosthesis bending or interface pressure. No such auxiliary sensor data were collected in the current human-subject experiment, so an independent check of the soft-body deformation component cannot be added without new experiments. revision: no

standing simulated objections not resolved
  • Independent validation of soft-body deformation (prosthesis bending or interface pressures) is not possible because the original experiment did not record these quantities.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The hybrid-link formulation implicitly relies on standard rigid-body dynamics assumptions plus a soft-rod deformation model whose stiffness and damping parameters are not specified; no new entities are postulated.

free parameters (1)
  • prosthesis stiffness and damping coefficients
    Required to close the soft-rod deformation equations; values not reported in abstract and therefore treated as fitted or chosen parameters.
axioms (1)
  • domain assumption The interaction force between rigid skeleton and flexible prosthesis can be captured by a single unified hybrid-link formulation without additional contact constraints.
    Invoked when the authors state that the hybrid system integrates rigid and soft bodies for inverse dynamics.

pith-pipeline@v0.9.0 · 5484 in / 1394 out tokens · 38211 ms · 2026-05-14T18:16:42.991145+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    The fastest sprinter in 2068 has an artificial limb? ; 2015

    Hobara H, Kobayashi Y, Heldoorn TA, et al. The fastest sprinter in 2068 has an artificial limb? ; 2015

  2. [2]

    Somatosensory computation for man-machine interface from motion-capture data and musculoskeletal human model

    Nakamura Y, Yamane K, Fujita Y, et al. Somatosensory computation for man-machine interface from motion-capture data and musculoskeletal human model. IEEE Transactions on Robotics. 2005;21(1):58–66

  3. [3]

    Opensim: open-source software to create and analyze dynamic simulations of movement

    Delp SL, Anderson FC, Arnold AS, et al. Opensim: open-source software to create and analyze dynamic simulations of movement. IEEE transactions on biomedical engineering. 2007;54(11):1940–1950

  4. [4]

    A prosthesis-specific multi-link segment model of lower- limb amputee sprinting

    Rigney SM, Simmons A, Kark L. A prosthesis-specific multi-link segment model of lower- limb amputee sprinting. Journal of Biomechanics. 2016;49(14):3185–3193

  5. [5]

    Inverse optimal control based enhancement of sprinting motion analysis with and without running-specific prostheses

    Kleesattel ALEN, Mombaur K. Inverse optimal control based enhancement of sprinting motion analysis with and without running-specific prostheses. In: 2018 7th IEEE Inter- national Conference on Biomedical Robotics and Biomechatronics (Biorob); IEEE; 2018. p. 556–562

  6. [6]

    Can forward dynamics simulation with sim- ple model estimate complex phenomena?: Case study on sprinting using running-specific prosthesis

    Murai A, Hobara H, Hashizume S, et al. Can forward dynamics simulation with sim- ple model estimate complex phenomena?: Case study on sprinting using running-specific prosthesis. Robomech Journal. 2018;5(1):1–8

  7. [7]

    Control of elastic soft robots based on real-time finite element method

    Duriez C. Control of elastic soft robots based on real-time finite element method. In: 2013 IEEE international conference on robotics and automation; IEEE; 2013. p. 3982–3987

  8. [8]

    A system for acquisition and modelling of ice- hockey stick shape deformation from player shot videos

    Mendhurwar K, Handa G, Zhu L, et al. A system for acquisition and modelling of ice- hockey stick shape deformation from player shot videos. In: Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition Workshops; 2020. p. 890–891

  9. [9]

    Effect of ice hockey stick stiffness on performance

    Kays BT, Smith L V. Effect of ice hockey stick stiffness on performance. Sports Engineer- ing. 2017;20:245–254

  10. [10]

    A study on the dynamic structural behavior of olympic sabres

    Vivaldi F, Cortese L, Coppola T, et al. A study on the dynamic structural behavior of olympic sabres. Procedia Structural Integrity. 2018;8:345–353

  11. [11]

    A framework for the automated design and modelling of soft robotic systems

    Runge G, Raatz A. A framework for the automated design and modelling of soft robotic systems. CIRP Annals. 2017;66(1):9–12

  12. [12]

    Review on generic methods for mechanical modeling, simulation and control of soft robots

    Schegg P, Duriez C. Review on generic methods for mechanical modeling, simulation and control of soft robots. Plos one. 2022;17(1):e0251059

  13. [13]

    Nonlinear problems of elasticity

    Antman SS. Nonlinear problems of elasticity. New York: Springer; 2005

  14. [14]

    Discrete cosserat approach for multisection soft manip- ulator dynamics

    Renda F, Boyer F, Dias J, et al. Discrete cosserat approach for multisection soft manip- ulator dynamics. IEEE Transactions on Robotics. 2018;34(6):1518–1533

  15. [15]

    A geometric and unified approach for modeling soft-rigid multi- body systems with lumped and distributed degrees of freedom

    Renda F, Seneviratne L. A geometric and unified approach for modeling soft-rigid multi- body systems with lumped and distributed degrees of freedom. In: 2018 IEEE Interna- tional Conference on Robotics and Automation (ICRA); IEEE; 2018. p. 1567–1574

  16. [16]

    Dynamics computation of a hybrid multi-link humanoid robot integrating rigid and soft bodies

    Ishigaki T, Yamamoto K. Dynamics computation of a hybrid multi-link humanoid robot integrating rigid and soft bodies. In: 2021 IEEE/RSJ International Conference on Intel- ligent Robots and Systems (IROS); IEEE; 2021. p. 2816–2821

  17. [17]

    Sorosim: A matlab toolbox for hybrid rigid-soft 21 robots based on the geometric variable-strain approach

    Mathew AT, Hmida IB, Armanini C, et al. Sorosim: A matlab toolbox for hybrid rigid-soft 21 robots based on the geometric variable-strain approach. IEEE Robotics and Automation Magazine. 2023 9;30:106–122

  18. [18]

    Inverse kinematics of hybrid multi-link system and its application to motion capture for athlete wearing sports prosthesis

    Kim S, Ishigaki T, Shimane Y, et al. Inverse kinematics of hybrid multi-link system and its application to motion capture for athlete wearing sports prosthesis. In: 2022 IEEE- RAS 21st International Conference on Humanoid Robots (Humanoids); IEEE; 2022. p. 837–842

  19. [19]

    Rigid body dynamics algorithms

    Featherstone R. Rigid body dynamics algorithms. Springer; 2008

  20. [20]

    The heat of shortening and the dynamic constants of muscle

    Hill A V. The heat of shortening and the dynamic constants of muscle. Proceedings of the Royal Society of London Series B-Biological Sciences. 1938;126(843):136–195

  21. [21]

    Myodesis or myoplasty in trans-femoral amputations

    Geertzen JH, van der Schans SM, Jutte PC, et al. Myodesis or myoplasty in trans-femoral amputations. what is the best option? an explorative study. Medical Hypotheses. 2019; 124:7–12

  22. [22]

    Outcomes after 294 transtibial amputations with the posterior myocutaneous flap

    Brown BJ, Iorio ML, Klement M, et al. Outcomes after 294 transtibial amputations with the posterior myocutaneous flap. The international journal of lower extremity wounds. 2014;13(1):33–40

  23. [23]

    Application of piece-wise constant strain model to flexible deformation calculation of sports prosthesis and stiffness estimation

    Shimane Y, Ishigaki T, Kim S, et al. Application of piece-wise constant strain model to flexible deformation calculation of sports prosthesis and stiffness estimation. In: 2022 IEEE Int. Conf. Intell. Robots Syst.; (To be published)

  24. [24]

    Prediction of ground reaction forces during gait based on kinematics and a neural network model

    Oh SE, Choi A, Mun JH. Prediction of ground reaction forces during gait based on kinematics and a neural network model. Journal of biomechanics. 2013;46(14):2372–2380

  25. [25]

    Impedance characteristics of a neuromusculoskeletal model of the human arm i

    Stroeve S. Impedance characteristics of a neuromusculoskeletal model of the human arm i. posture control. Biological cybernetics. 1999;81(5-6):475–494

  26. [26]

    Analysis of fundamental human movement patterns through the use of in-depth antagonistic muscle models

    Winters JM, Stark L. Analysis of fundamental human movement patterns through the use of in-depth antagonistic muscle models. IEEE transactions on biomedical engineering. 1985;(10):826–839. 22