Experimental study of cancellous bone under large strains and a constitutive probabilistic model

Effect of geometrical design variables on implantation configuration and fixation stiffness of titling bone anchors: A parametric finite element study

The mechanical interaction of a tilting anchor and cancellous bones of various densities was simulated using finite element modeling. The model enjoyed a sophisticated representation of the bone, as an elasto-plastic material with large deformation capability. The anchor's tilting action during implantation phase, as well as its fixation stiffness during pull-out test, were predicted by the model and a parametric study was performed to investigate the effects of the anchor's distal width and corner fillet radius, on these measures. The model predictions were validated against the results of an experimental test on ovine humerus specimens. The model could reasonably reproduce the tilting action of the anchor during the implantation phase. Comparison of the model predictions with the experimental results revealed similar trends during both the implantation and the pull-out phases, but smaller displacement magnitudes (end points: 1.4 vs. 2.1 mm and 4.6 vs. 5.2 mm, respectively). The results of the parametric study indicated substantial increase in the fixation stiffness with increasing bone density. Reducing the distal width and increasing the fillet radius improved the anchor's implantation configuration and fixation stiffness in low-density bones. For high-density bone applications, however, a larger distal width was favored for improving the fixation stiffness.

Gaussian random field-based characterization and reconstruction of cancellous bone microstructure considering the constraint of correlation structure

The macro scale physical properties of cancellous bone materials are governed by the microstructural features, which is of great significance for the multi-scale research of cancellous bone and the inverse design of bone-mimicking materials. Therefore, it is essential to characterize the natural cancellous bone samples, and reconstruct the microstructures with the biomimetic osteointegration and mechanical properties. In this research, a novel approach for the characterization and reconstruction of cancellous bone was proposed, based on the medical image analysis and anisotropic three-dimensional Gaussian random field (GRF). The geometric similarity, i.e. the interface curvature distribution (ISD), was meticulously studied, which is important to the osteointegration ability. And the mechanical properties were validated by the stress-strain curves under the large compressive strain simulated by the smoothed particle hydrodynamic (SPH) method. In addition, the effects of the generation parameters of GRF-based biomimetic microstructures on the apparent properties were analyzed. The ISD results demonstrated that both GRF and micro-CT groups had the similar columnar morphological properties, while the latter had more hyperbolic features. And it was found that the GRF-based biomimetic microstructures and the natural bone samples based on micro-CT (MCT) had the similar failure mode. The concordance correlation coefficient between MCT and GRF pairs was 0.8685, with a Pearson ρ value of 0.8804, and significance level p<0.0001. The Bland-Altman LoA was 0.1647 MPa with 95 % (1.96SD) lower and upper bound value between -0.2892 and 0.6185 MPa. The two groups had almost the same elastic modulus with the mean absolute percentage error (MAPE) of 7.84 %. While the yield stress and total conversion energy of the GRF-based samples were lower than those of the natural bone samples, and the MAPE were 16.99 % and 16.27 %, respectively. Although it meant the lower structural efficiency, the huge design space of this approach and advanced 3D printing technology can provide great potential for the design of orthopedic implants.

Smoothed particle hydrodynamics implementation to enhance vertebral fracture finite element model in a cervical spine segment under compression

Spinal cord injuries (SCIs) can arise from compression loading when a vertebra fractures and bone fragments are pushed into the spinal canal. Experimental studies have demonstrated the importance of both fracture initiation and post-fracture response in the investigation of vertebral fractures and spinal canal occlusion resulting from compression. Finite element models, such as the Global Human Body Models Consortium (GHBMC) model, focused on predicting the initiation location of fractures using element erosion to model hard tissue fracture. However, the element erosion method resulted in a loss of material and structural support during compression, which limited the ability of the model to predict the post-fracture response. The current study aimed to improve the post-fracture response by combining strain-based element erosion with smoothed particle hydrodynamics (SPH) to preserve the volume of the trabecular bone during compression fracture. The proposed implementation was evaluated using a model comprising two functional spinal units (FSUs) (C5-C6-C7) extracted from the GHBMC 50th percentile male model, and loaded under central compression. The original and enhanced models were compared to experimental force-displacement data and measured occlusion of the spinal canal. The enhanced model with SPH improved the shape and magnitude of the force-displacement response to be in good agreement with the experimental data. In contrast to the original model, the enhanced SPH model demonstrated occlusion on the same order of magnitude as reported in the experiments. The SPH implementation improved the post-fracture response by representing the damaged material post-fracture, providing structural support throughout compression loading and material flow leading to occlusion.

A Damage Model to Trabecular Bone and Similar Materials: Residual Resource, Effective Elasticity Modulus, and Effective Stress under Uniaxial Compression

Experimental research of bone strength remains costly and limited for ethical and technical reasons. Therefore, to predict the mechanical state of bone tissue, as well as similar materials, it is desirable to use computer technology and mathematical modeling. Yet, bone tissue as a bio-mechanical object with a hierarchical structure is difficult to analyze for strength and rigidity; therefore, empirical models are often used, the disadvantage of which is their limited application scope. The use of new analytical solutions overcomes the limitations of empirical models and significantly improves the way engineering problems are solved. Aim of the paper: the development of analytical solutions for computer models of the mechanical state of bone and similar materials. Object of research: a model of trabecular bone tissue as a quasi-brittle material under uniaxial compression (or tension). The new ideas of the fracture mechanics, as well as the methods of mathematical modeling and the biomechanics of bone tissues were used in the work. Compression and tension are considered as asymmetric mechanical states of the material. Results: a new nonlinear function that simulates both tension and compression is justified, analytical solutions for determining the effective and apparent elastic modulus are developed, the residual resource function and the damage function are justified, and the dependences of the initial and effective stresses on strain are obtained. Using the energy criterion, it is proven that the effective stress continuously increases both before and after the extremum point on the load-displacement plot. It is noted that the destruction of bone material is more likely at the inflection point of the load-displacement curve. The model adequacy is explained by the use of the energy criterion of material degradation. The results are consistent with the experimental data available in the literature.

Influence of the fixation technique on the mechanical properties of human cancellous bone of the femoral head

BACKGROUND

The femoral head is of central importance for the force transmission from the suprapelvic body mass to the lower extremity. However, the condition of the subcortical bone and its mechanical properties in case of pathological changes due to coxarthrosis or femoral head necrosis differ from the healthy condition.

METHODS

Fresh femoral heads were gathered during hip total endoprosthesis surgeries and cylindrical cancellous bone samples were extracted with a hollow drill. By means of a uniaxial tensile-compression test system, the compressive strength was determined for two different specimen types (fresh and 24 h storage in acetone). Exemplary tests on an exceptionally large femoral head were performed to compare properties of fresh, fresh-deep-frozen and acetone-stored samples.

FINDINGS

The deformation behaviour and the material parameters determined were very heterogeneous. For most of the specimens, a destructive material test was successfully carried out, i.e. the compressive strength was determined. The average strength of fresh specimens was slightly higher than that of acetone specimens. On the other hand, the average Young's modulus of the acetone specimens was higher than that of the fresh specimens.

INTERPRETATION

The lower Young's moodulus of the fresh samples compared to the acetone samples could indicate a causal effect of the degreasing influence of the acetone. The partly considerable individual differences in compressive strength and failure compression can have patient-specific influencing factors such as constitution and physical fitness as well as causes in the initial pathological condition.

An Analytical Approach to Study the Intraoperative Fractures of Femoral Shaft During Total Hip Arthroplasty

An analytical approach which is popular in micromechanical studies has been extended to the solution for the interference fit problem of the femoral stem in cementless total hip arthroplasty (THA). The multiple inhomogeneity problem of THA in transverse plane, including an elliptical stem, a cortical wall, and a cancellous layer interface, was formulated using the equivalent inclusion method (EIM) to obtain the induced interference elastic fields. Results indicated a maximum interference fit of about 210 μm before bone fracture, predicted based on the Drucker–Prager criterion for a partially reamed section. The cancellous layer had a significant effect on reducing the hoop stresses in the cortical wall; the maximum press fit increased to as high as 480 μm for a 2 mm thick cancellous. The increase of the thickness and the mechanical quality, i.e., stiffness and strength, of the cortical wall also increased the maximum interference fit before fracture significantly. No considerable effect was found for the implant material on the maximum allowable interference fit. It was concluded that while larger interference fits could be adapted for younger patients, care must be taken when dealing with the elderly and those suffering from osteoporosis. A conservative reaming procedure is beneficial for such patients; however, in order to ensure sufficient primary stability without risking bone fracture, a preoperative analysis might be necessary.

Porous low modulus Ti40Nb compacts with electrodeposited hydroxyapatite coating for biomedical applications

Porous ß-type non-toxic Ti40Nb alloy was prepared by compaction of mechanically alloyed powder mixed with NaCl or Mg particles as space-holder material. The compacts with porosity of 36-80% demonstrated a very low Young's modulus of ~1.5-3 GPa and compression strength of ~10-35 MPa, which is suitable for potential implant material application. Porous samples were electrochemically covered with hydroxyapatite. The influence of the deposition time and of the electrolyte concentrations on the morphology of the hydroxyapatite coating was studied. It is demonstrated that a homogenous coating of hydroxyapatite crystals with different shape and size can be obtained on the surface of the porous samples.