First-order model for the analysis of stump stresses for below-knee amputees

This study extends what had been a purely numerical model that used influence-factor matrices to relate the stump stresses to prosthesis loads for unilateral, below/knee amputees. Previously published influence-factor matrices are now factored into a coefficient matrix times the inverse of a stump geometry matrix. Using actual stump parameters, new information is learned about how the resistive moment of the stump balances the flexion-extension moment of the prosthesis and why certain normal stresses reach a maximum during a specific portion of the stance phase of the prosthesis.     

A comparison of passive flexion-extension to normal gait in the ovine stifle joint

Obtaining accurate values of joint tissue loads in human subjects and animals in vivo requires exact 3D-reproduction of joint kinematics and comparisons of in vivo motions between subjects and animals, and also necessitates an accurate reference position. For the knee, passive flexion-extension of isolated joints by hand has been assumed to produce bony motions similar to those of normal gait. We hypothesized that passive flexion-extension kinematics would not accurately reproduce in vivo gait, and, further, that such kinematics would vary significantly between testers. In vivo gait motions of four ovine stifle joints were measured in six degrees of freedom, as were passive flexion-extension motions after sacrifice. Passive flexion-extension motions were performed by three testers on the same stifle joints used in vitro. Results showed statistically significant differences in all degrees of freedom, with the largest differences in the proximal-distal and internal-external directions. Differences induced by muscle loads and kinetic factors in vivo were most evident during stance and hoof-off phases of gait. The in vitro passive paths generated by hand created motions with large variability both between and within individual testers. The user dependence and "area" of motion of passive flexion-extension indicates that passive flexion-extension is contained in a volume of motion, rather than constrained to a unique path. The assumption that the passive path has relevance to precise bone positions during normal in vivo gait is not supported by these results. Thus, using passive flexion-extension as a reference between joints may introduce large motion variability in the observed outcome, and large potential errors in determining joint tissue loads.     

A finite element study of the effect of diametral interface gaps on the contact areas and pressures in uncemented cylindrical femoral total hip components

Uncemented femoral total hip components rely entirely on contact with the prepared femur for their initial fixation. The contact areas and stresses between a straight tubular bone and a metal cylindrical prosthesis 12.5 cm long and 13 mm in diameter were calculated in a finite element model which includes uniform diametral gaps varying from 20 to 500 microns, using transverse loads from 100 to 2000 N. Frictionless three-dimensional contact elements were used between the bone and the prosthesis. Contact stresses were high and irregular in all cases, and the contact areas were small. Two regions of contact were apparent for lower loads and larger gaps. A third region of contact occurred near the distal tip of the implant at higher loads. This region of contact markedly increased the contact stresses at the distal tip of the prosthesis. A 20 microns overlap between bone and implant was modelled to assess a slight interference fit. The contact stress distribution in this case was markedly different from the stress distribution with a 20 microns diametral gap. The data collectively indicates that gaps of less than 20 microns between bone and implant can substantially change contact stress distributions.     

A three-dimensional non-linear finite element study of the effect of cement-prosthesis debonding in cemented femoral total hip components.

A three-dimensional non-linear finite element analysis of a cemented femoral component in which the component was partially debonded from the cement mantle was used to assess the effects of debonding on stresses in the cement. Three cases of partial cement-metal debonding were modelled with debonding of the proximal portion of the implant down to a horizontal plane which was 35, 62.5, or 82.5 mm below the prosthesis collar. Each situation was studied under loads simulating both gait and stairclimbing. Also, complete debonding between the implant and the surrounding cement mantle was modeled for loads simulating gait. Under stair climbing loads with partial cement-mental debonding, hoop stresses of 13-18 MPa were observed in the cement at the cement-metal interface at the proximal postero-medial corner of the implant. Similarly, in stair climbing, the maximum principal stresses in the cement were also adjacent to the proximal postero-medial region of the implant. These stresses were compressive and increased from 15 MPa with fully bonded interfaces to 48 MPa with debonding down to 82.5 mm below the prosthesis collar. Under gait loads, complete debonding caused high compressive stresses up to 34.9 MPa in the cement distal to the prosthesis tip. Thus, cement failure subsequent to prosthesis debonding is likely in the proximal region in a partially debonded implant due to stair climbing loads and is likely below the prosthesis tip in a fully debonded implant due to gait loading.     

Three-dimensional finite element analysis of glenoid replacement prostheses: a comparison of keeled and pegged anchorage systems

Glenoid component loosening is the dominant cause of failure in total shoulder arthroplasty. It is presumed that loosening in the glenoid is caused by high stresses in the cement layer. Several anchorage systems have been designed with the aim of reducing the loosening rate, the two major categories being "keeled" fixation and "pegged" fixation. However, no three-dimensional finite element analysis has been performed to quantify the stresses in the cement or to compare the different glenoid prosthesis anchorage systems. The objective of this study was to determine the stresses in the cement layer and surrounding bone for glenoid replacement components. A three-dimensional model of the scapula was generated using CT data for geometry and material property definition. Keeled and pegged designs were inserted into the glenoid, surrounded by a 1-mm layer of bone cement. A 90 deg arm abduction load with a full muscle and joint load was applied, following van der Helm (1994). Deformations of the prosthesis, stresses in the cement, and stresses in the bone were calculated. Stresses were also calculated for a simulated case of rheumatoid arthritis (RA) in which bone properties were modified to reflect that condition. A maximum principal stress-based failure model was used to predict what quantity of the cement is at risk of failure at the levels of stress computed. The prediction is that 94 percent (pegged prosthesis) and 68 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival in normal bone. In RA bone, however, the situation is reversed where 86 percent (pegged prosthesis) and 99 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival. Bone stresses are shown to be not much affected by the prosthesis design, except at the tip of the central peg or keel. It is concluded that a "pegged" anchorage system is superior for normal bone, whereas a "keeled" anchorage system is superior for RA bone.     

Stress analysis and failure prediction in the proximal femur before and after total hip replacement

An investigation was performed to determine the effects of the presence of two lengths of proximal Müller prosthesis on predicted failure loads, as compared to those for an intact femur. Three-dimensional stresses in a bone/cement/prosthesis system were determined using finite element methods, with both isotropic and transversely isotropic material properties used for the diaphyseal cortex. Significant increases in prosthesis stem stresses were found when the transversely isotropic material properties were employed in the diaphyseal cortex. This leads to the conclusion that accurate anisotropic material properties for bone are essential for precise stress determination and optimum design in prosthetic implants. Failure loads were also predicted for vertical compression and axial torque, similar to available experimental conditions, and were within the range of the experimental failure data found in the literature. The technique developed herein can be used to systematically assess existing as well as future implant designs, taking into account the complex three-dimensional interaction effects of the overall bone/cement/prosthesis system.     

Analytical description of minimum energy expenditure surfaces

Mechanical energy expenditure during level walking was evaluated and graphed for two unilateral, below-knee amputees over time and a range of adjustments of the flexion-extension alignment angle. The resulting mechanical energy surfaces were then least-squared fitted with an analytical function that was linear in time and quadratic in flexion-extension alignment angle. The least-squares analysis showed that there was a flexion-extension adjustment that minimized the mechanical energy expenditure and that this optimal adjustment was very close to the design point set by certified prosthetists.     

Load transfer mechanics between trans-tibial prosthetic socket and residual limb--dynamic effects

The effects of inertial loads on the interface stresses between trans-tibial residual limb and prosthetic socket were investigated. The motion of the limb and prosthesis was monitored using a Vicon motion analysis system and the ground reaction force was measured by a force platform. Equivalent loads at the knee joint during walking were calculated in two cases with and without consideration of the material inertia. A 3D nonlinear finite element (FE) model based on the actual geometry of residual limb, internal bones and socket liner was developed to study the mechanical interaction between socket and residual limb during walking. To simulate the friction/slip boundary conditions between the skin and liner, automated surface-to-surface contact was used. The prediction results indicated that interface pressure and shear stress had the similar double-peaked waveform shape in stance phase. The average difference in interface stresses between the two cases with and without consideration of inertial forces was 8.4% in stance phase and 20.1% in swing phase. The maximum difference during stance phase is up to 19%. This suggests that it is preferable to consider the material inertia effect in a fully dynamic FE model.     

Contact stresses and fatigue life in a knee prosthesis: comparison between in vitro measurements and computational simulations

The evaluation of contact areas and pressures in total knee prosthesis is a key issue to prevent early failure. The first part of this study is based on the hypothesis that the patterns of contact stresses on the tibial insert of a knee prosthesis at different stages of the gait cycle could be an indicator of the wear performances of a knee prosthesis. Contact stresses were calculated for a mobile bearing knee prosthesis by means of finite element method (FEM). Contact areas and stresses were also measured through in vitro tests using Fuji Prescale film in order to support the FEM findings.The second part of this study addresses the long-term structural integrity of metal tibial components in terms of fatigue life by means of experimental tests and FEM simulations. Fatigue experimental evaluations were performed on Cr-Co alloy tibial tray, based on ISO standards. FEM models were used to calculate the stress patterns. The failure risk was estimated with a standard fatigue criterion on the basis of the results obtained from the FEM calculations. Experimental and computational results showed a positive matching.     

FE-analysis of surface stresses for the tribological system in total hip prostheses]

The implantation of a total hip prosthesis is an operation which is performed frequently due to advanced hip joint damage both in humans and in veterinary medicine in dogs. The long-term result of a hip prosthesis is mainly determined by aseptic loosening of the prosthesis; among other causes, abrasion particles of the tribological pairing are responsible for the loosening. For the analysis of the surface stresses with different tribological pairings, a finite element model was generated which was based on the CAD data of a commercial total hip prosthesis. After transmission of a physiological force in the components of the three tribological pairings ceramic/polyethylene, ceramic/ceramic and metal/polyethylene, stresses were calculated. Stresses in the ceramic/ceramic tribological pairings were conspicuously higher than in the other material pairings. In the future adapted prostheses have to be developed that ensure optimal friction and absorption characteristics of the components.     

Effects of shear stress on articular chondrocyte metabolism

The articular cartilage of diarthrodial joints experiences a variety of stresses, strains and pressures that result from normal activities of daily living. In normal cartilage, the extracellular matrix exists as a highly organized composite of specialized macromolecules that distributes loads at the bony ends. The chondrocyte response to mechanical loading is recognized as an integral component in the maintenance of articular cartilage matrix homeostasis. With inappropriate mechanical loading of the joint, as occurs with traumatic injury, ligament instability, bony malalignment or excessive weight bearing, the cartilage exhibits manifestations characteristic of osteoarthritis. Breakdown of cartilage in osteoarthritis involves degradation of the extracellular matrix macromolecules and decreased expression of chondrocyte proteins necessary for normal joint function. Osteoarthritic cartilage often exhibits increased amounts of type I collagen and synthesis of proteoglycans characteristic of immature cartilage. The shift in cartilage phenotype in response to altered load yields a matrix that fails to support normal joint function. Mathematical modeling and experimental studies in animal models confirm an association between altered loading of diarthrotic joints and arthritic changes. Both types of studies implicate shear forces as a critical component in the destructive profile. The severity of cartilage destruction in response to altered loads appears linked to expression of biological factors influencing matrix integrity and cellular metabolism. Determining how shear stress alters chondrocyte metabolism is fundamental to understanding how to limit matrix destruction and stimulate cartilage repair and regeneration. At present, the precise biochemical and molecular mechanisms by which shear forces alter chondrocyte metabolism from a normal to a degenerative phenotype remain unclear. The results presented here address the hypothesis that articular chondrocyte metabolism is modulated by direct effects of shear forces that act on the cell through mechanotransduction processes. The purpose of this work is to develop critical knowledge regarding the basic mechanisms by which mechanical loading modulates cartilage metabolism in health and disease. This presentation will describe the effects of using fluid induced shear stress as a model system for stimulation of articular chondrocytes in vitro. The fluid induced shear stress was applied using a cone viscometer system to stimulate all the cells uniformly under conditions of minimal turbulence. The experiments were carried using high-density primary monolayer cultures of normal and osteoarthritic human and normal bovine articular chondrocytes. The analysis of the cellular response included quantification of cytokine release, matrix metalloproteinase expression and activation of intracellular signaling pathways. The data presented here show that articular chondrocytes exhibit a dose- and time-dependent response to shear stress that results in the release of soluble mediators and extracellular matrix macromolecules. The data suggest that the chondrocyte response to mechanical stimulation contributes to the maintenance of articular cartilage homeostasis in vivo.     

A new transducer for facet force measurement in the lumbar spine: benchmark and in vitro test results.

A new transducer capable of direct measurement of time-dependent loads in human lumbar facet joints was developed and tested. The transducer was comprised of a force-sensitive resistor (FSR) in series with a pressure-sensitive film. A wide range of experiments revealed the performance attributes and limitations of the FSR. The output signal of the FSR is actually sensitive to both force and area of contact independently. Therefore, a pressure-sensitive film was used to quantify the contact area. At least two transformation equations were calculated for each FSR corresponding to known contact areas. Each equation was a linearization of the log of the FSR output vs the log of the applied ramp loads. Coefficients of determination (CD) were calculated for small (21 mm2) and large (32 mm2) contact areas, and were found to exceed 0.900 for all data. The average of nine cycles was nearly linear for some FSRs (CD of 0.999). FSR output signal and contact area were recorded in cadaveric lumbar facets under ramp load. The appropriate transformation equation was determined by a linear interpolation between benchmark equations based on the contact area measured in vitro. Facet force measurements compared well with those of other researchers. The transducer was found to be quite easy to use.     

Biomechanics of the C5-C6 Spinal Unit Before and After Placement of a disc prosthesis

The study consists of a biomechanical comparison between the intact C5–C6 spinal segment and the same segment implanted with the Bryan^TM artificial disc prosthesis (Medtronic Ltd., Memphis, TN, USA), by the use of the finite element (FE) method. Our target is the prediction of the influence of prosthesis placement on the resulting mechanics of the C5–C6 spine unit. A FE model of the intact C5–C6 segment was built, employing realistic models of the vertebrae, disc and ligaments. Simulations were conducted imposing a compression preload combined to a flexion/extension moment, a pure lateral bending moment and a pure torsion moment, and the calculated results were compared to data from literature. The model was then modified to include the Bryan^TM cervical disc prosthesis, and the simulations were repeated. The location of the instantaneous center of rotation (ICR) of C5 with respect to C6 throughout flexion/extension was calculated in both models. In general, the moment–rotation curves obtained from the disc prosthesis-implanted model were comparable to the curves obtained from the intact model, except for a slightly greater stiffness induced by the artificial disc. The position of the calculated ICRs was rather stable throughout flexion-extension and was generally confined to a small area, qualitatively matching the corresponding physiological region, in both models. These results imply that the Bryan^TM disc prosthesis allows to correctly reproduce a physiological flexion/extension at the implanted level. The results of this study have quantified aspects that may assist in optimizing cervical disc replacement primarily from a biomechanical point of view.     

Simulation of fretting wear at orthopaedic implant interfaces

Osteolysis due to wear debris is a primary cause of failure of total joint replacements. Although debris produced by the joint articulating surfaces has been studied and simulated extensively, fretting wear debris, produced at nonarticulating surfaces, has not received adequate attention. We developed a three-station fretting wear simulator to reproduce in vivo motion and stresses at the interfaces of total joint replacements. The simulator is based on the beam bending theory and is capable of producing cyclic displacement from 3 to 1000 microns, under varying magnitudes of contact stresses. The simulator offers three potential advantages over previous studies: The ability to control the displacement by load, the ability to produce very small displacements, and dynamic normal loads as opposed to static. A pilot study was designed to test the functionality of the simulator, and verify that calculated displacements and loads produced the predicted differences between two commonly used porous ingrowth titanium alloy surfaces fretting against cortical bone. After 1.5 million cycles, the simulator functioned as designed, producing greater wear of bone against the rougher plasma-sprayed surface compared to the fiber-mesh surface, as predicted. A novel pin-on-disk apparatus for simulating fretting wear at orthopaedic implant interfaces due to micromotion is introduced. The test parameters measured with the fretting wear simulator were as predicted by design calculations, and were sufficient to measure differences in the height and weight of cortical bone pins rubbing against two porous ingrowth surfaces, plasma-sprayed titanium and titanium fiber mesh.     

Patient-Specific Design and Biomechanical Evaluation of a Novel Bipolar Femoral Hemi-Knee Prosthesis

While total knee replacement is successful, hemiarthroplasty is necessary for some young, obese and active patients who are especially not suitable for unicompartmental or total knee prostheses. Hemiarthroplasty also provides an opportunity for children with bone tumors. The design ofhemiarthroplasty should be patient-specific to reduce contact stress and friction as well as instability, compared to conventional hemi-knee prosthesis. A novel bipolar hemi-knee prosthesis with two flexion stages was developed according to a healthy male＇s knee morphological profile. The motion mode of the bipolar hemi-knee prosthesis was observed through roentgenoscopy in vitro experiment. The biomechanical properties in one gait cycle were evaluated though finite element simulation. The bipolar hemi-knee prosthesis was found to produce knee flexion at two stages through X-ray images. The first stage is the motion from upright posture to a specified 60~ flexion, followed by the second stage of motion subsequently to deep flexion. The finite element simulation results also show that the designed hemi-knee prosthesis has the ability to reduce stresses on the joint contact surfaces. Therefore, it is possible for the bipolar hemi-knee prosthesis to provide better biotribological performances because it can reduce stresses and potentially wear on the opposing contacting surface during a gait cycle, orovidin~ a t~romisin~ treatment strate~v in future Joint renair znd renlneement     

On the estimation of joint kinematics during gait.

In gait analysis, the concepts of Euler and helical (screw) angles are used to define the three-dimensional relative joint angular motion of lower extremities. Reliable estimation of joint angular motion depends on the accurate definition and construction of embedded axes within each body segment. In this paper, using sensitivity analysis, we quantify the effects of uncertainties in the definition and construction of embedded axes on the estimation of joint angular motion during gait. Using representative hip and knee motion data from normal subjects and cerebral palsy patients, the flexion-extension axis is analytically perturbed +/- 15 degrees in 5 degrees steps from a reference position, and the joint angles are recomputed for both Euler and helical angle definitions. For the Euler model, hip and knee flexion angles are relatively unaffected while the ab/adduction and rotation angles are significantly affected throughout the gait cycle. An error of 15 degrees in the definition of flexion-extension axis gives rise to maximum errors of 8 and 12 degrees for the ab/adduction angle, and 10-15 degrees for the rotation angles at the hip and knee, respectively. Furthermore, the magnitude of errors in ab/adduction and rotation angles are a function of the flexion angle. The errors for the ab/adduction angles increase with increasing flexion angle and for the rotation angle, decrease with increasing flexion angle. In cerebral palsy patients with flexed knee pattern of gait, this will result in distorted estimation of ab/adduction and rotation. For the helical model, similar results are obtained for the helical angle and associated direction cosines.(ABSTRACT TRUNCATED AT 250 WORDS)     

Microstructural Inhomogeneity of Electrical Conductivity in Subcutaneous Fat Tissue

Microscopic peculiarities stemming from a temperature increase in subcutaneous adipose tissue (sWAT) after applying a radio-frequency (RF) current, must be strongly dependent on the type of sWAT. This effect is connected with different electrical conductivities of pathways inside (triglycerides in adipocytes) and outside (extra-cellular matrix) the cells and to the different weighting of these pathways in hypertrophic and hyperplastic types of sWAT. The application of the RF current to hypertrophic sWAT, which normally has a strongly developed extracellular matrix with high concentrations of hyaluronan and collagen in a peri-cellular space of adipocytes, can produce, micro-structurally, a highly inhomogeneous temperature distribution, characterized by strong temperature gradients between the peri-cellular sheath of the extra-cellular matrix around the hypertrophic adipocytes and their volumes. In addition to normal temperature effects, which are generally considered in body contouring, these temperature gradients can produce thermo-mechanical stresses on the cells’ surfaces. Whereas these stresses are relatively small under normal conditions and cannot cause any direct fracturing or damage of the cell structure, these stresses can, under some supportive conditions, be theoretically increased by several orders of magnitude, causing the thermo-mechanical cell damage. This effect cannot be realized in sWAT of normal or hyperplastic types where the peri-cellular structures are under-developed. It is concluded that the results of RF application in body contouring procedures must be strongly dependent on the morphological structure of sWAT.     

Pelvic stresses in vitro--II. A study of the efficacy of metal-backed acetabular prostheses

Six embalmed human cadaveric hemi-pelves with their associated proximal femurs have been tested in vitro using 25 strain-gauge rosettes on each hemi-pelvis. Loads were applied up to 2.5 kN and principal stresses were computed from the principal strain data. Acetabular prostheses, cemented in place upon a cartilage-devoid but intact subchondral bone-plate, showed little change in stress-patterns when compared with the normal data, regardless of whether or not the component employed metal-backing. The use of 30 anchoring holes of 6.4 mm diameter, in the intact subchondral bone-plate, had little effect upon the stress-patterns, regardless of whether metal-backing was employed upon the prosthesis. When the subchondral bone-plate was removed, there were notable changes in the stress-pattern in the periacetabular region and on the medial wall of pelvis. The metal-backed prosthesis did not produce such notable changes as its plastic counterpart, when the subchondral bone-plate was removed. The use of a plastic prosthesis cemented in a Protrusio ring, in an acetabulum devoid of subchondral bone, produced notable changes in the stress-patterns in the whole periacetabular region and on the medial wall.     

Sagittal plane kinematics of the adult hyoid bone

The hyoid bone is a unique bone in the skeleton not articulated to any other bone. The hyoid muscles, which attach to the hyoid bone, may play a role in neck mechanics, but analysis of their function requires quantifying hyoid bone mechanics. The goal of this study was to obtain the detailed kinematics of the hyoid bone over a large range of flexion-extension motion using radiographs at 5 postures. The position of the hyoid bone in the sagittal plane was characterized with respect to head, jaw, and vertebral movements. Sex differences in hyoid kinematics were also investigated. We hypothesized that (1) the position of the hyoid bone in the sagittal plane is linearly correlated with motion of the head, jaw, and vertebrae, and (2) the hyoid position, size, and kinematics are sex-specific. We found that the hyoid bone X, Y, and angular position generally had strong linear correlations with the positions of the head, jaw, and the cervical vertebrae C1-C4. Hyoid X and angular position was also correlated to C5. Sex differences were found in some regressions of the hyoid bone with respect to C1-C5. The angular and linear measurements of the hyoid bone showed sex differences in absolute values, which were not evident after normalization by posture or neck size. Incorporating these results to neck models would enable accurate modeling of the hyoid muscles. This may have implications for analyzing the mechanics of the cervical spine, including loads on neck structures and implants.