Geometric sensitivity of patient-specific finite element models of the spine to variability in user-selected anatomical landmarks

Software to create individualised finite element (FE) models of the osseoligamentous spine using pre-operative computed tomography (CT) data-sets for spinal surgery patients has recently been developed. This study presents a geometric sensitivity analysis of this software to assess the effect of intra-observer variability in user-selected anatomical landmarks. User-selected landmarks on the osseous anatomy were defined from CT data-sets for three scoliosis patients and these landmarks were used to reconstruct patient-specific anatomy of the spine and ribcage using parametric descriptions. The intra-observer errors in landmark co-ordinates for these anatomical landmarks were calculated. FE models of the spine and ribcage were created using the reconstructed anatomy for each patient and these models were analysed for a loadcase simulating clinical flexibility assessment. The intra-observer error in the anatomical measurements was low in comparison to the initial dimensions, with the exception of the angular measurements for disc wedge and zygapophyseal joint (z-joint) orientation and disc height. This variability suggested that CT resolution may influence such angular measurements, particularly for small anatomical features, such as the z-joints, and may also affect disc height. The results of the FE analysis showed low variation in the model predictions for spinal curvature with the mean intra-observer variability substantially less than the accepted error in clinical measurement. These findings demonstrate that intra-observer variability in landmark point selection has minimal effect on the subsequent FE predictions for a clinical loadcase.

Analysis of strain-rate dependent mechanical behavior of single chondrocyte : a finite element study

Various studies have been conducted to investigate the effects of impact loading on cartilage damage and chondrocyte death. These have shown that the rate and magnitude of the applied strain significantly influence chondrocyte death, and that cell death occurred mostly in the superficial zone of cartilage suggesting the need to further understand the fundamental mechanisms underlying the chondrocytes death induced at certain levels of strain-rate. To date there is no comprehensive study providing insight on this phenomenon. The aim of this study is to examine the strain-rate dependent behavior of a single chondrocyte using a computational approach based on Finite Element Method (FEM). An FEM model was developed using various mechanical models, which were Standard Neo-Hookean Solid (SnHS), porohyperelastic (PHE) and poroviscohyperelastic (PVHE) to simulate Atomic Force Microscopy (AFM) experiments of chondrocyte. The PVHE showed, it can capture both relaxation and loading rate dependent behaviors of chondrocytes, accurately compared to other models.

Ability of recombinant human bone morphogenetic protein 2 to enhance bone healing in the presence of tobramycin : evaluation in a rat segmental defect model

Objective To determine whether locally applied tobramycin influences the ability of recombinant human bone morphogenetic protein 2 (rhBMP-2) to heal a segmental defect in the rat femur. Methods The influence of tobramycin on the osteogenic differentiation of mesenchymal stem cells was first evaluated in vitro. For the subsequent, in vivo experiments, a 5-mm segmental defect was created in the right femur of each of 25 Sprague-Dawley rats and stabilized with an external fixator and four Kirschner wires. Rats were divided in four groups: empty control, tobramycin (11 mg)/absorbable collagen sponge, rhBMP-2 (11 μg)/absorbable collagen sponge, and rhBMP-2/absorbable collagen sponge with tobramycin. Bone healing was monitored by radiography at 3 and 8 weeks. Animals were euthanized at 8 weeks and the properties of the defect were compared with the intact contralateral femur. Bone formation in the defect region was assessed by dual-energy x-ray absorptiometry, microcomputed tomography, histology, and mechanical testing. Results Tobramycin exerted a dose-dependent inhibition of alkaline phosphatase induction and calcium deposition by mesenchymal stem cells cultured under osteogenic conditions. The inhibition was reversed in the presence of 500 ng/mL of rhBMP-2. Segmental defects in the rat femora failed to heal in the absence of rhBMP-2. Tobramycin exerted no inhibitory effects on the ability of rhBMP-2 to heal these defects and increased the bone area of the defects treated with rhBMP-2. Data obtained from all other parameters of healing, including dual-energy x-ray absorptiometry, microcomputed tomography, histology, and mechanical testing, were unaffected by tobramycin. Conclusions Although our in vitro results suggested that tobramycin inhibits the osteogenic differentiation of mesenchymal stem cells, this could be overcome by rhBMP-2. Tobramycin did not impair the ability of rhBMP-2 to heal critical-sized femoral defects in rats. Indeed, bone area was increased by nearly 20% in the rhBMP-2 group treated with tobramycin. This study shows that locally applied tobramycin can be used in conjunction with rhBMP-2 to enhance bone formation at fracture sites.

Second-order elastic finite element analysis of structures using a single element per member

Finite element frame analysis programs targeted for design office application necessitate algorithms which can deliver reliable numerical convergence in a practical timeframe with comparable degrees of accuracy, and a highly desirable attribute is the use of a single element per member to reduce computational storage, as well as data preparation and the interpretation of the results. To this end, a higher-order finite element method including geometric non-linearity is addressed in the paper for the analysis of elastic frames for which a single element is used to model each member. The geometric non-linearity in the structure is handled using an updated Lagrangian formulation, which takes the effects of the large translations and rotations that occur at the joints into consideration by accumulating their nodal coordinates. Rigid body movements are eliminated from the local member load-displacement relationship for which the total secant stiffness is formulated for evaluating the large member deformations of an element. The influences of the axial force on the member stiffness and the changes in the member chord length are taken into account using a modified bowing function which is formulated in the total secant stiffness relationship, for which the coupling of the axial strain and flexural bowing is included.

Non-linear refined finite element elastic analysis of steel frames with generalized transverse member loads

In the finite element modelling of steel frames, external loads usually act along the members rather than at the nodes only. Conventionally, when a member is subjected to these transverse loads, they are converted to nodal forces which act at the ends of the elements into which the member is discretised by either lumping or consistent nodal load approaches. For a contemporary geometrically non-linear analysis in which the axial force in the member is large, accurate solutions are achieved by discretising the member into many elements, which can produce unfavourable consequences on the efficacy of the method for analysing large steel frames. Herein, a numerical technique to include the transverse loading in the non-linear stiffness formulation for a single element is proposed, and which is able to predict the structural responses of steel frames involving the effects of first-order member loads as well as the second-order coupling effect between the transverse load and the axial force in the member. This allows for a minimal discretisation of a frame for second-order analysis. For those conventional analyses which do include transverse member loading, prescribed stiffness matrices must be used for the plethora of specific loading patterns encountered. This paper shows, however, that the principle of superposition can be applied to the equilibrium condition, so that the form of the stiffness matrix remains unchanged with only the magnitude of the loading being needed to be changed in the stiffness formulation. This novelty allows for a very useful generalised stiffness formulation for a single higher-order element with arbitrary transverse loading patterns to be formulated. The results are verified using analytical stability function studies, as well as with numerical results reported by independent researchers on several simple structural frames.

Phosphorus-based functional groups as hydrogen bonding templates for rotaxane formation

We report on the use of the hydrogen bond acceptor properties of some phosphorus-containing functional groups for the assembly of a series of [2]rotaxanes. Phosphinamides, and the homologous thio– and selenophosphinamides, act as hydrogen bond acceptors that, in conjunction with an appropriately positioned amide group on the thread, direct the assembly of amide-based macrocycles around the axle to form rotaxanes in up to 60% yields. Employing solely phosphorus-based functional groups as the hydrogen bond accepting groups on the thread, a bis(phosphinamide) template and a phosphine oxide-phosphinamide template afforded the corresponding rotaxanes in 18 and 15 % yields, respectively. X-Ray crystallography of the rotaxanes shows the presence of up to four intercomponent hydrogen bonds between the amide groups of the macrocycle and various hydrogen bond accepting groups on the thread, including rare examples of amide-to-phosphonamide, -thiophosphinamide and -selenophosphinamide groups. With a phosphine oxide-phosphinamide thread, the solid state structure of the rotaxane is remarkable, featuring no direct intercomponent hydrogen bonds but rather a hydrogen bond network involving water molecules that bridge the H-bonding groups of the macrocycle and thread through bifurcated hydrogen bonds. The incorporation of phosphorus-based functional groups into rotaxanes may prove useful for the development of molecular shuttles in which the macrocycle can be used to hinder or expose binding ligating sites for metal-based catalysts.

Self assembly of plasmonic core-satellite nano-assemblies mediated by hyperbranched polymer linkers

The morphology of plasmonic nano-assemblies has a direct influence on optical properties, such as localised surface plasmon resonance (LSPR) and surface enhanced Raman scattering (SERS) intensity. Assemblies with core-satellite morphologies are of particular interest, because this morphology has a high density of hot-spots, while constraining the overall size. Herein, a simple method is reported for the self-assembly of gold NPs nano-assemblies with a core-satellite morphology, which was mediated by hyperbranched polymer (HBP) linkers. The HBP linkers have repeat units that do not interact strongly with gold NPs, but have multiple end-groups that specifically interact with the gold NPs and act as anchoring points resulting in nano-assemblies with a large (~48 nm) core surrounded by smaller (~15 nm) satellites. It was possible to control the number of satellites in an assembly which allowed optical parameters such as SPR maxima and the SERS intensity to be tuned. These results were found to be consistent with finite-difference time domain (FDTD) simulations. Furthermore, the multiplexing of the nano-assemblies with a series of Raman tag molecules was demonstrated, without an observable signal arising from the HBP linker after tagging. Such plasmonic nano-assemblies could potentially serve as efficient SERS based diagnostics or biomedical imaging agents in nanomedicine.

The analysis of large scale data taken from the world groundnut (Arachis hypogaea L.) germplasm collection I. Two-way quantitative data

Data associated with germplasm collections are typically large and multivariate with a considerable number of descriptors measured on each of many accessions. Pattern analysis methods of clustering and ordination have been identified as techniques for statistically evaluating the available diversity in germplasm data. While used in many studies, the approaches have not dealt explicitly with the computational consequences of large data sets (i.e. greater than 5000 accessions). To consider the application of these techniques to germplasm evaluation data, 11328 accessions of groundnut (Arachis hypogaea L) from the International Research Institute for the Semi-Arid Tropics, Andhra Pradesh, India were examined. Data for nine quantitative descriptors measured in the rainy and post-rainy growing seasons were used. The ordination technique of principal component analysis was used to reduce the dimensionality of the germplasm data. The identification of phenotypically similar groups of accessions within large scale data via the computationally intensive hierarchical clustering techniques was not feasible and non-hierarchical techniques had to be used. Finite mixture models that maximise the likelihood of an accession belonging to a cluster were used to cluster the accessions in this collection. The patterns of response for the different growing seasons were found to be highly correlated. However, in relating the results to passport and other characterisation and evaluation descriptors, the observed patterns did not appear to be related to taxonomy or any other well known characteristics of groundnut.

Experimental and finite element analysis to identify the source of vibration of a coach

This paper presents an investigation on the cause of severe vibration problem of a coach with four-cylinder engine running at an idle state using vibration and impact hammer modal experiments to obtain the main vibration frequency components and the natural characteristics of the coach. The vibration results indicate that the main vibration component comes from the vibration transmitted from the engine to the chassis frame, which is closely related with the engine idle speed. Based on structural simulation analysis of the coach’s chassis frame and comparison with modal testing, the coach severe vibration problem was due to coupling resonance between the engine idle frequency and the fourth natural frequency of the chassis frame. The solution to eliminate the vibration problem is provided by changing the local structure stiffness of the chassis frame. The contribution of this paper lies in providing a solution to solve similar problems.

A finite volume scheme with preconditioned Lanczos method for two-dimensional space-fractional reaction–diffusion equations

Fractional differential equations have been increasingly used as a powerful tool to model the non-locality and spatial heterogeneity inherent in many real-world problems. However, a constant challenge faced by researchers in this area is the high computational expense of obtaining numerical solutions of these fractional models, owing to the non-local nature of fractional derivatives. In this paper, we introduce a finite volume scheme with preconditioned Lanczos method as an attractive and high-efficiency approach for solving two-dimensional space-fractional reaction–diffusion equations. The computational heart of this approach is the efficient computation of a matrix-function-vector product f(A)bf(A)b, where A A is the matrix representation of the Laplacian obtained from the finite volume method and is non-symmetric. A key aspect of our proposed approach is that the popular Lanczos method for symmetric matrices is applied to this non-symmetric problem, after a suitable transformation. Furthermore, the convergence of the Lanczos method is greatly improved by incorporating a preconditioner. Our approach is show-cased by solving the fractional Fisher equation including a validation of the solution and an analysis of the behaviour of the model.

An examination of microfinance self-help groups and the poorest of the poor women in Tamil Nadu, India

This study examined the formation and operation of women's microfinance self-help groups in southern India and investigated whether or not the poorest of the poor women were accepted as members of those groups. The study found that caste was used as a selection criterion. Many eligible women excluded themselves from joining the self-help group due to their own lack of education, age, poor health, poverty and lack of trust in the system. The research revealed that self-help groups enhanced women's income and education, improved village infrastructure, and reduced household conflict. Factors that might prevent inclusion of the poorest of the poor in future microfinance programs were identified.

Finite element analysis of the in-plane shear behaviour of masonry panels confined with reinforced grouted cores

A combined experimental and numerical program was conducted to study the in-plane shear behaviour of hollow concrete masonry panels containing reinforced grout cores. This paper is focused on the numerical program. A two dimensional macromodelling strategy was used to simulate the behaviour of the confined masonry (CM) shear panels. Both the unreinforced masonry and the confining element were modelled using macromasonry properties and the steel reinforcement was modelled as an embedded truss element located within the grout using perfectly bonded constraint. The FE model reproduced key behaviours observed in the experiments, including the shear strength, the deformation and the crack patterns of the unconfined and confined masonry panels. The predictions of the validated model were used to evaluate the existing in-plane shear expressions available in the national masonry standards and research publications.

Swelling induced finite strain flexure in a rectangular block of an isotropic elastic material

The deformation of a rectangular block into an annular wedge is studied with respect to the state of swelling interior to the block. Nonuniform swelling fields are shown to generate these flexure deformations in the absence of resultant forces and bending moments. Analytical expressions for the deformation fields demonstrate these effects for both incompressible and compressible generalizations of conventional hyperelastic materials. Existing results in the absence of a swelling agent are recovered as special cases.

Using patient specific finite element models to understand the biomechanics of paediatric spinal deformity surgery

Adolescent idiopathic scoliosis (AIS) is a spinal deformity, which may require surgical correction by attaching rods to the patient’s spine using screws inserted into the vertebrae. Complication rates for deformity correction surgery are unacceptably high. Determining an achievable correction without overloading the adjacent spinal tissues or implants requires an understanding of the mechanical interaction between these components. We have developed novel patient specific modelling software to create individualized finite element models (FEM) representing the thoracolumbar spine and ribcage of scoliosis patients. We are using these models to better understand the biomechanics of spinal deformity correction.

Keyword field-free conjunctive keyword searches on encrypted data and extension for dynamic groups

We consider the following problem: a user stores encrypted documents on an untrusted server, and wishes to retrieve all documents containing some keywords without any loss of data confidentiality. Conjunctive keyword searches on encrypted data have been studied by numerous researchers over the past few years, and all existing schemes use keyword fields as compulsory information. This however is impractical for many applications. In this paper, we propose a scheme of keyword field-free conjunctive keyword searches on encrypted data, which affirmatively answers an open problem asked by Golle et al. at ACNS 2004. Furthermore, the proposed scheme is extended to the dynamic group setting. Security analysis of our constructions is given in the paper.