The application of topology optimization on the quantitative description of the external shape of bone structure

Computational model of articular cartilage regeneration induced by scaffold implantation in vivo

Computational models allow to explain phenomena that cannot be observed through an animal model, such as the strain and stress states which can highly influence regeneration of the tissue. For this purpose, we have developed a simulation tool to determine the mechanical conditions provided by the polymeric scaffold. The computational model considered the articular cartilage, the subchondral bone, and the scaffold. All materials were modeled as poroelastic, and the cartilage had linear-elastic oriented collagen fibers. This model was able to explain the remodeling process that subchondral bone goes through, and how the scaffold allowed the conditions for cartilage regeneration. These results suggest that the use of scaffolds might lead the cartilaginous tissue growth in vivo by providing a better mechanical environment. Moreover, the developed computational model demonstrated to be useful as a tool prior experimental in vivo studies, by predicting the possible outcome of newly proposed treatments allowing to discard approaches that might not bring good results.

A tool for solving bone growth related problems using finite elements adaptive meshes

Long bones geometry changes in response to longitudinal growth in the epiphyseal plates and hydroxyapatite apposition in the periosteum. Due to its relevance for growth modulation and orthotics performance, researchers have extensively modeled these phenomena, using the finite elements method for it almost since the introduction of modern computers. This is a rather complex task that, besides the inherent difficulty of solving the models equations, requires considering a moving boundary. Here, the development of a new computational tool for its resolution is described. A generalized formulation of these problems is established based on the most common approaches taken in the literature and a novel finite elements algorithm is proposed for its resolution. The later allows a significant reduction of the spatial discretization requirements, the computational cost and the numerical errors associated with more classical approaches. The potentiality of the method is demonstrated by its application to three cases of practical interest, namely, hemiepiphysiodesis treatment, growth in the distal femur and bone remodeling around hip prosthesis. Eight relevant cases of study and an open source implementation of the proposed algorithm are also provided as supplementary material.

Three-dimensional topology optimization model to simulate the external shapes of bone

Elucidation of the mechanism by which the shape of bones is formed is essential for understanding vertebrate development. Bones support the body of vertebrates by withstanding external loads, such as those imposed by gravity and muscle tension. Many studies have reported that bone formation varies in response to external loads. An increased external load induces bone synthesis, whereas a decreased external load induces bone resorption. This relationship led to the hypothesis that bone shape adapts to external load. In fact, by simulating this relationship through topology optimization, the internal trabecular structure of bones can be successfully reproduced, thereby facilitating the study of bone diseases. In contrast, there have been few attempts to simulate the external structure of bones, which determines vertebrate morphology. However, the external shape of bones may be reproduced through topology optimization because cells of the same type form both the internal and external structures of bones. Here, we constructed a three-dimensional topology optimization model to attempt the reproduction of the external shape of teleost vertebrae. In teleosts, the internal structure of the vertebral bodies is invariable, exhibiting an hourglass shape, whereas the lateral structure supporting the internal structure differs among species. Based on the anatomical observations, we applied different external loads to the hourglass-shaped part. The simulations produced a variety of three-dimensional structures, some of which exhibited several structural features similar to those of actual teleost vertebrae. In addition, by adjusting the geometric parameters, such as the width of the hourglass shape, we reproduced the variation in the teleost vertebrae shapes. These results suggest that a simulation using topology optimization can successfully reproduce the external shapes of teleost vertebrae. By applying our topology optimization model to various bones of vertebrates, we can understand how the external shape of bones adapts to external loads.

In this paper, we developed a computational method to investigate the relationship between three-dimensional bone shape and external loads imposed on bones. Many studies report that bone formation varies in response to external loads. An increased external load induces bone synthesis, whereas a decreased external load induces bone resorption. This relationship led to the hypothesis that the shape of bones adapts to external load. However, it remains unclear whether this hypothesis can explain the shape of bones. Here, we constructed a three-dimensional mathematical model that imitates the cellular activities of bone formation to attempt the reproduction of the shape of teleost vertebrae. In teleosts, the shape of the vertebrae differs among the species. We set the multiple types of external load conditions in the simulations and compared the simulation results with different teleost vertebrae. The produced structures that can resist the deformation of the surrounding tissues exhibited multiple structural features similar to the vertebrae of several teleost species. This result shows that the formation of bone shape can be explained by the adaptation to external load.

Exploring conditions that make cortical bone geometry optimal for physiological loading

While physiological loading on lower long bones changes during bone development, the bone cross section either remains circular or slowly changes from nearly circular to other shapes such as oval and roughly triangular. Bone is said to be an optimal structure, where strength is maximized using the optimal distribution of bone mass (also called Wolff's law). One of the most appropriate mathematical validations of this law would be a structural optimization-based formulation where total strain energy is minimized against a mass and a space constraint. Assuming that the change in cross section during bone development and homeostasis after adulthood is direct result of the change in physiological loading, this work investigates what optimization problem formulation (collectively, design variables, objective function, constraints, loading conditions, etc.) results in mathematically optimal solutions that resemble bones under actual physiological loading. For this purpose, an advanced structural optimization-based computational model for cortical bone development and defect repair is presented. In the optimization problem, overall bone stiffness is maximized first against a mass constraint, and then also against a polar first moment of area constraint that simultaneously constrains both mass and space. The investigation is completed in two stages. The first stage is developmental stage when physiological loading on lower long bones (tibia) is a random combination of axial, bending and torsion. The topology optimization applied to this case with the area moment constraint results into circular and elliptical cross sections similar to that found in growing mouse or human. The second investigation stage is bone homeostasis reached in adulthood when the physiological loading has a fixed pattern. A drill hole defect is applied to the adult mouse bone, which would disrupt the homeostasis. The optimization applied after the defect interestingly brings the damaged section back to the original intact geometry. The results, however, show that cortical bone geometry is optimal for the physiological loading only when there is also a constraint on polar moment of area. Further numerical experiments show that application of torsion along with the gait-analysis-based physiological loading improves the results, which seems to indicate that the cortical bone geometry is optimal for some amount of torsion in addition to the gait-based physiological loading. This work has a potential to be extended to bone growth/development models and fracture healing models, where topology optimization and polar moment of area constraint have not been introduced earlier.

Design of complex bone internal structure using topology optimization with perimeter control

Large facial bone loss usually requires patient-specific bone implants to restore the structural integrity and functionality that also affects the appearance of each patient. Titanium alloys (e.g., Ti-6Al-4V) are typically used in the interfacial porous coatings between the implant and the surrounding bone to promote stability. There exists a property mismatch between the two that in general leads to complications such as stress-shielding. This biomechanical discrepancy is a hurdle in the design of bone replacements. To alleviate the mismatch, the internal structure of the bone replacements should match that of the bone. Topology optimization has proven to be a good technique for designing bone replacements. However, the complex internal structure of the bone is difficult to mimic using conventional topology optimization methods without additional restrictions. In this work, the complex bone internal structure is recovered using a perimeter control based topology optimization approach. By restricting the solution space by means of the perimeter, the intricate design complexity of bones can be achieved. Three different bone regions with well-known physiological loadings are selected to illustrate the method. Additionally, we found that the target perimeter value and the pattern of the initial distribution play a vital role in obtaining the natural curvatures in the bone internal structures as well as avoiding excessive island patterns.

Simulation on the internal structure of three-dimensional proximal tibia under different mechanical environments

Background

Bone can adjust its morphological structure to adapt to the changes of mechanical environment, i.e. the bone structure change is related to mechanical loading. This implies that osteoarthritis may be closely associated with knee joint deformity. The purposes of this paper were to simulate the internal bone mineral density (BMD) change in three-dimensional (3D) proximal tibia under different mechanical environments, as well as to explore the relationship between mechanical environment and bone morphological abnormity.

Methods

The right proximal tibia was scanned with CT to reconstruct a 3D proximal tibia model in MIMICS, then it was imported to finite element software ANSYS to establish 3D finite element model. The internal structure of 3D proximal tibia of young normal people was simulated using quantitative bone remodeling theory in combination with finite element method, then based on the changing pattern of joint contact force on the tibial plateau in valgus knees, the mechanical loading was changed, and the simulated normal tibia structure was used as initial structure to simulate the internal structure of 3D proximal tibia for old people with 6° valgus deformity. Four regions of interest (ROIs) were selected in the proximal tibia to quantitatively analyze BMD and compare with the clinical measurements.

Results

The simulation results showed that the BMD distribution in 3D proximal tibia was consistent with clinical measurements in normal knees and that in valgus knees was consistent with the measurement of patients with osteoarthritis in clinics.

Conclusions

It is shown that the change of mechanical environment is the main cause for the change of subchondral bone structure, and being under abnormal mechanical environment for a long time may lead to osteoarthritis. Besides, the simulation method adopted in this paper can more accurately simulate the internal structure of 3D proximal tibia under different mechanical environments. It helps to better understand the mechanism of osteoarthritis and provides theoretical basis and computational method for the prevention and treatment of osteoarthritis. It can also serve as basis for further study on periprosthetic BMD changes after total knee arthroplasty, and provide a theoretical basis for optimization design of prosthesis.

Influence of mastication and edentulism on mandibular bone density

The aim of this study was to demonstrate that external loading due to daily activities, including mastication, speech and involuntary open-close cycles of the jaw contributes to the internal architecture of the mandible. A bone remodelling algorithm that regulates the bone density as a function of stress and loading cycles is incorporated into finite element analysis. A three-dimensional computational model is constructed on the basis of computerised tomography (CT) images of a human mandible. Masticatory muscle activation involved during clenching is modelled by static analysis using linear optimisation. Other loading conditions are approximated by imposing mandibular flexure. The simulations predict that mandibular bone density distribution results in a tubular structure similar to what is observed in the CT images. Such bone architecture is known to provide the bone optimum strength to resist bending and torsion during mastication while reducing the bone mass. The remodelling algorithm is used to simulate the influence of edentulism on mandibular bone loss. It is shown that depending on the location and number of missing teeth, up to one-third of the mandibular bone mass can be lost due to lack of adequate mechanical stimulation.

Parametric investigation of load-induced structure remodeling in the proximal femur

The process of adaptive bone remodeling can be simulated with a self-optimizing finite element method. The basic remodeling rules attempt to obtain a constant value for the strain energy per unit bone mass, by adapting density. The precise solution is dependent on the loads, initial conditions, and the parameters of the remodeling rule. While there are several investigations on developing algorithms leading to the bone density distribution in the proximal femur, these algorithms often require a large number of iterations. The aim of this study was to develop a more efficient adaptive bone remodeling algorithm, and to identify how the bone density distribution of the proximal femur was affected by parameters that govern the remodeling process. The forces at different phases of the gait cycle were applied as boundary conditions. The bone density distributions from these forces were averaged to estimate the density distribution in the proximal femur. The effect of varying the initial bone density, spatial influence function, non-linear order of the adaptive algorithm, and the influence range on the converged solution were investigated. The proposed procedure was shown to converge in a fewer number of iterations and requiring less computational time, while still generating a realistic bone density distribution. It was also shown that varying the identified parameters within reasonable upper and lower bounds had very little impact on the qualitative form of the converged solution. In contrast, the convergence rate was affected to a greater degree by variation of these parameters. In all cases, the solutions obtained are comparable with the actual density in the proximal femur, as measured by Dual-energy X-ray absorptiometry (DEXA) scans.

Simulated evolution of the vertebral body based on basic multicellular unit activities

A numerical model based on the theory of bone remodeling is proposed to predict the evolution of trabecular bone architecture within the vertebral body and to investigate the process of degeneration in vertebral bone. In this study, particular attention is paid on the description of microstructure changes during the aging process. To take into account the effect of basic multicellular units (BMUs), a set of computational algorithms has been developed. It is assumed that BMU activation probability depends on the state of damaged bone tissue (damage accumulation, ω), which is evaluated according to previous research concerning bone fatigue damage. Combining these algorithms with the finite-element method (FEM), the microstructure of vertebral bone has been predicted for up to 8 simulated years. Moreover, biomechanical material properties have been monitored to investigate the changes of vertebral bone with age. This study shows that the simulation based on BMU activities has the potential to define and predict the morphological evolution of the vertebral body. It can be concluded that the novel algorithms incorporating the coupled effects of both adaptive remodeling and microdamage remodeling could be utilized to gain greater insight into the mechanism of bone loss in the elderly population.

Prediction of shape and internal structure of the proximal femur using a modified level set method for structural topology optimisation

A computational framework was developed to simulate the bone remodelling process as a structural topology optimisation problem. The mathematical formulation of the Level Set technique was extended and then implemented into a coronal plane model of the proximal femur to simulate the remodelling of internal structure and external geometry of bone into the optimal state. Results indicated that the proposed approach could reasonably mimic the major geometrical and material features of the natural bone. Simulation of the internal bone remodelling on the typical gross shape of the proximal femur, resulted in a density distribution pattern with good consistency with that of the natural bone. When both external and internal bone remodelling were simulated simultaneously, the initial rectangular design domain with a regularly distributed mass reduced gradually to an optimal state with external shape and internal structure similar to those of the natural proximal femur.

Computational simulation of simultaneous cortical and trabecular bone change in human proximal femur during bone remodeling

In this study, we developed a numerical framework that computationally determines simultaneous and interactive structural changes of cortical and trabecular bone types during bone remodeling, and we investigated the structural correlation between the two bone types in human proximal femur. We implemented a surface remodeling technique that performs bone remodeling in the exterior layer of the cortical bone while keeping its interior area unchanged. A micro-finite element (microFE) model was constructed that represents the entire cortical bone and full trabecular architecture in human proximal femur. This study simulated and compared the bone adaptation processes of two different structures: (1) femoral bone that has normal cortical bone shape and (2) perturbed femoral bone that has an artificial bone lump in the inferomedial cortex. Using the proposed numerical method in conjunction with design space optimization, we successfully obtained numerical results that resemble actual human proximal femur. The results revealed that actual cortical bone, as well as the trabecular bone, in human proximal femur has structurally optimal shapes, and it was also shown that a bone abnormality that has little contribution to bone structural integrity tends to disappear. This study also quantitatively determined the structural contribution of each bone: when the trabecular adaptation was complete, the trabecular bone supported 54% of the total load in the human proximal femur while the cortical bone carried 46%.

The effects of the spatial influence function on orthotropic femur remodelling

The morphology and internal structure of bone are modulated by the mechanical stimulus. The osteocytes can sense the stimulus signals from the adjacent regions and respond to them through bone growth or bone absorption. This mechanism can be modelled as the spatial influence function (SIF) in bone adaptation algorithm. In this paper, the remodelling process was simulated in human femurs using an adaptation algorithm with and without SIF, and the trabecular bone was assumed to be orthotropic. A different influence radius and weighting factor were adopted to study the effects of the SIF on the bone density distribution and trabecular alignment. The results have shown that the mean density and L—T ratio (the ratio of longitudinal modulus to transverse modulus) had an excellent linear relationship with the weighting factor when the influence radius was small. The characteristics of density distribution and L—T ratio accorded with the actual observation or measurement when a small weighting factor was used. The large influence radius and weighting factor led to unrealistic results. In contrast, the SIF hardly affected the trabecular alignment, as the mean variation angles of principal axes were less than 1.0 degree for any influence radius and weighting factor

Modelling external bone adaptation using evolutionary structural optimisation

External remodelling is significant in the bone healing process, and it is essential to predict the bone external shape in the design of artificial bone grafts. This paper demonstrates the effectiveness of the evolutionary structural optimisation (ESO) method for the simulation of bone morphology. A two-dimensional ESO strategy is developed which is capable of finding the modified bone topology beginning with any geometry under any loading conditions. The morphology of bone structure is described by the quantitative bone adaptation theory, which is integrated with the finite element method. The evolutionary topology optimisation process is introduced to find the bone shape. A rectangle, which occupies a larger space than the external shape of the bone structure, is specified as a design domain; the evolutionary process iteratively eliminates and redistributes material throughout the domain to obtain an optimum arrangement of bone materials. The technique has been tested on a wide range of examples. In this paper, the formation of trabecular bone architecture around an implant is studied; as another example, the growth of the coronal section of a vertebral body is predicted. The examples support the assertion that the external shape of bone structure can be successfully predicted by the proposed ESO procedure.

The numerical simulation of osteophyte formation on the edge of the vertebral body using quantitative bone remodeling theory

OBJECTIVES

To extend the quantitative prediction of the external shape of bone structure to the simulation of osteophyte formation on the edge of vertebral body.

METHODS

The high-order nonlinear equation of bone remodeling was used to control osteophyte formation process. The idea of topology optimization in engineering was imported to allow the outgrowth of osteophyte. Osteophyte on the edge of a vertebral body in its mid-sagittal plane was simulated numerically.

RESULTS

Osteophyte formation is an adaptive bone remodeling process in response to the progressive changes of mechanical environment, which were mainly caused by intervertebral disc degeneration.

CONCLUSION

The numerical simulation in this paper extended the structural simulation based on the quantitative bone remodeling theory to bone morphological abnormity in orthopaedics, which can help to better understand the relationship between bone morphological abnormity and the mechanical environment.