Modeling of inflicted head injury by shaking trauma in children: what can we learn? : Update to parts I&II: A systematic review of animal, mathematical and physical models.

Inflicted shaking trauma can cause injury in infants, but exact injury mechanisms remain unclear. Controversy exists, particularly in courts, whether additional causes such as impact are required to produce injuries found in cases of (suspected) shaking. Publication rates of studies on animal and biomechanical models of inflicted head injury by shaking trauma (IHI-ST) in infants continue rising. Dissention on the topic, combined with its legal relevance, makes maintaining an up-to-date, clear and accessible overview of the current knowledge-base on IHI-ST essential. The current work reviews recent (2017-2023) studies using models of IHI-ST, serving as an update to two previously published reviews. A systematic review was conducted in Scopus and PubMed for articles using animal, physical and mathematical models for IHI-ST. Using the PRISMA methodology, two researchers independently screened the publications. Two, five, and ten publications were included on animal, physical, and mathematical models of IHI-ST, respectively. Both animal model studies used rodents. It is unknown to what degree these can accurately represent IHI-ST. Physical models were used mostly to investigate gross head-kinematics during shaking. Most mathematical models were used to study local effects on the eye and the head's internal structures. All injury thresholds and material properties used were based on scaled adult or animal data. Shaking motions used as inputs for animal, physical and mathematical models were mostly greatly simplified. Future research should focus on using more accurate shaking inputs for models, and on developing or and validating accurate injury thresholds applicable for shaking.

OpenSense: An open-source toolbox for inertial-measurement-unit-based measurement of lower extremity kinematics over long durations.

BACKGROUND: The ability to measure joint kinematics in natural environments over long durations using inertial measurement units (IMUs) could enable at-home monitoring and personalized treatment of neurological and musculoskeletal disorders. However, drift, or the accumulation of error over time, inhibits the accurate measurement of movement over long durations. We sought to develop an open-source workflow to estimate lower extremity joint kinematics from IMU data that was accurate and capable of assessing and mitigating drift. METHODS: We computed IMU-based estimates of kinematics using sensor fusion and an inverse kinematics approach with a constrained biomechanical model. We measured kinematics for 11 subjects as they performed two 10-min trials: walking and a repeated sequence of varied lower-extremity movements. To validate the approach, we compared the joint angles computed with IMU orientations to the joint angles computed from optical motion capture using root mean square (RMS) difference and Pearson correlations, and estimated drift using a linear regression on each subject's RMS differences over time. RESULTS: IMU-based kinematic estimates agreed with optical motion capture; median RMS differences over all subjects and all minutes were between 3 and 6 degrees for all joint angles except hip rotation and correlation coefficients were moderate to strong (r = 0.60-0.87). We observed minimal drift in the RMS differences over 10 min; the average slopes of the linear fits to these data were near zero (- 0.14-0.17 deg/min). CONCLUSIONS: Our workflow produced joint kinematics consistent with those estimated by optical motion capture, and could mitigate kinematic drift even in the trials of continuous walking without rest, which may obviate the need for explicit sensor recalibration (e.g. sitting or standing still for a few seconds or zero-velocity updates) used in current drift-mitigation approaches when studying similar activities. This could enable long-duration measurements, bringing the field one step closer to estimating kinematics in natural environments.

Towards Single Camera Human 3D-Kinematics.

Markerless estimation of 3D Kinematics has the great potential to clinically diagnose and monitor movement disorders without referrals to expensive motion capture labs; however, current approaches are limited by performing multiple de-coupled steps to estimate the kinematics of a person from videos. Most current techniques work in a multi-step approach by first detecting the pose of the body and then fitting a musculoskeletal model to the data for accurate kinematic estimation. Errors in training data of the pose detection algorithms, model scaling, as well the requirement of multiple cameras limit the use of these techniques in a clinical setting. Our goal is to pave the way toward fast, easily applicable and accurate 3D kinematic estimation. To this end, we propose a novel approach for direct 3D human kinematic estimation D3KE from videos using deep neural networks. Our experiments demonstrate that the proposed end-to-end training is robust and outperforms 2D and 3D markerless motion capture based kinematic estimation pipelines in terms of joint angles error by a large margin (35% from 5.44 to 3.54 degrees). We show that D3KE is superior to the multi-step approach and can run at video framerate speeds. This technology shows the potential for clinical analysis from mobile devices in the future.

OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement.

Movement is fundamental to human and animal life, emerging through interaction of complex neural, muscular, and skeletal systems. Study of movement draws from and contributes to diverse fields, including biology, neuroscience, mechanics, and robotics. OpenSim unites methods from these fields to create fast and accurate simulations of movement, enabling two fundamental tasks. First, the software can calculate variables that are difficult to measure experimentally, such as the forces generated by muscles and the stretch and recoil of tendons during movement. Second, OpenSim can predict novel movements from models of motor control, such as kinematic adaptations of human gait during loaded or inclined walking. Changes in musculoskeletal dynamics following surgery or due to human-device interaction can also be simulated; these simulations have played a vital role in several applications, including the design of implantable mechanical devices to improve human grasping in individuals with paralysis. OpenSim is an extensible and user-friendly software package built on decades of knowledge about computational modeling and simulation of biomechanical systems. OpenSim's design enables computational scientists to create new state-of-the-art software tools and empowers others to use these tools in research and clinical applications. OpenSim supports a large and growing community of biomechanics and rehabilitation researchers, facilitating exchange of models and simulations for reproducing and extending discoveries. Examples, tutorials, documentation, and an active user forum support this community. The OpenSim software is covered by the Apache License 2.0, which permits its use for any purpose including both nonprofit and commercial applications. The source code is freely and anonymously accessible on GitHub, where the community is welcomed to make contributions. Platform-specific installers of OpenSim include a GUI and are available on simtk.org.

Is my model good enough? Best practices for verification and validation of musculoskeletal models and simulations of movement.

Computational modeling and simulation of neuromusculoskeletal (NMS) systems enables researchers and clinicians to study the complex dynamics underlying human and animal movement. NMS models use equations derived from physical laws and biology to help solve challenging real-world problems, from designing prosthetics that maximize running speed to developing exoskeletal devices that enable walking after a stroke. NMS modeling and simulation has proliferated in the biomechanics research community over the past 25 years, but the lack of verification and validation standards remains a major barrier to wider adoption and impact. The goal of this paper is to establish practical guidelines for verification and validation of NMS models and simulations that researchers, clinicians, reviewers, and others can adopt to evaluate the accuracy and credibility of modeling studies. In particular, we review a general process for verification and validation applied to NMS models and simulations, including careful formulation of a research question and methods, traditional verification and validation steps, and documentation and sharing of results for use and testing by other researchers. Modeling the NMS system and simulating its motion involves methods to represent neural control, musculoskeletal geometry, muscle-tendon dynamics, contact forces, and multibody dynamics. For each of these components, we review modeling choices and software verification guidelines; discuss variability, errors, uncertainty, and sensitivity relationships; and provide recommendations for verification and validation by comparing experimental data and testing robustness. We present a series of case studies to illustrate key principles. In closing, we discuss challenges the community must overcome to ensure that modeling and simulation are successfully used to solve the broad spectrum of problems that limit human mobility.

Minimization of metabolic cost of transport predicts changes in gait mechanics over a range of ankle-foot orthosis stiffnesses in individuals with bilateral plantar flexor weakness.

Neuromuscular disorders often lead to ankle plantar flexor muscle weakness, which impairs ankle push-off power and forward propulsion during gait. To improve walking speed and reduce metabolic cost of transport (mCoT), patients with plantar flexor weakness are provided dorsal-leaf spring ankle-foot orthoses (AFOs). It is widely believed that mCoT during gait depends on the AFO stiffness and an optimal AFO stiffness that minimizes mCoT exists. The biomechanics behind why and how an optimal stiffness exists and benefits individuals with plantar flexor weakness are not well understood. We hypothesized that the AFO would reduce the required support moment and, hence, metabolic cost contributions of the ankle plantar flexor and knee extensor muscles during stance, and reduce hip flexor metabolic cost to initiate swing. To test these hypotheses, we generated neuromusculoskeletal simulations to represent gait of an individual with bilateral plantar flexor weakness wearing an AFO with varying stiffness. Predictions were based on the objective of minimizing mCoT, loading rates at impact and head accelerations at each stiffness level, and the motor patterns were determined via dynamic optimization. The predictive gait simulation results were compared to experimental data from subjects with bilateral plantar flexor weakness walking with varying AFO stiffness. Our simulations demonstrated that reductions in mCoT with increasing stiffness were attributed to reductions in quadriceps metabolic cost during midstance. Increases in mCoT above optimum stiffness were attributed to the increasing metabolic cost of both hip flexor and hamstrings muscles. The insights gained from our predictive gait simulations could inform clinicians on the prescription of personalized AFOs. With further model individualization, simulations based on mCoT minimization may sufficiently predict adaptations to an AFO in individuals with plantar flexor weakness.

How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds.

The lengths and velocities of muscle fibers have a dramatic effect on muscle force generation. It is unknown, however, whether the lengths and velocities of lower limb muscle fibers substantially affect the ability of muscles to generate force during walking and running. We examined this issue by developing simulations of muscle-tendon dynamics to calculate the lengths and velocities of muscle fibers from electromyographic recordings of 11 lower limb muscles and kinematic measurements of the hip, knee and ankle made as five subjects walked at speeds of 1.0-1.75 m s(-1) and ran at speeds of 2.0-5.0 m s(-1). We analyzed the simulated fiber lengths, fiber velocities and forces to evaluate the influence of force-length and force-velocity properties on force generation at different walking and running speeds. The simulations revealed that force generation ability (i.e. the force generated per unit of activation) of eight of the 11 muscles was significantly affected by walking or running speed. Soleus force generation ability decreased with increasing walking speed, but the transition from walking to running increased the force generation ability by reducing fiber velocities. Our results demonstrate the influence of soleus muscle architecture on the walk-to-run transition and the effects of muscle-tendon compliance on the plantarflexors' ability to generate ankle moment and power. The study presents data that permit lower limb muscles to be studied in unprecedented detail by relating muscle fiber dynamics and force generation to the mechanical demands of walking and running.

Flexing computational muscle: modeling and simulation of musculotendon dynamics.

Muscle-driven simulations of human and animal motion are widely used to complement physical experiments for studying movement dynamics. Musculotendon models are an essential component of muscle-driven simulations, yet neither the computational speed nor the biological accuracy of the simulated forces has been adequately evaluated. Here we compare the speed and accuracy of three musculotendon models: two with an elastic tendon (an equilibrium model and a damped equilibrium model) and one with a rigid tendon. Our simulation benchmarks demonstrate that the equilibrium and damped equilibrium models produce similar force profiles but have different computational speeds. At low activation, the damped equilibrium model is 29 times faster than the equilibrium model when using an explicit integrator and 3 times faster when using an implicit integrator; at high activation, the two models have similar simulation speeds. In the special case of simulating a muscle with a short tendon, the rigid-tendon model produces forces that match those generated by the elastic-tendon models, but simulates 2-54 times faster when an explicit integrator is used and 6-31 times faster when an implicit integrator is used. The equilibrium, damped equilibrium, and rigid-tendon models reproduce forces generated by maximally-activated biological muscle with mean absolute errors less than 8.9%, 8.9%, and 20.9% of the maximum isometric muscle force, respectively. When compared to forces generated by submaximally-activated biological muscle, the forces produced by the equilibrium, damped equilibrium, and rigid-tendon models have mean absolute errors less than 16.2%, 16.4%, and 18.5%, respectively. To encourage further development of musculotendon models, we provide implementations of each of these models in OpenSim version 3.1 and benchmark data online, enabling others to reproduce our results and test their models of musculotendon dynamics.

Does enforcing glenohumeral joint stability matter? A new rapid muscle redundancy solver highlights the importance of non-superficial shoulder muscles.

The complexity of the human shoulder girdle enables the large mobility of the upper extremity, but also introduces instability of the glenohumeral (GH) joint. Shoulder movements are generated by coordinating large superficial and deeper stabilizing muscles spanning numerous degrees-of-freedom. How shoulder muscles are coordinated to stabilize the movement of the GH joint remains widely unknown. Musculoskeletal simulations are powerful tools to gain insights into the actions of individual muscles and particularly of those that are difficult to measure. In this study, we analyze how enforcement of GH joint stability in a musculoskeletal model affects the estimates of individual muscle activity during shoulder movements. To estimate both muscle activity and GH stability from recorded shoulder movements, we developed a Rapid Muscle Redundancy (RMR) solver to include constraints on joint reaction forces (JRFs) from a musculoskeletal model. The RMR solver yields muscle activations and joint forces by minimizing the weighted sum of squared-activations, while matching experimental motion. We implemented three new features: first, computed muscle forces include active and passive fiber contributions; second, muscle activation rates are enforced to be physiological, and third, JRFs are efficiently formulated as linear functions of activations. Muscle activity from the RMR solver without GH stability was not different from the computed muscle control (CMC) algorithm and electromyography of superficial muscles. The efficiency of the solver enabled us to test over 3600 trials sampled within the uncertainty of the experimental movements to test the differences in muscle activity with and without GH joint stability enforced. We found that enforcing GH stability significantly increases the estimated activity of the rotator cuff muscles but not of most superficial muscles. Therefore, a comparison of shoulder model muscle activity to EMG measurements of superficial muscles alone is insufficient to validate the activity of rotator cuff muscles estimated from musculoskeletal models.

Personalisation of Plantarflexor Musculotendon Model Parameters in Children with Cerebral Palsy.

Neuromusculoskeletal models can be used to evaluate aberrant muscle function in cerebral palsy (CP), for example by estimating muscle and joint contact forces during gait. However, to be accurate, models should include representative musculotendon parameters. We aimed to estimate personalised parameters that capture the mechanical behaviour of the plantarflexors in children with CP and typically developing (TD) children. Ankle angle (using motion capture), torque (using a load-cell), and medial gastrocnemius fascicle lengths (using ultrasound) were measured during slow passive ankle dorsiflexion rotation for thirteen children with spastic CP and thirteen TD children. Per subject, the measured rotation was input to a scaled OpenSim model to simulate the torque and fascicle length output. Musculotendon model parameters were personalised by the best match between simulated and experimental torque-angle and fascicle length-angle curves according to a least-squares fit. Personalised tendon slack lengths were significantly longer and optimal fibre lengths significantly shorter in CP than model defaults and than in TD. Personalised tendon compliance was substantially higher in both groups compared to the model default. The presented method to personalise musculotendon parameters will likely yield more accurate simulations of subject-specific muscle mechanics, to help us understand the effects of altered musculotendon properties in CP.

Conclusion or Illusion: Quantifying Uncertainty in Inverse Analyses From Marker-Based Motion Capture due to Errors in Marker Registration and Model Scaling.

Estimating kinematics from optical motion capture with skin-mounted markers, referred to as an inverse kinematic (IK) calculation, is the most common experimental technique in human motion analysis. Kinematics are often used to diagnose movement disorders and plan treatment strategies. In many such applications, small differences in joint angles can be clinically significant. Kinematics are also used to estimate joint powers, muscle forces, and other quantities of interest that cannot typically be measured directly. Thus, the accuracy and reproducibility of IK calculations are critical. In this work, we isolate and quantify the uncertainty in joint angles, moments, and powers due to two sources of error during IK analyses: errors in the placement of markers on the model (marker registration) and errors in the dimensions of the model's body segments (model scaling). We demonstrate that IK solutions are best presented as a distribution of equally probable trajectories when these sources of modeling uncertainty are considered. Notably, a substantial amount of uncertainty exists in the computed kinematics and kinetics even if low marker tracking errors are achieved. For example, considering only 2 cm of marker registration uncertainty, peak ankle plantarflexion angle varied by 15.9°, peak ankle plantarflexion moment varied by 26.6 N⋅m, and peak ankle power at push off varied by 75.9 W during healthy gait. This uncertainty can directly impact the classification of patient movements and the evaluation of training or device effectiveness, such as calculations of push-off power. We provide scripts in OpenSim so that others can reproduce our results and quantify the effect of modeling uncertainty in their own studies.

Muscle coordination retraining inspired by musculoskeletal simulations reduces knee contact force.

Humans typically coordinate their muscles to meet movement objectives like minimizing energy expenditure. In the presence of pathology, new objectives gain importance, like reducing loading in an osteoarthritic joint, but people often do not change their muscle coordination patterns to meet these new objectives. Here we use musculoskeletal simulations to identify simple changes in coordination that can be taught using electromyographic biofeedback, achieving the therapeutic goal of reducing joint loading. Our simulations predicted that changing the relative activation of two redundant ankle plantarflexor muscles-the gastrocnemius and soleus-could reduce knee contact force during walking, but it was unclear whether humans could re-coordinate redundant muscles during a complex task like walking. Our experiments showed that after a single session of walking with biofeedback of summary measures of plantarflexor muscle activation, healthy individuals reduced the ratio of gastrocnemius-to-soleus muscle activation by 25 ± 15% (p = 0.004, paired t test, n = 10). Participants who walked with this "gastrocnemius avoidance" gait pattern reduced late-stance knee contact force by 12 ± 12% (p = 0.029, paired t test, n = 8). Simulation-informed coordination retraining could be a promising treatment for knee osteoarthritis and a powerful tool for optimizing coordination for a variety of rehabilitation and performance applications.

A marker registration method to improve joint angles computed by constrained inverse kinematics.

Accurate computation of joint angles from optical marker data using inverse kinematics methods requires that the locations of markers on a model match the locations of experimental markers on participants. Marker registration is the process of positioning the model markers so that they match the locations of the experimental markers. Markers are typically registered using a graphical user interface (GUI), but this method is subjective and may introduce errors and uncertainty to the calculated joint angles and moments. In this investigation, we use OpenSim to isolate and quantify marker registration-based error from other sources of error by analyzing the gait of a bipedal humanoid robot for which segment geometry, mass properties, and joint angles are known. We then propose a marker registration method that is informed by the orientation of anatomical reference frames derived from surface-mounted optical markers as an alternative to user registration using a GUI. The proposed orientation registration method reduced the average root-mean-square error in both joint angles and joint moments by 67% compared to the user registration method, and eliminated variability among users. Our results show that a systematic method for marker registration that reduces subjective user input can make marker registration more accurate and repeatable.

Correction: A marker registration method to improve joint angles computed by constrained inverse kinematics.

[This corrects the article DOI: 10.1371/journal.pone.0252425.].

Minimal formulation of joint motion for biomechanisms.

Biomechanical systems share many properties with mechanically engineered systems, and researchers have successfully employed mechanical engineering simulation software to investigate the mechanical behavior of diverse biological mechanisms, ranging from biomolecules to human joints. Unlike their man-made counterparts, however, biomechanisms rarely exhibit the simple, uncoupled, pure-axial motion that is engineered into mechanical joints such as sliders, pins, and ball-and-socket joints. Current mechanical modeling software based on internal-coordinate multibody dynamics can formulate engineered joints directly in minimal coordinates, but requires additional coordinates restricted by constraints to model more complex motions. This approach can be inefficient, inaccurate, and difficult for biomechanists to customize. Since complex motion is the rule rather than the exception in biomechanisms, the benefits of minimal coordinate modeling are not fully realized in biomedical research. Here we introduce a practical implementation for empirically-defined internal-coordinate joints, which we call "mobilizers." A mobilizer encapsulates the observations, measurement frame, and modeling requirements into a hinge specification of the permissible-motion manifold for a minimal set of internal coordinates. Mobilizers support nonlinear mappings that are mathematically equivalent to constraint manifolds but have the advantages of fewer coordinates, no constraints, and exact representation of the biomechanical motion-space-the benefits long enjoyed for internal-coordinate models of mechanical joints. Hinge matrices within the mobilizer are easily specified by user-supplied functions, and provide a direct means of mapping permissible motion derived from empirical data. We present computational results showing substantial performance and accuracy gains for mobilizers versus equivalent joints implemented with constraints. Examples of mobilizers for joints from human biomechanics and molecular dynamics are given. All methods and examples were implemented in Simbody™-an open source multibody-dynamics solver available at https://Simtk.org.

Muscle Contributions to Upper-Extremity Movement and Work From a Musculoskeletal Model of the Human Shoulder.

Musculoskeletal models enable movement scientists to examine muscle function by computing the mechanical work done by muscles during motor tasks. To estimate muscle work accurately requires a model that is physiologically plausible. Previous models of the human shoulder have coupled scapula movement to humeral movement. While coupled movement produces a stereotypical scapulohumeral rhythm, it cannot model shrugging or independent movement of the scapula and humerus. The artificial coupling of humeral elevation to scapular rotation permits muscles that cross the glenohumeral joint, such as the rotator-cuff muscles and deltoids, to do implausible work to elevate and rotate the scapula. In reality, the motion of the scapula is controlled by thoracoscapular muscles, yet the roles of these muscles in shoulder function remains unclear. To elucidate the roles of the thoracoscapular muscles, we developed a shoulder model with an accurate scapulothoracic joint and includes scapular muscles to drive its motion. We used the model to compute the work done by the thoracoscapular muscles during shrugging and arm elevation. We found that the bulk of the work done in upper-extremity tasks is performed by the largest muscles of the shoulder: trapezius, deltoids, pectoralis major, and serratus-anterior. Trapezius and serratus anterior prove to be important synergists in performing upward-rotation of the scapula. We show that the large thoracoscapular muscles do more work than glenohumeral muscles during arm-elevation tasks. The model, experimental data and simulation results are freely available on SimTK.org to enable anyone to explore our results and to perform further studies in OpenSim 4.0.

A neuromusculoskeletal tracking method for estimating individual muscle forces in human movement.

A neuromusculoskeletal tracking (NMT) method was developed to estimate muscle forces from observed motion data. The NMT method combines skeletal motion tracking and optimal neuromuscular tracking to produce forward simulations of human movement quickly and accurately. The skeletal motion tracker calculates the joint torques needed to actuate a skeletal model and track observed segment angles and ground forces in a forward simulation of the motor task. The optimal neuromuscular tracker resolves the muscle redundancy problem dynamically and finds the muscle excitations (and muscle forces) needed to produce the joint torques calculated by the skeletal motion tracker. To evaluate the accuracy of the NMT method, kinematics and ground forces obtained from an optimal control (parameter optimization) solution for maximum-height jumping were contaminated with both random and systematic noise. These data served as input observations to the NMT method as well as an inverse dynamics analysis. The NMT solution was compared to the input observations, the original optimal solution, and a simulation driven by the inverse dynamics torques. The results show that, in contrast to inverse dynamics, the NMT method is able to produce an accurate forward simulation consistent with the optimal control solution. The NMT method also requires 3 orders-of-magnitude less CPU time than parameter optimization. The speed and accuracy of the NMT method make it a promising new tool for estimating muscle forces using experimentally obtained kinematics and ground force data.

A Biomechanical Model of the Scapulothoracic Joint to Accurately Capture Scapular Kinematics during Shoulder Movements.

The complexity of shoulder mechanics combined with the movement of skin relative to the scapula makes it difficult to measure shoulder kinematics with sufficient accuracy to distinguish between symptomatic and asymptomatic individuals. Multibody skeletal models can improve motion capture accuracy by reducing the space of possible joint movements, and models are used widely to improve measurement of lower limb kinematics. In this study, we developed a rigid-body model of a scapulothoracic joint to describe the kinematics of the scapula relative to the thorax. This model describes scapular kinematics with four degrees of freedom: 1) elevation and 2) abduction of the scapula on an ellipsoidal thoracic surface, 3) upward rotation of the scapula normal to the thoracic surface, and 4) internal rotation of the scapula to lift the medial border of the scapula off the surface of the thorax. The surface dimensions and joint axes can be customized to match an individual's anthropometry. We compared the model to "gold standard" bone-pin kinematics collected during three shoulder tasks and found modeled scapular kinematics to be accurate to within 2 mm root-mean-squared error for individual bone-pin markers across all markers and movement tasks. As an additional test, we added random and systematic noise to the bone-pin marker data and found that the model reduced kinematic variability due to noise by 65% compared to Euler angles computed without the model. Our scapulothoracic joint model can be used for inverse and forward dynamics analyses and to compute joint reaction loads. The computational performance of the scapulothoracic joint model is well suited for real-time applications; it is freely available for use with OpenSim 3.2, and is customizable and usable with other OpenSim models.

Rectus femoris transfer surgery affects balance recovery in children with cerebral palsy: A computer simulation study.

Stiff-knee gait is a troublesome movement disorder among children with cerebral palsy (CP), where peak swing phase knee flexion is diminished due to over-activity of the rectus femoris muscle. A common treatment for stiff-knee gait, rectus femoris transfer surgery, moves the muscle's distal tendon from the patella to the sartorius insertion on the tibia. As a biarticular muscle, rectus femoris may play a role in motor control and have unrecognized benefits for maintaining balance. We used musculoskeletal modeling, neuromuscular control, and forward dynamic simulation to investigate the role of rectus femoris tendon transfer surgery on balance recovery after support-surface perturbations for children with CP adopting two different crouched postures. We combined both high-level supraspinal and low-level spinal signals to generate 92 muscle excitations for tracking experimental whole body center of mass positions and velocities. Stability during balance recovery was evaluated by the minimum distance between the extrapolated center of mass and base of support boundary (bmin) and the minimum time to reach the boundary (TtBmin). The balance recovery of pre-surgical simulations (bmin=2.3+1.1cm, TtBmin=0.2+0.1s) were different (p=0.02), on average, than post-surgical simulations (bmin=-4.9+11.4cm, TtBmin=-0.1+0.3s) of rectus femoris transfers. The moderate crouch simulations (bmin=2.4+0.4cm, TtBmin=0.2+0.03s) were more stable than the mild crouch simulations (bmin=1.2+0.3cm, TtBmin=0.1+0.02s) following anterior translations of the support surface. These findings suggest that tendon transfer of rectus femoris affects balance recovery in children with CP.

Simulating Ideal Assistive Devices to Reduce the Metabolic Cost of Running.

Tools have been used for millions of years to augment the capabilities of the human body, allowing us to accomplish tasks that would otherwise be difficult or impossible. Powered exoskeletons and other assistive devices are sophisticated modern tools that have restored bipedal locomotion in individuals with paraplegia and have endowed unimpaired individuals with superhuman strength. Despite these successes, designing assistive devices that reduce energy consumption during running remains a substantial challenge, in part because these devices disrupt the dynamics of a complex, finely tuned biological system. Furthermore, designers have hitherto relied primarily on experiments, which cannot report muscle-level energy consumption and are fraught with practical challenges. In this study, we use OpenSim to generate muscle-driven simulations of 10 human subjects running at 2 and 5 m/s. We then add ideal, massless assistive devices to our simulations and examine the predicted changes in muscle recruitment patterns and metabolic power consumption. Our simulations suggest that an assistive device should not necessarily apply the net joint moment generated by muscles during unassisted running, and an assistive device can reduce the activity of muscles that do not cross the assisted joint. Our results corroborate and suggest biomechanical explanations for similar effects observed by experimentalists, and can be used to form hypotheses for future experimental studies. The models, simulations, and software used in this study are freely available at simtk.org and can provide insight into assistive device design that complements experimental approaches.