Search

Implant malalignment has been reported to be a primary reason for revision total knee arthroplasty (TKA). In addition, altered muscle coordination patterns are commonly observed in TKA patients, which is thought to alter knee contact loads. A comprehensive understanding of the influence of surgical implantation and muscle recruitment strategies on joint contact mechanics is crucial to improve surgical techniques, increase implant longevity, and inform rehabilitation protocols. In this study, a detailed musculoskeletal model with a 12 degrees of freedom knee was developed to represent a TKA subject from the CAMS-Knee datasets. Using motion capture and ground reaction force data, a level walking cycle was simulated and the joint movement and loading patterns were estimated using a novel technique for concurrent optimization of muscle activations and joint kinematics. In addition, over 12′000 Monte Carlo simulations were performed to predict knee contact mechanics during walking, considering numerous combinations of implant alignment and muscle activation scenarios. Validation of our baseline simulation showed good agreement between the model kinematics and loading patterns against the in vivo data. Our analyses reveal a considerable impact of implant alignment on the joint kinematics, while variation in muscle activation strategies mainly affects knee contact loading. Moreover, our results indicate that high knee compressive forces do not necessarily originate from extreme kinematics and vice versa. This study provides an improved understanding of the complex inter-relationships between loading and movement patterns resulting from different surgical implantation and muscle coordination strategies and presents a validated framework towards population-based modelling in TKA. [ABSTRACT FROM AUTHOR]

Representation of realistic muscle geometries is needed for systematic biomechanical simulation of musculoskeletal systems. Most of the previous musculoskeletal models are based on multibody dynamics simulation with muscles simplified as one-dimensional (1D) line-segments without accounting for the large muscle attachment areas, spatial fibre alignment within muscles and contact and wrapping between muscles and surrounding tissues. In previous musculoskeletal models with three-dimensional (3D) muscles, contractions of muscles were among the inputs rather than calculated, which hampers the predictive capability of these models. To address these issues, a finite element musculoskeletal model with the ability to predict contractions of 3D muscles was developed. Muscles with realistic 3D geometry, spatial muscle fibre alignment and muscle-muscle and muscle-bone interactions were accounted for. Active contractile stresses of the 3D muscles were determined through an efficient optimization approach based on the measured kinematics of the lower extremity and ground force during gait. This model also provided stresses and strains of muscles and contact mechanics of the muscle-muscle and muscle-bone interactions. The total contact force of the knee predicted by the model corresponded well to the in vivo measurement. Contact and wrapping between muscles and surrounding tissues were evident, demonstrating the need to consider 3D contact models of muscles. This modelling framework serves as the methodological basis for developing musculoskeletal modelling systems in finite element method incorporating 3D deformable contact models of muscles, joints, ligaments and bones. [ABSTRACT FROM AUTHOR]

Outcomes of total knee arthroplasty (TKA) are dependent on surgical technique, patient variability, and implant design. Non-optimal design or alignment choices may result in undesirable contact mechanics and joint kinematics, including poor joint alignment, instability, and reduced range of motion. Implant design and surgical alignment are modifiable factors with potential to improve patient outcomes, and there is a need for robust implant designs that can accommodate patient variability. Our objective was to develop a statistical shape-function model (SFM) of a posterior stabilized implanted knee to instantaneously predict joint mechanics in an efficient manner. Finite element methods were combined with Latin hypercube sampling and regression analyses to produce modeling equations relating nine implant design and six surgical alignment parameters to tibiofemoral (TF) joint mechanics outcomes during a deep knee bend. A SFM was developed and TF contact mechanics, kinematics, and soft tissue loads were instantaneously predicted from the model. Average normalized root-mean-square error predictions were between 2.79% and 9.42%, depending on the number of parameters included in the model. The statistical shape-function model generated instantaneous joint mechanics predictions using a maximum of 130 training simulations, making it ideally suited for integration into a patient-specific design and alignment optimization pipeline. Such a tool may be used to optimize kinematic function to achieve more natural motion or minimize implant wear, and may aid the engineering and clinical communities in improving patient satisfaction and surgical outcomes. [ABSTRACT FROM AUTHOR]

Abstract Modeling joint contact in OpenSim is not well understood. This study systematically investigated the variables associated with the elastic foundation contact model within OpenSim by performing a series of controlled benchtop experiment and concomitant simulations. Four metal-on-plastic interactions were modeled, including a model of a total knee replacement (TKR). Load-displacement curves were recorded during cyclic loading between 100 and 750 N. Geometries were imported and into OpenSim and contact mechanics were modeled with the on-board elastic foundation algorithm. A hybrid optimization algorithm determined that stiffness and dissipation coefficients for TKR implants were 1.52 × 1010 N/m and 57.7 Ns/m, respectively. Estimations of contact forces were 10.2% of blinded experimental data (average root mean square error: 76.82 ± 11.47 N). In the second portion of this study, freely available eTibia TKR renderings were used to test the ubiquity of the tuning parameters. They were also used to perform a sensitivity analysis of material stiffness and mesh density with regard to penetration depth and computational time. When a stiffness of 1 × 1010 was applied to an eTibia model with 5000 faces, a 100 kg load caused 0.259 mm of penetration. Under the same conditions, the tuned model experienced 0.300 mm of penetration. Material stiffnesses between 1 × 1013 and 1 × 1015 N/m increased computation time by factors of 12–23. This study provides much needed clarity regarding the use of the OpenSim EF algorithm. It also demonstrates the utility of OpenSim to model deformable materials and complex geometries, and this approach can be adapted to make reasonable estimations for both natural and surgically modified joints. [ABSTRACT FROM AUTHOR]

It has been suggested that the repetitive nature of altered joint tissue loading which occurs after anterior cruciate ligament (ACL) rupture can contribute to the development of osteoarthritis (OA). However, changes in dynamic knee joint contact stresses after ACL rupture have not been quantified for activities of daily living. Our objective was to characterize changes in dynamic contact stress profiles that occur across the tibial plateau immediately after ACL transection. By subjecting sensor-augmented cadaveric knees to simulated gait, and analyzing the resulting contact stress profiles using a normalized cross-correlation algorithm, we tested the hypothesis that common changes in dynamic contact stress profiles exist after ACL injury. Three common profiles were identified in intact knees, occurring on the: (I) posterior lateral plateau, (II) posterior medial plateau, and (III) central region of the medial plateau. In ACL-transected knees, the magnitude and shape of the common dynamic stress profiles did not change, but their locations on the tibial plateau and the number of knees identified for each profile changed. Furthermore, in the ACL transected knees, a unique common contact stress profile was identified in the posterior region of the lateral plateau near the tibial spine. This framework can be used to understand the regional and temporal changes in joint mechanics after injury. [ABSTRACT FROM AUTHOR]

The critical role of cell mechanics in tissue health has led to the development of many in vitro methods that measure the elasticity of the cytoskeleton and whole cells, yet the connection between these local cell properties and bulk measurements of tissue mechanics remains unclear. To help bridge this gap, we have developed a monolayer indentation technique for measuring multi-cellular mechanics in vitro . Here, we measure the elasticity of cell monolayers and uncover the role of fluid permeability in these multi-cellular systems, finding that the resistance of fluid transport through cells controls their force–response at long times. [ABSTRACT FROM AUTHOR]

The posterior medial meniscal root (PMMR) experiences variable and multiaxial forces during loading. Current methods to measure these forces are limited and fail to adequately characterize the loads in all three dimensions at the root. Our novel technique resolved these limitations with the installation of a 3-axis sensing construct that we hypothesized would not affect contact mechanics, would not impart extraneous loads onto the PMMR, would accurately measure forces, and would not deflect under joint loads. Six cadaveric specimens were dissected to the joint capsule and a sagittal-plane, femoral condyle osteotomy was performed to gain access to the root. The load sensor was placed below the PMMR and was validated across four tests. The contact mechanics test demonstrated a contact area precision of 44 mm2 and a contact pressure precision of 5.0 MPa between the pre-installation and post-installation states. The tibial displacement test indicated an average bone plug displacement of < 1 mm in all directions. The load validation test exhibited average precision values of 0.7 N in compression, 0.5 N in tension, 0.3 N in anterior-posterior shear, and 0.3 N in medial–lateral shear load. The bone plug deflection test confirmed < 2 mm of displacement in any direction when placed under a load. This is the first study to successfully validate a technique for measuring both magnitude and direction of forces experienced at the PMMR. This validated method has applications for improving surgical repair techniques and developing safer rehabilitation and postoperative protocols that decrease root loads. [ABSTRACT FROM AUTHOR]

Femoroacetabular impingement syndrome (FAIS) is a motion-related pathology of the hip characterized by pain, morphological abnormalities of the proximal femur, and an elevated risk of joint deterioration and hip osteoarthritis. Activities that require deep flexion are understood to induce impingement in cam FAIS patients, however, less demanding activities such as walking and pivoting may induce pain as well as alterations in kinematics and joint stability. Still, the paucity of quantitative descriptions of cam FAIS has hindered understanding underlying hip joint mechanics during such activities. Previous in silico studies have employed generalized model geometry or kinematics to simulate impingement between the femur and acetabulum, which may not accurately capture the interplay between morphology and motion. In this study, we utilized models with participant-specific bone and articular soft tissue anatomy and kinematics measured by dual-fluoroscopy to compare hip contact mechanics of cam FAIS patients to controls during four activities of daily living (internal/external pivoting and level/incline walking). Averaged across the gait cycle during incline walking, patients displayed increased strain in the anterior joint (labrum strain: p-value = 0.038, patients: 11.7 ± 6.7 %, controls: 5.0 ± 3.6 %; cartilage strain: p-value = 0.029, patients: 9.1 ± 3.3 %, controls: 4.2 ± 2.3). Patients also exhibited increased average anterior cartilage strains during external pivoting (p-value = 0.039; patients: 13.0 ± 9.2 %, controls: 3.9 ± 3.2 %]). No significant differences between patient and control contact area and strain were found for level walking and internal pivoting. Our study provides new insights into the biomechanics of cam FAIS, including spatiotemporal hip joint contact mechanics during activities of daily living. [ABSTRACT FROM AUTHOR]

The aim of this study was to develop a finite element (FE) hip model with subject-specific geometry and biphasic cartilage properties. Different levels of detail in the representation of fibre reinforcement were considered to evaluate the feasibility to simplify the complex depth-dependent fibre pattern in the native hip joint. A FE model of a cadaveric hip with subject-specific geometry was constructed through micro-computed-tomography (µCT) imaging. The cartilage was assumed to be biphasic and fibre-reinforced with different levels of detail in the fibre representation. Simulations were performed for heel-strike, mid-stance and toe-off during walking and one-leg-stance over 1500 s. It was found that the required level of detail in fibre representation depends on the parameter of interest. The contact stress of the native hip joint could be realistically predicted by simplifying the fibre representation to being orthogonally reinforced across the whole thickness. To predict the fluid pressure, depth-dependent fibre organisation is needed but specific split-line pattern on the surface of cartilage is not necessary. Both depth-dependent and specific surface fibre orientations are required to simulate the strains. [ABSTRACT FROM AUTHOR]

Although it is easier to generate finite element discretizations with tetrahedral elements, trilinear hexahedral (HEX8) elements are more often used in simulations of articular contact mechanics. This is due to numerical shortcomings of linear tetrahedral (TET4) elements, limited availability of quadratic tetrahedron elements in combination with effective contact algorithms, and the perceived increased computational expense of quadratic finite elements. In this study we implemented both ten-node (TET10) and fifteen-node (TET15) quadratic tetrahedral elements in FEBio ( www.febio.org ) and compared their accuracy, robustness in terms of convergence behavior and computational cost for simulations relevant to articular contact mechanics. Suitable volume integration and surface integration rules were determined by comparing the results of several benchmark contact problems. The results demonstrated that the surface integration rule used to evaluate the contact integrals for quadratic elements affected both convergence behavior and accuracy of predicted stresses. The computational expense and robustness of both quadratic tetrahedral formulations compared favorably to the HEX8 models. Of note, the TET15 element demonstrated superior convergence behavior and lower computational cost than both the TET10 and HEX8 elements for meshes with similar numbers of degrees of freedom in the contact problems that we examined. Finally, the excellent accuracy and relative efficiency of these quadratic tetrahedral elements was illustrated by comparing their predictions with those for a HEX8 mesh for simulation of articular contact in a fully validated model of the hip. These results demonstrate that TET10 and TET15 elements provide viable alternatives to HEX8 elements for simulation of articular contact mechanics. [ABSTRACT FROM AUTHOR]

Our objective was to quantify the effect of ACL transection on dynamic knee joint contact force distributions during simulated gait. Given the prevalence of medial compartment osteoarthritis in un-reconstructed ACL ruptured knees, we hypothesized that changes in contact mechanics after ACL transection would be most prevalent in the medial compartment. Twelve human cadaveric knees were tested using a dynamic knee gait simulator which was programmed to mimic a clinical Lachman exam and gait. An electronic pressure sensor was placed on the medial and lateral tibial plateaus under the menisci to quantify dynamic contact forces before and after ACL transection. Tibial translations and rotations, medial and lateral plateau peak contact stress, and position and velocity of the Weighted Center of Contact (WCoC) were computed. After ACL transection, the tibia translated more anteriorly in the Lachman examination and at heel strike during gait. Changes in contact mechanics across the medial tibial plateau during simulated gait were: an increase in the velocity of WCoC and a posterior shift in the WCoC, both of which occurred at heel strike; increased peak contact forces in the posterior-peripheral quadrant of the tibial plateau at 45% of the gait cycle; and an additional posterior shift in WCoC from 25 to 55% of the gait cycle. The only change in contact mechanics in the lateral plateau was a decrease in WCoC velocity in late stance. This data is suggested to further the study of biomechanical pathways (biomechanical biomarkers) in the relationship between altered knee contact mechanics and chondrocyte metabolic responses after ACL transection. [ABSTRACT FROM AUTHOR]

Erosion of articular cartilage is a concern following distal humeral hemiarthroplasty, because native cartilage surfaces are placed in contact with stiff metallic implant components, which causes decreases in contact area and increases in contact stresses. Recently, reverse engineered implants have been proposed which are intended to promote more natural contact mechanics by reproducing the native bone or cartilage shape. In this study, finite element modeling is used in order to calculate changes in cartilage contact areas and stresses following distal humeral hemiarthroplasty with commercially available and reverse engineered implant designs. At the ulna, decreases in contact area were −34±3% ( p =0.002), −27±1% ( p <0.001) and −14±2% ( p =0.008) using commercially available, bone reverse engineered and cartilage reverse engineered designs, respectively. Peak contact stresses increased by 461±57% ( p =0.008), 387±127% ( p =0.229) and 165±16% ( p =0.003). At the radius, decreases in contact area were −21±3% ( p =0.013), −13±2% ( p <0.006) and −6±1% ( p =0.020), and peak contact stresses increased by 75±52% ( p >0.999), 241±32% ( p =0.010) and 61±10% ( p =0.021). Between the three different implant designs, the cartilage reverse engineered design yielded the largest contact areas and lowest contact stresses, but was still unable to reproduce the contact mechanics of the native joint. These findings align with a growing body of evidence indicating that although reverse engineered hemiarthroplasty implants can provide small improvements in contact mechanics when compared with commercially available designs, further optimization of shape and material properties is required in order reproduce native joint contact mechanics. [ABSTRACT FROM AUTHOR]

The biomechanics of distance running are studied in relation to both understanding injury mechanisms and improving performance. Kinematic methods must be used to identify the stance phase of running when data are recorded during running on a standard treadmill or outside the laboratory. Recently, a focus on foot strike patterns has emerged in the field. Thus, there is a need for a kinematic method to identify foot contact that is equally effective for both rearfoot and non-rearfoot strike patterns. The purpose of this study was to determine whether a new kinematic method could accurately determine foot contact during running in both rearfoot and non-rearfoot strikers. Overground gait data were collected at on 22 runners, 11 with a rearfoot strike pattern and 11 with a non-rearfoot strike pattern. Data were processed to identify foot contact from: vertical ground reaction force, two previously published kinematic methods, and our new kinematic method. Limits of agreement were used to determine bias and random error of each kinematic method compared to ground reaction force onset. The new method had comparable random error at 200 Hz sampling frequency (5 ms per frame) to the previous methods (7 frames vs 6–9 frames) and produced the same offset for both strike patterns (3 frames), while the existing methods had different offsets for different strike patterns (4 or 7 frames). Study findings support use of this new method, as it can be applied to all running strike patterns without adjusting the frame offset, simplifying data processing. [ABSTRACT FROM AUTHOR]

Although the contact pressure increases during implantation of a wedge-shaped implant, friction coefficients tend to be measured under constant contact pressure, as endorsed in standard procedures. Abrasion and plastic deformation of the bone during implantation are rarely reported, although they define the effective interference, by reducing the nominal interference between implant and bone cavity. In this study radial forces were analysed during simulated implantation and explantation of angled porous and polished implant surfaces against trabecular bone specimens, to determine the corresponding friction coefficients. Permanent deformation was also analysed to determine the effective interference after implantation. For the most porous surface tested, the friction coefficient initially increased with increasing normal contact stress during implantation and then decreased at higher contact stresses. For a less porous surface, the friction coefficient increased continually with normal contact stress during implantation but did not reach the peak magnitude measured for the rougher surface. Friction coefficients for the polished surface were independent of normal contact stress and much lower than for the porous surfaces. Friction coefficients were slightly lower for pull-out than for push-in for the porous surfaces but not for the polished surface. The effective interference was as little as 30% of the nominal interference for the porous surfaces. The determined variation in friction coefficient with radial contact force, as well as the loss of interference during implantation will enable a more accurate representation of implant press-fitting for simulations. [ABSTRACT FROM AUTHOR]

Meniscal implants have been developed in an attempt to provide pain relief and prevent pathological degeneration of articular cartilage. However, as yet there has been no systematic and comprehensive analysis of the effects of the meniscal design variables on meniscal function across a wide patient population, and there are no clear design criteria to ensure the functional performance of candidate meniscal implants. Our aim was to develop a statistically-augmented, experimentally-validated, computational platform to assess the effect of meniscal properties and patient variables on knee joint contact mechanics during the activity of walking. Our analysis used Finite Element Models (FEMs) that represented the geometry, kinematics as based on simulated gait and contact mechanics of three laboratory tested human cadaveric knees. The FEMs were subsequently programmed to represent prescribed meniscal variables (circumferential and radial/axial moduli--Ecm, Erm, stiffness of the meniscal attachments--Slpma, Slamp) and patient variables (varus/valgus alignment--VVA, and articular cartilage modulus--Ec). The contact mechanics data generated from the FEM runs were used as training data to a statistical interpolator which estimated joint contact data for untested configurations of input variables. Our data suggested that while Ecm and Erm of a meniscus are critical in determining knee joint mechanics in early and late stance (peak 1 and peak 3 of the gait cycle), for some knees that have greater laxity in the mid-stance phase of gait, the stiffness of the articular cartilage, Ec, can influence force distribution across the tibial plateau. We found that the medial meniscus plays a dominant load-carrying role in the early stance phase and less so in late stance, while the lateral meniscus distributes load throughout gait. Joint contact mechanics in the medial compartment are more sensitive to Ecm than those in the lateral compartment. Finally, throughout stance, varus-valgus alignment can overwhelm these relationships while the stiffness of meniscal attachments in the range studied have minimal effects on the knee joint mechanics. In summary, our statistically-augmented, computational platform allowed us to study how meniscal implant design variables (which can be controlled at the time of manufacture or implantation) interact with patient variables (which can be set in FEMs but cannot be controlled in patient studies) to affect joint contact mechanics during the activity of simulated walking. [ABSTRACT FROM AUTHOR]

Since meniscal geometry affects the cartilage contact pressures, it is essential to carefully define the geometry of the synthetic meniscal implant that we developed. Recently, six independent modes of sizeand shape-related geometry variation were identified through 3D statistical shape modeling (SSM) of the medial meniscus. However, this model did not provide information on the functional importance of these geometry characteristics. Therefore, in this study finite element simulations were performed to determine the influence of anatomically-based meniscal implant size and shape variations on knee cartilage contact pressures. Finite element simulations of the knee joint were performed for a total medial meniscectomy, an allograft, the average implant geometry, six implant sizes and ten shape variations. The geometries of the allograft and all implant variations were based on the meniscus SSM. Cartilage contact pressures and implant tensile strains were evaluated in full extension under 1200 N of axial compression. The average implant induced cartilage peak pressures intermediate between the allograft and meniscectomy and also reduced the cartilage area subjected to pressures >5 MPa compared to the meniscectomy. The smaller implant sizes resulted in lower cartilage peak pressures and compressive strains than the allograft, yet high implant tensile strains were observed. Shape modes 2, 3 and 6 affected the cartilage contact stresses but to a lesser extent than the size variations. Shape modes 4 and 5 did not result in changes of the cartilage stress levels. The present study indicates that cartilage contact mechanics are more sensitive to implant size than to implant shape. Down-sizing the implant resulted in more favorable contact mechanics, but caused excessive material strains. Further evaluations are necessary to balance cartilage contact pressures and material strains to ensure cartilage protection and longevity of the implant. [ABSTRACT FROM AUTHOR]

The contact mechanics of artificial metal-on-polyethylene hip joints are believed to affect the lubrication, wear and friction of the articulating surfaces and may lead to the joint loosening. Finite element analysis has been widely used for contact mechanics studies and good agreements have been achieved with current experimental data; however, most studies were carried out with idealist spherical geometries of the hip prostheses rather than the realistic worn surfaces, either for simplification reason or lacking of worn surface profile. In this study, the worn surfaces of the samples from various stages of hip simulator testing (0 to 5 million cycles) were reconstructed as solid models and were applied in the contact mechanics study. The simulator testing results suggested that the center of the head has various departure value from that of the cup and the value of the departure varies with progressively increased wear. This finding was adopted into the finite element study for better evaluation accuracy. Results indicated that the realistic model provided different evaluation from that of the ideal spherical model. Moreover, with the progressively increased wear, large increase of the contact pressure (from 12 to 31 MPa) was predicted on the articulating surface, and the predicted maximum von Mises stress was increased from 7.47 to 13.26 MPa, indicating the marked effect of the worn surface profiles on the contact mechanics of the joint. This study seeks to emphasize the importance of realistic worn surface profile of the acetabular cup especially following large wear volume. [ABSTRACT FROM AUTHOR]

While correction of dysplastic acetabular deformity has been a focus of both clinical treatment and research, concurrent femoral deformities have only more recently received serious attention. The purpose of this study was to determine how including abnormalities in femoral head-neck offset and femoral version alter computationally derived contact stresses in patients with combined dysplasia and femoroacetabular impingement (FAI). Hip models with patient-specific bony anatomy were created from preoperative and postoperative CT scans of 20 hips treated with periacetabular osteotomy and femoral osteochondroplasty. To simulate performing only a PAO, a third model was created combining each patient's postoperative pelvis and preoperative femur geometry. These three models were initialized with the femur in two starting orientations: (1) standardized template orientation, and (2) using patient-specific anatomic landmarks. Hip contact stresses were computed in all 6 model sets during an average dysplastic gait cycle, an average FAI gait cycle, and an average stand-to-sit activity using discrete element analysis. No significant differences in peak contact stress (p = 0.190 to 1), mean contact stress (p = 0.273 to 1), or mean contact area (p = 0.050 to 1) were identified during any loading activity based on femoral alignment technique or inclusion of femoral osteochondroplasty. These findings suggest that presence of abnormal femoral version and/or head-neck offset deformities are not themselves predominant factors in intra-articular contact mechanics during gait and stand-to-sit activities. Inclusion of modified movement patterns caused by these femoral deformities may be necessary for models to adequately capture the mechanical effects of these clinically recognized risk factors for negative outcomes. [ABSTRACT FROM AUTHOR]

Edge loading can negatively impact the biomechanics and long-term performance of hip replacements. Although edge loading has been widely investigated for hard-on-hard articulations, limited work has been conducted for hard-on-soft combinations. The aim of the present study was to investigate edge loading and its effect on the contact mechanics of a modular metal-on-polyethylene (MoP) total hip replacement (THR). A three-dimensional finite element model was developed based on a modular MoP bearing. Different cup inclination angles and head lateral microseparation were modelled and their effect on the contact mechanics of the modular MoP hip replacement were examined. The results showed that lateral microseparation caused loading of the head on the rim of the cup, which produced substantial increases in the maximum von Mises stress in the polyethylene liner and the maximum contact pressure on both the articulating surface and backside surface of the liner. Plastic deformation of the liner was observed under both standard conditions and microseparation conditions, however, the maximum equivalent plastic strain in the liner under microseparation conditions of 2000 mm was predicted to be approximately six times that under standard conditions. The study has indicated that correct positioning the components to avoid edge loading is likely to be important clinically even for hard-on-soft bearings for THR. [ABSTRACT FROM AUTHOR]

The 16-station RandomPOD wear test system, previously validated for prosthetic hip wear, was used in the simulation of knee wear mechanisms with a ball-on-flat test configuration. This consisted of a CoCr pin with a ground and polished spherical bearing surface (radius 28 mm) against a conventional, gamma-sterilized UHMWPE disk in serum lubrication. The biaxial motion, consisting of x and y translations, and the load was non-cyclic. Relative to the disk, the center of contact wandered within a circle of 10 mm diameter, and the average sliding velocity was 15.5 mm/s (ranging from 0 to 31 mm/s). The load varied non-cyclically between 0 and 142 N (average 73 N). In the 60-day test with 16 similar wear couples, moderate adhesive wear, the principal wear mechanism of a well-functioning prosthetic knee, dominated. This showed as a burnished, circular wear mark (diameter 13.2 mm, area 137 mm² ). The wear factor was 2.04 ± 0.03 x 10-6mm³/Nm (mean ± 95 percent confidence limit). For the first time a truly multidirectional, realistic and uniform, large capacity pin-on-disk simulation of knee wear mechanisms was implemented. [ABSTRACT FROM AUTHOR]

Despite significant advances in scaffold design, manufacture, and development, it remains unclear what forces these scaffolds must withstand when implanted into the heavily loaded environment of the knee joint. The objective of this study was to fully quantify the dynamic contact mechanics across the tibial plateau of the human knee joint during gait and stair climbing. Our model consisted of a modified Stanmore knee simulator (to apply multi-directional dynamic forces), a two-camera motion capture system (to record joint kinematics), an electronic sensor (to record contact stresses on the tibial plateau), and a suite of post-processing algorithms. During gait, peak contact stresses on the medial plateau occurred in areas of cartilage-cartilage contact; while during stair climb, peak contact stresses were located in the posterior aspect of the plateau, under the meniscus. On the lateral plateau, during gait and in early stair-climb, peak contact stresses occurred under the meniscus, while in late stair-climb, peak contact stresses were experienced in the zone of cartilage-cartilage contact. At 45% of the gait cycle, and 20% and 48% of the stair-climb cycle, peak stresses were simultaneously experienced on both the medial and lateral compartment, suggesting that these phases of loading warrant particular consideration in any simulation intended to evaluate scaffold performance. Our study suggests that in order to design a scaffold capable of restoring 'normal' contact mechanics to the injured knees, the mechanics of the intended site of implantation should be taken into account in any pre-clinical testing regime. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

For nearly a century, articular cartilage has been known for its exceptional tribological properties. For nearly as long, there have been research efforts to elucidate the responsible mechanisms for application toward biomimetic bearing applications. It is now widely accepted that interstitial fluid pressurization is the primary mechanism responsible for the unusual lubrication and load bearing properties of cartilage. Although the biomechanics community has developed elegant mathematical theories describing the coupling of solid and fluid (biphasic) mechanics and its role in interstitial lubrication, quantitative gaps in our understanding of cartilage tribology have inhibited our ability to predict how tribological conditions and material properties impact tissue function. This paper presents an analytical model of the interstitial lubrication of biphasic materials under migrating contact conditions. Although finite element and other numerical models of cartilage mechanics exist, they typically neglect the important role of the collagen network and are limited to a specific set of input conditions, which limits general applicability. The simplified approach taken in this work aims to capture the broader underlying physics as a starting point for further model development. In agreement with existing literature, the model indicates that a large Peclet number, Pe, is necessary for effective interstitial lubrication. It also predicts that the tensile modulus must be large relative to the compressive modulus. This explains why hydrogels and other biphasic materials do not provide significant interstitial pressure under high Pe conditions. The model quantitatively agrees with in-situ measurements of interstitial load support and the results have interesting implications for tissue engineering and osteoarthritis problems. This paper suggests that a low tensile modulus (from chondromalacia or local collagen rupture after impact, for example) may disrupt interstitial pressurization, increase shear stresses, and activate a condition of progressive surface damage as a potential precursor of osteoarthritis. [ABSTRACT FROM AUTHOR]

This short communication provides details on customized Tekscan Analysis Programs (TAP) which extract comprehensive contact mechanics metrics from piezoelectric sensors in articulating joints across repeated loading cycles. The code provides functionality to identify regions of interest (ROI), compute contact mechanic metrics, and compare contact mechanics across multiple test conditions or knees. Further, the variability of identifying ROIs was quantified between seven different users and compared to an expert. Overall, the contribution of four variables were studied: two knee specimens; two points in the gait cycle; two averaging methods; and seven observers, to determine if variations in these values played a role in accurately quantifying the ROI. The relative error between the force ratio from each observer's ROI and the expert ROI was calculated as the output of interest. A multivariate linear mixed effects model was fit to the four variables for the relative error with an observer- and knee-specific random intercept. Results from the fitted model showed a statistically significant difference at the 0.05 level in the mean relative errors at the two gait points. Additionally, variability in the relative errors attributed to the observer, knee, and random errors was quantified. To reduce variability amongst users, by ensuring low inter-observer variability and increasing segmentation accuracy of knee contact mechanics, a training module and manual have been included as supplemental material. By sharing this code and training manual, we envisage that it can be used and modified to analyze outputs from a range of sensors, joints, and test conditions. [ABSTRACT FROM AUTHOR]

It is important to study joint contact mechanics to better understand the processes which lead to cartilage degradation. The purpose of this study was to develop and validate a finite element (FE) model of a human elbow capable of predicting joint contact area and stress. A cylindrical constrained elbow joint loading apparatus was used to measure the cartilage compression and contact area for a single cadaveric specimen. A computer model of the same joint was created based on computed tomography images of the specimen, and the same loading was simulated using FE contact analysis. The model-predicted joint compression and contact area corresponded closely with experiment-measured results (differences of -4.9% and +9.6%). A sensitivity analysis showed that the model results were sensitive to cartilage and bone material properties, as well as the cartilage thickness distribution. The results of this study underline the importance of using accurate material properties and physiological cartilage thickness distributions when simulating cartilage contact mechanics. [ABSTRACT FROM AUTHOR]

Computational models of the natural hip joint are needed to examine and optimise tissue sparing interventions where the natural cartilage remains part of the bearing surfaces. Although the importance of interstitial fluid pressurisation in the performance of cartilage has long been recognized, few studies have investigated the time dependent interstitial fluid pressurisation in a three dimensional natural hip joint model. The primary aim of this study was to develop a finite element model of the natural hip incorporating the biphasic cartilage layers that was capable of simulating the joint response over a prolonged physiological loading period. An initial set of sensitivity studies were also undertaken to investigate the influence of hip size, clearance, cartilage properties, thickness and hemiarthroplasty on the contact mechanics of the joint. The contact stress, contact area, fluid pressure and fluid support ratio were calculated and cross-compared between models with different parameters to evaluate their influence. It was found that the model predictions for the period soon after loading were sensitive to the hip size, clearance, cartilage aggregate modulus, thickness and hemiarthroplasty, while the time dependent behaviour over 3000 s was influenced by the hip clearance and cartilage aggregate modulus, permeability, thickness and hemiarthroplasty. The modelling methods developed in this study provide a basic platform for biphasic simulation of the whole hip joint onto which more sophisticated material models or other input parameters could be added in the future. [ABSTRACT FROM AUTHOR]

Abstract: Knee joint forces measured from instrumented implants provide important information for testing the validity of computational models that predict knee joint forces. The purpose of this study was to validate a parametric numerical model for predicting knee joint contact forces against measurements from four subjects with instrumented TKRs during the stance phase of gait. Model sensitivity to abnormal gait patterns was also investigated. The results demonstrated good agreement for three subjects with relatively normal gait patterns, where the difference between the mean measured and calculated forces ranged from 0.05 to 0.45 body weights, and the envelopes of measured and calculated forces (from three walking trials) overlapped. The fourth subject, who had a “quadriceps avoidance” external moment pattern, initially had little overlap between the measured and calculated force envelopes. When additional constraints were added, tailored to the subject’s gait pattern, the model predictions improved to complete force envelope overlap. Coefficient of multiple determination analysis indicated that the shape of the measured and calculated force waveforms were similar for all subjects (adjusted coefficient of multiple correlation values between 0.88 and 0.92). The parametric model was accurate in predicting both the magnitude and waveform of the contact force, and the accuracy of model predictions was affected by deviations from normal gait patterns. Equally important, the envelope of forces generated by the range of solutions substantially overlapped with the corresponding measured envelope from multiple gait trials for a given subject, suggesting that the variable strategic processes of in vivo force generation are covered by the solution range of this parametric model. [Copyright &y& Elsevier]

Abstract: Interstitial fluid pressurization, a consequence of a biphasic tissue structure, is essential to the load bearing and lubrication properties of articular cartilage. Focal tissue degradation may interfere with this protective mechanism, eventually leading to gross degeneration and osteoarthritis. Our long-term goal is to determine whether local contacts can be used as a means to probe local tissue integrity and functionality. In the present work, Hertzian rate-controlled microindentation was used as a model of the more complicated sliding system to directly determine the effects of contact radius and deformation rate on interstitial load support. During localized contact between a steel spherical probe and bovine articular cartilage, the equilibrium and non-equilibrium responses were well-fit by the Hertz model (R 2>0.998) with a mean equilibrium contact modulus of 0.93MPa. The effective contact modulus and fluid load fraction were independent of indentation depth, contact radius, and normal force; both increased monotonically with indentation rate. At 21μm/s indentation rate, the cartilage was effectively stiffened by 6-fold with the fluid pressure supporting 85% of the contact force. The results motivated a simple analytical model that directly links the tribomechanical response (including fluid load support) and the Peclet number to measurable material properties and controllable experimental variables. This paper demonstrates that tribological contacts can be used to probe local functional properties. Such measurements can add important insights into the roles of focal tissue damage and impaired local functionality in the pathogenesis of osteoarthritis. [Copyright &y& Elsevier]

Abstract: A finite element formulation of neutral solute transport across a contact interface between deformable porous media is implemented and validated against analytical solutions. By reducing the integral statements of external virtual work on the two contacting surfaces into a single contact integral, the algorithm automatically enforces continuity of solute molar flux across the contact interface, whereas continuity of the effective solute concentration (a measure of the solute mechano-chemical potential) is achieved using a penalty method. This novel formulation facilitates the analysis of problems in biomechanics where the transport of metabolites across contact interfaces of deformable tissues may be of interest. This contact algorithm is the first to address solute transport across deformable interfaces, and is made available in the public domain, open-source finite element code FEBio (http://www.febio.org). [Copyright &y& Elsevier]

Abstract: 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. [Copyright &y& Elsevier]

Abstract: The human hip joint withstands high contact forces during daily activity and is therefore susceptible to injury and structural deterioration over time. Knowledge of muscle-force contributions to hip joint loading may assist in the development of strategies to prevent and manage conditions such as osteoarthritis, femoro-acetabular impingement and fracture. The main aim of this study was to determine the contributions of individual muscles to hip contact force in normal walking. Muscle contributions to hip contact force were calculated based on a previously published dynamic optimization solution for normal walking, which provided the time histories of joint motion, ground reaction forces, and muscle forces during the stance and swing phases of gait. The force developed by each muscle plus its contribution to the ground reaction force were used to determine the muscle’s contribution to hip contact force. Muscles were the major contributors to hip contact force, with gravitational and centrifugal forces combined contributing less than 5% of the total contact force. Four muscles that span the hip – gluteus medius, gluteus maximus, iliopsoas, and hamstrings – contributed most significantly to the three components of the hip contact force and hip contact impulse (integral of hip contact force over time). Three muscles that do not span the hip – vasti, soleus, and gastrocnemius – also contributed substantially to hip joint loading. These results provide additional insight into lower-limb muscle function during walking and may also be relevant to studies of cartilage degeneration and bone remodelling at the hip. [Copyright &y& Elsevier]

Abstract: Falls initiated by slips and trips are a serious health hazard to older adults. Experimental studies have provided important descriptions of postural responses to slipping, but it is difficult to determine why some slips result in falls from experiments alone. Computational modeling and simulation techniques can complement experimental approaches by identifying causes of failed recovery attempts. The purpose of this study was to develop a method to determine the impact of a systematic reduction in the foot–floor friction coefficient (μ) on the kinematics of walking shortly after heel contact (∼200ms). A walking model that included foot–floor interactions was utilized to find the set of moments that best tracked the joint angles and measured ground reaction forces obtained from a non-slipping (dry) trial. A “passive” slip was simulated by driving the model with the joint-moments from the dry simulation and by reducing μ. Slip simulations with values of μ greater than the subject-specific peak required coefficient of friction (RCOF), an experimental measure of slip-resistant gait, resulted in only minor deviations in gait kinematics from the dry condition. In contrast, slip simulations run in environments characterized by μ

Abstract: Musculoskeletal models are currently the primary means for estimating in vivo muscle and contact forces in the knee during gait. These models typically couple a dynamic skeletal model with individual muscle models but rarely include articular contact models due to their high computational cost. This study evaluates a novel method for predicting muscle and contact forces simultaneously in the knee during gait. The method utilizes a 12 degree-of-freedom knee model (femur, tibia, and patella) combining muscle, articular contact, and dynamic skeletal models. Eight static optimization problems were formulated using two cost functions (one based on muscle activations and one based on contact forces) and four constraints sets (each composed of different combinations of inverse dynamic loads). The estimated muscle and contact forces were evaluated using in vivo tibial contact force data collected from a patient with a force-measuring knee implant. When the eight optimization problems were solved with added constraints to match the in vivo contact force measurements, root-mean-square errors in predicted contact forces were less than 10N. Furthermore, muscle and patellar contact forces predicted by the two cost functions became more similar as more inverse dynamic loads were used as constraints. When the contact force constraints were removed, estimated medial contact forces were similar and lateral contact forces lower in magnitude compared to measured contact forces, with estimated muscle forces being sensitive and estimated patellar contact forces relatively insensitive to the choice of cost function and constraint set. These results suggest that optimization problem formulation coupled with knee model complexity can significantly affect predicted muscle and contact forces in the knee during gait. Further research using a complete lower limb model is needed to assess the importance of this finding to the muscle and contact force estimation process. [Copyright &y& Elsevier]

Abstract: Diameter and diametral clearance of the bearing surfaces of metal-on-metal hip implants and structural supports have been recognised as key factors to reduce the dry contact and hydrodynamic pressures and improve lubrication performance. On the other hand, application of aspherical bearing surfaces can also significantly affect the contact mechanics and lubrication performance by changing the radius of the curvature of a bearing surface and consequently improving the conformity between the head and the cup. In this study, a novel metal-on-metal hip implant employing a specific aspherical bearing surface, Alpharabola, as the acetabular surface was investigated for both contact mechanics and elastohydrodynamic lubrication under steady-state conditions. When compared with conventional spherical bearing surfaces, a more uniform pressure distribution and a thicker lubricant film thickness within the loaded conjunction were predicted for this novel Alpharabola hip implant. The effects of the geometric parameters of this novel acetabular surface on the pressure distribution and lubricant thickness were investigated. A significant increase in the predicted lubricant film thickness and a significant decrease in the dry contact and hydrodynamic pressures were found with appropriate combinations of these geometric parameters, compared with the spherical bearing surface. [Copyright &y& Elsevier]

Abstract: The knowledge of articular cartilage contact biomechanics in the knee joint is important for understanding the joint function and cartilage pathology. However, the in vivo tibiofemoral articular cartilage contact biomechanics during gait remains unknown. The objective of this study was to determine the in vivo tibiofemoral cartilage contact biomechanics during the stance phase of treadmill gait. Eight healthy knees were magnetic resonance (MR) scanned and imaged with a dual fluoroscopic system during gait on a treadmill. The tibia, femur and associated cartilage were constructed from the MR images and combined with the dual fluoroscopic images to determine in vivo cartilage contact deformation during the stance phase of gait. Throughout the stance phase of gait, the magnitude of peak compartmental contact deformation ranged between 7% and 23% of the resting cartilage thickness and occurred at regions with thicker cartilage. Its excursions in the anteroposterior direction were greater in the medial tibiofemoral compartment as compared to those in the lateral compartment. The contact areas throughout the stance phase were greater in the medial compartment than in the lateral compartment. The information on in vivo tibiofemoral cartilage contact biomechanics during gait could be used to provide physiological boundaries for in vitro testing of cartilage. Also, the data on location and magnitude of deformation among non-diseased knees during gait could identify where loading and later injury might occur in diseased knees. [Copyright &y& Elsevier]

Abstract: Patellofemoral pain syndrome (PFPS) is a disorder of the patellofemoral (PF) joint in which abnormal tracking is often cited as a factor in pain development. PF tracking is partially dependent on passive stabilizers (ex: PF geometry). Relations amongst PFPS, PF tracking, and contact mechanics are poorly understood. In-vivo investigation of passive PF joint stabilizers including PF tracking, contact mechanics, cartilage thickness, and patellar shape will allow structural characterization of the PF joint and may highlight differences associated with PFPS. This study examined the role that passive stabilizers play in PFPS (n=10) versus healthy subjects (n=10). PF tracking (contact area centroid migration), cartilage thickness, shape, congruence, and contact patterns were quantified using magnetic resonance imaging during isometric loading at 15°, 30°, and 45° of knee flexion. Distinct relationships were identified between patellar shape and tracking and contact, particularly at low flexion (15–30°). Healthy subjects exhibited distinct PF tracking and contact patterns related to Type I patella shape (80%) with increasing total contact area (p<0.001) and proximal centroid migration (15–30°p=0.012; 30–45°p<0.001) for increasing knee angles. PFPS subjects deviated from these patterns at low flexion, demonstrating higher total contact area than healthy subjects (p=0.046 at 15°), lack of proximal centroid migration (15–30°), and more Type II (30%) and III (20%) patella shapes. This study highlights a new finding that patellar shape combined with low degrees of flexion (15–30°) may be important to consider, as this is where PFPS tracking and contact patterns deviate from healthy. [Copyright &y& Elsevier]

Abstract: Hemiarthroplasty is a potentially attractive alternative to knee replacement for young, active patients, as it allows preservation of more bone stock for potential revisions. However, there has been limited success with hemiarthroplasty or spacers to date. The wear and degradation of the biomaterial–cartilage interface is of paramount importance in the design and success of hemiarthroplasties. A comprehensive understanding of the tribological performance of hemiarthroplasty implants in the natural joint is required. The objective of this study was to investigate the tribological response of bovine medial compartmental knees, both natural and hemiarthroplasty replaced, under physiological loads and motion. The conformity of these metallic hemiarthroplasties was varied (conforming plates with radius of 50mm and radius of 100mm and a flat plate design), in order to examine the effects of conformity and contact stress, on the friction, friction shear stress and cartilage degeneration. With decreasing conformity of hemiarthroplasty bearings, an increase in contact stress was found, which resulted in elevated friction, elevated friction shear stress and increased cartilage degeneration. A strong correlation was found between contact stress and wear and between friction shear stress and wear. This new and unique in vitro tribological simulation has shown the direct elevation of friction, surface fibrillation and biomechanical wear of cartilage, upon replacing the tibia with a hemiarthroplasty, particularly when using low conformity hemiarthroplasty designs. [Copyright &y& Elsevier]

Abstract: The objective of this study was to provide a detailed experimental assessment of the two-dimensional cartilage strain distribution on the cross-section of the human patellofemoral joint (PFJ) subjected to physiological load magnitudes and rates. The medial side of six human PFJs sectioned along their mid-sagittal plane was loaded up to the equivalent of two body weights on a whole joint, and strain measurements obtained from digital image correlation are reported at 0.5s. Normal strains tangential to the articular surface and shear strains in the plane of the cross-section showed consistent patterns among all specimens, whereas normal strains perpendicular to the articular surface exhibited some variability that may be attributed to subject-specific variations in material properties through the depth of the articular layers. Elevated tensile and compressive principal normal strains were observed near the articular surface, around the center of the contact region, with additional locations of elevated compressive strains occurring at the bone–cartilage interface. Under an average contact stress of ∼3.3MPa, the peak compressive principal normal strains for the patella and femur averaged −0.158±0.072 and −0.118±0.051, respectively, magnitudes that are significantly greater than the relative changes in cartilage thickness, −0.090±0.030 and −0.072±0.038 (p<0.005). These experimental results provide a detailed description of the manner by which human PFJ articular layers deform in situ under physiological load conditions. [Copyright &y& Elsevier]

Abstract: 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. [Copyright &y& Elsevier]

Abstract: In vivo knee contact forces are difficult to determine using numerical methods because there are more unknown forces than equilibrium equations available. We developed parametric methods for computing contact forces across the knee joint during the stance phase of level walking. Three-dimensional contact forces were calculated at two points of contact between the tibia and the femur, one on the lateral aspect of the tibial plateau, and one on the medial side. Muscle activations were parametrically varied over their physiologic range resulting in a solution space of contact forces. The obtained solution space was reasonably small and the resulting force pattern compared well to a previous model from the literature for kinematics and external kinetics from the same patient. Peak forces of the parametric model and the previous model were similar for the first half of the stance phase, but differed for the second half. The previous model did not take into account the transverse external moment about the knee and could not calculate muscle activation levels. Ultimately, the parametric model will result in more accurate contact force inputs for total knee simulators, as current inputs are not generally based on kinematics and kinetics inputs from TKR patients. [Copyright &y& Elsevier]

Abstract: Vibrotactile thresholds at the fingertips are affected by a number of individual, environmental, and testing factors. In the current study, we theoretically analyzed the effects of the contact orientation of the probe on the fingertip and the static pre-indentation on the dynamic deformation of the soft tissues of the fingertip in the vibrotactile tests using a nonlinear finite element model. The fingertip considered in the 3D finite element model is the distal phalanx, the portion from the distal end to the distal interphalangeal (DIP) joint articulation. The fingertip is contacted by the probe at four different contact locations, which are regulated by contact angles (, , , and ), and three different pre-indentations (0.5, 1.0, and 1.5mm). The model predictions indicated that the average spatial summation of the vibration displacement (SVD) at the fingertip depends on the static pre-indentation and the probe/indentor contact orientation; although the resonance characteristics of the fingertip are not affected by either the pre-indentation or the contact location. The location-dependence of the vibration exposure factors at the fingertip was found to increase with increasing static pre-indentation. At a static indentation of 1.5mm, the test condition specified in the ISO-13091-1 standard, the values of the SVDs determined at different probe/fingertip contact orientations differ as much as 125%. Since the dynamic displacements of the soft tissues are believed to affect the vibrotactile threshold, the current results suggest that the contact orientation of the probe on the fingertip should be strictly defined and restricted to obtain reliable results in the vibrotactile perception threshold tests. [Copyright &y& Elsevier]

Abstract: Mammography is currently the most widely used screening and diagnostic tool for breast cancer. Because X-ray images are 2D projections of a 3D object, it is not trivial to localise features identified in mammogram pairs within the breast volume. Furthermore, mammograms represent highly deformed configurations of the breast under compression, thus the tumour localisation process relies on the clinician''s experience. Biomechanical models of the breast undergoing mammographic compressions have been developed to overcome this limitation. In this study, we present the development of a modelling framework that implements Coulomb''s frictional law with a finite element analysis using a -continuous Hermite mesh. We compared two methods of this contact mechanics implementation: the penalty method, and the augmented Lagrangian method, the latter of which is more accurate but computationally more expensive compared to the former. Simulation results were compared with experimental data from a soft silicon gel phantom in order to evaluate the modelling accuracy of each method. Both methods resulted in surface-deformation root-mean-square errors of less than 2mm, whilst the maximum internal marker prediction error was less than 3mm when simulating two mammographic-like compressions. Simulation results were confirmed using the augmented Lagrangian method, which provided similar accuracy. We conclude that contact mechanics on soft elastic materials using the penalty method with an appropriate choice of the penalty parameters provides sufficient accuracy (with contact constraints suitably enforced), and may thus be useful for tracking breast tumours between clinical images. [Copyright &y& Elsevier]

Abstract: Computational models have recently been developed to replicate experimental conditions present in the Stanmore knee wear simulator. These finite element (FE) models, which provide a virtual platform to evaluate total knee replacement (TKR) mechanics, were validated through comparisons with experimental data for a specific implant. As with any experiment, a small amount of variability is inherently present in component alignment, loading, and environmental conditions, but this variability has not been previously incorporated in the computational models. The objectives of the current research were to assess the impact of experimental variability on predicted TKR mechanics by determining the potential envelope of joint kinematics and contact mechanics present during wear simulator loading, and to evaluate the sensitivity of the joint mechanics to the experimental parameters. In this study, 8 component alignment and 4 experimental parameters were represented as distributions and used with probabilistic methods to assess the response of the system, including interaction effects. The probabilistic FE model evaluated two levels of parameter variability (with standard deviations of component alignment parameters up to 0.5mm and 1°) and predicted a variability of up to 226% (3.44mm) in resulting anterior–posterior (AP) translation, up to 169% (4.30°) in internal–external (IE) rotation, but less than 10% (1.66MPa) in peak contact pressure. The critical alignment parameters were the tilt of the tibial insert and the IE rotational alignment of the femoral component. The observed variability in kinematics and, to a lesser extent, contact pressure, has the potential to impact wear observed experimentally. [Copyright &y& Elsevier]

Joint kinematics and contact mechanics dictate the success of current total knee replacement (TKR) devices. Efficient computer models present an effective way of evaluating these characteristics. Predicted contact stress and area due to articulations at the tibio-femoral and patello-femoral interfaces indicate potential clinical performance. Previous finite element (FE) knee models have generally been used to predict contact stresses and/or areas during static or quasi-static loading conditions. Explicit dynamic FE analyses have recently been used to efficiently predict TKR kinematics and contact mechanics during dynamic loading conditions. The objective of this study was to develop and experimentally validate an explicit FE TKR model that incorporates tibio-femoral and patello-femoral articulations. For computational efficiency, we developed rigid body analyses that can reasonably reproduce the kinematics, contact pressure distribution, and contact area of a fully deformable system.Results from the deformable model showed that the patello-femoral and tibio-femoral kinematics were in good agreement with experimental knee simulator measurements. Kinematic results from the rigid body analyses were nearly identical to those from the fully deformable model, and the contact pressure and contact area correlation was acceptable given the great reduction in analysis time. Component mesh density studied had little effect on the predicted kinematics, particularly for the patellar component, and small effects on the predicted contact pressure and area. These analyses have shown that, at low computational cost, a force-control dynamic simulation of a gait cycle can yield useful and predictable results. [Copyright &y& Elsevier]