Summary Statement. The constraint behavior of total knee arthroplasty (TKA) prosthesis usually has to be physically tested. This study presents a computer simulation model using finite element analysis (FEA) and demonstrates its effectiveness in predicting the femorotibial constraint behavior of TKA implants. Introduction. TKA prostheses are semi-constrained artificial joints. A well-functioning TKA prosthesis should be designed with a good balance between stability and mobility, meaning the femorotibial constraint of the artificial joint cannot be excessive or too lax. To assess the constraint behavior of a TKA prosthesis, physical testing is usually required, and an industrial test standard has been developed for this purpose. Benefiting from technological advancement, computer simulation has become increasingly useful in many industries, including medical device research and development. FEA has been extensively used in stress analysis and structural evaluation of various orthopaedic implants. This study presented an FEA-based simulation to evaluate the femorotibial constraint behavior of TKA prosthesis, and demonstrated the effectiveness of the method by validating it through physical testing. Methods. A Cruciate Retaining (CR) TKA prosthesis design (Optetrak Logic CR, size 3, Exactech, FL, USA) was used in this study. The prosthesis system consists of a femoral component, a tibial insert, and a tibial baseplate. CAD models of the implants assembled at 0° of flexion were used for the simulation. Finite element models were generated using 10-node tetrahedral elements, with all materials considered linear elastic. Boundary conditions were set up according to the ASTM F1223 standard. The tibial baseplate was fixed distally. A constant compressive force (710 N) was applied on the femoral component. Nonlinear Surface-Surface-Contact was defined at the femorotibial articulating surfaces as well as between the tibial insert and tibial baseplate. A coefficient of friction of 0.2 determined from the physical test was input into the simulation. The femoral component was driven under a displacement-controlled scheme to slide along the anterior-posterior (AP) direction on the tibial insert. At each time step, constraint force occurring at the articulating surface was derived from the reaction force at the distal fixation of the tibial baseplate. The force-displacement curve was plotted by combining the results of all time steps to characterize the constraint behavior of the prosthesis. A nonlinear FEA solver (NX Nastran SOL601, Siemens, TX, USA) was used to solve the simulation. In addition, five samples of the prostheses were physically tested per ASTM F1223. Simulation results were compared to the physical testing. Results. The simulation successfully captured the movement of contact location and pressure along the movement of the femoral component. The force-displacement curve predicted by the simulation exhibited a very close hysteresis loop profile as the results of physical testing. Using the curve slope from 0 to 5 mm to characterise the constraint in the most relevant displacement range, the simulation predicted 45.7 N/mm anteriorly and 36.4 N/mm posteriorly, which are less than 10% different from the physical testing results (46.4 N/mm anteriorly and 39.6 N/mm posteriorly). Discussion/Conclusion. This study demonstrated that the simulation was able to closely predict the femorotibial constraint behavior of the TKA prosthesis under ASTM F1223 testing. The simulation results resembled the physical test results not only in the general profile of the curve but also in the magnitude of slope values. The increased difference at the far anterior region could be related to the fact that no material nonlinearity was considered in the current simulation, a factor that could be improved in future studies. A validated simulation method could be very useful in TKA prosthesis design. Since no physical prototypes are required, design evaluation and optimization can be achieved in a much easier and faster manner

Surgical education of fracture fixation biomechanics relies mainly on simplified illustrations to distill the essence of the underlying principles. These mostly consist of textbook drawings or hands-on exercises during courses, both with unique advantages such as broad availability and haptics, respectively. Computer simulations are suited to bridge these two approaches; however, the validity of such simulations must be guaranteed to teach the correct aspects. Therefore, the aim of this study was to validate finite element (FE) simulations of bone-plate constructs to be used in surgical education in terms of fracture gap movement and implant surface strain. The validation procedure was conducted in a systematic and hierarchical manner with increasing complexity. First, the material properties of the isolated implant components were determined via four-point bending of the plate and three-point bending of the screw. Second, stiffness of the screw-plate interface was evaluated by means of cantilever bending to determine the properties of the locking mechanism. Third, implant surface strain and fracture gap motion were measured by testing various configurations of entire fixation constructs on artificial bone (Canevasit) in axial compression. The determined properties of the materials and interfaces assessed in these experiments were then implemented into FE models of entire fixation constructs with different fracture width and screw configurations. The FE-predicted implant surface strains and fracture gap motions were compared with the experimental results. The simulated results of the different construct configurations correlated strongly with the experimentally measured fracture gap motions (R. 2. >0.99) and plate surface strains (R. 2. >0.95). In a systematic approach, FE model validation was achieved successfully in terms of fracture gap motion and implant deformation, confirming trustworthiness for surgical education. These validated models are used in a novel online education tool OSapp (. https://osapp.ch/. ) to illustrate and explain the biomechanical principles of fracture fixations in an interactive manner

This study aims to create a novel computational workflow for frontal plane laxity evaluation which combines a rigid body knee joint model with a non-linear implicit finite-element model wherein collateral ligaments are anisotropically modelled using subject-specific, experimentally calibrated Holzpfel-Gasser-Ogden (HGO) models. The framework was developed based on CT and MRI data of three cadaveric post-TKA knees. Bones were segmented from CT-scans and modelled as rigid bodies in a multibody dynamics simulation software (MSC Adams/view, MSC Software, USA). Medial collateral and lateral collateral ligaments were segmented based on MRI-scans and are modelled as finite elements using the HGO model in Abaqus (Simulia, USA). All specimens were submitted varus/valgus loading (0-10Nm) while being rigidly fixed on a testing bench to prevent knee flexion. In subsequent computer simulations of the experimental testing, rigid bodies kinematics and the associated soft-tissue force response were computed at each time step. Ligament properties were optimised using a gradient descent approach by minimising the error between the experimental and simulation-based kinematic response to the applied varus/valgus loads. For comparison, a second model was defined wherein collateral ligaments were modelled as nonlinear no-compression spring elements using the Blankevoort formulation. Models with subject-specific, experimentally calibrated HGO representations of the collateral ligaments demonstrated smaller root mean square errors in terms of kinematics (0.7900° +/− 0.4081°) than models integrating a Blankevoort representation (1.4704° +/− 0.8007°). A novel computational workflow integrating subject-specific, experimentally calibrated HGO predicted post-TKA frontal-plane knee joint laxity with clinically applicable accuracy. Generally, errors in terms of tibial rotation were higher and might be further reduced by increasing the interaction nodes between the rigid body model and the finite element software. Future work should investigate the accuracy of resulting models for simulating unseen activities of daily living

Abstract. Objectives. Catastrophic neck injuries in rugby tackling are rare (2 per 100,000 players per year) with 38% of these injuries occurring in the tackle. The aim of this study was to determine the primary mechanism of cervical spine injury during rugby tackling and to highlight the effect of tackling technique on intervertebral joint loads. Methods. In vivo and in vitro experimental data were integrated to generate realistic computer simulations representative of misdirected tackles. MRI images were used to inform the creation of a musculoskeletal model. In vivo kinematics and neck muscle excitations were collected during lab-based staged tackling of the player. Impact forces were collected in vitro using an instrumented anthropometric test device during experimental simulations of rugby collisions. Experimental kinematics and muscle excitations were prescribed to the model and impact forces applied to seven skull locations (three cranial and four lateral). To examine the effects of technique on intervertebral joint loads the model's neck angle was altered in steps of 5° about each rotational axis resulting in a total of 1,623 experimentally informed simulations of misdirected tackles. Results. Neck flexion angles and cranial impact locations had the largest effects on maximal compression, anterior shear and flexion moment loads. During posterior cranial impacts compression forces and flexion moments increased from 1500 to 3200 N and 30 to 60 Nm respectively between neck angles of 30° extension and 30° flexion. This was more evident at the C5-C6 and C6-C7 joints. Anterior shear loads remained stable throughout neck angle ranges however during anterior impacts they were directed posteriorly when the neck was flexed. Conclusions. The combination of estimated joint loads in the lower cervical spine support buckling as the primary injury mechanism of anterior bilateral facet dislocations observed in misdirected rugby tackles and highlights the importance of adopting a correct tackling technique. 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

Summary Statement. Femorotibial constraint is a key property of a total knee arthroplasty (TKA) prosthesis and should reflect the intended function of the device. With a validated simulation methodology, this study evaluated the constraint of two TKA prostheses designed for different intentions. Introduction. TKA prostheses are semi-constrained artificial joints. Femorotibial constraint level is a major property of a prosthesis and should be designed to match the device's intended function. Cruciate Retaining (CR) prostheses are usually indicated for patients with a functioning posterior cruciate ligament (PCL). For patients without a fully functioning PCL, CR-Constrained (CRC) prostheses with additional built-in constraint may be indicated. A CRC prosthesis usually consists of a CR femoral component and a tibial insert which has a more conforming sagittal profile to offer an increased femorotibial constraint. This study evaluated the anterior-posterior (AP) constraint behavior of two lines of prostheses (CR and CRC) from a same TKA product family. Using a validated computer simulation approach, multiple sizes of each product line were evaluated. Methods. Both the CR and CRC prostheses are from the same TKA product family (Optetrak Logic, Exactech, FL, USA) and share identical femoral components and tibial baseplates. The CRC tibial inserts have a more conforming sagittal profile than the CR tibial inserts, especially in the anterior aspect. Three sizes (sizes 1, 3, and 5) from each product line were included in this study. Computer simulations using finite element analysis (FEA) were performed to evaluate the femorotibial constraint of each prosthesis per ASTM F1223 standard [1]. The simulation has been validated by comparison with physical testing (more details submitted in a separate paper to CORS 2013). Briefly, FEA models were created using 10-node tetrahedral elements with all materials considered linear elastic. The tibial baseplate was distally fixed and a constant compressive force (710 N) was applied to the femoral component. Nonlinear Surface-Surface-Contact was established at the articulating surfaces, as well as between the tibial insert and the tibial baseplate. A coefficient of friction of 0.1 was assumed for all articulations [2]. The femoral component was driven under a displacement-controlled scheme to slide along AP direction on the tibial insert. Constraint force occurring at the articulation was derived from the reaction force at the distal fixation; thus, the force-displacement curve can be plotted to characterise the constraint behavior of the prosthesis. A nonlinear FEA solver (NX Nastran SOL601, Siemens, TX, USA) was used to solve the simulations. Results. The force-displacement curves predicted by the simulation exhibited the hysteresis loop appearance for both CR and CRC prostheses. The profile of the curves was generally consistent across different sizes for both product lines. The anterior constraint of the CRC prosthesis was significantly greater than the CR prosthesis. The posterior constraint of the CRC prosthesis was also slightly greater. Larger sizes exhibited reduced constraint compared to smaller sizes. Discussion/Conclusion. The increased constraint of the CRC prosthesis revealed in the study is consistent with the geometrical characteristics and the functional intent of the device. The CRC tibial insert is expected to provide significantly greater anterior constraint than the CR prosthesis to prevent paradoxical femoral translation when the patient's PCL is not fully functioning. The CRC tibial insert is also expected to provide slightly increased posterior constraint due to its elevated posterior lip. The observed hysteresis loop appearance is consistent with physical testing and the existence of friction. The reduced constraint on larger sizes is functionally desirable to offer proportional translation freedom. This study demonstrated the effectiveness of the simulation approach in quantifying the constraint behavior of different TKA prosthesis designs

Osteochondral (OC) grafting is one available method currently used to repair full thickness cartilage lesions with good results clinically when grafting occurs in patients with specific positive prognostic factors. However, there is poor understanding of the effect of individual patient and surgical factors. With limited tissue availability, development of Finite Element (FE) models taking into account these variations is essential. The aim of this study was to evaluate the effect of altering the material properties of OC grafts and their host environment through computer simulation. A generic FE model (ABAQUS CAE 2017) of a push-out test was developed as a press-fit bone cylinder (graft) sliding inside a bone ring (host tissue). Press-fit fixation was simulated using an interference fit. Overlap between host and graft (0.01mm–0.05mm) and coefficient of friction (0.3–0.7) were varied sequentially. Bone Young's moduli (YM) were varied individually between graft and host within the range of otherwise derived tissue moduli (46MPa, 82MPa, 123MPa). Increasing both overlap and frictional coefficient increased peak dislodging force independently (overlap: 490% & frictional coefficient: 176% across range tested). Increasing bone modulus also increased dislodging force, with host bone modulus (107%, 128%, and 140% increase across range, when Graft YM = 123MPa, 82 MPa, and 46MPa, respectively) having a greater influence than graft modulus (28%, 19% and 10% increase across range, when Host YM = 123 MPa, 82MPa and 46MPa, respectively). As anticipated increasing overlap and friction caused an increase in force necessary to dislodge the graft. Importantly, differentially changing the graft and host material properties changed the dislodging force indicating that difference between graft and host may be an important factor in the success or failure clinically of osteochondral grafting

The high risk and the associated high mortality of secondary, contralateral hip fractures [1,2] could justify internal, invasive prophylactic reinforcement of the osteoporotic proximal femur to avoid these injuries in case of a low energy fall. Previous studies have demonstrated high potential of augmentation approaches [3,4,5], but to date there has no ideal solution been found. The development of optimized reinforcement strategies can be aided with validated computer simulation tools that can be used to evaluate new ideas. A validated non-linear finite element (FE) simulation tool was used here to predict the yield and fracture load of twelve osteoporotic or osteopenic proximal femora in sideways fall based on high resolution CT images. Various augmentation strategies using bone cement or novel metal implants were developed, optimized and virtually performed on the bone models. The relative strengthening compared to the non-augmented state was evaluated using case-specific FE analyses. Strengthening effect of the cement-based augmentation was linearly proportional to cement volume and was significantly affected by cement location. With the clinically acceptable 12.6 ± 1.2 ml volume and optimized location of the cement cloud, compared to the non-augmented state, 71 ± 26% (42 – 134%) and 217 ± 166% (83 – 509%) increase in yield force and energy was reached, respectively. These were significantly higher than previously published experimental results using the “central” cement location [5], which could be well predicted by our FE models. The optimized metal implant could provide even higher strengthening effect: 140 ± 39% (76 – 194%) increase in yield force and +357 ± 177% (132 – 691%) increase in yield energy. However, for metal implants, a higher risk of subcapital fractures was indicated. For both cement and metal, the originally weaker bones were strengthened exponentially more compared to the stronger ones. The ideal solution for prophylactic augmentation should provide an appropriate balance between the requirements of being clinically feasible, ethically acceptable and mechanically sufficient. Even with the optimized location, the cement-based approach may not provide enough strengthening effect and adequate reproducibility of the identified optimal cement cloud position may not be achieved clinically. While the metal implant based strategy appears to be able to deliver the required strengthening effect, the ethical acceptance of this more invasive option is questionable. Further development is therefore required to identify the ideal, clinically relevant augmentation strategy. This may involve new cement materials, less invasive metal implants, or a combination of both. The FE simulation approach presented here could help to screen the potential ideas and highlight promising candidates for experimental evaluation

Total ankle replacement (TAR) is the main surgical option in case of severe joint osteoarthritis. The high failure rate of current TAR is often associated to inappropriate prosthetic articulating surfaces designed according to old biomechanical concepts such the fixed axis of rotation, thus resulting in non-physiological joint motion. A recent image-based 3D morphological study of the normal ankle (Siegler et al. 2014) has demonstrated that the ankle joint surfaces can be approximated by a saddle-shaped cone with its apex located laterally (SSCL). We aimed at comparing the kinematic effects of this original solution both with the intact joint and with the traditional prosthetic articulating surfaces via in-silico models and in-vitro measurements. Native 3D morphology of ten normal cadaver ankle specimens was reconstructed via MRI and CT images. Three custom-fit ankle joint models were then developed, according to the most common TAR designs: cylindrical, symmetrically-truncated medial apex cone (as in Inman's pioneering measures), and the novel lateral apex cone, i.e. SSCL. Bone-to-bone motion, surface-to-surface distance maps, and ligament forces and deformations were evaluated via computer simulation. Prototypes of corresponding prosthesis components were designed and manufactured via 3D-printing, both in polymer-like-carbon and in cobalt-chromium-molybdenum powders, for in-vitro tests on the cadaver specimens. A custom testing rig was used for application of external moments to the ankle joint in the three anatomical planes; a motion tracking system with trackers pinned into the bone was used to measure tibial, talar and calcaneal motion (Franci et al. 2009), represented then as tibiotalar, subtalar and ankle complex 3D joint rotations. Each ankle specimen was tested in the intact joint configuration and after replacement of the articulating surfaces according with the three joint models: cylindrical, medial apex cone and SSCL. Results. Small intra-specimen data variability in cycle-to-cycle joint kinematics was found in all cadaver ankles, the maximum standard deviation of all rotation patterns being smaller than 2.0 deg. In-silico ligament strain/stress analysis and in-vitro joint kinematic and load transfer measurements revealed that the novel SSCL surfaces reproduce more natural joint patterns than those with the most common surfaces used in current TAR. TAR based on a saddle-shaped skewed truncated cone with lateral apex is expected to restore more normal joint function. Additional tests are undergoing for further biomechanical validation. The present study has also demonstrated the feasibility and the quality of the full process of custom TAR design and production for any specific subject. This implies a thorough procedure, from medical imaging to the production of artificial surfaces via 3D printing, which is allowing for personalised implants to become the future standard in total joint replacement

Dislocation remains a major concern after total hip replacement, and is often attributed to malposition of the components. The optimum position for placement of the components remains uncertain. We have attempted to identify a relatively safe zone in which movement of the hip will occur without impingement, even if one component is positioned incorrectly. A three-dimensional computer model was designed to simulate impingement and used to examine 125 combinations of positioning of the components in order to allow maximum movement without impingement. Increase in acetabular and/or femoral anteversion allowed greater internal rotation before impingement occurred, but decreases the amount of external rotation. A decrease in abduction of the acetabular components increased internal rotation while decreasing external rotation. Although some correction for malposition was allowable on the opposite side of the joint, extreme degrees could not be corrected because of bony impingement.

We introduce the concept of combined component position, in which anteversion and abduction of the acetabular component, along with femoral anteversion, are all defined as critical elements for stability.