Mimicking Bone Anisotropic Structure with Modified Gyroid Scaffolds; A Finite Element Analysis

Abstract

The structure of the bone is very complex and heterogeneous; this causes different mechanical and biological properties in its longitudinal and transverse directions. For example, the modulus of elasticity and the permeability of the trabecular bone in a longitudinal and radial direction can vary up to several times. Therefore, implant design that matches these differences is necessary to maximize compliance with the host bone. Given that, in this study, a gyroid structure that generally is used in bone scaffolds was modified to design anisotropic scaffolds. Therefore, the gyroid triply periodic minimal surface trigonometric function was manipulated, and five different architectures were denoted as G(-50), G(-25), G(0), G(+25), and G(+50) with a constant porosity of 80% were developed. The effective elastic moduli of the models were calculated using finite element analysis. The results showed an anisotropicity rate of 0.21, 0.62, 1.50 and 2.23 in elastic moduli for G(-50), G(-25), G(+25) and G(+50) models respectively. As well, the permeability of the models was calculated using computational fluid dynamics (CFD) analysis. Anisotropic models showed different permeability in longitudinal and transverse directions. Longitudinal permeability to lateral direction rate were 0.67, 0.80, 1.25 and 1.47 for G(-50), G(-25), G(+25) and G(+50) models respectively.Kemiğin yapısı karmaşık ve heterojendir, bu da boylamasına ve enine yönlerinde farklı mekanik ve biyolojik özelliklere neden olur. Örneğin, trabeküler kemiğin uzunlamasına ve enine yönde elastik modülü ve geçirgenliği birkaç kata kadar değişebilir. Dolaysıyla, implantların tasarımında konuk kemikle uyum sağlaması için bu farklılıkları dikkate almak gerekir. Bu çalışmada, yaygın olarak kemik iskeleleri tasarımında kullanılan gyroid yapısı, izotropik olmayan iskeleler modellemek için modifiye edilmiştir. Bu nedenle, gyroid üçlü periyodik minimal yüzey trigonometrik fonksiyonu manipüle edilerek ve %80 sabit bir gözenekliliğe sahip beş farklı iskele G (-50), G (-25), G (0), G (+25) ve G (+50) modeli elde edilmiştir. Modellerin etkili elastik modülleri sonlu elemanlar analizi kullanılarak hesaplanmıştır. Analiz sonuçları G (-50), G (-25), G (+25) ve G (+50) modellerin boylarınca elastisite modülünün enlerine göre sırasıyla 0.21, 0.62, 1.50 ve 2.23 oranda olduğunu göstermiştir. Ayrıca modellerin geçirgenliği hesaplamalı akışkanlar dinamiği (CFD) analizi kullanılarak hesaplanmıştır. İzotropik olmayan modeller boyuna ve enine yönlerde farklı geçirgenlik göstermiştir. G (-50), G (-25), G (+25) ve G (+50) modellerin geçirgenliği boylarınca enlerine göre oranı sırasıyla 0.67, 0.80, 1.25 ve 1.47 olarak hesaplanmıştır.