Towards Functional Preoperative Planning in Orthopaedic Surgery

Las cirugíıas del aparato locomotor suponen más de 20 millones de intervencionesanuales para la corrección de lesiones que afectan a músculos, articulaciones,ligamentos, tendones, huesos o nervios; elementos que conforman el sistema musculoesquelético. Este tipo de afecciones de la biomecánica pueden tener diversos orígenes; siendo los principales los traumatismos, las lesiones degenerativas en huesos y tejidos blandos, los malos hábitos posturales o motores, y los de origen congénito.El uso de las tecnologías actuales en los procesos de corrección de estas afecciones forma parte del día a día en los quirófanos y en la monitorización de los pacientes.Sin embargo, el uso de técnicas computacionales que permitan la preparación de las intervenciones quirúrgicas antes de proceder con la cirugía están todavía lejos de formar parte del proceso de evaluación preoperatoria en este tipo de lesiones. Por este motivo, el objetivo principal de esta tesis consiste en demostrar la viabilidad del uso de herramientas computacionales en la planificación preoperatoria de diferentes cirugías ortopédicas.Entre los tipos de cirugías más comunes, la mayor parte de ellas se centran en las articulaciones del tren inferior de la anatomía humana. Por este motivo, este trabajo se centraría en el análisis de diferentes cirugías cuya finalidad es solucionar lesiones en las principales articulaciones del tren inferior: región sacrolumbar, cadera, rodilla y tobillo.Para poder realizar el análisis de estas cirugías se hizo uso de algunas de lasherramientas computacionales más usadas habitualmente y cuya capacidad en diversos ámbitos ha sido comprobada. Se ha utilizado la reconstrucción 3D para la obtención de modelos anatómicos sobre los que comprobar la viabilidad de las cirugías. Estas reconstrucciones se basan en las imágenes médicas obtenidas mediante Tomografia Axial Computerizada (TAC) o Resonancia Magnética (RM). Las imágenes procedentes de RM permiten diferenciar todos los tejidos de la anatomía, incluyendo los blandos tales como tendones o cartílagos; mientras que los TAC facilitan la diferenciación de los huesos. Esta última es la prueba más habitual en los diagnósticos.Para su análisis y reconstrucción se hizo uso de los software Mimics v 20.0 y3-matic 11.0 (Materialise NV, Leuven, Belgium). Como alternativa para la generación de los modelos cuando no se dispone de las imágenes necesarias para realizar la reconstrucción o cuando se requiere dotar de flexibilidad a estos modelos, se recurrió al modelado en el software de análisis por elementos finitos Abaqus/CAE v.6.14 (Dassault Syst`emes, Suresnes, France). Dicho software fue además utilizado para la simulación del efecto de las diferentes cirugías sobre la región de interés. Para resalizar las simulaciones, se incluyeron en los modelos aquellos parámetros, elementos y condiciones necesarios para poder representar las caraterísticas propias de cada cirugía. Finalmente, para aquellas situaciones que requerían del análisis de datos se hizo uso de tecnologías de machine learning. La solución seleccionada para estos casos fueron las redes neuronales artificiales (ANN). Dichas redes se desarrollaronhaciendo uso del software MATLAB R2018b (MathWorks, Massachusetts, USA).El estudio de la rodilla se centra en uno de los ligamentos clave en la estabilidad de la rótula y que, sin embargo, es uno de los menos analizados hasta ahora, el ligamento medial patelofemoral. La reconstrucción de este ligamento es la principal solución clínica para solventar esta inestabilidad y diferentes cirugías utilizadas para dicho fin han sido analizadas mediante el desarrollo de un modelo paramétrico en elementos finitos que permita su simulación. En este modelo es posible adaptar la geometría de la rodilla de forma que se puedan simular diferentes condiciones que pueden afectar a la estabilidad de la rótula, tales como la displasia troclear y la patella alta.El estudio de la región sacrolumbar se centra en el análisis de diferentes posibles configuraciones para las cirugías de fusión vertebral. El análisis se centró en la fijación con tornillos y la influencia del Polimetimetacrilato (PMMA) como elemento de fijación en las vértebras. Para ello, se reconstruyó el modelo óseo de diferentes pacientes que necesitaron este tipo de intervención. Sobre estos modelos se simularon mediante elementos finitos las diferentes configuraciones consideradas de forma que se pudiera comparar su comportamiento en diferentes casos.En el caso de la cadera, el estudio se centra en el análisis de la artroplastia total de cadera, que implica el reemplazo de la articulación anatómica por una prótesis habitualmente de titanio. Cuando este tipo de cirugías es realizado, es común que surjan posteriormente problemas derivados de la disposición de la prótesis y que pueden llevar al pinzamiento entre sus componentes y, en algunas ocasiones, su dislocación.Esto ocurre cuando el rango de movimiento de la articulación es reducido. Este tipo de sucesos son más comunes cuando se realizan los movimientos de extensión externa (EE) o de rotación interna (RI) de la extremidad. El estudio se desarrolló con el objetivo de elaborar una herramienta computacional capaz de predecir este choque y dislocación basándose en el diámetro de la cabeza del femur y de los ángulos de abducción y anteversión. Para ello, se recurrió al uso de redes neuronales artificales(ANN). Se configuró una red independiente para cada movimiento (EE y RI) y cada posible evento (pinzamiento y dislocación), de forma que se obtuvieron cuatro redes completamente independientes. Para el entrenamiento y primer testeo de las redes se recurrió a un modelo paramétrico en elementos finitos de la prótesis con el que se realizaron diferentes simulaciones determinando el rango de movimiento para cada caso. Finalmente, las redes fueron de nuevo validadas con el uso de datos procedentes de pacientes que sufrieron dislocación tras ser sometidos a este tipo de cirugías.Finalmente, el estudio de la región del tobillo se centró en la lesión de la sindesmosis del tobillo. Este tipo de lesiones implica la rotura de algunos de los ligamentos que unen los principales huesos de esta articulación (tibia, peroné y astrágalo) junto con parte de la membrana intraósea, que se extiende a lo largo de la tibia y el peroné ligando ambos huesos. Cuando se produce este tipo de lesiones, es necesario recurrir a la inclusión de elementos que fijen la articulación y prevengan la separación de los huesos. Los métodos más comunes y que centran este análisis comprenden la fijación con tornillos y la fijación mediante botón de sutura. Para poder realizar un análisis que permita comparar la efectividad y incidencia de este tipo de cirugías se recurrióa la reconstruccción 3D de la articulación de un paciente que sufrió este tipo de lesión. Con este modelo geométrico, se procedió al desarrollo de diferentes modelos en elementos finitos que incluyeran cada una de las alternativas consideradas. Las simulaciones de estos modelos junto a las situaciones anatómicas y lesionadas, permitió hacer una aproximación sobre la solución quirúrgica que mejor restablece el estado incial sano de la región afectada.Locomotor system surgeries represents more the 20 million interventions per year for the correction of injuries that affect muscles, joints, ligaments, tendons, bones or nerves; elements that form themusculoskeletal system. This kind of biomechanical affections may have several sources, being the main ones traumas, bones and soft tissues degenerative injuries, poor postural or motor habits and those of congenital source. The use of current technologies in the correction process for these injuries is part of the day-to-day in the operating rooms and the monitoring of patients. However, the use of computational tools that allow preoperative planning is still far from being part of the preoperative evaluation process in this kind of injuries. For this reason, the main goal of this thesis consists in demonstrating the viability of the use of computational tools in the preoperative planning of different orthopaedic surgeries. Among the most common surgeries, most of them focus in the lower body joints of the human anatomy. For this reason, this work will focus in the analysis of different surgeries whose purpose is to solve injuries in the main joints of the lower body: lumbosacral region, hip, knee and ankle. Some of the most commonly used computational tools, and whose capability in different fields has been widely proven, were used in order to be able of performing the analysis of these surgeries. 3D reconstruction has been used for obtaining anatomical models in which the viability of the surgeries could be verified. These reconstructions are based on the medical images obtained through Computerized Tomography (CT) or Magnetic Resonance Imaging (RMI). Images from RMI allow differentiating all the tissues of the anatomy, including soft ones such as tendons and cartilages; while CT scans make easier the bones differentiation. This last procedure is the most commonly used in diagnoses. For their analysis and reconstruction software Mimics v 20.0 and 3-Matic 11.0 (Materialise NV, Leuven, Belgium) were used. As alternative for the models generation when the necessary images for the reconstruction are not available or when flexibility is required for these models, modelling in the Finite Element Analysis software Abaqus/CAE v.6.14 (Dassault Syst‘emes, Suresnes, France) was used. This software was also used for the simulation of the effects of the different surgeries in the interest region. In order to perform the simulations, those parameters, elements and conditions necessary to represent the characteristics of each surgery were included. Finally, for those situations requiring data analysis, machine learning technologies were used. The selected solution for these cases were Artificial Neural Networks (ANN). These networks were developed using the software MATLAB R2018b (MathWorks, Massachusetts, USA). The study of the knee joint focuses in one of the key ligaments for the patellar stability and which, however, is one of the least analysed so far, the medial patellofemoral ligament. The reconstruction of this ligament is the main clinical solution for solving this instability and different surgeries used for that purpose have been analysed through the development of a finite element parametric model that allows their simulation. In this model adapting knee geometry is possible so that those conditions that can affect the stability of the patella, such as trochlear dysplasia or patella alta, can be simulated. The study of the lumbosacral region focuses in the analysis of different possible configurations for spine fusion surgeries. The analyses focused in the pedicle screws fixation and the influence of polymethyl methacrylate (PMMA) as fixation element in the vertebrae. To do this, osseous models for different patients that required this kind of intervention were reconstructed. The different configurations considered were simulated on these models through finite element analysis comparing their behaviour. In the case of the hip, the study focuses in the analysis of the total hip arthroplasty, which implies replacing the anatomical joint by a prosthesis, usually made of titanium. When this kind of surgery is performed, it is common for later issues arising from the arrangement of the prosthesis and which can lead to impingement between its components and, on some occasions, their dislocation. This happens when the range of movement of the joint is limited. This kind of events are more common when the external extension (EE) or internal rotation (IR) movements of the leg are performed. The study was developed with the goal of elaborating a computational tool able to predict the impingement and dislocation based on the diameter of the head of the femur and the anteversion and abduction angles. To do this, artificial neural networks (ANN) were used. An independent network was configured for each movement (EE and IR) and for each possible event (impingement and dislocation), so that four completely independent networks. For the training and the first testing of the networks, a parametric finite element model of the hip was used; with which different simulations were performed determining the range of movement for each case. Finally, the networks were validated again with the use of data proceeding from patients that suffered dislocation after going through this kind of surgery. Finally, the study of the ankle region focused in the ankle syndesmosis injury. This kind of injuries implies the tear of some ligaments that connect the main bones of the joint (tibia, fibula and talus) together with part of the intraosseous membrane, which extends along the tibia and fibula linking both bones. When this kind of injuries happens, it is necessary to resort to the inclusion of elements that fix the joint and prevent the bones distance. The most common methods, which focus this analysis, include the screws fixation and the suture button fixation. In order to carry out an analysis that allows comparing the effectiveness and incidence of this kind of surgeries, a 3D reconstruction of the joint from a patient that suffered this kind of injury was used. With this geometrical model, different finite element models including each of the considered alternatives were developed. The simulations of these models, together with the injured and anatomical situations, allowed an approximation of the surgical solution that better restores the initial healthy state of the affected region.<br /

Predicting Surgical Outcome in Patients with Recurrent Patellar Dislocation

Introduction Lateral dislocation of the patella is a common injury in active adolescents and young adults. Patients who are ultimately managed surgically have a significantly lower risk of recurrent dislocation. However, determining the optimal surgical treatment remains a challenge, with patients sometimes undergoing multiple surgeries prior to successful stabilization. The aim of this study is to computationally evaluate patients that have undergone multiple surgeries to correct for recurrent lateral patellar dislocation and predict their clinical outcome. Methods Our patient cohort consisted of 16 patients with patella dislocation. Patient-specific imaging were used to create three-dimensional (3D) finite element (FE) models of the knee joint and evaluate patellofemoral (PF) stability at multiple time points pre- and post-surgery for each patient. We applied these models to predict the clinical success or failure of each surgery. Specifically, the FE model simulated a knee extension activity while a tibia external torsion, a recognized cause of patellofemoral pain and instability, was applied to assess PF stability. A healthy control group of 12 participants was also included to assess the ability of the model to identify successful outcomes. In addition, five anatomic factors of risk were measured, and statistical analysis was performed to establish if significant differences exist among pre-surgery, post-surgery and healthy control groups. Lastly, a logistic regression model was implemented, trained with anatomic values, and used to classify subjects into likelihood of dislocation categories in order to differentiate between successful and unsuccessful surgical outcomes. Feature scaling and feature combination (namely, principal component analysis (PCA)) was applied to improve the predictive performance of the regression model. Results Of 12 control participants, 12 pre-surgery subjects (8 patients after an initial unsuccessful MRPLR and 4 without any), and 9 post-surgery subjects (5 after a successful trochleoplasty and 4 patients after MPFLR), the FE model correctly predicted 29 out of 33 surgery outcomes (87.9% accuracy). Post-surgery simulations predicted patellofemoral stability metrics similar to the healthy control group. Particularly, post-trochleoplasty subjects were associated with an increased ability to provide constraint force on the patella lateral facet, and a lower involvement of the medial patellofemoral ligament, particularly close to full extension. A one-way ANOVA showed that four out of five anatomic factors were significantly different between the pre-surgery and the control group, and three of them also between the pre- and post- surgery group, suggesting that the surgery was able to restore a physiological condition. Lastly, logistic regression classification performance demonstrated 72.2% and 78.9% accuracy before and after PCA, respectively. Conclusion The overall aim of this study is to provide surgeons with a useful and validated computational tool that can predict the likelihood of patellar dislocation and differentiate, prior to clinical intervention, between a successful versus unsuccessful surgery, to determine the optimal treatment pathways for individual patients. Preliminary results are promising, but an improvement of the model and a larger clinical dataset are necessary to improve accuracy and comprehensively validate model performance

Specimen-Specific Natural, Pathological, and Implanted Knee Mechanics Using Finite Element Modeling

There is an increasing incidence of knee pain and injury among the population, and increasing demand for higher knee function in total knee replacement designs. As a result, clinicians and implant manufacturers are interested in improving patient outcomes, and evaluation of knee mechanics is essential for better diagnosis and repair of knee pathologies. Common knee pathologies include osteoarthritis (degradation of the articulating surfaces), patellofemoral pain, and cruciate ligament injury and/or rupture. The complex behavior of knee motion presents unique challenges in the diagnosis of knee pathology and restoration of healthy knee function. Quantifying knee mechanics is essential for developing successful rehabilitation therapies and surgical treatments. Researchers have used in-vitro and in-vivo experiments to quantify joint kinematics and loading, but experiments can be costly and time-intensive, and contact and ligament mechanics can be difficult to measure directly. Computational modeling can complement experimental studies by providing cost-effective solutions for quantifying joint and soft tissue forces. Musculoskeletal models have been used to measure whole-body motion, and predict joint and muscle forces, but these models can lack detail and accuracy at the joint-level. Finite element modeling provides accurate solutions of the internal stress/strain behavior of bone and soft tissue using subject-specific geometry and complex contact and material representations. While previous FE modeling has been used to simulate injury and repair, models are commonly based on literature description or average knee behavior. The research presented in this dissertation focused on developing subject-specific representations of the TF and PF joints including calibration and validation to experimental data for healthy, pathological, and implanted knee conditions. A combination of in-vitro experiment and modeling was used to compare healthy and cruciate-deficient joint mechanics, and develop subject-specific computational representations. Insight from in-vitro testing supported in-vivo simulations of healthy and implanted subjects, in which PF mechanics were compared between two common patellar component designs and the impact of cruciate ligament variability on joint kinematics and loads was assessed. The suite of computational models developed in this dissertation can be used to investigate knee pathologies to better inform clinicians on the mechanisms surrounding injury, support the diagnosis of at-risk patients, explore rehabilitation and surgical techniques for repair, and support decision-making for new innovative implant designs

The use of rapid prototyped implants to simulate knee joint abnormalities for in vitro testing: a validation study with replica implants of the native trochlea

To investigate the biomechanical effect of skeletal knee joint abnormalities, the authors propose to implant pathologically shaped rapid prototyped implants in cadaver knee specimens. This new method was validated by replacing the native trochlea by a replica implant on four cadaver knees with the aid of cadaver-specific guiding instruments. The accuracy of the guiding instruments was assessed by measuring the rotational errors of the cutting planes (on average 3.01 degrees in extension and 1.18 degrees in external/internal rotation). During a squat and open chain simulation, the patella showed small differences in its articulation with the native trochlea and the replica trochlea, which could partially be explained by the rotational errors of the implants. This study concludes that this method is valid to investigate the effect of knee joint abnormalities with a replica implant as a control condition to account for the influence of material properties and rotational errors of the implant

Enhancing rotational placement of reconstruction prostheses of the distal femur after sarcoma resection

Introduction: Currently there are no accepted international guidelines for the correct placement of reconstruction prostheses in the axial plane of the femur after en bloc resection. The most accepted method is based on the linea aspera as an intraoperative landmark, indicating posterior. This study was conducted to address the reliability of the linea aspera as a landmark for rotational alignment. Material and methods: 50 CT angiographies of the right limb were used for this purpose. These 2D images were reconstructed into 3D models using proprietary software (materialize NV, Leuven, Belgium). The posterior condylar line was used as a reference axis. The orientation of the linea aspera was described as the angle between the perpendicular line to the PCL, through the center of the diaphysis, and the lateral (a) and medial labium (beta). Results: The linear mixed model shows that the alpha- and beta-angles are significantly associated with the distance from the joint line (p<0.001) and vary significantly between subjects (p<0.001). The alpha-angle has the lowest variance and approximates more closely true posterior, while the median beta-angle never overlaps true posterior. Discussion: When a surgeon would blindly rely on the linea aspera as a posterior landmark roughly 78% of the femoral implants would exceed the accepted +/- 3 degrees deviation around the surgical transepicondylar axis (sTEA) as defined in total knee replacement. The linea aspera is not a reliable landmark for axial rotation of femoral implants. The position is highly dependent on the osteotomy height and in addition differs between individual patients. Preoperative assessment of the linea aspera is advocated in order to reduce the risk of malrotation. As the height of the osteotomy cannot always be determined correctly preoperatively, a table was designed as a guideline for how much a deviation from the planned resection height will affect the rotation of the implant. (c) 2020 IPEM. Published by Elsevier Ltd. All rights reserved