Comparison of simplified bone-screw interface models in materially nonlinear µFE simulations.

Micro finite-element (µFE) simulations serve as a crucial research tool to assist laboratory experiments in the biomechanical assessment of screw anchorage in bone. However, accurately modelling the interface between bone and screw threads at the microscale poses a significant challenge. Currently, the gold-standard approach involves employing computationally intensive physical contact models to simulate this interface. This study compared nonlinear µFE predictions of deformations, whole-construct stiffness, maximum force and damage patterns of three different computationally efficient simplified interface approaches to the general contact interface in Abaqus Explicit, which was defined as gold-standard and reference model. The µCT images (resolution: 32.8 µm) of two human radii with varying bone volume fractions were utilized and a screw was virtually inserted up to 50% and 100% of the volar-dorsal cortex distance. Materially nonlinear µFE models were generated and loaded in tension, compression and shear. In a first step, the common simplification of using a fully-bonded interface was compared to the general contact interface, revealing overestimations of whole-construct stiffness (19% on average) and maximum force (26% on average), along with inaccurate damage pattern replications. To enhance predictions, two additional simplified interface models were compared: tensionally strained element deletion (TED) and a novel modification of TED (TED-M). TED deletes interface elements strained in tension based on a linear-elastic simulation before the actual simulation. TED-M extends the remaining contact interface of TED by incorporating neighboring elements to the contact area. Both TED and TED-M reduced the errors in whole-construct stiffness and maximum force and improved the replication of the damage distributions in comparison to the fully-bonded approach. TED was better in predicting whole-construct stiffness (average error of 1%), while TED-M showed lowest errors in maximum force (1% on average). In conclusion, both TED and TED-M offer computationally efficient alternatives to physical contact modelling, although the fully-bonded interface may deliver sufficiently accurate predictions for many applications.

Cortical and trabecular mechanical properties in the femoral neck vary differently with changes in bone mineral density.

Graphical Abstract.

Critical loss of primary implant stability in osteosynthesis locking screws under cyclic overloading.

Primary implant stability, which refers to the stability of the implant during the initial healing period is a crucial factor in determining the long-term success of the implant and lays the foundation for secondary implant stability achieved through osseointegration. Factors affecting primary stability include implant design, surgical technique, and patient-specific factors like bone quality and morphology. In vivo, the cyclic nature of anatomical loading puts osteosynthesis locking screws under dynamic loads, which can lead to the formation of micro cracks and defects that slowly degrade the mechanical connection between the bone and screw, thus compromising the initial stability and secondary stability of the implant. Monotonic quasi-static loading used for testing the holding capacity of implanted screws is not well suited to capture this behavior since it cannot capture the progressive deterioration of peri­implant bone at small displacements. In order to address this issue, this study aims to determine a critical point of loss of primary implant stability in osteosynthesis locking screws under cyclic overloading by investigating the evolution of damage, dissipated energy, and permanent deformation. A custom-made test setup was used to test implanted 2.5 mm locking screws under cyclic overloading test. For each loading cycle, maximum forces and displacement were recorded as well as initial and final cycle displacements and used to calculate damage and energy dissipation evolution. The results of this study demonstrate that for axial, shear, and mixed loading significant damage and energy dissipation can be observed at approximately 20 % of the failure force. Additionally, at this load level, permanent deformations on the screw-bone interface were found to be in the range of 50 to 150 mm which promotes osseointegration and secondary implant stability. This research can assist surgeons in making informed preoperative decisions by providing a better understanding of the critical point of loss of primary implant stability, thus improving the long-term success of the implant and overall patient satisfaction.

Development of a density-based topology optimization of homogenized lattice structures for individualized hip endoprostheses and validation using micro-FE.

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.

Three-Dimensional Osseointegration Patterns of Cementless Femoral Stems: An ex Vivo Study with High-Resolution Imaging and Histological Evaluation.

BACKGROUND: Osseointegration is essential for the long-term survival of cementless femoral stems and is dependent on periprosthetic bone quality and correct implantation technique. The aim of this study was to evaluate the 3-dimensional long-term fixation patterns of, and bone microarchitecture around, cementless hip stems. METHODS: Four specimens with varying degrees of bone quality and fixation characteristics from body donors who had received Alloclassic Zweymüller hip stems during their lifetime (mean time in situ at the time of death: 12.73 years) were evaluated with use of radiographs, high-resolution computed tomography (CT) scans, and hard-tissue histology. The CT voxel size was 85 µm, and the following parameters were calculated: total bone volume, total bone volume fraction, trabecular bone volume, trabecular bone volume fraction, cortical bone volume, cortical bone volume fraction, and cortical thickness. Bone-implant contact and canal fill index values for each Gruen zone of the specimens were calculated with use of histological samples. RESULTS: Femoral stems with apparently good cortical contact on clinical radiographs showed higher values for cortical bone volume, trabecular bone volume, and cortical thickness in the high-resolution CT analysis than femoral stems with apparently weak cortical contact on clinical radiographs. Based on the histological evaluation, the mean bone-implant contact ranged from 22.94% to 57.24% and the mean canal fill index ranged from 52.33% to 69.67% among the specimens. CONCLUSIONS: This study demonstrated different osseointegration patterns of cementless femoral stems on the basis of radiographs, high-resolution CT scans, and histological evaluation. Femora with high cortical bone volume and cortical thickness were associated with higher canal fill indices, whereas femora with low cortical bone volume and cortical thickness had lower canal fill indices and showed a characteristic corner-anchorage pattern. CLINICAL RELEVANCE: Osseointegration patterns and thus the long-term survival of cementless femoral stems are dependent on cortical bone volume and cortical thickness.

Trabecular architecture of the distal femur in extant hominids.

Extant great apes are characterized by a wide range of locomotor, postural and manipulative behaviours that each require the limbs to be used in different ways. In addition to external bone morphology, comparative investigation of trabecular bone, which (re-)models to reflect loads incurred during life, can provide novel insights into bone functional adaptation. Here, we use canonical holistic morphometric analysis (cHMA) to analyse the trabecular morphology in the distal femoral epiphysis of Homo sapiens (n = 26), Gorilla gorilla (n = 14), Pan troglodytes (n = 15) and Pongo sp. (n = 9). We test two predictions: (1) that differing locomotor behaviours will be reflected in differing trabecular architecture of the distal femur across Homo, Pan, Gorilla and Pongo; (2) that trabecular architecture will significantly differ between male and female Gorilla due to their different levels of arboreality but not between male and female Pan or Homo based on previous studies of locomotor behaviours. Results indicate that trabecular architecture differs among extant great apes based on their locomotor repertoires. The relative bone volume and degree of anisotropy patterns found reflect habitual use of extended knee postures during bipedalism in Homo, and habitual use of flexed knee posture during terrestrial and arboreal locomotion in Pan and Gorilla. Trabecular architecture in Pongo is consistent with a highly mobile knee joint that may vary in posture from extension to full flexion. Within Gorilla, trabecular architecture suggests a different loading of knee in extension/flexion between females and males, but no sex differences were found in Pan or Homo, supporting our predictions. Inter- and intra-specific variation in trabecular architecture of distal femur provides a comparative context to interpret knee postures and, in turn, locomotor behaviours in fossil hominins.

High-resolution local trabecular strain within trabecular structure under cyclic loading.

Trabecular bone structure is a complex microstructure consisting of rods and plates, which poses challenges for its mechanical characterization. Digital image correlation (DIC) offers the possibility to characterize the strain response on the surface of trabecular bone. This study employed DIC equipped with a telecentric lens to investigate the strain state of individual trabeculae within their trabecular structure by assessing the longitudinal strain of the trabeculae at both the middle and near the edges of the trabeculae. Due to the high-resolution of the used DIC system, local surface strain of trabeculae was analyzed too. Lastly, the correlation between longitudinal trabecular strain and the orientation and slenderness of the trabeculae was investigated. The results showed that the strain magnification close to the edge of the trabeculae was higher and reached up to 8-folds the strain along the middle of the trabeculae. On the contrary, no strain magnification was found for most of the trabeculae between the longitudinal trabecular strain along the middle of the trabeculae and the globally applied strain. High-resolution full-field strain maps were obtained on the surface of trabeculae showing heterogeneous strain distribution with increasing load. No significant correlation was found between longitudinal trabecular strain and its orientation or slenderness. These findings and the applied methodology can be used to broaden our understanding of the deformation mechanisms of trabeculae within the trabecular network.

The deep trabecular structure of first metacarpals in extant hominids.

OBJECTIVES: Recent studies have associated subarticular trabecular bone distribution in the extant hominid first metacarpal (Mc1) with observed thumb use, to infer fossil hominin thumb use. Here, we analyze the entire Mc1 to test for interspecific differences in: (1) the absolute volume of trabecular volume fraction, (2) the distribution of the deeper trabecular network, and (3) the distribution of trabeculae in the medullary cavity, especially beneath the Mc1 disto-radial flange. MATERIALS AND METHODS: Trabecular bone was imaged using micro-computed tomography in a sample of Homo sapiens (n = 11), Pan paniscus (n = 10), Pan troglodytes (n = 11), Gorilla gorilla (n = 10) and Pongo sp., (n = 7). Using Canonical Holistic Morphometric Analysis (cHMA), we tested for interspecific differences in the trabecular bone volume fraction (BV/TV) and its relative distribution (rBV/TV) throughout the Mc1, including within the head, medullary cavity, and base. RESULTS: P. paniscus had the highest, and H. sapiens the lowest, BV/TV relative to other species. rBV/TV distribution statistically distinguished the radial concentrations and lack of medullary trabecular bone in the H. sapiens Mc1 from all other hominids. H. sapiens and, to a lesser extent, G. gorilla also had a significantly higher trabecular volume beneath the disto-radial flange relative to other hominids. DISCUSSION: These results are consistent with differences in observed thumb use in these species and may also reflect systemic differences in bone volume fraction. The trabecular bone extension into the medullary cavity and concentrations beneath the disto-radial flange may represent crucial biomechanical signals that will aid in the inference of fossil hominin thumb use.

The 2-layer elasto-visco-plastic rheological model for the material parameter identification of bone tissue extended by a damage law.

The response of bone tissue to mechanical load is complex and includes plastic hardening, viscosity and damage. The quantification of these effects plays a mayor role in bone research and in biomechanical clinical trials as to better understand related diseases. In this study, the damage growth in individual wet human trabeculae subjected to cyclic overloading is quantified by inverse rheological modeling. Therefore, an already published rheological material model, that includes linear elasticity, plasticity and viscosity is extended by a damage law. The model is utilized in an optimization process to identify the corresponding material parameters and damage growth in single human trabeculae under tensile load. Results show that the damage model is leading to a better fit of the test data with an average root-mean-square-error (RMSE) of 2.52 MPa compared to the non-damage model with a RMSE of 3.03 MPa. Although this improvement is not significant, the damage model qualitatively better represents the data as it accounts for the visible stiffness reduction along the load history. It returns realistic stiffness values of 11.92 GPa for the instantaneous modulus and 5.73 GPa for the long term modulus of wet trabecular human bone. Further, the growth of damage in the tissue along the load history is substantial, with values above 0.8 close to failure. The relative loss of stiffness per cycle is in good agreement with comparable literature. Inverse rheological modeling proves to be a valuable tool for quantifying complex constitutive behavior from a single mechanical measurement. The evolution of damage in the tissue can be identified continuously over the load history and separated from other effects.

A mechanically validated open-source silicone model for the training of gastric perforation sewing.

BACKGROUND: Gastrointestinal perforation is commonly seen in emergency departments. The perforation of the stomach is an emergency situation that requires immediate surgical treatment. The necessary surgical skills require regular practical training. Owing to patient`s safety, in vivo training opportunities in medicine are restricted. Animal tissue especially porcine tissue, is commonly used for surgical training. Due to its limiting factors, artificial training models are often to be preferred. Many artificial models are on the market but to our knowledge, none that mimic the haptic- and sewing properties of a stomach wall at the same time. In this study, an open source silicone model of a gastric perforation for training of gastric sewing was developed that attempts to provide realistic haptic- and sewing behaviour. METHODS: To simulate the layered structure of the human stomach, different silicone materials were used to produce three different model layups. The production process was kept as simple as possible to make it easily reproducible. A needle penetration setup as well as a systematic haptic evaluation were developed to compare these silicone models to a real porcine stomach in order to identify the most realistic model. RESULTS: A silicone model consisting of three layers was identified as being the most promising and was tested by clinical surgeons. CONCLUSIONS: The presented model simulates the sewing characteristics of a human stomach wall, is easily reproducible at low-costs and can be used for practicing gastric suturing techniques. TRIAL REGISTRATIONS: Not applicable.

Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy.

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.

Accuracy of osseointegrated screw-bone construct stiffness and peri-implant loading predicted by homogenized FE models relative to micro-FE models.

Computational predictions of stiffness and peri-implant loading of screw-bone constructs are highly relevant to investigate and improve bone fracture fixations. Homogenized finite element (hFE) models have been used for this purpose in the past, but their accuracy has been questioned given the numerous simplifications, such as neglecting screw threads and modelling the trabecular bone structure as a continuum. This study aimed to investigate the accuracy of hFE models of an osseointegrated screw-bone construct when compared to micro-FE models considering the simplified screw geometry and different trabecular bone material models. Micro-FE and hFE models were created from 15 cylindrical bone samples with a virtually inserted, osseointegrated screw (fully bonded interface). Micro-FE models were created including the screw with threads (=reference models) and without threads to quantify the error due to screw geometry simplification. In the hFE models, the screws were modelled without threads and four different trabecular bone material models were used, including orthotropic and isotropic material derived from homogenization with kinematic uniform boundary conditions (KUBC), as well as from periodicity-compatible mixed uniform boundary conditions (PMUBC). Three load cases were simulated (pullout, shear in two directions) and errors in the construct stiffness and the volume average strain energy density (SED) in the peri-implant region were evaluated relative to the micro-FE model with a threaded screw. The pooled error caused by only omitting screw threads was low (max: 8.0%) compared to the pooled error additionally including homogenized trabecular bone material (max: 92.2%). Stiffness was predicted most accurately using PMUBC-derived orthotropic material (error: -0.7 ± 8.0%) and least accurately using KUBC-derived isotropic material (error: +23.1 ± 24.4%). Peri-implant SED averages were generally well correlated (R2 ≥ 0.76), but slightly over- or underestimated by the hFE models and SED distributions were qualitatively different between hFE and micro-FE models. This study suggests that osseointegrated screw-bone construct stiffness can be predicted accurately using hFE models when compared to micro-FE models and that volume average peri-implant SEDs are well correlated. However, the hFE models are highly sensitive to the choice of trabecular bone material properties. PMUBC-derived isotropic material properties represented the best trade-off between model accuracy and complexity in this study.

Comparison of linear and nonlinear stepwise µFE displacement predictions to digital volume correlation measurements of trabecular bone biopsies.

Digital volume correlation (DVC) enables to evaluate the ability of µFE models in predicting experimental results on the mesoscale. In this study predicted displacement fields of three different linear and materially nonlinear µFE simulation methods were compared to DVC measured displacement fields at specific load steps in the elastic regime (StepEl) and after yield (StepUlt). Five human trabecular bone biopsies from a previous study were compressed in several displacement steps until failure. At every compression step, µCT images (resolution: 36 µm) were recorded. A global DVC algorithm was applied to compute the displacement fields at all loading steps. The unloaded 3D images were then used to generate homogeneous, isotropic, linear and materially nonlinear µFE models. Three different µFE simulation methods were used: linear (L), nonlinear (NL), and nonlinear stepwise (NLS). Regarding L and NL, the boundary conditions were derived from the interpolated displacement fields at StepEl and StepUlt, while for the NLS method nonlinear changes of the boundary conditions of the experiments were captured using the DVC displacement field of every available load step until StepEl and StepUlt. The predicted displacement fields of all µFE simulation methods were in good agreement with the DVC measured displacement fields (individual specimens: R2>0.83 at StepEl and R2>0.59 at StepUlt; pooled data: R2>0.97 at StepEl and R2>0.92 at StepUlt). At StepEl, all three simulation methods showed similar intercepts, slopes, and coefficients of determination while the nonlinear µFE models improved the prediction of the displacement fields slightly in all Cartesian directions at StepUlt (individual specimens: L: R2>0.59 and NL, NLS: R2>0.68; pooled data: L: R2>0.92 and NL, NLS: R2>0.94). Damaged/overstrained elements in L, NL, and NLS occurred at similar locations but the number of overstrained elements was overestimated when using the L simulation method. Considering the increased solving time of the nonlinear µFE models as well as the acceptable performance in displacement prediction of the linear µFE models, one can conclude that for similar use cases linear µFE models represent the best compromise between computational effort and accuracy of the displacement field predictions.

Feasibility of aluminum phantom radiography for osteoporosis detection in postmenopausal women with a fragility fracture of the distal radius compared to DXA and HR-pQCT.

Recently, promising results have been reported for detection of osteoporosis with use of an aluminum phantom. Therefore, the aim of this study was to evaluate the feasibility of radiography-based bone mineral density (BMD) measurement using a graded aluminum phantom. This study included 27 postmenopausal women with a distal radius fracture. Aluminum phantom radiography of the healthy radius was conducted as well as high-resolution peripheral quantitative computed tomography (HR-pQCT) measurement of the ultradistal radius and dual energy X-ray absorptiometry (DXA) of the radius, spine, and hip. A strong correlation was observed between aluminum phantom radiography-based mean gray value (mGV) and DXA-derived BMD, especially for the ultradistal radius (ρ = 0.75; p < 0.001). A moderate correlation for the femoral neck (ρ = 0.61 and p < 0.001) between modalities was found. Radius mGV and HR-pQCT-derived BMD only showed a moderate correlation (ρ = 0.48; p < 0.09). Aluminum phantom radiography might serve as a cost efficient, highly available, low-radiation dose screening, and diagnostic method for osteoporosis additively to DXA measurements. Especially, an application in areas with constrained DXA availability and such as preoperative trauma settings would be beneficial. However, further investigation and assessment of specificity and sensitivity is needed.

The predictive ability of a QCT-FE model of the proximal femoral stiffness under multiple load cases is strongly influenced by experimental uncertainties.

Despite significant improvements in terms of the predictive ability of Quantitative Computed Tomography based Finite Element (QCT-FE) models in estimating femoral strength (fracture load and stiffness), no substantial clinical adoption of this method has taken place to date. Narrowing the wide variability of FE results by standardizing the methodology and validation protocols, as well as reducing the uncertainties in the FEA process have been proposed as routes towards improved reliability. The aim of this study was to: First, validate a QCT-FE model of proximal femoral stiffness in multiple stance load cases, and second, using a parametric approach, determine the influence of select experimental and modeling parameters on the predictive ability of our model. Ten fresh frozen human femoral samples were tested in neutral stance, 15° adducted and 15° abducted load cases. Voxel-based linear-elastic QCT-FE models of the samples were generated to predict the models' stiffness values in all load cases. The base FE models were validated against the experimental results using linear regression. Thirty six deviated models were created using the minimum and maximum values of experiment-based "plausible range" for 18 parameters in 4 categories of embedding, loading, material, and segmentation. The predictive ability of the models were compared in terms of the coefficient of determination (R2) of the linear regression between the measured and predicted stiffness values in all load cases. Our model was capable of capturing 90% of the variation in the experimental stiffness of the samples in neutral stance position (R2 = 0.9, concordance correlation coefficient (CCC) = 0.93, percent root mean squared error (RMSE%) = 8.4%, slope and intercept not significantly different from unity and zero, respectively). Embedding and loading categories strongly affected the predictive ability of the models with an average percent difference in R2 of 4.36% ± 2.77 and 2.96% ± 1.69 for the stance-neutral load case, respectively. The performance of the models were significantly different in adducted and abducted load cases with their R2 dropping to 71% and 70%, respectively. Similarly, off-axes load cases were affected by the parameters differently compared to the neutral load case, with the loading parameter category imposing more than 10% difference on their R2, larger than all other categories. We also showed that automatically selecting the best performing plausible value for each parameter and each sample would result in a perfectly linear correlation (R2> 0.99) between the "tuned" model's predicted stiffness and experimental results. Based on our results, high sensitivity of the model performance to experimental parameters requires extra diligence in modeling the embedding geometry and the loading angles since these sources of uncertainty could dwarf the effects of material modeling and image processing parameters. The results of this study could help in improving the robustness of the QCT-FE models of proximal femur by limiting the uncertainties in the experimental and modeling steps.

A Density-Dependent Target Stimulus for Inverse Bone (Re)modeling with Homogenized Finite Element Models.

Inverse 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.

The biomechanical behavior of 3D printed human femoral bones based on generic and patient-specific geometries.

BACKGROUND: Bone is a highly complex composite material which makes it hard to find appropriate artificial surrogates for patient-specific biomechanical testing. Despite various options of commercially available bones with generic geometries, these are either biomechanically not very realistic or rather expensive. METHODS: In this work, additive manufacturing was used for the fabrication of artificial femoral bones. These were based on CT images of four different commercially available femoral bone surrogates and three human bones with varying bone density. The models were 3D printed using a low-budget fused deposition modeling (FDM) 3D printer and PLA filament. The infill density was mechanically calibrated and varying cortical thickness was used. Compression tests of proximal femora simulating stance were performed and the biomechanical behavior concerning ultimate force, spring stiffness, and fracture pattern were evaluated as well as compared to the results of commercial and cadaveric bones. RESULTS: Regarding the ultimate forces and spring stiffness, the 3D printed analogs showed mechanical behavior closer to their real counterparts than the commercially available polyurethan-based surrogates. Furthermore, the increase in ultimate force with increasing bone density observed in human femoral bones could be reproduced well. Also, the fracture patterns observed match well with fracture patterns observed in human hip injuries. CONCLUSION: Consequently, the methods presented here show to be a promising alternative for artificial generic surrogates concerning femoral strength testing. The manufacturing is straightforward, cheap, and patient-specific geometries are possible.

Influence of aging on mechanical properties of the femoral neck using an inverse method.

Today, we are facing rapid aging of the world population, which increases the incidence of hip fractures. The gold standard of bone strength assessment in the laboratory is micro-computed finite element analysis (µFEA) based on micro-computed tomography (µCT) images. In clinics, the standard method to assess bone fracture risk is based on areal bone mineral density (aBMD), measured by dual-energy X-ray absorptiometry (DXA). In addition, homogenized finite element analysis (hFEA) constructed from quantitative computed tomography reconstructions (QCT) predicts clinical bone strength more accurately than DXA. Despite considerable evidence of degradation of bone material properties with age, in the past fifty years of finite element analysis to predict bone strength, bone material parameters remained independent of age. This study aims to assess the influence of age on apparent modulus, yield stress, and strength predictions of the human femoral neck made by laboratory-available bone volume fraction (BV/TV) and µFEA; and by clinically available DXA and hFEA. Using an inverse method, we test the hypothesis that FEA material parameters are independent of age. Eighty-six human femora were scanned with DXA (aBMD) and with QCT. The femoral necks were extracted and scanned at 16 µm resolution with µCT. The grayscale images were downscaled to 32 µm and 65 µm for linear and non-linear analyses, respectively, and segmented. The µFE solver ParOSolNL (non-linear) and a standard hFEA method were applied to the neck sections with the same material properties for all samples to compute apparent modulus, yield stress, and strength. Laboratory-available BV/TV was a good predictor of apparent modulus (R2 = 0.76), almost as good as µFEA (R2 = 0.79). However, yield stress and strength were better predicted by µFEA (R2 = 0.92, R2 = 0.86, resp.) than BV/TV (R2 = 0.76, R2 = 0.76, resp.). For clinically available variables, prediction of apparent modulus was better with hFEA than aBMD (R2 = 0.67, R2 = 0.58, resp.). hFEA outperformed aBMD for predictions of yield stress (R2 = 0.63 vs R2 = 0.34 for female and R2 = 0.55 for male) and strength (R2 = 0.48 vs R2 = 0.33 for female and R2 = 0.15 for male). The inclusion of age did not improve the multiple linear models for apparent modulus, yield stress, and strength. The resolution of the µFE meshes seems to account for most morphological changes induced by aging. The errors between the simulation and the experiment for apparent modulus, yield stress, and strength were age-independent, suggesting no rationale for correcting tissue material parameters in the current FE analysis of the aging femoral neck.

Fracture toughness determination of porcine muscle tissue based on AQLV model derived viscous dissipated energy.

The ability of soft collagenous tissue (SCT) to withstand propagation of a defect in the presence of a macroscopic crack is termed the 'fracture toughness parameter'. In soft tissues not undergoing significant plastic deformation, it is purported that a considerable amount of additional energy is dissipated during failure processes, due to viscoelasticity. Hence the total work, measured experimentally during failure, is the sum of fracture and viscoelastic energies. Previous authors have aimed to apply constitutive modeling to describe viscoelastic hysteresis for fracture toughness determination with a tendency of models to either over or underestimate the viscous energy. In this study, the fracture toughness of porcine muscle tissue is determined using two strategies. Firstly, it was determined experimentally by calculation of the difference in dissipated energy of notched and unnotched tissue specimens undergoing cyclic 'triangular wave' excitation with increasing strain levels in uniaxial tension. The second strategy involved the extension and use of the adaptive quasi-linear viscoelastic model (AQLV) to model cyclic loading (model parameters were obtained from a previous study) and sequentially the dissipated energy was calculated. The mean value of the dissipated energy based on the AQLV approach was then subtracted from the total dissipated energy of notched porcine muscle tissue samples to determine the fracture toughness. The mean experimental viscous dissipated energy ratio was 0.24 ± 0.04 in the experimental approach, compared to 0.28 ± 0.03 for the AQLV model. Fracture toughness determined experimentally yielded 0.84 ± 0.80 kJ/m2, and 0.71 ± 0.76 kJ/m2 for the AQLV model, without a significant difference (p = 0.87). Hence, the AQLV model enables a reasonable estimation of viscous dissipated energy in porcine muscle tissue with the advantage to perform tests only on notched specimens, instead of testing additional unnotched samples. Moreover, the AQLV model will help to better understand the constitutive viscoelastic behaviour of SCTs and might also serve as a basis for future fracture toughness determination with constitutive model simulations.

Articular contact vs. embedding: Effect of simplified boundary conditions on the stress distribution in the distal radius and volar plate implant loading.

Boundary conditions (BCs) are often simplified in experimental and numerical models simulating distal radius fractures and their treatments. The aim of this study was to investigate the effects of simplified BCs at the radiocarpal joint: (1) on the stress distribution in the intact distal radius, and (2) on the loading of a volar locking plate (VLP) used for distal radius fracture treatment. Finite element models of the distal radius with contact between carpals and cartilage were created as reference models for an intact bone and a fractured bone with VLP treatment. Four models with simplified BCs were compared to these reference models: One with embedding material instead of carpals, one with carpals tied to the radius; each loaded either uniaxially or with statically equivalent loading to the reference model. Differences in distal bone stress distributions and mechanical parameters of the VLP (fracture gap movement, plate peak stresses, distal screw loads) were generally largest for the uniaxially loaded, embedded model (up to 250 % in individual screw loads) and smallest for the model with tied carpals and statically equivalent loads (<25 % for all parameters). Differences were greatly reduced if statically equivalent loads were applied, but subchondral stress peaks were absent without carpals. In conclusion, implementing realistic resultant forces and moments is more important than the exact articular load distribution, but carpal bones should be included if subchondral bone stresses are analyzed. In this case, a tie constraint may replace articular contact modelling with acceptable accuracy if statically equivalent loading is applied.