1
|
Bahdasariants S, Barela AMF, Gritsenko V, Bacca O, Barela JA, Yakovenko S. Does joint impedance improve dynamic leg simulations with explicit and implicit solvers? PLoS One 2023; 18:e0282130. [PMID: 37399198 DOI: 10.1371/journal.pone.0282130] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 06/15/2023] [Indexed: 07/05/2023] Open
Abstract
The nervous system predicts and executes complex motion of body segments actuated by the coordinated action of muscles. When a stroke or other traumatic injury disrupts neural processing, the impeded behavior has not only kinematic but also kinetic attributes that require interpretation. Biomechanical models could allow medical specialists to observe these dynamic variables and instantaneously diagnose mobility issues that may otherwise remain unnoticed. However, the real-time and subject-specific dynamic computations necessitate the optimization these simulations. In this study, we explored the effects of intrinsic viscoelasticity, choice of numerical integration method, and decrease in sampling frequency on the accuracy and stability of the simulation. The bipedal model with 17 rotational degrees of freedom (DOF)-describing hip, knee, ankle, and standing foot contact-was instrumented with viscoelastic elements with a resting length in the middle of the DOF range of motion. The accumulation of numerical errors was evaluated in dynamic simulations using swing-phase experimental kinematics. The relationship between viscoelasticity, sampling rates, and the integrator type was evaluated. The optimal selection of these three factors resulted in an accurate reconstruction of joint kinematics (err < 1%) and kinetics (err < 5%) with increased simulation time steps. Notably, joint viscoelasticity reduced the integration errors of explicit methods and had minimal to no additional benefit for implicit methods. Gained insights have the potential to improve diagnostic tools and accurize real-time feedback simulations used in the functional recovery of neuromuscular diseases and intuitive control of modern prosthetic solutions.
Collapse
Affiliation(s)
- Serhii Bahdasariants
- Department of Human Performance, School of Medicine, West Virginia University, Morgantown, WV, United States of America
| | - Ana Maria Forti Barela
- Institute of Physical Activity and Sport Sciences, Cruzeiro do Sul University, São Paulo, SP, Brazil
| | - Valeriya Gritsenko
- Department of Human Performance, School of Medicine, West Virginia University, Morgantown, WV, United States of America
- Department of Neuroscience, School of Medicine, West Virginia University, Morgantown, WV, United States of America
| | - Odair Bacca
- Institute of Physical Activity and Sport Sciences, Cruzeiro do Sul University, São Paulo, SP, Brazil
| | - José Angelo Barela
- Department of Physical Education, São Paulo State University, Rio Claro, SP, Brazil
| | - Sergiy Yakovenko
- Department of Human Performance, School of Medicine, West Virginia University, Morgantown, WV, United States of America
- Department of Neuroscience, School of Medicine, West Virginia University, Morgantown, WV, United States of America
- Department of Mechanical and Aerospace Engineering, Benjamin M. Statler College of Engineering and Mineral Resources, West Virginia University, Morgantown, WV, United States of America
- Department of Chemical and Biomedical Engineering, B.M. Statler College of Engineering and Mineral Resources, West Virginia University, Morgantown, WV, United States of America
| |
Collapse
|
2
|
Bahdasariants S, Barela AMF, Gritsenko V, Bacca O, Barela JA, Yakovenko S. Does joint impedance improve dynamic leg simulations with explicit and implicit solvers? BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.09.527805. [PMID: 36798166 PMCID: PMC9934618 DOI: 10.1101/2023.02.09.527805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The nervous system predicts and executes complex motion of body segments actuated by the coordinated action of muscles. When a stroke or other traumatic injury disrupts neural processing, the impeded behavior has not only kinematic but also kinetic attributes that require interpretation. Biomechanical models could allow medical specialists to observe these dynamic variables and instantaneously diagnose mobility issues that may otherwise remain unnoticed. However, the real-time and subject-specific dynamic computations necessitate the optimization these simulations. In this study, we explored the effects of intrinsic viscoelasticity, choice of numerical integration method, and decrease in sampling frequency on the accuracy and stability of the simulation. The bipedal model with 17 rotational degrees of freedom (DOF)-describing hip, knee, ankle, and standing foot contact-was instrumented with viscoelastic elements with a resting length in the middle of the DOF range of motion. The accumulation of numerical errors was evaluated in dynamic simulations using swing-phase experimental kinematics. The relationship between viscoelasticity, sampling rates, and the integrator type was evaluated. The optimal selection of these three factors resulted in an accurate reconstruction of joint kinematics (err < 1%) and kinetics (err < 5%) with increased simulation time steps. Notably, joint viscoelasticity reduced the integration errors of explicit methods and had minimal to no additional benefit for implicit methods . Gained insights have the potential to improve diagnostic tools and accurize real-time feedback simulations used in the functional recovery of neuromuscular diseases and intuitive control of modern prosthetic solutions.
Collapse
Affiliation(s)
- Serhii Bahdasariants
- Department of Human Performance, School of Medicine, West Virginia University, Morgantown, WV, USA
| | - Ana Maria Forti Barela
- Institute of Physical Activity and Sport Sciences, Cruzeiro do Sul University, São Paulo, SP, Brazil
| | - Valeriya Gritsenko
- Department of Human Performance, School of Medicine, West Virginia University, Morgantown, WV, USA
- Department of Neuroscience, School of Medicine, West Virginia University, Morgantown, WV, USA
| | - Odair Bacca
- Institute of Physical Activity and Sport Sciences, Cruzeiro do Sul University, São Paulo, SP, Brazil
| | - José Angelo Barela
- Department of Physical Education, São Paulo State University, Rio Claro, SP, Brazil
| | - Sergiy Yakovenko
- Department of Human Performance, School of Medicine, West Virginia University, Morgantown, WV, USA
- Department of Neuroscience, School of Medicine, West Virginia University, Morgantown, WV, USA
- Department of Mechanical and Aerospace Engineering, Benjamin M. Statler College of Engineering and Mineral Resources, West Virginia University, Morgantown, WV, USA
- Department of Chemical and Biomedical Engineering, B.M. Statler College of Engineering and Mineral Resources, West Virginia University, Morgantown, WV, USA
| |
Collapse
|
3
|
Naemi R, Chatzistergos PE, Chockalingam N. A mathematical method for quantifying in vivo mechanical behaviour of heel pad under dynamic load. Med Biol Eng Comput 2015; 54:341-50. [DOI: 10.1007/s11517-015-1316-5] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2014] [Accepted: 05/18/2015] [Indexed: 11/27/2022]
|
4
|
Farzaneh Y, Akbarzadeh A, Akbari AA. Online bio-inspired trajectory generation of seven-link biped robot based on T–S fuzzy system. Appl Soft Comput 2014. [DOI: 10.1016/j.asoc.2013.05.013] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
5
|
NAEMI ROOZBEH, CHOCKALINGAM NACHIAPPAN. Mathematical Models to Assess Foot–Ground Interaction. Med Sci Sports Exerc 2013; 45:1524-33. [DOI: 10.1249/mss.0b013e31828be3a7] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
6
|
Kao PC, Lewis CL, Ferris DP. Joint kinetic response during unexpectedly reduced plantar flexor torque provided by a robotic ankle exoskeleton during walking. J Biomech 2010; 43:1401-7. [PMID: 20171638 DOI: 10.1016/j.jbiomech.2009.12.024] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2009] [Revised: 11/19/2009] [Accepted: 12/22/2009] [Indexed: 11/29/2022]
Abstract
During human walking, plantar flexor activation in late stance helps to generate a stable and economical gait pattern. Because plantar flexor activation is highly mediated by proprioceptive feedback, the nervous system must modulate reflex pathways to meet the mechanical requirements of gait. The purpose of this study was to quantify ankle joint mechanical output of the plantar flexor stretch reflex response during a novel unexpected gait perturbation. We used a robotic ankle exoskeleton to mechanically amplify the ankle torque output resulting from soleus muscle activation. We recorded lower-body kinematics, ground reaction forces, and electromyography during steady-state walking and during randomly perturbed steps when the exoskeleton assistance was unexpectedly turned off. We also measured soleus Hoffmann- (H-) reflexes at late stance during the two conditions. Subjects reacted to the unexpectedly decreased exoskeleton assistance by greatly increasing soleus muscle activity about 60ms after ankle angle deviated from the control condition (p<0.001). There were large differences in ankle kinematic and electromyography patterns for the perturbed and control steps, but the total ankle moment was almost identical for the two conditions (p=0.13). The ratio of soleus H-reflex amplitude to background electromyography was not significantly different between the two conditions (p=0.4). This is the first study to show that the nervous system chooses reflex responses during human walking such that invariant ankle joint moment patterns are maintained during perturbations. Our findings are particularly useful for the development of neuromusculoskeletal computer simulations of human walking that need to adjust reflex gains appropriately for biomechanical analyses.
Collapse
Affiliation(s)
- Pei-Chun Kao
- School of Kinesiology, University of Michigan, Ann Arbor, MI 48109-2214, USA.
| | | | | |
Collapse
|
7
|
Abstract
SUMMARYThis paper presents a two-level control strategy for bipedal walking mechanism that accounts for implicit control of push-off on the between-step control level and tracking of imposed holonomic constraints on kinematic variables via feedback control on within-step control level. The proposed control strategy was tested in a biologically inspired model with minimal set of segments that allows evolution of human-like push-off and power absorption. We investigated controller's stability characteristics by using Poincaré return map analysis in eight simulation cases and further evaluated the performance of the biped walking model in terms of how variations in torso position and gait velocity relate to push-off and power absorption. The results show that the proposed control strategy, with the same set of controller's gains, enables stable walking in a variety of chosen gait parameters and can accommodate to various trunk inclinations and gait velocities in a similar way as seen in humans.
Collapse
|
8
|
Jo S, Massaquoi SG. A model of cerebrocerebello-spinomuscular interaction in the sagittal control of human walking. BIOLOGICAL CYBERNETICS 2007; 96:279-307. [PMID: 17124602 DOI: 10.1007/s00422-006-0126-0] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2006] [Accepted: 09/05/2006] [Indexed: 05/12/2023]
Abstract
A computationally developed model of human upright balance control (Jo and Massaquoi on Biol cybern 91:188-202, 2004) has been enhanced to describe biped walking in the sagittal plane. The model incorporates (a) non-linear muscle mechanics having activation level -dependent impedance, (b) scheduled cerebrocerebellar interaction for control of center of mass position and trunk pitch angle, (c) rectangular pulse-like feedforward commands from a brainstem/ spinal pattern generator, and (d) segmental reflex modulation of muscular synergies to refine inter-joint coordination. The model can stand when muscles around the ankle are coactivated. When trigger signals activate, the model transitions from standing still to walking at 1.5 m/s. Simulated natural walking displays none of seven pathological gait features. The model can simulate different walking speeds by tuning the amplitude and frequency in spinal pattern generator. The walking is stable against forward and backward pushes of up to 70 and 75 N, respectively, and with sudden changes in trunk mass of up to 18%. The sensitivity of the model to changes in neural parameters and the predicted behavioral results of simulated neural system lesions are examined. The deficit gait simulations may be useful to support the functional and anatomical correspondences of the model. The model demonstrates that basic human-like walking can be achieved by a hierarchical structure of stabilized-long loop feedback and synergy-mediated feedforward controls. In particular, internal models of body dynamics are not required.
Collapse
Affiliation(s)
- Sungho Jo
- Department of Electrical Engineering and Computer Science, Computer Science and Artificial Intelligence Laboratory, Laboratory for Information and Decision Systems, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.
| | | |
Collapse
|
9
|
Zajac FE, Neptune RR, Kautz SA. Biomechanics and muscle coordination of human walking: part II: lessons from dynamical simulations and clinical implications. Gait Posture 2003; 17:1-17. [PMID: 12535721 DOI: 10.1016/s0966-6362(02)00069-3] [Citation(s) in RCA: 223] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Principles of muscle coordination in gait have been based largely on analyses of body motion, ground reaction force and EMG measurements. However, data from dynamical simulations provide a cause-effect framework for analyzing these measurements; for example, Part I (Gait Posture, in press) of this two-part review described how force generation in a muscle affects the acceleration and energy flow among the segments. This Part II reviews the mechanical and coordination concepts arising from analyses of simulations of walking. Simple models have elucidated the basic multisegmented ballistic and passive mechanics of walking. Dynamical models driven by net joint moments have provided clues about coordination in healthy and pathological gait. Simulations driven by muscle excitations have highlighted the partial stability afforded by muscles with their viscoelastic-like properties and the predictability of walking performance when minimization of metabolic energy per unit distance is assumed. When combined with neural control models for exciting motoneuronal pools, simulations have shown how the integrative properties of the neuro-musculo-skeletal systems maintain a stable gait. Other analyses of walking simulations have revealed how individual muscles contribute to trunk support and progression. Finally, we discuss how biomechanical models and simulations may enhance our understanding of the mechanics and muscle function of walking in individuals with gait impairments.
Collapse
Affiliation(s)
- Felix E Zajac
- Rehabilitation R&D Center (153), VA Palo Alto Health Care System and Department of Mechanical Engineering, Stanford University, CA, USA.
| | | | | |
Collapse
|
10
|
Abstract
Recent interest in using modeling and simulation to study movement is driven by the belief that this approach can provide insight into how the nervous system and muscles interact to produce coordinated motion of the body parts. With the computational resources available today, large-scale models of the body can be used to produce realistic simulations of movement that are an order of magnitude more complex than those produced just 10 years ago. This chapter reviews how the structure of the neuromusculoskeletal system is commonly represented in a multijoint model of movement, how modeling may be combined with optimization theory to simulate the dynamics of a motor task, and how model output can be analyzed to describe and explain muscle function. Some results obtained from simulations of jumping, pedaling, and walking are also reviewed to illustrate the approach.
Collapse
Affiliation(s)
- M G Pandy
- Department of Kinesiology, Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, USA.
| |
Collapse
|
11
|
Neptune RR, Wright IC, Van Den Bogert AJ. A Method for Numerical Simulation of Single Limb Ground Contact Events: Application to Heel-Toe Running. Comput Methods Biomech Biomed Engin 2001; 3:321-334. [PMID: 11264857 DOI: 10.1080/10255840008915275] [Citation(s) in RCA: 103] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
The objective of this work was to develop a method to simulate single-limb ground contact events, which may be applied to study musculoskeletal injuries associated with such movements. To achieve this objective, a three-dimensional musculoskeletal model was developed consisting of the equations of motion for the musculoskeletal system, and models for the muscle force generation and ground contact elements. An optimization framework and a weighted least-squares objective function were presented that generated muscle stimulation patterns that optimally reproduced subject-specific movement data. Experimental data were collected from a single subject to provide initial conditions for the simulation and tracking data for the optimization. As an example application, a simulation of the stance phase of running was generated. The results showed that the average difference between the simulation and subject's ground reaction force and joint angle data was less than two inter-trial standard deviations. Further, there was good agreement between the model's muscle excitation patterns and experimentally collected electromyography data. These results give confidence in the model to examine musculoskeletal loading during a variety of landing movements and to study the effects of various factors associated with injury. Limitations were examined and areas of improvement for the model were presented.
Collapse
Affiliation(s)
- R. R. Neptune
- Human Performance Laboratory, University of Calgary, Calgary, AB T2N 1N4
| | | | | |
Collapse
|
12
|
Whittlesey SN, van Emmerik RE, Hamill J. The swing phase of human walking is not a passive movement. Motor Control 2000; 4:273-92. [PMID: 10900056 DOI: 10.1123/mcj.4.3.273] [Citation(s) in RCA: 38] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Many studies have assumed that the swing phase of human walking at preferred velocity is largely passive and thus highly analogous to the swing of an unforced pendulum. In other words, while swing-phase joint moments are generally nonzero during swing, it was assumed that they were either zero or at least negligibly small compared to gravity. While neglect of joint moments does not invalidate a study by default, it remains that the limitations of such an assumption have not been explored thoroughly. This paper makes five arguments that the swing phase cannot be passive, using both original data and the literature: (1) Computer simulations of the swing phase require muscular control to be accurate. (2) Swing-phase joint moments, while smaller than those during stance, are still greater than those due to gravity. (3) Gravity accounts for a minority of the total kinetics of a swing phase. (4) The kinetics due to gravity do not have the pattern needed to develop a normal swing phase. (5) There is no correlation between pendular swing times and human walking periods in overground walking. The conclusion of this paper is that the swing phase must be an actively controlled process, and should be assumed to be passive only when a study does not require a quantitative result. This conclusion has significant implications for many areas of gait research, including clinical study, control theory, and mechanical modeling.
Collapse
Affiliation(s)
- S N Whittlesey
- Department of Exercise Science, Totman Building, University of Massachusetts Amherst, Amherst, MA 01003, USA
| | | | | |
Collapse
|
13
|
Abstract
Significant ground reaction forces exceeding body weight occur during the heel-strike phase of gait. The standard methods of analytical dynamics used to solve the impact problem do not accommodate well the heel-strike collision due to the persistent contact at the front foot and presence of contact at the back foot. These methods can cause a non-physical energy gain on the order of the total kinetic energy of the system at impact. Additionally, these standard techniques do not quantify the contact force, but the impulse over the impact. We present an energy-conserving impact algorithm based on the penalty method to solve for the ground reaction forces during gait. The rigid body assumptions are relaxed and the bodies are allowed to penetrate one another to a small degree. Associated with the deformation is a potential, from which the contact forces are derived. The empirical coefficient-of-restitution used in the standard approaches is replaced by two parameters to characterize the stiffness and the damping of the materials. We solve two simple heel-strike models to illustrate the shortcomings of a standard approach and the suitability of the proposed method for use with gait.
Collapse
Affiliation(s)
- M L Kaplan
- Department of Mechanical Engineering, Division of Biomechanical Engineering, Stanford University, CA 94305-4040, USA
| | | |
Collapse
|