Search

Mechanical performance of clavicle fracture fixation plates is often evaluated using finite element (FE) analysis. Typically, these studies use simplified loading conditions and assume a transversal fracture orientation. However, the loading conditions and fracture orientation influence how the fracture site and thus fixation plate is loaded. In this study, a musculoskeletal model that included the clavicle muscles and scapulohumeral rhythm was defined based on previously published models. The standard OpenSim workflow (inverse kinematics, inverse dynamics, static optimization, and joint reaction analysis) was used to calculate muscle and joint contact forces based on 3D marker data collected in three subjects during seven activities of daily living (ADL). These loading conditions were then applied to a 3D clavicle model with three different fracture orientations and the mean resulting moments on both fragments were calculated to assess fracture stability. Magnitude of glenohumeral contact forces showed good agreement with instrumented shoulder prosthesis data, whereas simulated muscle activations were comparable to experimental EMG data. An oblique fracture orienting from superomedial to inferolateral was the least self-stabilizing. The loading to which the clavicle is exposed during ADL tasks is more complex than the simplified loading conditions typically used as boundary conditions in FE analyses of clavicle fracture fixation plates. Additionally, transversal fractures did not represent the least self-stabilizing fracture orientation, and thus calculated stresses in the plate could be underestimated. Therefore, more complex loading conditions and evaluation of a midshaft fracture running from superomedial to inferolateral is more relevant in FE analyses of midshaft clavicle fracture fixation plates. ispartof: JOURNAL OF ORTHOPAEDIC RESEARCH vol:40 issue:9 pages:2179-2188 ispartof: location:United States status: published

BACKGROUND: Osteoarthritis is a highly prevalent disease affecting the hip and knee joint and is characterized by load-mediated pain and decreased quality of life. Dependent on involved joint, patients present antalgic movement compensations, aiming to decrease loading on the involved joint. However, the associated alterations in mechanical loading of the ipsi- and contra-lateral lower limb joints, are less documented. Here, we documented the biomechanical fingerprint of end-stage hip and knee osteoarthritis patients in terms of ipsilateral and contralateral hip and knee loading during walking and stair ambulation. METHODS: Three-dimensional motion-analysis was performed in 20 hip, 18 knee osteoarthritis patients and 12 controls during level walking and stair ambulation. Joint contact forces were calculated using a standard musculoskeletal modelling workflow in Opensim. Involved and contralateral hip and knee joint loading was compared against healthy controls using independent t-tests (p

OBJECTIVE: This study proposes a computationally efficient method to quantify the effect of surgical inaccuracies on ligament strain in total knee arthroplasty (TKA). More specifically, this study describes a framework to determine the implant position and required surgical accuracy that results in a ligament balanced post-operative outcome with a probability of 90%. METHODS: The response surface method is used to translate uncertainty in the implant position parameters to uncertainty in the ligament strain. The designed uncertainty quantification technique allows for an optimization with feasible computational cost towards the planned implant position and the tolerated surgical error for each of the twelve degrees of freedom of the implant position. RESULTS: It is shown that the error does not allow for a ligament balanced TKA with a probability of 90% using preoperative planning. Six critical implant position parameters can be identified, namely AP translation, PD translation, VV rotation, IE rotation for the femoral component and PD translation, VV rotation for the tibial component. CONCLUSION: We introduced an optimization process that allows for the computation of the required surgical accuracy for a ligament balanced postoperative outcome using preoperative planning with feasible computational cost. SIGNIFICANCE: Towards the research society, the proposed method allows for a computationally efficient uncertainty quantification on a complex model. Towards surgical technique developers, six critical implant position parameters were identified, which should be the focus when refining surgical accuracy of TKA, leveraging better patient satisfaction. ispartof: IEEE Transactions On Biomedical Engineering vol:68 issue:11 pages:3273-3280 ispartof: location:United States status: published

Current fixation plates used to operatively stabilize clavicular fractures are suboptimal, leading to reoperation rates up to 53%. Plate irritation, which can be caused by a poor geometric fit and plate thickness, has been found to be an important factor for reoperation. Moreover, muscle attachment sites (MAS) have to be detached occasionally. To improve current surgical clavicle fracture treatment with plate osteosynthesis, a change in plate design is required. The goal of this study was to design a patient-specific clavicle fracture fixation plate that includes geometrical optimization and stiffness determination. Its biomechanical performance has been compared with a commercial plate by examining the geometric fit, anatomical outline, stresses and interfragmentary motion using a finite element analysis with physiological loading and boundary conditions. Evaluation showed a better geometrical fit of the patient-specific plate as well as an improved fracture reduction. Displacements between fracture fragments were lower in case of the patient-specific plate, both when a fracture gap and no fracture gap were present. Results indicate a superior mechanical performance in terms of all investigated outcomes of the patient-specific plate compared to the commercial plate, while better aligning with the patient-specific geometry and without the need for MAS release. Due to the patient-specific geometry and reduced thickness, these fixation plates are expected to decrease the operation time, as intraoperative contouring will become irrelevant, and to decrease reoperation rates as implant irritation will be minimized. ispartof: JOURNAL OF ORTHOPAEDIC RESEARCH vol:40 issue:7 pages:1695-1706 ispartof: location:United States status: published

BACKGROUND: Musculoskeletal modelling is used to assess musculoskeletal loading during gait. Linear scaling methods are used to personalize generic models to each participant's anthropometry. This approach introduces simplifications, especially when used in paediatric and/or pathological populations. This study aimed to compare results from musculoskeletal simulations using various models ranging from linear scaled to highly subject-specific models, i.e., including the participant's musculoskeletal geometry and electromyography data. METHODS: Magnetic resonance images (MRI) and gait data of one typically developing child and three children with cerebral palsy were analysed. Musculoskeletal simulations were performed to calculate joint kinematics, joint kinetics, muscle forces and joint contact forces using four modelling frameworks: 1) Generic-scaled model with static optimization, 2) Generic-scaled model with an electromyography-informed approach, 3) MRI-based model with static optimization, and 4) MRI-based model with an electromyography-informed approach. FINDINGS: Root-mean-square-differences in joint kinematics and kinetics between generic-scaled and MRI-based models were below 5° and 0.15 Nm/kg, respectively. Root-mean-square-differences over all muscles was below 0.2 body weight for every participant. Root-mean-square-differences in joint contact forces between the different modelling frameworks were up to 2.2 body weight. Comparing the simulation results from the typically developing child with the results from the children with cerebral palsy showed similar root-mean-square-differences for all modelling frameworks. INTERPRETATION: In our participants, the impact of MRI-based models on joint contact forces was higher than the impact of including electromyography. Clinical reasoning based on overall root-mean-square-differences in musculoskeletal simulation results between healthy and pathological participants are unlikely to be affected by the modelling choice. ispartof: CLINICAL BIOMECHANICS vol:87 ispartof: location:England status: published

BACKGROUND: Children with cerebral palsy (CP) present with a pathological gait pattern due to musculoskeletal impairments, such as muscle weakness and altered bony geometry. However, the effect of these impairments on gait performance is still unknown. Research aim:This study aimed to explore the effect of hip muscle weakness and femoral deformities on the gait performance of CP and typical developing (TD) subjects. METHODS: 6400 musculoskeletal models were created by weakening the hip extensors, abductors, adductors and flexors from 0% to 75 % and increasing the femoral anteversion angle (FAA) and neck shaft angle (NSA) from 20° to 60° and 120° to 160°, respectively. One TD and five CP gait patterns were imposed to each model and muscle forces were calculated. The effect of weakness and bony deformities on the capability gap (CG) at the hip, i.e. the lack in hip moment generating capacity to perform the gait pattern, was investigated using regression analysis. RESULTS: The CG of apparent equinus, stiff knee gait, TD gait, jump gait and true equinus increased with 0.080, 0.038, 0.015, 0.023 and 0.005 Nm/kg per 10 percent hip abductor weakness increase, with 0.211, 0.130, 0.056, 0.045 and 0.011 Nm/kg per 10 degrees FAA increase and with 0.163, 0.080, 0.036, 0.043 and 0.011 Nm/kg per 10 degrees NSA increase, respectively. Combined weakness and bony deformities explained 96 %, 85 %, 82 %, 65 %, 40 % and 35 % of the variance in the CG of apparent equinus, TD gait, stiff knee gait, jump gait, true equinus and crouch gait, respectively. SIGNIFICANCE: The results suggest that surgical correction of femoral deformities is more likely to be effective than strength training of hip muscles in enhancing CP gait performance. Jump gait, true equinus and especially crouch were more robust, while apparent equinus and stiff knee gait were limited by hip weakness and femoral deformities. ispartof: GAIT & POSTURE vol:83 pages:280-286 ispartof: location:England status: published

Orthognathic surgery is performed to realign the jaws of a patient through several osteotomies. The state-of-the-art bone plates used to maintain the bone fragments in place are made of titanium. The presence of these non-degradable plates can have unwanted side effects on the long term (e.g. higher infection risk) if they are not removed. Using a biodegradable material such as magnesium may be a possible solution to this problem. However, biodegradation leads to a decrease of mechanical strength, therefore a time-dependent computational approach can help to evaluate the performance of such plates. In the present work, a computational framework has been developed to include biodegradation and bone healing algorithms in a finite element model. It includes bone plates and the mandible, which are submitted to masticatory loads during the early healing period (two months following the surgery). Two different bone plate designs with different stiffnesses have been tested. The stiff design exhibited good mechanical stability, with maximum Von Mises stress being less than 40% of the yield strength throughout the simulation. The flexible design shows high stresses when the bone healing has not started in the fracture gaps, indicating possible failure of the plate. However, this design leads to a higher bone healing quality after two months, as more cartilage is formed due to higher strains exerted in fracture gaps. We therefore conclude that in silico modelling can support tuning of the design parameters to ensure mechanical stability and while promoting bone healing.

BACKGROUND: The purpose of this study was to compare muscle forces and hip contact forces (HCFs) during squatting in patients with cam-type femoroacetabular impingement (cam-FAI) before and after hip corrective surgery and with healthy control participants. METHODS: Ten symptomatic male patients with cam-FAI performed deep squatting preoperatively and at 2 years postoperatively. Patients were matched by age and body mass index to 10 male control participants. Full-body kinematics and kinetics were computed, and muscle forces and HCFs were estimated using a musculoskeletal model and static optimization. Normalized squat cycle (%SC) trials were compared using statistical nonparametric mapping (SnPM). RESULTS: Postoperatively, patients with cam-FAI squatted down with higher anterior pelvic tilt, higher hip flexion, and greater hip extension moments than preoperatively. Preoperative patients demonstrated lower anterior pelvic tilt and lower hip flexion compared with the participants in the control group. Postoperative patients showed increased semimembranosus force compared with their preoperative values. Preoperative forces were lower than the control group for the adductor magnus, the psoas major, and the semimembranosus; however, the preoperative patients showed greater inferior gluteus maximus forces than the patients in the control group, whereas the postoperative patients did not differ from the control patients. Higher posterior, superior, and resultant HCF magnitudes were identified postoperatively in comparison with the preoperative values. Preoperative posterior HCF was lower than in the control group, whereas the postoperative posterior HCF did not differ from those in the control group. CONCLUSIONS: Higher postoperative anterior pelvic tilt was associated with an indication of return to closer to normal pelvic motion, which resembled data from the control group. Lower preoperative anterior pelvic tilt was associated with muscle force imbalance, indicated by decreased semimembranosus and increased gluteus maximus forces. The overall increased postoperative muscle forces were associated with improved pelvic mobility and increased HCFs that were comparable with the control-group standards. CLINICAL RELEVANCE: Muscle forces and HCFs may be indicative of postoperative joint health restoration and alleviated symptoms. ispartof: JOURNAL OF BONE AND JOINT SURGERY-AMERICAN VOLUME vol:102 issue:21S pages:34-42 ispartof: location:United States status: published

Custom reverse shoulder implants represent a valuable solution for patients with large bone defects. Since each implant has unique patient-specific features, finite element (FE) analysis has the potential to guide the design process by virtually comparing the stability of multiple configurations without the need of a mechanical test. The aim of this study was to develop an automated virtual bench test to evaluate the initial stability of custom shoulder implants during the design phase, by simulating a fixation experiment as defined by ASTM F2028-14. Three-dimensional (3D) FE models were generated to simulate the stability test and the predictions were compared to experimental measurements. Good agreement was found between the baseplate displacement measured experimentally and determined from the FE analysis (Spearman&rsquo, s rank test, p &lt, 0.05, correlation coefficient &rho, s = 0.81). Interface micromotion analysis predicted good initial fixation (micromotion &lt, 150 &micro, m, commonly used as bone ingrowth threshold). In conclusion, the finite element model presented in this study was able to replicate the mechanical condition of a standard test for a custom shoulder implants.

BACKGROUND: Many patients with cerebral palsy present a pathologic gait pattern, which presumably induces aberrant musculoskeletal loading that interferes with natural bone growth, causing bone deformations on the long term. Botulinum toxin interventions and single-event multilevel surgeries are used to restore the gait pattern, assuming that a normal gait pattern restores musculoskeletal loading and thus prevents further bone deformation. However, it is unknown if these interventions are able to restore musculoskeletal loading. Hence, we investigated the impact of botulinum toxin injections and single-event multilevel surgery on musculoskeletal loading. METHODS: Gait data collected in 93 children with bilateral cerebral palsy, which included pre- and post multi-level botulinum toxin (49 children) and single-event multilevel surgery (44 children) assessments, and 15 typically developing children were retrospectively processed using a musculoskeletal modelling workflow to calculate joint angles, moments, muscle and joint contact force magnitudes and orientations. Differences from the typically developing waveform were expressed by a root-mean square difference were compared using paired t-tests for each intervention separately (alpha

Background Apart from biomechanical alterations in movement patterns, it is known that movement limitations in persons with knee osteoarthritis (PwKOA) are related to an individual’s perception and belief regarding pain and disability. To gain more insights into the functional movement behaviour of PwKOA in a clinical setting, inertial sensor technology can be applied. This study first aims to evaluate the ability of inertial sensors to discriminate between healthy controls (HC) and PwKOA. Secondly, this study aims to determine the relationship between movement behaviour, pain-related factors and disability scores. Methods Twelve HC and 19 PwKOA were included. Five repetitions of six functional movement tasks (walking, forward lunge, sideward lunge, ascent and descent stairs, single leg squat and sit-to-stand) were simultaneously recorded by the inertial sensor system and a camera-based motion analysis system. Statistically significant differences in angular waveforms of the trunk, pelvis and lower limb joints between HC and PwKOA were determined using one-dimensional statistical parametric mapping (SPM1D). The Knee injury and Osteoarthritis Outcome Score and TAMPA scale for Kinesiophobia were used to evaluate the relationship between discriminating joint motion, pain-related factors and disability using spearman’s correlation coefficients. Results PwKOA had significantly less trunk rotation, internal pelvis rotation and knee flexion ROM during walking. Additionally, the reduced knee flexion (i.e. at the end of the stance phase and swing phase) was related to increased level of perceived pain. During the sideward lunge, PwKOA had significantly less knee flexion, ankle plantarflexion and hip abduction. This decreased hip abduction (i.e. during stance) was related to higher fear of movement. Finally, PwKOA had significantly less knee flexion during the forward lunge, single leg squat and during ascent and descent stairs. No significant correlations were observed with disability. Conclusions Inertial sensors were able to discriminate between movement characteristics of PwKOA and HC. Additionally, significant relationships were found between joint motion, perceived pain and fear of movement. Since inertial sensors can be used outside the laboratory setting, these results are promising as they indicate the ability to evaluate movement deviations. Further research is required to enable measurements of small movement deviations in clinically relevant tasks.

It is well founded that the mechanical environment may regulate bone regeneration in orthopedic applications. The purpose of this study is to investigate the mechanical contributions of the scaffold and the host to bone regeneration, in terms of subject specificity, implantation site and sensitivity to the mechanical environment. Using a computational approach to model mechano-driven regeneration, bone ingrowth in porous titanium scaffolds was simulated in the distal femur and proximal tibia of three goats and compared to experimental results. The results showed that bone ingrowth shifted from a homogeneous distribution pattern, when scaffolds were in contact with trabecular bone (max local ingrowth 12.47%), to a localized bone ingrowth when scaffolds were implanted in a diaphyseal location (max local ingrowth 20.64%). The bone formation dynamics revealed an apposition rate of 0.37±0.28%/day in the first three weeks after implantation, followed by limited increase in bone ingrowth until the end of the experiment (12 weeks). According to in vivo data, we identified one animal whose sensitivity to mechanical stimulation was higher than the other two. Moreover, we found that the stimulus initiating bone formation was consistently higher in the femur than in the tibia for all the individuals. Overall, the dependence of the osteogenic response on the host biomechanics means that, from a mechanical perspective, the regenerative potential depends on both the scaffold and the host environment. Therefore, this work provides insights on how the mechanical conditions of both the recipient and the scaffold contribute to meet patient and location-specific characteristics. ispartof: BONE vol:144 ispartof: location:United States status: published

Adequate deltoid and rotator cuff lengthening in total shoulder arthroplasty (TSA) is crucial to maximize the postoperative functional outcome and to avoid complications. Hence surgeons and patients could benefit from including muscle length information in preoperative planning software.Although different methods have been introduced to automatically indicate patient-specific muscle attachment and wrapping points, the definition of a fast and accurate workflow is still a challenge, due to the large variability in bone shapes. Therefore, the goal of this study is to develop and evaluate the accuracy of a novel method to automatically estimate muscle elongation, based on a statistical shape modelling (SSM) approach.Firstly, humerus and scapula SSMs were used to automatically indicate the attachment points of the main shoulder muscles: subscapularis, supraspinatus, infraspinatus, teres minor and deltoid. Secondly, a wrapping algorithm was applied to identify the points where muscles wrap around bones or potential implants. Finally, the accuracy of the automatically indicated landmarks and its effect on the muscle elongation were evaluated by comparing the manually indicated landmarks with the landmarks identified through the SSM for a set of 40 healthy shoulder CT-scans.The low errors on elongation values suggest that the presented automated workflow is a promising tool for allowing surgeons to evaluate patient-specific muscle elongations during preoperative planning. Although the evaluation was limited to healthy joints, this method allows to easily process large datasets and to potentially find a correlation between muscle elongations and postoperative outcome.

Total knee arthroplasty (TKA) is a common procedure that has become the standard of treatment for severe cases of knee osteoarthritis. Biomechanics and quality of movement similar to healthy were found to improve patient-reported outcomes.In this study, an evaluated musculoskeletal model predicted ligament, contact and muscle forces together with secondary tibiofemoral kinematics. An artificial neural network applied to the musculoskeletal model searched for the optimal implant position in a given range that will minimize the root-mean-square-error (RMSE) between post- TKA and native experimental tibiofemoral kinematics during a squat.We found that, using a cruciate-retaining implant, native kinematics could be accurately reproduced (average RMSE 1.47 mm (± 0.89 mm) for translations and 2.89° (± 2.83°) for rotations between native and optimal TKA alignment). The required implant positions changes maximally 2.96 mm and 2.40o. This suggests that when using pre- operative planning, off-the-shelf CR implants allow for reproducing native knee kinematics post-operatively.

Multi-scale simulations, combining muscle and joint contact force (JCF) from musculoskeletal simulations with adaptive mechanobiological finite element analysis, allow to estimate musculoskeletal loading and predict femoral growth in children. Generic linearly scaled musculoskeletal models are commonly used. This approach, however, neglects subject- and age-specific musculoskeletal geometry, e.g. femoral neck-shaft angle (NSA) and anteversion angle (AVA). This study aimed to evaluate the impact of proximal femoral geometry, i.e. altered NSA and AVA, on hip JCF and femoral growth simulations. Musculoskeletal models with NSA ranging from 120° to 150° and AVA ranging from 20° to 50° were created and used to calculate muscle and hip JCF based on the gait analysis data of a typically developing child. A finite element model of a paediatric femur was created from magnetic resonance images. The finite element model was morphed to the geometries of the different musculoskeletal models and used for mechanobiological finite element analysis to predict femoral growth trends. Our findings showed that hip JCF increase with increasing NSA and AVA. Furthermore, the orientation of the hip JCF followed the orientation of the femoral neck axis. Consequently, the osteogenic index, which is a function of cartilage stresses and defines the growth rate, barely changed with altered NSA and AVA. Nevertheless, growth predictions were sensitive to the femoral geometry due to changes in the predicted growth directions. Altered NSA had a bigger impact on the growth results than altered AVA. Growth simulations based on mechanobiological principles were in agreement with reported changes in paediatric populations. ispartof: Plos One vol:15 issue:7 ispartof: location:United States status: Published online

Objective Early detection of degenerative changes in the cartilage matrix composition is essential for evaluating early interventions that slow down osteoarthritis (OA) initiation. T1rho and T2 relaxation times were found to be effective for detecting early changes in proteoglycan and collagen content. To use these magnetic resonance imaging (MRI) methods, it is important to document the topographical variation in cartilage thickness, T1rho and T2 relaxation times in a healthy population. As OA is partially mechanically driven, the relation between these MRI-based parameters and localized mechanical loading during walking was investigated. Design MR images were acquired in 14 healthy adults and cartilage thickness and T1rho and T2 relaxation times were determined. Experimental gait data was collected and processed using musculoskeletal modeling to identify weight-bearing zones and estimate the contact force impulse during gait. Variation of the cartilage properties (i.e., thickness, T1rho, and T2) over the femoral cartilage was analyzed and compared between the weight-bearing and non-weight-bearing zone of the medial and lateral condyle as well as the trochlea. Results Medial condyle cartilage thickness was correlated to the contact force impulse ( r = 0.78). Lower T1rho, indicating increased proteoglycan content, was found in the medial weight-bearing zone. T2 was higher in all weight-bearing zones compared with the non-weight-bearing zones, indicating lower relative collagen content. Conclusions The current results suggest that medial condyle cartilage is adapted as a long-term protective response to localized loading during a frequently performed task and that the weight-bearing zone of the medial condyle has superior weight bearing capacities compared with the non-weight-bearing zones.

This study aims at defining gait pathomechanics in patients with hip osteoarthritis (OA) and their effect on hip joint loading by combining analyses of hip kinematics, kinetics and contact forces during gait. Twenty patients with hip OA and 17 healthy volunteers matched for age and BMI performed three-dimensional gait analysis. Hip OA level was evaluated based on plane radiographs using the Tonnis classification. Hip joint kinematics, kinetics as well as hip contact forces were calculated. Waveforms were time normalized and compared between groups using statistical parametric mapping analysis. Patients walked with reduced hip adduction angle and reduced hip abduction and external rotation moments. The work generated by the hip abductors during the stance phase of gait was largely decreased. These changes resulted in a decrease and a more vertical and anterior orientation of the hip contact forces compared to healthy controls. This study documents alterations in hip kinematics and kinetics resulting in decreased hip loading in patients with hip OA. The results suggested that patients altered their gait to increase medio-lateral stability, thereby decreasing demand on the hip abductors. These findings support discharge of abductor muscles that may bear clinical relevance of tailored rehabilitation targeting hip abductor muscles strengthening and gait retraining.

Model-based simulations of walking have the theoretical potential to support clinical decision making by predicting the functional outcome of treatments in terms of walking performance. Yet before using such simulations in clinical practice, their ability to identify the main treatment targets in specific patients needs to be demonstrated. In this study, we generated predictive simulations of walking with a medical imaging based neuro-musculoskeletal model of a child with cerebral palsy presenting crouch gait. We explored the influence of altered muscle-tendon properties, reduced neuromuscular control complexity, and spasticity on gait function in terms of joint kinematics, kinetics, muscle activity, and metabolic cost of transport. We modeled altered muscle-tendon properties by personalizing Hill-type muscle-tendon parameters based on data collected during functional movements, simpler neuromuscular control by reducing the number of independent muscle synergies, and spasticity through delayed muscle activity feedback from muscle force and force rate. Our simulations revealed that, in the presence of aberrant musculoskeletal geometries, altered muscle-tendon properties rather than reduced neuromuscular control complexity and spasticity were the primary cause of the crouch gait pattern observed for this child, which is in agreement with the clinical examination. These results suggest that muscle-tendon properties should be the primary target of interventions aiming to restore a more upright gait pattern for this child. This suggestion is in line with the gait analysis following muscle-tendon property and bone deformity corrections. The ability of our simulations to distinguish the contribution of different impairments on walking performance opens the door for identifying targeted treatment strategies with the aim of designing optimized interventions for neuro-musculoskeletal disorders.

Finite element modeling (FEM) can predict hip cartilage contact mechanics. This study investigated how subject-specific boundary conditions and joint geometry affect acetabular cartilage contact mechanics using a multi-scale workflow. For two healthy subjects, musculoskeletal models calculated subject-specific hip kinematics and loading, which were used as boundary conditions for FEM. Cartilage contact mechanics were predicted using either generic or subject-specific FEM and boundary conditions. A subject-specific mesh resulted in a more lateral contact. Effects of subject-specific boundary conditions varied between both subjects. Results highlight the complex interplay between loading and kinematics and their effect on cartilage contact mechanics. ispartof: COMPUTER METHODS IN BIOMECHANICS AND BIOMEDICAL ENGINEERING vol:22 issue:16 pages:1323-1333 ispartof: location:England status: published

This study aimed to find an optimal measurement protocol of elbow and forearm kinematics using different modelling and tracking methods. Kinematic data of elbow flexion/extension and forearm pronation/supination was acquired using optical motion capture from 12 healthy male volunteers. Segment coordinate systems for humerus, forearm, radius, ulna, and hand were defined. Different tracking methods, using anatomical markers or rigid or point maker clusters, were used to compute the three-dimensional rotations. Marker placement errors were assessed to evaluate the rigid body assumption. Multiple comparisons demonstrated statistical differences between tracking methods: compared to using only anatomical markers, tracking using clusters reduced the estimated range of pronation/supination by 14.9%-43.2%, while it estimated increased flexion/extension by 5.3%-9.1%. The study suggests using only anatomical markers exerts the optimal estimation of elbow and forearm kinematics. Modelling using the coordinate systems of the humerus and forearm and of the humerus and ulna, respectively, demonstrated good consistency with literature and are correspondingly regarded as the most appropriate approach for measuring pronation/supination and flexion/extension. The results are valuable in establishing a measurement protocol for analysing elbow and forearm kinematics, avoiding confusions and misinterpretations in communicating results from different methodologies. ispartof: PROCEEDINGS OF THE INSTITUTION OF MECHANICAL ENGINEERS PART H-JOURNAL OF ENGINEERING IN MEDICINE vol:233 issue:11 pages:1113-1121 ispartof: location:England status: published

BACKGROUND: Deviations in knee joint geometry and alignment were previously related to an increased risk for knee OA. These were hypothesized to influence the load distribution over the articular cartilage. Therefore, this study evaluated the effect of altered knee joint geometry and alignment in the coronal and transverse plane on the medial-lateral load distribution and ligament strain using a musculoskeletal modeling approach. METHODS: Joint kinematics during gait were measured in 15 healthy adults. Using different musculoskeletal models with altered geometry of the tibia plateau or knee joint malalignment in the coronal and transverse plane, the resulting muscle, ligament and contact forces were calculated. Next, the distribution of the load over the medial and lateral condyle was analyzed and compared to the reference loading distribution, with neutral geometry and alignment, using repeated-measures ANOVA and individual t-tests, with a Bonferroni-corrected alpha level. RESULTS: Coronal plane malalignment significantly affected the load distribution. Small changes in coronal tibial slope had less pronounced effects on the load distribution, but increased ligament strains. Transverse plane malalignment only minimally affected the load distribution. CONCLUSION: Coronal plane knee malalignment affected knee loading, with increased varus alignment resulting in increased medial loading. This confirms a causal relation between coronal malalignment and increased medial compartment loading and suggests a potential role of aberrant coronal plane alignment on OA initiation. Altered coronal tibial slope induced increased ligament strains, potentially contributing to a cascade of knee laxity and subsequently more extreme knee malalignment. ispartof: KNEE vol:26 issue:4 pages:813-823 ispartof: location:Netherlands status: published

BACKGROUND: Children with cerebral palsy (CP) present altered gait patterns and electromyography (EMG) activity compared to typically developing children. To temporarily reduce muscular activity and to correct the abnormal muscle force balance, Botulinum Toxin type A (BTX-A) injections are used. RESEARCH QUESTION: What is the effect of BTX-A injections on dynamic muscle forces during gait, when calculated using an EMG-constrained approach?. METHODS: Retrospective data of ten typically developing (TD) and fourteen children with spastic diplegic CP were used for musculoskeletal modeling and dynamic simulations of gait, before and after BTX-A treatment. Individual muscle forces were calculated using an EMG-constrained optimization, in which EMG of eight muscles was used as muscle excitation signal to constrain the muscle activation patterns. Paired t-tests were used to compare average modelled muscle forces in different phases of the gait cycle pre- and post-BTX-A, summarized in the muscle profile score. Two-sample t-tests were used to determine significant differences between TD and pre- and post-BTX-A modelled muscle forces. RESULTS: For most muscles, the force was decreased in CP compared to TD children in all phases of the gait cycle, both before and after BTX-A treatment. Differences in muscle forces before and after BTX-A treatment were limited, with only few significant differences between pre- and post-BTX-A. Compared to a standard static optimization approach, imposing the EMG activity increased modelled muscle forces for most muscles. SIGNIFICANCE: Our findings indicate that BTX-A treatment has a limited effect on the muscle balance in CP children. Besides that, the use of EMG-constrained optimization is recommended when studying muscle balance in children with CP. ispartof: GAIT & POSTURE vol:82 pages:54-60 ispartof: location:England status: published

BACKGROUND: Selective dorsal rhizotomy aims to reduce spasticity in children with cerebral palsy. Early investigations indicated postoperative weakness, whereas more recent studies showed that selective dorsal rhizotomy either does not change or improves muscle strength. All previous studies assessed muscle strength in a static position, which did not represent the walking situation. The aim of this study was to analyze the influence of selective dorsal rhizotomy on muscle forces during gait. METHODS: Motion capture data of 25 children with spastic cerebral palsy and 10 typically developing participants were collected. A musculoskeletal OpenSim model was used to calculate joint kinematics, joint kinetics and muscle forces during gait. Static optimization and an electromyography-informed approach to calculate muscle forces were compared. A Muscle-Force-Profile was introduced and used to compare the muscle forces during walking before and after a selective dorsal rhizotomy. FINDINGS: Independent of the approach used (electromyography-informed versus static optimization), selective dorsal rhizotomy significantly normalized forces in spastic muscles during walking and did not reduce the contribution of non-spastic muscles. INTERPRETATION: This study showed that selective dorsal rhizotomy improves dynamic muscle forces in children with cerebral palsy and leads to less gait pathology, as shown in the improvement in joint kinematics and joint kinetics. Individual muscle force analyses using the Muscle-Force-Profile extend standard joint kinematics and joint moment analyses, which might improve clinical-decision making in children with cerebral palsy in the future. The reference data of our participants and MATLAB code for the Muscle-Force-Profile are publicly available on simtk.org/projects/muscleprofile. ispartof: CLINICAL BIOMECHANICS vol:65 pages:26-33 ispartof: location:England status: published

To account for proximal femoral deformities in children with cerebral palsy (CP), subject-specific musculoskeletal models are needed. Non-rigid deformation (NRD) deforms generic onto personalized bone geometry and thereby transforms the muscle points. The goal of this study was to determine to what extent the models and simulation outcomes in CP patients differ when including subject-specific detail using NRD or Magnetic Resonance Imaging (MRI)-based models. The NRD models slightly overestimated hip contact forces compared to MRI models and differences in muscle point positions and moment arm lengths (MALs) remained, although differences were smaller than for the generic model. ispartof: Computer Methods In Biomechanics And Biomedical Engineering vol:22 issue:4 pages:376-385 ispartof: location:England status: published

Total knee arthroplasty (TKA) is a reliable surgical procedure, yet up to a fifth of primary implant patients remains unsatisfied. Musculoskeletal modeling (MSM) has the potential to explore the relationship between implant alignment and functional outcome [3]. Consequently, implant alignment can be quantitatively optimized to restore the pre- TKA joint behavior and, therefore, achieve the most favorable functional outcome for the specific patient. For this reason, we developed a method to optimize the implant alignment, with the aim of restoring the native kinematics and ligament elongations of the patient before undergoing TKA. Subject-specific optimization towards ligament elongations demonstrated to accurately emulate the pre-TKA ligament behavior, in contrast to the mechanically aligned approach. However, the values of the optimized implant positions resulting from the pre-TKA kinematic optimization were extreme in some cases. The presented modelling approach is a promising starting point for allowing surgeons to evaluate the patient-specific implant alignment and restore the patient- specific biomechanics.

Children with cerebral palsy (CP) often present aberrant hip geometry, more specifically increased femoral anteversion and neck-shaft angle. Furthermore, altered gait patterns are present within this population. This study analyzed the effect of aberrant femoral geometry, as present in subjects with CP, on muscle control of pelvis orientation, hip and knee joint kinematics. We ran dynamic simulations of CP and normal gait using two musculoskeletal models, one reflecting normal femoral geometry and one reflecting proximal femoral deformities. The results show that the combination of aberrant bone geometry and CP-specific gait characteristics beneficially increased the potentials of gluteus medius and maximus to control the hip and knee extension. In contrast, the potentials of the hamstrings to extend the hip decreased but to flex the knee increased. In conclusion, this study emphasizes the concomitant effect of the presence of proximal femoral deformity and CP gait characteristics on the muscle control of pelvis orientation and hip and knee joint kinematics during single stance. Not accounting for subject-specific geometry will affect the calculated muscles’ potential during gait. Therefore, the use of generic models to assess muscle function in the presence of femoral deformity and CP gait should be treated with caution. publisher: Elsevier articletitle: The role of altered proximal femoral geometry in impaired pelvis stability and hip control during CP gait: A simulation study journaltitle: Gait & Posture articlelink: http://dx.doi.org/10.1016/j.gaitpost.2015.11.010 content_type: article copyright: Copyright © 2015 Elsevier B.V. All rights reserved. ispartof: Gait & Posture vol:44 pages:61-67 ispartof: location:England status: published

Several approaches may be used for hip replacement surgery either in combination with conventional total hip arthroplasty (THA) or resurfacing hip arthroplasty (RHA). This study investigates the differences in hip loading during gait one year or more after surgery in three cohorts presenting different surgical procedures, more specific RHA placed using the direct lateral (RHA-DLA, n=8) and posterolateral (RHA-PLA, n=14) approach as well as THA placed using the direct anterior (THA-DAA, n=12) approach. For the DAA and control subjects, hip loading was also evaluated during stair ascent and descent to evaluate whether these motions can better discriminate between patients and controls compared to gait. Musculoskeletal modeling in Opensim was used to calculate in vivo joint loading. Results showed that for all operated patients, regardless the surgical procedure, hip loading was decreased compared to control subjects, while no differences were found between patient groups. This indicates that THA via DAA results in similar hip loading as a RHA via DLA or PLA. Stair climbing did not result in more distinct differences in hip contact force magnitude between patients and controls, although differences in orientation were more distinct. However, patients after hip surgery did adjust their motion pattern to decrease the magnitude of loading on the hip joint compared to control subjects. publisher: Elsevier articletitle: Does surgical approach or prosthesis type affect hip joint loading one year after surgery? journaltitle: Gait & Posture articlelink: http://dx.doi.org/10.1016/j.gaitpost.2015.11.009 content_type: article copyright: Copyright © 2015 Elsevier B.V. All rights reserved. ispartof: Gait & Posture vol:44 pages:74-82 ispartof: location:England status: published

Following metal-on-metal hip arthroplasty, edge loading (i.e., loading near the edge of a prosthesis cup) can increase wear and lead to early revision. The position and coverage angle of the prosthesis cup influence the risk of edge loading. This study investigates the effect of altered gait patterns, more specific hip, and pelvis kinematics, on the orientation of hip contact force and the consequent risk of antero-superior edge loading using muscle driven simulations of gait. With a cup orientation of 25° anteversion and 50° inclination and a coverage angle of 168°, many gait patterns presented risk of edge loading. Specifically at terminal double support, 189 out of 405 gait patterns indicated a risk of edge loading. At this time instant, the high hip contact forces and the proximity of the hip contact force to the edge of the cup indicated the likelihood of the occurrence of edge loading. Although the cup position contributed most to edge loading, altering kinematics considerably influenced the risk of edge loading. Increased hip abduction, resulting in decreasing hip contact force magnitude, and decreased hip extension, resulting in decreased risk on edge loading, are gait strategies that could prevent edge loading. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1069-1076, 2016.

BackgroundThe unipodal stance task is a clinical task that quantifies postural stability and alignment of the lower limb joints, while weight bearing on one leg. As persons with knee osteoarthritis (KOA) have poor postural and knee joint stability, objective assessment of this task might be useful.ObjectiveTo investigate the discriminant validity of three-dimensional joint kinematics and centre of mass displacement (COM) between healthy controls and persons with knee KOA, during unipodal stance using inertial sensors. Additionally, the reliability, agreement and construct validity are assessed to determine the reproducibility and accuracy of the discriminating parameters.MethodsTwenty healthy controls and 19 persons with unilateral severe KOA were included. Five repetitions of the unipodal stance task were simultaneously recorded by an inertial sensor system and a camera-based system (gold standard). Statistical significant differences in kinematic waveforms between healthy controls and persons with severe knee KOA were determined using one-dimensional statistical parametric mapping (SPM1D).ResultsPersons with severe knee KOA had more lateral trunk lean towards the contralateral leg, more hip flexion throughout the performance of the unipodal stance task, more pelvic obliquity and COM displacement towards the contralateral side. However, for the latter two parameters the minimum detectable change was greater than the difference between healthy controls and persons with severe knee KOA. The construct validity was good (coefficient of multiple correlation 0.75, 0.83 respectively) and the root mean squared error (RMSE) was low (RMSE

We present a simulation platform that will enable clinicians to evaluate the effect of different treatment options on gait performance in children with cerebral palsy (CP) in order to select the treatment with the highest potential to normalize the patient’s gait pattern. We present a case study to demonstrate the use of the platform. We created a neuro-musculoskeletal model of a 10-year old female child with mild spastic triplegic CP (GMFCS II) who was treated with single-event multilevel surgery based on medical imaging and motion capture data collected before the surgery. Based on this model, we predicted that the treatment would reduce the capability gap, i.e. the torque deficit of the patient with respect to the joint torques needed for normal walking. This prediction was in accordance with the closer-to-normal post-treatment gait kinetics of the child.

1.. Introduction: Gait analysis together with musculoskeletal modeling can be used to calculate muscle forces and assess pathological gait [1]. No generic, pediatric musculoskeletal models are available and, therefore, linear scaling methods are commonly used to personalize a generic, adult musculoskeletal model to the child’s anthropometry. 2. Research: How different are joint kinematics, joint kinetics and muscle force estimates of generic scaled models compared to medical-imaging based models in typically developing (TD) children? 3. Methods: 3D motion capture data and magnetic resonance images (MRI) of a TD boy (age: 8 years; height: 1.23 m; weight: 20.4 kg) were collected. Two musculoskeletal OpenSim models were created: (1) a scaled generic model (M_gen), and (2) a MRI-based model, which included subject-specific musculoskeletal geometry (M_mri) [2]. Joint kinematics, joint kinetics and muscle forces were calculated for each model using OpenSim 3.3 [3]. Joint kinematics, joint kinetics, muscle force waveforms, as well as femoral anteversion angle, neck-shaft angle and hip joint centre location were compared between both models. 4. Results: Joint kinematics and joint kinetics were surprisingly similar between the M_gen and M_mri with root-mean-square-differences below 2.8° and 0.05Nm/kg for joint angles and moments, respectively (Fig. 1, Fig. 2). Depending on the analyzed muscle, differences in muscle forces varied substantially (up to 230% difference) between the M_gen and M_mri (Fig. 3). Femoral anteversion and neck-shaft angles differed between M_gen and M_mri by 12 and 5 degrees, respectively. The hip joint centre position differed between both models by 5, 15 and 6 mm in the anterior/posterior, superior/inferior and medial/lateral direction, respectively.

ispartof: pages:41-42 ispartof: Gait Posture vol:65 Suppl 1 pages:41-42 ispartof: location:England status: published

Stair motion in the presence of hip osteoarthritis (OA) has received less attention than level walking. Its more strenuous aspect may shed the light on different locomotor strategies when compared to walking. We, therefore, aimed to define stair motion features associated to hip OA and to evaluate whether these specific features would differ from level walking and better characterize the hip pathological condition. Principal component and linear discriminant analyses were, respectively, used as data reduction and classification techniques. Our study highlighted that most of stair motion features associated to hip OA were similar to the ones of walking. Stair descent presented with the lowest misclassification error rate, ranging from 12% to 19% (estimated by cross-validation). But, features that may be considered as a mechanism to reduce demand on the hip abductors were found to be more important in the stair ascent condition. This was reflected by both, greater importance in the classification rule and variance compared with walking, that is, decreased hip internal rotation moment at mid-stance (72.50% vs. 57.63%) and increased trunk lateroflexion toward affected side (56.43% vs. 29.37%). This study emphasized the importance of investigating stair motion in hip osteoarthritic population by highlighting specific locomotor strategies. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res

Exercise plays a pivotal role in maximizing peak bone mass in adulthood and maintaining it through aging, by imposing mechanical loading on the bone that can trigger bone mineralization and growth. The optimal type and intensity of exercise that best enhances bone strength remains, however, poorly characterized, partly because the exact peak loading of the bone produced by the diverse types of exercises is not known. By means of integrated motion capture as an input to dynamic simulations, contact forces acting on the hip of 20 young healthy adults were calculated during walking and running at different speeds. During walking, hip contact forces (HCFs) have a two-peak profile whereby the first peak increases from 4.22 body weight (BW) to 5.41 BW and the second from 4.37 BW to 5.74 BW, by increasing speed from 3 to 6 km/h. During running, there is only one peak HCF that increases from 7.49 BW to 10.01 BW, by increasing speed from 6 to 12 km/h. Speed related profiles of peak HCFs and ground reaction forces (GRFs) reveal a different progression of the two peaks during walking. Speed has a stronger impact on peak HCFs rather than on peak GRFs during walking and running, suggesting an increasing influence of muscle activity on peak HCF with increased speed. Moreover, results show that the first peak of HCF during walking can be predicted best by hip adduction moment, and the second peak of HCF by hip extension moment. During running, peak HCF can be best predicted by hip adduction moment. The present study contributes hereby to a better understanding of musculoskeletal loading during walking and running in a wide range of speeds, offering valuable information to clinicians and scientists exploring bone loading as a possible nonpharmacological osteogenic stimulus. © 2015 American Society for Bone and Mineral Research.

Critical appraisal of the literature highlights that the discriminative power of gait-related features in patients with hip osteoarthritis (OA) has not been fully explored. We aimed to reduce the number of gait-related features and define the most discriminative ones comparing the three-dimensional gait analysis of 20 patients with hip osteoarthritis (OA) with those of 17 healthy peers. First, principal component analysis was used to reduce the high-dimensional gait data into a reduced set of interpretable variables for further analysis, including tests for group differences. These differences were indicative for the selection of the top 10 variables to be included into linear discriminant analysis models (LDA). Our findings demonstrated the successful data reduction of hip osteoarthritic-related gait features with a high discriminatory power. The combination of the top variables into LDA models clearly separated groups, with a maximum misclassification error rate of 19%, estimated by cross-validation. Decreased hip/knee extension, hip flexion and internal rotation moment were gait features with the highest discriminatory power. This study listed the most clinically relevant gait features characteristics of hip OA. Moreover, it will help clinicians and physiotherapists understand the movement pathomechanics related to hip OA useful in the management and design of rehabilitation intervention.

In musculoskeletal modelling, several optimization techniques are used to calculate muscle forces, which strongly influence resultant hip contact forces (HCF). The goal of this study was to calculate muscle forces using four different optimization techniques, i.e., two different static optimization techniques, computed muscle control (CMC) and the physiological inverse approach (PIA). We investigated their subsequent effects on HCFs during gait and sit to stand and found that at the first peak in gait at 15–20% of the gait cycle, CMC calculated the highest HCFs (median 3.9 times peak GRF (pGRF)). When comparing calculated HCFs to experimental HCFs reported in literature, the former were up to 238% larger. Both static optimization techniques produced lower HCFs (median 3.0 and 3.1 pGRF), while PIA included muscle dynamics without an excessive increase in HCF (median 3.2 pGRF). The increased HCFs in CMC were potentially caused by higher muscle forces resulting from co-contraction of agonists and antagonists around the hip. Alternatively, these higher HCFs may be caused by the slightly poorer tracking of the net joint moment by the muscle moments calculated by CMC. We conclude that the use of different optimization techniques affects calculated HCFs, and static optimization approached experimental values best.

Abnormal kinematics and kinetics have been reported in hip osteoarthritis (OA) patients before and after total hip arthroplasty (THA). These changes can affect the loading of the ipsilateral hip, as well as the contralateral hip and knee joint. As it is not clear how hip and knee loading evolves in THA patients during the first year after surgery, the goal of this study is to define how joint loading changes in patients before and at three evaluation times after THA surgery. Musculoskeletal modelling in combination with gait analysis data was used to calculate hip and knee contact forces in 14 patients before and 3-, 6- and 12-months after unilateral THA, as well as in 18 healthy controls. Results showed that bilateral hip and knee loading were decreased compared to controls, both before and after THA surgery. Loading symmetry was altered compared to controls at 3-months post-surgery for the hip and at all evaluation times, except for 6-months post-surgery, for the knee, with ipsilateral joint loading decreased compared to the contralateral side. To conclude, 12-months after THA joint loading was not normalized, with both hip and knee loading in patients decreased compared to controls. Therefore, no overloading of the ipsi- or contralateral hip and knee joint was found before and up to one year after unilateral THA.

Children with cerebral palsy (CP) often present aberrant hip geometry, specifically increased femoral anteversion and neck-shaft angle. Furthermore, altered gait patterns are present within this population. We analyzed the effect of aberrant femoral geometry, as present in CP subjects, on hip contact force (HCF) during pathological and normal gait. We ran dynamic simulations of CP-specific and normal gait using two musculoskeletal models (MSMs), one reflecting normal femoral geometry and one reflecting proximal femoral deformities. The combination of aberrant bone geometry and CP-specific gait characteristics reduced HCF compared to normal gait on a CP subject-specific MSM, but drastically changed the orientation of the HCF vector. The HCF was orientated more vertically and anteriorly than compared to HCF orientation during normal gait. Furthermore, subjects with more pronounced bony deformities encountered larger differences in resultant HCF and HCF orientation. When bone deformities were not accounted for in MSMs of pathologic gait, the HCF orientation was more similar to normal children. Thus, our results support a relation between aberrant femoral geometry and joint loading during pathological/normal gait and confirm a compensatory effect of altered gait kinematics on joint loading. © 2014 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 32:1406–1415, 2014.

ispartof: pages:48-49 ispartof: Gait and Posture vol:65 pages:48-49 status: published

ispartof: pages:110-111 ispartof: Gait Posture vol:65 Suppl 1 pages:110-111 ispartof: location:England status: published