A systematic comparison between FEBio and PolyFEM for biomechanical systems

Publikation: Bidrag til tidsskriftTidsskriftartikelForskningfagfællebedømt

Background and Objectives: Finite element simulations are widely employed as a non-invasive and cost-effective approach for predicting outcomes in biomechanical simulations. However, traditional finite element software, primarily designed for engineering materials, often encountered limitations in contact detection and enforcement, leading to simulation failure when dealing with complex biomechanical configurations. Currently, a lot of model tuning is required to get physically accurate finite element simulations without failures. This adds significant human interaction to each iteration of a biomechanical model. This study addressed these issues by introducing PolyFEM, a novel finite element solver that guarantees inversion- and intersection-free solutions with completely automatic collision detection. The objective of this research is to validate PolyFEM's capabilities by comparing its results with those obtained from a well-established finite element solver, FEBio. Methods: To achieve this goal, five comparison scenarios were formulated to assess and validate PolyFEM's performance. The simulations were reproduced using both PolyFEM and FEBio, and the final results were compared. The five comparison scenarios included: (1) reproducing simulations from the FEBio test suite, consisting of static, dynamic, and contact-driven simulations; (2) replicating simulations from the verification paper published alongside the original release of FEBio; (3) a biomechanically based contact problem; (4) creating a custom simulation involving high-energy collisions between soft materials to highlight the difference in collision methods between the two solvers; and (5) performing biomechanical simulations of biting and quasi-stance. Results: We found that PolyFEM was capable of replicating all simulations previously conducted in FEBio. Particularly noteworthy is PolyFEM's superiority in high-energy contact simulations, where FEBio fell short, unable to complete over half of the simulations in Scenario 4. Although some of the simulations required significantly more simulation time in PolyFEM compared to FEBio, it is important to highlight that PolyFEM achieved these results without the need for any additional model tuning or contact declaration. Discussion: Despite being in the early stages of development, PolyFEM currently provides verified solutions for hyperelastic materials that are consistent with FEBio, both in previously published workflows and novel finite element scenarios. PolyFEM exhibited the ability to tackle challenging biomechanical problems where other solvers fell short, thus offering the potential to enhance the accuracy and realism of future finite element analyses.

OriginalsprogEngelsk
Artikelnummer107938
TidsskriftComputer Methods and Programs in Biomedicine
Vol/bind244
Antal sider14
ISSN0169-2607
DOI
StatusUdgivet - 2024

Bibliografisk note

Funding Information:
We acknowledge funding from NSF 2053851 , NSF 1835712 , NSF 1908767 , NIH R01HD097187 , NIH R01-HD083383-06 , NSERC DGECR-2021-00461 , RGPIN 2021-0370 , Horizon 2020 MSCA grant No. 764644 , and NIH R01DK133328 .

Publisher Copyright:
© 2023 Elsevier B.V.

ID: 383885994