Bone Stress-Strain State Evaluation Using CT Based FEM

Frontiers in Mechanical Engineering

卷:7

Artur O. Fedyanin¹  Tatyana V. Baltina¹  Maxim E. Baltin¹  Filip O. Fadeev²  Rustem R. Islamov²  Oskar A. Sachenkov³  Oleg V. Gerasimov³  Nikita V. Kharin³  Pavel V. Bolshakov³  Evgeny O. Statsenko⁴ 

[1] Department of Human and Animal Physiology, Institute of Fundamental Medicine and Biology, Kazan Federal University, Kazan, Russia;

[2] Department of Medical Biology and Genetics, Kazan State Medical University, Kazan, Russia;

[3] Department of Theoretical Mechanics, Institute of Mathematics and Mechanics, Kazan Federal University, Kazan, Russia;

[4] Laboratory of X-ray Tomography, Institute of Geology and Petroleum Technologies, Kazan Federal University, Kazan, Russia;

关键词: bone imaging;    bone health;    fracture risk;    FEM;    CT;    CT/FEA;   

DOI  :  10.3389/fmech.2021.688474

来源: DOAJ

【 摘 要 】

Nowadays, the use of a digital prototype in numerical modeling is one of the main approaches to calculating the elements of an inhomogeneous structure under the influence of external forces. The article considers a finite element analysis method based on computed tomography data. The calculations used a three-dimensional isoparametric finite element of a continuous medium developed by the authors with a linear approximation, based on weighted integration of the local stiffness matrix. The purpose of this study is to describe a general algorithm for constructing a numerical model that allows static calculation of objects with a porous structure according to its computed tomography data. Numerical modeling was carried out using kinematic boundary conditions. To evaluate the results obtained, computational and postprocessor grids were introduced. The qualitative assessment of the modeling data was based on the normalized error. Three-point bending of bone specimens of the pig forelimbs was considered as a model problem. The numerical simulation results were compared with the data obtained from a physical experiment. The relative error ranged from 3 to 15%, and the crack location, determined by the physical experiment, corresponded to the area where the ultimate strength values were exceeded, determined by numerical modeling. The results obtained reflect not only the effectiveness of the proposed approach, but also the agreement with experimental data. This method turned out to be relatively non-resource-intensive and time-efficient.

【 授权许可】

Unknown