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.

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.

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.

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.

Exploration of the defect’s effect on the mechanical properties of different orientated nanowires

Molecular dynamics (MD) simulations have been carried out to investigate the defect’s effect on the mechanical properties of copper nanowire with different crystallographic orientations, under tensile deformation. Three different crystallographic orientations have been considered. The deformation mechanism has been carefully discussed. It is found that the Young’s modulus is insensitive to the defect, even when the nanowire’s crystallographic orientation is different. However, due to the defect’s effect, the yield strength and yield strain appear a large decrease. The defects have played a role of dislocation sources, the slips or stacking faults are first generated around the locations of the defects. The necking locations have also been affected by different defects. Due to the surface defect, the plastic deformation has received a large influence for the <001>/{110} and <110> orientated nanowires, and a relative small influence is seen for the <111> nanowire.

Advanced numerical characterization of mono-crystalline copper with defects

Based on the embedded atom method (EAM) and molecular dynamics (MD) method, the mono-crystalline copper with different defects is investigated through tension and nanoindentation simulation. The single-crystal copper nanowire with surface defects is firstly studied through tension. For validation, the tension simulations for nanowire without defect are carried out under different temperatures and strain rates. The defects on nanowires are then systematically studied in considering different defects orientation distribution. It is found that the Young’s modulus is insensitive of surface defects and centro-plane defects. However, the yield strength and yield point show a significant decrease due to the different defects. Specially, the 〖45〗^° defect in surface and in (200) plane exerts the biggest influence to the yield strength, about 34.20% and 51.45% decrease are observed, respectively. Different defects are observed to serve as a dislocation source and different necking positions of the nanowires during tension are found. During nanoindentation simulation, dislocation is found nucleating below the contact area, but no obvious dislocation is generated around the nano-cavity. Comparing with the perfect substrate during nanoindentation, the substrate with nano-cavities emerged less dislocations, it is supposed that the nano-cavity absorbed part of the indent energy, and less plastic deformation happened in the defected substrate.

Structural and thermal performance of cold-formed steel stud wall systems under fire conditions

Cold-formed steel stud walls are a major component of Light Steel Framing (LSF) building systems used in commercial, industrial and residential buildings. In the conventional LSF stud wall systems, thin steel studs are protected from fire by placing one or two layers of plasterboard on both sides with or without cavity insulation. However, there is very limited data about the structural and thermal performance of stud wall systems while past research showed contradicting results, for example, about the benefits of cavity insulation. This research was therefore conducted to improve the knowledge and understanding of the structural and thermal performance of cold-formed steel stud wall systems (both load bearing and non-load bearing) under fire conditions and to develop new improved stud wall systems including reliable and simple methods to predict their fire resistance rating. Full scale fire tests of cold-formed steel stud wall systems formed the basis of this research. This research proposed an innovative LSF stud wall system in which a composite panel made of two plasterboards with insulation between them was used to improve the fire rating. Hence fire tests included both conventional steel stud walls with and without the use of cavity insulation and the new composite panel system. A propane fired gas furnace was specially designed and constructed first. The furnace was designed to deliver heat in accordance with the standard time temperature curve as proposed by AS 1530.4 (SA, 2005). A compression loading frame capable of loading the individual studs of a full scale steel stud wall system was also designed and built for the load-bearing tests. Fire tests included comprehensive time-temperature measurements across the thickness and along the length of all the specimens using K type thermocouples. They also included the measurements of load-deformation characteristics of stud walls until failure. The first phase of fire tests included 15 small scale fire tests of gypsum plasterboards, and composite panels using different types of insulating material of varying thickness and density. Fire performance of single and multiple layers of gypsum plasterboards was assessed including the effect of interfaces between adjacent plasterboards on the thermal performance. Effects of insulations such as glass fibre, rock fibre and cellulose fibre were also determined while the tests provided important data relating to the temperature at which the fall off of external plasterboards occurred. In the second phase, nine small scale non-load bearing wall specimens were tested to investigate the thermal performance of conventional and innovative steel stud wall systems. Effects of single and multiple layers of plasterboards with and without vertical joints were investigated. The new composite panels were seen to offer greater thermal protection to the studs in comparison to the conventional panels. In the third phase of fire tests, nine full scale load bearing wall specimens were tested to study the thermal and structural performance of the load bearing wall assemblies. A full scale test was also conducted at ambient temperature. These tests showed that the use of cavity insulation led to inferior fire performance of walls, and provided good explanations and supporting research data to overcome the incorrect industry assumptions about cavity insulation. They demonstrated that the use of insulation externally in a composite panel enhanced the thermal and structural performance of stud walls and increased their fire resistance rating significantly. Hence this research recommends the use of the new composite panel system for cold-formed LSF walls. This research also included steady state tensile tests at ambient and elevated temperatures to address the lack of reliable mechanical properties for high grade cold-formed steels at elevated temperatures. Suitable predictive equations were developed for calculating the yield strength and elastic modulus at elevated temperatures. In summary, this research has developed comprehensive experimental thermal and structural performance data for both the conventional and the proposed non-load bearing and load bearing stud wall systems under fire conditions. Idealized hot flange temperature profiles have been developed for non-insulated, cavity insulated and externally insulated load bearing wall models along with suitable equations for predicting their failure times. A graphical method has also been proposed to predict the failure times (fire rating) of non-load bearing and load bearing walls under different load ratios. The results from this research are useful to both fire researchers and engineers working in this field. Most importantly, this research has significantly improved the knowledge and understanding of cold-formed LSF walls under fire conditions, and developed an innovative LSF wall system with increased fire rating. It has clearly demonstrated the detrimental effects of using cavity insulation, and has paved the way for Australian building industries to develop new wall panels with increased fire rating for commercial applications worldwide.

A nanoscratch method for measuring hardness of thin films

Thin solid films were extensively used in the making of solar cells, cutting tools, magnetic recording devices, etc. As a result, the accurate measurement of mechanical properties of the thin films, such as hardness and elastic modulus, was required. The thickness of thin films normally varies from tens of nanometers to several micrometers. It is thus challenging to measure their mechanical properties. In this study, a nanoscratch method was proposed for hardness measurement. A three-dimensional finite element method (3-D FEM) model was developed to validate the nanoscratch method and to understand the substrate effect during nanoscratch. Nanoindentation was also used for comparison. The nanoscratch method was demonstrated to be valuable for measuring hardness of thin solid films.

Atomistic exploration of deformation properties of copper nanowires with pre-existing defects

Based on the embedded atom method (EAM) and molecular dynamics (MD) method, in this paper, the tensile deformation properties of Cu nanowires (NWs) with different pre-existing defects, including single surface defects, surface bi-defects and single internal defects, are systematically studied. In-depth deformation mechanisms of NWs with pre-existing defects are also explored. It is found that Young's modulus is insensitive to different pre-existing defects, but yield strength shows an obvious decrease. Defects are observed influencing greatly on NWs' tensile deformation mechanisms, and playing a role of dislocation sources. Besides of the traditional deformation process dominated by the nucleation and propagation of partial dislocations, the generations of twins, grain boundaries, fivefold deformation twins, hexagonal close-packed (HCP) structure and phase transformation from face-centred cubic (FCC) structure to HCP structure have been triggered by pre-existing defects. It is found that surface defect intends to induce larger influence to yield strength than internal defect. Most importantly, the defect that lies on slip planes exerts larger influence than other defects. As expected, it is also found that the more or longer of the defect, the bigger influence will be induced.

Theoretical and numerical investigation of bending properties of Cu nanowires

Based on the AFM-bending experiments, a molecular dynamics (MD) bending simulation model is established which could accurately account for the full spectrum of the mechanical properties of NWs in a double clamped beam configuration, ranging from elasticity to plasticity and failure. It is found that, loading rate exerts significant influence to the mechanical behaviours of nanowires (NWs). Specifically, a loading rate lower than 10 m/s is found reasonable for a homogonous bending deformation. Both loading rate and potential between the tip and the NW are found to play an important role in the adhesive phenomenon. The force versus displacement (F-d) curve from MD simulation is highly consistent in shapes with that from experiments. Symmetrical F-d curves during loading and unloading processes are observed, which reveal the linear-elastic and non-elastic bending deformation of NWs. The typical bending induced tensile-compressive features are observed. Meanwhile, the simulation results are excellently fitted by the classical Euler-Bernoulli beam theory with axial effect. It is concluded that, axial tensile force becomes crucial in bending deformation when the beam size is down to nanoscale for double clamped NWs. In addition, we find shorter NWs will have an earlier yielding and a larger yielding force. Mechanical properties (Young’s modulus & yield strength) obtained from both bending and tensile deformations are found comparable with each other. Specifically, the modulus is essentially similar under these two loading methods, while the yield strength during bending is observed larger than that during tension.

Effective shear modulus approach for two dimensional solids and plate bending problems by meshless point collocation method

For the analysis of material nonlinearity, an effective shear modulus approach based on the strain control method is proposed in this paper by using point collocation method. Hencky’s total deformation theory is used to evaluate the effective shear modulus, Young’s modulus and Poisson’s ratio, which are treated as spatial field variables. These effective properties are obtained by the strain controlled projection method in an iterative manner. To evaluate the second order derivatives of shape function at the field point, the radial basis function (RBF) in the local support domain is used. Several numerical examples are presented to demonstrate the efficiency and accuracy of the proposed method and comparisons have been made with analytical solutions and the finite element method (ABAQUS).

MD investigations for mechanical properties of copper nanowires with and without surface defects

Based on the molecular dynamics (MD) method, the single-crystalline copper nanowire with different surface defects is investigated through tension simulation. For comparison, the MD tension simulations of perfect nanowire are firstly carried out under different temperatures, strain rates, and sizes. It has concluded that the surface-volume ratio significantly affects the mechanical properties of nanowire. The surface defects on nanowires are then systematically studied in considering different defect orientation and distribution. It is found that the Young’s modulus is insensitive of surface defects. However, the yield strength and yield point show a significant decrease due to the different defects. Different defects are observed to serve as a dislocation source.