The Use of Finite Element Method Analysis for Modeling Different Osteotomy Patterns and Biomechanical Analysis of Craniosynostosis Correction

Purpose: Several post-processing algorithms for 3D visualization of the skull in craniosynostosis with their specific advantages and disadvantages have been already described. The Finite Element Method (FEM) described herein can also be used to evaluate the efficacy of the cutting patterns with respect to an increase in the projected surface area under assumed uniform loading of the manipulated and cut bone segments. Methods: The FEM analysis was performed. Starting with the classic cranial osteotomies for bifrontal craniotomy and orbital bandeau a virtually mirroring of the unaffected triangular shaped frontal bone was performed to achieve a cup-shaped sphere of constant thickness of 2.5 mm with a radius of 65 mm. Mechanical properties required for the analysis were Young's modulus of 340 MPa and Poisson's ratio of 0.22. Four different cutting patterns from straight to curved geometries have been projected onto the inner surface of the sphere with a cutting depth set to 2/3rds of the shell thickness. The necessary force for the deformation, the resulting tensions and the volume loss due to the osteotomy pattern were measured. Results: Better outcomes were realized with pattern D. The necessary force was 73.6% smaller than the control group with 66N. Best stress distribution was achieved. Curved cutting patterns led to the highest peak of stress and thus to a higher risk of fracture. Straight bone cuts parallel to the corners or to the thighs of the sphere provided a better distribution of stresses with a small area with high stress. Additionally, also with pattern D a surface increase of 20.7% higher than reference was registered. Conclusion: As a proof of concept for different cutting geometries for skull molding in the correction of craniosynostosis, this computational model shows that depending of the cutting pattern different biomechanical behavior is achieved.