Study on FE Modeling Method of the 50th Percentile Chinese Shank with Optimization Methods under Multi-Loading Condition
Zhang Guanjun1, Wang Longliang1, Deng Xiaopeng1, Du Xianping1, Guan Fengjiao2*
1(State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha 410082,China) 2(College of Mechatronic Engineering and Automation, National University of Defense Technology, Changsha 410073, China)
Abstract：The lower extremity of the human body is one of the most vulnerable parts in traffic accidents, and the finite element (FE) model of the lower limb has become an important tool to study mechanisms and protective methods of lower limb injuries. In order to ensure the fidelity of the model, it is necessary to carry out a comprehensive validation under multi-loading conditions. The biomechanical test data are different in the mechanical response due to the sample size and material diversity, so it is difficult to use one model to meet the multiple test data. The differences in geometric dimensions between Chinese and occidental human body would result in differences of biomechanical response, so it is necessary to develop the FE model of Chinese human body. In this paper, Geometry model based on the CT and MRI data was scaled to the 50 percentile male according to scaling factor derived from key dimensions of the tibia. The scaling method was also used to scale biomechanical test data to the biomechanical response data corresponding to the size of the leg. The elastic modulus, stress-strain curve, failure strain of the tibia and fibula and physical bulk modulus of the muscle were selected as design variables, the optimization method was used to fit tibia, fibula and calf biomechanical responses under different load positions in the quasi-static and dynamic loading conditions, which solved the problem with diversity caused by the test sample. When the elastic modulus of tibia and fibula were 18.43 and 18.23GPa, respectively, the failure strains were 1.156% and 0.8%, respectively, and physical bulk modulus of the calf was 11.33MPa, the calf model could fit the biomechanical responses well in many groups of tests. The FE model verification methods based on the optimization of multi-loading conditionsmade the FE model acquire better biofidelity effectively.
张冠军, 王龙亮, 邓小朋, 杜现平, 官凤娇. 中国50th人体小腿有限元模型多工况优化建模方法[J]. 中国生物医学工程学报, 2017, 36(2): 195-204.
Zhang Guanjun, Wang Longliang, Deng Xiaopeng, Du Xianping, Guan Fengjiao. Study on FE Modeling Method of the 50th Percentile Chinese Shank with Optimization Methods under Multi-Loading Condition. Chinese Journal of Biomedical Engineering, 2017, 36(2): 195-204.
[1] Zhang Guanjun, Cao Libo, Hu Jingwen, et al. A field data analysis of risk factors affecting the injury risks in vehicle-to-pedestrian crashes [J]. Ann Adv Automot Med, 2008, 52: 199-214. [2] Kim YS, Choi HH, Cho YN, et al. Numerical investigations of interactions between the knee-thigh-hip complex with vehicle interior structures [J]. Stapp Car Crash J, 2005, 49: 85-115. [3] Yang Jikuang. Injury biomechanics in car-pedestrian collisions: development, validation and application of human-body mathematical models [D]. Gothenburg: Chalmers University of Technology, 1997. [4] Bermond F, Ramet M, Bouquet R, et al. A finite element model of the pedestrian knee joint in lateral impact [C] // Proceedings of the International Research Council on the Biomechanics of Injury Conference. Eindhoven: INRETS-Bron, 1993: 117-129. [5] Wykowski E, Sinnhuber R, Appel H. Finite element model of human lower extremities in a frontal impact [C]//Proceedings of the International Research Council on the Biomechanics of Injury Conference. Goteborg: IRCOBI Secretariat, 1998: 101-116. [6] Untaroiu CD. Development and validation of a finite element model of human lower limb: including detailed geometry, physical material properties, and component validations for pedestrian injuries [D]. Charlottesville: University of Virginia, 2005. [7] Untaroiu CD, Yue N, Shin J. A finite element model of the lower limb for simulating automotive impacts [J]. Annals of Biomedical Engineering, 2013, 41(3): 513-526. [8] 杨济匡, 方海峰. 人体下肢有限元动力学分析模型的建立和验证 [J]. 湖南大学学报(自然科学版), 2005, 32(5): 31-36. [9] 张冠军. 行人下肢的碰撞损伤特性及相关参数研究 [D].长沙: 湖南大学, 2009. [10] 李正东, 刘宁国, 黄平, 等. 下肢有限元模型的建立及损伤机制重建 [J]. 中国司法鉴定, 2012(6): 37-42. [11] 蒋小晴, 杨济匡, 王丙雨, 等. 乘员股骨在轴向压力-弯矩下的损伤生物力学机理研究 [J]. 力学学报, 2013, 46(3): 465-474. [12] 刘俊先, 张兴和. 中国正常人体测量值 [M]. 北京: 中国医药科技出版社, 1994: 62-63. [13] Novitskaya E, Chen PY, Hamed E, et al. Recent advances on the measurement and calculation of the elastic moduli of cortical and trabecular bone: a review [J]. Theoretical and Applied Mechanics, 2011, 38(3): 209-297. [14] Beillas P, Begeman PC, Yang KH, et al. Lower limb: advanced FE model and new experimental data [J]. Stapp Car Crash J, 2001, 45: 469-494. [15] Karimi A, Navidbakhsh M. Material properties in unconfined compression of gelatin hydrogel for skin tissue engineering applications [J]. Biomedical Engineering/ Biomedizinische Technik, 2014, 59(6): 479-486. [16] Snedeker JG, Muser MH, Walz FH. Assessment of pelvis and upper leg injury risk in car-pedestrian collisions: comparison of accident statistics, impactor tests and a human body finite element model [J]. Stapp Car Crash J, 2003, 47: 437-458. [17] Iwamoto M, Tamura A, Furusu K, et al. Development of a finite element model of the human lower extremity for analyses of automotive crash injuries [C] //SAE 2000 World Congress. Detroit: SAE International, 2000: 2000-01-0621. [18] Schuster PJ, Chou CC, Prasad P, et al. Development and validation of a pedestrian lower limb non-linear 3-d finite element model [J]. Stapp Car Crash J, 2000, 44: 315-334. [19] Takahashi Y, Kikuchi Y, Konosu A, et al. Development and validation of the finite element model for the human lower limb of pedestrians [J]. Stapp Car Crash J, 2000, 44: 335-355. [20] Chawls A, Mukherjee S, Mohan D, et al. Validation of lower extremity model in THUMS [C] //Proceedings of the International Research Council on the Biomechanics of Injury conference. Graz: IRCOBI Secretariat, 2004: 155-166. [21] Arnoux PJ, Cesari D, Behr M, et al. Pedestrian lower limb injury criteria evaluation: a finite element approach [J]. Traffic Injury Prevention, 2005, 6(3): 288-297. [22] Untaroiu C, Darvish K, Crandall J, et al. A finite element model of the lower limb for simulating pedestrian impacts [J]. Stapp Car Crash J, 2005, 49: 157-181. [23] Silvestri C, Ray M. Development of a finite element model of the knee-thigh-hip of a 50th percentile male including ligaments and muscles [J]. International Journal of Crashworthiness, 2009, 14(2): 215-229. [24] Neale M, Thomas R, Bateman H, et al. A finite element modelling investigation of lower leg injuries [C] //Proceedings of the 20th International Technical Conference on the Enhanced Safety of Vehicles (ESV). Lyon: NHTSA, 2007: 07-0077. [25] Yamada H. Strength of biological materials [M]. Baltimore: The Williams & Wilkins Company, 1970: 49-73. [26] Ehler E, L sche H. Die menschliche tibia unter biegebelastung [J]. Beitr. Orthop, 1970, 17(5): 291-304. [27] Asang E. Experimental biomechanics of the human leg: a basis for interpreting typical skiing injury mechanisms [J]. The Orthopedic Clinics of North America, 1976, 7(1): 63-73. [28] Kerrigan J, Ivansson B, Bose D, et al. Rate-sensitive constitutive and failure properties of human collateral knee ligaments [C]//Proceedings of the International Research Council on the Biomechanics of Injury conference. Lisbon: IRCOBI Secretariat, 2003: 177-190. [29] Begeman P, Paravasthu N. Static and dynamic compression loading of the lower leg [C]//Proceedings of the Seventh Injury Prevention Through Biomechanics Symposium. Detroit: Wayne State University, 1997. [30] Burstein AH, Reilly DT, Martens M. Aging of bone tissue: mechanical properties [J]. J Bone Joint Surg Am, 1976, 58(1): 82-86. [31] Martens M, van Audekercke R, Delport P, et al. The mechanical characteristics of cancellous bone at the upper femoral region [J]. Journal of Biomechanics, 1983, 16(12): 971-983. [32] Galbusera F, Freutel M, Dürselen L, et al. Material models and properties in the finite element analysis of knee ligaments: a literature review[J]. Frontiers in Bioengineering and Biotechnology, 2014, 2: 1-11. [33] 官凤娇. 冲击载荷下的生物组织材料参数反求及损伤研究 [D]. 长沙: 湖南大学, 2011. [34] Lessley D, Crandall J, Shaw G, et al. A normalization technique for developing corridors from individual subject responses [C] //SAE 2000 World Congress. Detroit: SAE International, 2004: 2004-01-0288. [35] Takahashi Y, Kikuchi Y, Mori F, et al. Advanced finite element lower limb model for pedestrians [C] //Proceedings of the 18th International Technical Conference on the Enhanced Safety of Vehicles (ESV). Lyon: NHTSA, 2003: No 218. [36] Kerrigan JR. A computationally efficientmathematicalmodelofthepedestrianlowerextremity [D]. Charlottesville: University of Virginia, 2008. [37] Rho JY, Ashman RB, Turner CH. Young's modulus of trabecular and cortical bone material: ultrasonic and microtensile measurements [J]. J Biomech, 1993, 26(2): 111-119. [38] Snyder SM, Schneider E. Estimation of mechanical properties of cortical bone by computed tomography [J].J Orthop Res, 1991, 9(3): 422-431. [39] Ruan J, El-Jawahri R, Chai L, et al. Prediction and analysis of human thoracic impact response and injuries in cadaver impacts using a full human body finite element model [J]. Stapp Car Crash J, 2003, 47: 299-321.
