the investigation of occlusal force integrates several domains of expertise and can facilitate an improved understanding of the mechanics of mastication , facial morphology , periodontal status , and temporomandibular joint diseases . especially during treatment planning for dental implants or during the maintenance phase for patients with periodontitis , evaluating occlusal force is a critical part of the clinical assessment . since loading is one of the factors that determine the outcome of dental implants and the maintenance of the periodontium , careful consideration of the biomechanics relating to loading is critical ; however , few , in any , studies have investigated the biomechanics of mastication . in the field of orthodontics anteroposterior disproportions have been categorized into class i ( normal ) , class ii ( retrusion of the mandible ) , and class iii ( protrusion of the mandible ) . a significantly altered pattern of occlusal force can significantly lower patients ' quality of life , especially that of periodontally compromised patients . a number of scientific studies have attempted to measure occlusal force in order to investigate the relationship between occlusal force and facial morphology . in order to study the effect of skeletal disproportions on masticatory force , casts and articulators have been routinely utilized to simulate masticatory movements . however , this methodology leads to inaccuracies because it is difficult to properly measure the exact occlusal force transferred from the teeth to the jaws . recently , finite element analysis ( fea ) has been introduced and used in various situations to assess the transfer of occlusal force , including onto implant prostheses and edentulous areas . , fea models were constructed using cone - beam computed tomography ( cbct ) images and dental scan images from three subjects with different molar occlusion relationships . the aim of this study was to evaluate the transfer of different occlusal forces in the three - dimensional ( 3d ) reconstructed images obtained from cbct and scanned dental models . three representative patients with classes i , ii , and iii malocclusions were selected from patients who visited the department of orthodontics at yonsei university . this study received the approval from the institutional review board , dental hospital , yonsei university ( reference number : 13 - 0070 ) . the skulls were scanned using a clinical cbct scanner ( rayscan symphony , ray co. , hwaseong , korea ) in the transverse plane with both a slice thickness and a scan increment of 2 mm , resulting in 20-slice images . the cbct scans were read into the visualization software ( dassault systems , waltham , ma , usa ) using a hp z800 workstation ( hewlett - packard co. , palo alto , ca , usa ) ; cpu , w5580 ( 3.2 ghz , turbo speed : 3.5 ghz , number of cores : 4 ) ; ram , 32 gb ; hdd , ssd , 128 gb ; vga , quadro fx4800 , 1.5 gb . bones and muscles were separated by uploading the cbct images into mimics ( materialise , leuven , belgium ) . two separate 3d files were exported to visualize the solid morphology of skeletal outlines without considering inner structures . this was performed to prevent calculation errors during the integration with the 3d scanned images of teeth . as well , the apical half of teeth during intercuspation was removed since it was not clearly distinguishable in the cbct images . the files created by the mimics software were transformed into the step file format using 3-matic ( materialise ) and finally analyzed with solidworks ( dassault systmes , velizy , france ) . the cast models were then scanned using a 3d optical laser scanner ( orapix kod-500 , seoul , korea ) , and imported into meshlab ( isti - cnr , pisa , italy ) . 2 ) . to simulate motion , two separate models were employed in the initial joint position . the initial stop points of teeth were designated as initial pressure points and colored in order to be utilized in structural simulations . for the structural simulation , connected joints were selected and joints were designated in the lower model to simulate initial pressure points onto colored areas . in this study , it was difficult to clearly visualize the muscles and ligaments , so the endpoints were intentionally defined at places where teeth were occluded . meshes were created with 677,852 total nodes and 474,152 total elements , each sized 1 mm with a tolerance of 0.05 mm ( fig . the following four weights were applied : 908 n ; 1,470 n ; 1,960 n ; and 2,205 n ( 100 kg , 150 kg , 200 kg , and 225 kg , respectively ) . 11.0 ( ansys inc . , canonsburg , pa , usa ) and processed on an hp xw6400 workstation ( hewlett - packard co. ) . the properties of the materials in the simulation were based on the properties of human bone determined when extracted teeth were physically fractured in a compression machine at 2,205 n. the specifications of the ceramic bone were as follows : tensile strength of 1.723410 n / m , compressive strength of 5.514910 n / m , young 's modulus of 2.205910 n / m , poisson 's ratio of 0.22 , mass density of 2,300 kg / m , and shear modulus of 9.040710 n / m . the 3d reconstruction of human skulls showed distinctive differences among the types of skeletal relationships ( fig . images from skulls with a class i relationship showed ideal anterior - posterior intermaxillary positions while images from skulls belonging to classes ii and iii showed skeletal malocclusions . complications presented by patients with each skeletal relationship were simulated . following the confirmation of mandible movement and projection , distortion measurements led to the following results : in class i , a 0.54-mm distortion was observed after loading 100 kg , a 0.80-mm distortion was observed after loading 150 kg , a 1.07-mm distortion was observed after loading 200 kg , and a 1.2-mm distortion was observed after loading 225 kg . in class ii , a 0.26-mm distortion was observed after loading 100 kg , a 0.38-mm distortion was observed after loading 150 kg , a 0.52 distortion was observed after loading 200 kg , and a 0.58-mm distortion was observed after loading 225 kg . in class iii , a 0.54-mm distortion was observed after loading 100 kg , a 0.80-mm distortion was observed after loading 150 kg , a 1.07-mm distortion was observed after loading 200 kg , and a 1.2-mm distortion was observed after loading 225 kg . the objective of this study was to evaluate the impact of different occlusal forces on the distortion of mandibles using fea of cbct images of the skeleton and 3d scanned dental images , which were later simulated by motion capture data . previous studies have shown that cbct images are appropriate for generating 3d models of teeth . however , the resolution is not always clear enough to distinguish all relevant details , although cbct images are widely used in various dental fields . the visualizations of dental surfaces and occlusion are not ideal because cbct images are subject to scattering from artifacts at the occlusal level . in the present study , dental information was manually acquired from impression models of each patient which were later scanned by a 3d scanner . the integration of digital dental casts into cbct scans has been attempted by various authors , and high - quality results have proven to be readily achievable . the technique of this integration may be of clinical relevance because it may provide important information for planning dental treatment and investigating oral disease . in particular , joint movement was simulated in the present study to fully mimic the movement of the temporomandibular joint and occlusion patterns . mechanically , bite force is substantianlly affected by the geometry of the lever system of the mandible , including the orientation of the muscles . since mechanical factors such as the ratio of anterior / posterior facial height , mandibular inclination , and the gonial angle are directly connected to the maximum bite force , more vertical and more acute gonial angles result in greater occlusal force . a number of studies have measured occlusal forces using strain gauges or miniature bite forks . however , these methods have certain limitations since they cause an open bite of at least 5 mm . therefore , indirect measurement methods have been introduced , involving the use sound transmission during mastication . recently , improved measurement methods have been introduced , involving pressure sensitive sheets ( dental prescale 50h type r , fuji film co. , tokyo , japan ) and charge - coupled device cameras . in the present study , since we used skeletal and dental information from korean patients , the results from this study can provide important information for further research . previous studies have shown that greater gonial angles provide an inferior level system for the elevator muscles and smaller perpendicular vectors to the occlusal table . as well , patients with no abnormalities experience greater occlusal forces than patients with dolichofacial abnormalities , whereas patients with brachyfacial abnormalities have even greater bite forces . this information is useful in the planning of implant treatment or rehabilitation after periodontal therapy . dental impressions are also simultaneously taken prior to the procedure . therefore , using such dental data , clinicians can predict the occlusal pattern , occlusal force , and the highly or minimally loaded areas using 3d simulation as proposed by the current authors in the present study . such procedures may become an essential part of total treatment planning in the near future . despite its limitations , the current study shows that the simulation of cbct images and dental scans can be successfully integrated and the occlusal forces can be clearly simulated in different skeletal malocclusions . further studies are warranted to fully evaluate the impact of skeletal malocclusion on masticatory performance using information about muscle attachment and 3d temporomandibular joint movement .
purposethe aim of this study was to evaluate the transfer of different occlusal forces in various skeletal malocclusions using finite element analysis ( fea).methodsthree representative human cone - beam computed tomography ( cbct ) images of three skeletal malocclusions were obtained from the department of orthodontics , yonsei university dental hospital , seoul , south korea . 
 the cbct scans were read into the visualization software after separating bones and muscles by uploading the cbct images into mimics ( materialise ) . 
 two separate three - dimensional ( 3d ) files were exported to visualize the solid morphology of skeletal outlines without considering the inner structures . 
 individual dental impressions were taken and stone models were scanned with a 3d scanner . 
 these images were integrated and occlusal motions were simulated . 
 displacement and von mises stress were measured at the nodes of the fea models . 
 the displacement and stress distribution were analyzed . 
 fea was performed to obtain the 3d deformation of the mandibles under loads of 100 , 150 , 200 , and 225 kg.resultsthe distortion in all three skeletal malocclusions was comparable . 
 greater forces resulted in observing more distortion in fea.conclusionsfurther studies are warranted to fully evaluate the impact of skeletal malocclusion on masticatory performance using information on muscle attachment and 3d temporomandibular joint movements.graphical abstract