Eftaxiopoulos, Dimitrios

Stamou, Aikaterini

Trikkas, Panayiotis

English

Coventry University Pure Portal

FE analysis of the proximal femurFinite element analysis of the proximal femurP. Trikkas, A. Stamou and D. A. EftaxiopoulosSchool of Applied Mathematical and Physical SciencesNational Technical University of AthensGreeceFreeFem++ Workshop 2009FE analysis of the proximal femurSegmentationComputer tomography scans of the hip, from the VisibleWoman Project, were used.The open source medical image tool Slicer, was used forthe segmentation of the cortical and the trabecular part ofthe proximal femur.Automatic simple region growing and manualsegmentation were combined.Two .stl files were created, one for each part of the bone,after surface smoothing was performed. The surface nodeswere coincident on the interface between the two volumes.FE analysis of the proximal femurSlicer coronal viewFE analysis of the proximal femurThe two .stl filesFE analysis of the proximal femurMeshing, bonding and boundary condition surfacesThe open source FE discretization tool Gmsh was used formeshing the volumes of the cortical and the spongy bone.Around 400000 tetrahedra were created.The bonding (welding) of the two volumes was achievedvia a code in FreeFem++.New surface labels for the assignment of boundaryconditions were created within Gmsh.FE analysis of the proximal femurThe bonding codeload "msh3"load "tetgen"mesh3 Th1("floioskat.mesh");mesh3 Th2("spogkat.mesh");int[int] r1=[3,5], r2=[1,7];mesh3 Th3=tetgreconstruction(Th2,switch="rYYCVVV",refface=r2, reftet=r1);mesh3 Th=Th1+Th3;savemesh (Th,"miktokat.mesh");FE analysis of the proximal femurThe boundary condition surfacesFE analysis of the proximal femurThe boundary condition volumeFE analysis of the proximal femurThe material properties assignment codeint floios=reg(-141.642, -11.5074, 6.68015);Ep=600.+17000.*(region==floios);nu=0.3-0.02*(region==floios);FE analysis of the proximal femurThe Poisson’s ratio distributionFE analysis of the proximal femurThe applied forcesFE analysis of the proximal femurThe FE analysis codesolve Elasticity([u1,u2,u3],[v1,v2,v3],solver=CG)=int3d(Th)((Ep*nu/((1.+nu)*(1.-2.*nu)))*div(u1,u2,u3)*div(v1,v2,v3)+2.*(Ep/(2.*(1.+nu)))*(epsilon(u1,u2,u3)'*epsilon(v1,v2,v3)))-int3d(Th,55)(((-0.54/volume55)*836)*v1)-......-int2d(Th,70)(((0.58/area70)*836)*v1)FE analysis of the proximal femurThe displacement magnitudeFE analysis of the proximal femurThe maximum principal strainFE analysis of the proximal femurThe maximum shear strainFE analysis of the proximal femurThe maximum principal stressFE analysis of the proximal femurThe maximum shear stressFE analysis of the proximal femurThe Von Mises stressFE analysis of the proximal femurThe Von Mises stressFE analysis of the proximal femurThe strain energy densityFE analysis of the proximal femurConclusionsA linear finite element analysis of the proximal femur wasperformed.The maximum displacement magnitude occurs within thefemoral head where the hip contact force is applied.The maximum tensile fracture risk in the cortex seems tobe on the lateral part of the upper diaphysis and on theupper part of the femoral neck.

Finite Element Analysis of the proximal femur

https://pure.coventry.ac.uk/ws/files/30261087/Finite_Element_Analysis_of_the_proximal_femur..pdf

Finite Element Analysis of the proximal femur

Abstract

Similar works

Full text

Available Versions

Coventry University Pure Portal