A bridging transition technique for the combination of meshfree method with finite element method in 2D solids and structures

For certain continuum problems, it is desirable and beneficial to combine two different methods together in order to exploit their advantages while evading their disadvantages. In this paper, a bridging transition algorithm is developed for the combination of the meshfree method (MM) with the finite element method (FEM). In this coupled method, the meshfree method is used in the sub-domain where the MM is required to obtain high accuracy, and the finite element method is employed in other sub-domains where FEM is required to improve the computational efficiency. The MM domain and the FEM domain are connected by a transition (bridging) region. A modified variational formulation and the Lagrange multiplier method are used to ensure the compatibility of displacements and their gradients. To improve the computational efficiency and reduce the meshing cost in the transition region, regularly distributed transition particles, which are independent of either the meshfree nodes or the FE nodes, can be inserted into the transition region. The newly developed coupled method is applied to the stress analysis of 2D solids and structures in order to investigate its’ performance and study parameters. Numerical results show that the present coupled method is convergent, accurate and stable. The coupled method has a promising potential for practical applications, because it can take advantages of both the meshfree method and FEM when overcome their shortcomings.

Comparison of 3D finite element analysis derived stiffness and BMD to determine the failure load of the excised proximal femur

Introduction: Bone mineral density (BMD) is currently the preferred surrogate for bone strength in clinical practice. Finite element analysis (FEA) is a computer simulation technique that can predict the deformation of a structure when a load is applied, providing a measure of stiffness (Nmm−1). Finite element analysis of X-ray images (3D-FEXI) is a FEA technique whose analysis is derived froma single 2D radiographic image. Methods: 18 excised human femora had previously been quantitative computed tomography scanned, from which 2D BMD-equivalent radiographic images were derived, and mechanically tested to failure in a stance-loading configuration. A 3D proximal femur shape was generated from each 2D radiographic image and used to construct 3D-FEA models. Results: The coefficient of determination (R2%) to predict failure load was 54.5% for BMD and 80.4% for 3D-FEXI. Conclusions: This ex vivo study demonstrates that 3D-FEXI derived from a conventional 2D radiographic image has the potential to significantly increase the accuracy of failure load assessment of the proximal femur compared with that currently achieved with BMD. This approach may be readily extended to routine clinical BMD images derived by dual energy X-ray absorptiometry. Crown Copyright © 2009 Published by Elsevier Ltd on behalf of IPEM. All rights reserved

Reducing MRI gradient coil vibration with rib stiffeners

This work investigates the effect of rib stiffeners on the free and forced vibration of a gradient coil in a Magnetic Resonance Imaging (MRI) scanner. Several reinforcement schemes are studied in this paper. One scheme utilizes the existing holes in the gradient coil structure (typically reserved for magnetic shims) to produce the reinforcement. Non-ferrous, non-magnetic carbon fibre rib stiffeners are employed to fill these holes in several ways to strengthen a gradient coil. Another scheme replaces the inner half of the gradient coil material with a grid of interconnected axial and circumferential rib stiffeners. It is found that the structural stiffness of the gradient coil increases substantially when the coil is reinforced by carbon fibre rib stiffeners. The reinforcement affects the noise and vibration response of the gradient coil structure in the following ways. It increases the frequency range of forced response of the gradient coil at low frequencies due to the increased resonant frequency of the fundamental mode of the coil. Secondly, it reduces the forced response amplitude of the coil structure (which is governed by the structural stiffness of the coil). Thirdly, it reduces the number of natural modes in the low and medium frequency range and therefore lessens the chance of the coil structure being excited resonantly by magnetic resonance signal acquisition sequences. It is shown that gradient coils modelled by solid finite element models have higher stiffness along the coil’s circumference and lower stiffness in the axial direction than those using shell finite element models.

Comparison of 3D finite element analysis derived stiffness and BMD to determine the failure load of the excised proximal femur

Bone mineral density (BMD) is currently the preferred surrogate for bone strength in clinical practice. Finite element analysis (FEA) is a computer simulation technique that can predict the deformation of a structure when a load is applied, providing a measure of stiffness (N mm− 1). Finite element analysis of X-ray images (3D-FEXI) is a FEA technique whose analysis is derived from a single 2D radiographic image. This ex-vivo study demonstrates that 3D-FEXI derived from a conventional 2D radiographic image has the potential to significantly increase the accuracy of failure load assessment of the proximal femur compared with that currently achieved with BMD.

Stability and convergence of a new explicit finite-difference approximation for the variable-order nonlinear fractional diffusion equation

In this paper, we consider the variable-order nonlinear fractional diffusion equation View the MathML source where xRα(x,t) is a generalized Riesz fractional derivative of variable order View the MathML source and the nonlinear reaction term f(u,x,t) satisfies the Lipschitz condition |f(u1,x,t)-f(u2,x,t)|less-than-or-equals, slantL|u1-u2|. A new explicit finite-difference approximation is introduced. The convergence and stability of this approximation are proved. Finally, some numerical examples are provided to show that this method is computationally efficient. The proposed method and techniques are applicable to other variable-order nonlinear fractional differential equations.

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.

An experimental and finite element investigation of the biomechanics of vertebral compression fractures

Osteoporosis is a disease characterized by low bone mass and micro-architectural deterioration of bone tissue, with a consequent increase in bone fragility and susceptibility to fracture. Osteoporosis affects over 200 million people worldwide, with an estimated 1.5 million fractures annually in the United States alone, and with attendant costs exceeding $10 billion dollars per annum. Osteoporosis reduces bone density through a series of structural changes to the honeycomb-like trabecular bone structure (micro-structure). The reduced bone density, coupled with the microstructural changes, results in significant loss of bone strength and increased fracture risk. Vertebral compression fractures are the most common type of osteoporotic fracture and are associated with pain, increased thoracic curvature, reduced mobility, and difficulty with self care. Surgical interventions, such as kyphoplasty or vertebroplasty, are used to treat osteoporotic vertebral fractures by restoring vertebral stability and alleviating pain. These minimally invasive procedures involve injecting bone cement into the fractured vertebrae. The techniques are still relatively new and while initial results are promising, with the procedures relieving pain in 70-95% of cases, medium-term investigations are now indicating an increased risk of adjacent level fracture following the procedure. With the aging population, understanding and treatment of osteoporosis is an increasingly important public health issue in developed Western countries. The aim of this study was to investigate the biomechanics of spinal osteoporosis and osteoporotic vertebral compression fractures by developing multi-scale computational, Finite Element (FE) models of both healthy and osteoporotic vertebral bodies. The multi-scale approach included the overall vertebral body anatomy, as well as a detailed representation of the internal trabecular microstructure. This novel, multi-scale approach overcame limitations of previous investigations by allowing simultaneous investigation of the mechanics of the trabecular micro-structure as well as overall vertebral body mechanics. The models were used to simulate the progression of osteoporosis, the effect of different loading conditions on vertebral strength and stiffness, and the effects of vertebroplasty on vertebral and trabecular mechanics. The model development process began with the development of an individual trabecular strut model using 3D beam elements, which was used as the building block for lattice-type, structural trabecular bone models, which were in turn incorporated into the vertebral body models. At each stage of model development, model predictions were compared to analytical solutions and in-vitro data from existing literature. The incremental process provided confidence in the predictions of each model before incorporation into the overall vertebral body model. The trabecular bone model, vertebral body model and vertebroplasty models were validated against in-vitro data from a series of compression tests performed using human cadaveric vertebral bodies. Firstly, trabecular bone samples were acquired and morphological parameters for each sample were measured using high resolution micro-computed tomography (CT). Apparent mechanical properties for each sample were then determined using uni-axial compression tests. Bone tissue properties were inversely determined using voxel-based FE models based on the micro-CT data. Specimen specific trabecular bone models were developed and the predicted apparent stiffness and strength were compared to the experimentally measured apparent stiffness and strength of the corresponding specimen. Following the trabecular specimen tests, a series of 12 whole cadaveric vertebrae were then divided into treated and non-treated groups and vertebroplasty performed on the specimens of the treated group. The vertebrae in both groups underwent clinical-CT scanning and destructive uniaxial compression testing. Specimen specific FE vertebral body models were developed and the predicted mechanical response compared to the experimentally measured responses. The validation process demonstrated that the multi-scale FE models comprising a lattice network of beam elements were able to accurately capture the failure mechanics of trabecular bone; and a trabecular core represented with beam elements enclosed in a layer of shell elements to represent the cortical shell was able to adequately represent the failure mechanics of intact vertebral bodies with varying degrees of osteoporosis. Following model development and validation, the models were used to investigate the effects of progressive osteoporosis on vertebral body mechanics and trabecular bone mechanics. These simulations showed that overall failure of the osteoporotic vertebral body is initiated by failure of the trabecular core, and the failure mechanism of the trabeculae varies with the progression of osteoporosis; from tissue yield in healthy trabecular bone, to failure due to instability (buckling) in osteoporotic bone with its thinner trabecular struts. The mechanical response of the vertebral body under load is highly dependent on the ability of the endplates to deform to transmit the load to the underlying trabecular bone. The ability of the endplate to evenly transfer the load through the core diminishes with osteoporosis. Investigation into the effect of different loading conditions on the vertebral body found that, because the trabecular bone structural changes which occur in osteoporosis result in a structure that is highly aligned with the loading direction, the vertebral body is consequently less able to withstand non-uniform loading states such as occurs in forward flexion. Changes in vertebral body loading due to disc degeneration were simulated, but proved to have little effect on osteoporotic vertebra mechanics. Conversely, differences in vertebral body loading between simulated invivo (uniform endplate pressure) and in-vitro conditions (where the vertebral endplates are rigidly cemented) had a dramatic effect on the predicted vertebral mechanics. This investigation suggested that in-vitro loading using bone cement potting of both endplates has major limitations in its ability to represent vertebral body mechanics in-vivo. And lastly, FE investigation into the biomechanical effect of vertebroplasty was performed. The results of this investigation demonstrated that the effect of vertebroplasty on overall vertebra mechanics is strongly governed by the cement distribution achieved within the trabecular core. In agreement with a recent study, the models predicted that vertebroplasty cement distributions which do not form one continuous mass which contacts both endplates have little effect on vertebral body stiffness or strength. In summary, this work presents the development of a novel, multi-scale Finite Element model of the osteoporotic vertebral body, which provides a powerful new tool for investigating the mechanics of osteoporotic vertebral compression fractures at the trabecular bone micro-structural level, and at the vertebral body level.

Doped poly(vinylidene fluoride) as solid polymer electrolyte in nanocrystalline dye-sensitized solar cell

A new solid composite polymer electrolyte was reported by incorporating Azino-bis-(3-ethyl benzo thiazoline-6-sulphonate) ion [ABTS] as dopant in poly(vinylidene flouride) along with redox couple (1-/13-). Under certain conditions, the electrolyte composition forms brush like nano-rods while it is doped with Azino-bis-(3-ethly) benzo thiazoline-6-sulphonate) ion [ABTS], a pi-electron donor. The polymer electrolyte forms nanoscale interpenetrating network with the crystalline order of the polymer electrolyte that seems to be a desirable architecture for the active layer of the photoelectrochemical cell. With this new polymer electrolyte, dye-sensitized solar cell was fabricated using N3 dye absorbed over Ti02- nonoparticles (photoanode) and conducting carbon cement coated on the conducting press (FTO, photocathode). This polymer composite has been successfully used as a promising candidate as solid polymer electrolyte in nanocrystalline dye-sensitized solar cell.

Solid-state nanocrystalline dye-sensitized solar cell using 2-mercapto benzimidazole doped poly(ethylene oxide) as a new polymer electrolyte

New composite doped poly (ethylene oxide) polymer electrolyte was developed using 2-mercapto benzimidazole as plasticizer and iodide/triiodide as redox couple. The fabrication of the cell involves Poly(ethylene oxide)/ 2-mercapto benzimidazole / iodide/triiodide as polymer electrolyte in dye-sensitized solar cell fabricated with N3 dye and TiO2 nanoparticles as the photoanode and Platinum coated FTO (fluorine doped SnO2) as counter electrode. The current-volatage characteristics under simulated sunlight AM1.5 shows a short circuit current Isc of 8.7mA and open circuit photovoltage 508 mV. The conductivity measurements for the new polymer electrolyte and the photoelectrochemical measurments were carried out systematically. In 2-mercapto benzimidazole the electron rich sulphur and nitrogen atoms, act as pi-electron donors that form good interaction with iodine which plays a vital role in the performance of the fabricated dye-sensitized solar cells. The resonance effect increases the stability of the cell to a considerable extent. These results suggest that the new composite polymer electrolyte performs as a promising new doped polymer-electrolyte.

Evaluation of piezoelectric PVDF polymers for use in space environments. III. Comparison of the effects of vacuum UV and gamma radiation

Films of piezoelectric PVDF and P(VDF-TrFE) were exposed to vacuum UV (115-300 nm VUV) and -radiation to investigate how these two forms of radiation affect the chemical, morphological, and piezoelectric properties of the polymers. The extent of crosslinking was almost identical in both polymers after -irradiation, but surprisingly, was significantly higher for the TrFE copolymer after VUV-irradiation. Changes in the melting behavior were also more significant in the TrFE copolymer after VUV-irradiation due to both surface and bulk crosslinking, compared with only surface crosslinking for the PVDF films. The piezoelectric properties (measured using d33 piezoelectric coefficients and D-E hysteresis loops) were unchanged in the PVDF homopolymer, while the TrFE copolymer exhibited more narrow D-E loops after exposure to either - or VUV-radiation. The more severe damage to the TrFE copolymer in comparison with the PVDF homopolymer after VUV-irradiation is explained by different energy deposition characteristics. The short wavelength, highly energetic photons are undoubtedly absorbed in the surface layers of both polymers, and we propose that while the longer wavelength components of the VUV-radiation are absorbed by the bulk of the TrFE copolymer causing crosslinking, they are transmitted harmlessly in the PVDF homopolymer.

Evaluation of piezoelectric PVDF polymers for use in space environments. II. Effects of atomic oxygen and vacuum UV exposure

The effects of atomic oxygen (AO) and vacuum UV radiation simulating low Earth orbit conditions on two commercially available piezoelectric polymer films, poly(vinylidene fluoride) (PVDF) and poly(vinylidene fluoride-trifluoroethylene) P(VDF-TrFE), have been studied. Surface erosion and pattern development are significant for both polymers. Erosion yields were determined as 2.8 � 10�24 cm3/atom for PVDF and 2.5 � 10�24 cm3/atom for P(VDF-TrFE). The piezoelectric properties of the residual material of both polymers were largely unchanged after exposure, although a slight shift in the Curie transition of the P(VDF-TrFE) was observed. A lightly cross-linked network was formed in the copolymer presumably because of penetrating vacuum ultraviolet (VUV) radiation, while the homopolymer remained uncross-linked. These differences were attributed to varying degrees of crystallinity and potentially greater absorption, and hence damage, of VUV radiation in P(VDFTrFE) compared with PVDF.

Evaluation of piezoelectric poly(vinylidene fluoride) polymers for use in space environments. I. Temperature limitations

Smart materials, such as thin-film piezoelectric polymers, are interesting for potential applications on Gossamer spacecraft. This investigation aims to predict the performance and long-term stability of the piezoelectric properties of poly(vinylidene fluoride) (PVDF) and its copolymers under conditions simulating the low-Earthorbit environment. To examine the effects of temperature on the piezoelectric properties of PVDF, poly(vinylidenefluoride-co-trifluoroethylene), and poly(vinylidenefluoride-cohexafluoropropylene), the d33 piezoelectric coefficients were measured up to 160 8C, and the electric displacement/electric field (D–E) hysteresis loops were measured from �80 to þ110 8C. The room-temperature d33 coefficient of PVDF homopolymer films, annealed at 50, 80, and 125 8C, dropped rapidly within a few days of thermal exposure and then remained unchanged. In contrast, the TrFE copolymer exhibited greater thermal stability than the homopolymer, with d33 remaining almost unchanged up to 125 8C. The HFP copolymer exhibited poor retention of d33 at temperatures above 80 8C. In situ D–E loop measurements from �80 to þ110 8C showed that the remanent polarization of the TrFE copolymer was more stable than that of the PVDF homopolymer. D–E hysteresis loop and d33 results were also compared with the deflection of the PVDF homopolymer and TrFE copolymer bimorphs tested over a wide temperature range.

Selection and optimization of piezoelectric polyvinylidene fluoride polymers for adaptive optics in space environments

Various piezoelectric polymers based on polyvinylidene fluoride (PVDF) are of interest for large aperture space-based telescopes. Dimensional adjustments of adaptive polymer films depend on charge deposition and require a detailed understanding of the piezoelectric material responses which are expected to deteriorate owing to strong vacuum UV, � -, X-ray, energetic particles and atomic oxygen exposure. We have investigated the degradation of PVDF and its copolymers under various stress environments detrimental to reliable operation in space. Initial radiation aging studies have shown complex material changes with lowered Curie temperatures, complex material changes with lowered melting points, morphological transformations and significant crosslinking, but little influence on piezoelectric d33 constants. Complex aging processes have also been observed in accelerated temperature environments inducing annealing phenomena and cyclic stresses. The results suggest that poling and chain orientation are negatively affected by radiation and temperature exposure. A framework for dealing with these complex material qualification issues and overall system survivability predictions in low earth orbit conditions has been established. It allows for improved material selection, feedback for manufacturing and processing, material optimization/stabilization strategies and provides guidance on any alternative materials.