收稿日期: 2024-10-02
网络出版日期: 2025-01-25
基金资助
国家自然科学基金(52111530189)
版权
Hydrodynamic finite element analysis of biological scaffolds with different pore sizes for cell growth and osteogenic differentiation
Received date: 2024-10-02
Online published: 2025-01-25
Supported by
the National Natural Science Foundation of China(52111530189)
Copyright
目的: 应用流体力学有限元分析方法,建立相同孔隙率、不同孔径的三周期极小曲面(triply periodic minimal surfaces,TPMS)螺旋形(Gyroid)多孔支架模型,模拟体内微环境,通过比较不同孔径支架的流体流速、壁面剪切应力、渗透率等相关参数变化,分析不同孔径的多孔支架对细胞黏附、增殖和成骨分化的可能影响。方法: 利用nTopology软件建立3组孔径的Gyroid多孔支架模型,支架尺寸为10 mm×10 mm×10 mm,孔径大小分别为400、600、800 μm,内部结构为各向同性。利用ANSYS 2022R1软件划分为非结构四面体网格,网格总数300+万。设定边界条件,流场域入口速度分别为0.01、0.1、1 mm/s,出口压力为0 Pa。根据Navier-Stokes方程计算流体流经支架时的压力、流速、壁面剪切应力,根据达西定律(Darcy’ s law)计算渗透率,使用ANSYS 2022R1软件中的Static structural模块对上述3种孔径支架的结构模型进行抗压强度分析。结果: 当入口流速分别为0.01、0.1、1 mm/s时,壁面剪切应力与流体流速呈线性关系,流速增加会导致壁面剪切应力增大。在0.1 mm/s流速下,当流体流经孔径为400、600、800 μm的3组支架时,压力呈梯度分布并逐渐减小,入口端压力依次为0.272、0.083、0.079 Pa;平均流速依次为0.093、0.078、0.070 mm/s;平均壁面剪切应力依次为2.955、1.343、1.706 mPa;渗透率依次为0.54×10-8、1.80×10-8、1.89×10-8 m2。计算3组支架内部最适合细胞黏附、增殖与成骨分化的壁面剪切应力范围所在区域占比,其中600 μm孔径支架该剪切应力范围内的内表面面积占比最大(27.65%),其次是800 μm孔径支架(17.30%),400 μm孔径支架占比最小(1.95%)。400、600和800 μm孔径支架的抗压强度依次为23、26、34 MPa。结论: 3组孔径的Gyroid支架在压应力作用下,应力分布均匀;600和800 μm孔径的Gyroid支架渗透率明显高于400 μm组,600 μm孔径的支架平均壁面剪切应力最小,且内部适合细胞生长与成骨分化的壁面区域占比最大,可能更适合于细胞黏附、增殖与成骨分化。
胡轶博 , 吕伟佳 , 夏炜 , 刘亦洪 . 基于细胞生长与成骨分化的不同孔径生物支架流体力学有限元分析[J]. 北京大学学报(医学版), 2025 , 57(1) : 97 -105 . DOI: 10.19723/j.issn.1671-167X.2025.01.015
Objective: The triply periodic minimal surface (TPMS) Gyroid porous scaffolds were built with identical porosity while varying pore sizes were used by fluid mechanics finite element analysis (FEA) to simulate the in vivo microenvironment. The effects of scaffolds with different pore sizes on cell adhesion, proliferation, and osteogenic differentiation were evaluated through calculating fluid velocity, wall shear stress, and permeability in the scaffolds. Methods: Three types of gyroid porous scaffolds, with pore sizes of 400, 600 and 800 μm, were established by nTopology software. Each scaffold had dimensions of 10 mm × 10 mm × 10 mm and isotropic internal structures. The models were imported to the ANSYS 2022R1 software, and meshed into over 3 million unstructured tetrahedral elements. Boun- dary conditions were set with inlet flow velocities of 0.01, 0.1, and 1 mm/s, and outlet pressure of 0 Pa. Pressure, velocity, and wall shear stress were calculated as fluid flowed through the scaffolds using the Navier-Stokes equations. At the same time, permeability was determined based on Darcy' s law. The compressive strength of scaffolds with different pore sizes was evaluated by ANSYS 2022R1 Static structural analysis. Results: A linear relationship was observed between the wall shear stress and fluid velocity at inlet flow rates of 0.01, 0.1 and 1 mm/s, with increasing velocity leading to higher wall shear stress. At the flow velocity of 0.1 mm/s, the initial pressures of scaffolds with pore sizes of 400, 600 and 800 μm were 0.272, 0.083 and 0.079 Pa, respectively. The fluid pressures were gradually decreased across the scaffolds. The average flow velocities were 0.093, 0.078 and 0.070 mm/s, the average wall shear stresses 2.955, 1.343 and 1.706 mPa, permeabilities values 0.54×10-8 1.80×10-8 and 1.89×10-8 m2 in the scaffolds with pore sizes of 400, 600 and 800 μm. The scaffold surface area proportions according with optimal wall shear stress range for cell growth and osteogenic differentiation were calcula-ted, which was highest in the 600 μm scaffold (27.65%), followed by the 800 μm scaffold (17.30%) and the 400 μm scaffold (1.95%). The compressive strengths of the scaffolds were 23, 26 and 34 MPa for the 400, 600 and 800 μm pore sizes. Conclusion: The uniform stress distributions appeared in all gyroid scaffold types under compressive stress. The permeabilities of scaffolds with pore sizes of 600 and 800 μm were significantly higher than the 400 μm. The average wall shear stress in the scaffold of 600 μm was the lowest, and the scaffold surface area proportion for cell growth and osteogenic differentiation the largest, indicating that it might be the most favorable design for supporting these cellular activities.
| 1 | Ray S , Nandi SK , Dasgupta S . Enhanced bone regeneration using Antheraea mylitta silk fibroin and chitosan based scaffold: In-vivo and in-vitro study[J]. Biomed Mater, 2023, 18 (5): 10. |
| 2 | Li Z , Tang S , Shi Z , et al. Multi-scale cellular PLA-based bionic scaffold to promote bone regrowth and repair[J]. Int J Biol Macromol, 2023, 245, 125511. |
| 3 | Zhang J , Tong D , Song H , et al. Osteoimmunity-regulating biomimetically hierarchical scaffold for augmented bone regeneration[J]. Adv Mater, 2022, 34 (36): e2202044. |
| 4 | 邹运, 韩青, 徐晓麟, 等. 骨科和口腔颌面外科3D打印模型的精度验证和可靠性分析[J]. 吉林大学学报(医学版), 2017, 43 (5): 996- 1001. |
| 5 | 吴其右, 崔博宇, 夏炜, 等. 基于细胞黏附的不同微结构3D打印多孔生物支架流体力学有限元分析[J]. 组织工程与重建外科杂志, 2024, 20 (3): 293- 299. |
| 6 | Luan HQ , Wang LT , Ren WY , et al. The effect of pore size and porosity of Ti6Al4V scaffolds on MC3T3-E1 cells and tissue in rabbits[J]. Sci China Technol Sci, 2019, 62 (7): 9. |
| 7 | Ma S , Tang Q , Han X , et al. Manufacturability, mechanical properties, mass-transport properties and biocompatibility of triply periodic minimal surface (TPMS) porous scaffolds fabricated by selective laser melting[J]. Mater Des, 2020, 195, 109034. |
| 8 | Wu J , Zhang Y , Lyu Y , et al. On the various numerical techniques for the optimization of bone scaffold[J]. Materials (Basel), 2023, 16 (3): 974. |
| 9 | Karageorgiou V , Kaplan D . Porosity of 3D biomaterial scaffolds and osteogenesis[J]. Biomaterials, 2005, 26 (27): 5474- 5491. |
| 10 | Ouyang P , Dong H , He X , et al. Hydromechanical mechanism behind the effect of pore size of porous titanium scaffolds on osteoblast response and bone ingrowth[J]. Mater Des, 2019, 183, 108151. |
| 11 | Tsuruga E , Takita H , Itoh H , et al. Pore size of porous hydroxyapatite as the cell-substratum controls BMP-induced osteogenesis[J]. J Biochem, 1997, 121 (2): 317- 324. |
| 12 | 王林, 马真胜, 李涤尘, 等. 灌注培养促进人胚成骨细胞在大体积可控微结构支架内的均匀扩增[J]. 中华医学杂志, 2013, 93 (25): 1970- 1974. |
| 13 | 崔越. 3D打印高强度三周期极小曲面羟基磷灰石支架用于骨修复的研究[D]. 广州: 华南理工大学, 2021. |
| 14 | Ali D , Ozalp M , Blanquer SBG , et al. Permeability and fluid flow-induced wall shear stress in bone scaffolds with TPMS and lattice architectures: A CFD analysis[J]. Euromech Fluids, 2020, 79, 376- 385. |
| 15 | Prakoso AT , Basri H , Adanta D , et al. The effect of tortuosity on permeability of porous scaffold[J]. Biomedicines, 2023, 11 (2): 427. |
| 16 | Porter B , Zauel R , Stockman H , et al. 3D computational mode-ling of media flow through scaffolds in a perfusion bioreactor[J]. J Biomech, 2005, 38 (3): 543- 549. |
| 17 | Pires T , Santos J , Ruben RB , et al. Numerical-experimental analysis of the permeability-porosity relationship in triply periodic minimal surfaces scaffolds[J]. J Biomech, 2021, 117, 110263. |
| 18 | 王真. 羟基磷灰石多孔骨支架的光固化制备工艺及力学与生物学性能研究[D]. 济南: 山东大学, 2020. |
| 19 | Zhu T , Cui Y , Zhang M , et al. Engineered three-dimensional scaffolds for enhanced bone regeneration in osteonecrosis[J]. Bioact Mater, 2020, 5 (3): 584- 601. |
| 20 | 姜至秀, 季俣辰, 刘丹瑜, 等. Gyroid结构钛仿生骨支架修复下颌骨节段性缺损的生物力学性能[J]. 中国组织工程研究, 2025, 29 (22): 4621- 4628. |
| 21 | Chan SW, Jusoh N, Abdul SA. Effect of fluid properties on bone scaffold permeability[C/OL]// 4th International Conference for Innovation in Biomedical Engineering and Life Sciences, 2022. (2024-03-22)[2024-06-26]. https://doi.org/10.1007/978-3-031-56438-3_3. |
| 22 | Prakoso AT , Basri H , Adanta D , et al. The effect of tortuosity on permeability of porous scaffold[J]. Biomedicines, 2023, 11 (2): 427. |
| 23 | 张传辉, 李建军, 杨军. 动态压力对负载胰岛素样生长因子1基因兔脂肪间充质干细胞增殖能力和机械性能的影响[J]. 中国组织工程研究, 2021, 25 (13): 6. |
| 24 | 熊婉琦, 李振豪, 崔焱, 等. 生物力学作用对成骨细胞生物特性的影响[J]. 中国组织工程研究, 2024, 28 (21): 3407- 3412. |
| 25 | Chen X , Guo J , Yuan Y , et al. Cyclic compression stimulates osteoblast differentiation via activation of the Wnt/β-catenin signaling pathway[J]. Mol Med Rep, 2017, 15 (5): 2890- 2896. |
| 26 | Yu W , Qu H , Hu G , et al. A microfluidic-based multi-shear device for investigating the effects of low fluid-induced stresses on osteoblasts[J]. PLoS One, 2014, 9 (2): e89966. |
| 27 | Pfister C , Bozsak C , Wolf P , et al. Cell shape-dependent shear stress on adherent cells in a micro-physiologic system as revealed by FEM[J]. Physiol Meas, 2015, 36 (5): 955- 966. |
/
| 〈 |
|
〉 |