Assessment of the Effect of Mesh Density on the Material Property Discretisation Within QCT Based FE Models: A Practical Example Using the Implanted Proximal Tibia

A three-dimensional, quantitative computed tomography based finite element model of a proximal implanted tibia was analysed in order to assess the effect of mesh density on material property discretisation and the resulting influence on the predicted stress distribution. The mesh was refined on the contact surfaces (matched meshes) with element sizes of 3, 2, 1.4, 1 and 0.8 mm. The same loading conditions were used in all models (bi-condylar load: 60% medial, 40% lateral). Significant variations were observed in the modulus distributions between the coarsest and finest mesh densities. Poor discretisation of the material properties also resulted in poor correlations of the stresses and risk ratios between the coarsest and finest meshes. Little difference in Young's modulus, von Mises stress and risk ratio distributions were observed between the three finest models; hence, it was concluded that for this particular case an element size of 1.4 mm on the contact surfaces was enough to properly describe the stiffness, stress and risk ratio distributions within the bone. Poor convergence of the material property distribution occurred when the element size was significantly larger than the pixel size of the source CT data. It was concluded that unless there is convergence in the Young's modulus distribution, convergence of the stress field or of other parameters of interest will not occur either.     

Investigation of fixation screw pull-out strength on human spine

For the purpose of fixation and stabilizing the spine, the pull-out strength of the screw is one of the most important factors to be considered. The material properties of the bone, coupled with the principal dimensions of the screw such as major diameter, minor diameter, pitch and purchase length, may affect the pull-out strength of the fixation screw. In this study, the effects of various factors on the bone screw pull-out strength were studied using finite element method of analysis. A three-dimensional finite element model simulating the threaded connection of bone and surgical screw was constructed. The behavior of both the bone and the screw during screw pull-out were analyzed and discussed. The effect of the screw parameters on the screw pull-out strength was also investigated. The results showed that failure of the connection was due to bone shearing which occurred along a cylindrical surface determined by the outer perimeter of the screw and penetrated length. The distribution of the shear stress along the threaded length in the bone was nearly uniform. The effect of the major diameter on the screw pull-out strength was more significant than those of the minor diameter and the pitch. The minor diameter hardly affected the pull-out strength.     

Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1-L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.     

Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1–L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.     

An approximate model for cancellous bone screw fixation

This paper presents a finite element (FE) model to identify parameters that affect the performance of an improved cancellous bone screw fixation technique, and hence potentially improve fracture treatment. In cancellous bone of low apparent density, it can be difficult to achieve adequate screw fixation and hence provide stable fracture fixation that enables bone healing. Data from predictive FE models indicate that cements can have a significant potential to improve screw holding power in cancellous bone. These FE models are used to demonstrate the key parameters that determine pull-out strength in a variety of screw, bone and cement set-ups, and to compare the effectiveness of different configurations. The paper concludes that significant advantages, up to an order of magnitude, in screw pull-out strength in cancellous bone might be gained by the appropriate use of a currently approved calcium phosphate cement.     

Assessment of the effect of mesh density on the material property discretisation within QCT based FE models: a practical example using the implanted proximal tibia

A three-dimensional, quantitative computed tomography based finite element model of a proximal implanted tibia was analysed in order to assess the effect of mesh density on material property discretisation and the resulting influence on the predicted stress distribution. The mesh was refined on the contact surfaces (matched meshes) with element sizes of 3, 2, 1.4, 1 and 0.8 mm. The same loading conditions were used in all models (bi-condylar load: 60% medial, 40% lateral). Significant variations were observed in the modulus distributions between the coarsest and finest mesh densities. Poor discretisation of the material properties also resulted in poor correlations of the stresses and risk ratios between the coarsest and finest meshes. Little difference in Young's modulus, von Mises stress and risk ratio distributions were observed between the three finest models; hence, it was concluded that for this particular case an element size of 1.4 mm on the contact surfaces was enough to properly describe the stiffness, stress and risk ratio distributions within the bone. Poor convergence of the material property distribution occurred when the element size was significantly larger than the pixel size of the source CT data. It was concluded that unless there is convergence in the Young's modulus distribution, convergence of the stress field or of other parameters of interest will not occur either.     

Micro-CT and micro-FE analysis of pedicle screw fixation under different loading conditions

Anchorage of pedicle screw instrumentation in the elderly spine with poor bone quality remains challenging. In this study, micro finite element (µFE) models were used to assess the specific influence of screw design and the relative contribution of local bone density to fixation mechanics. These were created from micro computer tomography (µCT) scans of vertebras implanted with two types of pedicle screws, including a full region-or-interest of 10 mm radius around each screw, as well as submodels for the pedicle and inner trabecular bone of the vertebral body. The local bone volume fraction (BV/TV) calculated from the µCT scans around different regions of the screw (pedicle, inner trabecular region of the vertebral body) were then related to the predicted stiffness in simulated pull-out tests as well as to the experimental pull-out and torsional fixation properties mechanically measured on the corresponding specimens. Results show that predicted stiffness correlated excellently with experimental pull-out strength (R² > 0.92, p < .043), better than regional BV/TV alone (R2 = 0.79, p = .003). They also show that correlations between fixation properties and BV/TV were increased when accounting only for the pedicle zone (R² = 0.66–0.94, p ≤ .032), but with weaker correlations for torsional loads (R² < 0.10). Our analyses highlight the role of local density in the pedicle zone on the fixation stiffness and strength of pedicle screws when pull-out loads are involved, but that local apparent bone density alone may not be sufficient to explain resistance in torsion.     

Efficiency considerations for the purely tapered interference fit (TIF) abutments used in dental implants

A tapered interference fit provides a mechanically reliable retention mechanism for the implant-abutment interface in a dental implant. Understanding the mechanical properties of the tapered interface with or without a screw at the bottom has been the subject of a considerable amount of studies involving experiments and finite element (FE) analysis. In this paper, approximate closed-form formulas are developed to analyze the mechanics of a tapered interference fit. In particular, the insertion force, the efficiency, defined as the ratio of the pull-out force to insertion force, and the critical insertion depth, which causes the onset of plastic deformation, are analyzed. It is shown that the insertion force is a function of the taper angle, the contact length, the inner and outer radii of the implant, the static and the kinetic coefficients of friction, and the elastic modulii of the implant/abutment materials. The efficiency of the tapered interference fit, which is defined as the ratio of the pull-out force to insertion force, is found to be greater than one, for taper angles that are less than 6 deg when the friction coefficient is 0.3. A safe range of insertion forces has been shown to exist. The lower end of this range depends on the maximum pull-out force that may occur due to occlusion in the multiple tooth restorations and the efficiency of the system; and the upper end of this range depends on the plastic deformation of the abutment and the implant due to interference fit. It has been shown that using a small taper angle and a long contact length widens the safe range of insertion forces.     

Do Implicit Attitudes Predict Actual Voting Behavior Particularly for Undecided Voters?

The prediction of voting behavior of undecided voters poses a challenge to psychologists and pollsters. Recently, researchers argued that implicit attitudes would predict voting behavior particularly for undecided voters whereas explicit attitudes would predict voting behavior particularly for decided voters. We tested this assumption in two studies in two countries with distinct political systems in the context of real political elections. Results revealed that (a) explicit attitudes predicted voting behavior better than implicit attitudes for both decided and undecided voters, and (b) implicit attitudes predicted voting behavior better for decided than undecided voters. We propose that greater elaboration of attitudes produces stronger convergence between implicit and explicit attitudes resulting in better predictive validity of both, and less incremental validity of implicit over explicit attitudes for the prediction of voting behavior. However, greater incremental predictive validity of implicit over explicit attitudes may be associated with less elaboration.     

Convergence behavior of high-resolution finite element models of trabecular bone

The convergence behavior of finite element models depends on the size of elements used, the element polynomial order, and on the complexity of the applied loads. For high-resolution models of trabecular bone, changes in architecture and density may also be important. The goal of this study was to investigate the influence of these factors on the convergence behavior of high-resolution models of trabecular bone. Two human vertebral and two bovine tibial trabecular bone specimens were modeled at four resolutions ranging from 20 to 80 microns and subjected to both compressive and shear loading. Results indicated that convergence behavior depended on both loading mode (axial versus shear) and volume fraction of the specimen. Compared to the 20 microns resolution, the differences in apparent Young's modulus at 40 microns resolution were less than 5 percent for all specimens, and for apparent shear modulus were less than 7 percent. By contrast, differences at 80 microns resolution in apparent modulus were up to 41 percent, depending on the specimen tested and loading mode. Overall, differences in apparent properties were always less than 10 percent when the ratio of mean trabecular thickness to element size was greater than four. Use of higher order elements did not improve the results. Tissue level parameters such as maximum principal strain did not converge. Tissue level strains converged when considered relative to a threshold value, but only if the strains were evaluated at Gauss points rather than element centroids. These findings indicate that good convergence can be obtained with this modeling technique, although element size should be chosen based on factors such as loading mode, mean trabecular thickness, and the particular output parameter of interest.     

Simulation of a knee joint replacement during a gait cycle using explicit finite element analysis.

The stress distribution within the polyethylene insert of a total knee joint replacement is dependent on the kinematics, which in turn are dependent on the design of the articulating surfaces, the relative position of the components and the tension of the surrounding soft tissues. Implicit finite element analysis techniques have been used previously to examine the polyethylene stresses. However, these have essentially been static analyses and hence ignored the influence of the kinematics. The aim of this work was to use an explicit finite element approach to simulate both the kinematics and the internal stresses within a single analysis. A simulation of a total knee joint replacement subjected to a single gait cycle within a knee wear simulator was performed and the results were compared with experimental data.The predicted kinematics were in close agreement with the experimental data. Various solution-dependent parameters were found to have little influence on the predicted kinematics. The predicted stresses were found to be dependent on the mesh density. This study has shown that an explicit finite element approach is capable of predicting the kinematics and the stresses within a single analysis at relatively low computational cost.     

A new material mapping procedure for quantitative computed tomography-based, continuum finite element analyses of the vertebra

Inaccuracies in the estimation of material properties and errors in the assignment of these properties into finite element models limit the reliability, accuracy, and precision of quantitative computed tomography (QCT)-based finite element analyses of the vertebra. In this work, a new mesh-independent, material mapping procedure was developed to improve the quality of predictions of vertebral mechanical behavior from QCT-based finite element models. In this procedure, an intermediate step, called the material block model, was introduced to determine the distribution of material properties based on bone mineral density, and these properties were then mapped onto the finite element mesh. A sensitivity study was first conducted on a calibration phantom to understand the influence of the size of the material blocks on the computed bone mineral density. It was observed that varying the material block size produced only marginal changes in the predictions of mineral density. Finite element (FE) analyses were then conducted on a square column-shaped region of the vertebra and also on the entire vertebra in order to study the effect of material block size on the FE-derived outcomes. The predicted values of stiffness for the column and the vertebra decreased with decreasing block size. When these results were compared to those of a mesh convergence analysis, it was found that the influence of element size on vertebral stiffness was less than that of the material block size. This mapping procedure allows the material properties in a finite element study to be determined based on the block size required for an accurate representation of the material field, while the size of the finite elements can be selected independently and based on the required numerical accuracy of the finite element solution. The mesh-independent, material mapping procedure developed in this study could be particularly helpful in improving the accuracy of finite element analyses of vertebroplasty and spine metastases, as these analyses typically require mesh refinement at the interfaces between distinct materials. Moreover, the mapping procedure is not specific to the vertebra and could thus be applied to many other anatomic sites.     

Numerical analysis of an anastomotic device

The explicit dynamic finite element method was utilized to investigate the deformation behaviour of a woven wire mesh tubular device that is used in a side-to-side anastomotic procedure for achieving gastrointestinal anastomosis. The numerical model was initially verified by comparison to experimental results that were obtained using a specialized testing mechanism. Once validated, the finite element model (FEM) was parameterized to ascertain the influence of several device parameters on its deformation behaviour. The importance of these parameters, as related to its optimal design for use in minimally invasive surgery (MIS), was subsequently ascertained and discussed.     

A finite element model of the human knee joint for the study of tibio-femoral contact

As a step towards developing a finite element model of the knee that can be used to study how the variables associated with a meniscal replacement affect tibio-femoral contact, the goals of this study were 1) to develop a geometrically accurate three-dimensional solid model of the knee joint with special attention given to the menisci and articular cartilage, 2) to determine to what extent bony deformations affect contact behavior, and 3) to determine whether constraining rotations other than flexion/extension affects the contact behavior of the joint during compressive loading. The model included both the cortical and trabecular bone of the femur and tibia, articular cartilage of the femoral condyles and tibial plateau, both the medial and lateral menisci with their horn attachments, the transverse ligament, the anterior cruciate ligament, and the medial collateral ligament. The solid models for the menisci and articular cartilage were created from surface scans provided by a noncontacting, laser-based, three-dimensional coordinate digitizing system with an root mean squared error (RMSE) of less than 8 microns. Solid models of both the tibia and femur were created from CT images, except for the most proximal surface of the tibia and most distal surface of the femur which were created with the three-dimensional coordinate digitizing system. The constitutive relation of the menisci treated the tissue as transversely isotropic and linearly elastic. Under the application of an 800 N compressive load at 0 degrees of flexion, six contact variables in each compartment (ie., medial and lateral) were computed including maximum pressure, mean pressure, contact area, total contact force, and coordinates of the center of pressure. Convergence of the finite element solution was studied using three mesh sizes ranging from an average element size of 5 mm by 5 mm to 1 mm by 1 mm. The solution was considered converged for an average element size of 2 mm by 2 mm. Using this mesh size, finite element solutions for rigid versus deformable bones indicated that none of the contact variables changed by more than 2% when the femur and tibia were treated as rigid. However, differences in contact variables as large as 19% occurred when rotations other than flexion/extension were constrained. The largest difference was in the maximum pressure. Among the principal conclusions of the study are that accurate finite element solutions of tibio-femoral contact behavior can be obtained by treating the bones as rigid. However, unrealistic constraints on rotations other than flexion/extension can result in relatively large errors in contact variables.     

A Voxel-based Formulation for Contact Finite Element Analysis

To date, voxel-based finite element models have not been feasible for contact problems, owing to the inherent stair-step boundary discontinuities. New preprocessing techniques are reported herein to mesh these boundaries smoothly, for purposes of contact stress analysis. Further, new techniques are reported to concentrate the mesh resolution automatically near the articular surface, thus reducing the problem size to levels compatible with executing nonlinear problems on contemporary engineering workstations. Close approximations to Hertzian analytical solutions were obtained for spherical and cylindrical geometries meshed in this manner, and an illustrative anatomical contact problem of the human hip joint is presented.     

Numerical Analysis of An Anastomotic Device

The explicit dynamic finite element method was utilized to investigate the deformation behaviour of a woven wire mesh tubular device that is used in a side-to-side anastomotic procedure for achieving gastrointestinal anastomosis. The numerical model was initially verified by comparison to experimental results that were obtained using a specialized testing mechanism. Once validated, the finite element model (FEM) was parameterized to ascertain the influence of several device parameters on its deformation behaviour. The importance of these parameters, as related to its optimal design for use in minimally invasive surgery (MIS), was subsequently ascertained and discussed.     

A three-dimensional finite element study on the biomechanical behavior of an FGBM dental implant in surrounding bone

Dental implants made of functionally graded biomaterials (FGBM) have been receiving increasing attention due to their unique advantage of being able to simultaneously satisfy biocompatibility, strength, corrosion resistance, etc., which a single composition with a uniform structure cannot satisfy. This paper investigates the biomechanical behavior of a threaded FGBM dental implant/surrounding bone system under static and harmonic occlusal forces by using a three-dimensional finite element method. The implant is a mixture of a bioceramic and a biometal with a smooth gradient in both the material composition and properties in the longitudinal direction. The interaction of the implant and the supporting bone tissues is considered. Three contact conditions at the implant-bone interface are used to model different osseointegration stages. A comprehensive parametric study is conducted to highlight the influence of the material properties, the volume fraction index, the occlusal force orientation, and the osseointegration quality on the maximum von-Mises stress, deformation distribution, natural frequencies, and harmonic response.     

Contact analysis of a posterior cervical spine plate using a three-dimensional canine finite element model

The development of a three-dimensional finite element model of a posteriorly plated canine cervical spine (C3-C6) including contact nonlinearities is described. The model was created from axial CT scans and the material properties were derived from the literature. The model demonstrated sufficient accuracy from the results of a mesh convergence test. Significant steps were taken toward establishing model validation by comparison of plate surface strains with a posteriorly plated canine cervical spine under three-point bending. This model was developed to better characterize the contact pressures at the various interfaces under average physiologic canine loading. The analysis showed that the screw-plate interfaces had the highest values of all the mechanical parameters evaluated.     

Mechanics of the tapered interference fit in dental implants

In evaluation of the long-term success of a dental implant, the reliability and the stability of the implant-abutment interface plays a great role. Tapered interference fits provide a reliable connection method between the abutment and the implant. In this work, the mechanics of the tapered interference fits were analyzed using a closed-form formula and the finite element (FE) method. An analytical solution, which is used to predict the contact pressure in a straight interference, was modified to predict the contact pressure in the tapered implant-abutment interface. Elastic-plastic FE analysis was used to simulate the implant and abutment material behavior. The validity and the applicability of the analytical solution were investigated by comparisons with the FE model for a range of problem parameters. It was shown that the analytical solution could be used to determine the pull-out force and loosening-torque with 5-10% error. Detailed analysis of the stress distribution due to tapered interference fit, in a commercially available, abutment-implant system was carried out. This analysis shows that plastic deformation in the implant limits the increase in the pull-out force that would have been otherwise predicted by higher interference values.     

Sampling efficiency in explicit and implicit membrane environments studied by peptide folding simulations

Efficient exploration of the conformational space of peptides embedded in biological membranes is vital to extract converged thermodynamic and kinetic data from simulation studies. Recently developed implicit membrane models promise vast increases in sampling efficiency compared to explicit membrane simulations, allowing for ab initio structure prediction and functional studies. In this study, a previously developed implicit membrane model, based on the generalized Born method, is compared to an explicit di‐palmitoyl‐phosphatidyl‐choline lipid bilayer and an octane slab membrane mimic. The complete folding process of a synthetic 16‐residue peptide is compared using these three setups. Since the comparison requires the entire folding pathway to be captured, individual simulations ranged up to 3 μs of MD. A quantitative sampling comparison using a wide range of performance metrics reveals that the implicit membrane model is at least 2 orders of magnitude more efficient than the simplest explicit setups. Proteins 2009. © 2008 Wiley‐Liss, Inc.