Finite Element Modelling and Analysis of Trans-Tibial Prosthetic Socket

Table of contents

1. Introduction

he socket is a basic component for prosthetic performance. Below-knee amputees generally demonstrate some gait abnormalities such as lower walking speed [1], incresed energy cost [2], and asymmetries between legs of unilateral amputees in stance phase cycle, step length and maximum vertical force [3]. Successful fitment of prosthesis may be achieved by understanding the biomechanical structure of socket and its material, weight, thickness in particular to fulfill the desirable load distribution in soft tissues and bone of residual limb. Most commonly used socket design in developing countries is pattelar tendon bearing (PTB) socket developed following the World War II at the University of California, Berkeley in the late 1950 s [4,5]. The Finite Element Method (FEM) has been used widely in biomechanics to obtain stress, strain and deformation in complicated systems and have been identified as an important tool in analysing load transfer in prosthesis [6]. The finite element analysis (FEA) models have been used to study the effects of the inertial loads and contact conditions on the interface between prosthetic socket and stump of an amputee Author ? ? ?: M Tech, Mechanical Engineering Dept., ISM Dhanbad, India. e-mails: [email protected], [email protected], [email protected] during gait cycle [7,8].The finite element methode has been used as a tool for parametric study and evaluation of prosthetic socket [9,10].

It is common for amputees to experience pain and discomfort in the residual limb while wearing the prosthetic socket [11]. For a lower limb amputee, the comfortableness of wearing prosthesis depends on the distribution of stress at the interface of residual limb and prosthetic socket is either at the pressure-tolerant (PT) or pressure-relief (PR) areas. By employing the technology of computer-aided engineering, the quality uncertainty and labour intensity of traditional process of fabricating a prosthetic socket can be improved. Lower limb prosthesis allows ambulation and improves the performance of daily routine activities. However, poorfitted socket can lead to complications that have adverse effects on the activity level and gait cycle of people with lower limb amputation [12].

The interface between the stump of lower limb amputees and their prostheses is the prosthetic socket. The contact pressure at the residual limb and prosthetic socket interface is an essential index, and is considered as a promising measure towards good socket design. Therefore, the fundamental concern is to understand pressure distribution at the stump-socket interface. Although the use of pressure sensor is a direct experimental approach towards estimating interface pressure, the analytical approach is an alternating to the experimental one, and finite element modeling of the socket has been used to analyse the contact pressure. Although, the complex features of the soft limb tissues and of their interaction with the socket still remains difficult to model [13].

The variation of interface pressure between the stump and socket is an important factor in socket design and fit. Lower limb prosthetic socket users experience pressure between the stump and socket during daily routine activities. The underlying soft tissues and skin of the stump are not habitual to weight bearing; thus, there is the risk of degenerative tissue ulcer in the stump because of cyclic or constant peak pressure applied by the prosthetic socket [14]. The pressure also can lead to various skin deases such as follicular hyperkeratosis, allergic contact dermatitis, infection and veracious hyperplasia [15][16][17].

Despite significant scrutiny in the field of prosthetics in the previous decades, still many amputees experience pressure ulcers with the use of prostheses. Sometimes, skin problems lead to chronic infection, which may necessitate re-amputation. This will obviate the long-term use of prosthesis, which indicatively reduces the routine activities of prosthesis users and the quality of life [18]. Many studies have concentrated on interface pressure magnitude between the socket and stump during level walking [19][20].

2. a) Trans-Tibialprosthesis Description

The artificial limb consists of a foot-ankle unit which needs to be attached to the remainder of the amputee's natural leg or stump. The foot ankle unit is attached directly to the socket frame. The artificial shank can be attached to the foot ankle unit and then attached to the socket frame for a below-knee amputation. Today the sockets are roughly quadrilateral in shape. They attempt to have total contact between the stump and the socket.

3. II.

4. Finite Element Model

A frequently used numerical analysis technique in biomechanics is the finite element technique, a computational approach for interface stress or structural deformation calculation evoluted in engineering mechanics. It has been introduces as a useful tool to understand the load transfer mechanics between a residual limb and its prosthetic socket. The finite element technique is a full-field analysis for calculating the state of stress and elestic strain in the specific field. This technique is well suited for parametric analysis in the process of design. The previous finite element analyses showed the significance of considering prestress in predicting interface stresses at loading stage [21][22].

One left unilateral male trans-tibial amputee participated in this study. The volunteer was 45 years of age, 166 cm tall, 70 kg in mass, and the cause of his amputation was an accident. He has been an active amputee for five years, using his prosthetic limb for all his daily chores. The simplified geometry of his residual limb was modeled in Pro-engineer and then it is being imported (in IGES formate) and modified in ANSYS 14.0 Workbench.

The finite element model was kept as simple as possible in terms of material properties and boundary conditions. Different materials for 2 mm, 3 mm, 4 mm, 5 mm, and 6 mm unit volume layer thickness was used for creating the 3-D FE model. Also the three dimension finite element model is developed for varrying the length of the socket of 16 cm, 17 cm, 18 cm, 19 cm, 20 cm, 21 cm. The model was meshed with brick element solid 185 with fused tibia and fibula bones. A total of 41,073 elements and 20,438 nodes were used. IV.

5. Results

The results for total deformation, shear stress and equivalent von-Mises stress of developed and modified transtibial socket model were obtained by using (ANSYS Workbench v14.0) program. Figure (1) shows the meshed view of the finite element three dimension socket model, and figure (2) shows the maximum von-Mises stress developed for different loading conditions of stance phase like Initial On the meshed model fixed support is being applied at the distal end of the socket, distal end of the socket is further attached with the remaining parts of the prosthesis like shank, ankle foot. The different loading conditions as listed in table 2 were quasi-static approximations using experimentally obtained maximum vertical ground reaction for the prosthetic side of same subject while walking at a given speed using CGD gait cycle analyzer [24][25][26][27].

6. III.

7. Material Properties

In this analysis the different material used are composite, polypropylene, 90/10 PP/PE (90% polypropylene and 10% ethylene), high density polyethylene (HDPE) and low density polyethylene (LDPE). The mechanical properties of the socket material were assumed to be linearly elastic, isotropic and homogeneous. Socket were analyzed for different materials and their valus of Young's modulus, Ultimate strength, Poisson's ratio and density is listed in table 1.

8. Tsai-Hill Criterion, C TH

= ?? 1 2 ?? ?? 2 ? ?? 1 * ?? 2 ?? ?? 2 + ?? 2 2 ?? ?? 2 + ?? 2 ?? ??? 2

Where CTH is the Tsai-Hill failure coefficient, Sv, St and Ssh are the ultimate strengths of composite in the vertical, transverse and shear directions respectively listed in table 5 and ?1, ?2 and ? are the imposed stresses in the longitudinal, transverse, and shear planes. If the value of CTH is less than one than design is safe. The values of maximum von-Mises stress, shear stress and total deformation off all the material in different length were analyzed and shown in figures 11-13, and it is found that as the length of socket increses the values of stress and deformation decreases. The decrease in value of deformation as increase of length is higher in case of LDPE material.

9. d) Structural Behavior vs Thickness

The values of maximum von-Mises stress, shear stress and total deformation off all the material in different thickness were analyzed and shown in figures 8-10, and it is found that as the thickness of socket increses the values of stress and deformation decreases. The decrease in value of deformation as increase of thickness is higher in case of LDPE material.

10. Strength (in MPa)

11. Sv

St Ssh The variation of factor of safety as a function of Weight of the socket for a socket of different thickness is shown in figure 7, where factore of safety is being calculeted by dividing maximum von-Mises stress at a load of 620 N with the endurance limit (50% ultimate tensile strength value) [28]. During daily activities of an amputee the total load of knee joint in transtibial prosthesis passes on the prosthetic socket. During normal walking , the total joint reaction forces at knee joint is three to four times increases than the total body weight, during jumping and fast running load on knee joint increses more [29]. Therefore, six factore of safety is minimum desirable to withstand the loading of socket. The factor of safety is just below the level of five for LDPE and HDPE so, it can be suggested that LDPE and HDPE are note suitable for prosthetic socket design.

12. b) Case 2 : Analysis of failure

The finite element simulation result of rotation and displacement in different parts of socket validate the biomechanical requirement of structural integrity in patellar tendon bearing socket. Figure 7 shown below describes the variation of Tsai-Hill coefficient with tensile and compressive strength. The value of CTH coefficient in 2mm thick composite for tensile strength is 0.1864wich is only five times factore of safety but thickness between 3 mm (0.0724) to 4 mm (0.031) has a factore of safety more than twenty times. Therefore, the optimum solution of composite material of thickness 3 mm to 4 mm satisfied the Tsai-Hill criterion. In all materials it is found that the von-Mises stress, shear stress and total deformation is inversely proportional to thickness except for LDPE of the socket figure 8-10. However the stress and stress variation were higher in case of 2mm and 3 mm socket and it is relatively low in case of 4mm and 6mm. Thus 3 mm to 4 mm could be a optimal solution in terms of thickness of the socket for all materials where this much thickness is used. The variation of von-Mises stress, shear stress and total deformation for different thickness of prosthetic socket were shown in figure 8, 9 and 10 respectively. The value of total deformation in case of LDPE of thickness less than 3 mm goes higher and it may loss biomechanical load bearing ability. Thus the result indicates that the LDPE socket length is not suitable fore fabrication of PTB socket of below 4 mm thickness. In all materials it is found that the von-Mises stress, shear stress and total deformation is inversely proportional to length of the socket figure 11-13. However the stress and stress variation were higher in case of 16 and 17 cm length socket and it is relatively low in case of 19cm and 20cm. Thus 19 cm to 20cm could be a viable solution in terms of length of the socket for all materials where this much length is possible. The variation of von-Mises stress, shear stress and total deformation for different length of prosthetic socket were shown in figure 11, 12 and 13 respectively. The value of total deformation in case of LDPE of length less than 16 cm goes higher and it may loss biomechanical load bearing ability. Thus the result indicates that the LDPE socket length is not suitable fore fabrication of PTB socket of below 16cm length.

13. Conclusions

The results summarized that assimilating local submissive properties within socket wall can be an effective methods to distribute maximum stress areas and also to relief contact pressure between the socket and stump. Based on the results and the discussion, the composite material is cheap, excellent strength, widely available butit has high weight that make it only useful to be used for adult with higher weights. The results obtained from analysis can be used as a reference to choose socket material, thickness and its optimal length for manufacturing of socket in developing countries. The socket buildup of composite material gives the optimal solution for patellar-tendon bearing socket design. The study reconnoitered further future scope for parametric analysis, investigating the effects of liner, socket stiffness, rectification scheme, soft tissues, and materials for the socket/stump interface stress distribution. Review of secondary physical conditions associated with lower-limb amputation and long-term prosthesis use.

Figure 1.
Contact/Heel Strike (HS), Foot Flat/Loading Response (LR), Mid-Stance (MS), Terminal Stance/Heel Off (HO) and Pre-Swing/Toe Off (TO).
Figure 2. Figure 1 :Figure 3 :
13Figure 1 : Meshed view Figure 2 : Equivalent von-Mises stress at different stance Figure 3 : shows the shear stress and figure (4) shows the total deformation developed for different loading conditions like Initial Contact/Heel Strike (HS), Foot Flat/Loading Response (LR), Mid-Stance (MS), Terminal Stance/Heel Off (HO) and Pre-Swing/Toe Off (TO) of stance phaseThe maximum values of von-Mises stress, shear stress and total deformation occure in five phases of stance under the load of 1000 N are listed in the table 3 as shown below.
Figure 3. Figure 3 :
3Figure 3 : Shear Stress
Figure 4. Figure 4 :
4Figure 4 : Total deformation Different approaches have been used to analyzed for socket material and thickness optimization: a) Socket thickness vs Factor of safety For different thickness of composite, polypropylene (PP), 90/10 PP/PE, HDPE and LDPE materials the weight of the socket were calculated and listed in table 4. b) Tsai-Hill Criterion for socket failure Socket failure is analyzed by Tsai-Hill Criterion based on maximum distortion criterion and used in this analysis to compare socket failure [23].
Figure 5. Figure 5 :
5Figure 5 : Equivalent von-Mises stress at pressure tolerated areas in stance phase of gait cycle during normal walking on plane surface at a given speed
Figure 6. Table 1 :
1
Figure 7. Table 2 :
2
Year 2014
44
( ) A Vol ume XIV Issue IV Version I
Global Journal of Researches in Engineering
Figure 8. Table 3 :
3
Stance phase von- Shear Total
Mises stress deformation
stress (MPa) (mm)
(MPa)
Initial 6.55 1.02 0.55
Contact/Heel
Strike (HS)
Foot Flat/Loading 10.54 1.64 0.88
Response (LR)
Mid-Stance (MS) 7.45 1.16 0.63
Terminal 10.54 1.64 0.88
Stance/Heel Off
(HO)
Pre-Swing/Toe Off 8.54 1.33 0.72
(TO) of stance
phase
Figure 9. Table 4 :
4
Thickness Composite PP PP/P HDP LDP
(mm) E E E
2 140 106 102 111 107
3 209 160 155 168 174
4 280 212 208 214 222
5 350 266 261 278 269
6 420 320 315 334 323
Table 5 : Tensile and compressive strength of
composite [31]
c) Structural Behavior vs Length
Figure 10.
Tension Compression 584 803 43 187 Year 2014 44
ume XIV Issue IV Version I
1

Appendix A

  1. , ProsthetOrthot Int 1999. 23 p. .
  2. Clinical investigation of the interface pressure in the trans-tibial socket with Dermo and Seal-In X5 liner during walking and their effect on patient satisfaction. Abu Alis , N A Osman , N Mortaza , A Eshraghi , H Gholizadeh , Wan Abas , Wab . Clinical Biomechanics 2012. 27 p. .
  3. The patellar-tendonbearing below knee prosthesis. C W Radcliff , J Foort . Biomechanics laboratory 1961. University of California
  4. Experimental and numerical predictions of the ultimete strength of a low-cost composite transtibial prosthesis. J Hahl , M Taya . Journal of Rehabilitation Research and Development 2000. 37 p. .
  5. Biomechanical basis of human movement. 2 nd Edition William & Wilkins Publisher. J Hamill , Knutzen . Philadelphia 1961. 237.
  6. Accelerographic, temporal, and distance gait: factors in below-knee amputees. J L Robinson , G L Smidt , J S Arora . Phys Ther 1977. 57 p. .
  7. Changes in EMG power spectrum (high-to-low ratio) with force fatigue in humans. J Moxham , Rht Edwards , M Aubier , A Detroyer , G Farkas , P T Macklem , C Roussos . Journal of Applied Physiology 1982. 53 p. .
  8. Skin disorders in amputees. Lyoncc , J Kulkarni , E Zimerson , Van Rosse , M H Beck . J Am Acad Dermatol 2000. 42 p. .
  9. Parametric analysis using the finite element method to investigate prosthetic interface stresses for persons with trans-tibial amputation. M B Silver-Thorn , D S Childress . J Rehabil Res Dev 1996. 33 p. .
  10. Prediction and experimental verification of residual limb/prosthetic socket interface pressure for below knee amputees, M B Silver-Thron . 1991. Evanstoon. Northwestern University
  11. The quasistatic response of compliant prosthetic sockets for transtibial amputees using finite element methods. M C Faustini , R R Neptune , R H Crawford . Medical Engineering & Physics 2006. 28 p. .
  12. An iterative method for viscoelastic modeling of prosthetic feet. M D Geil . Journal of Biomechanics 2002. 35 p. .
  13. VerstraeteMC and Njus GO. Finite element analysis as a tool for parametric prosthetic foot design and evaluation. Technique development in the solid ankle cushioned heel (SACH). M M Saunders , E P Schwentker , D B Kay , G Bennett , C R Jacobs . Computer Methods in Biomechanics and Biomedical 2003. 6 p. .
  14. Energy/speed relation of below-knee amputees walking on a motor-driven treadmill. N H Molen . Int Z Angew Physiol 1973. 31 p. .
  15. Dermatologic conditions associated with use of lower-extremity prosthesis. N L Dudek , M B Marks , S C Marshall , J P Chardon . Arch Phys Med Rehabil 2005. 86 p. .
  16. Ambulation monitoring of transtibial amputation subjects with patient activity monitorversus pedometer. N L Dudek , O D Khan , E D Lemaire , M B Marks , L Saville , N Maffiulletti . J Rehabil Res 2008. 45 p. .
  17. Socket/stump interface dynamic pressure distributions recorded during the prosthetic stance phase of gait of a trans-tibial amputee wearing a hydrocast socket, P Convery , A Buis .
  18. Timing and relative intensity of hip extensor and abductor muscle action during level and stair ambulation: an EMG Study. Perry J Lyonsk , J Gronley , L Barnes , D Antonelli . Journal of the American Physical Therapy Association 1983. 63 p. .
  19. Analysis of transtibial prosthetic socket materials using finite element method. P K Lenka , A R Choudhury . J. Biomedical Science and Engineering 2011. 4 p. .
  20. Radcliff CW Functional considerations in the fitting of above-knee prostheses. Artificial Limb 1955. 2 p. .
  21. Phasic activity of intrinsic muscles of the foot. R A Mann , V T Inman . Journal of Bone and Joint Surgery 1964. 46 p. .
  22. Principles of composite materials mechanics, R F Gibson . 1994. New York: McGraw-Hill.
  23. , R Gailey , K Allen , J Cstles , J Kucharic , M Roeder . 2014.
  24. Energy cost of walking of amputees: the influence of level of amputation. R L Waters , J Perry , D Aatonelli , H Hislop . J Bone Joint 1976. 58 p. .
  25. ume XIV Issue IV Version I interface of above-knee prosthesis. Advances in Engineering Software 1997. 28 p. .
  26. Electromyography of idiopathic toe walking. V Kalen , N Adler , E E Bleck . Journal of Pediatric Orthopedics 1986. 6 p. .
  27. Finite element modeling of the contact interface between trans-tibial residual limb and prosthetic socket. Wcc Lee , M Zhang , Jia , Jtm Cheung . Medical Engineering and Physics 1961. 26 p. .
  28. A computation model for monolimb design. Wcc Lee , M Zhang , X H Jia , D A Boone . Proceedings of the International Society of Biomechanics, (the International Society of BiomechanicsDunedin, New Zealand
    ) 2003. p. 234.
  29. FE modeling of the load transfer between trans-tibial residual limb and prosthetic socket. Wcc Lee , M Zhang , X H Jia , Jtm Cheung . Med Engg Phys 2004. 26 p. .
  30. Regional differences in pain threshold and tolerance ofthe transtibial residual limb: Including the effects of age and interface material. Wcc Lee , M Zhang , A F Mak . Arch Phys Med Rehabil 2005. 86 p. .
  31. Regional differences in pain threshold and tolerance of the transtibial residual limb: Including the effects of age and interface material. Wcc Lee , M Zhang , A F Mak . Arch Phys Med Rehabil 2005. 86 p. .
  32. Load transfer mechanics between trans-tibial prosthetic socket and residual limb-dynamic effects. X Jia , M Zhang , Lee Wcc . Journal of Biomechanics 2004. 37 p. .
Notes
1
© 2014 Global Journals Inc. (US)
Date: 2014-01-15