Physical modeling with orthotropic material based on harmonic fields

OpenMandible: An open-source framework for highly realistic numerical modelling of lower mandible physiology

OBJECTIVE

Computer modeling of lower mandible physiology remains challenging because prescribing realistic material characteristics and boundary conditions from medical scans requires advanced equipment and skill sets. The objective of this study is to provide a framework that could reduce simplifications made and inconsistency (in terms of geometry, materials, and boundary conditions) among further studies on the topic.

METHODS

The OpenMandible framework offers: 1) the first publicly available multiscale model of the mandible developed by combining cone beam computerized tomography (CBCT) and μCT imaging modalities, and 2) a C++ software tool for the generation of simulation-ready models (tet4 and hex8 elements). In addition to the application of conventional (Neumann and Dirichlet) boundary conditions, OpenMandible introduces a novel geodesic wave propagation - based approach for incorporating orthotropic micromechanical characteristics of cortical bone, and a unique algorithm for modeling muscles as uniformly directed vectors. The base intact model includes the mandible (spongy and compact bone), 14 teeth (comprising dentin, enamel, periodontal ligament, and pulp), simplified temporomandibular joints, and masticatory muscles (masseter, temporalis, medial, and lateral pterygoid).

RESULTS

The complete source code, executables, showcases, and sample data are freely available on the public repository: https://github.com/ArsoVukicevic/OpenMandible. It has been demonstrated that by slightly editing the baseline model, one can study different "virtual" treatments or diseases, including tooth restoration, placement of implants, mandible bone degradation, and others.

SIGNIFICANCE

OpenMandible eases the community to undertake a broad range of studies on the topic, while increasing their consistency and reproducibility. At the same time, the needs for dedicated equipment and skills for developing realistic simulation models are significantly reduced.

Semi-automatic segmentation of femur based on harmonic barrier

BACKGROUND AND OBJECTIVE

Segmentation of the femur from the hip joint in computed tomography (CT) is an important preliminary step in hip surgery planning and simulation. However, this is a time-consuming and challenging task due to the weak boundary, the varying topology of the hip joint, and the extremely narrow or blurred space between the femoral head and the acetabulum. To address these problems, this study proposed a semi-automatic segmentation framework based on harmonic fields for accurate segmentation.

METHODS

The proposed method comprises three steps. First, with high-level information provided by the user, shape information provided by neighboring slices as well as the statistical information in the mask, a region selection method is proposed to effectively locate joint space for the harmonic field. Second, incorporated with an improved gradient, the harmonic field is used to adaptively extract a curve as the barrier that separates the femoral head from the acetabulum accurately. Third, a divide and conquer segmentation strategy based on the harmonic barrier is used to combine the femoral head part and body part as the final segmentation result.

RESULTS

We have tested 40 hips with considerately narrow or disappeared joint spaces. The experimental results are evaluated based on Jaccard, Dice, directional cut discrepancy (DCD) and receiver operating characteristic (ROC), and we achieve the higher Jaccard of 84.02%, Dice of 85.96%, area under curve (AUC) of 89.3%, and the lower error with DCD of 0.52mm. The effective ratio of our method is 79.1% even for cases with severe malformation. The results show that our method performs best in terms of effectiveness and accuracy on the whole data set.

CONCLUSIONS

The proposed method is efficient to segment femurs with narrow joint space. The accurate segmentation results can assist the physicians for osteoarthritis diagnosis in future.

Anisotropic Finite Element Modeling Based on a Harmonic Field for Patient-Specific Sclera

Purpose. This study examined the influence of anisotropic material for human sclera. Method. First, the individual geometry of patient-specific sclera was reproduced from a laser scan. Then, high quality finite element modeling of individual sclera was performed using a convenient automatic hexahedral mesh generator based on harmonic field and integrated with anisotropic material assignment function. Finally, comparison experiments were designed to investigate the effects of anisotropy on finite element modeling of sclera biomechanics. Results. The experimental results show that the presented approach can generate high quality anisotropic hexahedral mesh for patient-specific sclera. Conclusion. The anisotropy shows significant differences for stresses and strain distribution and careful consideration should be given to its use in biomechanical FE studies.

Performance Analysis of a Fiber Reinforced Plastic Oil Cooler Cover Considering the Anisotropic Behavior of the Fiber Reinforced PA66

In this paper, a simulation method based on an orthogonal anisotropic material is proposed. A numerical example using a simple plate is presented to show the difference in the static performance between the orthogonal anisotropic and the isotropic models. Comparing with the tested modal data of a diesel engine oil cooler cover made by glass fiber reinforced polyamide 66 (PA66), the proposed simulation method was confirmed to be much closer to reality than the general isotropic model. After that, a comprehensive performance comparison between the plastic oil cooler covers with the orthogonal anisotropic and the isotropic fiber orientations was carried out including a static deformation and stress analysis under a pressure-temperature coupled load, a forced response analysis, and an acoustic analysis under real operating conditions. The results show that the stress, the deformation, the peak vibration velocity, and the overall sound power level of the orthogonal anisotropic model are different from that obtained with the isotropic model. More importantly, the proposed method can provide a much more detailed frequency content compared to the isotropic model.

Automatic Tooth Segmentation of Dental Mesh Based on Harmonic Fields

An important preprocess in computer-aided orthodontics is to segment teeth from the dental models accurately, which should involve manual interactions as few as possible. But fully automatic partition of all teeth is not a trivial task, since teeth occur in different shapes and their arrangements vary substantially from one individual to another. The difficulty is exacerbated when severe teeth malocclusion and crowding problems occur, which is a common occurrence in clinical cases. Most published methods in this area either are inaccurate or require lots of manual interactions. Motivated by the state-of-the-art general mesh segmentation methods that adopted the theory of harmonic field to detect partition boundaries, this paper proposes a novel, dental-targeted segmentation framework for dental meshes. With a specially designed weighting scheme and a strategy of a priori knowledge to guide the assignment of harmonic constraints, this method can identify teeth partition boundaries effectively. Extensive experiments and quantitative analysis demonstrate that the proposed method is able to partition high-quality teeth automatically with robustness and efficiency.

Interactive tooth partition of dental mesh base on tooth-target harmonic field

The accurate tooth partition of dental mesh is a crucial step in computer-aided orthodontics. However, tooth boundary identification is not a trivial task for tooth partition, since different shapes and their arrangements vary substantially among common clinical cases. Though curvature field is traditionally used for identifying boundaries, it is normally not reliable enough. Other methods may improve the accuracy, but require intensive user interaction. Motivated by state-of-the-art general interactive mesh segmentation methods, this paper proposes a novel tooth-target partition framework that employs harmonic fields to partition teeth accurately and effectively. In addition, a refining strategy is introduced to successfully segment teeth from the complicated dental model with indistinctive tooth boundaries on its lingual side surface, addressing an issue that had not been solved properly before. To utilise high-level information provided by the user, smart and intuitive user interfaces are also proposed with minimum interaction. In fact, most published interactive methods specifically designed for tooth partition are lacking efficient user interfaces. Extensive experiments and quantitative analyses show that our tooth partition method outperforms the state-of-the-art approaches in terms of accuracy, robustness and efficiency.

Influence of orthotropy on biomechanics of peri-implant bone in complete mandible model with full dentition

Objective. The study was to investigate the impact of orthotropic material on the biomechanics of dental implant, based on a detailed mandible with high geometric and mechanical similarity. Materials and Methods. Multiple data sources were used to elaborate detailed biological structures and implant CAD models. In addition, an extended orthotropic material assignment methodology based on harmonic fields was used to handle the alveolar ridge region to generate compatible orthotropic fields. The influence of orthotropic material was compared with the commonly used isotropic model and simplified orthotropic model. Results. The simulation results showed that the values of stress and strain on the implant-bone interface almost increased in the orthotropic model compared to the isotropic case, especially for the cancellous bone. However, the local stress concentration was more obvious in the isotropic case compared to that in orthotropic case. The simple orthotropic model revealed irregular stress and strain distribution, compared to the isotropic model and the real orthotropic model. The influence of orthotropy was little on the implant, periodontal ligament, tooth enamel, and dentin. Conclusion. The orthotropic material has significant effect on stress and strain of implant-bone interface in the mandible, compared with the isotropic simulation. Real orthotropic mechanical properties of mandible should be emphasized in biomechanical studies of dental implants.