Analysis of relationship between loading condition and cranial cracking pattern using a three-dimensional finite element model

Yoshimori Kiriyama; Yudai Sato; Yota Muramatsu; Teppei Mano; Katsumasa Tanaka; Kotaro Oshio

doi:10.1186/s12891-022-05215-x

Analysis of relationship between loading condition and cranial cracking pattern using a three-dimensional finite element model

BMC Musculoskelet Disord. 2022 Mar 31;23(1):310. doi: 10.1186/s12891-022-05215-x.

Authors

Yoshimori Kiriyama¹, Yudai Sato², Yota Muramatsu², Teppei Mano³, Katsumasa Tanaka⁴, Kotaro Oshio⁵

Affiliations

¹ Department of Mechanical Systems Engineering, Kogakuin University, Tokyo, Japan. kiriyama@cc.kogakuin.ac.jp.
² Department of Mechanical Systems Engineering, Kogakuin University, Tokyo, Japan.
³ Course of Mechanical Engineering, Kogakuin University, Tokyo, Japan.
⁴ Kogakuin University, Tokyo, Japan.
⁵ Department of Neurosurgery, St. Marianna University, Kanagawa, Japan.

Abstract

Background: A hairline crack on the cranium can occur even under a small external load or impact and are thus often observed in patients who have experienced an accidental fall or collision. Typical finite element analysis is useful to analyze the stress concentration or the propagation of stress waves. However, a stress propagation model does not accurately reproduce the features of hairline cracks on the cranium. The objective in this study was to reproduce cranial hairline cracks.

Methods: A three-dimensional finite element model of the cranial bone was developed from a patient CT images. The model consists of the frontal, parietal, occipital and temporal bones, and the bones are connected with the sutures. Additionally, the model comprised three layers; the external and internal tables and the diploe. The model was analyzed using the extended finite element method (X-FEM), and a forming limit diagram (FLD) was embedded in the model. In this study, the model was symmetrized bilaterally using the model developed from the left side of the skull. The FLD in this study was assumed to be a relationship between the maximum and minimum strains when a fracture occurs. A total of 13 typical loadings were applied to the model: loading points on the top, left, and back of the cranium were considered, and at each loading point, loads were applied with four or five different directions, namely perpendicular to the cranium and inclined in the anterior, posterior, superior, or inferior at an angle of 45^∘.

Results: Under all loading conditions, many small cracks formed radially at the loading points. Moreover, some large cracks formed under the certain loading conditions. The crack shapes on the top and left side could be associated with the specific loading directions, whereas cracks on the back did not show distinguishing characteristics depending on the loading directions. The present model was reproduced anatomically and morphologically, and the results were similar to those obtained in previous cadaver experiments.

Conclusions: Through X-FEM analysis of the FE model embedded with an FLD, hairline cracks in the cranium were reproduced, and a few crack shapes were identified as potential markers for estimating the loading conditions.

Keywords: Extended finite element method; Forming limit diagram; Hairline cracks; Skull fracture.

MeSH terms

Accidental Falls
Finite Element Analysis
Fractures, Bone*
Humans
Skull* / diagnostic imaging