Femoral bone growth predictions based on personalized multi-scale simulations: validation and sensitivity analysis of a mechanobiological model

Musculoskeletal function is pivotal to long-term health. However, various patient groups develop torsional deformities, leading to clinical, functional problems. Understanding the interplay between movement pattern, bone loading and growth is crucial for improving the functional mobility of these patients and preserving long-term health. Multi-scale simulations in combination with a mechanobiological bone growth model have been used to estimate bone loads and predict femoral growth trends based on cross-sectional data. The lack of longitudinal data in the previous studies hindered refinements of the mechanobiological model and validation of subject-specific growth predictions, thereby limiting clinical applications. This study aimed to validate the growth predictions using magnetic resonance images and motion capture data-collected longitudinally-from ten growing children. Additionally, a sensitivity analysis was conducted to refine model parameters. A linear regression model based on physical activity information, anthropometric data and predictions from the refined mechanobiological model explained 70% of femoral anteversion development. Notably, the direction of femoral development was accurately predicted in 18 out of 20 femurs, suggesting that growth predictions could help to revolutionize treatment strategies for torsional deformities.

Facial expression recognition through muscle synergies and estimation of facial keypoint displacements through a skin-musculoskeletal model using facial sEMG signals

The development of facial expression recognition (FER) and facial expression generation (FEG) systems is essential to enhance human-robot interactions (HRI). The facial action coding system is widely used in FER and FEG tasks, as it offers a framework to relate the action of facial muscles and the resulting facial motions to the execution of facial expressions. However, most FER and FEG studies are based on measuring and analyzing facial motions, leaving the facial muscle component relatively unexplored. This study introduces a novel framework using surface electromyography (sEMG) signals from facial muscles to recognize facial expressions and estimate the displacement of facial keypoints during the execution of the expressions. For the facial expression recognition task, we studied the coordination patterns of seven muscles, expressed as three muscle synergies extracted through non-negative matrix factorization, during the execution of six basic facial expressions. Muscle synergies are groups of muscles that show coordinated patterns of activity, as measured by their sEMG signals, and are hypothesized to form the building blocks of human motor control. We then trained two classifiers for the facial expressions based on extracted features from the sEMG signals and the synergy activation coefficients of the extracted muscle synergies, respectively. The accuracy of both classifiers outperformed other systems that use sEMG to classify facial expressions, although the synergy-based classifier performed marginally worse than the sEMG-based one (classification accuracy: synergy-based 97.4%, sEMG-based 99.2%). However, the extracted muscle synergies revealed common coordination patterns between different facial expressions, allowing a low-dimensional quantitative visualization of the muscle control strategies involved in human facial expression generation. We also developed a skin-musculoskeletal model enhanced by linear regression (SMSM-LRM) to estimate the displacement of facial keypoints during the execution of a facial expression based on sEMG signals. Our proposed approach achieved a relatively high fidelity in estimating these displacements (NRMSE 0.067). We propose that the identified muscle synergies could be used in combination with the SMSM-LRM model to generate motor commands and trajectories for desired facial displacements, potentially enabling the generation of more natural facial expressions in social robotics and virtual reality.

Knowledge-Based Deep Learning for Time-Efficient Inverse Dynamics

Accurate understanding of muscle activation and muscle forces plays an essential role in neuro-rehabilitation and musculoskeletal disorder treatments. Computational musculoskeletal modeling has been widely used as a powerful non-invasive tool to estimate them through inverse dynamics using static optimization, but the inherent computational complexity results in time-consuming analysis. In this paper, we propose a knowledge-based deep learning framework for time-efficient inverse dynamic analysis, which can predict muscle activation and muscle forces from joint kinematic data directly while not requiring any label information during model training. The Bidirectional Gated Recurrent Unit (BiGRU) neural network is selected as the backbone of our model due to its proficient handling of time-series data. Prior physical knowledge from forward dynamics and pre-selected inverse dynamics based physiological criteria are integrated into the loss function to guide the training of neural networks. Experimental validations on two datasets, including one benchmark upper limb movement dataset and one self-collected lower limb movement dataset from six healthy subjects, are performed. The experimental results have shown that the selected BiGRU architecture outperforms other neural network models when trained using our specifically designed loss function, which illustrates the effectiveness and robustness of the proposed framework.

Medial and lateral knee contact forces and muscle forces during sit-to-stand in patients one year after unilateral total knee arthroplasty

Understanding how forces are transmitted through the knee after TKA is essential, as it may explain why many patients experience pain or functional limitations during various activities. This study compared knee muscle forces and knee contact forces (KCF) during sit-to-stand in patients one year after unilateral total knee arthroplasty (TKA) with either a medial ball-and-socket (MBS) or posterior stabilized (PS) implant and compared them to a group of similarly healthy aged controls (CTRL). A musculoskeletal model and static optimization estimated lower limb kinematics, knee kinetics, muscle forces, and KCFs. The normalized sit-to-stand cycle was compared among the groups using statistical nonparametric mapping, and peak between-limb differences were compared using discrete statistics. The PS group required greater forward lean during the sit-to-stand task, causing greater spine flexion, posterior pelvic tilt, and decreased hip flexion on the operated limb. PS and MBS groups favoured their non-operated limb, resulting in less range of motion throughout the lower limb, lower knee kinetics, muscle forces, and KCFs on the operated limb. Compared to the controls, the MBS and PS groups had reduced medial compartment KCF. The control group did favour their dominant limb over their non-dominant limb. Post-operative rehabilitation should continue to promote greater use of the operated knee to have more symmetrical loading between operated and non-operated limbs and improve strength and mobility at the hip and ankle joints. One year after surgery, TKA patients remain with reduced muscle forces and KCF on their operated limb during a sit-to-stand task, regardless whether they received an MBS or PS implant.

Validation of full-length radiograph based musculoskeletal modeling method to estimate medial and lateral knee contact forces

BACKGROUND

Anatomical parameters of the pelvis, femur, and tibia derived from the full-length radiograph can be used to create a more accurate musculoskeletal model compared to marker-based linear scaling method. However, whether this model leads to more accurate estimations of medial knee contact force (MCF) and lateral knee contact force (LCF) than marker-based linear scaling method is still unknown.

RESEARCH QUESTION

This main purpose of this study was to determine whether musculoskeletal model generated from full-length radiograph improves the estimations of MCF and LCF.

METHODS

An open-source dataset including marker trajectories, ground reaction forces, in vivo knee contact forces, and full-length radiograph was used to evaluate the accuracy of full-length radiograph musculoskeletal modeling method. Subject-specific musculoskeletal models were created using anatomical parameters derived from the full-length radiograph or marker-based linear scaling methods. MCF and LCF were estimated using musculoskeletal simulations of normal walking trails. The accuracy of modeling methods was determined by comparing the estimated and in vivo measured MCF and LCF.

RESULTS

Compared to the marker-based linear scaling approach, the full-length radiograph musculoskeletal modeling method exhibited decreases of 38.3 % and 41.3 % in root mean square error for MCF and LCF respectively, as well as reductions of 50.0 % and 49.3 % in mean peak errors for MCF and LCF respectively.

SIGNIFICANCE

The full-length radiograph musculoskeletal modeling method provides a more accurate way to estimate MCF and LCF compared to the traditional maker-based linear scaling approach, which may contribute to understand the initiation, progression, and treatment of OA.

Myobolica: A Stochastic Approach to Estimate Physiological Muscle Control Variability

The inherent redundancy of the musculoskeletal systems is traditionally solved by optimizing a cost function. This approach may not be correct to model non-adult or pathological populations likely to adopt a "non-optimal" motor control strategy. Over the years, various methods have been developed to address this limitation, such as the stochastic approach. A well-known implementation of this approach, Metabolica, samples a wide number of plausible solutions instead of searching for a single one, leveraging Bayesian statistics and Markov Chain Monte Carlo algorithm, yet allowing muscles to abruptly change their activation levels. To overcome this and other limitations, we developed a new implementation of the stochastic approach (Myobolica), adding constraints and parameters to ensure the identification of physiological solutions. The aim of this study was to evaluate Myobolica, and quantify the differences in terms of width of the solution band (muscle control variability) compared to Metabolica. To this end, both muscle forces and knee joint force solutions bands estimated by the two approaches were compared to one another, and against (i) the solution identified by static optimization and (ii) experimentally measured knee joint forces. The use of Myobolica led to a marked narrowing of the solution band compared to Metabolica. Furthermore, the Myobolica solutions well correlated with the experimental data (R 2 = 0.92 , RMSE = 0.3 BW), but not as much with the optimal solution (R 2 = 0.82 , RMSE = 0.63 BW). Additional analyses are required to confirm the findings and further improve this implementation.

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.

Hip contact forces can be predicted with a neural network using only synthesised key points and electromyography in people with hip osteoarthritis

OBJECTIVE

To develop and validate a neural network to estimate hip contact forces (HCF), and lower body kinematics and kinetics during walking in individuals with hip osteoarthritis (OA) using synthesised anatomical key points and electromyography. To assess the capability of the neural network to detect directional changes in HCF resulting from prescribed gait modifications.

DESIGN

A calibrated electromyography-informed neuromusculoskeletal model was used to compute lower body joint angles, moments, and HCF for 17 participants with mild-to-moderate hip OA. Anatomical key points (e.g., joint centres) were synthesised from marker trajectories and augmented with bias and noise expected from computer vision-based pose estimation systems. Temporal convolutional and long short-term memory neural networks (NN) were trained using leave-one-subject-out validation to predict neuromusculoskeletal modelling outputs from the synthesised key points and measured electromyography data from 5 hip-spanning muscles.

RESULTS

HCF was predicted with an average error of 13.4 ± 7.1% of peak force. Joint angles and moments were predicted with an average root-mean-square-error of 5.3 degrees and 0.10 Nm/kg, respectively. The NN could detect changes in peak HCF that occur due to gait modifications with good agreement with neuromusculoskeletal modelling (r2 = 0.72) and a minimum detectable change of 9.5%.

CONCLUSION

The developed neural network predicted HCF and lower body joint angles and moments in individuals with hip OA using noisy synthesised key point locations with acceptable errors. Changes in HCF magnitude due to gait modifications were predicted with high accuracy. These findings have important implications for implementation of load-modification based gait retraining interventions for people with hip OA in a natural environment (i.e., home, clinic).

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.

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.

A framework based on subject-specific musculoskeletal models and Monte Carlo simulations to personalize muscle coordination retraining

Excessive loads at lower limb joints can lead to pain and degenerative diseases. Altering joint loads with muscle coordination retraining might help to treat or prevent clinical symptoms in a non-invasive way. Knowing how much muscle coordination retraining can reduce joint loads and which muscles have the biggest impact on joint loads is crucial for personalized gait retraining. We introduced a simulation framework to quantify the potential of muscle coordination retraining to reduce joint loads for an individuum. Furthermore, the proposed framework enables to pinpoint muscles, which alterations have the highest likelihood to reduce joint loads. Simulations were performed based on three-dimensional motion capture data of five healthy adolescents (femoral torsion 10°-29°, tibial torsion 19°-38°) and five patients with idiopathic torsional deformities at the femur and/or tibia (femoral torsion 18°-52°, tibial torsion 3°-50°). For each participant, a musculoskeletal model was modified to match the femoral and tibial geometry obtained from magnetic resonance images. Each participant's model and the corresponding motion capture data were used as input for a Monte Carlo analysis to investigate how different muscle coordination strategies influence joint loads. OpenSim was used to run 10,000 simulations for each participant. Root-mean-square of muscle forces and peak joint contact forces were compared between simulations. Depending on the participant, altering muscle coordination led to a maximum reduction in hip, knee, patellofemoral and ankle joint loads between 5 and 18%, 4% and 45%, 16% and 36%, and 2% and 6%, respectively. In some but not all participants reducing joint loads at one joint increased joint loads at other joints. The required alteration in muscle forces to achieve a reduction in joint loads showed a large variability between participants. The potential of muscle coordination retraining to reduce joint loads depends on the person's musculoskeletal geometry and gait pattern and therefore showed a large variability between participants, which highlights the usefulness and importance of the proposed framework to personalize gait retraining.

Osteoarthritis year in review 2023: Biomechanics

Biomechanics plays a significant yet complex role in osteoarthritis (OA) onset and progression. Identifying alterations in biomechanical factors and their complex interactions is critical for gaining new insights into OA pathophysiology and identification of clearly defined and modifiable mechanical treatment targets. This review synthesized biomechanics studies from March 2022 to April 2023, from which three themes relating to human gait emerged: (1) new insights into the pathogenesis of OA using computational modeling and machine learning, (2) technology-enhanced biomechanical interventions for OA, and (3) out-of-lab biomechanical assessments of OA. We further highlighted future-focused areas which may continue to advance the field of biomechanics in OA, with a particular emphasis on exploiting technology to understand and treat biomechanical mechanisms of OA outside the laboratory. The breadth of studies included in this review highlights the complex role of biomechanics in OA and showcase numerous innovative and outstanding contributions to the field. Exciting cross-disciplinary efforts integrating computational modeling, mobile sensors, and machine learning methods show great promise for streamlining in vivo multi-scale biomechanics workflows and are expected to underpin future breakthroughs in the understanding and treatment of biomechanics in OA.

Reducing Knee Joint Loads During Stance Phase With a Rigid-Soft Hybrid Exoskeleton

High mechanical loads generated during walking may accelerate the wear of the knee joint. General knee exoskeletons mainly reduce knee joint load by applying pure torque assistance to reduce the force required by knee extensor muscles. This study aims to further reduce the knee joint load during the early and middle stance phases through a gait intervention strategy that combines torque assistance and vertical force assistance. Comprehensive experiments were conducted to verify the gait intervention strategy in reducing the knee joint load. The results demonstrated that the strategy significantly reduced the maximum and mean knee joint force during the early and middle stance phases. Our studies indicated that the strategy may have the potential to reduce the knee joint load from multiple factors, including muscles, ligaments, and the ground reaction force.

Influence of femoral anteversion angle and neck-shaft angle on muscle forces and joint loading during walking

Femoral deformities, e.g. increased or decreased femoral anteversion (AVA) and neck-shaft angle (NSA), can lead to pathological gait patterns, altered joint loads, and degenerative joint diseases. The mechanism how femoral geometry influences muscle forces and joint load during walking is still not fully understood. The objective of our study was to investigate the influence of femoral AVA and NSA on muscle forces and joint loads during walking. We conducted a comprehensive musculoskeletal modelling study based on three-dimensional motion capture data of a healthy person with a typical gait pattern. We created 25 musculoskeletal models with a variety of NSA (93°-153°) and AVA (-12°-48°). For each model we calculated moment arms, muscle forces, muscle moments, co-contraction indices and joint loads using OpenSim. Multiple regression analyses were used to predict muscle activations, muscle moments, co-contraction indices, and joint contact forces based on the femoral geometry. We found a significant increase in co-contraction of hip and knee joint spanning muscles in models with increasing AVA and NSA, which led to a substantial increase in hip and knee joint contact forces. Decreased AVA and NSA had a minor impact on muscle and joint contact forces. Large AVA lead to increases in both knee and hip contact forces. Large NSA (153°) combined with large AVA (48°) led to increases in hip joint contact forces by five times body weight. Low NSA (108° and 93°) combined with large AVA (48°) led to two-fold increases in the second peak of the knee contact forces. Increased joint contact forces in models with increased AVA and NSA were linked to changes in hip muscle moment arms and compensatory increases in hip and knee muscle forces. Knowing the influence of femoral geometry on muscle forces and joint loads can help clinicians to improve treatment strategies in patients with femoral deformities.

Modeling Human Suboptimal Control: A Review

This review paper provides an overview of the approaches to model neuromuscular control, focusing on methods to identify nonoptimal control strategies typical of populations with neuromuscular disorders or children. Where possible, the authors tightened the description of the methods to the mechanisms behind the underlying biomechanical and physiological rationale. They start by describing the first and most simplified approach, the reductionist approach, which splits the role of the nervous and musculoskeletal systems. Static optimization and dynamic optimization methods and electromyography-based approaches are summarized to highlight their limitations and understand (the need for) their developments over time. Then, the authors look at the more recent stochastic approach, introduced to explore the space of plausible neural solutions, thus implementing the uncontrolled manifold theory, according to which the central nervous system only controls specific motions and tasks to limit energy consumption while allowing for some degree of adaptability to perturbations. Finally, they explore the literature covering the explicit modeling of the coupling between the nervous system (acting as controller) and the musculoskeletal system (the actuator), which may be employed to overcome the split characterizing the reductionist approach.

Latarjet's muscular alterations increase glenohumeral joint stability: A theoretical study

The surgical Latarjet procedure aims to stabilise the glenohumeral joint following anterior dislocations. Despite restoring joint stability, the procedure introduces alterations of muscle paths which likely modify the shoulder dynamics. Currently, these altered muscular functions and their implications are unclear. Hence, this work aims to predict changes in muscle lever arms, muscle and joint forces following a Latarjet procedure by using a computational approach. Planar shoulder movements of ten participants were experimentally assessed. A validated upper-limb musculoskeletal model was utilised in two configurations, i.e., a baseline model, simulating normal joint, and a Latarjet model simulating its related muscular alterations. Muscle lever arms and differences in muscle and joint forces between models were derived from the experimental marker data and static optimisation technique. Lever arms of most altered muscles, hence their role, were substantially changed after Latarjet. Altered muscle forces varied by up to 15% of the body weight. Total glenohumeral joint force increased by up to 14% of the body weight after Latarjet, mostly due to increase in compression force. Our simulation indicated that the Latarjet muscular alterations lead to changes in the muscular recruitment and contribute to the stability of the glenohumeral joint by increasing compression force during planar motions.

Can a Musculoskeletal Model Adapted to Knee Implant Geometry Improve Prediction of 3D Contact Forces and Moments?

Tibiofemoral contact loads are crucial parameters in the onset and progression of osteoarthrosis. While contact loads are frequently estimated from musculoskeletal models, their customization is often limited to scaling musculoskeletal geometry or adapting muscle lines. Moreover, studies have usually focused on superior-inferior contact force without investigating three-dimensional contact loads. Using experimental data from six patients with instrumented total knee arthroplasty (TKA), this study customized a lower limb musculoskeletal model to consider the positioning and the geometry of the implant at knee level. Static optimization was performed to estimate tibiofemoral contact forces and contact moments as well as musculotendinous forces. Predictions from both a generic and a customized model were compared to the instrumented implant measurements. Both models accurately predict superior-inferior (SI) force and abduction-adduction (AA) moment. Notably, the customization improves prediction of medial-lateral (ML) force and flexion-extension (FE) moments. However, there is subject-dependent variability in the prediction of anterior-posterior (AP) force. The customized models presented here predict loads on all joint axes and in most cases improve prediction. Unexpectedly, this improvement was more limited for patients with more rotated implants, suggesting a need for further model adaptations such as muscle wrapping or redefinition of hip and ankle joint centers and axes.

Maintaining soldier musculoskeletal health using personalised digital humans, wearables and/or computer vision

OBJECTIVES

The physical demands of military service place soldiers at risk of musculoskeletal injuries and are major concerns for military capability. This paper outlines the development new training technologies to prevent and manage these injuries.

DESIGN

Narrative review.

METHODS

Technologies suitable for integration into next-generation training devices were examined. We considered the capability of technologies to target tissue level mechanics, provide appropriate real-time feedback, and their useability in-the-field.

RESULTS

Musculoskeletal tissues' health depends on their functional mechanical environment experienced in military activities, training and rehabilitation. These environments result from the interactions between tissue motion, loading, biology, and morphology. Maintaining health of and/or repairing joint tissues requires targeting the "ideal" in vivo tissue mechanics (i.e., loading and strain), which may be enabled by real-time biofeedback. Recent research has shown that these biofeedback technologies are possible by integrating a patient's personalised digital twin and wireless wearable devices. Personalised digital twins are personalised neuromusculoskeletal rigid body and finite element models that work in real-time by code optimisation and artificial intelligence. Model personalisation is crucial in obtaining physically and physiologically valid predictions.

CONCLUSIONS

Recent work has shown that laboratory-quality biomechanical measurements and modelling can be performed outside the laboratory with a small number of wearable sensors or computer vision methods. The next stage is to combine these technologies into well-designed easy to use products.

Gluteal Muscle Forces during Hip-Focused Injury Prevention and Rehabilitation Exercises

PURPOSE

This study aimed to compare and rank gluteal muscle forces in eight hip-focused exercises performed with and without external resistance and describe the underlying fiber lengths, velocities, and muscle activations.

METHODS

Motion capture, ground reaction forces, and electromyography (EMG) were used as input to an EMG-informed neuromusculoskeletal model to estimate gluteus maximus, medius, and minimus muscle forces. Participants were 14 female footballers (18-32 yr old) with at least 3 months of lower limb strength training experience. Each participant performed eight hip-focused exercises (single-leg squat, split squat, single-leg Romanian deadlift [RDL], single-leg hip thrust, banded side step, hip hike, side plank, and side-lying leg raise) with and without 12 repetition maximum (RM) resistance. For each muscle, exercises were ranked by peak muscle force, and k-means clustering separated exercises into four tiers.

RESULTS

The tier 1 exercises for gluteus maximus were loaded split squat (95% confidence interval [CI] = 495-688 N), loaded single-leg RDL (95% CI = 500-655 N), and loaded single-leg hip thrust (95% CI = 505-640 N). The tier 1 exercises for gluteus medius were body weight side plank (95% CI = 338-483 N), loaded single-leg squat (95% CI = 278-422 N), and loaded single-leg RDL (95% CI = 283-405 N). The tier 1 exercises for gluteus minimus were loaded single-leg RDL (95% CI = 267-389 N) and body weight side plank (95% CI = 272-382 N). Peak gluteal muscle forces increased by 28-150 N when exercises were performed with 12RM external resistance compared with body weight only. Peak muscle force coincided with maximum fiber length for most exercises.

CONCLUSIONS

Gluteal muscle forces were exercise specific, and peak muscle forces increased by varying amounts when adding a 12RM external resistance. These findings may inform exercise selection by facilitating the targeting of individual gluteal muscles and optimization of mechanical loads to match performance, injury prevention, or rehabilitation training goals.

Limiting the Use of Electromyography and Ground Reaction Force Data Changes the Magnitude and Ranking of Modelled Anterior Cruciate Ligament Forces

Neuromusculoskeletal models often require three-dimensional (3D) body motions, ground reaction forces (GRF), and electromyography (EMG) as input data. Acquiring these data in real-world settings is challenging, with barriers such as the cost of instruments, setup time, and operator skills to correctly acquire and interpret data. This study investigated the consequences of limiting EMG and GRF data on modelled anterior cruciate ligament (ACL) forces during a drop–land–jump task in late-/post-pubertal females. We compared ACL forces generated by a reference model (i.e., EMG-informed neural mode combined with 3D GRF) to those generated by an EMG-informed with only vertical GRF, static optimisation with 3D GRF, and static optimisation with only vertical GRF. Results indicated ACL force magnitude during landing (when ACL injury typically occurs) was significantly overestimated if only vertical GRF were used for either EMG-informed or static optimisation neural modes. If 3D GRF were used in combination with static optimisation, ACL force was marginally overestimated compared to the reference model. None of the alternative models maintained rank order of ACL loading magnitudes generated by the reference model. Finally, we observed substantial variability across the study sample in response to limiting EMG and GRF data, indicating need for methods incorporating subject-specific measures of muscle activation patterns and external loading when modelling ACL loading during dynamic motor tasks.

Physics-Informed Deep Learning for Musculoskeletal Modeling: Predicting Muscle Forces and Joint Kinematics From Surface EMG

Musculoskeletal models have been widely used for detailed biomechanical analysis to characterise various functional impairments given their ability to estimate movement variables (i.e., muscle forces and joint moments) which cannot be readily measured in vivo. Physics-based computational neuromusculoskeletal models can interpret the dynamic interaction between neural drive to muscles, muscle dynamics, body and joint kinematics and kinetics. Still, such set of solutions suffers from slowness, especially for the complex models, hindering the utility in real-time applications. In recent years, data-driven methods have emerged as a promising alternative due to the benefits in speedy and simple implementation, but they cannot reflect the underlying neuromechanical processes. This paper proposes a physics-informed deep learning framework for musculoskeletal modelling, where physics-based domain knowledge is brought into the data-driven model as soft constraints to penalise/regularise the data-driven model. We use the synchronous muscle forces and joint kinematics prediction from surface electromyogram (sEMG) as the exemplar to illustrate the proposed framework. Convolutional neural network (CNN) is employed as the deep neural network to implement the proposed framework. Simultaneously, the physics law between muscle forces and joint kinematics is used the soft constraint. Experimental validations on two groups of data, including one benchmark dataset and one self-collected dataset from six healthy subjects, are performed. The experimental results demonstrate the effectiveness and robustness of the proposed framework.

ECG Measurement Uncertainty Based on Monte Carlo Approach: An Effective Analysis for a Successful Cardiac Health Monitoring System

Measurement uncertainty is one of the widespread concepts applied in scientific works, particularly to estimate the accuracy of measurement results and to evaluate the conformity of products and processes. In this work, we propose a methodology to analyze the performance of measurement systems existing in the design phases, based on a probabilistic approach, by applying the Monte Carlo method (MCM). With this approach, it is feasible to identify the dominant contributing factors of imprecision in the evaluated system. In the design phase, this information can be used to identify where the most effective attention is required to improve the performance of equipment. This methodology was applied over a simulated electrocardiogram (ECG), for which a measurement uncertainty of the order of 3.54% of the measured value was estimated, with a confidence level of 95%. For this simulation, the ECG computational model was categorized into two modules: the preamplifier and the final stage. The outcomes of the analysis show that the preamplifier module had a greater influence on the measurement results over the final stage module, which indicates that interventions in the first module would promote more significant performance improvements in the system. Finally, it was identified that the main source of ECG measurement uncertainty is related to the measurand, focused towards the objective of better characterization of the metrological behavior of the measurements in the ECG.

Electromyography measurements of the deep hip muscles do not improve estimates of hip contact force

The deep hip muscles are often omitted in studies investigating hip contact forces using neuromusculoskeletal modelling methods. However, recent evidence indicates the deep hip muscles have potential to change the direction of hip contact force and could have relevance for hip contact loading estimates. Further, it is not known whether deep hip muscle excitation patterns can be accurately estimated using neuromusculoskeletal modelling or require measurement (through invasive and time-consuming methods) to inform models used to estimate hip contact forces. We calculated hip contact forces during walking, squatting, and squat-jumping for 17 participants using electromyography (EMG)-informed neuromusculoskeletal modelling with (informed) and without (synthesized) intramuscular EMG for the deep hip muscles (piriformis, obturator internus, quadratus femoris). Hip contact force magnitude and direction, calculated as the angle between hip contact force and vector from femoral head to acetabular center, were compared between configurations using a paired t-test deployed through statistical parametric mapping (P < 0.05). Additionally, root mean square error, correlation coefficients (R²), and timing accuracy between measured and modelled deep hip muscle excitation patterns were computed. No significant between-configuration differences in hip contact force magnitude or direction were found for any task. However, the synthesized method poorly predicted (R²-range 0.02-0.3) deep hip muscle excitation patterns for all tasks. Consequently, intramuscular EMG of the deep hip muscles may be unnecessary when estimating hip contact force magnitude or direction using EMG-informed neuromusculoskeletal modelling, though is likely essential for investigations of deep hip muscle function.