3-D finite element modelling of facial soft tissue and preliminary application in orthodontics

Prediction of soft tissue aesthetics is important for achieving an optimal outcome in orthodontic treatment planning. Previously, applicable procedures were mainly restricted to 2-D profile prediction. In this study, a generic 3-D finite element (FE) model of the craniofacial soft and hard tissue was constructed, and individualisation of the generic model based on cone beam CT data and mathematical transformation was investigated. The result indicated that patient-specific 3-D facial FE model including different layers of soft tissue could be obtained through mathematical model transformation. Average deviation between the transformed model and the real reconstructed one was 0.47 ± 0.77 mm and 0.75 ± 0.84 mm in soft and hard tissue, respectively. With boundary condition defined according to treatment plan, such FE model could be used to predict the result of orthodontic treatment on facial soft tissue.     

An eFTD-VP framework for efficiently generating patient-specific anatomically detailed facial soft tissue FE mesh for craniomaxillofacial surgery simulation

Accurate surgical planning and prediction of craniomaxillofacial surgery outcome requires simulation of soft tissue changes following osteotomy. This can only be achieved by using an anatomically detailed facial soft tissue model. The current state-of-the-art of model generation is not appropriate to clinical applications due to the time-intensive nature of manual segmentation and volumetric mesh generation. The conventional patient-specific finite element (FE) mesh generation methods are to deform a template FE mesh to match the shape of a patient based on registration. However, these methods commonly produce element distortion. Additionally, the mesh density for patients depends on that of the template model. It could not be adjusted to conduct mesh density sensitivity analysis. In this study, we propose a new framework of patient-specific facial soft tissue FE mesh generation. The goal of the developed method is to efficiently generate a high-quality patient-specific hexahedral FE mesh with adjustable mesh density while preserving the accuracy in anatomical structure correspondence. Our FE mesh is generated by eFace template deformation followed by volumetric parametrization. First, the patient-specific anatomically detailed facial soft tissue model (including skin, mucosa, and muscles) is generated by deforming an eFace template model. The adaptation of the eFace template model is achieved by using a hybrid landmark-based morphing and dense surface fitting approach followed by a thin-plate spline interpolation. Then, high-quality hexahedral mesh is constructed by using volumetric parameterization. The user can control the resolution of hexahedron mesh to best reflect clinicians’ need. Our approach was validated using 30 patient models and 4 visible human datasets. The generated patient-specific FE mesh showed high surface matching accuracy, element quality, and internal structure matching accuracy. They can be directly and effectively used for clinical simulation of facial soft tissue change.     

A novel method for prediction of dynamic smiling expressions after orthodontic treatment: a case report

Smile esthetics has become increasingly important for orthodontic patients, thus prediction of post-treatment smile is necessary for a perfect treatment plan. In this study, with a combination of three-dimensional craniofacial data from the cone beam computed tomography and color-encoded structured light system, a novel method for smile prediction was proposed based on facial expression transfer, in which dynamic facial expression was interpreted as a matrix of facial depth changes. Data extracted from the pre-treatment smile expression record were applied to the post-treatment static model to realize expression transfer. Therefore smile esthetics of the patient after treatment could be evaluated in pre-treatment planning procedure. The positive and negative mean values of error for prediction accuracy were 0.9 and ? 1.1 mm respectively, with the standard deviation of ± 1.5 mm, which is clinically acceptable. Further studies would be conducted to reduce the prediction error from both the static and dynamic sides as well as to explore automatically combined prediction from the two sides.     

Finite element-based force/moment-driven simulation of orthodontic tooth movement

The objectives of this study were to develop a numerically controlled experimental set-up to predict the movement caused by the force systems of orthodontic devices and to experimentally verify this system. The presented experimental set-up incorporated an artificial tooth fixed via a 3D force/moment sensor to a parallel kinematics robot. An algorithm determining the initial movement of the tooth in its elastic embedding controlled the set-up. The initial tooth movement was described by constant compliances. The constants were obtained prior to the experiment in a parameterised finite element (FE) study on the basis of a validated FE model of a human molar. The long-term tooth movement was assembled by adding up a multiple of incremental steps of initial tooth movements. A pure translational movement of the tooth of about 8 mm resulted for a moment to force ratio of ? 8.85 mm, corresponding to the distance between the bracket and the centre of resistance. The correct behaviour of this linear elastic model in its symmetry plane allows for simulating single tooth movement induced by orthodontic devices.     

面部软组织厚度测量及其在面貌复原中的应用

软组织厚度作为颅骨面貌复原的基础, 具有重要的应用价值。本文借助计算机技术对西安地区132例成年人颅面数据样本开展软组织测量、分析及应用研究, 结果表明, 1)通过分析特征点处软组织厚度和面部软组织分布图, 发现面部软组织分布具有一定的规律, 额头区域软组织厚度薄且样本间差异小, 脸颊区域软组织厚且样本间差异大; 2)通过比较不同年龄段男性软组织厚度的均值, 发现20-30岁阶段软组织厚度均值最小, 50-60岁阶段软组织厚度均值其次, 30-40岁阶段软组织厚度均值最大, 但30-40岁和40-50岁两个年龄段的软组织厚度近似; 通过比较不同年龄段女性软组织厚度的均值, 发现20-30岁阶段软组织厚度均值最小, 30-40岁阶段软组织厚度均值其次, 40-50岁阶段的软组织厚度均值最大; 3)特征点处软组织厚度标准差可以反映面貌体态的差异, 因此根据10个脸颊特征点的软组织厚度均值和标准差实现面貌体态分类; 4)根据不同性别、年龄、体态对应的软组织平均厚度, 应用计算机技术实现给定颅骨的三维面貌复原, 复原结果相比于传统手工复原的结果更加科学。     

Computational modelling and optimisation of soft tissue outcome in cranio-maxillofacial surgery planning

Cranio-maxillofacial (CMF) surgery operations are associated with rearrangement of facial hard and soft tissues, leading to dramatic changes in facial geometry. Often, correction of the aesthetical patient's appearance is the primary objective of the surgical intervention. Due to the complexity of the facial anatomy and the biomechanical behaviour of soft tissues, the result of the surgical impact cannot always be predicted on the basis of surgeon's intuition and experience alone. Computational modelling of soft tissue outcome using individual tomographic data and consistent numerical simulation of soft tissue mechanics can provide valuable information for surgeons during the planning stage. In this article, we present a general framework for computer-assisted planning of CMF surgery interventions that is based on the reconstruction of patient's anatomy from 3D computer tomography images and finite element analysis of soft tissue deformations. Examples from our clinical case studies that deal with the solution of direct and inverse surgical problems (i.e. soft tissue prediction, inverse implant shape design) demonstrate that the developed approach provides a useful tool for accurate prediction and optimisation of aesthetic surgery outcome.     

安氏Ⅰ类错合拔牙与非拔牙矫治对口唇形态影响的研究

目的：研究安氏Ⅰ类错合拔牙与非拔牙矫治对口唇形态的影响。方法：从直丝弓矫治的AngleⅠ类错合患者治疗前后的X线侧位片中随机选取拔除4个第一前磨牙患者15例（A组）,非拔牙矫治患者15例（B组）,经X线头影软组织测量分析比较矫治前后拔牙组与非拔牙组口唇形态的变化,对所得数据进行统计学处理。结果：拔牙矫治后上下唇的突度有明显改善,平均减少1.42和2.03mm;上下唇的长度也平均增加0.51和1.58mm;非拔牙矫治患者治疗后鼻唇角、下唇突度、上下唇长度均有增加,但矫治前后无统计学差异。结论：拔牙矫治有利于减小上下唇突度从而改善软组织侧貌。     

安氏I类错合拔牙与非拔牙矫治对口唇形态影响的研究

目的:研究安氏Ⅰ类错合拔牙与非拔牙矫治对口唇形态的影响.方法:从直丝弓矫治的Angle Ⅰ类错合患者治疗前后的X线侧位片中随机选取拔除4个第一前磨牙患者15例(A组),非拔牙矫治患者15例(B组),经X线头影软组织测量分析比较矫治前后拔牙组与非拔牙组口唇形态的变化,对所得数据进行统计学处理.结果:拔牙矫治后上下唇的突度有明显改善,平均减少1.42和2.03 mm;上下唇的长度也平均增加0.51和1.58 mm;非拔牙矫治患者治疗后鼻唇角、下唇突度、上下唇长度均有增加,但矫治前后无统计学差异.结论:拔牙矫治有利于减小上下唇突度从而改善软组织侧貌.

Abstract：

Objictive: To investigate the effect of Angle Class Ⅰ malocclusion after orthodontic treatment, with and without extractions on lip position changes. Methods: 30 patients with Angle Class Ⅰ malocclusion were chosen. 15 patients were treated by 4 first-premolars extraction (Group A) and 15 patients were treated without extraction (Group B). The soft tissue X-ray cephalometric of the patients were measured before and after the treatment and compared statistically. Results: After the extraction treatment, the upper and lower projecting lip reduced by 1.42 mm and 2.03, mmrespectively. The length of the upper and lower lips increased by 0.51mm and 1.58mm, respectively. For the group B, the nasolabial angle, the lower lip protrusion, the length of upper and lower lips had been increased, though there had no statistical significance before and after treatment. Conclusions: After extraction treatment the upper and lower projecting lips decreased. The patients with extractment treatment had the facial aesthelics.     

Simulation of bone remodelling in orthodontic treatment

Orthodontic treatments not only displace irregular teeth but also induce responses in surrounding bone tissues. Bone remodelling is regarded as the regulatory mechanism triggered by mechanical loading. This study was aimed at investigating the effect of orthodontic loading on both tooth movement and neighbouring bone density distribution. A set of computational algorithms incorporating both external and internal remodelling mechanisms was implemented into a patient-specific 3D finite element (FE) model to investigate and analyse orthodontic treatment under four typical modes of orthodontic loading. The consequence of orthodontic treatment was reproduced numerically by using this FE-based technique. The results indicated that the diverse modes of orthodontic loading would result in different magnitudes of tooth movement and particular morphology of bone density distribution. It is illuminated that the newly developed algorithms may replicate the clinical situation more closely compared with the previous proposed method.     

Comparison of methods to determine the centre of resistance of teeth

In orthodontic treatment, the locations of the centre of resistance (CR) of individual teeth and the applied load system are the major determinants for the type of tooth movement achieved. Currently, CR locations have only been specified for a relatively small number of tooth specimen for research purposes. Analysing cone beam computed tomography data samples from three upper central incisors, this study explores whether the effort to establish accurate CR estimates can be reduced by (i) morphing a pre-existing simplified finite element (FE) mesh to fit to the segmented 3D tooth-bone model, and (ii) individualizing a mean CR location according to a small parameter set characterising the morphology of the tooth and its embedding. The FE morphing approach and the semi-analytical approach led to CR estimates that differ in average only 0.04 and 0.12 mm respectively from those determined by very time-consuming individual FE modelling (standard method). Both approaches may help to estimate the movement of individual teeth during orthodontic treatment and, thus, increase the therapeutic efficacy.     

Orbital and maxillofacial computer aided surgery: patient-specific finite element models to predict surgical outcomes

This paper addresses an important issue raised for the clinical relevance of Computer-Assisted Surgical applications, namely the methodology used to automatically build patient-specific finite element (FE) models of anatomical structures. From this perspective, a method is proposed, based on a technique called the mesh-matching method, followed by a process that corrects mesh irregularities. The mesh-matching algorithm generates patient-specific volume meshes from an existing generic model. The mesh regularization process is based on the Jacobian matrix transform related to the FE reference element and the current element.This method for generating patient-specific FE models is first applied to computer-assisted maxillofacial surgery, and more precisely, to the FE elastic modelling of patient facial soft tissues. For each patient, the planned bone osteotomies (mandible, maxilla, chin) are used as boundary conditions to deform the FE face model, in order to predict the aesthetic outcome of the surgery. Seven FE patient-specific models were successfully generated by our method. For one patient, the prediction of the FE model is qualitatively compared with the patient's post-operative appearance, measured from a computer tomography scan. Then, our methodology is applied to computer-assisted orbital surgery. It is, therefore, evaluated for the generation of 11 patient-specific FE poroelastic models of the orbital soft tissues. These models are used to predict the consequences of the surgical decompression of the orbit. More precisely, an average law is extrapolated from the simulations carried out for each patient model. This law links the size of the osteotomy (i.e. the surgical gesture) and the backward displacement of the eyeball (the consequence of the surgical gesture).     

Computer-assisted development of epitheses after optical recording of facial defects]

A major drawback of conventional impression techniques used for customizing facial prostheses is the fact that pressure applied deforms soft tissue, thus reducing accuracy and causing patient discomfort. A possible solution is the use of optical 3-D coordinate measuring techniques, such as the fringe projection technique, which enables precise and contact-free recording of facial surfaces. The application of this method is demonstrated on a patient who lost his left eye and part of the jaw bone during surgery for cancer. 3-D CAD software that supports the construction of a facial prosthesis from the data obtained has been developed. For this purpose, spline functions are used to define border curves, and the intact half of the face is used as a model for the prosthetic surface. The resulting digital data are used to construct first a model made of synthetic resin, and then a final wax model with the aid of rapid prototyping techniques.     

Deformation and stress distribution of the human foot after plantar ligaments release: A cadaveric study and finite element analysis

The majority of foot deformities are related to arch collapse or instability, especially the longitudinal arch. Although the relationship between the plantar fascia and arch height has been previously investigated, the stress distribution remains unclear. The aim of this study was to explore the role of the plantar ligaments in foot arch biomechanics. We constructed a geometrical detailed three-dimensional (3-D) finite element (FE) model of the human foot and ankle from computer tomography images. The model comprised the majority of joints in the foot as well as bone segments, major ligaments, and plantar soft tissue. Release of the plantar fascia and other ligaments was simulated to evaluate the corresponding biomechanical effects on load distribution of the bony and ligamentous structures. These intrinsic ligaments of the foot arch were sectioned to simulate different pathologic situations of injury to the plantar ligaments, and to explore bone segment displacement and stress distribution. The validity of the 3-D FE model was verified by comparing results with experimentally measured data via the displacement and von Mise stress of each bone segment. Plantar fascia release decreased arch height, but did not cause total collapse of the foot arch. The longitudinal foot arch was lost when all the four major plantar ligaments were sectioned simultaneously. Plantar fascia release was compromised by increased strain applied to the plantar ligaments and intensified stress in the midfoot and metatarsal bones. Load redistribution among the centralized metatarsal bones and focal stress relief at the calcaneal insertion were predicted. The 3-D FE model indicated that plantar fascia release may provide relief of focal stress and associated heel pain. However, these operative procedures may pose a risk to arch stability and clinically may produce dorsolateral midfoot pain. The initial strategy for treating plantar fasciitis should be non-operative.     

Biologic basis for maxillary osteotomies

Adult Rhesus monkeys were used as experimental models to investigate revascularization and bone healing in different single-stage anterior, posterior and total maxillary osteotomy techniques. Microangiographic and histologic studies demonstrated that intraosseous and intrapulpal circulation to the mobilized maxillary segments were maintained by the experimental flap designs which maintained intact soft tissue; the fragments healed by osseous union within six weeks without immobilization of the mandible. The treatment of many severe dental-facial deformities is difficult and challenging. Functional and stable occlusions with facial balance and harmony have been attained in many adult patients by maxillary osteotomy techniques. The Rhesus monkey serves as an excellent experimental model to develop new biologically sound maxillary surgical orthodontic techniques.     

Three dimensional patient-specific collagen architecture modulates cartilage responses in the knee joint during gait

Site-specific variation of collagen fibril orientations can affect cartilage stresses in knee joints. However, this has not been confirmed by 3-D analyses. Therefore, we present a novel method for evaluation of the effect of patient-specific collagen architecture on time-dependent mechanical responses of knee joint cartilage during gait. 3-D finite element (FE) models of a human knee joint were created with the collagen architectures obtained from T₂ mapped MRI (patient-specific model) and from literature (literature model). The effect of accuracy of the implementation of collagen fibril architecture into the model was examined by using a submodel with denser FE mesh. Compared to the literature model, fibril strains and maximum principal stresses were reduced especially in the superficial/middle regions of medial tibial cartilage in the patient-specific model after the loading response of gait (up to ?413 and ?26%, respectively). Compared to the more coarsely meshed joint model, the patient-specific submodel demonstrated similar strain and stress distributions but increased values particularly in the superficial cartilage regions (especially stresses increased >60%). The results demonstrate that implementation of subject-specific collagen architecture of cartilage in 3-D modulates location- and time-dependent mechanical responses of human knee joint cartilage. Submodeling with more accurate implementation of collagen fibril architecture alters cartilage stresses particularly in the superficial/middle tissue.     

Applications of stem cells in orthodontics and dentofacial orthopedics:Current trends and future perspectives

A simple overview of daily orthodontic practice involves use of brackets, wires and elastomeric modules. However, investigating the underlying effect of orthodontic forces shows various molecular and cellular changes. Also, orthodontics is in close relation with dentofacial orthopedics which involves bone regeneration. In this review current and future applications of stem cells(SCs) in orthodontics and dentofacial orthopedics have been discussed. For craniofacial anomalies, SCs have been applied to regenerate hard tissue(such as treatment of alveolar cleft) and soft tissue(such as treatment of hemifacial macrosomia). Several attempts have been done to reconstruct impaired temporomandibular joint. Also, SCs with or without bone scaffolds and growth factors have been used to regenerate bone following distraction osteogenesis of mandibular bone or maxillary expansion. Current evidence shows that SCs also have potential to be used to regenerate infrabony alveolar defects and move the teeth into regenerated areas. Future application of SCs in orthodontics could involve accelerating tooth movement, regenerating resorbed roots and expanding tooth movement limitations. However, evidence supporting these roles is weak and further studies are required to evaluate the possibility of these ideas.     

Digital image correlation and finite element modelling as a method to determine mechanical properties of human soft tissue in vivo

The mechanical properties of human soft tissue are crucial for impact biomechanics, rehabilitation engineering, and surgical simulation. Validation of these constitutive models using human data remains challenging and often requires the use of non-invasive imaging and inverse finite element (FE) analysis. Post-processing data from imaging methods such as tagged magnetic resonance imaging (MRI) can be challenging. Digital image correlation (DIC), however, is a relatively straightforward imaging method. DIC has been used in the past to study the planar and superficial properties of soft tissue and excised soft tissue layers. However, DIC has not been used to non-invasive study of the bulk properties of human soft tissue in vivo. Thus, the goal of this study was to assess the use of DIC in combination with FE modelling to determine the bulk material properties of human soft tissue. Indentation experiments were performed on a silicone gel soft tissue phantom. A two camera DIC setup was then used to record the 3D surface deformation. The experiment was then simulated using a FE model. The gel was modelled as Neo-Hookean hyperelastic, and the material parameters were determined by minimising the error between the experimental and FE data. The iterative FE analysis determined material parameters (μ=1.80 kPa, K=2999 kPa) that were in close agreement with parameters derived independently from regression to uniaxial compression tests (μ=1.71 kPa, K=2857 kPa). Furthermore the FE model was capable of reproducing the experimental indentor force as well as the surface deformation found (R²=0.81). It was therefore concluded that a two camera DIC configuration combined with FE modelling can be used to determine the bulk mechanical properties of materials that can be represented using hyperelastic Neo-Hookean constitutive laws.     

Three-dimensional finite element simulations of the mechanical response of the fingertip to static and dynamic compressions

The analysis of the mechanics of the contact interactions of fingers/handle and the stress/strain distributions in the soft tissues in the fingertip is essential to optimize design of tools to reduce many occupation-related hand disorders. In the present study, a three-dimensional (3D) finite element (FE) model for the fingertip is proposed to simulate the nonlinear and time-dependent responses of a fingertip to static and dynamic loadings. The proposed FE model incorporates the essential anatomical structures of a finger: skin layers (outer and inner skins), subcutaneous tissue, bone and nail. The soft tissues (inner skin and subcutaneous tissue) are considered to be nonlinearly viscoelastic, while the hard tissues (outer skin, bone and nail) are considered to be linearly elastic. The proposed model has been used to simulate two loading scenarios: (a) the contact interactions between the fingertip and a flat surface and (b) the indentation of the fingerpad via a sharp wedge. For case (a), the predicted force/displacement relationships and time-dependent force responses are compared with the published experimental data; for case (b), the skin surface deflection profiles were predicted and compared with the published experimental observations. Furthermore, for both cases, the time-dependent stress/strain distributions within the tissues of the fingertip were calculated. The good agreement between the model predictions and the experimental observations indicates that the present model is capable of predicting realistic time-dependent force/displacement responses and stress/strain distributions in the soft tissues for dynamic loading conditions.