Amy B. Zavatsky

Experimental Validation of a Finite Element Model of a Human Cadaveric Tibia

Finite element (FE) models of long bones are widely used to analyze implant designs. Experimental validation has been used to examine the accuracy of FE models of cadaveric femurs; however, although convergence tests have been carried out, no FE models of an intact and implanted human cadaveric tibia have been validated using a range of experimental loading conditions. The aim of the current study was to create FE models of a human cadaveric tibia, both intact and implanted with a unicompartmental knee replacement, and to validate the models against results obtained from a comprehensive set of experiments. Seventeen strain rosettes were attached to a human cadaveric tibia. Surface strains and displacements were measured under 17 loading conditions, which consisted of axial, torsional, and bending loads. The tibia was tested both before and after implantation of the knee replacement. FE models were created based on computed tomography (CT) scans of the cadaveric tibia. The models consisted of ten-node tetrahedral elements and used 600 material properties derived from the CT scans. The experiments were simulated on the models and the results compared to experimental results. Experimental strain measurements were highly repeatable and the measured stiffnesses compared well to published results. For the intact tibia under axial loading, the regression line through a plot of strains predicted by the FE model versus experimentally measured strains had a slope of 1.15, an intercept of 5.5 microstrain, and an R(2) value of 0.98. For the implanted tibia, the comparable regression line had a slope of 1.25, an intercept of 12.3 microstrain, and an R(2) value of 0.97. The root mean square errors were 6.0% and 8.8% for the intact and implanted models under axial loads, respectively. The model produced by the current study provides a tool for simulating mechanical test conditions on a human tibia. This has considerable value in reducing the costs of physical testing by pre-selecting the most appropriate test conditions or most favorable prosthetic designs for final mechanical testing. It can also be used to gain insight into the results of physical testing, by allowing the prediction of those variables difficult or impossible to measure directly.

Determination of gait patterns in children with spastic diplegic cerebral palsy using principal components

This study developed an objective graphical classification method of spastic diplegic cerebral palsy (CP) gait patterns based on principal component analysis (PCA). Gait analyses of 20 healthy and 20 spastic diplegic CP children were examined to define gait characteristics. PCA was used to reduce the dimensionality of 27 parameters (26 selected kinematics variables and age of the children) for the 40 subjects in order to identify the dominant variability in the data. Fuzzy C-mean cluster analysis was performed plotting the first three principal components, which accounted for 61% of the total variability. Results indicated that only the healthy children formed a distinct cluster; however it was possible to recognise gait patterns in overlapping clusters in children with spastic diplegia. This study demonstrates that it is possible to quantitatively classify gait types in CP using PCA. Graphical classification of gait types could assist in clinical evaluation of the children and serve as a validation of clinical reports as well as aid treatment planning.

Ligament Forces at the Knee during Isometric Quadriceps Contractions

A mathematical model of the knee in the sagittal plane was used to investigate the ligament forces resulting when a posteriorly directed external force, applied to the tibia, resists extension of the knee under increasing isometric quadriceps contractions. The model is based on simple geometric representations of the bones, ligaments and muscles at the knee. An elementary mechanical analysis was used to predict which ligament, the anterior or posterior cruciate, was loaded at a given flexion angle and known line of action of the external force. Ligament force, as a proportion of the external force, was calculated first assuming the ligaments to be represented by single, inextensible lines. Modelling the ligaments as continuous arrays of extensible fibres then showed that tibio-femoral translations and ligament forces increased non-linearly with increasing muscle forces and approached asymptotic values which depended on flexion angle. In most positions of the joint, the calculated asymptotic ligament force values were less than the reported ultimate strength of human ligament, despite quadriceps forces of over three times body weight. The possibility of these asymptotic values of ligament force may explain why, at certain flexion angles, large forces can be developed by the muscles at the knee without ligament rupture.

Cruciate ligament loading during isometric muscle contractions. A theoretical basis for rehabilitation.

A model of the knee in the sagittal plane was used to investigate the ligament forces resulting when an ex ternal load applied to the tibia resisted either extension or flexion of the knee under increasing isometric quad riceps or hamstrings contractions, respectively. An el ementary mechanical analysis showed which ligament, the anterior or posterior cruciate, was loaded at a given flexion angle and known line of action of the external load. Ligament force, as a proportion of the external load, was also calculated. The results serve as guide lines for the design of injury-specific physical therapy techniques for use after cruciate ligament reconstruc tion.

Gait compensations caused by foot deformity in cerebral palsy

Cerebral palsy (CP) is a complex syndrome, with multiple interactions between joints and muscles. Abnormalities in movement patterns can be measured using motion capture techniques, however determining which abnormalities are primary, and which are secondary, is a difficult task. Deformity of the foot has anecdotally been reported to produce compensatory abnormalities in more proximal lower limb joints, as well as in the contralateral limb. However, the exact nature of these compensations is unclear. The aim of this paper was to provide clear and objective criteria for identifying compensatory mechanisms in children with spastic hemiplegic CP, in order to improve the prediction of the outcome of foot surgery, and to enhance treatment planning. Twelve children with CP were assessed using conventional gait analysis along with the Oxford Foot Model prior to and following surgery to correct foot deformity. Only those variables not directly influenced by foot surgery were assessed. Any that spontaneously corrected following foot surgery were identified as compensations. Pelvic rotation, internal rotation of the affected hip and external rotation of the non-affected hip tended to spontaneously correct. Increased hip flexion on the affected side, along with reduced hip extension on the non-affected side also appeared to be compensations. It is likely that forefoot supination occurs secondary to deviations of the hindfoot in the coronal plane. Abnormal activity in the tibialis anterior muscle may be consequent to tightness and overactivity of the plantarflexors. On the non-affected side, increased plantarflexion during stance also resolved following surgery to the affected side.

Simultaneous in vitro measurement of patellofemoral kinematics and forces following Oxford medial unicompartmental knee replacement.

The Oxford medial unicompartmental knee replacement was designed to reproduce normal mobility and forces in the knee, but its detailed effect on the patellofemoral joint has not been studied previously. We have examined the effect on patellofemoral mechanics of the knee by simultaneously measuring patellofemoral kinematics and forces in 11 cadaver knee specimens in a supine leg-extension rig. Comparison was made between the intact normal knee and sequential unicompartmental and total knee replacement. Following medial mobile-bearing unicompartmental replacement in 11 knees, patellofemoral kinematics and forces did not change significantly from those in the intact knee across any measured parameter. In contrast, following posterior cruciate ligament retaining total knee replacement in eight knees, there were significant changes in patellofemoral movement and forces. The Oxford device appears to produce near-normal patellofemoral mechanics, which may partly explain the low incidence of complications with the extensor mechanism associated with clinical use.

Simultaneous in vitro measurement of patellofemoral kinematics and forces.

This study involved the development and testing of a system for the simultaneous in vitro measurement of tibiofemoral kinematics and patellofemoral kinematics and forces. Knee motion was tracked using a Vicon 370, and patellofemoral force was measured using a six degree-of-freedom transducer based on the design of Singerman et al. Using this system, twelve knee specimens were tested in supine leg extension under a simulated quadriceps force. The comprehensive set of results corresponds well to the individual results of previous studies. The measurement system will be of value in assessing the effects of total knee arthroplasty on patellar function.

Identifying gait events without a force plate during running: a comparison of methods.

This paper presents a comparison of four different methods of identifying the times of foot-strike and toe-off during running based on gait marker trajectories. The event times predicted by the methods were compared to those identified using a force plate for both over-ground and treadmill running. The effect of using different threshold values for the detection of gait events using force plate data was also investigated, and as a result, all assessments of event detection accuracy were based on a cut-off value of 10N. The most accurate method of foot-strike detection depended on whether the runner landed with a rear- or a mid-foot strike. For rear-foot-strike running, the best method of identifying foot-strike used the vertical acceleration profile of the posterior heel distal marker and the vertical position profile of the hallux marker. For mid-foot-strike running, the best method of identifying foot-strike used the vertical velocity profile of the mean positions of the posterior heel distal marker and a marker midway between the second and third metatarsal heads. The most accurate method of identifying toe-off did not depend on type of foot-strike and was based on the vertical acceleration and position profiles of the hallux marker.

Correlation between lower limb bone morphology and gait characteristics in children with spastic diplegic cerebral palsy.

Background: Children with spastic diplegic cerebral palsy (CP) exhibit abnormal walking patterns and frequently develop lower limb, long bone deformities. It is important to determine if any relationship exists between bone morphology and movement of the lower limbs in children with CP. This is necessary to explain and possibly prevent the development of these deformities. Methods: This study investigated the relationship between bone morphology and gait characteristics in 10 healthy children (age range, 6-13 years; mean, 8 years 7 months; SD, ±2 years 7 months) and 9 children with spastic diplegic CP (age range, 6-12 years; mean, 9 years 2.5 months; SD, ±1 year 10.5 months) with no previous surgery. Three-dimensional magnetic resonance images were analyzed to define bone morphology. Morphological characteristics, such as the bicondylar angle, neck-shaft angle, anteversion angle, and tibial torsion, were measured. Gait analyses were performed to obtain kinematic characteristics of CP and normal childrens gait. Principal component analysis was used to reduce the dimensionality of 27 parameters (26 kinematics variables and age of the children) to 8 independent variables. Correlations between gait and bone morphology were determined for both groups of children. Results: Results indicated that in healthy children, hip adduction was correlated with neck-shaft and bicondylar angles. In CP children, pelvic obliquity correlated with neck-shaft angle, and foot rotation with bicondylar angle. In the transverse plane, hip and pelvic rotational kinematics were related to femoral anteversion in healthy children and to tibial torsion in CP children. Conclusion: Different development was observed in femoral and tibial morphology between CP and healthy children. The relationship between bone shape and dynamic gait patterns also varied between these populations. This needs to be taken into account, particularly when surgical treatment is planned. Clinical Relevance: Understanding the relationship between gait abnormality and bone deformity could eventually help in developing treatment regimens that will address gait deviations at the correct level and promote normal bone growth in children with CP.

Experimental validation of a finite element model of a composite tibia.

Abstract Composite bones are synthetic models made to simulate the mechanical behaviour of human bones. Finite element (FE) models of composite bone can be used to evaluate new and modified designs of joint prostheses and fixation devices. The aim of the current study was to create an FE model of a composite tibia and to validate it against results obtained from a comprehensive set of experiments. For this, 17 strain rosettes were attached to a composite tibia (model 3101, Pacific Research Laboratories, Vashon, Washington, USA). Surface strains and displacements were measured under 13 loading conditions. Two FE models were created on the basis of computed tomography scans. The models differed from each other in the mesh and material properties assigned. The experiments were simulated on them and the results compared with experimental results. The more accurate model was selected on the basis of regression analysis. In general, experimental strain measurements were highly repeatable and compared well with published results. The more accurate model, in which the inner elements representing the foam were assigned isotropic material properties and the elements representing the epoxy layer were assigned transversely isotropic material properties, was able to simulate the mechanical behaviour of the tibia with acceptable accuracy. The regression line for all axial loads combined had a slope of 0.999, an intercept of -6.24 microstrain, and an R2 value of 0.962. The root mean square error as a percentage was 5 per cent.