Subject-specific hip geometry affects predicted hip joint contact forces during gait

Hip loading affects bone remodeling and implant fixation. In this study, we have analyzed the effect of subject-specific modeling of hip geometry on muscle activation patterns and hip contact forces during gait, using musculoskeletal modeling, inverse dynamic analysis and static optimization. We first used sensitivity analysis to analyze the effect of isolated changes in femoral neck-length (NL) and neck-shaft angle (NSA) on calculated muscle activations and hip contact force during the stance phase of gait. A deformable generic musculoskeletal model was adjusted incrementally to adopt a physiological range of NL and NSA. In a second similar analysis, we adjusted hip geometry to the measurements from digitized radiographs of 20 subjects with primary hip osteoarthrosis. Finally, we studied the effect of hip abductor weakness on muscle activation patterns and hip contact force. This analysis showed that differences in NL (41-74 mm) and NSA (113-140 degrees ) affect the muscle activation of the hip abductors during stance phase and hence hip contact force by up to three times body weight. In conclusion, the results from both the sensitivity and subject-specific analysis showed that at the moment of peak contact force, altered NSA has only a minor effect on the loading configuration of the hip. Increased NL, however, results in an increase of the three hip contact-force components and a reduced vertical loading. The results of these analyses are essential to understand modified hip joint loading, and for planning hip surgery for patients with osteoarthrosis.     

Subject-specific hip geometry and hip joint centre location affects calculated contact forces at the hip during gait

Hip loading affects the development of hip osteoarthritis, bone remodelling and osseointegration of implants. In this study, we analyzed the effect of subject-specific modelling of hip geometry and hip joint centre (HJC) location on the quantification of hip joint moments, muscle moments and hip contact forces during gait, using musculoskeletal modelling, inverse dynamic analysis and static optimization. For 10 subjects, hip joint moments, muscle moments and hip loading in terms of magnitude and orientation were quantified using three different model types, each including a different amount of subject-specific detail: (1) a generic scaled musculoskeletal model, (2) a generic scaled musculoskeletal model with subject-specific hip geometry (femoral anteversion, neck-length and neck-shaft angle) and (3) a generic scaled musculoskeletal model with subject-specific hip geometry including HJC location. Subject-specific geometry and HJC location were derived from CT. Significant differences were found between the three model types in HJC location, hip flexion–extension moment and inclination angle of the total contact force in the frontal plane. No model agreement was found between the three model types for the calculation of contact forces in terms of magnitude and orientations, and muscle moments. Therefore, we suggest that personalized models with individualized hip joint geometry and HJC location should be used for the quantification of hip loading. For biomechanical analyses aiming to understand modified hip joint loading, and planning hip surgery in patients with osteoarthritis, the amount of subject-specific detail, related to bone geometry and joint centre location in the musculoskeletal models used, needs to be considered.     

Subject-specific geometrical detail rather than cost function formulation affects hip loading calculation*

This study assessed the relative importance of introducing an increasing level of medical image-based subject-specific detail in bone and muscle geometry in the musculoskeletal model, on calculated hip contact forces during gait. These forces were compared to introducing minimization of hip contact forces in the optimization criterion. With an increasing level of subject-specific detail, specifically MRI-based geometry and wrapping surfaces representing the hip capsule, hip contact forces decreased and were more comparable to contact forces measured using instrumented prostheses (average difference of 0.69 BW at the first peak compared to 1.04 BW for the generic model). Inclusion of subject-specific wrapping surfaces in the model had a greater effect than altering the cost function definition.     

Achilles tendon stress is more sensitive to subject-specific geometry than subject-specific material properties: A finite element analysis

This study used subject-specific measures of three-dimensional (3D) free Achilles tendon geometry in conjunction with a finite element method to investigate the effect of variation in subject-specific geometry and subject-specific material properties on tendon stress during submaximal isometric loading. Achilles tendons of eight participants (Aged 25–35 years) were scanned with freehand 3D ultrasound at rest and during a 70% maximum voluntary isometric contraction. Ultrasound images were segmented, volume rendered and transformed into subject-specific 3D finite element meshes. The mean (±SD) lengths, volumes and cross-sectional areas of the tendons at rest were 62 ± 13 mm, 3617 ± 984 mm³ and 58 ± 11 mm² respectively. The measured tendon strain at 70% MVIC was 5.9 ± 1.3%. Subject-specific material properties were obtained using an optimisation approach that minimised the difference between measured and modelled longitudinal free tendon strain. Generic geometry was represented by the average mesh and generic material properties were taken from the literature. Local stresses were subsequently computed for combinations of subject-specific and generic geometry and material properties. For a given geometry, changing from generic to subject-specific material properties had little effect on the stress distribution in the tendon. In contrast, changing from generic to subject-specific geometry had a 26-fold greater effect on tendon stress distribution. Overall, these findings indicate that the stress distribution experienced by the living free Achilles tendon of a young and healthy population during voluntary loading are more sensitive to variation in tendon geometry than variation in tendon material properties.     

A superellipsoid-plane model for simulating foot-ground contact during human gait

Musculoskeletal models and forward dynamics simulations of human movement often include foot–ground interactions, with the foot–ground contact forces often determined using a constitutive model that depends on material properties and contact kinematics. When using soft constraints to model the foot–ground interactions, the kinematics of the minimum distance between the foot and planar ground needs to be computed. Due to their geometric simplicity, a considerable number of studies have used point–plane elements to represent these interacting bodies, but few studies have provided comparisons between point contact elements and other geometrically based analytical solutions. The objective of this work was to develop a more general-purpose superellipsoid–plane contact model that can be used to determine the three-dimensional foot–ground contact forces. As an example application, the model was used in a forward dynamics simulation of human walking. Simulation results and execution times were compared with a point-like viscoelastic contact model. Both models produced realistic ground reaction forces and kinematics with similar computational efficiency. However, solving the equations of motion with the surface contact model was found to be more efficient (~18% faster), and on average numerically ~37% less stiff. The superellipsoid–plane elements are also more versatile than point-like elements in that they allow for volumetric contact during three-dimensional motions (e.g. rotating, rolling, and sliding). In addition, the superellipsoid–plane element is geometrically accurate and easily integrated within multibody simulation code. These advantages make the use of superellipsoid–plane contact models in musculoskeletal simulations an appealing alternative to point-like elements.     

Non-rigid deformation to include subject-specific detail in musculoskeletal models of CP children with proximal femoral deformity and its effect on muscle and contact forces during gait

To account for proximal femoral deformities in children with cerebral palsy (CP), subject-specific musculoskeletal models are needed. Non-rigid deformation (NRD) deforms generic onto personalized bone geometry and thereby transforms the muscle points. The goal of this study was to determine to what extent the models and simulation outcomes in CP patients differ when including subject-specific detail using NRD or Magnetic Resonance Imaging (MRI)-based models. The NRD models slightly overestimated hip contact forces compared to MRI models and differences in muscle point positions and moment arm lengths (MALs) remained, although differences were smaller than for the generic model.     

A novel computational method to determine subject-specific bite force and occlusal loading during mastication

The evaluation of three-dimensional occlusal loading during biting and chewing may assist in development of new dental materials, in designing effective and long-lasting restorations such as crowns and bridges, and for evaluating functional performance of prosthodontic components such as dental and/or maxillofacial implants. At present, little is known about the dynamic force and pressure distributions at the occlusal surface during mastication, as these quantities cannot be measured directly. The aim of this study was to evaluate subject-specific occlusal loading forces during mastication using accurate jaw motion measurements. Motion data was obtained from experiments in which an individual performed maximal effort dynamic chewing cycles on a rubber sample with known mechanical properties. A finite element model simulation of one recorded chewing cycle was then performed to evaluate the deformation of the rubber. This was achieved by imposing the measured jaw motions on a three-dimensional geometric surface model of the subject’s dental impressions. Based on the rubber’s deformation and its material behaviour, the simulation was used to compute the resulting stresses within the rubber as well as the contact pressures and forces on the occlusal surfaces. An advantage of this novel modelling approach is that dynamic occlusal pressure maps and biting forces may be predicted with high accuracy and resolution at each time step throughout the chewing cycle. Depending on the motion capture technique and the speed of simulation, the methodology may be automated in such a way that it can be performed chair-side. The present study demonstrates a novel modelling methodology for evaluating dynamic occlusal loading during biting or chewing.     

Finite Element Analysis of the Contact Mechanics of Ceramic-on-Ceramic Hip Resurfacing Prostheses

<正> Ceramics are good alternative to metal as bearing couple materials because of their better wear resistance. A Finite Element(FE) study was performed to investigate the contact mechanics and stress distribution of Ceramic-on-Ceramic (COC) hip resurfacingprostheses. It was focused in particular on a parametric study to examine the effects of radial clearance, loading,alumina coating on the implants, bone quality, and fixation of cup-bone interface. It was found that a reduction in the radialclearance had the most significant effect on the predicted contact pressure distribution among all of the parameters considered inthis study. It was determined that there was a significant influence of non-metallic materials, such as the bone underneath thebearing components, on the predicted contact mechanics. Stress shielding within the bone tissue was found to be a major concernwhen regarding the use of ceramic as an alternative to metallic resurfacing prostheses. Therefore, using alumina implantswith a metal backing was found to be the best design for ceramic resurfacing prostheses in this study. The loading, bone quality,and acetabular cup fixation conditions were found to have only minor effects on the predicted contact pressure distribution alongthe bearing surfaces.     

Analytical and numerical analysis of human dental occlusal contact

The knowledge of contact forces in teeth surfaces during mastication or para-functional movements can help to understand processes related to friction and wear of human dental enamel. The development of a numerical model for analysis of the occlusal contact between two antagonistic teeth is proposed, which includes three basic steps: the characterisation of the surface roughness, its homogenisation using an assumed distribution function and the numerical determination of the resulting forces. Finite element strain results for the main different asperities are statistically combined, deriving the predicted macroscopic behaviour of the interface. Axisymmetric and 3D numerical models with an elasto-plastic constitutive law are used to simulate micro-indentations and micro-contacts, respectively. The contact is allowed to occur locally in planes not necessarily parallel to the surface's mean plane, a problem for which there is no analytical solution. The three identified parameters, homogenised surface hardness (3.68 GPa), surface yield stress (3.08 GPa) and static friction coefficient (0.23), agree with the experimental values reported in the literature.     

A multibody knee model with discrete cartilage prediction of tibio-femoral contact mechanics

Combining musculoskeletal simulations with anatomical joint models capable of predicting cartilage contact mechanics would provide a valuable tool for studying the relationships between muscle force and cartilage loading. As a step towards producing multibody musculoskeletal models that include representation of cartilage tissue mechanics, this research developed a subject-specific multibody knee model that represented the tibia plateau cartilage as discrete rigid bodies that interacted with the femur through deformable contacts. Parameters for the compliant contact law were derived using three methods: (1) simplified Hertzian contact theory, (2) simplified elastic foundation contact theory and (3) parameter optimisation from a finite element (FE) solution. The contact parameters and contact friction were evaluated during a simulated walk in a virtual dynamic knee simulator, and the resulting kinematics were compared with measured in vitro kinematics. The effects on predicted contact pressures and cartilage–bone interface shear forces during the simulated walk were also evaluated. The compliant contact stiffness parameters had a statistically significant effect on predicted contact pressures as well as all tibio-femoral motions except flexion–extension. The contact friction was not statistically significant to contact pressures, but was statistically significant to medial–lateral translation and all rotations except flexion–extension. The magnitude of kinematic differences between model formulations was relatively small, but contact pressure predictions were sensitive to model formulation. The developed multibody knee model was computationally efficient and had a computation time 283 times faster than a FE simulation using the same geometries and boundary conditions.     

Extracellular matrix integrity affects the mechanical behaviour of in-situ chondrocytes under compression

Cartilage lesions change the microenvironment of cells and may accelerate cartilage degradation through catabolic responses from chondrocytes. In this study, we investigated the effects of structural integrity of the extracellular matrix (ECM) on chondrocytes by comparing the mechanics of cells surrounded by an intact ECM with cells close to a cartilage lesion using experimental and numerical methods. Experimentally, 15% nominal compression was applied to bovine cartilage tissues using a light-transmissible compression system. Target cells in the intact ECM and near lesions were imaged by dual-photon microscopy. Changes in cell morphology (N_cell=32 for both ECM conditions) were quantified. A two-scale (tissue level and cell level) Finite Element (FE) model was also developed. A 15% nominal compression was applied to a non-linear, biphasic tissue model with the corresponding cell level models studied at different radial locations from the centre of the sample in the transient phase and at steady state. We studied the Green-Lagrange strains in the tissue and cells. Experimental and theoretical results indicated that cells near lesions deform less axially than chondrocytes in the intact ECM at steady state. However, cells near lesions experienced large tensile strains in the principal height direction, which are likely associated with non-uniform tissue radial bulging. Previous experiments showed that tensile strains of high magnitude cause an up-regulation of digestive enzyme gene expressions. Therefore, we propose that cartilage degradation near tissue lesions may be due to the large tensile strains in the principal height direction applied to cells, thus leading to an up-regulation of catabolic factors.     

Finite element modeling for predicting the contact pressure between a foam mattress and the human body in a supine position

In this paper, we generated finite element (FE) models to predict the contact pressure between a foam mattress and the human body in a supine position. Twenty-year-old males were used for three-dimensional scanning to produce the FE human models, which was composed of skin and muscle tissue. A linear elastic isotropic material model was used for the skin, and the Mooney–Rivlin model was used for the muscle tissue because it can effectively represent the nonlinear behavior of muscle. The contact pressure between the human model and the mattress was predicted by numerical simulation. The human models were validated by comparing the body pressure distribution obtained from the same human subject when he was lying on two different mattress types. The experimental results showed that the slope of the lower part of the mattress caused a decrease in the contact pressure at the heels, and the effect of bone structure was most pronounced in the scapula. After inserting a simple structure to function as the scapula, the contact pressure predicted by the FE human models was consistent with the experimental body pressure distribution for all body parts. These results suggest that the models proposed in this paper will be useful to researchers and designers of products related to the prevention of pressure ulcers.     

Dynamic contact stress patterns on the tibial plateaus during simulated gait: A novel application of normalized cross correlation

The spatial distribution and pattern of local contact stresses within the knee joint during activities of daily living have not been fully investigated. The objective of this study was to determine if common contact stress patterns exist on the tibial plateaus of human knees during simulated gait. To test this hypothesis, we developed a novel normalized cross-correlation (NCC) algorithm and applied it to the contact stresses on the tibial plateaus of 12 human cadaveric knees subjected to multi-directional loads mimicking gait. The contact stress profiles at different locations on the tibial plateaus were compared, where regions with similar contact stress patterns were identified across specimens. Three consistent regional patterns were found, among them two most prominent contact stress patterns were shared by 9–12 of all the knees and the third pattern was shared by 6–8 knees. The first pattern was located at the posterior aspect of the medial tibial plateau and had a single peak stress that occurred during the early stance phase. The second pattern was located at the central-posterior aspects of the lateral plateau and consisted of two peak stresses coincident with the timing of peak axial force at early and late stance. The third pattern was found on the anterior aspect of cartilage-to-cartilage contact region on the medial plateau consisted of double peak stresses. The differences in the location and profile of the contact stress patterns suggest that the medial and lateral menisci function to carry load at different points in the gait cycle: with the posterior aspect of the medial meniscus consistently distributing load only during the early phase of stance, and the posterior aspect of the lateral meniscus consistently distributing load during both the early and late phases of stance. This novel approach can help identify abnormalities in knee contact mechanics and provide a better understanding of the mechanical pathways leading to post-traumatic osteoarthritis.     

Effect of prism orientation and loading direction on contact stresses in prismatic enamel of primates: implications for interpreting wear patterns

The ability of prisms to effectively dissipate contact stress at the surface will influence wear rates in teeth. The aim of this investigation was to begin to quantify the effect of prism orientation on surface stresses. Seven finite element models of enamel microstructure were created, each model differing in the angulation of prism orientation with regard to the wear surface. For validation purposes, the mechanical behavior of the model was compared with published experimental data. In order to test the enamel under lateral loads, a compressed food particle was dragged across the surface from the dentino-enamel junction (DEJ) towards the outer enamel surface (OES). Under these conditions, tensile stresses in the enamel model increased with increases in the coefficient of friction. More importantly, stresses were found to be lowest in models in which the prisms approach the surface at lower angles (i.e., more obliquely cut prisms), and highest when the prisms approached the surface at 60 degrees (i.e., less obliquely cut). Finally, the direction of travel of the simulated food particle was reversed, allowing comparison of the difference in behavior between trailing and leading edge enamels (i.e., when the food particle was dragged either towards or away from the DEJ). Stresses at the trailing edge were usually lower than stresses at the leading edge. Taken together with what is known about prism orientation in primate teeth, such findings imply greater wear resistance at the intercuspal region and less wear resistance at the lateral enamel at midcrown. Such findings appear to be supported by archeological evidence.     

Numerical simulation of the insertion process of an uncemented hip prosthesis in order to evaluate the influence of residual stress and contact distribution on the stem initial stability

The long-term success of a cementless total hip arthroplasty depends on the implant geometry and interface bonding characteristics (fit, coating and ingrowth) and on stem stiffness. This study evaluates the influence of stem geometry and fitting conditions on the evolution and distribution of the bone–stem contact, stress and strain during and after the hip stem insertion, by means of dynamic finite element techniques. Next, the influence of the mechanical state (bone–stem contact, stress and strain) resulted from the insertion process on the stem initial resistance to subsidence is investigated. In addition, a study on the influence of bone–stem interface conditions (friction) on the insertion process and on the initial stem stability under physiological loading is performed. The results indicate that for a stem with tapered shape the contact in the proximal part of the stem was improved, but contact in the calcar region was achieved only when extra press-fit conditions were considered. Changes in stem geometry towards a more tapered shape and extra press fit and variation in the bone–stem interface conditions (contact amount and high friction) led to a raise in the total insertion force. A direct positive relationship was found between the stem resistance to subsidence and stem geometry (tapering and press fit), bone–stem interface conditions (bone–stem contact and friction interface) and the mechanical status at the end of the insertion (residual stress and strain). Therefore, further studies on evaluating the initial performance of different stem types should consider the parameters describing the bone–stem interface conditions and the mechanical state resulted from the insertion process.     

A 3D lower limb musculoskeletal model for simultaneous estimation of musculo-tendon,joint contact,ligament and bone forces during gait

Musculo-tendon forces and joint reaction forces are typically estimated using a two-step method, computing first the musculo-tendon forces by a static optimization procedure and then deducing the joint reaction forces from the force equilibrium. However, this method does not allow studying the interactions between musculo-tendon forces and joint reaction forces in establishing this equilibrium and the joint reaction forces are usually overestimated. This study introduces a new 3D lower limb musculoskeletal model based on a one-step static optimization procedure allowing simultaneous musculo-tendon, joint contact, ligament and bone forces estimation during gait. It is postulated that this approach, by giving access to the forces transmitted by these musculoskeletal structures at hip, tibiofemoral, patellofemoral and ankle joints, modeled using anatomically consistent kinematic models, should ease the validation of the model using joint contact forces measured with instrumented prostheses. A blinded validation based on four datasets was made under two different minimization conditions (i.e., C1 – only musculo-tendon forces are minimized, and C2 – musculo-tendon, joint contact, ligament and bone forces are minimized while focusing more specifically on tibiofemoral joint contacts). The results show that the model is able to estimate in most cases the correct timing of musculo-tendon forces during normal gait (i.e., the mean coefficient of active/inactive state concordance between estimated musculo-tendon force and measured EMG envelopes was C1: 65.87% and C2: 60.46%). The results also showed that the model is potentially able to well estimate joint contact, ligament and bone forces and more specifically medial (i.e., the mean RMSE between estimated joint contact force and in vivo measurement was C1: 1.14BW and C2: 0.39BW) and lateral (i.e., C1: 0.65BW and C2: 0.28BW) tibiofemoral contact forces during normal gait. However, the results remain highly influenced by the optimization weights that can bring to somewhat aphysiological musculo-tendon forces.