Techniques for mitigating the checkerboard formation: application in bone remodeling simulations

Bone remodelling prediction using mechanical stimulus with bone connectivity theory in porous implants

Strain energy density (SED) is considered to be the primary remodelling stimulus influencing the process of bone growth into porous implants. A bone remodelling algorithm incorporating the concept of bone connectivity, that newly formed bone should only grow from existing bone, was developed to provide a more biologically realistic simulation of bone growth. Results showed that the new algorithm prevented the occurrence of unconnected mature bone within porous implants, an unrealistic phenomenon observed using conventional adaptive elasticity theories. The bone connectivity algorithm had minimal effect (0.67% difference) on the final bone density distribution for standard bending and torsional moment cases. For a porous implant model, both algorithms, with and without bone connectivity implementation, reached the same final stiffness, with a difference of less than 0.01%. The bone connectivity algorithm predicted a slower and more gradual bone remodelling process, requiring at least 50% additional time for full remodelling compared to the conventional adaptive elasticity algorithm, which should be accounted for in the planning of rehabilitation strategies. The developed modelling can be employed to improve porous implant designs to achieve better clinical outcomes.

Biomechanical analysis of printable functionally graded material (FGM) dental implants for different bone densities

The long-term success of a dental implant is related to the material and design of the implant, and bone density. Conventional implants cause stress-shielding due to a mismatch between the implant and bone stiffness. Functionally graded porous materials and designs are a great choice for the design of implants to control the local stiffness at a certain location to meet the biomechanical requirements. The purpose of this study is to analyze five designs of axial and radial functionally graded materials (FGM) implants besides the conventional implant and conical and cylindrical shapes that were simulated with five different bone densities. The results showed that strain in bone increased with a decrease in cancellous bone density. The shape of the implant did not play an important role in strain/stress distribution. Conventional implants showed optimal strain (1000-2240 με) in low-density (0.7-0.8 g/cm3) bone, however, FGM implants produced optimal strain (990-1280 με) in the high-density bone (0.9-1 g/cm3) as compared to conventional implants. The proposed designs of FGM implants have the potential to address the complications of conventional implants in high-density bone.

Framework of sampling the subject-specific static loads from the gait cycle of interindividual variation

BACKGROUND AND OBJECTIVE

Numerous techniques for bone remodeling simulation have been developed based on Wolff's law. However, most studies have been conducted with empirically determined static loads, which cannot reflect subject-specific characteristics. We recently proposed a new concept of representative static loads (RSLs) to efficiently consider the effect of cyclically repeated dynamic loads on bone remodeling simulation. Based on this concept, the goal of this study is to sample the subject-specific static loads (SSL) from a general gait cycle of interindividual variation.

METHODS

A total of 15 gait cycles (ten normal and five abnormal cycles) obtained from the public database were used in this study. Each gait cycle was applied to a femur FE model constructed from the clinical CT scan data to evaluate the strain energy distribution as a reference. Then, a natural coordinate was introduced to maintain the predefined locations of extreme points (i.e., two peaks and one valley) for both normal and abnormal gait cycles. To determine the RSLs in the natural coordinate, five out of ten normal gait cycles were used. Through an inverse transformation for each gait cycle, the RSLs in the natural coordinate were converted to the SSLs in the original coordinate. Topology optimization results with the proposed SSLs were compared with those with a single full gait cycle (reference). For comparison, topology optimization was also conducted with empirically determined loads (EDLs) which have been widely used in the literature.

RESULTS

For normal gait cycles, the proposed SSLs reduced the average computing cost by 95.86% while suppressing the errors of bone mass distribution and apparent stiffness below maximum 4.24% and 1.72%, respectively. Even for abnormal gait cycles, the errors of bone mass distribution and apparent stiffness were suppressed below maximum 9.49% and 2.12%, respectively. Conversely, the conventional EDLs (peak loads selected in this study) showed significantly larger errors of maximum 47.28% and 30.31% in bone mass distribution and apparent stiffness for normal gait cycles.

CONCLUSION

By virtue of using the coordinate transformation for each gait cycle, the proposed SSLs achieved a higher accuracy in the bone mass distribution and apparent stiffness than the previous RSLs and EDLs. Furthermore, this approach can be used for abnormal gait cycles which have higher interindividual variation.