A musculoskeletal foot model for clinical gait analysis

A valid novel ground reaction force distribution algorithm to determine midfoot kinetics of gait with a single force plate

BACKGROUND

Understanding midfoot joint kinetics is valuable for improved treatment of foot pathologies. Segmental foot kinetics cannot currently be obtained in a standard gait lab without the use of multiple force plates or a pedobarographic plate overlaid with a force plate due to the single ground reaction force (GRF) vector.

RESEARCH QUESTION

Can an algorithm be created to distribute the GRF into multiple segmental vectors that will allow for calculation of accurate midfoot and ankle moments?

METHODS

20 pediatric subjects (10 typically developing, 10 with foot pathology) underwent multi-segment foot gait analysis using the Milwaukee Foot Model. A novel force distribution algorithm (FDA) was developed to proportionally divide the GRF into segmental vectors based on the location of the center of pressure and these vectors were used to calculate midfoot and ankle moments with an inverse dynamic approach. FDA GRF segmental vectors and midfoot/ankle moments were compared to metrics obtained from a previously validated pedobarographic plate methodology using correlations, statistical parametric mapping, and effect size with 95 % confidence intervals.

RESULTS

All force distributions and midfoot/ankle moments waveforms were highly correlated with R> 0.99 for hindfoot and forefoot forces, R> 0.99 for sagittal and transverse plane moments, and R> 0.95 for coronal plane moments. No statistical differences were found during 2nd rocker where the FDA was applied.

SIGNIFICANCE

Midfoot and ankle moments can be accurately obtained with our algorithm using standard equipment utilized in clinical and research motion analysis labs without the requirement of additional trials, or targeted walking by patients.

Advanced subject-specific neck musculoskeletal modeling unveils sex differences in muscle moment arm and cervical spine loading

Neck pain and injuries are growing healthcare burdens with women having a higher incidence rate and poorer treatment outcomes than males. A better understanding of sex differences in neck biomechanics, foundational for more targeted, effective prevention or treatment strategies, calls for more advanced subject-specific musculoskeletal modeling. Current neck musculoskeletal models are based on generic anatomy, lack subject specificity beyond anthropometric scaling, and are unable to accurately reproduce neck strengths exhibited in vivo without arbitrary muscle force scaling factors or residual torque actuators. In this work, subject-specific neck musculoskeletal models of 23 individuals (11 male, 12 female) were constructed by integrating multi-modality imaging and biomechanical measurements. Each model simulated maximal voluntary neck static exertions in three postures: neck flexion in a neutral posture, flexion in a 40° extended posture, and extension in a 40° flexed posture. Quantitative model validation showed close agreement between model-predicted muscle activation and EMG measurement. The models unveiled that (1) males have greater moment arms in one flexor muscle group and five extensor muscle groups, (2) females exhibited higher cervical spinal compression per unit exertion force in the flexed posture, and (3) the variability of compression force was much greater in females in all three exertions but most notably in the extension with a flexed "dropped head" position. These insights illuminated a plausible pathway from sex differences in neck biomechanics to sex disparities in the risk and prevalence of neck pain.

Virtual reality skateboarding training for balance and functional performance in degenerative lumbar spine disease

BACKGROUND

Degenerative lumbar spine disease (DLD) is a prevalent condition in middle-aged and elderly individuals. DLD frequently results in pain, muscle weakness, and motor impairment, which affect postural stability and functional performance in daily activities. Simulated skateboarding training could enable patients with DLD to engage in exercise with less pain and focus on single-leg weight-bearing. The purpose of this study was to investigate the effects of virtual reality (VR) skateboarding training on balance and functional performance in patients with DLD.

METHODS

Fourteen patients with DLD and 21 age-matched healthy individuals completed a 6-week program of VR skateboarding training. The motion capture and force platform systems were synchronized to collect data during a single-leg stance test (SLST). Musculoskeletal simulation was utilized to calculate muscle force based on the data. Four functional performance tests were conducted to evaluate the improvement after the training. A Visual Analogue Scale (VAS) was also employed for pain assessment.

RESULTS

After the training, pain intensity significantly decreased in patients with DLD (p = 0.024). Before the training, patients with DLD took longer than healthy individuals on the five times sit-to-stand test (p = 0.024). After the training, no significant between-group differences were observed in any of the functional performance tests (p > 0.05). In balance, patients with DLD were similar to healthy individuals after the training, except that the mean frequency (p = 0.014) was higher. Patients with DLD initially had higher biceps femoris force demands (p = 0.028) but shifted to increased gluteus maximus demand after the training (p = 0.037). Gluteus medius strength significantly improved in patients with DLD (p = 0.039), while healthy individuals showed consistent muscle force (p > 0.05).

CONCLUSION

This is the first study to apply the novel VR skateboarding training to patients with DLD. VR skateboarding training enabled patients with DLD to achieve the training effects in a posture that relieves lumbar spine pressure. The results also emphasized the significant benefits to patients with DLD, such as reduced pain, enhanced balance, and improved muscle performance.

Development and validation of a novel ankle joint musculoskeletal model

An improved understanding of contact mechanics in the ankle joint is paramount for implant design and ankle disorder treatment. However, existing models generally simplify the ankle joint as a revolute joint that cannot predict contact characteristics. The current study aimed to develop a novel musculoskeletal ankle joint model that can predict contact in the ankle joint, together with muscle and joint reaction forces. We modelled the ankle joint as a multi-axial joint and simulated contact mechanics between the tibia, fibula and talus bones in OpenSim. The developed model was validated with results from experimental studies through passive stiffness and contact. Through this, we found a similar ankle moment-rotation relationship and contact pattern between our study and experimental studies. Next, the musculoskeletal ankle joint model was incorporated into a lower body model to simulate gait. The ankle joint contact characteristics, kinematics, and muscle forces were predicted and compared to the literature. Our results revealed a comparable peak contact force and the same muscle activation patterns in four major muscles. Good agreement was also found in ankle dorsi/plantar-flexion and inversion/eversion. Thus, the developed model was able to accurately model the ankle joint and can be used to predict contact characteristics in gait.

An investigation of tendon strains in jersey finger injury load cases using a finite element neuromuscular human body model

Introduction: A common hand injury in American football, rugby and basketball is the so-called jersey finger injury (JFI), in which an eccentric overextension of the distal interphalangeal joint leads to an avulsion of the connected musculus flexor digitorum profundus (FDP) tendon. In the field of automotive safety assessment, finite element (FE) neuromuscular human body models (NHBMs) have been validated and are employed to evaluate different injury types related to car crash scenarios. The goal of this study is to show, how such a model can be modified to assess JFIs by adapting the hand of an FE-NHBM for the computational analysis of tendon strains during a generalized JFI load case. Methods: A jersey finger injury criterion (JFIC) covering the injury mechanisms of tendon straining and avulsion was defined based on biomechanical experiments found in the literature. The hand of the Total Human Model for Safety (THUMS) version 3.0 was combined with the musculature of THUMS version 5.03 to create a model with appropriate finger mobility. Muscle routing paths of FDP and musculus flexor digitorum superficialis (FDS) as well as tendon material parameters were optimized using literature data. A simplified JFI load case was simulated as the gripping of a cylindrical rod with finger flexor activation levels between 0% and 100%, which was then retracted with the velocity of a sprinting college football player to forcefully open the closed hand. Results: The optimization of the muscle routing node positions and tendon material parameters yielded good results with minimum normalized mean absolute error values of 0.79% and 7.16% respectively. Tendon avulsion injuries were detected in the middle and little finger for muscle activation levels of 80% and above, while no tendon or muscle strain injuries of any kind occurred. Discussion: The presented work outlines the steps necessary to adapt the hand model of a FE-NHBM for the assessment of JFIs using a newly defined injury criterion called the JFIC. The injury assessment results are in good agreement with documented JFI symptoms. At the same time, the need to rethink commonly asserted paradigms concerning the choice of muscle material parameters is highlighted.

Newly compiled Tai Chi (Bafa Wubu) promotes lower extremity exercise: a preliminary cross sectional study

BACKGROUND

Tai Chi (Bafa Wubu) is a new type of simplified Tai Chi widely practiced by Tai Chi enthusiasts that has developed and perfected simplified Tai Chi movement and enriched Tai Chi practice methods. When practicing, Tai Chi athletes and enthusiasts can choose the Bafa Wubu movements to practice according to their physical conditions. The purpose of this article is to discuss the mechanism by which Bafa Wubu promotes lower extremity exercise from the perspective of exercise biomechanics.

OBJECTIVES

This article aims to explore the scientific training methods and technical characteristics of Bafa Wubu, and its contribution to comprehensive exercise of the lower extremities, by analyzing the biomechanical characteristics of the lower extremities of participants who practice Bafa Wubu at different levels and by comparing their ground reaction force, lower limb joints, and muscles during Bafa Wubu.

METHODS

A total of 16 male participants were recruited and divided into an amateur group (N = 8) and a professional group (N = 8). The data were collected by a BTS 3D infrared-based motion capture system, and Kistler 3D force plate. The lower extremity joint forces and muscle strength were calculated by anybody simulation software with inverse dynamics.

RESULTS

During elbowing and leaning sideways with steps sideways (ELS), the ground reaction force of the professional group was significantly higher than that of the amateur group in the sagittal, vertical, and frontal axes (P < 0.01). While stepping forward, backward, and sideways, the professional group's joints loading at the hip, knee, and ankle was always higher in the vertical direction (P < 0.01). Furthermore, during warding off with steps forward (WOF), laying with steps forward (LF), and rolling back with steps backward (RBB), hip joint loading increased in the med-lat direction. During actions with steps backward and sideways, the professional group's ankle flexion/extension torque and hip abduction/rotation torque were significantly larger than those of the amateur group (P < 0.01). Different actions in Bafa Wubu activate muscles to different degrees, whereas the iliacus is mainly responsible for stabilizing postures when practitioners perform standing knee lifting motions.

CONCLUSIONS

Professional groups who have been practicing Tai Chi (Bafa Wubu) for a long time have higher ground reaction force, and the force on the three joints of the lower extremities is different for various movements, which has positive significance for exercising the joints of the lower extremities. In addition, various motions activate muscles of different types at different levels. For amateurs to practice different movements to stimulate the muscles, targeted areas of practice promote the lower extremity muscles' synergistic force. In summary, the muscles and joints of the lower extremity can obtain comprehensive and balanced exercise through Bafa Wubu.

Study on Two Typical Progressive Motions in Tai Chi (Bafa Wubu) Promoting Lower Extremity Exercise

BACKGROUND

By comparatively investigating the joints, muscles and bones of the lower extremity during two progressive motions in Bafa Wubu and normal walking, this paper aims to enrich the diversity of walking exercise and scientifically provide theoretical guidance for primary practitioners. The scientific training methods and technical characteristics of Bafa Wubu, as well as its contribution to comprehensive exercise of the lower extremities, are further explored.

METHODS

A total of eight professional athletes of Tai Chi at the national level were recruited. The kinetic parameters of the lower extremity were calculated using AnyBody 7.2 musculoskeletal modeling. Stress analysis of the iliac bone was performed using an ANSYS 19.2 workbench.

RESULTS

In Bafa Wubu, the ground reaction force during two progressive motions was significantly smaller than that noted during normal walking. During warding off with steps forward and laying with steps forward, the load at the three joints of the lower extremity was significantly smaller than that during normal walking in the frontal plane, but significantly greater than that noted during normal walking in the vertical axis. In addition, the lower limb joint torque was higher than that of normal walking in both progressive movements, and lower limb muscle activation was higher. The iliac bone loads during the two progressive motions were larger than those during normal walking, and the maximum loading point differed.

CONCLUSIONS

This is the first study to demonstrate the biomechanical performance of Bafa Wubu in professional athletes of Tai Chi. Two progressive motions of Bafa Wubu require the lower extremity to be slowly controlled, thereby resulting in a smaller ground reaction force. In addition, the loads of the three joints at the lower extremity all increase in the vertical direction and decrease in the lateral direction, reducing the possibility of lateral injury to the joints. In addition, the two progressive motions significantly enhance the muscle strength of the plantar flexion muscles, dorsiflexor, and muscles around the thigh, and effectively stimulate the bones of the lower extremity. Therefore, progressive motion training contributes to improving the controlling and supporting capabilities of the lower extremities during normal walking.

Stair-climbing wheelchair proven to maintain user's body stability based on AnyBody musculoskeletal model and finite element analysis

The electric stair-climbing wheelchair is a beneficial mobile assistance device for older adults and disabled persons with poor walking ability, as it reduces the daily walking and climbing burden. In this paper, 11 older adults were tested when using a stair-climbing wheelchair in three environments: flat ground, slopes, and stairs. The kinematic and dynamic parameters of the lower limb joints were simulated by AnyBody 7.2 human model simulation software using Vicon 3D infrared motion capture, a 3D force table, and analyzed by ANSYS 19.2 Workbench. The joint force, joint moment, and muscle strength did not change significantly under the three environments when using the wheelchair. Through finite element analysis of the mechanical properties of the human body, when using the wheelchair, no significant differences in the overall stress distributions of the fifth lumbar spine, hip bone, or femur were found among the three environments, no significant differences in deformation and displacement were found, and the stress distribution was relatively stable. Therefore, the human body is stable enough to use the electric stair-climbing wheelchair in the three test environments, all of which will be commonly encountered in daily life.

Effect of seat tilting on biomechanics of L4 in tractor drivers

The tilting of the tractors' seat during deep tillage operations affects the drivers' sitting position, which can lead to lumbar spine injuries. To investigate the effect of seat tilting on the driver's L4 biomechanics, we built a tractor-driver musculoskeletal model in AnyBody™. The maximum activity of the driver's erector spinae at different tilted angles were measured by sEMG and compared with the simulation results to validate the model. The spatial position of the driver's spine at different tilted angles were obtained by 3 D motion capture. The model simulated the driver's spine posture during the actual tilt and investigated the effects of different tilted angles and vibration on the biomechanics of the driver's L4 . The results showed that as the tilt angle of the tractor increased, the load on the driver's L4 also increased, especially the shear force increased at a faster rate than the axial and normal forces, with the shear force on the driver's L4 increasing from 0 N to 138.7 N when the tractor was tilted from 0° to 15°. When vibration was applied to the musculoskeletal model, the maximum value of the shear force on the driver became progressively greater as the angle of tilt of the tractor increased. Overall, tilting the tractor can have a large impact on the biomechanics of the driver's L4, and tilting the tractor may be an important cause of lumbar spine injuries for tractor drivers.

Applications of Biomechanical Foot Models to Evaluate Dance Movements Using Three-Dimensional Motion Capture: A Review of the Literature

Dance movement requires excessive, repetitive range of motion (ROM) at the foot-ankle complex, possibly contributing to the high rate of injury among dancers. However, we know little about foot biomechanics during dance movements. Researchers are using three-dimensional (3D) motion capture systems to study the in vivo kinematics of joint segments more frequently in dance-medicine research, warranting a literature review and quality assessment evaluation. The purpose of this literature review was to identify and evaluate studies that used 3D motion capture to analyze in vivo biomechanics of the foot and ankle for a cohort of dancers during dance-specific movement. Three databases (PubMed, Ovid MEDLINE, CINAHL) were accessed along with hand searches of dance-specific journals to identify relevant articles through March 2020. Using specific selection criteria, 25 studies were identified. Fifteen studies used single-segment biomechanical foot models originally created to study gait, four used a novel two-segment model, and six utilized a multi-seg- ment foot model. Nine of the studies referenced common and frequently published gait marker sets and four used a dance-specific biomechanical model with purposefully designed foot segments to analyze the dancers' foot and ankle. Description of the biomechanical models varied, reducing the reproducibility of the models and protocols. Investigators concluded that there is little evidence that the extreme total, segmental, and inter-segmental foot and ankle ROM exerted by dancers are being evaluated during dance-specific movements using 3D motion capture. Findings suggest that 3D motion capture is a robust measurement tool that has the capability to assist researchers in evaluating the in vivo, inter-segmental motion of the foot and ankle to potentially discover many of the remaining significant factors predisposing dancers to injury. The literature review synthesis is presented with recommendations for consideration when evaluating results from studies that utilized a 3D biomechanical foot model to evaluate dance-specific movement.

Outcome of Gastrocnemius Soleus Facial Lengthening in Ambulatory Patients With Cerebral Palsy

BACKGROUND

Treatment of equinus contractures in children with cerebral palsy (CP) varies across centers. Existing literature utilizes mixed study populations with a variety of procedures. As such, there is limited knowledge regarding recurrence rates and efficacy of a single procedure performed on a homogenous cohort. Here we retrospectively evaluate outcomes from gastroc soleus fascial lengthenings (GSFL) performed at 2 centers with consistent approaches in both patient selection and operative technique.

METHODS

Subjects meeting inclusion criteria including CP diagnosis, ambulation status, and minimum follow-up criteria were identified. Revision rate was reported based on need for additional calf lengthening procedures. Functional outcomes were evaluated using physical exam measures and selected variables from computational gait analysis. Outcomes factors were identified by comparing revised subjects to unrevised. Longitudinal outcomes of index surgeries were assessed by comparing preoperative functional data to short-term, mid-term, and long-term data.

RESULTS

A total of 64 subjects with 87 limbs met inclusion criteria. In all, 25% of subjects and 21% of limbs went on to revision. Factors influencing revision were age at index surgery and gross motor function classification system (GMFCS) level. More than half of revised limbs had index surgery before age 7. Revision rates for subjects less than 7 were 44% compared with a 17% revision rate for ages 7 to 12, and a 4% revision rate on children older than 12. GMFCSIII subjects had significantly higher revision rates (43%) compared with GMFCSII (18%) and GMFCSI (11%) subjects. Ankle range of motion measures improved significantly with GSFL and most maintained improvements at all time periods. GSFL did not lead to significant calcaneal gait or crouch.

CONCLUSIONS

This study evaluates long term efficacy of GSFL to address equinus in ambulatory children with CP. Overall revision rates are similar to previous reports for GSFL and other calf lengthening procedures. This information may be useful in setting expectations and counselling families. Younger subjects and those with more severe involvement are more likely to need revision surgery, with these factors compounding the likelihood in the younger GMFCSIII child.

LEVEL OF EVIDENCE

Level III-retrospective comparative study.

Dynamics analysis for flexion and extension of elbow joint motion based on musculoskeletal model of Anybody

PURPOSE

Little is known about how biomechanics governs the biological nature for humeral motion dynamically. Elbow motion ought to be investigated based on a musculoskeletal model and evidence the physiologic principle of upper limbs.

METHOD

A humeral model was reconstructed by MIMICS after CT images input in *.dicom format, it was processed by Geomagic Studio for Surfaces, then gridded mesh and assigned materials by Hypermesh. On the other hand, a musculoskeletal model was built by AnyBody, physical motions were then simulated to export boundary condition and myodynamia during flexion and extension. Finally, all the humeral model and boundary were imported to Abaqus for finite element analysis.

RESULT

During the simulative motion of flexion, the primary muscles are brachii biceps, brachialis anticus and teretipronator, their myodynamia increased and then decreased gradually, and reached its peak value at 30°; During extension, the main muscles are triceps brachii and brachialis anticus, their myodynamia increased and then decreased gradually too, and reached peak at 50°; In these two cases, their strain and displacement distributed at the middle of humerus.

CONCLUSION

AnyBody is a novel modelling system to simulate physical motion, for example flexion and extension. Biceps brachii and brachialis anticus are functional for flexion, and triceps brachii plays a key role in extension critically. This simulation confirms the physiologic rule for sport event, humeral fixation and postoperative healing with clinical significance that minimizing joint forces from injury onset may promote pain-free ways.

Extrinsic foot muscle forces and joint contact forces in flexible flatfoot adult with foot orthosis: A parametric study of tibialis posterior muscle weakness

BACKGROUND

The posterior tibialis tendon dysfunction (PTTD) is typically associated with progressive flatfoot deformity, which could be alleviated with foot orthosis. However, the evaluation of tibialis posterior (TP) weakness on lower limb mechanics of flatfoot adults with foot orthoses is scarce and requires further investigation.

RESEARCH QUESTION

This study aimed to examine the effects of TP weakness on lower limb mechanics in flatfoot adults with foot orthosis through gait analysis and musculoskeletal modelling.

METHODS

Fifteen young adults with flatfoot were recruited from University to perform a gait experiment with and without foot orthoses. Data collected from the motion capture system were used to drive the musculoskeletal modelling for the estimation of the joint force and extrinsic muscle forces of the lower limb. A parametric analysis was conducted by adjusting the TP muscle strength from 40 % to 100 %. Two-way repeated measures ANOVA was used to compare the peak extrinsic foot muscle forces and joint forces among different levels of TP weakness and insole conditions.

RESULTS

TP weakness significantly increased ankle joint force superoinferiorly (F = 125.9, p < 0.001) and decreased anteroposteriorly (F = 125.9, p < 0.001), in addition to a significant increase in the muscle forces of flexor hallucis longus (p < 0.001) and flexor digitorum longus (p < 0.001). Besides, the foot orthosis significantly reduced most peak muscle forces whilst significantly reduced the second peak knee force and peak ankle force compared to the control condition (F = 8.79-30.9, p < 0.05).

SIGNIFICANCE

The increased extrinsic foot muscle forces (flexor hallucis longus and flexor digitorum longus) and ankle joint forces in the TP weakness condition indicated that TP weakness may induce compensatory muscle activation and attenuated joint load. The abnormal muscle and joint mechanics in flatfoot adults with TP weakness might be restored by the orthosis.

Anthropometric scaling of musculoskeletal models of the hand captures age-dependent differences in lateral pinch force

Musculoskeletal models and computer simulations enable non-invasive study of muscle function and contact forces. Hand models are useful for understanding the complexities of hand strength, precision movement, and the dexterity required during daily activities. Yet, generic models fail to accurately represent the entire scope of the population, while subject-specific models are labor-intensive to create. The objective of this study was to assess the efficacy of scaled generic models to represent the broad spectrum of strength profiles across the lifespan. We examined one hundred lateral pinch simulations using a generic model of the wrist and thumb anthropometrically scaled to represent the full range of heights reported for four ages across childhood, puberty, older adolescence, and adulthood. We evaluated maximum lateral pinch force produced, muscle control strategies, and the effect of linearly scaling the maximum isometric force. Our simulations demonstrated three main concepts. First, anthropometric scaling could capture age-dependent differences in pinch strength. Second, a generic muscle control strategy is not representative of all populations. Lastly, simulations do not employ optimal fiber length to complete a lateral pinch task. These results demonstrate the potential of anthropometrically-scaled models to study hand strength across the lifespan, while also highlighting that muscle control strategies may adapt as we age. The results also provide insight to the force-length relationship of thumb muscles during lateral pinch. We conclude that anthropometric scaling can accurately represent age characteristics of the population, but subject-specific models are still necessary to represent individuals.

The modified Shriners Hospitals for Children Greenville (mSHCG) multi-segment foot model provides clinically acceptable measurements of ankle and midfoot angles: A dual fluoroscopy study

BACKGROUND

Several multi-segment foot models have been developed to evaluate foot and ankle motion using skin-marker motion analysis. However, few multi-segment models have been evaluated against a reference standard to establish kinematic accuracy.

RESEARCH QUESTION

How accurately do skin-markers estimate foot and ankle motion for the modified Shriners Hospitals for Children Greenville (mSHCG) multi-segment foot model when compared against the reference standard, dual fluoroscopy (DF), during gait, in asymptomatic participants?

METHODS

Five participants walked overground as full-body skin-marker trajectory data and DF images of the foot and shank were simultaneously acquired. Using the mSHCG model, ankle and midfoot angles were calculated throughout stance for both motion analysis techniques. Statistical parametric mapping assessed differences in joint angles and marker positions between skin-marker and DF motion analysis techniques. Paired t tests, and linear regression models were used to compare joint angles and range of motion (ROM) calculated from the two techniques.

RESULTS

In the coronal plane, the skin-marker model significantly overestimated ROM (p = 0.028). Further, the DF model midfoot ROM was significantly positively related to differences between DF and skin-marker midfoot angles (p = 0.035, adjusted R² = 0.76). In the sagittal plane, skin-markers underestimated ankle angles by as much as 7.26°, while midfoot angles were overestimated by as much as 9.01°. However, DF and skin-marker joint angles were not significantly different over stance. Skin-markers on the tibia, calcaneus, and fifth metatarsal had significantly different positions than the DF markers along the direction of walking for isolated portions that were less than 10 % of stance. Euclidean distances between DF and skin-markers positions were less than 9.36 mm.

SIGNIFICANCE

As the accuracy of the mSHCG model was formerly unknown, the results of this study provide ranges of confidence for key angles calculated by this model.

Strength Training Characteristics of Different Loads Based on Acceleration Sensor and Finite Element Simulation

Deep squat, bench press and hard pull are important ways for people to improve their strength. The use of sensors to measure force is rare. Measuring strength with sensors is extremely valuable for people to master the intensity of exercise to scientifically effective exercise. To this end, in this paper, we used a real-time wireless motion capture and mechanical evaluation system of the wearable sensor to measure the dynamic characteristics of 30 young men performing deep squat, bench press and hard pull maneuvers. The data of tibia were simulated with AnyBody 5.2 and ANSYS 19.2 to verify the authenticity. The result demonstrated that the appropriate force of the deep squat elbow joint, the hip joint and the knee joint is 40% 1RM, the appropriate force of the bench press is 40% 1RM and the appropriate force of the hard pull is 80% 1RM. The external force is the main factor of bone change. The mechanical characteristics of knee joint can be simulated after the Finite Element Analysis and the simulation of AnyBody model are verified.

Relationships between the foot posture Index and static as well as dynamic rear foot and arch variables

Clinicians, podiatrists and researchers have been quantifying foot posture and movement in various speed conditions and populations. Common variables to assess foot posture/movement are the Foot Posture Index (FPI-6), Achilles tendon angle (β), rear foot angle (γ) and longitudinal arch angle (LAA). These variables were frequently used in clinical and biomechanical settings. This study aimed to determine the relationship between the biomechanical variables (β, γ & LAA) in static and dynamic conditions and the clinically used FPI-6 and their redundancy. Forty participants performed bipedal standing, over ground walking and running trials. Manual assessment data (FPI-6), kinematic data and ground reaction forces were collected. Discrete biomechanical variables (β, γ & LAA) were calculated at various time points (e.g. heel strike). A Principal Component Analysis (PCA) was performed to quantify the contribution of each variable to the overall variance in the data set. Spearman correlations were used to assess the relationship between the sub-measures of the FPI-6 and the biomechanical variables. Two major components were found that explained 85.2% of the overall variance, consisting of LAA and β variables, respectively. Only LAA variables showed significant, but moderate correlations (r < -0.6) with the fifth sub-measurement of the FPI-6. The LAA and β describe independent movements, which dominate foot posture/movement in static and dynamic conditions. The FPI-6 sub-measurements did not closely reflect static nor dynamic behavior of the rear and mid foot. The FPI-6 and biomechanical variables may not be used interchangeably for screening or grouping individuals according to their foot posture/movement.

Joint moments and contact forces in the foot during walking

The net force and moment of a joint have been widely used to understand joint disease in the foot. Meanwhile, it does not reflect the physiological forces on muscles and contact surfaces. The objective of the study is to estimate active moments by muscles, passive moments by connective tissues and joint contact forces in the foot joints during walking. Joint kinematics and external forces of ten healthy subjects (all males, 24.7 ± 1.2 years) were acquired during walking. The data were entered into the five-segment musculoskeletal foot model to calculate muscle forces and joint contact forces of the foot joints using an inverse dynamics-based optimization. Joint reaction forces and active, passive and net moments of each joint were calculated from muscle and ligament forces. The maximum joint reaction forces were 8.72, 4.31, 2.65, and 3.41 body weight (BW) for the ankle, Chopart's, Lisfranc and metatarsophalangeal joints, respectively. Active and passive moments along with net moments were also obtained. The maximum net moments were 8.6, 8.4, 5.4 and 0.8%BW∙HT, respectively. While the trend of net moment was very similar between the four joints, the magnitudes and directions of the active and passive moments varied between joints. The active and passive moments during walking could reveal the roles of muscles and ligaments in each of the foot joints, which was not obvious in the net moment. This method may help narrow down the source of joint problems if applied to clinical studies.

A methodological framework for detecting ulcers' risk in diabetic foot subjects by combining gait analysis, a new musculoskeletal foot model and a foot finite element model

Diabetic foot is one of the most debilitating complications of diabetes and may lead to plantar ulcers. In the last decade, gait analysis, musculoskeletal modelling (MSM) and finite element modelling (FEM) have shown their ability to contribute to diabetic foot prevention and suggested that the origin of the plantar ulcers is in deeper tissue layers rather than on the plantar surface. Hence the aim of the current work is to develop a methodology that improves FEM-derived foot internal stresses prediction, for diabetic foot prevention applications. A 3D foot FEM was combined with MSM derived force to predict the sites of excessive internal stresses on the foot. In vivo gait analysis data, and an MRI scan of a foot from a healthy subject were acquired and used to develop a six degrees of freedom (6 DOF) foot MSM and a 3D subject-specific foot FEM. Ankle kinematics were applied as boundary conditions to the FEM together with: 1. only Ground Reaction Forces (GRFs); 2. OpenSim derived extrinsic muscles forces estimated with a standard OpenSim MSM; 3. extrinsic muscle forces derived through the (6 DOF) foot MSM; 4. intrinsic and extrinsic muscles forces derived through the 6 DOF foot MSM. For model validation purposes, simulated peak pressures were extracted and compared with those measured experimentally. The importance of foot muscles in controlling plantar pressure distribution and internal stresses is confirmed by the improved accuracy in the estimation of the peak pressures obtained with the inclusion of intrinsic and extrinsic muscle forces.

Dynamic Flexion Stiffness of Foot Joints During Walking

BACKGROUND

Dynamic stiffness can be used for studying foot pathologic abnormalities and for developing prostheses and orthoses. Although previous works have studied the role of ankle joint stiffness during gait, other foot joints have not yet been analyzed. We sought to characterize the dynamic stiffness of the ankle, midtarsal, and metatarsophalangeal joints during normal walking.

METHODS

Kinematics and contact data from four healthy individuals during walking were registered with a three-dimensional motion analysis system and a pressure platform. Stance phases with flexion moment-angle linear relationships were identified, and dynamic stiffnesses were calculated from the slope of their linear regressions. Intraparticipant repeatability was analyzed using analyses of variance, and interparticipant variability was checked through the SD of averaged participant stiffnesses.

RESULTS

Flexion moment-angle linear relationships were identified (R(2) > 0.98) during the early and late midstance phases and the propulsion phase at the ankle (2.76, 5.23, and 3.42 N·m/kg/rad, respectively) and midtarsal (15.88, 3.90, and 4.64 N·m/kg/rad, respectively) joints. At the metatarsophalangeal joint, a linear relationship (R(2) > 0.96) occurred only during the propulsion phase (0.11 N·m/kg/rad). High dynamic stiffness variability was observed during the late and early midstance phases at the ankle and midtarsal joints, respectively.

CONCLUSIONS

These results may serve as a basis for future studies aimed at investigating the role of dynamic stiffness identified herein in different foot disorders. The importance of properly controlling the samples in such studies is highlighted. Study of the dynamic stiffnesses identified might be used in the design of prostheses, orthoses, and other assistive devices.

A Patient-Specific Foot Model for the Estimate of Ankle Joint Forces in Patients with Juvenile Idiopathic Arthritis

Juvenile idiopathic arthritis (JIA) is the leading cause of childhood disability from a musculoskeletal disorder. It generally affects large joints such as the knee and the ankle, often causing structural damage. Different factors contribute to the damage onset, including altered joint loading and other mechanical factors, associated with pain and inflammation. The prediction of patients’ joint loading can hence be a valuable tool in understanding the disease mechanisms involved in structural damage progression. A number of lower-limb musculoskeletal models have been proposed to analyse the hip and knee joints, but juvenile models of the foot are still lacking. This paper presents a modelling pipeline that allows the creation of juvenile patient-specific models starting from lower limb kinematics and foot and ankle MRI data. This pipeline has been applied to data from three children with JIA and the importance of patient-specific parameters and modelling assumptions has been tested in a sensitivity analysis focused on the variation of the joint reaction forces. This analysis highlighted the criticality of patient-specific definition of the ankle joint axes and location of the Achilles tendon insertions. Patient-specific detection of the Tibialis Anterior, Tibialis Posterior, and Peroneus Longus origins and insertions were also shown to be important.

CONTRIBUTION OF MUSCLE TENSION FORCE TO MEDIAL KNEE CONTACT FORCE AT FAST WALKING SPEED

Fast walking is considered as a factor that causes pain in patients suffering from knee disorders. This study examined the effect of walking speed on the medial knee contact force and identified contributions to the muscle tension on the medial knee contact force during fast walking using musculoskeletal simulation analysis. The muscle contribution to the medial knee contact force was calculated based on the joint angles and ground reaction force for the normal and fast walking experiments of seven subjects. The muscle force and joint reaction force were used to estimate the medial knee contact force. Results showed, in average, 70% increase in medial knee contact force at the first peak and 34% increase at the second peak with a fast walking speed, compared to when they walked at a normal walking speed. The remarkable increase in the first peak was mainly contributed by the increase in the quadriceps force resisting the external knee flexion moment. In contrast, the moderate increase of second peak was contributed by the increase in the gastrocnemius muscle force. These results suggest that the increase in medial knee contact force at fast walking speeds is caused by the increased muscle force.

Statistical shape and appearance models of bones

When applied to bones, statistical shape models (SSM) and statistical appearance models (SAM) respectively describe the mean shape and mean density distribution of bones within a certain population as well as the main modes of variations of shape and density distribution from their mean values. The availability of this quantitative information regarding the detailed anatomy of bones provides new opportunities for diagnosis, evaluation, and treatment of skeletal diseases. The potential of SSM and SAM has been recently recognized within the bone research community. For example, these models have been applied for studying the effects of bone shape on the etiology of osteoarthritis, improving the accuracy of clinical osteoporotic fracture prediction techniques, design of orthopedic implants, and surgery planning. This paper reviews the main concepts, methods, and applications of SSM and SAM as applied to bone.

Discrepancies in anthropometric parameters between different models affect intervertebral rotations when loading finite element models with muscle forces from inverse static analyses

In only a few published finite element (FE) simulations have muscle forces been applied to the spine. Recently, muscle forces determined using an inverse static (IS) model of the spine were transferred to a spinal FE model, and the effect of methodical parameters was investigated. However, the sensitivity of anthropometric differences between FE and IS models, such as body height and spinal orientation, was not considered. The aim of this sensitivity study was to determine the influence of those differences on the intervertebral rotations (IVRs) following the transfer of muscle forces from an IS model to a FE model. Muscle forces were estimated for 20° flexion and 10° extension of the upper body using an inverse static musculoskeletal model. These forces were subsequently transferred to a nonlinear FE model of the spino-pelvic complex, which includes 243 muscle fascicles. Deviations of body height (±10 cm), spinal orientation in the sagittal plane (±10°), and body weight (±10 kg) between both models were intentionally generated, and their influences on IVRs were determined. The changes in each factor relative to their corresponding reference value of the IS model were calculated. Deviations in body height, spinal orientation, and body weight resulted in maximum changes in the IVR of 19.2%, 26% and 4.2%, respectively, relative to T12-S1 IVR. When transferring muscle forces from an IS to a FE model, it is crucial that both models have the same spinal orientation and height. Additionally, the body weight should be equal in both models.

Can modular strategies simplify neural control of multidirectional human locomotion?

Each human lower limb contains over 50 muscles that are coordinated during locomotion. It has been hypothesized that the nervous system simplifies muscle control through modularity, using neural patterns to activate muscles in groups called synergies. Here we investigate how simple modular controllers based on invariant neural primitives (synergies or patterns) might generate muscle activity observed during multidirectional locomotion. We extracted neural primitives from unilateral electromyographic recordings of 25 lower limb muscles during five locomotor tasks: walking forward, backward, leftward and rightward, and stepping in place. A subset of subjects also performed five variations of forward (unidirectional) walking: self-selected cadence, fast cadence, slow cadence, tiptoe, and uphill (20% incline). We assessed the results in the context of dimensionality reduction, defined here as the number of neural signals needing to be controlled. For an individual task, we found that modular architectures could theoretically reduce dimensionality compared with independent muscle control, but we also found that modular strategies relying on neural primitives shared across different tasks were limited in their ability to account for muscle activations during multi- and unidirectional locomotion. The utility of shared primitives may thus depend on whether they can be adapted for specific task demands, for instance, by means of sensory feedback or by being embedded within a more complex sensorimotor controller. Our findings indicate the need for more sophisticated formulations of modular control or alternative motor control hypotheses in order to understand muscle coordination during locomotion.

A computational model for dynamic analysis of the human gait

Biomechanical models are important tools in the study of human motion. This work proposes a computational model to analyse the dynamics of lower limb motion using a kinematic chain to represent the body segments and rotational joints linked by viscoelastic elements. The model uses anthropometric parameters, ground reaction forces and joint Cardan angles from subjects to analyse lower limb motion during the gait. The model allows evaluating these data in each body plane. Six healthy subjects walked on a treadmill to record the kinematic and kinetic data. In addition, anthropometric parameters were recorded to construct the model. The viscoelastic parameter values were fitted for the model joints (hip, knee and ankle). The proposed model demonstrated that manipulating the viscoelastic parameters between the body segments could fit the amplitudes and frequencies of motion. The data collected in this work have viscoelastic parameter values that follow a normal distribution, indicating that these values are directly related to the gait pattern. To validate the model, we used the values of the joint angles to perform a comparison between the model results and previously published data. The model results show a same pattern and range of values found in the literature for the human gait motion.

Influence of joint models on lower-limb musculo-tendon forces and three-dimensional joint reaction forces during gait

Several three-dimensional (3D) lower-limb musculo-skeletal models have been developed for gait analysis and different hip, knee and ankle joint models have been considered in the literature. Conversely to the influence of the musculo-tendon geometry, the influence of the joint models - i.e. number of degrees of freedom and passive joint moments - on the estimated musculo-tendon forces and 3D joint reaction forces has not been extensively examined.

In this paper musculo-tendon forces and 3D joint reaction forces have been estimated for one subject and one gait cycle with nine variations of a musculoskeletal model and outputs have been compared to measured electromyographic signals and knee joint contact forces. The model outputs are generally in line with the measured signals. However, the 3D joint reaction forces were higher than published values and the contact forces measured for the subject. The results of this study show that, with more degrees of freedom in the model, the musculo-tendon forces and the 3D joint reaction forces tend to increase but with some redistribution between the muscles. In addition, when taking into account passive joint moments, the 3D joint reaction forces tend to decrease during the stance phase and increase during the swing phase.

Although further investigations are needed, a five-degree-of-freedom lower-limb musculo-skeletal model with some angle-dependent joint coupling and stiffness seems to provide satisfactory musculo-tendon forces and 3D joint reaction forces.

The jump shot - a biomechanical analysis focused on lateral ankle ligaments

Handball is one of the top four athletic games with highest injury risks. The jump shot is the most accomplished goal shot technique and the lower extremities are mostly injured. As a basis for ankle sprain simulation, the aim of this study was to extend the ankle region of an existing musculoskeletal full-body model through incorporation of three prominent lateral ankle ligaments: ligamentum fibulotalare anterius (LFTA), ligamentum fibulotalare posterius (LFTP), ligamentum fibulocalcaneare (LFC). The specific objective was to calculate and visualise ligament force scenarios during the jumping and landing phases of controlled jump shots. Recorded kinematic data of performed jump shots and the corresponding ground reaction forces were used to perform inverse dynamics. The calculated peak force of the LFTA (107 N) was found at maximum plantarflexion and of the LFTP (150 N) at maximum dorsiflexion. The peak force of the LFC (190 N) was observed at maximum dorsiflexion combined with maximum eversion. Within the performed jump shots, the LFTA showed a peak force (59 N to 69 N) during maximum plantarflexion in the final moment of the lift off. During landing, the force developed by the LFTA reached its peak value (61 N to 70 N) at the first contact with the floor. After that, the LFTP developed a peak force (70 N to 118 N). This model allows the calculation of forces in lateral ankle ligaments. The information obtained in this study can serve as a basis for future research on ankle sprain and ankle sprain simulation.

Generation of subject-specific, dynamic, multisegment ankle and foot models to improve orthotic design: a feasibility study

Background

Currently, custom foot and ankle orthosis prescription and design tend to be based on traditional techniques, which can result in devices which vary greatly between clinicians and repeat prescription. The use of computational models of the foot may give further insight in the biomechanical effects of these devices and allow a more standardised approach to be taken to their design, however due to the complexity of the foot the models must be highly detailed and dynamic.

Methods/Design

Functional and anatomical datasets will be collected in a multicentre study from 10 healthy participants and 15 patients requiring orthotic devices. The patient group will include individuals with metarsalgia, flexible flat foot and drop foot.

Each participant will undergo a clinical foot function assessment, 3D surface scans of the foot under different loading conditions, and detailed gait analysis including kinematic, kinetic, muscle activity and plantar pressure measurements in both barefoot and shod conditions. Following this each participant will undergo computed tomography (CT) imaging of their foot and ankle under a range of loads and positions while plantar pressures are recorded. A further subgroup of participants will undergo magnetic resonance imaging (MRI) of the foot and ankle.

Imaging data will be segmented to derive the geometry of the bones and the orientation of the joint axes. Insertion points of muscles and ligaments will be determined from the MRI and CT-scans and soft tissue material properties computed from the loaded CT data in combination with the plantar pressure measurements. Gait analysis data will be used to drive the models and in combination with the 3D surface scans for scaling purposes. Predicted plantar pressures and muscle activation patterns predicted from the models will be compared to determine the validity of the models.

Discussion

This protocol will lead to the generation of unique datasets which will be used to develop linked inverse dynamic and forward dynamic biomechanical foot models. These models may be beneficial in predicting the effect of and thus improving the efficacy of orthotic devices for the foot and ankle.