Modelling and Finite Element Analysis of Stress Distribution in Total Hip Replacement Using Different Material Combinations
Modelling and Finite Element Analysis of Stress Distribution in Total Hip Replacement Using Different Material Combinations
Modelling and Finite Element Analysis of Stress Distribution in Total Hip Replacement Using Different Material Combinations
net/publication/294407103
CITATIONS READS
0 732
5 authors, including:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Saravana Kumar Gurunathan on 14 February 2016.
Abstract:
The hip joint is one of the largest weight bearing structures in the human body. Depending on the type
of activity, the hip joint is subjected to load intensity fluctuating between zero and maximum of about 7
times the body weight. Because of this loading, mechanical failures of the hip joint and of the
neighboring bony structure, particularly the upper femoral region, constitute a large proportion of the
problems confronting the orthopedic surgeons. Numerical modeling techniques, such as the finite
element method, open the possibility of predicting hip joint biomechanics noninvasively and to
estimate the magnitude of the stresses in hip joint.
The paper discusses development of a CAD model for the artificial hip stem. The steps involved in this complex model
includes 3D laser scanning of the hip stem, repairing and cleaning the model using SolidWorks, meshing the model using
HyperMesh and finally importing the model to MSC.Marc for finite element analysis. The stresses in the hip stem were
analyzed using the commercial software MSC.Marc for various material combinations such as Titanium alloys and Co-Cr
alloys. The paper also discusses development of an integrated model where the artificial hip stem is inserted in femur and
the finite element analysis is carried out to study the stress distribution in the impalnt and bone. The importance of
mathematical tools to understand the complex biomechanics of the hip stem are highlighted.
Keywords: Total Hip Replacement, Reverse engineering, 3D scanning, Meshing, Finite Element Analysis.
M&MT 2011 01
Common causes of implant failure can be broadly
Acetabulum categorized as follows:
Cancellous
Head a) Aseptic loosening
b) Septic loosening
M&MT 2011 02
the implant. understanding of the biomechanical function of the hip
joint prosthesis and stress distrbution within bone and
2. Biocompatibility: Implant materials are expected to be
implant for different material of the implant is essential.
highly non toxic, non allergic, and non inflammatory.
Finite element method, one of the advanced simulation
3. High corrosion & wear resistance: Implant is expected
techniques has been used in orthopaedic biomechanics
not to release any toxic or allergic elements during
for many decades. It is an important tool used in the
debris formation over the years of usage.
design and analysis of total joint replacement and other
4. Osseo-integration: It is the process of implant surface orthopaedic devices. Xio - Wen Zheng developed a 3D
material integrating with adjacent bones & tissues. A finite element model to study stress distribution of the
layer of hydroxyl-apatite on the surface of implant is artificial femoral stem in a different cemented conditions
commonly used for this purpose which promotes and loading conditions [5]. El- Sheikh et al. studied the
fibrous growth of adjacent tissues resulting in an effect of prosthesis Young's modulus and cement Young's
improved fix. Inability to osseo-integration results in modulus on the stresses in the prosthesis stem and
loosening of implant. cement layer [6]. Pyburn et al. carried out FE analysis of
the hip implant for different cross sections of the stem [7].
David et al. carried out finite element analysis of different
hip stem designs. Subsequently, the stem that had
highest stress and displacement were optimized for
stress and displacement [8]. Oguz Kayabasi et al studied
the effect of static, dynamic and fatigue behaviour on
three dimensional shape optimization of hip prosthesis by
finite lement method [9]. Taylor et. al investigated
cancellous bone stresses surrounding femoral
component using finite element analysis [10]. The
resulting cancellous bone stresses were compared to that
of intact femur. Stress shielding of the femur is known to
be a principal factor in aseptic loosening of hip
replacements. Gross and Abel carried out finite element
Fig. 3. Elastic modulus of various biomedical alloys [4]
analysis of a hollow stemmed hip implant for reducing the
effects of stress shielding, while maintaining acceptably
Conventionally, Stainless Steel had been used with
low levels of stress in the cement [11]. Using finite element
polyethelene acetabular cup until a few decades ago
modelling, the stresses in the proximal femur using
when Ti6Al4V, a titanium alloy was being deployed for the
different shapes of hollow stem were compared with those
same purpose. Co-Cr alloys, on the other hand, have
produced using comparable sizes of solid stem with
been used recently for their very good wear resistant
different values of elastic modulus. Results showed that a
properties to make acetabular parts and femoral head.
reduction in stress shielding could be achieved with a
Thus, making a metal-on-metal joint contact. Different
hollow stem [11]. Oguz Kayabasi et al carried out static
materials with their elastic modulus used for impalnt are
and dynamic analysis of hip prosthesis for four different
shown in Fig. 3. Average service life of artificial hip joint is
stem shapes to investigate an optimum stem shape [12].
around 10 years. Middle aged patients and especially
Based on static and dynamic FE analysis results, safety
younger patients need to have arthoplasty replaced two
factors for fatigue life of the hip prosthesis were calculated
or more times within their expected life. A wide range of
[12]. Though extensive FE work has been carried out by
materials are available which may be suitable for use in
various researchers, using combination of materials for
artificial hip joint prosthesis. Young's modulus of
stem and neck were not attempted. In the present study,
prosthesis material is a critical design variable because
various Titanium alloys have been taken into
for a given stem geometry design, it largely determines
consideration, and their combination is analyzed for
how the load is transferred to the bone. So, thorough
M&MT 2011 03
stress distribution in implant & bone using finite element
method.
Methodology:
This paper describes the entire process of reverse
engineering of hip implant, followed by developing an Fig. 4 (a) Prototype
integrated CAD model of femur and implant, further
undergoing the process of meshing to be analyzed using
finite element analysis. The novel approach of work
focuses on applying different material properties for
different parts of implant in order to achieve an improved
stress distribution.
Reverse Engineering
The first part of the work deals with developing a CAD
model for an implant prototype available at Sandvik
Materials Technology R&D exhibition, Pune. Figure 4(a)
shows the hip stem. The hip stem is scanned using 3D
laser scanning facility at the Department of Engineering
Design, IIT Madras, at Chennai. This scanning was done
using a point scanner. But the point cloud data obtained
does not meet the desired quality. Figure 4(b) shows the
output from this scanning. Since the scanner was such
that laser emitting and receiving source was stationary,
and the object was kept on a rotary indexing table, a few
surfaces have not been traced by the laser at all. Also,
because of a complicated geometric shape of the
prototype, it was difficult to fix the object in desired
orientations. Therefore the scanning was again
performed using Perceptron ScanWorks manual
scanning system which uses the same principle of
triangulation with a resolution of 0.1mm. It is a hand held
line scanner which provides better flexibility in capturing
the shape and geometry of the implant. Since the scanner
is hand held, it provides enhanced flexibility to capture the Fig. 4(b) Scanned output from IIT Madras facility
geometry even at corners. Point cloud data is initially
extrapolated to form continuous geometry. The output appears good but not perfect. So, actual measurements
from second scanning is shown in Fig. 4(c). Though it were taken from the physical prototype and SolidWorks
M&MT 2011 04
The final CAD model of implant is integrated with 3D solid
model of left thigh femur downloaded from the repository
available in the website [13] which has been widely
referred in the literature. Inputs taken from Dr. Dastane at
Aditya Birla Memorial Hospital, Pune & Dr. Sancheti,
Sancheti Hospital, Pune provided information needed for
assembling femur & implant. This work considers
cementless hip replacement assembly. In a real scenario,
broaching is used to rasp the cancellous bone so that the
implant fits exactly into the void created. In 3D, it is
ensured that the implant lie as close to the centre at
different cross sections of the bone. A cut is made through
the neck of femur parallel to the shoulder of implant. This
process is called Osteoctomy during the real operation.
Figure 5 shows the final model.
M&MT 2011 05
this purpose. A large part of assembly is uniaxial and
relatively more uni-form than the rest. Thus, this part of the
assembly is converted into 3D solid 'hexahedron'
elements each with eight nodes. The process starts by
choosing one cross section, with creating a 2D mesh of
quad elements and refining them with the available tools.
This set of elements is duplicated and translated at the
other end of the part. The 2D planar elements are refined
again and projected as close to the geometric cross
section. Reference points and vectors are chosen at two
sets of element lie as close to the vertical axis, and 3D
hexahedron elements are created in the vertical direction 40
M&MT 2011 06
material properties were assumed to be homogeneous, Boundary conditions (Static)
isotropic and perfectly elastic except cancellous which Femur load
was assiumed to be isotropic and elastic –plastic and Both forces at 20
degrees from
cortical bone was assiumed to be orthotropic and elastic
vertical
–plastic.
Boundary Conditions
Abductor
Two different load cases were considered. In the first load muscle force
Y ~1.25kN
case, stress analysis of implant alone with different
X
material combination is carried out. A force of 3kN, 4KN
and 5KN is applied as a point load on the nodes over an Z
area on the implant head at an angle of 20 degrees from Fixed
the vertical. This is due to the fact that if force is applied on 18
Fig. 7. Boundary conditions for implant alone Fig. 9. Von Mises streses in neck of the implantg
M&MT 2011 07
observed in the neck when stem is made of lowest
stiffness material i.e. TiNbZr and neck is made of highest
stiffness material i.e. Ti67.
Fig. 10. von Mises stress contours on the medial side of Another interesting result is that for a given load, same
the implant peak Von Mises are generated when same material is
used for neck and stem but strain in higher stiffness
Effect of different material combinations for implant on the material are lower. Similar results were obtained when
stress distribution within the implant when 1/3 part of the load varied to 4 and 5 KN and full length of the stem was
stem is fixed is shown in Table 2. fixed.
Table 2: Peak von Mises stresses in the neck of the implant Further analyses are done with the integrated bone-
1/3 part is fixed implant assembly model. It is to be noted that the abductor
von Mises stresses (MPa) at 3KN muscle force of magnitude ~1.25kN was considered
Load acting on upper trochanter in the direction almost exactly
Neck\Stem TiNbZr Ti64 TiCP Ti67
opposite to the femoral load. Load is varied from 3kN to
5kN. Stress contours for different material combination
TiNbZr 335.2 316.5 314.3 313.1 are plotted and results are tabulated.
Ti64 359.8 335.2 332.5 331.1
TiCP 364.9 340.2 337.4 336
Ti67 366.2 339.5 336.6 335.2 Comparison of peak stresses (neck_stem)
356Mpa- TiNbZr_TiCP
From the Table 2, it is observed that for a given neck 378Mpa- Ti64_Ti64
material when stiffness of the stem increases, peak Von 405Mpa- Ti67_TiNbZr
the neck. Similarly for a given stem material when Eq. Von Mises Stresses
M&MT 2011 08
TiNbZr(70)
Ti64(110)
5kN
4kN
TiCP(114) 3kN
Ti67(120)
0 10 20 30 40 50 60 70
TiNbZr(70)
Ti64(110)
5kN
4kN
TiCP(114) 3kN
Ti67(120)
Material neck and stem with lower stiffness and neck higher
for Neck stiffness generate higher stiffness. Increasing neck
stiffness for a given Ti64 as a stem material, stresses in
neck increases.
M&MT 2011 09
stem is resting on cancellous bone. Next higher stresses 2. Manish Kadam - Research Manager, Sandvik
were observed on the medial side at proximal end. Material Technology R&D, Sandvik Asia Pvt Ltd.
Graphs in Fig. 11(c) indicate there is little variation with 3. Dr. Palla Sivaprasad, General Manager, Sandvik
different implant material combinations Materials Technology R&D, Sandvik Asia Pvt. Ltd
Detailed study of proper selection of prosthesis material References
from point of view of stress distrbution in the bone will
1. William J. Hozack, Richard H. Rothman, Robert E.
form the basis for our future work. Booth Jr., Richard A. Balderston, Joy C. Cohn, Gary
Conclusions T. Pickens, The Journal of Arthroplasty, Volume 5,
Issue 1, March 1990, pp. 41-47,
From the peak stress values obtained from implant alone
2. In Gwun Jang, Il Yong Kim, Journal of
analysis and complete THR assembly analysis, following
Biomechanics, Volume 41, Issue 11, 7 August 2008,
conclusions are drawn: pp. 2353-2361.
1. Von Mises stresses of the artificial femoral stem are 3. D. R. Sumner, T. M. Turner, R. Igloria, R. M. Urban,
not significantly affected by cemented condition. J. O. Galante Journal of Biomechanics, Vol 31,
2. Neck peak stress is sensitive to elastic modulus. Issue 10, pp. 909-917,
3. Increasing neck stiffness increases neck stress for a 4. M. Geetha, A.K. Singh, R. Asokamani, A.K. Gogia,
given stem material. Progress in Materials Science, Volume 54, Issue 3,
4. Increasing stem stiffness decreases neck stress for a May 2009, pp. 397-425,
given neck material 5. Xiao-Wen Zheng in Bioinformatics and biomedical
5. Peak stresses in the neck do not vary when same engineering, 2009 ICBBE 3Rd international
material is used for neck and stem. conference
6. Neck with lower elastic modulus combined with stem 6. H.F.El-Sheikh, B.J. Macdonald M.S.J.Hashmi,
of higher elastic modulus turns out to be a better Journal of material processing technlogy,122,2002,
combination in the case of static loading, as far as pp. 309-317,
peak stresses in prosthesis neck are concerned. 7. E.Pyburn, T.Goswami, Materials and Design, 25,
7. There is no significant change in stress distribution in 2004, pp. 705-713.
cortical and cancellous bones with different material 8. David Bennett, Tarun Goswami, Materials &
combinations. Design, Volume 29, Issue 1, 2008, pp. 45-60,.
8. In general, Ti-29Nb-13Ta-4.5Zr appears as one of 9. Oguz kayabasi and Bulent Ekici, Materials and
the most promising alloys for neck providing better Design, Vol. 28, 2007, pp 2269-2277.
solutions ahead. 10. M.Taylor, K.E. Tanner, M.A.R. Freeman, A.L.
Yettram, Med eng Physics Vol 17, No7, pp. 544-
Acknowledgement:
550, 1995.
We are grateful to Perceptron Non-contact Metrological 11. A finite element analysis of hollow stemmed hip
Solutions, Chennai and Sandvik Mining & Construction prosthesis as a means of reducing stress shielding
for providing their support, facilities and guidance of the femur, Journal of Biomechanics, 34, 2001 pp.
through the course of this project. We are also thankful 995–1003, S Gross, E.W.Abel.
to Dr. Manish Dastane, Aditya Birla Memorial hospital 12. Oguz Kayabasi, Fehmi Erzincanli, Advances in
Pune, for his valuable inputs. Engineering Software, Volume 37, Issue 7, July
2006, pp. 477-483.
Notes about Authors:
13. Marco Viceconti, Massimiliano Casali, Barbara
1. Manu Gupta is a final year dual degree (BTech & Massari, Luca Cristofolini, Sanzio Bassini, Aldo
MTech) engineering student at the Department of Toni, Journal of Biomechanics, vol. 20, No. 9, pp.
Engineering Design, Indian Institute of Technology, 1241 (1996); www.biomedtown.org
Chennai, India.
n
M&MT 2011 10