Comparison of eight published lumbar spine finite element models

Due to its ability to represent intricate systems with material nonlinearities as well as irregular loading, boundary, geometrical and material domains, the finite element (FE) method has been recognized as an important computational tool in spinal biomechanics. Current FE models generally account for a single distinct spinal geometry with one set of material properties despite inherently large inter-subject variability. The uncertainty and high variability in tissue material properties, geometry, loading and boundary conditions has cast doubt on the reliability of their predictions and comparability with reported in vitro and in vivo values. A multicenter study was undertaken to compare the results of eight well-established models of the lumbar spine that have been developed, validated and applied for many years. Models were subjected to pure and combined loading modes and their predictions were 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; 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 previously published median in vitro values. However, the ranges of predictions were larger and exceeded the in vitro ranges, especially for facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with in vivo values. The simulations yielded median facet joint forces of 0 N in flexion, 38 N in extension, 14 N in lateral bending and 60 N in axial rotation that could not be validated due to the paucity of in vivo facet joint forces. In light of high inter-subject variability, one must be cautious when generalizing predictions obtained from one deterministic model. This study demonstrates however that the predictive power increases when FE models are combined together. The median of individual numerical results can hence be used as an improved tool in order to estimate the response of the lumbar spine.

Soft tissue human thigh and buttock finite element model to simulate vehicle seat cushion indentation

Accurate modelling of automotive occupant posture is strongly related to the mechanical interaction between human body soft tissue and flexible seat components. This paper presents a finite-element study simulating the deflection of seat cushion foam and supportive seat structures, as well as human buttock and thigh soft tissue when seated. The thigh-buttock surface shell model was based on 95th percentile male subject scan data and made of two layers, covering thin to moderate thigh and buttock proportions. To replicate the effects of skin and fat, the neoprene rubber layer was modelled as a hyperelastic material with viscoelastic behaviour. The analytical seat model is based on a Ford production seat. The result of the finite-element indentation simulation is compared to a previous simulation of an indentation with a hard shell human model of equal geometry, and to the physical indentation result. We conclude that SAE composite buttock form and human-seat indentation of a suspended seat cushion can be validly simulated.

Effects of strain artefacts arising from a pre-defined callus domain in models of bone healing mechanobiology

Iterative computational models have been used to investigate the regulation of bone fracture healing by local mechanical conditions. Although their predictions replicate some mechanical responses and histological features, they do not typically reproduce the predominantly radial hard callus growth pattern observed in larger mammals. We hypothesised that this discrepancy results from an artefact of the models’ initial geometry. Using axisymmetric finite element models, we demonstrated that pre-defining a field of soft tissue in which callus may develop introduces high deviatoric strains in the periosteal region adjacent to the fracture. These bone-inhibiting strains are not present when the initial soft tissue is confined to a thin periosteal layer. As observed in previous healing models, tissue differentiation algorithms regulated by deviatoric strain predicted hard callus forming remotely and growing towards the fracture. While dilatational strain regulation allowed early bone formation closer to the fracture, hard callus still formed initially over a broad area, rather than expanding over time. Modelling callus growth from a thin periosteal layer successfully predicted the initiation of hard callus growth close to the fracture site. However, these models were still susceptible to elevated deviatoric strains in the soft tissues at the edge of the hard callus. Our study highlights the importance of the initial soft tissue geometry used for finite element models of fracture healing. If this cannot be defined accurately, alternative mechanisms for the prediction of early callus development should be investigated.

Microscale consolidation analysis of relaxation behavior of single living chondrocytes subjected to varying strain-rates

Besides the elastic stiffness, the relaxation behavior of single living cells is also of interest of various researchers when studying cell mechanics. It is hypothesized that the relaxation response of the cells is governed by both intrinsic viscoelasticity of the solid phase and fluid-solid interactions mechanisms. There are a number of mechanical models have been developed to investigate the relaxation behavior of single cells. However, there is lack of model enable to accurately capture both of the mechanisms. Therefore, in this study, the porohyperelastic (PHE) model, which is an extension of the consolidation theory, combined with inverse Finite Element Analysis (FEA) technique was used at the first time to investigate the relaxation response of living chondrocytes. This model was also utilized to study the dependence of relaxation behavior of the cells on strain-rates. The stress-relaxation experiments under the various strain-rates were conducted with the Atomic Force Microscopy (AFM). The results have demonstrated that the PHE model could effectively capture the stress-relaxation behavior of the living chondrocytes, especially at intermediate to high strain-rates. Although this model gave some errors at lower strain-rates, its performance was acceptable. Therefore, the PHE model is properly a promising model for single cell mechanics studies. Moreover, it has been found that the hydraulic permeability of living chondrocytes reduced with decreasing of strain-rates. It might be due to the intracellular fluid volume fraction and the fluid pore pressure gradients of chondrocytes were higher when higher strain-rates applied.

A finite element modelling method for thin layer mortared masonry systems

This paper deals with a finite element modelling method for thin layer mortared masonry systems. In this method, the mortar layers including the interfaces are represented using a zero thickness interface element and the masonry units are modelled using an elasto-plastic, damaging solid element. The interface element is formulated using two regimes; i) shear-tension and ii) shearcompression. In the shear-tension regime, the failure of joint is consiedered through an eliptical failure criteria and in shear-compression it is considered through Mohr Coulomb type failure criterion. An explicit integration scheme is used in an implicit finite element framework for the formulation of the interface element. The model is calibrated with an experimental dataset from thin layer mortared masonry prism subjected to uniaxial compression, a triplet subjected to shear loads a beam subjected to flexural loads and used to predict the response of thin layer mortared masonry wallettes under orthotropic loading. The model is found to simulate the behaviour of a thin layer mortated masonry shear wall tested under pre-compression and inplane shear quite adequately. The model is shown to reproduce the failure of masonry panels under uniform biaxial state of stresses.

Predicting surgical outcomes for paediatric spinal deformity patients using subject specific finite element models

INTRODUCTION 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. Our novel patient specific modelling software creates individualized finite element models (FEM) representing the thoracolumbar spine and ribcage of scoliosis patients. We have recently applied the model to investigate the influence of increasing magnitudes of surgically applied corrective force on predicted deformity correction...

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

A three-dimensional hybrid smoothed finite element method (H-SFEM) for nonlinear solid mechanics problems

This paper presents a novel three-dimensional hybrid smoothed finite element method (H-SFEM) for solid mechanics problems. In 3D H-SFEM, the strain field is assumed to be the weighted average between compatible strains from the finite element method (FEM) and smoothed strains from the node-based smoothed FEM with a parameter α equipped into H-SFEM. By adjusting α, the upper and lower bound solutions in the strain energy norm and eigenfrequencies can always be obtained. The optimized α value in 3D H-SFEM using a tetrahedron mesh possesses a close-to-exact stiffness of the continuous system, and produces ultra-accurate solutions in terms of displacement, strain energy and eigenfrequencies in the linear and nonlinear problems. The novel domain-based selective scheme is proposed leading to a combined selective H-SFEM model that is immune from volumetric locking and hence works well for nearly incompressible materials. The proposed 3D H-SFEM is an innovative and unique numerical method with its distinct features, which has great potential in the successful application for solid mechanics problems.

Mechanical analysis of foot plantar fascia during normal walking condition

Foot plantar fascia is an important foot tissue in stabilizing the longitudinal arch of human foot. Direct measurement to monitor the mechanical situation of plantar fascia at human locomotion is difficult. The purpose of this study was to construct a three-dimensional finite element model of the foot to calculate the internal stress/strain value of plantar fascia during different stage of gait. The simulated stress distribution of plantar fascia was the lowest at heel-strike, which concentrated on the medial side of calcaneal tubercle. The peak stress of plantar fascia was appeared at push-off, and the value is more than 5 times of the heel-strike position. Current FE model was able to explore the plantar fascia tension trend at the main sub-phases of foot. More detailed fascia model and intrinsic muscle forces could be developed in the further study.

Normalized wall index specific and MRI-based stress analysis of atherosclerotic carotid plaques: A study comparing acutely symptomatic and asymptomatic patients

Background: Biomechanical stresses play an important role in determining plaque stability. Quantification of these simulated stresses can be potentially used to assess plaque vulnerability and differentiate different patient groups. Methods and Results: 54 asymptomatic and 45 acutely symptomatic patients underwent in vivo multicontrast magnetic resonance imaging (MRI) of the carotid arteries. Plaque geometry used for finite element analysis was derived from in vivo MRI at the sites of maximum and minimum plaque burden. In total, 198 slices were used for the computational simulations. A pre-shrink technique was used to refine the simulation. Maximum principle stress at the vulnerable plaque sites (ie, critical stress) was extracted for the selected slices and a comparison was performed between the 2 groups. Critical stress in the slice with maximum plaque burden is significantly higher in acutely symptomatic patients as compared to asymptomatic patients (median, inter quartile range: 198.0 kPa (119.8-359.0 kPa) vs 138.4 kPa (83.8-242.6 kPa), P=0.04). No significant difference was found in the slice with minimum plaque burden between the 2 groups (196.7 kPa (133.3-282.7 kPa) vs 182.4 kPa (117.2-310.6 kPa), P=0.82). Conclusions: Acutely symptomatic carotid plaques have significantly high biomechanical stresses than asymptomatic plaques. This might be potentially useful for establishing a biomechanical risk stratification criteria based on plaque burden in future studies.

Stress analysis of carotid plaque based on in vivo MRI of acute symptomatic and asymptomatic patients

Stress analysis within carotid plaques based on in vivo MR imaging has shown to be useful for the identification of vulnerable atheroma. This study is to investigate whether magnetic resonance imaging (MRI) based-biomechanical stress analysis of carotid plaques can differentiate acute symptomatic and asymptomatic patients. 54 asymptomatic and 45 acute symptomatic patients underwent in vivo multi-contrast MRI of the carotid arteries. Plaque geometry used for finite element analysis was derived from in vivo MR images at the site of maximum and minimum plaque burden. In total 198 slices were used for the computational simulations. A pre shrink technique was used to refine the simulation. Maximum principle stress at the vulnerable plaque sites (i.e. critical stress) was extracted for the selected slices and a comparison was performed between the two groups. Critical stress at the site of maximum plaque burden is significantly higher in acute symptomatic patients as compared to asymptomatic patients [median: 198.0kPa (inter quartile range (IQR) = (119.8 - 359.0) vs. 138.4kPa (83.8, 242.6), p=0.04]. No significant difference was found at the minimum plaque burden site between the two groups [196.7kPa (133.3- 282.7) vs. 182.4kPa (117.2 - 310. 6), p=0.82). Stress analysis at the site of maximal plaque burden can be effectively used for differentiating acute symptomatic carotid plaques from asymptomatic plaques. This maybe potentially used for development of biomechanical risk stratification criteria based on plaque burden in future studies.

Structural analysis and magnetic resonance imaging predict plaque vulnerability: A study comparing symptomatic and asymptomatic individuals

Background: More than half of all cerebral ischemic events are the result of rupture of extracranial plaques. The clinical determination of carotid plaque vulnerability is currently based solely on luminal stenosis; however, it has been increasingly suggested that plaque morphology and biomechanical stress should also be considered. We used finite element analysis based on in vivo magnetic resonance imaging (MRI) to simulate the stress distributions within plaques of asymptomatic and symptomatic individuals. Methods: Thirty nonconsecutive subjects (15 symptomatic and 15 asymptomatic) underwent high-resolution multisequence in vivo MRI of the carotid bifurcation. Stress analysis was performed based on the geometry derived from in vivo MRI of the carotid artery at the point of maximal stenosis. The finite element analysis model considered plaque components to be hyperelastic. The peak stresses within the plaques of symptomatic and asymptomatic individuals were compared. Results: High stress concentrations were found at the shoulder regions of symptomatic plaques, and the maximal stresses predicted in this group were significantly higher than those in the asymptomatic group (508.2 ± 193.1 vs 269.6 ± 107.9 kPa; P = .004). Conclusions: Maximal predicted plaque stresses in symptomatic patients were higher than those predicted in asymptomatic patients by finite element analysis, suggesting the possibility that plaques with higher stresses may be more prone to be symptomatic and rupture. If further validated by large-scale longitudinal studies, biomechanical stress analysis based on high resolution in vivo MRI could potentially act as a useful tool for risk assessment of carotid atheroma. It may help in the identification of patients with asymptomatic carotid atheroma at greatest risk of developing symptoms or mild-to-moderate symptomatic stenoses, which currently fall outside current clinical guidelines for intervention.

Analysis Of Curved Laminated Beams Of Bimodulus Composite-Materials

Composite materials exhibiting different moduli in tension and in compression, commonly called as bimodular composites are being used in many engineering fields. A finite element analysis is carried out for small deflection static behavior of laminated curved beams of bi modulus materials for both solid and hollow circular cross-sections using an iterative procedure. The finite element has 16 d.o.f. and uses the displacement field in terms of first order Hermite in terpolation polynomials. The neutral surface, i.e. the locus of points having zero axial strain is found to vary drastically depending on the loading, lay up schemes and radius of curvature. As il lustrations, plots of the cross-sections of the ruled neutral-surface are presented for some of the investigated cases. Using this element a few problems of curved laminated beams of bimodulus materials are solved for both solid and hollow circular cross-sections.

Analysis of flange joints with elastic bolt

The contact zone and pressure distribution between two elastic plates joined by an elastic bolt and nut are estimated using finite element analysis. Smooth interfacial conditions are assumed in all the regions of contact. Eight node axisymmetric ring elements are used to model the structure. The matrix solution is obtained through frontal technique and this solution technique is shown to be very efficient for the iterative scheme adopted to determine the extent of contact. A parametric study is conducted varying the elastic properties of bolt and plate materials, bolt head diameter and thickness of the plates. The method of approach presented in this paper provides a solution with a realistic idealization of tension flange joints.