This article has Open Peer Review reports available.
Does facial soft tissue protect against zygomatic fractures? Results of a finite element analysis
© Huempfner-Hierl et al. 2015
Received: 17 February 2015
Accepted: 9 June 2015
Published: 16 June 2015
Zygomatic fractures form a major entity in craniomaxillofacial traumatology. Few studies have dealt with biomechanical basics and none with the role of the facial soft tissues. Therefore this study should investigate, whether facial soft tissue plays a protecting role in lateral midfacial trauma.
A head-to-head encounter was simulated by way of finite element analysis. In two scenarios this impact - with and without soft tissues - was investigated to demonstrate the potential protective effects. To achieve realism, a transient simulation was chosen, which considers temporal dynamics and realistic material parameters derived from CT grey values.
The simulation results presented a typical zygomatic fracture with all relevant fracture lines. Including soft tissues did not change the maximum bony stress pattern, but increased the time period from impact to maximal stresses by 1.3 msec.
Although this could have clinical implications, facial soft tissues may be disregarded in biomechanical simulations of the lateral midface, if only the bony structures are to be investigated. Soft tissue seems to act as a temporal buffer only.
Lateral midfacial zygomatic fractures are frequently encountered in craniomaxillofacial traumatology. Typical causes are assaults, traffic accidents, or sports incidents [1–3]. Here a frequent situation is a player versus player impact in team sports like association football or rugby. Depending on local cultural habits 13 to 30 % of all sport-sustained fractures in the head and neck area are located in the lateral midface [2–7]. Typical victims are males aged between 18 to the mid thirties. The causative blunt impact often results from a head-to-head encounter as two players try to hit the ball with their heads, one reaching the ball, the other one his opponent’s zygoma.
Concerning biomechanical studies about facial traumatology researchers will always be confronted with certain difficulties. Many experiments have been performed on cadavers. Evidently, only restricted conclusions can be made, as cadavers will have undergone postmortal alterations and, in most cases, will not have been of the typical age group of persons suffering from zygomatic fractures. Moreover cadaver specimen will be destroyed in these experiment so that they are not repeatable.
Attempts have been made with small and big animal models, but whereas the anatomy of a sheep tibia may be comparable to the human tibia in a certain extent , the human facial skull will not be really represented by any animal model.
Since about thirty years finite element analysis (FEA) has expanded from technical application into biomechanical and medical research. Finite element models (FE-models) have developed from rather simple models at the beginning to very sophisticated 3D-models with increasing computing capacity and improving methods of data acquisition [9–12].
The authors have shown that finite element analysis can reproduce a head collision leading to a typical fracture pattern in a previous study without the integration of midfacial soft tissue . Regarding further biomechanical literature on zygomatic trauma, only reports concentrating on the field of zygomatic fracture osteosynthesis or necessary impact forces have been published by now [13–15].
The question arises, whether this simplification is acceptable, and how simulation of biomechanical parameters of facial soft tissue and bone would alter fracture patterns and stress propagation in the simulation of zygomatic fractures.
To answer this question a biomechanical study based on finite element analysis was initiated to investigate the influence of facial soft tissue in protecting against zygomatic fracture. The null hypothesis was that the facial soft tissue envelope would protect the lateral midface and would change the fracture pattern in a typical head-to-head encounter.
Two scenarios of head-to-head impacts as forehead versus zygoma impacts were created in ANSYS Workbench (ANSYS Classic V12.0.1; ANSYS Inc. Canonsburg, PA, USA). The first consisted of finite element models of two skulls without any soft tissue whereas in the second scenario soft tissue parameters were included in the victim’s skull model. Besides presence of soft tissue, all other parameters were identical.
Soft tissue simulation
To achieve a realistic numerical calculation of stresses, especially at the bone surface, the Young’s modulus of the skull was modified by limiting the lowest value to 11 000 MPa according to experimental studies [11, 19]. Concerning boundary conditions, the nodes of the occipital condyles of the victims’ skull were fixed in all degrees of freedom. A coulomb friction model was created in which the solid body friction was divided into sticking and sliding friction and a coefficient of 0.4 was assigned . The assaulting head was not fixed, but velocity and vector were defined. For evaluation both scenarios were analyzed regarding the time-dependent propagation of stresses within soft tissue and bone. Moreover the final stress patterns were compared.
Von Mises stresses were evaluated for both scenarios. According to the studies of Nagasao et al.  a yield criterion of 153 MPa was defined, where material parameters change from elastic to plastic behaviour. Higher values will cause fractures represented by plastic material deformation. For the analysis of potential fractures a theory of mechanical engineering was utilized, which states that failure of the examined part will occur, if two stress gradients meet above the yield limit. This represents the point, at which the object will fail to resist the load and break [25–27].
According to the statutes of the local ethical review committee, no approval of this study had been necessary.
Finally, the maximum stress distribution was reached at 1.9 msec in the soft tissue set-up with an identical pattern as in the bone model (Fig. 8). In this specific set-up critical values above the yield criterion could be seen on the impact side. Lower, uncritical stresses were present on the opposite side in the region of the zygomatic buttress, the orbital floor, and lateral orbital rim. In both scenarios, a peak impact force of 45 N at the impact site could be determined.
So in both FEA-simulations – with and without bone covering facial soft tissue – the same maximum stress pattern had to be registered. In analysing time dependency a delay of 1.3 msec for reaching the full stress pattern had been found. Comparison of stress curves as a function of time are displayed in Fig. 8. Here it is easily discernable that in the model with simulation of soft tissue bone stresses will occur at a time, when they will already have finished in the bone only model. But the resulting maximum von Mises stresses are almost identical.
In a biomechanical investigation based on finite element analysis the role of the facial soft tissue in blunt zygomatic trauma was investigated. In two scenarios a head-to-head encounter like seen in sports’ accidents was modelled, one with and the other one without simulating the soft tissue coverage. The results showed that the midfacial soft tissue led to a delay of impact stresses of 1.3 msec compared to the pure bone scenario. Soft tissue did not change the peak stresses or the slope of stress build-up. So the postulated null hypothesis which supposed a protective effect of the soft tissue envelope in this trauma scenario was rejected.
The following aspects are of interest for the final analysis: the tool chain in itself, the computational results, and the influence of the modeling of soft tissue for the obtained results.
The tool chain permitted the creation of highly detailed bone volume meshes and allowed the use of individual material parameters for the investigated bone. As explained in previous studies [12, 21] and by Szwedowski et al.  using patient-specific Young’s moduli for bone derived from the individual CT-scan should increase the quality of the study and is a step towards the individualization of craniofacial biomechanical simulation.
However, regarding the soft tissue two simplifications were made: First, the scalp of the impacting head was not modeled because the scalp is a uniform thin layer that is hardly compressible in contrast to the three-to-sevenfold thicker and more pliable cheek. Our measurement of 3 to 4 mm scalp thickness derived from an automatic software-based calculation corresponds well to data taken from literature and measurements performed on the CT-scan [30, 31]. Additionally, facial soft tissues were modelled as one material and was not divided into skin, fat, and muscle. This seems reasonable according to the work of Zachow et al. , who reported on the use of uniform material parameters for maxillofacial surgery simulation. In literature, varying material parameters for Young’s moduli and Poisson ratio have been reported for soft tissue. In a survey Choi and Zheng  stated values for Poisson ratios from 0.3 to 0.5. Here the suggestions of Zachow et al.  were followed, who suggested a value between 0.43 and 0.45. Regarding material parameters soft tissue shows the peculiarity that they depend on the extent of deformation . In our study a Young’s modulus of 0.5 MPa was chosen, which lies in the range of the results of Choi and Zheng, who found 0.59 to 0.6 MPa, and Zachow et al. . Two further conditions of the scenario have to be mentioned. First, the impact was modelled as a head butt without full body impact like seen in American football, in which full body mass instead of the head mass should be taken for the encounter . Secondly, neck flexion in impactor and victim as well as hair were not considered. Modelling neck flexion is important, if neck strain is to be evaluated, but as this was not intended in this investigation, it was omitted.
The second point of discussion is the accordance of the simulated fracture pattern with clinical findings. Here the model simulates all relevant fracture lines ranging from the lateral orbital wall, infraorbital rim, zygomatic buttress, and zygomatic arch. Even potential comminution of the zygomatic buttress is resembled and equals typical clinical CT scans (Fig. 5). Thus it can be stated that our model created a clinically correct copy of a typical lateral midfacial fracture.
The last issue is the influence of facial soft tissue regarding a blunt impact, in this specific case a head-butt as might be encountered in a sports event. The resulting stress patterns were almost identical for the chosen set-up and equalled a typical zygomatic fracture. Thus, under our chosen model specifications, facial soft tissue would not change the resulting fracture type. Whether there could be changes regarding the amount of comminution cannot be judged from our results. The question remains why no decrease in peak stresses could by seen in the soft tissue scenario. A possible answer could lie in the fact that the cheek is pliable and shows a lateral displacement if hit by a blunt impactor. As there are no obstacles for this soft tissue shift, the impactor would hit the remaining thin layer of skin and remaining fat after 1.3 msec delay. It can be concluded by the results that this lateral soft tissue displacement would not decrease the impact force substantially as identical peak stress values and slope of stress build-up were encountered in the zygoma.
As demonstrated, soft tissue seems to act as a temporal buffer increasing time period to maximum stresses in bone from 0.6 to 1.9 msec. This delay of 1.3 msec could have clinical implications as soft tissue could decrease the impact and consecutive acceleration of the concomitant brain tissues. As investigations on brain trauma biomechanics focus on linear and angular acceleration and not on impact force the effect can not be calculated [34, 35]. Therefore the potential effect of reducing the probability of subdural hematoma via bridge vein rupture or focal brain concussion will not be discussed here. A further possible effect of impact delay could lie in allowing the victim to evade from the impactor or perform defensive moves.
On the other hand it is questionable, whether a delay of only 1.3 msec for reaching the full stress pattern in the soft tissue scenario in comparison to the bone only scenario will make any noteworthy difference concerning the clinical result.
It is difficult to discuss our findings with preexisting literature as this investigation is according to our knowledge the first study on the influence of soft tissue on craniomaxillofacial trauma patterns. The only result which may be sufficiently discussed is the question whether the resulting stress pattern resembles a typical zygomatic fracture. As mentioned above the pattern displayed in Fig. 4 is in accordance with the findings demonstrated in existing studies [36, 37].
Within the limitations of this study, it can be concluded that facial soft tissue has only little influence on maximum stresses in bone and fracture pattern in case of a blunt head-to-head encounter.
The simulation of a head impact to the lateral midface with simulation of soft tissue results in protracting the maximum stress pattern by 1.3 msec.
Whether this acts as buffer cannot really be answered. Our results suggest that facial soft tissue does not really play a decisive role for protecting the facial skull.
Furthermore they suggest that soft tissue simulation may be disregarded in the lateral midface in instances, in which only stress patterns would be examined as it complicates simulation. So the results concerning influence of soft tissue support and justify existing FEA-models, in which soft tissue is neglected.
Facial soft tissue simulation will be a necessary step in special clinical questions like the effects of protective devices in sports which will be addressed in future investigations. The results demonstrate that finite element simulation is an appropriate means to perform biomechanical investigations as they correlate well with clinical findings.
We acknowledge support from the German Research Foundation (DFG) and Leipzig University within the program of Open Access Publishing.
The development of FAT software was supported by German Federal Ministry of Economics and Technology (Grant ZIM KF2036708SSo).
- Kraft A, Abermann E, Stigler R, Zsifkovits C, Pedross F, Kloss F, et al. Craniomaxillofacial trauma: synopsis of 14,654 cases with 35,129 injuries in 15 years. Craniomaxillofac Trauma Reconstruction. 2012;5:41–50.View ArticleGoogle Scholar
- Elhammali N, Bremerich A, Rustemeyer J. Demographical and clinical aspects of sports-related maxillofacial and skull base fractures in hospitalized patients. Int J Oral Maxillofac Surg. 2010;39:857–62.View ArticlePubMedGoogle Scholar
- Roccia F, Diaspro A, Nasi A, Berrone S. Management of sport-related maxillofacial injuries. J Craniofac Surg. 2008;19:377–82.View ArticlePubMedGoogle Scholar
- Antoun J, Lee KH. Sports-related maxillofacial fractures over an 11-year period. J Oral Maxillofac Surg. 2008;66:504–8.View ArticlePubMedGoogle Scholar
- Carroll SM, Jawad MA, West M, O’Connor TPF. One hundred and ten sports related facial fractures. Br J Sports Med. 1995;29:194–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Delibasi C, Yamazawa M, Nomura K, Iida S, Kogo M. Maxillofacial fractures sustained during sports played with a ball. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2004;97:23–7.View ArticleGoogle Scholar
- Frenguelli A, Ruscito P, Bicciolo G, Rizzo S, Massarelli M. Head and Neck Trauma in Sporting Activities. J Craniomaxillofac Surg. 1991;19:178–81.View ArticlePubMedGoogle Scholar
- den Boer FC, Patka P, Bakker FC, Wippermann BW, van Lingen A, Vink GQ, et al. New Segmental long bone defect model in sheep: quantitative analysis of healing with dual energy x-ray absorptiometry. J Orthop Res. 1999;17:654–60.View ArticlePubMedGoogle Scholar
- Takizawa H, Suigiura K, Baba M, Tachisawa T, Kadoyama S, Kabayama T, et al. Structural Mechanics of the Blowout Fracture: Numerical Computer Simulation of Orbital Deformation by the Finite Element Method. Neurosurgery. 1988;22:1053–5.View ArticlePubMedGoogle Scholar
- Voo L, Kumaresan S, Pintar FA, Yoganandan N, Sances A. Finite-element models of the human head. Med Biol Eng Comput. 1996;34:375–81.View ArticlePubMedGoogle Scholar
- Nagasao T, Miyamoto J, Nagasao M, Ogata H, Kaneko T, Tamaki T, et al. The Effect of Striking Angle on the Buckling Mechanism of Blowout Fracture. Plast Reconstr Surg. 2006;117:2373–80.View ArticlePubMedGoogle Scholar
- Schaller A, Voigt C, Huempfner-Hierl H, Hemprich A, Hierl T. Transient finite element analysis of a traumatic fracture oft the zygomatic bone caused by a head collision. Int J Oral Maxillofac Surg. 2012;41:66–73.View ArticlePubMedGoogle Scholar
- Rinehart GC, Marsh JL, Hemmer KM, Bresina S. Internal fixation of malar fractures: an experimental biophysical study. Plast Reconstr Surg. 1989;84:21–5.View ArticlePubMedGoogle Scholar
- Garza JR, Baratta RV, Odinet K, Metzinger S, Bailey D, Best R, et al. Impact tolerances of the rigidly fixated maxillofacial skeleton. Ann Plast Surg. 1993;30:212–6.View ArticlePubMedGoogle Scholar
- Samolczyk-Wanyura D, Wroblewski G. Numerical study of the upper part of the face displacements using finite element method. Czas Stomatol. 2007;8:503–18.Google Scholar
- Taddei F, Pancanti A, Viceconti M. An improved method for the automatic mapping of computed tomography numbers onto finite element models. Med Eng Phys. 2004;26:61–9.View ArticlePubMedGoogle Scholar
- Morgan EF, Bayraktar HH, Keaveny TM. Trabecular bone modulus-density relationship depend on anatomic site. J Biomech. 2003;36:897–904.View ArticlePubMedGoogle Scholar
- Huiskes R. Finite element analysis of acetabular reconstruction. Acta Orthop. 1987;58:620–5.View ArticleGoogle Scholar
- Dechow PC, Nail GA, Schwartz-Dabney CL, Ashman RB. Elastic properties of human supraorbital and mandibular bone. Am J Phys Anthropol. 1993;90:291–306.View ArticlePubMedGoogle Scholar
- Arampatzis A, Brüggemann GP, Metzler V. The effect of speed on leg stiffness and joint kinetics in human running. J Biomech. 1999;32:1349–53.View ArticlePubMedGoogle Scholar
- Schaller A, Huempfner-Hierl H, Hemprich A, Hierl T. Biomechanical mechanisms of orbital wall fractures – a transient finite element analysis. J Craniomaxillofac Surg. 2013;41:710–7.View ArticlePubMedGoogle Scholar
- Zachow S, Hierl T, Erdmann B. Predictability of tissue changes after osteotomy planning in maxillofacial surgery: A comparison with postoperative results. Int Cong Ser. 2004;1268:648–53.View ArticleGoogle Scholar
- Hierl T, Hümpfner-Hierl H, Kruber D, Gäbler T, Hemprich A, Wollny G. Requirements for a universal image analysis tool in dentistry and oral and maxillofacial surgery. In: Daskalaki A, editor. Dental Computing and Applications. Advanced techniques for clinical dentistry. Hershey, London: IGI Global; 2009. p. 79–89.View ArticleGoogle Scholar
- Shockey JS, Von Frauenhofer JA, Seligson D. A measurment of the coefficient of static friction of human long bones. Surf Technol. 1985;25:167–73.View ArticleGoogle Scholar
- Becker A, Needleman A, Richmond O, Tvergaard V. Void growth and failure in notched bars. J Mech Phys Solids. 1988;36:317–51.View ArticleGoogle Scholar
- Needleman A, Tvergaard V. An analysis of ductile rupture in notched bars. J Mech Phys Solids. 1984;32:461–90.View ArticleGoogle Scholar
- Vecchio FJ. Nonlinear finite element analysis of reinforced concrete membranes. ACI Struct J. 1989;86:26–35.Google Scholar
- Pena Coto N, Cruz Meira JB, Brito e Dias R, Driemeier L, de Oliveira Roveri G, Noritomi PY. Assesment of nose protector for sport activities: finite element analysis. Dent Traumatol. 2012;28:108–13.View ArticleGoogle Scholar
- Szwedowski TD, Whyne CM, Fialkov JA. Towards characterization of craniofacial biomechanics. J Craniofac Surg. 2010;21:202–7.View ArticlePubMedGoogle Scholar
- Delye H, Verschueren P, Depreitere B, Van Lierde C, Verpoest I, Berckmans D, et al. Biomechanics of frontal skull fracture. In: Gilchrist MD, editor. IUTAM Symposium on Impact Biomechanics: From Fundamental Insights to Applications. Dordrecht: Springer; 2005. p. 185–93.View ArticleGoogle Scholar
- Doorly MC, Horgan TJ, Gilchrist MD. Reconstruction of head injury cases arising from falls using the UCD brain trauma model. In: Gilchrist MD, editor. IUTAM Symposium on Impact Biomechanics: From Fundamental Insights to Applications. New York Dordrecht: Springer; 2005. p. 431–40.View ArticleGoogle Scholar
- Choi P, Ling H, Zheng Y. Indentation technique for simultaneous estimation of Young’s modulus and Poisson’s ratio of soft tissue. In: Leondes CT, editor. Biomechanical Systems Technology. Muscular Skeletal Systems. Singapore: World Scientific Co; 2009. p. 227–44.View ArticleGoogle Scholar
- Viano DC. Head impact biomechanics in sport. In: Gilchrist MD, editor. IUTAM Symposium on Impact Biomechanics: From Fundamental Insights to Applications. Dordrecht: Springer; 2005. p. 121–30.View ArticleGoogle Scholar
- Schmitt K-U, Niederer PF, Cronin DS, Muser MH, Walz F. Trauma Biomechanics. An Introduction to Injury Biomechanics. Heidelberg New York: Springer; 2014. p. 55–80.Google Scholar
- Yoganandan N, Pintar FA, Zhang J, Gennarelli TA, Beuse N. Biomechanical Aspects of Blunt and Penetrating Head Injuries. In: Gilchrist MD, editor. IUTAM Symposium on Impact Biomechanics: From Fundamental Insights to Applications. Dordrecht: Springer; 2005. p. 173–84.View ArticleGoogle Scholar
- Lee EI, Mohan K, Koshy JC, Hollier LH. Optimizing the surgical management of zygomaticomaxillary complex fractures. Sem Plast Surg. 2010;4:389–97.View ArticleGoogle Scholar
- Ogata H, Sakamoto Y, Kishi K: A new classification of zygomatic fracture featuring zygomaticofrontal suture: injury mechanism and a guide to treatment. Plastic Surgery Int J 2013, ID 383486 DOI:10.5171/2013.383486Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.