Summary Statement. Computational models are the primary tools for efficient design-phase exploration of knee replacement concepts before in vitro testing. To improve design-phase efficiency, a subject-specific computational platform was developed that allows designers to assess devices in realistic conditions by directly integrating subject-specific experimental data in these models. Introduction. Early in the design-phase of new implant design, numerous in vitro tests would be desirable to assess the influence of design parameters or component alignment on the performance of the device. However, cadaveric testing of knee replacement devices is a costly and time-consuming procedure, requiring manufacture of parts, preparation of cadaveric specimens, and personnel to carry of the experiments. Validated computational models are ideally suited for pre-clinical, high-volume design evaluation. Initial development of these models requires substantial time and expertise; once developed, however, computational simulations may be applied for comparative evaluation of devices in an extremely efficient manner [Baldwin et al. 2012]. Still, computational models are complementary of experimental testing and for this reason, computational models tuned with subject-specific experimental data, e.g. soft tissue parameters, could bring even more efficiency in the design phase. The objective of the current study was to develop a platform of tools that easily allows for subject-specific knee simulations. The system integrates with commercially available medical imaging and finite element software to allow for direct, efficient comparison of designs and surgical alignment under a host of different boundary conditions. Patients & Methods. MRI image was acquired, and 3D bone models were generated using the Mimics Innovation Suite® (Materialise NV, Leuven, Belgium). The two models (1) tibiofemoral (TF) joint laxity including ligamentous constraint and (2) whole joint (TF and patellofemoral (PF)) mechanics during dynamic activities of daily living (e.g. gait, squat, chair-rise), developed in Abaqus/Explicit (SIMULIA, Providence, RI), were then be adapted with integrated subject-specific attachment sites. Results. The suite of tools provides a platform for baseline evaluation of design factors, comparison of new implant designs with predicate devices, and assessment of robustness to surgical alignment. This platform is currently capable of taking into account subject-specific factors in order to provide realistic results in relation with experimental data. Implant material properties, ligament properties and initial conditions can be varied, and results compared, to evaluate the influence of a host of design and surgical factors on implant performance. The interface allows users without complex finite element expertise to setup, analyze and compare devices and interpret results. Discussion/Conclusion. A platform which allows implant designers to evaluate their design ideas in realistic conditions integrating subject-specific parameters and to compare with predicate devices has the potential to substantially decrease the development time for new devices. Designers can perform iterative modification to their devices to focus on an optimal design solution prior to in vitro testing, reducing the number of pre-clinical cadaveric experiments that may be required, and ultimately improving TKR mechanics in the patient population

Periprosthetic femoral fractures can occur as a complication of total hip arthroplasty and are often challenging to treat as the mechanical scenario is influenced by the presence of the metal prosthesis within the bone. This research focuses on finding the optimum fixation for transverse, Vancouver type B1 periprosthetic fractures, stabilised using locking plates and secured using screws. The aim of this study was to experimentally validate a computer model of a human femur, develop that model to represent a periprosthetic femoral fracture fixation and show how the model could be used to indicate differences between plating techniques.

In the first development stage, both a laboratory model and a finite element model were developed to evaluate the mechanical behaviour of an intact composite femur under axial loading. Axial strains were recorded along the medial length of the femur in both cases and compared to provide validation for the computational model predications. The computational intact femur model was then modified to include a cemented total hip replacement, and further adapted to include a periprosthetic fracture stabilised using a locking plate, with unicortical screws above, and bicortical screws below the transverse fracture.

For the intact femur case, the experimental and computational strain patterns correlated well with an average difference of 16%. Following the inclusion of the stem, there was a reduction in the strain in the region of the prosthesis reducing by an average of 45%. There was also a large increase in bulk stiffness with the introduction of the prosthesis. When the fracture and plate fixation were included, there was little difference in the proximal strain where the stem dominated, and the strains in the distal region were found to be highly sensitive to the distribution of the screws.

The results of this study indicate that screw configuration is an important factor in periprosthetic fracture fixation. A laboratory model of the periprosthetic facture case is now under development to further validate the computational models and the two approaches will then be used to determine optimum fixation methods for a range of clinical scenarios.

Introduction

Surgeons commonly resect additional distal femur during primary total knee arthroplasty (TKA) to correct a flexion contracture to restore range of motion and knee function. However, the effect of joint line elevation on the resulting TKA kinematics including frontal plane laxity is unclear. Thus, our goal was to quantify the effect of additional distal femoral resection on passive extension and mid-flexion laxity.

Introduction

Surgeons commonly resect additional distal femur during primary total knee arthroplasty (TKA) to correct a flexion contracture. However, the effect of joint line proximalization on TKA kinematics is unclear. Thus, our goal was to quantify the effect of additional distal femoral resection on knee extension and mid-flexion laxity.

Introduction: The complexity of the spine has made a complete understanding of its mechanical function difficult. As a consequence, biomechanical models have been used to describe the behaviour of the spine and its various components. A comprehensive mathematical model of the muscles of the lumbar spine and trunk is presented to enable computation of the forces and moments experienced by the lumbar intervertebral joints during physiological activities.

Methods: The model includes the nine major muscles crossing the region and concentrates on improving the estimated line of action for the muscles. The muscles are considered to consist of numerous fascicles, each with its own force producing potential based on size and line of action. The model respects the physical constraints imposed by the skeletal structure by ensuring that muscles maintain their anatomical position in various spinal postures. Validation was performed by comparing model predictions of maximum moments to published data from maximum isometric exertions in male volunteers. To highlight the potential novel uses of the model, three examples of muscle injury caused by surgical procedures were investigated; posterior lumbar surgery, impairment of abdominal muscles from anterior surgery and removal of the psoas major unilaterally during total hip replacement.

Results: The validation indicated that the model predicted forces similar to those measured in normal volunteers. The biomechanical changes resulting from the muscle injuries during the surgical procedures share several common features: decreased spinal compression and production of asymmetric moments during symmetric tasks.

Discussion: The results suggest that interference with muscles crossing or attaching to the lumbar spine can have a significant impact on its function.

INTRODUCTION: The complexity of the spine has made a complete understanding of its mechanical function difficult. As a consequence, biomechanical models have been used to describe the behaviour of the spine and its various components. A comprehensive mathematical model of the muscles of the lumbar spine and trunk is presented to enable computation of the forces and moments experienced by the lumbar intervertebral joints during physiological activities.

METHODS: The model includes the nine major muscles crossing the region and concentrates on improving the estimated line of action for the muscles. The muscles are considered to consist of numerous fascicles, each with its own force producing potential based on size and line of action. The model respects the physical constraints imposed by the skeletal structure by ensuring that muscles maintain their anatomical position in various spinal postures. Validation was performed by comparing model predictions of maximum moments to published data from maximum isometric exertions in male volunteers. To highlight the potential novel uses of the model, three examples of muscle injury caused by surgical procedures were investigated; posterior lumbar surgery, impairment of abdominal muscles from anterior surgery and removal of the psoas major unilaterally during total hip replacement.

RESULTS: The validation indicated that the model predicted forces similar to those measured in normal volunteers. The biomechanical changes resulting from the muscle injuries during the surgical procedures share several common features: decreased spinal compression and production of asymmetric moments during symmetric tasks.

DISCUSSION: The results suggest that interference with muscles crossing or attaching to the lumbar spine can have a significant impact on its function.

INTRODUCTION

Porous metallic materials, due to their capability of tailoring their mechanical properties to those of bone, have been suggested to be utilized in prosthesis to avoid the stress shielding phenomenon¹, believed to increase the risk of implant loosening².

The aim of this work is to obtain the most simplified model possible to simulate the mechanical behavior of a Ti6Al4V porous structure. For this purpose, a beam element model was analyzed and the results were then compared to a 3D-solid model.

The objective of this study was to compute the in vivo dynamic tibiofemoral contact forces for normal alignment, and then evaluate the change in contact forces and pressures with increasing varus-valgus and internal-external rotational malalignment of the femoral component. A three-dimensional computational model of the lower limb during deep knee bend was created using Kane’s method of dynamics. The change in forces from normal with malalignment of up to 10° valgus, 10° varus, 10° internal axial femoral rotation, and 10° internal axial femoral rotation were determined. In this study, varus-valgus malalignment had the greatest effect on medial-lateral pattelofemoral contact forces, with a maximum increase of 2.25 times body weight for 10° valgus malalignment. Axial malalignment had the greatest influence on tibiofemoral contact forces.

Numerous dynamic studies have evaluated the tibiofemoral contact pressures that follow total knee arthroplasty (TKA), and several static studies utilizing finite elements and pressure sensitive film have evaluated malalignment. The objective of this study was to compute the in vivo dynamic tibiofemoral contact forces for normal alignment and evaluate the change in contact pressure with increasing malalignment of the femoral component.

A three-dimensional computational model of the lower limb during deep flexion was created using Kane’s method of dynamics. A hybrid approach was used to determine the boundary conditions of the model. The motions of a total knee arthroplasty patient were measured using fluoroscopy. The motions of the patient were varied from the normal motions to simulate malalignment of the femoral component. The change in forces with malalignments of up to 10° valgus, 10° varus, 10° internal rotation, and 10° internal rotation were determined.

An increase in the axial tibiofemoral contact force from 2.44 times body weight (BW) to 2.62 BW and a decrease in the quadriceps force from 6.8 to 5.65 BW were observed with varus malalignment. The medial-lateral patellofemoral contact force decreased from 0.95 BW to 0.1 BW with 10° varus positioning of the femur and increased to 2.2 BW with 10° valgus positioning of the femur and a decrease in the patellar ligament forces from 1.70 to 1.63 BW was observed.

Changes in the tibiofemoral and patellofemoral forces of 1–2 BW were observed as the femur was malaligned with respect to the tibia. The most significant of these changes was the medial-lateral patellofemoral contact force. The implications of these findings are that malalignment could result in increased patellar subluxation or increased wear of the polyethylene component. Concerns were raised that this initial subject evaluated may not have had optimum alignment, thus leading to more optimal bearing surface stress conditions with varus malalignment. Future studies will be evaluated for subjects having the joint line restored to conditions for non-implanted knees.

Aims

Surgeons commonly resect additional distal femur during primary total knee arthroplasty (TKA) to correct a flexion contracture, which leads to femoral joint line elevation. There is a paucity of data describing the effect of joint line elevation on mid-flexion stability and knee kinematics. Thus, the goal of this study was to quantify the effect of joint line elevation on mid-flexion laxity.

Falls in adults are a major problem and can lead to injuries and death. In order to better understand falls and successful recoveries, identifying kinematics, kinetics, and muscle forces during recovery from loss of balance is crucial. To obtain reactive gait patterns, participants must be subjected to unexpected perturbations such as trips and slips. Previous researchers have reported kinetics recovery data following stumbling; however, the muscle force recovery patterns remain unknown. To better target exercises to reduce the risk of falls, we must first understand which muscles, their magnitude, and their coordination patterns, play a role in a successful recovery from a trip and a slip. Additionally, knowing the successful patterns of lower limb function can help with the diagnosis of faulty movements. A total of 20 healthy adults in their twenties with similar athletic backgrounds were perturbed on a split-belt treadmill using Computer-Assisted Rehabilitation Environment (Motkforce Link) at a preset speed of 1.1m/s. Two kinds of perturbations were administered: slip and trip. Slips were simulated by accelerating one belt, whereas trips were simulated by decelerating one belt. Both perturbations had similar intensity and only differed in the direction. Computational modeling was used to obtain lower-limb function during the compensatory step. SPM paired t-test was used to compare differences in recovery strategies between slip and trip through magnitude and patterns of joints. There were no significant differences in joint angles post tripping vs post-slipping. Results of net joint moments showed that compensating for the loss of balance due to tripping required a higher ankle plantarflexion moment than slipping (at 22-52%; 1.2± 0.3vs0.4±0.2, p<0.001). Additionally, larger gluteus maximus (at 40-50%;8.7±3.8vs2.7±1.1N/kg, p=0.001), gluteus medius (at23~33%; 22.6±5.7vs6.8±3.6N/kg, p<0.001) were generated than post-slipping, respectively. These findings suggested that greater GMAX and GMED forces are required post-trip recovery than slip. Future analysis of trip recovery showed the importance of ankle joint in recovering from forward and backward fall. These results can be used as references in remote diagnosis of joint and muscle weakness and assessment of the risk of falls with the use of accelerometers

Abstract. Objectives. A fibril reinforced multiphasic cartilage model was developed to improve the understanding of the depth-dependent cartilage internal structure and its through thickness biomechanical response. The heterogeneous model of cartilage was validated against full-field strain measurement obtained via Digital Image Correlation (DIC) during free swelling experiments. Methods. Hemi-cylindrical cartilage cores of 5mm diameter were obtained from porcine femoral condyles and humeral heads. The full field behaviour of these samples was monitored using DIC during an osmotic free swelling experiment performed following a standardised protocol [1]. Computational models were created in FEBio (version 2.8, . febio.org. ). The cartilage, submerged in saline solution was represented by a 1×1mm cube [2] with geometry and constrains set up to mimic the experimental conditions. Cartilage was modelled as a multiphasic material represented by one inhomogeneous layer with depth-dependent Young's modulus [3], zonally varied water content and zonally oriented collagen fibrils [4]. Experimental and predicted strain maps were compared to each other both qualitatively and quantitatively. Results. The numerical strain map showed high strain localisation close to the cartilage surface, with strain in this region reaching 40% and 12% for femoral and humeral samples respectively, this finding was confirmed in our experimental results. Strain magnitude gradually decreased with depth, reaching near-zero at around 200μm. This behaviour also matched experimental observations. Conclusions. Both sets of computational strain results exhibited very good agreement with experimental data, both in terms of cartilage through-thickness swelling behaviour and strain magnitude. Our results show the importance of including cartilage structural inhomogeneities and inclusions of collagen fibrils when simulating through-thickness cartilage swelling. These findings highlight the crucial role of collagen fibrils on both tissue solute transport properties as well as the overall biomechanical response of cartilage. Declaration of Interest. (b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project

Bone is a connective tissue that undergoes constant remodeling. Any disturbances during this process may result in undesired pathological conditions. A single nucleotide substitution (596T-A) in exon eight which leads to a M199K mutation in human RANKL was found to cause osteoclast-poor autosomal recessive osteopetrosis (ARO). Patients with ARO cannot be cured by hematopoietic stem cell transplantation and, without proper treatments, will die in their early age. To date, how this mutation alters RANKL function has not been characterized. We thus hypothesized that hRANKL M199 residue is a structural determinant for normal RANKL-RANK interaction and osteoclast differentiation. By sharing our findings, we aim to achieve an improved clinical outcome in treating bone-related diseases such as osteoporosis, ARO and osteoarthritis. Site-directed mutagenesis was employed to create three rat RANKL mutants, replacing the methionine 200 (human M199 equivalent residue) with either lysine (M200K), alanine (M200A) or glutamic acid (M200E). Recombinant proteins were subsequently purified through affinity chromatography and visualized by Coomassie blue staining and western blot. MTS was carried out before osteoclastogenesis assay in vitro to measure the cellular toxicity. Bone resorption pit assay, immuno-fluorescent staining, luciferase reporter assay, RT-PCR, western blot and calcium oscillation detection were also conducted to explore the biological effect of rRANKL mutants. Computational modeling, thermal Shift Assay, western blot and protein binding affinity experiments were later carried out for structural analyses. rRANKL mutants M200K/A/E showed a drastically reduced ability to induce osteoclast formation and did not demonstrate features of competitive inhibition against wild-type rRANKL. These mutants are all incapable of supporting osteoclastic polarization and bone resorption or activating RANKL-induced osteoclast marker gene transcription. Consistently, they were unable to induce calcium flux, and also showed a diminished induction of IκBa degradation and activation of NF-kB and NFATc1 transcriptional activity. Furthermore, the transcriptional activation of the antioxidant response element (ARE) crucial in modulating oxidative stress and providing cytoprotection was also unresponsive to stimulation with rM200s. Structural analyses showed that rM200 is located in a hydrophobic pocket critical for protein folding. Thermal shift and western blot assays suggested that rM200 mutants formed unstructured proteins, with disturbed trimerisation and the loss of affinity to its intrinsic receptors RANK and OPG. Taken together, we first demonstrates the underlying cause of M199-meidated ARO in a cellular and molecular level by establishing a phenotype in BMMs similar to observed in human samples. Further investigation hints the structural significance of a hydrophobic pocket within the TNF-like region. Combined with pharmaceutical studies on small-molecule drugs, this finding may represent a therapeutic target motif for future development of anti-resorptive treatments

Objectives. Secondary fracture healing is strongly influenced by the stiffness of the bone-fixator system. Biomechanical tests are extensively used to investigate stiffness and strength of fixation devices. The stiffness values reported in the literature for locked plating, however, vary by three orders of magnitude. The aim of this study was to examine the influence that the method of restraint and load application has on the stiffness produced, the strain distribution within the bone, and the stresses in the implant for locking plate constructs. Methods. Synthetic composite bones were used to evaluate experimentally the influence of four different methods of loading and restraining specimens, all used in recent previous studies. Two plate types and three screw arrangements were also evaluated for each loading scenario. Computational models were also developed and validated using the experimental tests. Results. The method of loading was found to affect the gap stiffness strongly (by up to six times) but also the magnitude of the plate stress and the location and magnitude of strains at the bone-screw interface. Conclusions. This study demonstrates that the method of loading is responsible for much of the difference in reported stiffness values in the literature. It also shows that previous contradictory findings, such as the influence of working length and very large differences in failure loads, can be readily explained by the choice of loading condition. Cite this article: A. MacLeod, A. H. R. W. Simpson, P. Pankaj. Experimental and numerical investigation into the influence of loading conditions in biomechanical testing of locking plate fracture fixation devices. Bone Joint Res 2018;7:111–120. DOI: 10.1302/2046-3758.71.BJR-2017-0074.R2

Background. Comminuted fractures involving the tibia are associated with a high level of complications including delayed healing and non-union, in conjunction with dramatically increased healthcare costs. Certain clinicians utilise a Pixel Value Ratio (PRV) of 1 to indicate such fracture healing. The subjectivity of this method has led to mixed outcomes including regenerate fracture. The poor prognosis of complex load bearing fractures is accentuated by the fact that no quantitative gold standard currently exists to which clinicians can reference regarding the definition of a healed fracture. The aim of the current study is to use patient specific finite element analysis of complex tibial fractures treated with Ilizarov frames to demonstrate callus maturation and to determine the optimum frame removal time. Methods. 3 patients (2 males, 1 female) were analysed following presentation with complex tibial fractures treated with Ilizarov frames. Patient specific computational analysis was performed according to radiographic data, incorporating maturing material properties to analyse the callus response to weight bearing over the healing timeframe. Computational results were compared to the PVR method to evaluate its efficacy in determining the optimum Ilizarov frame removal time. Results. All fractures were observed to clinically heal at a mean of 25.4 (±2.404) weeks. Following computational analysis however, the mean optimum Ilizarov frame removal time was seen to be 23.5 (±2.323) weeks. When compared with the PVR method, the suggested removal time presented a mean PVR of 1.025 (±0.017). Conclusion. Computational models of patient specific tibial fractures has shown promising correlations with the PVR method and has shown efficacy in predicting callus strength and subsequent optimum frame removal time. Level of Evidence. Level 4

Introduction. Whether anterior referencing (AR) or posterior referencing (PR) are optimal to position and size the femoral component in Total Knee Arthroplasty (TKA) remains controversial. This controversy stems, in part, from a lack of understanding of whether one technique more consistently balances the medial/lateral collateral ligaments (MCL & LCL) in flexion and extension. Therefore, our goal was to compare AR and PR in terms of: (1) maximum MCL and LCL forces in passive flexion, and (2) medial and lateral gaps at full extension and 90‖ of flexion. In addition, we identified geometric landmarks that could help predict the ligament forces during flexion. Methods. Computational models of six knees were virtually implanted with TKAs based on our previously-developed framework. AR and PR were simulated in each of the six models. A Posterior Stabilized implant was utilized. Standard AR and PR cuts and component positioning were simulated with the femoral component aligned parallel to the transepicondylar axis. In both AR and PR models, the distal femoral cut and the proximal tibial cut were perpendicular to the femoral and tibial mechanical axis, respectively. The amount of posterior bone resected with AR knees ranged from 4.2 to 10.8 mm, and with PR knees ranged from 4.2 to 8 mm. Ligament properties were standardized to reflect a balanced knee at full extension. Passive flexion under 500 N of compression was applied and the MCL and LCL forces were predicted. A new measure, the MCL ratio, that incorporated the femoral insertion of the anterior fiber of MCL relative to the posterior and distal femoral cuts was estimated (Fig. 1). A varus/valgus moment of 6 Nm was applied at full extension and 90‖ of flexion, and the corresponding lateral and medial gaps were measured. Results. In passive flexion, the maximum MCL force ranged from 2 to 87 N in AR and from 17 to 127 N in PR (Fig. 2). The LCL forces decreased to zero before 25‖ of flexion in all knees. The MCL ratio corresponded to the MCL force; the larger the MCL ratio, the larger was the maximum MCL force (Fig. 2). At full extension, AR and PR knees were balanced with a maximum difference in medial-lateral gap < 1 mm. However, in flexion, only two out of the six AR and PR knees produced a difference in medial-lateral gap < 2 mm. Conclusion. Neither AR nor PR consistently produced higher or lower maximum MCL force in flexion despite being well balanced in extension and aligned with the transepicondylar axis. Rather, the more the posterior bone resection, independent of AR or PR, the less was the maximum MCL force in flexion. Knees that produced symmetrical gaps at full extension and 90° of flexion were the ones with the lowest maximum MCL force. Therefore, our findings suggest that less MCL force in flexion promotes a more balanced knee. The MCL ratio corresponded to the variations in maximum MCL force in flexion; it may be used to help produce a more well-balanced knee joint. For any figures or tables, please contact authors directly

There are a number of periprosthetic femoral fracture (PFF) fixation failures. In several cases the effect of fracture configuration on the performance of the chosen fixation method has been underestimated. As a result, fracture movement within the window that seems to promote callus formation has not been achieved and fixations ultimately failed. This study tested the hypothesis that: PFF configuration and the choice of plate fixation method can be detrimental to healing. A series of computational models were developed, corroborated against measurements from a series of instrumented laboratory models and in vivo case studies. The models were used to investigate the fixation of different fracture configurations and plate fixation parameters. Surface strain and fracture movement were compared between the constructs. A strong correlation between the computational and experimental models was found. Computational models showed that unstable fracture configurations increase the stress on the plate fixation. It was found that bridging length plays a pivotal role in the fracture movement. Rigid fixations, where there is clinical evidence of failure, showed low fracture movement in the models (<0.05mm); this could be increased with different screw and plate configurations to promote healing. In summary our results highlighted the role of fracture configuration in PFF fixations and showed that rigid fixations that suppress fracture movement could be detrimental to healing

Background. Wear and fatigue damage to polyethylene components remain major factors leading to complications after total knee and unicompartmental arthroplasty. A number of wear simulations have been reported using mechanical test equipment as well as computer models. Computational models of knee wear have generally not replicated experimental wear under diverse conditions. This is partly because of the complexity of quantifying the effect of cross-shear at the articular interface and partly because the results of pin-on-disk experiments cannot be extrapolated to total knee arthroplasty wear. Our premise is that diverse experimental knee wear simulation studies are needed to generate validated computational models. We combined five experimental wear simulation studies to develop and validate a finite-element model that accurately predicted polyethylene wear in high and low crosslinked polyethylene, mobile and fixed bearing, and unicompartmental (UKA) and tricompartmental knee arthroplasty (TKA). Methods. Low crosslinked polyethylene (PE). A finite element analysis (FEA) of two different experimental wear simulations involving TKA components of low crosslinked polyethylene inserts, with two different loading patterns and knee kinematics conducted in an AMTI knee wear simulator: a low intensity and a high intensity. Wear coefficients incorporating contact pressure, sliding distance, and cross-shear were generated by inverse FEA using the experimentally measured volume of wear loss as the target outcome measure. The FE models and wear coefficients were validated by predicting wear in a mobile bearing UKA design. Highly crosslinked polyethylene (XLPE). Two FEA models were constructed involving TKA and UKA XLPE inserts with different loading patterns and knee kinematics conducted in an AMTI knee wear simulator. Wear coefficients were generated by inverse FEA. Results. Predicted wear rates were within 5% of experimental wear rates during validation tests. Unicompartmental mobile bearing back-side wear accounted for 46% of the total wear in the mobile bearing. Wear during the swing phase was 38% to 44% of total wear. Discussion & Conclusions. Crosslinking polyethylene primarily decreased (by nearly 10-fold) the wear generated by cross-shear. This result can be explained by the reduced propensity of crosslinked polyethylene molecules to orient in the dominant direction of sliding. A highly crosslinked fixed-bearing polyethylene insert can provide high wear performance without the increased risk for mobile bearing dislocation. Finite element analysis can be a robust and efficient method for predicting experimental wear. The value of this model is in rapidly conducting screening studies for design development, assessing the effect of varying patient activity, and assessing newer biomaterials. This FEA model was experimentally validated but requires clinical validation

Introduction. Joint mechanics and implant performance have been shown to be sensitive to ligament properties [1]. Computational models have helped establish this understanding, where optimization is typically used to estimate ligament properties for recreation of physically measured specimen-specific kinematics [2]. If available, contact metrics from physical tests could be used to improve the robustness and validity of these predictions. Understanding specimen-specific relationships between joint kinematics, contact metrics, and ligament properties could further highlight factors affecting implant survivorship and patient satisfaction. Instrumented knee implants offer a means to measure joint contact data both in-vivo and intra-operatively, and can also be used in a controlled experimental environment. This study extends on previous work presented at ISTA [3], and the purpose here was to evaluate the use of instrumented implant contact metrics during optimization of ligament properties for two specimens. The overarching goal of this work is to inform clinical joint balancing techniques and identify factors that are critical to implant performance. Methods. Total knee arthroplasties were performed on 4 (two specimens modeled) cadeveric specimens by an experienced orthopaedic surgeon. An instrumented trial implant (VERASENSE, OrthoSensor, Inc., Dania Beach, FL) was used in place of a standard insert. Experimentation was performed using a simVITROTM controlled robotic musculoskeletal simulator (Cleveland Clinic, Cleveland, OH) to apply intra-operative style loading and measure tibiofemoral kinematics. Three successive laxity style tests were performed at 10° knee flexion: anterior-posterior force (±100 N), varus-valgus moment (±5 Nm), and internal-external moment (±3 Nm). Tibiofemoral kinematics and instrumented implant contact metrics were measured throughout testing (Fig. 1). Specimen-specific finite element models were developed for two of the tested specimens and solved using Abaqus/Explicit (Dassault Systèmes). Relevant ligaments and rigid bone geometries were defined using specimen-specific MRIs. Virtual implantation was achieved using registration and each ligament was modeled as a set of nonlinear elastic springs (Fig. 1). Stiffness values were adopted from the literature [2] while the ligament slack lengths served as control variables during optimization. The objective was to minimize the root mean square difference between VERASENSE measured tibiofemoral contact metrics and the corresponding model results (Fig. 1). Results and Discussion. The models for both specimens successfully recreated joint kinematics with average errors less than 4° in rotations, and 3 mm in translations (not shown). Minus a systematic offset in θ for specimen 3, AFD and θ contact kinematics also realized good agreement for both specimens (Fig. 2). Contact forces were generally over-predicted, though both specimens recreated the experimental trends (Fig. 2). The present work shows continued progress towards simulation based tools that can be used for both research and to support the clinical decision making process. A separate ISTA submission presents assessment of these model's predictive capacity, while future work will evaluate additional specimens, and explore the sensitivity to uncertainties in experimental and modeling parameters. Acknowledgements. This work was supported by Orthosensor Inc. For any figures or tables, please contact authors directly (see Info & Metrics tab above).

Computational modelling has the potential of becoming a useful tool in assessing revision risk on a patient-specific basis. However, there are many difficulties encountered in generating subject-specific computational models that have unknown influences on such predictions, e.g. accuracy of the anatomical geometry and material properties of the patient. This study compares the influence of these two patient-specific parameters on predictions of revision risk due to aseptic loosening. First, X-rays from seventeen patients were processed using previously developed technique utilising rigid scaling of a generic femur to match selected dimensions from each patient’s post-operative X-ray and, then, the same set of 3D models was obtained by using an automated technique that generates 3D extra-cortical geometries from planar X-rays using a combination of 2D contour extraction and 3D warping of a generic model to match the extracted contour. A cement and cement-metal interfacial damage accumulation algorithm developed previously was used. For each geometric set two types of simulations were performed. First, constant cortical and cancellous bone apparent Young’s moduli were assumed. A second set of simulations used age-dependent Young’s moduli for each bone type. Walking and stair-climbing activities were simulated. Resultant migration of the prostheses was used to indicate revision risk. Factorial analysis has shown that the geometry has a larger influence on resultant migration magnitude for each case; however, unexpectedly, using more realistic geometry weakened the strength of predictions. This is most likely to be due ongoing mesh-induced contact problems