Roger Zauel

Validation of a new model-based tracking technique for measuring three-dimensional, in vivo glenohumeral joint kinematics.

Shoulder motion is complex and significant research efforts have focused on measuring glenohumeral joint motion. Unfortunately, conventional motion measurement techniques are unable to measure glenohumeral joint kinematics during dynamic shoulder motion to clinically significant levels of accuracy. The purpose of this study was to validate the accuracy of a new model-based tracking technique for measuring three-dimensional, in vivo glenohumeral joint kinematics. We have developed a model-based tracking technique for accurately measuring in vivo joint motion from biplane radiographic images that tracks the position of bones based on their three-dimensional shape and texture. To validate this technique, we implanted tantalum beads into the humerus and scapula of both shoulders from three cadaver specimens and then recorded biplane radiographic images of the shoulder while manually moving each specimens arm. The position of the humerus and scapula were measured using the model-based tracking system and with a previously validated dynamic radiostereometric analysis (RSA) technique. Accuracy was reported in terms of measurement bias, measurement precision, and overall dynamic accuracy by comparing the model-based tracking results to the dynamic RSA results. The model-based tracking technique produced results that were in excellent agreement with the RSA technique. Measurement bias ranged from -0.126 to 0.199 mm for the scapula and ranged from -0.022 to 0.079 mm for the humerus. Dynamic measurement precision was better than 0.130 mm for the scapula and 0.095 mm for the humerus. Overall dynamic accuracy indicated that rms errors in any one direction were less than 0.385 mm for the scapula and less than 0.374 mm for the humerus. These errors correspond to rotational inaccuracies of approximately 0.25 deg for the scapula and 0.47 deg for the humerus. This new model-based tracking approach represents a non-invasive technique for accurately measuring dynamic glenohumeral joint motion under in vivo conditions. The model-based technique achieves accuracy levels that far surpass all previously reported non-invasive techniques for measuring in vivo glenohumeral joint motion. This technique is supported by a rigorous validation study that provides a realistic simulation of in vivo conditions and we fully expect to achieve these levels of accuracy with in vivo human testing. Future research will use this technique to analyze shoulder motion under a variety of testing conditions and to investigate the effects of conservative and surgical treatment of rotator cuff tears on dynamic joint stability.

Accuracy of biplane x-ray imaging combined with model-based tracking for measuring in-vivo patellofemoral joint motion.

BackgroundAccurately measuring in-vivo motion of the knees patellofemoral (PF) joint is challenging. Conventional measurement techniques have largely been unable to accurately measure three-dimensional, in-vivo motion of the patella during dynamic activities. The purpose of this study was to assess the accuracy of a new model-based technique for measuring PF joint motion.MethodsTo assess the accuracy of this technique, we implanted tantalum beads into the femur and patella of three cadaveric knee specimens and then recorded dynamic biplane radiographic images while manually flexing and extending the specimen. The position of the femur and patella were measured from the biplane images using both the model-based tracking system and a validated dynamic radiostereometric analysis (RSA) technique. Model-based tracking was compared to dynamic RSA by computing measures of bias, precision, and overall dynamic accuracy of four clinically-relevant kinematic parameters (patellar shift, flexion, tilt, and rotation).ResultsThe model-based tracking technique results were in excellent agreement with the RSA technique. Overall dynamic accuracy indicated errors of less than 0.395 mm for patellar shift, 0.875° for flexion, 0.863° for tilt, and 0.877° for rotation.ConclusionThis model-based tracking technique is a non-invasive method for accurately measuring dynamic PF joint motion under in-vivo conditions. The technique is sufficiently accurate in measuring clinically relevant changes in PF joint motion following conservative or surgical treatment.

Comparison of the linear finite element prediction of deformation and strain of human cancellous bone to 3D digital volume correlation measurements.

The mechanical properties of cancellous bone and the biological response of the tissue to mechanical loading are related to deformation and strain in the trabeculae during function. Due to the small size of trabeculae, their motion is difficult to measure. To avoid the need to measure trabecular motions during loading the finite element method has been used to estimate trabecular level mechanical deformation. This analytical approach has been empirically successful in that the analytical models are solvable and their results correlate with the macroscopically measured stiffness and strength of bones. The present work is a direct comparison of finite element predictions to measurements of the deformation and strain at near trabecular level. Using the method of digital volume correlation, we measured the deformation and calculated the strain at a resolution approaching the trabecular level for cancellous bone specimens loaded in uniaxial compression. Smoothed results from linearly elastic finite element models of the same mechanical tests were correlated to the empirical three-dimensional (3D) deformation in the direction of loading with a coefficient of determination as high as 97% and a slope of the prediction near one. However, real deformations in the directions perpendicular to the loading direction were not as well predicted by the analytical models. Our results show, that the finite element modeling of the internal deformation and strain in cancellous bone can be accurate in one direction but that this does not ensure accuracy for all deformations and strains.

Variability of trabecular microstructure is age-, gender-, race- and anatomic site-dependent and affects stiffness and stress distribution properties of human vertebral cancellous bone☆

Cancellous bone microstructure is an important determinant of the mechanical integrity of vertebrae. The numerous microstructural parameters that have been studied extensively are generally represented as a single value obtained as an average over a sample. The range of the intra-sample variability of cancellous microstructure and its effect on the mechanical properties of bone are less well-understood. The objectives of this study were to investigate the extent to which human cancellous bone microstructure within a vertebra i) is related to bone modulus and stress distribution properties and ii) changes along with age, gender and locations thoracic 12 (T12) vs lumbar 1 (L1). Vertebrae were collected from 15 male (66±15 years) and 25 female (54±16 years) cadavers. Three dimensional finite element models were constructed using microcomputed tomography images of cylindrical specimens. Linear finite element models were used to estimate apparent modulus and stress in the cylinders during uniaxial compression. The intra-specimen mean, standard deviation (SD) and coefficient of variation (CV) of microstructural variables were calculated. Mixed model statistical analysis of the results demonstrated that increases in the intra-specimen variability of the microstructure contribute to increases in the variability of trabecular stresses and decreases in bone stiffness. These effects were independent from the contribution from intra-specimen average of the microstructure. Further, the effects of microstructural variability on bone stiffness and stress variability were not accounted for by connectivity and anisotropy. Microstructural variability properties (SD, CV) generally increased with age, were greater in females than in males and in T12 than in L1. Significant interactions were found between age, gender, vertebra and race. These interactions suggest that microstructural variability properties varied with age differently between genders, races and vertebral levels. The current results collectively demonstrate that microstructural variability has a significant effect on mechanical properties and tissue stress of human vertebral cancellous bone. Considering microstructural variability could improve the understanding of bone fragility and improve assessment of vertebral fracture risk.

Glenoid inclination: In vivo measures in rotator cuff tear patients and associations with superior glenohumeral joint translation

Glenoid inclination has been associated with rotator cuff tears and superior humeral translation, but the relationship between glenoid inclination and superior humeral translation has not been assessed in vivo. This study compared glenoid inclination between repaired and contralateral shoulders in 21 unilateral rotator cuff repair patients. As a secondary analysis, we assessed the relationship between glenoid inclination and in vivo superior humeral translation. Glenoid inclination was measured from patient-specific, computed tomography-based bone models. Glenohumeral joint motion was measured from biplane radiographs collected during coronal-plane abductions. Glenoid inclination was significantly lower for the rotator cuff tear shoulders (90.7 degrees ) than the asymptomatic, contralateral shoulders (92.3 degrees , P = .04). No significant correlation existed between increased glenoid inclination and superior-inferior translation of the uninjured shoulder (P > .30). This study failed to support the theory that glenoid inclination is responsible for superior humeral translation and the development of subacromial impingement.

Accuracy and Feasibility of Dual Fluoroscopy and Model-Based Tracking to Quantify In Vivo Hip Kinematics During Clinical Exams

Accurate measurements of in-vivo hip kinematics may elucidate the mechanisms responsible for impaired function and chondrolabral damage in hips with femoroacetabular impingement (FAI). The objectives of this study were to quantify the accuracy and demonstrate the feasibility of using dual fluoroscopy to measure in-vivo hip kinematics during clinical exams used in the assessment of FAI. Steel beads were implanted into the pelvis and femur of two cadavers. Specimens were imaged under dual fluoroscopy during the impingement exam, FABER test, and rotational profile. Bead locations measured with model-based tracking were compared with those measured using dynamic radiostereometric analysis. Error was quantified by bias and precision, defined as the average and standard deviation of the differences between tracking methods, respectively. A normal male volunteer was also imaged during clinical exams. Bias and precision along a single axis did not exceed 0.17 and 0.21 mm, respectively. Comparing kinematics, positional error was less than 0.48 mm and rotational error was less than 0.58°. For the volunteer, kinematics were reported as joint angles and bone-bone distance. These results demonstrate that dual fluoroscopy and model-based tracking can accurately measure hip kinematics in living subjects during clinical exams of the hip.

In vivo measurement of glenohumeral joint contact patterns

The objectives of this study were to describe a technique for measuring in-vivo glenohumeral joint contact patterns during dynamic activities and to demonstrate application of this technique. The experimental technique calculated joint contact patterns by combining CT-based 3D bone models with joint motion data that were accurately measured from biplane x-ray images. Joint contact patterns were calculated for the repaired and contralateral shoulders of 20 patients who had undergone rotator cuff repair. Significant differences in joint contact patterns were detected due to abduction angle and shoulder condition (i.e., repaired versus contralateral). Abduction angle had a significant effect on the superior/inferior contact center position, with the average joint contact center of the repaired shoulder 12.1% higher on the glenoid than the contralateral shoulder. This technique provides clinically relevant information by calculating in-vivo joint contact patterns during dynamic conditions and overcomes many limitations associated with conventional techniques for quantifying joint mechanics.

Stiffness of the endplate boundary layer and endplate surface topography are associated with brittleness of human whole vertebral bodies

Stress magnitude and variability as estimated from large scale finite element (FE) analyses have been associated with compressive strength of human vertebral cancellous cores but these relationships have not been explored for whole vertebral bodies. In this study, the objectives were to investigate the relationship of FE-calculated stress distribution parameters with experimentally determined strength, stiffness, and displacement based ductility measures in human whole vertebral bodies, investigate the effect of endplate loading conditions on vertebral stiffness, strength, and ductility and test the hypothesis that endplate topography affects vertebral ductility and stress distributions. Eighteen vertebral bodies (T6-L3 levels; 4 female and 5 male cadavers, aged 40-98 years) were scanned using a flat-panel CT system and followed with axial compression testing with Woods metal as filler material to maintain flat boundaries between load plates and specimens. FE models were constructed using reconstructed CT images and filler material was added digitally. Two different FE models with different filler material modulus simulating Woods metal and intervertebral disc (W-layer and D-layer models) were used. Element material modulus to cancellous bone was based on image gray value. Average, standard deviation, and coefficient of variation of von Mises stress in vertebral bone for W-layer and D-layer models and also the ratios of FE parameters from the two models (W/D) were calculated. Inferior and superior endplate surface topographical distribution parameters were calculated. Experimental stiffness, maximum load and work to fracture had the highest correlation with FE-calculated stiffness while experimental ductility measures had highest correlations with FE-calculated average von Mises stress and W-layer to D-layer stiffness ratio. Endplate topography of the vertebra was also associated with its structural ductility and the distribution parameter that best explained this association was kurtosis of inferior endplate topography. Our results indicate that endplate topography variations may provide insight into the mechanisms responsible for vertebral fractures.

Cancellous bone lamellae strongly affect microcrack propagation and apparent mechanical properties: Separation of patients with osteoporotic fracture from normal controls using a 2D nonlinear finite element method (biomechanical stereology)

Biomechanical stereology is proposed as a two-dimensional (2D) finite element (FE) method to estimate the ability of bone tissue to sustain damage and to separate patients with osteoporotic fracture from normal controls. Briefly, 2D nonlinear compact tension FE models were created from quantitative back scattered electron images taken of iliac crest bone specimens collected from the individuals with or without osteoporotic fracture history. The effects of bone mineral microstructure on predicted bone fracture toughness and microcrack propagation were examined. The 2D FE models were used as surrogates for the real bone tissues. The calculated microcrack propagation results and bone mechanical properties were examined as surrogates for measurements from mechanical testing of actual specimens. The results for the 2D FE simulation separated patients with osteoporotic fracture from normal controls even though only the variability in tissue mineral microstructure was used to build the models. The models were deliberately created to ignore all differences in mean mineralization. Hence, the current results support the following hypotheses: (1) that material heterogeneity is important to the separation of patients with osteoporotic fracture from normal controls; and (2) that 2D nonlinear finite element modeling can produce surrogate mechanical parameters that separate patients with fracture from normal controls.

Postfailure modulus strongly affects microcracking and mechanical property change in human iliac cancellous bone: a study using a 2D nonlinear finite element method

A two-dimensional (2D) finite element (FE) method was used to estimate the ability of bone tissue to sustain damage as a function of postfailure modulus. Briefly, 2D nonlinear compact-tension FE models were created from quantitative back-scattered electron images taken of human iliac crest bone specimens. The effects of different postfailure moduli on predicted microcrack propagation were examined. The 2D FE models were used as surrogates for real bone tissues. The crack number was larger in models with higher postfailure modulus, while mean crack length and area were smaller in these models. The rate of stiffness reduction was greater in the models with lower postfailure modulus. Hence, the current results supported the hypothesis that hard tissue postfailure properties have strong effects on bone microdamage morphology and the rate of change in apparent mechanical properties.