Stochastic failure analysis of proximal femur using an isogeometric analysis based nonlocal gradient-enhanced damage model

Prediction of the critical energy release rate for rat femoral cortical bone structure under different failure conditions

BACKGROUND AND OBJECTIVE

Critical energy release rate is a global fracture parameter that could be measured during the failing process, and its value may change under different failure conditions even in the same bone structure. The aim of this study was to propose an approach that combined the experimental test and finite element analysis to predict the critical energy release rates in the femoral cortical bone structures under compression and three-point bending loads.

METHODS

Three-point bending and compression experiments and the corresponding fracture simulations were performed on the rat femoral cortical bone structures. Different values of energy release rate were repeatedly assigned to the finite element models to perform fracture simulations, and then the load-displacement curves predicted in each simulation were compared with the experimental data to back-calculate the critical energy release rate.

RESULTS

The predicted data were similar to the experimental results when the calibrated energy release rate was suitable. The results showed that the cortical bone structure occurred shear open failure under compression load, and the predicted critical energy release rate was 0.12 N/mm. The same cortical bone structure occurred tensile open failure under three-point bending load, and the predicted critical energy release rate was 0.16 N/mm.

CONCLUSIONS

The critical energy release rates were different under various failure conditions in one cortical bone structure. A comprehensive analysis from the perspectives of material mechanical properties, failure mode, and damage fracture mechanism was conducted to reveal the reasons for the differences in the critical energy release rate in the cortical bone structure, which provided a theoretical basis for the measurement of the critical energy release rate and the accurate fracture simulation.

Experimental and numerical analysis of the influence of intramedullary nail position on the cut-out phenomenon

BACKGROUND AND OBJECTIVE

Proximal femur fractures, colloquially known as hip fractures, are a common pathology with increasing incidence in the last years due to the enhanced ageing population. Regarding the extracapsular fracture, the treatment for this pathology consists of a fixation of the fragments using an osteosynthesis device, mainly the intramedullary nail. This repairing method implies several complications, which may include the failure of the fixation device, frequently occurring due to the "cut-out" mechanism. The present work focuses on the study of how the position of the cephalic screw, which should be fixed during surgery, affects the cut-out risk. Through experimental tests and numerical models some variables that can be critical for the cut-out phenomenon are analysed.

METHODS

This study has been carried out through a numerical model based on the finite element method and experimental tests. The digital image correlation technique has been used in experimental tests to measure displacements on the femoral surface with the objective of numerical model validation. Some basic daily activities with different intramedullary nail positions have been analysed through the numerical model, considering variables that can induce the cut-out complication.

RESULTS

The results show how the intramedullary nail position clearly influences the cut-out risk, showing that displacements in the upper, anterior and posterior direction increase the cut-out risk, while displacement in the lower direction endangers the intramedullary nail itself. Thus, the centred position is the one which reduces the cut-out risk.

CONCLUSIONS

This work supposes an improvement in the knowledge of the cut-out phenomenon thanks to the combination of experimental testing and validated numerical models. The effects of different intramedullary nail positions in the femoral head are studied, including a novelty variable as torque, which is critical for the structural integrity of the fixation. The main conclusion of the work is the determination of the central intramedullary nail position as the most favourable one for decreasing the cut-out risk.

Fast NURBS-based parametric modeling of human calves with high-accuracy for personalized design of graduated compression stockings

BACKGROUND AND OBJECTIVES

Accurate human body models are increasingly demanded by high-quality human-centered ergonomic applications, especially the design and manufacturing of compressive functional apparels. However, existing parametric models in related works are not capable to accurately describe detailed local shape features of human.

METHODS

In this work, a high-accuracy parametric modeling approach for human limb was proposed. 3D Scans of human calves were studied. Key data points of the scanned human calves were identified according to human anatomy, forming a quasi-triangular mesh of feature points. Then, non-uniform rational B-splines (NURBS) method was implemented. Control points were calculated from the key data points, with which the human calf shapes can be reconstructed by the smooth NURBS surface, giving rise to a new parametric model of human calves. Error between the scanned and reconstructed calf shapes were analyzed to verify the effectiveness of this model.

RESULTS

Error analysis showed that, this proposed method delivers a high-efficiency and high-accuracy parametric shape modeling approach with averaged error observed as only 0.37% for all the 260 subjects, much less compared to previous relative works (around 5%). For tentative application, customized medical compression stockings were designed based on this model and proved as valid to exert desired gradient compression on the according calf mannequin.

CONCLUSIONS

By introducing the non-uniform rational B-splines method, a parametric model capable of characterizing human limbs with high-accuracy was proposed. Using very small amount of data, this model is expected to highly facilitate remote customized design and provide 3D shape references for design of compressive garments. Moreover, the proposed methods can inspire developments of other mixed modeling methods for high-accuracy applications.

Investigation on the Differences in the Failure Processes of the Cortical Bone under Different Loading Conditions

Cortical bone is a transversely isotropic material, and the mechanical properties may be related to the loading direction on the osteon. Therefore, analyzing the differences in the failure processes of cortical bone under different loading conditions is necessary to explore the measures for reducing the incidence of fracture. In this study, to investigate the effects of different loading directions on the fracture performance in the cortical bone, a numerical method that could simultaneously simulate the failure processes in the cortical bone structure under compression and bending loads was established based on continuum damage mechanics theory. The prediction accuracy and feasibility of the numerical method were first verified by comparing with the corresponding experimental results. Then, the differences in the failure process and fracture performance of the same cortical bone structure under compression and bending loads were investigated. The simulation results indicated that for the same structure, the slip-open failure mode appeared under compression load, and the crack propagated along a certain angle to the loading direction; the tension-open failure mode appeared under bending load, and the crack propagated along the direction perpendicular to the loading direction. Meanwhile, the fracture load was greater and the fracture time was later in the compression than in the bending condition. These phenomena stated that discrepant failure processes and fracture patterns occurred in the same cortical bone structure under different loading conditions. The main reason may be related to the tension–compression asymmetry and transversely isotropic characteristics in the cortical bone material. The fracture simulations in the cortical bone under different loading conditions could improve the prediction accuracy in bone biomechanics and provide the prevention method for cortical bone damage and fracture.