Mechanical behavior of bovine nasal cartilage under static and dynamic loading

Abnormal mechanical loading may trigger cartilage degeneration associated with osteoarthritis. Tissue response to load has been the subject of several in vitro studies. However, simple stimuli were often applied, not fully mimicking the complex in vivo conditions. Therefore, a rolling/plowing explant test system (RPETS) was developed to replicate the combined in vivo loading patterns. In this work we investigated the mechanical behavior of bovine nasal septum (BNS) cartilage, selected as tissue approximation for experiments with RPETS, under static and dynamic loading. Biphasic material properties were determined and compared with those of other cartilaginous tissues. Furthermore, dynamic loading in plowing modality was performed to determine dynamic response and experimental results were compared with analytical models and Finite Elements (FE) computations. Results showed that BNS cartilage can be modeled as a biphasic material with Young's modulus E=2.03±0.7 MPa, aggregate modulus H_A=2.35±0.7 MPa, Poisson's ratio ν=0.24±0.07, and constant hydraulic permeability k₀=3.0±1.3×10⁻¹⁵ m⁴ (N s)⁻¹. Furthermore, dynamic analysis showed that plowing induces macroscopic reactions in the tissue, proportionally to the applied loading force. The comparison among analytical, FE analysis and experimental results showed that predicted tangential forces and sample deformation lay in the range of variation of experimental results for one specific experimental condition. In conclusion, mechanical properties of BNS cartilage under both static and dynamic compression were assessed, showing that this tissue behave as a biphasic material and has a viscoelastic response to dynamic forces.     

Evaluation of an electromagnetic position tracking device for measuring in vivo, dynamic joint kinematics

An electromagnetic position tracking device was evaluated to determine its static and dynamic accuracy and reliability for applications related to measuring in vivo joint kinematics. The device detected the position and orientation of small coiled sensors, maintained in an electromagnetic field. System output was measured against known translations or rotations throughout the measurement volume. Average translational errors during static testing were 0.1 +/- 0.04, 0.2 +/- 0.17, and 0.8 +/- 0.81 mm (mean+/-SD) for sensors 50, 300, and 550 mm away from the field generator, respectively. Average rotational errors were 0.4 +/- 0.31 degrees, 0.4 +/- 0.21 degrees, and 0.9 +/- 0.85 degrees (mean +/- SD) for sensors located at the same distances. Since we intended to use this system in an animal walking on a treadmill, we incrementally moved the sensors under various treadmill conditions. The effects of treadmill operation on translational accuracy were found to be negligible. The effects of dynamic motions on sensor-to-sensor distance were also assessed for future data collection in the animal. Sensor-to-sensor distance showed standard deviations of 2.6 mm and a range of 13 mm for the highest frequency tested (0.23 Hz). We conclude that this system is useful for static or slow dynamic motions, but is of limited use for obtaining gait kinematics at higher speeds.     

Comparison of different camera calibration approaches for underwater applications

The purpose of this study was to compare three camera calibration approaches applied to underwater applications: (1) static control points with nonlinear DLT; (2) moving wand with nonlinear camera model and bundle adjustment; (3) moving plate with nonlinear camera model. The DVideo kinematic analysis system was used for underwater data acquisition. The system consisted of two gen-locked Basler cameras working at 100 Hz, with wide angle lenses that were enclosed in housings. The accuracy of the methods was compared in a dynamic rigid bar test (acquisition volume-4.5×1×1.5 m(3)). The mean absolute errors were 6.19 mm for the nonlinear DLT, 1.16 mm for the wand calibration, 1.20 mm for the 2D plate calibration using 8 control points and 0.73 mm for the 2D plane calibration using 16 control points. The results of the wand and 2D plate camera calibration methods were less associated to the rigid body position in the working volume and provided better accuracy than the nonlinear DLT. Wand and 2D plate camera calibration methods presented similar and highly accurate results, being alternatives for underwater 3D motion analysis.     

Accuracy of determining the point of force application with piezoelectric force plates

The accuracy of determining the point of force application with piezoelectric force plates, as specified by the manufacturer, is lower than needed for certain applications. The purpose of this study was to evaluate the accuracy of a commonly used plate (KISTLER type 9287) and to improve it by proposing a correction algorithm. Forces were applied to a wooden board, supported in one corner by a stylus that rested on the force plate. To determine the influence of position and magnitude of the force vector, the stylus was placed on 117 different locations, and calibrated masses were used to exert vertical forces between 0 and 2000 N. To determine the influence of loading rate, dynamic tests were performed in which a subject ran across the board. In static tests at a given stylus position with actual coordinates x (short axis) and y (long axis), it was found that the calculated coordinates x and y of the point of force application had virtually constant values at forces above 1000 N. In dynamic tests, oscillations could occur in x and y with an amplitude of more than 20 mm. When these were avoided or removed by filtering, static and dynamic tests at a given stylus position showed the same values for x and y at forces above 1000 N. Across stylus positions, the errors x-x and y-y (measured at 1600 N) ranged from -20 to +20 mm. The average over 117 points of the absolute errors magnitude of x-x and magnitude of y-y amounted to 3.5 and 6.3 mm, respectively (mean values of three plates).(ABSTRACT TRUNCATED AT 250 WORDS)     

The static accuracy and calibration of inertial measurement units for 3D orientation

Inertial measurement units (IMUs) are integrated electronic devices that contain accelerometers, magnetometers and gyroscopes. Wearable motion capture systems based on IMUs have been advertised as alternatives to optical motion capture. In this paper, the accuracy of five different IMUs of the same type in measuring 3D orientation in static situations, as well as the calibration of the accelerometers and magnetometers within the IMUs, has been investigated. The maximum absolute static orientation error was 5.2°, higher than the 1° claimed by the vendor. If the IMUs are re-calibrated at the time of measurement with the re-calibration procedure described in this paper, it is possible to obtain an error of less than 1°, in agreement with the vendor's specifications (XSens Technologies B.V. 2005. Motion tracker technical documentation Mtx-B. Version 1.03. Available from: www.xsens.com).

The new calibration appears to be valid for at least 22 days providing the sensor is not exposed to high impacts. However, if several sensors are ‘daisy chained’ together changes to the magnetometer bias can cause heading errors of up to 15°. The results demonstrate the non-linear relationship between the vendor's orthogonality claim of < 0.1° and the accuracy of 3D orientation obtained from factory calibrated IMUs in static situations. The authors hypothesise that the high magnetic dip (64°) in our laboratory may have exacerbated the errors reported. For biomechanical research, small relative movements of a body segment from a calibrated position are likely to be more accurate than large scale global motion that may have an error of up to 9.8°.     

A motion analysis marker-based method of determining centre of pressure during two-legged hopping

The fixed position of force plates has led researchers to pursue alternative methods of determining centre of pressure (CoP) location. To date, errors reported using alternative methods to the force plate during dynamic tasks have been high. The aim of this study was to investigate the accuracy of a motion analysis marker-based system to determine CoP during a two-legged hopping task. Five markers were attached to the left and right feet of eight healthy adults (5 females, 3 males, age: 25.0±2.8 years, height: 1.75±0.07 m, mass: 71.3±11.3 kg). Multivariate forward stepwise and forced entry linear regression was used with data from five participants to determine CoP position during quiet standing and hopping at various frequencies. Maximum standard error of the estimate of CoP position was 12 mm in the anteroposterior direction and 8 mm in the mediolateral. Cross-validation was performed using the remaining 3 participants. Maximum root mean square difference between the force plate and marker method was 14 mm for mediolateral CoP and 20 mm for anteroposterior CoP during 1.5 Hz hopping. Differences reduced to a maximum of 7 mm (mediolateral) and 14 mm (anteroposterior) for the other frequencies. The smallest difference in calculated sagittal plane ankle moment and timing of maximum moment was during 3.0 Hz hopping, and largest at 1.5 Hz. Results indicate the marker-based method of determining CoP may be a suitable alternative to a force plate to determine CoP position during a two-legged hopping task at frequencies greater than 1.5 Hz.     

The static accuracy and calibration of inertial measurement units for 3D orientation

Inertial measurement units (IMUs) are integrated electronic devices that contain accelerometers, magnetometers and gyroscopes. Wearable motion capture systems based on IMUs have been advertised as alternatives to optical motion capture. In this paper, the accuracy of five different IMUs of the same type in measuring 3D orientation in static situations, as well as the calibration of the accelerometers and magnetometers within the IMUs, has been investigated. The maximum absolute static orientation error was 5.2 degrees , higher than the 1 degrees claimed by the vendor. If the IMUs are re-calibrated at the time of measurement with the re-calibration procedure described in this paper, it is possible to obtain an error of less than 1 degrees , in agreement with the vendor's specifications (XSens Technologies B.V. 2005. Motion tracker technical documentation Mtx-B. Version 1.03. Available from: www.xsens.com). The new calibration appears to be valid for at least 22 days providing the sensor is not exposed to high impacts. However, if several sensors are 'daisy chained' together changes to the magnetometer bias can cause heading errors of up to 15 degrees . The results demonstrate the non-linear relationship between the vendor's orthogonality claim of < 0.1 degrees and the accuracy of 3D orientation obtained from factory calibrated IMUs in static situations. The authors hypothesise that the high magnetic dip (64 degrees ) in our laboratory may have exacerbated the errors reported. For biomechanical research, small relative movements of a body segment from a calibrated position are likely to be more accurate than large scale global motion that may have an error of up to 9.8 degrees .     

Design and demonstration of a new instrumented spatial linkage for use in a dynamic environment: application to measurement of ankle rotations during snowboarding

Joint injuries during sporting activities might be reduced by understanding the extent of the dynamic motion of joints prone to injury during maneuvers performed in the field. Because instrumented spatial linkages (ISLs) have been widely used to measure joint motion, it would be useful to extend the functionality of an ISL to measure joint motion in a dynamic environment. The objectives of the work reported by this paper were to (i) design and construct an ISL that will measure dynamic joint motion in a field environment, (ii) calibrate the ISL and quantify its static measurement error, (iii) quantify dynamic measurement error due to external acceleration, and (iv) measure ankle joint complex rotation during snowboarding maneuvers performed on a snow slope. An "elbow-type" ISL was designed to measure ankle joint complex rotation throughout its range (+/-30 deg for flexion/extension, +/-15 deg for internal/external rotation, and +/-15 deg for inversion/eversion). The ISL was calibrated with a custom six degree-of-freedom calibration device generally useful for calibrating ISLs, and static measurement errors of the ISL also were evaluated. Root-mean-squared errors (RMSEs) were 0.59 deg for orientation (1.7% full scale) and 1.00 mm for position (1.7% full scale). A custom dynamic fixture allowed external accelerations (5 g, 0-50 Hz) to be applied to the ISL in each of three linear directions. Maximum measurement deviations due to external acceleration were 0.05 deg in orientation and 0.10 mm in position, which were negligible in comparison to the static errors. The full functionality of the ISL for measuring joint motion in a field environment was demonstrated by measuring rotations of the ankle joint complex during snowboarding maneuvers performed on a snow slope.     

Developing simulations to reproduce in vivo fluoroscopy kinematics in total knee replacement patients

For clinically predictive testing and design-phase evaluation of prospective total knee replacement (TKR) implants, devices should ideally be evaluated under physiological loading conditions which incorporate population-level variability. A challenge exists for experimental and computational researchers in determining appropriate loading conditions for wear and kinematic knee simulators which reflect in vivo joint loading conditions. There is a great deal of kinematic data available from fluoroscopy studies. The purpose of this work was to develop computational methods to derive anterior–posterior (A–P) and internal–external (I–E) tibiofemoral (TF) joint loading conditions from in vivo kinematic data. Two computational models were developed, a simple TF model, and a more complex lower limb model. These models were driven through external loads applied to the tibia and femur in the TF model, and applied to the hip, ankle and muscles in the lower limb model. A custom feedback controller was integrated with the finite element environment and used to determine the external loads required to reproduce target kinematics at the TF joint. The computational platform was evaluated using in vivo kinematic data from four fluoroscopy patients, and reproduced in vivo A–P and I–E motions and compressive force with a root-mean-square (RMS) accuracy of less than 1 mm, 0.1°, and 40 N in the TF model and in vivo A–P and I–E motions, TF flexion, and compressive loads with a RMS accuracy of less than 1 mm, 0.1°, 1.4°, and 48 N in the lower limb model. The external loading conditions derived from these models can ultimately be used to establish population variability in loading conditions, for eventual use in computational as well as experimental activity simulations.     

Geometric mouse variation: Implications to the axial ulnar loading protocol and animal specific calibration

Large variations in axial ulnar load strain calibration results suggest that animal-specific calibrations may be necessary. However, the optimal set of geometric measures for performing an animal-specific calibration are not known, potentially as a result of confounding effects associated with experimentally introduced variation. The purpose of this study was to characterize the inherent variability of ulnar geometric measures known to influence periosteal midshaft strain during axial ulnar exogenous loading, and to further quantify the relationship between the variance of those geometric measures and periosteal strain during axial loading. Thirty-nine right mouse forelimbs were scanned with microCT. Seven geometric measures that influence periosteal strain resulting from a combined axial and bending loading were computed and used to estimate animal-specific strains on the periosteal midshaft. Animal specific strains were estimated using a theoretical model based on the generalized flexure formula. The predicted mean and standard deviation of the simulated midshaft strain gauge measurement resulting from the inter-animal geometric differences was −985±148 με/N. The complete beam bending term associated with bending about the I_min axis accounted for 89% of the variance and reduced the residual RMSE to 50.4 με. Eccentricity associated with the axial loading contributed a substantial portion of variation to the computed strain suggesting that calibration procedures to account for animal differences should incorporate that variable. The method developed in this study provides a relatively simple procedure for computing animal-specific strains using microCT scan data, without the need of a load/strain calibration study or computationally intensive finite element models.     

Evaluation and calibration of an electromagnetic tracking device for biomechanical analysis of lifting tasks.

Electromagnetic motion tracking devices are increasingly used as a kinematic measuring tool. The aim of this study was to evaluate a long-range transmitter in an environment with a conventional force plate present in order to assess its suitability for further biomechanical applications. Using a calibration apparatus developed in our lab and Optotrack measurements, the performances of the Motion Star were evaluated. Positions and orientations were measured in a 140 x 80 x 120 cm(3) space centered on the force plate. Using a mathematical model developed at Queen's University, these data were calibrated. Errors on position and orientation were less than 150 mm and 10 degrees before calibration of the Motion Star, and less than 20mm and 2 degrees after calibration, with no differences between data collected with the force plate switched on/off. These errors did not depend on sensor orientation. Variability of the signal was small indicating minimal noise. Field distortion was the largest source of measurement error, which increased with the distance between the transmitter and the sensor and the proximity of the sensor to the force plate. Before its use for biomechanical analysis of lifting tasks and validation of dynamic models using force plate data, the data from electromagnetic motion tracking devices must be calibrated to decrease the errors due to electromagnetic field distortion.     

Validation of imaging-based quantification of glenohumeral joint kinematics using an unmodified clinical biplane fluoroscopy system

Model-based tracking, using CT and biplane fluoroscopy, allows highly accurate quantification of glenohumeral motion and changes in the subacromial space. Previous investigators have used custom-built biplane fluoroscopes designed specifically for kinematic applications, which are available at few institutions and require FDA approval prior to clinical use. The aim of this study was to demonstrate the utility of an off-the-shelf clinical biplane fluoroscope for kinematic applications by validating model-based tracking for measurement of glenohumeral motion using an unmodified clinical system. Biplane images of each shoulder of a cadaver torso were acquired at various joint positions and during simulated movements along anatomical planes of motion. The pose of each humerus and scapula was determined using model-based tracking and compared to a bead-based gold standard. Error due to a temporal-offset between corresponding biplane images, characteristic of clinical biplane systems, was determined by comparison of measured and known relative position of 2 bead clusters of a phantom that was imaged while moved throughout the fluoroscopy image volume. Model-based tracking had global kinematic mean absolute errors of 0.27 mm and 0.29° (static), and 0.22–0.32 mm and 0.12–0.45° (dynamic). Glenohumeral mean absolute errors were 0.39 mm and 0.45° (static), and 0.36–0.42 mm and 0.41–0.48° (dynamic). The temporal-offset was predicted to add errors of 0.06–0.85 mm and 0.05–0.28° for cadaveric trials for the speeds examined. For defined speeds, sub-millimeter and sub-degree kinematic accuracy and precision were achieved using an unmodified clinical biplane fluoroscope for quantification of glenohumeral motion.     

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 specimen's 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.     

Fracture mechanics of the femoral neck in a composite bone model: Effects of platen geometry

Load applicator (platen) geometry used for axial load to failure testing of the femoral neck varies between studies and the biomechanical consequences are unknown. The purpose of this study was to determine if load application with a flat versus a conical platen results in differing fracture mechanics. Femurs were aligned in 25° of adduction and an axial compressive force was applied to the femoral heads at a rate of 6 mm/min until failure. Load application with the conical platen resulted in an average ultimate failure load, stiffness, and energy to failure of 9067 N, 4033 N/mm, and 12.12 J, respectively. Load application with the flat platen resulted in a significant (p<0.05) reduction in ultimate failure load (7620 N) and stiffness (2924 N/mm). Energy to failure (12.30 J) was not significantly different (p=0.893). Different fracture patterns were observed for the two platens and the conical platen produced fractures more similar to clinical observations. Use of a flat platen underestimates the strength and stiffness of the femoral neck and inaccurately predicts the associated fracture pattern. These findings must be considered when interpreting the results of prior biomechanical studies on femoral neck fracture and for the development of future femoral neck fracture models.     

Sensitivity deviation: instrumental linearity errors that influence concentration analyses and kinetic evaluation of biomolecular interactions.

Many scientific instruments utilise multiple element detectors, e.g. CCD's or photodiode arrays, to monitor the change in a position of an optical pattern. For example. instruments for affinity biosensing based on surface plasmon resonance (SPR) or resonant mirror are equipped with such detectors. An important and desired property of these bioanalytical instruments is that the calculation of the movement or change in shape follows the true change. This is often not the case and it may lead to linearity errors, and to sensitivity errors. The sensitivity is normally defined as the slope of the calibration curve. A new parameter is introduced to account for the linearity errors, the sensitivity deviation, defined as the deviation from the undistorted slope of the calibration curve. The linearity error and the sensitivity deviation are intimately related and the sensitivity deviation may lead to misinterpretation of kinetic data, mass transport limitations and concentration analyses. Because the linearity errors are small (e.g. 10 pg/mm2 of biomolecules on the sensor surface) with regard to the dynamic range (e.g. 30,000 pg/mm2), they can be difficult to discover. However, the linearity errors are often not negligible with regard to a typical response (e.g. 0-100 pg/mm2). and may therefore cause serious problems. A method for detecting linearity errors is outlined. Further on, this paper demonstrates how integral linearity errors of less than 1% can result in a sensitivity deviation of 10%, a value that in our opinion cannot be ignored in biospecific interaction analysis (BIA). It should also be stressed out that this phenomenon also occurs in other instruments using array detectors.     

A neural network model rapidly learning gains and gating of reflexes necessary to adapt to an arm's dynamics

Effects of dynamic coupling, gravity, inertia and the mechanical impedances of the segments of a multi-jointed arm are shown to be neutralizable through a reflex-like operating three layer static feedforward network. The network requires the proprioceptively mediated actual state variables (here angular velocity and position) of each arm segment. Added neural integrators (and/or differentiators) can make the network exhibit dynamic properties. Then, actual feedback is not necessary and the network can operate in a pure feedforward fashion. Feedforward of an additional load can easily be implemented into the network using descendent gating, and a negative feedback control loop added to the feedforward control reduces errors due to external noise. A training, which combines a least squared error based simultaneous learning rule (LSQ-rule) with a self-imitation algorithm based on direct inverse modeling, enables the network to acquire the whole inverse dynamics, limb parameters included, during one short training movement. The considerations presented also hold for multi-jointed manipulators.     

Dynamic compression can inhibit chondrogenesis of mesenchymal stem cells

The objective of this study was to investigate the influence of dynamic compressive loading on chondrogenesis of mesenchymal stem cells (MSCs) in the presence of TGF-β3. Isolated porcine MSCs were suspended in 2% agarose and subjected to intermittent dynamic compression (10% strain) for a period of 42 days in a dynamic compression bioreactor. After 42 days in culture, the free-swelling specimens exhibited more intense alcian blue staining for proteoglycans, while immunohistochemical analysis revealed increased collagen type II immunoreactivity. Glycosaminoglycan (GAG) content increased with time for both free-swelling and dynamically loaded constructs, and by day 42 it was significantly higher in both the core (2.5 ± 0.21%w/w vs. 0.94 ± 0.03%w/w) and annulus (1.09 ± 0.09%w/w vs. 0.59 ± 0.08%w/w) of free-swelling constructs compared to dynamically loaded constructs. This result suggests that further optimization is required in controlling the biomechanical and/or the biochemical environment if such stimuli are to have beneficial effects in generating functional cartilaginous tissue.     

Design and demonstration of a dynamometric horseshoe for measuring ground reaction loads of horses during racing conditions

Because musculoskeletal injuries to racehorses are common, instrumentation for the study of factors (e.g. track surface), which affect the ground reaction loads in horses during racing conditions, would be useful. The objectives of the work reported by this paper were to (1) design and construct a novel dynamometric horseshoe that is capable of measuring the complete ground reaction loading during racing conditions, (2) characterize static and dynamic measurement errors, and (3) demonstrate the usefulness of the instrument by collecting example data during the walk, trot, canter, and gallop for a single subject. Using electrical resistance strain gages, a dynamometric horseshoe was designed and constructed to measure the complete ground reaction force and moment vectors and the center of pressure. To mimic the load transfer surface of the hoof, the shape of the surface contacting the ground was similar to that of the solar surface of the hoof. Following static calibration, the measurement accuracy was determined. The root mean squared errors (RMSE) were 3% of full scale for the force component normal to the hoof and 9% for force components in the plane of the hoof. The dynamic calibration determined that the natural frequency with the full weight of a typical horse was 1744 Hz. Example data were collected during walking on a ground surface and during trotting, cantering, and galloping on a treadmill. The instrument successfully measured the complete ground reaction load during all four gaits. Consequently the dynamometric horseshoe is useful for studying factors, which affect ground reaction loads during racing conditions.     

Prediction of the structural response of the femoral shaft under dynamic loading using subject-specific finite element models

The goal of this study was to predict the structural response of the femoral shaft under dynamic loading conditions using subject-specific finite element (SS-FE) models and to evaluate the prediction accuracy of the models in relation to the model complexity. In total, SS-FE models of 31 femur specimens were developed. Using those models, dynamic three-point bending and combined loading tests (bending with four different levels of axial compression) of bare femurs were simulated, and the prediction capabilities of five different levels of model complexity were evaluated based on the impact force time histories: baseline, mass-based scaled, structure-based scaled, geometric SS-FE, and heterogenized SS-FE models. Among the five levels of model complexity, the geometric SS-FE and the heterogenized SS-FE models showed statistically significant improvement on response prediction capability compared to the other model formulations whereas the difference between two SS-FE models was negligible. This result indicated the geometric SS-FE models, containing detailed geometric information from CT images with homogeneous linear isotropic elastic material properties, would be an optimal model complexity for prediction of structural response of the femoral shafts under the dynamic loading conditions. The average and the standard deviation of the RMS errors of the geometric SS-FE models for all the 31 cases was 0.46 kN and 0.66 kN, respectively. This study highlights the contribution of geometric variability on the structural response variation of the femoral shafts subjected to dynamic loading condition and the potential of geometric SS-FE models to capture the structural response variation of the femoral shafts.     

Accuracy of femur reconstruction from sparse geometric data using a statistical shape model

Sparse geometric information from limited field-of-view medical images is often used to reconstruct the femur in biomechanical models of the hip and knee. However, the full femur geometry is needed to establish boundary conditions such as muscle attachment sites and joint axes which define the orientation of joint loads. Statistical shape models have been used to estimate the geometry of the full femur from varying amounts of sparse geometric information. However, the effect that different amounts of sparse data have on reconstruction accuracy has not been systematically assessed. In this study, we compared shape model and linear scaling reconstruction of the full femur surface from varying proportions of proximal and distal partial femur geometry in combination with morphometric and landmark data. We quantified reconstruction error in terms of surface-to-surface error as well as deviations in the reconstructed femur’s anatomical coordinate system which is important for biomechanical models. Using a partial proximal femur surface, mean shape model-based reconstruction surface error was 1.8 mm with 0.15° or less anatomic axis error, compared to 19.1 mm and 2.7–5.6° for linear scaling. Similar results were found when using a partial distal surface. However, varying amounts of proximal or distal partial surface data had a negligible effect on reconstruction accuracy. Our results show that given an appropriate set of sparse geometric data, a shape model can reconstruct full femur geometry with far greater accuracy than simple scaling.