Radiographic and non-invasive determination of the hip joint center location: effect on hip joint moments

Alternative Models for Pelvic Marker Occlusion in Cycling

Bike fitting aims to optimize riders' positions to improve their performance and reduce the risk of injury. To calculate joint angles, the location of the joint centers of the lower limbs needs to be identified. However, one of the greatest difficulties is the location of the hip joint center due to the frequent occlusion of the anterior superior iliac spine markers. Therefore, the objective of this study was to validate a biomechanical model adapted to cycling (modified pelvic model, MPM), based on the traditional pelvic model (TPM) with an additional lateral technical marker placed on the iliac crests. MPM was also compared with a widely used model in cycling, trochanter model (TM). Thirty-one recreational cyclists pedaled on a roller bike while the movement was captured with a 7-camera VICON system. The position of the hip joint center and knee angle were calculated and compared with the TPM continuously (along 10 pedaling cycles) and discreetly at 90° and 180° crank positions. No significant differences were found in the position of the hip joint center or in the knee flexion/extension angle between the TPM and the MPM. However, there are differences between TPM and TM (variations between 4.1° and 6.9° in favor of the TM at 90° and 180°; P < .001). Bland-Altman graphs comparing the models show an average difference or bias close to 0° (limits of agreement [0.2 to -8.5]) between TPM and MPM in both lower limbs and a mean difference of between -4° and -7° (limits of agreement [-0.6 to -13.2]) when comparing TPM and TM. Given the results, the new cycling pelvic model has proven to be valid compared with the TPM when performing bike fitting studies, with the advantage that the occluded markers are avoided. Despite its simplicity, the TM presents measurement errors that may be relevant when making diagnoses, which makes its usefulness questionable.

Conventional video recordings dependably quantify whole-body lifting strategy using the Stoop-Squat-Index: A methods comparison against motion capture and a reliability study

Whole-body lifting strategies could be derived from conventional video recordings using the Stoop-Squat-Index, which quantifies the ratio between trunk forward lean and lower extremity joint flexion from 0 (full squat) to 100 (full stoop). The purpose of this study was to compare Stoop-Squat-Indices derived from conventional video recordings to those from a three-dimensional marker-based motion capture system and to evaluate interrater and intrarater reliability of the video-based approach. Thirty healthy participants lifted a 5-kg box under different conditions (freestyle, squat, stoop). Kinematic data were recorded using a Vicon motion capture system (serving as reference standard) and an iPad camera. Stoop-Squat-Indices over the entire lifting cycle were derived separately from both approaches. Agreement was assessed using mean differences (video minus motion capture) and limits of agreement. Reliability was investigated by calculating intraclass correlation coefficients (ICC) and minimal detectable changes (MDC) over the course of the lifting cycle. Systematic errors were identified with Statistical Parametric Mapping-based T-tests. Systematic errors between the video-based and the motion capture-based approach were observed among all conditions. Mean differences in Stoop-Squat-Indices over the lifting cycle ranged from -6.9 to 3.2 (freestyle), from -1.8 to 5.3 (squat) and from -2.8 to -1.1 (stoop). Limits of agreement were lower when the box was close to the floor, and higher towards upright standing. Reliability of the video-based approach was excellent for most of the lifting cycle, with ICC above 0.995 and MDC below 3.5. These findings support using a video-based assessment of Stoop-Squat-Indices to quantify whole-body lifting strategy in field.

Development of a simulation system for femoroacetabular impingement detection based on 3D images

Image-based criteria have been adopted to diagnose femoroacetabular impingement (FAI). However, the overlapping property of the two-dimensional X-ray outlines and static and supine posture of taking computed tomography (CT) and magnetic resonance imaging images potentially affect the accuracy of the criteria. This study developed a CT image-based dynamic criterion to effectively simulate FAI, thereby providing a basis for physicians to perform pre-operative planning for arthroscopic surgery. Post-operative CT images of 20 patients with satisfactory surgical results were collected, and 10 sets of models were used to define the flexion rotation centre (FRC) of the three-dimensional FAI model. First, let these 10 groups of models simulate the FAI detection action and find the best centre offset, and then FRC is the result of averaging these 10 groups of best displacements. The model was validated in 10 additional patients. Finally, through the adjustment basis of FRC, the remaining 10 sets of models can find out the potential position of FAI during the dynamic simulation process. Rotational collisions detected using FRC indicate that the patient's post-operative flexion angle may reach 120° or greater, which is close to the actual result. The recommended surgical range of the diagnostic system (average length of 6.4 mm, width of 4.1 mm and depth of 3.2 mm) is smaller than the actual surgical results, which prevents the doctor from performing excessive resection operations, which may preserve more bones. The FRC diagnostic system detects the distribution of FAI in a simple manner. It can be used as a pre-operative diagnosis reference for clinicians, hoping to improve the effect and accuracy of debridement surgery.

AddBiomechanics: Automating model scaling, inverse kinematics, and inverse dynamics from human motion data through sequential optimization

Creating large-scale public datasets of human motion biomechanics could unlock data-driven breakthroughs in our understanding of human motion, neuromuscular diseases, and assistive devices. However, the manual effort currently required to process motion capture data and quantify the kinematics and dynamics of movement is costly and limits the collection and sharing of large-scale biomechanical datasets. We present a method, called AddBiomechanics, to automate and standardize the quantification of human movement dynamics from motion capture data. We use linear methods followed by a non-convex bilevel optimization to scale the body segments of a musculoskeletal model, register the locations of optical markers placed on an experimental subject to the markers on a musculoskeletal model, and compute body segment kinematics given trajectories of experimental markers during a motion. We then apply a linear method followed by another non-convex optimization to find body segment masses and fine tune kinematics to minimize residual forces given corresponding trajectories of ground reaction forces. The optimization approach requires approximately 3-5 minutes to determine a subject's skeleton dimensions and motion kinematics, and less than 30 minutes of computation to also determine dynamically consistent skeleton inertia properties and fine-tuned kinematics and kinetics, compared with about one day of manual work for a human expert. We used AddBiomechanics to automatically reconstruct joint angle and torque trajectories from previously published multi-activity datasets, achieving close correspondence to expert-calculated values, marker root-mean-square errors less than 2 cm, and residual force magnitudes smaller than 2% of peak external force. Finally, we confirmed that AddBiomechanics accurately reproduced joint kinematics and kinetics from synthetic walking data with low marker error and residual loads. We have published the algorithm as an open source cloud service at AddBiomechanics.org, which is available at no cost and asks that users agree to share processed and de-identified data with the community. As of this writing, hundreds of researchers have used the prototype tool to process and share about ten thousand motion files from about one thousand experimental subjects. Reducing the barriers to processing and sharing high-quality human motion biomechanics data will enable more people to use state-of-the-art biomechanical analysis, do so at lower cost, and share larger and more accurate datasets.

AddBiomechanics: Automating model scaling, inverse kinematics, and inverse dynamics from human motion data through sequential optimization

Creating large-scale public datasets of human motion biomechanics could unlock data-driven breakthroughs in our understanding of human motion, neuromuscular diseases, and assistive devices. However, the manual effort currently required to process motion capture data and quantify the kinematics and dynamics of movement is costly and limits the collection and sharing of large-scale biomechanical datasets. We present a method, called AddBiomechanics, to automate and standardize the quantification of human movement dynamics from motion capture data. We use linear methods followed by a non-convex bilevel optimization to scale the body segments of a musculoskeletal model, register the locations of optical markers placed on an experimental subject to the markers on a musculoskeletal model, and compute body segment kinematics given trajectories of experimental markers during a motion. We then apply a linear method followed by another non-convex optimization to find body segment masses and fine tune kinematics to minimize residual forces given corresponding trajectories of ground reaction forces. The optimization approach requires approximately 3-5 minutes to determine a subjecťs skeleton dimensions and motion kinematics, and less than 30 minutes of computation to also determine dynamically consistent skeleton inertia properties and fine-tuned kinematics and kinetics, compared with about one day of manual work for a human expert. We used AddBiomechanics to automatically reconstruct joint angle and torque trajectories from previously published multi-activity datasets, achieving close correspondence to expert-calculated values, marker root-mean-square errors less than 2 c m , and residual force magnitudes smaller than 2 % of peak external force. Finally, we confirmed that AddBiomechanics accurately reproduced joint kinematics and kinetics from synthetic walking data with low marker error and residual loads. We have published the algorithm as an open source cloud service at AddBiomechanics.org, which is available at no cost and asks that users agree to share processed and de-identified data with the community. As of this writing, hundreds of researchers have used the prototype tool to process and share about ten thousand motion files from about one thousand experimental subjects. Reducing the barriers to processing and sharing high-quality human motion biomechanics data will enable more people to use state-of-the-art biomechanical analysis, do so at lower cost, and share larger and more accurate datasets.

Restoring Rotation Center in Total Hip Arthroplasty for Developmental Dysplasia of the Hip with the Assistance of Three Dimensional Printing Technology: A Pilot Study

Objective

To develop a new method to restore hip rotation center exactly and rapidly in total hip arthroplasty (THA) with the assistance of three dimensional (3D) printing technology and evaluate its clinical and radiological outcomes.

Methods

From March 2014 to July 2018, a total of 17 patients (five hips of four men and 16 hips of 13 women) with end‐stage osteoarthritis secondary to developmental dysplasia of the hip who underwent THA were analyzed and followed up retrospectively. The average age is 58.00 ± 8.12 years (range from 45 to 71 years). Simulated operations were performed on 3D printed hip models for preoperative planning. The morphology of Harris fossa and acetabular notches were recognized and restored to locate the acetabular center. The size of bone defect was measured by the bone wax method. The agreement on the size of acetabular cup and bone defect between simulated operations and actual operations were analyzed. Harris Hip Score (HHS) was used to evaluate the recovery of hip joint function. The vertical distance and horizontal distance of the rotation center on the pelvis plain radiograph were measured, which were used to assess the efficacy of restoring hip rotation center and acetabular cup migration.

Results

The mean sizes of bone defect in simulated operations and THA were 4.58 ± 2.47 cm² and 4.55 ± 2.57 cm² respectively. There was no significant difference statistically between the sizes of bone defect in simulated operations and the actual sizes of bone defect in THA (t = 0.03, P = 0.97). The sizes of the acetabular cup of simulated operations on 3D print models showed a high rate of coincidence with the actual sizes in the operations (ICC = 0.93). All 17 patients were available for clinical and radiological follow‐up. The average follow‐up time was 18.35 ± 6.86 months (range, 12–36 months. The average HHS of the patients was improved from (38.33 ± 6.07) preoperatively to the last follow‐up (88.61 ± 3.44) postoperatively. The mean vertical and horizontal distances of hip rotation center on the pelvic radiographs were restored to 15.12 ± 1.25 mm and 32.49 ± 2.83 mm respectively. No case presented dislocation or radiological signs of loosening until last follow‐up.

Conclusions

The application of 3D printing technology facilitates orthopedists to recognize the morphology of Harris fossa and acetabular notches, locate the acetabular center and restore the hip rotation center rapidly and accurately.

A Normative Database of Hip and Knee Joint Biomechanics During Dynamic Tasks Using Four Functional Methods With Three Functional Calibration Tasks

Although predicted hip joint center (HJC) locations are known to vary widely between functional methods, no previous investigation has detailed functional method-dependent hip and knee biomechanics. The purpose of this study was to define a normative database of hip joint biomechanics during dynamic movements based upon functional HJC methods and calibration tasks. Thirty healthy young adults performed arc, star arc, and two-sided calibration tasks. Motion capture and ground reaction forces were collected during walking, running, and single-leg landings (SLLs). Two sphere-fit (geometric and algebraic) and two coordinate transformation techniques were implemented using each calibration (12 total method-calibration combinations). Surprisingly, the geometric fit-two-sided model placed the HJC at the midline of the pelvis and above the iliac spines, and thus was removed from analyses. A database of triplanar hip and knee kinematics and hip moments and powers was constructed using the mean of all subjects for the eleven method-calibration combinations. A nested analysis of variance approach compared calibration [method] peak hip kinematics and kinetics. Most method differences existed between geometric fit and coordinate transformations (58 of 84 total). No arc-star arc differences were found. Thirty-two differences were found between the two-sided and arc/star arc calibrations. This database of functional method based hip and knee biomechanics serves as an important reference point for interstudy comparisons. Overall, this study illustrates that functional HJC method can dramatically impact hip biomechanics and should be explicitly detailed in future work.

Can the hip joint center be estimated from pelvic dimensions in dysplastic hips?

BACKGROUND

We aimed to determine correlations between the hip joint center position and pelvic dimensions and whether the three-dimensional position of the original hip joint center could be estimated from pelvic landmarks in dysplastic and normal hips.

METHODS

We reviewed the pelvic CT scans of 70 patients (70 hips) with hip dysplasia. Seventy-seven normal hips were used as controls. The hip joint center coordinates (Cx, Cy, and Cz) and pelvic dimensions were measured with reference to the anterior pelvic plane coordinate system. Multiple regression formulas were used to estimate the original hip joint center.

RESULTS

The hip center for both dysplastic and normal hip was highly correlated with the distance between the anterior superior iliac spine (ASIS) in the coronal plane (r = 0.76 and 0.84), the distance from the ASIS to the pubic tubercle in the sagittal plane (r = 0.81 and 0.76), and distance from the pubic tubercle to the most posterior point of the ischium on the transverse plane (r = 0.76 and 0.78). The hip joint center could be estimated within a 5-mm error for more than 80% of hips on their respective axes in both dysplastic and normal hips.

CONCLUSIONS

The three-dimensional position of the original hip joint center was correlated with pelvic dimensions, and can be estimated with substantial accuracy using pelvic landmarks as references. Although these results are preliminary, this estimation method may be useful for surgeons planning total hip arthroplasties.

Methodological factors affecting joint moments estimation in clinical gait analysis: a systematic review

Quantitative gait analysis can provide a description of joint kinematics and dynamics, and it is recognized as a clinically useful tool for functional assessment, diagnosis and intervention planning. Clinically interpretable parameters are estimated from quantitative measures (i.e. ground reaction forces, skin marker trajectories, etc.) through biomechanical modelling. In particular, the estimation of joint moments during motion is grounded on several modelling assumptions: (1) body segmental and joint kinematics is derived from the trajectories of markers and by modelling the human body as a kinematic chain; (2) joint resultant (net) loads are, usually, derived from force plate measurements through a model of segmental dynamics. Therefore, both measurement errors and modelling assumptions can affect the results, to an extent that also depends on the characteristics of the motor task analysed (i.e. gait speed). Errors affecting the trajectories of joint centres, the orientation of joint functional axes, the joint angular velocities, the accuracy of inertial parameters and force measurements (concurring to the definition of the dynamic model), can weigh differently in the estimation of clinically interpretable joint moments. Numerous studies addressed all these methodological aspects separately, but a critical analysis of how these aspects may affect the clinical interpretation of joint dynamics is still missing. This article aims at filling this gap through a systematic review of the literature, conducted on Web of Science, Scopus and PubMed. The final objective is hence to provide clear take-home messages to guide laboratories in the estimation of joint moments for the clinical practice.

Prediction of Polyethylene Wear Rates from Gait Biomechanics and Implant Positioning in Total Hip Replacement

BACKGROUND

Patient-specific gait and surgical variables are known to play an important role in wear of total hip replacements (THR). However a rigorous model, capable of predicting wear rate based on a comprehensive set of subject-specific gait and component-positioning variables, has to our knowledge, not been reported.

QUESTIONS/PURPOSE

(1) Are there any differences between patients with high, moderate, and low wear rate in terms of gait and/or positioning variables? (2) Can we design a model to predict the wear rate based on gait and positioning variables? (3) Which group of wear factors (gait or positioning) contributes more to the wear rate?

PATIENTS AND METHODS

Data on patients undergoing primary unilateral THR who performed a postoperative gait test were screened for inclusion. We included patients with a 28-mm metal head and a hip cup made of noncrosslinked polyethylene (GUR 415 and 1050) from a single manufacturer (Zimmer, Inc). To calculate wear rates from radiographs, inclusion called for patients with a series of standing radiographs taken more than 1 year after surgery. Further, exclusion criteria were established to obtain reasonably reliable and homogeneous wear readings. Seventy-three (83% of included) patients met all criteria, and the final dataset consisted of 43 males and 30 females, 69 ± 10 years old, with a BMI of 27.3 ± 4.7 kg/m². Wear rates of these patients were determined based on the relative displacement of the femoral head with regard to the cup using a validated computer-assisted X-ray wear-analysis suite. Three groups with low (< 0.1 mm/year), moderate (0.1 to 0.2 mm/year), and high (> 0.2 mm/year) wear were established. Wear prediction followed a two-step process: (1) linear discriminant analysis to estimate the level of wear (low, moderate, or high), and (2) multiple linear and nonlinear regression modeling to predict the exact wear rate from gait and implant-positioning variables for each level of wear.

RESULTS

There were no group differences for positioning and gait suggesting that wear differences are caused by a combination of wear factors rather than single variables. The linear discriminant analysis model correctly predicted the level of wear in 80% of patients with low wear, 87% of subjects with moderate wear, and 73% of subjects with high wear based on a combination of gait and positioning variables. For every wear level, multiple linear and nonlinear regression showed strong associations between gait biomechanics, implant positioning, and wear rate, with the nonlinear model having a higher prediction accuracy. Flexion-extension ROM and hip moments in the sagittal and transverse planes explained 42% to 60% of wear rate while positioning factors, (such as cup medialization and cup inclination angle) explained only 10% to 33%.

CONCLUSION

Patient-specific wear rates are associated with patients' gait patterns. Gait pattern has a greater influence on wear than component positioning for traditional metal-on-polyethylene bearings.

CLINICAL RELEVANCE

The consideration of individual gait bears potential to further reduce implant wear in THR. In the future, a predictive wear model may identify individual, modifiable wear factors for modern materials.

Bilateral, Misalignment-Compensating, Full-DOF Hip Exoskeleton: Design and Kinematic Validation

A shared design goal for most robotic lower limb exoskeletons is to reduce the metabolic cost of locomotion for the user. Despite this, only a limited amount of devices was able to actually reduce user metabolic consumption. Preservation of the natural motion kinematics was defined as an important requirement for a device to be metabolically beneficial. This requires the inclusion of all human degrees of freedom (DOF) in a design, as well as perfect alignment of the rotation axes. As perfect alignment is impossible, compensation for misalignment effects should be provided. A misalignment compensation mechanism for a 3-DOF system is presented in this paper. It is validated by the implementation in a bilateral hip exoskeleton, resulting in a compact and lightweight device that can be donned fast and autonomously, with a minimum of required adaptations. Extensive testing of the prototype has shown that hip range of motion of the user is maintained while wearing the device and this for all three hip DOFs. This allowed the users to maintain their natural motion patterns when they are walking with the novel hip exoskeleton.

In-vivo quantification of dynamic hip joint center errors and soft tissue artifact

Hip joint center (HJC) measurement error can adversely affect predictions from biomechanical models. Soft tissue artifact (STA) may exacerbate HJC errors during dynamic motions. We quantified HJC error and the effect of STA in 11 young, asymptomatic adults during six activities. Subjects were imaged simultaneously with reflective skin markers (SM) and dual fluoroscopy (DF), an x-ray based technique with submillimeter accuracy that does not suffer from STA. Five HJCs were defined from locations of SM using three predictive (i.e., based on regression) and two functional methods; these calculations were repeated using the DF solutions. Hip joint center motion was analyzed during six degrees-of-freedom (default) and three degrees-of-freedom hip joint kinematics. The position of the DF-measured femoral head center (FHC), served as the reference to calculate HJC error. The effect of STA was quantified with mean absolute deviation. HJC errors were (mean±SD) 16.6±8.4mm and 11.7±11.0mm using SM and DF solutions, respectively. HJC errors from SM measurements were all significantly different from the FHC in at least one anatomical direction during multiple activities. The mean absolute deviation of SM-based HJCs was 2.8±0.7mm, which was greater than that for the FHC (0.6±0.1mm), suggesting that STA caused approximately 2.2mm of spurious HJC motion. Constraining the hip joint to three degrees-of-freedom led to approximately 3.1mm of spurious HJC motion. Our results indicate that STA-induced motion of the HJC contributes to the overall error, but inaccuracies inherent with predictive and functional methods appear to be a larger source of error.

Sex-based differences in lifting technique under increasing load conditions: A principal component analysis

The objective of the present study was to determine if there is a sex-based difference in lifting technique across increasing-load conditions. Eleven male and 14 female participants (n = 25) with no previous history of low back disorder participated in the study. Participants completed freestyle, symmetric lifts of a box with handles from the floor to a table positioned at 50% of their height for five trials under three load conditions (10%, 20%, and 30% of their individual maximum isometric back strength). Joint kinematic data for the ankle, knee, hip, and lumbar and thoracic spine were collected using a two-camera Optotrak motion capture system. Joint angles were calculated using a three-dimensional Euler rotation sequence. Principal component analysis (PCA) and single component reconstruction were applied to assess differences in lifting technique across the entire waveforms. Thirty-two PCs were retained from the five joints and three axes in accordance with the 90% trace criterion. Repeated-measures ANOVA with a mixed design revealed no significant effect of sex for any of the PCs. This is contrary to previous research that used discrete points on the lifting curve to analyze sex-based differences, but agrees with more recent research using more complex analysis techniques. There was a significant effect of load on lifting technique for five PCs of the lower limb (PC1 of ankle flexion, knee flexion, and knee adduction, as well as PC2 and PC3 of hip flexion) (p < 0.005). However, there was no significant effect of load on the thoracic and lumbar spine. It was concluded that when load is standardized to individual back strength characteristics, males and females adopted a similar lifting technique. In addition, as load increased male and female participants changed their lifting technique in a similar manner.

Evaluation of Anatomical and Functional Hip Joint Center Methods: The Effects of Activity Type, Gender, and Proximal Reference Segment

Accurate hip joint center (HJC) location is critical when studying hip joint biomechanics. The HJC is often determined from anatomical methods, but functional methods are becoming increasingly popular. Several studies have examined these methods using simulations and in vivo gait data, but none has studied high-range of motion activities, such a chair rise, nor has HJC prediction been compared between males and females. Furthermore, anterior superior iliac spine (ASIS) marker visibility during chair rise can be problematic, requiring a sacral cluster as an alternative proximal segment; but functional HJC has not been explored using this approach. For this study, the quality of HJC measurement was based on the joint gap error (JGE), which is the difference in global HJC between proximal and distal reference segments. The aims of the present study were to: (1) determine if JGE varies between pelvic and sacral referenced HJC for functional and anatomical methods, (2) investigate which functional calibration motion results in the lowest JGE and if the JGE varies depending on movement type (gait versus chair rise) and gender, and (3) assess whether the functional HJC calibration results in lower JGE than commonly used anatomical approaches and if it varies with movement type and gender. Data were collected on 39 healthy adults (19 males and 20 females) aged 14–50 yr old. Participants performed four hip “calibration” tests (arc, cross, star, and star-arc), as well as gait and chair rise (activities of daily living (ADL)). Two common anatomical methods were used to estimate HJC and were compared to HJC computed using a published functional method with the calibration motions above, when using pelvis or sacral cluster as the proximal reference. For ADL trials, functional methods resulted in lower JGE (12–19 mm) compared to anatomical methods (13–34 mm). It was also found that women had significantly higher JGE compared to men and JGE was significantly higher for chair rise compared to gait, across all methods. JGE for sacrum referenced HJC was consistently higher than for the pelvis, but only by 2.5 mm. The results indicate that dynamic hip range of movement and gender are significant factors in HJC quality. The findings also suggest that a rigid sacral cluster for HJC estimation is an acceptable alternative for relying solely on traditional pelvis markers.

Accuracy evaluation of a new stereophotogrammetry-based functional method for joint kinematic analysis in biomechanics

The human joint kinematics is an interesting topic in biomechanics and turns to be useful for the analysis of human movement in several fields. A crucial issue regards the assessment of joint parameters, like axes and centers of rotation, due to the direct influence on human motion patterns. A proper accuracy in the estimation of these parameters is hence required. On the whole, stereophotogrammetry-based predictive methods and, as an alternative, functional ones can be used to this end. This article presents a new functional algorithm for the assessment of knee joint parameters, based on a polycentric hinge model for the knee flexion-extension. The proposed algorithm is discussed, identifying its fields of application and its limits. The techniques for estimating the joint parameters from the metrological point of view are analyzed, so as to lay the groundwork for enhancing and eventually replacing predictive methods, currently used in the laboratories of human movement analysis. This article also presents an assessment of the accuracy associated with the whole process of measurement and joint parameters estimation. To this end, the presented functional method is tested through both computer simulations and a series of experimental laboratory tests in which swing motions were imposed to a polycentric mechanical analogue and a stereophotogrammetric system was used to record them.

Estimation of the hip joint centre in human motion analysis: a systematic review

BACKGROUND

Inaccuracies in locating the three-dimensional position of the hip joint centre affect the calculated hip and knee kinematics, force- and moment-generating capacity of muscles and hip joint mechanics, which can lead to incorrect interpretations and recommendations in gait analysis. Several functional and predictive methods have been developed to estimate the hip joint centre location, and the International Society of Biomechanics recommends a functional approach for use with participants that have adequate range of motion at the hip, and predictive methods in those with insufficient range of motion. The purpose of the current systematic review was to substantiate the International Society of Biomechanics recommendations. This included identifying the most accurate functional and predictive methods, and defining 'adequate' range of motion.

METHODS

A systematic search with broad search terms was performed including five databases.

FINDINGS

The systematic search yielded to 801 articles, of which 34 papers were included. Eleven different predictive and 13 different functional methods were identified. The results showed that the geometric sphere fit method and Harrington equations are the most accurate functional and predictive approaches respectively that have been evaluated in vivo.

INTERPRETATION

In regard to the International Society of Biomechanics recommendations, the geometric sphere fit method should be used in people with sufficient active hip range of motion and the Harrington equations should be used in patients without sufficient hip range of motion. Multi-plane movement trials with at least 60° of flexion-extension and 30° of ab-adduction range of motion are suggested when using functional methods.

The clinical impact of hip joint centre regression equation error on kinematics and kinetics during paediatric gait

Regression equations based on pelvic anatomy are routinely used to estimate the hip joint centre during gait analysis. While the associated errors have been well documented, the clinical significance of these errors has not been reported. This study investigated the clinical agreement of three commonly used regression equation sets (Bell et al., Davis et al. and Orthotrak software) against the equations of Harrington et al. Full 3-dimensional gait analysis was performed on 18 healthy paediatric subjects. Kinematic and kinetic data were calculated using each set of regression equations and compared to Harrington et al. In addition, the Gait Profile Score and GDI-Kinetic were used to assess clinical significance. Bell et al. was the best performing set with differences in Gait Profile Score (0.13°) and GDI-Kinetic (0.84 points) falling below the clinical significance threshold. Small deviations were present for the Orthotrak set for hip abduction moment (0.1 Nm/kg), however differences in Gait Profile Score (0.27°) and GDI-Kinetic (2.26 points) remained below the clinical threshold. Davis et al. showed least agreement with a clinically significant difference in GDI-Kinetic score (4.36 points). It is proposed that Harrington et al. or Bell et al. regression equation sets are used during gait analysis especially where inverse dynamic data are calculated. Orthotrak is a clinically acceptable alternative however clinicians must be aware of the effects of error on hip abduction moment. The Davis et al. set should be used with caution for inverse dynamic analysis as error could be considered clinically meaningful.

Displacement of the hip center of rotation after arthroplasty of Crowe III and IV dysplasia: a radiological and biomechanical study

To study the direction and biomechanical consequences of hip center of rotation (HCOR) migration in Crowe type III and VI hips after total hip arthroplasty, post-operative radiographs and CT scans of several unilaterally affected hips were evaluated. Using a three-dimensional model of the human hip, the HCOR was moved in all directions, and joint reaction force (JRF) and abductor muscle force (AMF) were calculated for single-leg stance configuration. Comparing to the normal side, HCOR had displaced medially and inferiorly by an average of 23.4% and 20.8%, respectively, of the normal femoral head diameter. Significant decreases in JRF (13%) and AMF (46.13%) were observed in a presumptive case with that amount of displacement. Isolated inferior displacement had a small, increasing effect on these forces. In Crowe type III and IV hips, the HCOR migrates inferiorly and medially after THA, resulting in a decrease in JRF, AMF, and abductor muscle contraction force.

Hip joint center localisation: A biomechanical application to hip arthroplasty population

AIM: To determine hip joint center (HJC) location on hip arthroplasty population comparing predictive and functional approaches with radiographic measurements.

METHODS: The distance between the HJC and the mid-pelvis was calculated and compared between the three approaches. The localisation error between the predictive and functional approach was compared using the radiographic measurements as the reference. The operated leg was compared to the non-operated leg.

RESULTS: A significant difference was found for the distance between the HJC and the mid-pelvis when comparing the predictive and functional method. The functional method leads to fewer errors. A statistical difference was found for the localization error between the predictive and functional method. The functional method is twice more precise.

CONCLUSION: Although being more individualized, the functional method improves HJC localization and should be used in three-dimensional gait analysis.

A QUASI-LINEAR VISCOELASTIC MODEL FOR THE PASSIVE PROPERTIES OF THE HUMAN HIP JOINT

Properties of passive elastic structures constituting the human hip joint can be exploited to increase efficiency of human locomotion. As studies estimating the passive contributions to the net joint moment often disregard damping properties of the joint such contributions overestimate the energy gained during leg retraction within swing and stance phase. We built an experimental apparatus to measure moment-angle-relations during motor guided cyclic movements over a wide range of angular velocities and step-like changes in hip angle. On the basis of the experimentally gained data set the objective of this study was to model the elastic as well as the damping characteristics of the joint in the sagittal plane utilizing the Quasi-Linear Viscoelastic theory (QLV). A double exponential function was conveniently employed to describe the elastic response. The dependency of the hip joint stiffness on biarticular muscles was incorporated by repeating the measurement protocol for different knee angles. Due to the fact that the stiffness characteristics of the elastic response were merely shifted over knee angles we introduced an equilibrium angle at the hip joint as exponential function of the knee angle eventually yielding an elastic response as a function of hip and knee angle. In order to cover the damping characteristics the reduced relaxation function comprising a continuous spectrum of relaxation was utilized. We exemplify the applicability of the QLV model on published kinematic data on human walking and estimated that approximately 27% of the energy passively stored at the hip dissipates during the gait cycle.

The personal lift-assist device and lifting technique: a principal component analysis

The personal lift-assist device (PLAD) is a non-motorised, on-body device that acts as an external force generator using the concept of stored elastic energy. In this study, the effect of the PLAD on the lifting kinematics of male and female lifters was investigated using principal component analysis. Joint kinematic data of 15 males and 15 females were collected using an opto-electronic system during a freestyle, symmetrical-lifting protocol with and without wearing the PLAD. Of the 31 Principal Components (PCs) retained in the models, eight scores were significantly different between the PLAD and no-PLAD conditions. There were no main effects for gender and no significant interactions. Results indicated that the PLAD similarly affected the lifting kinematics of males and females; demonstrating significantly less lumbar and thoracic flexion and significantly greater hip and ankle flexion when wearing the PLAD. These findings add to the body of work that suggest the PLAD may be a safe and effective ergonomic aid. STATEMENT OF RELEVANCE: The PLAD is an ergonomic aid that has been shown to be effective at reducing low back demands during manual materials handling tasks. This body of work establishes that the PLAD encourages safe lifting practices without adversely affecting lifting technique.

Effect of frontal plane tibiofemoral angle on the stress and strain at the knee cartilage during the stance phase of gait

Subject-specific three-dimensional finite element models of the knee joint were created and used to study the effect of the frontal plane tibiofemoral angle on the stress and strain distribution in the knee cartilage during the stance phase of the gait cycle. Knee models of three subjects with different tibiofemoral angle and body weight were created based on magnetic resonance imaging of the knee. Loading and boundary conditions were determined from motion analysis and force platform data, in conjunction with the muscle-force reduction method. During the stance phase of walking, all subjects exhibited a valgus-varus-valgus knee moment pattern with the maximum compressive load and varus knee moment occurring at approximately 25% of the stance phase of the gait cycle. Our results demonstrated that the subject with varus alignment had the largest stresses at the medial compartment of the knee compared to the subjects with normal alignment and valgus alignment, suggesting that this subject might be most susceptible to developing medial compartment osteoarthritis (OA). In addition, the magnitude of stress and strain on the lateral cartilage of the subject with valgus alignment were found to be larger compared to subjects with normal alignment and varus alignment, suggesting that this subject might be most susceptible to developing lateral compartment knee OA.

Validation of a Morphometric Reconstruction Technique Applied to a Juvenile Pelvis

Three-dimensional reconstructions of bone geometry from microCT (computed tomography) data are frequently used in biomechanical and finite element analyses. Digitization of bone models is usually a simple process for specimens with a complete geometry, but in instances of damage or disarticulation it can be very challenging. Subsequent to digitization, further imaging techniques are often required to estimate the geometry of missing bone or connecting cartilage. This paper presents an innovative approach to the reconstruction of incomplete scan data, to reproduce proper anatomical arrangements of bones, including absent connecting cartilaginous elements. Utilizing geometric morphometric tools, the reconstruction technique is validated through comparison of a reconstructed 9 year old pelvis, to the original CT data. A principal component analysis and an overlay of the two pelves provide a measure of the accuracy of the reconstructed model. Future work aims to investigate the biomechanical effects of any minor positional error on the bone's predicted structural properties through the use of finite element analysis.

A method to calculate the centre of the ankle joint: a comparison with the Vicon Plug-in-Gait model

BACKGROUND

In gait analysis, calculation of the ankle joint centre is a difficult task. The conventional way to calculate the ankle joint centre is using the Vicon Plug-in-Gait model. The present study proposes a new model, which calculates the joint centre from two markers positioned over the medial and lateral malleoli (i.e. Two-marker-model).

METHODS

In order to compare the proposed model with Plug-in-Gait model, gait data from healthy and patient subjects were captured using a motion capture system. The ankle joint centres were calculated by the two models. A test-retest experiment was carried out to check reliability and repeatability for Two-marker-model.

FINDINGS

Two ankle joint centres produced by two models were significantly different. The distances between two ankle joint centres were approximately 16.8 (mm), and the differences in the posterior-anterior, medial-lateral and inferior-superior directions were approximately 6.3, 7.7 and 8.2 (mm). Further error analysis highlighted that the probability of producing errors in Two-marker-model is lower than that in Plug-in-Gait model due to the Two-marker-model's simple and reliable marker positioning. The reliability and repeatability coefficients for the new model were greater than 0.9.

INTERPRETATION

In principle, the Plug-in-Gait model is more likely to produce errors than the Two-marker-model, because the former employs multiple markers from the pelvis to calf to define the ankle joint centre with marker positions being very user-dependent. The results suggest that the Two-marker-model can be considered an alternative to Plug-in-Gait model for calculating ankle joint centre.

Marker-Based Reconstruction of the Kinematics of a Chain of Segments: A New Method That Incorporates Joint Kinematic Constraints

The aim of skin-marker-based motion analysis is to reconstruct the motion of a kinematical model from noisy measured motion of skin markers. Existing kinematic models for reconstruction of chains of segments can be divided into two categories: analytical methods that do not take joint constraints into account and numerical global optimization methods that do take joint constraints into account but require numerical optimization of a large number of degrees of freedom, especially when the number of segments increases. In this study, a new and largely analytical method for a chain of rigid bodies is presented, interconnected in spherical joints (chain-method). In this method, the number of generalized coordinates to be determined through numerical optimization is three, irrespective of the number of segments. This new method is compared with the analytical method of Veldpaus et al. [1988, “A Least-Squares Algorithm for the Equiform Transformation From Spatial Marker Co-Ordinates,” J. Biomech., 21, pp. 45–54] (Veldpaus-method, a method of the first category) and the numerical global optimization method of Lu and O’Connor [1999, “Bone Position Estimation From Skin-Marker Co-Ordinates Using Global Optimization With Joint Constraints,” J. Biomech., 32, pp. 129–134] (Lu-method, a method of the second category) regarding the effects of continuous noise simulating skin movement artifacts and regarding systematic errors in joint constraints. The study is based on simulated data to allow a comparison of the results of the different algorithms with true (noise- and error-free) marker locations. Results indicate a clear trend that accuracy for the chain-method is higher than the Veldpaus-method and similar to the Lu-method. Because large parts of the equations in the chain-method can be solved analytically, the speed of convergence in this method is substantially higher than in the Lu-method. With only three segments, the average number of required iterations with the chain-method is 3.0±0.2 times lower than with the Lu-method when skin movement artifacts are simulated by applying a continuous noise model. When simulating systematic errors in joint constraints, the number of iterations for the chain-method was almost a factor 5 lower than the number of iterations for the Lu-method. However, the Lu-method performs slightly better than the chain-method. The RMSD value between the reconstructed and actual marker positions is approximately 57% of the systematic error on the joint center positions for the Lu-method compared with 59% for the chain-method.

Evaluation of formal methods in hip joint center assessment: an in vitro analysis

BACKGROUND

The hip joint center is a fundamental landmark in the identification of lower limb mechanical axis; errors in its location lead to substantial inaccuracies both in joint reconstruction and in gait analysis. Actually in Computer Aided Surgery functional non-invasive procedures have been tested in identifying this landmark, but an anatomical validation is scarcely discussed.

METHODS

A navigation system was used to acquire data on eight cadaveric hips. Pivoting functional maneuver and hip joint anatomy were analyzed. Two functional methods - both with and without using the pelvic tracker - were evaluated: specifically a sphere fit method and a transformation techniques. The positions of the estimated centers with respect to the anatomical center of the femoral head, the influence of this deviation on the kinematic assessment and on the identification of femoral mechanical axis were analyzed.

FINDINGS

We found that the implemented transformation technique was the most reliable estimation of hip joint center, introducing a - Mean (SD) - difference of 1.6 (2.7) mm from the anatomical center with the pelvic tracker, whereas sphere fit method without it demonstrated the lowest accuracy with 25.2 (18.9) mm of deviation. Otherwise both the methods reported similar accuracy (<3mm of deviation).

INTERPRETATION

The functional estimations resulted in the best case to be in an average of less than 2mm from the anatomical center, which corresponds to angular deviations of the femoral mechanical axis smaller than 1.7 (1.3) degrees and negligible errors in kinematic assessment of angular displacements.

Gait and stair function in total and resurfacing hip arthroplasty: a pilot study

Standard total hip arthroplasty (THA) is the established surgical treatment for patients older than 65 years with progressive osteoarthritis but survivorship curves wane in patients younger than 50. Resurfacing hip arthroplasty (RHA) is an alternative for younger, active patients reportedly providing superior range of motion. Quantitative investigation of functional recovery following arthroplasty may elucidate limitations that aid in device selection. Although limited long-term kinematic data are available, the early rate of recovery and gait compensations are not well described. This information may aid in refining rehabilitation protocols based on limitations specific to the implant. We presumed hip motion and forces for subjects receiving RHA are more similar to age-matched controls during physically demanding tasks, such as stair negotiation, at early time points than those for THA. In a pilot study, we quantified walking and stair negotiation preoperatively and 3 months postoperatively for seven patients with RHA (mean age, 49 years), seven patients with standard THA (mean age, 52 years), and seven age-matched control subjects (mean age, 56 years). Although both treatment groups demonstrated trends toward functional recovery, the RHA group had greater improvements in hip extension and abduction moment indicating typical loading of the hip. Further investigation is needed to determine if differences persist long term or are clinically meaningful.

Subject-specific hip geometry and hip joint centre location affects calculated contact forces at the hip during gait

Hip loading affects the development of hip osteoarthritis, bone remodelling and osseointegration of implants. In this study, we analyzed the effect of subject-specific modelling of hip geometry and hip joint centre (HJC) location on the quantification of hip joint moments, muscle moments and hip contact forces during gait, using musculoskeletal modelling, inverse dynamic analysis and static optimization. For 10 subjects, hip joint moments, muscle moments and hip loading in terms of magnitude and orientation were quantified using three different model types, each including a different amount of subject-specific detail: (1) a generic scaled musculoskeletal model, (2) a generic scaled musculoskeletal model with subject-specific hip geometry (femoral anteversion, neck-length and neck-shaft angle) and (3) a generic scaled musculoskeletal model with subject-specific hip geometry including HJC location. Subject-specific geometry and HJC location were derived from CT. Significant differences were found between the three model types in HJC location, hip flexion-extension moment and inclination angle of the total contact force in the frontal plane. No model agreement was found between the three model types for the calculation of contact forces in terms of magnitude and orientations, and muscle moments. Therefore, we suggest that personalized models with individualized hip joint geometry and HJC location should be used for the quantification of hip loading. For biomechanical analyses aiming to understand modified hip joint loading, and planning hip surgery in patients with osteoarthritis, the amount of subject-specific detail, related to bone geometry and joint centre location in the musculoskeletal models used, needs to be considered.

Can the acetabular position be derived from a pelvic frame of reference?

Acetabular center positioning has an effect on hip function. However, reported clinical and plain radiographic methods are inaccurate and unreliable for ascertaining acetabular implant location. In an exploratory study we asked whether the normal acetabular position can be derived from simple radiographically measurable pelvic dimensions. We analyzed computed tomographic scans of 37 normal hips using a pelvic frame of reference centered on the ipsilateral anterior-superior iliac spine. We defined the x-, y-, and z-coordinates of the hip center (C(x),C(y),C(z)) as a percentage of the corresponding pelvic dimensions (D(x),D(y),D(z)). C(x)/D(x) averaged 9%, C(y)/D(y) 34%, and C(z)/D(z) 37%. These ratios had narrow distributions with small confidence intervals. Interobserver agreement tests showed a mean intraclass correlation coefficient of 0.95. We observed gender differences in the ratios of as much as 4%, which correspond to differences of as much as 9 mm in the hip center position. The ratios provide a simple and reliable way of deriving the normal position of the hip center from the pelvic dimensions alone. This gives the surgeon a simple way of planning where the hip center should be and may be particularly helpful in revision hip arthroplasty or in cases involving extensive osteophytes, dysplasia, or protrusio.

Frontal knee alignment: three-dimensional marker positions and clinical assessment

We assessed the validity of the hip-knee-ankle angle measured statically during three-dimensional (3-D) gait analysis and the tibial angle using an inclinometer compared with the mechanical axis on radiographs. Eleven individuals (20 knees) with radiographic knee osteoarthritis (OA) participated in this study. We determined the following: the lower-limb mechanical axis using weightbearing long-leg radiographs; hip-knee-ankle angle using the techniques of 3-D gait analysis in a static standing position; and tibial alignment using an inclinometer. The mean mechanical axis (+/- standard deviation) for this cohort was 0.7 degrees +/- 7.2 degrees (range, -13 degrees-16 degrees). The tibial alignment and hip-knee-ankle angle correlated with the mechanical axis but the correlation between the mechanical axis and the hip-knee-ankle angle was stronger. Our data suggest the inclinometer and 3-D gait analysis are valid ways to estimate mechanical alignment of the knee.

A simple, anatomically based correction to the conventional ankle joint center

BACKGROUND

Conventional motion analysis studies define the ankle joint center as the midpoint between the most medial and lateral aspects of the malleoli, yet research points toward a more distal joint center location. The purpose of this study was to develop and evaluate an anatomically based correction that would move the conventional ankle joint center to a more accurate location.

METHODS

Lower extremity radiographs from 30 pediatric patients were analyzed retrospectively. An offset between the conventional and more accurate ankle joint centers was measured and correlated to other common anatomical measures based on conventional skin mounted marker positions. The best correlated measure was used to define a simple correction factor, which was subsequently evaluated by its effect on six degree-of-freedom ankle joint translations during normal gait (n=8).

FINDINGS

Shank length was found to have the highest bivariate linear correlation (r=0.89) with the offset. Adjusting the ankle joint center using a percentage of shank length (2.7%) was also as accurate as the regression equation in predicting offset (mean error 0.6mm, or 6% offset). Adjusting the ankle joint center using this simple percentage resulted in a 25% reduction in mean ankle joint translations during normal gait.

INTERPRETATION

The accuracy of the ankle joint center can be increased through a simple, anatomically based correction. This correction may prove beneficial in some kinematic and kinetic applications requiring increased anatomical fidelity.

Mechanical factors relate to pain in knee osteoarthritis

BACKGROUND

Pain experienced by people with knee osteoarthritis is related to psychosocial factors and damage to articular tissues and/or the pain pathway itself. Mechanical factors have been speculated to trigger this pain experience; yet mechanics have not been identified as a source of pain in this population. The purpose of this study was to identify whether mechanics could explain variance in pain intensity in people with knee osteoarthritis.

METHODS

Data from 53 participants with physician-diagnosed knee osteoarthritis (mean age=68.5 years; standard deviation=8.6 years) were analyzed. Pain intensity was reported on the Western Ontario and McMaster Universities Osteoarthritis Index. Mechanical measures included weight-bearing varus-valgus alignment, body mass index and isokinetic quadriceps torque. Gait analysis captured the range of adduction-abduction angle, range of flexion-extension angle and external knee adduction moment during level walking.

FINDINGS

Pain intensity was significantly related to the dynamic range of flexion-extension during gait and body mass index. A total of 29% of the variance in pain intensity was explained by mechanical variables. The range of flexion-extension explained 18% of variance in pain intensity. Body mass index added 11% to the model. The knee adduction moment was unrelated to pain intensity.

INTERPRETATION

The findings support that mechanical factors are related to knee osteoarthritis pain. Because limitations in flexion-extension range of motion and body size are modifiable factors, future research could examine whether interventions targeting these mechanics would facilitate pain management.

Hip joint center location by fitting conchoid shape to the acetabular rim region of MR images

This paper proposes a hip joint center (HJC) location method based on the approximating of the acetabulum with a conchoid shape to acetabular rim in MR images. As the human hip joint is not a perfect sphere but it is close to a conchoid shape, the accurate location of the HJC cannot be computed as the center of a sphere. By approximating the acetabulum with a conchoid shape, it is possible to compensate a possible hip joint center location error due to the inaccuracy of 3D surface models for using functional method to calculate HJC location. For the hip joint surgery planning, it is necessary to assess patients' hip range of motions based on the HJC, thus the accurate HJC location is important.

A framework for the functional identification of joint centers using markerless motion capture, validation for the hip joint

The objective of the study was to develop a framework for the accurate identification of joint centers to be used for the calculation of human body kinematics and kinetics. The present work introduces a method for the functional identification of joint centers using markerless motion capture (MMC). The MMC system used 8 color VGA cameras. An automatic segmentation-registration algorithm was developed to identify the optimal joint center in a least-square sense. The method was applied to the hip joint center with a validation study conducted in a virtual environment. The results had an accuracy (6mm mean absolute error) below the current MMC system resolution (1cm voxel resolution). Direct experimental comparison with marker-based methods was carried out showing mean absolute deviations over the three anatomical directions of 11.9 and 15.3mm if compared with either a full leg or only thigh markers protocol, respectively. Those experimental results were presented only in terms of deviations between the two systems (marker-based and markerless) as no real gold standard was available. The methods presented in this paper provide an important enabling step towards the biomechanical and clinical applications of markerless motion capture.