The feasibility of vibration analysis as a technique to detect osseointegration of transfemoral implants

One of the main causes of above knee or transfemoral amputation (TFA) in the developed world is trauma to the limb. The number of people undergoing TFA due to limb trauma, particularly due to war injuries, has been increasing. Typically the trauma amputee population, including war-related amputees, are otherwise healthy, active and desire to return to employment and their usual lifestyle. Consequently there is a growing need to restore long-term mobility and limb function to this population. Traditionally transfemoral amputees are provided with an artificial or prosthetic leg that consists of a fabricated socket, knee joint mechanism and a prosthetic foot. Amputees have reported several problems related to the socket of their prosthetic limb. These include pain in the residual limb, poor socket fit, discomfort and poor mobility. Removing the socket from the prosthetic limb could eliminate or reduce these problems. A solution to this is the direct attachment of the prosthesis to the residual bone (femur) inside the residual limb. This technique has been used on a small population of transfemoral amputees since 1990. A threaded titanium implant is screwed in to the shaft of the femur and a second component connects between the implant and the prosthesis. A period of time is required to allow the implant to become fully attached to the bone, called osseointegration (OI), and be able to withstand applied load; then the prosthesis can be attached. The advantages of transfemoral osseointegration (TFOI) over conventional prosthetic sockets include better hip mobility, sitting comfort and prosthetic retention and fewer skin problems on the residual limb. However, due to the length of time required for OI to progress and to complete the rehabilitation exercises, it can take up to twelve months after implant insertion for an amputee to be able to load bear and to walk unaided. The long rehabilitation time is a significant disadvantage of TFOI and may be impeding the wider adoption of the technique. There is a need for a non-invasive method of assessing the degree of osseointegration between the bone and the implant. If such a method was capable of determining the progression of TFOI and assessing when the implant was able to withstand physiological load it could reduce the overall rehabilitation time. Vibration analysis has been suggested as a potential technique: it is a non destructive method of assessing the dynamic properties of a structure. Changes in the physical properties of a structure can be identified from changes in its dynamic properties. Consequently vibration analysis, both experimental and computational, has been used to assess bone fracture healing, prosthetic hip loosening and dental implant OI with varying degrees of success. More recently experimental vibration analysis has been used in TFOI. However further work is needed to assess the potential of the technique and fully characterise the femur-implant system. The overall aim of this study was to develop physical and computational models of the TFOI femur-implant system and use these models to investigate the feasibility of vibration analysis to detect the process of OI. Femur-implant physical models were developed and manufactured using synthetic materials to represent four key stages of OI development (identified from a physiological model), simulated using different interface conditions between the implant and femur. Experimental vibration analysis (modal analysis) was then conducted using the physical models. The femur-implant models, representing stage one to stage four of OI development, were excited and the modal parameters obtained over the range 0-5kHz. The results indicated the technique had limited capability in distinguishing between different interface conditions. The fundamental bending mode did not alter with interfacial changes. However higher modes were able to track chronological changes in interface condition by the change in natural frequency, although no one modal parameter could uniquely distinguish between each interface condition. The importance of the model boundary condition (how the model is constrained) was the key finding; variations in the boundary condition altered the modal parameters obtained. Therefore the boundary conditions need to be held constant between tests in order for the detected modal parameter changes to be attributed to interface condition changes. A three dimensional Finite Element (FE) model of the femur-implant model was then developed and used to explore the sensitivity of the modal parameters to more subtle interfacial and boundary condition changes. The FE model was created using the synthetic femur geometry and an approximation of the implant geometry. The natural frequencies of the FE model were found to match the experimental frequencies within 20% and the FE and experimental mode shapes were similar. Therefore the FE model was shown to successfully capture the dynamic response of the physical system. As was found with the experimental modal analysis, the fundamental bending mode of the FE model did not alter due to changes in interface elastic modulus. Axial and torsional modes were identified by the FE model that were not detected experimentally; the torsional mode exhibited the largest frequency change due to interfacial changes (103% between the lower and upper limits of the interface modulus range). Therefore the FE model provided additional information on the dynamic response of the system and was complementary to the experimental model. The small changes in natural frequency over a large range of interface region elastic moduli indicated the method may only be able to distinguish between early and late OI progression. The boundary conditions applied to the FE model influenced the modal parameters to a far greater extent than the interface condition variations. Therefore the FE model, as well as the experimental modal analysis, indicated that the boundary conditions need to be held constant between tests in order for the detected changes in modal parameters to be attributed to interface condition changes alone. The results of this study suggest that in a clinical setting it is unlikely that the in vivo boundary conditions of the amputated femur could be adequately controlled or replicated over time and consequently it is unlikely that any longitudinal change in frequency detected by the modal analysis technique could be attributed exclusively to changes at the femur-implant interface. Therefore further development of the modal analysis technique would require significant consideration of the clinical boundary conditions and investigation of modes other than the bending modes.

Effects of opening in the cross girders of the flat bottom rail wagons to the load transferring mechanisms when used as road bridge deck

A Flat Bed Rail Wagon (FBRW) has been proposed as an alternative solution for replacing bridges on low traffic volume roads. The subject matter for this paper is to investigate the impediment to load transfer from cross girders to main girder, through visually identifiable structural flaws. Namely, the effect of having large openings at close proximity to the connection of the main girder to the cross girder of a FBRW was examined. It was clear that openings locally reduce the section modulus of the secondary members; however it was unclear how these reductions would affect the load transfer to the main girder. The results are presented through modeling grillage action for which the loads applied onto the FBRW were distributed through cross girders to the main girder.

Atomistic numerical investigation of single-crystal copper nanowire with surface defects

Based on the embedded atom method (EAM), a molecular dynamics (MD) simulation is performed to study the single-crystal copper nanowire with surface defects through tension. The tension simulations for nanowire without defect are first carried out under different temperatures, strain rates and time steps and then surface defect effects for nanowire are investigated. The stress-strain curves obtained by the MD simulations of various strain rates show a rate below 1 x 10(9) s-1 will exert less effect on the yield strength and yield point, and the Young's modulus is independent of strain rate. a time step below 5 fs is recommend for the atomic model during the MD simulation. It is observed that high temperature leads to low Young's modulus, as well as the yield strength. The surface defects on nanowires are systematically studied in considering different defect orientations. It is found that the surface defect serves as a dislocation source, and the yield strength shows 34.20% decresse with 45 degree surface defect. Both yield strength and yield point are significantly influenced by the surface defects, except the Young's modulus.

A coupled meshless technique/molecular dynamics approach for deformation characterization of monocrystalline metal

This paper presents a multiscale study using the coupled Meshless technique/Molecular Dynamics (M2) for exploring the deformation mechanism of mono-crystalline metal (focus on copper) under uniaxial tension. In M2, an advanced transition algorithm using transition particles is employed to ensure the compatibility of both displacements and their gradients, and an effective local quasi-continuum approach is also applied to obtain the equivalent continuum strain energy density based on the atomistic poentials and Cauchy-Born rule. The key parameters used in M2 are firstly investigated using a benchmark problem. Then M2 is applied to the multiscale simulation for a mono-crystalline copper bar. It has found that the mono-crystalline copper has very good elongation property, and the ultimate strength and Young's modulus are much higher than those obtained in macro-scale.

Plantar enthesopathy : thickening of the enthesis is correlated with energy dissipation of the plantar fat pad during walking

Background: The enthesis of the plantar fascia is thought to play an important role in stress dissipation. However, the potential link between entheseal thickening characteristic of enthesopathy and the stress-dissipating properties of the intervening plantar fat pad have not been investigated. Purpose: This study was conducted to identify whether plantar fat pad mechanics explain variance in the thickness of the fascial enthesis in individuals with and without plantar enthesopathy. Study Design: Case-control study; Level of evidence, 3. Methods: The study population consisted of 9 patients with unilateral plantar enthesopathy and 9 asymptomatic, individually matched controls. The thickness of the enthesis of the symptomatic, asymptomatic, and a matched control limb was acquired using high-resolution ultrasound. The compressive strain of the plantar fat pad during walking was estimated from dynamic lateral radiographs acquired with a multifunction fluoroscopy unit. Peak compressive stress was simultaneously acquired via a pressure platform. Principal viscoelastic parameters were estimated from subsequent stress-strain curves. Results: The symptomatic fascial enthesis (6.7 ± 2.0 mm) was significantly thicker than the asymptomatic enthesis (4.2 ± 0.4 mm), which in turn was thicker than the enthesis (3.3 ± 0.4 mm) of control limbs (P < .05). There was no significant difference in the mean thickness, peak stress, peak strain, or secant modulus of the plantar fat pad between limbs. However, the energy dissipated by the fat pad during loading and unloading was significantly lower in the symptomatic limb (0.55 ± 0.17) when compared with asymptomatic (0.69 ± 0.13) and control (0.70 ± 0.09) limbs (P < .05). The sonographic thickness of the enthesis was correlated with the energy dissipation ratio of the plantar fat pad (r = .72, P < .05), but only in the symptomatic limb. Conclusion: The energy-dissipating properties of the plantar fat pad are associated with the sonograpic appearance of the enthesis in symptomatic limbs, providing a previously unidentified link between the mechanical behavior of the plantar fat pad and enthesopathy.

Reinforced hydrogels for silicone copolymer delivery for scar remediation

Hypertrophic scars are formed by collagen overproduction in wounded areas and often occur in victims of severe burns. There are several methods for hypertrophic scar remediation and silicone gel therapy is one of the more successful methods. Research by others has shown that the activity of these gels may be due to migration of amphiphilic silicone oligomers from the gel and into the dermis, down-regulating production of collagen by fibroblasts. Normal silicone oil (PDMS) does not produce the same effect on fibroblasts. The main purpose of this project is the introduction of a particular amphiphilic silicone rake copolymer into an appropriate network which can absorb and release the silicone copolymer on the scarred area. Hydrogels are polymeric crosslinked networks which can swell in water or a drug solution, and gradually release the drug when applied to the skin. The application of gel enhances the effectiveness of the therapy, reduces the period of treatment and can be comfortable for patients to use. Polyethylene glycol (PEG) based networks have been applied in this research, because the amphiphilic silicone rake copolymer to be used as a therapy has polyethylene oxide (PEO) as a side chain. These PEO side chains have very similar chemical structure to a PEG gel chain so enhancing both the compatibility and the diffusion of the amphiphilic silicone rake copolymer into and out of the gel. Synthesis of PEG-based networks has been performed by two methods: in situ silsesquioxane formation as crosslink with a sol-gel reaction under different conditions and UV curing. PEG networks have low mechanical properties which is a fundamental limitation of the polymer backbone. For mechanical properties enhancement, composite networks were synthesized using nano-silica with different surface modification. The chemical structure of in situ silsesquioxane in the dry network has been examined by Solid State NMR, Differential Scanning Calorimetry (DSC) and swelling measurements in water. Mechanical properties of dry networks were tested by Dynamic Mechanical Thermal Analysis (DMTA) to determine modulus and interfacial interaction between silica and the network. In this way a family of self-reinforced networks has been produced that have been shown to absorb and deliver the active amphiphilic silicone- PEO rake copolymer.

The mechanical heterogeneity of the hard callus influences local tissue strains during bone healing : a finite element study based on sheep experiments

During secondary fracture healing, various tissue types including new bone are formed. The local mechanical strains play an important role in tissue proliferation and differentiation. To further our mechanobiological understanding of fracture healing, a precise assessment of local strains is mandatory. Until now, static analyses using Finite Elements (FE) have assumed homogenous material properties. With the recent quantification of both the spatial tissue patterns (Vetter et al., 2010) and the development of elastic modulus of newly formed bone during healing (Manjubala et al., 2009), it is now possible to incorporate this heterogeneity. Therefore, the aim of this study is to investigate the effect of this heterogeneity on the strain patterns at six successive healing stages. The input data of the present work stemmed from a comprehensive cross-sectional study of sheep with a tibial osteotomy (Epari et al., 2006). In our FE model, each element containing bone was described by a bulk elastic modulus, which depended on both the local area fraction and the local elastic modulus of the bone material. The obtained strains were compared with the results of hypothetical FE models assuming homogeneous material properties. The differences in the spatial distributions of the strains between the heterogeneous and homogeneous FE models were interpreted using a current mechanobiological theory (Isakson et al., 2006). This interpretation showed that considering the heterogeneity of the hard callus is most important at the intermediate stages of healing, when cartilage transforms to bone via endochondral ossification.

Biaxial cell stimulation : a mechanical validation

To analyse mechanotransduction resulting from tensile loading under defined conditions, various devices for in vitro cell stimulation have been developed. This work aimed to determine the strain distribution on the membrane of a commercially available device and its consistency with rising cycle numbers, as well as the amount of strain transferred to adherent cells. The strains and their behaviour within the stimulation device were determined using digital image correlation (DIC). The strain transferred to cells was measured on eGFP-transfected bone marrow-derived cells imaged with a fluorescence microscope. The analysis was performed by determining the coordinates of prominent positions on the cells, calculating vectors between the coordinates and their length changes with increasing applied tensile strain. The stimulation device was found to apply homogeneous (mean of standard deviations approx. 2% of mean strain) and reproducible strains in the central well area. However, on average, only half of the applied strain was transferred to the bone marrow-derived cells. Furthermore, the strain measured within the device increased significantly with an increasing number of cycles while the membrane's Young's modulus decreased, indicating permanent changes in the material during extended use. Thus, strain magnitudes do not match the system readout and results require careful interpretation, especially at high cycle numbers.

Spatial and temporal variations of mechanical properties and mineral content of the external callus during bone healing

After bone fracture, various cellular activities lead to the formation of different tissue types, which form the basis for the process of secondary bone healing. Although these tissues have been quantified by histology, their material properties are not well understood. Thus, the aim of this study is to correlate the spatial and temporal variations in the mineral content and the nanoindentation modulus of the callus formed via intramembranous ossification over the course of bone healing. Midshaft tibial samples from a sheep osteotomy model at time points of 2, 3, 6 and 9 weeks were employed. PMMA embedded blocks were used for quantitative back scattered electron imaging and nanoindentation of the newly formed periosteal callus near the cortex. The resulting indentation modulus maps show the heterogeneity in the modulus in the selected regions of the callus. The indentation modulus of the embedded callus is about 6 GPa at the early stage. At later stages of mineralization, the average indentation modulus reaches 14 GPa. There is a slight decrease in average indentation modulus in regions distant to the cortex, probably due to remodelling of the peripheral callus. The spatial and temporal distribution of mineral content in the callus tissue also illustrates the ongoing remodelling process observed from histological analysis. Most interestingly the average indentation modulus, even at 9 weeks, remains as low as 13 GPa, which is roughly 60% of that for cortical sheep bone. The decreased indentation modulus in the callus compared to cortex is due to the lower average mineral content and may be perhaps also due to the properties of the organic matrix which might be different from normal bone.

Mechanical Behavior of Articular Cartilage After Osteochondral Autograft Transfer in an Ovine Model

BACKGROUND: Grafting of autologous hyaline cartilage and bone for articular cartilage repair is a well-accepted technique. Although encouraging midterm clinical results have been reported, no information on the mechanical competence of the transplanted joint surface is available. HYPOTHESIS: The mechanical competence of osteochondral autografts is maintained after transplantation. STUDY DESIGN: Controlled laboratory study. METHODS: Osteochondral defects were filled with autografts (7.45 mm in diameter) in one femoral condyle in 12 mature sheep. The ipsilateral femoral condyle served as the donor site, and the resulting defect (8.3 mm in diameter) was left empty. The repair response was examined after 3 and 6 months with mechanical and histologic assessment and histomorphometric techniques. RESULTS: Good surface congruity and plug placement was achieved. The Young modulus of the grafted cartilage significantly dropped to 57.5% of healthy tissue after 3 months (P < .05) but then recovered to 82.2% after 6 months. The aggregate and dynamic moduli behaved similarly. The graft edges showed fibrillation and, in some cases (4 of 6), hypercellularity and chondrocyte clustering. Subchondral bone sclerosis was observed in 8 of 12 cases, and the amount of mineralized bone in the graft area increased from 40% to 61%. CONCLUSIONS: The mechanical quality of transplanted cartilage varies considerably over a short period of time, potentially reflecting both degenerative and regenerative processes, while histologically signs of both cartilage and bone degeneration occur. CLINICAL RELEVANCE: Both the mechanically degenerative and restorative processes illustrate the complex progression of regeneration after osteochondral transplantation. The histologic evidence raises doubts as to the long-term durability of the osteochondral repair.

Long-term effects of hydrogel properties on human chondrocyte behavior

Hydrogels provide a 3-dimensional network for embedded cells and offer promise for cartilage tissue engineering applications. Nature-derived hydrogels, including alginate, have been shown to enhance the chondrocyte phenotype but are variable and not entirely controllable. Synthetic hydrogels, including polyethylene glycol (PEG)-based matrices, have the advantage of repeatability and modularity; mechanical stiffness, cell adhesion, and degradability can be altered independently. In this study, we compared the long-term in vitro effects of different hydrogels (alginate and Factor XIIIa-cross-linked MMP-sensitive PEG at two stiffness levels) on the behavior of expanded human chondrocytes and the development of construct properties. Monolayer-expanded human chondrocytes remained viable throughout culture, but morphology varied greatly in different hydrogels. Chondrocytes were characteristically round in alginate but mostly spread in PEG gels at both concentrations. Chondrogenic gene (COL2A1, aggrecan) expression increased in all hydrogels, but alginate constructs had much higher expression levels of these genes (up to 90-fold for COL2A1), as well as proteoglycan 4, a functional marker of the superficial zone. Also, chondrocytes expressed COL1A1 and COL10A1, indicative of de-differentiation and hypertrophy. After 12 weeks, constructs with lower polymer content were stiffer than similar constructs with higher polymer content, with the highest compressive modulus measured in 2.5% PEG gels. Different materials and polymer concentrations have markedly different potency to affect chondrocyte behavior. While synthetic hydrogels offer many advantages over natural materials such as alginate, they must be further optimized to elicit desired chondrocyte responses for use as cartilage models and for development of functional tissue-engineered articular cartilage.

Cell density alters matrix accumulation in two distinct fractions and the mechanical integrity of alginate-chondrocyte constructs

Chondrocyte density in articular cartilage is known to change with the development and growth of the tissue and may play an important role in the formation of a functional extracellular matrix (ECM). The objective of this study was to determine how initial chondrocyte density in an alginate hydrogel affects the matrix composition, its distribution between the cell-associated (CM) and further removed matrix (FRM) fractions, and the tensile mechanical properties of the developing engineered cartilage. Alginate constructs containing primary bovine chondrocytes at densities of 0, 4, 16, and 64 million cells/ml were fabricated and cultured for 1 or 2 weeks, at which time structural, biochemical, and mechanical properties were analyzed. Both matrix content and distribution varied with the initial cell density. Increasing cell density resulted in an increasing content of collagen and sulfated-glycosaminoglycan (GAG) and an increasing proportion of these molecules localized in the CM. While the equilibrium tensile modulus of cell-free alginate did not change with time in culture, the constructs with highest cell density were 116% stiffer than cell-free controls after 2 weeks of culture. The equilibrium tensile modulus was positively correlated with total collagen (r2 = 0.47, p < 0.001) and GAG content (r2 = 0.68, p < 0.001), and these relationships were enhanced when analyzing only those matrix molecules in the CM fraction (r2 = 0.60 and 0.72 for collagen and GAG, respectively, each p < 0.001). Overall, the results of this study indicate that initial cell density has a considerable effect on the developing composition, structure, and function of alginate–chondrocyte constructs.

Depth-dependent biomechanical and biochemical properties of fetal, newborn, and tissue-engineered articular cartilage

Adult articular cartilage has depth-dependent mechanical and biochemical properties which contribute to zone-specific functions. The compressive moduli of immature cartilage and tissue-engineered cartilage are known to be lower than those of adult cartilage. The objective of this study was to determine if such tissues exhibit depth-dependent compressive properties, and how these depth-varying properties were correlated with cell and matrix composition of the tissue. The compressive moduli of fetal and newborn bovine articular cartilage increased with depth (p < 0.05) by a factor of 4-5 from the top 0.1 mm (28 +/- 13 kPa, 141 +/- 10 kPa, respectively) to 1 mm deep into the tissue. Likewise, the glycosaminoglycan and collagen content increased with depth (both p < 0.001), and correlated with the modulus (both p < 0.01). In contrast, tissue-engineered cartilage formed by either layering or mixing cells from the superficial and middle zone of articular cartilage exhibited similarly soft regions at both construct surfaces, as exemplified by large equilibrium strains. The properties of immature cartilage may provide a template for developing tissue-engineered cartilage which aims to repair cartilage defects by recapitulating the natural development and growth processes. These results suggest that while depth-dependent properties may be important to engineer into cartilage constructs, issues other than cell heterogeneity must be addressed to generate such tissues. (c) 2005 Elsevier Ltd. All rights reserved.

The effect of friction on indenter force and pile-up in numerical simulations of bone nanoindentation

Nanoindentation is a useful technique for probing the mechanical properties of bone, and finite element (FE) modeling of the indentation allows inverse determination of elasto-plastic constitutive properties. However, all but one FE study to date have assumed frictionless contact between indenter and bone. The aim of this study was to explore the effect of friction in simulations of bone nanoindentation. Two dimensional axisymmetric FE simulations were performed using a spheroconical indenter of tip radius 0.6 m and angle 90°. The coefficient of friction between indenter and bone was varied between 0.0 (frictionless) and 0.3. Isotropic linear elasticity was used in all simulations, with bone elastic modulus E=13.56GPa and Poisson‟s ratio f 0.3. Plasticity was incorporated using both Drucker-Prager and von Mises yield surfaces. Friction had a modest effect on the predicted force-indentation curve for both von Mises and Drucker-Prager plasticity, reducing maximum indenter displacement by 10% and 20% respectively as friction coefficient was increased from zero to 0.3 (at a maximum indenter force of 5mN). However, friction has a much greater effect on predicted pile-up after indentation, reducing predicted pile-up from 0.27 to 0.11 m with a von Mises model, and from 0.09 to 0.02 m with Drucker-Prager plasticity. We conclude that it is potentially important to include friction in nanoindentation simulations of bone if pile-up is used to compare simulation results with experiment.