Bone remodeling effects is a significant issue in predicting long term stability of hip arthroplasty. It has been frequently observed around the femoral components especially with the implantation of prosthesis stem. Presence of the stiffer materials into the femur has altering the stress distribution and induces changes in the architecture of the bone. Phenomenon of bone resorption and bone thickening are the common reaction in total hip arthroplasty (THA) which leading to stem loosening and instability. The objectives of this study are (i) to develop inhomogeneous model of lower limbs with hip osteoarthritis and THA and (ii) to predict the bone resorption behavior of lower limbs for both cases.

Biomechanical evaluations of lower limbs are established using the finite element method in predicting bone remodeling process. Lower limbs CT-based data of 79 years old female with hip osteoarthritis (OA) are used in constructing three dimensional inhomogenous models. The FE model of lower limbs was consisted of sacrum, left and right ilium and both femur shaft. Bond between cartilage, acetabulum and femoral head, sacrum and ilium were assumed to be rigidly connected. The inhomogeneous material properties of the bone are determined from the Hounsfield unit of the CT image using commercial biomedical software. A load case of 60kg body weight was considered and fixed at the distal cut of femoral shaft. For THA lower limbs model, the left femur which suffering for hip OA was cut off and implanted with prosthesis stem. THA implant is designed to be Titanium alloy and Alumina for stem and femoral ball, respectively.

Distribution of young modulus of cross-sectional inhomogeneous model is presented in Fig. 2 while model of THA lower limbs also shown in Fig. 2. Higher values of young modulus at the outer part indicate hard or cortical bone. Prediction of bone resorption is discussed with the respect of bone mineral density (BMD). Changes in BMD at initial age to 5 years projection were simulated for hip OA and THA lower limbs models.

The results show different pattern of stress distribution and bone mineral density between hip OA lower limbs and THA lower limbs. Stress is defined to be dominant at prosthesis stem while femur experienced less stress and leading to bone resorption. Projection for 5 years follow up shows that the density around the greater tronchanter appears to decrease significantly.

Effectiveness and long term stability of hip resurfacing and total hip arthroplasty for osteoarthritis patients are still debated nowadays. Several clinical and biomechanical issues have to be considered, including pain relief, return to function, femoral neck fractures, impingement and prosthesis loosening. Normally, patients with hip arthroplasties are facing gait adaptation and at risk of fall. Sudden impact loading and twisting during sideway falls may lead to femoral fractures and joint failures. The purposes of this study are (i) to investigate the stress behavior of hip resurfacing and total hip arthroplasty, and (ii) to predict pattern of femoral fractures during sideway falls and twisting configurations.

Computed tomography (CT) based images of a 54-year old male were used in developing a 3D femoral model. The femur model was designed to be inhomogeneous material as defined by Hounsfield Unit of the CT images. CAD data of hip arthroplasties were imported and aligned to represent RHA and THA femur modelas shown in Fig.1. Prosthesis stem is modeled as Ti-6Al-4V material while femoral ball as Alumina properties. Meanwhile, RHA implant is assigned as Co-Cr-Mo material. Four types of loading and boundary conditions were assigned to demonstrate different falling (FC) and twisting (TC) configurations (see Fig.2). Finite element analysis combined with a damage mechanics model was then performed to predict bone fractures in both arthroplasty models. Different loading magnitudes up to 4BW were applied to extrapolate the fracture patterns.

Prediction of femoral fracture for RHA and THA femurs are discussed in corresponding to maximum principal stress and damage formation criterion. The load bearing strain was set to 3000micron, the physiological bone loading that leads to bone formation. The test strength was wet to 80% of the yield strength determined from the CT images. Different locations of fracture are predicted in each configuration due to different loading direction and boundary conditions as shown in Fig.3. For falling configurations, fractures were projected at trochanteric region for intact and RHA femur, while THA femurs experience fracture at inner proximal region of bone. Differs to twisting configurations, both arthroplasties were predicted to fracture at the distal end of femurs.