Predictive structural models resulting in a trabecular bone topology closely resembling real bone would be a step toward 3D printing of sympathetic prosthetics. This study modifies an established trabecular bone structural adaptation approach, with the objective of achieving an improved adapted topology, specifically connectivity, compared to CT imaging studies; whilst retaining continuum level mechanical properties consistent with those reported in experimental studies. Strain driven structural adaptation models successfully identify trabecular trajectories, although tend to overpredict connectivity and skew trabecular radii distribution towards the smallest radius included in the adaptation. Radius adaptation of each trabecula is driven by a mechanostat approach with a target strain (1250 µɛ) below which radius is decreased (resorption), and above which radius is increased (apposition). Simulations include a lazy zone, in which neither resorption nor apposition takes place (1000 to 1500 µɛ); and a dead zone (<250 µɛ) in which complete resorption of trabeculae with the smallest included radius takes place. This study assesses the impact of increasing the dead zone threshold from <250 µɛ to <1000 µɛ, the lower limit of the lazy zone. In-silico structural models with an initial connectivity (number of trabeculae connecting at each joint) of 14 were generated using a nearest neighbour approach applied to a random cloud of points. Trabeculae were modelled using circular beams whose radii were adapted in response to normal strains caused by the axial force and bending moments due to a vertical pressure of 1 MPa applied to the top of the lattice, with the bottom of the lattice fixed in the vertical direction. Lattices in which nodes are either able (rigid jointed) or unable (pin jointed) to transmit bending moments were considered. Five virtual samples of each lattice type were used, and each simulation repeated twice: with a dead zone of either <250 µɛ or <1000 µɛ.Abstract
1.0 Objectives
2.0 Methods
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. 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, Abstract
Objectives
Methods
The screw fastening torque applied during bone fracture fixation has a decisive influence on subsequent bone healing. Insufficient screw tightness can result in device/construct instability; conversely, excessive torques risk damaging the bone causing premature fixation failure. This effect is even more prominent in osteoporotic bone, a condition associated annually with almost 9 million fractures worldwide. During fracture fixation, screw tightening torque is applied using subjective feel. This approach may not be optimal for patient”s recovery, increasing risk of fixation failure, particularly in osteoporotic bone, and potentially require revision surgical interventions. Besides bone density, various factors influence the performance of screw fixation. These factors include bone geometry, cortical thickness and time-dependant relaxation behaviour of the bone. If the influence of screw fastening torque on the bone and relationships between these factors was better understood, the surgical technique could be optimised to reduce the risk of complications. Within this study, we developed an axisymmetric finite element (FE) model of bone screw tightening incorporating viscoelastic behaviour of the cortical bone such as creep and stress relaxation. The model anticipated time-dependent behaviour of the bone for different bone thickness and density after a typical bone fixation screw had been inserted. The idealised model has been developed based on CT scans of bones with varying densities and inserted screws. The model was validated through a series of experiments involving bovine tibiae (4–5 months) to evaluate the evolution of surface strains with time (Ncorr v1.2). Stress distribution was assessed in photoelastic experiments using acrylic analogues. Relaxation tests have been performed in aqueous environment for up to 48 hours to ensure the relaxation would be complete. The creep behaviour (maximum principal strain) was compared against computational predictions. Our early simulations predicted relaxation strains on the surface of the bone to be 1.1% within 24 hours comparing favourably to 1.3% measured experimentally. Stress distribution patterns were in agreement with photoelastic results. Using experimentally derived viscoelastic properties, the model has the potential to predict creep and stress relaxation patterns after screw insertion with different fastening torques for bones with varying density and geometry. We aim to develop this into a planning tool providing guidance to surgeons for optimal tightening when using screw fixation, particularly in reduced quality bone.
The results of the original mobile bearing Oxford unicompartmental knee replacement (UKR) in the lateral compartment have been disappointing because of high dislocation rates (11%). This original implant used a flat bearing articulation on the tibial tray. To address the issue of dislocation a new implant (domed tibia with biconcave bearing to increase entrapment) was introduced with a modified surgical technique. The aim of this study was to compare the risk of dislocation between a domed and flat lateral UKR. Separate geometric computer models of an Oxford mobile bearing lateral UKR were generated for the two types of articulation between the tibial component and the meniscal bearing: Flat-on-flat (flat) and Concave-on-convex (domed). Each type of mobile bearing was used to investigate three distinct dislocation modes observed clinically: lateral to medial dislocation, with the bearing resting on the tray wall (L-M-Wall); medial to lateral dislocation, out of the joint space (M-L); anterior to posterior dislocation, out of the joint space (A-P). A size C tray and a medium femoral component and bearing were used in all models. The femoral component, tibial tray and bearing were first aligned in a neutral position. For each dislocation the tibial tray was restrained in all degrees of freedom. The femoral component was restrained from moving in the anterior-posterior directions and in the medial-lateral directions. The femoral component was also restrained from rotating about the anterior-posterior, medial-lateral and superior-inferior directions. This meant that the femoral component was only able to move in the superior-inferior direction. Different bearing sizes were inserted into the model and the effect that moving the femoral component medially and laterally had on the amount of distraction required to cause bearing dislocation was investigated.Introduction
Methods
Patients with cancer and bone metastases can have an increased risk of fracturing their femur. Treatment is based on the impending fracture risk: patients with a high fracture risk are considered for prophylactic surgery, whereas low fracture risk patients are treated conservatively with radiotherapy to decrease pain. Current clinical guidelines suggest to determine fracture risk based on axial cortical involvement of the lesion on conventional radiographs, but that appears to be difficult. Therefore, we developed a patient-specific finite element (FE) computer model that has shown to be able to predict fracture risk in an experimental setting and in patients. The goal of this study was to determine whether patient-specific finite element (FE) computer models are better at predicting fracture risk for femoral bone metastases compared to clinical assessments based on axial cortical involvement on conventional radiographs, as described in current clinical guidelines. 45 patients (50 affected femurs) affected with predominantly lytic bone metastases who were treated with palliative radiotherapy for pain were included. CT scans were made and patients were followed for six months to determine whether or not they fractured their femur. Non-linear isotropic FE models were created with the patient-specific geometry and bone density obtained from the CT scans. Subsequently, an axial load was simulated on the models mimicking stance. Failure loads normalized for bodyweight (BW) were calculated for each femur. High and low fracture risks were determined using a failure load of 7.5 × BW as a threshold. Experienced assessors measured axial cortical involvement on conventional radiographs. Following clinical guidelines, patients with lesions larger than 30 mm were identified as having a high fracture risk. FE predictions were compared to clinical assessments by means of diagnostic accuracy values (sensitivity, specificity and positive (PPV) and negative predictive values (NPV)). Seven femurs (14%) fractured during follow-up. Median time to fracture was 8 weeks. FE models were better at predicting fracture risk in comparison to clinical assessments based on axial cortical involvement (sensitivity 100% vs. 86%, specificity 74% vs. 42%, PPV 39% vs. 19%, and NPV 100% vs. 95%, for the FE computer model vs. axial cortical involvement, respectively). We concluded that patient-specific FE computer models improve fracture risk predictions of femoral bone metastases in advanced cancer patients compared to clinical assessments based on axial cortical involvement, which is currently used in clinical guidelines. Therefore, we are initiating a pilot for clinical implementation of the FE model.
Total shoulder arthroplasty is a well-tested procedure that offers pain relief and restores the joint function. However, failure rate is still high, and glenoid loosening is pointed as the main reason in orthopedic registers. In order to understand the principles of failure, the principal strain distributions after implantation with Comprehensive® Total Shoulder System of Biomet® were experimental and numerically studied to predict bone behavior. Fourth generation composite left humerus and scapula from Sawbones® were used. These were implanted with Comprehensive® Total Shoulder System (Biomet®) with a modular Hybrid® glenoid base and Regenerex® glenoid and placed in situ by an experienced surgeon. The structures were placed in order to simulate 90º abduction, including principal muscular actions. Muscle forces used were as follows: Deltoideus 300N, Infraspinatus 120N, Supraspinatus 90N, Subscapularis 225N. All bone structures were modeled considering cortical and the trabecular bone of the scapula. The components of prosthesis were placed in the same positions than those in the in vitro models. Geometries were meshed with tetrahedral linear elements, with material properties as follows: Elastic modulus of cortical bone equal to 16 GPa, elastic modulus of trabecular bone equal to 0.155 GPa, polyethylene equal to 1GPa and titanium equal to 110 GPa. The assumed Poisson's ratio was 0.3 in all except for polyethylene where we assumed a value of 0.4. The prosthesis was considered as glued to the adjacent bone. The finite element model was composed of 336 024 elements. At the glenoid cavity, the major influence of the strain distributions was observed at the posterior-superior region, in both cortical and trabecular bone structures. The system presents critical region around holes of fixation in glenoid component. At the trabecular bone, the maximum principal strains at the posterior-superior region ranged from 2250 µε to 3000 µε. While at the cortical bone, the maximum principal strains were 300 µε to 400 µε. The results observed evidence some critical regions of concern and the effect of implant in the bone strains mainly at the posterior-superior region of the glenoid cavity is pronounced. This indicates that this region is more affected by the implant if bone remodeling is a concern and it is due to the strain-shielding effect, which has been connected with loosening of the glenoid component.
Polyethylene wear of joint replacements can cause severe clinical complications, including; osteolysis, implant loosening, inflammation and pain. Wear simulator testing is often used to assess new designs, but it is expensive and time consuming. It is possible to predict the volume of polyethylene implant wear from finite element models using a modification of Archard's classic wear law [1–2]. Typically, linear elastic isotropic, or elasto-plastic material models are used to represent the polyethylene. The purpose of this study was to investigate whether use of a viscoelastic material model would significantly alter the predicted volumetric wear of a mobile-bearing unicompartmental knee replacement. Tensile creep-recovery experiments were performed to characterise the creep and relaxation behaviour of the polyethylene (moulded GUR 4150 samples machined to 180×20×1 mm). Samples were loaded to 3 MPa stress in 4 minutes, and then held for 6 hours, the tensile stress was removed and samples were left to relax for 6 hours. The mechanical test data was used fit to a validated three–dimensional fractional Maxwell viscoelastic constitutive material model [3]. An explicit finite element model of a mobile–bearing unicompartmental knee replacement was created, which has been described previously [4]. The medial knee replacement was loaded to 1200 N over a period of 0.2 s. The bearing was meshed using quadratic tetrahedral elements (1.5 mm seeding size based on results of a mesh convergence study), and the femoral component was represented as an analytical rigid body. Wear predictions were made from the contact stress and sliding distance using Archard's law, as has been described in the literature [1–2]. A wear factor of 5.24×10−11 was used based upon the work by Netter et al. [2]. All models were created and solved using ABAQUS finite element software (version 6.14, Simulia, Dassault Systemes). The fractional viscoelastic material model predicted almost twice as much wear (0.119 mm3/million cycles) compared to the elasto-plastic model (0.069 mm3/million cycles). The higher wear prediction was due to both an increased sliding distance and higher contact pressures in the viscoelastic model. These preliminary findings indicate the simplified elasto-plastic polyethylene material representation can underestimate wear predictions from numerical simulations. Polyethylene is known to be a viscoelastic material which undergoes creep clinically, and it is not surprising that it is necessary to represent that viscoelastic behaviour to accurately predict implant wear. However, it does increase the complexity and run time of such computational studies, which may be prohibitive.
Summary Statement. A coupled finite element - analytical model is presented to predict and to elucidate a clinical healing scenario where bone regenerates in a critical-sized femoral defect, bounded by periosteum or a periosteum substitute implant and stabilised via an intramedullary nail. Introduction. Bone regeneration and maintenance processes are intrinsically linked to mechanical environment. However, the cellular and subcellular mechanisms of mechanically-modulated bone (re-) generation are not fully understood. Recent studies with periosteum osteoprogenitor cells exhibit their mechanosensitivity in vitro and in situ. In addtion, while a variety of growth factors are implicated in bone healing processes, bone morphogenetic protein-2 (BMP-2) is recognised to be involved in all stages of bone regeneration. Furthermore, periosteal injuries heal predominantly via endochondral ossification mechanisms. With this background in mind, the current study aims to understand the role of mechanical environment on BMP-2 production and periosteally-mediated bone regeneration. The one-stage bone transport model [1] provides a clinically relevant experimental platform on which to model the mechanobiological process of periosteum-mediated bone regeneration in a critical-sized defect. Here we develop a model framework to study the cellular-, extracellular- and mechanically-modulated process of defect infilling, governed by the mechanically-modulated production of BMP-2 by osteoprogenitor cells located in the periosteum. Methods. Material properties of the healing callus and periosteum contribute to the strain stimulus sensed by osteoprogenitor cells therein. Using a mechanical finite element model, periosteal surface strains are first predicted as a function of callus properties. Strains are then input to a mechanistic mathematical model, where mechanical regulation of BMP-2 production mediates rates of cellular proliferation, differentiation and extracellular matrix (ECM) production, to predict healing outcomes. A parametric approach enables the spatial and temporal prediction of tissue regeneration via endochondral ossification.
