A meshfree-based local Galerkin method with condensation of degree of freedom for elastic dynamic analysis

Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dynamic analysis. In the present method, scattered nodes without connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approximation. Then local discrete equations can be simplified by condensation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by assembling all local discrete equations and are solved by using the standard implicit Newmark’s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is implemented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.

Porous scaffold architecture guides tissue formation

Critical-sized bone defect regeneration is a remaining clinical concern. Numerous scaffold-based strategies are currently being investigated to enable in vivo bone defect healing. However, a deeper understanding of how a scaffold influences the tissue formation process and how this compares to endogenous bone formation or to regular fracture healing is missing. It is hypothesized that the porous scaffold architecture can serve as a guiding substrate to enable the formation of a structured fibrous network as a prerequirement for later bone formation. An ovine, tibial, 30-mm critical-sized defect is used as a model system to better understand the effect of the scaffold architecture on cell organization, fibrous tissue, and mineralized tissue formation mechanisms in vivo. Tissue regeneration patterns within two geometrically distinct macroscopic regions of a specific scaffold design, the scaffold wall and the endosteal cavity, are compared with tissue formation in an empty defect (negative control) and with cortical bone (positive control). Histology, backscattered electron imaging, scanning small-angle X-ray scattering, and nanoindentation are used to assess the morphology of fibrous and mineralized tissue, to measure the average mineral particle thickness and the degree of alignment, and to map the local elastic indentation modulus. The scaffold proves to function as a guiding substrate to the tissue formation process. It enables the arrangement of a structured fibrous tissue across the entire defect, which acts as a secondary supporting network for cells. Mineralization can then initiate along the fibrous network, resulting in bone ingrowth into a critical-sized defect, although not in complete bridging of the defect. The fibrous network morphology, which in turn is guided by the scaffold architecture, influences the microstructure of the newly formed bone. These results allow a deeper understanding of the mode of mineral tissue formation and the way this is influenced by the scaffold architecture. Copyright © 2012 American Society for Bone and Mineral Research.

Hyaluronic acid enhances the mechanical properties of tissue-engineered cartilage constructs

There is a need for materials that are well suited for cartilage tissue engineering. Hydrogels have emerged as promising biomaterials for cartilage repair, since, like cartilage, they have high water content, and they allow cells to be encapsulated within the material in a genuinely three-dimensional microenvironment. In this study, we investigated the mechanical properties of tissue-engineered cartilage constructs using in vitro culture models incorporating human chondrocytes from osteoarthritis patients. We evaluated hydrogels formed from mixtures of photocrosslinkable gelatin-methacrylamide (Gel-MA) and varying concentrations (0–2%) of hyaluronic acid methacrylate (HA-MA). Initially, only small differences in the stiffness of each hydrogel existed. After 4 weeks of culture, and to a greater extent 8 weeks of culture, HA-MA had striking and concentration dependent impact on the changes in mechanical properties. For example, the initial compressive moduli of cell-laden constructs with 0 and 1% HA-MA were 29 and 41 kPa, respectively. After 8 weeks of culture, the moduli of these constructs had increased to 66 and 147 kPa respectively, representing a net improvement of 69 kPa for gels with 1% HA-MA. Similarly the equilibrium modulus, dynamic modulus, failure strength and failure strain were all improved in constructs containing HA-MA. Differences in mechanical properties did not correlate with glycosaminoglycan content, which did not vary greatly between groups, yet there were clear differences in aggrecan intensity and distribution as assessed using immunostaining. Based on the functional development with time in culture using human chondrocytes, mixtures of Gel-MA and HA-MA are promising candidates for cartilage tissue-engineering applications.

In vivo non-invasive high resolution MR-based method for the determination of the elastic modulus of arterial vessels

The mechanical properties of arterial walls have long been recognized to play an essential role in the development and progression of cardiovascular disease (CVD). Early detection of variations in the elastic modulus of arteries would help in monitoring patients at high cardiovascular risk stratifying them according to risk. An in vivo, non-invasive, high resolution MR-phase-contrast based method for the estimation of the time-dependent elastic modulus of healthy arteries was developed, validated in vitro by means of a thin walled silicon rubber tube integrated into an existing MR-compatible flow simulator and used on healthy volunteers. A comparison of the elastic modulus of the silicon tube measured from the MRI-based technique with direct measurements confirmed the method's capability. The repeatability of the method was assessed. Viscoelastic and inertial effects characterizing the dynamic response of arteries in vivo emerged from the comparison of the pressure waveform and the area variation curve over a period. For all the volunteers who took part in the study the elastic modulus was found to be in the range 50-250 kPa, to increase during the rising part of the cycle, and to decrease with decreasing pressure during the downstroke of systole and subsequent diastole.

Determination of elastic modulus of thin coating materials using nanoindentation and finite element analysis

A deconvolution method that combines nanoindentation and finite element analysis was developed to determine elastic modulus of thin coating layer in a coating-substrate bilayer system. In this method, the nanoindentation experiments were conducted to obtain the modulus of both the bilayer system and the substrate. The finite element analysis was then applied to deconvolve the elastic modulus of the coating. The results demonstrated that the elastic modulus obtained using the developed method was in good agreement with that reported in literature.

Comparison of mechanical and ultrasound elastic modulus of ovine tibial cortical bone

Finite element models of bones can be created by deriving geometry from anx-ray CT scan. Material properties such as the elastic modulus can then be applied using either a single or set of homogeneous values, or individual elements can have local values mapped onto them. Values for the elastic modulus can be derived from the CT density values using an elasticityversus density relationship. Many elasticity–density relationships have been reported in the literature for human bone. However, while ovine in vivo models are common in orthopaedic research, no work has been done to date on creating FE models of ovine bones. To create these models and apply relevant material properties, an ovine elasticity-density relationship needs to be determined. Using fresh frozen ovine tibias the apparent density of regions of interest was determined from a clinical CT scan. The bones were the sectioned into cuboid samples of cortical bone from the regions of interest. Ultrasound was used to determine the elastic modulus in each of three directions – longitudinally, radially and tangentially. Samples then underwent traditional compression testing in each direction. The relationships between apparent density and both ultrasound, and compression modulus in each directionwere determined. Ultrasound testing was found to be a highly repeatable non-destructive method of calculating the elastic modulus, particularly suited to samples of this size. The elasticity-density relationships determined in the longitudinal direction were very similar between the compression and ultrasound data over the density range examined.A clear difference was seen in the elastic modulus between the longitudinal and transverse directions of the bone samples, and a transverse elasticity-density relationship is also reported.

In vivo mri-based estimation of time-dependent elastic modulus in healthy arteries

The elastic properties of the arterial wall have been the subject of physiological, clinical and biomedical research for many years. There is convincing evidence that the elastic properties of the large arteries are seriously impaired in the presence of cardiovascular disease (CVD), due to alterations in the intrinsic structural and functional characteristics of vessels [1]. Early detection of changes in the elastic modulus of arteries would provide a powerful tool for both monitoring patients at high cardiovascular risk and testing the effects of pharmaceuticals aimed at stabilizing existing plaques by stiffening them or lowering the lipids.

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.

Interactive axial shortening of columns and walls in high rise buildings

Concrete is commonly used as a primary construction material for tall building construction. Load bearing components such as columns and walls in concrete buildings are subjected to instantaneous and long term axial shortening caused by the time dependent effects of "shrinkage", "creep" and "elastic" deformations. Reinforcing steel content, variable concrete modulus, volume to surface area ratio of the elements and environmental conditions govern axial shortening. The impact of differential axial shortening among columns and core shear walls escalate with increasing building height. Differential axial shortening of gravity loaded elements in geometrically complex and irregular buildings result in permanent distortion and deflection of the structural frame which have a significant impact on building envelopes, building services, secondary systems and the life time serviceability and performance of a building. Existing numerical methods commonly used in design to quantify axial shortening are mainly based on elastic analytical techniques and therefore unable to capture the complexity of non-linear time dependent effect. Ambient measurements of axial shortening using vibrating wire, external mechanical strain, and electronic strain gauges are methods that are available to verify pre-estimated values from the design stage. Installing these gauges permanently embedded in or on the surface of concrete components for continuous measurements during and after construction with adequate protection is uneconomical, inconvenient and unreliable. Therefore such methods are rarely if ever used in actual practice of building construction. This research project has developed a rigorous numerical procedure that encompasses linear and non-linear time dependent phenomena for prediction of axial shortening of reinforced concrete structural components at design stage. This procedure takes into consideration (i) construction sequence, (ii) time varying values of Young's Modulus of reinforced concrete and (iii) creep and shrinkage models that account for variability resulting from environmental effects. The capabilities of the procedure are illustrated through examples. In order to update previous predictions of axial shortening during the construction and service stages of the building, this research has also developed a vibration based procedure using ambient measurements. This procedure takes into consideration the changes in vibration characteristic of structure during and after construction. The application of this procedure is illustrated through numerical examples which also highlight the features. The vibration based procedure can also be used as a tool to assess structural health/performance of key structural components in the building during construction and service life.

Beat phenomena in metal nanowires, and their implications for resonance-based elastic property measurements

The elastic properties of 1D nanostructures such as nanowires are often measured experimentally through actuation of the nanowire at its resonance frequency, and then relating the resonance frequency to the elastic stiffness using elementary beam theory. In the present work, we utilize large scale molecular dynamics simulations to report a novel beat phenomenon in [110]oriented Ag nanowires. The beat phenomenon is found to arise from the asymmetry of the lattice spacing in the orthogonal elementary directions of the [110] nanowire, i.e. the [-110] and [001] directions, which results in two different principal moments of inertia. Because of this, actuations imposed along any other direction are found to decompose into two orthogonal vibrational components based on the actuation angle relative to these two elementary directions, with this phenomenon being generalizable to <110> FCC nanowires of different materials (Cu, Au, Ni, Pd and Pt). The beat phenomenon is explained using a discrete moment of inertia model based on the hard sphere assumption, the model is utilized to show that surface effects enhance the beat phenomenon, while the effect is reduced with increasing nanowires cross-sectional size or aspect ratio. Most importantly, due to the existence of the beat phenomena, we demonstrate that in resonance experiments only a single frequency component is expected to be observed, particularly when the damping ratio is relatively large or very small. Furthermore, for a large range of actuation angles, the lower frequency is more likely to be detected than the higher one, which implies that experimental predictions of Young’s modulus obtained from resonance may in fact be under predictions. The present study therefore has significant implications for experimental interpretations of Young’s modulus as obtained via resonance testing.

Electrospinning and crosslinking of low-molecular-weight poly(trimethylene carbonate-co-L-lactide) as an elastomeric scaffold for vascular engineering

The growth of suitable tissue to replace natural blood vessels requires a degradable scaffold material that is processable into porous structures with appropriate mechanical and cell growth properties. This study investigates the fabrication of degradable, crosslinkable prepolymers of l-lactide-co-trimethylene carbonate into porous scaffolds by electrospinning. After crosslinking by γ-radiation, dimensionally stable scaffolds were obtained with up to 56% trimethylene carbonate incorporation. The fibrous mats showed Young’s moduli closely matching human arteries (0.4–0.8 MPa). Repeated cyclic extension yielded negligible change in mechanical properties, demonstrating the potential for use under dynamic physiological conditions. The scaffolds remained elastic and resilient at 30% strain after 84 days of degradation in phosphate buffer, while the modulus and ultimate stress and strain progressively decreased. The electrospun mats are mechanically superior to solid films of the same materials. In vitro, human mesenchymal stem cells adhered to and readily proliferated on the three-dimensional fiber network, demonstrating that these polymers may find use in growing artificial blood vessels in vivo.

The modulus of elasticity of steel - is it 200 GPa?

In cold-formed steel construction, the use of a range of thin, high strength steels (0.35 mm thickness and 550 MPa yield stress) has increased significantly in recent times. A good knowledge of the basic mechanical properties of these steels is needed for a satisfactory use of them. In relation to the modulus of elasticity, the current practice is to assume it to be about 200 GPa for all steel grades. However, tensile tests of these steels have consistently shown that the modulus of elasticity varies with grade of steel and thickness. It was found that it increases to values as high as 240 GPa for smaller thicknesses and higher grades of steel. This paper discusses this topic, presents the tensile test results for a number of steel grades and thicknesses, and attempts to develop a relationship between modulus of elasticity, yield stress and thickness for the steel grades considered in this investigation.

Non-invasive clinical measurement of the viscoelastic properties of tendon using acoustic wave transmission

Background. In isotropic materials, the speed of acoustic wave propagation is governed by the bulk modulus and density. For tendon, which is a structural composite of fluid and collagen, however, there is some anisotropy requiring an adjustment for Poisson's ratio. This paper explores these relationships using data collected, in vivo, on human Achilles tendon and then compares estimates of elastic modulus and hysteresis against published values from in vitro mechanical tests. Methods. Measurements using conventional B-model ultrasound imaging, inverse dynamics and acoustic transmission techniques were used to determine dimensions, loading conditions and longitudinal speed of sound in the Achilles tendon during a series of isometric plantar flexion exercises against body weight. Upper and lower bounds for speed of sound versus tensile stress in the tendon were then modelled and estimates of the elastic modulus and hysteresis of the Achilles tendon derived. Results. Axial speed of sound varied between 1850 and 2090 ms-1 with a non-linear, asymptotic dependency on the level of tensile stress (5-35 MPa) in the tendon. Estimates derived for the elastic modulus of the Achilles tendon ranged between 1-2 GPa. Hysteresis derived from models of the stress-strain relationship, ranged from 3-11%. Discussion. Estimates of elastic modulus agree closely with those previously reported from direct measurements obtained via mechanical tensile tests on major weight bearing tendons in vitro [1,2]. Hysteresis derived from models of the stress-strain relationship is consistent with direct measures from various mamalian tendon (7-10%) but is lower than previous estimates in human tendon (17-26%) [3]. This non-invasive method would appear suitable for monitoring changes in tendon properties during dynamic sporting activities.

Effect of vertical misalignment of adjacent sleepers on the increase in dynamic loads around rail joints

In spite of the extensive usage of continuous welded rails, a number of rail joints still exist in the track. Although a number of them exist as part of turnouts in the yards where the speed is not of concern, the Insultated Rail Joints (IRJs) that exist in ballasted tracks remain a source of significant impact loading. A portion of the dynamic load generated at the rail joints due to wheel passage is transmitted to the support system which leads to permanent settlements of the ballast layer with subsequent vertical misalignment of the sleepers around the rail joints. The vertical misalignment of the adjacent sleepers forms a source of high frequency dynamic load raisers causing significant maintenance work including localised grinding of railhead around the joint, re-alignment of the sleepers and/or ballast tamping or track component renewals/repairs. These localised maintenance activities often require manual inspections and disruptions to the train traffic loading to significant costs to the rail industry. Whilst a number of studies have modelled the effect of joints as dips, none have specifically attended to the effect of vertical misalignment of the sleepers on the dynamic response of rail joints. This paper presents a coupled finite element track model and rigid body track-vehicle interaction model through which the effects of vertical of sleepers on the increase in dynamic loads around the IRJ are studied. The finite element track model is employed to determine the generated dip from elastic deformations as well as the vertical displacement of sleepers around the joint. These data (dip and vertical misalignments) are then imported into the rigid body vehicle-track interaction model to calculate the dynamic loads.

Experimental measurement of the mechanical properties of carotid atherothrombotic plaque fibrous cap

Eleven carotid atherothrombotic plaque samples were harvested from patients. Three samples that were highly calcified were discarded, while eight yielded results. The elastic properties of the material were estimated by fitting the measured indentation response to finite element simulations. The methodology was refined and its accuracy quantified using a synthetic rubber. The neo-Hookean form of the material model gave a good fit to the measured response of the tissue. The inferred shear modulus μ was found to be in the range 7-100 kPa, with a median value of 11 kPa. A review of published materials data showed a wide range of material properties for human atherothrombotic tissue. The effects of anisotropy and time dependency in these published results were highlighted. The present measurements were comparable to the static radial compression tests of Lee et al, 1991 [Structure-dependent dynamic behaviour of fibrous caps from human atherosclerotic plaques. Circulation 83, 1764-1770].