CS-FEM在膝关节有限元模型中的应用

Application of CS-FEM to finite element modelling of knee joints

摘要

目的 通过使用光滑有限元法（cell-based smoothed finite element method， CS-FEM）来提高膝关节有限元模型中网格的质量与效率。 方法 在传统有限元法的基础上引入梯度光滑技术，采用基于单元的光滑有限元法（CS-FEM）对膝关节有限元模型进行网格划分，对比传统有限元法与CS-FEM实验数据的精确性。 结果 CS-FEM较传统有限元法更加精确且减少了网格划分的复杂度，同时降低了单元数量，使得求解过程更加高效。 结论 传统有限元法（finite element method， FEM）方法由于其应用广泛和成熟性，在处理相对简单的几何结构和线性材料时具有较高的计算精度和稳定性。 然而，在复杂结构和非线性材料的生物力学模拟中，CS-FEM通过优化网格划分和减少应力集中，提供了更高的计算效率和精度。 因此，在未来膝关节和其他复杂生物结构的建模中，CS-FEM有望得到更广泛的应用。

Abstract

Objective To improve the quality and efficiency of the mesh in the finite element model of the knee joint by using the cell-based smoothed finite element method (CS-FEM). Methods The gradient smoothing technique is introduced on the basis of the traditional finite element method, and CS-FEM is used to mesh the finite element model of the knee joint, and to compare which of the two methods (traditional finite element method and CS-FEM) is more accurate for the experimental data of the previous study. Results CS-FEM is more accurate and reduces the complexity of meshing than the traditional finite element method, and reduces the number of cells, making the solution process more efficient. Conclusion The conventional FEM methods offer high computational accuracy and stability when dealing with relatively simple geometries and linear materials due to their wide application and maturity. However, in biomechanical simulation of complex structures and nonlinear materials, CS-FEM provides higher computational efficiency and accuracy by optimizing meshing and reducing stress concentration. Therefore, CS-FEM is expected to be more widely used in the modelling of knee joints and other complex biological structures in the future.

关键词

膝关节 / 有限元模型 / CS-FEM网格划分技术

Key words

knee joint / finite element modelling / CS-FEM meshing technique

参考文献

[1] ESTELL E G, MURPHY L A, GANGI L R, et al. Attachment of cartilage wear particles to the synovium negatively impacts friction properties[J]. Journal of biomechanics, 2021, 127: 110668.

[2] 陈彦飞, 鲁超, 赵勇, 等. 基于CT影像动态膝关节有限元模型的构建及仿真力学分析[J]. 中国骨伤, 2020, 33(5): 479-484.

[3] CHEN Yanfei, LU Chao, ZHAO Yong, et al. Construction and simulation mechanical analysis of dynamic knee joint finite element model based on CT image[J]. China journal of orthopaedics and traumatology, 2020, 33(5): 479-484(in Chinese).

[4] RAO H, BAKKER R, MCLACHLIN S, et al. Computational study of extrinsic factors affecting ACL strain during single-leg jump landing[J]. BMC musculoskeletal disorders, 2024, 25(1): 318.

[5] 鲍春雨, 郭宝川, 孟庆华. 人体膝关节有限元模型建立及其有效性验证[J]. 应用力学学报, 2017, 34(3): 559-563.

[6] BAO Chunyu, GUO Baochuan, MENG Qinghua. A finite element model of human knee joint and its validation[J]. Chinese journal of applied mechanics, 2017, 34(3): 559-563(in Chinese).

[7] 姜亚琼, 路坦, 徐彪, 等. 内侧副韧带不同程度损伤时膝关节的应力变化[J]. 中国组织工程研究, 2024, 28(33): 5270-5275.

[8] JIANG Yaqiong, LU Tan, XU Biao, et al. Stress changes of knee joint with different degrees of medial collateral ligament injury[J]. Chinese journal of tissue engineering research, 2024, 28(33): 5270-5275(in Chinese).

[9] 张楠, 孟庆华, 鲍春雨, 等. 脑卒中患者运动过程中动力学特征的智能预测[J]. 医用生物力学, 2024, 39(3): 489-496.

[10] ZHANG Nan, MENG Qinghua, BAO Chunyu, et al. Intelligent prediction for dynamic characteristics of stroke patients during exercise[J]. Journal of medical biomechanics, 2024, 39(3): 489-496(in Chinese).

[11] RAHIMIZADEH A, NOURMOHAMMADI Z, ARABNEJAD S, et al. Porous architected biomaterial for a tibial-knee implant with minimum bone resorption and bone-implant interface micromotion[J]. Journal of the mechanical behavior of biomedical materials, 2018, 78: 465-479.

[12] ALI A A, HARRIS M D, SHALHOUB S, et al. Combined measurement and modeling of specimen-specific knee mechanics for healthy and ACL-deficient conditions[J]. Journal of biomechanics, 2017, 57: 117-124.

[13] FAN Y H, ZHOU G Y, ZHAN G T, et al. Comparative study on the mechanical behavior of the interface between natural cartilage and artificial cartilage[J]. Soft materials, 2021, 19(4): 400-419.

[14] LIU G R, NGUYEN T T, DAI K Y, et al. Theoretical aspects of the smoothed finite element method (SFEM)[J]. International journal for numerical methods in engineering, 2007, 71(8): 902-930.

[15] LIU G R. A G space theory and a weakened weak (W²) form for a unified formulation of compatible and incompatible methods: part I theory [J]. International journal for numerical methods in engineering, 2010, 81(9): 1093-1126.

[16] CHEN J S, WU C T, YOON S, et al. A stabilized conforming nodal integration for Galerkin mesh-free methods[J]. International journal for numerical methods in engineering, 2001, 50(2): 435-466.

[17] LIU G R. Meshfree methods: moving beyond the finite element method[M]. 2nd ed. Boca Raton: CRC Press, 2009.

[18] ALMEIDA E S, SPILKE R L. Mixed and penalty finite element models for the nonlinear behavior of biphasic soft tissues in finite deformation: part Ⅱ - nonlinear examples[J]. Computer Methods in Biomechanics and Biomedical Engineering, 1997, 1(2): 151-170.

[19] DAI K Y, LIU G R, NGUYEN T T. An n-sided polygonal smoothed finite element method (nSFEM) for solid mechanics[J]. Finite elements in analysis and design, 2007, 43(11/12): 847-860.

[20] LIU G R, NGUYEN-THOI T, NGUYEN-XUAN H, et al. A node-based smoothed finite element method (NS-FEM) for upper bound solutions to solid mechanics problems[J]. Computers & structures, 2009, 87(1/2): 14-26.

[21] LIU G R, NGUYEN-THOI T, LAMK Y. A novel alpha finite element method (αFEM) for exact solution to mechanics problems using triangular and tetrahedral elements[J]. Computer methods in applied mechanics and engineering, 2008, 197(45/46/47/48): 3883-3897.

[22] 叶楠, 程建宇, 冯世哲. 基于CS-FEM的某制导炮弹弹体动力特性分析[J]. 机械设计, 2024, 41(4): 14-19.

[23] YE Nan, CHENG Jianyu, FENG Shizhe. Analysis on dynamic characteristics of guided projectile body based on CS-FEM[J]. Journal of machine design, 2024, 41(4): 14-19(in Chinese).

[24] 任爽, 时会娟, 张家豪, 等. 前交叉韧带重建术后移植物应力的有限元分析[J]. 北京大学学报(医学版), 2021, 53(5): 865-870.

[25] REN Shuang, SHI Huijuan, ZHANG Jiahao, et al. Finite element analysis of the graft stresses after anterior cruciate ligament reconstruction[J]. Journal of Peking University (health sciences), 2021, 53(5): 865-870(in Chinese).

[26] 杨骏良, 路坦, 徐彪, 等. 前交叉韧带部分断裂对膝关节应力影响的三维有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1347-1353.

[27] YANG Junliang, LU Tan, XU Biao, et al. Three-dimensional finite element analysis of effects of partial anterior cruciate ligament rupture on knee joint stress[J]. Chinese journal of tissue engineering research, 2024, 28(9): 1347-1353(in Chinese).

[28] PEŃA E, CALVO B, MARTÍNEZ MA, et al. A three-dimensional finite element analysis of the combined behavior of ligaments and menisci in the healthy human knee joint[J]. Journal of biomechanics, 2006, 39(9): 1686-1701.

[29] 刘晓敏, 王学锋, 吕劲, 等. 膝关节后外侧复合体有限元模拟重建的研究[J]. 中国修复重建外科杂志, 2012, 26(11): 1310-1314.

[30] LIU Xiaomin, WANG Xuefeng, LÜ Jin, et al. Analog reconstruction of posterolateral complex by the finite element[J]. Chinese journal of reparative and reconstructive surgery, 2012, 26(11): 1310-1314(in Chinese).

[31] GABRIEL MT, WONG E K, WOO S L Y, et al. Distribution of in situ forces in the anterior cruciate ligament in response to rotary loads[J]. Journal of orthopaedic research, 2004, 22(1): 85-89.

[32] SONG Y H, DEBSKI R E, MUSAHL V, et al. A three-dimensional finite element model of the human anterior cruciate ligament: a computational analysis with experimental validation[J]. Journal of biomechanics, 2004, 37(3): 383-390.

[33] GUAN W Y, BHOWMICK S, GAO G J, et al. A phase-field modelling for 3D fracture in elasto-plastic solids based on the cell-based smoothed finite element method[J]. Engineering fracture mechanics, 2021, 254: 107920.

[34] LI E, CHEN J N, ZHANG Z P, et al. Smoothed finite element method for analysis of multi-layered systems-applications in biomaterials[J]. Computers & structures, 2016, 168: 16-29.

[35] LIU G R, TRUNG N T. Smoothed finite element methods[M]. Boca Raton: CRC Press, 2010.

[36] LIU G R, NOURBAKHSHNIA N, CHEN L, et al. A novel general formulation for singular stress field using the es-fem method for the analysis of mixed-mode cracks[J]. International journal of computational methods, 2010, 7(1): 191-214.

[37] HOFACKER M, MIEHE C. A phase field model of dynamic fracture: robust field updates for the analysis of complex crack patterns[J]. International journal for numerical methods in engineering, 2013, 93(3): 276-301.

[38] GABRIEL A A, LI D, CHIOCCHETTI S, et al. A unified first-order hyperbolic model for nonlinear dynamic rupture processes in diffuse fracture zones[J]. Philosophical transactions of the royal society a: mathematical, physical and engineering sciences, 2021, 379(2196): 20200130.

[39] 常先坤, 张建新. 不同鞋底厚度的运动鞋对高中生步态参数和爆发力的影响[J]. 应用力学学报, 2024, 41(3): 719-728.

[40] CHANG Xiankun, ZHANG Jianxin. Gait parameters and explosive power of sports shoes with different sole thickness of high school students[J]. Chinese journal of applied mechanics, 2024, 41(3): 719-728(in Chinese).

[41] RAINA A, MIEHE C. A phase-field model for fracture in biological tissues[J]. Biomechanics and modeling in mechanobiology, 2016, 15(3): 479-496.

[42] MIEHE C, HOFACKER M, SCHÄNZEL L M, et al. Phase field modeling of fracture in multi-physics problems. Part II. Coupled brittle-to-ductile failure criteria and crack propagation in thermo-elastic-plastic solids[J]. Computer methods in applied mechanics and engineering, 2015, 294: 486-522.

[43] ULMER H, HOFACKER M, MIEHE C. Phase field modeling of brittle and ductile fracture[J]. PAMM, 2013, 13(1): 533-536.

[44] ALESSI R, MARIGO J J, MAURINI C, et al. Coupling damage and plasticity for a phase-field regularisation of brittle, cohesive and ductile fracture: one-dimensional examples[J]. International journal of mechanical sciences, 2018, 149: 559-576.

[45] 杜娟, 夏晓宇, 孟琳, 等. 基于肌骨动力学与有限元耦合模型的脊柱生物力学研究应用进展[J]. 应用力学学报, 2024, 41(1): 233-240.

[46] DU Juan, XIA Xiaoyu, MENG Lin, et al. A review on applications of computational models combining finite element methods and musculoskeletal model[J]. Chinese journal of applied mechanics, 2024, 41(1): 233-240(in Chinese).

[47] DUDA F P, CIARBONETTI A, SÁNCHEZ P J, et al. A phase-field/gradient damage model for brittle fracture in elastic-plastic solids[J]. International journal of plasticity, 2015, 65: 269-296.

[48] BORDEN MJ, HUGHES T J R, LANDIS C M, et al. A phase-field formulation for fracture in ductile materials: Finite deformation balance law derivation, plastic degradation, and stress triaxiality effects[J]. Computer methods in applied mechanics and engineering, 2016, 312: 130-166.

[49] AMBATI M, GERASIMOV T, DE LORENZIS L. Phase-field modeling of ductile fracture[J]. Computational mechanics, 2015, 55(5): 1017-1040.

[50] LIU G W, LI Q B, MSEKH MA, et al. Abaqus implementation of monolithic and staggered schemes for quasi-static and dynamic fracture phase-field model[J]. Computational materials science, 2016, 121: 35-47.

基金资助

国家社科基金资助项目(23BTY121)