Individual-specific multi-scale finite element simulation of cortical bone of human proximal femur

• Published in: Journal of computational physics Vol. 244; pp. 298 - 311
• Main Authors: Ascenzi, Maria-Grazia, Kawas, Neal P., Lutz, Andre, Kardas, Dieter, Nackenhorst, Udo, Keyak, Joyce H.
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 01.07.2013
• Subjects: 60 APPLIED LIFE SCIENCES; APATITES; Bone; BONE FRACTURES; BONE MINERAL DENSITY; BONE TISSUES; Bones; CAT SCANNING; COLLAGEN; Computation; Computer programs; FEMUR; FINITE ELEMENT METHOD; Finite element methods; Mathematical analysis; Micro-structure; MICROSTRUCTURE; Multi-scale model; NMR IMAGING; Physiological loading; RADIOLOGY AND NUCLEAR MEDICINE; Remodeling; Secondary osteons; SIMULATION; Strain; STRAINS; Stresses; SURGERY; Multi-scale model; Secondary osteons; QCT; MRI; Micro-structure; Finite element methods; ROI; 2PM; 3D; CT-scan; BMD; Bone; Physiological loading; 2AL; 1 M; FE
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2012.05.027

We present an innovative method to perform multi-scale finite element analyses of the cortical component of the femur using the individual’s (1) computed tomography scan; and (2) a bone specimen obtained in conjunction with orthopedic surgery. The method enables study of micro-structural characteristics regulating strains and stresses under physiological loading conditions. The analysis of the micro-structural scenarios that cause variation of strain and stress is the first step in understanding the elevated strains and stresses in bone tissue, which are indicative of higher likelihood of micro-crack formation in bone, implicated in consequent remodeling or macroscopic bone fracture. Evidence that micro-structure varies with clinical history and contributes in significant, but poorly understood, ways to bone function, motivates the method’s development, as does need for software tools to investigate relationships between macroscopic loading and micro-structure. Three applications – varying region of interest, bone mineral density, and orientation of collagen type I, illustrate the method. We show, in comparison between physiological loading and simple compression of a patient’s femur, that strains computed at the multi-scale model’s micro-level: (i) differ; and (ii) depend on local collagen-apatite orientation and degree of calcification. Our findings confirm the strain concentration role of osteocyte lacunae, important for mechano-transduction. We hypothesize occurrence of micro-crack formation, leading either to remodeling or macroscopic fracture, when the computed strains exceed the elastic range observed in micro-structural testing.