Higher coactivations of lower limb muscles increase stability during walking on slippery ground in forward dynamics musculoskeletal simulation

The energy efficiency theory of human bipedal locomotion has been widely accepted as a neuro-musculoskeletal control method. However, coactivation of agonist and antagonist muscles in the lower limb has been observed during various limb movements, including walking. The emergence of this coactivation cannot be explained solely by the energy efficiency theory and remains a subject of debate. To shed light on this, we investigated the role of muscle coactivations in walking stability using a forward dynamics musculoskeletal simulation combined with neural-network-based gait controllers. Our study revealed that a gait controller with minimal muscle activations had a high probability of falls under challenging gait conditions such as slippery ground and uneven terrain. Lower limb muscle coactivations emerged in the process of gait controller training on slippery ground. Controllers with physiological coactivation levels demonstrated a significantly reduced probability of falls. Our results suggest that achieving stable walking requires muscle coactivations beyond the minimal level of muscle energy. This study implies that coactivations likely emerge to maintain gait stability under challenging conditions, and both coactivation and energy optimization of lower limb muscles should be considered when exploring the foundational control mechanisms of human walking.

Impact forces in backward falls: Subject-specific video-based rigid body simulation of backward falls

A critical missing component in the study of real-world falls is the ability to accurately determine impact forces resulting from the fall. Subject-specific rigid body dynamic (RBD) models calibrated to video captured falls can quantify impact forces and provide additional insights into injury risk factors. RBD models were developed based on five backward falls captured on surveillance video in long-term care facilities in British Columbia, Canada. Model joint stiffness and initial velocities were calibrated to match the kinematics of the fall and contact forces were calculated. The effect of joint stiffnesses (neck, lumbar spine, hip, and knee joint) on head contact forces were determined by modifying the calibrated stiffness values ±25%. Fall duration, fall trajectories, and maximum velocities showed a close match between fall events and simulations. The maximum value of pelvic velocity difference between Kinovea (an open-source software 2D digitization software) and Madymo multibody modeling was found to be 6% ± 21.58%. Our results demonstrate that neck and hip stiffness values have a non-significant yet large effect on head contact force (t(3) = 1, p = 0.387 and t(3) = 2, p = 0.138), while lower effects were observed for knee stiffness, and the effect of lumbar spine stiffness was negligible. The subject-specific fall simulations constructed from real world video captured falls allow for direct quantification of force outcomes of falls and may have applications in improving the assessment of fall-induced injury risks and injury prevention methods.

On challenges in clinical assessment of hip fracture risk using image-based biomechanical modelling: a critical review

INTRODUCTION

Hip fracture is a common health risk among elderly people, due to the prevalence of osteoporosis and accidental fall in the population. Accurate assessment of fracture risk is a crucial step for clinicians to consider patient-by-patient optimal treatments for effective prevention of fractures. Image-based biomechanical modeling has shown promising progress in assessment of fracture risk, and there is still a great possibility for improvement. The purpose of this paper is to identify key issues that need be addressed to improve image-based biomechanical modeling.

MATERIALS AND METHODS

We critically examined issues in consideration and determination of the four biomechanical variables, i.e., risk of fall, fall-induced impact force, bone geometry and bone material quality, which are essential for prediction of hip fracture risk. We closely inspected: limitations introduced by assumptions that are adopted in existing models; deficiencies in methods for construction of biomechanical models, especially for determination of bone material properties from bone images; problems caused by separate use of the variables in clinical study of hip fracture risk; availability of clinical information that are required for validation of biomechanical models.

RESULTS AND CONCLUSIONS

A number of critical issues and gaps were identified. Strategies for effectively addressing the issues were discussed.

The effect of body configuration on the strain magnitude and distribution within the acetabulum during sideways falls: A finite element approach

While the incidence of hip fractures has declined during the last decades, the incidence of acetabular fractures resulting from low-energy sideways falls has increased, and the mechanisms responsible for this trend remain unknown. Previous studies have suggested that body configuration during the impact plays an important role in a hip fracture. Thus, the aim of this study was to investigate the effect of body configuration angles (trunk tilt angle, trunk flexion angle, femur horizontal rotation angle, and femur diaphysis angle) on low-energy acetabular fractures via a parametric analysis. A computed tomography-based (CT) finite element model of the ground-proximal femur-pelvis complex was created, and strain magnitude, time-history response, and distribution within the acetabulum were evaluated. Results showed that while the trunk tilt angle and femur diaphysis angle have the greatest effect on strain magnitude, the direction of the fall (lateral vs. posterolateral) contributes to strain distribution within the acetabulum. The results also suggest that strain level and distribution within the proximal femur and acetabulum resulting from a sideways fall are not similar and, in some cases, even opposite. Taken together, our simulations suggest that a more horizontal trunk and femoral shaft at the impact phase can increase the risk of low-energy acetabular fractures.

Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey

This paper provides a starting point for researchers and practitioners from biology, medicine, physics and engineering who can benefit from an up-to-date literature survey on patient-specific bone fracture modelling, simulation and risk analysis. This survey hints at a framework for devising realistic patient-specific bone fracture simulations. This paper has 18 sections: Section 1 presents the main interested parties; Section 2 explains the organzation of the text; Section 3 motivates further work on patient-specific bone fracture simulation; Section 4 motivates this survey; Section 5 concerns the collection of bibliographical references; Section 6 motivates the physico-mathematical approach to bone fracture; Section 7 presents the modelling of bone as a continuum; Section 8 categorizes the surveyed literature into a continuum mechanics framework; Section 9 concerns the computational modelling of bone geometry; Section 10 concerns the estimation of bone mechanical properties; Section 11 concerns the selection of boundary conditions representative of bone trauma; Section 12 concerns bone fracture simulation; Section 13 presents the multiscale structure of bone; Section 14 concerns the multiscale mathematical modelling of bone; Section 15 concerns the experimental validation of bone fracture simulations; Section 16 concerns bone fracture risk assessment. Lastly, glossaries for symbols, acronyms, and physico-mathematical terms are provided.

Explicit Finite Element Models Accurately Predict Subject-Specific and Velocity-Dependent Kinetics of Sideways Fall Impact

The majority of hip fractures in the elderly are the result of a fall from standing or from a lower height. Current injury models focus mostly on femur strength while neglecting subject-specific loading. This article presents an injury modeling strategy for hip fractures related to sideways falls that takes subject-specific impact loading into account. Finite element models (FEMs) of the human body were used to predict the experienced load and the femoral strength in a single model. We validated these models for their predicted peak force, effective pelvic stiffness, and fracture status against matching ex vivo sideways fall impacts (n = 11) with a trochanter velocity of 3.1 m/s. Furthermore, they were compared to sideways impacts of volunteers with lower impact velocities that were previously conducted by other groups. Good agreement was found between the ex vivo experiments and the FEMs with respect to peak force (root mean square error [RMSE] = 10.7%, R² = 0.85) and effective pelvic stiffness (R² = 0.92, RMSE = 12.9%). The FEMs were predictive of the fracture status for 10 out of 11 specimens. Compared to the volunteer experiments from low height, the FEMs overestimated the peak force by 25% for low BMI subjects and 8% for high BMI subjects. The effective pelvic stiffness values that were derived from the FEMs were comparable to those derived from impacts with volunteers. The force attenuation from the impact surface to the femur ranged between 27% and 54% and was highly dependent on soft tissue thickness (R² = 0.86). The energy balance in the FEMS showed that at the time of peak force 79% to 93% of the total energy is either kinetic or was transformed to soft tissue deformation. The presented FEMs allow for direct discrimination between fracture and nonfracture outcome for sideways falls and bridge the gap between impact testing with volunteers and impact conditions representative of real life falls. © 2019 American Society for Bone and Mineral Research.

Subject-specific ex vivo simulations for hip fracture risk assessment in sideways falls

The risk of hip fracture of a patient due to a fall can be described from a mechanical perspective as the capacity of the femur to withstand the force that it experiences in the event of a fall. So far, impact forces acting on the lateral aspect of the pelvic region and femur strength have been investigated separately. This study used inertia-driven cadaveric impact experiments that mimic falls to the side from standing in order to evaluate the subject-specific force applied to the hip during impact and the fracture outcome in the same experimental model. Eleven fresh-frozen pelvis-femur constructs (6 female, 5 male, age = 77 years (SD = 13 years), BMI = 22.8 kg/m² (SD = 7.8 kg/m²), total hip aBMD = 0.734 g/cm² (SD = 0.149 g/cm²)), were embedded into soft tissue surrogate material that matched subject-specific mass and body shape. The specimens were attached to metallic lower-limb constructions with subject-specific masses and subjected to an inverted pendulum motion. Impact forces were recorded with a 6-axis force plate at 10,000 Hz and three dimensional deflections in the pelvic region were tracked with two high-speed cameras at 5000 Hz. Of the 11 specimens, 5 femur fractures and 3 pelvis fractures were observed. Three specimens did not fracture. aBMD alone did not reliably separate femur fractures from non-fractures. However, a mechanical risk ratio, which was calculated as the impact force divided by aBMD, classified specimens reliably into femur fractures and non-fractures. Single degree of freedom models, based on specimen kinetics, were able to predict subject-specific peak impact forces (RMSE = 2.55% for non-fractures). This study provides direct evidence relating subject-specific impact forces and subject-specific strength estimates and improves the assessment of the mechanical risk of hip fracture for a specific femur/pelvis combination in a sideways fall.

A novel sideways fall simulator to study hip fractures ex vivo

Falls to the side are the leading cause of hip fractures in the elderly. The load that a person experiences during a fall cannot be measured with volunteers for ethical reasons. To evaluate injurious loads, while considering relevant energy input and body posture for a sideways fall, a subject-specific cadaveric impact experiment was developed. Full cadaveric femur-pelvis constructs (N = 2) were embedded in surrogate soft tissue material and attached to metallic surrogate lower limbs. The specimens were then subjected to an inverted pendulum motion, simulating a fall to the side with an impact to the greater trochanter. The load at the ground and the deformation of the pelvis were evaluated using a 6-axis force transducer and two high-speed cameras. Post-test, a trauma surgeon (PG) evaluated specimen injuries. Peak ground contact forces were 7132 N and 5641 N for the fractured and non-fractured specimen, respectively. We observed a cervical fracture of the femur in one specimen and no injuries in a second specimen, showing that the developed protocol can be used to differentiate between specimens at high and low fracture risk.

Sideways fall-induced impact force and its effect on hip fracture risk: a review

Osteoporotic hip fracture, mostly induced in falls among the elderly, is a major health burden over the world. The impact force applied to the hip is an important factor in determining the risk of hip fracture. However, biomechanical researches have yielded conflicting conclusions about whether the fall-induced impact force can be accurately predicted by the available models. It also has been debated whether or not the effect of impact force has been considered appropriately in hip fracture risk assessment tools. This study aimed to provide a state-of-the-art review of the available methods for predicting the impact force, investigate their strengths/limitations, and suggest further improvements in modeling of human body falling.

METHODS

We divided the effective parameters on impact force to two categories: (1) the parameters that can be determined subject-specifically and (2) the parameters that may significantly vary from fall to fall for an individual and cannot be considered subject-specifically.

RESULTS

The parameters in the first category can be investigated in human body fall experiments. Video capture of real-life falls was reported as a valuable method to investigate the parameters in the second category that significantly affect the impact force and cannot be determined in human body fall experiments.

CONCLUSIONS

The analysis of the gathered data revealed that there is a need to develop modified biomechanical models for more accurate prediction of the impact force and appropriately adopt them in hip fracture risk assessment tools in order to achieve a better precision in identifying high-risk patients. Graphical abstract Impact force to the hip induced in sideways falls is affected by many parameters and may remarkably vary from subject to subject.

Study of sex differences in the association between hip fracture risk and body parameters by DXA-based biomechanical modeling

There is controversy about whether or not body parameters affect hip fracture in men and women in the same way. In addition, although bone mineral density (BMD) is currently the most important single discriminator of hip fracture, it is unclear if BMD alone is equally effective for men and women. The objective of this study was to quantify and compare the associations of hip fracture risk with BMD and body parameters in men and women using our recently developed two-level biomechanical model that combines a whole-body dynamics model with a proximal-femur finite element model. Sideways fall induced impact force of 130 Chinese clinical cases, including 50 males and 80 females, were determined by subject-specific dynamics modeling. Then, a DXA-based finite element model was used to simulate the femur bone under the fall-induced loading conditions and calculate the hip fracture risk. Body weight, body height, body mass index, trochanteric soft tissue thickness, and hip bone mineral density were determined for each subject and their associations with impact force and hip fracture risk were quantified. Results showed that the association between impact force and hip fracture risk was not strong enough in both men (r=-0.31,p<0.05) and women (r=0.42,p<0.001) to consider the force as a sole indicator of hip fracture risk. The correlation between hip BMD and hip fracture risk in men (r=-0.83,p<0.001) was notably stronger than that in women (r=-0.68,p<0.001). Increased body mass index was not a protective factor against hip fracture in men (r=-0.13,p>0.05), but it can be considered as a protective factor among women (r=-0.28,p<0.05). In contrast to men, trochanteric soft tissue thickness can be considered as a protective factor against hip fracture in women (r=-0.50,p<0.001). This study suggested that the biomechanical risk/protective factors for hip fracture are sex-specific. Therefore, the effect of body parameters should be considered differently for men and women in hip fracture risk assessment tools. These findings support further exploration of sex-specific preventive and protective measurements to reduce the incidence of hip fractures.

A biomechanical sorting of clinical risk factors affecting osteoporotic hip fracture

Osteoporotic fracture has been found associated with many clinical risk factors, and the associations have been explored dominantly by evidence-based and case-control approaches. The major challenges emerging from the studies are the large number of the risk factors, the difficulty in quantification, the incomplete list, and the interdependence of the risk factors. A biomechanical sorting of the risk factors may shed lights on resolving the above issues. Based on the definition of load-strength ratio (LSR), we first identified the four biomechanical variables determining fracture risk, i.e., the risk of fall, impact force, bone quality, and bone geometry. Then, we explored the links between the FRAX clinical risk factors and the biomechanical variables by looking for evidences in the literature. To accurately assess fracture risk, none of the four biomechanical variables can be ignored and their values must be subject-specific. A clinical risk factor contributes to osteoporotic fracture by affecting one or more of the biomechanical variables. A biomechanical variable represents the integral effect from all the clinical risk factors linked to the variable. The clinical risk factors in FRAX mostly stand for bone quality. The other three biomechanical variables are not adequately represented by the clinical risk factors. From the biomechanical viewpoint, most clinical risk factors are interdependent to each other as they affect the same biomechanical variable(s). As biomechanical variables must be expressed in numbers before their use in calculating LSR, the numerical value of a biomechanical variable can be used as a gauge of the linked clinical risk factors to measure their integral effect on fracture risk, which may be more efficient than to study each individual risk factor.

A two-level subject-specific biomechanical model for improving prediction of hip fracture risk

BACKGROUND

Sideways fall-induced hip fracture is a major worldwide health problem among the elderly population. However, all existing biomechanical models for predicting hip fracture mainly consider the femur related parameters. Their accuracy is limited as hip fracture is significantly affected by loading conditions as well. The objective of this study was to develop a biomechanical model for improving assessment of hip fracture risk by subject-specific prediction of fall-induced loading conditions.

METHOD

All information required to construct the models was extracted from the subject's whole-body and hip medical image in order to make the models subject-specific. Fall-induced hip fracture risk for eighty clinical cases was calculated under two sets of loading conditions: subject-specific determined by the proposed model, and non-subject-specific obtained from empirical functions. The predicted hip fracture risk indices were then compared with clinical observations.

FINDINGS

It was found that the subject-specific prediction of fall-induced loading conditions significantly improves the hip fracture risk assessment. Consistent to the clinical observations, the fracture risk predicted by the proposed model suggested that obesity is a protective factor for hip fracture and underweight subjects are more likely to experience a hip fracture.

INTERPRETATIONS

This study shows that hip fracture risk is affected by a number of factors, including body weight, body height, impact force, body mass index, hip soft tissue thickness, and bone quality. The proposed model provides a comprehensive, fast, accurate, and non-expensive method for prediction of hip fracture risk which should lead to more effective prevention of hip fractures.