1
|
Müller P, Synek A, Stauß T, Steinnagel C, Ehlers T, Gembarski PC, Pahr D, Lachmayer R. Development of a density-based topology optimization of homogenized lattice structures for individualized hip endoprostheses and validation using micro-FE. Sci Rep 2024; 14:5719. [PMID: 38459092 PMCID: PMC10923877 DOI: 10.1038/s41598-024-56327-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2023] [Accepted: 03/05/2024] [Indexed: 03/10/2024] Open
Abstract
Prosthetic implants, particularly hip endoprostheses, often lead to stress shielding because of a mismatch in compliance between the bone and the implant material, adversely affecting the implant's longevity and effectiveness. Therefore, this work aimed to demonstrate a computationally efficient method for density-based topology optimization of homogenized lattice structures in a patient-specific hip endoprosthesis. Thus, the root mean square error (RMSE) of the stress deviations between the physiological femur model and the optimized total hip arthroplasty (THA) model compared to an unoptimized-THA model could be reduced by 81 % and 66 % in Gruen zone (GZ) 6 and 7. However, the method relies on homogenized finite element (FE) models that only use a simplified representation of the microstructural geometry of the bone and implant. The topology-optimized hip endoprosthesis with graded lattice structures was synthesized using algorithmic design and analyzed in a virtual implanted state using micro-finite element (micro-FE) analysis to validate the optimization method. Homogenized FE and micro-FE models were compared based on averaged von Mises stresses in multiple regions of interest. A strong correlation (CCC > 0.97) was observed, indicating that optimizing homogenized lattice structures yields reliable outcomes. The graded implant was additively manufactured to ensure the topology-optimized result's feasibility.
Collapse
Affiliation(s)
- Patrik Müller
- Institute of Product Development, Leibniz University of Hannover, Garbsen, 30823, Germany.
| | - Alexander Synek
- TU Wien, Institute for Lightweight Design and Structural Biomechanics, Vienna, 1060, Austria
| | - Timo Stauß
- Institute of Product Development, Leibniz University of Hannover, Garbsen, 30823, Germany
| | - Carl Steinnagel
- Institute of Product Development, Leibniz University of Hannover, Garbsen, 30823, Germany
| | - Tobias Ehlers
- Institute of Product Development, Leibniz University of Hannover, Garbsen, 30823, Germany
| | | | - Dieter Pahr
- TU Wien, Institute for Lightweight Design and Structural Biomechanics, Vienna, 1060, Austria
- Division Biomechanics, Karl Landsteiner University of Health Sciences, Krems, 3500, Austria
| | - Roland Lachmayer
- Institute of Product Development, Leibniz University of Hannover, Garbsen, 30823, Germany
| |
Collapse
|
2
|
Bachmann S, Pahr DH, Synek A. Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 236:107549. [PMID: 37084528 DOI: 10.1016/j.cmpb.2023.107549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2022] [Revised: 03/23/2023] [Accepted: 04/12/2023] [Indexed: 05/03/2023]
Abstract
BACKGROUND AND OBJECTIVE Measuring physiological loading conditions in vivo can be challenging, as methods are invasive or pose a high modeling effort. However, the physiological loading of bones is also imprinted in the bone microstructure due to bone (re)modeling. This information can be retrieved by inverse bone remodeling (IBR). Recently, an IBR method based on micro-finite-element (µFE) modeling was translated to homogenized-FE (hFE) to decrease computational effort and tested on the distal radius. However, this bone has a relatively simple geometry and homogeneous microstructure. Therefore, the objective of this study was to assess the agreement of hFE-based IBR with µFE-based IBR to predict hip joint loading from the head of the femur; a bone with more complex loading as well as more heterogeneous microstructure. METHODS hFE-based IBR was applied to a set of 19 femoral heads using four different material mapping laws. One model with a single homogeneous material for both trabecular and cortical volume and three models with a separated cortex and either homogeneous, density-dependent inhomogeneous, or density and fabric-dependent orthotropic material. Three different evaluation regions (full bone, trabecular bone only, head region only) were defined, in which IBR was applied. µFE models were created for the same bones, and the agreement of the predicted hip joint loading history obtained from hFE and µFE models was evaluated. The loading history was discretized using four unit load cases. RESULTS The computational time for FE solving was decreased on average from 500 h to under 1 min (CPU time) when using hFE models instead of µFE models. Using more information in the material model in the hFE models led to a better prediction of hip joint loading history. Inhomogeneous and inhomogeneous orthotropic models gave the best agreement to µFE-based IBR (RMSE% <14%). The evaluation region only played a minor role. CONCLUSIONS hFE-based IBR was able to reconstruct the dominant joint loading of the femoral head in agreement with µFE-based IBR and required considerably lower computational effort. Results indicate that cortical and trabecular bone should be modeled separately and at least density-dependent inhomogeneous material properties should be used with hFE models of the femoral head to predict joint loading.
Collapse
Affiliation(s)
- Sebastian Bachmann
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Gumpendorfer Straße 7, Vienna 1060, Austria.
| | - Dieter H Pahr
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Gumpendorfer Straße 7, Vienna 1060, Austria; Division Biomechanics, Karl Landsteiner University of Health Sciences, Dr. Karl-Dorrek-Straße 30, Krems 3500, Austria
| | - Alexander Synek
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Gumpendorfer Straße 7, Vienna 1060, Austria
| |
Collapse
|
3
|
Wang Y, Yang D, Zhao L, Xiao B, Ma Q, Dong L, Li X, Zhang Y, Zhang Z, Luo Y. Finite Element Analysis of Mechanical Characteristics of Internal Fixation for Treatment of Proximal Femoral Osteolytic Lesions in Children. Orthop Surg 2023; 15:628-638. [PMID: 36597674 PMCID: PMC9891963 DOI: 10.1111/os.13591] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/29/2022] [Revised: 10/13/2022] [Accepted: 10/19/2022] [Indexed: 01/05/2023] Open
Abstract
OBJECTIVES Clinically, it is very difficult to prevent pathological fracture caused by high recurrence rate of osteolytic disease of proximal femur in children. At present, there is no consensus in clinical studies of which internal fixation method can significantly reduce the probability of recurrence of pathological fracture. The study aims to research the mechanical properties of different internal fixations in the treatment of osteolytic lesions of proximal femur in children by finite element analysis, and to find out the optimal treatment. METHODS Based on finite element analysis, the osteolytic disease models of the femoral neck and intertrochanter in a child (8-year-old, boy) were established respectively, and different internal fixation models (plate and titanium elastic intramedullary nails, TENs) were assembled. For the osteolytic lesion of the femoral neck: model A1 was assembled with a plate; model A2 with two TENs crossing the physis; model A3 with two TENs without crossing the physis. And for pertrochanteric osteolytic lesion: model B1 was assembled with a plate, model B2 with two TENs crossing the physis and model B3 with two TENs without crossing the physis. The Eccentric bearing load, torsional restraintal restraint of calcar femorale and composite load were analyzed for each models. RESULTS When the yield strain of each model is reached, the stress concentration points are located in the proximal and distal femoral calcar. In the model of femoral neck lesions, the failure load of model A1 and model A2 are the same (1250 N), and the failure load of model A3 (980 N) is significantly lower than that of the former two; in the model of intertrochanteric lesions, the failure load of model B2 is the largest (1350 N), and the failure load of model B1 (1220 N) is lower than that of model B3 (1260 N), but both are smaller than that of model B2. CONCLUSION Through finite element analysis, TENs through the epiphyseal plate, is found to be the better internal fixation method for femoral neck lesions and intertrochanteric lesions under two different working conditions. The results of clinical correlation study provide new biomechanical information for orthopedic doctors to consider different treatment options for osteolytic lesions of proximal femur.
Collapse
Affiliation(s)
- Yi‐chen Wang
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Dan Yang
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Li‐hua Zhao
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Bo Xiao
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Qi‐chao Ma
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Liang‐chao Dong
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Xiang Li
- School of Mechanical EngineeringShanghai Jiao Tong UniversityShanghaiChina
| | - Yu‐guo Zhang
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Zi‐ming Zhang
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| | - Yi Luo
- Department of Orthopedics, Shanghai Children's Hospital, School of medicineShanghai Jiao Tong UniversityShanghaiChina
| |
Collapse
|
4
|
A Density-Dependent Target Stimulus for Inverse Bone (Re)modeling with Homogenized Finite Element Models. Ann Biomed Eng 2022; 51:925-937. [PMID: 36418745 PMCID: PMC10122636 DOI: 10.1007/s10439-022-03104-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Accepted: 10/17/2022] [Indexed: 11/25/2022]
Abstract
AbstractInverse bone (re)modeling (IBR) can infer physiological loading conditions from the bone microstructure. IBR scales unit loads, imposed on finite element (FE) models of a bone, such that the trabecular microstructure is homogeneously loaded and the difference to a target stimulus is minimized. Micro-FE (µFE) analyses are typically used to model the microstructure, but computationally more efficient, homogenized FE (hFE) models, where the microstructure is replaced by an equivalent continuum, could be used instead. However, also the target stimulus has to be translated from the tissue to the continuum level. In this study, a new continuum-level target stimulus relating relative bone density and strain energy density is proposed. It was applied using different types of hFE models to predict the physiological loading of 21 distal radii sections, which was subsequently compared to µFE-based IBR. The hFE models were able to correctly identify the dominant load direction and showed a high correlation of the predicted forces, but mean magnitude errors ranged from − 14.7 to 26.6% even for the best models. While µFE-based IBR can still be regarded as a gold standard, hFE-based IBR enables faster predictions, the usage of more sophisticated boundary conditions, and the usage of clinical images.
Collapse
|
5
|
A novel specimen shape for measurement of linear strain fields by means of digital image correlation. Sci Rep 2021; 11:17515. [PMID: 34471200 PMCID: PMC8410939 DOI: 10.1038/s41598-021-97085-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2021] [Accepted: 08/13/2021] [Indexed: 11/08/2022] Open
Abstract
Strains on the surface of engineering structures or biological tissues are non-homogeneous. These strain fields can be captured by means of Digital Image Correlation (DIC). However, DIC strain field measurements are prone to noise and filtering of these fields influences measured strain gradients. This study aims to design a novel tensile test specimen showing two linear gradients, to measure full-field linear strain measurements on the surface of test specimens, and to investigate the accuracy of DIC strain measurements globally (full-field) and locally (strain gauges' positions), with and without filtering of the DIC strain fields. Three materials were employed for this study: aluminium, polymer, and bovine bone. Normalized strain gradients were introduced that are load independent and evaluated at two local positions showing 3.6 and 6.9% strain change per mm. Such levels are typically found in human bones. At these two positions, two strain gauges were applied to check the experimental strain magnitudes. A third strain gauge was applied to measure the strain in a neutral position showing no gradient. The accuracy of the DIC field measurement was evaluated at two deformation stages (at [Formula: see text] 500 and 1750 μstrain) using the root mean square error (RMSE). The RMSE over the two linear strain fields was less than 500 μstrain for both deformation stages and all materials. Gaussian low-pass filter (LPF) reduced the DIC noise between 25% and 64% on average. As well, filtering improved the accuracy of the local normalized strain gradients measurements with relative difference less than 20% and 12% for the high- and low-gradient, respectively. In summary, a novel specimen shape and methodological approach are presented which are useful for evaluating and improving the accuracy of the DIC measurement where non-homogeneous strain fields are expected such as on bone tissue due to their hierarchical structure.
Collapse
|
6
|
Synek A, Dunmore CJ, Kivell TL, Skinner MM, Pahr DH. Inverse remodelling algorithm identifies habitual manual activities of primates based on metacarpal bone architecture. Biomech Model Mechanobiol 2019; 18:399-410. [PMID: 30413983 PMCID: PMC6418057 DOI: 10.1007/s10237-018-1091-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2018] [Accepted: 10/29/2018] [Indexed: 12/30/2022]
Abstract
Previously, a micro-finite element (micro-FE)-based inverse remodelling method was presented in the literature that reconstructs the loading history of a bone based on its architecture alone. Despite promising preliminary results, it remains unclear whether this method is sensitive enough to detect differences of bone loading related to pathologies or habitual activities. The goal of this study was to test the sensitivity of the inverse remodelling method by predicting joint loading histories of metacarpal bones of species with similar anatomy but clearly distinct habitual hand use. Three groups of habitual hand use were defined using the most representative primate species: manipulation (human), suspensory locomotion (orangutan), and knuckle-walking locomotion (bonobo, chimpanzee, gorilla). Nine to ten micro-computed tomography scans of each species ([Formula: see text] in total) were used to create micro-FE models of the metacarpal head region. The most probable joint loading history was predicted by optimally scaling six load cases representing joint postures ranging from [Formula: see text] (extension) to [Formula: see text] (flexion). Predicted mean joint load directions were significantly different between knuckle-walking and non-knuckle-walking groups ([Formula: see text]) and in line with expected primary hand postures. Mean joint load magnitudes tended to be larger in species using their hands for locomotion compared to species using them for manipulation. In conclusion, this study shows that the micro-FE-based inverse remodelling method is sensitive enough to detect differences of joint loading related to habitual manual activities of primates and might, therefore, be useful for palaeoanthropologists to reconstruct the behaviour of extinct species and for biomedical applications such as detecting pathological joint loading.
Collapse
Affiliation(s)
- Alexander Synek
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9/BE, Vienna, Austria.
| | - Christopher J Dunmore
- Animal Postcranial Evolution Lab, Skeletal Biology Research Centre, School of Anthropology and Conservation, University of Kent, Canterbury, UK
| | - Tracy L Kivell
- Animal Postcranial Evolution Lab, Skeletal Biology Research Centre, School of Anthropology and Conservation, University of Kent, Canterbury, UK
- Department of Human Evolution, Max Plank Institute for Evolutionary Anthropology, Leipzig, Germany
| | - Matthew M Skinner
- Animal Postcranial Evolution Lab, Skeletal Biology Research Centre, School of Anthropology and Conservation, University of Kent, Canterbury, UK
- Department of Human Evolution, Max Plank Institute for Evolutionary Anthropology, Leipzig, Germany
| | - Dieter H Pahr
- Institute of Lightweight Design and Structural Biomechanics, TU Wien, Getreidemarkt 9/BE, Vienna, Austria
- Department of Anatomy and Biomechanics, Karl Landsteiner Private University of Health Sciences, Krems an der Donau, Austria
| |
Collapse
|