Joint contact forces during semi-recumbent seated cycling.

Semi-recumbent cycling performed from a wheelchair is a popular rehabilitation exercise following spinal cord injury (SCI) and is often paired with functional electrical stimulation. However, biomechanical assessment of this cycling modality is lacking, even in unimpaired populations, hindering the development of personalised and safe rehabilitation programs for those with SCI. This study developed a computational pipeline to determine lower limb kinematics, kinetics, and joint contact forces (JCF) in 11 unimpaired participants during voluntary semi-recumbent cycling using a rehabilitation ergometer. Two cadences (40 and 60 revolutions per minute) and three crank powers (15 W, 30 W, and 45 W) were assessed. A rigid body model of a rehabilitation ergometer was combined with a calibrated electromyogram-informed neuromusculoskeletal model to determine JCF at the hip, knee, and ankle. Joint excursions remained consistent across all cadence and powers, but joint moments and JCF differed between 40 and 60 revolutions per minute, with peak JCF force significantly greater at 40 compared to 60 revolutions per minute for all crank powers. Poor correlations were found between mean crank power and peak JCF across all joints. This study provides foundation data and computational methods to enable further evaluation and optimisation of semi-recumbent cycling for application in rehabilitation after SCI and other neurological disorders.

Muscle synergy-informed neuromusculoskeletal modelling to estimate knee contact forces in children with cerebral palsy.

Cerebral palsy (CP) includes a group of neurological conditions caused by damage to the developing brain, resulting in maladaptive alterations of muscle coordination and movement. Estimates of joint moments and contact forces during locomotion are important to establish the trajectory of disease progression and plan appropriate surgical interventions in children with CP. Joint moments and contact forces can be estimated using electromyogram (EMG)-informed neuromusculoskeletal models, but a reduced number of EMG sensors would facilitate translation of these computational methods to clinics. This study developed and evaluated a muscle synergy-informed neuromusculoskeletal modelling approach using EMG recordings from three to four muscles to estimate joint moments and knee contact forces of children with CP and typically developing (TD) children during walking. Using only three to four experimental EMG sensors attached to a single leg and leveraging an EMG database of walking data of TD children, the synergy-informed approach estimated total knee contact forces comparable to those estimated by EMG-assisted approaches that used 13 EMG sensors (children with CP, n = 3, R2 = 0.95 ± 0.01, RMSE = 0.40 ± 0.14 BW; TD controls, n = 3, R2 = 0.93 ± 0.07, RMSE = 0.19 ± 0.05 BW). The proposed synergy-informed neuromusculoskeletal modelling approach could enable rapid evaluation of joint biomechanics in children with unimpaired and impaired motor control within a clinical environment.

Prediction of Achilles Tendon Force During Common Motor Tasks From Markerless Video.

Remodeling of the Achilles tendon (AT) is partly driven by its mechanical environment. AT force can be estimated with neuromusculoskeletal (NMSK) modeling; however, the complex experimental setup required to perform the analyses confines use to the laboratory. We developed task-specific long short-term memory (LSTM) neural networks that employ markerless video data to predict the AT force during walking, running, countermovement jump, single-leg landing, and single-leg heel rise. The task-specific LSTM models were trained on pose estimation keypoints and corresponding AT force data from 16 subjects, calculated via an established NMSK modeling pipeline, and cross-validated using a leave-one-subject-out approach. As proof-of-concept, new motion data of one participant was collected with two smartphones and used to predict AT forces. The task-specific LSTM models predicted the time-series AT force using synthesized pose estimation data with root mean square error (RMSE) ≤ 526 N, normalized RMSE (nRMSE) ≤ 0.21 , R 2 ≥ 0.81 . Walking task resulted the most accurate with RMSE = 189±62 N; nRMSE = 0.11±0.03 , R 2 = 0.92±0.04 . AT force predicted with smartphones video data was physiologically plausible, agreeing in timing and magnitude with established force profiles. This study demonstrated the feasibility of using low-cost solutions to deploy complex biomechanical analyses outside the laboratory.

Sagittal plane knee kinematics can be measured during activities of daily living following total knee arthroplasty with two IMU.

Knee function is rarely measured objectively during functional tasks following total knee arthroplasty. Inertial measurement units (IMU) can measure knee kinematics and range of motion (ROM) during dynamic activities and offer an easy-to-use system for knee function assessment post total knee arthroplasty. However, IMU must be validated against gold standard three-dimensional optical motion capture systems (OMC) across a range of tasks if they are to see widespread uptake. We computed knee rotations and ROM from commercial IMU sensor measurements during walking, squatting, sit-to-stand, stair ascent, and stair descent in 21 patients one-year post total knee arthroplasty using two methods: direct computation using segment orientations (r_IMU), and an IMU-driven iCloud-based interactive lower limb model (m_IMU). This cross-sectional study compared computed knee angles and ROM to a gold-standard OMC and inverse kinematics method using Pearson's correlation coefficient (R) and root-mean-square-differences (RMSD). The r_IMU and m_IMU methods estimated sagittal plane knee angles with excellent correlation (>0.95) compared to OMC for walking, squatting, sit-to-stand, and stair-ascent, and very good correlation (>0.90) for stair descent. For squatting, sit-to-stand, and walking, the mean RMSD for r_IMU and m_IMU compared to OMC were <4 degrees, < 5 degrees, and <6 degrees, respectively but higher for stair ascent and descent (~12 degrees). Frontal and transverse plane knee kinematics estimated using r_IMU and m_IMU showed poor to moderate correlation compared to OMC. There were no differences in ROM measurements during squatting, sit-to-stand, and walking across the two methods. Thus, IMUs can measure sagittal plane knee angles and ROM with high accuracy for a variety of tasks and may be a useful in-clinic tool for objective assessment of knee function following total knee arthroplasty.

Inclusion of a skeletal model partly improves the reliability of lower limb joint angles derived from a markerless depth camera.

A single depth camera provides a fast and easy approach to performing biomechanical assessments in a clinical setting; however, there are currently no established methods to reliably determine joint angles from these devices. The primary aim of this study was to compare joint angles as well as the between-day reliability of direct kinematics to model-constrained inverse kinematics recorded using a single markerless depth camera during a range of clinical and athletic movement assessments.A secondary aim was to determine the minimum number of trials required to maximize reliability. Eighteen healthy participants attended two testing sessions one week apart. Tasks included treadmill walking, treadmill running, single-leg squats, single-leg countermovement jumps, bilateral countermovement jumps, and drop vertical jumps. Keypoint data were processed using direct kinematics as well as in OpenSim using a full-body musculoskeletal model and inverse kinematics. Kinematic methods were compared using statistical parametric mapping and between-day reliability was calculated using intraclass correlation coefficients, mean absolute error, and minimal detectable change. Keypoint-derived inverse kinematics resulted in significantly smaller hip flexion (range = -9 to -2°), hip abduction (range = -3 to -2°), knee flexion (range = -5° to -2°), and greater dorsiflexion angles (range = 6-15°) than direct kinematics. Both markerless kinematic methods had high between-day reliability (inverse kinematics ICC 95 %CI = 0.83-0.90; direct kinematics ICC 95 %CI = 0.80-0.93). For certain tasks and joints, keypoint-derived inverse kinematics resulted in greater reliability (up to 0.47 ICC) and smaller minimal detectable changes (up to 13°) than direct kinematics. Performing 2-4 trials was sufficient to maximize reliability for most tasks. A single markerless depth camera can reliably measure lower limb joint angles, and skeletal model-constrained inverse kinematics improves lower limb joint angle reliability for certain tasks and joints.

Effect of different constraining boundary conditions on simulated femoral stresses and strains during gait.

Finite element analysis (FEA) is commonly used in orthopaedic research to estimate localised tissue stresses and strains. A variety of boundary conditions have been proposed for isolated femur analysis, but it remains unclear how these assumed constraints influence FEA predictions of bone biomechanics. This study compared the femoral head deflection (FHD), stresses, and strains elicited under four commonly used boundary conditions (fixed knee, mid-shaft constraint, springs, and isostatic methods) and benchmarked these mechanics against the gold standard inertia relief method for normal and pathological femurs (extreme anteversion and retroversion, coxa vara, and coxa valga). Simulations were performed for the stance phase of walking with the applied femoral loading determined from patient-specific neuromusculoskeletal models. Due to unrealistic biomechanics observed for the commonly used boundary conditions, we propose a novel biomechanical constraint method to generate physiological femur biomechanics. The biomechanical method yielded FHD (< 1 mm), strains (approaching 1000 µÎµ), and stresses (< 60 MPa), which were consistent with physiological observations and similar to predictions from the inertia relief method (average coefficient of determination = 0.97, average normalized root mean square error = 0.17). Our results highlight the superior performance of the biomechanical method compared to current methods of constraint for  both healthy and pathological femurs.