tensile and compressive forces of cervical tooth area lead to bonds interrupting among hydroxylapatite crystals of enamel and dentin , allowing penetration of small molecules , which are preventing re - establishment of bonds . abfraction lesions are characterized by a wedge - shaped defect of tissue , bordered with sharp internal and external edges [ figure 1 ] . abfraction lesions   each tooth tissue has different mechanical properties and morphology , which makes analysis of distribution process and stress concentration alongside these structures very complex . describe stress distribution alongside the enamel  as a cap around the tooth ,  which is finally concentrated at the cement - enamel junction . the aim of this study was to explore stress distribution inside tooth structures and abfraction lesions with different morphologies under various loading conditions . mandibular first premolar was scanned by micro computed tomography ( ct ) scanner ( skyscan 1076 kontich , belgium ) and the obtained images were reconstructed into transaxial sections using nrecon and ctan program software ( skyscan ) . approximately , 570 horizontal sections ( resolution up to 758  758 [ px ] ) were chosen for the reconstruction of three - dimensional ( 3d ) models [ figure 2 ] . micro computed tomography images of the mandibular first premolar a volumetric 3d cad tooth model was created using program packages matlab ( mathworks , inc . , the resulting model consisted of enamel , dentin , pulp , periodontal ligament , and reconstructed segment of alveolar bone [ figure 3 ] . in comparison to adjacent structures , young 's modulus of pulp is negligibly small , therefore pulp was modeled as an empty space . three - dimensional volumetric tooth tissue models reconstructed from micro computed tomography images ; enamel , dentin , pulp , periodontal ligament , and reconstructed segment of the alveolar bone to compare the stress distribution in lesions with different geometry , two types of cervical lesions were modeled ; the wedge - shaped lesion v lesion [ figure 4 ] and saucer - shaped lesion  u lesion [ figure 5 ] . v - shaped lesion , appearance on the three - dimensional tooth model , and lesion dimensions u - shaped lesion , appearance on the three - dimensional tooth model , and lesion dimensions the models were exported to finite element analysis software ansys workbench ( 14.0 ) and a finite element mesh of the models is made [ figure 6 ] . three - dimensional meshed models tooth tissue properties since the intensity of masticatory force and teeth contact surface is extremely variable , two occlusal contact types were chosen to show the effects of a favorable and unfavorable situation for the tooth and surrounding tissues . for the purpose of presenting the action of functional occlusal load , occlusal contacts are simulated on regions that imitate teeth contacts in central occlusion [ figure 7a ] . in addition to this , nonfunctional load is simulated on the external buccal cusp side , at 40 angle [ figure 7b ] . type of simulated occlusal contacts in ( a ) central occlusion ( b ) laterotrusion following steps were processing ( finite elements method [ fem ] calculation of tooth tissue straining ) and postprocessing ( result analysis by stress distribution criteria ) . to record the complex stress measured in ( mpa ) , von mises stress ( vms ) - the hypothesis of the highest distortion energy was used . vms determines the total of resulting stress in every point of the observed object  tooth element . results are presented in the form of images with the stress distribution presented as a color scale and numeric values obtained in ( mpa ) . the calculated vms on the sound tooth in the cervical region under functional load measured up to 12 ( mpa ) [ figure 8a ] , while the stresses in the same area under eccentric load are significantly higher and were over 50 ( mpa ) [ figure 8b ] . distribution of von mises stress at the cervical part of the tooth under ( a ) functional and ( b ) nonfunctional load of 200 ( n ) on the sagittal section of the cervical tooth area , the highest stress was observed along the dentino - enamel junction ( dej ) . the values of calculated vms in the region of sub - superficial enamel under nonfunctional load of 200 ( n ) range up to 60 ( mpa ) , which is 5 times higher in comparison to the values measured in superficial enamel [ figure 9 ] . distribution of von mises stress on the sagittal section of the superficial and sub - superficial enamel of sound tooth model under nonfunctional load of 200 ( n ) the sound model showed higher stress distribution pattern for both loading types compared with models with lesions . on the model with v - shaped lesions , discontinuity of the enamel and dentin resulted in a strain increase , as well as stress concentration around the apex / bottom of the lesion . values of the vms on the apex / bottom of the wedge lesion are extremely high and measured up to 266.13 ( mpa ) under nonfunctional load [ figures 10b and 11a ] and 93.24 ( mpa ) under functional load [ figures 10a ] . distribution of von mises stress in the v - shaped lesions under ( a ) functional and ( b ) nonfunctional load of 200 ( n ) distribution of von mises stress in the u - shaped lesions under ( a ) functional and ( b ) nonfunctional load of 200 ( n ) values of vms on the bottom of u - shaped lesion were up to 16 ( mpa ) under functional load [ figure 12a ] and 55 ( mpa ) under nonfunctional load [ figures 11b and 12b ] . higher values of strain up to 55 mpa under functional load and 180 ( mpa ) under nonfunctional load were measured with the saucer lesion only on the mesial and distal side of the lesion , where the enamel inserts with its sharp edges . distribution of von mises stress on the sagittal section of ( a ) v - shaped lesion and ( b ) u - shaped lesion and surrounding tooth tissues under nonfunctional load of 200 ( n ) the biomechanics of complex systems , such as a tooth , requires application of complex research methods . the fem used in this research is an adequate method for the biomechanical behavior prediction of the tooth under load . the models in this research are complex and accurate , not only because of the 3d structure , but also due to the large number of elements and nodes of which they are composed , as well as the fact that they are gained from a ct image that allowed a reproduction of the smallest morphological details [ figure 2 ] . the research tested the action of two types of occlusal loading on a model of sound tooth , model with wedge - shaped cervical lesion , i.e. , v lesion [ figure 4 ] and model with saucer - shaped cervical lesion , i.e. , u lesion [ figure 5 ] . this research showed that under eccentric forces , stress values on tested models were higher in all tooth tissues . occlusal stress is commonly reflected in the cervical tooth region , which is especially noticeable with the action of nonfunctional occlusal forces . calculated vms in the cervical tooth area under axial load measured up to 12 ( mpa ) , [ figure 8a ] , whereas the stresses in the same part of the tooth under paraxial load are significantly higher and measured up to 60 ( mpa ) [ figure 8b ] . yaman et al . , and rees and hammadeh reported that tensile stress values on the buccal aspect of the cervical region range from 60 to 90 ( mpa ) . rees and hammadeh analyzed the mechanisms of forming abfraction lesions and presented a theory of a possible undermining of the enamel due to stress concentration on the dej . results of our research confirmed the aforementioned theory , since the values of calculated stress in the cervical part of the tooth were 5 times higher in the sub - superficial compared to the superficial region of the enamel [ figure 9 ] . this indicates that breaking of the bonds between enamel prisms might occur in these layers exactly . the results of this study demonstrated that the geometric shape of the lesions is very important in the distribution of internal stress in the tooth . on the model with v - shaped lesion , the enamel and dentin discontinuity resulted in strain increase , as well as in stress concentration around the apex of the lesion . stress values on the apex / bottom of the v lesion are extremely high and measures up to 266.13 ( mpa ) under paraxial load [ figures 10b and 12a ] and 93.24 [ mpa ] under axial load [ figures 10a and 12b ] . the results of analyses showed that the bottom of cervical lesions concentrated stress under load , which clearly indicates that their further exposure to stress would lead to its deepening . stress values on the bottom of the u - shaped lesion under nonfunctional load measures up to 55 ( mpa ) [ figures 11b and 12b ] , which is 5 times lower comparing with v - shaped lesions . higher values of strain , up to 180 ( mpa ) , were calculated with the saucer lesion only in the areas where the enamel inserts with its sharp edges . according to these results , it can be concluded that the geometric shape of the lesion plays an important role in the distribution of the internal stress of the tooth . lesions with emphasized geometric discontinuity ( v - shaped lesions ) result in high strain concentration , whereas lesions with rounded contours ( u - shaped lesions ) concentrate less strain . the results obtained in this study can be important in clinical practice because they open a possibility of considering lesions remodeling with extremely sharp geometry into rounded lesions during their restoration . eccentric occlusal forces in all tooth tissues caused higher stress values on the tested modelsoccurrence of significant stress in the cervical part of the sound tooth model is caused by eccentric occlusal load . stress is almost 5 times higher in the sub - superficial layer of the cervical enamel in comparison to the superficial enamelgeometric shape of the lesion is very important in the distribution of internal stress in hard dental tissues . compared to the u - shaped lesions , v - shaped lesions show significantly higher stress concentrations under load . therefore , it could be concluded that their further exposure to stress would lead to its progressionthis research provides a good insight into stress distribution in one moment of tooth loading . limitation of this study is that the acting force is assumed as a static , not as a dynamic process , what it is in reality . eccentric occlusal forces in all tooth tissues caused higher stress values on the tested models occurrence of significant stress in the cervical part of the sound tooth model is caused by eccentric occlusal load . stress is almost 5 times higher in the sub - superficial layer of the cervical enamel in comparison to the superficial enamel geometric shape of the lesion is very important in the distribution of internal stress in hard dental tissues . compared to the u - shaped lesions , v - shaped lesions show significantly higher stress concentrations under load . therefore , it could be concluded that their further exposure to stress would lead to its progression this research provides a good insight into stress distribution in one moment of tooth loading . limitation of this study is that the acting force is assumed as a static , not as a dynamic process , what it is in reality .
objective : the present study aims to investigate the influence of presence and shape of cervical lesions on biomechanical behavior of mandibular first premolar , subjected to two types of occlusal loading using three - dimensional ( 3d ) finite element method ( fem).materials and methods:3d models of the mandibular premolar are created from a micro computed tomography x - ray image : model of sound mandibular premolar , model with the wedge - shaped cervical lesion ( v lesion ) , and model with saucer - shaped cervical lesion ( u lesion ) . by fem , 
 straining of the tooth tissues under functional and nonfunctional occlusal loading of 200 ( n ) is analyzed . 
 for the analysis , the following software was used : ctan program 1.10 and ansys workbench ( version 14.0 ) . 
 the results are presented in von mises stress.results:values of calculated stress in all tooth structures are higher under nonfunctional occlusal loading , while the functional loading is resulted in homogeneous stress distribution . 
 nonfunctional load in the cervical area of sound tooth model as well as in the sub - superficial layer of the enamel resulted with a significant stress ( over 50 [ mpa ] ) . 
 the highest stress concentration on models with lesions is noticed on the apex of the v - shaped lesion , while stress in saucer u lesion is significantly lower and distributed over wider area.conclusion:the type of the occlusal teeth loading has the biggest influence on cervical stress intensity . 
 geometric shape of the existing lesion is very important in the distribution of internal stress . compared to the u - shaped lesions , 
 v - shaped lesions show significantly higher stress concentrations under load . 
 exposure to stress would lead to its progression .