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.

Seasonal amplitude of hemorrhagic fever with renal syndrome in China : a call for attention to neglected regions

Hemorrhagic fever with renal syndrome (HFRS), a rodent-borne viral disease characterized by fever, hemorrhagic, kidney damage and hypotension, is caused by different species of hantaviruses [1]. Every year, HFRS affects thousands of people in Asia, and more than 90% of these cases are reported in China [2, 3]. Due to its high fatality, HFRS has attracted considerable research attention, and prior studies have predominantly focused on quantifying HFRS morbidity [4], identifying high risk areas [5] and populations [6], or exploring peak time of HFRS occurrence [3]. To date, no study has assessed the seasonal amplitude of HFRS in China, even though it reveals the seasonal fluctuation and thus may provide pivotal information on the possibility of HFRS outbreaks.

Endurance test-rig for diagnostics and prognostics of rolling element bearings

Detection of faults in roller element bearing is a topic widely discussed in the scientific field. Bearings diagnostics is usually performed by analyzing experimental signals, almost always vibration signals, measured during operation. A number of signal processing techniques have been proposed and applied to measured vibrations. The diagnostic effectiveness of the techniques depends on their capacities and on the environmental conditions (i.e. environmental noise). The current trend, especially from an industrial point of view, is to couple the prognostics to the diagnostics. The realization of a prognostic procedure require the definition of parameters able to describe the bearing condition during its operation. Monitoring the values of these parameters during time allows to define their trends depending on the progress of the wear. In this way, a relation between the variation of the selected parameters and the wear progress, useful for diagnostics and prognostics of bearings in real industrial applications, can be established. In this paper, a laboratory test-rig designed to perform endurance tests on roller element bearing is presented. Since the test-rig has operated for a short time, only some preliminary available results are discussed.

Tuning second order cyclostationarity parameters for the diagnostics of rolling element bearings

Failures on rolling element bearings usually originate from cracks that are detectable even in their early stage of propogation by properly analyzing vibration signals measured in the proximity of the bearing. Due to micro-slipping in the roller-races contact, damage-induced vibration signals belong to the family of quasi-periodic signals with a strong second order cyclostationary component. Cyclic coherence and its integrated form are widely considered as the most suitable tools for bearing fault diagnostics and their theoretical bases have been already consolidated. This paper presents how to correctly set the parameters of the cyclostationary analysis tool to be implemented in an automatable algorithm. In the first part of the paper some general guidelines are provided for the specific application. These considerations are further verified, applying cyclostationary tools to data collected in an experimental campaign on a specific test-rig.

Diagnostics of rolling element bearings for the traction system of high speed trains : experimental evidences

Monitoring of the integrity of rolling element bearings in the traction system of high speed trains is a fundamental operation in order to avoid catastrophic failures and to implement effective condition-based maintenance strategies. Diagnostics of rolling element bearings is usually based on vibration signal analysis by means of suitable signal processing techniques. The experimental validation of such techniques has been traditionally performed by means of laboratory tests on artificially damaged bearings, while their actual effectiveness in industrial applications, particularly in the field of rail transport, remains scarcely investigated. This paper will address the diagnostics of bearings taken from the service after a long term operation on a high speed train. These worn bearings have been installed on a test-rig, consisting of a complete full-scale traction system of a high speed train, able to reproduce the effects of wheel-track interaction and bogie-wheelset dynamics. The results of the experimental campaign show that suitable signal processing techniques are able to diagnose bearing failures even in this harsh and noisy application. Moreover, the most suitable location of the sensors on the traction system is also proposed.

Determination of strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes using atomic force microscopy and inverse finite element analysis

The aim of this paper is to determine the strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes, in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of these single cells at four different strain-rates. These results were then employed in inverse finite element analysis (FEA) using Modified Standard neo-Hookean Solid (MSnHS) idealization of these cells to determine their mechanical properties. In addition, a FEA model with a newly developed spring element was employed to accurately simulate AFM evaluation in this study. We report that both cytoskeleton (CSK) and intracellular fluid govern the strain-rate-dependent mechanical property of living cells whereas intracellular fluid plays a predominant role on fixed cells’ behavior. In addition, through the comparisons, it can be concluded that osteocytes are stiffer than chondrocytes at all strain-rates tested indicating that the cells could be the biomarker of their tissue origin. Finally, we report that MSnHS is able to capture the strain-rate-dependent mechanical behavior of osteocyte and chondrocyte for both living and fixed cells. Therefore, we concluded that the MSnHS is a good model for exploration of mechanical deformation responses of single osteocytes and chondrocytes. This study could open a new avenue for analysis of mechanical behavior of osteocytes and chondrocytes as well as other similar types of cells.

Towards a multi-element simulation of heavily over-consolidated bagasse using a material subroutine

A better understanding of the behaviour of prepared cane and bagasse, and the ability to model the mechanical behaviour of bagasse as it is squeezed in a milling unit to extract juice, would help identify how to improve the current process. There are opportunities to decrease bagasse moisture from a milling unit. The behaviour of bagasse in chutes is poorly understood. Previous investigations have shown that juice flow through bagasse obeys Darcy’s permeability law, that the grip of the rough surface of the grooves on the bagasse can be represented by the Mohr-Coulomb failure criterion for soils, and that the internal mechanical behaviour of the bagasse is critical state behaviour similar to that for sand and clay. Progress has been made in the last 11 years towards implementing a mechanical model for bagasse in finite element software. The objective is to be able to correctly simulate various simple mechanical loading conditions measured in the laboratory. Combining these behaviours together is thought to have a high probability of reproducing the complicated stress conditions in a milling unit. This paper reports on progress made towards modelling the fifth and final (and most challenging) of the simple loading conditions: the shearing of heavily over-consolidated bagasse, using a specific model for bagasse in a multi-element simulation.

Computational strategies for surface fitting using thin plate spline finite element methods

Thin plate spline finite element methods are used to fit a surface to an irregularly scattered dataset [S. Roberts, M. Hegland, and I. Altas. Approximation of a Thin Plate Spline Smoother using Continuous Piecewise Polynomial Functions. SIAM, 1:208--234, 2003]. The computational bottleneck for this algorithm is the solution of large, ill-conditioned systems of linear equations at each step of a generalised cross validation algorithm. Preconditioning techniques are investigated to accelerate the convergence of the solution of these systems using Krylov subspace methods. The preconditioners under consideration are block diagonal, block triangular and constraint preconditioners [M. Benzi, G. H. Golub, and J. Liesen. Numerical solution of saddle point problems. Acta Numer., 14:1--137, 2005]. The effectiveness of each of these preconditioners is examined on a sample dataset taken from a known surface. From our numerical investigation, constraint preconditioners appear to provide improved convergence for this surface fitting problem compared to block preconditioners.

Scaling laws for localised states in a nonlocal amplitude equation

It is well known that, although a uniform magnetic field inhibits the onset of small amplitude thermal convection in a layer of fluid heated from below, isolated convection cells may persist if the fluid motion within them is sufficiently vigorous to expel magnetic flux. Such fully nonlinear(‘‘convecton’’) solutions for magnetoconvection have been investigated by several authors. Here we explore a model amplitude equation describing this separation of a fluid layer into a vigorously convecting part and a magnetically-dominated part at rest. Our analysis elucidates the origin of the scaling laws observed numerically to form the boundaries in parameter space of the region of existence of these localised states, and importantly, for the lowest thermal forcing required to sustain them.

Stabilization of experimental bioretention basins during intermittent wetting and drying

Stormwater bioretention basins are subjected to spontaneous intermittent wetting and drying, unlike water treatment filter systems that are subjected to continuous feed. Drinking water filters when constructed new or after back-wash, are subjected to a phase of stabilization. Experiments show that bioretention basins are similarly impacted by intermittent wetting and drying. The common parameter monitored in the stabilisation of filters is the concentration of total solids in the outflow. Filter media in bioretention basins however, consists of a mix of particulate organic matter and fine sand. Organic carbon and solids are therefore needed to be monitored. Four Perspex bioretention filter columns of 94 mm (ID) were packed with a filter layer (800 mm), transition layer and a gravel layer and operated with synthetic stormwater in the laboratory. The filter layer contained 8% organic material by weight. A free board of 350 mm provided detention storage and head to facilitate infiltration. Synthetic stormwater was prepared by adding NH4NO3 (ammonium nitrate) and C2H5NO2 (glycine) and a mixture of kaolinite and montmorillonite clay, to tapwater. The columns were fed with synthetic stormwater with different Antecedent Dry Days (ADD) (0 – 25 day) and constant inflow concentration (2 ppm: nitrate-nitrogen, 1.5 ppm: ammonium-nitrogen, 2.5 ppm: organic-nitrogen 100 ppm: total suspended solids and 7 ppm: organic carbon) at a feed rate of 100mL.min (85.7cm/h). Samples were collected from the outflow at different time intervals between 2 – 150 min from the start of outflow and were tested for Total Suspended Solids (TSS) and Total Organic Carbon (TOC). Both TSS and TOC concentrations in the outflow were observed to be much higher than the concentration of both the parameters in the inflow during the stabilisation period indicating a phase of wash-off (first flush) which lasted for approximately 30 min for both parameters at the beginning of each storm event. The wash-off of TSS and TOC were found to be highly variable depending on the age of the filter and the number of antecedent dry days. The duration of stabilisation phase in the experiments is significant compared with many of the stormwater events. A computational analysis on total mass of each pollutant further affirmed the significance of the first flush of an event on removal of these pollutants. Therefore, the kinetics of the first flush in the stabilisation phase needs to be considered in the performance analysis of the systems.

Comparison of eight published static finite element models of the intact lumbar spine : predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centres around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.

A non-linear interface element model for thin layer high adhesive mortared masonry

A nonlinear interface element modelling method is formulated for the prediction of deformation and failure of high adhesive thin layer polymer mortared masonry exhibiting failure of units and mortar. Plastic flow vectors are explicitly integrated within the implicit finite element framework instead of relying on predictor–corrector like approaches. The method is calibrated using experimental data from uniaxial compression, shear triplet and flexural beam tests. The model is validated using a thin layer mortared masonry shear wall, whose experimental datasets are reported in the literature and is used to examine the behaviour of thin layer mortared masonry under biaxial loading.