Gongfa Chen

A new approach for assigning bone material properties from CT images into finite element models

Generation of subject-specific finite element (FE) models from computed tomography (CT) datasets is of significance for application of the FE analysis to bone structures. A great challenge that remains is the automatic assignment of bone material properties from CT Hounsfield Units into finite element models. This paper proposes a new assignment approach, in which material properties are directly assigned to each integration point. Instead of modifying the dataset of FE models, the proposed approach divides the assignment procedure into two steps: generating the data file of the image intensity of a bone in a MATLAB program and reading the file into ABAQUS via user subroutines. Its accuracy has been validated by assigning the density of a bone phantom into a FE model. The proposed approach has been applied to the FE model of a sheep tibia and its applicability tested on a variety of element types. The proposed assignment approach is simple and illustrative. It can be easily modified to fit users situations.

Computational investigations of mechanical failures of internal plate fixation.

Abstract This paper investigated the biomechanics of two clinical cases of bone fracture treatments. Both fractures were treated with the same locking compression plate but with different numbers of screws as well as different plate materials. The fracture treated with 12 screws (rigid fixation) failed at 7 weeks with the plate breaking; the fracture with six screws (flexible fixation) endured the entire healing process. It was hypothesized that the plate failure in the unsuccessful case was due to the material fatigue induced by stress concentration in the plate. As the two clinical cases had different fracture locations and different plate materials, finite element simulations were undertaken for each fractured bone fixed by both a rigid and a flexible method. This enabled comparisons to be made between the rigid and flexible fixation methods. The fatigue life was assessed for each fixation method. The results showed that the stress in the rigid fixation methods could be significantly higher than that in flexible fixation methods. The fatigue analyses showed that, with the stress level in flexible fixation (i.e. with fewer screws), the plate was able to endure 2000 days, and that the plate in rigid fixation could fail by fatigue fracture in 20 days. The paper concludes that the rigid fixation method resulted in serious stress concentrations in the plate, which induced fatigue failure. The flexible fixation gave sufficient stability and was better for fracture healing.

Tooth fracture risk analysis based on a new finite element dental structure models using micro-CT data

The finite element (FE) analysis is an effective method to study the strength and predict the fracture risk of endodontically-treated teeth. This paper presents a rapid method developed to generate a comprehensive tooth FE model using data retrieved from micro-computed tomography (μCT). With this method, the inhomogeneity of material properties of teeth was included into the model without dividing the tooth model into different regions. The material properties of the tooth were assumed to be related to the mineral density. The fracture risk at different tooth portions was assessed for root canal treatments. The micro-CT images of a tooth were processed by a Matlab software programme and the CT numbers were retrieved. The tooth contours were obtained with thresholding segmentation using Amira. The inner and outer surfaces of the tooth were imported into Solidworks and a three-dimensional (3D) tooth model was constructed. An assembly of the tooth model with the periodontal ligament (PDL) layer and surrounding bone was imported into ABAQUS. The material properties of the tooth were calculated from the retrieved CT numbers via ABAQUS users subroutines. Three root canal geometries (original and two enlargements) were investigated. The proposed method in this study can generate detailed 3D finite element models of a tooth with different root canal enlargements and filling materials, and would be very useful for the assessment of the fracture risk at different tooth portions after root canal treatments.

Computational approach to localization using global energy minimization

Damage localization induced by strain softening can be predicted by the direct minimization of a global energy function. This article concerns the computational strategy for implementing this principle for softening materials such as concrete. Instead of using heuristic global optimization techniques, our strategies are a hybrid of local optimization methods with a path-finding approach to ensure a global optimum. With admissible nodal displacements being independent variables, it is easy to deal with the geometric (mesh) constraint conditions. The direct search optimization methods recover the localized solutions for a range of softening lattice models which are representative of quasi-brittle structures

Mechanobiology of Bone Development and Computational Simulations

Bone is important because it provides the skeleton structural integrity and enables movement and locomotion. Its development and morphology follow its function. It adapts to changes of mechanical loading and has the ability to repair itself after damage or fracture. The processes of bone development, bone adaptation, and bone regeneration in fracture healing are regulated, in part, by mechanical stimuli that result when the bone is loaded.

Parametric Investigations on Columns with Un-bonded Reinforcing Bars for Damage Control

Influential factors on concrete cracking in column ends are studied by a non-linear finite element model. Following the previous strategy developed in the Building Research Institute in Japan, un-bonded zones are introduced into column ends and additional re-bars added around the un-bonded zones. The parametric studies of the length of the un-bonded zone and the amount of the additional re-bars are conducted. The load-displacement curves and cracking patterns in concrete are obtained and compared for all of the considered designs. The parametric studies show that the length of the un-bonded zone has significant influences on concrete cracking in column ends; the longer the un-bonded zone is, the better the concrete cracking is controlled in column ends. The studies also show that the additional re-bars slightly influence the load-displacement curves and cracking patterns. The influences of the bond-slip on the numerical analysis are also investigated and the proper parameters of the bond-slip are determined and recommended by the parametric studies.

PREDICTING THE FATIGUE LIFE OF INTERNAL FRACTURE FIXATION PLATES

Failures of fracture fixation plates, often related to fatigue fractures of the implants, have been reported (Banovetz et al, 1996). While metallurgical defects can usually be excluded, many of these fractures can be explained with the biomechanical situation. This study investigated the biomechanics of two clinical cases, both of which used a 14-hole locking compression plate. In the first case, a titanium plate was used in a rigid configuration with 12 screws resulting in plate breakage after 7 weeks (Sommer et al, 2003). In the second case, a stainless steel plate, which endured the entire healing process, was used in a flexible application with only 6 screws.

One-Dimensional Nonlinear Model for Prediction of Crack Spacing in Concrete Pavements

This paper proposes a one-dimensional non-linear model to predict the minimum and maximum crack spacings due to shrinkage in concrete pavements. The proposed model consists of two cohesive cracks and an elastic bar restrained by distributed elastic springs. The cohesive crack is characterized by an exponential softening constitutive relation. A set of non-linear equilibrium conditions are obtained. By varying the length of the elastic bar of the proposed model, the tensile forces acting on the cohesive cracks and the energy profiles are investigated. It is demonstrated that the cracking pattern varies with the length of the elastic bar (i.e. the spacing between the two possible cracks), from which the minimum and maximum crack spacings are obtained. Numerical analyses are made of a model pavement and the results indicate that it is the energy minimization principle that governs the cracking pattern. The proposed model provides physical insight into the mechanism of crack spacing in concrete pavements.