髋关节假体脱位分析软件的设计与开发

目的研发了髋关节假体脱位分析软件。软件基于ADAMS/VIEW软件进行二次开发,能够可视化的构建三维参数化的髋关节假体模型,模型能够模拟假体的六种运动。应用此软件可以为病人更好的选择和植入假体,评估各种髋关节假体的安全活动范围,用来设计新的关节假体。     

Hip endoprosthesis for in vivo measurement of joint force and temperature.

Friction between the prosthetic head and acetabular cup increases the temperature in hip implants during activities like walking. A hip endoprosthesis was instrumented with sensors to measure the joint contact forces and the temperature distribution along the entire length of the titanium implant. Sensors and two inductively powered telemetry units are placed inside the hip implant and hermetically sealed against body fluids. Each telemetry unit contains an integrated 8-channel telemetry chip and a radio frequency transmitter. Force, temperature and power supply data are transmitted at different frequencies by two antennas to an external twin receiver. The inductive power supply is controlled by a personal computer. Force and temperature are monitored in real time and all data are stored on a video tape together with the patient's images. This paper describes the design and accuracy of the instrumented implant and the principal function of the external system components.     

Determination of patient-specific multi-joint kinematic models through two-level optimization

Dynamic patient-specific musculoskeletal models have great potential for addressing clinical problems in orthopedics and rehabilitation. However, their predictive capability is limited by how well the underlying kinematic model matches the patient's structure. This study presents a general two-level optimization procedure for tuning any multi-joint kinematic model to a patient's experimental movement data. An outer level optimization modifies the model's parameters (joint position and orientations) while repeated inner level optimizations modify the model's degrees of freedom given the current parameters, with the goal of minimizing errors between model and experimental marker trajectories. The approach is demonstrated by fitting a 27 parameter, three-dimensional, 12 degree-of-freedom lower-extremity kinematic model to synthetic and experimental movement data for isolated joint (hip, knee, and ankle) and gait (full leg) motions. For noiseless synthetic data, the approach successfully recovered the known joint parameters to within an arbitrarily tight tolerance. When noise was added to the synthetic data, root-mean-square (RMS) errors between known and recovered joint parameters were within 10.4 degrees and 10 mm. For experimental data, RMS marker distance errors were reduced by up to 62% compared to methods that estimate joint parameters from anatomical landmarks. Optimized joint parameters found using a loaded full-leg gait motion differed significantly from those found using unloaded individual joint motions. In the future, this approach may facilitate the creation of dynamic patient-specific musculoskeletal models for predictive clinical applications.     

Multi-body simulation of a canine hind limb: model development, experimental validation and calculation of ground reaction forces

Background  

Among other causes the long-term result of hip prostheses in dogs is determined by aseptic loosening. A prevention of prosthesis complications can be achieved by an optimization of the tribological system which finally results in improved implant duration. In this context a computerized model for the calculation of hip joint loadings during different motions would be of benefit. In a first step in the development of such an inverse dynamic multi-body simulation (MBS-) model we here present the setup of a canine hind limb model applicable for the calculation of ground reaction forces.     

A new X-ray-free measurement method for postoperative 3D-position analysis of navigated inserted implants]

In this article a new X-ray-free measurement method for postoperative 3D-position analysis of dental implants is proposed. The advantage of this new method is the possibility of precise 3D comparison to preoperative planning data without the use of X-rays. Standard methods for postoperative implant analysis are performed manually on the basis of orthopantomography or computer tomography scans. The proposed method uses a navigation system and a specially designed measurement device. Analysis of the measurement data shows a mean position error of 0.2 mm +/-0.1 mm and a mean depth error of 0.4 mm +/- 0.3 mm. Thus, this method is suitable for postoperative analysis of the horizontal implant position and for comparison to preoperative planning     

MRI vs CT-based 2D-3D auto-registration accuracy for quantifying shoulder motion using biplane video-radiography

Biplane 2D-3D registration approaches have been used for measuring 3D, in vivo glenohumeral (GH) joint kinematics. Computed tomography (CT) has become the gold standard for reconstructing 3D bone models, as it provides high geometric accuracy and similar tissue contrast to video-radiography. Alternatively, magnetic resonance imaging (MRI) would not expose subjects to radiation and provides the ability to add cartilage and other soft tissues to the models. However, the accuracy of MRI-based 2D-3D registration for quantifying glenohumeral kinematics is unknown. We developed an automatic 2D-3D registration program that works with both CT- and MRI-based image volumes for quantifying joint motions. The purpose of this study was to use the proposed 2D-3D auto-registration algorithm to describe the humerus and scapula tracking accuracy of CT- and MRI-based registration relative to radiostereometric analysis (RSA) during dynamic biplanar video-radiography. The GH kinematic accuracy (RMS error) was 0.6–1.0 mm and 0.6–2.2° for the CT-based registration and 1.4–2.2 mm and 1.2–2.6° for MRI-based registration. Higher kinematic accuracy of CT-based registration was expected as MRI provides lower spatial resolution and bone contrast as compared to CT and suffers from spatial distortions. However, the MRI-based registration is within an acceptable accuracy for many clinical research questions.     

Accuracy of the functional method of hip joint center location: effects of limited motion and varied implementation.

Accurate location of the hip joint center is essential for computation of hip kinematics and kinetics as well as for determination of the moment arms of muscles crossing the hip. The functional method of hip joint center location involves fitting a pelvis-fixed sphere to the path traced by a thigh-fixed point while a subject performs hip motions; the center of this sphere is the hip joint center. The aim of the present study was to evaluate the potential accuracy of the functional method and the dependence of its accuracy on variations in its implementation and the amount of available hip motion. The motions of a mechanical linkage were studied to isolate the factors of interest, removing errors due to skin movement and the palpation of bony landmarks that are always present in human studies. It was found that reducing the range of hip motion from 30 degrees to 15 degrees did significantly increase hip joint center location errors, but that restricting motion to a single plane did not. The magnitudes of these errors, however, even in the least accurate cases, were smaller than those previously reported for either the functional method or other methods based on pelvis measurements of living subjects and cadaver specimens. Neither increasing the number of motion data observations nor analyzing the motion of a single thigh marker (rather than the centroid of multiple markers) was found to significantly increase error. The results of this study (1) imply that the limited range of motion that is often evident in subjects with hip pathology does not preclude accurate determination of the hip joint center when the functional method is used; and (2) provide guidelines for the use of the functional method in human subjects.     

A Novel Approach for Dynamic Testing of Total Hip Dislocation under Physiological Conditions

Constant high rates of dislocation-related complications of total hip replacements (THRs) show that contributing factors like implant position and design, soft tissue condition and dynamics of physiological motions have not yet been fully understood. As in vivo measurements of excessive motions are not possible due to ethical objections, a comprehensive approach is proposed which is capable of testing THR stability under dynamic, reproducible and physiological conditions. The approach is based on a hardware-in-the-loop (HiL) simulation where a robotic physical setup interacts with a computational musculoskeletal model based on inverse dynamics. A major objective of this work was the validation of the HiL test system against in vivo data derived from patients with instrumented THRs. Moreover, the impact of certain test conditions, such as joint lubrication, implant position, load level in terms of body mass and removal of muscle structures, was evaluated within several HiL simulations. The outcomes for a normal sitting down and standing up maneuver revealed good agreement in trend and magnitude compared with in vivo measured hip joint forces. For a deep maneuver with femoral adduction, lubrication was shown to cause less friction torques than under dry conditions. Similarly, it could be demonstrated that less cup anteversion and inclination lead to earlier impingement in flexion motion including pelvic tilt for selected combinations of cup and stem positions. Reducing body mass did not influence impingement-free range of motion and dislocation behavior; however, higher resisting torques were observed under higher loads. Muscle removal emulating a posterior surgical approach indicated alterations in THR loading and the instability process in contrast to a reference case with intact musculature. Based on the presented data, it can be concluded that the HiL test system is able to reproduce comparable joint dynamics as present in THR patients.     

Fusion of radiostereometric analysis data into computed tomography space: application to the elbow joint

Improvement of joint prostheses is dependent upon information concerning the biomechanical properties of the joint. Radiostereometric analysis (RSA) and electromagnetic techniques have been applied in previous cadaver and in vivo studies on the elbow joint to provide valuable information concerning joint motion axes. However, such information is limited to mathematically calculated positions of the axes according to an orthogonal coordinate system and is difficult to relate to individual skeletal anatomy. The aim of this study was to evaluate the in vivo application of a new fusion method to provide three-dimensional (3D) visualization of flexion axes according to bony landmarks. In vivo RSA data of the elbow joint's flexion axes was combined with data obtained by 3D computed tomography (CT). Results were obtained from five healthy subjects after one was excluded due to an instable RSA marker. The median error between imported and transformed RSA marker coordinates and those obtained in the CT volume was 0.22 mm. Median maximal rotation error after transformation of the rigid RSA body to the CT volume was 0.003 degrees . Points of interception with a plane calculated in the RSA orthogonal coordinate system were imported into the CT volume, facilitating the 3D visualization of the flexion axes. This study demonstrates a successful fusion of RSA and CT data, without significant loss of RSA accuracy. The method could be used for relating individual motion axes to a 3D representation of relevant joint anatomy, thus providing important information for clinical applications such as the development of joint prostheses.     

Investigation of the XCAT phantom as a validation tool in cardiac MRI tracking algorithms

PurposeTo describe our magnetic resonance imaging (MRI) simulated implementation of the 4D digital extended cardio torso (XCAT) phantom to validate our previously developed cardiac tracking techniques. Real-time tracking will play an important role in the non-invasive treatment of atrial fibrillation with MRI-guided radiosurgery. In addition, to show how quantifiable measures of tracking accuracy and patient-specific physiology could influence MRI tracking algorithm design.MethodsTwenty virtual patients were subjected to simulated MRI scans that closely model the proposed real-world scenario to allow verification of the tracking technique’s algorithm. The generated phantoms provide ground-truth motions which were compared to the target motions output from our tracking algorithm. The patient-specific tracking error, e_p, was the 3D difference (vector length) between the ground-truth and algorithm trajectories. The tracking errors of two combinations of new tracking algorithm functions that were anticipated to improve tracking accuracy were studied. Additionally, the correlation of key physiological parameters with tracking accuracy was investigated.ResultsOur original cardiac tracking algorithm resulted in a mean tracking error of 3.7 ± 0.6 mm over all virtual patients. The two combinations of tracking functions demonstrated comparable mean tracking errors however indicating that the optimal tracking algorithm may be patient-specific.ConclusionsCurrent and future MRI tracking strategies are likely to benefit from this virtual validation method since no time-resolved 4D ground-truth signal can currently be derived from purely image-based studies.     

A soft tissue artefact model driven by proximal and distal joint kinematics

When analysing human movement through stereophotogrammetry, skin-markers are used. Their movement relative to the underlying bone is known as a soft tissue artefact (STA). A mathematical model to estimate subject- and marker-specific STAs generated during a given motor task, is required for both skeletal kinematic estimators and comparative assessment using simulation. This study devises and assesses such a mathematical model using the paradigmatic case of thigh STAs. The model was based on two hypotheses: (1) that the artefact mostly depends on skin sliding, and thus on the angles of hip and knee; (2) that the relevant relationship is linear. These hypotheses were tested using data obtained from passive hip and knee movements in non-obese specimens and from running volunteers endowed with both skin- and pin-markers.     

Superparamagnetic Iron Oxide Nanoparticles Function as a Long-Term,Multi-Modal Imaging Label for Non-Invasive Tracking of Implanted Progenitor Cells

The purpose of this study was to determine the ability of superparamagnetic iron oxide (SPIO) nanoparticles to function as a long-term tracking label for multi-modal imaging of implanted engineered tissues containing muscle-derived progenitor cells using magnetic resonance imaging (MRI) and X-ray micro-computed tomography (μCT). SPIO-labeled primary myoblasts were embedded in fibrin sealant and imaged to obtain intensity data by MRI or radio-opacity information by μCT. Each imaging modality displayed a detection gradient that matched increasing SPIO concentrations. Labeled cells were then incorporated in fibrin sealant, injected into the atrioventricular groove of rat hearts, and imaged in vivo and ex vivo for up to 1 year. Transplanted cells were identified in intact animals and isolated hearts using both imaging modalities. MRI was better able to detect minuscule amounts of SPIO nanoparticles, while μCT more precisely identified the location of heavily-labeled cells. Histological analyses confirmed that iron oxide particles were confined to viable, skeletal muscle-derived cells in the implant at the expected location based on MRI and μCT. These analyses showed no evidence of phagocytosis of labeled cells by macrophages or release of nanoparticles from transplanted cells. In conclusion, we established that SPIO nanoparticles function as a sensitive and specific long-term label for MRI and μCT, respectively. Our findings will enable investigators interested in regenerative therapies to non-invasively and serially acquire complementary, high-resolution images of transplanted cells for one year using a single label.     

A musculo-skeletal model of rat ankle motion and its experimental test on rat

A biomechanical musculo-skeletal model of functional electrical stimulation (FES)-induced rat ankle motion was implemented and tested in rat experiments. The muscle model is a new Hill-based model which includes established physiological relations of force-velocity and force-length-frequency. However, the series-elastic component and the activation component of previous Hill-based models are replaced by a new component which accounts for dynamic time delays and recruitment that occur in real muscle force generation during limb movements. The skeletal model includes gravity and dynamic forces that occur in real rat ankle motions. In computer simulations, various FES patterns were applied to the tibialis anterior (TA) and soleus (SO) model muscles to produce walk-like ankle motions. In lab experiments, the same stimulation patterns were applied by epimysial electrodes implanted in the TA and SO muscles of live rats cordotomized at level T7. The resulting rat motions were recorded by video camera. Video data was converted to ankle angle-vs-time files for comparison with corresponding model angle-vs-time files. Over a physiologically significant range of ankle motions, model parameters were adjustable to yield model motions that agreed with rat motions to within 2 degrees (root mean square differences of rat and model ankle angles). This is shown in plots of model and rat motions presented here for representative cases of FES. The accuracy of our model in reproducing real ankle motions supports the hypothesis that our new muscle model generates correct muscle forces over a useful range of limb motions. It suggests that the model may be useful in the design of FES neural prostheses.     

Brief communication: Three‐dimensional motion analysis of hindlimb during brachiation in a white‐handed gibbon (Hylobates lar)

In brachiating gibbons, it is thought that there is little movement in the hindlimb joints and that lateral body movement is quite limited. These hypotheses are based on naked‐eye observations, and no quantitative motion analyses of the hindlimbs have been reported. This study quantitatively describes the three‐dimensional movements of the lower trunk and distal thigh during continuous‐contact brachiation in a white‐handed gibbon (Hylobates lar) to evaluate the roles of the trunk and hindlimb. The results revealed that the lower trunk moved both laterally and vertically. The lateral movement of the lower trunk resulted from the lateral inclination of the trunk by gravity. The vertical movement of the trunk was converted into forward velocity, indicating an exchange between potential and kinetic energy. We also observed flexion and extension of the hip, although the excursion was within a small range. In addition, the lateral movement of the hindlimb in thedirection opposite to that of trunk movement helped to reduce the lateral sway of the body. These results suggest that during continuous‐contact brachiation a gibbon uses hip flexion and extension motions to increase the kinetic energy in the swing. In addition, fine motions of the hip may restrict the lateral sway of the center of body mass. Am J Phys Anthropol 142:650–654, 2010. © 2010 Wiley‐Liss, Inc.     

Development of a statistical shape-function model of the implanted knee for real-time prediction of joint mechanics

Outcomes of total knee arthroplasty (TKA) are dependent on surgical technique, patient variability, and implant design. Non-optimal design or alignment choices may result in undesirable contact mechanics and joint kinematics, including poor joint alignment, instability, and reduced range of motion. Implant design and surgical alignment are modifiable factors with potential to improve patient outcomes, and there is a need for robust implant designs that can accommodate patient variability. Our objective was to develop a statistical shape-function model (SFM) of a posterior stabilized implanted knee to instantaneously predict joint mechanics in an efficient manner. Finite element methods were combined with Latin hypercube sampling and regression analyses to produce modeling equations relating nine implant design and six surgical alignment parameters to tibiofemoral (TF) joint mechanics outcomes during a deep knee bend. A SFM was developed and TF contact mechanics, kinematics, and soft tissue loads were instantaneously predicted from the model. Average normalized root-mean-square error predictions were between 2.79% and 9.42%, depending on the number of parameters included in the model. The statistical shape-function model generated instantaneous joint mechanics predictions using a maximum of 130 training simulations, making it ideally suited for integration into a patient-specific design and alignment optimization pipeline. Such a tool may be used to optimize kinematic function to achieve more natural motion or minimize implant wear, and may aid the engineering and clinical communities in improving patient satisfaction and surgical outcomes.     

A novel methodology for generating 3D finite element models of the hip from 2D radiographs

Finite element (FE) modelling has been proposed as a tool for estimating fracture risk and patient-specific FE models are commonly based on computed tomography (CT). Here, we present a novel method to automatically create personalised 3D models from standard 2D hip radiographs. A set of geometrical parameters of the femur were determined from seven a–p hip radiographs and compared to the 3D femoral shape obtained from CT as training material; the error in reconstructing the 3D model from the 2D radiographs was assessed. Using the geometry parameters as the input, the 3D shape of another 21 femora was built and meshed, separating a cortical and trabecular compartment. The material properties were derived from the homogeneity index assessed by texture analysis of the radiographs, with focus on the principal tensile and compressive trabecular systems. The ability of these FE models to predict failure load as determined by experimental biomechanical testing was evaluated and compared to the predictive ability of DXA. The average reconstruction error of the 3D models was 1.77 mm (±1.17 mm), with the error being smallest in the femoral head and neck, and greatest in the trochanter. The correlation of the FE predicted failure load with the experimental failure load was r²=64% for the reconstruction FE model, which was significantly better (p<0.05) than that for DXA (r²=24%). This novel method for automatically constructing a patient-specific 3D finite element model from standard 2D radiographs shows encouraging results in estimating patient-specific failure loads.     

Evaluation of a subject-specific female gymnast model and simulation of an uneven parallel bar swing

A gymnast model and forward dynamics simulation of a dismount preparation swing on the uneven parallel bars were evaluated by comparing experimental and predicted joint positions throughout the maneuver. The bar model was a linearly elastic spring with a frictional bar/hand interface, and the gymnast model consisted of torso/head, arm and two leg segments. The hips were frictionless balls and sockets, and shoulder movement was planar with passive compliant structures approximated by a parallel spring and damper. Subject-specific body segment moments of inertia, and shoulder compliance were estimated. Muscles crossing the shoulder and hip were represented as torque generators, and experiments quantified maximum instantaneous torques as functions of joint angle and angular velocity. Maximum torques were scaled by joint torque activations as functions of time to produce realistic motions. The downhill simplex method optimized activations and simulation initial conditions to minimize the difference between experimental and predicted bar-center, shoulder, hip, and ankle positions. Comparing experimental and simulated performances allowed evaluation of bar, shoulder compliance, joint torque, and gymnast models. Errors in all except the gymnast model are random, zero mean, and uncorrelated, verifying that all essential system features are represented. Although the swing simulation using the gymnast model matched experimental joint positions with a 2.15cm root-mean-squared error, errors are correlated. Correlated errors indicate that the gymnast model is not complex enough to exactly reproduce the experimental motion. Possible model improvements including a nonlinear shoulder model with active translational control and a two-segment torso would not have been identified if the objective function did not evaluate the entire system configuration throughout the motion. The model and parameters presented in this study can be effectively used to understand and improve an uneven parallel bar swing, although in the future there may be circumstances where a more complex model is needed.     

Changes in Muscle Volume and Composition After Treatment of Hip Dysplasia with Periacetabular Osteotomy

BackgroundPeriacetabular osteotomy (PAO) is a common treatment for pre-arthritic hip dysplasia in young adults. The purpose of this study was to better understand changes in muscle volume and composition after PAO visualized using magnetic resonance imaging (MRI).MethodsA prospectively collected series of individuals that underwent PAO for hip dysplasia were reviewed to identify subjects with pre- and postoperative MRI. In our practice, MRI was obtained preoperatively and greater than 6 months after PAO for persistent hip pain. MRI sequences were selected to optimize visualization of the muscle volume, fatty infiltration, and hip joint cartilage. MRI images were selected at predetermined bony landmarks and analyzed using 3D Slicer (©2021, www.slicer.org) software to measure muscle diameter and calculate muscle cross-sectional area (CSA) in 17 individual muscles surrounding the hip. Muscle atrophy was graded using the Goutallier classification for fatty infiltration and acetabular cartilage condition was graded using the Outerbridge classification. We compared pre- and postoperative muscle area and composition as well as cartilage for each case.ResultsA series of six female patients met our inclusion criteria. Mean age was 26 years at time of surgery. All cases had MRI sequences adequate for muscle volume measurements. Fatty infiltration and cartilage changes were recorded in four subjects with appropriate MRI sequences. Separating muscle groups, external rotators underwent the largest volume increase. Hip flexors demonstrated mild volume decrease. CSA change among external rotators averaged +12%, hip flexors -9.3%, and hip abductors -9.2% after PAO. All muscles had either the same or increased fatty infiltration after surgery, with gluteus medius and iliacus undergoing the most average increase. Similarly, cartilage condition worsened by a small margin in this series.ConclusionOur results provide preliminary indication that PAO may have noticeable effects on muscle characteristics and cartilage in the early postoperative period. This was a limited case series of subjects with adequate pre- and post-operative MRI imaging.Level of Evidence: IV     

Monitoring hip posture in total hip arthroplasty using an inertial measurement unit-based hip smart trial system: An in vitro validation experiment using a fixed pelvis model

Intraoperative measurement of hip posture is the basis for assessing hip range of motion (ROM) and predicting postoperative functional limits allowable for activities of daily living. Although computer navigation for total hip arthroplasty (THA) has improved the accuracy of intraoperative ROM evaluation, it has not gained widespread popularity due to its complex and time-consuming protocol. We therefore developed an inertial measurement unit-based hip smart trial system (IMUHST) for intraoperative monitoring of hip posture. An in vitro validation experiment was conducted using bone models with a three-dimensional measurement model as the reference standard. The absolute mean error, Bland – Altman analysis and intra-class correlation coefficient demonstrated that the validity and reliability of this system meets the requirement for clinical application. Given that monitoring posture is the basis for evaluating the direction(s) of potential impingement, subluxation and dislocation, the IMUHST is a promising development direction of computer assisted surgery in THA.     

Validation of a low-dose hybrid RSA and fluoroscopy technique: Determination of accuracy, bias and precision

Analyzing skeletal kinematics with radiostereometric analysis (RSA) following corrective orthopedic surgery allows the quantitative comparison of different implant designs. The purpose of this study was to validate a technique for dynamically estimating the relative position and orientation of skeletal segments using RSA and single plane X-ray fluoroscopy. Two micrometer-based in vitro phantom models of the skeletal segments in the hip and knee joints were used. The spatial positions of tantalum markers that were implanted into each skeletal segment were reconstructed using RSA. The position and orientation of each segment were determined in fluoroscopy images by minimizing the difference between the markers measured and projected in the image plane. Accuracy was determined in terms of bias and precision by analyzing the deviation between the applied displacement protocol and measured pose estimates. Measured translational accuracy was less than 100 microm parallel to the image plane and less than 700 microm in the direction orthogonal to the image plane. The measured rotational error was less than 1 degrees . Measured translational and rotational bias was not statistically significant at the 95% level of confidence. The technique allows real-time kinematic skeletal measurements to be performed on human subjects implanted with tantalum markers for quantitatively measuring the motion of normal joints and different implant designs.