A Three-Dimensional Finite Element Model for Biomechanical Analysis of the Hip

The objective of this study was to construct a three-dimensional (3D) finite element model of the hip. The images of the hip were obtained from Chinese visible human dataset. The hip model includes acetabular bone, cartilage, labrum, and bone. The cartilage of femoral head was constructed using the AutoCAD and Solidworks software. The hip model was imported into ABAQUS analysis system. The contact surface of the hip joint was meshed. To verify the model, the single leg peak force was loaded, and contact area of the cartilage and labrum of the hip and pressure distribution in these structures were observed. The constructed 3D hip model reflected the real hip anatomy. Further, this model reflected biomechanical behavior similar to previous studies. In conclusion, this 3D finite element hip model avoids the disadvantages of other construction methods, such as imprecision of cartilage construction and the absence of labrum. Further, it provides basic data critical for accurately modeling normal and abnormal loads, and the effects of abnormal loads on the hip.     

The influence of the acetabular labrum on hip joint cartilage consolidation: a poroelastic finite element model

The goal of this study was to investigate the influence of the acetabular labrum on the consolidation, and hence the solid matrix strains and stresses, of the cartilage layers of the hip joint. A plane-strain finite element model was developed, which represented a coronal slice through the acetabular and femoral cartilage layers and the acetabular labrum. Elements with poroelastic properties were used to account for the biphasic solid/fluid nature of the cartilage and labrum. The response of the joint over an extended period of loading (10,000s) was examined to simulate the nominal compressive load that the joint is subjected to throughout the day. The model demonstrated that the labrum adds an important resistance in the flow path of the fluid being expressed from the cartilage layers of the joint. Cartilage layer consolidation was up to 40% quicker in the absence of the labrum. Following removal of the labrum from the model, the solid-on-solid contact stresses between the femoral and acetabular cartilage layers were greatly increased (up to 92% higher), which would increase the friction between the joint surfaces. In the absence of the labrum, the centre of contact shifted towards the acetabular rim. Subsurface strains and stresses were much higher without the labrum, which could contribute to fatigue damage of the cartilage layers. Finally, the labrum provided some structural resistance to lateral motion of the femoral head within the acetabulum, enhancing joint stability and preserving joint congruity.     

Influence of cartilage cap modelling on FEM simulation of femoral head stress]

The present study reports on the finite element analysis (FEA) of the femoral head in a process of preparation for a program for the realistic simulation of correctional osteotomies of the proximal femur. While the material properties have been studied extensively, only few publications consider the influence of the cartilage layer geometry on FE stimulation of the hip joint. Various models of the femoral head with and without the cartilage layer are generated and analysed. On looking at the maximum surface stresses, we found a strong influence of the cartilage layer and the subchondral osseous layer on the magnitude of the von Mises equivalent stress. The model with an anatomically realistic cartilage layer and compact bone shows stresses of between 4 and 5.5 MPa, depending on the position of the joint, while the model with a concentric cartilage layer has a maximum von Mises stress of 0.8 MPa. Only on simulation of a "realistic" cartilage layer, with a maximum thickness at the "pole" and minimum thickness at the "equator" do the changes in stress distribution--determined by changes in the position of the femoral head--become visible. Owing to major artefacts and the inability to create a realistic cartilage layer, voxel-based models of the femur are not suitable for the simulation of the femoral head surface.     

Hemiarthroplasty of hip joint: An experimental validation using porcine acetabulum

Biphasic properties of articular cartilage allow it to be an excellent bearing material and have been studied through several simplified experiments as well as finite element modelling. However, three-dimensional biphasic finite element (FE) models of the whole joint are rare. The current study was carried out to experimentally validate FE methodology for modelling hemiarthroplasty. Material properties such as equilibrium elastic modulus and permeability of porcine acetabular cartilage were initially derived by curve-fitting an experimental deformation curve with that obtained using FE. These properties were then used in the hemiarthroplasty hip joint modelling. Each porcine acetabular cup was loaded with 400N using a 34mm diameter CoCr femoral head. A specimen-specific FE model of each acetabular cup was created using μCT and a series of software processes. Each model was analysed under conditions similar to those tested experimentally. Contact stresses and contact areas predicted by the model, immediately after loading, were then compared with the corresponding experimentally measured values. Very high peak contact stresses (maximum experimental: 14.09MPa) were recorded. A maximum difference of 12.42% was found in peak contact stresses. The corresponding error for contact area was 20.69%. Due to a fairly good agreement in predicted and measured values of contact stresses and contact areas, the integrated methodology developed in this study can be used as a basis for future work. In addition, FE predicted total fluid load support was around 80% immediately after loading. This was lower than that observed in conforming contact problems involving biphasic cartilage and was due to a smaller local contact area and variable clearance making fluid exudation easier.     

Hip joint degeneration due to cam impingement: a finite element analysis

The goal of this study was to investigate the impact of cam impingement, a biomechanical risk factor, on hip joint degeneration and ultimately coxarthrosis. 3D finite element solid models of a healthy and a pathologic hip were developed based on clinical reports. The biphasic characteristics of cartilaginous tissues were considered to identify localised solid matrix overloading during normal walking and sitting down (SD). Localised femoral intrusion at the anterior-superior pelvic horn was revealed in the pathologic hip during SD, where the radial and meridional solid stresses in the acetabular cartilage and circumferential solid stresses within the acetabular labrum increased by 3.7, 1.5 and 2.7 times, respectively. The increased solid-on-solid stresses, reduction in fluid-load support and associated higher friction during articulation may result in joint wear and other degenerative changes in the hip.     

Role of the acetabular labrum in load support across the hip joint

The relatively high incidence of labral tears among patients presenting with hip pain suggests that the acetabular labrum is often subjected to injurious loading in vivo. However, it is unclear whether the labrum participates in load transfer across the joint during activities of daily living. This study examined the role of the acetabular labrum in load transfer for hips with normal acetabular geometry and acetabular dysplasia using subject-specific finite element analysis. Models were generated from volumetric CT data and analyzed with and without the labrum during activities of daily living. The labrum in the dysplastic model supported 4-11% of the total load transferred across the joint, while the labrum in the normal model supported only 1-2% of the total load. Despite the increased load transferred to the acetabular cartilage in simulations without the labrum, there were minimal differences in cartilage contact stresses. This was because the load supported by the cartilage correlated with the cartilage contact area. A higher percentage of load was transferred to the labrum in the dysplastic model because the femoral head achieved equilibrium near the lateral edge of the acetabulum. The results of this study suggest that the labrum plays a larger role in load transfer and joint stability in hips with acetabular dysplasia than in hips with normal acetabular geometry.     

Finite element and experimental cortex strains of the intact and implanted tibia

Finite Element (FE) models for the simulation of intact and implanted bone find their main purpose in accurately reproducing the associated mechanical behavior. FE models can be used for preclinical testing of joint replacement implants, where some biomechanical aspects are difficult, if not possible, to simulate and investigate in vitro. To predict mechanical failure or damage, the models should accurately predict stresses and strains. Commercially available synthetic femur models have been extensively used to validate finite element models, but despite the vast literature available on the characteristics of synthetic tibia, numerical and experimental validation of the intact and implant assemblies of tibia are very limited or lacking. In the current study, four FE models of synthetic tibia, intact and reconstructed, were compared against experimental bone strain data, and an overall agreement within 10% between experimental and FE strains was obtained. Finite element and experimental (strain gauge) models of intact and implanted synthetic tibia were validated based on the comparison of cortex bone strains. The study also includes the analysis carried out on standard tibial components with cemented and noncemented stems of the P.F.C Sigma Modular Knee System. The overall agreement within 10% previously established was achieved, indicating that FE models could be successfully validated. The obtained results include a statistical analysis where the root-mean-square-error values were always <10%. FE models can successfully reproduce bone strains under most relevant acting loads upon the condylar surface of the tibia. Moreover, FE models, once properly validated, can be used for preclinical testing of tibial knee replacement, including misalignment of the implants in the proximal tibia after surgery, simulation of long-term failure according to the damage accumulation failure scenario, and other related biomechanical aspects.     

The effect of boundary condition on the biomechanics of a human pelvic joint under an axial compressive load: a three-dimensional finite element model

The finite element (FE) model of the pelvic joint is helpful for clinical diagnosis and treatment of pelvic injuries. However, the effect of an FE model boundary condition on the biomechanical behavior of a pelvic joint has not been well studied. The objective of this study was to study the effect of boundary condition on the pelvic biomechanics predictions. A 3D FE model of a pelvis using subject-specific estimates of intact bone structures, main ligaments and bone material anisotropy by computed tomography (CT) gray value was developed and validated by bone surface strains obtained from rosette strain gauges in an in vitro pelvic experiment. Then three FE pelvic models were constructed to analyze the effect of boundary condition, corresponding to an intact pelvic joint, a pelvic joint without sacroiliac ligaments and a pelvic joint without proximal femurs, respectively. Vertical load was applied to the same pelvis with a fixed prosthetic femoral stem and the same load was simulated in the FE model. A strong correlation coefficient (R(2)=0.9657) was calculated, which indicated a strong correlation between the FE analysis and experimental results. The effect of boundary condition changes on the biomechanical response depended on the anatomical location and structure of the pelvic joint. It was found that acetabulum fixed in all directions with the femur removed can increase the stress distribution on the acetabular inner plate (approximately double the original values) and decrease that on the superior of pubis (from 7 MPa to 0.6 MPa). Taking sacrum and ilium as a whole, instead of sacroiliac and iliolumber ligaments, can influence the stress distribution on ilium and pubis bone vastly. These findings suggest pelvic biomechanics is very dependent on the boundary condition in the FE model.     

The biomechanical effect of the collar of a femoral stem on total hip arthroplasty

To investigate the biomechanical effect of collars, finite element analyses are carried out through two hip joints that are implanted using collared and collarless stems, respectively, and an intact hip joint model. For the analyses, the sacrum, coxal bone, and the cancellous and cortical bones of a femur are modelled using finite elements based on X-ray computed tomographic images taken from a 27-year-old woman. From the results, it is found that a collar with perfect calcar contact prevents stem subsidence and decreases the proximal-lateral gap and the lateral stem tilting. Therefore, it can impart reasonable biomechanical stability for total hip arthroplasty. However, its low load transmission ability and increased stem tilting effect due to the imperfect contact between the collar and the calcar are found to be serious problems that need to be solved. Results of clinical follow-up are presented for supporting the computational results.     

The biomechanical effect of the collar of a femoral stem on total hip arthroplasty

To investigate the biomechanical effect of collars, finite element analyses are carried out through two hip joints that are implanted using collared and collarless stems, respectively, and an intact hip joint model. For the analyses, the sacrum, coxal bone, and the cancellous and cortical bones of a femur are modelled using finite elements based on X-ray computed tomographic images taken from a 27-year-old woman. From the results, it is found that a collar with perfect calcar contact prevents stem subsidence and decreases the proximal–lateral gap and the lateral stem tilting. Therefore, it can impart reasonable biomechanical stability for total hip arthroplasty. However, its low load transmission ability and increased stem tilting effect due to the imperfect contact between the collar and the calcar are found to be serious problems that need to be solved. Results of clinical follow-up are presented for supporting the computational results.     

Normal hip joint contact pressure distribution in single-leg standing--effect of gender and anatomic parameters.

A practical and easy-to-use analysis technique that can study the patient's hip joint contact force/pressure distribution would be useful to assess the effect of abnormal biomechanical conditions and anatomical deformities on joint contact stress for treatment planning purpose. This technique can also help to establish the normative database on hip joint contact pressure distribution in men and women in different age groups. Twelve anatomic parameters and seven biomechanical parameters of the hip joint in a normal population (41 females, 15 males) were calculated. The inter-parameter correlations were investigated. The pressure distribution in the hip joint was calculated using a three-dimensional discrete element analysis (DEA) technique. The 3D contact geometry of the hip joint was estimated from a 2D radiograph by assuming that the femoral head and the acetabular surface were spherical in shape. The head-trochanter ratio (HT), femoral head radius, pelvic height, the joint contact area, the normalized peak contact pressure, abductor force, and the joint contact force were significantly different between men and women. The normalized peak contact pressure was correlated both with acetabular coverage and head-trochanter ratio. Change of abductor force direction within normal variation did not affect the joint peak contact pressure. However, in simulated dysplastic conditions when the CE angle is small or negative, abductor muscle direction becomes very sensitive in joint contact pressure estimation. The models and the results presented can be used as the reference base in computer simulation for preoperative planning in pelvic or femoral osteotomy.     

人工半骨盆假体置换联合腰椎椎弓根螺钉固定术后生物力学的有限元分析

目的:建立人工半骨盆假体置换与联合腰椎椎弓根螺钉固定后的三维有限元模型,评价腰骶段生物力学改变后半骨盆假体力学结构的特点。方法:采用CT薄层扫描采集原始数据,分别建立正常骨盆、半骨盆假体置换术后以及半骨盆假体置换联合腰椎椎弓根螺钉固定术后骨盆的三维有限元模型,分别在第4腰椎上终板平面施以500 N的垂直纵向载荷,分析不同骨盆模型的应力分布特点。结果:与正常骨盆有限元模型相比,半骨盆假体置换术后健侧骨盆应力分布以骶髂关节、髋臼窝及耻骨为主,置换侧半骨盆假体以耻骨连接棒、髋臼杯及髂骨座为主,最大应力出现在耻骨连接棒,应力峰值为65.62 MPa。联合腰椎椎弓根螺钉固定后健侧应力相对减小,置换侧髂骨固定座与骶骨固定处应力相对减小,应力分布以腰椎椎弓根钉棒、耻骨连接棒及髋臼杯为主,最大应力出现在椎弓根螺钉,应力峰值为107 MPa。结论:半骨盆假体置换联合腰椎椎弓根螺钉固定后钉棒分担了半骨盆置换后健侧骨盆及置换侧髂骨固定座与骶骨固定处附近的部分应力,缓解应力集中现象,降低术后骨盆破坏风险,一定程度上增加了半骨盆置换后骨盆的稳定性。     

An in vitro investigation of the acetabular labral seal in hip joint mechanics

Labrum pathology may contribute to early joint degeneration through the alteration of load transfer between, and the stresses within, the cartilage layers of the hip. We hypothesize that the labrum seals the hip joint, creating a hydrostatic fluid pressure in the intra-articular space, and limiting the rate of cartilage layer consolidation. The overall cartilage creep consolidation of six human hip joints was measured during the application of a constant load of 0.75 times bodyweight, or a cyclic sinusoidal load of 0.75+/-0.25 times bodyweight, before and after total labrum resection. The fluid pressure within the acetabular was measured. Following labrum resection, the initial consolidation rate was 22% greater (p=0.02) and the final consolidation displacement was 21% greater (p=0.02). There was no significant difference in the final consolidation rate. Loading type (constant vs. cyclic) had no significant effect on the measured consolidation behaviour. Fluid pressurisation was observed in three of the six hips. The average pressures measured were: for constant loading, 541+/-61kPa in the intact joint and 216+/-165kPa following labrum resection, for cyclic loading, 550+/-56kPa in the intact joint and 195+/-145kPa following labrum resection. The trends observed in this experiment support the predictions of previous finite element analyses. Hydrostatic fluid pressurisation within the intra-articular space is greater with the labrum than without, which may enhance joint lubrication. Cartilage consolidation is quicker without the labrum than with, as the labrum adds an extra resistance to the flow path for interstitial fluid expression. However, both sealing mechanisms are dependent on the fit of the labrum against the femoral head.     

应用三维有限元分析髋臼骨结核软骨下骨塌陷的风险

目的建立人体髋臼骨结核三维有限元模型,探讨不同部位髋臼骨结核软骨下骨塌陷的风险。方法通过正常髋关节CT数据,利用Mimics软件和ANSYS有限元软件,建立正常髋关节三维有限元模型（模型A）、髋臼顶部骨结核（模型B）、髋臼中心部骨结核（模型C）、髋臼前部骨结核（模型D）、髋臼后部骨结核（模型E）三维有限元模型,模拟人体单脚站立进行加载,分析髋臼软骨下骨峰值Von Mises应力和初始微动值。结果建立了正常髋关节和不同部位髋臼骨结核三维有限元模型,各模型含节点269284,三维四面体单元184786个。通过加载分析结果显示：与正常髋关节相比,峰值Von Mises应力,依次增加84%、3%、21%、67%;髋臼软骨下骨初始微动值依次增加66%、11%、17%、29%。结论髋臼顶部结核软骨下骨峰值Von Mises应力和初始微动值最大,塌陷的风险最大。     

The low permeability of healthy meniscus and labrum limit articular cartilage consolidation and maintain fluid load support in the knee and hip

The knee meniscus and hip labrum appear to be important for joint health, but the mechanisms by which these structures perform their functions are not fully understood. The fluid phase of articular cartilage provides compressive stiffness and aids in maintaining a low friction articulation. Healthy fibrocartilage, the tissue of meniscus and labrum, has a lower fluid permeability than articular cartilage. In this study we hypothesized that an important function of the knee meniscus and the hip labrum is to augment fluid retention in the articular cartilage of a mechanically loaded joint. Axisymmetric hyperporoelastic finite element models were analyzed for an idealized knee and an idealized hip. The results indicate that the meniscus maintained fluid pressure and inhibited fluid exudation in knee articular cartilage. Similar, but smaller, effects were seen with the labrum in the hip. Increasing the fibrocartilage permeability relative to that of articular cartilage gave a consolidation rate and loss of fluid load support comparable to that predicted by meniscectomy or labrectomy. The reduced articular cartilage fluid pressure that was calculated for the joint periphery is consistent with patterns of endochondral ossification and osteophyte formation in knee and hip osteoarthritis. High articular central strains and loss of fluid load support after meniscectomy could lead to fibrillation. An intact low-permeability fibrocartilage is important for limiting fluid exudation from articular cartilage in the hip and knee. This may be an important aspect of the role of fibrocartilage in protecting these joints from osteoarthritis.     

The influence and biomechanical role of cartilage split line pattern on tibiofemoral cartilage stress distribution during the stance phase of gait

This study presents an evaluation of the role that cartilage fibre ‘split line’ orientation plays in informing femoral cartilage stress patterns. A two-stage model is presented consisting of a whole knee joint coupled to a tissue-level cartilage model for computational efficiency. The whole joint model may be easily customised to any MRI or CT geometry using free-form deformation. Three ‘split line’ patterns (medial–lateral, anterior–posterior and random) were implemented in a finite element model with constitutive properties referring to this ‘split line’ orientation as a finite element fibre field. The medial–lateral orientation was similar to anatomy and was derived from imaging studies. Model predictions showed that ‘split lines’ are formed along the line of maximum principal strains and may have a biomechanical role of protecting the cartilage by limiting the cartilage deformation to the area of higher cartilage thickness.     

Articular cartilage friction increases in hip joints after the removal of acetabular labrum

The acetabular labrum is believed to have a sealing function. However, a torn labrum may not effectively prevent joint fluid from escaping a compressed joint, resulting in impaired lubrication. We aimed to understand the role of the acetabular labrum in maintaining a low friction environment in the hip joint. We did this by measuring the resistance to rotation (RTR) of the hip, which reflects the friction of the articular cartilage surface, following focal and complete labrectomy. Five cadaveric hips without evidence of osteoarthritis and impingement were tested. We measured resistance to rotation of the hip joint during 0.5, 1, 2, and 3 times body weight (BW) cyclic loading in the intact hip, and after focal and complete labrectomy. Resistance to rotation, which reflects articular cartilage friction in an intact hip was significantly increased following focal labrectomy at 1-3 BW loading, and following complete labrectomy at all load levels. The acetabular labrum appears to maintain a low friction environment, possibly by sealing the joint from fluid exudation. Even focal labrectomy may result in increased joint friction, a condition that may be detrimental to articular cartilage and lead to osteoarthritis.     

髋骨术后仿真生物力学分析与临床研究

目的:探讨坐耻骨肿瘤切除术后髋骨生物力学变化及其与恢复期并发症产生的关系,指导术后骨盆重建.方法:构建单侧完整髋骨及坐耻骨肿瘤切除术后髋骨有限元模型,在相同约束和负载条件下,计算髋骨相同部位位移、应力及应变值,比较与分析模拟手术前后力学变化,结合临床评价意义.结果:手术前后髋骨节点位移、应力及应变分布区域相似,手术前后骶髂关节节点位移及应变、髋臼顶节点应力及应变有显著性差异;手术前后骶髂关节节点应力、坐骨大切迹应力及应变无显著性差异.结论:坐耻骨肿瘤切除术后主要影响同侧骶髂关节、髋臼顶的生物力学状态,恢复期并发症发生的原因由骶髂关节、髋臼顶生物力学改变及耻骨联合与耻骨上支连接与稳定作用的消失所致,三者相互作用,相互影响.     

The mechanics of focal chondral defects in the hip

There is a mean incidence of osteoarthritis (OA) of the hip in 8% of the overall population. In the presence of focal chondral defects, defined as localized damage to the articular cartilage, there is an increased risk of symptomatic progression toward OA. This relationship between chondral defects and subsequent development of OA has led to substantial efforts to develop effective procedures for surgical cartilage repair. This study examined the effects of chondral defects and labral delamination on cartilage mechanics in the dysplastic hip during the gait cycle using subject-specific finite element analysis. Models were generated from volumetric CT data and analyzed with simulated chondral defects at the chondrolabral junction on the posterior acetabulum during five distinct points in the gait cycle. Focal chondral defects increased maximum shear stress on the osteochondral surface of the acetabular cartilage, when compared to the intact case. This effect was amplified with labral delamination. Additionally, chondral defects increased the first principal Lagrange strain on the articular surface of the acetabular cartilage and labrum. Labral delamination relieved some of this tensile strain. As defect size was increased, contact stress increased in the medial zone of the acetabulum, while it decreased anteriorly. The results suggest that in the presence of chondral defects and labral delamination the cartilage experiences elevated tensile strains and shear and contact stress, which could lead to further damage of the cartilage, and subsequent arthritic progression. The framework presented here will serve as the procedure for future finite element studies on cartilage mechanics in hips with varying disease states with simulated chondral defects and labral tears.     

Specimen-specific predictions of contact stress under physiological loading in the human hip: validation and sensitivity studies

Hip osteoarthritis may be initiated and advanced by abnormal cartilage contact mechanics, and finite element (FE) modeling provides an approach with the potential to allow the study of this process. Previous FE models of the human hip have been limited by single specimen validation and the use of quasi-linear or linear elastic constitutive models of articular cartilage. The effects of the latter assumptions on model predictions are unknown, partially because data for the instantaneous behavior of healthy human hip cartilage are unavailable. The aims of this study were to develop and validate a series of specimen-specific FE models, to characterize the regional instantaneous response of healthy human hip cartilage in compression, and to assess the effects of material nonlinearity, inhomogeneity and specimen-specific material coefficients on FE predictions of cartilage contact stress and contact area. Five cadaveric specimens underwent experimental loading, cartilage material characterization and specimen-specific FE modeling. Cartilage in the FE models was represented by average neo-Hookean, average Veronda Westmann and specimen- and region-specific Veronda Westmann hyperelastic constitutive models. Experimental measurements and FE predictions compared well for all three cartilage representations, which was reflected in average RMS errors in contact stress of less than 25 %. The instantaneous material behavior of healthy human hip cartilage varied spatially, with stiffer acetabular cartilage than femoral cartilage and stiffer cartilage in lateral regions than in medial regions. The Veronda Westmann constitutive model with average material coefficients accurately predicted peak contact stress, average contact stress, contact area and contact patterns. The use of subject- and region-specific material coefficients did not increase the accuracy of FE model predictions. The neo-Hookean constitutive model underpredicted peak contact stress in areas of high stress. The results of this study support the use of average cartilage material coefficients in predictions of cartilage contact stress and contact area in the normal hip. The regional characterization of cartilage material behavior provides the necessary inputs for future computational studies, to investigate other mechanical parameters that may be correlated with OA and cartilage damage in the human hip. In the future, the results of this study can be applied to subject-specific models to better understand how abnormal hip contact stress and contact area contribute to OA.